经典力学与量子力学中的一维谐振子[摘要]一维谐振动是一种最简单的振动形式，许多复杂的运动都可分析为一维谐振动。

本文以一维谐振子为研究对象，首先讨论经典力学与量子力学中的一维谐振子的运动方程和能量特征，然后分析坐标表象以及粒子数表象下的一维谐振子，最后讨论经典力学与量子力学中的一维谐振子的区别与联系。

[关键词]谐振子经典力学量子力学运动方程能量分布1 前言所谓谐振，在运动学中就是简谐振动。

一个劲度系数为k的轻质弹簧的一端固定，另一端固结一个可以自由运动的质量为m的物体，就构成一个弹簧振子[1]。

该振子是在一个位置（即平衡位置）附近做往复运动。

在这种振动形式下，物体受力的大小总是和它偏离平衡位置的距离成正比，并且受力方向总是指向平衡位置。

这种情况即为一维谐振子。

一维谐振子在应用上有很大价值,因为经典力学告诉我们只要选择适当的坐标,任意粒子体系的微小振动都可以认为是一些相互独立的振子的运动的集合。

普朗克在他的辐射理论中将辐射物质的中心当作一些谐振子,从而得到和实验相符合的结果。

在分子光谱中,我们可以把分子的振动近似地当作谐振子的波函数。

另外在量子场论中电磁场的问题也能归结成谐振子的形式。

因此在量子力学中,谐振子问题的地位较经典物理中来得重要。

应用线性谐振子模型可以解决许多量子力学中的实际问题。

本文将以一维谐振子为研究对象，首先分别讨论经典力学与量子力学中一维谐振子的运动方程和能量特征，然后讨论坐标表象以及粒子数表象下的一维谐振子，最后分析经典力学与量子力学中的一维谐振子的区别与联系并简要讨论经典力学与量子力学的过渡问题。

从而帮助我们更加深入的理解一维谐振子的物理实质，充分认识微观粒子的波粒二象性。

2 经典力学中的一维谐振子在经典力学中基本方程以牛顿定律为基础，研究质点位移随时间变化的规律，反映质点特征的是运动方程和能量。

因此我们可以从运动方程和能量这两方面出发讨论一维谐振子的运动特征。

一个劲度系数为k 的轻质弹簧的一端固定，另一端固结一个可以自由运动的质量为m 的物体，就构成一个弹簧振子[1]，如图2.1。

当弹簧处于自然长度时，物体处于平衡位置，取作坐标原点，以O 表示。

沿弹簧长度方向（取作x 轴方向）拉动物体然后释放，则物体将在O 点两侧作往复运动。

图2.1 弹簧振子2.1 一维谐振子的运动方程图2.1中的物体可视为一个质点。

设x 代表质点相对于平衡位置的位移，则质点所受的力kx F -=，其中k 为劲度系数。

负号表示F 与位移方向相反，因而总是指向平衡位置。

由牛顿第二定律，谐振子的运动微分方程为：kx xm -= 即 02=+x x ω （2.1.1）这是一个二阶的常系数线性微分方程。

令mk =ω （2.1.2） ω即简谐运动的角频率，由振动系统本身的性质嗦决定。

将（2.1.2）式代入（2.1.1）式，则可求出（2.1.1）式的通解：iwt iwt Ne Me t x -+=)(（2.1.3） )sin(ϕω+=t A这就是谐振子的运动方程[2]。

其中M 和N 是任意常数，由质点的初位置和初速度确定。

A 是振幅，ϕ是初相位。

（2.1.3）式表明质点应作简谐振动[2]。

2.2 一维谐振子的能量在谐振子问题中，振子的总能量可以反映出振子的运动特征。

因此我们可以从谐振子的动能和势能出发，求解谐振子的总能量，进而帮助我们分析振子的运动特征。

由（2.1.3）式可知，振子的速度为：)cos(ϕωω+==t A dtdx v 振子的动能为：)(cos 21)(212122222ϕωω+===t C m dt dx m mv E k 由（2.1.2）式，有： )(cos 2122ϕω+=t kA E k （2.2.1） 由（2.2.1）式可知，振子的动能变化频率为ω2。

振子的势能（以平衡位置的势能为零）为：2021kx Fdx E xp =-=⎰ 即为：)(sin 2121222ϕω+==t kA kx E p （2.2.2） 由（2.2.2）式可知，振子的势能变化频率也为ω2。

因此，由（2.2.1）式和（2.2.3）式可得，振子的总能量为：221kA E E E p k =+= （2.2.3） 由（2.2.3）式可知：谐振子的总能量不随时间改变，即其机械能守恒[3]。

（2.2.3）式还说明：对于一定的振子（m 和k 给定，因而ω给定），总能量与振幅的平方成正比[3]。

振幅不仅给出了简谐运动的运动范围，而且还反映了振动系统总能量的大小，或者说反映了振动的强度。

3 量子力学中的一维谐振子在量子力学中，粒子状态用波函数表示，为了描述微观粒子状态随时间变化的规律，就需要找出波函数所满足的运动方程，即薛定谔方程。

因此下面将从谐振子的哈密顿算符出发，求解振子的定态薛定谔方程，进而分析量子力学中一维谐振子的运动特征。

3.1 用运动方程求解的一维谐振子我们可以从谐振子的势能函数出发，写出谐振子的哈密顿算符及薛定谔方程，并求谐振子的能量和定态波函数的解，进而讨论能量分布特点。

取谐振子的平衡位置0r 为坐标原点，并选原点为势能的零点，则有0)(0=r E p 。

仅考虑一维情况。

由于kz j y i x r ˆˆˆ+-= 在x 轴方向分振动的谐振子在x 处的势能可以表示为：221)(kx x E p =（3.1.1） 势能曲线是一条定点在原点的抛物线，如图3.1所示：图3.1 一维谐振子的势能一维谐振子的经典哈密顿函数为：22212kx m p H +=设振子的原子质量为μ，则振子的频率为：m k=ω振子的哈密顿算符可以写为：22222ˆ212x m dx d m H ω+-=相应的定态薛定谔方程)(ψψE H = 为：)()()212(22222x E x x m dx d m ψψω=+- （3.1.1）这是一个二阶线性微分方程。

如果振子的运动不受限制，x 的变化范围为+∞<<∞-x 。

当∞→x 时，（3.1.1）式的解一般为无穷大，表示振子在无穷远处的几率为无穷大。

这不符合物理要求。

但若振子的能量E 取下列特殊值[4]：ω )21(+=n E )2,1,0( =n （3.1.2） 其中 为普朗克常数，ω为经典力学中谐振子的频率。

则对每一个n 值，方程（3.1.1）都有一个在全区间+∞<<∞-x 中有界的解。

而且当∞→x 时，这个解趋于零。

这显然符合对谐振子问题的物理要求。

与（3.1.2）式的能量值相应的关于定态波函数的解为： )()(2x H Ae x n x αψα-= （3.1.3） 其中0>=hkm α，)(ξn H 是E 的一个n 次多项式，称为厄米多项式[4]。

其前四项为： 1)(0=ξHξξ2)(1=H24)(22-=ξξHξξξ128)(33-=H 由于)(x H n α是x 的n 次多项式，且0>α。

因此，当∞→x 时，（3.1.3）式趋于零。

由（3.1.2）式知，在量子力学中谐振子的能量是分立的，与振幅无关，只依赖于振子的固有特性---振子的本征频率ω。

（3.1.2）式还表明[5]，频率为ω的振子其能量的改变只能是能量单元ω 的整数倍。

这一点同普朗克的能量子假设是一致的。

但量子力学中振子的最低能级（基态能量）不再是零而是ω 21，称为零点能[5]。

它充分体现了粒子具有波粒二象性。

3.2 坐标表象中的一维谐振子粒子系统的状态用以空间坐标为自变量，以时间为参量的波函数),(t r 来描述，这种表示形式称为坐标表象[6]。

下面从谐振子的哈密顿算符出发，求解谐振子的能量本征值和定态波函数，并对谐振子在量子力学与经典力学中的几率密度进行比较，给出量子谐振子向经典谐振子过渡的条件。

一维谐振子的哈密顿算符为：22222212x m dx d m H ω+-= 坐标表象中，振子的定态薛定谔方程为：)()()212(22222x E x x m dx d m ψψω=+- 引入没有量纲的变量ξ代替x ，它们的关系是x x m αωξ≡≡ ωαm = ωλ E 2=（3.2.1） 以ω2乘以式（3.1.1），利用式（3.1.2）和式（3.1.3），薛定谔方程可改写为 0)(222=-+ψξλψξd d （3.2.2） 这是一个变系数的二阶常微分方程。

当ξ很大时，λ与2ξ相比可以略去，因而在±∞→ξ时，（3.2.2）式可以写为：ψξξψ222=d d 它的解是ψ～22ξ±e 。

因为波函数的标准条件要求当±∞→ξ时，ψ应为有限，所以对波函数只取指数上的负号：ψ～22ξ-e 。

根据上面的讨论，我们把ψ写成如下形式来求（3.2.2）式的解：)()(22ξξψξH e -= （3.2.3）（3.2.3）式代入（3.2.2）式可得)(ψH 满足方程：0)1(222=-+-H d dH d H d λξξξ（3.2.3） 用级数解法，把H 展开成ξ的幂级数，来求着方程的解。

这个级数必须只含有限项，才能在±∞→ξ时使)(ξψ为有限；而级数只含有限项的条件是λ为奇数：12+=n λ， 2,1,0=n代入（3.2.1）式即可得谐振子的能级为：)21(+=n E n ω ， 2,1,0=n （3.2.4） 可见,谐振子的能量只能取分立值，两相邻能级间隔均为ω ，即：ω =-+n n E E 1在坐标表象中可以明显看出：是描述粒子波动性的波函数)(x ψ受到势能场的约束使能谱分立。

与（3.2.4）式定态能量对应的定态波函数：)()(222x H eN x n x n n αψα-= 式中2121)!2(n N n n πα=是归一化常数，它由归一化条件1)()(*=⎰∞∞-dx x x n n ψψ确定。

图3.2中画出了3,2,1,0=n 的几率密度2n ψ（图中实线），图中虚线是经典谐振子的平均位置密度。

从图3.2可见，经典谐振子不能进入A x >的区域。

而量子谐振子能进入这种区域，但进入以后指数衰减。

可见，量子振子和经典振子完全不同，但当n 增大时，nE E ∆减小，量子振子向经典振子过渡。

图3.2 一维谐振子的位置几率密度分布 3.3 粒子数表象中的谐振子以粒子数算符的本征矢|n 〉为基矢的表象称为占有数表象[7]，又叫粒子数表象。

我们可以通过引入升降算符求解谐振子，求出谐振子的能量本征值以及坐标算符xˆ和动量算符p ˆ的矩阵元。

一维谐振子的经典哈密顿函数为：22212kx p H +=μ （3.3.1） 在量子力学里，谐振子的哈密顿算符具有同一形式：222ˆ21ˆ21x p H μωμ+-= （3.3.2） 将经典泊松括号换成量子泊松括号：[]1ˆ,ˆ-=PB p x → []1ˆ,ˆ-=p x i由xˆ与p ˆ的对易关系： [] i x p p x p x =-=ˆˆˆˆˆ,ˆ （3.3.3）定义两个非厄米算符aˆ和+a ˆ： )ˆ1ˆ(2ˆx pi x a μωμω+=)ˆ1ˆ(2ˆx pi x a μωμω-=+ （3.3.4） 这两个非厄米算符满足如下基本对易关系：][1ˆ,ˆ=+a a （3.3.5） 则（3.3.4）式的逆变换关系为： )ˆˆ(2ˆa a x +=+μω（3.3.6） )ˆˆ(2ˆa ap -=+μω 利用（3.3.6）式，代入（3.3.2）式，并考虑对易关系（3.3.5），哈密顿算符又可表示为：)21ˆˆ(ˆ+=+a a H ω （3.3.7） 由于Hˆ与算符a a ˆˆ+仅仅相差一个常数矩阵，所以只需求解a a ˆˆ+得本征值问题。