---
title: Heavy-Element Paleodetectors in Ancient Minerals
url: https://www.emergentmind.com/topics/heavy-element-paleodetectors
type: topic
---

# Heavy-Element Paleodetectors in Ancient Minerals

Heavy-element paleodetectors are paleo-detectors implemented in ancient minerals whose target nuclei are sufficiently heavy to enhance coherent scattering and, in inelastic scenarios, to overcome kinematic thresholds for scattering. In the broad paleodetector program, insulators or poor semiconductors preserve latent damage tracks from nuclear recoils over geological timescales; heavy-element variants emphasize minerals containing high-mass nuclei such as lead, while related work on ultra-basic rocks such as olivine establishes the microphysics of track formation and readout in candidate solids [2606.05299] [1811.06844]. The resulting detectors are time-integrating solid-state track detectors rather than real-time instruments, with sensitivity controlled by recoil kinematics, stopping powers, mineral radiopurity, annealing stability, and nanoscale three-dimensional microscopy [2106.06559].

## 1. Concept and detection principle

Heavy-element paleodetectors inherit the general paleo-detector concept: a recoiling nucleus produced by dark matter or neutrino scattering deposits energy in a solid and creates a permanent damage track of broken bonds, provided the local energy deposition exceeds a material-dependent threshold and the mineral’s annealing time far exceeds its age [1811.06844]. In this framework, the detector is the mineral itself, and the measurement is performed by post hoc imaging of the accumulated damage record rather than by instrumenting a large target mass in real time.

The basic attraction of the method is the combination of long integration times and nanometric readout. Paleo-detectors use small samples of naturally occurring rocks that have been deep underground, typically for $\mathcal{O}(1)$ Gyr, and modern microscopy techniques promise nanometer-resolution readout in macroscopic samples [2106.06559]. In heavy-element proposals, integration over geological times $T\sim10^6$–$10^9$ yr builds up $O(10^3$–$10^6)$ t·yr exposure in even $\sim100$ g samples, while conventional forecasts for the broader paleo-detector program reach keV recoil thresholds and 100 kilotonne-yr exposures [2606.05299] [2106.06559].

The preference for heavy nuclei follows directly from the scattering kinematics and rates. For elastic spin-independent scattering, the sensitivity scales roughly as $A_T^2$, so large mass number nuclei are preferred [1811.06844]. For inelastic dark matter, the dependence is stronger still: the maximum accessible mass splitting is
$$
\delta_{\max,A}=\tfrac12\,\mu_A\,v_{\max}^2,
$$
so heavy nuclei increase the kinematic reach because $\mu_A$ grows with target mass [2606.05299].

## 2. Stopping powers, latent-track formation, and morphology

Track formation in insulating crystals is governed by the competition between electronic stopping and nuclear stopping. An energetic ion loses energy through ionization and excitation of target electrons,
$$
S_e=\frac{dE}{dx}_e,
$$
and through elastic collisions with target nuclei,
$$
S_n=\frac{dE}{dx}_n.
$$
A useful diagnostic is the fractional nuclear stopping
$$
f_n(E)\equiv \frac{S_n}{S_e+S_n},
$$
which parameterizes the transition from electronically dominated energy loss to collision-dominated cascade damage [2604.09732].

In the electronic, thermal-spike regime, track radius is commonly parameterized by the Szenes scaling law,
$$
R^2\propto
\begin{cases}
\ln(S_e/S_{e,\mathrm{th}}), & 0\le \ln(S_e/S_{e,\mathrm{th}})\le 1,\\
S_e/S_{e,\mathrm{th}}, & \ln(S_e/S_{e,\mathrm{th}})\ge 1,
\end{cases}
$$
where $S_{e,\mathrm{th}}$ is the electronic stopping threshold for continuous track formation [2604.09732]. In the nuclear-dominated regime, by contrast, damage is described as discontinuous “islands” of amorphous or vacancy-rich zones rather than a smooth cylinder, and empirically tracks become spotty when $S_n/S_{\rm total}\gtrsim0.5$ [2604.09732].

The clearest experimental benchmark for this morphology in a candidate paleo-detector mineral is the transmission electron microscopy study of Au-irradiated olivine. In that work, natural olivine Mg$_{1.8}$Fe$_{0.2}$SiO$_4$ was irradiated with 15 MeV Au$^{+5}$ and examined by focused ion-beam sectioning and scanning transmission electron microscopy at target depths from 423 to 2920 nm, without etching [2604.09732]. The measured widths changed only modestly across the energy range studied: approximate weighted means for linear tracks were 6.0, 6.1, 6.3, 6.8, 7.1, 7.4, 7.9, 7.7, and 7.0 nm from 423 to 2920 nm, while circular-track estimates ranged from 6.5 to 9.2 nm over the same depths [2604.09732]. The paper reports a significant change in track continuity across 0.4–12.9 MeV, with smooth, continuous tracks at shallow depths and spotty, discontinuous tracks at larger depths, consistent with the transition from electronic to nuclear stopping dominance predicted by SRIM [2604.09732].

This result is significant for heavy-element paleodetectors for two reasons. First, it demonstrates that latent tracks in a paleo-detector candidate can be resolved directly at the few-nanometer scale without etching. Second, it anchors the broader claim that track morphology is not determined by recoil energy alone but also by how that energy partitions between $S_e$ and $S_n$. A plausible implication is that heavy-element target optimization must account for morphology changes, not only total scattering rate.

## 3. Minerals, radiopurity, and geological provenance

Material choice is constrained simultaneously by scattering physics, track retention, and backgrounds. The principal criteria stated in the paleo-detector literature are large mass number for coherence enhancement, electrical resistivity $\gtrsim 2\times10^3\,\Omega\,\mathrm{cm}$ to allow track recording, high melting or annealing temperature for long track retention, ultra-low concentrations of radioactive contaminants such as $^{238}$U and $^{232}$Th, and in some cases the presence of hydrogen for neutron moderation [1811.06844].

| Geological setting or class | Examples | Noted properties |
|---|---|---|
| Marine evaporites | Halite, gypsum, epsomite | Typically $C^{238}\sim10^{-11}$ in weight; some contain H |
| Ultra-basic rocks | Olivine, phlogopite, nickelbischofite | Typically $C^{238}\sim10^{-10}$ in weight |
| Brine precipitates from deep geothermal aquifers | Laurionite, cottunite, paralaurionite, galena | Deep $(>4\ \mathrm{km})$ sources; Pb-bearing; extremely low U is possible |

The explicit heavy-element proposal centers on ancient, radiopure minerals containing lead. Brine precipitates from deep geothermal aquifers are identified as a possible source because hot $(>100\,^\circ\mathrm{C})$, reducing conditions render UO$_2$ and USiO$_4$ insoluble, leading to measured uranium concentrations $C^{238}_{\rm BP}\lesssim3\times10^{-12}$ g U/g in Gulf-Coast aquifers and a theoretical minimum $\lesssim1\times10^{-13}$ g/g [2606.05299]. These aquifers lie at 4–5 km depth, where cosmogenic neutrons are negligible relative to radiogenic ones, and their ages can range from Oligocene values of 23–34 Myr to older horizons up to $\sim1$ Gyr [2606.05299].

Laurionite is the benchmark Pb-bearing mineral in this literature. Its chemical formula is PbClOH, its Pb mass fraction is $w_{\rm Pb}=0.80$, its electrical resistivity is $\rho_r\sim10^8$–$10^{10}\ \Omega\cdot\mathrm{cm}$, and its hydrogen content helps moderate fast-neutron backgrounds [2606.05299]. Other Pb-bearing minerals explicitly noted are cottunite PbCl$_2$, paralaurionite PbOHCl, and galena PbS [2606.05299]. By contrast, olivine appears as an ultra-basic rock target with Fe and Mg, important both as a general paleo-detector material and as the mineral in which track-formation microphysics has been directly characterized [1811.06844] [2604.09732].

## 4. Readout modalities and the observable track spectrum

The observable in paleodetection is usually the damage-track length distribution rather than the recoil-energy spectrum itself. For a recoiling nucleus of energy $E_R$, the track length is approximated by its range,
$$
x(E_R)=\int_0^{E_R}\frac{dE'}{|dE/dx(E')|},
$$
with stopping powers obtained in practice from SRIM or TRIM [2106.06559]. The differential track-length spectrum is then written as
$$
\frac{dR}{dx}(x)=\frac{dR}{dE_R}\bigl(E_R(x)\bigr)\left|\frac{dE}{dx}\bigl(E_R(x)\bigr)\right|,
$$
or, more generally, as a convolution over the conditional distribution of track lengths at fixed recoil energy when stochastic cascade effects are included [2504.08885].

Two readout architectures dominate the literature. High-resolution Helium-Ion-Beam Microscopy is associated with $\sigma_L=1$ nm and $M=10$ mg, corresponding to an effective threshold $L_{\min}\simeq0.5$ nm, approximately 1 keV in heavy-element minerals [2106.06559]. Small-Angle X-ray Scattering tomography is associated with $\sigma_L=15$ nm and $M=100$ g, corresponding to $L_{\min}\simeq7.5$ nm, approximately a few keV [2106.06559]. In the heavy-element Higgsino study, the corresponding processable volumes are given as $\sim6$ mm$^3$ for HIBM and $\sim60$ cm$^3$ for SAXS [2606.05299].

Microscopy-based microphysical validation has also been demonstrated directly in olivine. The 2026 STEM study used a Thermo-Fisher Spectra 300 at 300 keV, bright-field imaging for track measurements, pixel sizes of 0.25 nm or 0.36 nm, and focused ion-beam staircase lift-outs at 10 depths [2604.09732]. Track extraction employed ridge detection plus transverse Gaussian fits for oblique tracks and ParticleAnalyzer plus circular-Gaussian fits for near-normal tracks [2604.09732]. This establishes that few-nanometer tracks can be measured without etching, but the scalable readout challenge remains the transition from $\mu$m$^3$-scale microscopy to macroscopic sample volumes.

A central revision introduced by TRIM-based sensitivity studies is that recoil energy and track length are not in one-to-one correspondence. The full distribution is
$$
\frac{dR}{dx}(x)=\sum_i \int_0^\infty dE_R\,\frac{dR_i}{dE_R}\,P_{i,\mathrm{track}}(E_R)\,P_i(x|E_R),
$$
where $P_i(x|E_R)$ is the conditional track-length distribution and $P_{i,\mathrm{track}}(E_R)$ is the probability that a recoil produces any visible track at all [2504.08885]. This refinement becomes critical near threshold.

## 5. Backgrounds, statistical inference, and revised sensitivities

Background control in paleodetectors is a mineralogical and geological problem as much as a detector-physics problem. Cosmogenic backgrounds are negligible if samples spend $\gtrsim5$ km-water-equivalent underground, or at depths $\gtrsim5$ km rock, where the cosmogenic neutron flux is negligible relative to the relevant radiogenic and neutrino backgrounds [2106.06559] [1811.06844]. Radiogenic backgrounds are dominated by the $^{238}$U chain, spontaneous fission, and $(\alpha,n)$ neutrons; one particularly important feature is the isolated U-234 recoil from the first $\alpha$-decay of U-238, which produces a monoenergetic 72 keV heavy recoil and a characteristic track-length signature near 200 nm in analyses that do not observe low-$Z$ tracks [1811.06844].

Hydrogen-bearing minerals are repeatedly emphasized because they moderate radiogenic neutrons efficiently. This is why gypsum and sinjarite are favorable in the general dark-matter forecasts, and why laurionite is attractive among Pb-bearing candidates [2106.06559] [2606.05299]. Fiducial uranium concentrations are $C_{^{238}}\simeq10^{-11}$ g/g for marine evaporites and $C_{^{238}}\simeq10^{-10}$ g/g for ultra-basic rocks, while the heavy-element Higgsino proposal targets $C^{238}\lesssim3\times10^{-12}$ g/g for ultimate reach [2106.06559] [2606.05299].

The modern sensitivity forecasts rely on binned spectral analyses with nuisance parameters for sample age, neutrino flux normalizations, and radioactive contamination. One formulation writes the expected counts in track-length bin $i$ as
$$
N_i=M\,t_{\rm age}\int dL'\,W(L';L_i^{\min},L_i^{\max})\frac{dR}{dL'},
$$
with a Gaussian window function of width $\sigma_L$, and uses a Poisson likelihood with Gaussian priors on nuisance parameters [2106.06559]. In the 2021 projections, the 90% CL exclusion criterion follows the profile-likelihood ratio with Asimov value $q=2.71$ [2106.06559]. Earlier work also formulated sensitivity through a signal-to-noise requirement
$$
\mathrm{SNR}=\frac{S}{\sqrt{\sum_i(B_i+\Sigma_i^2 B_i^2)}}\ge 3,\qquad S\ge 5,
$$
with $\Sigma_\nu=100\%$ for neutrinos and $\Sigma_{\rm radio}=1\%$ for radioactive backgrounds [1811.06844].

For elastic spin-independent dark matter, the updated 2021 projections give the following representative numbers. In a high-resolution scenario with 10 mg, $\sigma_L=1$ nm, and $t_{\rm age}=1$ Gyr, paleo-detectors can probe down to the conventional neutrino floor in a Xe-based direct detection experiment for masses $\lesssim10$ GeV/$c^2$; for $m_\chi=1$ GeV/$c^{-2}$ the sensitivity reaches $\sigma_p^{\rm SI}\sim10^{-43}$ cm$^2$, and for $m_\chi=5$ GeV/$c^{-2}$ it reaches $\sigma_p\sim10^{-44}$ cm$^2$ [2106.06559]. In a high-exposure scenario with 100 g and $\sigma_L=15$ nm, the projected sensitivity is $\sigma_p\sim10^{-46}$ cm$^2$ at $m_\chi=50$ GeV/$c^{-2}$ and $\sim3\times10^{-47}$ cm$^2$ at 100 GeV/$c^{-2}$ [2106.06559].

A major correction to earlier optimism comes from realistic TRIM-based track modeling. In olivine, below $E_R\sim1$ keV the track-formation probability falls below 50%; at $E_R=0.6$ keV only $\sim20$–30% of recoils yield any visible track, while by $E_R\sim10$ keV the yield is approximately 100% [2504.08885]. At the opposite end, once electronic stopping dominates above $E_R\sim50$–100 keV, the average track length saturates at $L_{\max}\sim200$ nm, producing a “track length barrier” [2504.08885]. The paper concludes that previous studies overestimated the number of tracks caused by weakly interacting particles, and that the realistic limits degrade by 0.5–1 decade in $\sigma_{\rm SI}$ at low $m_\chi$ and by $O(1)$ in $m_\chi$ relative to one-to-one energy–length treatments [2504.08885]. For heavy-element minerals, the same work notes that because electronic stopping scales roughly as $Z^2$, heavy-$Z$ targets exhibit this plateau at even shorter lengths and lower recoil energies [2504.08885].

## 6. Heavy nuclei for inelastic dark matter and halo-history sensitivity

The most distinctive role of heavy-element paleodetectors emerges in inelastic dark matter models such as the Higgsino, where existing paleo-detector targets lack sufficiently heavy nuclei to overcome the kinematic threshold for scattering [2606.05299]. In this scenario, dark matter $\chi$ up-scatters to a heavier state $\chi'$ with mass splitting $\delta$, and the minimum velocity for recoil energy $E_R$ is
$$
v_{\min}(E_R)=\frac{1}{\sqrt{2m_AE_R}}\left(\frac{m_A}{\mu_A}E_R+\delta\right).
$$
No scattering is kinematically allowed if $\delta>\delta_{\max,A}$, so target mass becomes the decisive figure of merit [2606.05299].

This is why Pb-bearing minerals are central. Since $\mu_A\simeq m_A$ for $m_\chi\gg m_A$, the ordering $\delta_{\max,\mathrm{Pb}}\gg\delta_{\max,\mathrm{Xe}}\gg\delta_{\max,\mathrm{Fe}}$ follows directly, and Pb with $A=207$ materially extends the accessible inelastic frontier [2606.05299]. The differential recoil rate used in the Higgsino benchmark is
$$
\frac{dR}{dE_R}(E_R)=\frac{\rho_\chi}{m_\chi}\,\frac{m_A}{2\mu_n^2}\,\sigma_{n\chi}\,A^2\,F_A^2(E_R)\,\eta(v_{\min}(E_R)),
$$
with benchmark parameters $\rho_\chi\approx0.3$ GeV/cm$^3$, $m_\chi\approx1.1$ TeV, and $\sigma_{n\chi}\simeq10^{-39}$ cm$^2$ [2606.05299].

The benchmark heavy-element forecast uses laurionite with $V=60$ cm$^3$, age $T=1$ Gyr, SAXS readout at $\sigma_x=15$ nm, $C^{238}=3\times10^{-12}$ g/g, and depth $\ge5$ km [2606.05299]. In the Standard Halo Model this yields sensitivity up to $\delta\approx560$ keV at $\sigma_{n\chi}=10^{-39}$ cm$^2$ [2606.05299]. When a Large Magellanic Cloud-induced high-velocity tail is included, the reach extends to $\delta\approx880$ keV for the present-day velocity distribution and to $\delta\approx960$ keV at pericenter; with the theoretical minimum uranium concentration, the ultimate reach is stated as $\delta\lesssim920$ keV [2606.05299].

The same study argues that heavy-element paleodetectors are uniquely sensitive to the history of the dark matter high-velocity tail. Simulations of Auriga halo 13 indicate that the Large Magellanic Cloud pericenter about 50 Myr ago injected a fast, unbound dark matter population, increasing the maximum Earth-frame speed from the Standard Halo Model value to present-day $v_{\max}\approx950$ km/s and pericenter $v_{\max}\approx990$ km/s [2606.05299]. For $\delta\gtrsim550$ keV, upscattering is kinematically forbidden until the LMC boost, so a Gyr-old sample would accumulate approximately 950 Myr of background before any signal; a younger $T\approx50$ Myr sample therefore optimizes signal-to-noise for large $\delta$ [2606.05299]. This temporal selectivity is a distinctive feature of time-integrating detectors and has no direct analogue in conventional real-time experiments.

Relaxed scenarios are also quantified. Even radio-impure samples with $C^{238}$ up to $3\times10^{-7}$ g/g from depths of only 2 km still probe $\delta\lesssim500$ keV, and a $10^6\times$ smaller volume, $V\sim60$ mm$^3$, remains competitive for $\delta\lesssim300$ keV [2606.05299]. This does not eliminate the challenges of radiopurity and depth, but it indicates that the unusually large Higgsino-nucleon cross section partially relaxes requirements that are otherwise stringent in the paleodetector program.

Taken together, the literature defines heavy-element paleodetectors as a specialized branch of paleo-detection in which mineralogy is tuned to high-mass nuclei, especially lead, while track-formation microphysics, background rejection, and nanoscale readout remain inherited from the broader ancient-mineral program. The central design tension is now clear: heavy nuclei increase coherent rate and inelastic reach, but stronger electronic stopping compresses the observable track spectrum and pushes the track-length barrier to shorter scales [2504.08885]. This suggests that the future development of heavy-element paleodetectors will depend on jointly optimizing target $Z$, hydrogen content, radiopurity, age, depth, and readout resolution rather than maximizing any single parameter alone.

Source: https://www.emergentmind.com/topics/heavy-element-paleodetectors