结合案例图分析,本文以这个表面催化决速步项目为例,讲 dimer 找鞍点和 ZPE 校正。这是一个表面催化客户做的决速步活化能搜索项目。客户研究电催化反应,需要知道最慢一步的能垒来估计反应速率。第一张案例图是 dimer 方法得到的反应能量曲线,沿反应坐标在约 0.5 处出现一个峰,标注 dimer TS,活化能 Ea=0.90 eV。第二张案例图是 dimer 迭代中最大力的收敛曲线。

方法上 dimer 的优势是不需要预知完整路径,只要一个初始猜测结构和反应坐标方向就能自动找最近鞍点,比 NEB 更适合路径复杂的表面反应。但我先用 NEB 粗略探路确认大致坐标范围,再用 dimer 精修——单独用 dimer 容易卡在局部极小。 我的设置:VASP 优化反应物和产物,线性插值得初始结构,dimer 优化到鞍点(EDIFFG=−0.05 eV/Å),最后做频率验证。客户最初只要”Ea 一个数”,我坚持给了过渡态结构和虚频,他后来做速率估算正好用上。
图里读出的数要落地。第一步确认 dimer 收敛——最大力 < 0.05 eV/Å 才算鞍点找到,否则峰值是假峰。我画过 dimer 最大力随迭代步的收敛曲线:初始 dimer 方向偏离真反应坐标时,前 5 步力在 2→1 eV/Å 徘徊,修正方向后第 8 步起快速下降,15 步到 0.03 eV/Å 平台——这条曲线证明鞍点已收敛。 第二步做频率验证:真过渡态有且仅有一个虚频(本例 −320 cm⁻¹)。我把”单虚频”作为交付硬判据:虚频在 −200~−500 cm⁻¹ 区间属典型过渡态,若虚频接近 0 说明鞍点太平、反应坐标没抓住。第三步把 Ea 含 ZPE 报给客户,并标注误差 ±0.05 eV。
做的过程中我踩过实打实的坑。第一坑是 dimer 对初始方向敏感:方向错会收敛到附近极小点而非真鞍点,我让 dimer 方向沿疑似反应坐标(某键伸长、某原子迁移方向),方向修正后收敛曲线才正常下降。 第二坑是 ZPE,含氢体系活化能要做零点能校正,可能改变 0.1–0.3 eV,第一版没做 ZPE 时 Ea=0.78 eV,校正后 0.90 eV。第三坑是只用 NEB 不精修,粗路径给的 Ea 误差可达 0.1 eV,必须 dimer 收尾。
作为可复用的项目经验,第一,dimer 建议先用 NEB 探路再用 dimer 精修,成功率更高。第二,活化能务必含 ZPE 校正并标注 ±0.05 eV。第三,报告要有结构图+虚频+能垒三件套,缺一不可。客户最终以 hollow 位点活化能作为 ORR 决速步标度。说到底,过渡态项目的交付标准就三条——有鞍点、有且仅有一虚频、能垒含 ZPE;缺任何一条,客户拿去算速率都是空中楼阁。
计算工程师团队:500+ | 成功计算案例:5W+ | 1W+ 助力顶刊发布经验
需要模拟代算服务可以联系我们,可免费评估项目的计算费用和免费出计算方案。
cp2k第一性原理计算
CP2K模拟计算
CP2K模拟计算
castep计算吸附能:Materials Studio里的表面建模与能量收敛复盘
CP2K计算能带:大体系电子结构模拟的混合基组方案
能带理论计算:固体能带结构的DFT模拟方法与工程应用
CP2K计算能带:混合基组DFT方法在周期性体系中的实战应用
CP2K分子动力学模拟详解:大体系加速策略与GPW方法实战