Forming Double Neutron Stars using Detailed Binary Evolution Models with POSYDON: Comparison to the Galactic Systems

本研究利用 POSYDON 双星群体合成代码揭示了银河系双中子星通过两种截然不同的共用包亚通道（Case B 和 Case C）形成，这要求特定的核心定义或高抛射效率以及低超新星爆发踢力，从而解释了观测到的轨道周期分裂并约束了这些系统的详细演化物理过程。

要点

大局观：一场宇宙侦探故事

想象一下，宇宙是一个巨大的、混乱的舞池。大多数时候，恒星都在独自起舞。但有时，两颗质量巨大的恒星会同时诞生，并手牵着手。在数百万年的时间里，它们会经历一场疯狂而戏剧性的舞蹈，最终有两种结局：要么永远漂离彼此，要么撞在一起并合并。

这篇论文是一场关于双中子星 (DNS) 的侦探故事。它们是宇宙舞池中的“重量级选手”——两颗密度极高的死寂恒星（中子星）正绕着彼此运行。在我们所在的银河系中，我们已知有 25 对这样的组合。作者们使用了一个名为 POSYDON 的超级计算机模拟程序，来弄清楚这些星对究竟是如何形成的，以及为什么它们呈现出现在的样子。

主要发现：通往同一目的地的两条不同路径

这篇论文中最令人惊讶的发现是，制造双中子星并不只有一种方法。作者发现，这些星对是通过**两个截然不同的“子通道”**形成的，就像是通往同一座城市的两条不同高速公路。

区别在于当这场大戏开始时，第二颗恒星的“年龄”和“成熟度”：

1. “年轻”路径（氦核通道）：

   • 场景： 想象一颗年轻的恒星，它还没有多少成长的机会。它的核心仍然是由氦组成的（就像一个还没读完高中的青少年）。
   • 戏剧性时刻： 当这颗年轻的恒星试图膨胀时，它会离它的伴侣（第一颗中子星）太近。它们会纠缠在一个巨大的、混乱的气体云中，被称为共同包层 (Common Envelope)。
   • 结果： 为了从这种纠缠中生存下来，系统需要大量的能量来踢开这团气体云。如果成功了，它们最终会形成一个非常紧密、快速运行的轨道。这些就是那些最终会撞在一起并产生引力波的星对（即“合并型”）。

2. “年老”路径（碳氧核通道）：

   • 场景： 现在想象一颗更年长、更成熟的恒星。它已经消耗掉了大部分燃料，核心由碳和氧组成（就像一位已经结束职业生涯的退休人员）。
   • 戏剧性时刻： 这颗年老的恒星也会与它的伴侣一起纠缠在气体云中。
   • 结果： 因为这颗年老恒星的气体云比较松散，更容易被推开，所以系统能更容易地从纠缠中幸存。然而，它们最终会处于一个更宽、更慢的轨道上。这些星对之间的距离太远，以至于在宇宙的寿命内都不会发生碰撞。

类比： 把这想象成两对情侣试图在拥挤的房间里平息一场争吵。

• 年轻的情侣（氦核）是在一个狭窄、拥挤的空间里吵架。要推开人群需要付出巨大的努力（能量），但如果他们成功了，最后会抱得非常紧。
• 年老的情侣（碳氧核）是在一个宽敞的房间里吵架。推开人群很容易，但最后他们会站得很远。

“踢”的问题：为什么舞池如此湿滑

当一颗恒星死亡时，它会作为超新星爆发。这种爆炸通常会给新的中子星一个“踢”（kick），就像台球中的母球被撞击一样。

• 旧理论： 科学家以前认为这些“踢”是非常巨大且随机的（就像用大锤猛击母球）。
• 新发现： 作者发现，如果这些“踢”真的那么大，大多数星对都会飞散，永远不会成为双中子星。
• 解决方案： 为了符合我们在天空中实际观察到的现象，“踢”必须非常温柔（就像轻轻的一拍）。论文表明，由于第二颗恒星在爆炸前已经被剥离了大部分重量，因此它获得的“踢”非常轻微。这使得这对星对能够保持在一起。

第一颗恒星的“回收利用”

在第二颗恒星爆炸之前，第一颗中子星（较老的那颗）会经历一次“大改造”。随着第二颗恒星脱落气体，第一颗恒星会吸收其中的一部分。这就像是一个宇宙级的饮食计划，实际上让这颗恒星旋转得更快、更亮，使其变成一颗“回收脉冲星”。

论文显示，“年轻路径”（氦核）为第一颗恒星提供了更好的“饮食”，使其旋转得更快。这解释了为什么我们在天空中看到的合并型星对（亚群 i）旋转得很快，而由于非合并型星对（亚群 ii）旋转得较慢。

那些“怪胎”呢？

论文指出，有一小组由三颗双中子星组成的特殊群体不属于上述任何一条路径。它们的轨道非常古怪，难以解释。作者承认他们的模型无法轻易解释这三个天体，这表明它们可能是通过完全不同的方式形成的（例如，是在拥挤的星团中通过碰撞形成的，而不是作为一对星对诞生的）。

结果总结

• 分裂： 双中子星通过两条主要路径形成：一条针对“年轻”捐赠者（导致合并型星对），另一条针对“年老”捐赠者（导致非合并型星对）。
• 能量： 要让“年轻路径”奏效，宇宙需要非常高效地抛弃气体云。
• “踢”： 创造这些恒星的爆炸必须是温和的，而不是暴力的，才能让这对星对保持在一起。
• 合并： 那些最终会碰撞并产生引力波（如著名的 GW170817 事件）的星对，几乎肯定来自“年轻路径”。

本论文没有说明的内容

• 它并没有预测下一次合并何时发生。
• 它并不声称解释了宇宙中每一个奇特的恒星（特别是上面提到的那三个“怪胎”）。
• 它没有包含一个关于恒星如何“自转减速”或逐渐暗淡的模型，这限制了它目前与现实世界观测结果的完美匹配程度。

简而言之，这篇论文通过详细的计算机模拟表明，宇宙有两种截然不同的制造双中子星的配方，而其中只有一种配方能产生那些最终会碰撞并震动时空结构的星对。

技术摘要

技术摘要：利用基于 POSYDON 的详细双星演化模型形成双中子星

问题陈述
双中子星（DNSs）是理解大质量双星演化、超新星机制以及引力波源的关键实验室。虽然银河系中已识别出二十多个双中子星，但其形成通道仍然复杂且存在争议。以往的群体合成研究主要依赖于快速双星群体合成（BPS）代码，这些代码利用拟合公式或单星演化轨迹库。这些方法往往缺乏对双星相互作用（特别是在共同包层（CE）阶段和核心坍缩超新星（SN）阶段）的自洽处理。此外，观测到的银河系双中子星群体在轨道周期和离心率上表现出明显的双分性，标准模型在没有精细调优的情况下难以重现这一现象。本文旨在通过更严谨、自洽的框架，解决约束导致双中子星形成的详细演化阶段的需求。

方法论
作者使用了 POSYDON（基于详细双星演化模拟的群体合成），这是一种最先进的双星群体合成代码，它整合了基于一维代码 MESA 的详细单星和双星演化模型。与快速 BPS 代码不同，POSYDON 采用预计算的恒星演化网格和即时计算来处理复杂的物理过程，包括磁制动、引力波辐射、共同包层（CE）演化和超新星反冲（kicks）。

该研究模拟了 $2 \times 10^6$ 个双星群体，涵盖 25 种不同的模型变体（汇总于表 2）。操纵的关键变量包括：

1. 共同包层（CE）处理： 变化包括 CE 抛射效率（$\alpha_{CE}$）、核心-包层边界的定义（氢质量分数分别为 1%、10% 和 30%），以及 post-CE 质量损失的处理方式（稳定质量转移 vs. 风损失 vs. 单相）。
2. 超新星（SN）模型： 不同的核心坍缩结果预设（Sukhbold 等人 2016；Patton & Sukhbold 2020；Fryer 等人 2012），区分铁核坍缩（CCSN）和电子捕获超新星（ECSN）。
3. 诞生反冲（Natal Kick）预设： 将标准的 Hobbs-Maxwellian 分布（$\sigma \approx 265$ km/s）与“基于物理机制”的模型进行比较：
   • AsymEj： 源自非对称喷发机制（Janka 2017）的反冲，其大小随喷发质量缩放。
   • Linear： 反冲随喷发质量线性缩放（Richards 等人 2023）。
   • Low Maxwellian： 降低的色散模型（$\sigma \approx 20$ km/s）。

模型假设为太阳金属丰度及恒定的恒星形成史。生成的合成群体在轨道周期（$P_{orb}$）和离心率（$e$）分布以及合并率方面，与观测到的 25 个银河系双中子星样本（表 1）进行定性比较。

主要贡献与结果

• 形成通道的双分性： 研究确定了 CE 形成通道的一个基本双分性，该双分性由捐赠者恒星在与第一颗诞生中子星开始洛希瓣溢流（RLO）时的演化状态决定：

  • 氦核通道（Case B RLO）： 捐赠者在拥有氦核时启动质量转移。这些系统通常具有较短的轨道周期和较高的包层结合能。
  • 碳/氧核通道（Case C RLO）： 捐赠者在 RLO 之前进一步演化以发育出碳-氧核。这些系统拥有结合较松的包层，并导致较宽的轨道周期。
  • 意义： 这种双分性自然地解释了银河系双中子星群体在短周期合并系统（子群体 i）与宽轨道非合并系统（子群体 ii）之间的分裂。

• 对共同包层（CE）物理机制的约束： 氦核通道（产生合并双中子星的通道）的形成对 CE 效率高度敏感。模型显示，该通道仅在以下条件下存活：

  • 使用较宽松的核心定义（H-分数 $\simeq 30\%$），或者
  • 采用较高的 CE 抛射效率（$\alpha_{CE} \gtrsim 1.2$）。
    若无上述条件，氦核系统将在 CE 阶段内发生合并。C/O 核通道则更为稳健，即使在较低的 $\alpha_{CE}$ 值下也能存活。

• 诞生反冲约束： 标准的 Hobbs-Maxwellian 反冲分布（$\sigma \approx 265$ km/s）被排除在外，因为它会在第一次超新星爆发期间破坏大多数双星，导致合并率比观测值低几个数量级。观测到的银河系双中子星群体最符合低速度反冲预设（$\lesssim 50$ km/s）。“基于物理机制”的 AsymEj 和 Linear 模型自然地为被剥离的前身星产生较低的反冲，成功重现了轨道分布和局部合并率密度（$O(10) - O(100)$ Gpc$^{-3}$ yr$^{-1}$）。

• 脉冲星回收与质量吸积： 研究发现，脉冲星回收（使第一颗诞生中子星旋转加速）主要发生在与氦星伴星的稳定 Case BB 质量转移阶段，而非 CE 阶段。通过氦核通道形成的系统经历了更高效的回收，这与观测到的较短自转周期（属于子群体 i）一致。

• 超新星类型与前身星：

  • 第一次超新星主要是铁核坍缩（CCSN）。
  • 第二次超新星通常涉及超剥离（ultra-stripped）前身星。大约 60% 的合并双中子星前身星在第二次超新星爆发时正处于洛希瓣溢流状态。
  • 合并双中子星诞生于前身星中值前超新星包层质量约为 $\sim 0.2 M_\odot$ 的系统。这表明很大一部分合并双中子星可能被观测为 Ib 型或 Ic 型超新星。

• 合并率与引力波源： 模型预测，只有通过氦核通道形成的双中子星会在哈勃时间内发生合并。因此，推断引力波事件如 GW170817 和 GW190425 源自该通道。基于物理反冲模型的预测局部合并率与 LIGO-Virgo-Kagra 以及银河系双中子星样本的约束相一致。

意义与主张
本文声称，详细的演化模型对于解决双中子星的形成历史至关重要。通过超越拟合公式并转向自洽的网格模型，作者证明了观测到的银河系双中子星的多样性是捐赠者恒星在 CE 阶段开始时的演化状态的自然结果。

作者谦虚地指出，由于目前 POSYDON 框架内缺乏脉冲星演化模型（用于考虑选择效应和自转减慢），其结论是基于定性比较得出的。他们也承认其结果是针对太阳金属丰度的。然而，该研究确立了以下几点：

1. 形成通道的双分性（氦核 vs. C/O 核）是双星演化的一个稳健特征，解释了观测到的轨道周期分裂。
2. 高效的 CE 抛射（通过宽松的核心定义或高 $\alpha_{CE}$）是形成合并群体所必需的。
3. 低诞生反冲是双中子星双星存活的必要条件，支持了近期从高速度反冲分布转向物理机制模型的理论转变。
4. 不稳定的 Case BB 质量转移在这些详细模型中不会产生合并双中子星，这对有关超快合并系统形成的假设提出了挑战。

这项工作为未来在集成脉冲星演化模型后进行定量贝叶斯分析奠定了基础。