过渡态计算

结合案例图分析，第一张案例图横轴是反应坐标（从0到1归一化），纵轴是能量。灰色圆点是反应物，绿色方块是产物，橙色圆点是NEB生成的7个image，蓝色曲线是插值后的能量面。反应物能量设在0 eV附近，产物能量比反应物低约0.25 eV，说明产物在热力学上更稳定。能量最高点出现在反应坐标约0.5的位置，对应过渡态，能量为0.85 eV，这就是正向反应的活化能E_a。从过渡态到产物一侧是下坡，到反应物一侧也是下坡，符合鞍点特征。

第二张案例图用三个椭圆示意了反应过程中的几何变化：反应物是ABC三个片段按一定顺序排列，过渡态中B和C发生相对滑移，产物中C的位置发生翻转。图下标注”NEB 7 images · CI-NEB收敛”，说明用的是攀爬图像弹性带方法（CI-NEB），让最高能量image自动爬到鞍点，并且用7个image已经收敛。

从能量路径图和几何示意图的组合来看，这是一个过渡态计算项目——反应是一个小分子在催化剂表面从左侧吸附态翻转到右侧吸附态的过程，客户关心的是这个重排能不能在室温下发生。能量路径给出了热力学和动力学判据，几何示意图给出了反应过程的直观理解。

项目计算设置上，我用VASP做NEB。初始和末态先用标准DFT优化到力收敛<0.02 eV·Å^{-1}；然后用nebmake.pl生成7个初始image，包括初态和末态共9个点。INCAR里设置IMAGES=7，SPRING=-5，LCLIMB=.TRUE.开启CI-NEB。交换关联用PBE，平面波截断能450 eV，K点按表面超胞取3×3×1。收敛判据是最大力<0.05 eV·Å^{-1}，实际跑了约60步后达到。

根据0.85 eV的活化能，我用过渡态理论估算室温下的速率常数。取尝试频率约10^{13} s^{-1}，k_BT在300 K约为0.026 eV，k≈10^{13}·exp(−0.85/0.026)≈10^{-1} s^{-1}。这意味着在室温下反应可以在秒量级发生，对催化循环来说是可以接受的动力学。如果E_a超过1.2 eV，室温速率会大幅下降到分钟甚至小时量级。

踩坑经验有三条。第一个是初始路径质量：如果nebmake.pl线性插值出来的中间结构出现原子重叠，NEB很容易发散。我习惯先手动调整几个关键原子的位置，或者先用IDPP方法做更平滑的初猜。第二个是spring常数：−5这个值对大多数表面反应合适，但如果能量面很平，image会滑向产物/反应物，需要增大到−10。第三个是收敛判据：对于0.85 eV这么高的活化能，0.05 eV·Å^{-1}的力收敛足够；但如果E_a<0.3 eV，需要更严的0.03 eV·Å^{-1}才能精确定位鞍点。

更多关于过渡态计算（NEB、DIMER、LST/QST）的参数设置与收敛技巧，我在过渡态计算栏目里有系统整理。