宇宙学基础与ΛCDM模型技术剖析

宇宙学作为探索宇宙整体结构、演化和基本组成的科学，其核心目标是通过精确的观测和理论模型，重构宇宙从极早期至今的完整历史。在现代宇宙学体系中，ΛCDM（Lambda Cold Dark Matter，含宇宙学常数的冷暗物质）模型已成为描述宇宙大尺度行为的标准范式。该模型基于广义相对论和宇宙学原理（即宇宙在大尺度上均匀且各向同性），将宇宙的能量组成归结为暗能量（以宇宙学常数Λ表征）、冷暗物质（Cold Dark Matter, CDM）和重子物质，并辅以微小的初始密度扰动。ΛCDM模型成功解释了宇宙微波背景辐射（CMB）的黑体谱和各向异性功率谱、大尺度结构的形成、Ia型超新星的加速膨胀现象以及轻元素的丰度等关键观测事实。然而，随着观测技术的飞速进步和数据的日益精确，该模型也面临着严峻的挑战，其中最突出的便是哈勃张力（Hubble Tension）——即早期宇宙（通过CMB推断）与晚期宇宙（通过局部距离阶梯测量）得到的哈勃常数H₀之间存在超过5σ的显著差异。这一张力可能暗示着标准模型之外的新物理，如早期暗能量、修改引力理论或额外的相对论性粒子成分。因此，深入剖析ΛCDM模型的数学基础、理解其核心方程的推导过程、掌握宇宙学参数的观测约束技术，并能够通过数值工具进行参数拟合，已成为现代天体物理学研究者的必备技能。本章将系统性地展开这些内容，首先从宇宙学原理和大爆炸理论的基本框架出发，严格推导弗里德曼方程，进而详细阐述宇宙学参数（如物质密度参数Ωₘ、暗能量密度参数Ω_Λ、哈勃常数H₀）的观测约束方法，包括普朗克卫星的CMB功率谱分析、Ia型超新星的测距技术以及重子声波振荡（BAO）的观测原理。随后，通过一个基于Python Astropy库的完整实践案例，展示如何整合多源数据拟合宇宙学参数。最后，将深入讨论ΛCDM模型当前面临的主要局限性，特别是哈勃张力的可能起源及其对宇宙学未来发展的深远影响。

宇宙学原理与大爆炸理论框架

现代宇宙学的理论基石是宇宙学原理，该原理假设宇宙在足够大的尺度上（通常大于100 Mpc）是均匀的（homogeneous）和各向同性的（isotropic）。均匀性意味着宇宙的物质分布在大尺度上没有 preferential 位置，各向同性则意味着从任意位置观测，宇宙在宏观上呈现相同的性质，没有 preferential 方向。这一原理并非显然成立，但得到了星系巡天（如SDSS、2dF）和CMB观测的强力支持：星系在百万秒差距尺度上形成“宇宙网”（纤维状结构和空洞），但平均密度趋于一致；CMB温度各向异性在扣除偶极和多极后，剩余涨落仅约10⁻⁵量级，体现了高度的各向同性。

在宇宙学原理的约束下，时空度规必须具有弗里德曼-勒梅特-罗伯逊-沃尔克（FLRW）形式。采用球坐标（r, θ, φ），度规可写为：

$$ 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) $$

其中：

• $c$ 为光速，
• $t$ 为宇宙时（共动坐标），
• $a(t)$ 为**尺度因子**（scale factor），描述宇宙随时间的膨胀，通常归一化使得当前时刻 $t_0$ 满足 $a(t_0) = 1$，
• $k$ 为**曲率常数**，取值 $k=+1$（封闭球面）、$k=0$（平坦欧氏）、$k=-1$（开放双曲面），反映了空间的全局几何。

FLRW度规的几何性质由尺度因子 $a(t)$ 和曲率 $k$ 完全决定。共动坐标（comoving coordinates）随宇宙膨胀而“固定”在物质上，实际物理距离由 $a(t)$ 缩放。例如，两个共动点之间的固有距离（proper distance）为 $d(t) = a(t) \Delta r$。宇宙的膨胀动力学由爱因斯坦场方程决定，结合物质能量-动量张量，即可导出描述 $a(t)$ 演化的核心方程。

大爆炸理论（Big Bang Theory）并非指宇宙从一个点爆炸开来，而是描述宇宙从一个极热、极密的状态开始膨胀和冷却的图景。其核心证据链包括： 1. 哈勃-勒梅特定律：1929年，埃德温·哈勃发现星系的红移 $z$ 与距离 $d$ 近似成正比，$z \approx H_0 d / c$（对小红移），表明宇宙正在膨胀。红移 $z$ 定义为 $1+z = a(t_0)/a(t_{\text{em}})$，即波长被宇宙膨胀拉长。 2. 宇宙微波背景辐射（CMB）：1965年，彭齐亚斯和威尔逊意外发现全天空各向同性的微波辐射，温度约2.7 K，黑体谱完美。这是宇宙在约38万年时（再复合时期）释放的光子，因宇宙膨胀而红移到微波波段，是早期宇宙的“快照”。 3. 轻元素丰度：大爆炸核合成（BBN）理论预测，在宇宙最初几分钟内，高温高密环境引发核反应，产生氢（~75%质量）、氦-4（~25%）、以及微量的氘、氦-3、锂-7。这些丰度与最古老恒星和星系中观测到的值高度吻合。

这三个支柱共同构成了大爆炸范式的坚实基础。然而，标准大爆炸模型本身无法解释宇宙的平坦性、视界问题以及初始扰动的来源，这些缺陷由暴胀理论（Inflation）弥补。暴胀提出宇宙在极早期（约10⁻³⁶秒）经历指数膨胀，解决了平坦性（使 $k \approx 0$）和视界问题（使可观测宇宙源自因果联系区域），并量子涨落放大为宇宙结构的种子。暴胀的细节虽未完全确定，但其预言（如标量谱指数 $n_s \approx 0.96$）已得到CMB观测的强有力支持。

弗里德曼方程的严格推导

弗里德曼方程是描述FLRW宇宙动力学的一对微分方程，直接源于爱因斯坦场方程。本节将展示从场方程到弗里德曼方程的完整推导过程，这是理解宇宙学模型数学结构的关键。

爱因斯坦场方程为：

$$ 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$ 是牛顿引力常数。

对于FLRW度规，度规张量分量为：

$$ g_{00} = -c^2, \quad g_{ij} = a(t)^2 \tilde{g}_{ij} $$

其中 $\tilde{g}_{ij}$ 是三维空间度规（曲率为 $k$ 的空间度规）。为简化，可引入共时坐标（proper time）$d\tau = dt$（对共动观测者），并令 $c=1$（自然单位），稍后恢复。

计算克里斯托费尔符号（Christoffel symbols）：

$$ \Gamma^\rho_{\mu\nu} = \frac{1}{2} g^{\rho\sigma} \left( \partial_\mu g_{\nu\sigma} + \partial_\nu g_{\mu\sigma} - \partial_\sigma g_{\mu\nu} \right) $$

非零分量包括：

• $\Gamma^0_{ij} = a \dot{a} \tilde{g}_{ij}$（因为 $g_{ij}=a^2\tilde{g}_{ij}$，对时间求导得 $\partial_0 g_{ij}=2a\dot{a}\tilde{g}_{ij}$）
• $\Gamma^i_{0j} = \Gamma^i_{j0} = \dot{a}/a \delta^i_j$
• $\Gamma^i_{jk}$ 由空间度规 $\tilde{g}_{ij}$ 决定，记为 $\tilde{\Gamma}^i_{jk}$。

接着计算里奇张量（Ricci tensor）$R_{\mu\nu} = \partial_\rho \Gamma^\rho_{\mu\nu} - \partial_\nu \Gamma^\rho_{\mu\rho} + \Gamma^\rho_{\rho\lambda} \Gamma^\lambda_{\mu\nu} - \Gamma^\rho_{\nu\lambda} \Gamma^\lambda_{\mu\rho}$。

• $R_{00}$ 分量：

$$ R_{00} = -\frac{3\ddot{a}}{a} $$ 推导：$\Gamma^\rho_{00}$ 非零仅当 $\rho=0$？实际上 $\Gamma^i_{00}=0$，$\Gamma^0_{00}=0$，但需计算 $\partial_i \Gamma^i_{00}$ 等。更系统：$R_{00} = \partial_\rho \Gamma^\rho_{00} - \partial_0 \Gamma^\rho_{\rho 0} + \cdots$。已知 $\Gamma^i_{0j} = (\dot{a}/a)\delta^i_j$，则 $\Gamma^\rho_{\rho 0} = \Gamma^0_{00} + \Gamma^i_{i0} = 0 + 3(\dot{a}/a)$。$\partial_0 \Gamma^\rho_{\rho 0} = 3(\ddot{a}/a - \dot{a}^2/a^2)$。$\partial_\rho \Gamma^\rho_{00}$：$\Gamma^i_{00}=0$，$\Gamma^0_{00}=0$，所以为0。交叉项：$\Gamma^\rho_{\rho\lambda}\Gamma^\lambda_{00} - \Gamma^\rho_{0\lambda}\Gamma^\lambda_{\rho 0}$。第一项：$\Gamma^\rho_{\rho\lambda}\Gamma^\lambda_{00}=0$（因 $\Gamma^\lambda_{00}=0$）。第二项：$\Gamma^\rho_{0\lambda}\Gamma^\lambda_{\rho 0}$。对 $\rho,\lambda$ 求和。$\Gamma^i_{0j}\Gamma^j_{i0} = (\dot{a}/a)^2 \delta^i_j \delta^j_i = 3(\dot{a}/a)^2$。其他组合为零。故 $R_{00} = 0 - 3(\ddot{a}/a - \dot{a}^2/a^2) - 3(\dot{a}/a)^2 = -3\ddot{a}/a$。

• $R_{ij}$ 分量：

$$ R_{ij} = \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) g_{ij} $$ 推导：利用空间曲率，$R_{ij} = \tilde{R}_{ij} - \frac{\ddot{a}}{a} g_{ij} - 2\frac{\dot{a}^2}{a^2} g_{ij}$？更直接：对于FLRW，$R_{ij} = \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) g_{ij}$，其中 $\tilde{R}_{ij} = 2k \tilde{g}_{ij}$ 是三维空间里奇张量（常曲率空间）。验证：$\tilde{R}_{ij} = (k) \tilde{g}_{ij}$ 对于三维球面/双曲面？实际上，对于度规 $d\ell^2 = dr^2/(1-kr^2) + r^2 d\Omega^2$，里奇标量 $\tilde{R}=6k$，且 $\tilde{R}_{ij} = k \tilde{g}_{ij}$？需查证：常曲率空间有 $\tilde{R}_{ij} = (n-1)k \tilde{g}_{ij}$，n=3时 $\tilde{R}_{ij}=2k \tilde{g}_{ij}$。是的，所以 $\tilde{R}_{ij} = 2k \tilde{g}_{ij} = (2k/a^2) g_{ij}$。然后 $R_{ij} = \tilde{R}_{ij} - \frac{\ddot{a}}{a} g_{ij} - 2\frac{\dot{a}^2}{a^2} g_{ij}$？从标准公式：$R_{0i0j} = -a\ddot{a} \tilde{g}_{ij}$，$R_{ijkl} = a^2(\dot{a}^2 + k) (\tilde{g}_{ik}\tilde{g}_{jl} - \tilde{g}_{il}\tilde{g}_{jk})$，然后收缩得 $R_{ij} = R^k_{ikj} = g^{kl}R_{kilj}$。计算得 $R_{ij} = \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) g_{ij}$。正确。

里奇标量：

$$ R = g^{\mu\nu} R_{\mu\nu} = g^{00}R_{00} + g^{ij}R_{ij} = (-\frac{1}{c^2})(-3\ddot{a}/a) + \frac{1}{a^2} \delta^{ij} \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) g_{ij} $$

注意 $g^{ij} = (1/a^2) \tilde{g}^{ij}$，且 $\tilde{g}^{ij}\tilde{g}_{ij}=3$，所以 $g^{ij}g_{ij}=3$。因此：

$$ R = \frac{3\ddot{a}}{c^2 a} + \frac{3}{a^2} \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) = \frac{3}{c^2}\frac{\ddot{a}}{a} + 3\frac{\ddot{a}}{a^3} + 6\frac{\dot{a}^2}{a^4} + \frac{6k}{a^4} $$

这看起来混乱，因为我们混合了 $c$。最好在推导中保持 $c$ 显式，或使用 $c=1$ 然后恢复。通常宇宙学中设 $c=1$，但为清晰，我们保留 $c$。

实际上，在FLRW度规中，$g_{00}=-c^2$，所以 $g^{00}=-1/c^2

$$ R = g^{\mu\nu} R_{\mu\nu} = g^{00} R_{00} + g^{ij} R_{ij} $$

在 $c=1$ 单位制下，$g^{00} = -1$，$g^{ij} = a^{-2} \tilde{g}^{ij}$，且 $\tilde{g}^{ij} \tilde{g}_{ij} = 3$。代入 $R_{00} = -3\ddot{a}/a$ 和 $R_{ij} = \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) g_{ij}$，得：

$$ R = (-1) \cdot \left(-\frac{3\ddot{a}}{a}\right) + a^{-2} \tilde{g}^{ij} \cdot \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) a^2 \tilde{g}_{ij} = \frac{3\ddot{a}}{a} + 3 \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) = 6\left( \frac{\ddot{a}}{a} + \frac{\dot{a}^2}{a^2} + \frac{k}{a^2} \right) $$

由此可得爱因斯坦张量分量：

• $G_{00} = R_{00} - \frac{1}{2} R g_{00} = -\frac{3\ddot{a}}{a} - \frac{1}{2} \cdot 6\left( \frac{\ddot{a}}{a} + \frac{\dot{a}^2}{a^2} + \frac{k}{a^2} \right) \cdot (-1) = -\frac{3\ddot{a}}{a} + 3\left( \frac{\ddot{a}}{a} + \frac{\dot{a}^2}{a^2} + \frac{k}{a^2} \right) = 3\left( \frac{\dot{a}^2}{a^2} + \frac{k}{a^2} \right)$
• $G_{ij} = R_{ij} - \frac{1}{2} R g_{ij} = \left( \frac{\ddot{a}}{a} + 2\frac{\dot{a}^2}{a^2} + \frac{2k}{a^2} \right) g_{ij} - \frac{1}{2} \cdot 6\left( \frac{\ddot{a}}{a} + \frac{\dot{a}^2}{a^2} + \frac{k}{a^2} \right) g_{ij} = \left( -2\frac{\ddot{a}}{a} - \frac{\dot{a}^2}{a^2} - \frac{k}{a^2} \right) g_{ij}$

对于完美流体，能量-动量张量为 $T_{\mu\nu} = (\rho + p) u_\mu u_\nu + p g_{\mu\nu}$，其中四速度 $u^\mu = (1,0,0,0)$（共动坐标），故 $u_\mu = g_{\mu\nu} u^\nu$ 给出 $u_0 = -1$，$u_i=0$。因此：

• $T_{00} = (\rho + p)(-1)(-1) + p g_{00} = \rho + p - p = \rho$（因为 $g_{00}=-1$）
• $T_{ij} = p g_{ij}$

将 $G_{\mu\nu}$ 与 $T_{\mu\nu}$ 代入场方程 $G_{\mu\nu} + \Lambda g_{\mu\nu} = 8\pi G T_{\mu\nu}$（$c=1$）：

00 分量： $$ 3\left( \frac{\dot{a}^2}{a^2} + \frac{k}{a^2} \right) + \Lambda (-1) = 8\pi G \rho \quad\Rightarrow\quad \frac{\dot{a}^2}{a^2}