"从一维到三维，不只是增加两个坐标轴——角动量出现了，自旋出现了，整个化学的基础出现了。"
氢原子的精确解是量子力学最伟大的胜利之一。

第4章三维量子力学：原子为什么不坍缩

4.1 故事：原子尺度的建筑设计师

22世纪，"纳构科技"公司专为客户设计原子尺度的电子器件。首席设计师苏然接到了一份看似不可能的任务：

"设计一个单原子晶体管。电子必须被精确约束在原子核附近，但又要有特定的能级跃迁来传输信号。"

苏然打开模拟软件，调出氢原子的三维波函数图像。屏幕上浮现出一个个概率云——有的像球（s轨道），有的像哑铃（p轨道），有的像更复杂的四叶草（d轨道）。

"经典物理告诉我们电子应该螺旋掉进原子核，"她对团队说，"但量子力学说电子会形成这些稳定的概率云。为什么呢？"

年轻的工程师小王指着屏幕："因为三维空间中的薛定谔方程有角动量这个守恒量。电子不只是有能量，还有角动量。角动量给电子一个’离心势垒’，阻止它无限接近原子核。"

"更关键的是，"苏然补充，"在三维中，能级不仅由主量子数 $n$ 决定，还有角量子数 $l$ 和磁量子数 $m$ 。这三个数一起，决定了电子云的形状、大小和空间取向。"

核心洞察：三维空间中的量子力学引入了角动量，它既是守恒量，也是保护原子免于坍缩的机制。

4.2 前置知识：球坐标系与球谐函数

4.2.1 球坐标系基础

在三维空间中描述一个点，除了直角坐标 $(x, y, z)$ ，还有另一种自然的选择——球坐标 $(r, \theta, \phi)$ ：

坐标	名称	定义	取值范围
$r$	径向距离	到原点的距离	$[0, \infty)$
$\theta$	极角（天顶角）	与 $z$ 轴的夹角	$[0, \pi]$
$\phi$	方位角	在 $xy$ 平面上的投影与 $x$ 轴的夹角	$[0, 2\pi)$

坐标变换关系：

$\begin{aligned} x &= r\sin\theta\cos\phi \\ y &= r\sin\theta\sin\phi \\ z &= r\cos\theta \end{aligned}$

体积元（由雅可比行列式导出）：

$dV = r^2 \sin\theta \, dr \, d\theta \, d\phi$

例题：计算半径为 $R$ 的球的体积，验证 $V = \frac{4}{3}\pi R^3$ 。

解答：

$V = \int_0^R r^2 dr \int_0^\pi \sin\theta \, d\theta \int_0^{2\pi} d\phi = \frac{R^3}{3} \cdot 2 \cdot 2\pi = \frac{4}{3}\pi R^3$

为什么用球坐标？ 因为原子核产生的库仑势是球对称的——从任何方向看，势能在相同距离上都一样。球坐标把"球对称"这个几何特征编码进了坐标系本身，让方程大大简化。这就像是面对一个圆形蛋糕，用极坐标描述比直角坐标自然得多。

4.2.2 球谐函数初步

球谐函数 $Y_l^m(\theta, \phi)$ 是角动量算符 $\hat{L}^2$ 和 $\hat{L}_z$ 的共同本征函数：

$\hat{L}^2 Y_l^m = \hbar^2 l(l+1) Y_l^m, \quad \hat{L}_z Y_l^m = \hbar m Y_l^m$

其中 $l = 0, 1, 2, ...$ 是角量子数， $m = -l, -l+1, ..., l$ 是磁量子数。

正交归一性：

$\int_0^\pi \sin\theta \, d\theta \int_0^{2\pi} d\phi \, Y_{l'}^{m'*} Y_l^m = \delta_{l'l} \delta_{m'm}$

宇称：在空间反演 $\vec{r} \to -\vec{r}$ （即 $r \to r$ , $\theta \to \pi - \theta$ , $\phi \to \phi + \pi$ ）下：

$Y_l^m(\pi - \theta, \phi + \pi) = (-1)^l Y_l^m(\theta, \phi)$

前几个球谐函数（显式表达式）：

$\begin{aligned} Y_0^0 &= \frac{1}{\sqrt{4\pi}} \\ Y_1^0 &= \sqrt{\frac{3}{4\pi}} \cos\theta \\ Y_1^{\pm 1} &= \mp \sqrt{\frac{3}{8\pi}} \sin\theta \, e^{\pm i\phi} \\ Y_2^0 &= \sqrt{\frac{5}{16\pi}} (3\cos^2\theta - 1) \end{aligned}$

物理直觉： $Y_0^0$ 是常数——s轨道是球对称的，电子在空间中任何方向的分布都一样。 $Y_1^0 \propto \cos\theta$ 在 $\theta = 0$ （北极）最大，在 $\theta = \pi/2$ （赤道）为零——这就是p轨道的哑铃形状。球谐函数就是"量子力学中的乐高积木"——它们构成了所有角向分布的完备基。

4.2.3 球谐函数的实数组合

化学中常用的轨道图像是实数组合：

$\begin{aligned} p_x &= \frac{1}{\sqrt{2}}(Y_1^{-1} - Y_1^{+1}) \propto \sin\theta\cos\phi \propto x/r \\ p_y &= \frac{i}{\sqrt{2}}(Y_1^{-1} + Y_1^{+1}) \propto \sin\theta\sin\phi \propto y/r \\ p_z &= Y_1^0 \propto \cos\theta \propto z/r \end{aligned}$

例题：验证 $p_x$ 轨道在 $x$ 轴方向（ $\theta = \pi/2$ , $\phi = 0$ ）取最大值，在 $y$ - $z$ 平面（ $x=0$ ）为零。

解答： $p_x \propto x/r$ ，在 $x$ 轴上 $x = r$ ，所以最大；在 $y$ - $z$ 平面上 $x = 0$ ，所以为零。这与化学中"p轨道沿x轴延伸"的图像完全一致。

4.3 三维薛定谔方程与球坐标

4.3.1 方程形式：从一维到三维的飞跃

三维时间依赖薛定谔方程：

$i\hbar \frac{\partial \Psi}{\partial t} = -\frac{\hbar^2}{2m}\nabla^2\Psi + V\Psi$

对于中心势场（ $V = V(r)$ ，只与到原点的距离有关），使用球坐标 $(r, \theta, \phi)$ 最方便。

拉普拉斯算符在球坐标中：

$\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}$

graph TD
    A[三维薛定谔方程] -->|"中心势场
V = V(r)"| B[球坐标分离变量]
    B --> C["r: 径向距离"]
    B --> D["θ: 极角"]
    B --> E["φ: 方位角"]
    C --> F["决定能级和
电子到核的距离"]
    D --> G["决定电子云的
南北极分布"]
    E --> H["决定电子云的
东西方向取向"]
    
    style A fill:#e3f2fd
    style B fill:#e8f5e9
    style F fill:#fff3e0
    style G fill:#fce4ec
    style H fill:#fce4ec

4.3.2 分离变量：把三维问题拆成三个一维问题

假设 $\psi(r, \theta, \phi) = R(r)Y(\theta, \phi)$ ，代入定态薛定谔方程后，方程分离为：

角向方程（球谐函数）：

$\frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\left(\sin\theta\frac{\partial Y}{\partial\theta}\right) + \frac{1}{\sin^2\theta}\frac{\partial^2 Y}{\partial\phi^2} = -l(l+1)Y$

解为球谐函数 $Y_l^m(\theta, \phi)$ ，其中 $l = 0, 1, 2, ...$ 是角量子数， $m = -l, -l+1, ..., l$ 是磁量子数。

径向方程：

$-\frac{\hbar^2}{2m}\frac{d^2 u}{dr^2} + \left[V(r) + \frac{\hbar^2}{2m}\frac{l(l+1)}{r^2}\right]u = Eu$

其中 $u(r) = rR(r)$ ，有效势多了一个离心势垒项 $\frac{\hbar^2 l(l+1)}{2mr^2}$ 。

物理直觉：离心势垒项是什么？想象你用绳子甩一个小球——球不会撞到你手上，因为旋转产生了"离心力"把它往外推。在量子力学中，角动量 $l$ 越大，离心势垒越强，电子越难接近原子核。这就是为什么s电子（ $l=0$ ）可以穿透到核附近，而d电子（ $l=2$ ）被排斥在外层。

graph TD
    A[三维薛定谔方程] -->|"球坐标 + 分离变量"| B[角向方程]
    A -->|"球坐标 + 分离变量"| C[径向方程]
    B --> D["球谐函数 Yₗᵐ"]
    C --> E["径向波函数 R(r)"]
    D --> F["决定形状
s, p, d, f..."]
    E --> G["决定大小
能级"]
    
    style D fill:#e3f2fd
    style E fill:#e8f5e9
    style F fill:#fff3e0
    style G fill:#fce4ec

4.4 角动量：旋转的量子化

4.4.1 经典角动量 vs 量子角动量

经典力学中 $\vec{L} = \vec{r} \times \vec{p}$ 。在量子力学中，角动量变为算符：

$\hat{L}_x = y\hat{p}_z - z\hat{p}_y = -i\hbar\left(y\frac{\partial}{\partial z} - z\frac{\partial}{\partial y}\right)$

类似定义 $\hat{L}_y$ 和 $\hat{L}_z$ 。

总角动量平方算符：

$\hat{L}^2 = \hat{L}_x^2 + \hat{L}_y^2 + \hat{L}_z^2$

与第3章的联系：角动量算符是厄米算符（可以验证），所以它们的本征值是实数，本征态可以构成完备基。这是第3章形式主义的直接应用。

graph LR
    A[经典角动量] -->|"L⃗ = r⃗ × p⃗"| B[量子角动量算符]
    B --> C["L̂ₓ, L̂ᵧ, L̂_z"]
    C --> D["L̂² = L̂ₓ² + L̂ᵧ² + L̂_z²"]
    D --> E["本征值: ℏ²l(l+1)"]
    C --> F["本征值: ℏm"]
    
    style B fill:#e3f2fd
    style E fill:#e8f5e9
    style F fill:#fff3e0

4.4.2 角动量的对易关系：量子力学最漂亮的代数结构

这是量子力学中最漂亮的代数结构之一：

$[\hat{L}_x, \hat{L}_y] = i\hbar\hat{L}_z, \quad [\hat{L}_y, \hat{L}_z] = i\hbar\hat{L}_x, \quad [\hat{L}_z, \hat{L}_x] = i\hbar\hat{L}_y$

三个分量互相不对易——这意味着你无法同时精确测量 $L_x$ 、 $L_y$ 和 $L_z$ 。最多只能精确测量一个分量（通常选 $L_z$ ）和总角动量平方 $L^2$ 。

物理直觉：为什么角动量的三个分量不能同时确定？想象一下：如果你精确知道 $L_z$ ，电子就在xy平面上做确定的圆周运动。但圆周运动意味着 $L_x$ 和 $L_y$ 在不断变化——就像旋转的陀螺，你知道它的转轴方向，但垂直于转轴的方向就在快速变化。这不是技术限制，是量子世界的内在属性。

完整推导：验证 $[\hat{L}_x, \hat{L}_y] = i\hbar\hat{L}_z$ 。

$\hat{L}_x = y\hat{p}_z - z\hat{p}_y, \quad \hat{L}_y = z\hat{p}_x - x\hat{p}_z$

$[\hat{L}_x, \hat{L}_y] = [y\hat{p}_z - z\hat{p}_y, z\hat{p}_x - x\hat{p}_z]$

展开四项，利用 $[z, \hat{p}_z] = i\hbar$ 和位置与位置、动量与动量都对易：

$= y\hat{p}_x[z, \hat{p}_z] - \hat{p}_y x[\hat{p}_z, z]$

$= y\hat{p}_x(i\hbar) - \hat{p}_y x(-i\hbar)$

$= i\hbar(y\hat{p}_x - x\hat{p}_y) = i\hbar\hat{L}_z$

升降算符：定义 $L_\pm = L_x \pm iL_y$ ，则有：

$[L_z, L_\pm] = \pm \hbar L_\pm, \quad [L^2, L_\pm] = 0$

$L_+$ 把 $m$ 增加1（"升"）， $L_-$ 把 $m$ 减少1（"降"）。利用这个代数结构，可以从一个本征态出发，生成全部 $2l+1$ 个态——这就是角动量本征值的完整推导方法。

4.4.3 本征值：角动量的量子化

总角动量平方： $\hat{L}^2 Y_l^m = \hbar^2 l(l+1) Y_l^m$ ， $l = 0, 1, 2, ...$
z分量： $\hat{L}_z Y_l^m = \hbar m Y_l^m$ ， $m = -l, -l+1, ..., l$

关键结论：

量子数	名称	允许值	物理意义
$l$	角量子数	$0, 1, 2, ...$	总角动量大小 $L = \sqrt{l(l+1)}\hbar$
$m$	磁量子数	$-l$ 到 $l$ 的整数	$L_z$ 分量，共 $2l+1$ 个取值

注意： $L_z$ 的最大值是 $l\hbar$ ，但总角动量 $L = \sqrt{l(l+1)}\hbar > l\hbar$ 。所以角动量矢量永远不能完全指向z轴方向——这是不确定性原理在角动量上的体现。

graph TD
    subgraph "l=2 的角动量取向（空间量子化）"
        A["L² = 6ℏ²"] --> B["m=2: Lz=2ℏ"]
        A --> C["m=1: Lz=1ℏ"]
        A --> D["m=0: Lz=0"]
        A --> E["m=-1: Lz=-1ℏ"]
        A --> F["m=-2: Lz=-2ℏ"]
    end
    
    style A fill:#e3f2fd
    style B fill:#fff3e0
    style C fill:#fff3e0
    style D fill:#fff3e0
    style E fill:#fff3e0
    style F fill:#fff3e0

历史背景：空间量子化的概念最早由索末菲（Sommerfeld）和玻尔提出。斯特恩-盖拉赫实验（1922年）首次直接观测到了空间量子化——银原子束在磁场中分裂成离散的两束，而不是连续分布。这个实验不仅证实了角动量的量子化，还意外发现了电子自旋。

4.5 氢原子：量子力学的皇冠明珠

4.5.1 库仑势：最简单的原子

库仑势：

$V(r) = -\frac{e^2}{4\pi\varepsilon_0}\frac{1}{r} = -\frac{k_e e^2}{r}$

这是一个吸引性的中心势，电子被原子核束缚。

物理直觉：库仑势的特点是"长程"和"奇点"。 $1/r$ 在 $r \to 0$ 时发散，在 $r \to \infty$ 时缓慢衰减。这种"既深又广"的势阱，恰好能产生离散的束缚态和连续的散射态。如果势阱太浅（比如屏蔽库仑势），可能就没有束缚态；如果势阱太深，能级结构又会完全不同。自然选择的 $1/r$ 形式，恰好给了我们丰富的原子光谱。

4.5.2 双粒子系统：约化质量

在之前的讨论中，我们假设原子核是静止的。实际上，电子和质子都绕着它们的质心运动。在质心坐标系中，问题约化为一个约化质量为 $\mu$ 的单粒子在库仑势中运动：

$\mu = \frac{m_e m_p}{m_e + m_p}$

对于氢原子：

$\mu = \frac{m_e m_p}{m_e + m_p} = \frac{9.109 \times 10^{-31} \times 1.673 \times 10^{-27}}{9.109 \times 10^{-31} + 1.673 \times 10^{-27}} \approx 9.104 \times 10^{-31} \text{ kg}$

约化质量比电子质量小约 $0.05\%$ 。这意味着氢原子的能级比用 $m_e$ 计算的结果稍微低一点。

例题：计算氢原子基态能量的约化质量修正。

解答：

$E_1 = -\frac{\mu}{2\hbar^2}\left(\frac{e^2}{4\pi\varepsilon_0}\right)^2 = -\frac{\mu}{m_e} \times 13.6 \text{ eV}$

$\frac{\mu}{m_e} = \frac{m_p}{m_e + m_p} = \frac{1.673 \times 10^{-27}}{1.674 \times 10^{-27}} \approx 0.999455$

$E_1 = -0.999455 \times 13.6 = -13.598 \text{ eV}$

即基态能量比无穷核质量近似低约 $0.007$ eV。这看起来很小，但在精密光谱学中非常重要——氢原子光谱的测量精度达到 $10^{-12}$ 量级，约化质量修正是必须考虑的。

4.5.3 径向方程的求解：级数法

氢原子的径向方程（令 $u(r) = rR(r)$ ）：

$-\frac{\hbar^2}{2\mu}\frac{d^2 u}{dr^2} + \left[-\frac{e^2}{4\pi\varepsilon_0 r} + \frac{\hbar^2 l(l+1)}{2\mu r^2}\right]u = Eu$

Griffiths在书中用级数法详细推导，步骤如下：

无量纲化：令 $\rho = \kappa r$ ， $\kappa = \sqrt{-2\mu E}/\hbar$ （对于束缚态 $E < 0$ ），方程变为标准形式。
渐进分析：
- $r \to \infty$ 时， $u \sim e^{-\rho}$ （指数衰减，束缚态要求）
- $r \to 0$ 时， $u \sim \rho^{l+1}$ （离心势垒要求波函数在原点为零）
设解的形式： $u(\rho) = \rho^{l+1} e^{-\rho} v(\rho)$ ，代入后得到 $v(\rho)$ 满足的微分方程。
幂级数展开： $v(\rho) = \sum_{j=0}^{\infty} c_j \rho^j$ 。递推关系要求级数必须在某处截断（否则在 $r \to \infty$ 时发散）。
截断条件：定义主量子数 $n = j_{max} + l + 1$ ，得到能量量子化条件。

4.5.4 求解结果：量子力学最优美的答案

通过求解径向方程，得到：

能量本征值（惊人的简单！）：

$E_n = -\frac{\mu}{2\hbar^2}\left(\frac{e^2}{4\pi\varepsilon_0}\right)^2\frac{1}{n^2} = -\frac{13.6 \text{ eV}}{n^2}, \quad n = 1, 2, 3, ...$

物理意义：

$n=1$ （基态）： $E_1 = -13.6$ eV，电子被紧紧束缚
$n=2$ ： $E_2 = -3.4$ eV
$n=3$ ： $E_3 = -1.51$ eV
$n \to \infty$ ： $E_\infty = 0$ ，电离阈值

主量子数 $n$ 决定了能量。但波函数还需要 $l$ 和 $m$ ：

允许的组合：对于给定的 $n$ ， $l$ 可以取 $0, 1, ..., n-1$ ；对于给定的 $l$ ， $m$ 可以取 $-l, ..., l$ 。

径向波函数 $R_{nl}(r)$ 与连带拉盖尔多项式相关：

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

其中 $\rho = 2r/(na_0)$ ，a_0 = 4\pi\varepsilon_0\hbar^2/(\mu e^2) \approx 0.529 \text{ Å} 是玻尔半径。

graph TD
    subgraph 氢原子能级
        A["n=1, -13.6eV
基态"] -->|"Lyman系
紫外"| B["n=2, -3.4eV"]
        B -->|"Balmer系
可见光"| C["n=3, -1.51eV"]
        C -->|"Paschen系
红外"| D["n=4, -0.85eV"]
        D --> E["n=∞, 0eV
电离"]
    end
    
    style A fill:#e3f2fd
    style B fill:#e8f5e9
    style C fill:#fff3e0
    style D fill:#fce4ec
    style E fill:#ffebee

简并度：对于给定的 $n$ ，总简并度为：

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

考虑自旋后，简并度为 $2n^2$ 。

物理直觉：为什么能量只依赖 $n$ 而不依赖 $l$ ？这是库仑势的特殊性质！对于一般的中心势，能量会同时依赖 $n$ 和 $l$ 。但 $1/r$ 势恰好有额外的对称性（SO(4)动力学对称性），导致 $l$ 的"偶然简并"。这就像一个巧合：自然选择的势能形式，让能级结构格外简单。

4.5.5 径向分布函数与最概然半径

电子在半径 $r$ 到 $r+dr$ 的球壳中出现的概率为：

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

例题：计算氢原子基态（ $1s$ ）和第一激发态（ $2s$ ）的最概然半径。

解答：

基态波函数： $R_{10}(r) = 2a_0^{-3/2} e^{-r/a_0}$

$P(r) = |R_{10}|^2 r^2 = \frac{4}{a_0^3} r^2 e^{-2r/a_0}$

令 $dP/dr = 0$ ：

$\frac{d}{dr}(r^2 e^{-2r/a_0}) = 2r e^{-2r/a_0} - \frac{2r^2}{a_0}e^{-2r/a_0} = 0$

r = a_0 \approx 0.529 \text{ Å}

基态的最概然半径恰好是玻尔半径！

对于 $2s$ 态： $R_{20}(r) = \frac{1}{\sqrt{2}a_0^{3/2}}(1 - \frac{r}{2a_0})e^{-r/(2a_0)}$

$P(r) = r^2 |R_{20}|^2$ ，求导后有两个极值点： $r = (3 \pm \sqrt{5})a_0$ ，即 $r \approx 0.764a_0$ 和 $r \approx 5.236a_0$ 。后者是全局最大值，所以 $2s$ 态的最概然半径约为 $5.24 a_0 \approx 2.77$ Å。

4.5.6 光谱学与量子跃迁：从理论到实验

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

E_{\text{光子}} = E_{n_i} - E_{n_f} = 13.6\left(\frac{1}{n_f^2} - \frac{1}{n_i^2}\right) \text{ eV}

这就是氢原子光谱，解释了19世纪末发现的巴尔末系、莱曼系等谱线。

例题：计算从 $n=3$ 到 $n=2$ 跃迁发射的光子波长（巴尔末系 $H_\alpha$ 线）。

解答：

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

\lambda = \frac{hc}{E} = \frac{4.14 \times 10^{-15} \text{ eV·s} \times 3 \times 10^8 \text{ m/s}}{1.89 \text{ eV}} \approx 656 \text{ nm}

这正是红色的 $H_\alpha$ 线！量子力学不仅"解释"了光谱，还精确预测了每一条谱线的位置。

4.5.7 电子云的形状：原子的"指纹"

记号	$n$	$l$	$m$	形状	名称
$1s$	1	0	0	球对称	基态
$2s$	2	0	0	球对称 + 径向节点
$2p$	2	1	-1,0,1	哑铃形	p轨道
$3d$	3	2	-2,…,2	四叶草/环形	d轨道

物理直觉：s轨道是球对称的，因为 $l=0$ 意味着没有角动量，电子可以从任何角度接近原子核。p轨道像哑铃，因为 $l=1$ 的电子有一个单位角动量，在某个平面内运动，但在垂直方向上被"推开"。d轨道更复杂，因为 $l=2$ 的电子在两个方向上都有角动量分量。

graph LR
    subgraph 轨道形状与化学
        A["l=0
s轨道"] --> B[球对称]
        C["l=1
p轨道"] --> D["哑铃形
三种取向"]
        E["l=2
d轨道"] --> F["四叶草
五种取向"]
        B --> G["σ键"]
        D --> H["π键"]
        F --> I[配位化学]
    end
    
    style B fill:#e3f2fd
    style D fill:#e8f5e9
    style F fill:#fff3e0
    style G fill:#fce4ec
    style H fill:#fce4ec
    style I fill:#fce4ec

4.6 自旋：内禀角动量

4.6.1 实验证据：斯特恩-盖拉赫实验

Stern-Gerlach 实验（1922年）发现：银原子束在非均匀磁场中分裂成两束，而不是连续分布。

这意味着电子有一种内禀的角动量，与空间运动无关。它只能取两个离散值："自旋向上"和"自旋向下"。

历史注记：斯特恩和盖拉赫做这个实验时，本想验证空间量子化理论。他们预期看到三条线（ $m = -1, 0, +1$ ），结果只看到两条。乌伦贝克（Uhlenbeck）和古德斯密特（Goudsmit）在1925年提出"电子自旋"假说，解释说额外的一条线来自电子的内禀角动量。泡利最初反对这个"经典旋转"的图像，但后来不得不接受数学结果。

graph TD
    A["Stern-Gerlach实验"] -->|银原子束| B[非均匀磁场]
    B --> C[分裂成两束]
    C --> D["电子自旋
s = ½"]
    D --> E["m_s = +½"]
    D --> F["m_s = -½"]
    
    style A fill:#e3f2fd
    style C fill:#e8f5e9
    style D fill:#fff3e0

4.6.2 自旋1/2的数学：泡利矩阵

自旋算符 $\vec{S}$ 满足与轨道角动量相同的对易关系，但本征值不同：

$S^2 |s, m_s\rangle = \hbar^2 s(s+1)|s, m_s\rangle$

对于电子， $s = 1/2$ （因此叫"自旋1/2"）， $m_s = \pm 1/2$ 。

泡利矩阵（自旋算符的矩阵表示）：

$S_x = \frac{\hbar}{2}\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}, \quad S_y = \frac{\hbar}{2}\begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}, \quad S_z = \frac{\hbar}{2}\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}$

自旋态用旋量（二分量列向量）表示：

$|+\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \quad |-\rangle = \begin{pmatrix} 0 \\ 1 \end{pmatrix}$

例题：验证 $|+\rangle$ 是 $S_z$ 的本征态，本征值为 $+\hbar/2$ 。

解答：

$S_z|+\rangle = \frac{\hbar}{2}\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}\begin{pmatrix} 1 \\ 0 \end{pmatrix} = \frac{\hbar}{2}\begin{pmatrix} 1 \\ 0 \end{pmatrix} = +\frac{\hbar}{2}|+\rangle$

4.6.3 自旋的物理意义：不是"自转"

自旋是量子力学的内禀属性，不是电子在"自转"。如果把电子想象成一个经典小球，要产生观测到的磁矩，它的表面转速需要超过光速——这是不可能的。

物理直觉：自旋更像是"标签"或"内部自由度"。每个电子自带这个标签，取值为"上"或"下"。它不是空间运动产生的，而是时空对称性（洛伦兹群）在量子力学中的体现。就像质量和电荷是粒子的内禀属性一样，自旋也是。

与第5章的联系：自旋的取值（整数或半整数）直接决定了粒子是玻色子还是费米子，进而决定了泡利不相容原理是否适用。自旋是"统计性质"的根源。

graph TD
    A[角动量] --> B["轨道角动量 L"]
    A --> C["自旋角动量 S"]
    B -->|空间运动| D["l = 0,1,2,..."]
    C -->|内禀属性| E["s = 0, ½, 1,..."]
    D --> F["电子: l任意"]
    E --> G["电子: s = ½"]
    G --> H[费米子]
    G --> I[泡利不相容原理]
    
    style B fill:#e3f2fd
    style C fill:#fff3e0
    style F fill:#e8f5e9
    style G fill:#e8f5e9
    style H fill:#ffebee
    style I fill:#ffebee

4.7 多电子原子：从氢到真实世界

4.7.1 屏蔽效应与有效核电荷

氢原子只有一个电子，问题简单。但多电子原子中，电子之间还有相互作用：

$V = -\frac{Ze^2}{4\pi\varepsilon_0 r} + \sum_{i<j} \frac{e^2}{4\pi\varepsilon_0 |\vec{r}_i - \vec{r}_j|}$

第二项是电子-电子排斥，让问题无法精确求解。

近似方法：把内层电子对核电荷的屏蔽效应用有效核电荷 $Z_{\text{eff}}$ 描述。外层电子感受到的势近似为：

$V_{\text{eff}}(r) = -\frac{Z_{\text{eff}}e^2}{4\pi\varepsilon_0 r}$

物理直觉：内层电子像一层"防护服"，把原子核的正电荷屏蔽了一部分。外层电子感受不到完整的 $Z$ 个正电荷，而是更少的 $Z_{\text{eff}}$ 。这就是为什么碱金属（如钠）的外层电子很容易失去——它感受到的有效核电荷很小，束缚很弱。

例题：用简单模型估算锂原子（ $Z=3$ ） $2s$ 电子的有效核电荷。

解答：

锂的基态电子排布是 $1s^2 2s^1$ 。 $1s$ 电子靠近原子核，几乎可以完全屏蔽两个正电荷。外层 $2s$ 电子感受到的有效核电荷约为：

$Z_{\text{eff}} = Z - \sigma = 3 - 2 \times 0.85 = 1.3$

其中 $\sigma = 0.85$ 是每个 $1s$ 电子对 $2s$ 电子的屏蔽常数（Slater规则）。由此估算 $2s$ 电子的能量：

$E_{2s} \approx -\frac{Z_{\text{eff}}^2}{n^2} \times 13.6 \text{ eV} = -\frac{1.3^2}{4} \times 13.6 \approx -5.75 \text{ eV}$

实验测得锂的电离能为 $5.39$ eV，与这个粗略估算同数量级。

4.7.2 洪德规则与电子排布

除了泡利不相容原理，电子排布还遵循洪德规则（Hund’s Rules）：

最大自旋多重度：电子优先以平行自旋占据不同的轨道
最小轨道角动量：对于给定自旋多重度，取最小 $L$
自旋-轨道耦合：对于未满半满壳层， $J = |L-S|$ ；对于超过半满的壳层， $J = L+S$

物理直觉：洪德第一条规则意味着电子"喜欢"自旋平行。为什么？因为平行自旋的波函数在空间上是反对称的（交换两个自旋平行的电子，自旋部分不变，空间部分必须变号），这降低了电子靠近的概率，从而减小了库仑排斥能。这是一个"量子力学为了减少能量"而自发产生的关联效应。

graph TD
    A[多电子原子] --> B[屏蔽效应]
    A --> C[泡利不相容原理]
    A --> D[洪德规则]
    B --> E["有效核电荷 Z_eff"]
    C --> F[电子分层排布]
    D --> G[优先平行自旋]
    E --> H[外层电子易失去]
    F --> I[元素周期表结构]
    G --> J[降低库仑排斥]
    
    style A fill:#e3f2fd
    style E fill:#e8f5e9
    style F fill:#fff3e0
    style G fill:#fce4ec

4.8 本章总结

graph TD
    A[第4章核心] --> B[三维薛定谔方程]
    A --> C[角动量量子化]
    A --> D[氢原子精确解]
    A --> E[自旋]
    A --> F[多电子原子]
    
    B -->|球坐标分离变量| G["径向 × 角向"]
    C -->|"L², Lz本征值"| H["l(l+1)ℏ², mℏ"]
    D -->|库仑势| I["Eₙ = -13.6/n² eV"]
    E -->|"Stern-Gerlach"| J["s = ½, 两个态"]
    F -->|屏蔽效应| K[有效核电荷]
    
    G --> L["离心势垒
阻止坍缩"]
    H --> M["空间量子化
2l+1个态"]
    I --> N["解释了
氢原子光谱"]
    J --> O["泡利矩阵
二分量旋量"]
    K --> P["元素周期表
化学性质"]
    
    L --> Q["第5章泡利原理
电子分层基础"]
    M --> R["原子光谱
塞曼效应"]
    N --> S["量子力学
最伟大的胜利"]
    O --> T["核磁共振
量子计算"]
    P --> U["化学键
材料科学"]
    
    style L fill:#e3f2fd
    style M fill:#e8f5e9
    style N fill:#fff3e0
    style O fill:#fce4ec
    style P fill:#fce4ec

带走的三句话：

角动量是三维量子力学的核心——它解释了为什么原子不坍缩，为什么能级有简并。 离心势垒是电子的"防护罩"。
氢原子的精确解是量子力学的胜利： $E_n \propto -1/n^2$ ，解释了全部原子光谱。 从紫外到可见光到红外，一切谱线尽在其中。
自旋是内禀属性，不是空间运动——它是量子世界的"新维度"。 自旋的取值决定了粒子统计，进而决定了整个物质世界。

"氢原子是量子力学的实验室。如果你理解了氢原子，你就理解了量子力学的一半。" —— Griffiths

4.9 练习与思考

球坐标系积分：用球坐标计算半径为 $R$ 的球内，满足 $r < R/2$ 且 $\theta < \pi/2$ （上半球内层）的体积。验证你的结果。
球谐函数计算：利用显式表达式，验证 $Y_1^1$ 和 $Y_1^{-1}$ 的正交归一性。计算 $|Y_1^1|^2$ 在 $\theta = \pi/2$ 的值——这对应p轨道在赤道平面的概率密度。
角动量对易子：验证 $[\hat{L}_x, \hat{L}_y] = i\hbar\hat{L}_z$ 。提示：用位置-动量的基本对易关系 $[x, p_x] = i\hbar$ 。然后思考：为什么这个对易关系意味着 $L_x$ 和 $L_y$ 不能同时精确测量？
约化质量修正：氘原子（一个电子绕氘核运动，氘核质量约为 $3.34 \times 10^{-27}$ kg）的约化质量是多少？与氢原子相比，氘原子的基态能量偏移了多少（用eV和百分比表示）？这在实验上可观测吗？
氢原子能级简并：对于 $n=3$ ，有多少个不同的量子态 $(n, l, m)$ ？考虑自旋后总数是多少？为什么化学中的"第三电子层"最多容纳18个电子？画出 $n=3$ 的所有可能态的能量简并图。
径向分布函数：画出 $P(r) = r^2 |R_{10}|^2$ 的草图（不需要精确绘制，标出峰值位置和衰减趋势）。为什么尽管 $|\psi_{1s}(0)|^2$ 在原点不为零，但最概然半径却在 $r = a_0$ 处？
自旋测量：一个电子处于自旋态 $|\chi\rangle = \frac{1}{\sqrt{2}}\begin{pmatrix}1 \\ 1\end{pmatrix}$ （这是 $S_x$ 的本征态）。测量 $S_z$ 可能得到什么结果？各自的概率是多少？测量 $S_x$ 呢？如果测量 $S_z$ 得到 $+\hbar/2$ ，再测量 $S_x$ ，结果会如何？这体现了什么量子力学原理？
有效核电荷估算：用Slater规则估算碳原子（ $Z=6$ ，电子排布 $1s^2 2s^2 2p^2$ ） $2p$ 电子的有效核电荷 $Z_{\text{eff}}$ ，并解释为什么碳的电离能（11.26 eV）远高于锂的电离能（5.39 eV）。

第4章 三维量子力学：原子为什么不坍缩