Neutrinoless Double Beta Decay (0νββ)

• Neutrinoless double beta decay (0νββ) is a hypothetical nuclear process where two neutrons convert into two protons without neutrino emission, signaling lepton number violation.
• The process is analyzed through various nuclear models and experimental techniques to determine effective Majorana mass and examine physics beyond the Standard Model.
• Experimental approaches focus on achieving precise energy resolution and ultra-low background conditions to detect half-lives exceeding 10^26 years and constrain new physics scenarios.

Neutrinoless double beta decay (0νββ) is a hypothesized second-order weak nuclear process in which two neutrons in a nucleus are simultaneously converted into two protons, emitting two electrons and no neutrinos: $(A, Z) \to (A, Z+2) + 2e^{-}$ This process is forbidden in the Standard Model by the conservation of total lepton number (ΔL = 2), but is generically allowed in theories where neutrinos are Majorana particles (i.e., they are their own antiparticles). Evidence for 0νββ would unambiguously demonstrate lepton-number violation, test the Majorana nature of neutrinos, and directly probe fundamental physics beyond the Standard Model. Today, experimental searches constrain the half-life for this decay to exceed $10^{25}$–$10^{26}$ years in favorable isotopes, setting crucial limits on the effective Majorana mass and the scale of lepton-number–violating new physics (Cardani, 2018); see also (Dolinski et al., 2019, Brugnera, 17 Jan 2025).

1. Theoretical Framework and Significance

0νββ is expected in extensions of the Standard Model where neutrinos acquire Majorana masses, as in seesaw scenarios. The “black-box theorem” (Schechter–Valle) establishes that any process inducing 0νββ also generates a Majorana mass term for the neutrino at some order (Cardani, 2018, Deppisch et al., 2012). The canonical mechanism (“mass mechanism”) involves the exchange of light Majorana neutrinos between two $\beta$-decay vertices, requiring a helicity flip proportional to the neutrino mass $m_i$:

$$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$

where $q \sim 100$ MeV is the typical virtual momentum, and the effective parameter is

$$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$

with $U_{ei}$ the elements of the PMNS mixing matrix and $m_i$ the neutrino mass eigenvalues (Cardani, 2018, Rodejohann, 2010, Bilenky et al., 2012).

Observation of 0νββ would have several profound implications:

• Demonstrate violation of lepton number by two units (ΔL = 2)
• Establish that neutrinos are Majorana fermions ($10^{25}$0)
• Provide access to the absolute neutrino mass scale and Majorana CP-violating phases
• Constrain or discriminate between normal and inverted neutrino mass hierarchies
• Impact scenarios for baryogenesis via leptogenesis, since lepton-number violation is a prerequisite for generating the matter–antimatter asymmetry of the Universe (Cardani, 2018, Dolinski et al., 2019)

2. Formalism and Rate Formulae

The inverse half-life for $10^{25}$1 via light-neutrino exchange is factorized as:

$10^{25}$2

where:

• $10^{25}$3 is the exactly calculable phase-space factor (yr$10^{25}$4), scaling as $10^{25}$5; typical values are $10^{25}$6–$10^{25}$7 yr$10^{25}$8 for $10^{25}$9–3 MeV
• $10^{26}$0 is the nuclear matrix element (NME), a dimensionless quantity encapsulating nuclear structure; values depend on the calculation method and typically lie in the range $10^{26}$1–7 (uncertainty factor 2–3)
• $10^{26}$2 is the effective Majorana mass, as above
• $10^{26}$3 MeV is the electron mass (Cardani, 2018, Faessler, 2011, Bilenky et al., 2012)

The structure of $10^{26}$4 is conventionally decomposed as:

$10^{26}$5

where $10^{26}$6 (Gamow–Teller), $10^{26}$7 (Fermi), and $10^{26}$8 (tensor) are evaluated using various many-body nuclear approaches (Grebe, 1 Apr 2025, Faessler, 2011, Rodejohann, 2010).

Oscillation data tightly constrain the mixing angles and mass-squared splittings, but the absolute mass scale and Majorana phases remain unconstrained. For the allowed ranges:

• Inverted mass ordering ($10^{26}$9): $\beta$0–50 meV
• Normal ordering ($\beta$1): $\beta$2 meV This corresponds to expected $\beta$3 half-lives $\beta$4–$\beta$5 yr (inverted hierarchy) or $\beta$6 yr (normal hierarchy), for typical $\beta$7 and $\beta$8 (Cardani, 2018, Dolinski et al., 2019).

3. Nuclear Matrix Elements and Theoretical Uncertainties

Computing $\beta$9 is the main theoretical challenge. The principal methods are:

• Nuclear Shell Model (SM): truncation to valence space, full correlations; typically gives lower $m_i$0 due to limited orbitals
• Quasiparticle Random Phase Approximation (QRPA): large single-particle space, includes more intermediate-state correlations; $m_i$1 parameter adjusted to reproduce $m_i$2
• Interacting Boson Model (IBM-2): maps pairs to bosons, good global trends
• Projected Hartree–Fock–Bogoliubov (PHFB) and Energy-Density Functional (EDF) approaches: include pairing, deformation, and multi-reference correlations

Representative matrix element ranges (not exhaustive) are: | Isotope | Shell Model | QRPA | IBM-2 | PHFB | |-----------------|------------|---------|--------|--------| | $m_i$3Ge | 2–3 | 2–6 | 2–6 | 3–5 | | $m_i$4Te | 2–5 | 2–5 | 2–5 | ... | | $m_i$5Xe | 1.5–4 | 1.6–3.5 | 2.5–4.5| ... |

Sources of NME uncertainty include model space size, treatment of short-range correlations, quenching of the axial coupling $m_i$6, nuclear deformation, and omitted two-body weak currents. Discrepancies among methods contribute a factor-of-2–3 systematic uncertainty in inferred $m_i$7 from a measured $m_i$8 (Cardani, 2018, Grebe, 1 Apr 2025, Faessler, 2011, Faessler, 2012).

Lattice QCD and nuclear effective field theory are emerging as complementary tools for systematically reducing such uncertainties, providing direct calculation of certain low-energy constants and contact terms relevant for $m_i$9 (Grebe, 1 Apr 2025, Detmold et al., 2020).

4. Non-Standard Mechanisms and Beyond–Standard-Model Physics

While the “mass mechanism” is canonical, other processes can induce $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$0:

• Heavy Majorana neutrino exchange: left–right symmetric models with heavy $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$1 states contribute short-range operators; amplitude scales as $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$2
• Right-handed currents: new $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$3 bosons or mixing can alter electron chirality and the angular spectrum
• Supersymmetric (SUSY) models: R-parity–violating couplings enable squark–gluino or slepton–neutralino mediated operators, often via short-range diagrams
• Leptoquarks, Higgs triplets, and other exotic mediators: each introduces higher-dimensional ($$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$4) operators with characteristic operator structures and nuclear responses (Deppisch et al., 2012, Dolinski et al., 2019, Banerjee et al., 9 Oct 2025)

The master half-life formula for multiple mechanisms is:

$$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$5

where $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$6 are particle-physics parameters (typically dimensionless couplings or mass ratios). The contributions may interfere constructively or destructively; CP-violating phases and the structure of nuclear operators determine observable signatures (Faessler, 2012, Banerjee et al., 9 Oct 2025).

Disentangling the dominant underlying mechanism can potentially be achieved by:

• Comparing $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$7 rates across multiple isotopes (different $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$8 and operator sensitivities)
• Analyzing event kinematics (electron angular and energy distributions) in tracking detectors
• Correlating with complementary high-energy searches (colliders, $$\mathcal{A}_{\text{light}} \propto G_F^2\,\frac{m_{\beta\beta}}{q^2}$$9, direct searches for HNLs or SUSY particles) (Deppisch et al., 2012, Faessler, 2012, Dolinski et al., 2019).

5. Experimental Searches and Constraints

Experimental searches exploit isotopes with favorable $q \sim 100$0-values ($q \sim 100$1 MeV) and long $q \sim 100$2 lifetimes, using technologies optimized for background rejection, energy resolution, and large masses. Leading techniques include:

• High-purity Ge diodes: GERDA, Majorana Demonstrator, LEGEND ($q \sim 100$3Ge)
• Liquid/gaseous Xe TPCs: EXO-200, KamLAND-Zen, nEXO, NEXT ($q \sim 100$4Xe)
• Bolometric calorimeters: CUORE, CUPID ($q \sim 100$5Te, $q \sim 100$6Mo, $q \sim 100$7Se)
• Large liquid scintillator detectors: SNO+ ($q \sim 100$8Te), KamLAND-Zen ($q \sim 100$9Xe)
• Tracking calorimeters: NEMO-3, SuperNEMO ($$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$0Se, $$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$1Mo)

Recent 90% C.L. half-life limits and corresponding $$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$2 constraints (using a range of $$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$3) are summarized below (Cardani, 2018, Dolinski et al., 2019, Brugnera, 17 Jan 2025):

Isotope	Experiment	$$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$4 [yr]	$$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$5 Bound [meV]
$$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$6Xe	KamLAND-Zen	$$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$7	$$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$8–$$m_{\beta\beta} = \left| \sum_{i=1}^3 U_{ei}^2 m_i \right|$$9
$U_{ei}$0Ge	GERDA Phase II	$U_{ei}$1	$U_{ei}$2–$U_{ei}$3
$U_{ei}$4Xe	EXO-200	$U_{ei}$5	$U_{ei}$6–$U_{ei}$7
$U_{ei}$8Te	CUORE	$U_{ei}$9	$m_i$0–$m_i$1
$m_i$2Mo	NEMO-3	$m_i$3	$m_i$4–$m_i$5

The leading current and next-generation experiments aim to cover the entire inverted ordering region ($m_i$6–50 meV), targeting sensitivities $m_i$7–$m_i$8 yr (Brugnera, 17 Jan 2025, Guinn et al., 2019, Agostini et al., 2018, Grebe, 1 Apr 2025).

Key experimental challenges include:

• Background suppression: exploiting deep underground laboratories, ultra-pure materials, active veto systems (liquid argon, scintillator), event topology (tracking, pulse-shape discrimination)
• Excellent energy resolution: crucial for distinguishing the monoenergetic 0νββ peak from the $m_i$9 spectrum and backgrounds (ranging from $10^{25}$00 FWHM in bolometers/HPGe to $10^{25}$01 in liquid scintillator detectors)
• Scaling to large isotope mass: hundreds of kg to tonne scale to reach inverted-hierarchy sensitivity (Cardani, 2018, Dolinski et al., 2019, Brugnera, 17 Jan 2025)

6. Implications and Future Prospects

A positive observation of $10^{25}$02 would establish:

• Lepton-number violation and thus a breakdown of Standard Model accidental symmetries
• The Majorana nature of neutrinos, confirming that neutrino mass arises at least partially via Majorana terms
• The absolute neutrino mass scale, and, given sufficient precision, provide constraints or measurement of the Majorana phases and the ordering of neutrino masses
• Evidence for B–L violation, with direct links to baryogenesis scenarios via leptogenesis (Cardani, 2018, Dolinski et al., 2019)

Conversely, null results at sensitivities corresponding to $10^{25}$03 meV would disfavor standard inverted ordering under the light-neutrino exchange scenario. They would also place stringent constraints on models of non-standard lepton-number-violating physics, such as TeV-scale left–right symmetric models or R-parity–violating SUSY (Banerjee et al., 9 Oct 2025).

Next-generation experiments—such as LEGEND-1000 ($10^{25}$04Ge), nEXO ($10^{25}$05Xe), CUPID ($10^{25}$06Te, $10^{25}$07Mo), SNO+ (high-loading $10^{25}$08Te)—aim to achieve sensitivities sufficient to probe the full inverted-hierarchy parameter space and, with further scaling and theoretical improvements in NME calculation, even approach the normal-hierarchy regime (Brugnera, 17 Jan 2025, Paton, 2019, Cardani, 2018, Grebe, 1 Apr 2025).

Ultimate interpretation will require:

• Multi-isotope, multi-technology confirmation
• Advances in nuclear theory to reduce NME uncertainties below $10^{25}$0910–20%
• Complementary information from cosmological sum-of-mass limits, single $10^{25}$10-decay experiments (e.g., KATRIN), and accelerator-based LNV searches (Grebe, 1 Apr 2025, Cardani, 2018, Rodejohann, 2010)

In summary, the search for 0νββ is entering a precision era in which meaningful conclusions regarding neutrino mass, new physics, and fundamental symmetries depend not only on experimental reach but also on detailed control of nuclear structure and theoretical interpretation. A discovery—or continued exclusion—will have far-reaching consequences for particle physics and cosmology.