Structure of the $^8$B and $^8$Li nuclei and the astrophysical $S_{17}(0)$-factor of the $^7$Be($p,\gamma$)$^8$B direct capture process within a three-body model

Using a three-body potential cluster model with the hyperspherical Lagrange-mesh method, this study calculates the structural properties and asymptotic normalization coefficients of $^8$B and $^8$Li nuclei to derive a precise zero-energy astrophysical $S_{17}(0)$-factor of $22.492 \pm 0.014$ eV b for the $^7$Be($p,\gamma$)$^8$B reaction, revealing that the spin-2 channel dominates the capture process.

The Gist

Imagine the atomic nucleus not as a solid marble, but as a tiny, chaotic dance floor where particles are constantly spinning and holding hands. This paper is a detailed study of two specific dancers on that floor: the Boron-8 (8B) and Lithium-8 (8Li) nuclei.

The authors, working from Uzbekistan, wanted to understand exactly how these nuclei are built and how they behave when they interact with other particles. Here is the breakdown of their work in simple terms.

1. The Setup: A Three-Person Dance

Most people think of a nucleus as a single blob, but the authors treat these specific nuclei as a three-body system.

• The Dancers: They imagine the nucleus as a group of three distinct parts: an Alpha particle (a tight cluster of 2 protons and 2 neutrons), a Helium-3 or Tritium nucleus (a smaller cluster), and a single proton or neutron.
• The Model: They used a mathematical "dance floor" called the hyperspherical Lagrange-mesh method. Think of this as a super-precise 3D grid that allows them to calculate exactly how these three parts move and hold onto each other without crashing into forbidden zones (a concept called the "Pauli exclusion principle," which is like a rule saying two dancers can't occupy the exact same spot at the same time).

2. The Goal: Measuring the "Grip" (ANC)

The main thing the researchers wanted to measure is something called the Asymptotic Normalization Coefficient (ANC).

• The Analogy: Imagine the nucleus is a magnet. The ANC measures how strong the magnetic pull is at the very edge of the magnet, just as a piece of iron is about to snap onto it.
• Why it matters: In the world of stars, nuclei are constantly trying to stick together to create energy. To know how likely they are to stick, you need to know exactly how strong that "edge grip" is. If the grip is too weak, they bounce off; if it's just right, they fuse.

The team calculated this "grip strength" for two different scenarios:

1. Boron-8: How tightly does a proton hold onto a Beryllium-7 core?
2. Lithium-8: How tightly does a neutron hold onto a Lithium-7 core?

They found that the "grip" is different depending on the spin of the particles (like whether the dancers are spinning clockwise or counter-clockwise). They calculated these values with high precision, ensuring their math converged (stopped changing) when they added enough detail to the model.

3. The Big Question: The Solar Thermostat

The ultimate reason for this study is to solve a mystery about the Sun.

• The Reaction: The Sun shines because of a chain reaction where Beryllium-7 grabs a proton to become Boron-8. This step is the "bottleneck" of the process.
• The Problem: We can't easily measure this reaction in a lab because the Sun's core is incredibly hot, but the reaction happens at very low energies where the electric repulsion between particles is like a massive wall.
• The Solution: By calculating the "grip strength" (ANC) perfectly in their model, they could predict the Astrophysical S-factor. Think of the S-factor as a "probability score" for how often this fusion happens.

4. The Results: A New Number for the Sun

The team calculated a specific number for this probability: 22.492 eV b.

Here is how their result compares to the "rulebooks" scientists use:

• Solar Fusion II (The Old Rulebook): Suggested a value around 20.8. The authors' result is a bit higher than this.
• Solar Fusion III (The Newer Rulebook): Suggested a value of 20.5. The authors' result is definitely higher than this.
• The "Best" Solar Model (BAR2M): Interestingly, the most successful modern model of the Sun currently uses a value of 22.4.

The Conclusion: The authors' calculation (22.49) is almost a perfect match for the value used in the most successful current model of the Sun (22.4). This suggests that their way of modeling the three-body dance is very accurate and supports the idea that the Sun's internal temperature and energy production might be slightly different than what the "Solar Fusion III" rulebook suggests.

Summary

In short, the authors built a highly detailed mathematical simulation of how Boron-8 and Lithium-8 nuclei are constructed. By measuring exactly how tightly their outer particles are held, they calculated a specific probability for a nuclear reaction that powers the Sun. Their number matches the most successful modern solar models, suggesting that our current understanding of the Sun's "engine" might need to be tweaked slightly to match their findings.

Technical Summary

1. Problem Statement

The radiative capture reaction $^{7}\text{Be}(p, \gamma)^{8}\text{B}$ is a critical bottleneck in the solar proton-proton (pp) chain, determining the flux of high-energy solar neutrinos. A major challenge in nuclear astrophysics is determining the zero-energy astrophysical S-factor, $S_{17}(0)$, with high precision.

• Experimental Limitations: Direct measurement at stellar energies is impossible due to the Coulomb barrier and extremely low cross-sections.
• Theoretical Discrepancies: Existing theoretical models and indirect experimental extractions (via transfer reactions or breakup) yield conflicting values for the Asymptotic Normalization Coefficients (ANCs) and the resulting $S_{17}(0)$.
  • Solar Fusion II recommends $S_{17}(0) \approx 20.8 \pm 0.7 \text{ (th)} \pm 1.4 \text{ (exp)}$ eV b.
  • Solar Fusion III recommends a slightly lower value of $20.5 \pm 0.70$ eV b.
  • However, the most successful modern Solar Model (BAR2M) relies on a value of 22.4 eV b.
• Goal: The authors aim to resolve these discrepancies by calculating the structure of $^{8}\text{B}$ and its mirror nucleus $^{8}\text{Li}$ using a self-consistent three-body model to derive precise ANC values and the corresponding $S_{17}(0)$.

2. Methodology

The authors employ a three-body potential cluster model treating the nuclei as:

• $^{8}\text{B}$: $\alpha + ^{3}\text{He} + p$
• $^{8}\text{Li}$: $\alpha + ^{3}\text{H} + n$

Key Technical Components:

• Hyperspherical Lagrange-Mesh Method: The Schrödinger equation is solved in hyperspherical coordinates. The wave function is expanded over a basis of hyperspherical harmonics.
  • Convergence was achieved with a maximum hypermomentum of $K_{\text{max}} = 22$ for ground states ($2^+$) and $K_{\text{max}} = 28$ for excited states ($1^+$).
• Interaction Potentials: Realistic two-body potentials were adopted from literature:
  • $\alpha$-Nucleon ($\alpha$-p, $\alpha$-n): Voronchev et al. potential (including spin-orbit and Pauli forbidden states).
  • $\alpha$-$^{3}\text{He}$ / $\alpha$-$^{3}\text{H}$: Gaussian potentials with Pauli forbidden states.
  • $p$-$^{3}\text{He}$ / $n$-$^{3}\text{H}$: Gaussian attractive + exponential repulsive potentials.
  • Note: Orthogonalizing Pseudopotentials (OPP) were used to eliminate Pauli-forbidden states.
• ANC Extraction: The ANC values for the virtual decays ($^{8}\text{B} \to ^{7}\text{Be} + p$ and $^{8}\text{Li} \to ^{7}\text{Li} + n$) were determined by matching the overlap integral of the three-body wave function with the asymptotic Whittaker function at large inter-cluster distances.
• Astrophysical Factor Calculation: The $S_{17}(0)$ factor was calculated using the asymptotic theory developed by D. Baye, which relates the S-factor to the S-wave scattering length and the squared ANC ($C^2$).

3. Key Contributions and Results

A. Nuclear Structure Properties

The model successfully reproduced experimental binding energies without requiring three-body forces:

• Binding Energies:
  • $^{8}\text{B}$ ($2^+$): Calculated $-1.718$ MeV vs. Exp. $-1.724$ MeV.
  • $^{8}\text{Li}$ ($2^+$): Calculated $-4.340$ MeV vs. Exp. $-4.499$ MeV.
• Matter Radii:
  • $^{8}\text{B}$: Calculated $2.76$ fm (Exp. $2.58 \pm 0.06$ fm). The model slightly overestimates the radius by ~5%, consistent with its halo nature.
  • $^{8}\text{Li}$: Calculated $2.56$ fm (Exp. $2.50 \pm 0.06$ fm), well within experimental error.

B. Asymptotic Normalization Coefficients (ANCs)

The study provided self-consistent ANC values for spin channels $S=1$ and $S=2$:

• $^{8}\text{B} \to ^{7}\text{Be} + p$:
  • $C_{S=1} = 0.211 \text{ fm}^{-1/2}$
  • $C_{S=2} = 0.739 \text{ fm}^{-1/2}$
  • Total Squared ANC: $C^2 = 0.591 \text{ fm}^{-1}$.
• $^{8}\text{Li} \to ^{7}\text{Li} + n$:
  • $C_{S=1} = 0.220 \text{ fm}^{-1/2}$
  • $C_{S=2} = 0.774 \text{ fm}^{-1/2}$
  • Total Squared ANC: $C^2 = 0.648 \text{ fm}^{-1}$.

Comparison: The calculated $C^2$ for $^{8}\text{B}$ ($0.591$) is higher than the ab-initio NCSM result ($0.509$) but aligns well with values extracted from the $^{7}\text{Be}(d,n)^{8}\text{B}$ reaction ($0.613 \pm 0.060$).

C. Astrophysical S-Factor ($S_{17}(0)$)

Using the derived ANCs and experimental scattering lengths, the zero-energy astrophysical factor was calculated:

• Spin 2 Channel (Dominant): $S^{(2)}_{17}(0) = 20.838 \pm 0.014$ eV b.
• Spin 1 Channel: $S^{(1)}_{17}(0) = 1.654 \pm 0.003$ eV b.
• Total Result: $S_{17}(0) = 22.492 \pm 0.014$ eV b.

4. Significance and Conclusion

• Resolution of Solar Model Discrepancy: The calculated value of 22.492 eV b is in excellent agreement with the value (22.4 eV b) used in the BAR2M solar model, which is currently considered the most successful model for solar neutrino predictions.
• Comparison with Reviews:
  • The result is consistent with the Solar Fusion II estimate ($20.8 \pm \dots$).
  • It is slightly higher than the Solar Fusion III recommendation ($20.5 \pm 0.70$ eV b), suggesting that the lower value in SFIII might underestimate the contribution of the spin-2 channel or the specific ANC magnitude derived from this rigorous three-body approach.
• Methodological Validation: The study demonstrates that the hyperspherical Lagrange-mesh method, when combined with realistic two-body potentials, provides a robust and self-consistent framework for describing light halo nuclei and extracting astrophysical factors without adjustable three-body forces.
• Implication for Neutrino Physics: The findings support the higher $S_{17}(0)$ value, which has direct implications for the predicted flux of $^{8}\text{B}$ solar neutrinos and the interpretation of solar neutrino data.

In summary, this paper provides a rigorous theoretical confirmation that the astrophysical S-factor for $^{7}\text{Be}(p, \gamma)^{8}\text{B}$ is likely closer to 22.5 eV b, bridging the gap between theoretical nuclear structure models and the requirements of modern solar modeling.