S. E. Woosley

The evolution of massive stars including mass loss - Presupernova models and explosion

The evolution of massive stars of 35, 40, 60, and 85 solar masses is followed through all stages of nuclear burning to the point of Fe core collapse. Critical nuclear reaction and mass-loss rates are varied. Efficient mass loss during the Wolf-Rayet (WR) stage is likely to lead to final masses as small as 4 solar masses. For a reasonable parameterization of the mass loss, there may be convergence of all WR stars, both single and in binaries, to a narrow band of small final masses. Our representative model, a 4.25 solar-mass WR presupernova derived from a 60 solar mass star, is followed through a simulated explosion, and its explosive nucleosynthesis and light curve are determined. Its properties are similar to those observed in Type Ib supernovae. The effects of the initial mass and mass loss on the presupernova structure of small mass WR models is also explored. Important properties of the presupernova star and its explosion can only be obtained by following the complete evolution starting on the main sequence.

3D simulations of Rayleigh-Taylor mixing in core-collapse SNe with CASTRO

We present multidimensional simulations of the post-explosion hydrodynamics in three different 15 solar mass supernova models with zero, 10^{-4} solar metallicity, and solar metallicities. We follow the growth of the Rayleigh-Taylor instability that mixes together the stellar layers in the wake of the explosion. Models are initialized with spherically symmetric explosions and perturbations are seeded by the grid. Calculations are performed in two-dimensional axisymmetric and three-dimensional Cartesian coordinates using the new Eulerian hydrodynamics code, CASTRO. We find as in previous work, that Rayleigh-Taylor perturbations initially grow faster in 3D than in 2D. As the Rayleigh-Taylor fingers interact with one another, mixing proceeds to a greater degree in 3D than in 2D, reducing the local Atwood number and slowing the growth rate of the instability in 3D relative to 2D. By the time mixing has stopped, the width of the mixed region is similar in 2D and 3D simulations provided the Rayleigh-Taylor fingers show significant interaction. Our results imply that 2D simulations of light curves and nucleosynthesis in supernovae (SNe) that die as red giants may capture the features of an initially spherically symmetric explosion in far less computational time than required by a full 3D simulation. However, capturing large departures from spherical symmetry requires a significantly perturbed explosion. Large scale asymmetries cannot develop through an inverse cascade of merging Rayleigh-Taylor structures; they must arise from asymmetries in the initial explosion.

Evolution and explosion of very massive primordial stars

While the modern stellar IMF shows a rapid decline with increasing mass, theoretical investigations suggest that very massive stars (>100 solar masses) may have been abundant in the early universe. Other calculations also indicate that, lacking metals, these same stars reach their late evolutionary stages without appreciable mass loss. After central helium burning, they encounter the electron-positron pair instability, collapse, and burn oxygen and silicon explosively. If sufficient energy is released by the burning, these stars explode as brilliant supernovae with energies up to 100 times that of an ordinary core collapse supernova. They also eject up to 50 solar masses of radioactive Ni56. Stars less massive than 140 solar masses or more massive than 260 solar masses should collapse into black holes instead of exploding, thus bounding the pair-creation supernovae with regions of stellar mass that are nucleosynthetically sterile. Pair-instability supernovae might be detectable in the near infrared out to redshifts of 20 or more and their ashes should leave a distinctive nucleosynthetic pattern.

The first direct measurement of ¹²C (¹²C,n) ²³Mg at stellar energies

Neutrons produced by the carbon fusion reaction (12)C((12)C,n)(23)Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and improved extrapolation technique based on experimental data from the mirror reaction (12)C((12)C,p)(23)Na. The new reaction rate has been determined with a well-defined uncertainty that exceeds the precision required by astrophysics models. Using our constrained rate, we find that (12)C((12)C,n)(23)Mg is crucial to the production of Na and Al in pop-III pair instability supernovae. It also plays a nonnegligible role in the production of weak s-process elements, as well as in the production of the important galactic γ-ray emitter (60)Fe.

First Direct Measurement of C12(C12,n)Mg23 at Stellar Energies

Neutrons produced by the carbon fusion reaction (12)C((12)C,n)(23)Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and improved extrapolation technique based on experimental data from the mirror reaction (12)C((12)C,p)(23)Na. The new reaction rate has been determined with a well-defined uncertainty that exceeds the precision required by astrophysics models. Using our constrained rate, we find that (12)C((12)C,n)(23)Mg is crucial to the production of Na and Al in pop-III pair instability supernovae. It also plays a nonnegligible role in the production of weak s-process elements, as well as in the production of the important galactic γ-ray emitter (60)Fe.

Fates of the most massive primordial stars

We present our results of numerical simulations of the most massive primordial stars. For the extremely massive non-rotating Pop III stars over 300M⊙, they would simply die as black holes. But the Pop III stars with initial masses 140 - 260M⊙ may have died as gigantic explosions called pair-instability supernovae (PSNe). We use a new radiation-hydrodynamics code CASTRO to study evolution of PSNe. Our models follow the entire explosive burning and the explosion until the shock breaks out from the stellar surface. In our simulations, we find that fluid instabilities occurred during the explosion. These instabilities are driven by both nuclear burning and hydrodynamical instability. In the red supergiant models, fluid instabilities can lead to significant mixing of supernova ejecta and alter the observational signature.