说明:本文华算科技主要介绍活性位点在材料计算中的定位对象、候选位点筛选、电子结构证据和工作状态校验。


活性位点是什么?
活性位点在 DFT文章里不是一个孤立原子的名字,而是 反应物、表面结构和局部电子态共同限定的反应中心。同一个金属原子在平整晶面、边缘、缺陷邻位、异质界面或孔道内壁中,配位数、近邻元素和局域电势都不同;吸附物靠近之后,真正参与成键的原子、相邻原子的极化、表面是否重构,都会改变活性判断。写成“Fe 位点”“Cu 位点”之前,模型里应已经有 晶面指数、终止方式、覆盖度、吸附物构型和自旋状态。

定位活性位点的第一步是给出可复查的结构对象。周期性 slab 模型里,候选位点常来自 顶位、桥位、空位邻位、台阶位、边缘位和界面位;分子筛或 MOF 模型里,候选对象还包括金属配位数、孔径限制和可弛豫区域。图1的蓝色矩形标出不同材料上氢的最低能吸附位置,读者能直接看到:吸附位置随表面原子排列变化,活性位点并不自动等同于元素种类。
很多误判来自模型对象没有固定。若只写“NiFe-LDH 的 Fe 是活性位”,Fe 处在本征层板、掺杂邻位还是电化学重构后的羟基环境中,判断含义并不相同。活性位点定位要把原子编号、近邻配位、表面终止和吸附物姿态写在同一个模型语境里,否则后面的吸附能、Bader 电荷和 PDOS 都会变成互相难以比较的数字。
在孔道催化体系里,金属中心、骨架酸性位和空间位阻会把活性位点限定成局部反应环境。图2中的 Au 交换 ZSM-5 和 MCM-22 簇模型把金属中心、骨架 T 位和局部可弛豫区域固定下来;CO2 或环氧化物进入孔道后,反应中心不只是 Au+,还包含邻近骨架 O、酸性位和空间位阻。活性位点是一个局域反应环境,原子坐标、配位饱和度和孔道几何共同决定反应物能否采用可反应的吸附姿态。



候选位点怎样缩小范围?
候选位点筛选通常从反应物或关键中间体吸附开始,但吸附能只回答一个局部问题:在给定模型、参考态和构型下,吸附物与表面相互作用有多强。用于定位活性位点时,吸附能要和键长、吸附姿态、覆盖度和后续反应步相连。过弱的吸附常见于反应物难以活化,过强的吸附会让中间体或产物停留时间过长;两侧都会拉低周转频率。
图3给出 Au-[M]-ZSM-5 和 Au-[M]-MCM-22 中环氧乙烷吸附态、过渡态和产物吸附态的优化构型。用于判断位点时,比较 C-O 键长、Au-O 接触距离、环张力释放和骨架 T 位邻近环境。同一反应物在不同孔道与酸性位环境中的几何响应,比单个吸附能数值更能提示哪类位点参与键活化。

吸附构型还要接受反应路径检验。某个位点也许能稳定 CO2 或环氧化物,却可能在成环过渡态、质子转移或产物脱附步骤产生较高能垒。图4把不同分子筛骨架和环氧化物的活化自由能放入同一相关图,位点差异转化为过渡态自由能的变化。

在电催化体系中,电位、pH、溶剂和离子环境会改变吸附自由能。CHE 近似、显式水层、隐式溶剂、表面带电模型给出的趋势可能不同,尤其是 *H、*OH、*OOH、*CO、*OCHO 这类带有明显极化的中间体。候选位点从“能吸附”进入“能完成关键步”时,才接近活性判断;模型没有覆盖反应路径时,定位结果只能称为吸附偏好。
同一表面上还可能存在覆盖度效应。低覆盖度下稳定的位点,在高覆盖度或共吸附条件下会因侧向排斥、氢键网络、表面偶极或位阻改变排序。对于 HER、OER、ORR、CO2RR 等反应,覆盖度和参考态必须跟自由能图保持一致,否则不同位点之间的能量差会混入模型尺寸误差和电荷补偿方式差异。


电子结构证据怎样限定真正参与反应的原子?
几何和能量能缩小候选范围,电子结构负责回答另一个问题:哪一个局域原子轨道在反应步中改变了成键关系。Bader 电荷给出原子盆电子数变化,差分电荷密度给出空间积累和耗尽区域,PDOS 反映费米能级附近轨道贡献,COHP 或 ICOHP 描述成键/反键占据。几类图谱回答的问题不同,数值方向一致时,活性位点判断才具有较强约束。
NiFe-LDH 是典型例子。图5把 Fe 周围原子结构、Bader 电荷、差分电荷和 COHP 纳入同一组计算:Ga3+ 掺入后,Fe 邻近的电子密度分布和 Fe-O 成键强度发生改变。这里的重点不是“电荷多一点就活性高”,而是 Fe 位点附近的电子抽取、轨道成键和反应中间体稳定化方向相互吻合。

PDOS 的用法也要具体到轨道和能区。费米能级附近 Fe 3d 态增强,可能提高与含氧中间体轨道相互作用的机会;若反键态被过度占据,吸附物反而可能难以转化或脱附。图6中的 Fe 位点 pDOS、AEM 路径和 Gibbs 自由能可用于判断电子态变化是否真的转化为 *OH、*O、*OOH 等中间体自由能的改变。

差分电荷密度常被误写成价态结论。正确处理方式是固定颜色含义、等值面大小和参考态:吸附后体系电荷密度减去干净表面与孤立吸附物电荷密度,得到的是空间重排,不是形式价态。形式价态仍要由 XPS、XANES、磁矩、配位环境和投影态共同限定;DFT 中的电荷分区还会受赝势、投影半径和算法影响。
COHP、ICOHP 或类似成键分析适合处理“哪个键在反应中被削弱或增强”。例如含氧中间体吸附在 Fe 位点时,Fe-O 反键态占据降低、ICOHP 绝对值增加,通常说明键相互作用增强;若后续脱附或 O-O 成键受阻,过强 Fe-O 相互作用又会形成动力学代价。活性位点的电子结构证据要与反应坐标中的键变化对应,而不是停留在单张静态电荷图。


工作状态会重新定义活性位点
很多催化材料在反应电位、气氛或溶液中会发生表面重构。预处理后的晶面、空气中表征到的价态、真空 DFT 模型里的终止方式,未必就是工作电位下的反应中心。对于 OER 中的羟基化层、CO2RR 中的氧化物衍生金属、HER 中的氢覆盖表面,活性位点判断应包含工作状态下的结构和吸附物。
δ-Bi2O3/La2O3 的 CO2 电还原案例说明了工作状态的影响。图7包含 原位 Raman、原位 FT-IR 和准原位 XPS 信号,能够跟踪电位变化下 Bi-O、吸附中间体和表面化学态的响应。若计算只采用反应前氧化物模型,谱学中出现的还原、缺陷生成或吸附物峰位变化会让位点归属发生偏移。

计算模型可以把工作状态转化为可比较的结构假设。图8给出 δ-Bi2O3 与 δ-Bi2O3/La2O3 异质结构的 DFT 优化模型、pCOHP 和 CO2RR/HER 路径。pCOHP 用于判断 Bi-O、Bi-La 等局域键相互作用,反应路径用于比较 *HCOO、*HCOOH 和 *H 的相对稳定范围。

把活性位点写成“某元素”时,电位、界面结构、中间体和竞争反应会从结论里消失。更精确的写法是“在 -1.15 V vs. RHE 条件下,La2O3 邻近的 δ-Bi2O3 界面 Bi 位点更有利于 *HCOO 路径”,其中包含电位、界面结构、中间体和竞争 HER。位点、反应物和工作电位同写在一句话里,结论的适用范围才不会被扩大到所有 Bi 位点。
形成能、表面能、相图、AIMD、声子谱和 NEB 能垒分别回答缺陷热力学、表面构型保持和过渡态高度。某缺陷是否热力学可形成,某表面是否在温度下保持构型,某反应步是否存在过高过渡态,对活性位点可信度有不同约束。活性位点定位完成后,仍要检查 位点在反应环境中是否保留、关键中间体是否可转化、产物是否能离开表面。文章写到这里就停在材料状态和图谱信号上,避免把局部计算趋势扩大成普遍催化结论。
