吸附能中的形变能怎么算？

引言

吸附能很负，并不等于吸附物与表面的直接相互作用同样强。把吸附能分成形变能与相互作用能，才能区分“成键获得多少能量”和“为了形成该构型付出多少结构代价”。分子弯曲、键长改变、表面原子抬升或局域重构所需的能量，已经包含在常规 E_ads 中。

这套分解适用于分子吸附、团簇负载、强重构表面和柔性二维材料。它不要求修改 VASP 源码，关键是从优化后的复合结构中提取冻结片段，并保持晶胞、原子坐标、POTCAR、ENCUT、k 点、泛函、色散修正和自旋设置一致。

图1 Cu(111)–有机分子界面的研究对象与吸附构型总览。分子骨架、吸附高度和表面响应共同进入吸附能分解。 DOI: 10.1002/cplu.202400771。

总吸附能为什么看不见形变代价？

常规吸附能比较复合体系、平衡表面和孤立吸附物三个状态。若定义吸附放热为负，E_ads = E_AB – E_A – E_B。E_AB 来自优化后的吸附体系；E_A 与 E_B 通常分别来自各自平衡构型。这个差值把结构形变与片段成键合并在一个数中。

E_ads = E_AB – E_A^eq – E_B^eq　　式(1)

分子从气相平衡构型变成吸附态几何时，内部键和二面角会改变。形变能对应两个片段各自从平衡构型走到吸附态冻结构型的能量升高。表面从干净 slab 变成吸附后的冻结几何时，顶层原子的位移也会进入这项能量。

吸附能分解示意把准备能、相互作用能、吸附高度和电荷转移放在同一个界面框架中。对平面芳香分子，色散可能带来稳定作用，同时 Pauli 排斥抬高电子相互作用项；对明显弯曲的化学吸附构型，分子准备能和表面准备能都会增加。

若两个位点的 E_ads 只差很小，却对应完全不同的表面重构，单看总能会掩盖差异。一个位点可能具有更强的片段相互作用，同时承担更大的表面形变；另一个位点相互作用较弱，但结构代价也较小。相近的总吸附能可以来自不同的能量补偿关系。

图2 吸附能分解及几何、电子界面量的关系示意。准备能对应吸附物和表面从平衡构型变到吸附构型的代价。 DOI: 10.1002/cplu.202400771。

四个冻结片段能量怎样定义？

从优化后的 AB 结构复制一份 POSCAR，删除 B 只保留 A，得到 A 在吸附态几何下的冻结片段 A*；再复制一份并删除 A，得到 B*。两个冻结片段都不做离子优化，只进行静态单点。A* 与 B* 的坐标、晶胞和原始 AB 结构完全对应。

吸附物形变能写成 E_def,A = E_A* – E_A^eq，表面形变能写成 E_def,B = E_B* – E_B^eq。这里的 eq 表示各片段独立优化后的平衡状态，星号表示从复合体系抽出的冻结几何。

E_def = (E_A*-E_A^eq) + (E_B*-E_B^eq)　　式(2)

片段相互作用能由 E_int = E_AB – E_A* – E_B* 给出。E_int 比较的是同一吸附态几何下“合在一起”和“拆成冻结片段”的能量差。片段从平衡构型发生形变的代价已经单独放入 E_def。

E_int = E_AB – E_A* – E_B*；　E_ads = E_int + E_def　　式(3)

式(3)也是最直接的闭合检查。若 E_int + E_def 与原始 E_ads 不一致，通常说明四个单点没有使用同一电子结构口径，或冻结片段的坐标、晶胞和电子数处理发生变化。数值差应落在输出精度与电子收敛误差允许范围内。

对带电吸附物或电化学界面，删除片段后总电荷怎样分配并不唯一。此时的冻结片段分解依赖人为选定的电荷态，分量只能在同一电荷拆分方案内比较。若无法给出物理上明确的片段电荷，优先报告中性片段分解并说明其边界。

图3 多种有机分子在 Cu(111) 上的吸附能及其色散、电子相互作用和准备能分量。总吸附能接近时，分量仍可能明显不同。 DOI: 10.1002/cplu.202400771。

VASP 单点怎样保持同一计算口径？

AB、A*、B* 三个静态单点要使用相同晶胞、ENCUT、KPOINTS、POTCAR 版本、交换相关泛函、IVDW、LASPH、自旋极化和收敛标准。孤立分子 A^eq 若放在单独盒子里计算，需要确认盒子尺寸与 Γ 点采样已经收敛；为减少基准差，也可以把分子放进与 slab 相同的晶胞中再做平移检查。

冻结片段计算设置 NSW = 0 或 IBRION = -1，只读取当前结构的静态能。形变能常是多个大数相减，任何输入口径变化都会直接污染分量。金属表面应统一 ISMEAR、SIGMA 和能量读取口径；若最终比较使用 energy without entropy E0，所有构型都读取最后一个完成静态步对应的 E0。

从 AB 删除吸附物后，B* 的表面原子不能回到干净 slab 的平衡位置；否则得到的是另一种表面状态，不再是表面形变能。类似地，A* 不能在真空中重新优化。冻结片段的任务就是保留吸附造成的几何变化，只让电子密度重新自洽。

自旋体系还要检查 MAGMOM 与最终磁矩。删除吸附物后，片段的电子态可能从 AB 中继承到另一局域极小值；可分别用合理的初始磁矩进行几次静态单点，选取与目标片段电荷和自旋定义一致的解。若 A*、B* 的自旋态不清楚，形变与相互作用分量也会带有电子态选择差异。

能量分量与吸附高度的关系说明，总 E_ads 的变化可来自色散、电子相互作用和准备能共同调整。VASP 的 DFT-D3 修正为加和项时，可进一步读取色散贡献；采用非局域 vdW 泛函时，不能直接照搬 D3 的分项口径。能量分解的层级应与所用泛函能够提供的物理量对应。

图4 吸附高度变化时吸附能及分量的响应。冻结几何单点可把结构形变代价与片段间相互作用分开。 DOI: 10.1002/cplu.202400771。

哪些对比最能说明形变来自哪里？

对柔性分子吸附，可把 E_def,A 与关键二面角、吸附高度和分子内键长变化对应。若某位点总吸附更稳定，却伴随更大的分子形变能，说明片段相互作用足以补偿构象代价；这与“分子越接近平衡构型，吸附越强”并不等价。

对氧化物或二维材料，表面形变能可能成为主要项。若吸附后金属原子抬升、配位数改变或局部晶格翘曲，应把 E_def,B 与位移、键角和层间距离一起比较。准备能分布说明，化学吸附体系往往伴随更明显的几何调整，但准备能大小不能单独定义成键类型。

覆盖度比较需要使用相同表面面积和一致的归一化单位。片段划分改变时，E_int 与 E_def 的含义也随之改变。一个大分子占据多个位点时，按每分子与按每活性位归一化会得到不同量级；共吸附体系还要把所有共吸附物纳入 A 片段，或明确拆成多个片段。

若吸附后发生断键、质子转移或表面原子被抽出，原来的 A/B 身份已经改变。此时仍能做冻结片段单点，但它回答的是指定碎片划分下的能量关系，不宜再称为完整分子的普通吸附形变能。应把反应终态按新的化学计量重新定义，并与反应能或形成能一起分析。

完整结果至少包含 E_ads、E_def,A、E_def,B、E_int、闭合误差、片段电荷自旋、主要位移和关键键长。当强相互作用与大形变相互抵消时，分解结果比单个吸附能更能解释位点差异；当形变项接近零时，总吸附能主要反映冻结片段间的相互作用。

参照构型还会影响形变能的物理含义。柔性分子可能有多个气相构象，表面也可能存在几个接近的重构态。若 A^eq 选全局最低气相构象，E_def,A 包含构象转变与局域弯曲；若选与吸附态同构的局部极小值，它只衡量成键附近的附加形变。两种参照都可计算，但不能混在同一系列位点排序中。

闭合误差应与电子收敛精度、k 点和截断能测试一起报告。若分量只有几百分之一电子伏特的差异，先把 EDIFF、精确度和静态网格收紧，再判断位点差别；若闭合误差随真空层或超胞变化，说明冻结片段受到周期镜像或偶极相互作用影响。这样的数值检查能把真正的结构代价与计算噪声分开。

图5 化学吸附与物理吸附体系的准备能分布。较大的准备能说明成键前伴随更明显的分子或表面形变。 DOI: 10.1002/cplu.202400771。

本章要点总结

要点：总吸附能同时包含片段相互作用和吸附物、表面的形变代价。

要点：冻结片段必须从同一个优化复合结构中提取，并保持同一晶胞和电子结构设置。

要点：Eads = Eint + Edef 可用于检查分解是否闭合，带电体系还需说明片段电荷口径。

要点：形变能要和键长、位移、二面角及表面重构共同解释。