Computational studies of atmospherically-relevant chemical reactions in water clusters and on liquid water and ice surfaces

Water-Catalyzed Formation of Reactive Oxygen Species from NO2 on a Weakly Hydrated Calcite Surface

The interfaces of weakly hydrated mineral substrates have been shown to serve as catalytic sites for chemical reactions that may not be accessible in the gas phase or under bulk conditions. Currently known mechanisms for the formation of reactive oxygen species (ROS) from nitrogen dioxide (NO2) involve NO2 dimerization. Here, we report the formation of the ROS HONO via a mechanism involving simple adsorption of a single NO2 molecule on a weakly hydrated calcite substrate. First-principles molecular dynamics simulations coupled with enhanced sampling techniques show how an adsorbed water sublayer can enhance NO2 adsorption on calcite compared to adsorption on a bare dry substrate. On the weakly hydrated calcite surface, an interfacial electric field facilitates proton extraction from water, thus allowing HONO formation from a single adsorbed NO2, i.e., without the need for the formation of a NO2 dimer precomplex. HONO formation on calcite is kinetically more favorable than that in the gas phase, with a reaction barrier of 14 kcal/mol on the weakly hydrated calcite surface compared to 27 kcal/mol in the gas phase. Further photocatalysed HONO production by visible light and HONO dissociation are hampered on calcite, unlike the process on silica. NO2 is a significant anthropogenic pollutant, and understanding its chemistry is crucial for explaining the high ROS levels and haze formation in polluted areas or prebiotic ROS generation. These findings emphasize how mineral substrates under water-restricted hydration conditions can trigger chemical pathways that are unexpected in the gas phase or under bulk conditions.

Molecular Insights into Chemical Reactions at Aqueous Aerosol Interfaces

Atmospheric aerosols facilitate reactions between ambient gases and dissolved species. Here, we review our efforts to interrogate the uptake of these gases and the mechanisms of their reactions both theoretically and experimentally. We highlight the fascinating behavior of N2O5 in solutions ranging from pure water to complex mixtures, chosen because its aerosol-mediated reactions significantly impact global ozone, hydroxyl, and methane concentrations. As a hydrophobic, weakly soluble, and highly reactive species, N2O5 is a sensitive probe of the chemical and physical properties of aerosol interfaces. We employ contemporary theory to disentangle the fate of N2O5 as it approaches pure and salty water, starting with adsorption and ending with hydrolysis to HNO3, chlorination to ClNO2, or evaporation. Flow reactor and gas-liquid scattering experiments probe even greater complexity as added ions, organic molecules, and surfactants alter the interfacial composition and reaction rates. Together, we reveal a new perspective on multiphase chemistry in the atmosphere.

Understanding the active role of water in laboratory chamber studies of reactions of the OH radical with alcohols of atmospheric relevance

In this work, we studied the reactions of three cyclic aliphatic alcohols with OH at room temperature, atmospheric pressure and different humidities in a Teflon reaction chamber. It was determined that the lower the solubility of the alcohol in water, the larger the effect of the humidity on the acceleration of the reaction. This experimental evidence allows suggesting that the acceleration is due to the reaction of the co-adsorbed reactants at the air-water interface of a thin water film deposited on the Teflon walls of the reaction chamber, instead of between co-reactants dissolved in the water film or due to gas phase catalysis as previously suggested. Therefore, formation of thin water films on different surfaces could have some implications on the tropospheric chemistry of these alcohols in the tropical regions of the planet with high humidity.

Formation of Chlorine in the Atmosphere by Reaction of Hypochlorous Acid with Seawater

The highly reactive dihalogens play a significant role in the oxidative chemistry of the troposphere. One of the main reservoirs of these halogens is hypohalous acids, HOX, which produce dihalogens in the presence of halides (Y-), where X, Y = Cl, Br, I. These reactions occur in and on aerosol particles and seawater surfaces and have been studied experimentally and by field observations. However, the mechanisms of these atmospheric reactions are still unknown. Here, we establish the atomistic mechanism of HOCl + Cl- → Cl2 + OH- at the surface of the water slab by performing ab initio molecular dynamics (AIMD) simulations. Main findings are (1) This reaction proceeds by halogen-bonded complexes of (HOCl)···(Cl-)aq surrounded with the neighboring water molecules. (2) The halogen bonded (HOCl)···(Cl-)aq complexes undergo charge transfer from Cl- to OH- to form transient Cl2 at neutral pH. (3) The addition of a proton to one proximal water greatly facilitates the Cl2 formation, which explains the enhanced rate at low pH.

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Significant influence of water molecules on the SO3 + HCl reaction in the gas phase and at the air-water interface

The products resulting from the reactions between atmospheric acids and SO3 have a catalytic effect on the formation of new particles in aerosols. However, the SO3 + HCl reaction in the gas-phase and at the air-water interface has not been considered. Herein, this reaction was explored exhaustively by using high-level quantum chemical calculations and Born Oppenheimer molecular dynamics (BOMD) simulations. The quantum calculations show that the gas-phase reaction of SO3 + HCl is highly unlikely to occur under atmospheric conditions with a high energy barrier of 22.6 kcal mol-1. H2O and (H2O)2 play obvious catalytic roles in reducing the energy barrier of the SO3 + HCl reaction by over 18.2 kcal mol-1. The atmospheric lifetimes of SO3 show that the (H2O)2-assisted reaction dominates over the H2O-assisted reaction within the altitude range of 0-5 km, whereas the H2O-assisted reaction is more favorable within an altitude range of 10-50 km. BOMD simulations show that H2O-induced formation of the ClSO3-⋯H3O+ ion pair and HCl-assisted formation of the HSO4-⋯H3O+ ion pair were identified at the air-water interface. These routes followed a stepwise reaction mechanism and proceeded at a picosecond time scale. Interestingly, the formed ClSO3H in the gas phase has a tendency to aggregate with sulfuric acids, ammonias, and water molecules to form stable clusters within 40 ns simulation time, while the interfacial ClSO3- and H3O+ can attract H2SO4, NH3, and HNO3 for particle formation from the gas phase to the water surface. Thus, this work will not only help in understanding the SO3 + HCl reaction driven by water molecules in the gas-phase and at the air-water interface, but it will also provide some potential routes of aerosol formation from the reaction between SO3 and inorganic acids.

Atmospheric Chemistry of NH2SO3H in Polluted Areas: An Unexpected Isomerization of NH2SO3H in Acid-Polluted Regions

NH2SO3H is an effective nucleation agent for the formation of atmospheric aerosols and cloud particles. So, the ammonolysis of SO3 to form NH2SO3H without and with neutral (H2O) and basic (NH3) trace gases has been extensively investigated. However, the acidic trace gas X (X = H2SO4 and CH3SO3H)-assisted ammonolysis of SO3 is still up for debate. In this work, a comprehensive theoretical investigation of X-assisted ammonolysis of SO3 and its reverse reaction (the isomerization of NH2SO3H to form SO3-···NH3+) was carried out in the gas phase and at the air-water interface. The gas-phase results show that X-assisted isomerization of NH2SO3H to form SO3-···NH3+ is more energetically and kinetically favorable than its reverse reaction and the isomerization of NH2SO3H in the presence of H2O and NH3. Such unexpected findings revealed that gas-phase NH2SO3H is highly reactive in the presence of acidic trace gas in contrast to the high stability of NH2SO3H in neutral and basic conditions. At the air-water interface, the X-assisted isomerization reaction of NH2SO3H involves multiple water molecules. The loop structure of the reaction center (X···NH2SO3H···3H2O) promotes the transfer of protons in the water molecules to form the SO3-···NH3+ ion pair, which can then interact with several interfacial water molecules to form ammonium bisulfate. These interfacial reaction channels follow a stepwise mechanism and proceed at the picosecond time-scale. The findings of this study will contribute to a better understanding of the atmospheric behavior of NH2SO3H in polluted acidic trace gases.

Distinctive Heterogeneous Reaction Mechanism of ClNO2 on the Air-Water Surface Containing Cl

The heterogeneous reaction of nitryl chloride (ClNO2) on the air-water surface plays a significant role in the chloride lifecycle. The air-water surface is ubiquitous on ice surfaces under supercooled conditions, affecting the uptake and heterogeneous reaction processes of trace gases. Previous studies suggest that ClNO2 is formed on Cl-doped ice surfaces following the N2O5 uptake. Herein, a distinctive heterogeneous reaction mechanism of ClNO2 is suggested on an air-water surface containing Cl under supercooled conditions using combined classic molecular dynamics (MD) and Born-Oppenheimer MD simulations. It is found that N2O5 dissociates into a NO2+ and NO3- ionic pair on the top air-water surface. In the top layer of the surface containing barely any Cl-, NO2+ proceeds through hydrolysis and produces H3O+ and HNO3. Thus, surface acidification appears because of H3O+ yields. With NO2+ diffusion to the deep layer of the surface, NO2+ reacts with Cl- and forms ClNO2. Note that ClNO2 formation competes with NO2+ hydrolysis, and the rate of ClNO2 formation is 27.7[Cl-] larger than that of NO2+ hydrolysis. Afterward, the reaction of ClNO2 with Cl- becomes barrierless with the catalysis by H3O+, which is not feasible on a neutral surface. Cl2 is thus generated and escapes into the atmosphere (low solubility of Cl2), contributing to the Cl radical. The proposed mechanism bolsters the current understanding of ClNO2's fate and its role in Cl chemistry in extremely cold environments like the Arctic and other high-latitude regions in wintertime.

Multiple evaluations of atmospheric behavior between Criegee intermediates and HCHO: Gas-phase and air-water interface reaction

Given the high abundance of water in the atmosphere, the reaction of Criegee intermediates (CIs) with (H₂O)₂ is considered to be the predominant removal pathway for CIs. However, recent experimental findings reported that the reactions of CIs with organic acids and carbonyls are faster than expected. At the same time, the interface behavior between CIs and carbonyls has not been reported so far. Here, the gas-phase and air-water interface behavior between Criegee intermediates and HCHO were explored by adopting high-level quantum chemical calculations and Born-Oppenheimer molecular dynamics (BOMD) simulations. Quantum chemical calculations evidence that the gas-phase reactions of CIs + HCHO are submerged energy or low energy barriers processes. The rate ratios speculate that the HCHO could be not only a significant tropospheric scavenger of CIs, but also an inhibitor in the oxidizing ability of CIs on SO_x in dry and highly polluted areas with abundant HCHO concentration. The reactions of CH₂OO with HCHO at the droplet's surface follow a loop structure mechanism to produce i) SOZ (), ii) BHMP (HOCH₂OOCH₂OH), and iii) HMHP (HOCH₂OOH). Considering the harsh reaction conditions between CIs and HCHO at the interface (i.e., the two molecules must be sufficiently close to each other), the hydration of CIs is still their main atmospheric loss pathway. These results could help us get a better interpretation of the underlying CIs-aldehydes chemical processes in the global polluted urban atmospheres.

Effect of Linker Structure and Functionalization on Secondary Gas Formation in Metal-Organic Frameworks

Rare-earth terephthalic acid (BDC)-based metal-organic frameworks (MOFs) are promising candidate materials for acid gas separation and adsorption from flue gas streams. However, previous simulations have shown that acid gases (H₂O, NO₂, and SO₂) react with the hydroxyl on the BDC linkers to form protonated acid gases as a potential degradation mechanism. Herein, gas-phase computational approaches were used to identify the formation energies of these secondary protonated acid gases across multiple BDC linker molecules. Formation energies for secondary protonated acid gases were evaluated using both density functional theory (DFT) and correlated wave function methods for varying BDC-gas reaction mechanisms. Upon validation of DFT to reproduce wave function calculation results, rotated conformational linkers and chemically functionalized BDC linkers with -OH, -NH₂, and -SH were investigated. The calculations show that the rotational conformation affects the molecule stability. Double-functionalized BDC linkers, where two functional groups are substituted onto BDC, showed varied reaction energies depending on whether the functional groups donate or withdraw electrons from the aromatic system. Based on these results, BDC linker design must balance adsorption performance with degradation via linker dehydrogenation for the design of stable MOFs for acid gas separations.

Theoretical Study on the Gas-Phase and Aqueous Interface Reaction Mechanism of Criegee Intermediates with 2-Methylglyceric Acid and the Nucleation of Products

Criegee intermediates (CIs) are important in the sink of many atmospheric substances, including alcohols, organic acids, amines, etc. In this work, the density functional theory (DFT) method was used to calculate the energy barriers for the reactions of CH3CHOO with 2-methyl glyceric acid (MGA) and to evaluate the interaction of the three functional groups of MGA. The results show that the reactions involving the COOH group of MGA are negligibly affected, and that hydrogen bonding can affect the reactions involving α-OH and β-OH groups. The water molecule has a negative effect on the reactions of the COOH group. It decreases the energy barriers of reactions involving the α-OH and β-OH groups as a catalyst. The Born-Oppenheimer molecular dynamic (BOMD) was applied to simulate the reactions of CH3CHOO with MGA at the gas-liquid interface. Water molecule plays the role of proton transfer in the reaction. Gas-phase calculations and gas-liquid interface simulations demonstrate that the reaction of CH3CHOO with the COOH group is the main pathway in the atmosphere. The molecular dynamic (MD) simulations suggest that the reaction products can form clusters in the atmosphere to participate in the formation of particles.

Effects of Microhydration on the Mechanisms of Hydrolysis and Cl- Substitution in Reactions of N2 O5 and Seawater

The reaction of N₂ O₅ at atmospheric interfaces has recently received considerable attention due to its importance in atmospheric chemistry. N₂ O₅ reacts preferentially with Cl^- to form ClNO₂ /NO₃ ^- (Cl^- substitution), but can also react with H₂ O to form 2HNO₃ (hydrolysis). In this paper, we explore these competing reactions in a theoretical study of the clusters N₂ O₅ /Cl^- /nH₂ O (n=2-5), resulting in the identification of three reaction motifs. First, we uncovered an S_N 2-type Cl^- substitution reaction of N₂ O₅ that occurs very quickly due to low barriers to reaction. Second, we found a low-lying pathway to hydrolysis via a ClNO₂ intermediate (two-step hydrolysis). Finally, we found a direct hydrolysis pathway where H₂ O attacks N₂ O₅ (one-step hydrolysis). We find that Cl^- substitution is the fastest reaction in every cluster. Between one-step and two-step hydrolysis, we find that one-step hydrolysis barriers are lower, making two-step hydrolysis (via ClNO₂ intermediate) likely only when concentrations of Cl^- are high.

Computational chemistry of cluster: Understanding the mechanism of atmospheric new particle formation at the molecular level

New particle formation (NPF), which exerts significant influence over human health and global climate, has been a hot topic and rapidly expands field of research in the environmental and atmospheric chemistry recent years. Generally, NPF contains two processes: formation of critical nucleus and further growth of the nucleus. However, due to the complexity of the atmospheric nucleation, which is a multicomponent process, formation of critical clusters as well as their growth is still connected to large uncertainties. Detection limits of instruments in measuring specific gaseous aerosol precursors and chemical compositions at the molecular level call for computational studies. Computational chemistry could effectively compensate the deficiency of laboratory experiments as well as observations and predict the nucleation mechanisms. We review the present theoretical literatures that discuss nucleation mechanism of atmospheric clusters. Focus of this review is on different nucleation systems involving sulfur-containing species, nitrogen-containing species and iodine-containing species. We hope this review will provide a deep insight for the molecular interaction of nucleation precursors and reveal nucleation mechanism at the molecular level.

Local Ice-like Structure at the Liquid Water Surface

Experiments and computer simulations have established that liquid water's surfaces can deviate in important ways from familiar bulk behavior. Even in the simplest case of an air-water interface, distinctive layering, orientational biases, and hydrogen bond arrangements have been reported, but an overarching picture of their origins and relationships has been incomplete. Here we show that a broad set of such observations can be understood through an analogy with the basal face of crystalline ice. Using simulations, we demonstrate a number of structural similarities between water and ice surfaces, suggesting the presence of domains at the air-water interface with ice-like features that persist over 2-3 molecular diameters. Most prominent is a shared characteristic layering of molecular density and orientation perpendicular to the interface. Lateral correlations of hydrogen bond network geometry point to structural similarities in the parallel direction as well. Our results bolster and significantly extend previous conceptions of ice-like structure at the liquid's boundary and suggest that the much-discussed quasi-liquid layer on ice evolves subtly above the melting point into a quasi-ice layer at the surface of liquid water.

Preparation and Characterization of Metastable trans-Dinitrogen Tetroxide

Under atmospheric conditions, NO₂ is in equilibrium with its dimers, N₂O₄, which can exist in the form of constitutional isomers and stereoisomers whose relative stabilities and reactivities are still being debated. Experimental limitations facing the spectroscopic characterization of the isomers of N₂O₄ prevent us from determining their relative contributions to reaction mechanisms possibly causing discrepancies in the reported reaction orders and rates. Using reflection-absorption infrared spectroscopy, molecular beam deposition, and matrix isolation techniques, it is shown that the relative abundances of NO₂ and its dimers can be controlled by heating or cooling the deposited gas. The comparison of spectra acquired from samples prepared using molecular beam deposition with those obtained using tube dosing deposition demonstrates how the N₂O₄ isomer distributions are sensitive to details of the experimental conditions and sample preparation protocols. These observations not only provide a better understanding of a possible source for the disagreements found in the literature, but also a methodology to control and quantify the chemical speciation in NO₂ vapors in terms of the relative abundances of NO₂ and of the various isomers of N₂O₄.

Chemistrees: Data-Driven Identification of Reaction Pathways via Machine Learning

We propose to analyze molecular dynamics (MD) output via a supervised machine learning (ML) algorithm, the decision tree. The approach aims to identify the predominant geometric features which correlate with trajectories that transition between two arbitrarily defined states. The data-driven algorithm aims to identify these features without the bias of human “chemical intuition”. We demonstrate the method by analyzing the proton exchange reactions in formic acid solvated in small water clusters. The simulations were performed with ab initio MD combined with a method to efficiently sample the rare event, path sampling. Our ML analysis identified relevant geometric variables involved in the proton transfer reaction and how they may change as the number of solvating water molecules changes.

Mechanistic Study of the Aqueous Reaction of Organic Peroxides with HSO3 - on the Surface of a Water Droplet

Aqueous reactions between organic peroxides and SO₂ are of intense interest in atmospheric science because of their ubiquitous implications for sulfate formation in secondary aerosols. However, the relative yields of the reaction products (inorganic vs. organic sulfates) remain controversial (i.e., 90 % vs. 40-70 % for inorganic sulfate) due in part to the lack of understanding of the underlying reaction mechanisms. Here, our computational results suggest that the reactions of HSO₃ ^- (dissolved SO₂ ) with organic peroxides are initiated on the surface of water nanodroplets and then proceed under two reaction pathways, in which the S atom of HSO₃ ^- attacks either the O1 or O2 atom of the peroxide group -O(O2)O(O1)H, leading to the formation of inorganic and organic sulfates, respectively. Notably, we find that thse reaction initiated by O1 atom exhibits a relatively low energy barrier and high reaction rate, which favours the formation of inorganic sulfate.

Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface

Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.

Weighted-Graph-Theoretic Methods for Many-Body Corrections within ONIOM: Smooth AIMD and the Role of High-Order Many-Body Terms

We present a weighted-graph-theoretic approach to adaptively compute contributions from many-body approximations for smooth and accurate post-Hartree-Fock (pHF) ab initio molecular dynamics (AIMD) of highly fluxional chemical systems. This approach is ONIOM-like, where the full system is treated at a computationally feasible quality of treatment (density functional theory (DFT) for the size of systems considered in this publication), which is then improved through a perturbative correction that captures local many-body interactions up to a certain order within a higher level of theory (post-Hartree-Fock in this publication) described through graph-theoretic techniques. Due to the fluxional and dynamical nature of the systems studied here, these graphical representations evolve during dynamics. As a result, energetic "hops" appear as the graphical representation deforms with the evolution of the chemical and physical properties of the system. In this paper, we introduce dynamically weighted, linear combinations of graphs, where the transition between graphical representations is smoothly achieved by considering a range of neighboring graphical representations at a given instant during dynamics. We compare these trajectories with those obtained from a set of trajectories where the range of local many-body interactions considered is increased, sometimes to the maximum available limit, which yields conservative trajectories as the order of interactions is increased. The weighted-graph approach presents improved dynamics trajectories while only using lower-order many-body interaction terms. The methods are compared by computing dynamical properties through time-correlation functions and structural distribution functions. In all cases, the weighted-graph approach provides accurate results at a lower cost.

Barrierless HONO and HOS(O)2-NO2 Formation via NH3-Promoted Oxidation of SO2 by NO2

In the troposphere, the knowledge about nitrous acid (HONO) sources is incomplete. The missing source of sulfate and fine particles cannot be explained during haze events. Air quality models cannot predict high levels of secondary fine-particle pollution. Despite extensive studies, one challenging issue in atmospheric chemistry is identifying the source of HONO. Here, we present direct ab initio molecular dynamics simulation evidence and typical air pollution events of the formation of gaseous HONO, nitrogen dioxide/hydrogen sulfite (HOS(O)2-NO₂ or NO₂-HSO₃) from nitrogen dioxide (NO₂), sulfur dioxide (SO₂), water (H₂O), and ammonia (NH₃) molecules in a proportion of 2:1:3:3. The reactions show a new mechanism for the formation of HONO and NO₂-HSO₃ in the troposphere, especially when the concentration of NO₂, SO₂, H₂O, and NH₃ is high (e.g., 2:1:3:3 or higher) in the air. Contrary to the proportion NO₂, SO₂, H₂O, and NH₃ equaling to 1:1:3:1 and 1:1:3:2, the proportion (2:1:3:3) enables barrierless reactions and weak interactions between molecules via the formation of HONO, NO₂-HSO₃, and NH₃/H₂O. In addition, field observations are carried out, and the measured data are summarized. Correlation analysis supported the conversion of NO₂ to HONO during observational studies. The weak interactions promote proton transfer, resulting in the generation of HONO, NO₂-HSO₃, and NH₃/H₂O pairs.

Molecular Mechanisms and Atmospheric Implications of Criegee Intermediate-Alcohol Chemistry in the Gas Phase and Aqueous Surface Environments

Criegee intermediates and alcohols are important species in the atmosphere. In this study, we use quantum chemistry and Born-Oppenheimer molecular dynamics (BOMD) simulations to investigate the reaction between methanol/ethanol and Criegee intermediates (anti- or syn-CH₃CHOO) in the gas phase and at the air-water interface. Reactions at the interface are found to be much faster than those in the gas phase. When water molecules are available, loop structures can be formed to facilitate the reaction. In addition, nonloop reaction pathways characterized by the formation of hydrated protons, although with a low possibility, are also identified at the air-water interface. Implications of our results on the fate of Criegee intermediates in the atmosphere are discussed, which deepen our understanding of Criegee intermediate-alcohol chemistry in humid environments.

Determination of the amine-catalyzed SO3 hydrolysis mechanism in the gas phase and at the air-water interface

New particle formation (NPF) involving amines in the atmosphere is considered an aggregation process, during which stable molecular clusters are formed from amines and sulfuric acid via hydrogen bond interaction. In this work, ab initio dynamics simulations of ammonium bisulfate formation from a series of amines, SO₃, and H₂O molecules were carried out in the gas phase and at the air-water interface. The results show that reactions between amines and hydrated SO₃ molecules in the gas phase are barrierless or nearly barrierless processes. The reaction rate is related to the basicity of gas-phase amines-the stronger the basicity, the faster the reaction. Furthermore, SO₃ hydrolysis catalyzed by amines occurs simultaneously with H₂SO₄-amine cluster formation. At the air-water interface, reactions between amines and SO₃ involve multiple water molecules. The reaction center's ring structure (amine-SO₃-nH₂O) promotes the transfer of protons in the water molecules. The formed ammonium cation (-RNH₃⁺) and the bisulfate anion (HSO₄^-) are present and stable by means of hydrogen bond interaction. The cluster formation mechanism provides new insights into NPF involving amines, which may play an important role in the formation of aerosols in some heavily polluted areas - e.g., those with a high amine concentration.

Photodecay of guaiacol is faster in ice, and even more rapid on ice, than in aqueous solution

Snowpacks contain a wide variety of inorganic and organic compounds, including some that absorb sunlight and undergo direct photoreactions. How the rates of these reactions in, and on, ice compare to rates in water is unclear: some studies report similar rates, while others find faster rates in/on ice. Further complicating our understanding, there is conflicting evidence whether chemicals react more quickly at the air-ice interface compared to in liquid-like regions (LLRs) within the ice. To address these questions, we measured the photodegradation rate of guaiacol (2-methoxyphenol) in various sample types, including in solution, in ice, and at the air-ice interface of nature-identical snow. Compared to aqueous solution, we find modest rate constant enhancements (increases of 3- to 6-fold) in ice LLRs, and much larger enhancements (of 17- to 77-fold) at the air-ice interface of nature-identical snow. Our computational modeling suggests the absorption spectrum for guaiacol red-shifts and increases on ice surfaces, leading to more light absorption, but these changes explain only a small portion (roughly 2 to 9%) of the observed rate constant enhancements in/on ice. This indicates that increases in the quantum yield are primarily responsible for the increased photoreactivity of guaiacol on ice; relative to solution, our results suggest that the quantum yield is larger by a factor of roughly 3-6 in liquid-like regions and 12-40 at the air-ice interface.

Molecular reactions at aqueous interfaces

This Review aims to critically analyse the emerging field of chemical reactivity at aqueous interfaces. The subject has evolved rapidly since the discovery of the so-called 'on-water catalysis', alluding to the dramatic acceleration of reactions at the surface of water or at its interface with hydrophobic media. We review critical experimental studies in the fields of atmospheric and synthetic organic chemistry, as well as related research exploring the origins of life, to showcase the importance of this phenomenon. The physico-chemical aspects of these processes, such as the structure, dynamics and thermodynamics of adsorption and solvation processes at aqueous interfaces, are also discussed. We also present the basic theories intended to explain interface catalysis, followed by the results of advanced ab initio molecular-dynamics simulations. Although some topics addressed here have already been the focus of previous reviews, we aim at highlighting their interconnection across diverse disciplines, providing a common perspective that would help us to identify the most fundamental issues still incompletely understood in this fast-moving field.

On the Theoretical Determination of Photolysis Properties for Atmospheric Volatile Organic Compounds

Volatile organic compounds (VOCs) are ubiquitous atmospheric molecules that generate a complex network of chemical reactions in the troposphere, often triggered by the absorption of sunlight. Understanding the VOC composition of the atmosphere relies on our ability to characterize all of their possible reaction pathways. When considering reactions of (transient) VOCs with sunlight, the availability of photolysis rate constants, utilized in general atmospheric models, is often out of experimental reach due to the unstable nature of these molecules. Here, we show how recent advances in computational photochemistry allow us to calculate in silico the different ingredients of a photolysis rate constant, namely, the photoabsorption cross-section and wavelength-dependent quantum yields. The rich photochemistry of tert-butyl hydroperoxide, for which experimental data are available, is employed to test our protocol and highlight the strengths and weaknesses of different levels of electronic structure and nonadiabatic molecular dynamics to study the photochemistry of (transient) VOCs.

Determination of reactions between Criegee intermediates and methanesulfonic acid at the air-water interface

In recent years, Criegee chemistry has become an important research focus due to its relevance in regulating concentrations of tropospheric OH radicals, hydroperoxides, sulfates, nitrates, and aerosols. However, to date, its interface behavior remains poorly understood. Thus, in this study, we used the Born-Oppenheimer molecular dynamics (BOMD) simulation method to explore the reaction mechanisms between Criegee intermediates (CIs) and methylsulfonic acid (MSA) at the air-water interface, then compared the observed behaviors with those in the gas phase. The addition of Criegee intermediates to MSA is nearly a barrierless reaction and follows a loop-structure mechanism in the gas phase. The high rate constants indicate that the Criegee intermediates and MSA reactions are the main acid removal channels. At the water's surface, the interaction of Criegee intermediates with MSA includes three main channels: 1) direct addition reaction, 2) H₂O-mediated hydroperoxide formation, and 3) MSA-mediated Criegee hydration. These reaction channels follow a loop-structure or a stepwise mechanism and proceed at the picosecond time-scale. The results of this work broaden our understanding of Criegee atmospheric behaviors in polluted urban and marine areas, which in turn will aid in developing more effective pollution control measures.

Tailoring electric field standing waves in reflection-absorption infrared spectroscopy to enhance absorbance from adsorbates on ice surfaces

The spectroscopic detection of molecules adsorbed onto ice surfaces at coverages similar to those encountered under typical environmental conditions requires high surface selectivity and sensitivity that few techniques can afford. An experimental methodology allowing a significant enhancement in the absorbance from adsorbed molecules is demonstrated herein. It exploits Electric Field Standing Wave (EFSW) effects intrinsic to grazing incidence Reflection-Absorption Infrared (RAIR) spectroscopy, where film thickness dependent optical interferences occur between the multiple reflections of the IR beam at the film-vacuum and the substrate-film interfaces. In this case study, CH₄ is used as a probe molecule and is deposited on a 20 ML coverage dense amorphous solid water film adsorbed onto solid Ar underlayers of various thicknesses. We observe that, at thicknesses where destructive interferences coincide with the absorption features from the CH stretching and HCH bending vibrational modes of methane, their intensity increases by a factor ranging from 10 to 25. Simulations of the RAIR spectra of the composite stratified films using a classical optics model reproduce the Ar underlayer coverage dependent enhancements of the absorbance features from CH₄ adsorbed onto the ice surface. They also reveal that the enhancements occur when the square modulus of the total electric field at the film's surface reaches its minimum value. Exploiting the EFSW effect allows the limit of detection to be reduced to a coverage of (0.2 ± 0.2) ML CH₄, which opens up interesting perspectives for spectroscopic studies of heterogeneous atmospheric chemistry at coverages that are more representative of those found in the natural environment.

Microscopic Mechanisms of N2O5 Hydrolysis on the Surface of Water Droplets

Reactions of N₂O₅, in particular heterogeneous hydrolysis, play a vital role in determining the chemistry of the atmosphere. The N₂O₅ heterogeneous hydrolysis reaction has been the subject of extensive research for decades, yet the physicochemical details of the mechanism have not been established. In this study, we show that this reaction can occur on the surface of a pure water droplet. We compute a relevant transition state for a nano-size model system and follow its evolution in time by means of ab initio molecular dynamics. This transition state, where N₂O₅ has a strong ion-pair character, leads directly to HNO₃. Both electrophilic and nucleophilic mechanisms take place. It is suggested that corresponding simulations for hydrolysis in the bulk are desirable.

Thermally Induced Hydrogen-Bond Rearrangements in Small Water Clusters and the Persistent Water Tetramer

Small water clusters absorb heat and catalyze pivotal atmospheric reactions. Yet, experiments produced conflicting results on water cluster distribution under atmospheric conditions. Additionally, it is unclear which "phase transitions" such clusters exhibit, at what temperatures, and what are their underlying molecular mechanisms. We find that logarithmically small tails in the radial probability densities of (H₂O) _n clusters (n = 2 - 6) provide direct testimony for such transitions. Using the best available water potential (MB-pol), an advanced thermostating algorithm (g-BAOAB), and sufficiently long trajectories, we map the "bifurcation", "melting", and (hitherto unexplored) "vaporization" transitions, finding that both melting and vaporization proceed via a "monomer on a ring" conformer, exhibiting huge distance fluctuations at the vaporization temperatures (T _v). T _v may play a role in determining the atmospheric cluster size distribution such that the dimer and tetramer, with their exceptionally low/high T _v values, are under/over-represented in these distributions, as indeed observed in nondestructive mass spectrometric measurements.

Modelling consortium for chemistry of indoor environments (MOCCIE): integrating chemical processes from molecular to room scales

We report on the development of a modelling consortium for chemistry in indoor environments that connects models over a range of spatial and temporal scales, from molecular to room scales and from sub-nanosecond to days, respectively. Our modeling approaches include molecular dynamics (MD) simulations, kinetic process modeling, gas-phase chemistry modeling, organic aerosol modeling, and computational fluid dynamics (CFD) simulations. These models are applied to investigate ozone reactions with skin and clothing, oxidation of volatile organic compounds and formation of secondary organic aerosols, and mass transport and partitioning of indoor species to surfaces. MD simulations provide molecular pictures of limonene adsorption on SiO2 and ozone interactions with the skin lipid squalene, providing kinetic parameters such as surface accommodation coefficient, desorption lifetime, and bulk diffusivity. These parameters then constrain kinetic process models, which resolve mass transport and chemical reactions in gas and condensed phases for analysis of experimental data. A detailed indoor chemical box model is applied to simulate α-pinene ozonolysis with improved representation of gas-particle partitioning. Application of 2D-volatility basis set reveals that OH-induced aging sometimes drives increases in indoor organic aerosol concentrations, due to organic mass functionalization and enhanced partitioning. CFD simulations show that concentrations of ozone and primary product change near the human surface rapidly, indicating non-uniform spatial distributions from the occupant surface to ambient air, while secondary ozone product is relatively well-mixed throughout the room. This development establishes a framework to integrate different modeling tools and experimental measurements, opening up an avenue for development of comprehensive and integrated models with representations of various chemistry in indoor environments.

Mechanisms and competition of halide substitution and hydrolysis in reactions of N2O5 with seawater

S_N2-type halide substitution and hydrolysis are two of the most ubiquitous reactions in chemistry. The interplay between these processes is fundamental in atmospheric chemistry through reactions of N₂O₅ and seawater. N₂O₅ plays a major role in regulating levels of O₃, OH, NO _x , and CH₄. While the reactions of N₂O₅ and seawater are of central importance, little is known about their mechanisms. Of interest is the activation of Cl in seawater by the formation of gaseous ClNO₂, which occurs despite the fact that hydrolysis (to HNO₃) is energetically more favorable. We determine key features of the reaction landscape that account for this behavior in a theoretical study of the cluster N₂O₅/Cl^-/H₂O. This was carried out using ab initio molecular dynamics to determine reaction pathways, structures, and time scales. While hydrolysis of N₂O₅ occurs in the absence of Cl^-, results here reveal that a low-lying pathway featuring halide substitution intermediates enhances hydrolysis.

Ion reactions in atmospherically-relevant clusters: mechanisms, dynamics and spectroscopic signatures

Exploring models of reactions of N₂O₄ with ions in water in order to provide molecular-level understanding of these processes.

Elucidating the molecular mechanisms of Criegee-amine chemistry in the gas phase and aqueous surface environments

Computational results suggest that the reactions ofantisubstituted Criegee intermediates with amine could lead to oligomers, which may play an important role in new particle formation and hydroxyl radical generation in the troposphere.

Single-Molecule Catalysis Revealed: Elucidating the Mechanistic Framework for the Formation and Growth of Atmospheric Iodine Oxide Aerosols in Gas-Phase and Aqueous Surface Environments

Iodine oxide aerosols are ubiquitous in many coastal atmospheric environments. However, the exact mechanism responsible for their homogeneous nucleation and subsequent cluster growth remains to be fully established. Using quantum chemical calculations, we propose a new mechanistic framework for the formation and subsequent growth of iodine oxide aerosols, which takes advantage of noncovalent interactions between iodine oxides (I₂O₅ and I₂O₄) and iodine acids (HIO₃ and HIO₂). Larger iodine oxide clusters are suggested to be formed in a facile manner and with enhanced exothermicity. The newly proposed mechanisms follow both concerted and stepwise pathways. In all these new chemistries, an O:I ratio of 2-2.5 is predicted, which satisfies an experimentally derived criterion recently proposed for identifying iodine oxides involved in atmospheric aerosol formation. Born-Oppenheimer molecular dynamics simulations at the air-water interface suggest that I₂O₅ and I₄O₁₀, which are two of the most common nucleating iodine oxides, react with interfacial water on the picosecond time scale and result in novel nucleating species such as H₂I₂O₆ and HI₄O₁₁^- or I₃O₈. An important implication of these simulation results is that aqueous surfaces, which are ubiquitous in the atmosphere, may activate iodine oxides to result in a new class of nucleating compounds, which can form mixed aerosol particles with potent precursors, such as HIO₃ or H₂SO₄, in marine air masses via typical acid-based interactions. Overall, these results give a better understanding of iodine-rich aerosols in diverse environments.

Formation of HONO from the NH3-promoted hydrolysis of NO2 dimers in the atmosphere

One challenging issue in atmospheric chemistry is identifying the source of nitrous acid (HONO), which is believed to be a primary source of atmospheric "detergent" OH radicals. Herein, we show a reaction route for the formation of HONO species from the NH₃-promoted hydrolysis of a NO₂ dimer (ONONO₂), which entails a low free-energy barrier of 0.5 kcal/mol at room temperature. Our systematic study of HONO formation based on NH₃ + ONONO₂ + nH₂O and water droplet systems with the metadynamics simulation method and a reaction pathway searching method reveals two distinct mechanisms: (i) In monohydrates (n = 1), tetrahydrates (n = 4), and water droplets, only one water molecule is directly involved in the reaction (denoted the single-water mechanism); and (ii) the splitting of two neighboring water molecules is seen in the dihydrates (n = 2) and trihydrates (n = 3) (denoted the dual-water mechanism). A comparison of the computed free-energy surface for NH₃-free and NH₃-containing systems indicates that gaseous NH₃ can markedly lower the free-energy barrier to HONO formation while stabilizing the product state, producing a more exergonic reaction, in contrast to the endergonic reaction for the NH₃-free system. More importantly, the water droplet reduces the free-energy barrier for HONO formation to 0.5 kcal/mol, which is negligible at room temperature. We show that the entropic contribution is important in the mechanism by which NH₃ promotes HONO formation. This study provides insight into the importance of fundamental HONO chemistry and its broader implication to aerosol and cloud processing chemistry at the air-water interface.

Barrierless Proton Transfer in the Weak C-H···O Hydrogen Bonded Methacrolein Dimer upon Nonresonant Multiphoton Ionization in the Gas Phase

Intermolecular proton transfer (IMPT) in a C-H···O hydrogen bonded dimer of an α,β-unsaturated aldehyde, methacrolein (MC), upon nonresonant multiphoton ionization by 532 nm laser pulses (10 ns), has been investigated using time-of-flight (TOF) mass spectrometry under supersonic cooling condition. The mass peaks corresponding to both the protonated molecular ion [(MC)H⁺] and intact dimer cation [(MC)₂]^•+ show up in the mass spectra, and the peak intensity of the former increases proportionately with the latter with betterment of the jet cooling conditions. The observations indicate that [(MC)₂]^•+ is the likely precursor of (MC)H⁺ and, on the basis of electronic structure calculations, IMPT in the dimer cation has been shown to be the key reaction for formation of the latter. Laser power dependences of ion yields indicate that at this wavelength the dimer is photoionized by means of 4-photon absorption process, and the total 4-photon energy is nearly the same as the predicted vertical ionization energy of the dimer. Electronic structure calculations reveal that the optimized structures of [(MC)₂]^•+ correspond to a proton transferred configuration wherein the aldehydic hydrogen is completely shifted to the carbonyl oxygen of the neighboring moiety. Potential energy scans along the C-H···O coordinate also show that the IMPT in [(MC)₂]^•+ is a barrierless process.

Insight into Chemistry on Cloud/Aerosol Water Surfaces

Cloud/aerosol water surfaces exert significant influence over atmospheric chemical processes. Atmospheric processes at the water surface are observed to follow mechanisms that are quite different from those in the gas phase. This Account summarizes our recent findings of new reaction pathways on the water surface. We have studied these surface reactions using Born-Oppenheimer molecular dynamics simulations. These studies provide useful information on the reaction time scale, the underlying mechanism of surface reactions, and the dynamic behavior of the product formed on the aqueous surface. According to these studies, the aerosol water surfaces confine the atmospheric species into a specific orientation depending on the hydrophilicity of atmospheric species or the hydrogen-bonding interactions between atmospheric species and interfacial water. As a result, atmospheric species are activated toward a particular reaction on the aerosol water surface. For example, the simplest Criegee intermediate (CH₂OO) exhibits high reactivity toward the interfacial water and hydrogen sulfide, with the reaction times being a few picoseconds, 2-3 orders of magnitude faster than that in the gas phase. The presence of interfacial water molecules induces proton-transfer-based stepwise pathways for these reactions, which are not possible in the gas phase. The strong hydrophobicity of methyl substituents in larger Criegee intermediates (>C1), such as CH₃CHOO and (CH₃)₂COO, blocks the formation of the necessary prereaction complexes for the Criegee-water reaction to occur at the water droplet surface, which lowers their proton-transfer ability and hampers the reaction. The aerosol water surface provides a solvent medium for acids (e.g., HNO₃ and HCOOH) to participate in reactions via mechanisms that are different from those in the gas and bulk aqueous phases. For example, the anti-CH₃CHOO-HNO₃ reaction in the gas phase follows a direct reaction between anti-CH₃CHOO and HNO₃, whereas on a water surface, the HNO₃-mediated stepwise hydration of anti-CH₃CHOO is dominantly observed. The high surface/volume ratio of interfacial water molecules at the aerosol water surface can significantly lower the energy barriers for the proton transfer reactions in the atmosphere. Such catalysis by the aerosol water surface is shown to cause the barrier-less formation of ammonium bisulfate from hydrated NH₃ and SO₃ molecules rather than from the reaction of H₂SO₄ with NH₃. Finally, an aerosol water droplet is a polar solvent, which would favorably interact with high polarity substrates. This can accelerate interconversion of different conformers (e.g., anti and syn) of atmospheric species, such as glyoxal, depending on their polarity. The results discussed here enable an improved understanding of atmospheric processes on the aerosol water surface.