黑洞碰撞計算器

計算雙黑洞系統的併合時間、引力波能量與最終黑洞性質。

輸入兩個黑洞的質量、初始軌道分離、偏心率與軌道傾角,即可使用 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 結果,應使用後牛頓或數值相對論波形模型。