計算新手一枚，最近用Vasp計算析氫反應(yīng)HER的吸附H吉普斯自由能變化△G。
使用GGA  PBE  DFT+D3, 模擬計算Pt-111對H的吸附，4X4X4 slab，固定底部兩層，使用公式：△G=E（H吸附）+ 0.24 eV。調(diào)整了很多輸入文件參數(shù)，但是，我計算的△G均在 -0.3 eV左右，這個結(jié)果和文獻報道不符合（-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 返回小木蟲查看更多