pith. machine review for the scientific record. sign in

arxiv: 2605.13659 · v1 · submitted 2026-05-13 · ✦ hep-ph · astro-ph.CO· hep-ex

Recognition: unknown

Calorimetric approach to paleo-detection of dark matter

Authors on Pith no claims yet

Pith reviewed 2026-05-14 17:54 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.COhep-ex
keywords paleo-detectorsdark matternuclear recoilslattice vacanciesbackground discriminationolivineSRIMWIMP sensitivity
0
0 comments X

The pith

Combining track lengths with lattice vacancy counts in mineral samples lets paleo-detectors distinguish dark matter recoils from neutron backgrounds and reach cross-section sensitivities of order 10^{-48} cm^{2}.

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

The paper introduces a calorimetric readout for paleo-detectors that counts stable lattice vacancies created by each nuclear recoil in addition to measuring track length. Simulations in olivine show that the vacancy count alone matches the sensitivity of existing track-only methods for a 100 gGyr exposure. When the two observables are combined they supply an event-by-event proxy for the energy-loss rate dE/dx, which identifies the recoiling nucleus. Because dark-matter scattering cross sections grow with nuclear mass squared while neutron scattering does not, this identification suppresses the dominant neutron background by more than an order of magnitude at moderate WIMP masses. The resulting sensitivity becomes comparable to that projected for future direct-detection experiments.

Core claim

Using full-cascade SRIM simulations the authors show that a vacancy-only readout already reaches sensitivities similar to track-only analyses; combining the two observables supplies a proxy for |dE/dx| and therefore for the recoiling nuclear species, suppressing the neutron background by more than an order of magnitude and yielding spin-independent dark-matter–nucleon cross sections of order 10^{-48} cm^{2} at WIMP masses of a few tens of GeV for a 100 gGyr exposure.

What carries the argument

Dual readout of track length plus number of stable lattice vacancies, serving as a per-event proxy for |dE/dx| to identify the recoiling nucleus.

Load-bearing premise

Full-cascade SRIM simulations correctly predict the number of stable lattice vacancies created by each recoil and these vacancies remain measurable as a reliable per-event observable in olivine.

What would settle it

An experimental calibration showing that the measured vacancy count in olivine does not scale with |dE/dx| as predicted by SRIM, or that the combined observables fail to suppress the neutron background in a controlled sample.

Figures

Figures reproduced from arXiv: 2605.13659 by Patrick Huber, Samuel Hedges.

Figure 1
Figure 1. Figure 1: Distribution of track length (left panel) and number of vacancies (right panel) as a function of primary recoil energy, [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Projection into the xy-plane of two events simulated with SRIM at approximately 10 keV primary recoil energy. Filled disks show the result in quick mode, open circles show the result in full cascade mode. Each colored arc corresponds to one of the three track-length definitions discussed in the text: Euclidean first-to-last distance, arc length along the time-ordered cascade, and diameter (largest pairwise… view at source ↗
Figure 3
Figure 3. Figure 3: Track length (left panel) and vacancy count (right panel) distributions for the dark matter signal at [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: 90% C.L. sensitivity limits on the spin-independent dark matter–nucleon cross section as a function of dark matter [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: 90% C.L. sensitivity limits on the spin-independent dark matter–nucleon cross section as a function of dark matter [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Joint two-dimensional histograms in the (vacancy count, track length) plane for the neutron background (blue) and [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: 90% C.L. sensitivity limits on the spin-independent dark matter–nucleon cross section as a function of dark matter [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Neutrino and neutron fog for the 100 g ·1 Gyr olivine benchmark, for the track-only (left panel), vacancy-only (middle panel), and combined (right panel) analyses. Colored regions show contours of the fog parameter n = (d ln σ/d ln N) −1 (Eq. 4), and the black lines mark the corresponding 90% C.L. sensitivities. The labels identify the dominant background behind each island: neutrons (“neut.”), solar neutr… view at source ↗
Figure 9
Figure 9. Figure 9: Mass-versus-resolution landscape for paleo-detector readout. SPIM (selective-plane illumination microscopy) reaches [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
read the original abstract

We present the first paleo-detector dark matter sensitivity analysis based on a calorimetric readout, in which the number of stable lattice vacancies produced by each nuclear recoil is used as a per-event observable complementary to the track length. Using full-cascade SRIM simulations in olivine, we compute the expected sensitivity for a 100 gGyr exposure. We find that a vacancy-only readout reaches a sensitivity envelope very similar to that of state-of-the-art track-only analyses. The combination of the two observables provides an event-by-event proxy for |dE/dx| and hence for the recoiling nuclear species. Since the neutron-nucleus cross section is approximately flat in nuclear mass while the dark-matter--nucleus cross section scales as $A^2$, this discrimination suppresses the dominant neutron background by more than an order of magnitude at moderate dark matter masses. The combined-analysis sensitivity reaches spin-independent dark-matter--nucleon cross sections of order $10^{-48}\,\mathrm{cm}^2$ at WIMP masses of a few tens of GeV, comparable to future direct detection experiments. A two-stage readout combining selective-plane illumination microscopy with scanning electron microscopy is identified as a path to making a 100 g-scale analysis plausible.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a calorimetric paleo-detection approach for dark matter in which the number of stable lattice vacancies produced by nuclear recoils in olivine serves as a per-event observable complementary to track length. Using full-cascade SRIM simulations, the authors compute sensitivities for a 100 gGyr exposure and show that combining the two observables provides an event-by-event proxy for dE/dx (and thus recoiling nuclear species), enabling neutron background suppression via the differing A-dependence of DM-nucleus versus neutron-nucleus cross sections. They project spin-independent DM-nucleon cross-section sensitivities of order 10^{-48} cm² at WIMP masses of a few tens of GeV, comparable to future direct-detection experiments, and identify a two-stage microscopy readout as a practical path forward.

Significance. If the SRIM-based vacancy predictions prove accurate, the combined observable would constitute a genuine advance in paleo-detection by adding background discrimination power without requiring larger exposures. The explicit identification of a feasible readout technology strengthens the practical relevance of the projected reach.

major comments (2)
  1. [SRIM simulation description and sensitivity calculation] The central sensitivity projections rest on SRIM predictions of stable lattice vacancy counts in olivine; no experimental calibration, cross-check, or uncertainty band for vacancy production is presented, leaving the calorimetric observable unvalidated at the level required for the claimed background rejection.
  2. [Background discrimination and combined-analysis results] The order-of-magnitude neutron suppression is derived from the A² scaling of the DM cross section versus the flatter neutron cross section, but the paper does not quantify the overlap between signal and background distributions in the combined (vacancy, track-length) plane or the resulting signal efficiency after cuts.
minor comments (2)
  1. [Results section] The abstract states that the vacancy-only readout reaches sensitivities 'very similar' to track-only analyses; the main text should include a direct side-by-side comparison plot or table with numerical values.
  2. [Method section] Notation for the vacancy observable (e.g., whether it is total vacancies or vacancies per unit track length) should be defined explicitly at first use and used consistently.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful review and constructive feedback. We address each major comment below and have revised the manuscript to strengthen the presentation of the SRIM-based results and the background discrimination analysis.

read point-by-point responses
  1. Referee: The central sensitivity projections rest on SRIM predictions of stable lattice vacancies in olivine; no experimental calibration, cross-check, or uncertainty band for vacancy production is presented, leaving the calorimetric observable unvalidated at the level required for the claimed background rejection.

    Authors: We agree that experimental calibration would provide the strongest validation of the vacancy observable. As this work is a theoretical sensitivity study, we have added a new subsection reviewing existing SRIM validation benchmarks for similar mineral targets in the literature, together with conservative uncertainty bands obtained by varying key simulation parameters such as displacement threshold energies. Within these bands the projected neutron suppression and sensitivity reach remain stable, although we explicitly note that dedicated experimental calibration would be a valuable next step. revision: partial

  2. Referee: The order-of-magnitude neutron suppression is derived from the A² scaling of the DM cross section versus the flatter neutron cross section, but the paper does not quantify the overlap between signal and background distributions in the combined (vacancy, track-length) plane or the resulting signal efficiency after cuts.

    Authors: We have performed the requested quantification and added new figures and text to the revised manuscript. The two-dimensional distributions of signal and neutron background events in the (vacancy count, track length) plane are now shown explicitly, together with the optimized cut contours. The resulting signal efficiency after cuts is 65–80 % across the relevant WIMP mass range while preserving the order-of-magnitude neutron suppression, confirming that the quoted sensitivity is achievable with realistic selection efficiencies. revision: yes

Circularity Check

0 steps flagged

No significant circularity; relies on external SRIM simulations

full rationale

The paper computes sensitivity using full-cascade SRIM simulations (external code) to predict stable lattice vacancies in olivine as a per-event observable, combined with standard nuclear recoil cross-section scalings (A^2 for DM-nucleus, approximately flat for neutrons). No derivation step reduces by the paper's own equations to a quantity fitted from its own data, and no self-citation is load-bearing for the central sensitivity claim. The result is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of SRIM modeling for vacancy production and on the known scaling of dark-matter versus neutron cross sections; no new free parameters are fitted and no new entities are postulated.

axioms (2)
  • domain assumption SRIM simulations accurately model the production of stable lattice vacancies by nuclear recoils in olivine
    Invoked to compute the expected vacancy number per recoil event for the sensitivity calculation.
  • domain assumption Neutron-nucleus cross section is approximately flat in nuclear mass while dark-matter--nucleus cross section scales as A²
    Used to argue that the combined observable suppresses neutron background by more than an order of magnitude.

pith-pipeline@v0.9.0 · 5510 in / 1455 out tokens · 34769 ms · 2026-05-14T17:54:03.389554+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages · 5 internal anchors

  1. [2]

    Fossil tracks of charged particles in mica and the age of minerals,

    P. B. Price and R. M. Walker, “Fossil tracks of charged particles in mica and the age of minerals,” Journal of Geophysical Research68, 4847–4862 (1963)

  2. [3]

    G. A. Wagner and P. Van den haute,Fission Track Dating(Ferdinand Enke Verlag, Stuttgart, 1992)

  3. [4]

    Fission track dating of phosphate minerals and the thermochronology of apatite,

    A. J. W. Gleadow, D. X. Belton, B. P. Kohn, and R. W. Brown, “Fission track dating of phosphate minerals and the thermochronology of apatite,” Reviews in Mineralogy and Geochemistry48, 579–630 (2002)

  4. [5]

    Search for magnetic monopoles in deep ocean deposits,

    R. L. Fleischer, Jr. Hart, H. R., I. S. Jacobs, P. B. Price, W. M. Schwarz, and F. Aumento, “Search for magnetic monopoles in deep ocean deposits,” Physical Review184, 1393–1397 (1969)

  5. [6]

    Search for supermassive magnetic monopoles using mica crystals,

    P. B. Price and M. H. Salamon, “Search for supermassive magnetic monopoles using mica crystals,” Physical Review Letters56, 1226–1229 (1986)

  6. [7]

    Compositional dependence of the formation of nuclear tracks in muscovite mica,

    D. P. Snowden-Ifft, M. K. Y. Chan, and R. Frenkel, “Compositional dependence of the formation of nuclear tracks in muscovite mica,” Phys. Rev. Lett.74, 4133–4136 (1995)

  7. [8]

    Paleo-detectors: Searching for Dark Matter with Ancient Minerals

    Andrzej K. Drukier, Sebastian Baum, Katherine Freese, Maciej G´ orski, and Patrick Stengel, “Paleo-detectors: Searching for Dark Matter with Ancient Minerals,” Phys. Rev. D99, 043014 (2019), arXiv:1811.06844 [astro-ph.CO]

  8. [9]

    Searching for Dark Matter with Paleo-Detectors,

    Sebastian Baum, Andrzej K. Drukier, Katherine Freese, Maciej G´ orski, and Patrick Stengel, “Searching for Dark Matter with Paleo-Detectors,” Phys. Lett. B803, 135325 (2020), arXiv:1806.05991 [astro-ph.CO]

  9. [10]

    Mineral detection of neutrinos and dark matter. A whitepaper,

    Sebastian Baumet al., “Mineral detection of neutrinos and dark matter. A whitepaper,” Phys. Dark Univ.41, 101245 (2023), arXiv:2301.07118 [astro-ph.IM]

  10. [11]

    DARWIN: towards the ultimate dark matter detector,

    J. Aalberset al.(DARWIN), “DARWIN: towards the ultimate dark matter detector,” JCAP11, 017 (2016)

  11. [12]

    The XLZD Design Book: Towards the Next-Generation Dark Matter Observatory,

    J. Aalberset al.(XLZD), “The XLZD Design Book: Towards the Next-Generation Dark Matter Observatory,” (2024)

  12. [13]

    Measuring Changes in the Atmospheric Neutrino Rate Over Gigayear Timescales,

    Johnathon R. Jordan, Sebastian Baum, Patrick Stengel, Alfredo Ferrari, Maria Cristina Morone, Paola Sala, and Joshua Spitz, “Measuring Changes in the Atmospheric Neutrino Rate Over Gigayear Timescales,” Phys. Rev. Lett.125, 231802 (2020), arXiv:2004.08394 [hep-ph]

  13. [14]

    Measuring solar neutrinos over gigayear timescales with paleo detectors,

    Natalia Tapia-Arellano and Shunsaku Horiuchi, “Measuring solar neutrinos over gigayear timescales with paleo detectors,” Phys. Rev. D103, 123016 (2021), arXiv:2102.01755 [hep-ph]

  14. [15]

    Rocks, water, and noble liquids: Unfolding the flavor contents of supernova neutrinos,

    Sebastian Baum, Francesco Capozzi, and Shunsaku Horiuchi, “Rocks, water, and noble liquids: Unfolding the flavor contents of supernova neutrinos,” Phys. Rev. D106, 123008 (2022), arXiv:2203.12696 [hep-ph]

  15. [16]

    Sedimentary rocks from Mediterranean drought in the Messinian age as a probe of the past cosmic ray flux,

    Lorenzo Caccianiga, Lorenzo Apollonio, Federico Maria Mariani, Paolo Magnani, Claudio Galelli, and Alessandro Veutro, “Sedimentary rocks from Mediterranean drought in the Messinian age as a probe of the past cosmic ray flux,” Phys. Rev. D110, L121301 (2024), arXiv:2405.04908 [astro-ph.HE]

  16. [17]

    A volcanic chronosequence as a time-resolved paleo-detector array to study the cosmic-ray flux in the late Pleistocene and Holocene,

    Claudio Galelli, Lorenzo Caccianiga, Lorenzo Apollonio, Paolo Magnani, and Vincent Breton, “A volcanic chronosequence as a time-resolved paleo-detector array to study the cosmic-ray flux in the late Pleistocene and Holocene,” JCAP04, 023 (2026), arXiv:2510.23126 [astro-ph.HE]

  17. [18]

    New Projections for Dark Matter Searches with Paleo-Detectors,

    S. Baum, T. D. P. Edwards, K. Freese, and P. Stengel, “New Projections for Dark Matter Searches with Paleo-Detectors,” Instruments5, 21 (2021), arXiv:2106.06559 [astro-ph.CO]

  18. [19]

    Refining the sensitivity of new physics searches with ancient minerals,

    A. Fung, T. Lucas, L. Balogh, M. Leybourne, and A. C. Vincent, “Refining the sensitivity of new physics searches with ancient minerals,” Phys. Rev. D112, 043040 (2025), arXiv:2504.08885 [hep-ph]

  19. [20]

    Passive Low-Energy Nuclear-Recoil Detection with Color Centers,

    B. K. Cogswell, A. Goel, and P. Huber, “Passive Low-Energy Nuclear-Recoil Detection with Color Centers,” Phys. Rev. Applied16, 064060 (2021), arXiv:2104.13926 [physics.ins-det]

  20. [21]

    Nuclear recoil detection with color centers in bulk lithium fluoride,

    G. R. Araujoet al.(PALEOCCENE), “Nuclear recoil detection with color centers in bulk lithium fluoride,” (2025), arXiv:2503.20732 [physics.ins-det]

  21. [22]

    Benchtop mesoSPIM: a next-generation open-source light-sheet microscope for cleared samples,

    N. Vladimirovet al., “Benchtop mesoSPIM: a next-generation open-source light-sheet microscope for cleared samples,” Nature Communications (2024), 10.1038/s41467-024-46770-2

  22. [23]

    The CMS experiment at the CERN LHC,

    S. Chatrchyanet al.(CMS), “The CMS experiment at the CERN LHC,” JINST3, S08004 (2008)

  23. [24]

    The ATLAS Experiment at the CERN Large Hadron Collider,

    G. Aadet al.(ATLAS), “The ATLAS Experiment at the CERN Large Hadron Collider,” JINST3, S08003 (2008)

  24. [25]

    SRIM – The Stopping and Range of Ions in Matter (2010),

    J. F. Ziegler, M. D. Ziegler, and J. P. Biersack, “SRIM – The Stopping and Range of Ions in Matter (2010),” Nucl. Instrum. Meth. B268, 1818–1823 (2010). 13

  25. [26]

    On the use of SRIM for computing radiation damage exposure,

    S. Agarwal, Y. Lin, C. Li, R. E. Stoller, and S. J. Zinkle, “On the use of SRIM for computing radiation damage exposure,” Nucl. Instrum. Meth. B503, 11–29 (2021)

  26. [27]

    Implication of neutrino backgrounds on the reach of next generation dark matter direct detection experiments

    J. Billard, L. Strigari, and E. Figueroa-Feliciano, “Implication of neutrino backgrounds on the reach of next generation dark matter direct detection experiments,” Phys. Rev. D89, 023524 (2014), arXiv:1307.5458 [hep-ph]

  27. [28]

    O’Hare,New Definition of the Neutrino Floor for Direct Dark Matter Searches,Phys

    Ciaran A. J. O’Hare, “New definition of the neutrino floor for direct dark matter searches,” Phys. Rev. Lett.127, 251802 (2021), arXiv:2109.03116 [hep-ph]

  28. [29]

    Comparison of ion-beam irradiation effects in xsub 2yosub 4 compounds,

    L Wang, W Gong, S Wang, and R C Ewing, “Comparison of ion-beam irradiation effects in xsub 2yosub 4 compounds,” Journal of the American Ceramic Society82(1999)

  29. [30]

    Sputter- ing of grains in c-type shocks,

    P. W. May, G. Pineau des Forˆ ets, D. R. Flower, D. Field, N. L. Allan, and J. A. Purton, “Sputter- ing of grains in c-type shocks,” Monthly Notices of the Royal Astronomical Society318, 809–816 (2000), https://academic.oup.com/mnras/article-pdf/318/3/809/3474381/318-3-809.pdf

  30. [31]

    thesis, Institut d’Astrophysique Spatiale (2014)

    Marco Bocchio,Modelling Dust Processing and Evolution in Extreme Environments as seen by Herschel Space Observatory, Ph.D. thesis, Institut d’Astrophysique Spatiale (2014)

  31. [32]

    Olivine Mineral Data,

    Webmineral Mineralogy Database, “Olivine Mineral Data,”https://webmineral.com/data/Olivine.shtml, accessed: 2026-05-07

  32. [33]

    Optical properties of some F-aggregate centers in LiF,

    J. Nahum and D. A. Wiegand, “Optical properties of some F-aggregate centers in LiF,” Physical Review154, 817–830 (1967)

  33. [34]

    Optical bands of F 2 and F + 3 color centers in LiF,

    G. Baldacchini, M. Cremona, G. d’Auria, R. M. Montereali, and V. Kalinov, “Optical bands of F 2 and F + 3 color centers in LiF,” Physical Review B54, 17508–17516 (1996)

  34. [35]

    A. M. Zaitsev,Optical Properties of Diamond: A Data Handbook(Springer, Berlin, Heidelberg, 2001)

  35. [36]

    Silicon carbide color centers for quantum applications,

    S. Castelletto and A. Boretti, “Silicon carbide color centers for quantum applications,” Journal of Physics: Photonics2, 022001 (2020)

  36. [37]

    Inelastic and Elastic Scattering of 187-Mev Electrons from Selected Even-Even Nuclei,

    R. H. Helm, “Inelastic and Elastic Scattering of 187-Mev Electrons from Selected Even-Even Nuclei,” Phys. Rev.104, 1466–1475 (1956)

  37. [38]

    Review of mathematics, numerical factors, and corrections for dark matter experiments based on elastic nuclear recoil,

    J. D. Lewin and P. F. Smith, “Review of mathematics, numerical factors, and corrections for dark matter experiments based on elastic nuclear recoil,” Astropart. Phys.6, 87–112 (1996)

  38. [39]

    The RAVE Survey: Constraining the Local Galactic Escape Speed

    M. C. Smithet al., “The RAVE Survey: Constraining the Local Galactic Escape Speed,” Mon. Not. Roy. Astron. Soc. 379, 755–772 (2007), arXiv:astro-ph/0611671

  39. [40]

    The Local Dark Matter Density

    J. I. Read, “The Local Dark Matter Density,” J. Phys. G41, 063101 (2014), arXiv:1404.1938 [astro-ph.GA]

  40. [41]

    Dark Matter Search Results from 4.2 Tonne-Years of Exposure of the LUX-ZEPLIN (LZ) Experiment,

    J. Aalberset al.(LZ), “Dark Matter Search Results from 4.2 Tonne-Years of Exposure of the LUX-ZEPLIN (LZ) Experi- ment,” (2024), arXiv:2410.17036 [hep-ex]

  41. [42]

    The CMS Data Acquisition System for the Phase-2 Upgrade

    A. Zabiet al., “The CMS data acquisition system for the Phase-2 upgrade,” Proceedings of Science (Topical Workshop on Electronics for Particle Physics 2018) (2019), 10.22323/1.343.0129, arXiv:1806.08975

  42. [43]

    Imaging nuclear recoil damage in minerals with light-sheet microscopy,

    Samuel Hedges, “Imaging nuclear recoil damage in minerals with light-sheet microscopy,” Talk presented at the MDνDM 2026 – Mineral Detection of Neutrinos and Dark Matter workshop Karlsruhe, Germany (2026)

  43. [44]

    A petavoxel fragment of human cerebral cortex reconstructed at nanoscale resolution,

    A. Shapson-Coeet al., “A petavoxel fragment of human cerebral cortex reconstructed at nanoscale resolution,” Science 384, eadk4858 (2024). 14