This document provides a starting point for approaching molecular simulations, guiding beginning practitioners to what issues they need to know about before and while starting their first simulations, and why those issues are so critical. This document makes no claims to provide an adequate introduction to the subject on its own. Instead, our goal is to help people know what issues are critical before beginning, and to provide references to good resources on those topics. We also provide a checklist of key issues to consider before and while setting up molecular simulations which may serve as a foundation for other best practices documents.

An in-depth knowledge of the interaction of water with amorphous silica is critical to fundamental studies of interfacial hydration water, as well as to industrial processes such as catalysis, nanofabrication, and chromatography. Silica has a tunable surface comprising hydrophilic silanol groups and moderately hydrophobic siloxane groups that can be interchanged through thermal and chemical treatments. Despite extensive studies of silica surfaces, the influence of surface hydrophilicity and chemical topology on the molecular properties of interfacial water is not well understood. In this work, we controllably altered the surface silanol density, and measured surface water diffusivity using Overhauser dynamic nuclear polarization (ODNP) and complementary silica-silica interaction forces across water using a surface forces apparatus (SFA). The results show that increased silanol density generally leads to slower water diffusivity and stronger silica-silica repulsion at short aqueous separations (less than ∼4 nm). Both techniques show sharp changes in hydration properties at intermediate silanol densities (2.0-2.9 nm). Molecular dynamics simulations of model silica-water interfaces corroborate the increase in water diffusivity with silanol density, and furthermore show that even on a smooth and crystalline surface at a fixed silanol density, adjusting the spatial distribution of silanols results in a range of surface water diffusivities spanning ∼10%. We speculate that a critical silanol cluster size or connectivity parameter could explain the sharp transition in our results, and can modulate wettability, colloidal interactions, and surface reactions, and thus is a phenomenon worth further investigation on silica and chemically heterogeneous surfaces.

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

The interactions of water with solid surfaces govern their apparent hydrophobicity/hydrophilicity, influenced at the molecular scale by surface coverage of chemical groups of varied nonpolar/polar character. Recently, it has become clear that the precise patterning of surface groups, and not simply average surface coverage, has a significant impact on the structure and thermodynamics of hydration layer water, and, in turn, on macroscopic interfacial properties. Here we show that patterning also controls the dynamics of hydration water, a behavior frequently thought to be leveraged by biomolecules to influence functional dynamics, but yet to be generalized. To uncover the role of surface heterogeneities, we couple a genetic algorithm to iterative molecular dynamics simulations to design the patterning of surface functional groups, at fixed coverage, to either minimize or maximize proximal water diffusivity. Optimized surface configurations reveal that clustering of hydrophilic groups increases hydration water mobility, while dispersing them decreases it, but only if hydrophilic moieties interact with water through directional, hydrogen-bonding interactions. Remarkably, we find that, across different surfaces, coverages, and patterns, both the chemical potential for inserting a methane-sized hydrophobe near the interface and, in particular, the hydration water orientational entropy serve as strong predictors for hydration water diffusivity, suggesting that these simple thermodynamic quantities encode the way surfaces control water dynamics. These results suggest a deep and intriguing connection between hydration water thermodynamics and dynamics, demonstrating that subnanometer chemical surface patterning is an important design parameter for engineering solid-water interfaces with applications spanning separations to catalysis.

A tetrahedral structure resulting from hydrogen bonding is a hallmark of liquid water and plays a significant role in determining its unique thermophysical properties. This water feature has helped understand anomalous properties and physically interpret and model hydrophobic solvation thermodynamics. Tetrahedrality is well described by the geometric relationship of any central water molecule with two of its nearest neighbors in the first coordination shell, as defined by the corresponding “three-body” angle. While order parameters and even full water models have been developed using specific or average features of the three-body angle distribution, here we examine the distribution holistically, tracking its response to changes in temperature, density, and the presence of model solutes. Surprisingly, we find that the three-body distribution responds by varying primarily along a single degree of freedom, suggesting a remarkably simplified view of water structure. We characterize three-body angle distributions across temperature and density space and identify principal components of the variations with state conditions. We show that these principal components embed physical significance and trace out transitions between tetrahedral and simple-fluid-like behavior. Moreover, we find that the ways three-body angles vary within the hydration shells of model colloids of different types and sizes are nearly identical to the variations seen in bulk water across density and temperature. Importantly, through the principal directions of these variations, we find that perturbations to the hydration-water distributions well predict the thermodynamics associated with colloid solvation, in particular, the relative entropy of this process that captures indirect, solvent-mediated contributions to the hydration free energy.