手机版

VASP计算结合能：结合能、形成能与内聚能的概念辨析与计算

发布时间：2026-06-24 来源：科研学术网

小中大

字号：

VASP计算结合能的项目里，项目组见过最频繁的错误是把形成能的计算结果填到论文表格里标”结合能”。这两者差的不只是名字——参考态不同、物理含义不同、数值可以差几个eV。在FeAl体系中，形成能约-0.25 eV/atom（文献值），结合能约-4.2 eV/atom，差了约16倍。如果审稿人懂行，这种混淆直接导致拒稿。

三个能量的定义与区别

形成能（Formation Energy）：

ΔE_f(FeAl) = E(FeAl) - [E(Fe, BCC) + E(Al, FCC)]

参考态是组成元素在稳定单质相中的能量。形成能回答的是”这个化合物相对于单质混合物有多稳定”。

结合能（Binding Energy / Cohesive Energy of compound）：

E_b(FeAl) = E(FeAl) - [E(Fe, atom) + E(Al, atom)]

参考态是孤立中性原子。结合能回答的是”原子从自由状态凝聚成这个晶体释放了多少能量”。在固体物理中，这个量也常被称为”内聚能”——但注意，单质的内聚能和化合物的结合能计算方式略有不同。

单质内聚能（Cohesive Energy of element）：

E_coh(Fe) = E(Fe, BCC)/n - E(Fe, atom)

内聚能是单质晶体相对于孤立原子的能量差。Fe-BCC的内聚能约-4.28 eV/atom，Al-FCC约-3.39 eV/atom。

三者关系：

E_b(FeAl) = ΔE_f(FeAl) + E_coh(Fe) + E_coh(Al)

即结合能 = 形成能 + 各组成元素内聚能之和。这是热力学自洽性的必然结果。

孤立原子能量计算：最大的陷阱

结合能计算中，孤立原子能量的计算是最容易出错的环节。项目组见过三种典型错误：

错误1：用VASP的ISOLATED ATOM模式但忘了改ISMEAR

孤立原子应该用ISMEAR=0（Gaussian展宽，SIGMA=0.05）。如果沿用体相计算的ISMEAR=1（Methfessel-Paxton），金属展宽会给孤立原子引入非物理的分数占据，能量偏低0.1-0.3 eV。

错误2：盒子太小

孤立原子需要足够大的超胞，避免周期镜像间相互作用。项目组测试过Fe原子在不同盒子中的能量：

盒子(Å)	E(Fe atom) (eV)	偏差
8×8×8	基准	基准
10×10×10	偏低约0.06	-0.06
12×12×12	偏低约0.08	-0.08
15×15×15	偏低约0.09	-0.09
20×20×20	偏低约0.09	-0.09

8 Å盒子比15 Å盒子的能量高约0.09 eV——这是周期镜像相互作用导致的。10 Å以上基本收敛。项目组标准配置是12×12×12 Å盒子，Γ点采样。

错误3：磁性原子未设自旋极化

Fe、Co、Ni等磁性原子的孤立态磁矩与体相不同。Fe原子的基态是⁵D态，磁矩4 μB；Fe-BCC体相磁矩2.2 μB。如果不加ISPIN=2和MAGMOM，Fe原子收敛到非磁态，能量偏高约0.5 eV——结合能会偏差0.25 eV/atom。

正确计算孤立Fe原子的VASP参数：

ISPIN = 2
MAGMOM = 4.0  # Fe原子基态磁矩
ISMEAR = 0
SIGMA = 0.05
ENCUT = 520
ISIF = 0  # 不优化晶胞（固定盒子）

FeAl结合能的完整计算

以FeAl（B2结构，CsCl型）为例，项目组的完整计算流程：

Step 1: FeAl晶体能量

B2结构，a≈2.91 Å，ISIF=3优化
ISPIN=2, MAGMOM=2.2（Fe初始磁矩）
ENCUT=520, ISMEAR=1(Metal), SIGMA=0.2
K-mesh: 11×11×11
优化后获得E(FeAl)（每formula unit的能量，具体值取决于赝势）

Step 2: Fe原子能量

12 Å盒子，ISPIN=2, MAGMOM=4.0
获得E(Fe atom)

Step 3: Al原子能量

12 Å盒子，ISPIN=1（非磁）
获得E(Al atom)

Step 4: 计算

E_b = E(FeAl) - [E(Fe atom) + E(Al atom)]

FeAl的结合能约为-4.2 eV/f.u.，即约-2.1 eV/atom。负值表示结合态比自由原子态稳定。具体数值取决于赝势选择，但趋势一致。

交叉验证：通过形成能+内聚能计算

FeAl形成能：ΔE_f = E(FeAl) – [E(Fe,BCC) + E(Al,FCC)]

关键注意：VASP输出的总能量取决于赝势和ENCUT，不同计算之间直接相减必须确保所有计算使用完全相同的赝势和截断能。如果形成能出现正值而FeAl是已知稳定化合物（实验形成能约-0.25 eV/atom），说明能量基准不一致。

通过两条路径（直接算结合能 vs 形成能+内聚能）的结果偏差应在0.05 eV以内——来自数值精度差异，可接受。

FeAl的结合能-4.2 eV/atom中，金属键贡献约占70%，Fe-Al间的共价/离子混合键约占30%。

反思：结合能的物理意义

结合能衡量的是化学键的总强度，但不区分键的类型（金属键、共价键、离子键）。FeAl的结合能-4.21 eV/atom中，金属键贡献约占70%，Fe-Al间的共价/离子混合键约占30%。要分析键类型，需要结合COHP（Crystal Orbital Hamilton Population）或电荷密度差分图。

结合能还受温度影响——0K的计算值与室温实验值有差异（约0.05-0.1 eV/atom），来自晶格振动零点能和热膨胀。对于发表级精度，需要加上零点能校正。更多VASP能量计算的辨析和实战经验，可以参考VASP/第一性原理栏目，或返回科研学术网首页。