第二章薛定谔方程

第二章薛定谔方程

本章介绍：本章将系统介绍波动力学。波函数统计解释和态叠加原理是量子力学的两个基本假设。薛定谔方程是波动力学的核心。在一定的边界条件和初始条件下求解薛定谔方程，可以给出许多能与实验直接比较的结果。

§2.1 波函数的统计解释

§2.1.1波动—粒子两重性矛盾的分析按照德布罗意的观点，和每个粒子相联系的都有一个波。怎样理解粒子性和波动性之间的联系，这是量子力学首先遇到的根本问题。b5E2RGbCAP

2.1.1波动—粒子两重性矛盾的分析能否认为波是由粒子组成？

粒子的单缝和双缝实验表明，如减小入射粒子强度，让粒子近似的一个一个从粒子源射出，实验发现，虽然开始时底片上的感光点是无规则的，但只要时间足够长，感光点足够多，底片上仍然会出现衍射条纹。如果波是由粒子做成，那末，波的干涉、衍射必然依赖于粒子间的相互作用。这和上述实验结果相矛盾，实际上，单个粒子也具有波动性的。p1EanqFDPw

能否认为粒子是由波组成？

比如说，电子是三维空间的物质波包，波包的大小即电子的大小，波包的速度即电子的速度，但物质波包是色散的，即使原来的物质波包很小，但经过一段时间后，也会扩散到很大的空间去，或者形象地说，随着时间的推移，粒子将越来越“胖”，这与实验相矛盾

DXDiTa9E3d经典物理对自然界所形成的基本物理图像中有两类物理体系：

◆一类是实物粒子

◆另一类是相互作用场<波）经典粒子是以同时确定的坐标和动量来描述其运动状态，粒子的运动遵从经典力学规律，在运动过程中具有确定严格的轨道。粒子的能量，动量在粒子限度的空间小区域集中；当其与其它物理体系作用时，只与粒子所在处附近的粒子相互作用，并遵从能量、动量的单个交换传递过程，其经典物理过程是粒子的碰撞；“定域”是粒子运动的特征。RTCrpUDGiT经典波动则是以场量<振幅、相位等）来描述其运动状态，遵从经典波动方程，波的能量和动量周期性分布于波所传播的空间而不是集中在空间一点，即波的能量、动量是空间广延的。波与其他物质体系相互作用时，可同时与波所在广延空间内的所有物理体系相互作用，其能量可连续变化，波满足叠加原理，“非定域”是波动性运动的特性。5PCzVD7HxA◆◆在经典物理中，粒子和波各为一类宏观体系的呈现，反映着两类对象，两种物质形态，其运动特点是不相容的，即具有粒子性运动的物质不会具有波动性；反之具有波动性运动的物质不会具有粒子性。jLBHrnAILg综上所述，微观粒子既不是经典的粒子又不是经典的波，或者说它既是量子概念的粒子又是量子概念的波。其量子概念中的粒子性表示他们是具有一定的能量、动量和质量等粒子的属性，但不具有确定的运动轨道，运动规律不遵从牛顿定律；其量子概念中的波动性并不是指某个实在物理量在空间的

波动，而是指用波函数的模的平方表示在空间某处粒子被发现的概率。xHAQX74J0X◆现在被物理学家们普遍接受的波函数解释是玻恩提出的统计解释。他认为，粒子在衍射或干涉实验中所揭示的波动性质，既可以看成是大量粒子在同一实验中的统计结果，也可以认为是单个粒子在多次相同实验中显示的统计结果。LDAYtRyKfE

◆玻恩的统计解释：波函数在某一时刻在空间的强度，即其振幅绝对值的平方与在这一点找到粒子的几率成正比，和粒子联系的波是概率波Zzz6ZB2Ltk

§2.1.2波函数统计解释波函数的的特点：1.由于给出在

时刻，粒子在处出现的几率密度，因此原则上可由统计平均公式：

求出力学量的平均值。在这种意义下，波函数描述了微观粒子的运动状态，微观粒子的运动状态叫量子态。波函数应该是的单值、有界、连续函数。3.不确定性：

a.常数因子的不确定性：若为常数，则和描述同一

个物理状态。

b.相角的不确定性：由于与的模相同，因此不定。

4.可归一化：

5、容易将波函数统计解释推广到多粒子体系。

6.描述粒子微观运动的波函数与可以用其他量<如动量）为自变量。

，

薛定谔 薛定谔(Schroding,1897-1961>奥地利人,因发现原子理论的有效的新形式一波动力学与狄拉克(Dirac,1902-1984>因创立相对论性的波动方程一狄拉克方程,共同分享了1933年度诺贝尔物理学奖dvzfvkwMI1

玻恩M.玻恩，(Max Born 1882～1970>德国理论物理学家，量子力学的奠基人之一。主要成就是创立矩阵力学和对波函数作出统计解释。1954年因波函数的统计解释荣获诺贝尔物理学奖。rqyn14ZNXI §2.2态叠加原理

态叠加原理是量子力学中一个很重要的原理，这一节先作一些初步介绍，随着学习量子力学内容的不断深入，会不断加深对态叠加原理的理解。EmxvxOtOco

态叠加原理：

如果是体系可能的状态，则它们的线性叠加所得出的波函数

也是体系的一个可能状态；当体系处于态时，出现的概率是

，可以是有限的，也可以是无限的。

几点讨论：

I.测量力学量得出的是一些可能值但这些可能值的相对概率，或者说每个可能态的相对权重，是完全确定的。

II.态叠加原理中所谓的叠加，是波函数的叠加，或者说是概率幅的叠加，而不是概率的叠加。因而它必然会出现干涉、衍射等现象。SixE2yXPq5

III.在量子力学中，对于概率波而言，波的干涉是描述粒子运动状态的概率波本身的干涉，而不是粒子之间的干涉。6ewMyirQFL

IV.一般来说，依赖于时间，是的函数，因此态叠加原理不仅对某一时刻成立，而且随时间的变化，态叠加原理仍然成立。

kavU42VRUs

§2.3薛定谔方程

经典力学中，体系运动状态随时间的变化遵循牛顿力学。和经典力学类似，我们也应建立一个决定波函数随时间变化规律的方程式。从物理上，这个方程式必须满足下述条件：y6v3ALoS89

I.由于波函数满足态叠加原理，而态叠加原理对任何时间都成立，因此描述波函数随时间变化的方程应该是线性方程。M2ub6vSTnPII.方程的系数仅含有质量、电荷等内禀量，不应含有和个别粒子运动状态特定性质有关的量，如动量。

III.因为波函数的自变量是坐标和时间，因此它必然是关于坐标和时间的偏微分方程。

IV.由于经典力学是量子力学的极限情况，因此这个方程必须满足对应原理，当取经典极限时，它能过渡到牛顿方程。0YujCfmUCw

V.对于自由粒子，这个方程的解应该是单色平面波的波函数。

方程的建立

对平面波式

分别对坐标和时间求微商后得：，

由上两式可以看出能量与动量作用在波函数上的结果与算符及作用在波函数上的结果相同，即存在对应关系：

eUts8ZQVRd?1926年，薛定谔推广上述规则到一般情况，建立了描述波函数演化规律的薛定谔方程，得到薛定谔方程：

薛定谔方程式量子力学的基本假设之一，但必须指出，我们并未建立薛定谔方程，因为只知道微分方程的解是不足以建立微分方程的。sQsAEJkW5TB.以上对应关系式<2.3.3）式，只是在直角坐标系中的对应关系，在其他坐标系中不一定成立。

下面我们讨论一下定态情况：若势能不显含时间，则薛定谔方程可用分离变量法求解，此时可令：将上式代入薛

定谔方程并用遍除等式两边，可得：，

此即定态薛定谔方程。

方程<2.3.5）的解可直接给出为代入<2.3.4）并将吸

收入中去，并有归一化条件来确定，有，按照德布罗意关系，就是体系出于这个波函数所描写的状态时的能量。由此可见，体系出于上述波函数所描述的状态时，能量具有确定值，这种状态称为定态。波函数称为定态波函数。GMsIasNXkA以

表示体系的能量算符的第个本征值，是与相应的波函数，则

体系的第个定态波函数是

含时的薛定谔方程的一般解，可以写成这些定态波函数的线性叠加：

§2.4概率流密度与概率流守恒定律?本节我们将进一步讨论粒子在一定区域内出现的几率将怎样随时间变化。

?设描述粒子状态的波函数是，在时刻、在点周围单位体积内粒子出现的几率是

?几率密度随时间的变化率为由薛定谔方程及其

共轭：，可得：

令：

称为概率流密度，由<2.4.1）式得：<2.4.2）式就是概率流守恒定律对上式两边同时对任意空间体积积分

这是概率流守恒定律的积分表示。此式表明，在空间某体积内发现粒子的概率在单位时间内的增量，必定等于在同一时间内通过

的边界流入体积的概率。TIrRGchYzg

若以粒子的质量乘和，则有：是在时刻在

点的质量密度。是质量流密度，满足：即量子力学中的质量守恒定律。

B.同样，以粒子电荷乘和后，得到是电荷密度，

是电流密度，

方程

是量子力学中的电荷守恒定律。

§2.5一维方势阱本节将以一维定态为例，求解已知势场的定态薛定谔方程。了解怎样确定定态的能量E，从而看出能量量子化是薛定谔方程的自然结果。7EqZcWLZNX

§2.5.1一维无限深方势阱已知粒子所处的势场为：粒子在势阱势能为零，在阱外势能为无穷大，在阱壁上受极大的斥力。称为一维无限深方势阱。其定态薛定谔方程为：

，当时，根据波

函数的连续性和有限性条件得：，

令：

则薛定谔方程可简写为：，

它的解是：，

利用边界条件及，得

，

解为：，，，为奇数）lzq7IGf02E，，，为偶数） zvpgeqJ1hk带入<2.5.1）得体系的能级：

显然，一维无限深方势阱的能谱是分立谱，这个分离的能谱就是量子化了的能级。

由图可以看出，在不同能级上粒子出现的概率密度是不同的。在基态，粒子出现的概率在阱区中部为最大，而越靠近阱壁概率越小，阱壁上概率为零。在激发态，粒子在阱内出现的概率是起伏变化的，随着量子数的增大，起伏变化越来频繁。NrpoJac3v1而在经典物理中，粒子在阱内各处出现的概率是相等的。 由图可以推断，只有当量子数很大时，粒子在阱内各处的概率才趋于均匀。粒子的最低能量状态称为基态，就是的状态，基态能量为

此本征值能量称为零点能，是束缚在无限深方势阱内粒子所具有的最低能量。归一化以后的波函数为：，

。我们把粒子只能束缚在空间的有限区域，在无穷远处

波函数为零的状态称为束缚态。

§2.5.2一维有限深方势阱

求解势场为的薛定谔方程。讨论的情况：在区，相应的薛定谔方程是，

在时，有界的解是：，

在区，薛定谔方程是：，

其解为

1.在区，取，解取有偶宇称的情况

利用处波函数对数微商的连续条件都可得引入

可将<2.5.3）是改写为

另外，又<2.5.1）和<2.5.2）有可得

联立<2.5.5）--<2.5.6）式，解出，再由<2.5.4）可给出能谱。

2.在区，取，解取有奇宇称的情况

同样，利用波函数对数微商在连续条件得：

同样，联立<2.5.6）--<2.5.7）式，解出，再由<2.5.4）可给出能谱。

<2.5.5）--<2.5.7）都是超越方程，可用图解法求出能谱。

在平面中分别就<2.5.5）与<2.5.6）式作相应的曲线，曲线的交

点表示具有偶宇称是相应的能谱。如上图。1nowfTG4KI由以上图可见，对于偶宇称态，由于曲线经过原点，因此无论多么小，两条曲线总有交点，这意味着至少有一个束缚态，且相应的宇称为偶。fjnFLDa5Zo

同样，作<2.5.6）和<2.5.7）式相应曲线，他们的交点表示波函数其宇称时相应的能谱。所得结果见上图。tfnNhnE6e5

由以上图可见，对于奇宇称态，当且仅当时，即当时，曲线才有交点，才出现奇宇称态解。显然，一维无限深势阱的结果可作为一维方势阱的特例得出。当时，可得这正是阱宽为的一维无限深势阱的能谱公式。

§2.6 一维方势垒前面讨论了束缚态，这一节我们讨论散射态首先讨论一维方势垒问题。

势垒为设能量为的粒子从势垒的左方向右方运动，下面分别就与来讨论。

1.的情形

此时，满足的薛定谔方程为

为方便起见，令方程可改为：

其解分别为

利用在和处波函数连续性和波函数微商连续性条件

可得出与关系

由概率流密度公式可得入射波的概率流密度为透射波的概率流密度为：反射波的概率流密度为：反射系数为：透射系数为：

由上两式可见，一般情况下，

透射系数，反射系数，而这之和为1。这表明，在量子力学中，即是粒子的能量大于势垒高度，仍有部分被反射回来。这正是微观粒子具有波动性的体现。HbmVN777sL由第二式可见，一般情况下透射系数，当的特定情况下，其透射系数，这种情形下的透射现象叫做共振透射。

2. 的情形

此时，为虚数。但若令，则系数关系变为

反射系数和透射系数为：

由此可见，反射系数和透射系数，且二者之和等于1，这表明，在量子力学中，即使粒子的能量小于势垒的高度，粒子仍有一部分透射过去。V7l4jRB8Hs 这种粒子在其能量小于势垒高度时，仍然会有部分粒子穿过势垒的现象叫

隧道效应，又叫隧穿效应隧道效应的应用：1.扫描隧道显微镜

扫描隧道显微镜可以显示表面原子台阶和原子排布的表面三维图案。

在表面物理、材料科学和生命科学等诸多领域中，扫描隧道显微镜都能提供十分有价值的信息。

2.隧道二极管是一种利用隧道效应的半导体器件，也是隧道效应的重要应用之一。

由于隧道效应而使其伏安特性曲线出现负阳区，因而隧道二级管具有高频、低噪声的特点。

隧道二级管是低频放大器、低频噪声振荡器和超高速开关电路中的重要器件。

§2.7 一维谐振子本节我们来讨论一维谐振子问题。

一维谐振子的哈密顿量为：满足的定态定谔方程为：为方便求解，引入系数：

则方程可改写为：

这是一个变系数的二阶常微分方程，当很大时，，上式中的

可略去。从而，得到上式的渐进方程

其解就是原方程的解，又由于波函数在时的有限性

条件，得

为了求出在整个空间都合适的解，可以将系数视为的某一特定函

数，假设方程的解为

在有限时应该有限，在时它的行为也必须保证波函数

有限。代回薛定谔方程，得到待定系数满足的方程

对作泰勒展开

可由

得

当时，的渐进行为是

与的渐进行为相同。若为无穷级数时，在时将趋向无穷大。为了在时，波函数仍有限，必须断为多项

式。因为如果是多项式，当时，它趋于无穷的行为永远

比趋于零慢，从而保证了在是有限。由<2.7.2）可知方程(2.7.1>有解的条件为83lcPA59W9此时，

有

这是厄M方程，其解为厄M多项式。厄M多项式有三种重要表示：1.级数表示：

式中为偶数时。为奇数。

2．积分表示： 3.微分表示：

厄M多项式具有如下性质：

递推关系微分性质正交归一性：完备性：

式中的展开系数为：

由式<2.7.1）即可得能量本征值为：

叫振动量子数。从而其波函数为：式中归一

化常数为：由<2.7.2）可见，一维谐振子的能量也是量子化的，并且能量间隔相等，为。一维谐振子基态能量：

叫零点能。mZkklkzaaP 经典与量子的比较

按经典力学的结论，一维谐振子的能量如图谐振子只能处于的范围内,

的区域则是经典禁区。

而在量子力学中，由于隧道效应，粒子可以到达经典进取，也就是说在所谓经典禁区内发现粒子的概率不为零。按经典力学的规律，在处振子的速度最大停留时间最短，在处振子的速度为零停留时间最长。将这一规律应用于微观粒子，自然会得到在

处粒子出现的概率最小，而在处粒子出现的概率最大。而实际情况如何呢？由时的波函数及概率密度的图：AVktR43bpw 可以看出，在量子数较小的时候，粒子位置的概率密度的分布与经典结论明显不同。

可以推断，随着量子数的增大，概率密度的平均值将越来越接近经典结论。

§2.8一维周期场

设空间周期为，考虑到势场的周期性条件：

,为任意整数。则在晶格周期势场中运动粒子的薛定

谔方程为

显然，都满足上述薛定谔方程，并且具有同一本征值，从而可以得到：

其中为一常数。

由波函数的有界性，当时，若为实数，必使。所以，必为复相因子，令则：

所以，粒子在空间呈现的几率也是周期性的。在某一个周期内，定态薛定谔方程为

下面就和两种情况分别讨论：1.的情况

令

则方程的解为：

同理下一个周期中的解为：

由在处的连续性条件

在处的连续性条件

稍加整理，有：

具有非零解的条件是其

系数行列式为零

展开并整理，再除以得

所以

于是

上式为一超越方程，为简单起见，只讨论

只讨论的极限情形，此时，有

则：同时，此时有：，

所以，=由于

则

定态薛定谔方程讲义

定态薛定谔方程 一、定态Schr?dinger 方程 2 2(,)[()](,)2i r t V r r t t m ψψ?=-?+? (1) 在一般情况下，从初始状态ψ(r,0)求 ψ(r,t)是不容易的。以下，我们考虑一个很重要的特殊情形——假设势场V 不显含时间 t （在经典力学中，在这种势场中运动的粒子，其机械能守恒），此时薛定谔方程（1）可以用分离变量数法求其特解。 ()V r 与t 无关时，可以分离变量 令(,)()()r t r f t ψψ= 代入(1)式 2 2()1[()]()()()2i df t V r r f t dt r m ψψ=-?+ E = 其中E 是即不依赖于t ，也不依赖于r 的常量，这样 ()()df t i Ef t dt = (2) 2 2[()]()()2V r r E r ψψμ -?+= (3) ——定态薛定谔方程 由(2)解得 Et i ce t f -=)( 其中c 为任意常数。把常数c 放到()E r ψ 里面去，则 (,)()i Et E r t r e ψψ-= (4) 这个波函数与时间的关系是正弦式的，其角频率是ω=Ε/?按照德布罗意关系E=h ν=?ω，E 就是该体系处于这个波函数所描写状态时的能量。由此可见，当体系处于（4）式所描写状态时，能量具有确定值E ，所以这种状态称为定态，波函数ψ(r,t)称为定态波函数。 定态有两个含义：1、(,)()i Et E r t r e ψψ-= ；2、E 具有确定值；（判断是否为定态的依 据） 空间波函数()E r ψ 可由方程 2 2[()]()()2E E V r r E r m ψψ-?+= 和具体问题()E r ψ 应满足的边界条件得出。方程(3)称为定态Schr?dinger 方程，()E r ψ 也可

最新薛定谔方程及其解法

关于薛定谔方程 一．定义及重要性 薛定谔方程（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) 其中，是质量，是位置，是相依于时间的波函数，是约化普朗克常数，是位势。类似地，在三维空间里，一个单独粒子运动于位势中的含时薛定谔方程为 。(2) 假若，系统有个粒子，则波函数是定义于 -位形空间，所有可能的粒子位置空间。用方程表达， 。 其中，波函数的第个参数是第个粒子的位置。所以，第个粒子的位置是。 不含时薛定谔方程 不含时薛定谔方程不相依于时间，又称为本征能量薛定谔方程，或定态薛定谔方程。顾名思义，本征能量薛定谔方程，可以用来计算粒子的本征能量与其它相关的量子性质。 应用分离变量法，猜想的函数形式为 ； 其中，是分离常数，是对应于的函数．稍回儿，我们会察觉就是能量． 代入这猜想解，经过一番运算，含时薛定谔方程 (1) 会变为不含时薛定谔方程： 。 类似地，方程 (2) 变为

。 历史背景与发展 爱因斯坦诠释普朗克的量子为光子，光波的粒子；也就是说，光波具有粒子的性质，一种很奇奥的波粒二象性。他建议光子的能量与频率成正比。在相对论里，能量与动量之间的关系跟频率与波数之间的关系相同，所以，连带地，光子的动量与波数成正比。 1924年，路易·德布罗意提出一个惊人的假设，每一种粒子都具有波粒二象性。电子也有这种性质。电子是一种波动，是电子波。电子的能量与动量决定了它的物质波的频率与波数。1927年，克林顿·戴维和雷斯特·革末将缓慢移动的电子射击于镍晶体标靶。然后，测量反射的强度，侦测结果与X射线根据布拉格定律 (Bragg's law) 计算的衍射图案相同。戴维森-革末实验彻底的证明了德布罗意假说。 薛定谔夜以继日地思考这些先进理论，既然粒子具有波粒二象性，应该会有一个反应这特性的波动方程，能够正确地描述粒子的量子行为。于是，薛定谔试着寻找一个波动方程。哈密顿先前的研究引导著薛定谔的思路，在牛顿力学与光学之间，有一种类比，隐蔽地暗藏于一个察觉里。这察觉就是，在零波长极限，实际光学系统趋向几何光学系统；也就是说，光射线的轨道会变成明确的路径，遵守最小作用量原理。哈密顿相信，在零波长极限，波传播会变为明确的运动。可是，他并没有设计出一个方程来描述这波行为。这也是薛定谔所成就的。他很清楚，经典力学的哈密顿原理，广为学术界所知地，对应于光学的费马原理。借着哈密顿-雅可比方程，他成功地创建了薛定谔方程。薛定谔用自己设计的方程来计算氢原子的谱线，得到了与用玻尔模型计算出的能级相同的答案。 但是，薛定谔对这结果并不满足，因为，索末菲似乎已经正确地计算出氢原子光谱线精细结构常数的相对论性的修正。薛定谔试着用相对论的能量动量关系式，来寻找一个相对论性方程（现今称为克莱因-高登方程），可以描述电子在库仑位势的量子行为。薛定谔计算出这方程的定态波函数。可是，相对论性的修正与索末菲的公式有分歧。虽然如此，他认为先前非相对论性的部分，仍旧含有足够的新结果。因此，决定暂时不发表相对论性的修正，只把他的波动方程与氢原子光谱分析结果，写为一篇论文。1926年，正式发表于物理学界[2]。从此，给予了量子力学一个新的发展平台。 薛定谔方程漂亮地解释了的行为，但并没有解释的意义。薛定谔曾尝试解释代表电荷的密度，但却失败了。1926年，就在薛定谔第四篇的论文发表之后几天，马克斯·玻恩提出概率幅的概念，成功地解释了的物理意义[3]。可是，薛定谔本人一直不承认这种统计或概率的表示方法，和它所伴随的非连续性波函数坍缩。就像爱因斯坦的认为量子力学是基本为确定性理论的统计近似，薛定谔永远无法接受哥本哈根诠释。在他有生最后一年，他写给马克斯·玻恩的一封信，薛定谔清楚地表明了这看法。 含时薛定谔方程导引

实验三 定态薛定谔方程的矩阵解法

实验三 定态薛定谔方程的矩阵解法 一．实验目的 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．实验结果 程序设计：

薛定谔方程

第一章 薛定谔方程 §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 ( -?=ψ

第二章 薛定谔方程

第二章 薛定谔方程 本章介绍：本章将系统介绍波动力学。波函数统计解释和态叠加原理是量子力学的两个基本假设。薛定谔方程是波动力学的核心。在一定的边界条件和初始条件下求解薛定谔方程，可以给出许多能与实验直接比较的结果。 §2.1 波函数的统计解释 §2.1.1 波动—粒子两重性矛盾的分析按照德布罗意的观点，和每个粒子相联系的都有一个波。怎样理解粒子性和波动性之间的联系，这是量子力学首先遇到的根本问题。 2.1.1 波动—粒子两重性矛盾的分析能否认为波是由粒子组成？ 粒子的单缝和双缝实验表明，如减小入射粒子强度，让粒子近似的一个一个从粒子源射出，实验发现，虽然开始时底片上的感光点是无规则的，但只要时间足够长，感光点足够多，底片上仍然会出现衍射条纹。如果波是由粒子做成，那末，波的干涉、衍射必然依赖于粒子间的相互作用。这和上述实验结果相矛盾，实际上，单个粒子也具有波动性的。 能否认为粒子是由波组成？ 比如说，电子是三维空间的物质波包，波包的大小即电子的大小，波包的速度即电子的速度，但物质波包是色散的，即使原来的物质波包很小，但经过一段时间后，也会扩散到很大的空间去，或者形象地说，随着时间的推移，粒子将越来越“胖”，这与实验相矛盾 经典物理对自然界所形成的基本物理图像中有两类物理体系： ◆一类是实物粒子 ◆另一类是相互作用场（波）经典粒子是以同时确定的坐标和动量来描述其运动状态，粒子的运动遵从经典力学规律，在运动过程中具有确定严格的轨道。粒子的能量，动量在粒子限度的空间小区域集中；当其与其它物理体系作用时，只与粒子所在处附近的粒子相互作用，并遵从能量、动量的单个交换传递过程，其经典物理过程是粒子的碰撞；“定域”是粒子运动的特征。经典波动则是以场量（振幅、相位等）来描述其运动状态，遵从经典波动方程，波的能量和动量周期性分布于波所传播的空间而不是集中在空间一点，即波的能量、动量是空间广延的。波与其他物质体系相互作用时，可同时与波所在广延空间内的所有物理体系相互作用，其能量可连续变化，波满足叠加原理，“非定域”是波动性运动的特性。◆◆在经典物理中，粒子和波各为一类宏观体系的呈现，反映着两类对象，两种物质形态，其运动特点是不相容的，即具有粒子性运动的物质不会具有波动性；反之具有波动性运动的物质不会具有粒子性。综上所述，微观粒子既不是经典的粒子又不是经典的波，或者说它既是量子概念的粒子又是量子概念的波。其量子概念中的粒子性表示他们是具有一定的能量、动量和质量等粒子的属性，但不具有确定的运动轨道，运动规律不遵从牛顿定律；其量子概念中的波动性并不是指某个实在物理量在空间的波动，而是指用波函数的模的平方表示在空间某处粒子被发现的概率。◆现在被物理学家们普遍接受的波函数解释是玻恩提出的统计解释。他认为，粒子在衍射或干涉实验中所揭示的波动性质，既可以看成是大量粒子在同一实验中的统计结果，也可以认为是单个粒子在多次相同实验中显示的统计结果。 ◆玻恩的统计解释：波函数在某一时刻在空间的强度，即其振幅绝对值的平方与在这一点找到粒子的几率成正比，和粒子联系的波是概率波 §2.1.2 波函数统计解释 波函数的的特点：1.由于 2 |),(|t r ψ给出在 t 时刻，粒子在 r 处出现的几率密度，因此原 则上可由统计平均公式：? ?>= <)(r f 。在这种意义下，波函数),(t r ψ描述了微观粒子的运

非线性薛定谔方程数值解的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年开始应用的，这种方法己经被证明是相同精度下数值求解非线性薛定愕方程最快的方法，部分原因是它采用了快速傅立叶变换算法(Fast Fourier Transform Algorithm)。基于MATLAB 科学计算软件以及MATLAB 强大的符号计算功能，完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 一般情况下，光脉冲沿光纤传播时受到色散和非线性效应的共同作用，假设当传输距离 很小的时候，两者相互独立作用，那么，根据这种思想可建立如下分步傅立叶数值算法的数 学模型： 把待求解的非线性薛定谔方程写成以下形式: ??()U D N U z ?=+? （I ） （II ）

§4.3薛定谔方程

§4.3 薛定谔方程 在这一节，我们讨论态随时间变化的规律问题。大家知道，在经典力学中，当质点的初始状态为已知时，由其运动方程就可以知道以后任一时刻的运动状态。在量子力学中的情况也是这样的，即当粒子在初始时刻的态为已知时，在以后任一时刻的态也要由一个相应的方程来决定。所不同的是：在经典力学中，质点的状态用质点的坐标和速度描写，质点的运动方程就是我们所熟知的牛顿运动方程。而在量子力学中，微观粒子的状态则用波函数来描写，决定粒子状态变化的方程不再是牛顿运动方程，而是下面我们要建立的薛定谔方程。从物理上，这个方程式 必须满足下述条件： 一、在非相对论条件下，薛定谔方程应该满足的条件 1、在粒子的速度v c 时，质量为m 的粒子的总能量为：22p E U m =+ 2、方程是线性的 由于波函数满足态叠加原理，而态叠加原理对任何时间都成立，因此描述波函数随时间变化的方程应该是线性方程。即如果1ψ和2ψ是方程的解，那么它们的线性迭加 1122c c ψ+ψ也是方程的解。 3、方程的系数仅含有质量、电荷等内禀量，不应含有和个别粒子运动状态特定性质有关的 量，如动量、能量等。 4. 方程应当是波函数 (,)r t ψ 对时间的一阶微分方程 因为我们所要建立的是波函数(,)r t ψ 随时间变化的运动方程，而波函数完全描述态，因此方程必须波函数 ),(t r ψ对时间的一阶微分方程。也就是说方程必然包 含(,)r t t ?ψ? ，但方程不包含22 (,)r t t ?ψ? ，否则需要利用两个初始条件(,0)r ψ 和0(,)|t r t t =?ψ? 才能确定),(t r ψ，这就意味着体系的初始状态不能由波函数(,0)r ψ 完全描述，违反了波函数完全描述态体系运动状态的基本假设。 二、自由粒子波函数所满足的微分方程 下面，就以自由粒子为例，来建立满足上述条件的运动方程。自由粒子的波函数就是德布罗意平面波函数 ()()·,i p r Et r t Ae -ψ= （1）

薛定谔方程

薛定谔方程(Schr?dinger equation)是一个由奥地利物理学家薛定谔在1926年描述量子力学中波函数的运动方程[1]，被认为是量子力学的奠基理论之一。 薛定谔方程主要分为含时薛定谔方程与不含时薛定谔方程。含时薛定谔方程相依于时间，专门用来计算一个量子系统的波函数，怎样随着时间演变。不含时薛定谔方程不相依于时间，可以计算一个定态量子系统，对应于某本征能量的本征波函数。波函数又可以用来计算，在量子系统里，某个事件发生的概率幅。而概率幅的绝对值的平方，就是事件发生的概率密度。 薛定谔方程的解答，清楚地描述量子系统里，量子尺寸粒子的统计性量子行为。量子尺寸的粒子包括基本粒子，像电子、质子、正电子、等等，与一组相同或不相同的粒子，像原子核。 薛定谔方程可以转换为海森堡的矩阵力学，或费曼的路径积分表述 (path integral formulation) 。薛定谔方程是个非相对论性的方程，不能够用于相对论性理论。海森堡表述比较没有这么严重的问题；而费曼的路径积分表述则完全没有这方面的问题。 目录 [隐藏] ? 1 含时薛定谔方程 ? 2 不含时薛定谔方程 ? 3 历史背景与发展 ? 4 含时薛定谔方程导引 o 4.1 启发式导引 ? 4.1.1 假设 ? 4.1.2 波函数以复值平面波来表达波函数 o 4.2 薛定谔的导引 ? 5 特性 o 5.1 线性方程 ? 5.1.1 证明 o 5.2 实值的本征态 o 5.3 幺正性 ? 5.3.1 证明 o 5.4 完备基底 ? 6 相对论性薛定谔方程 ?7 解析方法 ?8 实例 o8.1 自由粒子 o8.2 一维谐振子 o8.3 球对称位势 ?8.3.1 角部分解答

非线性薛定谔方程的孤子解和怪波解

非线性薛定谔方程的孤子解和怪波解 摘要：光纤中光波的传输模型一直是当前研究的热点理论模型之一，从非线性薛定谔方程到金格堡-朗道方程，都试图对其进行更好的阐释，其次对于非线性动力学系统中，非线性薛定谔方程的解有呈现出非常多有趣的特征，对于其中特定解的研究能够让我们了解脉冲演化的本质，所以本文主要从孤子解的传输入手，并且简单介绍了怪波解的解形式。 薛定谔方程又称薛定谔波动方程，是量子力学的一个基本方程，同时又是量子力学的基本假设之一，由奥地利物理学家薛定谔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 =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα

大学物理-一维定态薛定谔方程的应用

一维定态薛定谔方程 的应用 授课人： 物理科学与技术学院

势 阱 日常生活中的各种井（阱） 物理学中研究微观粒子运动状态时常用的模型，因其势能函数曲线的形状如同井而得名 水井 窨井 陷阱 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 =???，

薛定谔方程与它的基本意义

薛定谔方程 维基百科，自由的百科全书 跳转到：导航, 搜索 汉漢▼ 量子力学 不确定性原理 入门·数学表述显示▼背景 经典力学·旧量子论·干涉 哈密顿量·狄拉克符号 显示▼基本概念 量子态·波函数·态矢量 态叠加原理·波粒二象性 量子测量·不确定性原理 泡利不相容原理·量子缠结 量子脱散·量子隧穿效应 埃伦费斯特定理 显示▼实验 双缝实验·薛定谔的猫 戴维孙-革末实验 施特恩-格拉赫实验 贝尔不等式实验 波普尔实验·量子擦除器 显示▼构想

薛定谔绘景·海森堡绘景 相互作用绘景·矩阵力学 求和的历史 显示▼方程 薛定谔方程·泡利方程 克莱因-高登方程 狄拉克方程 显示▼量子力学诠释 哥本哈根诠释·Ensemble 隐变量·交易诠释 多世界诠释·一致性历史 系综诠释·量子逻辑 显示▼进阶理论 量子场论·量子引力 万有理论 显示▼科学家 普朗克、玻尔、薛定谔、海森堡 泡利、德布罗意、埃伦费斯特、玻姆 玻恩、爱因斯坦、冯?诺伊曼 费曼、狄拉克、维恩、埃弗里特 索末菲、其他 本模板：查看? 讨论? 编辑? 历史 薛定谔方程是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程[1]，也是量子力学的一个基本假定，其正确性只能靠实验来检验。就好像牛顿定律在经典力学的地位，薛定谔方程在量子力学里占有中心的地位。 薛定谔方程主要分为含时薛定谔方程与不含时薛定谔方程。含时薛定谔方程相依于时间，专门用来计算一个量子系统的波函数，怎样随着时间演变。不含时薛定谔方程不相依于时间，可以计算一个定态量子系统，对应于某本征能量的本征波函数。波函数又可以用来计算，在量子系统里，某个事件发生的几率幅。而几率幅的绝对值的平方，就是事件发生的几率密度。薛定谔方程的解答，清楚地描述量子系统里，量子尺寸粒子的统计性量子行为。量子尺寸的

§16.3 一维定态薛定谔方程的建立和求解举例

§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年版.

薛定谔方程及其解法

一． 定义及重要性 薛定谔方程（Schrdinger equation ）是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程，也是量子力学的一个基本假定，其正确性只能靠实验来检验。是将物质波的概念和波动方程相结合建立的二阶偏微分方程，可描述微观粒子的运动，每个微观系统都有一个相应的薛定谔方程式，通过解方程可得到波函数的具体形式以及对应的能量，从而了解微观系统的性质。 薛定谔方程是量子力学最基本的方程，亦是量子力学的一个基本假定，它的正确性只能靠实验来检验。 二． 表达式 三． 定态方程 ()()2 22V r E r m ηψψ+??-?=???? 所谓势场，就是粒子在其中会有势能的场，比如电场就是一个带电粒子的势场；所谓定态，就是假设波函数不随时间变化。 其中，E 是粒子本身的能量；v(x ，y ，z)是描述势场的函数，假设不随时间变化。 可化为 薛定谔方程的解法 一． 初值解法；欧拉法，龙格库塔法

二．边值解法；差分法，打靶法，有限元法 龙格库塔法（对欧拉法的完善） 给定初值问题 有限元方法 有限元的概念早在几个世纪前就已产生并得到了应用，例如用多边形（有限个直线单元）逼近圆来求得圆的周长，但作为一种方法而被提出，则是最近的事。有限元法最初被称为矩阵近似方法，应用于航空器的结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。 有限元分析的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件，从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。不同于求解（往往是困难的）满足整个定义域边界条件的函数的Rayleigh Ritz法，有限元法将函数定义在简单几何形状（如二维问题中的三角形或任意四边形）的单元域上（分片函数），且不考虑整个定义域的复杂边界条件，这是有限元法优于其他近似方法的原因之一。

非线性薛定谔方程数值解的MATLAB仿真

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 ）

第二章薛定谔方程习题

第二章 薛定谔方程 习题 (课本44页) 证明在定态中，概率流密度与时间无关。 证明：当一个系统处于定态时，其波函数),(t r ?可以写作， ?? ? ??-=Et i r t r ex p )(),(φ? 于是便有， ?? ? ??=Et i r t r ex p )(),(**φ? 根据概率流密度的定义式有， ????? ????? ??-??? ??-???? ????? ??-= ??????????????? ??-???? ??-????????? ?????? ??-=?-?≡????????ψψψψt iE t iE t iE t iE m i t iE t iE t iE t iE m i m i J exp exp exp exp 2exp exp exp exp 2) (2* ***** 即有，)(2)(2****φφφφ?????-?=?-?=m i m i J 显然，在定态中概率流密度与时间无关。从某种意义上说明上述波函数称为定态波函数是名副其实的。 由下列两定态波函数计算概率流密度：⑴ )exp(11ikr r = ?，⑵ )exp(1 2ikr r -=?。 从所得结果说明1?表示向外传播的球面波，2?表示向内(即向原点)传播的球面波。 解：在解本题之前，首先给出一个函数f 的梯度在球坐标系下的表达式，即 ? θθ? θ??+??+??=?f r e f r e r f e f r sin 1?1?? ⑴ 首先求解函数1?的概率流密度 r ikr ikr r ikr ikr ikr r ikr e mr k r ike r e e r e r ike r e e r e m i r ikr r ikr r ikr r ikr m i m i J ???2)exp()exp()exp()exp(2) (22221*1*111 =?????????? ??+--???? ??-+-=?? ? ???---?=?-?=---???? 可见，概率流密度1J 与r 同号，这便意味着1J 的指向是向外的，即1?表示向外传播的球面波。

第二章薛定谔方程

第二章薛定谔方程 本章介绍：本章将系统介绍波动力学。波函数统计解释和态叠加原理是量子力学的两个基本假设。薛定谔方程是波动力学的核心。在一定的边界条件和初始条件下求解薛定谔方程，可以给出许多能与实验直接比较的结果。 §2.1 波函数的统计解释 §2.1.1波动—粒子两重性矛盾的分析按照德布罗意的观点，和每个粒子相联系的都有一个波。怎样理解粒子性和波动性之间的联系，这是量子力学首先遇到的根本问题。b5E2RGbCAP 2.1.1波动—粒子两重性矛盾的分析能否认为波是由粒子组成？ 粒子的单缝和双缝实验表明，如减小入射粒子强度，让粒子近似的一个一个从粒子源射出，实验发现，虽然开始时底片上的感光点是无规则的，但只要时间足够长，感光点足够多，底片上仍然会出现衍射条纹。如果波是由粒子做成，那末，波的干涉、衍射必然依赖于粒子间的相互作用。这和上述实验结果相矛盾，实际上，单个粒子也具有波动性的。p1EanqFDPw 能否认为粒子是由波组成？ 比如说，电子是三维空间的物质波包，波包的大小即电子的大小，波包的速度即电子的速度，但物质波包是色散的，即使原来的物质波包很小，但经过一段时间后，也会扩散到很大的空间去，或者形象地说，随着时间的推移，粒子将越来越“胖”，这与实验相矛盾

DXDiTa9E3d经典物理对自然界所形成的基本物理图像中有两类物理体系： ◆一类是实物粒子 ◆另一类是相互作用场<波）经典粒子是以同时确定的坐标和动量来描述其运动状态，粒子的运动遵从经典力学规律，在运动过程中具有确定严格的轨道。粒子的能量，动量在粒子限度的空间小区域集中；当其与其它物理体系作用时，只与粒子所在处附近的粒子相互作用，并遵从能量、动量的单个交换传递过程，其经典物理过程是粒子的碰撞；“定域”是粒子运动的特征。RTCrpUDGiT经典波动则是以场量<振幅、相位等）来描述其运动状态，遵从经典波动方程，波的能量和动量周期性分布于波所传播的空间而不是集中在空间一点，即波的能量、动量是空间广延的。波与其他物质体系相互作用时，可同时与波所在广延空间内的所有物理体系相互作用，其能量可连续变化，波满足叠加原理，“非定域”是波动性运动的特性。5PCzVD7HxA◆◆在经典物理中，粒子和波各为一类宏观体系的呈现，反映着两类对象，两种物质形态，其运动特点是不相容的，即具有粒子性运动的物质不会具有波动性；反之具有波动性运动的物质不会具有粒子性。jLBHrnAILg综上所述，微观粒子既不是经典的粒子又不是经典的波，或者说它既是量子概念的粒子又是量子概念的波。其量子概念中的粒子性表示他们是具有一定的能量、动量和质量等粒子的属性，但不具有确定的运动轨道，运动规律不遵从牛顿定律；其量子概念中的波动性并不是指某个实在物理量在空间的

波函数和薛定谔方程-力学量算符

波函数和薛定谔方程-力学量算符1．一维运动的粒子处在 的状态，其中，求： （1）粒子动量的几率分布函数； （2）粒子动量的平均值。 [解]首先将归一化，求归一化系数A。 （1）动量的几率分布函数是 注意到中的时间只起参数作用，对几率分布无影响，因此可有 令 代入上式得

（2） 动量p的平均值的结果从物理上看是显然的，因为对本题说来，粒子动量是和是的几率是相同的。讨论： ①一维的傅里叶变换的系数是而不是。 ②傅里叶变换式中的t可看成参变量。因此，当原来坐标空间的波函数不含时间变量时， 即相当于的情况，变换式的形式保持不变。 ③不难证明，若是归一化的，则经傅里叶变换得到也是归一化的。 2．设在时，粒子的状态为 求粒子动量的平均值和粒子动能的平均值。 [解]方法一：根据态迭加原理和波函数的统计解释。任意状态总可以分解为单色 平面波的线性和，即，展开式的系数表示粒子的动量为p时的几率。知道了几率分布函数后，就可按照 求平均值。

在时，动量有一定值的函数，即单色德布罗意平面波为，与的展开式比较可知，处在状态的粒子动量可以取 ，而， 粒子动量的平均值为 A可由归一化条件确定 故 粒子动能的平均值为 。 方法二：直接积分法

根据函数的性质，只有当函数的宗量等于零时，函数方不为零，故的可能值有 而 则有及。 讨论：①由于单色德布罗意平面波当时不趋于零，因此的归一化积分是发散的，故采用动量几率分布的概念来求归一化系数。 ②本题的不是平方可积的函数，因此不能作傅氏积分展开，只能作傅氏级数展 开，即这时对应于波函数的是分立谱而不是连续谱，因此计算积分， 得到函数。 ③在连续谱函数还未熟练以前，建议教学时只引导学生按方法一做，在第三章函 数讲授后再用函数做一遍，对比一下，熟悉一下函数的运算。 3．一维谐振子处在 的状态，求： （1）势能的平均值； （2）动量的几率分布函数； （3）动能的平均值 [解]先检验是否归一化。 是归一化的。 （1）