Bridging electrostatic properties between nanoscopic and microscopic highly charged 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.

Real-Time Monitoring of the Evaporation and Fission of Electrospray-Ionized Polystyrene Beads and Bacterial Pellets at Elevated Temperatures

This study uses real-time monitoring, at microsecond time scales, with a charge-sensing particle detector to investigate the evaporation and fission processes of methanol/micrometer-sized polystyrene beads (PS beads) droplets and bacterial particles droplets generated via electrospray ionization (ESI) under elevated temperatures. By incrementally raising capillary temperatures, the solvent, such as methanol on 0.75 μm PS beads, experiences partial evaporation. Further temperature increase induces fission, and methanol molecules continue to evaporate until PS ions are detected after this range. Similar partial evaporation is observed on 3 μm PS beads. However, the shorter period of the fission temperature range is necessary compared to 0.75 μm PS beads. For the spherical-shaped bacterium, Staphylococcus aureus, the desolvation process shows a similar fission period as compared to 0.75 μm PS beads. Comparably, the rod-shaped bacteria, Escherichia coli EC11303, and E. coli strain W have shorter fission periods than S. aureus. This research provides insights into the evaporation and fission mechanisms of ESI droplets containing different sizes and shapes of micrometer-sized particles, contributing to a better understanding of gaseous macroion formation.

Spiers Memorial Lecture: Water at interfaces

In this article we discuss current issues in the context of the four chosen subtopics for the meeting: dynamics and nano-rheology of interfacial water, electrified/charged aqueous interfaces, ice interfaces, and soft matter/water interfaces. We emphasize current advances in both theory and experiment, as well as important practical manifestations and areas of unresolved controversy.

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.

Toward a molecular understanding of the surface composition of atmospherically relevant organic particles

Many mass spectrometry methods using various ionization sources provide bulk composition of airborne particles, but little is known about the surface species that play a major role in determining their physicochemical properties that impact air quality, climate, and health. The present work shows that the composition of surface layers of atmospherically relevant submicron organic particles can be probed without the use of an external ionization source. Solid dicarboxylic acid particles are used as models, with glutaric acid being the most efficient at generating ions. Coating with small diacids or products from α-pinene ozonolysis demonstrates that ions are ejected from the surface, providing surface molecular characterization of organic particles on the fly. This unique approach provides a path forward for elucidating the role of the surface in determining chemical and physical properties of particles, including heterogeneous reactions, particle growth, water uptake, and interactions with biological systems.

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.

Microdroplets can act as electrochemical cells

A water microdroplet in air or oil typically possesses an electric double layer (EDL) from the preferential adsorption of surface-bound ions at the periphery. We present the calculations of the ion gradients within a microdroplet at equilibrium, including systems containing buffers and water autoionization. These ion gradients are used to calculate the potential energy stored within the microdroplet. We consider how this stored potential energy can be utilized to drive chemical reactions, much like an electrochemical cell. Effective voltages as high as 111 mV are found for microdroplets having a low surface charge density (0.01 ions per nm2). Two sources of potential energy are investigated: (1) the electrostatic energy of the EDL of the microdroplet and (2) shifts in other chemical equilibria coupled to the main reaction through the EDL. A particularly important example of the latter is water autoionization, wherein the reaction of interest causes a flattening of the [H+] gradient within the EDL, resulting in a net recombination of H+ and OH- throughout the microdroplet. Numerical calculations are performed using a continuum model consisting of a balance between the electromigration and diffusion of ions throughout the microdroplet. Our treatment accounts for the autoionization of water and any chemical equilibrium of buffers present. The results are presented for uncharged water microdroplets with low amounts of salts and simple buffers in them. However, the calculational method presented here can be applied to microdroplets of any net charge, composed of any solvent, containing ions of any valence, and containing complex mixtures of chemical equilibria.

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.

Molecular Characterization of the Surface Excess Charge Layer in Droplets

The surface excess charge layer (SECL) in droplets has often been associated with distinct chemistry. We examine the effect of the nature of ions in the composition and structure of SECL by using molecular dynamics. We find that in the presence of simple ions the thickness of SECL is invariant not only with respect to droplet size but also with respect to the nature of the ions. In the presence of simple ions, this layer has a thickness of ∼1.5-1.7 nm but in the presence of macroions it may extend to ∼2.0 nm. The proportion of ions contained in SECL depends on the nature of the ions and the droplet size. For the same droplet size, I^- and model H₃O⁺ ions show considerably higher concentration than Na⁺ and Cl^- ions. We identify the maximum ion concentration region, which, in nanodrops, may partially overlap with SECL. As the relative shape fluctuations decrease when microdrop size is approached, the overlap between SECL and maximum ion concentration region increases. We suggest the extension of the bilayer droplet structure assumed in the equilibrium partitioning model of Enke to include the maximum ion concentration region that may not coincide with SECL in nanodrops. We compute the ion concentrations in SECL, which are those that should enter the kinetic equation in the ion-evaporation mechanism, instead of the overall drop ion concentration that has been used.

Restricted intramolecular rotation of fluorescent molecular rotors at the periphery of aqueous microdroplets in oil

Fluorescent molecular rotor dyes, including Cy3, Cy5, and Alexa Fluor 555, dissolved in micron-sized aqueous droplets (microdroplets) in oil were excited, and the fluorescence intensity was recorded as function of time. We observed lengthening of the fluorescence lifetime of these dyes at the water-oil periphery, which extended several microns inward. This behavior shows that intramolecular rotation is restricted at and near the microdroplet interface. Lengthened lifetimes were observed in water microdroplets but not in microdroplets composed of organic solvents. This lifetime change was relatively insensitive to added glycerol up to 60%, suggesting that solution viscosity is not the dominant mechanism. These restricted intramolecular rotations at and near the microdroplet periphery are consistent with the reduced entropy observed in chemical reactions in microdroplets compared to the same reaction conditions in bulk solution and helps us further understand why microdroplet chemistry differs so markedly from bulk-phase chemistry.

Strong Electric Field Observed at the Interface of Aqueous Microdroplets

Chemical reactions in aqueous microdroplets often exhibit unusual kinetic and thermodynamic properties not observed in bulk solution. While an electric field has been implicated at the water interface, there has been no direct measurement in aqueous microdroplets, largely due to the lack of proper measurement tools. Herein, we employ newly developed stimulated Raman excited fluorescence microscopy to measure the electric field at the water-oil interface of microdroplets. As determined by the vibrational Stark effect of a nitrile-bearing fluorescent probe, the strength of the electric field is found to be on the order of 10⁷ V/cm. This strong electric field aligns probe dipoles with respect to the interface. The formation of the electric field likely arises from charge separation caused by the adsorption of negative ions at the water-oil interface of microdroplets. We suggest that this strong electric field might account in part for the unique properties of chemical reactions reported in microdroplets.

Simple model for the electric field and spatial distribution of ions in a microdroplet

It is well established that the chemistry in microdroplets has been found to be radically different from reactions in bulk, particularly in the case of water. It has also been established that there is a threshold size for microdroplets to behave differently than droplets near the 10 µm diameter range. We present a three-dimensional electrostatic treatment in the spirit of the Gouy-Chapman model for double layers at interfaces. Our treatment predicts a strong concentration of charged molecules toward the surface of the droplet. As the droplet size deceases, the majority of the volume of the liquid experiences a large DC electric field. Such electric fields are highly unusual in a conducting fluid such as water. We believe that this unique environment helps to explain the reaction rate acceleration and new chemistry that have been observed in microdroplets compared to bulk phase.