过渡态计算

项目过程：H₂ 在 Cu(111) 表面的解离能垒

这个项目源于客户做铜基催化加氢，需要同时了解两个关键动力学参数：气相 H₂ 在 Cu(111) 上解离成吸附态 H* 的活化能，以及 H* 在 Cu(111) 表面扩散的能垒受应变如何影响。前者决定反应能否启动，后者决定活性位点能否被及时供给。

我分两步做计算。第一步用 CI-NEB 搜索 H₂ 解离过渡态：模型是 4×4 的 Cu(111) 超胞，4 层原子，底部两层固定；H₂ 初始放在表面上方约 2.5 Å，解离态是两个 H 原子分别吸附在 fcc hollow 位；NEB 取 8 个 image，用 climbing image 定位鞍点。第二步用应变扫描研究 H 原子迁移：对 Cu(111) 衬底施加 −3% 到 +3% 的双轴应变，重新计算 H 扩散的最小能量路径。

案例分析：解离能垒 0.45 eV，H 扩散能垒对应变敏感

第一张案例图是 H₂ 解离的 CI-NEB 最小能量路径。初态 H₂+Cu(111) 能量设为零，随着反应坐标推进，能量在第 5 个 image 处达到鞍点，对应解离活化能 E_a=0.45 eV；末态 2H* 的能量比初态低 0.30 eV，说明解离过程放热。0.45 eV 的能垒意味着室温下 H₂ 在 Cu(111) 上的解离速率较慢，需要一定温度或应变调控来活化。

第二张案例图是 H 原子迁移能垒随双轴应变 ε 的变化。未应变时迁移能垒约 2.35 eV；压缩应变（ε<0）使能垒升高，拉伸应变（ε>0）使能垒降低，整体呈线性关系，灵敏度约为 0.18 eV/%（图中标注为 ∂E_a/∂ε 的绝对值）。这说明晶格拉伸不仅能促进 H₂ 解离，也会加速 H* 在表面的扩散，有利于催化循环。

总结教训：不要把不同过程的能垒混为一谈

第一，解离能垒和扩散能垒物理意义不同。0.45 eV 是 H–H 键断裂的能垒，2.35 eV 是 H 原子在表面 hopping 的能垒，两者差了近 2 eV，写文章时必须明确区分。

第二，NEB 的 image 数量要匹配过程复杂度。H₂ 解离用 8 个 image 足够；H 扩散路径更短，5 个 image 就够了，过多反而浪费计算量。

第三，应变计算中衬底层数很关键。4 层是底线，6 层以上迁移能垒差异小于 0.03 eV，但计算量增加约 60%。

项目经验：多步骤能垒要综合看催化活性

只看 H₂ 解离能垒 0.45 eV，会觉得 Cu(111) 活性还行；但加上 H 扩散能垒 2.35 eV，就会发现表面 H* 的迁移可能是限速步骤。做催化机理分析时，最好把吸附、解离、扩散、脱附四步的能垒都列出来，找出真正的瓶颈步骤。这个项目让我认识到：给客户讲过渡态时，要把“哪一步最难”说清楚，而不是只给一个数字。

这篇文章把我在这个项目里用到的计算流程、参数取舍和图上读数做了完整复盘。如果你也在做类似方向，可以参考我的思路：过渡态计算。