第11-12章 对称性与角动量:守恒律的量子版本

📑 目录

第11-12章 对称性与角动量:守恒律的量子版本

故事引入:镜像世界的密码

2089年,"对称号"空间站发现了一个诡异的现象。宇航员陈默在进行量子态传输实验时,每次将物体从舱A传送到舱B,镜像反射后的版本总是失败。经过三个月排查,他在日志中写道:"宇宙似乎 handedness——它有左右之分。"更深层的问题随之浮现:为什么空间平移不变性保证动量守恒,旋转不变性保证角动量守恒?这些经典力学中的"巧合"在量子世界中有着怎样的数学根基?当陈默将实验数据发给地球的理论组时,一位老教授回了一封只有一行公式的邮件:[H,Lz]=0dLz/dt=0[H, L_z] = 0 \Rightarrow d\langle L_z \rangle/dt = 0。"去读懂它,"教授说,"这就是宇宙对称性的全部密码。"而这,正是Shankar在第11-12章要带我们穿越的旅程——从平移和旋转的数学定义,到角动量算符的代数结构,再到氢原子解中自然涌现的球谐函数。对称性不是装饰,它是量子力学的骨架。


11.0 前置知识:群论基础

对称性的数学语言是群论。Shankar在第11章用到了群论的基本概念,但通常以物理学家的"实用主义"方式呈现——先使用后严谨定义。让我们在这里补充必要的群论基础,作为阅读本章的数学钥匙。

11.0.1 群的定义

一个 (G,)(G, \cdot) 是一个集合 GG 配上一个二元运算 \cdot,满足四条公理:

  1. 封闭性a,bG,abG\forall a, b \in G, \quad a \cdot b \in G
  2. 结合律a,b,cG,(ab)c=a(bc)\forall a, b, c \in G, \quad (a \cdot b) \cdot c = a \cdot (b \cdot c)
  3. 单位元eG,aG,ea=ae=a\exists e \in G, \quad \forall a \in G, \quad e \cdot a = a \cdot e = a
  4. 逆元aG,a1G,aa1=a1a=e\forall a \in G, \quad \exists a^{-1} \in G, \quad a \cdot a^{-1} = a^{-1} \cdot a = e

例子

  • 三维旋转群 SO(3)SO(3):所有三维空间中的旋转矩阵,运算为矩阵乘法
  • 平移群:所有空间平移变换,运算为变换的复合
  • 宇称群:{I,Π}\{I, \Pi\},其中 Π2=I\Pi^2 = I,这是一个只有两个元素的群

11.0.2 子群与生成元

一个子群 HGH \subset GGG 的子集,自身也构成群。

生成元:对于连续群(Lie群),无穷小变换的集合由生成元决定。例如 SO(3)SO(3) 的三个生成元就是角动量分量 Lx,Ly,LzL_x, L_y, L_z。任何有限旋转都可以通过对生成元的指数映射得到:

R(n,θ)=exp(iθnL)R(\mathbf{n}, \theta) = \exp\left(-\frac{i\theta \mathbf{n} \cdot \mathbf{L}}{\hbar}\right)

11.0.3 群的表示

群的表示是将群的每个元素映射为线性空间上的算符(或矩阵),并保持群乘法结构:

D(g1g2)=D(g1)D(g2)D(g_1 \cdot g_2) = D(g_1) D(g_2)

表示的维度是目标线性空间的维度。SO(3)SO(3) 的不可约表示正是角动量代数提供的 (2l+1)(2l+1) 维空间——这就是为什么 l=0l=0(1维)、l=1l=1(3维)、l=2l=2(5维)等在原子光谱中自然出现。

graph TD
    A[群论基础] --> B["群G: 封闭性+结合律+单位元+逆元"]
    B --> C["连续群/Lie群"]
    C --> D["生成元: 无穷小变换"]
    D --> E["SO(3)生成元: L_x, L_y, L_z"]
    B --> F[群的表示]
    F --> G[将群元映射为算符]
    G --> H[不可约表示]
    H --> I["SO(3)不可约表示维度 = 2l+1"]
    I --> J["l=0: 1维
s轨道"] I --> K["l=1: 3维
p轨道"] I --> L["l=2: 5维
d轨道"] style B fill:#e3f2fd style I fill:#fff8e1

11.1 对称变换与守恒律:Noether定理的量子版本

11.1.1 主动变换 vs 被动变换

在深入数学之前,必须澄清一个概念陷阱。当我们说"把系统向右平移1米"时,有两种截然不同的含义:

  • 主动变换:保持坐标系不动,将物理系统(波函数、粒子)向右移动
  • 被动变换:保持系统不动,将坐标系原点向左移动1米

两种视角描述的是同一个物理状态,只是参考系不同。在量子力学中,主动变换由幺正算符 UU 实现:

ψ=Uψ|\psi'\rangle = U|\psi\rangle

11.1.2 连续对称性与守恒量

考虑一个无穷小变换 U(ϵ)=IiϵG/+O(ϵ2)U(\epsilon) = I - i\epsilon G/\hbar + O(\epsilon^2),其中 GG 是厄米算符(称为变换的生成元)。如果Hamiltonian在该变换下不变:

UHU=H[G,H]=0UHU^\dagger = H \quad \Rightarrow \quad [G, H] = 0

则由Heisenberg运动方程:

dGdt=i[H,G]+Gt=0\frac{dG}{dt} = \frac{i}{\hbar}[H, G] + \frac{\partial G}{\partial t} = 0

GG 是守恒量! 这就是量子版本的Noether定理:每一个连续对称性对应一个守恒量。

graph TD
    A[连续对称变换] --> B[无穷小变换]
    B --> C["U(ε) = I - iεG/ℏ"]
    C --> D["G: 厄米生成元"]
    D --> E{"[G,H] = 0?"}
    E -->|是| F[对称性存在]
    F --> G["dG/dt = 0"]
    G --> H[G是守恒量]
    E -->|否| I[对称性破缺]
    H --> J["平移→动量守恒"]
    H --> K["旋转→角动量守恒"]
    H --> L["时间平移→能量守恒"]
    style F fill:#e8f5e9
    style I fill:#ffebee

11.1.3 平移不变性 → 动量守恒

无穷小平移 xx+ϵx \to x + \epsilon 作用于波函数:

ψ(x)ψ(xϵ)ψ(x)ϵdψdx=(1iϵP)ψ(x)\psi(x) \to \psi(x - \epsilon) \approx \psi(x) - \epsilon \frac{d\psi}{dx} = \left(1 - \frac{i\epsilon P}{\hbar}\right)\psi(x)

生成元正是动量算符 P=iddxP = -i\hbar \frac{d}{dx}!如果 [P,H]=0[P, H] = 0,动量守恒。

推广到三维:P=(Px,Py,Pz)\mathbf{P} = (P_x, P_y, P_z),平移不变性意味着空间均匀性,即"宇宙没有特殊位置"。

11.1.4 时间平移不变性 → 能量守恒

时间平移 tt+τt \to t + \tau 的生成元是Hamiltonian本身。如果 HH 不显含时间(H/t=0\partial H/\partial t = 0),则:

[H, H] = 0 \quad \text{(平凡成立)}

能量自然守恒。如果 HH 显含时间(例如外加交变电场),则能量不守恒——系统与外界交换能量。

11.1.5 空间反演与宇称

宇称变换(Parity)是一个分立对称性(不是连续的):

Π:rr\Pi: \mathbf{r} \to -\mathbf{r}

由于 Π2=I\Pi^2 = I,宇称算符的本征值只能是 ±1\pm 1

Πψ=±ψ\Pi|\psi\rangle = \pm |\psi\rangle

  • η=+1\eta = +1偶宇称(even parity),如 ψ(x)=+ψ(x)\psi(-x) = +\psi(x)
  • η=1\eta = -1奇宇称(odd parity),如 ψ(x)=ψ(x)\psi(-x) = -\psi(x)

宇称守恒曾是"神圣"的守恒律,直到1956年李政道和杨振宁提出弱相互作用中宇称不守恒,吴健雄的 60^{60}Co实验证实了这一惊人预言。在强相互作用和电磁相互作用中,宇称仍然守恒;在弱相互作用中,大自然显示出了"手性偏好"。

graph LR
    subgraph 连续对称
        A1[平移]
        B1[旋转]
        C1[时间平移]
        D1["生成元: P, L, H"]
    end
    E[Noether定理]
    subgraph 分立对称
        A2["宇称 Π"]
        B2[时间反演 T]
        C2[电荷共轭 C]
        D2["本征值: ±1"]
    end
    A1 --> E
    B1 --> E
    C1 --> E
    A2 --> F[分立守恒律]
    B2 --> F
    E --> G["动量/角动量/能量守恒"]
    F --> H["η = ±1"]
    style E fill:#e3f2fd
    style F fill:#fff3e0

11.2 旋转对称性与角动量

11.2.1 二维旋转:单参数群

在二维平面中,绕原点旋转角度 ϕ\phi 的变换矩阵为:

R(ϕ)=(cosϕsinϕsinϕcosϕ)R(\phi) = \begin{pmatrix} \cos\phi & -\sin\phi \\ \sin\phi & \cos\phi \end{pmatrix}

无穷小旋转 ϕ=ϵ1\phi = \epsilon \ll 1

R(ϵ)(1ϵϵ1)=Iϵ(0110)R(\epsilon) \approx \begin{pmatrix} 1 & -\epsilon \\ \epsilon & 1 \end{pmatrix} = I - \epsilon \begin{pmatrix} 0 & 1 \\ -1 & 0 \end{pmatrix}

在量子力学中,绕z轴的无穷小旋转算符为:

U(Rz(ϵ))=IiϵLzU(R_z(\epsilon)) = I - \frac{i\epsilon L_z}{\hbar}

其中 Lz=XPyYPx=iϕL_z = XP_y - YP_x = -i\hbar \frac{\partial}{\partial\phi}(在球坐标中)。

11.2.2 三维旋转:非对易群

三维旋转有三个独立的生成元LxL_xLyL_yLzL_z。它们的对易关系构成了角动量代数的核心:

\boxed{[L_i, L_j] = i\hbar \varepsilon_{ijk} L_k}

用分量写出来:

  • [Lx,Ly]=iLz[L_x, L_y] = i\hbar L_z
  • [Ly,Lz]=iLx[L_y, L_z] = i\hbar L_x
  • [Lz,Lx]=iLy[L_z, L_x] = i\hbar L_y

这个代数结构的非凡之处在于:它不依赖于具体的坐标表示。只要三个算符满足上述对易关系,它们就是"角动量"——无论它们是轨道角动量 L=r×pL = r \times p,还是自旋角动量 SS,或者两者之和 J=L+SJ = L + S

graph TD
    A["三维旋转群 SO(3)"] --> B["三个生成元 L_x, L_y, L_z"]
    B --> C[对易关系]
    C --> D["[L_i, L_j] = iℏ ε_ijk L_k"]
    D --> E["李代数 so(3)"]
    E --> F[不可约表示]
    F --> G["表示空间: 2l+1维"]
    G --> H["l = 0, 1, 2, ..."]
    H --> I["s轨道, p轨道, d轨道..."]
    D --> J[与具体实现无关]
    J --> K[轨道角动量 L]
    J --> L[自旋角动量 S]
    J --> M["总角动量 J = L+S"]
    style D fill:#ffebee
    style E fill:#e3f2fd

12.0 前置知识:三维旋转的几何

在进入角动量的量子代数之前,让我们先在经典几何中建立三维旋转的直觉。量子力学中的旋转算符继承经典旋转的结构,但作用在Hilbert空间而非欧几里得空间。

12.0.1 Euler角与旋转分解

任何三维旋转都可以分解为三个绕不同轴的旋转(Euler角):

R(α,β,γ)=Rz(α)Ry(β)Rz(γ)R(\alpha, \beta, \gamma) = R_z(\alpha) R_y(\beta) R_z(\gamma)

其中 α,β,γ\alpha, \beta, \gamma 分别是绕z轴、y轴、z轴的旋转角。这种分解不是唯一的,但它展示了三维旋转的三个自由度。

关键的几何事实:三维旋转不可对易。先绕x轴转90°再绕y轴转90°,与先绕y轴转90°再绕x轴转90°,结果完全不同。这正是角动量对易关系 [Lx,Ly]=iLz0[L_x, L_y] = i\hbar L_z \neq 0 的经典起源。

12.0.2 旋转的矩阵表示

绕z轴旋转角度 ϕ\phi 的矩阵:

Rz(ϕ)=(cosϕsinϕ0sinϕcosϕ0001)R_z(\phi) = \begin{pmatrix} \cos\phi & -\sin\phi & 0 \\ \sin\phi & \cos\phi & 0 \\ 0 & 0 & 1 \end{pmatrix}

绕y轴:

Ry(ϕ)=(cosϕ0sinϕ010sinϕ0cosϕ)R_y(\phi) = \begin{pmatrix} \cos\phi & 0 & \sin\phi \\ 0 & 1 & 0 \\ -\sin\phi & 0 & \cos\phi \end{pmatrix}

验证 Rx(ϵ)Ry(ϵ)Ry(ϵ)Rx(ϵ)R_x(\epsilon) R_y(\epsilon) \neq R_y(\epsilon) R_x(\epsilon)O(ϵ2)O(\epsilon^2) 阶,即可得到对易子的几何解释。

12.0.3 从经典到量子:几何→代数

经典旋转作用于三维向量 r=(x,y,z)\mathbf{r} = (x, y, z)。量子旋转作用于波函数 ψ(r)\psi(\mathbf{r}),或更抽象地作用于态向量 ψ|\psi\rangle。但两者的数学结构(群和对易关系)是同一的。这就是为什么Shankar可以从纯代数对易关系出发,推导出角动量的全部量子化性质——几何直觉给了我们信心,但代数推导给出了严格性。


12.1 LzL_z 的本征值问题

12.1.1 环上的量子力学

考虑一个粒子被限制在半径为 rr 的圆环上运动(例如在平面极坐标中,固定 rr)。Hamiltonian简化为:

H=Lz22I,I=mr2H = \frac{L_z^2}{2I}, \quad I = mr^2

LzL_z 的本征方程:

Lzm=mmL_z |m\rangle = m\hbar |m\rangle

在坐标表象中:

iddϕψm(ϕ)=mψm(ϕ)-i\hbar \frac{d}{d\phi} \psi_m(\phi) = m\hbar \psi_m(\phi)

解为:

ψm(ϕ)=12πeimϕ,m=0,±1,±2,\psi_m(\phi) = \frac{1}{\sqrt{2\pi}} e^{im\phi}, \quad m = 0, \pm 1, \pm 2, \ldots

周期性边界条件 ψ(ϕ+2π)=ψ(ϕ)\psi(\phi + 2\pi) = \psi(\phi) 要求 mm 必须是整数!这是角动量量子化的来源——不是人为假设,而是空间拓扑结构的直接后果。如果粒子生活在一个"有缝"的环上(ϕ\phi 范围 [0,2π)[0, 2\pi) 但端点不连接),mm 可以是任意实数,因为边界条件改变了。

12.1.2 从离散到连续:平面波的类比

这个结构与无限深势阱中的动量本征态完全类似。在无限直线上,PP 的本征值 pp 连续;在周期环上,LzL_z 的本征值 mm\hbar 离散。拓扑约束决定了谱的性质。

graph LR
    subgraph 无限直线
        A1["X ∈ (-∞, +∞)"]
        B1["P = -iℏ d/dx"]
        C1["本征值: p连续"]
        D1["ψ_p(x) = e^{ipx/ℏ}"]
    end
    subgraph 周期圆环
        A2[φ ∈ ["0, 2π)"]
        B2["L_z = -iℏ d/dφ"]
        C2["本征值: mℏ, m∈ℤ"]
        D2["ψ_m(φ) = e^{imφ}/√(2π)"]
    end
    E[拓扑决定谱结构]
    F[周期性边界条件]
    A1 --> E
    A2 --> E
    A2 --> F
    F --> C2
    style E fill:#e3f2fd

12.2 角动量代数与升降算符

12.2.1 构造阶梯:L+L_+LL_-

面对对易关系 [Lz,L±]=±L±[L_z, L_\pm] = \pm \hbar L_\pm,其中 L±=Lx±iLyL_\pm = L_x \pm i L_y,Shankar展示了角动量代数最优美的结构之一——阶梯算符(Ladder operators)。

l,m|l, m\rangleL2L^2LzL_z 的共同本征态:

L2l,m=2l(l+1)l,mL^2 |l, m\rangle = \hbar^2 l(l+1) |l, m\rangle

Lzl,m=ml,mL_z |l, m\rangle = m\hbar |l, m\rangle

L+L_+ 作用于 l,m|l, m\rangle 会提升 mm 值:

Lz(L+l,m)=(L+Lz+L+)l,m=(m+1)(L+l,m)L_z (L_+ |l, m\rangle) = (L_+ L_z + \hbar L_+) |l, m\rangle = (m+1)\hbar (L_+ |l, m\rangle)

所以 L+l,ml,m+1L_+ |l, m\rangle \propto |l, m+1\rangle。类似地,LL_- 降低 mm 值。

12.2.2 本征值的边界与量子化

关键问题:这个梯子能无限延伸吗?由于 L2Lz2=Lx2+Ly20L^2 - L_z^2 = L_x^2 + L_y^2 \geq 0(两个厄米算符平方和),我们有:

2l(l+1)m220l(l+1)ml(l+1)\hbar^2 l(l+1) - m^2 \hbar^2 \geq 0 \quad \Rightarrow \quad -\sqrt{l(l+1)} \leq m \leq \sqrt{l(l+1)}

L+L_+ 每次增加 mm 1个单位,所以 mm 必须有上界。设最大值为 mmaxm_{\max},则:

L+l,mmax=0L_+ |l, m_{\max}\rangle = 0

利用 LL+=L2Lz2LzL_- L_+ = L^2 - L_z^2 - \hbar L_z

l,mmaxLL+l,mmax=2[l(l+1)mmax(mmax+1)]=0\langle l, m_{\max} | L_- L_+ | l, m_{\max} \rangle = \hbar^2 [l(l+1) - m_{\max}(m_{\max}+1)] = 0

解得 mmax=lm_{\max} = l。类似地,mmin=lm_{\min} = -l。由于 mm 每次变化1,从 l-lll 共有 2l+12l+1 个值,所以 2l+12l+1 必须是整数 \Rightarrow ll 必须是整数或半整数

\boxed{l = 0, \frac{1}{2}, 1, \frac{3}{2}, 2, \ldots, \quad m = -l, -l+1, \ldots, l-1, l}

12.2.3 代数推导的普适性

这个推导完全没有用到 L=r×pL = r \times p 的具体形式!它只依赖对易关系 [Li,Lj]=iεijkLk[L_i, L_j] = i\hbar \varepsilon_{ijk} L_k。这意味着:

  • 轨道角动量ll 必须是整数(因为波函数必须在 2π2\pi 旋转下单值),对应 l=0,1,2,l = 0, 1, 2, \ldots
  • 自旋角动量ss 可以是半整数(如 s=1/2s = 1/2),因为自旋没有经典对应,不需要单值性条件
graph TD
    A["角动量代数
[L_i, L_j] = iℏ ε_ijk L_k"] --> B["定义升降算符
L± = L_x ± iL_y"] B --> C["L±改变m ±1"] C --> D{"L² - L_z² ≥ 0"} D --> E[m有上下界] E --> F["m_max = l
m_min = -l"] F --> G["l必须是
整数或半整数"] G --> H["2l+1个态
= 表示维度"] H --> I["轨道: l=0,1,2,..."] H --> J["自旋: s=1/2,3/2,..."] I --> K["s, p, d, f轨道"] J --> L["电子自旋
光子自旋"] style G fill:#fff8e1 style B fill:#e8f5e9

12.2.4 数值例子:升降算符的矩阵元计算

让我们计算 l=1l=1 系统中升降算符的具体矩阵元。l=1l=1m=1,0,+1m = -1, 0, +1,基矢为 1,1,1,0,1,+1|1, -1\rangle, |1, 0\rangle, |1, +1\rangle

归一化公式

L±l,m=l(l+1)m(m±1)l,m±1L_\pm |l, m\rangle = \hbar \sqrt{l(l+1) - m(m\pm 1)} |l, m\pm 1\rangle

计算各矩阵元

  1. L+1,1=2(1)(0)1,0=21,0L_+ |1, -1\rangle = \hbar \sqrt{2 - (-1)(0)} |1, 0\rangle = \hbar\sqrt{2} |1, 0\rangle

  2. L+1,0=201,+1=21,+1L_+ |1, 0\rangle = \hbar \sqrt{2 - 0} |1, +1\rangle = \hbar\sqrt{2} |1, +1\rangle

  3. L+1,+1=221,+2=0L_+ |1, +1\rangle = \hbar \sqrt{2 - 2} |1, +2\rangle = 0(上界!)

  4. L1,+1=2101,0=21,0L_- |1, +1\rangle = \hbar \sqrt{2 - 1 \cdot 0} |1, 0\rangle = \hbar\sqrt{2} |1, 0\rangle

  5. L1,0=201,1=21,1L_- |1, 0\rangle = \hbar \sqrt{2 - 0} |1, -1\rangle = \hbar\sqrt{2} |1, -1\rangle

  6. L1,1=0L_- |1, -1\rangle = 0(下界!)

矩阵形式(以 1,+1,1,0,1,1|1,+1\rangle, |1,0\rangle, |1,-1\rangle 为基):

L+=(020002000),L=(000200020)L_+ = \hbar \begin{pmatrix} 0 & \sqrt{2} & 0 \\ 0 & 0 & \sqrt{2} \\ 0 & 0 & 0 \end{pmatrix}, \quad L_- = \hbar \begin{pmatrix} 0 & 0 & 0 \\ \sqrt{2} & 0 & 0 \\ 0 & \sqrt{2} & 0 \end{pmatrix}

验证 LxL_xLyL_y

Lx=L++L2=2(010101010)L_x = \frac{L_+ + L_-}{2} = \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix}

Ly=L+L2i=2(0i0i0i0i0)L_y = \frac{L_+ - L_-}{2i} = \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & -i & 0 \\ i & 0 & -i \\ 0 & i & 0 \end{pmatrix}

验证对易关系 [Lx,Ly]=iLz[L_x, L_y] = i\hbar L_z

其中 Lz=(100000001)L_z = \hbar \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & -1 \end{pmatrix}

直接矩阵乘法可验证对易关系成立(留给读者作为练习)。这个例子展示了角动量代数的具体实现——从抽象对易关系到可触摸的3×3矩阵。

12.2.5 数值例子:Clebsch-Gordan系数的数值计算

考虑两个自旋 s1=s2=1/2s_1 = s_2 = 1/2 的耦合(如两个电子)。单粒子基矢:=1/2,+1/2|\uparrow\rangle = |1/2, +1/2\rangle=1/2,1/2|\downarrow\rangle = |1/2, -1/2\rangle

总自旋 S=S1+S2S = S_1 + S_2 的可能取值:s=1s = 1(三重态)或 s=0s = 0(单态)。

s=1s = 1(三重态,对称)

1,+1=|1, +1\rangle = |\uparrow\uparrow\rangle

1,0=12(+)|1, 0\rangle = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle + |\downarrow\uparrow\rangle)

1,1=|1, -1\rangle = |\downarrow\downarrow\rangle

s=0s = 0(单态,反对称)

0,0=12()|0, 0\rangle = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle - |\downarrow\uparrow\rangle)

CG系数表

| m1m_1 | m2m_2 | 1,+1|1, +1\rangle | 1,0|1, 0\rangle | 1,1|1, -1\rangle | 0,0|0, 0\rangle |
|-------|-------|---------------|---------------|----------------|---------------|
| +1/2 | +1/2 | 1 | 0 | 0 | 0 |
| +1/2 | -1/2 | 0 | 1/21/\sqrt{2} | 0 | 1/21/\sqrt{2} |
| -1/2 | +1/2 | 0 | 1/21/\sqrt{2} | 0 | 1/2-1/\sqrt{2} |
| -1/2 | -1/2 | 0 | 0 | 1 | 0 |

这些系数决定了氢原子基态的超精细结构(质子与电子自旋耦合)和氦原子的自旋排列。两个电子处于单态时空间波函数对称,可以靠近;处于三重态时空间波函数反对称,彼此"回避"。


12.3 球谐函数与三维旋转

12.3.1 L2L^2LzL_z 的共同本征函数

在球坐标中:

L2=2[1sinθθ(sinθθ)+1sin2θ2ϕ2]L^2 = -\hbar^2 \left[\frac{1}{\sin\theta} \frac{\partial}{\partial\theta}\left(\sin\theta \frac{\partial}{\partial\theta}\right) + \frac{1}{\sin^2\theta} \frac{\partial^2}{\partial\phi^2}\right]

Lz=iϕL_z = -i\hbar \frac{\partial}{\partial\phi}

LzL_z 的本征函数自动包含 eimϕe^{im\phi} 因子。代入 L2L^2 的本征方程,令 x=cosθx = \cos\theta,得到连带Legendre方程

ddx[(1x2)dPdx]+[l(l+1)m21x2]P=0\frac{d}{dx}\left[(1-x^2)\frac{dP}{dx}\right] + \left[l(l+1) - \frac{m^2}{1-x^2}\right]P = 0

其解为连带Legendre函数 Plm(cosθ)P_l^m(\cos\theta)

归一化的完整本征函数是球谐函数

\boxed{Y_l^m(\theta, \phi) = (-1)^m \sqrt{\frac{(2l+1)(l-m)!}{4\pi(l+m)!}} P_l^m(\cos\theta) e^{im\phi}}

其中 l=0,1,2,l = 0, 1, 2, \ldotsm=l,l+1,,lm = -l, -l+1, \ldots, l

12.3.2 球谐函数的物理图像与表

球谐函数是原子轨道角向部分的数学描述:

ll名称mm物理形状
0s0球对称
1p-1, 0, +1哑铃形(沿x, y, z)
2d-2, -1, 0, +1, +2四叶草、环形等
3f-3, …, +3更复杂的形状

YlmY_l^m 的模平方 Ylm2|Y_l^m|^2 给出在方向 (θ,ϕ)(\theta, \phi) 找到电子的概率密度。例如:

  • Y00=1/4πY_0^0 = 1/\sqrt{4\pi}:完全球对称(s轨道)
  • Y10=3/(4π)cosθY_1^0 = \sqrt{3/(4\pi)} \cos\theta:沿z轴的哑铃形(p_z轨道)
  • Y1±1=3/(8π)sinθe±iϕY_1^{\pm 1} = \mp \sqrt{3/(8\pi)} \sin\theta e^{\pm i\phi}:在xy平面内旋转的p轨道
graph TD
    A["球坐标分离变量
ψ("\\"r,θ,φ\\"") = R(r)Y("\\"θ,φ\\"")"] --> B["径向方程
依赖于V(r)"] A --> C["角向方程
普适的"] C --> D["L² Y = ℏ²l("\\"l+1\\"")Y"] C --> E["L_z Y = mℏ Y"] D --> F[连带Legendre方程] E --> G["e^{imφ}因子"] F --> H["P_l^m(cosθ)"] G --> I["Y_l^m("\\"θ,φ\\"")"] H --> I I --> J["s: l=0"] I --> K["p: l=1
m=-1,0,+1"] I --> L["d: l=2
m=-2,...,+2"] I --> M["f: l=3, ..."] J --> N[球对称] K --> O[哑铃形] L --> P[四叶草] style I fill:#e3f2fd style C fill:#f3e5f5

12.3.3 宇称与球谐函数

空间反演 rr\mathbf{r} \to -\mathbf{r} 在球坐标中对应 (r,θ,ϕ)(r,πθ,ϕ+π)(r, \theta, \phi) \to (r, \pi - \theta, \phi + \pi)。球谐函数在此变换下的行为:

Ylm(πθ,ϕ+π)=(1)lYlm(θ,ϕ)Y_l^m(\pi - \theta, \phi + \pi) = (-1)^l Y_l^m(\theta, \phi)

宇称为 (1)l(-1)^l

  • ll 为偶数(s, d, g…):偶宇称
  • ll 为奇数(p, f, h…):奇宇称

这在选择定则中起关键作用:电磁跃迁通常要求宇称变化 Δl=±1\Delta l = \pm 1,因此s态和p态之间可以跃迁,但s态到s态是禁戒的(电偶极近似下)。

12.3.4 角动量的叠加:Clebsch-Gordan系数

当两个角动量耦合(例如 J=L+SJ = L + S),总角动量的本征态 j,mj|j, m_j\rangle 可以用单粒子态的直积展开:

j,mj=ml,msl,ml;s,msj,mjl,mls,ms|j, m_j\rangle = \sum_{m_l, m_s} \langle l, m_l; s, m_s | j, m_j \rangle |l, m_l\rangle \otimes |s, m_s\rangle

展开系数 l,ml;s,msj,mj\langle l, m_l; s, m_s | j, m_j \rangle 称为Clebsch-Gordan系数。例如自旋-轨道耦合:

j=l+1/2,mj=l+mj+1/22l+1l,mj1/2+lmj+1/22l+1l,mj+1/2|j = l + 1/2, m_j\rangle = \sqrt{\frac{l + m_j + 1/2}{2l+1}} |l, m_j - 1/2\rangle |\uparrow\rangle + \sqrt{\frac{l - m_j + 1/2}{2l+1}} |l, m_j + 1/2\rangle |\downarrow\rangle

这些系数决定了原子光谱的精细结构,也是核物理和粒子物理中角动量分析的基本工具。

12.3.5 数值例子:球谐函数的具体计算

让我们计算几个低阶球谐函数的具体数值。

Y00Y_0^0(s轨道)

Y00(θ,ϕ)=14π0.282Y_0^0(\theta, \phi) = \frac{1}{\sqrt{4\pi}} \approx 0.282

在所有方向都相同,完全球对称。

Y10Y_1^0(p_z轨道)

Y10(θ,ϕ)=34πcosθY_1^0(\theta, \phi) = \sqrt{\frac{3}{4\pi}} \cos\theta

  • θ=0\theta = 0(北极,z轴正方向):Y10=3/(4π)0.488Y_1^0 = \sqrt{3/(4\pi)} \approx 0.488
  • θ=π/2\theta = \pi/2(赤道面):Y10=0Y_1^0 = 0
  • θ=π\theta = \pi(南极):Y10=3/(4π)0.488Y_1^0 = -\sqrt{3/(4\pi)} \approx -0.488

概率密度 Y102=34πcos2θ|Y_1^0|^2 = \frac{3}{4\pi} \cos^2\theta

  • 在z轴方向最大
  • 在xy平面为零——电子永远不会出现在赤道面上!这就是p轨道的"哑铃形"节点结构。

Y1+1Y_1^{+1}(p_{+1}轨道)

Y1+1(θ,ϕ)=38πsinθeiϕY_1^{+1}(\theta, \phi) = -\sqrt{\frac{3}{8\pi}} \sin\theta \, e^{i\phi}

概率密度:

Y1+12=38πsin2θ|Y_1^{+1}|^2 = \frac{3}{8\pi} \sin^2\theta

  • 在z轴方向(θ=0\theta = 0):为零
  • 在xy平面(θ=π/2\theta = \pi/2):最大,Y1+12=38π0.119|Y_1^{+1}|^2 = \frac{3}{8\pi} \approx 0.119

通过线性组合 Y1+1Y_1^{+1}Y11Y_1^{-1},可以得到实函数的p_x和p_y轨道:

pxsinθcosϕ,pysinθsinϕp_x \propto \sin\theta \cos\phi, \quad p_y \propto \sin\theta \sin\phi

这些实函数在化学中更常用,但复数形式的 YlmY_l^mLzL_z 的本征态,在量子力学中更基本。


本章总结

对称性不是物理学的装饰,而是它的骨架。平移不变性保证动量守恒,旋转不变性保证角动量守恒,时间平移不变性保证能量守恒——这些不是巧合,而是Noether定理在不同对称群下的具体表现。量子力学中的对称变换由幺正算符实现,生成元就是守恒的可观测量。宇称这种分立对称性则揭示了大自然的手性偏好(尤其在弱相互作用中),打破了经典物理中"左右完全等价"的预期。

角动量代数 [Li,Lj]=iεijkLk[L_i, L_j] = i\hbar \varepsilon_{ijk} L_k 是量子力学中最优美的数学结构之一。仅从这个对易关系出发,升降算符方法就导出了本征值的完整量子化:ll 为整数或半整数,mml-lll 以整数步长变化。轨道角动量要求 ll 为整数,源于波函数在 2π2\pi 旋转下的单值性;自旋角动量允许半整数,则是纯粹的量子效应,没有经典对应。

球谐函数 Ylm(θ,ϕ)Y_l^m(\theta, \phi) 作为 L2L^2LzL_z 的共同本征函数,是三维旋转群不可约表示的基矢。它们描述了原子轨道的角向分布,从球对称的s轨道到哑铃形的p轨道,再到更复杂的d和f轨道。宇称 (1)l(-1)^l 在选择定则中扮演关键角色,决定了哪些光谱跃迁是被允许的。

graph TD
    A[对称性] --> B[连续对称]
    A --> C[分立对称]
    B --> D["平移 → P守恒"]
    B --> E["旋转 → L守恒"]
    B --> F["时间平移 → E守恒"]
    C --> G["宇称 Π = ±1"]
    C --> H[时间反演]
    E --> I[角动量代数]
    I --> J["[L_i,L_j] = iℏε_ijk L_k"]
    J --> K[升降算符]
    K --> L["l, m量子化"]
    L --> M["球谐函数 Y_l^m"]
    M --> N["原子轨道
光谱结构"] G --> O[选择定则] style E fill:#e8f5e9 style I fill:#e3f2fd style M fill:#f3e5f5

练习与思考

  1. 宇称守恒的检验:证明如果Hamiltonian在空间反演下不变([Π,H]=0[\Pi, H] = 0),那么非简并的能量本征态必须有确定的宇称。提示:考虑 Πn\Pi |n\ranglen|n\rangle 的关系。

  2. 角动量代数推导:仅使用对易关系 [Li,Lj]=iεijkLk[L_i, L_j] = i\hbar \varepsilon_{ijk} L_k,证明 [L2,Lz]=0[L^2, L_z] = 0。这说明 L2L^2LzL_z 可以有共同本征态。然后验证 L±l,mL_\pm |l, m\rangle 仍然是 L2L^2 的本征态(本征值不变),但 LzL_z 的本征值改变 ±\pm \hbar

  3. 球谐函数的宇称:利用球谐函数的显式表达式,证明 Ylm(πθ,ϕ+π)=(1)lYlm(θ,ϕ)Y_l^m(\pi - \theta, \phi + \pi) = (-1)^l Y_l^m(\theta, \phi)。解释为什么p轨道(l=1l=1)具有奇宇称,而s轨道(l=0l=0)具有偶宇称,并讨论这对原子跃迁选择定则 Δl=±1\Delta l = \pm 1 的物理意义。

  4. 升降算符的数值验证:对于 l=1l=1 系统,利用12.2.4节给出的 Lx,Ly,LzL_x, L_y, L_z 矩阵,直接计算 [Lx,Ly][L_x, L_y] 并验证等于 iLzi\hbar L_z。然后计算 L2=Lx2+Ly2+Lz2L^2 = L_x^2 + L_y^2 + L_z^2,验证它在 1,m|1, m\rangle 基下是对角的,本征值为 222\hbar^2(因为 l(l+1)2=22l(l+1)\hbar^2 = 2\hbar^2)。