Searches for Continuous Gravitational Waves from Supernova Remnants in the first part of the LIGO-Virgo-KAGRA Fourth Observing run

本文利用 LIGO-Virgo-KAGRA 第四观测期前八个月的数据，通过五种分析流程对 15 个邻近超新星遗迹进行了无先验信息的定向连续引力波搜索，未发现信号证据，并针对最邻近的 Vela Jr. 源设定了目前最灵敏的应变振幅上限及相应的中子星椭圆度和 r 模振幅限制。

要点

这是一篇关于寻找“宇宙心跳”的科学报告。为了让你轻松理解，我们可以把这篇论文想象成一群天文学家在“听诊”宇宙的故事。

🌌 故事背景：我们在找什么？

想象一下，宇宙中有一些极其致密、旋转极快的“死星”，叫做中子星。它们就像宇宙中的超级陀螺，每秒能转几百圈。

如果这些陀螺长得不够圆（哪怕只有一点点像土豆一样的凸起），当它们疯狂旋转时，就会像洗衣机脱水时那样，搅动时空，产生一种看不见的涟漪，叫做**“连续引力波”**。

• 之前的发现：以前我们主要听到的是“砰”的一声（比如两个黑洞撞在一起），那是瞬态引力波。
• 现在的目标：这次我们要找的是持续不断的“嗡嗡”声（连续引力波）。这就像是在嘈杂的摇滚音乐会上，试图听清远处一只蚊子持续发出的微弱嗡嗡声。

🎯 我们的任务：监听 15 个“嫌疑犯”

这篇论文记录了科学家们在2023 年 5 月到 2024 年 1 月期间（也就是 LIGO-Virgo-KAGRA 第四次观测运行的前八个月），利用地球上最灵敏的“耳朵”（引力波探测器），专门盯着15 个年轻的超新星遗迹（也就是恒星爆炸后留下的“废墟”）。

为什么选这些地方？因为恒星爆炸后，中心往往会留下一颗年轻、转得飞快、可能长得有点歪的中子星。它们是产生“嗡嗡声”的最佳嫌疑人。

🛠️ 我们用了什么工具？（5 种不同的“听诊器”）

为了不漏掉任何声音，科学家们动用了5 种不同的算法（管道），就像用 5 种不同的方法去过滤噪音：

1. BSD (带采样数据)：像是一个**“频率跳跃侦探”**，它把数据切成小块，一块一块地找规律。
2. PyStoch (随机相关)：像一个**“交叉验证员”**，它比较两个探测器（一个在华盛顿，一个在路易斯安那）的数据，看是否有同步的微弱信号。
3. Viterbi (单谐波)：像一个**“追踪猎犬”**，它假设信号可能会因为中子星转得越来越慢而稍微变调，然后一路追踪下去。
4. Viterbi (双谐波)：这是猎犬的**“双胞胎版”**，它假设信号可能同时以两种频率出现（就像一个人同时发出两个音高）。
5. Weave (编织者)：像一个**“精密织布机”**，它用极其复杂的模板去覆盖所有可能的信号，确保没有死角。

🔍 结果如何？

很遗憾，我们什么也没听到。

• 没有发现新信号：在这 15 个目标中，没有发现任何确凿的“嗡嗡声”。
• 全是噪音：所有看起来像信号的“候选者”，经过仔细检查后，都被证明是仪器本身的杂音或者是宇宙背景噪音的巧合。
  • 比喻：就像你在森林里听到一声响动，以为是老虎，但走近一看，只是风吹树叶的声音。

📉 虽然没有抓到，但我们得到了什么？

虽然没有抓到“真凶”，但这篇论文非常有价值，因为它划定了“禁区”：

1. 设定了“最安静”的标准：
   科学家告诉宇宙：“如果那里有中子星在转，它发出的声音必须比我们设定的这个界限还要轻，否则我们早就听到了。”

   • 对于最近的Vela Jr.（船帆座 X-1 遗迹），我们设定的界限非常严格：如果那里有凸起，它的变形程度必须小于一亿分之一（$10^{-7}$）。这意味着那里的中子星必须圆得像一颗完美的台球，稍微有点歪都不行。

2. 刷新了纪录：
   这次搜索比以前的任何一次都要灵敏。我们的“耳朵”比以前更灵了，能听到的声音更小。这就像是从用普通听诊器升级到了用核磁共振级别的听诊器。

3. 排除了某些理论：
   如果某些理论预测中子星会有很大的变形（像土豆一样），那么我们的结果就告诉那些理论家：“不对，中子星没那么歪。”

🚀 总结

这就好比我们在大海里撒了一张巨大的网，想捞起一条特定的鱼。虽然这次没捞到鱼，但我们证明了：

1. 我们的网比以前更密、更结实了。
2. 如果那条鱼真的存在，它一定比我们要找的还要小、还要隐蔽。
3. 我们排除了很多关于鱼长什么样的错误猜测。

未来的展望：
随着探测器越来越灵敏，数据越来越多，也许在未来的某一天，我们真的能听到那来自宇宙深处的、持续不断的“心跳声”，从而揭开中子星内部神秘的面纱。

---

一句话总结：这是一次极其精密的“宇宙静默”搜索，虽然没抓到“幽灵”（引力波信号），但我们成功地把宇宙中可能存在的“幽灵”藏身之处，压缩到了前所未有的微小范围。

技术摘要

这是一份关于利用 LIGO-Virgo-KAGRA (LVK) 第四次观测运行（O4）早期数据，针对超新星遗迹（SNRs）进行连续引力波（CW）搜索的论文技术总结。

1. 研究背景与问题 (Problem)

• 科学目标：连续引力波（CWs）被认为是由非轴对称、快速旋转的中子星（NS）发出的长时标、准单色信号。尽管瞬态引力波（如双星并合）已被多次探测到，但连续引力波尚未被直接探测到。
• 物理意义：探测 CWs 将提供关于中子星内部结构、物态方程、地壳形变、磁场构型以及能量耗散机制（如引力波辐射或 r 模式振荡）的关键信息。
• 目标源：年轻的超新星遗迹（SNRs）被认为是年轻中子星的理想宿主。由于这些中子星尚未被电磁波（如射电脉冲）探测到（即“隐藏”的中子星或中心致密天体 CCO），通过引力波直接搜索是发现它们的有效手段。
• 挑战：CW 信号极其微弱，且中子星的自转频率和自转减慢率（spin-down）通常未知，且可能受到“自转游荡”（spin wandering）或 glitches 的影响，使得全相干搜索计算成本过高。

2. 方法论 (Methodology)

该研究使用了 LVK 合作组 O4 运行前八个月（O4a，2023 年 5 月 24 日至 2024 年 1 月 16 日）的数据，针对 15 个 邻近的超新星遗迹（包括著名的 Cas A、Vela Jr.、SNR 1987A 等）进行了定向搜索。

研究采用了 5 种互补的搜索管道（Pipelines），以覆盖不同的信号模型和参数空间：

1. Band-Sampled-Data (BSD)：
   • 基于频率 - 霍夫（Frequency-Hough）变换的半相干方法。
   • 搜索范围：10 Hz - 600 Hz。
   • 特点：允许一阶自转减慢（$\dot{f}$）为正或负（考虑自转加速的可能性），忽略二阶项。
2. PyStoch：
   • 基于互相关（Cross-correlation）的辐射计方法，通常用于随机引力波背景，此处用于定向搜索。
   • 搜索范围：20 Hz - 1726 Hz。
   • 特点：无模型（unmodeled）搜索，不假设具体的自转演化模型，通过自适应频率分箱来应对信号展宽。
3. Viterbi 单谐波 (SHV)：
   • 基于隐马尔可夫模型（HMM）的半相干跟踪算法。
   • 特点：专门设计用于跟踪具有长期自转减慢和随机自转游荡的信号。假设信号频率为自转频率 $f_*$ 的两倍（$2f_*$）。
4. Viterbi 双谐波 (DHV)：
   • 同样基于 HMM，但同时跟踪两个频率分量：$f_*$ 和 $2f_*$。
   • 特点：针对可能同时发射基频和谐波的特殊信号模型（如双谐波发射场景），提高了对此类信号的灵敏度。
5. Weave：
   • 基于模板的半相干搜索，使用 F-统计量（F-statistic）。
   • 搜索范围：20 Hz - 475 Hz（避开仪器噪声严重的区域）。
   • 特点：假设功率律自转减慢（$\dot{f} \propto -f^n$，制动指数 $n$ 在 2-7 之间），覆盖较宽的参数空间。

数据处理：

• 使用了 LIGO Hanford (H1) 和 LIGO Livingston (L1) 的校准应变数据。
• 实施了严格的去噪程序，包括剔除已知仪器线（lines）、门控（gating）处理以及针对仪器伪影的剔除。
• 所有候选体均经过后续验证（Follow-up），包括增加相干时间、5-矢量统计量检验等，以排除噪声涨落。

3. 主要贡献 (Key Contributions)

• 最广泛的搜索：这是迄今为止对超新星遗迹进行的最灵敏、最宽频带的定向连续引力波搜索之一，覆盖了 15 个目标源。
• 多管道协同：首次在同一数据集上综合应用了 5 种不同的先进搜索管道，每种管道针对不同的物理假设（如自转游荡、双谐波发射、制动指数变化等），确保了参数空间的最大覆盖。
• 针对 SNR 1987A 的更新：将 SNR 1987A 重新纳入搜索（因其已进入遗迹阶段并显示出致密天体存在的证据），并更新了其距离和年龄参数。
• 对 Ming et al. (2025) 候选体的验证：针对之前 Einstein@Home 项目报道的 G347.3-0.5 的一个候选体（假警报概率约 27%），本研究使用不同的管道配置进行了独立分析，未发现兼容的候选体，进一步确认了该信号可能源于噪声。

4. 研究结果 (Results)

• 无探测发现：在所有 15 个目标源中，未发现任何具有统计显著性的连续引力波信号。所有幸存的候选体在经过后续验证后，均被确认为仪器噪声或统计涨落。
• 上限设定：研究设定了 95% 置信水平的应变振幅（$h_0$）上限。
  • 最严格限制：对于邻近源 Vela Jr. (G266.2-1.2)，在 300 Hz 附近，最严格的应变上限达到 $\sim 4 \times 10^{-26}$。
  • 其他显著结果：
    • BSD 管道在 Vela Jr. 处达到 $h_0 \approx 5.94 \times 10^{-26}$ (250 Hz)。
    • PyStoch 在 Vela Jr. 处达到 $h_0 \approx 9.02 \times 10^{-26}$ (320 Hz)。
    • SHV 管道在 G65.7+1.2 处达到 $h_0 \approx 6.14 \times 10^{-26}$ (164 Hz)。
    • DHV 管道在 $f_* \approx 213$ Hz 处达到 $h_0 \approx 1.82 \times 10^{-25}$。
• 物理参数约束：
  • 椭圆率 ($\epsilon$)：将应变上限转化为中子星的椭圆率限制。对于 Vela Jr.，椭圆率限制低至 $\epsilon \lesssim 10^{-7}$（在高频段），这已经接近或超过了普通中子星的理论上限，表明这些年轻中子星非常接近球形，或者其自转减慢并非完全由引力波主导。
  • r 模式振幅 ($\alpha$)：对于 r 模式发射模型，在 400 Hz 以上，r 模式振幅限制达到 $\alpha \lesssim 10^{-5}$，远小于非线性饱和机制预期的 $\alpha \sim 10^{-3}$。
• 灵敏度提升：与 O3 运行相比，由于观测时间更长（8 个月 vs 6 个月）和探测器灵敏度提高（至少 1.5 倍），应变灵敏度显著提升，部分目标的限制比间接自转减慢极限（spin-down limit）还要严格。

5. 意义与结论 (Significance & Conclusions)

• 排除特定模型：虽然未探测到信号，但这些结果排除了许多年轻中子星具有极大形变（高椭圆率）或强烈 r 模式振荡的模型。
• 技术验证：证明了多管道联合搜索策略在覆盖复杂参数空间（如未知的自转游荡、双谐波发射）方面的有效性。
• 未来展望：随着 LVK 第五、第六次观测运行（O5, O6）的到来，探测器灵敏度将进一步提升，结合更长的数据积累和多信使天文学的进展，探测连续引力波的希望将大大增加。
• 总结：这项工作代表了目前对超新星遗迹中年轻中子星进行连续引力波搜索的最先进水平，为理解致密天体的物理性质设定了新的基准。