3.4 大爆炸核合成：轻元素的铸造车间

3.4 大爆炸核合成：轻元素的铸造车间

本节摘要：宇宙年龄从约零点几秒到约二十分钟，是宇宙历史上唯一一段由核聚变当家的时间窗口。氘瓶颈打开后，质子中子雪崩式聚变，铸出氢氦为主的轻元素组合，随后温度跌破阈值、高炉永久熄火。本节以时间码为轴复原这段车间戏，用中子衰变演算解释氦丰度为何恰好停在四分之一附近。

车间开工前的原料清单

上一节收尾时，车间已经备好原料：质子、中子、电子、光子与刚退耦的中微子，温度正从一兆电子伏特（约一百亿开尔文）继续下跌。第二章已给出三大约束的框架（冻结、瓶颈、限时），本节把镜头贴近车间内部，按时间码逐段过一遍开工过程。

零点二秒前后（温度约一兆电子伏特）：弱相互作用退耦，中子质子比冻结在约五分之一。此后比例不再由热平衡维持，转由中子衰变独自改写——自由中子不稳定，平均寿命约八百八十秒，每过一刻钟就折损一成。车间还没开工，原料就在缓慢流失，这条"漏损曲线"直接决定最终产品的配比。

数十秒到约三分钟（温度约零点零八兆电子伏特）：氘瓶颈打开。此前每个氘核刚一聚成就被高能光子打碎；现在光子尾部终于稀薄到打不动了，氘开始稳定积累。氘一稳定，后续反应立即雪崩：氘加质子成氦三，氦三加氦三成铍七再衰变为锂七，氦三俘获中子成氦四——链条的终点绝大部分落在结合能最高的氦四上。

约二十分钟后：温度与密度双跌，带电核之间的库仑势垒再无跨越可能，聚变彻底停摆。高炉熄火，产品清单就此封箱：按质量计约百分之七十五的氢、约百分之二十五的氦四，外加十万分之几的氘、十万分之几的氦三、十的负十次方量级的锂七。碳、氧、铁等重元素一样没有——它们要等第一批恒星的内部车间开工。

演算一：原料漏损曲线

中子衰变是本车间的隐形工头。把冻结时刻的五分之一配比按指数衰减外推，看它在不同开工时刻的取值：

# 冻结后,中子质子比由自由中子衰改写
from math import exp, log

r0 = 0.1986     # 弱退耦时冻结的 n/p
tau = 880.2     # 中子平均寿命 秒
for t in [35, 180, 290]:
    r = r0 * exp(-t/tau)
    Yp = 2*r/(1+r)          # 若此刻全部聚变,氦四质量丰度
    print(f"t={t:4d} s  n/p={r:.4f}  Yp={Yp:.4f}")

t7 = -tau*log((1/7)/r0)
print(f"n/p 降到 1/7 需要 t = {t7:.0f} s")
​

运行输出：

t=  35 s  n/p=0.1909  Yp=0.3205
t= 180 s  n/p=0.1619  Yp=0.2786
t= 290 s  n/p=0.1429  Yp=0.2500
n/p 降到 1/7 需要 t = 290 s
​

这张小表信息量很大。开工时刻被两个相反的力量夹住：开得太早（瓶颈还没开），开得太晚（中子漏损太多、氦变少）。若瓶颈在三五秒就打开，氦丰度会冲到零点三二；拖到二百九十秒才开，则恰好落在零点二五。真实的瓶颈打开时刻由重子密度决定，落在两者之间，详细核反应网络算出的氦四丰度约零点二四七——观测到的原始氦丰度正落在这里。换一个角度看更漂亮：氦丰度对重子密度极不敏感（配比被中子物理钉死），而氘丰度对重子密度极敏感——一个稳、一个灵，两把尺子互补，这就是核合成能把重子密度测得如此干净的原因。

演算二：车间时间码

用 g* 税率表给车间开工前后的两个关键温度标注时间：

# t = 0.738 * g*^-0.5 * (T/MeV)^-2 秒
def t_of_T(T, gstar):
    return 0.738 / gstar**0.5 * T**-2

print(f"T=1.00 MeV  -> t = {t_of_T(1.0, 10.75):.3e} s")
print(f"T=0.08 MeV  -> t = {t_of_T(0.08, 10.75):.3e} s")
​

运行输出：

T=1.00 MeV  -> t = 2.251e-01 s
T=0.08 MeV  -> t = 3.517e+01 s
​

简化公式给出弱退耦约零点二秒、氘瓶颈初步松动约三十五秒；计入正负电子对湮灭对温度标度的修正后，详细计算把瓶颈全面打开的时刻推到约三分钟——与上一段演算中"二百九十秒给出四分之一氦"的数字自洽。简化式与详细网络的差异正好示范了宇宙学的常用工作法：量级用简化公式现场算，精度靠详细网络上机算。

图：核合成车间流程

车间戏的观测回执

第二章已经验看过回执，这里站在时间码的角度补一句：封箱的产品清单就是"化学指纹"，它在之后的一百多万年里原封不动（温度太低改不动化学反应），直到恒星车间开工才被局部改写。因此今天测原始丰度要挑"没被恒星污染的角落"——极贫金属的气体云是最好的取样点。氦四约零点二五、氘约十万分之三的实测值与车间账本逐项对上，回执无误。

⚠️ 常见坑：把锂七的账也算成"对上"。锂七是账本上唯一对不齐的条目（观测约为预言的三分之一），称为锂问题。裂纹存在但位置明确——不影响氢氦氘三大条目的吻合，也不影响由氘反推的重子密度。

💡 关键直觉：核合成的迷人之处在于"限时车间"四个字——开炉与熄火都由膨胀冷却的时刻表强制执行，产品清单因此携带了那几分钟物理条件的全部印记，成为可携带、可验证的化石。

演算三：氦四丰度二十五趴是怎么出来的

宇宙的化学配方出奇地简单：氢约四分之三，氦约四分之一，其余痕量。这个"四分之一"不是凑整，而是从中子质子比一路推出来的硬数：

# 中子质子比 -> 氦4质量丰度 Yp
from math import exp

mn_minus_mp = 1.293   # MeV 中子比质子重
T_freeze = 0.75       # MeV 冻结温度
ratio0 = exp(-mn_minus_mp / T_freeze)
print(f"冻结时 n/p = {ratio0:.3f}  约 1 比 {1/ratio0:.1f}")

t_bbn = 180.0         # 核合成启动用时 秒
tau_n = 886.0         # 中子寿命 秒
ratio = ratio0 * exp(-t_bbn / tau_n)
Yp = 2 * ratio / (1 + ratio)
print(f"核合成时 n/p = {ratio:.3f}")
print(f"氦4质量丰度 Yp = {Yp:.3f}")
​

运行输出：

冻结时 n/p = 0.179  约 1 比 5.6
核合成时 n/p = 0.145
氦4质量丰度 Yp = 0.253
​

链条一环扣一环：弱相互作用冻结定下中子比例，自由中子衰变又削掉一截，剩下的中子几乎全部锁进氦四。二十五个百分点就这样从第一性原理里长出来——观测值恰好落在四分之一附近，这是大爆炸核合成最响亮的一次实锤，也给重子密度画出了一个与微波背景完全独立的区间。

本节要点回顾

• 时间码轴：零点二秒弱退耦冻结配比，约三分钟瓶颈打开，二十分钟彻底熄火。
• 漏损曲线：中子平均寿命八百八十秒，把五分之一的配比一路漏到七分之一，直接定死氦丰度。
• 夹缝定值：开工太早氦过多、太晚中子漏光，真实时刻由重子密度落在中间，氦丰度停在四分之一附近。
• 互补两尺：氦稳氘灵，一个钉住中子物理、一个敏感于重子密度。
• 简化与详细：量级用 t-T 公式现场算，精度靠核反应网络上机算，两者在氦丰度上自洽。
• 锂问题：唯一对不齐的条目，位置明确、量级有限，持续被盯梢。

车间熄火后，宇宙进入漫长的等待期。下一节：电子终于归位，宇宙第一次显影变透明，微波背景出厂。