请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 677|回复: 0

[VASP] 解决VASP中SCF不收敛的方法总结

[复制链接]

19

帖子

0

威望

233

eV
积分
252

Level 3 能力者

发表于 2018-10-17 17:01:49 | 显示全部楼层 |阅读模式
达到步数上限时SCF仍没有收敛是实际研究中极为容易碰见的问题,VASP中经常会碰到电子不收敛的情况,也叫做SCF不收敛。经常会在电子结构比较复杂的体系中碰到。以下是我提供的可能解决方案,有不妥的地方,欢迎讨论。同时欢迎加入理论计算化学研习社qq群:905699144



依次尝试以下方法:
(1)先检查几何结构是不是合理,非常离谱的初始结构会导致SCF难收敛。
(2)依次检查已经设置的参数,对照ppt的建议,是不是选择合理,对于+U和ICHARGE=11的任务,添加LMAXMIN=4 (对于d区体系),添加LMAXMIN=6 (对于f区体系)
(3)检查是不是ISTART=1读取了不合理的波函数,如果是,rm WAVECAR CHGCAR重新跑。
(4)如果ALGO=Fast或者VeryFast,换成Normal。
(5)尝试使用更大的SIGMA值,先粗略收敛,再读取CHGCAR和WAVECAR用小的SIGMA计算。
(5)对于非磁性体系(闭壳层ISPIN=1)添加:(注意AMIX和BMIX对收敛有很大影响,可以自己调试)
AMIX = 0.2
BMIX = 0.0001 #almost zero, but 0 will crash some versions
(6)对于磁性体系(自旋极化,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
(7)  尝试更换不同的ISMEAR。
(8)  检查体系是不是特殊的磁性排列,即MAGMOM设置是否合理。
(9)  提高积分精度,PREC=Accurate 。
(10) 提高格点精度,ADDGRID = .TRUE. 。
(11) 先用1 1 1 K点计算收敛,再读取CHGCAR,用高K点计算。
(12) 尝试不同的ALGO 比如:ALGO=Conjugate。
(13) 如果在结构优化或者MD过程中,某一步突然不收敛,使用MAXMIX = 50
(14) 尝试用更小的ENCUT或者更大的ENCUT的预收敛。
(15) 换更小的赝势或者更soft的赝势。

(16) 最后给出一个VASP官方教程里解决的不收敛的方法:
1.用ALGO=N (是否收敛N to 2,Y to 6)
2. ICHARG=12 (no charge update,非自洽计算,N to5, Y to 3)
3. ICHARG=2 AMIX=0.1 ; BMIX=0.01(N to 4, Y to 6)
4. increase BMIXBMIX=3.0 ; AMIN=0.01(N to 5, Y to 6)
5. Bug report

6. Use this setting

评分

参与人数 6eV +30 收起 理由
lsing + 5 赞!
Marisa + 5 谢谢
obaica + 5 谢谢
wuhuawujiu + 5 赞!
archer + 5 赞!
Lanthanide + 5

查看全部评分

您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-11-17 05:13 , Processed in 0.103567 second(s), 25 queries .

快速回复 返回顶部 返回列表