Astronomy with Radioactivities: Chapter 9, Nuclear Reactions

Nuclear reaction rates determine the abundances of isotopes in stellar burning processes. A multitude of reactions determine the reaction flow pattern which is described in terms of reaction network simulations. The reaction rates are determined by laboratory experiments supplemented by nuclear reaction and structure theory. We will discuss the experimental approach as well as the theoretical tools for obtaining the stellar reaction rates. A detailed analysis of a reaction is only possible for a few selected cases which will be highlighted in this section. The bulk of nuclear reaction processes is however described in terms of a statistical model approach, which relies on global nuclear structure and reaction parameters such as level density and mass and barrier penetration, respectively. We will discuss a variety of experimental facilities and techniques used in the field, this includes low energy stable beam experiments, measurements at radioactive beam accelerators, and neutron beam facilities.

💡 Research Summary

This chapter provides a comprehensive overview of how nuclear reaction rates shape the isotopic abundances produced during stellar burning and, consequently, influence stellar evolution, supernova explosions, and the synthesis of heavy elements. It begins by defining the stellar reaction rate, N_A⟨σv⟩, as the Maxwell‑Boltzmann‑averaged product of the nuclear cross‑section σ(E) and the relative velocity v of reacting particles. Because stellar interiors operate at temperatures ranging from a few million to several billion kelvin and densities up to 10⁸ g cm⁻³, the reacting nuclei possess relatively low kinetic energies, making quantum tunnelling through the Coulomb barrier the dominant mechanism. The probability of tunnelling is highly sensitive to the barrier height, the reduced mass of the system, and the energy distribution of the particles, which together determine the temperature dependence of the reaction rate.

The experimental determination of reaction rates is divided into three principal approaches. First, low‑energy stable‑beam experiments employ conventional accelerators, high‑purity germanium or LaBr₃ γ‑detectors, and recoil separators to directly measure cross‑sections at energies as close as possible to the astrophysical Gamow window. Second, radioactive‑ion beam (RIB) facilities enable the study of reactions involving short‑lived nuclei that cannot be produced as stable targets. Techniques such as inverse kinematics, Coulomb dissociation, the Trojan‑Horse method, and the extraction of Asymptotic Normalization Coefficients (ANCs) are used to infer the desired (α,γ), (p,γ), or (n,γ) rates. Third, neutron‑beam facilities provide high‑flux, pulsed neutron sources combined with time‑of‑flight spectrometers to measure (n,γ) cross‑sections crucial for the slow neutron‑capture (s‑process) path. The chapter highlights the importance of high‑intensity facilities such as FRIB (USA), RIKEN (Japan), FAIR (Germany), and the upcoming SPARC neutron source for extending measurements to exotic isotopes far from stability.

On the theoretical side, the Hauser–Feshbach statistical model forms the backbone of large‑scale reaction‑rate calculations. When the level density of the compound nucleus is sufficiently high, individual resonances overlap and the reaction can be described by averaged transmission coefficients. The model requires global inputs: nuclear mass models (e.g., FRDM, HFB‑21), level‑density prescriptions (constant‑temperature, back‑shifted Fermi‑gas, microscopic combinatorial), optical‑model potentials for neutrons, protons, and α‑particles, and γ‑strength functions. These inputs are calibrated against the growing body of experimental data using Bayesian inference or machine‑learning techniques, allowing the generation of extensive reaction libraries such as REACLIB and JINA’s STARLIB. The chapter emphasizes that while the statistical approach works well for the majority of reactions, a limited set of “key reactions” (e.g., ¹²C(α,γ)¹⁶O, ³He(α,γ)⁷Be, ⁷Be(p,γ)⁸B) dominate the nucleosynthetic flow and therefore demand high‑precision, often direct, measurements.

The integration of experimental data and theoretical rates into reaction‑network simulations is discussed in detail. Networks solve a coupled set of differential equations that track the time evolution of thousands of isotopic abundances under prescribed temperature–density trajectories. Sensitivity studies identify which reactions most strongly affect final abundances, guiding experimental priorities. Uncertainty quantification, typically performed via Monte‑Carlo sampling of reaction‑rate probability density functions, reveals that a handful of reactions contribute disproportionately to the overall error budget in predictions of, for example, solar neutrino fluxes or the r‑process abundance peaks.

Finally, the chapter looks ahead to emerging technologies. Next‑generation RIB facilities will deliver higher beam intensities and broader isotope coverage, enabling direct measurements of reactions previously accessible only through indirect methods. Advanced detector arrays with γ‑ray tracking (e.g., GRETA, AGATA) and active target time‑projection chambers promise improved angular resolution and efficiency, crucial for low‑cross‑section experiments. Coupled with ever‑more sophisticated statistical models and high‑performance computing, these developments are expected to reduce the uncertainties in stellar reaction rates dramatically, thereby sharpening our understanding of the chemical evolution of the universe.