根据有限差分法原理,将求解范围划分为一系列等间距的离散节点后,一维薛定谔方程转化为可以用一个矩阵方程表示的节点线性方程组。利用MATLAB提供的矩阵左除命令,即可得到各未知节点的函数近似值。该方法概念简单,使用方便,不需要花费较多精力编程即可求解大型线性方程组。
文章编号: 17—7521)3 020 6288 ( 00— 4—5 0 0
通过有限差分和 MA L B矩阵运算 T A
直接求解一维薛定谔方程
王忆锋唐利斌
(明物理研究所,云南昆明 6 0 2 )昆 5 2 3
摘要:根据有限差分法原理,求解范围划分为一系列等间距的离散节点后,一维薛将定谔方程转化为可以用一个矩阵方程表示的节点线性方程组。利用 MA L B提供的矩 TA 阵左除命令,即可得到各未知节点的函数近似值。该方法概念简单,使用方便,不需要花费较多精力编程即可求解大型线性方程组。 关键词:半导体;量子力学;薛定谔方程;有限差分法; MA L B TA中图分类号: O 1. 431文献标识码: A DO: 1 . 6/ in17—75 00 3 0 I 0 99js . 2 8. 1. . 9 3 .s 6 8 2 0 0
D i e t So ut o f One di e i r c l i n o - m nsona c Sdi e ua i l S hr ng r Eq ton
t r ug nie D i e e e a ATLA B a r x Com put t o h o h Fi t f r nc nd M M t i a in
WANG Yi e g, TANG ib n — n f k— i
( n ig nt ue f hs s K n n un n6 02, h a Ku m n s tt o yi, u mi Y n a 523 C i ) I i P c g n
A bs r t Ac o d n o t e f ie di e e c rn i l, n— m e so a c r di g r e ua i n c n b t ac: c r i g t h n t f r n e p i c p e a o e di n i n lS h S n e q to a e i c n e t d i t e fn d ll e r e u to x e s d i a rx e u to fe h p c s d v d d o v r e n o a s t o o a i a q a i ns e pr s e n a m t i q a i n a t r t e s a e i i i e n
i t e i s o ic e e n d s wi h a q a n e v 1 n o a s r e f d s r t o e t n e u li t r a .Th a rx lf i ii n c m m a d o e e n t e e m t i e t d v so o n
f r d i h
M ATL AB o s,wa e c n b s d t e i e t e f n to p r x ma i n o a h u kn wn n d lf n to r a e u e o d rv h u c i n a p o i to fe c n o o a u c i n.
Th s me h d i sm p e i o c p c n e i n n o e a i n a a o v a g i e re u to s wiho t i t o s i l n c n e t o v n e t i p r t o nd c n s l e l r e l a q a i n t u j n mo e e o t n p o r mm i . r f r si r g a ng Ke y wor: s m io d c o; u n u me ha i s S h 6 i g re u to; n t i e e e m e h d; AT— ds e c n u t r q a t m c n c; c r d n e q a i n f ie d f r nc t o M i
T. AR
1引言
薛定谔方程是量子力学中最基本的方程,
MA L B应用的简捷性不仅体现于特殊函 TA
数的计算]也反映在部分数学物理方程的求,解过程中[引。利用 MA L B强大的矩阵计算功 TA能,可以大幅减少求解一维薛定谔方程所需的
工作量。
其地位有如经典力学中的牛顿第三定律、电磁学中的麦克斯韦方程。碲镉汞等半导体能带结构的计算就需要求解薛定谔方程。一般情况下,薛定谔方程没有解析解,需要对其做数值计算,其求解也一直作为—个典型的数学物理问题被反复分析讨论,且形成了各种解决方案。些方案这无一例外地要解决两个基本问题,一个是算法
2一维薛定谔方程的求解过程分析
一
维定态薛定谔波动方程如下:
dF 2
设计,另一个是通过编写程序来实现算法。
+ )=0
( 1 )
收稿日期: 091—3 20 01 作者简介:王忆锋 ( 6一,男,湖南零陵人,工学士,高级工程师,主要从事器件仿真研究。 Em i w ny 13 9 ) - a: agi l ̄
ng 3 s n c r 6@ i a. o n
IF A E ( O T L)V L3, o3 N R R D M N HY/ O . N . 1,MA 00 R2 1