计算化学公社

标题: VASP几何优化中电子步数的设置问题 [打印本页]

作者
Author: jiangning198511 时间: 2017-3-9 14:28
标题: VASP几何优化中电子步数的设置问题
各位好，请教二个问题
VASP程序中电子步数默认为60步，并且手册中说明这个参数一般不建议修改，在几何优化(EDIFF=1e-5,EDIFFG=-0.01)过程中经常出现SCF在60步内不收敛，VASP程序会继续优化不停止，就像高斯中的IOP(5/13=1)。既然高斯中的IOP(5/13=1)应该禁止使用，那么在VASP中是否也应该在SCF收敛的条件(加大步数，加大SIGMA值)下再优化结构？还有，在最初的几步几何优化中SCF在60步内不收敛，优化若干步后SCF就能收敛了，这样条件下得到的几何结构是否合理可信？

作者
Author: mylifecg 时间: 2017-3-9 15:23
如果在60步不收敛，将IBRION=2改为1，cp CONTCAR POSCAR继续优化，如此反复，最终至结构收敛为止！

作者
Author: jiangning198511 时间: 2017-3-9 15:48

mylifecg 发表于 2017-3-9 15:23
如果在60步不收敛，将IBRION=2改为1，cp CONTCAR POSCAR继续优化，如此反复，最终至结构收敛为止！

SCF收敛性不用管了吗？

作者
Author: 卡开发发 时间: 2017-3-9 16:36

jiangning198511 发表于 2017-3-9 15:48
SCF收敛性不用管了吗？

楼上可能误认为是离子步难收敛的情形了。

一般还是先检查结构和磁性以及和截断能、k点等参数，若不收敛可以修改ALGO或者MIX的方案，尽量还是不要将SIGMA修改过大。

若60步当中出现SCF收敛，结果很可能是不可靠的，实际再进行自洽很有可能能量与离子弛豫有差别。有时候通过cp CONTCAR再计算（注意不要读取波函数和电荷密度）能起作用，但有时候即便这样做，新的结构做出的SCF还是会出现不收敛，不太提倡这样的做法。

作者
Author: jiangning198511 时间: 2017-3-9 16:46

卡开发发发表于 2017-3-9 16:36
楼上可能误认为是离子步难收敛的情形了。

一般还是先检查结构和磁性以及和截断能、k点等参数，若不收 ...

我碰到的是优化初始几步SCF不收敛，等优化若干步，SCF在60步能就收敛了，这样的结果是否可信

作者
Author: 卡开发发 时间: 2017-3-9 19:29

jiangning198511 发表于 2017-3-9 16:46
我碰到的是优化初始几步SCF不收敛，等优化若干步，SCF在60步能就收敛了，这样的结果是否可信

cp CONTCAR POSCAR之后再算一次自洽确保SCF收敛以及得到的energy、force、stress仍然落在收敛范围内。

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