说明:本文华算科技主要介绍异质界面在催化计算中怎样改变活性位、电子结构、吸附中间体和反应能量。


异质界面的计算对象是什么?
异质界面在催化计算中指两种组成、晶相或局域配位不同的区域共用一段接触面,常见对象包括金属/氧化物、硫化物/硫化物、氧化物/碳载体和金属团簇/配位碳。DFT模型通常把接触区域写入同一个周期 slab,真空层、晶格匹配、表面终止、接触面原子排布和吸附物覆盖度共同限定 界面位点 的化学环境。催化活性讨论的对象由单个裸露原子扩展为吸附物、邻近配位原子、载体和局域电荷共同形成的 反应微区。

FeMoO4@Ds-AC 这类氧化物/碳复合材料提供了一个直观样本:显微形貌、晶格条纹和元素分布确认氧化物颗粒与碳基底存在紧密接触,DFT 中再把 FeMoO4 片段与碳片写成可优化的界面结构。真实颗粒接触 与 周期界面模型 属于两类对象,前者来自 TEM、EDS、XPS 等表征,后者承担电子态和反应能量计算;两者一致的部分通常是接触元素、近邻配位和可吸附区域。
模型建立后,界面讨论才进入可计算变量:Fe-C 或 Mo-C 距离、接触区键级、费米能级附近态密度、平面平均电势、差分电荷密度和吸附构型能量。界面本征活性 不由某个材料名称决定,而由接触面上原子配位、轨道能级和吸附物几何共同决定。若只把 FeMoO4 和碳分别计算,再把两个结果相加,界面处的 轨道杂化、局域极化和电荷重排会从模型中消失。

异质界面的活性改变常从结构项开始,因为两相晶格常数不同,接触面会出现面内应变、错配位和局部键长偏移。slab 中固定哪几层、放松哪几层,会改变 Fe-C、Mo-C、Pd-O、Ru-O 等界面键长;界面形成能 可用于比较不同堆垛或终止面,吸附前后键长变化 则把结构扰动连接到中间体稳定程度。
周期模型还会把界面沿面内方向无限复制,低覆盖度吸附物在相邻超胞间可能产生人为相互作用。扩大表面超胞、比较不同吸附位和检查偶极修正,可以减少电荷重排在真空方向的串扰;slab 厚度、真空层 与 k 点采样 改变后,吸附能和费米能级附近态密度若保持同一趋势,界面结论才不容易被模型尺寸主导。


界面为什么会改变吸附中间体?
催化提升在计算结果里常表现为 *H、*OH、*O、*OOH、CO* 等中间体的吸附自由能进入较合适的区间。异质界面改变吸附态的直接原因通常有三类:接触面电荷重新分布改变吸附原子的局域电势,金属 d 轨道与吸附物 p 轨道发生 轨道杂化,低配位或混合配位原子形成不同于单相表面的 配位不饱和位点。
Pd-CeO2 界面的 HOR 计算把这种关系写得很清楚:Pd 团簇向 CeO2 侧转移电子后,Pd 3d 态的 d 带中心下移,*OH 与 HBE 等吸附指标同步变化。Bader 电荷 在这里只给出净电荷分配,差分电荷密度 给出空间重排区域,PDOS 才能判断吸附物轨道与金属态的能量重叠;三个量对应同一界面,却回答不同问题。

RuO2/CoMoO4 异质结把 OER 和 HER 相关中间体放在同一组计算中比较,平面平均差分电荷沿 z 方向给出电子积累和耗尽区域,OER 的 *OH、*O、*OOH 自由能与 *H 吸附自由能对应不同反应分支。界面产生的 局域电场 和 电子富集区 会改变吸附物成键强度,但吸附能过强会拖慢中间体转化,吸附能过弱又会削弱初始活化;有利区间由具体反应坐标给出。

差分电荷与吸附自由能对应到同一组吸附构型后,界面效应才从电子重排进入催化判断。Pd-CeO2 中 Pd 侧电子耗尽会削弱过强的 *OH 键合,RuO2/CoMoO4 中不同金属-氧配位会改变 *O 与 *OOH 的相对能量;吸附位点编号、覆盖度 和 参考态 必须随自由能数值一同保留。
PDOS 重叠只反映态密度在能量轴上的相遇,成键强弱还可由 COHP、ICOHP 或 ELF 辅助判断。若界面让金属 d 态与吸附物反键轨道的占据降低,吸附键可能变弱;若吸附物 p 态与界面金属态在费米能级以下形成更明显的成键贡献,吸附态稳定程度会提高。成键态/反键态占据 比单独的峰高更接近化学键强弱的计算含义。


反应自由能和能垒怎样随界面改变?
异质界面提升催化并不只体现在某个吸附构型能量降低,反应坐标上的每一步都可能重新排序。对 HER 来说,Volmer 步的 H2O 解离、*H 形成、Heyrovsky 或 Tafel 成氢各自有不同能量项;对 OER 来说,AEM 路线中的 *OH、*O、*OOH 转化与 LOM 路线中的晶格氧参与给出不同能量轮廓。自由能图 的作用在于定位高能步骤,而 NEB 能垒 或水解离能用于补足动力学信息。
NiCoSx@Ni3S2 的计算同时列出功函数差、DOS、d 带中心、HER 的 Volmer-Tafel/Volmer-Heyrovsky 曲线以及 OER 的 AEM/LOM 路线。界面处 Ni 位和 Co 位的 *H、*OH、*O、*OOH 稳定程度不同,反应优势来自 多位点分工 和 步骤能量匹配。若只比较单相 Ni3S2 的 ΔGH*,水解离和 OER 中间体转化的能量变化会被遗漏。

Ru/RuO2 金属/氧化物界面提供了另一个 OER 样本。差分电荷和 Bader 结果显示金属 Ru 与 RuO2 接触后电荷分布发生偏移,Ru 4d 与 O 2p 的 PDOS 反映 t2g/eg 相关态变化,OER 自由能曲线把 *OH 到 *OOH 的能量差列入同一坐标。这里的 d 轨道调制 作为孤立描述符时信息不足,只有与 OER 中间体自由能 同时出现时,才构成对活性差异的计算支撑。

自由能曲线还区分热力学高能步骤和动力学能垒。若 *OH 形成自由能降低,但水解离能或 Heyrovsky 过渡态仍然偏高,HER 速率仍可能受水活化限制;若 *OOH 形成能降低,但产物 O2 脱附或晶格氧恢复代价上升,OER 的表观电流也可能受后续步骤限制。界面协同 在计算文本中应写成“哪一个位点稳定哪一个中间体、哪一步能量下降、哪一步仍是高能步骤”,而不是只给出总过电位差。


异质界面结论怎样限定计算条件?
界面计算给出的结论依赖表面终止、覆盖度、电位参照、溶剂处理、自旋态和反应后结构。用于解释 OER 的 Ru/RuO2 slab 若只代表反应前结构,所得 PDOS 和自由能适合描述初始界面;在高电位和酸性介质下,Ru 价态、Ru-O 配位和晶格氧参与程度可能改变,fresh/spent 谱学 与同位素 DEMS 信号就成为限定计算结论范围的实验约束。

把异质界面写成催化提升原因时,正文应给出材料体系、接触面原子、反应类型和计算参照态,例如 Pd-CeO2 用于 HOR,RuO2/CoMoO4 和 Ru/RuO2 用于 OER/HER,FeMoO4@Ds-AC 用于氧化物/碳接触。同一差分电荷图 只能描述电子密度空间重排,同一 d 带中心 只能描述金属 d 态位置,同一自由能图 只覆盖给定反应路线;把它们与吸附构型、位点配位和反应后相态配套校验,异质界面提升催化的判断才有明确适用范围。
电位和溶剂处理会改写同一界面的吸附序列,因为带电 slab、隐式溶剂、显式水层和 CHE 参照给出的自由能项并不相同。碱性 HER 中 H2O 解离和 OH– 脱附占有较大权重,酸性 OER 中 Ru-O 键稳定和晶格氧交换更敏感;pH 条件、电位标尺、溶剂模型 与 自旋态 改变时,界面活性排序可能随之改变。
界面形成稳定性和反应中保留状态也要分开处理。形成能较低只表示给定原子排布在热力学上较有利,AIMD 或反应后谱学才对应工作条件下的结构保留;若催化过程中生成羟基化层、氧空位或金属溶出,初始异质界面 与 真实工作态界面 的活性位可能已经不同。
当界面位点使过强吸附变弱、过弱吸附变强,或把水解离、质子/电子转移和产物脱附分配到不同近邻位点,反应自由能曲线会出现较低的高能步骤;当反应后相态改变,原本的静态界面模型则只保留为初始结构参照。异质界面的价值 来自可计算、可表征、可限定的局域化学环境,区别于材料名称相加后的经验判断。
