Molecular Dynamics Simulations of Atmospheric Oxidants at the Air−Water Interface: Solvation and Accommodation of OH and O3

Atmospheric Intermediates at the Air-Water Interface

The air-water interface (AWI) is a ubiquitous reaction field different from the bulk phase where unexpected reactions and physical processes often occur. The AWI is a region where air contacts cloud droplets, aerosol particles, the ocean surface, and biological surfaces such as fluids that line human epithelia. In Earth's atmosphere, short-lived intermediates are expected to be generated at the AWI during multiphase reactions. Recent experimental developments have enabled the direct detection of atmospherically relevant, short-lived intermediates at the AWI. For example, spray ionization mass spectrometric analysis of water microjets exposed to a gaseous mixture of ozone and water vapor combined with a 266 nm laser flash photolysis system (LFP-SIMS) has been used to directly probe organic peroxyl radicals (RO2·) produced by interfacial hydroxyl radicals (OH·) + organic compound reactions. OH· emitted immediately after the laser flash photolysis of carboxylic acid at the gas-liquid interface have been directly detected by time-resolved, laser-induced florescence techniques that can be used to study atmospheric multiphase photoreactions. In this Featured Article, we show some recent experimental advances in the detection of atmospherically important intermediates at the AWI and the associated reaction mechanisms. We also discuss current challenges and future prospects for atmospheric multiphase chemistry.

Chemical Kinetics in Microdroplets

Micrometer-sized compartments play significant roles in driving heterogeneous transformations within atmospheric and biochemical systems as well as providing vehicles for drug delivery and novel reaction environments for the synthesis of industrial chemicals. Many reports now indicate that reaction kinetics are accelerated under microconfinement, for example, in sprays, thin films, droplets, aerosols, and emulsions. These observations are dramatic, posing a challenge to our understanding of chemical reaction mechanisms with potentially significant practical consequences for predicting the complex chemistry in natural systems. Here we introduce the idea of kinetic confinement, which is intended to provide a conceptual backdrop for understanding when and why microdroplet reaction kinetics differ from their macroscale analogs.

On the Statistical Mechanics of Mass Accommodation at Liquid-Vapor Interfaces

We propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

Iodide oxidation by ozone at the surface of aqueous microdroplets

The oxidation of iodide by ozone occurs at the sea-surface and within sea spray aerosol, influencing the overall ozone budget in the marine boundary layer and leading to the emission of reactive halogen gases. A detailed account of the surface mechanism has proven elusive, however, due to the difficulty in quantifying multiphase kinetics. To obtain a clearer understanding of this reaction mechanism at the air-water interface, we report pH-dependent oxidation kinetics of I- in single levitated microdroplets as a function of [O3] using a quadrupole electrodynamic trap and an open port sampling interface for mass spectrometry. A kinetic model, constrained by molecular simulations of O3 dynamics at the air-water interface, is used to understand the coupled diffusive, reactive, and evaporative pathways at the microdroplet surface, which exhibit a strong dependence on bulk solution pH. Under acidic conditions, the surface reaction is limited by O3 diffusion in the gas phase, whereas under basic conditions the reaction becomes rate limited on the surface. The pH dependence also suggests the existence of a reactive intermediate IOOO- as has previously been observed in the Br- + O3 reaction. Expressions for steady-state surface concentrations of reactants are derived and utilized to directly compute uptake coefficients for this system, allowing for an exploration of uptake dependence on reactant concentration. In the present experiments, reactive uptake coefficients of O3 scale weakly with bulk solution pH, increasing from 4 × 10-4 to 2 × 10-3 with decreasing solution pH from pH 13 to pH 3.

Mass Transfer through Vapor-Liquid Interfaces Studied by Non-Stationary Molecular Dynamics Simulations

Molecular dynamics (MD) simulations are highly attractive for studying the influence of interfacial effects, such as the enrichment of components, on the mass transfer through the interface. In a recent work, we have presented a steady-state MD simulation method for investigating this phenomenon and tested it using model mixtures with and without interfacial enrichment. The present study extends this work by introducing a non-stationary MD simulation method. A rectangular simulation box that contains a mixture of two components 1 + 2 with a vapor phase in the middle and two liquid phases on both sides is used. Starting from a vapor-liquid equilibrium state, a non-stationary molar flux of component 2 is induced by inserting particles of component 2 into the center of the vapor phase in a pulse-like manner. During the isothermal relaxation process, particles of component 2 pass through the vapor phase, cross the vapor-liquid interface, and enter the liquid phase. The system thereby relaxes into a new vapor-liquid equilibrium state. During the relaxation process, spatially resolved responses for the component densities, fluxes, and pressure are sampled. To reduce the noise and provide measures for the uncertainty of the observables, a set of replicas of simulations is carried out. The new simulation method was applied to study mass transfer in two binary Lennard-Jones mixtures: one that exhibits a strong enrichment of the low-boiling component 2 at the vapor-liquid interface and one that shows no enrichment. Even though both mixtures have similar transport coefficients in the bulk phases, the results for mass transfer differ significantly, indicating that the interfacial enrichment influences the mass transfer.

Phase State and Relative Humidity Regulate the Heterogeneous Oxidation Kinetics and Pathways of Organic-Inorganic Mixed Aerosols

Inorganic species always coexist with organic materials in atmospheric particles and may influence the heterogeneous oxidation of organic aerosols. However, very limited studies have explored the role of the inorganics in the chemical evolution of organic species in mixed aerosols. This study examines the heterogeneous oxidation of glutaric acid-ammonium sulfate and 1,2,6-hexanetriol-ammonium sulfate aerosols by hydroxyl radicals (OH) under varied organic mass fractions (f_org) and relative humidity in a flow tube reactor. Coupling the oxidation kinetics and product measurements with kinetic model simulations, we found that under both low relative humidity (RH, 30-35%) and high RH conditions (85%), the decreased f_org from 0.7 to 0.2 accelerates the oxidation of the organic materials by a factor of up to 11. We suggest that the faster oxidation kinetics under low-RH conditions is due to full or partial phase separation, with the organics greatly enriched at the particle outer region, while enhanced "salting-out" of the organics and OH adsorption caused by higher inorganics could explain the observations under high-RH conditions. Analysis of the oxidation products reveals that the dilution of organics by the inorganic salts and corresponding water uptake under high-RH conditions will favor alkoxy radical fragmentation by a factor of 3-4 and inhibit its secondary chain propagation chemistry. Our results suggest that atmospheric organic aerosol oxidation lifetime and composition are strongly impacted by the coexistent inorganic salts.

Spontaneous Formation of Hydrogen Peroxide in Water Microdroplets

There is accumulating evidence that many chemical reactions are accelerated by several orders of magnitude in micrometer-sized aqueous or organic liquid droplets compared to their corresponding bulk liquid phase. However, the molecular origin of the enhanced rates remains unclear as in the case of spontaneous appearance of 1 μM hydrogen peroxide in water microdroplets. In this Letter, we consider the range of ionization energies and whether interfacial electric fields of a microdroplet can feasibly overcome the high energy step from hydroxide ions (OH^-) to hydroxyl radicals (OH^•) in a primary H₂O₂ mechanism. We find that the vertical ionization energies (VIEs) of partially solvated OH^- ions are greatly lowered relative to the average VIE in the bulk liquid, unlike the case of the Cl^- anion which shows no reduction in the VIEs regardless of solvation environment. Overall reduced hydrogen-bonding and undercoordination of OH^- are structural features that are more readily present at the air-water interface, where the energy scale for ionization can be matched by statistically probable electric field values.

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O₃. However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. To examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients (trans-aconitic acid (C₆H₆O₆), maleic acid (C₄H₄O₄), and sodium nitrite (NaNO₂)). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O₃ partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O₃ and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O₃, compared to gas and bulk phases is rapidly maintained and remains constant even as O₃ is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O₃ loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O₃ arrival in the bulk by diffusion. As a result, bulk-phase O₃ becomes depleted from its Henry's law solubility. This bulk-phase O₃ depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k_rxn ≈ 10³-10⁴ M^-1 s^-1), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k_rxn ≈ 10⁵ M^-1 s^-1) where surface reactions strongly impact the observed decay kinetics.

In situ analysis of the bulk and surface chemical compositions of organic aerosol particles

Understanding the chemical and physical properties of particles is an important scientific, engineering, and medical issue that is crucial to air quality, human health, and environmental chemistry. Of special interest are aerosol particles floating in the air for both indoor virus transmission and outdoor atmospheric chemistry. The growth of bio- and organic-aerosol particles in the air is intimately correlated with chemical structures and their reactions in the gas phase at aerosol particle surfaces and in-particle phases. However, direct measurements of chemical structures at aerosol particle surfaces in the air are lacking. Here we demonstrate in situ surface-specific vibrational sum frequency scattering (VSFS) to directly identify chemical structures of molecules at aerosol particle surfaces. Furthermore, our setup allows us to simultaneously probe hyper-Raman scattering (HRS) spectra in the particle phase. We examined polarized VSFS spectra of propionic acid at aerosol particle surfaces and in particle bulk. More importantly, the surface adsorption free energy of propionic acid onto aerosol particles was found to be less negative than that at the air/water interface. These results challenge the long-standing hypothesis that molecular behaviors at the air/water interface are the same as those at aerosol particle surfaces. Our approach opens a new avenue in revealing surface compositions and chemical aging in the formation of secondary organic aerosols in the atmosphere as well as chemical analysis of indoor and outdoor viral aerosol particles.

Representation of Multiphase OH Oxidation of Amorphous Organic Aerosol for Tropospheric Conditions

Organic aerosol (OA) is ubiquitous in the atmosphere and, during transport, can experience chemical transformation with consequences for air quality and climate. Prediction of the chemical evolution of OA depends on its reactivity with atmospheric oxidants such as the OH radical. OA particles undergo amorphous phase transitions from liquid to solid (glassy) states in response to temperature changes, which, in turn, will impact its reactivity toward OH oxidation. To improve the predictability of OA reactivity toward OH oxidation, the reactive uptake coefficients (γ) of OH radicals reacting with triacontane and squalane serving as amorphous OA surrogates were measured at temperatures from 213-293 K. γ increases strongest with temperature when the organic species is in the liquid phase, compared to when being in the semisolid or solid phase. The resistor model is applied, accounting for the amorphous phase state changes using the organic species' glass transition temperature and fragility, to evaluate the physicochemical parameters of the temperature dependent OH uptake process. This allows for the derivation of a semiempirical formula, applicable to models, to predict the degree of oxidation and chemical lifetime of the condensed-phase organic species for typical tropospheric temperature and humidity when OA particle viscosity is known.

The mobility and solvation structure of a hydroxyl radical in a water nanodroplet: a Born-Oppenheimer molecular dynamics study

Hydroxyl radicals (OH*) play a crucial role in atmospheric chemistry and biological processes. In this study, Born-Oppenheimer molecular dynamics simulations are performed under ambient conditions for a hydroxyl radical in a water nanodroplet containing 191 water molecules. Density functional theory calculations are performed at the BLYP-D3 level with some test calculations at the B3LYP-D3 level. In two 150 ps trajectories, either with OH* initially located in the interior region or at the surface of the water nanodroplet, the OH* radical ends up in the subsurface layer of the nanodroplet, which is different from the "surface preference" predicted from previous empirical force field simulations. The solvation structure of OH* contains fluctuating hydrogen bonds, as well as a two-center three-electron hemibond in some cases. The mobility of OH* is enhanced by hydrogen transfer, which has a free energy barrier of ∼4.6 kcal mol-1. The results presented in this study deepen our understanding of the structure and dynamics of OH* in aqueous solutions, especially around the air-water interface.

Reactive Uptake of Ozone to Simulated Seawater: Evidence for Iodide Depletion

The reaction of ozone with iodide in the ocean is a major ozone dry deposition pathway, as well as an important source of reactive iodine to the marine troposphere. Few prior laboratory experiments have been conducted with environmentally relevant ozone mixing ratios and iodide concentrations, leading to uncertainties in the rate of the reaction under marine boundary layer conditions. As well, there remains disagreement in the literature assessment of the relative contributions of an interfacial reaction via ozone adsorbed to the ocean surface versus a bulk reaction with dissolved ozone. In this study, we measure the uptake coefficient of ozone over a buffered, pH 8 salt solution replicating the concentrations of iodide, bromide, and chloride in the ocean over an ozone mixing ratio of 60-500 ppb. Due to iodide depletion in the solution, the measured ozone uptake coefficient is dependent on the exposure time of the solution to ozone and its mixing ratio. A kinetic multilayer model confirms that iodide depletion is occurring not only within ozone's reactodiffusive depth, which is on the order of microns for environmental conditions, but also deeper into the solution as well. Best model-measurement agreement arises when some degree of nondiffusive mixing is occurring in the solution, transporting iodide from deeper in the solution to a thin, diffusively mixed upper layer. If such mixing occurs rapidly in the environment, iodide depletion is unlikely to reduce ozone dry deposition rates. Unrealistically high bulk-to-interface partitioning of iodide is required for the model to predict a substantial interfacial component to the reaction, indicating that the Langmuir-Hinshelwood mechanism is not dominant under environmental conditions.

Photoinduced Oxidation Reactions at the Air-Water Interface

Chemistry on water is a fascinating area of research. The surface of water and the interfaces between water and air or hydrophobic media represent asymmetric environments with unique properties that lead to unexpected solvation effects on chemical and photochemical processes. Indeed, the features of interfacial reactions differ, often drastically, from those of bulk-phase reactions. In this Perspective, we focus on photoinduced oxidation reactions, which have attracted enormous interest in recent years because of their implications in many areas of chemistry, including atmospheric and environmental chemistry, biology, electrochemistry, and solar energy conversion. We have chosen a few representative examples of photoinduced oxidation reactions to focus on in this Perspective. Although most of these examples are taken from the field of atmospheric chemistry, they were selected because of their broad relevance to other areas. First, we outline a series of processes whose photochemistry generates hydroxyl radicals. These OH precursors include reactive oxygen species, reactive nitrogen species, and sulfur dioxide. Second, we discuss processes involving the photooxidation of organic species, either directly or via photosensitization. The photochemistry of pyruvic acid and fatty acid, two examples that demonstrate the complexity and versatility of this kind of chemistry, is described. Finally, we discuss the physicochemical factors that can be invoked to explain the kinetics and thermodynamics of photoinduced oxidation reactions at aqueous interfaces and analyze a number of challenges that need to be addressed in future studies.

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.

Effect of Bulk Composition on the Heterogeneous Oxidation of Semi-Solid Atmospheric Aerosols

The OH-initiated heterogeneous oxidation of semi-solid saccharide particles with varying bulk compositions was investigated in an atmospheric pressure flow tube at 30% relative humidity. Reactive uptake coefficients were determined from the rate loss of the saccharide reactants measured by mass spectrometry at different monosaccharide (methyl-β-d-glucopyranoside, C7H14O6) and disaccharide (lactose, C12H22O11) molar ratios. The reactive uptake for the monosaccharide was found to decrease from 0.53 ± 0.10 to 0.05 ± 0.06 as the mono-to-disaccharide molar ratio changed from 8:1 to 1:1. A reaction–diffusion model was developed in order to determine the effect of chemical composition on the reactive uptake. The observed decays can be reproduced using a Vignes relationship to predict the composition dependence of the reactant diffusion coefficients. The experimental data and model results suggest that the addition of the disaccharide significantly increases the particle viscosity leading to slower mass transport phenomena from the bulk to the particle surface and to a decreased reactivity. These findings illustrate the impact of bulk composition on reactant bulk diffusivity which determines the rate-limiting step during the chemical transformation of semi-solid particles in the atmosphere.

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.

Interfacial vs Bulk Ozonolysis of Nerolidol

Ozone readily reacts with olefins with the formation of more reactive Criegee intermediates (CIs). The transient CIs impact HO _x cycles, and they play a role in new particle formation in the troposphere. Oxidation by O₃ occurs both in the gas-phase, in the liquid phase, and at air-water and air-aerosol interfaces. In light of the importance of O₃ in environmental and engineered chemical transformations, we have investigated the ozonolysis mechanisms of a triolefin C₁₅-alcohol, nerolidol (Nero, a biogenic sesquiterpene), at the air-water interface in the presence of acetonitrile. Surface-sensitive pneumatic ionization mass spectrometric detection of α-hydroxy-hydroperoxides and functionalized carboxylates, generated by the hydration and isomerization of CIs, respectively, enables us to evaluate the relative reactivity of each C=C toward O₃. In addition, we compare bulk-phase ozonolysis chemistry to similar reactions taking place at the air-water interface. Our experimental results show that O₃ reacts primarily with the (CH₃)₂C=CH- and -(CH₃)C=CH- moieties (>∼98%), while the O₃ attack on the terminal -HC=CH₂ site (<∼2%) is a minor pathway during both interfacial and bulk ozonolysis. The presence of functionalized-carboxylates on interfaces but not in bulk-phase reactions with O₃ indicates that the isomerization of the CIs is not hindered at the air-water interface due to the lower availability of water .

Mass Transfer Thermodynamics through a Gas-Liquid Interface

Molecular level information about thermodynamic variations (enthalpy, entropy, and free energy) of a gas molecule as it crosses a gas-liquid interface is strongly lacking from an experimental perspective under equilibrium conditions. Herein, we perform in situ measurements of water interacting with the ionic liquid (IL) 1-butyl-3-methylimidazolium acetate, [C₄mim][Ace], using ambient pressure X-ray photoelectron spectroscopy in order to assess the interfacial uptake of water quantitatively as a function of temperature, pressure, and water mole fraction ( x_w). The surface spectroscopy results are compared to existing bulk water absorption experiments, showing that the amount of water in the interfacial region is consistently greater than that in the bulk. The enthalpy and entropy of water sorption vary significantly between the gas-liquid interface and the bulk as a function of x_w, with a crossover that occurs near x_w = 0.6 where the water-IL mixture converts from being homogeneous ( x_w < 0.6) to nanostructured ( x_w > 0.6). Free energy results reveal that water at the gas-IL interface is thermodynamically more favorable than that in the bulk, consistent with the enhanced water concentration in the interfacial region. The results herein show that the efficacy for an ionic liquid to absorb a gas phase molecule is not merely a function of bulk solvation parameters but also is significantly influenced by the thermodynamics occurring across the gas-IL interface during the mass transfer process.

OH-Radical Oxidation of Lung Surfactant Protein B on Aqueous Surfaces

Air pollutants generate reactive oxygen species on lung surfaces. Here we report how hydroxyl radicals (·OH) injected on the surface of water react with SP-B_1-25, a 25-residue polypeptide surrogate of human lung surfactant protein B. Our experiments consist of intersecting microjets of aqueous SP-B_1-25 solutions with O₃/O₂/H₂O/N₂(g) gas streams that are photolyzed into ·OH(g) in situ by 266 nm laser nanosecond pulses. Surface-sensitive mass spectrometry enables us to monitor the prompt (<10 μs) and simultaneous formation of primary O _n -containing products/intermediates (n≤5) triggered by the reaction of ·OH with interfacial SP-B_1-25. We found that O-atoms from both O₃ and ·OH are incorporated into the reactive cysteine Cys₈ and Cys₁₁ and tryptophan Trp₉ components of the hydrophobic N-terminus of SP-B_1-25 that lies at the topmost layers of the air-liquid interface. Remarkably, these processes are initiated by ·OH additions rather than by H-atom abstractions from S-H, C-H, or N-H groups. By increasing the hydrophilicity of the N-terminus region of SP-B_1-25, these transformations will impair its role as a surfactant.

An Overview of Dynamic Heterogeneous Oxidations in the Troposphere

Due to the adverse effect of atmospheric aerosols on public health and their ability to affect climate, extensive research has been undertaken in recent decades to understand their sources and sinks, as well as to study their physical and chemical properties. Atmospheric aerosols are important players in the Earth’s radiative budget, affecting incoming and outgoing solar radiation through absorption and scattering by direct and indirect means. While the cooling properties of pure inorganic aerosols are relatively well understood, the impact of organic aerosols on the radiative budget is unclear. Additionally, organic aerosols are transformed through chemical reactions during atmospheric transport. The resulting complex mixture of organic aerosol has variable physical and chemical properties that contribute further to the uncertainty of these species modifying the radiative budget. Correlations between oxidative processing and increased absorptivity, hygroscopicity, and cloud condensation nuclei activity have been observed, but the mechanisms behind these phenomena have remained unexplored. Herein, we review environmentally relevant heterogeneous mechanisms occurring on interfaces that contribute to the processing of aerosols. Recent laboratory studies exploring processes at the aerosol–air interface are highlighted as capable of generating the complexity observed in the environment. Furthermore, a variety of laboratory methods developed specifically to study these processes under environmentally relevant conditions are introduced. Remarkably, the heterogeneous mechanisms presented might neither be feasible in the gas phase nor in the bulk particle phase of aerosols at the fast rates enabled on interfaces. In conclusion, these surface mechanisms are important to better understand how organic aerosols are transformed in the atmosphere affecting the environment.

Kinetics Study of OH Uptake onto Deliquesced NaCl Particles by Combining Laser Photolysis and Laser-Induced Fluorescence

Despite the role of hydroxyl radical (OH) uptake onto sea-salt particles as a daytime chlorine source, affecting the chemical processes in the marine boundary layer, its uptake coefficient has not yet been confirmed by direct measurement methods. This study reports the application of a combination technique of laser flash photolysis generation and laser-induced fluorescence detection for the direct kinetic measurement of OH uptake onto deliquesced NaCl particles. The uptake coefficient was not constant and inversely depended on the initial OH concentration, indicating that the first uptake step is Langmuir-type adsorption. The resistance model, including surface processes, well reproduced the observed uptake coefficient. The model predicted an uptake coefficient for the atmospheric relevant OH concentration within the range from 0.77 to 0.95. Such values may lead to emissions of Cl₂ higher than those predicted in previous studies based on other values. Hence, the proposed value may provide more reliable estimations of ozone formation, oxidation of volatile organic compounds, secondary organic aerosol formation, and lifetime of methane and elemental mercury in the marine boundary layer.

Rapid Water Transport through Organic Layers on Ice

Processes involving atmospheric aerosol and cloud particles are affected by condensation of organic compounds that are omnipresent in the atmosphere. On ice particles, organic compounds with hydrophilic functional groups form hydrogen bonds with the ice and orient their hydrophobic groups away from the surface. The organic layer has been expected to constitute a barrier to gas uptake, but recent experimental studies suggest that the accommodation of water molecules on ice is only weakly affected by condensed short-chain alcohol layers. Here, we employ molecular dynamics simulations to study the water interactions with n-butanol covered ice at 200 K and show that the small effect of the condensed layer is due to efficient diffusion of water molecules along the surface plane while seeking appropriate sites to penetrate, followed by penetration driven by the combined attractive forces from butanol OH groups and water molecules within the ice. The water molecules that penetrate through the n-butanol layer become strongly bonded by approximately three hydrogen bonds at the butanol-ice interface. The obtained accommodation coefficient (0.81 ± 0.03) is in excellent agreement with results from previous environmental molecular beam experiments, leading to a picture where an adsorbed n-butanol layer does not alter the apparent accommodation coefficient but dramatically changes the detailed molecular dynamics and kinetics.

Diffusive confinement of free radical intermediates in the OH radical oxidation of semisolid aerosols

Multiphase chemical reactions (gas + solid/liquid) involve a complex interplay between bulk and interface chemistry, diffusion, evaporation, and condensation. Reactions of atmospheric aerosols are an important example of this type of chemistry: the rich array of particle phase states and multiphase transformation pathways produce diverse but poorly understood interactions between chemistry and transport. Their chemistry is of intrinsic interest because of their role in controlling climate. Their characteristics also make them useful models for the study of principles of reactivity of condensed materials under confined conditions. In previous work, we have reported a computational study of the oxidation chemistry of a liquid aliphatic aerosol. In this study, we extend the calculations to investigate nearly the same reactions at a semisolid gas-aerosol interface. A reaction-diffusion model for heterogeneous oxidation of triacontane by hydroxyl radicals (OH) is described, and its predictions are compared to measurements of aerosol size and composition, which evolve continuously during oxidation. These results are also explicitly compared to those obtained for the corresponding liquid system, squalane, to pinpoint salient elements controlling reactivity. The diffusive confinement of the free radical intermediates at the interface results in enhanced importance of a few specific chemical processes such as the involvement of aldehydes in fragmentation and evaporation, and a significant role of radical-radical reactions in product formation. The simulations show that under typical laboratory conditions semisolid aerosols have highly oxidized nanometer-scale interfaces that encapsulate an unreacted core and may confer distinct optical properties and enhanced hygroscopicity. This highly oxidized layer dynamically evolves with reaction, which we propose to result in plasticization. The validated model is used to predict chemistry under atmospheric conditions, where the OH radical concentration is much lower. The oxidation reactions are more strongly influenced by diffusion in the particle, resulting in a more liquid-like character.

N2O5at water surfaces: binding forces, charge separation, energy accommodation and atmospheric implications

Studying the interactions between N₂O₅and water in nano-sized clusters, in bulk and on the surface of water.

Heterogeneous OH Oxidation, Shielding Effects, and Implications for the Atmospheric Fate of Terbuthylazine and Other Pesticides

Terbuthylazine (TBA) is a widely used herbicide, and its heterogeneous reaction with OH radicals is important for assessing its potential to undergo atmospheric long-range transport and to affect the environment and public health. The apparent reaction rate coefficients obtained in different experimental investigations, however, vary by orders of magnitude depending on the applied experimental techniques and conditions. In this study, we used a kinetic multilayer model of aerosol chemistry with reversible surface adsorption and bulk diffusion (KM-SUB) in combination with a Monte Carlo genetic algorithm to simulate the measured decay rates of TBA. Two experimental data sets available from different studies can be described with a consistent set of kinetic parameters resolving the interplay of chemical reaction, mass transport, and shielding effects. Our study suggests that mass transport and shielding effects can substantially extend the atmospheric lifetime of reactive pesticides from a few days to weeks, with strong implications for long-range transport and potential health effects of these substances.

Molecular dynamics simulation of the local concentration and structure in multicomponent aerosol nanoparticles under atmospheric conditions

MD simulations predicted core–shell or partially engulfed morphologies (depending on the type of the organic compound present) in multicomponent aerosol nanoparticles.

A theoretical study on the reaction of ozone with aqueous iodide

Atmospheric iodine chemistry plays a key role in tropospheric ozone catalytic destruction, new particle formation, and as one of the possible sinks of gaseous polar elemental mercury. Moreover, it has been recently proposed that reaction of ozone with iodide on the sea surface could be the major contributor to the chemical loss of atmospheric ozone. However, the mechanism of the reaction between aqueous iodide and ozone is not well known. The aim of this paper is to improve the understanding of such a mechanism. In this paper, an ab initio study of the reaction of aqueous iodide and ozone is presented, evaluating thermodynamic data of the different reactions proposed in previous experimental studies. In addition, the structures, energetics and possible evolution of the key IOOO(-) intermediate are discussed for the first time.

"Sizing" Heterogeneous Chemistry in the Conversion of Gaseous Dimethyl Sulfide to Atmospheric Particles

The oxidation of biogenic dimethyl sulfide (DMS) emissions is a global source of cloud condensation nuclei. The amounts of the nucleating H2SO4(g) species produced in such process, however, remain uncertain. Hydrophobic DMS is mostly oxidized in the gas phase into H2SO4(g) + DMSO(g) (dimethyl sulfoxide), whereas water-soluble DMSO is oxidized into H2SO4(g) in the gas phase and into SO4(2-) + MeSO3(-) (methanesulfonate) on water surfaces. R = MeSO3(-)/(non-sea-salt SO4(2-)) ratios would therefore gauge both the strength of DMS sources and the extent of DMSO heterogeneous oxidation if Rhet = MeSO3(-)/SO4(2-) for DMSO(aq) + ·OH(g) were known. Here, we report that Rhet = 2.7, a value obtained from online electrospray mass spectra of DMSO(aq) + ·OH(g) reaction products that quantifies the MeSO3(-) produced in DMSO heterogeneous oxidation on aqueous aerosols for the first time. On this basis, the inverse R dependence on particle radius in size-segregated aerosol collected over Syowa station and Southern oceans is shown to be consistent with the competition between DMSO gas-phase oxidation and its mass accommodation followed by oxidation on aqueous droplets. Geographical R variations are thus associated with variable contributions of the heterogeneous pathway to DMSO atmospheric oxidation, which increase with the specific surface area of local aerosols.

Super-Maxwellian helium evaporation from pure and salty water

Helium atoms evaporate from pure water and salty solutions in super-Maxwellian speed distributions, as observed experimentally and modeled theoretically. The experiments are performed by monitoring the velocities of dissolved He atoms that evaporate from microjets of pure water at 252 K and 4-8.5 molal LiCl and LiBr at 232-252 K. The average He atom energies exceed the flux-weighted Maxwell-Boltzmann average of 2RT by 30% for pure water and 70% for 8.5m LiBr. Classical molecular dynamics simulations closely reproduce the observed speed distributions and provide microscopic insight into the forces that eject the He atoms from solution. Comparisons of the density profile and He kinetic energies across the water-vacuum interface indicate that the He atoms are accelerated by He-water collisions within the top 1-2 layers of the liquid. We also find that the average He atom kinetic energy scales with the free energy of solvation of this sparingly soluble gas. This free-energy difference reflects the steeply decreasing potential of mean force on the He atoms in the interfacial region, whose gradient is the repulsive force that tends to expel the atoms. The accompanying sharp decrease in water density suppresses the He-water collisions that would otherwise maintain a Maxwell-Boltzmann distribution, allowing the He atom to escape at high energies. Helium is especially affected by this reduction in collisions because its weak interactions make energy transfer inefficient.

Ozone uptake on glassy, semi-solid and liquid organic matter and the role of reactive oxygen intermediates in atmospheric aerosol chemistry

Humidity-induced phase transition and formation of reactive oxygen intermediates are important processes in the heterogeneous ozonolysis of unsaturated organic compounds in the atmosphere.

The influence of the surface composition of mixed monolayer films on the evaporation coefficient of water

The influence of mixed component organic surface films on the evaporation rate of water from an aqueous droplet is reported.

OH-Radical Specific Addition to Glutathione S-Atom at the Air-Water Interface: Relevance to the Redox Balance of the Lung Epithelial Lining Fluid

Antioxidants in epithelial lining fluids (ELF) prevent inhaled air pollutants from reaching lung tissue. This process, however, may upset ELF's redox balance, which is deemed to be expressed by the ratio of the major antioxidant glutathione (GSH) to its putative oxidation product GSSG. Previously, we found that at physiological pH O3(g) rapidly oxidizes GS(2-)(aq) (but not GSH(-)) to GSO3(-) rather than GSSG. Here, we report that in moderately acidic pH ≤ 5 media ·OH(g) oxidizes GSH(-)(aq) to sulfenic GSOH(-), sulfinic GSO2(-), and sulfonic GSO3(-) acids via ·OH specific additions to reduced S-atoms. The remarkable specificity of ·OH on water versus its lack of selectivity in bulk water implicates an unprecedented steering process during [OH···GSH] interfacial encounters. Thus, both O3 and ·OH oxidize GSH to GSOH(-) under most conditions, and since GSOH(-) is reduced back to GSH in vivo by NADPH, redox balance may be in fact signaled by GSH/GSOH ratios.

Near-unity mass accommodation coefficient of organic molecules of varying structure

Atmospheric aerosol particles have a significant effect on global climate, air quality, and consequently human health. Condensation of organic vapors is a key process in the growth of nanometer-sized particles to climate relevant sizes. This growth is very sensitive to the mass accommodation coefficient α, a quantity describing the vapor uptake ability of the particles, but knowledge on α of atmospheric organics is lacking. In this work, we have determined α for four organic molecules with diverse structural properties: adipic acid, succinic acid, naphthalene, and nonane. The coefficients are studied using molecular dynamics simulations, complemented with expansion chamber measurements. Our results are consistent with α = 1 (indicating nearly perfect accommodation), regardless of the molecular structural properties, the phase state of the bulk condensed phase, or surface curvature. The results highlight the need for experimental techniques capable of resolving the internal structure of nanoparticles to better constrain the accommodation of atmospheric organics.

Spectroscopic signatures of ozone at the air-water interface and photochemistry implications

First-principles simulations suggest that additional OH formation in the troposphere can result from ozone interactions with the surface of cloud droplets. Ozone exhibits an affinity for the air-water interface, which modifies its UV and visible light spectroscopic signatures and photolytic rate constant in the troposphere. Ozone cross sections on the red side of the Hartley band (290- to 350-nm region) and in the Chappuis band (450-700 nm) are increased due to electronic ozone-water interactions. This effect, combined with the potential contribution of the O3 + hν → O((3)P) + O2(X(3)Σg(-)) photolytic channel at the interface, leads to an enhancement of the OH radical formation rate by four orders of magnitude. This finding suggests that clouds can influence the overall oxidizing capacity of the troposphere on a global scale by stimulating the production of OH radicals through ozone photolysis by UV and visible light at the air-water interface.