二维电极材料里的 Li 迁移势垒，经常被拿来判断倍率性能。真正计算时，难点并不只是把 VTST 的 NEB 参数写进 INCAR，而是先判断一个迁移事件有没有被定义清楚：初态和末态是不是同一种 Li 覆盖度、同一种表面、同一种超胞，迁移前后是否只改变 Li 的占位而没有改变材料骨架身份。NEB 算到的是两个已定义局域极小值之间的最小能量路径，不是任意两张结构之间的能量插值。 后续检查应回到端点和图像串本身。
端点怎样才算可比较？
初态和末态都要是局域稳定构型。对二维材料来说，常见端点包括 hollow 位、bridge 位、top 位、孔洞中心位、层间四面体位或八面体位。端点优化后，如果 Li 从预设位点滑到同一个低能盆地，说明这两个候选位点在当前模型下不是两个独立端点，继续做 NEB 只会得到很短的回落路径。
端点还要保持可比较的化学组成。单个 Li 在 4×4 单层表面扩散，就不要把末态改成另一种 Li 覆盖度；层间扩散也不要让末态多出局域层间滑移。端点一致性检查的核心，是确认能量差来自 Li 位置改变，而不是来自超胞、覆盖度、磁构型或骨架重构改变。 这一步通常比后处理曲线更早决定结果能否比较。
图1 二维硼化物负极材料中不同 Li 储存位点和扩散通道的结构示例。位点命名、覆盖度和迁移方向应先在同一个超胞中固定，再进入 NEB 计算。DOI: 10.3390/nano15120886。
一个实用做法是先做静态位点扫描：把 Li 放在所有高对称位点，分别优化，记录最终位点、相对能和最近邻配位。只有最终仍能分成不同局域极小值的位点，才适合作为 NEB 的候选端点。若多个初始位置都优化到同一位置，它们不是多个扩散终点，而是同一吸附盆地的不同初猜。
Eads = E(material + Li) – E(material) – E(Li reference) 式(1)
式(1)用于先筛 Li 的稳定占位。E(material + Li) 是含 Li 结构总能，E(material) 是同一超胞下材料总能，E(Li reference) 可按研究口径选金属 Li 或其他参照；迁移势垒计算前，应先确认端点位点在同一参照下确实是局域极小值。
扩散路径怎样对应相邻位点？
二维表面常见的错误，是把两个相距很远但能量较低的位点直接连成一条路径。线性插值会把 Li 拉过不合理的短键区域，最高能图像变成几何碰撞，而不是迁移鞍点。更合理的处理是先在结构图上找相邻位点网络，再把长距离扩散拆成若干个相邻跳跃。
Li 从 hollow A 到 hollow B 时，如果中间必须经过 bridge 附近的通道，NEB 端点可以设为 A 和 B，但初始图像要沿 A-bridge-B 分布。层间 Li 从一侧孔洞穿到另一侧孔洞时，如果中间存在亚稳中间位点，就应拆成 A→M 和 M→B 两段。迁移路径的物理身份由相邻位点网络决定，不能只由端点直线距离决定。 路径拆分后，正向和反向势垒也更容易对应具体跳跃。
图2 Li 在二维材料不同扩散方向上的能量剖面示例。不同通道的势垒差异说明迁移路径需要按晶向、位点连接和局域配位分别比较。DOI: 10.3390/nano15120886。
实际输入中，端点 POSCAR 应先用同一精度优化：同一 POTCAR、ENCUT、KPOINTS、ISMEAR/SIGMA、vdW 修正和磁性设置。对二维材料，真空层和偶极修正也要统一；对带电补偿或外加电场模型，端点之间还要保持同一背景设置。只要端点计算口径不同，后面的 NEB 曲线就很难解释。
图像编号代表什么？
NEB 的 IMAGES 不是越多越好。图像太少时，最高点可能跨过真正鞍点；图像过多时，如果初始路径差，图像之间会在低能盆地附近堆积，计算量增加却没有更清楚的鞍点。二维材料中一个相邻 Li 跳跃通常可以从 3–7 个中间图像开始，复杂穿孔或层间路径再增加图像数。
图像要按同一原子映射生成。Li 原子在 POSCAR 中的编号、材料骨架中可移动原子的编号都应一致；如果端点结构经过重新排序，插值会把错误原子连在一起。NEB 图像的横坐标是经过弹簧约束和切向投影修正后的路径坐标，不是普通结构优化步数。 因而图像编号不能直接解释为真实时间序列。
图3 端点对齐示意。端点优化前后的原子映射若发生漂移，插值得到的图像串可能连接到错误的迁移事件。DOI: 10.1016/j.mex.2026.103899。
在 VASP 输入里，常用组合是 IBRION=3、POTIM=0、SPRING=-5，结合 LCLIMB=.TRUE. 做 CI-NEB。初始阶段可以先不用 climbing image，把路径稳定到合理能带；随后打开 LCLIMB 精修最高能图像。力收敛阈值可从 EDIFFG=-0.03 eV/Å 或 -0.05 eV/Å 开始，最终用于发表的势垒需要结合体系尺寸和目标精度重新收敛。
初始路径为什么会改变势垒？
线性插值只适合结构变化简单、没有明显避让动作的迁移。Li 在二维孔道、缺陷附近或层间狭窄通道迁移时，中间图像很容易出现 Li-骨架原子距离过短。此时可以采用 IDPP/SIDPP 思路生成初始路径，或者手动把中间图像放在合理通道上，再进行短步预优化。
同一端点下可以分别用线性插值和 IDPP 路径起算，比较早期图像的最近邻距离、能量峰位置和最终 MEP。如果两个路径最后收敛到同一势垒，说明路径比较稳健；如果势垒差异很大，应回头检查是否存在多个迁移通道。多个初始路径收敛到不同能量剖面时，较低势垒路径才更可能对应动力学可达通道。 这一判断还需要构型连续性作支撑。
图4 SIDPP 路径生长示意。对复杂迁移事件，逐步增加图像并优化中间构型，可以降低插值路径落入错误局域结构的概率。DOI: 10.1016/j.mex.2026.103899。
还有一种常见情况是表面原子在 NEB 过程中发生明显重构。若端点未允许骨架充分弛豫，图像串会把表面松弛和 Li 迁移混在一起；若全部原子都自由移动，又可能出现端点附近骨架差异过大。可操作的折中是先比较固定底层、放松表层和全放松三组端点，确认迁移势垒对约束层数是否敏感。
正向和反向势垒怎样读？
NEB 曲线给出一串相对能量。若初态能量设为 0，最高能图像相对初态的能量就是正向势垒；最高能图像相对末态的能量是反向势垒。端点不等能时，不能只报一个数就结束，因为正向和反向扩散的动力学含义不同。
Em,forward = ETS – Einitial；Em,reverse = ETS – Efinal 式(2)
对于电极材料，低势垒通常说明 Li 扩散快，但这个判断要和终态稳定性一起看。若末态比初态高很多，正向势垒低并不代表 Li 会大量占据末态；若两个端点接近等能，势垒才更接近材料内部扩散的主要限制。迁移势垒描述跨越鞍点的代价，端点能差描述两个位点的热力学偏好。两者合起来才说明扩散方向和占位概率。
图5 一维能量路径示例。初态、鞍点和终态的相对能量共同决定正向势垒、反向势垒和反应能。DOI: 10.1016/j.mex.2026.103899。
路径图应同时标出位点编号、路径方向、相对能量零点和图像数。只给一张平滑曲线而不交代端点结构，会无法判断势垒对应哪一个迁移事件；只给结构图而不标能量，也无法判断该通道是否真正主导扩散。
鞍点身份怎样确认？
NEB 收敛不是只看 OUTCAR 末尾的 force。需要逐个打开 CONTCAR，检查 Li 是否沿预期通道移动，表面骨架是否发生非目标重构，最高能图像附近是否存在明显短键或不合理穿越。如果最高能图像只是几何碰撞，势垒会虚高；如果图像串塌到端点附近，势垒会虚低。
严格一点的过渡态确认，可以取最高能图像或经 CI-NEB 精修后的鞍点结构做频率分析。对单个 Li 迁移事件，过渡态附近应出现沿迁移方向的一条虚频，而不是多个骨架软模共同参与。势垒数值只有和鞍点构型、虚频方向、端点能差一起成立时，才适合用于扩散能力判断。 频率结果可作为鞍点身份的独立证据。
实际执行时，可以先做位点扫描，筛选局域极小值；再按相邻位点网络定义迁移事件；接着生成 3–7 个中间图像并做预检查；NEB 收敛后读取正向和反向势垒；最后用构型和必要的频率分析确认鞍点身份。这样的结果，才能和二维材料中的倍率性能、扩散各向异性和缺陷辅助迁移建立清楚联系。
要点：Li 迁移势垒计算要先定义同一覆盖度、同一超胞和同一材料骨架下的端点。
要点：扩散路径应来自相邻位点网络，不应只按端点直线距离生成插值。
要点：NEB 图像是反应坐标采样点，势垒读取需要同时报告正向势垒、反向势垒和端点能差。
要点：最高能图像还要结合构型连续性、力收敛和必要的虚频方向确认鞍点身份。
