密度泛函理论
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Local Density Approximation (LDA)
E LDA x
r
Cx
4/
3
r
dr
LDA x
r
Cx
1/
3
Local Spin Density Approximation (LSDA)
ELSDA x
21/3Cx
4 / 3
4/3
dr
LSDA x
1 2
Cx 1/3
BPW91 = B88 + PW91 B3PW91 = B3 + PW91
➢ 一般而言,GGA 比 LSDA 效果要好得多。GGA 的计算量与 HF 相仿,但构型和振 动频率的精确度一般要好于 MP2,与 CC 可比拟。DFT 方法总的来说对静电作用描 述得更好一些,而对于范德华作用描述得差一些。
➢ DFT算法的实现与 HF 相似。基本的DFT算法复杂度为 M 4 ,最新计算技术使 DFT 的计算量线性化。
2b
2x0 Q
tan1
Q 2x b
f 1 4/3 1 4/3 2 2 21/3 1
x rS X x x2 bx cQ 4c b2
GGA (见下)中的 PW91 修改了 VWN 的泛函形式:
PW91 c/a
x
2a
1 x2
ln 1
2a
1
1x 2 x2 3x3 4 x4
1
f 1 1 2/3 1 2/3 2
其中 c rS , 在 LSDA 部分已经给出。
Becke(B95):更好地满足一些基本的物理约束。
B95 c
c
c
c
c
1 a
x2 x2
1 PW 91,
c
c
1 bx2 2
D
DLDA
PW 91, c
DLDA
2 C 5/3
1
4/3 1
4/3
关联项
Vosko,Wilk,and Nusair (VWN)
VWN c
rS ,
c
rS , 0
a
rS
f f
''
0
1
4
c
rS ,1 c
rS , 0
f
4
c/a
x
A ln
x2
X x
2b Q
tan1
Q 2x b
பைடு நூலகம்
bx0
X x0
ln
x x0 2 X x
➢ DFT 的最大问题在于没有统一的理论方法系统地提高计算精度,即更复杂的泛函 形式不一定计算精度越高,而是与被研究体系密切相关。
5/3
F
➢ 混合方法
混合 HF 和 DFT 给出的能量项。
Becke 3 parameter functional (B3)
EB3 xc
1 a
E LSDA x
aExHF
bExB88
E LSDA c
cEcGGA
➢ 交换和关联项的组合应用
SVWN = LSDA + VWN BLYP = B88 + LYP BP86 = B88 + P86 B3LYP = B3 + LYP B3P86 = B3 + P86
2 2 eab ab eab a3 4b
e ab
8
a ab 2 b 2 2DD
N i
i
2
2
4
Perdew and Wang (PW91)
PW91
x
LDA x
xa1 sinh1 xa2 1 xa1 sinh1
a3 xa2
a4ebx2 a5x2
x2
x
4/3
PW c
91
H0 t, rS , H1 t, rS ,
H0
t, rS
,
b1
f
3
ln
1
t2 a 1 At2
At 4 A2t 4
H1
t,
rS
,
16
3 2
1/3 C c
f
t e 3 2 dx2 / f 2
t
192
2
1/ 6
2
f
7/6
a
exp
bc rS , / f 3
P86 c
LDA c
P86 c
P86 c
eC f 7/3
2
a
C C 7/6
f 21/3
1
2
5/3
1
2
5/3
C b1
b2 b3rS b4rS2 1 b5rS b6rS2 b7rS3
Perdew and Wang(PW91 or P91):改进 P86。
关联项
Lee, Yang, and Parr (LYP)
LYP c
a
1 d 1/3
ab 9
ec1/ 3 1 d 1/3
8/3 18
22/3
CF
8/3
8/3
18tW
2tW 2
2tW 2
2 1
2 2 2
tW
1 8
2
2
Perdew(P86):修正 LSDA 的梯度项。
虽然证明了电子密度和基态能量的一一对应关系是存在的,但是两者之间的泛函形 式未知。各种 DFT 的目的就是从不同的简化物理图象出发,给出近似的泛函形式。
E [ρ] = T [ρ] + Ene [ρ] + J [ρ] + K [ρ] Ene [ρ] (核子-电子势能)和 J [ρ] (库仑积分)已知,T [ρ](动能)和 K [ρ](交换积 分)未知。
➢ Gradient Corrected Methods
Gradient Corrected or Generalized Gradient Approximation (GGA): 泛函不仅决定 于电子密度,还决定于电子密度的梯度。
交换项
Perdew and Wang (PW86): 修正 LSDA 的泛函形式:加入高阶项。
➢ Kohn-Sham 理论
DFT 的 HF 理论。给定了未知泛函的形式后,类似 HF 方法,得到准本征态方程 Kohn-Sham 方程。HF的计算量,但是自动包括了电子关联的贡献。
➢ Local Density Methods
假设局域电子密度可以被认为是均匀电子气,或等效地说,电子密度是随空间缓慢 变化的函数。 交换项
PW 86 x
LDA x
1 ax2 bx4 cx6
1/15
x
4/3
Becke (B or B88): 正确的能量密度渐进行为。
B88 x
LDA x
B x
88
B88 x
1/3
1
6
x2 x sinh1
x
Becke and Roussel (BR): 加入轨道波函数的导数项。
BR x
密度泛函理论
➢ Hohenberg-Kohn 定理
1)基态系统的所有物理性质都由电子密度唯一决定,能量与电子密度为一一映射。
2)对应于电子密度的变分原理:任意近似电子密度所对应的能量值都大于等于基 态对应的真正密度所决定的能量值。
➢ 密度泛函理论 (Density Functional Theory, DFT)