[论文解读] 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的三维速度分布，用于定向暗物质探测，重建了径向速度分布，并揭示了角向与径向分量中的昼夜与年周期调制。关键结果表明，当事件统计量足够时，径向速度分布中的年周期调制可被探测到；简化后的平移麦克斯韦分布能很好地近似精确形式，从而为未来实验提供高效的分析手段。

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）的三维速度分布建立模型，以支持定向暗物质探测。
• 研究WIMP速度分布的角向与径向分量中的昼夜与年周期调制。
• 开发并验证贝叶斯拟合技术，用于从模拟数据中重建径向WIMP速度分布。
• 评估事件统计量与速度分布模型选择对径向分布重建精度的影响。
• 为多个主要地下探测站点提供实验室参考系中的角向速度模式参考。

提出的方法

• 在银河系、黄道、赤道、地球、水平及实验室坐标系中，使用蒙特卡洛方法模拟3D 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 $ 的中位数及1σ/2σ置信区间。
• 通过比较在银河系逃逸速度 $ v_{\rm esc} = 500 $ 与 $ 600 $ km/s 条件下精确形式与简化形式的速度分布，验证模型准确性。

实验结果

研究问题

• RQ1 昼夜与年周期调制如何影响不同参考系中3D WIMP速度的角向分布？
• RQ2 贝叶斯重建能否准确恢复来自模拟定向数据的径向WIMP速度分布？
• RQ3 速度分布模型选择（如 $ f_{1,{\rm sh},\Delta v} $）如何影响 $ v_0 $ 与 $ \Delta v $ 参数估计的精度？
• RQ4 事件数（每60天50个 vs. 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 时，精确形式与简化形式之间的形状差异消失，验证了在数据分析中使用简化模型的合理性。
• 在不同季节与实验室位置下，重建的速度分布与参数不确定性保持一致，证实了贝叶斯方法的鲁棒性。