说明:本文华算科技主要介绍催化剂中毒在材料计算中的含义、计算来源、常见图谱和适用范围。




催化剂中毒指毒物分子、原子或离子在活性位附近形成持续占据,使原本参与反应物吸附、键断裂、质子/电子转移或产物脱附的位置被改写。材料计算里讨论的中毒,通常从毒物吸附构型、吸附能和覆盖度开始;这些量把“位点被占”转成可比较的结构和能量结果。
毒物的来源可以是反应物中的杂质,例如 CO、H2S、SO2、Cl– 或含磷分子,也可以是反应过程中生成的强吸附中间体。中毒的计算对象是毒物-表面复合态,包括毒物原子与金属中心的键长、毒物覆盖度、相邻吸附物排布,以及表面原子是否发生明显弛豫。

DFT 处理这类复合态时,常把干净表面、毒物吸附态和反应物吸附态放在同一套泛函、同一超胞、同一覆盖度口径下比较。Eads 给出热力学占位倾向,结构弛豫给出局域配位改变,电子结构图谱给出轨道和电荷响应,三者共同限定毒物是否真正削弱反应位点。
“毒物”这个词并不要求材料完全失去所有活性。一个 CO 分子占据 Pt 位点时,其他位点仍可能吸附 H*、O2 或含氧中间体;但若毒物覆盖度升高,可用位点数、反应物进入几率和后续脱附窗口都会改变。计算文章讨论中毒时,应把毒物覆盖的局部模型写明,而不能把单个位点结果扩展成整颗催化剂的失效结论。
催化剂中毒还分可逆和持久两种情况。弱配位分子在温度、电位或气氛变化后可能脱附;S、Cl 或强配位 CO 在特定金属位点上形成深能谷时,清除过程会牵涉较高脱附能垒、氧化路径或表面重构。能量差只描述一个状态的相对高低,再生难度还要由脱附路径、氧化路径和工作态结构共同约束。




在 Pt 等金属位点附近,毒物分子进入吸附口袋后最先改写竞争吸附排序。以 Pt 位点为例,CO 的 C 端靠近金属中心时,Pt-C 相互作用通常较强,Eads 可明显低于许多反应物或弱吸附分子。强吸附毒物占据同一金属中心后,反应物分子可选的构型空间随之收缩。

CO、O2 和 CO2 等分子在同一 Pt-graphene 模型上的吸附能差异,能把“抢占位点”写成能量排序。CO 的 V’C-mode 构型中,C 端与 Pt 发生近距离相互作用;若同一位点原本承担 O2 活化或 H* 形成,毒物吸附态越深,反应物替换毒物所需的热力学驱动力越高。

竞争吸附还受到覆盖度约束。低覆盖度下,毒物只占据少数强吸附位;覆盖度升高后,毒物-毒物侧向相互作用、局部偶极和表面弛豫会改变每个新增毒物的增量吸附能。同一毒物在孤立吸附和高覆盖吸附下,对反应物的排斥程度可能不同,这个差异常被微观动力学模型转成位点可用率。
当 CO、S 或 Cl– 覆盖金属中心后,剩余反应物会转向邻位、桥位或缺陷位并形成新的构型偏好。CO 覆盖的 Pt 表面会把水、OH–、H* 或含氧中间体推向邻位、桥位或缺陷位;这些新位置的键长和配位数与干净表面不同。中毒后的活性下降常来自“可用位点减少”和“剩余位点能量重排”两部分,吸附能表格里应区分毒物位点、反应物位点和共吸附位点。
若毒物分子参与后续反应,吸附过深并不总是坏结果。例如 SO2 的裂解、CO 的氧化清除或含氯物种的电化学脱附,都可能经过毒物-表面键的转化。DFT 的任务是分辨占位态、转化态和脱附态的相对能量,再把这些状态与温度、电位和气相化学势对应起来。




毒物覆盖后,电子结构最先反映在表面局域态、诱导电荷和界面势分布中。CO 终止 Pt(111) 的计算显示,CO 层附近的诱导电荷峰位和裸 Pt 表面不同;毒物层改变的是界面响应空间位置,这会影响电化学界面中的电容、局部电场和吸附中间体稳定。

平面平均电荷密度适合描述界面法向上的电子响应,但它给出的是空间分布,和某个原子的形式价态不同。若要把电荷响应分配到 Pt、C、O 或吸附物区域,Bader 或 Hirshfeld 分区可以提供积分数值;分区电荷必须和差分电荷密度的空间峰位对应,否则数值变化容易被误写成孤立原子得失电子。

Pt-CO 界面的逆电容剖面说明,CO 层并非只是几何盖帽;金属区、吸附层、CO-溶剂间隙和扩散层各自贡献不同。对电催化中毒而言,毒物覆盖态会改变双电层附近的电势降和局域响应,吸附能之外的电场修正、显式水层和电位条件会改变结论范围。
S 与过渡金属纳米团簇形成 TM-S 键后,DFT 图谱常呈现强化学吸附特征,Eads、ΔEint、ΔEdist、配位数和平均 TM-TM 键长都会变化。若结构畸变能和 d 态中心同步移动,毒物吸附已经越过简单占位,进入金属骨架和电子态重排。

d 态中心 εd 在中毒分析中只是一个电子结构描述符。S 或 CO 吸附后,金属 d 态与毒物 p 轨道发生杂化,费米能级附近态密度和反键态占据都可能变化;PDOS、COHP/ICOHP 和键长变化可共同判断 TM-S、Pt-C 或 M-Cl 键是否削弱后续中间体吸附。
单张差分电荷密度图给出的判断范围有限,尤其在金属表面和小团簇中,电子离域会让“得电子”与“成键增强”分开。材料计算文章写催化剂中毒时,适合把毒物吸附能、毒物-金属键长、表面 d 态、反应中间体自由能和脱附路径并列列出,每一项对应不同的物理对象。




抗中毒判断不宜追求毒物吸附能一味变弱。毒物吸附过弱时,含硫或含氧分子可能无法被活化;吸附过深时,活性位会长期被覆盖。合理的计算目标是形成可转化、可清除、不过度占位的吸附态,并让目标反应中间体仍能保持合适的能量窗口。
SO2 在 Ti13、Zr13 和 Hf13 团簇上的结果说明,含硫分子可能从分子吸附转入解离吸附。解离态中 TM-S 与 TM-O 键同时出现,PDOS 中金属 d 态、S p 态和 O p 态发生重叠;这种毒化状态的再生难度由断键路径和氧化/还原环境共同限制。

判断抗中毒能力时,模型条件必须写实:slab 厚度、缺陷位、表面晶面、毒物覆盖度、共吸附水层、电位修正、自旋态和 U 值都会改变吸附能排序。同一毒物在台阶位、空位、单原子位和金属纳米团簇上的配位数不同,同一 Eads 数值不能跨模型直接比较。
再生过程的计算对象应从毒物吸附态转到路径态。CO 清除常涉及 OH* 辅助氧化或 CO2 释放,S 毒化可能涉及 SOx 形成、硫迁移或金属-硫键断裂,Cl– 相关中毒还受电位和溶剂化控制。脱附能、NEB 能垒、反应自由能和覆盖度演化分别约束热力学、动力学和位点占有率。
在热催化或电催化工作条件下,毒物覆盖度会随气氛、电位和温度变化而改变,同一材料可能在低覆盖度保留邻位反应,在高覆盖度转入整体失活。吸附等温线、覆盖度自由能和微观动力学可以把 Eads 转成位点占有率,避免把一个 0 K 构型当成工作态比例。
最终写成计算结论时,可以把催化剂中毒分成三类信号:毒物占位信号、电子结构失配信号和再生路径信号。毒物占位由 Eads 与构型给出,电子结构失配由 PDOS、COHP、差分电荷密度和 εd 给出,再生路径由脱附或氧化能垒给出。这三类信号分别对应位点数量、键合状态和工作态寿命。
