Accurate electron-recoil ionization factors for dark matter direct detection in xenon, krypton and argon

A. R. Caddell,V. V. Flambaum,B. M. Roberts

DOI: https://doi.org/10.1103/PhysRevD.108.083030

2023-05-09

Abstract:While most scintillation-based dark matter experiments search for Weakly Interacting Massive Particles (WIMPs), a sub-GeV WIMP-like particle may also be detectable in these experiments. While dark matter of this type and scale would not leave appreciable nuclear recoil signals, it may instead induce ionization of atomic electrons. Accurate modelling of the atomic wavefunctions is key to investigating this possibility, with incorrect treatment leading to a large suppression in the atomic excitation factors. We have calculated these atomic factors for argon, krypton and xenon and present the tabulated results for use with a range of dark matter models. This is made possible by the separability of the atomic and dark matter form factor, allowing the atomic factors to be calculated for general couplings; we include tables for vector, scalar, pseudovector, and pseudoscalar electron couplings. Additionally, we calculate electron impact total ionization cross sections for xenon using the tabulated results as a test of accuracy. Lastly, we provide an example calculation of the event rate for dark matter scattering on electrons in XENON1T and show that these calculations depend heavily on how the low-energy response of the detector is modelled.

High Energy Physics - Phenomenology,Cosmology and Nongalactic Astrophysics,Atomic Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to accurately calculate the ionization factors of atomic electrons for weakly interacting massive particles (WIMP - like particles) on the sub - GeV scale when directly detecting dark matter in xenon, krypton and argon. Such light dark matter does not produce significant nuclear recoil signals, but may cause ionization of atomic electrons. Therefore, accurately modeling the atomic wave function is the key, and incorrect treatment will lead to a significant reduction in the atomic excitation factor. The author solves this problem by calculating these atomic factors and providing the results in tabular form for use in combination with various dark matter models.

Specifically, the paper mainly focuses on the following points:
1. **Calculation of atomic ionization factors**: In order to study the possibility of interaction between sub - GeV WIMP - like particles and atomic electrons, it is necessary to accurately calculate the atomic ionization factors. These factors are calculated for different types of electron couplings (vector, scalar, pseudo - vector, pseudo - scalar) and are provided in tabular form.
2. **Calculation of total ionization cross - section for electron collisions**: As an accuracy test, the author also calculates the total ionization cross - section for electron collisions in xenon and compares it with experimental data.
3. **Calculation of event rates**: The author provides an example calculation, showing the event rate of dark matter - electron scattering in the XENON1T experiment, and emphasizes the influence of low - energy response modeling on the calculation results.

### Main contributions
- **Accurate calculation of atomic ionization factors**: The author uses the relativistic Hartree - Fock method, taking into account important many - body effects, ensuring accuracy in different energy deposition and momentum transfer ranges.
- **Verification of total ionization cross - section for electron collisions**: The accuracy of the calculation method is verified by comparison with high - precision experimental data.
- **Example of event rate calculation**: A specific calculation example is provided, showing how to combine atomic ionization factors and dark matter models to calculate event rates, taking into account the response characteristics of the detector.

### Formula summary
- **Effective Yukawa potential**:
\[
V=\hbar c\alpha_{\chi}\frac{e^{-\mu r}}{r}
\]
where \(\alpha_{\chi}\) is the dark matter - electron coupling strength, and \(\mu\) is the reciprocal of the interaction length scale.

- **Differential scattering cross - section**:
\[
\frac{d\sigma_{n\kappa}}{dE}=\frac{8\pi\alpha_{\chi}^2}{cv^2}\int_{q_-}^{q_+}\frac{q dq}{(q^2+\mu^2)^2}K_{n\kappa}(E,q)\frac{E}{E_H}
\]
where \(E\) is the energy deposition, \(v\) is the fixed dark matter velocity, \(q\) is the momentum transfer, \(E_H = m_e c^2\alpha^2\approx27.2\ eV\) is the Hartree energy unit, and \(K\) is the dimensionless atomic excitation factor.

- **Atomic excitation factor**:
\[
K_{n\kappa}\equiv E_H\sum_m\sum_{\kappa'm'}\left|\langle\epsilon\kappa'm'|e^{i\mathbf{q}\cdot\mathbf{r}}|n\kappa m\rangle\right|^2
\]

- **Velocity - averaged differential scattering cross - section**:
\[
\left\langle\frac{d\sigma}{dE}\right\rangle=\frac{\bar{\sigma}_e c^2}{2m_e c^2}\int dv\frac{f(v)}{v/c}\int_{q_-}^{q_+}\frac{a_0^2q dq}{q^3}|F_{\chi}(q)|^2K(E,q)
\]
where \(\bar{\sigma}\)