Nuclear Reaction Rate Uncertainties

• Nuclear reaction rate uncertainties are systematic errors in astrophysics that affect stellar evolution, nucleosynthesis, and age dating.
• They are quantified through experimental measurements and Monte Carlo analyses that reveal key sensitivity coefficients and error propagation.
• Reducing these uncertainties via targeted experiments and improved models leads to more accurate predictions of stellar ages and chemical yields.

Nuclear reaction rate uncertainties are a fundamental source of systematic error in nearly all applications of nuclear astrophysics, including stellar evolution, nucleosynthesis, and chronometry. These uncertainties propagate through stellar models to impact predictions of stellar ages, chemical evolution, and nucleosynthetic yields, often dominating the theoretical error budgets that underpin interpretation of high-precision observations. The mechanisms by which nuclear physics uncertainties enter, their propagation and statistical combination, and the methodology to reduce and interpret these uncertainties are central to modern stellar and galactic modeling.

1. Formalism of Nuclear Reaction Rates and Uncertainties

The astrophysical thermonuclear reaction rate per particle pair for a reaction $a + b \rightarrow c + d$ at center-of-mass energy $E$ is typically expressed as

$$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$

where $S(E)$ is the astrophysical $S$-factor encoding nuclear matrix elements and wave-function overlap, $v = \sqrt{2E/\mu}$ is the relative velocity, and $\mu$ is the reduced mass. The rate at stellar temperatures is the Maxwellian average,

$$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$

Uncertainties in $S(E)$, resonance energies, or partial widths propagate into $N_A\langle \sigma v \rangle$.

Formally, for a quantity of interest $E$0 (e.g., rate or abundance), sensitivity coefficients are defined as

$E$1

A fractional uncertainty $E$2 in input $E$3 then induces $E$4 (Rauscher, 2012). For reactions proceeding via isolated resonances, propagation is dominated by the uncertainties in resonance strengths and energies; in statistical model rates, uncertainties in nuclear level densities, optical model potentials, and photon strength functions become dominant.

2. Methods of Quantifying and Propagating Uncertainties

2.1 Laboratory, Solar, and Experimental Constraints

Crucial reactions, such as $E$5H$E$6H in the solar core, have uncertainties constrained by direct experiment, but laboratory measurements at solar energies are hampered by extremely low cross sections. Current best experimental and theoretical studies yield, for the solar pp-fusion $E$7,

$E$8

but this 1% systematic uncertainty induces a $E$9 systematic uncertainty in main-sequence turn-off ages ($$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$0) for low-mass stars (2206.13570). Similar fractional uncertainties propagate proportionally in other pp-chain and CNO-cycle S-factors.

2.2 Helioseismic and Solar Neutrino Constraints

Advanced methods combine helioseismic frequency separation ratios (e.g., $$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$1, $$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$2) with solar neutrino flux measurements to constrain reaction rates in-situ. The sensitivity of separation ratios to core structure allows $$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$3 minimization techniques to yield a new best uncertainty,

$$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$4

a factor of seven improvement over laboratory priors, by integrating over the covariance of observables and using Markov Chain Monte Carlo (MCMC) methods. Solar neutrino fluxes ($$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$5) provide independent constraints, with flux scaling as power-laws in S-factors (e.g., $$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$6) (2206.13570).

2.3 Monte Carlo and Sensitivity Analysis Frameworks

State-of-the-art uncertainty quantification employs large-scale Monte Carlo (MC) sampling, in which thousands of reaction rates are simultaneously varied within their prescribed uncertainty intervals, typically following factor-uncertainty (f.u.) distributions or uniform/log-normal samplings in the logarithm of the variation factor (Fields et al., 2017, Nishimura et al., 2019).

A key advance is simultaneous MC variation, which accounts for interdependence and nonlinear response not captured by varying one rate at a time (Nishimura et al., 2019). For each MC realization, final isotopic abundances are computed, and the resulting distribution quantifies propagated nuclear uncertainties. Key rates are identified using Pearson or Spearman correlation coefficients; a Level-1 key is $$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$7 for abundance $$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$8 and rate $$\sigma(E) = \frac{1}{E}\, S(E) \exp[-2\pi \eta(E)], \quad \eta(E) = Z_a Z_b e^2 / (\hbar v)$$9. Removal or fixing of such rates allows iterative identification of sub-dominant contributors.

3. Impact on Astrophysical Models and Observables

3.1 Stellar Evolution and Ages

Stellar models demonstrate near-linear propagation of S-factor uncertainties into observable quantities such as stellar turn-off ages. For $S(E)$0 in solar-type stars, a 1% rate uncertainty induces a $S(E)$1 error in $S(E)$2; for $S(E)$3 (solar-inferred), the corresponding age uncertainty is $S(E)$4 ($S(E)$5 Myr) for a $S(E)$6 star, an order-of-magnitude reduction over previous estimates (2206.13570).

3.2 Core-Collapse Supernovae and Explosive Nucleosynthesis

In explosive environments (e.g., core-collapse supernovae), temperatures exceed $S(E)$7 GK and chemical equilibria (NSE) dominate. Here, the equilibrium composition depends solely on nuclear partition functions and binding energies, making most reaction-rate uncertainties inconsequential for final iron-group yields, except for isotopes synthesized or frozen out during freeze-out from NSE. A handful of $S(E)$8- and proton-induced rates (e.g., $S(E)$9Ti($S$0,p)$S$1V, $S$2Ca($S$3,$S$4)$S$5Ti) emerge as bottlenecks impacting key radioactive isotopes (e.g., $S$6Ti), with abundance uncertainties of 10–40% even after $S$7 MC samples (Nishimura et al., 3 Nov 2025).

3.3 X-ray Burst Models

Type I X-ray bursts are exquisitely sensitive to rates along the $S$8p- and $S$9-process path. MC and one-at-a-time sensitivity studies consistently identify $v = \sqrt{2E/\mu}$0O$v = \sqrt{2E/\mu}$1, $v = \sqrt{2E/\mu}$2Al$v = \sqrt{2E/\mu}$3, $v = \sqrt{2E/\mu}$4Cu$v = \sqrt{2E/\mu}$5, and $v = \sqrt{2E/\mu}$6Ga$v = \sqrt{2E/\mu}$7 as reactions whose current factor-level uncertainties (up to $v = \sqrt{2E/\mu}$8) fundamentally limit precision constraints on neutron star mass-radius inference from burst light curves and ashes (Sultana et al., 17 Oct 2025, Meisel et al., 2018, Cyburt et al., 2016). Carbon survival post-burst is governed by a set of $v = \sqrt{2E/\mu}$9 and $\mu$0 rates, directly affecting the triggering of superbursts.

3.4 Nucleosynthesis Beyond Iron and Heavy-Element Production

In the i-process, Hauser–Feshbach rates for neutron-rich unstable nuclei carry intrinsic and extrapolation errors up to factors of 10–30 (Denissenkov et al., 2016, Martinet et al., 2023). Propagation leads to 0.5–1.0 dex (i.e., factors of 3–10) uncertainty in surface abundances of $\mu$1 elements in i-process models. Odd-Z elements and actinides (e.g., Th, U) have even higher uncertainty due to single-isotope leverage. Monte Carlo correlation analysis robustly identifies $\mu$2 key (n,$\mu$3) reactions that dominate error budgets in traces such as La, Eu, and Th (Martinet et al., 2023).

4. Statistical and Systematic Aspects of Uncertainties

4.1 Parameter and Model Uncertainties

Uncertainties are categorized as parameter (statistical) or model (systematic). Parameter uncertainties arise from experimental errors in resonance energies, spectroscopic factors, or widths; model uncertainties reflect nuclear reaction models, e.g., choices of nuclear level density or photon strength function (Martinet et al., 2023). Both are assessed by varying input parameters in global fits (e.g., TALYS code) and propagating rate envelopes through stellar or network calculations. For neutron capture on i-process nuclei, maximum-to-minimum MACS uncertainty bands reach factors of 2–60 for $\mu$4 (Martinet et al., 2023).

4.2 Correlated and Uncorrelated Error Structure

Accounting for nontrivial parameter correlations is essential. For example, resonance strengths normalized to a single reference resonance inherit correlated uncertainties; MC sampling requires proper treatment of the covariance matrix, e.g., via Cholesky or eigenvalue decomposition (Longland, 2017, Longland et al., 2020). In complex resonance-dominated rates, ignoring such correlations can misestimate the rate uncertainty by factors of 2–5, particularly in high level-density (many overlapping resonances) regions.

4.3 Sensitivity Maps and Ground-State Contributions

The fraction of the stellar reaction rate traced to ground-state transitions, $\mu$5, and the sensitivity coefficients $\mu$6 as a function of input $\mu$7 and temperature $\mu$8, are essential for targeting experimental efforts (Rauscher, 2012). For (p,$\mu$9) and (n,$$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$0) on stable nuclei, $$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$1 is often 0.6–0.9, but drops sharply for neutron-rich species, rendering laboratory constraints less effective at reducing rate uncertainty.

5. Strategies for Reducing Uncertainties and Prioritizing Measurements

5.1 Experimental Prioritization

Key reactions for uncertainty reduction are identified through MC sensitivity analysis:

• Solar-age determinations: $$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$2H$$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$3H ($$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$4), $$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$5He($$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$6He,2p)$$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$7He, $$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$8Be($$N_A \langle \sigma v \rangle (T) = N_A \sqrt{\frac{8}{\pi\mu}}(kT)^{-3/2} \int_0^{\infty} S(E) \exp\Big[-\frac{E}{kT}-2\pi\eta(E)\Big] dE$$9,$S(E)$0)$S(E)$1B (2206.13570).
• X-ray bursts: $S(E)$2O$S(E)$3, $S(E)$4Al$S(E)$5, $S(E)$6Cu$S(E)$7, $S(E)$8Ga$S(E)$9, $N_A\langle \sigma v \rangle$0O$N_A\langle \sigma v \rangle$1, $N_A\langle \sigma v \rangle$2Ne$N_A\langle \sigma v \rangle$3 (Sultana et al., 17 Oct 2025, Cyburt et al., 2016, Meisel et al., 2018).
• Iron-group production in CCSNe: $N_A\langle \sigma v \rangle$4Ti($N_A\langle \sigma v \rangle$5)$N_A\langle \sigma v \rangle$6V, $N_A\langle \sigma v \rangle$7Ca($N_A\langle \sigma v \rangle$8)$N_A\langle \sigma v \rangle$9Ti, $E$00Ni($E$01)$E$02Cu (Nishimura et al., 3 Nov 2025).
• i-process: $E$03Ba$E$04, $E$05Sm$E$06, $E$07Bi$E$08, $E$09Cs$E$10 (Martinet et al., 2023).

Recommended strategies emphasize direct cross section measurements for high-$E$11 reactions, and indirect or surrogate methods (e.g., transfer reactions, $E$12-strength and level-density measurements) for those dominated by excited-state contributions or with poor ground-state sensitivity.

5.2 Theoretical and Model Developments

Reduction of model uncertainties requires more accurate nuclear-level density and strength function models, improved optical potentials for charged and neutron projectiles, and incorporation of experimental systematics within reaction models. Routine reporting of sensitivity coefficients and model envelopes is necessary to support robust network and model error propagation (Rauscher, 2012).

5.3 Statistical Frameworks and Data Reporting

A unified statistical treatment of uncertainties, explicitly incorporating both parameter and model covariances, is recommended for all future reaction-rate evaluations (Longland, 2017, Longland et al., 2020). Monte Carlo rate libraries must provide not only recommended rates and uncertainty intervals but also covariance matrices or sufficient information to allow correlated sampling.

6. Consequences for Astrophysical Inference and Outlook

Reduction of nuclear reaction rate uncertainties has direct implications for the precision of stellar age dating, nucleosynthetic yield predictions, and the ability to match high-resolution chemical evolution and presolar grain data. For example, the correlation slope $E$13 of $E$14 vs. $E$15 in presolar SiC grains, previously in tension with GCE models, can be fully accounted for by current nuclear uncertainties—particularly in $E$16Si$E$17—demonstrating the astrophysical significance of rate envelope narrowing (Fok et al., 2024).

Future progress will be contingent on targeted experimental campaigns aimed at the most impactful reactions—particularly those on unstable nuclei along nucleosynthetic flow paths—and on continued refinement of statistical and model-based uncertainty quantification tools. As nuclear astrophysics models incorporate a broader range of observables and higher-precision data, the systematic reduction and transparent propagation of nuclear reaction rate uncertainties will remain an essential, enabling component.