The Dramatic Effect of Water Structure on Hydration Forces and the Electrical Double Layer

Abstract

Forces between hydrophilic surfaces mediated by water are important in various systems from lipid membranes and solid surfaces to colloids and macromolecules, first discovered as a significant addition to DLVO forces at the nanoscale. These “hydration forces” have been studied in great detail experimentally using osmotic stress measurements, surface force apparatus, and AFM, and they have also been the subject of multiple theories and simulations. One spectacular feature observed in experimental and simulation studies was the nonmonotonic, oscillatory decay in the forces between atomically smooth surfaces. Forces between “rougher” surfaces exhibit only quasi-exponential, monotonic decay. Here we revisit this hydration force problem by exploring the consequences of an extended phenomenological Landau–Ginzburg approach that describes nonlocal correlations in water, linking them with the key features of the wave-number k-dependent nonlocal dielectric function of water. With corresponding boundary conditions, this theory predicts the observed oscillatory decay in hydration force between ideally flat surfaces, the oscillatory mode disappearing with just a tiny roughness of the surfaces (of mean height ca. of the size of a water molecule). This study also brings an important side message. Explanation of these observations appears only possible under an assumption of two modes of polarization in water, consistent with the behavior of the response function, i.e., Lorentzian at small k and resonance-like at higher k. This resolves the “force oscillation–non-oscillation” paradigm, which is a strong, although indirect indication of the existence of these two modes. We also consider other important subjects, such as how the distribution of ions near a charged surface reacts to the propensity for overscreening oscillations due to polarized water. This is important not only for the interactions between charged surfaces but also for the fundamental understanding of the structure of the electrical double layer at electrochemical interfaces. We show that even in dilute electrolytes, the distribution of ions in the vicinity of the polarized interface follows, although not literally, preferential positions corresponding to the potential wells caused by “resonance” water layering. For a sharp interface, the theory predicts that the decaying spatial oscillation profiles extend over a 1 to 2 nm distance from the interface. With the smearing of the interface and the corresponding suppression of the resonance water layering, oscillations in the spatial distribution of ions subside, resulting in a familiar Gouy–Chapman–Stern picture. At longer distances from the interface, whether smeared or not, the ion distribution profiles become Gouy–Chapman-like. The effect of the boundary conditions on water polarization at the interface goes beyond a trivial shift of the potential of zero charge. We show that they can dramatically affect the ion distribution near the charged surface. Last, but not least, we study how the interfacial water layerin...

Orientational Ordering in Nano-confined Polar Liquids

“…It has also recently been shown that the smearing of the surface (and of the boundary condition) may dramatically affect the polarization profiles; a strong and nontrivial effect is also expected to come from the lateral inhomogeneity of the surfaces . In the present study we have considered ideally smooth, laterally homogeneous, and sharp surfaces, whereas in reality in many systems they are not as such.…”

Multiscale Modeling of Aqueous Electric Double Layers

“…Third, density variations have a direct effect on the adsorption of ions due to steric exclusion . In particular, FF MD simulations show that ions at charged hydrophobic surfaces accumulate between the peaks in the oscillating water density profile, an effect expected to vanish for surfaces exhibiting a roughness of as little as the size of a molecule . In these simulated density profiles, the traditional interfacial layers can be clearly identified, viz., partially dehydrated ions between the surface and the first density peak form the inner Helmholtz layer of adsorbed ions, fully hydrated ions between the first two density peaks form the outer Helmholtz layer, and the rest of the fluid forms the diffuse layer.…”

“…Water-mediated hydration forces are caused not only by molecular density fluctuations. The coupling between polarization and molecular density fluctuations gives rise to a resonance peak in the wave-vector-dependent nonlocal dielectric response function, which in turn causes an oscillatory decay in the water-mediated forces between smooth surfaces . Yet to show that the oscillatory structure is present also in the absence of the AFM tip, the interfacial density of water oxygen can be extracted from a fit to X-ray spectroscopy data.…”

“…In particular, the hydrophilic diamond surface features protruding hydroxide groups, and the silica surface consists of alternating silicon and oxygen atoms as well as hydroxide groups. Classical density functional theory has shown that smearing over the size of about 0.25 nm is already sufficient to suppress the oscillations in the ion density and hydration force resulting from the nonlocal dielectric response …”

Properties of aqueous electrolyte solutions at carbon electrodes: effects of concentration and surface charge on solution structure, ion clustering and thermodynamics in the electric double layer

“…More sophisticated mean-field models of the EDL were developed to address some of these effects. 15–21…”