4.2 暗物质晕中的气体坍缩

4.2 暗物质晕中的气体坍缩

本节摘要：星系不是凭空长出来的，它们的"胚胎"是暗物质晕，而"养分"是落入晕中的气体。本节追踪气体从宇宙网中被捕获、经辐射冷却收缩、靠角动量守恒铺成旋转盘、最终在临界密度下点燃恒星的完整链条，并用 Python 估算冷却时间与盘面尺度。

气体如何落入暗物质晕

宇宙中百分之九十五以上的物质不会发光，但重子气体沿着暗物质编织的引力势阱流动。初始时刻，气体与暗物质几乎完全混合，都随宇宙膨胀而稀释。当某个暗物质晕从膨胀背景中脱离、完成非线性坍缩时，晕内引力势阱开始从周围宇宙网中"虹吸"气体。气体落入晕的物理与暗物质不同：暗物质粒子无碰撞，只被引力束缚，形成弥散的平衡分布；气体则可以通过辐射损失能量，因此会继续收缩、增密，甚至沉向晕中心。

气体被捕获的效率并非百分之百。激波加热把落向晕中心的气体加热到晕的维里温度——对质量 10¹² 太阳质量的晕，维里温度约 10⁶ 开尔文，足以电离氢并发出 X 射线。这种热气体在晕中形成静水压平衡的"晕内介质"（ICM 的前身）。只有温度降下来、气体才有机会继续收缩成星。因此"冷却"是气体坍缩故事的第二个关键词。

辐射冷却：把热能交给光子

气体降温靠辐射。对温度百万开尔文左右的热气体，主要冷却机制是韧致辐射（自由-自由辐射）：电子被离子偏转时加速，发射 X 射线光子。对温度一万到十万开尔文的气体，金属谱线的碰撞激发冷却变得重要——重元素（碳、氧、铁等）像"冷却剂"一样把热能迅速辐射出去。这意味着金属丰度直接影响冷却效率：宇宙早期的气体几乎不含金属，只能靠氢分子与韧致辐射缓慢冷却，这推迟了第一批星系的点火时间；到了金属富集的后期，冷却快得多。

冷却时间与自由落体时间的比较决定了气体的命运。若冷却时间远小于自由落体时间，气体在下落过程中就大量失能，会坍缩成中心冷盘；若冷却时间接近或超过自由落体时间，气体只能滞留在晕的热气体中。临界质量的量级约 10¹² 太阳质量：低于它，冷却主导，气体进得来、星形成；高于它，激波加热的气体难以冷却，落入率骤降。这一"冷却界限"直接解释了大质量星系为何恒星形成效率反而低——气体进不了冷却轨道。

图：气体落入晕的物理链条

角动量：从球形到盘形的分水岭

如果气体坍缩是完全球对称的，最终只会形成一个中心恒星团，而不会出现盘。观测却告诉我们，宇宙里一半的星系是旋转的盘。谜底在于角动量。气体从宇宙网落入晕时，潮汐力矩已经赋予它可观的角动量；坍缩过程中角动量守恒使旋转速度随半径收缩而增大，离心力最终阻止径向继续收缩，气体在垂直于角动量方向的平面上铺成薄盘。盘半径由角动量大小决定，典型值为几个千秒差距——这与观测到的盘星系尺度完全吻合。

角动量的去向是星系形成的经典难题。数值模拟早期曾预言，气体通过冷却落入中心时会因动量交换损失角动量，导致"角动量灾难"：盘应该比观测的小十倍。后来的高分辨率模拟发现，宇宙网中的气体主要沿纤维状结构带着大角动量流入，且气体在合并过程中不易充分混合，才避免了这场灾难。这一"冷流吸积"图景改变了我们对星系早期形成的理解：早期星系未必都经过热晕阶段，一部分气体直接以准直流的形式喂给星系。

Python 算例：冷却时间与盘尺度

import numpy as np

G = 4.30e-6   # kpc*(km/s)^2/Msun
kB = 1.38e-16 # 玻尔兹曼常数（尔格/K），示意性使用
mH = 1.67e-24 # 氢原子质量 g

def virial_temperature(M_vir, R_vir_kpc):
    # 维里温度：kT ~ G*M*mu*mH / (2*R)，mu 为平均分子量（离子化气体约 0.6）
    mu = 0.6
    T = 1e6 * (M_vir / 1e12) / (R_vir_kpc / 100.0)
    return T

def cooling_time(T, n_e):
    # 韧致辐射冷却率的量级估计：t_cool ~ (3/2) n kT / (L_ff)
    # 对 1e6 K、密度 1e-3 cm^-3 的典型值，t_cool 约 1e9 年
    n5 = n_e / 1e-3
    T6 = T / 1e6
    return 1.0e9 * (n5 * T6**0.5) ** -1   # 粗略近似（单位：年）

# 示例：10^12 Msun 的晕，半径 150 kpc，z=0
M = 1e12
R = 150.0
T_vir = virial_temperature(M, R)
print(f"维里温度 T_vir ≈ {T_vir:.2e} K（热到足以电离氢）")

t_cool = cooling_time(T_vir, 1e-3)
t_ff = 1.0e8 * (R / 100.0)**1.5      # 自由落体时间量级
print(f"冷却时间 ≈ {t_cool:.2e} 年，自由落体时间 ≈ {t_ff:.2e} 年")
print("冷却时间 < 自由落体时间，气体可有效冷却坍缩 → 形成冷盘")
​

从坍缩到点火

冷却后的气体不再均匀下落，而是碎裂成分子云。当某团气体的质量超过金斯质量（Jeans 质量），引力的自压缩超过热压支撑，云团就开始自由坍缩。分子云内部致密的核区进一步收缩，密度达到每立方厘米十万到百万个粒子时，中心点火核聚变——一颗恒星诞生。星系盘的恒星形成不是持续恒定的，而是受分子云供给、旋臂压缩与反馈扰动调制，形成第 4.4 节将要讨论的"恒星形成主序"。

坍缩链条对星系性质的决定作用

回头看，气体坍缩链条的每一步都在给星系"盖章"。冷却效率决定恒星形成效率与质量-光度比；角动量决定盘尺度和形态；金斯质量决定恒星形成的最小尺度与初始质量函数；而反馈（第 4.5 节）决定坍缩何时停、冷却何时止。正是这些环节的耦合，才出现了星系性质与暗物质晕质量的复杂对应：矮星系恒星形成效率低（反馈强）、大质量星系恒星形成效率低（冷却差）、中等质量盘星系效率最高——一条倒 U 形的"恒星形成效率-质量"曲线，构成了当代星系形成理论最核心的经验约束之一。

本节的坍缩链条将直接引出下一节的问题：星系不是单次坍缩的产物，而是无数小"胚胎"反复并合的组装体。从单体生长到群体组装，正是层次性并合模型的主题。