[论文解读] Detection of Gravitational Waves through Observations of a Group of Pulsars

А. Е. Родин|arXiv (Cornell University)|Dec 30, 2010

Statistical and numerical algorithms参考文献 12被引用 7

一句话总结

本文提出了一种利用毫秒脉冲星计时残差检测随机引力波背景的新方法，通过奇异谱分析（SSA）将信号分解为具有不同光谱特征的分量。该方法在第8和第9个分量中识别出具有统计显著性的角相关性（ρ = 0.92 ± 0.10），表明存在能量密度为 Ω₉h² ~ 10⁻⁸ 的引力波背景，超过以往的上限。

ABSTRACT

We suggest a new approach to the detection of gravitational waves using observations of a group of millisecond pulsars. In contrast to the usual method, based on increasing the accuracy of the arrival times of pulses by excluding possible distorting factors, our method supposes that the additive phase noise that is inevitably present even in the most accurate observational data has various spectral components, which have characteristic amplitudes and begin to appear on different time scales. We use the "Caterpillar" (Singular Spectral Analysis, SSA) method to decompose the signal into its components. Our initial data are the residuals of the pulse arrival times for six millisecond pulsars. We constructed the angular correlation function for components of the decomposition of a given number, whose theoretical form for the case of an isotropic and homogeneous gravitational-wave background is known. The individual decomposition components show a statistically significant agreement with the theoretical expectations (correlation coefficient $ρ=0.92\pm 0.10$).

研究动机与目标

• 开发一种利用毫秒脉冲星独立于极高计时精度的随机引力波背景检测替代方法。
• 解决在脉冲星计时数据中区分引力波信号与固有噪声及系统误差的挑战。
• 检验对计时残差进行谱分解是否能揭示与引力波背景一致的关联分量。
• 评估从分解分量中得出的角相关函数的统计显著性。

提出的方法

• 该方法使用奇异谱分析（SSA），一种数据驱动的分解技术，将脉冲星计时残差分解为具有不同光谱特征的分量。
• SSA直接从时间序列构建正交基函数，无需对噪声结构作先验假设。
• 对每个分解分量计算角相关函数 ζ(θ) = (3/2)x log(x) − x/4 + 1/2 + 1/2 δ(x)，其中 x = (1 − cosθ)/2，以检验是否存在类似引力波的相关性。
• 分析聚焦于第8和第9个分量，其与各向同性引力波背景的理论预期高度一致。
• 该方法考虑了仪器误差，并使用线性插值填补40天平均时间序列中的间隙。
• 通过相关系数 ρ 和 p 值估计（Pr(ρ=0.92, N=19) < 10⁻⁷）评估统计显著性。

实验结果

研究问题

• RQ1 对脉冲星计时残差进行谱分解是否能揭示与随机引力波背景一致的角相关性分量？
• RQ2 计时残差的分量是否表现出光谱特征（如红噪声、蓝噪声）以指示不同的物理起源？
• RQ3 特定分量中观测到的角相关性是否具有统计显著性，并与引力波背景的理论预测一致？
• RQ4 即使尚未被证实，该方法是否能在高于以往上限的水平上探测到引力波背景？

主要发现

• 第8和第9个SSA分量的角相关函数与理论Hellings-Downs曲线高度一致，相关系数为 ρ = 0.92 ± 0.10。
• 该相关性由偶然导致的概率小于 10⁻⁷，表明具有强烈的统计显著性。
• 第8和第9个分量的关联振幅为 0.5 μs，对应于六年的观测跨度下引力波背景能量密度 Ω₉h² ~ 10⁻⁸。
• 该能量密度估计值显著高于以往的上限，后者通常为 Ω₉h² < 10⁻¹⁰。
• 分量1–5表现出红噪声谱，分量6–9具有中频峰值，分量10–12则呈现蓝噪声特征，表明残差中存在多样的物理起源。
• 重构分量的误差占仪器误差的19%至32%，表明分解具有稳健性。