Theoretical methodology for prediction of tropospheric oxidation of dimethyl phosphonate and dimethyl methylphosphonate

Rate constants for the reactions of OH radicals with dimethyl phosphonate [DMHP, (CH₃O)₂P(O)H] and dimethyl methylphosphonate [DMMP, (CH₃O)₂P(O)CH₃] have been calculated by ab initio structural methods and semiclassical dynamics modeling and compared with experimental measurements over the temperature range 250-350 K. The structure and energetics of reactants and transition structures are determined for all hydrogen atom abstraction pathways that initiate the atmospheric oxidation mechanism. Structures are obtained at the CCSD/6-31++G* level of chemical theory, and the height of the activation barrier is determined by a variant of the G2MP2 method. A Transfer Hamiltonian is used to compute the minimum energy path in the neighborhood of the transition state (TS). This calculation provides information about the curvature of the potential energy surface in the neighborhood of the TS, as well as the internal forces that are needed by the semiclassical flux-flux autocorrelation function (SCFFAF) dynamics model used to compute the temperature-dependent reaction rate constants for the various possible abstraction pathways. The computed temperature-dependent rate curves frequently lie within the experimental error bars.