纯干货：DFT+U | VASP零基础实用教程

密度泛函理论通常不能很好地应用于含有局域的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的带宽比较