计算化学公社

标题: NEB计算结果远低于文献报道值的可能原因？ [打印本页]

作者
Author: Daisy_zju 时间: 2022-5-31 16:16
标题: NEB计算结果远低于文献报道值的可能原因？
本帖最后由 Daisy_zju 于 2022-6-7 14:55 编辑

求助各位老师，最近我在使用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
(, 下载次数 Times of downloads: 2)

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和文献描述一致。

(, 下载次数 Times of downloads: 3)

(, 下载次数 Times of downloads: 2)

(, 下载次数 Times of downloads: 2)

5、下图是计算结果图，能垒为0.06 eV
(, 下载次数 Times of downloads: 3)
原子最大受力的变化很正常,收敛标准设置为-0.03 eV/AA
(, 下载次数 Times of downloads: 2)
过渡态（03）的切向力也很正常，在03处由负转正且接近0
(, 下载次数 Times of downloads: 4)
验证过渡态只有一个虚频
(, 下载次数 Times of downloads: 4)

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

作者
Author: yuexiatx 时间: 2022-6-1 14:37
vtst有问题吧

作者
Author: Daisy_zju 时间: 2022-6-1 21:21

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

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

作者
Author: Daisy_zju 时间: 2022-6-1 21:22

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

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

作者
Author: Daisy_zju 时间: 2022-6-7 14:58

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

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

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

作者
Author: Daisy_zju 时间: 2022-6-10 15:42

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

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

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