Interfacial water on crystalline silica: a comparative molecular dynamics simulation study

Understanding the properties of interfacial water at solid-liquid interfaces is important in a wide range of applications. Molecular dynamics is becoming a widespread tool for this purpose. Unfortunately, however, the results of such studies are known to strongly depend on the selection of force fields. It is, therefore, of interest to assess the extent by which the implemented force fields can affect the predicted properties of interfacial water. Two silica surfaces, with low and high surface hydroxyl density, respectively, were simulated implementing four force fields. These force fields yield different orientation and flexibility of surface hydrogen atoms, and also different interaction potentials with water molecules. The properties for interfacial water were quantified by calculating contact angles, atomic density profiles, surface density distributions, hydrogen bond density profiles and residence times for water near the solid substrates. We found that at low surface density of hydroxyl groups, the force field strongly affects the predicted contact angle, while at high density of hydroxyl groups, water wets all surfaces considered. From a molecular-level point of view, our results show that the position and intensity of peaks observed from oxygen and hydrogen atomic density profiles are quite different when different force fields are implemented, even when the simulated contact angles are similar. Particularly, the surfaces simulated by the CLAYFF force field appear to attract water more strongly than those simulated by the Brodka and Zerda force field. It was found that the surface density distributions for water strongly depend on the orientation of surface hydrogen atoms. In all cases, we found an elevated number of hydrogen bonds formed between interfacial water molecules. The hydrogen bond density profile does not depend strongly on the force field implemented to simulate the substrate, suggesting that interfacial water assumes the necessary orientation to maximise the number of water-water hydrogen bonds irrespectively of surface properties. Conversely, the residence time for water molecules near the interface strongly depends on the force field and on the flexibility of surface hydroxyl groups. Specifically, water molecules reside for longer times at contact with rigid substrates with high density of hydroxyl groups. These results should be considered when comparisons between simulated and experimental data are attempted.