依次尝试以下方法: (1)如果力在接近收敛的区域震荡,很有可能是电子步精度不够高,尝试EDIFF=1E-6或者EDIFF = 1E-7 (2)如果力在接近收敛的区域震荡尝试IBRION=1 (3)如果初始结构不好,或者在优化的前几步结构就崩了,IBRION=2,POTIM缩小到0.2或0.1。 (4)如果在极小点附近出现震荡,可以缩小到0.05 (5)调整初始结构。 (6)如果体系结构复杂,原子多导致结构收敛慢,力收敛趋势和能量收敛趋势向下,尝试增加收敛步数,NSW=500,不建议取大于500的值。 (7)尝试IBRION=3,阻尼md的方法往往对超大体系收敛效果好。 (8)调整晶体的对称性,先使用ISYM=0优化,如果优化一开始的初猜结构具有较高对称性,那么优化过程中这个对称性多数情况下会一直保持。如果实际上体系本身势能面极小点的结构并没有那么高的对称性,则很容易造成优化不能收敛,而且即便收敛了也会有虚频。 (9)改变收敛标准:几何优化的一般EDIFFG= -0.02,过渡态EDIFFG=-0.03. 对于体系较大难以收敛的体系,见到文献里有用EDIFFG= -0.05的。但是不建议使用更大的收敛标准了。 (10)使用更精确的积分格点精度,PREC=Accurate等。
|