Anthony Mezzacappa

Bio: Anthony Mezzacappa is an academic researcher from University of Tennessee. The author has contributed to research in topics: Neutrino & Supernova. The author has an hindex of 52, co-authored 192 publications receiving 9652 citations. Previous affiliations of Anthony Mezzacappa include University of California, Santa Barbara & Joint Institute for Nuclear Research.

Stability of Standing Accretion Shocks, with an Eye toward Core-Collapse Supernovae

Abstract: We examine the stability of standing, spherical accretion shocks. Accretion shocks arise in core-collapse supernovae (the focus of this paper), star formation, and accreting white dwarfs and neutron stars. We present a simple analytic model and use time-dependent hydrodynamics simulations to show that this solution is stable to radial perturbations. In two dimensions we show that small perturbations to a spherical shock front can lead to rapid growth of turbulence behind the shock, driven by the injection of vorticity from the now nonspherical shock. We discuss the ramifications this instability may have for the supernova mechanism.

Protoneutron star evolution and the neutrino-driven wind in general relativistic neutrino radiation hydrodynamics simulations

Abstract: Massive stars end their lives in explosions with kinetic energies on the order of 10 51 erg. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta, which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy, and the material accelerates to an outflow that is known as the neutrino-driven wind. For the first time we simulate the collapse, bounce, explosion, and the neutrino-driven wind phases consistently over more than 20 s. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three-flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino-driven explosions for low- and intermediatemass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the 8. 8M � O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post-explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino-driven wind. On the other hand, we generally find lower neutrino luminosities and mean neutrino energies, as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino-driven wind is proton-rich for more than 10 s and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino-driven wind studies. Despite the moderately high entropies of about 100 kB/baryon and the fast expansion timescales, the conditions found in our models are unlikely to favor r-process nucleosynthesis. The simulations are carried out until the neutrino-driven wind settles down to a quasi-stationary state. About 5 s after the bounce, the peak temperature inside the PNS already starts to decrease because of the continued deleptonization. This moment determines the beginning of a cooling phase dominated by the emission of neutrinos. We discuss the physical conditions of the quasi-static PNS evolution and take the effects of deleptonization and mass accretion from early fallback into account.

Probing the gravitational well: no supernova explosion in spherical symmetry with general relativistic boltzmann neutrino transport

Abstract: We report on the stellar core collapse, bounce, and postbounce evolution of a $13 {\mathrm{M}}_{\ensuremath{\bigodot}}$ star in a self-consistent general relativistic spherically symmetric simulation based on Boltzmann neutrino transport. We conclude that approximations to exact neutrino transport and the omission of general relativistic effects were not alone responsible for the failure of numerous preceding attempts to model supernova explosions in spherical symmetry. Compared to simulations in Newtonian gravity, the general relativistic simulation results in a smaller shock radius. We however argue that the higher neutrino luminosities and rms energies in the general relativistic case could lead to a larger supernova explosion energy.

Supernova simulations with Boltzmann neutrino transport: a comparison of methods

Abstract: Accurate neutrino transport has been built into spherically symmetric simulations of stellar core collapse and postbounce evolution. The results of such simulations agree that spherically symmetric models with standard microphysical input fail to explode by the delayed, neutrino-driven mechanism. Independent groups implemented fundamentally different numerical methods to tackle the Boltzmann neutrino transport equation. Here we present a direct and detailed comparison of such neutrino radiation-hydrodynamics simulations for two codes, AGILE-BOLTZTRAN of the Oak Ridge-Basel group and VERTEX of the Garching group. The former solves the Boltzmann equation directly by an implicit, general relativistic discrete-angle method on the adaptive grid of a conservative implicit hydrodynamics code with second-order TVD advection. In contrast, the latter couples a variable Eddington factor technique with an explicit, moving-grid, conservative high-order Riemann solver with important relativistic effects treated by an effective gravitational potential. The presented study is meant to test our neutrino radiation-hydrodynamics implementations and to provide a data basis for comparisons and verifications of supernova codes to be developed in the future. Results are discussed for simulations of the core collapse and postbounce evolution of a 13 M☉ star with Newtonian gravity and a 15 M☉ star with relativistic gravity.

Three-dimensional Core-collapse Supernova Simulated Using a 15 M⊙ Progenitor

Abstract: We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 M⊙ progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions; however, the onset of explosion (shock revival) is delayed by ~100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. In the ~100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. Finally, we posit that shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas such plumes form more rapidly in 2D, permitting earlier explosions.

First year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Determination of cosmological parameters

Abstract: WMAP precision data enable accurate testing of cosmological models. We find that the emerging standard model of cosmology, a flat � -dominated universe seeded by a nearly scale-invariant adiabatic Gaussian fluctuations, fits the WMAP data. For the WMAP data only, the best-fit parameters are h ¼ 0:72 � 0:05, � bh 2 ¼ 0:024 � 0:001, � mh 2 ¼ 0:14 � 0:02, � ¼ 0:166 þ0:076 � 0:071 , ns ¼ 0:99 � 0:04, and � 8 ¼ 0:9 � 0:1. With parameters fixed only by WMAP data, we can fit finer scale cosmic microwave background (CMB) measure- ments and measurements of large-scale structure (galaxy surveys and the Lyforest). This simple model is also consistent with a host of other astronomical measurements: its inferred age of the universe is consistent with stellar ages, the baryon/photon ratio is consistent with measurements of the (D/H) ratio, and the inferred Hubble constant is consistent with local observations of the expansion rate. We then fit the model parameters to a combination of WMAP data with other finer scale CMB experiments (ACBAR and CBI), 2dFGRS measurements, and Lyforest data to find the model's best-fit cosmological parameters: h ¼ 0:71 þ0:04 � 0:03 , � bh 2 ¼ 0:0224 � 0:0009, � mh 2 ¼ 0:135 þ0:008 � 0:009 , � ¼ 0:17 � 0:06, ns(0.05 Mpc � 1 )=0 :93 � 0:03, and � 8 ¼ 0:84 � 0:04. WMAP's best determination of � ¼ 0:17 � 0:04 arises directly from the temperature- polarization (TE) data and not from this model fit, but they are consistent. These parameters imply that the age of the universe is 13:7 � 0:2 Gyr. With the Lyforest data, the model favors but does not require a slowly varying spectral index. The significance of this running index is sensitive to the uncertainties in the Ly� forest. By combining WMAP data with other astronomical data, we constrain the geometry of the universe, � tot ¼ 1:02 � 0:02, and the equation of state of the dark energy, w < � 0:78 (95% confidence limit assuming w �� 1). The combination of WMAP and 2dFGRS data constrains the energy density in stable neutrinos: � � h 2 < 0:0072 (95% confidence limit). For three degenerate neutrino species, this limit implies that their mass is less than 0.23 eV (95% confidence limit). The WMAP detection of early reionization rules out warm dark matter. Subject headings: cosmic microwave background — cosmological parameters — cosmology: observations — early universe On-line material: color figure

The Observation of Gravitational Waves from a Binary Black Hole Merger

Abstract: On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160) Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

Linear and nonlinear waves, by G. B. Whitham. Pp.636. £50. 1999. ISBN 0 471 35942 4 (Wiley).

Abstract: to be done in this area. Face recognition is a problem that scarcely clamoured for attention before the computer age but, having surfaced, has involved a wide range of techniques and has attracted the attention of some fine minds (David Mumford was a Fields Medallist in 1974). This singular application of mathematical modelling to a messy applied problem of obvious utility and importance but with no unique solution is a pretty one to share with students: perhaps, returning to the source of our opening quotation, we may invert Duncan's earlier observation, 'There is an art to find the mind's construction in the face!'.

Collapsars: Gamma-ray bursts and explosions in 'failed supernovae'

Abstract: Using a two-dimensional hydrodynamics code (PROMETHEUS), we explore the continued evolution of rotating helium stars, Mα 10 M☉, in which iron-core collapse does not produce a successful outgoing shock but instead forms a black hole of 2-3 M☉. The model explored in greatest detail is the 14 M☉ helium core of a 35 M☉ main-sequence star. The outcome is sensitive to the angular momentum. For j16 ≡ j/(1016 cm2 s-1) 3, material falls into the black hole almost uninhibited. No outflows are expected. For j16 20, the infalling matter is halted by centrifugal force outside 1000 km where neutrino losses are negligible. The equatorial accretion rate is very low, and explosive oxygen burning may power a weak equatorial explosion. For 3 j16 20, however, a reasonable value for such stars, a compact disk forms at a radius at which the gravitational binding energy can be efficiently radiated as neutrinos or converted to beamed outflow by magnetohydrodynamical (MHD) processes. These are the best candidates for producing gamma-ray bursts (GRBs). Here we study the formation of such a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter α ≈ 0.001 and 0.1. Infall along the rotational axis is initially uninhibited, and an evacuated channel opens during the first few seconds. Meanwhile the black hole is spun up by the accretion (to a ≈ 0.9), and energy is dissipated in the disk by MHD processes and radiated by neutrinos. For the α = 0.1 model, appreciable energetic outflows develop between polar angles of 30° and 45°. These outflows, powered by viscous dissipation in the disk, have an energy of up to a few times 1051 ergs and a mass ~1 M☉ and are rich in 56Ni. They constitute a supernova-like explosion by themselves. Meanwhile accretion through the disk is maintained for approximately 10-20 s but is time variable (±30%) because of hydrodynamical instabilities at the outer edge in a region where nuclei are experiencing photodisintegration. Because the efficiency of neutrino energy deposition is sensitive to the accretion rate, this instability leads to highly variable energy deposition in the polar regions. Some of this variability, which has significant power at 50 ms and overtones, may persist in the time structure of the burst. During the time followed, the average accretion rate for the standard α = 0.1 and j16 = 10 model is 0.07 M☉ s-1. The total energy deposited along the rotational axes by neutrino annihilation is (1-14) × 1051 ergs, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of energy in the polar regions, at a constant rate of 5 × 1050 ergs s-1 per pole, results in strong relativistic outflow jets beamed to about 1% of the sky. These jets may be additionally modulated by instabilities in the sides of the "nozzle" through which they flow. The jets blow aside the accreting material, remain highly focused, and are capable of penetrating the star in ~10 s. After the jet breaks through the surface of the star, highly relativistic flow can emerge. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible, ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that collapsars will always make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to produce in this model and may require merging neutron stars and black holes for their explanation.