Quantitative kinetics reveal that reactions of HO2 are a significant sink for aldehydes in the atmosphere and may initiate the formation of highly oxygenated molecules via autoxidation

Large aldehydes are widespread in the atmosphere and their oxidation leads to secondary organic aerosols. The current understanding of their chemical transformation processes is limited to hydroxyl radical (OH) oxidation during daytime and nitrate radical (NO3) oxidation during nighttime. Here, we report quantitative kinetics calculations of the reactions of hexanal (C5H11CHO), pentanal (C4H9CHO), and butanal (C3H7CHO) with hydroperoxyl radical (HO2) at atmospheric temperatures and pressures. We find that neither tunneling nor multistructural torsion anharmonicity should be neglected in computing these rate constants; strong anharmonicity at the transition states is also important. We find rate constants for the three reactions in the range 3.2-7.7 × 10-14 cm3 molecule-1 s-1 at 298 K and 1 atm, showing that the HO2 reactions can be competitive with OH and NO3 oxidation under some conditions relevant to the atmosphere. Our findings reveal that HO2-initiated oxidation of large aldehydes may be responsible for the formation of highly oxygenated molecules via autoxidation. We augment the theoretic studies with laboratory flow-tube experiments using an iodide-adduct time-of-flight chemical ionization mass spectrometer to confirm the theoretical predictions of peroxy radicals and the autoxidation pathway. We find that the adduct from HO2 + C5H11CHO undergoes a fast unimolecular 1,7-hydrogen shift with a rate constant of 0.45 s-1. We suggest that the HO2 reactions make significant contributions to the sink of aldehydes.

Nature or number of species in a transition state: the key role of catalytically active molecules in hydrogen transfer stages in atmospheric aldehyde reactions

For the first time, the two factors (the number of sites in the transition state and the nature of the catalytically active species) that affect the energy barriers (Ea and ΔG‡) in atmospheric aldehyde reactions are proposed. The contribution of each factor to the energy barriers of the ammonization and amination stages, dehydration, and intramolecular hydrogen transfer is studied using the example of the acetaldehyde and glyoxal interactions with ammonia in aqueous solution. A regular decrease in energy barriers is observed in a series of 4-, 6-, and 8-membered transition states (TSs) regardless of the nature of the catalytically active species and their numbers. The 8-membered TSs of ammonization, amination, and dehydration reactions are the most efficient catalytic systems. The role of the nature of catalytically active species is secondary and is expressed in different cases through the influence of entropy and different acidity/basicity of catalytically active species and their structures. The regularities for the stage of intramolecular hydrogen transfer stand out from those for the ammonization, amination, and dehydration stages. The intramolecular hydrogen transfer is organized by three atoms in TSs without the participation of catalytically active species, while the 5- and 7-membered TSs are formed with the participation of such species. A proportional decrease in energy barrier with a sequential increase in the number of TS sites (3-, 5-, and 7-) is not observed. A sharp decrease in the barriers occurs only during the formation of the 7-membered TSs, while the 5-membered structures lie above the 3-membered catalytically inactive structures on the potential energy surface (PES) regardless of the nature of the species forming these structures.

New Reactions for the Formation of Organic Nitrate in the Atmosphere

Organic nitrates make an important contribution to the formation of secondary organic aerosols, but the formation mechanisms of organic nitrates are not fully understood at the molecular level. In the present work, we explore a new route for the formation of organic nitrates in the reaction of formaldehyde (HCHO) with nitric acid (HNO₃) catalyzed by water (H₂O), ammonia (NH₃), and dimethylamine ((CH₃)₂NH) using theoretical methods. The present results using CCSD(T)-F12a/cc-pVTZ-F12//M06-2X/MG3S unravel that dimethylamine has a stronger catalytic ability in the reaction of HCHO with HNO₃, reducing the barrier by 21.97 kcal/mol, while water and ammonia only decrease the energy barrier by 7.35 and 13.56 kcal/mol, respectively. In addition, the calculated kinetics combined with the corresponding concentrations of these species show that the HCHO + HNO₃ + (CH₃)₂NH reaction can compete well with the naked HCHO + HNO₃ reaction at 200-240 K, which may make certain contributions to the formation of organic nitrates under some atmospheric conditions.

Quantitative Kinetics of HO2 Reactions with Aldehydes in the Atmosphere: High-Order Dynamic Correlation, Anharmonicity, and Falloff Effects Are All Important

Kinetics provides the fundamental parameters for elucidating sources and sinks of key atmospheric species and for atmospheric modeling more generally. Obtaining quantitative kinetics in the laboratory for the full range of atmospheric temperatures and pressures is quite difficult. Here, we use computational chemistry to obtain quantitative rate constants for the reactions of HO₂ with HCHO, CH₃CHO, and CF₃CHO. First, we calculate the high-pressure-limit rate constants by using a dual-level strategy that combines conventional transition state theory using a high level of electronic structure wave function theory with canonical variational transition state theory including small-curvature tunneling using density functional theory. The wave-function level is beyond-CCSD(T) for HCHO and CCSD(T)-F12a (Level-A) for XCHO (X = CH₃, CF₃), and the density functional (Level-B) is specifically validated for these reactions. Then, we calculate the pressure-dependent rate constants by using system-specific quantum RRK theory (SS-QRRK) and also by an energy-grained master equation. The two treatments of the pressure dependence agree well. We find that the Level-A//Level-B method gives good agreement with CCSDTQ(P)/CBS. We also find that anharmonicity is an important factor that increases the rate constants of all three reactions. We find that the HO₂ + HCHO reaction has a significant dependence on pressure, but the HO₂ + CF₃CHO reaction is almost independent of pressure. Our findings show that the HO₂ + HCHO reaction makes important contribution to the sink for HCHO, and the HO₂ + CF₃CHO reaction is the dominant sink for CF₃CHO in the atmosphere.

High time-resolved variations of proteins in PM2.5 during haze pollution periods in Xi'an, China

Proteinaceous matter is an important component of PM_2.5, which can cause adverse health effects and also influence the air quality and climate change. However, there is little attention to high time-resolved variations and potential role of aerosol proteins during haze pollution periods. In this study, PM_2.5 samples were first collected by a medium flow sampler in autumn and winter in Xi'an, China. Then three high time-resolved monitoring campaigns during haze pollution periods were conducted to determine the evolving characteristics of total protein concentration and explore the interactive relationship between protein and other chemical compositions. The results showed that the average protein concentration in PM_2.5 in Xi'an (5.46 ± 3.32 μg m^-3) was higher than those in most cities of China, and varied by seasons and air pollution conditions. In particular, the protein concentration in PM_2.5 increased with the increase of air quality index (AQI). The continuous variations of aerosol proteins during the haze pollution periods further showed that PM_2.5, atmospheric humidity and long-distance air mass transport exerted the significant impacts on the protein components in aerosols. Based on the present observation, it is suggested that aerosol proteins might affect the generation of secondary aerosols under haze weather conditions. The present results may provide a new possible insight into the variations and the role of aerosol proteinaceous matter during the formation and development of haze pollution.

Temperature-dependent kinetics of the atmospheric reaction between CH2OO and acetone

Criegee intermediates are important oxidants produced in the ozonolysis of alkenes in the atmosphere. Quantitative kinetics of the reactions of Criegee intermediates are required for atmospheric modeling. However, the experimental studies do not cover the full relevant range of temperature and pressure. Here we report the quantitative kinetics of CH₂OO + CH₃C(O)CH₃ by using our recently developed dual strategy that combines coupled cluster theory with high excitation levels for conventional transition state theory and well validated levels of density functional theory for direct dynamics calculations using canonical variational transition theory including tunneling. We find that the W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ electronic structure method can be used to obtain quantitative kinetics of the CH₂OO + CH₃C(O)CH₃ reaction. Whereas previous investigations considered a one-step mechanistic pathway, we find that the CH₂OO + CH₃C(O)CH₃ reaction occurs in a stepwise manner. This has implications for the modeling of Criegee-intermediate reactions with other ketones and with aldehydes. In the kinetics calculations, we show that recrossing effects of the conventional transition state are negligible for determining the rate constant of CH₂OO + CH₃C(O)CH₃. The present findings reveal that the rate ratio between CH₂OO + CH₃C(O)CH₃ and OH + CH₃C(O)CH₃ has a significant negative dependence on temperature such that the CH₂OO + CH₃C(O)CH₃ reaction can contribute as a significant sink for atmospheric CH₃C(O)CH₃ at low temperature. The present findings should have broad implications in understanding the reactions of Criegee intermediates with carbonyl compounds and ketones in the atmosphere.

Large Pressure Effects Caused by Internal Rotation in the s-cis-syn-Acrolein Stabilized Criegee Intermediate at Tropospheric Temperature and Pressure

Criegee intermediates are important atmospheric oxidants, and quantitative kinetics for stabilized Criegee intermediates are key parameters for atmospheric modeling but are still limited. Here we report barriers and rate constants for unimolecular reactions of s-cis-syn-acrolein oxide (scsAO), in which the vinyl group makes it a prototype for Criegee intermediates produced in the ozonolysis of isoprene. We find that the MN15-L and M06-2X density functionals have CCSD(T)/CBS accuracy for the unimolecular cyclization and stereoisomerization of scsAO. We calculated high-pressure-limit rate constants by the dual-level strategy that combines (a) high-level wave function-based conventional transition-state theory (which includes coupled-cluster calculations with quasiperturbative inclusion of quadruple excitations because of the strongly multiconfigurational character of the electronic wave function) and (b) canonical variational transition-state theory with small-curvature tunneling based on a validated density functional. We calculated pressure-dependent rate constants both by system-specific quantum Rice-Ramsperger-Kassel theory and by solving the master equation. We report rate constants for unimolecular reactions of scsAO over the full range of atmospheric temperature and pressure. We found that the unimolecular reaction rates of this larger-than-previously studied Criegee intermediate depend significantly on pressure. Particularly, we found that falloff effects decrease the effective unimolecular cyclization rate constant of scsAO by about a factor of 3, but the unimolecular reaction is still the dominant atmospheric sink for scsAO at low altitudes. The large falloff caused by the inclusion of the stereoisomerization channel in the master equation calculations has broad implications for mechanistic analysis of reactions with competitive internal rotations that can produce stable rotamers.

Possible atmospheric source of NH2SO3H: the hydrolysis of HNSO2 in the presence of neutral, basic, and acidic catalysts

NH₂SO₃H can directly participate in H₂SO₄-(CH₃)₂NH-based cluster formation, and thereby substantially enhance the cluster formation rate. Herein, the reaction mechanisms and kinetics for the formation of NH₂SO₃H from the hydrolysis of HNSO₂ without and with neutral (H₂O, (H₂O)₂, and (H₂O)₃), basic (NH₃ and CH₃NH₂), and acidic (HCOOH, H₂SO₄, H₂SO₄⋯H₂O, and (H₂SO₄)₂) catalysts were studied theoretically at the CCSD(T)-F12/cc-pVDZ-F12//M06-2X/6-311+G(2df,2pd) level. The calculated results showed that neutral, basic, and acidic catalysts decrease the energy barrier by over 18.1 kcal mol^-1; meanwhile, the product formation of NH₂SO₃H was more strongly bonded to neutral, basic, and acidic catalysts than to the reactants HNSO₂ and H₂O. This reveals that the reported neutral, basic, and acidic catalysts promote the formation of NH₂SO₃H from the hydrolysis of HNSO₂ both kinetically and thermodynamically. Kinetic calculations using the master equation showed that (H₂O)₂ (100% RH) dominate over the other catalysts within the range of 0-10 km altitudes and 230-320 K with its rate ratio larger by at least 2.98 times, whereas HCOOH (3.2 × 10⁹ molecules cm^-3) is the most favorable catalysts at 15 km altitude in the troposphere. Overall, the present results will provide a definitive example that neutral, basic, and acidic catalysts have important influences on atmospheric reactions.

A computational study of the HO2 + SO3 → HOSO2 + 3O2 reaction catalyzed by water monomer, water dimer and small clusters of sulfuric acid: kinetics and atmospheric implications

Herein, the reaction mechanisms and kinetics for the HO2 + SO3 → HOSO2 + 3O2 reaction catalyzed by water monomer, water dimer and small clusters of sulfuric acid have been...

Atmospheric Ring-Closure and Dehydration Reactions of 1,4-Hydroxycarbonyls in the Gas Phase: The Impact of Catalysts

1,4-Hydroxycarbonyls can potentially undergo sequential reactions involving cyclization followed by dehydration to form dihydrofurans. As dihydrofurans contain a double bond, they are highly reactive toward atmospheric oxidants such as OH, O₃, and NO₃. In the present study, we use ab initio calculations to examine the impact of various atmospheric catalysts on the energetics and kinetics of the gas-phase cyclization and dehydration reaction steps associated with 4-hydroxybutanal, a prototypical 1,4-hydroxycarbonyl molecule. The cyclization step transforms 4-hydroxybutanal into 2-hydroxytetrahydrofuran, which can subsequently undergo dehydration to form 2,3-dihydrofuran. As the barriers associated with the cyclization and dehydration steps for 4-hydroxybutanal are, respectively, 34.8 and 63.0 kcal/mol in the absence of a catalyst, both reaction steps are inaccessible under atmospheric conditions in the gas phase. However, the presence of a suitable catalyst can significantly reduce the reaction barriers, and we have examined the impact of a single molecule of H₂O, HO₂ radical, HC(O)OH, HNO₃, and H₂SO₄ on these reactions. We find that H₂SO₄ reduces the reaction barriers the greatest, with the barrier for the cyclization step being reduced to -13.1 kcal/mol and that for the dehydration step going down to 9.2 kcal/mol, measured relative to their respective separated starting reactants. Interestingly, our kinetic study shows that HNO₃ gives the fastest rate due to the combined effects of a larger atmospheric concentration and a reduced barrier. Thus, our study suggests that, with acid catalysis, the cyclization reaction step can readily occur for 1,4-hydroxycarbonyls in the gas phase. Because the dehydration step exhibits a significant barrier even with acid catalysis, the 2-hydroxytetrahydrofuran products, once formed, are likely lost through their reaction with OH radicals in the atmosphere. We have investigated the reaction pathways and the rate constant for this bimolecular reaction in the presence of excess molecular oxygen (³O₂), as it would occur under tropospheric conditions, using computational chemistry over the 200-300 K temperature range. We find that the main products from these OH-initiated oxidation reactions are succinaldehyde + HO₂ and 2,3-dihydro-2-furanol + HO₂.

Atmospheric Kinetics: Bimolecular Reactions of Carbonyl Oxide by a Triple-Level Strategy

Criegee intermediates in the atmosphere serve as oxidizing agents to initiate aerosol formation, which are particularly important for atmospheric modeling, and understanding their kinetics is one of the current outstanding challenges in climate change modeling. Because experimental kinetics are still limited, we must rely on theory for the complete picture, but obtaining absolute rates from theory is a formidable task. Here, we report the bimolecular reaction kinetics of carbonyl oxide with ammonia, hydrogen sulfide, formaldehyde, and water dimer by designing a triple-level strategy that combines (i) benchmark results close to the complete-basis limit of coupled-cluster theory with the single, double, triple, and quadruple excitations (CCSDTQ/CBS), (ii) a new hybrid meta density functional (M06CR) specifically optimized for reactions of Criegee intermediates, and (iii) variational transition-state theory with both variable rection coordinates and optimized reaction paths, with multidimensional tunneling, and with pressure effects. For (i) we have found that quadruple excitations are required to obtain quantitative reaction barriers, and we designed new composite methods and strategies to reach CCSDTQ/CBS accuracy. The present findings show that (i) the CH₂OO + HCHO reaction can make an important contribution to the sink of HCHO under wide atmospheric conditions in the gas phase and that (ii) CH₂OO + (H₂O)₂ dominates over the CH₂OO + H₂O reaction below 10 km.

Atmospheric Chemistry of CH2OO: The Hydrolysis of CH2OO in Small Clusters of Sulfuric Acid

The hydrolysis of CH₂OO is not only a dominant sink for the CH₂OO intermediate in the atmosphere but also a key process in the formation of aerosols. Herein, the reaction mechanism and kinetics for the hydrolysis of CH₂OO catalyzed by the precursors of atmospheric aerosols, including H₂SO₄, H₂SO₄···H₂O, and (H₂SO₄)₂, have been studied theoretically at the CCSD(T)-F12a/cc-pVDZ-F12//B3LYP/6-311+G(2df,2pd) level. The calculated results show that the three catalysts decrease the energy barrier by over 10.3 kcal·mol^-1; at the same time, the product formation of HOCH₂OOH is more strongly bonded to the three catalysts than to the reactants CH₂OO and H₂O, revealing that small clusters of sulfuric acid promote the hydrolysis of CH₂OO both kinetically and thermodynamically. Kinetic simulations show that the H₂SO₄-assisted reaction is more favorable than the H₂SO₄···H₂O- (the pseudo-first-order rate constant being 27.9-11.5 times larger) and (H₂SO₄)₂- (between 2.8 × 10⁴ and 3.4 × 10⁵ times larger) catalyzed reactions. Additionally, due to relatively lower concentration of H₂SO₄, the hydrolysis of CH₂OO with H₂SO₄ cannot compete with the CH₂OO + H₂O or (H₂O)₂ reaction within the temperature range of 280-320 K, since its pseudo-first-order rate ratio is smaller by 4-7 or 6-8 orders of magnitude, respectively. However, the present results provide a good example of how small clusters of sulfuric acid catalyze the hydrolysis of an important atmospheric species.

New mechanistic pathways for the formation of organosulfates catalyzed by ammonia and carbinolamine formation catalyzed by sulfuric acid in the atmosphere

Sulfuric acid exerts a remarkable catalytic role in the H₂SO₄ + HCHO + NH₃ reaction that leads to the formation of carbinolamine.

Effect of NH3 and HCOOH on the H2O2 + HO → HO2 + H2O reaction in the troposphere: competition between the one-step and stepwise mechanisms

The H₂O₂ + HO → HO₂ + H₂O reaction is an important reservoir for both radicals of HO and HO₂ catalyzing the destruction of O₃. Here, this reaction assisted by NH₃ and HCOOH catalysts was explored using the CCSD(T)-F12a/cc-pVDZ-F12//M06-2X/aug-cc-pVTZ method and canonical variational transition state theory with small curvature tunneling. Two possible sets of mechanisms, (i) one-step routes and (ii) stepwise processes, are possible. Our results show that in the presence of both NH₃ and HCOOH catalysts under relevant atmospheric temperature, mechanism (i) is favored both energetically and kinetically than the corresponding mechanism (ii). At 298 K, the relative rate for mechanism (i) in the presence of NH₃ (10, 2900 ppbv) and HCOOH (10 ppbv) is respectively 3–5 and 2–4 orders of magnitude lower than that of the water-catalyzed reaction. This is due to a comparatively lower concentration of NH₃ and HCOOH than H₂O which indicates the positive water effect under atmospheric conditions. Although NH₃ and HCOOH catalysts play a negligible role in the reservoir for both radicals of HO and HO₂ catalyzing the destruction of O₃, the current study provides a comprehensive example of how acidic and basic catalysts assisted the gas-phase reactions.

The H₂O₂ + HO → HO₂ + H₂O reaction is an important reservoir for both radicals of HO and HO₂ catalyzing the destruction of O₃.

Unexpected quenching effect on new particle formation from the atmospheric reaction of methanol with SO3

Despite the high abundance in the atmosphere, alcohols in general and methanol in particular are believed to play a small role in atmospheric new particle formation (NPF) largely due to the weak binding abilities of alcohols with the major nucleation precursors, e.g., sulfuric acid (SA) and dimethylamine (DMA). Herein, we identify a catalytic reaction that was previously overlooked, namely, the reaction between methanol and SO₃, catalyzed by SA, DMA, or water. We found that alcohols can have unexpected quenching effects on the NPF process, particularly in dry and highly polluted regions with high concentrations of alcohols. Specifically, the catalytic reaction between methanol and SO₃ can convert methanol into a less-volatile species--methyl hydrogen sulfate (MHS). The latter was initially thought to be a good nucleation agent for NPF. However, our simulation results suggest that the formation of MHS consumes an appreciable amount of atmospheric SO₃, disfavoring further reactions of SO₃ with H₂O. Indeed, we found that MHS formation can cause a reduction of SA concentration up to 87%, whereas the nucleation ability of MHS toward new particles is not as good as that of SA. Hence, a high abundance of methanol in the atmosphere can lower the particle nucleation rate by as much as two orders of magnitude. Such a quenching effect suggests that the recently identified catalytic reactions between alcohols and SO₃ need to be considered in atmospheric modeling in order to predict SA concentration from SO₂, while also account for their potentially negative effect on NPF.

Hydrolysis of Formyl Fluoride Catalyzed by Sulfuric Acid and Formic Acid in the Atmosphere

Formyl fluoride (HFCO) is an important atmospheric molecule, and its reaction with the OH radical is an important pathway when degradation of HFCO is considered in earth's troposphere. Here, we study the hydrolysis of formyl fluoride (HFCO + H₂O) with sulfuric acid (H₂SO₄) and formic acid (HCOOH) acting as catalysts by utilizing M06-2X, CCSD(T)-F12a, and conventional transitional state theory with Eckart tunneling to explore the atmospheric impact of the above-said hydrolysis reactions. Our calculated results show that H₂SO₄ has a remarkably catalytic role in the gas-phase hydrolysis of HFCO, as the energy barriers of the HFCO + H₂O reaction are reduced from 39.22 and 41.19 to 0.26 and -0.63 kcal/mol with respect to the separate reactants, respectively. In addition, we also find that H₂SO₄ can significantly accelerate the decomposition of FCH(OH)₂ into hydrogen fluoride (HF) and HCOOH. This is because while the barrier height for the unimolecular decomposition of FCH(OH)₂ into HF and HCOOH is 31.63 kcal/mol, the barrier height for the FCH(OH)₂ + H₂SO₄ reaction is predicted to be -5.99 kcal/mol with respect to separate reactants. Nevertheless, the comparative relative rate analysis shows that the reaction between HFCO and the OH radical is still the most dominant pathway when the tropospheric degradation of HFCO is taken into account and that the gas-phase hydrolysis of HFCO may only occur with the help of H₂SO₄ when the atmospheric concentration of OH is about 10¹ molecules cm^-3 or less. Having an understanding from the present study that the gas-phase hydrolysis of HFCO in the presence of H₂SO₄ has very limited role possibly in the absence of sunlight, we also prefer here to emphasize that the HFCO + H₂O + H₂SO₄ reaction may occur on the surface of secondary organic aerosols for the formation of HCOOH.

Ammonolysis as an important loss process of acetaldehyde in the troposphere: energetics and kinetics of water and formic acid catalyzed reactions

The reaction of ammonia with acetaldehyde as a potential source of 1,1-aminoethanol in the troposphere has been investigated by electronic structure and chemical kinetics calculations.

Unimolecular reaction of acetone oxide and its reaction with water in the atmosphere

Criegee intermediates (i.e., carbonyl oxides with two radical sites) are known to be important atmospheric reagents; however, our knowledge of their reaction kinetics is still limited. Although experimental methods have been developed to directly measure the reaction rate constants of stabilized Criegee intermediates, the experimental results cover limited temperature ranges and do not completely agree well with one another. Here we investigate the unimolecular reaction of acetone oxide [(CH3)2COO] and its bimolecular reaction with H2O to obtain rate constants with quantitative accuracy comparable to experimental accuracy. We do this by using CCSDT(Q)/CBS//CCSD(T)-F12a/DZ-F12 benchmark results to select and validate exchange-correlation functionals, which are then used for direct dynamics calculations by variational transition state theory with small-curvature tunneling and torsional and high-frequency anharmonicity. We find that tunneling is very significant in the unimolecular reaction of (CH3)2COO and its bimolecular reaction with H2O. We show that the atmospheric lifetimes of (CH3)2COO depend on temperature and that the unimolecular reaction of (CH3)2COO is the dominant decay mode above 240 K, while the (CH3)2COO + SO2 reaction can compete with the corresponding unimolecular reaction below 240 K when the SO2 concentration is 9 × 1010 molecules per cubic centimeter. We also find that experimental results may not be sufficiently accurate for the unimolecular reaction of (CH3)2COO above 310 K. Not only does the present investigation provide insights into the decay of (CH3)2COO in the atmosphere, but it also provides an illustration of how to use theoretical methods to predict quantitative rate constants of medium-sized Criegee intermediates.