Energetics and dynamics of the reactions of O(3P) with dimethyl methylphosphonate and sarin

Electronic structure and molecular dynamics calculations were performed on the reaction systems O((3)P) + sarin and O((3)P) + dimethyl methylphosphonate (DMMP), a sarin simulant. Transition state geometries, energies, and heats of reaction for the major reaction pathways were determined at several levels of theory, including AM1, B3LYP/6-311+G(d,p), and CBS-QB3. The major reaction pathways for both systems are similar and include H-atom abstraction, H-atom elimination, and methyl elimination, in rough order from low to high energy. The H-atom abstraction channels have fairly low barriers (approximately 10 kcal mol(-1)) and are close to thermoneutral, while the other channels have relatively high energy barriers (>40 kcal mol(-1)) and a wide range of reaction enthalpies. We have also found a two-step pathway leading to methyl elimination through O-atom attack on the phosphorus atom for DMMP and sarin. For sarin, the two-step methyl elimination pathway is significantly lower in energy than the single-step pathway. We also present results of O((3)P) + sarin and O((3)P) + DMMP reaction cross sections over a broad range of collision energies (2-10 km s(-1) collision velocities) obtained using the direct dynamics method with an AM1 semiempirical potential. These excitation functions are intended as an approximate guide to future hyperthermal measurements, which to our knowledge have not yet examined either of these systems. The reaction barriers, reaction enthalpies, transition state structures, and excitation functions are generally similar for DMMP and sarin, with some moderate differences for methyl elimination energetics, which indicates DMMP will likely be a good substitute for sarin in many O((3)P) chemical investigations.