金属单原子、双原子和小团簇经常被放在同一张图里比较，但 M1、M2、M4并不是“金属数量逐渐增加”的简单序列。单原子位点通常没有金属—金属键，双原子开始出现邻位协同，小团簇则可能表现出金属纳米颗粒的雏形。位点尺度一变，吸附物看到的成键对象就变了。
因此，M1、M2、M4的吸附能不能直接混着排。它们的参照态、有效位点数、金属利用率、构型自由度和电荷分布都不同。若不先分清这些口径，吸附能更负可能只是因为团簇提供了更多金属—吸附物接触，而不是单个位点本身更活泼。
单原子位点为什么不是最小团簇？
单原子通常被载体配位固定，例如 M—N₄、M—C₃ 或氧化物缺陷位。它的电子结构强烈依赖载体，金属中心常表现为孤立配位中心。团簇则有金属—金属键，吸附物可以同时接触多个金属原子，局域电子态也更接近小尺寸金属颗粒。
这一区别会直接影响 VASP 建模方式。M1模型通常要先确认金属是否被缺陷、N 配位或空位稳定住；M4模型还要确认团簇内部构型，例如平面四原子、四面体、链状或贴附在载体上的畸变结构。单原子模型问的是一个配位中心，团簇模型问的是一个多原子界面。
因此，M1不能被简单理解为 M4的“少放三个原子”。M1的能级、磁矩和电荷更多由载体配位决定；M4的电子态还会受到金属—金属键、团簇形状和接触面积影响。若用同一套语言解释两者，容易把载体调控和团簇协同混成一件事。
图1 Fe5、Co5、Pt5等金属团簇负载在石墨烯上的结构示例，说明团簇位点具有金属—金属键和多原子吸附界面。DOI：10.1021/acsami.5c08441
双原子和四原子为什么会多出协同吸附？
M2可以让 *OOH、*O 或 *CO 形成桥式吸附，M4则可能出现顶位、桥位、三重位和团簇边位。吸附构型一多，最低能结构就不再代表“某一个金属原子”的性质，而是代表团簇整体的几何与电子协同。团簇吸附能常常包含多中心成键贡献，不能直接拿来和单原子顶位吸附比强弱。
多中心协同有时会带来好处，例如一个金属原子抓住 O 端，另一个金属原子稳定 H 端或另一个 O 端，中间体就不再依赖单个金属位点完成全部成键。但它也可能带来过强吸附，尤其是 *O、*OH 或 *OOH 这类含氧中间体，容易被多个金属同时拉住，导致产物脱附或后续加氢变慢。
做构型枚举时，M2至少要放端位和桥位；M4则要把团簇顶面、边位、团簇—载体界面位都试一遍。若最低能结构总是落在团簇边缘，而不是你想比较的中心位点，就要按最终位点重新命名，不能继续用初始位点描述结果。
图2 单原子、双原子和三原子位点在 ORR 中的 O—O 活化方式示例，展示多金属位点会引入不同电子接受/反馈通道。DOI：10.1021/jacsau.3c00432
参照态要不要随金属数量改变？
比较 M1和 M4时，首先要区分两个问题：一个是“这个结构稳定吗”，另一个是“它吸附某个中间体强不强”。结构稳定性应看金属结合能、团聚倾向或形成能；吸附强度才看 Eads。如果把金属加入载体的能量和吸附物结合能混在一起，结论会失去可比性。
同样是吸附能，M1和 M4的 Esite含义也不同。M1的空位点通常只有一个金属中心和载体配位；M4的空位点包含团簇内部键合、团簇—载体界面和可能的重构空间。若吸附后 M4团簇发生形变，Eads里还会混入一部分团簇重排能。
Eads= Esite+ads− Esite− Eadsorbate
这个式子里，Esite必须对应同一个位点尺度。M1的 Esite是单原子结构，M4的 Esite是团簇结构。它们可以分别计算吸附能，但不能只按数值大小说明 M4的“单个位点”更强。
更清楚的做法，是把表格拆成两层：第一层列结构稳定性，例如金属结合能、最短 M—M 距离、团簇是否保持形状；第二层列吸附结果，例如最低能位点、Eads、关键键长和自由能修正。这样 M1、M2、M4可以被放在同一套研究里，但不会被误当成同一种位点。
图3 不同金属团簇位点的 ORR 自由能路径示例，说明团簇位点常通过多中心吸附改变反应台阶。DOI：10.1021/acsami.5c08441
同一个吸附物为什么在团簇上会有多个身份？
以 *OOH 为例，它在单原子位点上多为端位吸附；在 M2或 M4上可能出现端位、桥位、侧向吸附，甚至 O—O 键被过度拉长。对于 *CO、*OH 和 *O 也类似。吸附构型不同，吸附能反映的是不同成键方式，不宜只看数字。比较不同尺度位点时，必须同时比较吸附位点、键长和中间体取向。
例如 *CO 在单原子上常是 C 端顶位吸附，到了小团簇上可能变成桥式 CO；*OH 在 M1上是 M—O 单键，在 M4上可能被两个金属共同稳定。此时如果只写“*OH 在 M4上吸附更强”，信息是不完整的，因为强吸附可能来自成键数增加，而不是单个金属中心电子结构更合适。
自由能路径也会因此改变。单原子位点常适合逐步比较 *OH、*O、*OOH 的端位吸附；团簇位点可能出现桥式 *O、侧向 *OOH 或 O—O 断裂。若反应路径中的中间体身份已经改变，就不能把 M1和 M4的台阶图简单重叠。
图4 单/双原子位点的轨道与电荷转移示例，说明多原子位点可能通过不同电子通道改变吸附。DOI：10.1021/jacsau.3c00432
先分别确认稳定性，再做吸附
建议先为 M1、M2、M4分别计算结构稳定性：金属结合能、相邻金属距离、是否迁移、是否塌陷到缺陷位。只有稳定结构确认后，再放入 *H、*OH、*OOH、*CO 或目标中间体。结构都没站住时，吸附能再漂亮也很难说明位点可用。
结构优化后要保存几何指标。M1看金属—载体键长和金属是否离开锚定位；M2看 M—M 距离是否保留；M4看团簇形状是否从初始构型变成另一个异构体。若吸附前后团簇发生大幅重排，吸附能最好单独标注，因为它不再只是吸附物和固定团簇之间的相互作用。
如果计算资源允许，可以为 M4准备两个以上初始异构体。小团簇经常存在多个能量接近的构型，吸附物可能稳定其中某一个。只用一个初始团簇构型，有时会把“团簇异构体选择”误写成“吸附物偏好某个位点”。
# ======== M1/M2/M4 吸附优化参考示例 ========
SYSTEM = M4_support_OOH_relax
ENCUT = 520
EDIFF = 1E-5
EDIFFG = -0.03
ISMEAR = 0
SIGMA = 0.05
IBRION = 2
NSW = 150
ISIF = 2
# ======== 静态单点参考示例 ========
SYSTEM = M4_support_OOH_static
ENCUT = 520
EDIFF = 1E-6
IBRION = -1
NSW = 0
LORBIT = 11
LCHARG = .TRUE.
结果应该按什么单位归一化？
如果讨论原子利用率，可以把活性按金属原子数归一化；如果讨论单位位点活性，应明确最低能吸附构型对应的是一个金属中心、两个金属中心还是团簇边位。对催化路径来说，最可靠的写法是按“位点类型”分组，而不是把所有 Eads混成一个排行榜。
归一化还要看你想回答的问题。若关心贵金属利用率，M1的优势可能来自每个金属原子都暴露；若关心单个位点的反应能力，M4的多中心吸附可能更有利；若关心真实催化剂，团簇稳定性和生长倾向又会成为关键。同一组吸附能，按金属原子、按位点或按表面积归一化，得到的结论可能不同。
论文中可以把 M1、M2、M4分别列出“最低能构型”和“参与成键的金属数”。这样别人一眼能看出，M4的强吸附到底来自更强的单点电子结构，还是来自两个以上金属同时抓住中间体。
图5 金属团簇 ORR 活性示意，提示团簇位点的活性来自结构尺度、金属种类和吸附路径共同作用。DOI：10.1021/acsami.5c08441
如果 M1、M2、M4都用于回答“哪种结构更可能形成”，可以比较形成能或结合能；如果用于回答“哪类位点更适合吸附 *OOH”，可以比较各自最低能吸附构型和自由能，但要同时报告吸附方式。若要证明 M4比 M1更好，还需要说明这种团簇稳定存在、位点数量足够，并且不会在反应中继续长大。
可以直接比较的情况，是三者使用同一载体、同一超胞、同一吸附物参照，并且都完成了充分构型枚举。即便如此，结论也应该写成“在各自最低能构型下，M4对 *OOH 的稳定更强”，而不是写成“每个金属原子都更活泼”。后者需要额外的归一化和位点数量证据。
换句话说，M1、M2、M4可以放在同一项研究里，但不能被写成同一类位点的线性升级。它们更像三种不同结构家族：孤立中心、双位点协同和小团簇界面。
只能并列讨论的情况也很常见：比如 M1是 M—N₄，M4是金属团簇贴在石墨烯上，两者锚定位、局域电荷和吸附位点都不同。此时更适合说它们代表两类候选结构，而不是用一个吸附能数字宣布谁更优。
要点一：M₁、M₂、M₄ 的差异不只是金属数量，而是成键对象和吸附方式改变。
要点二：吸附能可以分别算，但不能不看构型就混排行。
要点三：结构稳定性和吸附强度要分开讨论，金属结合能不能替代吸附能。
要点四：不同尺度位点最好按位点家族分组，而不是直接写成谁一定更强。
