一维量子谐振子的概率分布

一维量子谐振子的概率分布

摘要：线性谐振子问题作为一种普遍的模型，所以在经典力学中和量子力学中都受到很大关注。并且谐振子包括很多类型，我们就先研究量子谐振子的问题。量子谐振子是很多复杂物理模型的基础，量子谐振子在前几个量子态时，概率密度与经典情况相差较多，随着量子数的增加，随之相似性也会增加。可以通过使用数学软件将量子谐振子的概率分布绘制成图像，从而得出一维量子谐振子的概率分布。

关键词：经典谐振子 一维量子谐振子 波函数 量子谐振子概率分布

1.引言：

谐振子的振动是一种很常见的物理模型，它在很多方面得到应用。谐振子大体可分为经典力学和量子力学两部分，谐振在运动学就是简谐振动，这样的振动是物体在某一位置附近往复偏离该振动中心位置，在这样的振动方式下，物体所受到的力的大小总是与它偏离平衡位置的大小成正比关系，并且物体总是受到指向平衡位置的力。谐振子具有周期运动的物理特征，一些复杂的物理基础可以运用谐振子运动来解决。

通过对经典谐振子的研究，得到经典谐振子的函数关系式。再利用量子力学中的不确定关系得到量子谐振子的能量最低点，即平衡位置，最后得到谐振子的波函数，从而得到了谐振子的概率。随着量子数的增加，利用软件Mathematica绘制一维量子谐振子的概率分布。再和经典的线性谐振子来作比较，得到经典谐振子的关系。 2.经典一维谐振子：

首先让我们谐振子在物理中是非常常见的模型，我们很早就已经接触过 ，并且有了一定的了解。下面来讨论一维弹性力的一维简谐振子。例如：质量为m的物体放在光滑的桌面上，在其水平的方向上受到一个弹簧作用，在某一位置处质点所受力的大小为零，则把这一点叫做平衡位置。弹簧的劲度系数为k，物体m在弹簧弹性力的作用下沿弹簧方向运动，作用于质点的力和质点距离平衡位置的位移成正比，这样受力的质点就是一个典型的一维简谐振子。大家都知道，质量为m的质点在做简谐振动的过程中用x来表示质点便偏移平衡位置的距离，也就是质点的位置，也是弹簧的伸长或压缩的量。当x很小时，质点受力为F，则力F和x之间的线性关系为F??kx，并且可知弹簧的弹性力是线性回复力，弹簧振子

1

作简谐振动，再根据牛顿第二定律：

d2xm??kx， dt所以得运动微分方程为：

x''??在此中??2kx???2x， mk（决定于弹簧的劲度系数和滑块的质量），由此可以得到一个比较普遍md2x2??x?0的形式，也可以将它的解为的定义。比如质点运动的动力学方程为2dtx?Acos(?t??)。

相

应

的

，

质

点

的

动

量

的

表

达

式

为

p?mx?i?pxi??mA?sin(?t??)i??Bsin(?t??)i，这里出现的B?Akm 。这个

质点的动量的方向为x方向的分量。

还有一种典型的谐振子，那就是单摆，在研究单摆模型中，要用长度不可变的轻线悬挂一个小球，我们把小球看作为质点，质点所受的的合力为质点的重力和悬挂线拉力的合力，使得质点在竖直的平面内沿圆弧摆动，并且要求摆动相对于悬线竖直位置的夹角很小，把这个夹角记作? 。现在分析质点在沿运动方向所受的力为弧线的切线方向，记作F 。质点的

质量为m ，切向力F的大小为mgsin?，且当

?=0这个位置时，

sin?????23!??55!?··· ，所以当?角很小时，就可以略去级数展开式中的高次项，即

，从公式中就可以看出来，切向力和角位

sin??? ，这样切向力就可写成F??mg?

移反号，使得质点总要返回平衡位置，已知F力是线性回复力，所以会做简谐振动，可以得出单摆的动力学方程，假设线长为l ，所以：

gd2(l?)d2?g2??m??mg???? ， ，令 22lldtdt

即，可得到：

2

d2?2????0 2dt这样得到的结果和弹簧振子得到的结果是一样的，所以单摆得到的也是一个简谐振动。 3.量子谐振子

比如在一维的系统内粒子的势能

1m?2x2，其中?是常量，这种形式的称为线性谐2?Uk?0，这样就可以写成U=U0?(x?a)2，?x2振子。例如，两原子势能与x的关系，其中在两原子间距中有一个稳定平衡点，把这一点记作a，在x=a处，势能U具有极小值，即

式中的k和U0是常数，这就是一维线性谐振子的势能，通常情况下，一个体系平衡位置附近的运动都可以用线性谐振子来表示。谐振子的势能为1m?2x2 ，且坐标和时间可表示

2为x?asin(?t??)，我们把a作为振幅，?是初相。 我们选择适合的坐标系，领粒子势能为

1m?2x2 ，为了更加方便，我们引入了无量纲的变量?来代替x，所以该体系的 薛2定谔方程为：

?d2?m?22?(E?x)??0

2mdx22其中的关系为??

m?x??x ?，其中??m? ?。 令??2E ??，以

2乘方???d2?m?22?(E?x)??0， 由 ??程

2mdx22定谔方程变化为：

m?x??x和 ???m?式， 薛?d2?2?(???)??0 2d?2我们当?在????时的渐近时，也就是?非常大时，可得?和?相比可略去， 所以在

????

时，方程可表示为：

d2?2?(???)??0 2d? 3

d?可以写成为??2? ，该解为 ??e2d?的渐近解。

2??22d2? ，因此这就是方程?(???2)??02d?

波函数标准的条件是，当???? 时，?应该为有限的，因此取指数的负号，即

??22??ed2? 。所以就可以把?写成以下形式， 从而求得方程?(???2)??0的解为：2d???22?(?)?eH(?) ，从上式中求得的函数H(?)在自变量为有限时应有限的，而且当

????时，这样就的让H(?)必须保证?(?)为有限。将

??22?(?)?e 代入方程

H(?)

?d2?m?22?(E?x)??0 22mdx2先求出?(?)?e??22H(?)式的二级微商：

d?dH?2?(??H?)e d?d?22?2

?22

d?dHdH2?(?H?2???H?)e22D?d?d??

d2?将上式代入?(???2)??0式中，可得到H(?)所能满足的方程：2d?dHd2H(??1)H?2???0 2d?d?把H展成关于?的幂级数，而且级数只含有限的项的条件是?为奇数：??2n?1 ，且n=0,1,2,3··· 将式代入2公式??2E1后，就可得到线性谐振子的能级为E???(n?) ，??2n=0,1,2,3··· 因此，两个相邻能级的间隔均为?? ，线性谐振子的基态能量为

4

E0?1?? （n=0），这个被称之为零点能。这个在量子力学中是特有的，在旧量子论2中没有。 4．能量的最小值

根据不确定关系（测不准关系）可以的到动量和坐标的关系， 因为

??xpx?px?i? ，k??

于是得到：

?2(?x)(?px)?

422这个就是坐标和动量的关系，(?x)和(?px)2不能同时为零，其中坐标x的均方偏差越小，那么它的共轭的动量的均方偏差就越大。对于线性谐振子的零点能，我们可以利用不确定关系：

2?2(?x)(?px)?

422其中振子的平均能量是：

p21E??m?2x2 ，

2m2得到坐标的期望值是：

??x?Nn2???e??2x2Hn2(?x)xdx ，

该式中积分号下的函数是关于x的奇函数， 动量的期望值：

p?分步积分后：

?2eNn???i????2x2Hn2d?(?x)[edx?2x22HN(?X)]dx

5

p??均方差公式：

?2eNn???i????2x2d?Hn2(?x)[edx?2x22HN(?X)]dx??p

(?F)2?(F?F)2?F2?2FF?F2?F2?F2

因为x?0，p?0 ，得到(?x)2?x2 ，(?px)2?p2 ，将其代入

p21E??m?2x2

2m2得线性谐振子的能量期望值：

(?p)21E??m?2(?x)2

2m2又由于不确定关系使(?x)和(?px)2不能同时为零，所以E的最小值也就不能为零，

2?2并且必须是有限的正值。为了求得E的最小值，使得(?x)(?px)?取等号，即可得

422到：

?2 (?px)?24(?x)2 将该式代入

(?p)21E??m?2(?x)2式，

2m2得到

?21122E??m?(?x) 28m(?x)2将此式对(?x)求导就得到了E的最小值，又由(?x)22?? ，得到了E的最小

2m? 6

值为

1m? 。从以上关系可以看出，线性谐振子的基态能量是由不确定关系所求得的2最小值。

5．线性谐振子波函数

对于??2n?1 ，且n=0,1,2,3··· 中不同的n或者不同的? ，在方程

dHd2H(??1)H?2???0中有着不同的解Hn(?) ，我们把Hn(?)称之为厄2d?d?米多项式，也可以用Hn(?)n?(?1)ne?2dn??2end? 来表示。其中的Hn(?)的最高

次幂是n ，并且它的系数是2 ，所以就可以利用

Hn(?)?(?1)e推出

n?2dn??2e nd?dHn?2nHn?1(?) d?Hn?1(?)?2?Hn(?)?2nHn?1(?)?0

Hn(?)?(2?)n?n(n?1)(2z)n?2n?2[]n!2?...?(?1)(2z)

n[]2n2n?nn(n?1)n??2[式中的]?? ，n为偶数时为 ，当n为奇数时为。

222?(n?1)??2从而得到了几个厄米多项式：

H0?1 H1?2? H2?4?2?2 H3?8?3?12? H4?16?4?48?2?12

7

H5?32?5?160?3?120? H6?64?6?480?4?720?2?120 H7?128?7?1344?5?780?3?384?

H8?256?8?3584?6?8280?4?10848?2?1680 H9?512?9?9216?7?38064?5?34176?3?9504?H10?1024?10?23040?8?140640?6?217392?4?214272?2?30240

由?(?)?e??22H(?)得到的关于能量En的波函数，即

?n(?)?Nne或?n(x)??22Hn(?)

?Nne??22x2Hn(?x)

式中的Nn是归一化因子，所以利用归一化条件:

??*??n(x)?n?(x)dx??n,n?

??（其中归一化因子Nn?()?2nn!12?12） 。由以上关系可以得到

?n(?x)?(?1)n?n(x)的关系。

从而得到谐振子的波函数如下

?0(x)???14e??2x2/2

?1(x)??2(x)?2???14?xe??2x2/222?14(4?2x2?2)e??2x2/2

8

?3(x)??4(x)??3?14(2?x?3)e22??2x2/2?24?14(4?4x4?12?2x2?3)e??2x2/2?5(x)??60?14(4?x?20?x?15?x)e4422??2x2/2?6(x)??125?14(8?6x6?60?4x4?90?2x2?15)e??22x/2?7(x)??8(x)??670??14(8?x?84?x?210?x?105)e664422??2x2/22470?14(16?8x8?224?6x6?840?4x4?840?2x2?105)e??22x/2?9(x)???7235?14(16?8x8?1512?4x4?2520?2x2?945)e??22x/2

?10(x)?7207?14(32?10x10?720?8x8?5040?6x6?12600?4x4?9450?2x2?945)e??22x/2

6.量子谐振子的概率分布：

所以从谐振子的波函数来获得线性谐振子的概率分布

即，概率密度为?(?)??(?) ，n=0,1,2,3··· ，其中???x 所以就可以得到谐振子的概率：

2W(0)?1?12?1e???

2W(1)?e??2??2

W(2)?3?13?e???(2?2?1)2

2W(3)?e??2??2(2?2?3)2

9

W(4)?124?160?1e???(4?4?12?2?3)2

2W(5)?e????2(4?4?20?2?15)2

2W(6)?720?1e???(8?6?60?4?90?2?15)2

2W(7)?2520?1e??2??2(8?6?84?4?210?2?105)2

W(8)?40320?1181440?1e???(16?8?224?6?840?4?840?2?105)2

2W(9)?e??2??2(16?8?288?6?1512?4?2520?2?945)2

W(10)?3628800?e???(32?10?720?8?5040?6?12600?4?9450?2?945)2

2利用Mathematica画出谐振子的概率分布，图分别如下： 图一：

2.01.51.00.5?1.0?0.50.51.0 基态（n=0）是的概率密度（6——1）

图二：

10

0.350.300.250.200.150.100.05?2?112 n=1时波函数的概率密度 （6——2）

图三：

15010050?4?224 n=4时波函数的概率密度 （6——3）

11

图四：

15000010000050000?6?4?2246 n=6时波函数的概率密度 （6——4）

图五：

5?1084?1083?1082?1081?108?55 n=8时波函数的概率密度 （6——5）

图六：

12

3?10122?10121?1012?6?4?2246 n=10时波函数的概率密度 （6——6）

总结：

从上图可得到，量子谐振子的概率分布随着量子态的增加而呈现规律，当量子态很小的时候，这时的概率密度与经典的谐振子几乎没有什么相似处，但是随着n的增大，它们之间的相似性会随着增加，当达到n=10的时候，这时的量子谐振子和经典谐振子的情况几乎一样，其中的一点差别就是量子谐振子的?(?)会迅速的震荡而已。

2参考文献：

?1? 曾谨言 .【量子力学导论】第二版 ?M? . 北京大学出版社，2004

?2? 周世勋 .【量子力学教程】第二版 ?M? . 高等教育出版社，2009 ?3? 漆安慎 杜婵英 【.力学普通物理教程】第二版 ?M? . 高等教育出版社，2005 ?4? 陈世民 .【理论力学简明教程】第二版 ?M? . 高等教育出版社，2008

】 第四版 ?M? . 科学出版社，2007 ?5? 曾谨言 .【量子力学（卷?）

?6? 王正行 .【量子力学原理】 ?M?. 北京大学出版社，2003

?7? 刘明 .【经典谐振子与量子谐振子】 ?J? 湖北教育学院 物理与电子工程系，

武汉 430205

?8? 徐秀玮，孔祥伟，田丽杰，柳盛典，朱友亮，赵旭光 【量子谐振子与经典谐振子的比较】 ?J?

?9?Baym.G .Lecteres on Quantum Mechanics Benjam in New York [M] ,1987.

?10?May R M . Nolinear Problems in ecology and resource management [M] . North

Holland Publ Co,1983

13

[11]宋鹤山 . 量子力学 . [M] 大连理工大学出版社，2004 [12]姚玉洁 . 量子力学 . [M] 吉林大学出版社，1988 [13]刘全慧 . 等全波包与一维谐振子[J] . 大学物理，2002

[14]王可，陶才德 . 一维线性谐振子量子束缚态能级的几种计算方法 [J] . 绵阳师范学院学报，2005

[15]文焕邦，刘敬乾 . 量子力学 [M] . 四川科学技术出版社，1986

One-dimensional quantum harmonic oscillator probability distribution Abstract: Linear Resonance problem as a general model，Therefore, in classical mechanics and quantum mechanics are greatly concerned. And resonator includes many types, We first study the quantum harmonic oscillator problem. Quantum harmonic oscillator is the basis of many complex physical model. When the first few quantum harmonic oscillator quantum state, Classical probability density difference between the more, With the increase of quantum numbers, Along with similarity will increase. By using mathematical software quantum harmonic oscillator probability distributions graphed, To arrive at the one-dimensional quantum harmonic oscillator probability distribution.

14

致 谢

在写论文的过程中，由于自己学习的不够扎实，遇到了很多困难，经过老师不厌其烦的讲解和指导，最终才将其完成。在此过程中深深的感觉到作为一个合格的老师的不容易，也深深的被老师扎实的物理功底所折服。在向解老师学习知识的同时，也应该向他学习对于学术那种严谨和认真的态度，无疑对我以后是一种激励。还有就是曾经教导我的老师们，是老师们一点一点的努力和汗水让我们懂得越来越多，没有他们曾经的努力也就没有现在的我们。

我的论还有很多不足，还需要老师们的指导和纠正，我会虚心学习和认真听取老师的建议，感谢各位老师，您辛苦了！

15

本文来源：https://www.bwwdw.com/article/d2dp.html