Hydration forces between silica surfaces: experimental data and predictions from different theories

Silica is a very interesting system that has been thoroughly studied in the last decades. One of the most outstanding characteristics of silica suspensions is their stability in solutions at high salt concentrations. In addition to that, measurements of direct-interaction forces between silica surfaces, obtained by different authors by means of surface force apparatus or atomic force microscope (AFM), reveal the existence of a strong repulsive interaction at short distances (below 2 nm) that decays exponentially. These results cannot be explained in terms of the classical Derjaguin, Landau, Verwey, and Overbeek (DLVO) theory, which only considers two types of forces: the electrical double-layer repulsion and the London-van der Waals attraction. Although there is a controversy about the origin of the short-range repulsive force, the existence of a structured layer of water molecules at the silica surface is the most accepted explanation for it. The overlap of structured water layers of different surfaces leads to repulsive forces, which are known as hydration forces. This assumption is based on the very hydrophilic nature of silica. Different theories have been developed in order to reproduce the exponentially decaying behavior (as a function of the separation distance) of the hydration forces. Different mechanisms for the formation of the structured water layer around the silica surfaces are considered by each theory. By the aid of an AFM and the colloid probe technique, the interaction forces between silica surfaces have been measured directly at different pH values and salt concentrations. The results confirm the presence of the short-range repulsion at any experimental condition (even at high salt concentration). A comparison between the experimental data and theoretical fits obtained from different theories has been performed in order to elucidate the nature of this non-DLVO repulsive force.     

Forces between surfaces in liquids

The results and implications of direct force measurements between molecularly smooth mica surfaces in liquids are reviewed. These discussions include four interactions fundamental to colloid science: van der Waals forces, double layer forces, adhesion forces and structural or solvation forces (e.g. hydration forces). Also considered are the effects of preferential surface adsorption of solute molecules on these interactions, e.g. surfactant adsorptions from aqueous solutions and water condensation from non-aqueous solvents.In aqueous media it is apparent that the DLVO theory is valid at all surface separations down to the “force barrier”, but that under certain conditions hydration forces can become significant at distances below 30 Å.The measured adhesion force between two solid surfaces can be simply related to their surface energies and where meniscus forces are also present due to “capillary condensation” from vapor solvent, their effect on adhesion can be understood in terms of straightforward bulk thermodynamic principles. Here, too, it is concluded that structural forces cannot be ignored.Our results suggest that structural forces may either very monotonically with distance or be oscillatory with a periodicity equal to the molecular size. Their origin, nature, mode of action and importance for particle interactions will no doubt take many years to sort out.     

Approaches to hydration, old and new: Insights through Hofmeister effects

Hydration effects in colloidal interactions or problems involving electrolytes are usually taken care of by effective electrostatic potentials that subsume notions like hydrated ion size, Gurney potentials, soft and hard, chaotropic and cosmotropic ions, and inner and outer Helmholtz planes. Quantum fluctuation (dispersion) forces between ions and between ions and surfaces are missing from classical theories, at least not explicit in standard approaches to hydration. This paper outlines an evolving back-to-basics approach that allows these ion specific forces to be included in theories quantitatively. In this approach ab initio quantum mechanics is used to calculate dynamic polarisabilities of ions and to quantify bare ion radii. The ionic dispersion potentials between ions, and between ions and surfaces in water can then be given explicit analytic form from an extension of Lifshitz theory. They are included in the theory along with electrostatic potentials. In a first stage the primitive (continuum solvent) model provides a skeletal theory on which to build in hydration. Extension of the ab initio calculations to include “dressed” ions, i.e. water hydration shells for cosmotropic ions, quadrupolar and octupolar polarisability contributions and; for colloids, allowance for a surface hydration layer, permits quantification of Hofmeister effects and Gurney potentials. With these extensions, primary hydration forces (short range repulsion) arise due to an interplay between surface hydration layers and specific ion interactions. Apparent longer range “secondary hydration forces” are shown to be a consequence of ion-surface dispersion interactions and are not true “hydration forces”.     

Structure and dynamics of a fluid phase bilayer on a solid support as observed by a molecular dynamics computer simulation

Simulations of a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine lipid bilayer interacting with a solid surface of hydroxylated nanoporous amorphous silica have been carried out over a range of lipid-solid substrate distances. The porous solid surface allowed the water layer to dynamically adjust its thickness, maintaining equal pressures above and below the membrane bilayer. Qualitative estimates of the force between the surfaces leads to an estimated lipid-silicon distance in very good agreement with the results of neutron scattering experiments. Detailed analysis of the simulation at the separation suggested by experiment shows that for this type of solid support the water layer between surfaces is very narrow, consisting only of bound waters hydrating the lipid head groups and hydrophilic silica surface. The reduced hydration, however, has only minor effects on the head group hydration, the orientation of water molecules at the interface, and the membrane dipole potential. Whereas these structural properties were not sensitive to the presence of the solid substrate, the calculated diffusion coefficient for translation of the lipid molecules was altered significantly by the silica surface.     

Formation of ice-like water structure on the surface of an antifreeze protein

Antifreeze proteins (AFPs) are found in different species from polar, alpine, and subarctic regions where they serve to inhibit ice crystal growth by adsorption to ice surfaces. Computational methods have the power to investigate the antifreeze mechanism in atomic detail. Molecular dynamics simulations of water under different conditions have been carried out to test our water model for simulations of biological macromolecules in extreme conditions: very low temperatures (200 K) and at the ice/liquid water interface. We show that the flexible F3C water model reproduces properties of water in the solid phase (ice I(h)), the supercooled liquid phase, and at the ice/liquid water interface. Additionally, the hydration of the type III AFP from ocean pout was studied as a function of temperature. Hydration waters on the ice-binding surface of the AFP were less distorted and more tetrahedral than elsewhere on the surface. More ice-like hydrating water structures formed on the ice-binding surface of the protein such that it created an ice-like structure in water within its first hydration layer but not beyond, suggesting that this portion of the protein has high affinity for ice surfaces.     

Structured water in partially dehydrated yeast cells and at partially hydrophobized fumed silica surface

Nonfreezable water structured due to interaction with endocellular objects in yeast cells (endocellular water) or with partially hydrophobic fumed silica (interfacial water) was studied by means of (1)H NMR spectroscopy with layer-by-layer freezing-out of bulk water and quantum chemical methods applied to water clusters in the gas and liquid (chloroform and cyclohexene) phases and adsorbed on silylated silica. Variation in cell hydration as well as in amount of water adsorbed on modified fumed silica leads to changes in the ratio between contributions of two water states characterized by NMR chemical shifts at delta(H)=1.1-1.7 and 4-5 ppm. Lowering of hydration and temperature results in an enhancement of the first signal, and the opposite result is observed for the second signal. These effects may be explained by structured water distribution in the form of relatively large nanodroplets (delta(H)=4-5 ppm is close to that for bulk water) and small clusters of the 2D structure, in which the fraction of hydrogen atoms out of the hydrogen bonds (delta(H)=1.1-1.7 ppm) is larger than that in nanodroplets.     

Supercrystal structures of polyhedral PbS nanocrystals

We analyzed the interaction between chemically grafted polysaccharide layers in aqueous solutions. To fabricate such layers, an end-terminated dextran silane coupling agent was synthesized and the polydextran was grafted to oxidized silicon wafers and to silica particles. This resulted in the formation of a 28 nm thick layer (in air) and a grafted amount of 40 mg/m(2) as determined by ellipsometry. The physical properties of the grafted layer were investigated in aqueous solutions by atomic force microscope imaging and colloidal probe force measurements. Surface and friction forces were measured between one bare and one polydextran coated silica surface. A notable feature was a bridging attraction due to affinity between dextran and the silica surface. Surface interactions and friction forces were also investigated between two surfaces coated with grafted polydextran. Repulsive forces were predominant, but nevertheless a high friction force was observed. The repulsive forces were enhanced by addition of sodium dodecyl sulfate (SDS) that associates with the tethered polydextran layers. SDS also decreased the friction force. Our data suggests that energy dissipation due to shear-induced structural changes within the grafted layer is of prime importance for the high friction forces observed, in particular deformation of protrusions in the surface layer.     

The polarization model for hydration/double layer interactions: the role of the electrolyte ions

The interactions between hydrophilic surfaces in water cannot be always explained on the basis of the traditional Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, and an additional repulsion, the "hydration force" is required to accommodate the experimental data. While this force is in general associated with the organization of water in the vicinity of the surface, different models for the hydration were typically required to explain different experiments. In this article, it is shown that the polarization-model for the double layer/hydration proposed by the authors can explain both (i) the repulsion between neutral lipid bilayers, with a short decay length ( approximately 2 A), which is almost independent of the electrolyte concentration, and, at the same time, (ii) the repulsion between weakly charged mica surfaces, with a longer decay length ( approximately 10 A), exhibiting not only a dependence on the ionic strength, but also strong ion-specific effects. The model, which was previously employed to explain the restabilization of protein-covered latex particles at high ionic strengths and the existence of a long-range repulsion between the apoferritin molecules at moderate ionic strengths, is extended to account for the additional interactions between ions and surfaces, not included in the mean field electrical potential. The effect of the disorder in the water structure on the dipole correlation length is examined and the conditions under which the results of the polarization model are qualitatively similar to those obtained by the traditional theory via parameter fitting are emphasized. However, there are conditions under which the polarization model predicts results that cannot be recovered by the traditional theory via parameter fitting.     

Hofmeister effects in the restabilization of IgG--latex particles: testing Ruckenstein's theory

The hydration interaction is responsible for the colloidal stability observed in protein-coated particles at high ionic strengths. The origin of this non-DLVO interaction is related not only to the local structure of the water molecules located at the surface but also to the structure of those molecules involved in the hydration of the ions that surround the colloidal particles. Ruckenstein and co-workers have recently developed a new theory based on the coupling of double-layer and hydration interactions. Its validity was contrasted by their fitting of experimental data obtained with IgG-latex particles restabilized at high salt concentration. The theory details the important role played by the counterions in the stability at high salt concentrations by proposing an ion pair reaction forming surface dipoles. These surface dipoles are responsible of repulsive interactions between two approaching surfaces. This paper checks the theory with recent data where some ions associated with the Hofmeister series (NO(3)(-), SCN(-) and Ca(2+)) restabilize the same kind of IgG-latex systems by means of hydration forces. Surprisingly, these ions induce stability acting even as co-ions, likely by modifying the water structure at the surface, but not forming surface ion pairs. Therefore, this experimental evidence would question Ruckenstein's theory based on the surface dipole formation for explaining the observed restabilization phenomena.     

Probing colloidal forces between a Si3N4 AFM tip and single nanoparticles of silica and alumina

The atomic force microscope (AFM) has been used to measure surface forces between silicon nitride AFM tips and individual nanoparticles deposited on substrates in 10(-4) and 10(-2) M KCl solutions. Silica nanoparticles (10 nm diameter) were deposited on an alumina substrate and alumina particles (5 to 80 nm diameter) were deposited on a mica substrate using aqueous suspensions. Ionic concentrations and pH were used to manage attractive substrate-particle electrostatic forces. The AFM tip was located on deposited nanoparticles using an operator controlled offset to achieve stepwise tip movements. Nanoparticles were found to have a negligible effect on long-range tip-substrate interactions, however, the forces between the tip and nanoparticle were detectable at small separations. Exponentially increasing short-range repulsive forces, attributed to the hydration forces, were observed for silica nanoparticles. The effective range of hydration forces was found to be 2-3 nm with the decay length of 0.8-1.3 nm. These parameters are in a good agreement with the results reported for macroscopic surfaces of silica obtained using the surface force apparatus suggesting that hydration forces for the silica nanoparticles are similar to those for flat silica surfaces. Hydration forces were not observed for either alumina substrates or alumina nanoparticles in both 10(-4) M KCl solution at pH 6.5 and 10(-2) M KCl at pH 10.2. Instead, strong attractive forces between the silicon nitride tip and the alumina (nanoparticles and substrate) were observed.     

Study of hydration layers near a hydrophilic surface in water through AFM imaging

In this work, a hydrophilic silica plate exposed in air, and immersed in an aqueous solution was studied through atomic force microscopy (AFM) imaging in contact‐ and tapping‐mode operations. It was experimentally found that the tapping‐mode AFM images of the silica surface were different when it was immersed in an aqueous solution from those when it was exposed in air. The former showed fewer topographic features. However, the contact‐mode AFM images of the silica surface were almost uninfluenced by the medium in which the surface was placed. This phenomenon might be attributed to the existence of hydration layers near the silica surface in the aqueous solution. The layers are like a large sheet on the surface that hides the details, so that an AFM tip in the tapping mode can read only the hydration layer and therefore image only the rough outline of the surface. This result might suggest the existence of hydration layers near a hydrophilic surface immersed in water. Copyright © 2006 John Wiley & Sons, Ltd.     

On the relationship between charge distribution, surface hydration, and the structure of the interface of metal hydroxides

The double layer structure of metal (hydr)oxides is discussed. Charge separation may exist between the minimum distance of approach of electrolyte ions and the DDL domain. The corresponding capacitance value of the outer Stern layer is similar to the capacitance value of the inner Stern layer. The extended Stern model implicitly supports a hydration structure at the near-surface with some discrete layering of water and electrolyte ions. The significance of dipole orientation is analyzed theoretically. Dipole theory in combination with a calculated ion charge distribution is compared with the experimental overall charge distribution. Ion charge distribution for various oxyanions has been calculated applying the Brown bond valence concept to the geometry of surface complexes that have been optimized with MO/DFT calculations. The comparison is done in detail for silicic acid adsorption on goethite. In addition, results are discussed for arsenite, carbonate, sulfate, and phosphate, using the same approach. The dipole correction depends on the charge introduced in a neutral surface by ion adsorption, which differs for the various ions studied. The fractional correction factor phi derived for the experimental data agrees with the theoretical value phi(m)=0.17+/-0.02. On an absolute scale, the dipole corrections are usually limited to the range about 0-0.15 v.u. The CD values calculated with MO/DFT are not particularly sensitive (approximately 0.03 v.u.) to the precise Fe-octahedral geometry, which suggests that a calculated CD is a reasonable approximation in ion adsorption modeling for ill-defined Fe-oxides like HFO and natural Fe oxide materials of soils.     

Molecular dynamics simulation of nanocolloidal amorphous silica particles: Part I

Explicit molecular dynamics simulations were applied to a pair of amorphous silica nanoparticles in aqueous solution, with diameter of 4.4 nm and with four different background electrolyte concentrations, to extract the mean force acting between the two silica nanoparticles. Dependences of the interparticle forces on the separation and the background electrolyte concentration were demonstrated. The nature of the interaction of the counterions with charged silica surface sites (deprotonated silanols) was investigated. A "patchy" double layer of adsorbed sodium counterions was observed. Dependences of the interparticle potential of mean force on the separation and the background electrolyte concentration were demonstrated. Direct evidence of the solvation forces is presented in terms of changes of the water ordering at the surfaces of the isolated and double nanoparticles. The nature of the interaction of the counterions with charged silica surface sites (deprotonated silanols) was investigated in terms of quantifying the effects of the number of water molecules separately inside each pair of nanoparticles by defining an impermeability measure. A direct correlation was found between the impermeability (related to the silica surface "hairiness") and the disruption of water ordering. Differences in the impermeability between the two nanoparticles are attributed to differences in the calculated electric dipole moment.     

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.     

From simple surface models to lipid membranes: Universal aspects of the hydration interaction from solvent-explicit simulations

A review of atomistic simulation approaches including explicit water for the study of hydration forces between polar surfaces is presented. In particular, we discuss different methods for keeping the chemical potential of water constant and compare advantages and limitations of each method. It turns out that modifications of hydration forces due to surface softness can be accounted for by a convolution over the surface shape profile. Universal aspects of the hydration interaction observed in simulations of different surface chemistries are highlighted, while special attention is given to hydration forces between self-assembled phospholipid membranes.     

Hydration structure of human lysozyme investigated by molecular dynamics simulation and cryogenic X-ray crystal structure analyses: on the correlation between crystal water sites,solvent density,and solvent dipole

The hydration structure of human lysozyme was studied with cryogenic X-ray diffraction experiment and molecular dynamics simulations. The crystal structure analysis at a resolution of 1.4 A provided 405 crystal water molecules around the enzyme. In the simulations at 300 K, the crystal structure was immersed in explicit water molecules. We examined correlations between crystal water sites and two physical quantities calculated from the 1-ns simulation trajectories: the solvent density reflecting the time-averaged distribution of water molecules, and the solvent dipole measuring the orientational ordering of water molecules around the enzyme. The local high solvent density sites were consistent with the crystal water sites, and better correlation was observed around surface residues with smaller conformational fluctuations during the simulations. Solvent dipoles around those sites exhibited coherent and persistent ordering, indicating that the hydration water molecules at the crystal water sites were highly oriented through the interactions with hydrophilic residues. Those water molecules restrained the orientational motions of adjoining water molecules and induced a solvent dipole field, which was persistent during the simulations around the enzyme. The coherent ordering was particularly prominent in and around the active site cleft of the enzyme. Because the ordering was significant up to the third to fourth solvent layer region from the enzyme surface, the coherently ordered solvent dipoles likely contributed to the molecular recognition of the enzyme in a long-distance range. The present work may provide a new approach combining computational and the experimental studies to understand protein hydration.     

The hydration of hydrophobic substances

A model of the hydration of hydrophobic substances in water is suggested. The models of fluctuation formation of empty cavities in water as a stage of hydration extensively used in the literature were shown to be at variance with experiment. The fundamental role played by the interphase boundary surface was emphasized. On this surface, the successive addition of water molecules with the formation of capsules around hydrophobic molecules occurred. The physical meaning of the Ostwald equation was revealed. This equation characterized the distribution of hydrophobic volatile substances between the gas and aqueous phases. The method of optical probes (hydrophobic aromatic molecules) was used to reveal the synergistic character of autocorrelation of dispersion interactions between water and hydrophobic substance molecules. This synergism was at variance with the Lennard-Jones potential. The synergism (superadditivity) of dispersion attraction forces, which strengthened their directional character, caused the self-organization and enhanced stability of hydration capsules with encapsulated hydrophobic molecules. Computer models were used to show that the spatially directional character of dispersion interactions necessary for the self-organization of hydrated aggregates could be simulated by the molecular mechanics method on the basis of orientational correlation of water molecules and hydrophobic substances in the starting system.     

Hydrogen bond lifetime dynamics at the interface of a surfactant monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.     

Protein driven patterning of self-assembled cubosomic nanostructures: long oriented nanoridges

Self-assembly lipid/protein cubosomic nanostructures are generated at high hydration level (dispersion of 5% lipid only) and examined by freeze-fracture electron microscopy (FF-EM) and synchrotron X-ray diffraction (XRD). The fracture surface of the three-dimensional (3D) soft-matter membranous assembly reveals starlike nanopatterns of oriented 100-nm-long cubosomic nanoridges with lateral periodicity defined by their 21-nm diameters. The average water channel radius in these liquid crystalline cubosomic nanoarchitectures, determined by high-resolution FF-EM and XRD, is 18.0 Angstrom. The protein-directed fragmentation of a diamond-type lipid cubic phase at high hydration can induce 3D patterns of oriented nanoporous building blocks, which are a unique example of tertiary organization of functionalized fluid lipid/water interfaces.     

Water density in the electric double layer at the insulator/electrolyte solution interface

I studied the spatial structure of the thick transition region between n-hexane and a colloidal solution of 7-nm silica particles by X-ray reflectivity and grazing incidence small-angle scattering. The interfacial structure is discussed in terms of a semiquantitative interface model wherein the potential gradient at the n-hexane/sol interface reflects the difference in the potentials of "image forces" between the cationic Na(+) and anions (nanoparticles) and the specific adsorption of surface charge at the interface between the adsorbed layer and the solution, as well as at the interface between the adsorbed layer and n-hexane. The X-ray scattering data revealed that the average density of water in the field approximately 10(9)-10(10) V/m of the electrical double layer at the hexane/silica sol interface is the same as, or only few percent higher (1-7%) than, its density under normal conditions.