Two-Dimensional Pulsar Distance Inference from Nanohertz Gravitational Waves

技术摘要：基于纳赫兹引力波的二维脉冲星距离推断

问题陈述
脉冲星计时阵列（PTA）目前在定位纳赫兹连续引力波（CGW）源（如单个超大质量黑洞双星，SMBHB）方面受到限制。主要的瓶颈在于脉冲星距离的不确定性较大。由于脉冲星项（pulsar term）包含有关天球定位的有价值信息，而当前的距离不确定性（通常为数十至数百秒差距）超过了引力波波长（在 10 nHz 时约为 1 秒差距），因此脉冲星项通常被视为一种自噪声而非约束条件。现有的利用 CGW 推断脉冲星距离的方法依赖于来自单个脉冲星的一维（1D）边缘化后验分布，或结合来自多个源的一维后验分布。这些方法未能充分利用包含在联合参数空间中的信息，往往会因为脉冲星距离与引力波源参数（如啁啾质量和频率）之间的简并性，而抹除掉距离后验分布中的多峰结构。

方法论
作者提出了一种二维（2D）联合后验方法，利用多个 CGW 源来推断脉冲星距离。该方法的核心包括：

1. 相位-距离映射： 利用脉冲星项初始相位（$\Phi_p$）与脉冲星距离（$L_p$）之间的解析关系。由于 $\Phi_p$ 是对 $2\pi$ 取模的，单个源会产生具有退化峰值的周期性、多峰距离后验分布。
2. 二维联合构建： 作者并非通过对其他脉冲星进行边缘化或结合一维后验分布，而是利用 $N$ 个 CGW 源的数据，为脉冲星对 $(p, q)$ 构建联合概率密度函数（PDF）。该联合 PDF 定义为：
   $$\text{PDF}_{pq}(L_p, L_q) \propto \pi_{pq}(L_p, L_q) \prod_{n=1}^{N} P_{\Phi pq,n}(L_p, L_q)$$
   其中 $\pi_{pq}$ 是基于视差的联合高斯先验，而 $P_{\Phi pq,n}$ 代表第 $n$ 个源的贡献，将两个脉冲星的相位后验映射到 $(L_p, L_q)$ 平面上。
3. 破缺简并性： 不同的引力波源会根据其频率和天球位置，在距离后验分布上施加截然不同的周期性结构。在二维平面中，这些源会产生“简并带”。虽然单个简并带是平行且具有歧义的，但来自具有不匹配周期性的多个源的简并带之交集会抑制伪峰，从而隔离出真实的距离解。
4. 边缘化与选择： 对二维后验分布关于 $L_q$ 进行边缘化，以获得 $L_p$ 的一维后验分布。为了最大限度地保留信息，作者选择使 68% 置信区间半宽度最小的配对 $(p, q)$ 作为脉冲星 $p$ 的最终约束。

模拟与设置
该方法使用模拟的平方公里阵列（SKA）时代 PTA 进行了测试：

• 数据集： 85 颗脉冲星，其天球位置取自当前数据集，标称距离取自澳大利亚国家望远镜中心脉冲星目录。
• 观测参数： 20 年的时间跨度，2 周的采样间隔，计时噪声建模为白高斯噪声（$\sigma_n = 50$ ns 和 $100$ ns）。
• 引力波源： 一个由 100 个各向同性分布的 SMBHB 组成的目录，具有固定的啁啾质量（$M = 5 \times 10^9 M_\odot$）和频率范围（3–15 nHz）。
• 先验： 基于预测的 SKA 计时视差精度（对于 $\sim 1$ kpc 的脉冲星约为 $\sim 10$ pc）在脉冲星距离上设定高斯先验。
• 推断： 使用 Eryn 系综 MCMC 采样器进行贝叶斯推断，假设引力波源的天球位置是固定的（假设通过电磁对应体进行了外部识别）。

关键结果

• 精度提升： 对于距离约为 $\sim 1$ kpc 的脉冲星，二维方法仅需使用少量 CGW 源即可实现亚秒差距级的距离精度（$\Delta L_p < 1$ pc）。
  • 对于计时噪声 $\sigma_n = 50$ ns，当 $N=4$ 时，实现亚秒差距精度的实现比例上升至 $\sim 72\%$；当 $N=5$ 时，上升至 $\sim 87\%$。
  • 即使在 $\sigma_n = 100$ ns 时，当 $N=5$ 时，仍有 $\sim 47\%$ 的实现达到了该精度。
• 与一维方法的比较： 二维方法显著优于现有的一维边缘化方法。在说明性案例中，二维方法将距离不确定性从 6.5 pc（一维）降低到了 0.4 pc（二维）。一维方法之所以失效，是因为在对第二个脉冲星进行边缘化时，会将距离-距离简并性投影到单一轴上，从而模糊了多峰结构。
• 全阵列影响： 对于包含 85 颗脉冲星的全阵列，结合视差先验与 $N=5$ 个 CGW 源，可将具有亚秒差距精度的脉冲星数量从 4 颗（仅靠视差）增加到中位数 47 颗（约 55%），前提是 $\sigma_n = 50$ ns。
• 对源属性的依赖性： 该方法的有效性取决于源频率和天球位置的多样性。处于“啁啾质量主导”机制（较高频率）下的源会产生具有不同斜率的简并带，这对于打破脉冲星距离之间的相关性至关重要。

意义与主张
论文声称，这种二维联合后验方法能够利用少量 CGW 源实现稳健、高精度的脉冲星距离测量，克服了全多脉冲星联合推断在计算上的不可行性。作者指出，实现亚秒差距精度具有两个主要意义：

1. 宿主星系识别： 它显著改善了 CGW 源的天球定位精度。利用解析缩放关系，作者展示了二维方法可以增加在 10 年观测时间内精确测量的脉冲星数量（超过衍射极限 $N_{dist} \approx 9$），而一维方法通常无法达到这一阈值。这种定位区域的缩小（对于 $N=5$，缩小倍数高达 $\sim 24$）对于即将到来的时域巡天（如 LSST、Roman）中识别宿主星系至关重要。
2. 多信使天文学： 精确的定位有助于识别电磁对应体，从而实现多信使研究，例如通过标准引力波源测量哈勃常数。

作者指出，其分析假设了一个简化的 PTA 环境（白噪声、单个 CGW 源），未来的工作需要量化重叠信号、随机引力波背景（SGWB）以及更真实的噪声模型的影响。他们还强调，该方法依赖于对 CGW 源天球位置的外部确定，这可能通过定向搜索或宿主星系过滤来获得。