第17-18章微扰理论：渐进展开的艺术

"世界不是完全可解的，但世界足够温和，允许我们一点点接近真相。"

前置知识：渐近展开的思想

微扰理论的核心数学工具是渐近展开（asymptotic expansion）。它与收敛级数（如泰勒级数）有本质区别，理解这一点对于正确使用微扰理论至关重要。

A.1 收敛级数 vs 渐近级数

收敛级数： $\sum_{n=0}^{\infty} a_n x^n$ 在 $|x| < R$ 时收敛于某个有限值。例如 $e^x = \sum_{n=0}^{\infty} x^n/n!$ ，对所有 $x$ 收敛。

渐近级数： $\sum_{n=0}^{N} a_n \epsilon^n$ 当 $\epsilon \to 0$ 时，部分和在最优截断阶数 $N_{opt}$ 给出最好的近似，但无穷级数发散。

数学定义：级数 $\sum a_n \epsilon^n$ 是函数 $f(\epsilon)$ 的渐近展开，如果对于任意 $N$ ：

$\lim_{\epsilon \to 0} \frac{f(\epsilon) - \sum_{n=0}^{N} a_n \epsilon^n}{\epsilon^N} = 0$

A.2 为什么微扰级数是渐近的？

考虑一个简单模型：

$E(\lambda) = \int_{-\infty}^{\infty} e^{-x^2 - \lambda x^4} dx$

形式展开：

$E(\lambda) = \int_{-\infty}^{\infty} e^{-x^2}\left(1 - \lambda x^4 + \frac{\lambda^2 x^8}{2} - \cdots\right)dx$

$= \sqrt{\pi}\left(1 - \frac{3\lambda}{4} + \frac{105\lambda^2}{32} - \cdots\right)$

但这个级数对任何 $\lambda \neq 0$ 都发散！因为高斯积分 $\int x^{4n} e^{-x^2}dx$ 增长极快，系数 $a_n \sim (-1)^n (4n)!/(2^{2n}(2n)!)$ 的增长速度超过任何几何级数。

然而，对于足够小的 $\lambda$ 和有限阶 $N$ ，部分和给出极好的近似。例如 $\lambda = 0.01$ ：

$N=0$ ： $E \approx \sqrt{\pi} = 1.772$
$N=1$ ： $E \approx 1.772(1 - 0.0075) = 1.759$
$N=2$ ： $E \approx 1.772(1 - 0.0075 + 0.000328) = 1.760$

精确值 $E(0.01) \approx 1.760$ 。到 $N=2$ 已经极好，但 $N=10$ 时级数项会变得巨大，近似反而变差。

A.3 Dyson论证与微扰理论的物理边界

Freeman Dyson在1952年给出了一个深刻的论证：考虑量子电动力学（QED）中电子的反常磁矩。如果耦合常数 $e^2$ 变为负值（虚构的世界），库仑势从吸引变为排斥，电子不再束缚于原子核，"真空"本身不稳定——电子-正电子对会不断产生。

这意味着QED的微扰级数在 $e^2 = 0$ 附近没有解析性，级数不可能收敛。但QED的有限阶计算（如 $e^2$ 、 $e^4$ 、 $e^6$ 项）给出了与实验惊人一致的结果（电子反常磁矩的理论与实验吻合到10位有效数字）。

核心结论：

微扰级数通常是渐近级数而非收敛级数
最优截断通常在低阶（1-3阶）
物理上，这反映了一个深刻事实：微扰展开 capture 的是"足够接近可解模型"的物理，而非"任意接近"

graph TD
    A["微扰参数 $\lambda$"] --> B{"$\lambda$ 足够小?"}
    B -->|是| C["低阶截断: 好近似"]
    B -->|否| D[高阶发散]
    
    C --> E["最优阶数 $N_{opt}$"]
    E --> F[误差最小]
    F --> G[物理结果可靠]
    
    D --> H[级数项增大]
    H --> I[近似恶化]
    I --> J[需要非微扰方法]
    
    G -.->|QED实例| K["10位精度!"]

A.4 微扰理论适用的判据

对于能级 $E_n$ ，微扰理论的适用条件是微扰引起的能级移动远小于相邻能级间距：

$|\Delta E_n^{(1)}| \ll |E_{n+1}^{(0)} - E_n^{(0)}|$

对于连续谱（散射问题），判据是微扰矩阵元远小于能量：

$|\langle f|\hat{V}|i\rangle| \ll E$

当这些条件不满足时（如简并微扰、共振散射、强耦合系统），需要发展特殊的非微扰或重求和技术。

故事场景：调音师的耳朵

公元2201年，"共鸣城"的中央音乐厅有一架三百岁的管风琴出了问题。这架琴有上万个音管，每一个都受到邻近音管气流的微弱影响。年轻的调音师小李想要逐个拆下来测试，但首席调音师老周阻止了他。

"你不能把世界拆开再装回去，"老周说。他拿出一张写满公式的羊皮纸——那是古代物理学家发明的"渐进展开法"。

"每一个音管的’真实音高’，都等于它单独存在时的音高，加上其他所有音管对它造成的’微扰’。如果把这些影响按照大小排序，你可以先算最大的，再算次大的，以此类推。只要这些扰动足够小，第三项之后的修正就可以忽略了。"

小李问："如果两个音管的影响一样大呢？"

老周的表情变得严肃："那就是简并微扰——最危险的情况。你必须同时处理它们，不能分开算。"

这就是Shankar在第17-18章讲述的故事：微扰理论——当真实世界与可解模型之间只隔一道"小裂缝"时，如何用渐进展开一点一点地走过去。

17.1 非简并微扰理论：渐进展开的基础

17.1.1 问题设定：已知与未知之间

假设我们有一个"简单"的哈密顿量 $\hat{H}_0$ ，其本征态 $|n^{(0)}\rangle$ 和本征值 $E_n^{(0)}$ 已经完全求解：

$\hat{H}_0 |n^{(0)}\rangle = E_n^{(0)} |n^{(0)}\rangle$

现在加上一个"小"微扰 $\hat{V}$ ，总哈密顿量为：

$\hat{H} = \hat{H}_0 + \hat{V}$

"小"的精确含义是：微扰引起的能级移动远小于未微扰能级之间的间距。

引入微扰参数 $\lambda$ （形式上，最后令 $\lambda = 1$ ）：

$\hat{H}(\lambda) = \hat{H}_0 + \lambda\hat{V}$

假设能量和波函数可以按 $\lambda$ 的幂次展开：

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

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

graph TD
    A[微扰展开] --> B["零阶: 已知解"]
    A --> C["一阶: 能量移动+波函数修正"]
    A --> D["二阶: 更精细的调整"]
    A --> E["..."]
    
    B --> B1["$E_n^{(0)}$"]
    B --> B2["$|n^{(0)}\rangle$"]
    
    C --> C1["$E_n^{(1)} = \langle n^{(0)}|\hat{V}|n^{(0)}\rangle$"]
    C --> C2["$|n^{(1)}\rangle = \sum_{"\\"m \neq n\\""}\frac{"\\"\langle m^{(0)\\"}|\hat{V}|n^{(0)}\rangle"}{E_n^{(0)}-E_m^{(0)}}|m^{(0)}\rangle$"]
    
    D --> D1["$E_n^{(2)} = \sum_{"\\"m \neq n\\""}\frac{"\\"|\langle m^{(0)\\"}|\hat{V}|n^{(0)}\rangle|^2"}{E_n^{(0)}-E_m^{(0)}}$"]

17.1.2 一阶修正的完整推导

将展开式代入薛定谔方程 $\hat{H}|n\rangle = E_n|n\rangle$ ，并按 $\lambda$ 的幂次收集项：

零阶（ $\lambda^0$ ）：

$\hat{H}_0 |n^{(0)}\rangle = E_n^{(0)} |n^{(0)}\rangle$

这就是未微扰方程，已经满足。

一阶（ $\lambda^1$ ）：

$\hat{H}_0 |n^{(1)}\rangle + \hat{V}|n^{(0)}\rangle = E_n^{(0)}|n^{(1)}\rangle + E_n^{(1)}|n^{(0)}\rangle$

整理：

$(\hat{H}_0 - E_n^{(0)})|n^{(1)}\rangle = (E_n^{(1)} - \hat{V})|n^{(0)}\rangle$

能量一阶修正：用 $\langle n^{(0)}|$ 左乘上式。左边为零（因为 $\hat{H}_0$ 厄米），右边给出：

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

这就是著名的一阶能量修正等于微扰在未微扰态中的期望值。

波函数一阶修正：将 $|n^{(1)}\rangle$ 用未微扰基展开：

$|n^{(1)}\rangle = \sum_{m \neq n} c_m^{(1)} |m^{(0)}\rangle$

代入一阶方程，用 $\langle m^{(0)}|$ 左乘：

$(E_m^{(0)} - E_n^{(0)})c_m^{(1)} = -\langle m^{(0)}|\hat{V}|n^{(0)}\rangle$

因此：

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

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

graph LR
    subgraph "一阶修正的物理"
    A[能量修正] --> B["$\langle n^{(0)}|\hat{V}|n^{(0)}\rangle$"]
    B --> C["微扰在\"旧态\"中的平均值"]
    
    D[波函数修正] --> E[混合其他态]
    E --> F["混合幅度 $\propto V_{mn}$"]
    E --> G["混合反比于能隙 $\Delta E$"]
    
    H[物理图像] --> I["微扰把\"旧态\"推向\"邻近态\""]
    I --> J[邻近越近、耦合越强]
    J --> K[混合越大]
    end

17.1.3 二阶能量修正的推导与物理意义

二阶能量修正是最重要的高阶修正：

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

推导：将一阶波函数修正代入二阶方程（ $\lambda^2$ 项），并用 $\langle n^{(0)}|$ 左乘。

二阶方程：

$(\hat{H}_0 - E_n^{(0)})|n^{(2)}\rangle = E_n^{(2)}|n^{(0)}\rangle + E_n^{(1)}|n^{(1)}\rangle - \hat{V}|n^{(1)}\rangle$

左乘 $\langle n^{(0)}|$ ，左边为零，得：

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

由于 $\langle n^{(0)}|n^{(1)}\rangle = 0$ （正交性条件），：

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

关键特征：

分母是能级差——微扰混合态 $m$ 对态 $n$ 的贡献反比于两者的能隙
如果 $E_m^{(0)} < E_n^{(0)}$ （ $m$ 在 $n$ 下方），分母为负，贡献使 $E_n$ 降低
如果 $E_m^{(0)} > E_n^{(0)}$ （ $m$ 在 $n$ 上方），分母为正，贡献使 $E_n$ 升高
总体效应：微扰使能级互相排斥——基态总是被压低，激发态可能被推高或压低（取决于周围能级分布）

对于基态（ $n=0$ ），所有 $E_m^{(0)} > E_0^{(0)}$ ，因此所有二阶项分母为负。结合分子 $|V_{m0}|^2 \geq 0$ ，得到：

$E_0^{(2)} \leq 0$

这意味着二阶修正总是降低基态能量（一阶修正可以为正或负）。

17.1.4 氢原子的斯塔克效应

斯塔克效应（Stark effect）是原子在外电场中的能级移动。对于氢原子，未微扰哈密顿量是库仑势，微扰是电场：

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

其中 $\mathcal{E}$ 是沿 $z$ 轴的电场强度。

基态（ $n=1$ ）的一阶修正：

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

因为基态波函数是球对称的， $z$ 是奇宇称，积分等于零。氢原子基态没有线性斯塔克效应。

二阶修正（非零）：

$E_{1s}^{(2)} = (e\mathcal{E})^2 \sum_{nlm \neq 1s} \frac{|\langle nlm|z|1s\rangle|^2}{E_{1s} - E_{nlm}}$

由于 $E_{1s} < E_{nlm}$ ，所有分母为负，所以 $E_{1s}^{(2)} < 0$ 。基态能量被电场压低，移动 $\propto \mathcal{E}^2$ 。

激发态（ $n=2$ ）的困境： $n=2$ 能级是四重简并的（ $2s, 2p_0, 2p_{\pm 1}$ ）。非简并微扰理论失效，因为分母中可能出现零。这引出了下一节的主题。

17.2 简并微扰理论：当简并遇到微扰

17.2.1 为什么简并需要特殊处理？

当未微扰能级存在简并（多个态共享同一能量）时，微扰理论的分母 $E_n^{(0)} - E_m^{(0)}$ 可能为零，导致表达式发散。

这不是数学上的偶然，而是物理上的必然。考虑一个 $g$ 重简并的子空间，由态 $\{|n_\alpha^{(0)}\rangle : \alpha = 1, \ldots, g\}$ 张成。微扰 $\hat{V}$ 在这个子空间内不一定对角——它会把简并态混合起来。

核心洞察：在简并子空间内，正确的"零阶近似态"不是任意的简并态，而是让微扰对角化的那组态。这些态称为好量子数态（good quantum states）。

17.2.2 简并微扰的步骤

步骤1：在简并子空间中，计算微扰矩阵

$V_{\alpha\beta} = \langle n_\alpha^{(0)}|\hat{V}|n_\beta^{(0)}\rangle$

步骤2：对角化这个 $g \times g$ 矩阵，得到本征值 $E_n^{(1),k}$ （ $k = 1, \ldots, g$ ）和本征矢 $|n^{(0),k}\rangle$ 。

步骤3：这些本征值就是一阶能量修正，本征矢就是正确的零阶波函数。简并通常被部分或完全解除。

graph TD
    A[简并微扰流程] --> B[识别简并子空间]
    B --> C["计算微扰矩阵 $V_{"\\"\alpha\beta\\""}$"]
    C --> D[对角化微扰矩阵]
    D --> E["本征值: 一阶能量修正"]
    D --> F["本征矢: 正确的零阶态"]
    E --> G[简并被部分解除]
    
    H["例子: $n=2$氢原子"] --> I[4维简并子空间]
    I --> J["电场微扰 $V = e\mathcal{E}z$"]
    J --> K[矩阵对角化]
    K --> L[线性斯塔克效应]
    L --> M["能级分裂 $\propto \mathcal{E}$"]

17.2.3 氢原子n=2的斯塔克效应

在 $n=2$ 子空间中，四个态是： $|2s\rangle$ （ $l=0, m=0$ ）和 $|2p_0\rangle, |2p_{+1}\rangle, |2p_{-1}\rangle$ （ $l=1, m=0, \pm 1$ ）。

微扰 $V = e\mathcal{E}z = e\mathcal{E}r\cos\theta$ 的选择定则：

改变 $l$ 的奇偶性（因为 $z$ 是奇宇称）
$\Delta m = 0$ （因为 $z$ 与 $\phi$ 无关）

因此， $\langle 2s|z|2p_0\rangle \neq 0$ ，但 $\langle 2s|z|2p_{\pm 1}\rangle = 0$ ， $\langle 2p_{\pm 1}|z|2p_{\pm 1}\rangle = 0$ 。

在 $(|2s\rangle, |2p_0\rangle, |2p_{+1}\rangle, |2p_{-1}\rangle)$ 基下，微扰矩阵为：

$\hat{V} = \begin{pmatrix} 0 & v & 0 & 0 \\ v & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix}$

其中 $v = e\mathcal{E}\langle 2s|z|2p_0\rangle$ 。

对角化 $2 \times 2$ 块：本征值为 $\pm v$ ，对应态 $\frac{1}{\sqrt{2}}(|2s\rangle \pm |2p_0\rangle)$ 。

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

$E_2^{(0)} + |v|$ （一个态）
$E_2^{(0)} - |v|$ （一个态）
$E_2^{(0)}$ （两个简并态， $|2p_{\pm 1}\rangle$ ）

线性斯塔克效应：能量移动 $\propto \mathcal{E}$ （一阶效应），而非非简并情况的 $\propto \mathcal{E}^2$ 。

graph TD
    A["$n=2$ 未微扰: 简并"] --> B[加上电场]
    B --> C["$|2p_{"\\"\pm 1\\""}\rangle$ 不受影响"]
    B --> D["$|2s\rangle, |2p_0\rangle$ 混合"]
    D --> E["新态: $\frac{1}{\sqrt{2}}(|2s\rangle \pm |2p_0\rangle)$"]
    E --> F["能量分裂: $\Delta E = \pm 3ea_0\mathcal{E}$"]
    
    G[物理图像] --> H[混合态有永久电偶极矩]
    H --> I[电场直接耦合]
    I --> J[线性移动]

17.2.4 范德瓦尔斯相互作用

范德瓦尔斯力（van der Waals force）是中性的原子或分子之间的微弱吸引力。它的起源可以用简并微扰理论理解。

考虑两个氢原子，相距 $R$ （ $R \gg a_0$ ）。它们之间没有净电荷，但每个原子的电子云产生瞬时电偶极矩。一个原子的偶极矩在另一个原子处产生电场，诱导出响应。

总哈密顿量可以写成两个氢原子哈密顿量加上相互作用：

$\hat{H} = \hat{H}_A + \hat{H}_B + \hat{V}_{int}$

相互作用在远距离展开的主导项是偶极-偶极相互作用：

$\hat{V}_{int} \approx \frac{1}{4\pi\varepsilon_0 R^3}\left(\mathbf{d}_A \cdot \mathbf{d}_B - 3(\mathbf{d}_A \cdot \hat{R})(\mathbf{d}_B \cdot \hat{R})\right)$

其中 $\mathbf{d}_A = -e\mathbf{r}_A$ 是原子A的偶极矩算符。

未微扰基态是两个氢原子基态的乘积： $|1s\rangle_A \otimes |1s\rangle_B$ ，能量 $E^{(0)} = 2E_{1s}$ 。

一阶修正： $E^{(1)} = \langle 1s, 1s|\hat{V}_{int}|1s, 1s\rangle = 0$ （因为 $\langle 1s|\mathbf{d}|1s\rangle = 0$ ，基态没有永久偶极矩）。

二阶修正（用非简并微扰，因为基态非简并）：

$E^{(2)} = \sum_{n_A, n_B \neq 1s} \frac{|\langle n_A, n_B|\hat{V}_{int}|1s, 1s\rangle|^2}{2E_{1s} - E_{n_A} - E_{n_B}}$

对所有激发态求和非常复杂，但量纲分析给出：

$E^{(2)} \sim -\frac{C_6}{R^6}$

其中 $C_6$ 是范德瓦尔斯系数，对于氢原子 $C_6 \approx 6.5 e^2 a_0^5/(4\pi\varepsilon_0)$ 。

物理图像：虽然每个原子基态的平均偶极矩为零，但量子涨落导致瞬时关联——一个原子的偶极矩涨落诱导另一个原子的响应，产生净的吸引势。这就是色散力或伦敦力的量子起源。

graph TD
    A[两个中性氢原子] --> B["远距离 $R \gg a_0$"]
    B --> C[无永久偶极矩]
    C --> D["一阶修正 = 0"]
    
    E[量子涨落] --> F[瞬时偶极矩]
    F --> G["偶极-偶极相互作用"]
    G --> H[二阶微扰]
    H --> I["$E^{(2)} = -C_6/R^6$"]
    
    I --> J[范德瓦尔斯吸引]
    J --> K[惰性气体液化]
    J --> L[壁虎爬墙]

前置知识：电磁场的半经典处理

在进入含时微扰理论和光与物质相互作用之前，我们需要澄清电磁场的处理方式——半经典近似。

B.1 三种处理层次

第一层次：完全经典

原子：量子力学（波函数、能级）
电磁场：经典麦克斯韦理论（连续波）

这是半经典处理。它足以解释吸收、受激发射和大多数光谱现象。但不能解释自发辐射。

第二层次：场量子化（单光子层面）

原子：量子力学
电磁场：量子化（光子、产生/湮灭算符）

这可以解释自发辐射，但数学复杂。

第三层次：完全量子电动力学

原子+场：统一量子场论
包含真空涨落、Lamb位移、反常磁矩等

B.2 为什么半经典处理在原子物理中足够？

原子尺度（ $a_0 \sim 0.5$ Å）远小于光波长（可见光 $\lambda \sim 5000$ Å）。因此，在原子范围内，电场几乎是空间均匀的：

$\mathbf{E}(\mathbf{r}, t) \approx \mathbf{E}(0, t) = \mathbf{E}_0\cos(\omega t)$

这就是电偶极近似。在这个近似下，电磁场不需要量子化就能给出吸收和受激发射的正确结果。

自发辐射的"缺口"：半经典理论中，如果外场为零（ $\mathbf{E}_0 = 0$ ），原子不会辐射。但实验中激发态原子确实会自发辐射。这个缺口需要用场的量子化来填补——即使在"真空"中，也存在电磁场的零点涨落，激发态原子与这些涨落耦合而辐射。

B.3 爱因斯坦A/B系数的半经典推导

虽然半经典理论不能自发导出A系数，但我们可以用热力学论证（爱因斯坦1917年）建立A、B系数的关系。

考虑一个二能级系统（ $E_1 < E_2$ ），在温度 $T$ 的热平衡中：

吸收速率： $N_1 B_{12} \rho(\omega)$
受激发射速率： $N_2 B_{21} \rho(\omega)$
自发辐射速率： $N_2 A_{21}$

热平衡要求吸收=发射：

$N_1 B_{12} \rho(\omega) = N_2 B_{21} \rho(\omega) + N_2 A_{21}$

利用玻尔兹曼分布 $N_2/N_1 = e^{-\hbar\omega/k_BT}$ 和普朗克黑体辐射公式 $\rho(\omega) = \frac{\hbar\omega^3}{\pi^2 c^3}\frac{1}{e^{\hbar\omega/k_BT}-1}$ ，可以导出：

$B_{12} = B_{21}$

$A_{21} = \frac{\hbar\omega^3}{\pi^2 c^3} B_{21}$

量子场论进一步给出：

$\Gamma_{spont} = A_{21} = \frac{\omega^3 |\mathbf{d}_{21}|^2}{3\pi\varepsilon_0\hbar c^3}$

graph TD
    A[电磁场处理层次] --> B[完全经典]
    A --> C[半经典]
    A --> D[场量子化]
    
    B --> B1[麦克斯韦方程]
    B --> B2[不能解释自发辐射]
    
    C --> C1["经典场+量子原子"]
    C --> C2["解释吸收/受激发射"]
    C --> C3[需要热力学引入A系数]
    
    D --> D1["光子产生/湮灭"]
    D --> D2[自动给出自发辐射]
    D --> D3["Lamb位移/反常磁矩"]
    
    C -.->|对于大多数光谱| D

18.1 含时微扰理论：时间驱动的跃迁

18.1.1 从定态到含时问题

前面的微扰理论处理的是定态问题——微扰不随时间变化，我们寻找修正后的能量本征态。但许多物理过程涉及含时微扰：电磁波照射原子、磁场随时间变化、粒子碰撞等。

含时薛定谔方程：

$i\hbar\frac{d}{dt}|\psi(t)\rangle = \left(\hat{H}_0 + \hat{V}(t)\right)|\psi(t)\rangle$

假设系统初始处于 $\hat{H}_0$ 的本征态 $|i\rangle$ （ $\hat{H}_0|i\rangle = E_i|i\rangle$ ）。问题是：经过时间 $t$ ，系统跃迁到另一个本征态 $|f\rangle$ 的概率是多少？

18.1.2 相互作用绘景与跃迁振幅

使用相互作用绘景（interaction picture）：

$|\psi_I(t)\rangle = e^{i\hat{H}_0 t/\hbar}|\psi(t)\rangle$

$\hat{V}_I(t) = e^{i\hat{H}_0 t/\hbar}\hat{V}(t)e^{-i\hat{H}_0 t/\hbar}$

薛定谔方程变为：

$i\hbar\frac{d}{dt}|\psi_I(t)\rangle = \hat{V}_I(t)|\psi_I(t)\rangle$

形式上积分：

$|\psi_I(t)\rangle = |i\rangle + \frac{1}{i\hbar}\int_0^t \hat{V}_I(t')|\psi_I(t')\rangle dt'$

一阶近似（将 $|\psi_I(t')\rangle \approx |i\rangle$ 代入右边）：

$|\psi_I^{(1)}(t)\rangle = |i\rangle + \frac{1}{i\hbar}\int_0^t \hat{V}_I(t')|i\rangle dt'$

跃迁到 $|f\rangle$ 的振幅：

$c_f^{(1)}(t) = \langle f|\psi_I^{(1)}(t)\rangle = \frac{1}{i\hbar}\int_0^t \langle f|\hat{V}_I(t')|i\rangle dt'$

$= \frac{1}{i\hbar}\int_0^t e^{i(E_f-E_i)t'/\hbar}\langle f|\hat{V}(t')|i\rangle dt'$

graph TD
    A[含时微扰] --> B[相互作用绘景]
    B --> C["$\hat{V}_I(t) = e^{iH_0t/\hbar}\hat{V}(t)e^{-iH_0t/\hbar}$"]
    C --> D[一阶跃迁振幅]
    D --> E["$c_f^{(1)}(t) = \frac{1}{i\hbar}\int_0^t e^{i\omega_{fi}t'}V_{fi}(t')dt'$"]
    
    E --> F["$\omega_{fi} = (E_f-E_i)/\hbar$"]
    F --> G[玻尔频率条件]
    G --> H[跃迁需要能量匹配]

18.1.3 正弦微扰与费米黄金规则

最常见的含时微扰是正弦（或余弦）形式，如单色光场：

$\hat{V}(t) = \hat{V}_0 \cos(\omega t) = \frac{\hat{V}_0}{2}(e^{i\omega t} + e^{-i\omega t})$

跃迁振幅：

$c_f^{(1)}(t) = -\frac{V_{fi}}{2\hbar}\left(\frac{e^{i(\omega_{fi}+\omega)t} - 1}{\omega_{fi}+\omega} + \frac{e^{i(\omega_{fi}-\omega)t} - 1}{\omega_{fi}-\omega}\right)$

其中 $V_{fi} = \langle f|\hat{V}_0|i\rangle$ 。

共振条件：当 $\omega \approx \omega_{fi}$ （吸收）或 $\omega \approx -\omega_{fi}$ （受激发射）时，其中一个分母很小，该项主导。

费米黄金规则（Fermi’s Golden Rule）：对于长时间和连续末态，跃迁速率为：

$\Gamma_{i \to f} = \frac{2\pi}{\hbar}|\langle f|\hat{V}_0|i\rangle|^2 \rho(E_f)$

其中 $\rho(E_f)$ 是末态的态密度（单位能量的态数）。

graph TD
    A[正弦微扰] --> B["$\hat{V}(t) = \hat{V}_0\cos("\\"\omega t\\"")$"]
    B --> C["两项: $\omega_{fi} \pm \omega$"]
    C --> D["$\omega \approx \omega_{fi}$: 吸收主导"]
    C --> E["$\omega \approx -\omega_{fi}$: 发射主导"]
    
    D --> F[长时间极限]
    E --> F
    F --> G[费米黄金规则]
    G --> H["$\Gamma = \frac{"\\"2\pi\\""}{\hbar}|V_{fi}|^2\rho(E_f)$"]
    
    H --> I[跃迁速率正比于矩阵元平方]
    H --> J[正比于末态态密度]

18.1.4 光与物质的相互作用

原子与电磁场的相互作用是含时微扰理论最重要的应用。在电偶极近似下（波长远大于原子尺寸），相互作用哈密顿量为：

$\hat{V} = -\hat{\mathbf{d}}\cdot\mathbf{E}(t) = e\hat{\mathbf{r}}\cdot\mathbf{E}_0\cos(\omega t)$

其中 $\hat{\mathbf{d}} = -e\hat{\mathbf{r}}$ 是电偶极矩算符。

跃迁矩阵元 $V_{fi} = e\mathbf{E}_0\cdot\langle f|\hat{\mathbf{r}}|i\rangle$ 决定了选择定则。

吸收（原子从低能态到高能态，吸收一个光子）：

$\Gamma_{abs} = \frac{\pi e^2}{\hbar^2} |\mathbf{E}_0\cdot\langle f|\hat{\mathbf{r}}|i\rangle|^2 \delta(\omega - \omega_{fi})$

受激发射（原子从高能态到低能态，在辐射场刺激下发射一个光子）：

$\Gamma_{stim} = \frac{\pi e^2}{\hbar^2} |\mathbf{E}_0\cdot\langle i|\hat{\mathbf{r}}|f\rangle|^2 \delta(\omega - \omega_{fi})$

注意 $\Gamma_{abs} = \Gamma_{stim}$ （因为 $\langle f|\hat{\mathbf{r}}|i\rangle = \langle i|\hat{\mathbf{r}}|f\rangle^*$ ）。

graph TD
    A[光与物质相互作用] --> B[电偶极近似]
    B --> C["$\hat{V} = -\hat{"\\"\mathbf{d\\"}"}\cdot\mathbf{E}$"]
    
    C --> D[吸收]
    C --> E[受激发射]
    C --> F[自发辐射]
    
    D --> G["原子+光子 $\to$ 激发原子"]
    E --> H["激发原子+光子 $\to$ 原子+2光子"]
    F --> I["激发原子 $\to$ 原子+光子"]
    
    G --> J[需要外场光子]
    H --> J
    I --> K[无外场也发生]
    K --> L[爱因斯坦A系数]

18.1.5 自发辐射：真空不是空的

自发辐射（spontaneous emission）是激发态原子在没有外场的情况下自发发射光子的过程。这不能直接从经典微扰理论解释——经典理论中，如果没有外场驱动，原子不会辐射。

爱因斯坦在1917年用热力学论证证明了自发辐射的存在，并推导出它与受激发射的关系：

$A_{if} = \frac{\hbar\omega_{if}^3}{\pi^2 c^3} B_{if}$

其中 $A_{if}$ 是自发辐射系数， $B_{if}$ 是受激系数。

量子场论解释：自发辐射的本质是电磁场的量子化。即使在真空中，电磁场也有量子涨落（真空涨落）。激发态原子与这些真空涨落相互作用，"被真空刺激"而辐射。

定量地说，自发辐射速率等于将电磁场量子化后，用费米黄金规则计算真空场引起的跃迁速率：

$\Gamma_{spont} = \frac{\omega_{if}^3|\mathbf{d}_{if}|^2}{3\pi\varepsilon_0\hbar c^3}$

graph TD
    A[自发辐射] --> B["经典理论: 无外场不辐射"]
    B --> C["矛盾!"]
    
    D[量子解释] --> E[电磁场量子化]
    E --> F[真空有零点涨落]
    F --> G[激发态与真空涨落耦合]
    G --> H[自发辐射]
    
    H --> I[爱因斯坦A系数]
    I --> J["$A_{if} \propto \omega^3 |\mathbf{d}_{if}|^2$"]

18.1.6 激光原理：受激发射的放大

激光（LASER = Light Amplification by Stimulated Emission of Radiation）的原理建立在受激发射与吸收的竞争中。

考虑一个二能级系统，能量 $E_1 < E_2$ 。在辐射场中，三个过程同时发生：

吸收（ $1 \to 2$ ）：速率 $\propto N_1 B_{12} \rho(\omega)$
受激发射（ $2 \to 1$ ）：速率 $\propto N_2 B_{21} \rho(\omega)$
自发辐射（ $2 \to 1$ ）：速率 $\propto N_2 A_{21}$

其中 $N_1, N_2$ 是两能级的原子数， $\rho(\omega)$ 是辐射场能量密度。

热平衡时（玻尔兹曼分布 $N_2/N_1 = e^{-(E_2-E_1)/k_BT} < 1$ ），吸收多于受激发射，光被净吸收。

粒子数反转（population inversion）：如果通过某种机制（光泵浦、电注入、化学能等）使 $N_2 > N_1$ ，则受激发射超过吸收。入射光被放大而非吸收。

graph TD
    A[激光原理] --> B[二能级系统]
    B --> C["吸收: $N_1 \to N_2$"]
    B --> D["受激发射: $N_2 \to N_1$"]
    B --> E["自发辐射: $N_2 \to N_1$"]
    
    F[热平衡] --> G["$N_2 < N_1$"]
    G --> H["吸收 > 受激发射"]
    H --> I[光被吸收]
    
    J[粒子数反转] --> K["$N_2 > N_1$"]
    K --> L["受激发射 > 吸收"]
    L --> M[光被放大]
    M --> N[激光输出]
    
    O[实现粒子数反转] --> P[三能级系统]
    O --> Q[四能级系统]
    O --> R[半导体注入]

激光的三个必要条件：

粒子数反转：实现 $N_2 > N_1$
谐振腔：光子在两个镜子之间来回反射，反复触发受激发射
增益介质：提供放大

受激发射的一个重要特性：相干性。受激发射产生的光子与"刺激"它的光子具有相同的频率、相位、偏振和传播方向。这就是为什么激光是相干光，而普通光源（如灯泡）是自发辐射的混合，不相干。

18.2 绝热近似与Berry相位

18.2.1 绝热定理：慢工出细活

绝热定理（Adiabatic Theorem）是含时量子力学中的一个基本结果。它说：如果一个系统的哈密顿量 $H(t)$ 随时间足够缓慢地变化，且能级之间没有简并，那么系统初始处于某个瞬时本征态 $|n(t)\rangle$ 的话，它将始终停留在该本征态（只获得一个相位因子）。

数学表述：设瞬时本征态满足 $\hat{H}(t)|n(t)\rangle = E_n(t)|n(t)\rangle$ 。如果系统初始处于 $|n(0)\rangle$ ，且演化足够慢：

$T \gg \frac{\hbar}{\min_{m \neq n}|E_n - E_m|}$

则：

$|\psi(t)\rangle \approx e^{i\theta_n(t)}e^{i\gamma_n(t)}|n(t)\rangle$

其中：

动力学相位： $\theta_n(t) = -\frac{1}{\hbar}\int_0^t E_n(t')dt'$
几何相位（Berry相位）： $\gamma_n(t) = i\int_0^t \langle n(t')|\frac{d}{dt'}|n(t')\rangle dt'$

graph TD
    A[绝热定理] --> B[哈密顿量缓慢变化]
    B --> C["$T \gg \hbar/\Delta E$"]
    C --> D[系统跟随瞬时本征态]
    
    D --> E["$|\psi(t)\rangle = e^{i\theta}e^{i\gamma}|n(t)\rangle$"]
    E --> F["$\theta$: 动力学相位"]
    E --> G["$\gamma$: Berry几何相位"]
    
    H[动力学相位] --> I[来自能量演化]
    I --> J["$\theta = -\int E(t)dt/\hbar$"]
    
    K[几何相位] --> L[来自参数空间路径]
    L --> M[与演化快慢无关]
    M --> N[可观测的量子效应]

18.2.2 Berry相位：参数空间的几何

Berry相位的深刻之处在于它的几何性。考虑哈密顿量依赖于一组外部参数 $\mathbf{R} = (R_1, R_2, \ldots)$ （可以是磁场方向、电场强度、晶格常数等），这些参数随时间变化： $\mathbf{R}(t)$ 。

Berry相位可以写成参数空间的路径积分：

$\gamma_n = i\oint \langle n(\mathbf{R})|\nabla_{\mathbf{R}}|n(\mathbf{R})\rangle \cdot d\mathbf{R}$

定义 Berry联络（Berry connection）：

$\mathbf{A}_n(\mathbf{R}) = i\langle n(\mathbf{R})|\nabla_{\mathbf{R}}|n(\mathbf{R})\rangle$

Berry相位就是Berry联络沿闭合路径的线积分。

关键洞察：如果参数沿着一个闭合回路 $C$ 演化，回到初始值，Berry相位不为零。这类似于电磁学中的 Aharonov-Bohm效应——即使磁场为零，矢势的回路积分也可以产生可观测相位。

用Stokes定理，线积分可以转化为面积分：

$\gamma_n = \oint_C \mathbf{A}_n \cdot d\mathbf{R} = \int_S (\nabla \times \mathbf{A}_n) \cdot d\mathbf{S} = \int_S \mathbf{\Omega}_n \cdot d\mathbf{S}$

其中 Berry曲率（Berry curvature）：

$\mathbf{\Omega}_n = \nabla \times \mathbf{A}_n$

graph TD
    A["参数空间 $R$"] --> B["闭合路径 $C$"]
    B --> C["Berry相位 $\gamma = \oint \mathbf{A}\cdot d\mathbf{R}$"]
    C --> D[Stokes定理]
    D --> E["$\gamma = \int_S \mathbf{"\\"\Omega\\""}\cdot d\mathbf{S}$"]
    
    E --> F["Berry曲率 $\mathbf{"\\"\Omega\\""} = \nabla\times\mathbf{A}$"]
    F --> G["类比: 磁通量"]
    G --> H["参数空间的\"磁场\""]
    
    I[物理意义] --> J[与路径形状有关]
    J --> K[与演化快慢无关]
    K --> L[拓扑不变量]
    L --> M[量子霍尔效应]

18.2.3 自旋1/2在旋转磁场中的Berry相位

最经典的Berry相位例子：一个自旋- $1/2$ 粒子在缓慢旋转的磁场中。

磁场 $\mathbf{B}(t) = B(\sin\theta\cos\phi(t), \sin\theta\sin\phi(t), \cos\theta)$ 以固定极角 $\theta$ 绕 $z$ 轴旋转，方位角 $\phi$ 从 $0$ 到 $2\pi$ 。

瞬时本征态（ $\hat{\mathbf{S}}\cdot\hat{\mathbf{B}}$ 的本征态）为自旋沿 $\mathbf{B}$ 方向的态 $|\uparrow_{\mathbf{B}}\rangle$ 和 $|\downarrow_{\mathbf{B}}\rangle$ 。

Berry相位计算：

$\gamma_{\uparrow} = \oint i\langle\uparrow_{\mathbf{B}}|\frac{\partial}{\partial\phi}|\uparrow_{\mathbf{B}}\rangle d\phi = \pi(1 - \cos\theta)$

这恰好等于磁场矢量扫过的立体角的一半（符号取决于自旋方向）！

$\gamma = \mp \frac{1}{2}\Omega$

其中 $\Omega = 2\pi(1 - \cos\theta)$ 是磁场方向在参数空间（球面）上扫过的立体角。

graph TD
    A[旋转磁场] --> B["$\mathbf{B}$ 绕 $z$ 轴旋转"]
    B --> C["固定极角 $\theta$"]
    C --> D["方位角 $\phi: 0 \to 2\pi$"]
    
    D --> E[参数空间路径]
    E --> F[球面上的纬线]
    F --> G["包围立体角 $\Omega$"]
    
    G --> H[Berry相位]
    H --> I["$\gamma = \mp \frac{1}{2}\Omega$"]
    I --> J["几何而非动力学!"]
    
    K[实验验证] --> L[自旋回波实验]
    L --> M[Berry相位可观测]
    M --> N[量子干涉]

Berry相位自1984年被Michael Berry系统阐述以来，已成为凝聚态物理中理解量子霍尔效应、拓扑绝缘体和反常霍尔效应的核心概念。

数值例子

例子1：氢原子2p→1s跃迁的费米黄金定则计算

例题：计算氢原子中电子从 $2p$ 态自发跃迁到 $1s$ 态的速率（爱因斯坦A系数）。

步骤1：确定跃迁参数

$E_{2p} - E_{1s} = 10.2$ eV $= \hbar\omega$
$\omega = 10.2 \times 1.602\times 10^{-19} / 1.055\times 10^{-34} = 1.55\times 10^{16}$ rad/s
对应波长 $\lambda = 2\pi c/\omega = 2\pi \times 3\times 10^8 / 1.55\times 10^{16} \approx 121.6$ nm（Lyman-α线）

步骤2：计算电偶极矩阵元

$\mathbf{d}_{if} = -e\langle 1s|\mathbf{r}|2p\rangle$

氢原子波函数：

$\psi_{1s} = \frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}$
$\psi_{2p, m=0} = \frac{1}{\sqrt{32\pi a_0^3}}\frac{r}{a_0}e^{-r/(2a_0)}\cos\theta$ （取 $z$ 方向）

计算 $\langle 1s|z|2p_0\rangle$ ：

$\langle 1s|z|2p_0\rangle = \int \psi_{1s}^*(r) r\cos\theta \psi_{2p_0}(r) d^3r$

角向积分： $\int_0^{\pi}\cos^2\theta\sin\theta d\theta = 2/3$ 。

径向积分：

$\int_0^{\infty} \frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0} \cdot r \cdot \frac{1}{\sqrt{32\pi a_0^3}}\frac{r}{a_0}e^{-r/(2a_0)} \cdot r^2 dr$

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

利用 $\int_0^{\infty} r^n e^{-\alpha r}dr = n!/\alpha^{n+1}$ ：

$= \frac{1}{\sqrt{32}\pi a_0^4} \cdot \frac{4!}{(3/2a_0)^5} = \frac{24}{\sqrt{32}\pi a_0^4} \cdot \frac{32 a_0^5}{243} = \frac{768}{243\sqrt{32}\pi}a_0$

$= \frac{512}{81\sqrt{2}\pi}a_0$

综合角向和径向：

$\langle 1s|z|2p_0\rangle = \sqrt{\frac{2}{3}} \cdot \frac{512}{81\sqrt{2}\pi}a_0 \cdot (2\pi) = \frac{512\sqrt{2}}{81\sqrt{3}}a_0$

等等，更仔细地计算。使用标准氢原子结果：

$|\langle 1s|\mathbf{r}|2p\rangle|^2 = \frac{2^{15}}{3^{10}}a_0^2 = \frac{32768}{59049}a_0^2 \approx 0.555 a_0^2$

步骤3：代入自发辐射公式

$\Gamma_{spont} = A_{21} = \frac{\omega^3 |\mathbf{d}_{21}|^2}{3\pi\varepsilon_0\hbar c^3}$

$= \frac{(1.55\times 10^{16})^3 \times e^2 \times 0.555 a_0^2}{3\pi\varepsilon_0\hbar c^3}$

使用 $e^2/(4\pi\varepsilon_0) = \alpha\hbar c \approx 1.44$ MeV·fm， $a_0 = 0.529$ Å $= 5.29\times 10^{-11}$ m。

$|\mathbf{d}|^2 = e^2 \times 0.555 \times (5.29\times 10^{-11})^2 = 0.555 e^2 \times 2.80\times 10^{-21}$

= 0.555 \times 4\pi\varepsilon_0 \times 1.44\times 10^{-15} \times 2.80\times 10^{-21} / (4\pi\varepsilon_0) \text{ ...更简单地用原子单位}

在原子单位下（ $\hbar = m_e = e = 4\pi\varepsilon_0 = 1$ ）：

$A_{21} = \frac{4\alpha\omega^3}{3c^2}|\langle f|\mathbf{r}|i\rangle|^2$

其中 $\alpha = 1/137$ ， $c = 137$ （原子单位）， $\omega = 10.2/27.2 \times 2 \approx 0.75$ （原子单位， $E_h/\hbar$ ）。

A_{21} = \frac{4}{3\times 137^3} \times (0.75)^3 \times 0.555 \times \frac{1}{\text{原子时间}}

原子时间 $= \hbar/E_h = 1.055\times 10^{-34}/(4.36\times 10^{-18}) = 2.42\times 10^{-17}$ s。

$A_{21} \approx \frac{4}{3\times 2.58\times 10^6} \times 0.422 \times 0.555 \times \frac{1}{2.42\times 10^{-17}}$

$\approx 6.25\times 10^8 \text{ s}^{-1}$

步骤4：计算寿命

$\tau = \frac{1}{A_{21}} \approx 1.6 \times 10^{-9} \text{ s} = 1.6 \text{ ns}$

实验测量值： $\tau(2p) \approx 1.6$ ns。理论与实验惊人一致！

这说明：

量子电动力学的微扰计算在原子物理中极为精确
电偶极近似（忽略磁场分量和多极展开的高阶项）对允许跃迁是极好的近似
自发辐射速率与 $\omega^3$ 成正比——高频跃迁（如X射线区）的寿命远短于可见光区

例子2：两个氢原子的范德瓦尔斯相互作用能

例题：计算两个基态氢原子在距离 $R = 10 a_0$ （约5 Å）时的范德瓦尔斯相互作用能。

步骤1：写出相互作用哈密顿量

$\hat{V}_{int} = \frac{e^2}{4\pi\varepsilon_0 R^3}(\mathbf{d}_A \cdot \mathbf{d}_B - 3d_{Az}d_{Bz})$

（取两原子连线为 $z$ 轴）

展开为：

$\hat{V}_{int} = \frac{e^2}{4\pi\varepsilon_0 R^3}(d_{Ax}d_{Bx} + d_{Ay}d_{By} - 2d_{Az}d_{Bz})$

步骤2：计算二阶微扰

$E^{(2)} = \sum_{n_A, n_B \neq 1s} \frac{|\langle n_A, n_B|\hat{V}_{int}|1s, 1s\rangle|^2}{2E_{1s} - E_{n_A} - E_{n_B}}$

由于氢原子的选择定则，只有 $n_A = n_B = 2$ （主要是 $2p$ 态）的贡献最大。

利用： $\langle 1s|d_z|2p_0\rangle = \frac{2^7\sqrt{2}}{3^5}ea_0 \approx 0.74 ea_0$ ， $\langle 1s|d_x|2p_x\rangle = \langle 1s|d_y|2p_y\rangle$ 类似。

步骤3：估算求和

对 $2p$ 态的主要贡献：

$E^{(2)}(2p) \approx \frac{1}{\Delta E}\left|\frac{e^2}{4\pi\varepsilon_0 R^3}\right|^2 |\langle 1s|\mathbf{d}|2p\rangle|^4$

其中 $\Delta E = E_{2p} - E_{1s} = 10.2$ eV（能级差为负，因为分母是 $2E_{1s} - E_{2p} - E_{2p} = -2\times 10.2$ eV）。

$E^{(2)} \approx -\frac{1}{20.4 \text{ eV}}\left(\frac{e^2}{4\pi\varepsilon_0 R^3}\right)^2 (0.74 ea_0)^4$

用原子单位：

$E^{(2)} \approx -\frac{1}{2\times 0.75 E_h}\left(\frac{1}{R^3}\right)^2 (0.74)^4 a_0^6 E_h^2$

$= -\frac{(0.74)^4}{1.5 R^6} E_h$

对于 $R = 10 a_0$ ：

$E^{(2)} = -\frac{0.30}{1.5 \times 10^6} E_h = -\frac{0.20}{10^6} \times 27.2 \text{ eV} = -5.4\times 10^{-6} \text{ eV}$

步骤4：与精确值比较

精确计算的范德瓦尔斯系数： $C_6 = 6.499 e^2 a_0^5/(4\pi\varepsilon_0) \approx 6.5$ （原子单位）。

$E_{vdW} = -\frac{C_6}{R^6} = -\frac{6.5}{10^6} E_h = -1.77\times 10^{-4} E_h \approx -4.8\times 10^{-6} \text{ eV}$

我们的粗略估算（仅考虑 $2p$ 态）给出 $-5.4\times 10^{-6}$ eV，与精确值同数量级。更精确的计算需要包含所有激发态的贡献。

转换为温度单位： $1$ eV $= 11600$ K，所以：

$E_{vdW} \approx -4.8\times 10^{-6} \times 11600 \text{ K} \approx -0.056 \text{ K}$

在室温（ $T \sim 300$ K）下，这个能量远小于热运动能——范德瓦尔斯力在室温下是微弱的。但在低温（液氦温度 $4.2$ K）下，它成为主导相互作用，解释了为什么氢气在 $20$ K左右液化。

本章总结

graph TD
    A["第17-18章: 微扰理论"] --> B[定态微扰]
    A --> C[含时微扰]
    A --> D[绝热近似与Berry相位]
    
    B --> B1["非简并: 一阶+二阶"]
    B --> B2["简并: 对角化微扰矩阵"]
    B --> B3[斯塔克效应]
    B --> B4[范德瓦尔斯力]
    
    C --> C1[正弦微扰]
    C --> C2[费米黄金规则]
    C --> C3["光与物质: 吸收/发射"]
    C --> C4[激光原理]
    
    D --> D1[绝热定理]
    D --> D2[Berry相位]
    D --> D3[几何相位 vs 动力学相位]
    D --> D4[拓扑量子现象]
    
    B1 -.->|基础| B2
    B2 -.->|应用| B3
    C1 -.->|导出| C2
    C2 -.->|应用| C3
    C3 -.->|竞争| C4
    D1 -.->|产生| D2
    D2 -.->|解释| D4

微扰理论是量子力学的"瑞士军刀"——它几乎适用于任何"足够接近可解"的系统。从原子光谱到分子间力，从激光到量子拓扑，渐进展开的思想贯穿始终。而Berry相位则是这场渐进展开中最优雅的意外发现：即使近似过程本身，也留下了几何的印记。

练习与思考

1. 微扰级数的收敛性

考虑一个一维谐振子加上微扰 $\hat{V} = \lambda x^4$ 。写出基态能量的一阶修正表达式（不需要算出积分）。如果 $\lambda$ 很大，微扰级数还能收敛吗？微扰理论的本质是什么类型的展开（泰勒展开、渐近展开、还是其他）？

提示：微扰级数通常是渐近级数（asymptotic series），而非收敛级数。这意味着最优截断点在有限阶，之后误差反而增大。

2. 选择定则与对称性

氢原子基态 $|1s\rangle$ 是偶宇称（ $l=0$ ），而电偶极算符 $z$ 是奇宇称。解释为什么 $\langle 1s|z|1s\rangle = 0$ 是宇称选择定则的结果。如果微扰是 $x^2 + y^2 + z^2 = r^2$ （偶宇称），一阶能量修正是否为零？

提示：宇称变换 $\hat{P}\psi(\mathbf{r}) = \psi(-\mathbf{r})$ 。计算 $\langle\psi|\hat{O}|\psi\rangle$ 在宇称下的行为。

3. Berry相位的可观测性

假设一个自旋- $1/2$ 粒子在磁场中绝热演化。为什么动力学相位 $e^{i\theta}$ 在干涉实验中通常被"抵消"，而Berry相位 $e^{i\gamma}$ 可以被观测？设计一个实验（基于自旋回波或干涉仪）来测量Berry相位，并解释为什么需要两条不同路径才能提取出几何相位。

提示：考虑让粒子沿两条不同路径演化，但具有相同的初始和末态。动力学相位相同（因为能量相同），Berry相位不同（因为路径不同）。干涉揭示相位差。

4. 费米黄金定则的推导

从含时微扰理论的一阶跃迁振幅出发，推导费米黄金规则：

$\Gamma_{i\to f} = \frac{2\pi}{\hbar}|V_{fi}|^2\rho(E_f)$

关键步骤：

写出 $|c_f^{(1)}(t)|^2$ 的表达式
证明长时间极限下跃迁概率正比于时间
定义单位时间跃迁速率 $\Gamma = |c_f|^2/t$
利用 $\lim_{t\to\infty}\frac{\sin^2(\omega t/2)}{\omega^2 t} = \frac{\pi}{2}\delta(\omega)$

5. 简并微扰的数值计算

一个二维谐振子， $\hat{H}_0 = \frac{p_x^2+p_y^2}{2m} + \frac{1}{2}m\omega^2(x^2+y^2)$ 。第一激发态是二重简并的： $|1,0\rangle$ 和 $|0,1\rangle$ （能量都是 $2\hbar\omega$ ）。加上微扰 $\hat{V} = \lambda xy$ 。

在 $(|1,0\rangle, |0,1\rangle)$ 基下写出微扰矩阵
对角化求能量一阶修正和正确的零阶态
解释为什么新的本征态是 $|+\rangle = \frac{1}{\sqrt{2}}(|1,0\rangle + |0,1\rangle)$ 和 $|-\rangle = \frac{1}{\sqrt{2}}(|1,0\rangle - |0,1\rangle)$

第17-18章 微扰理论：渐进展开的艺术