宇宙学标准模型：框架、参数与挑战

现代宇宙学在二十世纪末叶完成了一次深刻的范式转移，从基于推测的理论思辨，演变为一门以精确观测为基石、以数学模型为框架的精密科学。这一转变的核心成果，便是宇宙学标准模型——即Λ冷暗物质模型（Lambda Cold Dark Matter, ΛCDM）的确立。该模型以广义相对论为理论基础，以弗里德曼-勒梅特-罗伯逊-沃尔克（FLRW）度规描述均匀各向同性的宇宙几何，并引入两种神秘的成分——冷暗物质（Cold Dark Matter, CDM）和宇宙学常数Λ（代表暗能量）——来解释从宇宙微波背景辐射（CMB）的各向异性到星系大尺度结构的形成等一系列观测事实。其简洁而强大的解释力，使其成为描述宇宙从暴胀后演化至今的“标准模型”。然而，这个模型的成功本身就孕育着深刻的物理学危机：宇宙常数问题所揭示的真空能量灾难、哈勃常数（$H_0$）的局部与全局测量之间的显著张力、以及暗物质与暗能量本质的完全未知，共同构成了对现有物理图景的严峻挑战。本章旨在系统回顾ΛCDM模型的数学与物理基础，梳理其关键参数的定义与观测确定方法，并批判性地审视其在解释观测数据方面的成功与内在的理论困境。这不仅是对现代宇宙学基石的一次严谨检阅，更是为后续深入探讨暗物质与暗能量之本质所必需的理论铺垫。

宇宙学标准模型的数学框架建立在爱因斯坦场方程对均匀各向同性时空的应用之上。广义相对论将时空的几何与物质的能量-动量分布联系起来，其核心方程为： $$G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$$ 其中$G_{\mu\nu}$是爱因斯坦张量，$\Lambda$是宇宙学常数，$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$是归一化曲率常数，取值为+1（闭合）、0（平坦）或-1（开放）。将FLRW度规代入爱因斯坦场方程，并假设宇宙内容为理想流体（$T^{\mu}_{\nu} = \text{diag}(-\rho c^2, p, p, p)$），即可推导出描述宇宙动力学演化的两个弗里德曼方程：

第一个弗里德曼方程（膨胀方程）： $$H^2(t) \equiv \left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho(t) - \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}$$

这两个方程是ΛCDM模型的动力学核心。它们表明，宇宙的膨胀速率$H$由三部分贡献：物质-辐射的能量密度$\rho$（导致减速）、空间曲率$k$（正曲率减速，负曲率加速）以及宇宙学常数$\Lambda$（导致加速膨胀）。方程还隐含了能量守恒定律：$\dot{\rho} + 3H(\rho + p/c^2) = 0$。对于不同成分的宇宙组分，其状态方程$p = w\rho c^2$不同：辐射$w=1/3$，物质（包括重子和冷暗物质）$w=0$，宇宙学常数$w=-1$。由此，各组分的密度随标度因子演化规律为$\rho \propto a^{-3(1+w)}$。因此，物质密度$\rho_m \propto a^{-3}$，辐射密度$\rho_r \propto a^{-4}$，而暗能量密度$\rho_\Lambda = \text{常数}$。这意味着在宇宙早期，辐射主导；随后物质主导；而在最近约50亿年，暗能量主导，导致宇宙加速膨胀。

在FLRW框架下，定义关键的宇宙学密度参数是分析宇宙几何与命运的核心。对于任意组分$i$，其密度参数定义为该组分当前密度与临界密度的比值： $$\Omega_i = \frac{\rho_{i,0}}{\rho_{c,0}}, \quad \text{其中} \quad \rho_{c,0} = \frac{3H_0^2}{8\pi G}$$ 临界密度$\rho_{c,0}$是将宇宙几何弯曲为平坦（$k=0$）所需的恰好密度。当前总密度参数$\Omega_{\text{tot}} = \sum_i \Omega_i$与曲率$k$的关系由第一个弗里德曼方程在$t=t_0$时给出： $$\Omega_{\text{tot}} - 1 = \frac{kc^2}{a_0^2 H_0^2}$$ 因此，$\Omega_{\text{tot}} = 1$对应平坦宇宙（$k=0$）；$\Omega_{\text{tot}} > 1$对应闭合宇宙（$k=+1$）；$\Omega_{\text{tot}} < 1$对应开放宇宙（$k=-1$）。ΛCDM模型假设宇宙由重子物质（baryons）、冷暗物质（CDM）、辐射（包括光子和中微子）以及暗能量（以$\Lambda$形式）组成。因此，总密度参数可分解为： $$\Omega_{\text{tot}} = \Omega_b + \Omega_{\text{cdm}} + \Omega_r + \Omega_\Lambda$$ 其中，$\Omega_b$是重子密度参数，$\Omega_{\text{cdm}}$是冷暗物质密度参数，$\Omega_r$是辐射密度参数（主要包括CMB光子和大爆炸核合成遗迹的中微子），$\Omega_\Lambda$是暗能量密度参数。在标准模型中，辐射密度在当今极为微小（$\Omega_r \sim 10^{-4}$），通常可忽略。模型的核心预言是宇宙在几何上高度接近平坦，即$\Omega_{\text{tot}} \approx 1$，这要求$\Omega_b + \Omega_{\text{cdm}} + \Omega_\Lambda \approx 1$。这一“平坦性”预言得到了CMB观测的强力支持，并成为模型的重要支柱。

宇宙学参数的精确测定是检验ΛCDM模型的关键。这些参数包括：

• **哈勃常数 $H_0$**：描述当前宇宙膨胀速率，单位为 km/s/Mpc。其倒数即为哈勃时间 $t_H = 1/H_0$，是宇宙年龄的粗略量级估计。
• **物质密度参数 $\Omega_m$**：通常指总物质（重子+冷暗物质）密度参数，$\Omega_m = \Omega_b + \Omega_{\text{cdm}}$。
• **重子密度参数 $\Omega_b h^2$**：其中$h = H_0 / (100 \text{ km/s/Mpc})$是归一化哈勃参数。使用$\Omega_b h^2$而非$\Omega_b$可消除$H_0$的强相关性，是CMB分析中的标准表述。
• **暗能量密度参数 $\Omega_\Lambda$**：在平坦宇宙假设下，$\Omega_\Lambda = 1 - \Omega_m$。
• **物质功率谱标度 $A_s$ 与谱指数 $n_s$**：描述原初密度扰动的幅度与尺度依赖性。
• **再复合红移 $z_*$**：CMB最后散射面的红移，约$z_* \approx 1100$。

这些参数的确定依赖于多种独立的宇宙学观测手段，它们从不同角度、不同物理过程对模型进行交叉验证。其中，宇宙微波背景辐射（CMB）各向异性的精密测量、重子声波振荡（BAO）作为“标准尺”的应用，以及局域距离阶梯法测量哈勃常数，构成了当前参数测定的三大支柱。

宇宙微波背景辐射是大爆炸的余晖，其温度各向异性（约百万分之一量级）携带了宇宙在再复合时期（约38万年）的初始条件与物理信息。CMB温度功率谱$C_\ell$（作为多极矩$\ell$的函数）呈现出清晰的声学峰结构，这是光子-重子流体在引力势阱中振荡的直接印记。峰的位置、高度和间距对宇宙学参数极为敏感：

• **第一峰位置**：主要取决于宇宙的几何曲率。在平坦宇宙中，第一峰位于$\ell \approx 220$。其精确位置直接约束$\Omega_{\text{tot}}$。Planck卫星的测量给出了$\ell_1 \approx 220.0 \pm 0.5$，强烈支持$\Omega_{\text{tot}} = 1.0007 \pm 0.0019$（Planck 2018结果），表明宇宙在1%精度内平坦。
• **峰高比**：第一峰与第二峰的高度比主要敏感于物质密度$\Omega_m h^2$与重子密度$\Omega_b h^2$的比值。更高的重子密度增强声波振荡的压缩相（第一峰），而物质密度影响势阱深度。精确测量峰高比可独立确定$\Omega_b h^2$和$\Omega_m h^2$。
• **阻尼尾部**：小角尺度（高$\ell$）功率的指数衰减（Silk阻尼）由光子扩散长度决定，对$\Omega_b h^2$和再复合物理敏感。

通过将理论计算的CMB功率谱（基于Boltzmann方程求解，如CAMB或CLASS代码）与观测数据进行马尔可夫链蒙特卡洛（MCMC）拟合，可以同时约束多个参数。欧洲空间局的Planck卫星任务（2013, 2015, 2018数据释放）提供了迄今最精确的全天CMB温偏振测量。其2018年结果给出的基准ΛCDM参数为（TT,TE,EE+lowE+ lensing）：

• $H_0 = 67.4 \pm 0.5$ km/s/Mpc
• $\Omega_b h^2 = 0.02237 \pm 0.00015$
• $\Omega_c h^2 = 0.1200 \pm 0.0012$（$\Omega_c$即$\Omega_{\text{cdm}}$）
• $\tau = 0.054 \pm 0.007$（光学深度）
• $n_s = 0.9649 \pm 0.0042$
• $\ln(10^{10} A_s) = 3.044 \pm 0.014$

这些参数以极高的精度（通常优于1%）定义了标准宇宙学模型。值得注意的是，Planck结果隐含$\Omega_\Lambda \approx 0.685$，$\Omega_m \approx 0.315$，且宇宙平坦（$\Omega_{\text{tot}} \approx 1$）。

重子声波振荡（BAO）提供了另一种独立的几何探针，其原理源于早期宇宙中光子-重子流体内的声波传播。在再复合之前，光子与重子强耦合，形成一种有效的“声学流体”。密度扰动在此流体中引发向外传播的声波。当电子与质子复合形成中性原子后，光子退耦，声波传播停止，在物质分布中留下一个特征尺度——声视界（sound horizon）$r_s$，即声波在光子退耦前能传播的最大距离。这个尺度在宇宙中表现为物质分布（如星系、类星体）的相关函数中出现的一个微弱的过量信号，即BAO特征。$r_s$在红移$z$处的角直径距离$D_A(z)$和径向哈勃距离$D_H(z) = c/H(z)$中表现为一个标准的“刻度”： $$\frac{r_s}{D_V(z)} \propto \text{BAO信号位置}$$ 其中$D_V(z) = \left[ (1+z)^2 D_A^2(z) \cdot \frac{cz}{H(z)} \right]^{1/3}$是一个组合距离。通过测量不同红移处星系巡天（如SDSS/BOSS/eBOSS）中BAO信号的峰值位置，可以确定$D_V(z)/r_s$。由于$r_s$由早期宇宙物理（由CMB精确测定）决定，BAO测量直接给出了宇宙膨胀历史$H(z)$和距离-红移关系$D_A(z)$，从而对$\Omega_m$、$\Omega_\Lambda$和$H_0$提供约束。BAO测量与CMB结合，能有效打破参数间的简并，特别是对曲率和暗能量状态的约束。例如，BAO数据确认了在低红移处宇宙的加速膨胀，并支持平坦ΛCDM模型。

除了几何探针，大尺度结构（LSS）的形成与演化也为ΛCDM模型提供了动力学检验。在物质主导时期，密度扰动在引力作用下增长。在线性微扰理论中，物质功率谱$P(k)$的演化由增长因子$D(z)$描述：$P(k, z) = P_{\text{prim}}(k) \cdot T^2(k) \cdot D^2(z)$，其中$P_{\text{prim}}(k) \propto k^{n_s}$是原初谱，$T(k)$是转移函数（编码了辐射-物质相等前后的物理），$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 = 0$$ 其中$\Omega_m(a) = \Omega_m a^{-3} / (\Omega_m a^{-3} + \Omega_\Lambda + \Omega_r a^{-4})$。在平坦ΛCDM中，增长因子可解析近似。测量星系红移空间 distortions（RSD）可推断$f\sigma_8(z)$，其中$f = d\ln D/d\ln a$是增长速率，$\sigma_8$是当前$8 h^{-1}$Mpc球内密度涨落的RMS。$f\sigma_8$对$\Omega_m$和引力理论敏感。当前LSS观测（如BOSS、KiDS、DES）给出的$f\sigma_8$测量与Planck ΛCDM预言在1-2σ内一致，但部分弱引力透镜测量稍低，存在轻微张力。

ΛCDM模型在解释观测数据方面取得了巨大成功，这使其成为“标准模型”。其成功主要体现在以下几个方面：

CMB功率谱的精确拟合：如前所述，Planck卫星观测到的CMB温度与偏振功率谱（从低$\ell$的 Sachs-Wolfe plateau到高$\ell$的声学峰系列及阻尼尾部）与基于6个基本ΛCDM参数的理论预言在统计上完美吻合。图1.1（概念图）展示了这一拟合。模型成功预测了：

• 声学峰的位置序列：第一峰在$\ell \approx 220$，后续峰等间距分布（间距约$\Delta \ell \approx 300$，对应声视界角尺度）。峰的位置精确对应平坦几何。
• 峰高的交替模式：奇偶峰高比由重子密度$\Omega_b h^2$决定，模型拟合值$\Omega_b h^2 \approx 0.0224$与大爆炸核合成（BBN）独立测得的轻元素丰度（如$^4\text{He}$、D）高度一致，这是宇宙学一个非凡的交叉验证。
• 极化谱的TE和EE相关：TE交叉谱中负峰的位置与温度谱正峰对齐，EE谱中峰高比与温度谱不同，这些特征均被模型准确预言。

BAO尺度的确认：从SDSS等巡天中测量的BAO信号位置，与基于CMB确定的$r_s$所预言的$D_V(z)/r_s$值高度一致。例如，在红移$z \approx 0.38, 0.51, 0.61$处的测量值与Planck ΛCDM预言吻合。这独立证实了CMB所给出的宇宙膨胀历史，并支持了平坦宇宙中暗能量主导的结论。

宇宙年龄的一致性：由最佳拟合参数积分弗里德曼方程得到的宇宙年龄$t_0 = \int_0^\infty \frac{dz}{(1+z)H(z)} \approx 13.80 \pm 0.02$ Gyr（Planck 2018），与最古老球状星团（约13.4 Gyr）和放射性衰变测年法（如铀-铅测年）的结果兼容，没有矛盾。

结构形成的定性成功：在ΛCDM框架下，通过N体数值模拟（如Millennium Simulation、IllustrisTNG），从 nearly scale-invariant 的原初扰动出发，能够自洽地生长出从星系到星系团、纤维状结构、宇宙空洞的宇宙大尺度结构网。模拟预言的结构统计量（如两点相关函数、质量函数）与巡天观测（如SDSS、2dF）定性吻合，尽管在最小尺度（矮星系数量、核心-陡坡问题）存在潜在矛盾，可能涉及重子物理反馈的复杂性。

然而，ΛCDM模型的“成功”表象下，隐藏着深刻的物理学困境。这些困境并非简单的观测偏差，而是指向模型基础假设的根本性质疑，主要集中在这一领域：

宇宙常数问题（真空能量灾难）：这是理论物理学中最 notorious 的不匹配。在量子场论中，真空并非空无一物，而是充斥着量子涨落，每种场都有其零点能。将这些零点能贡献求和（需截断于普朗克尺度），得到的真空能量密度估计值比观测到的暗能量密度$\rho_\Lambda \approx 6.91 \times 10^{-27} \text{ kg/m}^3$（对应$\Omega_\Lambda \approx 0.685$）大120个数量级（$10^{120}$倍）。即使考虑某种未知的对称性机制将真空能量精确抵消至零，为何剩余值恰好等于当前观测值（约$(2.3 \times 10^{-3} \text{ eV})^4$）？这个值在宇宙演化中并非恒定：在辐射主导时期，$\rho_\Lambda$远小于$\rho_r$；在物质主导时期，$\rho_\Lambda$开始接近$\rho_m$；恰在当下，二者可比，导致加速膨胀。这种“巧合”被称为“宇宙 coincidence problem”。ΛCDM模型将$\Lambda$视为一个常数，但并未解释其数值的微小与特殊 timing，这暗示可能需要对引力理论或真空本质有全新的理解。

哈勃常数（$H_0$）张力：这是当前最突出的观测层面矛盾。利用早期宇宙探针（CMB）在平坦ΛCDM假设下反推的$H_0$，与利用晚期宇宙局域距离阶梯法直接测量的$H_0$之间存在显著且持续的偏离。

• **CMB 推断值**：Planck 2018给出$H_0 = 67.4 \pm 0.5$ km/s/Mpc。此值依赖于ΛCDM模型框架，通过将观测到的CMB功率谱角尺度（由$\theta_* = r_s / D_A(z_*)$给出）与理论声视界$r_s$（由$\Omega_b h^2$, $\Omega_c h^2$等决定）比较来约束$D_A(z_*)$，进而结合其他参数得到$H_0$。
• **局域距离阶梯测量**：以SH0ES团队（Supernovae and H0 for the Equation of State）为代表，通过校准造父变星（Cepheids）距离，再用Ia型超新星（SNe Ia）将距离尺度延伸到遥远星系，最终测量哈勃常数。其最新结果为$H_0 = 73.0 \pm 1.0$ km/s/Mpc（Riess et al. 2022