贡献者： _Eden_; addis
下文中所采用的电磁单位制为高斯单位制或洛伦兹-亥维赛单位制（Lorentz–Heaviside Units）1。
泡利方程可以写成以下的形式
\begin{equation}
\boxed{
i\hbar \frac{\partial}{\partial{t}} \varphi = \left[\frac{(-i\hbar\nabla - (e/c) \boldsymbol{\mathbf{A}} )^2}{2m}+e\phi-\frac{e}{mc} \boldsymbol{\mathbf{S}} \cdot \boldsymbol{\mathbf{B}} \right] \varphi}~.
\end{equation}
其中 $\varphi$ 是一个二分量的波函数，$e=-|q_e|$ 代表电子的电荷量。它可以有效地描述的非相对论情形下电磁场中自旋 $1/2$ 粒子的行为。
1. 从狄拉克方程到泡利方程
预备知识 2 狄拉克方程的非相对论近似
，
电磁场中的狄拉克方程
从式 10 出发：
\begin{equation}
i\hbar \frac{\partial}{\partial{t}} \varphi = \frac{1}{2m} \left[ \boldsymbol{\mathbf{\sigma}} \cdot \left( \boldsymbol{\mathbf{P}} - \frac{q}{c} \boldsymbol{\mathbf{A}} \right) \right] ^2 \varphi + e\phi \varphi~.
\end{equation}
并利用
式 5 开始的一系列推导对它进行化简，可以直接得到泡利方程。这种推导方式简单粗暴，直接取相对论性狄拉克方程的非相对论极限得到最终的解。
其中 $\frac{e}{mc} \boldsymbol{\mathbf{S}} \cdot \boldsymbol{\mathbf{B}} $ 是自旋与磁场耦合项。代表电子磁矩与外磁场之间的相互作用。如果将这一项舍去将可以得到非相对论情形下电磁场中无自旋带点粒子的薛定谔方程 式 1 。
2. 从自旋 $1/2$ 粒子的非相对论方程到泡利方程
或者我们也可以换一个角度，这种思路直接从非相对论的 Galilean 变换出发 [1]，依据非相对论情形下的时空对称性得到量子力学的算符以及电子波函数所满足的非相对论方程。虽然这种推导看起来拐弯抹角，但其实和狄拉克方程的推导是完全一致的。大部分的推导已经在文章中展现过，我们这里直接采用 中的结果。对于自旋 $1/2$ 粒子，描述它的非相对论性波函数具有 $4$ 个分量，可以表示为 $\psi= \begin{pmatrix}\varphi\\\chi\end{pmatrix} $，其中 $\varphi,\chi$ 为双分量的波函数，且满足方程式 11 ：
\begin{equation}
\begin{aligned}
& \left[\hat H-\frac{( \boldsymbol{\mathbf{\sigma}} \cdot \hat{ \boldsymbol{\mathbf{P}} })^2}{2m} \right] \varphi=0~,\\
&
\left[\hat H-\frac{( \boldsymbol{\mathbf{\sigma}} \cdot \hat{ \boldsymbol{\mathbf{P}} })^2}{2m} \right] \chi=0~.
\end{aligned}
\end{equation}
在电磁场中，需要对 $\hat H,\hat{ \boldsymbol{\mathbf{P}} }$ 作如下的替换：
\begin{equation}
\hat H\rightarrow \hat H-e\phi, \hat{ \boldsymbol{\mathbf{P}} }\rightarrow \hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} ~.
\end{equation}
其中 $\phi$ 为电势，$ \boldsymbol{\mathbf{A}} $ 为磁矢势。那么
\begin{equation}
\left[ \boldsymbol{\mathbf{\sigma}} \cdot(\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} ) \right] ^2=(\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} )^2+i \boldsymbol{\mathbf{\sigma}} \cdot[(\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} )\times (\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} )]~.
\end{equation}
注意这里的 $\hat \sigma,\hat{ \boldsymbol{\mathbf{P}} }, \boldsymbol{\mathbf{A}} $ 都应当被视为作用于 Hilbert 空间上的算符，这个 Hilbert 空间中的 “矢量” 是如同 $\phi,\chi$ 那样的二分量波函数。因此算符间不一定满足交换律，所以右边的两个算符的 “叉积” 不为 $0$。计算可得
\begin{equation}
\begin{aligned}
\left[ \boldsymbol{\mathbf{\sigma}} \cdot (\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} ) \right] ^2
&=(\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} )^2+i \boldsymbol{\mathbf{\sigma}} \cdot[-\frac{e}{c}(\hat{ \boldsymbol{\mathbf{P}} }\times \boldsymbol{\mathbf{A}} - \boldsymbol{\mathbf{A}} \times \hat{ \boldsymbol{\mathbf{P}} })]\\
&=(\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} )^2+i \boldsymbol{\mathbf{\sigma}} \cdot[\frac{e}{c}i\hbar \nabla \times \boldsymbol{\mathbf{A}} ]\\
&=(\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} )^2-\frac{e}{c}\hbar \boldsymbol{\mathbf{\sigma}} \cdot \boldsymbol{\mathbf{B}} ~.
\end{aligned}
\end{equation}
于是有泡利方程
\begin{equation}
\left[(\hat H-e\phi)-\frac{1}{2m}(\hat{ \boldsymbol{\mathbf{P}} }-\frac{e}{c} \boldsymbol{\mathbf{A}} )^2+\frac{e}{2mc}\hbar \boldsymbol{\mathbf{\sigma}} \cdot \boldsymbol{\mathbf{B}} \right] \varphi=0~,
\end{equation}
或者写为
\begin{equation}
i\hbar \frac{\partial}{\partial{t}} \varphi = \left[\frac{(-i\hbar\nabla - (e/c) \boldsymbol{\mathbf{A}} )^2}{2m}+e\phi-\frac{e}{mc} \boldsymbol{\mathbf{S}} \cdot \boldsymbol{\mathbf{B}} \right] \varphi~,
\end{equation}
$ \boldsymbol{\mathbf{S}} $ 为自旋算符。注意到自旋 $1/2$ 的粒子所满足的泡利方程比自旋为 $0$ 的粒子在电磁场中的薛定谔方程要多出一项 $- \boldsymbol{\mathbf{\mu}} \cdot \boldsymbol{\mathbf{B}} $，其中 $ \boldsymbol{\mathbf{\mu}} =e \boldsymbol{\mathbf{S}} /m$ 为粒子的自旋磁矩。
1. ^ 参考：维基百科
[1] ^ 刘绝平. 量子力学 北京高等教育出版社, 2012.8