Joshua W. Allen

RPMDrate: Bimolecular chemical reaction rates from ring polymer molecular dynamics ☆

United States. Dept. of Energy (Office of Basic Energy Sciences under the Energy Frontier Research Center for Combustion Science (Grant No. DE-SC0001198))

Communication: Full dimensional quantum rate coefficients and kinetic isotope effects from ring polymer molecular dynamics for a seven-atom reaction OH + CH4 → CH3 + H2O

The kinetic isotope effect (KIE) of the seven-atom reactions OH + CH4 → CH3 + H2O and OH + CD4 → CD3 + HDO over the temperature range 200-1000 K is investigated using ring polymer molecular dynamics (RPMD) on a full-dimensional potential energy surface. A comparison of RPMD with previous theoretical results obtained using transition state theory shows that RPMD is a more reliable theoretical approach for systems with more than 6 atoms, which provides a predictable level of accuracy. We show that the success of RPMD is a direct result of its independence of the choice of transition state dividing surface, a feature that is not shared by any of the transition state theory-based methods. Our results demonstrate that RPMD is a prospective method for studies of KIEs for polyatomic reactions for which rigorous quantum mechanical calculations are currently impossible.

Chemically Activated Formation of Organic Acids in Reactions of the Criegee Intermediate with Aldehydes and Ketones

Reactions of the Criegee intermediate (CI, ˙CH2OO˙) are important in atmospheric ozonolysis models. In this work, we compute the rates for reactions between ˙CH2OO˙ and HCHO, CH3CHO and CH3COCH3 leading to the formation of secondary ozonides (SOZ) and organic acids. Relative to infinitely separated reactants, the SOZ in all three cases is found to be 48-51 kcal mol(-1) lower in energy, formed via 1,3-cycloaddition of ˙CH2OO˙ across the C=O bond. The lowest energy pathway found for SOZ decomposition is intramolecular disproportionation of the singlet biradical intermediate formed from cleavage of the O-O bond to form hydroxyalkyl esters. These hydroxyalkyl esters undergo concerted decomposition providing a low energy pathway from SOZ to acids. Geometries and frequencies of all stationary points were obtained using the B3LYP/MG3S DFT model chemistry, and energies were refined using RCCSD(T)-F12a/cc-pVTZ-F12 single-point calculations. RRKM calculations were used to obtain microcanonical rate coefficients (k(E)) and the reservoir state method was used to obtain temperature and pressure dependent rate coefficients (k(T, P)) and product branching ratios. At atmospheric pressure, the yield of collisionally stabilized SOZ was found to increase in the order HCHO < CH3CHO < CH3COCH3 (the highest yield being 10(-4) times lower than the initial ˙CH2OO˙ concentration). At low pressures, chemically activated formation of organic acids (formic acid in the case of HCHO and CH3COCH3, formic and acetic acid in the case of CH3CHO) was found to be the major product channel in agreement with recent direct measurements. Collisional energy transfer parameters and the barrier heights for SOZ reactions were found to be the most sensitive parameters determining SOZ and organic acid yield.

Reply to comment on “Automatic estimation of pressure-dependent rate coefficients” (J. W. Allen, C. F. Goldsmith, and W. H. Green, Phys. Chem. Chem. Phys., 2011, 14, 1131–1155)

Sources of error in determining pressure-dependent rate coefficients are discussed in the context of automatic mechanism generation. Often the approximations made in converting the master equations into phenomenological rate coefficients are not the largest source of error. However, in some cases some of the methods are seriously inaccurate. Automation of Miller and Klippenstein’s accurate CSE method, including automatic identification of equilibrated species, would be very helpful.