REVIEW 2 major objections 6 minor 17 references
Search for the in-situ production of $^{77}$Ge in the GERDA neutrinoless double-beta decay experiment
T0 review · 2 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read The full GERDA Phase II data set shows zero 77Ge delayed-coincidence candidates, bounding the in-situ production rate below 0.216 nuclei per kilogram-year at 90% CL.
desk verdict Solid GERDA background measurement with a wrong abstract phrase; worth reviewing after a small fix. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing mechanism is the delayed-coincidence signature of the $^{77}$Ge decay chain. Beta decay of $^{77}$Ge ($T_{1/2}=11.21$ h) populates the $9/2^+$ isomeric state of $^{77}$As at 475 keV ($T_{1/2}=114$ µs) in about 33% of decays, and the delayed electromagnetic de-excitation emits gammas of 211, 264, or 475 keV. Because the 114 µs half-life is much longer than the roughly 1.5 µs charge-collection time of the germanium detectors, the delayed pulse arrives after the prompt pulse in the same detector and appears as pile-up on the waveform's exponential tail. New digital signal processing applies trapezoidal filters of length matched to the trigger separation, extracts the heights of both pulses, reconstructs the time difference, and is sensitive down to 4.5 µs separation; Monte Carlo simulation then fixes the combined selection efficiency at 32.4%. A unified small-signal counting procedure converts the observed zero candidates into the upper limit, with the branching-ratio uncertainties treated as nuisance parameters.
What would settle it
Reanalysing the GERDA Phase II waveforms with a time-difference threshold below 4.5 µs and a delayed-energy acceptance extended into the continuum: if this finds a cluster of candidates whose time differences follow the 114 µs exponential of the $^{77\mathrm{m}}$As isomer, the paper's $N_\mathrm{obs}=0$ result and its quoted production-rate limit would be contradicted; if no such cluster appears, the null is strengthened.
Extended reading notes
Core claim
The central result is a null observation with a sharp bound. In the full GERDA Phase II exposure of 103.7 kg·yr, after applying the delayed-coincidence selection built from the $^{77}$Ge decay chain, the observed number of candidates is $N_\mathrm{obs}=0$ against an expected random-coincidence background of $N_\mathrm{rc}=0.04$. A unified small-signal upper-limit calculation gives an upper limit of $<2.4$ delayed-coincidence counts at 90% CL, which translates into a $^{77}$Ge production rate of $<0.216$ nuc/(kg·yr) and, assuming equal production of the ground and isomeric states and a ground-state population probability of $(59.5\pm8.1)\%$, a combined $^{77}$Ge and $^{77\mathrm{m}}$Ge production rate of $<0.38$ nuc/(kg·yr). The paper reads this as corroboration of the earlier Monte Carlo estimate of $(0.21\pm0.07)$ nuc/(kg·yr), and a Bayesian update centers the rate at $0.18$ nuc/(kg·yr) with a 1$\sigma$ credibility interval of $[0.106, 0.251]$ nuc/(kg·yr). It further claims that tagging the $^{77\mathrm{m}}$As isomer can reject about 62% of the $^{77}$Ge background events that survive other cuts, and that this keeps the in-situ cosmogenic background in LEGEND-1000 below $10^{-6}$ cts/(keV·kg·yr), at most 10% of the experiment's background budget.
Load-bearing premise
The limit assumes the Monte Carlo simulation of the germanium detectors correctly predicts how often a real $^{77}$Ge decay chain would survive the selection, because the quoted $<0.5\%$ efficiency uncertainty covers only the statistical spread of the generated samples and not the full systematic accuracy of the detector simulation; if the simulated efficiency is too high, the production-rate limit is too strong.
Editorial extensions
If this is right
- The combined $^{77}$Ge and $^{77\mathrm{m}}$Ge in-situ production rate at the GERDA site is bounded below $0.38$ nuc/(kg·yr) at 90% CL, about a factor of two above the Monte Carlo prediction and roughly an order of magnitude tighter than the earlier GERDA-based bound of $<4.1$ nuc/(kg·yr).
- A Bayesian update of the simulation prior by the new likelihood centers the production rate at $0.18$ nuc/(kg·yr) with a 1$\sigma$ credibility interval of $[0.106, 0.251]$ nuc/(kg·yr), a scaling factor of $0.85^{+0.35}_{-0.34}$ relative to the original estimate.
- Tagging the $^{77\mathrm{m}}$As delayed de-excitation rejects about 62% of the $^{77}$Ge background events that survive the other analysis cuts, lowering the $^{77}$Ge/$^{77\mathrm{m}}$Ge background-index contribution to $(1.50\pm0.07\pm0.67)\times10^{-6}$ cts/(keV·kg·yr).
- For LEGEND-1000, the in-situ cosmogenic background contribution is estimated below $10^{-6}$ cts/(keV·kg·yr), no more than 10% of the target background index of $<10^{-5}$ cts/(keV·kg·yr).
- After one year of LEGEND-200 running at its final detector mass, the sensitivity of this delayed-coincidence search doubles, so the limit would tighten by a factor of two or a signal could appear.
Reading between the lines
- An editorial inference: because the expected random-coincidence rate is only 0.04 counts, applying the $^{77\mathrm{m}}$As tag as an active veto in LEGEND-200 costs almost no signal efficiency, so the same delayed-coincidence logic could be run online as a live background tag rather than offline.
- An editorial inference: the factor-of-two headroom between the simulation and the limit is small next to the simulation's 35% systematic uncertainty, so the agreement is weaker than it might appear; a larger exposure is needed to say whether the true rate sits near the prior or near the limit.
- An editorial inference: the eight observed continuum candidates with delayed energies between 200 and 500 keV do not follow the 114 µs lifetime, but a future analysis that includes the continuum rather than only the three gamma peaks could gain additional sensitivity.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a search for in-situ production of 77Ge in the full GERDA Phase II data set (103.7 kg·yr) using a delayed coincidence method that exploits the 114 µs isomeric state of 77As. No candidate events are found (N_obs = 0), with an expected random-coincidence background of N_rc = 0.04, leading to a Feldman-Cousins upper limit of <2.4 counts and a 77Ge production rate of <0.216 nuc/(kg·yr) at 90% CL. Using a neutron-capture model with a ground-state population fraction ε_g = (59.5±8.1)%, the total 77(m)Ge production rate is derived as <0.38 nuc/(kg·yr). A Bayesian update of a previous Monte Carlo prediction gives a posterior central value of 0.18 nuc/(kg·yr). The paper further estimates that the 77(m)Ge background contribution in LEGEND-1000 would be sub-dominant, below 10^-6 cts/(keV·kg·yr).
Significance. If the result holds, it is the strongest experimental constraint on in-situ 77(m)Ge production at LNGS and an important validation of cosmogenic background models for LEGEND-1000. The analysis uses the full GERDA Phase II exposure, a novel digital signal processing approach for pile-up reconstruction, and a carefully described delayed-coincidence selection with cross-checks against calibration data. The null observation, the internal consistency of the limit derivation, and the explicit treatment of nuisance parameters are strengths. The main concerns are a misleading statement in the abstract and the absence of a systematic uncertainty for the simulation-based efficiencies; both are addressable.
major comments (2)
- [Abstract vs. Section 3] The abstract states that the total 77(m)Ge production rate limit of <0.38 nuc/(kg·yr) is derived "assuming equal production rates". However, Section 3 derives this number by converting the 77Ge ground-state limit using ε_g = (ε_d + (1−ε_d)·ε_IT) = (59.5±8.1)%, with ε_d = (50±10)% and ε_IT = (19±2)%. A genuine equal-production assumption (ε_g = 0.5) would give a total limit of <0.216/0.5 = <0.43 nuc/(kg·yr), not <0.38. The quoted 0.38 follows from the neutron-capture model, not from equal production rates. This misstates the basis of the headline result and should be corrected in the abstract, for example by stating the adopted model assumption explicitly.
- [Section 2.3, Table 1] The uncertainty quoted for the selection efficiencies (<0.5%) is the statistical precision of the generated Monte Carlo samples, not a systematic uncertainty of the MaGe/GERDA detector simulation. The total efficiency directly enters the conversion from counts to production rate; an unquantified systematic bias in the simulated energy/multiplicity, time-difference, or pile-up acceptance would propagate directly into the quoted limit. The paper should either quantify systematic uncertainties of the simulation-based efficiencies (e.g., by varying simulation parameters or comparing with calibration data) or explicitly discuss the expected size of such effects and justify that they are negligible relative to the statistical precision of the null measurement.
minor comments (6)
- [Abstract] The sentence "an upper limit on the production rate of was set at <0.216 nuc/(kg·yr)" is missing the isotope; it should refer to 77Ge.
- [Figures 5 and 6] The axis label "dT [s]" appears in both figures, but the plotted values (0–20 and 0–800) are in microseconds, not seconds. Please correct the units to µs.
- [Section 2.3 and Figure 6 caption] The text states that the FWTM contains 96% of the peak area, while the caption of Figure 6 states 97%. These values should be reconciled.
- [Figure 7] The x-axis label contains a typo: "produdction" should be "production".
- [Introduction] The unit "nuc" (nuclei) is not defined at first use; please define it explicitly.
- [Section 3] The statement that the previous prediction is "further corroborated" should be qualified, since the 90% CL upper limit of 0.38 nuc/(kg·yr) is about a factor of two above the predicted central value of 0.21 nuc/(kg·yr); the result is consistent with but does not strongly confirm the prediction.
Circularity Check
No circular derivation: the primary 77Ge upper limit comes directly from GERDA data and is independent of the previous MC prediction; only the secondary total-rate conversion relies on a same-collaboration epsilon_g parameter, and the abstract's 'equal production rates' wording is internally inconsistent rather than circular.
full rationale
The primary upper limit is not circular. The paper records N_obs = 0 candidate delayed coincidences with N_rc = 0.04 expected random coincidences, applies Feldman-Cousins to obtain <2.4 counts, and converts to a 77Ge production rate using independently documented decay-scheme branching (33±1)% from ENSDF, exposure, and simulation-based selection efficiencies (total 32.4%) that are checked against calibration data. No parameter of the previous Monte Carlo prediction is fitted into this limit. The total 77(m)Ge rate is a secondary conversion: it divides the ground-state limit by epsilon_g = (epsilon_d + (1 - epsilon_d) * epsilon_IT) = (59.5±8.1)%, where epsilon_d is taken from the collaboration's earlier work [4] and from statistical-model calculations [15]. This is a same-group citation, but it is not a circular reduction because the measured 77Ge limit is not derived from epsilon_g and the paper does not define epsilon_g in terms of the total rate it is used to produce. The Bayesian update is transparently a prior (the earlier simulation) updated by the new likelihood, so it is disclosed model combination rather than hidden reuse. The abstract's phrase 'assuming equal production rates' conflicts with the body's epsilon_g = 59.5% ground-state fraction; a literal 50% equal-rate assumption would give about <0.43 nuc/(kg·yr) rather than <0.38. That is an internal inconsistency to be corrected, not a circularity. The LEGEND-1000 projection is explicitly a scaled estimate from Refs. [4,17] and is not presented as an independent measurement.
Assumptions & free parameters
free parameters (1)
- epsilon_d, direct ground-state population fraction after neutron capture on 76Ge =
(50 ± 10)%
assumptions (4)
- domain assumption MaGe/Geant4 simulation accurately models the GERDA detector response, including energy deposition, multiplicity, trigger thresholds, and pile-up reconstruction.
- domain assumption Cosmogenic neutron capture on 76Ge is the only significant in-situ source of 77Ge; radiogenic neutrons from outside are absorbed by the water shield.
- standard math Nuclear decay data (77Ge half-life, 77mAs half-life and energies, branching ratios) from ENSDF are accurate.
- standard math Feldman-Cousins and Cousins-Highland prescriptions correctly handle the Poisson upper limit with nuisance parameters.
Cite this review
Pith. "Pith review of Search for the in-situ production of $^{77}$Ge in the GERDA neutrinoless double-beta decay experiment." pith.science (2026). https://pith.science/paper/AU6KZDBW
@misc{pith2026250616955,
author = {Pith},
title = {Pith review of: Search for the in-situ production of $^77$Ge in the GERDA neutrinoless double-beta decay experiment},
year = {2026},
howpublished = {\url{https://pith.science/paper/AU6KZDBW}},
note = {Machine review of arXiv:2506.16955}
}
abstract
The beta decay of $^{77}$Ge and $^{77\mathrm{m}}$Ge, both produced by neutron capture on $^{76}$Ge, is a potential background for Germanium based neutrinoless double-beta decay search experiments such as GERDA or the LEGEND experiment. In this work we present a search for $^{77}$Ge decays in the full GERDA Phase II data set. A delayed coincidence method was employed to identify the decay of $^{77}$Ge via the isomeric state of $^{77}$As (9/2$^+$, 475 keV, ${T_{1/2} = 114}\,\mu $s, $^{77\mathrm{m}}$As). New digital signal processing methods were employed to select and analyze pile-up signals. No signal was observed, and an upper limit on the production rate of was set at $<0.216$ nuc/(kg$\cdot$yr) (90% CL). This corresponds to a total production rate of $^{77}$Ge and $^{77\mathrm{m}}$Ge of $<0.38$ nuc/(kg$\cdot$ yr) (90% CL), assuming equal production rates. A previous Monte Carlo study predicted a value for in-situ $^{77}$Ge and $^{77\mathrm{m}}$Ge production of (0.21$\pm$0.07) nuc/(kg$\cdot$yr), a prediction that is now further corroborated by our experimental limit. Moreover, tagging the isomeric state of $^{77\mathrm{m}}$As can be utilised to further suppress the $^{77}$Ge background. Considering the similar experimental configurations of LEGEND-1000 and GERDA, the cosmogenic background in LEGEND-1000 at LNGS is estimated to remain at a sub-dominant level.
Reference graph
Works this paper leans on
-
[1]
Agostiniet al.(GerdaCollab.), Upgrade for Phase II of theGerdaexperiment , Eur
M. Agostiniet al.(GerdaCollab.), Upgrade for Phase II of theGerdaexperiment , Eur. Phys. J. C 78, 388 (2018) https://doi.org/10.1140/epjc/s10052-018-5812-2
-
[2]
M. Agostiniet al.(GerdaCollab.), Final Results ofGerda on the Search for Neutrinoless Double-βDecay, Phys. Rev. Lett. 125, 252502 (2020)https://doi.org/10.1103/ PhysRevLett.125.252502
work page 2020
-
[3]
L. Pandolaet al., Monte Carlo evaluation of the muon- induced background in theGerdadouble beta decay ex- periment, Nucl. Instrum. Methods A 570, 149–158 (2007) https://doi.org/10.1016/j.nima.2006.10.103
-
[4]
C. Wiesinger, L. Pandola and S. Sch¨ onert, Virtual depth by active background suppression: Revisiting the cos- mic muon induced background ofGerdaPhase II, Eur. Phys. J. C 78, 597 (2018)https://doi.org/10.1140/epjc/ s10052-018-6079-3
doi:10.1140/epjc/ 2018
-
[5]
N. Abgrallet al.(LEGEND Collab.), The Large En- riched Germanium Experiment for NeutrinolessββDecay: LEGEND-1000 Preconceptual Design Report, arXiv (2021) https://doi.org/10.48550/arXiv.2107.11462
-
[6]
I. J. Arnquistet al., (MajoranaCollab.), Signatures of muonic activation in theMajoranaDemonstrator, Phys. Rev. C 105, 014617 (2022)https://doi.org/10.1103/ PhysRevC.105.014617
work page 2022
-
[7]
From ENSDF database as of February 21st, 2024. Ver- sion available athttp://www.nndc.bnl.gov/ensarchivals/ Full Evaluation Balraj Singh ENSDF 30-Sep-2020
work page 2024
-
[8]
M. Boswellet al., MaGe – a Geant4-based Monte Carlo application framework for low-background germanium ex- periments. IEEE Trans. Nucl. Sci. 58, 1212 (2011)https: //doi.org/10.1109/TNS.2011.2144619
arXiv 2011
Show all 17 references
-
[9]
A. Lazzaro, Signal processing and event classification for a background free neutrinoless double beta decay search with theGerdaexperiment, PhD Thesis, Technische Universit¨ at M¨ unchen (2019)https://mediatum.ub.tum.de/1507626
2019
-
[10]
Agostiniet al., The MGDO software library for data analysis in Ge neutrinoless double-beta decay experiments, J
M. Agostiniet al., The MGDO software library for data analysis in Ge neutrinoless double-beta decay experiments, J. Phys. Conf. Ser. 375 042027 (2012)https://doi.org/10. 1088/1742-6596/375/1/042027
2012
-
[11]
Agostiniet al.(GerdaCollab.), Calibration of the Gerdaexperiment, Eur
M. Agostiniet al.(GerdaCollab.), Calibration of the Gerdaexperiment, Eur. J. Phys. C 81 (2021) 682https: //doi.org/10.1140/epjc/s10052-021-09403-2
2021 doi
-
[12]
G. J. Feldman and R. D. Cousins, A Unified approach to the classical statistical analysis of small signals, Phys. Rev. D 57 (1998) 3873-3889https://doi.org/10.1103/PhysRevD. 57.3873 11
1998 doi
-
[13]
R. D. Cousins and V. Highland, Incorporating systematic uncertainties into an upper limit, Nucl. Instrum. Meth. A 320 (1992) 331-33https://doi.org/10.1016/0168-9002(92) 90794-5
1992 doi
-
[14]
R. L. Workmanet al.(Particle Data Group), Prog. Theor. Exp. Phys. 2022, 083C01 (2022)
2022
-
[15]
Grabmayr, Cross sections and gamma cascades in 77Ge needed for background reduction in 0νββexperi- ments on 76Ge, Eur
P. Grabmayr, Cross sections and gamma cascades in 77Ge needed for background reduction in 0νββexperi- ments on 76Ge, Eur. Phys. J. A (2024) 60, 115,https: //doi.org/10.1140/epja/s10050-024-01336-0
2024 doi
-
[16]
L. Vanhoefer, Limitations of Rare Event HPGe Exper- iments due to Muon-Induced Neutron Background, PhD Thesis, Technische Universit¨ at M¨ unchen (2018)https:// mediatum.ub.tum.de/1446626
2018
-
[17]
Neuberger, M
M. Neuberger, M. Morellaet al., Strategies for cosmo- genic 77(m)Ge reduction for LEGEND-1000 Experiment, Poster, Neutrino 2022https://doi.org/10.5281/zenodo. 6804443
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.