单原子锚定位点怎么选？从结构构建到稳定性评估再到反应匹配

说明：本文华算科技主要介绍单原子锚定位点的结构来源、稳定性指标和反应适用条件。

单原子锚定位点的计算对象是什么

单原子催化剂里的“锚定位点”，不是普通吸附物在表面挑一个顶位、桥位或空位。它指的是载体上能够固定孤立金属原子的局域结构：邻近原子给出配位，缺陷或孔洞给出几何约束，载体能带和费米能级附近态密度参与成键。DFT 模型里，常见写法包括 M-N_x、M-C_x、M-O_x、金属@空位和金属@边缘位，M 是单个金属原子，x 对应最近邻配位原子数。

锚定位点的选择从三个量开始：单原子与载体的结合能 E_bind，孤立金属或金属体相的参考能，载体本身形成该缺陷的代价。E_bind越负，单个金属原子和该局域结构之间的相互作用越强；但若载体缺陷形成能很高，实验中该位点浓度可能很低。另一个约束来自金属的内聚能 E_coh，当金属-载体结合强度弱于金属-金属聚集倾向时，模型里的单原子在热处理或电化学条件下可能转向团簇。

用气相金属原子作参考时，E_bind偏向描述“单个原子放到载体上”的能量收益；用金属体相或金属团簇参考态时，结果更接近抗团聚能力。超胞尺寸也会改变数值：同一个 M-N₄位点在小超胞里等效覆盖度较高，金属中心之间可能产生镜像相互作用。参考态、超胞尺寸和覆盖度要和 E_bind一起记录，否则不同论文里的锚定能很难横向比较。

单空位和双空位给出两类常见锚定位点，单空位周围常形成三配位口袋，金属原子可能明显高出石墨烯平面；双空位能让四个邻近原子参与配位，部分金属会进入更接近平面的结构。面外位移、M-X 键长和配位数会影响后续中间体的吸附几何，若两种位点只按空位数量命名，后续吸附几何差异会变成一个标签。

图2. 单空位和双空位石墨烯中，过渡金属原子形成三配位或四配位锚定结构。DOI：10.1039/d6cp00932h

哪些局域结构更容易固定单原子

碳载体里常见的锚定位点来自吡啶 N、石墨 N 邻近空位、双空位、边缘缺口和含 B、P、S、O 的杂原子配位。N 原子电负性较高，容易降低邻近碳位的电子密度；B 或 P 会改变孔洞附近的键长和局域应变。对于同一个 Fe、Co 或 Ni 原子，TM-3N₃、TM-4N₄、TM-2N2C不是名字差异，它们给出的配位场、d 轨道分裂和自旋态可能完全不同。

同一类空位还存在掺杂原子排列差异。双空位里两个 N 原子可以呈 trans、cis(5-ring) 或 cis(6-ring) 排布，邻近配位原子的相对位置会改变金属中心的对称性。配位原子种类和孔洞大小共同决定金属是否进入平面、是否形成四个短 M-X 键，以及费米能级附近是否出现局域 d 态。

氧化物载体里的锚定位点常来自氧空位和低配位金属阳离子、台阶边缘和表面羟基。氧空位可让 Au、Pt、Rh、Ru 等金属原子靠近还原态阳离子，形成 M-O-Ce、M-O-Ti 或 M-O-In 等局域结构。此时计算对象包含氧化学势、自旋态、U 值和表面电荷补偿方式；这些设置会改变空位形成能和金属原子电荷分布。

氮化物、硫化物和硼氮化物载体还会给出不同的锚定图像。Ni-TiN 例子里，Ni 原子可位于氮空位、完整 N-top 位点，也可形成小团簇；三种模型的局域配位和金属-载体接触面积不同，后续 CH₄ 裂解路径也不同。用于选位点时，缺陷位、完整位和团簇模型应并列计算，否则可能把建模最方便的位点当成真实反应中心。

锚定强弱怎样进入稳定性筛选

稳定性筛选可比较同一载体上不同金属、不同配位口袋的 E_bind、M-X 键长和金属面外位移。以 TM-4NDG 为例，Ti、V、Cr、Mn、Fe、Co、Ni、Cu、Zn 都能放入四氮掺杂石墨烯空位，但键能、TM-N 距离和 Δz 并不相同。Fe、Co、Ni 接近平面时，四配位口袋对金属原子的几何约束更强；Ti、V 的面外位移则提示金属半径和轨道占据会改变局域结构。

结合能给出总量，能量分解能看出差异来源。TM-4NDG 的相互作用可分为静电吸引、Pauli 排斥和轨道相互作用。若某个金属的 ΔE_oi更负，通常对应 d 轨道与 N 掺杂石墨烯片段轨道之间的杂化增强；若 Pauli 排斥上升，短键长未必带来更强有效稳定。E_bind、键长、Δz、PDOS 或 COHP同步给出趋势时，锚定位点筛选才有物理约束。

静态结构仍需接受迁移检查。金属原子从强锚定位点迁移到邻近弱位点时，NEB 能垒 E_diff决定热处理或反应温度下的可动性。Ni-TiN 里，Ni 从氮空位向 N-top 位点迁移时经过两个过渡态，能垒分别出现在局域配位改变的位置。若 E_diff偏低，即使初始 E_bind较负，金属原子仍可能在有限温度下寻找新的配位环境或参与团聚。

高温合成、还原气氛和电化学电位会改变单原子周围的配体。AIMD 可以检查几皮秒到几十皮秒内的配位保持情况，团簇竞争模型可比较 M₁、M₂、M₄在同一缺陷附近的能量。若 M₂或 M₄结构能量明显更低，单原子模型需要配合低负载、空间隔离或配体饱和条件，否则孤立金属中心在真实样品中会转向多核结构。

实验谱学能给稳定性筛选增加约束。HAADF-STEM 的孤立亮点只能确认重元素分散状态，EXAFS 中 M-N 或 M-O 第一配位峰、XANES 吸收边位置、CO-DRIFTS 中线性吸附峰才更接近局域配位。若计算给出M₂、M₄竞争结构能量较低，而谱学仍呈现单配位壳层，样品制备中的空间隔离或配体限制就要进入模型假设。

锚定位点还要怎样服务反应判断

单原子固定以后，锚定位点仍要与目标反应匹配。HER 中可比较 ΔG_H*，CO₂RR 中可比较 *COOH、*CO、*OCHO，ORR 中可比较 *OOH、*O、*OH。若金属-载体相互作用过弱，单原子稳定性不足；若相互作用过强，金属 d 态可能远离费米能级，中间体吸附和产物脱附都会改变。筛选锚定位点时，稳定性和活性不是两个独立表格，同一局域配位会同步改写 E_bind、d 带中心、Bader 电荷和中间体自由能。

NGA-COF@Pt 的计算给出一个清楚例子：不同 Pt 配位环境下，HER 自由能曲线、TDOS 和 d 带中心发生同步变化。Pt-N₄、Pt-2H、Pt-2Cl、Pt-NP 等模型中，H* 吸附强度和电子态位置不同，最适合固定 Pt 的结构未必就是 ΔG_H*最接近 0 eV 的结构。对于 CO₂RR 或 OER 体系，也可把锚定模型推进到关键中间体，比较自由能图、差分电荷密度和反应态几何。

电催化模型还会受到电位、pH、显式水分子和表面覆盖度影响。M-N₄在真空 slab 中保持平面配位，并不保证在含 *OH 或 *CO 覆盖时仍维持相同价态；氧化物上的 M-O_x位点在还原电位下可能产生新的氧空位。恒电位计算、溶剂修正和覆盖度扫描可以把锚定稳定性与工作条件靠近，0 K 优化结构只能作为初始构型证据。

带电中间体会引入界面偶极，周期性 slab 里的真空层、电荷补偿和偶极修正都会影响吸附自由能。对于 CO₂RR，*COOH 与 *OCHO 的偶极方向不同，锚定位点周围的局域电场会改变二者相对稳定性；对于 ORR，*OOH 的端位吸附和侧向吸附会给出不同 O-O 键长。吸附构型、偶极修正和显式水氢键应和锚定模型一并保留。

实际建模中，一个可用的锚定位点至少要同时满足四类信号：载体形成该缺陷的热力学代价可接受，M-X 键长和配位数在结构优化后保持稳定，E_bind相对 E_coh能抵消金属团聚倾向，E_diff不低到让单原子在工作温度下快速迁移。之后才进入反应自由能、NEB 反应能垒、溶剂或电位修正。锚定位点筛选停在一个负的结合能数值时，只回答“单原子能否暂时放在这里”；迁移能垒、轨道杂化和中间体自由能加入后，模型才开始回答“这个局域结构能否作为反应位点”。