宇宙学标准模型与隐匿组分的提出背景

1929年，埃德温·哈勃（Edwin Hubble）通过分析24个星系的光谱红移与造父变星校准的距离，首次发现了哈勃定律（Hubble's law）：星系的退行速度$v$与距离$d$成正比，比例系数即哈勃常数$H_0$，这一发现直接颠覆了此前静态宇宙的认知，为大爆炸宇宙学（Big Bang cosmology）奠定了观测基础[^1]。此后近一个世纪的观测进展不断验证并完善了这一框架，却也揭示出标准模型中无法回避的组分缺失矛盾，最终引出了暗物质（dark matter）与暗能量（dark energy）这两类占据宇宙总能量95%的隐匿组分。

1.1 大爆炸宇宙学的核心框架与弗里德曼方程

现代宇宙学的理论基石建立在两个基本假设之上：其一是宇宙学原理（cosmological principle），即宇宙在大尺度上（>100 Mpc）满足均匀性（homogeneity）与各向同性（isotropy），不存在特殊的中心或方向；其二是广义相对论（general relativity） 作为引力的正确描述，其核心是爱因斯坦场方程（Einstein field equations）： $$G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$$ 其中$G_{\mu\nu}$为爱因斯坦张量，描述时空的曲率；$\Lambda$为宇宙学常数；$g_{\mu\nu}$为时空度规；$T_{\mu\nu}$为物质能量动量张量；$G$为万有引力常数，$c$为光速。

在宇宙学原理的约束下，时空度规可简化为弗里德曼-勒梅特-罗伯逊-沃尔克度规（FLRW metric）： $$ds^2 = -c^2 dt^2 + a^2(t) \left[ \frac{dr^2}{1-kr^2} + r^2 (d\theta^2 + \sin^2\theta d\phi^2) \right]$$ 其中$a(t)$为标度因子（scale factor），描述宇宙的相对膨胀尺度；$k$为曲率参数，取值为$+1$（闭合宇宙）、$0$（平坦宇宙）或$-1$（开放宇宙），对应不同的时空几何。

将FLRW度规代入爱因斯坦场方程，可推导出两个独立的动力学方程，即弗里德曼方程（Friedmann equations）： :::info 📖 纸书侧边栏：弗里德曼方程的推导要点 1. 计算FLRW度规下的克里斯托费尔符号（Christoffel symbols） 2. 求解里奇张量（Ricci tensor）$R_{\mu\nu}$与里奇标量（Ricci scalar）$R$ 3. 代入爱因斯坦张量$G_{\mu\nu} = R_{\mu\nu} - \frac{1}{2} R g_{\mu\nu}$ 4. 结合理想流体的能量动量张量$T_{\mu\nu} = (\rho + p/c^2) u_\mu u_\nu + p g_{\mu\nu}$，得到时空与物质耦合的演化方程 ::: 第一个方程为哈勃参数演化方程： $$\left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3}$$ 其中$\dot{a} = da/dt$为标度因子对宇宙时$t$的导数，$\dot{a}/a$定义为哈勃参数$H(t)$，描述宇宙的膨胀速率；$\rho$为宇宙总能量密度；$p$为总压强。该方程表明，宇宙的膨胀速率由物质密度、曲率与宇宙学常数共同决定。

第二个方程为加速度方程： $$\frac{\ddot{a}}{a} = -\frac{4\pi G}{3} \left( \rho + \frac{3p}{c^2} \right) + \frac{\Lambda c^2}{3}$$ 其中$\ddot{a} = d^2a/dt^2$为标度因子的二阶导数。该方程揭示了宇宙膨胀的加速或减速机制：当物质与压强的贡献$\rho + 3p/c^2$为正时，引力会抑制膨胀，导致$\ddot{a}<0$；当宇宙学常数或负压强组分的贡献超过物质项时，会出现$\ddot{a}>0$的加速膨胀现象。

为定量描述不同组分对宇宙总能量密度的贡献，引入密度参数（density parameter） 的定义： $$\Omega_i = \frac{\rho_i}{\rho_c}, \quad \rho_c = \frac{3H^2}{8\pi G}$$ 其中$\rho_c$为临界密度（critical density），对应宇宙恰好平坦（$k=0$）时的总能量密度，当前（$H=H_0$）的临界密度约为$\rho_{c0} \approx 9.2 \times 10^{-27} \text{ kg m}^{-3}$，相当于每立方米空间中分布约5个氢原子的质量。总密度参数定义为$\Omega_{\text{total}} = \sum_i \Omega_i + \Omega_k$，其中$\Omega_k = -kc^2/(a^2 H^2)$为曲率密度参数，反映时空的几何偏离平坦的程度。

不同宇宙组分的密度参数随标度因子$a$的演化规律由其状态方程$p = w \rho c^2$决定：将连续方程$\dot{\rho} + 3H(\rho + p/c^2) = 0$与弗里德曼方程联立，可得$\Omega_i(a) \propto a^{-3(1+w)}$。对于非相对论性物质（重子、暗物质），$w \approx 0$，因此$\rho_m \propto a^{-3}$；对于辐射（光子、中微子），$w = 1/3$，因此$\rho_r \propto a^{-4}$（额外的$a^{-1}$因子来自光子波长的拉伸）；对于宇宙学常数（或更一般的精质quintessence），$w = -1$，因此$\rho_\Lambda = \text{const}$，不随宇宙膨胀而稀释。这一演化规律是理解宇宙组分随时间的核心依据：在宇宙诞生后的前约4.7万年，辐射主导宇宙膨胀；此后至约98亿年前，物质（重子+暗物质）主导膨胀；最近的约50亿年，暗能量主导，导致宇宙加速膨胀。

1.2 20世纪宇宙学观测的关键里程碑

大爆炸宇宙学的理论框架在20世纪通过一系列关键观测得到了反复验证，这些观测不仅确立了标准模型的地位，也为后续暗物质与暗能量的发现埋下了伏笔。

1.2.1 哈勃膨胀的精确测量

哈勃1929年的原始结果仅基于24个近邻星系，测量误差较大，甚至曾因距离校准的错误导致哈勃常数被高估了近10倍。此后数十年间，天文学家通过造父变星周光关系、Ia型超新星（Type Ia supernova）、星系表面亮度起伏（surface brightness fluctuation）等多种测距方法，逐步扩展了可观测的宇宙深度。到20世纪末，已有数十万个星系的距离与红移被测量，哈勃定律在红移$z \ll 1$的线性区域得到了极高精度的验证，当前的哈勃常数测量精度已达到1%量级：根据普朗克卫星（Planck satellite）的宇宙微波背景辐射测量结果，$H_0 = 67.4 \pm 0.5 \text{ km s}^{-1} \text{ Mpc}^{-1}$；而基于近邻造父变星与Ia型超新星的本地测量结果为$H_0 = 73.2 \pm 1.3 \text{ km s}^{-1} \text{ Mpc}^{-1}$，二者之间的“哈勃张力（Hubble tension）”是当前宇宙学的前沿争议之一，可能预示着标准模型之外的新物理[^2]。

:::info 📱 电子版交互元素 点击查看[普朗克卫星官网的哈勃常数测量结果](https://www.cosmos.esa.int/web/planck/publications)，或观看[哈勃定律发现历程的科普视频](https://www.youtube.com/watch?v=tiJcTuqQmSQ)。 :::

1.2.2 原初核合成对重子密度的约束

大爆炸核合成（Big Bang Nucleosynthesis, BBN）发生在大爆炸后1秒至3分钟的时间窗口内，此时宇宙温度从$10^{10}\text{ K}$降至$10^9\text{ K}$，高能光子与中微子退耦，质子和中子开始结合形成轻原子核。BBN的核心预言是：宇宙中轻元素（氘、氦-3、氦-4、锂-7）的丰度仅由宇宙早期的重子-光子比$\eta = n_b / n_\gamma$决定，其中$n_b$为重子数密度，$n_\gamma$为光子数密度。

当前的天文观测对BBN预言的验证达到了极高的精度：观测到的氦-4质量分数为$Y_p = 0.245 \pm 0.003$，氘丰度$(D/H) = (2.527 \pm 0.030) \times 10^{-5}$，锂-7丰度$(^7\text{Li}/\text{H}) = (1.6 \pm 0.3) \times 10^{-10}$，与BBN的理论计算高度吻合[^3]。通过匹配观测丰度与理论预言，可得到重子密度参数的精确约束： $$\Omega_b h^2 = 0.0224 \pm 0.0001$$ 其中$h = H_0 / (100 \text{ km s}^{-1} \text{ Mpc}^{-1})$为无量纲哈勃常数，当前$h \approx 0.674$，因此$\Omega_b \approx 0.049$，即重子物质仅占宇宙总能量密度的约5%。这一结果直接表明，宇宙中绝大部分物质并非由质子、中子等重子组成，为暗物质的存在提供了第一个间接证据。

1.2.3 宇宙微波背景辐射的各向异性测量

1965年，阿尔诺·彭齐亚斯（Arno Penzias）与罗伯特·威尔逊（Robert Wilson）在调试射电望远镜时，意外发现了存在于全天空、温度为$2.725 \pm 0.001\text{ K}$的各向同性微波辐射，即宇宙微波背景辐射（Cosmic Microwave Background, CMB），这是大爆炸宇宙学最直接的证据之一：该辐射是宇宙诞生后约38万年，电子与质子复合为中性氢，光子不再被散射而自由传播留下的“余晖”。

1992年，宇宙背景探测者卫星（COBE）首次测量到CMB的各向异性，其温度涨落幅度仅为$10^{-5}$量级，说明早期宇宙的密度分布高度均匀，为后续星系、星系团等大尺度结构的形成提供了初始条件。2001年发射的威尔金森微波各向异性探测器（WMAP）和2009年发射的普朗克卫星，对CMB的温度与偏振各向异性进行了全天空高精度扫描，测量到了角功率谱上的多组声学峰（acoustic peaks）。

CMB角功率谱的测量可同时约束宇宙的几何平坦性、重子密度与暗物质密度：第一个声学峰的位置对应宇宙的总密度参数$\Omega_{\text{total}} = 1.0007 \pm 0.0019$，与平坦宇宙的预言完全一致，说明时空曲率可忽略；声学峰的相对振幅则由重子与暗物质的密度共同决定：重子物质会增强第一个峰的振幅，而暗物质会抑制第二个及后续峰的振幅。结合BBN的重子密度约束，CMB测量给出了暗物质密度的精确值$\Omega_{\text{DM}} h^2 = 0.120 \pm 0.001$，对应$\Omega_{\text{DM}} \approx 0.264$。

:::info 📖 纸书侧边栏：CMB声学峰的物理起源 光子与重子流体在早期宇宙中耦合，形成以声速传播的声波：在光子脱耦时刻（$z \approx 1100$），声波传播

在光子脱耦时刻（$z \approx 1100$），声波传播的最大距离对应第一个声学峰的角尺度，其位置对宇宙总密度高度敏感；而重子与暗物质的相对密度决定了声学振荡的阻尼程度与峰振幅的比值。普朗克卫星的高精度测量显示，第一个声学峰的位置对应角尺度$\theta_* \approx 0.596^\circ$，与平坦宇宙的预言完全一致，验证了$\Omega_{\text{total}} \approx 1$的结论；同时，第二个声学峰的振幅被显著抑制，直接证明了宇宙中存在不与光子耦合的冷暗物质（cold dark matter, CDM）组分——若仅有重子物质，第二个峰的振幅会比观测值高约30%，这一差异在普朗克数据中以超过$50\sigma$的显著性被确认，是暗物质存在的最直接证据之一。

```mermaid xychart-beta title "宇宙组分密度参数随红移演化" x-axis "红移 z" [0, 0.5, 1, 2, 5, 10, 1100] y-axis "密度参数 $\Omega_i h^2$" 0 --> 0.5 line "重子物质" [0.022, 0.055, 0.088, 0.176, 0.44, 0.88, 44] line "冷暗物质" [0.12, 0.30, 0.48, 0.96, 2.4, 4.8, 240] line "暗能量" [0.68, 0.645, 0.432, 0.204, 0.056, 0.016, ~0] ```

上述三类里程碑观测共同构建了大爆炸宇宙学的标准框架，但仅用重子物质无法解释一系列跨尺度的引力效应观测结果，这些矛盾最终引出了暗物质与暗能量的假说。

1.3 标准宇宙学模型的观测矛盾

1.3.1 星系旋转曲线异常

星系旋转曲线（galaxy rotation curve）是指星系中恒星、电离氢区等发光体的旋转速度$v(R)$随到星系中心距离$R$的变化关系。在经典力学框架下，若星系的总质量$M(R)$集中在发光区域（盘面与核球），则根据开普勒第三定律，$v(R) = \sqrt{GM(R)/R}$，在发光半径$R_{\text{eff}}$之外，$M(R)$近似为常数，因此$v(R) \propto R^{-1/2}$，即旋转曲线应随距离

在发光半径$R_{\text{eff}}$之外，若质量仅来自可见的重子物质，则$M(R)$近似为常数，因此$v(R) \propto R^{-1/2}$，即旋转曲线应随距离增加而下降。然而，1933年天文学家弗里茨·兹威基（Fritz Zwicky）在研究后发座星系团（Coma Cluster）中星系的运动时，首次注意到引力效应与可见物质质量的严重偏离，但这一发现直到1970年代才因维拉·鲁宾（Vera Rubin）对仙女座星系（M31）及数十个螺旋星系的旋转曲线系统测量而得到广泛关注。鲁宾团队使用基特峰国家天文台（KPNO）的4米望远镜，测量了M31中电离氢区（HII区）的视向速度，发现即使在距离星系中心30 kpc的盘面外侧，旋转速度仍保持在约250 km/s，而非预期的下降至100 km/s以下。这一平缓旋转曲线（flat rotation curve）现象在后续的数十个螺旋星系、矮不规则星系、透镜星系中均被观测到，其普遍性直接否定了“质量仅集中于可见盘面”的假设。

根据牛顿力学，旋转速度的平方与 enclosed 质量的关系为$v^2(R) = \frac{GM(R)}{R}$，若$v(R)$为常数$v_0$，则$M(R) = \frac{v_0^2 R}{G}$，即 enclosed 质量随半径线性增加，远超过可见盘面的质量贡献。以银河系为例，其可见盘面的质量约为$6 \times 10^{10} M_\odot$，而根据旋转曲线推算的半径50 kpc处的 enclosed 质量约为$10^{12} M_\odot$，说明存在一个延伸至数百千秒差距的不可见质量晕，即暗物质晕（dark matter halo）。目前最广泛使用的暗物质晕密度剖面为纳瓦罗-弗伦克-怀特（Navarro-Frenk-White, NFW）剖面，由冷暗物质（Cold Dark Matter, CDM）的数值模拟得到： $$\rho_{\text{NFW}}(r) = \frac{\rho_s}{\frac{r}{r_s}\left(1 + \frac{r}{r_s}\right)^2}$$ 其中$r_s$为特征尺度半径，$\rho_s$为特征密度。NFW剖面的核心区域$\rho \propto r^{-1}$，外侧区域$\rho \propto r^{-3}$，能够很好地拟合大多数螺旋星系的旋转曲线，但也存在一些挑战：矮星系的观测显示其核心区域密度更接近常数（即“核心剖面”，core profile），而非NFW的尖点（cusp）剖面，这一**核心-尖点问题