反应路径势能模拟计算

反应路径势能模拟计算用以定位反应物、过渡态、中间体与产物，并给出能垒与反应热，是理解催化机理、酶催化与材料相变的计算骨架。本项目承接路径计算时，不迷信单一方法，而是用”初猜（NEB/动力学）+ 精修（TS 优化）+ 验证（频率/IRC）”三级流程确保过渡态真实。下文围绕路径计算的完整链路展开。

路径计算的三类目标

1. 最小能量路径（MEP）：反应沿最低能垒走的那条线，对应实际最可能机理。
2. 过渡态（TS）：一阶鞍点，唯一虚频指向反应坐标，其能量即能垒。
3. 反应热 / 自由能变：产物与反应物能量差，决定热力学倾向。

初猜：NEB 与简谐过渡态

• NEB（微动弹性带）：用一串 Images 插值，弹簧力把路径拉直，本项目对固体表面反应、团簇重排默认用 CI-NEB（爬坡修正）找鞍点。
• Dimer / 自适应对少自由度体系更快。
• 初猜路径质量决定收敛速度，本项目先用粗 NEB（少量 images）探路，再加密精修。

过渡态精修与验证

找到候选 TS 后必须验证，否则可能是假鞍点：

4. 频率分析：恰好一个虚频（负本征值），其余为正。本项目用虚频振动模式确认它确实沿反应坐标。
5. IRC（内禀反应坐标）：从 TS 向两侧跑，确认一端连反应物、一端连产物。连错说明找到的是别的路径。
6. 能垒交叉验证：TS 能量应高于两端；若低于反应物，是极小值不是鞍点。

方法层级

本项目对能垒敏感结论用两种泛函交叉验证（如 PBE vs r2SCAN），避免单泛函系统误差。

案例：表面 C–H 活化

某金属表面甲烷 C–H 活化，CI-NEB 给出 TS 能垒 0.82 eV，频率确认唯一虚频对应 C–H 拉伸+金属接近，IRC 顺利连到吸附 CH₃ + H 产物。PBE 与 r2SCAN 能垒差 0.07 eV，结论稳健。该值解释了实验观测的低温活性阈值。

坦诚的局限

NEB 给出的是最小能量路径，但真实反应可能走更高能但熵更有利的路径；多通道反应需逐个排查。本项目对关键机理会并列 2–3 条候选路径比较能垒与自由能，而不是只报一条。反应路径与第一性原理计算的衔接可查本站的第一性原理计算栏目。更多科研计算服务见科研学术网首页。