Water Confined in Cylindrical Pores: A Molecular Dynamics Study

Molecular dynamics simulations of water confined in two hydrophilic cylindrical pores—PH and PL—that differ in their silanol surface concentration (7.6 and 3.0 nm⁻², respectively) have been performed at 300 K. A strong interaction of interfacial water molecules with the pore was systematically found and gives rise to a layering effect, a significant distortion of both the hydrogen bond network (HBN) and the tetrahedral structure of these water molecules, and a corresponding subdiffusive mean square displacement along the main axis of the pores. By contrast, water molecules in the inner part of both pores were found to behave similarly to bulk water. The HBN and the tetrahedral configuration of water were more gradually distorted in the PL pore given the larger heterogeneity and rugosity of the surface, and the number of water–pore hydrogen bonds did not scale linearly with the silanol surface concentration of the pores, in part because of the close proximity between silanols in the PH pore. With the PL pore, the dynamic slowing down of water was found consistent with the experiment, suggesting that it provides a better model for the cylindrical MCM-41 mesopores. The structural and dynamical properties of water vary little with the silica force field.