跳轉到

2015 - LIGO 重力波偵測（GW150914）#

案件摘要#

2015 年 9 月 14 日 09:50:45 UTC，LIGO 的 Hanford（華盛頓）與 Livingston（路易斯安那）兩座干涉儀先後捕捉到同一個訊號：GW150914。這是兩個質量約 36 與 29 個太陽質量的黑洞合併成 62 個太陽質量黑洞的瞬間，约 3 個太陽質量的能量以重力波形式在 0.2 秒內輻射出去。人類第一次「聽見」時空的漣漪，開啟重力波天文學時代。Weiss、Barish、Thorne 獲 2017 年諾貝爾物理學獎。

前因 -- 為什麼會有這個案子#

  • 1916–1918：愛因斯坦從線性化場方程導出重力波，以光速 \(c\) 傳播、具兩個極化自由度。但他本人一度懷疑它只是座標變換的假象。
  • 1957 Chapel Hill 會議：Feynman 的「黏珠論證」等思想實驗確立重力波攜帶真實能量，偵測在原則上可行。
  • 1960s–70s：Weber 以共振棒（鋁圓柱）嘗試偵測，靈敏度不足且結果無法重現；同時 Hulse–Taylor 雙脈衝星提供了間接證據，激勵社群繼續尋找直接偵測之路。
  • 方案：Rainer Weiss（MIT）提出雷射干涉儀方案——讓重力波拉伸/壓縮空間，改變干涉儀兩臂的長度差，用雷射干涉讀出。Kip Thorne（Caltech）提供理論預測（波形、合併事件率），Ronald Drever 建構關鍵光學技術。
  • 工程難題：應變 \(h \sim 10^{-21}\) 意味著 4 km 臂長上的長度變化 \(\Delta L = hL \sim 4\times 10^{-18}\) m——比質子半徑（ \(10^{-15}\) m）還小三個數量級。地震、地震雜訊、熱雜訊、散粒雜訊全都足以掩蓋訊號，必須一一排除。

線索與推理 -- 數學式、程式、理論#

線索一：干涉儀原理#

Michelson 干涉儀中，雷射分為兩束分別往返兩臂。當重力波（ \(+\) 極化）沿垂直方向入射時，兩臂長度做差動振盪：

\[ L_x(t) = L\left(1 + \frac{h}{2}\right),\qquad L_y(t) = L\left(1 - \frac{h}{2}\right),\qquad \Delta L(t) = h\,L \]

往返光程差轉為干涉條紋的相位差 \(\Delta\phi = \dfrac{4\pi}{\lambda}\Delta L\) 。以 \(\lambda = 1064\) nm、 \(\Delta L \sim 4\times 10^{-18}\) m， \(\Delta\phi \sim 10^{-8}\) rad——只能靠法布里–珀羅共振腔（光在 4 km 臂內往返約 300 次，等效臂長放大至 1200 km）與功率回收、訊號回收技術讀出。

線索二：靈敏度與雜訊預算#

GW150914 的應變峰值 \(h \approx 1.0\times 10^{-21}\) ：

\[ \Delta L = h\,L = 10^{-21} \times 4\times 10^3\ \text{m} = 4\times 10^{-18}\ \text{m} \]

主要雜訊來源與頻段：低頻（ \(< 10\) Hz）為地震與牛頓引力梯度雜訊；中頻為懸掛鏡的熱雜訊；高頻（ \(> 200\) Hz）為光的散粒雜訊（ \(\propto 1/\sqrt{P_{\rm laser}}\) ，故需 100 W 量級雷射）。鏡子以四重擺懸掛隔振，測試質量 40 kg，真空管維持 \(10^{-9}\) Pa。

線索三：兩站同時觸發#

宇宙射線或本地擾動不可能同時在相距 3002 km 的兩站觸發。兩站訊號時間差 \(\Delta t = 6.9^{+0.5}_{-0.4}\) ms，恰好落在重力波以光速橫越兩站的最大時間差（ \(\le 10\) ms + 地球曲率修正）內——方向鎖定在南天球。這是「同案兩處現場」的鐵證。

理論：啁啾訊號與質量反演#

雙星系統的軌道頻率隨能量流失而升高，重力波頻率（兩倍軌道頻率）跟著升高——振幅越來越大、頻率越來越快的「啁啾」（chirp）。由廣義相對論後牛頓展開的領頭項：

\[ f_{\rm GW}(t) \simeq \frac{5}{256\pi}\left(\frac{G\mathcal{M}_c}{c^3}\right)^{-5/3}(t_{\rm merge}-t)^{-3/8},\qquad h(t) \propto \frac{(G\mathcal{M}_c)^{5/3}}{c^4 D_L} f_{\rm GW}^{2/3} \]

其中啁啾質量（chirp mass）：

\[ \mathcal{M}_c = \frac{(m_1 m_2)^{3/5}}{(m_1+m_2)^{1/5}} \]

從訊號可直接讀出 \(\mathcal{M}_c = 30^{+5}_{-4}\ M_\odot\) ；再由合併前的頻率上限（ \(f \approx 150\) Hz，更高頻時已進入最後不穩定軌道）推得 \(m_1 \approx 36\ M_\odot\) 、 \(m_2 \approx 29\ M_\odot\) 、光度距離 \(D_L \approx 410\) Mpc（紅移 \(z \approx 0.09\) ）。合併後由 ringdown（餘振）頻率與阻尼時間驗證：這是一個質量 62 \(M_\odot\) 、自旋 0.67 的克爾黑洞——與數值相對論模擬的波形完全吻合。輻射能量：

\[ E_{\rm GW} = (36+29-62)\,M_\odot c^2 \approx 3\,M_\odot c^2 \approx 5.4\times 10^{47}\ \text{J} \]

瞬時峰值光度約 \(3.6\times 10^{49}\) W——超過可觀測宇宙全部恆星光度之和。

Python 模擬：啁啾波形示意#

import numpy as np
import matplotlib.pyplot as plt

# 啁啾波形（領頭階後牛頓近似）
G, c = 6.674e-11, 2.998e8
Msun = 1.989e30
m1, m2 = 36*Msun, 29*Msun
Mc = (m1*m2)**0.6 / (m1+m2)**0.2        # chirp mass
DL = 410e6 * 3.086e16                    # 光度距離 (m)

# 頻率-時間關係： df/dt ∝ f^(11/3)，反解 f(t)
tc = 0.0
t = np.linspace(-0.45, -0.005, 4000)     # 合併前 0.45 秒
tau = tc - t
f = (5/256) * (G*Mc/c**3)**(-5/3) * tau**(-3/8) / np.pi

# 截斷在最後穩定軌道頻率
f_isco = c**3 / (6**1.5 * np.pi * G * (m1+m2))
f = f[f < f_isco]; tau = tau[:len(f)]

# 應變振幅 ∝ f^(2/3) / DL
h_amp = 4 * (G*Mc)**(5/3) * (np.pi*f)**(2/3) / (c**4 * DL)
# 高斯衰減包絡模擬 ringdown 結尾
h = h_amp * np.cos(2*np.pi*f*tau) * np.exp(-((t[:len(f)]+0.02)/0.03)**2)

fig, ax = plt.subplots(2, 1, figsize=(9, 6), sharex=True)
ax[0].plot(t*1e3, h*1e21, 'b-', lw=0.8)
ax[0].set_ylabel('strain $h$ [$\\times 10^{-21}$]')
ax[0].set_title('GW150914 chirp (schematic): 36+29 $M_\\odot$ -> 62 $M_\\odot$')

# 瞬時頻率（啁啾）
ax[1].plot(t[:len(f)]*1e3, f, 'r-', lw=1.5)
ax[1].set_xlabel('time to merger [ms]'); ax[1].set_ylabel('$f_{GW}$ [Hz]')
ax[1].set_yscale('log')
plt.tight_layout(); plt.show()

輸出顯示經典的啁啾形貌：振幅隨頻率攀升而增大，最後在約 150 Hz 處到達最後穩定軌道並合併——這正是 Hanford/Livingston 記錄到的圖形（經匹配濾波後的「漂亮的波形」）。

結案 -- 後果與影響#

  • 2016 年 2 月 11 日：LIGO 公布 GW150914，論文發表於 PRL，作者超過 1000 人（LIGO 科學合作組織 + Virgo）。
  • 2017 諾貝爾物理學獎：Rainer Weiss（干涉儀構想與領導）、Barry Barish（把 LIGO 從原型變成大科學工程）、Kip Thorne（理論與波形預測）。
  • 重力波天文學時代：
  • GW170817（2017）：雙中子星合併，同時被 Fermi 衛星觀測到伽瑪射線暴 GRB 170817A——多信使天文學誕生；重力波與電磁波聯合定位，確認千新星（kilonova）重元素來源，並以 \(f_{\rm GW}\) 與電磁紅移直接檢驗宇宙膨脹（ \(H_0\) 獨立測量）。
  • O1/O2/O3 觀測季累積偵測上百次黑洞合併；2023 年起進入 O4，事件率達每週數次。
  • 未來：太空干涉儀 LISA（2030s）瞄準毫赫茲頻段，將偵測超大型黑洞雙星。

關鍵人物與文獻#

  • Rainer Weiss（1932–2025）：干涉儀設計者，2017 諾獎。
  • Barry Barish（1935–2025）：LIGO 前主任，2017 諾獎。
  • Kip Thorne（1940–）：理論家，2017 諾獎。
  • Ronald Drever（1931–2017）：關鍵光學技術（錯失諾獎，因其辭世早於頒獎）。
  • 文獻：
  • Abbott et al. (LIGO/Virgo), PRL 116, 061102 (2016) —— GW150914 偵測論文。
  • Abbott et al., ApJL 848, L12 (2017) —— GW170817 多信使觀測。
  • Abbott et al., Nature 551, 85 (2017) —— 重力波標準汽笛測 \(H_0\) 。