北京大学学报(自然科学版)
,第
39
卷,第
2
期,
2003
年
3
月
Ac
ta
Scientiarurn Naturaliurn
Universitatis
Pekinens
凹,
Vo1
.39 ,
No.2
(陆
r
,
2003)
近似相对论与非相对论密度泛函理论
自精度计算方法和程序。
王繁胡向前洪功义
2)
黎乐民
3)
(北京大学化学与分子工程学院,稀土材料化学及应用国家重点实验室,北京,
100871
)
摘
要
在己有的高精度非相对论密度泛函理论计算程序中编入新的方法进行近似相对论密度
泛函计算、产生对称性群轨道和实现解析能量梯度计算及几何构型自动优化功能的程序。改进了
原有程序的数值积分方案并且在
MPICH
环境下实现了程序的并行化。程序能满足对各种体系,特
别是含重元素的体系进行量子化学计算的一般要求。对扩展后程序的功能和结构做了简要介绍。
关键词
密度泛函理论,相对论效应,解析能量梯度,计算程序
中图分类号
0641
。寻|
~,~
仁
1
相对论效应对含重元素体系的化学和物理性质有显著影响
[IJ
在对这类体系的理论研究
中必须予以考虑。直接精确求解Di
rac
方程的计算量很大目前只能对很小的体系进行。Di
rac
算符的规则近似展开(
regular
approximation
expansion)
方法
[2
J
被认为是目前精度高的近似相对论
计算方法之一,得到的方程可以用变分法处理,而且算符形式简单,容易实现自治叠代计算,计
算量较小,结果精度满足一般化学理论研究的要求
[3
,
4
J
。
在基于
Koh
n-
Sham
方程的密度泛函理论
[5J
中交换能和相关能用电子密度的泛函近似表
示。虽然泛函的形式通常比较复杂需要用数值积分方法计算矩阵元但与
Ha
此
ree-
Fock
方法
相比计算量较小而计算精度更高。因此相关的计算方法被广泛用于各种体系的理论研究中。
作者所在科研组己编写了高精度非相对论密度泛函计算程序
[6
,
7
J
该程序的主要特点是:
( 1
)可以采用原子数值基、
STO
基或者两者混合使用;
(2)
所有矩阵元都用高精度的数值积分
方法计算;
(3)
采用电荷密度分割和多极展开的方法求解
Poisson
方程计算库仑势,可以达到
需要的任何精度;
(
4)
采用过渡态方法
[8J
计算体系的总能量,减少了数值计算误差。具体细节
见参考文献
[6
,
7]
。
在近似相对论计算方法的研究中作者提出过
ZORA
(Ze
roth-
order
regular
approximation)
方程
的简化求解方法
[9
,
10
J
、空间限域势函数
(SLF)
ZORA
方法
[IIJ
和排除奇点的相对论功能算符近似
1)
国家自然科学基金资助项目
(29892163
,29928002)
2)
现在地址
llipartment
of
Chemistry
,
University
of
Southern
California
,
L
四An
geles
,
CA
90089-1062
3)
通讯联系人
收稿日期
200
去。
1-14