Numerov算法求解一维薛定谔方程研究

定态薛定谔方程讲义

定态薛定谔方程 一、定态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ω/()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 =???，

薛定谔方程

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

非线性薛定谔方程数值解的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 强大的符号计算功能，完全可以实现分步傅立叶数值算法来对脉冲形状和频谱进行仿真。 一般情况下，光脉冲沿光纤传播时受到色散和非线性效应的共同作用，假设当传输距离 很小的时候，两者相互独立作用，那么，根据这种思想可建立如下分步傅立叶数值算法的数 学模型： （I ）

§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)是描述势场的函数，假设不随时间变化。

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 =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα

薛定谔方程及其解法

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

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

第二章薛定谔方程习题

第二章 薛定谔方程 习题 (课本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?表示向外传播的球面波。

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

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

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

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

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

量子笔记1 —— 一维薛定谔方程

量子笔记1 —— 一维薛定谔方程 给出某种一维势，求解一维薛定谔方程的束缚定态解及其能级的题目是常见的量子力学的题型之一，这种题型的求解虽有其固有模式，但具体处理过程中也牵涉到很多技巧和要注意之处。下面我通过两个例子来试图对其解题模式和某些解题过程中的常见技巧和经验作出一个概括性的总结，作为量子力学复习的第一阶段的一个阶段性小结。 例一. 质量为μ的粒子在一维势场 ()()???><=''+-=0,V 0 0V V V 0 x x x x ，，αδ 中运动，其中α与0V 均为实数。（1）试给出存在束缚态的条件，并给出其能量本征值和相应的本征函数。（2）给出粒子处于0>x 区域中的概率，它是大于1/2，还是小于1/2，为什么？ 解：在此势场中的束缚定态能量E<0， 令202E V 2,E 2 ）（）（-=-=μγμβ （1） 不包括x=0点的定态方程为 0),()(0),()(222222>=<=x x dx x d x x dx x d ψγψψβψ （2） )(x ψ满足条件 0)(),0()0(=±∞=-+ψψψ （3） )0(2)0()0(2ψμαψψ - ='-'-+ （4） 在势阱中，能量要小于势阱边缘的势能，否则就不能形成束缚态，而是游离态了。这里的势阱是一个δ势，所以要有束缚态，并需要求E<0。 本来一维定态薛定谔方程是ψψE V dx d m =+-)2(22 2 ，用熟了就可以直接像上面那样写了。特别的，当粒子限定在圆周上运动时，或者写自由转子的薛定谔方程时，只要将上式中的m 换成I ，x 换成?就可以了。其中2R I μ=是粒子的转动惯量。即)()()(R 222 22?ψ?ψ??μE V d d =??? ? ??+-

薛定谔方程及其解法 - 副本

薛定谔方程 一． 定义及重要性 薛定谔方程（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 =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα

mathematic求解薛定谔方程程序

5.编程求解薛定谔方程： BeginPackage["QuantumWell`"] Clear[PsiSym,PsiASym,Spectrum] PsiSym::usage= "PsiSym[x_ ,k ,a ] determines the symmetric eigenfunction for a potential well of depth -V0. The input parameter k fixes the energy and 2a the width of the well. Psisym is useful for a numerical representation of eigenfunctions." PsiASym::usage= "PsiASym[x_ ,k_ ,a_] determines the antisymmetric eigenfunction for a potential well of depth-V0. The input parameter k fixes the energy and 2a the width of the well. PsiASym is useful for a numerical representation of eigenfunctions." Spectrum::usage="Spectrum[V0_ ,a_]calculates the negative eigenvalues in a potential well. V0 is the potential depth and 2a the width of the well.The eigenvalues are returend as a list and are available in the variables lsymand lasym as replacement rules. The corresponding plots of eigenfunctions are stored in the variables Plsym and Plasym. The determining equation for the eieenvalues is plotted." (*-一define global variables-一*) Plsym::usage= "Variables containing the symmetric plots of the eigenfunctions." Plasym::usage ="Variables containing the antisymmetric plots of the eigenfunctions." lsym::usage= "List of symmetric eigenvalues." lasym::usage ="List of antisymmetric eigenvalues." k::usage ="Eigenvalue." Begin["`Private`"] (*一symmetric eigenfunctions一*) Psisym[x ,k ,a ]:=Module[{kapa, Al]},Kapa = k Tan [k a]; (*一normalization constant一*) A1=1/Sqrt[a Exp[-2 a kapa] (1+1/(kapa a)+kapa/(k^2 a)+Kapa^2/k^2)]; (*一define the three domains of solution一一*) Which[-Infinity

第二章薛定谔方程

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

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