重力波天文學進階：應變、啁啾與配對濾波的物理

為什麼一個比質子還小的位移，能告訴我們十三億光年外發生了什麼？

你已經知道重力波是時空的漣漪，由加速運動的質量產生，並由 LIGO 這類干涉儀偵測。但這裡有個讓人不安的數字：2015 年那場著名的 GW150914 事件，在 LIGO 的 4 公里長手臂上造成的長度變化，只有約 $10^{-18}$ 公尺——大約是質子直徑的千分之一。

這引出一連串入門篇通常略過、卻是重力波天文學真正核心的問題：人類究竟如何「聽見」這麼微弱的訊號？為什麼我們能從一段只有零點幾秒的波形，反推出兩個黑洞的質量、自旋、距離與所在方向？又為什麼有些重力波我們現在的儀器永遠聽不到，得靠「把整個銀河系變成一台偵測器」才行？

這篇進階篇要做的，是把入門篇留下的「黑盒子」一個個打開：應變的精確數學、波形如何編碼物理、配對濾波（matched filtering）如何從雜訊中撈出訊號，以及為什麼重力波頻譜橫跨二十多個數量級、需要完全不同的偵測策略。

應變 h：把時空的拉伸寫成一個數字

入門篇說重力波會「拉伸與擠壓」空間。進階一點，我們需要一個量來精確描述這件事——這就是應變（strain）$h$。它的定義非常直觀：

$$ h = \frac{\Delta L}{L} $$

其中 $L$ 是兩個自由漂浮的測試質量之間的距離，$\Delta L$ 是重力波經過時造成的距離變化。應變是無量綱的，因為它是長度除以長度。

關鍵在於：重力波是個潮汐效應。它不是把所有東西往同一個方向推，而是在垂直於傳播方向的平面上，一個軸拉長、另一個軸縮短，半個週期後互換。這就是著名的「加號偏振」（plus polarization, $h_+$）與「叉號偏振」（cross polarization, $h_\times$）。一圈測試質點會在波經過時交替變成橫橢圓與直橢圓。

對 GW150914，峰值應變約 $h \approx 10^{-21}$。把這個數字放進 LIGO 的 4 公里手臂：

$$ \Delta L = h \cdot L = 10^{-21} \times 4 \times 10^{3}\ \text{m} = 4 \times 10^{-18}\ \text{m} $$

這就是那個小於質子直徑的位移。能測到它，靠的是干涉測量把光程差轉成光的相位差，再加上 Fabry-Pérot 共振腔讓光在手臂裡來回數百次、等效手臂長度放大到上千公里，以及功率回收（power recycling）把雷射有效功率推到數百千瓦。雜訊壓制（地震隔離、量子壓縮光 squeezed light）才是真正的工程奇蹟，但物理核心就是這個無量綱的 $h$。

為什麼振幅隨距離只衰減 1/r，而不是 1/r²？

這是重力波天文學一個極其重要、卻常被誤解的點。你熟悉的光、聲音，其能量通量（intensity）隨距離以 $1/r^2$ 衰減——這是球面波能量守恆的必然結果。重力波的能量通量也是 $1/r^2$。

但重力波偵測器測的不是能量通量，而是應變振幅 $h$。而振幅與通量的關係是通量正比於振幅平方（$F \propto h^2$）。因此：

$$ F \propto \frac{1}{r^2} \quad\Longrightarrow\quad h \propto \frac{1}{r} $$

這個 $h \propto 1/r$ 的線性衰減，是重力波天文學能「看得很遠」的根本原因。對比一下：如果你把偵測器的靈敏度（最小可測 $h$）提升兩倍，你能探測的距離也增加兩倍，而能探測的體積（正比於 $r^3$）增加八倍。換句話說，靈敏度每改善一個因子，可觀測的事件率就以三次方放大。這就是為什麼每一代干涉儀升級都帶來偵測數量的爆炸性成長：從 O1 的零星幾起，到 O3 之後累積上百起併合事件。

啁啾：一段波形如何把整個系統的物理寫進去

雙緻密星（黑洞或中子星）併合的波形有個招牌特徵：頻率與振幅都隨時間快速上升，直到併合瞬間達到峰值——這在頻譜圖上是一條向上掃的曲線，聽起來像鳥叫，所以叫啁啾（chirp）。

啁啾不是裝飾，它編碼了系統最重要的物理量。在系統還相隔很遠、可用牛頓加廣義相對論微擾（post-Newtonian）描述的「旋近（inspiral）」階段，重力波頻率 $f$ 的演化由一個組合質量主導，稱為啁啾質量（chirp mass）$\mathcal{M}$：

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

頻率隨時間的變化率直接給出啁啾質量：

$$ \dot{f} = \frac{96}{5}\,\pi^{8/3}\left(\frac{G\mathcal{M}}{c^3}\right)^{5/3} f^{11/3} $$

這個式子很美：只要量測波形頻率上升的速度，就能解出 $\mathcal{M}$，完全不需要事先知道距離或方向。啁啾質量是雙星質量資訊裡最容易、最精準被測到的量；個別質量 $m_1$、$m_2$ 與質量比要靠更高階的相位修正才能拆開，誤差通常較大。

波形可分成三段，各自由不同的物理主宰：

• 旋近（inspiral）：兩星互繞、軌道收縮，可用解析的 post-Newtonian 展開描述。
• 併合（merger）：兩個視界融合，時空強烈非線性，只能靠數值相對論（numerical relativity）用超級電腦解愛因斯坦方程式得到。
• 鈴宕（ringdown）：併合後形成的單一黑洞像被敲響的鐘一樣振盪、輻射掉殘餘形變，回到平滑的克爾（Kerr）黑洞。這段是衰減的正弦波，其頻率與衰減時間由最終黑洞的質量與自旋唯一決定——這是檢驗「無毛定理（no-hair theorem）」的黃金測試場。

動手算一下：兩個黑洞併合前最後一刻在做什麼？

讓我們估算 GW150914 那兩個約 $30\,M_\odot$ 黑洞，在訊號剛進入 LIGO 頻段（重力波頻率約 $f_{\rm GW} = 35$ Hz）時的軌道狀態。

重力波頻率是軌道頻率的兩倍（因為質量分布的「四極矩」每半圈重複一次），所以軌道頻率 $f_{\rm orb} = f_{\rm GW}/2 = 17.5$ Hz——它們每秒互繞約 17.5 圈。

用克卜勒第三定律估軌道半長軸（取總質量 $M = 60\,M_\odot \approx 1.2 \times 10^{32}$ kg）：

$$ a = \left(\frac{G M}{4\pi^2 f_{\rm orb}^2}\right)^{1/3} $$

代入 $G = 6.67\times10^{-11}$、$f_{\rm orb}=17.5$ Hz：

$$ a \approx \left(\frac{6.67\times10^{-11} \times 1.2\times10^{32}}{4\pi^2 \times 17.5^2}\right)^{1/3} \approx 6.6 \times 10^{5}\ \text{m} \approx 660\ \text{km} $$

兩個各重 30 倍太陽的黑洞，此刻相距只有約 660 公里——比台灣南北還短——卻每秒互繞 17.5 圈，繞行線速度已達光速的相當比例。再過約 0.2 秒，它們會以接近光速的一半相撞併合，在不到一秒內輻射出約 $3\,M_\odot$ 的質量等效能量，那一刻的瞬時功率超過可觀測宇宙中所有恆星光度的總和。這就是為什麼如此遙遠的事件，仍能在地球上留下可測的痕跡。

配對濾波：在雜訊汪洋裡撈出已知形狀的針

GW150914 的應變只有 $10^{-21}$，而 LIGO 的雜訊在很多頻段同樣甚至更大。那訊號怎麼會「看得見」？答案是配對濾波（matched filtering）——重力波資料分析的核心技術。

直覺是這樣：雜訊是隨機的，沒有固定形狀；但重力波訊號的形狀我們事先就算得出來（透過 post-Newtonian 與數值相對論）。如果你拿一個已知的理論波形 $template$，沿著資料時間軸滑動、做相關運算，當 template 與藏在雜訊裡的真實訊號對齊時，相關值會出現一個尖峰；而雜訊因為隨機，平均下來會互相抵消。

更精確地說，配對濾波計算資料 $d(t)$ 與模板 $h(t)$ 的雜訊加權內積：

$$ \langle d \mid h \rangle = 4\,\mathrm{Re}\int_0^\infty \frac{\tilde{d}(f)\,\tilde{h}^*(f)}{S_n(f)}\,df $$

其中 $\tilde{d}$、$\tilde{h}$ 是傅立葉變換，$S_n(f)$ 是偵測器的功率譜密度（power spectral density）——它描述每個頻率的雜訊有多強。注意分母 $S_n(f)$：它自動讓演算法在雜訊小的頻段「相信」資料、在雜訊大的頻段「打折」資料。這正是配對濾波在高斯雜訊下被證明為最佳線性偵測器的原因。

實務上，分析團隊得準備一個模板庫（template bank），涵蓋數十萬種不同質量與自旋組合的波形，因為你不知道下一個訊號是哪種系統。每個模板都跟資料做相關，找出超過門檻的訊噪比（SNR）尖峰，再用多座偵測器的時間一致性與假警報率（false alarm rate）統計確認它不是雜訊巧合。GW150914 的偵測顯著性超過 $5.1\sigma$，相當於每二十萬年才會出現一次的雜訊巧合——這就是為什麼我們能如此有信心地說「這是真的」。

重力波頻譜：為什麼一台偵測器永遠不夠

電磁波天文學早就懂這件事：你不能用同一台望遠鏡看無線電波和伽馬射線。重力波也一樣，而且頻段跨度更誇張——從 LIGO 的數百赫茲，到脈衝星計時陣列追的奈赫茲（nHz），橫跨超過十個數量級的頻率。

不同頻段對應完全不同的天體源，需要完全不同的偵測技術：

• 地面干涉儀（10 Hz – 數 kHz）：LIGO、Virgo、KAGRA。看恆星級黑洞與中子星併合。受限於地表震動與儀器尺寸，低頻測不到。
• 太空干涉儀（mHz）：規劃中的 LISA，三顆衛星組成邊長 250 萬公里的三角形干涉儀，在太空中不受地震干擾。目標是超大質量黑洞（百萬到億倍太陽質量）併合、極端質量比旋近（EMRI）。
• 脈衝星計時陣列（nHz）：NANOGrav、EPTA 等。把銀河系裡數十顆毫秒脈衝星當成宇宙級的時鐘網路，監測它們脈衝到達時間的微小、相關的偏移。2023 年這些團隊宣布了與重力波背景一致的證據，很可能來自宇宙中無數對超大質量黑洞的疊加哼鳴。

為什麼頻段決定源？因為雙星併合的特徵頻率反比於系統大小，而系統大小正比於質量。質量越大，最內穩定圓軌道（ISCO）越大、軌道頻率越低，重力波頻率也越低。恆星級黑洞在數百赫茲尖叫，億倍太陽質量的黑洞則以奈赫茲低吟——慢到一個週期要好幾年。

重點回顧

• 應變 $h = \Delta L / L$ 是描述重力波強度的無量綱量；它是潮汐性的，有 $h_+$ 與 $h_\times$ 兩種偏振，對 LIGO 而言峰值約 $10^{-21}$。
• 因為偵測器測的是振幅而非能量通量，重力波振幅隨距離以 $h \propto 1/r$ 線性衰減，使靈敏度每改善一倍、可觀測體積增加八倍。
• 啁啾質量 $\mathcal{M}$ 由頻率上升速率直接決定，是雙星系統裡最容易精確測得的量；波形分旋近、併合、鈴宕三段，分別由微擾理論、數值相對論、克爾黑洞準正規模主導。
• 配對濾波用已知理論模板與資料做雜訊加權相關，是高斯雜訊下的最佳偵測器，靠模板庫掃描數十萬種波形組合。
• 重力波頻譜橫跨十多個數量級，地面干涉儀、太空 LISA、脈衝星計時陣列各看不同質量尺度的天體，一台偵測器看不到整個宇宙。

深入探討（研究所視角）

若你想再往下走，以下是幾個把上述圖像建立在嚴謹理論之上的方向。

橫向無跡規範與四極公式。 重力波在弱場下是度規對平直時空的微擾 $g_{\mu\nu} = \eta_{\mu\nu} + h_{\mu\nu}$，$|h_{\mu\nu}| \ll 1$。選用橫向無跡規範（transverse-traceless gauge）後，真空中的擾動滿足波動方程 $\Box \bar{h}_{\mu\nu} = 0$，只剩兩個物理自由度（即 $h_+$、$h_\times$）。輻射的領導項來自源的質量四極矩二階時間導數：

$$ h_{ij}^{TT} = \frac{2G}{c^4 r}\,\ddot{I}_{ij}^{TT} $$

注意這裡沒有單極（質量守恆）也沒有偶極（動量守恆）輻射——這是重力波與電磁波的根本差異，也解釋了為何球對稱的塌縮（如理想化的超新星）不輻射重力波。係數裡的 $G/c^4 \approx 10^{-44}$ 是個極小的數，這正是重力波如此微弱、需要天體級質量加速到相對論性速度才能產生可測訊號的原因。

ISCO 與強場印記。 post-Newtonian 展開在旋近末期失效，因為速度趨近光速、場強趨近強場。對非自旋黑洞，史瓦西時空的最內穩定圓軌道在 $r_{\rm ISCO} = 6GM/c^2$，標誌旋近結束、動態併合開始。黑洞自旋會移動 ISCO 位置（順行自旋讓它更靠近視界），因此自旋會在波形相位裡留下印記——這是我們能從波形反推自旋的物理基礎，也是測量有效自旋參數 $\chi_{\rm eff}$ 的依據。

鈴宕與黑洞光譜學。 鈴宕階段的擾動可展開成克爾黑洞的準正規模（quasi-normal modes, QNM）疊加，其複數頻率（實部為振盪頻率、虛部為衰減率）僅由最終黑洞質量與自旋決定。若能在單一事件中量到兩個以上獨立模態並驗證它們互相一致，就能檢驗無毛定理、甚至探測對廣義相對論的偏離——這個新興領域稱為黑洞光譜學（black hole spectroscopy），是次世代偵測器（如愛因斯坦望遠鏡、宇宙探索者）的核心科學目標之一。

參數估計的貝氏框架。 真實分析不是只報「測到一個訊號」，而是對整個參數空間（質量、自旋、距離、傾角、天空位置、合併相位等十多個維度）做貝氏後驗推斷，用隨機抽樣方法（如 MCMC 或 nested sampling）建構後驗分布。距離與傾角高度簡併（兩者都影響振幅），天空定位則靠多座偵測器間的訊號到達時間差三角定位，這也是為什麼偵測器網路越大、定位越準——對中子星併合的電磁對應體後續追蹤（如 GW170817 的千新星）至關重要。

如果你對其中任一條線索著迷，建議從 Maggiore 的《Gravitational Waves》或 Creighton & Anderson 的資料分析專書入手，並試著自己用公開的 GWOSC（Gravitational Wave Open Science Center）資料，跑一次配對濾波，親手把 GW150914 從雜訊裡撈出來。那一刻，你就真正聽見了時空的聲音。