1. 狄拉克方程的原始形式
\[ i\hbar\frac{\partial\psi}{\partial t} = H_D\psi \]
其中
\[ H_D = c,\boldsymbol{\alpha}\cdot\mathbf{p} + \beta mc^2 + V \]
这里我加上了外势 \(V\)（比如氢原子中的库仑势），因为自旋-轨道耦合项正是从 \(V\) 来的。
在狄拉克表示中：
\[ \boldsymbol{\alpha} = \begin{pmatrix} 0 & \boldsymbol{\sigma}\ \boldsymbol{\sigma} & 0 \end{pmatrix}, \quad \beta = \begin{pmatrix} I & 0\ 0 & -I \end{pmatrix} \]
把 4 分量旋量写成两个 2 分量旋量：
\[ \psi = \begin{pmatrix} \phi\ \chi \end{pmatrix} \]
其中 \(\phi\) 是大分量（非相对论极限下主导），\(\chi\) 是小分量。
2. 分块写出方程
代入 \(H_D\)：
\[ i\hbar\frac{\partial}{\partial t} \begin{pmatrix}\phi\\chi\end{pmatrix}
\begin{pmatrix} mc^2+V & c,\boldsymbol{\sigma}\cdot\mathbf{p}\ c,\boldsymbol{\sigma}\cdot\mathbf{p} & -mc^2+V \end{pmatrix} \begin{pmatrix}\phi\\chi\end{pmatrix} \]
得到两个耦合方程：
\[ i\hbar\dot{\phi} = (mc^2+V)\phi + c,\boldsymbol{\sigma}\cdot\mathbf{p},\chi \tag{1} \]
\[ i\hbar\dot{\chi} = c,\boldsymbol{\sigma}\cdot\mathbf{p},\phi + (-mc^2+V)\chi \tag{2} \]
3. 非相对论极限：提取静止质量能量
在非相对论极限下，能量主要来自静止质量：
\[ E \approx mc^2 + E_{\text{NR}} \]
其中 \(E_{\text{NR}}\) 远小于 \(mc^2\)。为此做变换：
\[ \begin{pmatrix}\phi\\chi\end{pmatrix}
e^{-imc^2t/\hbar} \begin{pmatrix}\tilde{\phi}\\tilde{\chi}\end{pmatrix} \]
这样 \(\tilde{\phi},\tilde{\chi}\) 随时间变化缓慢。代入后：
\[ i\hbar\dot{\tilde{\phi}} = V\tilde{\phi} + c,\boldsymbol{\sigma}\cdot\mathbf{p},\tilde{\chi} \tag{3} \]
\[ i\hbar\dot{\tilde{\chi}} = c,\boldsymbol{\sigma}\cdot\mathbf{p},\tilde{\phi} + (-2mc^2+V)\tilde{\chi} \tag{4} \]
4. 消去小分量
在非相对论极限下，\(\tilde{\chi}\) 随时间变化很慢，且 \(|V|\ll mc^2\)，所以 (4) 式中：
\[ i\hbar\dot{\tilde{\chi}} \ll 2mc^2\tilde{\chi} \]
可以近似忽略左边，得到：
\[ 0 \approx c,\boldsymbol{\sigma}\cdot\mathbf{p},\tilde{\phi} - 2mc^2\tilde{\chi} \]
即：
\[ \tilde{\chi} \approx \frac{\boldsymbol{\sigma}\cdot\mathbf{p}}{2mc},\tilde{\phi} \]
这是小分量用大分量表示的近似关系。
5. 得到大分量的有效方程
把 \(\tilde{\chi}\) 代回 (3)：
\[ i\hbar\dot{\tilde{\phi}} = V\tilde{\phi} + c,\boldsymbol{\sigma}\cdot\mathbf{p} \cdot \frac{\boldsymbol{\sigma}\cdot\mathbf{p}}{2mc},\tilde{\phi} \]
即：
\[ i\hbar\dot{\tilde{\phi}} = \left[ V + \frac{(\boldsymbol{\sigma}\cdot\mathbf{p})^2}{2m} \right]\tilde{\phi} \]
利用泡利矩阵恒等式：
\[ (\boldsymbol{\sigma}\cdot\mathbf{A})(\boldsymbol{\sigma}\cdot\mathbf{B}) = \mathbf{A}\cdot\mathbf{B} + i\boldsymbol{\sigma}\cdot(\mathbf{A}\times\mathbf{B}) \]
取 \(\mathbf{A}=\mathbf{B}=\mathbf{p}\)，则 \(\mathbf{p}\times\mathbf{p}=0\)，所以：
\[ (\boldsymbol{\sigma}\cdot\mathbf{p})^2 = \mathbf{p}^2 \]
于是：
\[ i\hbar\dot{\tilde{\phi}} = \left[ \frac{\mathbf{p}^2}{2m} + V \right]\tilde{\phi} \]
这就是薛定谔方程。到此为止，自旋还没有显现，因为 \(\boldsymbol{\sigma}\) 被平方消掉了。
6. 如何得到精细结构项？
上面的推导太粗糙了：忽略了 \(\dot{\tilde{\chi}}\)、忽略了 \(V\) 对小分量的影响、没有考虑高阶修正。要得到 \(H_{\text{NR}}\) 中那些修正项，需要更系统地做Foldy-Wouthuysen 变换。
6.1 Foldy-Wouthuysen 变换的思路
FW 变换的目标是：通过一系列幺正变换，把哈密顿量中奇算子（把大分量和小分量耦合的项，如 \(\boldsymbol{\alpha}\cdot\mathbf{p}\)）消去，使哈密顿量变成块对角形式：
\[ H_{\text{FW}} = \begin{pmatrix} H_+ & 0\ 0 & H_- \end{pmatrix} \]
其中 \(H_+\) 描述正能态（电子），\(H_-\) 描述负能态（正电子）。然后对 \(H_+\) 做非相对论展开，就得到 \(H_{\text{NR}}\)。
6.2 展开到 \(1/c^2\) 阶的结果
经过 FW 变换并展开到 \(1/c^2\) 阶，正能态的有效哈密顿量为：
\[ H_{\text{NR}} = \frac{\mathbf{p}^2}{2m} + V
- \frac{\mathbf{p}^4}{8m^3c^2}
- \frac{\hbar^2}{8m^2c^2}\nabla^2V
- \frac{\hbar}{4m^2c^2}\boldsymbol{\sigma}\cdot(\nabla V\times\mathbf{p})
- \cdots \]
下面逐项解释来源。
7. 各项的物理来源
7.1 \(\frac{\mathbf{p}^2}{2m}+V\)：零阶项
就是普通的非相对论哈密顿量。
7.2 \(-\frac{\mathbf{p}^4}{8m^3c^2}\)：相对论动能修正
来自相对论能量展开：
\[ E = \sqrt{p^2c^2+m^2c^4} = mc^2 + \frac{p^2}{2m} - \frac{p^4}{8m^3c^2} + \cdots \]
减去静止质量 \(mc^2\) 后，第二项就是它。它导致原子能级的质量速度修正。
7.3 \(\frac{\hbar^2}{8m^2c^2}\nabla^2V\)：达尔文项
这一项来自狄拉克方程中小分量在势场中的“抖动”（Zitterbewegung）。它的物理意义是：电子并不是一个点，而是在康普顿波长尺度内快速抖动，感受到的势是空间平均后的势。对库仑势：
\[ V = -\frac{Ze^2}{4\pi\epsilon_0 r} \]
有：
\[ \nabla^2 V = \frac{Ze^2}{\epsilon_0}\delta^3(\mathbf{r}) \]
所以达尔文项只对 \(s\) 态（\(|\psi(0)|^2\neq 0\)）有贡献，导致 \(s\) 态能级上移。
7.4 \(\frac{\hbar}{4m^2c^2}\boldsymbol{\sigma}\cdot(\nabla V\times\mathbf{p})\)：自旋-轨道耦合
这一项是自旋-轨道耦合。对中心势：
\[ \nabla V = \frac{1}{r}\frac{dV}{dr}\mathbf{r} \]
所以：
\[ \nabla V\times\mathbf{p} = \frac{1}{r}\frac{dV}{dr}\mathbf{r}\times\mathbf{p} = \frac{1}{r}\frac{dV}{dr}\mathbf{L} \]
于是：
\[ H_{\text{SO}} = \frac{\hbar}{4m^2c^2}\frac{1}{r}\frac{dV}{dr}\boldsymbol{\sigma}\cdot\mathbf{L} \]
对库仑势 \(V=-Ze^2/(4\pi\epsilon_0 r)\)：
\[ \frac{1}{r}\frac{dV}{dr} = \frac{Ze^2}{4\pi\epsilon_0 r^3} \]
所以：
\[ H_{\text{SO}} = \frac{Ze^2\hbar}{4m^2c^2(4\pi\epsilon_0)} \frac{\boldsymbol{\sigma}\cdot\mathbf{L}}{r^3} \]
注意：这里有一个因子 2 的“托马斯进动”问题。严格的相对论推导给出的是 \(\frac{1}{2m^2c^2}\) 而不是 \(\frac{1}{4m^2c^2}\)，差一个 1/2 来自电子自旋的托马斯进动。上面 FW 变换直接得到的是 \(\frac{\hbar}{4m^2c^2}\)，但通常教材中会说明需要额外考虑托马斯进动，最终系数为 \(\frac{1}{2m^2c^2}\)。
7.5 更高阶项
继续展开会得到更多项，如：
- 外加磁场下的塞曼项
- 电场梯度导致的高阶修正
- 等等
8. 总结：\(H_{\text{NR}}\) 的完整来源
\[ \boxed{ H_{\text{NR}} = \frac{\mathbf{p}^2}{2m} + V
- \frac{\mathbf{p}^4}{8m^3c^2}
- \frac{\hbar^2}{8m^2c^2}\nabla^2V
- \frac{\hbar}{4m^2c^2}\boldsymbol{\sigma}\cdot(\nabla V\times\mathbf{p})
- \cdots } \]
| 项 | 名称 | 来源 |
|---|---|---|
| \(\frac{p^2}{2m}+V\) | 非相对论哈密顿量 | 零阶近似 |
| \(-\frac{p^4}{8m^3c^2}\) | 相对论动能修正 | 相对论能量展开 |
| \(\frac{\hbar^2}{8m^2c^2}\nabla^2V\) | 达尔文项 | 小分量抖动 / Zitterbewegung |
| \(\frac{\hbar}{4m^2c^2}\boldsymbol{\sigma}\cdot(\nabla V\times\mathbf{p})\) | 自旋-轨道耦合 | 大分量与小分量耦合 + 势场 |
推导方法：从狄拉克方程出发，分离大分量 \(\phi\) 和小分量 \(\chi\)，消去小分量得到大分量的有效方程，再做系统展开（或 FW 变换）到 \(1/c^2\) 阶，就得到 \(H_{\text{NR}}\)。
关键点：自旋-轨道耦合项中的 \(\boldsymbol{\sigma}\) 来自狄拉克矩阵结构，\(\mathbf{L}\) 来自 \(\nabla V\times\mathbf{p}\)，两者结合正是自旋与轨道角动量的耦合。