One Year and One Night to 1% in $\mathbf{H_0}$: Efficient Spectroscopic Strategy for Dark Siren Cosmology

技术摘要：一年与一晚实现 $H_0$ 的 1% 精度的研究

问题陈述
哈勃常数 ($H_0$) 仍是宇宙学中一个显著的张力问题，来自宇宙微波背景（CMB）的测量值与局部距离阶梯的测量值之间存在约 $5\sigma$ 的差异。引力波（GW）标准汽笛为解决这一张力提供了一条独立的途径。虽然“明汽笛”（具有电磁对应体的事件）提供了直接约束，但大多数紧凑双星并合（CBC），特别是双黑洞（BBH）并合，缺乏可探测的电磁对应体。这些“暗汽笛”需要与引力波定位体积内的红移星系目录进行统计关联。

目前的暗汽笛分析依赖于光度星系目录（例如 GLADE+、DES），这些目录在深度或天区覆盖方面存在不完整性，且具有较大的光度红移不确定性 ($\sigma_z \sim 10^{-2}$)。虽然光谱红移 ($\sigma_z \sim 10^{-3}–10^{-4}$) 提供了更高的精度，但对于体积限制型星系目录而言，最优深度以及为实现良好定位事件的完整光谱覆盖所付出的观测成本尚未得到系统性的量化。具体而言，目前尚不清楚浅层目录（例如 $r \sim 19$）是否会引入偏差或低估不确定性（相比于深层目录 $r \sim 24$），以及是否存在一个可行的光谱观测任务，以实现亚 2% 的 $H_0$ 精确度。

方法论
作者采用贝叶斯统计框架，利用一组 $N_{GW}$ 个引力波事件和一个星系目录（CAT）来推断 $H_0$。后验概率由引力波应变数据的联合似然函数与星系识别概率构建而成。

1. 模拟与注入：

   • 使用 Uchuu-$\nu^2$GC 暗物质 N 体模拟构建了一个模拟星系目录，并通过半解析模型分配星系属性。
   • 该模拟生成了一个延伸至红移 $z=0.5$ 的光锥，包含 $i$ 波段星等最高达 24.1 的星系。
   • 抽取了 $245,760$ 个星系，并根据恒星质量进行加权，以优先选择作为潜在 BBH 宿主的质量较大的星系。
   • BBH 并合属性（质量、自旋）采用 Broken Power Law + 两个峰值和高斯分量自旋模型分配。并合率遵循 Madau-Dickinson 星形成历史。
   • 使用 Fisher 矩阵（GWBENCH）估计 GW 参数，并通过 IMRPhenomXPHM 波形族通过贝叶斯推断（Bilby）进行精细化。

2. 网络与事件选择：

   • 分析了两种探测器网络配置：
     • HLI#：LIGO Hanford、LIGO Livingston 和位于 A# 灵敏度的 LIGO-India。
     • HLV+：LIGO Hanford、LIGO Livingston 和 Virgo，处于 A+ 灵敏度。
   • 从 10 年的模拟观测中，选择了定位最好的事件（由天区面积 $\Delta\Omega_{90}$ 定义）。本研究重点关注单年观测中定位最好的前 10 个事件。

3. 目录深度分析：

   • 作者将体积限制型目录的表观星等极限 ($r_{app}$) 从 $r=17.5$ 变化到 $r=24$。
   • 他们评估了不同事件数量（3、5、10、15）下的分数不确定性 ($\sigma_{H_0}/H_0$) 和统计偏差（后验中位数与注入值的偏移）。

4. 观测成本估算：

   • 论文评估了获取 10 个定位最好事件的光谱红移所需的望远镜时间。
   • 对多个设施进行了建模：DESI、Subaru/PFS、HET/VIRUS、AAT/2dF、4MOST 和 Magellan/IMACS。
   • 基于仪器规格、目标密度及光纤/狭缝容量计算了曝光时间和开销。

主要贡献与结果

1. 最优目录深度：

   • 研究表明，对于定位良好的事件，一个表观星等极限为 $r \ge 19$ 的浅层体积限制型目录，在 $H_0$ 不确定性方面并不会比深层目录（$r \sim 24$）产生低估。
   • 对于前 10 个定位最好的事件，中位精度从 $r=17.5$ 时的 $0.92\%$ 仅略微提升至 $r=24$ 时的 $1.0\%$。
   • 然而，浅于 $r=19$ 的目录表现出显著的统计不稳定性和潜在偏差，因为某些事件可能匹配不到任何星系，从而导致无信息的后验分布。因此，$r=19$ 被确定为深度与稳定性之间的最优折衷点。

2. 精度与偏差：

   • A# 网络 (HLI#)： 结合一年内 10 个定位最好的事件，可获得 $H_0$ 的无偏测量，中位精度为 $0.98^{+0.56}_{-0.32}\%$。精度分布的 84 百分位值为 $1.54\%$，安全低于 2% 的目标。该测量在统计上是无偏的，中位偏移量与零一致（$-0.03^{+0.15}_{-0.45}\%$）。
   • A+ 网络 (HLV+)： 使用相同的 10 事件策略，HLV+ 网络得到的 $H_0$ 中位精度为 $2.63^{+1.87}_{-1.53}\%$。只有 32% 的实现达到了 2% 的目标，尽管 88% 的实现达到了 5% 或更低。

3. 观测可行性：

   • 针对 10 个定位最好事件的专门光谱随访任务，在理想条件下至少需要一个夜晚的望远镜时间。
   • 在保守假设下（考虑调度、开销和变化的观测目标密度），该任务需要大约五个夜晚（总计约 $20.5$ 小时）。
   • 诸如 DESI、Subaru/PFS 和 Magellan/IMACS 等仪器能够在单个夜晚内完成此项调查，而 HET/VIRUS 虽然速度较慢，但在缺乏现有光度目录的区域具有独特优势。

意义与主张
本文声称确定了一种高效且可行的策略，通过暗汽笛实现亚 1% 的哈勃常数测量。其主要意义在于证明了：

• 成本效益： 实现高精度宇宙学并不需要深层的全天区光谱调查。针对每年仅 10 个定位最好的事件进行浅层（$r \sim 19$）靶向光谱随访是足够的。
• 可行性： 所需的观测时间（1–5 个夜晚）完全在标准望远镜分配限制之内，这使得该策略对于即将到来的 A# 和 A+ 观测运行具有即时可操作性。
• 鲁棒性： 结果对于目录深度（只要 $r \ge 19$）以及特定的宿主星系恒星质量权重是稳健的，为在等待下一代引力波探测器或大规模光谱调查之前解决哈勃张力提供了一条切实路径。

作者指出，本征速度会引入微小的额外红移不确定性，并且大约 20% 的定位良好事件落在避日带（zone of avoidance）内，这将按比例增加所需的观测时间。然而，核心结论仍然是，“一年与一晚”（乐观情况下）或“一年与五晚”（保守情况下）的行动可以交付一个精确且无偏的 $H_0$ 测量。