No Arabic abstract
The impact of nuclear physics uncertainties on nucleosynthesis in thermonuclear supernovae has not been fully explored using comprehensive and systematic studies with multiple models. To better constrain predictions of yields from these phenomena, we have performed a sensitivity study by post-processing thermodynamic histories from two different hydrodynamic, Chandrasekhar-mass explosion models. We have individually varied all input reaction and, for the first time, weak interaction rates by a factor of ten and compared the yields in each case to yields using standard rates. Of the 2305 nuclear reactions in our network, we find that the rates of only 53 reactions affect the yield of any species with an abundance of at least 10^-8 M_sun by at least a factor of two, in either model. The rates of the 12C(a,g), 12C+12C, 20Ne(a,p), 20Ne(a,g) and 30Si(p,g) reactions are among those that modify the most yields when varied by a factor of ten. From the individual variation of 658 weak interaction rates in our network by a factor of ten, only the stellar 28Si(b+)28Al, 32S(b+)32P and 36Ar(b+)36Cl rates significantly affect the yields of species in a model. Additional tests reveal that reaction rate changes over temperatures T > 1.5 GK have the greatest impact, and that ratios of radionuclides that may be used as explosion diagnostics change by a factor of less than two from the variation of individual rates by a factor of 10. Nucleosynthesis in the two adopted models is relatively robust to variations in individual nuclear reaction and weak interaction rates. Laboratory measurements of a limited number of reactions would help to further constrain predictions. As well, we confirm the need for a consistent treatment for relevant stellar weak interaction rates since simultaneous variation of these rates (as opposed to individual variation) has a significant effect on yields in our models.
The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and alpha-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucleosynthesis. Temperature-dependent rate uncertainties combining realistic experimental and theoretical uncertainties are used. From detailed statistical analyses uncertainties in the final abundances are derived as probability density distributions. Furthermore, based on rate and abundance correlations an automated procedure identifies the most important reactions in complex flow patterns from superposition of many zones or tracers. The method so far was already applied to a number of nucleosynthesis processes. Here we focus on the production of p-nuclei in white dwarfs exploding as thermonuclear (type Ia) supernovae. We find generally small uncertainties in the final abundances despite of the dominance of theoretical nuclear uncertainties. A separate analysis of low- and high-density regions indicates that the total uncertainties are dominated by the high-density regions.
The explosion energy of thermonuclear (Type Ia) supernovae is derived from the difference in nuclear binding energy liberated in the explosive fusion of light fuel nuclei, predominantly carbon and oxygen, into more tightly bound nuclear ash dominated by iron and silicon group elements. The very same explosive thermonuclear fusion event is also one of the major processes contributing to the nucleosynthesis of the heavy elements, in particular the iron-group elements. For example, most of the iron and manganese in the sun and its planetary system were produced in thermonuclear supernovae. Here, we review the physics of explosive thermonuclear burning in carbon-oxygen white dwarf material and the methodologies utilized in calculating predicted nucleosynthesis from hydrodynamic explosion models. While the dominant explosion scenario remains unclear, many aspects of the nuclear combustion and nucleosynthesis are common to all models and must occur in some form in order to produce the observed yields. We summarize the predicted nucleosynthetic yields for existing explosion models, placing particular emphasis on characteristic differences in the nucleosynthetic signatures of the different suggested scenarios leading to Type Ia supernovae. Following this, we discuss how these signatures compare with observations of several individual supernovae, remnants, and the composition of material in our galaxy and galaxy clusters.
Thermonuclear supernovae originating from the explosion of a white dwarf accreting mass from a companion star have been suggested as a site for the production of $p$ nuclides. Such nuclei are produced during the explosion, in layers enriched with seed nuclei coming from prior strong $s$ processing. These seeds are transformed to proton-richer isotopes mainly by photodisintegration reactions. Several thousand trajectories from a 2D explosion model were used in a Monte Carlo approach. Temperature-dependent uncertainties were assigned individually to thousands of rates varied simultaneously in post-processing in an extended nuclear reaction network. The uncertainties in the final nuclear abundances originating from uncertainties in the astrophysical reaction rates were determined. In addition to the 35 classical $p$ nuclides, abundance uncertainties were also determined for the radioactive nuclides $^{92}$Nb, $^{97,98}$Tc, $^{146}$Sm, and for the abundance ratios $Y$(${}^{92}$Mo)/$Y$(${}^{94}$Mo), $Y$(${}^{92}$Nb)/$Y$(${}^{92}$Mo), $Y$(${}^{97}$Tc)/$Y$(${}^{98}$Ru), $Y$(${}^{98}$Tc)/$Y$(${}^{98}$Ru), and $Y$(${}^{146}$Sm)/$Y$(${}^{144}$Sm), important for Galactic Chemical Evolution studies. Uncertainties found were generally lower than a factor of two, although most nucleosynthesis flows mainly involve predicted rates with larger uncertainties. The main contribution to the total uncertainties comes from a group of trajectories with high peak density originating from the interior of the exploding white dwarf. The distinction between low-density and high-density trajectories allows more general conclusions to be drawn, also applicable to other simulations of white dwarf explosions.
We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to evaluate the reactions that impact more significantly the final abundances. We found that beta-decay rate uncertainties affect typically nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are approximately 50%.
The propagation of uncertainties in reaction cross sections and rates of neutron-, proton-, and $alpha$-induced reactions into the final isotopic abundances obtained in nucleosynthesis models is an important issue in studies of nucleosynthesis and Galactic Chemical Evolution. We developed a Monte Carlo method to allow large-scale postprocessing studies of the impact of nuclear uncertainties on nucleosynthesis. Temperature-dependent rate uncertainties combining realistic experimental and theoretical uncertainties are used. The importance of contributions of cross sections of reactions on excited states of the nuclear targets, which have weights different from from the thermal Boltzmann population factors, is explained. From detailed statistical analyses of the Monte Carlo data uncertainties in the final abundances are derived as probability density distributions. Furthermore, based on rate and abundance correlations an automated procedure identifies the most important reactions in complex flow patterns from superposition of many zones or tracers. The method already has been applied to a number of nucleosynthesis processes.