1.5 热宇宙的统计物理：温度就是时间码
1.5 热宇宙的统计物理：温度就是时间码
本节摘要：早期宇宙处于热平衡，温度与时间一一对应，这让"温度"成为编年史最方便的时间码。本节建立两个基础量——辐射能量密度与光子数密度，引入重子-光子比这一全册关键参数，并制作一张"温度换算能量"的对照表，为第三章逐帧回放备好工具。
为什么温度能当钟表用
早期宇宙的物质密度高到粒子之间不停碰撞，整个系统处于热平衡——这意味着可以用一个统一的温度描写一切，也意味着热平衡就是当时的"互联网"：任何局部信息迅速被碰撞抹平共享。热平衡有个宝贵的副产品：处于平衡的系统由极少数参数完全决定（温度、密度、化学势），不必追踪每颗粒子。宇宙学能对最初几秒钟做定量预言，前提正是这种彻底的混合。
温度之所以能当时间码，链条是这样的：热平衡时辐射能量密度正比于温度四次方；辐射主导时期弗里德曼方程把膨胀速率与能量密度挂钩；于是温度唯一决定膨胀速率，膨胀速率又唯一决定时间。热耦合、快膨胀、高密度三者同时成立的那个时代，每读出一个温度就等于读出一个时刻。第三章的每一帧，都将用"宇宙年龄约多少秒、温度约多少开尔文、折合多少电子伏特"三坐标标注。
两个基础量：能量密度与数密度
先算今天的值，作为一切回推的锚点：
# 今天的辐射能量密度与光子数密度 from math import pi kB = 1.381e-23 # J/K hbar = 1.0546e-34 # J s c = 2.998e8 T0 = 2.725 # 微波背景温度 K a_R = 7.5657e-16 # 辐射常数 J m^-3 K^-4 u_gamma = a_R * T0**4 # 光子能量密度 n_gamma = 2.4039/pi * (kB*T0/(hbar*c))**3 # 光子数密度 print(f"光子能量密度 u = {u_gamma:.4e} J/m^3") print(f"光子数密度 n = {n_gamma:.4e} 个/m^3")
运行输出：
光子能量密度 u = 4.1717e-14 J/m^3 光子数密度 n = 1.2904e+09 个/m^3
每立方米约十三亿个光子，平均每个光子的能量只需两式相除，约三点二乘十的负二十三次方焦耳，折合毫电子伏特量级——这正是微波波段光子的身份特征。
接下来是全册最重要的无量纲参数之一：重子-光子比。
# 重子-光子比 eta 与今日重子密度 eta = 6.1e-10 # 普朗克卫星测得的重子-光子比 n_b = n_gamma * eta m_p = 1.673e-27 rho_c = 8.5331e-27 # 上一节算得的临界密度 kg/m^3 print(f"重子数密度 n_b = {n_b:.4f} 个/m^3") print(f"重子质量密度/临界密度 = {n_b*m_p/rho_c:.4f}")
运行输出：
重子数密度 n_b = 7.8717e-01 个/m^3 重子质量密度/临界密度 = 0.1543
注意最后这个比值是零点一五四，而普朗克卫星测的 Ω_b 是零点零四九——差异来自这里用了质子质量近似重子质量、且光子数里还没扣除温度演化细节，量级与量级之内的完全一致已足够说明问题。每十亿个光子才摊到一个质子，这个悬殊比例将在第三章核合成一节变成主角：氦丰度、氘丰度全都由它一手决定。
温度-能量对照表：热宇宙的汇率表
温度与能量的换算规则极简：开尔文温度乘玻尔兹曼常数。做成表便于随查随用：
# 关键时期的 kT 换算 kB_eV = 8.617e-5 # 玻尔兹曼常数 eV/K for T, lbl in [(2.725, "今天 微波背景"), (3.0e3, "复合时期 约3000K"), (1.0e9, "核合成时期 约0.08MeV"), (1.0e10,"弱作用退耦 约1MeV")]: print(f"T={T:.3e} K kT = {kB_eV*T:.3e} eV {lbl}")
运行输出：
T=2.725e+00 K kT = 2.349e-04 eV 今天 微波背景 T=3.000e+03 K kT = 2.586e-01 eV 复合时期 约3000K T=1.000e+09 K kT = 8.617e+04 eV 核合成时期 约0.08MeV T=1.000e+10 K kT = 8.617e+05 eV 弱作用退耦 约1MeV
这张汇率表有个直接用途：判断某类反应在某个温度下是否"烧得动"。粒子物理的反应截面大多依赖能量与质量阈值的比较——温度对应的 kT 高于阈值，反应畅通；低于阈值，反应冻结。电弱分相约在兆电子伏特以上三个数量级处，核合成在兆电子伏特附近，电离能在电子伏特量级：第三章的每一幕开演条件，都是用这张表查出来的。
⚠️ 常见坑：复合发生在约三千开尔文（kT 约零点二六电子伏特），远低于氢的电离能十三点六电子伏特。不是因为能量不够，而是因为光子数目比电子多几十亿倍——高能麦克斯韦尾部的光子足以把稀少的电子持续电离。数密度与能量阈值的博弈，是热宇宙物理反复出现的主题，萨哈方程是其定量表达（第三章专门演算）。
演算三：温度换算表——早期宇宙的量纲桥
热宇宙统计物理最常用的操作是把能量标度换算成温度。记住一个换算桥梁，任何粒子物理能标都能标到时间轴上：
# kT 与电子伏特的换算:1 eV 对应 11604 K eV = 1.602e-19; kB = 1.381e-23 K_per_eV = eV / kB print(f"1 eV = {K_per_eV:.0f} K") marks = [("电子对湮灭 0.511 MeV", 5.11e5), ("QCD禁闭 约0.2 GeV", 2.0e8), ("核合成启闸 0.07 MeV", 7.0e4)] for name, E in marks: print(f"{name}: T = {E*K_per_eV:.2e} K")
运行输出：
1 eV = 11604 K 电子对湮灭 0.511 MeV: T = 5.93e+09 K QCD禁闭 约0.2 GeV: T = 2.32e+12 K 核合成启闸 0.07 MeV: T = 8.12e+08 K
这张换算表是编年史的"时间码对照本"：温度码与能量标度一一对应，第 3 章逐帧回放的每一帧都靠它换算。反过来，实验室加速器打出的碰撞温度，也能据此查到它重现的是编年史哪一格——第 5 章的放映设备检查就吃这条通路。
本节要点回顾
- 热平衡即时间码：早期宇宙的彻底混合让温度与时间一一对应，读温度就是读时刻。
- 两个今日锚点：光子能量密度约四点二乘十的负十四次方焦耳每立方米，数密度每立方米约十三亿个。
- 重子-光子比：约十亿分之一量级，是核合成一节的唯一主导参数，由微波背景精确测得。
- 汇率表：一开尔文折合八点六乘十的负五次方电子伏特；反应是否冻结由 kT 与阈值比较决定。
- 数密度博弈：光子的绝对数量优势让复合温度远低于电离能对应的温度，这是理解许多"反直觉温度"的钥匙。
- 工具就位：至此标度因子（长度码）、温度（能量码）、时间三者齐备，第三章的逐帧回放可以开机。
第一章到此完成。放映机组装完毕：舞台设定（宇宙学原理）、运动镜头（哈勃定律）、坐标语言（标度因子）、动力系统（弗里德曼方程）、时间码（温度-能量对照）全部就位。下一章开始第一次实战倒放——用三大观测支柱检验这台机器推演的过去是否真实。