外观
Coleman QFT Lecture 20: Dirac方程 II:求解
约 5414 字大约 18 分钟
总结
先对上一章做一个简单的回顾:
- Lorentz群的表示 D(21,0),D(0,21) 是旋量表示,其表示空间中生活着Weyl旋量 u±,两者在旋转下的变换相同,boost下的变换不同,可以通过宇称相互转换。通过Pauli矩阵可以构造作为矢量分量的双线性型:
Vμ=(u+†u+,u+†σiu+),Wμ=(W0,Wi)=(u−†u−,−u−†σiu−)
- 矢量的散度是标量,可以用于构造Lagrangian:
L+=±i(u+†∂0u++u+†σ⋅∇u+)L−=±i(u−†∂0u−−u−†σ⋅∇u−)
其变分可以得到Weyl方程。两种旋量场分别代表了右手与左手(螺旋度 ±1/2 )的无质量粒子。Weyl理论不具有宇称对称性,两种粒子在宇称下相互转化。
- 将两种旋量场都放进Lagrangian中,可以构造出含质量项的宇称不变的Lagrangian:
L=±[i(u+†∂0u++u+†σ⋅∇u++u−†∂0u−−u−†σ⋅∇u−)−m(u+†u++u−†u+)]
通过定义Weyl基底下的Dirac旋量与Dirac矩阵:
ψ≡(u+u−),α=(σ00−σ),β=(0110)
可以得到更简单的Lagrangian:
L=±[iψ†(∂0+α⋅∇)ψ−mψ†βψ]
以及变分得到的Dirac方程:
i∂0ψ+iα⋅∇ψ=mβψ
20.1 Dirac 基底
Dirac矩阵是厄米的,其满足称为Clifford代数的反对易关系:
{αi,αj}=2δij,{αi,β}=0,{β,β}=2
这可以很容易由Pauli矩阵的反对易性质推出:
{σi,σj}=2δij
Dirac旋量是两个Weyl旋量直接放在一起得到的旋量,其对应的Lorentz群表示是两个不可约旋量表示的直和 D(21,0)⊕D(0,21)。我们知道对两个旋量表示,生成元的表示为:
L=21σ,M=±2iσ
因此其直和为:
L=21(σ00σ)≡21Σ,M=2i(σ00−σ)=2iα
宇称对Dirac旋量的作用是交换 u±,即
P:ψ(x,t)→βψ(−x,t)
以上,我们都是在Weyl基底下进行讨论,其称为Dirac方程的Weyl表示。但当时Dirac用的并不是Weyl基底,他选择的是:
ψ=21(u++u−u+−u−)
这种基底称为Dirac表示或者标准表示(standard representation)。可以发现 iψ†∂0ψ=iu+†∂0u++iu−†∂0u− 不变。而Dirac矩阵变为:
α=(0σσ0),β=(100−1)
此时的生成元表示仍满足:
L=21Σ,M=2iα
提示
Dirac表示与Weyl表示只是基底变换而已,其通过下面的相似变换联系:
ψD=TψW,αD=TαWT−1,βD=TβWT−1,T=21(111−1)
20.2 平面波解
接下来我们来实际求解Dirac方程。不论在何种基底下,Dirac方程形式上总可写为:
i∂0ψ+iα⋅∇ψ=βmψ
由于其描述的是自由的质量为 m 的粒子,我们自然会猜测平面波解。首先,我们猜测指数上带负号的正频率解:
\psi(\vec x,t)=u_\vec p\,\mathrm{e}^{-ip\cdot x}
注意不要混淆其与Coleman QFT - Lecture 19. Dirac方程 I:构造拉格朗日量中求解Weyl方程时的 u_\vec p,现在的 u_\vec p 是四分量的列向量。后面我们会说明,Dirac方程可以推出质量为 m 的Klein-Gordan方程,因此平面波的波矢应当处在质壳上,即
Ep=p2+m2
则Dirac方程变为:
[E_p-\vec \alpha\cdot\vec p\,]u_\vec p=\beta mu_\vec p
首先,我们考虑静止粒子( p=0 ),则
u0=βu0
因此其即 β 的本征矢。具体考虑标准表示,此时 β=diag(1,−1)。因此本征矢为:
u0(1)=2m1000,u0(2)=2m0100
这里前面的系数 2m 的作用在后面会显现出来。那么,动量不为零的解又是什么呢?我们不需要直接求解Dirac方程,而是对上面的解施加boost即可。
u_\vec p^{(r)}=\mathrm{e}^{-i\hat a\cdot\vec M\phi}u_0^{(r)}=\mathrm{e}^{\frac12\vec\alpha\cdot\hat a\phi}u_0^{(r)}
其中 $$\hat a=\dfrac{\vec p}{|\vec p|},\quad \cosh\phi=\dfrac{E_p}{m}$$为了简单起见,考虑沿 z 轴的动量,由于
e21αzϕ=cosh21ϕ+αzsinh21ϕ=cosh2ϕ0sinh2ϕ00cosh2ϕ0−sinh2ϕsinh2ϕ0cosh2ϕ00−sinh2ϕ0cosh2ϕ
其中
cosh2ϕ=2coshϕ+1=2mEp+m,sinh2ϕ=2mEp−m
因此
up(1)=Ep+m0Ep−m0,up(2)=0Ep+m0−Ep−m,
我们可以做一个简单的验证。前面知道
Vμ=(u+†u+,u+†σiu+),Wμ=(W0,Wi)=(u−†u−,−u−†σiu−)
构成两个矢量。我们也可以将其合成一个由Dirac旋量构造的矢量,即
Uμ=(ψ†ψ,ψ†αψ)
对于 p=0,有下面的关系:
(u0(r)†u0(s),u0(r)†αu0(s))=(2m,0)δrs
而上面的解满足
(up(r)†up(s),up(r)†αup(s))=2(Ep,p)δrs
这正是能动量的Lorentz变换。 接下来,我们来看另一部分,关于负频率的平面波解,即
\psi(\vec x,t)=v_\vec p\mathrm{e}^{+ip\cdot x},\quad p^0=E_p=\sqrt{\vec p^2+m^2}
代入Dirac方程得到:
[E_p-\vec\alpha\cdot\vec p\,]v_\vec p=-\beta mv_\vec p
应用相同的流程。对于零动量 p=0 情形:
v0=−βv0⇒v0(1)=2m0001,v0(2)=2m0010
同样用boost得到沿 z 方向运动的解:
vp(1)=0−Ep−m0Ep+m,vp(2)=Ep−m0Ep+m0
以及正交归一关系:
(vp(r)†vp(s),vp(r)†αvp(s))=2(Ep,p)δrs
负频率似乎代表着负能解,这一点当时引起了很多人的困扰。毕竟如果有负能解的存在,意味着能量无下界,会导致物理系统的不稳定。Dirac的提议是: 之所以带正能的电子不会掉落到负能态,是因为“除了会有少量速度较小的态, 所有的负能态已经被占满了。”随后,这种空穴解释被正电子的图像所取代:负能解代表带正电的正电子。当然,这样解释可以避免能量不正定的原因需要依靠量子化Dirac理论完成,这将在下一章进行。
最后,我们来说明平面波解也是角动量的本征态。由于我们考虑的是 z 方向的boost,因此可以仅考虑 z 方向的角动量 Lz,其矩阵表示为:
Lz=21(σz00σz)=2110000−1000010000−1
对于零动量解,易验证:
Lzu0(1)=21u0(1),Lzv0(2)=21v0(2),Lzu0(2)=−21u0(2),Lzv0(1)=−21v0(1)
而对于 pz=0 的粒子,由于沿 z 方向的boost和 Lz 对易,故上面的本征关系不变。因此,角标 1,2 标记了角动量本征态,或者说螺旋度本征态。
20.3 Pauli 定理
接下来,我们稍微偏题一下,进行一些数学上的讨论。现在,我已经介绍了Dirac表示与Weyl表示,两者之间可以通过相似变换相联系。但Dirac矩阵的许多性质是不依赖于具体基底的,而是源于其代数结构。而这里说的代数结构,就是指前面提到的Clifford algebra:
{αi,αj}=2δij,{αi,β}=0,{β,β}=2
这一关系显然在相似变换下不变。那么,Clifford代数能否将Dirac矩阵限制在至多差一个相似变换的范围内呢?答案是可以的。这就是Pauli定理,表述如下:
Pauli theorem
任何一组满足Dirac代数的 4×4 矩阵表示都等价于Weyl表示。
此定理暗示了下面的结果:所有只包含自旋1/2粒子且宇称对称的Lorentz群的表示都是等价的。也就是说,Dirac理论是唯一可以构建出的有宇称对称性的自旋1/2粒子的理论。下面就来证明该定理: 首先,定义
Mi=2iαi,[Mi,Mj]=−iϵijkLk
我们来证明这样定义出的 L 与 M 满足Lorentz群生成元的对易关系。会用到的引理:
ϵijkϵijl=2δkl,[AB,C]=A{B,C}−{A,C}B
则
−iϵijkLk=−41[αi,αj]=−21αiαj⇒Lk=−4iϵijkαiαj
[Lk,Ml]=81ϵijk[αiαj,αl]=41ϵijk(αiδjl−δilαj)=41(ϵilkαi−ϵljkαj)=−21ϵkliαi=iϵkliMi
[Li,Lj]=−161ϵiklϵjmn[αkαl,αmαn]=−161ϵiklϵjmn([αkαl,αm]αn+αm[αkαl,αn])=−81ϵiklϵjmn(αkαnδlm−αlαnδkm+αmαkδln−αmαlδkn)=−81(αjαi−δij−δij+αjαi+δij−αiαj−αiαj+δij)=41(αiαj−αjαi)=iϵijkLk
这样就验证了生成元的对易关系,因此构造出的 L,M 构成了Lorentz群的一个4维表示。那么是哪一个呢? 首先,不可约表示只可能是 D(23,0),D(21,21),D(0,23),但其不可能,因为
Lz2=−41αxαyαxαy=41αx2αy2=41
因此 L 的表示矩阵只有本征值 ±1/2。而现在,这三个不可约表示都存在超过 l=21 的角动量部分,因此是不可能的。换句话说,我们找的表示应该只含有自旋1/2表示,因此只能是旋量表示的直和,有三种可能:
D(21,0)⊕D(21,0),D(21,0)⊕D(0,21),D(0,21)⊕D(0,21)
接下来使用 β,由恒等式:
βαβ−1=βαβ=−α
因此
βMβ−1=−M,βLβ−1=L
我们注意到,此变换的作用与宇称变换完全相同。因此 β 就是宇称变换下的相似变换矩阵,也就是说此表示有宇称不变性。故我们只剩下一个选择,即表示为 D(21,0)⊕D(0,21)。 以上所有的过程都只用到了Clifford代数的条件,因此对于任意的满足Clifford代数的 4×4 矩阵 {α,β},其都会生成相同的表示 D(21,0)⊕D(0,21)。由于 M 正比于 α,故存在相似变换矩阵 T,使得:
TαT−1=αW=(σ00−σ)
由于 β 与 α 反对易且模方为1,其形式一定为:
TβT−1=β′=(0λ−1⋅1λ⋅10)
现在再做一个相似变换:
S=(λ⋅1001):αW→SαWS−1=αW,β′→Sβ′S−1=(0110)=βW
这样就证明了任意一组满足Clifford代数的 4×4 矩阵都与Weyl表示等价。
20.4 γ 矩阵
下面我将介绍关于Dirac矩阵的快速运算方法。首先假设Dirac矩阵式是厄米的:
α†=α,β†=β
Weyl表示与Dirac表示都是满足条件的。但这个性质的确依赖于基底。现在定义旋量的“伴随”:
ψˉ=ψ†β
ψˉ 称为Dirac伴随。矩阵(线性变换)的Dirac伴随定义为:
Aˉ≡βA†β
由于 β 的逆即为自身,Dirac伴随与厄米共轭有许多相同的性质。比如
D(Λ)ψ=ψD(Λ)
由于 ψ†βψ=ψˉψ 为Lorentz标量,因此
Λ:ψˉψ→ψD(Λ)D(Λ)ψ=ψˉψ
因此 D(Λ)D(Λ)=1,也就是说Lorentz群表示虽然不是厄米的,但在Dirac伴随下是等距的。这称为“Dirac幺正性”。
我们前面提到,用Dirac旋量构造的矢量形式为:
Uμ=(ψ†ψ,ψ†αψ)
我们想把所有的 ψ† 改造成 ψˉ 的形式,因此在其中插入 β2=I:
Uμ=(ψ†β2ψ,ψ†β2αψ)=(ψˉβψ,ψˉβαψ)=ψˉγμψ
其中定义了gamma矩阵
\gamma^\mu=(\gamma^0,\vec\gamma)\equiv(\beta,\beta\vec\alpha)$$这样,原先的厄米共轭全部用Dirac伴随替代,Dirac矩阵也用新的gamma矩阵替代。这样,构造矢量就可以直接看角标了:$\gamma^\mu$ 自带一个矢量指标,因此只用取其关于 $\psi,\bar\psi$ 的二次型即可得到矢量 $U^\mu=\bar\psi\gamma^\mu\psi$。 在标准表示下,gamma矩阵的显式矩阵为:
\gamma^0=\beta=\begin{pmatrix} 0&\mathbb 1\ \mathbb 1& 0 \end{pmatrix},\quad \gammai=\beta\alphai=\begin{pmatrix} 0&\sigma^i\ -\sigma^i& 0 \end
gamma矩阵是Lorentz变换下的矢量,因此:
\overline {D(\Lambda)},\gamma^\mu D(\Lambda)=\Lambda\mu_\nu\gamma\nu
Clifford代数现在可以用gamma矩阵简化为:
{\gamma\mu,\gamma\nu}=2,\eta^{\mu\nu}\mathbb{1}_{4\times 4}
这里 $\eta^{\mu\nu}=\mathrm{diag}(1,-1,-1,-1)$ 为闵氏度规。 拉格朗日量也可以被简化:
\mathcal L=\pm\left[i\psi^\dagger(\partial_0+\vec \alpha\cdot\nabla)\psi-m\psi\dagger\beta\psi\right]=\pm[i\bar\psi\beta\partial_0\psi+i\bar\psi\beta\vec\alpha\cdot\nabla\psi-m\bar\psi\psi]=\pm\bar\psi(i\gamma\mu\partial_\mu-m)\psi
这个形式已经非常简单了,而且一眼就能看出这是一个标量。但是Feynman不满足于此,他观察到很多地方需要将gamma矩阵与别的矢量做内积,因此引入了Feynman斜线符号:
\gamma^\mu a_\mu\equiv a \mkern -9.5 mu /
其有一些性质:
(a \mkern -9.5 mu /)^2=a_\mu a_\nu\gamma\mu\gamma\nu=\dfrac12 a_\mu a_\nu{\gamma\mu,\gamma\nu}=\eta^{\mu\nu}a_\mu a_\nu=a^2
{a \mkern -9.5 mu /,b\mkern -9.5mu /}=2a^\mu b_\mu=2a\cdot b
现在Lagrangian可以写为:
\mathcal L=\pm\bar\psi(i\partial\mkern-9.5mu/-m)\psi
Dirac方程即为:
(i\partial\mkern-9.5mu/-m)\psi=0
这个形式真是非常简洁了。最后证明一下由Dirac方程可以推出KG方程:左作用 $(i\partial\mkern-9.5mu/+m)$ 可得:
(i\partial\mkern-9.5mu/+m)(i\partial\mkern-9.5mu/-m)\psi=-[\partial\mkern-9.5mu/2+m2]=-(,\square2+m2)\psi=0
正是质量为 $m$ 的Klein-Gordan方程。 ## 20.5 旋量的双线性型 我们看到,$\bar\psi\psi$ 是Lorentz标量,而 $\bar\psi\gamma^\mu\psi$ 是Lorentz矢量。那么,其他的双线性型呢?我们将会看到,16个线性独立的双线性型可以分类为:标量(1),极矢量(4),反对称二阶张量(6),赝矢量(4)、赝标量(1)。 首先来考虑比较好处理的**宇称**,Dirac旋量在宇称作用下的变换关系为:
P:\psi(\vec x,t)\to\beta\psi(-\vec x,t)=\gamma^0\psi(-\vec x,t)
gamma矩阵的Dirac伴随性质为:
\overline{\gamma0}=\beta\beta\dagger\beta=\gamma^0,\quad \overline{\gammai}=\beta\vec\alpha\dagger\beta\dagger\beta=\beta\vec\alpha=\gammai
因此gamma矩阵是∗∗Dirac自伴∗∗的。故
\overline{\gamma^0\psi(-\vec x,t)}=\bar\psi(-\vec x,t)\gamma^0
因此
P:\bar\psi\psi(\vec x,t)\to\bar\psi\gamma0\gamma0\psi(-\vec x,t)=\bar\psi\psi(-\vec x,t)
因此这是一个(真)标量(scalar)。另一方面:
P:\bar\psi\gamma^\mu\psi(\vec x,t)\to\bar\psi\gamma0\gamma\mu\gamma^0\psi(-\vec x,t)=\begin{cases} \bar\psi\gamma^0\psi(-\vec x,t),\quad \mu=0\ -\bar\psi\gamma^i\psi(-\vec x,t),\quad\mu\ne 0 \end
其中用到了 $\{\gamma^\mu,\gamma^0\}=2\delta^{\mu 0}$。因此宇称会反转矢量的空间部分,故这是一个**极矢量** (polar vector)。 接下来,我们考虑二阶张量 $\bar\psi\gamma^\mu\gamma^\nu\psi$,其可分解为对称与反对称部分:
\bar\psi\gamma\mu\gamma\nu\psi=\dfrac12\bar\psi{\gamma\mu,\gamma\nu}\psi+\dfrac12\bar\psi[\gamma\mu,\gamma\nu]\psi=\eta{\mu\nu}\bar\psi\psi+\dfrac12\bar\psi[\gamma\mu,\gamma^\nu]\psi
我们按习惯定义反对称部分:
\sigma^{\mu\nu}\equiv\dfrac i2 [\gamma\mu,\gamma\nu]$$ σμν/2 其实就是在Coleman QFT - Lecture 18. Lorentz群的表示中说明的一般的Lorentz变换的生成元 Jμν 的表示。这源于 γμ 的矢量性:
D(Λ)γμD(Λ)=D−1(Λ)γμD(Λ)=Λνμγν
由
D(Λ(ϵμν))=exp(−2iϵμνJμν)
应用无穷小变换 Λνμ=I+ϵνμ 得到:
[Jρσ,γμ]=i(ημργσ−ημσγρ)
进行计算:
[[γρ,γσ],γμ]=γρ{γσ,γμ}−{γρ,γμ}γσ−γσ{γρ,γμ}+{γσ,γμ}γρ=4(ημργσ−ημσγρ)
因此 Jμν=σμν/2 满足条件。不管怎样,ψˉσμνψ 是一个Lorentz变换下的反对称二阶张量。 那如果有三个gamma矩阵呢?其会不会组合出三阶张量呢?这是不会的,或者说你的确构造出了一个三阶张量,但是其可以分解为更低阶的张量。这源于四阶全反对称张量 ϵμνρσ 是Lorentz变换下的不变量。因此我们可以定义下面的对象:
γ5≡iγ0γ1γ2γ3=(0110)
其性质为
(γ5)2=1,γ5†=γ5=−γ5,{γ5,γμ}=0
iψˉγ5ψ 是一个Lorentz标量,但是是一个赝标量 (pseudoscalar)。因为
P:iψˉγ5ψ→iψˉγ0γ5γ0ψ=−iψˉγ5ψ
因此我们说明了四个gamma矩阵凑出的四阶张量实际上会分解出一个赝标量。回到三个gamma矩阵的情形,其一定可以写为
ψˉγμγ5ψ
这是Lorentz变换下的矢量。其在宇称下的变换为:
P:ψˉγμγ5ψ→ψˉγ0γμγ5γ0ψ={−ψˉγ0γ5ψ,μ=0ψˉγiγ5ψ,μ=0
反而是时间分量反号而空间分量不变,因此这是一个赝矢量 (pseudovector)。现在,我们完成了对双线性型的全部分类,列表如下:

任何一个双线性型均可以分解为上面这五种类型的线性组合。
20.6 正交性与完备性
最后,我们来看Dirac方程的平面波解具有的正交性与完备性,这在后一章经常会用到。首先,我们在20.2节构造了下面的矢量:
(up(r)†up(s),up(r)†αup(s))=(vp(r)†vp(s),vp(r)†αvp(s))=2(Ep,p)δrs
这可以用gamma矩阵表示为
\bar u_\vec p^{(r)}\gamma^\mu u_\vec p^{(s)}=\bar v_\vec p^{(r)}\gamma^\mu v_\vec p^{(s)}=2p^\mu\delta^{rs}
将其与 pμ 做内积得到:
\bar u_\vec p^{(r)}p\mkern-8.5mu/ u_\vec p^{(s)}=\bar v_\vec p^{(r)}p\mkern-8.5mu/ v_\vec p^{(s)}=2m^2\delta^{rs}
利用Dirac方程:
(p\mkern-8.5mu/-m)u_\vec p^{(r)}=0,\quad (p\mkern-8.5mu/+m)v_\vec p^{(s)}=0
得到
\bar u_\vec p^{(r)}u_\vec p^{(s)}=2m\delta^{rs},\quad \bar v_\vec p^{(r)}v_\vec p^{(s)}=-2m\delta^{rs}
那混合积 \bar v_\vec p^{(r)}u_\vec p^{(s)} 和 \bar u_\vec p^{(r)}v_\vec p^{(s)} 呢?利用
m\bar u_\vec p^{(r)}v_\vec p^{(s)}=\bar u_\vec p^{(r)}p\mkern-8.5mu /v_\vec p^{(s)}=\bar u_\vec p^{(r)}(-m)v_\vec p^{(s)}\Rightarrow \bar u_\vec p^{(r)}v_\vec p^{(s)}=0
因此
\bar u_\vec p^{(r)}v_\vec p^{(s)}=\bar u_\vec p^{(r)}\gamma^\mu v_\vec p^{(s)}=\bar v_\vec p^{(r)}u_\vec p^{(s)}=\bar v_\vec p^{(r)}\gamma^\mu u_\vec p^{(s)}=0
这样,我们知道四个平面波解是Dirac正交的。接下来研究完备性,考虑下面的定义:
P_u\equiv\dfrac{1}{2m}\sum_{r=1}^2u_\vec p^{(r)}\bar u_\vec p^{(r)},\quad P_v\equiv\dfrac{1}{2m}\sum_{r=1}^2v_\vec p^{(r)}\bar v_\vec p^{(r)}
可验证:
P_u u_\vec p^{(s)}=\dfrac{1}{2m}\sum_{r=1}^2u_\vec p^{(r)}\underbrace{\bar u_\vec p^{(r)}u_\vec p^{(s)}}_{2m\delta^{rs}}=u_\vec p^{(s)},\quad P_u v_\vec p^{(s)}=0,\quad P_v u_{\vec p}^{(s)}=0,P_v v_\vec p^{(s)}=v_\vec p^{(s)}
因此 Pu,Pv 称为旋量投影算符。具体计算表明:
Pu=2m1(p/+m),Pv=−2m1(p/−m)
这时上面的投影关系就很显然了,只需代入Dirac方程即可看出。做一些性质的检验:
Pu+Pv=2m1(p/+m−p/+m)=1
Pu2=4m21(p/2+m2+2mp/)=4m21(2m2+2mp/)=2m1(p/+m)=Pu
Pv2=4m21(p/2+m2−2mp/)=4m21(2m2−2mp/)=−2m1(p/−m)=Pv