First-principles–based reaction kinetics from reactive molecular dynamics simulations: Application to hydrogen peroxide decomposition

This paper presents our vision of how to use in silico approaches to extract the reaction mechanisms and kinetic parameters for complex condensed-phase chemical processes that underlie important technologies ranging from combustion to chemical vapor deposition. The goal is to provide an analytic description of the detailed evolution of a complex chemical system from reactants through various intermediates to products, so that one could optimize the efficiency of the reactive processes to produce the desired products and avoid unwanted side products. We could start with quantum mechanics (QM) to ensure an accurate description; however, to obtain useful kinetics we need to average over ∼10-nm spatial scales for ∼1 ns, which is prohibitively impractical with QM. Instead, we use the reactive force field (ReaxFF) trained to fit QM to carry out the reactive molecular dynamics (RMD). We focus here on showing that it is practical to extract from such RMD the reaction mechanisms and kinetics information needed to describe the reactions analytically. This analytic description can then be used to incorporate the correct reaction chemistry from the QM/ReaxFF atomistic description into larger-scale simulations of ∼10 nm to micrometers to millimeters to meters using analytic approaches of computational fluid dynamics and/or continuum chemical dynamics. In the paper we lay out the strategy to extract the mechanisms and rate parameters automatically without the necessity of knowing any details of the chemistry. We consider this to be a proof of concept. We refer to the process as RMD2Kin (reactive molecular dynamics to kinetics) for the general approach and as ReaxMD2Kin (ReaxFF molecular dynamics to kinetics) for QM-ReaxFF–based reaction kinetics.