5.2 星系形成的脚手架：从晕到发光体
5.2 星系形成的脚手架：从晕到发光体
本节摘要：引力不稳定性不挑食——越密的区域越早越过坍缩线，小扰动先成小晕，小晕合并成大晕，结构按"先小后大"的次序等级式生长，与热暗物质预言的"先大后小"相反，这是冷暗物质模型的早期关键胜绩。本节计算三类质量晕的维里半径与特征速度、团级晕的维里温度与自由落时标，把"星系为什么住在这种晕里、星系团为什么发不了光"讲成定量问题。
坍缩的次序：先小后大
线性增长越过 δ≈1 后，扰动区脱离膨胀、转入坍缩，弛豫成维里化的束缚体——暗物质晕。次序由密度门槛决定：坍缩红移 z_col 正比于该尺度的初始扰动幅度，而初始功率谱在小尺度上幅度更高（加上小尺度更早进入视界吃过辐射期的亏之后，净效果仍保持小幅领先），于是小晕先成，大晕后成。几十亿年里，晕像滚雪球一样并合：十亿级的晕并成百亿级，百亿级并成万亿级的星系团，暗物质骨架上逐渐挂满发光的星系。
这次序是判决性的。热暗物质（如 30 eV 中微子）自由流到超星系团尺度才停下，只能"先大后小"——先形成薄饼状超结构再碎裂成星系。1980 年代的星系巡天与红移测量已经能分辨成团次序：观测到的星系团里有年老星系、星系里没有比星系团更年轻的层次结构，"先小后大"胜诉，热暗物质出局。这是暗物质"冷"属性从画像到判例的完整闭环。
动手验算之一：三类晕的身材表
晕的身材由一个质量数决定。用 3.3 节的定义（平均密度 200 倍临界密度）反解半径，再得特征速度：
import math G, MS, KPC, MPC = 6.674e-11, 1.989e30, 3.0857e19, 3.0857e22 H0 = 67.36e3 / MPC rho_c = 3 * H0 ** 2 / (8 * math.pi * G) def r200(M): # M 单位 Msun，返回米 return (3 * M * MS / (4 * math.pi * 200 * rho_c)) ** (1 / 3) print("（M200，r200，v200）：") for M in [1e10, 1e12, 1e15]: r2 = r200(M) v2 = math.sqrt(G * M * MS / r2) print(" M=%.0e Msun：r200=%.0f kpc，v200=%.0f km/s" % (M, r2 / KPC, v2 / 1e3))
输出：
（M200，r200，v200）： M=1e+10 Msun：r200=46 kpc，v200=31 km/s M=1e+12 Msun：r200=212 kpc，v200=143 km/s M=1e+15 Msun：r200=2117 kpc，v200=1426 km/s
矮星系晕三万公里每秒的十到三十倍……准确说：矮星系级 31 km/s、银河系级 143 km/s、星系团级 1426 km/s——半径与速度都按质量的立方根标定（r ∝ M^⅓，v ∝ M^⅓）。这张身材表与观测直接对表：矮星系测到的转速约 10 到 50 km/s，银河系 200 出头（盘内含重子略高），星系团速度弥散上千——模型尺寸合适。
动手验算之二：温度与时标，谁发不了光
晕维里化后，掉进去的气体会被激波加热到维里温度，能不能发光、能不能冷却成恒星，取决于温度与冷却时标的赛跑。
import math G, MS, MPC, KPC, YR = 6.674e-11, 1.989e30, 3.0857e22, 3.0857e19, 3.156e7 k_B, KEV = 1.381e-23, 1.602e-16 H0 = 67.36e3 / MPC # 接上块的常数与函数 rho_c = 3 * H0 ** 2 / (8 * math.pi * G) def r200(M): return (3 * M * MS / (4 * math.pi * 200 * rho_c)) ** (1 / 3) v_cl = math.sqrt(G * 1e15 * MS / r200(1e15)) T_cl = 0.6 * 1.673e-27 * v_cl ** 2 / (3 * k_B) print("团级晕：v=%.0f km/s -> T=%.1e K（kT=%.1f keV，与富团X射线同量级）" % (v_cl / 1e3, T_cl, T_cl * k_B / KEV)) rho_h = 1e12 * MS / (4 * math.pi / 3 * r200(1e12) ** 3) t_ff = math.sqrt(3 * math.pi / (32 * G * rho_h)) print("1e12 晕平均密度对应的自由落时标 t_ff = %.1f Gyr" % (t_ff / YR / 1e9))
输出：
团级晕：v=1426 km/s -> T=4.9e+07 K（kT=4.2 keV，与富团X射线同量级） 1e12 晕平均密度对应的自由落时标 t_ff = 1.6 Gyr
第一笔账接上第2章：星系团晕的维里温度正落在 X 射线波段的 keV 区间——X 射线卫星看到的团内热气体，就是掉进团级势阱被激波加热的重子，账目严丝合缝。第二笔账是星系形成的节拍器：自由落时标十亿年上下，冷却时标若比它短，气体一层层冷掉、恒星点亮，星系长成；若比它长，气体就一直泡在热相——星系团正是这种"冷却破产"的大晕，所以团里星系占比极小，大部分重子以百万度的 X 射线气体形式终老。半人马团、后发团的弥散 X 射线辉光，就是没冷却下来的重子库。
脚手架与住户：晕-星系的匹配
暗物质晕是脚手架，星系是住户，两者的匹配规则叫晕占领模型：晕要够大才点得起灯（气体冷却成恒星的效率在 10¹² Msun 附近达到峰值），太大反而点不亮（冷却破产）。于是有了今日的图景：每颗亮星系背后有一个晕，但晕不都住着亮星系——大晕住一窝（星系团住几百个）、中晕住一两个、最小的晕（10⁸ 到 10¹⁰ Msun）基本空置，只有超暗矮星系级别的住户。"有房无人"与"一房多户"并存，这个分配规则把旋转曲线、卫星计数、成团统计串成一张网，也是 5.4 节"卫星缺失"疑点的定量舞台。
💡 关键直觉：星系质量函数（星系按质量的数目分布）比晕质量函数陡得多——不是所有晕都发光，重子冷却这道闸门决定了"住户规则"，暗物质只管盖房。
变式：并合仍在进行
等级成团不是历史剧，是进行时。银河系正在吞并人马座矮星系，未来与仙女座的并合将把本星系群滚成一个椭圆星系；LIGO 听到的黑洞并合是恒星级"脚手架"的余响。观测并合率与模拟对比，是 ΛCDM 动力学剧本的又一重验尸级检验——目前对得上。
本节要点回顾
- 次序即判例：先小后大的等级成团与观测一致，热暗物质的先大后小出局。
- 身材立方根律：10¹⁰/10¹²/10¹⁵ Msun 晕的半径 46/212/2117 kpc，速度 31/143/1426 km/s。
- 冷却是闸门：团级晕维里温度数千万度（keV 级 X 射线），冷却破产所以不点灯；星系级晕自由落时标约 1.6 Gyr，来得及冷却。
- 住户规则：晕占领模型串起星系-晕匹配，"空置小晕"是下一节疑点的伏笔。
脚手架讲完了，下一节进机房——看十亿粒子的引力 N 体模拟怎么把以上一切真正演出来。