WKB近似模拟器 返回
量子力学

WKB近似模拟器 — 量子隧穿概率

一个体验"量子隧穿效应"的工具——粒子在经典力学中无法跨越的势垒,在量子力学中却能穿透。改变粒子能量、势垒高度与宽度、质量等参数,实时查看WKB近似计算的透射概率 T 与势垒内衰减常数 κ。

参数设置
粒子能量 E
eV
入射粒子的动能
势垒高度 V0
eV
矩形势垒的势能
势垒宽度 L
nm
粒子所要穿透的势垒厚度
粒子质量 m
mₑ
以电子质量为1的单位。质量越重隧穿越困难
计算结果
隧穿概率 T
衰减常数 κ (nm⁻¹)
透射深度 1/κ (nm)
E/V0 比值
势差 V0−E (eV)
经典判定
势垒与波函数 ψ(x)

蓝线为能量 E,灰色矩形为势垒 V0。波函数 ψ(x) 在势垒前振荡,势垒内部指数衰减,通过后作为透射波呈现小振幅振荡。

隧穿概率 vs 势垒宽度 L
隧穿概率 vs 粒子能量 E
理论与主要公式

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

WKB透射概率通用公式。κ 是势垒内波函数的衰减率,x₁·x₂ 是经典转折点(E=V的位置)。积分值越大,透射概率越小。

$$T = e^{-2\kappa L}$$

对于高度V0、宽度L的矩形势垒,κ为常数,积分简化为κL,透射概率变成这个简洁的指数函数。透射深度1/κ是波函数衰减到1/e的距离。

WKB近似在势能相对波长变化缓慢时有效。矩形势垒严格上属于急变势,但作为学习衰减本质的教学材料被广泛使用。

WKB近似原理

🙋
"量子隧穿效应"是说粒子能穿过墙壁,对吧?这真的在现实中会发生吗?太科幻了。
🎓
是的,确实发生。而且是真实的物理现象。简单地说,古典力学中你把球往墙扔,如果能量不足就100%反弹回来。但电子这样的微观粒子既有波的性质,波函数在势垒内部不会突然消失,而是呈指数衰减地渗透。如果势垒足够薄,衰减还没完全消失它就已经穿到另一边去了——这就是隧穿效应。
🙋
明白了。那么这个"穿过的量"怎么计算呢?薛定谔方程好像很复杂啊。
🎓
这时候"WKB近似"就派上用场了。名字来自三位科学家 Wentzel、Kramers、Brillouin。与其精确求解薛定谔方程,不如利用一个巧妙的想法:在势能相对波长变化缓慢的区域,波函数可以近似为 exp(±∫k dx) 的指数形式。这样隧穿概率就能简化成 T ≈ exp(−2∫κ dx) 这样简单的积分。也叫半经典近似。
🙋
我看左边的图,当势垒宽度 L 增大时,隧穿概率 T 掉得特别快,即使在对数坐标上也是直线。
🎓
很好的观察!对于矩形势垒,T = exp(−2κL)。指数项包含了L,所以增大L时T呈指数下降。对数坐标上呈直线正是这个原因。默认参数下T ≈ 2×10⁻⁸,大概是说"被撞一亿次才成功一次",但注意不是零!如果把势垒宽度减半,T会增加四倍左右。这就是为什么改变势垒厚度会产生巨大效应。
🙋
质量 m 变重时也很难隧穿。为什么呢?
🎓
因为衰减常数 κ 与 √(m(V0−E)) 成正比。质量越大,κ就越大,波函数在势垒内衰减得越快。所以电子隧穿比较容易,但质子(电子质量的约1836倍)就很难。α衰变之所以具有极长的半衰期,正是因为α粒子很重。如果把能量E靠近V0,那么(V0−E)变小,κ下降,T就会急剧增大。下面的能量图可以看到这个效果。
🙋
我把能量 E 调成大于 V0 了,判定词变成"经典上也能透过"。这就不是隧穿了,对吧?
🎓
正确!当 E ≥ V0 时,粒子能从经典的角度"翻过"势垒,已经不是隧穿效应了。此时(V0−E)为负,√里头不成立,κ=0,隧穿公式 T=exp(−2κL) 就不适用。本工具用判定语句明确区分了两种情况。隧穿效应只在 E < V0 的"经典上绝对无法通过"的情况下才有意义。

常见问题

WKB近似(Wentzel–Kramers–Brillouin近似)是薛定谔方程的半经典近似求解方法。在势能相对于德布罗意波长变化缓慢的区域,波函数可近似为指数形式 exp(±∫k dx)。利用这一近似,即使对于隧穿效应这样原本很难精确求解的问题,透射概率也可简化为积分 T ≈ exp(−2∫κ dx)。本工具对矩形势垒应用了WKB公式。
在经典力学中,能量E低于势垒高度V0的粒子绝对无法穿过势垒,必然100%反射。但在量子力学中,粒子具有波的性质,波函数在势垒内部不会突然变为零,而是指数衰减地渗入势垒。如果势垒宽度有限,衰减还未完全消失时波函数就已到达另一侧,这部分就能够逃逸——这就是隧穿效应。透射概率取决于势垒内部的衰减量。
矩形势垒的WKB透射概率为 T = exp(−2κL),其中κ是势垒内的衰减常数,L是势垒宽度。κ = √(2m(V0−E))/ℏ,因此T对势垒宽度L呈指数递减,对粒子质量m与势差(V0−E)的平方根也呈指数递减。也就是说,将势垒宽度减半或降低高度,T会大幅增加。反之,对于重粒子或厚而高的势垒,T会急剧减小。
WKB近似适用于势能相对粒子波长变化足够缓慢的情况。反之,边界陡峭的矩形势垒或势能急剧变化的点(经典转折点附近)会产生较大误差。本工具的矩形势垒在严格意义上属于急变势,但作为理解衰减本质的教学材料,仍采用WKB公式 T=exp(−2κL)。在α衰变或扫描隧道显微镜这样的实际缓变势问题中,WKB近似给出很高的精度。

实际应用

扫描隧道显微镜(STM):探针与样品表面间的真空间隙极小,电子可通过隧穿在间隙中流动。隧穿电流遵从 T = exp(−2κL),间隙L改变0.1nm时电流就变化10倍。这种极端灵敏性正是STM实现原子级分辨率的基础。试试用本工具改变L参数,你会看到T瞬间变化多个数量级。

α衰变与原子核物理:放射性原子核释放的α粒子是通过隧穿库仑势垒逃逸的。1928年,Gamow正是用WKB近似解释了α衰变半衰期。由于隧穿概率对势垒呈指数关系,半衰期可以跨越20多个数量级,从10⁻⁷秒到10¹⁷年。这就是Geiger-Nuttal规律的物理基础。

半导体器件:隧穿二极管、闪存的擦写、超薄栅氧化膜的漏电流等现代电子学现象都与隧穿息息相关。随着器件微型化,栅绝缘膜薄至几纳米,原本不想要的隧穿漏电流成为芯片设计不得不考虑的因素。WKB近似仍是快速估算这种漏电流的第一近似。

核聚变与恒星能量:太阳内部即使温度极高,经典上两个质子还是撞不出足够能量融合。质子必须通过隧穿库仑势垒才能发生核聚变,恒星因此才会发光。隧穿概率对温度极其敏感,中心温度微小的变化就会导致融合速率剧烈改变。WKB近似用于计算这个"Gamov因子"。

常见误解与注意

最常见的误解是"隧穿粒子在势垒内部借了能量"。实际上隧穿过程中粒子能量守恒,势垒内的波函数与入射时拥有同样的能量E。"不确定性原理允许粒子临时借能量翻过势垒"这种说法看上去直观但不严格。真实情况是波函数在势垒内从振荡变成指数衰减(衰减波),这种渗透到另一侧,并无能量借贷。本工具的波函数图也清晰显示势垒内部的纯衰减特征。

第二个误解是"WKB近似对矩形势垒完全精确"。WKB近似的前提是"势能相对波长缓慢变化",而矩形势垒的边界处势能跳跃不连续,违反了这个前提。教学上常用它是因为 T=exp(−2κL) 形式简洁,但严格的矩形势垒透射率需要加前置因子,尤其在E接近V0或势垒很薄时,WKB的指数项可能有数倍偏差。在α衰变或缓变势的实问题上,WKB才真正发挥威力。

第三个误解是"隧穿概率那么小可以忽略"。本工具默认参数T ≈ 2×10⁻⁸看起来极小,但在试验次数极多的场景下就不容忽视。栅氧化膜中每秒有天文数字的碰撞尝试,10⁻⁸的概率叠加后就成为可观的漏电流,足以影响芯片设计。α衰变也是"单个粒子概率微小,但原子数目庞大"。应该看期望值,不是概率绝对值。

使用指南

  1. 设置粒子能量eNum(eV)与扫描范围eRange,输入势垒高度vNum(eV)与宽度wNum(nm)
  2. 选择质量mNum(电子质量me单位),模拟器自动计算WKB积分∫κ(x)dx
  3. 实时读取衰减常数κ=√(2m(V0−E))/ℏ与透射概率T≈exp(−2κW)
  4. 在E/V0比值0.3至0.9的范围扫描参数,观察隧穿概率的指数递减现象

具体计算案例

α粒子(m=4 me、E=5 MeV)穿过金核库仑势垒(V0=8 MeV、W=2 fm)时,κ≈0.48 fm⁻¹,透射概率T≈10⁻⁴⁰。相比之下,电子(m=1 me、E=2 eV)穿过有限势垒(V0=5 eV、W=1 nm)时,κ≈1.03 nm⁻¹,T≈0.12,体现出有意义的隧穿概率,适用于半导体隧穿接点与STM原理。

实务中的注意事项