A simple analytic model for astrophysical S-factors

We propose a physically transparent analytic model of astrophysical S-factors as a function of a center-of-mass energy E of colliding nuclei (below and above the Coulomb barrier) for non-resonant fusion reactions. For any given reaction, the S(E)-model contains four parameters [two of which approximate the barrier potential, U(r)]. They are easily interpolated along many reactions involving isotopes of the same elements; they give accurate practical expressions for S(E) with only several input parameters for many reactions. The model reproduces the suppression of S(E) at low energies (of astrophysical importance) due to the shape of the low-r wing of U(r). The model can be used to reconstruct U(r) from computed or measured S(E). For illustration, we parameterize our recent calculations of S(E) (using the Sao Paulo potential and the barrier penetration formalism) for 946 reactions involving stable and unstable isotopes of C, O, Ne, and Mg (with 9 parameters for all reactions involving many isotopes of the same elements, e.g., C+O). In addition, we analyze astrophysically important 12C+12C reaction, compare theoretical models with experimental data, and discuss the problem of interpolating reliably known S(E) values to low energies (E <= 2-3 MeV).

💡 Research Summary

This paper presents a novel, physically transparent analytic model for the astrophysical S-factor, a key quantity in nuclear astrophysics that describes the probability of fusion reactions in stellar environments after removing the dominant exponential dependence due to Coulomb barrier penetration. The model aims to provide a simple yet accurate parameterization of S(E) as a function of center-of-mass energy E for non-resonant reactions, applicable both below and above the Coulomb barrier.

The core innovation lies in the adoption of a simplified, piecewise effective nucleus-nucleus potential U(r). For large separations, U(r) is the pure Coulomb potential (α/r). For small separations (inside the barrier), it is approximated by an inverted parabola centered at the barrier radius. This specific form allows for the analytic evaluation of the semi-classical barrier penetration integral, which is central to calculating the reaction cross-section at sub-barrier energies. The resulting S-factor is expressed through an explicit formula containing only four parameters with clear physical meanings: S0 (a scale pre-factor), EC (the Coulomb barrier height), δ (a parameter controlling the width of the barrier peak, related to the low-r wing of U(r)), and ξ (a parameter adjusting the behavior above the barrier).

A major strength of the model is the ease of interpolating these parameters across many reactions involving isotopes of the same elements. To demonstrate its practical utility, the authors apply it to parameterize a large dataset of S-factors they previously calculated for 946 fusion reactions involving stable and unstable isotopes of carbon, oxygen, neon, and magnesium. Remarkably, all reactions for a given element pair (e.g., all C+O reactions) can be described with only about 9 parameters, a drastic reduction from the thousands of parameters required by a naive fit, showcasing the model’s efficiency for large-scale astrophysical reaction network calculations.

The model successfully reproduces the suppression (or in some cases, enhancement) of S(E) at very low energies, a phenomenon of paramount astrophysical importance for determining stellar reaction rates. This behavior is directly linked to the shape of U(r) at small internuclear distances, encapsulated by the parameter δ. The paper also highlights that the model can be used in reverse: to constrain the effective interaction potential U(r) from computed or measured S(E) data.

In the final part, the authors focus on the astrophysically critical 12C+12C fusion reaction. They compare predictions from various theoretical models, including their own, with available experimental data. A central discussion point is the significant challenge and inherent uncertainty in extrapolating the S-factor from the measurable energy range (down to about 2-3 MeV) to the even lower energies relevant for carbon burning in stars. The proposed analytic model provides a physically grounded framework for performing such extrapolations more reliably, as its parameters are directly tied to the properties of the potential barrier. In conclusion, this work offers a versatile and intuitive tool for parameterizing, analyzing, and extrapolating astrophysical S-factors, with potential applications in stellar modeling and nucleosynthesis studies.