Systematic Enumeration of Elementary Reaction Steps in Surface Catalysis

The direct synthesis of complex chemicals from simple precursors (such as syngas) is one of the main objectives of current research in heterogeneous catalysis. To rationally design catalytic materials for this purpose, it is essential to identify the critical elementary reaction steps that ultimately determine a catalyst’s activity and selectivity with respect to a desired product. Unfortunately, the number of potentially relevant elementary steps is in the thousands, even for relatively simple target species like ethanol. The challenge of identifying the critical steps is thus akin to finding the proverbial needle in a haystack. Recently, a model-reduction scheme has been proposed, which tackles this problem by prescreening the barriers of all potential reactions with computationally inexpensive approximations. Although this route appears highly promising, it raises the question of how the starting point of the model-reduction process can be determined. In this contribution, we present a systematic method for enumerating all intermediates and elementary reactions relevant to a chemical process of interest. Using this approach, we construct reaction networks for C,H,O-containing systems consisting of up to four non-hydrogen atoms (more than 1 million reactions). Importantly, the scheme goes beyond simple bond-breaking reactions and allows considering rearrangement and transfer reactions as well. The presented reaction networks thus cover the chemistry of syngas-based processes (and beyond) to an unprecedented scale.

[1]  Andrew J. Medford,et al.  Analyzing the Case for Bifunctional Catalysis. , 2016, Angewandte Chemie.

[2]  B. Grzybowski,et al.  The core and most useful molecules in organic chemistry. , 2006, Angewandte Chemie.

[3]  Ture R. Munter,et al.  Scaling properties of adsorption energies for hydrogen-containing molecules on transition-metal surfaces. , 2007, Physical review letters.

[4]  J. Reymond,et al.  Exploring chemical space for drug discovery using the chemical universe database. , 2012, ACS chemical neuroscience.

[5]  Kevin Van Geem,et al.  Automatic reaction network generation using RMG for steam cracking of n‐hexane , 2006 .

[6]  Sebastian Matera,et al.  Adlayer inhomogeneity without lateral interactions: rationalizing correlation effects in CO oxidation at RuO2(110) with first-principles kinetic Monte Carlo. , 2011, The Journal of chemical physics.

[7]  C. Campbell The Degree of Rate Control: A Powerful Tool for Catalysis Research , 2017 .

[8]  Satoshi Maeda,et al.  Systematic exploration of the mechanism of chemical reactions: the global reaction route mapping (GRRM) strategy using the ADDF and AFIR methods. , 2013, Physical chemistry chemical physics : PCCP.

[9]  Stefan Grimme,et al.  Towards first principles calculation of electron impact mass spectra of molecules. , 2013, Angewandte Chemie.

[10]  Zachary W. Ulissi,et al.  To address surface reaction network complexity using scaling relations machine learning and DFT calculations , 2017, Nature Communications.

[11]  Pavlo O. Dral,et al.  Quantum chemistry structures and properties of 134 kilo molecules , 2014, Scientific Data.

[12]  Shawn M. Kathmann,et al.  Ethanol synthesis from syngas over Rh-based/SiO2 catalysts: A combined experimental and theoretical modeling study , 2010 .

[13]  J. Nørskov,et al.  Why gold is the noblest of all the metals , 1995, Nature.

[14]  Thomas Bligaard,et al.  The Brønsted–Evans–Polanyi relation and the volcano curve in heterogeneous catalysis , 2004 .

[15]  Dionisios G. Vlachos,et al.  Microkinetic Modeling for Water-Promoted CO Oxidation, Water−Gas Shift, and Preferential Oxidation of CO on Pt , 2004 .

[16]  Thomas Bligaard,et al.  Assessing the reliability of calculated catalytic ammonia synthesis rates , 2014, Science.

[17]  Woo Youn Kim,et al.  Efficient prediction of reaction paths through molecular graph and reaction network analysis† †Electronic supplementary information (ESI) available: Detailed information on reaction networks and pathways for two example reactions, Cartesian coordinates of molecules in the reaction networks obtained , 2017, Chemical science.

[18]  Brian A. Rohr,et al.  Computational Design of Active Site Structures with Improved Transition-State Scaling for Ammonia Synthesis , 2018 .

[19]  Sebastian Matera,et al.  A practical approach to the sensitivity analysis for kinetic Monte Carlo simulation of heterogeneous catalysis. , 2016, The Journal of chemical physics.

[20]  J. Nørskov,et al.  Towards the computational design of solid catalysts. , 2009, Nature chemistry.

[21]  R. West,et al.  Automatic Generation of Microkinetic Mechanisms for Heterogeneous Catalysis , 2017 .

[22]  Linda J. Broadbelt,et al.  Computer Generated Pyrolysis Modeling: On-the-Fly Generation of Species, Reactions, and Rates , 1994 .

[23]  H. Eyring The Activated Complex in Chemical Reactions , 1935 .

[24]  J. Nørskov,et al.  Using scaling relations to understand trends in the catalytic activity of transition metals , 2008, Journal of physics. Condensed matter : an Institute of Physics journal.

[25]  Eric Vanden-Eijnden,et al.  Simplified and improved string method for computing the minimum energy paths in barrier-crossing events. , 2007, The Journal of chemical physics.

[26]  Karsten Reuter,et al.  Perspective: On the active site model in computational catalyst screening. , 2017, The Journal of chemical physics.

[27]  Thomas Bligaard,et al.  Density functional theory in surface chemistry and catalysis , 2011, Proceedings of the National Academy of Sciences.

[28]  J. Nørskov,et al.  Universal Brønsted-Evans-Polanyi Relations for C–C, C–O, C–N, N–O, N–N, and O–O Dissociation Reactions , 2011 .

[29]  Markus Reiher,et al.  Context-Driven Exploration of Complex Chemical Reaction Networks. , 2017, Journal of chemical theory and computation.

[30]  Sebastian Matera,et al.  Addressing global uncertainty and sensitivity in first-principles based microkinetic models by an adaptive sparse grid approach. , 2018, The Journal of chemical physics.

[31]  Markus Reiher,et al.  Comprehensive Analysis of the Neglect of Diatomic Differential Overlap Approximation. , 2018, Journal of chemical theory and computation.

[32]  Ping Liu Water-gas shift reaction on oxide∕Cu(111): Rational catalyst screening from density functional theory. , 2010, The Journal of chemical physics.

[33]  B. Grzybowski,et al.  The 'wired' universe of organic chemistry. , 2009, Nature chemistry.

[34]  K. Reuter,et al.  Structure sensitivity in oxide catalysis: First-principles kinetic Monte Carlo simulations for CO oxidation at RuO2(111). , 2015, The Journal of chemical physics.

[35]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[36]  Jonathan E. Sutton,et al.  Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO2 , 2018 .

[37]  Thomas Bligaard,et al.  Discovery of technical methanation catalysts based on computational screening , 2007 .

[38]  B S Clausen,et al.  Catalyst design by interpolation in the periodic table: bimetallic ammonia synthesis catalysts. , 2001, Journal of the American Chemical Society.

[39]  Manos Mavrikakis,et al.  On the mechanism of low-temperature water gas shift reaction on copper. , 2008, Journal of the American Chemical Society.

[40]  Kevin Van Geem,et al.  Comprehensive reaction mechanism for n-butanol pyrolysis and combustion , 2011 .

[41]  Steven T. Evans,et al.  Mechanism of the Water Gas Shift Reaction on Pt: First Principles, Experiments, and Microkinetic Modeling , 2008 .

[42]  Michail Stamatakis,et al.  Unraveling the Complexity of Catalytic Reactions via Kinetic Monte Carlo Simulation: Current Status and Frontiers , 2012 .

[43]  Clausen,et al.  Design of a surface alloy catalyst for steam reforming , 1998, Science.

[44]  David Weininger,et al.  SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules , 1988, J. Chem. Inf. Comput. Sci..

[45]  Jean-Louis Reymond,et al.  Enumeration of 166 Billion Organic Small Molecules in the Chemical Universe Database GDB-17 , 2012, J. Chem. Inf. Model..

[46]  Young K. Park,et al.  Construction and optimization of complex surface‐reaction mechanisms , 2000 .

[47]  William H. Green,et al.  Reaction Mechanism Generator: Automatic construction of chemical kinetic mechanisms , 2016, Comput. Phys. Commun..

[48]  J. Nørskov,et al.  Ammonia synthesis and decomposition on a Ru-based catalyst modeled by first-principles , 2009 .

[49]  Johannes T. Margraf,et al.  Automatic generation of reaction energy databases from highly accurate atomization energy benchmark sets. , 2017, Physical chemistry chemical physics : PCCP.

[50]  Kurt Stokbro,et al.  Improved initial guess for minimum energy path calculations. , 2014, The Journal of chemical physics.

[51]  M. Fiałkowski,et al.  Architecture and evolution of organic chemistry. , 2005, Angewandte Chemie.

[52]  R. McGibbon,et al.  Discovering chemistry with an ab initio nanoreactor , 2014, Nature chemistry.

[53]  J. Nørskov,et al.  Ammonia Synthesis from First-Principles Calculations , 2005, Science.

[54]  Quantum chemistry: Quadruply bonded carbon. , 2012, Nature chemistry.

[55]  G. Frenking,et al.  Critical comments on "One molecule, two atoms, three views, four bonds?". , 2013, Angewandte Chemie.

[56]  Sebastian Matera,et al.  Examination of the concept of degree of rate control by first-principles kinetic Monte Carlo simulations , 2009 .

[57]  Ping Liu,et al.  Mechanism of ethanol synthesis from syngas on Rh(111). , 2009, Journal of the American Chemical Society.

[58]  K. Reuter Ab Initio Thermodynamics and First-Principles Microkinetics for Surface Catalysis , 2016, Catalysis Letters.