关键词： 双原子催化剂，CO PROX，微反应动力学，高通量计算

以清洁能源 H ₂ 为燃料的质子交换膜燃料电池应用前景广阔，而经蒸汽重整、水煤气变换等制备 H ₂ 残留的微量 CO 将严重毒化 Pt 电极。富氢氛围中 CO 优先氧化（ PROX ）是进气路上在线去除痕量 CO 的理想方式，氧化铈等金属氧化物负载高分散过渡金属体系被认为是极具前景的 PROX 催化剂。

此前，华中科技大学单斌教授团队（微纳材料设计与制造研究中心）基于选择性 ALD 方法制备的 Pt ₁ /CoO _x 、 Pt-Cu/CeO _x 等精细纳米结构表现出优异的 CO 氧化活性和选择性（ J. Mater. Chem. A , 2020 ,8, 10180; Nat. Commun. , 2020 , 11, 4240 ）， DFT 计算结合微反应动力学分析阐明了 CO 与 H ₂ 竞争氧化机制（ ACS Appl. Mater. Interfaces , 2022 , 14, 27762 ）。与单原子催化剂（ SAC ）相比，双原子催化剂（ DAC ）不仅具备高原子利用率和高选择性等优势，同时有望通过协同效应和电子调控效应等打破线性比例关系极限，实现 PROX 性能的进一步提升。对此，单斌教授团队与云南省贵金属新材料集团有限公司张爱敏首席研究员团队合作，在 ACS Catalysis 期刊发表题为“ Synergistic Dual-Atom Catalysts on Ceria for Enhanced CO Preferential Oxidation: Insights from High-Throughput First- Principles Microkinetics ”的研究文章 （点击文末「阅读原文」，直达链接） 。

图 1. 负载型双原子 PROX 催化剂高通量搜索示意图

该项工作建立了包含了双原子模型构建、稳定性评估、反应网络 DFT 计算、微反应动力学分析和描述符识别等模块的高通量工作流。首先，以 CeO ₂ (110) 为载体构建了包含 3d, 4d, 5d 过渡金属和 5 种氧配位环境组合的 2325 种负载双原子构型，如图 1 所示。进而通过不同机器学习方法建立稳定性预测模型，其中基于压缩感知的符号回归方法（ Sure Independence Screening and Sparsifying Operator, SISSO ）表现出良好的预测效果，如图 2 所示。 SISSO 推荐的 3D 描述符表明，双金属原子的 D 电子数之和、电负性倒数之和以及氧配位数是影响稳定性的关键，且 D 电子数更少、电负性更低以及氧配位数适中对应更稳定的 DACs 。

图 2. 负载过渡金属双原子构型的稳定性分析

进而面向所有同质 DACs （ M _A =M _B ），如图 3 所示，高通量计算 CO 直接氧化、 H ₂ 直接氧化、 OCOO 、 COOH 、氧迁移和 H （ OH ）辅助等可能发生路径中的关键基元步骤。结果表明，以适中氧配位数（ N _O =2 ， * ）为原点的包含 OCOO 、 O 迁移和直接氧化等路径构成的 MvK 催化循环被认为是主导机制，即 O ₂ 预解离的 MvK 机理。为验证该机理并量化分析影响 CO 氧化活性和 CO ₂ 选择性的关键因素，以实验报道的 Ir ₁ -Ir ₁ DAC 作为代表，基于 CO PROX 全反应网络的 DFT 计算能量信息，在实验条件下进行微反应动力学分析。

图 3. 包含主要路径的 CO PROX 简化反应网络示意图

为实现对 DACs PROX 性能的整体评估，基于团队提出的“态到态”方法建立了复合微反应动力学模型。受负载双原子影响而显著区别于干净 CeO ₂ 表面位点的区域被视为双原子限域位点，并对该区域进行整体化处理，即将限域位点中所有局域环境不同的可能共吸附组合视为相互独立的反应中间态，将基元反应视为相邻反应态之间的跳转。其中，氧迁移路径通过含氧物种的扩散沟通了双原子限域位点和干净 CeO ₂ 表面位点（平均场近似）。如图 4 所示，速率控制因子和选择性控制因子的计算结果表明， CO 直接氧化（ * → OO _s * ）和 OCOO （ OO _s * → O _V * ）路径主导了低温 CO 氧化活性， H ₂ 直接氧化（ * → OO _s * ）和氧迁移路径（ OO _s * → O _V * ）是导致高温选择性下降的关键。进一步结合速率贡献度分析，将主导路径涉及的所有基元反应组合为简化 PROX 网络。与全反应网络求解结果的对比分析表明，简化网络能够实现对 CO 氧化速率、 H ₂ 氧化速率和 CO ₂ 选择性在 PROX 温度窗口内的近似描述。

图 4. Ir ₁ -Ir ₁ DAC 的 PROX 全反应网络的微反应动力学分析

随后，针对所有同质 DACs 进行简化网络的批量 DFT 计算， SISSO 方法推荐的 PROX 性能描述符指向解离 O ₂ 的氧去除能（或氧空位形成能， E _Ov ）、 CO 和 H 在金属位点的吸附能，其与 CO 氧化活性和选择性的关系如图 5 所示。 H 吸附能适中（ -0.6~0.0eV ）时， CO 吸附能介于 -2.0~-1.5eV 、氧空位形成能介于 2.6~3.2eV 是兼具良好 CO 氧化活性和 CO ₂ 选择性的理想范围，且该区间可随 H 吸附的减弱而进一步扩大。 Ir ₁ -Ir ₁ 和 Ru ₁ -Ru ₁ 等贵金属 DACs 在同质 DACs 中接近活性火山图的顶端，同时也展现出良好的选择性，与实验报道相符，但 PROX 性能仍具提升空间。

图 5. H 原子在金属位点弱结合和强结合下的 CO 氧化活性和选择性

而后，面向所有 DACs 针对 PROX 性能描述符进行高通量计算，从而预测其 CO 氧化活性和 CO ₂ 选择性，如图 6 所示。部分金基和贵金属基 DACs ，如 Fe-Au 、 Cu-Au 和 Zn-Ru 等，展现出优异的低温 PROX 性能。其中， Fe-Au 、 Fe-Ru 等 6 种 DACs 在整个工作温度窗口均保持 TOF>0.1s ^-1 且选择性 >50% ，被认为是 DACs/CeO ₂ 中最具潜力的 PROX 催化体系。最后以 Fe-Au 为例，分别针对 Fe-Fe 、 Au-Au 和 Fe-Au 进行微反应动力学对比验证和电子结构分析。结果表明 Fe-Au 兼具了 Au-Au 高活性和 Fe-Fe 高选择性的优势，这源于 Au 、 Fe 对活性氧结合强度以及氢吸附强度的有效协同调控。

图 6. 基于描述符的 PROX 性能预测与微反应动力学验证

综上所述，本项研究基于团队提出的态到态复合方法，面向金属氧化物负载型过渡金属双原子体系的 CO PROX 过程，开展了第一原理计算结合微反应动力学分析的机器学习辅助高通筛选工作，发现了数种潜在高效的双原子 PROX 催化剂，为负载型 DACs 的理性设计和性能调控提供了理论指导准则。

单斌教授 ： 华中科技大学材料科学与工程学院教授、博士生导师、兼任中科院宁波技术与工程材料所客座研究员。教育部新世纪优秀人才支持计划获得者，湖北省特聘教授，湖北省杰出青年基金获得者，中国稀土学会催化专业委员会委员。主持和参与了包括国家自然科学基金重点项目、面上项目，国家重大科学研究计划项目在内的多项国家级及省部级项目。在 Science , Nat. Commun. , Angew. Chem. Int. Ed. , Phys. Rev. Lett. 等国际权威杂志发表论文 150 余篇，他引超 8000 余次；获发明专利授权 40 余项。获日内瓦国际发明展特别嘉许金奖、湖北省技术发明一等奖、中国稀土学会科学技术二等奖、湖北省专利银奖等奖项。畅销教材《计算材料学 - 从算法原理到代码实现》第一主编。

主要研究方向：计算材料学、先进催化与原子层沉积的跨尺度模拟、人工智能在材料中的应用。

张爱敏研究员： 云南省贵金属新材料控股集团有限公司首席研究员，昆明理工大学博士生导师，云南省政府特殊津贴获得者。现任九三学社云南省第九届委员会新社员工作委员会副主任，九三学社昆明贵金属研究所支社副主委。先后主持国家高技术研究发展计划（ 863 计划）、国家重点研发计划、国家工信部科技成果转化、国家自然科学基金、云南省重大科技专项在内的课题 16 项，获省部级科技奖 7 项，获国家专利授权 6 项，获软件著作权 5 件，合作出版著作 3 部，以第一或通讯作者在国内外重要刊物发表论文 40 余篇。

主要研究方向：工业催化剂研发与产业化、材料基因工程建设及关键技术研究。

刘璋博士研究生 ： 华中科技大学环境科学与工程学院 2020 级在读博士研究生。目前在 ACS Catal. , Appl. Catal. B-Environ. , Chem. Mater. , J. Mater. Chem. A 等期刊发表论文 20 余篇。

主要研究方向：面向催化反应和原子层沉积过程的跨尺度模拟。

更多科研作图、软件使用、表征分析、SCI 写作、名师介绍等干货知识请进入后台自主查询。