计算新手一枚，最近用Vasp计算析氢反应HER的吸附H吉普斯自由能变化△G。
使用GGA  PBE  DFT+D3, 模拟计算Pt-111对H的吸附，4X4X4 slab，固定底部两层，使用公式：△G=E（H吸附）+ 0.24 eV。调整了很多输入文件参数，但是，我计算的△G均在 -0.3 eV左右，这个结果和文献报道不符合（-0.09 eV）。这个问题头疼了很久，找不到原因，请各位虫友大牛指点迷津。谢谢。

附INCAR

INCAR-basic
   PREC   = Normal
   ENCUT  =450 eV
NELM = 300
   LREAL = Auto
   IVDW = 11
   ALGO = Fast
   EDIFF  = 1E-5
  ISMEAR = 1
   ISPIN = 1
   SIGMA = 0.2
   NPAR = 4
Ionic Relaxation
   ICHARG = 2
   EDIFFG = -0.03
   NSW   = 300   
   IBRION = 2
   POTIM = 0.5
   ISIF   = 2   
   LCHARG = .FALSE.
   LWAVE  =  F 返回小木虫查看更多