晶格平移

晶格平移是一种离散对称操作,在固体物理学中具有极其重要的应用

一维周期势

单粒子在晶格势场中的哈密顿量为

H=p22m+V(q),V(q+a)=V(q)H=\frac{p^2}{2m}+V(q),\quad V(q+a)=V(q)

定义平移算符 T(na)=e−iℏp(na)T(na)=e^{-\frac{i}{\hbar}p(na)},其对势能的作用为

T†(na)V(q)T(na)=eiℏp(na)V(q)e−iℏp(na)=V(q+na)=V(q)T^\dagger(na)V(q)T(na)=e^{\frac{i}{\hbar}p(na)}V(q)e^{-\frac{i}{\hbar}p(na)}=V(q+na)=V(q)

因此有

[T(na),V(q)]=0⇒[T(na),H]=0[T(na),V(q)]=0\Rightarrow[T(na),H]=0

这说明HH 和 T(na)T(na) 可以同时对角化。又因为 T(na)T(na) 是幺正算符,其本征值的模为 11,那么就可以设其本征值的形式为 t(na)=eiθ(na)t(na)=e^{i\theta(na)},其中 θ\theta 为 nana 的函数。

定义 ∣m⟩|m\rangle 为粒子局域在第 mm 个格点上的态

格点处的态

那么有

T(na)∣m⟩=∣m+n⟩T(na)|m\rangle=|m+n\rangle

这说明局域态 {m}\{m\} 不是 T(na)T(\mathbf{n}\mathbf{a}) 的本征态。

H和T(na)的共同本征态

我们做紧束缚近似:将单体问题形象地看成是粒子跳跃问题,即粒子只能在相邻的格点间跳跃。设粒子在第 jj 个格点的能量为 ε0\varepsilon_0,在相邻格点间跳跃的能量为 −J-J,则哈密顿量可写成

H=ε0∑j=−∞∞∣j⟩⟨j∣−J∑j=−∞∞(∣j⟩⟨j+1∣+∣j+1⟩⟨j∣),J>0H=\varepsilon_0\sum_{j=-\infty}^\infty|j\rangle\langle j|-J\sum_{j=-\infty}^\infty(|j\rangle\langle j+1|+|j+1\rangle\langle j|),\quad J>0

式中第一项为在位能,第二项为最近邻粒子跃迁强度,粒子只能在相邻格点间跳跃。我们现在考虑还是无穷长的晶格,所以具有平移不变性 T†(na)HT(na)=HT^\dagger(na)HT(na)=H,若是有限长晶格则不具有平移不变性。

现将紧束缚的哈密顿量作用在 ∣n⟩\ket{n} 态上,有

H∣n⟩=ε0∑j∣j⟩⟨j∣n⟩−J∑j(∣j⟩⟨j+1∣n⟩+∣j+1⟩⟨j∣n⟩)=ε0∑j∣j⟩δjn−J∑j(∣j⟩δj+1,n+∣j+1⟩δjn)=ε0∣n⟩−J(∣n−1⟩+∣n+1⟩)\begin{aligned} H|n\rangle&=\varepsilon_0\sum_j|j\rangle\langle j|n\rangle-J\sum_j(|j\rangle\langle j+1|n\rangle+|j+1\rangle\langle j|n\rangle) \\[1em] &=\varepsilon_0\sum_j|j\rangle\delta_{jn}-J\sum_j(|j\rangle\delta_{j+1,n}+|j+1\rangle\delta_{jn}) \\[1em] &=\varepsilon_0|n\rangle-J(|n-1\rangle+|n+1\rangle) \end{aligned}

所以 ∣n⟩\ket{n} 也不是 HH 的本征态。但我们要的是 HH 和 T(na)T(na) 的共同本征态,怎么办?线性展开!

构建所有局域态的线性组合

∣k⟩=∑n=−∞∞Cn(k)∣n⟩|k\rangle=\sum_{n=-\infty}^\infty C_n(k)|n\rangle

其中 kk 是波矢,用于标记一个特定的本征态。我们刚刚已经提到要求

H∣k⟩=Ek∣k⟩T(a)∣k⟩=eiθk(a)∣k⟩H|k\rangle=E_k|k\rangle \qquad T(a)|k\rangle=e^{i\theta_k(a)}|k\rangle

利用 HH 和 T(a)T(a) 作用在 ∣n⟩|n\rangle 上的结果,经过整理可以分别得到 Cn(k)C_n(k) 需满足的条件为

ε0Cn(k)−J[Cn+1(k)+Cn−1(k)]=EkCn(k)\varepsilon_0C_n(k)-J[C_{n+1}(k)+C_{n-1}(k)]=E_kC_n(k)

Cn−1(k)=eiθk(a)Cn(k)C_{n-1}(k)=e^{i\theta_k(a)}C_n(k)

上两式中第二式的解显然是

Cn(k)=Ae−inθk(a)C_n(k)=Ae^{-in\theta_k(a)}

我们要求 θk(a)\theta_k(a) 是无量纲的,那么最简单的选择就是 θk(a)=ka\theta_k(a)=ka,其中 kk 为波数。再结合第一个方程,可以得到能谱为

Ek=ε0−J(e−ika+eika)=ε0−2Jcos⁡kaE_k=\varepsilon_0-J\left(e^{-ika}+e^{ika}\right)=\varepsilon_0-2J\cos ka

  • 本征态 ∣k⟩|k\rangle 是局域态 {∣j⟩}\{|j\rangle\} 的傅里叶变换,称为布洛赫态,由波数 kk 标记。

  • 这就相当于从局域态的表象经幺正变换到能量本征态的表象,又因为希尔伯特空间的维数不因表象变换而变化,所以 kk 也是连续变量,理论上可以取任意实数值,但是由于 kk 出现在 cos⁡ka\cos ka 中,本征能量连续分布在 ε0−2J\varepsilon_0 - 2J 到 ε0+2J\varepsilon_0 + 2J 之间,所以 kk 没必要取全部实数值,这就是布里渊区的概念

  • kk 通常取在 [−πa,πa)[-\dfrac{\pi}{a}, \dfrac{\pi}{a}) 范围内,称为第一布里渊区。

能谱的宽度正比于 JJ,当 J=0J=0 时,晶格势垒变得无限高,粒子不会跳跃,哈密顿量只剩下在位能

J=0势垒无限高

H0=ε0∑j=−∞∞∣j⟩⟨j∣,H0∣j⟩=ε0∣j⟩H_0=\varepsilon_0\sum_{j=-\infty}^\infty|j\rangle\langle j|,\qquad H_0|j\rangle=\varepsilon_0|j\rangle

每一个局域态都是能量本征态,本征能量为在位能,此时本征态为无穷维简并。

一维晶格能带

最近邻跃迁 JJ 解除了无穷重简并,形成了一个连续的能带。

J=0J=0 时的坐标表象下的 mm 和能量表象下的 kk 态的本征能量虽然都是 ε0\varepsilon_0,但它们并不相同,前者是局域态,后者是完全波动的布洛赫态。

布洛赫定理

布洛赫态 ∣k⟩|k\rangle 的波函数为:⟨q′∣k⟩=A∑ne−ikna⟨q′∣n⟩\langle q^{\prime}|k\rangle=A\sum_ne^{-ikna}\langle q^{\prime}|n\rangle

替换为 q′+aq'+a 仍然成立:⟨q′+a∣k⟩=A∑ne−ikna⟨q′+a∣n⟩\langle q^{\prime}+a|k\rangle=A\sum_ne^{-ikna}\langle q^{\prime}+a|n\rangle

又因为 Wannier function\text{Wannier function},∣n⟩\ket{n} 在 q′+maq'+ma 处的波函数和 ∣n−m⟩\ket{n-m} 在 q′q' 处的波函数相同,从数学上来讲就是

⟨q′+ma∣n⟩=⟨q′∣T†(ma)∣n⟩=⟨q′∣T(−ma)∣n⟩=⟨q′∣n−m⟩\langle q^{\prime}+ma|n\rangle=\bra{q^{\prime}}T^\dagger(ma)|n\rangle=\bra{q^{\prime}}T(-ma)|n\rangle=\langle q^{\prime}|n-m\rangle

那么有

⟨q′+a∣k⟩=A∑ne−ikna⟨q′∣n−1⟩=A∑ne−ik(n+1)a⟨q′∣n⟩=e−ikaA∑ne−ikna⟨q′∣n⟩=e−ika⟨q′∣k⟩\langle q^{\prime}+a|k\rangle=A\sum_ne^{-ikna}\langle q^{\prime}|n-1\rangle=A\sum_ne^{-ik(n+1)a}\langle q^{\prime}|n\rangle \\[1em] =e^{-ika}A\sum_ne^{-ikna}\langle q^{\prime}|n\rangle=e^{-ika}\langle q^{\prime}|k\rangle

再在上式两边同时乘以 eik(q′+a)e^{ik\left(q^{\prime}+a\right)},得到

eik(q′+a)⟨q′+a∣k⟩=eik(q′+a)e−ika⟨q′∣k⟩=eikq′⟨q′∣k⟩e^{ik(q^{\prime}+a)}\langle q^{\prime}+a|k\rangle=e^{ik(q^{\prime}+a)}e^{-ika}\langle q^{\prime}|k\rangle=e^{ikq^{\prime}}\langle q^{\prime}|k\rangle

上两式对比说明,布洛赫态的波函数满足本是不满足周期性,但是当布洛赫态乘以因子 eikqe^{ikq} 后就满足周期性,所以我们定义
uk(q′)≡eikq′⟨q′∣k⟩u_k(q^{\prime})\equiv e^{ikq^{\prime}}\langle q^{\prime}|k\rangle,则有

uk(q′)=uk(q′+a)u_k(q^{\prime})=u_k(q^{\prime}+a)

那么波函数 ⟨q′∣k⟩≡ψk(q′)\langle q^{\prime}|k\rangle\equiv\psi_k(q^{\prime}) 有以下形式

ψk(q′)=e−ikq′uk(q′)\psi_k(q^{\prime})=e^{-ikq^{\prime}}u_k(q^{\prime})

这就是著名的布洛赫定理。

环上的有限尺寸问题

不止无限长晶格有平移对称性,有限长晶格在首尾相连成环后也有平移对称性。

我们取有限个格点 NN 个,首尾相连成环

H=ε0∑j=1N∣j⟩⟨j∣−J∑j=1N(∣j⟩⟨j+1∣+Н.с.),J>0H=\varepsilon_0\sum_{j=1}^N|j\rangle\langle j|-J\sum_{j=1}^N(|j\rangle\langle j+1|+\mathrm{Н.с.),}\qquad J>0

其中为简单起见假设 NN 为偶整数。指标 jj 从 11 取到 NN,因为假设世界是一个长度为 NN 的环,使得第 (N+m)(N + m) 个点就是第 mm 个点,也就是说满足周期性边界条件

TN(a)∣m⟩=∣N+m⟩=∣m⟩T^N(a)|m\rangle=|N+m\rangle=|m\rangle

可以证明,环上的有限元系统仍然满足平移不变性

T†(a)HT(a)=HT^\dagger(a)HT(a)=H

其布洛赫态的构造与无限长晶格类似,只能取到 NN,这时可以进行归一化

∣k⟩≡1NΣn=1Ne−ikna∣n⟩|k\rangle\equiv\frac{1}{\sqrt{N}}\Sigma_{n=1}^Ne^{-ikna}|n\rangle

我们还是要求 ∣k⟩|k\rangle 是 平移算符 TT 的本征态,但此时这一要求就会导致 kk 只能取离散值

T(a)∣k⟩=T(a)1N∑n=1Ne−ikna∣n⟩=1N∑n=1Ne−ikna∣n+1⟩=1N∑n=2N+1e−ik(n−1)a∣n⟩=1N∑n=2Ne−ik(n−1)a∣n⟩+1Ne−ikNa∣N+1⟩=eika1N∑n=2Ne−ikna∣n⟩+1Ne−ikNa∣1⟩\begin{aligned} T(a)|k\rangle &= T(a) \frac{1}{\sqrt{N}} \sum_{n=1}^N e^{-ikna} |n\rangle = \frac{1}{\sqrt{N}} \sum_{n=1}^N e^{-ikna} |n+1\rangle \\[1em] &= \frac{1}{\sqrt{N}} \sum_{n=2}^{\textcolor{red}{N+1}} e^{-ik(n-1)a} |n\rangle = \frac{1}{\sqrt{N}} \sum_{n=2}^{\textcolor{red}{N}} e^{-ik(n-1)a} |n\rangle + \textcolor{red}{\frac{1}{\sqrt{N}} e^{-ikNa} |N+1\rangle} \\[1em] &= \textcolor{orange}{e^{ika} \frac{1}{\sqrt{N}} \sum_{n=2}^N e^{-ikna} |n\rangle} + \textcolor{teal}{\frac{1}{\sqrt{N}} e^{-ikNa} |1\rangle} \end{aligned}

又由于

T(a)∣k⟩=eika∣k⟩=eika1N∑n=2Ne−ikna∣n⟩+eika1Ne−ika∣1⟩T(a)|k\rangle=e^{ika}|k\rangle=e^{ika}\frac{1}{\sqrt{N}}\sum_{n=2}^{N}e^{-ikna}|n\rangle+e^{ika}\frac{1}{\sqrt{N}}e^{-ika}|1\rangle

对比上两式得到

e−ikNa=1⟹k=2πmNa,m∈Ze^{-ikNa}=1\quad\Longrightarrow\quad k=\frac{2\pi m}{Na},m\in\mathbb{Z}

也就是说 kk 只能取离散值,但我们可以将 kk 限制在第一布里渊区内

k∈KN=even={−πa,−πa+2πNa,⋯ ,−2πNa,0,2πNa,⋯ ,πa−2πNa}≡{k1,k2,⋯ ,kN}k\in K_{N=\mathrm{even}}=\left\{-\frac{\pi}{a},-\frac{\pi}{a}+\frac{2\pi}{Na},\cdots,-\frac{2\pi}{Na},0,\frac{2\pi}{Na},\cdots,\frac{\pi}{a}-\frac{2\pi}{Na}\right\}\equiv\{k_1,k_2,\cdots,k_N\}

kα≡−πa+2(α−1)πNak_\alpha\equiv-\frac{\pi}{a}+\frac{2(\alpha-1)\pi}{Na}

因为有限情况下态的个数是有限的,所以 kk 的个数也应是有限的。


这里给出两个有用的等式

  1. kk 空间正交性:

    ⟨kα∣kβ⟩=1N∑n=1Nei(kα−kβ)na=δαβ,α,β=1,...,N\langle k_\alpha | k_\beta \rangle = \frac{1}{N} \sum_{n=1}^{N} e^{i(k_\alpha - k_\beta)na} = \delta_{\alpha\beta}, \quad \alpha, \beta = 1, ..., N

  2. 实空间正交性:(n,mn,m 表示格点编号)

    1N∑α=1Neikα(n−m)a=δnm\frac{1}{N} \sum_{\alpha=1}^{N} e^{ik_\alpha(n-m)a} = \delta_{nm}

这两个式子将会用来证明有限晶格下的哈密顿量在 kk 空间的对角化,为此我们先找到布洛赫态 ∣k⟩\ket{k} 到局域态 ∣n⟩\ket{n} 的逆变换

∣n⟩≡1NΣα=1Neikαna∣kα⟩|n\rangle\equiv\frac{1}{\sqrt{N}}\Sigma_{\alpha=1}^Ne^{ik_\alpha\mathrm{n}a}|k_\alpha\rangle

kk 空间和实空间的傅里叶变换

由于哈密顿量最开始是在局域态下给出的,那么我们将上式直接代入哈密顿量,先看跃迁项

∑j=1N∣j⟩⟨j+n∣=∑j=1N1N∑α=1Neikαja∣kα⟩1N∑β=1Ne−ikβ(j+n)a⟨kβ∣=∑α,βe−ikβna[1N∑j=1Nei(kα−kβ)ja]∣kα⟩⟨kβ∣=∑α,βe−ikβnaδαβ∣kα⟩⟨kβ∣=∑αe−ikαna∣kα⟩⟨kα∣\begin{aligned} \sum_{j=1}^N|j\rangle\langle j+n|&=\sum_{j=1}^N\frac{1}{\sqrt{N}}\sum_{\alpha=1}^Ne^{ik_\alpha ja}|k_\alpha\rangle\frac{1}{\sqrt{N}}\sum_{\beta=1}^Ne^{-ik_\beta(j+n)a}\langle k_\beta| \\[1em] &=\sum_{\alpha,\beta}e^{-ik_\beta na}[\frac{1}{N}\sum_{j=1}^Ne^{i(k_\alpha-k_\beta)ja}]|k_\alpha\rangle\langle k_\beta|=\sum_{\alpha,\beta}e^{-ik_\beta na}\delta_{\alpha\beta}|k_\alpha\rangle\langle k_\beta| \\[1em] &=\sum_\alpha e^{-ik_\alpha\mathrm{n}a}|k_\alpha\rangle\langle k_\alpha| \end{aligned}

上式将跃迁项 ∑∣j⟩⟨j+n∣\sum |j\rangle\langle j+n| 变换到 kk 空间,且是对角的。利用上式直接计算紧束缚模型的哈密顿量

H=ε0∑j=1N∣j⟩⟨j∣−J∑j=1N(∣j⟩⟨j+1∣+H.c.)=ε0∑αe−ikα⋅0⋅a∣kα⟩⟨kα∣−J(∑αe−ikα⋅1⋅a∣kα⟩⟨kα∣+H.c.)=∑α(ε0−2Jcos⁡kαa)∣kα⟩⟨kα∣\begin{aligned} \mathrm{H} & =\varepsilon_0\sum_{j=1}^N|j\rangle\langle j|-J\sum_{j=1}^N(|j\rangle\langle j+1|+\mathrm{H.c.}) \\[1em] & =\varepsilon_0\sum_\alpha e^{-ik_\alpha\cdot0\cdot a}|k_\alpha\rangle\langle k_\alpha|-J(\sum_\alpha e^{-ik_\alpha\cdot1\cdot a}|k_\alpha\rangle\langle k_\alpha|+\mathrm{H.c.}) \\[1em] & =\sum_\alpha(\varepsilon_0-2J\cos k_\alpha a)|k_\alpha\rangle\langle k_\alpha| \end{aligned}

可以看出有限尺寸晶格的哈密顿量在布洛赫态下也是对角化的,其本征能量形式与无限长晶格时相同,只不过此时的 kk 只能取离散值。

空间反演与宇称

空间反演

不同维度坐标的空间反演和镜像分别为

空间反演 镜像

三维空间的反演和镜像关系为

空间反演=绕 −y^ 轴旋转 180∘×关于 x0z 平面的镜面反射\begin{aligned} \text{空间反演} &= \textcolor{green}{\text{绕 } -\hat{y} \text{ 轴旋转 } 180^\circ} \times \textcolor{brown}{\text{关于 } x0z \text{ 平面的镜面反射}} \end{aligned}

反演和镜像的关系
  • 二维的空间反演不改变坐标轴的手性,而三维的空间反演会改变坐标轴的手性。
  • 镜像操作会改变坐标轴的手性
  • 相同的手性可通过旋转联系起来

我们来看一个反直觉的事实:电流方向在空间反演下不改变,而在镜像下改变,这和我们对矢量的理解是相反的,下面我们将会讲到磁场不是真正的矢量。

赝矢量的特性

磁场:由电流(或运动电荷)产生 ⇒\Rightarrow 在空间反演下不变。
极矢量:在空间反演下变号的矢量(例如 x,p,E\boldsymbol{x}, \boldsymbol{p}, \boldsymbol{E})。
轴矢量:在空间反演下不变号的矢量(例如 S,L,x×p,B\boldsymbol{S}, \boldsymbol{L}, \boldsymbol{x} \times \boldsymbol{p}, \boldsymbol{B})。
赝标量:在旋转下不变但在空间反演下变号的量,由轴矢量和极矢量的乘积组合而成(例如 E⋅B\boldsymbol{E} \cdot \boldsymbol{B} 或 S⋅P\boldsymbol{S} \cdot \boldsymbol{P})。

量子力学中的空间反演

在量子力学中,空间反演是通过一个酉算符 π\pi,称为宇称算符来实现的,此时 π\pi 不是圆周率而是一个算符。对于给定态 ∣ψ⟩|\psi\rangle,空间反演态是 π∣ψ⟩\pi|\psi\rangle。

在量子力学中我们无法得到其确定的坐标,但是可以谈论其平均值,我们要求平均值在空间反演下变号:

(⟨ψ∣π†)x(π∣ψ⟩)=−⟨ψ∣x∣ψ⟩  ⟹  π†xπ=−x(\langle\psi|\pi^\dagger)\boldsymbol{x}(\pi|\psi\rangle) = -\langle\psi|\boldsymbol{x}|\psi\rangle \implies \pi^\dagger \boldsymbol{x} \pi = -\boldsymbol{x}

π\pi 是幺正的,那么 π−1=π†  ⟹  π−1xπ=−x  ⟹  \pi^{-1} = \pi^\dagger \implies \pi^{-1}\boldsymbol{x}\pi = -\boldsymbol{x} \implies {π,x}=0\{\pi, \boldsymbol{x}\} = 0,即宇称算符与位置算符反对易。

连续两次空间反演完全不产生变化:ππ∣ψ⟩=eiθ∣ψ⟩\pi\pi|\psi\rangle = e^{i\theta}|\psi\rangle。为方便起见通常选取 θ=0\theta = 0,即 π=π−1\pi = \pi^{-1}。因此,π\pi 既是幺正的也是厄米的:π=π−1=π†\pi = \pi^{-1} = \pi^\dagger。既然 π2=1\pi^2 = 1,那么算符 π\pi 就满足一个代数方程,其本征值只能是 +1+1 或 −1-1。

π∣x′⟩\pi|\boldsymbol{x}'\rangle 是 x\boldsymbol{x} 的本征值为 −x′-\boldsymbol{x}' 的本征态,因此 π∣x′⟩=eiθ∣−x′⟩\pi|\boldsymbol{x}'\rangle = e^{i\theta}|-\boldsymbol{x}'\rangle,再次按惯例选取 θ=0\theta = 0,我们得到 π∣x′⟩=∣−x′⟩\pi|\boldsymbol{x}'\rangle = |-\boldsymbol{x}'\rangle。


我们考虑 π\pi 算符对波函数的作用。令 π^\hat{\pi} 为 π\pi 在 x\boldsymbol{x} 基底下的表示:⟨x′∣π∣ψ⟩=π^⟨x′∣ψ⟩=π^ψ(x′)\langle \boldsymbol{x}'|\pi|\psi\rangle = \hat{\pi}\langle \boldsymbol{x}'|\psi\rangle = \hat{\pi}\psi(\boldsymbol{x}'),由 π\pi 的厄米性 ⟨x′∣π=⟨−x′∣\langle \boldsymbol{x}'|\pi = \langle -\boldsymbol{x}'|

⟨x′∣π∣ψ⟩=⟨−x′∣ψ⟩=ψ(−x′)\langle \boldsymbol{x}'|\pi|\psi\rangle = \langle -\boldsymbol{x}'|\psi\rangle = \psi(-\boldsymbol{x}')

因此 π\pi 对波函数的作用就是改变坐标的符号

π^ψ(x′)=ψ(−x′)\hat{\pi}\psi(\boldsymbol{x}') = \psi(-\boldsymbol{x}')

所有 x′\boldsymbol{x}' 的偶(奇)函数都是 π^\hat{\pi} 的本征函数,其本征值为 ±1\pm 1:

π^ψe/o(x′)=±ψe/o(x′)\hat{\pi}\psi_{e/o}(\boldsymbol{x}') = \pm\psi_{e/o}(\boldsymbol{x}')

假设有两个态 ∣ϕ+⟩|\phi_+\rangle 和 ∣ϕ−⟩|\phi_-\rangle,分别是 π\pi 的两个本征态 π∣ϕ±⟩=±∣ϕ±⟩\pi|\phi_\pm\rangle = \pm|\phi_\pm\rangle,于是

⟨x′∣π∣ϕ±⟩=±⟨x′∣ϕ±⟩  ⟹  ϕ±(−x′)=±ϕ±(x′)\langle \boldsymbol{x}'|\pi|\phi_\pm\rangle = \pm\langle \boldsymbol{x}'|\phi_\pm\rangle \implies \phi_\pm(-\boldsymbol{x}') = \pm\phi_\pm(\boldsymbol{x}')

这意味着本征值为 11 的态的波函数是偶函数,而本征值为 −1-1 的态的波函数是奇函数。偶(奇)波函数或其对应的态被称为具有宇称 +1(−1)+1(-1),或称为偶(奇)宇称。

本征值为 ±1\pm1 的态 ⟺\Longleftrightarrow 偶(奇)宇称态

可以证明 π\pi 对动量算符的作用同样为

π†pπ=−p\pi^\dagger \boldsymbol{p} \pi = -\boldsymbol{p}

那么角动量 L=x×pL=x\times p 在 π\pi 作用是不变的: π†Lπ=L\pi^\dagger L\pi=L,我们希望相应的自旋角动量 SS 也不变号 π†Sπ=S\pi^\dagger S\pi=S,因此总角动量 J=L+SJ=L+S 也不变号 π†Jπ=J\pi^\dagger J\pi=J。

可以证明平移算符和宇称算符不对易,而转动算符和宇称算符对易。

假设系统的哈密顿量在空间反演下不变,即 π†Hπ=H\pi^\dagger H \pi = H,且 ∣n⟩|n\rangle 是 HH 的一个非简并本征态,其本征值为 EnE_n:H∣n⟩=En∣n⟩H|n\rangle = E_n|n\rangle,那么 ∣n⟩|n\rangle 也是 π\pi 的本征态。

对易算符具有相同的本征态

宇称不守恒

如果 [π,H]=0[\pi, H] = 0,那么我们就说宇称本征值 ±1\pm 1 是一个守恒量。长期以来,人们一直认为自然界的基本定律在空间反演(或镜像)下是不变的,简单地说就是自然界不区分左手性和右手性。“吴氏实验”(1956):表明自然界不遵守这种对称性。

放射性原子核 2760Co{}_{27}^{60}\text{Co} 发生 β\beta 衰变:

2760Co→2860Ni+e−+νˉe+2γ{}_{27}^{60}\text{Co} \to {}_{28}^{60}\text{Ni} + e^- + \bar{\nu}_e + 2\gamma

2760Co{}_{27}^{60}\text{Co} 原子核具有非零自旋 S\boldsymbol{S} 和磁矩,因此在低温下可以通过磁场使其自旋定向排列。

宇称不守恒实验

图中左边的情况是预期中宇称守恒的情况,在镜像世界中,电子应该更倾向于沿着自旋的正方向发射;而右边的情况是实际实验结果,电子的发射方向与原子核自旋方向存在明确的关联关系:电子始终更倾向于沿着自旋的反方向发射。这一现象表明,在 β\beta 衰变过程中,宇称并不是守恒的。

想象一下你骑一辆自行车经过一片水洼,往车轮两边溅出的水不是一样多的,这是一个非常违背直觉的现象,人们原本认为左右是相对的、不可区分的,但是现在可以利用多少来区分了。

我们可以做定性分析:令 S\boldsymbol{S} 为核自旋算符,P\boldsymbol{P} 为电子动量算符。对于物理态 ∣ψ⟩|\psi\rangle(在实验室中),观察到 ⟨ψ∣S⋅P∣ψ⟩<0\langle\psi|\boldsymbol{S} \cdot \boldsymbol{P}|\psi\rangle < 0。然而,对于反演态 ∣ψ′⟩=π∣ψ⟩|\psi'\rangle = \pi|\psi\rangle,我们有

⟨ψ′∣S⋅P∣ψ′⟩=⟨ψ∣π†S⋅Pπ∣ψ⟩=−⟨ψ∣S⋅P∣ψ⟩>0\langle\psi'|\boldsymbol{S} \cdot \boldsymbol{P}|\psi'\rangle = \langle\psi|\pi^\dagger \boldsymbol{S} \cdot \boldsymbol{P}\pi|\psi\rangle = -\langle\psi|\boldsymbol{S} \cdot \boldsymbol{P}|\psi\rangle > 0

现在,假设 β\beta 衰变的哈密顿量与 π\pi 对易 πH=Hπ\pi H = H\pi,这就是宇称守恒,那么

H∣ψ⟩=E∣ψ⟩  ⟹  Hπ∣ψ⟩=πH∣ψ⟩=Eπ∣ψ⟩  ⟹  H∣ψ′⟩=E∣ψ′⟩H|\psi\rangle = E|\psi\rangle \quad \implies \quad H\pi|\psi\rangle = \pi H|\psi\rangle = E\pi|\psi\rangle \quad \implies \quad H|\psi'\rangle = E|\psi'\rangle

这就会导致 ∣ψ⟩|\psi\rangle 和 ∣ψ′⟩|\psi'\rangle 要么是简并态,要么是同一个态(非简并)。但是它们又不可能是同一个态,否则 ⟨ψ∣S⋅P∣ψ⟩=⟨ψ′∣S⋅P∣ψ′⟩=0\langle\psi|\boldsymbol{S} \cdot \boldsymbol{P}|\psi\rangle =\langle\psi'|\boldsymbol{S} \cdot \boldsymbol{P}|\psi'\rangle= 0。因此,只有当 2760Co{}_{27}^{60}\text{Co} 原子核的自旋极化态是简并的,才有可能维持空间反演对称性,但是进一步的核结构理论不支持这一假设。

自发对称性破缺

假设哈密顿量 HH 具有某种对称性(例如平移、转动以及空间反演对称性),而其基态在该对称性下却是变化的,我们就说对称性发生了自发破缺。

这说明基态不是该对称算符的本征态,根据定理,对称性自发破缺下的基态必然是简并的。

这个现象的经典解释如下图:若是一个单阱谐振子,基态是唯一且对称的,但若是一个双阱势,则会处于 x=−ax = -a 或 x=ax = a,这两个互为空间反演。

对称性自发破缺经典解释

成立的条件是:存在超过 11 个基态,且这些态并不具有对称不变性,而是通过对称性操作互相映射。


自发对称性破缺的量子示例为:无限势垒模型

无限势垒模型

在基底 {∣a⟩,∣−a⟩}\{|a\rangle, |-a\rangle\} 下,哈密顿量和宇称算符分别为

H=(E000E0)π=(0110)H = \begin{pmatrix} E_0 & 0 \\[0.5em] 0 & E_0 \end{pmatrix} \qquad \pi = \begin{pmatrix} 0 & 1 \\[0.5em] 1 & 0 \end{pmatrix}

由于势垒无限高,两个能级之间不会发生跃迁,哈密顿量的非对角项为零。因为 [H,π]=0[H, \pi] = 0,所以 HH 和 π\pi 具有共同本征态,我们构造出两个本征态为

∣S⟩=12(∣a⟩+∣−a⟩),∣A⟩=12(∣a⟩−∣−a⟩)|S\rangle = \frac{1}{\sqrt{2}}(|a\rangle + |-a\rangle), |A\rangle = \frac{1}{\sqrt{2}}(|a\rangle - |-a\rangle)

  ⟹  π∣S⟩=∣S⟩,π∣A⟩=−∣A⟩\implies \pi|S\rangle = |S\rangle, \quad \pi|A\rangle = -|A\rangle

但是上式中的对称态和反对称态不是物理存在的。假设在测量过程中粒子被发现位于势垒的一侧,那么由于无限势垒的阻隔,绝无可能实现 ∣S/A⟩|S/A\rangle 态。粒子无法感知到另一个与其自身所在环境完全一样的“平行宇宙”的存在。在此类问题中,现实物理世界将不再是对称的,且 ∣S/A⟩|S/A\rangle 代表了无法实现的物理情形 —— 即对称性发生了自发破缺。


若是有限深势垒,那么哈密顿量中的非对角项不为零

H=(E0−V−VE0)=E0−Vσx  ⟹  {H∣S⟩=(E0−V)∣S⟩H∣A⟩=(E0+V)∣A⟩H = \begin{pmatrix} E_0 & -V \\[0.5em] -V & E_0 \end{pmatrix} = E_0 - V\sigma_x \implies \begin{cases} H|S\rangle = (E_0 - V)|S\rangle \\[0.5em] H|A\rangle = (E_0 + V)|A\rangle \end{cases}

此时简并解除,∣S⟩|S\rangle 是唯一的基态,所以有限势垒问题中不会发生自发对称性破缺,这就是通过隧穿实现的对称性恢复。

我们再来计算从 ∣a⟩|a\rangle 态演化至到 ∣−a⟩|-a\rangle 态所需的时间。如果我们从定域态 ∣a⟩=12(∣S⟩+∣A⟩)|a\rangle = \dfrac{1}{\sqrt{2}}(|S\rangle + |A\rangle) 开始演化,那么经过时间 tt 后的态为

∣ψ(t)⟩=H∣a⟩=12(e−iℏ(E0−V)t∣S⟩+e−iℏ(E0+V)t∣A⟩)=12e−iℏ(E0−V)t(∣S⟩+e−iℏ2Vt∣A⟩)|\psi(t)\rangle =H\ket{a}= \frac{1}{\sqrt{2}}(e^{-\frac{i}{\hbar}(E_0-V)t}|S\rangle + e^{-\frac{i}{\hbar}(E_0+V)t}|A\rangle) = \frac{1}{\sqrt{2}}e^{-\frac{i}{\hbar}(E_0-V)t}(|S\rangle + e^{-\frac{i}{\hbar}2Vt}|A\rangle)

该态会以频率 ω=2V/ℏ\omega = 2V/\hbar 在 ∣a⟩|a\rangle 和 ∣−a⟩|-a\rangle 之间来回振荡。若是势垒有限高但是又足够大,那么 V→0⇒T=πℏ2V→∞V \to 0 \Rightarrow T = \dfrac{\pi\hbar}{2V} \to \infty,势垒越高,粒子从 ∣a⟩|a\rangle 态隧穿到 ∣−a⟩|-a\rangle 态所需的时间就越长,在势垒无限高时,粒子永远无法从 ∣a⟩|a\rangle 态隧穿到 ∣−a⟩|-a\rangle 态,从而实现了自发对称性破缺。

因此,我们观察到的宇称不对称性,例如人体心脏在左边,可以通过非对称的初始状态来解释,也就是一开始就处于 ∣−a⟩\ket{-a} 或 ∣a⟩\ket{a} 态,而不需要用宇称不守恒来解释。

时间反演

时间反演更准确的说法应该是运动反演,但时间反演这个表述已经被广泛接受,所以我们接下来都使用时间反演这个术语。

经典力学中的时间反演

若 x(t)\boldsymbol{x}(t) 是牛顿方程 mx¨(t)=−∇V(x(t))m\ddot{\boldsymbol{x}}(t) = -\nabla V(\boldsymbol{x}(t)) 的解(例如钟摆),那么 x(−t)\boldsymbol{x}(-t) 也应是一个解:

mdd(−t)dd(−t)x(−t)=−∇V(x(−t))  ⟹  mddtddtx(−t)=−∇V(x(−t))m \frac{d}{d(-t)}\frac{d}{d(-t)}\boldsymbol{x}(-t) = -\nabla V(\boldsymbol{x}(-t)) \quad \implies \quad m \frac{d}{dt}\frac{d}{dt}\boldsymbol{x}(-t) = -\nabla V(\boldsymbol{x}(-t))

若是满足上式则具有时间反演对称性。时间反演对称性破坏的两个例子:

  1. 引入摩擦力:mx¨(t)=−∇V(x(t))−γx˙m\ddot{\boldsymbol{x}}(t) = -\nabla V(\boldsymbol{x}(t)) - \gamma \dot{\boldsymbol{x}},其中包含速度 x˙\dot{\boldsymbol{x}},若是将运动过程(例如钟摆运动)拍成视频,倒放和正放是可以区分的。

  2. mx¨(t)=qx˙(t)×Bm\ddot{\boldsymbol{x}}(t) = q\dot{\boldsymbol{x}}(t) \times \boldsymbol{B},其中也包含速度 x˙\dot{\boldsymbol{x}},如下图所示,正向运动和反向运动分别代表时间向前和时间向后的过程,可以看出两个不同方向的电子会得到两个不同的轨迹,因此不存在时间反演对称性

    磁场中的时间反演

    但是如果你同时反转 B\boldsymbol{B},那么 x(−t)\boldsymbol{x}(-t) 仍然是下列方程的一个解:

    mx¨=qx˙×(−B)m\ddot{\boldsymbol{x}} = q\dot{\boldsymbol{x}} \times (-\boldsymbol{B})

    我们说磁场的存在破坏了时间反演对称性,那是因为我们将磁场视为独立于所考虑物理系统的背景场。实际上,磁场是由运动电荷产生的,而时间反演下电荷会反向运动,因此磁场也会随之反向,从而恢复时间反演对称性。

量子力学中的时间反演

量子力学中的时间反演和经典中有所不同。量子中使用波函数描述状态而不是 x,p\boldsymbol{x},\boldsymbol{p},使用算符描述力学量而不是 x,p\boldsymbol{x},\boldsymbol{p},经典和量子之间通过平均值联系起来,量子力学的平均值应回到经典力学的结论。在这里我们要求量子力学平均值的时间反演对称性与经典力学一致。

定义时间反演态(或运动反演态)为 T∣ψ⟩\mathcal{T}|\psi\rangle。如果时间反演对称性能够保持,如下图所示,我们希望演化后的时间反演,再经过相同时间的演化,等于初始态的时间反演态

量子态的时间反演

T∣ψ⟩=e−iℏHt(Te−iℏHt∣ψ⟩)\mathcal{T}|\psi\rangle = e^{-\frac{i}{\hbar}Ht}(\mathcal{T}e^{-\frac{i}{\hbar}Ht}|\psi\rangle)

  ⟹  eiℏHtT∣ψ⟩=Te−iℏHt∣ψ⟩\implies e^{\frac{i}{\hbar}Ht}\mathcal{T}|\psi\rangle = \mathcal{T}e^{-\frac{i}{\hbar}Ht}|\psi\rangle

由于时间 tt 是任意的,我们先考虑无穷小时间微元 t=dtt = dt,则有

(1+iℏHdt)T∣ψ⟩=T(1−iℏHdt)∣ψ⟩\left(1 + \frac{i}{\hbar}Hdt\right)\mathcal{T}|\psi\rangle = \mathcal{T}\left(1 - \frac{i}{\hbar}Hdt\right)|\psi\rangle

最后得到 T\mathcal{T} 需要满足的条件

iHT=−TiHiH\mathcal{T} = -\mathcal{T}iH

你可能会问,上式中的常数 ii 为什么不能消去?不能!因为 T\mathcal{T} 不是线性算符,实际上,T\mathcal{T} 对 ii 也有作用。

那么你又问了,时间反演算符 T\mathcal{T} 可以是幺正算符吗?我们可以先假设 T\mathcal{T} 是幺正算符,那么有

{HT=−THH∣n⟩=En∣n⟩\begin{cases} H\mathcal{T} = -\mathcal{T}H \\ H|n\rangle = E_n|n\rangle \end{cases}

  ⟹  HT∣n⟩=−TH∣n⟩=(−En)T∣n⟩\implies H\mathcal{T}|n\rangle = -\mathcal{T}H|n\rangle = (-E_n)\mathcal{T}|n\rangle

即反演变换态 T∣n⟩\mathcal{T}|n\rangle 是 HH 本征能量为 −En-E_n的本征态。但这即使在自由粒子的基本情况下也不能成立,因为我们要求能谱必须要存在下限。这说明,时间反演不是幺正算符,我们不允许消去 iHT=−TiHiH\mathcal{T} = -\mathcal{T}iH 中的 ii。实际上,时间反演算符与共轭操作有关。

反线性算符

反线性算符 AA 满足以下关系

A(c1∣ψ1⟩+c2∣ψ2⟩)=c1∗A∣ψ1⟩+c2∗A∣ψ2⟩A(c_1|\psi_1\rangle + c_2|\psi_2\rangle) = c_1^* A|\psi_1\rangle + c_2^* A|\psi_2\rangle

即对常数取复共轭。两个反线性算符的乘积是线性的

A2A1c∣ψ⟩=A2c∗A1∣ψ⟩=cA2A1∣ψ⟩A_2 A_1 c|\psi\rangle = A_2 c^* A_1 |\psi\rangle = c A_2 A_1 |\psi\rangle

我们一般约定:反线性算符只向右作用,而不向左作用。这导致与线性算符不同,反线性算符无法像线性算符那样完美地适配狄拉克符号。

与线性算符不同,反线性操作依赖于基底的选择。例如有 A^[ci∣ψ⟩]=−c∗iA^∣ψ⟩\hat{A}[ci|\psi\rangle] = -c^*i \hat{A}|\psi\rangle,但是如果我们令 ∣ψ′⟩=i∣ψ⟩|\psi'\rangle = i|\psi\rangle,那么 A^[c∣ψ′⟩]=c∗A^∣ψ′⟩=c∗iA^∣ψ⟩\hat{A}[c|\psi'\rangle] = c^* \hat{A}|\psi'\rangle= c^* i \hat{A}|\psi\rangle。因为基矢选取的不同,同样的操作得到了不同的结果。复共轭算符是最简单的反线性算符,其作用就仅是简单的取共轭,我们来看它在不同基底下的作用

  1. 正交归一基底集 {∣n⟩}\{|n\rangle\}。对于任意态 ∣ψ⟩|\psi\rangle,在该 nn-基底下的复共轭算符 K(n)K_{(n)} 定义为:

    K(n)∣ψ⟩=K(n)∑n⟨n∣ψ⟩∣n⟩=∑n⟨ψ∣n⟩K(n)∣n⟩=∑n⟨ψ∣n⟩∣n⟩K_{(n)}|\psi\rangle = K_{(n)} \sum_n \langle n|\psi\rangle|n\rangle = \sum_n \langle \psi|n\rangle K_{(n)}|n\rangle = \sum_n \langle \psi|n\rangle|n\rangle

    其中 ∣n⟩|n\rangle 应被理解为基底列向量 (0,⋯ ,1,⋯ ,0)T(0, \cdots, 1, \cdots, 0)^T。由于概率幅取复共轭不影响概率,因此我们认为:

    K(n)∣n⟩=∣n⟩K_{(n)}|n\rangle = |n\rangle

  2. 正交归一基底矢量集 {∣ν⟩}\{|\nu\rangle\}。按照第一种基底的结论,在这个 ν\nu-基底下,复共轭算符的作用为 K(ν)∣ψ⟩=∑ν⟨ψ∣ν⟩∣ν⟩K_{(\nu)}|\psi\rangle = \sum_\nu \langle\psi|\nu\rangle|\nu\rangle

    那我又要问了:这两个不同基矢下的复共轭算符等价吗?分别作用在同一个态 ∣ψ⟩|\psi\rangle 上

    K(ν)∣ψ⟩=∑ν⟨ψ∣ν⟩∣ν⟩=∑ν⟨ψ∣ν⟩∑n⟨n∣ν⟩∣n⟩=∑n(∑ν⟨ψ∣ν⟩⟨n∣ν⟩)∣n⟩K(n)∣ψ⟩=∑n⟨ψ∣n⟩∣n⟩=∑n(∑ν⟨ψ∣ν⟩⟨ν∣n⟩)∣n⟩\begin{aligned} K_{(\nu)}|\psi\rangle &= \sum_\nu \langle\psi|\nu\rangle|\nu\rangle = \sum_\nu \langle\psi|\nu\rangle \sum_n \langle n|\nu\rangle|n\rangle = \sum_n \left(\sum_\nu \langle\psi|\nu\rangle \langle n|\nu\rangle\right)|n\rangle \\[1em] K_{(n)}|\psi\rangle &= \sum_n \langle\psi|n\rangle|n\rangle = \sum_n \left(\sum_\nu \langle\psi|\nu\rangle \langle \nu|n\rangle\right)|n\rangle \end{aligned}

    只有当内积 ⟨n∣ν⟩\langle n|\nu\rangle 对所有 nn 和 ν\nu 均为实数时,才有 K(n)=K(ν)K_{(n)} = K_{(\nu)},两组基下的复共轭算符才是等价的。

    例如,考虑自旋 1/2 粒子的态 ∣↑⟩y^=12(∣↑⟩+i∣↓⟩)|\uparrow\rangle_{\hat{y}} = \dfrac{1}{\sqrt{2}}(|\uparrow\rangle + i|\downarrow\rangle)。
    在 σz\sigma_z-基底下:∣↑⟩y^≐12(1i)  ⟹  K(z)∣↑⟩y^=12(1−i)=∣↓⟩y^|\uparrow\rangle_{\hat{y}} \doteq \dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ i \end{pmatrix} \implies K_{(z)}|\uparrow\rangle_{\hat{y}} = \dfrac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ -i \end{pmatrix} = |\downarrow\rangle_{\hat{y}}
    在 σy\sigma_y-基底下:∣↑⟩y^≐(10)  ⟹  K(y)∣↑⟩y^=(10)=∣↑⟩y^|\uparrow\rangle_{\hat{y}} \doteq \begin{pmatrix} 1 \\ 0 \end{pmatrix} \implies K_{(y)}|\uparrow\rangle_{\hat{y}} = \begin{pmatrix} 1 \\ 0 \end{pmatrix} = |\uparrow\rangle_{\hat{y}}

    σz\sigma_z 和 σy\sigma_y 基底下的复共轭算符不等价:
    K(z)≠K(y)K_{(z)} \neq K_{(y)}。但是可以验证,K(x)=K(z)K_{(x)} = K_{(z)},因为 σx\sigma_x 和 σz\sigma_z 基底之间的变换矩阵为实数矩阵。

反幺正算符

如果算符 AA 满足以下条件,则称为反幺正算符:

  1. 反线性

  2. 存在逆算符 A−1A^{-1}

  3. 不改变作用前后态的模:对于任意 ∣ψ⟩|\psi\rangle,∣∣∣ψ⟩∣∣=∣∣A∣ψ⟩∣∣|||\psi\rangle|| = ||A|\psi\rangle||。

我们可以证明,一个反幺正算符 AA 总是可以写为一个幺正算符和复共轭算符的乘积:

A=UKA = UK

由于复共轭算符 KK 是依赖于基底的,因此反幺正算符 AA 也是依赖于基底的。


我们已经约定反线性算符只向右作用,而不向左作用,那么若是真出现向左作用的情况该怎么办呢?

假设 AA 是一个反幺正算符,而 LL 是一个线性算符。定义 ∣ψ′⟩≡A∣ψ⟩|\psi'\rangle \equiv A|\psi\rangle 和 ∣ϕ′⟩≡A∣ϕ⟩|\phi'\rangle \equiv A|\phi\rangle,那么有

⟨ψ∣L∣ϕ⟩=⟨ϕ′∣AL†A−1∣ψ′⟩\langle\psi|L|\phi\rangle = \langle\phi'|AL^\dagger A^{-1}|\psi'\rangle

每当我们遇到复共轭算符时,我们总是选取某个基展开得到系数,就能消去复共轭算符,从而进行常规计算;而左矢则需要先行展开计算,再代入矩阵元计算式。

推论:令 L=1L = 1,我们得到:

⟨ψ∣ϕ⟩=⟨ϕ′∣ψ′⟩=⟨ψ′∣ϕ′⟩∗\langle\psi|\phi\rangle = \langle\phi'|\psi'\rangle = \langle\psi'|\phi'\rangle^*

这表明,两个时间反演前后态的内积是相互共轭的。

若 FF 是厄米算符,则:

⟨ψ∣F∣ϕ⟩=⟨ϕ′∣AFA−1∣ψ′⟩  ⟹  ⟨ψ∣F∣ψ⟩=⟨ψ′∣AFA−1∣ψ′⟩\langle\psi|F|\phi\rangle = \langle\phi'|AFA^{-1}|\psi'\rangle \quad \implies \quad \langle\psi|F|\psi\rangle = \langle\psi'|AFA^{-1}|\psi'\rangle

时间反演算符是反幺正的

在介绍反线性和反幺正算符后,我们就可以确定时间反演算符的性质了。我们知道,若是变换前后态的模不变,那么这个变换只能是幺正的或反幺正的。

如前所述,我们假定 T\mathcal{T} 是一个反幺正算符是合理的。因为平均值与经典相联系,所以对于时间反演态 ∣ψ′⟩=T∣ψ⟩|\psi'\rangle = \mathcal{T}|\psi\rangle,我们预期动量期望值会反号(也就是运动反演)⟨ψ′∣p∣ψ′⟩=−⟨ψ∣p∣ψ⟩\langle\psi'|\boldsymbol{p}|\psi'\rangle = -\langle\psi|\boldsymbol{p}|\psi\rangle。利用前述结论 ⟨ψ∣F∣ψ⟩=⟨ψ′∣AFA−1∣ψ′⟩\langle\psi|F|\psi\rangle = \langle\psi'|AFA^{-1}|\psi'\rangle,我们得到:

−⟨ψ′∣p∣ψ′⟩=⟨ψ∣p∣ψ⟩=⟨ψ′∣TpT−1∣ψ′⟩-\langle\psi'|\boldsymbol{p}|\psi'\rangle = \langle\psi|\boldsymbol{p}|\psi\rangle = \langle\psi'|\mathcal{T}\boldsymbol{p}\mathcal{T}^{-1}|\psi'\rangle

对比得到

TpT−1=−p\mathcal{T}\boldsymbol{p}\mathcal{T}^{-1} = -\boldsymbol{p}

同理,位置算符 x\boldsymbol{x} 在时间反演下应当保持不变,由此可导出轨道角动量 L\boldsymbol{L} 的变换性质:

TxT−1=x  ⟹  TLT−1=−L\mathcal{T}\boldsymbol{x}\mathcal{T}^{-1} = \boldsymbol{x} \implies \mathcal{T}\boldsymbol{L}\mathcal{T}^{-1} = -\boldsymbol{L}

对于自旋 S\boldsymbol{S} 和总角动量 J\boldsymbol{J},我们同样希望他们与轨道角动量具有相同的变换性质:

TST−1=−S,TJT−1=−J\mathcal{T}\boldsymbol{S}\mathcal{T}^{-1} = -\boldsymbol{S}, \quad \mathcal{T}\boldsymbol{J}\mathcal{T}^{-1} = -\boldsymbol{J}

现在考察对易关系 [xj,pj]=iℏ[x_j, p_j] = i\hbar 在时间反演下的变换:

TiℏT−1=T[xj,pj]T−1=TxjT−1TpjT−1−TpjT−1TxjT−1=xj(−pj)−(−pj)xj=−[xj,pj]=−iℏ  ⟹  TiT−1=−i\begin{aligned} \mathcal{T}i\hbar\mathcal{T}^{-1} &= \mathcal{T}[x_j, p_j]\mathcal{T}^{-1} = \mathcal{T}x_j\mathcal{T}^{-1}\mathcal{T}p_j\mathcal{T}^{-1} - \mathcal{T}p_j\mathcal{T}^{-1}\mathcal{T}x_j\mathcal{T}^{-1} \\[0.5em] &= x_j(-p_j) - (-p_j)x_j = -[x_j, p_j] = -i\hbar \end{aligned} \\[1em] \implies \mathcal{T}i\mathcal{T}^{-1} = -i

所以为了满足理论的自洽性,T\mathcal{T} 必须是一个反线性算符。

薛定谔方程的时间反演

我们之前提到过,tt 在量子力学中仅仅是一个参数,其不能被算符直接影响,在薛定谔方程两边同时作用时间反演算符 T\mathcal{T}

T[iℏ∂t∣ψ(t)⟩]=T[H∣ψ(t)⟩]  ⟹  −iℏ∂tT∣ψ(t)⟩=THT−1T∣ψ(t)⟩\mathcal{T}[i\hbar\partial_t|\psi(t)\rangle] = \mathcal{T}[H|\psi(t)\rangle] \implies -i\hbar\partial_t\mathcal{T}|\psi(t)\rangle = \mathcal{T}H\mathcal{T}^{-1}\mathcal{T}|\psi(t)\rangle

假设 HH 在时间反演下是不变的,即 THT−1=H\mathcal{T}H\mathcal{T}^{-1} = H,则有

−iℏ∂tT∣ψ(t)⟩=HT∣ψ(t)⟩→t→−tiℏ∂tT∣ψ(−t)⟩=HT∣ψ(−t)⟩-i\hbar\partial_t\mathcal{T}|\psi(t)\rangle = H\mathcal{T}|\psi(t)\rangle \xrightarrow{t \to -t} i\hbar\partial_t\mathcal{T}|\psi(-t)\rangle = H\mathcal{T}|\psi(-t)\rangle

上式说明 T∣ψ(−t)⟩\mathcal{T}|\psi(-t)\rangle 也是原来薛定谔方程的一个解,我们称该系统具有时间反演对称性。一般来说,对称性对应守恒量,但是哈密顿量 HH 在反线性时间反演变换下的不变性并不对应任何守恒量。因为哈密顿量的守恒来自于时间平移不变性,且时间反演算符并非厄米的。

取而代之的是,[T,H]=0[\mathcal{T}, H] = 0 说明薛定谔方程的解是成对出现的,即 ∣ψ(t)⟩|\psi(t)\rangle 和 T∣ψ(−t)⟩\mathcal{T}|\psi(-t)\rangle。坐标基底下的薛定谔方程:

iℏ∂tψ(x′,t)=[−ℏ2∇22m+V(x′)]ψ(x′,t)i\hbar\partial_t\psi(\boldsymbol{x}', t) = \left[-\frac{\hbar^2\nabla^2}{2m} + V(\boldsymbol{x}')\right]\psi(\boldsymbol{x}', t)

取复共轭(时间反演)

−iℏ∂tψ∗(x′,t)=[−ℏ2∇22m+V∗(x′)]ψ∗(x′,t)-i\hbar\partial_t\psi^*(\boldsymbol{x}', t) = \left[-\frac{\hbar^2\nabla^2}{2m} + V^*(\boldsymbol{x}')\right]\psi^*(\boldsymbol{x}', t)

取 t→−tt \to -t,那么有

iℏ∂tψ∗(x′,−t)=[−ℏ2∇22m+V∗(x′)]ψ∗(x′,−t)i\hbar\partial_t\psi^*(\boldsymbol{x}', -t) = \left[-\frac{\hbar^2\nabla^2}{2m} + V^*(\boldsymbol{x}')\right]\psi^*(\boldsymbol{x}', -t)

如果 V∗(x′)=V(x′)V^*(\boldsymbol{x}') = V(\boldsymbol{x}'),那么 ψ∗(x′,−t)\psi^*(\boldsymbol{x}', -t) 也是原薛定谔方程的解。注意这里的 −t-t 只是相对于原时间参数 tt 的一个符号变化,表示反向演化,但是 ψ∗(x′,−t)\psi^*(\boldsymbol{x}', -t) 中的每个时刻和对应 ψ∗(x′,t)\psi^*(\boldsymbol{x}', t) 中的时刻遵循同样的物理规律。

以波包为例,时间反演算符的作用就是将波包的动量反向,也就是改变运动方向

在不考虑自旋的情况下,我们可以将坐标表象下的时间反演算符 T^\hat{\mathcal{T}} 等同于复共轭算符 K0K_0:T^=K0\hat{\mathcal{T}} = K_0,即 T^ψ(x,t)=ψ∗(x,t)\hat{\mathcal{T}}\psi(\boldsymbol{x}, t) = \psi^*(\boldsymbol{x}, t)。

这一点也可以通过以下事实得到印证:

  1. 动量算符在坐标表象下为 p^=−iℏ∇\hat{\boldsymbol{p}} = -i\hbar\nabla,复共轭操作会改变其符号,因此 T^p^T^−1=−p^\hat{\mathcal{T}}\hat{\boldsymbol{p}}\hat{\mathcal{T}}^{-1} = -\hat{\boldsymbol{p}}。

  2. 磁场会破坏时间反演对称性,因为 −μL⋅B∝−L^⋅B-\boldsymbol{\mu}_L \cdot \boldsymbol{B} \propto -\hat{\boldsymbol{L}} \cdot \boldsymbol{B} 在时间反演下会变化

我们可以证明 T∣x′⟩=∣x′⟩\mathcal{T}|\boldsymbol{x}'\rangle = |\boldsymbol{x}'\rangle 和 T∣p′⟩=∣−p′⟩\mathcal{T}|\boldsymbol{p}'\rangle = |-\boldsymbol{p}'\rangle。这与我们的经典认识相符。

若 [T,H]=0[\mathcal{T}, H] = 0 且本征态 ∣n⟩|n\rangle 是非简并的,则波函数 ⟨x′∣n⟩\langle x'|n\rangle 必为实函数和一个与 x′x' 无关的相位因子之积

时间反演的平方

连续两次施加时间反演变换应当使物理情形保持不变,即:

T2∣ψ⟩=eiθ∣ψ⟩\mathcal{T}^2|\psi\rangle = e^{i\theta}|\psi\rangle

那么我们如何确定相位因子 eiθe^{i\theta} 的可能取值?

方法 11:首先考察三次时间反演操作

T3∣ψ⟩=TT2∣ψ⟩=Teiθ∣ψ⟩=e−iθT∣ψ⟩\mathcal{T}^3|\psi\rangle = \mathcal{T}\mathcal{T}^2|\psi\rangle = \mathcal{T}e^{i\theta}|\psi\rangle = e^{-i\theta}\mathcal{T}|\psi\rangle

上式利用了 T\mathcal{T} 的反线性性质,将复数常数移出时取了共轭。现在考虑叠加态 ∣ψ⟩+T∣ψ⟩|\psi\rangle + \mathcal{T}|\psi\rangle。对于这个态,施加 T2\mathcal{T}^2 也应产生一个(可能不同的)相位因子 eiθ′e^{i\theta'}:

T2(∣ψ⟩+T∣ψ⟩)=eiθ′(∣ψ⟩+T∣ψ⟩)\mathcal{T}^2(|\psi\rangle + \mathcal{T}|\psi\rangle) = e^{i\theta'}(|\psi\rangle + \mathcal{T}|\psi\rangle)

将两边展开

T2∣ψ⟩+T3∣ψ⟩=eiθ∣ψ⟩+e−iθT∣ψ⟩\mathcal{T}^2|\psi\rangle + \mathcal{T}^3|\psi\rangle = e^{i\theta}|\psi\rangle + e^{-i\theta}\mathcal{T}|\psi\rangle

比较左右两边的系数:

eiθ∣ψ⟩+e−iθT∣ψ⟩=eiθ′∣ψ⟩+eiθ′T∣ψ⟩e^{i\theta}|\psi\rangle + e^{-i\theta}\mathcal{T}|\psi\rangle = e^{i\theta'}|\psi\rangle + e^{i\theta'}\mathcal{T}|\psi\rangle

  ⟹  eiθ=eiθ′=e−iθ  ⟹  e2iθ=1  ⟹  eiθ=±1\implies e^{i\theta} = e^{i\theta'} = e^{-i\theta} \implies e^{2i\theta} = 1 \implies e^{i\theta} = \pm 1

所以有

T2∣ψ⟩=±∣ψ⟩\mathcal{T}^2|\psi\rangle = \pm|\psi\rangle

尽管连续两次时间反演操作看似是一个平凡变换,但其对应的算符 T2\mathcal{T}^2 并不一定是恒等算符。


我们还可以利用作用两次时间反演算符但取具体表象的方法

方法 22:我们在某个特定基底(例如 nn-基底)下考察方程 T2∣ψ⟩=eiθ∣ψ⟩\mathcal{T}^2|\psi\rangle = e^{i\theta}|\psi\rangle

T2∣ψ⟩=U(n)K(n)U(n)K(n)∑n⟨n∣ψ⟩∣n⟩\mathcal{T}^2|\psi\rangle = U_{(n)}K_{(n)}U_{(n)}K_{(n)} \sum_n \langle n|\psi\rangle|n\rangle

其中 ∣n⟩=(0,⋯ ,1,⋯ ,0)T|n\rangle = (0, \cdots, 1, \cdots, 0)^T,且 U(n)U_{(n)} 为幺正矩阵。于是有:

T2∣ψ⟩=U(n)K(n)U(n)∑n⟨ψ∣n⟩∣n⟩=U(n)U(n)∗K(n)∑n⟨ψ∣n⟩∣n⟩=U(n)U(n)∗∑n⟨n∣ψ⟩∣n⟩=U(n)U(n)∗∣ψ⟩\begin{aligned} \mathcal{T}^2|\psi\rangle &= U_{(n)}K_{(n)}U_{(n)} \sum_n \langle \psi|n\rangle|n\rangle = U_{(n)}U_{(n)}^* K_{(n)} \sum_n \langle \psi|n\rangle|n\rangle \\ &= U_{(n)}U_{(n)}^* \sum_n \langle n|\psi\rangle|n\rangle = U_{(n)}U_{(n)}^* |\psi\rangle \end{aligned}

由此可知 U(n)U(n)∗=eiθ1=eiθU(n)U(n)†U_{(n)}U_{(n)}^* = e^{i\theta} \mathbb{1} = e^{i\theta} U_{(n)} U_{(n)}^\dagger。消去左侧的 U(n)U_{(n)} 可得:

U(n)∗=eiθU(n)†→取复共轭U(n)=e−iθU(n)T→取转置U(n)T=e−iθU(n)U_{(n)}^* = e^{i\theta} U_{(n)}^\dagger \xrightarrow{\text{取复共轭}} U_{(n)} = e^{-i\theta} U_{(n)}^T \xrightarrow{\text{取转置}} U_{(n)}^T = e^{-i\theta} U_{(n)}

将 U(n)U_{(n)} 的表达式代回

U(n)T=e−iθ(e−iθU(n)T)=e−2iθU(n)TU_{(n)}^T = e^{-i\theta} (e^{-i\theta} U_{(n)}^T) = e^{-2i\theta} U_{(n)}^T

  ⟹  e−2iθ=1  ⟹  eiθ=±1\implies e^{-2i\theta} = 1 \implies e^{i\theta} = \pm 1

因此 U(n)U_{(n)} 是对称或者反对称矩阵

U(n)=±U(n)TU_{(n)} = \pm U_{(n)}^T

分别对应于 T2∣ψ⟩=±∣ψ⟩]\mathcal{T}^2|\psi\rangle = \pm|\psi\rangle]。

时间反演和自旋

自旋下的时间反演算符

两次时间反演变换并不一定是个恒等算符,我们允许其具有相位因子 −1-1,那么这个因子什么时候会出现呢?就是自旋存在时。

一般来说,对于轨道变量选用坐标表象,而对于自旋算符则选用标准 SzS_z 基底。那么在这两个表象下,时间反演算符 T\mathcal{T} 可以直接等同于复共轭算符 K0K_0 吗?

我们在标准基底下考察泡利矩阵的性质(SyS_y 为纯虚矩阵,Sx,SzS_x, S_z 为实矩阵)

K0SyK0=Sy∗K0K0=−Sy,K0SxK0=Sx∗K0K0=Sx,K0SzK0=Sz∗K0K0=SzK_0 S_y K_0 = S_y^* K_0 K_0 = -S_y, \quad K_0 S_x K_0 = S_x^* K_0 K_0 = S_x, \quad K_0 S_z K_0 = S_z^* K_0 K_0 = S_z

显然,仅用 K0K_0 无法满足时间反演对自旋的要求 TST−1\mathcal{T}\boldsymbol{S}\mathcal{T}^{-1}(即所有分量都应反号)。由于 T=UK0\mathcal{T} = U K_0,为了得到正确的变换性质 TST−1=−S\mathcal{T}\boldsymbol{S}\mathcal{T}^{-1} = -\boldsymbol{S},我们要凑出幺正算符 UU,其必须具有以下作用:

−Sy=UK0SyK0U†=USy∗K0K0U†=−USyU†−Sx=UK0SxK0U†=USx∗K0K0U†=USxU†−Sz=UK0SzK0U†=USz∗K0K0U†=USzU†\begin{aligned} -S_y &= U K_0 S_y K_0 U^\dagger &= U S_y^* K_0 K_0 U^\dagger &= -U S_y U^\dagger \\[0.5em] -S_x &= U K_0 S_x K_0 U^\dagger &= U S_x^* K_0 K_0 U^\dagger &= U S_x U^\dagger \\[0.5em] -S_z &= U K_0 S_z K_0 U^\dagger &= U S_z^* K_0 K_0 U^\dagger &= U S_z U^\dagger \end{aligned}

观察上述方程可知,UU 算符的作用是翻转 SxS_x 和 SzS_z 的符号,同时保持 SyS_y 不变。这就对应于自旋空间中绕 yy 轴旋转 π\pi 角度的转动

U=e−iℏSyπU = e^{-\frac{i}{\hbar}S_y\pi}

自旋下的时间反演算符

因此,包含自旋的时间反演算符专门选定 yy 方向的转动

T=e−iℏSyπK0\mathcal{T} = e^{-\frac{i}{\hbar}S_y\pi} K_0

值得注意的是,在此基底下矩阵 U=e−iℏSyπU = e^{-\frac{i}{\hbar}S_y\pi} 是实正交矩阵,即满足 U†=UTU^\dagger = U^T


对于 S=1/2S = 1/2 的情形,我们有:

T(αβ)=e−i2σyπK0(αβ)=(cos⁡π2−iσysin⁡π2)(α∗β∗)=−iσy(α∗β∗)=(0−110)(α∗β∗)=(−β∗α∗)\begin{aligned} \mathcal{T} \begin{pmatrix} \alpha \\[0.5em] \beta \end{pmatrix} &= e^{-\frac{i}{2}\sigma_y\pi} K_0 \begin{pmatrix} \alpha \\[0.5em] \beta \end{pmatrix} = \left(\cos\frac{\pi}{2} - i\sigma_y \sin\frac{\pi}{2}\right) \begin{pmatrix} \alpha^* \\[0.5em] \beta^* \end{pmatrix} \\[0.5em] &= -i\sigma_y \begin{pmatrix} \alpha^* \\[0.5em] \beta^* \end{pmatrix} = \begin{pmatrix} 0 & -1 \\[0.5em] 1 & 0 \end{pmatrix} \begin{pmatrix} \alpha^* \\[0.5em] \beta^* \end{pmatrix} = \begin{pmatrix} -\beta^* \\[0.5em] \alpha^* \end{pmatrix} \end{aligned}

于是时间反演的平方 T2\mathcal{T}^2 为

T2=e−iℏSyπK0e−iℏSyπK0=e−iℏSyπe−iℏSyπK0K0=e−iℏSy2π\mathcal{T}^2 = e^{-\frac{i}{\hbar}S_y\pi} K_0 e^{-\frac{i}{\hbar}S_y\pi} K_0 = e^{-\frac{i}{\hbar}S_y\pi} e^{-\frac{i}{\hbar}S_y\pi} K_0 K_0 = e^{-\frac{i}{\hbar}S_y 2\pi}

这里利用了 e−iℏSyπe^{-\frac{i}{\hbar}S_y\pi} 是实矩阵,且考虑到轨道部分的贡献 e−iℏLy2π=1e^{-\frac{i}{\hbar}L_y 2\pi} = 1,最终得到:

T2=e−iℏJy2π\mathcal{T}^2 = e^{-\frac{i}{\hbar}J_y 2\pi}

T2\mathcal{T}^2 等价于绕 yy 轴旋转一整周的旋转算符

T2={1对于整数j 的任意态−1对于半奇整数j 的任意态\mathcal{T}^2 = \begin{cases} 1 & \text{对于整数} j \text{ 的任意态} \\[1em] -1 & \text{对于半奇整数} j \text{ 的任意态} \end{cases}

由此可推导幺正算符 UU 的性质:

对于 S=S = 整数:e−iℏSy2π=1  ⟹  e−iℏSyπ=eiℏSyπ  ⟹  U=U†=UTe^{-\frac{i}{\hbar}S_y 2\pi} = 1 \implies e^{-\frac{i}{\hbar}S_y\pi} = e^{\frac{i}{\hbar}S_y\pi} \implies U = U^\dagger = U^T;
对于 S=S = 半奇整数:e−iℏSy2π=−1  ⟹  e−iℏSyπ=−eiℏSyπ  ⟹  U=−U†=−UTe^{-\frac{i}{\hbar}S_y 2\pi} = -1 \implies e^{-\frac{i}{\hbar}S_y\pi} = -e^{\frac{i}{\hbar}S_y\pi} \implies U = -U^\dagger = -U^T。

总结以上两点,可以得到通式:

e−iℏSyπ=(−1)2SeiℏSyπe^{-\frac{i}{\hbar}S_y\pi} = (-1)^{2S} e^{\frac{i}{\hbar}S_y\pi}

自旋的时间反演实例

小 d 矩阵就是绕 yy 轴转动的生成元,我们来考察 Wigner dd-矩阵在 β=π\beta = \pi(即绕 yy 轴旋转 180∘180^\circ)时的具体数值。对于自旋 j=1/2j=1/2:

d(1/2)(β)=(cos⁡β2−sin⁡β2sin⁡β2cos⁡β2)  ⟹  d(1/2)(π)=(0−110)d^{(1/2)}(\beta) = \begin{pmatrix} \cos\dfrac{\beta}{2} & -\sin\dfrac{\beta}{2} \\[0.5em] \sin\dfrac{\beta}{2} & \cos\dfrac{\beta}{2} \end{pmatrix} \implies d^{(1/2)}(\pi) = \begin{pmatrix} 0 & -1 \\[0.5em] 1 & 0 \end{pmatrix}

这正是我们之前推导出的幺正算符 UU 的形式。对于自旋 j=1j=1:

d(1)(β)=(12(1+cos⁡β)−sin⁡β212(1−cos⁡β)sin⁡β2cos⁡β−sin⁡β212(1−cos⁡β)sin⁡β212(1+cos⁡β))  ⟹  d(1)(π)=(0010−10100)d^{(1)}(\beta) = \begin{pmatrix} \dfrac{1}{2}(1+\cos\beta) & -\dfrac{\sin\beta}{\sqrt{2}} & \dfrac{1}{2}(1-\cos\beta) \\[1em] \dfrac{\sin\beta}{\sqrt{2}} & \cos\beta & -\dfrac{\sin\beta}{\sqrt{2}} \\[1em] \dfrac{1}{2}(1-\cos\beta) & \dfrac{\sin\beta}{\sqrt{2}} & \dfrac{1}{2}(1+\cos\beta) \end{pmatrix} \implies d^{(1)}(\pi) = \begin{pmatrix} 0 & 0 & 1 \\[1em] 0 & -1 & 0 \\[1em] 1 & 0 & 0 \end{pmatrix}

由此可见,一般规律为 dm′m(j)(π)∝δm′,−md_{m'm}^{(j)}(\pi) \propto \delta_{m',-m},即该旋转操作将磁量子数 mm 翻转为 −m-m。


NN 粒子系统的时间反演

考虑一个由 NN 个粒子组成的系统,若粒子均为零自旋,则对任意粒子数 NN,均有 T2=1\mathcal{T}^2 = 1。若粒子均为自旋 1/21/2,系统的时间反演算符是各个粒子时间反演算符的直积:

T=e−iℏS1yπK0⊗e−iℏS2yπK0⊗⋯⊗e−iℏSNyπK0\mathcal{T} = e^{-\frac{i}{\hbar}S_1^y\pi}K_0 \otimes e^{-\frac{i}{\hbar}S_2^y\pi}K_0 \otimes \cdots \otimes e^{-\frac{i}{\hbar}S_N^y\pi}K_0

由于每个费米子贡献一个 −1-1 因子(即 Ti2=−1\mathcal{T}_i^2 = -1),总的效果为

T2=(−1)⊗(−1)⊗⋯⊗(−1)=(−1)N\mathcal{T}^2 = (-1) \otimes (-1) \otimes \cdots \otimes (-1) = (-1)^N

这表明,当 NN 为偶数时,T2=1\mathcal{T}^2 = 1(表现类似玻色子);当 NN 为奇数时,T2=−1\mathcal{T}^2 = -1(表现类似费米子)。

克莱默定理:如果一个由奇数个自旋 1/21/2 粒子组成的系统,其哈密顿量在时间反演下保持不变,那么它的所有定态都是简并的。

证明:设 ∣E,N⟩|E, N\rangle 为哈密顿量 HH 的本征态,本征值为 EE

H∣E,N⟩=E∣E,N⟩H|E, N\rangle = E|E, N\rangle

由于系统具有时间反演不变性 TH=HT\mathcal{T}H = H\mathcal{T},我们有

HT∣E,N⟩=TH∣E,N⟩=ET∣E,N⟩H\mathcal{T}|E, N\rangle = \mathcal{T}H|E, N\rangle = E\mathcal{T}|E, N\rangle

这表明 T∣E,N⟩\mathcal{T}|E, N\rangle 也是 HH 的本征态,且能量同样为 EE。假设 ∣E,N⟩|E, N\rangle 是非简并的,则 T∣E,N⟩\mathcal{T}|E, N\rangle 必定与 ∣E,N⟩|E, N\rangle 共线(仅相差一个相位因子):

T∣E,N⟩=eiθ∣E,N⟩\mathcal{T}|E, N\rangle = e^{i\theta}|E, N\rangle

对该式再次作用 T\mathcal{T}:

T2∣E,N⟩=T(eiθ∣E,N⟩)=e−iθT∣E,N⟩=e−iθeiθ∣E,N⟩=∣E,N⟩\mathcal{T}^2|E, N\rangle = \mathcal{T}(e^{i\theta}|E, N\rangle) = e^{-i\theta}\mathcal{T}|E, N\rangle = e^{-i\theta}e^{i\theta}|E, N\rangle = |E, N\rangle

然而,对于奇数个自旋 1/21/2 粒子(NN 为奇数),我们已知 T2=−1\mathcal{T}^2 = -1,即应有 T2∣E,N⟩=−∣E,N⟩\mathcal{T}^2|E, N\rangle = -|E, N\rangle。这与上述推导矛盾,因此原假设错误,∣E,N⟩|E, N\rangle 必定是简并的。

一般的简并都是由于系统的对称性引起的,而克莱默定理表明,时间反演对称也能造成简并。对于处于非对称电场中且存在自旋-轨道耦合 (L⋅S)(\boldsymbol{L} \cdot \boldsymbol{S}) 的系统,该定理并非平凡结论。在这种奇数个粒子的系统中,自旋和轨道角动量均不守恒。克莱默定理指出,因为电场附加能也具有时间反演对称性,所以纯电场无法消除克莱默简并。但是外部磁场可以打破这种简并,因为磁场与磁矩耦合,在哈密顿量中引入了形如 γS⋅B\gamma \boldsymbol{S} \cdot \boldsymbol{B} 的项,而该项不具有时间反演不变性。