Low-energy $^{3}$He($α,γ$)$^{7}$Be reaction within the Skyrme potential framework

This paper employs a microscopic Skyrme Hartree-Fock-based folding potential model to simultaneously describe low-energy elastic scattering and calculate the astrophysical $S$ factor for the $^{3}$He($\alpha,\gamma$)$^{7}$Be reaction, yielding a recommended zero-energy value of $S_{34}(0) = 0.610 \pm 0.024$ keV b that agrees well with experimental data.

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 helium in balloons or the oxygen in our lungs), these chefs need to mix ingredients together. One of the most important recipes in the "solar kitchen" is a reaction where a tiny particle called Helium-3 (a light version of helium) bumps into a regular Helium-4 nucleus to create Beryllium-7.

This recipe is tricky. It's like trying to push two strong magnets together when their same poles are facing each other—they naturally repel. In physics, this is called the Coulomb barrier. Because they repel so strongly, this reaction happens very rarely, and when it does, it releases very little energy. This makes it incredibly hard for scientists to measure in a lab, and even harder to predict with math.

Here is what this paper does, explained simply:

1. The Problem: A Missing Recipe Card

Scientists know this reaction is crucial for how the Sun shines and how the early universe formed its first elements. However, because the reaction is so rare, the "recipe cards" (experimental data) we have are fuzzy and sometimes contradict each other. We need a better way to predict exactly how often this happens, especially at the low energies found inside stars.

2. The Solution: Building a Better Map

The authors decided to stop guessing and start building a microscopic map of the forces at play. Instead of treating the particles as simple, featureless balls, they used a sophisticated mathematical framework called the Skyrme Hartree-Fock model.

Think of this like this:

• Old Way: Trying to predict how two cars will crash by just looking at their speed and weight.
• This Paper's Way: Looking at the engine, the tires, the suspension, and the driver's behavior to understand exactly how they interact.

They started by mapping out how a single proton interacts with a Helium-4 nucleus. They found that their "map" was slightly off, so they applied a tiny "zoom" or scaling factor (like adjusting the focus on a camera lens) to make the map match real-world observations perfectly.

3. The Folding Technique: Making a Sandwich

Once they had the perfect map for a single proton, they needed to figure out how a whole Helium-3 nucleus (which is made of 2 protons and 1 neutron) interacts with Helium-4.

They used a technique called folding. Imagine you have a spread of peanut butter (the force field of Helium-4) and you want to know how a slice of bread (Helium-3) sits on top of it. You don't just look at the center of the bread; you look at how every crumb of the bread interacts with the peanut butter. By "folding" the single-proton map over the shape of the Helium-3 nucleus, they created a complete interaction map for the whole system.

They tested two different "shapes" for the Helium-3 bread (based on different scientific theories about how its particles are arranged) and found that one shape fit the data much better than the other.

4. The Result: A Clearer Picture of the Sun

Using this new, highly detailed map, they calculated the Astrophysical S-factor. Think of the S-factor as the "efficiency rating" of the recipe. How many times does the reaction happen per second inside the Sun?

• The Outcome: Their calculation matched existing experimental data very well.
• The Recommendation: They provided a new, highly precise number for this efficiency: 0.610. This number helps astronomers refine their models of how the Sun produces energy and how much neutrino radiation (ghostly particles) it sends to Earth.

Why This Matters

This paper is like upgrading from a hand-drawn sketch of a city to a high-definition GPS.

• For the Sun: It helps us understand exactly how much energy our star is generating and why it shines the way it does.
• For the Universe: It helps us understand how the first elements were forged in the Big Bang.
• For Science: It proves that we can use fundamental laws of physics (nuclear forces) to predict complex cosmic events without needing to rely solely on guesswork or messy experiments.

In short, the authors built a better microscope to look at the smallest building blocks of the universe, allowing us to see the "recipe" of the stars with much greater clarity.

Technical Summary

1. Problem Statement

The $^3$He($\alpha, \gamma$)$^7$Be reaction is a critical process in the proton-proton (pp) chain of stellar hydrogen burning and Big Bang nucleosynthesis (BBN). Its reaction rate directly influences solar neutrino fluxes and primordial element abundances.

• Experimental Challenge: At astrophysical energies (center-of-mass energy $E < 500$ keV), the reaction cross-section is extremely low due to the strong Coulomb barrier. This leads to significant uncertainties in experimental data, with existing measurements showing inconsistencies.
• Theoretical Challenge: While various models (R-matrix, cluster models, shell models) exist, there is a need for a unified, microscopic framework that can simultaneously describe elastic scattering and radiative capture with minimal phenomenological adjustments, thereby improving predictive power for astrophysical extrapolations.

2. Methodology

The authors employ a microscopic potential-model approach derived from the Skyrme Hartree-Fock (HF) framework, extended to the continuum. The methodology proceeds in three main stages:

A. Construction of the Nucleon-Target Potential

• Base Calculation: Self-consistent Skyrme HF calculations are performed using the SLy4 interaction to generate single-particle states and the mean-field potential for the $^4$He nucleus.
• Continuum Extension: The HF potential is extended to the continuum to describe scattering states.
• Scaling: A central scaling factor ($\lambda_0$) is introduced to slightly adjust the central potential to reproduce experimental $p + \alpha$ elastic scattering phase shifts. The spin-orbit strength is also scaled ($\lambda_{ls}$) where necessary.

B. Construction of the $^3$He + $\alpha$ Potential (Folding Model)

• Folding Procedure: The nucleon-$\alpha$ potential is folded with the density distribution of the $^3$He projectile to construct the nucleus-nucleus potential:
  $$V_{^3\text{He}+\alpha}(E, r) = \sum_{q=n,p} \int \rho_q^{(^3\text{He})}(\vec{r}_a) V_q(E, s) d^3r_a$$
• Density Distributions: Two different two-parameter Fermi (2pF) density distributions for $^3$He are tested to assess sensitivity:
  1. Schechter and Canto (SC) density: RMS radius $r_{rms} = 2.316$ fm.
  2. Ngô and Ngô (Ngo) density: RMS radius $r_{rms} = 2.182$ fm (closer to the experimental value of 2.023 fm).

C. Radiative Capture Calculation

• Transition Mechanism: The calculation focuses on the Electric Dipole (E1) transition, which dominates at low astrophysical energies.
• Wave Functions: Radial wave functions for scattering and bound states are obtained by solving the Schrödinger equation using the folded HF potentials.
• Cross Section: The astrophysical S factor, $S(E)$, is calculated using the E1 transition matrix elements, which depend on the overlap integral between scattering and bound state wave functions.

3. Key Contributions

• Unified Framework: This is the first application of the Skyrme HF-based folding model to describe both low-energy nucleus-nucleus elastic scattering ($^3$He + $\alpha$) and radiative capture ($^3$He($\alpha, \gamma$)$^7$Be) within a single microscopic framework.
• Parameter Reduction: The model minimizes free parameters. The scaling factors are constrained primarily by elastic scattering phase shifts, rather than fitting the capture data directly.
• Density Sensitivity Analysis: The study systematically investigates the impact of the projectile density distribution (SC vs. Ngo) on the calculated S factor and Asymptotic Normalization Coefficients (ANCs).

4. Results

Elastic Scattering ($p + \alpha$ and $^3$He + $\alpha$)

• $p + \alpha$: A scaling factor of $\lambda_0 = 0.88$ successfully reproduces experimental $s_{1/2}$ phase shifts without modifying the spin-orbit term.
• $^3$He + $\alpha$:
  • The folded potentials reproduce experimental $s$-wave phase shifts well.
  • Optimal scaling factors were found to be $\lambda_0 \approx 0.85$ (Ngo) and $0.88$ (SC) for scattering states.
  • The spin-orbit strength required for $^3$He + $\alpha$ ($\lambda_{ls} = 0.14$) is roughly an order of magnitude weaker than in the $p + \alpha$ case, suggesting a specific structural feature of the $^3$He cluster.
  • Angular distributions calculated with these scaled potentials agree well with experimental data up to 3.0 MeV.

Radiative Capture and Astrophysical S Factor

• Binding Energies: To reproduce the binding energies of the $^7$Be ground state ($3/2^-$) and first excited state ($1/2^-$), different scaling factors were required for bound states ($\lambda_0 \approx 1.04 - 1.07$) compared to scattering states.
• Asymptotic Normalization Coefficients (ANCs): The ANCs extracted using the Ngo density showed better agreement with empirical values (within 3% of Tursunmahatov et al.) compared to the SC density.
• S Factor ($S_{34}$):
  • The calculated S factor shows good agreement with experimental data from Refs. [5–13].
  • The Ngo density yields a better overall fit to the data, attributed to its more realistic spatial distribution (smaller RMS radius).
  • The branching ratio between capture to the ground state and the first excited state is calculated as $R = 0.43$, consistent with previous findings and insensitive to the density choice.

Recommended Value

The extrapolated zero-energy S factor ($S_{34}(0)$) values are:

• SC Density: $0.666 \pm 0.026$ keV b
• Ngo Density: $0.610 \pm 0.024$ keV b

The authors recommend the Ngo density result ($0.610 \pm 0.024$ keV b) as the final value. This aligns excellently with previous evaluations and sits near the upper limit of the recent Solar Fusion III evaluation.

5. Significance

• Predictive Power: The study demonstrates that a microscopic potential derived from the Skyrme HF framework, with minimal scaling, can accurately predict low-energy nuclear reactions without relying on purely phenomenological fits to capture data.
• Astrophysical Impact: The recommended $S_{34}(0)$ value provides a robust input for solar models and BBN calculations, helping to resolve discrepancies in solar neutrino flux predictions.
• Methodological Validation: The consistency between the scaling factors derived from scattering data and those required for capture calculations validates the unified nature of the approach, suggesting that the underlying nuclear interaction is correctly modeled.
• Density Dependence: The work highlights the importance of accurate projectile density distributions in folding models, showing that the Ngo density provides a superior description for the $^3$He + $\alpha$ system compared to the SC density.