单原子模型建好以后，很多人第一反应是算一个结合能：越负越容易锚定，似乎就能证明金属原子站得住。真正麻烦的地方在于，结合能的“参照物”有好几种：孤立金属原子、体相金属、金属前驱体、干净载体、缺陷载体。参照能量不同，结合能回答的问题就不同。
如果只想看一个金属原子是否愿意落到某个锚定位点，可以用孤立原子参照；如果要讨论单原子会不会向金属团簇转化，就必须把体相金属或团簇参照放进来；如果载体本身带空位、N 配位孔或羟基，还要先把载体参照写清楚。
图1 单原子、纳米颗粒以及动态聚集/分散过程的自由能示意。DOI: 10.1002/advs.202308046。
孤立原子参照看锚定强弱
最常见的单原子结合能，是把金属原子从孤立态放到载体位点上，看体系能量降低多少。这个定义适合比较同一载体上不同金属，或者同一金属在不同锚定位点上的相对偏好。前提是孤立原子自旋、多重度、盒子尺寸和赝势口径都一致。
Ebind = EM/S − ES − EM(atom)
这里 EM/S 是金属单原子负载后的总能，ES 是对应载体的总能，EM(atom) 是孤立金属原子的能量。按这个写法，Ebind 越负，说明从孤立原子到锚定位点的放热越强。这个数值说明锚定过程有多强，不等于已经排除了团聚趋势。
图2 金属-杂原子掺杂碳载体单原子位点的结构矩阵与筛选流程。DOI: 10.3389/fchem.2022.873609。
为什么同一个符号会被写成不同含义？
有些论文把上式叫 binding energy，有些叫 adsorption energy，也有些把金属化学势换成体相金属后称为 formation energy。符号不统一并不一定错，真正需要看的是减掉了什么。若减的是孤立原子，讨论的是原子落位；若减的是体相金属每原子能量，讨论的是相对金属聚集态的分散偏好；若减的是金属盐或前驱体化学势，则更接近合成热力学口径。
因此，不能只截取一个 Ebind 数字去跨论文比较。不同泛函、不同 PAW 赝势、不同孤立原子自旋、不同载体缺陷都会改变绝对值。更可用的做法，是在同一套 VASP 设置下，把一组金属或一组位点放进同一张比较里。
μM 决定形成能的物理口径
如果希望把单原子锚定写得更一般，可以用金属化学势 μM。它可以取孤立原子能量，也可以取体相金属每原子能量，还可以按实验环境近似为前驱体中的金属化学势。不同选择对应不同问题，不应混在一列里排序。
Eform = EM/S − ES − μM
当 μM = EM(atom) 时，Eform 与孤立原子锚定能等价；当 μM = EM(bulk) 时，它更接近“从金属相中取出一个原子放到载体上是否有利”。金属参照从孤立原子换成体相金属后，数值会整体抬高，排序也可能改变。
图3 过渡金属单原子嵌入杂原子掺杂碳结构后的形成能对比。DOI: 10.3389/fchem.2022.873609。
孤立原子计算要留意什么？
孤立金属原子通常放在足够大的盒子里，只用 Γ 点，并打开自旋极化。过渡金属原子有多个可能自旋态，MAGMOM 初值会影响收敛到的能量。若这一步没有统一，后面所有 Ebind 都会带着系统偏差。对 Fe、Co、Ni、Cu 这类元素，建议保留孤立原子的磁矩、能量和 INCAR 设置，确保整组金属用同一口径。
体相金属参照也要单独收敛。晶格常数、k 点密度、ISIF 设置和磁性状态都会影响 EM(bulk)。比较多个金属时，不能把数据库值、不同泛函值和自己算的单原子体系混用；否则看似是团聚趋势，实际可能只是参照能量来源不一致。
干净载体和缺陷载体不能随手替换
单原子并不是总落在完美表面。常见锚定位点包括 N4 孔、N3 缺陷、氧空位、硫空位、羟基邻位和台阶位。若金属嵌入前需要先挖孔或造空位，ES 就应该是“已经挖好孔、但还没有放金属”的载体，而不是完美载体。载体参照要跟锚定位点保持同一结构身份。
Eanchor = EM/S* − ES* − EM(atom)
这里 S* 表示已经具备锚定位点的载体，例如带 N4 孔的石墨烯、带氧空位的氧化物或含缺陷边缘的二维材料。这样写，Eanchor 只衡量金属落入这个位点的锚定强弱；若直接用完美载体作参照，空位形成能、掺杂形成能和金属锚定能会被揉在一起。
什么时候要把缺陷形成能也算进去？
如果问题是“这种单原子位点整体是否容易形成”，就不能只看 Eanchor。例如 M-N4 位点要先形成 N 配位孔，氧化物单原子位点要先生成氧空位，这些步骤本身也有能量成本。此时应把缺陷形成能和金属锚定能拆开报告，再根据研究对象决定是否合并。
一个更清楚的做法是并列写两句话：第一句说明缺陷或配位孔的形成代价，第二句说明金属原子进入该位点后的锚定收益。若锚定很强，但缺陷形成能很高，那么该位点在热力学上并不一定容易大量形成；若缺陷已经由实验处理、等离子体或热退火引入，锚定能就更适合作为后续分散指标。
体相金属参照能把团聚倾向放进来
单原子是否会团聚，取决于金属原子更愿意留在载体位点，还是回到金属-金属相互作用更强的聚集态。孤立原子参照往往会让锚定看起来很强，因为金属原子从气相落到任何表面都可能显著放热。要讨论团聚，必须把金属内聚能或体相金属每原子能量作为对照。
ΔEagg = EM/S* − ES* − EM(bulk)
若 ΔEagg 小于 0，说明相对于体相金属每原子能量，单原子位点仍有热力学优势；若 ΔEagg 大于 0，则金属原子在热力学上可能更偏向金属聚集态。这个判断还需要迁移能垒、工作温度、覆盖度和反应气氛配合。团聚趋势不是只看一个结合能，而是看锚定收益能否抵消金属-金属内聚驱动力。
图4 内聚能差异驱动多金属原子分散或选择性聚集的机制示意。DOI: 10.1038/s41467-024-52520-1。
为什么还要看迁移和覆盖度？
热力学上偏向团聚，不代表计算模型一优化就一定形成团簇；热力学上偏向单原子，也不代表高温反应下完全不会迁移。单原子在载体上的迁移能垒、相邻金属原子距离、载体缺陷密度、反应中间体覆盖度都会改变真实分散状态。尤其在氧化物和碳载体上，吸附的 O、OH、CO、N2 等中间体可能改变金属位点的配位环境，使团聚趋势随反应条件变化。
所以更可靠的比较不是只排一个 Ebind，而是同时给出三类量：孤立原子参照的锚定强弱、体相金属参照的抗团聚趋势，以及必要时的迁移能垒。若三者都支持单原子分散，结论才更完整；若只有锚定能很负，而 ΔEagg 为正，就要把“可锚定”和“抗团聚”分开写。
图5 不同金属元素内聚能与原子分散/聚集构型之间的相关关系。DOI: 10.1038/s41467-024-52520-1。
先固定同一组计算边界
同一批单原子结合能比较中，泛函、PAW 赝势、截断能、k 点、展宽、ISPIN、MAGMOM、vdW 修正、DFT+U 口径都要保持一致。载体模型也要统一超胞尺寸、缺陷浓度、真空层、固定层和偶极修正。否则不同金属之间的差异，可能来自计算边界，而不是金属-载体相互作用。
比较推荐的任务组合是：先优化 S*，再优化 M/S*，然后分别计算孤立金属原子和体相金属每原子能量。若需要讨论工作态，还可以在相同 M/S* 位点上加入 *H、*OH、*CO 或 *N2，观察吸附后金属原子是否明显偏离锚定位点。
结果表述要把三个参照写在一起
结果部分可以把同一金属位点写成三行或三列：Ebind 对孤立原子，Eform 对选定 μM，ΔEagg 对体相金属。再配合 M-N、M-O 或 M-C 键长、Bader 电荷、PDOS 和差分电荷密度，说明锚定来自配位键、缺陷电荷补偿还是金属-载体杂化。只有参照能量、几何结构和电子结构同时交代清楚，单原子分散能力比较才有可比性。
写结论时可以用更具体的句子：在相同 N4 载体和相同 VASP 设置下，M 原子相对孤立原子具有较负 Ebind，说明锚定位点对 M 有强吸引；相对体相金属的 ΔEagg 仍为负，说明该位点在热力学上有利于保持单原子分散；Bader 电荷和 M-N 键长进一步支持 M 与配位 N 之间存在明显相互作用。
要点一：孤立金属原子参照主要回答锚定强弱，不能直接等同于抗团聚能力。
要点二：μM 可以取孤立原子、体相金属或前驱体化学势，不同选择不能混排。
要点三：载体参照要与锚定位点一致，N4 孔、氧空位和完美载体不能随手替换。
要点四：ΔEagg 使用体相金属每原子能量，更适合讨论金属团聚趋势。
要点五：完整比较应同时报告参照定义、几何配位、Bader 电荷、PDOS 或差分电荷密度。
