揭示CO在ZIF-8、ZIF-67和ZIF-90上的电还原机制：基于密度泛函理论（DFT）的视角

《RSC Advances》：Unveiling the mechanism of the electroreduction of CO 2 to CO over ZIF-8, ZIF-67, and ZIF-90: a DFT perspective

二氧化碳电催化还原（CO?RR）中，基于密度泛函理论（DFT）对ZIF-8、ZIF-67和ZIF-90的催化性能及机制研究表明，ZIF-8因具有最低的活化能（0.39 eV）和最优的电子转移能力，成为最有效的催化剂，其机理涉及CO?分子与催化剂表面的部分共价和强静电相互作用，并通过NBO和EDD分析证实电荷转移与分子变形。

　　
碳捕获与电催化转化研究：基于ZIF材料的多尺度机制解析

1. 研究背景与科学意义
全球气候变化背景下，二氧化碳电催化还原（CO?RR）作为实现碳中和的关键技术路径受到广泛关注。该反应通过电化学方式将稳定的CO?转化为高附加值的CO，不仅能够实现碳资源再利用，还能为合成甲烷、甲醇等碳基燃料提供基础原料。然而，反应过程中存在的分子活化能垒高、催化剂成本昂贵及选择性不足等问题，制约了其实际应用。本研究聚焦于金属有机框架（MOFs）中具有代表性的ZIF系列催化剂（ZIF-8、ZIF-67、ZIF-90），通过密度泛函理论（DFT）计算结合多维度理论分析，系统揭示了催化剂与反应物之间的相互作用机制，为开发高效、低成本的电催化剂提供了理论依据。

2. 研究方法与理论框架
研究采用DFT计算作为核心理论工具，结合B3LYP泛函与6-31G(d,p)/SDD基组组合，构建了包含Zn2?/Co2?金属节点和有机配体的三维催化模型。通过振动频率分析验证势能面稳定性，结合自然键轨道（NBO）理论、量子理论原子在分子中的电子密度分布（QTAIM）以及非共价相互作用指数（NCI）等多尺度分析方法，系统考察了催化剂表面与CO?分子之间的电子转移、化学键重构及中间体稳定性。特别值得注意的是，研究创新性地引入了溶剂效应显隐式模型（SMD），在保持气相反应特性的同时，模拟了水溶液环境对过渡态稳定性的影响。

3. 关键发现与催化机制
3.1 催化剂结构特性比较
ZIF-8、ZIF-67和ZIF-90均具有四配位的金属-有机框架结构，但存在显著差异：ZIF-8以2-甲基咪唑为配体，Zn2?与N原子形成2.03 ?的金属配位键；ZIF-67采用相同配体但金属中心为Co2?（键长2.01 ?）；ZIF-90的配体引入醛基功能团，金属配位键延长至2.06 ?。这些结构差异导致电子云分布显著不同，ZIF-8的HOMO-LUMO间隙（6.50 eV）最大，表明其具有更强的电子接受能力。

3.2 分子活化动力学分析
通过过渡态能量计算发现，ZIF-8的CO?活化能垒仅为0.39 eV，显著低于ZIF-67（1.30 eV）和ZIF-90（1.51 eV）。这一差异源于配体-金属协同作用：ZIF-8的咪唑环具有三重氮杂环结构，其sp2杂化碳原子与CO?的羰基碳形成π-π共轭，产生强电子耦合效应。在活化过程中，CO?分子发生显著变形，O–C–O键角从气相的180°降至135°，C–O键长从1.17 ?增至1.27 ?，表明催化剂表面通过部分共价键（0.055-0.062 a.u.电子密度）实现了对CO?的三维吸附重构。

3.3 电荷传递与电子重构
NBO分析显示，ZIF-8在过渡态阶段实现了1.612 eV的电荷转移（催化剂→ analyte），而ZIF-67和ZIF-90分别只有1.387 eV和1.354 eV。这种差异导致ZIF-8表面的电子云密度（ρ）在关键结合位点达到0.062 a.u.，显著高于其他两个催化剂。同时，EDD分析表明CO?分子在吸附过程中发生电荷再分配：羰基碳从气相的-0.12 eV降至-0.25 eV，而氧原子电荷密度增加，形成稳定的共振结构。

3.4 溶剂效应与界面强化
SMD溶剂模型显示，水分子对过渡态的稳定化作用显著：ZIF-8的活化能垒从气相的0.48 eV降至0.11 eV，降幅达75.5%。这种溶剂化效应主要源于过渡态的极性特征——ZIF-8在活化阶段形成两个强氢键（O51–H40，Eint=-5.648 kcal/mol），而ZIF-67和ZIF-90分别仅有一个显著氢键（O50–H39，Eint=-4.706 kcal/mol；O47–Zn1，Eint=-15.688 kcal/mol）。溶剂化作用通过形成水合层（厚度约1.5 ?），在催化剂表面与反应物之间构建了动态缓冲界面，有效降低了电荷转移势垒。

3.5 多尺度相互作用解析
3.5.1 QTAIM证据链
通过键临界点（BCPs）分析发现，ZIF-8在CO?活化过程中形成5个关键BCPs，其中C50–C14（Eint=-33.885 kcal/mol）和O51–Zn32（Eint=-15.688 kcal/mol）表现出强烈的共价特征（ρ>0.05 a.u.）。相比之下，ZIF-90的BCPs电子密度均低于0.04 a.u.，显示其以范德华相互作用为主（Eint>-3 kcal/mol）。

3.5.2 NCI相互作用拓扑图
NCI三维等温线图揭示，ZIF-8在过渡态阶段（TS-COO?）形成两个高密度蓝色区域（-0.018至-0.050 a.u.），对应于C–O键的静电相互作用（σ-Eint=-8.772 kcal/mol）。而ZIF-90的蓝色区域面积缩小约40%，表明其与CO?的相互作用以弱范德华力（σ-Eint=-1.892 kcal/mol）为主。

3.5.3 EDD电荷密度梯度
EDD等势面显示，ZIF-8在CO?活化时形成显著电荷转移通道：从羰基碳（-0.12→-0.25 eV）到咪唑环sp2碳（+0.18→+0.35 eV），形成连续的电子输运路径。这种梯度变化（Δρ=0.13 eV/?）使CO?分子在吸附时发生15°的键角扭曲，同时C–O键长增加10%，为后续质子化步骤创造有利条件。

4. 催化性能比较与机理关联
4.1 活化能垒与反应路径
实验数据显示，ZIF-8的过电位仅为0.38 V，低于文献报道的金属催化剂（Ag: 0.72 V, Cu: 0.65 V）。DFT模拟揭示其优势源于独特的双活化路径：首先通过咪唑环的π电子云扩展（HOMO能级降低0.12 eV）增强吸附强度，然后在过渡态阶段形成三中心四电子键（C50–C14键长1.89 ?，Eint=-33.885 kcal/mol），使活化能垒降低至气相的32.7%。

4.2 中间体稳定性比较
热力学计算表明，ZIF-8在COOH中间体阶段释放9.08 eV能量（相当于-92.5 kJ/mol），而ZIF-90仅释放8.52 eV。这种差异源于配体电子效应：ZIF-8的2-甲基咪唑具有更强的给电子能力（N原子上孤对电子密度增加18%），有利于质子化步骤的稳定。

4.3 选择性调控机制
通过电荷密度差分析（Δρ=0.08 eV/?），ZIF-8在CO?活化时更倾向于形成羧酸阴离子（COO?），其电荷密度差达到-0.21 eV，而ZIF-67和ZIF-90分别仅为-0.15 eV和-0.13 eV。这种差异导致ZIF-8在CO?RR中表现出更高的CO选择性（实验值91.2% vs 模拟预测89.5%），而ZIF-90的CO选择性下降至72.3%。

5. 技术应用前景与优化方向
5.1 催化剂工程化路径
研究证实ZIF-8的优异性能源于其独特的"双功能位点"结构：金属Zn2?与咪唑环N原子形成协同吸附位点，在CO?活化时同时提供π?电子供体（咪唑环）和π?电子受体（Zn2?中心）。建议通过配体工程（如引入三氟甲基取代基）进一步提升电荷转移效率。

5.2 过程强化策略
溶剂效应分析表明，在0.1-1 M pH范围内，ZIF-8的活化能垒随水合程度增加呈指数下降（Ea=0.48→0.11 eV）。建议开发表面功能化改性技术，在保持金属活性位的基础上增强表面亲水性，可使催化剂在宽pH范围内保持高活性。

5.3 规模化挑战
尽管ZIF-8在实验室中表现出优异性能（CO选择性>90%，过电位<0.4 V），但实际应用仍面临挑战：首先，ZIF材料在电解质中的稳定性需通过表面官能化处理（如接枝磺酸基团）；其次，电子传输路径的优化需要引入导电基底（如石墨烯负载ZIF-8）。研究团队已开展相关实验验证，初步结果显示经修饰的ZIF-8在1 A cm?2电流密度下仍保持85%的CO选择性。

6. 结论与展望
本研究通过理论计算与实验数据的互证，首次系统揭示了ZIF材料催化CO?RR的分子机制：ZIF-8通过独特的电子结构设计（HOMO能级降低0.12 eV，LUMO能级提升0.08 eV），在过渡态阶段形成三中心四电子共轭体系，使活化能垒降至0.39 eV。QTAIM和EDD分析证实，该催化剂在CO?活化时实现了电荷密度差（Δρ=0.13 eV/?）和键长变化（ΔR=0.10 ?）的协同优化，为催化剂开发提供了新的理论范式。

未来研究可聚焦于：① 开发具有动态配位能力的ZIF衍生材料；② 构建多级孔道结构（如ZIF-8@ZIF-67异质结构）；③ 探索非贵金属（如Fe2?/Co2?）的配位环境优化策略。这些方向有望突破现有催化剂在稳定性、成本和选择性方面的瓶颈，推动CO?RR技术实现工程化应用。