从文献中找到一个U值，填进INCAR，计算也顺利收敛了，是否就能继续分析结果？还需要看文献采用哪种DFT+U方法、U加在哪组轨道上，以及结构和磁态是否与自己的模型一致。尤其是文献中的U与Ueff，含义并不总相同。设置时先把这些条件对应起来，再检查带隙、磁矩和结构变化，才有依据判断参数是否合适。

　　一、VASP怎么设置DFT+U参数

　　DFT+U通过INCAR标签设置，没有统一的“推荐U值”适用于所有含同种元素的材料。下面以常用的Dudarev方法为例说明。

　　1、先确定方法和参数定义

　　①在INCAR中设置LDAU=.TRUE.，开启DFT+U。

　　②设置LDAUTYPE=2，采用Dudarev方法。

　　③核实参考资料给出的是U、J，还是有效参数Ueff，再填写LDAUU和LDAUJ。

　　该方法中起作用的是Ueff=U−J。例如，需要表达Ueff为4 eV时，可以令U为4 eV、J为0。这里的数值只用于说明参数关系，不代表某种材料的推荐设置。LDAUTYPE=1采用不同形式，不能直接照搬这组对应关系。

　　2、按元素类型填写轨道和数值

　　①核对POSCAR的元素类型顺序，以及POTCAR的拼接顺序。

　　②按相同顺序填写LDAUL：d轨道用2，f轨道用3，不加U的类型用−1。

　　③为每个类型填写LDAUU、LDAUJ，数值单位为eV，不加U的位置也保留对应项。

　　例如，类型顺序为Fe、O，仅对Fe的d轨道处理时，LDAUL写为“2−1”。数组按类型填写，不是按原子总数填写；如果人为拆分了同元素的不同类型，也需要分别对应。

　　3、补齐相关计算条件

　　①对需要考虑磁性的体系，设置合适的自旋计算方式及MAGMOM初值。

　　②涉及d、f轨道时，分别检查LMAXMIX=4或6，尤其是在读取CHGCAR进行固定电荷密度计算前。

　　③对相关过渡金属或f电子体系，检查LASPH=.TRUE.等设置是否与计算方法一致。

　　LMAXMIX涉及CHGCAR中保存的在位占据信息。准备继续计算能带时，应在生成电荷密度阶段就设置妥当，不能只在后续INCAR中补上。

　　4、确认程序实际采用了哪些设置

　　①运行后检查OUTCAR中的DFT+U参数，核对方法、轨道及数值。

　　②需要分析轨道占据时，设置LDAUPRINT=1，在OUTCAR中查看在位占据矩阵。

　　③检查电子步和结构优化的收敛情况，保存输入文件及最终结果。

　　二、VASP中U值选取不合理如何评估

　　带隙偏大、磁矩异常或结构变化明显，都可以成为检查线索，但不能单凭其中一项断定U值错误。计算精度、磁态和比较对象也可能影响结论。

　　1、先把数值误差和磁态问题查清楚

　　①检查截断能、k点、电子收敛标准及展宽设置，确认目标性质已有足够的数值精度。

　　②核对优化后的结构、原子磁矩和预期磁序，避免把错误磁态的结果用于比较。

　　③在同一U值下测试有物理依据的不同磁性初态，并比较最终状态及能量。

　　DFT+U计算可能落入不同局域极小值。两次运行都收敛，仍可能得到不同磁态或轨道占据。此时先查初态和最终解，比立即增减U更有针对性。

　　2、在有依据的区间内比较性质

　　①根据同体系文献或理论计算确定候选区间，为不同U值建立独立目录。

　　②保持泛函、赝势及数值精度一致，先在共同结构上比较电子性质，区分直接的U效应与结构变化。

　　③如果研究平衡结构，再按一致的优化标准分别弛豫，比较晶格、键长及磁矩。

　　④将带隙、投影态密度和相关结构指标一起整理，观察结论是否随U明显改变。

　　带隙接近实验，只能说明这一指标吻合，还要看带边轨道成分、磁性及结构是否合理。如果为了匹配带隙而使其他关键性质明显偏离，就需要重新审视参数或方法。

　　3、比较能量时不要跨U挑最低值

　　①比较不同磁态或结构时，在每一套固定的U、J及其他计算条件下计算相应能量差。

　　②分别观察各组能量差随U的变化，检查稳定性排序是否改变。

　　不能把U为3 eV和5 eV时的总能量直接放在一起，选择数值更低的一组作为“最佳U”。不同参数改变了能量泛函，这样的比较没有选参意义。如果相序在候选区间内反转，应说明结论对U敏感，并补充选参依据。

　　三、VASP的U值从哪里获得更有依据

　　对照文献时，元素名称相同只是起点。氧化态、配位环境、所用泛函和局域轨道定义，都关系到参数能否迁移。

　　1、把文献中的计算条件一并核对

　　①查明结构、价态、磁序、泛函、赝势及DFT+U形式。

　　②确认作者报告的是U还是Ueff，以及参数通过什么性质或方法确定。

　　③将文献值作为候选起点，用自己的模型检验相关性质。

　　同一种材料的体相参数，移到表面、缺陷或吸附模型后，也需要检查适用性。局域环境已经变化时，引用原数值应同时交代它的来源和使用范围。

　　2、必要时计算U，而非只靠试值

　　①缺少可靠参考时，可考虑线性响应或约束随机相近似等方法。

　　②采用VASP线性响应流程时，按相应步骤建立扰动、计算占据响应，并检查超胞及扰动范围的收敛。

　　LDAUTYPE=3用于线性响应计算，此时LDAUU、LDAUJ表示施加给不同自旋通道的扰动势强度，不能当成普通DFT+U中的最终U、J直接使用。

　　总结

　　选U时，需要回答的是它能否合理描述研究对象，而不只是某次计算是否收敛、带隙是否接近实验。先核对标签含义和类型顺序，再检查磁态、结构与电子性质，可以避免不少看似属于U值、实际来自输入或收敛的问题。若研究结论随U明显变化，也应把这种敏感性如实保留。希望以上说明能帮助大家建立更清楚的选参依据，让后续分析有据可查。对于VASP相关需求，欢迎联系我们进一步交流。