Image charges and dispersion forces in electric double layers: the dependence of wall-wall interactions on salt concentration and surface charge density

The interaction pressure between two planar charged walls is calculated for a range of conditions. The diffuse electric double layers between the two wall surfaces are treated with ion-wall dispersion forces and ionic image charge interactions taken into account. Both these interactions are due to dielectric discontinuities at the surfaces. Ion-ion and ion-image charge correlations are explicitly included. The ion-wall dispersion interactions can give rise to appreciable ion specific effects, which are particularly strong when the counterions to the surfaces are highly polarizable. The mechanisms of these effects are investigated, and their influence on the net interaction pressure between the walls is studied for a range of surface charge densities, strengths of the anion-wall dispersion interaction and bulk electrolyte concentrations. When the strength of the anion-wall dipersion interaction is increased, the pressure generally becomes less repulsive (or more attractive) for positive surfaces. The opposite happens in general for negative surfaces but to a much lesser extent. The effects are largest for large surface charge densities and high electrolyte concentrations. The image charge interactions give rise to a considerable depletion attraction between the walls for low surface charge densities.     

Electrolyte solutions between uncharged walls

The anisotropic hypernetted chain approximation is applied to calculate the ion concentration profiles, the pair distribution functions and the surface-surface interaction between two uncharged slabs immersed in a primitive model electrolyte solution. The dielectric constant of the slabs is much smaller than that of the solvent. The full complexity of this problem is considered without further approximations. At short surface separations the ions are largely expelled from the space between the slabs due to repulsive image charge interactions. The pressure results clearly demonstrate the interdependence of the image charge interactions and the dispersion interactions between the slabs.     

Charge transport in confined concentrated solutions: A minireview

This Minireview summarizes several recent experiments clouding over prevailing theoretical understanding of charge transport behaviors in electrochemical systems; they are nonlinear concentration dependence of ionic conductivity, ultra-long Debye length in ionic liquids, nonmonotonic double layer charging behavior, and anomalous increase in area specific capacitance with decreasing nanopore size. Theoretical activities reveal that nanoconfinement and high concentration exert strong influence on charge distribution and transport via strong ion-ion correlations and ion-wall interactions. By exemplifying where and why classical theories of charge transport fail, we defy the popular point of view that theoretical electrochemistry is well-established and we are left with applications of these theories only.     

Modification of Surface Forces by Metal Ion Adsorption

Abstract

Sorption of ions may lead to variations in interparticle forces and, thus, changes in the stability of colloidal particles. Chemical interactions between metal ions and colloidal particles modify the molecular structure of the surface, the surface charge, and the electrical potential between colloidal particles. These modifications to the surface and to the electrical double layer due to metal ion sorption are reflected in the interaction force between a particle and another surface, which is measured in this study by atomic force microscopy (AFM). Specifically, AFM is used to investigate the sorption of copper ions from aqueous solutions by silica particles. The influence of metal ion concentration and solution ionic strength on surface forces is studied under transient conditions. Results show that as the metal ion concentration is decreased, charge reversal occurs and a longer period of time is required for the system to reach equilibrium. The ionic strength has no significant effect on sorption kinetics. Furthermore, neither metal concentration nor ionic strength exhibits any effect on sorption equilibria, indicating that for the experimental conditions used in this study, the surface sites of the silica particle are fully occupied by copper ions.     

Dispersion and Hydrogen Bonding Rule: Why the Vaporization Enthalpies of Aprotic Ionic Liquids Are Significantly Larger than those of Protic Ionic liquids

It is well known that gas‐phase experiments and computational methods point to the dominance of dispersion forces in the molecular association of hydrocarbons. Estimates or even quantification of these weak forces are complicated due to solvent effects in solution. The dissection of interaction energies and quantification of dispersion interactions is particularly challenging for polar systems such as ionic liquids (ILs) which are characterized by a subtle balance between Coulomb interactions, hydrogen bonding, and dispersion forces. Here, we have used vaporization enthalpies, far‐infrared spectroscopy, and dispersion‐corrected calculations to dissect the interaction energies between cations and anions in aprotic (AILs), and protic (PILs) ionic liquids. It was found that the higher total interaction energy in PILs results from the strong and directional hydrogen bonds between cation and anion, whereas the larger vaporization enthalpies of AILs clearly arise from increasing dispersion forces between ion pairs.     

Hofmeister effects: interplay of hydration, nonelectrostatic potentials, and ion size

The classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloids, and corresponding theories of electrolytes, are unable to explain ion specific forces between colloidal particles quantitatively. The same is true generally, for surfactant aggregates, lipids, proteins, for zeta and membrane potentials and in adsorption phenomena. Even with fitting parameters the theory is not predictive. The classical theories of interactions begin with continuum solvent electrostatic (double layer) forces. Extensions to include surface hydration are taken care of with concepts like inner and outer Helmholtz planes, and "dressed" ion sizes. The opposing quantum mechanical attractive forces (variously termed van der Waals, Hamaker, Lifshitz, dispersion, nonelectrostatic forces) are treated separately from electrostatic forces. The ansatz that separates electrostatic and quantum forces can be shown to be thermodynamically inconsistent. Hofmeister or specific ion effects usually show up above ≈10(-2) molar salt. Parameters to accommodate these in terms of hydration and ion size had to be invoked, specific to each case. Ionic dispersion forces, between ions and solvent, for ion-ion and ion-surface interactions are not explicit in classical theories that use "effective" potentials. It can be shown that the missing ionic quantum fluctuation forces have a large role to play in specific ion effects, and in hydration. In a consistent predictive theory they have to be included at the same level as the nonlinear electrostatic forces that form the skeletal framework of standard theory. This poses a challenge. The challenges go further than academic theory and have implications for the interpretation and meaning of concepts like pH, buffers and membrane potentials, and for their experimental interpretation. In this article we overview recent quantitative developments in our evolving understanding of the theoretical origins of specific ion, or Hofmeister effects. These are demonstrated through an analysis that incorporates nonelectrostatic ion-surface and ion-ion dispersion interactions. This is based on ab initio ionic polarisabilities, and finite ion sizes quantified through recent ab initio work. We underline the central role of ionic polarisabilities and of ion size in the nonelectrostatic interactions that involve ions, solvent molecules and interfaces. Examples of mechanisms through which they operate are discussed in detail. An ab initio hydration model that accounts for polarisabilities of the tightly held hydration shell of "cosmotropic" ions is introduced. It is shown how Hofmeister effects depend on an interplay between specific surface chemistry, surface charge density, pH, buffer, and counterion with polarisabilities and ion size. We also discuss how the most recent theories on surface hydration combined with hydrated nonelectrostatic potentials may predict experimental zeta potentials and hydration forces.     

Mechanisms behind concentration profiles illustrated by charge and concentration distributions around ions in double layers

The mechanisms behind the behaviour of concentration and charge density profiles in diffuse electric double layers are investigated quantitatively for 1:1 and 2:2 electrolytes. This is done by analysing various contributions to the mean force that acts on each ion. The forces are obtained from the calculated ionic charge and concentration distributions around individual ions at various positions in the double layer. These distributions are presented graphically which allows an immediate visual illustration of the mechanisms in action. Some features studied are charge inversion in double layers for divalent aqueous electrolytes, overcompensation of surface charge due to large amounts of physisorbed counterions, ion size effects in the double layer structure and various mechanisms that cause deviations from the predictions of the Poisson–Boltzmann approximation. A major objective of the paper is to present the results in a visual form and explain aspects of modern double layer theory in a simple manner.     

Free energy decomposition analysis of bonding and nonbonding interactions in solution

A free energy decomposition analysis algorithm for bonding and nonbonding interactions in various solvated environments, named energy decomposition analysis-polarizable continuum model (EDA-PCM), is implemented based on the localized molecular orbital-energy decomposition analysis (LMO-EDA) method, which is recently developed for interaction analysis in gas phase [P. F. Su and H. Li, J. Chem. Phys. 130, 074109 (2009)]. For single determinant wave functions, the EDA-PCM method divides the interaction energy into electrostatic, exchange, repulsion, polarization, desolvation, and dispersion terms. In the EDA-PCM scheme, the homogeneous solvated environment can be treated by the integral equation formulation of PCM (IEFPCM) or conductor-like polarizable continuum model (CPCM) method, while the heterogeneous solvated environment is handled by the Het-CPCM method. The EDA-PCM is able to obtain physically meaningful interaction analysis in different dielectric environments along the whole potential energy surfaces. Test calculations by MP2 and DFT functionals with homogeneous and heterogeneous solvation, involving hydrogen bonding, vdW interaction, metal-ligand binding, cation-π, and ionic interaction, show the robustness and adaptability of the EDA-PCM method. The computational results stress the importance of solvation effects to the intermolecular interactions in solvated environments.     

Atomic force microscopy study of cellulose surface interaction controlled by cellulose binding domains

Colloidal probe microscopy has been used to study the interaction between model cellulose surfaces and the role of cellulose binding domain (CBD), peptides specifically binding to cellulose, in interfacial interaction of cellulose surfaces modified with CBDs. The interaction between pure cellulose surfaces in aqueous electrolyte solution is dominated by double layer repulsive forces with the range and magnitude of the net force dependent on electrolyte concentration. AFM imaging reveals agglomeration of CBD adsorbed on cellulose surface. Despite an increase in surface charge owing to CBD binding to cellulose surface, force profiles are less repulsive for interactions involving, at least, one modified surface. Such changes are attributed to irregularity of the topography of protein surface and non-uniform distribution of surface charges on the surface of modified cellulose. Binding double CBD hybrid protein to cellulose surfaces causes adhesive forces at retraction, whereas separation curves obtained with cellulose modified with single CBD show small adhesion only at high ionic strength. This is possibly caused by the formation of the cross-links between cellulose surfaces in the case of double CBD.     

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.     

Attraction of iodide ions by the free water surface, revealed by simulations with a polarizable force field based on Drude oscillators

Recent theoretical and experimental studies have shown that polarizable anions, such as iodide and bromide, preferentially accumulate close to the surface of electrolyte solutions. This finding is in sharp contrast to the previously prevailing idea that salts are dielectrically excluded from the free water surface and opens up new avenues for research in specific salt effects. In this work, we have verified the ability of a recently introduced polarizable water model, SWM4-DP, to reproduce this behavior, by simulations of a NaI/water slab, corresponding to a 1.2 M solution. The water and ion polarizabilities are modeled by classical Drude oscillator particles. As revealed by the simulations, a double layer is formed close to the free water surface, with the iodide ions located closer to the interface and the sodium ions at a neighboring, interior layer. Near the surface, all solution species acquire an induced dipole moment, that is perpendicular to the surface and points toward the exterior. The double charge layer causes ordering of water at a subsurface region. Simulations with a simpler system of a single iodide ion in a water slab show that the surface position is stabilized by induced charge interactions; in contrast, the charge-dipole interactions between the iodide permanent charge and the water permanent dipole moment favor the bulk position. Thus, the polarizabilities of ion and water are essential for explaining the increased preference of iodide for the air-water interface, in accordance with other studies.     

The Double‐Faced Nature of Hydrogen Bonding in Hydroxy‐Functionalized Ionic Liquids Shown by Neutron Diffraction and Molecular Dynamics Simulations

We characterize the double‐faced nature of hydrogen bonding in hydroxy‐functionalized ionic liquids by means of neutron diffraction with isotopic substitution (NDIS), molecular dynamics (MD) simulations, and quantum chemical calculations. NDIS data are fit using the empirical potential structure refinement technique (EPSR) to elucidate the nearest neighbor H???O and O???O pair distribution functions for hydrogen bonds between ions of opposite charge and the same charge. Despite the presence of repulsive Coulomb forces, the cation–cation interaction is stronger than the cation–anion interaction. We compare the hydrogen‐bond geometries of both “doubly charged hydrogen bonds” with those reported for molecular liquids, such as water and alcohols. In combination, the NDIS measurements and MD simulations reveal the subtle balance between the two types of hydrogen bonds: The small transition enthalpy suggests that the elusive like‐charge attraction is almost competitive with conventional ion‐pair formation.     

Effect of ionic van der Waals forces on the diffuse potential of model colloids

In this work, the role of ionic dispersion forces in specific ion effects is evaluated through Monte Carlo simulations in the primitive model. More specifically, we assess the effect of such forces on the diffuse potential, since this property is essential to understand the electrokinetic behavior of colloids. In this way, ion specificity arises naturally since ion dispersion forces depend on ionic polarizability, which differs for ions with the same valence. This property is included in the primitive model by means of the Lifshitz theory. The results for different ions are summarized with the help of a nondimensional parameter characterizing the relative weight of ionic van der Waals interactions. Our data reveal that for small ions with high polarizability the ionic dispersion forces can considerably contribute to the specificity of the diffuse potential. In any case, the specific ion effects due to ion polarizability are strongly influenced by ion size.     

Van der Waals effects on structure and optical properties in organic photovoltaics

We have recently reported that dispersion forces can have drastic effects on the structure and charge separation in organic photovoltaics (OPVs) (Martinez et al., J. Phys. Chem. C 121, 20 134 [2017]). Here we investigate dimer complexes formed by the polymer P3HT and the fullerene derivative PCBM. We show how van der Waals (vdW) interactions affect the geometrical structure, which has strong effects on the electronic structure and UV-Vis absorption spectrum. Time-dependent density functional theory calculations demonstrate that the experimentally observed blue-shift of the absorption maximum of P3HT/PCBM OPV cells with respect to pure P3HT results mainly from distortions in P3HT due to vdW interactions between donor and acceptor fragments. Reduced absorption in the red region of the UV-Vis spectrum results from distortions of P3HT and small charge transfer between P3HT and PCBM. These results are in qualitative agreement with experiments and recent theoretical results on the corresponding solid-state films.     

Effects of ionic strength and surface charge on protein adsorption at PEGylated surfaces

PEGylated Nb2O5 surfaces were obtained by the adsorption of poly(L-lysine)-g-poly(ethylene glycol) (PLL-g-PEG) copolymers, allowing control of the PEG surface density, as well as the surface charge. PEG (MW 2 kDa) surface densities between 0 and 0.5 nm(-2) were obtained by changing the PEG to lysine-mer ratio in the PLL-g-PEG polymer, resulting in net positive, negative and neutral surfaces. Colloid probe atomic force microscopy (AFM) was used to characterize the interfacial forces associated with the different surfaces. The AFM force analysis revealed interplay between electrical double layer and steric interactions, thus providing information on the surface charge and on the PEG layer thickness as a function of copolymer architecture. Adsorption of the model proteins lysozyme, alpha-lactalbumin, and myoglobin onto the various PEGylated surfaces was performed to investigate the effect of protein charge. In addition, adsorption experiments were performed over a range of ionic strengths, to study the role of electrostatic forces between surface charges and proteins acting through the PEG layer. The adsorbed mass of protein, measured by optical waveguide lightmode spectroscopy (OWLS), was shown to depend on a combination of surface charge, protein charge, PEG thickness, and grafting density. At high grafting density and high ionic strength, the steric barrier properties of PEG determine the net interfacial force. At low ionic strength, however, the electrical double layer thickness exceeds the thickness of the PEG layer, and surface charges "shining through" the PEG layer contribute to protein interactions with PLL-g-PEG coated surfaces. The combination of AFM surface force measurements and protein adsorption experiments provides insights into the interfacial forces associated with various PEGylated surfaces and the mechanisms of protein resistance.     

The Effect of Surface Charge on Adsorption of a Cationic Polyelectrolyte

The regularities of adsorption of a cationic polyelectrolyte, poly(diallyldimethylammonium chloride), on the surface of fused quartz are studied at different values of solution pH by capillary electrokinetics. It is shown that the polyelectrolyte adsorption on a negatively charged surface depends on the value of the surface charge and increases with its growth. At a low charge value (pH 3.8), the polyelectrolyte adsorption increases the quartz surface charge. The driving forces of the adsorption are both electrostatic interaction and forces of nonelectrostatic nature, probably hydrophobic interactions and a change in entropy due to the displacement of counterions from a double layer. The adsorption of poly(diallyldimethylammonium chloride) on quartz from alkaline and neutral solutions is irreversible, which indicates the key role of the electrostatic interaction. At low values of the surface charge, the nonelectrostatic interactions play the main role, thereby resulting in polyelectrolyte desorption.     

Monte Carlo simulation of electrical double-layer formation from mixtures of electrolytes inside nanopores

The formation of the electrical double layer (EDL) in the presence of trivalent and monovalent ions inside a slit-type nanopore was simulated via the canonical Monte Carlo method using a primitive model. In large pores, the distribution of ionic species is similar to that observed in an isolated planar double layer. Screening of surface charge is determined by the competitive effects between ion size and charge asymmetry of the counterions. On the other hand, as the pore size approaches the dimension of the ionic species, phenomena such as EDL overlapping become enhanced by ion-size effects. Simulation results demonstrate that EDL overlapping is not only a function of such parameters as ionic strength and surface charge density, but also a function of the properties of the ionic species involved in the EDL. Furthermore, charge inversion can be observed under certain conditions when dealing with mixtures of asymmetric electrolytes. This phenomenon results from strong ion-ion correlation effects and the asymmetries in size and charge of ionic species, and is most significant in the case of trivalent counterions with larger diameters. The simulation results provide insights into the fundamental mechanisms behind the formation of EDL within nanopores as determined by pore size and by the properties of ionic species present in solution. The findings of this work are relevant to ion sorption and transport within nanostructured materials.     

Gradients of the polarization energy in the effective fragment potential method

The effective fragment potential (EFP) method is an ab initio based polarizable classical method in which the intermolecular interaction parameters are obtained from preparative ab initio calculations on isolated molecules. The polarization energy in the EFP method is modeled with asymmetric anisotropic dipole polarizability tensors located at the centroids of localized bond and lone pair orbitals of the molecules. Analytic expressions for the translational and rotational gradients (forces and torques) of the EFP polarization energy have been derived and implemented. Periodic boundary conditions (the minimum image convention) and switching functions have also been implemented for the polarization energy, as well as for other EFP interaction terms. With these improvements, molecular dynamics simulations can be performed with the EFP method for various chemical systems.     

Electrocapillary at contact: potential-dependent adhesion between a gold electrode and a mica surface

Using the electrochemical surface forces apparatus, we investigated adhesion (from pull-off measurements) between gold and mica as the potential of the gold surface was changed externally. Measurements were performed at different concentrations of KClO(4) in a potential window where the gold electrode is ideally polarizable. At applied potentials where the gold-mica interactions are repulsive, we obtain double layer forces that are predictable by the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloid stability but deviate from the theory at short range. At applied potentials where the gold-mica interactions are attractive, we observed a very strong dependence of adhesion on the applied potential, a result that cannot be directly related to DLVO theory. We show, however, that an approach based on electrocapillary thermodynamics can be employed to model the potential dependence of adhesion seen in our measurements. This electrocapillary approach presents evidence of charging at the gold-mica interface and stresses the relation between the charge within and outside of the contact area.