第2章 定态薛定谔方程：量子化的自然结果

"如果你把粒子关在一个盒子里，它不会像经典小球那样在壁上反弹——它会形成驻波，只允许特定的能量存在。"

这就是量子化的起源。不是人为假设，而是薛定谔方程的自然结果。

---

第2章 定态薛定谔方程：盒子里的粒子会唱歌

0. 前置知识：常微分方程基础

定态薛定谔方程本质上是一个二阶常微分方程（ODE）。在深入求解之前，我们需要复习几个关键技巧。

0.1 分离变量法

如果方程可以写成：

$\frac{d^2\psi}{dx^2} = f(x)\psi$

我们尝试找到形如 $\psi(x) = X(x)$ 的解，其中 $X$ 满足特定边界条件。

核心思想：将偏微分方程拆成多个常微分方程。第2章的核心技巧就是分离变量——将时间依赖的薛定谔方程拆成空间部分 $\psi(x)$ 和时间部分 $\phi(t)$。

0.2 边界条件与本征值问题

二阶ODE需要两个边界条件才能唯一确定解。在量子力学中，边界条件不是任意的，而是由物理要求决定的：

1. 波函数连续性：$\psi$ 必须在空间连续（势能有限处）
2. 导数连续性：$d\psi/dx$ 必须在空间连续（势能有限处）
3. 归一化条件：$\int |\psi|^2 dx = 1$（束缚态）
4. 渐近行为：$x \to \pm\infty$ 时 $\psi \to 0$（束缚态）

本征值问题：当边界条件"太严格"时，方程只对特定的参数值（本征值 $E_n$）才有非零解。这就是量子化的数学根源——不是物理假设，而是边界条件的必然结果。

类比：一根琴弦，两端固定（边界条件）。只有特定频率的振动模式才能存在——这些频率就是"本征值"，对应的振动模式就是"本征态"。

0.3 数值例题：简单ODE的边界条件

问题：方程 $\frac{d^2\psi}{dx^2} = -k^2\psi$，边界条件 $\psi(0) = \psi(L) = 0$。求允许的非零解。

解答：
通解为 $\psi(x) = A\sin(kx) + B\cos(kx)$。

由 $\psi(0) = 0$：$B = 0$。

由 $\psi(L) = 0$：$A\sin(kL) = 0$。

若 $A \neq 0$（非零解），则 $\sin(kL) = 0$，即 $kL = n\pi$，$n = 1, 2, 3, ...$

所以 $k_n = \frac{n\pi}{L}$，解为 $\psi_n(x) = A\sin(\frac{n\pi x}{L})$。

这就是无限深势阱的雏形——边界条件直接导致了 $k$（从而能量）的量子化。

---

2.1 故事：量子琴师的调音难题

2099年，"量子音乐厅"是全球最热门的演出场地。这里的乐器不是传统的弦乐，而是量子阱——将单个电子囚禁在纳米尺度的势阱中，通过激光激发，电子在不同能级间跃迁时发出纯净的单色光。

但今晚，首席量子琴师陈默遇到了麻烦。

"这台量子钢琴跑调了。"她盯着频谱仪，"理论上，把势阱宽度从2nm改成3nm，发射波长应该从520nm移到……不对，实际测出来完全不是线性关系！"

她的助手——一台古董AI——慢悠悠地回答："因为电子在势阱里不是’自由移动’，而是形成驻波。只有半波长的整数倍等于阱宽的模式才被允许。改变阱宽，允许的模式会按平方关系变化，不是线性的。"

"什么意思？"

"意思就是：量子世界的音阶不是连续的，而是离散的。你不能任意调音，只能跳到特定的’音符’上。"

陈默愣了几秒，然后笑了："所以量子钢琴不是钢琴，而是……竖琴？"

"可以这么说。只不过琴弦是势阱壁，音符是能量本征值。"

核心洞察：束缚导致离散。粒子被限制在空间中时，能量只能取特定值——这就是量子化的本质。

---

2.2 定态：薛定谔方程的可分离解

第1章的时间依赖薛定谔方程是：

$i\hbar \frac{\partial \Psi}{\partial t} = \hat{H}\Psi$

如果势能 $V(x)$ 不随时间变化，我们可以用分离变量法：假设 $\Psi(x,t) = \psi(x)\phi(t)$。

2.2.1 分离变量的结果

代入后得到两个方程：

时间部分：

$\phi(t) = e^{-iEt/\hbar}$

空间部分（定态薛定谔方程）：

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

或者写成哈密顿算符的形式：

$\hat{H}\psi = E\psi$

这就是量子力学中最重要的方程之一：能量本征值方程。它的解 $\psi_n(x)$ 称为能量本征态，对应的 $E_n$ 是能量本征值。

历史背景：从波动方程到本征值问题

分离变量法不是薛定谔发明的——它是19世纪数学家求解波动方程（如琴弦振动、热传导）的标准工具。但薛定谔的天才在于：他将这一经典数学工具用于原子中的"物质波"，从而将量子化问题转化为本征值问题。

在经典物理中，本征值问题描述的是琴弦的固有频率、鼓膜的振动模式。在量子力学中，本征值描述的是原子允许存在的能量。数学形式相同，物理内涵完全不同。

graph TD
    subgraph 分离变量法的统一性
        A[经典琴弦振动] --> B[波动方程]
        C[量子物质波] --> B
        B --> D[分离变量]
        D --> E[空间方程]
        D --> F[时间方程]
        E --> G["经典: 固有频率"]
        E --> H["量子: 允许能量"]
    end
    
    style G fill:#e8f5e9
    style H fill:#fff3e0

2.2.2 为什么叫"定态"？

定态的完整波函数是：

$\Psi_n(x,t) = \psi_n(x)e^{-iE_nt/\hbar}$

注意概率密度：

$|\Psi_n(x,t)|^2 = |\psi_n(x)|^2$

概率密度不随时间变化！ 粒子的"云"静止在空间中，这就是"定态"的含义。

graph TD
    A[时间依赖薛定谔方程] -->|"V(x)不含时"| B[分离变量]
    B --> C["空间部分: Hψ = Eψ"]
    B --> D["时间部分: e^(-iEt/ℏ)"]
    C --> E["能量本征态 ψₙ"]
    D --> E
    E --> F["定态: |Ψ|² = |ψ|² 不随时间变化"]
    
    style F fill:#e8f5e9

物理直觉：定态不是"静止不动"

定态并不意味着粒子"凝固"在空间中。恰恰相反，定态是粒子所有可能运动方式的稳定模式——就像琴弦的驻波不是弦不动，而是弦以固定模式振动。

在定态中，概率分布不变，但波函数本身随时间演化（乘以相位因子 $e^{-iE_nt/\hbar}$）。这个相位因子不影响任何物理测量结果（因为概率只依赖模平方），但它在叠加态中至关重要——不同能量的相位以不同速率旋转，导致干涉图案随时间变化。

---

2.3 无限深方势阱：最简单的"量子牢房"

2.3.1 问题设置

势能函数：

$V(x) = \begin{cases} 0, & 0 \leq x \leq a \\ \infty, & \text{otherwise} \end{cases}$

电子被严格限制在 $[0,a]$ 区间内，外面势能无穷大，粒子不可能逃逸。

2.3.2 求解

在阱内 ($0 < x < a$)，定态薛定谔方程变为：

$-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} = E\psi$

这是标准的波动方程。解为正弦和余弦函数的组合。加上边界条件 $\psi(0) = \psi(a) = 0$（波函数必须在壁处为零），我们得到：

$\psi_n(x) = \sqrt{\frac{2}{a}}\sin\left(\frac{n\pi x}{a}\right), \quad n = 1, 2, 3, ...$

能量本征值：

$E_n = \frac{n^2\pi^2\hbar^2}{2ma^2}$

graph TD
    subgraph 无限深方势阱的波函数
        A["n=1"] --> B["sin(πx/a)"]
        C["n=2"] --> D["sin(2πx/a)"]
        E["n=3"] --> F["sin(3πx/a)"]
    end
    
    style B fill:#e3f2fd
    style D fill:#e8f5e9
    style F fill:#fff3e0

例题详解：归一化常数的确定

问题：证明无限深势阱波函数的归一化常数是 $\sqrt{2/a}$。

分步推导：

第一步：写出归一化条件

$\int_0^a |\psi_n|^2 dx = \int_0^a A^2 \sin^2\left(\frac{n\pi x}{a}\right) dx = 1$

第二步：利用三角恒等式

$\sin^2\theta = \frac{1 - \cos(2\theta)}{2}$

$\int_0^a A^2 \cdot \frac{1 - \cos(2n\pi x/a)}{2} dx = 1$

第三步：拆分积分

$\frac{A^2}{2}\left[\int_0^a dx - \int_0^a \cos\left(\frac{2n\pi x}{a}\right) dx\right] = 1$

第四步：第一项是 $a$；第二项是完整周期积分，结果为0

$\frac{A^2}{2} \cdot a = 1$

第五步：解出 $A$

$A = \sqrt{\frac{2}{a}}$

物理意义：势阱越宽，波函数整体幅度越小——因为粒子"摊"在更大的空间内，每个位置的概率密度必须降低以保持总概率为1。

2.3.3 关键特征

特征	经典物理	量子力学
能量	可以任意连续取值	离散：$E_n \propto n^2$
最低能量	可以为零（静止在底部）	$E_1 > 0$，存在零点能
阱内分布	均匀分布	有节点，概率分布不均匀
节点数	无	$n-1$ 个节点

零点能的存在是海森堡不确定性原理的直接后果：如果粒子能量为零，动量也为零，那 $\sigma_p = 0$，意味着 $\sigma_x \to \infty$——但粒子被限制在阱宽 $a$ 内，$\sigma_x$ 最大也只有 $a$。所以能量不能为零。

graph TD
    subgraph 零点能的物理起源
        A[粒子被限制在阱内] --> B["位置不确定度有限
σₓ ~ a"]
        B --> C[不确定性原理]
        C --> D["动量不确定度非零
σₚ ≥ ℏ/(2σₓ)"]
        D --> E["动能非零
E = p²/2m > 0"]
        E --> F[零点能]
    end
    
    style F fill:#ffebee

2.3.4 故事续集：音阶的秘密

陈默回到实验室，在白板上推导：

"对于宽度 $a = 2$ nm 的势阱，基态能量是：

$E_1 = \frac{\pi^2\hbar^2}{2ma^2} \approx 0.094 \text{ eV}$

第一激发态 $E_2 = 4E_1$，第二激发态 $E_3 = 9E_1$……"

她停住了："所以能级差不是等间距的！$E_2 - E_1 = 3E_1$，$E_3 - E_2 = 5E_1$……"

"对。"AI说，"这就是为什么量子钢琴的音阶不是等距的。而且当你把阱宽从2nm改成3nm，基态能量会按 $1/a^2$ 缩放，变成原来的 $4/9$。"

"所以波长变长了，颜色从绿光移向红光……"

"正是如此。"

例题详解：能级的数值估算

问题：一个电子被限制在宽度 $a = 1$ nm 的无限深势阱中。计算前三个能级的能量（以eV为单位），并求从 $n=2$ 跃迁到 $n=1$ 时发射光子的波长。

分步推导：

第一步：代入常数

\hbar = 1.055 \times 10^{-34} \text{ J·s}

$m = 9.11 \times 10^{-31} \text{ kg}$

$a = 10^{-9} \text{ m}$

第二步：计算 $E_1$

$E_1 = \frac{\pi^2 (1.055 \times 10^{-34})^2}{2(9.11 \times 10^{-31})(10^{-9})^2}$

$= \frac{\pi^2 \times 1.113 \times 10^{-68}}{1.822 \times 10^{-48}}$

$\approx 6.02 \times 10^{-20} \text{ J}$

转换为eV（$1 \text{ eV} = 1.602 \times 10^{-19} \text{ J}$）：

$E_1 \approx 0.376 \text{ eV}$

第三步：计算前三个能级

$E_1 \approx 0.376 \text{ eV}$

$E_2 = 4E_1 \approx 1.504 \text{ eV}$

$E_3 = 9E_1 \approx 3.384 \text{ eV}$

第四步：计算跃迁能量和波长

$\Delta E = E_2 - E_1 = 3E_1 \approx 1.128 \text{ eV}$

光子波长：

\lambda = \frac{hc}{\Delta E} = \frac{(4.136 \times 10^{-15} \text{ eV·s})(2.998 \times 10^8 \text{ m/s})}{1.128 \text{ eV}}

$\approx 1.10 \times 10^{-6} \text{ m} = 1100 \text{ nm}$

结论：这个跃迁发射的是红外光（刚好在近红外波段）。如果势阱更窄（比如 $a = 0.5$ nm），能量会增大4倍，波长缩短到约275 nm（紫外）。

---

2.4 谐振子：量子世界的弹簧

2.4.1 问题设置

势能：

$V(x) = \frac{1}{2}m\omega^2 x^2$

这是量子力学中最重要的模型之一，因为任何系统在平衡点附近的微小振动都可以近似为谐振子。

graph TD
    subgraph 谐振子的普适性
        A[分子振动] --> B[谐振子近似]
        C[晶格振动] --> B
        D[电磁场模式] --> B
        E[量子场论激发] --> B
        B --> F[统一描述]
    end
    
    style B fill:#e8f5e9
    style F fill:#fff3e0

2.4.2 求解与结果

通过巧妙的变量替换和级数解法（Griffiths在书中详细展示了推导），得到：

能量本征值：

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

关键区别：谐振子的能级是等间距的！每级相差 $\hbar\omega$。这与势阱的 $n^2$ 依赖完全不同。

graph LR
    subgraph 能级对比
        A[无限深方势阱] --> B["Eₙ ∝ n²"]
        C[谐振子] --> D["Eₙ ∝ n+½"]
    end
    
    style B fill:#ffebee
    style D fill:#e8f5e9

波函数：用厄米多项式 $H_n(y)$ 表示，其中 $y = x\sqrt{m\omega/\hbar}$：

$\psi_n(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\frac{1}{\sqrt{2^n n!}} H_n(y) e^{-y^2/2}$

基态 ($n=0$) 是一个高斯波包——粒子最可能出现在平衡位置，但也有一定概率出现在其他地方。

物理直觉：为什么能级等间距？

势能 $V(x) = \frac{1}{2}m\omega^2 x^2$ 是二次型。无论粒子能量多高，它总是在"同样形状"的势阱里运动。能量的增加只是让粒子"爬得更高"，但每爬升相同的高度 $\hbar\omega$，就会遇到一个新的允许模式。

相比之下，无限深势阱的"壁"是垂直的——能量越高，粒子运动越快，波长越短，允许的节点数越多。能级不是等间距的，因为每个新增模式需要"挤进"更小的空间。

graph TD
    subgraph "等间距 vs 不等间距"
        A[谐振子] --> B["势能: 二次型
V ∝ x²"]
        B --> C["等间距能级
ΔE = ℏω"]
        
        D[无限深势阱] --> E["势能: 垂直壁
V = 0或∞"]
        E --> F["不等间距能级
ΔE ∝ n"]
    end
    
    style C fill:#e8f5e9
    style F fill:#fff3e0

2.4.3 物理意义

谐振子的等间距能级解释了：

• 黑体辐射谱（光子能量量子化）
• 晶格振动（声子）
• 量子场论的基本激发模式

历史背景：普朗克的黑体辐射

1900年，普朗克为了解释黑体辐射曲线，被迫假设能量是量子化的：$E = n\hbar\omega$。当时他并不清楚这个假设的深层原因——只是数学上"凑"出了正确的公式。

现在我们知道，黑体辐射腔中的电磁场模式本质上就是谐振子。每个模式的能量量子化是谐振子能级的直接结果。普朗克无意中触及了量子力学的核心。

graph TD
    subgraph 从普朗克到谐振子
        A[普朗克假设] --> B["E = nℏω"]
        B --> C[解释黑体辐射]
        D[薛定谔谐振子] --> E[精确导出]
        E --> F["Eₙ = (n+½)ℏω"]
        F --> G["零点能 ℏω/2"]
        G --> H[真空不空]
    end
    
    style F fill:#e8f5e9
    style H fill:#fff3e0

深入：谐振子的级数解法

Griffiths 在书中展示了完整的级数解法。这里我们概述核心步骤：

第一步：无量纲化变量

令 $\xi = \sqrt{\frac{m\omega}{\hbar}} x$，方程变为：

$\frac{d^2\psi}{d\xi^2} = \left(\xi^2 - K\right)\psi$

其中 $K = \frac{2E}{\hbar\omega}$。

第二步：渐近行为分析

当 $\xi \to \infty$，方程近似为 $\frac{d^2\psi}{d\xi^2} \approx \xi^2\psi$，解为 $\psi \sim e^{\pm\xi^2/2}$。只有负指数 $e^{-\xi^2/2}$ 可归一化。

所以设 $\psi(\xi) = h(\xi) e^{-\xi^2/2}$，其中 $h(\xi)$ 是一个"温和"的多项式。

第三步：递推关系

代入后得到 $h(\xi)$ 满足 Hermite 方程：

$\frac{d^2h}{d\xi^2} - 2\xi\frac{dh}{d\xi} + (K-1)h = 0$

用幂级数 $h(\xi) = \sum_{j=0}^{\infty} a_j \xi^j$ 代入，得到递推关系：

$a_{j+2} = \frac{2j + 1 - K}{(j+1)(j+2)} a_j$

第四步：截断条件

如果级数不截断，$h(\xi)$ 在 $\xi \to \infty$ 时以 $e^{\xi^2}$ 增长，导致 $\psi$ 不可归一化。因此级数必须在某一项截断——即存在某个 $n$ 使得 $a_{n+2} = 0$。

这要求：

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

结论：能级量子化不是假设，而是可归一化条件的数学必然结果。

---

深入：升降算符法

级数解法虽然直接，但代数上很繁琐。升降算符（Ladder Operators）提供了一条优雅的捷径。

定义：

$\hat{a}_+ = \frac{1}{\sqrt{2m\hbar\omega}}(-i\hat{p} + m\omega\hat{x}) = \frac{1}{\sqrt{2}}\left(\xi - \frac{d}{d\xi}\right)$

$\hat{a}_- = \frac{1}{\sqrt{2m\hbar\omega}}(i\hat{p} + m\omega\hat{x}) = \frac{1}{\sqrt{2}}\left(\xi + \frac{d}{d\xi}\right)$

关键性质：

$\hat{H} = \left(\hat{a}_+\hat{a}_- + \frac{1}{2}\right)\hbar\omega = \left(\hat{a}_-\hat{a}_+ - \frac{1}{2}\right)\hbar\omega$

升降算符的作用：如果 $\psi_n$ 是能量为 $E_n$ 的本征态，那么：

$\hat{a}_+\psi_n \propto \psi_{n+1}, \quad \hat{a}_-\psi_n \propto \psi_{n-1}$

即 $\hat{a}_+$ 把能量"升高"一个 $\hbar\omega$，$\hat{a}_-$ 把能量"降低"一个 $\hbar\omega$。

为什么不能无限降低？

因为能量必须大于零（谐振子势能为正，总能量不可能为负）。所以存在一个"最低台阶"——基态 $\psi_0$，满足：

$\hat{a}_-\psi_0 = 0$

解这个方程：

$\frac{1}{\sqrt{2}}\left(\xi + \frac{d}{d\xi}\right)\psi_0 = 0 \Rightarrow \frac{d\psi_0}{d\xi} = -\xi\psi_0 \Rightarrow \psi_0 \propto e^{-\xi^2/2}$

这正是高斯波包！基态能量：

$E_0 = \frac{1}{2}\hbar\omega$

然后每用一次 $\hat{a}_+$，能量增加 $\hbar\omega$：

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

物理直觉：升降算符就像量子楼梯。你能往上爬（$\hat{a}_+$），也能往下走（$\hat{a}_-$），但不能跌到地面以下——地面就是零点能 $\frac{1}{2}\hbar\omega$。这个"地面"不是空无一物，而是量子涨落的最低限度。

---

数值例题：谐振子的零点能与光子能量

问题：一个质量 $m = 10^{-26}$ kg（典型分子质量）的粒子，在谐振子势中振动，角频率 $\omega = 10^{14}$ rad/s（红外光频率）。计算：

1. 零点能 $E_0$
2. 相邻能级差 $\Delta E$
3. 对应的跃迁光子的波长

解答：
第一步：计算零点能

$E_0 = \frac{1}{2}\hbar\omega = \frac{1}{2}(1.055 \times 10^{-34})(10^{14}) = 5.28 \times 10^{-21} \text{ J}$

转换为 eV：

$E_0 = \frac{5.28 \times 10^{-21}}{1.602 \times 10^{-19}} \approx 0.033 \text{ eV}$

第二步：能级差

$\Delta E = \hbar\omega = 2E_0 \approx 0.066 \text{ eV}$

第三步：光子波长

\lambda = \frac{hc}{\Delta E} = \frac{(4.136 \times 10^{-15})(2.998 \times 10^8)}{0.066} \approx 1.88 \times 10^{-5} \text{ m} = 18.8 \text{ μm}

结论：这是典型的红外振动光谱。分子振动的能级差在红外波段，这正是红外光谱学能探测分子振动的原因。零点能约0.033 eV——即使在绝对零度，分子也不能停止振动。

---

2.5 自由粒子与δ函数势

2.5.1 自由粒子

当 $V(x) = 0$ everywhere，定态薛定谔方程的解是平面波：

$\psi(x) = Ae^{ikx} + Be^{-ikx}$

对应能量 $E = \frac{\hbar^2 k^2}{2m}$。

但自由粒子有个微妙的问题：平面波不可归一化（$|\psi|^2 = |A|^2$ 在全空间积分发散）。这意味着自由粒子不存在真正的定态。Griffiths在第2章末尾用波包的概念解决了这个问题——定态只是数学工具，真实的自由粒子是波包。

物理直觉：为什么自由粒子没有定态？

定态意味着概率分布不随时间变化。但自由粒子以恒定速度运动——它"现在在这里，下一秒在那里"。一个真正的"不动"的自由粒子需要同时以所有速度向所有方向运动，这对应的就是平面波（不可归一化）。

真实的粒子是波包——许多平面波的叠加，位置局域化，速度近似确定。波包会随时间扩散（展宽），这正是因为它由不同动量的分量组成，每个分量以不同速度运动。

graph TD
    subgraph 自由粒子的描述
        A[定态平面波] --> B["动量精确确定
Δp = 0"]
        B --> C["位置完全不确定
Δx → ∞"]
        D[波包] --> E["动量近似确定
Δp 小"]
        E --> F["位置近似确定
Δx 小"]
        F --> G["满足不确定性原理
ΔxΔp ≥ ℏ/2"]
    end
    
    style C fill:#ffebee
    style G fill:#e8f5e9

2.5.2 δ函数势

势能 $V(x) = -\alpha\delta(x)$，一个无限尖锐的吸引势。

束缚态（$E < 0$）：只有一个！

$\psi(x) = \frac{\sqrt{m\alpha}}{\hbar}e^{-m\alpha|x|/\hbar^2}, \quad E = -\frac{m\alpha^2}{2\hbar^2}$

散射态（$E > 0$）：计算透射系数，会发现一个有趣的现象——即使粒子能量高于势垒，透射率也不一定是100%。

例题详解：δ函数势阱的束缚态求解

问题：求解 $V(x) = -\alpha\delta(x)$ 的束缚态。

分步推导：

第一步：分区求解

对于 $x \neq 0$，$V(x) = 0$，方程为：

$-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} = E\psi$

对于束缚态 $E < 0$，令 $\kappa = \sqrt{-2mE}/\hbar$，解为指数函数：

$\psi(x) = \begin{cases} Ae^{-\kappa x}, & x > 0 \\ Ae^{\kappa x}, & x < 0 \end{cases}$

（对称性要求 $\psi(-x) = \psi(x)$）

第二步：边界条件

在 $x = 0$ 处，波函数连续：$\psi(0^+) = \psi(0^-) = A$

导数不连续（因为有δ函数势）。对薛定谔方程从 $-\epsilon$ 到 $+\epsilon$ 积分：

$-\frac{\hbar^2}{2m}\left[\frac{d\psi}{dx}\right]_{0^-}^{0^+} - \alpha\psi(0) = 0$

第三步：计算导数跃变

$\frac{d\psi}{dx}\bigg|_{0^+} = -A\kappa, \quad \frac{d\psi}{dx}\bigg|_{0^-} = A\kappa$

$\left[\frac{d\psi}{dx}\right]_{0^-}^{0^+} = -A\kappa - A\kappa = -2A\kappa$

代入边界条件：

$-\frac{\hbar^2}{2m}(-2A\kappa) = \alpha A$

$\frac{\hbar^2 \kappa}{m} = \alpha$

第四步：解出能量

$\kappa = \frac{m\alpha}{\hbar^2}$

$E = -\frac{\hbar^2 \kappa^2}{2m} = -\frac{m\alpha^2}{2\hbar^2}$

第五步：归一化

$\int_{-\infty}^{\infty} |\psi|^2 dx = 2|A|^2 \int_0^{\infty} e^{-2\kappa x} dx = \frac{|A|^2}{\kappa} = 1$

$A = \sqrt{\kappa} = \sqrt{\frac{m\alpha}{\hbar^2}}$

结论：δ函数势阱只有一个束缚态，波函数指数衰减，能量与 $\alpha^2$ 成正比。

graph TD
    subgraph δ势阱的特征
        A[无限尖锐的势阱] --> B[只有一个束缚态]
        B --> C[指数衰减波函数]
        C --> D["无论α多小
总有束缚态"]
        D --> E[一维势的普遍性定理]
    end
    
    style B fill:#e8f5e9
    style D fill:#fff3e0

---

2.6 有限方势阱与隧穿效应

2.6.1 阱外有波！

与无限深势阱不同，有限深势阱 ($V(x) = -V_0$ for $|x| < a$，否则为0) 允许粒子以一定概率出现在阱外。

在阱外 ($|x| > a$)，$E < 0$ 时薛定谔方程的解是指数衰减：

$\psi(x) \propto e^{-\kappa|x|}$

其中 $\kappa = \sqrt{-2mE}/\hbar$。

graph TD
    subgraph 有限势阱的束缚态波函数
        A["阱内: 正弦/余弦振荡"] --> B[阱壁]
        B --> C["阱外: 指数衰减尾巴"]
    end
    
    style A fill:#e3f2fd
    style C fill:#fff3e0

这意味着：粒子有一定概率出现在经典物理不允许的区域！

物理直觉：为什么阱外有波？

经典物理中，粒子能量必须大于势能才能存在——否则动能会为负，速度为虚数。但量子力学中，波函数在"禁区"指数衰减，不代表粒子"闯入"了禁区然后被弹回，而是粒子在禁区内"试探"了一段距离后才返回。

这就像你试图翻过一座山：经典力学说要么你有足够能量翻过去，要么你根本到不了另一边。量子力学说，即使你没有足够能量，你也有一定概率"隧穿"山体到达另一边——不是因为你找到了隧道，而是因为你的"物质波"衰减得不够快，在另一侧仍有残余振幅。

深入：有限深势阱的超越方程

Griffiths 在书中用图解法求解了有限深势阱。这里我们概述对称势阱的完整求解过程：

势阱设置：$V(x) = -V_0$ for $|x| < a$，$V(x) = 0$ for $|x| > a$。

偶宇称解（$\psi(-x) = \psi(x)$）：

阱内 ($|x| < a$)：$\psi(x) = A\cos(lx)$，其中 $l = \sqrt{2m(E+V_0)}/\hbar$

阱外 ($|x| > a$)：$\psi(x) = Be^{-\kappa|x|}$，其中 $\kappa = \sqrt{-2mE}/\hbar$

边界条件（在 $x = a$ 处）：

1. 波函数连续：$A\cos(la) = Be^{-\kappa a}$
2. 导数连续：$-Al\sin(la) = -B\kappa e^{-\kappa a}$

两式相除，消去 $A, B$：

$l\tan(la) = \kappa$

令 $z = la$，$z_0 = \frac{a}{\hbar}\sqrt{2mV_0}$，则：

$\tan(z) = \sqrt{\frac{z_0^2}{z^2} - 1}$

这是一个超越方程，只能用图解或数值法求解。

图解法的物理意义：画 $y = \tan(z)$ 和 $y = \sqrt{(z_0/z)^2 - 1}$，交点就是允许的能级。

• $z_0$ 越大（势阱越深越宽），交点越多——束缚态越多
• 对于一维对称势阱，至少存在一个偶宇称束缚态，无论势阱多浅
• 第一个奇宇称态出现需要 $z_0 > \pi/2$

数值例题：一个电子在 $V_0 = 5$ eV，$a = 0.5$ nm 的有限深势阱中。估算 $z_0$ 并判断有几个束缚态。

解答：

$z_0 = \frac{a}{\hbar}\sqrt{2mV_0} = \frac{0.5 \times 10^{-9}}{1.055 \times 10^{-34}}\sqrt{2(9.11\times 10^{-31})(5 \times 1.602\times 10^{-19})}$

$\approx \frac{0.5 \times 10^{-9}}{1.055 \times 10^{-34}} \times 3.82 \times 10^{-24} \approx 1.81$

由于 $z_0 \approx 1.81 > \pi/2 \approx 1.57$，所以至少有一个偶宇称态，且可能有一个奇宇称态。

精确地说，偶宇称解的超越方程在 $z \in (0, \pi/2)$ 内总有一个解（因为 $\tan(z)$ 从0到 $\infty$，而右边从 $\infty$ 到有限值）。所以至少一个束缚态。如果 $z_0 > \pi/2$，则可能有第二个束缚态。

结论：这个势阱大约有 1-2 个束缚态。具体数目需要数值求解超越方程。

---

2.6.2 量子隧穿

如果势能是势垒而非势阱 ($V(x) = V_0 > 0$ for $|x| < a$)，粒子能量 $E < V_0$：

经典物理：粒子不可能穿过。
量子力学：波函数在势垒内指数衰减，但在另一侧仍有非零振幅——粒子有概率隧穿过去！

透射概率近似为：

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

graph LR
    A[入射波] --> B["势垒区
指数衰减"]
    B --> C["透射波
T = e^(-2κa)"]
    B --> D["反射波
R = 1-T"]
    
    style C fill:#e8f5e9
    style D fill:#ffebee

例题详解：隧穿概率的数值估算

问题：一个电子（质量 $9.11\times 10^{-31}$ kg）面对一个高度 $V_0 = 5$ eV、宽度 $a = 1$ nm的势垒，能量 $E = 4$ eV。估算隧穿概率。如果换成质子（质量是电子的1836倍），概率会如何变化？

分步推导：

第一步：计算参数

$V_0 - E = 1 \text{ eV} = 1.602 \times 10^{-19} \text{ J}$

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

$= \frac{\sqrt{2.92 \times 10^{-49}}}{1.055 \times 10^{-34}}$

$\approx \frac{5.40 \times 10^{-25}}{1.055 \times 10^{-34}}$

$\approx 5.12 \times 10^9 \text{ m}^{-1}$

第二步：计算指数

$2\kappa a = 2(5.12 \times 10^9)(10^{-9}) = 10.24$

第三步：计算隧穿概率

$T = e^{-10.24} \approx 3.6 \times 10^{-5}$

大约每10万个电子中有3-4个能隧穿过去。

第四步：质子的情况

质子质量 $m_p = 1836 m_e$。由于 $\kappa \propto \sqrt{m}$：

$\kappa_p = \sqrt{1836} \kappa_e \approx 42.8 \kappa_e$

$2\kappa_p a = 42.8 \times 10.24 \approx 438$

$T_p = e^{-438} \approx 10^{-190}$

物理意义：质子的隧穿概率小到可以忽略不计。这就是为什么核聚变需要极高的温度和压力——质子必须"热"到足以经典地克服库仑势垒，或者等待极长的时间才能隧穿（太阳核心的质子平均等待约 $10^9$ 年才隧穿一次，但由于粒子数庞大，总反应速率仍然可观）。

graph TD
    subgraph 隧穿的现实意义
        A[扫描隧道显微镜] --> B["电子隧穿过真空
探测表面形貌"]
        C[太阳核聚变] --> D[质子隧穿库仑势垒]
        D --> E["极低概率 × 极大数量
= 可观反应率"]
        F[闪存存储器] --> G["电子隧穿氧化层
写入/擦除数据"]
    end
    
    style B fill:#e8f5e9
    style E fill:#e8f5e9
    style G fill:#e8f5e9

---

2.7 本章总结

graph TD
    A[第2章核心] --> B[定态薛定谔方程]
    A --> C[无限深势阱]
    A --> D[谐振子]
    A --> E["自由粒子/δ势/势垒"]
    
    B -->|求解| F["Hψ = Eψ
分离变量法"]
    C -->|结果| G["离散能级
Eₙ ∝ n²"]
    D -->|结果| H["等间距能级
Eₙ ∝ n+½"]
    E -->|结果| I["隧穿效应
经典禁区有概率"]
    
    style G fill:#e3f2fd
    style H fill:#e8f5e9
    style I fill:#fff3e0

带走的三句话：

1. 束缚导致量子化：粒子被限制时，能量只能取离散值。
2. 零点能不可避免：被束缚的粒子不能静止，这是不确定性原理的要求。
3. 隧穿是真实的：粒子可以穿过经典物理认为"不可能"的区域。

与其他章节的联系

graph LR
    A[第2章] --> B["第1章
时间依赖方程的特例"]
    A --> C["第3章
本征值问题的形式化"]
    A --> D["第4章
三维求解技巧"]
    A --> E["第6章
微扰理论修正"]
    B --> F[定态是时间依赖方程的可分离解]
    C --> G["本征态 = 希尔伯特空间的基"]
    D --> H["三维势阱、氢原子径向方程"]
    E --> I[近似求解复杂势场]
    
    style A fill:#e3f2fd
    style B fill:#e8f5e9
    style C fill:#e8f5e9
    style D fill:#e8f5e9
    style E fill:#fff3e0

"量子力学最深刻的真理之一：限制产生离散性。" —— Griffiths

---

2.8 练习与思考

1. 势阱尺度：证明无限深方势阱中，基态能量可以写成 $E_1 = \frac{h^2}{8ma^2}$（用 $h$ 而非 $\hbar$）。这说明量子效应在质量小、空间窄时最显著。估算一个质量 $m = 1$ g 的宏观粒子在 $a = 1$ cm 势阱中的基态能量，并与室温热能 $k_B T \approx 0.025$ eV 比较。

2. 节点定理：无限深势阱的第 $n$ 个本征态有 $n-1$ 个节点。证明波函数在相邻节点之间的区域符号相反，并解释为什么这保证了本征态的正交性。（提示：利用 $\int_0^a \psi_m \psi_n dx = 0$ for $m \neq n$）

3. 隧穿估算：一个电子（质量 $9.11\times 10^{-31}$ kg）面对一个高度 $V_0 = 5$ eV、宽度 $a = 1$ nm的势垒，能量 $E = 4$ eV。精确计算（而非近似）透射系数 $T$，并与近似公式 $T \approx e^{-2\kappa a}$ 比较。如果势垒宽度加倍到2 nm，透射概率如何变化？