Narrowband searches for continuous gravitational waves from known pulsars in the first two parts of the fourth LIGO--Virgo--KAGRA observing run

本文利用 LIGO-Virgo-KAGRA 第四次观测运行的数据，首次对包含双星系统的 34 颗已知脉冲星进行了涵盖频率二阶导数扩展范围的窄带连续引力波搜索，未发现信号并设定了严格的应变振幅上限，其中蟹状脉冲星的约束达到了其自转减慢极限的 2% 以下。

要点

这是一篇关于寻找“宇宙心跳”中微弱杂音的科学报告。为了让你轻松理解，我们可以把这篇论文想象成一次**“在超级嘈杂的宇宙派对中，寻找特定歌手走调瞬间”**的大行动。

1. 背景：什么是“连续引力波”？

想象一下，宇宙中有一种特殊的恒星叫脉冲星（Pulsar）。它们就像宇宙中的灯塔，旋转得极快，并且非常稳定地发出“滴答、滴答”的电磁波信号（就像完美的节拍器）。

根据爱因斯坦的理论，如果这些脉冲星长得不够圆（比如表面有个小凸起，或者内部有点变形），它们在旋转时就会像洗衣机里甩干不平衡的衣物一样，搅动时空，产生一种持续的、微弱的涟漪，这就是连续引力波（CW）。

• 难点：这种涟漪太微弱了，就像试图在摇滚音乐节的震耳欲聋的噪音中，听到一根针掉在地上的声音。
• 目标：科学家想捕捉这种声音，以此了解脉冲星内部到底是什么做的（是像果冻一样软，还是像钻石一样硬？）。

2. 这次行动做了什么？（O4 观测期）

这篇论文记录了科学家利用LIGO、Virgo 和 KAGRA（三个巨大的“引力波听诊器”）在 2023 年到 2025 年期间（被称为 O4 观测期）进行的一次大搜索。

• 搜索对象：他们盯着34 颗已知的脉冲星。这是历史上针对“窄带搜索”（一种特定的搜索方法）数量最多的一次。
• 新玩法：
  • 以前：科学家假设脉冲星发出的引力波和它发出的无线电波是完美同步的（就像节拍器和鼓点完全重合）。
  • 这次：科学家想得更周全。他们允许引力波和无线电波有一点点**“时间差”或“频率偏差”**。
  • 比喻：以前是只找“完美跟唱”的歌手；这次是找那些可能因为感冒（ glitches，脉冲星突变）或者稍微有点跑调，但依然在唱歌的歌手。他们甚至第一次尝试去听那些在双星系统（像双人舞一样互相绕圈）中跳舞的脉冲星。

3. 他们找到了什么？

很遗憾，没有直接抓到“凶手”（没有发现确凿的引力波信号）。

• 排除干扰：在搜索过程中，他们发现了一些看似像信号的“异常点”。但经过仔细检查，发现这些其实是仪器本身的“打嗝”或环境噪音（比如地球上的震动或电磁干扰），就像在录音时突然有人咳嗽了一声，被误以为是歌手的走调。
• 最终结果：虽然没有抓到信号，但科学家设定了“最严限制”。

4. 没抓到，那有什么意义？（上限与约束）

虽然没有抓到鬼，但这次行动告诉我们要**“鬼”不可能有多大**。

• 比喻：想象你在找一只藏在房间里的老鼠。虽然你没看见它，但你把房间的缝隙都堵死了，并且测量了房间的大小。你可以自信地说：“如果这里有老鼠，它一定比一只仓鼠还小，否则我早就看见了。”
• 具体成就：
  • 对于其中20 颗脉冲星，科学家设定的限制比理论预测的“最大可能值”还要严格。
  • 最牛的案例：著名的蟹状星云脉冲星（Crab Pulsar）。科学家发现，它发出的引力波能量，最多只占它旋转能量损失的 0.04%。
  • 这意味着：蟹状星云几乎是一个完美的球体，或者它的“变形”非常非常小。这极大地缩小了我们对中子星内部结构的猜测范围。

5. 总结：为什么这很重要？

这就好比我们在探索**“宇宙物质的终极形态”**。

• 中子星是宇宙中最稠密的物体，一茶匙的物质就有几亿吨重。
• 通过这种“听音辨位”的搜索，即使没听到声音，我们也知道了这些恒星**“有多硬”、“有多圆”**。
• 这次搜索使用了更聪明的算法（允许一点点偏差），并且第一次尝试了更复杂的“双人舞”（双星系统）场景，为未来真正抓到那个微弱的“宇宙心跳”打下了坚实的基础。

一句话总结：
科学家拿着超级灵敏的耳朵，在宇宙中听了很久，虽然没听到预期的“走调声”，但他们成功证明了那些恒星比想象中更“圆润”，并且排除了很多关于它们内部结构的错误猜想。

技术摘要

这是一份关于利用 LIGO-Virgo-KAGRA (LVK) 第四轮观测运行（O4）前两个阶段（O4a 和 O4b）数据，对已知脉冲星进行连续引力波（CW）窄带搜索的论文技术总结。

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

• 连续引力波 (CW) 的探测挑战：快速旋转且质量分布非轴对称的中子星（脉冲星）是连续引力波的潜在源。然而，CW 信号极其微弱，远低于目前探测到的瞬态引力波信号（如黑洞或中子星并合）。
• 相位失配问题：传统的“全相干”搜索假设引力波（GW）的相位演化与电磁波（EM）观测完全锁定。然而，由于脉冲星计时噪声、自旋突变（Glitches）或引力与电磁辐射机制的不匹配，GW 频率可能偏离 EM 预测值。
• 现有方法的局限：全相干搜索对相位失配非常敏感，一旦失配会导致灵敏度急剧下降。因此，需要一种既能保持高灵敏度，又能容忍小范围频率和自旋参数失配的搜索策略。
• 研究目标：利用 O4 观测数据，对已知脉冲星进行窄带搜索 (Narrowband Search)，以探测那些可能因上述原因导致频率偏移的 CW 信号，并首次将此类搜索扩展到双星系统中的脉冲星。

2. 方法论 (Methodology)

• 数据源：
  • 引力波数据：LIGO 汉福德 (H1) 和利文斯顿 (L1) 探测器在 O4a (2023.05-2024.01) 和 O4b (2024.04-2025.01) 期间的数据。Virgo 和 KAGRA 因在 O4 后期才加入，未包含在此分析中。
  • 电磁数据：来自 Fermi-LAT、Nancay、Jodrell Bank、MeerKAT、Chandra、NICER 等望远镜的脉冲星计时解（Ephemerides），用于提供脉冲星的自转参数（$f, \dot{f}, \ddot{f}$）和轨道参数。
• 搜索算法：
  • 采用 5n-vector 窄带流水线，基于频域匹配滤波（Matched Filtering）。
  • 参数空间扩展：
    • 在频率 ($f$) 和自旋一阶导数 ($\dot{f}$) 空间围绕 EM 预测值进行窄带扫描（通常范围 $\delta = 10^{-3}$）。
    • 创新点：首次探索了自旋二阶导数 ($\ddot{f}$) 空间的参数范围，以应对更复杂的自旋演化。
    • 双星系统：首次在该类窄带搜索中针对双星系统脉冲星，通过修正轨道多普勒效应来搜索信号。
  • 数据处理：
    • 使用 SFDB (Short Fourier Data Base) 框架。
    • 对发生 Glitch 的脉冲星，将数据在 Glitch 发生前后分割成独立段进行分析，排除 Glitch 期间的数据。
    • 对双星系统，根据轨道参数不确定性确定最大相干积分时间，部分长轨道周期脉冲星被排除或分段处理。
• 统计推断：
  • 计算检测统计量 $S$，通过设定误报概率 ($p_{fa} = 1\%$) 确定阈值。
  • 若未发现显著信号，则通过向真实数据中注入模拟信号来设定 95% 置信度 (CL) 的上限 (Upper Limits, ULs)。

3. 关键贡献 (Key Contributions)

1. 最大规模的窄带搜索：这是先进探测器时代以来，针对窄带搜索目标数量最多的一次研究，共分析了 34 颗 已知脉冲星。
2. 双星系统的首次应用：首次将窄带搜索方法应用于双星系统中的脉冲星，扩展了搜索的科学范围。
3. 高阶自旋参数探索：首次在窄带搜索中引入对 $\ddot{f}$（频率二阶导数）的扫描，提高了对自旋演化复杂脉冲星的探测鲁棒性。
4. 修正与验证：在实施新功能时，发现并修复了之前分析（Abac et al. 2025a）中自旋修正的一个 Bug，重新运行后确认之前的上限未受影响，且无新候选体。
5. 灵敏度提升：相比 O1-O3 及 O4a 单独的分析，由于观测时间加倍和探测器升级，灵敏度提升了约 $\sqrt{2}$ 倍。

4. 主要结果 (Results)

• 探测情况：未发现任何连续引力波信号。
• 异常值处理：在 PSR J0117+5914 和 PSR J1826-1334 中发现了统计显著的异常值（Outliers），但经光谱分析确认，这些信号是由 L1 探测器的仪器噪声瞬态（与 3.4 Hz 和 16.3 Hz 校准线相互作用有关）引起的，因此被剔除。
• 上限设定：
  • 为所有 34 颗脉冲星设定了应变振幅 ($h_0$) 的 95% 置信度上限。
  • 突破自旋下降极限 (Spin-down Limit)：在 20 次 分析中（包括对 Glitch 脉冲星的分段分析），设定的上限低于理论上的自旋下降极限。这意味着如果这些脉冲星仅通过引力波辐射损失能量，其辐射功率将小于观测到的自旋下降功率。
  • 最佳约束：对蟹状星云脉冲星 (PSR J0534+2200) 的约束最为严格。其应变上限约为自旋下降极限的 2% ($h_0^{95\%} / h_{sd} \approx 0.02$)。
  • 能量辐射比例：对于蟹状星云脉冲星，这意味着其自旋下降能量中通过引力波通道辐射的比例小于 0.04%。
• 椭圆率限制：
  • 对于年轻脉冲星 ($f \lesssim 100$ Hz)，椭圆率 ($\epsilon$) 上限主要在 $10^{-5}$ 到 $10^{-3}$ 之间。
  • 对于毫秒脉冲星 ($f \gtrsim 100$ Hz)，上限在 $10^{-8}$ 到 $10^{-7}$ 之间。
  • 对于 PSR J0711-6830，上限接近自旋下降极限（比值约 1.15），若考虑中子星转动惯量可能更大（$I_{zz} \approx 3 \times 10^{38} \text{ kg m}^2$），该结果具有潜在的探测意义。

5. 科学意义 (Significance)

• 中子星内部物理：结果排除了中子星存在巨大“山脉”（由地壳应力或强磁场支撑）导致其通过引力波快速损失能量的可能性。特别是对于蟹状星云脉冲星，极低的能量辐射比例限制了其内部非轴对称形变的程度。
• 物态方程约束：虽然目前的上限尚未完全达到某些理论模型（如夸克物质核心或 CFL 相）预测的极端椭圆率，但不断收紧的上限正在逐步排除部分极端物态方程模型。
• 方法学进步：证明了窄带搜索在处理 Glitch 脉冲星和双星系统方面的有效性，为未来更长时间的观测（O5 及以后）提供了成熟的技术框架。
• 多信使天文学：通过结合高精度的电磁计时数据和引力波搜索，进一步加深了对脉冲星自转机制和引力波辐射效率的理解。

总结：该研究代表了连续引力波窄带搜索的最新进展，通过扩大目标样本、引入双星搜索和高阶自旋参数，显著提升了探测能力。虽然未直接探测到信号，但设定的严格上限极大地限制了中子星内部形变和引力波辐射效率的理论空间，特别是为蟹状星云脉冲星提供了目前最严格的约束。