Kinetics of Thermal Unimolecular Decomposition of Acetic Anhydride

Apr 6, 2017 - (8). The master equation (ME) has been shown to be a powerful ... valence cc-pVXZ basis sets,(20, 21) and A, B, α and β are constants ...
0 downloads 0 Views 2MB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article

Kinetics of Thermal Unimolecular Decomposition of Acetic Anhydride - An Integrated Deterministic and Stochastic Model Tam V.-T. Mai, Minh van Duong, Hieu Thanh Nguyen, Kuang C. Lin, and Lam K Huynh J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b00015 • Publication Date (Web): 06 Apr 2017 Downloaded from http://pubs.acs.org on April 7, 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 31

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 of Thermal Unimolecular Decomposition of Acetic Anhydride An Integrated Deterministic and Stochastic Model

Tam V.-T. Mai,1 Minh v. Duong,1 Hieu T. Nguyen,1 Kuang C. Lin2 and Lam K. Huynh3,* 1

Institute for Computational Science and Technology, SBI Building, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam.

2

Department of Mechanical and Electromechanical Engineering, National Sun Yat-Sen University, Kaohsiung, 80424, Taiwan.

3

International University, Vietnam National University - HCMC, Quarter 6, Linh Trung Ward, Thu Duc District, Ho Chi Minh City, Vietnam.

* Corresponding authors. Email address: [email protected] | [email protected] (LKH) Tel: (84-8) 2211.4046 (Ext. 3233) Fax: (84-8) 3724.4271

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 2 of 31

Abstract An integrated deterministic and stochastic model, within the master equation/Rice–Ramsperger– Kassel–Marcus (ME/RRKM) framework, was first used to characterize temperature- and pressuredependent behaviors of thermal decomposition of acetic anhydride in a wide range of conditions (i.e., 300–1500 K & 0.001–100 atm). Particularly, using potential energy surface and molecular properties obtained from high-level electronic structure calculations at CCSD(T)/CBS, macroscopic thermodynamic properties and rate coefficients of the title reaction were derived with corrections for hindered internal rotation and tunneling treatments. Being in excellent agreement with the scattered experimental data, the results from deterministic and stochastic frameworks confirmed and complemented each other to reveal the main decomposition pathway proceeds via a 6-membered-ring transition state with the 0-K barrier of 35.2 kcal⋅mol-1. This observation was further understood and confirmed by the sensitivity analysis on the time-resolved species profiles and the derived rate coefficients with respect to the ab initio barriers. Such an agreement suggests the integrated model can be confidently used for a wide range of conditions as a powerful post-facto and predictive tool in detailed chemical kinetic modeling and simulation for the title reaction and thus can be extended to complex chemical reactions.

2

ACS Paragon Plus Environment

Page 3 of 31

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

1. Introduction Anhydride acetic, (CH3CO)2O (referred to as AH), can be considered as an acetylating reagent1,2 that thermally decomposes to give ketene and acetic acid, (CH3CO)2O → CH2CO + CH3COOH (Rxn. 1)3-7. This reaction has truly received considerable interest for the purpose to produce ketene (CH2CO)5, an important intermediate in the combustion of hydrocarbon fuels. Despite its importance, only a few kinetic investigations concerning the title reaction have been studied both experimentally3-5 and theoretically5. Using a flow reactor, Szwark and Murawski3 measured a rate constant of k(T) = 1.00×1012×exp(-34500 cal⋅mol-1/RT) s-1 within the temperature range of 553 – 646 K and the total pressure range of 6 – 13 torr. By analyzing the reaction products, they found the decomposition took place in the unimolecular process to yield equal molarity of ketene and acetic acid with no other products observed. Blake and Speis4 conducted this reaction in a static and flow apparatus at the temperature of 470 – 643 K and over a wide pressure range of 7 – 232 torr, leading to the rate constant k(T) = 1.48×1011×exp(-32194 cal⋅mol-1/RT) s-1. They also found that the only products were ketene and acetic acid. More recently, Akao et al.5 reported a rate constant expression as k(T) = 6.31×1011×exp(-32983 cal⋅mol-1/RT) s-1 by using the shock tube method to evaluate the Arrhenius parameters of this composition at higher temperature of 750 – 980 K and pressure range of 19 – 80 torr. They also performed ab initio molecular orbital calculations by identifying the conceivable transition states (TSs) and discussed the reaction mechanism. To our best knowledge, there is no detailed computational study on the mechanism and kinetics of the title reaction even though a 6membered-ring transition state has been suggested for this reaction by analogy with ester elimination8.

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

The master equation (ME) has been shown to be a powerful framework for modelling kinetic behaviors of a complex gas-phase chemical system for a wide range of conditions (e.g., temperature and pressure). The solutions of the ME equation, i.e., the species profiles as well as the rate coefficients, can be obtained by either deterministic9,10 or stochastic11,12 approaches. In principles, it is challenging for a deterministic model to provide reliable results if there is no clear separation between the two time scales, namely internal energy relaxation eigenvalues (IEREs) and chemically significant eigenvalues (CSEs) (e.g., the time-scale difference is less than one order of magnitude which normally occurs at high temperature and low pressure). Alternatively, a stochastic model can reveal the real physics/chemistry only with a statistically–sufficient number of trials, which normally requires much computational time and effort to correctly describe minor events or channels/species. Therefore, the use of both approaches will help to complement each other, and thus gain more confidence on the calculated numbers for a wide range of different conditions, even at the extreme cases as mentioned above. In such a case, they certainly will help to gain more confidence in the calculated data as well as to shed more light to the nature of complicated processes. In this study, we explored the mechanism of the title reaction by carrying out accurate ab initio calculations whose results are then used to characterize its pressure and temperaturedependent behaviors in a wide range of conditions using the combined deterministic and stochastic statistical rate model. It is expected to provide an integrated picture to bridge the gap between fundamental chemistry and reaction engineering for this significant issue in atmospheric chemistry and combustion.

4

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

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

2. Computational Details Structures of reactants, products and TSs of the title reaction were optimized using B3LYP level of theory13,14 and 6-311G(2d,d,p) basis set denoted as CBSB7. The calculated B3LYP/CBSB7 vibrational frequencies with the scaling factor of 0.9915 were used to calculate the zero-point energy (ZPE) corrections. Electronic energies were calculated with the B3LYP geometries by using the coupled cluster theory with single and double excitations16-18 and perturbative estimate of triple excitations [CCSD(T)]19.

Extrapolation for the calculated

electronic energies was further performed to complete basis set (CBS) using cc-pVDZ and ccpVTZ20,21 basis sets suggested by Truhlar22. The total single point energy is a sum of the HartreeFock (HF) energy, EHF and the correlation energy, Ecorr: E∞tot = E∞HF + E∞corr. The HF energy is assumed to approach its CBS limit by power laws: ExHF = E∞HF + AX-α and Excorr = E∞corr + BX-β, where X = 2, 3 for D, T extrapolation with the Dunning’s correlation consistent polarized valance cc-pVXZ basis sets20,21, and A, B, α and β are constants (α = 3.4 and β = 2.4 for CCSD and CCSD(T) calculations)22. The reported root-mean-square (RMS) deviation of the CCSD(T)/CBS procedure was from 1.3 to 2.4 kcal/mol for the total energies of the considered test cases22. The energetic profiles extrapolated from CCSD(T)/CBS calculations, were then compared with other accurate composite methods such CBS-QB315, CBS-APNO23, G324 and G425 levels of theory. The CCSD(T)/CBS energies, accompanied with the B3LYP/CBSB7 geometries and frequencies information, were used to carry out kinetic analysis. All these electronic structure calculations were calculated using the GAUSSIAN 09 package26. Thermodynamic and kinetic analyses, using both deterministic and stochastic models, were carried out using the Multi-Species Multi-Channels (MSMC) code27. The hindered internal rotation (HIR) correction28 and tunneling effect29 were taken into account. For the HIR treatment, 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

the hindrance potentials, V(θ) , as a function of torsional angle, θ, along the single bonds (i.e., C−C and C−O, cf. supplementary Figure S1) were explicitly obtained at the B3LYP/6-31G(d) level via relaxed surface scans with the step size of 10o for the dihedral angles corresponding to the rotations. The HIR parameters were determined with the use of MSMC-GUI30. The calculated m‡ σ external rate constants included the reaction path degeneracy, L = × ‡ (where m‡ and m are the m σ external #

‡ number of optical isomers of transition states and reactant, and σ exter nal and σ external are the external

rotation symmetry numbers of the transition state and reactant, respectively). Note that the internal symmetry number σ internal (e.g., σ internal = 3 for the rotation of the methyl CH3 group along the C-C bond) was explicitly included in the HIR treatment. The energy-transfer process was computed on the basis of the exponential-down model with Edown = 500 cm-1 for Ar as the bath gas31. The Lennard-Jones (L-J) parameters Ar32 while

ε / kB = 93.3 K and σ = 3.542 Å were used for

ε / kB = 454.0 K and σ = 5.413 Å were taken from the data of C4H7O333 to represent

C4H6O3. The density/sum of states were directly counted using the Beyer-Swinehart algorithm34 (for HO modes) in combination/convolution with the modified Stein-Rabinovitch algorithm35 (for external/hindered rotational modes) with the energy step of 1 cm-1. The maximum energy level was determined by the convergence tolerance ε = 1.0×10-6 in the sum of the partition function, which is temperature dependent. For the title system, maximum energy level at T = 1500 K is roughly 40,000 cm-1. For kinetic calculations, the energy grain size of 200 cm-1 was used.

6

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

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

3. Results and Discussion Potential Energy Surface Acetic anhydride has two non-planar minima, labeled AH[sp,sp] and AH[sp,ac], related to each other by large amplitude motions36 via the transition state TS1 of the rotation of the (O=C−O-C) dihedral angle (cf. Figures 1 & 2). The tittle molecule AH can decompose via two 4membered-ring transition states (TS2 and TS4 from AH[sp,sp] and AH[sp,ac], respectively) and a 6-membered-ring transition state (TS3 from AH[sp,ac]). Figure 2 presented the optimized geometries of the reactants, products and transition states of the title reaction at the B3LYP/CBSB7 level of theory and compared with literature data. The calculated bond lengths and angles are consistent with those of the literature data with maximum errors of ~ 0.1 Å and ~ 1o in bonds and angles, respectively. In particular, the geometries of the TS2, TS3 and TS4 obtained at the B3LYP/CBSB7 level in this study are very similar to the HF/3-21G values reported by Akao et al.5, indicating that the geometries are not sensitive to the calculation methods. The corresponding harmonic vibrational frequencies provided in the Supplementary Table S1 are slightly higher than those available in the literature. (e.g., the average differences are 88 and 49 cm-1 for CH3COOH and CH2CO, respectively).

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

Figure 1. Possible pathways for the unimolecular decomposition of (CH3CO)2O to form CH3COOH and CH2CO. (TS1 is in the conformational process between AH[sp,sp] and AH[sp,ac]; TS2 and TS4 are TSs for the pyrolysis via the 4-membered ring from AH[sp,sp] and AH[sp,ac] conformers, respectively; TS3 occurs in the unimolecular decomposition via the 6membered ring from AH[sp,ac] conformer). Details of the species and transition states involved were presented in Figure 2.

8

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

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

Figure 2. Optimized geometries (distances in Å and angles in degrees) of reactants (AH[sp,sp] and AH[sp,ac]), products (CH3COOH and CH2CO) and transition states (TS1, TS2, TS3 and TS4) at B3LYP/CBSB7 level of theory. The numbers in the parentheses are the reference values from the work of a Atkinson et al.6; b Hellwege et al.37; c Kuchitsu38; d Akao et al.5.

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

Page 10 of 31

Figure 3 presents the energy obtained from single-point energy CCSD(T)/CBS calculations at the optimized B3LYP/CBSB7 geometries (denoted as CCSD(T)/CBS//B3LYP/CBSB7) as extrapolated benchmark values39. The structure of TS3 is similar to those of the elimination of ethyl esters to form carboxylic acids and ethylene40-43; that is, a hydrogen atom in the terminal methyl group migrates to the oxygen atom of carboxyl group [-C(=O)-] to produce acetic acid and ketene. As a result, the pathway through the 6-membered-ring TS is expected to be more favourable than the higher barrier channels via the 4-membered-ring TS2 and TS4. Besides this reaction path, we found three dissociation channels where C-C, C-O and C-H bonds of AH can be broken via barrierless reactions with very high endothermicity (e.g., 89.1, 90.7, and 97.6 kcal/mol above the reactant energy, respectively, cf. Figure S5). Therefore, their contribution to the total rate of decomposition of AH can reasonably be neglected, at least for the kinetic purpose. The energies computed at CCSD(T)/CBS//B3LYP/CBSB7 match very well with those of the other methods except for the calculated data of Akao et al.5 who employed a very small basis set, namely 3-21G (cf. Table 1). Our analysis herein suggests the CCSD(T)/CBS energetic profiles accompanied with the geometry and frequency information from B3LYP/CBSB7 can be considered as the most accurate data available for detailed kinetic analysis. In order to examine the validity of the single-reference electron-correlation CCSD(T) method, T1 diagnostics for the species involved were calculated and presented in the supplementary Table S5.

The T1

diagnostic values of less than 0.02 and 0.04 for closed shell and open-shell species, respectively, suggest the nondynamic correlation energy is small44-48; thus, there is no need to consider highorder methods for the title systems.

10

ACS Paragon Plus Environment

Page 11 of 31

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

Figure 3. ZPE-included potential energy surface of the unimolecular decomposition reaction, (CH3CO)2O → CH3COOH + CH2CO at the CCSD(T)/CBS//B3LYP/CBSB7 level of theory. Values in the parentheses are the reference numbers from the work of Akao et al.5. All numbers are in kcal⋅mol-1. Table 1. Relative energies (in kcal⋅mol-1) of the species involved and the corresponding transition states at different levels of theory at 0 K. The zero-point energy (ZPE) correction was included. AH AH TS1 [sp,sp] [sp,ac] CCSD(T)/CBS[a],[b] 0.0 0.2 0.3 CBS-QB3[a] 0.0 0.3 0.3 CBS-APNO[a] 0.0 0.1 0.8 [a] G3 0.0 0.0 0.3 G4[a] 0.0 0.4 0.3 MP4(SDTQ)//HF[c] 0.0 0.4 [a] This work; [b] CCSD(T)/CBS//B3LYP/CBSB7; [c] The MP4(SDTQ)/3-21G//HF/3-21G. Method

TS2

TS3

TS4

Products

59.0 35.2 58.9 22.6 59.1 36.6 58.8 23.5 34.5 23.6 34.9 23.0 36.3 58.3 22.8 61.3 37.4 60.3 work of Akao et al.5, calculated at

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

Because

the

energy

is

~0.2

kcal·mol-1

higher

than

Page 12 of 31

that

of

AH[sp,sp]

at

CCSD(T)/CBS//B3LYP/CBSB7 level (the experimental value36 is ~ 0.1 kcal⋅mol-1), AH[sp,ac] was included in the hindered internal rotation treatment along the CH3COO−COCH3 axis for AH[sp,sp] (cf. Figure 4). The similar treatment was also used for TS4 in which TS2 was included in the HIR treatment (cf. Supplementary Figure S1). Therefore, the unimolecular decomposition of AH[sp,sp] can proceed via either the 6-membered-ring TS (TS3) or the 4membered-ring one (TS4) with the corresponding HIR treatments, where all possible conformers were implicitly included, as described above.

AH[sp,ac] AH[sp,sp]

AH[sp,sp]

Figure 4. Hindrance potentials for AH[sp,sp] at B3LYP/6-31G(d) level of theory with the step size of 10 degrees. AH[sp,ac] is located at the dihedral angle of 120o or 180o due to the rotation of the C-O bonds.

12

ACS Paragon Plus Environment

Page 13 of 31

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

The calculated thermodynamic properties including heat of formation (∆Hf298) and entropy (S298) for the species involved were compared to the literature data (cf. Table 2) in an attempt to evaluate the reliability of our computed numbers. In general, the calculated values are consistent with the literature data (e.g., when compared to the ATcT data, the differences are less than 0.5 kcal⋅mol-1 and 2.0 cal⋅mol-1⋅K-1 for ∆Hf298 and S298, respectively). Therefore, more confidence has been gained in our calculated data which were derived from a solid framework as presented previously.

Table 2. Comparison of calculated thermodynamic properties of species explicitly considered in the kinetic analysis with literature data (NIST= Webbook NIST, webbook.nist.gov, ATcT = Active Thermochemical Tables49,50,[a], GA = Group Additivity51). Unit: △Hf

298 K

in kcal.mol-1,

S298 K in cal.mol−1·K−1 (1 kcal = 4.184 kJ). Method ∆Hf298 K S298 K Ab initio[b] -139.4 91.5 NIST N/A 93.252 AH[sp,sp] GA -137.1 93.2 [b] Ab initio -103.3 87.0 TS3 Ab initio[b] -80.7 95.9 TS4 Ab initio[b] -103.9 68.3 ATcT[a] -103.3 67.7 CH3COOH 53 NIST -103.5 ± 0.6 67.654 GA -103.4 67.5 -11.4 57.8 Ab initio[b] [a] ATcT -11.6 60.1 CH2CO NIST -11.4 ± 0.455 N/A GA N/A N/A [a] Values collected from Burcat’s online database, http://garfield.chem.elte.hu/Burcat/burcat.html (access date: May 2016); [b] This work at CBS-QB3 level. △Hf was calculated by atomization method. Species

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

Page 14 of 31

Kinetic Analysis The pressure- and temperature-dependent analyses were carried out using the potential energy surface and molecular properties calculated at the CCSD(T)/CBS//B3LYP/CBSB7. In particular, both deterministic and stochastic simulations were independently performed within the ME/RRKM framework to predict the species profiles as well as macroscopic rate coefficients for engineering purposes. Representatively, Figure 5 shows the normalized time-resolved profiles for acetic anhydride and the products (CH3COOH and CH2CO) at the condition of 500 K and 1 atm. The resutls for other conditions can be found in the supplementary material. It is found that the high-temperature chemistry is not so different from the low temperature chemistry, thus it was not discussed explicitly here.

Figure 5. Time-resolved species profiles for the unimolecular decomposition of acetic anhydride at 500 K and 1 atm using both deterministic and stochastic (with different numbers of trials: 500, 1000 and 10000) approaches (1 atm = 0.1 MPa).

14

ACS Paragon Plus Environment

Page 15 of 31

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

Within the deterministic model, our eigenpair analysis showed that a mixing in eigenvalue spectrum, between the fastest chemically significant eigenvalues (CSEs) and the slowest internal energy relaxation eigenvalues (IEREs) where the difference is less than 1 order of magnitude, occurs at T ≥ 1000 K with P = 1.32×10-6 atm (c.f. Figure 6). In such a case, the derived phenomenological rate coefficients using the CSE approach, thus, can only be used with careful evaluation and the stochastic simulation can be alternatively carried out but with a sufficient number of trials. As the pressure increases, the merging regime occurs at a slightly higher temperature (e.g., T ≥ 1200 K at P = 1.32×10-3 atm and T ≥ 1400 K at P = 1.32×10-1 atm, cf. Figure S2); thus the mixing issue is less noticeable at high pressure. IEREs

CSEs

Figure 6. Eigenvalue spectra as a function of temperature for (CH3CO)2CO → CH2CO + CH3COOH (via both TS3 and TS4) at 1.32×10-6 atm. Only the eigenvalues with the smallest magnitude were presented: dashed and solid lines represent chemically significant eigenvalues (CSEs) and internal energy relaxation eigenvalues (IEREs) modes, respectively.

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 31

It is seen in Figure 5 that the stochastic method with a trial number increased from 500 to 10000 is able to predict the profiles which are identical to those obtained from the deterministic framework; therefore, the number of 10000 trials can be concluded as statistically adequate at this condition, and this statistically adequate number is found to increase with the temperature. The computer time was found to increase with the temperature due to the increase of accessible energy levels. Details of the deterministic as well as stochastic simulations at different conditions were presented in Supplementary Figures S3. In short, we have demonstrated for the first time that the two approaches were systematically carried out within the solid ME/RRKM framework to complement each other in an attempt to confidently capture the correct chemistry/physics of this reaction. Rate coefficients k(T, P) as a function of pressure and temperature, calculated for temperature ranging between 300 and 1500 K and for pressures between 0.001 and 100 atm, were depicted in Figures 7 & 8. It can be seen that k(T, P) approaches the high-pressure limit at lower pressure when the temperature decreases. In other words, the fall-off region is shifted to the high pressure with temperature and this region is more noticeable at high temperature regime (e.g., T > 1000 K). The kinetic and thermodynamic data in the Chemkin format were provided in the supplementary Tables S2 & S3 in order to facilitate reactor modeling/simulation using this subkinetic model. As seen in Figure 9, our calculated rate coefficients, k(T, P), are in good agreement (i.e., within the experimental uncertainty) with the experimental data measured by Akao et al. (1996)5, Blake et al. (1971)56 and Blake et al. (1976)57 for both forward and reverse reactions at different pressures as a function of temperatures. The agreement in both directions implicitly reflects the validity of the calculated thermodynamic data for the species involved in this study. Note that the calculated forward rate coefficients also excellently matched with the

16

ACS Paragon Plus Environment

Page 17 of 31

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

older experimental data by Blake 19714 and Szwarc 19513 (cf. Figure 10) and the calculations were carried out using the full PES described in Figure 2 with no adjustment.

Figure 7. Calculated rate coefficients for unimolecular decomposition of (CH3CO)2O as a function of pressure at different temperatures (i.e., 300, 500, 700, 800, 900, 1000, 1200 and 1500 K). The calculations were carried out using the full PES described in Figure 3.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

8 P = 0.001 atm P = 0.01 atm P = 0.1 atm P = 1 atm P = 10 atm P = 100 atm high-P limit

6 4

logk (1/s)

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 31

2 0 -2 -4 -6 -8 -10 -12 -14 0.6

1.0

1.4

1.8 2.2 1000/T (1/K)

2.6

3.0

3.4

Figure 8. Calculated rate coefficients for unimolecular decomposition of (CH3CO)2O as a function of temperature at different pressures (i.e, 0.001, 0.01, 0.1, 1, 10, 100 atm and highpressure limit). The calculations were carried out using the full PES described in Figure 3.

18

ACS Paragon Plus Environment

Page 19 of 31

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

10-2 10-1

10-2

Figure 9. Comparison between calculated and experimental rate coefficients as a function of temperature at different pressures: (a) forward reaction, (CH3CO)2CO → CH2CO + CH3COOH; (b) reverse reaction, CH2CO + CH3COOH → (CH3CO)2CO. Experimental data are from the work of Akao and coworkers (Akao 19965) and Blake and coworkers (Blake 197156 and Blake 197657). 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

(a)

1.0

(CH3CO)2CO → CH2CO + CH3COOH (P =9.21×10-3 - 3.05×10-1 atm & T = 470 - 643 K)

logk (1/s)

0.0 -1.0 -2.0 expt'l (Blake 1971 - fitting) This work (P = 9.21× 10-3 atm)

-3.0

This work (P = 3.05×10-1 atm) -4.0 1.50

1.65

1.80

1.95

2.10

1000/T (1/K)

(b)

5.0

(CH3CO)2CO → CH2CO + CH3COOH (P = 1.68×10-2 - 1.79 ×10-2 atm & T = 553 - 923 K)

4.0 3.0 logk (1/s)

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 31

2.0 1.0 0.0

expt'l (Szwarc 1951 - fitting)

-1.0

This work (P = 1.68×10-2 atm) This work (P = 1.79×10-2 atm)

-2.0 1.00

1.15

1.30

1.45

1.60

1.75

1000/T (1/K) Figure 10. Comparison between calculated and the earlier experimental rate constants (Blake 19714 and Szwarc 19513) for forward reaction, (CH3CO)2CO → CH2CO + CH3COOH as a function of temperature at different pressures. The calculations were carried out using the full PES described in Figure 3.

20

ACS Paragon Plus Environment

Page 21 of 31

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

The branching ratio of the two decomposition channels via two TSs (TS3:TS4) slightly decreases with temperature (e.g., ~ 999 at P = 1 atm and T = 300 – 1500 K, cf. Figure S4 (b)), reflecting the channel via a 6-membered-ring TS (TS3) dominates at the considered conditions. It is worth noting that if these two channels are treated independently (i.e., two separate simulations with each channel), the corresponding ratio is 24 (about 40 times smaller); thus, it is crucial to consider the full PES (as presented in Figure 2) in the kinetic analysis. The result also confirms the hypotheses in the previous experiment5 that the 4-membered-ring TS reaction pathway cannot compete or play any role even at high temperature when comparing with the 6-membered-ring TS pathway. The TS3:TS4 ratio is also found to decrease with pressure (e.g., 999 and 76 at 0.001 and 100 atm, respectively, at 1500 K; cf. Figure S4 (a) & (c)). Sensitivity analyses were also carried out to quantify the dependence of the ME solutions on the calculated ab initio input parameters in the model. Specifically, we examined the effects of the calculated barrier ∆V ‡ of the TS3 and TS4 channels and energy transfer ∆Edown on the timeresolved species mole fraction and the temperature-dependent rate coefficients. The normalized sensitivity coefficients for the species with respect to the barrier change of the TS3 channel (i.e., d(ln[X])/d(ln( ∆ V ‡ (T S 3) ) in Figure 11a) are negative for the reactants and TS3-via products (the products via the TS3 channel) but positive with larger magnitude for the TS4-via products. In other words, the increase of ∆ V ‡ (T S 3) decreased the mole fraction of the reactants and TS3-via products while it increased that of TS4-via products (with a larger normalized effect). Similarly, the increase of ∆ V ‡ (T S 4) basically had no effect on the TS3-via products (the coefficient is close to zero) while it decreased the mole fraction of the reactants and TS4-via products. Since there are two competing channels (i.e., TS3 vs. TS4), the sign as well as the magnitude of the sensitivity coefficients suggest that the TS3 pathway is the dominant channel. The sensitivity 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

Page 22 of 31

analysis was also carried out for the calculated rate constants (cf. Figure 11b). It was observed that the change of the TS3 barrier, ∆∆ V ‡ (T S3) , has an opposite effect on the TS3 rate constants while ∆∆ V ‡ (T S4) had no effect at all. It is interesting to find that both changes caused significant fluctuation of the rate constants for the TS4 channel, the minor channel. This observation can be explained by the competition between excitation/de-excitation and dissociation via different pathways which are time-dependent, thus it makes the kinetics of this system even non-equilibrium. Note that the perturbation has larger effect on the minor channel. It is also partially due to the large uncertainties of the corresponding calculated data which were calculated using the stochastic approach as discussed previously. The sensitivity analysis was also carried out with respected to the energy transfer Edown for the species mole faction as well as rate coefficients (provided in Figure S6). The results are consistent with those observed with the PES analysis about the detailed mechanism. In particular, the change of ∆E has negative effect on the reactant profile (i.e., negative value, cf. Figure S6-a) while it is more complicated for the product channels (both negative and positive, depending on the reaction time). Since the product channel via TS4 is a minor one, the sensitivity coefficient for this channel is larger and more fluctuating than that via TS3 for the temperature range of 3001500 K (cf Figure S6-b).

22

ACS Paragon Plus Environment

Page 23 of 31

(a)

T = 750 K, P = 0.025 atm 3.2E+06 2.8E+06 CH3COOH + CH2CO (via TS4)

d(ln[X])/d(ln∆ ∆V‡)

2.4E+06

+ Solid lines: ∆∆V‡(TS3)= + 1 kcal/mol + Dashed lines: ∆∆V‡(TS4)= + 1 kcal/mol

2.0E+06 1.6E+06

(CH3CO)2O

1.2E+06 8.0E+05

CH3COOH + CH2CO (via TS3)

4.0E+05 0.0E+00 -4.0E+05 0.E+00

2.E-02

(b)

4.E-02 Time (s)

6.E-02

8.E-02

T = 750 - 1000 K, P = 0.025 atm 8.E+05 6.E+05

d(lnk)/d(ln∆ ∆V‡)

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

via TS4

4.E+05 2.E+05

via TS3

0.E+00 -2.E+05 -4.E+05 750

800

850

900

950

1000

T (K) Figure 11. Normalized sensitivity coefficients for the species mole fraction [X] (a) and rate coefficients (b) with respect to the barrier height ∆ V ‡ of the TS3 (solid lines) and TS4 (dashed lines) channels. Note that the analysis was carried out with the full PES described in Figure 3.

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

Page 24 of 31

4. Conclusions In this study, the unimolecular decomposition of acetic anhydride, (CH3CO)2CO → CH2CO + CH3COOH, was investigated by using the high-level quantum chemical and state-ofthe-art integrated deterministic-stochastic statistical rate theory calculations. The results show that the reaction dominantly takes place via a 6-membered-ring transition state over the other two 4-membered-ring TSs in the considered T-P range. The temperature-pressure dependent rate coefficient calculations based on the RRKM/ME solutions were found to be in good agreement with scattered literature experimental data for both directions. The calculated thermodynamic data and rate coefficients, therefore, can be confidently used for various applications in the wide range of conditions (T = 300 – 1500 K and P = 0.001 – 100 atm).

Associated content *Supporting Information 0K

Table S1: The optimized geometries, electronic energies at 0 K ( Eelec ), zero-point energy (ZPE) corrections

and

vibrational

frequencies

of

the

species

involved,

calculated

at

CCSD(T)/CBS//B3LYP/CBSB7 level of theory for the title reaction. Table S2: Parameters obtained from the modified Arrhenius expression of the calculated values for the pressure-dependent rate coefficients over the temperature range of 500 – 1200 K. The calculation was carried out with the channel via TS3. Table S3: Detailed kinetic submechanism in Chemkin format for (CH3CO)2O → CH3COOH + CH2CO. Table S4: Rate constants and uncertainty for (CH3CO)2O → CH3COOH + CH2CO via both TS3 and TS4 at T = 750 – 980 K and P = 0.025 atm, using stochastic approach with 107 trials. Table S5: T1 diagnostic values for the decomposition of acetic anhydride computed at CCSD(T)/cc-pVXZ (X = D, T), using geometries at B3LYP/CBSB7.

24

ACS Paragon Plus Environment

Page 25 of 31

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

Figure S1: Hindrance potentials for AH[sp,sp], TS3, TS4 and CH3COOH at B3LYP/6-31G(d) level of theory. Figure S2: Eigenvalue spectra as a function of temperature for (CH3CO)2CO → CH2CO + CH3COOH (via both TS3 and TS4) at different pressures: (a) 1.32×10-3 atm and (b) 1.32×10-1 atm. Figure S3: Time-resolved species profiles for unimolecular decomposition of acetic anhydride (via TS3) at 1000 K (a) and 1500 K (b) at 1 atm using both stochastic and deterministic (with different numbers of trials: 500, 1000 and 10000) approaches. Figure S4: Branching ratios of acetic anhydride decomposition via two TSs, TS3 and TS4 at 0.001 atm (a) and 100 atm (b). The calculations were carried out using the full PES described in Figure 3. Figure S5: All possible pathways for the thermal decomposition of acetic anhydride calculated at CCSD(T)/CBS//B3LYP/CBSB7 level of theory (@ 0 K and including ZPE corrections). Figure S6: Normalized sensitivity coefficients for the species mole fraction [X] (a) and rate coefficients (b) with respect to the energy transfer ∆ = +100 cm-1. Note that the analysis was carried out with the full PES described in Figure 3. Text S1: Procedure of thermochemical calculations. This material is available free of charge via the Internet at http://pubs.acs.org. Author information Corresponding

Author:

Lam

K.

Huynh.

E-mail

address:

[email protected],

[email protected]. Tel: (84-8) 2211.4046 (ext 3233). Fax: (84-8)3724.4271.

Acknowledgements Computing resources and financial support provided by the Institute for Computational Science and Technology – Ho Chi Minh City and International University, VNU-HCM are gratefully 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

Page 26 of 31

acknowledged. LKH also acknowledges the Department of Science and Technology – Ho Chi Minh City for funding (Grant No. 1166/QĐ-SKHCN). KCL would like to express his gratitude to Ministry of Science and Technology in Taiwan for the financial support under Contract No. MOST 103-2221-E-110-073-MY2.

26

ACS Paragon Plus Environment

Page 27 of 31

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

References (1) Dodson, J. R.; Leite, T.; Pontes, N. S.; Peres Pinto, B.; Mota, C. J., Green acetylation of solketal and glycerol formal by heterogeneous acid catalysts to form a biodiesel fuel additive. ChemSusChem 2014, 7 (9), 2728-2734. (2) Noga, M. J.; Asperger, A.; Silberring, J., N-terminal H3/D3-acetylation for improved highthroughput peptide sequencing by matrix-assisted laser desorption/ionization mass spectrometry with a time-of-flight/time-of-flight analyzer. Rapid. Commun. Mass Spectrom. 2006, 20 (12), 1823-1827. (3) Szwarc, M.; Murawski, J., The kinetics of the thermal decomposition of acetic anhydride. Trans. Faraday Soc. 1951, 47, 269-274. (4) Blake, P. G.; Speis, A., Kinetics of thermal decomposition of acetic anhydride in flow and static systems. J. Chem. Soc. B 1971, 1877-1878. (5) Akao, M.; Saito, K.; Okada, K.; Takahashi, O.; Tabayashi, K., Thermal unimolecular decomposition of acetic anhydride in the gas phase. Ber. Bunsen-Ges. Phys. Chem. 1996, 100 (7), 1237-1241. (6) Atkinson, S. J.; Noble-Eddy, R.; Masters, S. L., Gas-Phase Structures of Ketene and Acetic Acid from Acetic Anhydride Using Very-High-Temperature Gas Electron Diffraction. J. Phys. Chem. A 2016, 120 (12), 2041-2048. (7) Voronova, K.; Mozaffari Easter, C. M.; Torma, K. G.; Bodi, A.; Hemberger, P.; Sztaray, B., Bifurcated dissociative photoionization mechanism of acetic acid anhydride revealed by imaging photoelectron photoion coincidence spectroscopy. Phys. Chem. Chem. Phys. 2016, 18 (36), 25161-25168. (8) Richardson, W. H.; O'Neal, H. E., Comprehensive Chemical Kinetics. Elsevier: Amsterdam, 1972; Vol. 5. (9) Miller, J. A.; Klippenstein, S. J., Master equation methods in gas phase chemical kinetics. J. Phys. Chem. A 2006, 110 (36), 10528-10544. (10) Robertson, S. H.; Pilling, M. J.; Jitariu, L. C.; Hillier, I. H., Master equation methods for multiple well systems: application to the 1-,2-pentyl system. Phys. Chem. Chem. Phys. 2007, 9 (31), 4085-4097. (11) Gillespie, D. T., A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J. Comput. Phys. 1976, 22 (4), 403-434. (12) Gillespie, D. T.; Hellander, A.; Petzold, L. R., Perspective: Stochastic algorithms for chemical kinetics. J. Chem. Phys. 2013, 138 (17), 170901. (13) Becke, A. D., Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98 (7), 5648-5652. (14) Lee, C.; Yang, W.; Parr, R. G., Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37 (2), 785-789. (15) Montgomery-Jr., J. A.; Frisch, M. J.; Ochterski, J. W.; Petersson, G. A., A complete basis set model chemistry. VI. Use of density functional geometries and frequencies. J. Chem. Phys. 1999, 110 (6), 2822-2827. (16) Purvis, G. D.; Bartlett, R. J., A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples. J. Chem. Phys. 1982, 76 (4), 1910-1918.

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 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 31

(17) Scuseria, G. E.; Janssen, C. L.; Schaefer, H. F., An efficient reformulation of the closedshell coupled cluster single and double excitation (CCSD) equations. J. Chem. Phys. 1988, 89 (12), 7382-7387. (18) Scuseria, G. E.; Schaefer, H. F., Is coupled cluster singles and doubles (CCSD) more computationally intensive than quadratic configuration interaction (QCISD)? J. Chem. Phys. 1989, 90 (7), 3700-3703. (19) Masgrau, L.; González-Lafont, A. n.; Lluch, J. M., The reactions CHnD4−n+OH→P and CH4+OD→CH3+HOD as a test of current direct dynamics multicoefficient methods to determine variational transition state rate constants. II. J. Chem. Phys. 2001, 115 (10), 4515-4526. (20) Dunning, T. H., Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90 (2), 1007-1023. (21) Kendall, R. A.; Dunning, T. H.; Harrison, R. J., Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96 (9), 6796-6806. (22) Truhlar, D. G., Basis-set extrapolation. Chem. Phys. Lett. 1998, 294 (1-3), 45-48. (23) Ochterski, J. W.; Petersson, G. A.; Montgomery, J. A., A complete basis set model chemistry. V. Extensions to six or more heavy atoms. J. Chem. Phys. 1996, 104 (7), 2598-2619. (24) Curtiss, L. A.; Raghavachari, K.; Redfern, P. C.; Rassolov, V.; Pople, J. A., Gaussian-3 (G3) theory for molecules containing first and second-row atoms. J. Chem. Phys. 1998, 109 (18), 7764-7776. (25) Curtiss, L. A.; Redfern, P. C.; Raghavachari, K., Gaussian-4 theory. J. Chem. Phys. 2007, 126 (8), 084108. (26) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A., et al. Gaussian 09, Revision A.1, Gaussian, Inc.: Wallingford CT, 2009. (27) Duong, M. v.; Nguyen, H. T.; Truong, N.; Le, T. N. M.; Huynh, L. K., Multi-Species Multi-Channel (MSMC): An Ab Initio-based Parallel Thermodynamic and Kinetic Code for Complex Chemical Systems. Int. J. Chem. Kinet. 2015, 47 (9), 564-575. (28) Mai, T. V. T.; Duong, M. v.; Le, X. T.; Huynh, L. K.; Ratkiewicz, A., Direct ab initio dynamics calculations of thermal rate constants for the CH4 + O2 = CH3 + HO2 reaction. Struct. Chem. 2014, 25 (5), 1495-1503. (29) Eckart, C., The Penetration of a Potential Barrier by Electrons. Phys. Rev. 1930, 35 (11), 1303-1309. (30) Le, T. H. M.; Do, S. T.; Huynh, L. K., Algorithm for auto-generation of hindered internal rotation parameters for complex chemical systems. Comput. Theor. Chem. 2017, 1100, 61-69. (31) Alquaity, A. B.; Giri, B. R.; Lo, J. M.; Farooq, A., High-Temperature Experimental and Theoretical Study of the Unimolecular Dissociation of 1,3,5-Trioxane. J. Phys. Chem. A 2015, 119 (25), 6594-6601. (32) Poling, B. E.; Prausnitz, J. M.; O’Connel, J. P., The Properties of Gases and Liquids. 5th ed., McGraw-Hill Education, New York, Chicago, San Francisco, Athens, London, Madrid, Mexico City, Milan, New Delhi, Singapore, Sydney, Toronto: 2004. (33) Kaiser, E. W.; Wallington, T. J.; Hurley, M. D.; Platz, J.; Curran, H. J.; Pitz, W. J.; Westbrook, C. K., Experimental and Modeling Study of Premixed Atmospheric-Pressure Dimethyl Ether−Air Flames. J. Phys. Chem. A 2000, 104 (35), 8194-8206. (34) Beyer, T.; Swinehart, D. F., Algorithm 448: number of multiply-restricted partitions. Communications of the ACM 1973, 16 (6), 379.

28

ACS Paragon Plus Environment

Page 29 of 31

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

(35) Stein, S. E.; Rabinovitch, B. S., Accurate evaluation of internal energy level sums and densities including anharmonic oscillators and hindered rotors. J. Chem. Phys. 1973, 58 (6), 2438-2445. (36) Wu, G.; Van Alsenoy, C.; Geise, H. J.; Sluyts, E.; Van der Veken, B. J.; Shishkov, I. F.; Khristenko, Acetic Anhydride in the Gas Phase, Studied by Electron Diffraction and Infrared Spectroscopy, Supplemented With ab Initio Calculations of Geometries and Force Fields. J. Phys. Chem. A. 2000, 104 (7), 1576-1587. (37) Hellwege, K.; Hellwege, A., Landolt-Bornstein: Group II: Atomic and Molecular Physics Structure Data of Free Polyatomic Molecules. Springer-Verlag: Berlin, 1976; Vol. 7. (38) Kuchitsu, K., Structure of Free Polyatomic Molecules - Basic Data. Springer: Berlin, 1998. (39) Ranasinghe, D. S.; Petersson, G. A., CCSD(T)/CBS atomic and molecular benchmarks for H through Ar. J. Chem. Phys. 2013, 138 (14), 144104. (40) Saito, K.; Sasaki, T.; Yoshinobu, I.; Imamura, A., Thermal decomposition of ethyl acetate. Branching ratio of the competing paths in the pyrolysis of the produced acetic acid. Chem. Phys. Lett. 1990, 170 (4), 385-388. (41) Saheb, V.; Hosseini, S. M. A., Theoretical studies on the kinetics and mechanism of multichannel gas-phase unimolecular reaction of ethyl acetate. Comput. Theor. Chem. 2013, 1009, 4349. (42) Giri, B. R.; AlAbbad, M.; Farooq, A., High-temperature unimolecular decomposition of ethyl propionate. Chem. Phys. Lett. 2016, 664, 184-190. (43) AlAbbad, M.; Giri, B. R.; Szőri, M.; Farooq, A., On the high-temperature unimolecular decomposition of ethyl levulinate. Proc. Combust. Inst. 2017, 36 (1), 187-193. (44) Lee, T. J.; Taylor, P. R., A diagnostic for determining the quality of single-reference electron correlation methods. Int. J. Quantum Chem. 1989, 36 (S23), 199-207. (45) Jayatilaka, D.; Lee, T. J., Open‐shell coupled‐cluster theory. J. Chem. Phys. 1993, 98 (12), 9734-9747. (46) Rienstra-Kiracofe, J. C.; Allen, W. D.; Schaefer, H. F., The C2H5+ O2 Reaction Mechanism: High-Level ab Initio Characterizations. J. Phys. Chem. A 2000, 104 (44), 98239840. (47) Ali, M. A.; Sonk, J. A.; Barker, J. R., Predicted Chemical Activation Rate Constants for HO2 + CH2NH: The Dominant Role of a Hydrogen-Bonded Pre-reactive Complex. J. Phys. Chem. A 2016, 120 (36), 7060-7070. (48) Monge-Palacios, M.; Rafatijo, H., On the role of the termolecular reactions 2O2 + H2 --> 2HO2 and 2O2 + H2 --> H + HO2 + O2 in formation of the first radicals in hydrogen combustion: ab initio predictions of energy barriers. Phys. Chem. Chem. Phys. 2017, 19 (3), 2175-2185. (49) Ruscic, B.; Pinzon, R. E.; Laszewski, G. v.; Kodeboyina, D.; Burcat, A.; Leahy, D.; Montoya, D.; Wagner, A. F., J. Phys.: Conf. Ser. 2005, 16, 561-570. (50) Ruscic, B.; Pinzon, R. E.; Morton, M. L.; Srinivasan, N. K.; Su, M. C.; Sutherland, J. W.; Michael, J. V., Active Thermochemical Tables: accurate enthalpy of formation of hydroperoxyl radical, HO2. J. Phys. Chem. A 2006, 110 (21), 6592-6601. (51) Cohen, N.; Benson, S. W., Estimation of heats of formation of organic compounds by additivity methods. Chem. Rev. 1993, 93 (7), 2419-2438. (52) Stull D.R., J., The Chemical Thermodynamics of Organic Compounds. Wiley: New York, 1969. (53) Webbook NIST. The heat of formation of acetic acid, CH3COOH. http://webbook.nist.gov/cgi/cbook.cgi?ID=C64197&Units=CAL&Mask=1 (accessed July 2016). 29

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 30 of 31

(54) Weltner, W., The Vibrational Spectrum, Associative and Thermodynamic Properties of Acetic Acid Vapor. J. Am. Chem. Soc. 1955, 77 (15), 3941-3950. (55) Nuttall, R. L.; Laufer, A. H.; Kilday, M. V., The enthalpy of formation of ketene. J. Chem. Thermodyn. 1971, 3 (2), 167-174. (56) Blake, P. G.; Davies, H. H., Reactions of keten. Part I. Kinetics of the gas-phase reaction with acetic acid. J. Chem. Soc. B 1971, 1727-1728. (57) Blake, P. G.; Vayjooee, M. B., Reactions of keten. Part V. Gas-phase reactions with carboxylic acids. Effect of alkyl substituents. J. Chem. Soc., Perkin Trans. 2 1976, (9), 988-990.

30

ACS Paragon Plus Environment

Page 31 of 31

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

31

ACS Paragon Plus Environment