3.3.5 结构优化计算不收敛怎么办#
对于原子多的大体系,常出现的问题是结构优化不收敛。此时要分清楚两种情况:
(1)某一离子步中的静态计算不收敛,使得原子受力和压强不正常,最终结构演化到奇葩的情况。该情况,请参照“自洽计算不收敛”问题,先解决自洽计算不收敛。
(2)(静态计算收敛了)原子受力始终没法降低到EDIFFG以下。有时候查看relax过程中的原子受力,会发现达到最大离子步NSW时,原子受力仍然在0.1附近震荡,没法降低到0.01(假如EDIFFG=-0.01)。此时可能是INCAR中POTIM太大(默认0.5),体系一直在局域极小值附近震荡没法达到极小值,此时需要减小POTIM,比如改成0.2或者0.1,再复制CONTCAR为POSCAR,继续优化。或者提高EDIFF的精度,从1E-6提高到1E-8;或者修改IBRION算法,1和2都尝试,个人经验,2比1更好。或者网上搜索“结构优化不收敛”的相关经验。