No Arabic abstract
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 $ u p$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a large reaction network. As the detailed astrophysical conditions of the $ u p$ process still are unknown, a parameter study was performed, with 23 trajectories covering a large range of entropies and $Y_mathrm{e}$. The resulting abundance uncertainties are given for each trajectory. The $ u p$ process has been speculated to contribute to the light $p$ nuclides but it was not possible so far to reproduce the solar isotope ratios. It is found that it is possible to reproduce the solar $^{92}$Mo/$^{94}$Mo abundance ratio within nuclear uncertainties, even within a single trajectory. The solar values of the abundances in the Kr-Sr region relative to the Mo region, however, cannot be achieved within a single trajectory. They may still be obtained from a weighted superposition of different trajectories, though, depending on the actual conditions in the production site. For a stronger constraint of the required conditions, it would be necessary to reduce the uncertainties in the 3$alpha$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK.
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.
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.
Asymptotic giant branch (AGB) stars with low initial mass (1 - 3 Msun) are responsible for the production of neutron-capture elements through the main s-process (main slow neutron capture process). The major neutron source is 13C(alpha, n)16O, which burns radiatively during the interpulse periods at about 8 keV and produces a rather low neutron density (10^7 n/cm^3). The second neutron source 22Ne(alpha, n)25Mg, partially activated during the convective thermal pulses when the energy reaches about 23 keV, gives rise to a small neutron exposure but a peaked neutron density (Nn(peak) > 10^11 n/cm^3). At metallicities close to solar, it does not substantially change the final s-process abundances, but mainly affects the isotopic ratios near s-path branchings sensitive to the neutron density. We examine the effect of the present uncertainties of the two neutron sources operating in AGB stars, as well as the competition with the 22Ne(alpha, gamma)26Mg reaction. The analysis is carried out on AGB the main-s process component (reproduced by an average between M(AGB; ini) = 1.5 and 3 Msun at half solar metallicity, see Arlandini et al. 1999), using a set of updated nucleosynthesis models. Major effects are seen close to the branching points. In particular, 13C(alpha, n)16O mainly affects 86Kr and 87Rb owing to the branching at 85Kr, while small variations are shown for heavy isotopes by decreasing or increasing our adopted rate by a factor of 2 - 3. By changing our 22Ne(alpha, n)25Mg rate within a factor of 2, a plausible reproduction of solar s-only isotopes is still obtained. We provide a general overview of the major consequences of these variations on the s-path. A complete description of each branching will be presented in Bisterzo et al., in preparation.
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.