大尺度结构形成

概述#

今天的宇宙在大尺度上（$\gtrsim$ 几百 Mpc）呈现出令人惊叹的泡沫状网格结构——宇宙网（cosmic web）：巨大的空洞（voids）被由星系构成的纤维（filaments）和墙（sheets）分隔，纤维的交叉节点处是密度最高的星系团（galaxy clusters）和超星系团。

这一宏大结构并非偶然。它是宇宙早期——宇宙微波背景（CMB）时代——极微弱的量子密度起伏（振幅约 $10^{-5}$）在引力作用下，历经约 138 亿年增长放大的结果。理解大尺度结构的形成，就是理解从量子世界到宇宙网的这段史诗旅程。

破除误解：结构不是"从无到有长出来"，是对比度被放大#

"宇宙从均匀变成有结构"这句话容易读成：原本什么都没有，后来凭空长出星系。准确的说法是：起伏一开始就在，只是极其微弱，引力把它的对比度放大了几万倍。

给个量：CMB 显示宇宙 38 万岁时的密度起伏幅度只有约 $10^{-5}$——十万分之一。今天一个星系团内部的密度是宇宙平均值的几百倍以上，一个巨洞里则不到十分之一。从 $10^{-5}$ 到 $10^2$，跨了七个数量级，全靠引力不稳定性。

机制本身很朴素，是一个正反馈：某处密度略高 → 引力略强 → 吸引更多物质 → 密度更高。稍稀疏处则相反，物质持续流失。所以引力不创造起伏，它只加剧已有的不平等。

由此有两个重要推论：

第一，结构形成的速度是宇宙学参数的探针。放大得多快，取决于物质密度有多高（引力多强）与膨胀有多快（引力多难赢）。所以"今天结构长到什么程度"直接约束 $\Omega_m$ 与暗能量——这就是 $S_8$ 参数的含义。

第二，必须有暗物质。普通物质在复合之前与光子紧密耦合，辐射压会把起伏抹平，无法提前坍缩。暗物质不与光子作用，因此能在复合之前就开始成团，提前挖好引力势阱，等普通物质一解放就落进去。如果只有普通物质，138 亿年不够长出今天的结构。 这是暗物质最有力的证据之一，且完全独立于星系旋转曲线。

现场：模拟成了检验理论的主要手段#

引力不稳定性一旦进入非线性阶段（密度对比超过 1），方程就没有解析解了。于是这门学科出现了一个不寻常的局面：理论预言主要靠数值模拟给出，而不是靠公式。流程是这样的：把 CMB 测到的初始功率谱作为初始条件，在盒子里撒下几十亿到上千亿个粒子代表暗物质，让它们只在引力下相互作用，按宇宙学模型演化到今天，然后把结果"观测"一遍——统计成团程度、数星系团、量丝状结构——与真实巡天对比。

这条路线给出了 ΛCDM 最漂亮的成功：模拟自然产生了红移巡天里看到的丝状网络与巨洞，包括统计意义上的尺度分布。但它也暴露了最顽固的问题，都集中在小尺度上：

• 卫星星系数目：模拟预言银河系周围应有大量小暗物质晕，观测到的矮星系却少得多（"missing satellites"）。
• 核心-尖点问题：模拟给出的暗物质晕中心密度是陡峭的尖点，许多矮星系的观测却指向平坦的核心。

这些张力的解释路线泾渭分明：一派认为问题在重子物理没算好（超新星反馈会把矮星系里的气体吹走、也会把暗物质晕中心"加热"变平），另一派认为问题在暗物质本身（可能是温暗物质或自相互作用暗物质）。

值得注意的是这个格局：大尺度上 ΛCDM 精确到几个百分点，小尺度上争议不断。 而恰恰是小尺度最难算——因为那里重子物理最复杂。所以这场争论的本质，是"我们不确定这是新物理，还是我们的气体动力学不够好"。

初始条件：来自暴胀的量子起伏#

宇宙暴胀（Cosmic Inflation）模型预言，暴胀场（inflaton field）的量子涨落被暴胀放大，产生近似高斯分布、接近尺度无关的原初密度起伏谱：

$$P(k) \propto k^{n_s}$$

其中 $n_s$ 是标量谱指数（scalar spectral index）。Planck 2018 测量 $n_s = 0.9649 \pm 0.0042$（Planck Collaboration 2020，A&A 641，A10），轻微偏离尺度无关（$n_s = 1$），与简单暴胀模型预言一致。CMB 的温度各向异性（约 $\delta T/T \sim 10^{-5}$）直接映射了宇宙复合时代（$z \approx 1100$）的密度起伏——这就是大尺度结构的"初始条件"快照。

线性增长：引力不稳定性#

物质密度涨落的增长#

在物质主导时代，对于尺度远大于 Jeans 长度的密度起伏 $\delta(\mathbf{x}, t) = (\rho - \bar\rho)/\bar\rho$，在线性近似下，密度对比的增长方程为：

$$\ddot{\delta} + 2H\dot{\delta} - \frac{3}{2} H_0^2 \Omega_m (1+z)^3 \delta = 0$$

其解为 $\delta(t) = D_+(t) \delta_0$，其中 $D_+(t)$ 是线性增长因子（linear growth factor）。在爱因斯坦-德西特宇宙（$\Omega_m = 1$，$\Omega_\Lambda = 0$）中，$D_+ \propto a \propto t^{2/3}$（即密度起伏随膨胀因子线性增长）。在 $\Lambda$CDM 宇宙中，暗能量主导之后 $D_+$ 增长减缓。

功率谱演化#

物质功率谱 $P(k, z) = D_+^2(z) T^2(k) P_0(k)$ 描述了不同空间尺度（波数 $k$）上密度涨落的振幅。$T(k)$ 是转移函数（transfer function），描述在辐射-物质相等时代之前，各波数模式的进化修正。CDM 的转移函数在小尺度（$k \gg k_\mathrm{eq}$）有 $\sim k^{-2}$ 的压缩，反映了在辐射主导时代进入视界的模式被梅塞特效应（Meszaros effect）压制。

非线性结构形成：从壳层到网#

当 $\delta \sim 1$（密度涨落与背景相当）时，线性理论失效，结构进入非线性坍缩阶段。

Zel'dovich 近似#

Zel'dovich（1970）给出了从线性到弱非线性的半解析近似，描述物质在引力下从均匀分布向"煎饼"结构（pancakes）坍缩的过程：

$$\mathbf{r}(\mathbf{q}, t) = \mathbf{q} - D_+(t) \nabla \Phi_0(\mathbf{q})$$

其中 $\mathbf{q}$ 是初始（拉格朗日）坐标，$\Phi_0$ 是初始引力势。Zel'dovich 近似预言，坍缩沿密度场的最大曲率方向优先发生，形成薄的二维"壳层"——这正是宇宙网纤维和墙状结构的理论起源。

晕的形成：Press-Schechter 理论#

Press & Schechter（1974，ApJ 187，425）给出了第一个解析预测不同质量的暗物质晕数量密度的框架。假设密度场为高斯随机场，对高斯随机场进行平滑（对应不同质量尺度 $M$），密度起伏超过临界值（$\delta_c \approx 1.686$，由球对称坍缩计算）的区域坍缩为束缚暗物质晕：

$$\frac{dn}{d\ln M} = \sqrt{\frac{2}{\pi}} \frac{\bar\rho}{M} \frac{\delta_c}{\sigma(M)} \left|\frac{d\ln\sigma}{d\ln M}\right| \exp\left(-\frac{\delta_c^2}{2\sigma^2(M)}\right)$$

这一晕质量函数（halo mass function，HMF）的定性形式与 N 体模拟高度吻合，现代精确版本（如 Sheth-Tormen 或 Tinker 函数）将误差控制在约 5%–10%（Tinker et al. 2008，ApJ 688，709）。HMF 对 $\Omega_m$ 和 $\sigma_8$ 极为敏感，是利用星系团数量约束宇宙学参数的理论基础。

N 体模拟：数值宇宙#

现代大规模 N 体（或流体力学）数值模拟是研究非线性大尺度结构形成的主要工具：

Springel 等（2005，Nature 435，629）发表的 Millennium 模拟再现了从大尺度纤维网到单个星系的多尺度结构，是迄今引用最多的宇宙学模拟之一。流体力学模拟（如 IllustrisTNG、EAGLE、SIMBA）在重力模拟的基础上加入气体、恒星形成、超新星反馈和 AGN 反馈，从第一性原理重现了星系的颜色双模分布、质量函数和金属丰度，是现代星系形成理论的核心工具。

宇宙网的观测#

星系分布的丝状结构#

二维投影下可直接看到宇宙网的壮观图像：CfA2 "长城"（Geller & Huchra 1989，Science 246，897）是最早被识别的大尺度结构之一，长约 5 亿光年（约 150 Mpc）、宽约 2 亿光年（约 60 Mpc），当时是已知最大的宇宙结构。SDSS 的三维星系分布图（$z < 0.2$，约10万个星系）清晰展示了纤维、空洞和超星系团的层次结构（Tegmark et al. 2004，ApJ 606，702）。

空洞统计#

宇宙空洞（voids）占宇宙体积约 60%–80%，是宇宙网中最大的单一结构。利用 SDSS 和 2dFGRS 数据识别的空洞直径约 10–100 Mpc，平均密度约为宇宙平均的 20%–30%。空洞的形态和大小分布对暗能量态方程 $w$ 和中微子质量敏感（Pisani et al. 2019，BAAS 51，40）。

大尺度结构与宇宙学参数#

大尺度结构是约束宇宙学参数的多功能工具：

• 物质功率谱形状：约束 $\Omega_m h$（物质密度与哈勃常数之积）
• BAO 峰位置：约束暗能量（详见 重子声学振荡）
• 增长率 $f\sigma_8$（通过红移空间畸变）：检验广义相对论与修改引力理论
• 星系团数量：约束 $\sigma_8$ 和 $\Omega_m$
• 弱引力透镜（详见 弱引力透镜与宇宙剪切）：约束 $S_8$

中微子质量与大尺度结构#

中微子是宇宙中已知质量最小的粒子之一，但其总质量（三代之和 $\sum m_\nu$）对宇宙大尺度结构有可测量的影响：中微子在宇宙早期以相对论速度自由流动（free streaming），抹平了自由流长度 $\lambda_\mathrm{fs}$ 以下的密度起伏。三代中微子的总质量约束：

$$\lambda_\mathrm{fs} \approx 340\ \text{Mpc} \cdot \left(\frac{m_\nu}{1\ \text{eV}}\right)^{-1}$$

在 $k > k_\mathrm{fs}$ 尺度，物质功率谱的抑制量约为 $\Delta P/P \approx -8 f_\nu$，其中 $f_\nu = \Omega_\nu/\Omega_m$ 是中微子对物质密度的贡献份额。Planck CMB + BAO 联合约束 $\sum m_\nu < 0.12$ eV（95% 置信，Planck 2018），而 DESI 2024 进一步给出 $\sum m_\nu < 0.072$ eV（结合CMB+BAO），接近粒子物理中微子振荡实验给出的下限（约 0.06 eV）。大尺度结构（尤其是 BAO + WL + 星系功率谱）未来有望首次测量而非仅约束中微子绝对质量。

跨域连接#

• 宇宙网与纤维状结构：坍缩沿密度场曲率最大的方向优先发生，因此结构先压成面、再收成丝、最后在交点成团。这条顺序是纤维状形态的理论来源，也给出一个可检验的推论：不同形态成分各占可预测的体积与质量份额，而不是随机分布。
• 动力系统：密度对比小于一时线性理论足够，超过一之后方程失去解析解，弱非线性还能用摄动展开，强非线性只剩数值。这条分界不是技术上的懒惰，而是结构性的：同一个正反馈在小振幅下可线性化，在大振幅下必然产生尺度之间的耦合。
• 相变与临界现象：近似尺度无关的初始谱加上引力的自相似，给出层级式成团与幂律形式的质量函数；坍缩阈值扮演临界值的角色，超过阈值的稀有涨落被指数压制。因此最大质量结构的数目对涨落振幅极其敏感，这正是用团数约束宇宙学参数的依据。
• 分治法：直接两两求引力是粒子数平方的代价，千亿粒子根本算不动。把远处质量按树结构或网格合并成等效项，复杂度降到近似线性对数。推论：模拟的分辨率上限由算法而不是机器决定，而分辨率不足会人为抑制小尺度结构。
• 气候模型：两边共有同一个困境：小尺度过程必须参数化，而结论对参数化敏感。于是"这是新物理还是子网格物理不够好"成了同一种争论——矮星系数目与暗物质晕内部密度轮廓的分歧，与云反馈之争在方法论上完全同构。

参考文献#

• Press, W. H. & Schechter, P. (1974). Formation of galaxies and clusters of galaxies by self-similar gravitational condensation. Astrophysical Journal, 187, 425–438.
• Springel, V. et al. (2005). Simulations of the formation, evolution and clustering of galaxies and quasars. Nature, 435, 629–636.
• Geller, M. J. & Huchra, J. P. (1989). Mapping the universe. Science, 246, 897–903.
• Zeldovich, Ya. B. (1970). Gravitational instability: an approximate theory for large density perturbations. Astronomy & Astrophysics, 5, 84–89.
• Planck Collaboration (2020). Planck 2018 results. X. Constraints on inflation. Astronomy & Astrophysics, 641, A10.
• Pillepich, A. et al. (2018). Simulating galaxy formation with the IllustrisTNG model. Monthly Notices of the Royal Astronomical Society, 473, 4077–4106.