Correlated Uncertainties in Monte Carlo Reaction Rate Calculations

Abstract

Context. Monte Carlo methods have enabled nuclear reaction rates from uncertain inputs to be presented in a statistically meaningful manner. However, these uncertainties are currently computed assuming no correlations between the physical quantities that enter those calculations. This is not always an appropriate assumption. Astrophysically important reactions are often dominated by resonances, whose properties are usually normalized to a well-known reference resonance. This insight provides a basis from which to develop a flexible framework for including correlations in Monte Carlo reaction rate calculations. Aims. The aim of this work is to develop and test a method for including correlations in Monte Carlo reaction rate calculations when the input has been normalized to a common reference. Methods. A mathematical framework is developed for including correlations between input parameters in Monte Carlo reaction rate calculations. The magnitude of those correlations is calculated from the uncertainties typically reported in experimental papers, where full correlation information is not available. The method is applied to four illustrative examples: a fictional 3-resonance reaction, ${}^{27}$Al(p,$\gamma$)${}^{28}$Si, ${}^{23}$Na(p,$\alpha$)${}^{20}$Ne, and ${}^{23}$Na($\alpha$,p)${}^{26}$Al. Results. Reaction rates at low temperatures that are dominated by a few isolated resonances are found to minimally impacted by correlation effects. However, reaction rates determined from many overlapping resonances can be significantly affected. Uncertainties in the ${}^{23}$Na($\alpha$,p)${}^{26}$Al reaction, for example, increase by up to a factor of 5. This highlights the need to take correlation effects into account in reaction rate calculations, and provide insight into which cases are expected to be most affected by them. The impact of correlation effects on nucleosynthesis is also investigated.

9 pages, 7 figures, accepted to Astronomy and Astrophysics. Typographical errors corrected

References in corpus (7)

• Fast rotating massive stars and the origin of the abundance patterns in galactic globular clusters
• The nucleosynthesis of Al26 and Fe60 in solar metallicity stars extending in mass from 11 to 120 Msun: the hydrostatic and explosive contributions
• On the self-enrichment scenario of galactic globular clusters: Constraints on the IMF
• Statistical Methods for Thermonuclear Reaction Rates and Nucleosynthesis Simulations
• Study of the astrophysically important ${}^{23}\mathrm{Na}(\mathbf{\alpha },\mathbf{p}{)}^{26}\mathrm{Mg}$ and ${}^{23}\mathrm{Na}(\mathbf{\alpha },\mathbf{n}{)}^{26}\mathrm{Al}$ reactions
• Thermonuclear reaction rate of ${}^{18}$Ne($\alpha$,$p$)${}^{21}$Na from Monte-Carlo calculations
• Correlated Uncertainties in Monte Carlo Reaction Rate Calculations

Cited by in corpus (5)

• Correlated Uncertainties in Monte Carlo Reaction Rate Calculations
• Evaluation of the ${}^{13}$N($\alpha$,p)${}^{16}$O thermonuclear reaction rate and its impact on the isotopic composition of supernova grains
• Thermonuclear reaction rate of ${}^{29}$Si(p,$\gamma$)${}^{30}$P
• Correlated Energy Uncertainties in Reaction Rate Calculations
• The 2025 Evaluation of Experimental Thermonuclear Reaction Rates (ETR25)