计算化学公社

标题: 求助vasp 6.3.0优化CeO2结构电子步迭代能量无法收敛的问题，恳请路过的大神指点 [打印本页]

作者
Author: 熊啊熊 时间: 2023-7-15 14:14
标题: 求助vasp 6.3.0优化CeO2结构电子步迭代能量无法收敛的问题，恳请路过的大神指点
用新搭建的一台服务器优化一个CeO2结构，优化过程使用60个核，Material Studio上下载的cif文件，切了1 1 1表面，总共有32个Ce原子、64个O原子，计算前固定了底部三层原子。
优化过程中修改了很多参数，电子步能量始终无法收敛，并且dE降低到10的-2次方后还会跳回到10的4次方及以上，然后计算报错：Error EDDDAV: Cal1 to ZHEGy failed. Returncode = 9 1 10

下面是我的INCAR：其中AlGO试过Normal、Fast、veryfast等，ENCUT最高加到了520 eV，AMIX增加到0.2或者降低到0.05均不能收敛。

SYSTEM = opt-CeO2
NCORE = 6
ISTART = 1
ICHARG = 1
LWAVE = .TRUE.
LCHARG = .TRUE.
LVTOT = .FALSE.
LVHAR = .FALES.
#LORBIT = 11
#NEDOS = 2000

ENCUT = 450
ISMEAR = 0
SIGMA = 0.05
EDIFF = 1E-4
NELMIN = 5
NELM = 300
GGA = PE
LREAL = Auto
PREC = Normal
#LDIPOL = .TRUE.
#IDIPOL = 3

EDIFFG = -0.03
IBRION = 2
POTIM = 0.2
NSW = 300
ISYM = 0
ISIF = 3

ISPIN = 2
ALGO = Normal
MAGMOM = 32*1.5 64*0

LDAU = .TRUE.
LDAUTYPE = 2
LMAXMIX= 6
LDAUL = 3 -1
LDAUU = 5.0 0.0
LDAUJ = 0.0 0.0
AMIX = 0.1
BMIX = 0.0001
AMIX_MAG = 0.8
BMIX_MAG = 0.0001

IVDW = 12

恳请大神指路！！！

作者
Author: 淦饭青年 时间: 2023-7-15 14:36
把IVDW = 12删掉试试

作者
Author: wypkdhd 时间: 2023-7-15 18:01
简单，ALGO = Fast，注释掉NCORE = 6，使用59个核计算。

作者
Author: 熊啊熊 时间: 2023-7-16 14:34
本帖最后由熊啊熊于 2023-7-16 14:45 编辑

淦饭青年发表于 2023-7-15 14:36
把IVDW = 12删掉试试

厉害啊，虽然还是不太稳定，有几个电子步300步后仍然无法正常收敛，但是整体上可以收敛了，这是什么原理呢？
后面计算有机小分子的吸附能还是需要考虑范德华矫正，一会看看吸附小分子之后加上这个参数能不能收敛。

作者
Author: 熊啊熊 时间: 2023-7-16 14:38

wypkdhd 发表于 2023-7-15 18:01
简单，ALGO = Fast，注释掉NCORE = 6，使用59个核计算。

我试试，大佬可以说说原理吗

作者
Author: 淦饭青年 时间: 2023-7-18 18:18

熊啊熊发表于 2023-7-16 14:34
厉害啊，虽然还是不太稳定，有几个电子步300步后仍然无法正常收敛，但是整体上可以收敛了，这是什么原理 ...

我看你只优化了CeO2，没有加吸附，不用加范德华矫正，要是算吸附的话，可以加上，可以参照这个帖子http://bbs.keinsci.com/thread-15521-1-1.html

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