"当德布罗意波长远小于势能变化的尺度时，量子力学看起来像经典力学——但带有量子修正。"
WKB近似是连接经典世界与量子世界的桥梁。

第8章 WKB近似：半经典的量子航行

8.1 故事：量子航海家的导航图

24世纪，"量子航海家"驾驶着微型探测器穿越复杂的势场地形。他们的导航图上，有些区域标注为"经典海域"——波涛规律，可以用经典力学预测；有些区域是"量子迷雾"——波函数干涉严重，一切都不确定。

资深航海家老王指着前方的一片渐变浅滩："这里势能变化很缓慢，德布罗意波长相对于势的变化尺度很短。这种情况下，我们可以用WKB近似——把波函数想象成局部平面波，振幅和相位随位置缓慢变化。"

"像经典的粒子，但带着量子相位？"

"对。WKB是’半经典’的：粒子的运动轨迹像经典粒子，但波的相位累积产生量子效应——比如干涉和隧穿。"

年轻的副手小李问："那遇到势垒呢？"

"经典物理说’此路不通’。但WKB告诉我们，如果势垒不是无穷陡峭，波函数会在势垒内指数衰减，然后在另一侧重新出现——就像船只穿过薄雾，看似消失了，却在雾的另一端重新浮现。"

核心洞察：WKB近似在势能变化缓慢的区域有效，它将波函数表示为局部平面波，是理解隧穿和量子化条件的强大工具。

8.2 前置知识：渐近展开与半经典方法

8.2.1 什么是渐近展开？

在数学物理中，许多函数没有简单的闭式表达式，但可以在某个极限下展开为级数。渐近展开（Asymptotic Expansion）是一种特殊的级数展开：

$f(x) \sim \sum_{n=0}^{\infty} \frac{a_n}{x^n} \quad (x \to \infty)$

与泰勒级数不同，渐近级数通常发散——但它的前几项可以给出极好的近似。

例子：误差函数的渐近展开

$\text{erfc}(x) \sim \frac{e^{-x^2}}{x\sqrt{\pi}}\left(1 - \frac{1}{2x^2} + \frac{3}{4x^4} - \cdots\right)$

当 $x$ 很大时，只取第一项就能得到很好的近似。

8.2.2 半经典方法的思想

半经典方法（Semiclassical Method）的核心思想是：在普朗克常数 $\hbar$ "很小"的极限下（即作用量远大于 $\hbar$ ），量子力学应该接近经典力学。

形式化地，我们将波函数按 $\hbar$ 的幂次展开：

$\psi(x) = \exp\left(\frac{i}{\hbar}S(x)\right) = \exp\left(\frac{i}{\hbar}S_0(x) + iS_1(x) + i\hbar S_2(x) + \cdots\right)$

$S_0$ ：经典作用量（零阶，对应经典轨迹）
$S_1$ ：量子修正（一阶，给出振幅）
$S_2, ...$ ：更高阶量子修正

对应原理：当 $\hbar \to 0$ 时， $S_0$ 主导，量子力学退化为经典力学。

8.2.3 WKB近似的数学基础

WKB近似本质上是快速振荡积分的驻相近似在微分方程中的应用。当相位快速变化时，只有相位稳定的点附近对积分有贡献——这与经典轨迹的概念密切相关。

8.3 WKB波函数的严格推导

8.3.1 基本假设

假设势能 $V(x)$ 变化缓慢，即在一个德布罗意波长 $\lambda = h/p$ 内，势能变化很小：

$\lambda \left|\frac{dV}{dx}\right| \ll |E - V(x)|$

或者等价地：

$\frac{d\lambda}{dx} \ll 1$

这意味着波函数在局部看起来像平面波，但振幅和波数随位置缓慢变化。

8.3.2 WKB波函数的推导

设 $\psi(x) = A(x)e^{i\phi(x)}$ ，代入定态薛定谔方程：

$-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} + V(x)\psi = E\psi$

令 $\psi(x) = \exp(iS(x)/\hbar)$ ，其中 $S(x)$ 是复函数。代入后得到：

$i\hbar S'' - (S')^2 + p^2(x) = 0$

其中 $p(x) = \sqrt{2m(E-V(x))}$ 是局部动量。

将 $S(x)$ 按 $\hbar$ 展开： $S = S_0 + \hbar S_1 + \hbar^2 S_2 + ...$

零阶（ $\hbar^0$ ）：

$-(S_0')^2 + p^2 = 0 \Rightarrow S_0' = \pm p(x)$

$S_0 = \pm \int p(x) dx$

一阶（ $\hbar^1$ ）：

$iS_0'' - 2S_0'S_1' = 0 \Rightarrow S_1' = \frac{iS_0''}{2S_0'} = \frac{ip'}{2p}$

$S_1 = \frac{i}{2}\ln p + C$

因此：

$e^{iS_1} = \frac{C'}{\sqrt{p(x)}}$

WKB波函数：

经典区域（ $E > V(x)$ ，动量为实数）：

$\psi(x) \approx \frac{C}{\sqrt{p(x)}} \exp\left(\pm \frac{i}{\hbar}\int p(x) dx\right)$

非经典区域（ $E < V(x)$ ，动量为虚数）：

$\psi(x) \approx \frac{C}{\sqrt{|p(x)|}} \exp\left(\pm \frac{1}{\hbar}\int |p(x)| dx\right)$

graph TD
    A[定态薛定谔方程] -->|"假设 λ|dV/dx| ≪ |E-V|"| B[WKB近似]
    B -->|"E > V"| C["经典区域
振荡波函数"]
    B -->|"E < V"| D["非经典区域
指数衰减/增长"]
    
    style C fill:#e8f5e9
    style D fill:#ffebee

物理解释：

分母 $\sqrt{p(x)}$ ：粒子速度慢的地方（ $p$ 小），停留时间长，概率密度大
相位积分：累积的相位产生干涉效应

8.3.3 WKB近似的阶数

WKB近似可以系统地展开到更高阶：

零阶 $S_0$ ：给出经典轨迹和局部动量
一阶 $S_1$ ：给出振幅修正 $\propto 1/\sqrt{p}$
高阶：更精细的量子修正

在 $\hbar \to 0$ 的极限下，WKB恢复经典力学——这正是对应原理的体现。

8.4 转折点与连接公式

8.4.1 转折点的问题

在 $E = V(x)$ 处（转折点）， $p(x) = 0$ ，WKB波函数的分母为零——发散！

这是WKB近似的边界。在转折点附近，势能不能近似为常数，WKB失效。

graph TD
    A["经典区域 E>V"] -->|转折点| B["p(x) = 0"]
    B -->|WKB发散| C[失效区域]
    D["非经典区域 E|转折点| B
    
    style B fill:#ffebee
    style C fill:#fff3e0

8.4.2 Airy函数与线性近似

解决方法是：在转折点附近用线性近似 $V(x) \approx V(x_0) + V'(x_0)(x - x_0)$ 精确求解薛定谔方程。

这个方程的解是艾里函数（Airy Function） $\text{Ai}(\xi)$ ，其中 $\xi = (x - x_0)/\alpha$ ， $\alpha = (\hbar^2/2mV')^{1/3}$ 。

Airy函数有两个重要渐近行为：

$\xi \to +\infty$ ：指数衰减（对应非经典区域）
$\xi \to -\infty$ ：振荡（对应经典区域）

物理图像：Airy函数就像"量子转折处的标准件"——它平滑地连接了两侧的WKB波函数。

8.4.3 连接公式

核心结果：

对于左侧经典、右侧非经典的转折点（势垒）：

经典区域 ( $x < x_0$ ) 的振荡波函数：

$\psi \approx \frac{2C}{\sqrt{p}} \sin\left(\frac{1}{\hbar}\int_x^{x_0} p dx + \frac{\pi}{4}\right)$

连接到非经典区域 ( $x > x_0$ ) 的衰减波函数：

$\psi \approx \frac{C}{\sqrt{|p|}} \exp\left(-\frac{1}{\hbar}\int_{x_0}^x |p| dx\right)$

相位 $\pi/4$ 的来源：这正是Airy函数渐近展开中的Stokes相位。

graph TD
    subgraph 转折点连接
        A["经典区域
E > V"] -->|转折点| B["线性近似区
Airy函数"]
        B -->|连接公式| C["非经典区域
E < V"]
    end
    
    style A fill:#e8f5e9
    style B fill:#fff3e0
    style C fill:#ffebee

8.5 量子化条件的推导

8.5.1 势阱中的WKB量子化

考虑一个势阱，有两个转折点 $x_1$ 和 $x_2$ 。在经典区域 $x_1 < x < x_2$ 内，波函数必须是驻波——从两边反射回来的波相位匹配。

应用连接公式，要求波函数在势阱内自洽，得到Bohr-Sommerfeld量子化条件：

$\int_{x_1}^{x_2} p(x) dx = \left(n + \frac{1}{2}\right)\pi\hbar, \quad n = 0, 1, 2, ...$

物理意义：经典轨道的作用量必须是 $\hbar$ 的半整数倍。

为什么是 $n + \frac{1}{2}$ 而不是 $n$ ？

旧量子论给出的是 $\oint p dx = n h$ ，但WKB近似告诉我们，每个转折点贡献一个 $\pi/2$ 的相位偏移，两个转折点总共贡献 $\pi$ ，即半个波长的等效作用量。

graph TD
    A[势阱两个转折点] -->|连接公式相位匹配| B["每个转折点
贡献 π/2 相位"]
    B -->|"总相位 π"| C["n + ½ 量子化"]
    C -->|作用量量子化| D["∮p dx = (n+½)h"]
    
    style C fill:#e8f5e9
    style D fill:#fff3e0

8.5.2 数值例子：谐振子的WKB验证

问题：用WKB量子化条件推导一维谐振子能级。

解答：

对于谐振子 $V(x) = \frac{1}{2}m\omega^2 x^2$ ，转折点在 $x = \pm a$ ，其中 $a = \sqrt{2E/m\omega^2}$ 。

计算积分：

$\int_{-a}^{a} p(x) dx = \int_{-a}^{a} \sqrt{2m\left(E - \frac{1}{2}m\omega^2 x^2\right)} dx$

令 $x = a\sin\theta$ ，则 $dx = a\cos\theta d\theta$ ，积分限变为 $-\pi/2$ 到 $\pi/2$ ：

$p(x) = \sqrt{2mE - m^2\omega^2 x^2} = m\omega a \cos\theta$

$\int_{-a}^{a} p(x) dx = m\omega a^2 \int_{-\pi/2}^{\pi/2} \cos^2\theta d\theta = m\omega a^2 \cdot \frac{\pi}{2} = \frac{\pi E}{\omega}$

应用量子化条件：

$\frac{\pi E}{\omega} = \left(n + \frac{1}{2}\right)\pi\hbar$

$E_n = \left(n + \frac{1}{2}\right)\hbar\omega$

WKB给出了精确结果！ 对于谐振子，WKB近似恰好是精确的。

8.6 隧穿概率的WKB公式

8.6.1 势垒穿透

对于宽度为 $a$ 、高度为 $V_0$ 的方势垒， $E < V_0$ ：

经典计算透射系数（WKB近似下的主导项）：

$T \approx \exp\left(-2\int_{x_1}^{x_2} \frac{|p(x)|}{\hbar} dx\right) = \exp\left(-\frac{2}{\hbar}\int_{x_1}^{x_2} \sqrt{2m(V(x) - E)} dx\right)$

对于矩形势垒：

$T \approx \exp\left(-2\kappa a\right), \quad \kappa = \frac{\sqrt{2m(V_0 - E)}}{\hbar}$

8.6.2 数值例子：方势垒隧穿

问题：电子（质量 $m = 9.11 \times 10^{-31}$ kg）遇到高度 $V_0 = 5$ eV、宽度 $a = 1$ nm = $10^{-9}$ m 的势垒，入射能量 $E = 2$ eV。估算隧穿概率。

解答：

步骤1：计算 $\kappa$

$V_0 - E = 3 \text{ eV} = 3 \times 1.6 \times 10^{-19} \text{ J} = 4.8 \times 10^{-19} \text{ J}$

$\kappa = \frac{\sqrt{2m(V_0-E)}}{\hbar} = \frac{\sqrt{2 \times 9.11 \times 10^{-31} \times 4.8 \times 10^{-19}}}{1.055 \times 10^{-34}}$

$= \frac{\sqrt{8.75 \times 10^{-49}}}{1.055 \times 10^{-34}} = \frac{9.35 \times 10^{-25}}{1.055 \times 10^{-34}} \approx 8.86 \times 10^{9} \text{ m}^{-1}$

步骤2：计算指数

$2\kappa a = 2 \times 8.86 \times 10^{9} \times 10^{-9} = 17.72$

步骤3：隧穿概率

$T \approx e^{-17.72} \approx 2.0 \times 10^{-8}$

结果：隧穿概率约为 $2 \times 10^{-8}$ ，即约两千万分之一。

物理意义：即使是纳米尺度的势垒，电子的隧穿概率也非常小——这正是经典物理直觉所预期的。但在原子尺度（ $a \sim 1$ Å），隧穿概率会显著增大。

8.7 α衰变：量子隧穿的经典应用

8.7.1 物理图像

放射性α衰变是量子隧穿的经典例子。α粒子（氦核）在原子核内被强核力束缚，但核外有库仑势垒。

经典物理：α粒子能量 $E$ 远低于势垒峰值，不可能逃出。
量子力学：α粒子以一定概率隧穿库仑势垒。

8.7.2 伽莫夫的WKB计算

伽莫夫（Gamow）的模型：

α粒子在核内自由运动（核力作用区），以频率 $f$ 撞击势垒壁。每次撞击有概率 $T$ 隧穿而出。

衰变常数：

$\lambda = f \cdot T$

其中 $f \approx v/2R \approx 10^{21}$ s $^{-1}$ （α粒子在核内的运动频率）。

8.7.3 数值例子：U-238的α衰变

问题：估算铀-238的α衰变半衰期。

已知数据：

α粒子能量： $E = 4.27$ MeV $= 4.27 \times 10^6 \times 1.6 \times 10^{-19}$ J $= 6.83 \times 10^{-13}$ J
铀核电荷： $Z = 92$
α粒子质量： $m = 6.64 \times 10^{-27}$ kg
核半径： $R \approx 1.2 \times A^{1/3}$ fm $= 1.2 \times 238^{1/3} \approx 7.44$ fm $= 7.44 \times 10^{-15}$ m

解答：

步骤1：确定势垒和转折点

库仑势： $V(r) = \frac{2Ze^2}{4\pi\varepsilon_0 r}$

经典转折点（ $V(r_2) = E$ ）：

$r_2 = \frac{2Ze^2}{4\pi\varepsilon_0 E} = \frac{2 \times 92 \times (1.6 \times 10^{-19})^2}{4\pi \times 8.85 \times 10^{-12} \times 6.83 \times 10^{-13}}$

$= \frac{4.71 \times 10^{-36}}{7.58 \times 10^{-23}} \approx 6.21 \times 10^{-14} \text{ m} = 62.1 \text{ fm}$

内转折点 $r_1 = R \approx 7.44$ fm。

步骤2：计算势垒积分

$\gamma = \int_{r_1}^{r_2} \sqrt{2m(V(r)-E)} dr$

令 $V(r) = E \cdot \frac{r_2}{r}$ ，则：

$\gamma = \sqrt{2mE} \int_{r_1}^{r_2} \sqrt{\frac{r_2}{r} - 1} dr$

这个积分可以用变量替换 $r = r_2 \sin^2\theta$ 求解：

$\gamma = \sqrt{2mE} \cdot r_2 \left[\arccos\sqrt{\frac{r_1}{r_2}} - \sqrt{\frac{r_1}{r_2}\left(1-\frac{r_1}{r_2}\right)}\right]$

代入数值：

$r_1/r_2 = 7.44/62.1 \approx 0.12$
$\arccos(\sqrt{0.12}) = \arccos(0.346) \approx 1.21$ rad
$\sqrt{0.12 \times 0.88} \approx 0.325$

$\gamma = \sqrt{2 \times 6.64 \times 10^{-27} \times 6.83 \times 10^{-13}} \times 6.21 \times 10^{-14} \times (1.21 - 0.325)$

= 9.52 \times 10^{-20} \times 6.21 \times 10^{-14} \times 0.885 = 5.23 \times 10^{-33} \text{（SI单位）}

转换为无量纲形式（除以 $\hbar$ ）：

$\frac{\gamma}{\hbar} = \frac{5.23 \times 10^{-33}}{1.055 \times 10^{-34}} \approx 49.6$

步骤3：计算隧穿概率和半衰期

$T \approx e^{-2\gamma/\hbar} = e^{-99.2} \approx 1.3 \times 10^{-43}$

衰变常数：

$\lambda = f \cdot T \approx 10^{21} \times 1.3 \times 10^{-43} = 1.3 \times 10^{-22} \text{ s}^{-1}$

半衰期：

T_{1/2} = \frac{\ln 2}{\lambda} = \frac{0.693}{1.3 \times 10^{-22}} \approx 5.3 \times 10^{21} \text{ s} \approx 1.7 \times 10^{14} \text{ 年}

实验值： $T_{1/2} \approx 4.5 \times 10^{9}$ 年

说明：这个粗略估算（忽略了核势的细节、角动量等因素）给出的数量级与实验值在同一区间（都是 $10^9$ 年量级），但具体数值有偏差。更精确的计算需要考虑：

核表面势的精确形状
α粒子的角动量（ $l \neq 0$ 时势垒更高）
核内预形成因子

8.7.4 Geiger-Nuttall定律

Geiger-Nuttall定律：

$\ln T_{1/2} = a + \frac{b}{\sqrt{E}}$

半衰期的对数与粒子能量的平方根成反比。WKB近似完美地解释了这一经验规律——因为隧穿概率的对数正比于势垒积分 $\gamma$ ，而 $\gamma$ 与 $1/\sqrt{E}$ 相关。

graph TD
    A["α粒子在核内"] -->|库仑势垒| B[经典禁戒区]
    B -->|量子隧穿| C[衰减波函数]
    C -->|"透射概率 T"| D["α衰变"]
    D -->|WKB解释| E["Geiger-Nuttall定律
ln T₁/₂ ∝ 1/√E"]
    
    style C fill:#ffebee
    style D fill:#e8f5e9
    style E fill:#fff3e0

8.8 扫描隧道显微镜（STM）

扫描隧道显微镜是WKB隧穿公式最直接的技术应用。

原理：金属探针靠近导电样品表面（间距约几个Å）。探针和样品之间是真空——经典物理说电子不可能穿过。但量子隧穿允许电子穿过这个"真空势垒"。

隧穿电流：

$I \propto e^{-2\gamma d}$

其中 $d$ 是探针-样品间距， $\gamma = \sqrt{2m\phi}/\hbar$ ， $\phi$ 是功函数。

关键洞察：电流对间距极度敏感—— $d$ 增加1Å，电流衰减约一个数量级。这种指数敏感性让STM获得了原子级分辨率。

1986年，宾尼格和罗雷尔因发明STM获得诺贝尔物理学奖。

graph TD
    A[金属探针] -->|"真空势垒
几个Å"| B[样品表面]
    A -->|电子隧穿| C[隧穿电流]
    C -->|"I ∝ e^(-2γd)"| D[原子级分辨率]
    D -->|扫描成像| E["STM图像
看见单原子"]
    
    style D fill:#e8f5e9
    style E fill:#fff3e0

8.9 共振隧穿与双势垒结构

共振隧穿是WKB近似的一个精彩应用。考虑两个势垒之间夹着一个势阱（双势垒结构）。

经典物理：粒子必须克服两个势垒才能通过。
量子力学：如果粒子的能量恰好等于势阱内的某个"准束缚态"能级，隧穿概率会急剧增大——甚至接近1！

WKB分析：双势垒的透射系数在相位共振条件下可以接近1。

技术应用：共振隧穿二极管（RTD）是超高速电子器件的核心。

graph TD
    A[入射粒子] -->|势垒1| B["势阱
准束缚态"]
    B -->|势垒2| C[透射]
    B -->|"能量匹配
共振条件"| D["透射率 ≈ 1"]
    D -->|应用| E["共振隧穿二极管
RTD"]
    
    style D fill:#e8f5e9
    style E fill:#fff3e0

8.10 本章总结

带走的三句话：

WKB近似适用于势能变化缓慢的区域：把波函数看作局部平面波。
转折点是经典与量子的边界：WKB失效，但连接公式可以桥接两侧。
WKB量子化条件解释了为什么经典轨道的作用量必须是ℏ的半整数倍。

"WKB是量子力学中最接近直觉的近似。它告诉你：在大多数情况下，粒子确实像经典一样运动——但波的相位累积会带来量子奇迹。" —— Griffiths

8.11 练习与思考

谐振子WKB验证：用WKB量子化条件推导谐振子能级 $E_n = (n + 1/2)\hbar\omega$ 。计算积分 $\int_{-a}^{a} \sqrt{2m(E - \frac{1}{2}m\omega^2 x^2)} dx$ ，其中 $a = \sqrt{2E/m\omega^2}$ 。
方势垒隧穿：电子能量 $E = 2$ eV，遇到高度 $V_0 = 5$ eV、宽度 $a = 1$ nm 的矩形势垒。计算 $\kappa = \sqrt{2m(V_0-E)}/\hbar$ 和隧穿概率 $T \approx e^{-2\kappa a}$ 。
α衰变估算：假设α粒子（质量 $m \approx 6.6 \times 10^{-27}$ kg，能量 $E \approx 5$ MeV）面对一个库仑势垒，峰值约 $V_{max} \approx 30$ MeV，有效宽度 $a \approx 10$ fm ( $10^{-14}$ m)。用WKB公式估算隧穿概率 $T$ 的数量级。
WKB失效条件：WKB近似在什么情况下失效？列举三种情况，并说明为什么。提示：考虑转折点附近的性质、势能突变、以及 $E \approx V$ 的区域。