谐振子-相图_153601860
- 格式:pdf
- 大小:391.51 KB
- 文档页数:9
各向异性谐振子的能级简并刘永宏 指导教师:焦志莲 (太原师范学院物理系,太原030031)【摘 要】 给出了二维、三维各向异性谐振子的能级及波函数,并讨论各种情况下二维、三维各向异性谐振子的能级简并问题。
【关键词】各向异性谐振子,能级,波函数,能级简并。
0. 引 言各向同性谐振子的能级简并问题,很多量子力学教材都进行了讨论,譬如:曾谨言写的《量子力学导论》就对各向同性谐振子作了详细而深刻的分析。
但是对于各向异性谐振子的问题,则很少有教材中进行专门的讨论。
各向异性谐振子有其独特的能级简并和对称性,且在一定的近似条件下,可转变为各向同性谐振子来处理。
所以对于各向异性谐振子的能级简并研究,既能进一步加深对各向同性谐振子的理解和应用,同时又能为学习和探究更深层次的各向异性谐振子奠定基础。
本文先给出二维,三维各向异性谐振子的能级及波函数,然后讨论相应各向异性谐振子的能级简并度问题。
1. 各向异性谐振子的能级及波函数1.1 二维各向异性谐振子的能量及波函数当各向异性谐振子为二维情况时,体系哈密顿量在oxy 坐标系中可以表示为222222(,)112222y x x y x y P P H x y μωμωμμ=+++ (1) 令 22222211,2222yx x x y y P P H x H y μωμωμμ=+=+ (2)求解哈密顿本征值方程,可以得体系能量及波函数的表示为,11()()22x y n n x x y y E n n ωω=+++ (3),(,)()()x y x y n n n x n y x y x y ααψ=ψψ (4)其中,各维波函数为221()exp()();0,1,2,2x n x x x x x x x x N x H x n ααααψ=-== (5)221()exp()();0,1,2,2y n y y y y y y y y N y H y n ααααψ=-== (6)1.2 三维各向异性谐振子的能量及波函数在三维空间o xyz -中,三维谐振子的哈密顿量为222222222(,,)111222222y x z x y z x y z P P P H x y z μωμωμωμμμ=+++++ (7) 令 222222222111,,222222y x z x x y y z z P P P H x H y H z μωμωμωμμμ=+=+=+ (8)由三维谐振子体系哈密顿量的本征值方程,可以求出的体系哈密顿量的本征值及相应的本征值函数为,,111()()()222x y z n n n x x y y z z E n n n ωωω=+++++ (9),,(,,)()()()x y zx y z nn n n x n y n z x y z x y z αααψ=ψψψ (10)其中,()()x yn x n y x y ααψψ、的具体表示与(5)、(6)式完全相同,z 方向的波函数为221()exp()();0,1,2,2z n z z z z z z z z N z H z n ααααψ=-== (11)2各向异性谐振子的能级简并一般情况下,各向异性谐振子的能级并不简并。
§ 2.4 一维谐振子一、能量本征方程 二、级数解法三、本征值和本征波函数平衡位置附近的微振动可近似认为是简谐振动。
例如原子核内质子和中子的振动、原子和分子的振动、固体晶格离子的振动等。
一、能量本征方程取振子的平衡位置为坐标原点22222212ˆx m x m H ω+-=d d)()(21222222x E x x m x m ψ=ψ⎥⎥⎦⎤⎢⎢⎣⎡+-ωd d因为0min =V ,∞→min out V ,所以∞<<E 0,谐振子只有束缚态,0)(lim =ψ±∞→x x 。
设ωαm =引入无量纲量 ⎪⎭⎫⎝⎛==ωλαξ 21,E x能量本征值问题转化成如下定解问题0)()()(222=ψ-+ψξξλξξd d)(lim =ψ±∞→ξξ下面会看到,束缚态条件要求λ只能取特定值,2,1,0,12=+=n n λ这导致能量的量子化。
首先把上述方程转化成可以用级数求解的形式。
考虑±∞→ξ的渐近解。
这时系数为λ的项可以忽略,方程趋近于0222=ψ-ψξξd d渐近通解为2222eeξξ-+≈ψB A ,(±∞→ξ)但因22ξe不满足束缚态的条件,所以渐近解取为22~ξ-ψe把波函数写成)(2ξξu -=ψe代入方程 0)(222=ψ-+ψξλξd d 后,求解ψ的问题则转化成求解u 的方程)1(222=-+-u uu λξξξd d d d这个方程称为Hermite 方程,可以用级数求解。
二、级数解法在原点0=ξ附近,用幂级数kk k a u ξξ∑∞==0)(代入Hermite 方程,得0)1(2)1(01122=-+--∑∑∑∞=-∞=-∞=k k kk k k k k k a ka a k k ξλξξξ把前两项的求和序号改为从0开始0)1(2)1)(2(02=-+-++∑∑∑∞=∞=∞=+k k kk k k k k k a ka a k k ξλξξ由此得到展开系数ka 的递推关系,2,1,0,)1)(2()1(22=++--=+k a k k k a k k λ只要给定0a 或者1a ,就可以把)(ξu 分成只含偶次项和只含奇次项的级数+++=+++=553312442201)()(ξξξξξξξa a a u a a a u而波函数为⎪⎩⎪⎨⎧=ψ--)()()(221222ξξξξξu u e e当∞→k 时)(1ξu 的相邻后项对前项的系数比值的极限为m k k k k a a k k 12)1)(2()1(22=→++--=+λ, ,2,1=m这与2e ξ的幂级数相邻项系数比值11+m 的极限相同。
文献综述题目:线性谐振子相图研究姓名:学号:系别:物理与电子信息工程系专业:物理学年级:指导教师:2009年2月7 日文献综述一、前言线性谐振子是量子力学中可以精确求解的有限几个事例之一[1],其中最简单的线性谐振子是简谐振子。
自然界中任何一个力学系统,只要某一个物理量在其稳定平衡点附近作微小振动,便可以用简谐振子模型来描述,例如:复摆的振动、分子的振动、晶格的振动、原子核表面振动以及辐射场的振动等。
在选择适当的坐标系之后,复杂的运动往往可以分解成若干彼此独立的一维简谐振动(simple harmonic vibration )。
简谐振动作为一种最简单最基本的振动,往往还是复杂运动的初步近似,是研究振动的基础。
因此研究它在理论上和应用上都有重大的意义。
其中从相空间的角度来研究振动系统的力学问题如今已经成为一个研究趋势。
因为相图里包含着完整的力学系统的全部信息,无须去解复杂的运动方程[2]。
计算机技术软硬件的飞速发展,为此研究趋势提供了现实条件。
本论文从简谐振子的基本定义出发,在Fortran 90条件下进行数值模拟并在Origin75 软件下获得简谐振子的相图。
二、主体2.1简谐振动的定义定义一: 物体只在弹性力或准弹性 (线性回复力)作用下发生的运动,即动力学方程为的运动为简谐振动[2]。
定义二: 在无外来强迫力作用下, 物体相对于平衡点的位移随时间按余弦(或正弦)规律变化即 则称物体作简谐振动式即简谐振动的表达式[3]。
—振幅;—角频率;—相位;—初相位。
位移随时间的变化曲线称为振动曲线。
广义定义:某个物理量随时间的变化是按正弦或余弦规律,则可称该物理量做简谐动,可用表示 。
自然界中任何一个力学系统中,只要某一个物理量在其稳定平衡点附近作微小振动,便可以用这种简谐振子模型来描述,例如:复摆的振动、分子的振动、晶格的振动,原子核表面振动、辐射场的振动以及电磁场振动等等。
2.2简谐振动的基本特征及动力学特征简谐振动位移随时间的变化 cos()x A t ωφ=+222d d xx o tω+=()cos()x t A t ωφ=+()cos()x t A t ωφ=+物体作简谐振动时,速度为:物体作简谐振动时,加速度为:可见物体做简谐振动时,其速度、加速度都以同样的角频率作简谐振动,相位依次超前π/2。
第三章 谐振子一 内容提要1 一维线性谐振子的能级与波函数2221)(x x V μω= 222212ˆˆx p Hμω+= ,3,2,1)21(=ω+=n n E n)()(2221x H eN x n x n n α-=ψ [其中 !2n N n n πα=μω=α ] 2 谐振子的升降算符 [1] 升降算符)ˆˆ(2ˆp i x aμω-μω=+ )ˆˆ(21p ix μω-α= )ˆˆ(2ˆp i x aμω+μω= )ˆˆ(21pix μω+α= 则 )ˆˆ(2ˆ++μω=a ax)ˆˆ(2ˆ+-μω-=a a i p [2] 升降算符的性质11ˆ++ψ+=ψn n n a1ˆ-ψ=ψn n n a1]ˆ,ˆ[=+a a二 例题讲解1 一维谐振子如果考虑非谐振微扰项4'ˆx Hλ=,求体系能级的一级修正。
解:>+<μωλ>=<λ>==<+n a an n x n n Hn E n 424')1()ˆˆ()2(ˆ 可以导出 )122(3)ˆˆ(24++>=+<+n n n a an 那么 =)1(n E )122()(4322++μωλn n2 已知单摆在重力作用下能在竖直平面内摆动。
求:[1] 小角度近似下,体系的能量本征值及归一化本征函数。
[2] 由于小角度近似而引起的体系基态能级的一级近似。
解:摆球平衡位置作为势能零点 摆球重力势能为)cos 1(θ-==mgl mgh V (1)[1] 由公式 -θ+θ-=θ42!41!211c o s(2)得在小角度近似下的二级修正势能为:2221))211(1(θ=θ--≈mgl mgl V (3)体系Hanmilton 为V L IV mr mv r V mv H z +=+⨯=+=ˆ21)(2121ˆ222 即:22221)(21ˆθ+θ=mgl d d i ml H(4) 当 θ≈θ=→θl l x sin 0设 lg =ω (4)可以变为22222212ˆx m dx d m H ω+= (5) (5)与一维谐振子类似,则(5)的解为:,3,2,1)21(=ω+=n n E n)()(2221x H eN x n x n n α-=ψ [其中 !2n N n n πα=μω=α ] (6) [2] )cos 1()(21ˆ22θ-+θ=mgl d d i ml H(7) 则微扰项20'21)cos 1(ˆˆˆθ-θ-=-=mgl mgl H H H (8) 以(2)式取前三项代入(8)得434'241!41ˆmgx l mgl H-=θ-= (9) 利用上题可以得到=)1(n E )122())(241(43223++ω-n n m mg l )122()(321223++ω=n n m mg l3 质量为m 的粒子处于一维谐振子势场)0(21)(21>=k kx x V 的基态[1] 如果弹性系数k 突然变为k 2,即势场变为)0()(22>=k kxx V ,随即测量粒子的能量,求发现粒子处于新势场)(2x V 的基态的概率;[2] 势场突然由)(1x V 变成)(2x V 后,不进行测量,经过一段时间τ后,势场又恢复成)(1x V ,问τ取什么值时粒子仍恢复到原来)(1x V 场的基态(概率100%)?解:[1] 粒子的波函数),(t x ψ随时间变化应满足dinger o Schreq ψ+∂ψ∂-=ψ∂∂V xm t i 2222 当V 突然改变(由)(1x V →)(2x V ),但变化量有限时ψ仍然是t 的连续函数,即V 突变时ψ不变。