Implementing the Blowers-Masel Approximation to Scale Activation Energy Based on Reaction Enthalpy in Mean-Field Microkinetic Modeling for Catalytic Methane Partial Oxidation

Mean-field microkinetic modeling is a powerful tool for catalyst design and the simulation of catalytic processes. The reaction enthalpies in a microkinetic model often need to be adjusted when changing species’ binding energies to model different catalysts, when performing thermodynamic sensitivity analyses, and when fitting experimental data. When altering reaction enthalpies, the activation energies should also be reasonably altered to ensure realistic reaction rates. The Blowers-Masel approximation (BMA) relates the reaction barrier to the reaction enthalpy. Unlike the Brønsted-Evans-Polani relationship, the BMA requires less data because only one parameter, the intrinsic activation energy, needs to be determined. We validate this application of BMA relations to model surface reactions by comparing against density functional theory data taken from the literature. By incorporating the BMA rate description into the open-source Cantera software, we enable a new workflow, demonstrated herein, allowing rapid screening of catalysts using linear scaling relationships and BMA kinetics within the process simulation software. For demonstration purposes, a catalyst screening for catalytic methane partial oxidation on 81 hypothetical metals is conducted. We compared the results with and without BMA-corrected rates. The heat maps of various descriptors (e.g., CH₄ conversion, syngas yield) show that using BMA rates instead of Arrhenius rates (with constant activation energies) changes which metals are most active. Heat maps of sensitivity analyses can help identify which reactions or species are the most influential in shaping the descriptor map patterns. Our findings indicate that while using BMA-adjusted rates did not markedly affect the most sensitive reactions, it did change the most influential species.