Folding model approach to the elastic $p+^{12,13}$C scattering at low energies and radiative capture $^{12,13}$C$(p,γ)$ reactions

This paper employs a consistent folding model with a realistic density-dependent nucleon-nucleon interaction to simultaneously describe both elastic $p+^{12,13}$C scattering and the astrophysical $S$ factors of the radiative capture $^{12,13}$C$(p,\gamma)$ reactions at low energies.

The Gist

The Big Picture: Cooking in the Stars

Imagine stars as giant cosmic kitchens. The Sun is a relatively simple kitchen, but heavier stars are like high-pressure, high-temperature pressure cookers. To keep these massive stars from collapsing, they need to generate huge amounts of energy. They do this by fusing hydrogen into helium, but they need a special "starter" or "catalyst" to make the process efficient.

This catalyst is Carbon. The paper focuses on two specific recipes in the star's kitchen:

1. Recipe A: A proton (a hydrogen nucleus) hits a Carbon-12 atom to make Nitrogen-13.
2. Recipe B: A proton hits a Carbon-13 atom to make Nitrogen-14.

These reactions are the first steps in the CNO cycle (Carbon-Nitrogen-Oxygen cycle), which is the main way massive stars burn fuel. If we don't know exactly how fast these reactions happen, we can't accurately predict how stars live, die, or how heavy elements are created in the universe.

The Problem: The "Black Box" of Nuclear Physics

Scientists want to calculate the speed of these reactions. However, inside an atom, things are messy. You have protons and neutrons zipping around, interacting in complex ways. Trying to calculate every single interaction between every particle is like trying to predict the exact path of every single grain of sand in a hurricane. It's too complicated.

So, scientists use a model. Think of it like a map. You don't need to know the texture of every single tree in a forest to navigate it; you just need a good map of the roads. In nuclear physics, this "map" is called a potential. It's a mathematical description of how the proton "feels" the pull of the Carbon nucleus as it approaches.

The Old Way vs. The New Way

The Old Way (Phenomenological Models):
Previously, scientists built their maps by guessing the shape of the road and then tweaking the numbers until the map matched the data they had. It worked, but it was a bit like drawing a map based on a blurry photo and then erasing and redrawing lines until it looked right. It had a lot of "free parameters" (adjustable knobs) that didn't necessarily reflect the real physics.

The New Way (The Folding Model):
The authors of this paper used a method called the Folding Model.

• The Analogy: Imagine you have a bag of marbles (the Carbon nucleus) and you want to know how a single new marble (the proton) will roll over them. Instead of guessing the shape of the pile, you take the actual shape of the bag of marbles and "fold" the force of the new marble over it. You calculate the total pull based on the actual density of the marbles inside.
• The Benefit: This method uses real physics (how protons and neutrons actually interact) and has very few adjustable knobs. It's a "microscopic" approach, meaning it looks at the tiny details to build the big picture.

What Did They Do?

The team took this "Folding Model" map and tested it in two ways:

1. The Bounce Test (Elastic Scattering):
   First, they fired protons at Carbon atoms and watched them bounce off (like billiard balls). They compared the real bounce data with their Folding Model map.

   • Result: The map worked great! It predicted exactly how the protons would bounce, provided they tweaked the "strength" of the map by about 20-30%. This proved their map was accurate.

2. The Capture Test (Radiative Capture):
   Next, they used the same map to predict what happens when the proton doesn't bounce, but instead gets stuck to the Carbon atom, releasing a flash of light (a gamma ray). This is the actual "cooking" reaction.

   • Result: The map worked well for the Carbon-12 reaction. It predicted the energy output and the rate of the reaction very accurately.
   • The Hiccup: For the Carbon-13 reaction, the map got the general trend right but missed a specific, very sharp "spike" in the data (a resonance). It's like their map showed a smooth hill, but the real data showed a sudden, sharp cliff.

Why the Hiccup?

The authors realized that the Carbon-13 nucleus is a bit "spinny" (it has a magnetic-like property called spin). Their model treated the interaction as a simple push-and-pull, ignoring the complex "twisting" forces caused by this spin.

• The Analogy: Imagine trying to predict how a spinning top interacts with a wall. If you only look at the speed, you might miss how the spin makes it wobble and bounce differently. The authors suspect that adding a "spin-spin" interaction term to their math would fix the sharp cliff issue, but that's a job for a future paper.

The Takeaway

This paper is a success story for consistency.

• They used one single, realistic map (the Folding Model) to describe both the bouncing protons and the captured protons.
• They didn't need to invent different rules for different situations.
• Even though they missed one tiny detail (the sharp spike in Carbon-13), their approach is much more reliable than the old "guess-and-check" methods.

In short: They built a better, more physics-based map of the atomic world. This map helps us understand how stars cook their fuel, which helps us understand how the universe is made. It's a solid step forward in the quest to decode the secrets of the stars.

Technical Summary

1. Problem Statement

The paper addresses the theoretical description of two critical nuclear astrophysics processes:

• Elastic Scattering: Proton scattering off $^{12}$C and $^{13}$C nuclei at energies around the Coulomb barrier.
• Radiative Capture: The $^{12}$C(p, $\gamma$)$^{13}$N and $^{13}$C(p, $\gamma$)$^{14}$N reactions, which are the first steps in the CNO cycle and crucial for stellar energy production and nucleosynthesis.

The Challenge:
Current models often rely on phenomenological potentials (like Woods-Saxon) or the R-matrix method. While phenomenological potentials are flexible, they introduce many free parameters and lack a microscopic foundation. The R-matrix method fits data well but struggles to predict reactions involving unstable nuclei without experimental input. There is a need for a consistent, microscopic mean-field approach that can describe both the scattering (continuum) and bound states using the same underlying nuclear potential, thereby reducing free parameters and improving predictive power for astrophysical environments.

2. Methodology

The authors employ a microscopic folding model based on the mean-field approximation to generate the proton-nucleus optical potential (OP).

• The Folding Model:

  • The nuclear potential $V(r)$ is calculated by folding a realistic effective nucleon-nucleon (NN) interaction with the ground-state densities of the target nucleus.
  • Interaction: They use the CDM3Y3 interaction, a density-dependent version of the M3Y-Paris interaction. Crucially, this version includes the rearrangement term (RT) derived from the Hugenholtz-van Hove theorem, which is essential for accurately describing elastic scattering at low energies.
  • Densities: Target densities ($^{12}$C and $^{13}$C) are generated using the Independent Particle Model (IPM) with Woods-Saxon single-particle wavefunctions.
  • Components: The total potential includes the nuclear mean-field (central), spin-orbit, and Coulomb terms. The exchange term is treated using the WKB approximation to localize the non-local potential.

• Consistency Strategy:

  • The same folded potential is used for:
    1. Elastic Scattering: Analyzed via the Optical Model (OM).
    2. Radiative Capture: Used to generate both the initial scattering wavefunctions ($\chi$) and the final bound-state wavefunctions ($\phi$) in the potential model calculation of the $(p, \gamma)$ cross-section.
  • Renormalization: To account for higher-order effects (beyond mean-field) and missing interactions (like spin-spin), the folded potentials are renormalized by factors $N_s$ (for scattering/resonance) and $N_b$ (for bound states).

• Reaction Calculation:

  • The radiative capture cross-section is calculated using the Fermi Golden Rule for electric dipole (E1) transitions.
  • The total cross-section is the sum of resonant (via specific excited states) and non-resonant (direct capture) contributions.
  • Spectroscopic factors ($S_F$) and Asymptotic Normalization Coefficients (ANC) are derived or fitted to match experimental data.

3. Key Contributions

1. Unified Framework: The paper demonstrates that a single, microscopic folding model (CDM3Y3 with rearrangement terms) can consistently describe both elastic scattering data and astrophysical S-factors for radiative capture, bridging the gap between scattering and reaction theory.
2. Validation of Density Dependence: It confirms that the specific density-dependent NN interaction (CDM3Y3) with the rearrangement term is necessary to reproduce elastic scattering data at low energies without invoking complex imaginary potentials (as the energy is below inelastic thresholds).
3. Microscopic ANC Prediction: The study provides microscopic calculations of Asymptotic Normalization Coefficients (ANC) for the final states, comparing them favorably with values extracted from transfer reactions and R-matrix analyses.
4. Analysis of Spin-Spin Limitations: The authors explicitly identify the limitations of the current model regarding the $^{13}$C(p, $\gamma$) reaction, attributing discrepancies in reproducing sharp resonance peaks to the neglect of spin-spin interaction terms in the potential.

4. Results

A. Elastic Scattering ($p + ^{12,13}$C)

• The folded potentials successfully describe elastic scattering data at energies of ~1.5–4.6 MeV.
• Renormalization: The real folded potential required renormalization factors of $N_s \approx 1.32\text{--}1.33$ for $^{12}$C and $N_s \approx 1.20\text{--}1.25$ for $^{13}$C to fit the data.
• The imaginary potential was found to be negligible at these sub-barrier energies, validating the use of a purely real microscopic potential.

B. $^{12}$C(p, $\gamma$)$^{13}$N Reaction

• Resonance: The dominant $1/2^+$ resonance at $E \approx 0.422$ MeV was reproduced with $N_s \approx 1.317$.
• Direct Capture: The non-resonant direct capture contribution was found to be negligible compared to the resonance.
• Spectroscopic Factor: A best-fit $S_F \approx 0.40$ was required for the resonance, which is lower than shell-model predictions ($\sim 0.6$), suggesting the potential model effectively spreads the E1 strength.
• Transition Strength: The calculated reduced transition probability $B(E1) \approx 20.31$ e$^2$fm$^2$ agrees well with the experimental value ($22.77 \pm 0.25$ e$^2$fm$^2$).
• ANC: The calculated final state ANC ($A_F^2 \approx 2.06$ fm$^{-1}$) matches well with R-matrix and transfer reaction results.

C. $^{13}$C(p, $\gamma$)$^{14}$N Reaction

• Resonance Challenge: The sharp $1^-$ resonance at $E \approx 0.51$ MeV is difficult to reproduce with standard mean-field potentials.
  • The model required $N_s \approx 1.20$ for the resonance and $N_s \approx 1.08$ for the non-resonant background.
  • The best-fit spectroscopic factor for the resonance was $S_F \approx 0.23$.
• Comparison with Phenomenology: Previous studies using phenomenological Woods-Saxon potentials achieved a fit with $S_F \approx 0.15$ by assuming a mixture of $1^-$ and $0^-$ states. The authors argue this is unrealistic because the $0^-$ state is a higher energy excitation.
• Potential Depth Issue: To reproduce the sharp peak, phenomenological models require unrealistically shallow potentials with extremely small diffuseness ($a \approx 0.04$ fm) or deep Gaussian potentials. The authors' folding model (with realistic diffuseness $a \approx 0.65$ fm) cannot reproduce the sharpness of this peak without including spin-spin interactions, which are currently omitted.
• Excited States: The model successfully described capture to excited states of $^{14}$N ($0^+$ at 2.31 MeV and $1^+$ at 3.95 MeV) with reasonable ANC values.

5. Significance and Conclusion

• Astrophysical Relevance: The study provides a reliable, parameter-reduced method for calculating reaction rates for the CNO cycle, essential for stellar evolution models.
• Model Validity: The work validates the folding model with density-dependent interactions as a robust tool for nuclear astrophysics, capable of describing both scattering and capture consistently.
• Future Directions: The authors conclude that while the mean-field approach is highly successful, the inability to reproduce the sharp $^{13}$C(p, $\gamma$) resonance peak highlights the need to include spin-spin interactions and potentially non-local potentials (like the Perey-Buck form) in future refinements. The consistency of the approach suggests it is a viable path for studying reactions involving unstable nuclei where experimental data is scarce.