Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.
Citation: K. Chenoweth, A.C.T. van Duin, and W.A. Goddard (2008), "ReaxFF Reactive Force Field for Molecular Dynamics Simulations of Hydrocarbon Oxidation", The Journal of Physical Chemistry A, 112(5), 1040-1053. DOI: 10.1021/jp709896w.
Abstract: To investigate the initial chemical events associated with high-temperature gas-phase oxidation of hydrocarbons, we have expanded the ReaxFF reactive force field training set to include additional transition states and chemical reactivity of systems relevant to these reactions and optimized the force field parameters against a quantum mechanics (QM)-based training set. To validate the ReaxFF potential obtained after parameter optimization, we performed a range of NVT−MD simulations on various hydrocarbon/O2 systems. From simulations on methane/O2, o-xylene/O2, propene/O2, and benzene/O2 mixtures, we found that ReaxFF obtains the correct reactivity trend (propene > o-xylene > methane > benzene), following the trend in the C−H bond strength in these hydrocarbons. We also tracked in detail the reactions during a complete oxidation of isolated methane, propene, and o-xylene to a CO/CO2/H2O mixture and found that the pathways predicted by ReaxFF are in agreement with chemical intuition and our QM results. We observed that the predominant initiation reaction for oxidation of methane, propene, and o-xylene under fuel lean conditions involved hydrogen abstraction of the methyl hydrogen by molecular oxygen forming hydroperoxyl and hydrocarbon radical species. While under fuel rich conditions with a mixture of these hydrocarbons, we observed different chemistry compared with the oxidation of isolated hydrocarbons including a change in the type of initiation reactions, which involved both decomposition of the hydrocarbon or attack by other radicals in the system. Since ReaxFF is capable of simulating complicated reaction pathways without any preconditioning, we believe that atomistic modeling with ReaxFF provides a useful method for determining the initial events of oxidation of hydrocarbons under extreme conditions and can enhance existing combustion models.
See Computed Properties Notes: The file "ffield.reax.CHO_2008" was provided by Adri van Duin. From Prof. van Duin: "The ffield-file contains the force field parameters; this file is readable by LAMMPS." The ReaxFF manual (including file formatting information) was obtained from http://www.wag.caltech.edu/home/duin/manual.html. All files were posted with Prof. van Duin's approval. The standalone ReaxFF program is available without charge for academic users by emailing him. File(s):