Appearance
氢原子
一、问题的提出
1.1 氢原子的重要性
氢原子是自然界中最简单的原子——它仅由一个质子和一个电子组成，在量子力学中恰好构成一个精确可解的二体问题。正因为这种简洁性，氢原子成为了量子力学最重要的试金石之一。任何量子力学理论的正确性，首先需要经受氢原子光谱的检验。
历史地位
氢原子光谱的巴尔末公式（1885年）是通向现代量子力学的第一块铺路石。玻尔正是在解释氢原子光谱的尝试中提出了他的原子模型，而薛定谔方程的建立也以氢原子为第一个完整的应用案例。
氢原子之所以重要，可以归纳为以下几个核心原因：
- 精确可解性：它是少数几个可以精确求解薛定谔方程的三维量子系统之一，为理解量子力学的基本原理提供了一个完美的范例。
- 理论检验：氢原子光谱的精细结构、超精细结构、塞曼效应和斯塔克效应等，分别检验了量子力学中相对论修正、自旋理论、电磁相互作用等各个层面的正确性。
- 概念基础：氢原子的求解过程自然地引入了主量子数、角量子数、磁量子数以及自旋量子数等核心概念，这些概念构成了整个原子物理学的语言基础。
- 多电子原子的出发点：理解氢原子是理解多电子原子的前提。类氢离子的波函数常被用作多电子原子波函数的零级近似。
1.2 薛定谔方程
氢原子中，电子在质子产生的库仑势场中运动。将质子和电子视为二体问题，经过质心坐标变换后，问题简化为一个约化质量为 $\mu = \frac{m_e m_p}{m_e + m_p} \approx m_e$ 的粒子在固定中心力场中的运动。
库仑势的形式为：
$$ V(r) = -\frac{e^2}{4\pi\varepsilon_0 r} $$
其中 $e$ 是基本电荷，$\varepsilon_0$ 是真空介电常数，$r$ 是电子到核的距离。负号表示吸引力。
三维定态薛定谔方程写作：
$$ -\frac{\hbar^2}{2m}\nabla^2\psi(\mathbf{r}) - \frac{e^2}{4\pi\varepsilon_0 r}\psi(\mathbf{r}) = E\psi(\mathbf{r}) $$
符号约定
在量子力学文献中，有时会使用高斯单位制，此时库仑势写作 $V(r) = -e^2/r$。本文采用国际单位制（SI），初学者在阅读不同教材时需要注意单位的转换。
这个方程由于势能 $V(r)$仅依赖于径向距离$r$，具有球对称性，因此我们自然选择球坐标 $(r, \theta, \phi)$ 来求解。球对称性意味着系统的哈密顿量与角动量算符对易，角动量是守恒量，这为分离变量提供了物理基础。
二、分离变量
2.1 球坐标下的薛定谔方程
在球坐标中，拉普拉斯算符 $\nabla^2$ 的形式为：
$$ \nabla^2 = \frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial}{\partial r}\right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{r^2\sin^2\theta}\frac{\partial^2}{\partial\phi^2} $$
将此代入薛定谔方程，并利用角动量平方算符在球坐标中的表达式：
$$ \hat{L}^2 = -\hbar^2\left[\frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\phi^2}\right] $$
我们可以将哈密顿量重写为：
$$ \hat{H} = -\frac{\hbar^2}{2m}\frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial}{\partial r}\right) + \frac{\hat{L}^2}{2mr^2} - \frac{e^2}{4\pi\varepsilon_0 r} $$
由于 $\hat{L}^2$仅与$\theta$和$\phi$有关，而势能仅与$r$ 有关，我们设波函数为径向部分和角向部分的乘积：
$$ \psi(r,\theta,\phi) = R(r)\thinspace{}Y(\theta,\phi) $$
2.2 角向方程
将分离变量形式的波函数代入薛定谔方程后，角向部分满足本征方程：
$$ \hat{L}^2Y(\theta,\phi) = \hbar^2\lambda\thinspace{}Y(\theta,\phi) $$
这就是轨道角动量的本征方程。进一步分离变量 $Y(\theta,\phi) = \Theta(\theta)\Phi(\phi)$，可以解得：
- $\Phi(\phi) = \frac{1}{\sqrt{2\pi}}e^{im\phi}$，其中 $m = 0, \pm 1, \pm 2, \ldots$
- $\Theta(\theta)$由连带勒让德方程描述，仅当$\lambda = l(l+1)$且$|m| \leq l$ 时才有物理上可接受的解
最终得到的角向函数是球谐函数 $Y_{lm}(\theta,\phi)$：
$$ Y_{lm}(\theta,\phi) = (-1)^m\sqrt{\frac{(2l+1)}{4\pi}\frac{(l-m)!}{(l+m)!}}\thinspace{}P_l^m(\cos\theta)\thinspace{}e^{im\phi} $$
其中 $P_l^m$ 是连带勒让德多项式。球谐函数满足正交归一性：
$$ \int_0^{2\pi}\int_0^\pi Y_{l'm'}^*(\theta,\phi)Y_{lm}(\theta,\phi)\sin\theta\thinspace{}d\theta\thinspace{}d\phi = \delta_{ll'}\delta_{mm'} $$
与角动量的关联
这部分内容与角动量与自旋密切相关。球谐函数 $Y_{lm}$是轨道角动量算符$\hat{L}^2$和$\hat{L}_z$ 的共同本征函数： $$ \hat{L}^2 Y_{lm} = \hbar^2 l(l+1) Y_{lm}, \quad \hat{L}_z Y_{lm} = \hbar m Y_{lm} $$ 其中 $l$ 称为角量子数，$m$ 称为磁量子数。$l$ 必须是非负整数。
2.3 径向方程
将角向方程的结果代入薛定谔方程，径向部分 $R(r)$ 满足：
$$ -\frac{\hbar^2}{2m}\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{dR}{dr}\right) + \left[\frac{\hbar^2 l(l+1)}{2mr^2} - \frac{e^2}{4\pi\varepsilon_0 r}\right]R(r) = ER(r) $$
引入 $u(r) = rR(r)$，方程简化为一个一维形式：
$$ -\frac{\hbar^2}{2m}\frac{d^2u}{dr^2} + V_{\text{eff}}(r)\thinspace{}u(r) = E\thinspace{}u(r) $$
其中有效势为：
$$ V_{\text{eff}}(r) = -\frac{e^2}{4\pi\varepsilon_0 r} + \frac{\hbar^2 l(l+1)}{2mr^2} $$
离心势垒的物理意义
有效势中的第二项 $\frac{\hbar^2 l(l+1)}{2mr^2}$称为离心势垒。它来源于角动量，表现为一个排斥性的势。经典力学中，角动量不为零的粒子无法到达$r=0$ 处，因为离心力会将其推开。量子力学中，$l=0$ 的 s 态电子不受离心势垒的影响，可以出现在核的位置，这对于理解如核磁共振等现象具有重要意义。
对于 $l > 0$的态，有效势在$r_0 = \frac{4\pi\varepsilon_0\hbar^2 l(l+1)}{me^2}$ 处有一个极小值，这对应于经典力学中稳定的圆轨道。但这个经典图像的类比仅具有启发意义，量子力学中电子的分布由概率密度描述。
三、径向方程的求解
3.1 无量纲化
为了简化方程，我们引入无量纲变量。首先定义玻尔半径：
$$ a_0 = \frac{4\pi\varepsilon_0\hbar^2}{me^2} \approx 0.529\thinspace\text{Å} $$
然后引入无量纲径向变量 $\rho = r/a_0$。同时定义无量纲能量参数 $\kappa$：
$$ \kappa = \sqrt{-\frac{2mE}{\hbar^2}}\thinspace{}a_0 $$
对于束缚态（$E < 0$），$\kappa$ 是正实数。径向方程变为：
$$ \frac{d^2u}{d\rho^2} = \left[1 - \frac{2}{\kappa\rho} + \frac{l(l+1)}{\rho^2}\right]u $$
3.2 渐近行为
分析径向方程在 $\rho \to \infty$和$\rho \to 0$ 两个极限下的行为，对于理解解的结构至关重要。
$\rho \to \infty$的渐近行为：当$\rho$很大时，方程中$1/\rho$和$1/\rho^2$ 项可以忽略，得：
$$ \frac{d^2u}{d\rho^2} \approx u \quad\Rightarrow\quad u(\rho) \sim e^{-\rho} $$
考虑到归一化（波函数必须在无穷远处趋于零），我们取衰减解 $u \sim e^{-\rho}$。更一般地，最终解会显示 $u \sim e^{-\rho/n}$，其中 $n$ 是主量子数。
$\rho \to 0$的渐近行为：当$\rho$很小时，离心势垒项$l(l+1)/\rho^2$ 占主导，方程近似为：
$$ \frac{d^2u}{d\rho^2} \approx \frac{l(l+1)}{\rho^2}u $$
设 $u \sim \rho^s$，代入得 $s(s-1) = l(l+1)$，解得 $s = l+1$或$s = -l$。由于波函数在原点必须有限（$u(0) = 0$），我们取 $s = l+1$，即：
$$ u(\rho) \sim \rho^{l+1}, \quad R(\rho) \sim \rho^l \quad (\rho \to 0) $$
物理直觉
$l$越大，电子在核附近出现的概率越小。这是因为离心势垒将电子"推开"了核。对于$l=0$ 的 s 态，$R(0) \neq 0$，电子在核处有非零的概率密度，这一事实在微扰理论中处理费米接触相互作用时至关重要。
3.3 级数解法
将渐近行为提取出来，设径向函数的形式为：
$$ u(\rho) = \rho^{l+1}e^{-\rho/n}\thinspace{}F(\rho) $$
将这个形式代入径向方程，$F(\rho)$ 满足的微分方程可以化为合流超几何方程（Kummer 方程）：
$$ \rho\frac{d^2F}{d\rho^2} + (2l+2-\rho)\frac{dF}{d\rho} + \left(\frac{1}{\kappa}-l-1\right)F = 0 $$
设 $F(\rho)$为幂级数$F(\rho) = \sum_{k=0}^{\infty} c_k\rho^k$，代入后得到递推关系：
$$ c_{k+1} = \frac{k + l + 1 - 1/\kappa}{(k+1)(k+2l+2)}\thinspace{}c_k $$
关键的量子化条件：为了使波函数在 $\rho \to \infty$时可归一化，级数必须在某有限项处截断，即$F(\rho)$必须是一个多项式。这意味着存在某个整数$k_{\max}$使得$c_{k_{\max}+1} = 0$，从而：
$$ k_{\max} + l + 1 = \frac{1}{\kappa} $$
令 $n = k_{\max} + l + 1$，则 $n$是一个正整数，且$n \geq l+1$。由此我们得到：
$$ \kappa = \frac{1}{n} $$
这就是量子化条件的来源——束缚态的要求迫使级数截断，从而能量只能取离散值。
3.4 能级公式
由 $\kappa = 1/n$和$\kappa$ 的定义，立即得到氢原子的能级公式：
$$ E_n = -\frac{me^4}{8\varepsilon_0^2 h^2}\frac{1}{n^2} = -\frac{13.6\thinspace\text{eV}}{n^2} $$
其中 $n = 1, 2, 3, \ldots$ 称为主量子数。
量子数的取值范围
从级数截断条件 $n = k_{\max} + l + 1$可知，对于给定的$n$：
- 角量子数 $l = 0, 1, 2, \ldots, n-1$（共 $n$ 个取值）
- 磁量子数 $m = -l, -l+1, \ldots, l$（共 $2l+1$ 个取值）
因此第 $n$ 个能级的简并度为： $$ \sum_{l=0}^{n-1}(2l+1) = n^2 $$ 这个 $n^2$ 简并度是库仑势特有的性质，来源于库仑势的动力学对称性——拉普拉斯-龙格-楞次矢量守恒。
3.5 径向波函数
径向波函数 $R_{nl}(r)$ 的具体形式为：
$$ R_{nl}(r) = \sqrt{\left(\frac{2}{na_0}\right)^3\frac{(n-l-1)!}{2n[(n+l)!]}}\thinspace\left(\frac{2r}{na_0}\right)^l\thinspace{}e^{-r/(na_0)}\thinspace{}L_{n-l-1}^{2l+1}\negthinspace\left(\frac{2r}{na_0}\right) $$
其中 $L_{q}^{p}(x)$ 是缔合拉盖尔多项式（associated Laguerre polynomials）。径向波函数满足归一化条件：
$$ \int_0^\infty |R_{nl}(r)|^2\thinspace{}r^2 dr = 1 $$
前几个径向波函数：
- $R_{10}(r) = 2a_0^{-3/2}\thinspace{}e^{-r/a_0}$ —— 基态（1s）
- $R_{20}(r) = \frac{1}{\sqrt{2}}a_0^{-3/2}\left(1 - \frac{r}{2a_0}\right)e^{-r/(2a_0)}$ —— 2s 态
- $R_{21}(r) = \frac{1}{2\sqrt{6}}a_0^{-3/2}\thinspace\frac{r}{a_0}\thinspace{}e^{-r/(2a_0)}$ —— 2p 态
完整的氢原子波函数为：
$$ \psi_{nlm}(r,\theta,\phi) = R_{nl}(r)\thinspace{}Y_{lm}(\theta,\phi) $$
这构成了量子力学基本假设中态矢量的坐标表示，也是波函数与薛定谔方程的一个经典应用实例。
四、量子数与原子结构
4.1 三个量子数
氢原子的量子态由三个量子数完全确定：
| 量子数 | 符号 | 取值范围 | 物理意义 |
|---|---|---|---|
| 主量子数 | $n$ | $1, 2, 3, \ldots$ | 决定能量$E_n = -13.6/n^2$ eV |
| 角量子数 | $l$ | $0, 1, \ldots, n-1$ | 决定轨道角动量$\Vert\mathbf{L}\Vert = \hbar\sqrt{l(l+1)}$ |
| 磁量子数 | $m$ | $-l, \ldots, +l$ | 决定角动量$z$分量$L_z = m\hbar$ |
常见误解
玻尔模型中角动量为 $n\hbar$，但量子力学给出的是 $\hbar\sqrt{l(l+1)}$。对于基态（$n=1, l=0$），量子力学预言角动量为零，而玻尔模型预言为 $\hbar$。实验（如基态氢原子没有轨道磁矩）支持量子力学的结论。
这三个量子数自然地出现在求解薛定谔方程的过程中，它们分别对应于径向方程（$n$）、角向方程（$l$和$m$）的可解条件。正如量子力学发展史所展示的，这些量子数的发现经历了一个从经验规则（巴尔末公式、里德伯公式）到理论推导（薛定谔方程）的漫长过程。
4.2 电子云
量子力学中，电子不再沿着经典轨道运动，而是以概率云的形式分布在核周围。概率密度为：
$$ |\psi_{nlm}(r,\theta,\phi)|^2 = |R_{nl}(r)|^2\thinspace|Y_{lm}(\theta,\phi)|^2 $$
径向概率分布：在 $r$到$r+dr$ 的球壳内找到电子的概率为：
$$ P_{nl}(r)\thinspace{}dr = r^2|R_{nl}(r)|^2\thinspace{}dr $$
$r^2|R_{nl}(r)|^2$称为径向概率密度。对于基态$n=1, l=0$：
$$ P_{10}(r) = \frac{4}{a_0^3}r^2 e^{-2r/a_0} $$
该函数在 $r = a_0$（玻尔半径）处取最大值，与玻尔模型中的基态轨道半径一致。但这只是最大概率位置，而非唯一的轨道半径——电子有一定的概率出现在任何 $r$ 处。
角向概率分布：$|Y_{lm}(\theta,\phi)|^2$ 描述了电子概率分布的角度依赖性。例如：
- $Y_{00}$ 是球形对称的（s 轨道）
- $Y_{10} \propto \cos\theta$ 呈哑铃形（p 轨道）
- $Y_{20} \propto (3\cos^2\theta - 1)$ 呈四瓣形（d 轨道）
轨道形状的直观理解
当我们在化学课本中看到 s 轨道是球形、p 轨道是哑铃形、d 轨道是四瓣梅花形时，这些形状其实就是 $|Y_{lm}(\theta,\phi)|^2$ 在三维空间中的极坐标图（polar plot）。这些形状直接决定了化学键的方向性和分子的几何构型。
4.3 原子轨道
不同 $(n,l)$ 组合对应的轨道及其光谱学记号：
| $n$ | $l=0$(s) | $l=1$(p) | $l=2$(d) | $l=3$ (f) |
|---|---|---|---|---|
| 1 | 1s | — | — | — |
| 2 | 2s | 2p | — | — |
| 3 | 3s | 3p | 3d | — |
| 4 | 4s | 4p | 4d | 4f |
光谱学记号的历史来源：s = sharp（锐线系），p = principal（主线系），d = diffuse（漫线系），f = fundamental（基线系）。字母 g, h, i... 在 $l=4,5,6...$ 时按字母顺序继续使用。
每个轨道都可以容纳两个自旋相反的电子（泡利不相容原理），这正是元素周期表结构的量子力学基础。
五、玻尔模型的对应
5.1 玻尔半径
玻尔半径 $a_0$ 是原子尺度的自然单位：
$$ a_0 = \frac{4\pi\varepsilon_0\hbar^2}{me^2} \approx 0.529\thinspace\text{Å} = 5.29 \times 10^{-11}\thinspace\text{m} $$
$a_0$ 是氢原子基态电子出现概率最大的径向距离。有趣的是，量子力学中玻尔半径作为一个自然长度尺度出现，但并非像玻尔模型那样被假定为固定的轨道半径，而是从理论和基本常数中自然推导出来的。
历史视角
玻尔在1913年通过量子化角动量 $L = n\hbar$的先验假设得到了$a_0$。而薛定谔在1926年从波动方程出发，无需任何人为量子化条件，自然地导出了 $a_0$ 作为氢原子问题的特征长度。这是量子力学发展史中的一个重要转折点。
5.2 能级公式的一致性
量子力学给出的能级公式 $E_n = -13.6\thinspace\text{eV}/n^2$ 与玻尔模型完全一致。这并非巧合，而是因为库仑势具有特殊的动力学对称性。
然而，两者在角动量预言上有根本差异：
| 物理量 | 玻尔模型 | 量子力学 |
|---|---|---|
| 角动量 | $L = n\hbar$ | $L = \hbar\sqrt{l(l+1)}$ |
| 基态角动量 | $L = \hbar$ | $L = 0$ |
| 允许的 $l$值 | $l = n$（隐含） | $l = 0, 1, \ldots, n-1$ |
实验证据明确支持量子力学：基态氢原子没有轨道磁矩（斯特恩-盖拉赫实验未观测到基态氢原子的轨道磁矩偏转），这与 $L=0$ 的预言一致。
5.3 电子在核的概率
一个玻尔模型完全无法解释的现象是：s 态电子（$l=0$）在核处有非零的概率密度。
对于 1s 态： $$ |\psi_{100}(0)|^2 = \frac{1}{\pi a_0^3} $$
电子在核处的非零概率密度导致了一些重要的物理效应，包括：
- 费米接触相互作用：电子与核磁矩的接触超精细相互作用
- 电子俘获：核内质子俘获 K 层电子的过程
- 正电子素（positronium）的湮灭：电子和正电子在 $r=0$ 处重合时发生湮灭
六、精细结构
6.1 相对论修正
前面的讨论基于非相对论薛定谔方程。然而，实际的氢原子光谱线在更高分辨率下会显示出分裂——这就是精细结构。精细结构主要来源于三个相对论效应：
（1）动能的相对论修正：电子的动能应为 $T = \sqrt{p^2c^2 + m^2c^4} - mc^2$，展开到 $p^4$ 阶：
$$ \hat{H}_{\text{rel}} = -\frac{\hat{p}^4}{8m^3c^2} $$
这个修正项的数量级约为 $\alpha^2 E_n \approx 10^{-4}E_n$，其中 $\alpha$ 是精细结构常数。
（2）达尔文项（Darwin term）：仅对 $l=0$ 的 s 态有贡献，来源于电子的颤动运动（Zitterbewegung）：
$$ \hat{H}_{\text{Darwin}} = \frac{\hbar^2}{8m^2c^2}\nabla^2 V = \frac{\hbar^2 e^2}{8m^2c^2\varepsilon_0}\cdot 4\pi\delta^3(\mathbf{r}) $$
（3）自旋-轨道耦合：见下节。
6.2 自旋-轨道耦合
电子具有内禀自旋角动量 $\mathbf{S}$，其自旋磁矩为 $\boldsymbol{\mu}_s = -\frac{e}{m}\mathbf{S}$。在电子静止系中，核绕电子运动产生磁场，该磁场与电子自旋磁矩相互作用产生自旋-轨道耦合能：
$$ \hat{H}_{\text{SO}} = \frac{1}{2m^2c^2}\frac{1}{r}\frac{dV}{dr}\thinspace\hat{\mathbf{L}}\cdot\hat{\mathbf{S}} $$
对于库仑势 $V(r) = -e^2/(4\pi\varepsilon_0 r)$：
$$ \hat{H}_{\text{SO}} = \frac{e^2}{8\pi\varepsilon_0 m^2c^2}\frac{1}{r^3}\thinspace\hat{\mathbf{L}}\cdot\hat{\mathbf{S}} $$
总角动量
当存在自旋-轨道耦合时，$\hat{L}_z$和$\hat{S}_z$不再分别与哈密顿量对易，但总角动量$\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$仍然是守恒量。因此我们使用量子数$j$（总角动量量子数）和 $m_j$ 来标记态。这与角动量与自旋中讨论的角动量耦合理论完全一致。
$\hat{\mathbf{L}}\cdot\hat{\mathbf{S}}$ 的本征值为： $$ \langle\hat{\mathbf{L}}\cdot\hat{\mathbf{S}}\rangle = \frac{\hbar^2}{2}[j(j+1) - l(l+1) - s(s+1)] $$
对于单电子，$s = 1/2$，$j = l \pm 1/2$（当 $l > 0$）。
6.3 精细结构常数
精细结构常数 $\alpha$ 是物理学中最重要的无量纲常数之一：
$$ \alpha = \frac{e^2}{4\pi\varepsilon_0\hbar c} \approx \frac{1}{137.036} $$
$\alpha$ 表征了电磁相互作用的强度。精细结构常数的小数值（$\alpha \ll 1$）意味着量子电动力学（QED）的微扰展开是可行的，这也是为什么我们可以用微扰理论来处理精细结构修正。
6.4 精细结构能级
综合考虑三种相对论修正，氢原子能级的精细结构公式为（狄拉克方程的精确解）：
$$ E_{nj} = m c^2\left[\left(1 + \left(\frac{\alpha}{n - (j+1/2) + \sqrt{(j+1/2)^2 - \alpha^2}}\right)^2\right)^{-1/2} - 1\right] $$
展开到 $\alpha^4$ 阶：
$$ E_{nj} \approx -\frac{13.6\thinspace\text{eV}}{n^2}\left[1 + \frac{\alpha^2}{n}\left(\frac{1}{j+1/2} - \frac{3}{4n}\right)\right] $$
对于 $n=2$ 能级：
- $2S_{1/2}$（$l=0, j=1/2$）和 $2P_{1/2}$（$l=1, j=1/2$）具有相同的能量（狄拉克方程预言）
- $2P_{3/2}$（$l=1, j=3/2$）能量略高
$2P_{3/2}$与$2P_{1/2}$之间的能量差对应氢原子$\mathrm{H}_\alpha$线（656.3 nm）的双重结构，分裂间距约为$4.53 \times 10^{-5}\thinspace\text{eV}$。
兰姆位移
1947年，兰姆和雷瑟福实验发现 $2S_{1/2}$和$2P_{1/2}$ 并非完全简并，两者之间存在约 1057 MHz 的能量差，这就是兰姆位移。兰姆位移是量子电动力学（QED）的重要实验验证，来源于电子与真空涨落的相互作用——这是狄拉克方程无法预言的纯量子场论效应。
七、超精细结构
7.1 核自旋与电子磁矩的相互作用
原子核具有内禀自旋 $\mathbf{I}$和核磁矩$\boldsymbol{\mu}_I = g_I\frac{e}{2m_p}\mathbf{I}$。核磁矩与电子产生的磁场之间的相互作用导致能级的进一步分裂——超精细结构。
超精细相互作用的哈密顿量为：
$$ \hat{H}_{\text{HFS}} = \frac{\mu_0}{4\pi}g_I\frac{e}{2m_p}\frac{e}{m_e}\frac{\hat{\mathbf{I}}\cdot\hat{\mathbf{S}}}{r^3} $$
对于 s 态电子，还需要加上费米接触项：
$$ \hat{H}_{\text{contact}} = \frac{2\mu_0}{3}g_I\frac{e}{2m_p}\frac{e}{m_e}|\psi(0)|^2\thinspace\hat{\mathbf{I}}\cdot\hat{\mathbf{S}} $$
超精细相互作用能通常用总自旋 $\hat{\mathbf{F}} = \hat{\mathbf{I}} + \hat{\mathbf{J}}$ 来描述，对应的能量为：
$$ \Delta E_{\text{HFS}} = \frac{A}{2}[F(F+1) - I(I+1) - J(J+1)] $$
其中 $A$ 是超精细常数。
7.2 氢原子21cm线
氢原子超精细结构最著名的例子是基态（$n=1, l=0$）的超精细分裂。质子的核自旋 $I = 1/2$，电子自旋 $s = 1/2$，总自旋 $\mathbf{F} = \mathbf{I} + \mathbf{S}$可以取$F = 0$（单态，自旋反平行）和 $F = 1$（三重态，自旋平行）。
$F=1$和$F=0$之间的能量差对应于频率$\nu = 1420.4057517667\thinspace\text{MHz}$，即波长 $\lambda = 21.106\thinspace\text{cm}$。这就是著名的氢原子21厘米线。
天体物理意义
21cm线在天体物理学中具有极其重要的地位：
- 中性氢分布探测：银河系和河外星系中中性氢原子（HI 区）的分布可以通过21cm线来测绘，这是绘制银河系螺旋结构的主要手段。
- 宇宙学红移测量：21cm线可用于测量宇宙大尺度结构，探测宇宙再电离时期。
- SETI 搜索：21cm线被认为是一个"宇宙水洞"——地外文明可能使用这个频率进行星际通信，因为它是宇宙中最自然的频率标准之一。
- 暗物质探测：通过21cm线的宇宙学观测可以间接探测暗物质的性质。
八、塞曼效应
8.1 正常塞曼效应
当氢原子置于外磁场 $\mathbf{B} = B\hat{z}$ 中时，原子能级发生分裂，这就是塞曼效应。电子轨道磁矩与外磁场的相互作用为：
$$ \hat{H}_Z = \frac{e}{2m}(\hat{\mathbf{L}} + 2\hat{\mathbf{S}})\cdot\mathbf{B} $$
如果忽略自旋（即考虑正常塞曼效应），则：
$$ \hat{H}_Z = \frac{eB}{2m}\hat{L}_z $$
能级分裂为等间距的 $2l+1$ 个子能级：
$$ \Delta E = m_l\thinspace\mu_B B $$
其中 $\mu_B = e\hbar/(2m) \approx 5.79 \times 10^{-5}\thinspace\text{eV/T}$ 是玻尔磁子。
选择定则 $\Delta m_l = 0, \pm 1$导致光谱线分裂为三条：原频率$\omega_0$和$\omega_0 \pm \omega_L$（其中 $\omega_L = eB/(2m)$ 是拉莫尔频率）。
8.2 反常塞曼效应
当考虑自旋时，出现反常塞曼效应。在弱磁场中（$B$较小，精细结构未被破坏），使用$|n, l, j, m_j\rangle$ 基矢，塞曼能量为：
$$ \Delta E = g_J\thinspace{}m_j\thinspace\mu_B B $$
其中朗德 g 因子为：
$$ g_J = 1 + \frac{j(j+1) + s(s+1) - l(l+1)}{2j(j+1)} $$
$g_J$因子是$m_j$各子能级分裂间距的决定因素。对于$l=0$ 的 s 态，$g_J = 2$，分裂完全来自电子自旋。
8.3 与角动量与自旋的关系
塞曼效应是角动量理论的直接实验验证。$g_J$ 因子的公式来源于角动量的矢量耦合模型，是角动量与自旋中讨论的角动量代数（CG 系数、Wigner-Eckart 定理）的自然结果。朗德 g 因子的实验测量值与理论预言的高度一致，为角动量量子理论提供了强有力的支持。
九、斯塔克效应
9.1 外电场中的能级移动
当氢原子置于外电场 $\mathbf{E} = E\thinspace\hat{z}$ 中时，电子与电场的相互作用为：
$$ \hat{H}_S = eE\thinspace{}z = eE\thinspace{}r\cos\theta $$
由于氢原子能级（除 $n=1$ 外）具有简并性，外电场会在简并态之间产生非零的矩阵元，需要用简并微扰理论来处理。参见微扰理论中的相关讨论。
9.2 线性斯塔克效应
氢原子（及类氢离子）表现出线性斯塔克效应——能级移动与电场强度 $E$成正比。这是氢原子独有的特征，来源于$n$相同而$l$ 不同的态之间的简并。
对于 $n=2$ 能级，$2s$和$2p$ 态之间的微扰矩阵元为：
$$ \langle 200|\hat{H}_S|210\rangle = 3eE a_0 $$
对 $2 \times 2$ 矩阵对角化后，$n=2$ 能级分裂为三个子能级：
$$ E = E_2^{(0)},\thickspace E_2^{(0)} \pm 3eE a_0 $$
为什么其他原子没有线性斯塔克效应？
在非氢原子中，$n$相同但$l$不同的态通常不具有相同的能量（因为多电子屏蔽效应破坏了库仑势的简并性），因此$\hat{H}_S$ 在非简并态之间的对角元为零，线性斯塔克效应消失，仅留下二次斯塔克效应（$\propto E^2$）。
十、类氢离子
10.1 核电荷数 Z > 1
类氢离子（如 $\mathrm{He}^+$、$\mathrm{Li}^{2+}$、$\mathrm{Be}^{3+}$等）是仅有一个电子但核电荷数为$Z > 1$的离子。其薛定谔方程与氢原子相同，只需将库仑势中的$e^2$替换为$Ze^2$：
$$ V(r) = -\frac{Ze^2}{4\pi\varepsilon_0 r} $$
类氢离子的能级为：
$$ E_n = -\frac{Z^2 \cdot 13.6\thinspace\text{eV}}{n^2} $$
- $\mathrm{He}^+$（$Z=2$）：$E_1 = -54.4\thinspace\text{eV}$
- $\mathrm{Li}^{2+}$（$Z=3$）：$E_1 = -122.4\thinspace\text{eV}$
类氢离子的波函数也可以通过标度变换得到：$a_0 \to a_0/Z$，$R_{nl}(r) \to Z^{3/2}R_{nl}(Zr)$。
10.2 多电子原子的基础
氢原子和类氢离子的波函数为多电子原子的计算提供了基础。在微扰理论和变分法中，类氢轨道常被用作零级近似波函数。此外，氢原子求解过程中引入的量子数概念（$n, l, m$）和原子轨道（s, p, d, f...）构成了整个原子物理学和量子化学的语言体系。
从氢原子到氦原子，再到更复杂的多电子原子，量子力学基本假设提供了统一的理论框架，而波函数与薛定谔方程给出了求解实际问题的数学工具。氢原子的精确解作为理论的出发点，其重要性怎么强调都不为过。