§2-3 波函数及其物理意义经典物理中,某一种波是用相应的物理量随时间和空间的变化的关系式来描述的。频率为ν,波长为λ,沿x 方向的平面波可表示为:
[]
cos 2()A r n t Φ=?πλν i n K 对三维情况,沿方向传播的平面波,其函数形
式可表示为:
[]
cos 2()A x t πλνΦ=?[]
exp 2()A i r n t Φ=?πλν i 通常写成复指数形式
为了表示微观粒子的波粒二象性,可以用平面波来描写自由粒子
[]
exp 2()A i r n t Ψ=?πλν i νE h λh p ==此平面波的频率和波长与自由粒子的能量和动量由德布罗意关系联系起来。
exp ()i A p r Et ??Ψ=????? i 2h π
= 其中此式称为自由粒子的波函数。
如果粒子处于随时间和位置变化的力场中运动,他的动量和能量不再是常量(或不同时为常量),粒子的状态就不能用平面波描写,而必须用较复杂的波描写,一般记为:
()()Ψ,,,Ψx y z t r,t K 或描写粒子状态的波函数,它通常是一个复函数。
电子的衍射实验电子源感
光
屏
P
P O Q Q O 玻恩对波函数的统计解释——几率波
结论:衍射实验所揭示的电子的波动性是:许多电子在同一个实验中的统计结果,或者是一个电子在许多次相同实验中的统计结果。
1.入射电子流强度大,很快显示衍射图样。
在衍射极大值得地方,波的强度(波函数振幅的平方)大,在衍射极小值的地方,波的强度很小或等于零。即电子波衍射的强弱正比于波函数振幅的平方()()()
2Ψ,Ψ,Ψ,r t r t r t ?=K K K 2.入射电子流强度小,开始
显示电子的微粒性,没有
分布规则,长时间亦显示
相同的衍射图样。
在衍射极大值得地方,电子出现的几率大,在衍射极小值的地方,电子出现的几率小或等于零。即电子波衍射的强弱正比于电子出现的几率。
波函数意义的统计解释:
波函数在空间某点的强度(振幅绝对值的平方)和在这点找到粒子的几率成比例。描写粒子的波是几率波,波函数Ψ(r ,t )有时也称为几率幅。这就是首先由玻恩提出的波函数的几率解释,它是量子力学的基本原理。
根据波函数的几率解释,假如知道了描写微观体系的波函数,由波函数振幅绝对值的平方,就可以得到粒子在空间任意一点出现的几率,也就可以得到体系的各种性质,因此我们说波函数描写体系的状态,所以波函数也称为态函数。
设波函数描写粒子的状态,在空间一点和时刻t ,波的强度为()Φ,r t K 2ΦΦΦ
?=r K (),dP r t K 用表示微观粒子在空间处体积元内出现的几率,则
r K d τ()()22,Φ,dP r t A r t d τ
=K K
微观粒子在时刻t ,空间处单位体积内出现的几
率为r K ()()()22,,Φ,dP r t P r t A r t d τ
==K K K 几率密度
()()22,Φ,V V P r t d τA r t d τ=∫
∫K K 在体积V 内,t 时刻找到粒子的几率为:
由于粒子在空间总要出现,所以在全空间找到粒子的几率应为一。()()22,Φ,1P r t d τA r t d τ∞∞
==∫∫K K ()12
21Φ,A r t d τ∞????=????
∫K 可得令()()Φ,Ψ,A r t r t =K K ()()()()2,Ψ,Ψ,Ψ,1P r t d τr t d τr t r t d τ?∞∞∞
==∫∫∫K K K K =()Ψ,r t K 称为归一化波函数,A 为归一化常数。
波函数的归
一化条件
()()2,Ψ,P r t r t =K K 几率密度为:
由于粒子在全空间出现的几率等于一,所以粒子在空间各点出现的几率只取决于波函数在空间各点强度的相对比例,而不取决于强度的绝对大小。在t 时刻,空间任意两点r 1和r 2处找到粒子的相对几率之比
222
111222Φ(,)Φ(,)Ψ(,)Φ(,)Φ(,)Ψ(,)r t A r t r t r t A r t r t ==G G G G G G 所描写状态的相对几率是相同的,所以描述的是同一状态。
()()()Φ,Ψ,Φ,r t r t A r t K K K 和=因而,将波函数乘上一个常数后,所描写的粒子状态不变。
这与经典波不同。经典波波幅增大一倍(原来的2 倍),则相应的波动能量将为原来的4 倍,因而代表完全不同的波动状态。经典波无归一化问题。
解:令已归一化波函数为Ψ(),Ψ()Φ()
x x A x =28A a
=22A a =2
22012Ψ()cos 4a πx x dx A dx
a +∞?∞=∫∫22041cos 1114242
a πx a a A dx A +===∫)
,0(a x ∈a x x π2cos 2
1)(=Φ将其归一化。例题:设波函数所以归一化波函数为()22Ψcos πx x a a
=课堂练习:
设,为常数,求归一化系数A 。()222x x Ae
αψ?=α解:()22221
x x dx A e dx αψ∞
∞??∞?∞==∫∫2221
x A e dx
α∞
??∞=∫2y e dy π∞??∞=∫12A απ??=????2
A απ=
综上所述,的物理意义为微观粒子在t 时刻出现在空间处单位体积内
的几率,微观粒子的波动性表现在几率密度的波动。()()()2*Ψ,Ψ,Ψ,r t r t r t =K K K r K ()
2Ψ,r t K 利用波函数求力学量平均值
已知波函数,可得微观粒子在t 时刻出现在空间处单位体积内的几率,即可利用各个物理量在值求得该物理量在t 时刻的平均
值,因此可以完全描述粒子的状态。
()Ψ,r t K r K ()2Ψ,r t K r K
在统计学中知道,
当可能值为离散值时: 一个物理量的平均值等于物理量出现的各种可能值乘上相应的几率求和;当可能值为连续取值时:一个物理量出现的各种可能值乘上相应的几率密度求积分。
(1)坐标平均值
()2*(,)|(,)|(,)(,)x t xP r t d x r t d r t x r t d ∞∞?∞?∞∞
?∞=
=Ψ=ΨΨ∫∫∫τττ
(2)力学量f (r )平均值
()(,)(,)f r t f r r t d τ?
=ΨΨ∫∫∫
3/22 0() 0 0
x xe x x x ??≥=?
例题:设粒子波函数为(2)求x 2的平均值
解:(1)3320()()4λx x ψx x ψx dx λx e dx ∞∞
?
??∞==∫∫(2)22
342()()4λx x ψx x ψx dx λx e dx ∞∞???∞?∞==∫∫32λ=23λ=10!n αx n n x e dx α
∞
?+=∫0301(,,)r a r e a ?=ψθ?π(1)求r 的平均值
例题:氢原子处于状态(2)求势能的平均值2e r ?解:(1)02223000
0(,,)(,,)1sin r ππa r ψr θφr ψr θφd τre r θdrd θd φπa ∞
??∞∞?==∫∫∫∫023043300003443!2
2r a a r e dr a a a ∞?===??????
∫10!n αx n n x e dx α∞?+=∫
(2)0022222300002
2230000
(,,)(,,)1sin sin r ππa r ππa e U ψr θφψr θφd τr e e r θdrd θd φπa r e e r θdrd θd φπa ∞??∞∞?∞???=?????
=?=?∫∫∫∫∫∫∫02222233000004412r a e e e re dr a a a a ∞?=?=?=???????
∫10!n αx n n x e dx α
∞?+=∫
§2-4 薛定谔方程
微观粒子量子状态用波函数完全描述,波函数确定之后,粒子的任何一个力学量的平均值及其测量的可能值和相应的几率分布也都被完全确定,波函数完全描写微观粒子的状态。因此量子力学最核心的问题就是要解决以下两个问题:
(1)在各种情况下,找出描述系统的各种可能的波函数;
(2)波函数如何随时间演化。
这些问题在1926年薛定谔提出了波动方程之后得到了圆满解决。
1、方程的引入1、方程的引入
(一)自由粒子满足的方程
Ψexp ()exp ()x y z i A p r Et i A p x p y p z Et ??=??????
??=++?????
G G ==描写自由粒子波函数:
将上式对t 微商,得:
1i E i E t t ?Ψ?=?Ψ→Ψ=Ψ?? ()将Ψ对坐标二次微商,得:
()ΨΨx y z i p x p y p z Et x i Ae p x x ++???==??==
2222ΨΨx p x ?=??=
22222222ΨΨΨ,Ψy z p p y z ??=?=???==
同理有Ψ++?=?Ψ?+?Ψ?+?Ψ?][12222222222z y x p p p z y x
引入拉普拉斯算符2222222x y z
????=++???Ψ++?=?Ψ?+?Ψ?+?Ψ?][12222222222z y x p p p z y x
2
221ΨΨp ?=?=上式写为22ΨΨ2i t μ?=???==()
22
2ΨΨ222p μμ
??==上式两边乘以得:22μ?=22p E T μ==对于自由粒子22ΨΨ2E μ
??==即为自由粒子满足的方程
()
ΨΨ1i E t ?=?=
(二)势场U (r ) 中运动的粒子
若粒子处于势场中运动,则2(,)2p E U r t μ
=+G (,)U r t G ()
22
2ΨΨ222p μμ??==()2
2ΨΨ,Ψ2i U r t t μ?=??+?=K =()222ΨΨΨ,Ψ22p i E U r t t μμ???+?=?=????
?=K =整理得该方程称为薛定谔方程,也常称为波动方程。
()ΨΨ1i E t ?=?=两式相减得:2、波函数的标准条件2、波函数的标准条件
2、根据波函数的统计解释,要求波函数的模的平方是单值的,从而保证几率密度在如何时刻都是确定的。
1、归一化条件要求,在空间任何有限体积元中找到粒子的几率为有限值。
()()22,Φ,V
V P r t d τA r t d τ==∫∫K K 有限值
3、波函数对坐标的二阶偏微商要求波函数Ψ(r, t)及其各阶微商的连续性。
概括之,波函数在全空间每一点通常应满足单值、有限、连续三个条件,该条件称为波函数的标准条件。
3、态叠加原理
3、态叠加原理
同光学中波的叠加原理一样,量子力学中也存在波叠加原理。因为量子力学中的波函数决定体系的状态,称为状态波函数,所以量子力学的波叠加原理称为态叠加原理。
一般情况下,如果Ψ1和Ψ2是体系的可能状态,那末它们的线性叠加Ψ= C 1Ψ1+ C 2Ψ2 也是该体系的一个可能状态,其中C 1和C 2是复常数,这就是量子力学的态叠加原理。
态叠加原理一般表述:
若Ψ1,Ψ2 ,..., Ψn ,...是体系的一系列可能的状态,则这些态的线性叠加Ψ= C 1Ψ1+ C 2Ψ2+ ...+ C n Ψn + ... (其中C 1, C 2 ,...,C n ,...为复常数)。也是体系的一个可能状态。
处于Ψ态的体系,部分的处于Ψ1态,部分的处于Ψ2态...,部分的处于Ψn ,...
2
2Ψ(,)[()]Ψ(,)2i r t U r r t t μ?=??+?G G G ==)
()(),(t f r t r ψ=Ψ现在让我们讨论有外场情况下,且与不显含时间的薛定谔方程:令:U (r )与t 无关时,可以分离变量
代入22()()()[]()2d i ψr f t f t U ψr dt
μ=??+G G ==4、定态薛定谔方程4、定态薛定谔方程
()()
ψr f t G 两边同除
等式两边是相互无关的物理量,故应等于与t , r 无关的常数,设为E 。
2211()[]()()()2d i f t U ψr f t dt ψr μ
=??+G ==G 22()()[]()()2d i f t E f t d t U ψr E ψr μ
?=??????+=??=G G = /~)(iEt e t f ?()()(,)()i Et r t r f t r e
ψψ?Ψ== 于是:也应适用于自由粒子的波函数,对比自由粒子的波
函数,由德布罗意关系可
知:()()(,)()i Et r t r f t r e
ψψ?Ψ== exp ()i A p r Et ??Ψ=?????
i E 就是体系处于波函数Ψ(r ,t )所描写的状态时的能量。也就是说,此时体系能量有确定的值,所以这种状态称为定态,波函数Ψ(r ,t )称为定态波函数。
该方程称为定态薛定谔方程,ψ(r )也可称为定态波函数,或可看作是t =0时刻ψ(r ,0)的定态波函数。
22[]()()2U ψr E ψr μ
??+=G G =空间波函数ψ(r )可由方程
和具体问题波函数应满足的条件得出。
对于微观粒子的三维情况,定态薛定谔方程的直角坐标形式为:
2
222222[]()()2U ψr E ψr μx
y z ??????+++=???????G G =相应的球坐标形式为:
22222222111[sin 2sin sin ]()()
r θμr r r r θθθr θφU ψr E ψr ????????????++?????????????????+==G G
(,)ΨΨn n n
P r t ?=G 定态在空间几率密度与时间无关
()()[exp(/)][exp(/)]n n n n ψr iE t ψr iE t ?=??K K ==()()exp(/)exp(/)
n n n n ψr iE t ψr iE t ?=?K K ==)
()(r r n n ψψ?=综上所述,当Ψ满足下列三个等价条件中的任何一个时,Ψ就是定态波函数:
1. Ψ描述的状态其能量有确定的值;
2. Ψ满足定态薛定谔方程;
3. |Ψ|2 与t 无关。
关于薛定谔方程 一.定义及重要性 薛定谔方程(Schrdinger equation)是由奥地利物理学家薛定谔提 出的量子力学中的一个基本方程,也是量子力学的一个基本假定, 其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合 建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都 有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式 以及对应的能量,从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程,亦是量子力学的一个基 本假定,它的正确性只能靠实验来检验。 二.表达式 三.定态方程 ()() 2 2 2 V r E r m η ψψ + ?? -?= ?? ?? 所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。 其中,E是粒子本身的能量;v(x,y,z)是描述势场的函数,假设不随时间变化。
2 2 22222 z y x ??????++=? 可化为 d 0)(222 =-+ψψv E h m dx 薛定谔方程的解法 一. 初值解法;欧拉法,龙格库塔法 二. 边值解法;差分法,打靶法,有限元法 龙格库塔法(对欧拉法的完善) 给定初值问题 ). ()()((3) ) ,(),()( ,,(2) )(),( 311212 2111021h O t y t y hk y h t f k y t f k k c k c h y y y c c a y b t a y t f dt dy i i i i i i i i =-???????++==++==?????=≤≤=++的局部截断误差使以下数值解法的值及确定常数ββα βα
波函数和薛定谔方程-力学量算符 1.一维运动的粒子处在 的状态,其中,求: (1)粒子动量的几率分布函数; (2)粒子动量的平均值。 [解]首先将归一化,求归一化系数A。 (1)动量的几率分布函数是 注意到中的时间只起参数作用,对几率分布无影响,因此可有 令 代入上式得 (2) 动量p的平均值的结果从物理上看是显然的,因为对本题说来,粒子动量是和是的几率是相同的。讨论: ①一维的傅里叶变换的系数是而不是。 ②傅里叶变换式中的t可看成参变量。因此,当原来坐标空间的波函数不含时间变量时, 即相当于的情况,变换式的形式保持不变。
③不难证明,若是归一化的,则经傅里叶变换得到也是归一化的。 2.设在时,粒子的状态为 求粒子动量的平均值和粒子动能的平均值。 [解]方法一:根据态迭加原理和波函数的统计解释。任意状态总可以分解为单色平面波的线性和,即,展开式的系数表示粒子的动量为p时的几率。知道了几率分布函数后,就可按照 求平均值。 在时,动量有一定值的函数,即单色德布罗意平面波为,与的展开式比较可知,处在状态的粒子动量可以取 ,而,粒子动量的平均值为 A可由归一化条件确定 故 粒子动能的平均值为 。 方法二:直接积分法
根据函数的性质,只有当函数的宗量等于零时,函数方不为零,故的可能值有 而 则有及。 讨论:①由于单色德布罗意平面波当时不趋于零,因此的归一化积分是发散的,故采用动量几率分布的概念来求归一化系数。 ②本题的不是平方可积的函数,因此不能作傅氏积分展开,只能作傅氏级数展开,即 这时对应于波函数的是分立谱而不是连续谱,因此计算积分,得到函数。 ③在连续谱函数还未熟练以前,建议教学时只引导学生按方法一做,在第三章函数讲 授后再用函数做一遍,对比一下,熟悉一下函数的运算。 3.一维谐振子处在 的状态,求: (1)势能的平均值; (2)动量的几率分布函数; (3)动能的平均值 [解]先检验是否归一化。 是归一化的。 (1) . 其中应用及 (2)由于是平方可积的,因此可作傅氏变换求动量几率分布函数
实验三 定态薛定谔方程的矩阵解法 一.实验目的 1.掌握定态薛定谔方程的矩阵解法。 2.掌握几种矩阵特征值问题数值解法的原理,会调用相应的子程序求解具体问题。 二.实验内容 1.问题描述 以/2ω/()m ω为长度单位,一维谐振子的哈密顿量为 2 202d H x dx =-+, 其本征值为21n E n =+,本证波函数为 2 /2)()n n x H x ?=-, 其中()n H x 为厄米多项式,满足递推关系 11()2()2()n n n H x xH x nH x +-=-。 用矩阵方法求 2 22d H x x dx =-++ 的本证能量和相应的波函数。 2.问题分析 H E ψψ= 0()|j j j t c ψ?∞ ==>∑ 0||i i j i j i j c E c x Ec ??∞ =+<>=∑ 11|j j j x ???-+>=>>
11||||j j j j x x ????-+<>= <>= 0010010 112111,211,11,1 n n n n n n n n n n n n E x c c x E x c c E x E x c c x E c c -------?????????????????????????=??????????????????????? ? 3.程序编写 子程序及调用方法见《FORTRAN 常用算法程序集(第二版)》第三章 徐士良,P97 4.实验要求 ◆用恰当的算法求解以上实对称三对角矩阵的特征值问题。 ◆取n=8,给出H 的全部特征值和相应的特征向量。 5.实验步骤 ● 启动软件开发环境Microsoft Developer Studio 。 ● 创建新工作区shiyan03。 ● 创建新项目xm3。 ● 创建源程序文件xm3.f90,编辑输入源程序文本。 ● 编译、构建、运行、调试程序。 6.实验结果 程序设计:
一维定态薛定谔方程 的应用 授课人: 物理科学与技术学院
势 阱 日常生活中的各种井(阱) 物理学中研究微观粒子运动状态时常用的模型,因其势能函数曲线的形状如同井而得名 水井 窨井 陷阱 U x O a U
() U x x O a ∞ ∞00()0 , x a U x x x a ≤≤?=?∞<>? 这是一个理想化的物理模型, 应用定态薛定谔方程求解波函数, 有利于进一步理解在微观系统中 能量量子化和概率密度等概念 这样的势能函数称为 一维无限深势阱
建立定态薛定谔方程并求解 假设微观粒子质量为 ,由 m 22 2d ()()()2d U x x E x m x ψψ??-+=???? x a U x 0()0≤≤=阱内( ) : 22 2d ()()2d x E x m x ψψ-= x x a U x 0 , ()<>→∞ 阱外( ): 令: 2 22mE k =得通解: ()sin() x A kx ψ?=+ 微观粒子的能量不可能达到 无穷大,所以粒子不可能在阱外出现,或者说粒子在阱外出现的概率为零。 ()0 x ψ≡222 d 0d k x ψψ+=
利用标准条件确定 和 k ?因 在整个 轴上必须连续 x ()x ψsin() 0()0 0 0 A kx x a x x x ?ψ+≤≤?=? <>?,(0)sin 0 A ψ?== a A ka ()sin()0 ψ?=+=求归一化的波函数 一维无限深势阱中 微观粒子的波函数 2220π()d sin d a n x x A x x a ψ+∞-∞=??221 A a =?= 2A a = n a x x a x a x x a π2sin 0()00 , ψ? ≤≤?=??<>?() π ()sin 1,2,3n x A x n a ψ==??, 0?=π n k a =()1,2,3n =???,
波函数和薛定谔方程 一、波函数的统计解释、叠加原理和双缝干涉实验 微观粒子具有波粒二象性<德布罗意假设); 德布罗意关系<将描述粒子和波的物理量联系在一起) 物质波<微观粒子—实物粒子) 引入波函数<概率波幅)—描述微观粒子运动状态 对于微观粒子来说,如果不考虑“自旋”一类的“内禀”态,单值波函数是其物理状态的最详尽描述。至少在目前量子力学框架中,我们不能获得比波函数更多的物理信息。b5E2RGbCAP 微观粒子的状态用波函数完全描述 ——量子力学中的一条基本原理 该原理包含三方面内容:粒子的状态用波函数表示、波函数的统计解释和对波函数性质的要求。 要明确“完全”的含义是什么。按着波函数的统计解释,波函数统计性的描述体系的量子态,若已知单粒子<不考虑自旋)波函数,则不仅可以确定粒子的位置概率分布,而且如动量等粒子的其它力学量的概率分布也均可通过波函数而完全确定。由此可见,只要已知体系的波函数,便可获得该体系的一切物理信息。从这个意义上说,有关体系的全部信息已包含在波函数中,所以说微观粒子的状态用波函数完全描述。p1EanqFDPw 必须强调指出,波函数给出的有关粒子的“信息”本质上是统计性质的。例如,在适当条件下制备动量为p的粒子,然后测量其空间位置,我们根本无法预言测量的结果,我们只能知道获得各种可能结果的概率。DXDiTa9E3d
很自然,人们会提出这样的疑问:既然量子力学只能给出统计结果,那就只需引入一个概率分布函数<象经典统计力学那样),何必假定一个复值波函数呢?RTCrpUDGiT 事实上,引入复值波函数的物理基础,乃是量子力学中的又一条基本原理——叠加原理。 这条原理告诉我们,两种状态的叠加,绝不是概率相加,而是带有相位的复值波函数的叠加<数学求和)。正因如此,在双缝干涉实验中,我们才能看见屏上的干涉花纹。5PCzVD7HxA 实物粒子双缝干涉实验分析 我们首先只打开一条狭缝,根据粒子的波动性,可以预言屏上将显示波 长<为粒子动量)的单缝衍射花纹。但是,根据粒子的微粒性, 它们将是一个一个打上去的,怎样将这两种性质的描述调和起来呢?为此,我们想象将入射粒子束强度降低,直到只一个粒子通过狭缝,这时屏上会出现很微弱的衍射花纹吗?当然不会!单个粒子只能作为一个不可分割的整体打到屏上的一个点,从而出现一个小斑点。如果让这种微弱的粒子束<几乎让粒子一个一个地通过狭缝)长时间照射狭缝<相当于一个粒子的多次行为),结果发现,屏上一个一个斑点逐渐增加,最后形成一种接近连续的分布,它恰恰就是单缝衍射花纹!<单个粒子具有波动性的有力证明)jLBHrnAILg 这提示:粒子的波动性只是一种“概率波”,或者干脆说只是一种概率分布而已。这种看法对吗?这种说法容易造成误解,因为它忽略了叠加原理的要求。xHAQX74J0X 为了说明这一点,我们继续分析双缝干涉实验。
薛定谔方程 在一维空间里,一个单独粒子运动于位势中的含时薛定谔方程为 ;(1) 其中,是质量,是位置,是相依于时间的波函数,是约化普朗克常数,是位势。类似地,在三维空间里,一个单独粒子运动于位势中的含时薛定谔方程为 。(2) 假若,系统有个粒子,则波函数是定义于 -位形空间,所有可能的粒子位置空间。用方程表达, 。 其中,波函数的第个参数是第个粒子的位置。所以,第个粒子的位置是。 不含时薛定谔方程 不含时薛定谔方程不相依于时间,又称为本征能量薛定谔方程,或定态薛定谔方程。顾名思义,本征能量薛定谔方程,可以用来计算粒子的本征能量与其它相关的量子性质。 应用分离变量法,猜想的函数形式为 ; 其中,是分离常数,是对应于的函数.稍回儿,我们会察觉就是能量. 代入这猜想解,经过一番运算,含时薛定谔方程 (1) 会变为不含时薛定谔方程: 。 类似地,方程 (2) 变为
。 历史背景与发展 爱因斯坦诠释普朗克的量子为光子,光波的粒子;也就是说,光波具有粒子的性质,一种很奇奥的波粒二象性。他建议光子的能量与频率成正比。在相对论里,能量与动量之间的关系跟频率与波数之间的关系相同,所以,连带地,光子的动量与波数成正比。 1924年,路易·德布罗意提出一个惊人的假设,每一种粒子都具有波粒二象性。电子也有这种性质。电子是一种波动,是电子波。电子的能量与动量决定了它的物质波的频率与波数。1927年,克林顿·戴维和雷斯特·革末将缓慢移动的电子射击于镍晶体标靶。然后,测量反射的强度,侦测结果与X射线根据布拉格定律 (Bragg's law) 计算的衍射图案相同。戴维森-革末实验彻底的证明了德布罗意假说。 薛定谔夜以继日地思考这些先进理论,既然粒子具有波粒二象性,应该会有一个反应这特性的波动方程,能够正确地描述粒子的量子行为。于是,薛定谔试着寻找一个波动方程。哈密顿先前的研究引导著薛定谔的思路,在牛顿力学与光学之间,有一种类比,隐蔽地暗藏于一个察觉里。这察觉就是,在零波长极限,实际光学系统趋向几何光学系统;也就是说,光射线的轨道会变成明确的路径,遵守最小作用量原理。哈密顿相信,在零波长极限,波传播会变为明确的运动。可是,他并没有设计出一个方程来描述这波行为。这也是薛定谔所成就的。他很清楚,经典力学的哈密顿原理,广为学术界所知地,对应于光学的费马原理。借着哈密顿-雅可比方程,他成功地创建了薛定谔方程。薛定谔用自己设计的方程来计算氢原子的谱线,得到了与用玻尔模型计算出的能级相同的答案。 但是,薛定谔对这结果并不满足,因为,索末菲似乎已经正确地计算出氢原子光谱线精细结构常数的相对论性的修正。薛定谔试着用相对论的能量动量关系式,来寻找一个相对论性方程(现今称为克莱因-高登方程),可以描述电子在库仑位势的量子行为。薛定谔计算出这方程的定态波函数。可是,相对论性的修正与索末菲的公式有分歧。虽然如此,他认为先前非相对论性的部分,仍旧含有足够的新结果。因此,决定暂时不发表相对论性的修正,只把他的波动方程与氢原子光谱分析结果,写为一篇论文。1926年,正式发表于物理学界[2]。从此,给予了量子力学一个新的发展平台。 薛定谔方程漂亮地解释了的行为,但并没有解释的意义。薛定谔曾尝试解释代表电荷的密度,但却失败了。1926年,就在薛定谔第四篇的论文发表之后几天,马克斯·玻恩提出概率幅的概念,成功地解释了的物理意义[3]。可是,薛定谔本人一直不承认这种统计或概率的表示方法,和它所伴随的非连续性波函数坍缩。就像爱因斯坦的认为量子力学是基本为确定性理论的统计近似,薛定谔永远无法接受哥本哈根诠释。在他有生最后一年,他写给马克斯·玻恩的一封信,薛定谔清楚地表明了这看法。 含时薛定谔方程导引
第一章 薛定谔方程 §1.1.波函数及其物理意义 1. 波函数: 用波函数描述微观客体的运动状态。 例:一维自由粒子的波函数 推广 :三维自由粒子波函数 2. 波函数的强度——模的平方 3. 波函数的统计解释 用光栅衍射与电子衍射对比的方式理解波函数的统计解释。 t 时刻,出现在空间(x,y,z )点附近单位体积内的粒子数与总粒子数之比。 t 时刻,粒子出现在空间(x,y,z )点附近单位体积内的概率。 t 时刻,粒子在空间分布的概率密度 4、 波函数的归一化条件和标准条件 归一化条件 粒子在整个空间出现的概率为1 标准条件:一般情况下, 有关特殊情况波函数所满足的条件参看曾谨言教程。 对微观客体的数学描述: 脱离日常生活经验,避免借用经典语言引起的表观矛盾 §1.2. 薛定谔方程 是量子力学的基本假设之一,只能建立,不能推导,其正确性由实验检验。 1. 建立 (简单→复杂, 特殊→一般) 一维自由粒子的振幅方程 非相对论考虑 2. 一维定态薛定谔方程 2 |),,,(|t z y x ψ1d d d d d ||2===?=ψ???N N N N V V N N V V V . 是单值、有限、连续的ψ0)(2d )(d 222=ψ+ψx mE x x 0)()(2d )(d 222=ψ-+ψx U E m x x
3. 三维定态薛定谔方程 4. 一般形式薛定谔方程 5. 多粒子体系的薛定谔方程 讨论: 1、薛定谔方程也称波动方程,描述在势场U 中粒子状态随时间的变化规律。 2 、建立方程而不是推导方程,正确性由实验验证。薛定谔方程实质上是一种基本假设,不能从其他更基本原理或方程推导出来,它的正确性由它解出的结果是否符合实验来检验。 3、薛定谔方程是线性方程。是微观粒子的基本方程,相当于牛顿方程。 4、自由粒子波函数必须是复数形式,否则不满足自由粒子薛定谔方程。 5、薛定谔方程是非相对论的方程。 量子力学的中心任务就是求解薛定谔方程。 求解问题的思路: 1. 写出具体问题中势函数U (r )的形式代入方程 2. 用分离变量法求解 3. 用归一化条件和标准条件确定积分常数 4. 讨论解的物理意义, 薛定谔的另一伟大科学贡献 《What is life ?》 薛定谔(Schroding,1897-1961)奥地利人,因发现原子理论的有效的新形式一波动力学与狄拉克(Dirac,1902-1984)因创立相对论性的波动方程一狄拉克方程,共同分享了1933年度诺贝尔物理学奖 定态薛定谔方程 一.定态薛定谔方程条件:V (r,t )=V(r), 与t 无关。用分离变量法, 令Ψ=φ(r)f(t),代入薛定谔方程,得两个方程: 此称定态薛定谔方程 整个定态波函数形式: ),,,(),,,()],,,(2[),,,(2121212221t r r t r r V t r r m t r r t i i i i ψ+ψ?-=ψ??∑)t (Ef t )t (f i =?? Et i ce )t (f -=)r (E )r ()r (V )r (m ?=?+??-222Et i e )r ( -?=ψ
薛定谔方程应用举例II---原子系统
? 氢原子 ? 电子自旋 ? 多电子原子
1
氢原子的定态薛定谔方程
?原子由一个原子核和核外电子构成,属于多粒子体系。多粒 子体系的总能量等于每一个粒子的能量与粒子间相互作用能量 之和。
?氢原子包括一个原子核和电子,库仑场是各向同性的,哈密 顿量可记作(绝热近似):
H?
=
?
h2 2me
?2
+
qeU(r)
me为电子质量,qe是电子电荷。U(r)为原子核静电场中的库 仑势,记作:
U(r) = ? Zqe = ? Z h2
4πε0r a1meqer
Z为核的电荷数,a1 = 4πε0?2/(meqe2) = 0.529?,为氢原子的第
一波尔轨道半径。
2
??? ?
h2 2me
?2
?
Zh 2 a1meqer
??ψ
?
(r)
=
E
?ψ
(r)
中心力场问题,采用球坐标,薛定谔方程为:
? ?? ??
h2 2me
?
????
1 r2
? ?r
r2
? ?r
?
L?2 r2
???? ?
Zh2
?
?ψ (r,?,θ ) =
a1mer ??
E ?ψ (r,?,θ )
用分离变量法求解,令:
ψ (r,θ ,φ) = R(r) ?Y (?,θ )
分别求解径向波函数R(r)和角向波函数Y(?,θ)。
3
[键入作者姓名] [键入文档标题] ——利用分步快速傅里叶变换对光纤中光信号的传输方程进行数值求解
1、非线性薛定谔方程 非线性薛定谔方程(nonlinear Schrodinger equation ,NLSE)是奥地利物理学家薛定谔于1926 年提出的,应用在量子力学系统中。由于量子力学主要研究粒子的动力学运动状态,所以不能运用牛顿力学公式来表示。通常在量子力学中,研究系统的状态一般通过波函数(x ,t)来表示。而对波函数的研究主要是求解非线性薛定谔方程。本文主要研究光脉冲在光纤中传输状态下的演变。 一般情况下,光脉冲信号在光纤中传输时,同时受到光纤的色散和非线性效应的影响。通过Maxwell 方程,考虑到光纤的色散和非线性效应,可以推导出光信号在光纤中的传输方程,即非线性薛定谔方程。NLSE 是非线性偏微分方程,一般很难直接求出解析解,于是通过数值方法进行求解。具体分为两大类:(1)分布有限差分法(split-step finite differencemethod ,SSFD);(2)分步傅里叶变换法(split-step Fourier transform method ,SSFT)。一般情况,在达到相同精度,由于分步傅里叶变换法采用运算速度快的快速傅里叶变换,所以相比较有限差分法运算速度快一到两个数量级。于是本文介绍分步傅里叶变换法来对光纤中光信号的传输方程,即非线性薛定谔方程进行数值求解。并通过MATLAB 软件对结果数值仿真。 非线性薛定谔方程的基本形式为: 22||t xx iu u u u =+ 其中u 是未知的复值函数. 目前,采用分步傅立叶算法(Split step Fourier Method)求解非线性薛定谔方程的数值解应用比较多。分步傅立叶方法最早是在1937年开始应用的,这种方法己经被证明是相同精度下数值求解非线性薛定愕方程最快的方法,部分原因是它采用了快速傅立叶变换算法(Fast Fourier Transform Algorithm)。基于MATLAB 科学计算软件以及MATLAB 强大的符号计算功能,完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 一般情况下,光脉冲沿光纤传播时受到色散和非线性效应的共同作用,假设当传输距离 很小的时候,两者相互独立作用,那么,根据这种思想可建立如下分步傅立叶数值算法的数 学模型: 把待求解的非线性薛定谔方程写成以下形式: ??()U D N U z ?=+? (I ) (II )
非线性薛定谔方程的孤子解和怪波解 摘要:光纤中光波的传输模型一直是当前研究的热点理论模型之一,从非线性薛定谔方程到金格堡-朗道方程,都试图对其进行更好的阐释,其次对于非线性动力学系统中,非线性薛定谔方程的解有呈现出非常多有趣的特征,对于其中特定解的研究能够让我们了解脉冲演化的本质,所以本文主要从孤子解的传输入手,并且简单介绍了怪波解的解形式。 薛定谔方程又称薛定谔波动方程,是量子力学的一个基本方程,同时又是量子力学的基本假设之一,由奥地利物理学家薛定谔1926年在《量子化就是本征值问题》中提出的,它在量子力学中的地位非常重要,相当于牛顿定律对于经典力学一样。 随着人们对世界的不断探索,非线性现象逐渐走进人们的视野,这种现象一般大都用非线性偏微分方程的数学模型来描述,显然线性方程已经不能满足人们的需求。 1973年,Hasegawa从含有非线性项的色散方程中推导出了非线性薛定谔方程。非线性薛定谔方程(NLS)是普适性很强的一个基本方程,最简单的形式是: 其中为常数。因为这个方程在几乎所有的物理分支及其他科学领域得到了广泛的应用,如超导,光孤子在光纤中传播,光波导,等离子体中的Langnui波等,所以许多学者对此方程的研究投入了很大的热情,至今还在生机勃勃的向前发展着。 1 分步傅里叶法计算演化过程 对于处理非线性性薛定谔方程,常用的数值仿真方式为分步傅里叶方法,为了简单起见,只考虑二阶色散和自相位调制,不考虑高阶色散、自陡以及四波混频等高阶非线性效应。上述方程中做 2 β为二阶色散,γ表示Kerr效应系数,g和α分别代表光纤中的增益和损耗。对上述方程转化到频域,先不考虑增益和损耗。可以得到 2 k k k k k dA i A i a a dz βγ =?+F. 其中2 2 2 k i β β ?=Ω 令() exp k k A B i z β =?可以得到 () 2exp k k k k dB i a a i z dz γβ =-? F 以上方程可以用四阶龙格库塔直接求解,但是速度较慢,所以我们需要做差分处理。 ()() ()()() 2 exp k k k k k B z z B z i a z a z i z z γβ +?- =-? ? F 再利用() exp k k A B i z β =?可以得到 ()()()() ()()() 2 2 exp exp exp k k k k k k k k A z z A i a z a z z i z a z i a z z i z γβ γβ ?? +?=+??? ?? ?? ?? ≈????? ?? F F 然后做傅里叶反变换就可以得到最终的结果 ()()()() 2 1exp exp - k k k k a z z a z i a z z i z γβ ?? +?=????? ?? F F
关于薛定谔方程 一. 定义及重要性 薛定谔方程(Schrdinger equation )是由奥地利物理 学家薛定谔提出的量子力学中的一个基本方程,也是量子力学的一个基本假定,其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式以及对应的能量,从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程,亦是量子力学的一个基本假定,它的正确性只能靠实验来检验。 二. 表达式 三. 定态方程 ()()2 22V r E r m ηψψ+??-?=???? 所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。 其中,E 是粒子本身的能量;v(x ,y ,z)是描述势场的函数,假设不随时间变化。
2 2 22222z y x ?? ????++=? 可化为d 0)(222=-+ψψ v E h m dx 薛定谔方程的解法 一. 初值解法;欧拉法,龙格库塔法 二. 边值解法;差分法,打靶法,有限元法 龙格库塔法(对欧拉法的完善) 给定初值问题 ).()()((3) ) ,() ,() ( ,,(2) )() ,( 3112122111021h O t y t y hk y h t f k y t f k k c k c h y y y c c a y b t a y t f dt dy i i i i i i i i =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα
第二章波函数和薛定谔方程 ●§2.1 波函数的统计解释 ●§2.2 态叠加原理 ●§2.3 薛定谔方程 ●§2.4 粒子流密度和粒子数守恒定律●§2.5 定态薛定谔方程 ●§2.6 一维无限深势阱 ●§2.7 线性谐振子 ●§2.8势垒贯穿
本章主要介绍了波函数的统计解释、薛定谔方程的建立过程、用定态薛定方程处理势阱问题和线性谐振子问题。
§2.1 波函数的统计解释(一)波函数 (二)波函数的解释 (三)波函数的性质
?? ????-?=ψ)(exp Et r p i A ?3个问题? 描写自由粒子的 平面波 ),(t r ψ?如果粒子处于随时间和位置变化的力场中运动,他的动量和能量不再是常量(或不同时为常量)粒子的状态就不能用平面波 描写,而必须用较复杂的波描写,一般记为: 描写粒子状态的 波函数,它通常 是一个复函数。 称为de Broglie 波。此式称为自由粒子的 波函数。 (1) ψ是怎样描述粒子的状态呢? (2) ψ如何体现波粒二象性的? (3) ψ描写的是什么样的波呢? (一)波函数
电子源感 光 屏(1)两种错误的看法 1. 波由粒子组成 如水波,声波,由分子密度疏密变化而形成的一种分布。 这种看法是与实验矛盾的,它不能解释长时间单个电子衍射实验。 电子一个一个的通过小孔,但只要时间足够长,底片上增 加呈现出衍射花纹。这说明电子的波动性并不是许多电子在空间聚集在一起时才有的现象,单个电子就具有波动性。 波由粒子组成的看法夸大了粒子性的一面,而抹杀 了粒子的波动性的一面,具有片面性。 P P O Q Q O 事实上,正是由于单个电子具有波动性,才能理解氢原子 (只含一个电子!)中电子运动的稳定性以及能量量子化这样一些量子现象。
(16.4.4) (16.4.5) (图16.4a )球极坐标 薛定谔方程对氢原子的应用 (一)氢原子的薛定谔方程 前一节讨论一维运动自由粒子的薛定谔方程及 其定态解.本节要讨论氢原子中电子的运动,这与 前一节有两点不同: (1)氢原子电子作三维空间运动,因此,薛定 谔方程(16.3.3)中的波函数ψ(x,t )应换成ψ(x,y,z,t ) 或ψ(r ,t ),而22x ??应换成=??+??+??222222z y x ▽2.此▽2称为拉普拉斯算符或拉氏算符. ??????<<的薛定谔方程三维运动自由粒子)c (v 222222222z y x )m 2/(t i ??+??+??=?=?ψ?-=?ψ? (16.4.1) (2)氢原子的电子不是自由粒子,它受到氢核的库仑力,此力的作用可用它们的电势能E p 表示.因此,氢原子电子的薛定谔方程可表示如下??,见〔附录16D 〕. ??????<<的薛定谔方程氢原子电子)c (v p 2p k p 22E )m 2/p (E E E E )m 2/(t i +=+=ψ+ψ?-=?ψ? (16.4.2) *(二)氢原子的定态薛定谔方程 定态解是解决氢原子各种问题的基础.参照(16.3.4)至(16.3.6)式,可把(16.4.2)式中的波函数ψ(r ,t )分离为空间部分u (r )和时间部分f (t ),并参照(16.3.10)式写出氢原子的定态薛定谔方程,见〔附录16E 〕. ψ(r ,t )=u (r )f (t ), f (t )=C /iEt e - (16.4.3) ??????<<的定态薛定谔方程氢原子电子)c (v r 4e E 0u )E E )(/m 2(u 02p p 22πε-==-+? 氢核的质量比电子的大得多,可认为氢核不动,电子绕核转动.其电势能可表成E p =-e 2/4πε0r .此势能E p 只与电子至氢核的距离r 有关,而与方向无关,即具有球对称性,应用球极坐标较为方便.如(图16.4a ),O 表氢核,e 表电子,r 为e 至O 的距离.θ为r 与z 轴的夹角,θ称天顶角或极角.?为r 在xOy 平面的投影与x 轴的夹角.故有 x=rsin θcos ?; y=rsin θsin ?; z=rcos θ (16.4.6) 拉氏算符 2222222z y x ??+??+??=?改用球坐标(r,θ,?)表示如下:?? ()() 22222222sin r 1sin sin r 1r r r r 1???θ+θ??θθ ??θ+????=?(16.4.7) 将此▽2算符代入(16.4.4)式,便得到以球坐标表示的氢原子定态薛定谔方程. ? 郭敦仁《量子力学初步》18—19,34—35页,1978年版. ? 程守洙、江之永编,王志符、朱讠永春等修订《普通物理学》第3册177—180页,1982年修订本. ? 郭敦仁《量子力学初步》35—45页,1978年版. ? 周世勋编《量子力学》59—72页,1961年版.
§16.3 一维定态薛定谔方程的建立和求解举例 (一)一维运动自由粒子的薛定谔方程 波函数随时间和空间而变化的基本方程,是薛定谔于1926年提出的,称为薛定谔波动方程,简称波动方程或薛定谔方程,它成为量子力学的基本方程. 将(16.2.14)式分别对t 和x 求导,然后从这两式消去E 、p 、和ψ,便可得到一维运动自由粒子的薛定谔方程: ψ-=?ψ?)/iE (t 即ψ=?ψ?E t i (16.3.1) ψ=?ψ ?22)/ip (x 2 ψ=ψ ?-2222p ????? ?????<<的薛定谔方程自由粒子轴运动的沿)c x (v 方程(16.3.3)中不含有能量E 和动量p ,表明此方程是不受E 和p 的数值限制的普遍方程. 请同学们自己试一试,如果上述波函数不用复数表式(16.2.14),改用类似于(16.2.1)式的余弦函数或正弦函数表式,就不会得到合乎要求的薛定谔方程(16.3.3)式?. 这薛定谔方程不是根据直接实验结果归纳而得,也不是由经典波动理论或其他理论推导出来的,它是在物质波假设的基础上,参照经典波动方程而建立起来的.薛定谔方程在微观领域中得到广泛的应用,它推导出来的结果,都与相关实验结果符合得很好,这才是薛定谔方程正确反映微观领域客观规律的最有力的证明. (二)一维运动自由粒子的定态薛定谔方程?? 上述薛定谔方程(16.3.3)是偏微分方程,从此方程可解出波函数ψ(x ,t ).在量子力学中最重要的解,是可把波函数ψ(x,t )分离成空间部分u (x )和时间部分f (t )两函数的乘积的特解,即 〔一维运动自由粒子的定态波函数〕 ψ(x,t )=u (x )f (t )(16.3.4) 将此式代入(16.3.3)式得: 22 2dx u d )t (f )m 2/(dt df )x (u i -= 两边除以ψ=uf 得: 22 2dx u d u 1)m 2/(dt df f 1i -= 此式左边是时间t 的函数,右边是坐标x 的函数.已知t 与x 是互相独立的自变量,左右两边相等,必须是两边都等于同一常量E ,即 ? 郭敦仁《量子力学初步》16—17页,人民教育出版社1978年版. ? 郭敦仁《量子力学初步》21—22页,人民教育出版社1978年版. ? 周世勋编《量子力学》32—33页,上海科学技术出版社1961年版.
求解非线性薛定谔方程的一类数值解法 张艳敏,刘明鼎 (青岛理工大学琴岛学院,山东青岛266106) 摘要:利用非标准有限差分方法构造了求解非线性薛定谔方程的两个非标准有限差分格式。对于离散后的差分格式,把关于时间和空间的步长函数作为分母逼近导数项。对于非线性项,通过非局部的离散方法计算了这两个非标准有限差分格式的局部截断误差。数值实验结果验证了非标准有限差分格式的有效性。关键词:非线性薛定谔方程;局部截断误差;数值解法中图分类号:O241.82 文献标识码:A 文章编号:2095-7726(2019)03-0008-03 薛定谔方程是物理学中量子力学的一个重要方程,可以用于研究深水波浪理论。柱(球)非线性薛定谔方程常用于描述单色波的一维自调适、光学的自陷现象、固体中的热脉冲传播和等离子体中的Langnui 波[1–5],因此对于此类方程的研究具有非常重要的意义。 薛定谔方程有线性和非线性两种,在本文中,我们研究的是非线性薛定谔方程。非线性薛定谔方程解的解析表达式是很难得到的,因此求解此类方程最常用的就是数值解法。求非线性薛定谔方程数值解的方法主要有差分方法、配置谱方法[6]、有限元方法[7]和平均离散梯度方法[8]等。在本文中,我们利用非标准有限差分方法研究了非线性薛定谔方程的数值解,这种方法已在求解偏微分方程中得到了广泛的应用[9],其优点是对非线性项作非局部离散,对导数项作离散后用步长函数作分母,这样不仅能保持差分方程的数值解与原方程的解析解具有相同的正性,而且能保持较好的数值稳定性。 1非标准有限差分格式的构造 现在我们利用文献[10-12]给出的方法构造非线 性薛定谔方程的两种非标准有限差分格式,要考虑的非线性薛定谔方程为 (1)相应的初边值条件为 其中:为虚数单位;、、和均为连续函数;和均为正数。 为了得到非线性薛定谔方程的差分格式,需要对式 (1)进行离散。首先,需要利用网格对区域进行分割,取空间步长时间步长其 次,在网格点处,定义数值解其中,且下面将分别构造式(1)的两种非标准有限差分格式。 1.1第一种非标准有限差分格式的构造 为了构造式(1)的第一种非标准有限差分格式,我们利用R.E.Mickens 提出的构造非标准有限差分格式的原理[10]和文献[13-14]中提到的方法,并利用给定的记号,对式(1)进行离散。离散后的差分方程为 其中,和为分母函数,且,且分母是通过步长函数逼近得到的。 从式(4)可以看出,和分别取代了和分母函数的选择依据了薛定谔方程解的性质[4]。 记对式(4)进行整理,可得第36卷第3期Vol.36No.3 新乡学院学报 Journal of Xinxiang University 2019年3月Mar.2019 收稿日期:2018-12-21 基金项目:山东省高校科技计划项目(J17KB053);青岛理工大学琴岛学院教育教学研究重点项目(2018003A)作者简介:张艳敏(1981—),女,山东东营人,副教授,硕士,研究方向:偏微分方程数值分析。通信作者:刘明鼎(1982—),男,辽宁大连人,副教授,硕士,研究方向:偏微分方程数值分析。 222 (,)(,) (,)(,)(,),i u x t u x t u x t u x t g x t t x ??=++??(,0)(), u x f x =(2) 01(0,)(), (,)()u t p t u L t p t =ìí =?。 (3)0,0;x H t T <£<£i (,)g x t ()f x 0()p x 1()p x H T [0,][0,]H T ′,h H M =Δt T N =。(,)m n x t (,),n m m n u u x t =(0,1,2,,),m x mh m M ==L Δ(0,1,2,,),n t n t n N ==L ,M N ++??Z Z 。111212 2(),i n n n n n m m m m m n n n m m m u u u u u u u g D D ++---+=++(4) 1D 2D 12exp(Δ)1,D t D =-=24sin ()2 h 1D 2D Δt 2,h 11122 ,,D R D R D ==
一. 定义及重要性 薛定谔方程(Schrdinger equation )是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程,也是量子力学的一个基本假定,其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式以及对应的能量,从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程,亦是量子力学的一个基本假定,它的正确性只能靠实验来检验。 二. 表达式 三. 定态方程 ()()2 22V r E r m ηψψ+??-?=???? 所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。 其中,E 是粒子本身的能量;v(x ,y ,z)是描述势场的函数,假设不随时间变化。 可化为 薛定谔方程的解法 一. 初值解法;欧拉法,龙格库塔法
二.边值解法;差分法,打靶法,有限元法 龙格库塔法(对欧拉法的完善) 给定初值问题 有限元方法 有限元的概念早在几个世纪前就已产生并得到了应用,例如用多边形(有限个直线单元)逼近圆来求得圆的周长,但作为一种方法而被提出,则是最近的事。有限元法最初被称为矩阵近似方法,应用于航空器的结构强度计算,并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力,随着计算机技术的快速发展和普及,有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域,成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。 有限元分析的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的(较简单的)近似解,然后推导求解这个域总的满足条件,从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。不同于求解(往往是困难的)满足整个定义域边界条件的函数的Rayleigh Ritz法,有限元法将函数定义在简单几何形状(如二维问题中的三角形或任意四边形)的单元域上(分片函数),且不考虑整个定义域的复杂边界条件,这是有限元法优于其他近似方法的原因之一。
admin [非线性薛定谔方程数值解的MATLAB仿真]——利用分步快速傅里叶变换对光纤中光信号的传输方程进行数值求解
1、非线性薛定谔方程 非线性薛定谔方程(nonlinear Schrodinger equation ,NLSE)是奥地利物理学家薛定谔于1926 年提出的,应用在量子力学系统中。由于量子力学主要研究粒子的动力学运动状态,所以不能运用牛顿力学公式来表示。通常在量子力学中,研究系统的状态一般通过波函数(x ,t)来表示。而对波函数的研究主要是求解非线性薛定谔方程。本文主要研究光脉冲在光纤中传输状态下的演变。 一般情况下,光脉冲信号在光纤中传输时,同时受到光纤的色散和非线性效应的影响。通过Maxwell 方程,考虑到光纤的色散和非线性效应,可以推导出光信号在光纤中的传输方程,即非线性薛定谔方程。NLSE 是非线性偏微分方程,一般很难直接求出解析解,于是通过数值方法进行求解。具体分为两大类:(1)分布有限差分法(split-step finite differencemethod ,SSFD);(2)分步傅里叶变换法(split-step Fourier transform method ,SSFT)。一般情况,在达到相同精度,由于分步傅里叶变换法采用运算速度快的快速傅里叶变换,所以相比较有限差分法运算速度快一到两个数量级。于是本文介绍分步傅里叶变换法来对光纤中光信号的传输方程,即非线性薛定谔方程进行数值求解。并通过MATLAB 软件对结果数值仿真。 非线性薛定谔方程的基本形式为: 22||t xx iu u u u =+ 其中u 是未知的复值函数. 目前,采用分步傅立叶算法(Split step Fourier Method)求解非线性薛定谔方程的数值解应用比较多。分步傅立叶方法最早是在1937年开始应用的,这种方法己经被证明是相同精度下数值求解非线性薛定愕方程最快的方法,部分原因是它采用了快速傅立叶变换算法(FastFourier Transform Algorithm)。基于MATLAB 科学计算软件以及MATLAB 强大的符号计算功能,完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 一般情况下,光脉冲沿光纤传播时受到色散和非线性效应的共同作用,假设当传输距离 很小的时候,两者相互独立作用,那么,根据这种思想可建立如下分步傅立叶数值算法的数 学模型: 把待求解的非线性薛定谔方程写成以下形式: ??()U D N U z ?=+? (I ) (II )
薛定谔方程 一. 定义及重要性 薛定谔方程(Schrdinger equation )是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程,也是量子力学的一个基本假定,其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合建立的二阶偏微分方程,可描述微观粒子的运动,每个微观系统都有一个相应的薛定谔方程式,通过解方程可得到波函数的具体形式以及对应的能量,从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程,亦是量子力学的一个基本假定,它的正确性只能靠实验来检验。 二. 表达式 三. 定态方程 ()()2 22V r E r m ηψψ+??-?=???? 所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。 其中,E 是粒子本身的能量;v(x ,y ,z)是描述势场的函数,假设不随时间变化。
2 2 22222z y x ?? ????++=? 可化为d 0)(222=-+ψψ v E h m dx 薛定谔方程的解法 一. 初值解法;欧拉法,龙格库塔法 二. 边值解法;差分法,打靶法,有限元法 龙格库塔法(对欧拉法的完善) 给定初值问题 ).()()((3) ) ,() ,() ( ,,(2) )(),( 3112122111021h O t y t y hk y h t f k y t f k k c k c h y y y c c a y b t a y t f dt dy i i i i i i i i =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα