宇宙的组成：可见与不可见

当我们仰望星空，看到的恒星、星系和星云只是宇宙的冰山一角。根据欧洲空间局普朗克卫星2018年发布的最新数据，宇宙的总能量密度组成中，我们熟悉的普通重子物质（Baryonic Matter）仅占约5%，而神秘的暗物质（Dark Matter）约占27%，暗能量（Dark Energy）则占据了剩余的68%。这意味着，宇宙中超过95%的成分是我们无法直接观测到的“暗”宇宙。这一发现彻底改变了我们对宇宙本质的理解，也开启了现代宇宙学最激动人心的探索篇章。本章将基于精确的观测数据，建立宇宙物质-能量组成的框架，并通过对比不同尺度结构的质量-光度比，直观展示“缺失质量”问题，最终引入关键概念：宇宙临界密度、哈勃常数争议以及宇宙微波背景（Cosmic Microwave Background, CMB）各向异性功率谱。我们还将用Python代码示例演示如何从WMAP/Planck数据反推宇宙组分，让抽象的理论变得触手可及。

宇宙的精确成分：普朗克2018年的答案

现代宇宙学进入精确测量时代，得益于卫星任务的精确观测。其中，欧洲空间局的普朗克卫星（Planck）于2009年发射，专门用于绘制宇宙微波背景（CMB）的精细图样。CMB是大爆炸后约38万年时，宇宙首次透明时释放的光子，如今红移至微波波段，均匀充满整个天空，温度约为2.725 K。这些光子携带了早期宇宙的“化石信息”，其微小的温度涨落（各向异性）反映了原初的密度扰动，进而决定了今天宇宙的大尺度结构。通过测量CMB温度涨落的角功率谱， cosmologists可以精确推断宇宙的几何、物质含量和初始条件。

普朗克卫星在2013年、2015年和2018年陆续发布数据，精度不断提升。2018年发布的最终结果（Planck 2018 results. VI. Cosmological parameters）基于完整的卫星任务数据，结合其他观测（如重子声波振荡Baryon Acoustic Oscillations, BAO），给出了宇宙学参数的黄金标准。在标准ΛCDM（Lambda Cold Dark Matter）模型下，宇宙的总能量密度由以下组分构成：

• **普通重子物质（Baryons）**：构成原子、恒星、行星和我们的物质，约占宇宙总能量密度的 **4.9%**（具体为Ω_b h² = 0.02237 ± 0.00015，其中h = H₀/100 km/s/Mpc）。
• **冷暗物质（Cold Dark Matter, CDM）**：不参与电磁相互作用，仅通过引力影响可见物质，约占 **26.8%**（Ω_c h² = 0.1200 ± 0.0012）。
• **暗能量（Dark Energy）**：表现为宇宙加速膨胀的驱动力，通常用宇宙学常数Λ描述，约占 **68.3%**（Ω_Λ = 0.6847 ± 0.0073）。
• **辐射（Photons and Neutrinos）**：包括CMB光子和中微子，占比极小，约0.01%，在当前宇宙演化中可忽略。

这些比例基于宇宙的总密度参数Ω_total = Ω_b + Ω_c + Ω_Λ ≈ 1（更精确为1.0007 ± 0.0019），表明宇宙在几何上是平坦的（Flat），符合暴胀理论的预言。密度参数定义为组分密度与临界密度（Critical Density）之比。临界密度是宇宙恰好平坦所需的平均密度，由哈勃常数H₀决定：

$$ \rho_c = \frac{3H_0^2}{8\pi G} $$

其中G为引力常数。当前H₀的最佳估计值存在争议：普朗克CMB数据给出H₀ = 67.36 ± 0.54 km/s/Mpc（假设ΛCDM模型），而局部宇宙测量（如造父变星和Ia型超新星）得到H₀ ≈ 73 km/s/Mpc。这一“哈勃张力”（Hubble Tension）可能暗示新物理，如早期暗能量或修改引力理论，是当前研究前沿。

从CMB功率谱到宇宙组分：物理图像

CMB的温度涨落通常用角功率谱C_ℓ描述，它表示不同角尺度（多极矩ℓ）上的涨落方差。功率谱呈现一系列峰值，其位置和高度对宇宙参数极为敏感：

• **第一峰（ℓ ≈ 220）**：对应约1°角尺度，由声波在光子-重子流体中的振荡形成。其位置主要敏感于**宇宙几何**：平坦宇宙中，第一峰出现在ℓ ≈ 220；若宇宙闭合（正曲率），峰向左移；开放（负曲率）则向右移。普朗克数据精确命中ℓ ≈ 220，强有力支持平坦宇宙。
• **第二峰和第三峰**：高度由**重子密度**Ω_b h²和**暗物质密度**Ω_c h²控制。重子增加流体惯性，抑制后续振荡，使第二峰相对第一峰降低；暗物质不参与声波，但提供额外引力势，影响振荡幅度。通过拟合峰值高度比，可以精确分离Ω_b和Ω_c。
• **阻尼尾部（ℓ > 1000）**：小尺度涨落因光子扩散（ Silk Damping）而衰减，其形状提供额外约束。

图1.1（概念图）展示了普朗克2018的TT（温度-温度）功率谱，误差棒极小，理论与观测完美吻合。这种精确匹配是ΛCDM模型的巨大成功，但也留下谜题：为何暗物质和暗能量的比例如此巧合？这引出人择原理等哲学讨论。

密度参数的计算示例

基于普朗克2018的最佳拟合参数（base_plikHM_TTTEEE_lowl_lowE_lensing），我们可以计算各组分的密度参数Ω（相对于临界密度）。注意，Planck直接给出Ω_b h²和Ω_c h²，需除以h²得到Ω_b和Ω_c，其中h = H₀/100。使用H₀ = 67.36 km/s/Mpc，则h = 0.6736。

计算过程：

• Ω_b = (Ω_b h²) / h² = 0.02237 / (0.6736²) ≈ 0.0493（4.93%）
• Ω_c = (Ω_c h²) / h² = 0.1200 / (0.6736²) ≈ 0.264（26.4%）
• 假设平坦宇宙Ω_total = 1，则Ω_Λ = 1 - Ω_b - Ω_c ≈ 0.686（68.6%）

与官方给出的0.6847非常接近，微小差异源于参数后验分布的细节。这些数字已成为现代宇宙学的“常识”，但背后是数十年观测技术的积累。

可见与不可见：质量-光度比的阶梯式揭示

宇宙的“暗”并非凭空猜测，而是从多个独立观测尺度上，通过质量-光度比（Mass-to-Light Ratio, M/L）的异常逐步揭示的。M/L表示单位光度对应的质量，太阳的M/L约为1（以太阳单位）。如果宇宙仅由可见恒星和气体组成，M/L应接近常数。然而，从太阳系到宇宙网，M/L急剧上升，暗示大量不可见质量的存在。

太阳系：M/L接近恒星值

在太阳系内，行星、小行星等天器的质量可通过轨道动力学精确测量（如利用探测器跟踪数据）。太阳系的总光度主要来自太阳，M/L ≈ 1-2 M_⊙/L_⊙，与单颗恒星一致。没有证据表明太阳系内存在大量暗物质，因为暗物质如果存在，应形成弥漫的晕，但其密度极低（本地暗物质密度约0.3 GeV/cm³），对太阳系动力学影响远小于行星摄动，目前无法直接探测。这提示暗物质分布是宇宙大尺度的，在星系尺度上才显著。

银河系：旋转曲线揭示暗物质晕

20世纪30年代，弗里茨·兹威基（Fritz Zwicky）在后发座星系团中首次提出“暗物质”概念，但直接证据来自星系尺度。1970年代，薇拉·鲁宾（Vera Rubin）和肯特·福特（Kent Ford）利用光谱仪测量螺旋星系旋转曲线（Rotation Curve）。根据牛顿力学，如果星系质量集中在可见盘内，旋转速度v应随半径r增加而下降（v ∝ 1/√r）。但观测显示，在盘外区域（r > 10 kpc），旋转速度几乎恒定（图1.2概念图），形成“平坦旋转曲线”。这意味着质量分布远超出可见光范围，必须存在一个巨大的、不可见的暗物质晕（Dark Matter Halo），其质量可达可见盘的10倍以上。

数学上，平坦旋转曲线v ≈ constant意味着 enclosed mass M(r) ∝ r，因此密度ρ ∝ 1/r²，形成等温晕 profile。典型螺旋星系的M/L可达10-100 M_⊙/L_⊙，远高于恒星值。例如，银河系的 rotation curve 在太阳位置（r ≈ 8 kpc）v ≈ 220 km/s，推算的总质量（包括晕）约10¹² M_⊙，而总光度约2×10¹⁰ L_⊙，M/L ≈ 50。这一异常无法用恒星质量-光度比变化解释，因为恒星M/L一般不超过10。

星系团：多重证据汇聚

星系团是宇宙中最大的引力束缚结构，包含数百至数千星系、热气体（X射线源）和暗物质。这里，暗物质证据来自三种独立方法，且结果一致：

1. 动力学质量：通过星系速度弥散σ估算维里质量。兹威基1933年测量后发座星系团，发现从红移得到的速度弥散对应的维里质量比从光度估计的可见质量大400倍，首次提出“dunkle Materie”（暗物质）。 2. X射线热气体：星系团内介质（Intracluster Medium, ICM）温度达10⁷-10⁸ K，发出X射线。根据流体静力平衡，气体压力梯度与引力平衡，可推算总质量分布。X射线质量通常与动力学质量一致，且远大于可见星系质量。 3. 引力透镜：大质量星系团扭曲背景星系的光，形成弧形或多重像。通过透镜模型可直接映射质量分布，无需假设动力学状态。著名案例子弹星系团（1E0657-558） 中，X射线气体（普通物质）因碰撞而滞后，而引力透镜显示的质量中心与星系分布重合，直接证明大部分质量不与可见物质耦合，是暗物质的“ smoking gun”。

图1.3展示了子弹星系团的多波段合成图：X射线（粉红）显示热气体，光学（红蓝）显示星系，引力透镜质量 contours（蓝绿）显示质量峰值与星系一致，而非气体。这成为暗物质存在的最直观证据。

宇宙大尺度结构：数值模拟的必然要求

在更大尺度上，宇宙呈现“宇宙网”结构：星系聚集在丝状结构和节点，中间是巨大空洞。这种结构的形成需要用冷暗物质（Cold Dark Matter, CDM）模型解释。在ΛCDM框架下，暗物质在早期宇宙由于微小密度扰动而开始坍缩，形成暗物质晕的“骨架”，普通重子物质随后落入势阱，形成星系和星系团。纯重子模型无法在宇宙年龄内形成观测到的结构，因为辐射压力会阻止坍缩。著名的Millennium Simulation（2005）使用CDM模型，成功再现了宇宙大尺度结构的统计性质（如两点相关函数、功率谱），而如果去掉暗物质，模拟结果与观测格格不入。

因此，从星系旋转曲线到宇宙网，M/L的阶梯式上升——太阳系~1，银河系~50，星系团~300，宇宙平均~300——构成了暗物质存在的多层次证据链。暗物质不是假设，而是当前宇宙学标准模型的基石。

暗能量的意外发现：宇宙加速膨胀

如果说暗物质的发现是“缺失质量”的渐进揭示，暗能量的出现则完全出乎意料。1998年，两个独立团队——超新星宇宙学计划（Supernova Cosmology Project, SCP）和高红移超新星搜索队（High-z Supernova Search Team）——通过观测遥远Ia型超新星，发现宇宙膨胀在加速，从而推断存在一种排斥性能量成分，即暗能量。

Ia型超新星：标准烛光

Ia型超新星源于白矮星吸积伴星物质达到钱德拉塞卡极限（约1.4 M_⊙）时的热核爆炸。其峰值光度高度一致（绝对星等M_B ≈ -19.3），因此可作为“标准烛光”测量宇宙距离。通过测量超新星的红移z和视亮度，可以绘制哈勃图（距离-红移关系）。在物质主导的减速膨胀宇宙中，遥远超新星应比预期更亮（因为膨胀在過去更快，距离更近）；而在加速膨胀宇宙中，它们会比预期更暗（因为膨胀在過去更慢，距离更远）。

1998年，两个团队分别观测到高红移（z ≈ 0.5-1）超新星比减速模型预测暗约0.25星等，统计显著性超过5σ。这一结果震惊学界，因为当时普遍预期宇宙在物质引力作用下应减速膨胀。加速膨胀意味着存在负压成分，根据广义相对论，这需要一种具有负方程态参数w = p/ρ < -1/3的能量。最简单的解释是宇宙学常数Λ，对应w = -1，其能量密度

暗能量的本质：从宇宙学常数到精质场

宇宙加速膨胀的发现直接指向一种具有负压的神秘能量成分，其方程态参数 $w = p/\rho$ 必须小于 $-1/3$ 才能产生排斥引力。最简单的候选者是宇宙学常数（Cosmological Constant），记为 $\Lambda$，对应 $w = -1$，其能量密度 $\rho_\Lambda$ 在时空均匀分布且不随宇宙膨胀而稀释（$\rho_\Lambda \propto a^0$）。在爱因斯坦场方程中，$\Lambda$ 作为几何项加入：

$$R_{\mu\nu} - \frac{1}{2}R g_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$$

这等价于在能量-动量张量中添加一项 $T_{\mu\nu}^{(\Lambda)} = -\frac{\Lambda c^4}{8\pi G} g_{\mu\nu}$，即真空能量具有负压 $p_\Lambda = -\rho_\Lambda c^2$。

然而，将 $\Lambda$ 解释为真空零点能时，量子场论预测的真空能量密度比观测值大 $10^{120}$ 倍，这是物理学中最严重的理论不匹配，称为宇宙学常数问题（Cosmological Constant Problem）。为解决此问题，物理学家提出精质场（Quintessence）模型：一种动态的、极轻的标量场 $\phi$，其势能 $V(\phi)$ 缓慢变化，产生有效的 $w \approx -1$ 但可能随时间演化。更激进的模型包括幻能量（Phantom Energy，$w < -1$）和修改引力理论（如 $f(R)$ 引力），但当前观测（如Planck 2018结合BAO和SNe）支持 $w = -1.03 \pm 0.03$，与宇宙学常数一致，未发现显著演化。

暗能量的均匀性也带来结构形成难题：如果暗能量过早主导，会抑制星系和星系团的形成。观测显示，暗能量在红移 $z \approx 0.7$ 后才开始主导，此前物质密度更高，结构得以生长。这要求 $\Omega_\Lambda$ 与 $\Omega_m$ 的比例恰好允许我们存在——这又引出了人择原理的讨论：如果 $\Omega_\Lambda$ 过大，宇宙会过早加速，星系无法形成；过小则宇宙可能已坍缩。我们的观测值处于允许星系形成的狭窄窗口内。

哈勃常数争议：局部宇宙与早期宇宙的裂痕

哈勃常数 $H_0$ 描述当前宇宙膨胀速率，单位 km/s/Mpc。其值直接影响宇宙年龄（$t_0 \approx 1/H_0$ 对于平坦宇宙）和临界密度 $\rho_c \propto H_0^2$。精确测定 $H_0$ 是宇宙学的核心任务之一，但两种主流方法给出了显著不同的结果，形成了哈勃张力（Hubble Tension）。

方法一：早期宇宙推断（CMB各向异性） 基于普朗克卫星对CMB温度涨落的测量，在假设ΛCDM模型正确的前提下，通过拟合功率谱得到 $H_0$。CMB在红移 $z \approx 1100$ 时释放，其涨落模式受早期宇宙物理（如声学尺度 $\theta_*$）约束。声学尺度由重子-光子流体的声速和扩散阻尼决定，观测值 $\theta_*^{\text{obs}}$ 与理论预测 $\theta_*^{\text{th}}$ 匹配，但 $\theta_*^{\text{th}}$ 依赖于 $\Omega_m h^2$ 和 $\Omega_b h^2$，而 $h = H_0/100$。因此，给定 $\Omega_m$ 和 $\Omega_b$，$\theta_*$ 固定了 $H_0$。普朗克2018给出 $H_0 = 67.36 \pm 0.54$ km/s/Mpc（TT,TE,EE+lowE+lensing），宇宙年龄约137.9亿年。

方法二：局部距离阶梯（Local Distance Ladder） 通过测量近距天体的距离和红移直接确定 $H_0$。经典方法：造父变星（Cepheids）作为第一级标准烛光，其周期-光度关系校准；再以Ia型超新星作为第二级，延伸至更远距离。哈勃空间望远镜（HST）的SH0ES团队（Supernovae and H0 for the Equation of State）精化此链，得到 $H_0 = 73.04 \pm 1.04$ km/s/Mpc（2022年更新）。另一独立方法使用引力透镜时间延迟：类星体被前景星系团透镜成多像，光路不同导致时间延迟，结合质量模型可直接给出 $H_0$。H0LiCOW合作组结果 $H_0 = 73.3^{+1.7}_{-1.8}$ km/s/Mpc，与SH0ES一致。

张力统计：两者相差约 $5\sigma$，远超随机误差。如果这不是系统误差，可能暗示ΛCDM模型需要扩展。 proposed solutions 包括：

• **早期暗能量**（Early Dark Energy, EDE）：在CMB释放前短暂贡献约10%能量密度，改变声学尺度，从而在保持CMB拟合的同时提高 $H_0$。
• **额外 relativistic 物种**（如额外中微子或暗辐射）：增加早期辐射密度，影响膨胀历史。
• **修改引力**：改变早期宇宙的扰动增长。
• **系统误差**：如造父变星金属丰度校准、星系团质量模型偏差。但多方独立验证使系统误差解释难度增大。

这一张力是当前宇宙学最活跃的研究领域，可能指向新物理。

宇宙微波背景各向异性功率谱：宇宙的指纹

CMB温度各向异性 $\Delta T(\hat{n})/T$ 在天空球面上的展开为：

$$\frac{\Delta T}{T}(\hat{n}) = \sum_{\ell=0}^{\infty}\sum_{m=-\ell}^{\ell} a_{\ell m} Y_{\ell m}(\hat{n})$$

其中 $Y_{\ell m}$ 是球谐函数，$a_{\ell m}$ 是系数。各向异性功率谱定义为：

$$C_\ell = \langle |a_{\ell m}|^2 \rangle$$

（对于各向同性天空，$C_\ell$ 与 $m$ 无关）。$C_\ell$ 表征角尺度 $\theta \approx 180^\circ/\ell$ 上的涨落方差。图1.4（概念图）展示普朗克2018的TT功率谱，从低ℓ（大尺度）到高ℓ（小尺度）呈现复杂