[論文レビュー] Fully ab-initio all-electron calculation of dark matter--electron scattering in crystals with evaluation of systematic uncertainties

Cyrus E. Dreyer, Rouven Essig+3 more|arXiv (Cornell University)|Jun 26, 2023

Dark Matter and Cosmic Phenomena被引用数 7

ひとこと要約

この研究は、原子中心ガウス基底を用いた全電子計算から結晶中の DM–電子散乱を第一原理で算出し、QCDark を導入し、Si および Ge の結晶因子とレートに関する系統的不確かさを分析する。

ABSTRACT

We calculate target-material responses for dark matter--electron scattering at the extit{ab-initio} all-electron level using atom-centered gaussian basis sets. The all-electron effects enhance the material response at high momentum transfers from dark matter to electrons, $q\gtrsim \mathcal{O}\left({10\ αm_e} ight)$, compared to calculations using conventional plane wave methods, including those used in QEDark; this enhances the expected event rates at energy transfers $E \gtrsim 10$~eV, especially when scattering through heavy mediators. We carefully test a range of systematic uncertainties in the theory calculation, including those arising from the choice of basis set, exchange-correlation functional, number of unit cells in the Bloch sum, $\mathbf{k}$-mesh, and neglect of scatters with very high momentum transfers. We provide state-of-the-art crystal form factors, focusing on silicon and germanium. Our code and results are made publicly available as a new tool, called Quantum Chemistry Dark (``QCDark'').

研究の動機と目的

• 全電子ガウス基底集合を用いて結晶中の DM–電子散乱を第一原理で計算するためのフレームワークを開発する。
• 高運動量移行での全電子効果が結晶の form factor および散乱レートに与える影響を定量化する。
• 基底集合、交換相関汎関数、k点格子、単位胞数、q-max 切断による系統的不確かさを評価する。
• 公開可能なコード（QCDark）を提供し、既存アプローチとベンチマークを行う。

提案手法

• 原子中心ガウス基底集合を用いた Kohn–Sham DFT 内で DM–electron 散乱を定式化する。
• 占有状態と非占有 KS 状態間の行列要素から結晶 form factors f_crystal(q,E_e) を計算する（Eq. 9）。
• 擬ポテンシャルを用いず全電子効果を組み込み、重い媒介子に対して高-q 精度を実現する。
• 重い媒介/軽い媒介の極限と誘電 screenings を考慮して DM–electron 断面をモデル化する（Eq. 10 および Eq. 11）。
• Monkhorst–Pack グリッドで Brillouin zone を離散化し、Eq. 10 に基づくレートを計算し、q- および E_e-bin によるビニングを行う（Eq. 14）。
• Si および Ge の基底集合（TZP, def2-TZVP）と汎関数（PBE0 など）をベンチマークし、実験的バンドギャップに合わせてスコシア補正を適用する。

実験結果

リサーチクエスチョン

• RQ1 全電子・ガウス基底計算は、結晶中の DM–電子散乱 form factor を平面波アプローチと比較してどのように変えるか？
• RQ2 Si および Ge に対する first-principles の DM–電子散乱レート計算で支配的な系統的不確かさは何で、予測イベントレートにどのように影響するか？
• RQ3 基底集合の選択、交換相関汎関数、k点サンプリング、q-max 切断が計算された結晶 form factor およびレートにどう影響するか？
• RQ4 結晶中のダイ dielectric screening は特に重い媒介に対して DM–電子散乱にどのような影響を与えるか？
• RQ5 予測されるレートは、シリコンおよびゲルマニウム標的でのサブ GeV DM の潜在的検出感度にどう翻訳されるか？

主な発見

• 全電子効果は高い運動量移行 q>約10 alpha m_e で材料の応答を高め、E>約10 eV のイベントレートを増加させ、特に重い媒介で顕著。
• 結晶 form factor は、従来の平面波結果と比較して q>約4 alpha m_e で明確な寄与（コア部および高周波成分の明示的寄与）により拡張された。
• 基底集合、XC汎関数、単位胞数、k-メッシュ、最大動量移動の系統的不確かさを定量化。適切な選択（例：TZP/def2-TZVP と PBE0）で収束する。
• シリコンとゲルマニウムは最新の結晶 form factor を提供。Ge は約30 eV 周辺の半コア 3d 遷移に対する感度を示し、収束にはより密な k-grid が必要。
• QCDark コードは以前のアプローチ（例：EXCEED-DM PAW 再構成）を再現・拡張し、より広く公開されている使用を可能にしている。