一维势阱是量子力学中最简单最基本的模型,在量子力学教学和科研过程中起着基础性作用[1]。其理论结果在许多实际系统中也得到了很好的应用,比如在低维量子系统(如量子点、量子面、量子线等)中[2]。但作为量子力学中的基本方程,薛定谔方程[3]的求解却并不简单,在《量子力学》课程的教学过程中,除了一维无限深势阱、一维谐振子势等特殊一维情形有解析解外,一般情形很难求解。即使对于处在一维有限深势阱中的运动的粒子,当其处于束缚态时,由于确定其能级的是超越方法,无法具体给出它们的能级的解析表达式和归一化波函数[4]。
薛定谔方程的求解可分为解析法和数值法两种,解析法主要有WKB法、变分法等,数值法有打靶法[5]和有限元(FEM)法等。作为一种有益的探索,本文以处理多层光波导问题的转移矩阵(TM)法为算法,求解粒子在不同一维有限深势阱中运动的能量本征值,并利用Matlab作为计算工具,绘制出对应的波函数图像。采用TM方法计算,通过不同势阱的设置,可以得到定性的准解析规律性结论,为解释实验现象提供重要的理论依据。而且,在《量子力学》课程的学习中,学生可借助不同势阱的计算分析将理论学习进行拓展,更直观更透彻地理解一维势阱相关方面的知识点,加深对相关物理概念的理解。
1 转移矩阵法
1.1 算法应用依据
转移矩阵法起源于光学,用于计算多层薄膜介质的反射率和透射率。在光波导技术中,转移矩阵法用来建立不同介质层间的场分布联系。而经过对比发现:不含时的薛定谔方程与平板波导的波动方程在形式上是一致的;同时一维势阱模型与平板波导模型相类似,因而转移矩阵法也能够用于解决一维有限深势阱问题[6]。 平板波导的波动方程为[7]
![d243225dd006cbb2fa8711d6f3e14db0.png](https://i-blog.csdnimg.cn/blog_migrate/aafdb19cf78a3fd5d74bca4362da5126.jpeg)
(1) 其中k0=2π/λ,是光在真空中的传播常数,Ey为y方向的电场,β为电磁场沿z方向的传播常数,nj对应于波导中不同层的折射率,j=0,1,2。 不含时的一维薛定谔方程形式为
![f890b902c2f8cb9538c3ce22098dfc6b.png](https://i-blog.csdnimg.cn/blog_migrate/58d544f3cb5bd9f5894ec7ee099746b2.jpeg)
其中V(x)表示势能,m为粒子质量,E和ψ分别代表粒子的能量本征值和本征函数,上式经过变形,有:
![c4d2c2e6030875c0e391aa03d965dd8f.png](https://i-blog.csdnimg.cn/blog_migrate/9c1b7893e231b577cdcd5e09bc78e87f.jpeg)
(2) 对比方程(2)与方程(1)即可发现:这两个方程在形式上是一致的,方程(2)中的本征函数ψ即对应方程(1)中的电场强度Ey,因此两个方程可采用相同的方法求解。 按照波动光学理论,光波是电磁横波,光波在空间任意位置的电磁场强度及所在介质性能之间的联系都是通过Maxwell方程组和物质方程建立,电磁波在两种介质形成的界面上反射和透射时的振幅反射因数和透射因数均可以由菲涅耳公式确定。若是一个多界面的薄膜系统,则根据以下两点[8]:(1)在每一界面处运用电磁场边界条件,将同一界面两侧的场分布联系起来;(2)利用与电磁场传播相伴随的相位差,将同一膜层上下两界面内侧的场分布联系起来,可将多界面系统看作是入射介质与薄膜、基底形成的等效介质之间的界面。将以上思路应用到任何一个复杂的薄膜系统,此时系统中光反射率和透射率问题,都可以通过其等效界面对应的等效介质进行计算。等效介质的光谱特性可以采用一个特征矩阵来表示,该矩阵也称之为转移矩阵。 非对称平板波导的模型如图1所示。通常来说折射率有n1>n0>n2,导波光可被约束在导波层中进行传播。将导波光类比为势场中运动的粒子,粒子会被约束在势能较低处运动,则导波层可看作势阱中势能低的地方,衬底层和覆盖层可看作势能相对较高的地方。从这一角度看,势阱中粒子运动的模型和波导中导波光的传播模型是一致的。基于以上分析,可将转移矩阵法应用到不含时薛定谔方程的求解。而在实际应用中,例如超晶格量子阱,虽然由能隙不同的材料组成可以形成复杂的势阱,但系统的基本特征完全可以通过有限深势阱问题获得很好的理解[2]。
![3fc0580073740c8517650870006f3272.png](https://i-blog.csdnimg.cn/blog_migrate/c695f184489fada0a7d3770b66af9a5c.jpeg)
图1 非对称平板波导模型
1.2 建模与求解
对方程(1)的求解,首先选取两个特解E1(x)和E2(x),使之满足E1(0)=E′2(0)=1和E′1(0)=E2(0)=0。则方程(1)的通解可设为E1(x)和E2(x)的线性叠加,即 Ey(x)=C1E1(x)+C2E2(x) (3) 在区间(0,d)界面上,根据(3)式可给出场分布及其导数,