Pith. sign in

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 →

arxiv 2502.03017 v3 pith:XOKTXE4P submitted 2025-02-05 nucl-ex

classification nucl-ex
keywords doublebetadecay136XeexcitedstatePandaX-4Thalf-lifelimitmulti-siteeventsstainlesssteelplatformbackgroundnaturalxenondetector
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper searches for a rare nuclear decay: $^{136}$Xe emits two electrons and two antineutrinos and lands in the first excited $0^+_1$ state of $^{136}$Ba, which then releases gamma rays. Such transitions are standard-model allowed but much rarer than ground-state double $\beta$ decay, and their rates constrain the nuclear matrix elements needed for neutrinoless double $\beta$ decay searches. Using the first 94.9 days of PandaX-4T data, the analysis reconstructs multi-site events at MeV energies and fits the single-site and multi-site spectra together. No signal appears, and the paper quotes a lower limit on the half-life of $7.5\times10^{22}$ years at 90% confidence. This is the first such limit from a natural-xenon detector.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

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)
  1. [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.
  2. [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)
  1. [Sec. 4] There is a typo in 'shilding water tank'; it should read 'shielding water tank'.
  2. [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.
  3. [Sec. 6.1, bullet list] The word 'nuissance' should be 'nuisance'.
  4. [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.
  5. [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

0 steps flagged · score 0.0 of 10

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 5 free parameters · 6 assumptions · 0 invented entities

The central result rests on Monte Carlo spectral shapes, a uniform-geometry stainless steel platform model, and the collaboration's own previous 2nu beta beta half-life measurement. No new particles or entities are invented. Five groups of fitted parameters are used, four of them nuisance or background normalizations and one being the signal count of interest.

free parameters (5)
  • SSP 232Th activity = 319 ± 24 mBq/kg (averaged)
    Floated in the simultaneous MS/SS spectral fit; this is the dominant new stainless steel platform background component.
  • SSP 238U activity = 463 ± 92 mBq/kg (averaged)
    Floated in the same fit; the SSP contribution is partly degenerate with detector-material backgrounds and with the signal.
  • Signal event count n_s = 220 ± 191 after SSP systematic
    Parameter of interest in the likelihood; no evidence for nonzero signal with p=0.13 and converted into the 7.5e22 yr half-life limit.
  • Detector material component counts = MS 2826/5919/8901/2671 and SS 343/733/729/263 for 60Co, 40K, 232Th, 238U
    Gaussian-constrained nuisance parameters in Eq. (5.1); post-fit values set the subtraction of Category 2 backgrounds.
  • MS/SS energy response parameters a,b,c,d,e = Eq. 3.2 and Eq. 3.3 with covariance matrix
    Five parameters per spectrum, pre-fit in the outer detector region and used as priors; they shape all templates in the FV fit.
assumptions (6)
  • domain assumption BambooMC/Geant4 and Decay0 provide reliable spectral shapes and SS/MS fractions for signal and backgrounds at MeV energies.
    Used in Sec. 3.1 and Sec. 4 to derive a signal MS fraction of 98.8% and all template shapes; no data-driven signal template is provided.
  • 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.
    Introduced in Sec. 4 and Sec. 6.1; the only robustness test moves sources 20 to 100 cm and assigns a 145-count systematic, without direct assay or full geometry variation.
  • 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.
    Sec. 4 Category 1; self-cited from a prior PandaX paper, not derived in this work, but an external input to this likelihood.
  • 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.
    Sec. 3.1; the discrimination is validated only against 232Th calibration data and BambooMC, not by a first-principles waveform model.
  • domain assumption The fiducial volume, 136Xe abundance (8.58 ± 0.11)%, and density give a correct exposure of 12.57 kg·yr.
    Sec. 3.2; an error in any of these scales the half-life limit linearly.
  • standard math The binned Poisson likelihood with Gaussian nuisance penalties in Eq. (5.1) yields valid confidence intervals.
    Standard statistical construction; no formal proof or coverage study is given.

how reviews work

0 comments
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.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Sensitivity of neutrinoless double beta decays from a combined analysis of ground and excited states

    hep-ph 2025-08 unverdicted novelty 5.0 of 10

    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

38 extracted references · 17 canonical work pages · cited by 1 Pith paper

  1. [1]

    Goeppert-Mayer, Double Beta-disintegration, Phys

    M. Goeppert-Mayer, Double Beta-disintegration, Phys. Rev. 48 (1935) 512

  2. [2]

    Saakyan, Two-Neutrino Double-Beta Decay, Ann

    R. Saakyan, Two-Neutrino Double-Beta Decay, Ann. Rev. Nucl. Part. Sci.63 (2013) 503

  3. [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 –

  4. [4]

    A vignone, III, S.R

    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]

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

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

  7. [7]

    PandaX collaboration, Determination of Double Beta Decay Half-Life of136Xe with the PandaX-4T Natural Xenon Detector, Research 2022 (2022) 9798721 [ 2205.12809]

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

Show all 38 references
  1. [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]

  2. [10]

    CUORE collaboration, Search for Majorana neutrinos exploiting millikelvin cryogenics with CUORE, Nature 604 (2022) 53

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

  4. [12]

    GERDA collaboration, Latest results from GERDA Phase II experiment on76Ge double-beta decay and exotic decay searches, PoS EPS-HEP2023 (2024) 158

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

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

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

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

  9. [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

  10. [18]

    NEMO-3 collaboration, Investigation of double beta decay of100Mo to excited states of 100Ru, Nucl. Phys. A925 (2014) 25 [ 1402.7196]

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

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

  13. [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 –

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

  15. [23]

    PandaX-4T collaboration, Dark Matter Search Results from the PandaX-4T Commissioning Run, Phys. Rev. Lett.127 (2021) 261802 [ 2107.13438]

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

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

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

  19. [27]

    PandaX collaboration, Search for Dark-Matter–Nucleon Interactions with a Dark Mediator in PandaX-4T, Phys. Rev. Lett.131 (2023) 191002 [ 2308.01540]

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

  21. [29]

    PandaX-4T collaboration, Low radioactive material screening and background control for the PandaX-4T experiment, JHEP 06 (2022) 147 [ 2112.02892]

  22. [30]

    PandaX collaboration, Searching for MeV-scale Axion-like Particles and Dark Photons with PandaX-4T, 2409.00773

  23. [31]

    PandaX collaboration, Measurement of two-neutrino double electron capture half-life of 124Xe with PandaX-4T, 2411.14355

  24. [32]

    GEANT4 collaboration, GEANT4–a simulation toolkit, Nucl. Instrum. Meth. A506 (2003) 250

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

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

  27. [35]

    PANDA-X collaboration, Horizontal position reconstruction in PandaX-II, JINST 16 (2021) P11040 [ 2106.08380]

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

  29. [37]

    https://www.hidenanalytical.com/products/residual-gas-analysis/

  30. [38]

    PandaX collaboration, Searching for Neutrinoless Double-Beta Decay of136Xe with PandaX-4T, 2412.13979. – 17 –

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.