催化表面电子效应：电荷转移、d带中心、EMSI、轨道势模型全解析

说明：本文华算科技主要介绍催化表面电子效应的计算含义，以及电荷转移、d 带中心、电子金属—载体相互作用（EMSI）和轨道势模型各自描述的物理过程。

金属位点与吸附分子形成化学键时，表面原子的价电子会重新分布。合金中的邻近原子、氧化物载体、局部应变和吸附物都能改变位点附近的电荷密度、能级位置与轨道占据；由此产生的吸附能或反应能差异，通常称为电子效应。它描述电子态如何响应结构与环境变化；数值比较须指定表面位点与吸附构型。

以氧化物负载金属为例，孤立金属颗粒和贴附在载体上的颗粒拥有不同的电子结构。接触面会形成金属—氧键，界面附近的电荷可能重新分配；颗粒外表面的反应位点又与埋在接触面内的金属原子不同。位点位置和颗粒尺寸均发生改变时，吸附能差异可能含有配位数和电子态两种贡献。

静态 DFT 通常从同一几何与电荷参照下的差分电荷密度、局域态密度和投影态密度提取电子响应。吸附能则反映给定构型的总能变化；两者属于不同观测量。对电催化表面，电位、溶剂和覆盖度还会改变界面电荷与吸附物稳定程度，真空中的中性 slab 结果只对应其自身模型。

Cu 掺杂后的晶格应变会改变 Cu—Cu 距离，也会改变 Cu 的 d 态，几何与电子效应常有共同来源。保持比较对象的晶面、吸附位和覆盖度一致，再观察电子态与吸附能的联动，才能将结构变化带来的贡献限制在可讨论的范围内。

金属与载体接触后，费米能级平衡和界面成键可使接触面附近的电子密度重新分布。计算中常用差分电荷密度 Δρ=ρ_整体−ρ_金属−ρ_载体 表示相同原子位置与计算网格下的密度差。正负等值面标出空间积累与耗尽；Bader 分区再给出指定原子盆内的电子数变化。差分图中的局部极化并不要求整个颗粒转移整数个电子。

Pt/CeO₂ 是一个可量化的例子。颗粒增大时，单颗粒转移到 CeO₂ 的总电子数上升，每个 Pt 原子的平均转移量却在约 30—70 个原子的颗粒附近达到峰值，随后下降。接触面积、颗粒表面原子比例与金属内部屏蔽共同影响这个趋势；“颗粒越大、每个原子转移越多”不符合该组模型。

在同一 Pt₈/CeO₂ 模型里，改变 Pt 向载体转移的电子数后，CO 与 H₂O 的吸附能明显变化，O 的数值变化较小。吸附物轨道与金属态的耦合方式决定电荷变化如何传到吸附键，因而电荷转移量相近的两个体系也可能具有不同的反应结果。

电荷转移有静态电荷分区的含义，也有实际反应中的电子传输含义。Bader 数值和差分密度描述给定构型的电子分布，电子传输速率涉及态间耦合、势垒、电位及时间尺度。把静态计算所得“转移了多少电子”直接写为“电子转移更快”，会越过这些动力学条件。

过渡金属表面的 d 带中心 ε_d 是选定金属 d 态 PDOS 的能量加权平均：ε_d=∫Eρ_d(E)dE/∫ρ_d(E)dE。比较不同位点时，投影原子、能量参照、积分窗口和自旋分支须保持一致。数值上移表示 d 态平均能量更接近所选参照；它本身没有写入 d 带宽度、轨道对称性或具体吸附分子的能级。

金属 d 态与吸附物轨道耦合后形成成键态和反键态。能级间距、耦合强度及反键态占据共同控制吸附键。对于一些同族过渡金属表面，d 带中心随金属种类移动能追踪吸附强弱；合金或单原子位点的轨道形状改变时，同一中心值可以对应不同的吸附能。

硝酸根还原为 NH₃ 的反应同时涉及桥式双齿吸附的 *NO₃ 和空位吸附的 *N。Cu(100) 与 Cu(111) 上各中间体的构型及自由能不同，原研究以*NO₃ 和 *N 的吸附能绘制活性火山图。若两种吸附物都随 ε_d 同向增强，单靠移动 d 带中心难以单独调节其中一个中间体。

Cu(100) 的吸附物轨道投影态密度显示，*N 的 p 轨道有成键和反键特征，*NO₃ 的前线轨道则呈现另一种能级分布。表面金属 d 态与不同吸附物轨道的重叠决定能级移动如何变成吸附能变化；相同的 d 带中心位移，不保证两种中间体获得相同的能量增量。

图中的模型扰动把 d 带中心上下移动时，*NO₃ 与 *N 吸附能大体同步改变。B2 有序金属间化合物的 (100) 位点还存在另一项：较短层间距使近表层金属 d 轨道与空位 *N 的 p 轨道发生更强的泡利排斥，*NO₃ 的桥式构型受其影响较小。CuPd 的预测位置由两种吸附能共同确定，说明位点几何和轨道排斥可使吸附能偏离简单的线性联动。

d 带中心适合比较结构接近、成键轨道相似的金属表面。氧化物、强共价位点或显著自旋极化的单原子中心，还需核对参与吸附的具体轨道与反键态占据。把电催化活性与真空 PDOS 直接对应时，还要记录电位、溶剂和覆盖度；这些条件会改变自由能与界面电荷。

电子金属—载体相互作用（EMSI）指载体通过界面成键、电荷重排和轨道耦合改变金属位点的电子态。金属—载体作用还可表现为颗粒形貌改变、金属包覆或氧迁移；这些结构与化学变化的证据应分别记录。对纳米颗粒，金属内部的电荷屏蔽使接触面附近的电子扰动最显著，暴露于反应物的颗粒外表面未必具有同等电荷变化。

Pd₁₁₉ 与 MgO(100) 接触后，界面第一层出现明显电子密度极化，远离载体的 Pd 原子变化很小；H、O 和 CO 在颗粒周边位点的吸附响应又各不相同。周边位点的吸附能与内侧位点分开计算，才能解释为何局部界面作用会改变 CO 氧化的失活起点。

在 Pt/B₄C 与 Pt/C 催化剂中，无外加电位的 XPS 测得 Pt 4f 与 C 1s 相对结合能差异；在相同外加电位下，原位 Pt L 边 XANES 却显示 Pt/B₄C 有较多未占据 d 态。真空中的电荷平衡和电化学双电层中的电荷状态采用不同参照，不能拿两个实验的“富电子/缺电子”标签直接互换。

当金属颗粒的粒径、晶面或覆盖度也随载体改变时，催化电流的差异含有多种来源。Pt/B₄C 的比较控制了部分颗粒尺寸与分散条件，并记录了谱学、电化学活性和稳定性；将载体作用指定为电子贡献，还需让金属位点的结构与反应环境具有可比性。

“轨道势”对应原研究提出的 orbital potential（G），研究对象是 H₂O₂ 在 TiO₂ 表面的化学吸附。模型把表面原子的价轨道参与体相能带的份额记作 f_B，参与局域表面配位键的份额记作 f_S，并以 f_B+f_S=1 描述两种分配。G 用来表征轨道由体相成键转向表面成键的竞争，不是电极电势，也不是 DFT 输出文件里直接列出的单一能量。

在该模型的粒径关系中，颗粒缩小会改变体相成键份额与表面成键份额，因而改变表面反应性。图中的蓝线、红线分别代表模型内的轨道势和表面反应性随粒径的趋势；它们是特定假设下的理论函数，并非所有氧化物或所有吸附物共享的实验曲线。纳米颗粒暴露晶面、缺陷数量和表面配位数若均随粒径变化，测得的趋势会含有活性位数量和局部配位变化的贡献。

把该模型用于另一种氧化物表面时，至少需检验尺寸系列的晶相、暴露面、缺陷浓度与 H₂O₂ 覆盖度。对于 TiO₂，同一种粒径变化若还引起氧空位增多，吸附量上升便可能含有缺陷位点增加的贡献；固定表面结构后的计算或可比实验，才有机会检验 f_S 与粒径的关系。