[论文解读] Search for Continuous Gravitational Wave Signals in Pulsar Timing Residuals: A New Scalable Approach with Diffusive Nested Sampling

Yu-Yang Songsheng, Yi-Qian Qian+6 more|arXiv (Cornell University)|Sep 1, 2021

Pulsars and Gravitational Waves Research参考文献 81被引用 8

一句话总结

本文提出扩散嵌套采样（DNest）作为一种可扩展的贝叶斯方法，用于在脉冲星计时残差中检测单个超大质量黑洞双星的连续引力波信号。结果表明，DNest能高效处理高维参数空间——尤其在存在未知脉冲星相位参数的情况下——即使在超过100颗脉冲星的情形下，也能实现高精度和强鲁棒性的源检测，并成功恢复源参数，包括天球位置、振幅、频率和并合时间，涵盖轨道演化的双星系统。

ABSTRACT

Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semi-analytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling shown the capability of coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including $\mathcal{O}(10^2)$ pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also high beneficial to include DNest in current toolboxes of PTA analysis.

研究动机与目标

• 解决在脉冲星时标阵列（PTA）中连续引力波检测的贝叶斯推断所面临的高维、多峰参数空间挑战。
• 克服现有方法（如并行退火MCMC和似然最大化/边缘化）在可扩展性及对信号模型假设方面的局限性。
• 开发并验证一种通用的、可扩展的贝叶斯框架，采用扩散嵌套采样（DNest），可在无需简化假设的前提下处理复杂信号模型。
• 展示该方法在检测和表征单个连续引力波源方面的能力，包括具有轨道演化的双星系统和未知脉冲星相位参数的源。
• 实现在大规模PTA中同时检测多个源的能力，这对未来如FAST和SKA等新一代射电望远镜至关重要。

提出的方法

• 将扩散嵌套采样（DNest）这一贝叶斯推断技术应用于脉冲星计时残差分析中的高维、多峰后验分布探索。
• 使用时域信号模型对来自单个超大质量黑洞双星的连续引力波信号进行建模，该模型包含脉冲星相位参数和演化的轨道参数。
• 利用信号模型计算作为时间函数的计时残差，包含由于脉冲星距离引起的延迟的地球项和脉冲星项。
• 通过依赖引力波驱动的轨道衰减推导出时间相关的轨道频率和相位，实现双星轨道的时间演化建模。
• 在关键参数（如啁啾质量、光度距离和并合时间）上采用非信息性、对数均匀先验，并对脉冲星距离使用具有20%不确定性的高斯先验。
• 通过使用DNest对完整后验分布进行采样实现盲源检测，该方法通过带有扩散更新的嵌套采样算法，高效探索复杂、高维的参数空间。

实验结果

研究问题

• RQ1 DNest能否在高维、多峰参数空间的脉冲星计时残差中有效检测来自单个超大质量黑洞双星的连续引力波信号？
• RQ2 在大规模PTA设置下，DNest在精度、鲁棒性和效率方面与传统MCMC方法（如并行退火）相比表现如何？
• RQ3 在不依赖简化假设的前提下，DNest在处理复杂信号模型（如具有轨道演化的双星系统和未知脉冲星相位参数）方面的能力有多强？
• RQ4 DNest能否在单次分析中同时检测多个源，其可扩展性与现有方法相比如何？
• RQ5 当双星处于非恒定频率状态时，DNest能否准确恢复关键源参数（如天球位置、振幅、频率和并合时间）？

主要发现

• DNest在包含超过100颗脉冲星的模拟数据中成功检测到信噪比（SNR）约为120的单个连续引力波源。
• 源的天球位置（赤经和赤纬）后验分布将真实源定位在约3 deg²的95%可信区域内，且与后验中心的分离为0.7°。
• 尽管与轨道倾角存在强相关性，计时残差的振幅仍被约束在约0.1 dex的不确定度内。
• 引力波的角频率被精确约束，而并合时间即使超过观测时长，也因地球项与脉冲星项之间频率差异而被高精度确定。
• DNest准确恢复了短并合时间（≲10 yr）源的参数，证明其具备在识别电磁对应体时探测标准烛光以用于宇宙学研究的能力。
• 在建模轨道演化双星系统时，该方法保持鲁棒且高效，表明DNest可直接推广至包含时变轨道参数的模型，无需额外近似。