|
|
体系为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
|
|