Probing Dark Matter with Strongly Lensed Binary Black Hole Mergers: Prospects in the Near Future

问题陈述
暗物质（DM）的本质仍然是宇宙学中最重大的开放性问题之一。虽然标准 $\Lambda$CDM 模型成功解释了大规模观测结果，但它缺乏对冷暗物质（CDM）粒子的直接探测。温暗物质（WDM）提供了一个引人注目的替代方案，由热残留粒子（例如惰性中微子）组成，这些粒子在相对论性阶段脱耦。WDM 引入了一个自由流向尺度，与 CDM 相比，它会抑制低质量暗物质晕的形成。这种抑制改变了晕的质量函数，导致更浅的引力势和更少的低质量结构。

目前对 WDM 粒子质量（$m_{\text{wdm}}$）的约束高度依赖于电磁（EM）探针，如 Lyman-$\alpha$ 森林、高红移类星体吸收光谱以及强引力透镜通量比，这些方法给出的界限通常在 $m_{\text{wdm}} \gtrsim 6\text{--}10$ keV 的范围内。然而，这些约束往往依赖于复杂的天体物理建模（例如恒星形成效率、重子物理学）。本文研究了通过引力波（GW）观测到的强透镜双黑洞（BBH）并合，是否可以提供一种互补的、独立的、且纯粹基于引力的暗物质本质探测手段。具体而言，作者旨在通过分析透镜事件的时间延迟分布，并严格考虑探测器选择效应，来预测即将到来的引力波观测运行（O5、O6 及以后）对 $m_{\text{wdm}}$ 的约束能力。

方法论
作者采用贝叶斯推断框架，利用模拟的透镜 BBH 并合群体来预测对 WDM 反质量（$m^{-1}_{\text{wdm}}$）的约束。该方法包含以下几个关键步骤：

1. 透镜群体建模： 作者使用 Schneider 等人 (2012) 的重标度处方，对 WDM 情景下相对于 CDM 的晕质量函数（HMF）的抑制进行建模。透镜被近似为奇异等温球（SIS），由速度弥散 $\sigma$ 参数化。透镜的数密度源自 HMF，其中较低的 $m_{\text{wdm}}$ 值会抑制低 $\sigma$（低质量）端的分布。
2. 强引力透镜概率： 通过对红移和速度弥散下的透镜数密度进行积分，计算强透镜的光学深度。使用泊松统计得出源至少被透镜一次的概率。
3. 探测器选择效应： 不同于以往假设对所有并合事件都敏感的下一代（XG）探测器的研究，本工作显式地为当前及近未来探测器（从 O4 到 Voyager/XG）建立了选择函数模型。作者通过将内在并合率（使用 Dominik 等人和 Madau & Dickinson 模型）与探测器灵敏度和强透镜概率进行卷积，计算了可探测的透镜比例（$u_{\text{det}}$）和时间延迟分布（$dP_{\text{det}}/d\Delta t$）。
4. 时间延迟特征： 核心物理特征在于，WDM 抑制了低质量晕的丰度，而这些晕正是产生短时间延迟（$\Delta t$）的原因。因此，与 CDM 相比，WDM 时间延迟分布在短 $\Delta t$ 尾部表现出缺失。作者推导出一个近似标度关系，表明抑制阈值 $\Delta t_{\text{min}}$ 与 $m_{\text{wdm}}^{-4.44}$ 成比例。
5. 贝叶斯推断： 作者构建了一个结合两个观测量的似然函数：观测到的透镜比例（$u_{\text{det}}$）和测得的时间延迟集合 $\{\Delta t_i\}$。他们假设对 $m^{-1}_{\text{wdm}}$ 采用均匀先验，并对宇宙学参数（$H_0, \Omega_m, \sigma_8$）进行边缘化处理，以测试是否存在简并性。

主要贡献

• 纳入选择效应： 本文对当前代（O5）和近未来（O6）网络进行了严谨的选择效应处理，证明了这些效对于准确预测至关重要，这与理想化的 XG 情景不同。
• 双观测量约束： 分析结合了透镜比例和时间延迟分布，表明虽然仅靠比例本身只能提供适度的约束，但纳入时间延迟分布后会显著收紧界限。
• 简并性分析： 作者明确展示了 WDM 特征（短 $\Delta t$ 尾部的抑制）与改变宇宙学参数的效果是截然不同的，后者会对整个 $\Delta t$ 范围内的分布产生扭曲。这证实了 $m_{\text{wdm}}$ 可以独立于 $H_0, \Omega_m$ 和 $\sigma_8$ 进行约束。
• 针对近期运行的预测： 该研究为 LIGO-Virgo-KAGRA (LVK) 网络升级提供了具体的预测，超越了理论上的 XG 投影，转向了实际的近未来预期。

结果

• O5 运行： 由于预期的透镜对数量较低，第五次观测运行预计仅能提供适度的约束（$m^{-1}_{\text{wdm}} \lesssim 0.30 \text{ keV}^{-1}$，或 $m_{\text{wdm}} \gtrsim 3.3 \text{ keV}$）。
• O6 运行： 到第六次观测运行，预测的界限将收紧至 $m^{-1}_{\text{wdm}} \lesssim 0.1\text{--}0.2 \text{ keV}^{-1}$（即 $m_{\text{wdm}} \gtrsim 5\text{--}10 \text{ keV}$）。这与目前通过电磁探针获得的最严苛的天体物理约束相当。
• 未来运行（Voyager/XG）： 在 O6 之后的网络预计将把约束提高约一个数量级，达到 $m^{-1}_{\text{wdm}} \lesssim 0.036 \text{ keV}^{-1}$（即 $m_{\text{wdm}} \gtrsim 28 \text{ keV}$），在许多情景下将超越当前的电磁极限。
• 稳健性： 结果显示，其对于 BBH 并合率模型（Dominik 与 MD-nodelay）以及速率定标的变化具有稳健性。此外，对宇宙学参数进行边缘化处理不会显著降低对 $m_{\text{wdm}}$ 的约束。
• 质量测量： 如果暗物质确实是温暗物质（例如 $m_{\text{wdm}} = 8.6 \text{ keV}$），作者表明未来的运行不仅可以设定上限，而且随着样本量的增加，还可以通过提高精度来测量粒子质量。

意义
本文声称，强透镜引力波提供了“一种互补、竞争且独立的暗物质探测手段”。通过仅依赖引力相互作用，该方法绕过了困扰电磁探针的各种天体物理不确定性（如重子反馈和恒星形成效率）。作者强调，到 O6 观测运行为止，引力波观测将达到与目前最佳电磁界限相竞争的精度水平，而未来的运行将提供更严苛的约束。这确立了强透镜引力波作为在不久的将来测试暗物质基本性质的有效且强大的工具。