密度泛函理论通常不能很好地应用于含有局域的d或f电子的强关联体系,例如,Mott绝缘体体系(包含一些3d过渡金属氧化物)。在Mott绝缘体中,当电子跃迁到另一个原子位时,若该原子位已存在电子,那跃迁中就必须克服一个库伦相互作用。如果该能量大于宽带,即使能带未填满,电子也不能完成自由运输,导致系统表现绝缘性质。对于这种情况,通常需要在密度泛函理论基础上加
DFT+U理论
入一个Hubbard模型中的原子占据位(on-site)库仑排斥项,这就是通常说的DFT+U方法,它是一种类似与Hartree-Fock平均场的方法。DFT+U方法是Anisimov等人在1991年建立的。Anisimov等人发现在传统的L(S)DA中,只包括了由Hund规则对应的交换参数J,而在Mott绝缘体体系中起决定作用的应当是Hubbard参数U,U通常比J大1个数量级。这致使了L(S)DA处理Mott绝缘体体系的失败。他们通过在原来的L(S)DA能量泛函中加入Hubbard参数U对应的一项,建立了DFT+U方法(或称为LDA+U方法)。
DFT+U,包括GGA+U,LDA+U;而GGA+U又可分为很多种,包括PBE+U, PW+U, PBEsol+U,取决于GGA的形式;+U也有很多方法,比较常见的一种就是给出U和J值,U是库伦排斥能,J是洪德耦合参数。
在VASP计算中DFT+U的设置
LDAU = .TRUE. # 控制计算中是否考虑在位库伦校正相
LDAUTYPE = 2 # +U的类型, 1|2|4; 2-默认值; U-J具有实际物理意义
LDAUL = -1 3 # 控制具体的原子轨道上加U; -1, 不加U; 1-p轨道; 2-d轨道; 3-f轨道
LDAUU = 0.0 5.5 # 电子库伦相互作用项( on-site Coulomb interaction),U值由LDAUU-LDAUJ确定
LDAUJ = 0.0 0.5 # 电子交换相互作用项( on-site exchange interaction)LDAUTYPE= 2 # +U的类型, 1|2|4; 2-默认值; U-J具有实际物理意义
LMAXMIX = 6 # 对+U体系, 设置线性混合参数, f轨道设为6
+U与未加U的带宽比较