6.3 新式证词：引力波警报器与21厘米信号

6.3 新式证词：引力波警报器与21厘米信号

本节摘要：前两节的装备仍在电磁波的账本里记账，本节介绍两路完全绕开电磁波的证词：双致密星并合的引力波自带距离（标准警报器），可以独立裁决哈勃常数与距离阶梯；氢原子的 21 厘米谱线把可测的膨胀史推到恒星点亮之前。本节演算啁啾质量反演与频率红移换算，并用引力波测哈勃常数的完整案例说明新证词的成色与边界。

第二章的五路证词有一个共同底座：全都经由光。光传播会被途中的物质弯折、消光、散射，证词因此总有"中途被动过手脚"的质疑。本节的两路新证词直接换底座——引力波几乎不与物质作用，21 厘米信号则来自宇宙最黑暗年代的氢原子云。对暗能量案卷而言，它们不是第五路、第六路的简单加号，而是能把 1.2 节哈勃常数读数独立复核的"编制外证人"。

引力波标准警报器：波形自带距离

双中子星或黑洞并合的引力波频率上升过程（啁啾信号）由一个组合量决定——啁啾质量。频率变化率把啁啾质量反演出来，波形的振幅又给出光度距离：距离不经过任何亮度校准，这就是"标准警报器"名字的由来。先练反演。

# 演算 1：啁啾质量反演与并合倒计时
import math
def mchirp(m1, m2):
    # 啁啾质量 = (m1*m2)^(3/5) / (m1+m2)^(1/5)
    return (m1*m2)**0.6/(m1+m2)**0.2
for pair in ((1.4, 1.4), (1.4, 10.0), (30.0, 30.0)):
    print("m1=%.1f m2=%.1f -> 啁啾质量 = %.3f 太阳质量" % (pair[0], pair[1], mchirp(*pair)))
# 双中子星（1.4+1.4）在 10 与 100 赫兹处的并合倒计时
# t(f) = 5/256 × (G*Mc/c^3)^(-5/3) × (pi*f)^(-8/3)
G = 6.674e-11; c = 2.998e8; msun = 1.989e30
Mc = mchirp(1.4, 1.4)*msun
for f in (10.0, 100.0):
    t = 5.0/256*(G*Mc/c**3)**(-5/3.0)/(math.pi*f)**(8/3.0)
    print("f = %5.0f Hz -> 距并合约 %.0f 秒" % (f, t))
# 输出:
# m1=1.4 m2=1.4 -> 啁啾质量 = 1.219 太阳质量
# m1=1.4 m2=10.0 -> 啁啾质量 = 2.994 太阳质量
# m1=30.0 m2=30.0 -> 啁啾质量 = 26.117 太阳质量
# f =    10 Hz -> 距并合约 1002 秒
# f =   100 Hz -> 距并合约 2 秒
​

信号读法有两层。啁啾质量只依赖波形形状，与宇宙学参数无关——它是一把免检的内在标尺；倒计时则说明警报器的工作窗口：10 赫兹处还剩一刻钟，100 赫兹处只剩两秒，探测网络要提前预警就需要低频端足够灵敏。距离反演出来之后，剩下的步骤与 2.1 的超新星读法镜像对称：距离有了，红移从并合宿主星系的光谱里取，两者相除就是哈勃常数。

案例：一记警报器能校出什么

背景：某次双中子星并合被三台探测器同时捕获，宿主星系的光学对应体在几小时内找到，红移约零点零一。操作：由波形振幅反演光度距离（数十兆秒差距量级），由宿主光谱取红移，两值相除得到哈勃常数；同时检查波形里是否有偏白的谱偏差——那可能是宇宙加速度在波形频移上的指纹（红移依赖的距离-频率关系）。结果：单次事件的哈勃常数初值落在七十上下，误差带宽达两位数——不是精度不足，而是单次事件的本动速度干扰（几十兆秒差距上 peculiar 速度占红移的比重可观）与振幅反演的取向简并。解读：警报器的价值在"编制外"——它不依赖造父变星与超新星链条（1.2 节的梯子），因而能仲裁哈勃张力（7.3 节）哪一端出了问题。变式：事件数累积后误差按样本数反比收缩（距离主导型误差），约百次独立并合即可把哈勃常数压进百分之一量级——这是把 1.2 的读数从"案发记录"升级为"复核基准"的路线图。

21 厘米：把测距史推回黑暗年代

第二路新证词来自最普通的氢原子：基态超精细结构跃迁发出 1420.4 兆赫的谱线，对应 21 厘米波长。第一代恒星光把它周围的氢电离，剩下的中性氢像一片片暗影刻在背景辐射上——每个红移处的 21 厘米信号都是一张"当时的中性氢分布图"。它的宇宙学用法与 BAO 同源：氢云里同样冻结着 147 兆秒差距的声学尺，且红移可以推到远超星系巡天极限的黑暗年代。

# 演算 2：21 厘米信号的观测频率（ν = 1420.40575 MHz / (1+z)）
for z in (10.0, 50.0, 100.0, 200.0):
    print("z = %5.0f -> 21cm 频率 = %7.2f MHz" % (z, 1420.40575/(1+z)))
# 输出:
# z =    10 -> 21cm 频率 =  129.13 MHz
# z =    50 -> 21cm 频率 =   27.85 MHz
# z =   100 -> 21cm 频率 =   14.06 MHz
# z =   200 -> 21cm 频率 =    7.07 MHz
​

频率表划出了战场：红移 10 到 200 的信号落在 7 到 129 兆赫之间，全部低于调频广播的频率下限附近或与其重叠。这正是 21 厘米证词的天价成本——银河系与地球电波的前向污染比信号强一万到十万倍，百兆赫以下还有电离层反射挡路，严肃的探测方案已经指向月球背面。但回报同样明确：红移 50 以上的膨胀史是全案唯一没有电磁证词的空白页，任何能读出那页声学尺的设备，都直接给 4.2 节的高红移测谎曲线补上最深的一段。

易错点：警报器不是蜡烛，暗警报也另有用法

标准警报器与标准烛光的区别值得单独立档。烛光的路数是"先标定绝对亮度、再量流量"——校准链就是它的软肋（1.2 节整段在讲校准）；警报器是"波形自持"——距离从广义相对论的波形方程里解出来，没有可漂移的标定环节。它的软肋换了一处：红移必须从宿主星系取，若电磁对应体找不到（黑洞并合多数如此），距离有了、红移没有，哈勃常数照样开不了口。这类"暗警报"并非废票：把大量无对应体事件与星系巡天目录做统计配对，红移信息以概率形式补齐，单次精度差些、样本却能翻倍。另一处常见误读是把警报器当"更准的蜡烛"——它的单次精度其实不如 Ia 型超新星，价值在误差来源正交：亮度校准、尘埃、消光这些蜡烛案的陈年干扰项，在波形里根本不存在。顺着这条记档标准回看本章四台设备，会发现同一逻辑贯穿全章——每台装备的价值排序，独立性与校准纪律永远排在灵敏度前面。

新证词的边界与全章收线

两路新证词的短板要如实入档。引力波警报器的事件率与定位精度仍受探测器灵敏度限制，一次并合给不出一个哈勃常数；21 厘米在前景污染面前尚未完成第一次干净的探测，目前只有上限级别的结论。它们的地位因此是"预备证人"而非"在案证人"——可一旦作证，独立性远超现有五路。装备升级一章至此收线：6.1 的光谱阵列、6.2 的空间哨所、本节的两路新底座，四件装备对应第四章判据的四个收紧方向（w 读数、w 演化、增长指数、哈勃常数）。装备只负责把读数测准，读数背后的动机问题——真空能为何这么小、时机为何这么巧——移交第七章的结案陈词。