[论文解读] New algorithms for radio pulsar search

Kendrick M. Smith|arXiv (Cornell University)|Oct 21, 2016

Pulsars and Gravitational Waves Research被引用 3

一句话总结

本文提出了四种新算法——最优恒周期搜索、相干树搜索、半相干搜索和分层搜索——通过利用递归树结构和相位相干信号组合，显著降低了射电脉冲星搜索的计算成本。这些方法在每个模型上实现 O(1) 的计算成本，接近最优检测性能，使大规模时间序列的搜索复杂度降低数个数量级。

ABSTRACT

The computational cost of searching for new pulsars is a limiting factor for upcoming radio telescopes such as SKA. We introduce four new algorithms: an optimal constant-period search, a coherent tree search which permits optimal searching with O(1) cost per model, a semicoherent search which combines information from coherent subsearches while preserving as much phase information as possible, and a hierarchical search which interpolates between the coherent and semicoherent limits. Taken together, these algorithms improve the computational cost of pulsar search by several orders of magnitude. In this paper, we consider the simple case of a constant-acceleration phase model, but our methods should generalize to more complex search spaces.

研究动机与目标

• 解决下一代射电望远镜（如 SKA）中脉冲星搜索日益增长的计算负担问题，其中搜索空间大小随时间序列长度的增加而表现极差的可扩展性。
• 克服传统功率谱折叠方法的局限性，该方法丢弃相位信息，相比相干方法性能次优。
• 开发更快、接近最优的搜索算法，在降低计算成本的同时保留相位相干性，使大规模时间序列的计算成本从 O(S) 降低至 O(T)。
• 通过引入在相干与非相干搜索策略之间插值的可扩展方法，实现对微弱脉冲星的实际检测。
• 提供一个可扩展的框架，适用于复杂模型（如多项式相位模型或双星脉冲星系统），突破恒定加速度假设的限制。

提出的方法

• 使用递归树算法，以每个模型 O(1) 的计算成本计算恒定加速度模型下的最优相干重叠统计量，达到计算成本的理论下限。
• 通过将数据划分为短的相干段，对每段计算相干统计量，并利用保持相位信息的递归方法组合，实现半相干搜索。
• 定义一种分层搜索，逐步增加相干时间，从而检测在长 timescale 上偏离恒定加速度模型的脉冲星。
• 对半相干统计量应用递归关系，尽管理论上搜索空间呈指数增长，仍可实现 O(T) 的评估时间。
• 通过调节相干时间 Tc，在相干与半相干模式之间插值，实现灵敏度与计算成本之间的自适应权衡。
• 利用分层结构标记“信号丢失”——即在短 Tc 下显著但在长 Tc 下显著性下降的候选体——指示潜在的双星或非恒定加速度脉冲星。

实验结果

研究问题

• RQ1 我们能否设计一种搜索算法，在实现接近最优灵敏度的同时，将恒定加速度模型下每个模型的计算成本从 O(S) 降低至 O(1)？
• RQ2 如何在多个相干段之间保留相位信息，以实现与时间序列长度呈线性关系的计算复杂度的半相干搜索？
• RQ3 分层搜索策略在多大程度上能够检测出对恒定加速度或多项式相位模型拟合较差的脉冲星？
• RQ4 所提出的方法能否扩展至高阶相位模型（如 3–4 次多项式），与完整相干搜索相比显著降低计算成本？
• RQ5 与传统的功率谱折叠或非相干叠加相比，新算法在性能和对 RFI 及假阳性信号的鲁棒性方面表现如何？

主要发现

• 递归树算法以每个模型 O(1) 的计算成本计算恒定加速度模型下的最优相干重叠统计量，达到计算成本的理论下限。
• 半相干搜索方法实现了对加速度缓慢变化模型的最优检测，尽管理论上搜索空间呈指数增长，计算成本仍保持线性 O(T)。
• 分层搜索成功识别出偏离恒定加速度模型的脉冲星，通过检测在相干时间增加时显著性下降的统计信号实现。
• 在分层搜索中，完整的相位相干信噪比（Ŝ）得以保留，从而可通过与相干检测阈值比较，实现对假阳性的稳健排除。
• 这些方法对弱、未掩蔽的 RFI 具有鲁棒性，因为 RFI 通常以非相干方式贡献，难以模拟真实相干脉冲星信号的尖锐似然峰。
• 所提出的算法预计可将计算成本降低数个数量级，使下一代望远镜（如 SKA 和 CHIME）大规模时间序列的最优脉冲星搜索成为现实。