Uncertainty Quantification of the ⁷⁶Ge Neutrinoless Double-Beta Decay Nuclear Matrix Element
Pith reviewed 2026-05-22 08:24 UTC · model grok-4.3
The pith
The 0νββ nuclear matrix element for 76Ge is determined to be 2.46 with a standard deviation of 0.25 through Bayesian analysis.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By integrating simulated variations of two-body matrix elements into a Bayesian Model Averaging framework and benchmarking against empirical spectroscopic data, the analysis yields a constrained probability distribution for the 0νββ NME with central value 2.46 and standard deviation 0.25.
What carries the argument
Bayesian Model Averaging of systematic bounded fluctuations applied to the two-body matrix elements of the effective interaction in the interacting shell model.
If this is right
- The uncertainty in the NME is quantified at the level of 0.25 around the central value of 2.46.
- Correlation analysis reveals structural dependencies among nuclear observables that can benchmark effective interactions.
- This distribution provides a more reliable input for deriving limits on the effective Majorana neutrino mass from experimental half-life bounds.
- Similar statistical protocols can be extended to other candidate isotopes for 0νββ decay.
Where Pith is reading between the lines
- If the derived distribution holds, it suggests that shell-model calculations with this level of uncertainty can be compared across different isotopes to identify systematic trends.
- Future refinements of effective interactions could be guided by minimizing discrepancies in the correlated observables identified here.
- An independent ab initio calculation falling outside the 2.21-2.71 range would indicate missing physics in the fluctuation model.
Load-bearing premise
The assumption that systematic bounded fluctuations applied to the two-body matrix elements of an effective interaction previously used for 82Se, when propagated and averaged via Bayesian methods against low-energy spectroscopic data, fully capture the dominant theoretical uncertainty in the NME for 76Ge.
What would settle it
An independent calculation or future experimental constraint on the NME or related observables that yields a value outside the interval from 2.21 to 2.71 would falsify the reported probability distribution.
Figures
read the original abstract
The experimental pursuit of neutrinoless double-beta decay ($0\nu\beta\beta$) constitutes one of the most compelling avenues for probing lepton-number violation and exploring physics beyond the Standard Model. Within this landscape, $^{76}$Ge has consistently ranked among the most promising isotopes for current and next-generation bolometric and liquid-scintillator experiments, notably GERDA and LEGEND. In the present work, we adapt a rigorous statistical protocol previously established for $^{48}$Ca~\cite{Horoi-prc22} and $^{136}$Xe~\cite{Horoi-Xe-2023} to the $^{76}$Ge system, utilizing a valence configuration that aligns with our recent investigation of $^{82}$Se~\cite{Neacsu-Symmetry-2024}. Our methodology introduces systematic, bounded fluctuations to the two-body matrix elements of established effective interactions, subsequently monitoring how these perturbations propagate through a suite of low-energy nuclear observables. Special emphasis is placed on the $0\nu\beta\beta$ nuclear matrix element (NME), whose theoretical uncertainty currently dominates the interpretation of experimental half-life limits. By integrating these simulated variations into a Bayesian Model Averaging framework and benchmarking against empirical spectroscopic data, we derive a constrained probability distribution for the NME. The resulting analysis yields a central value of 2.46 with an associated standard deviation of 0.25, thereby quantifying the intrinsic theoretical spread within the interacting shell model approach. Furthermore, we perform a comprehensive correlation analysis across all computed observables to evaluate internal consistency, identify non-trivial structural dependencies, and establish benchmarks that may guide the refinement of future effective interactions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript adapts a Bayesian Model Averaging protocol, previously applied to 48Ca and 136Xe, to 76Ge. It introduces systematic bounded fluctuations to the two-body matrix elements of an effective interaction previously tuned for 82Se, propagates these through low-energy nuclear observables, and benchmarks against spectroscopic data to derive a probability distribution for the 0νββ nuclear matrix element with central value 2.46 and standard deviation 0.25.
Significance. If the fluctuations and averaging procedure are shown to capture the dominant sources of uncertainty, the result supplies a statistically grounded theoretical error bar on the NME that can be directly used to interpret half-life limits from GERDA and LEGEND. The accompanying correlation analysis across observables is a useful addition for identifying structural dependencies within the shell-model framework.
major comments (2)
- [§3.2] §3.2 (fluctuation protocol): the bounded variations are applied only to the TBMEs of the 82Se interaction while the valence-space truncation, single-particle energies, and short-range form of the 0νββ transition operator remain fixed; because the NME is sensitive to high-momentum components not constrained by low-energy spectroscopy, the quoted spread of 0.25 may not represent the full theoretical uncertainty.
- [§4] §4 (Bayesian averaging): the weighting of models by agreement with spectroscopic data assumes that the chosen low-energy observables sufficiently constrain the degrees of freedom relevant to the NME; no cross-validation against an independent set of observables or against results from a different model space is presented to test this assumption.
minor comments (2)
- [Eq. (7)] The definition of the fluctuation bounds (Eq. (7)) should be accompanied by a brief justification of why the chosen range is appropriate for 76Ge rather than simply inherited from the 82Se study.
- [Figure 5] Figure 5 (correlation matrix) would benefit from explicit labeling of which observables are spectroscopic and which are the NME itself to improve readability.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We respond to each major comment below, proposing revisions to clarify the scope of our uncertainty quantification.
read point-by-point responses
-
Referee: [§3.2] §3.2 (fluctuation protocol): the bounded variations are applied only to the TBMEs of the 82Se interaction while the valence-space truncation, single-particle energies, and short-range form of the 0νββ transition operator remain fixed; because the NME is sensitive to high-momentum components not constrained by low-energy spectroscopy, the quoted spread of 0.25 may not represent the full theoretical uncertainty.
Authors: We concur that the fluctuations are confined to the TBMEs of the 82Se interaction, with the valence-space truncation, single-particle energies, and the short-range form of the 0νββ operator held fixed. This protocol is intended to quantify the uncertainty stemming from variations in the effective two-body interaction while remaining consistent with the low-energy spectroscopic data used for benchmarking. Although the NME receives contributions from high-momentum components that are not fully constrained by low-energy observables, our approach provides a well-defined statistical measure of the spread within the interacting shell-model framework adopted here, in line with our earlier studies. We will revise §3.2 and the conclusions to explicitly delineate the boundaries of this uncertainty estimate and note that a more comprehensive assessment would require exploring additional model variations. revision: partial
-
Referee: [§4] §4 (Bayesian averaging): the weighting of models by agreement with spectroscopic data assumes that the chosen low-energy observables sufficiently constrain the degrees of freedom relevant to the NME; no cross-validation against an independent set of observables or against results from a different model space is presented to test this assumption.
Authors: The weighting in the Bayesian model averaging is based on the agreement with a standard set of low-energy spectroscopic observables, as successfully applied in our previous works on 48Ca and 136Xe. The correlation analysis presented in the manuscript demonstrates the internal consistency and highlights how variations in these observables relate to the NME. We recognize that an explicit cross-validation with a held-out set of observables or comparisons to calculations in alternative model spaces is not performed in this study. This would provide additional validation but lies beyond the present scope. We will add a brief discussion in §4 acknowledging this aspect of the methodology and referencing the supporting evidence from the correlation study. revision: partial
Circularity Check
Self-cited protocol and prior effective interaction yield NME distribution with moderate dependence on fitted inputs
specific steps
-
self citation load bearing
[Abstract]
"we adapt a rigorous statistical protocol previously established for ^{48}Ca~[Horoi-prc22] and ^{136}Xe~[Horoi-Xe-2023] to the ^{76}Ge system, utilizing a valence configuration that aligns with our recent investigation of ^{82}Se~[Neacsu-Symmetry-2024]"
The load-bearing methodology (systematic bounded TBME fluctuations + Bayesian Model Averaging against spectroscopic data) and the valence space are justified solely by citations to the authors' own prior papers rather than re-derived or independently validated in the present work; the final NME probability distribution therefore inherits its structure and bounds from those self-cited frameworks.
full rationale
The derivation applies Bayesian Model Averaging to bounded TBME fluctuations drawn from an effective interaction previously tuned for 82Se, with weights from low-energy spectroscopic data. The central NME value and spread (2.46 ± 0.25) are computed for 76Ge using this imported framework and chosen fluctuation bounds. While the specific numerical result for 76Ge is not identical to prior outputs by construction, the load-bearing statistical protocol, valence configuration, and interaction choice reduce to self-citations without re-derivation or external validation here. This creates moderate circularity but leaves independent computational content in the propagation and averaging steps.
Axiom & Free-Parameter Ledger
free parameters (1)
- fluctuation bounds on two-body matrix elements
axioms (1)
- domain assumption The valence configuration and effective interaction calibrated on 82Se are appropriate for 76Ge low-energy structure.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
introduces systematic, bounded fluctuations to the two-body matrix elements... Bayesian Model Averaging framework... central value of 2.46 with an associated standard deviation of 0.25
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
- [4]
- [5]
- [6]
- [7]
- [8]
- [9]
- [10]
- [11]
-
[12]
I. Avignone, F.T., S. Elliott, and J. Engel, Rev. Mod. Phys.80, 481 (2008)
work page 2008
-
[13]
J. Vergados, H. Ejiri, and F. Simkovic, Rep. Prog. Phys. 75, 106301 (2012)
work page 2012
- [14]
-
[15]
M. Agostini, G. Araujo, A. Bakalyarov, M. Balata, I. Barabanov, L. Baudis, C. Bauer, E. Bellotti, S. Be- logurov, A. Bettini,et al., Phys. Rev. Lett.125, 252502 (2020)
work page 2020
-
[16]
I. Arnquist, F. Avignone, A. Barabash, C. Barton, P. Barton, K. Bhimani, E. Blalock, B. Bos, M. Busch, M. Buuck,et al., Phys. Rev. Lett.130, 062501 (2023)
work page 2023
- [17]
-
[18]
S. Abe, S. Asami, M. Eizuka, S. Futagi, A. Gando, Y. Gando, T. Gima, A. Goto, T. Hachiya, K. Hata,et al., Phys. Rev. Lett.130, 051801 (2023)
work page 2023
-
[19]
M. Doi, T. Kotani, H. Nishiura, and E. Takasugi, Progr. Theor. Exp. Phys.69, 602 (1983)
work page 1983
-
[20]
M. Doi, T. Kotani, and E. Takasugi, Prog. Theor. Phys. Suppl.83, 1 (1985)
work page 1985
- [21]
-
[22]
O. Nit ¸escu, S. Ghinescu, V. Sevestrean, M. Horoi, F. ˇSimkovic, and S. Stoica, arXiv (2024), 2402.13784
-
[23]
O. Nit ¸escu, S. Ghinescu, S. Stoica, and F.ˇSimkovic, Uni- verse10, 98 (2024)
work page 2024
-
[24]
O. Nit ¸escu, R. Dvornick´ y, and F.ˇSimkovic, Phys. Rev. C109, 025501 (2024)
work page 2024
-
[25]
O. Nit ¸escu, R. Dvornick´ y, S. Stoica, and F. ˇSimkovic, Universe7, 147 (2021)
work page 2021
-
[26]
J. Nabi, M. Ishfaq, O. Nit ¸escu, M. Mirea, and S. Stoica, Universe6, 5 (2019)
work page 2019
- [27]
- [28]
- [29]
- [30]
-
[31]
E. Caurier, F. Nowacki, A. Poves, and J. Retamosa, Phys. Rev. Lett.77, 1954 (1996)
work page 1954
-
[32]
E. Caurier, G. Martinez-Pinedo, F. Nowacki, A. Poves, and A. Zuker, Rev. Mod. Phys.77, 427 (2005)
work page 2005
- [33]
- [34]
- [35]
- [36]
- [37]
- [38]
- [39]
- [40]
-
[41]
F. Simkovic, G. Pantis, J. Vergados, and A. Faessler, Phys. Rev. C60, 055502 (1999)
work page 1999
- [42]
- [43]
- [44]
- [45]
-
[46]
F. Simkovic, V. Rodin, A. Faessler, and P. Vogel, Phys. 10 Rev. C87, 045501 (2013)
work page 2013
- [47]
- [48]
-
[49]
T. Rodriguez and G. Martinez-Pinedo, Phys. Rev. Lett. 105, 252503 (2010)
work page 2010
-
[50]
P. Rath, R. Chandra, K. Chaturvedi, P. Lohani, P. Raina, and J. Hirsch, Phys. Rev. C88, 064322 (2013)
work page 2013
-
[51]
S. Novario, P. Gysbers, J. Engel, G. Hagen, G. Jansen, T. Morris, P. Navr´ atil, T. Papenbrock, and S. Quaglioni, Phys. Rev. Lett.126, 182502 (2021)
work page 2021
-
[52]
J. Yao, B. Bally, J. Engel, R. Wirth, T. Rodr´ ıguez, and H. Hergert, Phys. Rev. Lett.124, 232501 (2020)
work page 2020
- [53]
- [54]
-
[55]
M. Dolinski, A. Poon, and W. Rodejohann, Annu. Rev. Nucl. Part. Sci.69, 219 (2019)
work page 2019
- [56]
- [57]
-
[58]
B. A. Brown and W. A. Richter, Phys. Rev. C74, 034315 (2006)
work page 2006
- [59]
-
[60]
J. Men´ endez, A. Poves, E. Caurier, and F. Nowacki, Nucl. Phys. A818, 139 (2009)
work page 2009
- [61]
- [62]
- [63]
-
[64]
G. A. Miller and J. E. Spencer, Annals of Physics100, 562 (1976)
work page 1976
- [65]
- [66]
- [67]
- [68]
-
[69]
D. Zinatulina, V. Brudanin, V. Egorov, C. Petitjean, M. Shirchenko, J. Suhonen, and I. Yutlandov, Phys. Rev. C99, 024327 (2019)
work page 2019
- [70]
-
[71]
L. Collaboration, N. Abgrall, I. Abt, M. Agostini, A. Alexander, C. Andreoiu, G. Araujo, I. Avignone, F.T., W. Bae, A. Bakalyarov,et al., arXiv (2021), 2107.11462
-
[72]
G. Adhikari, S. Kharusi, E. Angelico, G. Anton, I. Arn- quist, I. Badhrees, J. Bane, V. Belov, E. Bernard, T. Bhatta,et al., J. Phys. G Nucl. Part. Phys.49, 015104 (2021)
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.