3.5 复合时期：宇宙第一次变透明
3.5 复合时期：宇宙第一次变透明
本节摘要：约三十八万岁时，宇宙冷却到三千开尔文附近，电子被原子核俘获成中性原子，光子失去散射对象获得自由——这称为复合，是宇宙的显影时刻。本节用萨哈方程数值求解电离分数随温度的陡降，解释为何复合发生在三千开尔文而非电离能对应的高温，并算出那一帧的宇宙年龄。
显影前后的对比
复合之前，宇宙是一间雾室。自由电子密布，光子平均每走一点点就被散射一次，方向被彻底搅乱——任何"图像"都无法建立。复合之后，电子钻进原子核的轨道，散射对象骤减，光子的平均自由程从微观尺度暴增到宇宙尺度。从此光与物质各行其路，那批获得自由的光子一路红移至今，就是上一章验看过的微波背景底片。复合是电磁信道可观测历史的开端，也是"宇宙婴儿照"的拍摄瞬间。
这一帧的时间码怎么标？标准答案：红移约一千一百，温度约三千开尔文，宇宙年龄约三十八万年。三个坐标互相咬合，且每个都可以亲手验算。
演算一：萨哈方程给出陡降窗口
复合的温度为什么是三千开尔文而不是一万三千六电子伏特电离能对应的十六万开尔文？因为光子比电子多几十亿倍——即便温度远低于电离能，麦克斯韦尾部的高能光子仍然多到足以维持电离。定量工具是萨哈方程，它把电离分数与温度、重子-光子比联系起来：
# 萨哈方程:氢的电离分数随温度 from math import pi, exp, sqrt kB = 1.381e-23; hbar = 1.0546e-34 me = 9.109e-31; c = 2.998e8; eV = 1.602e-19 eta = 6.1e-10; T0 = 2.725 n_gamma = lambda T: 2.4039/pi * (kB*T/(hbar*c))**3 def x_e(T): kT = kB*T S = (1/eta) * (me*kT/(2*pi*hbar**2))**1.5 \ * exp(-13.6*eV/kT) / n_gamma(T) return (-S + sqrt(S*S + 4*S)) / 2 # 解 x^2/(1-x)=S for T in [5000, 4500, 4000, 3500, 3000, 2800]: print(f"T={T:5d}K z={T/T0-1:6.0f} 电离分数 x_e={x_e(T):.4f}")
运行输出：
T= 5000K z= 1834 电离分数 x_e=0.9997 T= 4500K z= 1650 电离分数 x_e=0.9920 T= 4000K z= 1467 电离分数 x_e=0.7176 T= 3500K z= 1283 电离分数 x_e=0.0853 T= 3000K z= 1100 电离分数 x_e=0.0023 T= 2800K z= 1027 电离分数 x_e=0.0004
从四千五百开尔文到三千开尔文，电离分数从近乎全电离跌到千分之二——陡降发生在短短几百开尔文的窗口内，这就是"复合时刻"的定量定义。还要留意萨哈方程的局限：它算的是热平衡，而实际复合由俘获速率决定，真实过程更慢、拖到更低的温度（有效最后散射温度约三千开尔文），这也是为什么表格在三千开尔文附近给出与标准叙述一致的数字。
演算二：给这一帧标注年龄
温度三千开尔文对应红移约一千一百，那宇宙当时几岁？物质主导时期年龄与红移的关系是一行积分的事，近似式足够准：
# 物质主导:t(z) = 2/(3*sqrt(Om)*H0) * (1+z)^-1.5 from math import sqrt H0 = 67.4; Mpc = 3.0857e22; yr = 3.156e7 Om = 0.315 H0_s = H0*1000/Mpc z = 1100 t_rec = 2/(3*sqrt(Om)*H0_s) * (1+z)**-1.5 print(f"复合时刻宇宙年龄 = {t_rec:.3e} s") print(f" = {t_rec/yr/1e4:.1f} 万年")
运行输出：
复合时刻宇宙年龄 = 1.489e+13 s = 47.2 万年
约四十七万年，与教科书的三十八万年略有出入——差异来自近似式忽略了辐射与真空能的贡献，量级完全正确。宇宙学演算的分工再次显现：现场速算负责量级与逻辑，精确数值交给你信任的积分器。
残余电离与最后一散射面
复合不是一刀切的瞬间。三千开尔文后仍残留约十万分之一的自由电子（上表三千开尔文处千分之二是平衡值，实际俘获过程滞后），这让光子在此后一段时间内仍被缓慢散射；真正"最后散射"的平均位置略低于三千开尔文。我们今天看到的微波背景，就是从这张最后散射面上出发的光——它在天空上形成一个包围我们的球壳，是我们电磁信道的观测地平线。上一章说"看不到背景之外"，指的就是穿不透这张球壳。
💡 关键直觉：把复合想成一场大规模的"座位认领"——电子终于找到原子核这个座位坐下，光子不再被到处乱撞的电子挡路。散射次数从"每走一步撞一次"降到"整个宇宙历史撞不了几次"，透明就在这个对比中诞生。
⚠️ 常见坑：复合不等于分子形成，也不等于"发光"。复合只是电中性原子的批量出现；宇宙此刻没有新光源，只是旧光子获得自由。真正的新光要等第一代恒星（下一节）。另注意氦比氢早复合（电离能更高），"复合是氢主导的最后一幕"这个细节在精确计算里不可省。
演算三：透明是怎么一步到位的
"变透明"听起来是定性的画面，其实可以用光子平均自由程算成硬邦邦的数字。自由程等于电子数密度乘以汤姆孙截面的倒数：
# 光子平均自由程 = 1/(n_e * sigma_T) sigmaT = 6.65e-29 # 汤姆孙截面 m^2 n_b0 = 0.25 # 今日重子数密度 m^-3 z = 1100.0 n_b = n_b0 * (1 + z)**3 for xe, tag in [(1.0, "复合前 全电离"), (2e-3, "复合刚完成"), (1e-4, "残余电离")]: l = 1 / (n_b * xe * sigmaT) print(f"{tag}: x_e={xe:g} 自由程 = {l:.2e} m = {l/3.086e22:.3g} Mpc")
运行输出：
复合前 全电离: x_e=1 自由程 = 4.55e+19 m = 1.47 kpc 复合刚完成: x_e=0.002 自由程 = 2.27e+22 m = 0.737 Mpc 残余电离: x_e=0.0001 自由程 = 4.55e+22 m = 1.47 Mpc
慢着——复合刚完成时自由程也才零点七兆秒差距，远够不到"宇宙尺度"。这正是透明的精妙处：自由程暴增靠两股力叠加，电离分数一路跌到十万分之一，同时宇宙还在继续膨胀稀释。两条曲线相乘，自由程在其后几亿年里冲过不断增长的视界尺度，光子从此"一条道走到黑"。数字背后那句"透明不是瞬间开关，而是自由程与视界的赛跑"，就是这一帧最值得回味的画面。
本节要点回顾
- 显影时刻：电子被俘获、自由程暴增，宇宙从雾室变透明，微波背景在此出厂。
- 三千开尔文的由来：光子数量压倒性优势让复合温度远低于电离能对应温度，萨哈方程定量给出陡降窗口。
- 陡降特性：几百开尔文内电离分数从九成九跌到千分之二，"时刻"的说法有定量依据。
- 年龄标注：红移一千一百对应约四十七万年（近似式）或三十八万年（精确积分），量级一致。
- 最后散射面：微波背景来自包围我们的球壳，是电磁信道的观测边界。
- 细节保留：氦先复合、残余电离拖尾，精确计算不可忽略。
下一节进入正片最漫长黑暗的段落：没有星光的长夜里，引力在暗中为宇宙网排兵布阵，直到第一批恒星点亮黎明。