Kinetics of the Strongly Correlated CH3O + O2 Reaction: The

5 days ago - methods is challenging because the transition states and the .... It is a difficult challenge to use theoretical methods to model the CH3...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Rhode Island | University Libraries

Article

Kinetics of the Strongly Correlated CH3O + O2 Reaction: The Importance of Quadruple Excitations in Atmospheric and Combustion Chemistry Bo Long, Junwei Lucas Bao, and Donald G. Truhlar J. Am. Chem. Soc., Just Accepted Manuscript • DOI: 10.1021/jacs.8b11766 • Publication Date (Web): 13 Dec 2018 Downloaded from http://pubs.acs.org on December 13, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society Dec. 11, 2018

Kinetics of the Strongly Correlated CH3O + O2 Reaction: The Importance of Quadruple Excitations in Atmospheric and Combustion Chemistry Bo Long, a , b * Junwei Lucas Bao,b and Donald G. Truhlarb∗ a

College of Materials Science and Engineering, Guizhou Minzu University, Guiyang, 550025, China

b

Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA ABSTRACT: Kinetics measurements on radical–radical reactions are often unavailable experimentally, and obtaining quantitative rate constants for such reactions by theoretical methods is challenging because the transition states and the reactants are often strongly correlated. Treating strongly correlated systems by coupled cluster theory limited to single, double, and triple connected excitations is often inadequate. We therefore use a new method, called GMM(P), for extrapolation to the complete configuration interaction limit to go beyond triple excitations and in particular to approximate the CCSDTQ(P)/CBS limit. Here, we present this method and use it to investigate the CH3O + O2 reaction. The contribution of connected quadruple excitations to the barrier height energy is found to be –3.13 kcal/mol, and adding a quasiperturbative calculation of the effect of connected pentuple excitations brings the post-connected-triples contributions to –3.44 kcal/mol, which corresponds to Boltzmann factors that increase calculated rate constants by factors of 1.0 × 103, 3.3 × 102, and 18 at 250 K, 298 K, and 600 K, respectively. We present rate constants for temperatures from 250 K to 2000 K, and we find that the Arrhenius activation energy increases from 0.58 to 9.68 kcal/mol over this range. We also find reasonably good accuracy for the barrier height with the MN15-L exchange-correlation functional, and we calculate rate constants by a combination of GMM(P) and MN15-L electronic structure calculations and conventional and variational transition state theory, in particular canonical variational theory with small-curvature tunneling.. The present findings have broad implications for obtaining quantitative rate constants for complex reaction systems in atmospheric and combustion chemistry.

∗ To whom correspondence should be addressed. E-mail: [email protected] (Bo Long),

[email protected] (Donald G. Truhlar)

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 20 2

I. INTRODUCTION Alkoxy radicals (RO), which are formed via the oxidation of volatile organic compounds, are key intermediates in the atmosphere and in combustion.1,2 Alkoxy radicals react with O2 and undergo unimolecular decomposition and isomerization in the atmosphere, and these reactions have a significant impact on the formation of ozone and secondary organic aerosols.3,4,5,6,7 Due to the importance of alkoxy radicals in the atmosphere, obtaining the absolute rate constants of the RO + O2 reactions is very important for fully estimate the atmospheric fates of alkoxy radicals, but obtaining these rates for the larger and functionalized alkoxy radicals from experiment is very difficult.8 Therefore we turn to computational chemistry to predict the gas-phase reaction kinetics, and our goal is to obtain experimental accuracy or better. The simplest alkoxy radical is the methoxy radical (CH3O). The dominant sink of CH3O is reaction with O2,9 which has been extensively studied both experimentally1,2,3,10,11,12,13,14,15,16 and theoretically.17,18,19,20,21 For this reaction of a small radical there are experimental kinetic data from different sources that agree with one another within a factor of 1.5–2. The CH3O + O2 reaction that dominates at low temperature occurs by abstraction of a hydrogen atom to yield HCHO + HO2; this occurs via two transition states;21 these transition states are labeled TS1 and TS2, and they are shown – as reoptimized in the present work–in Figure 1. Calculations in Table S1 with the W3X-L method22 at a geometry optimized by RCCSD(T)-F12a/jun-cc-pVTZ show a difference in classical barrier heights of 9.26 kcal/mol for the two transition states and that the enthalpy of activation at 0 K for TS1 is lower than that of TS2 by 7.49 kcal/mol. (Tables with the prefix S are in supporting information, which should also be consulted for references for basis sets (see Table S2) and explanations of coupled cluster abbreviations.) The difference in classical barrier heights given by the W3X-L calculation is even larger than an earlier estimate19 of 7.5 kcal/mol for the difference in classical barrier heights. Thus TS2 makes a negligible contribution to the CH3O + O2 reaction, and it will not be considered further. The rest of this paper concerns the hydrogen abstraction reaction proceeding through TS1. It is a difficult challenge to use theoretical methods to model the CH3O + O2 reaction because the transition state of the CH3O + O2 reaction is strongly correlated. Strongly correlated systems are systems for which a single configuration state function does not provide a good

ACS Paragon Plus Environment

Page 3 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 3

zero-order description; they are often called multireference systems because their treatment with a low level of excitations (up to only double or triple excitations) often requires a multi-configuration reference function, but they may also be treated with single-configuration reference functions if one includes higher excitations. In the present work we use both single-reference coupled cluster theory and multireference complete active space second-order perturbation theory (CASPT2). Hu and Dibble21 found that calculations by the single-reference CCSD(T) method predict a barrier that is several kcal/mol too high. Setokuchi et al20 investigated the reaction using the G2M(RCC1) composite electronic structure method and obtained a rate constant in reasonable agreement with experiment based on a theoretical method containing no component at a higher level than triple excitations with a modest basis set plus a non-size-consistent empirical high-level correction, but the reliability of these methods in general is not validated. Thus, further investigations are needed to develop and validate theoretical methods that can reliably produce rate constants with experimental accuracy. In this report, we demonstrate the level of theory required to calculate the barrier height for the CH3O + O2 reaction, and we calculate the rate constant by canonical variational transition state theory with small-curvature tunneling (CVT/SCT). The present results not only address a specific gas-phase reaction, but also have broad implications for the treatment of radical–radical reactions with strongly correlated transition states such as the reactions is the classes RO2 + HO223,24,25 and RO2 + OH.26,27

II. METHODS One can often use lower-lower methods for optimizing geometries than for calculating reliable energetics, but nevertheless we calculated geometries at a high level.28 For open-shell doublet (CH3O, HO2, and transition states) and triplet (O2) states, we optimized structures and calculated frequencies with unrestricted explicitly correlated coupled cluster theory with single and double excitations and noniterative triple excitations based on restricted open-shell Hartree– Fock (ROHF) orbitals (RCCSD(T)-F12a)29,30 with the basis sets jun-cc-pVTZ and cc-pVDZ-F12. We used restricted RCCSD(T)-F12a based on based on closed-shell Hartree– Fock orbitals for singlet formaldehyde. To calculate energetics, we designed a new composite scheme to do single-point energy

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 20 4

calculations at an accuracy close to frozen-core CCSDTQ(P)/CBS, where frozen-core means that core orbitals are calculated by Hartree-Fock theory and are doubly occupied in all configurations. The frozen-core approximation is expected to be good for relative energies (barrier heights and energies of reaction), and thus the relative energies should be close to FCI/CBS, where FCI denotes full configuration interaction; FCI/CBS is also called complete configuration interaction, and it is equivalent to solving the many-electron Schrödinger equation. The highest-level of the new composite scheme, formulated to approximate the unrestricted coupled cluster theory with single, double, triple, and quadruple excitations and noniterative quintuple excitations at the complete basis set limit (CCSDTQ(P)/CBS), is called GMM(P). The energy is given by 𝐸 GMM(5) = 𝐸!"#!!"# + ∆𝐸!! ! + ∆𝐸(!)!! + ∆𝐸!!(!) + ∆𝐸(!)!! + ESO

(1)

where MW2-F12 is a method originally developed for a single-point CCSD(T)/CBS energy of benzene in the recently proposed W4-F12 approach31 and was subsequently used for studying the unimolecular reaction of (CH3)2COO;32 the other terms are defined as follows: ∆𝐸!! ! = E(RCCSDT/CBS) – E[RCCSD(T)/CBS]

(2)

∆𝐸(!)!! = E[CCSDT(Q)/CBS] – E(CCSDT/CBS)

(3)

∆𝐸!!(!) = E(CCSDTQ/VDZ(d)] – E[CCSDT(Q)/VDZ(d)]

(4)

∆𝐸(!)!! = E[CCSDTQ(P)/cc-VDZ] – E[CCSDTQ/cc-VDZ]

(5)

and ESO is the lowering of the energy by spin-orbit coupling. In eqs 2 and 3, CBS denotes extrapolation to the complete basis set limit, and this was accomplished by31,33 ∆𝐸! = ∆𝐸!"# +

! !!

(6)

where L is the highest angular momentum in the basis set. In GMM(P), we extrapolate using jun-cc-pVDZ and jun-cc-pVTZ for which L equals 2 and 3, respectively. Coupled cluster calculations based on ROHF orbitals are denoted RCC…, and coupled cluster calculations based on UHF orbitals are denoted CC…. The W3X-L procedure provides a cost-efficient approximation to the RCCSDT(Q)/CBS energy, where CBS denotes extrapolation to a complete basis set. Classical barrier height is the difference in potential energy of the saddle point and the reactants; enthalpy at 0 K is the potential energy plus the zero point vibrational energy; potential energy is the potential energy for nuclear motion, and it is the electronic energy

ACS Paragon Plus Environment

Page 5 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 5

including nuclear repulsion. All calculations include spin-orbit coupling for CH3O (spin-orbit coupling is zero to first order for other species). The spin-orbit coupling ESO of CH3O is estimated to be –0.09 kcal/mol, based on previous experimental34 work (see also refs. 35,36,37). This raised the barrier heights and energies of all reactions by 0.09 kcal/mol. We also define GMMQ in which we do not include ∆𝐸(!)!! and GMM(Q) in which we do not include ∆𝐸!!(!) and ∆𝐸(!)!! .

We found that special care must be exercised to be sure that the reference functions used for the coupled cluster calculations are converged to the physically correct solution. To facilitate this, we used the method of Vaucher and Reiher38 to find the best SCF orbitals by mixing pairs of randomly selected occupied−unoccupied molecular orbitals to search the multiple solutions of the unrestricted Hartree-Fock calculations. Vibrational frequencies for a given calculation (which are needed to calculate finite-temperature enthalpies and free energies) were calculated by the same method as used to optimize the geometry and were then scaled by our standard procedure39 to obtain more accurate zero point energies. The electronic structure and frequency calculations were carried out using the Gaussian 1640 and MN-GFM41 programs for density functional theory and using Molpro 201542 and MRCC43,44 for wave function theory. Rate constants were calculated using the Polyrate 2017-C45 and Gaussrate 2017-B46 dynamics codes. Structures. Tables S3 and S4 show that the transition state structures of TS1 optimized at various levels agree with one another. We conclude that using RCCSD(T)-F12a/jun-cc-pVTZ geometries should contribute negligible error; for brevity we will call this the jun method in the text and tables.

III. CONVERGENCE OF GMM(P) Details of the component calculations of GMM(P) are provided in Tables S5 and S6. We used multiple basis sets to validate that the components of the GMM(P) calculation are converged to the complete basis set limit. Table S5 shows only slight differences (0.13, 0.10 kcal/mol) for the RCCSDT – RCCSD(T) and CCSDT(Q) – CCSDT components between using cc-pV(D, T)Z basis sets to extrapolate to the CBS limit and using jun-cc-pV(D,T)Z basis sets to

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 20 6

extrapolate to the CBS limit; this shows that these two components are well converged to the CBS limit. Table S5 also tells us that the CCSDT(Q) – CCSDT component makes the main contribution to the beyond-CCSD(T) calculations for the CH3O + O2 reaction, and it shows that beyond-CCSD(T) calculations at least up to CCSDTQ(P) are required to obtain quantitative results for the CH3O + O2 reaction. Based on our previous paper,32 the convergence check in Table S5, and experience in the literature with higher-order excitations,47 it is reasonable to expect that relative energies computed by our GMM(P) best estimate are in most cases within a few tenths of a kcal/mol of the frozen-core complete configuration interaction limit.

IV. ENERGIES Zero Point Energies. CH3O is a complex molecule, with Jahn–Teller (JT) and significant spin–orbit coupling.48 CH3O has been extensively investigated by both experimental and theoretical methods.48,49,50,51 Here, we only need to consider the zero-point vibrational energy of CH3O because it is required to compute the enthalpy of activation. Table S7 shows that the scaled zero-point vibrational energy of CH3O is 22.7 kcal/mol by the jun method and is the same by RCCSD(T)-F12a/jun-cc-pVTZ. This value is in excellent agreement with the experimental value of 22.6 kcal/mol.49,50,51 Reaction Barriers. We first consider calculations at the jun geometry. The W2X22 and MW2-F1231,32 procedures provide cost-efficient approximations to the RCCSD(T)/CBS energy, with the latter having a larger basis set. Table 1 shows that these methods give barrier heights and 0 K energies of activation 2.23–2.74 kcal/mol higher than our best estimate. MW3X-L is a new method introduced here that combines MW2-F12 with beyond-CCSD(T) calculations from W3X-L. Table 1 shows that including higher-than-CCSD(T) contributions by MW3X-L lowers the deviations from our best estimate to 0.99 kcal/mol, and including beyond-CCSD(T) contributions by the more expensive W3X-L method lowers this further to 0.48 kcal/mol. The difference in Table 1 between the beyond-CCSD(T) result of GMM(Q) and the CCSD(T) result of MW2-F12 is very large, in particular 3.01 kcal/mol; this also shows that the CH3O + O2 reaction is much more strongly correlated than the OH + SO2,52 CH3O + HF,53 Criegee intermediates + H2O/NH3,32,54,55 unimolecular reactions of Criegee intermediates,32,54,56 and HCl + OH57 reactions because the beyond-CCSD(T) corrections up to CCSDT(Q) in these reactions are less than 1 kcal/mol. In addition, we note that the T1, D1, and T1/D1 diagnostics

ACS Paragon Plus Environment

Page 7 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 7

do not reliably indicate the character of the CH3O + O2 reaction. For example, the T1 diagnostic values are less than 0.044 for CH3O, O2, and TS1 in Table S8; these diagnostics would be interpreted using customary criteria to mean that the CH3O + O2 reaction has only minor multireference character. However, our investigations have shown that the CH3O + O2 reaction is not treated accurately by the methods that usually give good results for single-reference systems. The difference between MW2-F12//RCCSD(T)-F12a/jun-cc-pVTZ and W2X//RCCSD(T)-F12a/jun-cc-pVTZ is about 0.5 kcal/mol, which shows that W2X does not give a good approximation to the complete basis set limit for CCSD(T), although our previous investigation showed that W2X agrees well with MW2-F12 for the unimolecular reaction of (CH3)2COO.32 But some systems are harder and require higher levels than others. The contribution of connected quadruple excitations to the barrier height energy is found to be –3.13 kcal/mol, and adding a quasiperturbative calculation of the effect of connected pentuple excitations brings the post-connected-triples contributions to –3.44 kcal/mol, which corresponds to Boltzmann factors that increase calculated rate constants by factors of 1.0 × 103, 47, and 6 at 250 K, 298 K, and 1000 K, respectively. The GMM(P) benchmark calculations were used to identify a density functional method that gives a semiquantitatively correct description of the reaction. We found that Kohn-Sham density functional theory with the MN15-L exchange–correlation functional58 and the MG3S basis set has a mean unsigned error of only 0.84 kcal/mol (Table 1), which is about 2 kcal/mol lower than that of MW2-F12. This shows that density functional theory with a modern functional is capable of much higher accuracy than CCSD(T) for some strongly correlated systems, a finding that is consistent with our previous results.32,54 This result is especially encouraging because the MN15-L density functional has unique properties. Most exchange-correlation functionals that give good barrier heights have some percentage of Hartree–Fock exchange, which is known to bring in static correlation error, making the functional poorly suited for strongly correlated systems.59 MN15-L is a local functional (no Hartree-Fock exchange) and so it does not have this problem, and yet a recent survey60 found that its mean unsigned error on the 76 diverse barrier heights is only 1.66 kcal/mol, the best performance of 40 local functionals tested (the other local functionals tested have an average

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 20 8

mean unsigned error for these barrier heights of 6.80 kcal/mol). We conclude that MN15-L might be very useful for modeling reaction dynamics of strongly correlated systems. In an alternative strategy, strongly correlated systems are studied using multireference methods such as CASPT2,61 where the accuracy depends in part on active space and basis set.62 This method is quite expensive for large systems with the affordable limit often being 16 active electrons in 16 active orbitals with conventional solvers. However, CASPT2 is affordable for the present small reaction where we used an active space that has 11 electrons in 9 orbitals and is composed of four nominally doubly occupied orbitals [σ(C-H), σ(O-O), and two π(O-O)], the

p(CO) orbital that contains the unpaired electron, and four nominally unoccupied orbitals [σ*(C-H), σ*(O-O), and two π*(O-O)]. Therefore we can test it for the present reaction. Table 1 compares the energy of CASPT2(11,9)/aug-cc-pVTZ to our best estimate and shows high accuracy for the classical barrier height but poor accuracy when zero point energies are included. Tables S7 and S9-S10 show that latter is caused mainly by an inaccurate zero-point vibrational energy of CH3O in the CASPT2 calculation. We also studied the strategy of using CASPT2 geometries and frequencies for coupled cluster calculations. Again the results are of poor quality due to the inaccurate zero-point vibrational energy of CH3O. In particular, Table 1 shows a remarkable difference of 0.81 kcal/mol in enthalpy of activation at 0 K between W3X-L calculations with the RCCSD(T)-F12a/cc-pVTZ-F12 geometry and frequencies and those with the CASPT2 geometry and frequencies. Thus, the strategy if using multireference perturbation theory is not accurate in the present case.

V. RATE CONSTANTS A dual-level strategy, as used in our previous investigations,32,53,63 is employed to account for the difference between the benchmark barrier height and the density functional barrier height. First, the rate constant is calculated by conventional transition state theory (TST) without tunneling using GMM(P). Then MN15-L/MG3S was utilized to do direct dynamics calculations of the reaction rate using canonical variational transition state theory with small-curvature TST CVT/SCT tunneling.64,65 These rate constants are labeled respectively as kHL and kLL , where HL

and LL denote higher level and lower level, respectively. We can factor the latter rate constant

ACS Paragon Plus Environment

Page 9 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 9

as CVT/SCT CVT SCT TST SCT kLL = kLL κ LL = kLL ΓLLκ LL

(7)

CVT TST where kLL and kLL are the CVT and conventional TST rate constants (both without SCT is the tunneling) at level LL, ΓLL is the recrossing transmission coefficient, and κ LL

tunneling transmission coefficient. Then the final rate constant is calculated as TST SCT CVT/SCT ‡ ‡ k = kHL ΓLLκ LL = exp ⎡⎢− ΔGHL − ΔGLL / RT ⎤⎥ kLL ⎣ ⎦

(

)

(8)

The frequencies needed for zero point energies and partition coefficients are always calculated with the method used to optimize the geometry; thus they are calculated by the jun method for HL and by MN15-L/MG3S for LL. Then scale factors (Table S11) were used to scale all the calculated frequencies, by a method explained previously.38 The calculated rate constants are listed in Table 2 and Table S12 in the temperature range from 250 to 2000 K. The theoretical rate constant is fitted to a functional form proposed previously,66 and this yields 𝑘 = 1.22 ×10!!! cm3 molecule−1 s−1 (

!! !! !.!" ) exp !""

[−

[!.!" !"#$/!"#](!! !! ) !(! ! !!!! )

]

(9)

where R is the gas constant, T is temperature, and T0 is 434 K. The fit reproduces the calculations within 5%. The experimental rate constants from several groups over the temperature range from 250 K to 610 K have been fit to16 k = 3.81×10-21 T 2.40 exp(-208/T) cm3 molecule-1 s-1

250-610 K

(11)

and we compare to this data over the range 250 K to 450 K. We do not compare to experiment above 450 K because Wantuck et al.14 concluded that the measured rate constant above 450 K might contain significant contributions from other processes (collision-induced isomerization and collision-induced decomposition) as well as abstraction. Equation 11 agrees with the recommendation of an earlier review1 at temperatures 298-450 K within 10% (results below 298 K were not available at the time of the earlier review), and the various experimental data agree with each other within a factor of about 1½–2 over this temperature range. The final calculated rate constants in Table 2 agree with experimental results fitted by eq 11 within 34% over the 250-450 K temperature. Considering that quadruple and higher

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 20 10

excitations affect the rate constant in this temperature range by a factor of 47 at 450 K and a factor of 1.0×103 at 250 K, the quantitative agreement with experiment is remarkable (and perhaps even fortuitously better than it ought to be), but the main take-home message is that the predicted rate constant would be orders of magnitude off from experiment if one stopped with CCSDT.

VI. ACTIVATION ENERGIES The temperature-dependent activation energies in the CH3O + O2 reaction were computed from the fit as67 Ea = −R

d ln k d(1/ T )

(10)

that is, Ea is local slope of the Arrhenius plot and would be a constant for a linear Arrhenius plot. It has been known for a long time68,69 that deviations for constant Ea can be large when one considers a wide temperature range (that is, Arrhenius plots are far from linear), but the results in Table 2 are a particularly striking example, with Ea increasing from 0.58 to 9.86 kcal/mol. The dependence of Ea on temperature is large even above 800 K, where the variational effects and tunneling are negligible.

VII. CONCLUDING REMARKS Obtaining quantitative rate constants by theoretical methods requires obtaining quantitative energy barriers, which are determined by three factors: optimized geometries, zero-point vibrational energies, and single-point energies. Here, we show how to design computational strategies and methods to obtain quantitative energy barriers for strongly correlated systems in a practical way by electronic structure methods. We investigated the CH3O + O2 reaction acted as a prototype reaction with a strongly correlated electronic structure, and we calculated quantitative rate constants. We used coupled cluster theory, where the accuracy of the method depends on excitation level and basis set. We showed that RCCSD(T)-F12a/jun-cc-pVTZ can yield a reliable geometrical optimization and zero-point vibrational energy, whereas CASPT2(11,9)/aug-cc-pVTZ cannot provide reliable zero-point vibrational energies. We included connected excitations up to pentuple, and we used two choices of basis set to extrapolate key components to the complete basis set limit to show whether each coupled cluster component in the scheme is well converged to the complete basis

ACS Paragon Plus Environment

Page 11 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 11

set limit. The combination of connected excitations up to pentuple with complete basis set extrapolation is judged to give a best estimate of the barrier height within a few tenths of a kcal/mol of the frozen-core complete configuration interaction limit. We then found that Kohn-Sham density functional theory with the MN15-L exchange-correlation functional has remarkable accuracy, which is much higher than that of CCSD(T)/CBS. These findings in this work are broadly relevant to kinetics processes involves strongly correlated states, which includes many processes in atmospheric and combustion chemistry.

! ASSOCIATED CONTENT SSupporting Information The Supporting Information is available free of charge on the ACS Publications website DOI: 10.1021/

at

References for basis sets and number of contracted basis functions for transition state calculations; computational details of orbitals and coupled cluster calculations, W2X, MW2-F12, MW3X-L, and W3X-L composite methods, electronic partition functions, symmetry numbers, and spin-orbit coupling; scale factors for the vibrational frequencies; activation enthalpies at 0 K and barrier heights for TS1 and TS2; components of barrier height increments (kcal/mol) for TS1; selected internuclear distances and bond angles of TS1; frequencies and zero-point vibrational energies of CH3O, O2, and TS1; rate constants, transmission coefficients, and activation energies; and coordinates and absolute energies of optimized structures. (PDF ) ! AUTHOR INFORMATION Corresponding Authors *E-mail: [email protected], , [email protected] Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. ORCID Bo Long: 0000-0002-9358-2585 Junwei Lucas Bao: 0000-0002-4967-663X

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 20 12

Donald G. Truhlar: 0000-0002-7742-7294 Notes The authors declare no competing financial interest.

! ACKNOWLEDGMENT

This work was supported in part by the National Natural Science Foundation of China (41775125), the Science and Technology Foundation of Guizhou Province, China, ([2018]1080), the Science and Technology Foundation of Guizhou Provincial Department of Education, China ([2015]350), and the U.S. Department of Energy (DE-SC0015997). Computations were performed using resources of Minnesota Supercomputing Institute and the National Energy Research Scientific Computing Centre.

ACS Paragon Plus Environment

Page 13 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 13

Table 1. The enthalpy of activation (Δ𝐻!‡ , 0 Κ) and classical barrier height (𝑉 ‡ ) of the transition state TS1 (including spin-orbit coupling) and the mean unsigned error (MUE) of these quantities (in kcal/mol).a Method Geometry 𝑉‡ MUE Δ𝐻!‡ b GMM(P) (best estimate) jun 2.09 1.79 0.00 GMMQ ” 2.40 2.10 0.31 GMM(Q) ” 1.82 1.52 0.27 MW2-F12 ” 4.83 4.53 2.74 jun ” 4.21 3.91 2.12 c CASPT2 Consistent 1.41 1.77 0.35 Weizmann and modified Weizmann methods: W3X-L RCCSD(T)-F12a/cc-pVTZ-F12 ” jun ” RCCSD(T)-F12a/cc-pVDZ-F12 ” RQCISD/MG3S ” CASPT2 MW3X-L jun W2X ” ” RCCSD(T)-F12a/cc-pVDZ-F12

2.62 2.57 2.37 2.32 1.81 3.08 4.32 4.10

2.29 2.27 2.29 1.58 2.17 2.78 4.02 4.02

0.52 0.48 0.39 0.22 0.33 0.99 2.23 2.12

Density functional theory: MN15-L/maug-cc-pVTZ MN15-L/MG3S

3.31 3.19

2.45 2.36

0.94 0.84

Consistent Consistent

a

All reactant energies used to compute barrier heights and enthalpies include a spin-orbit contribution of –0.09 kcal/mol, and this raised all reported barrier heights and enthalpies of activation by 0.09 kcal/mol. b As explained in the text, jun denotes RCCSD(T)-F12a/jun-cc-pVTZ. c “Consistent” denotes that the energy and the geometry optimization were performed by the same method.

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 20 14

Table 2. The calculated rate constant, factors in the rate constant, and activation energies. from !"! a T k(theory)a k(experiment)b 𝐸! (theory)c 𝑘!" к!"# !! Г !! (cm3 molecule−1 s−1) (cm3 molecule−1 s−1) (cm3 molecule−1 s−1) (K) (kcal/mol) −16 −15 −16 7.56 × 10 1.31 × 10 9.4 × 10 250 3.28 0.90 0.58 −15 −15 −15 0.85 1.42 × 10 1.65 × 10 1.6 × 10 298 2.17 0.93 −15 −15 −15 0.87 1.45 × 10 1.67 × 10 1.7 × 10 300 2.14 0.93 −15 −15 −15 1.99 × 10 1.97 × 10 2.2 × 10 330 1,84 0.95 1.06 −15 −15 −15 3.60 × 10 2.89 × 10 4.0 × 10 400 1.48 0.97 1.53 −15 −15 −15 5.07 × 10 3.71 × 10 5.6 × 10 450 1.36 0.98 1.88 −14 −15 1.12 × 10 7.10 × 10 600 1.18 1.00 2.87 −14 −14 2.44 × 10 1.41 × 10 800 1.10 1.00 4.05 −14 −14 4.45 × 10 2.46 × 10 1000 1.06 1.00 5.09 −14 −14 7.28 × 10 3.92 × 10 1200 1.04 1.00 6.08 1.34 × 10−13 7.02 × 10−14 1500 1.03 0.99 7.51 −13 −14 1.88 × 10 9.78 × 10 1700 1.02 0.99 8.45 −13 −13 2.93 × 10 1.51 × 10 2000 1.01 0.99 9.86 a from eq 8 b from eq 11 c from eq 10

ACS Paragon Plus Environment

Page 15 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 15

Figure 1. Optimized geometries of the transition structures for the reaction of CH3O with O2, as calculated by RCCSD(T)-F12a/jun-cc-pVTZ.

! REFERENCES

1 Orlando, J. J.; Tyndall, G. S.; Wallington, T. J. The Atmospheric Chemistry of Alkoxy Radicals. Chem. Rev. 2003, 103, 4657-4690. 2 Atkinson, R. Rate Constants For The Atmospheric Reactions of Alkoxy Radicals: An Updated Estimation Method. Atmos. Environ. 2007, 41, 8468-8485. 3 Atkinson, R. Atmospheric Reactions of Alkoxy and Β-Hydroxyalkoxy Radicals. Int. J. Chem. Kinet. 1997, 29, 99-111. 4 Jenkin, M. E.; Clemitshaw, K. C. Ozone and Other Secondary Photochemical Pollutants: Chemical Processes Governing Their Formation in The Planetary Boundary Layer. Atmos. Environ. 2000, 34, 2499-2527. 5 Lim, Y. B.; Ziemann, P. J. Products and Mechanism of Secondary Organic Aerosol Formation from Reactions of n-Alkanes with OH Radicals in the Presence of NOx. Environ. Sci. Technol. 2005, 39, 9229-9236. 6 Davis, A. C.; Francisco, J. S. Reactivity Trends within Alkoxy Radical Reactions Responsible for Chain Branching. J. Am. Chem. Soc. 2011, 133, 18208-18219. 7 Carrasquillo, A. J.; Hunter, J. F.; Daumit, K. E.; Kroll, J. H. Secondary Organic Aerosol Formation via the Isolation of Individual Reactive Intermediates: Role of Alkoxy Radical Structure. J. Phys. Chem. A 2014, 118, 8807-8816. 8 Dibble, T. S.; Chai, J. Critical Review of Atmospheric Chemistry of Alkoxy Radicals. In Advances in Atmospheric Chemistry; World Scientific: 2017, vol. 1, pp 185-269. 9 Dibble, T. S. Comment on ”Isomerization of the Methoxy Radical Revisited: The Impact of Water Dimers” by Kumar, P.; Biswas, P.; Bandyopadhyay, B. Phys. Chem. Chem. Phys., 2016, 18, 27728 and “Isomerization of methoxy radical in the troposphere: competition between acidic, neutral and basic catalysts” by Sarkar, S.; Mallick, S.; Deepak; Kumar, P.;

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 20 16

Bandyopadhyay, B . Phys. Chem. Chem. Phys., 2017, 19, 278. Phys. Chem. Chem. Phys. 2018, 20 , 11481-11482. 10 Wiebe, H.; Villa, A.; Hellman, T.; Heicklen, J. Photolysis of Methyl Nitrite. J. Am. Chem. Soc. 1973, 95, 7−13. 11 Cox, R.; Derwent, R.; Kearsey, S.; Batt, L.; Patrick, K. Photolysis of Methyl Nitrite: Kinetics of the Reactions of the Methoxy Radical with O2. J. Photochem. 1980, 13, 149−163. 12 Gutman, D.; Sanders, N.; Butler, J. E. Kinetics of the Reactions of Methoxy and Ethoxy Radicals With Oxygen. J. Phys. Chem. 1982, 86, 66-70. 13 Lorenz, K.; Rhäsa, D.; Zellner, R.; Fritz, B. Laser Photolysis-LIF Kinetic Studies of the Reactions of CH3O and CH2CHO with O2 Between 300 and 500 K. Ber. Bunsen-Ges. Phys. Chem. 1985, 89, 341−342. 14 Wantuck, P. J.; Oldenborg, R. C.; Baughcum, S. L.; Winn, K. R. Removal Rate Constant Measurements For Methoxy Radical By Oxygen over the 298-973 K Range. J. Phys. Chem. 1987, 91, 4653-4655. 15 Atkinson, R.; Baulch, D. L.; Cox, R. A.; Crowley, J. N.; Hampson, R. F.; Hynes, R. G.; Jenkin, M. E.; Rossi, M. J.; Troe, J. Subcommittee, I. Evaluated Kinetic and Photochemical Data for Atmospheric Chemistry: Volume II – Gas Phase Reactions of Organic Species. Atmos. Chem. Phys. 2006, 6, 3625-4055. 16 Chai, J.; Hu, H.; Dibble, T. S.; Tyndall, G. S.; Orlando, J. J. Rate Constants and Kinetic Isotope Effects for Methoxy Radical Reacting with NO2 and O2. J. Phys. Chem. A 2014, 118, 3552-3563. 17 Jungkamp, T. P. W.; Seinfeld, J. H. The Enthalpy of Formation of Trioxy Radicals ROOO (R = H, CH3, C2H5). An ab initio study. Chem. Phys. Lett. 1996, 257, 15-22. 18 Jungkamp, T. P. W.; Seinfeld, J. H. The Mechanism of Methoxy Radical Oxidation: Hydrogen Abstraction Versus Trioxy Radical Formation. Chem. Phys. Lett. 1996, 263, 371-378. 19 Bofill, J. M.; Olivella, S.; Solé, A.; Anglada, J. M., The Mechanism of Methoxy Radical Oxidation by O2 in the Gas Phase. Computational Evidence for Direct H Atom Transfer Assisted by an Intermolecular Noncovalent O···O Bonding Interaction. J. Am. Chem. Soc. 1999, 121, 1337-1347. 20 Setokuchi, O.; Sato, M., Direct Dynamics of an Alkoxy Radical (CH3O, C2H5O, and i-C3H7O) Reaction with an Oxygen Molecule. J. Phys. Chem. A 2002, 106, 8124-8132. 21 Hu, H.; Dibble, T. S. Quantum Chemistry, Reaction Kinetics, and Tunneling Effects in the Reaction of Methoxy Radicals with O2. J. Phys. Chem. A 2013, 117, 14230-14242. 22 Chan, B.; Radom, L. W2X and W3X-L: Cost-Effective Approximations to W2 and W4 with kJ mol–1 Accuracy. J. Chem. Theory. Comput. 2015, 11, 2109-2119. 23 Orlando, J. J.; Tyndall, G. S. Laboratory Studies of Organic Peroxy Radical Chemistry: An Overview With Emphasis on Recent Issues of Atmospheric Significance. Chem. Soc. Rev. 2012, 41, 6294-6317. 24 Winiberg, F. A. F.; Dillon, T. J.; Orr, S. C.; Groß, C. B. M.; Bejan, I.; Brumby, C. A.; Evans, M. J.; Smith, S. C.; Heard, D. E.; Seakins, P. W. Direct Measurements of OH and Other Product Yields from the HO2  + CH3C(O)O2 Reaction. Atmos. Chem. Phys. 2016, 16 (6), 4023-4042. 25 Teng, A. P.; Crounse, J. D.; Wennberg, P. O. Isoprene Peroxy Radical Dynamics. J. Am. Chem. Soc. 2017, 139, 5367-5377. 26 Fittschen, C.; Whalley, L. K.; Heard, D. E. The Reaction of CH3O2 Radicals with OH

ACS Paragon Plus Environment

Page 17 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 17

Radicals: A Neglected Sink for CH3O2 in the Remote Atmosphere. Environ. Sci. Technol. 2014, 48, 7700-7701. 27 Caravan, R. L.; Khan, M. A. H.; Zádor, J.; Sheps, L.; Antonov, I. O.; Rotavera, B.; Ramasesha, K.; Au, K.; Chen, M.-W.; Rösch, D.; Osborn, D. L.; Fittschen, C.; Schoemaecker, C.; Duncianu, M.; Grira, A.; Dusanter, S.; Tomas, A.; Percival, C. J.; Shallcross, D. E.; Taatjes, C. A. The Reaction of Hydroxyl and Methylperoxy Radicals Is Not a Major Source of Atmospheric Methanol. Nature Communications 2018, 9, 4343. 28 Vereecken, L.; Glowacki, D. R.; Pilling, M. J. Theoretical Chemical Kinetics in Tropospheric Chemistry: Methodologies and Applications. Chem. Rev. 2015, 115, 4063-4114. 29 Adler, T. B.; Knizia, G.; Werner, H.-J. A Simple and Efficient CCSD(T)-F12 Approximation. J. Chem. Phys. 2007, 127, 221106. 30 Knizia, G.; Adler, T. B.; Werner, H.-J. Simplified CCSD(T)-F12 methods: Theory and Benchmarks. J. Chem. Phys. 2009, 130, 054104. 31 Sylvetsky, N.; Peterson, K. A.; Karton, A.; Martin, J. M. L. Toward a W4-F12 approach: Can Explicitly Correlated and Orbital-Based Ab Initio CCSD(T) Limits Be Reconciled? J. Chem. Phys. 2016, 144, 214101. 32 Long, B.; Bao, J. L.; Truhlar, D. G. Unimolecular Reaction of Acetone Oxide and its Reaction with Water in the Atmosphere. Proc. Natl. Acad. Sci. USA 2018, 115, 6135-6140. 33 Helgaker, T.; Klopper, W. Koch, H.; Nago, J. Basis-Set Convergence of Correlated Calculations on Water. J. Chem. Phys. 1997, 106, 9639-9646. 34 Foster, S. C.; Misra, P.; Lin, T. Y. D.; Damo, C. P.; Carter, C. C.; Miller, T. A. Free Jet-Cooled Laser-Induced Fluorescence Spectrum of Methoxy. 1. Vibronic Analysis of the A! and X! States. J. Phys. Chem. 1988, 92, 5914-5921. 35 Brossard, S. D.; Carrick, P. G.; Chappell, E. L.; Hulegaard, S. C.; Engelking, P. C. The 2 A1–2E Electronic Spectrum of Methoxy, CH3O: Vibrational Progressions, and Spin-Orbit and Jahn–Teller Couplings. J. Chem. Phys. 1986, 84, 2459-2465. 36 Misra, P.; Zhu, X.; Hsueh, C.-Y.; Halpern, J. B. Laser Excitation and Emission Spectroscopy of the Methoxy Radical in a Supersonic Jet. Chem. Phys. 1993, 178, 377-385. 37 Lee, Y. Y.; Wann, G. H.; Lee, Y. P. Vibronic Analysis of the Ã→X̃ Laser-Induced Fluorescence of Jet-Cooled Methoxy (CH3O) Radical. J. Chem. Phys. 1993, 99, 9465-9471. 38 Vaucher, A. C.; Reiher, M. Steering Orbital Optimization out of Local Minima and Saddle Points Toward Lower Energy. J. Chem. Theory. Comput. 2017, 13, 1219-1228. 39 Alecu, I. M.; Zheng, J.; Zhao, Y.; Truhlar, D. G. Computational Thermochemistry: Scale Factor Databases and Scale Factors for Vibrational Frequencies Obtained from Electronic Model Chemistries. J. Chem. Theory. Comput. 2010, 6, 2872-2887. 40 Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.;

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 20 18

Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian 16, Revision B.01, Gaussian, Inc., Wallingford CT, 2016. 41 Zhao, Y.; Peverati, R.; Tang, K.; Luo, S.; Yu, H. S.; He, X.; Truhlar, D. G. MN-GFM 6.7; Department of Chemistry, University of Minnesota: Minneapolis, MN 55455-0431, 2015. 42 Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G.; Shamasundar, K. R.; Adler, T. B.; Amos, R. D.; Bernhardsson, A.; Berning, A.; Cooper, D. L.; Deegan, M. J. O.; Dobbyn, A. J.; Eckert, F.; Goll, E.; Hampel, C.; Hesselmann, A.; Hetzer, G.; Hrenar, T.; Jansen, G.; Köppl, C.; and, Y. L.; Lloyd, A. W.; Mata, R. A.; May, A. J.; McNicholas, S. J.; Meyer, W.; Mura, M. E.; Nicklass, A.; O'Neill, D. P.; Palmieri, P.; Peng, D.; Pflüger, K.; Pitzer, R.; Reiher, M.; Shiozaki, T.; Stoll, H.; Stone, A. J.; Tarroni, R.; Thorsteinsson, T.; Wang, M. Molpro, version 2015.1, a package of ab initio programs; 43 M. Kállay, Z. Rolik, J. Csontos, I. Ladjánszki, L. Szegedy, B. Ladóczki, and G. Samu. See also Z. Rolik, L. Szegedy, I. Ladjánszki, B. Ladóczki, and M. Kállay. An Efficient Linear-Scaling Ccsd(T) Method Based on Local Natural Orbitals. J. Chem. Phys. 2013, 139, 094105. 44 Kállay, M.; Rolik, Z.; Ladjánszki, I.; Szegedy, L.; Ladóczki, B.; Csontos, J.; Kornis, B. MRCC, A Quantum Chemical Program Suite, http://www.mrcc.hu/. 45 Zheng, J.; Bao, L. J. W.; Meana-Pañeda, R.; Zhang, S.; Lynch, B. J.; Corchado, J. C.; Chuang, Y.-Y.; Fast, P. L.; Hu, W.-P.; Liu, Y.-P.; Lynch, G. C.; Nguyen, K. A.;Jackels, C. F.; Fernandez-Ramos, A.; Ellingson, B. A.; Melissas, V. S.;Villa, J.; Rossi, I.; Coitino, L.; Pu, J.; Albu, T. V.; Steckler, R.; Garrett, B. C.; Issacson, A. D.; Truhlar, D. G. Polyrate – version 2017-C; University of Minnesota: Minneapolis, 2018. 46 Zheng, J.; Bao, L. J. W.; Zhang, S.; Corchado, J. C.; Chuang, Y.-Y.; Coitiño, E. L.; Ellingson, B. A.; Truhlar, D. G. Gaussrate – version 2017-B; University of Minnesota: Minneapolis, 2018. 47 Karton, A. A computational chemist's guide to accurate thermochemistry for organic molecules. WIREs Comput. Mol. Sci. 2016, 6, 292–310. 48 Marenich, A. V.; Boggs, J. E. The molecular structure, spin–vibronic energy levels, and thermochemistry of CH3O. J. Mol. Struct. 2006, 780-781, 163-170. 49 Marenich, A. V.; Boggs, J. E. A model spin-vibronic hamiltonian for twofold degenerate electron systems: a variational ab initio study of X̃  2E CH3O•. J. Chem. Phys. 2004, 122, 024308. 50 Meana-Pañeda, R.; Xu, X.; Ma, H.; Truhlar, D. G. Computational Kinetics by Variational Transition-State Theory with Semiclassical Multidimensional Tunneling: Direct Dynamics Rate Constants for the Abstraction of H from CH3OH by Triplet Oxygen Atoms. J. Phys. Chem. A 2017, 121 , 1693-1707. 51 Jacox, M. E. Vibrational and Electronic Energy Levels of Polyatomic Transient Molecules, P.J. Linstrom and W.G. Mallard, in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. National Institute of Standards and Technology, Gaithersburg MD, 2018. 52 Long, B.; Bao, J. L.; Truhlar, D. G. Reaction of SO2 with OH in the Atmosphere. Phys. Chem. Chem. Phys. 2017, 19, 8091-8100. 53 Wei, M.-L.; Tan, X.-F.; Long, Z.-W.; Long, B. Atmospheric Chemistry of CH3O: Its Unimolecular Reaction and Reactions with H2O, NH3, and HF. RSC Adv. 2017, 7, 56211-56219.

ACS Paragon Plus Environment

Page 19 of 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society 19

54 Long, B.; Bao, J. L.; Truhlar, D. G. Atmospheric Chemistry of Criegee Intermediates: Unimolecular Reactions and Reactions with Water. J. Am. Chem. Soc. 2016, 138, 14409-14422. 55 Misiewicz, J. P.; Elliott, S. N.; Moore, K. B.; Schaefer, H. F. Re-examining Ammonia Addition to the Criegee Intermediate: Converging to Chemical Accuracy. Phys. Chem. Chem. Phys. 2018, 20, 7479-7491. 56 Fang, Y.; Liu, F.; Barber, V. P.; Klippenstein, S. J.; McCoy, A. B.; Lester, M. I. Communication: Real Time Observation of Unimolecular Decay of Criegee Intermediates to OH Radical Products. J. Chem. Phys. 2016, 144, 061102. 57 Mallick, S.; Kumar, P. Impact of Beyond-CCSD(T) Corrections on Reaction Energetics and Rate Constants of the OH• + HCl Reaction. J. Phys. Chem. A 2018, 122, 7151-7159. 58 Yu, H. S.; He, X.; Truhlar, D. G. MN15-L: A New Local Exchange-Correlation Functional for Kohn–Sham Density Functional Theory with Broad Accuracy for Atoms, Molecules, and Solids. J. Chem. Theory Comput. 2016, 12, 1280-1293. 59 Yu, H. S.; Li, S. L.; Truhlar, D. G. Perspective: Kohn-Sham Density Functional Theory Descending a Staircase. J. Chem. Phys. 2016, 145, 130901. 60 Wang, Y.; Jin, X.; Yu, H. S.; Truhlar, D. G.; He, X. Revised M06-L Functional for Improved Accuracy on Chemical Reaction Barrier Heights, Noncovalent Interactions, and Solid-State Physics. Proc. Nat. Acad. Sci. USA 2017, 114, 8487-8492. 61 Andersson, K.; Malmqvist, P. Å.; Roos, B. O. Second-Order Perturbation Theory with a Complete Active Space Self‐Consistent Field Reference Function. J. Chem. Phys. 1992, 96, 1218-1226. 62 Tishchenko, O.; Zheng, J.; Truhlar, D. G. Multireference Model Chemistries for Thermochemical Kinetics. J. Chem. Theory Comput. 2008, 4, 1208-1219. 63 Tan, X.-F.; Long, B.; Ren, D.-S.; Zhang, W.-J.; Long, Z.-W.; Mitchell, E. Atmospheric Chemistry of CH3CHO: The Hydrolysis of CH3CHO Catalyzed by H2SO4. Phys. Chem. Chem. Phys. 2018, 20, 7701-7709. 64 Liu, Y. P.; Lynch, G. C.; Truong, T. N.; Lu, D. H.; Truhlar, D. G.; Garrett, B. C. Molecular Modeling of The Kinetic Isotope Effect for the [1,5]-Sigmatropic Rearrangement of cis-1,3-Pentadiene. J. Am. Chem. Soc. 1993, 115, 2408-2415. 65 Bao, J. L.; Truhlar, D. G. Variational Transition State Theory: Theoretical Framework and Recent Developments. Chem. Soc. Rev. 2017, 46, 7548-7596. 66 Zheng, J.; Truhlar, D. G. Multi-Path Variational Transition State Theory For Chemical Reaction Rates of Complex Polyatomic Species: Ethanol + OH Reactions. Faraday Discuss. 2012, 157, 59-88. 67 Truhlar, D. G. Interpretation of the Activation Energy. J. Chem. Educ. 1978, 55, 309. 68 Blais, N.C.; Truhlar, D. G.; Garrett, B. C. Dynamical Calculation of the Temperature Dependence of the Activation Energy for a Chemical Reaction from 444 K to 2400 K. J. Phys. Chem. 1981, 85, 1094-1096. 69 Garrett, B. C.; Truhlar, D. G.; Bowman, J. M.; Wagner, A. F.; Robie, D.; Arepalli, S.; Presser, N.; Gordon, R. J. Ab initio Predictions and Experimental Confirmation of Large Tunneling Contributions to Rate Constants and Kinetic Isotope Effects For Hydrogen-Atom Transfer Reactions. J. Am. Chem. Soc. 1986, 108, 3515-3516.

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 20 20

TOC Graphic

ACS Paragon Plus Environment