解决VASP中SCF不收敛的方法总结
依次尝试以下方法:
- 先检查几何结构是不是合理,非常离谱的初始结构会导致SCF难收敛。
- 依次检查已经设置的参数,对照ppt的建议,是不是选择合理,对于+U和ICHARGE=11的任务,添加LMAXMIN=4 (对于d区体系),添加LMAXMIN=6 (对于f区体系)。
- 检查是不是ISTART=1读取了不合理的波函数,如果是,rm WAVECAR CHGCAR重新跑。
- 如果ALGO=Fast或者VeryFast,换成Normal。
- 尝试使用更大的SIGMA值,先粗略收敛,再读取CHGCAR和WAVECAR用小的SIGMA计算。
- 对于非磁性体系(闭壳层ISPIN=1)添加:(注意AMIX和BMIX对收敛有很大影响,可以自己调试)
AMIX = 0.2
BMIX = 0.0001 #almost zero, but 0 will crash some versions - 对于磁性体系(自旋极化,ISPIN=2)添加:
AMIX = 0.2
BMIX = 0.0001 #almost zero, but 0 will crash some versions
AMIX MAG = 0.8
BMIX MAG = 0.0001 #almost zero, but 0 will crash some versions - 尝试更换不同的ISMEAR。
- 检查体系是不是特殊的磁性排列,即MAGMOM设置是否合理。
- 提高积分精度,PREC = Accurate 。
- 提高格点精度,ADDGRID = .TRUE. 。
- 先用1 1 1 K点计算收敛,再读取CHGCAR,用高K点计算。
- 尝试不同的ALGO 比如:ALGO = Conjugate。
- 如果在结构优化或者MD过程中,某一步突然不收敛,使用MAXMIX = 50
- 尝试用更小的ENCUT或者更大的ENCUT的预收敛。
- 换更小的赝势或者更soft的赝势。
最后给出一个VASP官方教程里解决的不收敛的方法:
- 用ALGO=N (是否收敛N to 2,Y to 6)
- ICHARG=12 (no charge update,非自洽计算,N to5, Y to 3)
- ICHARG=2 AMIX=0.1 ; BMIX=0.01(N to 4, Y to 6)
- increase BMIXBMIX=3.0 ; AMIN=0.01(N to 5, Y to 6)
- Bug report
- Use this setting
解决 VASP 中 SCF 不收敛的方法总结 http://bbs.keinsci.com/thread-11126-1-1.html (出处: 计算化学公社)
转载请注明来源,欢迎对文章中的引用来源进行考证,欢迎指出任何有错误或不够清晰的表达。