New Windows on Heavy Dark Matter: Mineral Melt Modelling and X-Ray Readout for Muscovite Mica

Technical Summary: New Windows on Heavy Dark Matter: Mineral Melt Modelling and X-Ray Readout for Muscovite Mica

Problem Statement
Conventional direct detection experiments lose sensitivity to heavy composite dark matter (DM) candidates with masses significantly above a microgram ($\sim 10^{18}$ GeV) because the galactic flux through meter-scale detectors drops below one particle per year. While ancient minerals offer a "paleodetector" alternative with gigaton-year exposures, previous searches using muscovite mica relied on chemical etching and optical microscopy. These methods are inefficient for detecting large composites (radii from nanometers to microns) and lack systematic calibration of the damage mechanisms. Furthermore, prior constraints derived from recasting monopole searches in mica suffer from unverified assumptions regarding track retention (specifically for alpha-recoil damage) and sample selection biases that may have excluded high-damage events. This work addresses the need for a quantitative framework to model energy deposition by heavy composites, a calibrated readout method for micron-scale damage, and a robust geological validation protocol.

Methodology
The authors develop a multi-faceted approach combining theoretical modeling, numerical simulation, and experimental readout:

1. Theoretical Framework (Energy Deposition):

   • Two interaction models are considered: the Opaque (Geometric) Limit, where the composite reflects all incident nuclei within its cross-section, and the Diffuse (Constituent) Limit, where loosely bound constituents interact individually.
   • Thermal Spike Modeling: Using a Sedov-Taylor formalism, the authors derive analytic scalings for the melt track radius ($R_{melt}$) as a function of the composite radius ($R_D$). They treat energy deposition as an instantaneous adiabatic injection relative to thermal diffusion timescales.
   • SRIM/TRIM Calibration: Numerical simulations of nuclear recoil cascades are used to validate the analytic models at sub-micron scales and to calibrate the phonon efficiency factor ($\eta$), which determines the fraction of deposited energy available for local heating. The simulations yield $\eta \approx 0.75$.

2. Experimental Readout (XRF Transmission):

   • A novel, non-destructive readout method is demonstrated using X-ray Fluorescence (XRF) mapping with a copper backing contrast technique.
   • A thin copper sheet is placed beneath the cleaved mica. Intact mica attenuates X-rays, suppressing the Cu signal. Damage tracks (melt voids or bore holes) expose the copper, creating localized enhancements in Cu K$\alpha$ fluorescence.
   • Calibration is performed using laser-ablated melt regions (50 $\mu$m and 150 $\mu$m) to establish a minimum detectable feature radius ($R_{min} = 25$ $\mu$m) under standard operating conditions (Bruker M6 Jetstream).

3. Geological Validation:

   • The effective exposure time ($t_{exp}$) is constrained by a dual-age determination: primary crystallization age (via $^{87}$Rb/$^{86}$Sr or U/Pb geochronology) and track retention age (via in-situ $^{238}$U spontaneous fission track counting).
   • The authors argue that spontaneous fission tracks (retention temperature $\sim 325^\circ$C) serve as a conservative proxy for the retention of larger, hydrodynamic melt tracks produced by heavy DM, whereas alpha-recoil tracks (retention $\sim 30^\circ$C) are unreliable over gigayear timescales.

4. Sensitivity Projections:

   • The authors integrate the energy deposition models, XRF calibration, and geological constraints with the Standard Halo Model (SHM) to project 90% C.L. exclusion contours.
   • They account for overburden attenuation (energy loss as DM traverses the Earth) and identify a "boring regime" for large opaque composites that are decelerated below the melt threshold but retain enough energy to physically bore a cylindrical puncture through the lattice.

Key Contributions and Results

• Analytic and Numerical Models: The paper provides the first quantitative derivation of melt track radii for heavy composite DM in both opaque and diffuse regimes. It validates these models against SRIM simulations for $R_D < 1$ nm and establishes a calibrated phonon efficiency ($\eta \approx 0.75$).
• Readout Demonstration: The authors successfully demonstrate the copper-backed XRF method, identifying laser-ablated features as small as 50 $\mu$m in diameter (corresponding to a 25 $\mu$m radius) with a measured contrast of 14%. This establishes a practical detection threshold for micron-scale damage.
• New Detection Channels: The work identifies a "sub-melt" or "boring" detection mode for large opaque composites ($R_D \ge 25$ $\mu$m) that have been significantly slowed by the Earth's overburden. These composites create clean cylindrical voids rather than melt halos, which are detectable via XRF transmission.
• Re-evaluation of Prior Constraints: The authors critically revisit previous dark matter exclusions derived from the Price–Salamon monopole search. They identify that these limits are compromised by:
  1. The assumption that alpha-recoil tracks (which anneal on $\sim$Myr timescales) are reliable proxies for Gyr-scale DM tracks.
  2. The "optically clean" sample selection criterion, which likely rejected samples containing the macroscopic melt features the new method seeks to detect.
  • Consequently, the paper delineates a specific "recastable" band of parameter space (between the fission-track retention threshold and the onset of macroscopic melt features) where previous limits remain robust, and a complementary high-cross-section band where they do not.
• Projected Sensitivities: For a benchmark exposure of $1 \text{ m}^2 \times 10^9$ years, the paper presents projected sensitivity contours for both opaque and diffuse composite DM. These projections extend into regions of high mass and cross-section previously inaccessible to direct detection, particularly for large composites where the boring regime applies.

Significance
The paper claims to establish a new, robust framework for using muscovite mica as a paleodetector for heavy composite dark matter. Its significance lies in:

1. Extending Sensitivity: It opens a detection window for composites with radii from nanometers to microns, a regime where traditional etching methods are inefficient.
2. Methodological Rigor: It replaces unverified assumptions about track retention with a dual-geochronological validation strategy and provides the first systematic calibration of the damage mechanism (via SRIM and laser ablation).
3. Non-Destructive Readout: It introduces a rapid, large-area XRF readout method that avoids the chemical etching and potential sample destruction of previous searches.
4. Correcting Historical Limits: By identifying the shortcomings of prior monopole-based recasts, the paper clarifies the true exclusion power of existing mica data and defines the parameter space where new searches are necessary.

The authors conclude that scaling this methodology to square-meter exposures is within reach of existing instrumentation, offering a promising avenue to probe the fundamental nature of heavy composite dark matter candidates that have remained undetected by conventional means.