活性位点：反应含义、原子组成、数量与周转，以及工作条件如何改变它的身份

说明：本文华算科技主要介绍活性位点的反应含义、局部原子组成、数量与周转的关系，以及工作条件对位点身份的改变。

活性位点指的是怎样的局部结构？

活性位点是催化循环中直接容纳反应物、促使化学键变化，并能在产物离开后继续承担反应的局部原子集合。这个集合可能只有一个配位不饱和金属原子，或由相邻金属原子、缺陷周围的晶格原子、界面两侧的原子和吸附配体组成。“位点”描述的是一个可发生具体基元步骤的局部化学环境，元素名称只是其中一项信息。

催化循环还要求位点能够再生。反应物在局部结构上形成键并转化，产物释放后留下可继续接收反应物的位置。若某处只把中间体束缚在深能阱中，后续转化或脱附无法发生，该结构会积累覆盖物并压低周转。完成关键成键、断键和位点再生，才构成一个可持续的催化循环。

吸附位点、活性物种和活性相与它有关，含义却不相同。吸附位点只说明吸附物能够停留的局部能量极小构型，其中一部分只储存反应物或旁观物种；活性物种强调反应时的价态、配体或吸附态，例如 Ir–O；活性相指承担反应的物相，例如金属 Cu 或重构后的羟基氧化物。把三类对象合并成一句判断，容易把“能够吸附”误作“能够完成催化循环”。

图1中的氧化物衍生铜经历脱氧与分子动力学弛豫后，所得表面保留 Cu(111)，并生成方形排列、台阶和其他缺陷，呈现台面、晶界、空位与弯曲区域并存的结构。同一块铜表面因此包含多种局部原子集合；给整个颗粒贴上“Cu 活性位点”的标签，还没有指明反应实际发生在哪里。

活性位点还必须绑定某个基元步骤。CO₂ 电还原生成多碳产物时，两个相邻的 *CO 需要在足够近且取向合适的位置形成 C–C 键。图2把 155 个表面局部结构的筛选结果收敛到若干方形 Cu 原子集合，再用周期性 DFT 模型比较 *CO 偶联反应能和能垒。这里的反应中心由数个 Cu 原子共同构成，双分子靠近所需的间距和取向无法用其中某个 Cu 原子的元素标签表达。

为什么活性位点要写明反应步骤和产物？

一个局部结构降低某一步的能垒后，后续步骤仍可能转移到另一类局部结构。氧化物衍生铜上的平面方形、凸起方形和台阶方形位点都能促进 *CO 偶联；反应进入 *CH₂CHO 之后，C–O 断裂与继续加氢开始竞争，三类位点对中间体的配位方式随之分开。活性位点应写成“针对某个反应步骤或产物分支的位点”，否则不同论文里的“高活性”可能对应不同化学事件。

图3显示，平面方形和凸起方形位点更容易走向乙烯，台阶方形位点更利于醇类生成。台阶位的广义配位数较低，对 *CH₃CHO 结合更强；平面和凸起方形结构能够给 C–O 断裂后的 *O 留出四配位位置。“C–C 偶联位点”与“乙醇选择位点”对应不同的局部结构集合，后续分支仍由中间体键合与邻近空位共同决定。

当水解离、含碳中间体稳定或酸碱异裂分别由不同原子承担时，反应需要相邻位点协同。一个位置负责吸附并解离 H₂O，另一个位置稳定含碳中间体；酸位与碱位共同完成异裂；金属原子与邻近晶格氧共同形成 O–O 键。此时应报告原子间距、覆盖度和相邻空位概率。把每个吸附原子分别计作独立位点，会漏掉双位点或多位点反应对几何邻接的要求。

位点数量怎样与催化速率相连？

宏观速率由多类位点共同累积。对第 j 类位点，可以用位点数 N_j、处于可反应状态的比例 θ_j 和单位位点周转频率 TOF_j 表示其贡献；总速率近似写成 r = ΣN_jθ_jTOF_j。高总活性可能来自位点更多，也可能来自单个位点转得更快。质量活性、几何面积电流或总产率没有自动消除这两种来源。

计算给出的单位位点速率与实验电流比较时，分母尤其关键。若把全部金属含量当作位点数，埋在体相或被覆盖的原子会拉低表观 TOF；若只选最活跃的一类稀有缺陷作分母，又可能高估整个催化剂的平均周转。位点密度、可接近比例和单位位点速率需要分别报告，才能判断性能变化发生在哪个乘数上。

图4利用退火逐步抹平台阶和缺陷。三类方形位点的比例下降时，C₂₊ 产物电流同步下降；平面方形与凸起方形位点的占比对应乙烯电流，台阶方形位点的占比对应醇类电流。结构数量与产物速率按位点类型分别变化，使原子模型与宏观电流连接起来。不过线性趋势仍依赖位点统计方法、表面面积归一化以及退火未同时改变其他速率变量等条件。

统计到的位点中，只有能够接触反应物且未被竞争物种占据的部分可以进入周转。孔道堵塞、强吸附物覆盖、双电层中的离子占据和传质限制都会改变 θ_j；颗粒内部被包埋的原子也不参与表面反应。实验上用化学吸附、电化学探针或滴定估算位点密度时，探针的可达范围与目标中间体的几何要求需要同时检查。

图5中的 Ir 原子分散在重构后的 NiFeO 表面，显微图、Ir–O 配位和结构拟合共同说明这些原子的存在方式。结构表征确定“有多少、长什么样”，也限定了 DFT 模型中的近邻原子；单位位点速率由反应态覆盖和能垒决定。

工作条件为什么会改变位点身份？

真空优化得到的清洁表面只是一个参考状态。温度、分压、电位、pH、溶剂和覆盖度能够引起氧化还原、吸附诱导重构、配体交换与表面溶出。工作态活性位点应包含当时的价态、近邻配位和吸附物。空气中测得的前驱体结构若在反应开始后转成另一种表面，前驱体原子坐标不能直接代表催化循环中的反应中心。

图6记录 Ir–NiFeO 在开路、1.45 V 和 1.55 V（相对于 RHE）下的 XANES 与 EXAFS。开路时 OH⁻ 和 H₂O 吸附增加 Ir–O 配位；进入析氧电位后，Ir 的氧化态与局部配位继续变化；更高电位下发生 Ir–OH 脱质子，形成能参与水进攻和 O–O 偶联的 Ir–O 组合。位点身份由“孤立 Ir 原子”细化成了指定电位下的 Ir–O 协同结构，名称中多出的 O 正是反应态变化的一部分。

为了让计算对象匹配工作态，DFT 需要对这些可能存在的化学状态分别建模。图7比较 Ni、Fe、Ir 和 Ir–O 位点的 OER 自由能，Ir–O 模型中形成 *OOH 的自由能代价最低；PDOS 用于检查引入 Ir 后近费米能级电子态的变化。自由能路径检验局部结构能否完成连续基元步骤，PDOS 只解释其中的电子态变化；位点归属需要检查关键过渡态、覆盖度、竞争反应以及该结构在目标电位和化学势下能否维持。

在温度较高或表面容易重构的体系中，实验速率还可能来自一组持续互变的结构。此时计算应比较多个局部极小值的自由能、相互转化时间尺度和各自反应速率，不能只保留 0 K 下总能最低的构型。占比不高但周转很快的结构也可能贡献可观通量。

材料计算中的位点结论应写到可检验的粒度：具体晶面或界面、原子编号与近邻配位、承担的基元步骤、吸附覆盖、温度或电位，以及支持判断的能垒或自由能。若模型没有包含重构、显式溶剂、恒电位效应或竞争覆盖，结论只适用于所用静态结构和参考条件，不能外推为材料在所有工况下唯一的活性中心。