5.1 密度扰动与增长：暗物质的先发优势

5.1 密度扰动与增长：暗物质的先发优势

本节摘要：今天的星系要从十万分之一的原初涨落长出来，增长纪元的账目必须精确。辐射时代，光子-重子流体的扰动被压成驻波（对数级增长），纯重子宇宙到复合时攒不出足够的结构。暗物质不与光子耦合，从进入视界起持续积累，"先发优势"三个量级以上；复合后重子落进已就位的暗物质势阱，一夜之间补齐进度。本节用完整演算链量化这段历史，并算出重子 Jeans 质量这道门槛。

增长的两条交规

密度扰动 δ = δρ/ρ 的增长速率由两个纪元的两条"交规"决定。辐射主导期（z ≳ 3400）：宇宙膨胀太快，引力增长追不上，视界内扰动只能对数级缓增（Mészáros 效应）——更糟的是，重子与光子耦合成的流体在压缩时辐射压反弹，扰动呈声波式驻振荡，净增长几乎为零。物质主导期（z ≲ 3400）：引力赢得膨胀，线性扰动 δ ∝ a（尺度因子）增长，指数级追赶。

复合（z≈1090）是重子的解禁日。此后重子脱离光子，立刻掉进早已就位的暗物质势阱——CMB 测到的重子扰动幅度只有 10⁻⁵ 量级，若没有暗物质预先挖好的阱，这点涨落到今天也只长到 δ≈0.03，连一个星系都坍不出来。先发优势不是修辞，是结构存在的先决条件。

图 两条增长曲线的分岔

动手验算：先发优势的完整链条

把一条星系尺度扰动的成长史编成程序。设它在红移 10⁸ 进入视界，初幅 3×10⁻⁵；辐射期对数增长因子为 (1 + ln z_ent/z_eq)；物质期按尺度因子线性放大。

import math
d0, z_ent, z_eq, z_rec = 3e-5, 1e8, 3400, 1090
grow_rad = 1 + math.log(z_ent / z_eq)     # 辐射期对数增长
d_eq = d0 * grow_rad
d_rec = d_eq * z_eq / z_rec               # 物质期线性放大至复合
d_today = d_rec * z_rec                   # 复合至今继续线性
print("CDM 增长链：初幅 d=%.0e（视界进入时）" % d0)
print("  辐射期对数增长 ×%.1f -> d(z_eq)=%.1e" % (grow_rad, d_eq))
print("  物质期至复合 ×%.2f -> d(z_rec)=%.1e" % (z_eq / z_rec, d_rec))
print("  至今 ×%d -> d~%.2f（已进入非线性成团）" % (z_rec, d_today))
print("  对照：纯重子 d(z_rec)=%.0e，至今仅 %.3f -> 无法成团"
      % (d0, d0 * z_rec))
​

输出：

CDM 增长链：初幅 d=3e-05（视界进入时）
  辐射期对数增长 ×11.3 -> d(z_eq)=3.4e-04
  物质期至复合 ×3.12 -> d(z_rec)=1.1e-03
  至今 ×1090 -> d~1.15（已进入非线性成团）
  对照：纯重子 d(z_rec)=3e-05，至今仅 0.033 -> 无法成团
​

账目一目了然：同样的初始涨落，带上暗物质到今天长成 δ≈1（坍缩成结构），不带则停在 0.03（永远的稀薄纹路）。先发优势在复合时刻约为 36 倍（1.1×10⁻³ 对 3×10⁻⁵），这份差距线性放大到今天。顺带解释了一个 CMB 反直觉：声学峰十万分之一的涨落本该让宇宙至今平淡无奇，正因为看到的是"被押着的重子"，而真正的推手暗物质在 CMB 上只留暗记（引力势的早期等势演化），结构史才可能如此波澜。

动手验算：重子解禁日的门槛

复合后的重子也不是立刻就能任意聚团——它自带一道 Jeans 门槛：尺度小于声学视界对应的 Jeans 质量，气体压强顶住引力。拿复合后早期（z=1000、声速 5 km/s 的中性气体）算这道门槛：

import math
G, MS, KPC = 6.674e-11, 1.989e30, 3.0857e19
rho_c = 8.52e-27
om, z, cs = 0.3138, 1000, 5e3
rho_m = om * rho_c * (1 + z) ** 3
lam_J = cs * math.sqrt(math.pi / (G * rho_m))       # Jeans 波长（物理）
M_J = rho_m * (4 * math.pi / 3) * (lam_J / 2) ** 3 / MS
print("z=%d，c_s=%.0f km/s：" % (z, cs / 1e3))
print("  lam_J = %.1f kpc（共动），M_J = %.1e Msun" % (lam_J / KPC * 1e3, M_J))
​

输出：

z=1000，c_s=5 km/s：
  lam_J = 21.5 kpc（共动），M_J = 2.1e+05 Msun
​

Jeans 质量约两万倍太阳质量——比球状星团略小、比任何星系都小。也就是说重子解禁后，凡是有暗物质势阱的地方它都能掉进去，门槛低到不构成障碍；反过来，这个尺度以下的"微结构"只能靠暗物质单独支撑（这批最小暗物质微晕至今仍是理论预言的悬赏目标）。

解读：先发优势的观测回执

这套增长剧本不是自说自话，它有两张观测回执。其一，CMB 峰高：重子被押期间扰动横盘，CMB 上重子贡献的峰幅度偏低——若重子当年自由增长，峰高组合对不上。其二，星系红移分布：最远的星系已在 z≈10 到 14 被韦布望远镜抓到，若纯重子宇宙，那时的 δ 才 0.1 量级，坍缩无从谈起。数字环环相扣，先发优势是冷暗物质模型最不可替代的贡献。

⚠️ 常见坑：把"扰动增长"理解成所有尺度同步放大。视界内的尺度按上述规则走，超视界扰动由广义相对论性处理（绝热守恒），两套账不能混。本节的演算针对星系尺度（早进视界），对星系团以上尺度需改用完整传递函数。

两问两答

问：为什么辐射期暗物质也只对数缓增，而不是完全不涨？
答：它不受光子压强，但广义相对论不讲情面——辐射主导期宇宙膨胀太快（a ∝ t^½），引力拉拢物质的效率被膨胀稀释，亚视界扰动只能按 ln a 缓慢爬升。"对数增长"是暗物质与光子压强无关、却仍逃不脱膨胀稀释的双重结果；重子则是干脆被光子按住纹丝不动。

问：35 倍的先发优势听起来不多，值得这么大书特书吗？
答：倍数不吓人，"缺它就不成"才吓人。35 倍的差距发生在指数放大的起点——δ 终值 1.15 对 0.03，一个是坍缩成星系的门槛之上，一个是永远安静的纹路。用赛跑打比方：这不是领先 35 米，而是其中一个选手根本没被允许起跑。宇宙学里"差一个数量级"常常就是"存在与不存在"的差别。

本节要点回顾

• 两条交规：辐射期对数缓增（重子更是驻波横盘），物质期线性猛增。
• 先发优势量化：复合时 36 倍，今日 35 倍上下；没有它 δ 终值只有 0.03。
• Jeans 门槛低：z=1000 重子 Jeans 质量约 2×10⁵ Msun，远小于星系，势阱照单全收。
• 双重回执：CMB 峰高组合与高红移星系计数，都为这段增长史背书。

增长引擎就位。下一节看坍缩如何按质量分级排队——先小后大的等级成团，以及星系为何住进今天这样的晕里。