转动算符的表示
转动算符的矩阵元是 Wigner 函数
我们在 Dj 空间中定义 D(R)(其中 R 由 n^ 和 α 指定)的矩阵元为:
Dm′m(j)(R)=⟨jm′∣D^(R)∣jm⟩=⟨jm′∣e−ℏiJ⋅n^α∣jm⟩
这就是 Wigner 函数,D(j)(R) 是一个 m=2j+1 维的矩阵。你可能会想问,为什么不同的态只需要改变一个 m 指标而不需要同时改变 jm 指标呢,例如 ∣jm′⟩ 而不是 ∣j′m′⟩?这是因为转动算符不会改变角动量的值 j,它只会改变角动量在某一方向上的分量 m:
J2[D(R)∣jm⟩]=D(R)J2∣jm⟩=j(j+1)ℏ2[D(R)∣jm⟩]
因此,转动算符只能在给定的 j 子空间内起作用。由上式还可以看出,J2 在 ∣jm⟩ 下是 2j+1 重简并的。
在文献中,由 Dm′mj(R) 构成的 (2j+1)×(2j+1) 矩阵被称为旋转算符 D(R) 的 2j+1 维不可约表示。这意味着在这个 2j+1 维的子空间 Dj 内,没有任何更小的子空间能在所有旋转操作下保持“封闭”。换句话说,旋转操作会将 Dj 中的任意一个态混合成该空间内所有其他态的线性组合,而不会跑出这个空间,也不会局限在更小的部分里。
可以证明:矩阵 {Dj(R)} 构成了 {D(R)} 的一个表示。这意味着这些矩阵保留了群的乘法结构:即旋转算符的乘积对应于矩阵的乘积:
D(R′)D(R)=D(R′R)
Dj(R1)Dj(R2)=Dj(R1R2)
欧拉转动的矩阵表示
类比经典欧拉转动 R(α,β,γ)=Rz(α)Ry(β)Rz(γ),只需直接将经典转动中的每个转动替换为相应的量子转动算符:
D(α,β,γ)≡Dz^(α)Dy^(β)Dz^(γ)=e−ℏiJzαe−ℏiJyβe−ℏiJzγ
与之对应的 Wigner 函数为
Dm′m(j)(α,β,γ)=⟨jm′∣e−ℏiJzαe−ℏiJyβe−ℏiJzγ∣jm⟩=e−i(m′α+mγ)⟨jm′∣e−ℏiJyβ∣jm⟩
唯一非平凡的部分是 e−ℏiJyβ,由于 Jy 是 J+ 和 J− 的线性组合,其将不同的 m 值混合在一起,因此我们定义
dm′m(j)(β)≡⟨jm′∣e−ℏiJyβ∣jm⟩
称为 Wigner 小 d 矩阵。
自旋 1/2 的转动算符
自旋 1/2
自然界利用某些关系的最低维(维数=2)实现并不是先验的,但大量的实验(如原子光谱学、核磁共振等)证明事实确实如此。
我们在 j=1/2 情况下,将之前理论中的 J 替换为 S,m 只能取两个值:m=+1/2 和 m=−1/2。
SO(3) 和 SU(2)
根据转动算符 Dn^(α)=e−ℏiS⋅n^α,有
Dn^(1/2)(α)=e−ℏiS⋅n^α=e−2iσ⋅n^α=cos2α−isin2ασ⋅n^=cos2α−isin2αcosθsinθeiϕsinθe−iϕ−cosθ=cos2α−isin2αcosθ−isin2αsinθeiϕ−isin2αsinθe−iϕcos2α+isin2αcosθ≡A−B∗BA∗
其中 A=cos2α−isin2αcosθ,B=−isin2αsinθe−iϕ。可以验证 Dn^(1/2)(α) 的行列式为 1:
dete−2iσ⋅n^α=∣A∣2+∣B∣2=1
且由于 σ 为厄米算符,所以 e−2iσ⋅n^α 是幺正矩阵。
我们定义 SU(2) 群:如果一个 2×2 复矩阵 U 满足
- 幺正性:UU†=U†U=1
- 幺模性:detU=1
那么我们将 U 称为 SU(2) 群(2维特殊幺正群)的一个元素。满足这两个条件的所有矩阵构成了 SU(2) 群。U 的一般形式为
U(a,b)=(a−b∗ba∗),∣a∣2+∣b∣2=1
a 和 b 各自由两个实参数,加上 detU=1 的约束条件,SU(2) 群有三个自由参数,我们可以验证
U(a,b)U†(a,b)=(a−b∗ba∗)(a∗b∗−ba)=(aa∗+bb∗−b∗a∗+a∗b∗−ab+bab∗b+a∗a)=1
任意的转动算符 Dn^(21)(α)=e−2iσ⋅n^α 都属于 SU(2) 群。
我们再回到经典转动 R^(α)=e−iαI⋅n^,这是一个实正交的矩阵,且 detR(α)=1。我们称 Rn^(α) 为 SO(3) 群(特殊正交群)的一个元素。SO(3) 群有三个自由参数,与 SU(2) 群一样。
任意经典转动在自旋 1/2 的态空间中都有对应的量子转动,具有如下对应关系
Rn^(α)=e−iαI⋅n^∈SO(3)⟺Dn^(21)(α)=e−2iσ⋅n^α∈SU(2)
但是这不是一一对应的
Rn^(α+2π)Dn^(21)(α+2π)=Rn^(α)Rn^(2π)=Rn^(α)=Dn^(21)(α)Dn^(21)(2π)=−Dn^(21)(α)
对于所有的半奇数都有:Dn^(j)(α+2π)=Dn^(j)(α)Dn^(j)(2π)=−Dn^(j)(α)
因此,SU(2) 群是 SO(3) 群的二重覆盖群。换句话说,SU(2) 群中的每一个元素都对应 SO(3) 群中的一个元素,但 SO(3) 群中的每一个元素对应 SU(2) 群中的两个元素,分别相差一个负号。
由转动算符得到自旋分量算符的本征态
我们之前通过本征值问题求得了自旋分量算符的本征态,现在我们从一个动态的转动观点直接得到这些本征态。
我们想要找到一个联系 ∣ψ⟩=∣↑⟩ 和 ∣ψR⟩=∣↑⟩n^ 的旋转算符 D(R)。由前面的内容可知,σ 的平均值就是 n^ 方向的单位向量
⟨σ⟩R=n^⟨↑∣σ∣↑⟩n^=n^
这个式子对所有的 n^ 都成立,当选择 n^=z^ 时,有
⟨σ⟩=⟨↑∣σ∣↑⟩=z^
又因为 σ 与 S 成正比,所以同样满足
D†(R)σD(R)=Rσ⇒⟨σ⟩R=R⟨σ⟩
其中 ⟨σ⟩R 和 ⟨σ⟩ 分别是转动后和转动前在 ∣↑⟩ 上的平均值。
任意轴都可以由 z^ 通过转动得到:Rz^(ϕ)Ry^(θ)z^=n^
所以相应的 R 和 D(R) 分别为
R=Rz^(ϕ)Ry^(θ)⟶D(R)=Dz^(ϕ)Dy^(θ)
现在问题就简单了,我们不用再求本征值了,只需将 ∣↑⟩ 通过转动算符 D(R) 转动到 ∣↑⟩n^ 方向:
∣↑⟩n=D(R)∣↑⟩=(cos2ϕ−isin2ϕσ3)(cos2θ−isin2θσ2)(01)=(e−iϕ/200eiϕ/2)cos2θsin2θ−sin2θcos2θ(10)=(e−iϕ/200eiϕ/2)cos2θsin2θ=e−iϕ/2cos2θeiϕ/2sin2θ
为了找到两个态之间的对应,我们先根据σ 的平均值就是 n^ 方向的单位向量、平均值在转动情况下表现得与经典向量一致这两个性质,再从平均值出发,找到对应的算符的转动。
自旋共振:拉比问题
拉比公式
我们知道,只要将一个磁矩 μ 放在一个磁场 B 中,它就会进动。电子的反常磁矩为
μ=−g2meeS
这个磁矩与磁场的相互作用哈密顿量为
H=−μ⋅B=2megeS⋅B=4megeℏσ⋅B
其中 g 被称为朗德 g 因子。虽然狄拉克方程预测 g 严格等于 2,但实验值只是非常接近 2,量子电动力学修正了狄拉克的预测,这些修正可以按精细结构常数 α 的幂级数进行计算,该差异 a=21(g−2) 被称为粒子的反常磁矩,电子的单圈贡献为 ae=2πα。
下面,我们简单地令 g=2,那么
H=2meeℏσ⋅B=μBσ⋅B
其中 μB=2meeℏ 是玻尔磁子。
拉比问题:量子力学解法
拉比问题是少数可以精确求解的含时问题。考虑一个由静态磁场 B 和一个频率为 ωr 的旋转磁场 Brot 组成的磁场:
B=Brot(cosωrte1+sinωrte2)+Be3
一个自旋粒子在其中的哈密顿量为
H(t)=mee[Brot(cosωrtS1+sinωrtS2)+BS3]≡ωrot(cosωrtS1+sinωrtS2)+ω0S3
其中 ωrot=meeBrot 和 ω0=meeB。我们之前讲过含时幺正变换,其能够将含时哈密顿量转化为一个新的不含时哈密顿量,选取幺正变换 W(t)=eℏiS3ωrt,则
H′(t)≡W(t)HW†(t)+iℏdtdW(t)W†(t)=ωrot[cosωrt(S1cosωrt−S2sinωrt)+sinωrt(S1sinωrt+S2cosωrt)]+(ω0−ωr)S3=ωrotS1+(ω0−ωr)S3
为什么在转动坐标系下哈密顿量不含时了?因为我们选择的幺正变换 W(t) 实际上是一个绕 S3 轴以角速度 ωr 旋转的参考系,在这个参考系中,原本随时间变化的磁场和参考系是同步转动的,所以看起来是静止的,因此哈密顿量也变得不含时了。最后再使用 W†(t) 将解变回实验室参考系即可。
现在,不含时问题就可以直接做时间演化了。时间演化算符为
U′(t)=e−ℏiH′t=e−ℏi[ωrotS1+(ω0−ωr)S3]t=e−2i[ωrotσ1+(ω0−ωr)σ3]t
若定义
Ω=21(ω0−ωr)2+ωrot2n^=[ωrote1+(ω0−ωr)e3]/(2Ω)
则时间演化算符为
U′(t)=e−2i[ωrotσ1+(ω0−ωr)σ3]/(2Ω)2Ωt=e−2iσ⋅n^(2Ωt)=cosΩt−isinΩtσ⋅n^=cosΩt−isinΩtσ⋅n^=cosΩt−isinΩt[ωrotσ1+(ω0−ωr)σ3]/(2Ω)
设初态为 ∣ψ′(0)⟩=W(t)∣ψ(0)⟩=∣↑⟩,则在转动参考系下的时间演化为
∣ψ′(t)⟩=U′(t)∣ψ′(0)⟩=cosΩt∣↑⟩−isinΩt[ωrotσ1∣↑⟩+(ω0−ωr)σ3∣↑⟩]/(2Ω)=cosΩt∣↑⟩−isinΩt[ωrot∣↓⟩+(ω0−ωr)∣↑⟩]/(2Ω)=[cosΩt−isinΩt(ω0−ωr)/(2Ω)]∣↑⟩−isinΩtωrot/(2Ω)∣↓⟩=2Ω1e−iΩt2ω0−ωr+Ω2ωrot−eiΩt2ω0−ωr−Ω2ωrot
最后,将解变回实验室参考系(薛定谔绘景):
∣ψ(t)⟩=W†(t)∣ψ′(t)⟩=[cosΩt−isinΩt(ω0−ωr)/(2Ω)]e−2iωrt∣↑⟩−isinΩtωrot/(2Ω)e2iωrt∣↓⟩
我们计算初态为 ∣↑⟩ 时,自旋在任意时刻 t 处于 ∣↓⟩ 态的概率(翻转概率):
P↓(t)=∣⟨↓∣ψ(t)⟩∣2=∣−isinΩtωrot/(2Ω)e2iωrt∣2=8Ω21ωrot2(1−cos2Ωt)
这就是拉比公式。上式中转动变换产生的影响 e2iωrt 并没有起实际作用,说明转动参考系下的 Sz 的平均值和实验室参考系下的平均值是相同的。
总结:拉比公式
自旋磁矩在振荡磁场中的哈密顿量为
H(t)=mee[Brot(cosωrtS1+sinωrtS2)+BS3]
若定义
ωrot≡meeBrot,ω0≡meeB,Ω=21(ω0−ωr)2+ωrot2
则初态为 ∣↑⟩ 时,自旋在任意时刻 t 处于 ∣↓⟩ 态的概率为
P↓(t)=8Ω21ωrot2(1−cos2Ωt)
拉比振荡
可以看出,在 t=2Ωπ 时,翻转概率达到最大值
P↓max=P↓(t=2Ωπ)=4Ω21ωrot2=(ω0−ωr)2+ωrot2ωrot2=(1−ωr/ω0)2+(ωrot/ω0)2(ωrot/ω0)2
当外加磁场的频率 ωr 接近自旋在静态磁场 B 中的拉莫尔进动频率 ω0 时,翻转概率最大值能够达到 1,这就是拉比共振现象。此时,系统以频率 ωrot 在两个自旋态之间振荡,这种现象被称为拉比振荡。下图展示了不同失谐情况下 ωr 下的最大翻转概率
转动坐标系下拉比振荡时间演化的物理图像如下
在转动参考系下,转动磁场变为一个静态磁场 Brot=Brotex,而静态磁场在转动小减小但方向不变,根据哈密顿量可知其变为有效磁场 Beff=(B−emeωr)ez=(emeω0−emeωr)ez。因此,自旋在转动参考系下的总磁场为 Btot=Brot+Beff,自旋围绕 Btot 进动;而在实验室参考系下,自旋围绕 Btot 进动的同时,Brot 还在以角速度 ωr 绕 ez 旋转,这就解释了为什么转动参考系中 Sz 的平均值和实验室参考系下的平均值是相同的,因为自旋的平均值在 z 方向上的投影并不取决于 Brot 是否绕 z 轴旋转。
在共振情况下 ωr=ω0,有效磁场 Beff=0,转动磁场 Brot 不再是斜的,而是完全在 x 方向上,此时自旋围绕 x 轴进动,如下图所示。
此时若初态为 ∣↑⟩,则自旋是能够完全达到相反位置的 ∣↓⟩ 的,因此最大值为 1,而在失谐情况下,自旋围绕斜着的 Btot 进动,无法完全到达 ∣↓⟩ 态。
共振情形
在共振情况下 ωr=ω0,Ω=2ωrot=2meeBrot,初态为 ∣↑⟩ 时,翻转概率为
P↓(t)=21(1−cosωrott)
这可以用于描述二能级系统在外加谐场驱动下的动力学行为
其会在两个能级之间以频率 ωrot 来回振荡,这就是拉比振荡。
角动量的合成
角动量的合成并不神秘,仅仅是一个表象变换而已。其在原子光谱学、核物理与粒子物理以及量子磁学等方面有着重要的应用,且能更好地阐述基底变换的概念。
两种合成
例1:一个自旋 1/2 粒子:∣x′,↑/↓⟩=∣x′⟩⊗∣↑/↓⟩
因为轨道运动和自旋运动是独立的,故在由 {∣x′⟩} 张成的空间中的任意算符,都与由 {∣↑⟩,∣↓⟩} 张成的二维空间中的任何算符对易。特别地,轨道角动量 L 与自旋算符 S 对易。总角动量是这两者的和
J=L+S=L⊗1S+1L⊗S
很容易验证 J 也是一个角动量算符,其三个分量满足角动量算符的基本对易关系
[Jx,Jy]=[Lx+Sx,Ly+Sy]=[Lx,Ly]+[Sx,Sy]=iℏLz+iℏSz=iℏJz
对于轨道角动量,空间部分的基矢为
L2∣nlm⟩=l(l+1)ℏ2∣nlm⟩,Lz∣nlm⟩=mℏ∣nlm⟩
对于自旋角动量,自旋部分的基矢为
S2∣↑/↓⟩=43ℏ2∣↑/↓⟩,Sz∣↑/↓⟩=±2ℏ∣↑/↓⟩
我们已经说过角动量的合成就是基矢的变换,那么需要找到一组完备的基矢,使得它们是两个对易子的共同本征态。我们选取两套不同的相互对易的算符,这些可观测量分别为
Set I:L2,LZ,S2,SZSet II:J2,JZ,L2,S2
例2:两个自旋 1/2 粒子:
自旋算符应为
S=S1+S2=S1⊗12+11⊗S2
其中 S1 和 S2 各自满足各自的角动量分量基本对易关系
[Siα,Sjβ]=iℏδijεαβγSiγ(i,j=1,2且α,β=x,y,z)
同样可以验证总角动量 S 也是一个角动量算符,其满足角动量分量的基本对易关系
[Sx,Sy]=[S1x+S2x,S1y+S2y]=[S1x,S1y]+[S2x,S2y]=iℏS1z+iℏS2z=iℏSz
在这个例子中,两套不同的相互对易的算符分别为
Set I:S12,S1z,S22,S2zSet II:S2,Sz,S12,S22
设 set I 部分的基矢为
∣↑↑⟩,∣↑↓⟩,∣↓↑⟩,∣↓↓⟩
set II 部分的基矢为 ∣Sm⟩
S2∣Sm⟩=S(S+1)ℏ2∣Sm⟩Sz∣Sm⟩=mℏ∣Sm⟩
其中 S 和 m 分别为总自旋和总自旋在 z 方向的分量的量子数。接下来我们计算 S2 的矩阵表示,S2 可写为
S2=S12+2S1⋅S2+S22=2(S1xS2x+S1yS2y+S1zS2z)+23ℏ2=(S1+S2−+S1−S2++2S1zS2z)+23ℏ2
那么有
S2∣↑↑⟩S2∣↑↓⟩S2∣↓↑⟩S2∣↓↓⟩=ℏ2(0+0+2∣↑↑⟩)+23ℏ2∣↑↑⟩=2ℏ2∣↑↑⟩=ℏ2(0+∣↓↑⟩−2∣↑↓⟩)+23ℏ2∣↑↓⟩=ℏ2∣(↓↑⟩+↑↓⟩)=ℏ2(∣↑↓⟩+0−2∣↓↑⟩)+23ℏ2∣↓↑⟩=ℏ2(∣↑↓⟩+∣↓↑⟩)=ℏ2(0+0+2∣↓↓⟩)+23ℏ2∣↓↓⟩=2ℏ2∣↓↓⟩
因此,S2 在 {∣↑↑⟩,∣↑↓⟩,∣↓↑⟩,∣↓↓⟩} 基底下的矩阵表示为
S2≐ℏ22000011001100002
由于四个基底都是 Sz 的本征态,故 Sz 在该基底下的对角矩阵为
Sz≐ℏ100000000000000−1
为了将 S2 对角化,我们找到一组新的基
{∣↑↑⟩,2∣↑↓⟩+∣↓↑⟩,2∣↑↓⟩−∣↓↑⟩,∣↓↓⟩}
在这组新基底下,S2 的矩阵表示对角化,Sz 的矩阵不变
S2≐ℏ22000020000000002Sz≐ℏ100000000000000−1
我们又知道,S2 的本征值为 S(S+1)ℏ2,根据 S2 矩阵的对角元,S(S+1)=2 或 S(S+1)=0,解得 S=1 或 S=0。因此,总自旋量子数 S 的可能取值为 1 和 0。其中当 S=1 时,分别对应基矢中的第一个、第二个和第四个;当 S=0 时,对应基矢中的第三个。
接下来我们直接从矩阵角度获得 S2 和 Sz 的矩阵表示,定义单个自旋 1/2 的基底为
∣↑⟩≐(01)∣↓⟩≐(10)
两个自旋 1/2 的基底为它们的张量积
∣↑↑⟩≐(01)⊗(01)=1000,∣↓↑⟩≐(10)⊗(01)=0010,∣↑↓⟩≐(01)⊗(10)=0100∣↓↓⟩≐(10)⊗(10)=0001
根据 S1⊗12 和 11⊗S2 可得
S1z≐2ℏ(100−1)⊗(1001)=2ℏ1000010000−10000−1S2z≐(1001)⊗2ℏ(100−1)=2ℏ10000−1000010000−1
因此,z 分量总自旋算符为
Sz=S1z+S2z≐ℏ100000000000000−1
综上,两个自旋 1/2 的角动量耦合后会得到两个不同的总自旋量子数 S,分别为 1 和 0,对应的基矢分别为
∣S=1,m=1⟩=∣↑↑⟩∣S=1,m=0⟩=21(∣↑↓⟩+∣↓↑⟩)∣S=1,m=−1⟩=∣↓↓⟩⎭⎬⎫(三重态)∣S=0,sz=0⟩=21(∣↑↓⟩−∣↓↑⟩)(单重态)
这个结果可以形象地表示为
21⊗21=1⊕0
两个自旋 1/2 希尔伯特空间的直积,是一个自旋 1 空间与一个自旋 0 空间的直和,是可约的。
角动量的合成问题本质上是一个基底变换问题,我们将乘积态转换至总自旋态。
角动量的叠加
考虑两个属于独立自由度的角动量 J1 和 J2,量子数分别为 j1 和 j2,它们的对易子为零:
[Jiα,Jjβ]=iℏδijεαβγJiγ(i,j=1,2α,β=x,y,z)
它们各自对应一个无穷小转动,并直积出一个总的转动
(1−ℏiJ1⋅n^ε)⊗(1−ℏiJ2⋅n^ε)=1−ℏi(J1⊗1+1⊗J2)⋅n^ε≡1−ℏiJ⋅n^ε
其中 J≡J1⊗1+1⊗J2=J1+J2,这从无穷小转动的角度证明了我们之前对 J 的定义是合理的。我们之前证明了 J 也是一个角动量算符,满足基本对易关系:
[Jα,Jβ]=[J1α+J2α,J1β+J2β]=[J1α,J1β]+[J2α,J2β]=iℏεαβγJ1γ+iℏεαβγJ2γ=iℏεαβγ(J1γ+J2γ)=iℏεαβγJγ
容易看出 J 是整个系统的生成元,那么总的转动算符为
Dn^(θ)=e−ℏiJ⋅n^θ=e−ℏi(J1+J2)⋅n^θ=e−ℏiJ1⋅n^θe−ℏiJ2⋅n^θ
由于 J1 和 J2 都是标量函数,且 J1⋅J2 也是标量函数(因为只要两个三维空间中的算符具有相同转动规则,它们的点积就是一个标量函数),那么它们在总的转动下是不变的
⎩⎨⎧Dn^†(θ)J12Dn^(θ)=J12,Dn^†(θ)J22Dn^(θ)=J22,Dn^†(θ)J1⋅J2Dn^(θ)=J1⋅J2
因此我们猜测总角动量的每一个分量都跟 J12、J22 和 J1⋅J2 对易,事实验证确实如此,因此我们注意到
J2=J12+J22+2J1⋅J2=J12+J22+(J1+J2−+J1−J2++2J1zJ2z)
里面的三项是两两对易的
[J12,J22]=[J12,J1⋅J2]=[J22,J1⋅J2]=0
所以总角动量跟各分量之间也是对易的 [J2,J12]=[J2,J22]=0。
两种基矢选择
第一种:不合成(直积态),用各自的 Jz投影标定, J12,J1z,J22,J2z 的共同本征矢为:∣j1j2;m1m2⟩≡∣j1m1⟩⊗∣j2m2⟩
Ji2∣j1j2;m1m2⟩=ji(ji+1)ℏ2∣j1j2;m1m2⟩Jiz∣j1j2;m1m2⟩=miℏ∣j1j2;m1m2⟩
第二种:合成后用总的 J2 和 Jz 标定,J2,Jz,J12,J22 的共同本征矢为:∣j1j2;jm⟩
J2∣j1j2;jm⟩=j(j+1)ℏ2∣j1j2;jm⟩Jz∣j1j2;jm⟩=mℏ∣j1j2;jm⟩Ji2∣j1j2;jm⟩=ji(ji+1)ℏ2∣j1j2;jm⟩
对于给定的 j1 和 j2,J12 和 J22 这两个算符在 Dj1⊗Dj2 中都正比于单位算符
Ji2≐ji(ji+1)ℏ21j1⊗1j2
两种基矢之间的幺正变换:克莱布什-戈登系数
在右矢空间 Dj1⊗Dj2 中,我们有完备关系
m1m2∑∣j1j2;m1m2⟩⟨j1j2;m1m2∣=1
我们插入这个完备性关系得到
∣j1j2;jm⟩=m1m2∑∣j1j2;m1m2⟩⟨j1j2;m1m2∣j1j2;jm⟩
上式右边的内积就是两个基矢之间的转换系数,称为克莱布什-戈登系数,其具有以下性质
-
对m 的限制:只有当 m=m1+m2 时,克莱布什-戈登系数才非零:
⟨j1j2;m1m2∣j1j2;jm⟩=0,m=m1+m2
-
对 j 的限制:只有当 ∣j1−j2∣≤j≤j1+j2 时,克莱布什-戈登系数才非零:
⟨j1j2;m1m2∣j1j2;jm⟩=0,j<∣j1−j2∣orj>j1+j2
由于直和运算具有维数相加的性质,因此我们可以将这个范围写作
j1⊗j2=∣j1−j2∣⊕(∣j1−j2∣+1)⊕⋯⊕(j1+j2)
在不失一般性的情况下,我们总可以选择 j1≥j2,那么绝对值符号可以去掉,那么就有 2j2+1 个允许的 j 值。
对于上图最右上角的点 (j=m=j1+j2),它只能对应 (m1=j1,m2=j2),因此
∣j1j2;j1+j2,j1+j2⟩=∣j1j2;j1j2⟩
等式右边的第一个 j1+j2 表示 jmax,第二个 j1+j2 表示 mmax,等式右边的 j1 和 j2 分别表示 m1max 和 m2max。
对于再向下的一条直线 m=mmax−1=j1+j2−1 上的点
j≥m=j1+j2−1⇒j=j1+j2−1orj=j1+j2
那么在 m=j1+j2−1 的前提下,j=j1+j2 或 j=j1+j2−1 对应的基矢表示分别为
∣j1j2;j1+j2,j1+j2−1⟩=a∣j1j2;j1,j2−1⟩+b∣j1j2;j1−1,j2⟩∣j1j2;j1+j2−1,j1+j2−1⟩=c∣j1j2;j1,j2−1⟩+d∣j1j2;j1−1,j2⟩
为什么两式左边会是相同乘积态的线性组合呢?因为这两个红点在同一条直线上,而一条直线上有多少点就表示该直线代表的子空间是多少维的。且同时也代表在固定 m 下 j 的个数。
接下来就要确定系数 a,b,c,d 了。利用已知的最右上角的红点 ∣j1j2;j1+j2,j1+j2⟩=∣j1j2;j1j2⟩,两端分别左乘 J−,得到
J−∣j1j2;j1+j2,j1+j2⟩=(J1−+J2−)∣j1j2;j1j2⟩⟹2(j1+j2)∣j1j2;j1+j2,j1+j2−1⟩=2j1∣j1j2;j1−1,j2⟩+2j2∣j1j2;j1,j2−1⟩
最后就能对比得到 a 和 b 的值
∣j1j2;j1+j2,j1+j2−1⟩=j1+j2j2∣j1j2;j1,j2−1⟩+j1+j2j1∣j1j2;j1−1,j2⟩
为了继续确定 c,d,我们先做如下约定:
-
C-G 系数总是可以被选为实数(因为按照约定,我们已经选取 J± 的矩阵元为实数)。
-
对于 (m1,m2) 态,具有较大 m1 的那一个乘积态的系数被选为正值。
利用上述两个约定,∣j1j2;j1+j2−1,j1+j2−1⟩ 中的系数 c 和 d 可以通过正交归一关系被唯一确定
⟨j1j2;j1+j2−1,j1+j2−1∣j1j2;j1+j2−1,j1+j2−1⟩=1⟨j1j2;j1+j2−1,j1+j2−1∣j1j2;j1+j2,j1+j2−1⟩=0
解得 c 和 d 的值为
c=j1+j2j1,d=−j1+j2j2
将 c 和 d 代入之前的表达式中,得到
∣j1j2;j1+j2−1,j1+j2−1⟩=j1+j2j1∣j1j2;j1,j2−1⟩−j1+j2j2∣j1j2;j1−1,j2⟩
对于 m=mmax=j1+j2−2 的直线上的点,就会存在三个态:允许的 j 值为 j=j1+j2、j=j1+j2−1 和 j=j1+j2−2。由于 J− 只会在子空间 Dj 内作用,只改变 m 的值而不改变 j 的值,因此其中的两个态可直接通过 J− 作用在 m=j1+j2−1 的两个态上得到
∣j1j2;j1+j2,j1+j2−2⟩∼J−∣j1j2;j1+j2,j1+j2−1⟩∣j1j2;j1+j2−1,j1+j2−2⟩∼J−∣j1j2;j1+j2−1,j1+j2−1⟩
类似地,∣j1j2;j1+j2−2,j1+j2−2⟩ 中的 3 个展开系数,可以通过解 1 个归一化条件和 2 个正交关系确定。原则上,我们可以通过类似的方法,生成对应于 m=j1+j2−3,j1+j2−4,⋯ 的 J2 本征态,如下图所示
通过上图我们可以看出,对于固定的 m=m1+m2,我们总可以利用 C-G 系数将 J2 的本征态表示为乘积态的线性组合
∣j1j2;jm⟩=m1+m2=m∑∣j1j2;m1m2⟩⟨j1j2;m1m2∣j1j2;jm⟩
其中 m1+m2=m 共有 2j+1 项,从而共有 2j+1 个乘积态基矢,又因为不同 j 下的 ∣jm⟩ 是两两正交的,所以对应的 j 值也只能有 2j+1 项,因此 j 的最小值为 jmin=(j1+j2)−(2j)=j1−j2。
因为 C-G 系数是两个完备正交归一基底之间的幺正变换系数,且都为实数,所以能构成一个实正交矩阵。显然我们可以准备好所有的 j1,j2 的 C-G 系数表格,然后通过查表法获得任意需要的 C-G 系数,如下图所示
同理我们还可以反过来,将乘积态表示为总角动量态的线性组合
∣j1j2;m1m2⟩=j∑∣j1j2;jm⟩⟨j1j2;jm∣j1j2;m1m2⟩
此时只需对同一 m 对应的 C-G 系数进行转置读取即可。也就是说,若是 ∣j1j2;m1m2⟩→∣j1j2;jm⟩ 的变换那么就竖着读,若是 ∣j1j2;jm⟩→∣j1j2;m1m2⟩ 的变换那么就横着读。
L-S耦合
在原子物理中,电子的总角动量 J 由轨道角动量 L 和自旋角动量 S 组成
J=L+S
此时 j=l±21≡j±,mmax=l+21,mmin=−(l+21)。示意图如下
我们想要得到固定 m 下的一般态的 C-G 系数,那么一般对应两个 m 相同而 j 不同的态
l21;l+21,m⟩l21;l−21,m⟩=al21;m+21,−21⟩+bl21;m−21,21⟩=a′l21;m+21,−21⟩+b′l21;m−21,21⟩
我们有两种方法来确定 C-G 系数 a,b,a′,b′
①定义 m±=m±21,又因为 j±=l±21,结合我们之前得到的 J− 对 ∣jm⟩ 的降低作用
∣jm⟩=(2j)!(j−m)!(j+m)!(J−)j−m∣jj⟩(J−)j−m∣jj⟩=(j+m)!(2j)!(j−m)!∣jm⟩
我们可以得到
l21;j+m⟩=(2j+)!(j+−m)!(j++m)!(J−)j+−ml21;j+j+⟩
又因为最高的 ∣j+j+⟩ 态等于最高的 ∣l21;l21⟩ 态 l21;j+j+⟩=l21;l21⟩,再利用 J−=L−+S−,并利用二项式定理展开 (J−)j+−m
(J−)j+−m=(L−+S−)j+−m=n=0∑j−mCj+−mn(L−)n(S−)j+−m−n
这个二项式定理看起来很多项数,实际上只有两项不为零,因为 S− 只能作用一次就把自旋 1/2 的态降到最低态了,所以当 j+−m−n≥2 时,(S−)j+−m−nl21;l21⟩=0,因此我们只需要考虑 n=j+−m 和 n=j+−m−1 两项,计算下式
l21;j+m⟩=(2j+)!(j+−m)!(j++m)!n=0∑j−mCj+−mn(L−)n(S−)j+−m−nl21;l21⟩
可以得到
l21;j+m⟩=2l+11[l+m+l21;m−,21⟩+l−m−l21;m+,−21⟩]
可以利用同样的方法或者正交归一关系得到
l21;j−m⟩=2l+11[−l−m−l21;m−,21⟩+l+m+l21;m+,−21⟩]
②利用两个态的正交归一关系,得到
a2+b2=a′2+b′2=1aa′+bb′=0
但这只是三个条件,而我们有四个未知数,因此还需要一个额外条件,例如比值 a/b。将之前推导出的
J1⋅J2∣j1j2;jm⟩=2ℏ2[j(j+1)−j1(j1+1)−j2(j2+1)]∣j1j2;jm⟩
应用至 L−S 耦合基矢变换的等式两端,分别得到
2L⋅Sl21;j+m⟩=ℏ2[j+(j++1)−l(l+1)−43]l21;j+m⟩=lℏ2l21;j+m⟩2L⋅Sl21;m+,−21=(L+S−+L−S++2LzSz)l21;m+,−21=ℏ2((l+m+)(l−m++1)l21;m−,21−m+l21;m+,−21)2L⋅Sl21;m−,21⟩=(L+S−+L−S++2LzSz)l21;m−,21⟩=ℏ2((l−m−)(l+m−+1)l21;m+,−21+m−l21;m−,21)
将上面三式代入 l21;j±m⟩=al21;m+,−21⟩+bl21;m−,21⟩ 中,比较系数有
⎩⎨⎧la=b(l−m−)(l+m+)−am+lb=a(l+m+)(l−m−)+bm−
得到了 a,b 的比例关系
ba=l+m+l−m−
凑齐了四个方程,最后也能得到
l21;j+m⟩=2l+11[l+m+l21;m−,21⟩+l−m−l21;m+,−21⟩]
升降算符方向的投影算符
我们已经得到了 J± 作用后的结果分别为
ℏ22L⋅Sl21;j±m⟩=[j±(j±+1)−l(l+1)−43]l21;j±m⟩=⎩⎨⎧ll21;j+m⟩−(l+1)l21;j−m⟩
上式中可以看出,L⋅S 本征值分别为 l 和 −(l+1),对应的本征态分别为 l21;j+m⟩ 和 l21;j−m⟩。因此我们可以定义两个投影算符,将任意态投影到这两个子空间中
我们在之前证明了:一个具有有限个本征值的可观测量满足一个代数方程。事实上,可以证明 L⋅S 在 Dl⊗D1/2 中满足下列二次方程:
(ℏ22L⋅S−l)(ℏ22L⋅S+l+1)=0
之前的定理就是这时的逆定理:一个满足最简代数方程 ϕ(ξ)=∏i=1n(ξ−ci)=0 的厄米算符ξ,拥有 n 个本征值,由 {ci} 给出。但是注意,虽然本征值只有两个,这两个本征值是简并的。
我们将 ((2l+1)⋅2) 个 (j,m) 态按序排成一列
由于这些态都是 L⋅S 的本征态,因此我们可以将 L⋅S 在这个基底下表示为对角阵
ℏ22L⋅S=diag(l,⋯,l,−l−1,...,−l−1)
我们定义两个投影算符
℘+≡m=−j+∑j+l21;j+m⟩⟨l21;j+m=diag(1,⋯,1,0,…,0)℘−≡m=−j−∑j−l21;j−m⟩⟨l21;j−m=diag(0,⋯,0,1,…,1)
它们分别将 j=j+ 和 j=j− 的投影基矢叠加,分别形成 j+ 和 j− 子空间的投影算符,分别可以将 Dl⊗D1/2 中的任意态投影到 j+ 或 j− 子空间中去。这两个算符可以用 L⋅S 来表示
℘+=2l+1ℏ22L⋅S+l+1℘−=2l+1l−ℏ22L⋅S