[论文解读] A forward-modelling method to infer the dark matter particle mass from strong gravitational lenses

Qiuhan He, Andrew Robertson+8 more|arXiv (Cornell University)|Oct 25, 2020

Galaxies: Formation, Evolution, Phenomena参考文献 109被引用 33

一句话总结

本文提出一种利用强引力透镜残差约束暗物质粒子质量的正向建模方法，特别针对暖暗物质（WDM）模型进行检验。通过模拟不同暗物质晕群的透镜图像，并利用近似贝叶斯计算比较残差的功率谱，该方法在暗物质晕质量截止为10⁷ M⊙的CDM类模型中，实现了𝑚DM > 4.10 keV的中位2σ约束；在真实𝑚DM = 2.22 keV的WDM宇宙中，成功恢复出𝑚DM ≈ 1.43–3.21 keV的估计值。

ABSTRACT

A fundamental prediction of the cold dark matter (CDM) model of structure formation is the existence of a vast population of dark matter haloes extending to subsolar masses. By contrast, other dark matter models, such as a warm thermal relic (WDM), predict a cutoff in the mass function at a mass which, for popular models, lies approximately between $10^7$ and $10^{10}~{ m M}_\odot$. We use mock observations to demonstrate the viability of a forward modelling approach to extract information about low-mass dark haloes lying along the line-of-sight to galaxy-galaxy strong lenses. This can be used to constrain the mass of a thermal relic dark matter particle, $m_\mathrm{DM}$. With 50 strong lenses at Hubble Space Telescope resolution and a maximum pixel signal-to-noise ratio of $\sim50$, the expected median 2$\sigma$ constraint for a CDM-like model (with a halo mass cutoff at $10^{7}~{ m M}_\odot$) is $m_\mathrm{DM} > 4.10 \, \mathrm{keV}$ (50% chance of constraining $m_{ m DM}$ to be better than 4.10 keV). If, however, the dark matter is a warm particle of $m_\mathrm{DM}=2.2 \, \mathrm{keV}$, our 'Approximate Bayesian Computation' method would result in a median estimate of $m_\mathrm{DM}$ between 1.43 and 3.21 keV. Our method can be extended to the large samples of strong lenses that will be observed by future telescopes, and could potentially rule out the standard CDM model of cosmogony. To aid future survey design, we quantify how these constraints will depend on data quality (spatial resolution and integration time) as well as on the lensing geometry (source and lens redshifts).

研究动机与目标

• 开发一种统计方法，从强引力透镜残差中提取关于低质量暗物质晕的信息。
• 检验图像残差的功率谱是否能够约束由于暖暗物质导致的晕质量函数截止的半模式质量（𝑀hf）。
• 评估该方法对数据质量（空间分辨率、信噪比）和透镜几何结构（源红移与透镜红移）的敏感性。
• 评估利用统计透镜信号排除或约束标准ΛCDM冷暗物质模型的可行性。
• 为未来利用欧几里得与CSST等空间望远镜的大型强透镜样本提供约束框架。

提出的方法

• 使用正向建模方法生成模拟的强透镜图像，模拟来自平滑宏观模型的透镜弧，并叠加低质量暗物质晕引起的扰动。
• 改变暗物质粒子质量（𝑚DM）以产生不同的晕质量函数截止，参数化为半模式质量（𝑀hf）。
• 对观测数据和模拟的正向模型应用相同的透镜建模与残差分析流程，提取图像残差。
• 计算观测数据与模拟数据中残差的功率谱，以量化视线方向晕的统计特征。
• 使用近似贝叶斯计算（ABC）比较观测与模拟残差的功率谱，生成𝑀hf的后验分布。
• 结合50个透镜系统的分析结果，提升统计约束精度，同时考虑晕实现中的内在离散性。

实验结果

研究问题

• RQ1 强引力透镜弧中图像残差的功率谱能否在统计上约束暗物质晕质量函数的半模式质量（𝑀hf）？
• RQ2 使用50个哈勃空间望远镜质量的强透镜系统（信噪比上限约50），对暗物质粒子质量（𝑚DM）的预期约束为何？
• RQ3 对𝑚DM的约束如何依赖于透镜几何结构，如源红移和爱因斯坦环的空间尺度？
• RQ4 数据质量（特别是空间分辨率与积分时间）如何影响对低质量暗物质晕的敏感度？
• RQ5 该方法能否区分冷暗物质（CDM）模型（晕质量截止为10⁷ M⊙）与暖暗物质（WDM）模型（𝑚DM = 2.22 keV）？

主要发现

• 在50个哈勃空间望远镜分辨率、信噪比约50的强透镜系统下，该方法在晕质量截止为10⁷ M⊙的CDM类模型中，实现了𝑚DM > 4.10 keV的中位2σ约束。
• 在真实𝑚DM = 2.22 keV（对应𝑀hf = 10⁹ M⊙）的WDM宇宙中，该方法的中位2σ估计值为𝑚DM ∈ [1.43 keV, 3.21 keV]。
• 随着源红移升高，约束显著收紧；当真实𝑚DM = 8.86 keV时，源红移𝑧 = 2.5可提供𝑚DM的2σ上限为5.20 keV。
• 数据质量提升（如更高分辨率或更长积分时间）可带来更紧的约束，即使在欧几里得/CSST质量数据下，对𝑀hf的约束仍具有实际意义。
• 即使存在晕实现中的内在离散性，该方法仍具鲁棒性，但离散性限制了精度，因此需要大样本量才能实现紧密约束。
• 该方法可扩展至未来大规模透镜巡天，为未来可能排除标准ΛCDM模型提供可行路径，若发现低质量晕缺失或数量不足。