第1章：宇宙的组成与未知 / Chapter 1: The Composition of the Universe and the Unknown

宇宙的绝大部分是隐形的。当我们用望远镜凝视星空，看到的恒星、星系、星云——所有这些发光物质——加起来仅占宇宙总质能的5%。剩下的95%由两种神秘成分主导：暗物质（Dark Matter）约占27%，暗宇宙学常数（Dark Energy）约占68%。这一惊人结论并非科幻猜想，而是基于过去三十年精密宇宙学观测的坚实推论，其核心支柱是宇宙微波背景辐射（Cosmic Microwave Background, CMB）的测量，特别是欧洲空间局普朗克（Planck）卫星2018年发布的数据。这些数字构成了现代宇宙学的标准模型——ΛCDM（Lambda Cold Dark Matter），它像一张宇宙的“成分清单”，却揭示了我们对宇宙本质的认知存在巨大空白。本章将深入剖析这张清单如何被绘制出来，从CMB的微小温度涨落到宇宙的大尺度结构，展示如何将天空中的光点转化为精确的宇宙学参数，并回顾人类如何从赫歇尔（William Herschel）的星云图一步步走到今天的宇宙学精确时代。

1.1 宇宙学标准模型ΛCDM：一个简洁而成功的框架

ΛCDM模型是当前描述宇宙演化的最简有效理论。它假设宇宙在大尺度上均匀各向同性，物质成分包括普通重子物质（Baryons）、冷暗物质（Cold Dark Matter）和暗能量（表现为宇宙学常数Λ），空间几何由爱因斯坦场方程决定。模型的六个基本参数通常取为：重子物理密度Ω_b h²、冷暗物质物理密度Ω_c h²、再复合时期的声学角度θ_*、再电离光学深度τ、标量谱指数n_s以及标量振幅A_s。其中h = H₀/100 km/s/Mpc是归一化哈勃常数。通过普朗克2018年数据，这些参数被精确测定，进而推导出今天的组分比例：

普通物质（重子）密度：Ω_b = Ω_b h² / h² ≈ 0.049 冷暗物质密度：Ω_c = Ω_c h² / h² ≈ 0.265 暗能量密度（Λ）：Ω_Λ ≈ 0.686 曲率密度：Ω_k = 1 - (Ω_b + Ω_c + Ω_Λ) ≈ 0（宇宙平坦）

这些比例之和严格满足宇宙几何闭合条件：Ω_m + Ω_Λ + Ω_k = 1，其中Ω_m = Ω_b + Ω_c ≈ 0.314。普朗克2018年给出的最佳拟合值（基于plik_lite组合）为：

$$ Ω_b h^2 = 0.02237 \pm 0.00015, \quad Ω_c h^2 = 0.1200 \pm 0.0012, \quad 100θ_* = 1.04110 \pm 0.00030, \quad τ = 0.0544 \pm 0.0073, \quad \ln(10^{10} A_s) = 3.044 \pm 0.014, \quad n_s = 0.9649 \pm 0.0042 $$

结合H₀ = 67.4 ± 0.5 km/s/Mpc（普朗克TT,TE,EE+lowE+lensing），计算得Ω_b ≈ 0.0486，Ω_c ≈ 0.2589，Ω_Λ ≈ 0.6925。注意不同数据组合（如加入BAO）会略有调整，但核心结论稳定：普通物质仅约5%，暗物质约27%，暗能量约68%。这张“成分表”的权威性来自多个独立观测的交叉验证，而CMB是最关键的“定锚点”。

1.2 宇宙微波背景辐射：大爆炸的余晖与信息宝库

CMB是大爆炸后约38万年时释放的光子。此时宇宙温度降至约3000K，电子与质子结合成中性氢原子（再复合时期），光子得以自由传播，不再频繁散射。这些光子至今仍充满宇宙，但由于宇宙膨胀约1100倍，其波长被拉伸到微波波段，黑体谱峰值对应于约2.725K的温度。1965年，阿诺·彭齐亚斯（Arno Penzias）和罗伯特·威尔逊（Robert Wilson）意外发现这种弥漫各向的微波噪声，为宇宙学大爆炸理论提供了首个直接证据。1992年，COBE卫星首次精确测量CMB黑体谱并发现各向异性，获得诺贝尔奖。随后WMAP（2001-2010）和普朗克（2009-2013）将测量精度推向极致。

CMB的温度涨落ΔT/T极小，仅约10^{-5}量级。这些涨落源于早期宇宙的量子涨落在暴胀时期被放大，成为后来星系和大尺度结构的“种子”。涨落在天空中的分布是各向异性的，但统计上具有高斯性和各向同性，因此完全由角功率谱C_l描述。定义球谐系数a_{lm}为：

$$ \frac{ΔT}{T}(\theta, \phi) = \sum_{l=0}^{\infty} \sum_{m=-l}^{l} a_{lm} Y_{lm}(\theta, \phi) $$

则角功率谱为：

$$ C_l = \frac{1}{2l+1} \sum_{m=-l}^{l} |a_{lm}|^2 $$

C_l随l的变化曲线（通常绘图以l(l+1)C_l/(2π)对l）呈现一系列峰值和谷值，称为声学峰。这些峰值的物理起源是光子-重子流体在引力势阱中的声波振荡。在再复合前后，流体被引力拉入势阱，辐射压力又将其推出，形成振荡。不同l对应不同角度尺度：l ≈ 180°/θ，因此低l（大尺度）对应超视界涨落，高l（小尺度）对应小角度涨落。第一个峰值（l≈220）对应声波在再复合时期完成半个振荡的最大尺度，其位置对宇宙几何极度敏感：如果宇宙平坦，该峰值出现在l≈220；如果空间正曲率（闭合），峰值左移；负曲率（开放）则右移。普朗克数据精确显示峰值在l≈220，证实Ω_k ≈ 0，宇宙平坦。

峰值的高度和相对幅度则携带物质成分信息。重子物质（Ω_b）增加流体惯性，增强声波振幅，使第一个峰值相对第二个更高。暗物质（Ω_c）不参与电磁耦合，只提供额外引力势阱，影响振荡的初始条件和演化，从而改变峰值序列。通过比较观测C_l与理论计算（基于Boltzmann方程求解光子-重子-暗物质流体动力学），可以同时约束Ω_b、Ω_c、Ω_Λ等参数。理论计算需用代码如CAMB或CLASS，输入基本参数，输出理论功率谱。然后通过贝叶斯马尔可夫链蒙特卡洛（MCMC）方法（如CosmoMC或MontePython）扫描参数空间，找到使似然函数最大的参数组合。

1.3 多信使验证：大尺度结构与超新星的独立证据

CMB提供的是宇宙38万年时的“快照”，而宇宙组分比例需要与今天观测到的结构形成一致。大尺度结构（Large-Scale Structure, LSS）巡天，如斯隆数字化巡天（SDSS）、BOSS（Baryon Oscillation Spectroscopic Survey）和eBOSS，测量了数千万星系的空间分布。这些分布也存在声学特征，称为重子声学振荡（Baryon Acoustic Oscillation, BAO）。BAO是CMB声学峰在物质功率谱中的遗迹，其特征尺度（约150 Mpc）是标准尺，可用于测量宇宙膨胀历史。LSS功率谱P(k)的形状对Ω_m和Ω_b敏感：小尺度（高k）功率受重子拖曳抑制，大尺度（低k）形状反映物质密度。将CMB预测的初始功率谱（通过转移函数演化到今天）与LSS观测比较，可以独立检验Ω_m和Ω_b。结果与CMB高度一致：例如，BAO测量给出的Ω_m ≈ 0.31，与普朗克CMB结果完美吻合。

超新星Ia（Type Ia Supernovae）作为“标准烛光”，通过测量其视亮度与红移关系，给出宇宙距离-红移关系。1998年，两个团队（Supernova Cosmology Project和High-z Supernova Search Team）发现高红移超新星比预期暗，表明宇宙膨胀在加速。加速膨胀需要负压成分，即暗能量。超新星数据直接约束暗能量状态方程参数w = p/ρ，若w = -1则对应宇宙学常数Λ。结合CMB和BAO，超新星数据确认Ω_Λ ≈ 0.7且w ≈ -1。这三类观测（CMB、LSS/BAO、SNe Ia）分别对应宇宙不同 epoch（早期、中间、晚期），却收敛于同一套参数，这是ΛCDM模型最有力的支撑。

1.4 历史视角：从星云到精确宇宙学

人类对宇宙组成的认知经历了漫长演变。1785年，威廉·赫歇尔通过望远镜计数恒星，首次尝试绘制银河系结构，他误将银河系视为整个宇宙，但开启了“星海”观念。18世纪末，伊曼努尔·康德（Immanuel Kant）和约翰·赫歇尔（John Herschel）提出“宇宙岛”假说，认为某些星云是遥远的类似银河系的星系。1920年代，哈勃（Edwin Hubble）利用威尔逊山天文台的胡克望远镜，证实仙女座大星云是独立星系，并发现星系退行速度与距离成正比（哈勃定律），为宇宙膨胀奠定基础。1927年，勒梅特（Georges Lemaître）提出“原始原子”假说，认为宇宙从高密度状态开始膨胀；伽莫夫（George Gamow）、阿尔弗（Ralph Alpher）和赫尔曼（Robert Herman）在1948年预言了大爆炸核合成和CMB的存在。1965年CMB的发现使大爆炸理论成为主流。

暗物质的概念更早浮现。1933年，兹威基（Fritz Zwicky）研究后发座星系团时，利用维里定理发现星系速度弥散远大于可见物质引力所能束缚，提出“暗物质”存在。1970年代，薇拉·鲁宾（Vera Rubin）和肯特·福特（Kent Ford）测量螺旋星系旋转曲线，发现外层恒星旋转速度不随半径下降（与开普勒定律不符），而是近乎平坦，暗示星系被巨大暗物质晕包围。此后，引力透镜、星系团X射线观测等进一步证实暗物质普遍存在。暗能量则更晚：1998年前，宇宙学界普遍预期膨胀会因物质引力而减速，超新星结果颠覆了这一认知，迫使引入加速膨胀成分，最初称为“暗能量”，后与宇宙学常数关联。

从定性猜测（“缺失质量”）到定量精确测量，得益于技术飞跃。CMB测量从COBE（角分辨率7°，ΔT/T精度30 μK）到WMAP（角分辨率13角分，ΔT/T精度约15 μK）再到普朗克（角分辨率5角分，ΔT/T精度约5 μK）。普朗克卫星的低温辐射计在九个频率通道（30–857 GHz）工作，系统误差控制达前所未有的水平。这些数据通过复杂的前景去除（银河系同步辐射、尘埃辐射、点源）和校准，生成多张CMB成分图（如SMICA、NILC、SEVEM），其一致性本身即是对系统误差的检验。

1.5 技术深度：从ΔT/T到物质密度的转换

CMB温度涨落ΔT/T的测量精度达10^{-5}量级，但如何转化为物质密度？关键在于声学振荡的物理。在光子-重子耦合时期，密度涨落δ = δρ/ρ激发声波，角频率ω = k c_s / √3，其中c_s = c/√(3(1+R))是声速，R = 3ρ_b/(4ρ_γ)是重子-光子动量密度比。振荡方程近似为：

$$ \ddot{Θ} + \frac{\dot{a}}{a} \dot{Θ} + k^2 c_s^2 Θ = \frac{k^2}{3} Φ $$

其中Θ是温度扰动（约等于光子密度扰动），Φ是引力势（由暗物质主导）。在辐射主导时期，Φ衰减；物质主导时期，Φ恒定。再复合时，声波相位φ = k r_s，其中r_s是声学尺度（声波传播距离）：

$$ r_s = \int_0^{η_*} c_s dη $$

η是共形时间。观测到的C_l峰值对应k r_s ≈ nπ（n为整数）。因此，第一个峰值位置l_* ≈ π d_A / r_s，其中d_A是最后散射面角直径距离。d_A依赖于Ω_m和Ω_Λ，而r_s依赖于Ω_b h²和Ω_m h²。因此，l_*的测量同时约束几何（通过d_A）和物质成分（通过r_s）。第一个峰值高度主要取决于重子密度：Ω_b h²越大，R越大，c_s越小，声波振荡阻尼越强，但第一个峰值相对幅度增加。第二个峰值高度则对Ω_c h²敏感：暗物质增加引力势Φ，增强驱动，使第二个峰值相对第一个更高。第三个峰值进一步受光子扩散阻尼（Silk damping）影响，该阻尼尺度依赖于Ω_b h²和Ω_m h²。

通过计算理论C_l并拟合观测，可以解出这些参数。普朗克团队使用plik似然函数，直接比较观测多频温度（TT）、偏振（EE）及交叉（TE）功率谱与理论，在ΛCDM框架下进行MCMC采样。采样链收敛后，给出参数后验分布。例如，Ω_b h²与Ω_c h²存在轻微负相关，因为两者都影响峰值高度序列；但Ω_b h²主要控制第一与第二峰值比，Ω_c h²控制第二与第三峰值比，因此可区分。

1.6 实践案例：使用Planck公开数据计算宇宙组成

普朗克合作组公开了大量数据产品，使任何具备基础编程能力的研究者都能复现核心结果。以下是一个简化实践流程，使用Python和公开数据计算宇宙组分比例。

步骤1：获取数据 访问Planck Legacy Archive (PLA, https://pla.esac.esa.int)，下载：

• CMB温度功率谱：`COM_PowerSpect_CMB-base-plik-lite-l3.cls`（plik_lite组合的TT+TE+EE）
• 理论模板（可选）：CAMB或CLASS输出的理论谱（用于比较）
• 参数最佳拟合值文件：`base_plikHM_TTTEEE_lowl_lensing_latedata.paramnames`和`.mcmc`

步骤2：理解数据格式 `.cls`文件为文本，列依次为：l, D_l = l(l+1)C_l/(2π) [μK²] for TT, TE, EE, 以及对应误差（通常为对角线协方差）。注意D_l单位是μK²，需转换为C_l：C_l = 2π D_l / [l(l+1)]。

步骤3：计算观测功率谱（简化） 实际上，plik_lite已经提供直接可用的C_l，但为演示，可读取并转换为C_l：

```python import numpy as np import matplotlib.pyplot as plt

读取plik_lite数据

data = np.loadtxt('COM_PowerSpect_CMB-base-plik-lite-l3.cls') l = data

宇宙的绝大部分是隐形的。当我们用望远镜凝视星空，看到的恒星、星系、星云——所有这些发光物质——加起来仅占宇宙总质能的5%。剩下的95%由两种神秘成分主导：暗物质（Dark Matter）约占27%，暗宇宙学常数（Dark Energy）约占68%。这一惊人结论并非科幻猜想，而是基于过去三十年精密宇宙学观测的坚实推论，其核心支柱是宇宙微波背景辐射（Cosmic Microwave Background, CMB）的测量，特别是欧洲空间局普朗克（Planck）卫星2018年发布的数据。这些数字构成了现代宇宙学的标准模型——ΛCDM（Lambda Cold Dark Matter），它像一张宇宙的“成分清单”，却揭示了我们对宇宙本质的认知存在巨大空白。本章将深入剖析这张清单如何被绘制出来，从CMB的微小温度涨落到宇宙的大尺度结构，展示如何将天空中的光点转化为精确的宇宙学参数，并回顾人类如何从赫歇尔（William Herschel）的星云图一步步走到今天的宇宙学精确时代。

1.1 宇宙学标准模型ΛCDM：一个简洁而成功的框架

ΛCDM模型是当前描述宇宙演化的最简有效理论。它假设宇宙在大尺度上均匀各向同性，物质成分包括普通重子物质（Baryons）、冷暗物质（Cold Dark Matter）和暗能量（表现为宇宙学常数Λ），空间几何由爱因斯坦场方程决定。模型的六个基本参数通常取为：重子物理密度Ω_b h²、冷暗物质物理密度Ω_c h²、再复合时期的声学角度θ_*、再电离光学深度τ、标量谱指数n_s以及标量振幅A_s。其中h = H₀/100 km/s/Mpc是归一化哈勃常数。通过普朗克2018年数据，这些参数被精确测定，进而推导出今天的组分比例：

普通物质（重子）密度：Ω_b = Ω_b h² / h² ≈ 0.049 冷暗物质密度：Ω_c = Ω_c h² / h² ≈ 0.265 暗能量密度（Λ）：Ω_Λ ≈ 0.686 曲率密度：Ω_k = 1 - (Ω_b + Ω_c + Ω_Λ) ≈ 0（宇宙平坦）

这些比例之和严格满足宇宙几何闭合条件：Ω_m + Ω_Λ + Ω_k = 1，其中Ω_m = Ω_b + Ω_c ≈ 0.314。普朗克2018年给出的最佳拟合值（基于plik_lite组合）为：

$$ Ω_b h^2 = 0.02237 \pm 0.00015, \quad Ω_c h^2 = 0.1200 \pm 0.0012, \quad 100θ_* = 1.04110 \pm 0.00030, \quad τ = 0.0544 \pm 0.0073, \quad \ln(10^{10} A_s) = 3.044 \pm 0.014, \quad n_s = 0.9649 \pm 0.0042 $$

结合H₀ = 67.4 ± 0.5 km/s/Mpc（普朗克TT,TE,EE+lowE+lensing），计算得Ω_b ≈ 0.0486，Ω_c ≈ 0.2589，Ω_Λ ≈ 0.6925。注意不同数据组合（如加入BAO）会略有调整，但核心结论稳定：普通物质仅约5%，暗物质约27%，暗能量约68%。这张“成分表”的权威性来自多个独立观测的交叉验证，而CMB是最关键的“定锚点”。

1.1.1 为什么使用物理密度Ωh²？参数退化与h的困境

初学者常困惑：为何ΛCDM常用Ω_b h²和Ω_c h²而非直接Ω_b、Ω_c？这涉及参数间的强相关性（退化）。CMB功率谱主要敏感于物理密度（即乘以h²），因为声学尺度r_s ∝ 1/√(Ω_m h²)，而角直径距离d_A ∝ 1/H₀ ∝ 1/h。因此，l_* ∝ d_A/r_s 同时依赖h和Ω_m h²，导致Ω_m h²和h存在负相关：若h较大，则需较小Ω_m h²以保持l_*不变。直接使用Ω_m和h会导致拟合不确定度剧增。采用物理密度Ω_m h²可部分解耦，使约束更紧。但最终要得到今天的组分比例Ω_b、Ω_c，仍需独立测量H₀。这正是宇宙学“H₀张力”的根源：普朗克CMB给出的H₀（67.4 km/s/Mpc）与局部距离阶梯法（如SH0ES团队用造父变星校准超新星）给出的H₀（约73 km/s/Mpc）存在超过4σ的偏差。若采用局部H₀，则从相同CMB数据反推的Ω_Λ会略有降低，但95%未知的整体比例依然成立——普通物质仍只占约5%，暗物质和暗能量仍是主体。H₀张力可能暗示新物理（如早期暗能量），但未动摇ΛCDM的基本成分划分。

1.2 宇宙微波背景辐射：大爆炸的余晖与信息宝库

CMB是大爆炸后约38万年时释放的光子。此时宇宙温度降至约3000K，电子与质子结合成中性氢原子（再复合时期），光子得以自由传播，不再频繁散射。这些光子至今仍充满宇宙，但由于宇宙膨胀约1100倍，其波长被拉伸到微波波段，黑体谱峰值对应于约2.725K的温度。1965年，阿诺·彭齐亚斯（Arno Penzias）和罗伯特·威尔逊（Robe