ARTICLE pubs.acs.org/JPCA
Atmospheric Chemistry of Two Biodiesel Model Compounds: Methyl Propionate and Ethyl Acetate Vibeke F. Andersen,† Tesfaye A. Berhanu,† Elna J. K. Nilsson,† Solvejg Jørgensen,† Ole John Nielsen,† Timothy J. Wallington,‡ and Matthew S. Johnson*,† †
Copenhagen Center for Atmospheric Research (CCAR), Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 København Ø, Denmark ‡ Systems Analytics and Environmental Sciences Department, Research and Advanced Engineering, Ford Motor Company, Mail Drop RIC-2122, Dearborn, Michigan 48121-2053, United States
bS Supporting Information ABSTRACT: The atmospheric chemistry of two C4H8O2 isomers (methyl propionate and ethyl acetate) was investigated. With relative rate techniques in 980 mbar of air at 293 K the following rate constants were determined: k(C2H5C(O)OCH3 + Cl) = (1.57 ( 0.23) 1011, k(C2H5C(O)OCH3 + OH) = (9.25 ( 1.27) 1013, k(CH3C(O)OC2H5 + Cl) = (1.76 ( 0.22) 1011, and k(CH3C(O)OC2H5 + OH) = (1.54 ( 0.22) 1012 cm3 molecule1 s1. The chlorine atom initiated oxidation of methyl propionate in 930 mbar of N2/O2 diluent (with, and without, NOx) gave methyl pyruvate, propionic acid, acetaldehyde, formic acid, and formaldehyde as products. In experiments conducted in N2 diluent the formation of CH3CHClC(O)OCH3 and CH3CCl2C(O)OCH3 was observed. From the observed product yields we conclude that the branching ratios for reaction of chlorine atoms with the CH3, CH2, and OCH3 groups are 9 ( 2%, respectively. The chlorine atom initiated oxidation of ethyl acetate in N2/O2 diluent gave acetic acid, acetic acid anhydride, acetic formic anhydride, formaldehyde, and, in the presence of NOx, PAN. From the yield of these products we conclude that at least 41 ( 6% of the reaction of chlorine atoms with ethyl acetate occurs at the CH2 group. The rate constants and branching ratios for reactions of OH radicals with methyl propionate and ethyl acetate were investigated theoretically using transition state theory. The stationary points along the oxidation pathways were optimized at the CCSD(T)/cc-pVTZ//BHandHLYP/aug-cc-pVTZ level of theory. The reaction of OH radicals with ethyl acetate was computed to occur essentially exclusively (∼99%) at the CH2 group. In contrast, both methyl groups and the CH2 group contribute appreciably in the reaction of OH with methyl propionate. Decomposition via the R-ester rearrangement (to give C2H5C(O)OH and a HCO radical) and reaction with O2 (to give CH3CH2C(O)OC(O)H) are competing atmospheric fates of the alkoxy radical CH3CH2C(O)OCH2O. Chemical activation of CH3CH2C(O)OCH2O radicals formed in the reaction of the corresponding peroxy radical with NO favors the R-ester rearrangement.
1. INTRODUCTION Concerns related to climate change and energy security combined with a desire to provide support for rural communities has led to a substantial increase in the global production of biofuels over the past decade. Ethanol produced from corn and sugar cane is the dominant biofuel in the US and in Brazil, respectively, and is typically blended into gasoline before use (e.g., in E10, E25, E85). Fatty acid methyl esters (FAMEs) made from transesterification of plant oils or animal fats are the dominant biofuel in the EU and are typically blended into petroleum diesel before use (e.g., in B2, B20). FAME is usually also referred to as “biodiesel” as distinct from hydrotreated vegetable oils which are referred to as “renewable diesel”. The r 2011 American Chemical Society
majority of global biodiesel is produced in Europe, with rapeseed and sunflower being the two major feedstocks.13 The widespread use of FAMEs in diesel blends will lead to their release into the atmosphere. The atmospheric oxidation of esters is initiated by reaction with OH radicals and, to a lesser extent, chlorine atoms. Detailed kinetic and mechanistic data concerning these reactions are needed as inputs for atmospheric chemistry models to assess the environmental impacts associated with release of esters into the air. The available database concerning the atmospheric chemistry of esters is Received: May 24, 2011 Revised: July 9, 2011 Published: July 11, 2011 8906
dx.doi.org/10.1021/jp204819d | J. Phys. Chem. A 2011, 115, 8906–8919
The Journal of Physical Chemistry A
ARTICLE
Table 1. Rate Constants (Units of 1012 cm3 molecule1 s1) for Reactions of Chlorine Atoms and OH Radicals with Methyl Propionate and Ethyl Acetate k
reaction (1) C2H5C(O)OCH3 + Cl f products
method
reference
15.7 ( 2.3
relative rate 293 ( 0.5 K, 980 mbar air
this study
19.8 ( 2.6 15.1 ( 2.2
absolute rate 298 ( 2 K, 19.578 mbar He relative rate 296 ( 2 K, 990 mbar air
Notario et al.6 Cavalli et al.7
0.925 ( 0.127
relative rate 293 ( 0.5 K, 980 mbar air
this study
0.38
ab initioa Transition state theory, 298 K
this study
0.929 ( 0.113
relative rate 296 K, 990 mbar air
Cavalli et al.7
1.03 ( 0.04
absolute rate 296 K, 3370 mbar Ar
Wallington et al.8
0.83 ( 0.09
absolute rate 298 K, 133 mbar He
Le Calve et al.9
(3) CH3CHClC(O)OCH3 + Cl f products
2.71 ( 0.39
relative rate 293 ( 0.5 K, 980 mbar N2
this study
(4) CH3C(O)OC2H5 + Cl f products
17.6 ( 2.2 20.1 ( 2.5
relative rate 293 ( 0.5 K, 980 mbar air absolute rate 298 ( 2 K, 2080 mbar air
this study Notario et al.6
13.7 ( 2.0
absolute rate 264.5380 K, 26270 mbar He
Cuevas et al.10
17.6 ( 1.1
absolute rate 298 ( 2 K, 3930 mbar N2
Xing et al.11
17.6 ( 2.6
relative rate 296 ( 1 K, 3933 mbar N2
Xing et al.11
1.54 ( 0.22
relative rate 293 ( 0.5 K, 980 mbar air
this study
2.00
ab initioa transition state theory, 298 K
this study
1.51 ( 0.14
absolute rate 296 K, 3370 mbar Ar
Wallington et al.8
1.67 ( 0.22 1.73 ( 0.20
absolute rate 298 K, 1133 mbar He relative rate 298 ( 4 K, 1013 mbar air
El Boudali et al.13 Picquet-Varrault et al.14
(2) C2H5C(O)OCH3 + OH f products
(5) CH3C(O)OC2H5 + OH f products
a
CCSD(T)/cc-pVTZ//BHandHLYP/aug-cc-pVTZ.
limited,4,5 and the present work was undertaken to improve our understanding of the chemistry of this class of oxygenated organic compound. The FAMEs present in biodiesel (e.g., methyl oleate, methyl linoleate) have low vapor pressures and are difficult to study in smog chamber systems. In the present study we investigate, using experimental and computational methods, the chemistry of two small esters as models for the larger esters present in biodiesel. We conducted an experimental investigation of the kinetics and products of the OH radical and chlorine atom initiated oxidation of methyl propionate and ethyl acetate in 930980 mbar of N2/O2 diluent with, and without, added NO at 293296 K. In addition we conducted a computational study of the kinetics of reactions of OH radicals with methyl propionate and ethyl acetate. Two previous studies of k(Cl + methyl propionate)6,7 and three studies of k(OH + methyl propionate)79 have been reported. A study of the chlorine atom initiated oxidation of methyl propionate in 740 Torr air, at 296 K in the presence of NO was performed by Cavalli et al.7 There are no reported studies of the oxygen and NO dependence of the oxidation mechanism of methyl propionate. No theoretical studies of the OH-initiated oxidation of methyl propionate have been reported. Three previous studies of k(Cl + ethyl acetate)6,10,11 and three studies of k(OH + ethyl acetate)8,12,13 have been reported. The previous kinetic studies are summarized in Table 1. The products of the OH radical initiated oxidation of ethyl acetate in air at ambient temperature and pressure in the presence of NO have been studied by Tuazon et al.14 and Picquet-Varrault et al.15 The oxygen and temperature dependence of the products of chlorine initiated oxidation of ethyl acetate have been described by Orlando and Tyndall.16 The present work is the first experimental study of the NOx dependence of the product branching ratios of the chlorine atom initiated oxidation of ethyl acetate and the first theoretical study of the kinetics of the reactions of OH radicals with ethyl acetate.
2. EXPERIMENTAL SETUP Experiments were performed in the photochemical reactors at CCAR in Copenhagen17 and at Ford Motor Company in Dearborn, MI.18 The system in Copenhagen consists of a 100-L quartz reaction chamber sealed at both ends with stainless steel flanges and surrounded by 8 UV-A (325380 nm), 16 UV-C (254 nm), and 12 broad band sun lamps (325 nm to infrared). The system is enclosed in an insulating box equipped with a temperature control system to ensure stable (within 0.5 K) temperatures during experiments. Gas mixtures in the chamber were monitored with a Bruker IFS 66v/s FTIR spectrometer using an IR path length of 72 m and resolution of 0.125 cm1. The system at Ford Motor Company consists of a 140-L Pyrex reactor interfaced to a Mattson Sirius 100 FTIR spectrometer. The reactor is surrounded by 22 fluorescent blacklamps (GE F40T12BLB) which were used to initiate the experiments. Concentrations of reactants and products were monitored by FTIR spectroscopy. IR spectra were derived from 32 coadded interferograms with a spectral resolution of 0.25 cm1 and an analytical path length of 27.1 m. The experiments were performed at 296 ( 2 K. 3. MATERIALS AND METHOD Relative rate methods were used to measure rate constants for the following reactions
8907
C2 H5 CðOÞOCH3 þ Cl f products
ð1Þ
C2 H5 CðOÞOCH3 þ OH f products
ð2Þ
CH3 CHClCðOÞOCH3 þ Cl f products
ð3Þ
CH3 CðOÞOC2 H5 þ Cl f products
ð4Þ
CH3 CðOÞOC2 H5 þ OH f products
ð5Þ
dx.doi.org/10.1021/jp204819d |J. Phys. Chem. A 2011, 115, 8906–8919
The Journal of Physical Chemistry A
ARTICLE
Table 2. Reference Rate Constants Used in This Work rate constant (cm3 molecule1 s1)
reference
(6) C2H5Cl + Cl f products
8.04 1012
37
(7) C2H4 + Cl f products
9.29 1011
38
reaction
11
(8) C2H6 + Cl f products
5.9 10
(9) C3H8 + Cl f products
1.40 1010
39
(10) C2H6 + OH f products
2.40 1013
39
(11) c-C6H12 + OH f products
6.37 1012
40
(12) C3H8 + OH f products
1.10 1012
39
39
Reactions 612 were employed as the reference reactions. The reference rate constants used are listed in Table 2. C2 H5 Cl þ Cl f products
ð6Þ
C2 H4 þ Cl f products
ð7Þ
C2 H6 þ Cl f products
ð8Þ
C3 H8 þ Cl f products
ð9Þ
C2 H6 þ OH f products
ð10Þ
c-C6 H12 þ OH f products
ð11Þ
C3 H8 þ OH f products
ð12Þ
k1 was determined relative to k6 and k7, k2 was determined relative to k10 and k11, k3 was determined relative to k6 and k8, k4 was determined relative to k6, k7, and k9, and k5 was determined relative to k11 and k12. All chemicals were obtained from commercial sources: C2H5C(O)OCH3 (g99.0%, Fluka), CH3CHClC(O)OCH3 (97%, Aldrich), CH3C(O)OC2H5 (99.8%, Sigma-Aldrich), C2H5Cl (g99.7%, Gerling Holz + CO), C2H6 (g99.95%, Fluka), C2H4 (g99.5%, Aldrich), cC6H12 (99.5%, Aldrich), and C3H8 (98%, Aldrich). Relative rate experiments were performed in 980 mbar of air, or N2, diluent. Product studies were performed in 930 mbar of N2/O2 diluent with O2 partial pressures of 70930 mbar in the presence and absence of NO. The photochemical chain chlorination of methyl propionate was studied in 980 mbar of N2 diluent. Experimental conditions are shown in Tables S1 and S2 in the Supporting Information. Chlorine atoms were generated by photolysis of Cl2 in mixtures containing methyl propionate or ethyl acetate and reference compounds. OH radicals were generated by the photolysis of O3 in the presence of H2 O. O3 was generated from O2 with an ozone generator from O3 Technology and condensed on silica gel that was cooled with ethanol and dry ice to approximately 70 °C. The initial concentration of ozone and water vapor was approximately 5 times the total concentration of reactants. Photolysis intervals were typically between 5 s and 10 min, and an IR spectrum was recorded after each photolysis step. The IR absorption spectra of methyl propionate and ethyl acetate are provided in the Supporting Information. The reaction rates were determined using the relative rate method. Assuming that the reactant and reference compounds are only lost via the reaction of interest and that neither the reactant nor the reference
is formed in the system, the following relation is valid ! ! ½reference0 ½reactant0 kreference ln ¼ ln ½referencet kreactant ½reactantt
ðIÞ
where [reactant]0 and [reference]0 are the initial concentrations of reactant and reference, respectively, and [reactant]t and [reference]t are the concentrations of reactant and reference at time t. Plotting the loss of reactant versus reference compound gives a straight line with slope kreactant/kreference. Unless stated otherwise, all quoted uncertainties are 2 standard deviations from least-squares regressions and include uncertainties in the analysis of the IR spectra (typically 510%).
4. COMPUTATIONAL DETAILS The methodology used in the present work was the same as that employed in our recent study of the reaction of OH radicals with methyl acetate.19 The geometries of the stationary points along the reaction path were optimized using density functional theory BHandHLYP20,21 with the basis set aug-cc-pVTZ.22 The calculated frequencies were used to characterize the stationary points (a minimum or a saddle point) and to estimate the zero point vibrational energy (ZPE). A single point energy optimization with the coupled cluster single and double excitation method including a perturbative estimate of triples (CCSD(T))23 and the cc-pVTZ basis set was used to obtain reliable energies at each stationary point. All calculations were performed using Gaussian 03.24 One problem with applying unrestricted calculations to these systems is the potential for spin contamination from higher lying spin states. The expectation value of the total spin was less than 0.78 before and 0.75 after spin annihilation. Therefore, spin contamination is not considered to be severe. In general, spin contamination for the transition state structures and the product complexes was larger than that for the reactant complexes. We observe prereaction complexes with lower energy than the separated reactants. If we assume a fast equilibrium between the reactant and the prereactive complex, we can calculate the rate constant at 298 K using standard transition state theory (TST)25 kb T QTS Ea k ¼ kðTÞσ exp ðIIÞ h QR QOH kb T where k(T) is the quantum tunneling coefficient, σ the symmetry factor, i.e., the number of identical reaction paths, kb Boltzmann’s constant, h Planck’s constant, T temperature, Q the total partition function, and finally Ea is the activation energy at 0 K calculated as Ea = ETS (ER + EOH). The quantum tunneling coefficient was considered using the one-dimensional Wigner transmission coefficient26 !2 1 hυq kðTÞ ¼ 1 þ ðIIIÞ 24 kb T where υq is the imaginary frequency associated with the reaction coordinate.
5. RELATIVE RATE RESULTS To test for photolysis, mixtures containing methyl propionate, ethyl acetate and the reference compounds (but not Cl2 or O3/H2O) in air were irradiated using conditions (number of lamps and time span) similar to those used in the relative rate experiments. To test for 8908
dx.doi.org/10.1021/jp204819d |J. Phys. Chem. A 2011, 115, 8906–8919
The Journal of Physical Chemistry A
ARTICLE
Figure 1. Relative rate plots for the reaction of methyl propionate with Cl (left) and OH (right) in 980 mbar air at 293 ( 0.5 K.
Table 3. Results for Relative Rate Studies of the Cl and OH Initiated Oxidations of Methyl Propionate, Methyl-2-chloropropionate and Ethyl Acetatea rate constant reaction
a
kester/kreference
(1012 cm3 molecule1 s1)
(1) relative to (6)
1.81 ( 0.43
14.6 ( 3.4
(1) relative to (7)
0.18 ( 0.03
16.7 ( 3.1
(2) relative to (10)
3.65 ( 0.72
0.876 ( 0.174
(2) relative to (11)
0.15 ( 0.03
0.974 ( 0.185
(3) relative to (6) (3) relative to (8)
0.328 ( 0.064 0.047 ( 0.010
2.64 ( 0.52 2.78 ( 0.59
(4) relative to (6)
1.93 ( 0.38
15.5 ( 3.1
(4) relative to (7)
0.184 ( 0.036
17.1 ( 3.3
(4) relative to (9)
0.145 ( 0.033
20.3 ( 4.6
(5) relative to (11)
0.234 ( 0.044
1.49 ( 0.28
(5) relative to (12)
1.43 ( 0.30
1.58 ( 0.33
Experiments were performed at 293 ( 0.5 K.
unwanted heterogeneous reactions, mixtures (including Cl2 or O3/H2O) were left for 1030 min in the chamber in the dark. There was no discernible loss (