Molecular Characterization of the Surface Excess Charge Layer in Droplets

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.

Limitations of Atomistic Molecular Dynamics to Reveal Ejection of Proteins from Charged Nanodroplets

Atomistic molecular dynamics (MD) is frequently used to unravel the mechanisms of macroion release from electrosprayed droplets. However, atomistic MD is currently feasible for only the smallest window of droplet sizes appearing at the end steps of a droplet's lifetime. The relevance of the observations made to the actual droplet evolution, which is much longer than the simulated sizes, has not been addressed yet in the literature. Here, we perform a systematic study of the desolvation mechanisms of poly(ethylene glycol) (PEG), protonated peptides of different compositions, and proteins, to (a) obtain insight into the charging mechanism of macromolecules in larger droplets than those that are currently amenable to atomistic MD and (b) examine whether currently used atomistic MD modeling can establish the extrusion mechanism of proteins from droplets. To mimic larger droplets that are not amenable to MD modeling, we scale down the systems, by simulating a large droplet size relative to the macromolecule. MD of PEG charging reveals that, above a critical droplet size, ions are available near the backbone of the macromolecule, but charging occurs only transiently by transfer of ions from the solvent to the macroion, while below the critical size, the capture of the ion from PEG has a lifetime sufficiently long for the extrusion of a charged PEG from the aqueous droplet. This is the first report of the role of droplet curvature in the relation between macroion conformation and charging. Simulations of protonated peptides with a high degree of hydrophobicity show that partial extrusion of a peptide from the droplet surface is rare relative to desolvation by drying-out. Different from what has been presented in the literature, we argue that atomistic MD simulations have not sufficiently established the extrusion mechanism of proteins from droplets and their charging mechanism. We also argue that release of highly charged proteins can occur at an earlier stage of a droplet's lifetime than predicted by atomistic MD. In this earlier stage, we emphasize the key role of jets emanating from a droplet at the point of charge-induced instability in the release of proteins.

A Molecular Dynamic Study of the Effects of Surface Partitioning on the OH Radical Interactions with Solutes in Multicomponent Aqueous Aerosols

The surface-bulk partitioning of small saccharide and amide molecules in aqueous droplets was investigated using molecular dynamics. The air-particle interface was modeled using a 80 Å cubic water box containing a series of organic molecules and surrounded by gaseous OH radicals. The properties of the organic solutes within the interface and the water bulk were examined at a molecular level using density profiles and radial pair distribution functions. Molecules containing only polar functional groups such as urea and glucose are found predominantly in the water bulk, forming an exclusion layer near the water surface. Substitution of a single polar group by an alkyl group in sugars and amides leads to the migration of the molecule toward the interface. Within the first 2 nm from the water surface, surface-active solutes lose their rotational freedom and adopt a preferred orientation with the alkyl group pointing toward the surface. The different packing within the interface leads to different solvation shell structures and enhanced interaction between the organic molecules and absorbed OH radicals. The simulations provide quantitative information about the dimension, composition, and organization of the air-water interface as well as about the nonreactive interaction of the OH radicals with the organic solutes. It suggests that increased concentrations, preferred orientations, and decreased solvation near the air-water surface may lead to differences in reactivities between surface-active and surface-inactive molecules. The results are important to explain how heterogeneous oxidation mechanisms and kinetics within interfaces may differ from those of the bulk.

Salt Enrichment and Dynamics in the Interface of Supercooled Aqueous Droplets

The interconversion reaction of NaCl between the contact-ion pair (CIP) and the solvent-separated ion pair (SSIP) as well as the free-ion state in cold droplets has not yet been investigated. We report direct computational evidence that the lower is the temperature, the closer to the surface the ion interconversion reaction takes place. In supercooled droplets the enrichment of the subsurface in salt becomes more evident. The stability of the SSIP relative to the CIP increases as the ion-pairing is transferred toward the droplet's outer layers. In the free-ion state, where the ions diffuse independently in the solution, the number density of Cl^- shows a broad maximum in the interior in addition to the well-known maximum in the surface. In the study of the reaction dynamics, we find a weak coupling between the interionic NaCl distance reaction coordinate and the solvent degrees of freedom, which contrasts with the diffusive crossing of the free energy barrier found in bulk solution modeling. The H₂O self-diffusion coefficient is found to be at least an order of magnitude larger than that in the bulk solution. We propose to exploit the enhanced surface ion concentration at low temperature to eliminate salts from droplets in native mass spectrometry ionization methods.

Conical Shape Fluctuations Determine the Rate of Ion Evaporation and the Emitted Cluster Size Distribution from Multicharged Droplets

The ion evaporation mechanism (IEM) is perceived to be a major pathway for disintegration of multi-ion charged droplets found in atmospheric and sprayed aerosols. However, the precise mechanism of IEM and the effect of the nature of the ions in the emitted cluster size distribution have not yet been established despite its broad use in mass spectrometry and atmospheric chemistry over the past half century. Here, we present a systematic study of the emitted ion cluster distribution in relation to their spatial distribution in the parent droplet using atomistic modeling. It is found that in the parent droplet, multiple kosmotropic and weakly polarizable chaotropic ions (Cs⁺) are buried deeper within the droplet than polarizable chaotropic ions (Cl^-, I^-). This differentiation in the ion location is only captured by a polarizable model. It is demonstrated that the emitted cluster size distribution is determined by dynamic conical deformations and not by the equilibrium ion depth within the parent droplet as the IEM models assume. Critical factors that determine the cluster size distribution such as the charge sign asymmetry that have not been considered in models and in experiments are presented. We argue that the existing IEM analytical models do not establish a clear difference between IEM and Rayleigh fission. We propose a shift in the existing view for IEM from the equilibrium properties of the parent droplet to the chemistry in the conical shape fluctuations that serve as the centers for ion emission. Consequently, chemistry in the conical fluctuations may also be a key element to explain charge states of macromolecules in mass spectrometry and may have potential applications in catalysis due to the electric field in the conical region.

Grotthuss Molecular Dynamics Simulations for Modeling Proton Hopping in Electrosprayed Water Droplets

Excess protons in water exhibit unique transport properties because they can rapidly hop along H-bonded water wires. Considerable progress has been made in unraveling this Grotthuss diffusion mechanism using quantum mechanical-based computational techniques. Unfortunately, high computational cost tends to restrict those techniques to small systems and short times. Molecular dynamics (MD) simulations can be applied to much larger systems and longer time windows. However, standard MD methods do not permit the dissociation/formation of covalent bonds, such that Grotthuss diffusion cannot be captured. Here, we bridge this gap by combining atomistic MD simulations (using Gromacs and TIP4P/2005 water) with proton hopping. Excess protons are modeled as hydronium ions that undergo H₃O⁺ + H₂O → H₂O + H₃O⁺ transitions. In accordance with ab initio MD data, these Grotthuss hopping events are executed in "bursts" with quasi-instantaneous hopping across one or more waters. The bursts are separated by regular MD periods during which H₃O⁺ ions undergo Brownian diffusion. The resulting proton diffusion coefficient agrees with the literature value. We apply this Grotthuss MD technique to highly charged water droplets that are in a size regime encountered during electrospray ionization (5 nm radius, ∼17,000 H₂O). The droplets undergo rapid solvent evaporation and occasional H₃O⁺ ejection, keeping them at ca. 81% of the Rayleigh limit. The simulated behavior is consistent with phase Doppler anemometry data. The Grotthuss MD technique developed here should be useful for modeling the behavior of various proton-containing systems that are too large for high-level computational approaches. In particular, we envision future applications related to electrospray processes, where earlier simulations used metal cations while in reality excess protons dominate.

Electron Ionization of Size-Selected Positively and Negatively Charged Helium Droplets

A beam of size-selected charged helium droplets was crossed with an electron beam, and the ion efficiency curves for the product droplets in all different charge states were recorded. We estimate that the selected helium droplets on their passage through the electron beam are hit by several hundred electrons which can interact with the individual He atoms of the droplets. Reaction channels corresponding to the removal or capture of up to eight electrons were identified, and in all cases, inelastic scattering and the formation of metastable helium played a significant role.

Low Density Interior in Supercooled Aqueous Nanodroplets Expels Ions to the Subsurface

The interaction between water and ions within droplets plays a key role in the chemical reactivity of atmospheric and man-made aerosols. Here we report direct computational evidence that in supercooled aqueous nanodroplets a lower density core of tetrahedrally coordinated water expels the cosmotropic ions to the denser and more disordered subsurface. In contrast, at room temperature, depending on the nature of the ion, the radial distribution in the droplet core is nearly uniform or elevated toward the center. We analyze the spatial distribution of a single ion in terms of a reference electrostatic model. The energy of the system in the analytical model is expressed as the sum of the electrostatic and surface energy of a deformable droplet. The model predicts that the ion is subject to a harmonic potential centered at the droplet's center of mass. We name this effect "electrostatic confinement". The model's predictions are consistent with the simulation findings for a single ion at room temperature but not at supercooling. We anticipate this study to be the starting point for investigating the structure of supercooled (electro)sprayed droplets that are used to preserve the conformations of macromolecules originating from the bulk solution.