Uncertainty Quantification of the $^{76}$Ge Neutrinoless Double-Beta Decay Nuclear Matrix Element

This paper quantifies the theoretical uncertainty of the neutrinoless double-beta decay nuclear matrix element for $^{76}$Ge by applying a rigorous statistical protocol with bounded fluctuations to effective interactions and Bayesian Model Averaging, yielding a central value of 2.46 with a standard deviation of 0.25.

The Gist

Imagine the universe is a giant, complex puzzle, and one of the most mysterious pieces is the neutrino. Scientists suspect that neutrinos might be their own antiparticles (like a mirror image that is actually the same person). To prove this, they are looking for a very rare event called neutrinoless double-beta decay. It's like watching two people in a room suddenly swap places without anyone else entering or leaving—a violation of the usual rules of physics.

The paper you provided is about 76Germanium (76Ge), a specific type of atom that is a prime candidate for this experiment. However, there's a problem: while the experiments are getting better at looking for this decay, the math used to predict how likely it is to happen is full of guesswork.

Here is a simple breakdown of what the authors did, using some everyday analogies:

1. The Problem: The "Recipe" is Uncertain

Think of the atom (76Ge) as a complex cake. To predict how the cake will taste (or in this case, how likely the decay is), scientists use a "recipe" called a Nuclear Matrix Element (NME).

• The Issue: Different scientists have slightly different recipes. Some say the cake will be light and fluffy; others say it will be dense and heavy. Because we don't know which recipe is perfect, we don't know how to interpret the experimental results. If the experiment says "we didn't find it," is it because the decay doesn't exist, or because our recipe was wrong?

2. The Solution: The "Taste-Test" Simulation

Instead of guessing which recipe is right, the authors decided to run a massive simulation.

• The Analogy: Imagine you have three master bakers (three different mathematical models called Hamiltonians: JUN45, GCN2850, and JJ44b). Instead of just baking one cake with each, they decided to bake 200 slightly different versions of each cake.
• The Method: They took the original recipes and made tiny, random adjustments to the ingredients (the "two-body matrix elements"). They changed the amounts by about 10%—enough to see how sensitive the cake is to a pinch of salt or a splash of milk, but not enough to ruin the cake entirely.
• The Goal: They baked thousands of these "what-if" cakes to see how much the final result (the NME) wiggles around. This creates a safety margin or a "confidence zone" for the answer.

3. The Results: Finding the Sweet Spot

After running all these simulations, they looked at the data:

• The Average: They found that the most likely value for the NME is 2.46.
• The Uncertainty: They calculated that the answer is likely between 2.21 and 2.71 (give or take 0.25).
• The "Vibe Check": They didn't just look at the decay number. They also checked other things about the atom, like how much energy it takes to make it vibrate (excitation energies) or how it spins. They found that if the "recipe" predicts the decay rate correctly, it also predicts these other physical properties correctly. It's like checking if a cake rises properly; if it does, you can trust the recipe.

4. The Conclusion: A Better Map for the Future

The authors combined their three different bakeries into one super-recipe using a statistical method called "Bayesian Model Averaging."

• What this means: They didn't pick just one winner. Instead, they blended the three best guesses together to create a single, highly reliable probability map.
• Why it matters: This map tells experimentalists (the people building the detectors) exactly how much "wiggle room" they have in their calculations. It stops them from panicking if their numbers don't match a single, rigid prediction.

Summary

In short, this paper is like a quality control audit for the math used to hunt for neutrinoless double-beta decay. The authors didn't discover the decay itself; instead, they built a statistical safety net. They showed that even if we tweak the ingredients of our nuclear models slightly, the answer stays surprisingly stable. This gives scientists a much clearer, more honest picture of where they stand in the search for new physics.

Technical Summary

Technical Summary: Uncertainty Quantification of the $^{76}$Ge Neutrinoless Double-Beta Decay Nuclear Matrix Element

Problem Statement
The search for neutrinoless double-beta decay ($0\nu\beta\beta$) is a primary avenue for probing lepton-number violation and physics beyond the Standard Model. While experimental collaborations like GERDA and LEGEND have established increasingly stringent half-life limits for the $^{76}$Ge isotope, the extraction of the effective Majorana neutrino mass is currently limited by theoretical uncertainties. Specifically, the nuclear matrix element (NME), which encapsulates the complex many-body dynamics of the parent and daughter nuclei, represents the dominant source of theoretical error. Existing calculations using various nuclear structure methodologies (e.g., pn-QRPA, IBA, ISM) yield a broad spread of NME values. While the Interacting Shell Model (ISM) is noted for its systematic treatment of configuration mixing and minimal reliance on phenomenological tuning, a rigorous quantification of its intrinsic theoretical uncertainty for $^{76}$Ge within a specific valence space has been required to guide experimental interpretation.

Methodology
The authors adapt a rigorous statistical protocol previously established for $^{48}$Ca and $^{136}$Xe to the $^{76}$Ge system. The study utilizes the $jj44$ (or $f_5pg_9$) valence configuration built upon a $^{56}$Ni inert core, with active nucleons occupying the $1p_{3/2}$, $1p_{1/2}$, $0f_{5/2}$, and $0g_{9/2}$ orbitals.

The core of the methodology involves:

1. Ensemble Generation: Starting from three established effective Hamiltonians (JUN45, GCN2850, and JJ44b), the authors generate perturbed ensembles. The two-body matrix elements (TBME) of these interactions are varied uniformly within a $\pm 10\%$ window around their original values. This amplitude is chosen to remain within empirically motivated bounds while avoiding unphysical overfitting. Single-particle energies (SPE) are held fixed.
2. Sampling: For each of the three baseline interactions, 200 perturbed Hamiltonians are generated. Full diagonalizations are performed without model-space truncation to capture all relevant configuration mixing.
3. Observable Calculation: For each perturbed Hamiltonian, a suite of 12 low-energy observables is computed. These include the $0\nu\beta\beta$ NME ($M^{0\nu}$), the two-neutrino double-beta decay NME ($M^{2\nu}$), Gamow-Teller transition probabilities ($P_{GT}$), $B(E2)\uparrow$ transition strengths, and excitation energies for the $2^+$, $4^+$, and $6^+$ states in both $^{76}$Ge and $^{76}$Se.
   • For $2\nu\beta\beta$ and Gamow-Teller transitions, a phenomenological quenching factor $q=0.65$ is applied.
   • The $0\nu\beta\beta$ NME is evaluated without quenching using the closure approximation and Jastrow short-range correlations.
4. Bayesian Model Averaging: The results from the three Hamiltonian ensembles are synthesized into a unified probability distribution for $M^{0\nu}$. While initial evidence integrals favored GCN2850, a compromise weighting scheme ($W_{GCN2850} = 4/6$, $W_{JUN45} = W_{JJ44b} = 1/6$) is adopted to avoid over-reliance on a single interaction.

Key Results

• NME Distribution: The analysis yields a constrained probability distribution for the $^{76}$Ge $0\nu\beta\beta$ NME with a central value of 2.46 and a standard deviation of 0.25. This corresponds to a 90% confidence interval of $[1.89, 3.07]$.
• Stability: The ensemble averages closely track the unperturbed Hamiltonian predictions, and standard deviations remain modest across all observables. This indicates an absence of pathological sensitivity to specific TBME fluctuations and confirms the robustness of shell-model predictions in the $jj44$ space.
• Correlations: A comprehensive correlation analysis reveals:
  • A robust positive correlation ($\rho > 0.8$) between $M^{0\nu}$ and $M^{2\nu}$, despite their distinct operator structures.
  • Strong interdependencies among excited-state energies and between NMEs and quadrupole transition strengths.
  • Notable anti-correlations involving parent $B(E2)\uparrow$ values, reflecting underlying structural trade-offs in configuration mixing.
• Hamiltonian Dependence: The weighted sum of the distributions closely mirrors the GCN2850 kernel density estimate. The authors note that within the $jj44$ space, the choice of baseline interaction exerts less influence on NME dispersion than in other model spaces, contrasting sharply with the parameter sensitivity often observed in pn-QRPA approaches.
• Limitations: The study acknowledges that canonical effective charges do not fully reproduce experimental $B(E2)$ values, a known limitation of the $jj44$ space due to the absence of $f_{7/2}$ and $g_{7/2}$ orbitals which limits the capture of quadrupole collectivity. However, the authors verify that this does not alter the correlation structures or the NME conclusions.

Significance and Claims
The paper claims to provide the first rigorous quantification of the intrinsic theoretical spread for the $^{76}$Ge NME within the interacting shell model approach using the $jj44$ configuration. By integrating simulated variations into a Bayesian framework and benchmarking against empirical spectroscopic data, the authors establish a constrained probability distribution that can serve as a rigorous theoretical prior for extracting neutrino mass limits from experimental half-life bounds.

The authors emphasize that this distribution allows for more accurate isotope-mass scaling and resource allocation for next-generation detectors. Furthermore, the methodology offers a benchmark for future ab initio and many-body methods, potentially guiding the refinement of effective interactions and the development of consistent effective operators derived from first principles. The work bridges the gap between nuclear structure theory and precision neutrino physics by providing a controlled statistical ensemble to quantify theoretical uncertainties, rather than relying on ad hoc error bars.