现代宇宙学标准模型：ΛCDM框架的确立与观测基石

现代宇宙学标准模型——Λ冷暗物质（Lambda Cold Dark Matter, ΛCDM）模型——是描述宇宙从暴胀后至今演化的最成功理论框架。该模型以广义相对论为基础，结合宇宙学原理（假设宇宙在大尺度上均匀且各向同性），将宇宙组分归结为普通重子物质、冷暗物质以及宇宙学常数Λ（代表暗能量）三大部分。其核心动力学由弗里德曼方程描述，而模型的可靠性则建立在多重独立观测的坚实基石之上。本章将严格推导ΛCDM模型的数学基础，系统梳理支撑该模型的六大核心观测支柱，深入分析这些观测如何共同约束宇宙的关键参数，并探讨当前模型内部存在的显著张力，为后续深入探讨暗物质与暗能量的物理机制设定精确的起点。

1.1 弗里德曼方程：ΛCDM模型的数学基石

ΛCDM模型的理论起点是爱因斯坦场方程： $$R_{\mu\nu} - \frac{1}{2}R g_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$$ 在宇宙学原理的约束下，度规取弗里德曼-勒梅特-罗伯逊-沃尔克（FLRW）形式： $$ds^2 = -c^2 dt^2 + a(t)^2 \left[ \frac{dr^2}{1-kr^2} + r^2(d\theta^2 + \sin^2\theta d\phi^2) \right]$$ 其中 $a(t)$ 为尺度因子，$k$ 为曲率常数（$k=0, \pm1$）。假设宇宙内容为理想流体，能量动量张量为 $T^{\mu}_{\ \nu} = \text{diag}(-\rho c^2, p, p, p)$，其中 $\rho$ 为能量密度，$p$ 为压强。将FLRW度规代入场方程，可推导出两个独立方程：

第一弗里德曼方程（动力学方程）： $$H^2(t) \equiv \left(\frac{\dot{a}}{a}\right)^2 = \frac{8\pi G}{3}\rho - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3}$$ 第二弗里德曼方程（加速度方程）： $$\frac{\ddot{a}}{a} = -\frac{4\pi G}{3}\left(\rho + \frac{3p}{c^2}\right) + \frac{\Lambda c^2}{3}$$

定义临界密度 $\rho_c(t) = \frac{3H^2(t)}{8\pi G}$，则各组分可引入密度参数 $\Omega_i(t) = \rho_i(t)/\rho_c(t)$。对于平坦宇宙（$k=0$），第一方程简化为： $$\sum_i \Omega_i(t) = 1$$ 其中求和包含所有组分。在ΛCDM模型中，主要组分为：

• **重子物质**（Baryons, b）：$p \approx 0$，$\Omega_b$
• **冷暗物质**（Cold Dark Matter, c）：$p \approx 0$，$\Omega_c$
• **辐射**（Photons + 中微子, r）：$p = \rho c^2/3$，$\Omega_r$（ today 极小，但早期主导）
• **宇宙学常数**（Λ）：$p = -\rho c^2$，$\Omega_\Lambda$

因此，当代宇宙的总密度参数为： $$\Omega_{\text{tot}} = \Omega_b + \Omega_c + \Omega_\Lambda + \Omega_r \approx 1$$ 观测强烈支持 $\Omega_{\text{tot}}$ 非常接近1，意味着宇宙空间几何近乎平坦（$k \approx 0$）。这一平坦性可通过CMB各向异性功率谱的第一声学峰位置精确测定。

关键参数定义：

• $H_0$：当前哈勃常数，单位 km/s/Mpc。常用值约67-74。
• $h$：无量纲哈勃参数，$H_0 = 100h$ km/s/Mpc。
• $\Omega_b h^2$：重子物理密度，CMB直接测量。
• $\Omega_c h^2$：冷暗物质物理密度，CMB直接测量。
• $\Omega_\Lambda$：暗能量密度参数，由 $\Omega_{\text{tot}}$ 和物质密度推导。
• $n_s$：标量谱指数，描述原初扰动谱的倾斜。
• $\tau$：光学深度，表征宇宙再电离程度。

实践案例：从Planck数据反推参数 Planck卫星2018年释放的最终数据给出了基于ΛCDM模型的精确参数估计。例如，使用TT,TE,EE+lowE+lensing数据组合，得到： $$\Omega_b h^2 = 0.02237 \pm 0.00015,\quad \Omega_c h^2 = 0.1200 \pm 0.0012,\quad H_0 = 67.4 \pm 0.5\ \text{km/s/Mpc},\quad n_s = 0.9649 \pm 0.0042,\quad \tau = 0.0544 \pm 0.0073$$ 这些数值并非直接测量，而是通过马尔可夫链蒙特卡洛（MCMC）方法对理论功率谱与观测数据进行拟合后得到的后验分布统计量。其不确定性包含统计误差和系统误差的联合贡献。

1.2 六大核心观测支柱：从独立证据到联合约束

ΛCDM模型的成立并非依赖单一观测，而是由多种独立、互补的观测手段共同验证。每一根“支柱”从不同角度测量宇宙的几何、组分和演化历史，它们的一致性构成了现代宇宙学的黄金标准。

1.2.1 宇宙微波背景辐射（CMB）各向异性：早期宇宙的快照

CMB是宇宙再复合时期（约38万年）遗留下来的热辐射，其温度各向异性（约$\Delta T/T \sim 10^{-5}$）携带了早期宇宙物理条件的丰富信息。WMAP和Planck卫星已以前所未有的精度绘制了全天空CMB温度与偏振图。

物理机制：原初量子涨落在暴胀期间被放大为密度扰动。在光子-重子流体中，这些扰动激发声波，导致压缩与稀疏交替，形成特征性的声学峰。峰的位置、高度和宽度对宇宙参数高度敏感：

• **第一峰位置**：主要决定空间曲率。对于平坦宇宙，多极矩 $\ell_1 \approx 220$（对应角度约1°），与观测完美吻合。
• **峰高比**：第一峰与第二峰高度比敏感于重子密度 $\Omega_b h^2$（重子增加流体惯性，抑制峰值）。
• **阻尼尾部**：小尺度（高$\ell$）衰减受光子扩散长度影响，约束 $\Omega_b h^2$ 和 $n_s$。
• **偏振谱（EE, TE）**：提供独立于温度的约束，尤其对再电离光学深度 $\tau$ 极为敏感。

Planck数据深度分析：Planck 2018释放的TT功率谱覆盖 $2 \leq \ell \leq 2500$，TE和EE谱至 $\ell \approx 2000$。图1.1（建议电子版插入动态图表，展示Planck TT+TE+EE最佳拟合曲线与数据点）显示，理论曲线与观测数据在绝大多数范围内吻合极佳，残差随机分布，表明ΛCDM模型对CMB的描述极为成功。通过MCMC拟合，CMB独立给出 $\Omega_b h^2$、$\Omega_c h^2$、$n_s$、$\tau$ 的高精度估计，并间接推断 $\Omega_\Lambda$（需结合平坦性假设）。

技术细节：CMB数据分析涉及复杂的束函数校正、前景分离（银河系同步辐射、尘埃）、仪器噪声建模以及协方差矩阵估计。Planck合作组使用`CamSpec`和`Plik`等 likelihood 代码，在参数空间进行大规模MCMC采样（如`CosmoMC`或`MontePython`），产生后验分布。图1.2（建议电子版插入参数约束三角形图）展示主要参数的两两联合分布，清晰显示参数间的强相关性，如 $\Omega_c h^2$ 与 $n_s$ 的负相关。

1.2.2 Ia型超新星：宇宙加速膨胀的发现者

Ia型超新星源于白矮星的热核爆炸，其峰值光度高度一致（经光变曲线形状校准后），可作为“标准烛光”测量共动距离。1998年，两个团队通过观测高红移（$z \sim 0.5-1$）Ia型超新星发现，其距离比预期更远，意味着宇宙膨胀在约50亿年前开始加速，直接揭示了暗能量的存在。

距离模数与红移关系：观测到的视亮度 $m$ 与绝对亮度 $M$ 之差为距离模数 $\mu = m - M = 5 \log_{10}(d_L/\text{Mpc}) + 25$，其中 $d_L$ 为光度距离。在FLRW宇宙中： $$d_L(z) = (1+z) \frac{c}{H_0} \frac{1}{\sqrt{|\Omega_k|}} \sinh\left( \sqrt{|\Omega_k|} \int_0^z \frac{dz'}{E(z')} \right)$$ 其中 $E(z) = H(z)/H_0$，对于ΛCDM： $$E(z) = \sqrt{\Omega_m (1+z)^3 + \Omega_\Lambda + \Omega_r (1+z)^4 + \Omega_k (1+z)^2}$$ $\Omega_m = \Omega_b + \Omega_c$。通过测量 $\mu(z)$，可约束 $d_L(z)$，进而限制 $\Omega_\Lambda$ 和 $\Omega_m$。

实践案例：Pantheon+数据集：当前最全面的Ia型超新星汇编为Pantheon+，包含约1700个超新星，红移范围 $0.01 < z < 2.3$。图1.3（电子版可插入 $\mu$ vs $z$ 对比图）显示，数据点与不同宇宙学模型的预测曲线。在平坦ΛCDM下，联合CMB和BAO可给出 $\Omega_\Lambda \approx 0.69$，与超新星数据完美契合。然而，系统误差是主要挑战：尘埃消光、星族演化、校准链（从局部超新星到Cepheid变星再到几何距离）均需仔细处理。SH0ES团队通过HST观测Cepheid校准近邻超新星，得到 $H_0 = 73.04 \pm 1.04$ km/s/Mpc，与CMB推断值形成显著张力（见1.4节）。

1.2.3 大尺度结构（LSS）成团性：物质分布的指纹

星系巡天（如SDSS、2dF、BOSS、eBOSS、DESI）绘制了数百万星系的三维分布图，揭示宇宙在百Mpc尺度上呈现“宇宙网”结构：星系团、纤维状结构和空洞。这种成团性源于原初密度扰动的引力增长。

统计工具：两点相关函数与功率谱：

• **相关函数** $\xi(r)$：描述超过随机分布的概率，$\xi(r) = \langle \delta(\mathbf{x})\delta(\mathbf{x+r}) \rangle$，其中 $\delta = (\rho - \bar{\rho})/\bar{\rho}$。
• **功率谱** $P(k)$：$\xi(r)$ 的傅里叶变换，$P(k) = \langle |\delta_k|^2 \rangle$，描述不同尺度 $k$ 的扰动幅度。

在ΛCDM中，原初标量扰动功率谱为 $P_{\text{prim}}(k) = A_s (k/k_0)^{n_s-1}$，经转移函数 $T(k)$ 修正（反映辐射-物质相等前后的演化），得到 today 的物质功率谱： $$P(k,z=0) = \frac{2\pi^2}{k^3} A_s \left(\frac{k}{k_0}\right)^{n_s-1} T^2(k) D^2(z=0)$$ 其中 $D(z)$ 为增长因子，满足微分方程： $$\frac{d^2 D}{d \ln a^2} + \left[2 + \frac{1}{2}\frac{d\ln H^2}{d \ln a}\right] \frac{d D}{d \ln a} = \frac{3}{2} \Omega_m(a) D$$ 增长因子对 $\Omega_m$ 和 $\Omega_\Lambda$ 敏感，因为暗能量抑制结构增长。

观测约束：SDSS主星系巡天（MGS）和BOSS的CMASS样本测量了 $P(k)$ 和 $\xi(r)$，显示特征性的重子声学振荡（BAO）峰（见下节）。同时，功率谱的形状和幅度约束 $\Omega_m h$ 和 $\sigma_8$（8 Mpc/h球内密度涨落RMS）。LSS数据与CMB结合，能检验增长历史是否与广义相对论预言一致，目前未发现显著偏离。

案例：BOSS DR12分析：通过测量CMASS和LOWZ星系的BAO信号，结合红移空间扭曲（RSD）测量 $f\sigma_8$（增长速率），联合Planck数据，得到 $\Omega_m = 0.310 \pm 0.006$，$\sigma_8 = 0.800 \pm 0.014$。这些结果与CMB外推至 today 的预言一致，但轻微存在“$S_8$ 张力”（见1.4节）。

1.2.4 重子声波振荡（BAO）：宇宙的标准尺

在光子-重子流体中，原初扰动引发的声波在光子退耦前传播约一半波长，在物质功率谱中留下特征尺度——声视界 $r_d$，即声波在再复合前能传播的最大距离： $$r_d = \int_{z_d}^{\infty} \frac{c_s(z)}{H(z)} dz \approx 147\ \text{Mpc}$$ 其中 $c_s = c/\sqrt{3(1+R)}$ 为声速，$R = 3\rho_b/(4\rho_\gamma)$，$z_d \approx 1059$ 为退耦红移。这一尺度在CMB温度功率谱中表现为第一峰位置，在LSS中表现为相关函数在约150 Mpc处的隆起。

测量方法：通过星系巡天测量两点相关函数 $\xi(r)$ 或功率谱 $P(k)$，在径向上（沿视线方向）和横向上（垂直于视线）分别提取BAO尺度 $D_V(z)/r_d$、$D_M(z)/r_d$ 和 $H(z) r_d$，其中 $D_M$ 为角直径距离，$D_V = [z D_M^2 c / H(z)]^{1/3}$ 为复合距离。这些测量提供绝对距离标尺，无需校准。

关键巡天与结果：

• **SDSS-III BOSS**：在 $z \approx 0.38, 0.51, 0.61$ 测量BAO，精度达1%量级。
• **6dFGS**：低红移（$z \approx 0.106$）BAO测量。
• **eBOSS**：扩展至类星体（$z \approx 1.5$）和Lyα森林（$z \approx 2.3$）。
• **DESI**（进行中）：目标测量数千万星系和类星体的BAO，精度将达0.5%以下。

宇宙学约束：BAO测量直接给出 $D_M(z)/r_d$ 和 $H(z) r_d$。若 $r_d$ 已知（由CMB精确测定），即可约束 $H_0$ 和 $\Omega_\Lambda$。联合CMB和BAO可打破参数简并，例如显著改善 $\Omega_m$ 和 $H_0$ 的约束。当前BAO数据与CMB结合，在平坦ΛCDM下给出 $H_0 = 67.4 \pm 0.5$ km/s/Mpc，与CMB自洽，但与局部测量存在分歧。

1.2.5 宇宙年龄测定：一致性检验

宇宙年龄 $t_0$ 可从哈勃参数和密度参数积分得到： $$t_0 = \int_0^\infty \frac{dz}{(1+z)H(z)}$$ 在ΛCDM中，给定 $H_0$、$\Omega_m$、$\Omega_\Lambda$，可计算 $t_0$。同时，独立测定最古老天体的年龄提供下限。两者一致性是模型的重要检验。

最老星体年龄：球状星团中贫金属恒星通过核合成年代学（如Th/U比值）或恒星演化模型测定年龄。目前最老星团年龄估计为 $13.4 \pm 0.8$ Gyr。白矮星冷却序列也给出类似下限。

CMB推断年龄：Planck数据在ΛCDM下给出 $t_0 = 13.797 \pm 0.023$ Gyr，与星团年龄在1σ内一致。这一一致性支持模型从早期宇宙到 today 的演化图景。

1.2.6 大爆炸核合成（BBN）：轻元素丰度的预言

在宇宙最初几分钟（$T \sim 0.1-1$ MeV），质子和中子通过弱相互作用保持化学平衡，随后中子衰变和核合成形成轻元素：D, $^3$He, $^4$He, $^7$Li。BBN理论高度依赖重子-光子比 $\eta = n_b/n_\gamma$，该参数与 $\Omega_b h^2$ 直接相关： $$\Omega_b h^2 = \eta \frac{m_p c^2 \rho_c}{\rho_\gamma} \approx 3.65 \times 10^{-3} \eta$$ 其中 $\rho_\gamma$ 为 today 光子密度。

预言与观测：

• **氘（D）**：丰度对 $\eta$ 极为敏感（$D/H \propto \eta^{-1.6}$）。高红移类星体吸收系统（DLA）测量 $D/H = (2.527 \pm 0.030) \times 10^{-5}$，对应 $\Omega_b h^2 = 0.02233 \pm 0.00015$，与CMB独立结果惊人一致。
• **氦-4（$^4$He）**：质量分数 $Y_p$ 对 $\eta$ 敏感度较低，但受中子寿命影响。观测金属贫乏星系 $Y_p = 0.2449 \pm 0.0040$，与预言吻合。
• **锂-7（$^7$Li）**：存在显著矛盾：预言丰度比观测高约3倍（“锂问题”），可能涉及恒星物理或新物理。

BBN成功预言D和$^4$He丰度，为 $\Omega_b$ 提供了独立于CMB的验证，且覆盖宇宙演化两个极端时期（分钟 vs. 38万年），是ΛCDM模型的关键支柱。

1.3 参数约束：协方差、MCMC与联合分析

单一观测无法唯一确定所有参数，需联合分析多个数据集以打破简并。例如，CMB精确测定 $\Omega_b h^2$、$\Omega_c h^2$、$n_s$、$\tau$，但对 $\Omega_\Lambda$ 和 $H_0$ 的约束依赖于平坦性假设；BAO提供绝对距离，直接约束 $\Omega_m$ 和 $H_0$；超新星提供 $d_L(z)$，约束 $\Omega_\Lambda$。联合后，参数不确定性大幅缩小。

1.3.1 似然函数与协方差矩阵

给定理论模型参数 $\theta = (\Omega_b h^2, \Omega_c h^2, H_0, \tau, n_s, \dots)$，观测数据 $d$ 的似然函数为： $$\mathcal{L}(\theta) \propto \exp\left