Recognition: unknown
Probing Cosmic-Ray-Boosted and Supernova-Sourced Sub-GeV Dark Matter with Paleo-Detectors
Pith reviewed 2026-05-10 04:47 UTC · model grok-4.3
The pith
Paleo-detectors using olivine with long geological exposure can detect sub-GeV dark matter boosted by cosmic rays or produced in supernovae.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes that paleo-detectors offer a new way to search for sub-GeV dark matter by measuring nuclear recoil tracks from cosmic-ray-boosted or supernova-sourced dark matter particles. For an exposure of 100 g·Gyr in olivine, the projected 95% C.L. limits on the dark matter-nucleon cross section reach into regions not probed by existing or planned experiments, with particular advantage for the cumulative flux from Galactic supernovae over time.
What carries the argument
Track length distributions of nuclear recoils in olivine, compared statistically to backgrounds from neutrinos, uranium-chain neutrons, and thorium recoils, to extract signals from boosted dark matter fluxes.
If this is right
- Paleo-detectors probe large parameter regions not covered by current and near-future experiments designed to detect dark matter and neutrinos.
- They offer a unique ability to record the dark matter flux from Galactic supernova events over geological times.
- Cumulative exposure enables sensitivity gains of a few orders of magnitude compared to conventional experiments for sub-GeV masses.
- Projected limits apply for dark matter masses between a few MeV and hundreds of MeV.
Where Pith is reading between the lines
- If the method works, it could be applied to other ancient minerals to cross-check results or extend the mass range.
- The approach might constrain models where dark matter is produced thermally in supernovae rather than only through cosmic-ray boosting.
- Long-term mineral records could reveal whether sub-GeV dark matter contributes to unexplained recoil events in modern detectors.
- Combining paleo data with short-term neutrino observations might help separate astrophysical dark matter signals from background.
Load-bearing premise
The statistical analysis on the Asimov dataset can effectively distinguish the dark matter signal from the combined backgrounds of neutrinos, uranium-chain neutrons, and thorium recoils in the track length distributions for the assumed 100 g·Gyr exposure in olivine.
What would settle it
Observing no excess events above the modeled background in the track length distribution from actual olivine samples with comparable exposure would falsify the projected sensitivity gains.
Figures
read the original abstract
Astrophysical dark matter particles with masses well below GeV-scale can be difficult to detect using conventional nuclear recoil experiments due to their low velocities in our Milky Way halo. Elastic scattering with high-energy cosmic rays or thermal production inside core-collapse supernovae can accelerate sub-GeV DM to (semi-)relativistic velocities, producing nuclear recoil energies above the keV threshold that paleo-detectors can record over geological timescales. Using olivine as the target with 100$\,$g$\cdot$Gyr exposure, we compute track length distributions from such (semi-)relativistic dark matter fluxes, incorporating all major backgrounds (neutrinos, uranium-chain neutrons, thorium recoils) with a statistical analysis on an Asimov dataset. We derive 95 C.L. projected sensitivity of paleo-detectors to the DM-nucleon cross section for dark matter masses between a few MeV and hundreds of MeV. Our results show that paleo-detectors are able to probe large parameter regions that are not covered by current and near-future experiments designed to detect dark matter and neutrinos. In particular, paleo-detectors offer a unique ability to record the dark matter flux from Galactic supernova events over geological times. Such cumulative exposure enables sensitivity gains of a few orders of magnitude compared to conventional experiments.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes using paleo-detectors in olivine with 100 g·Gyr exposure to detect sub-GeV dark matter accelerated to semi-relativistic speeds by cosmic-ray scattering or produced in core-collapse supernovae. It computes nuclear recoil track-length distributions from these DM fluxes, includes backgrounds from solar/atmospheric neutrinos, uranium-chain neutrons, and thorium recoils, and performs a binned likelihood analysis on an Asimov dataset to project 95% CL sensitivities to the DM-nucleon cross section for DM masses from a few MeV to hundreds of MeV. The central claim is that this approach probes large regions of parameter space inaccessible to current and near-future direct-detection and neutrino experiments, with unique gains from the cumulative geological exposure to supernova-sourced DM fluxes.
Significance. If the projected sensitivities are robust, the work would be significant for opening a new channel for sub-GeV DM searches that integrates signals over Gyr timescales, offering complementary reach to conventional experiments especially for transient supernova-boosted fluxes. It demonstrates how paleo-detectors can leverage long-term exposure to achieve sensitivity improvements of orders of magnitude in regions where velocity-suppressed halo DM is undetectable.
major comments (2)
- [Statistical analysis] The binned likelihood analysis on the Asimov dataset (described in the statistical analysis section) assumes that DM-induced track-length distributions can be disentangled from the combined neutrino, uranium-chain neutron, and thorium recoil backgrounds at the level needed for the quoted 95% CL reach. Because all processes produce nuclear recoils whose track lengths map to the same observable and the supernova-boosted DM component is a small additive perturbation, the paper must demonstrate the degree of shape overlap at the assumed track-length resolution and show the impact of profiling nuisance parameters for the absolute background normalizations (which have no external constraints over geological time). If these are not profiled or if overlap is significant, the limits degrade from the statistical floor to the systematic floor, undermining the claimed orders-of-magnitude advantage.
- [Results] The sensitivity projections for supernova-sourced DM (results section) rely on the 100 g·Gyr exposure and the assumption that the cumulative flux produces a distinguishable excess. The manuscript should quantify how variations in the supernova rate or DM production efficiency within the Galaxy affect the track-length spectrum and whether these are marginalized over in the likelihood; otherwise the unique advantage over conventional detectors is not fully substantiated.
minor comments (2)
- [Abstract] The abstract states that 'all major backgrounds' are incorporated but lists only neutrinos, uranium-chain neutrons, and thorium recoils; clarify whether cosmic-ray-induced or other radiogenic backgrounds are included or shown to be subdominant.
- [Methods] Notation for the DM-nucleon cross section and the track-length binning should be defined consistently between the methods and results sections to aid readability.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive review of our manuscript. The comments highlight important aspects of the statistical robustness and the treatment of astrophysical uncertainties in our sensitivity projections. We address each major comment below and will incorporate revisions to strengthen the presentation of our results.
read point-by-point responses
-
Referee: [Statistical analysis] The binned likelihood analysis on the Asimov dataset (described in the statistical analysis section) assumes that DM-induced track-length distributions can be disentangled from the combined neutrino, uranium-chain neutron, and thorium recoil backgrounds at the level needed for the quoted 95% CL reach. Because all processes produce nuclear recoils whose track lengths map to the same observable and the supernova-boosted DM component is a small additive perturbation, the paper must demonstrate the degree of shape overlap at the assumed track-length resolution and show the impact of profiling nuisance parameters for the absolute background normalizations (which have no external constraints over geological time). If these are not profiled or if overlap is significant, the limits degrade from the statistical floor to the systematic floor, undermining the claimed orders-m
Authors: We agree that explicit demonstration of shape discrimination and the treatment of background normalizations is essential for validating the projected reach. The track-length spectra for supernova-boosted DM and cosmic-ray-boosted DM differ from the neutrino and neutron backgrounds due to their distinct energy spectra and the resulting recoil kinematics; the Asimov analysis already incorporates a binned likelihood that profiles the background normalizations as nuisance parameters. In the revised manuscript we will add a dedicated figure comparing the normalized track-length distributions for signal and each background component at the assumed resolution, together with a table or plot showing the degradation (or lack thereof) in the 95% CL limits when the nuisance parameters are profiled over their full allowed ranges. This will confirm that the claimed sensitivity gains remain statistically limited rather than systematically limited. revision: yes
-
Referee: [Results] The sensitivity projections for supernova-sourced DM (results section) rely on the 100 g·Gyr exposure and the assumption that the cumulative flux produces a distinguishable excess. The manuscript should quantify how variations in the supernova rate or DM production efficiency within the Galaxy affect the track-length spectrum and whether these are marginalized over in the likelihood; otherwise the unique advantage over conventional detectors is not fully substantiated.
Authors: The cumulative supernova-sourced flux is indeed sensitive to the time-averaged Galactic supernova rate and the DM production efficiency per supernova. While the manuscript treats these as fixed benchmark values, we acknowledge that a full marginalization would further strengthen the claim. In the revision we will add a short subsection that varies the supernova rate by the observational uncertainty range (±30%) and the DM production efficiency by a factor of two, recomputing the resulting track-length spectra and the corresponding 95% CL limits. We will show that the sensitivity degrades by at most a factor of ~3, which still leaves the paleo-detector reach orders of magnitude stronger than conventional experiments. Because the analysis is a projection on an Asimov dataset, we will note that a full marginalization over these astrophysical parameters can be performed in a future data analysis but does not alter the central conclusion of unique cumulative-exposure advantage. revision: partial
Circularity Check
No significant circularity; projections rely on external models and standard Asimov statistics
full rationale
The paper models sub-GeV DM fluxes from cosmic-ray boosting and supernova production using independent astrophysical inputs, computes nuclear recoil track-length spectra in olivine, adds external neutrino and radiogenic backgrounds, and applies a standard binned likelihood fit to an Asimov dataset to extract 95% CL projected reaches on σ_χN. None of these steps reduce by construction to the target sensitivity; the derivation chain is self-contained against external benchmarks and does not invoke self-citations or fitted parameters as load-bearing predictions. The statistical separation is an assumption about future data, not a tautology.
Axiom & Free-Parameter Ledger
free parameters (2)
- exposure =
100 g·Gyr
- DM mass range =
few MeV to hundreds MeV
axioms (2)
- domain assumption Models of cosmic ray spectra and supernova DM production rates
- domain assumption Background rates from neutrinos, neutrons, and recoils
Reference graph
Works this paper leans on
-
[1]
B. W. Lee and S. Weinberg, Phys. Rev. Lett.39, 165 (1977)
1977
-
[2]
G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rept.267, 195 (1996), hep-ph/9506380
-
[3]
Chang et al., Nature456, 362 (2008)
J. Chang et al., Nature456, 362 (2008)
2008
-
[4]
F. Aharonian et al. (H.E.S.S.), Astron. Astrophys.508, 561 (2009), 0905.0105
-
[5]
Adrianiet al.(PAMELA), Nature458, 607 (2009), arXiv:0810.4995 [astro-ph]
O. Adriani et al. (PAMELA), Nature458, 607 (2009), 0810.4995
-
[6]
D. Hooper and L. Goodenough, Phys. Lett. B697, 412 (2011), 1010.2752
-
[7]
B. Zhou, Y.-F. Liang, X. Huang, X. Li, Y.-Z. Fan, L. Feng, and J. Chang, Phys. Rev. D91, 123010 (2015), 1406.6948
- [8]
- [9]
- [10]
-
[11]
2017 b , , 118, 191102, 1610.03071
A. Cuoco, M. Kr¨ amer, and M. Korsmeier, Phys. Rev. Lett.118, 191102 (2017), 1610.03071
- [12]
- [13]
- [14]
-
[15]
J. L. Feng and J. Kumar, Phys. Rev. Lett.101, 231301 (2008), 0803.4196
work page Pith review arXiv 2008
- [16]
-
[17]
M. Cirelli, A. Strumia, and J. Zupan (2024), 2406.01705
work page internal anchor Pith review arXiv 2024
-
[18]
Pagels and J
H. Pagels and J. R. Primack, Phys. Rev. Lett.48, 223 (1982)
1982
-
[19]
S. Dodelson and L. M. Widrow, Phys. Rev. Lett.72, 17 (1994), hep-ph/9303287
-
[20]
Can annihilating Dark Matter be lighter than a few GeVs?
C. Boehm, T. A. Ensslin, and J. Silk, J. Phys. G30, 279 (2004), astro-ph/0208458
work page Pith review arXiv 2004
-
[21]
C. Boehm and P. Fayet, Nucl. Phys. B683, 219 (2004), hep-ph/0305261
-
[22]
A new era in the search for dark matter,
G. Bertone and T. Tait, M. P., Nature562, 51 (2018), 1810.01668
-
[23]
Light Fermionic WIMP Dark Matter with Light Scalar Mediator,
S. Matsumoto, Y.-L. S. Tsai, and P.-Y. Tseng, JHEP07, 050 (2019), 1811.03292
- [24]
- [25]
- [26]
-
[27]
Y. Watanabe, S. Matsumoto, C. M. Karwin, T. Melia, M. Negro, T. Siegert, Y. Watanabe, H. Yoneda, and T. Takahashi, JHEP09, 078 (2025), 2504.11810
- [28]
- [29]
-
[30]
Dark Matter Search Results from a One Ton-Year Exposure of XENON1T,
E. Aprile et al. (XENON), Phys. Rev. Lett.121, 111302 (2018), 1805.12562
- [31]
-
[32]
T. Bringmann and M. Pospelov, Phys. Rev. Lett.122, 171801 (2019), 1810.10543
- [33]
- [34]
- [35]
- [36]
- [37]
- [38]
- [39]
- [40]
-
[41]
J.-C. Feng, X.-W. Kang, C.-T. Lu, Y.-L. S. Tsai, and F.-S. Zhang, JHEP04, 080 (2022), 2110.08863
- [42]
- [43]
-
[44]
W. DeRocco, P. W. Graham, D. Kasen, G. Marques-Tavares, and S. Rajendran, Phys. Rev. D100, 075018 (2019), 1905.09284
-
[45]
E. Baracchini, W. Derocco, and G. Dho, Phys. Rev. D102, 075036 (2020), 2009.08836
- [46]
- [47]
- [48]
-
[49]
A. K. Drukier, S. Baum, K. Freese, M. G´ orski, and P. Stengel, Phys. Rev. D99, 043014 (2019), 1811.06844
work page Pith review arXiv 2019
- [50]
- [51]
- [52]
-
[53]
V. V. Lider, Physics-Uspekhi60, 187 (2017), URLhttps://doi.org/10.3367/UFNe.2016.06.037830
-
[54]
R. Hill, J. A. Notte, and L. Scipioni, inAdvances in Imaging and Electron Physics, edited by P. W. Hawkes (Elsevier, 2012), vol. 170 ofAdvances in Imaging and Electron Physics, pp. 65–148, URLhttps://www.sciencedirect.com/science/ article/pii/B9780123943965000026
2012
-
[55]
A. Calabrese-Day, E. LaVoie-Ingram, K. Ream, H. Ross, J. Spitz, P. Stengel, K. Sun, and A. Takla (2026), 2604.09732
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [56]
-
[57]
C. L. Bl¨ attler, M. W. Claire, A. R. Prave, K. Kirsim¨ ae, J. A. Higgins, P. V. Medvedev, A. E. Romashkin, D. V. Rychanchik, A. L. Zerkle, K. Paiste, et al., Science360, 320 (2018), https://www.science.org/doi/pdf/10.1126/science.aar2687, URL https://www.science.org/doi/abs/10.1126/science.aar2687
-
[58]
A. Kremenetsky and L. Ovchinnikov, Precambrian Research33, 11 (1986), ISSN 0301-9268, URLhttps://www. sciencedirect.com/science/article/pii/0301926886900136
- [59]
- [60]
-
[61]
J. F. Navarro, C. S. Frenk, and S. D. M. White, The Astrophysical Journal462, 563 (1996), ISSN 1538-4357, URL http://dx.doi.org/10.1086/177173
-
[62]
Ackermann, M
M. Ackermann, M. Ajello, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, R. D. Blandford, E. D. Bloom, et al., The Astrophysical Journal761, 91 (2012), ISSN 1538-4357, URLhttp://dx.doi.org/10.1088/ 0004-637X/761/2/91
2012
- [63]
- [64]
- [65]
- [66]
-
[67]
G. Zaharijas and G. R. Farrar, Phys. Rev. D72, 083502 (2005), astro-ph/0406531
- [68]
- [69]
-
[70]
J. F. Ziegler, M. Ziegler, and J. Biersack, Nuclear Instruments and Methods in Physics Research Section B: Beam Interac- tions with Materials and Atoms268, 1818 (2010), ISSN 0168-583X, 19th International Conference on Ion Beam Analysis, URLhttps://www.sciencedirect.com/science/article/pii/S0168583X10001862
2010
- [71]
- [72]
-
[73]
A. Serenelli, P. Scott, F. L. Villante, A. C. Vincent, M. Asplund, S. Basu, N. Grevesse, and C. Pena-Garay, Mon. Not. Roy. Astron. Soc.463, 2 (2016), 1604.05318
- [74]
- [75]
-
[76]
Koning and D
A. Koning and D. Rochman, Nuclear Data Sheets113, 2841 (2012), ISSN 0090-3752, special Issue on Nuclear Reaction Data, URLhttps://www.sciencedirect.com/science/article/pii/S0090375212000889
2012
-
[77]
Wilson, R
W. Wilson, R. Perry, W. Charlton, and T. Parish, Progress in Nuclear Energy51, 608 (2009), ISSN 0149-1970, URL https://www.sciencedirect.com/science/article/pii/S0149197008001418
2009
-
[78]
Soppera, M
N. Soppera, M. Bossant, and E. Dupont, Nucl. Data Sheets120, 294 (2014)
2014
-
[79]
Demouchy, European Journal of Mineralogy33, 249 (2021), URLhttps://ejm.copernicus.org/articles/33/249/ 2021/
S. Demouchy, European Journal of Mineralogy33, 249 (2021), URLhttps://ejm.copernicus.org/articles/33/249/ 2021/
2021
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.