No abstract

Surface reaction networks involving hydrocarbons exhibit enormous complexity with thousands of species and reactions for all but the very simplest of chemistries. We present a framework for optimization under uncertainty for heterogeneous catalysis reaction networks using surrogate models that are trained on the fly. The surrogate model is constructed by teaching a Gaussian process adsorption energies based on group additivity fingerprints, combined with transition-state scaling relations and a simple classifier for determining the rate-limiting step. The surrogate model is iteratively used to predict the most important reaction step to be calculated explicitly with computationally demanding electronic structure theory. Applying these methods to the reaction of syngas on rhodium(111), we identify the most likely reaction mechanism. Propagating uncertainty throughout this process yields the likelihood that the final mechanism is complete given measurements on only a subset of the entire network and uncertainty in the underlying density functional theory calculations.

Over the last century, the industrialization of agriculture and the consumption of fossil fuels have resulted in a significant imbalance and redistribution in nitrogen-containing resources. This has sparked an interest in developing more sustainable and resilient approaches for producing nitrogencontaining commodities such as fertilizers and fuels. One largely neglected but emerging approach is photocatalytic nitrogen fixation. There is significant evidence that this process occurs spontaneously in terrestrial settings, and it has been demonstrated in numerous engineered systems. Yet many questions still remain unanswered regarding the rates, mechanisms, and impacts of photocatalytically producing fixed nitrogen "out of thin air". This work reviews the fascinating history of the reaction and examines current progress toward understanding and improving photofixation of nitrogen. This is supplemented by a quantitative review of the thermodynamic considerations and limitations for various reaction mechanisms. Finally, future prospects and preliminary performance targets for photocatalytic nitrogen fixation are discussed.

We introduce a general method for estimating the uncertainty in calculated materials properties based on density functional theory calculations. We illustrate the approach for a calculation of the catalytic rate of ammonia synthesis over a range of transition-metal catalysts. The correlation between errors in density functional theory calculations is shown to play an important role in reducing the predicted error on calculated rates. Uncertainties depend strongly on reaction conditions and catalyst material, and the relative rates between different catalysts are considerably better described than the absolute rates. We introduce an approach for incorporating uncertainty when searching for improved catalysts by evaluating the probability that a given catalyst is better than a known standard.

No abstract

Descriptor-based analysis is a powerful tool for understanding the trends across various catalysts. In general, the rate of a reaction over a given catalyst is a function of many parameters-reaction energies, activation barriers, thermodynamic conditions, etc. The high dimensionality of this problem makes it very difficult and expensive to solve completely, and even a full solution would not give much insight into the rational design of new catalysts. The descriptor-based approach seeks to determine a few ''descriptors'' upon which the other parameters are dependent. By doing this it is possible to reduce the dimensionality of the problem-preferably to 1 or 2 descriptors-thus greatly reducing computational efforts and simultaneously increasing the understanding of trends in catalysis. The ''CatMAP'' Python module seeks to standardize and automate many of the mathematical routines necessary to move from ''descriptor space'' to reaction rates for heterogeneous (electro) catalysts. The module is designed to be both flexible and powerful, and is available for free online. A ''reaction model'' can be fully defined by a configuration file, thus no new programming is necessary to change the complexity or assumptions of a model. Furthermore, various steps in the process of moving from descriptors to reaction rates have been abstracted into separate Python classes, making it easy to change the methods used or add new functionality. This work discusses the structure of the code and presents the underlying algorithms and mathematical expressions both generally and via an example for the CO oxidation reaction. Graphical Abstract

Synthesis gas (CO + H2) conversion is a promising route to converting coal, natural gas, or biomass into synthetic liquid fuels. Rhodium has long been studied as it is the only elemental catalyst that has demonstrated selectivity to ethanol and other C2+ oxygenates. However, the fundamentals of syngas conversion over rhodium are still debated. In this work a microkinetic model is developed for conversion of CO and H2 into methane, ethanol, and acetaldehyde on the Rh (211) and (111) surfaces, chosen to describe steps and close-packed facets on catalyst particles. The model is based on DFT calculations using the BEEF-vdW functional. The mean-field kinetic model includes lateral adsorbate-adsorbate interactions, and the BEEF-vdW error estimation ensemble is used to propagate error from the DFT calculations to the predicted rates. The model shows the Rh(211) surface to be ∼6 orders of magnitude more active than the Rh(111) surface, but highly selective toward methane, while the Rh(111) surface is intrinsically selective toward acetaldehyde. A variety of Rh/SiO2 catalysts are synthesized, tested for catalytic oxygenate production, and characterized using TEM. The experimental results indicate that the Rh(111) surface is intrinsically selective toward acetaldehyde, and a strong inverse correlation between catalytic activity and oxygenate selectivity is observed. Furthermore, iron impurities are shown to play a key role in modulating the selectivity of Rh/SiO2 catalysts toward ethanol. The experimental observations are consistent with the structure-sensitivity predicted from theory. This work provides an improved atomic-scale understanding and new insight into the mechanism, active site, and intrinsic selectivity of syngas conversion over rhodium catalysts and may also guide rational design of alloy catalysts made from more abundant elements.

No abstract