Kinetics and Thermodynamics of Reactions Involving Criegee Intermediates: An Assessment of Density Functional Theory and Ab Initio Methods Through Comparison with CCSDT(Q)/CBS Data
Cameron D. Smith,Amir Karton
DOI: https://doi.org/10.1002/jcc.26106
2019-11-21
Journal of Computational Chemistry
Abstract:<p>Reactions involving Criegee intermediates (CIs, R<sub>1</sub>R<sub>2</sub>COO) are important in atmospheric ozonolysis models. In recent years, density functional theory (DFT) and CCSD(T)‐based ab initio methods are increasingly being used for modeling reaction profiles involving CIs. We obtain highly accurate CCSDT(Q)/CBS reaction energies and barrier heights for ring‐closing reactions involving atmospherically important CIs (R<sub>1</sub>/R<sub>2</sub> = H, Me, OH, OMe, F, CN, cyclopropene, ethylene, acetaldehyde, and acrolein). We use this benchmark data to evaluate the performance of DFT, double‐hybrid DFT (DHDFT), and ab initio methods for the kinetics and thermodynamics of these reactions. We find that reaction energies are more challenging for approximate theoretical procedures than barrier heights. Overall, taking both reaction energies and barrier heights into account, only one of the 58 considered DFT methods (the meta‐GGA MN12‐L) attains near chemical accuracy, with root‐mean‐square deviations (RMSDs) of 3.5 (barrier heights) and 4.7 (reaction energies) kJ mol<sup>−1</sup>. Therefore, MN12‐L is recommended for investigations where CCSD(T)‐based methods are not computationally feasible. For reaction barrier heights performance does not strictly follow Jacob's Ladder, for example, DHDFT methods do not perform better than conventional DFT methods. Of the ab initio methods, the cost‐effective CCSD(T)/CBS(MP2) approach gives the best performance for both reaction energies and barrier heights, with RMSDs of 1.7 and 1.4 kJ mol<sup>−1</sup>, respectively. All the considered Gaussian‐<i>n</i> methods show good performance with RMSDs below the threshold of chemical accuracy for both reaction energies and barrier heights, where G4(MP2) shows the best overall performance with RMSDs of 2.9 and 1.5 kJ mol<sup>−1</sup>, respectively. © 2019 Wiley Periodicals, Inc.</p>
chemistry, multidisciplinary