|
|
[求助]
计算贵金属的NEB出现问题
小虫在算Ir催化剂上苯的衍生物加氢过程中,发现做NEB的过程中出现了以前没有遇到过得奇怪情况。中间几个image的能量比始末态的能量高大概20 eV。当把最靠近末态的image做scf计算时,能量又回归到和末态相近的水平。整个NEB过程很合理,可以排除这个因素。INCAR如下:
######################
# Start parameters #
######################
ISTART = 0
INIWAV = 1
ICHARG = 2
NWRITE = 2
LCHARG = T ! write out "CHGCAR" file?
GGA = RE
PREC = Normal
ISPIN = 1
#############################
# Electronic relaxation 1 #
#############################
ENCUT = 400
NELM = 700
NELMIN = 2
NELMDL = -1
EDIFF = 1E-05
EDIFFG = -5E-02
######################
# Ionic relaxation #
######################
NSW = 40
IBRION = 3
POTIM = 0
ISIF = 0
ISYM = 0
ISMEAR = 0
SIGMA = 0.1
LREAL = AUTO
IALGO = 38
#####################
######neb######
######################
IMAGES = 5
SPRING = -5
EFIRST = -936.916918
ELAST = -936.642450
LCLIMB = TRUE
IOPT = 2
NPAR=2
NEB如下图所示:
请大家帮忙看看到底哪里出错误了。。。
![计算贵金属的NEB出现问题]()
NEB.jpg |
|