Low-energy $^7$Li($n,γ$)$^8$Li and $^7$Be($p,γ$)$^8$B radiative capture reactions within the Skyrme Hartree-Fock approach

This study employs the Skyrme Hartree-Fock potential model to simultaneously analyze low-energy radiative capture reactions in $^7$Be($p,γ$)$^8$B and $^7$Li($n,γ$)$^8$Li, successfully describing electric dipole transitions with minimal adjustment and determining the astrophysical $S_{17}(0)$ factor for the former reaction to be 22.3 eV b.

The Gist

Imagine the universe as a giant, cosmic kitchen where stars are the chefs. To cook the elements that make up everything around us (like the carbon in your DNA or the iron in your blood), these stellar chefs need to fuse tiny particles together. Two specific recipes are crucial for this cosmic cooking:

1. The Solar Recipe: A Lithium-7 nucleus catches a proton (a hydrogen nucleus) to become Boron-8. This happens in our Sun and creates high-energy neutrinos (ghostly particles that stream through us).
2. The Big Bang Recipe: A Lithium-7 nucleus catches a neutron to become Lithium-8. This happened in the early universe and helped form heavier elements.

The problem is that these reactions happen at incredibly low energies, making them very hard to measure in a lab. It's like trying to catch a specific, shy butterfly that only lands when the wind is perfectly calm. Because we can't measure them perfectly, there is a lot of uncertainty in our "cookbook" for how stars work and how the universe began.

The Paper's Mission: A New Way to Look at the Butterfly

This paper, written by physicists Nguyen Le Anh and Bui Minh Loc, tries to solve this uncertainty using a clever mathematical tool called the Skyrme Hartree-Fock approach.

Here is how they did it, using some everyday analogies:

1. The "Microscope" vs. The "Map"

Most scientists look at these reactions using "phenomenological models." Think of this like drawing a map based on where people say they walked. You adjust the lines until they match the story.

These authors, however, used a microscopic approach. Imagine instead of a map, they built a detailed 3D simulation of the terrain itself. They calculated exactly how the protons and neutrons behave inside the nucleus, like simulating the movement of every single grain of sand on a beach rather than just guessing where the tide goes.

2. The "Double-Check" Strategy

The two reactions (Lithium catching a proton vs. Lithium catching a neutron) are "twins" in the world of physics. They are so similar that if you understand one, you understand the other.

• The Neutron Twin: Lithium-7 + Neutron. This is easier to study because neutrons don't have an electric charge, so they don't get repelled by the nucleus.
• The Proton Twin: Lithium-7 + Proton. This is the hard one because protons are positively charged and repel each other (like trying to push two strong magnets together north-to-north).

The authors used their "microscope" to study the easy twin (the neutron reaction) first. Once they got the settings right there, they applied the exact same settings to the hard twin (the proton reaction). This cross-checking gives them much more confidence in their results.

3. Tuning the Radio

In their simulation, they had to adjust a few "knobs" (parameters) to make the math match reality.

• The "Depth" Knob: They adjusted how deep the "well" is that holds the particles together. They tuned this so the energy levels matched what we see in experiments.
• The "Spectroscopic Factor" Knob: This is a measure of how "pure" the particle arrangement is. Think of it like a recipe. If a cake recipe calls for 1 cup of flour, but your cake is only 80% flour and 20% mystery ingredients, you need a factor to correct the math. They adjusted this factor until their simulation perfectly matched the experimental data from previous years.

The Big Result: A Precise Number

The ultimate goal of this paper was to find a specific number called S17(0).

• What is it? It's the "efficiency rating" of the Sun's proton-catching reaction at zero energy.
• Why does it matter? This number tells astrophysicists exactly how much high-energy neutrino the Sun produces. If this number is wrong, our models of how the Sun shines and how old it is might be wrong.

The Verdict:
After running their complex simulations, the authors calculated the S17(0) factor to be 22.3 eV b.

This is a very precise number that sits comfortably among other recent measurements. It suggests that the "Skyrme Hartree-Fock" method (their microscopic simulation) is a very powerful tool. It proves that you don't need to guess the rules of the game; you can calculate them from the ground up and get a result that matches the real world.

In Summary

The authors built a high-tech, microscopic simulation of atomic nuclei. They used the easier version of a nuclear reaction (catching a neutron) to calibrate their machine, then used that same machine to predict the harder version (catching a proton). Their result gives us a clearer, more confident picture of how our Sun generates energy and how the universe created its elements.

Technical Summary

1. Problem Statement

The paper addresses the theoretical description of low-energy radiative capture reactions, specifically $^7\text{Li}(n, \gamma)^8\text{Li}$ and $^7\text{Be}(p, \gamma)^8\text{B}$. These reactions are critical in nuclear astrophysics:

• $^7\text{Be}(p, \gamma)^8\text{B}$: Crucial for the production of high-energy solar neutrinos and the nucleosynthesis of low-mass stars. Determining its cross-section at zero energy (the astrophysical $S_{17}(0)$ factor) is challenging because measurements must be extrapolated from extremely low energies where the cross-section is vanishingly small.
• $^7\text{Li}(n, \gamma)^8\text{Li}$: Important for the nucleosynthesis of heavy nuclei in inhomogeneous Big Bang models.

Due to isospin symmetry, these two reactions are often studied simultaneously. However, existing theoretical approaches (phenomenological potentials, shell models, $R$-matrix, and ab initio calculations) often struggle to consistently describe both bound and scattering states or require significant phenomenological adjustments. The authors aim to provide a microscopic, consistent approach using the Skyrme Hartree-Fock (SHF) formalism to calculate both bound and scattering wave functions simultaneously, minimizing free parameters while accurately reproducing experimental data.

2. Methodology

The authors employ a potential model where the radiative capture is treated as an electromagnetic transition between a continuum (scattering) state and a bound state. The core innovation lies in how the wave functions for these states are generated.

• Skyrme Hartree-Fock (SHF) Framework:

  • Instead of using phenomenological Woods-Saxon potentials, the authors use the SHF calculation (using the SLy4 interaction) to generate the single-particle wave functions.
  • Simultaneous Calculation: Both the bound state wave function ($\phi_{n\ell'j'}$) and the scattering state wave function ($\chi_{\ell j}$) are derived from the same microscopic mean-field potential.
  • Equations: The radial Hartree-Fock equations are solved, incorporating central, spin-orbit, and Coulomb potentials. The effective mass $m^*(r)$ is included to account for non-locality.
  • Corrections:
    • Bound States: The SHF approximation typically yields tightly bound states. The authors adjust a parameter $N_b$ (well-depth scaling) to ensure the single-particle energy matches the experimental nucleon separation energy ($Q$-value).
    • Scattering States: The mean-field approach does not naturally reproduce low-lying resonances. A parameter $N_s$ is adjusted to shift the scattering phase shifts and reproduce the exact positions of observed resonances.

• Transition Operators:

  • The study considers both Electric Dipole (E1) and Magnetic Dipole (M1) transitions.
  • E1 Transitions: Dominant for non-resonant capture. The cross-section depends on the radial overlap integral of the scattering and bound states.
  • M1 Transitions: Responsible for low-lying resonances (e.g., at 633 keV and 2184 keV for $^7\text{Be}+p$). The M1 operator includes contributions from orbital motion, spin, and the core magnetic moment.

• Spectroscopic Factors ($S_F$):

  • Since the SHF model is a mean-field approach, it lacks complex many-body correlations. A spectroscopic factor $S_F$ is introduced to account for missing configurations.
  • Strategy: $N_b$ and $N_s$ are constrained by experimental separation energies and resonance positions. $S_F$ is then fine-tuned to reproduce experimental cross-section data at specific energies (hundreds of keV).

3. Key Contributions

1. Microscopic Consistency: The work demonstrates that a single microscopic framework (Skyrme Hartree-Fock) can consistently describe both bound and scattering states for light nuclei ($A=7, 8$), which are often considered too light for mean-field approximations.
2. Simultaneous Analysis: The study successfully analyzes both $^7\text{Li}(n, \gamma)^8\text{Li}$ and $^7\text{Be}(p, \gamma)^8\text{B}$ side-by-side, leveraging isospin symmetry to validate the approach.
3. Resonance Handling: The authors explicitly model low-lying resonant M1 transitions (at 633 keV, 2184 keV, and 222 keV) by adjusting the scattering potential depth ($N_s$) to match resonance energies, rather than treating them purely phenomenologically.
4. Parameter Economy: The model requires minimal adjustment. For non-resonant E1 transitions, only the bound state depth ($N_b$) is adjusted. The spectroscopic factor is the primary parameter fitted to data, but its value is physically constrained by the resonance analysis.

4. Results

• Non-Resonant E1 Transitions:

  • The calculated E1 cross-sections for both reactions reproduce experimental data well.
  • For $^7\text{Li}(n, \gamma)^8\text{Li}$, the calculation matches data from Refs. [8, 28–30] with $N_b = 0.72$ and $S_F = 1.0$.
  • For $^7\text{Be}(p, \gamma)^8\text{B}$, the calculation matches data from Refs. [31–39] with $N_b = 0.70$ and $S_F = 0.8$.
  • The $S_{17}(0)$ value is found to be 22.3 eV b when $S_F = 0.8$ is used. This value is relatively insensitive to the specific choice of $N_b$ (e.g., changing $N_b$ from 0.70 to 1.00 only shifts $S_{17}(0)$ to 21.2 eV b, though it requires an unphysical $S_F$).

• Resonant M1 Transitions:

  • $^7\text{Be}(p, \gamma)^8\text{B}$:
    • The resonance at 633 keV (p-wave capture to $1p_{3/2}$) is reproduced with $N_s = 0.67$ and $S_F = 1.0$.
    • The resonance at 2184 keV (dominated by d-wave capture to $1d_{5/2}$) is reproduced with $N_s = 1.27$ and a small $S_F = 0.1$, indicating significant interference from other configurations at this energy.
  • $^7\text{Li}(n, \gamma)^8\text{Li}$:
    • The resonance at 222 keV is reproduced. The authors note a discrepancy in experimental data between Refs. [28] and [29]. To match Ref. [29], they use $S_F = 0.05$, which is consistent with the small $S_F$ found for the 2184 keV resonance in the mirror nucleus.

• Astrophysical S-Factor:

  • The final derived value for the $S_{17}(0)$ factor is 22.3 eV b. This falls within the range of recent experimental values (e.g., 20.6–21.5 eV b) and theoretical predictions (19.1–20.5 eV b), providing a robust theoretical constraint.

5. Significance

• Validation of Mean-Field for Light Nuclei: The paper challenges the notion that mean-field models like Skyrme Hartree-Fock are unsuitable for light nuclei ($A=7, 8$). It shows that with appropriate corrections for asymptotic behavior and resonance positions, SHF provides a reliable microscopic basis for radiative capture.
• Constraint on Solar Neutrinos: By providing a consistent theoretical determination of $S_{17}(0)$, the work contributes to reducing uncertainties in solar neutrino flux predictions, which are vital for understanding solar physics and neutrino oscillations.
• Methodological Framework: The approach of using SHF to generate wave functions and then applying a potential model for the transition matrix elements offers a computationally efficient alternative to complex ab initio or multi-cluster calculations, while maintaining a high degree of physical consistency between bound and continuum states.
• Isospin Symmetry: The successful simultaneous description of the neutron and proton capture reactions reinforces the utility of isospin symmetry in constraining nuclear potentials and spectroscopic factors.