VASP 输出里同时出现 free energy TOTEN、energy without entropy E0 和 entropy T*S，吸附能结果常在这些量之间发生混用。吸附能比较首先是能量口径比较，不能仅取 OUTCAR 或 OSZICAR 末尾出现的某个能量项直接相减。对金属表面，ISMEAR 和 SIGMA 会引入电子熵项；对分子、半导体和绝缘体，展宽项通常较小，但仍然要确认同一组结构采用同一读取方式。
OSZICAR 里有哪些能量？
OSZICAR 每个电子步会给出 F、E0 和 d E 等量。F 对应包含电子熵修正的自由能，OUTCAR 中常写为 free energy TOTEN；E0 是去掉电子熵项后的能量，常写作 energy without entropy。金属体系使用有限展宽后，F 和 E0 的差异会随 SIGMA、态密度和收敛状态变化。
吸附能常写为 Eads = Eslab+ads – Eslab – Eadsorbate。同一条公式里的三项必须来自同一种能量定义。若吸附态取 TOTEN，清洁表面取 E0，孤立分子又取另一个静态输出，最终数值会混入电子熵项差异。图1 的金属表面吸附能基准显示，不同吸附物和表面之间本来就有较大能量差，读取口径混乱会让误差来源难以分辨。
图1 Ni(111) 与 Pd(111) 上 MeOH、CH₃O 和 CO 吸附能的计算值与实验值对照。该图用于说明吸附能比较必须先统一能量定义。 DOI: 10.1038/s41467-022-34507-y。
电子熵项来自哪里？
金属表面附近存在连续能带，有限温度展宽用于帮助电子占据收敛。VASP 中 entropy T*S 反映展宽占据带来的熵项，F 与 E0 的差可以理解为电子自由能与去熵能量之间的差别。SIGMA 越大，或费米能级附近态密度越高，二者差异通常越明显。
计算吸附能时，电子熵项属于展宽占据带来的能量组成，其影响取决于不同端点之间能否近似抵消。若清洁表面、吸附态和分子参照的 entropy T*S 差异处在吸附能差值的同一量级，直接比较 TOTEN 会把展宽贡献带进吸附强弱排序。金属表面可同步记录 F、E0 和 entropy T*S，并用较小 SIGMA 或 tetrahedron 静态对照检查吸附能排序。
金属表面优先检查什么？
金属 slab、合金表面和强吸附体系常用 ISMEAR = 1 或 2 做结构优化，静态能量可再用更合适的展宽设置检查。若目标是比较同一金属表面上 CO 顶位、桥位和三重位，三种构型的 k 点、ENCUT、ISMEAR、SIGMA、LDIPOL、固定层和收敛阈值都应一致。
若研究对象是合金表面，还要核对各构型的表层元素排列和近邻配位是否保持同一比较口径。表面偏析、吸附诱导重构和磁矩翻转都会改变费米能级附近态密度，进而改变电子熵项。结果记录同步写入表层组成、近邻配位、磁矩状态和所用能量项，参数影响与结构差异才能分开判断。
图2 和图3 给出吸附能误差在不同泛函与吸附类型中的差异。这类对比显示，能量读取口径和计算设置会共同传递到吸附能误差。若某两种构型的 Eads 只差几 meV 到十几 meV，而电子熵项差异也在同一量级，结论应回到更严格的静态对照，而不是只报告一个排序。
图2 PBE+D3 吸附能误差与修正结果对照。图中化学吸附和 vdW 主导吸附分区提醒能量差会受模型与泛函共同影响。 DOI: 10.1038/s41467-022-34507-y。
图3 不同周期泛函及修正方案的吸附能平均绝对误差。该图说明吸附能口径和泛函选择会一起影响最终排序。 DOI: 10.1038/s41467-022-34507-y。
分子参照怎样处理？
孤立分子通常放在足够大的盒子中，k 点取 Γ 点，展宽可用适合分子体系的设置。表面吸附态则是周期 slab，k 点和展宽与金属态有关。由于孤立分子和金属表面电子结构差异很大，吸附能公式中的参照项更要明确能量类型。
气相分子还要单独检查自旋和盒子尺寸。O₂、NO、NO₂ 这类开壳层分子若自旋态设错，参照能会直接传入吸附能；CO、H₂O、NH₃ 等闭壳层分子也要避免分子之间周期镜像相互作用。若分子参照和表面端点采用不同能量口径，后续 ZPE 或熵修正也会失去清楚的起点。
结果记录可保留 EF、E0、T*S 和 TOTEN 四类量。对金属吸附能，若优化阶段使用较大 SIGMA，静态阶段可降低 SIGMA 并读取同一口径的 E0；若体系为半导体或绝缘体，仍应检查 F 与 E0 差值是否可以忽略。可比性来自同一规则下的端点能量，不是来自单个结构的总能精度。
SIGMA 改变排序怎么办？
若 SIGMA 从 0.20 eV 改为 0.05 eV 后，两个吸附构型排序发生变化，说明排序对电子展宽较敏感。此时先检查电子步收敛、k 点密度和磁性初值，再比较 E0 口径下的吸附能。对于费米能级附近存在尖峰 DOS 的合金或缺陷金属表面，k 点加密往往比单纯缩小 SIGMA 更有意义。
还有一类情况来自吸附前后金属性改变。若分子吸附使表面从半金属态变成明显开隙态，F 和 E0 差异在清洁表面和吸附态之间不会相同。此时应把 DOS、占据数和电子熵项一起列出，并说明吸附能使用的读取口径。图4 的周期模型与簇模型对应关系说明，能量差还会受到模型大小和截取方式影响，不能把所有误差都归到电子熵项。
图4 O/Pt(111) 周期模型与簇模型的对应关系。该图用于说明总能读取只是一环，模型截取和参照态同样影响吸附能。 DOI: 10.1038/s41467-022-34507-y。
小能量差如何报告？
吸附能差很小时，结构热振动、覆盖度、真空层、偶极修正和色散修正都可能改变排序。若 CO 在两个邻近位点的 Eads 只差 0.02 eV，而频率校正或展宽口径也在相近范围，结果更适合作为“能量接近”的判断。
TOTEN 与 E0 的选择要服务于同一组比较。优化看收敛过程时可以关注 TOTEN 和力，最终吸附能结果应说明读取 E0 还是 TOTEN，并提供展宽参数。图5 中 CO 吸附能随簇尺寸变化，说明小差值排序还会受模型尺度影响；对于周期 slab，类似问题体现为层厚、k 点、真空层和展宽的联合收敛。
图5 CO 在 Cu 簇上的吸附能随簇尺寸变化。该图用于提示局域模型、周期模型和能量项定义之间的可比性边界。 DOI: 10.1038/s41467-022-34507-y。
结果记录保留哪些量？
每个端点记录结构名、ISMEAR、SIGMA、ENCUT、KPOINTS、电子步收敛、TOTEN、E0、entropy T*S 和最终磁矩。若做频率或自由能校正，再另列 ZPE、TS 与 G。吸附能差异来自结构本身、电子熵项、磁性状态还是后处理口径，可通过这些量逐项排除。
计算时先用统一参数优化清洁表面、吸附态和孤立分子，再对代表性构型做静态能量对照，最后用同一能量口径代入 Eads。金属表面可测试 ISMEAR = 1 或 2 的优化与更小 SIGMA 静态结果，半导体表面可测试 ISMEAR = 0 或 tetrahedron 类设置。具体数值阈值以体系收敛测试为准。
结果边界怎样说明？
当 F 和 E0 差异远小于吸附能差，并且改变 SIGMA 不改变排序时，可以说明该组吸附能对电子熵项不敏感。当 F 和 E0 差异接近排序差值时，应把展宽敏感性作为结果边界，并补充 DOS 或 k 点对照。
以 CO/Ni(111) 和 CO/Pd(111) 的比较为例，如果两套表面的金属性不同，entropy T*S 在清洁表面和吸附态中的抵消程度也可能不同。此时可以先固定同一 SIGMA 得到主结果，再用较小 SIGMA 做静态对照；若两次 Eads 差值保持在预设收敛范围内，吸附能排序才更适合归因于表面差异。
对于含磁性元素的表面，还要把磁矩变化纳入同一组结果记录。吸附物可能改变局域磁矩和费米能级附近态密度，进而改变电子熵项大小。若某个构型的磁矩状态与其他构型不同，先比较磁构型能量和局域磁矩，再进入 TOTEN 或 E0 读取口径。
吸附能的可靠性来自端点、参照和能量读取的统一。TOTEN 和 E0 没有脱离场景的唯一答案；对金属表面，E0 常用于最终能量比较，TOTEN 可用于观察电子自由能收敛。需要交代的对象包括同一批结构是否使用同一种口径、展宽项是否足够小、排序是否经静态对照保持一致。
要点：TOTEN 包含电子熵项，E0 对应去掉电子熵项后的能量，二者不能在同一吸附能公式中混用。
要点：金属表面要同时记录 ISMEAR、SIGMA、F、E0 和 entropy T*S，并用静态对照检查排序。
要点：孤立分子、清洁 slab 和吸附态应使用统一的能量读取口径。
要点：当吸附能差接近电子熵项差异时，排序应回到 k 点、SIGMA、DOS 和磁性状态共同检查。
