⚛️ nuclear theory

Exact Calculation of Two-neutrino Double Beta Decay Rate

This paper proposes a rigorous theoretical framework that incorporates the full interdependence between nuclear structure and lepton kinematics to calculate two-neutrino double-beta decay rates, revealing deviations from traditional approximation methods for isotopes like $^{82}$ Se and $^{136}$ Xe.

Original authors: Stefan-Alexandru Ghinescu, Andrei Neacsu, Sabin Stoica

Published 2026-07-01

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Stefan-Alexandru Ghinescu, Andrei Neacsu, Sabin Stoica

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a nuclear physics experiment as a high-stakes dance. In this dance, a heavy atom (the "parent") wants to transform into a lighter one (the "child"). To do this, it has to throw off two electrons and two invisible particles called neutrinos. This specific dance move is called two-neutrino double-beta decay.

For decades, scientists trying to predict how fast this dance happens have had to use a "shortcut." They treated the two main parts of the dance as if they were happening in separate rooms:

The Nuclear Room: Where the heavy atoms shuffle their internal parts.
The Lepton Room: Where the flying electrons and neutrinos move around.

The Old Way (The Shortcut)
Traditionally, researchers calculated the "Nuclear Room" moves and the "Lepton Room" moves separately and then just multiplied the results together. They used approximations like "averaging" the energy of the middle steps or using a "Taylor expansion" (a mathematical way of guessing the curve of a line by looking at a few points).

Think of it like trying to predict the path of a ball thrown by a spinning dancer. The old method said, "Let's calculate how the dancer spins, then calculate how the ball flies, and assume they don't really affect each other." It was fast, but it introduced errors—sometimes up to 25% off.

The New Way (The Full Picture)
In this paper, the authors (Ghinescu, Neacsu, and Stoica) say, "Wait a minute. The dancer and the ball are connected. The way the dancer spins changes how the ball flies, and the ball's flight changes how the dancer spins."

They propose a new, "exact" calculation method. Instead of separating the rooms, they keep the whole dance floor open. They calculate the nuclear structure and the electron movement simultaneously, accounting for every tiny interaction between the heavy nucleus and the flying electrons.

What They Found
They tested this new method on two specific atoms: Selenium-82 and Xenon-136.

The Speed: Their new, more realistic calculation showed that these atoms decay 8% to 14% faster than the old methods predicted.
- Analogy: If the old map said a trip would take 100 minutes, the new map says it actually takes 90 minutes. The atoms are "aging" faster than we thought.
The Energy: They also looked at the "music" of the dance—the energy spectrum of the electrons.
- For Selenium-82, the difference between the old and new methods was huge (up to 20% in some parts of the energy range). It's like the old map said the music was a slow waltz, but the new map reveals it's actually a fast-paced salsa in certain moments.
- For Xenon-136, the difference was much smaller. The old shortcuts worked okay here, but the new method is still more precise.

Why It Matters
The authors aren't claiming this fixes a disease or builds a new engine. They are simply saying: "Our math is now more honest."

Because the new method shows the atoms decay faster, it suggests that the "adjustment knobs" (called quenching factors) scientists use to make their nuclear models match reality might need to be turned. It also gives experimentalists a better "ruler" to measure the dance. If they can measure the energy of the electrons very precisely, they can now tell if their theories about how the nucleus is built are correct, because the new math provides a much sharper target to hit.

In a Nutshell
The authors stopped treating the nucleus and the electrons as strangers in separate rooms. By letting them interact in the calculation, they found that the decay happens faster and the energy patterns are slightly different than previously thought, particularly for Selenium-82. This provides a more accurate foundation for future experiments trying to understand the fundamental rules of the universe.

Technical Summary: Exact Calculation of Two-neutrino Double Beta Decay Rate

Problem Statement
The calculation of two-neutrino double-beta decay ( $2\nu\beta\beta$ ) rates has historically relied on approximations that decouple nuclear structure from lepton kinematics. Traditional methods, including the closure, non-closure, and Taylor expansion approximations, factorize the decay rate into independent Nuclear Matrix Elements (NMEs) and Phase-Space Factors (PSFs).

Closure Approximation: Replaces intermediate nuclear state energies ( $E_N$ ) with an average energy $\langle E_N \rangle$ and approximates lepton energies, introducing errors up to 25% in NME calculations.
Taylor Expansion: Retains lepton energies in the denominator up to the fourth order but remains computationally intensive and mathematically questionable for certain isotopes where the energy denominator approaches zero (e.g., when the first $1^+$ state is close to the ground state of the intermediate nucleus).
These approximations fail to capture the full interdependence between nuclear structure and atomic/lepton effects, leading to uncertainties in decay rate predictions and electron spectra.

Methodology
The authors propose a novel, exact approach that integrates NMEs and PSFs into a unified formulation, preserving the full coupling between nuclear structure and lepton kinematics.

Formalism: Instead of factorizing the decay rate, the method computes the integral for each intermediate $1^+$ state individually and sums the results. The decay rate $\Gamma$ is calculated by integrating over electron energies ( $\epsilon_1, \epsilon_2$ ) and neutrino energies ( $\omega_1, \omega_2$ ) without approximating the energy denominators in the Gamow-Teller matrix elements.
Implementation: The authors applied this method to the isotopes $^{82}\text{Se}$ $^{82} Se$ and $^{136}\text{Xe}$ $^{136} Xe$ .
- Nuclear structure inputs ( $1^+$ spectra and Gamow-Teller strengths) were generated using the KSHELL program with the jun45 interaction for $^{82}\text{Se}$ (summing 50 levels) and the svd interaction for $^{136}\text{Xe}$ (summing 250 levels).
- Atomic effects, including electron wave functions and screening, were computed using the RADIAL package with a modified self-consistent Dirac-Hartree-Fock-Slater treatment.
- Integrals were evaluated using the SPADES package.
Scope: The study focuses on the $2\nu\beta\beta$ decay mode, calculating both total decay rates and differential distributions (single-electron and summed electron energy spectra). Extensions to other isotopes and the neutrinoless mode ( $0\nu\beta\beta$ ) are noted as future work.

Key Contributions

Unified Framework: The paper introduces a computational framework that eliminates the artificial decoupling of nuclear and atomic degrees of freedom, offering a more rigorous treatment of the $2\nu\beta\beta$ process.
Exact Calculation: By avoiding closure and Taylor expansion approximations, the method accounts for the specific energy dependence of intermediate states and lepton kinematics simultaneously.
Differential Spectra: The authors provide exact calculations for single-electron and summed electron energy spectra, allowing for a direct comparison of spectral shapes between exact and approximate methods.

Results

Decay Rates: For both $^{82}\text{Se}$ and $^{136}\text{Xe}$ , the exact method yields decay rates that are 8–14% higher than those obtained via traditional approximate methods (Closure and Taylor expansion). This corresponds to predicted half-lives that are shorter than those derived from approximate models.
Electron Spectra:
- $^{82}\text{Se}$ : The single-electron spectra show differences of up to 20% between the exact method and approximate methods. These deviations fall within the experimentally accessible range.
- $^{136}\text{Xe}$ : The differences are significantly smaller. This is attributed to the fact that the energy term $E_N - \frac{1}{2}(E_I + E_F)$ lies approximately 680 keV above $Q_{\beta\beta}$ , making the non-closure and Taylor approximations inherently more accurate for this isotope.
Half-Life Comparison: Theoretical half-lives calculated via the exact method are consistently lower than experimental values and those predicted by closure/Taylor methods, suggesting that current quenching factors used in NME calculations may require adjustment when using this exact formalism.

Significance and Claims
The authors claim that their approach provides a more realistic description of the $2\nu\beta\beta$ process by explicitly accounting for the interdependence between nuclear and atomic effects.

Validation Tool: The sensitivity of the single-electron spectra to the calculation method, particularly for $^{82}\text{Se}$ , offers a potential experimental tool for validating nuclear structure models and testing hypotheses regarding Single-State Dominance (SSD) versus Higher-State Dominance (HSD).
Theoretical Implications: The observed discrepancies suggest that the quenching factors introduced in NME calculations to reproduce experimental Gamow-Teller strengths may need re-evaluation when applied within this exact framework.
Future Directions: The paper positions this framework as a foundation for future investigations into nuclear and atomic effects, with extensions to other isotopes and the $0\nu\beta\beta$ mode currently underway. The authors emphasize that while their method increases computational complexity, it remains manageable with modern computing resources.

Technical Summary: Exact Calculation of Two-neutrino Double Beta Decay Rate

More like this