天文学范式与现代宇宙观 / Astronomical Paradigms and the Modern Cosmos

现代宇宙学作为一门精确科学，其核心框架建立在几个里程碑式的观测发现之上。从哥白尼颠覆地心说的科学革命，到哈勃揭示宇宙膨胀的观测突破，再到普朗克卫星对宇宙微波背景辐射的精细测绘，人类对宇宙的理解经历了深刻的范式转换。本章将系统追溯这一思想史脉络，深入解析红移与哈勃定律的物理本质与观测基础，详细阐述宇宙距离阶梯的构建原理与关键技术（包括造父变星的周光关系、Ia型超新星的标准化校准），并基于普朗克数据讨论现代宇宙学参数的最佳估计值及其内在张力。最后，通过Python实践演示如何运用这些概念进行宇宙学距离计算与哈勃图构建，为后续章节的深入讨论奠定基础。

历史脉络：从地心说到宇宙膨胀

人类宇宙观的演变是一部科学革命史。古希腊学者阿里斯塔克斯（Aristarchus）曾提出日心说雏形，但托勒密（Ptolemy）的地心说模型凭借其数学精致性主导西方思想长达十四个世纪。1543年，哥白尼（Copernicus）在《天体运行论》中系统阐述日心说，开启科学革命序幕。伽利略（Galileo）以望远镜观测到金星相位、木星卫星等证据，直接挑战地心体系；开普勒（Kepler）基于第谷·布拉赫（Tycho Brahe）的观测数据，提出行星运动三大定律，为牛顿（Newton）万有引力定律铺平道路。牛顿力学成功解释了太阳系动力学，但当时人们认为“银河即宇宙”，太阳位于银河系中心。

18世纪末，威廉·赫歇尔（William Herschel）通过恒星计数首次尝试绘制银河系结构，误将太阳置于中心。20世纪初，哈罗·沙普利（Harlow Shapley）通过球状星团分布证明太阳远离银河系中心，银河系尺度被扩展至约10万光年。与此同时，埃德温·哈勃（Edwin Hubble）利用威尔逊山天文台的100英寸望远镜，在1923-1924年确认仙女座“星云”为河外星系（Andromeda Galaxy），揭示银河系仅是宇宙中千亿星系之一，彻底改变人类对宇宙尺度的认知。

1929年，哈勃发表里程碑论文，基于斯里弗（Vesto Slipher） earlier 的红移测量和自行距离估计，发现星系退行速度与距离近似成正比：$v \approx H_0 d$。这一关系被命名为哈勃定律（Hubble's Law），成为宇宙膨胀的直接证据。哈勃最初估计 $H_0 \approx 500$ km/s/Mpc，因距离校准误差而偏高；后续数十年间，随着距离阶梯改进，$H_0$ 值逐步下调至约70 km/s/Mpc量级。哈勃的发现催生了现代宇宙学，将宇宙从一个静态背景转变为动态演化实体。

红移的物理本质：多普勒效应与度规膨胀

红移（Redshift）定义为观测波长与静止波长之比减一：$z = (\lambda_{\text{obs}} - \lambda_{\text{emit}})/\lambda_{\text{emit}}$。在宇宙学中，红移主要源于两种物理机制：

1. 多普勒红移（Doppler Redshift）：由光源相对观测者的运动引起，适用于局部惯性系。对于非相对论速度 $v \ll c$，有 $z \approx v/c$；相对论情形下，$1+z = \sqrt{(1+v/c)/(1-v/c)}$。星系在本动速度（peculiar velocity）导致的多普勒红移，是哈勃定律观测中的主要噪声源。

2. 宇宙学红移（Cosmological Redshift）：由时空度规膨胀引起，是弗里德曼-勒梅特-罗伯逊-沃尔克（FLRW）度规的直接结果。在均匀各向同性宇宙中，标度因子 $a(t)$ 描述宇宙膨胀，光子波长随宇宙膨胀而拉伸：$1+z = a(t_{\text{obs}})/a(t_{\text{emit}})$。宇宙学红移并非传统多普勒效应，而是光子在膨胀时空中传播的累积效应，其速度-距离关系在广义相对论框架下由测地线方程描述。

对于小红移（$z \ll 1$），宇宙学红移与多普勒红移在数值上近似等价，$v_{\text{rec}} \approx c z$，其中 $v_{\text{rec}}$ 为哈勃流退行速度。但当 $z \gtrsim 0.1$ 时，必须考虑宇宙学效应：退行速度可能超过光速（不违反相对论，因非局部运动），且速度-距离关系非线性。精确的共动距离（comoving distance）计算需积分弗里德曼方程：

$$\chi(z) = c \int_0^z \frac{dz'}{H(z')}$$

其中哈勃参数 $H(z) = H_0 \sqrt{\Omega_m (1+z)^3 + \Omega_r (1+z)^4 + \Omega_\Lambda + \Omega_k (1+z)^2}$。光度距离（luminosity distance）与角直径距离（angular diameter distance）则通过 $d_L = (1+z) S_k(\chi)$ 和 $d_A = d_L/(1+z)^2$ 与共动距离关联，$S_k(\chi)$ 为空间曲率修正函数（$k=0$ 时为 $\chi$）。

哈勃定律的观测验证依赖于独立距离测量。哈勃最初使用造父变星校准仙女座星系距离，但受限于仪器和校准误差。现代哈勃常数测定需跨越宇宙距离阶梯，每一步都引入系统误差，导致当前“哈勃张力”（Hubble Tension）——局部测量值（$H_0 \approx 73$ km/s/Mpc）与CMB推断值（$H_0 \approx 67.4$ km/s/Mpc）存在超过 $4\sigma$ 的偏差，可能指向新物理或系统误差。

宇宙距离阶梯：技术原理与校准链

宇宙距离阶梯（Cosmic Distance Ladder）是一套分层测距体系，利用“标准烛光”（standard candles）或“标准尺”（standard rulers）逐步扩展可测距离范围。每级方法校准上一级，误差逐级累积。核心步骤如下：

3.1 三角视差法（Trigonometric Parallax）

适用于太阳系内及近距恒星（$d \lesssim 10$ kpc）。通过地球公转基线（2 AU）观测恒星周年视差角 $\pi$（角秒），距离 $d$（秒差距）满足 $d = 1/\pi$。依巴谷卫星（Hipparcos，1989-1993）精确测量约10万颗恒星， Gaia卫星（2013至今）将精度提升至微角秒级，视差精度 $\sigma_\pi/\pi \sim 0.02-0.05$ 对于 $G<15$ 等星，覆盖银河系大部分区域。视差法是距离阶梯的基石，直接校准造父变星等后续标准烛光。

3.2 造父变星（Cepheid Variables）

造父变星是脉动变星，其光度周期 $P$ 与绝对星等 $M$ 存在稳定周光关系（Period-Luminosity Relation）：$M = a \log P + b$。经典造父变星（Type I）属年轻星族I，金属丰度较高，位于旋臂；II型造父变星（W Vir stars）属年老星族II，金属丰度低，周光关系不同。 Henrietta Leavitt 1912年通过小麦哲伦云造父变星发现周光关系，因假设云内恒星距离相同而确立相对关系。1920年代，哈勃与沙普利利用银河系内造父变星（通过视差或星团拟合校准）确定绝对校准，扩展至仙女座星系。

现代造父变星测距使用近红外波段（如 $H$ 带）减少消光影响，并通过周期性光度曲线拟合提高精度。哈勃太空望远镜（HST）的造父变星观测（如 H0 Key Project）将距离阶梯延伸至约30 Mpc，校准Ia型超新星。系统误差主要来自周光关系斜率/截距的金属丰度依赖性、红化修正、光度曲线形状校正等。

3.3 Ia型超新星（Type Ia Supernovae）

Ia型超新星源于碳氧白矮星吸积伴星物质达到钱德拉塞卡极限（$\sim 1.4 M_\odot$）引发的热核爆炸。其峰值绝对星等均一（$M_B \approx -19.3$），是宇宙学距离的理想标准烛光，可达 $z \sim 1.5$。但亮度存在 $\sim 30\%$ 自然分散，需校正：

• **Phillips 关系**：光变曲线衰减率 $\Delta m_{15}(B)$（B波段亮度下降15天幅度）与峰值亮度相关，更慢衰减者更亮。
• **颜色校正**：红颜色暗示更多消光或本征红化。
• **宿主星系质量校正**：大质量宿主星系中的超新星略亮，可能与金属丰度相关。

校准链：通过近距Ia型超新星（$z<0.1$）与造父变星宿主星系距离比较，建立绝对光度标定。此过程由H0 Key Project和SH0ES团队推进，SH0ES团队利用HST观测宿主星系造父变星，报告 $H_0 = 73.0 \pm 1.0$ km/s/Mpc（2022）。系统误差包括超新星样本选择偏差、流量平均效应（Malmquist bias）、尘埃消光模型等。

3.4 距离阶梯的误差传递与哈勃张力

距离阶梯的系统误差沿校准链传递。例如，造父变星周光关系的零点误差直接影响超新星校准，进而影响 $H_0$。此外，阶梯方法在 $z>0.1$ 区域需假设宇宙学模型（如平坦ΛCDM）将红移转换为距离，可能引入模型依赖偏差。哈勃张力可能源于：局部速度场扰动（如本动星系影响）、阶梯系统误差未充分建模、或早期宇宙与晚期宇宙物理不一致（如早期暗能量、修改引力等新物理）。

```mermaid graph TD A[三角视差<br/>Gaia: d < 10 kpc] --> B[造父变星<br/>周光关系校准<br/>d ~ 30 Mpc] B --> C[Ia型超新星<br/>Phillips校正<br/>d > 100 Mpc] C --> D[哈勃常数 H0] D --> E[宇宙学参数拟合] E --> F[ΛCDM模型] style A fill:#e1f5fe style B fill:#f3e5f5 style C fill:#e8f5e9 style D fill:#fff3e0 ```

现代宇宙学参数：ΛCDM模型与普朗克数据

当前标准宇宙学模型为ΛCDM（Lambda Cold Dark Matter），包含六个基本参数，由宇宙微波背景（CMB）各向异性、大尺度结构、超新星等多探针联合约束。该模型假设：宇宙均匀各向同性（宇宙学原理），度规由FLRW描述，物质成分包括重子物质（baryons）、冷暗物质（cold dark matter）和宇宙学常数（Λ，对应暗能量），初始扰动为 adiabatic、尺度不变高斯随机场。

普朗克卫星（Planck，2009-2013）提供全天空CMB温度与偏振高精度图，测量角功率谱至多极 $\ell \sim 2500$。基于2018年数据释放（Planck 2018 results），ΛCDM最佳拟合参数为：

| 参数 | 符号 | 值 | 单位 | |------|------|-------|------| | 哈勃常数 | $H_0$ | $67.4 \pm 0.5$ | km/s/Mpc | | 重子密度参数 | $\Omega_b h^2$ | $0.0224 \pm 0.0001$ | — | | 冷暗物质密度参数 | $\Omega_c h^2$ | $0.120 \pm 0.001$ | — | | 再复合时期光学深度 | $\tau$ | $0.054 \pm 0.007$ | — | | 标量谱指数 | $n_s$ | $0.965 \pm 0.004$ | — | | 标量振幅 | $\ln(10^{10} A_s)$ | $3.044 \pm 0.014$ | — |

其中 $h = H_0 / (100 \text{ km/s/Mpc})$，当前物质密度参数 $\Omega_m = \Omega_b + \Omega_c \approx 0.315$，暗能量密度参数 $\Omega_\Lambda \approx 0.685$（平坦宇宙假设 $\Omega_{\text{tot}}=1$）。这些参数精确描述宇宙组成：约5%重子物质，27%冷暗物质，68%暗能量。

CMB功率谱的声学峰结构敏感于这些参数：第一峰位置决定空间曲率（Planck确认平坦性）；峰高比约束 $\Omega_b/\Omega_c$；再复合时期物理决定峰位置与振幅。普朗克数据与ΛCDM模型高度吻合，$\chi^2$/dof ≈ 1，但残差在低多极区（$\ell<30$）存在轻微异常，可能暗示宇宙学原理失效或新物理。

哈勃张力是当前最大挑战。局部距离阶梯法（如SH0ES）给出 $H_0 = 73.0 \pm 1.0$ km/s/Mpc，而Planck在ΛCDM框架下推断 $H_0 = 67.4 \pm 0.5$ km/s/Mpc，偏差达 $4.4\sigma$（2023年统计）。可能解释包括：未校正的系统误差（如造父变星金属丰度效应、超新星光污染）、局部宇宙速度场扰动（如 $z<0.1$ 区域存在 $\sim 300$ km/s 本动速度）、或早期宇宙物理修改（如早期暗能量、额外 relativistic 物种、修改膨胀历史）。解决此张力需独立测距方法，如引力波标准汽笛（GW170817给出 $H_0 = 70^{+12}_{-8}$ km/s/Mpc，精度待提高）、强引力透镜时间延迟（H0LiCOW项目 $H_0 = 73.3^{+1.7}_{-1.8}$ km/s/Mpc）、或泰勒序列校准的造父变星独立分析。

Python实践：宇宙学距离计算与哈勃图构建

本实践使用 `astropy` 库（版本5.0+）进行宇宙学距离计算与可视化。`astropy.cosmology` 提供预定义宇宙学模型（如 Planck18、WMAP9）及距离计算函数。以下代码演示光度距离、角直径距离计算，并绘制不同宇宙学模型的哈勃图（速度-距离关系）。

```python import numpy as np import matplotlib.pyplot as plt from astropy.cosmology import Planck18 as cosmo_planck from astropy.cosmology import FlatLambdaCDM import astropy.units as u

定义红移范围（小红移，哈勃定律近似有效）

z = np.linspace(0, 0.2, 100) # 红移0到0.2

计算光度距离（Mpc）

d_l_planck = cosmo_planck.luminosity_distance(z).value

自定义宇宙学模型：H0=73 km/s/Mpc, Ωm=0.3, ΩΛ=0.7

cosmo_local = FlatLambdaCDM(H0=73 * u.km / u.s / u.Mpc, Om0=0.3) d_l_local = cosmo_local.luminosity_distance(z).value

计算哈勃流速度 v = H0 * d（仅小红移近似）

v_planck = cosmo_planck.H0.value * d_l_planck # 注意：H0单位 km/s/Mpc，d_l单位Mpc，v单位km/s v_local = cosmo_local.H0.value * d_l_local

绘制哈勃图

plt.figure(figsize=(10, 6)) plt.plot(d_l_planck, v_planck, 'b-', label=f'Planck18: $H_0$={cosmo_planck.H0.value:.1f}') plt.plot(d_l_local, v_local, 'r--', label=f'Local: $H_0$={cosmo_local.H0.value:.1f}') plt.xlabel('Luminosity Distance $d_L$ (Mpc)') plt.ylabel('Recession Velocity $v$ (km/s)') plt.title('Hubble Diagram for Low Redshift ($z < 0.2$)') plt.legend() plt.grid(True, alpha=0.3) plt.show()

宇宙学原理与FLRW度规的数学表述

现代宇宙学建立在两个观测基石之上：宇宙学原理（Cosmological Principle）与广义相对论（General Relativity）。宇宙学原理假设在大尺度上（$>100$ Mpc），宇宙在空间上是均匀（homogeneous）与各向同性（isotropic）的。这一假设得到多重观测支持：宇宙微波背景（CMB）温度各向异性仅约 $10^{-5}$；大尺度星系巡天（如 Sloan Digital Sky Survey, SDSS）显示星系分布呈“宇宙网”结构，但统计平均后密度涨落随尺度增大而趋平；类星体分布亦无优先方向。尽管局部存在星系团、超星系团等结构，但宇宙学原理在大于 $100$ Mpc 尺度上成立，为度规构建提供对称性约束。

在广义相对论框架下，物质-能量分布决定时空几何。满足宇宙学原理的度规由弗里德曼-勒梅特-罗伯逊-沃尔克（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$ 为空间曲率常数（$k=+1,0,-1$ 对应闭合、平坦、开放宇宙）；$(r,\theta,\phi)$ 为共动坐标（comoving coordinates），随宇宙膨胀而拉伸。物理距离由 $a(t)$ 缩放：共动距离 $\chi$ 不变，proper distance $d_p = a(t) \chi$。度规的时空对称性导致爱因斯坦场方程简化为弗里德曼方程（Friedmann Equations）：

第一方程（膨胀率）： $$H(t)^2 \equiv \left( \frac{\dot{a}}{a} \right)^2 = \frac{8\pi G}{3} \rho - \frac{kc^2}{a^2} + \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}$$

其中