A possible solution to the gallium anomaly moving beyond the leptonic wave function factorization

Technical Summary: A Possible Solution to the Gallium Anomaly Moving Beyond the Leptonic Wave Function Factorization

The Problem
For over three decades, a persistent discrepancy known as the "gallium anomaly" has existed between measured and predicted neutrino capture rates on $^{71}\text{Ga}$. Experiments utilizing radioactive sources ($^{51}\text{Cr}$ and $^{37}\text{Ar}$), specifically GALLEX, SAGE, and more recently BEST, have observed a $\sim20\%$ deficit in the electron neutrino capture rate ($\nu_e + ^{71}\text{Ga} \to ^{71}\text{Ge} + e^-$). This deficit now exceeds $5\sigma$ significance. While this anomaly has prompted speculation regarding new physics, such as short-baseline active-sterile neutrino oscillations, these scenarios face significant tension with results from reactor antineutrino experiments, solar neutrino bounds, MicroBooNE, and KATRIN. Consequently, a rigorous reassessment of the theoretical assumptions underlying the Standard Model prediction for the inverse beta decay (IBD) cross-section is required.

Methodology
The authors critically re-examine the standard theoretical treatment of the IBD cross-section, specifically challenging the "factorization" of leptonic and nuclear matrix elements derived from the detailed-balance (db) principle.

1. Critique of Factorization: The standard approach assumes that leptonic wave functions ($\psi_{e,\nu}$) are spatially constant over the nuclear volume, allowing the transition amplitude to be factorized into a product of a leptonic factor and a nuclear matrix element ($M_{nuc}$). This permits the IBD nuclear matrix element to be inferred directly from the electron capture (EC) rate of the inverse process ($^{71}\text{Ge} \to ^{71}\text{Ga} + \nu_e$) via the detailed-balance relation. The authors argue this approximation is invalid when high precision is required, as the radial dependence of the leptonic wave functions cannot be separated from the nuclear integral.
2. New Formalism: The paper abandons the factorization scheme. Instead, it calculates the full transition amplitude by integrating the exact Dirac-Hartree-Fock-Slater (DHFS) radial solutions for the electron wave functions directly with the nuclear weak transition density, $\rho_{TD}(r)$. The transition density is defined as $\rho_{TD}(r) = \Psi^*_{71\text{Ge}}(r) \hat{H}_{GT} \Psi_{71\text{Ga}}(r)$, where $\hat{H}_{GT}$ is the Gamow-Teller Hamiltonian.
3. Phenomenological Parametrization: Since a first-principles calculation of $\rho_{TD}(r)$ with controlled uncertainties is currently challenging, the authors employ data-driven phenomenological parametrizations. They test several functional forms for $\rho_{TD}(r)$, including Single Gaussian (SG), Double Gaussian (DG), and modified versions (mDG, mTG) that incorporate polynomial powers of $r$ to mimic the nodal structures of shell-model wave functions.
4. Constraints: The models are constrained by two simultaneous conditions:
   • Reproducing the precisely measured half-life of $^{71}\text{Ge}$ ($t_{1/2} = 11.465 \pm 0.003$ d).
   • Minimizing the $\chi^2$ discrepancy between the theoretical IBD cross-section and the experimental values from GALLEX, SAGE, and BEST.

Key Results
The study demonstrates that the gallium anomaly can be resolved without invoking new physics, provided the transition density possesses specific radial structures.

• Resolution of the Anomaly: The authors find that transition densities with at least one node (sign change) can reduce the theoretical ground-state IBD cross-section by approximately 20%, bringing it into agreement with experimental measurements.
• Parametrization Performance:
  • The Single Gaussian (SG) model, lacking nodes, fails to resolve the tension.
  • The Double Gaussian (DG) model resolves the anomaly but requires an unphysically extended transition density.
  • The modified Double (mDG) and modified Triple Gaussian (mTG) models successfully resolve the anomaly while maintaining compact transition densities localized around the nuclear surface, consistent with standard nuclear structure expectations.
• Comparison with Standard Assumptions: The authors compare their results against a "two-parameter Fermi" (2pF) shape, which equates the weak transition density with the nuclear charge distribution (a common simplifying assumption). They show that this assumption yields a cross-section significantly larger than the experimental value, confirming that the charge distribution is an unjustified proxy for the Gamow-Teller transition density.
• Consistency: All successful models (DG, mDG, mTG) strictly satisfy the experimental constraints on the $^{71}\text{Ge}$ half-life, proving that the reduction in the IBD cross-section does not violate the known EC rate.

Significance and Claims
The paper claims to provide a "proof of principle" that the gallium anomaly is likely a consequence of theoretical approximations rather than new physics. Specifically:

1. Theoretical Correction: The work establishes that the standard factorization of leptonic and nuclear currents introduces a significant bias in the predicted cross-section. The detailed-balance principle, in its standard factorized form, is insufficient for high-precision IBD calculations.
2. Mechanism: The resolution relies on the interplay between exact lepton wave functions and a transition density with specific radial structure (nodes). The authors emphasize that this mechanism does not require anomalously long-ranged densities; compact, surface-localized densities with sign changes are sufficient.
3. Implications for New Physics: By offering a Standard Model-consistent explanation, the paper suggests that the need for sterile neutrino interpretations of the gallium anomaly is eliminated, aligning with recent null results from MicroBooNE and KATRIN.
4. Future Directions: The authors conclude that while their phenomenological models demonstrate the possibility of a solution, a definitive resolution requires the nuclear physics community to perform reliable, accurate microscopic calculations of the Gamow-Teller transition densities for the $^{71}\text{Ga} \leftrightarrow ^{71}\text{Ge}$ system. They note that if their scenario is correct, existing normalizations of excited-state contributions based on electron capture data may need to be revisited.