计算化学公社

标题: 优化结构能量差距巨大，如何解决 [打印本页]

作者
Author: archer 时间: 2020-9-27 16:01
标题: 优化结构能量差距巨大，如何解决
体系为TiO2的吸附，neb计算的中间点比初末态能量低50eV
0       0.009381    -378.551500       0.000000
1       0.680975    -436.350100    -57.798600
2       1.287765    -436.153500    -57.602000
3       0.307553    -436.736700    -58.185200
4       0.497363    -436.768700    -58.217200
5       0.315934    -435.757900    -57.206400
6       0.008585    -377.732800       0.818700

把中间的结构拿去优化，得到的能量是-300多eV，和过渡态插点的计算不同。
表面的能量为-358eV，除此之外，有的吸附态的能量也优化出了-430eV，减去分子能量吸附能也是几十eV。
这个体系出现如此巨大的能量差距，可能的原因是什么，收敛到错误的电子态？
结构优化的INCAR如下:

System=TiO2
ISTART=0
ICHARG=2
ENCUT=400
EDIFF=1E-6
IBRION=2
POTIM=0.5
NSW=4000
EDIFFG=-0.03
ISMEAR=0
SIGMA=0.1
ISIF=2
ISPIN=2
GGA = PE
#NCORE=16
#LREAL= Auto
#ALGO = Fast

LDIPOL = .TRUE.
IDIPOL = 3 # add dipole correction

LDAU=.TURE.
LDAUTYPE=2 # Only U-J is meaningful
LDAUL= 2 -1 -1 -1
LDAUU= 4.2 0 0 0#U-J = 4.2 for Ti
LDAUJ= 0 0 0 0
LMAXMIX =4 # d electron

作者
Author: tanwgwenjun 时间: 2022-6-29 14:39
您好，我刚刚接触CL-NEB，也遇到了相同的情况，请问您当时是怎么解决的呢？方便分享一下解决方法吗

作者
Author: archer 时间: 2022-6-29 15:26

tanwgwenjun 发表于 2022-6-29 14:39
您好，我刚刚接触CL-NEB，也遇到了相同的情况，请问您当时是怎么解决的呢？方便分享一下解决方法吗

这个我也没有解决，放弃这个了。

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