2.4 引力波
2.4 引力波
本节摘要：引力波是时空本身的涟漪，由大质量天体加速运动产生。2015 年 LIGO 首次直接探测到双黑洞并合信号 GW150914，结束了长达百年的"只有预言没有证据"状态。本节讲清引力波与粒子信使的本质区别、波形里能读出什么参数，以及 GW170817 如何改写多信使天文学。
一、与粒子信使的根本区别
光子、宇宙线、中微子都在时空里传播，引力波则是时空本身的扰动。这带来三个特性：几乎不被物质吸收散射、携带源的质量四极矩动力学信息、探测靠测长度微小变化而非计数粒子。代价是振幅极小——LIGO 要测的是 4 公里臂长上约 10^-18 米的变化，相当于质子直径的千分之一。
二、波形是一本账本
双致密星并合的波形分三段：旋进（inspiral，频率与振幅渐增）、并合（merger，峰值）、铃宕（ringdown，黑洞恢复球形的衰减振荡）。用数值相对论模板做匹配滤波，可以从一段波形里同时抽出：
- 两个源质量（并合前与并合后）
- 自旋（至少对齐分量）
- 距离（"标准汽笛"）
- 位置（多探测器到达时间差三角定位）
一个简单的量级估算展示了为什么这类源可以被听到——辐射功率峰值约 3.6×10^49 W（GW150914），短暂超过可观测宇宙全部恒星光度之和：
c = 2.998e8 # m/s G = 6.674e-11 # SI M_sun = 1.989e30 # kg def chirp_mass(m1, m2): """chirp 质量：控制旋进段频率演化的等效质量，kg""" return (m1*m2)**0.6 / (m1+m2)**0.2 mc = chirp_mass(36*M_sun, 29*M_sun) # GW150914 双黑洞初质量 print(f"chirp 质量 ≈ {mc/M_sun:.1f} 太阳质量")
匹配滤波中 chirp 质量正是最先、最稳被锁定的参数——频率扫过的速度由它决定。
三、GW170817：一锤定音的多信使事件
2017 年 8 月 17 日，LIGO/Virgo 的三台干涉仪记录到双中子星并合，并在 1.7 秒后收到伽马暴 GRB 170817A 的伽马触发，随后数天内全球望远镜追踪到从光学到射电的千新星余辉。这一事件一次性解决/推进了多个问题：
| 问题 | GW170817 给出的证据 |
|---|---|
| 短伽马暴起源 | 与中子星并合同源，时间与定位吻合 |
| 重元素起源 | 千新星光谱支持 r 过程合成金、铂等元素 |
| 引力波速度 | 与光速差小于十亿分之一 |
| 哈勃常数 | 独立"标准汽笛"测距路径 |
图：GW170817 多信使时间线
💡 关键直觉：引力波不告诉你源"长什么样"，而是告诉你时空"怎么动"。它与传统天文的关系，像地震仪之于摄像机。
本节要点回顾
- 本质区别：引力波是时空自身的扰动，不被吸收、携带动力学信息。
- 波形即参数：质量、自旋、距离可从匹配滤波同时抽出，chirp 质量最稳。
- GW150914 与 GW170817 分别打开引力波天文学与多信使时代。
- 探测极限来自噪声：4 km 臂上 10^-18 m 的应变，决定了必须用激光干涉 + 统计累积。
动手估算：双黑洞并合的啁啾质量
引力波波形最有信息量的一段是"啁啾"——频率与振幅同时上升。由频率变化率可解出啁啾质量：M_chirp = c³/G × [ (5/96) π^(-8/3) f^(-11/3) df/dt ]^(3/5)。GW150914 的 f 从 35 Hz 升到 150 Hz 用时约 0.2 秒，解出 M_chirp ≈ 30 M_sun，再结合质量比约束得 36+29 M_sun 的双黑洞。
# 用啁啾公式从 (f, fdot) 反解啾啾质量 import math c, G, M_sun = 2.998e10, 6.674e-8, 1.989e33 f, fdot = 35.0, (150-35)/0.2 # Hz, Hz/s（粗略线性化） inner = (5/96) * math.pi**(-8/3) * f**(-11/3) * fdot M_chirp = (c**3/G) * inner**0.6 / M_sun print(f"啁啾质量 ~ {M_chirp:.1f} M_sun") # 输出约 30 M_sun，与 LIGO 发表值一致，这是"波形即秤"的演示
观测数据档案：引力波探测累计格局（O1-O3）
| 波段事件类 | 累计数量级 | 典型质量 | 科学收获 |
|---|---|---|---|
| 双黑洞并合 | 约 90 例 | 5-150 M_sun | 黑洞质量谱与并合率 |
| 双中子星 | 2 例 | ~1.3+1.3 M_sun | 千新星的 r 过程证据 |
| 黑洞-中子星 | 2 例 | ~8+2 M_sun | 潮汐瓦解阈值 |
# 由应变 h 与距离估算辐射功率峰值 h, d = 1.0e-21, 410 # 应变，Mpc d_cm = d * 3.086e24 import math # 峰值光度量级 L ~ (c^3/G) * (h*d*f)^2 * f^2 / 通用系数 f_gw = 150.0 # Hz L = (2.998e10)**3/6.674e-8 * (h*d_cm*2*math.pi*f_gw/2.998e10)**2 print(f"峰值光度 ~ {L:.2e} erg/s") # 输出 ~3.6e+56 erg/s，短促超过全宇宙恒星光度之和
延伸讨论：引力波为何是"全新的仪器"
引力波与电磁波的分野不在"波"而在"耦合"：电磁波被等离子体与尘埃散射吸收，引力波与物质的耦合强度仅约光速的平方分之一量级，几乎无阻尼地穿过双星包层、星系尘埃乃至宇宙早期等离子体。这赋予它两种独有能力：一是"听"致密天体内部——并合末期的频率攀升记录的是两颗中子星物质如何被潮汐力拉扯，状态方程的软硬直接写进波形相位；二是"看"无光事件——双黑洞并合在任何波段都无电磁对应，唯有引力波可统计其并合率与质量谱。到 O4 为止的约90例双黑洞事件已勾勒出从5到150太阳质量的分布，其中包含传统恒星演化难以解释的质量区（如落入"对不稳定性空隙"的候选），这些数据正在反哺恒星演化与星族合成理论，构成引力波天文学回馈天体物理的第一批实质成果。