CINEB路径正确，反应物和产物的能量高于过渡态

red

金属离子在2D材料表面的扩散势垒A-B-A。
0         0.000000      -120.040600         0.000000
   1         0.028548      -120.607400        -0.566800
   2         0.038461      -120.572200        -0.531600
   3         0.029367      -120.607400        -0.566800
   4         0.004664      -120.642300        -0.601700
   5         0.029008      -120.607400        -0.566800
   6         0.037808      -120.572200        -0.531600
   7         0.028449      -120.607400        -0.566800
   8         0.000000      -120.040400         0.000200
中间出现负的势阱的情况。我之前5个原子时候得到了扩散势垒是正确的。
分别计算单原子在各个吸附点位的能量也都高于1号和8号点，但是CINEB算出来的中间态的能量就比1号和8号更低。
请教各位老师这是什么情况？
计算参数
SYSTEM=NEB
LREAL= Auto
IVDW=10
ISTART= 0
ICHARG=2
LWAVE= .FALSE.
LCHARG=.FALSE.
MAGMOM= 4*1 13*0
###CINEB###
LCHAIN=0
LCLIMB=.T.
POTIM=0
IOPT=1
SPRING=-5
IMAGES=7

###ELE###
ENCUT= 520
NELM = 200
ALGO=Fast
PREC=Accurate
ISMEAR=0
SIGMA=0.1
EDIFF=1E-6
ISYM=1

IBRION=3
EDIFFG= -0.05
ISIF= 2
NSW=500
ISPIN=2

red

另外，检查了路径，路径并没有发生任何偏移。

波云见天

你是计算neb之前比较之前，有没有计算ab两个位点他们各自吸附能的高低。

red

波云见天 发表于 2021-10-27 08:49
你是计算neb之前比较之前，有没有计算ab两个位点他们各自吸附能的高低。

有的。我的这个路径能量单独拿出来构型计算能量是A=B的，也就是说1号和8号和中间的4号能量相等。

波云见天

检查一下neb中的结构，outcar中的压强是否在同一水平

red

波云见天 发表于 2021-10-27 11:44
检查一下neb中的结构，outcar中的压强是否在同一水平

请问老师，是每个结构中的external presssure这一项么？
我提取了出来，结果如下
00 -0.29 kB
01 -1.11 kB
02 -1.14 kB
03 -1.14 kB
04 -1.11 kB
05 -0.65 kB
06 -0.23 kB
请问老师这种情况正常吗？

red

波云见天 发表于 2021-10-27 11:44
检查一下neb中的结构，outcar中的压强是否在同一水平

感谢老师的解答， 我检查了正确的NEB发现他们的压强应该为正数，才能得到对的结果。请教老师如何解决这个问题呢？

yyy

你好，请问之后是怎么解决的呢，我也出现了这个问题