6.1 引力波：时空的涟漪航讯

6.1 引力波：时空的涟漪航讯

本节摘要：把弱场近似的"静态"限制放开，线性化场方程成为波动方程——度规的微小起伏 h 以光速传播，这就是引力波。它有两个十字形偏振（加号与叉号），只会由不对称的加速质量（四极矩变化）产生；经过之处，一米长的尺子真的伸缩约 10⁻²¹ 米。本站推导波动解的骨架、算清赫尔斯-泰勒脉冲双星与 GW150914 的应变与能量账，说明引力波探测为何是全新的感官。

放开静态限制：度规的波

第 3.3 节做弱场近似时加了三条限制，其中一条是静态——场不随时间变。把这条剪掉，度规写成 g_μν = η_μν + h_μν 且 |h| 远小于 1，线性化后的场方程立刻显出本来面目：一个波动方程

□h̄_μν = −(16πG/c⁴) T_μν

它的结构与电磁理论惊人平行：电磁波是麦克斯韦方程的波解，引力波是场方程线性化的波解；波速同样是 c（度规起伏踩着光锥奔跑）。真空处 T 为零，方程说起伏自由传播；有物质处，物质的四极矩变化是波源。在场方程发表的同年（1916 年），爱因斯坦就注意到了这个解，还一度怀疑它是否携带能量（"坐标伪波"之争）——直到同行用理想探测器账本证明波的潮汐效应是坐标无关的真物理。

引力波经过时发生什么？想象一队测试粒子排成圆圈，波沿垂直方向穿过：圆圈交替地被拉长压扁——先横向拉伸纵向压缩，再反过来，循环往复。这不是尺子"看起来"变了：激光在臂里跑的距离真实地变化，干涉条纹真实地移动。应变的定义 h = ΔL/L，即"长度的相对变化"。

十字形偏振：时空的两种搓法

广义相对论只允许两种偏振模式（广义的修改引力理论可以更多，这是 6.5 节的检验点）：加号偏振沿参考轴拉伸压缩，叉号偏振沿着旋转 45 度的轴拉伸压缩。一列波通常是两者的叠加。

图 6-1 引力波的两种偏振：测试粒子圆圈的呼吸

为什么必须是四极：对称不发声

波源的条件比电磁波苛刻得多。电磁辐射的最低阶是偶极（电荷来回振荡即可）；引力没有负质量，单极辐射被质量守恒禁止，偶极辐射被动量守恒禁止——第一阶可用的是四极。翻译成工程语言：转得再快的完美球也不发声，必须"不对称地加速"。圆轨道双星恰好是天然的四极机器：两团质量互相绕转，四极矩每转一圈变化两次，连续发射引力波。

引力波带走能量，波源必然损失——这个"副作用"成了引力波存在的第一个证据。赫尔斯与泰勒 1974 年发现的脉冲双星 PSR B1913+16：两颗中子星每 7.75 小时互绕一周，轨道周期每年缩短 76.5 微秒。三十年累计的轨道收缩曲线与广义相对论的引力波辐射预言吻合到千分之二——还没等到直接探测，引力波已经靠这笔间接账拿了 1993 年诺贝尔奖。

能量本身的记账方式要单独记一笔（2.5 节埋过伏笔）：引力场能量没有局部的张量密度，只有整体的能量流动账（通过无穷远球面的能流），这与场方程的非线性一脉相承——几何自己参与携带能量，就无法在每一点上把"物质能"与"场能"切开。好在"远处收到多少能量"是定义良好的：它正比于时间导数平方的约定积分，对双星系统积分出来，正是脉冲双星轨道衰减的那条曲线。理论账本与天文观测在 2.4 乘 10 的负 12 次方每秒的轨道变化率上对齐——这份耐心，是引力波时代之前最精密的广义相对论检验。

数值对账：潮汐量级与能量洪流

# 应变速算：两个标志性波源
import math
G = 6.674e-11; c = 2.99792458e8; Msun = 1.989e30
Mpc = 3.086e22
def strain(Mchirp_solar, f_hz, dist_m):
    Mc = Mchirp_solar*Msun
    return 4*(G*Mc)**(5/3)*(math.pi*f_hz)**(2/3)/(c**4*dist_m)
# 赫尔斯-泰勒：啁啾质量约 1.22 太阳质量，轨道频率对应引力波频率约 7.3e-5 Hz
h_ht = strain(1.22, 7.3e-5, 6.4e3*3.086e16)
# GW150914：啁啾质量约 28 太阳质量，150 Hz，410 Mpc
h_gw = strain(28.0, 150.0, 410*Mpc)
print(f"赫尔斯-泰勒双星   h ≈ {h_ht:.1e}（频率过低，LIGO 频段外）")
print(f"GW150914 并合瞬间 h ≈ {h_gw:.1e}（恰好压在探测线之上）")
print("一米长的臂在并合瞬间的变化量：", f"{h_gw*1e18:.4f} 乘 10 的负 18 次方米")
​

账本读法：GW150914 在 410 兆秒差距外掠过地球时，应变量级 10⁻²¹——四公里的臂长变化约 10⁻¹⁸ 米，千分之一颗质子的直径。量到了"质子直径的千分之一"却确凿无疑，因为两台相隔三千公里的探测器听到了同一个波形、相差 6.9 毫秒（正好是光速穿越的时间）。

# 能量洪流：并合瞬间的引力波光度
G = 6.674e-11; c = 2.99792458e8; Msun = 1.989e30
L_peak = 3.6e49            # GW150914 峰值光度，瓦（文献值）
L_sun = 3.828e26
E_rad = 3.0*Msun*c**2      # 辐射掉约 3 太阳质量的能量
print(f"峰值光度：{L_peak:.1e} W = {L_peak/L_sun:.1e} 倍太阳光度")
print(f"可观测宇宙全部恒星光度约 2e49 W——并合瞬间超过它")
print(f"总辐射能量：{E_rad:.2e} J，相当于 3 个太阳质量彻底化为涟漪")
​

第二笔账更夸张：并合峰值光度 3.6×10⁴⁹ 瓦，瞬时超过可观测宇宙全部恒星的光度总和；三倍太阳质量的静能在一百多毫秒内化 pure粹为时空涟漪。宇宙里最猛烈的能量事件，光学望远镜几乎看不见——全部账目记在引力波波段。这就是 6.2 节要开的"新感官"的意义。

⚠️ 航行警告："引力波是时空的振动"这句话容易被听成修辞。准确版本：引力波经过时，两点间距离与波长同量级的相对变化是真实的、可由干涉仪测量的物理事件；不接受"只是坐标在晃"的解读。

本站要点

• 波动解：放开静态限制后线性化场方程成为波动方程，度规起伏以光速传播，与电磁波结构平行；
• 两种偏振：加号与叉号十字形搓动，测试粒子圈交替拉压；偏振数是检验修改引力的旋钮；
• 四极门槛：质量守恒禁单极、动量守恒禁偶极，第一阶可用波源是四极——对称不发声；
• 间接证据：赫尔斯-泰勒脉冲双星轨道衰减三十年吻合到千分之二，1993 年诺奖入账；
• 量级账：GW150914 应变 10⁻²¹、四公里臂变化千分之一质子直径；峰值光度超全宇宙恒星总和，3 倍太阳质量化为涟漪；
• 下一站：量这么小的东西如何测量——LIGO 的四公里激光干涉工程与并合事件的反演算法。