第一性原理计算的理论方法
随着科技的发展,计算机性能也得到了飞速的提高,人们对物理理论的认识也更加的深入,利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手段。这主要是因为在许多情况下计算机模拟比实验更快、更省,还得意于计算机模拟可以预测一些当前实验水平难以达到的情况。然而在众多的模拟方法中,第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员的青睐,成为计算材料学的重要基础和核心计算。本章将介绍第一性原理计算的理论基础,研究方法和ABINIT软件包。
1.1 第一性原理
第一性原理计算( 简称从头计算,the abinitio calculation)
,指
从所要研究的材料的原子组分出发,运用量子力学及其它物理规律,通过自洽计算来确定指定材料的几何结构、电子结构、热力学性质和光学性质等材料物性的方法。基本思想是将多原子构成的实际体系理解成为只有电子和原子核组成的多粒子系统,运用量子力学等最基本的物理原理最大限度的对问题进行”非经验”处理。【1】第一性原理计算就只需要用到五个最基本的物理常量即( m o.e.h.c.k b ) 和元素周期表中各组分元素的电子结构,就可以合理地预测材料的许多物理性质。用第一性原理计算的晶胞大小和实验值相比误差只有几个百分点,其他性质也和实验结果比较吻合,体现了该理论的正确性。
第一性原理计算按照如下三个基本假设把问题简化:
1.利用Born-Oppenheimer 绝热近似把包含原子核和电子的多粒子问题转化为多电子问题。
2.利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单电子方程。
3.利用自洽迭代法求解单电子方程得到系统基态和其他性质。以下我将简单介绍这些第一性原理计算的理论基础和实现方法:绝热近似、密度泛函理论、局域密度近似(LDA)和广义梯度近似(GGA)、平面波及赝势方法、密度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINIT。
1.2量子力学与Born-Oppenheimer 近似固体是由原子核和核外的电子组成的,在原子核与电子之间,电子与电子之间,原子核与原子核之间都存在着相互作用。从物理学的角度来看,固体是一个多体的量子力学体系【2】,相应的体系哈密顿量可以写成如下形式:
H (r,R) E H(r ,R) (1-1)
其中r,R 分别代表所有电子坐标的集合、所有原子核坐标的集合。在不计外场作用下,体系的哈密顿量日包括体系所有粒子( 原子核和电子) 的动能和粒子之间的相互作用能,即
H H e H N H e N (1-2) 其中,以是电子部分的哈密顿量,形式为:
22 1 e2
H e(r) r2i 1 e(1-3)
i 2m 2i i,i i|r i r i |
上式的前一项代表电子的动能,后一项表示电子.电子之间的库仑相互作用
能,m是电子的质量。
原子核部分的哈密顿量H N ,可以写成:
221
H N (R) 2R j1V N(R j R j )
N j2M j R j2 j,j N j j
jj
(1-4)
原子核与电子的相互作用项可以写成:
H e N (r,R) V e N (r i r j )
i,j
(1-5)
对于这样一个多粒子体系要对其实际精确求解是非常困难的,因此对其进行简化和近似是非常的必要。考虑到电子的质量比原子核的质量小很多( 约103个数量级) ,相对来说,电子的运动速度比核的运动速度要快近千倍。当电子在做高速运动时,原子核只在平衡位置附近缓慢振动,电子能够绝热于原子核的运动。因此,可以将上面的多体问题分成两部分考虑:当考虑电子运动时,原子核要处在它们的瞬时位置上;当考虑原子核运动时,就不需要考虑不电子在空间的具体分布。这就是波恩(M.Born) 和奥本海默(J.E.Oppenheimer) 提出的绝热近似,或称波恩.奥本海默近似【2】,即Born-Oppenheimer 绝热近似。此时系统的哈密顿量简化为:
(1-6)
1.3 Hartree-Fock 轨道近似
利用 Born-Oppenheimer 绝热近似就容易把包含原子核和电子的多 粒子问题转化为多电子问题。求解方程 (1-6) 的困难在于电子与电 子之间的库伦相互作 用项。假设不考虑电子之间的相互作用,就容易得到相互独立的单
电子近似哈密顿量。为了把多电子问题简化成单电子问题【 3】, 如果把其他电子对所考虑电子的瞬时作用平均化和球对称化,则
(1-7) 这样就可以把多电子问题转变成单单子问题。这时,整个系统
的波 函数就是每个电子波函数 i (r i ) 连乘积。单电子波函数应该满足单电 子的 Hartree 方程:
H i 2
V(r i )
d r i
| i (ri )|
2m e
i (i i) i
|r i r r |
(1-8)
其中V(r) 是该电子所受到的核的作用势。 Hartree 方程描述了每个坐 标r 处单电子在核作用势和其它电子的平均势中的运动, E 是单电子 的能量,简化后就可以从假设的一组 i (r i ) 出发,求解波函数时引入 自治场方法,则整个系统的能量可以写为: E |H i (r) H i (r) E i
H
r 2
i 1
i 2m
2
i i,i
i
|r i e
2
V e N (r i R j )
r i | i,j
V i (r i )
i (i i)
d
r i
| i (r i ) |2
|r r |
(1-9) 上式并没有考虑到波函数是电子交换反对称的,于是需要考虑
尸口 础不相容原理,即把波函数写成 ( 斯莱特)Slater 行列式。此时体系 的总能要增加一个由电子交换引起的交换项,体系的总能可改写 成:
1 i (r i ) i (r i , ) i , (r i ,) i , (r i )
E | H
dr i i (r i )H i i (r i ) dr i dr i . i i i i i i i i
i 2 i,i ,
i
|r i t i , |
(1-10)
对应的单电子方程为:
2 2 | (r ) |2
(r ) (r )
[ 2 V(r i )] i (r i ) d r i
i i
i (r i ) dr i
i i i i
i
(r i ) ii i (r i ) 2m i (i i) i
|r i r r |
i (i i)
|r i r i | i
(1-11) 这就是 Hartree-Fock 方程【 4】。
2.1密度泛函的理论基础
密度泛函理论( Density Functional Theoty, 简称 DFT )【 5】是从
量子力学的基本原理出发,考虑电子结构,用体系的粒子数密度函 数替代电子波函数来描述体系的理论。 也就是说,假定固体、 原子、 分子等系统的基态能量和物理性质可以用电子密度函数唯一的确 定。密度泛函理论是由于考虑了电子相关作用的 Thomas-Fermi 模型 【6、7】,并在 Hobenberg 以及 Kohn 等人的工作【 8】后发展成的, 在经过 Kohn 和Sham(沈吕九 ) 改进得到的电子密度泛函理论中的单 电子方程,即
Kohn-Sham 方程【 9】,最终才使密度泛函理论得到实际的应用。密
度泛函理论是研究多粒子系统基态的重要方法之一,它不但成功将多电子问题转化为简单的单电子方程理论,而且也成为计算分子、固体等的电子结构和总能的有效手段。
2.2Thomas-Fermi-Dirac 近似
在1927年,H.Thomas和E.Fermi 就已经提出来建立在均匀电子气基础上的Thomas-Fermi模型【6、7】。在这个均匀的电子气模型中,电子不受外力,电子与电子之间也没有相互作用,经过求解电子运动的波动方程和简单的推导,就能看出,体系的能量仅与电子密度的函数有关。在1930年,Dirac 考虑了电子的交换相互作用并推导出来在外势V ext (r)中的电子的能量泛函的表达式如下:
E T
F (n) C1 d 3
rn(r)
53
d
3
rV ext(r)n(r) C2 d
3
rn(r)
43 1
d
3
rd
3
r
n(r)n(r )
TF 1 ext 2
2 | r r | (
2-12)
上式从左到右各项表达式分别表示:动能的局域近似、外力能作用、交换关联相互作用、经典的经典作用能。由于Thomas-Fermi-Dirac 近似太粗略简单,没有考虑到物理、化学中的一些本质现象而没用得到广泛的应用f 鲫。
2. 3 Hobenberg-Kohn 定理
密度泛函理论的基本理论基础是Hobenberg和Kohn提出的非均匀电子气理论的第一、第二定理。
第一定理:处于外势V ext ( r )中的不计自旋的电子体系,不可能存在另外一个外势V ext (r )也有相同的密度函数,即其外势V ext (r ) 可由电子
密度唯一决定。此时系统的哈密顿量H=T+V+U,这里T表示电子动能,
V是外势,U为电子相互作用势。在不同体系的哈密顿量H中,外势V
是不一样的,而电子动能T和电子相互作用势U的表达式是相同的。
因此只要外势确定,体系的哈密顿量H也就确定了。根据公式
H E ,只要H是确定的,系统的波函数也确定,也可以说电子密度决定了系统波函数的所有性质。
第二定理:对于已定的外势,体系基态能量能于基态能量泛函
E(n(r)) 的极小值。对于不计自旋的全同电子体系,其能量泛函
E(n(r)) 可写为:
2
E(n (r)) V(r)n (r)dr T[n (r)] e Cn (r)drdr E xc[n(r)]
2 |r r |
(2-13) 其中,第一项是电子在外势场中的势能,第二项表示无相互作用电子气的动能,
第三项是电子间的库伦作用能,第四项是电子间的交换关联能。第二定理的基本点是在粒子数不变条件下求能量对密度函数的变分,就可以得到体系基态的能量E(n) 。但是Hobenberg-Kohn定理中还存在一些不足之处:
(1)电子密度分布函数n (r) 的具体形式不明确。
(2)无相互作用电子气的动能泛函T[ n (r) ]不知道。
(3)电子间的交换关联能泛函E xc[n(r)] 不清楚。针对前两个问题可以用Kohn-Sham方程解决。第三个问题,通常是采用各种近似得到电子间的交换关联能。·
2. 4有效单电子近似: Kohn-Sham 方程
1965年, Kohn 和Sham 提出了这样一个假设:体系的电荷密度可以
用 电子波函数构造。此时电荷密度
N
2
n(r) | i (r)|2
i1
(2-14)
这样前面遇到的问题就可以顺利解决。 将 i (r) 代到(2 .13)变形成; E[n(r)] T o [n(r)] n(r)V ext (r)dr E h [n(r)] E xc [n(r)] ( 2-15 )
2 N
2
其中
T o [n(r)]
i 2 i
,
2m e i 1
( 2-16 )
E [n(r)]
1 n(r)n(r )
drdr
E h [n(r)]
drdr
2 r r
(2-17)
虽然 E xc [n(r)] 与电子密度 n(r) 之间的函数表达式不知道,但是
Kohn 和Sham 成功的将多电子体系的薛定谔方程问题简单的归结为单
电 子在周期性势场中的运动的单电子方程。此时,只要求解在周期性 势场N 个无相互作用的单电子方程:
2
[ 2
V KS [n(r )] i (r) i i (r) 2m
2-19)
(2-18)
V KS V[n(r)]
E h [n(r)] n(r ) E XC [ n( r )
]
根据Kohn-Sham 的本征值 i ,体系的总能量可写成:
E i
1 n(r)n(r )
drdr V XC [n(r)]n(r)dr E XC [n(r)] i 2 r r
(2-20)
需要注意的是 Kohn-Sham 方程中本征值没有实际的物理意义。唯一 的例外是体系的最高占据轨道,它的本征值对应于体系的离子化能 【10】。
2.5 交换关联能近似
电子间的交换关联能泛函 E xc [n(r)] 表示的是所有其它多体项对总能 的贡献。它的物理意思是:当单电子在一个多电子体系运动中,由 于考虑电子之间的库伦排斥,电子与体系之间就有交互关联作用。 换句话说,就是在同一时刻两个电子不可能占据同一个位置,也就 产生了交换关联能 E xc [n(r)] 。在HoBenerg-Kohn-Sham 的理论框架下, 多电子体系基态的薛定谔方程问题转化成了有效的单电子方程问 题,这种形式的描述比胁舰 P 粕出方程更严密更简洁。但前提是要 处理好交换关联能后这个理论才有实际的应用价值。所以交换关联 能泛函在密度泛函理论中占有非常重要的地位。
2. 6 局域密度近似 (LDA)
1965 年 Kohn 和 Sham 所 提 出 了 局 域 密 度 近 似 (Local Density Approximation) 【 1l 】。局域密度近似的主要原理是
假设非均匀电 子体系的电荷密度的变化是相当的缓慢,可以将这个体系分成很多 很多个足够小的体积元,近似的认为每个小体积元中的电荷密度是 一个常数刀 n(r) ,则在这样一个小体积元中的电子气分布
是均匀的 并且没有相互作用,而对于整个非均匀的电子体系总体来说,各个 小体积元
的电荷密度只与它所处的空间位置 r 有关。因此,交换关联能可以 写成如下形式:
E x L
c DA
n(r) xc (n(r))dr
(2-21)
对应的交换关联势写为:
LDA
E x L
c DA [n]
xc
[n]
V xc LDA
[ n(r )] xc
xc
[n] n xc
nn
(2-22)
其中 xc (n) 特指均匀电子气中的交换关联能密度。 交换关联近似的形式多种多样,目前在 LDA 自洽从头算中用得最多 的交换关联势是
Ceperley-L.Alder 交换关联势,它是采用目前最精 确 的 量 子 Monte-Carlo 方 法 计 算 均 匀 电 子 气 的 结 果 , 并 由 T.P.Perdew 和A,zunger 参数化得到的交换关联函数。一般分为交
换 和关联两个部分:
xc
[n] x [n] c [n]
(2-23)
由 Dirac 给出的交换能可写为:
x
[n] C x n(r )
1 3
(2-24)
这里
13
3 3
13
4
(2-25)
关联能的精确值最早由D.M.Ceperley 和B.L.Alder 通过量子Monte-Carlo 方法计算获得【12】。而xc(n) 由
T.P.Perder 和A.Zunger
参数【13】得到。交换能表达式如下:
(2-26)
关联能形式如下:
c[r s] 0.2846 (1 1.0529 r s 0.3334r s
0.0960 0.0622 ln r s 0.00232r s 0.004r s
(2-27)
这里Weigner-seitz 半径,在均匀电子气模型中,表达式为:
31 3 r
s
4 n(r )
(2-28)
对于价电子r 的值通常是1~6之间;对于芯电子而言r s通常是小于l
LDA近似一般适用于电子密度变化比较平缓的体系,对于一些强关联系统如过渡金属和稀土金属等缺陷是很明显的。因此,需要对其进行一些适当的改进和修正。这就使得各种广义梯度近似(GGA)得到了发展的空间。
2.7 广义梯度近似(GGA)
0.9164
广义梯度近似就是在局域密度近似的基础上考虑了电荷密度的梯 度,换个说法是: 交换关联能密度不仅仅和该体积元内的局域电 荷密度有联系,还跟邻近小体积元的电荷密度有关,这时就要考虑 这个空间电荷密度的变化,考虑到电荷密度分布的不均匀性,就要 引入电荷密度梯度。此时
E xc [n] n(r) xc (n(r )) dr E x G
c GA
( n(r ) | n(r)|)
(2-29)
近年来发展起来的广义梯度近似 (GGA)已经有很多中样式,比较常 见的交换关联能有 Perdew-Wang(PW9【1) 14】Perdew-Burke-
Emerhof 【PBE)【15】和 BECKE8【8 16】
需要说明的是: GGA 和 LDA 两种交换关联能近似没有孰优孰劣之分, 只能由实际计算的体系来判定。
参考文献
【5l 】吴兴惠,项金钟.现代材料计算与设计教程.北京:电子工业 出版社, 2002,
p .173.
【 52】 Bom MHuang KDynamical Theory of Ctrstal
Lattices.Oxford :Clarendon , l954
【53】DRHartreeProcCamPhilSoc ,24:89,1928. 【54】v.Fock .Phys . Rev.B , 75:01240l ,2007. 【55】Chelikowsky J R,Louie S GQuantumth eory of real
materials
M】.Kluwer Academy Press ,1989:1-11 .
56】Tomas Proc L H .The calculation of atomic fields
【J】.Cambridge PhiloSophy Society l927 ,23:542-545.【57】Fermi E.An method statistic par la determination diaconal
proprietary ,dell attome .Accad.Naz.Lincei ,1 927,
6:602.605.【58】H0benberg,P,Kohn W.inhomogeneous
e1etron .Gas[J].Physical Review B ,1964,l 36 :
864·871.
【59】Kohn w,W,Sham L J.self-consisent equations
including exchange and correlation effcts[J]
.Physical Rcview A ,
1965,
140:1133~A1138.
【60】阚二军,中国科技大学博士学位论文,(2008)
【61】Kohn Wand Sham L J,Quanutm Density Oscillations in an Inhomgeneous ElecnonGas,PhyS.Rev.1965,137:A1697.A1705 【62】D.M..Caperley and B.J.A1der ,Phys .Rev Lett .1 980,45:566
【63】T.P.Perdew and A.Zung.Self-Interaction
Correction to Density .Functional
Approximations for Many-Electron Systems.Phys.Rev.B.198l ,23:5048.
【64】P.Perdew,J.A.Chevary ,S .H.Vosko,K.A.Jackson ,M,R,PederSon,D.J.Singh,andCFiomais.Atoms,mo1ecules,so1ids and surfaces :Applications of the generalizedgradient approximation for exchange and correlation .Phys.ReV.B.1 992,46:6671.
【65】J.P.Perdew,K.Burke,M.Emzerhol£GeneraJizcd GradieIlt Approximation
Made Simple.PhyS.Rev.Lett .1996,77:3865.
【66】A.D.Becke.Phys.Rev.A,l988 ,38:3098.