2. 第一性原理计算基础回顾
在NEB（Nudged Elastic Band，微动弹性带）计算中，我们真正关心的量并不是某个构型的绝对总能量，而是初态与鞍点之间的能量差，也就是扩散势垒。这个量通常在0.3～2.0 eV量级，而DFT总能量本身却可以达到-1000 eV量级。也就是说，势垒是"两个大数相减"的结果，任何系统性误差都会在相减时被放大。因此，理解DFT总能量的构成、交换关联泛函、赝势、平面波截断能和k点采样这五个要素如何影响能量差，是判断NEB结果是否可信的前提。
一、DFT总能量：势垒是能量差，不是绝对值
Kohn-Sham框架下，体系总能量可写成：
E_total = T_s[ρ] + E_ne[ρ] + E_H[ρ] + E_xc[ρ] + E_ion-ion
= 动能 + 电子-离子吸引 + Hartree排斥 + 交换关联 + 离子间排斥其中每一项都依赖电子密度ρ，而ρ又由占据的Kohn-Sham轨道求和得到。NEB计算中，每个image都要独立做一次自洽场（SCF）收敛，得到该构型的E_total，再沿反应坐标构造能量曲线。因此，SCF收敛精度直接决定势垒的可重复性。工程上建议：
- 总能收敛阈值
EDIFF = 1E-6 ~ 1E-5 eV（VASP中对应EDIFF）； - 力收敛阈值
EDIFFG = -0.02 ~ -0.05 eV/Å，NEB的力收敛比能量收敛更关键； - 若EDIFF从1E-5放宽到1E-4，势垒可能漂移0.02～0.05 eV，对轻元素（H、Li）扩散尤其明显。
需要强调的是，NEB的中间image并不需要完全收敛到基态，但必须保证所有image使用相同的计算参数，否则能量曲线会出现非物理的"锯齿"。
二、交换关联泛函：决定势垒的系统偏差
交换关联泛函E_xc是DFT中唯一需要近似的项，也是势垒系统误差的主要来源。常见选择及其对扩散势垒的影响：
| 泛函 | 类型 | 典型影响 | 适用场景 |
|---|---|---|---|
| LDA (PZ) | 局域密度近似 | 势垒普遍偏低，约低估10%～20% | 金属体系快速估算 |
| PBE / PW91 | 广义梯度近似(GGA) | 金属、氧化物势垒较均衡，常用基准 | 大多数扩散计算 |
| RPBE | GGA修正 | 吸附能更准，势垒略高于PBE | 表面扩散、催化 |
| SCAN / r2SCAN | meta-GGA | 势垒精度提升，但成本高 | 高精度对比 |
| HSE06 | 杂化泛函 | 带隙修正，势垒变化可达0.1～0.3 eV | 半导体、绝缘体 |
实践建议：同一组NEB计算必须使用同一泛函。若文献报道某体系势垒为0.65 eV（PBE），而你用LDA得到0.52 eV，这并不一定是错误，而是泛函差异。做对比时，应明确标注泛函版本。
三、赝势：冻结芯电子带来的可移植性
赝势（Pseudopotential）将芯电子冻结，只显式处理价电子，从而大幅降低平面波数量。赝势对势垒的影响主要体现在：
- 价电子配置：例如Ti用3s3p4s3d还是4s3d，会影响过渡态电子结构；
- 截断半径：半径越大，赝势越"软"，所需截断能越低，但可移植性可能下降；
- PAW vs USPP：VASP中PAW（Projector Augmented Wave）通常比超软赝势精度更高，推荐用于NEB。
一个常见陷阱是：初态和鞍点使用了不同版本的赝势（例如从POTCAR库中混用），导致能量差出现0.05 eV以上的虚假偏移。因此，整个NEB路径必须使用完全相同的POTCAR。
四、平面波基组与截断能：收敛性检查不可跳过
平面波基组由截断能 ENCUT 控制。ENCUT不足会导致总能量偏高、力计算不准，进而使势垒出现非单调漂移。收敛测试的标准做法是：
# 示例：对同一扩散体系做ENCUT收敛测试
ENCUT = 300, 350, 400, 450, 500 eV
记录：E_initial, E_saddle, E_barrier = E_saddle - E_initial
判据：相邻ENCUT下势垒变化 < 0.02 eV 即认为收敛经验上，含3d过渡金属的氧化物体系，ENCUT通常需要400～520 eV；含轻元素（H、Li）的体系可适当降低，但不应低于350 eV。注意：ENCUT收敛测试应针对势垒，而不是仅针对总能量，因为势垒对基组误差更敏感。
五、k点采样：倒空间积分的精度控制
k点采样决定布里渊区积分的精度。对于扩散势垒，k点不足会导致：
- 总能收敛但力场出现噪声，NEB路径抖动；
- 鞍点位置偏移，势垒误差可达0.05～0.15 eV；
- 金属体系对k点更敏感，半导体/绝缘体相对宽松。
推荐使用Monkhorst-Pack网格，并做k点收敛测试。例如：
# 示例：k点收敛测试（以4x4x4超胞为例）
k-mesh = 2x2x2, 3x3x3, 4x4x4, 5x5x5
观察：E_barrier随k点加密的变化
判据：势垒变化 < 0.03 eV，且力收敛稳定对于大超胞（如5x5x5以上），1x1x1或2x2x2的Gamma-centered网格有时已足够，但必须通过测试确认。切忌直接套用文献k点而不做验证。
- 泛函：全路径统一，记录版本（PBE/SCAN/HSE06）；
- 赝势：全路径统一POTCAR，优先PAW；
- ENCUT：针对势垒做收敛测试，变化<0.02 eV；
- k点：针对势垒做收敛测试，变化<0.03 eV；
- EDIFF：1E-6～1E-5 eV；EDIFFG：-0.02～-0.05 eV/Å；
- 所有image参数完全一致，仅原子位置不同。
把上述五个要素串起来看，NEB势垒的精度控制本质上是一个误差预算问题：泛函决定系统偏差（通常0.05～0.2 eV），赝势和ENCUT决定基组误差（0.02～0.05 eV），k点决定倒空间积分误差（0.03～0.1 eV），SCF阈值决定数值噪声（0.01～0.03 eV）。这些误差不会自动抵消，而是可能叠加。因此，在报告一个扩散势垒时，至少应说明泛函、赝势类型、ENCUT、k点网格和力收敛标准，否则结果无法被复现或横向比较。
最后需要提醒：NEB计算中，初态和末态的结构优化应使用与NEB路径完全相同的参数。如果初态用ENCUT=500优化，而NEB用ENCUT=400，那么初态能量与路径起点能量不一致，势垒会被系统性高估或低估。参数一致性是NEB计算中最容易被忽视、却最影响结果可靠性的工程细节。