H-Atom Abstraction Reactions by Ground-State Ozone from Saturated

J. Würmel† and J. M. Simmie‡. † Galway Mayo Institute of Technology, Galway H91 T8NW, Ireland. ‡ School of Chemistry, National University of ...
0 downloads 10 Views 1MB Size
Subscriber access provided by University of Massachusetts Amherst Libraries

Article

H-Atom Abstraction Reactions by Ground State Ozone from Saturated Oxygenates Judith Wurmel, and John M Simmie J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b07760 • Publication Date (Web): 04 Oct 2017 Downloaded from http://pubs.acs.org on October 4, 2017

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry A 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 34

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

The Journal of Physical Chemistry

H-atom Abstraction Reactions by Ground State Ozone from Saturated Oxygenates J. Würmel∗,† and J. M. Simmie‡ †Galway Mayo Institute of Technology, Galway H91 T8NW, Ireland ‡School of Chemistry, National University of Ireland, Galway H91 TK33, Ireland E-mail: [email protected] Phone: +353-91-742391 Abstract Theoretical insights into H-abstraction by ozone from saturated species have been virtually nonexistent which is in sharp contrast to reactions with unsaturated species. Our computed rate constants at various levels of theory for reaction with tetrahydrofuran and its methyl derivatives shows that abstraction occurs primarily at the carbons situated beside the heterocyclic oxygen with least likely reaction from the methyl groups. All the methods tested are in broad relative agreement with this conclusion and with recent experiments although they do differ widely as regards their absolute values. Abstraction leads variously to the formation of hydrotrioxides, ROOOH, or more directly to R – OH + O2 (1 ∆g ). To understand some of the observed behaviour dimethyl and diethyl acyclic ethers, methyl and ethyl hydroperoxides and the alcohols methanol and ethanol were also included in this study. For the simplest system, CH3 OH + O3 −−→ H2 C – O + H2 O2 , an all-electron scalar relativistic CCSD(T)/CBS approach yields a barrier height which differs by > 5 kJ mol−1 when an additional multi-reference treatment up to CCSDT(Q) is considered. Draft: October 2, 2017

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Introduction There are ongoing efforts to ‘green’ manufacturing processes and one strand is the utilisation of green solvents, that is, a solvent manufactured from renewable resources with a low life cycle footprint. 1,2 In particular 2-methyltetrahydrofuran has found promise as an alternative solvent suitable for environmentally benign syntheses because of its remarkable stability, low miscibility with water, etc. 3 Concurrently there is a pressing need to produce high-energy density liquid fuels from lignocellulosic biomass, the so-called ‘next generation biofuels’, to reduce our dependency on fossil fuels and to lessen their negative environmental impact. 4–6 Whilst the disposal, accidental or otherwise, of such materials may impact on the water table — here we are concerned with emissions to the atmosphere and their fate thereof. In particular we focus on four oxolanes, more commonly known as tetrahydrofurans, oxolane (THF), both 2-methyl and 3-methyl oxolanes (2MTHF, 3MTHF) and 2,5-dimethyloxolane (25DMTHF). These are classified as heterocyclic compounds and more specifically as members of the cyclic ether family. The rates of reactions with the hydroxyl radical, for example, have been extensively studied. 7 We are concerned here with their reactivity towards ozone O3 (X1 A1 ) and although ozone is present in the troposphere the highest concentrations are to be found in heavily polluted urban areas where solvent spills might reasonably be expected to occur. Recent experimental determinations of the total rates of reaction of THF, 2MTHF and 25DMTHF with ozone in a photochemical reactor supply the essential verification that is needed. 8 Andersen and colleagues used both relative and absolute rate methods at 298 ± 3 K and report that the relative rates are as 1:2.9:7.1 for THF:2MTHF:25DMTHF. These are slow hydrogen-atom abstraction reactions in comparison say to reaction even with an unsaturated cyclic ether such as furan. Mechanistically, reaction with an alkene moiety C – C proceeds via a fast 1,3-dipolar addition reaction — an option not available to saturated compounds. In fact, in the over 1,800 records in a comprehensive kinetics database 9 there are very few entries featuring neutral saturated molecules, M + O3 −−→ products. 2 ACS Paragon Plus Environment

Page 2 of 34

Page 3 of 34

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

The Journal of Physical Chemistry

Cerkovnik et al. have shown that ozone inserts into the aldehydic C–H bond in benzaldehyde, C6 H5 CHO, leading to the formation of an organic hydrotrioxide which can then undergo a series of conformational changes ending up with the elimination of singlet oxygen, O2 (1 ∆g ), and benzoic acid. 10 Their experiments were backed up by B3LYP/6-31G(d,p) calculations which determined an activation energy of 46 kJ mol−1 and ∆r H of −240 kJ mol−1 , for the highly exothermic formation of the benzoyl hydrotrioxide, C6 H5 C(O)OOOH. More generally there are a number of studies of the mechanisms and kinetics of ozone reacting with unsaturated oxygenates, 11–13 acetaldehyde, 14 fluorocarbons, 15 methane, 16 formaldehyde, 17,18 unsaturated ketones, 19 alkyl hydroperoxides 20,21 and summaries in a recent review of atmospheric reaction mechanisms in the troposphere. 22 However, although there have been these experimental studies, the intrinsic multireference character of ozone 23 has probably accounted for the dearth of theoretical calculations for the reactions of systems larger than H2 CO + O3 . 17 On the application side the enhancement of combustion by ozone is of interest as, for example, seen in studies of the increase in burning velocities of methane/air flames and similar. 24–28 The aims and objectives of this study are therefore to determine the reactivity of the various sites in the cyclic ethers and thereby validate the experimental measurements, to contrast cyclic and acyclic ethers, and to understand how the reactivity of cognate hydroperoxides and alcohols compares.

Computational methodology The bulk of our computations used a recent Minnesota density functional, RMN12SX, 29 a screened-exchange hybrid functional coupled with a split-valence triple-ζ basis set with added diffuse and polarization functions, 6-311++G(d,p). Thus, zero-point corrected electronic energies, vibrational frequencies (scaled by 0.9755), moments of inertia and relaxed potential

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

energy scans were carried out with this model chemistry on a singlet surface with Gaussian16. 30 Wang and co-workers 17 explored both the singlet and triplet surfaces in H2 CO + O3 −−→ HCO• + HOOO• and found that the singlet barrier is considerably lower than that of the triplet. The H-atoms occupy non-equivalent positions on the THF ring and they can be labelled as axial and equatorial as shown in Fig 1 for just those on the vicinal carbon atoms to the heterocyclic oxygen.

Figure 1: Axial and equatorial sites on THF. Intrinsic reaction coordinate calculations were used to verify the connection between reactants, transition state and products. 31,32 The keyword Recorrect=Never was used to check that the end points of the IRC corresponded to the actual structures because of the variability of the products for apparently similar reactions. Very weakly-bonded prereaction complexes were found for which basis set superposition error corrections 33 were applied; given their very weak nature they are kinetically insignificant. The simplest transition states were checked in CCSD(T)-F12A/VDZ-f12//M06-2X/6311++G(d,p) molpro calculations 34–36 which showed T1 diagnostics well above the recommended limit, 37 for example, T1=0.045 for the C3 axial H-abstraction transition state from 3MTHF at CCSD(T)-F12a/VDZ-F12 — thereby necessitating multireference methods which are however unaffordable with the resources available. The sole exception to this limitation was for the small system CH3 OH + O3 −−→ CH2 O + HOOOH for which T1 > 0.06 and for which it was possible to employ much higher level composite methods, W2X and W3X-L. The former provides an accurate approximation to the all-electron scalarrelativistic CCSD(T)/CBS energy whilst the latter provides additionally a multi-reference 4 ACS Paragon Plus Environment

Page 4 of 34

Page 5 of 34

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

The Journal of Physical Chemistry

post-CCSD(T) treatment up to CCSDT(Q). 38 These model chemistries have been used previously, employing the applications molpro 34 and mrcc, 39 to compute formation enthalpies via an atomisation procedure, validated against the most credible database of formation enthalpies known, the Active Thermochemical Tables, 40 and shown to give excellent results. 41 The results of these calculations were then used to determine thermochemical parameters such as entropy, S(298.15K) and the specific heat at constant pressure, Cp (298.15K), and hence the kinetics of reaction via canonical transition state theory as framed in the thermo module of the application MultiWell. 42 Thermo includes a correction for unsymmetrical Eckart tunnelling and treats 1-dimensional hindered rotors by fitting Fourier series to both the potential energy, and the rotational constant, as a function of dihedral angle. 43 In the case of highly symmetric methyl rotors, determining the barrier height is sufficient to compute the contribution that these modes make with the more elaborate treatment being unnecessary.

Pseudo-rotational considerations Tetrahydrofuran, in common with most n-ring systems with n > 3, exhibits the phenomenon of pseudo-rotation with a recent study on the robustness of fitting molecular mechanics parameters labelling it as a ‘complex and subtle’ energy surface. 44 Thus, these authors have gone to great lengths to determine reliable parameters because of the ubiquity of the THF ring in furanose carbohydrates including ribose components in RNA and DNA. Experiment 45 has shown that the CS structure with a C–C–C–C dihedral angle of 0◦ is the global minimum; it is only at quite high levels of theory, MP2/aug-cc-pVTZ, that this is reflected by calculations which suggest that the C2 structure is the more stable. 46 However the structures are nearly isoenergetic with a difference of ≤ 0.2 kJ mol−1 . 47 And while there has been much work on aspects of the rotational spectroscopy not much has been done on the thermochemical implications, specifically on the nature of the ring puckering vibrations which in essence drive pseudo-rotation. Current codes do not readily lend themselves to a simple treatment of out-of-plane skeletal 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

modes. 48,49 In addition relaxed potential energy scans of the CCCC dihedral over a narrow range of values from −15 → +15◦ exhibit different behaviour at B3LYP, MN12SX, M06-2X and MP2 all at a common basis set of 6-311++G(d,p). Also Ghahremanpour et al. have argued that the pseudorotation inherent in five-membered rings is described relatively well by the harmonic approximation. 50 Consequently we have assumed that it is sufficient to consider ring puckering as pure vibrational modes although these low frequencies will have the greatest impact on the computed vibrational entropy.

Results The calculation of a reaction rate constant depends essentially on two quantities, (1) the difference in zero-point corrected electronic energies between the transition state, TS, and the reactant(s), A and B, more commonly known as the barrier height, E ‡ and (2) the ratio of partition functions, (QTS /QA QB ). The partition functions can be separated into individual contributions due to translational (T), rotational (R), vibrational (V) and electronic components. Of these the most critical are the vibrational terms which are very species dependent as opposed to translational or whole-molecule rotational partition functions. Unfortunately partition functions are not directly verifiable but their contribution to experimentally accessible quantities such as entropies and heat capacities enables validation.

Thermochemistry In the calculation of rate constants computed values for S and CP are tabulated and compared with the literature, Table 1. The calculated values include symmetries C2 for THF and DMTHF and C1 for 2MTHF and 3MTHF with the latter two being optically active. Only the methyl groups are treated as hindered rotors; the ring puckering modes are assumed to be pure vibrations except for THF itself for which a pseudo-rotational constant is known 51 and for which well characterised calculations are available. 52 The agreement is very good in

6 ACS Paragon Plus Environment

Page 6 of 34

Page 7 of 34

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

The Journal of Physical Chemistry

this particular case; a purely vibrational treatment gives S = 297.3 and CP = 80.79 J K−1 mol−1 which differ from the most accurate values. However since ring-opening, which would compromise the puckering modes, is not involved in these abstraction reactions we have assumed that neglecting the proper treatment will not materially affect the kinetic results which are based on the difference in the entropies of reactants and transition state. Table 1: Thermochemical data at 298.15 K S/J K−1 mol−1 calc. lit. THF 303.1 302.3 52 2MTHF 335.3 305.6 53 3MTHF 336.0 311.2 53 25DMTHF 357.5

Species

Cp /J K−1 mol−1 calc. lit. 76.7 77.2 52 105.9 103.2 53 106.4 95.1 53 132.4

For the transition states a hindered rotor analysis in Gaussian-16 shows that only rotation about the dihedral H···O – O – O needs to be considered and to a very good approximation a vibrational mode treatment is perfectly adequate. So for both reactants and transition states only methyl hindered rotors must be taken into account.

Energetics Although absolute enthalpies of formation are not required in determining the kinetics of reaction xTHF + O3 −−→ TS the following ∆f H(298.15 K) values were found for 2MTHF, 3MTHF and 25DMTHF of −225.1, −213.3 and −265.6 kJ mol−1 respectively, anchored by a value of −184.6 kJ mol−1 for THF. 54 A summary of the barrier heights at MN12SX/6-311++G(d,p), and the corresponding rate constants is given in Table 2 from which it can be seen that: 1. Abstraction from carbon atoms adjacent to the heterocyclic oxygen is the dominant site of attack 2. Followed by abstraction from the more remote carbon atoms and

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Table 2: Barriers / kJ mol−1 and rate constants / cm3 s−1 molecule−1 . E† k (298.15 K) tetrahydrofuran C2 ax 43.9 5.38 × 10−22 C2 eq 41.1 1.77 × 10−21 C3 ax 86.9 7.88 × 10−30 C3 eq 83.4 2.35 × 10−29 2-methyltetrahydrofuran C2 ax 27.4 1.49 × 10−19 CH3 103.0 1.95 × 10−33 C3 ax 82.8 9.06 × 10−30 C3 eq 88.4 5.53 × 10−31 C4 ax 77.4 6.06 × 10−29 C4 eq 88.0 1.61 × 10−30 C5 ax 34.1 5.78 × 10−21 C5 eq 43.4 1.85 × 10−22 3-methyltetrahydrofuran C2 ax 49.3 6.04 × 10−23 C2 eq 37.1 1.24 × 10−21 CH3 100.1 2.64 × 10−32 C3 ax 63.1 1.15 × 10−25 C4 ax 75.1 2.43 × 10−28 C4 eq 90.9 6.87 × 10−31 C5 ax 25.0 1.26 × 10−19 C5 eq 46.9 3.90 × 10−23 2,5-dimethyletrahydrofuran C2 ax 25.1 5.71 × 10−19 CH3 101.5 6.34 × 10−33 C3 ax 81.0 4.48 × 10−29 C3 eq 80.3 2.15 × 10−29 Site

3. Finally, least likely from the methyl groups (these always occupy an equatorial site). These observations are entirely consistent with known C – H bond dissociation energies which have C2 – H at 388–392, C3 – H at 412–416 and CH2 – H substantially higher at 431 kJ mol−1 and indeed they also parallel the barriers found for H-abstraction by H-atoms and CH3 radical. 54

8 ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

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

The Journal of Physical Chemistry

Kinetics There are very few quantitative measurements of the room temperature rates of reaction between saturated oxygenates and ozone. This is in contrast to the wealth of data for unsaturated oxygenates which typically react approximately 103 times faster. In this case by a different and well-known mechanism of addition across the double bond, ozonolysis, as elucidated by Criegee, which has synthetic significance. 55,56 As regards the total reactivity of each molecule with ozone comparisons can be made against very recent experimental measurements. The values recommended by Andersen and colleagues 8 are an average of determinations obtained from absolute measurements of the concentration of reactant as a function of time via FTIR in the presence of an excess of O3 and from relative experiments using ethyne as a standard. The rather large uncertainties in their final values are quite understandable given the difficulties in measuring the rates of these very slow reactions. Perhaps, somewhat fortuitously, the computed rate rate constant for THF lies within a factor of 3 of the experimental value. However in general the agreement is poor although the computed rate constants capture the relative reactivities of THF, 2MTHF and 25DMTHF in good agreement with experiment, albeit the absolute values are off by a large factor of up to ×10; note however that modest changes in barrier heights by ≈ 6 kJ mol−1 do bring these into very satisfactory agreement. These barriers do show large variabilities depending on the level of theory. Thus, for example, at M06-2X/6-311++G(d,p) the C2 axial and equatorial barriers for THF are 67.2 and 62.3 kJ mol−1 respectively, some 50% greater than the MN12SX barriers whereas the C3 sites only increase by about 20%. Even higher barriers emerge from ωB97X-D/aug-cc-pVTZ calculations. Computations of reactions of these oxolanes with ground state oxygen, 3O2 , show comparable patterns of reactivity at the same sites on the ring structures but at much slower rates as a result of barrier heights up to 5 fold higher than those calculated for reactions with ozone. 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Table 3: Experimental 8 and computed total rate constants at 298 K / cm3 s−1 molecule−1 Total rate constant, 10+21 × k reaction measured calc. THF 6.41 ± 2.90 2.3 2MTHF 18.7 ± 8.2 149 3MTHF — 126 25DMTHF 45.8 ± 21.8 571 Acyclic ethers In order to gain some insight into the cyclic ether, or THF results, abstraction from acyclic dimethyl and diethyl ethers was investigated. Judged purely on the barriers to reaction the cyclic ethers should be considerably more reactive than the acyclic ethers, Table 4. In the case of the highly reactive hydroxyl radical this distinction is blurred and the cyclic and acyclic ethers have comparable reactivities. 57,58 For diethyl ether abstraction of a secondary H-atom leads to the formation of 1-ethoxy ethanol, C2 H5 OCH(OH)CH3 , and O2 (1 ∆g ). The former was probably observed during a study on the atmospheric oxidation of diethyl ether initiated by chlorine atom but the author could not be certain. 59 However the unidentified IR bands match quite closely our calculated spectrum. Abstraction of a primary H-atom in dimethyl ether leads this time to the formation of methoxymethylhydrotrioxide, CH3 OCH2 OOOH. In neither case are pre-reaction complexes observed. Alkyl hydroperoxides Since total reaction rates for CH3 OOH + O3 and C2 H5 OOH + O3 of 7.3×10−21 and 8.4×10−21 cm3 molecule−1 s−1 respectively, were available in the literature 20,21 we determined barrier heights and IRCs for these reactions. The computed rate constants are not in good agreement with these literature values for reasons previously discussed. It is however somewhat surprising that ethyl hydroperoxide is not more reactive than experiment indicates since the barrier to abstraction faced by a secondary H-atom in ethyl hydroperoxide at 65 kJ mol−1 is

10 ACS Paragon Plus Environment

Page 10 of 34

Page 11 of 34

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

The Journal of Physical Chemistry

Table 4: Barriers / kJ mol−1 . Species THF C2 CH3 CH2 OC2 H5 CH3 OCH3 CH3 CH2 OOH CH3 CH2 OOH CH3 OOH CH3 CH2 OH CH3 CH2 OH CH3 OH

E† 41.1 55.5 58.3 65.3 98.8 76.2 56.0 80.5 60.9

quite lower than that by a primary H-atom in methyl hydroperoxide at 76 kJ mol−1 , Table 4. For methyl hydroperoxide abstraction of a primary H-atom leads to the formation of a αperoxymethylhydrotrioxide, HOOOCH2 OOH. The simpler α-hydroxymethylhydroperoxide, HOCH2 OOH, has been the focus of recent studies on its catalysed and un-catalysed decomposition in the atmosphere. 60 For ethylhydroperoxide C2 H5 OOH abstraction of a secondary H-atom forms 1-hydroperoxy ethanol and singlet oxygen, CH3 CH(OH)OOH + O2 . Alcohols In an attempt to select a system which might be computationally affordable with higher level methods we have found that H-abstraction from methanol, CH3 OH, exhibits a comparable barrier of 61 kJ mol−1 to that for dimethyl ether of 58 kJ mol−1 . W2X and W3X-L calculations at an MN12SX/6-311++G(d,p) geometry give barrier heights of 76.3 and 71.1 kJ mol−1 — a clear indication of the necessity of including a post-CCSD(T) treatment. However the nature of the transition state (cyclic) and the IRC indicates the direct formation of methanal and trioxidane H2 CO + HOOOH, effectively abstracting two H-atoms from CH3 OH in a single step, Fig. 2. By contrast for ethanol, abstraction of a secondary Hatom is more straightforward with reaction leading to the hydrotrioxide CH3 CH(OOOH)OH.

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 34

Figure 2: CH3 OH + O3 transition state. Formaldehyde CCSD(T)/6-311+G(d,p)//M05-2X/6-311+G(d,p) calculations by Voukides et al. 18 and BMCCCSD//BHandHLYP/6-311+G(d,p) by Wang and colleagues 17 showed that both addition and abstraction occurs in CH2 – O + O3 , Fig. 3. They both agree that abstraction faces lower barriers of 68 and 72 kJ mol−1 , respectively, but disagree on the magnitude of this difference. In this study a value of 54.3 kJ ml−1 was obtained for abstraction at MN12SX/6311++G(d,p) and of 76.3 kJ mol−1 for W2X — yet another illustration of the sensitivity of these reactions.

Figure 3: CH2 O + O3 transition states.

Products Although the products of reaction are of very little kinetic significance as regards these calculations, the nature of the products of xTHFs + O3 is very interesting. Intrinsic reaction 12 ACS Paragon Plus Environment

Page 13 of 34

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

The Journal of Physical Chemistry

coordinate computations show that as the initially formed hydrogen trioxide radical, HOOO• , departs it can re-attach to the newly formed radical carbon C• in one of two ways. If the terminal O-atom attaches then a hydrotrioxide is formed, ROOOH, whereas if the other ‘terminal’ O-atom is involved then a furanol will result concomitantly with the elimination of singlet dioxygen, R – OH + O2 (1 ∆g ). For 25DMTHF abstraction at a C3 equatorial site neatly illustrates the fine balance between these two possibilities, Fig 4, C···O(H) = 2.608 ≈ C···O(OOH) = 2.612 Å; as the IRC is traversed the hydrotrioxide is formed. For the case of abstraction at the C3 axial site the IRC goes on to eliminate O2 and form the furanol.

Figure 4: A snapshot along an IRC path This contrasting behaviour can be simply shown in plots of the C···H, C···O(H), C···O(OOH) and O···O(OH) bondlengths, Fig 5, as a function of energy along the IRCs. In the only other known case Cerkovnik and colleagues 10 found that a trioxide is formed, benzoyl hydrotrioxide, during the low temperature ozonation of benzaldehyde in preference to the formation of a cyclic tetraoxolane. However, the trioxide rapidly decomposes in an 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Figure 5: Bondlengths along IRC paths; ⊙ C3 axial,

14 ACS Paragon Plus Environment

Page 14 of 34

C3 equatorial

Page 15 of 34

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

The Journal of Physical Chemistry

intramolecular H-atom transfer to benzoic acid and singlet dioxygen and they concluded that ROOOH should not be observable by 1H,

13

C and

17

O NMR above −80 C in a number

of different organic solvents, contrary to previous reports.

This latter mechanism is not possible for our systems but the equivalent process, formation of an alcohol and elimination of dioxygen is. The prediction of hydrotrioxides as products from our theoretical calculations is somewhat surprising. The chemistry of organic hydrotrioxides is but sparsely known 61,62 and there are very few examples of their formation particularly via an H-atom abstraction reaction. Hydrotrioxides The only member of this family sufficiently well known to be included in the Active Thermochemical Tables 40 or to feature in a recent review 63 is trioxidane or dihydrogen trioxide, HOOOH, the trans conformer of which has a formation enthalpy, ∆f H(0 K), of −81.40±0.78 kJ mol−1 ; this is in good agreement with a W3X-L value of −82.1 kJ mol−1 . A recent communication considered the reaction between methyl peroxy and hydroxyl radicals from which an enthalpy of formation for methyl hydrotrioxide of −65.3±3.5 kJ mol−1 can be derived. 64 Earlier Jungkamp and Seinfeld 65 had derived values for ∆f H(298.15K) of −94 kJ mol−1 from G2M(RCC) calculations. None of these are in satisfactory accord with W2X and W3X-L calculations of −69.5 and −71.1 kJ mol−1 at 0 K and −84.7 and −86.9 kJ mol−1 at 298.15 K respectively, for the same conformer. These key values can now be used in isodesmic reactions with lower, more affordable, levels of theory to compute the formation enthalpies of larger hydrotrioxides. The bond dissociation energies of 3-substituted tetrahydrofuran hydrotrioxides, hydroperoxides and alcohols, calculated with the G4 model chemistry, 66 are shown in Table 5. Although most of these values are novel where comparisons are possible, for example, in the 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 34

case of tetrahydro-3-furanol, R–OH typically range from 390–400 kJ mol−1 and for RO–H 435-445 kJ mol−1 are not untypical. 67 The O–OH bond in tetrahydro-3-hydroperoxyfuran, RO–OH, at 177 kJ mol−1 is comparable to values of 180–190 kJ mol−1 to be found in a comprehensive tabulation of chemical bond energies. The only direct comparison that can be made is with methyl hydrotrioxide for which Lay and Bozzelli 68 calculated CH3 OO – OH of 143 kJ mol−1 and 124 kJ mol−1 for CH3 O – OOH using a variety of methods up to G2 model chemistry. Table 5: Bond dissociation energy, BDE(298.15 K) / kJ mol−1 ROOO–H ROO–OH RO–OOH R–OOOH

319 121 104 251

ROO–H RO–OH R–OOH

353 177 286

RO–H R–OH

440 390

Conclusions The reactivities of the various sites in the oxolanes THF, 2-MTHF, 3-MTHF and 25DMTHF is mapped out and it is shown that these are concentrated at those carbon atoms situated beside the heterocyclic oxygen with abstraction from the methyl groups being least likely. All the functionals and basis set combinations tested are in broad relative agreement with this conclusion although they do differ widely as regards their absolute values. Hence, absolute computed rate constants are not in good agreement with experiment almost certainly due to the limited theoretical treatment; however modest but consistent changes to the barrier heights does improve the situation. The need for the highest levels of theory to adequately determine the potential energy surface was a difficulty which could only be overcome by reducing the number of ‘heavy’ atoms to two in the reactant, methanol. The fact that the reaction mechanism then changes from the usual abstraction of a single H-atom renders a general finding somewhat less certain. However a 10 kJ mol−1 change in barrier height leads to a 50-fold change in rate constant 16 ACS Paragon Plus Environment

Page 17 of 34

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

The Journal of Physical Chemistry

at room temperature — an illustration of the extreme sensitivity encountered with these reactions. Moreover, the difference in barrier heights from W2X, from all-electron scalarrelativistic CCSD(T)/CBS energies, to W3X-L, which includes in addition a multi-reference treatment up to CCSDT(Q)/VDZ, leads to an 8-fold change in rate constant. In comparison to the cyclic ethers abstraction of a secondary hydrogen from an acyclic ether such as C2 H5 OC2 H5 exhibits a much higher barrier as does abstraction of secondary hydrogens from peroxides CH3 CH2 OOH and alcohols, CH3 CH2 OH, Table 4. Abstraction of primary hydrogens from acyclics, peroxides and alcohols are yet higher again reflecting the stronger C–H bonds. The nature of the products initially formed also varies with remote H-sites favouring the formation of a furanic hydrotrioxide and the vicinal sites favouring the formation of furanols. The properties of these novel hydrotrioxides are shown to correspond with some of the limited data available for the simplest alkyl hydrotrioxide CH3 OOOH. In this latter case we are able to provide an accurate value for its formation enthalpy using a methodology which is itself validated against the best collection of such data. 40

Acknowledgement The Irish Centre for High-End Computing, ICHEC, is thanked for the provision of computational resources to projects gmche001c and ngche041c. We also acknowledge advice from J.L Bao and D. G. Truhlar.

Supporting Information Available Cartesian coordinates, vibrational frequencies, symmetries and rotational constants for all reactive species are archived.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

References (1) Capello, C.; Fischer, U.; Hungerbuhler, K., What Is a Green Solvent? A Comprehensive Framework for the Environmental Assessment of Solvents. Green Chemistry 2007, 9, 927–934. (2) Slater, C. S.; Savelski, M. J.; Hitchcock, D.; Cavanagh, E. J., Environmental Analysis of the Life Cycle Emissions of 2-Methyl Tetrahydrofuran Solvent Manufactured from Renewable Resources. J. Environ. Sci. Health Part a-Toxic/Haz. Substances & Environ. Eng. 2016, 51, 487–494. (3) Pace, V.; Hoyos, P.; Castoldi, L.; de Maria, P. D.; Alcantara, A. R., 2Methyltetrahydrofuran (2-MeTHF): A Biomass-Derived Solvent with Broad Application in Organic Chemistry. ChemSusChem 2012, 5, 1369–1379. (4) Ulonska, K.; Voll, A.; Marquardt, W., Screening Pathways for the Production of Next Generation Biofuels. Energy & Fuels 2016, 30, 445–456. (5) Voll, A.; Marquardt, W., Benchmarking of Next-Generation Biofuels from a Process Perspective. Biofuels, Bioprod. Biorefin. 2012, 6, 292–301. (6) Janssen, A. J.; Kremer, F. W.; Baron, J. H.; Muether, M.; Pischinger, S.; Klankermayer, J., Tailor-Made Fuels from Biomass for Homogeneous Low-Temperature Diesel Combustion. Energy & Fuels 2011, 25, 4734–4744. (7) Moriarty, J.; Sidebottom, H.; Wenger, J.; Mellouki, A.; Le Bras, G. Kinetic Studies on the Reactions of Hydroxyl Radicals with Cyclic Ethers and Aliphatic Diethers J. Phys. Chem. A 2013, 107, 1499–1505. (8) Andersen, C.; Nielsen, O. J.; Oesterstroem, F. F.; Ausmeel, S.; Nilsson, E. J. K.; Sulbaek Andersen, M. P., Atmospheric Chemistry of Tetrahydrofuran, 2-

18 ACS Paragon Plus Environment

Page 18 of 34

Page 19 of 34

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

The Journal of Physical Chemistry

Methyltetrahydrofuran, and 2,5-Dimethyltetrahydrofuran: Kinetics of Reactions with Chlorine Atoms, OD Radicals, and Ozone. J. Phys. Chem. A 2016, 120, 7320–7326. (9) Manion, J. A.; Huie, R. E. ; Levin, R. D.; Burgess Jr., D. R. ; Orkin, V. L.; Tsang, W.; McGivern, W. S.; Hudgens, J. W.; Knyazev, V. D.; Atkinson, D. B.; Chai, E. ; Tereza, A. M.; Lin, C.-Y.; Allison, T. C.; Mallard, W. G.; Westley, F.; Herron, J. T.; Hampson, R. F.; Frizzell, D. H. NIST Chemical Kinetics Database, NIST Standard Reference Database 17, Version 7.0 (Web Version), Release 1.6.8, Data version 2015.12, National Institute of Standards and Technology, Gaithersburg, Maryland, 20899-8320. Web address: http://kinetics.nist.gov/ (10) Cerkovnik, J.; Plesničar, B.; Koller, J.; Tuttle, T., Hydrotrioxides Rather Than Cyclic Tetraoxides (Tetraoxolanes) as the Primary Reaction Intermediates in the LowTemperature Ozonation of Aldehydes. The Case of Benzaldehyde. J. Org. Chem. 2009, 74, 96–101. (11) Grosjean, D., Atmospheric Chemistry of Biogenic Hydrocarbons — Relevance To The Amazon. Química Nova 1995, 18, 184–201. (12) Grosjean, E.; Grosjean, D., Rate constants for the gas-phase reaction of ozone with unsaturated oxygenates. Int. J. Chem. Kinet. 1998, 30, 21–29. (13) Al Mulla, I.; Viera, L.; Morris, R.; Sidebottom, H.; Treacy, J.; Mellouki, A., Kinetics and mechanisms for the reactions of ozone with unsaturated oxygenated compounds. ChemPhysChem 2010, 11, 4069–4078. (14) Yang, J.; Li, Q. S.; Zhang, S. W., Direct Dynamics Study on the Reaction of Acetaldehyde with Ozone. J. Comput. Chem. 2008, 29, 247–255. (15) Yang, J.; Zhang, S. W.; Li, Q. S., Reaction Path Dynamics and Theoretical Rate Constants Calculation for CHn F(4 – n) + O3 −−→ HOOO + CH(n – 1) F(4 – n) (n=2,3) Reactions. Chem. J. Chinese Universities-Chinese 2007, 28, 1975–1977. 19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(16) Yang, J.; Li, Q. S.; Zhang, S. W., Reaction-Path Dynamics and Theoretical Rate Constants for the Reaction CH4 + O3 −−→ HOOO + CH3 . Int. J. Quantum Chem. 2007, 107, 1999–2005. (17) Wang, F.; Sun, H.; Sun, J. Y.; Jia, X. J.; Zhang, Y. J.; Tang, Y. Z.; Pan, X. M.; Su, Z. M.; Hao, L. Z.; Wang, R. S., Mechanistic and Kinetic Study of CH2 O + O3 Reaction. J. Phys. Chem. A 2010, 114, 3516–3522. (18) Voukides, A. C.; Konrad, K. M.; Johnson, R. P., Competing Mechanistic Channels in the Oxidation of Aldehydes by Ozone. J. Organic Chem. 2009, 74, 2108-2113. (19) Wang, J.; Zhou, L.; Wang, W.; Ge, M., Gas-phase reaction of two unsaturated ketones with atomic Cl and O3 : kinetics and products. Phys. Chem. Chem. Phys. 2015, 17, 12000–12012. (20) Chen, Z.; Wang, C., Rate constants of the gas-phase reactions of CH3 OOH with O3 and NOx at 293 K. Chem. Phys. Letts. 2006, 424, 233–238. (21) Wang, C.; Chen, Z., An experimental study for rate constants of the gas phase reactions of CH3 CH2 OOH with OH radicals, O3 , NO2 and NO. Atmos. Environ. 2008, 42, 6614– 6619. (22) Vereecken, L.; Francisco, J. S., Theoretical Studies of Atmospheric Reaction Mechanisms in the Troposphere. Chem. Soc. Revs. 2012, 41, 6259–6293. (23) Zhao, Y.; Lynch, B. J.; Truhlar, D. G., Multi-coefficient extrapolated density functional theory for thermochemistry and thermochemical kinetics. Phys. Chem. Chem. Phys. 2005, 7, 43–52. (24) Wang, Z. H.; Yang, L.; Li, B.; Li, Z. S.; Sun, Z. W.; Alden, M.; Cen, K. F.; Konnov, A. A., Investigation of Combustion Enhancement by Ozone Additive in CH4 /Air

20 ACS Paragon Plus Environment

Page 20 of 34

Page 21 of 34

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

The Journal of Physical Chemistry

Flames Using Direct Laminar Burning Velocity Measurements and Kinetic Simulations. Combust. Flame 2012, 159, 120–129. (25) Zhang, Y.; Zhu, M. M.; Zhang, Z. Z.; Shang, R. X.; Zhang, D. K., Ozone Effect on the Flammability Limit and near-Limit Combustion of Syngas/Air Flames with N2 , CO2 , and H2 O Dilutions. Fuel 2016, 186, 414–421. (26) Zhang, K. P.; Hu, G.; Liao, S. Y.; Zuo, Z. H.; Li, H. M.; Cheng, Q.; Xiang, C. X., Numerical Study on the Effects of Oxygen Enrichment on Methane/Air Flames. Fuel 2016, 176, 93–101. (27) Pinchak, M.; Ombrello, T.; Carter, C.; Gutmark, E.; Katta, V., The Effects of Hydrodynamic Stretch on the Flame Propagation Enhancement of Ethylene by Addition of Ozone. Phil. Trans. Royal Society A-Math. Phys. Eng. Sci. 2015, 373 (2048). (28) Gao, X.; Zhang, Y.; Adusumilli, S.; Seitzman, J.; Sun, W. T.; Ombrello, T.; Carter, C., The Effect of Ozone Addition on Laminar Flame Speed. Combust. Flame 2015, 162, 3914–3924. (29) Peverati, R.; Truhlar, D. G. Screened-Exchange Density Functionals with Broad Accuracy for Chemistry and Solid-State Physics Phys. Chem. Chem. Phys. 2012, 14, 16187–16191. (30) Gaussian 16, Revision A.03, 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.; 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

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.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian, Inc., Wallingford CT, 2016. (31) Fukui, K., The Path of Chemical Reactions — the IRC Approach. Acc. Chem. Res. 1981, 14, 363–368. (32) Hratchian, H. P.; Schlegel, H. B., Using Hessian Updating to Increase the Efficiency of a Hessian Based Predictor-Corrector Reaction Path Following Method. J. Chem. Theor. Comput. 2005, 1, 61–69. (33) Simon, S.; Duran, M.; Dannenberg, J. J., How Does Basis Set Superposition Error Change the Potential Surfaces for Hydrogen Bonded Dimers? J. Chem. Phys. 1996, 105, 11024–11031. (34) Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M. WIREs Comput Mol Sci. 2012, 2, 242–-253. MOLPRO, version 2012.1, a package of ab initio programs, Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G., et al. http://www.molpro.net. 16/04/2012 (35) Adler, T. B.; Knizia, G.; Werner, H. J., A Simple and Efficient CCSD(T)-F12 Approximation. J. Chem. Phys. 2007, 127, 221106. (36) Zhao, Y.; Truhlar, D. G., The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Accts 2008, 120, 215–241. 22 ACS Paragon Plus Environment

Page 22 of 34

Page 23 of 34

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

The Journal of Physical Chemistry

(37) Lee, T. J.; Taylor, P. R., A Diagnostic for Determining the Quality of Single-reference Electron Correlation Methods. Int. J. Quantum Chem. 1989, 36, 199–207. (38) Chan, B.; Radom, L., W2x and W3x-L: Cost-Effective Approximations to W2 and W4 with kJ mol−1 Accuracy. J. Chem. Theor. Comput. 2015, 11, 2109–2119. (39) Rolik, Z.; Szegedy, L.; Ladjánszki, I.; Ladóczki, B.; Kállay, M. J. Chem. Phys. 2013, 139, 094105. MRCC, a quantum chemical program suite written by Kállay,M.; Rolik, Z.; Csontos, J.; Ladjánszki, I.; Szegedy, L.; Ladóczki, B.; Samu, G. http://www.mrcc.hu 16-September-2015 (40) Ruscic, B. and Bross, D. H., Active Thermochemical Tables (ATcT) values based on ver. 1.122 of the Thermochemical Network (2016); available at http://atct.anl.gov/ ThermochemicalData/ accessed 25-May-2017. (41) Simmie, J. M.; Sheahan, J.N. Validation of a Database of Formation Enthalpies and of Mid-level Model Chemistries, J. Phys. Chem. A 2016, 120, 7370–7384. (42) Barker, J. R.; Nguyen, T. L.; Stanton, J. F.; Aieta, C.; Ceotto, M.; Gabas, F.; Kumar, T. J. D.; Li, C. G. L.; Lohr, L. L.; Maranzana, A., Ortiz, N. F.; Preses, J.M.; Simmie, J. M.; Sonk, J. A.; Stimac, P. J. Multiwell-2017 Software Suite, Ann Arbor, Michigan, USA, 2017. http://clasp-research.engin.umich.edu/multiwell/ (43) Somers, K. P. 2014, Ph. D. thesis National University of Ireland, Galway. (44) Vanommeslaeghe, K.; Yang, M.; MacKerell, Jr A. D., Robustness in the Fitting of Molecular Mechanics Parameters, J. Comp. Chem. 2015, 36, 1083–1101. (45) Yang, T. C.; Su, G. L.; Ning, C. G.; Deng, J. K.; Wang, F.; Zhang, S. F.; Ren, X. G.; Huang, Y. R., New Diagnostic of the Most Populated Conformer of Tetrahydrofuran in the Gas Phase. J. Phys. Chem. A 2007, 111, 4927–4933.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(46) Rayón, V. M.; Sordo, J. A., Pseudorotation Motion in Tetrahydrofuran: An Ab Initio Study. J. Chem. Phys. 2005, 122, 204303. (47) Ford, T. A., The Evolution of the Structural, Vibrational and Electronic Properties of the Cyclic Ethers — on Ring Size. An ab Initio Study. J. Mol. Struct. 2014, 1073, 125–133. (48) Harris, D. O.; Engerholm, G. G.; Tolman, C. A.; Luntz, A. C.; Keller, R. A.; Kim, H.; Gwinn, W. D., Ring Puckering in 5-Membered Rings .1. General Theory. J. Chem. Phys. 1969, 50, 2438–2445. (49) Engerholm, G. G.; Luntz, A. C.; Gwinn, W. D.; Harris, D. O., Ring Puckering in 5-Membered Rings .2. Microwave Spectrum Dipole Moment and Barrier to Pseudorotation in Tetrahydrofuran. J. Chem. Phys. 1969, 50, 2446–2457. (50) Ghahremanpour, M. M.; van Maaren, P. J.; Ditz, J. C.; Lindh, R.; van der Spoel, D., Large-scale Calculations of Gas Phase Thermochemistry: Enthalpy of Formation, Standard Entropy, and Heat Capacity. J. Chem. Phys. 2016, 145, 114305. (51) Melnik, D. G.; Gopalakrishnan, S.; Miller, T. A.; De Lucia, F. C., The Absorption Spectroscopy of the Lowest Pseudorotational States of Tetrahydrofuran. J. Chem. Phys. 2003, 118, 3589–3599. (52) Dorofeeva, O. V., Ideal-Gas Thermodynamic Properties of Oxygen HeterocyclicCompounds Part 1. 3-Membered, 4-Membered and 5-Membered Rings. Thermochim. Acta 1992, 194, 9–46. (53) Auzmendi-Murua, I.; Charaya, S.; Bozzelli, J. W., Thermochemical Properties of Methyl-Substituted Cyclic Alkyl Ethers and Radicals for Oxiranes, Oxetanes, and Oxolanes: C–H Bond Dissociation Enthalpy Trends with Ring Size and Ether Site. J. Phys. Chem. A 2013, 117, 378–392.

24 ACS Paragon Plus Environment

Page 24 of 34

Page 25 of 34

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

The Journal of Physical Chemistry

(54) Simmie, J. M. Kinetics and Thermochemistry of 2,5-Dimethyltetrahydrofuran and Related Oxolanes: Next Next-Generation Biofuels J. Phys. Chem. A 2012, 116, 4528– 4538. (55) Jones Jr., M. Organic Chemistry, 2nd ed., pg. 387, W. W. Norton & Co., New York, 2000. (56) Criegee, R., Mechanism of Ozonolysis. Angew. Chem-Int. Ed. in English 1975, 14, 745–752. (57) Mellouki, A.; Teton, S.; Lebras, G., Kinetics of OH Radical Reactions with a Series of Ethers. Int. J. Chem. Kinet. 1995, 27, 791–805. (58) Moriarty, J.; Sidebottom, H.; Wenger, J.; Mellouki, A.; Le Bras, G., Kinetic Studies on the Reactions of Hydroxyl Radicals with Cyclic Ethers and Aliphatic Diethers. J. Phys. Chem. A 2003, 107, 1499–1505. (59) Orlando, J. J. The Atmospheric Oxidation of Diethyl Ether: Chemistry of the C2 H5 – O – CH(O)CH3 Radical between 218 and 335 K. Phys. Chem. Chem. Phys. 2007, 9, 4189–4199. (60) Kumar, M.; Busch, D. H.; Subramaniam, B.; Thompson, W. H., Role of Tunable Acid Catalysis in Decomposition of α-Hydroxyalkyl Hydroperoxides and Mechanistic Implications for Tropospheric Chemistry. J. Phys. Chem. A 2014, 118, 9701–9711. (61) Shereshovets, V.; L Khursan, S.; Komissarov, V.; Tolstikov, G., Organic Hydrotrioxides. Russ. Chem. Revs. 2001, 70, 105–129. (62) Khalitova, L. R.; Grabovskiy, S. A.; Antipin, A. V.; Spirikhin, L. V.; Kabal’nova, N. N., Products of Ozone Oxidation of Some Saturated Cyclic Hydrocarbons. Russian J. Org. Chem. 2015, 51, 1710–1716.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

(63) Burgess, D. R., An Evaluation of Gas Phase Enthalpies of Formation for HydrogenOxygen (HxOy) Species. J. Research National Inst. Standards and Technol. 2016, 121, 108–138. (64) Müller, J.-F.; Liu, Z.; Nguyen, V. S.; Stavrakou, T.; Harvey, J. N.; Peeters, J., The Reaction of Methyl Peroxy and Hydroxyl Radicals as a Major Source of Atmospheric Methanol. Nature Comms. 2016, 7, 13213. (65) Jungkamp, T. P. W.; Seinfeld, J. H., The Enthalpy of Formation of Trioxy Radicals ROOO (R=H, CH3 , C2 H5 ). An ab initio Study. Chem. Phys. Letts. 1996, 257, 15–22. (66) Curtiss, L. A.; Redfern, P. C.; Raghavachari, K., Gaussian-4 Theory. J. Chem. Phys. 2007, 126, 084108. (67) Luo, Y.-R., Comprehensive Handbook of Chemical Bond Energies. CRC Press: Boca Raton, FL, 2007. (68) Lay, T. H.; Bozzelli, Enthalpies of Formation and Group Additivity of Alkyl Peroxides and Trioxides. J. W. J. Phys. Chem. A 1997, 101, 9505–9510.

26 ACS Paragon Plus Environment

Page 26 of 34

Page 27 of 34

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

The Journal of Physical Chemistry

TOC Graphic

Figure 6: TOC Graphic

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

Page 28 of 34

Page 29 of 34

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

Page 30 of 34

Page 31 of 34

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

Page 32 of 34

Page 33 of 34

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

ACS Paragon Plus Environment

Page 34 of 34