黑洞碰撞计算器

计算双黑洞系统的并合时间、引力波能量和最终黑洞性质。

输入两个黑洞的质量、初始轨道分离、偏心率和轨道倾角,即可使用 Peters 的广义相对论公式计算螺旋并合参数。

黑洞碰撞计算器
计算双黑洞系统的并合时间、引力波能量和最终黑洞性质。

关于黑洞碰撞计算器

当两个黑洞形成引力束缚的双星系统时,它们会通过辐射引力波逐渐损失能量和角动量,彼此螺旋靠近。这个过程称为并合前螺旋演化,最早由 Peter Peters 于 1964 年基于爱因斯坦广义相对论推导出数学框架。黑洞碰撞计算器采用 Peters 公式,估算双黑洞并合所需时间、辐射能量以及最终并合体的外观。 引力波天文学中最重要的派生量是啁啾质量,定义为 M_chirp = (m₁ m₂)^(3/5) / (m₁+m₂)^(1/5)。啁啾质量决定引力波频率上升的速率——也就是信号从低频扫向高频的典型特征。再结合对称质量比 η = m₁m₂/(m₁+m₂)²,啁啾质量就编码了计算一阶并合前螺旋动力学所需的全部信息。 Peters(1964)给出的圆轨道并合时间公式为 T = (5/256) × c⁵ a⁴ / (G³ m₁ m₂ M_total),其中 a 是初始半长轴,G 是引力常数,c 是光速。对于偏心轨道,并合时间会大约乘以 (1−e²)^(7/2) 的修正因子,这意味着高偏心双星比相同初始分离下的圆轨道更有效地辐射能量,也会更快并合。对于 e ≲ 0.6,这一近似会越来越准确。 并合前螺旋阶段辐射的引力波能量估计约为约化质量能量(μc²)的 5%,这与相近质量双星的数值相对论模拟相符。剩余质量形成最终并合黑洞,其史瓦西半径为 r_s = 2 G M_final / c²。最内稳态圆轨道(ISCO)的峰值引力波频率为 f_peak = c³ / (π × 6√6 × G × M_total),它标志着从并合前螺旋到急坠与振铃的过渡——也就是并合最响亮、最容易探测的时刻。 首个直接探测到的引力波事件 GW150914 由 LIGO 于 2015 年 9 月 14 日观测到。它来自两个质量约为 36 和 29 个太阳质量的黑洞,在 13 亿光年外完成并合。该事件在不到一秒内将约 3 个太阳质量的能量以引力波形式辐射出去,短暂地在引力波光度上压过了整个可观测宇宙。此后,LIGO–Virgo–KAGRA 合作组已探测到 90 多个双星并合事件,使引力波天文学成为一门精密科学。

黑洞碰撞示例

下表展示了具有代表性的双黑洞系统及其关键并合参数。

系统参数关键结果事件 / 背景
m₁=36 M☉, m₂=29 M☉, a=10,000,000 km, e=0T_merge ≈ 94.4 Myr, M_chirp ≈ 28.1 M☉, f_peak ≈ 67.6 Hz类似 GW150914(LIGO,2015)
m₁=1000 M☉, m₂=800 M☉, a=100,000,000 km, e=0.3T_merge ≈ 32.0 Myr, M_chirp ≈ 778 M☉, f_peak ≈ 2.44 Hz中等质量黑洞双星
m₁=20 M☉, m₂=20 M☉, a=5,000,000 km, e=0T_merge ≈ 25.0 Myr, M_chirp ≈ 17.4 M☉, f_peak ≈ 110 Hz等质量恒星级双星

如何使用黑洞碰撞计算器

  1. 输入每个黑洞的质量,单位为太阳质量 (M☉)。恒星级黑洞通常约为 3 到 100 M☉;超大质量黑洞可超过 10⁹ M☉。
  2. 输入初始轨道分离,单位为公里。这就是双星轨道的初始半长轴。
  3. 将轨道偏心率设为 0(圆轨道)到 0.99(近径向轨道)之间。大多数被 LIGO 探测到的事件在进入探测频段时已经接近圆轨道。
  4. 输入倾角,单位为度(0° = 面朝上,90° = 侧视)。这会影响地球上观测到的引力波幅度,但不会改变并合时间。
  5. 点击计算,即可查看啁啾质量、估计并合时间、引力波能量、最终黑洞质量、史瓦西半径和峰值引力波频率。

常见问题

什么是啁啾质量,它为什么重要?
啁啾质量 M_chirp = (m₁m₂)^(3/5)/(m₁+m₂)^(1/5) 是引力波探测中最重要的单一参数。它决定引力波频率上升的速率(即啁啾率),因此天文学家即使在不知道单个质量之前,也能从波形中非常精确地测出啁啾质量。
并合时间估算有多准确?
这里使用的 Peters 公式适用于早期并合前螺旋阶段,此时分离距离远大于史瓦西半径。偏心率修正 (1−e²)^(7/2) 是一个近似,在 e ≲ 0.6 时效果很好。对于高偏心轨道或更紧密的分离距离,需要数值相对论才能给出精确估计。
为什么偏心率越高,并合越快?
在偏心轨道最接近时刻(近星点),天体移动最快且彼此最接近,此时辐射的引力波功率会急剧增加。每一圈辐射的能量更多,因此轨道收缩更快,与平均分离距离相同的圆轨道相比,并合时间会更短。
什么是 ISCO,为什么它定义峰值引力波频率?
最内稳态圆轨道(ISCO)标记了 Schwarzschild(不自旋)黑洞周围不再存在稳定圆轨道的边界。当并合前螺旋到达这一位置时,小天体会迅速向内坠落。ISCO 处的轨道频率乘以 2 后就是引力波频率,它代表并合前螺旋信号的最高频率,也是并合振铃开始的标志。
有多少能量会以引力波形式辐射出去?
对于质量相近的黑洞并合,数值相对论模拟表明大约有总质量的 4–8% 会以引力波形式辐射出去。此计算器用约化质量能量的约 5% 作为粗略估计。对于 GW150914,大约 3 个太阳质量(约占总质量的 5%)在不到一秒内转化为引力波能量。
这个计算器可以用于中子星并合吗?
该并合前螺旋公式同样适用于中子星–中子星(BNS)和中子星–黑洞(NSBH)双星。不过对于 BNS 事件,潮汐瓦解和中子星状态方程会引入这里未包含的修正。你可以把它用于数量级估算;若要精确的 BNS 结果,应使用后牛顿或数值相对论波形模型。