Automatic Generation of Microkinetic Mechanisms for

Yongping Fu , Morgan T. Rea , Jie Chen , Darien J. Morrow , Matthew P. Hautzinger , Yuzhou Zhao , Dongxu Pan , Lydia H. Manger , John C. Wright , Rand...
0 downloads 0 Views 2MB Size
Subscriber access provided by University of Colorado Boulder

Article

Automatic Generation of Microkinetic Mechanisms for Heterogeneous Catalysis C. Franklin Goldsmith, and Richard Henry West J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b02133 • Publication Date (Web): 06 Apr 2017 Downloaded from http://pubs.acs.org on April 11, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Automatic Generation of Microkinetic Mechanisms for Heterogeneous Catalysis C. Franklin Goldsmith∗,† and Richard H. West∗,‡ †School of Engineering, Brown University ‡Department of Chemical Engineering, Northeastern University E‐mail:  [email protected][email protected] 

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract A novel approach is presented for generating microkinetic mechanisms in heterogeneous catalysis. The open-source software RMG-Cat automatically develops a detailed list of elementary surface reactions, including thermodynamic properties for the adsorbates and parameterized rate coefficients for the reactions. The software proposes numerous possible surface intermediates and reactions, but it only retains those species that have a sufficiently high rate of formation. RMG-Cat was tested on the dry reforming of methane on nickel. The software correctly found the same set of elementary reactions as in a previously compiled microkinetic mechanism, as well as a few missing reactions. These results demonstrate the potential of this approach for predicting the dominant pathways in heterogeneous catalysis.

2 ACS Paragon Plus Environment

Page 2 of 43

Page 3 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction A major goal in catalysis is to be able to predict the performance of novel materials. Accurate prediction of the reactivity and selectivity of catalytic materials requires a detailed list of elementary surface reactions, or a microkinetic mechanism. 1 The microkinetic mechanism is the critical bridge between the electronic potential energy surface and the catalyst’s performance under industrially relevant conditions. Unlike more traditional models used in kinetics (e.g. pseudo-steady state, quasi-equilibration, irreversibility), microkinetic models do not make a priori assumptions to simplify the chemistry. 2 The advent of accurate, reliable density functional theory (DFT) methods for computing the electronic energy for periodic systems has given scientists the ability to rationalize catalyst performance at the atomic scale. These developments have produced the significant insight, led by Dr. Jens Nørskov and others, that a high-dimensional potential energy surface can be reduced to a lower-dimensional surface. 3–10 These performance descriptors are often referred to as “linear scaling relations”, since, according to the models, the binding energy of an adsorbate scales linearly with the binding energy of the adatom through which it binds, and the activation energy scales linearly with the change in the binding energy upon reaction. These scaling relations not only reduce the time required to generate a microkinetic mechanism, they also explain key trends in catalytic performance for many industrial processes over a wide array of catalysts. 8,9,11–32 Although linear scaling relations have led to major insights into catalyst design, the process of going from these predictors to microkinetic models remains time consuming and error-prone. For a given reactant and material, the scientist must complete a minimum of four tasks: (i) assume/guess which crystalline facet(s) are significant; (ii) assume/guess which elementary reactions are important on that facet, such as adsorption, bond activation, and desorption; (iii) parameterize the corresponding rate coefficients, e.g. an Arrhenius equation for each reaction; and (iv) assemble and solve the coupled differential equations that govern the chemical dynamics. 33–35 To ensure that the model is thermodynamically consistent, the 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

scientist must either compute the reverse rate coefficients in terms of a linear expansion of basis-set reactions, 36 or provide the temperature-dependent free energy for every species in the model. Scaling relations dramatically accelerate step (iii), but the process still requires significant human intervention. An alternative procedure would be to automate the mechanism generation process: the user completes the first task by providing a list of the reactants and the material facets, and the computer completes the remaining three tasks: (i) determine which reactions are important for each surface lattice, (ii) obtain accurate parameterization of the thermodynamic properties and rate coefficients, and (iii) solve the governing equations – without subsequent human intervention. Automated mechanism generation is consistent and thorough: the software does not miss pathways, and the resulting mechanism is thermodynamically consistent. Software that can automatically generate elementary kinetic mechanisms are revolutionizing gas-phase chemistry by predicting the combustion properties and complex interactions of novel biofuels, fuel blends, and additives. 37–53 These codes accomplish in hours what took experienced researchers months or years to accomplish by hand. Automated mechanism generation involves many problems at the intersection of chemistry and computer science. At a minimum, the kineticist must teach the computer how to: 1. Recognize when two or more species in the mechanism are equivalent. 2. Predict all the possible elementary reactions for each species and pair of species. 3. Determine which of the possible reactions are actually important. 4. Estimate accurately all the necessary thermodynamic and kinetic parameters. 5. Ensure that the mechanism is thermodynamically consistent. 6. Include flexibility for new reactants on novel materials. 4 ACS Paragon Plus Environment

Page 4 of 43

Page 5 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

7. Accomplish all of the above in a bug-free manner more quickly than a human. A challenge with automatic mechanism generation is determining which reactions are going to be kinetically significant. Although the final microkinetic mechanism may be sensitive to only a few pathways, these kinetic bottlenecks are seldom known a priori. As the size of the reactant increases, the number of possible elementary reactions grows exponentially. For example, moving from adsorbed methanol to hexanol, the number of possible reactions increases from 50 to 500,000. 34 The number of intermediates and reactions to be considered is vast (e.g. millions for gas-phase combustion of JP-10 Jet Fuel 48 ), so the parameter estimation process must be very fast. If all of these reactions are included, then the resulting mechanism quickly becomes unsolvable and therefore useless. Although one could avoid the combinatorial explosion by truncating the number of expansion steps or iterations, this approach runs the risks of including reactions that are kinetically irrelevant and missing reactions that are kinetically significant. Thus, the central challenge in automatic mechanism generation is more than just teaching the computer how to apply some rules to predict reactions, it must include teaching the computer how to determine which of these reactions is going to be kinetically significant. This selection process should not be biased by the user’s expectation, since the pathways that are dominant on one material may be irrelevant on another material. To address this challenge, the authors have taken existing software that achieves these goals for gas-phase mechanism generation – RMG 53 – and modified it to work for chemical kinetics on surfaces. This offshoot of RMG will be referred to as RMG-Cat. The version described here is archived and freely available online. 54,55

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Methods A complete description of RMG is provided by Gao et al.. The salient features that have been ported over to RMG-Cat are described below, but for a complete description of the underlying software, the reader is referred to Ref. 53 and references therein. An automatic mechanism generator like RMG-Cat requires four main components to function: a unique representation of chemical species; templates and rules to generate reactions; methods to estimate all the thermodynamic and kinetic parameters in the mechanism; and an algorithm to choose which reactions to include in the final mechanism.

Species Representation: Chemical Graph Theory A crucial step in mechanism generation is the representation of chemical species, so as to ensure that the same adsorbate is not represented by two or more distinct names in the mechanism. To accomplish this task, RMG-Cat represents adsorbates using chemical graph theory (ChemGraph), whereby the atoms are represented by nodes, and the bonds are represented by edges. Efficient graph theory algorithms are used for (sub)graph isomorphism detection throughout the code, e.g. to determine when two adsorbates are chemically equivalent, and to decompose molecules into functional groups for parameter estimation. Each adsorbate is represented by an adjacency list, a compact notation to describe the molecular structure in which each row corresponds to an atom in the molecule. This notation can preserve information regarding the bonding, bond order, charge, unpaired electrons, and lone pairs. Presently, the adjacency lists cannot preserve stereochemistry.

6 ACS Paragon Plus Environment

Page 6 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

at om ele n m u un en mb pa t er lo ired ne e l fo pair ec. rm al ch ar ge

Page 7 of 43

H H H C3

1 Ni u0 p0 c0 {2,S} 2 O u0 p2 c0 {1,S}{3,S} 3 C u0 p0 c0 {2,S}{4,S}{5,S}{6,S} 4 H u0 p0 c0 {3,S} 5 H u0 p0 c0 {3,S} 6 H u0 p0 c0 {3,S}

O2 Ni

Ni

1

bonds

Ni

Figure 1: sample adjacency list for methoxy adsorbed on Ni. There is a single (S) bond between atoms 1 and 2, and single bonds between atom pairs (2,3), (3,4), (3,5), and (3,6)

For RMG-Cat, four new bond-types were created to represent the metal-adsorbate bond: single (ıe.g. CH3 *), double (CH2 *), triple (CH*), and quadruple (C*). Figure 1 illustrates the current ChemGraph notation for methoxy adsorbed on Ni(111). The current implementation of chemgraph in RMG-Cat does not distinguish between different binding sites on a given 1

facet; the code assumes that the lowest-energy binding site for a given adsorbate is chosen, which is consistent with the mean-field approach of the kinetics. That said, there is nothing to prevent the code from including binding site as a possible additional column in the adjacency list. The authors are currently working on adding a fifth bond type for physisorbed or van der Waals species (e.g. CH4 *), as well as the possibility of combining bond types to form multi-dentate species (e.g. di-sigma ethylene, CH2 CH2 **).

Thermodynamic Properties Not only must the software be able to represent each species uniquely, but it also must be able to predict the thermodynamic properties of that species, such as binding energies, enthalpy, entropy, and temperature-dependent heat capacity. RMG-Cat accomplishes this task by combining a database of thermochemistry for known species, as well as the ability to predict the thermochemistry for new species that are not in the database. The current database 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 43

includes 21 adsorbates in Ni(111), listed in Table 1, which were taken from Blaylock et al. for methane steam reforming. For each species, RMG-Cat provides: the enthalpy of formation of the adsorbate, ∆f H ◦ (298 K); the standard-state entropy, S ◦ (298 K); and a polynomial expansion of the heat capacity at constant pressure as a function of temperature, Cp (T ). The thermodynamic properties were computed using the binding energies and normal modes, as provided in Ref. 56, using standard statistical thermodynamic methods. 57 Table 1: List of species in the current thermodynamics database. Adsorbate H* H2 * C* CH* CH2 * O* CH3 * OH* CH4 * H2 O* CO* COH* HCO* CH2 O* CHOH* OO* CH3 O* CH2 OH* CO2 * COOH* CH3 OH*

Metal Ni(111)

In a typical mechanism generation process, RMG-Cat will consider thousands of different adsorbates, even if only 10-100 intermediate species ultimately are retained in the final mechanism. Consequently, RMG-Cat must be able to generate accurate estimates of the thermodynamic properties rapidly, without first pausing to perform an expensive DFT calculation. RMG-Cat accomplishes this task efficiently by combining the ChemGraph functionality described above with Group Additivity, according to which the contribution of each functional group in a molecule contributes a fixed amount to the enthalpy, entropy, and heat capacity. 58–64 When RMG-Cat proposes a new species and confirms that it is not a duplicate of existing species, it first looks in the thermodynamic database to see if it can find an exact match. If it cannot, then it proceeds to estimate it according to the following algorithm: 1. For each species in Table 1, replace any H-atom with a wild card, then determine which species in the table is the closest match for the adsorbate in question. 2. Use a pre-compiled Adsorption Correction for that species. 3. Break the bond between the adsorbate and the metal to form a gas-phase precursor. 8 ACS Paragon Plus Environment

Page 9 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

4. Look up the thermodynamic properties of the gas-phase precursor in a database or estimate them using group additivity (considering additional complications, such as resonantly stabilized isomers, and whether or not the ground electronic state is a singlet or triplet, etc). 5. Determine the adsorbate thermochemistry by adding the pre-compiled adsorption correction to the estimated gas-phase properties.

Xadsorbate = Xgas-phase + ∆Xadsorption ,

X = ∆f H ◦ (298 K) , S ◦ (298 K) , Cp (T )

(1)

The process is summarized in Equation (1). For example, suppose RMG-Cat proposes adsorbed ethyl, CH3 CH2 *. The code first checks to see if this species is in Table 1. Since it is not, it recognizes that CH3 CH2 * is a close match to adsorbed methyl, CH3 *, by treating one of the H-atoms in methyl as a wild card (i.e. the ‘R’ in Figure 2). It then reads from the database the corresponding adsorption correction, or the difference in enthalpy, entropy, and heat capacity between the gas-phase radical CH3 • and the surface adsorbed CH3 *. Next, it removes CH3 CH2 * from the surface and estimates the thermochemistry of the gas-phase radical CH3 CH2 • using group additivity (if it cannot find it in a database). Finally, it adds the adsorption correction to the group additivity value to obtain estimates of the enthalpy, entropy, and heat capacity of CH3 CH2 *. This process is built upon the methods developed by Goldsmith for estimating partition functions for adsorbates. 64 Briefly, the approach assumes that the partition function of the adsorbate can be separated into the conserved modes and the transitional modes. 65–69 The conserved modes correspond to the 3N − 6 internal degrees of freedom of the gas-phase precursor (assuming it is nonlinear), and the transitional modes are the 6 degrees of freedom that describe the motion of the adsorbate relative to the surface.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

qadsorbate ≈

3N −6 Y

qiconserved

i

6 Y

qitransitional

Page 10 of 43

(2)

i

where qiconserved are the partition functions for the vibrational frequencies of the gas-phase precursor, and qitransitional are the partition functions for the relative motion – typically either a 2D lattice or 2D gas model. 56,57 The model assumes that, although the vibrational frequencies of the precursor do change upon adsorption, the contribution of these changes to the thermodynamic properties is negligible compared to the loss of the 3 external translational and 3 external rotational degrees of freedom and the corresponding gain of 6 new surface modes. Accordingly, the adsorption correction, ∆Xadsorption , in Equation (1) is calculated by

∆Hadsorption = H 6 new modes − 3H translation − 3H rotation + ∆Ebinding

(3a)

∆Sadsorption = S 6 new modes − 3S translation − 3S rotation

(3b)

∆Cp,adsorption = Cp6 new modes − 3Cptranslation − 3Cprotation

(3c)

The corrections for enthalpy and entropy are calculated at 298 K; the corrections for the heat capacity are calculated at the same temperatures as the group additivity values: 300, 400, 500, 600, 800, 1000, and 1500 K. This method has been applied successfully to model high-temperature partial oxidation of methane on platinum. 70 The gas-phase properties that were used for the pre-compiled adsorption corrections for the species in Table 1 were taken from a previously compiled database of DFT frequencies. 71

10 ACS Paragon Plus Environment

Page 11 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

R M

H

C

O

N

etc

M

M

M

M

M

C*R4

C*R3

C*R2

C*R

M

M

M

M

RRR R R C C M

M

C* O*R2 M

M

O*R

O*

M

M

R C

RR C

R C

H C

CR C

OR C

H O

CR O

OR O

M

M

M

M

M

M

M

M

M

Figure 2: Hierarchical tree for functional groups for adsorbate thermochemistry. R represents any atom; M is the metal surface site. For each branch in the second level down, the number of vertical lines between M and R is the bond order; the dashed vertical line to the right is for physisorbed or van der Waals complexes. The final row at the bottom indicates that each R can be further refined

The adsorption corrections summarized in Equation (3) are, in effect, a form of group additivity for adsorbate binding groups. These group additive values are arranged in a hierarchical tree-structured database, illustrated in Figure 2, with the most general groups at the top and more specific instances as child nodes. When using the database to estimate thermochemistry for an unknown adsorbate, the tree is climbed from the most specific node until a node is found that has a group value assigned. This way, the best possible estimate is used, but even groups with no training data can be estimated approximately. Once the tree structure is organized, the group additive values for a given node are fitted by linear least squares regression against all the data that falls under that node and its descendants in the tree. 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 43

This hierarchical tree database of group additive values has been successfully implemented in RMG for gas-phase thermochemistry, 53 reaction kinetics, 72 solvation parameters, 73 transport properties, 74 and most recently transition state geometries. 75

Kinetics Chemical reactions are classified according to “reaction families”. Each family contains: a template to identify when the reaction can occur, a recipe to specify the changes to be made to the molecular graphs, and some rules to determine the reaction rate coefficient. After the graphs are manipulated according to the recipe, they are checked for isomorphism with previously generated molecules to ensure the products are unique and the reactions are not duplicated. RMG-Cat specifies a reaction in a single direction. It assumes that all reactions are reversible, and computes the reverse rate coefficient directly from the forward rate and the equilibrium constant. Consequently, the RMG-Cat mechanisms are always thermodynamically consistent.

H Bond fission recipe: • Break Bond: {*2,S,*3} • Create Bond: {*3,S,*4} • Change Bond: {*1,D,*2}

H H C3 O2

O2

M1 + M4

M1 + M 4

H H H C3

Figure 3: sample reaction recipe for the bond fission of adsorbed methoxy to O* + CH3 * on a generic metal. The reaction breaks the single C–O bond, forms a new single (S) bond between the carbon atom and the metal, and it changes the nature of the oxygen-metal bond to double (D). The code is written in Python and is highly accessible and human readable.

12 ACS Paragon Plus Environment

Page 13 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

H-abstraction recipe: • Break Bond: {*2,S,*3} • Create Bond: {*3,S,*5} • Change Bond: {*1,D,*2} • Change Bond: {*4,S,*5}

H

H H C2 3 O5 M1 + M4

H

H C2

H O5 3

M1 + M4

Figure 4: sample reaction recipe for the H-atom transfer between adsorbed methyl and an oxygen adatom to methylene and OH on a generic metal.

RMG contains over 40 distinct families for gas-phase reactions. RMG-Cat now contains three additional reaction families specific to surface chemistry: adsorption/desorption, bond-fission (e.g. Figure 3), and H-abstraction (e.g. Figure 4). Each reaction family is further subdivided. For example, the Adsorption family includes rules for both dissociative and non-dissociated adsorption, with specific rules for the formation of adsorbates with single, double, triple and quadruple bonds. The code represents the rate coefficients for adsorption in terms of sticking coefficients, with built-in subroutines to convert from sticking coefficients to adsorption rate constants, and from adsorption rate constants to desorption rate constants. 76 Although direct reactions between gas-phase species and adsorbates are not presently implemented, Eley-Rideal kinetics will be incorporated soon, as will be physisorption and multi-dentate adsorption. Because the reaction families in the gas-phase RMG are all elementary reactions, it only considers unimolecular and bimolecular reactions; modifications were necessary to allow for dissociative adsorption, since these reactions contain three reactants (a gas-phase molecule and two vacant surface sites). After a reaction has been proposed by a reaction family, RMG-Cat searches a database of reactions to see if it can find an exact match. If the proposed reaction is not in the database, the Arrhenius parameters will be estimated according to a set of rules specific to that reaction family. The rules are arranged in hierarchical tree-structured databases, similar to the thermochemistry described above. The estimation rules, formulated as Brønsted-Evans-Polanyi relationships, are derived automatically from the database of known rate coefficients. 4,8,35 Kinetics and transition states predictions typically use two trees, one for each reactant in 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 43

bimolecular reactions, and possibly a third for the linking group for intramolecular reactions; the training procedure considers all combinations of each tree. At present, the kinetics database is sparsely populated. As the software develops, more specific values will be added, thereby improving the estimates for novel reactions, and increasingly specific rate rules will be added to the kinetics trees, including coverage dependence. The default rate rule for bond fission is:

k

bond fission

 (T ) = A

T Tref

n exp



  Ea0 + α∆Hrxn /RT

(4)

  A = 1.0 × 1017 m2 /mol-s n = 0.0 [−] Ea0 = 1.92 [eV] (185.3 kJ/mol) α = 0.84 [−]

where the pre-exponential factor was taken from Ref. 77, and the BEP value for the activation energy was taken from Ref. 8. The default rate rule for abstracting a Hydrogen atom from an adsorbate is:

k

H-abstraction

 (T ) = A

T Tref

n exp [Ea /RT ]

(5)

  A = 1.0 × 1017 m2 /mol-s n = 0.0 [−] Ea = 0.41 [eV] (40 kJ/mol)

where the pre-exponential factor and activation energy were fit from similar values in Ref. 77. RMG-Cat also includes the possibility that functional groups other than a terminal H

14 ACS Paragon Plus Environment

Page 15 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

can be abstracted, but that these rate coefficients should be orders of magnitude smaller. For non-H abstraction the default rate rule is:

k

R-abstraction

 (T ) = A

T Tref

n exp [Ea /RT ]

(6)

  A = 1.0 × 1015 m2 /mol-s n = 0.0 [−] Ea = 0.83 [eV] (80 kJ/mol)

The default sticking coefficients for dissociative and non-dissociative adsorption are:

S0diss. adsorp

 (T ) = A

T Tref

n exp [Ea /RT ]

(7)

A = 1.0 × 10−2 [−] n = 0.0 [−] Ea = 0.43 [eV] (41.8 kJ/mol)

S0nondiss. adsorp

 (T ) = A

T Tref

n exp [Ea /RT ]

(8)

A = 5.0 × 10−2 [−] n = 0.0 [−] Ea = 0.10 [eV] (10 kJ/mol)

The rate rules have been normalized on a per equivalent reaction site formulation, and RMG-Cat automatically incorporates the reaction path degeneracy into the forward rate

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

coefficient. For example, the generic pre-exponential factor for H-abstraction is 1 × 1017 m2 /mol-s. For the reaction C* + CH3 CH2 * → CH* + CH3 CH*, the pre-exponential factor is multiplied by 2, to account for the two equivalent H-atoms on C1 of the adsorbed ethyl, as noted in Table 2. Finally, the user can provide their own set of preferential rate coefficients, or a reaction library. RMG-Cat automatically will use those values first, in lieu of the internal database and estimation process. The user can force RMG-Cat to include all the reactions in the reaction library (a.k.a. a “seed mechanism”), or allow RMG-Cat to determine whether or not a given reaction in the reaction library contributes significantly towards the flux of a species and include only those reactions with a significant contribution. Importantly, the reactions in the user-defined reaction library need not conform to one of the pre-defined reaction families. An example of this functionality will be described in the Demonstration Section below. Our combined decades of experience using automatic mechanism generation for combustion chemistry has taught us that mechanism generation is an iterative process. We anticipate it being used as follows: RMG-Cat predicts a mechanism; sensitivity (or degree or rate control) analysis determines which reactions are critical; 78–89 the rate coefficients for these reactions are computed from first principles, and the new values are both added to the database and used to retrain the rate rules; then the re-trained RMG-Cat is used to predict an improved mechanism.

Mechanism Generation RMG-Cat considers all possible reactions that can be expressed within the reaction family system, without limiting the number of iterations. To avoid an exponential explosion in the number of (potentially irrelevant) species, it includes only those reactions that have a sufficiently high rate of reaction. 38 RMG-Cat divides all the species it considers into two groups: the “core” and the “edge”. The core contains species that are essential to the kinetics 16 ACS Paragon Plus Environment

Page 16 of 43

Page 17 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of the system; the edge contains species that have been proposed by a reaction, but have not been proven to be kinetically relevant. The process begins with the selection of reactants forming the initial core, as well as the temperature, pressure, and species concentrations. The algorithm follows the fastest pathways, according to the conditions at which the model will be used, and ignores the slow pathways that do not contribute to the overall process. Generate reactions of all core species and estimate parameters

Starting species

Reaction conditions (T, P, Xi)

Run simulations, tracking creation rate of edge species, Rn

Termination tolerance

Are any Rn > Rcritical ?

Add species with max Rn to model core

yes

no Generated mechanism Figure 3: Flow chart to illustrate the rate-based algorithm expansion. The user specifies the reactants (top-left

arrow), temperature, and concentrations (center-left arrow), and the tolerance (lower-left arrow). If Figure 5: the Flow chart pressure, to illustrate the rate-based algorithm for mechanism expansion. any new species is produced at an appreciable rate, the integration routine stops, the species is added to the core, and the process starts over.

Althoughrate-based this algorithmexpansion works exceptionally well for mechanism development, the species The iterative, algorithm is gas-phase described in Figure 5. For each

transition to heterogeneous catalysis may present new scientific challenges, especially for adsorbate-adsorbate interactions, coverage-dependent reaction pathways, and poisoning. In in the mechanism (both core and edge), the net rate of formation is computed: anticipation of these problems, we will track the surface coverage of all edge species, and include their effect on the model once it exceeds a critical threshold, even if their rate of generation is low. RMG-Cat will be able to run multiple reactor models, with the final mechanism as the union of the individual simulations. RMG-Cat will work seamlessly with the open-source code Cantera dNi equations (ODE).87 The initial reaction for solving the system of governing ordinary Ridifferential = (9) dt model used in the integration routine is a zero-dimensional isothermal, isobaric batch reactor, but if an alternative reactor model is desired, then it can be built first in Cantera and then imported into RMG-Cat. where Ni isThis themethod totalofnumber of isspecies i in a given phase (e.g. gas modelingof themoles chemistry inherently a mean-field approximation. The or surface). A advantage of mean-field kinetics is that the system of governing equations can be solved in a characteristic rate is efficient determined the species in the core: in this proposal are at computationally manner,from and many of the systems of interest temperatures sufficiently high that the mean-field approximation is appropriate.4, 88 Furthermore, the implicit Backward Differential Formula (BDF) that is used to solve the stiff ODE’s lends itself naturally to perform local sensitivity analysis.89 Local sensitivity analysis is a numerically efficient procedure that provides insight into the reactions and/or intermediates that control the s 89-97 X in heterogeneous catalysis, this method is commonly referred observable kinetic phenomena; Rdegree Ri2 98-100 species i ∈ core (10) to as Campbell’s char =of rate control. i

17 ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 43

The user provides an error tolerance  ∈ (0, 1]. For each species j in the edge, if the net rate of formation exceeds the characteristic rate scaled by the error tolerance, i.e. Rj ≥ Rchar , then that species is moved from the edge to the core. When a species is moved from the edge to the core, the simulation starts over. This process continues until a user-specified termination criterion is met and no new species are moved from the edge to the core. The current version of RMG-Cat has two possible termination criteria: a reaction time, or a percent conversion of a reactant. If both are specified, then RMG-Cat will terminate when one of them is met. RMG-Cat currently uses a zero-dimensional, isothermal, isochoric, batch reactor in the integration routine. The state variables are the total number of moles for each species, including gas-phase and surface species. The governing equations are formulated as a system of ordinary differential equations (ODE):

dNi = f (N, t; λ) dt

(11)

where N is the vector of species mole numbers, t is time, and λ is the vector of parameters, such as rate coefficients and thermodynamic properties. To complete the system of equations, the model assumes a reactor volume determined by one mole of gaseous reactants at the specified temperature and pressure (e.g. 2.48 × 10−2 m3 at STP). The user provides a ratio of catalytically active surface area to volume (m−1 ), and a surface site density for the catalyst (mole/m2 ). The system of ODE’s are integrated using PyDAS, 90 which is a Python wrapper for the DASPK integrator for stiff ordinary differential equation. 91

18 ACS Paragon Plus Environment

Page 19 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Other Properties RMG was recently extended to include more heteroatoms, transport properties, phase behavior, and solvation effects. 92–96 These tools have been carried over into RMG-Cat, thereby allowing RMG-Cat to be extended to liquid-phase catalysis in the future. Moreover, although gas-phase reactions are not expected to be significant for the trial application below, the core gas-phase chemistry of RMG is still present in RMG-Cat, so it is possible to build homogeneous/heterogeneously coupled systems, such as oxidative coupling of methane or catalytic combustion.

Results and discussion Proof of Concept: Methane dry reforming on nickel To demonstrate the capabilities of automatic mechanism generation, we apply RMG-Cat to the problem of methane dry reforming on nickel. Our goal is to see if the code can produce a microkinetic mechanism that is consistent with a mechanism that was developed “by hand” by experts – in this case, the mechanism developed by Olaf Deutschmann and coworkers. 77 The original mechanism in Delgado et al. consists of 26 reversible reactions, with the forward and reverse rate coefficient expressed explicitly. For this study three sets of modifications were made to the original mechanism. First, for each of the 26 reversible reactions, the rate coefficient was imported in only one direction; the reverse direction was computed directly from the thermochemistry, so as to maintain thermodynamic consistency. 97 The second modification was to neglect coverage effects in the activation energy. These effects are important, and they will be incorporated into RMG-Cat in a systematic fashion in the near future. For demonstration purposes, however, they are not essential. The third modification was for the adsorption of three species. In the Deutschmann mechanism, the adsorption of CH4 , CO2 , and H2 O proceed through pre-adsorbed complexes – CH4 *, CO2 *, and H2 O*,

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

respectively – which then dissociate on the surface. Because the van der Waals functionality has not yet been fully implemented, RMG-Cat cannot represent these three species, and so these three reactions were removed. Instead, any reaction that forms CH4 *, CO2 *, or H2 O* as a product was modified to produce the gas-phase species instead. For example, the reaction OH* + H* * ) H2 O* + * becomes OH* + H* * ) H2 O + 2*. The rate parameters in the associative desorption direction were included in the mechanism. Once the physisorption reactions for CH4 , CO2 , and H2 O were removed, 23 reactions remained. These reactions were then analyzed and used to populate the kinetics library. Three reactions were used to populate the Surface Abstraction tree: CH2 * + OH* → CH3 * + O*, CH* + OH* → CH2 * + O*, and C* + OH* → CH* + O*. Six reactions were used to populate the Surface Dissociation tree: CH2 * + H* → CH3 * + *, CH* + H* → CH2 * + *, C* + H* → CH* + *, OH* + * → H* + O*, HCO* + * → CO* + H*, and COOH* + * → CO* + OH*. The remaining 14 reversible elementary reactions were treated as a reaction library. These reactions include the necessary adsorption/desorption reactions, as well as elementary reactions that do not fit within the three reaction families described above (e.g. O-atom insertion reactions, such as HCO* + HO* * ) COOH* + H*). If RMG-Cat proposes a reaction that is in the Deutschmann mechanism, then RMGCat will use the specified kinetic parameters for that reaction (either by finding it within the kinetics database or within the reaction library); if, on the other hand, RMG-Cat proposes a reaction that is not in the Deutschmann mechanism, then it will estimate it from the rate libraries, as predicted by Equations (4)-(8). For the sample runs, RMG-Cat was initialized with a molar composition of 10% CH4 , 10% O2 , and 80% of N2 at 1000 K and 1 bar. Additionally, the nickel surface was completely vacant, with a surface site density of 2.9 × 10−9 mole/cm2 , and the reactor was assumed to have a surface-to-volume ratio of 105 m−1 . The termination criterion was set to a CH4 conversion of 50% or 0.1 seconds, whichever comes first. The error tolerance was decreased gradually from  = 1 × 100 to  = 1 × 10−8 .

20 ACS Paragon Plus Environment

Page 20 of 43

Page 21 of 43

no. of species

10 3

no. of reactions

10 4

(a)

no. of core reactions

400

(b)

(c)

350

Abstraction Dissociation Adsorption Deutcschmann

300 10 3 ge

250

10 2

ed

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

g ed

10 2

co

e cor

200

e

150 100

re

50

10 1

10 1 10

0

10

-4

10

-8

0 10

error tolerance, ε

0

10

-4

10

-8

error tolerance, ε

10 0

10 -4

10 -8

error tolerance, ε

Figure 6: Results for adjusting the error tolerance. (a) and (b) demonstrate how the number of species and reactions, respectively, in the core and edge grows exponentially with increasing the tolerance tightness. (c) the total number of instances from each reaction family in the core, as a function of error tolerance. 98

Figure 6 demonstrates how the size of the mechanism increases as the error tolerance is tightened. At  = 1, the core contained 17 species and 30 reactions, and the edge contained 42 species and 61 reactions. By  = 1 × 10−8 , the core contained 68 species and 366 reactions, and the edge 740 species and 1505 reactions. An important aspect of the graph theory methodology is its computational efficiency. The total mechanism generation process was ∼1 minute for  = 1 × 10−1 and ∼5 minutes for  = 1 × 10−8 , on a single core of a 3.2 GHz Intel Core i5 (circa 2013 iMac).

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 43

Table 2: List of reactions in the RMG-Cat mechanism for an initial composition of 10% CH4 , 10% O2 , and 80% N2 at 1000 K and 1 bar. Units of A are [m2 /mol-s] (or dimensionless - see note b) and Ea are [kJ/mol]. Rxn

reaction

A

n

Ea

notes

a

1

CH3 * + H* * ) CH4 + 2*

1.44×1018

−0.087

63.4

c,g,h

10−2

2

CH2 * + H* * ) CH3 * + *

3.09×1019

−0.087

57.2

e, h

10−2

3

CH* + H* * ) CH2 * + *

9.77×1020

−0.087

81.0

e, h

10−2

4

CH2 * + CH2 * * ) CH* + CH3 *

2.00×1017

0.0

40.0

f

10−2

b,c,h

10−2

e,h

10−2

5

H2 + 2* * ) 2H*

3.00×10−2

0.0

5.0

6

H* + C* * ) CH* + *

1.70×1020

−0.5

157.9

7

CH3 * + C* * ) CH2 * + CH*

3.00×1017

0.0

40.0

f

10−2

8

CH2 * + C* * ) CH* + CH*

2.00×1017

0.0

40.0

f

10−2

COOH* + * * ) CO2 + H* + *

3.73×1016

0.475

33.6

c,g,h

10−2

2CO* * ) CO2 + * + C*

1.62×1010

0.5

241.7

c,g,h

10−2

9 10 11

C* + OH* * ) CO* + H*

3.88×1021

0.188

62.5

c,h

10−2

12

COOH* + * * ) CO* + OH*

1.46×1020

−0.213

54.3

e,h

10−2

2.98×1018

0.101

25.8

c, h

10−2

1.75×109

0.0

116.2

c,h

10−2

123.6

c,g,h

10−2

13

CH3 * + OH* * ) CH4 + O* + *

14

CO* + * * ) C* + O*

15

O* + CO* * ) CO2 + *

2.15×1015

0.0

16

OH* + * * ) H* + O*

2.25×1016

0.188

29.6

e,h

10−2

17

CH2 * + OH* * ) CH3 * + O*

1.39×1017

0.101

19.0

f,h

10−2

18

CH* + OH* * ) CH2 + O*

4.40×1018

0.101

42.4

f,h

10−2

19

C* + OH* * ) CH* + O*

2.43×1017

−0.312

118.9

f,h

10−2

5.00×10−1

0.0

0.0

b,c,h

10−2

109.9

c,h

10−2

15.9

c,h

10−2

e,h

10−2

f

10−2

20

CO + * * ) CO*

21

CH* + O* * ) HCO* + *

4.59×1016

0.0

22

HCO* + OH* * ) COOH* + H*

2.28×1016

0.263

23

HCO* + * * ) CO* + H*

3.71×1017

0.0

0.0

1.00×1017

0.0

40.0

24

HCO* + CH2 * * ) CO* + CH3 *

Continued on Next Page. . .

22 ACS Paragon Plus Environment

Page 23 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 2 – Continued Rxn

reaction

A

n

Ea

notes

a

25

HCO* + CH* * ) CO* + CH2 *

1.00×1017

0.0

40.0

f

10−2

26

HCO* + C* * ) CO* + CH*

1.00×1017

0.0

40.0

f

10−2

27

HCO* + O* * ) CO* + OH*

1.00×1017

0.0

40.0

f

10−2

28

CH3 CH2 * + * * ) CH2 * + CH3 *

1.00×1017

0.0

117.7

e

10−2

29

C2 H6 + 2* * ) 2 CH3 *

2.00×10−2

0.0

41.8

b,d

10−2

30

C2 H6 + 2* * ) CH3 CH2 * + H*

1.20×10−1

0.0

41.8

b,d

10−2

31

CH3 CH2 * + * * ) CH3 CH* + H*

2.00×1017

0.0

180.2

e

10−3

32

CH3 CH* + * * ) CH3 * + CH*

1.00×1017

0.0

148.7

e

10−3

33

CH3 CH* + OH* * ) CH3 CH2 * + O*

1.00×1017

0.0

40.0

f

10−3

34

CH3 CH* + HCO* * ) CH3 CH2 * + CO*

1.00×1017

0.0

40.0

f

10−3

35

CH2 * + CH3 CH2 * * ) CH3 * + CH3 CH*

2.00×1017

0.0

40.0

f

10−3

36

CH2 * +CH3 CH* * ) CH* + CH3 CH2 *

2.00×1017

0.0

40.0

f

10−3

37

C* + CH3 CH2 * * ) CH* + CH3 CH*

2.00×1017

0.0

40.0

f

10−3

38

C* + CH3 CH2 * * ) CH2 * + CH3 C*

1.00×1015

0.0

80.0

f

10−3

39

CH3 C* + * * ) CH3 * + C*

1.00×1017

0.0

604.3

e

10−3

40

CH3 CH* + * * ) CH3 C* + H*

1.00×1017

0.0

0.0

e

10−3

41

O* +CH3 CH* * ) OH* + CH3 C*

1.00×1017

0.0

40.0

f

10−3

42

CO* +CH3 CH* * ) HCO* + CH3 C*

1.00×1017

0.0

40.0

f

10−3

43

CH2 * +CH3 CH* * ) CH3 * + CH3 C*

1.00×1017

0.0

40.0

f

10−3

44

CH* +CH3 CH* * ) CH2 * + CH3 C*

1.00×1017

0.0

40.0

f

10−3

45

CH* +CH3 C* * ) C* + CH3 CH*

1.00×1017

0.0

40.0

f

10−3

1.00×1017

0.0

40.0

f

10−3

1.00×10−2

0.0

16.0

b,d

10−3

46 47

CH3 CH* + CH3 CH* * ) CH3 CH2 * + CH3 C* CH3 • + * * ) CH3 *

48

CH3 C*O + * * ) CH3 * + CO*

1.00×1017

0.0

316.8

e

10−3

49

CO* + CH3 CH2 * * ) CH3 C*O + CH2 *

1.00×1015

0.0

80.0

f

10−3

50

CO* + CH3 CH* * ) CH3 C*O + CH*

1.00×1015

0.0

80.0

f

10−3

Continued on Next Page. . .

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 43

Table 2 – Continued Rxn

reaction

A

n

Ea

notes

a

51

CH3 C*O + C* * ) CO* + CH3 C*

1.00×1015

0.0

80.0

f

10−3

52

H* + OH* * ) H2 O + 2*

1.85×1016

0.086

41.5

c,g,h

10−3

2.34×1016

0.274

92.3

c,g,h

10−3

b,d

10−4

OH* + OH* * ) H2 O + O* + *

53 54

C3 H8 + 2* * ) CH3 * + CH3 CH2 *

4.00×10−2

0.0

71.1

55

CH3 CH2 CH2 * + * * ) CH2 * + CH3 CH2 *

1.00×1017

0.0

215.7

e

10−4

1.20×10−1

0.0

41.8

d

10−4

56b

C3 H8 + 2* * ) H* + CH3 CH2 CH2 *

57

CH3 CH*CH3 + * * ) CH3 * + CH3 CH*

2.00×1017

0.0

220.8

e

10−4

58

CH3 CH2 * + CH3 CH* * ) CH2 * + CH3 CH*CH3

1.00×1015

0.0

80.0

f

10−4

59

CH3 CH* + CH3 CH* * ) CH* + CH3 CH*CH3

1.00×1015

0.0

80.0

f

10−4

60

C* + CH3 CH*CH3 * ) CH3 C* + CH3 CH*

2.00×1015

0.0

80.0

f

10−4

61

C3 H8 + 2* * ) H* + CH3 CH*CH3

4.00×10−2

0.0

41.8

b,d

10−4

62

CO* + CH3 CH2 * * ) CH3 C*O + CH2 *

1.00×1015

0.0

80.0

f

10−4

63

CO* + CH3 CH* * ) CH3 C*O + CH*

1.00×1015

0.0

80.0

f

10−4

64

C* + CH3 C*O * ) CH3 C* + CO*

1.00×1015

0.0

80.0

f

10−4

2.00×1015

0.0

80.0

f

10−4

65

CO* + CH3 CH*CH3 * ) CH3 C*O + CH3 CH* +

a

minimum value of the error tolerance  for RMG-Cat to find this reaction.

b

rate coefficient is a sticking coefficient; pre-exponential factor A is unitless.

c

reaction directly from Deutschmann library

d

adsorption

e

dissociation

f

abstraction

g

associative desorption (see discussion in text)

h

kinetic parameters from Deutschmann [77]

i

kinetic parameters estimated by RMG-Cat

24 ACS Paragon Plus Environment

Page 25 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In order to illustrate the rate-based algorithm in practice, Table 2 lists the reactions found by RMG-Cat for three values of error tolerance:  = 1×10−2 , 1×10−3 , and 1×10−4 . For  ≥ 1 × 10−2 , RMG-Cat returned 30 elementary reactions. 20 of those 30 reactions were direct hits from the Deutschmann mechanism, including the main pathways for adsorption of CH4 and CO2 and desorption of H2 and CO. When the tolerance was decreased to  = 1 × 10−3 , RMG-Cat found 23 additional reactions, including the 2 Deutschmann pathways that lead to small concentrations of H2 O. By  = 1 × 10−4 , the mechanism contains 65 reactions. Tightening the tolerance beyond 1 × 10−4 further increases the size and complexity of the mechanism, including the formation of trace gas-phase products not listed in Table 2, such as CH2 O. The only reaction in the original Deutschmann mechanism that RMG-Cat did not include under these conditions was the associative desorption of 2O* to O2 . Inspection of the edge confirms that RMG-Cat did propose this reaction, but it determined that the flux was insufficient to move it from the edge to the core. Under these conditions, O2 is not predicted to be formed at a significant rate; integration of the original literature mechanism yielded a peak mole fraction of y2 ≤ 1 × 10−20 , so the rate-based algorithm behaved properly. It is important to emphasize that 43 of the 65 rate coefficients in Table 2 were created by RMG-Cat based upon the rate rules, and as such have higher uncertainty than the 22 rate coefficients that were taken from Delgado et al.. Nonetheless, including these new proposed reactions in a mechanism is a useful starting point for subsequent sensitivity analysis. In time, we expect the accuracy of these predictions to improve considerably, as the database becomes populated with more values. Additionally, as illustrated in Figure 6c, the largest single percentage of the proposed reactions are abstraction reactions. This result is in contrast to many mechanisms in the literature, which assume that bond-fission pathways dominate. The higher-than-expected prevalance of abstraction reactions is due in part, no doubt, to the simplicity of the rate rule; with a more densely populated database, we could develop more accurate BEP-style rules for this family, and these better rules may end up predicting fewer reactions to be signifi-

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

cant. However, it is also worth emphasizing that these reactions could play an increasingly important role under high-pressure conditions for which there are fewer vacant sites. As detailed above, the RMG-Cat mechanism is larger than the original mechanism of Delgado et al.. To illustrate the impact of the new reactions, a new simulation was performed. To simplify the comparison, the same initial conditions were chosen as in Figure 6, but the termination criterion was set to 1 second, and the constant T,V batch reactor was integrated to 1 second. These results are illustrated in Figures 7 - 9. 0.10

0.20

0.006

0.05

moles

0.15

moles

moles

CH 4

0.10 0.05

0.00 0.5

1.0

0.002

H2 O 0.000

0.0

time (s)

0.5

1.0

time (s)

0.10

0.004

CO

0.00 0.0

0.0

0.5

1.0

time (s)

0.20

CO 2

RMG-Cat Delgado et al.

0.15

moles

moles

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 43

0.05

0.10 0.05

0.00

H2

0.00 0.0

0.5

time (s)

1.0

0.0

0.5

1.0

time (s)

Figure 7: Simulation results for methane dry reforming at 1000 K and 1 bar (see text for initial conditions). A comparison of major gas-phase species. Black dashes are the results from the Deutschmann mechanism, and the red lines are from the RMG-Cat mechanism. 98 Figure 7 presents the gas-phase species that were in the Deutschmann mechanism. The agreement among the reacants, CH4 + CO2 , and major products, CO + H2 , is excellent. The profiles for H2 O differ slightly, with the original mechanism predicting a peak H2 O that is 13% higher than that of the new mechanism.

26 ACS Paragon Plus Environment

Page 27 of 43

10 -3 10 -4 10 -5

moles

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

10 -6

C2 H6 10

-7

C3 H8

CH 2 O

10 -8

CH 3

10 -9 0.0

0.5

1.0

time (s)

Figure 8: Simulation results for methane dry reforming at 1000 K and 1 bar. Some gas-phase species predicted by RMG-Cat that were not in the Deutschmann mechanism. 98 Additionally, RMG-Cat found several gas-phase species there were absent from the Deutschmann mechanism. The most abundant of these new products are shown in Figure 8, and even these species are mostly present in sub-ppm levels. The concentration of methyl radicals, CH3 , peaks around 1 ms. Effectively, the adsorption/desorption reaction equilibrates, and the radical concentration decreases as the surface CH3 * is consumed by surface oxidation reactions.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry

10 0

vacant

10 -1

H*

CH* 10

site fraction

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 43

CO*

-2

C* CH 3 C*O

10 -3

CH 3 C* 10 -4

O*

10 -5

C 2 HO*

CH 2 *

O*

CHO 2 *

10 -6 10 -6

10 -5

10 -4

10 -3

10 -2

10 -1

10 0

time (s)

Figure 9: Simulation results for methane dry reforming at 1000 K and 1 bar. Surface site fractions. Solid lines are the RMG-Cat predictions, black dashes are predicted by the Deutschmann mechanism. 98

Figure 9 compares the main adsorbate site fractions. As with Figure 7, the agreement is excellent, with minor deviations for C* and O* in the µs region, and new adsorbates – such as CH3 C* and CH3 C*O – becoming significant in the ms region.

Future Directions RMG-Cat is open source, and it is intended to be a community–driven project. Our immediate goals are to populate the thermodynamic and kinetic databases with more data (including metals beyond nickel) and to refine the rate rules within the reaction families (e.g. further distinguish the bond fission rule based upon both the bond type and nearest neighbors, or refine the abstraction family to distinguish between abstracting a hydrogen atom versus some other moeity). In principle, reactor models other than a constant-T,V batch reaction could be used 28 ACS Paragon Plus Environment

Page 29 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

within the integration routine, such as plug flow reactors (PFR). The main challenge with a PFR is that the governing equations are differential algebraic equations (DAE) and thus are somewhat more difficult to implement numerically. Furthermore, there is nothing to prevent RMG-Cat from moving beyond mean-field models and apply the rate-based screening to stochastic methods, such as kinetic Monte Carlo (kMC), which may give a physically more realistic portrayal of the heterogeneous chemistry. 99,100

Conclusions A novel software was presented that can automatically generate microkinetic mechanisms for heterogeneous catalysis. The software iteratively proposes new reactions according to reactions templates, but it includes only the reactions that are kinetically significant. It combines a precompiled database of thermodynamic properties and rate coefficients for known species and reactions with the ability to estimate these parameters for new species and reactions. The software was tested on methane dry reforming on nickel. It successfully “discovered” the many of the same reactions as a previously compiled mechanism that was developed over many years by experts.

Acknowledgement The project was made possible by collaborative tools GitHub, Slack, FaceTime, Screen Sharing, and Blue Bottle. CFG gratefully acknowledges support from Brown University, and RHW the same from Northeastern University.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Supporting Information Available The following files are available free of charge. • Instructions, input files, output files, and post-processing scripts to generate figures 6–9, and the figures themselves, are at http://doi.org/10.6084/m9.figshare.4640098.v1 (ref. 98). • The source code for the version of RMG-Cat described here is archived at https://doi.org/10.5281/zenodo.290119 (ref. 54). • The database for the version of RMG-Cat described here is archived at https://doi.org/10.5281/zenodo.290120 (ref. 55). This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Dumesic, J. A. The Microkinetics of heterogeneous catalysis; American Chemical Society: Washington, DC, 1993. (2) Chorkendorff, I.; Niemantsverdriet, J. W. Concepts in Modern Catalysis and Kinetics; Wiley: Weinheim, 2007. (3) Hammer, B.; Nørskov, J. K. Electronic factors determining the reactivity of metal surfaces. Surface Science 1995, 343, 211–220. (4) Bligaard, T.; Norskov, J. K.; Dahl, S.; Matthiesen, J.; Christensen, C. H.; Sehested, J. The Bronsted-Evans-Polanyi relation and the volcano curve in heterogeneous catalysis. Journal of Catalysis 2004, 224, 206–217. (5) Nørskov, J. K. Universality in Heterogeneous Catalysis. Journal of Catalysis 2002, 209, 275–278. 30 ACS Paragon Plus Environment

Page 30 of 43

Page 31 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(6) Abild-Pedersen, F.; Greeley, J.; Studt, F.; Rossmeisl, J.; Munter, T.; Moses, P.; SkXXlason, E.; Bligaard, T.; Nørskov, J. K. Scaling Properties of Adsorption Energies for Hydrogen-Containing Molecules on Transition-Metal Surfaces. Physical Review Letters 2007, 99, 4–7. (7) Fernanez, E. M.; Moses, P. G.; Toftelund, A.; Hansen, H. A.; Martinez, J. I.; AbildPedersen, F.; Kleis, J.; Hinnemann, B.; Rossmeisl, J.; Bligaard, T. et al. Scaling relationships for adsorption energies on transition metal oxide, sulfide, and nitride surfaces. Angew. Chemie-International Ed. 2008, 47, 4683–4686. (8) Wang, S.; Temel, B.; Shen, J.; Jones, G.; Grabow, L. C.; Studt, F.; Bligaard, T.; Abild-Pedersen, F.; Christensen, C. H.; Nørskov, J. K. Universal Brønsted-EvansPolanyi Relations for C-C, C-O, C-N, N-O, N-N, and O-O Dissociation Reactions. Catalysis Letters 2011, 141, 370–373. (9) Jones, G.; Studt, F.; Abild-Pedersen, F.; Nørskov, J. K.; Bligaard, T. Scaling relationships for adsorption energies of C2 hydrocarbons on transition metal surfaces. Chemical Engineering Science 2011, 66, 6318–6323. (10) Medford, A. J.; Vojvodic, A.; Hummelshøj, J. S.; Voss, J.; Abild-Pedersen, F.; Studt, F.; Bligaard, T.; Nilsson, A.; Nørskov, J. K. From the Sabatier principle to a predictive theory of transition-metal heterogeneous catalysis. Journal of Catalysis 2015, 328, 36 – 42. (11) Schumacher, N.; Boisen, A.; Dahl, S.; Gokhale, A. A.; Kandoi, S.; Grabow, L. C.; Dumesic, J. A.; Mavrikakis, M.; Chorkendorff, I. Trends in low-temperature water-gas shift reactivity on transition metals. Journal of Catalysis 2005, 229, 265–275. (12) Andersson, M.; Bligaard, T.; Kustov, A.; Larsen, K.; Greeley, J.; Johannessen, T.; Christensen, C. H.; Nørskov, J. K. Toward computational screening in heterogeneous

31 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

catalysis: Pareto-optimal methanation catalysts. Journal of Catalysis 2006, 239, 501– 506. (13) Grabow, L. C.; Xu, Y.; Mavrikakis, M. Lattice strain effects on CO oxidation on Pt(111). Physical chemistry chemical physics : PCCP 2006, 8, 3369–74. (14) Falsig, H.; Bligaard, T.; Rass-Hansen, J.; Kustov, A. L.; Christensen, C. H.; Nørskov, J. K. Trends in catalytic NO decomposition over transition metal surfaces. Topics in Catalysis 2007, 45, 117–120. (15) Greeley, J.; Rossmeisl, J.; Hellman, A.; Nørskov, J. K. Theoretical Trends in Particle Size Effects. Zeitschrift für Physikalische Chemie 2007, 221, 1209–1220. (16) Grabow, L. C.; Gokhale, A. A.; Evans, S. T.; Dumesic, J. A.; Mavrikakis, M. Mechanism of the Water Gas Shift Reaction on Pt: First Principles, Experiments, and Microkinetic Modeling. Journal of Physical Chemistry C 2008, 112, 4608–4617. (17) Jones, G.; Jakobsen, J. G.; Shim, S. S.; Kleis, J.; Andersson, M. P.; Rossmeisl, J.; Abild-Pedersen, F.; Bligaard, T.; Helveg, S.; Hinnemann, B. et al. First principles calculations and experimental insight into methane steam reforming over transition metal catalysts. Journal of Catalysis 2008, 259, 147–160. (18) Studt, F.; Abild-Pedersen, F.; Bligaard, T.; Sørensen, R. Z.; Christensen, C. H.; Nørskov, J. K. On the role of surface modifications of palladium catalysts in the selective hydrogenation of acetylene. Angewandte Chemie (International ed. in English) 2008, 47, 9299–302. (19) Ferrin, P.; Simonetti, D.; Kandoi, S.; Kunkes, E.; Dumesic, J. A.; Nørskov, J. K.; Mavrikakis, M. Modeling ethanol decomposition on transition metals: a combined application of scaling and Brønsted-Evans-Polanyi relations. Journal of the American Chemical Society 2009, 131, 5809–15.

32 ACS Paragon Plus Environment

Page 32 of 43

Page 33 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(20) Grabow, L. C.; Uhlrich, J. J.; Kuech, T. F.; Mavrikakis, M. Effectiveness of in situ CH4 annealing treatments for the removal of oxygen from GaN surfaces. Surface Science 2009, 603, 387–399. (21) Grabow, L. C.; Hvolbæk, B.; Nørskov, J. K. Understanding trends in catalytic activity: The effect of adsorbate-adsorbate interactions for CO oxidation over transition metals. Topics in Catalysis 2010, 53, 298–310. (22) Hirvi, J. T.; Kinnunen, T.-J. J.; Suvanto, M.; Pakkanen, T. a.; Nørskov, J. K. CO oxidation on PdO surfaces. The Journal of Chemical Physics 2010, 133, 084704– 084704. (23) Durand, W. J.; Peterson, A. A.; Studt, F.; Abild-Pedersen, F.; Nørskov, J. K. Structure effects on the energetics of the electrochemical reduction of CO2 by copper surfaces. Surface Science 2011, 605, 1354–1359. (24) Grabow, L. C.; Mavrikakis, M. Mechanism of methanol synthesis on Cu through CO2 and CO hydrogenation. ACS Catal. 2011, 1, 365–384. (25) Grabow, L. C.; Studt, F.; Abild-Pedersen, F.; Petzold, V.; Kleis, J.; Bligaard, T.; Nørskov, J. K. Descriptor-based analysis applied to HCN synthesis from NH3 and CH4 . Angewandte Chemie (International ed. in English) 2011, 50, 4601–5. (26) Peterson, A. A.; Dreher, M.; Wambach, J.; Nachtegaal, M.; Dahl, S.; Nørskov, J. K.; Vogel, F. Evidence of Scrambling over Ruthenium-based Catalysts in Supercriticalwater Gasification. Chemcatchem 2012, 4, 1185–1189. (27) Peterson, A. A.; Nørskov, J. K. Activity Descriptors for CO2 Electroreduction to Methane on Transition-Metal Catalysts. Journal of Physical Chemistry Letters 2012, 3, 251–258.

33 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(28) Peterson, A. A.; Nørskov, J. K. Activity Descriptors for CO2 Electroreduction to Methane on Transition-Metal Catalysts. J. Phys. Chem. Lett. 2012, 3, 251–258. (29) Lausche, A. C.; Medford, A. J.; Khan, T. S.; Xu, Y.; Bligaard, T.; Abild-Pedersen, F.; Nørskov, J. K.; Studt, F. On the effect of coverage-dependent adsorbate-adsorbate interactions for CO methanation on transition metal surfaces. Journal of Catalysis 2013, 307, 275–282. (30) Medford, A. J.; Lausche, A. C.; Abild-Pedersen, F.; Temel, B.; Schjødt, N. C.; Nørskov, J. K.; Studt, F. Activity and Selectivity Trends in Synthesis Gas Conversion to Higher Alcohols. Topics in Catalysis 2013, 57, 135–142. (31) Montoya, J. H.; Peterson, A. A.; Norskov, J. K. Insights into CC Coupling in CO2 Electroreduction on Copper Electrodes. Chemcatchem 2013, 5, 737–742. (32) Yoo, J. S.; Khan, T. S.; Abild-Pedersen, F.; Nørskov, J. K.; Studt, F.; Road, S. H.; Park, M. On the role of the surface oxygen species during A-H (A = C, N, O) bond activation: a density functional theory study. Chemical Communications 2015, 51, 2621–2624. (33) Calle-Vallejo, F.; Loffreda, D.; Koper, M. T. M.; Sautet, P. Introducing structural sensitivity into adsorption-energy scaling relations by means of coordination numbers. Nature Chemistry 2015, 7, 403–410. (34) Sutton, J. E.; Vlachos, D. G. Building large microkinetic models with first-principles’ accuracy at reduced computational cost. Chemical Engineering Science 2015, 121, 190–199. (35) Medford, A. J.; Shi, C.; Hoffmann, M. J.; Lausche, A. C.; Fitzgibbon, S. R.; Bligaard, T.; Nørskov, J. K. CatMAP: A Software Package for Descriptor-Based Microkinetic Mapping of Catalytic Trends. Catalysis Letters 2015, 145, 794–807.

34 ACS Paragon Plus Environment

Page 34 of 43

Page 35 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(36) Mhadeshwar, A. B.; Kitchin, J. R.; Barteau, M. A.; Vlachos, D. G. The role of adsorbate-adsorbate interactions in the rate controlling step and the most abundant reaction intermediate of NH3 decomposition on RU. Catalysis Letters 2004, 96, 13–22. (37) Broadbelt, L. J.; Stark, S. M.; Klein, M. T. Computer-generated Pyrolysis Modelingon-the-fly Generation of Species, Reactions, and Rates. Industrial & Engineering Chemistry Research 1994, 33, 790–799. (38) Susnow, R. G.; Dean, A. M.; Green, W. H.; Peczak, P.; Broadbelt, L. J. Rate-based construction of kinetic models for complex systems. Journal of Physical Chemistry A 1997, 101, 3731–3740. (39) Broadbelt, L. J.; Pfaendtner, J. Lexicography of kinetic modeling of complex reaction networks. AIChE Journal 2005, 51, 2112–2121. (40) Mitsos, A.; Oxberry, G. M.; Barton, P. I.; Green, W. H. Optimal automatic reaction and species elimination in kinetic mechanisms. Combustion and Flame 2008, 155, 118–132. (41) Sharma, S.; Harper, M. R.; Green, W. H. Modeling of 1,3-hexadiene, 2,4-hexadiene and 1,4-hexadiene-doped methane flames: Flame modeling, benzene and styrene formation. Combustion and Flame 2010, 157, 1331–1345. (42) Van Geem, K. M.; Pyl, S. P.; Marin, G. B.; Harper, M. R.; Green, W. H. Accurate High-Temperature Reaction Networks for Alternative Fuels: Butanol Isomers. Industrial & Engineering Chemistry Research 2010, 49, 10399–10420. (43) Hansen, N.; Harper, M. R.; Green, W. H. High-temperature oxidation chemistry of nbutanol - experiments in low-pressure premixed flames and detailed kinetic modeling. Physical Chemistry Chemical Physics 2011, 13, 20262–20274.

35 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(44) Harper, M. R.; Van Geem, K. M.; Pyl, S. P.; Marin, G. B.; Green, W. H. Comprehensive reaction mechanism for n-butanol pyrolysis and combustion. Combustion and Flame 2011, 158, 16–41. (45) Battin-Leclerc, F.; Blurock, E.; Bounaceur, R.; Fournet, R.; Glaude, P.-A.; Herbinet, O.; Sirjean, B.; Warth, V. Towards cleaner combustion engines through groundbreaking detailed chemical kinetic models. Chemical Society reviews 2011, 40, 4762– 82. (46) Magoon, G. R.; Aguilera-Iparraguirre, J.; Green, W. H.; Lutz, J. J.; Piecuch, P.; Wong, H. W.; Oluwole, O. O. Detailed chemical kinetic modeling of JP-10 (exotetrahydrodicyclopentadiene) high-temperature oxidation: Exploring the role of biradical species in initial decomposition steps. International Journal of Chemical Kinetics 2012, 44, 179–193. (47) Hansen, N.; Merchant, S. S.; Harper, M. R.; Green, W. H. The predictive capability of an automatically generated combustion chemistry mechanism: Chemical structures of premixed iso-butanol flames. Combustion and Flame 2013, 160, 2343–2351. (48) Magoon, G. R.; Green, W. H. Design and implementation of a next-generation software interface for on-the-fly quantum and force field calculations in automated reaction mechanism generation. Computers & Chemical Engineering 2013, 52, 35–45. (49) Merchant, S. S.; Zanoelo, E. F.; Speth, R. L.; Harper, M. R.; Van Geem, K. M.; Green, W. H. Combustion and pyrolysis of iso-butanol: Experimental and chemical kinetic modeling study. Combustion and Flame 2013, 160, 1907–1929. (50) Allen, J. W.; Scheer, A. M.; Gao, C. W.; Merchant, S. S.; Vasu, S. S.; Welz, O.; Savee, J. D.; Osborn, D. L.; Lee, C.; Vranckx, S. et al. A coordinated investigation of the combustion chemistry of diisopropyl ketone, a prototype for biofuels produced by endophytic fungi. Combustion and Flame 2014, 161, 711–724. 36 ACS Paragon Plus Environment

Page 36 of 43

Page 37 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(51) Prozument, K.; Suleimanov, Y. V.; Buesser, B.; Oldham, J. M.; Green, W. H.; Suits, A. G.; Field, R. W. A Signature of Roaming Dynamics in the Thermal Decomposition of Ethyl Nitrite: Chirped-Pulse Rotational Spectroscopy and Kinetic Modeling. Journal of Physical Chemistry Letters 2014, 5, 3641–3648. (52) Vandewiele, N. M.; Magoon, G. R.; Van Geem, K. M.; Reyniers, M. F.; Green, W. H.; Marin, G. B. Experimental and Modeling Study on the Thermal Decomposition of Jet Propellant-10. Energy & Fuels 2014, 28, 4976–4985. (53) Gao, C. W.; Allen, J. W.; Green, W. H.; West, R. H. Reaction Mechanism Generator: automatic construction of chemical kinetic mechanisms. Computer Physics Communications 2016, 203, 212–225. (54) West, R.; Gao, C.; Allen, J.; Vandewiele, N.; gmagoon,; Buesser, B.; Bhoorasingh, P.; Han, K.; nyee,; Belinda, et al. rwest/RMG-Py: RMG-Cat version 1.0. 2017; https: //doi.org/10.5281/zenodo.290119. (55) Goldsmith, C.; West, R. rwest/RMG-database: Database for RMG-Cat version 1.0. 2017; https://doi.org/10.5281/zenodo.290120. (56) Blaylock, D. W.; Ogura, T.; Green, W. H.; Beran, G. J. O. Computational Investigation of Thermochemistry and Kinetics of Steam Methane Reforming on Ni(111) under Realistic Conditions. The Journal of Physical Chemistry C 2009, 113, 4898–4908. (57) Hill, T. L. An Introduction to Statistical Thermodynamics; Dover: New York, 1960. (58) Cohen, N.; Benson, S. W. Estimation of Heats of Formation or Organic-Compounds by Additivity Methods. Chemical Reviews 1993, 93, 2419–2438. (59) Kua, J.; Faglioni, F.; Goddard, W. A. Thermochemistry for Hydrocarbon Intermediates Chemisorbed on Metal Surfaces: CHn−m (CH3 )m with n = 1, 2, 3 and m ≤ n on

37 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Pt, Ir, Os, Pd, Rh, and Ru. Journal of the American Chemical Society 2000, 122, 203902321. (60) Salciccioli, M.; Chen, Y.; Vlachos, D. G. Density Functional Theory-Derived Group Additivity and Linear Scaling Methods for Prediction of Oxygenate Stability on Metal Catalysts: Adsorption of Open-Ring Alcohol and Polyol Dehydrogenation Intermediates on Pt-Based Metals. Journal of Physical Chemistry C 2010, 114, 20155–20166. (61) Salciccioli, M.; Edie, S. M.; Vlachos, D. G. Adsorption of Acid, Ester, and Ether Functional Groups on Pt: Fast Prediction of Thermochemical Properties of Adsorbed Oxygenates via DFT-Based Group Additivity Methods. The Journal of Physical Chemistry C 2012, 116, 1873–1886. (62) Vorotnikov, V.; Wang, S.; Vlachos, D. G. Group Additivity for Estimating Thermochemical Properties of Furanic Compounds on Pd(111). Industrial & Engineering Chemistry Research 2014, 53, 11929–11938. (63) Vorotnikov, V.; Vlachos, D. G. Group Additivity and Modified Linear Scaling Relations for Estimating Surface Thermochemistry on Transition Metal Surfaces: Application to Furanics. The Journal of Physical Chemistry C 2015, 119, 10417–10426. (64) Goldsmith, C. F. Estimating the Thermochemistry of Adsorbates Based Upon GasPhase Properties. Topics in Catalysis 2012, 55, 366–375. (65) Klippenstein, S. J.; Marcus, R. A. High pressure rate constants for unimolecular dissociation/free radical recombination: Determination of the quantum correction via quantum Monte Carlo path integration. The Journal of Chemical Physics 1987, 87, 3410. (66) Georgievskii, Y.; Klippenstein, S. J. Variable reaction coordinate transition state theory: Analytic results and application to the C2 H3 +H → C2 H3 reaction. The Journal of Chemical Physics 2003, 118, 5442–5442. 38 ACS Paragon Plus Environment

Page 38 of 43

Page 39 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(67) Goldsmith, C. F.; Green, W. H.; Klippenstein, S. J. Role of O2 + QOOH in lowtemperature ignition of propane. 1. Temperature and pressure dependent rate coefficients. The Journal of Physical Chemistry A 2012, 116, 3325–46. (68) Goldsmith, C. F.; Harding, L. B.; Georgievskii, Y.; Miller, J. A.; Klippenstein, S. J. Temperature and Pressure-Dependent Rate Coefficients for the Reaction of Vinyl Radical with Molecular Oxygen. Journal of Physical Chemistry A 2015, 119, 7766–79. (69) Goldsmith, C. F.; Burke, M. P.; Georgievskii, Y.; Klippenstein, S. J. Effect of nonthermal product energy distributions on ketohydroperoxide decomposition kinetics. Proceedings of the Combustion Institute 2015, 35, 283–290. (70) Korup, O.; Goldsmith, C. F.; Weinberg, G.; Geske, M.; Kandemir, T.; Schlogl, R.; Horn, R. Catalytic partial oxidation of methane on platinum investigated by spatial reactor profiles, spatially resolved spectroscopy, and microkinetic modeling. Journal of Catalysis 2013, 297, 1–16. (71) Goldsmith, C. F.; Magoon, G. R.; Green, W. H. Database of Small Molecule Thermochemistry for Combustion. The Journal of Physical Chemistry A 2012, (72) Allen, J. W. Predictive Chemical Kinetics: Enabling Automatic Mechanism Generation and Evaluation. Thesis, 2013. (73) Jalan, A.; West, R. H.; Green, W. H. An Extensible Framework for Capturing Solvent Effects in Computer Generated Kinetic Models. Journal of Physical Chemistry B 2013, 117, 2955–2970. (74) Harper, M. R. Automated reaction mechanism generation: data collaboration, Heteroatom implementation, and model validation. PhD Thesis, Massachusetts Institute of Technology, 2011.

39 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(75) Bhoorasingh, P. L.; West, R. H. Transition state geometry prediction using molecular group contributions. Physical Chemistry Chemical Physics 2015, 17, 32173–32182. (76) Kee, R. J.; Coltrin, M. E.; Glarborg, P. Chemically Reacting Flow: Theory and Practice, 1st ed.; John Wiley & Sons: Hoboken, New Jersey, 2003. (77) Delgado, K.; Maier, L.; Tischer, S.; Zellner, A.; Stotz, H.; Deutschmann, O. Surface Reaction Kinetics of Steam- and CO2 -Reforming as Well as Oxidation of Methane over Nickel-Based Catalysts. Catalysts 2015, 5, 871–904. (78) Dunker, A. M. The Decoupled Direct Method for Calculating Sensitivity Coefficients in Chemical-Kinetics. Journal of Chemical Physics 1984, 81, 2385–2393. (79) Coffee, T. P.; Heimerl, J. M. Sensitivity Analysis for Premixed, Laminar, Steady-State Flames. Combustion and Flame 1983, 50, 323–340. (80) Rabitz, H.; Kramer, M.; Dacol, D. Sensitivity Analysis in Chemical-Kinetics. Annual Review of Physical Chemistry 1983, 34, 419–461. (81) Yetter, R.; Eslava, L. A.; Dryer, F. L.; Rabitz, H. Elementary and Derived Sensitivity Information in Chemical-Kinetics. Journal of Physical Chemistry 1984, 88, 1497– 1507. (82) Caracotsios, M.; Stewart, W. E. Sensitivity Analysis of Initial-Value Problems with Mixed ODES and Algebraic Equations. Computers & Chemical Engineering 1985, 9, 359–365. (83) Leis, J. R.; Kramer, M. A. Sensitivity Analysis of Systems of Differential and Algebraic Equations. Computers & Chemical Engineering 1985, 9, 93–96. (84) Reuven, Y.; Smooke, M. D.; Rabitz, H. Sensitivity Analysis of Boundary-Value Problesm - Applications to Nonlinear Reaction Diffusion-Systems. Journal of Computational Physics 1986, 64, 27–55. 40 ACS Paragon Plus Environment

Page 40 of 43

Page 41 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(85) Smooke, M. D.; Rabitz, H.; Reuven, Y.; Dryer, F. L. Application of Sensitivity Analysis to Premixed Hydrogen-Air Flames. Combustion Science and Technology 1988, 59, 295–319. (86) Miller, J. A.; Kee, R. J.; Westbrook, C. K. Chemical-Kinetics and Combustion Modeling. Annual Review of Physical Chemistry 1990, 41, 345–387. (87) Campbell, C. T. Future Directions and Industrial Perspectives Micro- and macrokinetics: Their relationship in heterogeneous catalysis. Topics in Catalysis 1994, 1, 353–366. (88) Campbell, C. T. Finding the Rate-Determining Step in a Mechanism: Comparing DeDonder Relations with the “Degree of Rate Control”. J. Catal. 2001, 204, 520–524. (89) Stegelmann, C.; Andreasen, A.; Campbell, C. T. Degree of rate control: How much the energies of intermediates and transition states control rates. Journal of the American Chemical Society 2009, 131, 8077–8082. (90) Allen, J. W. PyDAS: A Python wrapper for the DASSL, DASPK, and DASKR differential algebraic system solvers. https://github.com/jwallen/PyDAS. (91) Petzold, L. DASSL and DASPK Differential Algebraic Equation Solvers. http://www. cs.ucsb.edu/%7Ecse/software.html. (92) Ashcraft, R. W.; Green, W. H. Thermochemical properties and group values for nitrogen-containing molecules. Journal of Physical Chemistry A 2008, 112, 9144– 9152. (93) Ashcraft, R. W.; Raman, S.; Green, W. H. Predicted reaction rates of HxNyOz intermediates in the oxidation of hydroxylamine by aqueous nitric acid. Journal of Physical Chemistry A 2008, 112, 7577–7593.

41 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(94) Jalan, A.; Ashcraft, R. W.; West, R. H.; Green, W. H. Predicting solvation energies for kinetic modeling. Annual Reports Section "C" (Physical Chemistry) 2010, 106, 211. (95) Timko, M. T.; Schmois, E.; Patwardhan, P.; Kida, Y.; Class, C. A.; Green, W. H.; Nelson, R. K.; Reddy, C. M. Response of Different Types of Sulfur Compounds to Oxidative Desulfurization of Jet Fuel. Energy & Fuels 2014, 28, 2977–2983. (96) Slakman, B.; West, R. H. Automatic Mechanism Generation for Liquid-Phase Reaction Kinetics. AIChE Annual Meeting. 2013. (97) Wehinger, G. D.; Kraume, M.; Berg, V.; Korup, O.; Mette, K.; Schlogl, R.; Behrens, M.; Horn, R. Investigating dry reforming of methane with spatial reactor profiles and particle-resolved CFD simulations. AIChE Journal 2016, 62, 4436–4452. (98) West, R.; Goldsmith, F. RMG-Cat version 1.0 examples for automated reaction mechanism generation of heterogeneous catalysis. 2017; http://doi.org/10.6084/m9. figshare.4640098.v1. (99) Temel, B.; Meskine, H.; Reuter, K.; Scheffler, M.; Metiu, H. Does phenomenological kinetics provide an adequate description of heterogeneous catalytic reactions? Journal of Chemical Physics 2007, 126, 204711–204723. (100) Meskine, H.; Matera, S.; Scheffler, M.; Reuter, K.; Metiu, H. Examination of the concept of degree of rate control by first-principles kinetic Monte Carlo simulations. Surface Science 2009, 603, 1724–1730.

42 ACS Paragon Plus Environment

Page 42 of 43

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Graphical TOC Entry 10 0

CH4 CO2

CO H2

vacant

10 -1

H*

CH*

CO*

10 -2

site fraction

Page 43 of 43

C* CH 3 C*O

10 -3

CH 3 C* 10 -4

O*

10 -5

C 2 HO*

CH 2 *

O*

CHO 2 *

10 -6 10 -6

10 -5

10 -4

10 -3

10 -2

time (s)

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX

43 ACS Paragon Plus Environment

10 -1

10 0