The structure of ionic aqueous solutions at interfaces: An intrinsic structure analysis

Effect of an Electric Field on the Structure and Stability of Atmospheric Clusters

We study the influence of an applied electric field on the structure and stability of some common bimolecular clusters that are found in the atmosphere. These clusters play an important role in new particle formation (NPF). For low values of the electric field (i.e., |E| ≤ 0.01 V Å-1), we demonstrate that the field response of the clusters can be predicted from simply calculating the dipole moment of the cluster and the dipole moments of the constituent molecules and that the influence on the association energy of the cluster is minimal (i.e., <0.5 kcal mol-1). For higher field strengths |E| > 0.2 V Å-1, there can be more dramatic effects on both structure and energetics, as the induced dipole, charge transfer, and geometric distortion play a larger role. Although such large fields are not very relevant in the atmosphere, they do exist in some situations of experimental interest, such as near interfaces and in intense laser fields.

Does the Sign of Charge Affect the Surface Affinity of Simple Ions?

The role the charge sign of simple ions plays in determining their surface affinity in aqueous solutions is investigated by computer simulation methods. For this purpose, the free surface of aqueous solutions of fictitious salts is simulated at finite concentration both with nonpolarizable point-charge and polarizable Gaussian-charge potential models. The salts consist of monovalent cations and anions that are, apart from the sign of their charge, identical to each other. In particular, we consider the small Na⁺ and the large I^- ions together with their charge-inverted counterparts. In an attempt to avoid the interference even between the behavior of cations and anions, we also simulate systems containing only one of the above ions, and determine the free energy profile of these ions across the liquid-vapor interface of water at infinite dilution by potential of mean force (PMF) calculations. The obtained results reveal that, in the case of small ions, the anion is hydrated considerably stronger than the cation due to the close approach of water H atoms, bearing a positive fractional charge. As a consequence, the surface affinity of a small anion is even smaller than that of its cationic counterpart. However, considering that small ions are effectively repelled from the water surface, the importance of this difference is negligible. Further, a change in the hydration energy trends of the two oppositely charged ions is observed with their increasing size. This change is largely attributed to the fact that, with increasing ion size, the factor of 2 in the magnitude of the fractional charge of the closely approaching water atoms (i.e., O around cations and H around anions) outweighs the closer approach of the H than the O atom in the hydration energy. Thus, for large ions, being already surface active themselves, the surface affinity of the anion is larger than that of its positively charged counterpart. Further, such a difference is seen even in the case when the sign of the surface potential favors the adsorption of cations.

Surface Affinity of Tetramethylammonium Iodide in Aqueous Solutions: A Combined Experimental and Computer Simulation Study

The surface affinity of tetramethylammonium iodide (TMAI) in aqueous solutions is investigated by surface tension measurements and molecular dynamics computer simulations. Experiments, performed in the entire composition range of solubility using the pendant drop method with two different setups, clearly reveal that TMAI is a weakly capillary active salt. Computer simulations performed with the AMBER force field reproduce the experimental data very well, while two other major force fields (i.e., CHARMM and OPLS) can still reproduce the experimental trend qualitatively; however, even qualitative reproduction of the experimental trend requires scaling down the ion charges according to the Leontyev-Stuchebrukhov correction. On the other hand, the GROMOS force field fails in reproducing the experimentally confirmed capillary activity of TMAI. Molecular dynamics simulation results show that, among the two ions, iodide has a clearly larger surface affinity than tetramethylammonium (TMA⁺). Further, the adsorption of the I^- anions is strictly limited to the first molecular layer beneath the liquid-vapor interface, which is followed by several layers of their depletion. On the other hand, the net negative charge of the surface layer, caused by the excess amount of I^- with respect to TMA⁺, is compensated by a diffuse layer of adsorbed TMA⁺ cations, extending to or beyond the fourth molecular layer beneath the liquid surface.

Single-Particle Dynamics at the Intrinsic Surface of Aqueous Alkali Halide Solutions

The distribution of ions in the proximity of the liquid-vapor interface of their aqueous solution has been the subject of an intense debate during the last decade. The effects of ionic polarizability have been one of its salient aspects. Much less has been said about the corresponding dynamical properties, which are substantially unexplored. Here, we investigate the single-particle dynamics at the liquid-vapor interface of several alkali halide solutions, using molecular dynamics simulations with polarizable and nonpolarizable force fields and intrinsic surface analysis. We analyze the diffusion coefficient, residence time, and velocity autocorrelation function of water and ions and investigate how these properties depend on the molecular layer where they reside. While anions are found in the first molecular layer for relatively long times, cations are only making quick excursions into it, thanks to thermal fluctuations. The in-layer residence time of ions and their molar fraction in the layer turned out to be linearly dependent on each other. We interpret this unexpected result using a simple two-state model. In addition, we found that, unlike water and other neat molecular liquids that show a different diffusion mechanism at the surface than in the bulk of their liquid phase, ions do not enjoy enhanced mobility in the surface layer of their aqueous solution. This result indicates that ions in the surface layer are shielded by their nearest water neighbors from being exposed to the vapor phase as much as possible. Such positions are available for the ions at the negatively curved troughs of the molecularly rugged liquid surface.

Surface Affinity of Alkali and Halide Ions in Their Aqueous Solution: Insight from Intrinsic Density Analysis

The surface tension of all aqueous alkali halide solutions is higher than that of pure water. According to the Gibbs adsorption equation, this indicates a net depletion of these ions in the interfacial region. However, simulations and experiments show that large, soft ions, such as I^-, can accumulate at the liquid/vapor interface. The presence of a loose hydration shell is usually considered to be the reason for this behavior. In this work, we perform computer simulations to characterize the liquid-vapor interface of aqueous alkali chloride and sodium halide solutions systematically, considering all ions from Li⁺ to Cs⁺ and from F^- to I^-. Using computational methods for the removal of surface fluctuations, we analyze the structure of the interface at a dramatically enhanced resolution, showing that the positive excess originates in the very first molecular layer and that the next 3-4 layers account for the net negative excess. With the help of a fictitious system with charge-inverted ion pairs, we also show that it is not possible to rationalize the surface affinity of ions in solutions in terms of the properties of anions and cations separately. Moreover, the surface excess is generally dominated by the smaller of the two ions.

Experimental Depth Profiles of Surfactants, Ions, and Solvent at the Angstrom Scale: Studies of Cationic and Anionic Surfactants and Their Salting Out

Neutral impact ion scattering spectroscopy (NICISS) is used to measure the depth profiles of ionic surfactants, counterions, and solvent molecules on the angstrom scale. The chosen surfactants are 0.010 m tetrahexylammonium bromide (THA⁺/Br^-) and 0.0050 m sodium dodecyl sulfate (Na⁺/DS^-) in the absence and presence of 0.30 m NaBr in liquid glycerol. NICISS determines the depth profiles of the elements C, O, Na, S, and Br through the loss in energy of 5 keV He atoms that travel into and out of the liquid, which is then converted into depth. In the absence of NaBr, we find that THA⁺ and its Br^- counterion segregate together because of charge attraction, forming a narrow double layer that is 10 Å wide and 150 times more concentrated than in the bulk. With the addition of NaBr, THA⁺ is "salted out" to the surface, increasing the interfacial Br^- concentration by 3-fold and spreading the anions over a ∼30 Å depth. Added NaBr similarly increases the interfacial concentration of DS^- ions and broadens their positions. Conversely, the dissolved Br^- ions are significantly depleted over a depth of 0-40 Å from the surface because of charge repulsion from DS^- ions within the interfacial region. These different interfacial Br^- propensities correlate with previously measured gas-liquid reactivities: gaseous Cl₂ readily reacts with Br^- ions in the presence of THA⁺ but drops 70-fold in the presence of DS^-, demonstrating that surfactant headgroup charge controls the reactivity of Br^- through changes in its depth profile.

Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide on the Structure of the Aqueous Solution-Air Interface

We present the results of ab initio molecular dynamics simulations of the solution–air interface of aqueous lithium bromide (LiBr). We find that, in agreement with the experimental data and previous simulation results with empirical polarizable force field models, Br^– anions prefer to accumulate just below the first molecular water layer near the interface, whereas Li⁺ cations remain deeply buried several molecular layers from the interface, even at very high concentration. The separation of ions has a profound effect on the average orientation of water molecules in the vicinity of the interface. We also find that the hydration number of Li⁺ cations in the center of the slab N_{c,Li⁺–H2O} ≈ 4.7 ± 0.3, regardless of the salt concentration. This estimate is consistent with the recent experimental neutron scattering data, confirming that results from nonpolarizable empirical models, which consistently predict tetrahedral coordination of Li⁺ to four solvent molecules, are incorrect. Consequently, disruption of the hydrogen bond network caused by Li⁺ may be overestimated in nonpolarizable empirical models. Overall, our results suggest that empirical models, in particular nonpolarizable models, may not capture all of the properties of the solution–air interface necessary to fully understand the interfacial chemistry.

Debye screening, overscreening and specific adsorption in solutions of organic ions

Tetrabutylammonium (TBA) and tetraphenylborate (TPB) ions dissolved in dichloroethane (DCE) are widely used in electrochemistry of liquid-liquid interfaces. Unlike alkali halide solutions in water, TBA-TPB-DCE solutions feature large organic ions and a solvent with a dielectric constant almost one order of magnitude lower than that of water. This is expected to dramatically amplify the impact of ionic correlations in the properties of the solution. Here we report atomistic simulations of TBA-TPB-DCE solutions and analyze ion correlations, clustering, and charge screening effects. We target concentrations in the range of 0.01-0.25 molal (m), hence exploring concentration regimes typical for many experimental investigations. We show that the transition from monotonic to oscillatory decay of the charge density, which signals the onset of strong ion correlations, takes place in this concentration interval, leading to overscreening effects. Furthermore, we investigate the distribution and adsorption of ions at the DCE-air interface. Unlike what is observed for small inorganic ions in water at similar concentrations, we find that TPB and TBA ions accumulate near the DCE surface, resulting in significant interfacial clustering and adsorption at concentrations ∼0.25 m. TPB ions adsorb more strongly leading to free energy wells of ∼1-2 kBT. The adsorption modifies significantly the electrostatic potential of the DCE-air interface, which undergoes a shift of 0.2 V in going from pure DCE to TBA-TPB-DCE solutions at 0.25 m.

Specific cation effects at aqueous solution-vapor interfaces: Surfactant-like behavior of Li+ revealed by experiments and simulations

It is now well established by numerous experimental and computational studies that the adsorption propensities of inorganic anions conform to the Hofmeister series. The adsorption propensities of inorganic cations, such as the alkali metal cations, have received relatively little attention. Here we use a combination of liquid-jet X-ray photoelectron experiments and molecular dynamics simulations to investigate the behavior of K⁺ and Li⁺ ions near the interfaces of their aqueous solutions with halide ions. Both the experiments and the simulations show that Li⁺ adsorbs to the aqueous solution-vapor interface, while K⁺ does not. Thus, we provide experimental validation of the "surfactant-like" behavior of Li⁺ predicted by previous simulation studies. Furthermore, we use our simulations to trace the difference in the adsorption of K⁺ and Li⁺ ions to a difference in the resilience of their hydration shells.

Sensitivity of Electron Transfer Mediated Decay to Ion Pairing

Ion pairing in electrolyte solutions remains a topic of discussion despite a long history of research. Very recently, nearest-neighbor mediated electronic de-excitation processes of core hole vacancies (electron transfer mediated decay, ETMD) were proposed to carry a spectral fingerprint of local solvation structure and in particular of contact ion pairs. Here, for the first time, we apply electron-electron coincidence detection to a liquid microjet, and record ETMD spectra of Li 1s vacancies in aqueous solutions of lithium chloride (LiCl) in direct comparison to lithium acetate (LiOAc). A change in the ETMD spectrum dependent on the electrolyte anion identity is observed for 4.5 M salt concentration. We discuss these findings within the framework of the formation and presence of contact ion pairs and the unique sensitivity of ETMD spectroscopy to ion pairing.

Reversible ultralow-voltage liquid-liquid electrowetting without a dielectric layer

Electrowetting-on-dielectric devices typically have operating voltages of 10-20 V. A reduction in the operating voltage could greatly reduce the energy consumption of these devices. Herein, fully reversible one-electrolyte electrowetting of a droplet on a solid metal surface is reported for the first time. A reversible change of 29° for an 800 mV step is achieved. The effects of surface roughness, electrolyte composition, electrolyte concentration and droplet composition are investigated. It was found that there is a dramatic dependence of the reversibility and hysteresis of the system on these parameters, contrary to theoretical predictions. When a 3-chloro-1-propanol droplet is used, a system with no hysteresis and a 40° change in angle are obtained.

Single Particle Dynamics at the Intrinsic Surface of Various Apolar, Aprotic Dipolar, and Hydrogen Bonding Liquids As Seen from Computer Simulations

We investigate the single molecule dynamics at the intrinsic liquid/vapor interface of five different molecular liquids (carbon tetrachloride, acetone, acetonitrile, methanol, and water). After assessing that the characteristic residence times in the surface layer are long enough for a meaningful definition of several transport properties within the layer itself, we characterize the dynamics of the individual molecules at the liquid surface by analyzing their normal and lateral mean-square displacements and lateral velocity autocorrelation functions and, in the case of the hydrogen bonding liquids (i.e., water and methanol), also the properties of the hydrogen bonds. Further, dynamical properties as well as the clustering of the molecules residing unusually long in the surface layer are also investigated. The global picture emerging from this analysis is that of a noticeably enhanced dynamics of the molecules at the liquid surface, with diffusion coefficients up to 4 times larger than in the bulk, and the disappearance of the caging effect at the surface of all liquids but water. The dynamics of water is dominated by the strong hydrogen bonding structure also at the liquid surface.

Gas-Microjet Reactive Scattering: Collisions of HCl and DCl with Cool Salty Water

Liquid microjets provide a powerful means to investigate reactions of gases with salty water in vacuum while minimizing gas-vapor collisions. We use this technique to explore the fate of gaseous HCl and DCl molecules impinging on 8 molal LiCl and LiBr solutions at 238 K. The experiments reveal that HCl or DCl evaporate infrequently if they become thermally accommodated at the surface of either solution. In particular, we observe minimal thermal desorption of HCl following HCl collisions and no distinct evidence for rapid, interfacial DCl→HCl exchange following DCl collisions. These results imply that surface thermal motions are not generally strong enough to propel momentarily trapped HCl or DCl back into the gas phase before they ionize and disappear into solution. Instead, only HCl and DCl molecules that scatter directly from the surface escape entry. These recoiling molecules transfer less energy upon collision to LiBr/H2O than to LiCl/H2O, reflecting the heavier mass of Br(-) than of Cl(-) in the interfacial region.

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.

Thermodynamics of iodide adsorption at the instantaneous air-water interface

We performed molecular dynamics simulations using both polarizable and non-polarizable force fields to study the adsorption of iodide to the air-water interface. A novel aspect of our analysis is that the progress of ion adsorption is measured as the distance from the instantaneous interface, which is defined by a coarse-graining scheme proposed recently by Willard and Chandler ["Instantaneous liquid interfaces," J. Phys. Chem. B 114, 1954-1958 (2010)]. Referring structural and thermodynamic quantities to the instantaneous interface unmasks molecular-scale details that are obscured by thermal fluctuations when the same quantities are referred to an average measure of the position of the interface, such as the Gibbs dividing surface. Our results suggest that an ion adsorbed at the interface resides primarily in the topmost water layer, and the interfacial location of the ion is favored by enthalpy and opposed by entropy.

Alkali halide solutions under thermal gradients: soret coefficients and heat transfer mechanisms

We report an extensive analysis of the non-equilibrium response of alkali halide aqueous solutions (Na(+)/K(+)-Cl(-)) to thermal gradients using state of the art non-equilibrium molecular dynamics simulations and thermal diffusion forced Rayleigh scattering experiments. The coupling between the thermal gradient and the resulting ionic salt mass flux is quantified through the Soret coefficient. We find the Soret coefficient is of the order of 10(-3) K(-1) for a wide range of concentrations. These relatively simple solutions feature a very rich behavior. The Soret coefficient decreases with concentration at high temperatures (higher than T ∼ 315 K), whereas it increases at lower temperatures. In agreement with previous experiments, we find evidence for sign inversion in the Soret coefficient of NaCl and KCl solutions. We use an atomistic non-equilibrium molecular dynamics approach to compute the Soret coefficients in a wide range of conditions and to attain further microscopic insight on the heat transport mechanism and the behavior of the Soret coefficient in aqueous solutions. The models employed in this work reproduce the magnitude of the Soret coefficient, and the general dependence of this coefficient with temperature and salt concentration. We use the computer simulations as a microscopic approach to establish a correlation between the sign and magnitude of the Soret coefficients and ionic solvation and hydrogen bond structure of the solutions. Finally, we report an analysis of heat transport in ionic solution by quantifying the solution thermal conductivity as a function of concentration. The simulations accurately reproduce the decrease of the thermal conductivity with increasing salt concentration that is observed in experiments. An explanation of this behavior is provided.