动态蒙特卡罗（DMC）模拟与响应面法（RSM）解析甲醇在Ni(100)上的双路径吸附机理并优化产氢率

《RSC Advances》：Resolving the dual-pathway adsorption mechanism of methanol on Ni(100) via dynamic Monte Carlo simulations and optimizing hydrogen yield with response surface methodologyOpen Access

　　

甲醇在Ni(100)上的吸附模式——是分子吸附还是解离吸附——四十余年来一直是表面化学中未解决的核心问题。研究人员将动态蒙特卡罗（Dynamic Monte Carlo, DMC）模拟与响应面法（Response Surface Methodology, RSM）结合，以厘清这一机理歧义并同时优化产氢率。在350×350网格上的晶格气模型模拟了甲醇分解的全部基元步骤，动力学参数经实验程序升温脱附（Temperature-Programmed Desorption, TPD）谱严格验证（R2=0.9897）及独立等温数据校验。模拟明确表明：仅分子吸附或仅解离吸附均不能复现实验H2与CO脱附曲线；只有同时包含两条通道的双路径模型方可匹配所有可观测量。该发现首次明确解决过渡金属表面化学中这一长期机理争议。采用中心复合RSM设计，研究人员考察温度、压力、时间对产氢率的影响，确定最优条件（379.7 K、3.50 μPa、107.8 s）下产氢率达95.5%，与模型预测（96.8%，95%置信区间内）高度一致。敏感性分析显示动力学参数±20%扰动下产氢率变化<1.5%，印证机理结论稳健。经验证的动力学网络为理性催化剂设计及分层微观动力学建模提供定量基础，推进了镍表面醇分解的基础认知与实用进程。

研究背景方面，氢作为高比能量密度的清洁载体在可持续能源与燃料电池场景中十分关键，甲醇因氢碳比高、易储运、可由生物质或CO₂加氢再生而受关注，其催化分解CH₃OH→2H₂+CO成为单晶与多晶镍基催化剂研究重点。Ni(100)是界定基元步骤的模型面，但甲醇在其上的初始吸附长期存在争议：部分TPD与振动光谱支持非解离分子吸附CH₃OH(ad)，另一派据低温即出现甲氧CH₃O与表面氢/氘主张O–H或O–D断裂的解离吸附CH₃O+H。密度泛函理论（Density Functional Theory, DFT）能给单步能垒却无法给出程序升温或等温下的宏观产物分布，实验TPD只给多路径集成结果而缺分子级分辨，因而四十年无定论。研究人员因此用DMC显式追踪晶格上吸附态时空演化，再用RSM做条件优化，目标为辨机理、建网络、找优参。

主要关键技术方法上，研究人员采用埃因霍温理工大学CARLOS程序在Linux下做DMC，以350×350周期平方晶格描述Ni(100)，含顶位、空位、两类桥位，用First Reaction Method解主方程，以CH₃OD匹配Johnson与Madix实验TPD（137 K吸附、27 K s^?1升温）；用中心复合设计（Central Composite Design, CCD）在300–400 K、1–5 μPa、60–120 s内以二次多项式拟合产氢率，软件为Design-Expert 12；以脱附H₂/HD/D₂中H/D原子数占吸附CH₃OD四个H/D原子比例定义产氢率；用分子-only、解离-only、双路径三类网络比对实验谱，并以RMSE、AIC、BIC与嵌套F检验防过拟合，再用独立等温数据交叉验证，参数经1000次bootstrap与±20%局部敏感性分析。

3.1 动力学模型验证与机理：研究人员系统测试三类候选机制，发现分子-only（R²≈0.51）与解离-only（R²≈0.60）均无法复现H₂、CO脱附峰形与峰位，双路径模型使H₂与CO的R²达0.989/0.991、RMSE降至0.38/0.42、AIC/BIC显著下降，嵌套F检验p<0.001，说明双路径非多余参数而是数据本身所需。该网络含CH₃OD分子通道与CH₃OD(g)→CH₃O+D解离通道，后续顺序脱氢CH₃O→CH₂O+H→CHO+H→CO+H及氘代支链CH₂OD→CHOD+H→COD+D→CO，终产物为H₂、D₂、HD、CO。控速步为甲氧C–H断裂（Step 3，E_a=82.0 kJ mol^?1），degree of rate control系数>0.7；H+H→H₂表观E_a=56.20 kJ mol^?1低于洁净面文献值，归因为CO等共吸附削弱Ni–H键的覆盖度效应。bootstrap给出拟合参数95% CI很窄、相关系数<0.65，±20%扰动产氢率最大变<1.5%，结论稳健。解离通道快速供D解释早期D₂/HD，分子通道储CH₂OD于高温分解解释宽H₂峰，两路竞争统一了文献矛盾。

3.2 响应面优化：CCD下二次模型极显著（F=194.7，p<0.0001，调整R²=0.9873，适配精度43.7），失拟不显著。编码方程R=92.54+2.24A+0.731B+0.265C?0.065AB?0.366B²?0.178BC表明温度正向最强、压力呈抛物线、时间过长与高压交互略降益。三维曲面显示约379.7 K、3.50 μPa、107.8 s最优；RSM预测96.8%（95% CI 94.2–99.4%），DMC实测95.5%落在区间内。Pareto分析温度占效应约86%、压力9%、时间1.2%。

3.3 对催化剂设计的意义：双路径结论可用于合金化引入亲氧金属强化O–H断裂以提升CH₃O*覆盖，或抑制分子通道以减少HD/D₂、提H₂纯度；针对甲氧C–H断裂用应变、配体、缺陷降垒可在更低温获更高产；验证网络可作双金属/掺杂镍筛选基准。

讨论部分中，与同面乙醇分解比，甲醇无C–C断裂副产物，380 K即得95.5%产氢，优于乙醇约450 K下82%，证明甲醇是更优低温氢载体。局限包括：1–5 μPa属超高真空，工业压力需并入平均场微观动力学；未含积碳/CO歧化，低温短时无碳信号；理想Ni(100)不含台阶扭结；横向作用仅近邻占位阻断层级。研究人员指出该DMC–RSM框架在进高压实验前可低成本标定本征动力窗，并计划扩至Ni(211)、加压力依赖侧作用、建分层微观动力学模型。

结论部分翻译：本工作通过证明甲醇在Ni(100)上经分子与解离双路径吸附，解决了表面化学中长期争议；单靠实验或计算此前均不能确立此点。研究人员用经TPD（R²=0.9897）与独立等温数据验证的DMC–RSM框架表明，仅双路径模型能复现全部脱附曲线，AIC/BIC证其非过拟合。机理上解离通道速供D成早期D₂/HD，分子通道储CH₂OD高温分解致宽H₂峰。二次RSM以温度为首要因子，于379.7 K、3.50 μPa、107.8 s获95.5%产氢；±20%参数扰动产氢率变<1.5%。该网络把甲氧C–H断裂（degree of rate control=0.72）定为控速步，指引合金/掺杂选择性促解离。研究推进了过渡金属醇反应认知，并为分层微观动力学建模提供可靠本征动力学网络，弥合表面科学与催化工程间鸿沟。论文发表于《RSC Advances》。