复现关于Ti2CO2表面CO2RR文献的时候出现问题，无法复现

CalebParadoX

各位老师好，我想复现一篇文献关于Ti2CO2表面的CO2RR历程。但是我在尝试的过程中第一步就出现了瓶颈，我的第一步PCET自由能变化达到了0.52eV，文献给出的是-0.08eV，不知道是结构还是计算方法出现问题了，希望各位老师能够帮忙看一下。
文献来源：Catal. Sci. Technol., 2022, 12, 2223
文献计算细节：（见图一）
我的INCAR：
Global Parameters
ISTART =  0            (Read existing wavefunction; if there)
NPAR = 2
ISPIN =  2             (Spin polarised DFT)
ICHARG =  2            (Non-self-consistent: GGA/LDA band structures)
LREAL  = Auto          (Projection operators: automatic)
ENCUT  =  450          (Cut-off energy for plane wave basis set, in eV)
PREC   =  Fine            (Precision level)
LWAVE  = .FALSE.        (Write WAVECAR or not)
LCHARG = .TRUE.        (Write CHGCAR or not)
ADDGRID= .TRUE.        (Increase grid; helps GGA convergence)
ISYM = 0
ISIF   =  2            (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions)

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

Ionic Relaxation
NSW    =  200          (Max electronic SCF steps)
IBRION =  2            (Algorithm: 0-MD; 1-Quasi-New; 2-CG)
EDIFFG = -0.01        (Ionic convergence; eV/AA)
IVDW = 11
LAECHG=.TRUE.
我的结构：第一个PCET（附件1）
CO2吸附（附件2）
计算KPOINTS：
K-Spacing Value to Generate K-Mesh: 0.040
0
Gamma
   7   7   1
0.0  0.0  0.0

小分子能量信息：
CO2：-22.99(E)+0.31(ZPE)-0.66(TS)=-23.34eV
H2：-6.76(E)+0.29(ZPE)-0.40(TS)=6.88eV
文献路径图：（见图2）