RMG Database for Chemical Property Prediction.

The Reaction Mechanism Generator (RMG) database for chemical property prediction is presented. The RMG database consists of curated datasets and estimators for accurately predicting the parameters necessary for constructing a wide variety of chemical kinetic mechanisms. These datasets and estimators are mostly published and enable prediction of thermodynamics, kinetics, solvation effects, and transport properties. For thermochemistry prediction, the RMG database contains 45 libraries of thermochemical parameters with a combination of 4564 entries and a group additivity scheme with 9 types of corrections including radical, polycyclic, and surface absorption corrections with 1580 total curated groups and parameters for a graph convolutional neural network trained using transfer learning from a set of >130 000 DFT calculations to 10 000 high-quality values. Correction schemes for solvent-solute effects, important for thermochemistry in the liquid phase, are available. They include tabulated values for 195 pure solvents and 152 common solutes and a group additivity scheme for predicting the properties of arbitrary solutes. For kinetics estimation, the database contains 92 libraries of kinetic parameters containing a combined 21 000 reactions and contains rate rule schemes for 87 reaction classes trained on 8655 curated training reactions. Additional libraries and estimators are available for transport properties. All of this information is easily accessible through the graphical user interface at https://rmg.mit.edu. Bulk or on-the-fly use can be facilitated by interfacing directly with the RMG Python package which can be installed from Anaconda. The RMG database provides kineticists with easy access to estimates of the many parameters they need to model and analyze kinetic systems. This helps to speed up and facilitate kinetic analysis by enabling easy hypothesis testing on pathways, by providing parameters for model construction, and by providing checks on kinetic parameters from other sources.

Reaction Mechanism Generator v3.0: Advances in Automatic Mechanism Generation.

In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.

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., CH4 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.

Quantifying the Impact of Parametric Uncertainty on Automatic Mechanism Generation for CO2 Hydrogenation on Ni(111).

Automatic mechanism generation is used to determine mechanisms for the CO2 hydrogenation on Ni(111) in a two-stage process while considering the correlated uncertainty in DFT-based energetic parameters systematically. In a coarse stage, all the possible chemistry is explored with gas-phase products down to the ppb level, while a refined stage discovers the core methanation submechanism. Five thousand unique mechanisms were generated, which contain minor perturbations in all parameters. Global uncertainty assessment, global sensitivity analysis, and degree of rate control analysis are performed to study the effect of this parametric uncertainty on the microkinetic model predictions. Comparison of the model predictions with experimental data on a Ni/SiO2 catalyst find a feasible set of microkinetic mechanisms within the correlated uncertainty space that are in quantitative agreement with the measured data, without relying on explicit parameter optimization. Global uncertainty and sensitivity analyses provide tools to determine the pathways and key factors that control the methanation activity within the parameter space. Together, these methods reveal that the degree of rate control approach can be misleading if parametric uncertainty is not considered. The procedure of considering uncertainties in the automated mechanism generation is not unique to CO2 methanation and can be easily extended to other challenging heterogeneously catalyzed reactions.

Engineering Bacillus megaterium Strains To Secrete Cellulases for Synergistic Cellulose Degradation in a Microbial Community.

Recent environmental concerns have intensified the need to develop systems to degrade waste biomass for use as an inexpensive carbon source for microbial chemical production. Current approaches to biomass utilization rely on pretreatment processes that include expensive enzymatic purification steps for the requisite cellulases. We aimed to engineer a synthetic microbial community to synergistically degrade cellulose by compartmentalizing the system with multiple specialized Bacillus megaterium strains. EGI1, an endoglucanase, and Cel9AT, a multimodular cellulase, were targeted for secretion from B. megaterium. A small library of signal peptides (SPs) with five amino acid linkers was selected to tag each cellulase for secretion from B. megaterium. Cellulase activity against amorphous cellulose was confirmed through a series of bioassays, and the most active SP constructs were identified as EGI1 with the LipA SP and Cel9AT with the YngK SP. The activity of the optimized cellulase secretion strains was characterized individually and in tandem to assess synergistic cellulolytic activity. The combination of EGI1 and Cel9AT yielded higher activity than either single cellulase. A coculture of EGI1 and Cel9AT secreting B. megaterium strains demonstrated synergistic behavior with higher activity than either monoculture. This cellulose degradation module can be further integrated with bioproduct synthesis modules to build complex systems for the production of high value molecules.