4. 有限深势阱的超越方程建立

无限深势阱的能级公式干净利落，但它假设阱壁无限高，波函数在边界处被强制归零。真实物理系统——半导体量子阱、金属纳米颗粒、原子核模型——阱壁高度都是有限的，波函数会渗入势垒区，边界条件从“强制为零”变成“连续且光滑”。这一变化使能级不再有解析闭式解，而必须由一个超越方程数值确定。本节就围绕这个超越方程的建立过程展开。

一、分区与坐标约定

取一维有限深方势阱，宽度为 2a，关于原点对称放置，势能定义为：

V(x) = 0      (|x| < a，阱内)
V(x) = V₀     (|x| ≥ a，势垒区，V₀ > 0)

约定粒子质量 m，总能量 E 满足 0 < E < V₀，即束缚态。由于势能关于 x = 0 对称，V(-x) = V(x)，哈密顿量具有空间反演对称性，这为后面按宇称分类解埋下伏笔。全文坐标原点取在势阱中心，阱壁位于 x = ±a，这一约定在推导和图示中保持一致。

二、分区域写出通解

三个区域分别写定态薛定谔方程。阱内 |x| < a：

-ħ²/(2m) ψ'' = Eψ   →   ψ'' + k²ψ = 0,   k = √(2mE)/ħ

通解为振荡形式 ψ(x) = A sin(kx) + B cos(kx)。阱外 |x| > a，令 κ = √(2m(V₀ - E))/ħ，方程变为 ψ'' - κ²ψ = 0，通解为指数形式。为保证波函数在 x → ±∞ 时平方可积，必须舍去发散项：

ψ(x) = C e^{-κx}     (x > a)
ψ(x) = D e^{κx}      (x < -a)

这里 k 与 κ 并非独立，二者通过能量联系起来：k² + κ² = 2mV₀/ħ²。这个关系式是后面把两个未知量化为一个的关键。

三、边界连续性与宇称分类

在 x = a 处要求 ψ 和 ψ' 连续，得到两个方程；在 x = -a 处同理再得两个方程。四个方程联立看似复杂，但利用势的对称性可以大幅简化：解必然可以取为偶宇称（ψ(-x) = ψ(x)）或奇宇称（ψ(-x) = -ψ(x)），两类解互不混合。

偶宇称解：阱内取 ψ = B cos(kx)，只需在 x = a 匹配：

B cos(ka) = C e^{-κa}
-kB sin(ka) = -κC e^{-κa}

两式相除消去常数，得 k tan(ka) = κ。

奇宇称解：阱内取 ψ = A sin(kx)，同样在 x = a 匹配：

A sin(ka) = C e^{-κa}
kA cos(ka) = -κC e^{-κa}

相除得 -k cot(ka) = κ，即 k cot(ka) = -κ。

这两个方程就是有限深方势阱的超越方程。所谓“超越”，指未知量 k 同时出现在代数因子和三角函数（或余切）的自变量中，无法用初等函数显式解出，只能数值求根。

四、化为无量纲形式

直接对 k tan(ka) = κ 求根不便，因为 k 与 κ 受约束。引入无量纲量 ξ = ka，η = κa，则约束变为：

ξ² + η² = 2mV₀a²/ħ² ≡ R²

其中 R = a√(2mV₀)/ħ 是刻画势阱“强度”的唯一参数。超越方程随之化为：

偶宇称：η = ξ tan ξ
奇宇称：η = -ξ cot ξ

几何意义非常清晰：解就是圆 ξ² + η² = R²（第一象限部分）与曲线 η = ξ tan ξ 或 η = -ξ cot ξ 的交点。每个交点对应一个束缚态能级，能量由 E = ħ²ξ²/(2ma²) 反算。

五、示例：数值求解一个具体势阱

设某半导体量子阱中电子有效质量 m* = 0.067 m₀，阱宽 2a = 10 nm（即 a = 5 nm），阱深 V₀ = 0.3 eV。先算强度参数：

R = a√(2m*V₀)/ħ
  ≈ 5×10⁻⁹ × √(2×0.067×9.11×10⁻³¹×0.3×1.6×10⁻¹⁹) / 1.055×10⁻³⁴
  ≈ 2.6

取 R ≈ 2.6，在 ξ ∈ (0, 2.6) 内找交点。偶宇称曲线 ξ tan ξ 在 ξ = π/2 ≈ 1.571 处发散，故第一个偶宇称根位于 (0, 1.571)；奇宇称曲线 -ξ cot ξ 在 ξ = π ≈ 3.14 处才发散，而圆半径只有 2.6，故奇宇称根位于 (1.571, 2.6) 区间内。用二分法或牛顿法可求得（以下为示例数值，供教学演示）：

可见该阱只容纳两个束缚态，且基态为偶宇称、第一激发态为奇宇称，能级交替出现。若把 V₀ 增大使 R 超过 π，就会出现第二个偶宇称解，依此类推。

六、与无限深势阱的对照

无限深势阱中 ka = nπ/2，能级严格按 n² 增长，波函数在阱壁处为零。有限深势阱中，ξ 根总略小于对应的 nπ/2，因为波函数渗入势垒使“有效宽度”变大，能级整体下移；且束缚态数目有限，R 越小态越少。当 V₀ → ∞ 时 R → ∞，η → ∞，超越方程退化为 tan ξ = ∞ 或 cot ξ = ∞，即 ξ = nπ/2，回到无限深结果。这一极限检验是验证推导正确性的重要手段。