Presupernova Evolution Research Articles

Hydrostatic 12C Burning in CO WDs: the Simmering Phase of SNe Ia Progenitors

AbstractAmong the possible progenitor scenarios considered for Type Ia supernovae, those that involve a carbon–oxygen white dwarf (CO WD) accreting stably towards the Chandrasekhar mass should undergo a phase of hydrostatic carbon burning under high densities and strong convection: the simmering or carbon–flash phase, which can extend for a few hundred years before explosion. During this phase the progenitor CO WD can burn a small fraction of its carbon hydrostatically, releasing energy and ashes that make the star convective and able to capture electrons from the degenerate plasma. In this work we present simplified pre–supernova evolution models of CO WDs growing towards the Chandrasekhar mass accreting matter stably and evolving through the simmering phase towards ignition in order to explore the effects of different initial masses and cooling times in the final chemical composition of the WD before explosion. Preliminary results show that, as expected, denser systems at the start of the simmering phase undergo stronger neutronization. The amount of neutronization is less than what is found in the one–zone models of Chamulak et al. (2008), about a third, and can vary by about a factor of two depending on the exact path to explosion.

Neutrino cooling rates due to 54,55,56Fe for presupernova evolution of massive stars

Accurate estimate of neutrino energy loss rates are needed for the study of the late stages of the stellar evolution, in particular for cooling of neutron stars and white dwarfs. The energy spectra of neutrinos and antineutrinos arriving at the Earth can also provide useful information on the primary neutrino fluxes as well as neutrino mixing scenario (it is to be noted that these supernova neutrinos are emitted from a much later stage in stellar evolution than that considered in this manuscript). Proton–neutron quasi-particle random phase approximation (pn-QRPA) theory has recently being used for a microscopic calculation of stellar weak interaction rates of iron isotopes with success. Here I present the detailed calculation of neutrino and antineutrino cooling rates due to key iron isotopes in stellar matter using the pn-QRPA theory. The rates are calculated on a fine grid of temperature-density scale suitable for core-collapse simulators. The calculated rates are compared against earlier calculations. The neutrino cooling rates due to isotopes of iron are in overall good agreement with the rates calculated using the large-scale shell model. During the presupernova evolution of massive stars, from oxygen shell burning till around end of convective core silicon burning phases, the calculated neutrino cooling rates due to 54Fe are three to four times larger than the corresponding shell model rates. The Brink’s hypothesis used in previous calculations can at times lead to erroneous results. The Brink’s hypothesis assumes that the Gamow-Teller strength distributions for all excited states are the same. It is, however, shown by the present calculation that both the centroid and total strength for excited states differ appreciably from the ground state distribution. These changes in the strength distributions of thermally populated excited states can alter the total weak interaction rates rather significantly. The calculated antineutrino cooling rates, due to positron capture and β-decay of iron isotopes, are orders of magnitude smaller than the corresponding neutrino cooling rates and can safely be neglected specially at low temperatures and high stellar densities.

β-DECAY OF KEY TITANIUM ISOTOPES IN STELLAR ENVIRONMENT

Amongst iron regime nuclei, β-decay rates on titanium isotopes are considered to be important during the late phases of evolution of massive stars. The key β-decay isotopes during presupernova evolution were searched from available literature and a microscopic calculation of the decay rates were performed using the proton–neutron quasiparticle random phase approximation (pn-QRPA) theory. As per earlier simulation results, electron capture and β-decay on certain isotopes of titanium are considered to be important for the presupernova evolution of massive stars. Earlier the stellar electron capture rates and neutrino energy loss rates due to relevant titanium isotopes were presented. In this paper we finally present the β-decay rates of key titanium isotopes in stellar environment. The results are also compared against previous calculations. The pn-QRPA β-decay rates are bigger at high stellar temperatures and smaller at high stellar densities compared to the large scale shell model results. This study can prove useful for the core-collapse simulators.

Ground and excited states Gamow-Teller strength distributions of iron isotopes and associated capture rates for core-collapse simulations

This paper reports on the microscopic calculation of ground and excited states Gamow-Teller (GT) strength distributions, both in the electron capture and electron decay direction, for 54,55,56Fe. The associated electron and positron capture rates for these isotopes of iron are also calculated in stellar matter. These calculations were recently introduced and this paper is a follow-up which discusses in detail the GT strength distributions and stellar capture rates of key iron isotopes. The calculations are performed within the framework of the proton-neutron quasiparticle random phase approximation (pn-QRPA) theory. The pn-QRPA theory allows a microscopic state-by-state calculation of GT strength functions and stellar capture rates which greatly increases the reliability of the results. For the first time experimental deformation of nuclei are taken into account. In the core of massive stars isotopes of iron, 54,55,56Fe, are considered to be key players in decreasing the electron-to-baryon ratio (Y e ) mainly via electron capture on these nuclide. The structure of the presupernova star is altered both by the changes in Y e and the entropy of the core material. Results are encouraging and are compared against measurements (where possible) and other calculations. The calculated electron capture rates are in overall good agreement with the shell model results. During the presupernova evolution of massive stars, from oxygen shell burning stages till around end of convective core silicon burning, the calculated electron capture rates on 54Fe are around three times bigger than the corresponding shell model rates. The calculated positron capture rates, however, are suppressed by two to five orders of magnitude.

Expanded calculation of weak-interaction-mediated neutrino cooling rates due to 56Ni in stellar matter

An accurate estimate of the neutrino cooling rates is required in order to study the various stages of stellar evolution of massive stars. Neutrino losses from proto-neutron stars play a crucial role in deciding whether these stars would be crushed into black holes or explode as supernovae. Both pure leptonic and weak-interaction processes contribute to the neutrino energy losses in stellar matter. At low temperatures and densities, the characteristics of the early phase of presupernova evolution, cooling through neutrinos produced via the weak interaction, are important. Proton–neutron quasi-particle random phase approximation (pn-QRPA) theory has recently been used with success for the calculation of stellar weak-interaction rates of fp-shell nuclide. The lepton-to-baryon ratio (Ye) during early phases of stellar evolution of massive stars changes substantially, mainly due to electron captures on 56Ni. The stellar matter is transparent to the neutrinos produced during the presupernova evolution of massive stars. These neutrinos escape the site and assist the stellar core in maintaining a lower entropy. Here, an expanded calculation of weak-interaction-mediated neutrino and antineutrino cooling rates due to 56Ni in stellar matter using the pn-QRPA theory is presented. This detailed scale is appropriate for interpolation purposes and is of greater utility for simulation codes. The calculated rates are compared with earlier calculations. During the relevant temperature and density regions of stellar matter the reported rates show few differences compared with the shell model rates and might contribute in fine-tuning of the lepton-to-baryon ratio during the presupernova phases of stellar evolution of massive stars.

Presupernova evolution and explosion of massive stars

We review our recent progresses on the presupernova evolution of massive stars in the range 11-120 M⊙ of solar metallicity. Special attention will be devoted to the effect of the mass loss rate during the Wolf-Rayet stages in determining the structure and the physical properties of the star prior the supernova explosion. We also discuss the explosive yields and the initial mass-remnant mass relation in the framework of the kinetic bomb induced explosion and hence the contribution of these stars to the global chemical enrichment of the interstellar medium

Weak-interaction-mediated rates on iron isotopes for presupernova evolution of massive stars

During the presupernova evolution of massive stars, the isotopes of iron, $^{54,55,56}$Fe, are advocated to play a key role inside the cores primarily decreasing the electron-to-baryon ratio ($Y_{e}$) mainly via electron capture processes thereby reducing the pressure support. Electron decay and positron capture on $^{55}$Fe, on the other hand, also has a consequential role in increasing the lepton ratio during the silicon burning phases of massive stars. The neutrinos and antineutrinos produced, as a result of these weak-interaction reactions, are transparent to the stellar matter and assist in cooling the core thereby reducing the entropy. The structure of the presupernova star is altered both by the changes in $Y_{e}$ and the entropy of the core material. The aim of this paper is to report the improved microscopic calculation of Gamow-Teller (GT$_{\pm}$) strength distributions of these key isotopes of iron using the pn-QRPA theory. The main improvement comes from the incorporation of experimental deformation values for these nuclei. Additionally six different weak-interaction rates, namely electron & positron capture, electron & positron decay, and, neutrino & antineutrino cooling rates, were also calculated in presupernova matter. The calculated electron capture and neutrino cooling rates due to isotopes of iron are in good agreement with the large-scale shell model (LSSM) results. The calculated beta decay rates, however, are suppressed by three to five orders of magnitude.

GAMMA-RAY BURST PRODUCTION AND SUPERNOVA SIGNATURES IN SLOWLY ROTATING COLLAPSARS

We consider accretion onto newborn black holes following the collapse of rotating massive stellar cores, at the threshold where a centrifugally supported disk gives way to nearly radial inflow for low angular momentum. For realistic initial conditions taken from pre-supernova (pre-SN) evolution calculations, the densities and temperatures involved require the use of a detailed equation of state and neutrino cooling processes, as well as a qualitative consideration of the effects of general relativity. Through two-dimensional dynamical calculations we show how the energy release is affected by the rotation rate and the strength of angular momentum transport, giving rise to qualitatively different solutions in limits of high and low angular momentum, each being capable of powering a gamma-ray burst (GRB). We explore the likelihood of producing Fe-group elements in the two regimes and suggest that while large and massive centrifugally supported disks are capable of driving strong outflows with a possible SN-like signature, quasi-radial flows lack such a feature and may produce a GRB without such an accompanying feature, as seen in GRB060505.

Neutrino energy loss rates and positron capture rates onCo55for presupernova and supernova physics

Proton-neutron quasi-particle random phase approximation (pn-QRPA) theory has recently being used for calculation of stellar weak interaction rates of $fp$-shell nuclide with success. Neutrino losses from proto-neutron stars play a pivotal role to decide if these stars would be crushed into black holes or explode as supernovae. The product of abundance and positron capture rates on $^{55}$Co is substantial and as such can play a role in fine tuning of input parameters of simulation codes specially in the presupernova evolution. Recently we introduced our calculation of capture rates on $^{55}$Co, in a luxurious model space of $7 \hbar \omega$, employing the pn-QRPA theory with a separable interaction. Simulators, however, may require these rates on a fine scale. Here we present for the first time an expanded calculation of the neutrino energy loss rates and positron capture rates on $^{55}$Co on an extensive temperature-density scale. These type of scale is appropriate for interpolation purposes and of greater utility for simulation codes. The pn-QRPA calculated neutrino energy loss rates are enhanced roughly up to two orders of magnitude compared with the large-scale shell model calculations and favor a lower entropy for the core of massive stars.

Final stages of massive stars. SN explosion and explosive nucleosynthesis

We review the presupernova evolution, the explosion and the explosive nucleosynthesis of massive stars in the range 11–120 M ๏ . We mainly focus on the advanced evolutionary phases prior to the explosion, i.e. from the core He exhaustion up to the onset of the iron core collapse.

On the Sensitivity of Massive Star Nucleosynthesis and Evolution to Solar Abundances and to Uncertainties in Helium‐Burning Reaction Rates

We explore the dependence of presupernova evolution and supernova nucleosynthesis yields on the uncertainties in helium-burning reaction rates. Using the revised solar abundances of Lodders for the initial stellar composition, instead of those of Anders and Grevesse, changes the supernova yields and limits the constraints that those yields place on the 12C(α,γ)16O reaction rate. The production factors of medium-weight elements (A = 16-40) were found to be in reasonable agreement with observed solar ratios within the current experimental uncertainties in the triple-α reaction rate. Simultaneous variations by the same amount in both reaction rates or in either of them separately, however, can induce significant changes in the central 12C abundance at core carbon ignition and in the mass of the supernova remnant. It therefore remains important to have experimental determinations of the helium-burning rates so that their ratio and absolute values are known with an accuracy of 10% or better.

A Massive Cluster of Red Supergiants at the Base of the Scutum‐Crux Arm

We report on the unprecedented red supergiant (RSG) population of a massive young cluster, located at the base of the Scutum-Crux Galactic arm. We identify candidate cluster RSGs based on 2MASS photometry and medium-resolution spectroscopy. With follow-up high-resolution spectroscopy, we use CO band-head equivalent width and high-precision radial velocity measurements to identify a core grouping of 26 physically associated RSGs—the largest such cluster known to date. Using the stars' velocity dispersion and their inferred luminosities in conjunction with evolutionary models, we argue that the cluster has an initial mass of ~40,000 M☉ and is therefore among the most massive in the galaxy. Further, the cluster is only a few hundred parsecs away from the cluster of 14 RSGs recently reported by Figer at al.. These two RSG clusters represent 20% of all known RSGs in the Galaxy, and now offer the unique opportunity to study the presupernova evolution of massive stars, and the blue- to red-supergiant ratio at uniform metallicity. We use GLIMPSE, MIPSGAL, and MAGPIS survey data to identify several objects in the field of the larger cluster which seem to be indicative of recent regionwide starburst activity at the point where the Scutum-Crux arm intercepts the Galactic bulge. Future abundance studies of these clusters will therefore permit the study of the chemical evolution and metallicity gradient of the Galaxy in the region where the disk meets the bulge.

Are the Models for Type Ia Supernova Progenitors Consistent with the Properties of Supernova Remnants?

We explore the relationship between the models for progenitor systems of Type Ia supernovae and the properties of the supernova remnants that evolve after the explosion. Most models for Type Ia progenitors in the single degenerate scenario predict substantial outflows during the presupernova evolution. Expanding on previous work, we estimate the imprint of these outflows on the structure of the circumstellar medium at the time of the supernova explosion, and the effect that this modified circumstellar medium has on the evolution of the ensuing supernova remnant. We compare our simulations with the observational properties of known Type Ia supernova remnants in the Galaxy (Kepler, Tycho, SN 1006), the Large Magellanic Cloud (0509-67.5, 0519-69.0, N103B), and M31 (SN 1885). We find that optically thick outflows from the white dwarf surface (sometimes known as accretion winds) with velocities above 200 km/s excavate large low-density cavities around the progenitors. Such large cavities are incompatible with the dynamics of the forward shock and the X-ray emission from the shocked ejecta in all the Type Ia remnants that we have examined.

Off‐Center Ignition in Type Ia Supernovae. I. Initial Evolution and Implications for Delayed Detonation

The explosion of a carbon-oxygen white dwarf as a Type Ia supernova is known to be sensitive to the manner in which the burning is ignited. Studies of the pre-supernova evolution suggest asymmetric, off-center ignition, and here we explore its consequences in two- and three-dimensional simulations. Compared with centrally ignited models, one-sided ignitions initially burn less and release less energy. For the distributions of ignition points studied, ignition within two hemispheres typically leads to the unbinding of the white dwarf, while ignition within a small fraction of one hemisphere does not. We also examine the spreading of the blast over the surface of the white dwarf that occurs as the first plumes of burning erupt from the star. In particular, our studies test whether the collision of strong compressional waves can trigger a detonation on the far side of the star as has been suggested by Plewa et al. (2004). The maximum temperature reached in these collisions is sensitive to how much burning and expansion has already gone on, and to the dimensionality of the calculation. Though detonations are sometimes observed in 2D models, none ever happens in the corresponding 3D calculations. Collisions between the expansion fronts of multiple bubbles also seem, in the usual case, unable to ignite a detonation. "Gravitationally confined detonation" is therefore not a robust mechanism for the explosion. Detonation may still be possible in these models however, either following a pulsation or by spontaneous detonation if the turbulent energy is high enough.

Beta decay rates of nuclei with 65 < A < 75 for pre-supernova and supernova evolution

The half-lives are calculated for the β− decay process for nuclei in the mass range ∼65–75 relevant for the core of a massive star at the late burning stage of stellar evolution and the collapse that leads to supernova explosion. These half-lives and rates are calculated by expressing the β− Gamow-Teller decay strengths in terms of smoothed bivariate strength densities. These strength densities are constructed in the framework of spectral averaging theory for two-body nuclear Hamiltonian in a large nuclear shell model space. The method has a natural extension to electron captures as well as weak interaction rates for r and rp-processes.

Discovery of a Nearby Twin of SN 1987A's Nebula around the Luminous Blue Variable HD 168625: Was Sk -69 202 an LBV?

Some core-collapse supernovae show clear signs of interaction with dense circumstellar material that often appears to be non-spherical. Circumstellar nebulae around supernova progenitors provide clues to the origin of that asymmetry in immediate pre-supernova evolution. Here I discuss outstanding questions about the formation of the ring nebula around SN1987A and some implications of similar ring nebulae around Galactic B supergiants. Several clues hint that SN1987A's nebula may have been ejected in an LBV-like event, rather than through interacting winds in a transition from a red supergiant to a blue supergiant.

Evolution of the First Stellar Generations

Although the theoretical study of very low metallicity (Z) and metal-free stars is not new, their importance has recently greatly increased since two related fields have been developing rapidly. The first is cosmological simulations of the formation of the first stars and of the reionisation period. The second is the observations of extremely metal poor stars. In this paper, we present pre-supernova evolution models of massive rotating stars at very low Z (Z=1e-8) and at Z=0. Rotation has a strong impact on mass loss and nucleosynthesis. Models reaching break-up velocities lose up to ten percents of their initial mass. In very low Z models, rotational and convective mixing enhances significantly the surface content in carbon, nitrogen and oxygen (CNO) when the star becomes a red supergiant. This induces a strong mass loss for stars more massive than about 60 solar masses. Our models predict type Ib,c supernovae and gamma-ray bursts at very low Z. Rotational mixing also induces a large production of CNO elements, in particular of primary nitrogen. The stellar wind chemical composition is compatible with the most metal-poor star know to date, HE 1327-2326, for CNO elements. Our models reproduce the early evolution of nitrogen in the Milky Way.

Core‐Collapse Very Massive Stars: Evolution, Explosion, and Nucleosynthesis of Population III 500–1000M⊙Stars

We calculate evolution, collapse, explosion, and nucleosynthesis of Population III very massive stars with 500 and 1000 M☉. Presupernova evolution is calculated in spherical symmetry. Collapse and explosion are calculated by a two-dimensional code, based on the bipolar jet models. We compare the results of nucleosynthesis with the abundance patterns of intracluster matter, hot gases in M82, and extremely metal-poor stars in the Galactic halo. It was found that both 500 and 1000 M☉ models enter the region of pair instability but continue to undergo core collapse. In the presupernova stage, silicon-burning regions occupy a large fraction, more than 20% of the total mass. For moderately aspherical explosions, the patterns of nucleosynthesis match the observational data of both the intracluster medium and M82. Our results suggest that explosions of Population III core-collapse very massive stars contribute significantly to the chemical evolution of gases in clusters of galaxies. For Galactic halo stars our [O/Fe] ratios are smaller than the observational abundances. However, our proposed scenario is naturally consistent with this outcome. The final black hole masses are ~230 and ~500 M☉ for the 500 and 1000 M☉ models, respectively. This result may support the view that Population III very massive stars are responsible for the origin of intermediate-mass black holes, which were recently reported to be discovered.

Abundance Ratios in the Galactic Bulge and Super Metal-Rich Type II Nucleosynthesis: Pitfalls of the Analysis

We present abundance results from our Keck/HIRES observations of giants in the Galactic Bulge. We confirm that the metallicity distribution of giants in the low-reddening bulge field Baade's Window can be well-fit by a closed-box enrichment model. We also confirm previous observations that find enhanced [Mg/Fe], [Si/Fe] and [Ca/Fe] for all bulge giants, including those at super-solar metallicities. However, we find that the [O/Fe] ratios of metal-rich bulge dwarfs decrease with increasing metallicity, contrary to what is expected if the enhancements of the other a-elements is due to Type II supernovae enrichment. We suggest that the decrease in oxygen production may be due to mass loss in the pre-supernova evolution of metal-rich progenitors.

Presupernova Evolution of Differentially Rotating Massive Stars Including Magnetic Fields

As a massive star evolves through multiple stages of nuclear burning on its way to becoming a supernova, a complex, differentially rotating structure is set up. Angular momentum is transported by a variety of classic instabilities and also by magnetic torques from fields generated by the differential rotation. We present the first stellar evolution calculations to follow the evolution of rotating massive stars including, at least approximately, all these effects, magnetic and nonmagnetic, from the zero-age main sequence until the onset of iron-core collapse. The evolution and action of the magnetic fields is as described by Spruit in 2002, and a range of uncertain parameters is explored. In general, we find that magnetic torques decrease the final rotation rate of the collapsing iron core by about a factor of 30-50 when compared with the nonmagnetic counterparts. Angular momentum in that part of the presupernova star destined to become a neutron star is an increasing function of main-sequence mass. That is, pulsars derived from more massive stars rotate faster and rotation plays a more important role in the star's explosion. The final angular momentum of the core has been determined—to within a factor of 2—by the time the star ignites carbon burning. For the lighter stars studied, around 15 M☉, we predict pulsar periods at birth near 15 ms, though a factor of 2 range is easily tolerated by the uncertainties. Several mechanisms for additional braking in a young neutron star, especially by fallback, are explored.