谐振子模型在量子系统本征值计算中的应用.docVIP

谐振子模型在量子系统本征值计算中的应用 　　【摘要】束缚量子系统的本征值精确可解模型稀少，如何利用数值方法得到系统的本征值是初等量子力学教学过程中不得不面对的问题之一。本文从线性谐振子模型出发，借助谐振子本征值态所建构的希尔伯特空间，利用Matlab提供的库函数求解任意束缚势量子系统的本征值和本征函数。数值方法的讲授将极大提高学生灵活应用量子力学知识解决实际问题的能力。 　　【关键词】谐振子 束缚态 薛定谔方程 Matlab 　　【中图分类号】O413 【文献标识码】A 【文章编号】2095-3089（2018）03-0048-02 　　引言 　　薛定谔方程束缚态离散谱的计算是初等量子力学所研究的基本问题之一。对于定态薛定谔方程本征值求解问题，由于可解的模型稀少，数值方法成为必不可少的解决方案。数值计算在量子力学教学方面如何取得突破是值得探讨的问题之一，其能有效帮助学生形成更加清晰的物理图像并能大幅提高学生初步掌握并解决实际问题的能力。本文不关注上述数值方法，而侧重于如何基于线性谐振子模型本征波函数所建构的正交完?浠?求解系统的本征值及其相应的本征函数。正交完备基所组成的希尔伯特相关性质在初等量子力学教学过程均会进行详细讲解，但是相关概念依旧相当抽象，学生难以掌握。本文所介绍的数值方法有助于学生理解相关概念。 　　Matlab是mathworks公司发布的集数值计算、矩阵计算、数据可视化以及系统建模仿真等功能的计算语言[1]。本文利用Matlab强大的矩阵运算和绘图功能实现定态薛定谔方程离散谱和相应本征函数的数值求解。文中数值计算所使用的主要Matlab代码已在正文中给出。 　　一、线性谐振子模型的本征值与本征函数 　　线性谐振子模型无论在经典力学还是在量子力学均占有重要地位，其应用领域非常广泛。一维线性谐振子所满足的定态薛定谔方程可写为：，其中和E分别为普朗克常数、系统哈密顿量、振子的有效质量、特征频率和能量。该系统的本征值和相应的本征波函数为[2] 　　（1） 　　和 　　（2） 　　其中n=0，1，2…，和。和分别为谐振子能级升降算符（或产生湮灭算符）。|n为谐振子本征能量标记为n的状态（即Fock态）。需要指出的是态|n具有正交归一性，即（克罗内克函数）。可证明的是坐标表象的本征波函数可构成完备的希尔伯特空间。上式中，为厄米多项式，其可显式写为： 　　（3） 　　其中，函数表示向下取整，Matlab标准库函数floor可实现该运算。由方程（2）和（3）可得到坐标表象下的本征函数，其可由以下自定义函数Harmonic实现： 　　function psi_n=Harmonic（alpha，x，n） 　　xt=sqrt（alpha）*x； 　　N_n=（alpha^2/pi）^0.25/sqrt（2.0^n*factorial（n））； 　　h = zeros（1，n+1）；n_fact=factorial（n）； 　　for m=0：floor（n/2） 　　h（2*m+1）=n_fact*（-1.0）^m/（factorial（m）... 　　*factorial（n-2.0*m））*2^（n-2.0*m）； 　　end 　　Herfun=polyval（h，xt）； 　　psi_n=N_n*Herfun.*exp（-1.0*xt.^2/2.0）； 　　需要指出的是，厄米多项式采用方程（3）计算，而不是Matlab提供的库函数hermiteH。该库函数运行十分缓慢，其原因在于该库函数根据厄米多项式递推关系生成高阶形式。 　　二、计算方法 　　接下来我们介绍如何利用线性谐振子本征态函数所建构的希尔伯特正交空间表示坐标和动量算符。由量子力学表象理论可知，坐标和动量算符具有矩阵形式，其矩阵元对应于直接内积运算，可表示为和（n和m=0，1，2…）。基于类似的矩阵形式，人们可以进一步的计算任意势函数的矩阵元形式，其形如。利用式子， 　　，和，上述矩阵元的值可表示为 　　（4） 　　和 　　（5） 　　相应的矩阵形式可写为： 　　（6） 　　和 　　（7） 　　需要指出的是实际的数值计算过程中，展开的基矢不可能取无穷多个，必须采取截断。此类截断将影响计算精度。如果指标n且其最大值为N（截断），那么上述坐标与动量算符将被截断成为的矩阵。在接下来的计算过程中，我们采用自然单位（）。 　　哈密顿量的矩阵表示都可以由方程（6）和（7）计算得到，即将坐标和动量算符直接由其矩阵形式替换，如。实际数值计算过程中，哈密顿量的矩阵表示也将被截断，本征值及其本征函数计算将被引入误差，且计算结果的收敛性与截断N大小有关[3]。如能得到哈密顿量矩阵的本征值，则该本征值对应该量子系统