Beyond-DFT $\textit{ab initio}$ Calculations for Accurate Prediction of Sub-GeV Dark Matter Experimental Reach

Elizabeth A. Peterson,Samuel L. Watkins,Christopher Lane,Jian-Xin Zhu

2023-09-30

Abstract:As the search space for light dark matter (DM) has shifted to sub-GeV DM candidate particles, increasing attention has turned to solid state detectors built from quantum materials. While traditional solid state detector targets (e.g. Si or Ge) have been utilized in searches for dark matter (DM) for decades, more complex, anisotropic materials with narrow band gaps are desirable for detecting sub-MeV dark matter through DM-electron scattering and absorption channels. In order to determine if a novel target material can expand the search space for light DM it is necessary to determine the projected reach of a dark matter search conducted with that material in the DM mass - DM-electron scattering cross-section parameter space. The DM-electron scattering rate can be calculated from first-principles with knowledge of the loss function, however the accuracy of these predictions is limited by the first-principles level of theory used to calculate the dielectric function. Here we perform a case study on silicon, a well-studied semiconducting material, to demonstrate that traditional Kohn-Sham density functional theory (DFT) calculations erroneously overestimate projected experimental reach. We show that for silicon this can be remedied by the incorporation of self-energy corrections as implemented in the GW approximation. Moreover, we emphasize the care that must taken in selecting the appropriate level of theory for predicting experimental reach of next-generation complex DM detector materials.

High Energy Physics - Phenomenology,Mesoscale and Nanoscale Physics,High Energy Physics - Experiment

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **Improve the prediction accuracy of sub - GeV DM (sub - GeV dark matter) detection experiments, especially when using new and complex materials as detectors**.

Specifically, as the focus of dark matter research gradually shifts to light dark matter candidate particles in the sub - GeV mass range, traditional semiconductor materials (such as Si and Ge) have limitations in detecting lighter dark matter. To expand the detection range, researchers begin to focus on quantum materials with characteristics such as narrow band gaps and anisotropy. These materials can detect sub - MeV - level dark matter through dark matter - electron scattering and absorption channels.

However, for these complex materials, traditional first - principles calculations based on density functional theory (DFT) may underestimate or overestimate the dark matter - electron scattering rate, thus leading to misjudgment of experimental sensitivity. Therefore, the focus of this paper is:

1. **Evaluate the limitations of traditional DFT methods in predicting dark matter - electron scattering**: Through a case study of silicon materials, the author shows that DFT may overestimate experimental sensitivity when calculating the dielectric function.
2. **Introduce a higher - precision calculation method**: In particular, the GW approximation (Green's function and Wannier functions) to correct the errors in DFT calculations, thereby more accurately predicting the experimental sensitivity of new detection materials.
3. **Emphasize the importance of choosing an appropriate theoretical level**: Ensure more reliable predictions for complex materials in the future and avoid being misled by low - precision calculations.

Through these improvements, researchers hope to provide more accurate theoretical support for future sub - GeV dark matter detection experiments and guide the selection and design of new materials.

### Summary of Key Formulas

- **Kohn - Sham Equation**:
\[
H_{\text{KS}}\psi_{\text{KS}}^{nk}(r)=\left[-\frac{1}{2}\nabla^{2}+v_{\text{ext}}(r)+\int\frac{\rho(r')}{|r - r'|}dr'+v_{\text{XC}}[\rho]\right]\psi_{\text{KS}}^{nk}(r)=\epsilon_{\text{KS}}^{nk}\psi_{\text{KS}}^{nk}(r)
\]

- **Correction Equation under the GW Approximation**:
\[
\left[-\frac{1}{2}\nabla^{2}+v_{\text{ext}}(r)+\int\frac{\rho(r')}{|r - r'|}dr'+\Sigma(\epsilon_{\text{GW}}^{nk})\right]\psi_{\text{GW}}^{nk}(r)=\epsilon_{\text{GW}}^{nk}\psi_{\text{GW}}^{nk}(r)
\]

- **Expression for Polarizability**:
\[
P_{GG'}(q,\omega)=\lim_{\eta\rightarrow0}\sum_{nn'k}\frac{\langle nk|e^{-i(q + G)\cdot r}|n'k + q\rangle\langle n'k + q|e^{i(q + G')\cdot r}|nk\rangle[f(\epsilon_{n'k + q})-f(\epsilon_{nk})]}{\epsilon_{n'k + q}-\epsilon_{nk}-\hbar\omega - i\hbar\eta}
\]

- **Scattering Rate Formula**:
\[
R(t)=\frac{\bar{\sigma}_e}{\rho_T\rho_\chi m_\chi\pi(2\pi)^{4}}\alpha\mu_{\chi e}^{2}\int d\omeg