6.2 引力波天文学：聆听双黑洞并合

6.2 引力波天文学：聆听双黑洞并合

本节摘要：LIGO 用互成直角的四公里激光臂把 10⁻²¹ 的应变变成可读的干涉条纹——GW150914 的啁啾在十分之一秒内扫过 35 到 250 赫兹，反演出 36 加 29 太阳质量黑洞并合为 62 太阳质量、3 倍太阳质量化为涟漪。本节拆解干涉仪的测距工程与噪声预算，用代码合成一段啁啾并从频率演化反演啁啾质量，最后盘点多信使天文学的新账本。

2015 年 9 月 14 日的九毫秒

工程试运行还没正式开始，汉福德与利文斯顿两台相隔 3002 公里的探测器在 6.9 毫秒内先后记录下同一段波形：频率从 35 赫兹爬升到 250 赫兹，振幅先增后骤然衰减——教科书级的"啁啾"，但发生在引力波频段。过滤器扫过上千个理论模板（6.5 节数值相对论的产物），信号显著性超过 5 西格玛：这不是仪器哼哼，是十三亿年前两个黑洞的最后三圈舞。

参数反演靠的是波形各阶段携带的不同信息。旋近段的频率与频率变化率锁定啁啾质量（两质量的一个特定组合）；并合后的阻尼振荡（铃宕）频率与衰减率锁定最终黑洞的质量与自转；振幅结合距离给出光度距离。GW150914 的账本：36 加 29 太阳质量并合为 62 太阳质量，差额 3 倍太阳质量全部化为引力波，距离 410 兆秒差距——那次并合发生时，地球上才刚刚出现多细胞生命。

测 10 的负 21 次方：干涉仪工程

引力波不改变光的波长，它改变的是臂的长度——所以正确的仪器是测距机，不是测光机。迈克耳孙干涉仪的几何恰好合适：一束激光在分束器分成两半，跑完两条互成直角的臂再汇合。引力波到来时（图 6-1 的十字搓动），一条臂拉伸时另一臂压缩，两束光的相位差被来回放大。

图 6-2 L 型干涉仪与噪声预算

噪声预算的读法：不同频段有不同的"对手"——10 赫兹以下地震与悬挂热噪声统治，中段是镜面涂层的热涨落，高频是光子计数的散粒噪声。工程对策一整套：臂藏在高真空钢管里、40 公斤镜子挂在四级隔振塔上、激光在臂内往返 280 次叠加相位、再注入"压缩光"压低量子噪声。两台探测器相距三千公里，用 6.9 毫秒的时间差做方向三角化、也做真伪互证（一台的 glitches 不能伪造两台的同款波形）。

数值复演：合成一段啁啾并反演质量

# 啁啾合成与反演：从频率演化读出啁啾质量
import math
G = 6.674e-11; c = 2.99792458e8; Msun = 1.989e30
def chirp_mass_from_f_fdot(f, fdot):
    # 领先阶公式：啁啾质量由频率与频率导数唯一决定
    return (c**3/G) * ((5/96) * math.pi**(-8/3) * f**(-11/3) * fdot)**(3/5)
def fdot_from_chirp(Mc, f):
    Mc_kg = Mc*Msun
    return (96/5) * math.pi**(8/3) * (G*Mc_kg/c**3)**(5/3) * f**(11/3)

Mc = 28.0                       # GW150914 级别的啁啾质量
f = 35.0
print("合成：从 35 Hz 出发的频率爬升")
for step in range(5):
    fdot = fdot_from_chirp(Mc, f)
    dt = 0.02/f                 # 每 0.02 圈更新一次
    f2 = f + fdot*dt
    print(f"  f = {f:6.1f} Hz -> {f2:6.1f} Hz（fdot = {fdot:.1f} Hz/s）")
    f = f2
print(f"反演验证：由 f=35 Hz 与初段 fdot 反演啁啾质量")
Mc_back = chirp_mass_from_f_fdot(35.0, fdot_from_chirp(Mc, 35.0))
print(f"  输入 {Mc} 太阳质量，反演得 {Mc_back:.1f} 太阳质量")
print("两个观测量锁定一个质量组合——波形即天平")
​

日志要点：频率每上一阶，频率导数按 11/3 次方飙升——35 赫兹时每秒只爬几赫兹，200 赫兹时每秒数千赫兹，这就是啁啾"先懒后狂"的数学面目；反演闭合，说明啁啾质量确实被 (f, df/dt) 两个数锁死。旋近段读完质量，铃宕段（阻尼振荡）读最终黑洞的质量与自转（频率约 250 赫兹、衰减毫秒级，与数值相对论模板对表），振幅给距离——一个事件，四重账目。

# 双站定向：6.9 毫秒时间差能定多准？
c = 2.99792458e8
d = 3002e3                     # 两台探测器间距，米
dt = 6.9e-3
theta = math.degrees(math.acos(c*dt/d))
print(f"时间差 {dt*1e3:.1f} ms 对应来向与基线夹角约 {theta:.0f} 度")
print("单环带宽度有限——两台只能定出数百平方度的环带；")
print("三台成网才能收窄到几十平方度，给望远镜指路")
​

定向账本：两台探测器只能把源定在一个环带上（GW150914 的 600 平方度），2017 年起三台成网（加上处女座）收窄到几十平方度——这正是 GW170817 能有光学对应体的前提。

多信使天文学的新账本

2017 年 8 月 17 日，三台探测器同收一段 90 余秒的长啁啾（频率更低、时长更长——中子星双星的签名），两秒后伽马卫星暴闪，十一小时后光学望远镜在定向区内捉到千新星。GW170817 一举开出四门新账：状态方程（潮汐形变写入波形，中子星不能太硬）；短伽马暴起源（双中子星并合坐实）；金银等重元素产地（千新星的放射性发光功率对上快中子俘获核合成）；标准警笛测哈勃常数（引力波振幅自带光度距离，配红移即可量宇宙，不依赖距离阶梯——6.3 节的伏笔）。引力波天文学从"能探测"走到"能测量"，用了不到两年。

航行警告：引力波事件的"定位"与天文台的习惯不是一个精度量级——几十平方度对光学巡天够用，对需要盯住单个星系的观测仍然太大；下一代（宇宙望远镜 LISA、爱因斯坦望远镜）把角分辨率再压两个量级，是这一行公认的下一站。

关键直觉：电磁天文学看的是"物质在发光"，引力波听的是"质量在加速"。前者的信使是光子，后者传递信息的是几何本身——不遮挡、不散射、穿整个宇宙如入无人之境，这既是它难测的原因，也是它一旦测到就信息纯粹的原因。