VASP 反应路径与过渡态入门
以H与H₂交换反应说明端点、插值、爬坡NEB和虚频检查。
阅读本篇目录
反应物和产物的能量差不能告诉我们反应有多难发生。寻找势垒需要连接两个端点，再找路径上的鞍点。本例实际计算H+H₂→H₂+H，用最小模型展示步骤。
1. 先优化端点#
在14 Å立方盒中放入三个H，保持相同原子顺序、晶胞和自旋设置。两端分别对应H1远离H2–H3、H3远离H1–H2。端点先独立优化到0.02 eV/Å。
本例使用PBE、550 eV、Gamma点、ISPIN=2、NUPDOWN=1。必须核对两个端点是想研究的态，不能仅凭文件名叫initial/final就开始算路径。
2. 再优化整条路径#
- 固定优化后的两个端点，在中间插入3张图像，总共5张。
- 保持元素与原子编号一一对应，使用IDPP改善初始插值；检查有没有原子碰撞或错误跨周期边界。
- 用ASE的爬坡NEB连接图像，弹簧常数0.1 eV/Ų，FIRE优化，目标NEB投影力0.02 eV/Å。
- 每次图像受力由VASP静态计算提供；使用Gamma专用程序计算本例单Gamma模型。
- 路径结束后画每张图像相对端点的能量，并检查最高点附近几何。
- 对过渡态做独立有限位移振动检查，确认不稳定模式对应预期成键/断键方向。
这条路线由ASE控制VASP图像计算，不要求把本例输入改成VTST多图像目录。资料包含最终五张图像和绘图数据；重复路径优化时应使用配套ASE控制脚本，不要把五个独立静态任务当成NEB优化。
3. 本次结果#
| 量 | 本次值 |
|---|---|
| 正向势垒 | 0.162532 eV |
| 逆向势垒 | 0.162533 eV |
| 最大NEB投影力 | 0.018531 eV/Å |
| 过渡态显著虚频 | 230.07i cm⁻¹ |
| 过渡态两段H–H距离 | 0.93630 / 0.93635 Å |
Gamma专用程序与标准程序的端点能量差均小于0.0001 eV，远小于本例势垒。两个端点等价，正反势垒也基本相同。振动分析中只有一个明显负曲率模式；其余接近零的模式需要结合整体平移/转动理解，不能把所有小负值都视为第二个反应方向。
4. 势垒和速率有什么区别#
本页是电子势能势垒。若要温度相关速率，还需要反应物与过渡态的振动、平动/转动或表面约束处理，构造相应活化自由能与标准态。轻氢还可能有隧穿效应。因此不能只把0.163 eV代进一个未说明的前因子就当成真实实验速率。
复杂表面扩散或催化反应可能有多条路径；增加图像、改变初始插值和检查相邻极小值都很重要。一个成功路径不能证明它是所有可能路径中最低的。
计算环境与资料下载#
本页数据来自鸢算科技集群实际计算。VASP 6.6.1 GPU版采用 NVHPC 25.3、CUDA 12.8、HPC-X 2.22.1 与 MKL 2023.2；PBE、PAW PBE 版本64、ENCUT=550 eV、EDIFF=1E-6。需要优化的自由原子采用0.02 eV/Å；静态位移点、外加应变点和动力学轨迹不按结构优化阈值判断。
下载输入、分析脚本和实算数据（ZIP） · 下载结果表（CSV）
解压后先读 README，选择本文对应目录，按POSCAR元素顺序准备本人有权使用的POTCAR，再运行 sbatch run.sh。POTCAR、WAVECAR及电荷密度大文件不随资料包分发。各目录的原始输入和实际结果表用于复现；更换体系后应重新提取自己的结果。