⚛️ nuclear theory

Efficient calculation of two-neutrino double-beta-decay nuclear matrix elements

This paper presents and validates an improved strength-function method based on Lanczos iterations that enables efficient and accurate calculation of two-neutrino double-beta-decay nuclear matrix elements without requiring full diagonalization, thereby facilitating the assessment of neutrinoless double-beta-decay rates and higher-order decay treatments.

Original authors: Mihai Horoi

Published 2026-06-29

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Mihai Horoi

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

The Big Picture: A Cosmic Puzzle

Imagine physicists are trying to solve a massive cosmic puzzle involving double-beta decay. This is a rare event where an atom's nucleus changes by turning two neutrons into two protons, releasing electrons in the process.

There are two versions of this event:

The "Allowed" Version (2νββ): This happens naturally and has been seen in labs. It's like a car driving down a known road.
The "Forbidden" Version (0νββ): This has never been seen. If found, it would prove that neutrinos are their own antiparticles and break a fundamental law of physics. This is like finding a car driving through a solid wall.

To find the "forbidden" version, scientists need to know exactly how the "allowed" version works. The paper focuses on calculating the Nuclear Matrix Element (NME) for the allowed version. Think of the NME as the "engine efficiency" of the atom. If we don't know the engine efficiency perfectly, we can't predict how fast the car should go, and we can't tell if a car driving through a wall is actually happening or just a glitch.

The Problem: Counting Too Many Stars

To calculate this "engine efficiency," scientists usually have to add up the contributions of every possible intermediate state the atom can pass through.

The Old Way: Imagine trying to count every single star in the night sky one by one to figure out how bright the sky is. For small neighborhoods (small atoms), you can do this. But for big cities (heavy atoms like Xenon or Tellurium), there are billions of stars. Counting them one by one takes forever, and your computer might crash before you finish.
The Convergence Issue: Sometimes, as you count, the numbers go up, then down, then up again. You might think you're done because the number looks stable, but then one more star changes everything. It's hard to know when you've truly reached the final answer.

The Solution: The "Shadow" Method

The author, M. Horoi, proposes a new, faster way to do this math. Instead of counting every single star, he uses a Lanczos strength-function method.

The Analogy:
Imagine you want to know the total weight of a giant, messy pile of sand.

The Old Way: You pick up every single grain of sand, weigh it, and add it to a list. This takes forever.
The New Way: You take a small, special scoop (the "doorway state") and dig into the pile. You don't need to see every grain. You just need to measure how the sand flows through your scoop and use a mathematical shortcut (the Lanczos matrix) to predict the total weight.

This method creates a "shadow" or a "map" of the whole pile based on a few key measurements. It avoids the need to list every single grain (state) individually.

What the Paper Actually Did

The author tested this new "scoop" method against the old "counting" method to see if it was accurate.

The Test Drive: He applied the method to several specific atoms (like Calcium-48, Germanium-76, and Xenon-136) that are popular targets for experiments.
The Comparison: For the smaller atoms (like Calcium-48), he could still count the grains one by one (calculate the full list of states). He compared his "scoop" result with the "counting" result.
- The Result: The "scoop" method matched the "counting" method perfectly. It reached the same final answer but did it much faster and with less memory.
The Heavy Lifting: For the bigger atoms (like Xenon-136), where counting every grain is impossible, the method still worked smoothly. It showed that the "shadow" method captures the most important parts of the pile without needing to see everything.
Bonus Feature: The author also showed that this same "scoop" can be used to calculate more complex, higher-order details (like the shape of the sand pile, not just the weight) which are needed for even more precise physics calculations.

Why This Matters

The paper claims that because the "allowed" decay (2νββ) and the "forbidden" decay (0νββ) are strongly linked, getting a fast and accurate calculation for the "allowed" one helps scientists guess the right settings for the "forbidden" one.

By using this efficient method, scientists can:

Test different theories (Hamiltonians) much faster.
Figure out the correct "quenching factor" (a correction number needed to make the math match reality).
Eventually, use these better numbers to set stricter limits on whether the "forbidden" decay exists.

In short: The paper introduces a mathematical shortcut that lets physicists calculate the behavior of complex atoms quickly and accurately, without needing supercomputers to count every single possibility. This helps them get closer to solving the mystery of the "forbidden" double-beta decay.

Technical Summary: Efficient Calculation of Two-Neutrino Double-Beta-Decay Nuclear Matrix Elements

Problem Statement
Reliable nuclear matrix elements (NMEs) are critical for interpreting double-beta decay ( $\beta\beta$ ) experiments and linking measured half-lives to fundamental weak-interaction physics. While the two-neutrino mode ( $2\nu\beta\beta$ ) is a Standard Model process observed in several nuclei, the neutrinoless mode ( $0\nu\beta\beta$ ) remains the primary signature for lepton-number violation. Recent statistical shell-model analyses indicate a strong correlation between $2\nu\beta\beta$ and $0\nu\beta\beta$ NMEs, suggesting that accurate calculations of the former can constrain the latter.

However, calculating $2\nu\beta\beta$ NMEs ( $M_{2\nu}$ ) presents a significant computational challenge. The standard approach requires an explicit summation over all intermediate $1^+$ states in the odd-odd nucleus. In large model spaces (e.g., heavy nuclei or expanded valence spaces like $sd-pf$), the number of contributing states becomes prohibitive, and the running sum may converge slowly or non-monotonically. Full diagonalization of the many-body Hamiltonian to obtain these states is often computationally infeasible, particularly for next-generation experiments aiming for sensitivities up to $10^{28}$ years.

Methodology
The paper presents an improved implementation of the strength-function method based on Lanczos iterations to calculate $M_{2\nu}$ without requiring full diagonalization. The approach follows the formalism of Ref. [24] and utilizes a numerical prescription detailed in Ref. [25].

Lanczos Iterations: The method initiates Lanczos iterations using Gamow-Teller (GT) "doorway" states constructed by applying the GT operator to the initial ( $0^+_i$ ) or final ( $0^+_f$ ) ground states.
Tridiagonal Matrix: Instead of storing full eigenvectors of the many-body Hamiltonian, the algorithm constructs a much smaller tridiagonal Lanczos matrix. The NME is approximated by summing over the eigenvalues and eigenvectors of this small matrix.
Coefficients: The expansion coefficients ( $g^\pm_m$ ) are derived from the Lanczos vectors and the energy denominators associated with the intermediate states.
Efficiency: This approach avoids the storage and processing costs of diagonalizing the full Hamiltonian. It tracks convergence iteratively using the Lanczos vectors themselves, offering a computational efficiency increase of approximately a factor of two compared to earlier strength-function implementations.
Extension to Higher Orders: The framework is extended to calculate higher-order NMEs ( $M_{2\nu-3}$ and $M_{2\nu-5}$ ) required for Taylor-expanded treatments of phase-space factors and related two-neutrino electron-capture observables. These terms involve higher powers of the energy denominator, making them more sensitive to low-lying states.

Key Contributions and Results
The author applies this method to several experimentally important isotopes ( $^{48}\text{Ca}$ , $^{76}\text{Ge}$ , $^{82}\text{Se}$ , $^{128,130}\text{Te}$ , and $^{136}\text{Xe}$ ) using standard effective Hamiltonians (e.g., GXPF1A, KB3G, JUN45, GCN2850, GCN5082, SVD) in corresponding valence spaces.

Benchmarking: The method is rigorously tested by comparing the "running sums" obtained from the Lanczos expansion against explicit sums over hundreds or thousands of intermediate $1^+$ $1^{+}$ states.
- For $^{48}\text{Ca}$ (pf-shell), the 25-iteration Lanczos result reproduces the direct sum over 1000 states, capturing the rapid rise, changes in strength regions, and the final stable plateau.
- For heavier nuclei like $^{82}\text{Se}$ and $^{136}\text{Xe}$ , where direct diagonalization of thousands of states is difficult, the Lanczos method reproduces the convergence behavior of direct sums (e.g., over 2000 states for $^{136}\text{Xe}$ ).
- The method successfully captures contributions from the GT giant resonance region and handles the non-trivial cancellations typical of $2\nu\beta\beta$ matrix elements.
Higher-Order Terms: Benchmarks for $M_{2\nu-3}$ and $M_{2\nu-5}$ in $^{48}\text{Ca}$ show excellent agreement with explicit sums. This confirms that the same Lanczos coefficients used for the leading NME can be generalized to higher inverse powers of the energy denominator.
Quenching Factors: The study determines optimal GT quenching factors ( $q_{opt}$ ) for various Hamiltonians to reproduce empirical matrix elements. The results highlight that these factors are phenomenological calibrations dependent on the specific interaction and valence space, rather than universal constants.

Significance and Claims
The paper claims that the proposed Lanczos strength-function method provides a practical and efficient solution for evaluating $M_{2\nu}$ in model spaces where direct summation is computationally expensive or impossible.

Accuracy: The method preserves the accuracy of explicit summation where benchmarks are possible, even in cases involving large cancellations and contributions from high-energy states.
Utility for $0\nu\beta\beta$ : Given the strong statistical correlation between $2\nu\beta\beta$ and $0\nu\beta\beta$ NMEs, reliable and efficient $2\nu\beta\beta$ calculations can help identify effective Hamiltonians and operator renormalizations that are favorable for $0\nu\beta\beta$ predictions.
Scalability: The efficiency of the method enables the exploration of enlarged model spaces (e.g., including $2\hbar\omega$ excitations) that were previously inaccessible via direct diagonalization.
Future Outlook: The author suggests that this framework can serve as a complementary large-scale benchmark for emerging ab initio calculations, providing statistical samples and phenomenological constraints derived from measured half-lives.

The work does not propose new experimental signatures or claim to solve the $0\nu\beta\beta$ problem directly but positions the efficient calculation of $2\nu\beta\beta$ NMEs as a necessary tool for reducing theoretical uncertainties in the broader search for neutrinoless decay.