宇宙学基础与标准模型

在探索暗宇宙的奥秘之前，我们必须建立坚实的宇宙学理论基础。现代宇宙学建立在两个核心支柱之上：爱因斯坦的广义相对论和宇宙学原理。本章将系统性地从宇宙学原理出发，推导出描述宇宙整体动力学的Friedmann方程，并定义Hubble参数、密度参数、曲率参数等关键概念。随后，我们将深入剖析当前的标准宇宙学模型——ΛCDM模型，详细阐释其三大组分——普通重子物质（约5%）、冷暗物质（约27%）和暗能量（约68%）——的物理本质及其对宇宙几何、膨胀历史和结构形成的决定性影响。通过梳理宇宙微波背景辐射（CMB）、重子声波振荡（BAO）等关键观测数据，我们将展示ΛCDM模型如何以惊人的精度拟合观测事实，并论证暗成分的引入并非理论臆测，而是描述我们宇宙的必然要求。最后，本章将为后续章节铺垫必要的数学物理工具，包括时空度规、能量动量张量、状态方程等，为深入探讨暗物质粒子物理和暗能量动力学扫清障碍。

1.1 宇宙学原理：现代宇宙学的基石

宇宙学原理（Cosmological Principle）是现代宇宙学的哲学与科学起点。它包含两个基本假设：均匀性（Homogeneity）——在足够大的尺度上（通常大于100 Mpc），宇宙的物质分布是均匀的，任意位置观测者看到的宇宙大尺度结构统计性质相同；各向同性（Isotropy）——在任意位置，宇宙在各个方向上的物理性质相同。这一原理的雏形可追溯至哥白尼的“日心说”所蕴含的“地球非中心”思想，但直到20世纪20年代，随着哈勃发现星系红移-距离关系（哈勃定律），宇宙学原理才被正式确立为构建宇宙模型的公理。

从数学上讲，均匀各向同性时空的度规必须具有Robertson-Walker（RW）度规的形式： $$ds^2 = -c^2 dt^2 + a(t)^2 \left[ \frac{dr^2}{1 - k r^2} + r^2 (d\theta^2 + \sin^2\theta d\phi^2) \right]$$ 其中 $a(t)$ 是标度因子（scale factor），描述宇宙随时间的膨胀；$k$ 是曲率常数，取值为 $+1, 0, -1$ 分别对应封闭、平坦、开放的宇宙几何；$r, \theta, \phi$ 是共动坐标（comoving coordinates），随宇宙膨胀而固定，实际物理距离由 $a(t)$ 缩放。RW度规的推导基于均匀各向同性群（空间旋转和三维平移对称性），是最大对称的三维空间形式，其唯一自由参数就是标度因子和曲率。

> 侧边栏：宇宙学原理的观测验证 > 宇宙学原理在 $z \lesssim 1$ 的局部宇宙中已被大量观测证实：星系巡天（如SDSS、2dF）显示，在尺度超过300 Mpc后，星系分布的涨落幅度低于1%；CMB温度各向异性在扣除偶极子后，涨落仅约 $10^{-5}$，且统计上高度各向同性。然而，在更大尺度（如超星系团尺度）或更高红移，均匀性假设仍需谨慎检验，这也是当前宇宙学前沿课题之一。

1.2 Friedmann方程：宇宙动力学的核心

将RW度规代入爱因斯坦场方程 $G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$，并假设宇宙物质可视为完美流体（perfect fluid），其能量动量张量为 $T_{\mu\nu} = (\rho + p/c^2) u_\mu u_\nu + p g_{\mu\nu}$，其中 $\rho$ 是能量密度，$p$ 是压强，$u^\mu$ 是四速度。通过计算Christoffel符号、Ricci张量并化简，可得到两个独立的动力学方程，称为Friedmann方程：

第一Friedmann方程（能量方程）： $$H^2(t) = \left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho(t) - \frac{k c^2}{a(t)^2} + \frac{\Lambda c^2}{3}$$ 第二Friedmann方程（加速度方程）： $$\frac{\ddot{a}}{a} = -\frac{4\pi G}{3} \left( \rho + \frac{3p}{c^2} \right) + \frac{\Lambda c^2}{3}$$

其中 $H(t)$ 是Hubble参数（Hubble parameter），表征宇宙膨胀速率。第一方程是能量守恒的体现，第二方程则直接反映引力与物质/能量成分的竞争。注意，若宇宙学常数 $\Lambda$ 被吸收进能量密度（定义 $\rho_\Lambda = \frac{\Lambda c^2}{8\pi G}$，$p_\Lambda = -\rho_\Lambda c^2$），则方程形式更统一。

由第一Friedmann方程，我们定义临界密度（critical density）： $$\rho_c(t) = \frac{3 H^2(t)}{8\pi G}$$ 这是使宇宙几何平坦（$k=0$）所需的密度。当前时刻 $t_0$ 的临界密度记为 $\rho_{c0}$，对应 $H_0$（Hubble常数）。观测值 $H_0 \approx 67.4$ km/s/Mpc（Planck 2018）给出 $\rho_{c0} \approx 8.5 \times 10^{-27}$ kg/m³，仅相当于每立方米几个氢原子。

进一步定义密度参数（density parameter）： $$\Omega_i(t) = \frac{\rho_i(t)}{\rho_c(t)}$$ 对于各组分 $i$（重子、暗物质、暗能量、辐射等），其总密度参数为 $\Omega_{\text{total}} = \sum_i \Omega_i$。将第一Friedmann方程两边除以 $H^2$，并引入曲率密度参数 $\Omega_k(t) = -\frac{k c^2}{a^2 H^2}$，可得宇宙学约束： $$\Omega_{\text{total}} + \Omega_k = 1$$ 这是一个几何恒等式，由RW度规和Einstein方程导出。若 $\Omega_{\text{total}} > 1$，则 $\Omega_k < 0$，$k=+1$，宇宙封闭；若 $\Omega_{\text{total}} < 1$，则 $\Omega_k > 0$，$k=-1$，宇宙开放；若 $\Omega_{\text{total}} = 1$，则 $\Omega_k=0$，$k=0$，宇宙平坦。

1.3 宇宙组分与状态方程

宇宙中不同组分的动力学行为由其状态方程（equation of state）决定：$p = w \rho c^2$，其中 $w$ 是无量纲状态方程参数。常见组分的 $w$ 值：

• **辐射**（光子、相对论性中微子）：$w = 1/3$，$\rho \propto a^{-4}$（因波长拉伸和数密度稀释）。
• **非相对论性物质**（重子、冷暗物质）：$w = 0$，$\rho \propto a^{-3}$（仅数密度稀释）。
• **宇宙学常数**（暗能量的一种形式）：$w = -1$，$\rho = \text{const}$。
• **一般暗能量**（如精质，Quintessence）：$w \approx -1$，可能随时间演化，$\rho \propto a^{-3(1+w)}$。

由连续性方程 $\dot{\rho} + 3H(\rho + p/c^2) = 0$ 可推导上述标度关系。状态方程 $w$ 直接影响宇宙膨胀历史：$w < -1/3$ 时，第二Friedmann方程中 $\ddot{a} > 0$，宇宙加速膨胀。

1.4 ΛCDM模型：标准宇宙学框架

ΛCDM（Lambda Cold Dark Matter）模型是当前描述宇宙演化的标准模型，其核心假设是：宇宙主要由冷暗物质（Cold Dark Matter, CDM）和宇宙学常数（Λ）主导，几何接近平坦，初始扰动为近尺度不变的绝热高斯涨落。该模型参数极少，却能以极高精度拟合从CMB到低红移的多种观测。

1.4.1 组分比例与物理本质

根据Planck 2018年CMB数据，ΛCDM模型的最佳拟合参数为（在 $H_0 = 67.4$ km/s/Mpc 下）：

• **重子物质**（Baryons, $\Omega_b h^2 \approx 0.0224$）：构成恒星、行星和我们的普通物质，占宇宙总质能约 **5%**。其物理本质已知，但重子丰度（$\Omega_b \approx 0.0486$）远低于总物质密度，表明存在大量非重子成分。
• **冷暗物质**（CDM, $\Omega_c h^2 \approx 0.120$）：占约 **27%**。所谓“冷”指其运动速度远低于光速（非相对论性），从而能在小尺度上形成引力势阱，驱动结构形成。其物理本质未知，候选者包括弱相互作用大质量粒子（WIMPs）、轴子（Axions）、惰性中微子等。CDM必须满足：不带电磁相互作用（故“暗”），与普通物质仅通过引力（可能还有弱相互作用）耦合，且在宇宙学时间尺度上运动缓慢。
• **暗能量**（Dark Energy, $\Omega_\Lambda \approx 0.6925$）：占约 **68%**。在ΛCDM中，暗能量被等效为宇宙学常数Λ，状态方程 $w = -1$，能量密度不随膨胀而稀释。其物理本质是最大谜团：可能是真空能（但理论预测值与观测值相差 $10^{120}$ 倍，即“宇宙学常数问题”），也可能是一种动态场（如精质）。观测表明 $w \approx -1$ 且随时间变化极小（$|w+1| < 0.05$），支持Λ的简单形式。

辐射组分（光子、中微子）在当今宇宙能量密度中可忽略（$\Omega_r \sim 10^{-4}$），但在早期宇宙（$z \gtrsim 3400$）占主导。

1.4.2 对宇宙几何与膨胀历史的决定性作用

总密度参数 $\Omega_{\text{total}} = \Omega_b + \Omega_c + \Omega_\Lambda + \Omega_r \approx 1.0007 \pm 0.0019$（Planck），表明宇宙极其接近平坦（$\Omega_k \approx 0$）。平坦性在暴胀宇宙学中自然产生，也是ΛCDM的关键预测之一。

不同组分的相对丰度直接决定了宇宙膨胀历史的分段行为： 1. 辐射主导期（$z \gtrsim 3400$）：$\rho_r \gg \rho_m, \rho_\Lambda$，$a(t) \propto t^{1/2}$。 2. 物质主导期（$3400 \gtrsim z \gtrsim 0.7$）：$\rho_m \gg \rho_r, \rho_\Lambda$，$a(t) \propto t^{2/3}$。此期间密度扰动增长，结构形成。 3. 暗能量主导期（$z \lesssim 0.7$）：$\rho_\Lambda \gg \rho_m, \rho_r$，$a(t) \propto e^{H_\Lambda t}$，其中 $H_\Lambda = H_0 \sqrt{\Omega_\Lambda} \approx 55$ km/s/Mpc，宇宙加速膨胀。

转折红移 $z_{\text{eq}}$（辐射-物质相等）和 $z_\Lambda$（物质-暗能量相等）由组分比例决定。例如，$z_\Lambda \approx 0.7$ 对应宇宙年龄约 $t \approx 0.8 t_0$（$t_0 \approx 13.8$ Gyr），即过去约60亿年宇宙进入加速期。

1.4.3 对结构形成的决定性作用

暗物质在结构形成中扮演“骨架”角色。其引力势阱在物质主导期增长，而普通物质（重子）因辐射压阻力，需在再复合（$z \approx 1100$）后才能落入势阱形成星系和星系团。若无冷暗物质：

• 仅靠重子物质，密度扰动增长受辐射压抑制，无法在观测时间内形成已见到的星系和星系团。
• CMB温度涨落功率谱中，第三峰之后的高度将急剧下降，与观测不符（观测显示多峰结构）。

ΛCDM中，暗物质密度 $\Omega_c$ 决定了物质功率谱的幅度和形状。重子物质 $\Omega_b$ 则调制声学振荡特征：在CMB功率谱中，峰与谷的相对高度比直接反映 $\Omega_b / \Omega_c$ 比例；在物质功率谱中，产生“重子声学振荡”（BAO）特征。

1.5 观测证据：ΛCDM的成功验证

ΛCDM模型并非凭空臆造，而是被多重独立观测以高精度证实。这些观测来自不同宇宙学探针，共同约束同一套参数，形成“交叉验证”的严密网络。

1.5.1 宇宙微波背景辐射（CMB）

CMB是再复合时期（$z \approx 1100$，$t \approx 380,000$ 年）遗留下来的黑体辐射，温度约2.725 K，各向异性仅约 $10^{-5}$。其温度涨落功率谱 $C_\ell$（角功率谱）对宇宙参数极其敏感：

• **第一峰位置**：主要决定宇宙曲率。观测峰值在 $\ell \approx 220$，对应角尺度约1°，精确符合平坦宇宙预测（$\Omega_{\text{total}} \approx 1$）。
• **峰高比**：第一峰与第二峰、第三峰的相对高度比，由 $\Omega_b h^2$ 和 $\Omega_c h^2$ 决定。重子物质增加第一峰高度并抑制后续峰；暗物质增加总物质密度，改变峰位置和整体幅度。
• **阻尼尾部**：小尺度（高 $\ell$）的衰减由光子扩散（Silk阻尼）造成，依赖于 $\Omega_b h^2$。

Planck卫星2018年发布的最终数据（基于TT, TE, EE, lowE组合）给出ΛCDM六参数拟合： | 参数 | 符号 | 值（Planck 2018） | |------|------|-------------------| | 重子密度 | $\Omega_b h^2$ | $0.02237 \pm 0.00015$ | | 冷暗物质密度 | $\Omega_c h^2$ | $0.1200 \pm 0.0012$ | | 角声学尺度 | $\theta_*$ | $1.04110 \pm 0.00030 \times 10^{-2}$ | | 再复合光学深度 | $\tau$ | $0.0544 \pm 0.0073$ | | 标量振幅 | $\ln(10^{10} A_s)$ | $3.044 \pm 0.014$ | | 标量谱指数 | $n_s$ | $0.9649 \pm 0.0042$ |

由此推导的 $\Omega_b \approx 0.0486$，$\Omega_c \approx 0.2589$，$\Omega_\Lambda \approx 0.6925$，$H_0 = 67.4 \pm 0.5$ km/s/Mpc。CMB的精度达到百分之一甚至千分之一水平，是ΛCDM最坚实的支柱。

1.5.2 重子声波振荡（BAO）

BAO是早期宇宙中声波在重子-光子流体中传播留下的特征尺度，约490 Mpc（共动），作为“标准尺”用于测量宇宙膨胀历史。通过测量星系或类星体分布中相关函数的峰值，可得到不同红移下的哈勃参数 $H(z)$ 和角直径距离 $D_A(z)$ 的组合（通常为 $D_V(z) = \left[ (1+z)^2 D_A^2(z) \frac{cz}{H(z)} \right]^{1/3}$）。

BAO测量（如BOSS、eBOSS、6dFGS）与CMB结合，能独立约束 $\Omega_m$ 和 $H_0$，且与CMB结果高度一致。例如，$z \approx 0.57$ 的BAO测量给出 $D_V(0.57)/r_d \approx 13.67$（$r_d$ 为声学半径），与ΛCDM预测吻合。BAO数据还支持平坦几何，且与超新星数据联合后，对暗能量状态方程 $w$ 的约束达 $w = -1.03 \pm 0.03$，与Λ完全一致。

1.5.3 超新星Ia（Type Ia Supernovae）

超新星Ia是来自白矮星热核爆炸的标准烛光，其峰值光度与光变曲线形状相关，可用于测量距离模数 $\mu = m - M$。1998年，两个团队（SCP、High-z SN）通过高红移（$z \sim 0.5-1$）超新星观测发现，宇宙在约50亿年前开始加速膨胀，直接证明了暗能量的存在。

在ΛCDM框架下，超新星数据约束 $\Omega_\Lambda$ 和 $\Omega_m$ 的degeneracy。结合CMB（固定 $\Omega_{\text{total}}$）后，超新星给出 $\Omega_\Lambda \approx 0.7$，$\Omega_m \approx 0.3$。当前大型项目（如Pantheon+，包含1700多颗超新星）将 $w$ 的约束提升至 $w = -1.013 \pm 0.025$（联合BAO和CMB），仍与 $w=-1$ 无显著偏差。

1.5.4 其他独立探针

• **引力透镜**：弱引力透镜测量物质分布（包括暗物质）的总量，与CMB预测的结构增长历史一致。
• **星系团丰度**：通过X射线、SZ效应或光学 richness 测量的星系团数密度随红移变化，对 $\Omega_m$ 和 $\sigma_8$（扰动幅度）敏感，与CMB一致。
• **局部 $H_0$ 测量**：造父变星+超新星校准给出 $H_0 \approx 73$ km/s/Mpc，与CMB推断值存在约 $4-6\sigma$ 张力。此“哈勃张力”可能暗示新物理（如早期暗能量），但尚未动摇ΛCDM整体框架，因其他探针（BAO、CMB）高度自洽。

1.6 暗成分的必然性：为何不能没有暗物质和暗能量？

1.6.1 暗物质的不可替代性

若假设宇宙仅含重子物质（$\Omega_b \approx 0.05$）和辐射，将无法解释以下观测： 1. 星系旋转曲线：旋涡星系外围恒星和气体的旋转速度远高于由可见物质引力预测值，需额外质量（暗物质晕）提供引力。 2. 星系团动力学：通过星系速度弥散或X射线热气体压力测量的星系团质量，远大于可见质量（维里定理）。 3. 引力透镜：大质量透镜（如子弹星系团）的引力场分布与可见物质分离，直接显示暗物质存在。 4. CMB功率谱：若 $\Omega_c = 0$，仅靠 $\Omega_b$ 无法产生观测到的第三峰高度；且第一峰位置要求 $\O