VASP计算结合能：从参数设置到结果验证的完整流程

VASP计算结合能在材料模拟领域几乎是评估结构稳定性的”标配”操作，但真正做起来远没有手册上写的那么轻松。这个项目曾在一组Ti-Al二元合金的结合能计算上卡了两周——结果反复波动，根源就藏在几个被忽略的参数细节里。

结合能的本质是把孤立原子拉拢成固体时释放的能量，数值越负意味着结构越稳定。计算公式并不复杂：

E_b = E_total – Σ E_atom

其中E_total是体系的总能，E_atom是各孤立原子的能量。看似直接，但问题在于：孤立原子能量的计算方式直接决定了最终结果的可靠性。项目组最初采用同一组ENCUT参数计算两者，却发现原子能量偏高约0.3 eV，导致结合能系统性偏正。调整后改用更高截断能单独计算原子能量，结果才收敛到文献参考值附近。

VASP提供了两类赝势，项目在这一步做了明确取舍。PAW（投影缀加波）方法在过渡金属体系中的精度更被认可，尤其是对d电子的描述更为准确。根据Physical Review B的系统性比较研究，PAW赝势在金属间化合物中的结合能误差通常在10 meV/atom以内，而Ultrasoft在某些含d轨道元素中偏差可达50 meV以上。

对于Ti-Al体系，项目选定了PAW_PBE版本的Ti_pv和Al赝势。_pv后缀意味着对Ti的半芯p态进行了显式处理，这对于涉及价电子跃迁的结合能计算尤为关键。选定之后不再反复更换——赝势的一致性是结果可比性的前提。

ENCUT截断能的收敛测试是绕不过去的一关。项目从400 eV起步，每50 eV递增，观察总能变化：

– 400 eV时，结合能为-4.12 eV/atom

– 450 eV时，结合能为-4.28 eV/atom

– 500 eV时，结合能为-4.31 eV/atom

– 550 eV时，结合能为-4.32 eV/atom

从500 eV到550 eV，变化已小于10 meV，项目最终采用520 eV作为计算截断能。差距不会说谎——低于450 eV时的结果偏差接近4%，足以改变稳定性排序的结论。

K点网格同样需要收敛。对于Ti3Al的六方原胞，项目从4×4×2逐步加密到12×12×4。在8×8×3以上，结合能变化进入5 meV以内的平台区。最终选用Gamma中心化的9×9×3网格，兼顾精度与计算效率。

结构弛豫并非”一键运行”那么简单。项目对比了两种方案：ISIF=3同时优化晶胞形状、体积和原子位置，与ISIF=2仅优化原子位置的固定晶胞计算。对于实验晶格常数已知的体系，ISIF=2的结果与实验值吻合更好，结合能偏差约2%。而ISIF=3在晶胞优化过度时，体积可能偏大1-2%，间接影响结合能。

EDIFFG的设置也需要注意。项目采用-0.01 eV/Å的力收敛判据，而非默认的-0.02。更严格的判据对含轻元素的体系尤为必要，Al原子的位置敏感性较高，松弛的判据可能让原子未能到达真正的力极小值点。

这个项目踩过几个值得警醒的坑。其一是自旋极化（ISPIN）的设置：Ti元素虽然基态非磁性，但在某些低配位环境中会显现局域磁矩。项目初期未开启ISPIN=1，结果在Ti-rich组成的计算中总能高出0.05 eV/atom，结合能偏正。开启后，磁矩收敛到接近零，总能稳定。对于含过渡金属的体系，默认开启自旋极化计算是更稳妥的策略。

其二是SIGMA值的选择。金属体系使用Methfessel-Paxton方案时，SIGMA=0.2是常见值，但对于带隙较小的半导体相，过大的SIGMA会导致占据数”抹平”，影响总能精度。项目在计算TiAl3（带隙接近零的金属间化合物）时，将SIGMA降至0.05并配合ISMEAR=0的高斯展宽，总能变化进入0.5 meV量级。

验证方面，项目参考了两条基准线：一是实验报道的Ti3Al形成焓约-72 kJ/mol（约-0.75 eV/atom），计算值为-0.73 eV/atom，吻合度在3%以内；二是与Materials Project数据库中同一结构的条目对比，偏差小于15 meV/atom。双交叉验证通过后，这组结合能数据才算真正站住脚。

回顾整个计算过程，参数收敛测试占据了将近60%的工时，但正是这些看似枯燥的逐点扫描保证了最终结果的可信度。VASP计算结合能从来不是跑通INCAR那么简单——每一个参数背后都对应着明确的物理判断，而结果的可重复性恰恰建立在这些判断之上。