VASP计算稀土元素不收敛怎么办？从CeO₂案例看磁性体系的自洽场策略

做过含稀土元素VASP计算的同行大概都有这个经历：跑了一天一夜，打开OSZICAR一看，电子步在10⁻³到10⁻⁴之间反复横跳，死活不收敛。这个项目在算CeO₂(111)面氧空位形成能时，就遇到了典型的”稀土魔咒”——Ce的4f电子态密度又窄又深，电子占据在f轨道之间反复转移，SCF像被按了循环键一样。

排查：真的是f电子的问题

不收敛很多原因，先排除硬件和参数层面的低级错误。K点密度、截断能、初始磁矩这些常规检查过一遍没有问题。把Ce的4f电子替换为赝核处理（即把f电子冻在芯里），SCF在18步内收敛到10⁻⁶——问题直接锁定在4f电子上。

不收敛机制是这样的：Ce³⁺（f¹）和Ce⁴⁺（f⁰）的电子构型能量接近，SCF迭代过程中电子在两种价态间来回跳跃。每一步的电荷密度和上一步差异偏大，导致Pulay mixer的混合效果变差——常规的AMIX=0.4、BMIX=1.0对这个体系来说”反应太慢”。

三招打通

第一招：切换到线性混合。INCAR里设置：

IMIX=0      # 关闭Pulay，启用线性混合
AMIX=0.05   # 大幅降低混合比例
AMIX_MAG=0.02  # 磁矩混合更保守

线性混合收敛慢但稳定，适合电荷密度剧烈震荡的场景。AMIX从默认0.4砍到0.05，相当于告诉VASP”每一步只信5%的新电荷密度”——慢是慢了，但不再震荡。

第二招：Smearing策略调整。CeO₂是绝缘体（实验带隙~3.2 eV），按理说用ISMEAR=0（Gaussian）就够了。但项目发现ISMEAR=-5（四面体方法+Blöchl修正）搭配SIGMA=0.05对这个体系收敛更友好——不是物理上更合理，而是四面体方法在窄带体系中对Fermi能附近的态密度积分更稳定，间接减少了占据数震荡。

第三招：磁矩初始化。Ce的初始磁矩设MAGMOM=1.0（对应Ce³⁺的f¹），而不是默认的0或随机值。给一个合理的初始磁矩猜测，等于告诉SCF”往这个方向走”，缩小搜索空间。项目在O 2p电子上额外加了MAGMOM=0.1，起辅助收敛作用。

三招组合下来，SCF在43步收敛到10⁻⁵，虽然比普通体系多了一倍步数，但至少不再震荡。如果需要更高精度，可以在此基础上切回Pulay混合跑到10⁻⁶——项目实际用了这个两步走策略。

事后复盘：什么体系容易出这种问题

f电子不收敛不是CeO₂的”专利”。含Pr、Nd、Sm、Eu、Gd等镧系元素的氧化物、氮化物、硫化物都可能中招。共同特征是f轨道局域化强、电子关联明显、不同价态能量接近。对这类体系，Pulay混合不是最好的选择，线性混合+超低AMIX是更稳妥的起手式。

如果你遇到的不收敛伴随着磁矩剧烈震荡（OSZICAR里mag列每步跳几百%），AMIX_MAG单独降到0.01比调AMIX更管用——磁矩震荡和电荷震荡有时是两码事，分开控制效果更好。

参考文献：Fabris et al., Phys. Rev. B, 2005, 72, 235215；Da Silva et al., Phys. Rev. B, 2007, 75, 045121.