Douwe Jan Bonthuis

Nanoscale Pumping of Water by AC Electric Fields

Using molecular dynamics simulations we demonstrate pumping of water through a carbon nanotube by time-dependent electric fields. The fields are generated by electrodes with oscillating charges in a broad gigahertz frequency range that are attached laterally to the tube. The key ingredient is a phase shift between the electrodes to break the spatiotemporal symmetry. A microscopic theory based on a polarization-dragging mechanism accounts quantitatively for our numerical findings.

Profile of the static permittivity tensor of water at interfaces: consequences for capacitance, hydration interaction and ion adsorption.

We derive the theoretical framework to calculate the dielectric response tensor and determine its components for water adjacent to hydrophilic and hydrophobic surfaces using molecular dynamics simulations. For the nonpolarizable water model used, linear response theory is found to be applicable up to an external perpendicular field strength of ∼2 V/nm, which is well beyond the experimental dielectric breakdown threshold. The dipole contribution dominates the dielectric response parallel to the interface, whereas for the perpendicular component it is essential to keep the quadrupole and octupole terms. Including the space-dependent dielectric function in a mean-field description of the ion distribution at a single charged interface, we reproduce experimental values of the interfacial capacitance. At the same time, the dielectric function decreases the electrostatic part of the disjoining pressure between two charged surfaces, unlike previously thought. The difference in interfacial polarizability between hydrophilic and hydrophobic surfaces can be quantized in terms of the dielectric dividing surface. Using the dielectric dividing surface and the Gibbs dividing surface positions to estimate the free energy of a single ion close to an interface, ion-specific adsorption effects are found to be more pronounced at hydrophobic surfaces than at hydrophilic surfaces, in agreement with experimental trends.

Beyond the continuum: how molecular solvent structure affects electrostatics and hydrodynamics at solid-electrolyte interfaces.

Standard continuum theory fails to predict several key experimental results of electrostatic and electrokinetic measurements at aqueous electrolyte interfaces. In order to extend the continuum theory to include the effects of molecular solvent structure, we generalize the equations for electrokinetic transport to incorporate a space dependent dielectric profile, viscosity profile, and non-electrostatic interaction potential. All necessary profiles are extracted from atomistic molecular dynamics (MD) simulations. We show that the MD results for the ion-specific distribution of counterions at charged hydrophilic and hydrophobic interfaces are accurately reproduced using the dielectric profile of pure water and a non-electrostatic repulsion in an extended Poisson-Boltzmann equation. The distributions of Na(+) at both surface types and Cl(-) at hydrophilic surfaces can be modeled using linear dielectric response theory, whereas for Cl(-) at hydrophobic surfaces it is necessary to apply nonlinear response theory. The extended Poisson-Boltzmann equation reproduces the experimental values of the double-layer capacitance for many different carbon-based surfaces. In conjunction with a generalized hydrodynamic theory that accounts for a space dependent viscosity, the model captures the experimentally observed saturation of the electrokinetic mobility as a function of the bare surface charge density and the so-called anomalous double-layer conductivity. The two-scale approach employed here-MD simulations and continuum theory-constitutes a successful modeling scheme, providing basic insight into the molecular origins of the static and kinetic properties of charged surfaces, and allowing quantitative modeling at low computational cost.

Unraveling the Combined Effects of Dielectric and Viscosity Profiles on Surface Capacitance, Electro-Osmotic Mobility, and Electric Surface Conductivity

We calculate the electro-osmotic mobility and surface conductivity at a solid-liquid interface from a modified Poisson-Boltzmann equation, including spatial variations of the dielectric function and the viscosity that where extracted previously from molecular dynamics simulations of aqueous interfaces. The low-dielectric region directly at the interface leads to a substantially reduced surface capacitance. At the same time, ions accumulate into a highly condensed interfacial layer, leading to the well-known saturation of the electro-osmotic mobility at large surface charge density regardless of the hydrodynamic boundary conditions. The experimentally well-established apparent excess surface conductivity follows from our model for all hydrodynamic boundary conditions without additional assumptions. Our theory fits multiple published sets of experimental data on hydrophilic and hydrophobic surfaces with striking accuracy, using the nonelectrostatic ion-surface interaction as the only fitting parameter.

Self-Energy-Limited Ion Transport in Subnanometer Channels

The current-voltage characteristics of the alpha-Hemolysin protein pore during the passage of single-stranded DNA under varying ionic strength C are studied experimentally. We observe strong blockage of the current, weak superlinear growth of the current as a function of voltage, and a minimum of the current as a function of C. These observations are interpreted as the result of the ion electrostatic self-energy barrier originating from the large difference in the dielectric constants of water and the lipid bilayer. The dependence of DNA capture rate on C also agrees with our model.

Theory and simulations of water flow through carbon nanotubes: prospects and pitfalls

We study water flow through carbon nanotubes using continuum theory and molecular dynamics simulations. The large slip length in carbon nanotubes greatly enhances the pumping and electrokinetic energy conversion efficiency. In the absence of mobile charges, however, the electro-osmotic flow vanishes. Uncharged nanotubes filled with pure water can therefore not be used as electric field-driven pumps, contrary to some recently ventured ideas. This is in agreement with results from a generalized hydrodynamic theory that includes the angular momentum of rotating dipolar molecules. The electro-osmotic flow observed in simulations of such carbon nanotubes is caused by an imprudent implementation of the Lennard-Jones cutoff. We also discuss the influence of other simulation parameters on the spurious electro-osmotic flow.

Electrokinetics at Aqueous Interfaces without Mobile Charges

We theoretically consider the possibility of using electric fields in aqueous channels of cylindrical and planar geometry to induce transport in the absence of mobile ionic charges. Using the Navier-Stokes equation, generalized to include the effects of water spinning, dipole orientation and relaxation, we show analytically that pumping of a dipolar liquid through an uncharged hydrophobic channel can be achieved by injecting torque into the liquid, based on the coupling between molecular spinning and fluid vorticity. This is possible using rotating electric fields and suitably chosen interfacial boundary conditions or transiently by suddenly switching on a homogeneous electric field. A static electric field, however, does not induce a steady state flow in channels, irrespective of the geometry. Using molecular dynamics (MD) simulations, we confirm that static fields do not lead to any pumping, in contrast to earlier publications. The pumping observed in MD simulations of carbon nanotubes and oil droplets in a static electric field is tracked down to an imprudent implementation of Lennard-Jones interaction truncation schemes.

Optimization of classical nonpolarizable force fields for OH− and H3O+

We optimize force fields for H3O(+) and OH(-) that reproduce the experimental solvation free energies and the activities of H3O(+) Cl(-) and Na(+) OH(-) solutions up to concentrations of 1.5 mol/l. The force fields are optimized with respect to the partial charge on the hydrogen atoms and the Lennard-Jones parameters of the oxygen atoms. Remarkably, the partial charge on the hydrogen atom of the optimized H3O(+) force field is 0.8 ± 0.1|e|--significantly higher than the value typically used for nonpolarizable water models and H3O(+) force fields. In contrast, the optimal partial charge on the hydrogen atom of OH(-) turns out to be zero. Standard combination rules can be used for H3O(+) Cl(-) solutions, while for Na(+) OH(-) solutions, we need to significantly increase the effective anion-cation Lennard-Jones radius. While highlighting the importance of intramolecular electrostatics, our results show that it is possible to generate thermodynamically consistent force fields without using atomic polarizability.

Nanoscale structure of the oil-water interface

X-ray reflectivity (XR) and atomistic molecular dynamics (MD) simulations, carried out to determine the structure of the oil-water interface, provide new insight into the simplest liquid-liquid interface. For several oils (hexane, dodecane, and hexadecane) the XR shows very good agreement with a monotonic interface-normal electron density profile (EDP) broadened only by capillary waves. Similar agreement is also found for an EDP including a sub-Å thick electron depletion layer separating the oil and the water. The XR and MD derived depletions are much smaller than reported for the interface between solid-supported hydrophobic monolayers and water.

Mechanosensitive channel activation by diffusio-osmotic force

For ion channel gating, the appearance of two distinct conformational states and the discrete transitions between them are essential, and therefore of crucial importance to all living organisms. We show that the physical interplay between two structural elements that are commonly present in bacterial mechanosensitive channels--namely, a charged vestibule and a hydrophobic constriction--creates two distinct conformational states, open and closed, as well as the gating between them. We solve the nonequilibrium Stokes-Poisson-Nernst-Planck equations, extended to include a molecular potential of mean force, and show that a first order transition between the closed and open states arises naturally from the diffusio-osmotic stress caused by the ions and the water inside the channel and the elastic restoring force from the membrane.