Modeling O(³P) and Ar scattering from the ionic liquid [emim] [NO₃] at 5 eV with hybrid QM/MM molecular dynamics

In this paper, we develop a hybrid QM/MM-MD direct dynamics method for the study of Ar and O scattering from room-temperature ionic liquids (RTIL) at hyperthermal energies, with the goal of providing an understanding of the reactivity of ionic liquids as potential hypergolic fuels. The RTIL is chosen to be 1 -ethyl3-methylimidazolium nitrate ([emim][N03]) as the bulk properties of this have been analyzed experimentally and theoretically, and a force field has been developed. Within our simulation time of 7.3 ps, Ar collisions with the liquid surface cause 2-3 proton transfers to occur from the imidazolium ring to the anion, and a small fraction of these events leads to the desorption of neutral "emim" or HNO₃ species. These chargetransfer processes also occur for O(³P) collisions, but in addition, there is a rich amount of chemical change at the surface leading to four main collision types, reaction with NO₃ to form species like [NO₂] and O₂, O addition to the imidazolium ring causing ring scission, H abstractions from the hydrocarbon chains (ethyl or methyl) or imidazolium ring generating OH, and O atom inelastic scattering. Within the MSINDO model, the major scattered products (all are neutral) from the O collisions, with their probabilities in parentheses, are O (0.13), emim (0.12), HNO₂ (0.05), HNO₃ (0.04), and OH (0.04). Of these products, only OH is unique to O atom collisions; the HNO₂ and HNO₃ products also arise when the incident atom is Ar or O.