一维线性谐振子薛定谔方程的简易揭发以及其波函数和几率的作图.docVIP

一维线性谐振子薛定谔方程的简易解法以及其波函数和几率密度的作图 林乐鑫 (华南师范大学物理与电信工程学院，广东广州，510006) 摘要：量子谐振子在量子力学理论中占有重要地位，物质的许多物理与化学性质都可以用线性谐振子模型解释。就教学而言，一维线性谐振子是教学中的难点，一般教材仅限于对Schrodinger方程的直接求解，本文将阐述一种简易的方法来求解一维线性谐振子薛定谔方程，并通过科学计算与模拟平台同时画出可调能级的波函数和几率密度曲线。在多媒体教学和学生的认识理解中起到了重要的作用。 关键词：谐振子，薛定谔方程，波函数，几率密度，科学计算机与模拟平台 无论是在经典力学还是在量子力学的理论中，谐振子都是重要的例子。它不仅在原理上十分典型和重要，而且应用非常广，物质的许多物理与化学性质都可以用线性谐振子模型解释。许多体系可近似地看作是线性谐振子。一般说来，任何一个体系在稳定平衡点附近振动都可以近似的用线性谐振子表示。 因此怎样在科学理论研究中处理相关问题时或者教学时迅速找到简单易行的方法来解这个薛定谔方程，显得尤为重要。在多媒体教学中，我运用科学计算与模拟平台，能够同时画出可调节能级的波函数和几率密度曲线，可以非常直观的进行图像的比较和讨论，不仅方便了多媒体教学，而且加深了学生对一维线性谐振子的理解。 1 解一维谐振子薛定谔方程 一维线性谐振子的薛定谔方程为 （1） 为了方便解题，我们引进新变量 （2） 可得 代人（1）并记 ，得 （3） 易知方程（3）的解为 （4） 根据波函数的边值条件，考虑到 在 时收敛,并且与相比 可以略去, 化简后得到方程： （5） 易知方程（5）的解 （6） 我们可以猜想方程解的形式为 （7） 其中 是一个未知函数，待求函数在 时也应收敛，只有这样才能满足波函数的标准条件。 易得到 满足方程 （8） 显然方程（8）是厄米特方程，记为。 前两个厄米特多项式为 再根据厄米特多项式递推关系可以求解： （9） 对应能量的波函数是 （10） 式中 是归一化常数，它由归一化条件 易得： （11） 综上所述，该波函数为： （12） 几率密度为： 2 绘图 用一般的计算机语言编程，是可以绘制出波函数和几率密度曲线，不过很麻烦，需要进行进一步推导，而且一次只能显示一个能级的图。 而科学计算与模拟平台是利用VC++和OPENGL的结合，二次开发创作的适合高中以上师生学习使用的软件平台，是一种全新的学习使用高级语言去定量计算和模拟自然科学过程的教学工具，它特别适用于培养学生研究式学习和教师进行创新教育的教学模式。现在我们利用这个平台来作图。 （1）先编辑一个子程序计算积分 double Jiecheng(int n) { int m; if(n==0) { m=1; return m; } if(n==1) { m=1; return m; } m=n*Jiecheng(n-1); return m; } （2）再编辑一个子程序计算厄米特多项式 double Hen(int n,double v) { double h; if(n==0) { h=1; return h; } if(n==1) { h=2*v; return h; } if(n=2) { h=2*v*Hen(n-1,v)-2*(n-1)*Hen(n-2,v); return h; } } （3）编辑主程序 void demoApp::RenderScene(int sceneIndex) { title.Show(-80,55, 0, false); title1.Show(-80,0, 0, false); title2.Show(-80,-50, 0, false); Point3f p[120],g[120]; Color4f color={0,0,0,1},color1={0.2,0.8,0.9,0.8}; int i,n; double u; n=P_radius; for(i=0;i120;i++) { g[i].x=-60+ i; g[i].z=0; g[i].y=14*sq