一维谐振子的数值方法研究
第３０卷第６期 佳木斯大学学报（自然科学版） ２０１２ 年１１月 ＪｏｕｒｎａｌｏｆＪｉａｍｕｓｉＵｎｉｖｅｒｓｉｔｙ（ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ）
文章编号：１００８－１４０２（２０１２）０６－０９５８－０３
Ｖｏｌ．３０Ｎｏ．６
Ｎｏｖ． ２０１２
一维谐振子的数值方法研究
陈海军
（陇东学院物理系，甘肃庆阳７４５０００）
①
摘 要： 利用数值计算方法研究了量子力学中的一维谐振子问题．简述了解析方法求解一维谐振子问题的理论过程，给出体系的能级及对应的定态波函数．利用数值计算方法求解定态薛定谔方程，得出了前三个能级和对应的波函数，结果与理论计算结果完全一致．最后利用有限差分法求解含时薛定谔方程，验证定态解的稳定性，表明定态解在整个演化过程中是稳定的．关键词： 谐振子；数值；定态；稳定性中图分类号： ＴＮ９７２．４４ 文献标识码： Ａ
¨式的定态Ｓｃｈｒｏｄｉｎｇｅｒ方程
２
２＋（λ－ζ）Ψ（ζ）＝０（２）２
２显然，当ζ→±∞时，Ψ（ζ）～ｅｘｐ（），
因此，令
２方程的解为
０ 引 言
谐振子问题是量子力学和经典力学中的基本问题．任何体系在平衡位置附近的振动，例如，分子的振动，晶格的振动，原子核表面振动以及辐射场的振动等等，在适当的坐标变换之后，都可以分解为若干彼此独立的谐振动．谐振动还可以作为任何复杂振动的初步近似．在量子力学中，谐振子问题的精确求解过程，不仅提供了一整套量子力学处理具体问题的方法，而且为处理其它量子力学问题提供了物理基础．
鉴于大多数量子力学教科书都提供了一维谐振子问题的解析求解过程，本文简述一维谐振子问题的理论求解过程，然后利用数值计算的方法求解其定态解，并研究其随时间的演化规律．这种数值求解过程是数值方法研究其它量子问题的基础．
图１ 解析解波函数及波函数的模方
１ 谐振子问题理论求解
根据量子力学教科书简述一维谐振子问题的理论求解过程．
¨描述体系的定态Ｓｃｈｒｏｄｉｎｇｅｒ方程是
２２
２２ωｘΨ＝ＥΨ（１）２２ｍ２
［１］
２
Ψ（ζ）＝Ｈ（ζ）ｅｘｐ（）
２
代入方程（２）得到
－＋（λ－１）Ｈ（ζ）＝０２
Ｈ（ζ）＝∑ａｖζ
ｖ＝０∞
ｖ
２
（３）
（４）
级数展开求解，令
选取无量纲单位，长度单位α＝能量单
位令ζ＝ｘ／α，λ＝Ｅ／（礹ω／２），得到无量纲形２
（５）
代入方程（４）得到递推公式
ａｖ＋２＝
ａｖ
（ｖ＋１）（ｖ＋２）
（６）
①
收稿日期：２０１２－１０－０１
作者简介：陈海军（１９８２－），男，甘肃镇原人，陇东学院物理系讲师．主要研究方向：量子力学教学法和冷原子体系的稳定性．
第６期
２
陈海军：一维谐振子的数值方法研究９５９
方程（６）表明级数Ｈ（ζ）和ｅｘｐ（ζ）具有相同的收敛性，是发散的，因此必须取有限项使得Ｈ（ζ）变为厄米多项式Ｈｎ（ζ），则有
ａｖ＝ｎ≠０，ａｖ＋（７）２＝ａｎ＋２＝０
得到λ＝２ｎ＋１，则能量
Ｅｎ＝（ｎω（８）
２
因此和Ｅｎ相应的归一化波函数是
Ψ）＝ＮｎＨｎ（ζ）ｅｘｐ（Ｎｎ＝ｎ（ζ２
２２
ｎ
ｎ！
（９）
其中Ｎｎ为归一化系数．
图２ 数值计算得到前三个能级所对应的波函数
图３ 定态随时间传播过程
表１ 前三个状态的能量及波函数
ｎ
λ
ＥｎΨｎ（ζ）
ｎ＝０１
礹ωα２４
（ζ２２ｎ＝１３
３礹ω２２ｅｘｐ（１
４
２α
ｎ＝２５
２
－１）ｅｘｐ（－
２
１
２４
２
表１给出前三个状态的能量及其波函数．图１给出前三个解析解对应的波函数模方，即粒子出现概率密度的空间分布．可以看出处于基态和第二激
９６０佳木斯大学学报（自然科学版）
ｉ（Ψｎ＋１
ｎ
２０１２年
发态（ｎ＝０，ｎ＝２）的粒子在势阱中心有一定概率
出现，其中，处于基态的粒子在势阱中心出现的概率最大，而处于第一激发态的粒子在势阱中心出现的概率为０．
－Ψ）ｎ＋１ｎ＋１ｎ＋１＝２ｊ＋１－２Ψ＋Ψｊｊ－１）
τ２ｈｎｎｎｎｎ
＋（Ψｊ＋１－２Ψｊ＋Ψｊ－１）］＋ＶｊΨ（１１）ｊ
２
２ 定态问题数值计算
量子力学中的大多数问题不可能求出解析解，
因此利用数值计算方法研究谐振子问题具有重要意义．该方法不仅可以求解谐振子问题本身，而且其中Ｖ＝ｘ是势能函数．
２
图３表示ｎ＝０，１，２所对应定态｜Ψ｜随时间的传播过程，即粒子出现概率空间分布随时间变化过程．图（ａ）表示基态概率分布随时间的传播过程，其实质是高斯波包在抛物型势阱中的传播过程，如果没有势阱的作用，波包由于色散作用，随时可以扩展到其它的量子问题上，只要把其中的势能函数替换为想要求解体系的势能函数即可．束缚定态问题的数值求解，在数值计算方法中归结为偏微分方程的边值问题求解．这类问题求解多采用打靶法，即取微分方程中某一个参数作为数值计算的试探参数，利用计算机快速计算的能力，在参数取值的物理区间内，按一定的步长改变参数值，直到得出的数值结果满足物理问题要求的边界条件即可．
本文利用Ｍａｔｈｅｍａｔｉｃａ软件中的ＮＤＳｏｌｖｅ命令
求解方程（２）的边值问题［２］
．
自变量的变化范围是－∞＜ζ＜∞，在数值计算时取－４≤ζ≤４．由于求解的是束缚定态问题，因此波函数Ψ（ζ）在边界上的取值应为０，数值计算时，采用｜Ψ（±４）｜＜１０－３
．另外，由于势能Ｖ（ζ）具有对称性，所以体系具有明确的宇称，在计算过程中偶宇称解和奇宇称解的求解分别进行．在计算过程中，取λ为试探参数，从０开始循环，步长为０．００１，直到定态结果出现为止．
计算结果如图２所示，图２（ａ）和（ｃ）表示偶宇称波函数（ｎ＝０，２），图２（ｂ）表示奇宇称波函数（ｎ＝１）．数值计算结果给出λ＝１，３，５，对应的能量是Ｅ＝礹ω／２，３礹ω／２，５礹ω／２，和理论结果一致．
３ 定态稳定性研究
Ｎｉｃｈｏｌｓｏｎ为了验证定态解的稳定性时间步长τ差分法求解含时和空间步长ｈ分别为Ｓｃｈｒ，我们利用Ｃｒａｎｋ－
¨ｏｄｉｎｇｅｒ０．００１和方程［３，４］
，
０．１，ｘ方向的网格数为２５０．把表１中的定态解作为求解含
时Ｓｃｈｒ¨化的规律ｏｄｉｎｇｅｒ．
方程的初始条件，研究其随时间变
含时Ｓｃｈｒ¨ｏｉｄｉｎｇｅｒ方程的无量纲形式为
２
矪ｔ＝２矪ｘ
２＋ｘΨ（１０）其中时间无量纲化单位是２／ω．差分格式是
间传播会扩散，而由于势阱的束缚作用，高斯波包在传播过程中保持初始时刻的波形，是稳定的传播过程．第一激发态的传播过程也是稳定的（图（ｂ））．图（ｃ）表示第二激发态随时间的传播过程，和前两个状态的传播过程不同，在传播时，两边的概率会向势阱的中心移动，但是并没有坍缩．另外，整个传播过程是呼吸振荡的．
４ 结 论
本文利用数值计算方法研究了量子力学中一维谐振子问题．包括定态和含时过程的求解，定态计算结果给出了能量和定态所对应的波函数，和理论分析结果一致．含时过程计算表明处于定态中的粒子出现概率空间分布是不随时间变化的，这是定态的基本性质之一．整个求解过程中涉及到定态问题和含时问题的数值计算，这两类问题的数值计算是研究量子问题的基础，尤其体系稳定性问题，由于能进行理论分析的体系很少，所以只能依靠数值求解含时薛定谔方程来研究．因此本文介绍的方法可以推广到其它量子体系的研究，甚至包括非线性薛定谔方程所描述的体系．参考文献：
［１］ 苏汝铿．量子力学（第二版）［Ｍ］．北京：高等教育出版社，
２００２．
［２］ 董键．Ｍａｔｈｅｍａｔｉｃａ与大学物理计算［Ｍ］．北京：清华大学出版
社，２０１０．
［３］ 彭芳麟．计算物理基础［Ｍ］．北京：高等教育出版社，２０１０．［４］ ＳａｄｈａｎＫＡｄｈｉｋａｒｉ，ＰａｕｌｓａｍｙＭｕｒｕｇａｎａｎｄａｍ．Ｂｏｓｅ－Ｅｉｎｓｔｅｉｎ
Ｃｏｎｄｅｎｓａｔｉｏｎ－ＰｉｔａｅｖｓｋｉｉＤｙｎａｍｉｃｓＥｑｕａｔｉｏｎ［ｆｒｏｍＪ］．ｔｈｅＪ．ＰｈｙｓＮｕｍｅｒｉｃａｌ．Ｂ：ＡｔＳｏｌｕｔｉｏｎ．Ｍｏｌ．Ｏｐｔｏｆｔｈｅ．ＰｈｙｓＧｒｏｓｓ．，
２００２，３５：２８３１－２８４３．