催化还原计算

第一次做催化还原计算，我拿着一篇文献照抄参数，算出来的 CO₂ 还原限速步能垒比人家低了 0.4 eV，高兴了半天，后来发现人家用了隐式溶剂而我没加。那次让我意识到，催化还原计算不是跑通 workflow 就完事，每一个建模假设都会直接写进能垒数字里。

一、为什么催化还原计算先选对表面模型

还原反应多发生在金属或氧化物表面。我一般先用 XRD / XPS 信息锁定暴露晶面，比如 Cu 常用 (111)，TiO₂ 常用 (101)。表面层数至少 4 层、底部两固定，真空层 15 Å 以上避免镜像相互作用。曾经用一个 3 层 slab 算氧还原，底固定层应力没释放，吸附能偏差 0.2 eV。做催化还原计算时，slab 的弛豫充分性比 k 点密度更该先确认。

二、核心原理：自由能台阶图怎么搭

催化还原计算的主线是算每个基元步骤的 ΔG。以 CO₂ → CO 为例，中间体 *COOH、*CO 的吸附自由能用 E_ads = E(slab+ads) − E(slab) − E(gas) 得到，再叠加上熵修正和零点能（ZPE）。我用的是 Nørskov 那套方法：气相分子自由能查表，表面态用 DFT 总能量加 ZPE 修正。关键是所有中间体要锚定在同一个参考（如 U = 0 V 的 RHE 标度），否则台阶图会整体平移，看不出谁是限速步。

三、关键技术要点：U 值和泛函不能乱填

过渡金属氧化物离不开 DFT+U。我做 Co₃O₄ 上的还原时，Co 的 U 从 3.0 扫到 4.0 eV，带隙从 0.8 eV 变到 1.6 eV，吸附能跟着飘 0.15 eV。最后我选 U=3.5 eV 是因为它和实验带隙对得最齐——催化还原计算里 U 不是拟合能垒的工具，而是先对标体相电子结构。泛函上我默认 PBE，做吸附偏弱就补 D3 色散，绝不在没依据时乱换杂化泛函。

四、实操流程：从建模到台阶图

我的标准动作：确定晶面与缺陷位点 → 搭建 slab 并充分弛豫 → 气相分子单独单点 → 逐个吸附构型优化（每个中间体至少试 2 个吸附位）→ 过渡态用 CI-NEB 插 5~7 个点搜 → 频率验证（反应物无虚频、TS 单一虚频）→ 汇总 ΔG 画台阶图。最容易被砍的是 TS 验证，但限速步能垒的可信度全靠它，我宁可多花半天也不跳过。

补一个具体的数字：在 Cu(211) 台阶位上算 CO₂ 到 *COOH，气相参考下 ΔG 是 0.85 eV，加上隐式溶剂和 U=0 V 修正后变成 0.62 eV，而实验报道的过电位对应约 0.58 eV。差 0.04 eV 在可接受范围，但如果不加溶剂，0.85 eV 会让你误判这个位点”不好”，直接排除掉一个其实最优的候选。这件事之后，我做任何催化还原计算的液相体系，溶剂和电势都是默认项，不再当作可选项。

五、常见踩坑：溶剂与电势被忽略

催化还原计算在电化学条件下跑，不加隐式溶剂（如 VASPsol）会系统性高估吸附，不加电位修正会让台阶图失去电压意义。我见过直接用气相机理解释液相还原的，能垒差出 0.3 eV 以上。另一个坑是只算最稳吸附位，漏掉动力学上更易到达的次稳位，过渡态搜索要两者都覆盖。

六、复盘：催化还原计算价值在”指路”

回过头看，催化还原计算最有用的是在实验筛选前告诉你”哪个位点、哪一步最难过”，把试错变成有方向的验证。我接这类项目时默认先要反应体系和目标产物，再决定算哪些中间体。把表面模型、U 值、溶剂电势三件事做对，台阶图才敢拿去和实验对标。被证明有用的，是愿意为一个 U 值反复扫参的扎实。下次有人拿一篇文献的能垒直接套到自己的体系，我会先问——他的溶剂和电位条件，和你的一样吗。