The multimessenger view of Pulsar Timing Array black holes with the Horizon-AGN simulation

问题陈述
超大质量黑洞（MBHs）的起源与演化仍不明确，特别是在局部宇宙中的合并过程。虽然脉冲星计时阵列（PTAs）最近探测到了与大质量黑洞双星（MBHBs）随机叠加一致的引力波背景（GWB）信号，但观测到的振幅与理论预测之间存在持续的张力，后者通常倾向于较低的振幅。此外，由于数据集的有限时长以及可能存在的、高于背景的单个可解析连续波（CW）源的存在，解释 GWB 光谱变得更加复杂。目前将理论模型与 PTA 数据进行比较的方法往往依赖于简化的假设，这可能会引入偏差。本文旨在通过利用 Horizon-AGN 余弦模拟来研究对 PTA 波段有贡献的 MBHBs 的性质，生成用于 GWB 和 CW 预测的现实种群，并批判性地评估从计时残差中推断 GWB 光谱属性的方法论。

方法论
本研究利用了 Horizon-AGN 流体动力学宇宙学模拟，该模拟对一个 $142 \, \text{cMpc}^3$ 的体积进行了标准 $\Lambda$CDM 宇宙学建模。作者通过追踪模拟中的 MBH ID 并识别数值合并过程，构建了一个 MBH 合并目录。为了补偿数值合并（分离度约为 $\sim 4 \, \text{kpc}$）与实际并合之间的物理延迟，作者应用了一个后处理流水线，其中包括：

1. 动力学摩擦： 基于等温球近似下的恒星动力学摩擦进行计算。
2. 双星硬化： 通过恒星硬化、环双星吸积盘中的粘性转矩以及引力波（GW）辐射进行建模，遵循既定方案（例如 Sesana & Khan 2015; Dotti et al. 2015）。
3. 偏心率演化： 通过损失锥散射、粘性阻尼和引力波辐射进行追踪，并根据之前的数值结果分配初始偏心率。

为了估计 GWB，作者采用了两种不同的方法：

• 高斯系综（GE）法： 一种确定性的特征应变谱计算，假设存在连续的源分布，这在文献中经常使用。
• 种群法： 一种随机方法，通过基于合并目录的泊松分布抽取旋进双星，生成 2,000 个“宇宙实现”。这考虑了宇宙方差和源的离散性质。

作者使用两种方法模拟了这些种群的 PTA 计时残差：

1. 高斯系综： 通过平均脉冲星响应将特征应变转换为计时残差功率谱密度（PSD）。
2. 种群法： 通过对所有实现中的所有双星的单个地球项 GW 应变贡献求和，显式地模拟有限的观测时长（$T_{\text{obs}}$）和频谱泄漏效应（通过 sinc 函数）。

他们将这些理论预测与代表当前 EPTA（25 颗脉冲星，10.33 年）和未来 SKA（80 颗脉冲星，10.33 年）灵敏度的“玩具”型 PTA 数据集进行了比较。他们还识别了潜在的 CW 候选体（信噪比 $S/N > 3$ 的单个双星），并通过模拟 AGN 发射和宿主星系性质分析了它们的电磁性质（光学、近红外、射电、X 射线）。

主要贡献与结果

• GWB 光谱推断与偏差： 研究强调了将理论模型与 PTA 数据进行比较时的一个显著方法论缺陷。标准的高斯系综法假设无限的观测时长，未能考虑到由有限 $T_{\text{obs}}$ 引起的频谱泄漏。当使用种群法结合现实的 PTA 分析（包括计时模型边际化）时，推断出的 GWB 光谱相对于真实的底层谱，表现出更平坦的斜率和在 $f = 1/\text{yr}$ 处更高的振幅。这种由低频源引起的频谱泄漏驱动的偏差，可能部分解释了理论预测与观测到的 PTA 振幅之间的张力。
• 随机性与源计数： 在低频（$\sim 3 \, \text{nHz}$）处，GWB 由数百到数千个双星主导，证实了随机背景近似的有效性。然而，在高频（$\sim 10\text{--}50 \, \text{nHz}$）处，贡献源的数量显著下降（降至几十个），导致随机性丧失并增加了宇宙方差。偏心种群比圆轨道种群在更高频率处保持随机性，这是由于其在高阶轨道谐波上的辐射。
• 源属性（背景 vs. CW）：
  • 背景： GWB 由红移 $z \in 0.05 - 1$ 且啁啾质量 $M_c \sim 10^{8.5} - 10^{9.5} \, M_\odot$ 的双星主导。它们主要位于质量较大的星系（$M_\star > 10^{11} \, M_\odot$）中，且处于质量约为 $\sim 10^{13} \, M_\odot$ 的晕内。
  • CW 候选体： 单个可解析源（CW 候选体）与一般的背景相比，往往具有更大的质量且位于较低的红移（$z < 0.4$）。它们主要位于大质量星系团和星系群（$M_{\text{halo}} > 10^{13} \, M_\odot$）的中心星系中。
• 电磁对应体：
  • 性质： 大多数 CW 候选体表现为活动星系核（AGN）而非类星体，其特征是低爱丁顿比（$f_{\text{Edd}} < 10^{-2}$）。
  • 可探测性： 在光学和近红外（NIR）波段，大多数 CW 候选体被其宿主星系掩盖，使得光变曲线搜索具有挑战性且存在偏差，研究倾向于寻找罕见的明亮系统。
  • 最佳波段： 射电和 X 射线波段提供了最佳的识别前景。对于 $f_{\text{Edd}} > 10^{-3}$ 的候选体，AGN 发射可以在这些波段中盖过宿主星系。即使对于低吸积系统，乐观估计表明高达 80%（射电）和 30%（X 射线）的候选体可以在宿主之上被探测到。
• 探测概率： 假设 10 年的基线，在 EPTA 级别的灵敏度下，探测到至少一个 $S/N > 3$ 的 CW 候选体的概率估计为 4%，而在 SKA 级别的灵敏度下为 20%。

意义
该论文声称其主要意义在于提供了一个统一的多信使框架，利用最先进的宇宙学模拟来理解 PTA 信号。至关重要的是，它挑战了将理论 GWB 模型与观测数据进行比较的标准方法论。通过证明有限的观测时长和频谱泄漏会人为地抬高推断出的 GWB 振幅并使谱斜率变平，作者指出，观测到的 PTA 现象与天体物理模型之间的表观张力，可能部分是分析技术造成的伪影，而非底层天体物理学的失败。此外，该研究为电磁后续观测提供了具体指导，建议搜索应优先考虑射电和 X 射线波段，并侧重于大质量星系群和星系团，而不是仅仅依赖于类星体系统的光学变率。