REVIEW 2 major objections 3 minor 1 cited by
Direct Measurement of the $^{39}$Ar Half-life from 3.4 Years of Data with the DEAP-3600 Detector
T0 review · 2 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read DEAP-3600 directly observed the $^{39}$Ar decay curve over 3.4 years and measured the half-life as $(302 \pm 8_{\mathrm{stat}} \pm 6_{\mathrm{sys}})$ years, about 12% longer than the accepted 269-year value.
desk verdict First direct 39Ar half-life measurement from a decay curve, but the light-yield correction as written is inverted and could shift the central value past the quoted systematic. 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 machine that carries the argument is the detector's ability to count $^{39}$Ar decays continuously and to model the trigger rate. DEAP-3600's liquid argon target is so rich in $^{39}$Ar that this isotope dominates triggers by two orders of magnitude; pulse-shape discrimination (the ratio of prompt to total charge, $F_{\mathrm{prompt}}$) separates the electron-recoil band from nuclear recoils, and a data-cleaning cut plus livetime calculation converts raw triggers into rates. The rate model explicitly includes Poisson probabilities for uncorrelated pile-up (one, two, or three $^{39}$Ar decays landing in the same 10 $\mu$s trigger window, and $^{39}$Ar plus Cherenkov light), so the exponential decay $R_{^{39}\mathrm{Ar}}(t)=R_{^{39}\mathrm{Ar}}e^{-t/\tau_{^{39}\mathrm{Ar}}}$ can be extracted with the mean lifetime $\tau$ as a free parameter. A light-yield correction, derived from daily calibrations with $^{40}$K, $^{208}$Tl, and the $^{39}$Ar $\beta$ spectrum, is applied to each event's PE value to remove time-dependent detector response changes.
What would settle it
Re-fit the trigger-rate curve with the light-yield correction omitted and again with its sign reversed (dividing by $Y_j/Y_0$ instead of multiplying). If either variant shifts the half-life by more than the quoted $\pm 6$-year systematic, the central result is controlled by the correction rather than by the decay; a second, independent check is to fit the two halves of the 3.4-year dataset separately and compare their $T_{1/2}$ values.
Extended reading notes
Core claim
The central discovery is a half-life value obtained from a direct decay-rate fit. Selecting clean electron-recoil events with pulse-shape discrimination in the 700–1200 PE window, the collaboration measured the $^{39}$Ar trigger rate as a function of time and fit it with a model that includes single decays, double and triple pile-up, $^{39}$Ar–Cherenkov pile-up, and a constant ERB $\gamma$ background. The fit yields $T_{1/2}=(302 \pm 8_{\mathrm{stat}} \pm 6_{\mathrm{sys}})$ years, and the paper reports a $p$-value of 0.008 for consistency with the NDS value, i.e. the two are in tension. When combined with the earlier DEAP specific-activity measurement, this corresponds to a $^{39}$Ar abundance in atmospheric argon of $(8.6 \pm 0.4)\times 10^{-16}$.
Load-bearing premise
The load-bearing assumption is that the PE-scale correction of Section 3.1 completely removes time-dependent light-yield variations, so that the remaining change in trigger rate is purely the $^{39}$Ar decay; if the correction leaves, or introduces, a drift, the measured half-life changes by more than the quoted systematic uncertainty.
Editorial extensions
If this is right
- If the 302-year value is right, the 269-year standard is about 12% too short, and the NDS evaluation (268 ± 8 years) disagrees at a p-value of 0.008.
- Combining the new half-life with DEAP's measured specific activity gives an atmospheric $^{39}$Ar abundance of $(8.6 \pm 0.4)\times10^{-16}$, a direct estimate that does not depend on the older half-life.
- Radiometric dating and tracer methods that use $^{39}$Ar would need to update their decay constants, shifting inferred ages by the corresponding factor.
- Future dark-matter experiments using atmospheric argon would see a slightly higher $^{39}$Ar background rate per unit mass if the half-life is longer, which affects their background-model inputs.
Reading between the lines
- Not stated in the paper, the sign of the light-yield correction in Section 3.1 is worth an independent check: the text says each event's PE value is corrected by multiplying by $Y_j/Y_0$, and if the intended correction was to divide by that ratio, the drift used to extract the half-life would move in the opposite direction and could shift the central value by more than the quoted $\pm 6$-year syst
- Because the light-yield correction dominates the systematic budget, an external test would be to fit the 700–900 PE and 900–1200 PE sub-regions separately; a significant disagreement between the two inferred half-lives would indicate residual spectral-shape sensitivity.
- The method of continuously counting a dominant endogenous decay and modeling pile-up statistics could be transplanted to other long-lived isotopes in large liquid detectors, such as $^{85}$Kr in liquid argon or $^{14}$C in organic scintillators, where no direct decay-curve half-life measurement exists.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports the first direct measurement of the 39Ar half-life by observing the decay curve with the DEAP-3600 liquid-argon detector over 3.4 years. The authors select electron-recoil events in a 700-1200 PE region of interest, model the ROI trigger rate as a combination of single, double, and triple 39Ar pile-up, 39Ar-Cherenkov pile-up, and a constant gamma background, and fit for the initial 39Ar rate and lifetime. The result is T1/2 = (302 ± 8stat ± 6sys) years, about 12% longer than the NDS value, with a reported p-value of 0.008 for disagreement. The systematic budget is dominated by light yield corrections.
Significance. If the result is correct, it is a valuable first direct decay-curve measurement of the 39Ar half-life, with implications for geochronology, radiochemistry, and argon-based dark matter detectors. The paper's strengths include a clean derivation of the fit model from Poisson statistics, a free-floating half-life (so no circularity in the central fit), very high statistics, explicit Monte Carlo and data cross-checks, and a reasonably complete systematic inventory. The main weakness is the light yield correction described in Section 3.1: as written, the direction of the correction appears inverted, and the dominant systematic is evaluated by perturbing that same factor. Because the expected decay signal over 3.4 years is only about 0.8% in rate, an error in this correction could change the central value well beyond the quoted uncertainties. The result is therefore conditional on resolving this point.
major comments (2)
- [Sec. 3.1 and Sec. 6.1] The stated correction factor appears inverted. Section 3.1 says 'The PE value of each event in a run is then corrected by a factor of Yj/Y0,' where Y is the light yield in PE/keVee. For a standard light yield definition, an event with a given energy has PE_j = Y_j * E, so the mapping to the reference run is PE_0 = (Y_0 / Y_j) * PE_j. The paper's formula multiplies by Y_j / Y_0, which would amplify, not remove, run-to-run light yield drift; Section 6.1 repeats this form as 'Yj/Y0 + ΔY'. Since the ROI is fixed in PE, an inverted correction changes the fraction of the 39Ar spectrum inside the ROI as a function of time, potentially mimicking a change in the decay rate. The quoted differential systematic (5.1 years) is obtained by perturbing this factor in magnitude, not by checking its sign, so it does not cover an inversion. The authors must clarify whether a non-standard definition of Y is being used, correct the formula if it is a typo, and show with a closure test (e.g., Monte Carlo injection of a known lifetime or a residual-rate versus light-yield check) that the implemented correction is in the direction that removes drift.
- [Sec. 6.1 and Fig. 6] The differential drift systematic is evaluated by adding a linear ramp of amplitude ΔY = 0.00076 to the correction factor and refitting. However, the calibration source light yields plotted in Fig. 6 show substantially larger, non-linear time variations (of order 0.1 PE/keVee over the dataset). While the nominal run-by-run correction may remove these variations if the calibration measurements perfectly track the detector response relevant to the ROI, the paper does not provide a direct closure check demonstrating that the corrected rates are insensitive to the observed non-linear light yield structure. Given that the physics signal is only about 0.8% in rate over the full dataset, a non-linear residual drift at the level visible in Fig. 6 could be comparable to the signal. The revision should include an explicit validation that the applied correction removes the observed time-dependent light yield changes, beyond the linear-ramp perturbation used for the systematic.
minor comments (3)
- [Title/header] The manuscript header contains '3.4 Y ears of Data'; the word 'Years' should be restored.
- [Sec. 3] The text 'Individual runtimes Trun vary in the range of [ O(1minute ), ∼2 days]' has spacing and notation issues; it should read 'O(1 minute)' or a conventional notation for order-of-magnitude ranges.
- [Table 2] In the 'Constraints' column, entries such as 'Run-dependent N/A' are unclear: the table should specify that the light yield correction factors are run-dependent parameters whose uncertainties are evaluated as described in Section 6.1, and what 'N/A' refers to.
Circularity Check
No significant circularity: the 39Ar half-life is a free fit parameter and no load-bearing step reduces to an input.
full rationale
The central claim is derived by fitting the observed time-dependent trigger rate with Eq. 13, where R39Ar(t) = R39Ar exp(-t/tau39Ar) and tau39Ar is explicitly a free-floating parameter. The auxiliary inputs are determined independently of the target half-life: RChv is measured from low-threshold data and matched by MC (Sec. 5.1); Rbg is obtained by subtracting the MC 39Ar spectrum from the total MC ERB sum (Sec. 5.3); the ROI fractions fi come from a toy MC based on the external beta-spectrum model of Ref. [28] and a detector response model; the efficiencies eps_i come from MC and data-side selection studies without reference to the decay rate change. Thus the half-life is not an input to the fit, and no equation in the paper defines the result in terms of itself. The self-citations (prior DEAP specific-activity measurement [7] and PSD analysis [25]) are not load-bearing for the half-life: [7] is used only in the conclusion to convert the measured half-life into an abundance, and [25] supports the established pulse-shape discrimination technique. There is no imported uniqueness theorem, no ansatz smuggled in via citation, and no renaming of a known empirical pattern. The reader-identified concern about the direction of the light-yield correction factor (Sec. 3.1) is a potential implementation/correctness issue, not a circularity: even if the factor were inverted, that would bias the fitted rate trend, but it would not make the half-life equivalent to an input by construction. The paper also honestly reports the light-yield correction as the dominant systematic (5.1 years), so no fitted parameter is being relabeled as a prediction. Overall, the derivation chain is self-contained against the measured decay curve, and no step reduces to its own inputs.
Assumptions & free parameters
free parameters (9)
- Rbg (ERB gamma background rate) =
1.65 ± 0.31 Hz
- RChv (Cherenkov trigger rate) =
538 ± 4 Hz
- f1 (single 39Ar spectrum fraction in ROI) =
0.21 ± 0.05
- f2 (double 39Ar pile-up spectrum fraction in ROI) =
0.20 ± 0.05
- f3 (triple 39Ar pile-up spectrum fraction in ROI) =
0.19 ± 0.05
- f4 (39Ar-Cherenkov pile-up spectrum fraction in ROI) =
0.21 ± 0.05
- ε2 (double 39Ar pile-up selection efficiency) =
0.9099 ± 0.0033
- ε3 (triple 39Ar pile-up selection efficiency) =
0.860 ± 0.039
- Light yield ratio Yj/Y0 =
Run-dependent, average ΔY = 0.00076
assumptions (6)
- standard math Poisson statistics for the occurrence of 39Ar decays in the 10 µs trigger window (Eqs. 5-12)
- domain assumption The 39Ar beta spectrum shape from Kostensalo, Suhonen, and Zuber (Ref. [28])
- domain assumption Detector response is a Gaussian smearing with mean <NDN> + Y*E (Eq. 14)
- domain assumption The ERB gamma background rate Rbg is constant across the dataset
- domain assumption The Cherenkov trigger rate RChv is constant in time
- domain assumption The data-cleaning cut and livetime calculation in Eq. 2 correctly account for deadtime
Cite this review
Pith. "Pith review of Direct Measurement of the $^{39}$Ar Half-life from 3.4 Years of Data with the DEAP-3600 Detector." pith.science (2026). https://pith.science/paper/HUYFEZYV
@misc{pith2026250113196,
author = {Pith},
title = {Pith review of: Direct Measurement of the $^39$Ar Half-life from 3.4 Years of Data with the DEAP-3600 Detector},
year = {2026},
howpublished = {\url{https://pith.science/paper/HUYFEZYV}},
note = {Machine review of arXiv:2501.13196}
}
abstract
The half-life of $^{39}$Ar is measured using the DEAP-3600 detector located 2 km underground at SNOLAB. Between 2016 and 2020, DEAP-3600 used a target mass of (3269 $\pm$ 24) kg of liquid argon distilled from the atmosphere in a direct-detection dark matter search. Such an argon mass also enables direct measurements of argon isotope properties. The decay of $^{39}$Ar in DEAP-3600 is the dominant source of triggers by two orders of magnitude, ensuring high statistics and making DEAP-3600 well-suited for measuring this isotope's half-life. Use of the pulse-shape discrimination technique in DEAP-3600 allows powerful discrimination between nuclear recoils and electron recoils, resulting in the selection of a clean sample of $^{39}$Ar decays. Observing over a period of 3.4 years, the $^{39}$Ar half-life is measured to be $(302 \pm 8_{\rm stat} \pm 6_{\rm sys})$ years. This new direct measurement suggests that the half-life of $^{39}$Ar is significantly longer than the accepted value, with potential implications for measurements using this isotope's half-life as input.
Forward citations
Cited by 1 Pith paper
-
Community Report from the 2025 SNOLAB Future Projects Workshop
A community report summarizes proposed future experiments and infrastructure needs for SNOLAB over the next 15 years.
Reference graph
Works this paper leans on
-
[1]
Böhlke, PAC 86(9), 1421 (2014)
J. Böhlke, PAC 86(9), 1421 (2014). DOI doi:10.1515/ pac-2013-0918. URL https://doi.org/10.1515/ pac-2013-0918
work page 2014
-
[2]
P. Benetti, et al. (W ARP Collaboration), NIM-A574(1), 83 (2007). DOI 10.1016/j.nima.2007.01.106. URL https://doi.org/10.1016/j.nima.2007.01.106
-
[3]
J. Calvo, et al. (ArDM Collaboration), JCAP 2018(12), 011 (2018). DOI 10.1088/1475-7516/2018/12/
-
[4]
Saldanha et al., PRC 100(2) (2019)
R. Saldanha et al., PRC 100(2) (2019). DOI 10.1103/ physrevc.100.024608. URL https://doi.org/10. 1103/physrevc.100.024608
work page 2019
-
[5]
H. Loosli, H. Oeschger, EPSL 5, 191 (1968). DOI 10.1016/S0012-821X(68)80039-1. URL https://doi. org/10.1016/S0012-821X(68)80039-1
-
[6]
A.R. Brosi, H. Zeldes, B.H. Ketelle, Phys. Rev. 79, 902 (1950). DOI 10.1103/PhysRev.79.902. URL https: //doi.org/10.1103/PhysRev.79.902
-
[7]
P. Adhikari, et al. (DEAP Collaboration), EPJ C 83(642) (2023). DOI 10.1140/epjc/s10052-023-11678-6. URL https://doi.org/10.1140/epjc/ s10052-023-11678-6
-
[8]
A.K. Baksi, Earth Syst. Sci. 131(157) (2022). DOI 10. 1007/s12040-022-01855-8. URL https://doi.org/ 10.1007/s12040-022-01855-8
Show all 32 references
-
[9]
C. Yu, R.C. Zhong, Y .L. Xie, W.B. Li, Minerals9(10) (2019). DOI 10.3390/min9100622. URL https://doi. org/10.3390/min9100622
2019 doi
-
[10]
Ritterbusch, et al., PNAS 119(40), e2200835119 (2022)
F. Ritterbusch, et al., PNAS 119(40), e2200835119 (2022). DOI 10.1073/pnas.2200835119. URL https: //doi.org/10.1073/pnas.2200835119
2022 doi
-
[11]
URL https://doi.org/10.1088/1475-7516/ 2018/12/011
2018 doi
-
[12]
Pool, J.M
M.L. Pool, J.M. Cork, R.L. Thornton, Phys. Rev.52, 239 (1937). DOI 10.1103/PhysRev.52.239. URL https: //doi.org/10.1103/PhysRev.52.239
1937 doi
-
[13]
Seaborg, I
G.T. Seaborg, I. Perlman, RMP 20, 585 (1948). DOI 10. 1103/RevModPhys.20.585. URL https://doi.org/ 10.1103/RevModPhys.20.585
1948 doi
-
[14]
Zeldes, et al., Phys
H. Zeldes, et al., Phys. Rev. 86, 811 (1952). DOI 10.1103/PhysRev.86.811. URL https://doi.org/10. 1103/PhysRev.86.811
1952 doi
-
[15]
Stoenner, O.A
R.W. Stoenner, O.A. Schaeffer, S. Katcoff, Science 148(3675), 1325 (1965). DOI 10.1126/science.148.3675
1965 doi
-
[16]
Baksi, D.A
A.K. Baksi, D.A. Archibald, E. Farrar, Canadian Journal of Physics 74(5-6), 263 (1996). DOI 10.1139/p96-041. URL https://doi.org/10.1139/p96-041
1996 doi
-
[17]
Stoenner, O.A
R.W. Stoenner, O.A. Schaeffer, R. Davis Jr., JGR65(8), 2521 (1960). DOI 10.1029/JZ065i010p03025. URL https://doi.org/10.1029/JZ065i010p03025
1960 doi
-
[18]
Amaudruz, et al
P.-A. Amaudruz, et al. (DEAP Collaboration), Astropar- ticle Physics 108, 1 (2019). DOI 10.1016/j.astropartphys. 2018.09.006. URL https://doi.org/10.1016/j. astropartphys.2018.09.006
2019 doi
-
[19]
Chen, Nuclear Data Sheets 149 (2018)
J. Chen, Nuclear Data Sheets 149 (2018). DOI 10.1016/j. nds.2018.03.001. URL https://doi.org/10.1016/ j.nds.2018.03.001
2018 doi
-
[20]
(DEAP Collaboration), NIM-A 922, 373 (2019)
P.-A Amaudruz, et al. (DEAP Collaboration), NIM-A 922, 373 (2019). DOI 10.1016/j.nima.2018.12.058. URL http://doi.org/10.1016/j.nima.2018.12.058
2019 doi
-
[21]
Francini, et al., JINST 8(09), C09010 (2013)
R. Francini, et al., JINST 8(09), C09010 (2013). DOI 10.1088/1748-0221/8/09/C09010. URL https://doi. org/10.1088/1748-0221/8/09/C09010
2013 doi
-
[22]
Seibert, et al
S. Seibert, et al. RAT (is an Analysis Tool) User’s Guide (2014). URL https://rat.readthedocs.io/ en/latest/
2014
-
[23]
Phys.: Conf
Thomas Lindner for the DEAP-3600 collaboration, J. Phys.: Conf. Ser. 664(8), 082026 (2015). DOI 10.1088/ 1742-6596/664/8/082026. URL https://doi.org/ 10.1088/1742-6596/664/8/082026
2015 doi
-
[24]
R. Brun, F. Rademakers, NIM-A389(1), 81 (1997). DOI 10.1016/S0168-9002(97)00048-X. URL https://doi. org/10.1016/S0168-9002(97)00048-X
1997 doi
-
[25]
Agostinelli, et al., NIM-A 506(3), 250 (2003)
S. Agostinelli, et al., NIM-A 506(3), 250 (2003). DOI 10.1016/S0168-9002(03)01368-8. URL https://doi. org/10.1016/S0168-9002(03)01368-8
2003 doi
-
[26]
Heindl, et al., EPL 91(6), 62002 (2010)
T. Heindl, et al., EPL 91(6), 62002 (2010). DOI 10. 1209/0295-5075/91/62002. URL https://doi.org/ 10.1209/0295-5075/91/62002
2010 doi
-
[27]
Adhikari, et al
P. Adhikari, et al. (DEAP Collaboration), EPJ C 81(823) (2021). DOI 10.1140/epjc/s10052-021-09514-w. URL https://doi.org/10.1140/epjc/ s10052-021-09514-w
2021 doi
-
[28]
Kostensalo, J
J. Kostensalo, J. Suhonen, K. Zuber, Journal of Physics G: Nuclear and Particle Physics 45(2), 025202 (2018). DOI 10.1088/1361-6471/aa958e. URL https://doi. org/10.1088/1361-6471/aa958e
2018 doi
-
[29]
Adhikari, et al
P. Adhikari, et al. (DEAP Collaboration), EPJ C 80(303) (2020). DOI 10.1140/epjc/s10052-020-7789-x. URL https://doi.org/10.1140/epjc/ s10052-020-7789-x 11
2020 doi
-
[30]
Ajaj, et al
R. Ajaj, et al. (DEAP Collaboration), PRD100(7) (2019). DOI 10.1103/physrevd.100.072009. URL http://doi. org/10.1103/PhysRevD.100.072009
2019 doi
-
[31]
Ajaj, et al
R. Ajaj, et al. (DEAP Collaboration), PRD 100, 022004 (2019). DOI 10.1103/PhysRevD.100.022004. URL https://doi.org/10.1103/PhysRevD.100. 022004
2019 doi
-
[1325]
148.3675.1325
URL https://doi.org/10.1126/science. 148.3675.1325
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.