在上一章《第十章：探索密度泛函理论之平面波基组详解！| 2026新版VASP基础教程》中，华算科技朱老师详细介绍了密度泛函理论基础—赝势与基组。VASP通过自洽场计算方法，通过密度泛函理论求解材料的性质。本章将正式引入这一理论核心，从物理图像出发说明自洽场计算在描述不同材料体系时的关键作用。
在 VASP 中，收敛指的是计算过程中的物理量（如总能量、原子力、应力张量等）在迭代过程中趋于稳定，变化幅度小于设定的阈值。
指在固定离子位置的情况下，求解 Kohn-Sham 方程得到的电子密度和总能量在迭代中趋于稳定。实质上是通过自洽循环，使得输入密度与输出密度一致。
指在电子自洽收敛的前提下，通过优化离子位置（或晶格参数），使得体系的总能量达到局部极小值，原子力趋于零或应力趋于零。
判断 VASP 是否收敛主要通过观察 OUTCAR 文件（或 OSZICAR）中的特定行。
关键指标：Free energy = 或 E0= 行后面的数值（能量）。
判据：连续两次迭代的自由能差值（dE）小于 EDIFF，即 |ΔE|
观察点：在 OUTCAR 中搜索 reached required accuracy 或 converged，若出现此类信息，则电子收敛成功。
关键指标：FORCES: 或 TOTAL-FORCE 行下的力值。
判据：所有原子力的最大值（max）小于 |EDIFFG|（注意 EDIFFG 为负值，取绝对值比较）。
观察点：在OUTCAR 中搜索 reached required accuracy – stopping structural energy minimisation，若出现此行信息，则离子弛豫收敛成功。
EDIFF not reached：电子收敛未达标。需检查 NELM 是否足够，或调整 ALGO、NELMDL、LMIX 参数。
POTIM or IBRION warnings：通常与离子优化步长过大或算法不匹配有关。
ZBRENT errors：通常与体积优化过程中的数值问题相关。
在实际计算中，特别是对于金属、半金属、磁性体系或包含稀土元素（镧系、锕系）的体系，往往会出现难以收敛的情况。以下是针对性解决方案：
问题：费米面附近电子态密度变化剧烈，导致 SCF 振荡。
解决方案：
将 ISMEAR 设为 1（Methfessel-Paxton）或 2（Gaussian），并适当调大 SIGMA（如 0.2-0.3 eV）。
采用 ALGO=All 或 ALGO=Fast。
调整混合参数：NELMDL=-6（线性混合）或 NELMDL=-4（Pulay 混合）。
问题：自旋极化导致电子结构复杂，初始磁矩设置不当。
解决方案：
确保 ISPIN=2 开启自旋极化。
在 MAGMOM 中设置合理的初始磁矩（如 MAGMOM = 5*0.6 表示 5 个原子初始磁矩为 0.6 μB）。
调整混合参数 LMAXMIX（对于含 d、f 元素的体系尤其重要）。
适当降低 EDIFF（如 1E-5）以提升收敛精度。
问题：f 电子局域化强，电子相关性高，导致 SCF 振荡。
解决方案：
将 LMAXMIX 设为 6（对应 f 电子）。
采用 ALGO=All 或 ALGO=Fast。
对于难收敛的体系，可以考虑使用 GGA+U（在 INCAR 中添加 LDAU=.TRUE.）来纠正电子相关性。
适当增大 NELM。
问题：晶格参数波动幅度大，导致能量或应力不收敛。
解决方案：
检查 ISIF 是否设置正确（如只优化位置应设为 2，优化体积和形状应设为 3 或 6）。
调整 POTIM（通常对体积优化设为 0.1-0.2，防止步长过大导致体系爆炸）。
增加 NSW（最大离子步数），以便系统有足够的迭代次数收敛。
自洽场理论基础：介绍了自洽场计算的理论依据与关键指标
收敛标准：介绍了电子弛豫与离子弛豫的收敛标准
收敛策略：介绍了金属、磁性、澜系体系的策略
下一章将正式引入本次教程的核心—PAW赝势原理。我们将从理论基础、核心构造、优势与局限性，详细介绍PAW赝势原理，以及它在处理不同材料体系时的独特优势的应用，敬请期待！
