[论文解读] Studying dark matter with MadDM 3.1: a short user guide

Chiara Arina, Jan Heisig+5 more|arXiv (Cornell University)|Dec 16, 2020

Dark Matter and Cosmic Phenomena参考文献 27被引用 4

一句话总结

MadDM 3.1 是一个集成在 MadGraph5_aMC@NLO 中的自动化数值工具，可为通用新物理模型计算暗物质可观测量，如遗迹密度、直接探测反冲率和间接探测能谱。它支持使用 Pythia 8 或 PPPC4DMID 对 2→n 湮灭过程进行全微分截面计算，实现对对撞机、直接探测和间接探测实验中暗物质候选者的精确现象学研究。

ABSTRACT

MadDM is an automated numerical tool for the computation of dark-matter observables for generic new physics models. We announce version 3.1 and summarize its features. Notably, the code goes beyond the mere cross-section computation for direct and indirect detection. For instance, it allows the user to compute the fully differential nuclear recoil rates as well as the energy spectra of photons, neutrinos and charged cosmic rays for arbitrary $2 o n$ annihilation processes. This short user guide equips researchers with all the relevant information required to readily perform comprehensive phenomenological studies of particle dark-matter models.

研究动机与目标

• 提供一个全面的自动化工具，用于计算直接探测、间接探测和遗迹密度搜索中的暗物质可观测量。
• 超越截面计算，支持对核反冲率和光子、中微子及宇宙射线的能量谱进行全微分计算。
• 通过事件级模拟或预计算库，支持通用新物理模型中任意 2→n 湮灭过程。
• 整合来自 LUX、Xenon1T 和 Fermi-LAT 的最新实验约束，实现稳健的现象学解释。
• 通过在 MG5_aMC@NLO 框架内采用插件架构，简化研究人员的工作流程，支持模型无关的研究。

提出的方法

• 与 MadGraph5_aMC@NLO 集成，为 UFO 格式中的任意模型生成矩阵元。
• 使用 Pythia 8 在精确模式下进行部分子喷注和强子化，实现对 2→n 湮灭过程的完整模拟。
• 在“快速”模式下利用 PPPC4DMID 进行快速预计算的湮灭能谱，降低计算成本。
• 通过标准玻尔兹曼方程计算遗迹密度，利用湮灭截面因子 ξ² 进行重标度。
• 对矮椭圆星系的伽马射线流量应用 Fermi-LAT 似然分析，以计算排除极限和 p 值。
• 生成包含解析可观测量（如通量、截面和似然）的输出文件，适用于单点计算和扫描运行。

实验结果

研究问题

• RQ1 如何为多种目标材料的直接探测实验计算全微分核反冲率？
• RQ2 2→n 暗物质湮灭过程中，快中子、中微子和宇宙射线的能量谱是什么？
• RQ3 来自 LUX、Xenon1T 和 Fermi-LAT 的实验约束如何影响暗物质模型的可行性？
• RQ4 使用 Pythia 8 与 PPPC4DMID 对间接探测能谱计算的精度和效率有何影响？
• RQ5 如何在模型无关的方式下，使用 MadDM 3.1 计算遗迹密度并与观测值进行比较？

主要发现

• MadDM 3.1 可计算作为能量、散射角和时间函数的全微分核反冲率，支持详细的直接探测模拟。
• 该代码计算了所有湮灭通道的平均速度截面 ⟨σv⟩₀（单位：cm³s⁻¹），其值在 1.67×10⁻³⁸ 至 1.76×10⁻³⁰ cm³s⁻¹ 之间，具体取决于通道。
• 对于同一参数点，地球处的快中子通量计算为 1.57×10⁻¹⁴ 个/(cm² sr)，中微子通量在 8.89×10⁻¹⁸ 至 9.71×10⁻¹⁸ 个/(cm² sr) 之间。
• Fermi-LAT 约束给出的排除极限为 1.11×10⁻²⁵ 至 2.85×10⁻²⁵ cm³s⁻¹，总截面极限通过完整似然分析计算得出。
• 该程序能正确识别并跳过截面为零的过程（如 xrxr 和 y0y0），确保计算效率。
• 输出文件如 MadDM_results.txt 和 scan_run_01.txt 包含结构化、机器可读的数据，包括似然、p 值和光谱通量，便于后续分析。