第VII章：微扰理论——当精确解成为奢望

"现实世界中，没有什么是完全孤立的。精确解是理想国的居民，微扰理论才是通往现实世界的护照。"

前置知识：渐近展开与微扰理论思想

在深入学习Dirac的微扰理论之前，我们需要理解一个更广泛的数学概念——渐近展开（asymptotic expansion），以及微扰理论在物理学中的深厚传统。

1. 从天体力学到量子力学

微扰理论并非量子力学的专利。它的起源可以追溯到18世纪的天体力学：

牛顿在《自然哲学的数学原理》中已经用迭代方法处理月球运动的微小扰动
欧拉、拉格朗日、拉普拉斯系统发展了行星轨道摄动理论
庞加莱在19世纪末证明了三体问题的不可积性，同时完善了渐近方法

天体力学的核心问题是：太阳-行星两体问题可精确求解，但第三个天体（如木星对地球轨道的摄动）引入了"微小"的修正。数学家发现，即使摄动很小，长期效应可能累积成显著结果——这就是长期项（secular terms）问题。

量子力学的微扰理论继承了这一整套思想，但有了新的数学结构：

经典摄动中，轨道是确定的轨迹，微扰改变轨迹参数
量子微扰中，态矢量是希尔伯特空间中的向量，微扰改变能级和态的成分

2. 渐近展开的本质

一个函数 $f(\epsilon)$ 的渐近展开是指形式幂级数：

$f(\epsilon) \sim \sum_{n=0}^\infty a_n \epsilon^n$

其中符号 $\sim$ （而不是 $=$ ）表示：对任意 $N$ ，

f(\epsilon) - \sum_{n=0}^N a_n \epsilon^n = o(\epsilon^N) \quad \text{当 } \epsilon \to 0

关键性质：

未必收敛：渐近级数通常是发散的！也就是说，如果取 $N \to \infty$ 对固定 $\epsilon$ ，级数可能不收敛
最优截断：存在一个"最优"的 $N_{opt}(\epsilon)$ ，使得部分和最接近真实值。超过此点，加入更多项反而使近似变差
指数小项不可见：形如 $e^{-1/\epsilon}$ 的项在任意阶渐近展开中都不出现，但对精确结果有贡献

graph TD
    A["渐近展开的核心悖论"] --> B["级数发散"]
    A --> C["截断后高度精确"]
    
    B --> D["对于固定 ε, N→∞ 不收敛"]
    C --> E["对于固定 N, ε→0 任意精确"]
    
    D --> F["例如: Stirling级数
    阶乘近似"]
    E --> G["例如: 量子微扰
    前几级给出极好结果"]
    
    A --> H["最优截断 N_opt ~ 1/ε"]
    H --> I["超过 N_opt:
    误差反而增大"]
    
    style A fill:#e3f2fd
    style C fill:#c8e6c9
    style H fill:#ffcc80
    style I fill:#ffcdd2

3. 量子微扰中的渐近性质

量子力学中的微扰展开在形式上是：

$E_n(\lambda) = E_n^{(0)} + \lambda E_n^{(1)} + \lambda^2 E_n^{(2)} + \cdots$

其中 $\lambda$ 是"小参数"。这个级数通常是渐近的而非收敛的。

著名的例子：

量子电动力学中的电子磁矩： $g_e = 2 + \frac{\alpha}{\pi} + O(\alpha^2)$ 。一级修正 $\alpha/\pi \approx 0.00232$ 给出与实验极好的吻合
反常塞曼效应：微扰级数在弱磁场中收敛良好，强磁场中失效（需要帕邢-巴克效应的重新组织）

物理直觉：微扰理论之所以有效，是因为"未微扰"系统已经捕获了物理的主导结构，微扰只是做小修正。当微扰与未微扰系统可比时，需要重新选择"未微扰"的参考点。

4. Dirac的哲学：变换作为微扰

Dirac在《量子力学原理》中对微扰的处理不仅是计算技巧，更是他变换理论的自然应用：

微扰 $V$ 被看作是对未微扰表象的一个连续变换
能级修正是变换的"生成元"的作用
跃迁概率是变换矩阵元的模方
定态微扰和含时微扰被统一到同一个数学结构之下

这种视角将微扰理论从"近似计算"提升为"结构理解"——微扰揭示了未微扰系统背后隐藏的变换结构。

引言：一场星际航行的意外

2147年，人类第一艘深空探测舰"奥伯龙号"在航行至比邻星b轨道的途中，遭遇了一场始料未及的危机。舰上工程师林然正在检查量子导航系统时，警报突然尖啸——飞船的主量子计算机原本设计用于在绝对真空和零磁场环境下运行，但一颗意外的微型陨石击穿了屏蔽层，引入了一个微弱的、随时间变化的电磁场。林然盯着屏幕上跳动的能级图：原本完美的、离散的量子态开始"渗漏"——电子从稳定的轨道跃迁到了其他能级。精确求解？不可能。飞船的导航算法需要在一小时内重新校准，而现有的薛定谔方程精确解完全无法描述这种"闯入者"的影响。林然想起祖父书房里那本泛黄的《量子力学原理》——微扰理论，那个在精确解不可得时，用最优雅的近似捕获真相的工具。他打开了第VII章。

一、微扰理论的核心思想

1.1 为什么要微扰？

在量子力学中，能够精确求解的哈密顿量屈指可数：谐振子、氢原子、自由粒子、自旋在磁场中的进动。然而，真实世界远比这些理想模型复杂。一个原子置于外电场中（斯塔克效应），一个电子在原子核附近感受到弱磁场（塞曼效应），或者一个原子系统与外界电磁场发生相互作用（辐射吸收与发射）——这些情景都无法用已知的精确解直接描述。

Dirac在本书中提出的微扰理论，其核心哲学是：将复杂的哈密顿量分解为可精确求解的部分和"微小"的扰动部分，然后通过系统的方法逐级逼近真实解。

设系统的总哈密顿量为：

$H = H_0 + V$

其中：

$H_0$ 是未微扰哈密顿量（unperturbed Hamiltonian），其本征值和本征态已知或可精确求解
$V$ 是微扰项（perturbation），满足 $|V| \ll |H_0|$ （在某种意义上"足够小"）

这个分解是整个微扰理论的基石。它不是权宜之计，而是量子力学处理真实世界问题的方法论宣言。

1.2 微扰的两种面目

Dirac敏锐地区分了微扰的两种情形：

与时间无关的微扰（Stationary Perturbation）：微扰 $V$ 不随时间变化，主要导致能级的移动和本征态的"混合"。

与时间有关的微扰（Time-Dependent Perturbation）：微扰 $V(t)$ 随时间变化，主要导致系统从一个定态跃迁到另一个定态。

这两种情形在物理上完全不同，但在Dirac的变换理论框架中，它们共享同一套数学基础——微扰引起的态矢量的连续变化。

graph TD
    A["总哈密顿量 H = H₀ + V"] --> B{"微扰类型?"}
    B -->|与时间无关| C[定态微扰理论]
    B -->|与时间有关| D[含时微扰理论]
    C --> C1["能级移动 Eₙ → Eₙ + ΔEₙ"]
    C --> C2["本征态混合 |n⟩ → |n⟩ + Σ cₘ|m⟩"]
    D --> D1["跃迁概率 |n⟩ → |m⟩"]
    D --> D2["吸收/发射辐射"]
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#0f3460,color:#fff

二、定态微扰理论：能级的精细移动

2.1 一级能量修正

假设未微扰系统满足：

$H_0 |n^{(0)}\rangle = E_n^{(0)} |n^{(0)}\rangle$

其中 $|n^{(0)}\rangle$ 构成一组完备的正交归一基。

当加入微小且恒定的微扰 $V$ 后，真实的本征态和本征值可以展开为微扰的幂级数：

$|n\rangle = |n^{(0)}\rangle + |n^{(1)}\rangle + |n^{(2)}\rangle + \cdots$

$E_n = E_n^{(0)} + E_n^{(1)} + E_n^{(2)} + \cdots$

其中上标表示微扰的阶数。

在一级近似下，能量修正为：

$E_n^{(1)} = \langle n^{(0)} | V | n^{(0)} \rangle$

这个结果的物理意义极其清晰：一级能量修正就是微扰算符在未微扰态中的期望值。它告诉我们，微扰如何"平均地"改变了系统的能量。

物理直觉：想象一个完美圆形的轨道上运行的行星（ $H_0$ ），突然引入了一个微小的"地形起伏"（ $V$ ）。行星能量的一级修正，就等于这个地形起伏在整个轨道上的平均值。

2.2 一级态矢量修正与混合

能量的修正只是故事的一半。微扰不仅改变能级，还会"扭曲"本征态：

$|n^{(1)}\rangle = \sum_{m \neq n} \frac{\langle m^{(0)} | V | n^{(0)} \rangle}{E_n^{(0)} - E_m^{(0)}} |m^{(0)}\rangle$

这个公式的分子 $\langle m^{(0)} | V | n^{(0)} \rangle = V_{mn}$ 是微扰在未微扰态之间的矩阵元，它衡量了微扰"耦合"不同本征态的强度。

分母 $E_n^{(0)} - E_m^{(0)}$ 则是两个未微扰能级的能量差。这里藏着微扰理论的一个关键洞察：

能级越接近，微扰引起的混合越强。

当两个能级简并（ $E_n^{(0)} = E_m^{(0)}$ ）时，上述公式失效——分母为零，发散。这正是简并微扰理论需要单独处理的原因。

graph LR
    subgraph 未微扰系统
        A1["|1⟩ E₁⁽⁰⁾"] ---|大能隙| B1["|2⟩ E₂⁽⁰⁾"]
        C1["|3⟩ E₃⁽⁰⁾"] ---|小能隙| D1["|4⟩ E₄⁽⁰⁾"]
    end
    subgraph 加入微扰V
        A2["|1⟩ E₁⁽⁰⁾+ΔE₁"] ---|微弱混合| B2["|2⟩ E₂⁽⁰⁾+ΔE₂"]
        C2["|3⟩ E₃⁽⁰⁾+ΔE₃"] ---|强烈混合| D2["|4⟩ E₄⁽⁰⁾+ΔE₄"]
    end
    style A1 fill:#1a1a2e,color:#fff
    style B1 fill:#1a1a2e,color:#fff
    style C1 fill:#16213e,color:#fff
    style D1 fill:#16213e,color:#fff
    style A2 fill:#0f3460,color:#fff
    style B2 fill:#0f3460,color:#fff
    style C2 fill:#e94560,color:#fff
    style D2 fill:#e94560,color:#fff

2.3 二级能量修正

二级能量修正捕获了微扰的更深层效应：

$E_n^{(2)} = \sum_{m \neq n} \frac{|\langle m^{(0)} | V | n^{(0)} \rangle|^2}{E_n^{(0)} - E_m^{(0)}}$

注意分子是模方（正数），但分母可正可负：

如果 $m$ 的能级高于 $n$ （ $E_m^{(0)} > E_n^{(0)}$ ），则分母为负，这一项降低基态能量
这与化学中的变分原理一致：微扰总是让基态能量更低（至少在一级近似中如此）

二级修正的物理图像是：微扰让系统"借用"了其他态的成分，从而进一步调整能量。

三、含时微扰理论与量子跃迁

3.1 相互作用绘景：舞台的转换

这是Dirac最伟大的贡献之一。在处理含时微扰时，Dirac不直接在薛定谔绘景中工作，而是引入相互作用绘景（Interaction Picture）——一种介于薛定谔绘景和海森堡绘景之间的中间舞台。

在相互作用绘景中，态矢量随时间的演化只由微扰 $V$ 驱动：

$i\hbar \frac{d}{dt} |\psi(t)\rangle_I = V_I(t) |\psi(t)\rangle_I$

其中：

$V_I(t) = e^{iH_0 t/\hbar} V e^{-iH_0 t/\hbar}$

这个看似简单的变换背后隐藏着深刻的物理：我们将已知的、由 $H_0$ 驱动的"自由演化"从问题中剥离，只关注微扰引起的"额外"演化。

用Dirac自己的话说，这是将"已知"和"未知"清晰分离的方法论。

graph TB
    subgraph 三种绘景
        A["薛定谔绘景
态矢量演化
算符静止"] 
        B["海森堡绘景
态矢量静止
算符演化"]
        C["相互作用绘景
态矢量由V驱动
算符由H₀驱动"]
    end
    A ---|"极端：算符不变"| C
    B ---|"极端：态不变"| C
    style A fill:#1a1a2e,color:#fff
    style B fill:#1a1a2e,color:#fff
    style C fill:#e94560,color:#fff

3.2 跃迁概率的诞生

假设系统在 $t=0$ 时处于 $H_0$ 的第 $n$ 个本征态：

$|\psi(0)\rangle = |n\rangle$

在微扰开启后，系统不再停留在 $|n\rangle$ ，而是获得了"游荡"到其他态的可能性。在时刻 $t$ ，系统处于态 $|m\rangle$ 的概率振幅为：

$c_m(t) = \langle m | \psi(t) \rangle$

在一级微扰近似下：

$c_m^{(1)}(t) = \frac{1}{i\hbar} \int_0^t \langle m | V_I(t') | n \rangle dt'$

或者写回薛定谔绘景：

$c_m^{(1)}(t) = \frac{1}{i\hbar} \int_0^t e^{i(E_m - E_n)t'/\hbar} \langle m | V(t') | n \rangle dt'$

这个公式的结构揭示了量子跃迁的核心机制：

跃迁振幅 = 微扰矩阵元 × 含时相位因子的积分

3.3 与时间无关的微扰引起的跃迁

当微扰 $V$ 不随时间变化时，上述积分可以精确完成：

$c_m^{(1)}(t) = \frac{\langle m | V | n \rangle}{i\hbar} \int_0^t e^{i\omega_{mn}t'} dt' = \frac{\langle m | V | n \rangle}{i\hbar} \cdot \frac{e^{i\omega_{mn}t} - 1}{i\omega_{mn}}$

其中 $\omega_{mn} = (E_m - E_n)/\hbar$ 。

跃迁概率：

$P_{n\to m}(t) = |c_m^{(1)}(t)|^2 = \frac{|\langle m | V | n \rangle|^2}{\hbar^2} \cdot \frac{\sin^2(\omega_{mn}t/2)}{(\omega_{mn}/2)^2}$

这个结果的物理内容极其丰富：

短时间行为：当 $t \ll 1/\omega_{mn}$ 时， $P \propto t^2$ ——概率随时间的平方增长
共振特性：函数 $\sin^2(xt/2)/(x/2)^2$ 在 $x = 0$ 处取最大值 $t^2$ ，随着 $|x|$ 增大快速振荡衰减
能量守恒：当 $t \to \infty$ 时，这个函数趋于 $\delta(\omega_{mn})$ ——跃迁只发生在能量相近的态之间

Dirac的洞察：即使微扰是时间无关的，如果它"突然开启"，系统会经历一个瞬态过程，在这个过程中不同能级之间会发生跃迁。最终，系统会弛豫到新的定态，但中间过程的跃迁概率由上述公式描述。

3.4 费米黄金定则

如果微扰是恒定的（ $V$ 不随时间变化），且在足够长时间后观测，上述积分会产生一个尖锐的频率选择特性：

$c_m(t) \propto \frac{\sin[(E_m - E_n)t/2\hbar]}{(E_m - E_n)/2\hbar}$

当 $t \to \infty$ 时，这个函数趋于 $\delta(E_m - E_n)$ ——能量守恒从跃迁概率中自然涌现！

对于连续谱或准连续谱，跃迁速率（单位时间的跃迁概率）成为更有意义的量。在一级微扰下，这就是著名的费米黄金定则（Fermi’s Golden Rule）：

$\Gamma_{n \to m} = \frac{2\pi}{\hbar} |\langle m | V | n \rangle|^2 \rho(E_m)$

其中：

$\Gamma_{n \to m}$ 是从 $|n\rangle$ 到 $|m\rangle$ 的跃迁速率
$\rho(E_m)$ 是末态的态密度（单位能量的态数目）
$|\langle m | V | n \rangle|^2$ 是微扰矩阵元的模方

flowchart LR
    A["初态 |n⟩
能量Eₙ"] -->|微扰V| B["末态 |m⟩
能量Eₘ"]
    B --> C{"Eₘ ≈ Eₙ?"}
    C -->|是| D["跃迁允许
费米黄金定则"]
    C -->|否| E["跃迁被抑制
能量不守恒"]
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#e94560,color:#fff
    style E fill:#533483,color:#fff

四、辐射与物质的相互作用：光子的诞生

4.1 电磁场作为微扰

当一个原子置于电磁场中时，哈密顿量变为：

$H = \frac{1}{2m}(\mathbf{p} - q\mathbf{A})^2 + q\phi + V_{\text{atomic}}$

对于弱场（ $|\mathbf{A}|$ 小），可以展开得到微扰项：

$V \approx -\frac{q}{m} \mathbf{p} \cdot \mathbf{A}$

这是最小耦合（minimal coupling）的结果。它告诉我们：电子通过动量算符与矢势耦合。

4.2 吸收与发射

考虑一个单模电磁场，频率为 $\omega$ ，原子有两个能级 $|g\rangle$ （基态）和 $|e\rangle$ （激发态），能级差 $E_e - E_g = \hbar \omega_0$ 。

微扰矩阵元为：

$\langle e | V | g \rangle \propto \langle e | \mathbf{p} \cdot \mathbf{A} | g \rangle$

当入射光的频率 $\omega \approx \omega_0$ 时，发生共振：原子吸收光子，从基态跃迁到激发态——这就是受激吸收。

反之，如果原子初始处于激发态，微扰同样可以诱导它跃迁到基态，释放光子——受激发射。

Dirac在1927年的论文中首次用这种方法推导出了受激发射系数，后来成为激光理论的基石。更令人惊叹的是，即使在没有任何光子的真空中，微扰理论也预言了自发辐射——激发态原子会自发地发射光子。这不是微扰理论的失败，而是电磁场量子化的必然结果：真空并非"空无一物"，而是充满了量子涨落。

sequenceDiagram
    participant Atom as 原子
    participant Field as 电磁场
    Note over Atom,Field: 共振条件 ℏω = Eₑ - E_g
    
    rect rgb("26,26,46")
        Note over Atom,Field: 受激吸收
        Field->>Atom: 入射光子 ℏω
        Atom->>Atom: |g⟩ → |e⟩
    end
    
    rect rgb("15,33,96")
        Note over Atom,Field: 受激发射
        Atom->>Field: 发射光子 ℏω
        Atom->>Atom: |e⟩ → |g⟩
        Field->>Field: 光子数 +1
相干放大
    end
    
    rect rgb("233,69,96")
        Note over Atom,Field: 自发辐射
        Atom->>Field: 发射光子 ℏω₀
        Atom->>Atom: |e⟩ → |g⟩
（真空涨落驱动）
    end

4.3 选择定则

并非所有跃迁都被允许。跃迁矩阵元 $\langle e | \mathbf{p} | g \rangle$ 要求初末态的角动量满足特定关系——这就是选择定则的起源：

电偶极跃迁： $\Delta l = \pm 1$ ， $\Delta m = 0, \pm 1$
如果矩阵元为零，该跃迁在一级近似下被"禁戒"（forbidden）

这些选择定则深刻影响了原子光谱的结构，也是化学中电子跃迁规则的理论来源。

五、反常塞曼效应：自旋的舞台

5.1 正常与反常

当原子置于弱磁场 $\mathbf{B}$ 中时，能级会发生分裂。如果忽略自旋（或自旋为零），简单的分析给出：

$\Delta E = \mu_B B m_l$

其中 $\mu_B = e\hbar/2m_e$ 是玻尔磁子， $m_l$ 是轨道角动量的磁量子数。每条能级分裂为 $2l+1$ 条——这就是正常塞曼效应。

然而，实验观测到的许多原子（尤其是碱金属）却表现出更复杂的分裂模式——反常塞曼效应。能级分裂的数目和间隔都无法用纯轨道角动量解释。

5.2 自旋的引入

Dirac的微扰理论为这个问题提供了完美的解决框架。当把自旋-轨道耦合和磁场的作用同时作为微扰处理时，能级修正为：

$\Delta E = \mu_B B (m_l + 2m_s)$

注意自旋的贡献系数是 2（准确的说是朗德g因子），而非轨道角动量的 1。这个"2"的自旋磁矩比，在Dirac的相对论量子力学中才能被完美解释（它来自Dirac方程的自然推导）。

反常塞曼效应的完整分析需要使用简并微扰理论：当自旋-轨道耦合与磁场相比可比拟时，原本的简并能级需要先被总角动量 $J = L + S$ 的本征态"预对角化"，然后再计算磁场引起的分裂。

5.3 朗德g因子

在LS耦合方案中，朗德g因子给出了反常塞曼分裂的有效磁矩：

$g_J = 1 + \frac{J(J+1) + S(S+1) - L(L+1)}{2J(J+1)}$

这个公式完美解释了为什么有些能级分裂是正常塞曼间隔的 $\frac{2}{3}$ ，有些是 $2$ 倍，有些是 $\frac{4}{3}$ 倍。它把整个原子光谱学的经验规则凝结为一个优雅的量子力学公式。

5.4 帕邢-巴克效应

在极强磁场中，自旋-轨道耦合被磁场压制，系统回到"正常塞曼效应"的行为——这是两个微扰竞争的经典例子。

graph TD
    A["无磁场
精细结构"] -->|弱磁场B| B["反常塞曼效应
E = E₀ + μ_Bg_Jm_JB"]
    B -->|中等磁场| C["帕邢-巴克效应
LS耦合被破坏"]
    C -->|强磁场| D["正常塞曼效应
L和S独立进动"]
    style A fill:#1a1a2e,color:#fff
    style B fill:#0f3460,color:#fff
    style C fill:#533483,color:#fff
    style D fill:#e94560,color:#fff

六、Dirac的数学视角：变换理论中的微扰

Dirac对微扰理论的处理方式比薛定谔和海森堡更加数学化和统一。在Dirac的变换理论（Transformation Theory）框架中：

微扰被看作是对未微扰表象的一个连续变换
能级修正是这个变换的"生成元"的作用
跃迁概率是变换矩阵元的模方

这种视角将定态微扰和含时微扰统一到同一个数学结构之下：

$U(t) = T \exp\left(-\frac{i}{\hbar} \int_0^t V_I(t') dt'\right)$

其中 $T$ 是时序算符（time-ordering operator），确保较早时刻的算符出现在右边（薛定谔绘景中）。这个指数展开：

$U(t) = 1 - \frac{i}{\hbar} \int_0^t V_I(t') dt' + \left(-\frac{i}{\hbar}\right)^2 \int_0^t dt_1 \int_0^{t_1} dt_2 V_I(t_1)V_I(t_2) + \cdots$

正是各级微扰的积分形式。第零级是"无微扰"，第一级是单光子过程（吸收或发射），第二级是双光子过程（散射、虚跃迁），依此类推。

深刻洞察：现代量子场论中的费曼图，本质上就是对这个时序指数展开的图形化表示。

graph LR
    subgraph Dyson级数展开
        A["U = 1"] --- B["+ V₁"]
        B --- C["+ V₂"]
        C --- D["+ V₃"]
        D --- E["..."]
    end
    subgraph 物理对应
        V1["一级
单光子过程"]
        V2["二级
双光子/虚跃迁"]
        V3["三级
三光子过程"]
    end
    B -.-> V1
    C -.-> V2
    D -.-> V3
    style A fill:#1a1a2e,color:#fff
    style B fill:#16213e,color:#fff
    style C fill:#0f3460,color:#fff
    style D fill:#533483,color:#fff
    style E fill:#e94560,color:#fff

数值例子

例1：一维谐振子的微扰能级移动

问题：一个带电粒子（电荷 $q = e$ ，质量 $m = m_e$ ）在一维谐振子势 $V_0(x) = \frac{1}{2}m\omega^2 x^2$ 中运动，其中 $\hbar\omega = 1$ eV。现在施加一个恒定弱电场 $\mathcal{E} = 10^4$ V/m，沿 $x$ 方向。计算基态和第一激发态的一级和二级能量修正。

解答：

微扰项：

$V = -q\mathcal{E}x = -e\mathcal{E}x$

参数：

$\hbar\omega = 1$ eV $= 1.602 \times 10^{-19}$ J
$\omega = \frac{\hbar\omega}{\hbar} = \frac{1.602 \times 10^{-19}}{1.055 \times 10^{-34}} = 1.519 \times 10^{15}$ rad/s
特征长度： $a = \sqrt{\frac{\hbar}{m\omega}} = \sqrt{\frac{1.055 \times 10^{-34}}{9.11 \times 10^{-31} \times 1.519 \times 10^{15}}} = \sqrt{7.62 \times 10^{-20}} = 8.73 \times 10^{-10}$ m $= 0.873$ Å

一级能量修正：

$E_n^{(1)} = \langle n | V | n \rangle = -e\mathcal{E} \langle n | x | n \rangle$

由于宇称： $|n\rangle$ 的宇称为 $(-1)^n$ ， $x$ 是奇宇称算符，所以 $\langle n | x | n \rangle = 0$ 。

因此： $E_0^{(1)} = 0$ ， $E_1^{(1)} = 0$

物理意义：线性微扰的一级修正对谐振子所有能级都为零，因为微扰（线性势）在谐振子本征态中的期望值为零。

二级能量修正：

$E_n^{(2)} = \sum_{m \neq n} \frac{|\langle m | V | n \rangle|^2}{E_n^{(0)} - E_m^{(0)}} = (e\mathcal{E})^2 \sum_{m \neq n} \frac{|\langle m | x | n \rangle|^2}{(n-m)\hbar\omega}$

利用升降算符： $x = \sqrt{\frac{\hbar}{2m\omega}}(\hat{a} + \hat{a}^\dagger) = \frac{a}{\sqrt{2}}(\hat{a} + \hat{a}^\dagger)$

矩阵元：

$\langle n+1 | x | n \rangle = \sqrt{\frac{\hbar(n+1)}{2m\omega}} = a\sqrt{\frac{n+1}{2}}$
$\langle n-1 | x | n \rangle = \sqrt{\frac{\hbar n}{2m\omega}} = a\sqrt{\frac{n}{2}}$
其他矩阵元为零

基态 ( $n=0$ ) 二级修正：

只有 $m=1$ 项贡献：

$E_0^{(2)} = \frac{(e\mathcal{E})^2 |\langle 1 | x | 0 \rangle|^2}{E_0^{(0)} - E_1^{(0)}} = \frac{(e\mathcal{E})^2 \cdot \frac{\hbar}{2m\omega}}{-\hbar\omega} = -\frac{(e\mathcal{E})^2}{2m\omega^2}$

代入数值：

$e = 1.602 \times 10^{-19}$ C
$\mathcal{E} = 10^4$ V/m
$m = 9.11 \times 10^{-31}$ kg
$\omega = 1.519 \times 10^{15}$ rad/s

$E_0^{(2)} = -\frac{(1.602 \times 10^{-19} \times 10^4)^2}{2 \times 9.11 \times 10^{-31} \times (1.519 \times 10^{15})^2}$

$= -\frac{(1.602 \times 10^{-15})^2}{2 \times 9.11 \times 10^{-31} \times 2.307 \times 10^{30}}$

$= -\frac{2.566 \times 10^{-30}}{4.204 \times 10^{0}} = -6.10 \times 10^{-31}\text{ J}$

转换为eV：

$E_0^{(2)} = -\frac{6.10 \times 10^{-31}}{1.602 \times 10^{-19}} = -3.81 \times 10^{-12}\text{ eV}$

第一激发态 ( $n=1$ ) 二级修正：

有两项贡献： $m=0$ 和 $m=2$

$E_1^{(2)} = \frac{(e\mathcal{E})^2|\langle 0|x|1\rangle|^2}{\hbar\omega - 0} + \frac{(e\mathcal{E})^2|\langle 2|x|1\rangle|^2}{\hbar\omega - 2\hbar\omega}$

$= \frac{(e\mathcal{E})^2 \frac{\hbar}{2m\omega}}{\hbar\omega} + \frac{(e\mathcal{E})^2 \frac{\hbar}{m\omega}}{-\hbar\omega}$

$= \frac{(e\mathcal{E})^2}{2m\omega^2} - \frac{(e\mathcal{E})^2}{m\omega^2} = -\frac{(e\mathcal{E})^2}{2m\omega^2}$

所以： $E_1^{(2)} = E_0^{(2)} = -3.81 \times 10^{-12}$ eV

物理讨论：

二级修正对所有能级相同（这是谐振子的特殊性质），能级间距不变，整体下移
修正量 $-3.81 \times 10^{-12}$ eV 相对于 $\hbar\omega = 1$ eV 是极小的（ $\sim 10^{-12}$ ），验证了微扰理论的适用性
这个结果可以通过精确求解验证：在外电场中，谐振子势能变为 $\frac{1}{2}m\omega^2 x^2 - e\mathcal{E}x = \frac{1}{2}m\omega^2(x-x_0)^2 - \frac{(e\mathcal{E})^2}{2m\omega^2}$ ，其中 $x_0 = \frac{e\mathcal{E}}{m\omega^2}$ 。精确解给出的能级移动正是 $-\frac{(e\mathcal{E})^2}{2m\omega^2}$ ，与二级微扰结果一致！

例2：氢原子的斯塔克效应（线性斯塔克效应）

问题：氢原子处于基态，施加沿 $z$ 方向的均匀电场 $\mathcal{E} = 10^5$ V/m。计算一级能量修正。如果氢原子处于 $n=2$ 能级，一级修正又如何？

解答：

微扰项：

$V = e\mathcal{E}z = e\mathcal{E}r\cos\theta$

基态 ( $n=1, l=0, m=0$ )：

$E_{1s}^{(1)} = \langle 1s | e\mathcal{E}z | 1s \rangle$

$|1s\rangle$ 是球对称的（ $l=0$ ）， $z = r\cos\theta$ 在 $\theta$ 积分中：

$\int_0^\pi \cos\theta \sin\theta d\theta = 0$

因此： $E_{1s}^{(1)} = 0$

物理意义：基态氢原子没有永久电偶极矩（球对称），所以线性斯塔克效应为零。二级效应存在但极小。

$n=2$ 能级：

$n=2$ 能级有 $l=0$ （ $2s$ ）和 $l=1$ （ $2p$ ， $m=-1,0,+1$ ）四个简并态。简并微扰理论要求在这四个态构成的子空间中 diagonalize $V$ 。

矩阵元 $\langle 2l'm' | e\mathcal{E}r\cos\theta | 2lm \rangle$ ：

利用选择定则（ $\Delta l = \pm 1, \Delta m = 0$ ），只有以下矩阵元非零：

$\langle 2s | V | 2p, m=0 \rangle$
$\langle 2p, m=0 | V | 2s \rangle$

计算 $\langle 2s | r\cos\theta | 2p_0 \rangle$ ：

$|2s\rangle = R_{20}(r)Y_0^0 = \frac{1}{\sqrt{2a_0^3}}\left(1-\frac{r}{2a_0}\right)e^{-r/(2a_0)} \cdot \frac{1}{\sqrt{4\pi}}$

$|2p_0\rangle = R_{21}(r)Y_1^0 = \frac{1}{2\sqrt{6a_0^3}}\frac{r}{a_0}e^{-r/(2a_0)} \cdot \sqrt{\frac{3}{4\pi}}\cos\theta$

矩阵元：

$\langle 2s | r\cos\theta | 2p_0 \rangle = \int_0^\infty r^2 dr \int d\Omega R_{20}(r)Y_0^0 \cdot r\cos\theta \cdot R_{21}(r)Y_1^0$

角向积分：

$\int Y_0^0 \cos\theta Y_1^0 d\Omega = \frac{1}{\sqrt{4\pi}} \cdot \sqrt{\frac{3}{4\pi}} \int \cos^2\theta d\Omega = \frac{\sqrt{3}}{4\pi} \cdot \frac{4\pi}{3} = \frac{1}{\sqrt{3}}$

径向积分：

$\int_0^\infty r^3 R_{20}(r) R_{21}(r) dr = \int_0^\infty r^3 \cdot \frac{1}{\sqrt{2}a_0^{3/2}}\left(1-\frac{r}{2a_0}\right) \cdot \frac{1}{2\sqrt{6}a_0^{3/2}}\frac{r}{a_0} e^{-r/a_0} dr$

$= \frac{1}{4\sqrt{3}a_0^4} \int_0^\infty r^4\left(1-\frac{r}{2a_0}\right)e^{-r/a_0} dr$

利用 $\int_0^\infty x^n e^{-ax} dx = n!/a^{n+1}$ ：

$= \frac{1}{4\sqrt{3}a_0^4} \left[\frac{4!}{(1/a_0)^5} - \frac{1}{2a_0} \frac{5!}{(1/a_0)^6}\right]$

$= \frac{1}{4\sqrt{3}a_0^4} \left[24a_0^5 - 60a_0^5\right] = \frac{1}{4\sqrt{3}a_0^4} (-36a_0^5) = -\frac{9a_0}{\sqrt{3}} = -3\sqrt{3}a_0$

因此：

$\langle 2s | r\cos\theta | 2p_0 \rangle = \frac{1}{\sqrt{3}} \times (-3\sqrt{3}a_0) = -3a_0$

在 $n=2$ 子空间中，微扰矩阵为（以 $|2s\rangle, |2p_{+1}\rangle, |2p_0\rangle, |2p_{-1}\rangle$ 为基）：

$V = e\mathcal{E}\begin{pmatrix} 0 & 0 & -3a_0 & 0 \\ 0 & 0 & 0 & 0 \\ -3a_0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix}$

这个矩阵的本征值为： $\pm 3e\mathcal{E}a_0, 0, 0$

因此 $n=2$ 能级分裂为三条：

$E = E_2 + 3e\mathcal{E}a_0$
$E = E_2$
$E = E_2 - 3e\mathcal{E}a_0$

代入数值：

$e = 1.602 \times 10^{-19}$ C
$\mathcal{E} = 10^5$ V/m
$a_0 = 0.529 \times 10^{-10}$ m

$3e\mathcal{E}a_0 = 3 \times 1.602 \times 10^{-19} \times 10^5 \times 0.529 \times 10^{-10}$

$= 3 \times 1.602 \times 0.529 \times 10^{-24} = 2.54 \times 10^{-24}\text{ J}$

转换为eV：

\Delta E = \frac{2.54 \times 10^{-24}}{1.602 \times 10^{-19}} = 1.59 \times 10^{-5}\text{ eV} = 15.9\text{ μeV}

物理讨论：

$n=2$ 能级分裂为三条，间隔约 15.9 μeV
这被称为线性斯塔克效应（能级移动正比于电场的一次方），是简并态微扰理论的结果
$2p_{\pm 1}$ 态不受一级微扰影响，因为它们不能与 $2s$ 或 $2p_0$ 通过 $z$ 方向的电偶极耦合（ $\Delta m \neq 0$ ）
对应的态是 $|2s\rangle \pm |2p_0\rangle$ 的叠加态——这些叠加态具有沿 $z$ 方向的永久电偶极矩 $\pm 3ea_0$

本章总结

Dirac的第VII章将微扰理论从一种计算技巧提升为量子力学理解世界的方法论。我们学习了：

渐近展开的性质：微扰级数通常是渐近而非收敛的，存在最优截断。微扰理论的有效性不在于"无穷级数求和"，而在于"前几项给出极好近似"。
定态微扰理论：当微扰不随时间变化时，系统能级发生精细移动，本征态发生混合。一级能量修正等于微扰的期望值，二级修正涉及与其他态的虚拟"borrowing"。
含时微扰理论：在相互作用绘景中，微扰驱动态矢量演化，产生量子跃迁。与时间无关的微扰如果突然开启，也会引发瞬态跃迁过程。费米黄金定则给出了跃迁速率的普适公式，能量守恒从长时间极限中自然涌现。
辐射理论：微扰理论统一描述了受激吸收、受激发射和自发辐射，揭示了原子与光场相互作用的全貌。
反常塞曼效应：通过微扰理论引入自旋，解释了正常与反常塞曼效应的过渡，朗德g因子成为原子光谱的精确预言工具。
变换理论的统一视角：Dirac将微扰视为连续变换，Dyson级数成为现代量子场论和费曼图技术的源头。

mindmap
  root("(#quot;微扰理论
第VII章#quot;)")
    渐近展开
      发散但有效
      最优截断
      指数小项
    定态微扰
      一级能量修正
⟨n|V|n⟩
      一级态修正
 borrowing
      二级能量修正
      简并微扰
    含时微扰
      相互作用绘景
      跃迁振幅
      时间无关微扰的瞬态
      费米黄金定则
      能量守恒
    辐射应用
      受激吸收
      受激发射
      自发辐射
      选择定则
    反常塞曼效应
      自旋-轨道耦合
      朗德g因子
      LS耦合
      帕邢-巴克效应
    变换理论视角
      Dyson级数
      时序算符
      费曼图起源

练习与思考

1. 二维谐振子的微扰

考虑一个二维各向同性谐振子，哈密顿量为 $H_0 = \frac{p_x^2 + p_y^2}{2m} + \frac{1}{2}m\omega^2(x^2 + y^2)$ 。加入微扰 $V = \lambda xy$ （其中 $\lambda \ll m\omega^2$ ）。用一级定态微扰理论计算基态和第一激发态的能量修正，并讨论简并态的混合。

2. 费米黄金定则的推导

从一阶含时微扰公式出发，假设微扰 $V$ 为常数，证明在 $t \to \infty$ 时，向连续末态的跃迁速率由费米黄金定则给出。说明为什么态密度 $\rho(E)$ 在这个过程中是不可或缺的。

3. 自发辐射的必然性

Dirac的微扰理论中，即使在真空中（无入射光子），激发态原子仍有非零的跃迁概率。这看起来像是微扰理论的"自发"结果。请从电磁场量子化的角度解释：为什么真空并非"空无一物"，以及自发辐射速率如何与真空涨落的强度相关联。这与爱因斯坦在1917年基于热力学平衡提出的A、B系数有什么关系？

4. 简并微扰的二维谐振子

考虑二维谐振子（ $\omega_x = \omega_y = \omega$ ）的 $n=2$ 能级（即总能量 $E = 2\hbar\omega$ ，不含零点能）。这个能级的简并度是多少？列出所有简并态。如果施加微扰 $V = \lambda x^2$ ，用简并微扰理论计算能级分裂。

"精确解是理想国的居民，微扰理论才是通往现实世界的护照。"

✍️ 本章记录于"奥伯龙号"事件后第3小时。林然仍在调试导航系统，但微扰理论的代码已经跑通了——在精确与近似之间，在理想与现实之间，量子力学给出了最优雅的妥协方案。❤️‍🔥