Quasiclassical Trajectory Studies of the O(3P)+CX4(vk=0, 1) → OX(v)+CX3(n1n2n3n4) [X = H and D] Reactions on an Ab Initio Potential Energy Surface

We report quasiclassical trajectory calculations of the integral and differential cross sections and the mode-specific product state distributions for the central-barrier O(P-3) + CH4/CD4(v(k) = 0, 1) [k = 1, 2, 3, 4] reactions using a full-dimensional ab initio potential energy surface. The mode-specific vibrational distributions for the polyatomic methyl products are obtained by doing a normal-mode analysis in the Eckart frame, followed by standard histogram binning (HB) and energy-based Gaussian binning (1GB). The reactant bending excitations slightly enhance the reactivity, whereas stretching excitations activate the reaction more efficiently. None of the reactant vibrational excitations is as efficient as an equivalent amount of translational energy to promote the reactions. The excitation functions without product zero-point energy (ZPE) constraint are in good agreement with previous 8-dimensional quantum mechanical (QM) results for the ground-state and stretching-excited O + CH4 reactions, whereas for the bending-excited reactions the soft ZPE constraint, which is applied to the sum of the product vibrational energies, provides better agreement with the QM cross sections. All angular distributions show the dominance of backward scattering indicating a direct rebound mechanism, in agreement with experiment. The title reactions produce mainly OH/OD(v = 0) products for all the initial states. HB significantly overestimates the populations of OH/OD(v = 1), especially in the energetic threshold regions, whereas 1GB provides physically correct results. The CH3/CD3 vibrational distributions show dominant populations for ground (v = 0), umbrella-excited (v(2) = 1, 2), in-plane-bending-excited (v(4) = 1), and v(2) + v(4) methyl product states. Neither translational energy nor reactant vibrational excitation transfers significantly into product vibrations.