时域有限差分法在势垒贯穿问题中的应用
第27卷第5期2007年9月
云南师范大学学报
JournalofYunnanNormalUniversity
V01.”No.5
Sep.2007
时域有限差分法在势垒贯穿问题中的应用‘
谷桂初,周庆
(云南大学物理系,云南昆明650091)
摘要:时域有限差分(Finite
Difference
TimeDom丑in,FDTD)法在量子力学问题中的应用越来越引起
人们的关注。系统介绍了该方法的基本原理、边界处理方法和初始试探波函数的选取。在此基础上.计算模拟了三种不同的势垒如方势垒、三角形势垒和谐振势垒下系统的能量变化、概率变化,再现了量子隧道效应这一重要物理现象。进一步表明了运用FuID法来解决量子力学阃题是方便可行的.并且以系统平均能量方法处理截断边界能得到比较理想的结果。对于解决其它量子力学问题具有一定的指导意义。
关键词:时域有限差分法;势垒贯穿;量子隧道效应中围分类号:0411.3。0413.1
文献标识码:A
文章编号:1007—9793(2007)05-0043—06
时域有限差分(Finite
DifferenceTimeDo-
互转化、概率变化过程和规律进行探讨分析。
FlYlD法对于解决其它量子力学问题具有一定的指导意义。
main,FDTD)法…在计算电磁学中得到了充分的应用和发展”J,同时在量子力学中的应用研究也
引起了人们的关注”oJ。在量子力学问题的计算
分析中,直接求解Schrijdinger方程往往比较困难。时域有限差分法的提出,为解决量子力学问
1
FDTD法
题提供了另一条相对简单方便的参考途径。其基本思想是:将含时SchrOdinger方程转换成有限差分形式。建立时间离散的递进序列,采用合适的截断边界条件和初始概率波函数在时域上进行求解。
势垒模型被广泛应用于物理、化学和分子生
物学”1。是量子力学中的一个经典模型之一,揭
1.1
SchtOdinger方程的有限差分形式
r
考虑一维含时Schrfdinger方程聃J
”
[一丽112嘉川圳帅’I)一11。0.(刈),11)
其中:m为微观粒子的静止质量;h为Planck常数,tI=1.05459x10一J・。;i刮‘二丁;以x)为势能函数;驴(*,t)为描述粒子运动的概率波函数。由
示了量子力学不同于经典物理学的一个重要现象——量子隧道效应。故对势垒问题的研究有着
重要的意义。
于廿(*,t)具有复数形式,可用其实部o(Ⅳ,1)和虚
部b(x,t)来表示,即
廿(*,I)=口(g,t)+i・b(*,£)
(2)
本文系统的介绍了时域有限差分法的基本原
理、边界处理方法和初始试探波函数的选取原则。并采用FDTD法对不同的势垒贯穿问题进行研究分析。对不同形状的势垒如方势垒、线性势垒(三角形)、谐振型势垒等的系统动能和势能的相
将式(2)代入式(1),则一维SchriJdinger方程变为
争(州)一丽11争(列)+古以小b(x,t)(孙)争∽牡嘉务∽小如小出力(36)
在足够大的有限区间内,将空间£和时间t
・收稿日期.*2007—04—04
作者简介:谷桂初(1982一),男,湖南省娄底市人,硕士生,主要从事时域有限差分法在计算电磁学和量子力学中
的应用研究.
’
通讯作者:周庆(1963一),女,贵州人,博士,教授,主要从事物理电子和光电子学的研究.
・“・
云南师范大学学报(自然科学版)第”卷
进行离散化,分别用^和n表示空间和时间的节
点。血和血表示空间和时间的差分增量,则
妒(*,I)=砂(t・缸,17,・△t)=妒(k,n)=驴4(k)(4)
根据有限差分近似m,有
警*击渺(㈨ⅢⅥ_1)一缈㈤]
掣一罄1旷Ⅵ)叫㈤]
(5a)(Sb)
故式(34)和式(3b)可近似写成
o一(%):n一似)一;鱼[6^(t+1)+扩(%一1)一Zm
26^(圳+_产以”叭^)
(6。)
6n+l(I)=64(&)+s‘.mlt[。4(^+1)+口“(k-1)一
2a“(I)]一!二—{÷业y(^)・n4(%)
(66)
其中”=畚。(6口)式和(66)式即为一维
SchriMinger方程的有限差分近似形式。1.2初始试探波函数的选取
初始试探波函数的选取,原则上要使其能够近似地描述粒子的初始运动。对于双原子分子振
动运动的初始试探波函数,一般采取谐振子的经
典波函数来处理。对粒子在自由空间中运动,其
概率一般满足高斯型分布,为此采取高斯波包型波函数为其初始试探波函数。其实部和虚部形式
分别如下:
口(£)=唧[一o.5・(X了--X0)2].c*
(兰等立)
m)
6(x)=exp
c-0.5..X-口Xo)2].血
(型竿生)
(76)
其中:^为其波长;口为高斯波包宽;‰为高斯脉冲的中心。在本文中,初始试探波函数需满足归一化条件
亡…刈∥出=l
(8)
即在运用FDTD法计算模拟时,初始渡函数需乘
以其归一化系数。1.3边界截断处理
在FDTD中,不可能也没必要计算无限大的区间,所以必须进行边界截断处理。在此采用系统平均能量的截断处理方法。对于概率波函数
一旦i争(圳)=E・班(州)
(9)
其中:E为系统的总能量。将式(2)代入(9)式,有
一孕鲁[n(刈)+i・6(刈)]=E・[。(¨)+f
・6(z,£)]
(10)
等式两边实郝和虚部分别对应相等,则
h骂竽=层.6(¨)
(1ln)
曲女唑掣:E.口(州)(1lb)
令边界端点(xmi。和K一)处的能量E“(k。)和E“
(k)为此时系统的平均能量E4,则由式(56)得
到边界条件的有限差分近似形式
∽cu圳cu一扣他,。2旬P“(&,。)=。4(Km.n)+铲1・6“(k)
ktr1(k)2口4(k)+iIuE・扩(k)
ck,圳ck,一扣耐ck,n2∞
其中:系统能量均值为
i=(吵㈣砂>=一T+一V
(13)
动能均值为
予砌I岫=(小差喜㈨
一鬟亡(班’磐)出
(14)
势能均值为P=(廿Ivl廿)
=亡y(*)I砂(钳)l“
(15)
第5期各桂初,等:时域有限差分法在势垒贯穿问题中的应用・45・
考虑单个电子与势垒的作用过程,如图1所
2势垒贯穿问题的计算分析
2.1计算模型与计算参数的选取
示。电子质量m=9.1x10—31蛔,初始动能KE
=lOeV,势垒最高vo=15eV.图l中P1、忍、P3分
别表示三个作用区间内找到电子的概率。
KElloeV
●—v、,o●
Pl
Kl
毳
K2
籍
I
c
P3
K3
K.
图1不同势垒(一方势垒、——三角形势垒、…谐振势垒)的作用模型
Fig.1
Diffvrentmodelsabouttheintemcllonoftheelectronanddifferent
potential嫡em
includingreeta“gul“potential—bar-
tier。deltoidpotential—barrierandsyntonicpotential—harrier.eapeclive]y
图1中各势垒函数为(单位:e矿)方势垒
(16)
作用过程。
图2给出了不同时刻的概率波函数。(^)随着空间节点的变化关系,其中横轴表示空间节点。可以看出,当时间步N=l步时,图中表示的是初始试探波函数随着空间节点k的变化关系,此时
,0.15K一45
墨≤^≤KoK≤t≤墨
(17)
电子还没有与方势垒作用;随着时间的推移,在Ⅳ=4500步左右,电子与方势垒开始作用;N=10000步时,其作用效果比较明显;在N=30000步的图形中,可以清楚的看出。在方势垒的右边有
y(k)=2-0.15k+75
1
0
k<%,I>墨
0015酽+1.2/(一2259,≤%≤墨
Hk)={。(z8)
【0k>玛,&(墨
r—n
。
一定的波形,说明电子有一定的概率可以穿透方势垒,这正是量子隧道效应的具体体现。
同时从量子力学观点来看,由于粒子的波动
在本文计算中,考察区问O≤z≤7rwz,空间增量△x=O.Olnm,系数5=1000sec・rn,一,总时间迭代步数取30000步。
对于初始波函数,令波长It=39Az;高斯波包
性,问题与波透过一层介质相似,有一部分波穿
过,一部分波被反射回去。故随着粒子与方势垒
的作用,图中左边区域(式。<k(置)概率波函数
的最大峰值减小,波形密集度增加,势垒区和右边
宽口=30A。;高斯脉冲的中心x0=150Az。
2.2计算结果和分析
采用FDTD法,计算模拟得到电子与势垒的
区域有一部分波。按照泼函数的统计诠释,粒子有一部分概率穿过势垒,有一定的概率被反射回
去。
・46・
云南师范大学学报(自然科学版)第”卷
图2不同时刻IN=1.4500.10000.300GO)的撬率波函数a(k)随空间节点k的变化关系.其中横轴表示空间节点
Fig.2
Theprobability
w州efunction口【k)with
respecttospace—stepk址differenttime—stepⅣ。whereⅣ;1-4500,10(】00
30000respec6vely.
一图3给出了电子与不同势垒的作用过程中,系统的动能期望值KE、势能期望值PE和各个区
在初始区域中的概率变化;P2表示的是粒子在势垒区的存在几率;1'3是粒子穿透势垒到达势垒右
间内概率(Pl、P2、P3及Psum)随时间的变化关
系,其中横轴表示时间步。
边的概率;Psum是计算区间总的概率,为Psum=
Pl+P2+P3。对于三个势垒中的概率变化,总有Psum=1,即不论各个区域概率如何变化,总的概
从系统的动能期望值与势能期望值的关系图
中可以看出,动能的减少,恰是系统势能的增加,系统总的能量不变。其中,与三角形势垒作用下,各能量变化幅度最大,而方势垒中的能量变化幅
率不随时间改变,恒等于1。这不仅是归一化条
件限制的结果,也正是概率守恒(或粒子数守恒)的表现。比较三个势垒中的概率曲线,发现三角
度最小。能量的变化反映了粒子与势垒的相互作
用。同时,在能量变化峰值的两侧,图形并不对
形势垒中的P3幅值最大,Pl和P2变化最大,而方势垒中的概率变化最小。也即三角形势垒的贯
穿系数最大。P3不为零,这正是量子隧道效应作用的结果。
称,图中方势垒和谐振势垒中比较明显。说明对
称型的势垒,能量的变化并不完全对称。这也反映
了微观粒子量子力学中的不确定性。
图3中的概率变化曲线中:Pl反映的是粒子
联系能量变化和概率变化来看,能量变化时,
各区域内的概率也相应发生变化。当粒子与势垒
第5期谷桂初,等:时域有限差分法在势垒贯穿问题中的应用
・47・
作用时,动能和势能的转化关系,恰恰与各区域内过程,得到了概率波函数、系统能量的变化规律,符合理论上的解析分析,说明了用系统平均能量
的方法来处理截断边界是可行的。
概率的变化相对应。能量变化和概率变化是对应
进行的。
图2和图3清楚的描述了粒子与势垒的作用
¨¨¨¨¨¨
●5●■w0■'卸¨2¥00e25No,∞¨
lH
¨
¨¨¨¨¨
●¥041
伸●呻
I,¨
兰●■●●
2_
符¨O
枷¨
●
郅∞¨●¨ts000
!●■,●●
2嘲●
25口一
埘¨
●580椰●伸|s000
IDtIHⅢH
I
x
2●¨●签0N’∞¨●
5蚺伸●∞
',时●26000∞O■3∞●●
圈3能量期望值和各区域概率髓时间的变化关系。其中横轴表示时间步
Fig.3
Theexpectationenergy
kineticandpotential
energy.-espoctivelyJofthesystemandtheprobabilityin
different吓㈣
withrespecttotime。wherethe
axismeanstime—step
・48・
云南师范大学学报(自然科学版)
pmbkrns
in
valvingMaxwell0
011
第27卷
equation
in
isotropie
3结论
media[J].IEEETramAntennasandPropagation.
本文系统地介绍了FDTD法的基本原理、边界处理和初始试探波函数的选取,发展了FDTD
[2]
1966,14(I):302—307.
葛德彪,闰玉渡.电磁渡时域有限差分法[M].西安:西安电子科技大学出版社,2002.
[3]黄整,陈晓敏,朱正和.时域有限差分法对双原子分
子振动光谱计算的应用[J].原子与分子物理学报,2000,17(4):576—582.[4]Chang
S,H.andTalloveA.,Finite—differencetime
8
法在势垒问题中的应用.以系统平均能量方法处
理截断边界.用此方法计算模拟了在不同势垒
(方势垒、三角形势垒、谐振势垒)下粒子波函数、
系统动能期望值、势能期望值、各区域内概率的变
化关系.再现了微观世界量子隧道效应这一重要
物理现象。进一步表昵,运用FlYrD法来解决量
—domainmodeloflasingactioninelectronatomic3827—3833.
four—leveltwo—Exp.,2004,12:
system[J],Optics
子力学问题是方便可行的,并且以系统平均能量
方法处理截断边界能得到比较理想的结果。对于
[5]范康年.物理化学(第二版)[M].北京:高等教育出
版社,2005.[6]
其它难于直接求解的量子力学问题的解决提供了
一种新的参考途径。
曾谨言.量子力学(卷I)[M].北京:科学出版社,
2000.
[7]ThomasJW.Numerical
参考文献
PartialDifferentialEquations:
Phb.
FiniteDi∞陆rence
Me&ads[M].Be硒ms:World
1ishingCorporation,1997.
[I]YeeKS.Numericalsolutionofinitialboundary
value
TheapplicationoftheFiniteDifferenceTime
Domainmethodintheproblemofpotentialbarrier
GUGui—chu,
(Depatimem
ZHOUQing
650091,China)
inquantum—mechahieal
ofPhysics,YunnanuniveBity,Kunming
Abstract:TheapplicationoftheFiniteDifferenceTime
Domain(F啪)method
problemsbecom∞moreandmoreinteresting.11lepostulate.theboundarycut—ogconditionandtheinitialtrialwavefunetiona"desefibed.Based
on
theFDTD
method.three
differentpotential—ballierproblemsin—
eludingrectangularpotential—barrier。deltoidpotential—barrierandsyntonicted.Theenergyvariationandtheprobabilityvariation
are
potential—barrier
important
are
investiga-
calculatedandmodeled.the
phenomenon
ofquantumtunneleffectreappeared.ItisconvenienttoapplytheFDTDmethodin
quantum—mechanical
problemsand
theresults
are
acceptablewhileusingtheboundarycut—o仃based
Oil
theexpectationenergyof
thesystem.Anditispotentiallyfurtber
extendedtootherquantum—mechanicalproblems.
Keywords:theFiniteDifferenceTimeDomainmethod,potential
barrier
penetration,quantum
tunneleffect
时域有限差分法在势垒贯穿问题中的应用
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
谷桂初, 周庆, GU Gui-chu, ZHOU Qing云南大学,物理系,云南,昆明,650091
云南师范大学学报(自然科学版)
JOURNAL OF YUNNAN NORMAL UNIVERSITY(NATURAL SCIENCES EDITION)2007,27(5)4次
参考文献(7条)
1.Yee K S Numerical solution of initial boundary value problems in volving Maxwell's equation inisotropic media 1966(01)
2.葛德彪;闫玉波 电磁波时域有限差分法 2002
3.黄整;陈晓敏;朱正和 时域有限差分法对双原子分子振动光谱计算的应用[期刊论文]-原子与分子物理学报2000(04)
4.Chang S H;Taflove A Finite-difference time-domain model of lasing action in a four-level two-electron atomic system 20045.范康年 物理化学 20056.曾谨言 量子力学 2000
7.Thomas J W Numerical Partial Differential Equations:Finite Difference Methods 1997
本文读者也读过(10条)
1. 刘岩.廖炽坚.LIU Yan.LIAO Chi-jian 有限时域差分法对量子隧道效应的研究[期刊论文]-甘肃联合大学学报(自然科学版)2010,24(3)
2. 邓玲玲.张雅男.DENG Ling-ling.ZHANG Ya-nan 基于时域有限差分法的势垒透射过程的蒙特卡罗模拟[期刊论文]-大学物理2010,29(1)
3. 李小珠.刘益民 在均匀磁场中薛定谔方程的求解[期刊论文]-韶关学院学报2003,24(3)
4. 李淑红.郑忠喜.高迎芳.刘文森 含时受迫谐振子薛定谔方程的精确解[期刊论文]-山西大学学报(自然科学版)2003,26(2)
5. 张红梅.刘德.ZHANG Hong-mei.LIU De 传递矩阵方法与矩形势垒的量子隧穿[期刊论文]-河北科技大学学报2006,27(3)
6. 田红.田旭 二维谐振子薛定谔方程的非定态解[期刊论文]-中南民族学院学报(自然科学版)2001,20(z1)7. 袁通全.YUAN Tong-quan 利用相空间中薛定谔方程的解构造分布函数[期刊论文]-河池学院学报2008,28(5)8. 胡文燕 一类非线性薛定谔方程解的爆破性质的研究综述[学位论文]2007
9. 张学龙.王明泉.Zhang Xuelong.Wang Mingquan 用路径积分方法研究势垒散射[期刊论文]-光子学报2000,29(8)10. 胡先权.许杰.罗光.郑瑞伦.HU Xian-quan.XU Jie.LUO Guang.ZHENG Rui-lun 幂函数叠加势的径向薛定谔方程的解析解[期刊论文]-原子与分子物理学报2007,24(5)
引证文献(4条)
1.刘岩.廖炽坚 有限时域差分法对量子隧道效应的研究[期刊论文]-甘肃联合大学学报(自然科学版) 2010(3)2.邓玲玲.张雅男 基于时域有限差分法的势垒透射过程的蒙特卡罗模拟[期刊论文]-大学物理 2010(1)
3.谷桂初.陈玖福.高利霄.周庆 时域有限差分法计算薛定谔方程中的数值稳定性条件[期刊论文]-云南师范大学学报(自然科学版) 2008(3)
4.陈玖福.GO Gui-chu.高利霄.ZHOU Qing 基于时域有限差分法研究含时受迫谐振子相干态问题[期刊论文]-云南师范大学学报(自然科学版) 2008(4)
本文链接:http://d.wanfangdata.com.cn/Periodical_ynsfdxxb200705011.aspx