5Li Resonance in α–p Elastic Scattering

• The paper demonstrates that the 5Li resonance elevates α–p elastic scattering cross sections by up to three orders of magnitude at resonant energies.
• Using advanced techniques like the ERF and Δ-method, the study accurately extracts resonance parameters and validates them against ab initio nuclear models.
• The findings lead to refined Geant4 simulations, thereby improving estimates of background events in neutrino detectors and enhancing precision measurements.

The $^5$Li resonance in $\alpha$–$p$ elastic scattering is a key nuclear phenomenon underlying a strongly enhanced cross section when an $\alpha$ particle scatters from a proton at center-of-mass energies coinciding with the formation of the unstable $^5$Li nucleus. This resonance, which manifests as sharp peaks in the elastic cross section, plays a critical role in diverse areas of nuclear structure, nuclear reaction modeling, and precision neutrino measurements, notably impacting the backgrounds in large-scale liquid scintillator neutrino detectors.

1. Physical Origin and Resonant Character

The $^5$Li resonance arises during $\alpha$–$p$ elastic scattering when the center-of-mass energy matches the excitation energy of the unbound $^5$Li nucleus. This system features a ground-state resonance with spin-parity $J^\pi=3/2^-$ and a first excited resonant state with $J^\pi=1/2^-$. The ground-state resonance is located approximately at $E_r \approx 1.4$–$1.7$ MeV above the $p+\alpha$ threshold and exhibits a total width in the range $\Gamma \approx 1.2$–$1.3$ MeV, while the $1/2^-$ resonance lies higher and is much broader ($E_r \approx 2.6$–$3.2$ MeV, $\Gamma \approx 4.5$–$6.6$ MeV) (Shirokov et al., 2018, Orlov, 2020). These resonant states are extremely short-lived (lifetime $\sim 4\times10^{-22}$ s for the $3/2^-$ state (Wang et al., 18 Jan 2026)), decaying back to $\alpha$ + $p$. At $\alpha$ energies above 5 MeV, the presence of these resonances elevates the elastic cross section by two to three orders of magnitude compared to the non-resonant Coulomb (Rutherford) scattering value (Wang et al., 18 Jan 2026).

2. Theoretical Formalism: Scattering, Cross Sections, and Resonance Extraction

The resonant scattering is described using the single-level, single-channel Breit–Wigner formula for the angle-integrated cross section,

$$\sigma_\mathrm{res}(E) = \frac{\pi}{k^2} \frac{2J+1}{(2s_\alpha+1)(2s_p+1)} \frac{\Gamma_\alpha(E)\,\Gamma_p(E)}{(E-E_R)^2 + (\Gamma/2)^2}$$

where $k$ is the center-of-mass momentum, $s_\alpha$ and $s_p$ are the $\alpha$ and proton spins, and $J$ is the resonance spin ($3/2$ for the ground state). The resonance energy $E_R$ and the partial/total widths $\Gamma_\alpha$, $\Gamma_p$, $\Gamma$ define the peak position and shape ($E_R \approx 8.8$ MeV in the $\alpha$–$p$ center-of-mass, corresponding to $1.4$–$1.7$ MeV above threshold depending on convention) (Wang et al., 18 Jan 2026, Shirokov et al., 2018, Orlov, 2020).

Background, non-resonant scattering is governed by the Rutherford formula. Differential cross sections combine Coulomb and resonant nuclear amplitudes, requiring R-matrix or effective-range analyses for accurate angular distributions. The strongly forward-peaked nature of the resonance is confirmed by both experimental data and theoretical fits, with differential cross sections at small angles ($\theta<20^\circ$) being enhanced by one to two orders of magnitude above the Coulomb expectation at resonance energies (Wang et al., 18 Jan 2026).

The extraction of resonance energies and widths from experimental phase-shift data can be accomplished via several methods:

• Effective-Range Function (ERF) Method: Fits the Coulomb-modified effective-range function $K_l(E)$ to measured phase shifts.
• “Delta” ($\Delta$) Method: Isolates the nuclear contribution ($\Delta_l(E)$) from the ERF, providing a more accurate determination of pole positions for charged systems (Orlov, 2020).

The $\Delta$-method has been shown to improve the self-consistency, particularly for large-charge systems, yielding more accurate $E_r$ and $\Gamma$ values and adjusted asymptotic normalization coefficients (ANCs).

3. Ab Initio and Many-Body Calculations

The structure and resonance properties of $^5$Li have been studied using the Single-State Harmonic Oscillator Representation of Scattering Equations (SS-HORSE) built on No-Core Shell Model (NCSM) calculations. By matching variational eigenstates to asymptotic forms in an oscillator basis and extracting complex-energy poles from phase shifts, realistic nuclear interactions such as JISP16 and Daejeon16 have reproduced the characteristic $3/2^-$ and $1/2^-$ resonance structure (Shirokov et al., 2018).

The SS-HORSE approach yields resonance parameters in good agreement with experimental data and with dispersion from different ab initio interactions. For example, Daejeon16 gives $E_{3/2^-}=1.52$ MeV ($\Gamma=1.05$ MeV), while JISP16 gives $1.84$ MeV ($1.80$ MeV), compared to experimental values of $1.69$ MeV ($1.23$ MeV) (Shirokov et al., 2018). The spin-orbit splitting $\Delta = E_r(1/2^-) - E_r(3/2^-)$ is reproduced at the 0.2 MeV level.

The s-wave ($\ell=0$) channel provides a slowly varying, non-resonant background: the effective-range expansion gives a scattering length $a_0 \sim -2$ fm and range $r_0 \sim 2$ fm, modestly enhancing the low-energy cross section (Shirokov et al., 2018).

4. Detector Simulations, Geant4 Modeling, and Experimental Relevance

In large liquid scintillator detectors (e.g., SNO+, JUNO, Daya Bay), cascade decays of $^{214}$Bi–$^{214}$Po produce $\alpha$ particles with energies overlapping the critical $^5$Li resonance region. Default Geant4 simulations, omitting the resonance, severely underestimate high-energy proton recoils and thus background rates from $\alpha$–$p$ events. Modifying Geant4 to include the full R-matrix-derived $\mathrm{d}\sigma/\mathrm{d}\Omega(E, \theta)$ reproduces the high-energy tail in deposited energy spectra, increasing the calculated misidentification probability by up to a factor of 10—from $1.0 \times 10^{-5}$ to $1.9 \times 10^{-4}$ for $7.68$ MeV $\alpha$ depositions in the window $2$–$2.5$ MeV (Wang et al., 18 Jan 2026).

This improved modeling reconciles simulated and observed background rates, particularly in the context of inverse beta decay (IBD) backgrounds in neutrino oscillation and geoneutrino measurements.

5. Resonance Parameters from Phase-Shift Analysis and Systematic Comparisons

The two dominant resonances in $^5$Li ($J^\pi = 3/2^-$ and $1/2^-$) possess well-established widths, resonance energies, and ANCs. Comparative studies using the ERF and Delta methods illustrate the following for the $3/2^-$ ground-state resonance:

• $E_r = 1.390 \pm 0.010$ MeV, $\Gamma = 1.301 \pm 0.020$ MeV, $|C_{3/2}| = 0.325 \pm 0.010$ fm$^{-1/2}$ ($\Delta$ method, preferred)
• For the $1/2^-$ excited-state: $E_r = 2.611 \pm 0.020$ MeV, $\Gamma = 4.534 \pm 0.030$ MeV, $|C_{1/2}| = 0.505 \pm 0.015$ fm$^{-1/2}$ (Orlov, 2020)

Width and energy estimates differ by up to 15–20% depending on method, with the $\Delta$ method generally yielding more robust extraction especially when significant Coulomb background is present.

State	$E_r$ (MeV, $\Delta$)	$\Gamma$ (MeV, $\Delta$)	$|C|$ (fm$^{-1/2}$, $\Delta$)
$3/2^-$	1.390	1.301	0.325
$1/2^-$	2.611	4.534	0.505

6. Impact on Precision Neutrino Measurements and Astrophysical Context

The $^5$Li resonance critically affects background estimates in measurements of neutrino oscillation parameters ($\theta_{12}$, $\theta_{13}$, $\Delta m^2_{21}$, $|\Delta m^2_{31}|$). In particular, the $\alpha$–$p$ resonance-induced enhancement causes correlated $^{214}$Bi–$^{214}$Po backgrounds to mimic IBD events, impacting the extracted geoneutrino and reactor antineutrino fluxes. For Daya Bay, the correct inclusion of the $^5$Li resonance increases the measured $\sin^2 2\theta_{13}$ by approximately $+0.012$, and the Particle Data Group's reported $\sin^2 \theta_{13}$ by $\sim0.006$ (1$\sigma$) (Wang et al., 18 Jan 2026). For geoneutrino fluxes, corrections to the low-energy tail shape alter the inferred $U/Th$ content and their ratio by up to $\mathcal{O}(1\%)$.

A plausible implication is that any future refinements in nuclear modeling or detector calibration must robustly account for $\alpha$–$p$ resonance cross sections above 5 MeV, as neglecting these effects may bias key neutrino and astrophysics observables.

7. Methodological Developments and Systematic Considerations

• The $\Delta$-method provides improved accuracy for extracting resonance parameters in charged systems, with the Coulomb background isolated and removed to yield more reliable nuclear resonance properties (Orlov, 2020).
• Ab initio NCSM/SS-HORSE calculations, with both JISP16 and Daejeon16 interactions, offer convergent, model-independent resonance characterizations and corroborate phase-shift and ANCs extracted from experimental data (Shirokov et al., 2018).
• Modifications to Monte Carlo transport codes like Geant4, informed by R-matrix theory, are required to simulate nuclear backgrounds in neutrino detectors at percent accuracy (Wang et al., 18 Jan 2026).

Accurate parameterization of $\alpha$–$p$ scattering including $^5$Li resonances is thus essential for nuclear reaction modeling, precise particle detector simulations, and correct interpretation of next-generation neutrino experiments.