反应路径势能模拟计算

这张图是本项目用 CI-NEB 算的一条反应路径：横轴是反应坐标（NEB image 编号），纵轴是势能。起点 IS（initial state）能量归到 0，经过一个中间体 IM，爬到最高点 TS（transition state），再落到终态 FS，整体放热 ΔE = −0.38 eV，活化能 Ea = 1.85 eV。做催化或反应机理研究时，这张图基本就是结论的一半：有没有能垒、能垒多高、反应是吸热还是放热，全在上面。

这张图在说什么

曲线从 IS 出发，先小涨到 IM（一个相对稳定的中间体），再爬到 TS。TS 处的红箭头标了 Ea = 1.85 eV，它是从 IS 到 TS 的能量差。整个反应终点 FS 比 IS 低 0.38 eV，说明热力学上有利。注意 IM 和 FS 之间还有一段起伏，这不是过渡态，而是路径上另一个局部极小或由于 NEB 图像分布造成的台阶——最终只认最高点的真正 TS。

CI-NEB 怎么跑

NEB（Nudged Elastic Band）是在反应物和产物之间插一堆中间构型，用弹性力把它们连起来找最低能量路径。普通 NEB 可能找不到准确的 TS，CI-NEB（Climbing Image NEB）会让能量最高的那个 image 爬坡，直到它刚好落在鞍点。

本项目常用设置：

| 参数 | 取值 | 说明 |

|——|——|——|

| IMAGES | 7–11 | 图像数太少路径不光滑，太多算不动 |

| IBRION | 1 / 2 | CI-NEB 用拟牛顿或共轭梯度 |

| SPRING | −5 | 弹簧力，负值表示用改进 NEB |

| IOPT | 1–3 | 优化算法，GN/DN/L-BFGS |

| ENCUT | ≥450 eV | 与结构优化一致 |

| PREC | Accurate | 能量差才可靠 |

插点方法

插点质量决定 NEB 能不能收敛到合理路径。本项目两种做法：

1. **线性插点**：直接对 IS 和 FS 的坐标做线性插值，简单但路径可能经过高能区。
2. **IDPP 插点**：考虑原子间距离做智能插值，生成更光滑的初始路径，收敛更快。

对复杂表面反应，本项目优先用 IDPP，再用线性插点做对比，确保 TS 不依赖初始路径。

TS 验证

CI-NEB 找到的最高点不一定是真正 TS，必须用振动频率验证：

– 只有 **一个虚频**（负频率），且振动方向沿反应坐标，才是真 TS。

– 虚频数 >1 说明你卡在更高阶鞍点或初始猜测有问题。

本项目习惯在 TS 结构上算 Hessian，并做 IRC（intrinsic reaction coordinate）或能量最小路径反向验证，确认虚频方向连接的是 IS 和 FS。

能垒怎么用

得到 Ea = 1.85 eV 后，可以用 Arrhenius 关系估算速率常数：

k = (k_B·T/h)·exp(−Ea/k_B·T)

室温下 1.85 eV 对应的速率非常慢，说明本反应在常温下不易发生，需要升温或换催化剂降垒。如果做催化剂筛选，目标就是让 Ea 降到 0.6–1.0 eV 区间。

常见坑

– **图像数不足**：5 个 image 常把 TS 抹平，能垒算低。

– **自由度冻结错误**：表面反应中底层原子是否固定会显著改变能垒，本项目通常弛豫表层 2–3 层。

– **自旋态不一致**：若反应涉及自旋翻转，IS/TS/FS 必须同自旋多重度，否则能量差无意义。

[反应路径势能模拟计算](https://www.keyanxueshu.com/category/dft/) 栏目有更多反应路径与过渡态搜索的实战经验；完整流程模板见科研学术网首页 https://www.keyanxueshu.com/。

局限

NEB/CI-NEB 给出的是 0 K 静态能垒；真实反应在有限温度下有零点能、熵、隧穿效应修正，催化体系还需考虑覆盖度、溶剂化和电场。