薛定谔方程的WKB近似与半经典转折点连接条件

求解一维定态薛定谔方程时，当粒子的动能远大于势能的变化（即量子效应很弱），或者粒子的德布罗意波长在势场中变化缓慢时，可以使用半经典近似方法。其中最经典、最实用的就是WKB近似。它的核心思想是：当势场变化足够缓慢时，粒子在局部区域内的行为近似于平面波，波函数的形式可以写为振幅乘以一个快速振荡的相位因子。

1. 建立标准WKB近似解

首先，写下一维定态薛定谔方程。

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

其中，$m$ 是粒子质量，$V(x)$ 是势能，$E$ 是总能量。

定义 波数函数 $k(x)$ 和衰减常数 $\kappa(x)$：

• 在经典允许区（$E > V(x)$）：$k(x) = \frac{1}{\hbar} \sqrt{2m(E - V(x))}$
• 在经典禁止区（$E < V(x)$）：$\kappa(x) = \frac{1}{\hbar} \sqrt{2m(V(x) - E)}$

应用 WKB近似条件。当势场满足以下条件时，近似成立：
$$ \left| \frac{d}{dx} \left( \frac{1}{k(x)} \right) \right| \ll 1 \quad \text{或等价地} \quad \left| \frac{m}{\hbar^2 k(x)^3} \frac{dV}{dx} \right| \ll 1 $$
这个条件的物理意义是：在一个德布罗意波长 $\lambda = 2\pi/k(x)$ 的距离内，势能 $V(x)$ 的变化远小于局部动能 $(E-V)$。

得到 标准WKB近似波函数。

• 在经典允许区 ($E > V$)：
  $$ \psi(x) \approx \frac{A}{\sqrt{k(x)}} \exp\left( i \int^x k(x') dx' \right) + \frac{B}{\sqrt{k(x)}} \exp\left( -i \int^x k(x') dx' \right) $$
• 在经典禁止区 ($E < V$)：
  $$ \psi(x) \approx \frac{C}{\sqrt{\kappa(x)}} \exp\left( \int^x \kappa(x') dx' \right) + \frac{D}{\sqrt{\kappa(x)}} \exp\left( -\int^x \kappa(x') dx' \right) $$
  其中，$A, B, C, D$ 是待定常数。因子 $1/\sqrt{k}$ 或 $1/\sqrt{\kappa}$ 是为了满足流密度守恒而引入的。

2. 识别转折点与连接问题的根源

WKB近似在势能满足 $E = V(x)$ 的点（称为转折点）附近失效。这是因为在这些点，$k(x)$ 或 $\kappa(x)$ 趋于零，WKB近似条件中的分母变为零，近似条件被破坏。

观察 转折点两侧的波函数行为差异巨大。假设有一个转折点 $x=a$，且 $x<a$ 是经典禁止区 ($E<V$)，$x>a$ 是经典允许区 ($E>V$)。

• 在经典禁止区 ($x<a$)，WKB解是指数衰减和增长的。物理上可接受的解通常只保留衰减项 $\exp(-\int \kappa dx)$。
• 在经典允许区 ($x>a$)，WKB解是振荡的行波或驻波。

理解 核心困难：我们无法直接用WKB波函数去匹配 $x=a$ 点两侧的解，因为WKB近似在该点完全失效。我们必须找到一种“连接公式”，能够跨越转折点，将两侧的WKB解（特别是它们的系数 $A, B, C, D$）联系起来。

3. 推导并应用连接公式

连接公式的严格推导需要在转折点 $x=a$ 附近，将势能 $V(x)$ 线性化（泰勒展开至一阶），即 $V(x) \approx V(a) + V'(a)(x-a) = E + F(x-a)$，其中 $F=V'(a)$ 是常数力。然后求解这个线性势下的精确解（Airy函数），并令其在远离转折点处与两侧的WKB渐进行为匹配，从而得到连接系数。

记忆并应用 最常用的两组连接公式。设转折点 $x=a$ 左侧 ($x<a$) 是经典禁止区，右侧 ($x>a$) 是经典允许区。

公式组A：指数衰减解 → 驻波解
如果粒子从左侧经典禁止区（取衰减解）入射，则连接公式为：
$$ \frac{2}{\sqrt{k}} \cos\left( \int_a^x k(x') dx' - \frac{\pi}{4} \right) \longleftrightarrow \frac{1}{\sqrt{\kappa}} \exp\left( -\int_x^a \kappa(x') dx' \right) $$

• 左侧：代表经典允许区的一个余弦驻波，其相位参考点在转折点 $x=a$，并有一个关键的 $-\pi/4$ 相位修正。
• 右侧：代表经典禁止区的一个衰减指数波。

公式组B：指数增长解 → 行波解
如果粒子从右侧经典允许区入射，则需要另一组公式。更通用的形式是使用复指数表示行波：
$$ \frac{1}{\sqrt{k}} \exp\left( i \int_a^x k(x') dx' - i\frac{\pi}{4} \right) \longleftrightarrow \frac{-i}{\sqrt{\kappa}} \exp\left( \int_x^a \kappa(x') dx' \right) $$
其复共轭关系也成立：
$$ \frac{1}{\sqrt{k}} \exp\left( -i \int_a^x k(x') dx' + i\frac{\pi}{4} \right) \longleftrightarrow \frac{i}{\sqrt{\kappa}} \exp\left( -\int_x^a \kappa(x') dx' \right) $$

• 左侧：代表向右传播的波（$e^{iS}$ 形式），同样带有 $-\pi/4$ 的相位修正。
• 右侧：代表指数增长的解。在物理问题中，通常需要与另一侧的衰减解组合，因此增长解是数学上必要的。

注意 积分方向：公式中积分 $\int_a^x$ 表示从转折点 $a$ 积分到点 $x$。

4. 解决实际问题的通用步骤

假设需要求解一个束缚态问题，势阱 $V(x)$ 有两个转折点 $x_1$ 和 $x_2$（$x_1 < x_2$），且 $V(x_1) = V(x_2) = E$。在 $x<x_1$ 和 $x>x_2$ 是经典禁止区，在 $x_1 < x < x_2$ 是经典允许区。

步骤1：写出各区域的WKB波函数形式。

• 区域I ($x < x_1$)：$\psi_I \approx \frac{D}{\sqrt{\kappa}} \exp\left( -\int_x^{x_1} \kappa dx' \right)$ （只保留物理上有限的衰减解）。
• 区域II ($x_1 < x < x_2$)：$\psi_{II} \approx \frac{A}{\sqrt{k}} \exp\left( i\int_{x_1}^x k dx' \right) + \frac{B}{\sqrt{k}} \exp\left( -i\int_{x_1}^x k dx' \right)$。
• 区域III ($x > x_2$)：$\psi_{III} \approx \frac{C}{\sqrt{\kappa}} \exp\left( -\int_{x_2}^x \kappa dx' \right)$ （只保留衰减解）。

步骤2：在左侧转折点 $x_1$ 应用连接公式。
观察区域I的解，它对应连接公式组A的右侧。应用 公式组A，将其转换为区域II左侧 ($x$ 略大于 $x_1$) 的波函数：
$$ \psi_{II}(x \gtrsim x_1) \leftrightarrow \frac{2D}{\sqrt{k}} \cos\left( \int_{x_1}^x k dx' - \frac{\pi}{4} \right) $$
展开 余弦项，得到：
$$ \psi_{II}(x \gtrsim x_1) = \frac{2D}{\sqrt{k}} \left[ \frac{1}{2} \exp\left( i\int_{x_1}^x k dx' - i\frac{\pi}{4} \right) + \frac{1}{2} \exp\left( -i\int_{x_1}^x k dx' + i\frac{\pi}{4} \right) \right] $$
比较 这个表达式与区域II的通解 $\psi_{II}$，确定 系数 $A$ 和 $B$ 与 $D$ 的关系：
$$ A = D e^{-i\pi/4}, \quad B = D e^{i\pi/4} $$
这意味着在 $x_1$ 点，向左传播的波 ($B$) 和向右传播的波 ($A$) 振幅相等，形成了纯驻波。

步骤3：在右侧转折点 $x_2$ 应用连接公式。
我们需要将区域II的通解 $\psi_{II}$ 在 $x_2$ 附近的形式，与区域III的衰减解 $\psi_{III}$ 连接起来。
首先，将 $\psi_{II}$ 改写为以 $x_2$ 为参考点的形式。利用 $\int_{x_1}^{x_2} k dx' = \Phi$ （定义为从 $x_1$ 到 $x_2$ 的“相位积分”），则：
$$ \psi_{II}(x) = \frac{1}{\sqrt{k}} \left[ A e^{i\Phi} \exp\left( i\int_{x_2}^x k dx' \right) + B e^{-i\Phi} \exp\left( -i\int_{x_2}^x k dx' \right) \right] $$
现在，关注 在 $x$ 略小于 $x_2$ 处，向左传播的波（第二项 $e^{-i\int_{x_2}^x k dx'}$）和向右传播的波（第一项 $e^{i\int_{x_2}^x k dx'}$）。应用 连接公式组B及其共轭，将区域III的衰减解与向左传播的波关联起来。关键在于，右侧的衰减解必须由一个纯粹的向左传播的波（相位为 $e^{-i(\dots - \pi/4)}$）连接而来。因此，要求 向右传播的波在 $x_2$ 处为零，即：
$$ A e^{i\Phi} = 0 $$
但这会导致平凡解。实际上，正确的匹配要求向右传播的波与向左传播的波满足特定的复数关系。进行 更细致的匹配后，可得到束缚态量子化条件：
$$ \int_{x_1}^{x_2} k(x) dx = \left( n + \frac{1}{2} \right) \pi, \quad n = 0, 1, 2, \dots $$
其中，$\int_{x_1}^{x_2} k(x) dx$ 就是上面定义的相位积分 $\Phi$。

步骤4：解释量子化条件的物理意义。
得到 最终的量子化条件。这个条件决定了能量 $E$ 只能取一系列离散的值 $E_n$。

• 条件的左边是粒子在经典允许区内运动一个完整周期（从 $x_1$ 到 $x_2$ 再返回）所积累的相位。
• 条件的右边表明，这个相位必须是 $\pi$ 的半整数倍。额外的 $\pi/2$ 来自于两个转折点各自贡献的 $-\pi/4$ 相位修正（$-\pi/4 - (-\pi/4) = -\pi/2$，取绝对值贡献 $\pi/2$）。

5. 总结关键操作要点

1. 判断 是否满足WKB近似条件：势场变化相对粒子波长是否足够缓慢。
2. 识别 问题的所有转折点 $E = V(x)$。
3. 划分 区域：在转折点之间是经典允许区（用振荡解），之外是经典禁止区（用指数衰减/增长解）。
4. 书写 各区域的WKB波函数通解形式。
5. 应用 连接公式：在每一个转折点处，使用正确的连接公式组（指数衰减解 ↔ 驻波解，指数增长解 ↔ 行波解），注意 $-\pi/4$ 的相位修正和积分方向。
6. 匹配 波函数在转折点两侧的系数，建立方程。
7. 求解 方程组。对于束缚态问题，最终通常会导出一个形如 $\int k dx = (n+1/2)\pi$ 的量子化条件，用以求解能量本征值 $E_n$。