pith. sign in

arxiv: 2605.21657 · v1 · pith:NV3YTPGQnew · submitted 2026-05-20 · ⚛️ nucl-th

Uncertainty Quantification of the ⁷⁶Ge Neutrinoless Double-Beta Decay Nuclear Matrix Element

Pith reviewed 2026-05-22 08:24 UTC · model grok-4.3

classification ⚛️ nucl-th
keywords neutrinoless double-beta decaynuclear matrix element76GeBayesian model averagingshell modeluncertainty quantificationeffective nuclear interaction
0
0 comments X

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.

This paper quantifies the theoretical uncertainty in the nuclear matrix element for neutrinoless double-beta decay in 76Ge. It does so by applying bounded fluctuations to two-body matrix elements of an effective interaction and using Bayesian model averaging to fit low-energy spectroscopic data. The resulting probability distribution has a central value of 2.46 and standard deviation 0.25. This constrained value helps in interpreting experimental limits on the decay half-life and probing physics beyond the Standard Model.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.21657 by Andrei Neacsu, Mihai Horoi.

Figure 1
Figure 1. Figure 1: The heatmap for all 12 observables when using the GCN2850 effective Hamiltonian. Very similar representations are [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Correlation matrix for observables that have correlation factor greater than 0.5, when using the GCN2850 Hamiltonian. [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Distributions based on experimental data (in red) compared with the KDE (in blue) obtained from the GCN2850 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Similar to Figure [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Similar to Figures 3 and 4, the distributions based on experimental data (in red) compared with the KDE (in blue) obtained from the JJ44b starting Hamiltonian [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: PDFs of the 0νββ NME distributions for the JUN45, jj44b and GCN2850 Hamiltonians and their weighted sum (red curve, see text for details). The substantial overlap among all distributions con￾firms that shell-model calculations yield highly consistent 0νββ NME predictions in the jj44 space, irrespective of baseline interaction choice or minor optimization varia￾tions. This contrasts sharply with pn-QRPA res… view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [§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.
  2. [§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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

1 steps flagged

Self-cited protocol and prior effective interaction yield NME distribution with moderate dependence on fitted inputs

specific steps
  1. 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

1 free parameters · 1 axioms · 0 invented entities

Assessment limited to abstract; the central result rests on the transferability of the fluctuation protocol and the adequacy of the chosen valence space and effective interaction.

free parameters (1)
  • fluctuation bounds on two-body matrix elements
    The paper introduces systematic bounded fluctuations whose specific numerical ranges are not stated in the abstract but are required to generate the quoted NME distribution.
axioms (1)
  • domain assumption The valence configuration and effective interaction calibrated on 82Se are appropriate for 76Ge low-energy structure.
    Invoked when adapting the protocol from prior work without re-deriving the interaction.

pith-pipeline@v0.9.0 · 5825 in / 1459 out tokens · 43508 ms · 2026-05-22T08:24:48.936228+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

72 extracted references · 72 canonical work pages

  1. [1]

    Horoi, A

    M. Horoi, A. Neacsu, and S. Stoica, Phys. Rev. C106, 054302 (2022)

  2. [2]

    Horoi, A

    M. Horoi, A. Neacsu, and S. Stoica, Phys. Rev. C107, 045501 (2023)

  3. [3]

    Neacsu and M

    A. Neacsu and M. Horoi, Symmetry16, 974 (2024)

  4. [4]

    Fukuda, T

    Y. Fukuda, T. Hayakawa, E. Ichihara, K. Inoue, K. Ishi- hara, H. Ishino, Y. Itow, T. Kajita, J. Kameda, S. Ka- suga,et al., Phys. Rev. Lett.81, 1562 (1998)

  5. [5]

    Ahmad, R

    Q. Ahmad, R. Allen, T. Andersen, J. Anglin, G. B¨ uhler, J. Barton, E. Beier, M. Bercovitch, J. Bigu, S. Biller, et al., Phys. Rev. Lett.87, 071301 (2001)

  6. [6]

    Pontecorvo, Sov

    B. Pontecorvo, Sov. J. Exp. Theor. Phys.6, 429 (1958)

  7. [7]

    Pontecorvo, Sov

    B. Pontecorvo, Sov. J. Exp. Theor. Phys.26, 984 (1968)

  8. [8]

    Taroni, Nat

    A. Taroni, Nat. Phys.11, 891 (2015)

  9. [9]

    Barabash, Universe6, 159 (2020)

    A. Barabash, Universe6, 159 (2020)

  10. [10]

    Schechter and J

    J. Schechter and J. Valle, Phys. Rev. D25, 2951 (1982)

  11. [11]

    Hirsch, S

    M. Hirsch, S. Kovalenko, and I. Schmidt, Phys. Lett. B 642, 106 (2006)

  12. [12]

    Avignone, F.T., S

    I. Avignone, F.T., S. Elliott, and J. Engel, Rev. Mod. Phys.80, 481 (2008)

  13. [13]

    Vergados, H

    J. Vergados, H. Ejiri, and F. Simkovic, Rep. Prog. Phys. 75, 106301 (2012)

  14. [14]

    Rodejohann, J

    W. Rodejohann, J. Phys. G39, 124008 (2012)

  15. [15]

    Agostini, G

    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)

  16. [16]

    Arnquist, F

    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)

  17. [17]

    Anton, I

    G. Anton, I. Badhrees, P. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. Cao, W. Cen, et al., Phys. Rev. Lett.123, 161802 (2019)

  18. [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)

  19. [19]

    M. Doi, T. Kotani, H. Nishiura, and E. Takasugi, Progr. Theor. Exp. Phys.69, 602 (1983)

  20. [20]

    M. Doi, T. Kotani, and E. Takasugi, Prog. Theor. Phys. Suppl.83, 1 (1985)

  21. [21]

    Suhonen and O

    J. Suhonen and O. Civitarese, Phys. Rep.300, 123 (1998)

  22. [22]

    Nit ¸escu, S

    O. Nit ¸escu, S. Ghinescu, V. Sevestrean, M. Horoi, F. ˇSimkovic, and S. Stoica, arXiv (2024), 2402.13784

  23. [23]

    Nit ¸escu, S

    O. Nit ¸escu, S. Ghinescu, S. Stoica, and F.ˇSimkovic, Uni- verse10, 98 (2024)

  24. [24]

    Nit ¸escu, R

    O. Nit ¸escu, R. Dvornick´ y, and F.ˇSimkovic, Phys. Rev. C109, 025501 (2024)

  25. [25]

    Nit ¸escu, R

    O. Nit ¸escu, R. Dvornick´ y, S. Stoica, and F. ˇSimkovic, Universe7, 147 (2021)

  26. [26]

    J. Nabi, M. Ishfaq, O. Nit ¸escu, M. Mirea, and S. Stoica, Universe6, 5 (2019)

  27. [27]

    Mirea, T

    M. Mirea, T. Pahomi, and S. Stoica, Rom. Rep. Phys. 67, 872 (2015)

  28. [28]

    Stoica and M

    S. Stoica and M. Mirea, Phys. Rev. C88, 037303 (2013)

  29. [29]

    Kotila and F

    J. Kotila and F. Iachello, Phys. Rev. C85, 034316 (2012)

  30. [30]

    Caurier, A

    E. Caurier, A. Poves, and A. Zuker, Phys. Lett. B252, 13 (1990)

  31. [31]

    Caurier, F

    E. Caurier, F. Nowacki, A. Poves, and J. Retamosa, Phys. Rev. Lett.77, 1954 (1996)

  32. [32]

    Caurier, G

    E. Caurier, G. Martinez-Pinedo, F. Nowacki, A. Poves, and A. Zuker, Rev. Mod. Phys.77, 427 (2005)

  33. [33]

    Horoi, S

    M. Horoi, S. Stoica, and B. Brown, Phys. Rev. C75, 034303 (2007)

  34. [34]

    Horoi and S

    M. Horoi and S. Stoica, Phys. Rev. C81, 024321 (2010)

  35. [35]

    Horoi, Phys

    M. Horoi, Phys. Rev. C87, 014320 (2013)

  36. [36]

    Horoi and B

    M. Horoi and B. Brown, Phys. Rev. Lett.110, 222502 (2013)

  37. [37]

    Sen’kov and M

    R. Sen’kov and M. Horoi, Phys. Rev. C90, 051301(R) (2014)

  38. [38]

    Neacsu and M

    A. Neacsu and M. Horoi, Phys. Rev. C91, 024309 (2015)

  39. [39]

    Horoi and A

    M. Horoi and A. Neacsu, Phys. Rev. C93, 024308 (2016)

  40. [40]

    Horoi and A

    M. Horoi and A. Neacsu, Phys. Rev. C98, 035502 (2018)

  41. [41]

    Simkovic, G

    F. Simkovic, G. Pantis, J. Vergados, and A. Faessler, Phys. Rev. C60, 055502 (1999)

  42. [42]

    Stoica and H

    S. Stoica and H. Klapdor-Kleingrothaus, Nucl. Phys. A 694, 269 (2001)

  43. [43]

    Rodin, A

    V. Rodin, A. Faessler, F. Simkovic, and P. Vogel, Nucl. Phys. A766, 107 (2006)

  44. [44]

    Kortelainen and J

    M. Kortelainen and J. Suhonen, Phys. Rev. C75, 051303(R) (2007)

  45. [45]

    Faessler, V

    A. Faessler, V. Rodin, and F. Simkovic, J. Phys. G39, 124006 (2012)

  46. [46]

    Simkovic, V

    F. Simkovic, V. Rodin, A. Faessler, and P. Vogel, Phys. 10 Rev. C87, 045501 (2013)

  47. [47]

    Barea and F

    J. Barea and F. Iachello, Phys. Rev. C79, 044301 (2009)

  48. [48]

    Barea, J

    J. Barea, J. Kotila, and F. Iachello, Phys. Rev. C87, 014315 (2013)

  49. [49]

    Rodriguez and G

    T. Rodriguez and G. Martinez-Pinedo, Phys. Rev. Lett. 105, 252503 (2010)

  50. [50]

    P. Rath, R. Chandra, K. Chaturvedi, P. Lohani, P. Raina, and J. Hirsch, Phys. Rev. C88, 064322 (2013)

  51. [51]

    Novario, P

    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)

  52. [52]

    J. Yao, B. Bally, J. Engel, R. Wirth, T. Rodr´ ıguez, and H. Hergert, Phys. Rev. Lett.124, 232501 (2020)

  53. [53]

    Belley, C

    A. Belley, C. Payne, S. Stroberg, T. Miyagi, and J. Holt, Phys. Rev. Lett.126, 042502 (2021)

  54. [54]

    Engel and J

    J. Engel and J. Men´ endez, Rep. Prog. Phys.80, 046301 (2017)

  55. [55]

    Dolinski, A

    M. Dolinski, A. Poon, and W. Rodejohann, Annu. Rev. Nucl. Part. Sci.69, 219 (2019)

  56. [56]

    Retamosa, E

    J. Retamosa, E. Caurier, and F. Nowacki, Phys. Rev. C 51, 371 (1995)

  57. [57]

    Balysh, A

    A. Balysh, A. DeSilva, V. Lebedev, K. Lou, M. Moe, M. Nelson, A. Piepke, A. Pronskiy, M. Vient, and P. Vo- gel, Phys. Rev. Lett.77, 5186 (1996)

  58. [58]

    B. A. Brown and W. A. Richter, Phys. Rev. C74, 034315 (2006)

  59. [59]

    Honma, T

    M. Honma, T. Otsuka, T. Mizusaki, and M. Hjorth- Jensen, Phys. Rev. C80, 064323 (2009)

  60. [60]

    Men´ endez, A

    J. Men´ endez, A. Poves, E. Caurier, and F. Nowacki, Nucl. Phys. A818, 139 (2009)

  61. [61]

    Cheal, E

    B. Cheal, E. Man´ e, J. Billowes, M. Bissell, K. Blaum, B. Brown, F. Charlwood, K. Flanagan, D. Forest, C. Geppert,et al., Phys. Rev. Lett.104, 252502 (2010)

  62. [62]

    Horoi, Universe10, 252 (2024)

    M. Horoi, Universe10, 252 (2024)

  63. [63]

    Qi and Z

    C. Qi and Z. Xu, Phys. Rev. C86, 044323 (2012)

  64. [64]

    G. A. Miller and J. E. Spencer, Annals of Physics100, 562 (1976)

  65. [65]

    Zuker, J

    A. Zuker, J. Retamosa, A. Poves, and E. Caurier, Phys. Rev. C52, R1741 (1995)

  66. [66]

    Sen’kov, M

    R. Sen’kov, M. Horoi, and B. Brown, Phys. Rev. C89, 054304 (2014)

  67. [67]

    Kortelainen and J

    M. Kortelainen and J. Suhonen, Europhys. Lett.58, 666 (2002)

  68. [68]

    Kortelainen and J

    M. Kortelainen and J. Suhonen, Nucl. Phys. A713, 501 (2003)

  69. [69]

    Zinatulina, V

    D. Zinatulina, V. Brudanin, V. Egorov, C. Petitjean, M. Shirchenko, J. Suhonen, and I. Yutlandov, Phys. Rev. C99, 024327 (2019)

  70. [70]

    Bajpai, L

    D. Bajpai, L. Baudis, V. Belov, E. Bossio, T. Cocolios, H. Ejiri, E. Sushenok, M. Fomina, I. Hashim, M. Heines, et al., arXiv (2024), 2404.12686

  71. [71]

    Collaboration, N

    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. [72]

    Adhikari, S

    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)