Effect of Surface Polarity on the Structure and Dynamics of Water in Nanoscale Confinement

We present a molecular dynamics simulation study of the structure and dynamics of water confined between silica surfaces using beta-cristobalite as a model template. We scale the surface Coulombic charges by means of a dimensionless number, k, ranging from 0 to 1, and thereby we can model systems ranging frorn hydrophobic apolar to hydrophilic, respectively. Both rotational and translational dynamics exhibit a nonmonotonic dependence on k characterized by a maximum in the in-plane diffusion coefficient, D-parallel to, at values between 0.6 and 0.8, and a minimum in the rotational relaxation time, tau(R), at k = 0.6. The slow dynamics observed in the proximity of the hydrophobic apolar surface are a consequence of beta-cristobalite templating an ice-like water layer. The fully hydrophilic surfaces (k = 1.0), on the other hand, result in slow interfacial dynamics due to the presence of dense but disordered water that forms strong hydrogen bonds with surface silanol groups. Confinement also induces decoupling between translational and rotational dynamics, as evidenced by the fact that TR attains values similar to that of the bulk, while D-parallel to is always lower than in the bulk. The decoupling is characterized by a more drastic reduction in the translational dynamics of water compared to rotational relaxation.