



NEB 计算中最让人头疼的现象之一,是中间图像优化着优化着就靠回初态、终态,或者一起落到某个中间低能结构里。表面上看像是路径“塌回去”,本质上却不是一个单一错误。它可能来自端点没有真正稳定,也可能来自插值结构太差、原子映射不一致,或反应本来就不是一步完成。NEB 图像塌回去,通常说明图像串没有被约束在目标反应坐标附近。
这类问题不能只靠增加 IMAGES 或重新提交解决。NEB 计算要同时满足四件事:初态和终态是同一个基本反应的两端,原子顺序和固定层完全一致,初始路径没有穿越不合理距离,优化后最高能图像能对应真实断键、成键或迁移动作。缺少任意一环,图像串都可能沿势能面滑向更容易到达的低能区域。

图1 端点最小化、原子对齐、SIDPP 路径生成和 CI-NEB-MMF 优化的流程。图中强调端点与初始路径质量会决定后续能否收敛到目标 MEP。DOI: 10.1016/j.mex.2026.103899。




NEB 的初态和终态必须先充分优化。若初态只是一个还没弛豫好的吸附构型,终态只是手动拉开的产物构型,图像串优化时会同时受到端点附近残余力和路径力的影响。中间图像可能不是朝过渡态移动,而是被吸回端点附近的低能盆地。
端点还必须表示同一个反应。比如 H 从 A 位迁到 B 位,初末态中只有 H 位置应该发生目标变化;如果末态同时改变了吸附物取向、表面重构和氢键网络,NEB 会试图一次连接多个变化。端点不清楚时,图像塌回去不是偶然,而是路径定义本身不成立。
VASP 的 NEB 依赖每个图像里完全相同的原子顺序。若初态和终态中两个氢原子、氧原子或表面同种原子被交换,插值会把错误的原子连成反应坐标,图像之间看似连续,实际对应的是原子身份错配。优化时,体系会自动回到更合理的局部结构,表现为图像收缩或乱跳。
有机分子、表面水层、多个 H* 共吸附和大分子断键最容易出现映射漂移。处理时要先固定原子编号,再用可视化工具逐帧检查目标键长、相邻原子距离和吸附位点是否连续。原子顺序一致,不等于化学映射已经正确。

图2 两阶段端点对齐过程示意。端点优化前后都需要处理原子映射漂移,否则插值图像可能连接到错误的反应坐标。DOI: 10.1016/j.mex.2026.103899。




线性插值只按坐标直线移动,不理解键长、排斥和表面几何。若反应涉及 H 绕过表面原子、O—H 键旋转、分子翻转或双位点迁移,简单插值可能让原子穿过表面,或让两个原子距离过近。优化一开始,体系先处理这些不合理近邻,而不是沿目标反应前进。
这时图像能量会出现异常高点,力也可能把图像推回初末态附近。更好的做法是使用 IDPP 或手工构造少数关键中间结构,让相邻图像之间的键长变化平滑。初始路径不是装饰,它决定 NEB 优化从哪个势能通道出发。
NEB 输入准备完成后,应逐帧检查目标键长和关键非键距离。例如水解离要看 O—H 断裂和 H—表面成键是否单调变化;H 迁移要看 H 到起点和终点位的距离是否平滑交替;C—H 断裂要看 C—H 拉长与 H—金属成键是否同时发生。
如果第 2 张图像就已经接近终态,或者中间某张图像出现不合理短距离,后续优化很难救回来。可以重新插值,也可以先把反应拆成“转向”和“断键”两段。图像串的几何连续性,要在正式计算前就确认。

图3 SIDPP 路径生长示意。逐步增加图像并反复优化中间图像,可以降低复杂反应插值时落入局部低能盆地的风险。DOI: 10.1016/j.mex.2026.103899。




CI-NEB 会把最高能图像推向鞍点。如果路径还没有基本成形,最高能图像可能只是一个坏插值结构,而不是靠近真实过渡态的结构。过早开启爬山,会把这个坏图像继续推高,导致路径扭曲,中间图像被拉向端点或被挤到同一低能区域。
更可靠的流程是先普通 NEB 预收敛,让图像串大致贴近最低能路径,再开启 LCLIMB。对于特别复杂的反应,可以先用较松 EDIFFG 得到路径形状,再对最高能附近图像做更精细收敛和频率验证。爬山图像应该用于精修鞍点,而不是替代前期路径整理。
Ea = ETS − EIS 式(1)
式(1)里的 Ea 是正向活化能,ETS 对应路径最高能附近的过渡态候选,EIS 是初态能量。若最高能图像来自坏插值或多步反应的中间扭曲结构,代入式(1)得到的数值就不是目标反应的能垒。
图像数太少时,一条弯曲路径会被粗糙地拉成几段直线,最高能区可能被跳过;图像数太多但初始路径差时,多个图像会挤到低能盆地,路径仍然不清楚。SPRING 设置也会影响图像沿路径的分布,过强或过弱都可能让图像间距失衡。
判断时不要只看 OUTCAR 是否收敛,还要画能量-图像序号曲线,看是否出现连续、单峰或合理多峰。如果多个图像能量几乎相同并堆在端点附近,说明路径没有均匀铺开。NEB 收敛的是一串结构,不只是一个最高能数字。

图4 CI-NEB 与最小模跟随方法衔接示意。先用 CI-NEB 获得全局路径,再对鞍点附近做局部精修,可以减少错误最高能图像带来的偏差。DOI: 10.1016/j.mex.2026.103899。




如果反应实际包含分子转向、质子迁移、断键和产物扩散多个动作,一条 NEB 可能会自动找到中间低能结构。图像落到这个结构附近,并不代表计算失败,而是体系告诉你:这里可能存在一个真实中间态。继续把它硬拉成单步路径,反而会让能垒解释混乱。
此时应把路径拆开:初态到中间态算一段,中间态到终态再算一段。每段只包含一个主要断键、成键或迁移动作。当图像反复落到同一个中间低能构型时,优先判断它是不是应该被升级为中间态。
不一定。如果图像回到端点,通常说明端点、映射或插值有问题;如果图像回到一个稳定中间构型,则可能说明反应需要分段。区别在于这个构型是否可独立优化、是否有清楚的化学身份、是否能作为下一段路径的端点。
对催化表面来说,H* 先迁移再重组、H₂O 先转向再断键、CO₂ 先弯曲再形成 *COOH,都可能需要分段。把多步变化拆开后,每段 NEB 的图像更容易保持连续,虚频验证也更容易对应目标反应坐标。

图5 HCN → HNC 异构化的一维能量路径示例。能量剖面中初态、鞍点和终态的相对位置,可用于理解图像串怎样逐步逼近收敛路径。DOI: 10.1016/j.mex.2026.103899。




第一步,重新优化初态和终态,并确认二者原子数、晶胞、选择性弛豫标记和原子顺序完全一致。第二步,检查目标键长或迁移距离是否沿图像序号平滑变化。第三步,先不开 LCLIMB 做预收敛,再开启 CI-NEB。第四步,对最高能图像做频率或 dimer 精修,确认唯一虚频方向对准目标反应。
FiNEB = −∇E⊥ + Fspring,∥ 式(2)
式(2)表示 NEB 中每张图像主要受到垂直于路径的真实力和沿路径的弹簧力共同作用。若路径初始方向错了,垂直投影和弹簧约束都会围绕错误路径工作,图像自然可能向低能端点收缩。
SYSTEM = neb_path_precheck
IMAGES = 5
SPRING = -5
LCLIMB = .FALSE. # 先普通 NEB 预收敛
IBRION = 3
POTIM = 0
EDIFF = 1E-5
EDIFFG = -0.05
NSW = 120
# 预收敛后再开启 LCLIMB = .TRUE.,并检查最高能图像的虚频方向。
排查时还要保存每一轮 CONTCAR 并做可视化。若图像 02、03、04 全部向初态靠拢,通常优先检查末态和插值;若所有图像向某个中间构型靠拢,优先检查是否应拆分路径;若最高能图像结构怪异,优先关闭 LCLIMB 重新预收敛。
dk = |Rk+1 − Rk| 式(3)
式(3)可以粗略表示相邻图像的结构距离。正式判断时不一定只用笛卡尔距离,更要看目标键长和关键吸附位点变化。若某两个相邻图像差别很大,而其他图像几乎重合,路径分布就不合理。
最后,能垒报告要写清楚路径身份:初态是什么、终态是什么、目标断键或迁移坐标是什么、是否出现中间态、最高能图像是否经过单虚频验证。这样即使某次 NEB 曾经“塌回去”,也能通过排查过程把问题变成可靠路径,而不是反复提交同一套输入。




要点:NEB 图像塌回去通常来自端点、原子映射、插值路径或多步反应定义问题。
要点:初态和终态必须充分优化,并且对应同一个基本反应。
要点:CI-NEB 应在路径基本成形后再用于鞍点精修。
要点:若图像反复落入同一个中间低能构型,应考虑把反应拆成两段。
