解决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

phoenix

我也想请问一下，我算的是一个金属表面加气体分子的体系，算一个特定结构的单点计算，电子步也是一直不收敛是为啥呀？我后面应该怎么改呀？

HydroKagami

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

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

xxxhhh

想请教一下，怎么判断scf是否收敛呢

wugang138

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

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

柠檬鬼lemon

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

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

北大-陶豫

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

发一下输入文件？

柠檬鬼lemon

请问我体系中只有一个原子在做单点能计算，想看一下单个原子的能量scf也不收敛是怎么回事呀，已经尝试过将ALGO=Normal，减小SIGMA都不行

澳大利亚挖土豆

感谢感谢