第14-15章 自旋与角动量耦合：SU(2)的世界

2026-05-16 科普 12.7k 字 65 分钟阅读
📑 目录
  1. 前置知识：SU(2)群与旋转
    1. A.1 什么是SU(2)群？
    2. A.2 SU(2)的泡利参数化
    3. A.3 SU(2)与SO(3)的2对1同态
    4. A.4 为什么SU(2)比SO(3)更根本？
  2. 故事场景：星际导航员的罗盘
  3. 14.1 自旋：一个纯粹的量子自由度
    1. 14.1.1 Stern-Gerlach实验与自旋的发现
    2. 14.1.2 自旋态空间与泡利矩阵
    3. 14.1.3 自旋的物理图像：不是旋转，但比旋转更根本
  4. 14.2 旋转与SU(2)群
    1. 14.2.1 旋转算符D®的显式构造
    2. 14.2.2 SU(2)与SO(3)的关系：双覆盖
    3. 14.2.3 旋转下自旋态的变换
  5. 14.3 磁共振：Rabi振荡
    1. 14.3.1 物理背景：核磁共振与自旋操控
    2. 14.3.2 Rabi模型与Rabi频率
    3. 14.3.3 Rabi振荡的解
    4. 14.3.4 数值例子：核磁共振中的Rabi频率
  6. 前置知识：Clebsch-Gordan系数表的使用方法
    1. B.1 CG系数的定义与符号
    2. B.2 查表方法
    3. B.3 CG系数的正交归一性
    4. B.4 对称性性质
  7. 15.1 角动量耦合：从单粒子到多粒子
    1. 15.1.1 为什么要耦合角动量？
    2. 15.1.2 两个自旋-1/2的耦合：最简单的例子
    3. 15.1.3 Clebsch-Gordan系数：一般情况与显式计算
    4. 15.1.4 数值例子：j1=1,j2=1/2j_1 = 1, j_2 = 1/2j​1​​=1,j​2​​=1/2 的CG系数显式计算
  8. 15.2 张量算符与Wigner-Eckart定理
    1. 15.2.1 什么是张量算符？
    2. 15.2.2 Wigner-Eckart定理：选择定则的源泉
    3. 15.2.3 应用：电偶极跃迁的选择定则
  9. 数值例子：CG系数在原子光谱中的应用
  10. 本章总结
  11. 练习与思考

第14-15章 自旋与角动量耦合：SU(2)的世界

"自旋不是旋转，但它比旋转更根本。"

前置知识：SU(2)群与旋转

在进入自旋的物理世界之前，我们需要先建立SU(2)群的数学直觉。这是理解自旋-1/21/2以及所有角动量理论的基石。

A.1 什么是SU(2)群？

SU(2)（Special Unitary group of degree 2，二维特殊酉群）的数学定义是：

SU(2)={U∈C2×2:U†U=I,detU=1}\text{SU}(2) = \{ U \in \mathbb{C}^{2\times 2} : U^\dagger U = I, \det U = 1 \}

即所有 2×22\times 2 复幺正矩阵中行列式为1的集合。任何SU(2)矩阵都可以显式写成：

U=(ab−b∗a∗),∣a∣2+∣b∣2=1U = \begin{pmatrix} a & b \\ -b^* & a^* \end{pmatrix}, \quad |a|^2 + |b|^2 = 1

这意味着SU(2)作为流形同构于三维球面 S3S^3——一个紧致、单连通的流形。

群的四大性质：

  1. 封闭性：若 U1,U2∈SU(2)U_1, U_2 \in \text{SU}(2)，则 U1U2∈SU(2)U_1 U_2 \in \text{SU}(2)
  2. 单位元：I=diag(1,1)I = \text{diag}(1,1)
  3. 逆元：U−1=U†U^{-1} = U^\dagger
  4. 结合律：矩阵乘法天然满足

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)\hat{\mathbf{n}} = (n_x, n_y, n_z) 是单位矢量，θ∈[0,4π)\theta \in [0, 4\pi)。注意：θ=2π\theta = 2\pi 时 U=−IU = -I，θ=4π\theta = 4\pi 时才回到 +I+I。

A.3 SU(2)与SO(3)的2对1同态

SO(3) 是三维空间旋转群，由 3×33\times 3 实正交矩阵构成（RTR=IR^T R = I, detR=1\det R = 1）。

从SU(2)到SO(3)的映射可以显式构造。对于 U∈SU(2)U \in \text{SU}(2)，定义：

Rij(U)=12Tr(σiUσjU†)R_{ij}(U) = \frac{1}{2}\text{Tr}(\sigma_i U \sigma_j U^\dagger)

可以验证：

  • RR 是实矩阵
  • RTR=IR^T R = I
  • detR=1\det R = 1
  • R(−U)=R(U)R(-U) = R(U)（2对1）

李代数同构：SU(2)和SO(3)的李代数是同构的。SU(2)的李代数由 {iσj/2}\{i\sigma_j/2\} 张成：

[iσi2,iσj2]=−ϵijkiσk2\left[\frac{i\sigma_i}{2}, \frac{i\sigma_j}{2}\right] = -\epsilon_{ijk}\frac{i\sigma_k}{2}

这与角动量对易关系 [Ji,Jj]=iℏϵijkJk[J_i, J_j] = i\hbar\epsilon_{ijk}J_k 在结构上完全一致（令 Ji=ℏσi/2J_i = \hbar\sigma_i/2）。

A.4 为什么SU(2)比SO(3)更根本？

从拓扑上看：

  • SO(3)是不可单连通的：π1(SO(3))=Z2\pi_1(\text{SO}(3)) = \mathbb{Z}_2
  • SU(2)是单连通的：π1(SU(2))={e}\pi_1(\text{SU}(2)) = \{e\}

SU(2)是SO(3)的万有覆盖群（universal covering group）。这意味着SU(2)的表示可以涵盖SO(3)的所有表示，还包括SO(3)没有的双值表示——这正是半整数自旋的来源。

在量子力学中，态矢量允许整体相位不确定性：∣ψ⟩∼eiϕ∣ψ⟩|\psi\rangle \sim e^{i\phi}|\psi\rangle。因此 −I-I 不是单位元，而是非平凡的群元素。旋转 2π2\pi 后波函数变号：

D^(n^,2π)∣ψ⟩=−∣ψ⟩\hat{D}(\hat{\mathbf{n}}, 2\pi)|\psi\rangle = -|\psi\rangle

这在经典物理中没有对应——SO(3)中旋转 2π2\pi 就是恒等操作。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∣↓⟩\hat{S}_z |\uparrow\rangle = +\frac{\hbar}{2}|\uparrow\rangle, \quad \hat{S}_z |\downarrow\rangle = -\frac{\hbar}{2}|\downarrow\rangle

"你们的陀螺仪不是坏了，"物理学家说，"它们只是第一次接触到了电子的自旋——一个纯粹的量子自由度。它不对应任何经典旋转，却比任何机械转动更根本。在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^\hat{\mathbf{L}}，但完全独立于空间运动
  • 量子化：对于自旋-ss粒子，SzS_z的可能取值为 msℏm_s\hbar，其中 ms=−s,−s+1,…,s−1,sm_s = -s, -s+1, \ldots, s-1, s
  • 电子自旋：s=1/2s = 1/2，只有两个态：∣↑⟩|\uparrow\rangle（ms=+1/2m_s=+1/2）和 ∣↓⟩|\downarrow\rangle（ms=−1/2m_s=-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/21/2粒子的态空间是二维复向量空间。我们选择 S^z\hat{S}_z 的本征态作为基：

∣↑⟩=(10),∣↓⟩=(01)|\uparrow\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \quad |\downarrow\rangle = \begin{pmatrix} 0 \\ 1 \end{pmatrix}

任意自旋态可以写成：

∣χ⟩=α∣↑⟩+β∣↓⟩=(αβ)|\chi\rangle = \alpha |\uparrow\rangle + \beta |\downarrow\rangle = \begin{pmatrix} \alpha \\ \beta \end{pmatrix}

其中 ∣α∣2+∣β∣2=1|\alpha|^2 + |\beta|^2 = 1。

自旋算符 S^=(S^x,S^y,S^z)\hat{\mathbf{S}} = (\hat{S}_x, \hat{S}_y, \hat{S}_z) 满足与角动量相同的对易关系：

[S^i,S^j]=iℏϵijkS^k[\hat{S}_i, \hat{S}_j] = i\hbar \epsilon_{ijk} \hat{S}_k

在 SzS_z 基下，自旋算符表示为 泡利矩阵（Pauli matrices）乘以 ℏ/2\hbar/2：

S^x=ℏ2σx=ℏ2(0110)\hat{S}_x = \frac{\hbar}{2}\sigma_x = \frac{\hbar}{2}\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}

S^y=ℏ2σy=ℏ2(0−ii0)\hat{S}_y = \frac{\hbar}{2}\sigma_y = \frac{\hbar}{2}\begin{pmatrix} 0 & -i \\ i & 0 \end{pmatrix}

S^z=ℏ2σz=ℏ2(100−1)\hat{S}_z = \frac{\hbar}{2}\sigma_z = \frac{\hbar}{2}\begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}

泡利矩阵的完整性质：

  1. 厄米性：σi†=σi\sigma_i^\dagger = \sigma_i（保证本征值为实数）
  2. 平方为单位矩阵：σi2=I\sigma_i^2 = I（本征值为 ±1\pm 1）
  3. 反对易关系：{σi,σj}=2δijI\{\sigma_i, \sigma_j\} = 2\delta_{ij}I
  4. 循环对易：[σi,σj]=2iϵijkσk[\sigma_i, \sigma_j] = 2i\epsilon_{ijk}\sigma_k
  5. 迹性质：Tr(σi)=0\text{Tr}(\sigma_i) = 0，Tr(σiσj)=2δij\text{Tr}(\sigma_i\sigma_j) = 2\delta_{ij}
  6. 行列式：det(σi)=−1\det(\sigma_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{"\\"\hbar\\""}{2}\sigma_i$"]
    E --> F["自旋本征值 $\pm \hbar/2$"]
    end

泡利矩阵之间还有一个优美的乘法关系：

σiσj=δijI+iϵijkσk\sigma_i \sigma_j = \delta_{ij} I + i\epsilon_{ijk}\sigma_k

这可以统一地写成：如果 a\mathbf{a} 和 b\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）的不可约表示分类。对于质量非零的粒子，这些表示由两个量子数标记：质量 mm 和 自旋 ss。自旋-1/21/2粒子的态空间在旋转下按照SU(2)群的二维表示变换。

14.2 旋转与SU(2)群

14.2.1 旋转算符D®的显式构造

在量子力学中，物理系统的主动旋转由作用在态上的旋转算符 D^(R)\hat{D}(R) 实现。对于绕单位矢量 n^\hat{\mathbf{n}} 转角度 θ\theta 的旋转，旋转算符为：

D^(n^,θ)=exp(−iℏθn^⋅J^)\hat{D}(\hat{\mathbf{n}}, \theta) = \exp\left(-\frac{i}{\hbar}\theta \, \hat{\mathbf{n}}\cdot\hat{\mathbf{J}}\right)

其中 J^\hat{\mathbf{J}} 是总角动量算符（可以是轨道 L^\hat{\mathbf{L}}、自旋 S^\hat{\mathbf{S}} 或两者之和）。

对于纯自旋-1/21/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（因为 n^\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})

写成矩阵形式。若 n^=(sinθ′cosϕ,sinθ′sinϕ,cosθ′)\hat{\mathbf{n}} = (\sin\theta'\cos\phi, \sin\theta'\sin\phi, \cos\theta')：

D^=(cosθ2−icosθ′sinθ2−isinθ′sinθ2e−iϕ−isinθ′sinθ2eiϕcosθ2+icosθ′sinθ2)\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×22\times 2 幺正矩阵，行列式为1。所有这样的矩阵构成 SU(2)群。

graph TD
    A["三维旋转群 SO(3)"] --> B["旋转 $R("\\"\hat{n}, \theta\\"")$"]
    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×33\times 3 实正交矩阵构成（RTR=IR^T R = I，detR=1\det R = 1）。

SU(2)：2×22\times 2 复幺正矩阵（U†U=IU^\dagger U = I，detU=1\det U = 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})

那么 UU 在自旋态上的作用对应于SO(3)中的一个旋转。但注意：

U(n^,θ+2π)=−U(n^,θ)U(\hat{\mathbf{n}}, \theta + 2\pi) = -U(\hat{\mathbf{n}}, \theta)

然而 U(n^,θ+4π)=+U(n^,θ)U(\hat{\mathbf{n}}, \theta + 4\pi) = +U(\hat{\mathbf{n}}, \theta)。

这意味着旋转 2π2\pi（360度）在SU(2)中表示为 −I-I，而旋转 4π4\pi（720度）才回到 +I+I。

物理意义：自旋-1/21/2粒子在旋转 2π2\pi 后，波函数变号！这不是可观测效应（因为概率 ∣ψ∣2|\psi|^2 不变），但在干涉实验中可以检测。这被称为自旋统计定理的一部分。

graph LR
    subgraph "SU(2) → SO(3) 双覆盖"
    A["$U("\\"\hat{n}, \theta\\"")$"] -->|映射| B["$R("\\"\hat{n}, \theta\\"")$"]
    C["$U("\\"\hat{n}, \theta + 2\pi\\"") = -U$"] -->|"同一SO(3)元素"| B
    D["$U("\\"\hat{n}, \theta + 4\pi\\"") = +U$"] -->|回到自身| A
    end
    
    style A fill:#ccffcc
    style C fill:#ffffcc

14.2.3 旋转下自旋态的变换

考虑一个自旋在 +z+z 方向的态 ∣↑⟩|\uparrow\rangle。如果我们绕 yy 轴旋转 θ\theta 角度，得到的新态为：

∣χ′⟩=D^(y^,θ)∣↑⟩=(cosθ2I−isinθ2σy)(10)|\chi'\rangle = \hat{D}(\hat{y}, \theta)|\uparrow\rangle = \left(\cos\frac{\theta}{2} I - i\sin\frac{\theta}{2}\sigma_y\right)\begin{pmatrix} 1 \\ 0 \end{pmatrix}

=cosθ2(10)−isinθ2(0i)=(cosθ2sinθ2)= \cos\frac{\theta}{2}\begin{pmatrix} 1 \\ 0 \end{pmatrix} - i\sin\frac{\theta}{2}\begin{pmatrix} 0 \\ i \end{pmatrix} = \begin{pmatrix} \cos\frac{\theta}{2} \\ \sin\frac{\theta}{2} \end{pmatrix}

这正是沿 (θ,ϕ=0)(\theta, \phi=0) 方向自旋的态。如果再做 S^z\hat{S}_z 测量，得到 ±ℏ/2\pm\hbar/2 的概率分别为 cos2(θ/2)\cos^2(\theta/2) 和 sin2(θ/2)\sin^2(\theta/2)。

14.3 磁共振：Rabi振荡

14.3.1 物理背景：核磁共振与自旋操控

磁共振（Magnetic Resonance）是现代物理学和医学中最重要应用之一。从核磁共振成像（MRI）到量子计算中的量子比特操控，其核心物理都是自旋在时变磁场中的演化。

考虑一个自旋-1/21/2粒子在磁场中的哈密顿量。经典磁场 B\mathbf{B} 与自旋的相互作用为：

\hat{H} = -\boldsymbol{\mu}\cdot\mathbf{B} = -\gamma \hat{\mathbf{S}}\cdot\mathbf{B}

其中 γ\gamma 是旋磁比（gyromagnetic ratio），对于电子 γe=−gee/(2me)≈−2μB/ℏ\gamma_e = -g_e e/(2m_e) \approx -2\mu_B/\hbar（μB\mu_B 是玻尔磁子，μB=eℏ/(2me)≈9.27×10−24\mu_B = e\hbar/(2m_e) \approx 9.27 \times 10^{-24} J/T）。

对于静止在 zz 方向的磁场 B0B_0：

H^0=−γB0S^z=ℏω02σz\hat{H}_0 = -\gamma B_0 \hat{S}_z = \frac{\hbar\omega_0}{2}\sigma_z

其中 拉莫尔频率 ω0=−γB0\omega_0 = -\gamma B_0。两个自旋态的能量差为 ℏω0\hbar\omega_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频率

现在加上一个垂直于 zz 轴的射频场（radio-frequency field），在 xyxy 平面内以频率 ω\omega 旋转：

B(t)=B0z^+B1(cosωtx^+sinωty^)\mathbf{B}(t) = B_0\hat{z} + B_1(\cos\omega t \, \hat{x} + \sin\omega t \, \hat{y})

总哈密顿量：

H^(t)=ℏω02σz+ℏω12(cosωtσx+sinωtσy)\hat{H}(t) = \frac{\hbar\omega_0}{2}\sigma_z + \frac{\hbar\omega_1}{2}(\cos\omega t \, \sigma_x + \sin\omega t \, \sigma_y)

其中 ω1=−γB1\omega_1 = -\gamma B_1 是Rabi频率（表征射频场的强度）。

这个含时问题可以通过旋转参考系（rotating frame）来简化。定义变换：

∣χ~(t)⟩=exp(iωt2σz)∣χ(t)⟩|\tilde{\chi}(t)\rangle = \exp\left(\frac{i\omega t}{2}\sigma_z\right)|\chi(t)\rangle

在新参考系中，哈密顿量变为时间无关的：

H~=ℏ2(ω0−ωω1ω1−(ω0−ω))=ℏΔ2σz+ℏω12σx\tilde{H} = \frac{\hbar}{2}\begin{pmatrix} \omega_0 - \omega & \omega_1 \\ \omega_1 & -(\omega_0 - \omega) \end{pmatrix} = \frac{\hbar\Delta}{2}\sigma_z + \frac{\hbar\omega_1}{2}\sigma_x

其中 失谐量 Δ=ω0−ω\Delta = \omega_0 - \omega。

graph LR
    subgraph "旋转参考系变换"
    A["实验室系: $\hat{H}(t)$ 含时"] --> B[旋转参考系]
    B --> C["$\tilde{H}$ 时间无关"]
    C --> D["$\tilde{H} = \frac{"\\"\hbar\Omega\\""}{2}(\cos\phi \, \sigma_z + \sin\phi \, \sigma_x)$"]
    D --> E["有效频率 $\Omega = \sqrt{"\\"\Delta^2 + \omega_1^2\\""}$"]
    end
    
    style C fill:#ccffcc

有效哈密顿量描述了一个有效磁场在 xzxz 平面内：

H~=ℏΩ2(cosϕσz+sinϕσx)\tilde{H} = \frac{\hbar\Omega}{2}(\cos\phi \, \sigma_z + \sin\phi \, \sigma_x)

其中 Ω=Δ2+ω12\Omega = \sqrt{\Delta^2 + \omega_1^2}，tanϕ=ω1/Δ\tan\phi = \omega_1/\Delta。

14.3.3 Rabi振荡的解

假设初始时自旋处于 ∣↓⟩|\downarrow\rangle（低能态）。在旋转参考系中，自旋围绕有效磁场以频率 Ω\Omega 进动。

回到实验室系，我们发现自旋在 ∣↑⟩|\uparrow\rangle 和 ∣↓⟩|\downarrow\rangle 之间周期性振荡：

P↑(t)=ω12Ω2sin2(Ωt2)P_{\uparrow}(t) = \frac{\omega_1^2}{\Omega^2}\sin^2\left(\frac{\Omega t}{2}\right)

P↓(t)=1−P↑(t)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$-脉冲: 完整周期"]

共振情况（ω=ω0\omega = \omega_0，即 Δ=0\Delta = 0）：

  • 翻转概率：P↑(t)=sin2(ω1t/2)P_{\uparrow}(t) = \sin^2(\omega_1 t/2)
  • π\pi-脉冲（ω1t=π\omega_1 t = \pi）：完全翻转 ∣↓⟩→∣↑⟩|\downarrow\rangle \to |\uparrow\rangle
  • π/2\pi/2-脉冲（ω1t=π/2\omega_1 t = \pi/2）：产生等权重叠加态 12(∣↑⟩+∣↓⟩)\frac{1}{\sqrt{2}}(|\uparrow\rangle + |\downarrow\rangle)

失谐情况（Δ≠0\Delta \neq 0）：

  • 最大翻转概率：P↑max=ω12/(Δ2+ω12)<1P_{\uparrow}^{max} = \omega_1^2/(\Delta^2 + \omega_1^2) < 1
  • 频率升高：Ω>ω1\Omega > \omega_1

Rabi振荡是量子调控的核心工具——从NMR到量子计算中的量子门，都依赖于精确控制这些振荡。

14.3.4 数值例子：核磁共振中的Rabi频率

例题：一个质子（自旋-1/21/2）在MRI扫描仪中，静磁场 B0=1.5B_0 = 1.5 T（特斯拉）。射频线圈产生一个垂直振荡磁场 B1=10−4B_1 = 10^{-4} T。计算：

  1. 拉莫尔频率 ω0\omega_0：
    质子的旋磁比 γp=2.675×108\gamma_p = 2.675 \times 10^8 rad/(s·T)

    ω0=γpB0=2.675×108×1.5=4.01×108 rad/s\omega_0 = \gamma_p B_0 = 2.675 \times 10^8 \times 1.5 = 4.01 \times 10^8 \text{ rad/s}

    对应频率 f0=ω0/(2π)≈63.9f_0 = \omega_0/(2\pi) \approx 63.9 MHz（这正是MRI中常用的射频频率）。

  2. Rabi频率 ω1\omega_1：

    ω1=γpB1=2.675×108×10−4=2.675×104 rad/s\omega_1 = \gamma_p B_1 = 2.675 \times 10^8 \times 10^{-4} = 2.675 \times 10^4 \text{ rad/s}

    对应频率 f1≈4.26f_1 \approx 4.26 kHz。

  3. π\pi-脉冲时间：

    tπ=πω1=π2.675×104≈1.17×10−4 s=117 μst_{\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. 共振时的振荡周期：

    TRabi=2πω1≈235 μsT_{Rabi} = \frac{2\pi}{\omega_1} \approx 235 \text{ }\mu\text{s}

如果射频场频率偏离 ω0\omega_0 仅 0.1%0.1\%（Δ=4×105\Delta = 4 \times 10^5 rad/s），则：

Ω=Δ2+ω12=(4×105)2+(2.675×104)2≈4.01×105 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↑max=ω12Ω2=(2.675×104)2(4.01×105)2≈0.0045P_{\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%0.45\%！这说明MRI中必须极其精确地调谐射频频率到拉莫尔频率。

前置知识：Clebsch-Gordan系数表的使用方法

在深入角动量耦合的数学结构之前，先学会"读表"是极其实用的技能。CG系数是角动量理论中最常用的数值工具之一。

B.1 CG系数的定义与符号

Clebsch-Gordan系数 ⟨j1,m1;j2,m2∣j,m⟩\langle j_1, m_1; j_2, m_2 | j, m \rangle 是未耦合基 ∣j1,m1⟩∣j2,m2⟩|j_1, m_1\rangle |j_2, m_2\rangle 到耦合基 ∣j,m⟩|j, m\rangle 的变换矩阵元。

Condon-Shortley相位约定：

  • 所有CG系数选为实数
  • 约定使得 ⟨j1,j1;j2,j2∣j1+j2,j1+j2⟩=1\langle j_1, j_1; j_2, j_2 | j_1+j_2, j_1+j_2 \rangle = 1
  • 升降算符的矩阵元为正

B.2 查表方法

CG系数表通常按 (j1,j2)(j_1, j_2) 对组织。例如：

两个自旋-1/21/2（j1=1/2,j2=1/2j_1 = 1/2, j_2 = 1/2）：

| 耦合态 ∣j,m⟩|j, m\rangle | 未耦合基展开 |
|—|—|
| ∥1,1⟩\|1, 1\rangle | ∥↑↑⟩\|\uparrow\uparrow\rangle |
| ∥1,0⟩\|1, 0\rangle | 12(∥↑↓⟩+∥↓↑⟩)\frac{1}{\sqrt{2}}(\|\uparrow\downarrow\rangle + \|\downarrow\uparrow\rangle) |
| ∥1,−1⟩\|1, -1\rangle | ∥↓↓⟩\|\downarrow\downarrow\rangle |
| ∥0,0⟩\|0, 0\rangle | 12(∥↑↓⟩−∥↓↑⟩)\frac{1}{\sqrt{2}}(\|\uparrow\downarrow\rangle - \|\downarrow\uparrow\rangle) |

j1=1,j2=1/2j_1 = 1, j_2 = 1/2（例如 L=1L=1 电子与 S=1/2S=1/2 自旋）：

总角动量 j=3/2j = 3/2 或 1/21/2。

∣3/2,3/2⟩=∣1,1⟩∣1/2,1/2⟩|3/2, 3/2\rangle = |1, 1\rangle |1/2, 1/2\rangle

∣3/2,1/2⟩=23∣1,0⟩∣1/2,1/2⟩+13∣1,1⟩∣1/2,−1/2⟩|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⟩=13∣1,0⟩∣1/2,1/2⟩−23∣1,1⟩∣1/2,−1/2⟩|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,mj, m，对不同 m1,m2m_1, m_2 求和

∑m1,m2⟨j1,m1;j2,m2∣j,m⟩⟨j1,m1;j2,m2∣j′,m′⟩=δjj′δmm′\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'}

列正交：固定 m1,m2m_1, m_2，对不同 j,mj, m 求和

∑j,m⟨j1,m1;j2,m2∣j,m⟩⟨j1,m1′;j2,m2′∣j,m⟩=δm1m1′δm2m2′\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系数在交换 j1↔j2j_1 \leftrightarrow j_2 时有：

⟨j2,m2;j1,m1∣j,m⟩=(−1)j1+j2−j⟨j1,m1;j2,m2∣j,m⟩\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→−mm \to -m）性质：

⟨j1,−m1;j2,−m2∣j,−m⟩=(−1)j1+j2−j⟨j1,m1;j2,m2∣j,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 为什么要耦合角动量？

在现实世界中，孤立的一个角动量几乎不存在。原子中的电子同时具有轨道角动量 L^\hat{\mathbf{L}}（来自绕核运动）和自旋角动量 S^\hat{\mathbf{S}}（内禀属性）。两者都与磁场耦合，形成总角动量 J^=L^+S^\hat{\mathbf{J}} = \hat{\mathbf{L}} + \hat{\mathbf{S}}。

更一般地，两个独立子系统各具有角动量 J^1\hat{\mathbf{J}}_1 和 J^2\hat{\mathbf{J}}_2，它们之和 J^=J^1+J^2\hat{\mathbf{J}} = \hat{\mathbf{J}}_1 + \hat{\mathbf{J}}_2 也满足角动量对易关系：

[J^i,J^j]=[J^1i+J^2i,J^1j+J^2j]=iℏϵijkJ^1k+iℏϵijkJ^2k=iℏϵijkJ^k[\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

因为不同子系统的算符对易：[J^1i,J^2j]=0[\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/21/2粒子（比如两个电子）。每个粒子有自旋态 ∣↑⟩|\uparrow\rangle 和 ∣↓⟩|\downarrow\rangle，总系统态空间是4维的张量积空间：

H=H1⊗H2\mathcal{H} = \mathcal{H}_1 \otimes \mathcal{H}_2

基矢为：

∣↑↑⟩=∣↑⟩1⊗∣↑⟩2=(1000)|\uparrow\uparrow\rangle = |\uparrow\rangle_1 \otimes |\uparrow\rangle_2 = \begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix}

∣↑↓⟩=(0100),∣↓↑⟩=(0010)|\uparrow\downarrow\rangle = \begin{pmatrix} 0 \\ 1 \\ 0 \\ 0 \end{pmatrix}, \quad |\downarrow\uparrow\rangle = \begin{pmatrix} 0 \\ 0 \\ 1 \\ 0 \end{pmatrix}

∣↓↓⟩=(0001)|\downarrow\downarrow\rangle = \begin{pmatrix} 0 \\ 0 \\ 0 \\ 1 \end{pmatrix}

总自旋算符 S^=S^1+S^2\hat{\mathbf{S}} = \hat{\mathbf{S}}_1 + \hat{\mathbf{S}}_2。计算 S^2\hat{S}^2 和 S^z\hat{S}_z 在这个基下的矩阵：

S^z=S^1z+S^2z=ℏ2(σz⊗I+I⊗σz)\hat{S}_z = \hat{S}_{1z} + \hat{S}_{2z} = \frac{\hbar}{2}(\sigma_z \otimes I + I \otimes \sigma_z)

S^2=S^12+S^22+2S^1⋅S^2=3ℏ24I+3ℏ24I+2S^1⋅S^2\hat{S}^2 = \hat{\mathbf{S}}_1^2 + \hat{\mathbf{S}}_2^2 + 2\hat{\mathbf{S}}_1\cdot\hat{\mathbf{S}}_2 = \frac{3\hbar^2}{4}I + \frac{3\hbar^2}{4}I + 2\hat{\mathbf{S}}_1\cdot\hat{\mathbf{S}}_2

总自旋点积可以写成：

S^1⋅S^2=S^1zS^2z+12(S^1+S^2−+S^1−S^2+)\hat{\mathbf{S}}_1\cdot\hat{\mathbf{S}}_2 = \hat{S}_{1z}\hat{S}_{2z} + \frac{1}{2}(\hat{S}_{1+}\hat{S}_{2-} + \hat{S}_{1-}\hat{S}_{2+})

通过直接计算（Shankar在第15章有详细推导），我们发现4个态可以重新组合为总自旋本征态：

三重态（triplet，s=1s=1）：

∣1,1⟩=∣↑↑⟩,ms=+1|1, 1\rangle = |\uparrow\uparrow\rangle, \quad m_s = +1

∣1,0⟩=12(∣↑↓⟩+∣↓↑⟩),ms=0|1, 0\rangle = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle + |\downarrow\uparrow\rangle), \quad m_s = 0

∣1,−1⟩=∣↓↓⟩,ms=−1|1, -1\rangle = |\downarrow\downarrow\rangle, \quad m_s = -1

单态（singlet，s=0s=0）：

∣0,0⟩=12(∣↑↓⟩−∣↓↑⟩),ms=0|0, 0\rangle = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle - |\downarrow\uparrow\rangle), \quad m_s = 0

验证：S^2∣1,ms⟩=1(1+1)ℏ2∣1,ms⟩=2ℏ2∣1,ms⟩\hat{S}^2|1, m_s\rangle = 1(1+1)\hbar^2|1, m_s\rangle = 2\hbar^2|1, m_s\rangle，S^2∣0,0⟩=0\hat{S}^2|0, 0\rangle = 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系数：一般情况与显式计算

对于两个任意角动量 j1j_1 和 j2j_2 的耦合，总角动量 jj 的取值范围为：

j=∣j1−j2∣,∣j1−j2∣+1,…,j1+j2j = |j_1 - j_2|, |j_1 - j_2| + 1, \ldots, j_1 + j_2

这是角动量加法规则（triangle rule），可以从经典矢量合成的图像直观理解。

总态数验证：

∑j=∣j1−j2∣j1+j2(2j+1)=(2j1+1)(2j2+1)\sum_{j=|j_1-j_2|}^{j_1+j_2}(2j+1) = (2j_1+1)(2j_2+1)

未耦合基的态 ∣j1,m1⟩⊗∣j2,m2⟩|j_1, m_1\rangle \otimes |j_2, m_2\rangle 到耦合基的 ∣j,m⟩|j, m\rangle 的变换由 Clebsch-Gordan系数 给出：

∣j,m⟩=∑m1,m2⟨j1,m1;j2,m2∣j,m⟩∣j1,m1⟩∣j2,m2⟩|j, m\rangle = \sum_{m_1, m_2} \langle j_1, m_1; j_2, m_2 | j, m\rangle \, |j_1, m_1\rangle |j_2, m_2\rangle

其中 m=m1+m2m = m_1 + m_2（因为 J^z=J^1z+J^2z\hat{J}_z = \hat{J}_{1z} + \hat{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:#ccffcc

CG系数有明确的解析表达式（由Racah公式给出），满足正交归一性：

∑m1,m2⟨j1,m1;j2,m2∣j,m⟩⟨j1,m1;j2,m2∣j′,m′⟩=δjj′δmm′\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'}

∑j,m⟨j1,m1;j2,m2∣j,m⟩⟨j1,m1′;j2,m2′∣j,m⟩=δm1m1′δm2m2′\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}

15.1.4 数值例子：j1=1,j2=1/2j_1 = 1, j_2 = 1/2 的CG系数显式计算

例题：计算 j1=1j_1 = 1（轨道角动量 L=1L=1）和 j2=1/2j_2 = 1/2（自旋 S=1/2S=1/2）耦合的CG系数。这在原子物理中极为常见——pp 轨道电子与自旋耦合形成总角动量 j=3/2j = 3/2 或 1/21/2。

步骤1：确定总角动量取值

j=∣1−1/2∣,∣1−1/2∣+1=1/2,3/2j = |1 - 1/2|, |1 - 1/2| + 1 = 1/2, 3/2

  • j=3/2j = 3/2：m=3/2,1/2,−1/2,−3/2m = 3/2, 1/2, -1/2, -3/2，共4个态
  • j=1/2j = 1/2：m=1/2,−1/2m = 1/2, -1/2，共2个态
  • 总计 4+2=6=(2×1+1)(2×1/2+1)4 + 2 = 6 = (2\times 1 + 1)(2\times 1/2 + 1) ✓

步骤2：用递推法构造最高权态

∣3/2,3/2⟩|3/2, 3/2\rangle 只能来自 m1=1,m2=1/2m_1 = 1, m_2 = 1/2：

∣3/2,3/2⟩=∣1,1⟩∣1/2,1/2⟩|3/2, 3/2\rangle = |1, 1\rangle |1/2, 1/2\rangle

即 ⟨1,1;1/2,1/2∣3/2,3/2⟩=1\langle 1, 1; 1/2, 1/2 | 3/2, 3/2 \rangle = 1。

步骤3：用降算符 J^−=J^1−+J^2−\hat{J}_- = \hat{J}_{1-} + \hat{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⟩\hat{J}_-|3/2, 3/2\rangle = \hbar\sqrt{(3/2+1/2)(3/2-1/2+1)}|3/2, 1/2\rangle = \hbar\sqrt{3}|3/2, 1/2\rangle

左边作用在未耦合态上：

J^−∣1,1⟩∣1/2,1/2⟩=J^1−∣1,1⟩∣1/2,1/2⟩+∣1,1⟩J^2−∣1/2,1/2⟩\hat{J}_-|1,1\rangle|1/2,1/2\rangle = \hat{J}_{1-}|1,1\rangle|1/2,1/2\rangle + |1,1\rangle\hat{J}_{2-}|1/2,1/2\rangle

=ℏ2∣1,0⟩∣1/2,1/2⟩+ℏ∣1,1⟩∣1/2,−1/2⟩= \hbar\sqrt{2}|1,0\rangle|1/2,1/2\rangle + \hbar|1,1\rangle|1/2,-1/2\rangle

因此：

∣3/2,1/2⟩=23∣1,0⟩∣1/2,1/2⟩+13∣1,1⟩∣1/2,−1/2⟩|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,0;1/2,1/2∣3/2,1/2⟩=2/3\langle 1, 0; 1/2, 1/2 | 3/2, 1/2 \rangle = \sqrt{2/3}
  • ⟨1,1;1/2,−1/2∣3/2,1/2⟩=1/3\langle 1, 1; 1/2, -1/2 | 3/2, 1/2 \rangle = \sqrt{1/3}

步骤4：用正交性构造 j=1/2j = 1/2 的态

∣1/2,1/2⟩|1/2, 1/2\rangle 必须与 ∣3/2,1/2⟩|3/2, 1/2\rangle 正交，且是 ∣1,0⟩∣1/2,1/2⟩|1,0\rangle|1/2,1/2\rangle 和 ∣1,1⟩∣1/2,−1/2⟩|1,1\rangle|1/2,-1/2\rangle 的线性组合：

∣1/2,1/2⟩=13∣1,0⟩∣1/2,1/2⟩−23∣1,1⟩∣1/2,−1/2⟩|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

即：

  • ⟨1,0;1/2,1/2∣1/2,1/2⟩=1/3\langle 1, 0; 1/2, 1/2 | 1/2, 1/2 \rangle = \sqrt{1/3}
  • ⟨1,1;1/2,−1/2∣1/2,1/2⟩=−2/3\langle 1, 1; 1/2, -1/2 | 1/2, 1/2 \rangle = -\sqrt{2/3}

步骤5：验证归一性和正交性

∣3/2,1/2⟩|3/2, 1/2\rangle 的系数平方和：2/3+1/3=12/3 + 1/3 = 1 ✓
∣1/2,1/2⟩|1/2, 1/2\rangle 的系数平方和：1/3+2/3=11/3 + 2/3 = 1 ✓
正交性：2/3⋅1/3+1/3⋅(−2/3)=0\sqrt{2/3}\cdot\sqrt{1/3} + \sqrt{1/3}\cdot(-\sqrt{2/3}) = 0 ✓

完整的CG系数表（j1=1,j2=1/2j_1=1, j_2=1/2）：

jjmmm1=1,m2=−1/2m_1=1, m_2=-1/2m1=1,m2=1/2m_1=1, m_2=1/2m1=0,m2=−1/2m_1=0, m_2=-1/2m1=0,m2=1/2m_1=0, m_2=1/2m1=−1,m2=−1/2m_1=-1, m_2=-1/2m1=−1,m2=1/2m_1=-1, m_2=1/2
3/23/2—1————
3/21/21/3\sqrt{1/3}—2/3\sqrt{2/3}———
3/2-1/2—2/3\sqrt{2/3}—1/3\sqrt{1/3}——
3/2-3/2———1/3\sqrt{1/3}—2/3\sqrt{2/3}
1/21/2−2/3-\sqrt{2/3}—1/3\sqrt{1/3}———
1/2-1/2—−1/3-\sqrt{1/3}—2/3\sqrt{2/3}——

这些系数在原子物理中直接决定光谱线的强度和偏振。例如，pp 电子的自旋-轨道耦合将 L=1L=1 能级分裂为 j=3/2j=3/2（四重态）和 j=1/2j=1/2（双重态），这正是碱金属原子光谱双线结构的来源。

15.2 张量算符与Wigner-Eckart定理

15.2.1 什么是张量算符？

在角动量理论中，物理量按其在旋转下的变换性质分类。一个不可约张量算符 T^(k)\hat{T}^{(k)} 是秩为 kk 的 (2k+1)(2k+1) 个算符的集合 {T^q(k):q=−k,…,k}\{\hat{T}^{(k)}_q : q = -k, \ldots, k\}，在旋转下按照角动量 kk 的表示变换：

D^(R)T^q(k)D^†(R)=∑q′Dq′q(k)(R)T^q′(k)\hat{D}(R)\hat{T}^{(k)}_q\hat{D}^\dagger(R) = \sum_{q'} D^{(k)}_{q'q}(R)\hat{T}^{(k)}_{q'}

例子：

  • 标量（k=0k=0）：T^0(0)\hat{T}^{(0)}_0，旋转不变。如 J^2\hat{\mathbf{J}}^2。
  • 矢量（k=1k=1）：三个分量 T^−1(1),T^0(1),T^+1(1)\hat{T}^{(1)}_{-1}, \hat{T}^{(1)}_0, \hat{T}^{(1)}_{+1}。如位置 r^\hat{\mathbf{r}} 或角动量 J^\hat{\mathbf{J}} 的球分量。
  • 球谐函数 Ylm(θ,ϕ)Y_{lm}(\theta,\phi)：本身就是秩 ll 的张量。

张量算符与角动量的对易关系（定义性质）：

[J^±,T^q(k)]=ℏk(k+1)−q(q±1)T^q±1(k)[\hat{J}_\pm, \hat{T}^{(k)}_q] = \hbar\sqrt{k(k+1)-q(q\pm1)} \hat{T}^{(k)}_{q\pm 1}

[J^z,T^q(k)]=ℏqT^q(k)[\hat{J}_z, \hat{T}^{(k)}_q] = \hbar q \hat{T}^{(k)}_q

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)\hat{T}^{(k)}_q 在角动量本征态之间的矩阵元可以分解为：

⟨j′,m′∣T^q(k)∣j,m⟩=⟨j,m;k,q∣j′,m′⟩⟨j′∣∣T^(k)∣∣j⟩2j′+1\langle j', m' | \hat{T}^{(k)}_q | j, m \rangle = \langle j, m; k, q | j', m' \rangle \frac{\langle j' || \hat{T}^{(k)} || j \rangle}{\sqrt{2j'+1}}

其中：

  • ⟨j,m;k,q∣j′,m′⟩\langle j, m; k, q | j', m' \rangle 是 Clebsch-Gordan系数
  • ⟨j′∣∣T^(k)∣∣j⟩\langle j' || \hat{T}^{(k)} || j \rangle 是 约化矩阵元（reduced matrix element），与 m,m′,qm, m', q 无关

物理意义：

  1. 选择定则：如果CG系数为零，矩阵元就为零。CG系数非零的条件是 m′=m+qm' = m + q 且 ∣j−k∣≤j′≤j+k|j-k| \leq j' \leq j+k。这就是原子光谱中的选择定则（selection rules）。
  2. 比例关系：不同 (m,q,m′)(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^\hat{\mathbf{d}} = -e\hat{\mathbf{r}} 是**秩1（矢量）**张量算符。

Wigner-Eckart定理告诉我们，跃迁矩阵元 ⟨n′,l′,m′∣r^q∣n,l,m⟩\langle n', l', m' | \hat{r}_q | n, l, m \rangle 非零的条件：

  • mm 选择定则：m′=m+qm' = m + q，其中 q=−1,0,+1q = -1, 0, +1（对应 x±iyx\pm iy 和 zz）
  • ll 选择定则：∣l−1∣≤l′≤l+1|l-1| \leq l' \leq l+1，即 l′=l−1,l,l+1l' = l-1, l, l+1。但由于宇称，实际上 l′=l±1l' = l \pm 1
  • 自旋：Δs=0\Delta s = 0（电偶极不改变自旋）

这些选择定则解释了为什么原子光谱是线状谱而非连续谱，以及为什么只有某些跃迁（allowed transitions）强度高。

数值例子：CG系数在原子光谱中的应用

例题：钠原子（Na）的 3p→3s3p \to 3s 跃迁（D线，589 nm）。3p3p 态有 l=1l=1，3s3s 态有 l=0l=0。考虑自旋-轨道耦合后：

  • 3p3/23p_{3/2}：j=3/2j = 3/2，四重态（mj=±3/2,±1/2m_j = \pm 3/2, \pm 1/2）
  • 3p1/23p_{1/2}：j=1/2j = 1/2，双重态（mj=±1/2m_j = \pm 1/2）

问题：从 3p3/23p_{3/2} 到 3s1/23s_{1/2}（j=1/2j=1/2）的跃迁中，计算CG系数决定的相对强度。

电偶极算符是秩1张量，k=1k=1。j=3/2→j′=1/2j=3/2 \to j'=1/2，k=1k=1。

CG系数 ⟨3/2,m;1,q∣1/2,m′⟩\langle 3/2, m; 1, q | 1/2, m' \rangle 非零的条件：

  • m′=m+qm' = m + q
  • ∣3/2−1∣≤1/2≤3/2+1|3/2 - 1| \leq 1/2 \leq 3/2 + 1，即 1/2≤1/2≤5/21/2 \leq 1/2 \leq 5/2 ✓

具体系数（查表或递推）：

对于 m′=1/2m' = 1/2：

  • m=3/2,q=−1m = 3/2, q = -1：⟨3/2,3/2;1,−1∣1/2,1/2⟩=−1/2\langle 3/2, 3/2; 1, -1 | 1/2, 1/2 \rangle = -\sqrt{1/2}
  • m=1/2,q=0m = 1/2, q = 0：⟨3/2,1/2;1,0∣1/2,1/2⟩=1/6\langle 3/2, 1/2; 1, 0 | 1/2, 1/2 \rangle = \sqrt{1/6}
  • m=−1/2,q=+1m = -1/2, q = +1：⟨3/2,−1/2;1,1∣1/2,1/2⟩=−1/3\langle 3/2, -1/2; 1, 1 | 1/2, 1/2 \rangle = -\sqrt{1/3}

相对跃迁强度正比于CG系数的平方：

| 初态 mm | 终态 m′m' | ∣CG∣2|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 |

对所有初态求平均（假设初态各 mjm_j 等概率占据）：

∣CG∣2‾=14(12+16+13+12)=14×3+1+2+36=924=38\overline{|CG|^2} = \frac{1}{4}\left(\frac{1}{2} + \frac{1}{6} + \frac{1}{3} + \frac{1}{2}\right) = \frac{1}{4} \times \frac{3+1+2+3}{6} = \frac{9}{24} = \frac{3}{8}

类似地，对于 3p1/2→3s1/23p_{1/2} \to 3s_{1/2}：

∣CG∣2‾=12(13+13)=13\overline{|CG|^2} = \frac{1}{2}\left(\frac{1}{3} + \frac{1}{3}\right) = \frac{1}{3}

两者强度比：

I(3p3/2→3s1/2)I(3p1/2→3s1/2)=(2j+1)∣CG∣2‾3/2(2j+1)∣CG∣2‾1/2=4×3/82×1/3=3/22/3=94=2.25\frac{I(3p_{3/2} \to 3s_{1/2})}{I(3p_{1/2} \to 3s_{1/2})} = \frac{(2j+1)\overline{|CG|^2}_{3/2}}{(2j+1)\overline{|CG|^2}_{1/2}} = \frac{4 \times 3/8}{2 \times 1/3} = \frac{3/2}{2/3} = \frac{9}{4} = 2.25

实验上，钠D线的 D2D_2（589.0 nm，3p3/2→3s1/23p_{3/2} \to 3s_{1/2}）与 D1D_1（589.6 nm，3p1/2→3s1/23p_{1/2} \to 3s_{1/2}）强度比约为 2:12: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("\\"\hat{n}, \theta\\"") = \exp("\\"-i\theta \hat{n}\cdot\sigma/2\\"")$"]
    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/21/2粒子的行为，也是理解所有角动量耦合、选择定则和量子调控的基础。从核磁共振到量子计算，从原子光谱到粒子物理，SU(2)的语言无处不在。

练习与思考

1. 泡利矩阵的恒等式

证明对于任意两个矢量 a\mathbf{a} 和 b\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})

并由此推导：如果 a\mathbf{a} 和 b\mathbf{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\sigma_i\sigma_j = \delta_{ij}I + i\epsilon_{ijk}\sigma_k。

2. SU(2)到SO(3)的映射验证

考虑 U=exp(−iθσy/2)=cos(θ/2)I−isin(θ/2)σyU = \exp(-i\theta \sigma_y/2) = \cos(\theta/2)I - i\sin(\theta/2)\sigma_y。证明这个 2×22\times 2 幺正矩阵作用于泡利矩阵的方式：

UσxU†=σxcosθ+σzsinθU\sigma_x U^\dagger = \sigma_x\cos\theta + \sigma_z\sin\theta

UσzU†=−σxsinθ+σzcosθU\sigma_z U^\dagger = -\sigma_x\sin\theta + \sigma_z\cos\theta

这与SO(3)中绕 yy 轴旋转 θ\theta 的变换一致。为什么这说明了SU(2)到SO(3)的2对1同态？

3. Rabi振荡的量子门应用

在量子计算中，自旋-1/21/2（或等效的两能级系统）是最简单的量子比特。证明：

  • 一个 π\pi-脉冲（ω1t=π\omega_1 t = \pi，共振）实现的是NOT门（∣↓⟩↔∣↑⟩|\downarrow\rangle \leftrightarrow |\uparrow\rangle）
  • 一个 π/2\pi/2-脉冲实现的是Hadamard门的类似操作（将基态变为等权重叠加态）
  • 如果在失谐条件下（Δ≫ω1\Delta \gg \omega_1），系统近似只经历相位演化，可以实现相位门

设计一个操作序列来实现从 ∣↓⟩|\downarrow\rangle 到 12(∣↓⟩+i∣↑⟩)\frac{1}{\sqrt{2}}(|\downarrow\rangle + i|\uparrow\rangle) 的变换。

4. CG系数显式计算

用递推法推导 j1=1,j2=1j_1 = 1, j_2 = 1 耦合的CG系数。总角动量 j=2,1,0j = 2, 1, 0。验证：

  • ∣2,2⟩=∣1,1⟩∣1,1⟩|2, 2\rangle = |1,1\rangle|1,1\rangle
  • ∣2,1⟩=12(∣1,0⟩∣1,1⟩+∣1,1⟩∣1,0⟩)|2, 1\rangle = \frac{1}{\sqrt{2}}(|1,0\rangle|1,1\rangle + |1,1\rangle|1,0\rangle)
  • ∣1,1⟩=12(∣1,0⟩∣1,1⟩−∣1,1⟩∣1,0⟩)|1, 1\rangle = \frac{1}{\sqrt{2}}(|1,0\rangle|1,1\rangle - |1,1\rangle|1,0\rangle)
  • ∣0,0⟩=13(∣1,1⟩∣1,−1⟩−∣1,0⟩∣1,0⟩+∣1,−1⟩∣1,1⟩)|0, 0\rangle = \frac{1}{\sqrt{3}}(|1,1\rangle|1,-1\rangle - |1,0\rangle|1,0\rangle + |1,-1\rangle|1,1\rangle)

5. Wigner-Eckart定理与选择定则

证明：对于电偶极跃迁（k=1k=1），如果初态和终态的宇称相同，则矩阵元必为零。这解释了为什么 Δl=0\Delta l = 0 的跃迁是禁戒的（电偶极算符是奇宇称，所以 ll 必须改变奇数个单位）。