1. 课程导入与量子隧穿现象概览

在经典力学中，一个能量为 E 的粒子遇到高度为 V₀ 的势垒时，只要 E < V₀，它就会被完全反射，绝不可能出现在势垒另一侧。然而在量子力学框架下，粒子的状态由波函数描述，波函数在势垒内部并不立即归零，而是以指数形式衰减；只要势垒宽度有限，波函数就能在另一侧以非零幅度"泄漏"出去。这种粒子穿过经典禁区、在另一侧被探测到的现象，就是量子隧穿效应。

隧穿并非纯粹的数学游戏。半导体中的隧道二极管、闪存单元的电荷注入与擦除、原子核的 α 衰变、扫描隧道显微镜（STM）的成像机制，都直接依赖这一效应。本课程的目标，是把"波函数如何穿过势垒"这一抽象过程，转化为可计算、可绘图、可交互观察的数值实验：用有限差分法离散薛定谔方程，用数值方法求解定态与含时问题，再把透射系数、概率密度、波包演化等结果可视化，从而建立"公式—代码—图像—物理判断"之间的完整链条。

从一维方势垒到可计算的模型

考虑一维定态薛定谔方程

-（ħ²/2m）·d²ψ/dx² + V(x)·ψ = E·ψ

取方势垒 V(x) = V₀（0 ≤ x ≤ a），其余区域为 0。在 E < V₀ 时，势垒内解为指数型，两侧为平面波，匹配边界条件后可得透射系数

T ≈ 16·(E/V₀)·(1-E/V₀)·exp(-2·κ·a)，  κ = √(2m(V₀-E))/ħ

这个近似式揭示了隧穿最核心的参数关系：透射概率对势垒宽度 a 和衰减常数 κ 呈指数敏感。以电子为例，若 V₀ − E = 1 eV，则 κ ≈ 5.12 nm⁻¹；势垒宽度从 0.5 nm 增至 1.0 nm，指数因子 exp(−2κa) 从约 6.0×10⁻³ 降到约 3.6×10⁻⁵，透射率下降两个数量级以上。这正是 STM 能实现原子级分辨率的物理根源：针尖与样品间距每变化约 0.1 nm，隧道电流就变化约一个数量级。

数值模拟与可视化的整体路径

本课程的学习路径围绕"建模—离散—求解—可视化—验证"展开，但各环节并非简单线性排列，而是相互反馈的闭环。下图给出课程的整体技术路线与各环节之间的依赖关系：

图中上半部分是计算主链：物理建模确定势垒与粒子参数，空间离散把连续方程变成矩阵或差分格式，数值求解得到波函数或时间演化，可视化把结果转成可读图像。下半部分则是这条主链服务的对象——真实的物理现象与工程应用。中间的验证环节不是可有可无的收尾，而是贯穿始终的质量控制：任何一条透射曲线、任何一帧波包动画，都应当能与解析近似、概率守恒或网格收敛性对上号，否则图像再漂亮也不可信。

本节需要建立的三个基本认识

第一，隧穿是波动性的必然结果，不是"粒子偷偷溜过去"。 波函数在势垒内按 exp(−κx) 衰减，衰减长度 1/κ 决定了势垒"看起来有多厚"。理解这一点，才能在数值结果中正确解读"势垒内概率密度不为零"这一现象，而不会误以为是计算错误。

第二，数值模拟的关键在于离散参数的选择。 网格步长 Δx 必须远小于波函数的最小特征长度（势垒内衰减长度 1/κ 和德布罗意波长都要考虑）。若 Δx 取得过大，指数衰减段只剩两三个格点，透射系数会出现严重误差。一般建议每个衰减长度内至少布置 10～20 个格点，并通过减半 Δx 观察 T 是否收敛来确认。

第三，可视化不是装饰，而是判断工具。 把 |ψ|² 沿 x 的分布画出来，能直观看到势垒内的指数衰减段、势垒右侧的透射波幅度，以及入射侧反射波形成的驻波条纹。把 T 对 E 或对 a 画成曲线，能立刻看出指数敏感区间在哪里。这些图像本身就是物理判读的依据。

本课程将交付什么

完成全部章节后，你将拥有一套可复用的隧穿模拟代码：能够输入势垒高度、宽度、粒子能量与质量，输出透射系数曲线、势垒内外的波函数分布图，以及含时波包穿过势垒的动画帧序列。更重要的是，你将掌握一套判断数值结果是否可信的方法——用解析近似做交叉验证，用概率守恒做全局检查，用网格收敛性做精度评估。这套方法不仅适用于方势垒，也能迁移到双势垒共振、任意形状势垒等更复杂的场景。

接下来的章节将从最基础的一维方势垒定态问题入手，逐步引入矩阵求解、含时演化与交互式可视化。建议在阅读时同步打开代码环境，把每一节的参数亲手改一改、图亲手画一画——隧穿的指数敏感性，只有亲眼看到透射率随势垒宽度骤降几个数量级，才会真正变成你自己的直觉。