[论文解读] General Relativistic Numerical Simulation on Coalescing Binary Neutron Stars and Gauge-Invariant Gravitational Wave Extraction

Mari Kawamura, K. Oohara+1 more|arXiv (Cornell University)|Jun 24, 2003

Pulsars and Gravitational Waves Research被引用 8

一句话总结

本文采用最大切片法与规范不变的引力波提取方法，对双中子星并合过程进行了三维广义相对论数值模拟。结果表明，合并后约25 mHz处可能存在准正规模振荡，但由于数值边界处的波反射，精度受到限制；结果与先前工作定性一致，但因初始条件非平衡态导致辐射能量较低。

ABSTRACT

We are developing 3 dimensional simulation codes for coalescing binary neutron stars. A code using the maximal slicing condition is obtained. To evaluate the gravitational radiation, we implemented a gauge-invariant wave extraction and compared the wave forms with the metric tensors at the wave zone. The energy spectrum of the waves was also evaluated to investigate the possibility that the excitation of the quasi-normal modes of the black hole, which may be formed after the merger of two stars, can be observed.

研究动机与目标

• 开发一种基于最大切片法的三维广义相对论模拟代码，用于双中子星并合，以避免共形切片法中出现的不稳定性。
• 实现规范不变的引力波提取，确保波形物理结果与坐标选择无关。
• 评估引力波能量谱，并研究并合后形成的黑洞是否可能激发准正规模。
• 与现有模型（特别是 Shibata 和 Uryu 的工作）进行比较，以验证代码，并评估辐射能量与角动量。

提出的方法

• 采用广义相对论的(3+1)形式，将爱因斯坦方程分解为12个演化方程和4个约束方程。
• 使用最大切片法（K=0）以提高长期数值稳定性，替代先前工作中不稳定的共形切片法。
• 通过扩展 Abrahams 等人的方法，实现规范不变的引力波提取，适用于非轴对称、偶极与奇极扰动。
• 通过求解 lapse 函数 α、位移矢量 βi 和共形因子 φ 的椭圆方程，以维持约束满足与坐标控制。
• 采用伪最小畸变条件以减少坐标漂移，时间上监测到 ∼2% 的 ∂j(γ̃ij) 误差。
• 使用 475×475×238 的网格，内存占用 80 GB，每次运行约需 100 个 CPU 小时，主要计算瓶颈为椭圆方程的求解。

实验结果

研究问题

• RQ1 采用最大切片法的三维广义相对论代码能否准确模拟双中子星并合过程，并提取规范不变的引力波形？
• RQ2 双中子星并合是否能激发所形成黑洞的准正规模？这一现象是否能在引力波谱中被观测到？
• RQ3 与准平衡模型相比，初始非平衡条件如何影响辐射能量与角动量？
• RQ4 数值边界反射在多大程度上扭曲了引力波信号？如何减轻这种影响？

主要发现

• 规范不变的波提取方法成功隔离了物理引力波信号，外部真空区域中的规范依赖模态被有效抑制。
• 在引力波能量谱中观测到约25 mHz处的峰值，与质量为2.8 M⊙ 的史瓦西黑洞的 l=2 基态准正规模一致。
• 辐射能量与角动量低于 Shibata 和 Uryu 的模型 H-1，主要原因是非平衡初始数据导致并合速度加快。
• ADM 质量在 1.7 ms 内保持在 10% 以内，但并合后因网格分辨率不足与黑洞形成而逐渐失守。
• 引力波在数值边界处的反射在约 2.5 ms 后变得显著，且由于边界距离过近，反射波持续增强。
• 代码性能受限于求解椭圆方程（α, βi, φ）的 CPU 时间，未来改进需依赖更快的求解器与更细的网格。