Molecular Dynamics Simulations on Nuclear Recoils in Silicon Crystals towards Single Electron-Hole Pair Ionization Yields

Abstract

We have developed a novel methodology utilizing molecular dynamics simulations to evaluate the ionization yields of nuclear recoils in crystalline silicon. This approach enables analytical exploration of atomic-scale transport within the lattice without necessitating parameterization. The quenching factors across the nuclear recoil energy range from 20 eV to 10 keV have been thoroughly investigated. A remarkable agreement with experimental data is achieved, particularly for the minimal energy regime conducted to date, reaching the level of a single electron-hole pair. This work presents evidence of a crucial and fundamental distribution of the quenching factor, which can be associated to the collisional interactions underlying the transport phenomena. The region below 4 keV of the quenching factor, where discrepancies have been observed with the Lindhard's model, is found to be significantly attributed to the lattice binding effect and the specific crystal structure. In contrast, a gradual functional relationship is identified below approximately 100 eV, indicating that the quenching factor is influenced by the crystallographic orientation of the target material. From a distributional perspective, our analysis allows for the determination of the minimum exclusion mass for the dark matter nucleon elastic scattering channel at 0.29 $\mathrm{GeV}/c^2$, thereby significantly enhancing sensitivity for the sub-$\mathrm{GeV}/c^2$ mass region.

7 pages, 5 figures, revised manuscript

References in corpus (27)

• Observation of Coherent Elastic Neutrino-Nucleus Scattering
• First results from the CRESST-III low-mass dark matter program
• Compatibility of DAMA/LIBRA dark matter detection with other searches
• Single-electron and single-photon sensitivity with a silicon Skipper CCD
• Limits on Light Weakly Interacting Massive Particles from the First 102.8 kg $\times$ day Data of the CDEX-10 Experiment
• Light Dark Matter Search with a High-Resolution Athermal Phonon Detector Operated Above Ground
• Search for Low-Mass Dark Matter with CDMSlite Using a Profile Likelihood Fit
• Measurement of the Coherent Elastic Neutrino-Nucleus Scattering Cross Section on CsI by COHERENT
• Results on Low-Mass Weakly Interacting Massive Particles from an 11 kg-day Target Exposure of DAMIC at SNOLAB
• SENSEI: First Direct-Detection Results on sub-GeV Dark Matter from SENSEI at SNOLAB
• The electron elevator: excitations across the band gap via a dynamical gap state?
• Searches for electron interactions induced by new physics in the EDELWEISS-III germanium bolometers
• Measurement of the ionization produced by sub-keV silicon nuclear recoils in a CCD dark matter detector
• Coherency in Neutrino-Nucleus Elastic Scattering
• Thermal detection of single e-h pairs in a biased silicon crystal detector
• Atomic limits in the search for galactic dark matter
• Results on sub-GeV Dark Matter from a 10 eV Threshold CRESST-III Silicon Detector
• Trajectory-dependent electronic excitations by light and heavy ions around and below the Bohr velocity
• Study of the ionization efficiency for nuclear recoils in pure crystals
• Investigating the sources of low-energy events in a SuperCDMS-HVeV detector
• A measurement of the ionization efficiency of nuclear recoils in silicon
• First measurement of the nuclear-recoil ionization yield in silicon at 100 eV
• Nuclear-recoil energy scale in CDMS II silicon dark-matter detectors
• Confirmation of the spectral excess in DAMIC at SNOLAB with skipper CCDs
• Ionization efficiency for nuclear recoils in silicon from $\sim 50$ eV to $3$ MeV
• Detection and Calibration of Low-Energy Nuclear Recoils for Dark Matter and Neutrino Scattering Experiments
• Intrinsic Fano factor of nuclear recoils for dark matter searches