[論文レビュー] Simulations of the 3-Dimensional Velocity Distribution of Halo Weakly Interacting Massive Particles for Directional Dark Matter Detection Experiments

Chung-Lin Shan|arXiv (Cornell University)|May 27, 2019

Dark Matter and Cosmic Phenomena参考文献 54被引用数 9

ひとこと要約

本稿は、方向性ダークマター検出を目的として、銀河ハロー内のWIMPの3次元速度分布をシミュレートし、ベイズ推論を用いて径方向速度分布を再構築することで、角度成分および径方向成分における昼夜および年間周期的変動を明らかにした。主な結果として、十分なイベント統計があれば、径方向速度分布における年間周期的変動が検出可能であり、簡略化されたシフトマクスウェル分布が正確な形に非常に近い近似を与えることから、将来的な実験における効率的なデータ解析が可能であることが示された。

ABSTRACT

In this paper, as a preparation of developing data analysis procedures for using 3-dimensional information offered by directional Dark Matter (DM) detection experiments, we study the patterns of the angular distribution of the Monte Carlo-generated 3-D velocity of halo Weakly Interacting Massive Particles (WIMPs) as well as apply the Bayesian fitting technique to reconstruct the radial distribution of the 3-D WIMP velocity. Besides the diurnal modulation of the angular WIMP velocity distribution, the so-called "directionality" of DM signals proposed in literature, we will also demonstrate possible "annual" modulations of both of the angular and the radial distributions of the 3-D WIMP velocity. Our Bayesian reconstruction results of (the annual modulation of) the radial WIMP velocity distribution will also be discussed in detail. For readers' reference, the angular distribution patterns of the 3-D WIMP velocity in the "laboratory (location)-dependent" reference frames of several underground laboratories are given in the Appendix.

研究の動機と目的

• 方向性ダークマター検出を目的として、銀河ハロー内における弱い力で相互作用する大質量粒子（WIMPs）の3次元速度分布をモデル化すること。
• WIMP速度分布の角度成分および径方向成分における昼夜および年間周期的変動を調査すること。
• シミュレートされたデータから径方向WIMP速度分布を再構築するためのベイズフィッティング技術の開発および検証すること。
• イベント統計および速度分布モデルの選択が、径方向分布再構築の精度に与える影響を評価すること。
• 主な地下検出施設における、実験装置依存のフレームにおける基準となる角度速度パターンを提供すること。

提案手法

• モンテカルロ法を用いて、銀河座標系、黄道座標系、赤道座標系、地球座標系、水平座標系、および実験室座標系の6つの座標系で3次元WIMP速度分布をシミュレートした。
• ベイズ推論を用いて径方向速度分布 $ f_{\rm r}(v) $ を再構築し、$ f_{1,{\rm sh},v_0}(v;v_0) $、$ f_{1,{\rm sh},v_{\rm e}}(v;v_{\rm e}) $、$ f_{1,{\rm sh},\Delta v}(v;v_0,\Delta v) $ といった複数のパラメトリック形式を用いた。
• 径方向速度の物理的モデルとしてシフトマクスウェル分布を採用し、$ f_{1,{\rm sh},{\rm vesc}}(v) $ および $ f_{1,{\rm sh}}(v) $ の解析的表現を導出した。
• 地球の公転運動をモデル化し、黄道座標系および赤道座標系において $ v_{\rm e}(t) $ を用いてWIMP速度分布に年間周期的変動を組み込んだ。
• 60日間の観測期間を対象に統計的フィッティングを実施し、$ v_0 $ および $ \Delta v $ の中央値および90%（1σ）・95%（2σ）信頼区間を推定した。
• 銀河脱出速度 $ v_{\rm esc} = 500 $ および $ 600 $ km/s の条件下で、正確な形と簡略化された形の速度分布を比較することで、モデルの正確性を検証した。

実験結果

リサーチクエスチョン

• RQ1 昼夜および年間周期的変動は、異なる参照フレームにおける3次元WIMP速度の角度分布にどのように影響を与えるか？
• RQ2 ベイズ再構築は、シミュレートされた方向性データから径方向WIMP速度分布を正確に回復できるか？
• RQ3 速度分布モデルの選択（例：$ f_{1,{\rm sh},\Delta v} $）が、$ v_0 $ および $ \Delta v $ のパラメータ推定の精度に与える影響は何か？
• RQ4 イベント数（60日間あたり50件対500件）の違いが、再構築された速度パラメータの統計的不確実性に与える影響は何か？
• RQ5 年間周期的変動および現実的な $ v_{\rm esc} $ 値下で、正確な形と簡略化された形のシフトマクスウェル速度分布はどのように比較されるか？

主な発見

• ベイズ再構築により、入力パラメータが良好に回復され、$ v_0 $ の中央値はイベント数およびモデルに応じて185–192 km/s、$ \Delta v $ は50–88 km/s程度に推定された。
• 60日間あたり500件のイベントの場合、$ v_0 $ の中央値は $ 185.0 \pm 9.0 $ km/s（1σ）、$ \Delta v $ は $ 83.7 \pm 17.0 $ km/s（1σ）と推定され、50件のケースに比べて精度が向上した。
• 径方向速度分布における年間周期的変動は検出可能であり、$ \Delta v $ は非対称な周期的変動を示す一方、$ v_0 $ は年平均で安定している。
• 簡略化されたシフトマクスウェル分布 $ f_{1,{\rm sh}}(v) $ は、正確な形 $ f_{1,{\rm sh},{\rm vesc}}(v) $ と極めて近似しており、特に $ v_{\rm esc} = 500 $ km/s 時に統計的不確実性と比較して無視できる差異しか示さない。
• $ v_{\rm esc} = 600 $ km/s 時に、正確な形と簡略化された形の間の形状差は消失し、簡略化モデルのデータ解析への適用が妥当であることが裏付けられた。
• 異なる季節および実験室位置においても、再構築された速度分布およびパラメータの不確実性が一貫しており、ベイズ手法の頑健性が確認された。