Automated Skeleton Network Generation for ReaxFF Molecular Dynamics Simulations of Hydrocarbon Fuel Pyrolysis and...

In this study, we present an automated approach of rate-based skeleton network generation for ReaxFF MD simulation (RxMD-SN) for deriving the reaction kinetic mechanism of large hydrocarbon fuels in pyrolysis and oxidation from large-scale ReaxFF MD simulations. The approach contains the statistical calculation of reaction rate constants and the generation of skeleton reaction networks using a rate-based algorithm. The RxMD-SN method takes advantage of reaction flux ranking at a small time interval in terms of temporal reaction rate to extract the core reaction networks, which allows for keeping the rare reaction events that may be dominant in a certain period of the reaction network. The kinetic models derived from ReaxFF MD simulation in CH₄ oxidation can reproduce what was obtained in the ReaxFF MD simulation, which demonstrates the capability of RxMD-SN in capturing the global reaction kinetics. An evaluation of reaction rate constants indicates that close kinetic parameters are shared for n-octane oxidation of similar reaction classes, shared oxidation reactions of CH₄ against n-heptane, and shared pyrolysis reactions of the RP-3 surrogate fuel against n-heptane. This capability of RxMD-SN is particularly beneficial in meeting the challenges in characterizing the oxidation reaction kinetics of large hydrocarbon molecules. RxMD-SN approach is potentially a general approach in chemical kinetics modeling on the basis of ReaxFF MD simulations.