Astrophysical Reaction Rates for Charged-Particle Induced Reactions on Proton-Rich Nuclides

This paper presents updated astrophysical reaction rates for proton- and alpha-induced reactions on proton-rich nuclei from Neon to Bismuth, calculated using an improved SMARAGD statistical model code, which offer a superior description of experimental data compared to previous rates.

The Gist

Imagine the universe as a giant, cosmic kitchen where stars are the chefs. To cook up the elements that make up our world (like the gold in your jewelry or the oxygen in your lungs), these chefs need to smash tiny particles together at incredibly high speeds and temperatures. This process is called nucleosynthesis.

The problem? We can't just walk into a star and measure how often these particles smash into each other. It's too hot, too far away, and the particles involved are often unstable and disappear before we can study them in a lab.

This paper is like a master recipe book written by a team of theoretical chefs (led by Thomas Rauscher) to help us understand how these cosmic cooking reactions happen, especially when the ingredients are "proton-rich" (a specific type of atomic ingredient).

Here is a breakdown of the paper using simple analogies:

1. The Problem: The "Missing Link" in the Recipe

In the past, scientists had to guess how fast these reactions happened. They used old recipes (models) that were okay, but sometimes they got the taste wrong.

• The Challenge: Some ingredients (nuclei) are so unstable they vanish instantly. Others are so heavy that they are hard to study.
• The Solution: The author used a super-computer program called SMARAGD. Think of SMARAGD as a high-tech "simulator" that predicts how particles will behave when they collide, even if we can't build them in a lab yet.

2. The Method: The "Statistical Crowd" vs. The "Individual Dancer"

To predict these reactions, the paper uses a method called the Hauser-Feshbach model.

• The Analogy: Imagine a crowded dance floor.
  • Low Energy (Direct Reaction): If there are only a few people on the floor, you can watch exactly who bumps into whom. This is easy to predict.
  • High Energy (Statistical Model): If the dance floor is packed with thousands of people spinning and jumping, you can't track individuals. Instead, you look at the average behavior of the crowd. You assume that if the energy is high enough, the particles mix so thoroughly that the specific details of one collision don't matter as much as the overall "vibe" of the crowd.
• Why it matters: For the heavy elements the paper focuses on, the "crowd" is so dense that the statistical approach is the only way to get a good estimate.

3. The New Ingredients: Better Potentials

The old recipes (like the "NON-SMOKER" model) used some rough estimates for how particles interact. The new SMARAGD version updated these estimates with better data:

• The "Alpha" Problem: One of the hardest ingredients to predict is the alpha particle (a chunk of helium nucleus). It's like trying to predict how a heavy bowling ball interacts with a wall when it's moving very slowly. The old models were a bit off.
• The Fix: The author used a new "map" (called the ATOMKI-V2 potential) that better describes how these alpha particles behave, especially when they are moving slowly near the "Coulomb barrier" (a force field that repels them, like trying to push two strong magnets together).
• The Result: The new map fits the few experimental data points we do have much better than the old maps did.

4. The "Ghost" Particles: Excited States

This is a crucial point in the paper.

• The Analogy: Imagine a target in a shooting range.
  • Ground State: The target is standing still and calm.
  • Excited State: The target is vibrating, shaking, and spinning because the room is so hot (like in a star).
• The Insight: In a lab, we usually shoot at the calm target. But inside a star, the targets are always vibrating. The paper explains that if you only look at the calm target, you miss the most important part of the reaction. The new rates account for these "vibrating targets," which changes the cooking speed significantly.

5. The Output: A Better Menu for Astronomers

The paper provides a massive list of numbers (reaction rates) for scientists to use in their simulations of stars and supernovas.

• Why it's better: It's more accurate, especially for reactions involving alpha particles.
• The Caveat: The author warns that while these numbers are the best theoretical guesses we have, they are still guesses. If we ever get better lab data, we should swap the theoretical numbers for real ones. Also, the model works best when the "dance floor" is crowded; if the energy is too low, the "crowd" isn't big enough, and the prediction might be off.

Summary

Think of this paper as updating the cosmic cookbook.

• Old Cookbook: Had some rough estimates and missed how "hot" the ingredients get in a star.
• New Cookbook (SMARAGD): Uses a better simulator, a more accurate map for how particles interact, and accounts for the fact that ingredients in a star are always "jittery."
• Result: We now have a much clearer picture of how the universe cooks up the heavy elements, making our understanding of stellar explosions and the origin of matter more precise.

Technical Summary

1. Problem Statement

Astrophysical simulations of nucleosynthesis (e.g., in the Early Universe, hydrostatic burning, and explosive events like supernovae) rely heavily on accurate nuclear reaction rates. However, direct experimental determination of these rates is often impossible because:

• Many target nuclei are unstable (proton-rich) or have short half-lives.
• Reaction cross-sections at astrophysical energies (often below the Coulomb barrier) are extremely small.
• In stellar plasmas, reactions occur not only on ground states but also on thermally excited states, which are difficult to measure in the lab.

Consequently, most rates used in large-scale astrophysical models are derived from theoretical statistical models (Hauser-Feshbach). Previous widely used compilations (e.g., NON-SMOKER, SMOKER, TALYS) rely on specific nuclear inputs and optical potentials that may not accurately reproduce experimental data, particularly for $\alpha$-induced reactions at low energies. There is a need for an updated, comprehensive rate set that incorporates modern nuclear data and improved theoretical treatments to better describe proton-rich nucleosynthesis (e.g., the $\gamma$-process).

2. Methodology

The author utilized an updated version of the SMARAGD (Statistical Model for Astrophysical Reactions And Global Direct reactions) code, a Hauser-Feshbach statistical model code. The methodology involved several key updates and specific implementations:

• Code Improvements:
  • Implemented a Fox-Goodwin algorithm to solve the Schrödinger equation for charged-particle transmission coefficients, combined with a stable routine for Coulomb wavefunctions. This improves accuracy below the Coulomb barrier.
  • Updated the treatment of parity distributions in the nuclear level density, making them dependent on excitation energy rather than assuming equi-distribution.
• Nuclear Input Data:
  • Masses: Taken from the 2020 Atomic Mass Evaluation (AME2020).
  • Spin/Parity: Ground and low-lying states taken from the 2021 ENSDF+Nubase snapshot (RIPL-3). Up to 40 levels were included where unambiguous assignments exist.
  • Level Density: A hybrid approach using a constant temperature formula for low excitation energies and a backshifted Fermi-gas formula for higher energies.
• Optical Model Potentials (Crucial Updates):
  • Protons: Used a modified version of the Lejeune (1980) potential with an increased depth for the imaginary part to better match experimental data.
  • $\alpha$-particles: Replaced the long-standing McF potential with the ATOMKI-V2 potential. This is a folding potential requiring nuclear matter density distributions, specifically designed to reproduce modern low-energy experimental data.
  • Neutrons: Continued use of the microscopic Lejeune (1980) potential.
• Photon Strength Functions:
  • E1 transitions (Giant Dipole Resonance) use a Lorentzian shape with an energy-dependent width modification at low energies to avoid overestimation.
  • M1 and E2 transitions are treated via single-particle and Giant Quadrupole Resonance models, respectively.
• Scope: Calculations cover target nuclides from Neon (Ne) to Bismuth (Bi), ranging from stability to the proton dripline. Rates are provided for proton- and $\alpha$-induced reactions (essential for proton-rich matter) and, for completeness, neutron-induced reactions on the same targets.

3. Key Contributions

• New Rate Set: The paper presents a comprehensive, machine-readable dataset of astrophysical reaction rates and cross-sections available via Zenodo.
• Improved $\alpha$-Reaction Modeling: By implementing the ATOMKI-V2 potential, the study addresses a long-standing issue where previous models failed to reproduce low-energy $\alpha$-induced reaction data.
• Stellar Cross-Section Formalism: The paper rigorously defines the stellar reaction cross-section ($\sigma^*$), which accounts for thermal excitation of target nuclei. It emphasizes that laboratory ground-state measurements often do not constrain the total astrophysical rate because excited-state contributions can dominate, especially at high temperatures.
• Sensitivity Analysis: The author provides detailed figures (Figs. 1–4) identifying which nuclear widths (particle vs. $\gamma$) dominate specific reaction rates. This helps experimentalists understand which nuclear properties are most critical to measure for specific isotopes.
• Fit Quality Assessment: The paper includes an automated fitting procedure (REACLIB parameterization) for the rates and explicitly maps the deviation between the fits and the calculated values, highlighting regions (e.g., spontaneous $\alpha$-emitters) where tabular values are preferred over fits due to large deviations (40–80%).

4. Results

• Comparison with Experiment: The new SMARAGD rates show significantly better agreement with available experimental data compared to NON-SMOKER and previous evaluations.
  • Proton Capture: The modified proton potential improved the reproduction of energy dependence and scale for experimental data (e.g., Khaliel et al., Harissopoulos et al.).
  • $\alpha$-Induced Reactions: The ATOMKI-V2 potential was essential for reproducing $(\alpha, n)$ data (e.g., on $^{86}$Kr and $^{96}$Zr). Since $(\alpha, n)$ cross-sections depend exclusively on the $\alpha$-width, this validates the new $\alpha$-potential.
• Comparison with Previous Codes:
  • $\alpha$-reactions: The largest differences between SMARAGD and NON-SMOKER arise from the change in the $\alpha$-optical potential. Differences are most pronounced for heavier nuclides where Coulomb barriers are higher.
  • Proton-reactions: Differences are smaller but present due to the modified proton potential and updated nuclear inputs.
  • Level Density: Including discrete experimental levels above the ground state resulted in moderate changes, validating the theoretical level density descriptions used for higher energies.
• Limitations: The statistical model is less applicable near shell closures and the driplines at low temperatures due to low level densities (few overlapping resonances). The paper provides temperature limits for applicability in the data tables.

5. Significance

This work provides the most up-to-date theoretical foundation for modeling nucleosynthesis in proton-rich environments (such as the $\gamma$-process in supernovae and X-ray bursts).

• Accuracy: The improved treatment of $\alpha$-transmission coefficients resolves discrepancies in previous models, leading to more reliable predictions of elemental abundances.
• Guidance for Experiments: By identifying which reaction widths dominate specific rates, the paper guides experimentalists on where to focus efforts to reduce uncertainties in astrophysical models.
• Standardization: The release of a unified, high-quality rate set (including fits and sensitivity data) allows for more consistent and accurate large-scale astrophysical simulations, replacing older, less accurate compilations like NON-SMOKER for charged-particle reactions.

In summary, Rauscher's paper represents a critical update in nuclear astrophysics, bridging the gap between theoretical statistical models and experimental reality through improved optical potentials and rigorous data handling, specifically enhancing the reliability of $\alpha$-induced reaction rates.