Atmospheric Oxidation Mechanism and Kinetic Studies for OH and NO₃ Radical-Initiated Reaction of Methyl Methacrylate

Abstract

The mechanism for OH and NO₃ radical-initiated oxidation reactions of methyl methacrylate (MMA) was investigated by using density functional theory (DFT) molecular orbital theory. Geometrical parameters of the reactants, intermediates, transition states, and products were fully optimized at the B3LYP/6-31G(d,p) level. Detailed oxidation pathways were presented and discussed. The rate constants were deduced by the canonical variational transition-state (CVT) theory with the small-curvature tunneling (SCT) correction and the multichannel Rice-Ramspergere-Kassele-Marcus (RRKM) theory, based on the potential energy surface profiles over the general atmospheric temperature range of 180–370 K. The calculated results were in reasonable agreement with experimental measurement.

1. Introduction

Methyl methacrylate (MMA, CH₃COOCH₂CH₃) is widely used in manufacture of resins and plastics [1–3]. Currently, over 2.04 × 10⁹ kg of MMA has been produced by industrial processes [4]. Since the increasing use of MMA, the emission into the atmosphere may greatly rise. The high vapor pressure (3.9 × 10³ Pa at 293 K) indicates that MMA exists mainly in the gas phase under the general atmospheric conditions [5,6]. Neurological symptoms have been reported in humans following acute exposure to methyl methacrylate [7,8]. Fetal abnormalities have been reported in animals exposed to methyl methacrylate by injection and inhalation [9]. Once released into the atmosphere, the unsaturated MMA may be oxidized by OH radicals during daytime, nitrate radicals (NO₃), and ozone molecules during the nighttime [10,11]. The most important oxidation degradation of MMA is initiated by reaction with OH radicals in the atmosphere. However, as the formation of OH mainly takes place in daytime, OH concentration rapidly decreases after sunset. As a consequence, reactions with OH radicals only occur during the day. NO₃ radical undergoes rapid photolysis upon absorption of radiation. During daytime, concentration of NO₃ radicals is very low. Therefore, daytime NO₃ chemistry is expected to be unimportant for the MMA. However, NO₃ has been identified and measured in the nighttime atmosphere, especially in some seriously polluted regions. The reactions with NO₃ will dominate the loss pathway of MMA during nighttime. Moreover, in certain locations during certain times of the year, reactions with Cl atoms may also be important. These atmospheric oxidation reactions may contribute to the formation of the secondary photooxidants and aerosols in the troposphere.

The large emission, high volatility, and toxicity of MMA make it a potential important source of environmental concern in the atmosphere. The previous studies on the atmospheric reaction of MMA mainly focus on obtaining kinetic parameters. Low pressure rate constant for the reaction of MMA with OH radicals has been determined using the discharge flow-laser induced fluorescence (DF/LIF) technique at the room temperature [12]. The rate constants for the reactions of MMA with OH radicals and Cl atoms were carried out using the relative rate methods at the ambient temperature and pressure conditions [13], giving k_OH = (4.15 ± 0.32) × 10⁻¹¹ cm³ molecule⁻¹ s⁻¹, and k_Cl = (2.82 ± 0.93) × 10⁻¹⁰ cm³ molecule⁻¹ s⁻¹. Arrhenius formula over the temperature range of 253–374 K was derived as k(T)(MMA + OH) = (2.5 ± 0.8 × 10⁻¹²) exp((825 ± 55)/T) cm³ molecule⁻¹ s⁻¹. The rate constants for the gas-phase reactions of a series of acrylate methyl (including MMA) with O₃ were determined using smog chamber techniques [14] at 1.01 × 10⁴ Pa and 294 ± 2 K. The rate constant of MMA was obtained as k_O3 = (6.7 ± 0.9) × 10⁻¹⁸ cm³ molecule⁻¹ s⁻¹. Recently, the reactions of MMA with OH radicals in the presence of NOx were studied in environmental chamber [15], revealing that the primary products are methyl pyruvate (92% ± 16%) and methanal (87% ± 12%). However, due to the scarcity of efficient detection schemes for radical intermediate species and commercially available standards, the atmospheric reaction mechanism of MMA is still unclear. In this work, quantum chemical and direct kinetic calculations were performed to elucidate the reaction mechanism of the OH and NO₃ radical-initiated atmospheric oxidation of MMA.

2. Results and Discussion

To confirm the reliability of calculated results, the geometries and vibrational frequencies of CH₃COOCH₂CH₃, CH₃CHO, and CH₃CH=CH₂ were calculated at the B3LYP/6-31G(d,p) level. The results are in reasonable accordance with the corresponding experimental values and the discrepancy remains within 1% for geometrical parameters and 4% for vibrational frequencies [16,17].

2.1. Reaction Mechanism

The reactions of MMA with OH and NO₃ radicals involve two kinds of pathways: H abstraction from MMA and addition of OH or NO₃ radical to the C=C bond. The reaction schemes embedded with the potential barriers and reaction enthalpies are depicted in Figures 1 and 2. For convenience of description, the C and H atoms in MMA are labeled, as depicted in Figures 1–6.

2.1.1. H Abstraction by OH and NO₃ Radicals

As shown in Figures 1 and 2, there are eight H atoms in the MMA molecule. Since the C²–C⁵ and O–C⁴ bonds are rotatable, three H atoms bonded to C⁴ are equivalent, and three H atoms bonded C⁵ are equivalent. So, there are four kinds of H atoms in the MMA molecule, indicating that four possible pathways can be identified: H abstraction from the C¹–H¹, C¹–H², C⁴–H, and C⁵–H bonds. The products of H abstraction by OH or NO₃ radicals are IM1–1, IM1–2, IM1–3, and IM1–4. Eight transition states were located (TS1, TS2, TS3, and TS4 for H abstraction by OH radicals and TS1′, TS2′, TS3′, and TS4′ for H abstraction by NO₃ radicals), and each one was identified with only one imaginary frequency. The reaction schemes embedded with the potential barriers and reaction heats are depicted in Figures 1 and 2. All the H abstraction processes are exothermic. The potential barriers of pathways 3 and 4 are lower than that of pathways 1 and 2. In addition, pathways 3 and 4 are more exothermic than pathways 1 and 2. Thus, the thermodynamically favored H abstraction pathways are pathways 3 and 4. The resulting IM1–3 and IM1–4 are the main H abstraction products. Similarly, pathways 7 and 8 are thermodynamically favored H abstraction pathways for the reaction of MMA with NO₃ radicals. Two same products, IM1–3 and IM1–4, can be yielded from pathways 7 and 8. Therefore, IM1–3 and IM1–4 are important intermediates produced in the oxidation process of MMA initiated by OH and NO₃ radicals.

IM1–3 and IM1–4 are activated radicals and can further react with the ubiquitous oxygen molecules in the atmosphere to form two organic peroxy radicals, IM1–5 and IM1–8. To evaluate the nature of the entrance channel for the formation of the organic peroxy radicals of IM1–5 and IM1–8, we examined the potential along the reaction coordinate, especially to determine whether there is a well-defined transition state or if the process proceeds via a loose transition state without a barrier, depicted in Figure 3. The profiles of the potential energy surface were scanned by varying the newly formed C⁴–O or C⁵–O bond. We found no energy exceeding the C⁴–O or C⁵–O bond dissociation threshold along the reaction coordinate. This shows that the reactions of IM1–3 and IM1–4 with O₂ proceed via a barrierless association. The processes are strongly exothermic by 32.19 and 18.29 kcal/mol.

In the troposphere, IM1–5 or IM1–8 will react immediately with ubiquitous NO. The entrance channel of the reactions is exoergic, leading to two vibrationally excited intermediates (denoted as IM1–6 and IM1–9), which promptly react via unimolecular decomposition. Unimolecular decompositions of IM1–6 and IM1–9 occur via cleavage of the O-O bond, forming NO₂ and two alkoxy radicals IM1–7 and IM1–10. The two unimolecular decomposition processes are endothermic with high potential barrier. IM1–7 can further react with O₂, followed by loss of HO₂, to form CH₂=C(CHO)COOCH₃ (P1) and HO₂, or via unimolecular decomposition and H migration to form CH₂=C(CH₃)C(O)OC(O)OH (P2). The decomposition of IM1–10 has a high barrier of 24.51 kcal/mol and cannot occur under the general atmospheric conditions. The reaction with O₂ is the primary removal pathway for IM1–10, and the products are CH₂=C(CHO)COOCH₃ (P3) and HO₂.

2.1.2. Addition of OH to MMA

Two carbon atoms in the C=C bond of MMA are not equivalent, thus, OH radicals can attack C¹ or C² atom to form two different adducts, IM1a and IM1b, as shown in Figure 1. The OH-MMA adducts, IM1a and IM1b, could easily further react with O₂/NO in the atmosphere and form intermediates IM1d and IM1h, respectively. As shown in Figure 4, the two addition steps are barrierless and strongly exothermic. Then, the energy-rich intermediates IM1d and IM1h can decompose to yield IM1e and IM1i. The two decomposition processes are endothermic with high potential barriers. Three possible unimolecular decomposition pathways from IM1e were identified. The first pathway occurs via H migration and breaking of C²–C⁵ bond to yield P3. The barrier height is 10.22 kcal/mol. The second and third decomposition pathways are cleavage of C¹–C² and C²–C³ bonds, with barriers of 17.07 and 15.86 kcal/mol. The results show that the first decomposition pathway is energetically more favorable. Unimolecular decomposition of IM1i only occurs via cleavage of C¹–C² bond to yield IM1k. In addition, IM1i can further react with O₂, followed by the loss of HO₂, to yield CH₃OC(O)C(CH₃)(OH)CHO (P6). The decomposition of IM1i needs to overcome a barrier of 19.23 kcal/mol, therefore, the reaction with O₂ is relatively more favorable under the general atmospheric conditions. The calculated results suggest that the main products from secondary reactions of OH-MMA adducts are P3, HO₂ and methanal.

2.1.3. Addition of NO₃ to MMA

Similar to OH addition to MMA, the addition of NO₃ to MMA can form two different NO₃-MMA adducts, IM2a and IM2b, as shown in Figure 2. The scanned profile of the potential energy surface shows that the two addition reactions are barrierless. The NO₃-MMA adducts can further react via decomposition or isomerization, as depicted in Figure 5, or be removed via reactions with O₂/NO, as depicted in Figure 6. As shown in Figure 5, in the decomposition processes, a three-membered ring product, P7, was generated from the intermediates IM2a and IM2b via loss of NO₂, with potential barriers of 19.89 and 15.61 kcal/mol; and a five-membered ring product, P8, was formed from the isomerization of IM2a and IM2b, with potential barriers of 21.60 and 20.41 kcal/mol. The results suggest that the decomposition and isomerization processes are competitive.

In the oxygen-rich atmosphere, the NO₃-MMA adducts could further react with O₂ to yield organic peroxy radicals. The detailed subsequent reactions are presented in Figure 6. The reactions of IM2a and IM2b with O₂ are barrierless and strongly exothermic by 49.19 and 61.94 kcal/mol to form IM2c and IM2f. In the troposphere, the two intermediates, IM2c and IM2f, will further react with ubiquitous NO immediately and yield two vibrationally excited intermediates, IM2d and IM2g, respectively. Unimolecular decompositions of IM2d and IM2g result in the formation of IM2e and IM2h via loss of NO₂. The alkoxy radical IM2e can further react via three possible unimolecular decomposition pathways, as depicted in Figure 6. Comparison of potential barriers of the three pathways shows that cleavage of the C¹–C² and C²–C³ bonds are energetically favorable, leading to the products of P3, P11 and methanal. IM2j produced from the decomposition of IM2h can further decompose to yield P3 via the loss of NO₂. This process is exothermic by 47.21 kcal/mol. IM2j can also further react with O₂, followed by the loss of HO₂ to produce P11. Calculations show that the reaction of IM2j with O₂ is energetically more favorable than the decomposition of IM2j. Comparison of the three kinds of the reaction pathways for NO₃-MMA adducts suggests that the reaction of the NO₃-MMA adduct with O₂/NO is thermodynamically and energetically more favorable than their unimolecular decomposition or isomerization.

2.2. Rate Constant Calculations

On the basis of the profile of the potential energy surface calculated by the CCSD(T)/6-31G(d) + CF//B3LYP/6-31G(d,p) method, the individual and overall rate constants for the OH and NO₃ radical-initiated reactions of MMA were calculated over the temperature range from 180 to 370 K, which is the typical temperature range of the troposphere. The rate constants of H abstraction from MMA were calculated by using the CVT/SCT method. The individual rate constants for the H abstraction pathways 1–8 are noted as k_abs¹–k_abs⁸, respectively. For the addition pathways, the rate constants were calculated by using the multichannel RRKM method. The rate constants for addition of OH to the C¹ and C² atoms of MMA are noted as k_OH¹ and k_OH², respectively; and the rate constants for addition of NO₃ to the C¹ and C² atoms are noted as k_NO3¹ and k_NO3², respectively. The overall rate constant for the reaction of MMA with OH radical is noted as k_OH, k_OH = k_OH¹ + k_OH² + k_abs¹ + k_abs² + k_abs³ + k_abs⁴; and the overall rate constant for the reaction with NO₃ radical is noted as k_NO3, k_NO3 = k_NO3¹ + k_NO3² + k_abs⁵ + k_abs⁶ + k_abs⁷ + k_abs⁸. Comparison of the rate constants of the H abstraction pathways and addition pathways shows that the H abstraction pathways can be negligible for their little contribution under the general atmospheric conditions.

At 298 K, the calculated overall rate constants for the reactions of MMA with OH and NO₃ radicals are 4.36 × 10⁻¹¹ and 3.64 × 10⁻¹⁵ cm³ molecule⁻¹ s⁻¹, excellently agreed with the available experimental values [13,18,19]. The good agreement infers that the calculated other rate constants are reasonable. The calculated overall rate constants are fitted over the temperature range of 180–370 K, and Arrhenius formulas are given in units of cm³ molecule⁻¹ s⁻¹:

$$k(T)(\text{MMA}+\text{OH})=(1.83\times {10}^{-12})\text{exp}(945.58/T)$$

(1)

$$k(T)(\text{MMA}+{\text{NO}}_3)=(6.75\times {10}^{-16})\text{exp}(502.48/T)$$

(2)

To assess the impact to the environment, it is critical to known the atmospheric lifetime of MMA. The global average concentration OH radical (c_OH) in daytime is 2 × 10⁶ molecule/cm³ [20], and the typical concentration of NO₃ radicals (c_NO3) in the continental boundary layer is 5 × 10⁸ molecule/cm³ [21]. According to the rate constants of the reactions of MMA with OH and NO₃ radicals, using the expression:

$${\tau }_{\text{OH}}=\frac{1}{k_{(\text{OH}+\text{MMA})}\times c_{\text{OH}}}$$

(3)

$${\tau }_{{NO}_3}=\frac{1}{k_{({\text{NO}}_3+\text{MMA})}\times c_{{\text{NO}}_3}}$$

(4)

The atmospheric lifetimes of MMA determined by OH and NO₃ radicals are 3 h and 6.5 days, respectively. The obtained lifetimes here can be comparable to 2–10 h for OH radicals [21] and 6 days for NO₃ radicals in the previous work [14]. Therefore, MMA is likely to be removed quickly by the reaction with OH radicals near to their emission sources. However, since the atmospheric lifetime of MMA were determined both by the rate constant of the reactions with radicals and the concentration of radicals, in some seriously polluted regions, the reaction with NO₃ radicals also may contribute to the removal of MMA in the atmosphere. Moreover, in some coastal areas where the concentration of Cl atoms can reach a peak value of 1 × 10⁵ molecule/cm³ [22,23], the MMA may be quikly removed by the reaction with Cl atoms [15].

3. Computational Methods

All the calculations were performed with the Gaussian 03 software package (Gaussian, Inc., Wallingford, CT, USA). Geometrical parameters of the reactants, intermediates, transition states, and products were fully optimized at the B3LYP level with a standard basis set 6-31G(d,p). The nature of stationary points, the zero-point energy (ZPE), and the thermal contribution to the free energy of activation were determined by the vibrational frequency calculation at the same level. For each transition state, the intrinsic reaction coordinate (IRC) calculation was performed using the same electronic structure theory to confirm it was connected to the right minima along the reaction pathway. The DFT geometries were then used in the single-point energy calculations at the frozen-core second-order Møller-Plesset perturbation theory (MP2) and the coupled-cluster theory with single and double excitations including perturbative corrections for the triple excitations (CCSD(T)) with several basis sets. Each single-point energy was further corrected using a factor, CF. This factor was determined from the energy difference between the MP2/6-311++G(d,p) and MP2/6-31G(d) levels, and the ranges of CF value are 0.001592 to 0.3571924 Hartree in the reactions of MMA with OH radical, and 0.010456 to 0.2470771 Hartree in the reactions of MMA with NO₃ radical. The values of calculated energies at the CCSD(T)/6-31G(d) level were then corrected by the CFs, corresponding to the CCSD(T)/6-31G(d) + CF level of theory [24,25]. Potential barriers (ΔE, ΔE = E_{transition state} − ∑E_reactants) and reaction heats (ΔH, ΔH = ∑E_products − ∑ _reactants) were determined in each channel.

Depending on different reaction types, two methods were used to calculate the rate constants. For the abstraction pathways, the widely used canonical variational transition-state theory with the small curvature tunneling (CVT/SCT) correction was adopted. The theoretical rate constants and their temperature dependence were calculated by using the Polyrate 9.3 program (University of Minnesota, Minneapolis, MN, USA). Actually, the CVT/SCT method has been successfully used in dealing with several bimolecular reactions [26]. For the addition pathways, the rate constants were calculated using multichannel Rice-Ramspergere-Kassele-Marcus (RRKM) theory. This method has been successfully used in the previous works for several addition reactions [27].

4. Conclusions

A theoretical study was performed on the mechanism of OH and NO₃ radical-initiated reactions of methyl methacrylate. The rate constants were calculated by using the CVT/SCT and multichannel RRKM method. Several specific conclusions can be drawn from this study:

(1)

The mechanism of OH and NO₃ radical-initiated reactions of MMA includes H abstraction pathways and the addition pathways, and the H abstraction pathways can be negligible because of their little contribution under the general atmospheric conditions.

(2)

The OH-MMA and NO₃-MMA adducts are open-shell activated radical intermediates, and can further react with O₂/NO in the atmosphere. For the OH radical-initiated reaction, the main products are CH₃C(O)C(O)OCH₃, HO₂ and methanal, consistent with the experimental results [15]. For the NO₃ radical-initiated oxidation, the reaction of the NO₃-MMA adduct with O₂/NO is thermodynamically and energetically more favorable than their unimolecular decomposition or isomerization under the general atmospheric conditions and the main products are CH₃C(O)C(O)OCH₃, CH₃C(O)CH₂NO₃ and methanal.

The calculated overall rate constants match well the available experimental values. The atmospheric life times of MMA determined by OH radicals and NO₃ radicals are 3 h and 6.5 days. NO₃-initiated oxidation reaction contributes little to the atmospheric losses of MMA except in polluted regions.