计算化学公社

标题: 解决VASP中SCF不收敛的方法总结 [打印本页]

作者
Author: 啦啦黑还黑 时间: 2018-10-17 17:01
标题: 解决VASP中SCF不收敛的方法总结

达到步数上限时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

作者
Author: 澳大利亚挖土豆 时间: 2023-6-26 18:47
感谢感谢

作者
Author: 柠檬鬼lemon 时间: 2023-11-8 21:05
请问我体系中只有一个原子在做单点能计算，想看一下单个原子的能量scf也不收敛是怎么回事呀，已经尝试过将ALGO=Normal，减小SIGMA都不行

作者
Author: 北大-陶豫 时间: 2023-11-9 10:31

柠檬鬼lemon 发表于 2023-11-8 21:05
请问我体系中只有一个原子在做单点能计算，想看一下单个原子的能量scf也不收敛是怎么回事呀，已经尝试过将A ...

发一下输入文件？

作者
Author: 柠檬鬼lemon 时间: 2023-11-9 11:18
本帖最后由柠檬鬼lemon 于 2023-11-9 16:59 编辑

北大-陶豫发表于 2023-11-9 10:31
发一下输入文件？

我算的是Pr （已经解决了，我加上NCORE=4就收敛了，很奇怪不知道原因是什么，最近碰上好多问题都是和核有关的）

作者
Author: wugang138 时间: 2024-3-6 15:30

柠檬鬼lemon 发表于 2023-11-9 11:18
我算的是Pr （已经解决了，我加上NCORE=4就收敛了，很奇怪不知道原因是什么，最近碰上好多问题都是和核 ...

单原子计算属于孤立系统计算，关闭对称性，破坏盒子的对称性。

作者
Author: xxxhhh 时间: 2024-5-6 13:45
想请教一下，怎么判断scf是否收敛呢

作者
Author: HydroKagami 时间: 2024-6-16 13:54

xxxhhh 发表于 2024-5-6 13:45
想请教一下，怎么判断scf是否收敛呢

看OSZICAR，假如电子步的dE一直在重复振荡（每一步都是正负差不多的值），那基本上就是不收敛了

作者
Author: phoenix 时间: 2024-7-1 16:29
我也想请问一下，我算的是一个金属表面加气体分子的体系，算一个特定结构的单点计算，电子步也是一直不收敛是为啥呀？我后面应该怎么改呀？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)