第14-15章 自旋与角动量耦合:SU(2)的世界
"自旋不是旋转,但它比旋转更根本。"
前置知识:SU(2)群与旋转
在进入自旋的物理世界之前,我们需要先建立SU(2)群的数学直觉。这是理解自旋-1/2以及所有角动量理论的基石。
A.1 什么是SU(2)群?
SU(2)(Special Unitary group of degree 2,二维特殊酉群)的数学定义是:
SU(2)={U∈C2×2:U†U=I,detU=1}
即所有 2×2 复幺正矩阵中行列式为1的集合。任何SU(2)矩阵都可以显式写成:
U=(a−b∗ba∗),∣a∣2+∣b∣2=1
这意味着SU(2)作为流形同构于三维球面 S3——一个紧致、单连通的流形。
群的四大性质:
- 封闭性:若 U1,U2∈SU(2),则 U1U2∈SU(2)
- 单位元:I=diag(1,1)
- 逆元:U−1=U†
- 结合律:矩阵乘法天然满足
A.2 SU(2)的泡利参数化
在物理中,最实用的参数化使用泡利矩阵:
U(\hat{\mathbf{n}}, \theta) = \exp\left(-i\frac{\theta}{2}\hat{\mathbf{n}}\cdot\boldsymbol{\sigma}\right) = \cos\frac{\theta}{2} \, I - i\sin\frac{\theta}{2} \, (\hat{\mathbf{n}}\cdot\boldsymbol{\sigma})
其中 n^=(nx,ny,nz) 是单位矢量,θ∈[0,4π)。注意:θ=2π 时 U=−I,θ=4π 时才回到 +I。
A.3 SU(2)与SO(3)的2对1同态
SO(3) 是三维空间旋转群,由 3×3 实正交矩阵构成(RTR=I, detR=1)。
从SU(2)到SO(3)的映射可以显式构造。对于 U∈SU(2),定义:
Rij(U)=21Tr(σiUσjU†)
可以验证:
- R 是实矩阵
- RTR=I
- detR=1
- R(−U)=R(U)(2对1)
李代数同构:SU(2)和SO(3)的李代数是同构的。SU(2)的李代数由 {iσj/2} 张成:
[2iσi,2iσj]=−ϵijk2iσk
这与角动量对易关系 [Ji,Jj]=iℏϵijkJk 在结构上完全一致(令 Ji=ℏσi/2)。
A.4 为什么SU(2)比SO(3)更根本?
从拓扑上看:
- SO(3)是不可单连通的:π1(SO(3))=Z2
- SU(2)是单连通的:π1(SU(2))={e}
SU(2)是SO(3)的万有覆盖群(universal covering group)。这意味着SU(2)的表示可以涵盖SO(3)的所有表示,还包括SO(3)没有的双值表示——这正是半整数自旋的来源。
在量子力学中,态矢量允许整体相位不确定性:∣ψ⟩∼eiϕ∣ψ⟩。因此 −I 不是单位元,而是非平凡的群元素。旋转 2π 后波函数变号:
D^(n^,2π)∣ψ⟩=−∣ψ⟩
这在经典物理中没有对应——SO(3)中旋转 2π 就是恒等操作。SU(2)才是量子力学中旋转的"正确"数学描述。
graph TD
A["SU(2): 单连通"] --> B["$\pi_1 = \{e\}$"]
B --> C[覆盖群]
C --> D[包含双值表示]
D --> E[半整数自旋]
F["SO(3): 非单连通"] --> G["$\pi_1 = \mathbb{Z}_2$"]
G --> H[基础群]
H --> I[仅单值表示]
I --> J[仅整数自旋]
C -.->|覆盖| H
E -.->|量子力学需要| D
故事场景:星际导航员的罗盘
公元2189年,深空探测船"羲和号"的导航员小林遇到了一个前所未有的困境。飞船穿越一片强磁场星云时,所有的经典陀螺仪全部失效——不是损坏,而是它们开始给出"不可能"的读数:同一个陀螺仪在不同时刻测量,竟会稳定地指向两个完全相反的方向。地球上的量子物理学家通过量子通讯链路查看了数据后,在屏幕上写下了一行公式:
S^z∣↑⟩=+2ℏ∣↑⟩,S^z∣↓⟩=−2ℏ∣↓⟩
"你们的陀螺仪不是坏了,"物理学家说,"它们只是第一次接触到了电子的自旋——一个纯粹的量子自由度。它不对应任何经典旋转,却比任何机械转动更根本。在SU(2)的世界里,'上’和’下’可以同时存在。"
小林盯着那个公式,突然意识到:人类花了三千年学会用罗盘辨别方向,却要到量子时代才真正理解"方向"本身的含义。
这就是Shankar第14-15章要带我们进入的世界——自旋与角动量耦合的SU(2)王国。
14.1 自旋:一个纯粹的量子自由度
14.1.1 Stern-Gerlach实验与自旋的发现
1922年,Otto Stern和Walther Gerlach让一束银原子通过非均匀磁场, expecting a continuous smear on the detection plate. Instead, they saw two distinct spots.
graph LR
subgraph "Stern-Gerlach实验"
A[银原子源] --> B[非均匀磁场]
B --> C{测量结果}
C -->|经典预期| D[连续分布]
C -->|实际观测| E[两个离散斑点]
E --> F["自旋向上 $m_s=+1/2$"]
E --> G["自旋向下 $m_s=-1/2$"]
end
style D fill:#ffcccc
style E fill:#ccffcc这个实验揭示了一个惊人的事实:银原子(最外层一个电子)的某种内禀属性只能取两个离散值,而不是经典物理预期的连续值。这种内禀属性被命名为自旋(spin),但它不是任何经典的旋转——电子被证明是点粒子,没有可旋转的"内部结构"。
自旋是纯粹的量子力学概念,具有以下特征:
- 内禀角动量:类似于轨道角动量 L^,但完全独立于空间运动
- 量子化:对于自旋-s粒子,Sz的可能取值为 msℏ,其中 ms=−s,−s+1,…,s−1,s
- 电子自旋:s=1/2,只有两个态:∣↑⟩(ms=+1/2)和 ∣↓⟩(ms=−1/2)
graph TD
A[经典旋转] --> B[任意角度可能]
A --> C[角动量连续取值]
D[量子自旋] --> E[仅离散值]
D --> F[内禀属性]
D --> G[非空间运动]
H["电子自旋1/2"] --> I[仅两个态]
H --> J["$m_s = \pm 1/2$"]
E -.->|对比| B
I -.->|特例| E 14.1.2 自旋态空间与泡利矩阵
自旋-1/2粒子的态空间是二维复向量空间。我们选择 S^z 的本征态作为基:
∣↑⟩=(10),∣↓⟩=(01)
任意自旋态可以写成:
∣χ⟩=α∣↑⟩+β∣↓⟩=(αβ)
其中 ∣α∣2+∣β∣2=1。
自旋算符 S^=(S^x,S^y,S^z) 满足与角动量相同的对易关系:
[S^i,S^j]=iℏϵijkS^k
在 Sz 基下,自旋算符表示为 泡利矩阵(Pauli matrices)乘以 ℏ/2:
S^x=2ℏσx=2ℏ(0110)
S^y=2ℏσy=2ℏ(0i−i0)
S^z=2ℏσz=2ℏ(100−1)
泡利矩阵的完整性质:
- 厄米性:σi†=σi(保证本征值为实数)
- 平方为单位矩阵:σi2=I(本征值为 ±1)
- 反对易关系:{σi,σj}=2δijI
- 循环对易:[σi,σj]=2iϵijkσk
- 迹性质:Tr(σi)=0,Tr(σiσj)=2δij
- 行列式:det(σi)=−1
graph LR
subgraph 泡利矩阵的性质
A["$\sigma_x = \begin{pmatrix}0&1\\1&0\end{pmatrix}$"] --> B["本征值 $\pm 1$"]
C["$\sigma_y = \begin{pmatrix}0&-i\\i&0\end{pmatrix}$"] --> B
D["$\sigma_z = \begin{pmatrix}1&0\\0&-1\end{pmatrix}$"] --> B
B --> E["$S_i = \frac{#quot;#quot;\hbar#quot;#quot;}{2}\sigma_i$"]
E --> F["自旋本征值 $\pm \hbar/2$"]
end泡利矩阵之间还有一个优美的乘法关系: $$\sigma_i \sigma_j = \delta_{ij} I + i\epsilon_{ijk}\sigma_k$$ 这可以统一地写成:如果 $\mathbf{a}$ 和 $\mathbf{b}$ 是两个三维矢量,那么 $$(\boldsymbol{\sigma}\cdot\mathbf{a})(\boldsymbol{\sigma}\cdot\mathbf{b}) = \mathbf{a}\cdot\mathbf{b} \, I + i\boldsymbol{\sigma}\cdot(\mathbf{a}\times\mathbf{b})$$ 其中 $\boldsymbol{\sigma} = (\sigma_x, \sigma_y, \sigma_z)$。 ### 14.1.3 自旋的物理图像:不是旋转,但比旋转更根本 为什么叫"自旋"如果它不是旋转?这个名字来自历史——早期物理学家试图用电子绕自身轴旋转来解释。但简单的计算表明:如果电子是一个半径约为经典电子半径的球体,要产生观测到的磁矩,其表面速度需要**超过光速**。 自旋的真正身份是在相对论性量子力学(狄拉克方程)中才完全揭示的:自旋是**洛伦兹群的表示**的自然组成部分。它不是"某种旋转",而是粒子在旋转群(更准确地说是SU(2))下的变换性质。
graph TD
A["经典图像: 小球自转"] --> B[表面超光速矛盾]
B --> C[经典图像失败]
D["量子图像: 内禀自由度"] --> E[与空间运动无关]
E --> F["Stern-Gerlach证明其存在"]
G["相对论图像: Dirac方程"] --> H[自旋是洛伦兹表示的一部分]
H --> I["自动出现,无需假设"]
C -.->|启发| D
D -.->|深化| G一个更深刻的理解来自Wigner的定理:粒子的内禀性质由**庞加莱群(Poincaré group)的不可约表示**分类。对于质量非零的粒子,这些表示由两个量子数标记:**质量** $m$ 和 **自旋** $s$。自旋-$1/2$粒子的态空间在旋转下按照SU(2)群的**二维表示**变换。 --- ## 14.2 旋转与SU(2)群 ### 14.2.1 旋转算符D(R)的显式构造 在量子力学中,物理系统的**主动旋转**由作用在态上的**旋转算符** $\hat{D}(R)$ 实现。对于绕单位矢量 $\hat{\mathbf{n}}$ 转角度 $\theta$ 的旋转,旋转算符为: $$\hat{D}(\hat{\mathbf{n}}, \theta) = \exp\left(-\frac{i}{\hbar}\theta \, \hat{\mathbf{n}}\cdot\hat{\mathbf{J}}\right)$$ 其中 $\hat{\mathbf{J}}$ 是总角动量算符(可以是轨道 $\hat{\mathbf{L}}$、自旋 $\hat{\mathbf{S}}$ 或两者之和)。 对于纯自旋-$1/2$系统,$\hat{\mathbf{J}} = \hat{\mathbf{S}} = \frac{\hbar}{2}\boldsymbol{\sigma}$,因此: $$\hat{D}(\hat{\mathbf{n}}, \theta) = \exp\left(-\frac{i\theta}{2}\hat{\mathbf{n}}\cdot\boldsymbol{\sigma}\right)$$ 展开这个指数。利用 $(\hat{\mathbf{n}}\cdot\boldsymbol{\sigma})^2 = I$(因为 $\hat{\mathbf{n}}$ 是单位矢量),我们得到**显式表达式**: $$\hat{D}(\hat{\mathbf{n}}, \theta) = \cos\frac{\theta}{2} \, I - i\sin\frac{\theta}{2} \, (\hat{\mathbf{n}}\cdot\boldsymbol{\sigma})$$ 写成矩阵形式。若 $\hat{\mathbf{n}} = (\sin\theta'\cos\phi, \sin\theta'\sin\phi, \cos\theta')$: $$\hat{D} = \begin{pmatrix} \cos\frac{\theta}{2} - i\cos\theta'\sin\frac{\theta}{2} & -i\sin\theta'\sin\frac{\theta}{2}e^{-i\phi} \\ -i\sin\theta'\sin\frac{\theta}{2}e^{i\phi} & \cos\frac{\theta}{2} + i\cos\theta'\sin\frac{\theta}{2} \end{pmatrix}$$ 这是一个 $2\times 2$ **幺正矩阵**,行列式为1。所有这样的矩阵构成 **SU(2)群**。
graph TD
A["三维旋转群 SO(3)"] --> B["旋转 $R(#quot;#quot;\hat{n}, \theta#quot;#quot;)$"]
B --> C["3×3正交矩阵"]
D["SU(2)群"] --> E["2×2幺正矩阵"]
E --> F["$U^\dagger U = I$"]
E --> G["$\det U = 1$"]
H["旋转算符 $\hat{D}$"] --> I["自旋1/2: SU(2)"]
H --> J["自旋1: SO(3)的3维表示"]
H --> K["自旋j: (2j+1)维表示"]
I -.->|同态| A
style I fill:#ccffcc
style A fill:#ccffcc 14.2.2 SU(2)与SO(3)的关系:双覆盖
这里有一个深刻而微妙的数学结构。SU(2)群和SO(3)群(三维空间旋转群)不是同一个群,但它们密切相关:
SO(3):三维空间中的旋转,由 3×3 实正交矩阵构成(RTR=I,detR=1)。
SU(2):2×2 复幺正矩阵(U†U=I,detU=1)。
关键关系是:每个SO(3)旋转对应两个SU(2)矩阵。
具体来说,如果我们定义从SU(2)到SO(3)的映射:
U(\hat{\mathbf{n}}, \theta) = \cos\frac{\theta}{2} I - i\sin\frac{\theta}{2} (\hat{\mathbf{n}}\cdot\boldsymbol{\sigma})
那么 U 在自旋态上的作用对应于SO(3)中的一个旋转。但注意:
U(n^,θ+2π)=−U(n^,θ)
然而 U(n^,θ+4π)=+U(n^,θ)。
这意味着旋转 2π(360度)在SU(2)中表示为 −I,而旋转 4π(720度)才回到 +I。
物理意义:自旋-1/2粒子在旋转 2π 后,波函数变号!这不是可观测效应(因为概率 ∣ψ∣2 不变),但在干涉实验中可以检测。这被称为自旋统计定理的一部分。
graph LR
subgraph "SU(2) → SO(3) 双覆盖"
A["$U(#quot;#quot;\hat{n}, \theta#quot;#quot;)$"] -->|映射| B["$R(#quot;#quot;\hat{n}, \theta#quot;#quot;)$"]
C["$U(#quot;#quot;\hat{n}, \theta + 2\pi#quot;#quot;) = -U$"] -->|"同一SO(3)元素"| B
D["$U(#quot;#quot;\hat{n}, \theta + 4\pi#quot;#quot;) = +U$"] -->|回到自身| A
end
style A fill:#ccffcc
style C fill:#ffffcc 14.2.3 旋转下自旋态的变换
考虑一个自旋在 +z 方向的态 ∣↑⟩。如果我们绕 y 轴旋转 θ 角度,得到的新态为:
∣χ′⟩=D^(y^,θ)∣↑⟩=(cos2θI−isin2θσy)(10)
=cos2θ(10)−isin2θ(0i)=⎝⎜⎛cos2θsin2θ⎠⎟⎞
这正是沿 (θ,ϕ=0) 方向自旋的态。如果再做 S^z 测量,得到 ±ℏ/2 的概率分别为 cos2(θ/2) 和 sin2(θ/2)。
14.3 磁共振:Rabi振荡
14.3.1 物理背景:核磁共振与自旋操控
磁共振(Magnetic Resonance)是现代物理学和医学中最重要应用之一。从核磁共振成像(MRI)到量子计算中的量子比特操控,其核心物理都是自旋在时变磁场中的演化。
考虑一个自旋-1/2粒子在磁场中的哈密顿量。经典磁场 B 与自旋的相互作用为:
\hat{H} = -\boldsymbol{\mu}\cdot\mathbf{B} = -\gamma \hat{\mathbf{S}}\cdot\mathbf{B}
其中 γ 是旋磁比(gyromagnetic ratio),对于电子 γe=−gee/(2me)≈−2μB/ℏ(μB 是玻尔磁子,μB=eℏ/(2me)≈9.27×10−24 J/T)。
对于静止在 z 方向的磁场 B0:
H^0=−γB0S^z=2ℏω0σz
其中 拉莫尔频率 ω0=−γB0。两个自旋态的能量差为 ℏω0。
graph TD
A["静磁场 $B_0$"] --> B["拉莫尔频率 $\omega_0$"]
B --> C["能级分裂 $\Delta E = \hbar\omega_0$"]
C --> D["基态: $|\downarrow\rangle$"]
C --> E["激发态: $|\uparrow\rangle$"]
F["射频场 $B_1(t)$"] --> G["共振条件 $\omega = \omega_0$"]
G --> H[自旋翻转]
H --> I[Rabi振荡] 14.3.2 Rabi模型与Rabi频率
现在加上一个垂直于 z 轴的射频场(radio-frequency field),在 xy 平面内以频率 ω 旋转:
B(t)=B0z^+B1(cosωtx^+sinωty^)
总哈密顿量:
H^(t)=2ℏω0σz+2ℏω1(cosωtσx+sinωtσy)
其中 ω1=−γB1 是Rabi频率(表征射频场的强度)。
这个含时问题可以通过旋转参考系(rotating frame)来简化。定义变换:
∣χ~(t)⟩=exp(2iωtσz)∣χ(t)⟩
在新参考系中,哈密顿量变为时间无关的:
H~=2ℏ(ω0−ωω1ω1−(ω0−ω))=2ℏΔσz+2ℏω1σx
其中 失谐量 Δ=ω0−ω。
graph LR
subgraph "旋转参考系变换"
A["实验室系: $\hat{H}(t)$ 含时"] --> B[旋转参考系]
B --> C["$\tilde{H}$ 时间无关"]
C --> D["$\tilde{H} = \frac{#quot;#quot;\hbar\Omega#quot;#quot;}{2}(\cos\phi \, \sigma_z + \sin\phi \, \sigma_x)$"]
D --> E["有效频率 $\Omega = \sqrt{#quot;#quot;\Delta^2 + \omega_1^2#quot;#quot;}$"]
end
style C fill:#ccffcc有效哈密顿量描述了一个**有效磁场**在 $xz$ 平面内: $$\tilde{H} = \frac{\hbar\Omega}{2}(\cos\phi \, \sigma_z + \sin\phi \, \sigma_x)$$ 其中 $\Omega = \sqrt{\Delta^2 + \omega_1^2}$,$\tan\phi = \omega_1/\Delta$。 ### 14.3.3 Rabi振荡的解 假设初始时自旋处于 $|\downarrow\rangle$(低能态)。在旋转参考系中,自旋围绕有效磁场以频率 $\Omega$ 进动。 回到实验室系,我们发现自旋在 $|\uparrow\rangle$ 和 $|\downarrow\rangle$ 之间**周期性振荡**: $$P_{\uparrow}(t) = \frac{\omega_1^2}{\Omega^2}\sin^2\left(\frac{\Omega t}{2}\right)$$ $$P_{\downarrow}(t) = 1 - P_{\uparrow}(t)$$
graph TD
A["初始: $|\downarrow\rangle$"] --> B[时间演化]
B --> C["$t = \pi/\Omega$: 最大混合"]
B --> D["$t = \pi/\omega_1$ (共振): $|\uparrow\rangle$"]
B --> E["$t = 2\pi/\omega_1$: 回到 $|\downarrow\rangle$"]
C --> F["概率: $P_\uparrow = \omega_1^2/\Omega^2$"]
D --> G["$\pi$-脉冲: 完全翻转"]
E --> H["$2\pi$-脉冲: 完整周期"]**共振情况**($\omega = \omega_0$,即 $\Delta = 0$): - 翻转概率:$P_{\uparrow}(t) = \sin^2(\omega_1 t/2)$ - **$\pi$-脉冲**($\omega_1 t = \pi$):完全翻转 $|\downarrow\rangle \to |\uparrow\rangle$ - **$\pi/2$-脉冲**($\omega_1 t = \pi/2$):产生等权重叠加态 $\frac{1}{\sqrt{2}}(|\uparrow\rangle + |\downarrow\rangle)$ **失谐情况**($\Delta \neq 0$): - 最大翻转概率:$P_{\uparrow}^{max} = \omega_1^2/(\Delta^2 + \omega_1^2) < 1$ - 频率升高:$\Omega > \omega_1$ Rabi振荡是量子调控的核心工具——从NMR到量子计算中的量子门,都依赖于精确控制这些振荡。 ### 14.3.4 数值例子:核磁共振中的Rabi频率 **例题**:一个质子(自旋-$1/2$)在MRI扫描仪中,静磁场 $B_0 = 1.5$ T(特斯拉)。射频线圈产生一个垂直振荡磁场 $B_1 = 10^{-4}$ T。计算: 1. **拉莫尔频率** $\omega_0$: 质子的旋磁比 $\gamma_p = 2.675 \times 10^8$ rad/(s·T) $$\omega_0 = \gamma_p B_0 = 2.675 \times 10^8 \times 1.5 = 4.01 \times 10^8 \text{ rad/s}$$ 对应频率 $f_0 = \omega_0/(2\pi) \approx 63.9$ MHz(这正是MRI中常用的射频频率)。 2. **Rabi频率** $\omega_1$: $$\omega_1 = \gamma_p B_1 = 2.675 \times 10^8 \times 10^{-4} = 2.675 \times 10^4 \text{ rad/s}$$ 对应频率 $f_1 \approx 4.26$ kHz。 3. **$\pi$-脉冲时间**: $$t_{\pi} = \frac{\pi}{\omega_1} = \frac{\pi}{2.675 \times 10^4} \approx 1.17 \times 10^{-4} \text{ s} = 117 \text{ }\mu\text{s}$$ 4. **共振时的振荡周期**: $$T_{Rabi} = \frac{2\pi}{\omega_1} \approx 235 \text{ }\mu\text{s}$$ 如果射频场频率偏离 $\omega_0$ 仅 $0.1\%$($\Delta = 4 \times 10^5$ rad/s),则: $$\Omega = \sqrt{\Delta^2 + \omega_1^2} = \sqrt{(4\times 10^5)^2 + (2.675\times 10^4)^2} \approx 4.01 \times 10^5 \text{ rad/s}$$ 最大翻转概率: $$P_{\uparrow}^{max} = \frac{\omega_1^2}{\Omega^2} = \frac{(2.675\times 10^4)^2}{(4.01\times 10^5)^2} \approx 0.0045$$ 仅为 $0.45\%$!这说明MRI中必须极其精确地调谐射频频率到拉莫尔频率。 --- ## 前置知识:Clebsch-Gordan系数表的使用方法 在深入角动量耦合的数学结构之前,先学会"读表"是极其实用的技能。CG系数是角动量理论中最常用的数值工具之一。 ### B.1 CG系数的定义与符号 Clebsch-Gordan系数 $\langle j_1, m_1; j_2, m_2 | j, m \rangle$ 是未耦合基 $|j_1, m_1\rangle |j_2, m_2\rangle$ 到耦合基 $|j, m\rangle$ 的变换矩阵元。 **Condon-Shortley相位约定**: - 所有CG系数选为实数 - 约定使得 $\langle j_1, j_1; j_2, j_2 | j_1+j_2, j_1+j_2 \rangle = 1$ - 升降算符的矩阵元为正 ### B.2 查表方法 CG系数表通常按 $(j_1, j_2)$ 对组织。例如: **两个自旋-$1/2$($j_1 = 1/2, j_2 = 1/2$)**: | 耦合态 $|j, m\rangle$ | 未耦合基展开 | |---|---| | $\|1, 1\rangle$ | $\|\uparrow\uparrow\rangle$ | | $\|1, 0\rangle$ | $\frac{1}{\sqrt{2}}(\|\uparrow\downarrow\rangle + \|\downarrow\uparrow\rangle)$ | | $\|1, -1\rangle$ | $\|\downarrow\downarrow\rangle$ | | $\|0, 0\rangle$ | $\frac{1}{\sqrt{2}}(\|\uparrow\downarrow\rangle - \|\downarrow\uparrow\rangle)$ | **$j_1 = 1, j_2 = 1/2$(例如 $L=1$ 电子与 $S=1/2$ 自旋)**: 总角动量 $j = 3/2$ 或 $1/2$。 $|3/2, 3/2\rangle = |1, 1\rangle |1/2, 1/2\rangle$ $|3/2, 1/2\rangle = \sqrt{\frac{2}{3}}|1, 0\rangle |1/2, 1/2\rangle + \sqrt{\frac{1}{3}}|1, 1\rangle |1/2, -1/2\rangle$ $|1/2, 1/2\rangle = \sqrt{\frac{1}{3}}|1, 0\rangle |1/2, 1/2\rangle - \sqrt{\frac{2}{3}}|1, 1\rangle |1/2, -1/2\rangle$ ### B.3 CG系数的正交归一性 CG系数满足两组正交关系: **行正交**:固定 $j, m$,对不同 $m_1, m_2$ 求和 $$\sum_{m_1,m_2}\langle j_1,m_1;j_2,m_2|j,m\rangle\langle j_1,m_1;j_2,m_2|j',m'\rangle = \delta_{jj'}\delta_{mm'}$$ **列正交**:固定 $m_1, m_2$,对不同 $j, m$ 求和 $$\sum_{j,m}\langle j_1,m_1;j_2,m_2|j,m\rangle\langle j_1,m'_1;j_2,m'_2|j,m\rangle = \delta_{m_1m'_1}\delta_{m_2m'_2}$$ ### B.4 对称性性质 CG系数在交换 $j_1 \leftrightarrow j_2$ 时有: $$\langle j_2, m_2; j_1, m_1 | j, m \rangle = (-1)^{j_1+j_2-j} \langle j_1, m_1; j_2, m_2 | j, m \rangle$$ 时间反演($m \to -m$)性质: $$\langle j_1, -m_1; j_2, -m_2 | j, -m \rangle = (-1)^{j_1+j_2-j} \langle j_1, m_1; j_2, m_2 | j, m \rangle$$ 这些对称性在原子物理和核物理的计算中极为有用,可以大幅减少需要记忆的系数数量。 --- ## 15.1 角动量耦合:从单粒子到多粒子 ### 15.1.1 为什么要耦合角动量? 在现实世界中,孤立的一个角动量几乎不存在。原子中的电子同时具有**轨道角动量** $\hat{\mathbf{L}}$(来自绕核运动)和**自旋角动量** $\hat{\mathbf{S}}$(内禀属性)。两者都与磁场耦合,形成**总角动量** $\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}$。 更一般地,两个独立子系统各具有角动量 $\hat{\mathbf{J}}_1$ 和 $\hat{\mathbf{J}}_2$,它们之和 $\hat{\mathbf{J}} = \hat{\mathbf{J}}_1 + \hat{\mathbf{J}}_2$ 也满足角动量对易关系: $$[\hat{J}_i, \hat{J}_j] = [\hat{J}_{1i} + \hat{J}_{2i}, \hat{J}_{1j} + \hat{J}_{2j}] = i\hbar\epsilon_{ijk}\hat{J}_{1k} + i\hbar\epsilon_{ijk}\hat{J}_{2k} = i\hbar\epsilon_{ijk}\hat{J}_k$$ 因为不同子系统的算符对易:$[\hat{J}_{1i}, \hat{J}_{2j}] = 0$。
graph TD
A[角动量耦合的场景] --> B["轨道+自旋: LS耦合"]
A --> C["两个电子: 总自旋"]
A --> D["核自旋+电子: 超精细结构"]
A --> E["双原子分子: 转动+振动"]
B --> F["$\hat{J} = \hat{L} + \hat{S}$"]
C --> G["$\hat{S}_{total} = \hat{S}_1 + \hat{S}_2$"]
D --> H["$\hat{F} = \hat{J} + \hat{I}$"]
F --> I[原子光谱精细结构]
G --> J[泡利不相容原理]
H --> K[21cm氢线] 15.1.2 两个自旋-1/2的耦合:最简单的例子
考虑两个自旋-1/2粒子(比如两个电子)。每个粒子有自旋态 ∣↑⟩ 和 ∣↓⟩,总系统态空间是4维的张量积空间:
H=H1⊗H2
基矢为:
∣↑↑⟩=∣↑⟩1⊗∣↑⟩2=⎝⎜⎜⎛1000⎠⎟⎟⎞
∣↑↓⟩=⎝⎜⎜⎛0100⎠⎟⎟⎞,∣↓↑⟩=⎝⎜⎜⎛0010⎠⎟⎟⎞
∣↓↓⟩=⎝⎜⎜⎛0001⎠⎟⎟⎞
总自旋算符 S^=S^1+S^2。计算 S^2 和 S^z 在这个基下的矩阵:
S^z=S^1z+S^2z=2ℏ(σz⊗I+I⊗σz)
S^2=S^12+S^22+2S^1⋅S^2=43ℏ2I+43ℏ2I+2S^1⋅S^2
总自旋点积可以写成:
S^1⋅S^2=S^1zS^2z+21(S^1+S^2−+S^1−S^2+)
通过直接计算(Shankar在第15章有详细推导),我们发现4个态可以重新组合为总自旋本征态:
三重态(triplet,s=1):
∣1,1⟩=∣↑↑⟩,ms=+1
∣1,0⟩=√21(∣↑↓⟩+∣↓↑⟩),ms=0
∣1,−1⟩=∣↓↓⟩,ms=−1
单态(singlet,s=0):
∣0,0⟩=√21(∣↑↓⟩−∣↓↑⟩),ms=0
验证:S^2∣1,ms⟩=1(1+1)ℏ2∣1,ms⟩=2ℏ2∣1,ms⟩,S^2∣0,0⟩=0。
graph TD
A["两个自旋1/2"] --> B[4维张量积空间]
B --> C["未耦合基: $|\uparrow\uparrow\rangle, |\uparrow\downarrow\rangle, |\downarrow\uparrow\rangle, |\downarrow\downarrow\rangle$"]
B --> D["耦合基: 总自旋本征态"]
D --> E["三重态 $s=1$: 3个态"]
D --> F["单态 $s=0$: 1个态"]
E --> G["$|1,1\rangle, |1,0\rangle, |1,-1\rangle$"]
F --> H["$|0,0\rangle$"]
C -.->|"Clebsch-Gordan变换"| D 15.1.3 Clebsch-Gordan系数:一般情况与显式计算
对于两个任意角动量 j1 和 j2 的耦合,总角动量 j 的取值范围为:
j=∣j1−j2∣,∣j1−j2∣+1,…,j1+j2
这是角动量加法规则(triangle rule),可以从经典矢量合成的图像直观理解。
总态数验证:
j=∣j1−j2∣∑j1+j2(2j+1)=(2j1+1)(2j2+1)
未耦合基的态 ∣j1,m1⟩⊗∣j2,m2⟩ 到耦合基的 ∣j,m⟩ 的变换由 Clebsch-Gordan系数 给出:
∣j,m⟩=m1,m2∑⟨j1,m1;j2,m2∣j,m⟩∣j1,m1⟩∣j2,m2⟩
其中 m=m1+m2(因为 J^z=J^1z+J^2z)。
graph LR
subgraph "Clebsch-Gordan变换"
A[未耦合基] -->|CG系数| B[耦合基]
C["$|j_1,m_1\rangle|j_2,m_2\rangle$"] --> D["$|j,m\rangle$"]
E["$(2j_1+1)(2j_2+1)$个态"] --> F[相同数量的态]
F --> G[按总j分块]
end
style B fill:#ccffccCG系数有明确的解析表达式(由Racah公式给出),满足正交归一性:
m1,m2∑⟨j1,m1;j2,m2∣j,m⟩⟨j1,m1;j2,m2∣j′,m′⟩=δjj′δmm′
j,m∑⟨j1,m1;j2,m2∣j,m⟩⟨j1,m1′;j2,m2′∣j,m⟩=δm1m1′δm2m2′
15.1.4 数值例子:j1=1,j2=1/2 的CG系数显式计算
例题:计算 j1=1(轨道角动量 L=1)和 j2=1/2(自旋 S=1/2)耦合的CG系数。这在原子物理中极为常见——p 轨道电子与自旋耦合形成总角动量 j=3/2 或 1/2。
步骤1:确定总角动量取值
j=∣1−1/2∣,∣1−1/2∣+1=1/2,3/2
- j=3/2:m=3/2,1/2,−1/2,−3/2,共4个态
- j=1/2:m=1/2,−1/2,共2个态
- 总计 4+2=6=(2×1+1)(2×1/2+1) ✓
步骤2:用递推法构造最高权态
∣3/2,3/2⟩ 只能来自 m1=1,m2=1/2:
∣3/2,3/2⟩=∣1,1⟩∣1/2,1/2⟩
即 ⟨1,1;1/2,1/2∣3/2,3/2⟩=1。
步骤3:用降算符 J^−=J^1−+J^2− 生成其他态
J^−∣3/2,3/2⟩=ℏ√(3/2+1/2)(3/2−1/2+1)∣3/2,1/2⟩=ℏ√3∣3/2,1/2⟩
左边作用在未耦合态上:
J^−∣1,1⟩∣1/2,1/2⟩=J^1−∣1,1⟩∣1/2,1/2⟩+∣1,1⟩J^2−∣1/2,1/2⟩
=ℏ√2∣1,0⟩∣1/2,1/2⟩+ℏ∣1,1⟩∣1/2,−1/2⟩
因此:
∣3/2,1/2⟩=√32∣1,0⟩∣1/2,1/2⟩+√31∣1,1⟩∣1/2,−1/2⟩
即:
- ⟨1,0;1/2,1/2∣3/2,1/2⟩=√2/3
- ⟨1,1;1/2,−1/2∣3/2,1/2⟩=√1/3
步骤4:用正交性构造 j=1/2 的态
∣1/2,1/2⟩ 必须与 ∣3/2,1/2⟩ 正交,且是 ∣1,0⟩∣1/2,1/2⟩ 和 ∣1,1⟩∣1/2,−1/2⟩ 的线性组合:
∣1/2,1/2⟩=√31∣1,0⟩∣1/2,1/2⟩−√32∣1,1⟩∣1/2,−1/2⟩
即:
- ⟨1,0;1/2,1/2∣1/2,1/2⟩=√1/3
- ⟨1,1;1/2,−1/2∣1/2,1/2⟩=−√2/3
步骤5:验证归一性和正交性
∣3/2,1/2⟩ 的系数平方和:2/3+1/3=1 ✓
∣1/2,1/2⟩ 的系数平方和:1/3+2/3=1 ✓
正交性:√2/3⋅√1/3+√1/3⋅(−√2/3)=0 ✓
完整的CG系数表(j1=1,j2=1/2):
| j | m | m1=1,m2=−1/2 | m1=1,m2=1/2 | m1=0,m2=−1/2 | m1=0,m2=1/2 | m1=−1,m2=−1/2 | m1=−1,m2=1/2 |
|---|
| 3/2 | 3/2 | — | 1 | — | — | — | — |
| 3/2 | 1/2 | √1/3 | — | √2/3 | — | — | — |
| 3/2 | -1/2 | — | √2/3 | — | √1/3 | — | — |
| 3/2 | -3/2 | — | — | — | √1/3 | — | √2/3 |
| 1/2 | 1/2 | −√2/3 | — | √1/3 | — | — | — |
| 1/2 | -1/2 | — | −√1/3 | — | √2/3 | — | — |
这些系数在原子物理中直接决定光谱线的强度和偏振。例如,p 电子的自旋-轨道耦合将 L=1 能级分裂为 j=3/2(四重态)和 j=1/2(双重态),这正是碱金属原子光谱双线结构的来源。
15.2 张量算符与Wigner-Eckart定理
15.2.1 什么是张量算符?
在角动量理论中,物理量按其在旋转下的变换性质分类。一个不可约张量算符 T^(k) 是秩为 k 的 (2k+1) 个算符的集合 {T^q(k):q=−k,…,k},在旋转下按照角动量 k 的表示变换:
D^(R)T^q(k)D^†(R)=q′∑Dq′q(k)(R)T^q′(k)
例子:
- 标量(k=0):T^0(0),旋转不变。如 J^2。
- 矢量(k=1):三个分量 T^−1(1),T^0(1),T^+1(1)。如位置 r^ 或角动量 J^ 的球分量。
- 球谐函数 Ylm(θ,ϕ):本身就是秩 l 的张量。
张量算符与角动量的对易关系(定义性质):
[J^±,T^q(k)]=ℏ√k(k+1)−q(q±1)T^q±1(k)
[J^z,T^q(k)]=ℏqT^q(k)
graph TD
A[物理量在旋转下的分类] --> B["标量: $k=0$"]
A --> C["矢量: $k=1$"]
A --> D["张量: $k=2,3,...$"]
B --> E[旋转不变]
C --> F[三个分量]
D --> G["$(2k+1)$个分量"]
E --> H["如 $\hat{J}^2$"]
F --> I["如 $\hat{r}, \hat{p}$"]
G --> J[如四极矩] 15.2.2 Wigner-Eckart定理:选择定则的源泉
Wigner-Eckart定理 是角动量理论中最强大的工具之一。它指出:张量算符 T^q(k) 在角动量本征态之间的矩阵元可以分解为:
⟨j′,m′∣T^q(k)∣j,m⟩=⟨j,m;k,q∣j′,m′⟩√2j′+1⟨j′∣∣T^(k)∣∣j⟩
其中:
- ⟨j,m;k,q∣j′,m′⟩ 是 Clebsch-Gordan系数
- ⟨j′∣∣T^(k)∣∣j⟩ 是 约化矩阵元(reduced matrix element),与 m,m′,q 无关
物理意义:
- 选择定则:如果CG系数为零,矩阵元就为零。CG系数非零的条件是 m′=m+q 且 ∣j−k∣≤j′≤j+k。这就是原子光谱中的选择定则(selection rules)。
- 比例关系:不同 (m,q,m′) 组合的矩阵元之比完全由CG系数决定,与具体物理系统无关。
graph TD
A["Wigner-Eckart定理"] --> B[矩阵元分解]
B --> C["CG系数: 几何部分"]
B --> D["约化矩阵元: 动力学部分"]
C --> E[仅依赖角动量量子数]
D --> F[依赖具体物理系统]
E --> G["选择定则 $\Delta j, \Delta m$"]
G --> H["如电偶极跃迁: $\Delta l = \pm 1$"] 15.2.3 应用:电偶极跃迁的选择定则
考虑原子中电子在电场作用下的跃迁。电偶极算符 d^=−er^ 是**秩1(矢量)**张量算符。
Wigner-Eckart定理告诉我们,跃迁矩阵元 ⟨n′,l′,m′∣r^q∣n,l,m⟩ 非零的条件:
- m 选择定则:m′=m+q,其中 q=−1,0,+1(对应 x±iy 和 z)
- l 选择定则:∣l−1∣≤l′≤l+1,即 l′=l−1,l,l+1。但由于宇称,实际上 l′=l±1
- 自旋:Δs=0(电偶极不改变自旋)
这些选择定则解释了为什么原子光谱是线状谱而非连续谱,以及为什么只有某些跃迁(allowed transitions)强度高。
数值例子:CG系数在原子光谱中的应用
例题:钠原子(Na)的 3p→3s 跃迁(D线,589 nm)。3p 态有 l=1,3s 态有 l=0。考虑自旋-轨道耦合后:
- 3p3/2:j=3/2,四重态(mj=±3/2,±1/2)
- 3p1/2:j=1/2,双重态(mj=±1/2)
问题:从 3p3/2 到 3s1/2(j=1/2)的跃迁中,计算CG系数决定的相对强度。
电偶极算符是秩1张量,k=1。j=3/2→j′=1/2,k=1。
CG系数 ⟨3/2,m;1,q∣1/2,m′⟩ 非零的条件:
- m′=m+q
- ∣3/2−1∣≤1/2≤3/2+1,即 1/2≤1/2≤5/2 ✓
具体系数(查表或递推):
对于 m′=1/2:
- m=3/2,q=−1:⟨3/2,3/2;1,−1∣1/2,1/2⟩=−√1/2
- m=1/2,q=0:⟨3/2,1/2;1,0∣1/2,1/2⟩=√1/6
- m=−1/2,q=+1:⟨3/2,−1/2;1,1∣1/2,1/2⟩=−√1/3
相对跃迁强度正比于CG系数的平方:
| 初态 m | 终态 m′ | ∣CG∣2 |
|—|—|—|
| 3/2 | 1/2 | 1/2 |
| 1/2 | 1/2 | 1/6 |
| -1/2 | 1/2 | 1/3 |
| -3/2 | -1/2 | 1/2 |
| -1/2 | -1/2 | 1/6 |
| 1/2 | -1/2 | 1/3 |
对所有初态求平均(假设初态各 mj 等概率占据):
∣CG∣2=41(21+61+31+21)=41×63+1+2+3=249=83
类似地,对于 3p1/2→3s1/2:
∣CG∣2=21(31+31)=31
两者强度比:
I(3p1/2→3s1/2)I(3p3/2→3s1/2)=(2j+1)∣CG∣21/2(2j+1)∣CG∣23/2=2×1/34×3/8=2/33/2=49=2.25
实验上,钠D线的 D2(589.0 nm,3p3/2→3s1/2)与 D1(589.6 nm,3p1/2→3s1/2)强度比约为 2:1,与上述理论估算一致。这个比例是CG系数和统计权重共同决定的结果,与具体的原子波函数细节(约化矩阵元)无关——这正是Wigner-Eckart定理的威力。
本章总结
graph TD
A["第14-15章: 自旋与角动量耦合"] --> B["自旋1/2的SU(2)世界"]
A --> C["角动量耦合: 从单粒子到多粒子"]
A --> D[张量算符与选择定则]
B --> B1["泡利矩阵 $\sigma_x, \sigma_y, \sigma_z$"]
B --> B2["旋转算符 $D(#quot;#quot;\hat{n}, \theta#quot;#quot;) = \exp(#quot;#quot;-i\theta \hat{n}\cdot\sigma/2#quot;#quot;)$"]
B --> B3["SU(2)双覆盖SO(3)"]
B --> B4["Rabi振荡: 自旋操控"]
C --> C1["两个自旋1/2: 三重态+单态"]
C --> C2["Clebsch-Gordan系数"]
C --> C3["总角动量 $j = |j_1-j_2|, ..., j_1+j_2$"]
D --> D1["不可约张量算符 $T^{(k)}_q$"]
D --> D2["Wigner-Eckart定理"]
D --> D3[电偶极跃迁选择定则]
B2 -.->|应用| B4
C1 -.->|推广| C2
C2 -.->|基础| D2这两章的核心洞见:自旋不是旋转,但它是比旋转更根本的数学结构。SU(2)群不仅描述了自旋-1/2粒子的行为,也是理解所有角动量耦合、选择定则和量子调控的基础。从核磁共振到量子计算,从原子光谱到粒子物理,SU(2)的语言无处不在。
练习与思考
1. 泡利矩阵的恒等式
证明对于任意两个矢量 a 和 b:
(\boldsymbol{\sigma}\cdot\mathbf{a})(\boldsymbol{\sigma}\cdot\mathbf{b}) = \mathbf{a}\cdot\mathbf{b} \, I + i\boldsymbol{\sigma}\cdot(\mathbf{a}\times\mathbf{b})
并由此推导:如果 a 和 b 平行,(\boldsymbol{\sigma}\cdot\mathbf{a})^2 = a^2 I;如果垂直,(\boldsymbol{\sigma}\cdot\mathbf{a})(\boldsymbol{\sigma}\cdot\mathbf{b}) = i\boldsymbol{\sigma}\cdot(\mathbf{a}\times\mathbf{b})。
提示:利用 σiσj=δijI+iϵijkσk。
2. SU(2)到SO(3)的映射验证
考虑 U=exp(−iθσy/2)=cos(θ/2)I−isin(θ/2)σy。证明这个 2×2 幺正矩阵作用于泡利矩阵的方式:
UσxU†=σxcosθ+σzsinθ
UσzU†=−σxsinθ+σzcosθ
这与SO(3)中绕 y 轴旋转 θ 的变换一致。为什么这说明了SU(2)到SO(3)的2对1同态?
3. Rabi振荡的量子门应用
在量子计算中,自旋-1/2(或等效的两能级系统)是最简单的量子比特。证明:
- 一个 π-脉冲(ω1t=π,共振)实现的是NOT门(∣↓⟩↔∣↑⟩)
- 一个 π/2-脉冲实现的是Hadamard门的类似操作(将基态变为等权重叠加态)
- 如果在失谐条件下(Δ≫ω1),系统近似只经历相位演化,可以实现相位门
设计一个操作序列来实现从 ∣↓⟩ 到 √21(∣↓⟩+i∣↑⟩) 的变换。
4. CG系数显式计算
用递推法推导 j1=1,j2=1 耦合的CG系数。总角动量 j=2,1,0。验证:
- ∣2,2⟩=∣1,1⟩∣1,1⟩
- ∣2,1⟩=√21(∣1,0⟩∣1,1⟩+∣1,1⟩∣1,0⟩)
- ∣1,1⟩=√21(∣1,0⟩∣1,1⟩−∣1,1⟩∣1,0⟩)
- ∣0,0⟩=√31(∣1,1⟩∣1,−1⟩−∣1,0⟩∣1,0⟩+∣1,−1⟩∣1,1⟩)
5. Wigner-Eckart定理与选择定则
证明:对于电偶极跃迁(k=1),如果初态和终态的宇称相同,则矩阵元必为零。这解释了为什么 Δl=0 的跃迁是禁戒的(电偶极算符是奇宇称,所以 l 必须改变奇数个单位)。