说明:本文华算科技主要介绍缺陷数量与催化活性的计量方式、缺陷间相互作用对单个位点反应的影响、DFT 超胞中的缺陷浓度含义,以及确定工作态缺陷区间所需的计算和实验信号。


缺陷数量与活性的基本含义是什么?
“缺陷数量”必须带有明确的计数基准。点缺陷可写成缺陷原子数占可占据格点数的比例 cdef = Ndef/Nsites,二维表面也常用单位面积缺陷数;非化学计量氧化物用 δ 表示平均缺氧程度;台阶、位错和孔洞则需要长度密度或面积分数。同一个百分数用于体相、表面和晶界时,代表的原子环境并不相同。
几何电流密度、质量活性、比活性和 TOF 使用不同的面积、质量或位点分母。几何电流密度包含电极面积、负载量、孔道利用率和传质贡献;质量活性除以催化剂质量;比活性通常按 BET 面积或电化学活性面积归一;周转频率 TOF 则要求估算可参与反应的位点数。用一个简化乘积表示,总速率约等于可接近位点数乘以单个位点速率。缺陷增多可以提高前一项,而单个位点速率可能下降。
Pt 纳米管经过不同时间的酸刻蚀后,裂纹和孔洞面积增加,甲醇氧化峰电流却在中等刻蚀程度达到较高值。横轴的相对缺陷密度与纵轴的峰电流密度采用各自的测量口径,图中的峰形关系适用于这组 Pt 纳米管和相同测试条件,不能替代其他材料中的缺陷计数。

真实材料中的缺陷含量通常由多种信号约束。LaxSr1−xCoO3−δ 系列里,平均氧缺失可由化学组成估算,原子分辨 ABF-STEM 给出局域 O 柱强度,EXAFS 拟合得到 Co—O 配位数,O K 边 EELS 记录未占据态变化。δ、局域 O 柱强度和 Co—O 配位数描述的是不同空间尺度;δ 随处理时间升高而 Co—O 配位数降低时,缺陷梯度才有明确的原子结构含义。

位于颗粒内部的空位、被吸附物占据的孔洞、已经聚并成孔边的缺陷,并不会贡献相同数量的可用位点。名义缺陷含量、表面可接近缺陷数和工作态活性位数应分别登记;用 XPS 原子比直接除以电流得到“单缺陷活性”,会把探测深度、表面污染和位点占用误差写进 TOF。


缺陷相互作用怎样改变反应速率?
孤立空位打断局域配位,使邻近金属原子出现配位不饱和,并改变吸附物与金属 d 轨道的杂化。相邻缺陷距离缩短后,局域应变场、缺陷态和吸附物覆盖开始相互耦合;新增缺陷不再复制原有位点,而会生成台阶簇、孔边或新的表面终止。
Pt 纳米管研究把甲醇脱氢生成 COads 的速率常数记为 k1,把 COads 氧化速率常数记为 k2。管中心在所测范围内随缺陷密度增加而逐步加快,初始缺陷已经较多的管端则出现先升后降。同一颗粒的中心与端部位于缺陷—速率曲线的不同区段,面积平均电流会掩盖这类空间差别。

缺陷带来的吸附增强若主要稳定 COads,甲醇脱氢可能加快,后续 CO 氧化却承受更高的表面占用。过量低配位 Pt 原子还会提高吸附诱导原子迁移的幅度,使工作表面持续重排。反应物活化与毒化中间体清除使用不同过渡态,单个位点很难让两项速率随缺陷密度保持同一变化方向。
DFT 模型中的 Pt(111)、Pt(211) 以及提高台阶缺陷比例的 Pt(211)-50% 和 Pt(211)-75%,分别给出甲醇脱氢与 CO 氧化的自由能序列。最大过渡态自由能随缺陷比例呈非单调变化;缺陷丰富的 Pt(211)-75% 并未让两条反应支路都取得最低值。

氧化物中还会出现另一类浓度效应:空位改变邻近金属—氧键、O 2p 与金属 3d 态的位置,并可能改换反应中参与供氧的原子。LSCO 数据覆盖多种 La/Sr 比和氧缺失量,BET 归一电流在同一氧缺失量附近仍有较大离散,AEM 与 LOM 区域也发生转换。缺陷种类相同并不保证活性位电子态相同,A 位组成和局域 Co—O 结构仍在调节反应。



DFT的缺陷浓度为何会改变计算结果?
周期性 DFT 中,一个含 N 个可占据格点的超胞删除一个原子,模型浓度就是 1/N;所有周期镜像都会重复同一缺陷。缩小超胞既提高缺陷比例,也缩短缺陷镜像间距,两项改动共同影响弹性弛豫、静电势、缺陷态展宽和吸附物相互作用。
单个缺陷的稀释极限适合研究近似孤立位点;高浓度模型描述的是相邻缺陷已经相互作用的材料状态。带电缺陷还会受到补偿背景、电荷修正和费米能级的影响。“一个空位”的计算结果没有脱离超胞尺寸的独立含义,同一泛函下至少应给出缺陷间最短距离和超胞中缺陷比例。
MoS2 基面模型从 0% 到 21.9% 硫空位,既改变单位晶胞表面自由能,也改变空位排列。连续生成的空位在该计算中更倾向于靠近已有空位,并形成锯齿状团簇;H* 吸附自由能随浓度进入接近热中性的区间。浓度序列包含数量变化,也包含构型变化,均匀分散空位与团簇空位应作为两个模型集合比较。

比较浓度效应时,总能必须换成同一参考态下的每缺陷形成能、每面积表面自由能或每位点吸附自由能。若高浓度模型多删除了两个原子,直接比较两个超胞总能差会把原子数变化写成“缺陷效应”。归一方式、化学势与电荷态必须保持一致,弛豫范围和磁性初态也应在各浓度模型中采用相同方案。
缺陷态在高浓度下可能由离散局域态扩展成缺陷带,邻近原子的 p/d 态峰位和占据也会重排。LSCO 的计算模型把 VO 浓度从约 1.8% 提高到 50%,Co 3d 与邻近 O 2p 的 PDOS 并非简单增高或降低;50% VO 已接近另一种配位网络。高浓度模型描述缺陷有序化或新化学计量时,结论不应冠以“孤立空位”。


怎样确定可工作的缺陷区间?
一组可信的浓度序列要在相同晶相、颗粒尺度、负载量和电解条件下改变缺陷含量。若球磨时间又改变比表面积和相组成,几何电流的上升无法分配给缺陷本身。缺陷含量与 BET/ECSA 归一活性应取自同一批次、同一反应前后状态,否则横轴和纵轴代表的是两份不同材料。
LSCO 的 36 组样品显示,氧缺失量升高后,若干组成的 OER 电流进入低值区。计算得到的 O 2p 带、Co 3d 轨道能级间隔、局域原子距离和 OER 反应能都随 VO 浓度变化。缺陷密度只能作为同一材料家族内的控制变量;La/Sr 比改变后,电子占据和 OER 机理也随之改变。

Pt 纳米管在连续甲醇氧化中,荧光产物热点由初始管端向中心迁移,最优缺陷区间随电位和反应时间变化;反应后的 HRTEM 在管端和中心都记录到新增无序区。反应过程中生成的缺陷会改变原先标定的缺陷密度,初始结构与工作结构不能共用一个横坐标值。

工作态缺陷区间可由四组具体结果相互约束:原位或反应后结构给出缺陷是否聚并、钝化或溶出;ECSA/TOF 区分位点数量与单个位点速率;DFT 浓度与构型序列比较吸附自由能、过渡态自由能和缺陷间相互作用;耐久测试记录活性峰值能维持多久。若实验峰值只出现在几何电流,而 TOF、结构和计算能量均无同向变化,结论应限定为该电极制备条件下的表观活性提升。
可报告的“最优缺陷浓度”应附带材料组成、缺陷类型、空间位置、反应电位、归一分母和测试时长。例如 9% 的基面硫空位、0.3 的体相缺氧 δ、每平方厘米的表面位点数,三者不可互换;同一材料在 H*、OH* 或 COads 覆盖下,也可能出现不同的峰值位置。
