过渡金属氧化物、稀土化合物、还有一些3d金属的硫族化合物——这些体系统统有一个共同特征：PBE算出来的带隙比实验值低得离谱，有时候甚至是金属态的预测结果而实验上是绝缘体。NiO是经典案例：实验带隙约4.3 eV，PBE只给出约0.8 eV，有的低对称结构甚至预测为金属。
问题出在d电子和f电子的强关联效应上。LDA和GGA对交换关联势的近似建立在均匀电子气基础上，在d/f轨道高度局域化的体系中系统性失败。DFT+U是修复这个问题的最经济方案——它不改变DFT的基本框架，只在特定轨道上施加一个类Hubbard的现场库仑修正。
U（Hubbard U）描述的是同一个原子上两个自旋相反的d/f电子之间的有效库仑排斥能。直观地理解：如果试图把两个电子同时放进同一个d轨道，它们会感受到额外的排斥能U。这个排斥能抑制了d轨道的部分占据，驱使体系向更局域化的电子态收敛。
为什么标准DFT算不出这个效应？因为LDA/GGA的交换关联泛函依赖于局域电子密度，而d电子之间的强关联是非局域的本质——同一个原子的两个d电子之间的排斥无法用只取决于该点密度的泛函来描述。DFT+U用现场库仑项U和洪特交换项J手动补上这部分缺失的物理。
VASP支持三种DFT+U的实现方式，对应LDAUTYPE参数：
| LDAUTYPE | 方案 | 双计数修正 | 适用场景 |
|---|---|---|---|
| 2 | Dudarev | 球对称平均 | 最常用，参数少（只需Ueff=U-J） |
| 4 | Liechtenstein | 全矩阵 | U和J独立可调，适合各向异性体系 |
| 1 | 简化旋转不变 | — | 旧版遗留，不推荐 |
Dudarev方案（LDAUTYPE=2）几乎成了标配。它只用一个有效参数Ueff = U – J，J值通常固定在0.5-1.0 eV（取决于元素），在INCAR中只需要设置LDAUU和LDAUJ。这种简洁性让它在高通量计算中特别受欢迎——调一个参数比调两个简单得多。
Liechtenstein方案（LDAUTYPE=4）允许分別指定U和J，对轨道各向异性（比如某些八面体配位中dxy和dz²轨道对U的响应不同）提供更精细的控制。但它的数值稳定性比Dudarev差——初值越远离收敛解，自洽迭代越容易发散。常规做法是用LDAUTYPE=2收敛后再切换到LDAUTYPE=4做精算。
这不是一个可以在计算中自洽决定的值。U是半经验参数，需要外部输入来确定。目前主流的确定方法有四条路线：
线性响应法（Cococcioni-de Gironcoli）。 在超胞中对特定原子施加微扰势，计算d轨道占据数对微扰的线性响应，反推U值。这是从头算最严谨的方法，但计算量大——需要一系列受约束的DFT计算。VASP从5.2版本开始集成了这个功能（LDAUTYPE=3），输出结果可以直接用作U值。
文献经验值。 对于常见过渡金属氧化物，大量文献已经标定了Ueff的推荐范围。这些是集合了大量实验对比后的经验值，在实际项目中比从头算的U值更可靠：
| 元素 | 推荐Ueff (eV) | 典型体系 |
|---|---|---|
| Ti (Ti³⁺) | 3.0-4.0 | TiO₂, SrTiO₃ |
| V | 3.0-3.5 | VO₂, V₂O₅ |
| Cr | 3.0-3.5 | Cr₂O₃ |
| Mn | 3.5-4.5 | MnO, LaMnO₃ |
| Fe | 4.0-5.0 | Fe₂O₃, Fe₃O₄ |
| Co | 3.0-4.0 | CoO, LiCoO₂ |
| Ni | 5.0-6.0 | NiO, LiNiO₂ |
| Cu | 4.0-6.0 | CuO, Cu₂O |
| Ce (4f) | 4.0-6.0 | CeO₂ |
| Eu (4f) | 6.0-8.0 | EuO |
拟合实验带隙。 调整U值直到计算带隙与实验带隙匹配。这个做法的争议很大——用实验结果”标定”计算参数等于放弃了预测性。但在材料筛选中，有时只需要正确排序不同材料的带隙大小而非预测绝对值，这样做的实用性高于理论洁癖。
与杂化泛函对比。 HSE06计算的带隙通常更接近实验。如果计算资源不够在全系列材料上跑HSE06，可以在少数代表性体系上跑HSE06→标定U值→在其余体系上用+U替代HSE06。这个”以HSE06校准+U”的做法在电池材料领域越来越普遍。
U值每增加1 eV，对计算结果的影响是系统性的：
带隙增大。 这是最直接的效果。U值将占据态下推、空态上推，扩大价带和导带之间的能量差距。以CoO为例，Ueff从0加到4 eV，带隙从0（金属）增到约2.5 eV。增大的幅度在不同体系中不同，取决于d轨道的局域化程度——d电子越局域，U效应越显著。
磁矩增大。 U值增强d轨道的自旋极化，驱使更多电子占据自旋向上的轨道。Fe₃O₄在U=0时净磁矩约3.5 μB/Fe₃O₄公式单位，U=4 eV时增到约4.0 μB。这个变化不是线性的，通常在U达到3-4 eV后趋于饱和。
d带中心下移。 占据的d态被U下推，d带中心向低能方向移动。这对催化计算的影响很大——d带中心是描述吸附能的关键描述符，U的取值可能改变催化活性的预测趋势。
能带色散变窄。 这是U值对电子结构最微妙的影响。更高的U使d轨道更局域化，d带宽度变窄（从PBE的约5-7 eV缩到约3-5 eV）。这个效果在角分辨光电子能谱（ARPES）中有直接可验证的实验对应。
不要对闭壳层d⁰和d¹⁰体系加U——Zn²⁺(d¹⁰)、Ti⁴⁺(d⁰)等没有部分占据d轨道的离子，U修正缺乏物理基础。这类体系的带隙问题更多来自自相互作用误差而非强关联效应，HSE06是更合适的选择。
检查OUTCAR中的d轨道占据矩阵。加U后d轨道占据数会在OUTCAR中输出，Ueff=4 eV时d轨道的总占据数变化应该小于0.1 e——如果变化显著，说明U值可能设得过大，或者体系本身不适合+U处理（比如近itinerant的3d电子）。
U对不同价态的影响不同——Mn²⁺和Mn⁴⁺在同一U值下的响应可能相差很大，因为它们d轨道占据数不同。混价体系中不同价态离子同时存在，用同一个U值可能对其中一种价态偏大或偏小——这是DFT+U方法固有的局限。
DFT+U是目前平衡计算成本和精度最实用的方案，但U值的选取仍然需要一定的经验判断。科研学术网（https://www.keyanxueshu.com）计算案例及推荐参数，覆盖电池正极、催化剂、磁性材料等多个方向。
