Black Hole-Neutron Star Binaries near Neutron Star Disruption Limit in the Mass Regime of Event GW230529

技术摘要：处于中子星瓦解极限质量区间的黑洞-中子星双星系统（以 GW230529 事件的质量范畴为例）

问题陈述
LIGO–Virgo–KAGRA (LVK) 合作组探测到的引力波事件 GW230529，标志着首次可能产生电磁（EM）对应物的黑洞-中子星（BHNS）并合事件。该事件涉及一个质量为 2.5–4.5 $M_\odot$ 的原初黑洞（BH）和一个质量为 1.2–2.0 $M_\odot$ 的次级致密天体，这使其处于一个中子星（NS）可能发生潮汐瓦解而非直接坠入黑洞的质量区间。潮汐瓦解是产生可观测电磁信号（如千新星和紧凑双星伽马射线暴 [cbGRBs]）的先决条件，因为这类信号需要吸积盘的形成以及富中子物质的抛射。

然而，预测特定的 BHNS 并合是否会导致瓦解，严重依赖于解析拟合公式，例如 FNH18 模型（Foucart 等人，2018），该模型用于估算并合后残余盘的质量。这些模型主要基于质量比 $q \gtrsim 3$ 且自旋远离瓦解临界阈值的模拟进行校准。因此，在接近 GW230529 所对应的质量与自旋参数的瓦解极限区域，缺乏高保真度的数值相对论（NR）数据。这一空白限制了准确预测电磁信号普遍程度的能力，也限制了通过电磁观测约束中子星状态方程（EOS）的能力。

方法论
作者利用谱爱因斯坦编码器（SpEC）对 BHNS 并合进行了高分辨率数值模拟。研究重点关注质量比 $q=3$ 且中子星质量为 $M_{NS} = 1.35 M_\odot$ 的情形，这些参数与 GW230529 事件一致。研究采用了两种不同的状态方程来测试不同程度的中子星紧凑度：

1. SFHo： 一种较软的状态方程，其半径 $R_{NS} \approx 11.9$ km。
2. DD2： 一种较硬的状态方程，其半径 $R_{NS} \approx 13.2$ km。

模拟通过改变黑洞自旋幅值（$\chi_{BH}$），来探测紧邻 FNH18 模型预测的瓦解临界自旋（$\chi_{crit}$）的区域。对于 SFHo，$\chi_{crit} \approx 0.037$；对于 DD2，$\chi_{crit} \approx -0.23$。研究包含了自旋略低于及略高于这些临界阈值的模拟情形。

关键方法论特征包括：

• 蒙特卡洛中微子输运： 这标志着首次在 BHNS 并合模拟中应用蒙特卡洛中微子辐射输运。与泄漏模型不同，该方法在位置和动量空间中追踪中微子轨迹，以精确模拟吸收、发射和散射。这对于确定并合后残余物及其流出的成分（电子分数，$Y_e$）至关重要，而这决定了 r-过程核合成的潜力。
• 分辨率研究： 对 SFHo 情形进行了三个分辨率的收敛性测试，以量化数值不确定性。
• 初始条件： 使用 Spells 求解器生成初始数据，并将轨道离心率降低至 $e \lesssim 0.009$。

关键结果

1. 盘质量与解析预测：

   • SFHo EOS： 数值模拟得到的总残余质量通常落在 FNH18 公式预期的建模误差范围内。然而，模拟结果始终显示出比模型预测略高的残余质量。
   • DD2 EOS： 观察到了显著的差异。在自旋接近瓦解极限时，FNH18 模型对残余盘质量的低估程度超过了其声明的误差范围。例如，在 $\chi_{BH} = -0.030$ 时，模型预测残余质量约为 $\sim 0.012 M_\odot$，而模拟结果约为 $\sim 0.049 M_\odot$。
   • 临界自旋阈值： 尽管存在质量差异，但触发瓦解所需的“临界自旋”值在解析模型与模拟之间是一致的（即残余质量从零转变为非零的转折点）。主要的区别在于残余质量曲线在阈值附近的斜率（$dM_{rem}/d\chi_{BH}$）：模拟显示其上升更为陡峭。

2. 低质量盘的出现：

   • 模拟显示，即使在自旋略低于预测的瓦解极限时（例如 SFHo 情形下 $\chi_{BH} = -0.0130$），也会形成非零的吸积盘。虽然这些盘的质量很低（$\sim 10^{-3} M_\odot$），但正如之前的研究（如 Gottlieb 等人）所暗示，它们足以驱动短促 cbGRBs。

3. 质量抛射与成分：

   • 动力学抛射： 并合过程中动力学抛射的质量非常低，通常 $\lesssim 0.004 M_\odot$。对于最低自旋配置，非束缚质量接近数值分辨率的极限，其不确定性可能达到 $\sim 100\%$。
   • 盘风： 研究表明，对于这些系统，通过后期盘风（粘性驱动和磁驱动）抛出的总质量可能将主导初始动力学抛射。
   • 电子分数 ($Y_e$)： 动力学抛射具有高度富中子特性（$Y_e < 0.1$）。然而，由于弱相互作用（中微子加热），流出物和盘物质表现出更宽的 $Y_e$ 范围（最高达 $\sim 0.3$）。DD2 EOS 比 SFHo 产生更富中子的抛射物。

4. 观测意义：

   • 预测这类类似于 GW230529 的并合所产生的任何千新星信号都将是暗弱的，可能由盘风产生的光学辐射主导，而非动力学抛射产生的红外辐射。
   • 用于 r-过程核合成的总抛射质量 $\lesssim 0.02 M_\odot$。

意义与主张
本文主张，虽然 FNH18 解析模型在预测瓦解阈值（临界自旋）的存在性方面仍大致可靠，但它系统性地低估了由此产生的残余盘质量，尤其是在靠近瓦解极限的较硬状态方程（DD2）情形下。

作者强调，他们的结果突显了当前校准在应用于接近瓦解极限的 BHNS 并合时的局限性。具体而言：

• 模型关于临界自旋以下残余质量为“零”的假设并不完全准确；即使在预测极限之下，也会形成低质量盘。
• 这些低质量盘在物理上具有重要意义，因为它们是短伽马射线暴的潜在前身，即使它们不会产生明亮的千新星。
• 研究结果表明，随着更多低质量比、近临界自旋领域模拟数据的出现，可能需要对拟合公式进行适度的重新校准。

该研究并未声称改变了 LVK 对 GW230529 为瓦解型并合概率的评估，指出观察到的差异不太可能从根本上改变 LVK 的分析结果。然而，它为理解此类特定且具有观测关键性的参数空间内的 BHNS 双星行为，提供了更严谨的数值基础。