Physicists Ditch Decades-Old Approximations in Double-Beta Decay Breakthrough

In a development that could reshape how nuclear physicists interpret some of the rarest events in the universe, a team of researchers has unveiled a new method for calculating the rate of two-neutrino double-beta decay that abandons the approximations physicists have relied on for decades. The study, published in The European Physical Journal C by S. A. Ghinescu, A. Neacsu, and S. Stoica of the International Center for Advanced Training and Research in Physics in Măgurele, Romania, together with colleagues at the Horia Hulubei National Institute of Physics and Nuclear Engineering, demonstrates that the traditional way of splitting the decay calculation into two independent pieces—nuclear structure on one side, lepton kinematics on the other—leaves measurable fingerprints in the predicted decay rates and, crucially, in the energy spectra of the emitted electrons.

Double-beta decay is a process in which an atomic nucleus converts two neutrons into two protons at once, emitting two electrons and two electron antineutrinos. It is extraordinarily rare, with half-lives that can exceed 10^19 years, yet it has become one of the most actively pursued topics in modern physics. The reason is that its hypothetical sibling, neutrinoless double-beta decay, in which no antineutrinos appear at all, would violate the rules of the Standard Model and prove that neutrinos are their own antiparticles—a discovery with profound implications for why the universe contains matter at all. The ordinary two-neutrino mode, which has been observed in several isotopes, serves as the indispensable benchmark: experiments measure its half-lives and electron spectra with ever-greater precision, and theorists must match those measurements with equally precise calculations to validate the nuclear structure models that underpin the search for the neutrinoless mode.

The problem, as the Romanian team points out, lies in a computational shortcut inherited from an era when the full calculation was simply impossible. The decay rate is, in principle, a single quantity that entangles the quantum structure of the decaying nucleus with the energies and momenta of the four emitted leptons. To make the mathematics tractable, theorists historically factorized it into a nuclear matrix element—a number encoding how easily the nucleus transitions through its intermediate states—multiplied by a phase-space factor describing the lepton kinematics. This factorization rests on approximations known as closure, non-closure, and Taylor expansion methods. In the closure approximation, the energies of the intermediate nuclear states are replaced by a single average energy, and the lepton energies are approximated as half of the total energy released by the decay.

That shortcut carries real risks. Because the emitted leptons carry only a few mega-electronvolts of energy, the decay rate is acutely sensitive to the exact positions of the lowest-lying intermediate states, particularly the first excited 1+ state of the intermediate nucleus. Replacing the true spectrum of these states with an average energy can therefore introduce significant errors into the nuclear matrix element and, by extension, into the predicted half-life. A more recent refinement, the Taylor expansion method, retained the lepton energies in the denominator of the matrix element up to the fourth order, improving the description and enabling new experimental analyses of electron spectra and the axial-vector coupling constant. But even this method keeps the nuclear and lepton degrees of freedom decoupled, is computationally intensive, and—critically—its mathematical validity becomes questionable for certain isotopes. In nuclei such as molybdenum-100 and selenium-82, where the first 1+ state of the intermediate nucleus lies at or very near the ground state, the expansion parameter approaches unity and the series may simply fail to converge in parts of the integration domain.

The new approach, which the authors call the Direct method, is conceptually disarmingly simple: it computes the decay rate exactly as the underlying theory defines it, without any approximation. For each intermediate 1+ state of the daughter nucleus, the team evaluates one integral in which the nuclear transition amplitudes and the lepton energies appear together, exactly as they do in nature, and then sums the contributions over all the intermediate states. The computational cost is higher, but with modern computing resources and summations over roughly 30 to 50 excited 1+ states, the task is entirely manageable. The result is a unified formula that preserves the full interdependence between nuclear structure and lepton kinematics—something no previous calculation of the electron spectra had achieved.

To test the method, the researchers applied it to two of the most important double-beta decay isotopes: selenium-82 and xenon-136. Both are staples of major experimental programs, including KamLAND-Zen and PandaX-4T for xenon-136, and both have precisely measured two-neutrino half-lives and spectra. The nuclear structure input was generated with the Antoine shell-model code, using the jun45 effective Hamiltonian in the jj44 model space for selenium-82 and the svd Hamiltonian in the jj55 space for xenon-136, with the summation over intermediate states truncated at 30 levels, a cutoff previously shown to be sufficient for convergence. The electron wave functions were computed with the RADIAL package using a self-consistent Dirac–Hartree–Fock–Slater treatment of the daughter atom, accounting for screening effects and the final nuclear size, while the multidimensional integrals were evaluated with the SPADES package. Notably, the team applied no quenching of the nuclear matrix elements or of the axial-vector coupling constant, and fitted no parameters to experimental data—a deliberate choice designed to isolate the differences between calculation methods rather than to chase agreement with measurements.

The results are striking. For both isotopes, the Direct method yields decay rates up to 16 percent higher than the approximate methods, corresponding to predicted half-lives that are up to 16 percent shorter. While the authors caution that a definitive interpretation requires extending the analysis to more isotopes and alternative nuclear models, the discrepancy plausibly stems from the nuclear-structure–lepton-kinematics interdependence that only the new method captures. Even more experimentally consequential are the spectral predictions. The single-electron energy spectrum—the distribution of energies of one of the two emitted electrons—shows differences of up to 20 percent between the Direct and approximate calculations for selenium-82, a deviation squarely within the reach of current experimental precision. The summed electron spectrum, by contrast, varies by less than about 10 percent between methods, and the angular correlation between the two emitted electrons by at most 5 percent. For xenon-136, where the relevant nuclear energy denominator sits roughly 680 kilo-electronvolts above the decay Q-value, the approximations perform much better and the differences shrink accordingly—exactly as the structure of the underlying formulas would predict.

This hierarchy of sensitivities is what makes the work potentially viral in its implications for the experimental community. Precision measurements of the two-neutrino electron spectrum have already been used by experiments such as KamLAND-Zen, CUPID-MO, and CUORE to constrain the quenching of nuclear matrix elements and even to probe the effective value of the axial-vector coupling constant. If the single-electron spectrum is genuinely more sensitive to the exact treatment of the nuclear-lepton coupling than previously appreciated, then existing and upcoming spectral measurements carry information that the standard analysis frameworks have been systematically missing. The authors also highlight that their framework is naturally suited to testing the Single-State Dominance and Higher-State Dominance hypotheses—competing pictures of how much each intermediate 1+ state contributes to the total decay amplitude—which remain contested in the literature and are difficult to discriminate experimentally without precisely calculated spectral shapes.

The team is already extending the framework to other isotopes and, more ambitiously, to the neutrinoless decay mode itself, where the closure approximation has traditionally been considered safe because the virtual neutrino carries momenta of order 100 mega-electronvolts, diluting the sensitivity to the intermediate-state energies. Whether the full interdependence between nuclear structure and lepton kinematics leaves detectable traces in the neutrinoless channel—where the stakes include nothing less than the nature of the neutrino and the origin of matter in the universe—remains an open and tantalizing question. For now, the message of the Direct method is clear: the approximations that made double-beta decay calculations possible for half a century are no longer necessary, and retiring them reveals measurable differences exactly where the next generation of experiments is looking.

Subject of Research: Direct calculation of two-neutrino double-beta decay rates with full nuclear-structure and lepton-kinematics coupling

Article Title: Direct calculation of the two-neutrino double-beta decay rate including the full lepton-energy dependence

Article References: Ghinescu, S. A., Neacsu, A., & Stoica, S. (2026). Direct calculation of the two-neutrino double-beta decay rate including the full lepton-energy dependence. The European Physical Journal C, 86(9), Article 1123. https://doi.org/10.1140/epjc/s10052-026-16356-x

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16356-x

Keywords: double-beta decay, two-neutrino decay, nuclear matrix elements, phase-space factors, closure approximation, Taylor expansion, electron spectra, selenium-82, xenon-136, nuclear shell model, neutrinoless decay, axial-vector coupling