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

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.

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.

A numerical compass for experiment design in chemical kinetics and molecular property estimation

Kinetic process models are widely applied in science and engineering, including atmospheric, physiological and technical chemistry, reactor design, or process optimization. These models rely on numerous kinetic parameters such as reaction rate, diffusion or partitioning coefficients. Determining these properties by experiments can be challenging, especially for multiphase systems, and researchers often face the task of intuitively selecting experimental conditions to obtain insightful results. We developed a numerical compass (NC) method that integrates computational models, global optimization, ensemble methods, and machine learning to identify experimental conditions with the greatest potential to constrain model parameters. The approach is based on the quantification of model output variance in an ensemble of solutions that agree with experimental data. The utility of the NC method is demonstrated for the parameters of a multi-layer model describing the heterogeneous ozonolysis of oleic acid aerosols. We show how neural network surrogate models of the multiphase chemical reaction system can be used to accelerate the application of the NC for a comprehensive mapping and analysis of experimental conditions. The NC can also be applied for uncertainty quantification of quantitative structure-activity relationship (QSAR) models. We show that the uncertainty calculated for molecules that are used to extend training data correlates with the reduction of QSAR model error. The code is openly available as the Julia package KineticCompass.

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.

pH jump kinetics in colliding microdroplets: accelerated synthesis of azamonardine from dopamine and resorcinol

Recent studies report the dramatic acceleration of chemical reactions in micron-sized compartments. In the majority of these studies the exact acceleration mechanism is unknown but the droplet interface is thought to play a significant role. Dopamine reacts with resorcinol to form a fluorescent product azamonardine and is used as a model system to examine how droplet interfaces accelerate reaction kinetics. The reaction is initiated by colliding two droplets levitated in a branched quadrupole trap, which allows the reaction to be observed in individual droplets where the size, concentration, and charge are carefully controlled. The collision of two droplets produces a pH jump and the reaction kinetics are quantified optically and in situ by measuring the formation of azamonardine. The reaction was observed to occur 1.5 to 7.4 times faster in 9-35 micron droplets compared to the same reaction conducted in a macroscale container. A kinetic model of the experimental results suggests that the acceleration mechanism arises from both the more rapid diffusion of oxygen into the droplet, as well as increased reagent concentrations at the air-water interface.

Optical trapping and light scattering in atmospheric aerosol science

Aerosol particles are ubiquitous in the atmosphere, and currently contribute a large uncertainty to climate models. Part of the endeavour to reduce this uncertainty takes the form of improving our understanding of aerosol at the microphysical level, thus enabling chemical and physical processes to be more accurately represented in larger scale models. In addition to modeling efforts, there is a need to develop new instruments and methodologies to interrogate the physicochemical properties of aerosol. This perspective presents the development, theory, and application of optical trapping, a powerful tool for single particle investigations of aerosol. After providing an overview of the role of aerosol in Earth's atmosphere and the microphysics of these particles, we present a brief history of optical trapping and a more detailed look at its application to aerosol particles. We also compare optical trapping to other single particle techniques. Understanding the interaction of light with single particles is essential for interpreting experimental measurements. In the final part of this perspective, we provide the relevant formalism for understanding both elastic and inelastic light scattering for single particles. The developments discussed here go beyond Mie theory and include both how particle and beam shape affect spectra. Throughout the entirety of this work, we highlight numerous references and examples, mostly from the last decade, of the application of optical trapping to systems that are relevant to the atmospheric aerosol.

Exploring Gas-Liquid Reactions with Microjets: Lessons We Are Learning

Liquid water is all around us: at the beach, in a cloud, from a faucet, inside a spray tower, covering our lungs. It is fascinating to imagine what happens to a reactive gas molecule as it approaches the surface of water in these examples. Some incoming molecules may first be deflected away after colliding with an evaporating water molecule. Those that do strike surface H₂O or other surface species may bounce directly off or become momentarily trapped through hydrogen bonding or other attractive forces. The adsorbed gas molecule can then desorb immediately or instead dissolve, passing through the interfacial region and into the bulk, perhaps diffusing back to the surface and evaporating before reacting. Alternatively, it may react with solute or water molecules in the interfacial or bulk regions, and a reaction intermediate or the final product may then desorb into the gas phase. Building a "blow by blow" picture of these pathways is challenging for vacuum-based techniques because of the high vapor pressure of water. In particular, collisions within the thick vapor cloud created by evaporating molecules just above the surface scramble the trajectories and internal states of the gaseous target molecules, hindering construction of gas-liquid reaction mechanisms at the atomic scale that we strive to map out.The introduction of the microjet in 1988 by Faubel, Schlemmer, and Toennies opened up entirely new possibilities. Their inspired solution seems so simple: narrow the end of a glass tube to a diameter smaller than the mean free path of the vapor molecules and then push the liquid through the tube at speeds of a car on a highway. The narrow liquid stream creates a sparse vapor cloud, with water molecules spaced far enough apart that they and the reactant gases interact, at most, weakly. Experimentalists, however, confront a host of challenges: nozzle clogging, unstable jetting, searching for vacuum-compatible solutions, measuring low signal levels, and teasing out artifacts because the slender jet is the smallest surface in the vacuum chamber. In this Account, we describe lessons that we are learning as we explore gases (DCl, (HCOOH)₂, N₂O₅) reacting with water molecules and solute ions in the near-interfacial region of these fast-flowing aqueous microjets.

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.