REVIEW 2 major objections 5 minor 1 cited by
Search for Double Beta Decay of $^{136}$Xe to the $0^+_1$ Excited State of $^{136}$Ba with PandaX-4T
T0 review · 2 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read The PandaX-4T experiment finds no evidence for the two-neutrino double beta decay of $^{136}$Xe to the first excited $0^+_1$ state of $^{136}$Ba, and sets a lower half-life limit of $7.5\times10^{22}$ years at 90% confidence — the first…
desk verdict A clean null result with a real but contained systematic caveat; the limit is plausible, though the stainless steel platform background deserves independent assay before the number is taken at face value. 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 working object is the decay signature itself: $2\nu\beta\beta$-$0^+_1$ events deposit two electron tracks plus gamma rays, either a 760.5 and 818.5 keV cascade or a single 1579.0 keV gamma, so 98.8% of signal events have multiple energy deposits. The analysis classifies events as single-site or multi-site by counting S2 peaks in the bottom photomultiplier waveform, reconstructs both spectra up to MeV energies, and fits them jointly with a binned likelihood. The likelihood treats the $^{232}$Th and $^{238}$U activities of the stainless steel platform as free parameters, while detector-response parameters enter as Gaussian-penalized nuisance parameters, and the overall MS/SS fraction carries a 3.4% and 13.3% systematic uncertainty from calibration data.
What would settle it
Measure the $^{232}$Th and $^{238}$U activities of the stainless steel platform by direct high-purity germanium assay of platform material; if they come out significantly below the fitted 319 and 463 mBq/kg, or if a detailed-geometry simulation changes the best-fit signal by more than the assigned 145-count systematic, then the no-signal conclusion and the $7.5\times10^{22}$ yr limit are not robust.
Extended reading notes
Core claim
On its own terms, the paper establishes that the two-neutrino double $\beta$ decay of $^{136}$Xe to the first excited $0^+_1$ state of $^{136}$Ba has a half-life longer than $7.5\times10^{22}$ years at 90% confidence. The evidence is a simultaneous fit of multi-site and single-site energy spectra in the 1100–2800 keV region, where the signal is expected to appear mostly as multi-site events because its de-excitation gammas scatter. The same fit identifies and, for the first time, evaluates the background from the stainless steel platform outside the cryostat, finding $^{232}$Th and $^{238}$U activities of $319\pm24$ and $463\pm92$ mBq/kg. The best-fit signal count is $220\pm191$, with a p-value of 0.13 against the no-signal hypothesis, so the paper reports no observation. It also validates that the platform background, while negligible for earlier single-site results below 1 MeV, is significant for MeV-scale multi-site analyses and will matter for future neutrinoless double $\beta$ decay searches.
Load-bearing premise
The limit rests on the assumption that the stainless steel platform background has the spectral shape given by a simplified uniform geometry with floating $^{232}$Th and $^{238}$U activities and zero $^{60}$Co and $^{40}$K; if that shape is wrong, the fitted signal of $220\pm191$ counts and hence the half-life limit could shift.
Editorial extensions
If this is right
- A natural-xenon detector without isotopic enrichment can set limits on this decay at the $10^{22}$-year scale, as this first limit demonstrates.
- The stainless steel platform background is negligible for the earlier single-site analyses below 1 MeV but contributes significantly to the MeV-scale multi-site spectrum and to future $^{136}$Xe neutrinoless double beta decay searches.
- Including the platform background in the fit changes the MS $\chi^2/\mathrm{NDF}$ from 938/425 to 543/423, showing that multi-site data strongly constrain external MeV gamma backgrounds.
- The anti-correlation between the fitted signal and the floating platform background means future measurements should reduce or constrain that background, for example with direct assay or the new active water veto.
Reading between the lines
- If direct assay of the platform steel finds $^{232}$Th or $^{238}$U activities far below the fitted 319 and 463 mBq/kg, the $220\pm191$ signal count would shift, and the same data could yield a stronger limit or a signal; the authors leave this re-interpretation implicit.
- The same MeV-scale multi-site reconstruction could be applied to other gamma-emitting rare processes in natural xenon, a generalization the paper does not explore.
- The paper's own fit shows an anti-correlation between signal and platform background, which implies that a detector with an actively vetoed or cleaner support structure could push this search well past $10^{23}$ years; this projection is not worked out in the paper.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a search for two-neutrino double beta decay of 136Xe to the 0^+_1 excited state of 136Ba using 94.9 days of PandaX-4T commissioning data with natural xenon. The analysis reconstructs both multi-site and single-site events, builds a simultaneous binned likelihood fit over the 1100--2800 keV region, and introduces an in-situ background component from the stainless steel platform surrounding the cryostat. No significant excess is found; the best-fit signal is 220 +/- 191 counts, and the paper quotes a 90% confidence level lower limit on the half-life of 7.5 x 10^22 yr, stated as the first such limit from a natural-xenon detector.
Significance. If accepted, this result demonstrates that a large natural-xenon dual-phase TPC can probe 2nu-beta-beta-0^+_1 decays at the 10^22 yr scale, complementing dedicated enriched-xenon searches. The paper's strengths are its use of MS/SS discrimination to extend the energy reach to the MeV scale, the substantial improvement in the background-only fit after adding the stainless steel platform component (chi2/NDF dropping from 938.4/425 to 543.3/423 in the MS spectrum), and the honest reporting of the anti-correlation between the signal and that component. The null conclusion itself is credible; however, the numerical half-life limit rests on a background template whose spectral shape is not independently validated, so the exact limit should be treated with caution until that template is better constrained.
major comments (2)
- [Sec. 4 and Sec. 6.2] The systematic uncertainty assigned to the stainless steel platform (SSP) background does not cover the most dangerous model errors. The relocation test in Sec. 6.2 moves 232Th and 238U sources from 20 to 100 cm from the cryostat while keeping the simplified uniform geometry and fixing 60Co and 40K to zero. It does not vary the source composition, the column/beam geometry, the possibility of surface or point-like contamination, or the spectral shape of the template. Because the paper states in Sec. 6.2 that the fitted signal and the SSP background are anti-correlated, a modest mis-modeling of the SSP spectral shape could shift the best-fit signal by more than the quoted 145 counts. Please add additional cross-checks, such as fits with alternative geometry or composition templates, an explicit scan over allowed 60Co/40K fractions, or a more conservative systematic treatment, before the numerical limit is presented as robust.
- [Sec. 6.2] The conversion from the fitted signal count (220 +/- 191) to the quoted 90% CL half-life limit of 7.5 x 10^22 yr is not documented. The text states that a lower limit is 'achieved' but does not specify how the upper bound on the signal count is derived from the likelihood, how the nuisance parameters and the 145-count systematic are incorporated into the limit, or how the exposure of 12.57 kg x yr and the 136Xe fraction enter the calculation. Please provide the limit-setting procedure, including the test statistic and the treatment of the systematic uncertainty, so that the central claim can be reproduced from the described inputs.
minor comments (5)
- [Sec. 4] There is a typo in 'shilding water tank'; it should read 'shielding water tank'.
- [Eq. (5.1) and Eq. (5.2)] The roles of the nuisance parameters eta_k_eff, eta_s, and eta_k_b are not fully defined in the text. Please specify which systematic in Table 1 each eta corresponds to, and state whether the MS and SS efficiency uncertainties are treated as uncorrelated in the fit.
- [Sec. 6.1, bullet list] The word 'nuissance' should be 'nuisance'.
- [Sec. 6.2 and Fig. 6] The figure labels 'chi2/NDF = 1.28' and 'chi2/NDF = 1.04' are difficult to read at reproduction size; please increase the font size and clarify which fit each value refers to.
- [References] Reference [37] is a bare URL; a formal citation with author, title, and year would be more appropriate for a journal submission.
Circularity Check
No significant circularity: the excited-state signal count is a free fit parameter, and the self-cited 2νββ half-life and screening results are independent measured inputs.
full rationale
The central result is obtained by a binned likelihood fit (Eq. 5.1) in which the 136Xe 2νββ-0+1 signal count n_s is a free parameter; it is not defined in terms of any input, and the signal spectral templates come from Decay0/BambooMC Monte Carlo, not from the data being fit. The self-citations that appear are used as inputs to the background model, not as the predicted quantity: the 136Xe 2νββ ground-state half-life from Ref. [7] normalizes Category 1 background, and the HPGe screening results of Ref. [29] anchor detector-material activities. Both are independently measured quantities with stated assumptions and are externally falsifiable, so their reuse is legitimate evidence rather than circularity. The stainless-steel platform background is floated in the fit, and the paper explicitly reports an anti-correlation with the signal; this is a statistical property of a real spectral decomposition, and the 20-100 cm relocation test provides a systematic uncertainty estimate rather than forcing the signal to equal a background. No uniqueness theorem, ansatz, or defining relation is imported from prior PandaX papers to make the limit follow by construction. The result is also benchmarked against external experimental limits from EXO-200 and KamLAND-Zen, so it is not a renaming of a known result. The remaining concern about the unassayed SSP spectral shape is a model-robustness/correctness issue, not a circularity issue, and the paper itself flags it as needing future direct assay.
Assumptions & free parameters
free parameters (5)
- SSP 232Th activity =
319 ± 24 mBq/kg (averaged)
- SSP 238U activity =
463 ± 92 mBq/kg (averaged)
- Signal event count n_s =
220 ± 191 after SSP systematic
- Detector material component counts =
MS 2826/5919/8901/2671 and SS 343/733/729/263 for 60Co, 40K, 232Th, 238U
- MS/SS energy response parameters a,b,c,d,e =
Eq. 3.2 and Eq. 3.3 with covariance matrix
assumptions (6)
- domain assumption BambooMC/Geant4 and Decay0 provide reliable spectral shapes and SS/MS fractions for signal and backgrounds at MeV energies.
- domain assumption The stainless steel platform background can be modeled as uniform, uniformly distributed 232Th and 238U sources in a simplified geometry, with 60Co and 40K set to zero.
- domain assumption The collaboration's own measured 136Xe 2nu beta beta ground-state half-life [7] is used for the dominant Category 1 background normalization.
- domain assumption SS/MS discrimination via S2B peak counting is valid at MeV energies and its calibration residual is covered by 3.4% (MS) and 13.3% (SS) uncertainties.
- domain assumption The fiducial volume, 136Xe abundance (8.58 ± 0.11)%, and density give a correct exposure of 12.57 kg·yr.
- standard math The binned Poisson likelihood with Gaussian nuisance penalties in Eq. (5.1) yields valid confidence intervals.
Cite this review
Pith. "Pith review of Search for Double Beta Decay of $^{136}$Xe to the $0^+_1$ Excited State of $^{136}$Ba with PandaX-4T." pith.science (2026). https://pith.science/paper/XOKTXE4P
@misc{pith2026250203017,
author = {Pith},
title = {Pith review of: Search for Double Beta Decay of $^136$Xe to the $0^+_1$ Excited State of $^136$Ba with PandaX-4T},
year = {2026},
howpublished = {\url{https://pith.science/paper/XOKTXE4P}},
note = {Machine review of arXiv:2502.03017}
}
abstract
We perform a search of double beta decay of $^{136}$Xe to the excited state, $0^+_1$, of $^{136}$Ba (2$\nu \beta \beta$-0$_1^+$), using the dual-phase xenon detector of PandaX-4T with the first 94.9-day commissioning data. The multi-site events are reconstructed up to the MeV energy scale, which helps to improve the background model significantly. The background contribution from the stainless steel platform outside PandaX-4T cryostat is evaluated for the first time. No significant evidence for 2$\nu \beta \beta$-$0_1^+$ is observed, resulting in a lower limit on half-life of $7.5 \times 10^{22}$ yr at the 90% confidence level. This is the first experimental limit on such a rare decay in a natural xenon-based detector.
Forward citations
Cited by 1 Pith paper
-
Sensitivity of neutrinoless double beta decays from a combined analysis of ground and excited states
Combined analysis of 0νββ decays to ground and excited states can significantly enhance experimental sensitivity depending on NME predictions.
Reference graph
Works this paper leans on
-
[1]
Goeppert-Mayer, Double Beta-disintegration, Phys
M. Goeppert-Mayer, Double Beta-disintegration, Phys. Rev. 48 (1935) 512
work page 1935
-
[2]
Saakyan, Two-Neutrino Double-Beta Decay, Ann
R. Saakyan, Two-Neutrino Double-Beta Decay, Ann. Rev. Nucl. Part. Sci.63 (2013) 503
work page 2013
-
[3]
Furry, On transition probabilities in double beta-disintegration, Phys
W.H. Furry, On transition probabilities in double beta-disintegration, Phys. Rev. 56 (1939) 1184. – 15 –
work page 1939
-
[4]
F.T. A vignone, III, S.R. Elliott and J. Engel, Double Beta Decay, Majorana Neutrinos, and Neutrino Mass, Rev. Mod. Phys.80 (2008) 481 [ 0708.1033]
arXiv 2008
-
[5]
KamLAND-Zen collaboration, Precision measurement of the136Xe two-neutrino ββ spectrum in KamLAND-Zen and its impact on the quenching of nuclear matrix elements, Phys. Rev. Lett.122 (2019) 192501 [ 1901.03871]
work page Pith review arXiv 2019
-
[6]
NEXT collaboration, Measurement of the136Xe two-neutrino double-β-decay half-life via direct background subtraction in NEXT, Phys. Rev. C 105 (2022) 055501 [ 2111.11091]
work page Pith review arXiv 2022
-
[7]
PandaX collaboration, Determination of Double Beta Decay Half-Life of136Xe with the PandaX-4T Natural Xenon Detector, Research 2022 (2022) 9798721 [ 2205.12809]
arXiv 2022
-
[8]
GERDA, (GERDA Collaboration)*collaboration, Final Results of GERDA on the Two-Neutrino Double-β Decay Half-Life of Ge76, Phys. Rev. Lett.131 (2023) 142501 [2308.09795]
work page Pith review arXiv 2023
Show all 38 references
-
[9]
Majorana collaboration, Final Result of the Majorana Demonstrator’s Search for Neutrinoless Double-β Decay in 76Ge, Phys. Rev. Lett.130 (2023) 062501 [ 2207.07638]
2023 arXiv
-
[10]
CUORE collaboration, Search for Majorana neutrinos exploiting millikelvin cryogenics with CUORE, Nature 604 (2022) 53
2022
-
[11]
KamLAND-Zen collaboration, Search for the Majorana Nature of Neutrinos in the Inverted Mass Ordering Region with KamLAND-Zen, Phys. Rev. Lett.130 (2023) 051801 [2203.02139]
2023
-
[12]
GERDA collaboration, Latest results from GERDA Phase II experiment on76Ge double-beta decay and exotic decay searches, PoS EPS-HEP2023 (2024) 158
2024
-
[13]
Jokiniemi, B
L. Jokiniemi, B. Romeo, C. Brase, J. Kotila, P. Soriano, A. Schwenk et al., Two-neutrino ββ decay of 136Xe to the first excited 0+ state in 136Ba, Phys. Lett. B838 (2023) 137689 [2211.03764]
2023 arXiv
-
[14]
Barabash, Double beta decay to the excited states: Review, AIP Conf
A.S. Barabash, Double beta decay to the excited states: Review, AIP Conf. Proc.1894 (2017) 020002 [ 1709.06890]
2017 arXiv
-
[15]
Barea, J
J. Barea, J. Kotila and F. Iachello, 0νββ and 2νββ nuclear matrix elements in the interacting boson model with isospin restoration, Phys. Rev. C 91 (2015) 034304 [ 1506.08530]
2015 arXiv
-
[16]
KamLAND-Zen collaboration, Search for double-beta decay of136Xe to excited states of 136Ba with the KamLAND-Zen experiment, Nucl. Phys. A946 (2016) 171 [ 1509.03724]
2016 arXiv
-
[17]
Barabash et al., Two neutrino double beta decay of Mo-100 to the first excited 0+ state in Ru-100, Phys
A.S. Barabash et al., Two neutrino double beta decay of Mo-100 to the first excited 0+ state in Ru-100, Phys. Lett. B345 (1995) 408
1995
-
[18]
NEMO-3 collaboration, Investigation of double beta decay of100Mo to excited states of 100Ru, Nucl. Phys. A925 (2014) 25 [ 1402.7196]
2014 arXiv
-
[19]
Kidd, J.H
M.F. Kidd, J.H. Esterline, S.W. Finch and W. Tornow, Two-neutrino double-β decay of 150Nd to excited final states in150Sm, Phys. Rev. C 90 (2014) 055501 [ 1411.3755]
2014 arXiv
-
[20]
CUORE collaboration, Search for double-beta decay of130Te to the 0+ states of 130Xe with CUORE, Eur. Phys. J. C81 (2021) 567 [ 2101.10702]
2021 arXiv
-
[21]
MAJORANA collaboration, Search for double-β decay of 76Ge to excited states of76Se with the MAJORANA DEMONSTRATOR, Phys. Rev. C 103 (2021) 015501 [ 2008.06014]. – 16 –
2021 arXiv
-
[22]
EXO-200 collaboration, Search for two-neutrino double-beta decay of136Xe to the excited state of 136Ba with the complete EXO-200 dataset*, Chin. Phys. C 47 (2023) 103001 [2303.01103]
2023 arXiv
-
[23]
PandaX-4T collaboration, Dark Matter Search Results from the PandaX-4T Commissioning Run, Phys. Rev. Lett.127 (2021) 261802 [ 2107.13438]
2021 arXiv
-
[24]
PandaX collaboration, Searching for Two-Neutrino and Neutrinoless Double Beta Decay of Xe134 with the PandaX-4T Experiment, Phys. Rev. Lett.132 (2024) 152502 [ 2312.15632]
2024 arXiv
-
[25]
PandaX collaboration, Search for Cosmic-Ray Boosted Sub-MeV Dark-Matter–Electron Scattering in PandaX-4T, Phys. Rev. Lett.133 (2024) 101805 [ 2403.08361]
2024 arXiv
-
[26]
(PandaX)„ PandaXcollaboration, Measurement of solar pp neutrino flux using electron recoil data from PandaX-4T commissioning run*, Chin. Phys. C 48 (2024) 091001 [2401.07045]
2024 arXiv
-
[27]
PandaX collaboration, Search for Dark-Matter–Nucleon Interactions with a Dark Mediator in PandaX-4T, Phys. Rev. Lett.131 (2023) 191002 [ 2308.01540]
2023 arXiv
-
[28]
J. Li, X. Ji, W. Haxton and J.S.Y. Wang, The second-phase development of the China JinPing underground Laboratory, Phys. Procedia61 (2015) 576 [ 1404.2651]
2015 arXiv
-
[29]
PandaX-4T collaboration, Low radioactive material screening and background control for the PandaX-4T experiment, JHEP 06 (2022) 147 [ 2112.02892]
2022 arXiv
-
[30]
PandaX collaboration, Searching for MeV-scale Axion-like Particles and Dark Photons with PandaX-4T, 2409.00773
-
[31]
PandaX collaboration, Measurement of two-neutrino double electron capture half-life of 124Xe with PandaX-4T, 2411.14355
-
[32]
GEANT4 collaboration, GEANT4–a simulation toolkit, Nucl. Instrum. Meth. A506 (2003) 250
2003
-
[33]
Chen et al., BambooMC — A Geant4-based simulation program for the PandaX experiments, JINST 16 (2021) T09004 [ 2107.05935]
X. Chen et al., BambooMC — A Geant4-based simulation program for the PandaX experiments, JINST 16 (2021) T09004 [ 2107.05935]
2021 arXiv
-
[34]
Ponkratenko, V.I
O.A. Ponkratenko, V.I. Tretyak and Y.G. Zdesenko, The Event generator DECAY4 for simulation of double beta processes and decay of radioactive nuclei, Phys. Atom. Nucl.63 (2000) 1282 [ nucl-ex/0104018]
2000 arXiv
-
[35]
PANDA-X collaboration, Horizontal position reconstruction in PandaX-II, JINST 16 (2021) P11040 [ 2106.08380]
2021 arXiv
-
[36]
Luo et al., Improvement on the linearity response of PandaX-4T with new photomultiplier tube bases, JINST 19 (2024) P05021 [ 2401.00373]
L. Luo et al., Improvement on the linearity response of PandaX-4T with new photomultiplier tube bases, JINST 19 (2024) P05021 [ 2401.00373]
2024 arXiv
-
[37]
https://www.hidenanalytical.com/products/residual-gas-analysis/
-
[38]
PandaX collaboration, Searching for Neutrinoless Double-Beta Decay of136Xe with PandaX-4T, 2412.13979. – 17 –
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.