Aqueous interfaces: Bridging the gap between continuum and molecular models

The particulars of the electrostatic and hydrodynamic interactions at the interface between charged objects and water are of immense importance for biological systems, as well as for a wide range of technological applications, such as gel electrophoresis and electrokinetic pumping. Despite the progress made in over a century of active research, the intricate structure of aqueous interfaces remain poorly understood. Using a two-scale modelling approach, we study the wide- ranging effects of the ordered and anisotropic molecular structure of interfacial water.

We extract molecular profiles of the dielectric function, the viscosity and the ionic free energy from atomistic molecular dynamics simulations, and incorporate this molecular information into continuum theory. The use of continuum theory enables us to reach length and time scales that are unattainable using only molecular dynamics simulations. At the same time, it provides us with the flexibility to easily vary system parameters such as geometry, salt concentration and surface charge density. One of the most significant advantages of this combination of molecular detail, flexibility of model parameters and sufficiently large scale is that it allows for direct comparison with experimental data.

Incorporating the dielectric function and the viscosity profile into the combined Stokes and Poisson-Boltzmann equations accurately captures the saturation of the electro-osmotic flow velocity as a function of the bare surface charge density, which is consistently found in experiments on hydrophobic and hydrophilic surfaces (Fig. 1b).

srf-covers
Cover image source: American Chemical Society.
srf-puzzles
Fig. 1. (a) An electric field parallel to a surface carrying a bare surface charge density σ0 produces an electro-osmotic flow profile, from which the electrokinetic surface charge density σek is calculated. The conductive surface charge density σc is calculated from the electric current I. (b) σek is experimentally found to saturate as a function of σ0; saturation occurs at higher values for hydrophobic surfaces than for hydrophilic surfaces. (c) σc exceeds σek at all surface types. (d) When the bulk salt concentration c0 is raised, σek increases at hydrophobic surfaces, and decreases at hydrophilic surfaces.
Moreover, the experimentally well-established excess surface conductivity (Fig. 1c), generally referred to as anomalous surface conductivity, follows directly from the same model. The latter result makes the assumption of electrical conductance through otherwise static water layers — which has been the common explanation of the experimental observations — superfluous. Finally, the dependence of the electro-osmotic flow velocity on the bulk salt concentration (Fig. 1d) is explained by the same combination of dielectic and hydrodynamic surface effects [1,2].

Our multiscale modeling approach thus provides us with the means for experimental verification of our theoretical models, and simultaneously offers new insight into the molecular origin of experimental observations.
[2] D. J. Bonthuis and R. R. Netz; Unraveling the combined effects of dielectric and viscosity profiles on surface capacitance, electro-osmotic mobility, and electric surface conductivity; Langmuir 28, 16049 (2012)
[1] D. J. Bonthuis and R. R. Netz; Beyond the continuum: How molecular solvent structure affects electrostatics and hydrodynamics at solid-electrolyte interfaces; J. Chem. Phys. B 117, 11397 (2013)