⚛️ nuclear theory

The renormalization of the shell-model neutrinoless double-beta decay operator starting from effective field theory (I)

This paper presents the first fully consistent shell-model calculation of neutrinoless double-beta decay matrix elements for 48Ca, 76Ge, and 82Se, deriving both the nuclear Hamiltonian and decay operators from chiral perturbation theory via many-body perturbation theory while also assessing convergence and theoretical uncertainties.

Original authors: L. Coraggio, G. De Gregorio, S. L. Lyu, N. Itaco

Published 2026-06-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: L. Coraggio, G. De Gregorio, S. L. Lyu, N. Itaco

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: Predicting a "Ghostly" Event

Imagine you are trying to predict the outcome of a very rare, almost magical event: a specific type of atomic nucleus (like a tiny solar system) spontaneously changing into a different one without emitting any neutrinos (ghostly particles). This is called neutrinoless double-beta decay.

If this happens, it proves that neutrinos are their own antiparticles and helps us understand why the universe has more matter than antimatter. However, this event is so rare that we can't just watch it happen in a lab yet. We have to calculate how likely it is to happen using complex math.

This paper is about building a better, more consistent calculator to predict this likelihood.

The Problem: The "Mismatched Toolkit"

For a long time, scientists trying to calculate this used two different toolkits that didn't quite fit together:

The Engine (The Nucleus): They used a set of rules to describe how protons and neutrons stick together.
The Trigger (The Decay): They used a different set of rules to describe how the decay happens.

It's like trying to drive a car where the engine was built by one mechanic using a specific blueprint, but the steering wheel was built by another mechanic using a completely different blueprint. The car might run, but it won't handle well, and you can't be sure if your calculations are accurate.

The Solution: The "Chiral" Blueprint

The authors of this paper decided to build both the engine and the steering wheel using the same master blueprint. This blueprint is called Chiral Perturbation Theory (ChPT).

Think of ChPT as a universal language for the strong force (the glue holding atoms together). It's based on the fundamental rules of the universe (Quantum Chromodynamics). By using this single language for both the nucleus and the decay process, the authors ensure that their "calculator" is perfectly consistent.

The Method: The "Shell Model" and "Renormalization"

To do the math, they used a method called the Shell Model.

The Analogy: Imagine the nucleus as a hotel with many floors (energy levels). Protons and neutrons are guests trying to find rooms. The "Shell Model" is a way to figure out exactly which guests are in which rooms and how they interact.
The Challenge: The hotel is huge, and there are too many rooms to check every single one. So, scientists pick a few "key floors" (a model space) to focus on.
The Fix (Renormalization): When you ignore the other floors, you lose some information. To fix this, the authors use a process called renormalization.
- Analogy: Imagine you are summarizing a 1,000-page book into a 10-page summary. You can't just delete the other 990 pages; you have to rewrite the 10 pages you keep so they feel like the whole book. The "renormalization" is the act of rewriting the rules for the 10 pages so they account for the influence of the 990 pages you left out.

What They Did

The team applied this consistent, "renormalized" method to three specific atomic nuclei: Calcium-48, Germanium-76, and Selenium-82. These are the top candidates scientists are currently looking at to find evidence of this decay.

They did two main things:

Tested the Engine: They first checked if their new, consistent rules could accurately predict the normal behavior of these atoms (like their energy levels and how they vibrate). They found that their "third-order" calculations (a high level of detail) matched real-world experiments very well.
Predicted the Decay: They used these validated rules to calculate the probability (the "Nuclear Matrix Element") of the neutrinoless double-beta decay happening.

The Results: A "Contact Term" Surprise

One of the most interesting findings is about a specific part of the math called the short-range contact term.

The Old View: Scientists used to think the decay happened only through a long-distance interaction (like two people shouting across a field).
The New View: The authors' consistent blueprint shows that there is also a "contact" interaction (like two people bumping into each other in a crowded room). This short-range bump is crucial for the math to work correctly. They had to include this "bump" to make the calculation stable.

Uncertainty: How Sure Are We?

In science, it's not enough to get a number; you need to know how much you can trust it.

The authors used a mathematical trick called a Padé approximant (think of it as a "best guess" curve that smooths out the bumps in their calculations) to estimate the error.
They found that while their numbers are solid, there is still some uncertainty because they are using a "perturbative" approach (building the answer step-by-step).
The Takeaway: Their calculated values for how likely the decay is are generally lower than some other recent studies, but they are consistent within their own error margins.

Summary

This paper is a major step forward because it stops using "mixed" toolkits. Instead, it builds the entire calculation for these rare atomic decays from the ground up using a single, consistent set of physical laws. They proved their method works by matching it to real-world data, and they provided a new, carefully calculated estimate of how likely these "ghostly" decays are to occur in Calcium, Germanium, and Selenium.

What they did NOT do:

They did not claim to have discovered the decay.
They did not claim to have solved the mystery of neutrino mass directly (they just provided a better tool to help solve it).
They did not apply this to medical or clinical uses.
They did not calculate results for other nuclei like Molybdenum-100 (though they mentioned they plan to do that in the future).

Technical Summary: Renormalization of the Shell-Model Neutrinoless Double- $\beta$ Decay Operator from Effective Field Theory

Problem Statement
The calculation of nuclear matrix elements ( $M^{0\nu}$ ) for neutrinoless double- $\beta$ ( $0\nu\beta\beta$ ) decay is critical for interpreting experimental limits on the effective Majorana neutrino mass. A persistent challenge in this field is the "quenching" of the axial coupling constant ( $g_A$ ) and the lack of a fully consistent framework where both the nuclear Hamiltonian and the decay operators are derived from the same underlying theory. Standard shell-model calculations often employ phenomenological effective Hamiltonians paired with standard decay operators, potentially missing short-range correlations and renormalization effects required by Quantum Chromodynamics (QCD). Specifically, the $0\nu\beta\beta$ operator derived within Chiral Perturbation Theory (ChPT) includes a short-range contact term at leading order (LO) that is absent in standard formulations, necessitating a consistent renormalization procedure.

Methodology
This work presents a fully consistent shell-model (SM) calculation for the $0\nu\beta\beta$ decay of $^{48}\text{Ca}$ , $^{76}\text{Ge}$ , and $^{82}\text{Se}$ . The approach integrates ChPT with many-body perturbation theory (MBPT) through the following steps:

Hamiltonian Construction: The effective SM Hamiltonian ( $H_{\text{eff}}$ ) is derived from a high-precision two-nucleon (2N) potential constructed by Entem and Machleidt at next-to-next-to-next-to-leading order (N3LO) in ChPT. This is supplemented by a three-nucleon (3N) force derived at N2LO. The low-energy constants (LECs) $c_D$ and $c_E$ for the 3N force are fixed to reproduce observables in $A=3$ systems ( $c_D = -1, c_E = -0.34$ ).
Operator Derivation: The $0\nu\beta\beta$ decay operator is constructed from ChPT currents at LO. This includes the standard long-range components (Fermi, Gamow-Teller, and Tensor) and a short-range contact term regulated by a Gaussian cutoff ( $\Lambda = 500$ MeV). The LEC for the contact term, $g_{NN}^\nu$ , is estimated as $-0.47 \text{ fm}^2$ based on charge-independence breaking considerations.
Renormalization via MBPT: Both $H_{\text{eff}}$ and the effective decay operators ( $\Theta_{\text{eff}}$ ) are derived using the Kuo-Lee-Ratcliff folded-diagram expansion. The $\hat{Q}$ -box and $\hat{\Theta}$ -box vertex functions are expanded up to the third order in perturbation theory. The model spaces are defined as the $0f1p$ shell (for $^{48}\text{Ca}$ ) and the $0f5/21p0g9/2$ shell (for $^{76}\text{Ge}$ and $^{82}\text{Se}$ ), built upon $^{40}\text{Ca}$ and $^{56}\text{Ni}$ cores, respectively.
Convergence and Uncertainty: The study investigates the perturbative convergence of spectroscopic properties, two-neutrino ( $2\nu\beta\beta$ ) matrix elements, and $0\nu\beta\beta$ matrix elements. Theoretical uncertainties arising from the truncation of the MBPT expansion are estimated using Padé approximants $[2|1]$ .

Key Contributions

First Consistent Framework: This is the first application of a shell-model calculation for $0\nu\beta\beta$ matrix elements where both the nuclear Hamiltonian and the decay operators are derived entirely from ChPT and renormalized consistently via MBPT.
Inclusion of Short-Range Terms: The work explicitly incorporates the LO short-range contact operator required for the renormalizability of the $0\nu\beta\beta$ process within the light Majorana neutrino exchange scenario.
Perturbative Analysis: The authors provide a detailed analysis of the order-by-order convergence of the matrix elements, demonstrating that while individual components (Gamow-Teller, Fermi, Tensor) converge well, the total $M^{0\nu}$ is sensitive to cancellations between components, particularly in $^{48}\text{Ca}$ .

Results

Spectroscopy: The calculated low-energy spectra and $B(E2)$ transition strengths for $^{48}\text{Ca}$ , $^{48}\text{Ti}$ , $^{76}\text{Ge}$ , $^{76}\text{Se}$ , $^{82}\text{Se}$ , and $^{82}\text{Kr}$ show satisfactory agreement with experimental data, particularly when using third-order $H_{\text{eff}}$ . The results validate the reliability of the nuclear wave functions generated by this approach.
$2\nu\beta\beta$ Decay: The calculated $M^{2\nu}$ matrix elements for ground-state to ground-state transitions agree well with experimental values when using third-order Hamiltonians, further supporting the validity of the wave functions.
$0\nu\beta\beta$ Matrix Elements: The calculated $M^{0\nu}$ $M^{0 ν}$ values (in units of dimensionless matrix elements) are:
- $^{48}\text{Ca} \to ^{48}\text{Ti}$ : $0.9$ (3rd order) vs. $0.8$ (Padé).
- $^{76}\text{Ge} \to ^{76}\text{Se}$ : $1.7$ (3rd order) vs. $1.5$ (Padé).
- $^{82}\text{Se} \to ^{82}\text{Kr}$ : $1.3$ (3rd order) vs. $1.1$ (Padé).
Uncertainty Estimation: The difference between the third-order result and the Padé approximant is used to estimate the theoretical uncertainty associated with the perturbative expansion. The study finds that the uncertainty is driven largely by the cancellation between the Gamow-Teller and other components.
Comparison: The results are compared with recent ab initio and shell-model studies. While the absolute values differ from some phenomenological calculations (e.g., those using empirical Hamiltonians), the authors note that their first-order results align more closely with the central values of studies using N2LO currents and empirical Hamiltonians, suggesting that the consistent renormalization significantly impacts the final matrix elements.

Significance and Claims
The authors claim that this work establishes a robust, theoretically consistent foundation for calculating $0\nu\beta\beta$ decay rates. By deriving both the Hamiltonian and the operators from ChPT, the study addresses the issue of operator renormalization and the quenching of $g_A$ within a unified framework. The authors emphasize that their results provide a baseline for estimating theoretical uncertainties specifically related to the many-body perturbation expansion. They modestly state that while their approach is consistent, a complete error budget would also require assessing uncertainties from the ChPT expansion of the Hamiltonian and currents themselves, which is a subject of ongoing work. The paper concludes by outlining future plans to include higher-order (N2LO) contributions to the $0\nu\beta\beta$ currents and to apply this framework to the more complex $^{100}\text{Mo}$ system.