计算化学公社

标题: 复现关于Ti2CO2表面CO2RR文献的时候出现问题，无法复现 [打印本页]

作者
Author: CalebParadoX 时间: 5 day ago
标题: 复现关于Ti2CO2表面CO2RR文献的时候出现问题，无法复现
本帖最后由 CalebParadoX 于 2025-8-8 19:16 编辑

各位老师好，我想复现一篇文献关于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）

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