NEB计算结果远低于文献报道值的可能原因？

Daisy_zju

求助各位老师，最近我在使用VASP计算CINEB复现一篇文献中Li在Li2Se表面的迁移能垒（https://doi.org/10.1021/acs.nanolett.1c02658），选取了和文献一样的参数计算之后发现吸附能为 -0.21 eV,
与文献报道的 -0.46 eV 不符，随之计算的扩散能垒(0.06 eV)也远低于文献的报道值(0.63 eV)，联系文章通讯作者未果，想请各位老师帮忙查看可能是哪里出了问题，下面是我的计算详情，附件是扩散轨迹图，如果有其他需要的计算文件我后续补上来。

1、下图是我的计算案例的扩散路径 ：triple hollow Li——bridge——triple hollow Se


2、我的INCAR文件内容：
IVDW = 11
Global Parameters
ISTART =  0            (Read existing wavefunction; if there)
NWRITE = 1
ICHARG =  1            (Non-self-consistent: GGA/LDA band structures)
LREAL  = Auto           (Projection operators: automatic)
# LDIPOL = .TRUE.
# IDIPOL = 3 # add dipole correction
ENCUT  =  400          (Cut-off energy for plane wave basis set, in eV)
PREC   =  Accurate     (Precision level)  
LWAVE = .TRUE. # whether the wavefunctions are written to the WAVECAR file
LCHARG = .TRUE. # LCHARG determines whether the charge densities (files CHGCAR and CHG) are written.
LVTOT = .FALSE. # whether the total local potential is written to the LOCPOT file
LVHAR = .FALSE. # whether the electrostatic potential is written to the LOCPOT file
LELF = .FALSE. # whether to create an ELFCAR file or not.
LCHARG = .FALSE.       (Write CHGCAR or not)
ADDGRID= .TRUE.        (Increase grid; helps GGA convergence)
LREAL = Auto
# LVTOT  = .TRUE.      (Write total electrostatic potential into LOCPOT or not)
# LVHAR  = .TRUE.      (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT =             (No. of electrons: charged cells; be careful)
# LPLANE = .TRUE.      (Real space distribution; supercells)
# NPAR   = 4           (Max is no. nodes; don't set for hybrids)
NCORE  = 4

Electronic Relaxation
  ISMEAR =  1            (Gaussian smearing; metals:1)
  SIGMA  =  0.05         (Smearing value in eV; metals:0.2)
  NELM   =  60           (Max electronic SCF steps)  
  NELMIN =  4            (Min electronic SCF steps)
  EDIFF  =  1E-07        (SCF energy convergence; in eV)
  GGA  =  PE             (PBEsol exchange-correlation)

Ionic Relaxation
  NELMIN =  6            (Min electronic SCF steps)
  NSW    =  300           (Max electronic SCF steps)
  IBRION =  3            (Algorithm: 0-MD; 1-Quasi-New; 2-CG)
  ISIF   =  2            (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions)
  EDIFFG = -0.02        (Ionic convergence; eV/AA)
  POTIM = 0            (POTIM sets the time step (MD) or step width scaling (ionic relaxations)
  ISYM =  0            (Symmetry: 0=none; 2=GGA; 3=hybrids)
  ICHAIN = 0
  LCLIMB = .TRUE.
  SPRING = -5
  IMAGES = 4
  IOPT = 1
（这里我觉得ENCUT取400也有点不大合理，但是暂时还是和文献保持了一致）

3、KPOINTS选取gamma-points, 5   5   1，赝势文件分别为PAW_PBE Li_sv和PAW_PBE Se。

4、文献的computational details，我的POSCAR和文献描述一致。








5、下图是计算结果图，能垒为0.06 eV

原子最大受力的变化很正常,收敛标准设置为-0.03 eV/AA

过渡态（03）的切向力也很正常，在03处由负转正且接近0

验证过渡态只有一个虚频




实在是联系不上通讯作者，个人觉得0.06 eV的能垒的确不合理，但是各方面自查都找不到原因，还请各位老师能帮忙一起看看！感激不尽！

yuexiatx

vtst有问题吧

Daisy_zju

yuexiatx 发表于 2022-6-1 14:37
vtst有问题吧

好的，我试试重新编译一下vtst，谢谢您！

Daisy_zju

yuexiatx 发表于 2022-6-1 14:37
vtst有问题吧

不知道第二个问题，初始距离大于5A的情况前辈遇到过或者了解过吗？

Daisy_zju

yuexiatx 发表于 2022-6-1 14:37
vtst有问题吧

前辈，我重新按您的建议编译了之后可以正常计算了，但是又出现了一些问题，上面的问题已经更新了，可以麻烦您再帮我看看吗？

yuexiatx

不太了解你的体系ISMEAR=1是否合适？ 看下几个结构OSZICAR里的mag值是否相同？是否收敛到不同的电子态上？
初末态过渡态结构 取ISMEAR=0，算下单点能对比看看？
PS 我也是zju的

Daisy_zju

yuexiatx 发表于 2022-6-10 14:35
不太了解你的体系ISMEAR=1是否合适？ 看下几个结构OSZICAR里的mag值是否相同？是否收敛到不同的电子态上？
...

的确是ISMEAR不合适，我更改之后熵的值比较合理了，但是能垒还是和文献报道区别很大，感觉NEB的随机性过强了