第一章：宇宙的组成 / Chapter 1: The Composition of the Universe

当我们仰望星空，看到的仅是宇宙的冰山一角。现代宇宙学揭示，我们熟悉的恒星、行星和气体仅占宇宙总质能的5%，而神秘的暗物质和暗能量分别贡献了27%和68%。这一惊人比例并非凭空猜测，而是基于宇宙微波背景辐射、大尺度结构和超新星观测的精密测量。本章将带你探索如何从这些观测数据中提取宇宙组成的密码，并建立描述宇宙演化的数学框架——宇宙学原理与弗里德曼方程。

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

原理的内容与哲学基础

宇宙学原理 (Cosmological Principle) 是现代宇宙学的核心假设，它断言：在足够大的尺度上（通常指数亿光年以上），宇宙是均匀 (homogeneous) 且各向同性 (isotropic) 的。均匀性意味着宇宙的物质分布不随位置改变；各向同性意味着从任何方向观察，宇宙的统计性质相同。这一原理并非凭空想象，而是对哥白尼原理 (Copernican Principle) 的推广——地球并非宇宙的特殊中心，进而推广到宇宙本身没有中心或边缘。

哲学上，宇宙学原理体现了“平庸原理” (Mediocrity Principle)：我们所在的时空位置不具有代表性。历史可追溯至牛顿的绝对空间概念被抛弃，爱因斯坦在建立广义相对论时，默认宇宙在大尺度上简单而规则。1922年，亚历山大·弗里德曼 (Alexander Friedmann) 首次在爱因斯坦场方程中引入这一假设，推导出动态宇宙解，为现代大爆炸理论铺平道路。

观测证据支持

尽管宇宙在小尺度上呈现显著结构（星系、星系团、空洞），但大尺度巡天提供了均匀性与各向同性的有力证据：

• **星系巡天**：如2度角场星系红移巡天 (2dFGRS) 和斯隆数字天空巡天 (SDSS) 显示，当尺度超过约1亿光年时，星系分布变得均匀。相关函数趋于零，功率谱在大尺度上平坦。
• **宇宙微波背景辐射 (CMB)**：COBE、WMAP和普朗克卫星测量到的CMB温度涨落仅约十万分之一（$\Delta T/T \sim 10^{-5}$），且涨落模式统计上各向同性，完美支持宇宙学原理。若存在全局各向异性，CMB功率谱将出现非物理的偶极或四极异常，但实际数据高度符合各向同性预测。

对宇宙学模型的约束

宇宙学原理直接导致弗里德曼-勒梅特-罗伯逊-沃尔克度规 (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=+1,0,-1$ 分别对应闭合、平坦、开放宇宙）。该度规的均匀各向同性性质，使得宇宙动力学完全由哈勃参数 $H(t) = \dot{a}/a$ 描述，极大简化了爱因斯坦场方程的求解。

> 侧边栏：宇宙学原理的局限性 > 近年研究探讨宇宙学原理是否在最大尺度上严格成立。例如，对CMB冷斑 (Cold Spot) 的分析、星系巡天中可能存在的超大尺度结构（如“宇宙墙”）引发争议。但主流观点认为，现有观测仍支持宇宙学原理在 $>100$ Mpc尺度上有效，任何偏离需更精确数据验证。

弗里德曼方程：宇宙动力学的数学描述

从爱因斯坦场方程到弗里德曼方程

在FLRW度规下，爱因斯坦场方程 $G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$ 简化为两个独立方程，其中第一个即弗里德曼方程 (Friedmann Equation)： $$ H^2 = \left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho - \frac{kc^2}{a^2} + \frac{\Lambda c^2}{3} $$ 这里 $\rho$ 是宇宙总能量密度（包括物质、辐射等），$k$ 是曲率参数，$\Lambda$ 是宇宙学常数。该方程将宇宙膨胀率 $H$ 与能量密度、曲率和暗能量联系起来，是宇宙学的心脏。

第二个方程是加速度方程： $$ \frac{\ddot{a}}{a} = -\frac{4\pi G}{3} \left( \rho + \frac{3p}{c^2} \right) + \frac{\Lambda c^2}{3} $$ 其中 $p$ 是压强。辐射主导时 $p = \rho c^2/3$，物质主导时 $p \approx 0$，暗能量主导时 $p \approx -\rho c^2$（宇宙学常数情形）。

关键参数：临界密度与密度参数

定义临界密度 (Critical Density) 为宇宙平坦 ($k=0$) 且无 $\Lambda$ 时的密度： $$ \rho_c = \frac{3H^2}{8\pi G} $$ 当前哈勃常数 $H_0 \approx 70$ km/s/Mpc（普朗克2018年给出 $H_0 = 67.4$ km/s/Mpc），计算得 $\rho_c \approx 8.5 \times 10^{-27}$ kg/m³（约5个氢原子每立方米）。引入密度参数 (Density Parameter)： $$ \Omega_i = \frac{\rho_i}{\rho_c}, \quad \Omega_\Lambda = \frac{\Lambda c^2}{3H_0^2}, \quad \Omega_k = 1 - \sum \Omega_i - \Omega_\Lambda $$ 其中 $\Omega_i$ 对应普通物质 (baryons)、暗物质 (cold dark matter) 等。弗里德曼方程可改写为： $$ \sum \Omega_i + \Omega_\Lambda + \Omega_k = 1 $$ 这直接给出宇宙的几何：$\Omega_k=0$ 为平坦，$\Omega_k>0$ 为开放，$\Omega_k<0$ 为闭合。

宇宙演化简史与成分主导

根据观测，宇宙演化经历不同主导成分阶段：

• **辐射主导** ($z \gg 3400$)：$\rho_r \propto a^{-4}$，$a \propto t^{1/2}$
• **物质主导** ($3400 > z > 0.5$)：$\rho_m \propto a^{-3}$，$a \propto t^{2/3}$
• **暗能量主导** ($z < 0.5$)：$\rho_\Lambda = \text{const}$，$a \propto e^{Ht}$（加速膨胀）

当前 $\Omega_\Lambda \approx 0.68$，$\Omega_m \approx 0.32$（其中普通物质 $\Omega_b \approx 0.05$，暗物质 $\Omega_c \approx 0.27$），故处于暗能量主导的加速膨胀阶段。

计算示例：从哈勃常数到临界密度

假设 $H_0 = 70$ km/s/Mpc，需转换为 SI 单位： $$ H_0 = 70 \times \frac{1000 \text{ m}}{1 \text{ km}} \times \frac{1 \text{ Mpc}}{3.086 \times 10^{22} \text{ m}} \approx 2.27 \times 10^{-18} \text{ s}^{-1} $$ 代入 $\rho_c = 3H_0^2/(8\pi G)$，其中 $G = 6.674 \times 10^{-11}$ m³ kg⁻¹ s⁻²： $$ \rho_c = \frac{3 \times (2.27 \times 10^{-18})^2}{8\pi \times 6.674 \times 10^{-11}} \approx 9.2 \times 10^{-27} \text{ kg/m}^3 $$ 若观测得普通物质密度 $\rho_b = 4.6 \times 10^{-28}$ kg/m³，则 $\Omega_b = \rho_b / \rho_c \approx 0.05$，符合现代值。

观测宇宙学：揭示宇宙组成的三大支柱

宇宙组成的精确比例并非来自单一观测，而是宇宙学标准模型 (ΛCDM模型) 对多种独立数据集的联合拟合。三大支柱——CMB、大尺度结构（BAO）和超新星——相互验证，形成严密证据链。

宇宙微波背景辐射：婴儿宇宙的 snapshot

黑体辐射与再复合

大爆炸后约38万年，宇宙冷却至约3000K，电子与质子结合形成中性氢（再复合），光子得以自由传播，形成今天观测到的CMB。其能谱完美符合温度为 $T_0 = 2.72548 \pm 0.00057$ K 的黑体辐射，这是大爆炸理论的 triumphs。

温度涨落与角功率谱

CMB并非完全均匀，存在约 $10^{-5}$ 量级的温度涨落 $\Delta T(\theta, \phi)$。这些涨落源于早期宇宙的量子涨落，经暴胀放大，成为星系形成的种子。涨落统计性质由角功率谱 $C_\ell$ 描述： $$ \langle a_{\ell m} a_{\ell' m'}^* \rangle = C_\ell \delta_{\ell\ell'} \delta_{mm'} $$ 其中 $a_{\ell m}$ 是球谐展开系数。$C_\ell$ 随 $\ell$（相当于角尺度 $\theta \sim 180^\circ/\ell$）变化，呈现一系列“声学峰”。

声学峰与参数约束

在光子-重子流体中，引力坍缩与辐射压力产生声波振荡，在再复合时刻“冻结”为温度涨落模式。第一峰位置对应声学尺度 $\theta_*$，与宇宙几何直接相关： $$ \theta_* = \frac{r_s(z_*)}{D_A(z_*)} $$ 其中 $r_s(z_*)$ 是再复合时的声学半径，$D_A(z_*)$ 是角直径距离。平坦宇宙下，$\theta_* \approx 0.6^\circ$，对应 $\ell \approx 200$。普朗克卫星测得第一峰在 $\ell \approx 220$，强烈支持宇宙平坦（$\Omega_k \approx 0$）。

峰的高度和序列敏感于：

• **普通物质密度 $\Omega_b h^2$**：增加重子质量增强声波惯性，提高奇数峰（压缩峰）相对偶数峰（稀疏峰）的高度。
• **冷暗物质密度 $\Omega_c h^2$**：影响引力势阱深度，改变峰的整体幅度和阻尼尾部。
• **光谱指数 $n_s$**：描述初始涨落谱形状。
• **光学深度 $\tau$**：再电离导致的大尺度功率抑制。

通过拟合整个 $C_\ell$ 谱（至 $\ell \sim 2500$），可精确提取参数。普朗克2018年结果（TT,TE,EE+lowE）： $$ \Omega_b h^2 = 0.02237 \pm 0.00015, \quad \Omega_c h^2 = 0.1200 \pm 0.0012 $$ 结合 $H_0$ 先验，得 $\Omega_b \approx 0.049$，$\Omega_c \approx 0.265$，$\Omega_\Lambda \approx 0.686$（假设平坦）。

> 交互提示：访问普朗克官网 (https://www.cosmos.esa.int/web/planck) 查看 $C_\ell$ 谱数据与拟合曲线，拖动参数滑块观察曲线变化。

大尺度结构：宇宙的骨架

星系巡天与相关函数

星系巡天（如SDSS、BOSS、eBOSS）绘制三维星系分布图，计算两点相关函数 $\xi(r)$ 或功率谱 $P(k)$。在约100 Mpc尺度上，出现重子声波振荡 (Baryon Acoustic Oscillation, BAO) 峰，即早期声波在物质分布中留下的“化石”标准尺。

BAO作为标准尺

BAO特征尺度即为声学半径 $r_d \approx 147$ Mpc（由CMB精确测定）。在星系关联中，沿和横越视线方向的峰位置给出： $$ \frac{r_d}{D_V(z)} = \frac{r_d}{\left[ (1+z)^2 D_A^2(z) \frac{cz}{H(z)} \right]^{1/3}} $$ 其中 $D_V$ 是组合距离。测量不同红移 $z$ 的BAO尺度，可约束 $H(z)$ 和 $D_A(z)$，进而限制 $\Omega_m$ 和 $\Omega_\Lambda$。

例如，SDSS-III BOSS在 $z=0.57$ 测量得 $D_V/r_d = 13.67 \pm 0.22$，与ΛCDM预测一致。BAO数据对 $\Omega_m$ 敏感，尤其当结合CMB时，能打破 $\Omega_m$-$\Omega_\Lambda$ 简并。

全功率谱分析

除BAO峰外，整个 $P(k)$ 形状也含丰富信息：

• 小尺度（$k \gg 0.1$ h/Mpc）受重子声学特征和星系形成物理影响。
• 大尺度（$k \ll 0.1$ h/Mpc）形状由转移函数决定，依赖于 $\Omega_m h$ 和 $\Omega_b h^2$。

通过将观测 $P(k)$ 与理论预言（使用Boltzmann代码如CAMB或CLASS）比较，可独立约束 $\Omega_m$ 和 $\Omega_b$，与CMB结果交叉验证。

超新星Ia：宇宙膨胀的加速器

作为标准烛光

Ia型超新星来自碳氧白矮星的热核爆炸，峰值亮度高度一致（经光变曲线形状校正后，亮度弥散约0.15星等），可作为标准烛光 (Standard Candle)。通过比较视亮度与红移，测量距离模数 $\mu = m - M$，进而得到 luminosity distance $d_L$： $$ \mu = 5 \log_{10}(d_L/\text{Mpc}) + 25 $$

发现加速膨胀

1998年，Supernova Cosmology Project 和 High-z Supernova Search Team 独立发现，高红移（$z \sim 0.5-1$）超新星比减速膨胀模型预言更暗，意味着宇宙在近期（$z<1$）膨胀加速。这要求存在负压强成分，即暗能量 (Dark Energy)，其状态方程参数 $w \equiv p/\rho c^2 \approx -1$（宇宙学常数情形）。

约束暗能量

超新星数据直接测量 $d_L(z)$，对 $\Omega_\Lambda$ 和 $w$ 敏感。单独超新星给出 $\Omega_\Lambda \approx 0.7$，但存在 $\Omega_m$-$\Omega_\Lambda$ 简并。联合CMB（给出 $\Omega_m h^2$）和BAO（给出 $D_A(z)$）可打破简并，精确确定 $\Omega_\Lambda$ 并限制 $w$。当前数据与 $w=-1$ 一致，未发现动力学暗能量（如精质场 quintessence）证据。

> 案例研究：Pantheon+数据集 > 包含1701颗Ia型超新星，覆盖 $0.01 < z < 2.3$。联合Planck CMB和BAO，得 $\Omega_\Lambda = 0.688 \pm 0.010$，$\Omega_m = 0.299 \pm 0.008$（假设平坦）。若允许 $w \neq -1$，$w = -1.03 \pm 0.04$，与宇宙学常数一致。

综合数据分析：从多重观测到精确比例

联合拟合与参数简并

单一观测存在参数简并：例如，CMB第一峰位置主要约束 $\Omega_k$，但对 $\Omega_\Lambda$ 不敏感（因 $\Omega_\Lambda$ 在 $z \sim 1100$ 可忽略）；超新星对 $\Omega_\Lambda$ 敏感但依赖 $\Omega_m$。因此，必须联合分析 多个数据集，利用不同红移探针的互补性。

标准流程： 1. 使用Boltzmann代码（如CAMB）计算给定 $\{\Omega_b h^2, \Omega_c h^2, \Omega_\Lambda, H_0, n_s, A_s\}$ 的理论预言（CMB功率谱、BAO距离、物质功率谱等）。 2. 构建似然函数 $\mathcal{L}(\text{数据}|\text{参数})$，比较理论与观测。 3. 用马尔可夫链蒙特卡洛 (MCMC) 方法采样参数空间，得到后验分布。

当前最佳估计：普朗克2018结果

基于Planck 2018释放数据，结合BAO和超新星，平坦ΛCDM模型的最佳拟合值（TT,TE,EE+lowE+BAO）： | 参数 | 符号 | 值 | 95% 置信区间 | |------|------|-----|--------------| | 普通物质密度 | $\Omega_b h^2$ | $0.02237 \pm 0.00015$ | [0.02207, 0.02267] | | 冷暗物质密度 | $\Omega_c h^2$ | $0.1200 \pm 0.0012$ | [0.1176, 0.1224] | | 哈勃常数 | $H_0$ (km/s/Mpc) | $67.4 \pm 0.5$ | [66.4, 68.4] | | 重子密度 | $\Omega_b$ | $0.0486 \pm 0.0006$ | [0.0474, 0.0498] | | 暗物质密度 | $\Omega_c$ | $0.2589 \pm 0.0057$ | [0.2475, 0.2703] | | 暗能量密度 | $\Omega_\Lambda$ | $0.6925 \pm 0.0073$ | [0.6779, 0.7071] | | 物质密度 | $\Omega_m$ | $0.315 \pm 0.007$ | [0.301, 0.329] | | 曲率参数 | $\Omega_k$ | $0.0007 \pm 0.0019$ | [-0.0029, 0.0043] |

近似取整：普通物质 5%，暗物质 27%，暗能量 68%，宇宙平坦（$\Omega_k \approx 0$）。

系统误差与争议

尽管ΛCDM模型极其成功，仍存在挑战：

• **哈勃张力**：CMB推导的 $H_0$（~67 km/s/Mpc）与局域测量（如SH0ES使用造父变星校准的Ia型超新星，$H_0