New Approach to Pulsar Timing Array Data Combination
FrankenStat I:一种脉冲星计时阵列数据组合的新方法
In 2023, after more than two decades of searching, pulsar timing array (PTA) collaborations around the world announced evidence for a stochastic gravitational wave background. It was quickly followed by work from the International Pulsar Timing Array (IPTA), demonstrating that the results of regional collaborations were consistent with each other. The combination of these datasets is still ongoing and represents a significant investment of time and expertise. In that IPTA comparison, authors of this Letter combined the separate datasets in the standard PTA optimal detection statistic for cross-correlations incoherently, that is, the data were combined without fitting a merged timing model across all PTA datasets, treating datasets of the same pulsar as independent, and neglecting the "same pulsar, different datasets" cross-correlations. This Letter refines that method by extending its core ideas beyond detection statistics and into a full, general data-combination method. We have demonstrated its efficacy and extreme efficiency on simulated data. This new method, "FrankenStat," is very similar in sensitivity and parameter-constraining power to traditional data-combination methods while completing the full data combination in just a few minutes.
展开 ▾将耗时数月的 IPTA 数据组合改为对象级拼接,可沿用现有 PTA 分析代码;仿真显示 SNR、p 值和灵敏度曲线与传统组合差异小于 1%。
此前 IPTA 对区域脉冲星计时阵列结果的比较 Agazie+ 2024 采用非相干组合并显示联合数据可显著提升信噪比;Lite 方法 Larsen+ 2025 则通过品质因数快速选择单阵列数据以缩短初步分析时间。FrankenStat 在这些思路基础上,借鉴早期单脉冲星 IPTA DR2 的 Borg 方法 Schult 2021,将拼接残差和分块对角计时模型扰动矩阵嵌入标准 PTA 似然框架 Johnson+ 2024,无需统一计时模型即可保留全阵列灵敏度。仿真表明该方法与完全组合的灵敏度差异小于 1%,且计算只需数分钟。未来结合 enterprise/discovery 等对象化分析接口以及自动微分和 GPU 加速,有望使 IPTA 数据组合与下游分析由数月缩短至数周,并推广到任意跨数据集信号搜索。
预印本 2025-12-16 · 刊出 2026-07-21 · 收录 2026-08-20