5.1 密度扰动与宇宙结构形成
5.1 密度扰动与宇宙结构形成
本节摘要：宇宙中所有宏大的结构——星系、星系团、宇宙网——都起源于暴胀时期微小的量子涨落。这些十万分之一量级的原初密度扰动在引力不稳定性驱动下逐渐增长，最终坍缩形成我们今天看到的天体。本节分析引力不稳定性的线性增长理论，区分暗物质和重子物质在结构形成中的不同角色，并讨论从线性扰动到非线性坍缩的转变。
学习目标
阅读完本节，你应当能够：
- 描述原初密度扰动功率谱的基本特征
- 推导引力不稳定性的线性增长因子
- 解释暗物质为何能在复合之前开始坍缩而重子不能
- 分析转移函数如何将原初功率谱转化为观测到的物质功率谱
- 理解从线性增长到非线性坍缩（球对称坍缩模型）的转变
从十万分之一到星系
CMB温度涨落仅为十万分之一，这意味着大爆炸后38万年时宇宙的物质分布极其均匀。但正是这些微小的密度不均匀性——密度略高的区域和密度略低的区域——在引力的驱动下逐渐放大，最终成长为今天的宇宙结构。
理解这一过程的关键是引力不稳定性：在均匀的物质分布中，密度略高的区域引力更强，会吸引周围物质向其汇聚，使密度对比度进一步增大。这就是金斯不稳定性——当扰动尺度大于金斯长度时，引力克服压强支撑，扰动增长。
原初功率谱
暴胀理论预言原初密度扰动几乎是标度不变的——不同尺度上的扰动幅度几乎相同。用功率谱 P(k) \propto k^{n_s} 来描述，其中 n_s 是谱指数。慢滚暴胀预言 n_s 略小于1（约0.965），与Planck卫星的精确测量高度一致。
原初功率谱是结构形成的"初始条件"。但观测到的物质功率谱并非标度不变的——它受到宇宙演化过程中多种物理过程的调制。
暗物质的"脚手架"作用
暗物质在结构形成中扮演着不可替代的"脚手架"角色。
在复合之前（38万年之前），宇宙中充满了光子，光子的压强会抵抗引力的坍缩效应。重子物质（电子和质子）与光子强烈耦合，因此重子扰动无法增长——光子的压强把它们"推开"了。
但暗物质不同——它不与光子发生电磁相互作用，不受光子压强的阻碍。因此，暗物质的密度扰动可以在复合之前就开始增长，提前形成引力势阱。
复合之后，光子脱耦，重子物质失去压强支撑，迅速落入暗物质已经形成的引力势阱中。这就是为什么冷暗物质模型（CDM）能如此成功地解释结构形成——暗物质提供了先行的引力骨架，重子物质随后在其上聚集。
线性增长理论
在扰动幅度较小的线性阶段，密度对比度 \delta = (\rho - \bar{\rho})/\bar{\rho} 的增长可以用线性增长因子 D(t) 来描述。
在物质主导时期，密度对比度线性增长：\delta(t) \propto D(t) \propto a(t)（即正比于尺度因子）。这意味着密度对比度在物质主导期间持续增大。
在线性阶段，不同k模式的增长速率相同，因此扰动的形状保持不变。但转移函数 T(k) 会调制不同尺度的增长幅度：小尺度（大k）的扰动增长受到 Silk 阻尼效应的抑制（光子在复合前通过随机游走抹平了小尺度结构），导致小尺度功率谱被压低。
转移函数的特征在物质功率谱上表现为一个峰值——对应的尺度就是复合时刻声波传播的特征距离（约150 Mpc），这就是BAO尺度的物理起源。
从线性到非线性坍缩
当密度对比度 \delta 接近1时，线性近似不再成立，扰动进入非线性阶段。此时必须 resort to 数值模拟（N体模拟）或近似解析模型。
球对称坍缩是最简单的非线性模型。考虑一个密度略高于平均值的球形区域，在引力作用下它最终会停止跟随宇宙膨胀而开始坍缩。坍缩过程中，球形区域内的物质密度不断增加，最终达到 virial 平衡——形成暗物质晕。
暗物质晕的质量函数（Press-Schechter形式）描述了不同质量晕的数量密度，是连接理论预测与观测的重要工具。
⚠️ 热暗物质问题：如果暗物质是"热"的（速度接近光速），它会抹平小尺度的密度扰动，导致小尺度结构无法形成。观测到的矮星系数量和星系内部密度分布要求暗物质是"冷"的——这就是选择冷暗物质（CDM）模型的主要原因。
💡 关键直觉：结构形成的核心故事是"暗物质先行，重子跟随"。暗物质在宇宙不透明时期就悄悄搭好了引力脚手架，等到宇宙变透明、重子物质可以自由运动时，只需"走"进已经建好的框架中。没有这个脚手架，138亿年远远不够让宇宙从均匀状态形成今天的结构。
一个数值例子：扰动如何长大
把理论变成直觉的最好方式是算一笔具体的账。CMB显示复合时刻（红移 z \approx 1089）密度扰动的典型幅度约为 \delta \approx 10^{-5}。在物质主导时期，线性增长因子正比于尺度因子 D(t) \propto a(t)，而 1+z = 1/a，因此扰动幅度随红移的关系是
把复合时刻的数值代进去：\delta_{1089} \approx 10^{-5}，意味着今天（z=0）的线性外推幅度约为 \delta_0 \approx 10^{-5} \times 1089 \approx 10^{-2}。也就是说，线性理论预言今天的密度扰动约1%——但观测到的星系和星系团密度对比度远大于1，达到10到100的量级。这告诉我们两个事实：
- 大尺度（> 100 Mpc）上的扰动确实仍接近线性，约百分之几，可以用线性理论描述；
- 小尺度上的扰动早已进入非线性阶段，必须借助N体模拟或球对称坍缩模型来研究。
球对称坍缩模型给出了一个经典结论：一个密度略高的球形区域，在线性外推值达到 \delta_c \approx 1.686 时会完成"塌缩"（形成束缚的暗物质晕）。而实际达到维里平衡时，该区域相对背景的平均过密度约为200倍——这就是为什么天文学家常用"半径内平均密度为宇宙平均密度200倍的区域"来定义暗物质晕和星系团的边界（R_{200}）。
这个数值框架还有一个实际用途：它让我们可以把观测到的结构数量与暗物质性质直接挂钩。如果暗物质是"热的"（比如质量很轻的中微子），它会通过自由流抹平小尺度扰动，导致矮星系数量远少于观测；只有"冷"的暗物质才能在小尺度保留足够多的坍缩种子。因此，矮星系计数本身就是一个对暗物质温度（质量）的观测约束。
一节小结
- 原初功率谱：暴胀预言几乎标度不变的扰动，n_s \approx 0.965
- 金斯不稳定性：尺度大于金斯长度时，引力克服压强，扰动增长
- 暗物质先行：不受光子压强阻碍，提前形成引力势阱
- 重子跟随：复合后落入暗物质势阱，结构形成加速
- 线性增长：物质主导期 \delta \propto a(t)，转移函数调制不同尺度
- 非线性坍缩：\delta \approx 1 时进入非线性，数值模拟必需