宇宙学标准模型与暗成分的引入
现代宇宙学标准模型(ΛCDM)的成功建立,是20世纪理论物理学与观测天文学协同发展的巅峰之作,其核心在于引入两种神秘的暗成分——暗物质与暗能量,以 reconcile 一系列看似矛盾的观测事实。该模型以极其简洁的框架——一个包含宇宙学常数(Λ)与冷暗物质(CDM)的平直弗里德曼-勒梅特-罗伯逊-沃尔克(FLRW)度规——精确描述了从宇宙最初时刻(约38万年)至今日(约138亿年)的演化历史,其预测与宇宙微波背景辐射(CMB)、Ia型超新星(SNe Ia)、重子声波振荡(BAO)及大尺度结构(LSS)等多类独立观测数据达到了前所未有的吻合程度。然而,这一“ concordance model ” 的辉煌背后,暗物质与暗能量的物理本质始终是笼罩在现代物理学上空的最大谜团,而近年来自不同观测团队对关键参数(如哈勃常数 $H_0$ 和物质涨落幅度 $S_8$)的测量结果之间暴露出的显著分歧,正剧烈动摇着模型看似稳固的基石,迫使学界重新审视其基础假设。本章将系统追溯ΛCDM模型的建立历程,从爱因斯坦场方程出发,阐释弗里德曼方程如何成为宇宙动力学的核心描述;通过对比分析CMB各向异性功率谱与SNe Ia距离-红移关系等关键数据集,揭示暗物质与暗能量被引入的历史必然性与观测必要性;最后,深入剖析当前模型所面临的核心参数张力,为后续章节对暗成分本质的深入探索奠定坚实的理论基础。
理论基石:FLRW度规与弗里德曼方程的宇宙动力学
现代宇宙学建立在广义相对论的基础之上。爱因斯坦场方程将时空几何与物质-能量分布联系起来: $$G_{\mu\nu} + \Lambda g_{\mu\nu} = \frac{8\pi G}{c^4} T_{\mu\nu}$$ 其中 $G_{\mu\nu}$ 是爱因斯坦张量,$g_{\mu\nu}$ 是度规张量,$T_{\mu\nu}$ 是能量-动量张量,$\Lambda$ 是宇宙学常数。为了描述整体均匀各向同性的宇宙,必须对度规施加强宇宙学原理约束,即在大尺度上,宇宙在空间每一点、朝每一个方向都统计相同。满足此条件的时空度规即FLRW度规,其线元可写为: $$ds^2 = -c^2 dt^2 + a^2(t) \left[ \frac{dr^2}{1 - k r^2} + r^2 (d\theta^2 + \sin^2\theta d\phi^2) \right]$$ 其中 $a(t)$ 是尺度因子,描述宇宙随时间的膨胀;$k$ 是归一化曲率常数,取值为 $+1, 0, -1$ 分别对应封闭、平直、开放的几何。将FLRW度规代入场方程,并假设宇宙内容为完美流体($T^{\mu}_{\;\nu} = \text{diag}(-\rho c^2, p, p, p)$),即可导出描述宇宙整体演化的核心方程——第一弗里德曼方程(或哈勃方程): $$H^2(t) \equiv \left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho(t) - \frac{k c^2}{a^2(t)} + \frac{\Lambda c^2}{3}$$ 第二弗里德曼方程(加速度方程)为: $$\frac{\ddot{a}}{a} = -\frac{4\pi G}{3} \left( \rho + \frac{3p}{c^2} \right) + \frac{\Lambda c^2}{3}$$ 这两方程是宇宙学标准模型的动力学骨架。$H(t)$ 是哈勃参数,其当前值 $H_0$ 为哈勃常数。方程右侧第一项代表物质(及相对论性成分)的引力效应(吸引),第二项是空间曲率贡献,第三项是宇宙学常数(排斥)。宇宙的膨胀历史完全由其中能量成分的密度随尺度因子的演化规律 $\rho(a)$ 所决定。
对于非相对论性物质(“冷”物质),其密度随体积膨胀而稀释,$\rho_m \propto a^{-3}$。对于相对论性成分(如辐射、中微子),$\rho_r \propto a^{-4}$(多出一个红移因子)。对于宇宙学常数 $\Lambda$,其能量密度 $\rho_\Lambda = \Lambda c^2 / (8\pi G)$ 为常数,不随膨胀变化,故其状态方程参数 $w = p/\rho c^2 = -1$。此外,若存在曲率,其贡献 $\rho_k \propto a^{-2}$。因此,当前时刻的哈勃参数可表示为: $$H_0^2 = H_0^2 \left[ \Omega_{m,0} a^{-3} + \Omega_{r,0} a^{-4} + \Omega_{k,0} a^{-2} + \Omega_{\Lambda,0} \right]$$ 其中 $\Omega_{i,0} = \rho_{i,0} / \rho_{\text{crit},0}$ 是各成分当前密度参数,$\rho_{\text{crit},0} = 3H_0^2/(8\pi G)$ 是当前临界密度。$\sum \Omega_{i,0} = 1$ 对应平直几何($k=0$),这是暴胀理论预言的典型结果,也是当前观测强烈支持的结论。
弗里德曼方程揭示了一个关键物理图景:宇宙的膨胀命运由总能量密度及其状态方程决定。若只有物质($\Omega_m=1$),膨胀将永远减速;若存在正 $\Lambda$($\Omega_\Lambda>0$),则当物质密度稀释到一定程度后,排斥性的 $\Lambda$ 将主导,导致膨胀加速。观测正是通过测量不同红移处的 $H(z)$ 或等效的距离-红移关系,来反推各 $\Omega_i$ 的值。
观测转折点:从静态宇宙到动态膨胀
哈勃定律与宇宙膨胀的发现
1929年,埃德温·哈勃利用威尔逊山天文台的胡克望远镜,系统测量了 nearby 星系的光谱红移与距离,发现红移 $z$ 与星系距离 $d$ 近似成正比:$v \approx H_0 d$,其中 $v \approx cz$(对于小红移)。这一哈勃定律首次提供了宇宙正在膨胀的观测证据,彻底否定了爱因斯坦曾坚持的静态宇宙模型(为此他后来称引入 $\Lambda$ 为“一生最大错误”,尽管其物理意义在后来被重新认识)。哈勃的原始 $H_0$ 值约为 $500 \text{ km/s/Mpc}$,远高于现代值,这主要受限于当时对造父变星周光关系定标的系统误差。此后数十年,通过改进距离阶梯(如利用Ia型超新星作为“标准烛光”),$H_0$ 的测量精度不断提高,但其值始终在 $50-100 \text{ km/s/Mpc}$ 区间内存在争议,直接影响了宇宙年龄的估算。
宇宙年龄危机与暗物质的早期暗示
在哈勃时代,若取 $H_0 \sim 50-100 \text{ km/s/Mpc}$ 且宇宙仅由物质主导($\Omega_m=1$),由弗里德曼方程积分可得的宇宙年龄 $t_0 \approx 2/(3H_0)$ 仅为 $60-130$ 亿年。然而,当时对恒星年龄(特别是球状星团)的估计已普遍超过 $100$ 亿年,甚至达 $150-200$ 亿年,这构成了严重的“宇宙年龄危机”。解决之道有二:要么降低 $H_0$,要么引入非零 $\Lambda$(使膨胀历史变慢,从而增加 $t_0$),要么承认存在大量不可见的物质(增加 $\Omega_m$ 但需调整 $H_0$)。实际上,早在1933年,弗里茨·兹威基在研究后发座星系团时,通过维里定理计算其动力学质量,发现远大于可见恒星质量之和,首次提出“暗物质”概念。但此工作因距离测量不准及对星系团动力学理解不足,在当时未受重视。直到1970年代,薇拉·鲁宾等人对旋涡星系旋转曲线的精确测量,提供了暗物质存在于星系尺度的强有力证据:星系外围恒星和气体的旋转速度并未像预期那样随半径增加而下降(开普勒定律),而是几乎保持平坦,这意味着星系质量分布远超出可见盘,存在巨大的暗物质晕。这些动力学异常,若仅用可见物质解释,则要求星系处于非平衡态或修改引力理论(如MOND),但前者与星系团观测矛盾,后者在宇宙学尺度上难以成功。暗物质假设以最简单的粒子物理模型(弱相互作用大质量粒子,WIMP)即可解释所有尺度上的引力效应,且与宇宙大尺度结构的形成理论自然契合。
暗物质在宇宙学标准模型中的必要性:结构形成的引擎
ΛCDM模型中的“CDM”即冷暗物质,指其运动速度远低于光速(非相对论性),从而能在早期宇宙中有效坍缩形成密度扰动种子。宇宙大爆炸后约38万年,光子与重子脱耦,形成我们观测到的CMB。此时,暗物质已因不受辐射压力影响而开始在其自身引力下增长微小的密度涨落,而重子物质仍被光子“粘滞”着,只能跟随暗物质势阱运动。脱耦后,重子才落入暗物质形成的势阱中,最终形成星系和星系团。这一“先暗后重”的层次结构形成图像,是ΛCDM的核心预言。
CMB各向异性功率谱为这一图像提供了最精确的验证。COBE卫星首次探测到CMB温度涨落($10^{-5}$ 量级),WMAP和普朗克卫星则给出了高精度全 sky 图谱。功率谱(温度涨落方差随角尺度的分布)呈现一系列峰值,其位置、高度和间隔对宇宙参数极为敏感。第一峰的位置主要对总物质密度(包括暗物质)和空间曲率敏感。观测显示第一峰位于多极矩 $\ell \approx 220$,对应角尺度约 $1^\circ$,这精确对应于平直几何($\Omega_{\text{total}} \approx 1$)的预期。第一峰与第二峰的高度比则对重子密度 $\Omega_b h^2$ 与物质密度 $\Omega_m h^2$ 的比值敏感:重子越多,声波振荡阻尼越强,第二峰相对第一峰越低。普朗克2018年结果给出 $\Omega_b h^2 \approx 0.0224$,$\Omega_c h^2 \approx 0.120$,其中 $\Omega_c$ 即为冷暗物质密度参数。若无暗物质($\Omega_c=0$),仅靠重子,则第一峰高度会显著降低,且峰值位置会因声速变化而移动,与观测严重不符。第三峰及以后峰值的相对高度则对暗物质总量和初始扰动谱指数 $n_s$ 敏感,进一步确认了暗物质的主导作用。
此外,CMB功率谱的低 $\ell$ 部分(大角尺度) 受到早期宇宙积分萨克斯-沃尔夫效应的影响,对 $\Omega_\Lambda$ 有间接约束,但精度有限。暗物质的存在不仅解释了CMB峰值结构,也是宇宙大尺度结构(星系巡天如SDSS、2dFGRS观测到的星系分布相关函数和功率谱)形成的唯一已知机制。数值模拟(如 Millennium Simulation)在ΛCDM框架下,能够重现观测到的星系相关函数、星系团质量函数、以及宇宙网(宇宙纤维状结构)的形态。若去除暗物质,仅靠重子,结构形成的时间尺度将远长于宇宙年龄,无法解释高红移处已存在的巨大星系团和类星体。
暗能量的发现:加速膨胀与宇宙学常数的回归
Ia型超新星作为标准烛光
Ia型超新星源于白矮星吸积伴星物质达到钱德拉塞卡极限(约 $1.4 M_\odot$)后的热核爆炸。其峰值光度在爆炸物理机制上具有高度一致性,经过光变曲线形状校正后,可作为“标准烛光”测量共动距离。1990年代中期,两个独立团队——超新星宇宙学计划(SCP)和高红移超新星搜索团队(HZT)——开始利用地面望远镜和哈勃太空望远镜,系统观测高红移($z \sim 0.5-1$)Ia型超新星,并与低红移($z<0.1$)样本比较,旨在测量宇宙膨胀历史的减速参数。1998年,两队几乎同时宣布震惊学界的发现:高红移超新星比预期更暗,意味着它们比在减速膨胀宇宙中预测的距离更远。在给定哈勃常数和物质密度下,这唯一合理的解释是宇宙膨胀在近期(约过去50-60亿年)从减速转为加速。根据加速度方程,加速要求 $\ddot{a}>0$,即 $\rho + 3p/c^2 < 0$。对于普通物质($w \geq 0$),其贡献为正,只能减速。因此,必须存在一种具有负压的组分,其状态方程 $w < -1/3$,且在当前能量密度中占主导。最简单的候选就是爱因斯坦的宇宙学常数 $\Lambda$,对应 $w=-1$,其能量密度恒定,压强 $p_\Lambda = -\rho_\Lambda c^2$。
其他独立证据的汇合
超新星的发现迅速被其他观测所支持: 1. 宇宙微波背景辐射:如前所述,CMB功率谱峰值位置对总密度 $\Omega_{\text{total}}$ 敏感。WMAP和普朗克的测量一致表明 $\Omega_{\text{total}} \approx 1$。结合大尺度结构(如星系巡天)对物质密度 $\Omega_m$ 的测量($\Omega_m \approx 0.3$),由 $\Omega_{\text{total}} = \Omega_m + \Omega_\Lambda + \Omega_k$ 且 $\Omega_k \approx 0$,立即推出 $\Omega_\Lambda \approx 0.7$。这是独立于超新星的、对暗能量存在的强有力佐证。 2. 重子声波振荡:BAO是早期宇宙中声波在光子-重子流体中传播留下的特征尺度(约 $150 \text{ Mpc}$),作为“标准尺”出现在星系分布的相关函数中。测量不同红移处BAO的角尺度或径向尺度,可独立测出 $D_A(z)$(角直径距离)和 $H(z)$。BAO数据与CMB、SNe结合,进一步精确了 $\Omega_\Lambda$ 并确认加速膨胀。 3. 大尺度结构与弱引力透镜:星系成团性和弱透镜测量对宇宙增长历史敏感。观测到的结构增长速率与ΛCDM预言一致,前提是存在暗能量抑制了晚期增长(因加速膨胀稀释了物质密度)。
这些独立数据集——CMB(早期宇宙快照)、SNe Ia(晚期几何)、BAO(中间红移几何)、LSS(增长历史)——在ΛCDM参数空间内实现了令人惊叹的交叉验证,构成了现代宇宙学标准模型的观测基石。
参数测量:ΛCDM模型的精确认定
当前,ΛCDM模型通常用6个基础参数描述(有时增加 $N_{\text{eff}}$、$\sum m_\nu$ 等):
- $\Omega_b h^2$:重子物理密度
- $\Omega_c h^2$:冷暗物质物理密度
- $\theta_*$:声视尺度(由CMB第一峰位置决定,与 $\Omega_{\text{total}}$ 相关)
- $\tau$:再电离光深
- $n_s$:标量扰动谱指数
- $\ln(10^{10} A_s)$:扰动幅度(或 $A_s$)
其中 $h = H_0 / (100 \text{ km/s/Mpc})$ 是无量纲哈勃常数。普朗克2018年联合CMB、BAO数据给出的最佳拟合值为(TT,TE,EE+lowE+lensing+BAO):
- $H_0 = 67.4 \pm 0.5 \text{ km/s/Mpc}$
- $\Omega_m = 0.315 \pm 0.007$
- $\Omega_\Lambda = 0.685 \pm 0.007$(由平直性约束)
- $\Omega_b h^2 = 0.02237 \pm 0.00015$
- $\Omega_c h^2 = 0.1200 \pm 0.0012$
- $n_s = 0.9649 \pm 0.0042$
- $\sigma_8 = 0.811 \pm 0.006$(线性理论 today 的密度涨落幅度)
这些参数以约 $0.1-1\%$ 的精度被确定,其一致性是ΛCDM作为“标准模型”的核心依据。例如,由CMB推演的 $H_0$ 与BAO数据结合独立测得的 $H_0$ 高度一致,共同支持了平直几何和 $\Omega_\Lambda \approx 0.7$。模型对CMB功率谱的拟合优度极高,$\chi^2$ 每自由度接近1,残差无显著系统偏差。
模型的内在张力:H₀危机与S₈危机
尽管ΛCDM在整体上表现卓越,近年来两类关键参数在不同测量技术间的显著分歧,已构成对模型基本假设的严峻挑战。
H₀ 危机:早期与晚期宇宙的“哈勃常数之争”
早期宇宙推断值:基于普朗克CMB数据,在ΛCDM框架下,通过将早期观测($z \sim 1100$)外推至 today,得到 $H_0^{\text{CMB}} = 67.4 \pm 0.5 \text{ km/s/Mpc}$。此值依赖于模型假设(特别是暗能量为 $w=-1$ 的 $\Lambda$,且无额外 relativistic 成分)。 晚期宇宙直接测量值:采用“距离阶梯”法,先校准 nearby 造父变星的周光关系(通过几何距离如 Gaia 卫星),再用这些造父变星标定 nearby Ia型超新星的绝对亮度,最后测量高红移Ia型超新星的距离以推算 $H_0$。最著名的SH0ES团队(Riess等)2022年最新结果为 $H_0^{\text{SH0ES}} = 73.04 \pm 1.04 \text{ km/s/Mpc}$。两者相差约 $9\%$,统计显著性超过 $5\sigma$($>6\sigma$ 的联合分析)。其他独立方法,如基于强引力透镜时间延迟(H0LiCOW/TDCOSMO)、基于晚期宇宙几何(如 megamaser 距离)的测量,大多支持 $H_0 \gtrsim 70$,与CMB推断值存在 $2-3\sigma$ 分歧。
这一张力若持续存在,可能指向: 1. 系统误差:两类测量中可能存在未识别的系统偏差。例如,造父变星周光关系的零点校准、Ia型超新星校正的演化效应、星际消光、引力透镜对超新星亮度的影响等。CMB分析中,对前景剔除、束校准、光度方程假设等也可能存疑。但目前双方系统误差评估均难以解释 $5\sigma$ 差异。 2. 新物理:在ΛCDM之外引入额外成分或修改早期宇宙物理。例如:
- **早期暗能量(EDE)**:在再电离前后($z \sim 5000$)存在一段短暂的暗能量主导期,其占比约 $10\%$,可增加早期膨胀速率,减小声视尺度 $\theta_*$,从而在相同 $\theta_*$ 下要求更大的 $H_0$ 以匹配观测。EDE模型能缓解 $H_0$ 张力,但可能引入其他问题(如与CMB极化或LSS冲突)。
- **额外 relativistic 成分**($\Delta N_{\text{eff}}$):增加有效中微子种类或某种轻粒子,可增加早期辐射密度,同样减小 $\theta_*$,从而提升 $H_0$。但普朗克对 $N_{\text{eff}}$ 的约束已很严($N_{\text{eff}} = 2.99 \pm 0.17$),空间有限。
- **修改暴胀或早期宇宙物理**:如改变再复合物理、引入非标准声速等。
- **修改引力理论**:在宇宙学尺度上偏离广义相对论,可能改变 $H(z)$ 演化。
目前,尚无一个模型能完美解决 $H_0$ 张力而不产生