迁移能垒理论计算：NEB/CI-NEB的初猜策略与鞍点定位

迁移能垒理论计算在锂离子电池固态电解质和催化反应机制研究中是绕不开的话题。CI-NEB方法看似成熟，但实际项目中image数量设少了鞍点跑偏、设多了计算量翻倍、弹簧常数调不对整条链散架——这些坑本质上都指向同一个问题：迁移路径的连续性如何用离散的image序列可靠地逼近。

NEB vs CI-NEB：一个弹簧常数之差，鞍点精度翻倍

基础NEB方法在每条反应路径上撒N个image，相邻image之间用虚拟弹簧连接，每个image在垂直于路径方向上不受力、在平行方向受弹簧力约束。这种并行优化的结果是整条路径的能量profile可以一次性算出来，但最高能量image不保证落在势能面的真实鞍点上——它只是N个离散image中碰巧最高的那个。

CI-NEB（Climbing Image NEB）在这个框架上加了一个关键改动：最高能量image不受弹簧力约束，而是沿势能面向上”爬”向鞍点。这意味着CI-NEB定位的鞍点精度只受势能面收敛标准控制，不再受image间距的污染。

一个在Li₃PS₄固态电解质中算Li+空位迁移的项目里，这一差异的体感极其清晰。用11个image的NEB，迁移能垒1.32 eV；切到CI-NEB同样11个image，能垒收敛到1.47 eV——0.15 eV的差距不来自路径本身的偏差，而是NEB的最高image偏离了真实鞍点位置约0.3埃。在Arrhenius公式里，0.15 eV的能垒差异在室温下对应约300倍的离子电导率差异，这个错误量级不是可以忽略的边缘效应。

image数量：太少漏鞍点，太多增噪声

image数的选择有一个简单的原则：相邻image之间的原子位移量控制在0.1-0.2埃范围内。以Li+从一个八面体位迁移到相邻四面体位为例，跳跃距离约2.5埃——7-11个image基本覆盖这个位移量。

遇到过两个极端情况：image太少（3-5个），最高能量image离真实鞍点可能差半个埃，整条路径的profile看起来是”梯子”而非平滑曲线——这是离散化错误，不是物理。image太多（15-20个），收敛变慢且弹簧常数需要重新标定，高估的弹簧力会拖后腿。

弹簧常数的默认值（5.0 eV/Ų）在大多数过渡金属氧化物体系中工作良好，但在含Li+的轻离子体系中可能会偏高——轻离子的共价作用弱，势能面更平坦，弹簧力稍强就会在过度约束下画出”比实际更直”的路径，低估能垒。实践中的策略：先用默认值跑一轮，看最高image的位置是否恰好处于键断裂/形成的”中点”附近——偏左或偏右超过10%路径长度，就要检查弹簧常数和image间距。

初猜路径：线性插值之外的选择

CI-NEB最常用的初猜是线性插值（初态和末态之间做直线内插），对于简单的空位跳跃迁移，这个初猜通常足够接近真实MEP。

但遇到两种情况时线性插值会出问题。第一种是协同迁移——多个离子同时移动时，线性插值可能产生非物理的原子重叠，导致初始image的能量异常高甚至SCF不收敛。解决方案是用IDPP（Image-Dependent Pair Potential）方法——它基于原子间距离而非笛卡尔坐标做插值，自然避开原子碰撞。

第二种是结构重排型迁移，如LiFePO₄中Li+沿[010]方向的曲折扩散——线性插值产生的路径穿过了不稳定的中间构型，NEB链在优化过程中可能”绕道”到另一条拓扑路径。此时需要分段处理：将初态-中间态-末态拆成两段独立的NEB，分别搜索MEP。

收敛诊断：力的阈值和路径一致性

CI-NEB的收敛判据通常设每个image上原子受力<0.03 eV/Å。但只靠受力阈值不够——还要检查能量profile的平滑性。如果最高image两边的斜率不连续，或者profile出现”平顶”（即多个image都在相近的能量），说明弹簧网络还没有达到真正的MEP。此时需要增加图像数量的同时适当减小弹簧常数。

一个反直觉的经验：如果跑完CI-NEB后最高image的力仍然在0.05 eV/Å以上但能量profile已显现明确的单峰形态，大概率是image间距不够——在鞍点区域加密image比继续降低力的收敛阈值更有效。追加2-4个image在原最高image附近重新跑，比无限压低力阈值更快收敛到正确的鞍点能量。

更多内容请访问 https://www.keyanxueshu.com/