在寂静的深空之中，两个巨大的天体相互旋绕，轨道不断收紧，直到发生剧烈的碰撞合并。当这些天体是黑洞时，它们的碰撞会向时空本身发送涟漪，即被称为引力波的波动。地球上的探测器极其灵敏，足以测量比单个原子还要微小的变化，它们时刻倾听着这些涟漪。2019年5月探测到的其中一次事件，与以往所见的任何事件都截然不同。它来自两个极其沉重的黑洞，它们的总质量接近太阳质量的150倍，并最终合并成了一个规模更大的单一黑洞。这次事件如此宏大，以至于它向地球发出的信号极其短暂，仅持续了不到一秒钟，探测器在信号消失前仅捕捉到了寥寥数个波形。由于信号如此短暂，科学家们面临着一个难题：这短短的数据包含了许多种可能的解释，而用于解读数据的标准工具却难以将它们区分开来。
一组研究人员现在回到了这个神秘的事件，试图利用一套旨在寻找特定且常被忽视的特征——轨道离心率（orbital eccentricity）——的新工具来解开谜团。在完美的平滑轨道中，两个物体绕彼此做圆周运动。然而，如果轨道被拉伸成椭圆形，则被称为离心轨道。几十年来，科学家一直假设在寂静孤立的空间中合并的黑洞拥有圆形的轨道。但如果黑洞形成于拥挤的恒星集群中，它们可能会在碰撞前被抛入椭圆形的路径。研究人员使用能够解释这些椭圆轨道的波形模型，以及由黑洞自旋引起的摆动运动，对2019年这一事件的数据进行了重新分析。他们发现，数据有力地支持了黑洞确实是在高度离心的椭圆轨道上运动的观点。分析表明，在信号开始的那一刻，轨道的拉伸程度为0.81，这一数值表明其是一个非常明显的椭圆形状而非圆形。
这项研究表明，此前对该事件的困惑是因为科学家使用的标准模型不允许观察这些高度椭圆的可能性。当研究人员将搜索范围限制在低离心率水平时，模型无法找到轨道的真实形状，导致数据看起来平坦且缺乏信息量。然而，一旦允许模型探索全范围的椭圆形状，数据便清晰地指向了一个高度离心的系统。这一发现还解决了第二个谜团：早期的研究曾暗示这些黑洞在向内螺旋的过程中在剧烈地摆动。新的分析表明，这种表象上的摆动实际上是一种错觉。之前的模型试图通过虚构一种摆动来解释椭圆形状的轨道。当椭圆形状被纳入计算后，对摆动的需求便消失了，数据也显示没有显著的自旋诱导进动（spin-induced precession）证据。
研究人员还检查了两颗黑洞的质量。他们的分析确认，较大的黑洞几乎肯定位于一个理论上的“质量间隙”内，即恒星被认为无法自然形成黑洞的一个重量区间。主黑洞落在该间隙内的概率至少为97%。较小黑洞的状态则不那么确定，取决于使用哪种具体的模型，但最稳健的模型表明，它也极有可能处于这一异常的质量范围内。为了确保他们的发现并非仅仅是数学计算的结果，团队进行了一系列测试。他们创建了模拟真实事件的伪造信号，将其注入分析软件，以观察模型是否能正确识别来源。当他们向软件输入一个具有高度椭圆轨道但没有摆动的伪造信号时，那些忽略椭圆形状的标准模型错误地发明了一种摆动来解释数据。相反，当他们输入一个具有真实摆动但轨道为圆形的信号时，模型正确识别了摆动，并未虚构出椭圆形状。这证实了模型运行正常，且真实的事件确实是高度离心的。
最后，团队将他们最先进的模型与已知是高度离心的黑洞碰撞计算机模拟进行了对比测试。他们发现，尽管这些模拟实验比该模型之前测试过的任何情况都要极端，但模型仍能成功重现模拟中的波形形状。这使他们对解释真实事件的结论充满信心。该研究得出结论，2019年的这一事件是两个处于高度拉伸的椭圆轨道中的黑洞的碰撞，很可能形成于拥挤的恒星环境中。这一发现改变了我们对这些巨大天体行为的理解，并解释了为何早期尝试描述该事件的方法会导致相互矛盾的结果。通过允许数据在不受强制要求呈圆形的情况下自我表达，研究人员揭示了一个更清晰的画面，描绘出了引力波天文学历史上最宏大、最剧烈的事件之一。
技术摘要：迈向对 GW190521 特性的全面理解
问题陈述
由 LIGO–Virgo–KAGRA (LVK) 合作组织探测到的引力波事件 GW190521，由于其极高的总质量，为参数估计带来了独特的挑战。该信号在探测器敏感频率范围内仅包含极少数可观测周期，使得推断过程对波形建模假设高度敏感。先前的分析得出了相互矛盾的结论，关于该事件的物理特性，特别是轨道离心率和自旋进动。虽然一些研究表明存在高度离心的源，但另一些研究则发现没有离心率证据，或倾向于自旋进动的类圆轨道配置。这些分歧归因于离心率与自旋进动之间的简并性，以及各种模型中使用的有限离心率先验范围。本文旨在通过使用包含离心率和自旋进动的尖端波形模型进行全面的重新分析，来解决这些争议。
方法论
作者使用 RIFT 参数估计算法对 GW190521 进行了贝叶斯重新分析。分析使用了来自 LIGO Hanford、LIGO Livingston 和 Virgo 天文台的 4 秒数据，在 20–224 Hz 的频率范围内评估似然度。
关键方法论组成部分包括：
- 波形模型： 研究采用了多种波形族来测试物理假设：
- 离心模型： TEOBResumS（包含通用自旋取向）、IMRPhenomTEHM 以及 SEOBNRv5EHM（限制为对齐自旋）。
- 自旋进动模型： TEOBResumS、IMRPhenomTPHM、IMRPhenomXPHM、NRSur7dq4 以及 SEOBNRv5PHM。
- 组合模型： 使用 TEOBResumS 来测试结合离心率和自旋进动的配置 (-E+P)。
- 先验： 对探测器坐标系下的质量、天球位置和亮度距离应用了均匀先验。对于离心率（e5 在 5 Hz 处），通常在 [0.0,0.9] 上设置均匀先验，尽管由于 SEOBNRv5EHM 在较高离心率下出现波形生成失败，其范围被限制在 [0.0,0.4]。
- 验证研究：
- 零噪声注入： 作者注入了基于 GW190 极值似然参数（既高度离心又高度自旋进动）的合成信号，以测试模型的恢复能力并识别系统误差来源。
- 数值相对论 (NR) 恢复： 注入了高度离心的 NR 波形（RIT eBBH 1637 和 1639），并使用 TEOBResumS-E+P 进行恢复，以验证模型在 GW190521 分析所青睐的高离心率机制下的性能。
主要贡献与结果
高离心率证据：
使用具有宽离心率先验（e5∈[0.0,0.9]）的 TEOBResumS 和 IMRPhenomTEHM 进行的分析，稳健地推断出一个高度离心的源，其 e5≈0.81−0.01+0.05。相比之下，SEOBNRv5EHM 和受限范围分析（e5≤0.4）产生的是无信息或平坦的后验分布。作者证明，先前研究中缺乏离心率证据的原因在于，将先验限制在高于高似然区域（e5>0.8）的范围之下，阻碍了后验分布获取真实的解。全范围分析通过对数贝叶斯因子 lnB≈4.9 被优于受限范围分析。
离心率–自旋进动简并：
当分析中包含离心率时（TEOBResumS-E+P），推断的有效进动自旋参数 (χp) 与类圆轨道自旋进动模型相比显著向较低值偏移。虽然类圆轨道模型倾向于中等到大的 χp (∼0.6)，但离心分析得到的 χp≈0.31。这表明，在类圆轨道分析中观察到的表观自旋进动，部分是由于模型试图解释未建模的离心特征而驱动的。数据显示，在妥善考虑离心率后，没有明显的自旋进动迹象。
质量比与 PISN 间隙：
所有波形模型都倾向于一个质量比相近的二元系统 (q≈0.7)。然而，某些模型（例如 IMRPhenomXPHM）在高度不对称质量处（q∼0.2）表现出次级模态。关于对撞不稳定性超新星 (PISN) 质量间隙（∼60–130M⊙）：
- 所有模型都一致认为，原初黑洞落在间隙内，且概率极高（≥97%）。
- 次级黑洞的概率取决于模型和离心率先验。全离心率范围分析将次级黑洞落在间隙内的概率提高到了 ∼92%，而受限范围或特定模型产生的概率较低。
波形系统误差的起源：
零噪声注入研究表明，在真实的 GW190521 分析中观察到的波形模型系统误差（质量、自旋和离心率后验的差异）主要源于类圆轨道模型对未建模离心率的响应。当使用类圆轨道模型分析高度离心的信号时，它们会错误地推断出显著的自旋进动，并表现出依赖于模型的质量比。相反，当使用真正具有自旋进动的类圆轨道信号时，不同模型会得出一致的结果。
模型验证：
使用 TEOBResumS-E+P 恢复高度离心的 NR 波形，证明了该模型能够重现整体波形形态并恢复注入的质量和自旋。然而，离心率的恢复并不总是唯一的；其中一个注入被正确识别为高度离心，而另一个则产生了包含零离心率支持的多模态后验。这突显了在短信号中唯一识别高离心率的难度。
意义与主张
本文声称，要完整理解 GW190521，需要明确包含轨道离心率的波形模型。作者得出结论，GW190521 最可能是一个高度离心的双黑洞并合事件。这种解释通过展示非离心结果源于对离心参数空间探索不足，从而解决了之前的争议。此外，本研究确立了此前推断的自旋进动很可能是离心率–自旋进动简并性的产物。
作者强调，限制波形模型的物理内容（例如，假设类圆轨道）会导致对短、高质量引力波信号产生定性的不同解释。虽然研究结果强烈支持离心解释，但作者指出，对于此类高离心率现象的稀有性，以及需要将贝叶斯证据与天体物理先验几率相结合，必须对天体物理起源的推断保持审慎。这项工作强调了扩大波形建模参数空间的必要性，以避免在未来的引力波天文学中产生系统偏差。
每周获取最佳 general relativity 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题，再试一次？
无垃圾邮件，随时退订。