宇宙学标准模型与观测基石 / The Cosmological Standard Model and Observational Pillars

当我们凝视星空，宇宙展现出一幅令人困惑的图景：95%的宇宙成分是不可见的暗物质与暗能量。这一惊人事实并非来自哲学思辨，而是建立在坚实的观测基石之上——宇宙学标准模型（ΛCDM model）及其关键参数的精确测定。本章将系统阐述这一模型的理论基础、核心方程推导，以及如何通过宇宙微波背景辐射（Cosmic Microwave Background, CMB）和重子声学振荡（Baryon Acoustic Oscillations, BAO）等观测手段，将抽象的宇宙学参数转化为可测量的物理量。我们将看到，现代宇宙学如何通过精密的卫星观测与大规模星系巡天，将宇宙的年龄、组成和演化历史锁定在令人惊叹的精度之内，为探索暗宇宙的本质铺平道路。

宇宙动力学的基本方程：从爱因斯坦到场方程简化

宇宙学标准模型的理论框架根植于爱因斯坦的广义相对论。其核心是爱因斯坦场方程： $$ G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu} $$ 其中 $G_{\mu\nu}$ 是爱因斯坦张量，$\Lambda$ 是宇宙学常数，$T_{\mu\nu}$ 是能量-动量张量。为了描述大尺度上均匀和各向同性的宇宙，我们必须对时空度规做出对称性假设。

FLRW度规与宇宙学原理

基于宇宙学原理（Cosmological Principle），即在足够大的尺度上，宇宙是均匀且各向同性的，时空度规必须具有弗里德曼-勒梅特-罗伯逊-沃尔克（Friedmann-Lemaître-Robertson-Walker, FLRW）形式： $$ 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$ 分别对应封闭、平坦、开放的宇宙几何。FLRW度规将问题简化为仅依赖于时间 $t$ 的尺度因子 $a(t)$，这是推导宇宙动力学方程的关键第一步。

能量-动量张量与完美流体假设

对于大尺度宇宙，物质和辐射可视为完美流体，其能量-动量张量为： $$ T^{\mu}_{\ \nu} = \text{diag}(-\rho c^2, p, p, p) $$ 其中 $\rho$ 是总能量密度（包括物质、辐射、暗能量等所有成分），$p$ 是压强。不同成分的状态方程（Equation of State）不同：

• **冷暗物质（CDM）与重子物质**：非相对论性，$p \approx 0$，故 $w \equiv p/(\rho c^2) = 0$。
• **辐射（光子、中微子）**：相对论性，$p = \rho c^2 / 3$，故 $w = 1/3$。
• **宇宙学常数（Λ）**：$p = -\rho c^2$，故 $w = -1$。
• **一般暗能量**：可能 $w \neq -1$，但标准模型取 $w = -1$。

总能量密度 $\rho = \sum_i \rho_i$，总压强 $p = \sum_i p_i$。在宇宙膨胀过程中，各成分的能量密度随尺度因子变化规律为： $$ \rho_i(a) = \rho_{i,0} \, a^{-3(1+w_i)} $$ 其中下标 $0$ 表示当前时刻值。因此：

• 物质（$w=0$）：$\rho_m \propto a^{-3}$（体积稀释）
• 辐射（$w=1/3$）：$\rho_r \propto a^{-4}$（体积稀释+红移）
• 宇宙常数（$w=-1$）：$\rho_\Lambda = \text{const.}$

弗里德曼方程的完整推导

将FLRW度规代入爱因斯坦场方程，经过计算克里斯托费尔符号、里奇张量 $R_{\mu\nu}$ 和曲率标量 $R$，可得非零的爱因斯坦张量分量： $$ G^0_{\ 0} = -3 \left( \frac{\dot{a}^2}{a^2} + \frac{k c^2}{a^2} \right) / c^2, \quad G^i_{\ j} = \left( -2 \frac{\ddot{a}}{a} - \frac{\dot{a}^2}{a^2} - \frac{k c^2}{a^2} \right) \delta^i_j / c^2 $$ 其中点表示对时间 $t$ 的导数。场方程 $G^0_{\ 0} + \Lambda = -\frac{8\pi G}{c^4} T^0_{\ 0} = \frac{8\pi G}{c^2} \rho$ 给出第一弗里德曼方程（哈勃方程）： $$ \left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho - \frac{k c^2}{a^2} + \frac{\Lambda c^2}{3} $$ 定义哈勃参数 $H(t) \equiv \dot{a}/a$，则： $$ H^2(t) = \frac{8\pi G}{3} \rho(t) - \frac{k c^2}{a(t)^2} + \frac{\Lambda c^2}{3} \tag{1} $$ 场方程的空间分量 $G^i_{\ j} = -\frac{8\pi G}{c^4} T^i_{\ j}$（$i \neq j$ 时自动满足，$i=j$ 时）给出第二弗里德曼方程（加速度方程）： $$ \frac{\ddot{a}}{a} = -\frac{4\pi G}{3} \left( \rho + \frac{3p}{c^2} \right) + \frac{\Lambda c^2}{3} \tag{2} $$ 方程（1）和（2）是宇宙动力学的核心。它们将宇宙的膨胀速率 $H$ 和加速度 $\ddot{a}$ 与宇宙的内容（$\rho$, $p$）、几何（$k$）和宇宙学常数 $\Lambda$ 联系起来。

密度参数与宇宙几何

定义临界密度（critical density）为恰好使 $k=0$ 的密度： $$ \rho_c(t) \equiv \frac{3 H(t)^2}{8\pi G} $$ 引入密度参数（density parameters）： $$ \Omega_i(t) \equiv \frac{\rho_i(t)}{\rho_c(t)}, \quad \Omega_k(t) \equiv -\frac{k c^2}{a(t)^2 H(t)^2}, \quad \Omega_\Lambda(t) \equiv \frac{\Lambda c^2}{3 H(t)^2} $$ 则第一弗里德曼方程（1）变为简洁的求和形式： $$ \sum_i \Omega_i(t) + \Omega_k(t) = 1 \tag{3} $$ 在当前时刻 $t_0$，有 $\Omega_{i,0} \equiv \Omega_i(t_0)$。观测表明宇宙非常接近平坦，即 $\Omega_k \approx 0$，故： $$ \Omega_{m,0} + \Omega_{r,0} + \Omega_{\Lambda,0} \approx 1 $$ 其中 $\Omega_m = \Omega_b + \Omega_c$，$\Omega_b$ 为重子密度，$\Omega_c$ 为冷暗物质密度。辐射密度 $\Omega_r$ 极小（$\sim 10^{-4}$），可忽略。

哈勃参数与红移的关系

在宇宙学中，红移 $z$ 定义为 $1+z = a_0 / a(t)$，其中 $a_0 \equiv a(t_0) = 1$ 是惯例。利用 $\rho_i(a) = \rho_{i,0} a^{-3(1+w_i)}$ 和 $H^2 = \frac{8\pi G}{3} \sum \rho_i - \frac{k c^2}{a^2}$，可写出哈勃参数作为红移的函数： $$ H(z)^2 = H_0^2 \left[ \Omega_{r,0} (1+z)^4 + \Omega_{m,0} (1+z)^3 + \Omega_{k,0} (1+z)^2 + \Omega_{\Lambda,0} \right] \tag{4} $$ 其中 $\Omega_{k,0} = 1 - \Omega_{r,0} - \Omega_{m,0} - \Omega_{\Lambda,0}$。这是宇宙学中最常用的公式之一，它将观测到的 $H(z)$ 与宇宙组分联系起来。

距离计算：共动距离、角直径距离和光度距离

观测宇宙学中，距离的定义依赖于观测方式。共动距离（comoving distance）沿测地线积分： $$ \chi(z) = c \int_0^z \frac{dz'}{H(z')} \tag{5} $$ 对于平坦宇宙（$\Omega_k=0$），角直径距离（angular diameter distance）和光度距离（luminosity distance）为： $$ D_A(z) = \frac{\chi(z)}{1+z}, \quad D_L(z) = (1+z) \chi(z) = D_A(z) (1+z)^2 $$ 这些距离公式是解释超新星、BAO等观测的基础。

宇宙微波背景辐射：早期宇宙的精确快照

CMB的基本性质与发现

宇宙微波背景辐射是大爆炸后约38万年（再复合时期）遗留下来的光子，其温度高度均匀，约为 $T_0 = 2.72548 \pm 0.00057$ K（COBE FIRAS测量），光谱为近乎完美的黑体谱。这强有力地支持了大爆炸理论。然而，温度各向异性（anisotropy）仅约 $\Delta T/T \sim 10^{-5}$，包含了宇宙初始条件、组成和几何的丰富信息。

各向异性的起源：暴胀与声学振荡

在暴胀（inflation）理论框架下，原初量子涨落在暴胀期间被放大到宇宙尺度，成为CMB各向异性的种子。这些初始曲率扰动具有近尺度不变的功率谱 $P_\mathcal{R}(k) \propto k^{n_s-1}$，其中 $n_s \approx 1$。

在光子-重子流体中，引力（由暗物质势阱驱动）与光子压力产生振荡，形成声学峰（acoustic peaks）。在再复合时刻（$z_* \approx 1100$），光子退耦，振荡被“冻结”，其相位信息编码在CMB温度涨落中。声学峰的物理尺度是声视界（sound horizon）： $$ r_s(z_*) = \int_{z_*}^{\infty} \frac{c_s(z)}{H(z)} dz \tag{6} $$ 其中 $c_s = c / \sqrt{3(1+R)}$ 是声速，$R = 3\rho_b/(4\rho_\gamma)$ 是重子-光子比。$r_s(z_*)$ 是早期宇宙中声波在光子退耦前能传播的最大距离，约150 Mpc（共动单位），是宇宙学中的“标准尺”。

CMB功率谱：从天空到 $C_\ell$

CMB各向异性在天球上展开为球谐函数： $$ \frac{\Delta T}{T}(\theta, \phi) = \sum_{\ell=0}^{\infty} \sum_{m=-\ell}^{\ell} a_{\ell m} Y_{\ell m}(\theta, \phi) $$ 假设各向同性，功率仅依赖于 $\ell$： $$ C_\ell \equiv \langle |a_{\ell m}|^2 \rangle $$ $C_\ell$ 是角功率谱，其形状包含丰富信息：

• **低 $\ell$（大角尺度）**： Sachs-Wolfe 效应，反映原初扰动和引力势。
• **第一个峰（$\ell \sim 200$）**：对应声视界角尺度 $\theta_s = r_s(z_*) / D_A(z_*)$。峰位 $\ell_1 \approx \pi / \theta_s$ 对空间曲率敏感：平坦宇宙 $\ell_1 \approx 220$，正曲率（封闭）峰移向更高 $\ell$，负曲率（开放）移向更低 $\ell$。
• **峰高与阻尼**：第一个峰高对重子密度 $\Omega_b h^2$ 敏感（重子增加惯性，提高峰高）；高 $\ell$ 的阻尼尾由光子扩散（Silk damping）造成，对 $\Omega_b h^2$ 和 $n_s$ 敏感。
• **相对峰高**：第二个峰相对于第一个的振幅对冷暗物质密度 $\Omega_c h^2$ 敏感（CDM影响势阱深度）。

WMAP与Planck：从探测到精确宇宙学

WMAP（2001-2010）首次提供了全天的CMB高精度地图，将宇宙学参数从“数量级”提升到“10%精度”。其关键结果支持平坦宇宙、$\Omega_\Lambda \approx 0.7$、$\Omega_m \approx 0.3$。

Planck（2009-2013，2015最终数据，2018释放）将精度推向极致。其功率谱测量到 $\ell \sim 2500$，系统误差极小。Planck 2018结果（基于高-ℓ TT, TE, EE + low-ℓ + lensing）给出（68% CL）： $$ \begin{aligned} \Omega_b h^2 &= 0.02237 \pm 0.00015 \\ \Omega_c h^2 &= 0.1200 \pm 0.0012 \\ \Omega_\Lambda &= 0.6847 \pm 0.0073 \\ H_0 &= 67.4 \pm 0.5 \ \text{km/s/Mpc} \\ n_s &= 0.9649 \pm 0.0042 \\ A_s &= (2.10 \pm 0.03) \times 10^{-9} \quad (\text{at } k=0.05 \ \text{Mpc}^{-1}) \\ \tau &= 0.0544 \pm 0.0073 \end{aligned} $$ 这些参数定义了所谓的“Planck基准模型”，成为现代宇宙学的黄金标准。注意 $H_0$ 与局部测量（如SH0ES的 $73.2 \pm 1.3$ km/s/Mpc）存在显著张力，这是当前研究热点。

参数拟合的统计学框架

CMB参数估计是典型的贝叶斯推断问题。给定数据 $d$（功率谱 $C_\ell^{\text{obs}}$ 及其协方差），模型参数 $\theta = (\Omega_b h^2, \Omega_c h^2, \dots)$ 的后验概率为： $$ P(\theta|d) \propto \mathcal{L}(d|\theta) \, P(\theta) $$ 其中似然函数 $\mathcal{L}$ 通常假设高斯分布（因多普勒峰测量精确）： $$ -2 \ln \mathcal{L} = (C_\ell^{\text{th}}(\theta) - C_\ell^{\text{obs}})^T \, \mathbf{Cov}^{-1} \, (C_\ell^{\text{th}}(\theta) - C_\ell^{\text{obs}}) $$ $C_\ell^{\text{th}}(\theta)$ 通过Boltzmann代码（如CAMB或CLASS）计算，求解光子-重子-暗物质系统的线性扰动方程。先验 $P(\theta)$ 通常取 wide flat priors。

采样采用马尔可夫链蒙特卡洛（MCMC）方法，如emcee或CosmoMC。收敛诊断（如Gelman-Rubin统计量）确保链已充分探索参数空间。后验分布给出参数最佳值和不确定性。

以下是一个简化的Python代码框架，演示如何用emcee拟合一个极简的CMB模型（仅含 $\Omega_b h^2$, $\Omega_c h^2$, $n_s$, $A_s$，固定其他参数）：

```python import numpy as np import emcee from scipy.linalg import cho_solve, cho_factor

假设我们有Planck-like的模拟数据：ell, C_ell_obs, Cov (协方差矩阵)

这里用占位符，实际需从文件加载或生成

ell = np.arange(2, 2500) # 多极矩 C_ell_obs = np.random.randn(len(ell)) * 1e-10 # 模拟观测值，实际有精确值 Cov = np.eye(len(ell)) * (1e-10)**2 # 对角协方差，简化

简化的理论计算函数（实际需调用CAMB/CLASS）

def compute_C_ell(theta, ell): """ 简化的CMB功率谱模型。实际中应使用Boltzmann代码。 theta = (ombh2, omch2, ns, As) 返回C_ell数组（单位：μK²） """ ombh2, omch2, ns, As = theta

这里仅用幂律近似：C_ell ∝ As * ell^{ns-1} * exp(-ell/ell_damp)

完全忽略声学峰结构！仅用于演示MCMC流程。

ell_damp = 1000 # 阻尼尺度，简化 C_ell = As * (ell / 100.0)**(ns-1) * np.exp(-ell/ell_damp)

转换为μK²：实际CAMB输出为(l(l+1)C_ell/(2π)) in μK²，这里简化

return C_ell * 1e10

对数似然函数

def log_likelihood(theta): C_th = compute_C_ell(theta, ell) residual = C_ell_obs - C_th

使用Cholesky分解求解协方差逆

try: factor, low = cho_factor(Cov) log_like = -0.5 * (residual @ cho_solve((factor, low), residual)) log_like -= 0.5 * len(residual) * np.log(2 * np.pi) # 常数项可省略 except np.linalg.LinAlgError: return -np.inf return log