薛定谔方程数值解.ppt

《薛定谔方程数值解.ppt》由会员分享，可在线阅读，更多相关《薛定谔方程数值解.ppt（33页珍藏版）》请在三一办公上搜索。

1、计算物理,http:/125.217.162.13/lesson/ComputationalPhysics,薛定谔方程数值解,薛定谔方程数值解,薛定谔方程定态方程的矩阵解法含时方程的解法非线性薛定谔方程解法薛定谔方程的有限元方法,薛定谔方程(1/1),薛定谔方程,单粒子,多粒子,定态薛定谔方程(势能不显含时间),一维的单粒子,定态方程的矩阵解法(1/9),实对称矩阵的对角化定理：如果 A 是实对称矩阵，那么存在正交矩阵 R，使得,雅可比方法：基于上述定理，用一系列简单的正交矩阵 RK，逐步将 A 对角化，即选择 RK，令,取 A0=A，使得当 K 时，AK diag(l1,l2,ln)本征值：

2、l1,l2,ln本征向量：,定态方程的矩阵解法(2/9),矩阵,对角化 22 实对称矩阵 A,定态方程的矩阵解法(3/9),对角化 nn 实对称矩阵 A,定态方程的矩阵解法(4/9),例：计算 33 实对称矩阵 A 的本征值和本征向量,A0=A，选 p=1,q=2,选 p=1,q=3,第 9 次,定态方程的矩阵解法(5/9),久期方程方法,例：计算实对称矩阵的本征值问题,久期方程和本征值,本征向量,定态方程的矩阵解法(6/9),定态薛定谔方程的矩阵解法有限差分法例：一维无限深势阱,定态薛定谔方程的差分格式,差分方程的实对称矩阵和本征值问题,定态方程的矩阵解法(7/9),波函数,有限差分法的步骤

3、将定态薛定谔方程转化为差分格式写出差分方程的实对称矩阵，并对角化,定态方程的矩阵解法(8/9),希耳伯特空间的方法例：一维无限深线性势阱,希耳伯特空间的基矢,哈密顿算符和矩阵元,定态方程的矩阵解法(9/9),对角化(以 为能量单位)，结果分析,步骤选择适当的表象(即基矢)，推导哈密顿矩阵元计算哈密顿矩阵，并对角化,含时方程的解法(1/10),非本征态的时间演化特点：初始态 系统的本征态,解法1：有限差分方法解多维扩散方程一维含时薛定谔方程的差分格式,利用 k 时的 y 值，求 k+1 时的 y 值要求解线性方程组隐式的,含时方程的解法(2/10),边界条件束缚态：y=0非束缚态：假设初始态是束

4、缚态，x 足够大，在 T 内，x 左右边界处的 y=0,方程组的矩阵形式：,含时方程的解法(3/10),含时方程的解法(4/10),含时方程的解法(5/10),例：先将一个粒子用谐振子势束缚在基态，然后放入无限深势阱的中央，求其波函数 Y 随时间的演化,初始态(即 k=0)：谐振子的基态,由 k=0 的差分方程组，求 k=1 时刻的波函数方程,未知量,已知量,含时方程的解法(6/10),含时方程的解法(7/10),含时方程的解法(8/10),由 k=1 的差分方程组，求 k=2 时刻的波函数方程,未知量,已知量,由 k=2 的差分方程组，求 k=3 时刻的波函数,n 维含时薛定谔方程的差分格式

5、：Nn Nn 的矩阵,含时方程的解法(9/10),解法2：(希耳伯特空间的方法)将初始态在基矢上展开势能函数不含时(Q 表象，基矢是 un),写成矩阵形式,，Y 和 H 都是矩阵,当 Q=H 时，un 是哈密顿量的本征态，H 是对角阵,含时方程的解法(10/10),当 Q H 时，un 不是哈密顿量的本征态，H 非对角,例：先将一个粒子用谐振子势束缚在基态，然后放入无限深势阱的中央，求其波函数 Y 随时间的演化,初始态：谐振子的基态,基矢：无限深势阱的能量本征函数,常数 cn,非线性薛定谔方程解法(1/10),非线性薛定谔方程(Gross-Pitaevskii方程),实质是非线性偏微分方程，一

6、般没有解析解希耳伯特空间的迭代法(定态)：非线性项线性化,非线性薛定谔方程解法(2/10),例：求解 f=x 2 的一维 G-P 方程(0 x l),方程：,基矢：,哈密顿矩阵元,非线性薛定谔方程解法(3/10),l=5,a=5,非线性薛定谔方程解法(4/10),有限差分的迭代法(一维定态)：非线性项线性化,非线性薛定谔方程解法(5/10),例：求解 f=x 2 的一维 G-P 方程(0 x l),迭代方程,系数矩阵,非线性薛定谔方程解法(6/10),有限差分法(一维含时)二元非线性方程组的迭代法：非线性问题线性化,(关键：Jocobi矩阵 J，初值要接近解),非线性薛定谔方程解法(7/10),多元非线性方程组,(关键：Jocobi矩阵 J，初值接近解),非线性薛定谔方程解法(8/10),例：解非线性方程,Jocobi矩阵,初值,迭代,非线性薛定谔方程解法(9/10),有限差分法(一维含时),非线性薛定谔方程解法(10/10),有限差分法(一维含时)(续),作业,PPT.22 用Fortran实现，要求将|Y(x,t)|2 输出到数据文件,