类型a

本研究由L. C. Grabow和M. Mavrikakis完成，两位作者均来自美国威斯康星大学麦迪逊分校（University of Wisconsin, Madison）化学与生物工程系。该研究于2011年3月4日在线发表于ACS Catalysis期刊，文章标题为“Mechanism of Methanol Synthesis on Cu through CO₂ and CO Hydrogenation”。研究得到美国化学会（American Chemical Society）出版支持，属于催化化学与表面科学领域的基础研究。

学术背景
甲醇合成是工业催化中极为重要的反应之一，长期以来以Cu/ZnO/Al₂O₃为催化剂，在230–280 °C、50–120 atm条件下由合成气（CO/CO₂/H₂）大规模生产甲醇。然而，尽管该技术已有近百年历史，关于甲醇合成的碳源究竟是CO还是CO₂、反应遵循何种机理路径、活性位是金属Cu还是Cu⁺等关键科学问题始终存在争议。早期研究认为CO加氢是主要路径，但Chinchen等人的同位素标记实验指出CO₂是主要碳源，随后大量研究倾向于支持CO₂加氢机理。与此同时，Cu/ZnO/Al₂O₃催化剂中ZnO载体是否对Cu有协同促进作用、反应是否具有结构敏感性等问题也未得到统一认识。已有多种动力学模型试图描述甲醇合成反应，但这些模型往往基于不同机理假设和速控步骤假定，彼此之间存在明显分歧。针对这一研究现状，作者旨在建立一个全面、系统的微观动力学模型，不预先假定速控步骤或反应路径，而是基于密度泛函理论（Density Functional Theory, DFT）系统考察CO和CO₂加氢以及水煤气变换（Water-Gas-Shift, WGS）反应网络，从而定量阐明工业条件下甲醇合成的主要反应机制。

研究流程与详细方法
本研究的工作流程可以划分为以下主要环节：DFT计算基础参数获取、反应网络构建、微观动力学模型开发、参数拟合与优化、机理分析与结果验证。

在DFT计算阶段，作者选用Cu(111)表面作为模型催化面，采用DACAPO总能计算程序包，以三层Cu原子平板模型（slab）周期性重复模拟表面结构，真空层厚度约10.6 Å。计算采用p(3×3)超胞，对应每物种覆盖度为1/9单层。使用超软赝势（ultrasoft Vanderbilt pseudopotentials）描述核电子相互作用，平面波截断能为25 Ry，表面布里渊区以54个k点采样，交换关联泛函选用PW91广义梯度近似（GGA-PW91）。通过迭代对角化Kohn-Sham哈密顿量获得自洽电子密度，并将总能外推至电子温度kBT=0 eV。在频率计算中，采用谐振子假设，通过内坐标下质量加权Hessian矩阵对角化分离振动模式与平动/转动模式。所有基元反应的最小能量路径均采用爬坡图像微动弹性能带方法（CI-NEB）计算，每个路径至少由七个中间构型离散化描述，过渡态经振动频率计算确认具有单一虚频。研究共计算了22种吸附物种和8种气相物种的性质以及49个基元步骤的反应能量学数据。

在反应网络构建方面，作者纳入了CO₂加氢、CO加氢、WGS反应以及副产物甲酸（HCOOH）、甲醛（CH₂O）和甲酸甲酯（HCOOCH₃）的生成路径。模型包含49个基元步骤，涉及的反应中间体包括甲酸根（HCOO*）、甲酸（HCOOH*）、羟甲氧基（CH₃O₂*）、甲醛（CH₂O*）、甲氧基（CH₃O*）、甲酰基（HCO*）、羧基（COOH*）、二氧亚甲基（H₂CO₂*）等。为描述中间体之间的构型与能量差异，作者对每个吸附物种在Cu(111)表面的最优吸附位、结合能（BE）和振动频率均进行了系统计算，并以此作为后续微观动力学模型的输入参数。

在微观动力学模型开发中，作者采用平均场（mean-field）近似，不考虑吸附质局部有序性和扩散限制，并假定所有物种占据一个表面位点，最大覆盖度限制为1单层。熵值由振动频率直接计算，通过Shomate方程拟合温度相关热容（Cp）和焓值。对于每个基元步骤，指前因子由初态与过渡态的熵差计算得到。为在对DFT计算值进行调整时不破坏势能面的基本特征，作者引入了一个参数ω，用以描述过渡态在反应坐标上相对于初态和终态的位置。前向活化能由公式Ef = Ef,DFT + ω(ΔH − ΔE_DFT)计算，其中Ef,DFT和ΔE_DFT来自DFT计算，ΔH为实际反应焓变。通过这种方式，可以在拟合实验数据时对活化能进行系统调节，同时保持与DFT计算得到的Brønsted-Evans-Polanyi线性关系基本一致。

在参数估计与模型优化阶段，作者采用Graaf等人发表的工业Cu/ZnO/Al₂O₃催化剂动力学实验数据对模型参数进行拟合。由于原始DFT参数预测的甲醇合成速率比实验值低约六个数量级，作者首先对WGS相关参数子集进行初步拟合，以Koryabkina等人的WGS动力学数据为参考，使WGS速率达到实验值同一数量级。随后对31个最敏感参数（包括若干中间体结合能的修正值δH_fit和若干基元步骤的ω值）进行最小二乘拟合，目标函数同时考虑了CH₃OH和H₂O的生成速率，两者权重相等。最终优化后的模型可以同时重现CH₃OH和H₂O生成速率，其R²值分别达到0.90和0.94，整体R²为0.92。

主要结果
DFT计算结果显示，在Cu(111)表面上，CO₂加氢生成HCOO的直接路径（CO₂ + H* → HCOO* + *）具有0.87 eV的活化能，且放热0.25 eV；而经由碳酸盐（CO₃*）中间体的路径因CO₂解离吸附的活化能高达1.77 eV而不可行。HCOO*进一步加氢时，在O原子上加氢生成HCOOH*的能垒为0.91 eV，显著低于在C原子上加氢生成H₂CO₂*的1.59 eV能垒，因此模型确定HCOOH*为HCOO*加氢的首选产物。HCOOH*继续加氢生成CH₃O₂*的能垒为1.04 eV，CH₃O₂*随后经C–O键断裂生成CH₂O*与OH*，该步热力学能垒为0.74 eV。CH₂O*加氢生成CH₃O*的能垒仅0.24 eV，而CH₃O*加氢生成CH₃OH*的能垒为1.17 eV，成为整个路径中能垒最高的步骤之一。

微观动力学模型拟合后的机理分析表明，在典型工业条件（T=499.3 K，P=29.9 atm，y_CO=0.053，y_CO₂=0.047，y_H₂=0.90）下，甲醇主要通过CO₂加氢路径生成，约占甲醇总产量的2/3，其余1/3来自CO加氢。CO₂加氢路径的中间体依次为HCOO*、HCOOH*、CH₃O₂*、CH₂O*和CH₃O*；CO加氢路径则经由HCO*、CH₂O*和CH₃O*。表面覆盖度计算显示，CH₃O*和H*是最主要的吸附物种，覆盖度分别为0.21单层和0.18单层，HCOO*覆盖度约0.04单层。速率控制度分析表明，CH₃O*加氢生成CH₃OH是CO和CO₂两条路径共同的慢步骤，而HCO + H* → CH₂O* + * 和HCOOH* + H* → CH₃O₂* + * 则分别控制CO和CO₂路径的相对贡献。WGS反应通过羧基（COOH）中间体进行，主要路径为CO + OH* → COOH* + ，随后COOH + OH* → CO₂* + H₂O*。WGS步骤对甲醇生成速率的速率控制度较小，说明WGS活性对甲醇合成速率影响有限。CO对甲醇合成的促进作用并非通过HCO*作为氢供体的辅助加氢路径实现，而是通过OH*的去除（即WGS反应）间接促进甲醇生成，但该作用远小于CO直接加氢对甲醇的贡献。此外，拟合后的中间体结合能普遍较DFT计算值稳定0.3–0.6 eV，这可能与Cu(111)模型表面与工业催化剂真实活性位之间的差异有关，也可能涉及范德华力或表面重构等因素。

结论与价值
本研究建立了一个基于DFT计算的全面的平均场微观动力学模型，首次同时定量描述了Cu催化甲醇合成中CO₂加氢、CO加氢和WGS反应之间的耦合关系。模型明确指出，在工业条件下CO₂加氢是甲醇的主要来源，约贡献2/3的甲醇产量，而CO加氢贡献其余1/3。HCOO*加氢首先生成HCOOH*而非H₂CO₂*，CH₃O₂*是连接甲酸路径与甲醛路径的关键中间体，CH₃O*加氢是共同速控步骤。该研究在方法论上通过引入ω参数和质量加权内坐标频率计算方法，在保持DFT计算结果基本特征的前提下实现了对复杂反应网络的参数拟合，为多相催化微观动力学建模提供了可借鉴的框架。在应用层面，该模型为优化工业甲醇合成催化剂的反应条件与进料组成提供了理论依据，也为进一步研究Cu基催化剂活性位结构与载体效应奠定了基础。

研究亮点
本研究的亮点主要体现在以下几个方面：其一，首次将HCOOH*和CH₃O₂*作为CO₂加氢路径中的关键中间体纳入微观动力学模型，并证明HCOOH*是HCOO*加氢的优选产物；其二，系统比较了直接加氢与CO辅助加氢路径，排除了HCO*作为氢供体的促进作用，阐明了CO促进效应的真正来源；其三，通过参数拟合将DFT计算与工业催化剂实验数据相连接，模型能够同时重现CH₃OH和H₂O生成速率；其四，分析了Cu(111)模型表面与工业催化剂活性位之间的可能差异，指出更开放的Cu表面（如Cu(110)、Cu(100)或部分氧化的Cu(211)）可能更接近真实活性位结构。