CI-NEB计算迁移能垒为负值可能是什么原因？

coldwater-fish

我计算的是单个Ni在晶体Si中的迁移能，参考了文献 https://doi.org/10.1103/PhysRevB.81.235213 和 https://doi.org/10.1063/1.4807799。Ni在Si中能量最低的位置是Tetrahedral interstitial site，所以初末态为相邻两个Tetrahedral interstitial site。我按照网上的教程，首先对初态和末态都进行了结构驰豫。
结构驰豫的INCAR为：
  SYSTEM = sini
  PREC=Accurate
  ISTART = 0
  ICHARG = 2
  ENCUT = 500
  NCORE = 4
  KPAR=2
  ALGO = Fast
  LREAL = Auto
  ISMEAR = -1; SIGMA = 0.02585093
  ISYM=0
  NSW = 1000
  IBRION = 2
  ISIF = 2
  EDIFF=1E-7
  EDIFFG=-0.001

过渡态计算的INCAR为：
  SYSTEM = sini
  PREC=Accurate
  ISTART = 0
  ICHARG = 2
  ENCUT = 500
  NCORE = 8
  KPAR=4
  ALGO = Fast
  LREAL = Auto
  ISMEAR = -1; SIGMA = 0.02585093
  ISYM=0
  NSW = 300
  IBRION = 3
  POTIM = 0
  ISIF = 2
  IOPT = 7
  ICHAIN = 0
  LCLIMB = .TRUE.
  IMAGES = 5
  SPRING = -5
  EDIFF=1E-7
  EDIFFG=-0.03


但收敛后得到的迁移势垒是负的。neb.dat结果如下：
  0     0.000000     0.000000     0.000000   0
  1     0.025110    -0.000011     0.000537   1
  2     0.595661    -0.061334     0.120075   2
  3     1.165837    -0.070324    -0.061625   3
  4     1.736184    -0.058790     0.051235   4
  5     2.304823    -0.057005    -0.127756   5
  6     2.867947     0.000001     0.000000   6


请问可能是哪里出了问题呢？

bingbingppk

您好，我出现了同样的问题，您的问题解决了吗？谢谢您！