MATLAB在量子力学中的应用二维方阱中粒子的运动.pdfVIP

MATLAB 在量子力学中的应用 ——二维方阱中粒子的运动 夏 静 (华中师范大学 物理科学与技术学院 武汉 430079) 摘要： 通过分离变量法得到二维无限势深方阱的解，然后采用强大的 matlab7.0 绘制波函数图形，建立直观形象的理论模型,并附上程序代码。 关键字：二维方阱；分离变量法； MATLAB ；波函数 一、引言 【1】 教科书 在薛定谔方程一章中有一个关于求一维无限深方势阱这 个习题，由此我们联想到二维情况，可以猜想，二维情况下的方阱可 以与狄拉克势阱类比。本文通过分离变量法求解二维无限势深方阱的 解，然后采用强大的matlab7.0 绘制波函数图形，这不仅有利于教师教 学，而且可以提高学生的兴趣并加深其对薛定谔方程的理解，对其有 很大的帮助。另外，本文同时附上代码，便于参考。 二、二维无限势深方阱中的波函数和能量 设无限深二维方势阱如图 1 ,势阱的长、宽分别为a 、b , 即一 个方形量子围栏。对于处在方形量子围栏中的粒子而言, 这围栏便成 【2】 了一个势垒, 它们被束缚在一个矩形势垒中,势场函数 为： 0 ， 0 x a 且 0 y b U (x , y ) = (1) ∞, x ≤ 0 或 x ≥ a 且 y ≤ 0 或 y ≥ b z y b a 图1 x −h2 2 [ ∇ +U ( x, y)]ψ( x, y) Eψ( x, y) 定态薛定谔方程： （2 ） 2m 2 2 2 ∂ ∂ 二维情况下，拉普拉斯算符表示为：∇ =+ （3 ） 2 2 ∂ x ∂ y 在阱外，即 U (x , y ) = ∞时，波函数 只能为0，否则失去意义。 ψ(x , y ) −h2 2 在阱内，即 U (x , y ) = 0 时， ∇ψ(x , y ) Eψ(x ,