一种量子化学的有限元并行算法.pdfVIP

2005高性能计算应用大会 226 一种量子化学的有限元并行算法 张海蕾 杨坤 李洪林 王希诚 (大连理工大学力学系，工业装备结构分析国家重点实验室大连116023) 采用并行的迁移式子空问迭代法求解其中的广义特征值问题，计算TLi。和BH分子的基态总能量，获得了 令人满意的结果。 关键词：量子化学，8节点等参元，迁移式子空间迭代法，并行计算 1引言 在量子化学计算中，利用基组展开方法求解 屹(，’)= 等西 分子轨道(M0)得到了广泛的应用，获得了很大的 成功。但基组展开方法的计算精度不高，通常是 通过增加基函数的数目来提高计算精度，这样做 大大增加了计算量，同时容易造成基函数局部过 完备从而线性相关。 其中，J，是分子轨道波函数，0是原子坐标， 有限元方法是一种实施简单的高精度数值 ，．是电子坐标。 计算方法 ，在工程计算中获得了广泛的应用。 库仑势圪(厂)同时满足=Poisson方程： 1975年，Askar首次应用有限元方法计算了氢原 AV，(，·)=一4印(，) (2) 子基态及第一激发态的能量(1】。有限元方法显示出 方程(1)，(2)中的r(r)表示分子的电荷密 其高精度的特点，并已取得了初步的成果。 度： 本文尝试采用8节点等参元和迁移式子空间迭 代法来求解双原子分子的基态总能量，为更好解 P(，。)=∑i妒i(，’)’9．(，’)· 决三维量子化学问题提供一条途径。 (胛i表示分子轨道上的电子占据数) 2Hartree-fock-s ater方程的有限元模 由于双原子分子具有轴对称性，采用柱坐标 型 变换将三维问题变成二维问题，在柱坐标系下分 子轨道的波函数为： 量子力学计算的一个主要内容就是研究多 (3) 妒(，7．二．目)=≯(，7：：)exp(i”7矽) 粒子体系的薛定谔方程的近似解法。量子力学 电子问的库仑势为： 中，分子轨道理论应用的最为广泛，它是用单 K(17，二．臼)=圪(，7．二) (4) 电子波函数来近似的表达全波函数。分子轨道 将方程(3)、(4)代入到(1)、(2)中整理得： m(x，Y，z)是电子笛卡儿坐标的(x，Y，Z)的函数，其平 方I(p(x，Y，Z)l2被解释为这一电子在空间的概率分 一圭刍(节雾卜三昙(印髻)4-12生7／声+聊痧=研≯(5) 布。 应用分子轨道理论，对于含有用个电子的分子 (6) 兰(j7Z彬T)+昙(]7§_．4印(弘) 体系，Hartree—fock—slater方程形式为：