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

Theoretical investigation of the OH-initiated atmospheric degradation mechanism of CX2CHX (X = H, F, Cl) by advanced quantum chemical and transition state theory methods

Halogenated olefins are anthropogenic compounds with many industrial applications but at the same time raising many environmental and health concerns. Gas-phase electrophilic addition of the OH radical to the olefinic CC bond represents the primary sink for these chemicals in the atmosphere, with the degree and type of halogenation playing a significant role in their overall reactivity. In this work, we present a theoretical investigation of the reaction mechanisms and kinetics for the reactions between the OH radical and CH2CH2 (ethylene, ETH), CF2CHF (trifluoroethylene, TFE) and CCl2CHCl (trichloroethylene, TCE), simulated by state-of-the-art protocols and methods, with the aim of providing a detailed interpretation of the available experimental results, as well as new data of relevance to tropospheric chemistry. Specifically, potential energy surfaces (PESs) are obtained using the jun-Cheap (jChS) composite scheme, whereas temperature and pressure dependent rate coefficients and product distributions in the 100-600 K temperature range are calculated within the Rice-Ramsperger-Kassel-Marcus/master equation (RRKM/ME) framework. The rates for barrierless channels are obtained from variable reaction coordinate-variational transition state theory (VRC-VTST) combined with the two transition state model. While the reactions with ETH and TFE proceed mainly via the formation of addition adducts at P = 1 atm and T = 298 K, the dominant channel for TCE is the Cl-elimination reaction. Global rate constants for the two halogenated olefins, TFE and TCE, are found to be pressure-independent, contrary to the case of ETH. The computed rate constants, as well as their temperature and pressure dependence, are in remarkable agreement with the available experimental data, and they are used to derive atmospheric lifetimes (τ) for both TFE and TCE as a function of altitude (h) in the atmosphere, by taking into account variations in the rate coefficients (k (T, P)) and [OH] concentration.

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.

Continuing benefits of the Montreal Protocol and protection of the stratospheric ozone layer for human health and the environment

The protection of Earth's stratospheric ozone (O3) is an ongoing process under the auspices of the universally ratified Montreal Protocol and its Amendments and adjustments. A critical part of this process is the assessment of the environmental issues related to changes in O3. The United Nations Environment Programme's Environmental Effects Assessment Panel provides annual scientific evaluations of some of the key issues arising in the recent collective knowledge base. This current update includes a comprehensive assessment of the incidence rates of skin cancer, cataract and other skin and eye diseases observed worldwide; the effects of UV radiation on tropospheric oxidants, and air and water quality; trends in breakdown products of fluorinated chemicals and recent information of their toxicity; and recent technological innovations of building materials for greater resistance to UV radiation. These issues span a wide range of topics, including both harmful and beneficial effects of exposure to UV radiation, and complex interactions with climate change. While the Montreal Protocol has succeeded in preventing large reductions in stratospheric O3, future changes may occur due to a number of natural and anthropogenic factors. Thus, frequent assessments of potential environmental impacts are essential to ensure that policies remain based on the best available scientific knowledge.

The role of aldehydes on sulfur based-new particle formation: a theoretical study

Aldehydes play a crucial role in the formation of atmospheric particles, attracting significant attention due to their environmental impact. However, the microscopic mechanisms underlying the formation of aldehyde-involved particles remain uncertain. In this study, through quantum chemical calculations and molecular dynamics (MD) simulations, we investigate the microscopic formation mechanisms of binary and ternary systems composed of three representative aldehydes, two sulfur-based acids, water, and two bases. Our research findings reveal that the most stable structures of acid-aldehyde clusters involve the connection of acids and aldehyde compounds through hydrogen bonds without involving proton transfer reactions, indicating relatively poor cluster stability. However, with the introduction of a third component, the stability of 18 clusters significantly increase. Among these, in ten systems, acids act as catalysts, facilitating reactions between aldehyde compounds and water or alkaline substances to generate glycols and amino alcohols. However, according to MD simulations conducted at 300 K, these acids readily dissociate from the resulting products. In the remaining eight systems, the most stable structural feature involves ion pairs formed by proton transfer reactions between acids and aldehyde compounds. These clusters exhibit remarkable thermodynamic stability. Furthermore, the acidity of the acid, the nature of nucleophilic agents, and the type of aldehyde all play significant roles in cluster stability and reactivity, and they have synergistic effects on the nucleation process. This study offers microscopic insights into the processes of new particle formation involving aldehydes, contributing to a deeper understanding of atmospheric chemistry at the molecular level.

On the Chemical Pathways Influencing the Effective Global Warming Potential of Commercial Hydrofluoroolefin Gases

The enforcement of a global hydrofluorocarbon (HFC) refrigerant phase down led to the introduction of hydrofluoroolefins (HFOs) as a low Global Warming Potential (GWP) substitute, given their low atmospheric lifetime. However, to this date it is not fully clear the long-term atmospheric fate of HFOs primary degradation products: trifluoro acetaldehyde (TFE), trifluoro acetyl fluoride (TFF), and trifluoroacetic acid (TFA). It particularly concerns the possibility of forming HFC-23, a potent global warming agent. Although the atmospheric reaction networks of TFE, TFF, and TFA have a fair level of complexity, the relevant atmospheric chemical pathways are well characterized in the literature, enabling a comprehensive hazard assessment of HFC-23 formation as a secondary HFO breakdown product in diverse scenarios. A lower bound of the HFOs effective GWP in a baseline scenario is found above regulatory thresholds. While further research is crucial to refine climate risk assessments, the existing evidence suggests a non-negligible climate hazard associated with HFOs.

Quantitative kinetics of the atmospheric reaction between isocyanic acid and hydroxyl radicals: post-CCSD(T) contribution, anharmonicity, recrossing effects, torsional anharmonicity, and tunneling

Hydroxyl radicals (OH) are the most important atmospheric oxidant, initiating atmospheric reactions for the chemical transformation of volatile organic compounds. Here, we choose the HNCO + OH reaction as a prototype reaction because it contains the fundamental reaction processes for OH radicals: H-abstraction reaction by OH and OH addition reaction. However, its kinetics are unknown under atmospheric conditions. We investigate the reaction of HNCO with OH by using the GMM(P).L method close to the accuracy of single, double, triple, and quadruple excitations and noniterative quintuple excitations with a complete basis set (CCSDTQ(P)/CBS) as benchmark results and a dual-level strategy for kinetics calculations. The calculated rate constant of HNCO + OH is in good agreement with the experimental data available at the temperatures between 620 and 2500 K. We find that the rate constant cannot be correctly obtained by using experimental data to extrapolate the atmospheric temperature ranges. We find that the post-CCSD(T) contribution is very large for the barrier height with the value of -0.85 kcal mol-1 for the H-abstraction reaction, while the previous investigations were done up to the CCSD(T) level. Moreover, we also find that recrossing effects, tunneling, torsional anharmonicity, and anharmonicity are important for obtaining quantitative kinetics in the OH + HNCO reaction.

Hydroxy esters atmospheric degradation: OH and Cl reactivity, products distribution and fate of the alkoxy radicals formed

Kinetic studies of the reaction of ethyl glycolate HOCH2C(O)OCH2CH3 with OH radicals (kOH) and Cl atoms (kCl) have been conducted by the relative method using a glass atmospheric reactor by "in situ" Fourier Transform Infrared (FTIR) and Gas Chromatography equipped with flame ionization detection by Solid Phase Micro Extraction (GC-FID/SPME) at room temperature and atmospheric pressure. The following relative rate coefficients were determined using several reference compounds and two different techniques: kEG + OH-FTIR = (4.36 ± 1.21) × 10-12; kEG + OH-GC-FID= (3.90 ± 0.74) × 10-12; and kEG + Cl-GC-FID= (6.40 ± 0.72) × 10-11 all values in units of cm3.molecule-1.s-1. Complementary product studies were performed under comparable conditions to the kinetic tests, in order to identify the reaction products and to postulate their tropospheric oxidation mechanisms. The reaction of OH radicals and Cl atoms with ethyl glycolate initiates via H-atom abstraction from alkyl groups of the molecule. Formic acid was positively identified as a reaction product by FTIR. On the other hand, formaldehyde, acetaldehyde, glycolic acid; and formic acid were identified by the GC-MS technique. The Structure-Activity Relationship, (SAR) calculations were also implemented to estimate the more favorable reaction pathways and compare them with the products identified. Tropospheric lifetimes of τOH = 34 h and τCl = 5.5 days were estimated to determine how these investigated reactions might affect the air quality. In this sense, average ozone production of [O3] = 0.75 and a Photochemical Ozone Creation Potential, POCP, of 38 were calculated for the hydroxyl ester studied.

Kinetics of Sulfur Trioxide Reaction with Water Vapor to Form Atmospheric Sulfuric Acid

Although experimental methods can be used to obtain the quantitative kinetics of atmospheric reactions, experimental data are often limited to a narrow temperature range. The reaction of SO3 with water vapor is important for elucidating the formation of sulfuric acid in the atmosphere; however, the kinetics is uncertain at low temperatures. Here, we calculate rate constants for reactions of sulfur trioxide with two water molecules. We consider two mechanisms: the SO3···H2O + H2O reaction and the SO3 + (H2O)2 reaction. We find that beyond-CCSD(T) contributions to the barrier heights are very large, and multidimensional tunneling, unusually large anharmonicity of high-frequency modes, and torsional anharmonicity are important for obtaining quantitative kinetics. We find that at lower temperatures, the formation of the termolecular precursor complexes, which is often neglected, is rate-limiting compared to passage through the tight transition states. Our calculations show that the SO3···H2O + H2O mechanism is more important than the SO3 + (H2O)2 mechanism at 5-50 km altitudes. We find that the rate ratio between SO3···H2O + H2O and SO3 + (H2O)2 is greater than 20 at altitudes between 10 and 35 km, where the concentration of SO3 is very high.

Unimolecular Reactions of E-Glycolaldehyde Oxide and Its Reactions with One and Two Water Molecules

The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH₂OH)CHOO], for its reactions with H₂O and (H₂O)₂, and for the reaction of the E-(CH₂OH)CHOO…H₂O complex with H₂O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH₂OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH₂OH)CHOO…H₂O + H₂O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H₂O)₂. In addition, we find that the atmospheric lifetime of E-(CH₂OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10^-6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH₂OH)CHOO enhances its reactivity.

The atmospheric relevance of primary alcohols and imidogen reactions

Organic alcohols as very volatile compounds play a crucial role in the air quality of the atmosphere. So, the removal processes of such compounds are an important atmospheric challenge. The main goal of this research is to discover the atmospheric relevance of degradation paths of linear alcohols by imidogen with the aid of simulation by quantum mechanical (QM) methods. To this end, we combine broad mechanistic and kinetic results to get more accurate information and to have a deeper insight into the behavior of the designed reactions. Thus, the main and necessary reaction pathways are explored by well-behaved QM methods for complete elucidation of the studying gaseous reactions. Moreover, the potential energy surfaces as  a main factor are computed for easier judging of the most probable pathways in the simulated reactions. Our attempt to find the occurrence of the considered reactions in the atmospheric conditions is completed by precisely evaluating the rate constants of all elementary reactions. All of the computed bimolecular rate constants have a positive dependency on both temperature and pressure. The kinetic results show that H-abstraction from the α carbon is dominant relative to the other sites. Finally, by the results of this study, we conclude that at moderate temperatures and pressures primary alcohols can degrade with imidogen, so they can get atmospheric relevance.

Barrier heights, reaction energies and bond dissociation energies for RH + HO2 reactions with coupled-cluster theory, density functional theory and diffusion quantum Monte Carlo methods

Hydrogen abstraction reactions by the HO₂ radical from hydrocarbon molecules are an important class of reactions in the autoignition of hydrocarbon fuels. Performance of DLPNO-CC and DFT methods using three hybrids and four double hybrids as well as FN-DMC with the single-Slater-Jastrow trial wavefunction on barrier heights and reaction energies of RH + HO₂ reactions as well as bond dissociation energies of the involved X-H molecules is evaluated by comparison with the highly accurate CCSD(T)-F12b/CBS results in this study. Our results show that the DLPNO-CCSD(T)-F12 method can achieve highly accurate barrier heights, reaction energies and X-H bond energies for RH + HO₂ reactions at a relatively low computational cost, and it is applicable to the H-abstraction reactions of larger molecules. Among all DFAs, MN15 and the employed double hybrids can achieve accurate barrier heights and reaction energies with MADs of less than or around 2 kJ mol^-1, but their error on X-H bond energies is more pronounced. Only DSD-BLYP and DSD-PBEB95 can provide X-H bond energies with MADs less than 4 kJ mol^-1. Considering dispersion correction in DFT calculations does not improve these barrier heights and reaction energies. The error of FN-DMC on barrier heights and reaction energies is slightly larger than that of MN15 and those of double hybrids, but it can achieve results within chemical accuracy for these reactions and the X-H bond energies.

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.