薛定谔方程
薛定谔方程（Schrodinger equation）是由奥地利物理学家薛定谔提出的量子力学中的一个基本方程，也是量子力学的一个基本假定，其正确性只能靠实验来检验。它是将物质波的概念和波动方程相结合建立的二阶偏微分方程，可描述微观粒子的运动，每个微观系统都有一个相应的薛定谔方程式，通过解方程可得到波函数的具体形式以及对应的能量，从而了解微观系统的性质。
1定义
薛定谔方程
薛定谔方程（Schrodinger equation）又称薛定谔波动方程（Schrodinger wave equation）在量子力学中，体系的状态不能用力学量（例如x）的值来确定，而是要用力学量的函数Ψ（x,t），即波函数（又称概率幅，态函数）来确定，因此波函数成为量子力学研究的主要对象。力学量取值的概率分布如何，这个分布随时间如何变化，这些问题都可以通过求解波函数的薛定谔方程得到解答。这个方程是奥地利物理学家薛定谔于1926年提出的，它是量子力学最基本的方程之一，在量子力学中的地位与牛顿方程在经典力学中的地位相当。
薛定谔方程是量子力学最基本的方程，亦是量子力学的一个基本假定，其正确性只能靠实验来确定。
2方程概述
量子力学中求解粒子问题常归结为解薛定谔方程或定态薛定谔方程。薛定谔方程广泛地用于原子物理、核物理和固体物理，对于原子、分子、核、固体等一系列问题中求解的结果都与实际符合得很好。
薛定谔方程仅适用于速度不太大的非相对论粒子，其中也没有包含关于粒子自旋的描述。当涉及相对论效应时，薛定谔方程由相对论量子力学方程所取代，其中自然包含了粒子的自旋。.薛定谔提出的量子力学基本方程。建立于1926年。它是一个非相对论的波动方程。它反映了描述微观粒子的状态随时间变化的规律，它在量子力学中的地位相当于牛顿定律对于经典力学一样，是量子力学的基本假设之一。设描述微观粒子状态的波函数为Ψ（r，t），质量为m的微观粒子在势场V（r，t）中运动的薛定谔方程为。在给定初始条件和边界条件以及波函数所满足的单值、有限、连续的条件下，可解出波函数Ψ（r，t）。由此可计算粒子的分布概率和任何可能实验的平均值（期望值）。当势函数V不依赖于时间t时，粒子具有确定的能量，粒子的状态称为定态。定态时的波函数可写成式中Ψ（r）称为定态波函数，满足定态薛定谔方程，这一方程在数学上称为本征方程，式中E为本征值，是定态能量，Ψ（r）又称为属于本征值E的本征函数。
薛定谔方程是量子力学的基本方程，它揭示了微观物理世界物质运动的基本规律，就像牛顿定律在经典力学中所起的作用一样，它是原子物理学中处理一切非相对论问题的有力工具，在原子、分子、固体物理、核物理、化学等领域中被广泛应用。
3提出历史
当法国物理学家德布罗意的“微观粒子也像光一样具有波粒二象性”的假说被美国物理学家
戴维逊和革末利用“电子的晶体粉末散射实验”证实后，薛定谔通过类比光谱公式成功地发现了可以描述微观粒子运动状态的方法——薛定谔方程
4作者简介
埃尔温·薛定谔（Erwin Schrodinger，1887—1961年）1887年8月12日出生于奥地利首都维也纳。1906年至1910年，他就学于维也纳大学物理系。1910年获得博士学位。毕业后，在维也纳大学第二物理研究所从事实验物理的工作。第一次世界大战期间，他应征服役于一个偏僻的炮兵要塞，利用闲暇时间研究理论物理。战后他仍回到第二物理研究所。1920年他到耶拿大学协助维恩工作。1921年薛定谔受聘到瑞士的苏黎世大学任数学物理教授，在那里工作了6年，薛定谔方程就是在这一期间提出的。
埃尔温·薛定谔
1927年薛定谔接替普朗克到柏林大学担任理论物理教授。1933年希特勒上台后，薛定谔对于纳粹政权迫害爱因斯坦等杰出科学家的法西斯行为深为愤慨，移居牛津，在马达伦学院任访问教授。同年他与狄拉克共同获得诺贝尔物理学奖。
1936年他回到奥地利任格拉茨大学理论物理教授。不到两年，奥地利被纳粹并吞后，他又陷入了逆境。1939年10月流亡到爱尔兰首府都柏林，就任都柏林高级研究所所长，从事理论物理研究。在此期间还进行了科学哲学、生物物理研究，颇有建树。出版了《生命是什么》一书，试图用量子物理阐明遗传结构的稳定性。1956年薛定谔回到了奥地利，被聘为维也纳大学理论物理教授，奥地利政府给予他极大的荣誉，设定了以薛定谔命名的国家奖金，由奥地利科学院授给。
5具体介绍
数学形式
一维薛定谔方程
三维薛定谔方程
定态薛定谔方程
单粒子薛定谔方程的数学表达形式
这是一个二阶线性偏微分方程，ψ(x，y，z)是待求函数，它是x,y,z三个变量的复数函数（就是说函数值不一定是实数，也可能是虚数）。式子最左边的倒三角是一个算符，意思是分别对ψ(x，y，z)的梯度求散度。
物理含义
这是一个描述一个粒子在三维势场中的定态薛定谔方程。所谓势场，就是粒子在其中会有势能的场，比如电场就是一个带电粒子的势场；所谓定态，就是假设波函数不随时间变化。其中，E是粒子本身的能量；U(x，y，z)是描述势场的函数，假设不随时间变化。薛定谔方程有一个很好的性质，就是时间和空间部分是相互分立的，求出定态波函数的空间部分后再乘上时间部分e^(-t*i*E*2π/h)以后就成了完整的波函数了。
薛定谔方程的解——波函数的性质
简单系统，如氢原子中电子的薛定谔方程才能求解，对于复杂系统必须近似求解。因为对于有Z 个电子的原子，其电子由于屏蔽效应相互作用势能会发生改变，所以只能近似求解。近似求解的方法主要有变分法和微扰法。
在束缚态边界条件下并不是E 值对应的所有解在物理上都是可以接受的。主量子数、角量子数、磁量子数都是薛定谔方程的解。要完整描述电子状态，必须要四个量子数。自旋磁量子数不是薛定谔方程的解，而是作为实验事实接受下来的。
主量子数n和能量有关的量子数。原子具有分立能级，能量只能取一系列值，每一个波函数都对应相应的能量。氢原子以及类氢原子的分立值为：
En=-1/n*2×2.18×10*(-18)J，n 越大能量越高电子层离核越远。主量子数决定了电子出现的最大几率的区域离核远近，决定了电子的能量。N=1，2，3，……；常用K、L、M、N……表示。
角量子数l和能量有关的量子数。电子在原子中具有确定的角动量L，它的取值不是任意的，只能取一系列分立值，称为角动量量子化。L=√l(l+1) ·(h/2π) ，l=0，1，2，……（n-1）。l 越大，角动量越大，能量越高，电子云的形状也不同。l=0，1，2，……常用s，p，d，f，g 表示，简单的说就是前面说的电子亚层。角量子数决定了轨道形状，所以也称未轨道形状量子数。s 为球型，p 为哑铃型，d 为花瓣，f 轨道更为复杂。
磁量子数m是和电子能量无关的量子数。原子中电子绕核运动的轨道角动量，在外磁场方向上的分量是量子化的，并由量子数m 决定，m 称为磁量子数。对于任意选定的外磁场方向Z，角动量L 在此方向上的分量LZ 只能取一系列分立值，这种现象称为空间量子化。LZ=m·h/2π，m=0，±1，±2……±l。磁量子数决定了原子轨道空间伸展方向，即原子轨道在空间的取向，s 轨道一个方向（球），p 轨道3 个方向，d 轨道5 个，f 轨道7 个……。l 相同，m 不同即形状相同空间取向不同的原子轨道能量是相同的。不同原子轨道具有相同能量的现象称为能量简并。
能量相同的原子轨道称为简并轨道，其数目称为简并度。如p 轨道有3 个简并轨道，简并度为3。简并轨道在外磁场作用下会产生能量差异，这就是线状谱在磁场下分裂的原因。ms
粒子的自旋也产生角动量，其大小取决于自旋磁量子数（ms)。电子自旋角动量是量子化的其值为Ls=√s(s+1) ·(h/2π) ，s= 1/2 ，s 为自旋量子数，自旋角动量的一个分量Lsz 应取下列分立值：Lsz= ms(h/2π), ms=±1/2。
原子光谱，在高分辨光谱仪下，每一条光线都是由两条非常接近的光谱线组成，为解释这一现象提出了粒子的自旋。电子的自旋表示电子的两种不同状态，这两种状态有不同的自旋角动量。
电子的自旋不是机械的自身旋转，是本身的内禀属性，是新的自由度。就像质量和电荷一样
是它的内在属性，电子的自旋角动量为：? /2。[1]
对应关系
希尔伯特空间与薛定锷方程
一般，物理上将物理状态与希尔伯特空间上的向量（vector），物理量与希尔伯特空间上的算符相对应。这种形式下的薛定锷方程如右图所示。
薛定锷方程
H为哈密顿算符。这个方程在这个形式下充分显示出了时间与空间的对应性（时间与能量相对应，正如空间与动量相对应，后述）。这种算符（物理量）不随时间变化而状态随时间变化的对自然现象的描述方法被称为薛定谔绘景。与之对应的是海森伯绘景。
空间坐标算符x与其对应的动量算符p满足以下交换关系：
交换关系
所谓的薛定锷表示就是将空间算符直接作为x，而动量算符为下面的包含微分的微分算符：
动量算符[2]
最新薛定谔方程及其解法
关于薛定谔方程 一．定义及重要性 薛定谔方程（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 =-???????++==++==?????=≤≤=++的局部截断误差使以下数值解法的值及确定常数ββα βα
定态薛定谔方程讲义
定态薛定谔方程 一、定态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 ψ 也可
实验三 定态薛定谔方程的矩阵解法
实验三 定态薛定谔方程的矩阵解法 一．实验目的 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) 其中，是质量，是位置，是相依于时间的波函数，是约化普朗克常数，是位势。类似地，在三维空间里，一个单独粒子运动于位势中的含时薛定谔方程为 。(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 ( -?=ψ
固体物理学 1-5-薛定谔方程应用举例II
薛定谔方程应用举例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
非线性薛定谔方程数值解的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 ）
非线性薛定谔方程的孤子解和怪波解
非线性薛定谔方程的孤子解和怪波解 摘要：光纤中光波的传输模型一直是当前研究的热点理论模型之一，从非线性薛定谔方程到金格堡-朗道方程，都试图对其进行更好的阐释，其次对于非线性动力学系统中，非线性薛定谔方程的解有呈现出非常多有趣的特征，对于其中特定解的研究能够让我们了解脉冲演化的本质，所以本文主要从孤子解的传输入手，并且简单介绍了怪波解的解形式。 薛定谔方程又称薛定谔波动方程，是量子力学的一个基本方程，同时又是量子力学的基本假设之一，由奥地利物理学家薛定谔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 =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα
薛定谔方程对氢原子的应用
（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 一维定态薛定谔方程的建立和求解举例
§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法，有限元法将函数定义在简单几何形状（如二维问题中的三角形或任意四边形）的单元域上（分片函数），且不考虑整个定义域的复杂边界条件，这是有限元法优于其他近似方法的原因之一。
非线性薛定谔方程数值解的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 ）
薛定谔方程及其解法 - 副本
薛定谔方程 一． 定义及重要性 薛定谔方程（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 =-???????++==++==????? =≤≤=++的局部截断误差使以下数值解法 的值及确定常数ββα βα
非线性薛定谔方程求解
CHAPTER IV NUMERICAL SOLUTIONS TO THE NONLINEAR SCHR?DINGER EQUATION 4.1Introduction In general,analytical solutions to the full Maxwell wave equation for a nonlinear optical system do not exist.Even numerical solutions to the wave equation are extremely difficult to implement due to the dimensionality of the problem.The vector form of the wave equation is a four-dimensional(three spatial,one temporal),second-order partial differential equation.Thus,approximations based on propagation conditions and experimental results are needed in order to solve an approximate scalar form of the wave equation,i.e.the nonlinear Schr?dinger equation.However,the approximations listed in the previous chapter do limit the generality and validity of the solutions.For example, the condition extreme nonlinearity,as for the case in supercontinuum generation,is a propagation regime where slowly varying envelope approximation may be violated. The purpose of this chapter is to provide an introduction to a very powerful method in numerically solving the NLSE,known as the split-step Fourier method (SSFM)[15].The chapter will begin with a list pointing the advantages of the SSFM
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