催化反应计算：从过渡态搜索到能垒精度的全链路方法

催化反应计算的目标，是还原催化剂表面真实发生的一切。这个项目聚焦于二氧化碳在铜基催化剂上的电还原反应，需要通过密度泛函理论计算完整的反应路径，确定速控步骤和活化能垒。项目组很快发现，从吸附构型到过渡态再到产物脱附，每一步都暗藏精度陷阱。

催化反应计算的第一步，是搭建可靠的反应路径。项目组采用CI-NEB（攀爬图像轻推弹性带）方法搜索CO2还原为CO的过渡态，初始路径设置7个中间图像。在Cu(111)表面上，CO2吸附能在PBE泛函下计算为-0.38 eV，与实验值-0.42 eV的偏差控制在0.04 eV以内，这个精度让项目组对后续计算有了信心。

然而，过渡态搜索的过程并不顺利。CI-NEB在收敛过程中反复出现图像点回弹现象，弹簧力与真实力之间的耦合导致路径在势能面上振荡。系统在这里卡了将近十天，项目组逐一排查参数后认定，弹簧常数设置为1.0 eV/Ų偏大，将其降至0.5 eV/Ų后，路径在187步内收敛。最终的活化能垒为0.89 eV，与Nørskov课题组在Nature Catalysis上报道的0.85 eV高度吻合。

差距不会说谎。项目组同时用LST/QST方法对同一反应步骤进行了独立验证，得到的过渡态构型与CI-NEB结果的RMSD仅为0.03 Å，但能垒值偏低0.07 eV。这种差异提示我们，LST/QST对复杂反应路径的描述能力有限，在多步反应中容易陷入局部极小值。

催化反应计算中，泛函的选择直接决定结果的物理合理性。项目组在对比PBE、RPBE和B3LYP三种泛函时发现，对于CO在Cu表面的吸附能，PBE给出-0.76 eV，RPBE给出-0.52 eV，而实验值为-0.63 eV。PBE存在过强结合的系统性偏差，RPBE则矫枉过正。项目组最终认定PBE-D3更加适合这个体系——加入Grimme色散校正后，吸附能修正至-0.61 eV，偏差缩小到3%。

色散校正在催化反应计算中的价值，在含大π共轭体系的反应中更为突出。项目组在计算苯分子在Pt(111)上的吸附时，纯PBE给出的吸附能为-0.95 eV，而PBE-D3给出-1.72 eV，后者与实验值-1.68 eV的吻合度显著优于前者。原因在于苯与金属表面的作用中，范德华贡献占据了总吸附能的40%以上，纯GGA泛函根本无法捕捉这部分能量。

活化能垒只是起点，催化反应计算真正有价值的一步，是将DFT结果转化为宏观动力学预测。项目组采用微观动力学模型，将各基元步骤的速率常数通过过渡态理论计算，耦合求解稳态覆盖度方程。在温度为298 K、CO2分压为1 atm的条件下，模型预测的CO生成速率为2.3×10⁻⁸ mol/cm²/s，与流动池实验测量值的2.1×10⁻⁸ mol/cm²/s在同一数量级。

这个结果令项目组兴奋，但也暴露了一个深层问题：模型中氢析出竞争反应的速率被高估了近一倍。根源在于DFT对水分子的吸附描述不够精确，PBE泛函下水分子的吸附能偏差达到0.12 eV，这一误差通过Boltzmann因子在动力学中被指数级放大。

催化反应计算必须正视自身的局限。项目组使用的平板模型截断于4层原子，底层2层固定，这种处理在描述表面重构效应时存在不足。Nørskov组的研究表明，动态表面重构可使活化能垒变化0.1-0.2 eV，而固定模型无法捕捉这一效应。此外，溶剂效应在本项目中仅通过隐式溶剂模型VASPsol近似处理，对于涉及质子耦合电子转移的步骤，显式水分子层的引入可能使能垒进一步降低0.15-0.30 eV。

回过头看，催化反应计算最容易被忽视的环节，恰恰是计算结果与实验数据的交叉验证。项目组在初期完全依赖DFT预测的反应路径设计实验验证方案，走了一段弯路后才意识到，缺少原位光谱数据的约束，计算结果的置信度无法独立建立。被证明有效的工作流程是：DFT计算与实验表征同步推进，用原位红外光谱数据约束反应中间体的吸附构型，再用计算预测的振动频率与实验光谱比对。这套方法在后续的氮还原反应研究中将机理判断的准确率提升到了90%以上。