Wide parameter-space O3 search for continuous gravitational waves from unknown neutron stars in binary systems

本文利用先进探测器，针对双星系统中未知中子星发出的连续引力波进行了首次宽参数空间搜索，覆盖频率高于 520 Hz 且轨道周期短于 3 天的范围，虽未探测到信号，但为这类源确立了迄今最严格的信号振幅、椭圆度及 r 模式振幅约束。

要点

想象宇宙中充满了一种恒定、低沉的嗡嗡声，就像一架永不停歇的巨大旋转陀螺发出的声音。这就是科学家所称的连续引力波。这些波是时空结构中的涟漪，由中子星——那些微小却极度致密、城市大小的恒星——因轻微不对称而产生。当它们旋转时，这种晃动会发出稳定的信号，就像灯塔的光束扫过海洋。

然而，寻找这些信号就像在飓风中试图听清一声耳语。大多数时候，我们并不确切知道该往哪里寻找、恒星旋转得多快，或者它是否正围绕着一个伴星（双星系统）起舞。

本文描述了一场由科学家利用先进LIGO探测器组织的、规模宏大的高科技“聆听派对”。以下是他们所做工作的简明解释：

1. 搜索：在宇宙干草堆中寻找一根针

科学家们决定扫描“频率图谱”中一大片未探索的区域。

• 新领域：以往的搜索主要关注较低音调（即较慢的旋转）。而这次团队将搜索推向了更高的音调，高达1,000赫兹。这就像终于调谐到一个此前从未有人检查过的高频电台。
• 双星挑战：许多中子星都有一颗伴星与之相互绕行。这增加了一层复杂性，就像试图听清一位同时在旋转木马上旋转的歌手。轨道运动会改变声音的音调（多普勒效应），使得探测更加困难。此次搜索针对的是那些轨道周期短至0.2天（少于5小时）的“旋转木马上的歌手”。

2. 方法：“筛子”策略

由于宇宙极其浩瀚且数据量巨大，他们无法以完美的专注力聆听每一秒的数据（那将需要超过现有能力的计算资源）。相反，他们采用了一种半相干策略：

• 粗略扫描：他们将数据分割成短片段（每段15分钟）并寻找模式。这就像用粗筛子来捕捉大石块。
• 精细过滤：当他们在粗略扫描中发现一块“石头”（潜在信号）时，便回到该特定位置，利用更长的数据片段以更高的精度进行审视。这就像用放大镜观察那块石头，看它究竟是钻石还是普通石块。

3. 结果：寂静，但这是非常重要的寂静

他们没有发现任何引力波。 没有发现新的中子星。

然而，在科学中，如果“零结果”能告诉我们一些重要信息，它依然是一项胜利。由于他们什么也没找到，现在可以以95%的置信度断言：

• “禁行区”：如果在距离地球100光年范围内存在自转频率高于495赫兹的中子星，那么它们的晃动幅度不足以被我们当前的技术探测到。
• 极限：他们为这些恒星可以有多“凹凸不平”设定了迄今为止最严格的标准。如果一颗恒星距离如此之近且旋转如此之快，其形状必须极其光滑（比煎饼还要平坦）。如果它再稍微凹凸一些，我们本应听到它的声音。

4. 为何这很重要

尽管他们没有发现信号，但这篇论文是一个重要的里程碑，因为：

• 我们突破了上限：他们成功搜索了比以往任何人高出两倍的频率。
• 我们覆盖了新领域：他们探索了以往从未用先进探测器搜索过的轨道周期（恒星相互绕行的速度）。
• 我们证明了技术可行：他们展示了其计算机方法能够处理搜索这些特定、高速、双星系统所带来的巨大复杂性。

简而言之：科学家们调高了他们的宇宙收音机的音量，扫描了一个全新的、高音调的频率范围，寻找成对起舞的恒星，结果一无所获。但通过证明那里什么都没有，他们绘制了一幅极其精确的地图，标明了这些恒星不可能存在的位置，从而为下一代的发现缩小了搜索范围。

技术摘要

技术摘要：针对双星系统中未知中子星连续引力波的大参数空间 O3 搜索

问题与动机
来自非对称自转中子星的连续引力波（CWs）尚未被探测到。虽然针对已知脉冲星的定向搜索和针对孤立未知天体的全天空搜索已经开展，但针对双星系统中未知中子星的搜索面临着独特的计算挑战。源的天体轨道运动引入了额外的参数（轨道周期、投影半长轴、升交点时刻等），显著扩大了参数空间。鉴于双星吸积是产生所需四极矩不对称性的自然机制，且已知自转频率高于 25 Hz 的脉冲星中约有半数处于双星系统中，将搜索覆盖范围扩展至这些系统至关重要。先前的搜索受限于计算约束，往往限制了频率范围和轨道周期。本文旨在探索参数空间中尚未涉足的区域，特别是利用先进探测器针对比以往研究更高的引力波频率和更短的轨道周期进行搜索。

方法论
作者利用先进 LIGO 探测器（H1 和 L1）第三次观测运行（O3）的公开数据进行了半相干搜索。数据已针对 60 Hz 谐波和校准线中的非高斯噪声进行了清洗，并通过门控处理移除了短时 glitches。

• 信号模型：搜索采用了 Jaranowski 等人（1998）的 CW 信号模型，包含内禀振幅参数（$h_0, \cos \iota, \psi, \phi_0$）和相位演化参数（频率 $f_0$、自转减慢 $f_1$、天球位置以及双星轨道参数：$P_{orb}, a_p, t_{asc}, \omega, e$）。
• 参数空间：搜索覆盖的引力波频率范围为 50 Hz 至 1,000 Hz。为了控制计算成本，轨道参数空间 $\{P_{orb}, a_p\}$ 被划分为五个不同的区域（A 至 E）。值得注意的是，E 区域（$P_{orb} \in [0.2, 0.4]$ 天，$a_p \in [0.03, 0.12]$ 光秒）在此项工作之前完全未被探索。搜索假设轨道偏心率较小（$e \leq 0.1$），且自转减慢限制与观测到的双星中子星群体一致。
• 搜索流程：分析采用了 BinarySkyHouF 半相干流程。
  • 初始阶段：数据被分割为 $N_{seg} = 31,675$ 个时长 $T_{seg} = 900$ 秒的片段。利用模板统计量的线性组合计算检测统计量（$2\hat{F}_{AB}$），这些模板的轨道参数为零，但通过加权以追踪包含非零轨道参数的时频模式。网格间距经过优化，使失配度保持在 2% 以下。
  • 后续跟进：根据显著性选择候选体，并对其进行 $T_{seg} = 2,700$ 秒的相干跟进。采用嵌套采样算法计算检测统计量和稳健性统计量（$\log_{10} \hat{B}_{1;S/GLtL}$），以区分天体物理信号与仪器谱线及噪声。
• 验证：将模拟测试信号注入数据中以校准检测阈值并估算上限。候选体针对硬件注入和非高斯干扰进行了审查。

主要贡献

1. 扩展频率范围：这是首次覆盖 520 Hz 以上直至 1,000 Hz 引力波频率的搜索。这一进展具有重要意义，因为在“回收”情景中，中子星通过吸积获得角动量，可能达到极高的自转频率。
2. 短轨道周期：这是首次利用先进探测器探索低于 3 天轨道周期的搜索，特别是针对 E 区域（$P_{orb} \approx 0.2-0.4$ 天）。
3. 参数空间广度：该搜索覆盖的参数空间比任何先前利用先进 LIGO 数据进行的调查大约四个数量级。
4. 计算策略：作者采用了“广度优先”的方法，与灵敏度优化的搜索相比，使用了更粗糙的模板网格和更短的相干基线，从而在固定的计算预算内覆盖了广阔的、此前未探索的区域。

结果

• 探测：未探测到任何天体物理连续引力波信号。
• 候选体：从初始约 127,500 个候选体中，有 15 个通过了跟进审查。然而，其中 12 个被确认为与硬件注入相关（特定频率处的孤立中子星信号），其余 3 个归因于非高斯干扰。没有一个是可辩护的 CW 信号。
• 上限：作者建立了迄今为止对这类信号振幅最严格的约束。
  • 灵敏度深度：平均灵敏度深度确定为 $\langle D_{95\%}^{low} \rangle = 18.1 \pm 1.1$ Hz$^{-1/2}$（50–500 Hz）和 $\langle D_{95\%}^{high} \rangle = 17.1 \pm 0.9$ Hz$^{-1/2}$（500–1000 Hz）。
  • 椭圆率约束：对于距离在 100 pc 以内且自转快于约 495 Hz 的中子星，在 95% 置信度下排除了大于 $5.2 \times 10^{-8}$ 的椭圆率（假设 $I_{zz} = 10^{38}$ kg m$^2$）。
  • r-模约束：对于相同的距离和自转速率（> 约 740 Hz），排除了大于 $1.5 \times 10^{-6}$ 的 r-模振幅。

意义
该论文声称，这项工作代表了迄今为止针对双星系统中未知中子星连续引力波最广泛的搜索。通过将频率极限推至 1,000 Hz 并探索低至 0.2 天的轨道周期，该搜索解决了参数空间中的关键空白，而这些区域正是吸积回收的中子星预期存在的所在。虽然未发现信号，但结果提供了迄今为止对这类源椭圆率和 r-模振幅最严格的约束。作者将此次搜索定位为调查新的、更敏感数据集的“蓝图”，展示了尽管存在计算扩展（$\propto f_0^5$），但在具有挑战性的高频区域中，最先进的技术所具备的能力。这项工作承认参数空间的广阔区域仍未被探索，但为未来更灵敏的搜索建立了框架。