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 parametrization. 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 a crucial and fundamental distribution of the quenching factor, which can be associated to the collisional interactions underlying the transport phenomena. The discrepancies observed with Lindhard’s model for the quenching factor at nuclear recoil energies below 4 keV are primarily attributed to lattice binding effects and the specific characteristics of the 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.

Close

PrevNext

View figure in article

• Dark matter
• Particle dark matter
• Molecular dynamics