A model-independent measurement of the Hubble constant from gravitational-wave standard sirens and electromagnetic observations

技术摘要：基于引力波标准汽笛与电磁观测的无模型哈勃常数测量

问题陈述
本文探讨了“哈勃张力”问题，即基于 $\Lambda$CDM 模型从宇宙微波背景（CMB）得出的哈勃常数测量值（$H_0 = 67.4 \pm 0.5$ km s$^{-1}$ Mpc$^{-1}$）与使用 Ia 型超新星（SNe Ia）及局部距离阶梯进行的晚期宇宙直接测量值（$H_0 = 73.04 \pm 1.04$ km s$^{-1}$ Mpc$^{-1}$）之间存在的显著差异，该差异约为 $5\sigma - 6\sigma$。传统的测量方法通常依赖于特定的宇宙学模型或关于声速标度（sound horizon scale）的假设。作者旨在提供一种完全无模型的 $H_0$ 确定方法，通过仅依赖观测数据并结合电磁（EM）和引力波（GW）数据集，来避免这些理论偏差。

方法论
研究采用了高斯过程（GP）回归，这是一种非参数贝叶斯方法，用于直接从观测数据中重建光度距离函数 $D_L(z)$ 及其一阶导数，而不假设预定义的函数形式（例如 $\Lambda$CDM）。

1. 框架： 作者使用了 GaPP Python 代码。他们将 $D_L(z)$ 建模为一个由均值函数和协方差核定义的高斯过程。由于平方指数核（squared exponential kernel）具有更好的稳定性且导数性质良好，因此选择了它而非 Matérn 核。
2. $H_0$ 的推导： 哈勃常数是通过利用重建函数进行解析推导得到的，其关系式为：
   $$H_0 = c \left[ \frac{d}{dz} \left( \frac{D_L(z)}{1+z} \right) \bigg|_{z=0} \right]^{-1}$$
   这要求在 $z=0$ 处对 $D_L(z)$ 及其导数进行精确重建。
3. 数据整合： 分析结合了以下数据：
   • 引力波： 来自 GWTC-3 星表的数据（90 个候选事件，$p_{astro} > 0.5$），提供了独立于宇宙距离阶梯的绝对光度距离。
   • 电磁探针：
     • 宇宙计时器（CC）：直接的 $H(z)$ 测量值，并转换为 $D_L(z)$。
     • DESI BAO DR2：重子声学振荡数据，其中声速标度 $r_d$ 被视为待定参数（nuisance parameter）并进行边缘化处理，以保持模型的无模型特性。
     • Ia 型超新星：三个汇编数据集（Pantheon+SH0ES, Union 3.0, 以及 DES-SN5YR）。
     • 伽马射线暴（GRBs）：作为高红移指示器，通过 Amati 关系使用。
4. 误差处理： 该方法通过在低红移（$z < 0.05$）处以平方和形式加入本征速度误差，来处理本征速度误差。至关重要的是，GP 框架显式地计算并传播了重建函数 $D_L(z)$ 与其导数 $D'_L(z)$ 之间的协方差，以确保统计一致的误差估计。

核心贡献

• 混合无模型方法： 本研究是首批将 GWTC-3 标准汽笛数据整合进 GP 框架，并与多样化的 EM 数据集共同用于约束 $H_0$，且不依赖于声速标度或特定暗能量模型的尝试之一。
• 量化 GW 的影响： 研究明确分离了 GW 数据对 $H_0$ 精度的贡献。通过对比包含和不包含 GW 数据的 GP 运行结果，作者证明添加 GWTC-3 数据使 $H_0$ 的 $1\sigma$ 不确定度降低了约 15–20%（在特定组合中从 $1.02$降至$0.58$ km s$^{-1}$ Mpc$^{-1}$）。
• 协方差处理： 在 GP 框架内实现解析微分，使得能够严格处理距离函数与其导数之间的相关性，这是在 $z=0$ 处最小化偏差的关键因素。

结果
无模型重建得到的 $H_0$ 值介于 Planck 和 SH0ES 结果之间，但表现出对 SH0ES 值（$H_0 = 73.04 \pm 1.04$ km s$^{-1}$ Mpc$^{-1}$）的强烈偏好。

• 对于大多数数据组合，所得 $H_0$ 与 SH0ES 结果之间的张力小于 $2\sigma$，表明不存在统计学上的显著差异。
• 与仅使用 EM 的分析相比，加入 GW 数据显著收紧了约束。
• 其中一个特定组合（GWTC3+CC+DESI+DESY5）略微倾向于 Planck 的测量值，但整体趋势与局部距离阶梯的结果一致。

意义与主张
作者审慎地指出，虽然这项研究并未最终解决哈勃张力，但它提供了一个稳健且一致的 $H_0$ 估计值，且仅使用观测输入便倾向于 SH0ES 的测定结果。其工作的意义在于：

1. 验证 GW 数据： 证明了即使在当前存在局限性（红移覆盖范围低、统计不确定度高）的情况下，GW 标准汽笛仍能提供绝对定标，当与 EM 数据结合时，能显著提高 $H_0$ 测量的精度。
2. 一致性检查： 在无模型框架下，确认了 GW 与 EM 测量结果之间的一致性。
3. 未来用途： 建立了一套可扩展的方法论，可应用于未来来自下一代探测器的更大规模 GW 星表以及更精确的 EM 数据集，为在没有理论偏差的情况下研究哈勃张力的起源提供路径。

论文承认了当前的局限性，例如对于没有电磁对应体的双黑洞（BBH）事件，在红移推断上依赖于群体层面的假设，并指出对这些系统误差进行严格处理仍是未来的研究方向。