Numerical simulations of black hole-neutron star mergers with equal and near-equal mass ratios

技术摘要：等质量与近等质量比黑洞-中子星并合的数值模拟

问题与动机
引力波事件 GW230529 的探测表明存在质量比 $q \approx 1/2.6$ 的黑洞-中子星（BHNS）并合，这暗示了比此前预期更为对称的 BHNS 系统可能存在。虽然标准的恒星演化模型预测质量比 $q \lesssim 1/3$ 且通常较小，但致密星团中的动力学形成通道或第二代并合过程可能会产生 $q \sim 1/2$ 甚至 $q \approx 1$ 的系统。理论上预计这些对称系统会产生高质量的吸积盘，从而驱动明亮的电磁瞬变现象（如千新星）。然而，数值相对论（NR）模拟在等质量和近等质量比参数空间内的覆盖非常稀疏。这种数据的缺乏限制了引力波波形模型的准确性以及残余物性质解析拟合公式的精度，可能导致对真实事件的误分类，并错失识别电磁对应体的机会。

方法论
作者使用 BAM 代码对十二组非自旋 BHNS 并合进行了数值相对论模拟。模拟涵盖了质量比 $q \in \{1, 1/2, 1/3\}$ 和中子星（NS）质量 $m_{NS} \in \{0.8, 1.2, 1.4, 2.2\} M_\odot$。采用了两种状态方程（EOS）：较软的 SLy 和较硬的 DD2。模拟采用分段多项式表示法处理 EOS，并使用共形 Z4 公式进行时空演化，使用 Valencia 公式进行流体动力学演化，模拟时长约为九个轨道周期。初始数据由 FUKK 代码构建，并将离心率降低至 $10^{-3}$ 以下。

本研究包括：

1. 引力波波形分析：在多个半径处提取波形，进行二阶外推，并对质心漂移进行修正。将其与 IMRPhenomNSBH 和 SEOBNRv4_ROM_NRTidalv2_NSBH 波形模型进行对比。
2. 残余物性质验证：计算动力学抛射物的质量和吸积盘的质量，以及并合后 10 ms 时残余黑洞（BH）的质量和自旋。将其与文献 [28, 31, 40, 41] 中的解析拟合公式进行对比。
3. 吸积盘动力学：分析具有显著盘质量配置的形态、旋转律和全局密度振荡。
4. 吸积建模：对吸积行为进行分类，并识别其受盘振荡调制的现象。
5. 千新星建模：基于动力学抛射物剖面并加入球对称风成分（盘质量的 30%），使用 POSSIS 辐射传输代码对千新星光变曲线进行建模。

关键结果

• 引力波波形：当前的 BHNS 波形模型（IMRPhenomNSBH 和 SEOBNRv4）在对称质量比区间内与 NR 波形表现出较差的一致性。大多数情况在并合时刻存在约 1 弧度的相位偏差。这些模型往往预测的并合时间较早，且具有比 NR 数据更长的铃宕阶段。值得注意的是，对于亚太阳质量中子星情况（$m_{NS}=0.8 M_\odot$），其潮汐变形性超出了现有模型的有效参数空间；在进行外推时，模型显示出显著的相位和振幅差异。
• 残余物性质：
  • 抛射物与吸积盘质量：用于动力学抛射物质量和吸积盘质量的解析拟合公式与模拟数据表现出良好的一致性。
  • 残余黑洞质量：解析公式能以极高的精度预测残余黑洞质量（误差通常 $< 1\%$）。
  • 残余黑洞自旋：更新后的拟合公式在预测残余自旋方面准确度较低，特别是对于具有大潮汐变形性的系统（DD2 EOS）。对于亚太阳中子星情况，自旋预测偏差接近 50%。作者认为这可能是由于过度拟合了具有较低潮汐变形性和更高质量不对称比的系统。
• 盘结构与动力学：
  • 形态：等质量（$q=1$）盘几乎立即圆化（并合后约 3 ms），而 $q=1/2$ 的盘则长时间受到螺旋密度波和回落物质的影响，在约 8 ms 后才趋于轴对称。
  • 旋转律：比角动量遵循幂律 $j \propto \Omega^{-n}$。指数 $n$ 随半径变化，在最内稳定圆轨道（ISCO）内部表现为恒定角动量（$n=0$），而在外盘则表现为亚开普勒值。这与长期轴对称模拟中常见的固定旋转律假设相矛盾。
  • 振荡：在质量 $> 10^{-2} M_\odot$ 的盘中识别出了全局俘获的 $g$-mode 密度振荡。其频率与残余黑洞的最大前行进进动频率一致。
  • 吸积：区分了三种吸积类型。B 型（$q=1/2$）表现出直接受全局 $g$-mode 振荡调制的吸积率变化。A 型（$q=1/2$）和 C 型（直接坠落）表现出不同的行为。某些系统的 $g$-mode 频率与在伽马射线暴中观察到的高频准周期振荡（QPOs）是兼容的。
• 千新星发射：模拟的千新星主要由盘风抛射物主导。由于较硬的 DD2 EOS 会产生更巨大的吸积盘，因此该类系统产生的千新星比较软的 SLy EOS 系统更亮。对于 200 Mpc 的距离，大多数系统将在四天内被薇拉·鲁宾天文台探测到，并在两天内被 DECam 在 $r$ 和 $i$ 波段探测到。

意义与主张
本文声称填补了对称 BHNS 并合在数值相对论模拟参数空间中的关键空白。其主要意义在于证明了当前的引力波波形模型在等质量及近等质量比情况下精度不足，这可能导致对未来探测事件（如 GW230529）的参数估计产生偏差。本研究验证了残余质量的解析拟合公式，但也强调了其在预测高潮汐变形系统残余自旋方面的局限性。此外，对盘动力学的详细分析揭示了复杂的旋转律和振荡驱动的吸积调制，表明 $g$-mode 可能在伽马射线暴光变曲线上留下可观测的印记。最后，千新星建模表明，对称 BHNS 是当前宽场巡天进行电磁随访的理想目标，与此前主导观测的不对称系统相比，其被探测到的概率更高。