时域薛定谔方程的FDTD计算.pdfVIP

中国西部科技 2010年06,q(中旬)第09卷第17期总第214期 时域薛定谔方程N,3Hm)计算 宋永 (四川省阿坝师范高等专科学校，四川阿坝611741) 摘要：在量子力学中，合时薛定谔方程的求解是比较困难的，往往需要做大量的数学变换，而FDTD算法的引入则可轻 松解决这个问题。本文主要介绍FDTD算法在量子力学中的应用及其算法设置。 关键词：薛定谔方程；FDTD；数值 671—6396．201 0．17．009 D01．1O．3969／j．issn．I l 时域薛定调方程 v(胁，，出～。(t)——(式5) 量子力学告诉我们时域薛定谔方程具有如下形式“一： 整理后的迭代格式为： m掣=去掣+脚㈤一(式D ’I，批)州峨k 其中，参数、I，“1)为状态变量，是位置和时间的函 ’|，r’o)州(k，磊h苦『o+1)一却孙)w孙一1)]一iAtvow抛)l 数，本身没有直接的物理含义，但是，所有相关的物理参 ——(式6) 数均可以由它决定，为复数，描述了微观粒子的状态随时 间变化的规律，是微观物理世界物质运动的基本规律，就 若知道，出时刻及周围函数值，就可以得到函数V在 像牛顿定律在经典力学中所起的作用一样。V(x)为势函O+l冷时刻的函数值。 3 激励源的设置 数，单位为eV，h(=1．054×104J—s)为Plank常数，m为 粒子质量。 激励源的设置是FDTD算法中必不可少的部分，但鼍子力 学中的激励源与电磁学和声学中的不同。由于波函数具有统 2 时域有限差分法——FD皿m町 计特性，任何时刻空间上的任何一位置均有粒子或体系存在 麦克斯韦方程组分别由两个旋度方程和两个散度方程 构成，其中两个旋度方程为： 的可能，即任何一空问位置上均有体系存在的几率，故量子 力学中的激励源设置其实是找到某-‘个时刻的量子态，得到 V×盲=署J 该时刻空间的几率分布。但往往体系的空间几率分布并不能 V五=一詈i——(式2) 其中E为电场强度：D为电通量密度：H为磁场强度； 确定，这就需要构造一个模型来模拟体系运动或者找到一个 近似量子态即试探波函数。该态的选择必须尽可能的接近体 B为磁通量密度；J为电流密度；J．为磁流密度。 系的真实态，通过体系的反复修正求解可以得到比较精确的 各向同性线性介质中存在如下本构关系： 真实波函数。常常选定计算初始时刻该体系的状态波函数来 D=￡E 作为开始激励，本文所采用粒子的初始状态均为宽度为的 i：届I 高斯源，即： J=oE l——(式3) Z=a．刮 卜。@o)=乒㈣娜e芈) 其中e、1．t分别为介电常数(F／m)和磁导系数 III，，@0)：乒∽咖(!警) (H／m)；0、0．分别表示电导率(S／m)和导磁率(Q 同时，迭代过程中，需要始终保持V的归一化，即： ／m)。 时域薛定谔方程可分裂成两部分