Microsolvation of the sodium and iodide ions and their ion pair in acetonitrile clusters: a theoretical study

The structural and thermodynamic properties of Na+(CH3CN)n, I-(CH3CN)n, and NaI(CH3CN)n clusters have been investigated by means of room-temperature Monte Carlo simulations with model potentials developed to reproduce the properties of small clusters predicted by quantum chemistry. Ions are found to adopt an interior solvation shell structure, with a first solvation shell containing approximately 6 and approximately 8 acetonitrile molecules for large Na+(CH3CN)n and I-(CH3CN)n clusters, respectively. Structural features of Na+(CH3CN)n are found to be similar to those of Na+(H2O)n clusters, but those of I-(CH3CN)n contrast with those of I-(H2O)n, for which "surface" solvation structures were observed. The potential of mean force calculations demonstrates that the NaI ion pair is thermodynamically stable with respect to ground-state ionic dissociation in acetonitrile clusters. The properties of NaI(CH3CN)n clusters exhibit some similarities with NaI(H2O)n clusters, with the existence of contact ion pair and solvent-separated ion pair structures, but, in contrast to water clusters, both types of ion pairs adopt a well-defined interior ionic solvation shell structure in acetonitrile clusters. Whereas contact ion pair species are thermodynamically favored in small clusters, solvent-separated ion pairs tend to become thermodynamically more stable above a cluster size of approximately 26. Hence, ground-state charge separation appears to occur at larger cluster sizes for acetonitrile clusters than for water clusters. We propose that the lack of a large Na+(CH3CN)n product signal in NaI(CH3CN)n multiphoton ionization experiments could arise from extensive stabilization of the ground ionic state by the solvent and possible inhibition of the photoexcitation mechanism, which may be less pronounced for NaI(H2O)n clusters because of surface solvation structures. Alternatively, increased solvent evaporation resulting from larger excess energies upon photoexcitation or major solvent reorganization on the ionized state could account for the observed solvent-selectivity in NaI cluster multiphoton ionization.     

Water Structure and Phase Transition Near a Surface

This article concerns the density and orientational structure of water near uncharged (hydrophobic) and charged (hydrophilic) surfaces and the surface-induced phase transition accompanying a powerful, attractive surface force. Previously reported studies on the water structure using the integral equation theories and computer simulations are critically reviewed, and some controversial aspects are reexamined by additional calculations. Effects of the surface curvature and surface charge density on the structure are analyzed in detail. A new analysis is performed on the phase transition for water containing a trace amount of a hydrophobic component at surfaces. It is shown that the transition occurs even at charged surfaces and the attraction induced between surfaces predominates over the strong Coulombic repulsion at sufficiently small separations. At highly charged surfaces, however, the transition does not occur.     

Surfaces affect ion pairing

In water, positive ions attract negative ions. That attraction can be modulated if a hydrophobic surface is present near the two ions in water. Using computer simulations with explicit and implicit water, we study how an ion embedded on a hydrophobic surface interacts with another nearby ion in water. Using hydrophobic surfaces with different curvatures, we find that the contact interaction between a positive and negative ion is strongly affected by the curvature of an adjacent surface, either stabilizing or destabilizing the ion pair. We also find that the solvent-separated ion pair (SSIP) can be made more stable than the contacting ion pair by the presence of a surface. This may account for why bridging waters are often found in protein crystal structures. We also note that implicit solvent models do not account for SSIPs. Finally, we find that there are charge asymmetries: an embedded positive charge attracting a negative ion is different than an embedded negative charge attracting a positive ion. Such asymmetries are also not predicted by implicit solvent models. These results may be useful for improving computational models of solvation in biology and chemistry.     

Interactions between nonpolar surfaces coated with the nonionic surfactant hexaoxyethylene dodecyl ether C12E6 and the origin of surface charges at the air/water interface

The interactions between nonpolar surfaces coated with the nonionic surfactant hexaoxyethylene dodecyl ether C12E6 were investigated using two techniques and three different types of surfaces. As nonpolar surfaces, the air/water interface, silanated negatively charged glass, and thiolated uncharged gold surfaces were chosen. The interactions between the air/water interfaces were measured with a thin film pressure balance in terms of disjoining pressure as a function of film thickness. The interactions between the solid/liquid interfaces were determined using a bimorph surface force apparatus. The influence of the nature of the surface on the interaction forces was investigated at surfactant concentrations below and above the cmc. The adsorption of the nonionic surfactant on the uncharged thiolated surface does not, as expected, lead to any buildup of a surface charge. On the other hand, adsorption of C12E6 on the charged silanated glass and the charged air/water interface results in a lowering of the surface charge density. The reduction of the surface charge density on the silanated glass surfaces is rationalized by changes in the dielectric permittivity around the charged silanol groups. The reason for the surface charge observed at the air/water interface as well as its decrease with increasing surfactant concentration is discussed and a new mechanism for generation of OH- ions at this particular interface is proposed.     

Electro-osmosis at inhomogeneous charged surfaces: hydrodynamic versus electric friction

Electrokinetic methods are efficient in probing the electrostatic surface properties of charged systems. However, anomalies observed in experiments indicate that the classical electrokinetic theory should be reconsidered. Using Green's function methods and hydrodynamic simulations, we investigate electro-osmosis driven by electric-field-induced ion motion near a charged planar substrate with smooth or rough boundary. First, a reformulation of electro-osmotic theory for planar charged surfaces employing Green's functions shows that the Helmholtz-Smoluchowski (HS) relation between electrostatic potential and solvent velocity is exact for smooth surfaces, even in the presence of ion correlations. Deviations from HS theory are caused by combined hydrodynamic and electric surface friction, as our hydrodynamic simulations of ions at smooth and corrugated charged surfaces in lateral electric fields demonstrate. Within the simulations, hydrodynamic interactions are treated in the continuum limit and the presence of a no-slip boundary condition at the surface is taken into account. While electrofriction is relevant in highly charged system and/or for multivalent ions, hydrodynamic friction is dominant in systems with moderate surface charge density and/or low ionic valency. We also derive the effective electrokinetic surface charge from the electro-osmotic solvent profiles, which is substantially reduced when compared with the bare value and shows qualitative agreement with the experimental tendency.     

How to distinguish energetic surface heterogeneity from electrostatic interactions in the case of hydrogen ion adsorption from solution onto oxides

It was shown that the adsorption of uncharged particles onto energetically heterogeneous surfaces with Gaussian-like distribution of adsorption energy can be described by the same adsorption isotherm as the adsorption of hydrogen ions from solution onto homogeneous oxide surfaces with surface potentials given by a quasi-Nernst formula. It exemplifies that the separation of electrostatic from energetic factors in the case of ion adsorption onto heterogeneous oxide surfaces is very difficult. Nevertheless, application of calorimetric data for proton adsorption makes it possible to estimate both factors.     

Molecular mechanism of transporting a polarizable iodide anion across the water-CCl4 liquid/liquid interface

The result of transferring a polarizable iodide anion across the H2O-CCl4 liquid/liquid interface was investigated in this study. The computed transfer-free energy profile or potential of mean force exhibits a minimum near the Gibbs dividing surface. These system characteristics are similar to those found in a corresponding study of iodide transfer across the H2O-vapor interface; however, the free energy minimum was lower at the H2O-vapor interface. Molecular dynamics simulations were also carried out to compare the concentrations of NaCl, NaBr, and NaI at the H2O-vapor and H2O-CCl4 interfaces. While the concentration of bromide and iodide ions were lower at the H2O-CCl4 interface when compared to the H2O-vapor interface, the chloride ion concentrations were similar at both interfaces. Analysis of the solvation structures of iodide and chloride ions revealed that the more polarizable iodide ion was less solvated than the chloride ion at the interface. This characteristic brought the iodide ion into greater contact with CCl4, resulting in repulsive interactions with CCl4 and reducing its tendency to move to the interface.     

Electrode screening by ionic liquids

In this work we are concerned with the short-range screening provided by the ionic liquid dimethylimidazolium chloride near a charged wall. We study the free energy profiles (or potentials of mean force) for charged and neutral solutes as a function of distance from a charged wall. Four different wall charge densities are used in addition to a wall with zero charge. The highest magnitude of the charge densities is ±1 e nm(-2) which is close to the maximum limit of charge densities accessible in experiments, while the intermediate charges ±0.5 e nm(-2) are in the range of densities typically used in most of the experimental studies. Positively and negatively charged solutes of approximately the size of a BF ion and a Cl(-) ion are used as probes. We find that the ionic liquid provides excellent electrostatic screening at a distance of 1-2 nm. The free energy profiles show minima which are due to layering in the ionic liquid near the electrodes. This indicates that the solute ions tend to displace ionic liquid ions in the layers when approaching the electrode. The important role of non-electrostatic forces is demonstrated by the oscillations in the free energy profiles of uncharged solutes as a function of distance from the wall.     

Combined depletion and electrostatic forces in polymer-induced membrane adhesion: a theoretical model

We develop a semi-quantitative analytical theory to describe adhesion between two identical planar charged surfaces embedded in a polymer-containing electrolyte solution. Polymer chains are uncharged and differ from the solvent by their lower dielectric permittivity. The solution mimics physiological fluids: It contains 0.1 M of monovalent ions and a small number of divalent cations that form tight bonds with the headgroups of charged lipids. The components have heterogeneous spatial distributions. The model was derived self-consistently by combining: (a) a Poisson-Boltzmann like equation for the charge densities, (b) a continuum mean-field theory for the polymer profile, (c) a solvation energy forcing the ions toward the polymer-poor regions, and (d) surface interactions of polymers and electrolytes. We validated the theory via extensive coarse-grained Molecular Dynamics (MD) simulations. The results confirm our analytical model and reveal interesting details not detected by the theory. At high surface charges, polymer chains are mainly excluded from the gap region, while the concentration of ions increases. The model shows a strong coupling between osmotic forces, surface potential and salting-out effects of the slightly polar polymer chains. It highlights some of the key differences in the behaviour of monomeric and polymeric mixed solvents and their responses to Coulomb interactions. Our main findings are: (a) the onset of long-ranged ion-induced polymer depletion force that increases with surface charge density and (b) a polymer-modified repulsive Coulomb force that increases with surface charge density. Overall, the system exhibits homeostatic behaviour, resulting in robustness against variations in the amount of charges. Applications and extensions of the model are briefly discussed.     

A method for treating the passage of a charged hard sphere ion as it passes through a sharp dielectric boundary

In the implicit solvent models of electrolytes (such as the primitive model (PM)), the ions are modeled as point charges in the centers of spheres (hard spheres in the case of the PM). The surfaces of the spheres are not polarizable which makes these models appropriate to use in computer simulations of electrolyte systems where these ions do not leave their host dielectrics. The same assumption makes them inappropriate in simulations where these ions cross dielectric boundaries because the interaction energy of the point charge with the polarization charge induced on the dielectric boundary diverges. In this paper, we propose a procedure to treat the passage of such ions through dielectric interfaces with an interpolation method. Inspired by the "bubble ion" model (in which the ion's surface is polarizable), we define a space-dependent effective dielectric coefficient, ε(eff)(r), for the ion that overlaps with the dielectric boundary. Then, we replace the "bubble ion" with a point charge that has an effective charge q/ε(eff)(r) and remove the portion of the dielectric boundary where the ion overlaps with it. We implement the interpolation procedure using the induced charge computation method [D. Boda, D. Gillespie, W. Nonner, D. Henderson, and B. Eisenberg, Phys. Rev. E 69, 046702 (2004)]. We analyze the various energy terms using a spherical ion passing through an infinite flat dielectric boundary as an example.     

Collisions of DCl with pure and salty glycerol: enhancement of interfacial D --> H exchange by dissolved NaI

The fate of DCl molecules striking pure glycerol and a 2.6 M NaI-glycerol solution is investigated using scattering, uptake, and residence time measurements. We find that dissolved Na+ and I- ions alter every gas-liquid pathway from the moment of contact of DCl with the surface to its eventual emergence as HCl. In particular, the salt enhances both trapping-desorption of DCl and interfacial DCl --> HCl exchange at the expense of DCl entry into the bulk solution. The reduced entry and enhanced desorption of thermalized DCl molecules are interpreted by assuming that Na+ and I- ions bind to interfacial OH groups and tie up surface sites that would otherwise capture incoming DCl molecules. These ion-glycerol interactions may also be responsible for enhancing interfacial D --> H exchange by disrupting the interfacial hydrogen bond network that carries the newly formed H+ ion away from its Cl- pair. This disruption may increase the fraction of interfacial Cl- and H+ that recombine and desorb immediately as HCl before the ions separate and diffuse deeply into the bulk.     

Electrostatic properties of aqueous salt solution interfaces: a comparison of polarizable and nonpolarizable ion models

The effects of ion force field polarizability on the interfacial electrostatic properties of approximately 1 M aqueous solutions of NaCl, CsCl, and NaI are investigated using molecular dynamics simulations employing both nonpolarizable and Drude-polarizable ion sets. Differences in computed depth-dependent orientational distributions, "permanent" and induced dipole and quadrupole moment profiles, and interfacial potentials are obtained for both ion sets to further elucidate how ion polarizability affects interfacial electrostatic properties among the various salts relative to pure water. We observe that the orientations and induced dipoles of water molecules are more strongly perturbed in the presence of polarizable ions via a stronger ionic double layer effect arising from greater charge separation. Both anions and cations exhibit enhanced induced dipole moments and strong z alignment in the vicinity of the Gibbs dividing surface (GDS) with the magnitude of the anion induced dipoles being nearly an order of magnitude larger than those of the cations and directed into the vapor phase. Depth-dependent profiles for the trace and z z components of the water molecular quadrupole moment tensors reveal 40% larger quadrupole moments in the bulk phase relative to the vapor which mimics a similar observed 40% increase in the average water dipole moment. Across the GDS, the water molecular quadrupole moments increase nonmonotonically (in contrast to the water dipoles) and exhibit a locally reduced contribution just below the surface due to both orientational and polarization effects. Computed interfacial potentials for the nonpolarizable salts yield values 20-60 mV more positive than pure water and increase by an additional 30-100 mV when ion polarizability is included. A rigorous decomposition of the total interfacial potential into ion monopole, water and ion dipole, and water quadrupole components reveals that a very strong, positive ion monopole contribution is offset by negative contributions from all other potential sources. Water quadrupole components modulated by the water density contribute significantly to the observed interfacial potential increments and almost entirely explain observed differences in the interfacial potentials for the two chloride salts. By lumping all remaining nonquadrupole interfacial potential contributions into a single "effective" dipole potential, we observe that the ratio of quadrupole to "effective" dipole contributions range from 2:1 in CsCl to 1:1.5 in NaI, suggesting that both contributions are comparably important in determining the interfacial potential increments. We also find that oscillations in the quadrupole potential in the double layer region are opposite in sign and partially cancel those of the "effective" dipole potential.     

Mass selective axial ion ejection from a linear quadrupole ion trap

The electric fields responsible for mass-selective axial ejection (MSAE) of ions trapped in a linear quadrupole ion trap have been studied using a combination of analytic theory and computer modeling. Axial ejection occurs as a consequence of the trapped ions' radial motion, which is characterized by extrema that are phase-synchronous with the local RF potential. As a result, the net axial electric field experienced by ions in the fringe region, over one RF cycle, is positive. This axial field depends strongly on both the axial and radial ion coordinates. The superposition of a repulsive potential applied to an exit lens with the diminishing quadrupole potential in the fringing region near the end of a quadrupole rod array can give rise to an approximately conical surface on which the net axial force experienced by an ion, averaged over one RF cycle, is zero. This conical surface has been named the cone of reflection because it divides the regions of ion reflection and ion ejection. Once an ion penetrates this surface, it feels a strong net positive axial force and is accelerated toward the exit lens. As a consequence of the strong dependence of the axial field on radial displacement, trapped thermalized ions can be ejected axially from a linear ion trap in a mass-selective way when their radial amplitude is increased through a resonant response to an auxiliary signal.     

Specific ion adsorption and surface forces in colloid science

Mean-field theories that include nonelectrostatic interactions acting on ions near interfaces have been found to accommodate many experimentally observed ion specific effects. However, it is clear that this approach does not fully account for the liquid molecular structure and hydration effects. This is now improved by using parametrized ionic potentials deduced from recent nonprimitive model molecular dynamics (MD) simulations in a generalized Poisson-Boltzmann equation. We investigate how ion distributions and double layer forces depend on the choice of background salt. There is a strong ion specific double layer force set up due to unequal ion specific short-range potentials acting between ions and surfaces.     

Dispersion self-free energies and interaction free energies of finite-sized ions in salt solutions

The role for many-body dipolar (dispersion) potentials in ion-solvent and ion-solvent-interface interactions is explored. Such many-body potentials, accessible in principle from measured dielectric data, are necessary in accounting for Hofmeister specific ion effects. Dispersion self-energy is the quantum electrodynamic analogue of the Born electrostatic self-energy of an ion. We here describe calculations of dispersion self-free energies of four different anions (OH-, Cl-, Br-, and I-) that take finite ion size into account. Three different examples of self-free energy calculations are presented. These are the self-free energy of transfer of an ion to bulk solution, which influences solubility; the dispersion potential acting between one ion and an air-water interface (important for surface tension calculations); and the dispersion potential acting between two ions (relevant to activity coefficient calculations). To illustrate the importance of dispersion self-free energies, we compare the Born and dispersion contributions to the free energy of ion transfer from water to air (oil). We have also calculated the change in interfacial tension with added salt for air (oil)-water interfaces. A new model is used that includes dispersion potentials acting on the ions near the interface, image potentials, and ions of finite size that are allowed to spill over the solution-air interface. It is shown that interfacial free energies require a knowledge of solvent profiles at the interface.     

Rectification of nanopores at surfaces

At the nanoscale, methods to measure surface charge can prove challenging. Herein we describe a general method to report surface charge through the measurement of ion current rectification of a nanopipette brought in close proximity to a charged substrate. This method is able to discriminate between charged cationic and anionic substrates when the nanopipette is brought within distances from ten to hundreds of nanometers from the surface. Further studies of the pH dependence on the observed rectification support a surface-induced mechanism and demonstrate the ability to further discriminate between cationic and nominally uncharged surfaces. This method could find application in measurement and mapping of heterogeneous surface charges and is particularly attractive for future biological measurements, where noninvasive, noncontact probing of surface charge will prove valuable.     

Electrokinetic actuation of an uncharged polarizable dielectric droplet in charged hydrogel medium

Electrokinetic transport of an uncharged nonconducting microsized liquid droplet in a charged hydrogel medium is studied. Dielectric polarization of the liquid drop under the action of an externally imposed electric field induces a non-homogeneous charge density at the droplet surface. The interactions of the induced surface charge of the droplet with the immobile charges of the hydrogel medium generates an electric force to the droplet, which actuates the drop through the charged hydrogel medium. A numerical study based on the first principle of electrokinetics is adopted. Dependence of the droplet velocity on its dielectric permittivity, bulk ionic concentration, and immobile charge density of the gel is analyzed. The surface conduction is significant in presence of charged gel, which creates a concentration polarization. The impact of the counterion saturation in the Debye layer due to the dielectric decrement of the medium is addressed. The modified Nernst–Planck equation for ion transport and the Poisson equation for the electric field is considered to take into account the dielectric polarization. A quadrupolar vortex around the uncharged droplet is observed when the gel medium is considered to be uncharged, which is similar to the induced charge electroosmosis around an uncharged dielectric colloid in free-solution. We find that the induced charge electrokinetic mechanism creates a strong recirculation of liquid within the droplet and the translational velocity of the droplet strongly depends on its size for the dielectric droplet embedded in a charged gel medium.     

Computational study of ion distributions at the air/liquid methanol interface

Molecular dynamic simulations with polarizable potentials were performed to systematically investigate the distribution of NaCl, NaBr, NaI, and SrCl(2) at the air/liquid methanol interface. The density profiles indicated that there is no substantial enhancement of anions at the interface for the NaX systems, in contrast to what was observed at the air/aqueous interface. The surfactant-like shape of the larger more polarizable halide anions, which is part of the reason they are driven to air/aqueous interfaces, was compensated by the surfactant nature of methanol itself. These halide anions had on average an induced dipole of moderate magnitude in bulk methanol. As a consequence, methanol hydroxy groups donated hydrogen bonds to anions where the negatively charged side of the anion induced dipole pointed, and methyl groups interacted with anions where the positively charged side of the anion-induced dipole pointed. Furthermore, salts were found to disrupt the surface structure of methanol. For the neat air/liquid methanol interface, there is relative enhancement of methyl groups at the outer edge of the air/liquid methanol interface in comparison with hydroxy groups, but with the addition of NaX this enhancement was reduced somewhat. Finally, with the additional of salts to methanol, the computed surface potentials decreased, which is in contrast to what is observed in corresponding aqueous systems, where the surface potential increases with the addition of salts. Both of these trends have been indirectly observed with experiments. The surface potential trends were found to be due to the greater propensity of anions for the air/water interface that is not present at the air/liquid methanol interface.