量子阱能级与波函数的MATLAB实现
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
哈尔滨师范大学
毕业论文
题目量子阱能级和波函数的MATLAB实现
学生王勇
指导教师孙文军教授
年级2010级
专业物理学(通用技术)
系别物理系
学院物理与电子工程学院
哈尔滨师范大学
2013年5月
量子阱能级与波函数的MATLAB 实现
王勇
摘 要:在量子力学中,通过对求解一维多量子阱束缚态能级满足的超越方程,我们可以求出量子阱的能级表达式。另外,通过薛定谔方程求出对氢原子波函数满足的方程,在MATLAB 中通过程序的编辑,可以方便地绘制量子阱能级和波函数的空间分布图,指出MATLAB 软件是解决量子力学可视化的有效工具。 关键词:MATLAB ;量子阱能级;波函数
量子力学作为物理类学生必修的一门专业课. 一直以来以教师难教、学生难学而著称,量子力学主要研究微观尺度下粒子的行为与相互作用/,例如量子力学中的许多概念如角动量理论,波函数等抽象难懂. 一些量子力学现象如隧道效应、势垒反射等与宏观现实不相符. 同时这些理论和现象无法直接用肉眼观察,使得这些理论和现象更加晦涩难懂. 因此,在量子力学教学中引入多媒体教学.,增加动态直观的演示,是进行教学改革、也是将晦涩难懂的抽象理论化为形象生动的图像理论。而选择一种操作简单、功能强大的多媒体制作软件是进行教育改革的前提,针对量子力学计算复杂、物理量具有矩阵形式等特点. 我们最终选择了具有强大计算能力同时又有极强图像处理能力的数学软件——MATLAB 。
一、MATLAB 软件介绍
MATLAB 是美国MATH WORK 公司从1982年开始推出的一套高性能的数值计算和可视化软件, 它集数值分析、矩阵运算和图形显示于一体. 构成了一个方便的、界面友好的用户环境,MATLAB 程序设计语言结构完整. 且具有良好的移植性,易学易用.此外. 它还提供了解决各类问题的工具箱,MATLAB 已成为应用学科、计算机辅助分析、设计、仿真和教学不可缺少的软件.
MATLAB 具有强大的图形绘制能力. 为科学研究提供了极大的方便,MATLAB 可以绘制二维、三维乃至四维图形.而且能对图形进行线型、立体、色彩、光线、视角等控制,用户只须指定绘图方式. 并提供充足的绘图数据. 用很少的程序指令就可得到直观、形象的图形结果.借助于MATLAB 的数值计算和图形处理技术.我们可以方便地绘制氢原子波函数的空间分布图. 直观感受微观状态下电子行为,加深对理论、概念的理解。
二、一维多量子阱能级和氢原子波函数模型的建立
(一)一维多量子阱能级
由量子理论知,一维多量子阱满足的薛定谔方程为
n n n E V dx d m h ψψ=⎥⎦
⎤⎢⎣⎡+-2
2
22
一维多量子阱共有n 个常数位势,其高度分别用n V ,V ,,V ,V ,V 1-n 321 来标志,选
21V V 和阶跃点的坐标为零(0x 1=),1j ,V +j V 阶跃点的坐标为j x ,第j 个位势的宽
度为1=-=j j j x x a 。通常将Vn V 1与称之为外区位势,两个外区位势的小者记为max V ,而j V 称为内区位势,内区位势的最小者为j V 。
设电子处于一维多量子阱中,当电子的能量满足max min V V < 个区域的波函数有限性要求的时,可以直接写出满足n V E V max min << ()x ik A x 1 e 1,21 -=ψ (1) ()x ik x ik e A e A x 2 2 2.22,12-+=ψ (2) ()3 3 3 .23.13ikx x ik e A e A x -+=ψ (3) ··· ··· ··· ()x ik n x ik n n m e A e A x 222,22,12n ------+=ψ (4) ()x ik x ik n n n n e A e A x 111n 21,11 ------+=,ψ (5) ()x ik n n n e A x ,1=ψ (6) 试中 ) E (2k J j j V m -= (j=1,2,3, ,n) (7) 其中,j k 是实数还是虚数应视E 与j V 的数值而定,一般情况下将j k 作为复数处理。 首先,由波函数在x=1x -n 处的连接条件 )x ()x (111---=n n n n ψψ (8) 可知 1111 1,11,21,1-----=+---n n n n n n x ik n x ik n x ik n e A e A e A (9) 再利用波函数一阶导数的连接条件 n n n n n n m m x )x ()(11 11----'= 'ψψ (10) 得到 11111,111 1,21-n 1A ---------= -n n n n n n x ik n n n n n x ik n x ik e A m k m k e A e , (11) 为简捷计,令 n n n n m k m k B 11 1-n --= (12) 将式(9)加上式(11)与减去式(11), 分别得到 11 1n ,11n 1,1)1(21 -----+= n n n n x ik x ik n e A B e A (13) 11 1n ,111n ,2)1(2 1 -------=n n n n x ik n x ik e A B e A (14) 其次,由波函数在2-=n x x 处的连接条件 ()()2122----=n n n n x x ψψ (15) 可知 2121222 21,21,12n 22n 1--------------+=+n n n n n n n n x ik n x ik n x ik x ik e A e A e A e A ,, (16) 再利用波函数一阶导数的连接条件 ()()1 212 22------'= 'n n n n n n m x m x ψψ (17) 得到 () 212122221,21n 122,221A ----------------=-n n n n n n n n x ik n x ik n x ik n x ik n e A e A B e A e ,, (18) 将式(16)加上式(18)与减去式(18),分别得到 ()()2121222,221,12n 2,112 1B 121 -------------++= n n n n n n x ik n n x ik n x ik n e A B e A e A (19) ()()2221221,221,122,212 1121 -------------++-= n n n n n n x ik n n x ik n n x ik n e A B e A B e A (20) 然后,用类似的方法做下去,可以得到x=j x 处的结果,即 ()()j j j j j j x ik j j ik j x ik j A A e A 11e B 121e B 121 1,2x !,1j ,1++-++-++= (21) ()()j j j j j j ik j j ik j j x ik j A B A B e A x 1,2x 1,1,211e 12 1e 121 ++-++-++-= (22)