The properties of ion-water clusters. II. Solvation structures of Na+, Cl−, and H+ clusters as a function of temperature

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.

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.

Quantum molecular simulations of micro-hydrated halogen anions

We report the results of a detailed and accurate investigation focused on structures and energetics of poly-hydrated halides employing first-principles polarizable halide-water potentials to describe the underlying forces. Following a bottom-up data-driven potential approach, we initially looked into the classical behavior of higher-order X^-(H₂O)_N clusters. We have located several low-lying energies, such as global and local minima, structures for each cluster, with various water molecules (up to N = 8) surrounding the halide anion (X^- = F^-, Cl^-, Br^-, I^-), employing an evolutionary programming method. It is found that the F^--water clusters exhibit different structural configurations than the heavier halides, however independently of the halide anion, all clusters show in general a selective growth with the anion preferring to be connected to the outer shell of the water molecule arrangements. In turn, path-integral molecular dynamics simulations are performed to incorporate explicitly nuclear quantum and thermal effects in describing the nature of halide ion microsolvation in such prototypical model systems. Our data reveal that at low finite temperatures, nuclear quantum effects affect certain structural properties, such as weakening hydrogen bonding between the halide anion and water molecules, with minor distortions in the water network beyond the first hydration shell, indicating local structure rearrangements. Such structural characteristics and the promising cluster size trends observed in the single-ion solvation energies motivated us to draw connections of small size cluster data to the limits of continuum bulk values, toward the investigation of the challenging computational modeling of bulk single ion hydration.

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.

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.

Ion Solvation and Transport in Narrow Carbon Nanotubes: Effects of Polarizability, Cation-π Interaction, and Confinement

Understanding ion solvation and transport under confinement is critical for a wide range of emerging technologies, including water desalination and energy storage. While molecular dynamics (MD) simulations have been widely used to study the behavior of confined ions, considerable deviations between simulation results depending on the specific treatment of intermolecular interactions remain. In the following, we present a systematic investigation of the structure and dynamics of two representative solutions, that is, KCl and LiCl, confined in narrow carbon nanotubes (CNTs) with a diameter of 1.1 and 1.5 nm, using a combination of first-principles and classical MD simulations. Our simulations show that the inclusion of both polarization and cation-π interactions is essential for the description of ion solvation under confinement, particularly for large ions with weak hydration energies. Beyond the variation in ion solvation, we find that cation-π interactions can significantly influence the transport properties of ions in CNTs, particularly for KCl, where our simulations point to a strong correlation between ion dehydration and diffusion. Our study highlights the complex interplay between nanoconfinement and specific intermolecular interactions that strongly control the solvation and transport properties of ions.

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.

Condensing water vapor to droplets generates hydrogen peroxide

It was previously shown [J. K. Lee et al., Proc. Natl. Acad. Sci. U.S.A, 116, 19294-19298 (2019)] that hydrogen peroxide (H2O2) is spontaneously produced in micrometer-sized water droplets (microdroplets), which are generated by atomizing bulk water using nebulization without the application of an external electric field. Here we report that H2O2 is spontaneously produced in water microdroplets formed by dropwise condensation of water vapor on low-temperature substrates. Because peroxide formation is induced by a strong electric field formed at the water-air interface of microdroplets, no catalysts or external electrical bias, as well as precursor chemicals, are necessary. Time-course observations of the H2O2 production in condensate microdroplets showed that H2O2 was generated from microdroplets with sizes typically less than ∼10 µm. The spontaneous production of H2O2 was commonly observed on various different substrates, including silicon, plastic, glass, and metal. Studies with substrates with different surface conditions showed that the nucleation and the growth processes of condensate water microdroplets govern H2O2 generation. We also found that the H2O2 production yield strongly depends on environmental conditions, including relative humidity and substrate temperature. These results show that the production of H2O2 occurs in water microdroplets formed by not only atomizing bulk water but also condensing water vapor, suggesting that spontaneous water oxidation to form H2O2 from water microdroplets is a general phenomenon. These findings provide innovative opportunities for green chemistry at heterogeneous interfaces, self-cleaning of surfaces, and safe and effective disinfection. They also may have important implications for prebiotic chemistry.

pH-Dependent Water Clusters in Photoacid Solution: Real-Time Observation by ToF-SIMS at a Submicropore Confined Liquid-Vacuum Interface

Water clusters are ubiquitously formed in aqueous solutions by hydrogen bonding, which is quite sensitive to various environment factors such as temperature, pressure, electrolytes, and pH. Investigation of how the environment has impact on water structure is important for further understanding of the nature of water and the interactions between water and solutes. In this work, pH-dependent water structure changes were studied by monitoring the changes for the size distribution of protonated water clusters by in-situ liquid ToF-SIMS. In combination with a light illumination system, in-situ liquid ToF-SIMS was used to real-time measure the changes of a light-activated organic photoacid under different light illumination conditions. Thus, the proton transfer and pH-mediated water cluster changes were analyzed in real-time. It was found that higher concentration of free protons could lead to a strengthened local hydrogen bonding network as well as relatively larger protonated water clusters in both organic acid and inorganic acid. Besides, the accumulation of protons at the liquid-vacuum interface under light illumination was observed owing to the affinity of organic molecules to the low-pressure gas phase. The application of in-situ liquid ToF-SIMS analysis in combination with in-situ light illumination system opened up an avenue to real-time investigate light-activated reactions. Besides, the results regarding water structure changes in acidic solutions showed important insights in related atmospheric and physiochemical processes.

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.

Effects of Temperature on Cs+(H2O)20 Clathrate Structure

Clusters consisting of 20 water molecules and a single cesium ion are especially stable due to their clathrate structure that is composed exclusively of three-coordinate water molecules. Clathrate stability was investigated using infrared photodissociation (IRPD) spectroscopy in the free-OH stretching region (∼3600-3800 cm^-1) at ion cell temperatures between 135 and 355 K. At 275 K and colder, IRPD spectra of Cs⁺(H₂O)₂₀ have just one acceptor-acceptor-donor band. At higher temperatures, a higher-energy acceptor-donor band emerges and grows in intensity. Non-clathrate Na⁺(H₂O)₂₀ structures contain both of these bands, which do not change significantly in intensity over the temperature range. These results indicate a rapid onset in the conversion from clathrate to non-clathrate structures with temperature and suggest that some clathrate population remains even at the highest temperatures investigated. These results provide new insights into the role of entropy in clathrate stability.

Propensity of Hydrated Excess Protons and Hydroxide Anions for the Air-Water Interface

Significant effort has been undertaken to better understand the molecular details governing the propensity of ions for the air-water interface. Facilitated by computationally efficient reactive molecular dynamics simulations, new and statistically conclusive molecular-scale results on the affinity of the hydrated excess proton and hydroxide anion for the air-water interface are presented. These simulations capture the dynamic bond breaking and formation processes (charge defect delocalization) that are important for correctly describing the solvation and transport of these complex species. The excess proton is found to be attracted to the interface, which is correlated with a favorable enthalpic contribution and consistent with reducing the disruption in the hydrogen bond network caused by the ion complex. However, a recent refinement of the underlying reactive potential energy function for the hydrated excess proton shows the interfacial attraction to be weaker, albeit nonzero, a result that is consistent with the experimental surface tension measurements. The influence of a weak hydrogen bond donated from water to the protonated oxygen, recently found to play an important role in excess hydrated proton transport in bulk water, is seen to also be important for this study. In contrast, the hydroxide ion is found to be repelled from the air-water interface. This repulsion is characterized by a reduction of the energetically favorable ion-water interactions, which creates an enthalpic penalty as the ion approaches the interface. Finally, we find that the fluctuation in the coordination number around water sheds new light on the observed entropic trends for both ions.

Molecular dynamics simulation of the cooperative effect by different force fields in monosodium glutamate aqueous solution

Different force fields (GROMOS and OPLS) in conjunction with different water (SPC, SPCE, TIP3P, TIP4P and TIP5P) were assessed using molecular dynamics simulations of monosodium glutamate (MSG) aqueous solution.

Properties of hydrated excess protons near phospholipid bilayers

The behavior of the hydrated excess proton near different lipid membranes is studied with the third generation of the multistate empirical valence bond (MS-EVB3) model [Wu, Y. J.; Chen, H. N.; Wang, F.; Paesani, F.; Voth, G. A. J. Phys. Chem. B 2008, 112, 467]. Dioleoylphosphatidylcholine (DOPC), dioleoylphosphatidylethanolamine (DOPE), and dioleoylphosphatidylglycerol (DOPG) are selected as example lipids. In spite of the differences of the head groups, the molecular dynamics simulations show that all the lipid membranes have a proton-collecting antenna effect with no free energy barrier between the bulk water and interface regions. By comparison with classical hydronium model simulations, it is found that an appropriate description of proton Grotthuss shuttling and associated charge defect delocalization are necessary to obtain the correct free energy profile for the hydrated excess proton. In addition, nanosecond time scale sampling is essential to evaluate the free energy profiles, because certain slow motions are needed to stabilize the excess proton in the deep membrane interface region. It is also found that the lateral diffusion coefficients are 1 order of magnitude smaller in the interface region than in bulk water for all the lipids. These coefficients are almost the same as those of the lipid head groups. Finally, since the lipid phosphates may possibly be protonated due to the proton antenna effect of the membrane, phosphate group protonation is investigated and discussed within the MS-EVB framework.

Micro-solvation of the Zn2+ ion-a case study

The structures and thermodynamic parameters of hydrated zinc ion clusters incorporating a single zinc ion and up to eighteen water molecules have been determined with a quantum mechanical hybrid density functional, namely B3LYP using cc-PVDZ basis functions for H and O and a split valence 6-31G (d, p) basis function for Zn. The geometries for all the zinc ion water clusters are optimized with several initial guess structures and without imposing any initial symmetry restriction. Zinc metal ion is found to be preferably four coordinated for smaller sizes of hydrated cluster but attains an octahedral coordination for larger sizes (n >or= 11) of hydrated cluster. Structures with seven or more than seven water molecules attached directly to the central zinc ion are not found. The calculated gas phase coordination number in the first solvation sphere of a large hydrated zinc metal ion is found to be six and the same is also confirmed in the force field based classical molecular dynamics simulation study for an aqueous zinc ion and thus confirms the experimental findings. The equilibrium zinc-oxygen distance of 2.09328 A at the present B3LYP level of study is in excellent agreement with the X-ray diffraction result of 2.093 +/- 0.002 A for a hexahydrated zinc cluster.

Calorimetric observation of the melting of free water nanoparticles at cryogenic temperatures

We present an experimental study of the thermodynamics of free, size-selected water cluster anions consisting of 48 and 118 molecules. The measured caloric curves of the clusters are bulklike at low temperatures but show a well-defined, particle-size specific transition at 93+/-3 K for (H2O)48- and 118+/-3 K for (H2O)118-. At the transition temperature the heat capacity strongly increases, which marks the onset of melting.

Spiers Memorial Lecture. Ions at aqueous interfaces

Studies of aqueous interfaces and of the behavior of ions therein have been profiting from a recent remarkable progress in surface selective spectroscopies, as well as from developments in molecular simulations. Here, we summarize and place in context our investigations of ions at aqueous interfaces employing molecular dynamics simulations and electronic structure methods, performed in close contact with experiment. For the simplest of these interfaces, i.e. the open water surface, we demonstrate that the traditional picture of an ion-free surface is not valid for large, soft (polarizable) ions such as the heavier halides. Both simulations and spectroscopic measurements indicate that these ions can be present and even enhanced at surface of water. In addition we show that the ionic product of water exhibits a peculiar surface behavior with hydronium but not hydroxide accumulating at the air/water and alkane/water interfaces. This result is supported by surface-selective spectroscopic experiments and surface tension measurements. However, it contradicts the interpretation of electrophoretic and titration experiments in terms of strong surface adsorption of hydroxide; an issue which is further discussed here. The applicability of the observed behavior of ions at the water surface to investigations of their affinity for the interface between proteins and aqueous solutions is explored. Simulations show that for alkali cations the dominant mechanism of specific interactions with the surface of hydrated proteins is via ion pairing with negatively charged amino acid residues and with the backbone amide groups. As far as halide anions are concerned, the lighter ones tend to pair with positively charged amino acid residues, while heavier halides exhibit affinity to the amide group and to non-polar protein patches, the latter resembling their behavior at the air/water interface. These findings, together with results for more complex molecular ions, allow us to formulate a local model of interactions of ions with proteins with the aim to rationalize at the molecular level ion-specific Hofmeister effects, e.g. the salting out of proteins.

Acid-base analysis: a critique of the Stewart and bicarbonate-centered approaches

When approaching the analysis of disorders of acid-base balance, physical chemists, physiologists, and clinicians, tend to focus on different aspects of the relevant phenomenology. The physical chemist focuses on a quantitative understanding of proton hydration and aqueous proton transfer reactions that alter the acidity of a given solution. The physiologist focuses on molecular, cellular, and whole organ transport processes that modulate the acidity of a given body fluid compartment. The clinician emphasizes the diagnosis, clinical causes, and most appropriate treatment of acid-base disturbances. Historically, two different conceptual frameworks have evolved among clinicians and physiologists for interpreting acid-base phenomena. The traditional or bicarbonate-centered framework relies quantitatively on the Henderson-Hasselbalch equation, whereas the Stewart or strong ion approach utilizes either the original Stewart equation or its simplified version derived by Constable. In this review, the concepts underlying the bicarbonate-centered and Stewart formulations are analyzed in detail, emphasizing the differences in how each approach characterizes acid-base phenomenology at the molecular level, tissue level, and in the clinical realm. A quantitative comparison of the equations that are currently used in the literature to calculate H+concentration ([H+]) is included to clear up some of the misconceptions that currently exist in this area. Our analysis demonstrates that while the principle of electroneutrality plays a central role in the strong ion formulation, electroneutrality mechanistically does not dictate a specific [H+], and the strong ion and bicarbonate-centered approaches are quantitatively identical even in the presence of nonbicarbonate buffers. Finally, our analysis indicates that the bicarbonate-centered approach utilizing the Henderson-Hasselbalch equation is a mechanistic formulation that reflects the underlying acid-base phenomenology.

Molecular dynamics simulation of LiTFSI-acetamide electrolytes: structural properties

The liquid structures of nonaqueous electrolytes composed of lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) and acetamide, with LiTFSI/acetamide molar ratios of 1:2, 1:4, and 1:6, were studied by molecular dynamics simulations. The simulations indicate that the Li+ cations prefer to be six-coordinate by the sulfonyl oxygen atoms of the TFSI- anions and the carbonyl oxygen atoms of the acetamide molecules, rather than by the most electronegative nitrogen atom of the TFSI- anion. Therefore, close Li+-TFSI- contact pairs exist in the system. The TFSI- anion prefers to provide only one of four possible oxygen atoms to coordinate to the same Li+ cation. Three conformations (cis, trans, and gauche) of the TFSI- anions were found to coexist in the liquid electrolyte. At high salt concentrations, the TFSI- anions mainly adopt the gauche conformation in order to provide more oxygen atoms to coordinate to different Li+ cations, while simultaneously reducing the repulsion among the Li+ cations. On the other hand, the fraction of TFSI- anions adopting the cis conformation is largest for the system with the molar ratio of 1:6, in which many clusters, mainly composed of the Li+ cations and the TFSI- anions, are immersed in the acetamide molecules. The size and charge distribution of clusters were also investigated. In the system with the molar ratio of 1:2, nearly all of the ions in the PBC (periodic boundary conditions) box aggregate into a bulky cluster that gradually disassembles into small clusters with decreasing salt concentration. The addition of acetamide molecules was found to effectively relax the liquid electrolyte structure, and the system with the molar ratio of 1:4 was found to exhibit a more homogeneous liquid structure than the other two electrolyte systems with molar ratios of 1:2 and 1:6.