REVIEW 3 major objections 6 minor 1 cited by
Absorption of Fermionic Dark Matter in the PICO-60 C$_{3}$F$_{8}$ Bubble Chamber
T0 review · 3 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read The PICO-60 bubble chamber, reanalyzed with a non-relativistic effective field theory for absorption, sets the leading spin-independent absorption limits below 23 MeV/c² and the first spin-dependent absorption limits.
desk verdict Real experimental limits on a plausible absorption channel, but the NREFT mapping is asserted rather than derived; the paper should be refereed, with the derivation or a clear caveat required before the spin-dependent limits are 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 carrying machinery is the NREFT operator basis for absorption, with the leading operators $\mathcal{O}_1 = 1$ (spin-independent) and $\mathcal{O}_4 = \vec S_\chi \cdot \vec S_N$ (spin-dependent). The key transformation is the operator substitution $\vec v \rightarrow \vec v_N^{\,\mathrm{avg}}$: replacing the relative velocity between dark matter and nucleus with the average nucleon velocity removes the dark-matter velocity from the interaction, altering the nuclear response functions and producing a peaked recoil signature at $E_R \approx m_\chi^2/(2m_N)$. On the statistical side, the limits come from a profile-likelihood-ratio test in which neutron-calibration efficiency curves are encoded into likelihood surfaces, convolved with thermodynamic-threshold uncertainties, and used to exclude cross sections at 90% confidence; the recoil-energy window runs from the thermodynamic bubble-nucleation threshold to 100 keV, the upper bound chosen because acoustic calibration above about 100 keV is lacking.
What would settle it
A first-principles nuclear-structure calculation of the $\mathcal{O}_4$ spin-dependent response in $^{19}\mathrm{F}$, with the dark-matter velocity contribution removed, would settle the mapping; if the calculated response differs from the code's, the reported spin-dependent absorption limits would not correspond to the claimed process. On the experimental side, a threshold scan that sees a mono-energetic recoil peak at $E_R \approx m_\chi^2/(2m_N)$ followed by an abrupt cutoff once the threshold passes that energy would confirm the absorption signal.
Extended reading notes
Core claim
The paper's central claim is that the non-relativistic effective field theory (NREFT) used for elastic dark-matter–nucleus scattering can be extended to the neutral-current absorption process $\chi + N \rightarrow \nu + N$. The extension consists of replacing the relative-velocity operator in the elastic operators with the average nucleon velocity operator, since the incoming dark matter's velocity is negligible for cold dark matter; this makes the dark-matter side of the interaction velocity-independent and changes the nuclear responses. The resulting recoil spectrum is peaked at $E_R \approx m_\chi^2/(2m_N)$, a signature that cuts off abruptly once the detector threshold exceeds that energy. Applying this prescription to PICO-60 C$_3$F$_8$ data with a profile-likelihood analysis, the paper reports 90% confidence level upper limits: the leading spin-independent absorption limits below 23 MeV/$c^2$ for the operator $\mathcal{O}_1 = 1$, and the first spin-dependent absorption limits for the operator $\mathcal{O}_4 = \vec S_\chi \cdot \vec S_N$, with the strongest $\mathcal{O}_4$ limits in the 12 to 27 MeV/$c^2$ range. The limits are cast as generic cross sections $\sigma_{NC} = m_\chi^2/(4\pi\Lambda^4)$ to permit direct comparison with existing literature.
Load-bearing premise
The load-bearing premise is that the NREFT operators for the absorption process $\chi N \rightarrow \nu N$ are correctly obtained from elastic scattering by replacing the relative velocity operator with the average nucleon velocity operator, and that the nuclear responses change in the way the modified code implements; the paper gives no derivation or independent benchmark for those changed responses, especially for $\mathcal{O}_4$.
Editorial extensions
If this is right
- Bubble chambers become competitive probes of MeV-scale fermionic dark-matter absorption, because the light C$_3$F$_8$ nucleus pushes the characteristic recoil energy higher than heavier targets would.
- The $\mathcal{O}_1$ limits are the strongest spin-independent absorption constraints reported for dark-matter masses below 23 MeV/$c^2$, so any model predicting such absorption in that mass range must now face these bounds.
- The $\mathcal{O}_4$ limits are the first constraints on spin-dependent absorption, opening the way to constrain magnetic-dipole-mediated dark matter that spin-independent searches cannot see.
- Because the signal peaks at $E_R \approx m_\chi^2/(2m_N)$ and cuts off abruptly once the threshold exceeds that energy, a threshold scan provides a discriminating signature that separates absorption from ordinary WIMP recoils.
- Applying the same framework to the latest dataset from a large xenon-based experiment projects sensitivity for both operators, showing that the method transfers to other detectors.
Reading between the lines
- If the velocity-operator substitution survives scrutiny, the same mapping should extend to other neutral-current absorption processes, so existing direct-detection data could be reinterpreted for sterile-neutrino or other light-fermion absorption models.
- The $\mathcal{O}_4$ limits rest on an unbenchmarked nuclear response; a first-principles calculation of the $^{19}\mathrm{F}$ spin-dependent response would either confirm the reported bounds or show they need rescaling, making it the sharpest test of the framework.
- The predicted mono-energetic peak at $E_R \approx m_\chi^2/(2m_N)$ suggests a multi-threshold search: if the same excess appears at thresholds below the peak and vanishes above it, that would be a distinctive absorption signature largely independent of the nuclear-response uncertainty.
- Because the peak recoil energy scales as $1/m_N$, even lighter targets such as deuterated liquids could extend absorption sensitivity to sub-MeV dark matter masses, a direction the paper does not explore.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This Letter extends the non-relativistic effective field theory (NREFT) used for elastic dark matter-nucleus scattering to the absorptive neutral-current process χN→νN. The authors state that the extension amounts to replacing the relative-velocity operator of the elastic theory with the average nucleon velocity, with the nuclear response functions modified accordingly, and they use a modified version of the WIMpy code to compute expected rates for the operators O1 (spin-independent) and O4 (spin-dependent). Reanalyzing the PICO-60 C3F8 bubble-chamber data with the profile-likelihood procedure established in earlier PICO analyses, the paper reports 90% C.L. upper limits on the absorption cross section σNC = mχ²/(4πΛ⁴): leading spin-independent limits for dark-matter masses below 23 MeV/c² and the first spin-dependent absorption limits (for O4) over roughly 12-27 MeV/c². The same framework is also applied to the most recent LZ dataset to produce sensitivity projections. The Letter concludes that bubble chambers are competitive probes of MeV-scale fermionic dark matter absorption.
Significance. If the theory input is sound, the result is significant for the direct-detection community: it opens a model-independent NREFT route to absorption searches, adds a spin-dependent channel to what has so far been a largely spin-independent program, and demonstrates the reach of bubble-chamber data at MeV-scale masses. The analysis has real strengths: it reuses the established PICO profile-likelihood machinery including efficiency surfaces and Seitz-threshold uncertainties; the O1 limits are shown in the same figure with CDEX-10, PandaX-4T, and Majorana Demonstrator limits, providing an external check of that channel; and no physics parameter is fitted to data and then re-labeled as a prediction. The weakness is that the new physics ingredient, the mapping from elastic to absorptive NREFT operators and the resulting nuclear responses, is asserted rather than derived, the modified code is not released, and the O4 channel has no external benchmark. The primary claims should therefore be regarded as plausible but not yet supported to the standard the journal requires.
major comments (3)
- [Absorption of Fermionic Dark Matter (Sec. II)] The central mapping from elastic NREFT operators to the absorptive process χN→νN is asserted, not derived. The text states that the only difference from the elastic case is the replacement of the velocity operator with the average nucleon velocity v_N^avg and that 'the nuclear responses provoked by each operator change,' but no explicit absorptive response functions are given anywhere in the Letter. This matters because for a light χ at rest the outgoing neutrino carries momentum p_ν ≈ m_χ (Eq. 3), so p_ν/m_χ is of order unity; a non-relativistic reduction of the leptonic current then produces spatial components that are not suppressed and are not captured by an operator basis built from 1, S_χ·S_N and v_N^avg ≈ q/(2m_N), which is tiny at these masses. If those O(1) terms feed into the nuclear response functions, the normalization of the O1 limits in Fig. 1 and the entire O4 limit in Fig. 2 would not correspond to the four-fermion interactions the paper claims to constrain. The authors cite the μ→e NREFT formalism of refs. [35,36], which was developed for exactly this kinematics, but they neither apply it nor compare with it. A derivation of the absorptive operators and of the resulting nuclear responses, or at least a quantitative argument showing that the neutrino-momentum terms are negligible, is required before the reported limits can be accepted.
- [Absorption of Fermionic Dark Matter and PICO-60 Data Analysis (Secs. II-III)] The modified WIMpy code and the nuclear response functions used for absorption are not made available, and there is no external benchmark for the O4 channel, so the central calculation is not reproducible as presented. The Letter states that the code 'has been modified to extend its applicability to absorptive interactions by incorporating the changes in velocity dependence and nuclear responses,' but the modified response functions are not listed and the code is not released. For O1 the comparison with CDEX-10, PandaX-4T, and MJD in Fig. 1 provides partial validation of the normalization, but for O4 there is no such check and, being a first limit, no consistency test exists. Given that the O4 limits are a headline result, the authors should provide the explicit response functions (at least for O1 and O4, including the fluorine spin-structure inputs) in an appendix or supplementary material, or release the modified code with documentation.
- [Results, Eq. (5)] The relation between the quoted cross section σNC = mχ²/(4πΛ⁴) and the event rate that enters the profile-likelihood analysis is not given, and the normalization conventions of the earlier experiments shown in Fig. 1 are not discussed. The Letter says the σ parameterization follows Dror et al. and other collaborations, but it does not state whether the CDEX-10, PandaX-4T, and MJD limits were recast to the same operator O1 and the same σ definition; if the conventions differ, the claim of leading limits below 23 MeV/c² cannot be assessed from the figure alone. The authors should state the rate formula used in the code (flux, nuclear response, and σ convention) and clarify the convention applied to each comparison dataset.
minor comments (6)
- [Absorption of Fermionic Dark Matter, Eq. (3)] Equation (3) is typeset ambiguously ('q = m_N m_χ + m_χ²/2 m_N + m_χ'); please rewrite it with clear parentheses or a displayed fraction, since the formula q ≈ m_χ for m_N ≫ m_χ is used to justify the kinematics of the whole analysis.
- [PICO-60 C3F8 Bubble Chamber and Data Analysis (Sec. III)] There is a grammatical typo: 'all enclosed within an stainless steel pressure vessel' should read 'a stainless steel pressure vessel.'
- [References] Reference [27] is listed as 'Zhang et al. (PandaX Collaboration)' without the lead author's initials; please supply the full author list or at least the first author's initials for consistency with the other entries.
- [Results and Conclusions] The sentence 'The reported limits remain consistent with the absence of a dark matter signal' is the only place where the observed event counts are connected to the limit-setting; a sentence recalling that eight single- and multi-bubble events were observed and taken as background under the maximally conservative assumption would make the section self-contained.
- [Fig. 2 caption] The LZ sensitivity projection in Fig. 2 is not described: please specify the dataset, exposure, energy threshold, and signal region used in the recast so that the comparison with the PICO-60 limit is meaningful.
- [Absorption of Fermionic Dark Matter (Sec. II)] The phrase 'a unique peaked signature at ER ≈ mχ²/(2mN)' is imprecise: for a cold χ the recoil spectrum is a mono-energetic line, broadened by the target momentum distribution and detector effects; consider stating this explicitly rather than 'peaked.'
Circularity Check
No significant circularity: limits are computed from observed PICO-60 event counts; the NREFT velocity substitution is an input assumption, not a quantity fitted to or derived from the data.
full rationale
The reported absorption limits are obtained by a standard profile-likelihood exclusion calculation on PICO-60's observed single- and multi-bubble events, with no dark-matter cross-section parameter fitted to the data and then re-labeled as a prediction. The NREFT operators O1 and O4 are introduced as theoretical inputs, and the paper's only adaptation is the replacement of the relative-velocity operator by the average nucleon velocity. That replacement is stated as an assumption ("since it is not expected that the velocity of the incoming dark matter particle contributes meaningfully to the scattering") rather than derived from the data, so any concern about its correctness is a physics-modeling risk, not a circular reduction of the analysis to its inputs. The citations to prior PICO analyses [18, 42, 38, 46] are methodological (limit-setting procedure, calibrations, WIMpy modification) and do not themselves carry the physics claim; the framework citations [34, 35, 36] are external to the collaboration. No equation in the paper defines the predicted rate in terms of the measured limit, nor fits a parameter that is then called a prediction. Accordingly, under the stated circularity criteria, no load-bearing step reduces to its own input.
Assumptions & free parameters
assumptions (4)
- domain assumption The non-relativistic expansion is valid for cold dark matter with mass comparable to or lighter than the nucleus, so the DM momentum is negligible and the neutrino is relativistic.
- domain assumption Only the leading operators O1 and O4 need be considered; in matching from four-fermion operators, one operator tends to dominate.
- ad hoc to paper The nuclear response functions for absorption are obtained from the elastic NREFT by replacing the relative velocity operator with the average nucleon velocity operator.
- domain assumption Standard halo model with local density 0.3 GeV/c2/cm3 and standard velocity distribution.
Cite this review
Pith. "Pith review of Absorption of Fermionic Dark Matter in the PICO-60 C$_{3}$F$_{8}$ Bubble Chamber." pith.science (2026). https://pith.science/paper/XJBL3BVA
@misc{pith2026250413089,
author = {Pith},
title = {Pith review of: Absorption of Fermionic Dark Matter in the PICO-60 C$_3$F$_8$ Bubble Chamber},
year = {2026},
howpublished = {\url{https://pith.science/paper/XJBL3BVA}},
note = {Machine review of arXiv:2504.13089}
}
abstract
Fermionic dark matter absorption on nuclear targets via neutral current interactions is explored using a non-relativistic effective field theory framework. An analysis of data from the PICO-60 C$_{3}$F$_{8}$ bubble chamber sets leading constraints on spin-independent absorption for dark matter masses below 23 MeV/$\textit{c}^2$ and establishes the first limits on spin-dependent absorptive interactions. These results demonstrate the sensitivity of bubble chambers to low-mass dark matter and underscore the importance of absorption searches in expanding the parameter space of direct detection experiments.
Figures
Forward citations
Cited by 1 Pith paper
-
Supernova cooling from neutrino-devouring dark matter
Supernova cooling excludes fermionic dark matter produced by neutrino scattering down to cross sections of about 10^-58 cm^2 for electrons and 10^-56 cm^2 for nucleons across the keV to 100 MeV mass range.
Reference graph
Works this paper leans on
-
[1]
Einasto, A
J. Einasto, A. Kaasik, and E. Saar, Nature250, 309 (1974)
1974
-
[2]
J. Ostriker, P. Peebles, and A. Yahil, The Astrophysical Journal193, L1 (1974)
work page 1974
-
[3]
J. P. Ostriker and P. J. E. Peebles, Astrophys. J.186, 467 (1973)
work page 1973
- [4]
-
[5]
N. A. Bahcall and X. Fan, Proc. Nat. Acad. Sci.95, 5956 (1998), arXiv:astro-ph/9804082
work page Pith review arXiv 1998
-
[6]
N. A. Bahcall and A. Kulier, Mon. Not. Roy. Astron. Soc. 439, 2505 (2014), arXiv:1310.0022 [astro-ph.CO]
work page Pith review arXiv 2014
-
[7]
E. S. Sheldon et al. (SDSS), Astrophys. J.703, 2232 (2009), arXiv:0709.1162 [astro-ph]
work page Pith review arXiv 2009
-
[8]
R. Mandelbaum, A. Slosar, T. Baldauf, U. Seljak, C. M. Hirata, R. Nakajima, R. Reyes, and R. E. Smith, Mon. Not.Roy.Astron.Soc.432,1544(2013),arXiv:1207.1120 [astro-ph.CO]
arXiv 2013
Show all 49 references
-
[9]
C. L. Bennettet al., The Astrophysical Journal Supple- ment Series208, 20 (2013)
2013
-
[10]
P. A. R. Adeet al. (Planck Collaboration), Astron. As- trophys.571, A1. 48 p (2013), arXiv:1303.5062
2013 arXiv
-
[11]
Menget al
Y. Menget al. (PandaX-4T Collaboration), Phys. Rev. Lett.127, 261802 (2021)
2021
-
[12]
121, 111302 (2018), arXiv:1805.12562 [astro-ph.CO]
E.Aprile etal.(XENONCollaboration),Phys.Rev.Lett. 121, 111302 (2018), arXiv:1805.12562 [astro-ph.CO]
2018 arXiv
-
[13]
Hornet al
M. Hornet al. (LUX), Nucl. Instrum. Meth. A784, 504 (2015)
2015
-
[14]
Adhikari et al
P. Adhikari et al. (DEAP Collaboration), Phys. Rev. Lett.128, 011801 (2022)
2022
-
[15]
Ahmed et al
Z. Ahmed et al. (CDMS-II), Science327, 1619 (2010), arXiv:0912.3592 [astro-ph.CO]
2010 arXiv
-
[16]
Kluck et al
H. Kluck et al. (CRESST Collaboration), Journal of Physics: Conference Series1468, 012038 (2020)
2020
- [17]
-
[18]
Amoleetal
C. Amoleetal. (PICO Collaboration), Phys. Rev. D100, 022001 (2019)
2019
-
[19]
S. E. Vahsen et al., (2020), arXiv:2008.12587 [physics.ins-det]
2020 arXiv
-
[20]
J. A. Dror, G. Elor, and R. Mcgehee, Phys. Rev. Lett. 124, 18 (2020), arXiv:1905.12635 [hep-ph]
2020 arXiv
-
[21]
J. A. Dror, G. Elor, and R. McGehee, Journal of High Energy Physics2020, 134 (2020)
2020
-
[22]
J. A. Dror, G. Elor, R. McGehee, and T.-T. Yu, Physical Review D103, 035001 (2021)
2021
-
[23]
Ge, X.-G
S.-F. Ge, X.-G. He, X.-D. Ma, and J. Sheng, Journal of High Energy Physics2022, 191 (2022)
2022
-
[24]
Gu et al
L. Gu et al. (PandaX Collaboration), Phys. Rev. Lett. 129, 161803 (2022)
2022
-
[25]
Lett.132, 041001 (2024)
I.J.Arnquist etal.(MajoranaCollaboration),Phys.Rev. Lett.132, 041001 (2024)
2024
-
[26]
W. H. Daiet al. (CDEX Collaboration), Phys. Rev. Lett. 129, 221802 (2022)
2022
-
[27]
(PandaX Collaboration), Phys
Zhang et al. (PandaX Collaboration), Phys. Rev. Lett. 129, 161804 (2022)
2022
-
[28]
P. D. l. T. Luque, S. Balaji, and J. Silk, Phys. Rev. Lett. 134, 101001 (2025)
2025
-
[29]
T. Li, J. Liao, and R.-J. Zhang, Journal of High Energy Physics2022, 71 (2022)
2022
-
[30]
Beltrán, P
R. Beltrán, P. D. Bolton, F. F. Deppisch, C. Hati, and M. Hirsch, Journal of High Energy Physics2024, 153 (2024)
2024
-
[31]
Zhang and W
Y. Zhang and W. Liu, Physical Review D107, 095031 (2023)
2023
-
[32]
Brdar, A
V. Brdar, A. De Gouvêa, Y.-Y. Li, and P. A. Machado, Physical Review D107, 073005 (2023)
2023
-
[33]
O. G. Miranda, D. K. Papoulias, O. Sanders, M. Tórtola, and J. W. F. Valle, Journal of High Energy Physics2021, 191 (2021)
2021
-
[34]
A. L. Fitzpatrick, W. Haxton, E. Katz, N. Lubbers, and Y. Xu, J. Cosmol. Astropart. Phys.2013, 004 (2013)
2013
-
[35]
Haxton, K
W. Haxton, K. McElvain, T. Menzo, E. Rule, and J. Zu- pan, JHEP11, 076 (2024), arXiv:2406.13818 [hep-ph]
2024 arXiv
-
[36]
W. C. Haxton, E. Rule, K. McElvain, and M. J. Ramsey-Musolf, Phys. Rev. C107, 035504 (2023), arXiv:2208.07945 [nucl-th]
2023 arXiv
-
[37]
Kavanagh, Physical Review D92(2015), 10.1103/PhysRevD.92.023513
B. Kavanagh, Physical Review D92(2015), 10.1103/PhysRevD.92.023513
2015 doi
-
[38]
Ali et al
B. Ali et al. (PICO Collaboration), Phys. Rev. D106, 042004 (2022)
2022
-
[39]
N. J. T. Smith, Eur. Phys. J. Plus127, 108 (2012)
2012
-
[40]
Amoleetal
C. Amoleetal. (PICO Collaboration), Phys. Rev. D100, 082006 (2019)
2019
-
[41]
Seitz, Physics of Fluids1, 2 (1958)
F. Seitz, Physics of Fluids1, 2 (1958)
1958
-
[42]
Amoleet al
C. Amoleet al. (PICO Collaboration), Phys. Rev. Lett. 118, 251301 (2017). 7
2017
-
[43]
Durnford and M.-C
D. Durnford and M.-C. Piro, Journal of Instrumentation 17, C01030 (2022)
2022
-
[44]
Cowan, K
G. Cowan, K. Cranmer, E. Gross, and O. Vitells, The European Physical Journal C71(2011), 10.1140/epjc/s10052-011-1554-0
2011 doi
-
[45]
Foreman-Mackey, D
D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Goodman, Publ. Astron. Soc. Pac.125, 306 (2013), arXiv:1202.3665 [astro-ph.IM]
2013 arXiv
-
[46]
Ali et al
B. Ali et al. (PICO Collaboration), Phys. Rev. D106, 122003 (2022)
2022
-
[47]
Lewin and P
J. Lewin and P. Smith, Astroparticle Physics6, 87 (1996)
1996
-
[48]
Website: https://www.alliancecan.ca
-
[133]
[40], respectively. The detectors are operated such that a bubble forms only when energy exceeding a ther- modynamic threshold is deposited within a critical ra- dius—this is known as the Seitz threshold [41]. Bubble nucleation is therefore a threshold process governed by the ...
2016
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.