结合能计算全流程：从能量差定义到参考碎片选择、几何口径、归一化处理与误差分析

说明：本文华算科技主要介绍结合能的能量差定义、参考碎片选择、几何口径、归一化方式和误差来源。

结合能的能量差怎样定义？

复合体系 AB 沿指定解离过程分成 A 与 B，复合态和分离态的电子总能差构成结合能。按放热过程记负值时，常用定义为 ΔE_bind=E_AB−E_A−E_B。ΔE_bind

另一套文献写法为 E_b=E_A+E_B−E_AB，稳定结合态给出正值。两套公式只差一个负号，物理口径相同。数值表若省略公式与符号约定，“−2.0 eV”和“+2.0 eV”可能描述同一解离过程；标题或列名里的 binding energy 无法替代公式。

分子复合物以两个单体为碎片 A、B，负载团簇以团簇与载体为碎片，二维异质结构以上下两个片层为碎片。原子、分子、晶格片层和带电碎片代表不同解离反应，即使元素组成相同，所得能量也属于不同物理过程。

距离扫描把分离态参考画成一条能量曲线：短距离区域由电子云重叠产生强排斥，曲线最低点给出平衡间距，远距离极限逐渐接近零。势阱深度与选定碎片的电子解离能相连，最低点位置则是该模型和计算方法下的平衡几何。

层间距扫描还包含两种几何协议。固定片层内部坐标、仅改变 d，所得曲线描述垂直分离坐标上的单点能；每个 d 都放松片层内部坐标，曲线会混入面内应变、起伏和局域重构的松弛能。平衡结合能应从预先声明的同一种扫描协议取得。

上述差值通常指 0 K 电子能。零点振动、有限温热修正和熵项尚未写入 ΔE_bind。讨论溶液缔合、气相二聚或温度依赖平衡时，应另算 ΔG_bind=ΔE_elec+ΔZPE+ΔH_thermal−TΔS，并写明标准态。

总体系/参考碎片保持同一计算口径？

E_AB、E_A 与 E_B 应采用同一泛函、赝势或基组、色散方案和数值收敛标准。元素赝势版本、U 值或相对论处理发生变化时，总能零点也随之变化。周期模型中的碎片计算宜保留与复合体系相同的超胞和面内晶格，借此抵消网格与周期尺寸误差。

周期平面波模型可从优化后的 AB 超胞复制两份：一份删除 B，只保留 A 在复合物中的坐标；另一份删除 A，只保留 B。晶格矢量、k 点网格、截断能和电子收敛阈值保持不变，三次单点能便构成冻结几何相互作用能。若目标改为自由碎片结合能，再分别优化 A、B 并计入形变项。

把 A、B 固定在复合物中的几何 R^AB，可写 ΔE_int=E_AB(R^AB)−E_A(R_A^AB)−E_B(R_B^AB)。该差值聚焦片段在既定几何下的电子相互作用，分子内键长改变、片层应变和载体重构未计入碎片松弛收益。

若 E_A、E_B 来自各自自由优化结构 R⁰，ΔE_bind 会包含形变代价。负值口径下有 ΔE_bind=ΔE_int+ΔE_def,A+ΔE_def,B，其中 ΔE_def,A=E_A(R_A^AB)−E_A(R_A⁰)。柔性分子、应变异质结和明显重构的载体常出现可观的形变项。

孤立原子参考应使用正确的基态自旋与占据，尤其是过渡金属、O、N 和自由基。带电复合物的碎片电荷之和应等于 AB 总电荷；不同电荷分配代表不同解离产物组合。周期带电体系还会引入均匀补偿背景和有限尺寸校正，相关处理必须在三个能量项中维持同一方案。

金属体系常使用电子展宽。三个能量项应读取同一类能量，例如均采用外推到 0 K 的 E₀，避免把含熵自由能 F 与 E₀ 混用。孤立分子或原子放在大真空盒中时，Γ 点、对称性、偶极修正和自旋极化也应明确；真空盒扩大会改变周期镜像之间的残余作用。

分子/团簇/界面的结合能怎样归一化？

分子二聚体通常报告 eV/复合物或 kJ·mol⁻¹；团簇常报告整个团簇的总结合能或每原子平均值；固体界面常报告 eV/界面、J·m⁻² 或 meV·Å⁻²。单位和除数共同限定数值含义，跨体系比较应统一到同一解离对象和归一化单位。

由 n 个同种原子组成的团簇采用正值口径时，平均结合能可写为 E_b(n)=[nE_atom−E_cluster(n)]/n。除以 n 后可比较不同尺寸的平均凝聚程度；总解离能则不除以 n。多组分团簇把 nE_atom 改为各元素孤立原子能之和。

每原子平均结合能随 n 增大可逐渐靠近体相内聚能，特定 n 的碎裂倾向则由 E_frag=E_n−m+E_m−E_n 描述。平均凝聚程度与指定碎裂反应属于两种能量问题，幻数团簇常在二阶能量差或碎裂能中出现局部异常。

界面负值口径可写为 γ_bind=(E_AB−E_A−E_B)/A。周期超胞若含两个等价界面，分母改为 2A。面内应变必须在参考片层中保留；把 A、B 各自放回自由晶格常数会额外计入应变释放，使界面电子作用与晶格失配能量叠加。

界面工程也常报告正值分离功 W_sep=(E_A+E_B−E_AB)/A。理想脆性分离且碎片几何冻结时，W_sep 与负值 γ_bind 只差符号。分离后表面重构、化学计量改变或原子转移会改写末态，数值转为指定断裂反应的功。

层状固体也常把结合能除以原子数。AB 堆叠石墨的纵轴采用 meV/atom，界面研究则常用 meV·Å⁻²。两种单位经过原胞面积和原子数换算后才可比较，原胞选取、层数与界面数都应写入换算关系。

哪些计算细节会改写结合能数值？

局域原子轨道基组中，片段 A 进入复合物后可借用 B 上的基函数，E_AB 获得额外的数值降低，结合能随之偏向更强。Boys–Bernardi counterpoise 让单体在完整复合物基组中计算：ΔE_int^CP=E_AB^AB−E_A^AB−E_B^AB，上标 AB 表示保留另一片段的 ghost 基函数。

CP 值仍受基组不完备误差约束，小基组中两类误差可能部分抵消。报告分子结合能时，基组、CP 口径和单体几何应成套出现。平面波基组没有同形式的片段借基问题，但截断能、真空盒大小和周期镜像仍会进入能量差。

层间作用和分子缔合含有长程色散，半局域泛函对非局域相关的描述有限。图2中 PBE 与 PBE+vdW-DF2 给出不同的势阱和远程衰减，表明色散方案属于结合能定义的一部分。不同修正方案之间的数值差异应通过同一结构集和同一收敛标准评估。

结合能是多个大总能相减后得到的小差值，总能各自“看似收敛”仍可能留下不可接受的差额。k 点、截断能、真空厚度和超胞尺寸的测试应直接监测 ΔE_bind 的变化量。例如单项总能变化 5 meV，而三个能量项的误差未相互抵消，所得结合能可出现同量级偏移。

孤立原子的自旋基态、界面面内应变、分子盒尺寸、slab 真空层、偶极修正、k 点和截断能都会移动 E_AB−E_A−E_B。结果表应同时标注公式符号、碎片几何、电荷与自旋、归一化单位、泛函及色散方案；其中任一计算条件改变，都可能让 E_b 产生可见偏移。