晶格平移

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

一维周期势

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

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

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

T(na)V(q)T(na)=eip(na)V(q)eip(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

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

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

格点处的态

那么有

T(na)m=m+nT(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=ε0j=jjJj=(jj+1+j+1j),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} 态上,有

Hn=ε0jjjnJj(jj+1n+j+1jn)=ε0jjδjnJj(jδj+1,n+j+1δjn)=ε0nJ(n1+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 的本征态。但我们要的是 HHT(na)T(na) 的共同本征态,怎么办?线性展开!

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

k=n=Cn(k)n|k\rangle=\sum_{n=-\infty}^\infty C_n(k)|n\rangle

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

Hk=EkkT(a)k=eiθk(a)kH|k\rangle=E_k|k\rangle \qquad T(a)|k\rangle=e^{i\theta_k(a)}|k\rangle

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

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

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

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

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

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

Ek=ε0J(eika+eika)=ε02JcoskaE_k=\varepsilon_0-J\left(e^{-ika}+e^{ika}\right)=\varepsilon_0-2J\cos ka

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

  • 这就相当于从局域态的表象经幺正变换到能量本征态的表象,又因为希尔伯特空间的维数不因表象变换而变化,所以 kk 也是连续变量,理论上可以取任意实数值,但是由于 kk 出现在 coska\cos ka 中,本征能量连续分布在 ε02J\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=ε0j=jj,H0j=ε0jH_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 的波函数为:qk=Aneiknaqn\langle q^{\prime}|k\rangle=A\sum_ne^{-ikna}\langle q^{\prime}|n\rangle

替换为 q+aq'+a 仍然成立:q+ak=Aneiknaq+an\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 处的波函数和 nm\ket{n-m}qq' 处的波函数相同,从数学上来讲就是

q+man=qT(ma)n=qT(ma)n=qnm\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+ak=Aneiknaqn1=Aneik(n+1)aqn=eikaAneiknaqn=eikaqk\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+ak=eik(q+a)eikaqk=eikqqke^{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)eikqqku_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)

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

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

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

环上的有限尺寸问题

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

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

H=ε0j=1NjjJj=1N(jj+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 为偶整数。指标 jj11 取到 NN,因为假设世界是一个长度为 NN 的环,使得第 (N+m)(N + m) 个点就是第 mm 个点,也就是说满足周期性边界条件

TN(a)m=N+m=mT^N(a)|m\rangle=|N+m\rangle=|m\rangle

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

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

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

k1NΣn=1Neiknan|k\rangle\equiv\frac{1}{\sqrt{N}}\Sigma_{n=1}^Ne^{-ikna}|n\rangle

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

T(a)k=T(a)1Nn=1Neiknan=1Nn=1Neiknan+1=1Nn=2N+1eik(n1)an=1Nn=2Neik(n1)an+1NeikNaN+1=eika1Nn=2Neiknan+1NeikNa1\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=eikak=eika1Nn=2Neiknan+eika1Neika1T(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

对比上两式得到

eikNa=1k=2πmNa,mZe^{-ikNa}=1\quad\Longrightarrow\quad k=\frac{2\pi m}{Na},m\in\mathbb{Z}

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

kKN=even={πa,πa+2πNa,,2πNa,0,2πNa,,πa2π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β=1Nn=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α(nm)a=δnm\frac{1}{N} \sum_{\alpha=1}^{N} e^{ik_\alpha(n-m)a} = \delta_{nm}

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

n1NΣα=1Neikαnakα|n\rangle\equiv\frac{1}{\sqrt{N}}\Sigma_{\alpha=1}^Ne^{ik_\alpha\mathrm{n}a}|k_\alpha\rangle

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

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

j=1Njj+n=j=1N1Nα=1Neikαjakα1Nβ=1Neikβ(j+n)akβ=α,βeikβna[1Nj=1Nei(kαkβ)ja]kαkβ=α,βeikβnaδαβkαkβ=αeikαnakα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}

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

H=ε0j=1NjjJj=1N(jj+1+H.c.)=ε0αeikα0akαkαJ(αeikα1akαkα+H.c.)=α(ε02Jcoskα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})。
赝标量:在旋转下不变但在空间反演下变号的量,由轴矢量和极矢量的乘积组合而成(例如 EB\boldsymbol{E} \cdot \boldsymbol{B}SP\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+11-1

πx\pi|\boldsymbol{x}'\ranglex\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}π\pix\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\rangleHH 的一个非简并本征态,其本征值为 EnE_nHn=EnnH|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 衰变:

2760Co2860Ni+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(在实验室中),观察到 ψSPψ<0\langle\psi|\boldsymbol{S} \cdot \boldsymbol{P}|\psi\rangle < 0。然而,对于反演态 ψ=πψ|\psi'\rangle = \pi|\psi\rangle,我们有

ψSPψ=ψπSPπψ=ψSPψ>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 要么是简并态,要么是同一个态(非简并)。但是它们又不可能是同一个态,否则 ψSPψ=ψSPψ=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 = -ax=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(aa)|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=(E0VVE0)=E0Vσx    {HS=(E0V)SHA=(E0+V)AH = \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)=Ha=12(ei(E0V)tS+ei(E0+V)tA)=12ei(E0V)t(S+ei2VtA)|\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/\hbara|a\ranglea|-a\rangle 之间来回振荡。若是势垒有限高但是又足够大,那么 V0T=π2VV \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ψ=eiHt(TeiHtψ)\mathcal{T}|\psi\rangle = e^{-\frac{i}{\hbar}Ht}(\mathcal{T}e^{-\frac{i}{\hbar}Ht}|\psi\rangle)

    eiHtTψ=TeiHtψ\implies e^{\frac{i}{\hbar}Ht}\mathcal{T}|\psi\rangle = \mathcal{T}e^{-\frac{i}{\hbar}Ht}|\psi\rangle

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

(1+iHdt)Tψ=T(1iHdt)ψ\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=THHn=Enn\begin{cases} H\mathcal{T} = -\mathcal{T}H \\ H|n\rangle = E_n|n\rangle \end{cases}

    HTn=THn=(En)Tn\implies H\mathcal{T}|n\rangle = -\mathcal{T}H|n\rangle = (-E_n)\mathcal{T}|n\rangle

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

反线性算符

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

A(c1ψ1+c2ψ2)=c1Aψ1+c2Aψ2A(c_1|\psi_1\rangle + c_2|\psi_2\rangle) = c_1^* A|\psi_1\rangle + c_2^* A|\psi_2\rangle

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

A2A1cψ=A2cA1ψ=cA2A1ψA_2 A_1 c|\psi\rangle = A_2 c^* A_1 |\psi\rangle = c A_2 A_1 |\psi\rangle

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

与线性算符不同,反线性操作依赖于基底的选择。例如有 A^[ciψ]=ciA^ψ\hat{A}[ci|\psi\rangle] = -c^*i \hat{A}|\psi\rangle,但是如果我们令 ψ=iψ|\psi'\rangle = i|\psi\rangle,那么 A^[cψ]=cA^ψ=ciA^ψ\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)nnψn=nψnK(n)n=nψnnK_{(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=nK_{(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(ν)ψ=νψνν=νψνnnνn=n(νψνnν)nK(n)ψ=nψnn=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(1i)=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. 存在逆算符 A1A^{-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ϕ=ϕALA1ψ\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ϕ=ϕAFA1ψ    ψFψ=ψAFA1ψ\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ψ=ψAFA1ψ\langle\psi|F|\psi\rangle = \langle\psi'|AFA^{-1}|\psi'\rangle,我们得到:

ψpψ=ψpψ=ψTpT1ψ-\langle\psi'|\boldsymbol{p}|\psi'\rangle = \langle\psi|\boldsymbol{p}|\psi\rangle = \langle\psi'|\mathcal{T}\boldsymbol{p}\mathcal{T}^{-1}|\psi'\rangle

对比得到

TpT1=p\mathcal{T}\boldsymbol{p}\mathcal{T}^{-1} = -\boldsymbol{p}

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

TxT1=x    TLT1=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},我们同样希望他们与轨道角动量具有相同的变换性质:

TST1=S,TJT1=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 在时间反演下的变换:

TiT1=T[xj,pj]T1=TxjT1TpjT1TpjT1TxjT1=xj(pj)(pj)xj=[xj,pj]=i    TiT1=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[itψ(t)]=T[Hψ(t)]    itTψ(t)=THT1Tψ(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 在时间反演下是不变的,即 THT1=H\mathcal{T}H\mathcal{T}^{-1} = H,则有

itTψ(t)=HTψ(t)ttitTψ(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)\rangleTψ(t)\mathcal{T}|\psi(-t)\rangle。坐标基底下的薛定谔方程:

itψ(x,t)=[222m+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)

取复共轭(时间反演)

itψ(x,t)=[222m+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)

ttt \to -t,那么有

itψ(x,t)=[222m+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_0T^=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. 磁场会破坏时间反演对称性,因为 μLBL^B-\boldsymbol{\mu}_L \cdot \boldsymbol{B} \propto -\hat{\boldsymbol{L}} \cdot \boldsymbol{B} 在时间反演下会变化

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

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

时间反演的平方

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

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

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

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

T3ψ=TT2ψ=Teiθψ=eiθ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θψ+eiθTψ\mathcal{T}^2|\psi\rangle + \mathcal{T}^3|\psi\rangle = e^{i\theta}|\psi\rangle + e^{-i\theta}\mathcal{T}|\psi\rangle

比较左右两边的系数:

eiθψ+eiθ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θ=eiθ    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)nnψ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ψnn=U(n)U(n)K(n)nψnn=U(n)U(n)nnψ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)=eiθU(n)T取转置U(n)T=eiθ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=eiθ(eiθU(n)T)=e2iθU(n)TU_{(n)}^T = e^{-i\theta} (e^{-i\theta} U_{(n)}^T) = e^{-2i\theta} U_{(n)}^T

    e2iθ=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=SyK0K0=Sy,K0SxK0=SxK0K0=Sx,K0SzK0=SzK0K0=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 无法满足时间反演对自旋的要求 TST1\mathcal{T}\boldsymbol{S}\mathcal{T}^{-1}(即所有分量都应反号)。由于 T=UK0\mathcal{T} = U K_0,为了得到正确的变换性质 TST1=S\mathcal{T}\boldsymbol{S}\mathcal{T}^{-1} = -\boldsymbol{S},我们要凑出幺正算符 UU,其必须具有以下作用:

Sy=UK0SyK0U=USyK0K0U=USyUSx=UK0SxK0U=USxK0K0U=USxUSz=UK0SzK0U=USzK0K0U=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_xSzS_z 的符号,同时保持 SyS_y 不变。这就对应于自旋空间中绕 yy 轴旋转 π\pi 角度的转动

U=eiSyπU = e^{-\frac{i}{\hbar}S_y\pi}

自旋下的时间反演算符

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

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

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


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

T(αβ)=ei2σyπK0(αβ)=(cosπ2iσysinπ2)(αβ)=iσy(αβ)=(0110)(αβ)=(βα)\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=eiSyπK0eiSyπK0=eiSyπeiSyπK0K0=eiSy2π\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}

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

T2=eiJy2π\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 = 整数:eiSy2π=1    eiSyπ=eiSyπ    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 = 半奇整数:eiSy2π=1    eiSyπ=eiSyπ    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

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

eiSyπ=(1)2SeiSyπe^{-\frac{i}{\hbar}S_y\pi} = (-1)^{2S} e^{\frac{i}{\hbar}S_y\pi}

自旋的时间反演实例

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

d(1/2)(β)=(cosβ2sinβ2sinβ2cosβ2)    d(1/2)(π)=(0110)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(1cosβ)sinβ2cosβsinβ212(1cosβ)sinβ212(1+cosβ))    d(1)(π)=(001010100)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}

由此可见,一般规律为 dmm(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=eiS1yπK0eiS2yπK0eiSNyπ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

HE,N=EE,NH|E, N\rangle = E|E, N\rangle

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

HTE,N=THE,N=ETE,NH\mathcal{T}|E, N\rangle = \mathcal{T}H|E, N\rangle = E\mathcal{T}|E, N\rangle

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

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

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

T2E,N=T(eiθE,N)=eiθTE,N=eiθ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,即应有 T2E,N=E,N\mathcal{T}^2|E, N\rangle = -|E, N\rangle。这与上述推导矛盾,因此原假设错误,E,N|E, N\rangle 必定是简并的。

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