New Two-Dimensional Models of Supernova Explosions by the Neutrino-Heating Mechanism: Evidence for Different Instability Regimes in Collapsing Stellar Cores

The neutrino-driven explosion mechanism for core-collapse supernovae in its modern flavor relies on the additional support of hydrodynamical instabilities in achieving shock revival. Two possible candidates, convection and the so-called standing accretion shock instability (SASI), have been proposed for this role. In this paper, we discuss new successful simulations of supernova explosions that shed light on the relative importance of these two instabilities. While convection has so far been observed to grow first in self-consistent hydrodynamical models with multi-group neutrino transport, we here present the first such simulation in which the SASI grows faster while the development of convection is initially inhibited. We illustrate the features of this SASI-dominated regime using an explosion model of a 27 solar mass progenitor, which is contrasted with a convectively-dominated model of an 8.1 solar mass progenitor with subsolar metallicity, whose early post-bounce behavior is more in line with previous 11.2 and 15 solar mass explosion models. We analyze the conditions discriminating between the two different regimes, showing that a high mass-accretion rate and a short advection time-scale are conducive for strong SASI activity. We also briefly discuss some important factors for capturing the SASI-driven regime, such as general relativity, the progenitor structure, a nuclear equation of state leading to a compact proto-neutron star, and the neutrino treatment. Finally, we evaluate possible implications of our findings for 2D and 3D supernova simulations. Our results show that a better understanding of the SASI and convection in the non-linear regime is required.Comment: 12 pages, 13 figures; revised version accepted for publication in Ap

Supernova Simulations from a 3D Progenitor Model -- Impact of Perturbations and Evolution of Explosion Properties

We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using three-dimensional (3D) multi-group neutrino hydrodynamics simulations of an 18 solar mass progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330ms after bounce whereas the shock is not revived in a model based on a spherically symmetric progenitor for at least another 300ms. We tentatively infer a reduction of the critical luminosity for shock revival by ~20% due to pre-collapse perturbations. This indicates that convective seed perturbations play an important role in the explosion mechanism in some progenitors. We follow the evolution of the 18 solar mass model into the explosion phase for more than 2s and find that the cycle of accretion and mass ejection is still ongoing at this stage. With a preliminary value of 0.77 Bethe for the diagnostic explosion energy, a baryonic neutron star mass of 1.85 solar masses, a neutron star kick of ~600km/s and a neutron star spin period of ~20ms at the end of the simulation, the explosion and remnant properties are slightly atypical, but still lie comfortably within the observed distribution. Although more refined simulations and a larger survey of progenitors are still called for, this suggests that a solution to the problem of shock revival and explosion energies in the ballpark of observations are within reach for neutrino-driven explosions in 3D.Comment: 23 pages, 22 figures, accepted for publication in MNRA

Presupernova Evolution of Rotating Massive Stars and the Rotation Rate of Pulsars

Rotation in massive stars has been studied on the main sequence and during helium burning for decades, but only recently have realistic numerical simulations followed the transport of angular momentum that occurs during more advanced stages of evolution. The results affect such interesting issues as whether rotation is important to the explosion mechanism, whether supernovae are strong sources of gravitational radiation, the star's nucleosynthesis, and the initial rotation rate of neutron stars and black holes. We find that when only hydrodynamic instabilities (shear, Eddington-Sweet, etc.) are included in the calculation, one obtains neutron stars spinning at close to critical rotation at their surface -- or even formally in excess of critical. When recent estimates of magnetic torques (Spruit 2002) are added, however, the evolved cores spin about an order of magnitude slower. This is still more angular momentum than observed in young pulsars, but too slow for the collapsar model for gamma-ray bursts.Comment: 10 pages, 2 figures, to appear in Proc. IAU 215 "Stellar Rotation

Simulations of Electron Capture and Low-Mass Iron Core Supernovae

The evolutionary pathways of core-collapse supernova progenitors at the low-mass end of the spectrum are beset with major uncertainties. In recent years, a variety of evolutionary channels has been discovered in addition to the classical electron capture supernova channel of super-AGB stars. The few available progenitor models at the low-mass end have been studied with great success in supernova simulations as the peculiar density structure makes for robust neutrino-driven explosions in this mass range. Detailed nucleosynthesis calculations have been conducted both for models of electron capture supernovae and low-mass iron core supernovae and revealed an interesting production of the lighter trans-iron elements (such as Zn, Sr, Y, Zr) as well as rare isotopes like Ca-48 and Fe-60. We stress the need to explore the low-mass end of the supernova spectrum further and link various observables to understand the diversity of explosions in this regime.Comment: 7 page, 3 figures, proceedings of the conference "The AGB-Supernova Mass Transition", to appear in Memorie della Societ\`a Astronomica Italian

Dependence of X-Ray Burst Models on Nuclear Reaction Rates

X-ray bursts are thermonuclear flashes on the surface of accreting neutron stars and reliable burst models are needed to interpret observations in terms of properties of the neutron star and the binary system. We investigate the dependence of X-ray burst models on uncertainties in (p,$\gamma$), ($\alpha$,$\gamma$), and ($\alpha$,p) nuclear reaction rates using fully self-consistent burst models that account for the feedbacks between changes in nuclear energy generation and changes in astrophysical conditions. A two-step approach first identified sensitive nuclear reaction rates in a single-zone model with ignition conditions chosen to match calculations with a state-of-the-art 1D multi-zone model based on the {\Kepler} stellar evolution code. All relevant reaction rates on neutron deficient isotopes up to mass 106 were individually varied by a factor of 100 up and down. Calculations of the 84 highest impact reaction rate changes were then repeated in the 1D multi-zone model. We find a number of uncertain reaction rates that affect predictions of light curves and burst ashes significantly. The results provide insights into the nuclear processes that shape X-ray burst observables and guidance for future nuclear physics work to reduce nuclear uncertainties in X-ray burst models.Comment: 24 pages, 13 figures, 4 tables, submitte

Nucleosynthesis in Massive Stars With Improved Nuclear and Stellar Physics

We present the first calculations to follow the evolution of all stable nuclei and their radioactive progenitors in stellar models computed from the onset of central hydrogen burning through explosion as Type II supernovae. Calculations are performed for Pop I stars of 15, 19, 20, 21, and 25 M_sun using the most recently available experimental and theoretical nuclear data, revised opacity tables, neutrino losses, and weak interaction rates, and taking into account mass loss due to stellar winds. A novel ``adaptive'' reaction network is employed with a variable number of nuclei (adjusted each time step) ranging from about 700 on the main sequence to more than 2200 during the explosion. The network includes, at any given time, all relevant isotopes from hydrogen through polonium (Z=84). Even the limited grid of stellar masses studied suggests that overall good agreement can be achieved with the solar abundances of nuclei between 16O and 90Zr. Interesting discrepancies are seen in the 20 M_sun model and, so far, only in that model, that are a consequence of the merging of the oxygen, neon, and carbon shells about a day prior to core collapse. We find that, in some stars, most of the ``p-process'' nuclei can be produced in the convective oxygen burning shell moments prior to collapse; in others, they are made only in the explosion. Serious deficiencies still exist in all cases for the p-process isotopes of Ru and Mo.Comment: 53 pages, 17 color figures (3 as separate GIF images), slightly extended discussion and references, accepted by Ap

Sensitivity of the C and O production on the 3-alpha rate

We investigate the dependence of the carbon and oxygen production in stars on the 3-alpha rate by varying the energy of the 02+-state of 12C and determine the resulting yields for a selection of low-mass, intermediate-mass, and massive stars. The yields are obtained using modern stellar evolution codes that follow the entire evolution of massive stars, including the supernova explosion, and consider in detail the 3rd dredge-up process during the thermally pulsating asymptotic giant branch of low-mass and intermediate-mass stars. Our results show that the C and O production in massive stars depends strongly on the initial mass, and that it is crucial to follow the entire evolution. A rather strong C production during the He-shell flashes compared to quiescent He burning leads to a lower sensitivity of the C and O production in low-mass and intermediate-mass stars on the 3-alpha-rate than predicted in our previous work. In particular, the C production of intermediate-mass stars seems to have a maximum close to the actual value of the 02+ energy level of 12C.Comment: Language improved; accepted for publication in Astrophysics and Space Scienc

Critical angular momentum distributions in collapsars: quiescent periods from accretion state transitions in long gamma-ray bursts

The rotation rate in pre-supernova cores is an important ingredient which can profoundly affect the post-collapse evolution and associated energy release in supernovae and long gamma ray bursts (LGRBs). Previous work has focused on whether the specific angular momentum is above or below the critical value required for the creation of a centrifugally supported disk around a black hole. Here, we explore the effect of the distribution of angular momentum with radius in the star, and show that qualitative transitions between high and low angular momentum flow, corresponding to high and low luminosity accretion states, can effectively be reflected in the energy output, leading to variability and the possibility of quiescent times in LGRBs.Comment: 22 pages, 6 figures, 2 Tables, accepted for publication in Ap