Prospects of resolving and localising individual supermassive black hole binaries with pulsar timing arrays: the host ranking challenge

技术摘要：利用脉冲星计时阵列解析与定位单个超大质量黑洞双星的前景

问题陈述
脉冲星计时阵列（PTAs）最近探测到了随机引力波（GW）背景，这很可能源于超大质量黑洞（SMBH）双星族群。下一个重大里程碑是解析出该背景中的单个双星，这将实现多信使天文学。然而，一个主要的挑战是 PTA 有限的天空定位能力。第一个被单独解析出的双星的 90% 置信度定位区域（$A_{90}$）预计将跨越数百平方度，其中可能包含数十万个候选宿主星系。此外，现有的电磁（EM）星系目录在恒星或尘埃污染严重的区域（如银河系平面）是不完整的，从而产生了“缺失”的宿主候选者。目前尚无标准化方法能够根据引力波后验概率和电磁属性高效地对这些海量候选者进行排序，以优先开展后续观测。

方法论
作者使用 holodeck 软件模拟了与 NANOGrav 15 年数据集（NG15）测量的引力波背景一致的现实超大质量黑洞双星族群。他们生成了 1,000 个双星族群的实现，并从每个族群中选择了在未来 PTA 配置下可探测但在 NG15 中不可探测的最强双星。

模拟了三种 PTA 配置以代表未来的灵敏度：

• IPTA_20： 20 年基线，包含 116 颗脉冲星（类似于即将到来的 IPTA 第三次数据发布）。
• IPTA_25（标准配置）： 25 年基线，包含 158 颗脉冲星（在 IPTA_20 的基础上增加了 42 颗新脉冲星）。
• IPTA_30： 30 年基线，包含 200 颗脉冲星。

作者将选定的双星注入到包含随机引力波背景的模拟计时数据中，并使用 QuickCW 流水线来约束双星参数（天空位置、啁啾质量、应变、频率）。他们将“可定位”源定义为信噪比（S/N）$\ge 8$ 的源，这是一个通常需要用于获得良好约束的天空位置的阈值。

为了识别宿主，作者将生成的 $A_{90}$ 定位图与两个全面的全天电磁目录进行交叉匹配：

1. ETG 目录： 包含约 1,940 万个静止型早期型星系（ETGs），具有光度红移和恒星质量（源自 2MPZ 和 WISExSCOS）。
2. AGN 目录： 包含约 450 万个活动星系核（AGN）候选者（来自 Assef 等人 2018 年的 R90 目录），具有 W1 波段视星等。

开发了一个宿主排序系统，为每个候选星系分配一个评分（$H$）。该评分是天空位置概率（$p_{sky}$）与星系电磁属性与引力波后验概率相匹配程度（$p_{cat}$）的乘积。

• 对于 ETG，$p_{cat}$ 是通过星系的红移（转换为光度距离）及其估计的 SMBH 质量（通过标度关系从恒星质量推导而来）得出的。
• 对于 AGN，$p_{cat}$ 是通过 W1 波段视星等得出的，该星等是总质量和光度距离的函数。
  作者通过估算 EM 目录未覆盖的 $A_{90}$ 面积比例，明确计算了“缺失”的宿主数量。

关键结果

• 可探测性与定位：

  • 在未来 0–5 年内（IPTA_20），解析双星（$S/N \ge 4$）的概率约为 21%，但实现良好约束定位（$S/N \ge 8$）的概率仅约为 0.3%。
  • 在未来 5–10 年内（IPTA_25），解析双星的概率上升至约 38%，实现定位的概率约为 3.8%。
  • 在未来 10 年内（IPTA_30），解析概率达到约 51%，实现定位的概率达到约 14.1%。
  • 可定位源预期相对较近（中值光度距离 $d_{lum} \approx 400$ Mpc）且质量极大（$M_{tot} \gtrsim 10^{9.8} M_{\odot}$）。

• 定位区域与宿主数量：

  • 对于标准配置（IPTA_25），38 个潜在可定位的双星中有 30 个满足 $A_{90} < 2,000$ deg$^2$。其 $A_{90}$ 中值为 408 deg$^2$。
  • 这些区域包含中值约为 190,000 个 ETG 和 40,000 个 AGN。
  • 在 IPTA_30 中，中值 $A_{90}$ 缩小至 212 deg$^2$，使中值宿主数量减少了约 40%。

• 覆盖不完整性：

  • EM 覆盖不完整是一个显著的限制。在标准配置中，缺失的 ETG 宿主中值约为 20,000 个，缺失的 AGN 宿主中值约为 6,000 个。
  • 对于某些源，$A_{90}$ 部分或完全位于目录足迹之外（例如，与银河系平面重叠），这使得在没有新巡天的情况下无法识别真实的宿主。

• 排序效率：

  • 排序方法成功排除了属性与引力波后验不符的星系。
  • ETG： 约 47%（IPTA_25）至 61%（IPTA_30）的 ETG 因红移或 SMBH 质量不兼容而被排除。其排序效率（$E_{95-50}$，即第 95 百分位数得分与中值之比）在 IPTA_25 中约为 12.4，在 IPTA_30 中约为 17.2。
  • AGN： 仅约 0.3%（IPTA_25）至 3.8%（IPTA_30）的 AGN 被排除。其排序效率显著较低（IPTA_25 约为 3.7），因为视星等是质量和距离的退化组合，且 PTA 对啁啾质量的约束较差。

意义与主张
本文声称提供了对第一个 PTA 解析的 SMBH 双星的候选宿主数量及不完整天空覆盖影响的首次现实量化。作者表明，虽然候选者数量庞大，但基于引力波后验和电磁属性的排序系统可以有效地缩小搜索范围，特别是对于具有红移和恒星质量信息的 ETG。

研究强调：

1. 多信使发现是可能的，但也具有挑战性： 在 10 年内定位源的概率是不可忽视的（$\approx 14\%$），但定位区域仍然很大。
2. 目录完整性至关重要： 由于巡天间隙（如银河系平面）导致的数千个缺失宿主的存在，构成了宿主识别中假阴性的风险。
3. 数据质量决定排序成功与否： 排序方法的效率高度依赖于独立红移和 SMBH 质量估计的可获得性。目前无法高效对 AGN 进行排序的原因在于 AGN 目录中缺乏红移/质量数据，以及 PTA 对啁啾质量的约束能力较弱。
4. 未来改进： 作者建议，未来通过改进脉冲星距离测量（以利用引力波信号中的脉冲星项）以及开发带有红移和质量估计的 AGN 目录，将显著增强定位和排序能力。

论文得出结论，虽然第一次多信使探测已近在咫尺，但其成功不仅取决于 PTA 的灵敏度，还取决于电磁目录的完整性和均匀性，以及用于优先开展后续观测的有效排序流水线的开发。