Neutrinoless Double Beta Decay Half-Life

Introduction

Neutrinos are among the most mysterious particles in modern physics. Experiments have already shown that they have mass, but an even deeper question remains open: are neutrinos distinct from their antiparticles, or are they actually the same object viewed in two roles? The second possibility is called the Majorana picture, and one of the clearest experimental tests of that idea is the search for neutrinoless double beta decay, often written as 0νββ. This calculator is designed to estimate the half-life associated with that hypothetical process when you supply three standard ingredients used in the literature: an effective Majorana mass, a phase-space factor, and a nuclear matrix element.

In ordinary double beta decay with two neutrinos, a nucleus changes charge by two units and emits two electrons together with two antineutrinos. That process is allowed in the Standard Model and has been measured in several isotopes. Neutrinoless double beta decay would be far more dramatic. In that case, the nucleus would still emit two electrons, but no neutrinos would appear in the final state. Symbolically, the transition is written in plain form as (Z + 2, A) -> (Z, A) + 2e-. If such a decay were observed, it would imply lepton-number violation and would strongly support the conclusion that neutrinos are Majorana particles.

The reason this matters is not just classification. A confirmed 0νββ signal would connect nuclear measurements to some of the biggest questions in particle physics and cosmology. It would provide evidence about the origin of neutrino mass, help constrain the neutrino mass ordering, and inform theories that try to explain why the universe contains more matter than antimatter. Because the expected half-lives are extraordinarily long, often beyond 10²⁶ years, even rough estimates are useful for understanding what present and future experiments can realistically probe.

This page therefore serves two purposes at once. First, it gives you a working calculator that turns common theoretical inputs into a predicted half-life. Second, it explains what those inputs mean in plain language so the result is easier to interpret. The tool is not a substitute for a full nuclear-structure calculation or an experimental sensitivity study, but it is a practical way to explore how strongly the predicted half-life changes when the assumed neutrino mass scale or nuclear parameters are varied.

Background and Physical Meaning

The standard factorized description of neutrinoless double beta decay separates the problem into three pieces: kinematics, nuclear structure, and particle physics. The kinematic part is summarized by the phase-space factor, usually denoted $G_{0\nu }$. This quantity depends on the isotope and on the available decay energy, and it carries units of inverse time. The nuclear-structure part is encoded in the nuclear matrix element $\left|M_{0\nu }\right|$, which summarizes how the initial and final nuclei are connected by the underlying transition operator. The particle-physics part enters through the effective Majorana mass $m_{\beta \beta }$, which combines neutrino masses, mixing angles, and Majorana phases.

Because these ingredients come from different areas of physics, the same half-life prediction can be discussed by nuclear theorists, neutrino phenomenologists, and experimental collaborations using slightly different conventions. This calculator follows the common light-Majorana-neutrino exchange expression and assumes that the values you enter are already in compatible units. In particular, the effective mass is entered in electronvolts, the phase-space factor is entered in inverse years, and the matrix element is treated as a dimensionless magnitude. The electron mass is built into the calculation as a fixed conversion scale.

One useful feature of the formula is that it makes the scaling behavior very clear. If the effective Majorana mass becomes smaller, the half-life grows rapidly. If the nuclear matrix element is larger, the decay is easier to observe and the half-life becomes shorter. If the phase-space factor is larger, the decay is also faster. This means that isotope choice matters, nuclear-model choice matters, and assumptions about neutrino masses matter. The calculator lets you see those dependencies directly instead of treating the half-life as an abstract number quoted in a paper.

How to Use

To use the calculator, enter one value in each of the three fields and then press the compute button. The first field is the effective Majorana mass $m_{\beta \beta }$ in eV. This is the particle-physics input and is often the quantity of greatest interest when people compare 0νββ results with neutrino-oscillation data. The second field is the phase-space factor $G_{0\nu }$ in yr⁻¹. You would normally take this from published tables for the isotope you care about. The third field is the magnitude of the nuclear matrix element $\left|M_{0\nu }\right|$, which is usually estimated from a nuclear model.

After you submit the form, the calculator returns a predicted half-life $T_{1/2}$ in years. The result is also accompanied by a short qualitative label indicating whether the predicted scale is currently testable, within next-generation sensitivity, or beyond that reach. That label is only a broad guide. It is not tied to one specific experiment, isotope, background model, or exposure assumption. It simply helps place the number in a familiar experimental context.

When choosing inputs, it is important to keep units consistent. The calculator expects $m_{\beta \beta }$ in eV, not meV or GeV. Likewise, the phase-space factor should be entered directly in inverse years. If you are copying values from a paper, check whether the authors include factors such as the axial coupling convention or other normalization choices in the quoted matrix element or phase-space factor. The calculator does not attempt to reconcile differing conventions automatically; it simply evaluates the expression using the numbers you provide.

For exploratory work, a good approach is to hold the isotope-dependent quantities fixed and vary the effective mass. That quickly shows how strongly the half-life responds to changes in neutrino-scale assumptions. You can also do the reverse: keep the effective mass fixed and compare different matrix-element choices to see how theoretical nuclear uncertainty propagates into the predicted half-life. This is especially useful when reading papers that quote a range rather than a single preferred nuclear matrix element.

Formula

The inverse half-life for neutrinoless double beta decay can be factorized into three main components, reflecting phase space, nuclear structure, and particle physics contributions:

Formula: T1_/2^− = G 0_ν | M 0_ν | 2(m_ββ / m_e) 2

${T{1}_{/2}}^{-}=G{0}_{\nu }\left|M{0}_{\nu }\right|{}^2\left(\frac{m_{\mathrm{\beta \beta }}}{m_e}\right){}^2$

Here, $G_{0\nu }$ is the phase-space factor with units of inverse time, computable from kinematics and Coulomb effects; $\left|M_{0\nu }\right|$ is the nuclear matrix element; $m_{\beta \beta }$ is the effective Majorana mass; and $m_e$ is the electron mass. The formula shows immediately that the decay rate is quadratic in both the matrix element and the effective mass ratio.

This calculator evaluates the rearranged form for the half-life itself:

Formula: T 1_/2 = 1 / (G 0_ν |^M (m_ββ/m_e)^2)

$T{1}_{/2}=\frac{1}{G{0}_{\nu }{\left|}^M{\left(\frac{m_{\mathrm{\beta \beta }}}{m_e}\right)}^2}$

Since $m\_e=0.511\text{MeV}=5.11\times {10}^5\text{eV}$, expressing all masses in electronvolts makes the numerical evaluation straightforward. In the page script, the electron mass is stored as 5.11e5 eV. The code then computes the denominator G × M² × mββ² and divides me² by that quantity to obtain the half-life in years.

The scaling is worth emphasizing because it helps interpret the output. If you double $m_{\beta \beta }$, the predicted half-life becomes four times shorter. If you halve the matrix element, the half-life becomes four times longer. This quadratic sensitivity is one reason why both neutrino-parameter uncertainty and nuclear-theory uncertainty have such a large impact on 0νββ forecasts.

Example

Suppose you want a quick estimate for a representative isotope and choose the following values: effective Majorana mass 0.05 eV, phase-space factor 1×10⁻¹⁴ yr⁻¹, and nuclear matrix element magnitude 5. These are the default values already loaded in the form. The calculator first squares the matrix element and the effective mass, then multiplies by the phase-space factor. Using the built-in electron mass scale, it returns a half-life of about 2.1×10²⁶ years.

That result is physically informative even before you compare it with a specific experiment. A half-life near 10²⁶ years sits in the range that current leading searches have begun to probe for some isotopes. If you keep the same phase-space factor and matrix element but reduce the effective mass to 0.02 eV, the half-life grows to roughly 1.3×10²⁷ years. If instead you increase the effective mass to 0.10 eV, the half-life drops to about 5.3×10²⁵ years. Those changes illustrate the strong inverse-square dependence on the effective mass.

The table below summarizes a few representative parameter choices, assuming a typical nuclear matrix element magnitude of 5:

A worked example like this is useful because it turns a symbolic formula into an intuition-building estimate. It also shows why experimental sensitivity targets are often discussed in decades of half-life rather than in small percentage changes. Moving from 10²⁶ to 10²⁷ years is not a minor improvement; it can correspond to a substantial shift in the effective mass region being tested.

Interpreting the Result

The number produced by the calculator is a predicted half-life, not a probability that a decay will happen during a short laboratory run. A half-life of 10²⁶ years means the process is extraordinarily rare for any single nucleus. Experiments compensate by observing enormous numbers of candidate nuclei over long periods while suppressing backgrounds as much as possible. As a result, the practical meaning of the output is comparative: it tells you whether your chosen parameter set points toward a decay rate that is already under pressure from existing limits or one that would require more ambitious future detectors.

The qualitative label shown below the numerical result is intentionally simple. “Currently testable” means the predicted half-life is below about 10²⁶ years in this page’s classification. “Within next-generation sensitivity” covers the region below about 10²⁷ years. “Beyond next-generation reach” indicates even longer half-lives. These thresholds are broad educational markers rather than rigorous discovery or exclusion criteria. Real sensitivity depends on isotope mass, enrichment, live time, energy resolution, background index, and analysis method.

It is also important to remember that the same measured half-life can imply different effective masses if different nuclear matrix elements are assumed. That is why papers often quote a band of inferred $m_{\beta \beta }$ values rather than a single number. This calculator can help you see that dependence directly by rerunning the estimate with several plausible matrix-element choices.

Limitations and Assumptions

This calculator is intentionally simple, so it is important to understand what it does not include. The formula implemented here assumes the standard light-Majorana-neutrino exchange mechanism. Other mechanisms for neutrinoless double beta decay exist in the theoretical literature, including heavy-particle exchange or right-handed current contributions, and those can modify the relationship between half-life and the effective parameters. If you are studying a nonstandard mechanism, this tool should be treated only as a rough analogy, not as a definitive model.

The calculator also assumes that the phase-space factor and nuclear matrix element you enter are already appropriate for the isotope and convention you want to use. In practice, published values can differ because of choices involving the axial coupling, short-range correlations, nuclear model space, or operator treatment. The tool does not validate those conventions or convert between them. It simply applies the numerical expression exactly as entered.

Another limitation is that the result is purely theoretical and does not include detector effects. Experimental sensitivity is not determined by half-life alone. Background rates, energy resolution, isotope abundance, enrichment fraction, exposure, and systematic uncertainties all matter. A predicted half-life that looks “testable” in the calculator may still be difficult for a real experiment if backgrounds are high or the isotope mass is limited. Conversely, a very long half-life may become accessible with enough exposure and sufficiently low background.

Finally, the effective Majorana mass itself is not a directly measured input in most contexts. It is a derived quantity that depends on neutrino masses, mixing parameters, and unknown Majorana phases. Because of possible cancellations, especially in the normal ordering, small values of $m_{\beta \beta }$ can arise even when individual neutrino masses are not tiny. That means the calculator is best used as a transparent exploration tool: it shows how assumptions map into half-life predictions, but it does not resolve the underlying neutrino-parameter ambiguities by itself.

Even with those caveats, the calculator is useful for teaching, quick comparisons, and order-of-magnitude checks. It captures the central relationship linking neutrino physics to a measurable nuclear observable and makes the extreme scale of the problem easier to appreciate. Whether you are reading a paper, preparing a lecture, or comparing isotopes for intuition, the tool provides a compact way to connect the theory inputs to the half-life scale that experiments seek to test.