第13章 氢原子：量子力学的皇冠明珠

故事引入：设计师的终极谜题

2156年，"原子铸造厂"的首席设计师艾琳接到了一个不可能的任务：用最简单的原料——一个带正电的质子核和一个带负电的电子——设计一个稳定的"量子建筑"。她的团队在三维设计空间中尝试了无数方案，但每次模拟都得出同一个结论：电子不会像经典行星那样盘旋坠落，而是形成一系列分立的、美丽的概率云。更令人震惊的是，这些概率云的能量层级只由一个整数 $n$ 决定，而形状则由两个量子数 $l$ 和 $m$ 描述——但它们竟然共享相同的能量！"大自然在设计原子时，"艾琳在最终报告中写道，"用最少的材料创造了最丰富的结构。一个质子，一个电子，无穷无尽的化学宇宙。"从 Bohr 的半经典模型到 Schrödinger 方程的精确解，氢原子是量子力学从玩具模型走向真实世界的转折点。它的求解不仅是数学上的胜利，更是人类理解物质结构的基石。

---

13.0 前置知识：中心力场问题的经典力学回顾

在用量子力学求解氢原子之前，让我们先在经典力学中理解中心力场问题的结构。经典图像为量子求解提供了坐标选择和分离变量的直觉。

13.0.1 开普勒问题：经典遗产

经典力学中的中心力场问题由牛顿引力或库仑力驱动。一个质量为 $m$ 的粒子在势能 $V(r) = -k/r$ 中运动（$k > 0$ 对应吸引力）。

关键守恒量：

1. 能量 $E = \frac{p^2}{2m} - \frac{k}{r}$
2. 角动量 $\mathbf{L} = \mathbf{r} \times \mathbf{p}$（向量守恒，方向和大小都守恒）
3. Runge-Lenz矢量 $\mathbf{A} = \frac{1}{mk}(\mathbf{p} \times \mathbf{L}) - \frac{\mathbf{r}}{r}$（仅对 $1/r$ 势存在！）

Runge-Lenz矢量的守恒意味着轨道是闭合的椭圆（没有进动），这是 $1/r$ 势的特殊性质。在量子力学中，Runge-Lenz矢量的存在将解释氢原子"意外的" $l$ 简并。

13.0.2 有效势与轨道形状

利用角动量守恒，运动可以约化到径向方程。引入有效势：

$V_{\text{eff}}(r) = \frac{L^2}{2mr^2} - \frac{k}{r}$

第一项是离心势（排斥性），第二项是库仑吸引。有效势的形状决定了轨道类型：

• $E < 0$：束缚态，椭圆轨道
• $E = 0$：抛物线轨道（刚好逃逸）
• $E > 0$：双曲线轨道（散射）

量子力学的求解将遵循同样的结构：分离变量为角向和径向，角向由 $L^2$ 决定，径向在有效势中运动。

13.0.3 从经典到量子： Bohr模型的启示

Bohr在1913年提出了半经典的氢原子模型：

1. 电子在圆轨道上运动，角动量量子化 $L = n\hbar$
2. 能级 $E_n = -\frac{m_e e^4}{2\hbar^2 n^2} = -\frac{13.6 \text{ eV}}{n^2}$

Bohr模型给出了正确的能级公式，但它是特设的——没有从第一原理推导。Schrödinger方程的求解将自然地导出相同能级，同时给出完整的波函数和量子数解释。理解Bohr模型的历史地位，有助于欣赏Schrödinger解的完备性。

graph TD
    A[中心力场经典力学] --> B["开普勒问题
V(r) ∝ -1/r"]
    B --> C["守恒量: E, L, A"]
    C --> D["Runge-Lenz矢量
轨道闭合"]
    B --> E["有效势
V_eff = L²/2mr² - k/r"]
    E --> F["E<0: 椭圆
束缚态"]
    E --> G["E=0: 抛物线
逃逸极限"]
    E --> H["E>0: 双曲线
散射态"]
    F --> I["Bohr模型
L=nℏ, E_n=-13.6/n² eV"]
    I --> J["Schrödinger方程
严格求解"]
    style C fill:#e3f2fd
    style I fill:#fff8e1

---

13.1 问题的设置与简化

13.1.1 两体问题到单体问题

氢原子由质子（质量 $m_p$）和电子（质量 $m_e$）组成，通过库仑力相互作用：

V(r) = -\frac{e^2}{4\pi\varepsilon_0 r} = -\frac{e^2}{r} \quad \text{(高斯单位制)}

利用第10章的质心坐标变换，两体Hamiltonian分离为：

H = \underbrace{\frac{P^2}{2M}}_{\text{质心运动}} + \underbrace{\left(\frac{p^2}{2\mu} - \frac{e^2}{r}\right)}_{\text{相对运动}}

其中约化质量：

$\mu = \frac{m_e m_p}{m_e + m_p} \approx m_e \left(1 - \frac{m_e}{m_p}\right) \approx 0.999455 \, m_e$

质心部分描述整个氢原子的自由运动（通常取为平面波或静止），相对运动部分描述电子相对于质子的内部结构。

13.1.2 球坐标下的分离变量

相对运动的Hamiltonian：

$H = -\frac{\hbar^2}{2\mu}\nabla^2 - \frac{e^2}{r}$

在球坐标 $(r, \theta, \phi)$ 中，Laplacian算符分解为径向和角向部分：

\nabla^2 = \underbrace{\frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2 \frac{\partial}{\partial r}\right)}_{\text{径向}} - \underbrace{\frac{L^2}{\hbar^2 r^2}}_{\text{角向}}

由于 $[H, L^2] = 0$ 和 $[H, L_z] = 0$（库仑势是球对称的！），我们可以寻找 $H$、$L^2$、$L_z$ 的共同本征态：

$\psi(r, \theta, \phi) = R_{nl}(r) Y_l^m(\theta, \phi)$

分离变量后得到两个方程：

角向方程（已解决，即球谐函数）：

$L^2 Y_l^m = \hbar^2 l(l+1) Y_l^m$

径向方程：

\left[-\frac{\hbar^2}{2\mu}\frac{1}{r^2}\frac{d}{dr}\left(r^2 \frac{d}{dr}\right) + \underbrace{\frac{\hbar^2 l(l+1)}{2\mu r^2}}_{\text{离心势}} - \frac{e^2}{r}\right] R_{nl}(r) = E_{nl} R_{nl}(r)

离心势项 $\frac{\hbar^2 l(l+1)}{2\mu r^2}$ 是角动量的"有效排斥势"——角动量越大的态越难接近原点。

graph TD
    A[氢原子两体问题] --> B[质心坐标变换]
    B --> C["质心运动: 平面波"]
    B --> D["相对运动: 单体问题"]
    D --> E[球坐标分离变量]
    E --> F["ψ = R(r)Y("\\"θ,φ\\"")"]
    F --> G[角向方程]
    F --> H[径向方程]
    G --> I["球谐函数 Y_l^m
已知解"]
    H --> J[含离心势的有效势]
    J --> K["V_eff = ℏ²l("\\"l+1\\"")/2μr² - e²/r"]
    K --> L[求解径向波函数 R_nl]
    style I fill:#e8f5e9
    style K fill:#fff3e0

---

13.2 径向方程的求解

13.2.1 无量纲化与特征尺度

引入Bohr半径（氢原子的特征长度尺度）：

\boxed{a_0 = \frac{\hbar^2}{\mu e^2} \approx 0.529 \times 10^{-10} \, \text{m} \approx 0.529 \, \text{Å}}

和Rydberg能量（特征能量尺度）：

\boxed{E_R = \frac{\mu e^4}{2\hbar^2} = \frac{e^2}{2a_0} \approx 13.6 \, \text{eV}}

无量纲变量：$\rho = r/a_0$，$\epsilon = E/E_R$。径向方程变为：

$\left[\frac{d^2}{d\rho^2} + \frac{2}{\rho}\frac{d}{d\rho} - \frac{l(l+1)}{\rho^2} + \frac{2}{\rho} + \epsilon\right] R(\rho) = 0$

13.2.2 渐近行为分析

在求解之前，先分析两个极限：

大距离极限 ($\rho \to \infty$)：方程近似为

$\frac{d^2 R}{d\rho^2} + \epsilon R \approx 0$

对于束缚态 ($\epsilon < 0$)，解必须指数衰减：$R \sim e^{-\sqrt{-\epsilon}\rho}$。定义 $n = 1/\sqrt{-\epsilon}$（最终将证明 $n$ 是主量子数），则渐近行为是 $e^{-\rho/n}$。

小距离极限 ($\rho \to 0$)：离心势主导

$\frac{d^2 R}{d\rho^2} + \frac{2}{\rho}\frac{dR}{d\rho} - \frac{l(l+1)}{\rho^2} R \approx 0$

设 $R \sim \rho^s$，得到指标方程 $s(s+1) = l(l+1)$，解为 $s = l$ 或 $s = -(l+1)$。后者在原点发散不可归一化，所以：

R(\rho) \xrightarrow{\rho \to 0} \rho^l

13.2.3 级数求解与量子化条件

综合渐近行为，设 $R(\rho) = \rho^l e^{-\rho/n} u(\rho)$，其中 $u(\rho)$ 是待求的多项式。代入径向方程得到 $u(\rho)$ 的方程，这是一个合流超几何方程。

级数展开 $u(\rho) = \sum_{k=0}^{\infty} c_k \rho^k$，递推关系要求级数必须在某阶截断为多项式，否则 $u(\rho)$ 指数增长破坏归一化。截断条件给出：

\boxed{n = n_r + l + 1, \quad n_r = 0, 1, 2, \ldots}

其中 $n_r$ 是径向节点数（波函数在 $r > 0$ 处的零点数）。由于 $n_r \geq 0$，我们得到主量子数：

\boxed{n = 1, 2, 3, \ldots}

而角量子数受限于：

\boxed{l = 0, 1, 2, \ldots, n-1}

13.2.4 能级公式：量子化的胜利

将 $n$ 代回能量表达式 $\epsilon = -1/n^2$：

\boxed{E_n = -\frac{E_R}{n^2} = -\frac{\mu e^4}{2\hbar^2 n^2} = -\frac{13.6 \, \text{eV}}{n^2}}

这是量子力学最伟大的公式之一。让我们逐项解读：

• 负号表示束缚态（将电子移到无穷远需要 $13.6/n^2$ eV的能量）
• $n = 1$（基态）：$E_1 = -13.6$ eV，这就是氢原子的电离能
• $n \to \infty$：$E \to 0$，对应电离极限
• 能级间距随 $n$ 增大而减小：$E_n - E_{n+1} \approx 2E_R/n^3 \to 0$，在高能处形成Rydberg态的密集系列

graph LR
    subgraph 能级图
        A1["电离极限: E=0"]
        B1["n=∞"]
        C1["n=4: -0.85 eV"]
        D1["n=3: -1.51 eV"]
        E1["n=2: -3.4 eV"]
        F1["n=1: -13.6 eV"]
    end
    G[光谱跃迁]
    H["Lyman系
n→1, 紫外"]
    I["Balmer系
n→2, 可见光"]
    J["Paschen系
n→3, 红外"]
    A1 --> B1
    B1 --> C1
    C1 --> D1
    D1 --> E1
    E1 --> F1
    E1 --> H
    D1 --> I
    C1 --> J
    style F1 fill:#ffebee
    style E1 fill:#fff8e1

---

13.3 波函数的完整结构与量子数

13.3.1 三个量子数的物理意义

氢原子的完整波函数由三个量子数标记：

$\psi_{nlm}(r, \theta, \phi) = R_{nl}(r) Y_l^m(\theta, \phi)$

量子数	名称	取值	物理意义
$n$	主量子数	$1, 2, 3, \ldots$	决定能量 $E_n = -E_R/n^2$，大致决定轨道尺度
$l$	角量子数	$0, 1, \ldots, n-1$	决定轨道角动量大小 $L = \sqrt{l(l+1)}\hbar$
$m$	磁量子数	$-l, -l+1, \ldots, l$	决定 $L_z$ 分量 $L_z = m\hbar$

命名约定（光谱学符号）：$l=0,1,2,3,4,\ldots$ 对应 $s, p, d, f, g, \ldots$

13.3.2 径向波函数 $R_{nl}(r)$

径向波函数包含一个指数衰减因子和一个关联Laguerre多项式：

\boxed{R_{nl}(r) = \sqrt{\left(\frac{2}{na_0}\right)^3 \frac{(n-l-1)!}{2n[(n+l)!]}} \, e^{-\rho/n} \left(\frac{2\rho}{n}\right)^l L_{n-l-1}^{2l+1}\left(\frac{2\rho}{n}\right)}

其中 $\rho = r/a_0$，$L_k^\alpha(x)$ 是关联Laguerre多项式。

前几个径向波函数：

• 1s ($n=1, l=0$): $R_{10}(r) = 2 a_0^{-3/2} e^{-r/a_0}$
• 2s ($n=2, l=0$): $R_{20}(r) = \frac{1}{\sqrt{2}} a_0^{-3/2} \left(1 - \frac{r}{2a_0}\right) e^{-r/(2a_0)}$
• 2p ($n=2, l=1$): $R_{21}(r) = \frac{1}{\sqrt{24}} a_0^{-3/2} \frac{r}{a_0} e^{-r/(2a_0)}$

graph TD
    A["径向波函数 R_nl(r)"] --> B[归一化因子]
    A --> C["指数衰减 e^{-r/na₀}"]
    A --> D["r^l 因子"]
    A --> E[关联Laguerre多项式]
    B --> F["∫|R|²r²dr = 1"]
    C --> G["束缚态特征
空间范围~n²a₀"]
    D --> H["原点行为
R ~ r^l"]
    E --> I["径向节点数
n_r = n-l-1"]
    I --> J["1s: 0个节点"]
    I --> K["2s: 1个节点"]
    I --> L["2p: 0个节点"]
    style C fill:#e3f2fd
    style I fill:#fff8e1

13.3.3 概率密度与电子云

径向概率密度（在 $r$ 到 $r+dr$ 球壳内找到电子的概率）：

$P(r) dr = |R_{nl}(r)|^2 r^2 dr$

重要区分：$|R_{nl}|^2$ 是径向概率密度（单位 $r$），而 $|\psi|^2$ 是体积概率密度（单位体积）。乘以 $r^2$ 是因为球壳体积 $4\pi r^2 dr$。

基态1s的最可几半径（概率密度最大处）：

$\frac{dP_{10}}{dr} = 0 \quad \Rightarrow \quad r = a_0$

这恰好是Bohr半径！但注意，Bohr模型中电子在固定轨道上运行，而量子力学中 $a_0$ 只是概率最大的位置——电子可以出现在任何地方（包括 $r > a_0$ 和 $r \to 0$ 处，尽管后者概率密度趋于零）。

---

13.4 简并度：隐藏的对称性

13.4.1 意外的简并

氢原子能级公式 $E_n = -E_R/n^2$ 只依赖于 $n$，与 $l$ 无关。这意味着对于给定的 $n$，所有 $l = 0, 1, \ldots, n-1$ 的态能量完全相同。

例如 $n=2$：

• $2s$ ($l=0$)：球对称，一个径向节点
• $2p$ ($l=1$)：哑铃形，无径向节点

这两种形状截然不同的轨道竟然有完全相同的能量！这不是偶然，而是源于库仑势的特殊对称性。

总简并度（考虑所有 $l$ 和 $m$）：

$g_n = \sum_{l=0}^{n-1} (2l+1) = n^2$

$n^2$ 的简并度是库仑势的特征。对于一般的中心势（例如屏蔽库仑势），简并度会被打破，只剩下 $m$ 的简并（$2l+1$ 重），因为球对称性保留但特殊的"动力学对称性"消失了。

13.4.2 Runge-Lenz矢量：隐藏的守恒量

库仑势 $V(r) \propto 1/r$ 拥有一个额外的守恒量——Runge-Lenz矢量：

$\mathbf{A} = \frac{1}{2\mu e^2}(\mathbf{L} \times \mathbf{p} - \mathbf{p} \times \mathbf{L}) + \frac{\mathbf{r}}{r}$

在经典力学中，Runge-Lenz矢量指向椭圆轨道的长轴方向，其守恒性说明轨道形状和取向不变（开普勒问题没有进动）。在量子力学中，$\mathbf{A}$ 与 $H$ 对易，生成了一个更大的对称群——$SO(4)$（四维旋转群），而非仅仅是三维旋转群 $SO(3)$。

$SO(4)$ 的不可约表示维度正是 $n^2$，解释了完整的简并结构。Shankar在这里展示了量子力学与高等数学的美妙连接：氢原子的"意外"简并实际上是更深层次对称性的自然表现。

graph TD
    A[氢原子能级简并] --> B["球对称SO(3)"]
    B --> C["l简并: 意外"]
    B --> D["m简并: 2l+1重
几何必然"]
    C --> E[库仑势特殊对称性]
    E --> F["Runge-Lenz矢量 A"]
    F --> G["A守恒: [A,H]=0"]
    G --> H["对称群扩大
SO(3) → SO(4)"]
    H --> I["SO(4)不可约表示
维度 = n²"]
    I --> J["总简并度 n²
1s:1, 2s+2p:4,
3s+3p+3d:9..."]
    E --> K["一般中心势
无此简并"]
    style F fill:#fff8e1
    style H fill:#e3f2fd
    style J fill:#f3e5f5

---

13.5 光谱与选择定则

13.5.1 光谱线系

电子从高能级 $n_i$ 跃迁到低能级 $n_f$ 时发射光子，能量：

$\hbar \omega = E_{n_i} - E_{n_f} = E_R \left(\frac{1}{n_f^2} - \frac{1}{n_i^2}\right)$

经典光谱线系：

• Lyman系 ($n_f = 1$)：紫外区，$\lambda = 121.6$ nm ($2\to 1$) 等
• Balmer系 ($n_f = 2$)：可见光区，$\lambda = 656.3$ nm ($3\to 2$, H$\alpha$线) 等
• Paschen系 ($n_f = 3$)：红外区

13.5.2 数值例子：氢原子光谱的精确波长计算

Lyman-α线（2→1跃迁）：

$\Delta E = E_2 - E_1 = -13.6 \left(\frac{1}{4} - 1\right) = 10.2 \text{ eV}$

波长：

\lambda = \frac{hc}{\Delta E} = \frac{4.136 \times 10^{-15} \text{ eV·s} \times 2.998 \times 10^8 \text{ m/s}}{10.2 \text{ eV}} \lambda = \frac{1.2398 \times 10^{-6} \text{ eV·m}}{10.2 \text{ eV}} \approx 1.216 \times 10^{-7} \text{ m} = 121.6 \text{ nm}

确实在紫外区。

Balmer-α线/Hα线（3→2跃迁）：

$\Delta E = 13.6 \left(\frac{1}{4} - \frac{1}{9}\right) = 13.6 \times \frac{5}{36} \approx 1.89 \text{ eV}$

$\lambda = \frac{1.2398 \times 10^{-6}}{1.89} \approx 6.56 \times 10^{-7} \text{ m} = 656 \text{ nm}$

这是可见红光！Hα线是宇宙中最丰富的光谱线之一，恒星和星云的发光主要来自此跃迁。

Balmer-β线/Hβ线（4→2跃迁）：

$\Delta E = 13.6 \left(\frac{1}{4} - \frac{1}{16}\right) = 13.6 \times \frac{3}{16} \approx 2.55 \text{ eV}$

$\lambda = \frac{1.2398 \times 10^{-6}}{2.55} \approx 4.86 \times 10^{-7} \text{ m} = 486 \text{ nm}$

这是蓝绿色光，肉眼可见。

13.5.3 电偶极跃迁的选择定则

并非所有能级之间的跃迁都是允许的。在电偶极近似下（光子波长 $\gg$ 原子尺度），跃迁振幅正比于矩阵元 $\langle n' l' m' | \mathbf{r} | n l m \rangle$。利用球谐函数的性质，可以证明非零矩阵元要求：

\boxed{\Delta l = \pm 1, \quad \Delta m = 0, \pm 1}

这是角动量守恒的直接结果：光子携带1单位角动量（$\hbar$），所以原子角动量必须改变 $\pm 1$。$\Delta m = 0$ 对应光子的线偏振（沿z轴），$\Delta m = \pm 1$ 对应圆偏振。

选择定则解释了为什么某些谱线缺失。例如 $2s \to 1s$ 是禁戒的（$\Delta l = 0$），实际上这个跃迁通过双光子发射进行（寿命约0.12秒，远长于普通的电偶极跃迁）。

graph TD
    A[电偶极跃迁] --> B["矩阵元 ⟨n'l'm'|r|nlm⟩"]
    B --> C{"是否非零?"}
    C -->|需要| D["宇称变化: (-1)^l → (-1)^{l+1}"]
    C -->|需要| E["角动量耦合: l ⊗ 1 = l±1"]
    D --> F["Δl = ±1"]
    E --> G["Δm = 0, ±1"]
    F --> H[选择定则]
    G --> H
    H --> I[允许跃迁]
    H --> J[禁戒跃迁]
    J --> K["2s→1s 双光子"]
    J --> L["高阶过程
磁偶极/电四极"]
    style H fill:#e8f5e9
    style J fill:#ffebee

---

13.6 超越Schrödinger：Lamb位移与精细结构

13.6.1 精细结构：相对论修正

Schrödinger方程是非相对论性的。Dirac方程（包含相对论效应）给出更精确的能级公式：

$E_{n,j} = \frac{\mu c^2}{\sqrt{1 + \frac{\alpha^2}{(n - \delta_j)^2}}} - \mu c^2$

其中 $\alpha = e^2/(\hbar c) \approx 1/137$ 是精细结构常数，$\delta_j = j + 1/2 - \sqrt{(j+1/2)^2 - \alpha^2}$ 是量子亏损。

展开到 $\alpha^4$ 阶：

$E_{n,j} \approx -\frac{E_R}{n^2}\left[1 + \frac{\alpha^2}{n}\left(\frac{1}{j+1/2} - \frac{3}{4n}\right)\right]$

关键修正：

• 动能相对论修正：$p^4$ 项，$\sim \alpha^4$
• 自旋-轨道耦合：$\mathbf{L} \cdot \mathbf{S}$ 项，$\sim \alpha^4$
• Darwin项：对s态的接触修正，$\sim \alpha^4$

这些修正将 $l$ 简并打破（因为 $j = l \pm 1/2$ 有不同的能量），但 $2s_{1/2}$ 和 $2p_{1/2}$ 在Dirac理论中仍然是简并的。

13.6.2 Lamb位移：QED的胜利

1947年，Lamb和Retherford用微波光谱技术精确测量了 $2s_{1/2}$ 和 $2p_{1/2}$ 的能级差：

$\Delta E_{\text{Lamb}} \approx 1058 \, \text{MHz} \approx 4.38 \times 10^{-6} \, \text{eV}$

这是一个极小的能量差（约 $10^{-6}$ 倍基态能量），但它在理论上意义巨大——它证明了Dirac方程还不够。这个分裂来源于量子电动力学效应：

1. 真空涨落：电子与"虚光子云"的持续相互作用，使电子的有效电荷分布略有展宽。s态电子更靠近原子核（$|\psi_{2s}(0)|^2 > |\psi_{2p}(0)|^2 = 0$），所以感受的真空涨落更强。

2. 电子自能：电子不断发射和吸收虚光子，获得额外的"电磁质量"。

Feynman、Schwinger和Tomonaga独立发展了计算这些效应的系统方法，分享了1965年诺贝尔物理学奖。Lamb位移是量子场论最成功的定量验证之一。

13.6.3 数值例子：Lamb位移的定量估算

Lamb位移的大小可以用Bethe的半定量公式估算：

$\Delta E_{\text{Lamb}}(2s_{1/2}) \sim \frac{\alpha^3}{\pi} \ln\left(\frac{1}{\alpha}\right) \cdot \frac{E_R}{n^3}$

对于 $n=2$：

$\Delta E_{\text{Lamb}} \sim \frac{(1/137)^3}{\pi} \ln(137) \cdot \frac{13.6}{8} \text{ eV}$

$\sim \frac{3.9 \times 10^{-7}}{3.14} \times 4.91 \times 1.7 \text{ eV}$

$\sim 1 \times 10^{-6} \text{ eV} \sim 300 \text{ MHz}$

（精确值为1058 MHz，这个估算在数量级上是正确的）

转换为波长：

$\lambda = \frac{c}{\nu} = \frac{3 \times 10^8}{1.058 \times 10^9} \approx 0.284 \text{ m} = 28.4 \text{ cm}$

这是微波波段！Lamb和Retherford正是用微波腔在28 cm附近扫描，发现了这个微小的能级移动。

13.6.4 超精细结构：核自旋的印记

质子和电子都有自旋 $1/2$，它们的磁矩相互作用产生超精细分裂：

$\Delta E_{\text{hfs}} \sim \frac{\mu_p \mu_e}{a_0^3} \sim 5.9 \times 10^{-6} \, \text{eV} \approx 1420 \, \text{MHz}$

著名的 21厘米线（波长21.1 cm，频率1420 MHz）就是氢原子基态超精细跃迁（自旋平行 $F=1$ 到反平行 $F=0$）发出的射电波。它是射电天文学中最普遍的观测谱线，也是SETI（搜寻地外文明）项目的标志性频率。

---

本章总结

氢原子是量子力学从数学结构走向物理现实的转折点。通过质心坐标变换，两体问题被优雅地简化为单体在库仑势中的运动。球坐标下的分离变量将问题拆分为已知的角向部分（球谐函数）和待解的径向部分，后者通过渐近分析和级数求解导出了量子化条件。

能级公式 $E_n = -13.6 \text{ eV}/n^2$ 是量子力学最辉煌的成就之一。它不仅精确解释了氢原子光谱的所有线系（Lyman、Balmer、Paschen等），更揭示了自然界量子化的深层结构。主量子数 $n$ 决定能量，角量子数 $l$ 决定轨道形状，磁量子数 $m$ 决定空间取向——三个量子数共同编织出原子世界的完整图景。

简并度 $n^2$ 背后隐藏着Runge-Lenz矢量守恒和 $SO(4)$ 动力学对称性，这是数学之美在物理中的体现。选择定则 $\Delta l = \pm 1$ 则是角动量守恒的直接推论，解释了为什么某些跃迁被禁戒而另一些则发射明亮的光谱线。

但Schrödinger的解只是开始。精细结构修正打破了 $l$ 简并，Lamb位移进一步揭示了QED的真空涨落，超精细结构则记录了原子核自旋的微弱印记。从13.6 eV到 $10^{-6}$ eV，跨越九个数量级的能级修正，展示了量子力学作为框架理论的包容性和开放性。

graph TD
    A[氢原子求解] --> B["两体→单体
质心系"]
    B --> C[球坐标分离变量]
    C --> D["角向: 球谐函数
Y_l^m"]
    C --> E["径向: 级数求解
关联Laguerre多项式"]
    E --> F["量子化条件
n = n_r + l + 1"]
    F --> G["能级公式
E_n = -E_R/n²"]
    G --> H["光谱线系
Lyman/Balmer/Paschen"]
    H --> I["选择定则
Δl=±1, Δm=0,±1"]
    G --> J["简并度 n²
Runge-Lenz矢量
SO(4)对称性"]
    G --> K["精细结构
Dirac方程"]
    K --> L["Lamb位移
QED真空涨落"]
    L --> M["超精细结构
核自旋: 21cm线"]
    style G fill:#e8f5e9
    style J fill:#e3f2fd
    style L fill:#fff8e1

---

练习与思考

1. Bohr半径的物理意义：证明氢原子基态（1s）中电子位置的最可几距离恰好等于Bohr半径 $a_0$。然后计算 $\langle r \rangle_{1s}$（期望值）并与最可几值比较。为什么这两个值不同？

2. 选择定则与对称性：利用宇称性质解释为什么 $2s \to 1s$ 的电偶极跃迁是禁戒的。提示：电偶极算符 $\mathbf{r}$ 是奇宇称（$\Pi \mathbf{r} \Pi = -\mathbf{r}$），而 $2s$ 和 $1s$ 都是偶宇称（$l=0$）。因此矩阵元 $\langle 1s | \mathbf{r} | 2s \rangle$ 的宇称如何？

3. Rydberg原子：考虑一个主量子数 $n = 100$ 的高激发态氢原子（Rydberg原子）。估算：(a) 这个态的轨道半径（用 $a_0$ 表示）；(b) 相邻能级间距 $E_{100} - E_{101}$（用eV表示）；© 将这个能量差转换为频率，判断它处于电磁波谱的哪个区域。讨论为什么Rydberg原子对电场扰动极其敏感。

4. 光谱线的数值验证：利用能级公式 $E_n = -13.6/n^2$ eV，计算 $4 \to 2$ 跃迁（Hβ线）的波长，并与实验值486.1 nm比较。若使用约化质量（$E_R = 13.598$ eV），波长偏移多少皮米？