Liquid-vapor oscillations of water in hydrophobic nanopores

Water plays a key role in biological membrane transport. In ion channels and water-conducting pores (aquaporins), one dimensional confinement in conjunction with strong surface effects changes the physical behavior of water. In molecular dynamics simulations of water in short (0.8 nm) hydrophobic pores the water density in the pore fluctuates on a nanosecond time scale. In long simulations (460 ns in total) at pore radii ranging from 0.35 nm to 1.0 nm we quantify the kinetics of oscillations between a liquid-filled and a vapor-filled pore. This behavior can be explained as capillary evaporation alternating with capillary condensation, driven by pressure fluctuations in the water outside the pore. The free energy difference between the two states depends linearly on the radius. The free energy landscape shows how a metastable liquid state gradually develops with increasing radius. For radii larger than ca. 0.55 nm it becomes the globally stable state and the vapor state vanishes. One dimensional confinement affects the dynamic behavior of the water molecules and increases the self diffusion by a factor of two to three compared to bulk water. Permeabilities for the narrow pores are of the same order of magnitude as for biological water pores. Water flow is not continuous but occurs in bursts. Our results suggest that simulations aimed at collective phenomena such as hydrophobic effects may require simulation times longer than 50 ns. For water in confined geometries, it is not possible to extrapolate from bulk or short time behavior to longer time scales.