[论文解读] A nested sampling code for targeted searches for continuous gravitational waves from pulsars

M. Pitkin, M. Isi+2 more|arXiv (Cornell University)|May 24, 2017

Pulsars and Gravitational Waves Research被引用 12

一句话总结

本文提出了一种基于嵌套采样法的代码 lalapps_pulsar_parameter_estimation_nested（lppen），用于利用地面探测器对已知脉冲星进行目标连续引力波搜索。该方法通过计算贝叶斯证据（边缘似然），实现稳健的贝叶斯参数估计与模型选择，能够将信号模型与噪声进行比较，其效率与一致性优于以往的MCMC和网格搜索方法，同时支持检测显著性评估。

ABSTRACT

This document describes a code to perform parameter estimation and model selection in targeted searches for continuous gravitational waves from known pulsars using data from ground-based gravitational wave detectors. We describe the general workings of the code and characterise it on simulated data containing both noise and simulated signals. We also show how it performs compared to a previous MCMC and grid-based approach to signal parameter estimation. Details how to run the code in a variety of cases are provided in Appendix A.

研究动机与目标

• 开发一种更高效、更稳健的方法，用于对已知脉冲星进行目标连续引力波搜索中的参数估计。
• 通过计算贝叶斯证据以比较信号模型与高斯噪声，实现模型选择。
• 克服以往MCMC和网格方法的局限性，包括效率低下和缺乏检测显著性判据。
• 与LALInference库集成，确保与现有LIGO/Virgo分析流程的兼容性。
• 支持在相位演化参数的小范围内进行搜索，当信号不完全匹配已知脉冲星星历表时亦可适用。

提出的方法

• 通过LALInference库使用嵌套采样法，计算信号模型的贝叶斯证据（边缘似然）。
• 采用针对脉冲星连续波信号参数空间（包括振幅、相位和自旋参数）量身定制的提议分布。
• 利用嵌套采样算法评估信号参数的后验概率分布。
• 通过比较信号模型与纯高斯噪声的证据，计算检测显著性的似然比。
• 在扫描相位演化参数时，采用降秩求积法以实现高效的似然评估。
• 使用已注入已知信号的模拟数据，与MCMC和网格方法进行结果验证。

实验结果

研究问题

• RQ1 嵌套采样能否为针对连续波搜索的参数估计提供比MCMC和网格方法更高效、更一致的替代方案？
• RQ2 该嵌套采样方法是否能可靠地计算信号模型的贝叶斯证据，并实现与噪声的模型选择？
• RQ3 与以往方法相比，该新代码在收敛性、准确性和计算成本方面的性能如何？
• RQ4 当引力波相位演化与已知脉冲星星历表存在轻微偏差时，该方法能否检测到信号？
• RQ5 相干信号模型与非相干信号模型之间的似然比，随信噪比（SNR）的标定行为如何？

主要发现

• lppen代码成功计算了信号模型的贝叶斯证据，实现了与噪声的模型选择，提供了正式的检测显著性度量。
• 与MCMC和网格方法相比，嵌套采样在收敛速度和一致性方面表现更优，且无需大量调参。
• 在模拟数据中，该代码能准确恢复注入的信号参数，展示了参数估计的鲁棒性与可靠性。
• 相干信号与噪声之间的似然比近似遵循 exp(ρ²) 的标定关系，与预期的灵敏度标定一致。
• 相干信号与非相干信号之间的似然比大致与 ρ 成正比，表明其对信噪比的依赖性较弱，与理论预期一致。
• 即使在低信噪比条件下，该方法仍保持良好性能，且似然比的解析近似与模拟结果高度吻合。