反应路径势能模拟计算

一、项目背景与科学问题

氢气分子在过渡金属表面的解离吸附是多相催化、储氢材料和燃料电池电极反应中的核心基元步骤。以 Cu(111) 表面为例，H₂ 从气相自由分子到两个解离吸附氢原子（2H*）的过程，需要跨越一个能量势垒，该势垒的高度直接决定了解离速率常数，进而影响催化活性与选择性。若能在原子尺度准确获得这一反应路径上的能量分布，就能预测表面在给定温度、压力条件下的催化行为，从而为催化剂设计提供定量依据。

本项目围绕”氢分子在 Cu(111) 表面解离吸附的反应路径势能面”这一科学问题，采用第一性原理密度泛函理论（DFT）结合爬坡弹性带（CI-NEB）方法，计算最小能量路径（MEP）上各个中间构型的能量，识别过渡态（TS）位置，并提取活化能与反应热两个关键热力学与动力学参数。

二、计算方法与理论背景

2.1 密度泛函理论

采用 VASP 平面波基组 DFT 程序，交换关联泛函选用 GGA-PBE 形式。Cu 的投影缀加波（PAW）赝势价电子构型取 3d¹⁰4s¹，平面波截断能设置为 400 eV，以保证总能量收敛到 1 meV/atom 以内。Cu(111) 表面采用周期性 slab 模型，真空层厚度大于 15 Å，用以消除相邻镜像层之间的相互作用。

2.2 NEB 与 CI-NEB 方法

弹性带（NEB）方法将反应路径离散为一组镜像结构（本计算取 7 个镜像），每个镜像与其相邻镜像之间通过弹簧力耦合，通过优化使所有镜像同时达到能量最低构型，从而勾勒出反应路径的能量剖面。爬坡弹性带（CI-NEB）在 NEB 基础上对能量最高的镜像施加力反转，使其收敛到精确的一阶鞍点，即过渡态。过渡态需通过振动分析验证，确保其仅有一个虚频。

三、计算设置与模型构建

• 表面模型：Cu(111) 4×4 超胞、4 层原子厚度，底层 2 层固定以模拟体相约束，顶层 2 层与吸附质自由弛豫。
• k 点采样：采用 Monkhorst–Pack 网格，表面 Brillouin 区内 k 点精度设置为标准水平（等效于 3×3×1 量级），满足收敛测试要求。
• 收敛判据：原子受力 < 0.02 eV/Å，能量收敛 < 1×10⁻⁵ eV。
• NEB 镜像数：7 个，初始路径采用线性插值生成，随后进行 CI-NEB 优化。

四、结果分析与案例图解读

图 1 的左面板给出 NEB 计算得到的一维能量剖面。横坐标为反应坐标（s），纵坐标为体系能量。从图中可以看到，路径起点为气相 H₂ 分子（能量定为 0 eV），终点为两个解离吸附的氢原子 2H*（能量为 −0.40 eV）。路径中能量最高点位于 s = 0.5 处，活化能 Eₐ = 0.50 eV，该位置即过渡态。计算得到的解离吸附反应能为 ΔE = −0.40 eV，属于放热反应，说明 Cu(111) 表面对氢的解离吸附在热力学上有利，但动力学上仍需克服约 0.5 eV 的势垒。这一数值与文献报道的 DFT 结果（0.45–0.55 eV）一致，验证了计算方案的可靠性。

图 1 的右面板绘制了二维势能面的等高线图，横纵坐标分别为 H–H 键长和 H–Cu 距离。图中星号标记的过渡态位于势垒鞍点位置，两个能量极小值分别对应气相 H₂（左侧）和吸附态 2H*（右侧）。二维势能面直观展示了反应通道的走向：体系从反应物区域沿最小能量通道越过鞍点到达产物区域。等高线密集的方向表明该方向上势能变化剧烈，对应 H–H 键断裂的化学过程。

五、结论与展望

本项目通过 DFT + CI-NEB 方法定量获得了 H₂ 在 Cu(111) 表面解离吸附的反应路径势能面，活化能 Eₐ = 0.50 eV、反应能 ΔE = −0.40 eV。该结果与文献值吻合，表明方法可靠。案例图以一维能量剖面配合二维势能面两种形式呈现，既给出关键动力学参数，又展示反应通道的全局特征，可支撑论文中”反应机理与能量学”章节的论述。

后续可在以下方向拓展：(1) 引入零点能（ZPE）修正，获得更精确的活化能；(2) 考虑覆盖度效应，研究高覆盖度下活化能的变化；(3) 将 NEB 获得的过渡态用于速率常数计算（如 TST/hTST），定量预测催化反应速率。