计算化学公社

标题: CI-NEB计算迁移能垒为负值可能是什么原因？ [打印本页]

作者
Author: coldwater-fish 时间: 2023-3-9 10:24
标题: CI-NEB计算迁移能垒为负值可能是什么原因？
我计算的是单个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

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

作者
Author: bingbingppk 时间: 2024-5-5 21:38
您好，我出现了同样的问题，您的问题解决了吗？谢谢您！

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