Search

The structure of liquid water in the proximity of an interface can deviate significantly from that of bulk water, with surface-induced structural perturbations typically converging to bulk values at about ~1 nm from the interface. While these structural changes are well established it is, in contrast, less clear how an interface perturbs the dynamics of water molecules within the liquid. Here, through an extensive set of molecular dynamics simulations of supercooled bulk and interfacial water films and nano-droplets, we observe the formation of persistent, spatially extended dynamical domains in which the average mobility varies as a function of the distance from the interface. This is in stark contrast with the dynamical heterogeneity observed in bulk water, where these domains average out spatially over time. We also find that the dynamical response of water to an interface depends critically on the nature of the interface and on the choice of interface definition. Overall these results reveal a richness in the dynamics of interfacial water that opens up the prospect of tuning the dynamical response of water through specific modifications of the interface structure or confining material.

When an ice crystal is born from liquid water two key changes occur: (i) the molecules order; and (ii) the mobility of the molecules drops as they adopt their lattice positions. Most research on ice nucleation (and crystallization in general) has focused on understanding the former with less attention paid to the latter. However, supercooled water exhibits fascinating and complex dynamical behavior, most notably dynamical heterogeneity (DH), a phenomenon where spatially separated domains of relatively mobile and immobile particles coexist. Strikingly, the microscopic connection between the DH of water and the nucleation of ice has yet to be unraveled directly at the molecular level. Here we tackle this issue via computer simulations which reveal that: (i) ice nucleation occurs in low-mobility regions of the liquid; (ii) there is a dynamical incubation period in which the mobility of the molecules drops prior to any ice-like ordering; and (iii) ice-like clusters cause arrested dynamics in surrounding water molecules. With this we establish a clear connection between dynamics and nucleation. We anticipate that our findings will pave the way for the examination of the role of dynamical heterogeneities in heterogeneous and solution-based nucleation.

Molecular adsorption on surfaces plays a central role in catalysis, corrosion, desalination, and many other processes of relevance to industry and the natural world. Few adsorption systems are more ubiquitous or of more widespread importance than those involving water and carbon, and for a molecular level understanding of such interfaces water monomer adsorption on graphene is a fundamental and representative system. This system is particularly interesting as it calls for an accurate treatment of electron correlation effects, as well as posing a practical challenge to experiments. Here, we employ many-body electronic structure methodologies that can be rigorously converged and thus provide faithful references for the molecule-surface interaction. In particular, we use diffusion Monte-Carlo (DMC), coupled cluster (CCSD(T)), as well as the random phase approximation (RPA) to calculate the strength of the interaction between water and an extended graphene surface. We establish excellent, sub-chemical, agreement between the complementary high-level methodologies, and an adsorption energy estimate in the most stable configuration of approximately -100\,meV is obtained. We also find that the adsorption energy is rather insensitive to the orientation of the water molecule on the surface, despite different binding motifs involving qualitatively different interfacial charge reorganisation. In producing the first demonstrably accurate adsorption energies for water on graphene this work also resolves discrepancies amongst previously reported values for this widely studied system. It also paves the way for more accurate and reliable studies of liquid water at carbon interfaces with cheaper computational methods, such as density functional theory and classical potentials.

Ice nucleation plays a significant role in a large number of natural and technological processes, but it is challenging to investigate experimentally because of the small time (ns) and short length scales (nm) involved. On the other hand, conventional molecular simulations struggle to cope with the relatively long timescale required for critical ice nuclei to form. One way to tackle this issue is to take advantage of free energy or path sampling techniques. Unfortunately, these are computationally costly. Seeded molecular dynamics is a much less demanding alternative that has been successfully applied already to study the homogeneous freezing of water. However, in the case of heterogeneous ice nucleation, nature's favourite route to form ice, an array of suitable interfaces between the ice seeds and the substrate of interest has to be built - and this is no trivial task. In this paper, we present a Heterogeneous SEEDing approach (HSEED) which harnesses a random structure search framework to tackle the ice-substrate challenge, thus enabling seeded molecular dynamics simulations of heterogeneous ice nucleation on crystalline surfaces. We validate the HSEED framework by investigating the nucleation of ice on: (i) model crystalline surfaces, using the coarse-grained mW model; and (ii) cholesterol crystals, employing the fully atomistic TIP4P/Ice water model. We show that the HSEED technique yields results in excellent agreement with both metadynamics and forward flux sampling simulations. Because of its computational efficiency, the HSEED method allows one to rapidly assess the ice nucleation ability of whole libraries of crystalline substrates - a long-awaited computational development in e.g. atmospheric science., Comment: 13 pages, 7 figures

Non-bonded potentials are included in most force fields and therefore widely used in classical molecular dynamics simulations of materials and interfacial phenomena. It is commonplace to truncate these potentials for computational efficiency based on the assumption that errors are negligible for reasonable cutoffs or compensated for by adjusting other interaction parameters. Arising from a metadynamics study of the wetting transition of water on a solid substrate, we find that the influence of the cutoff is unexpectedly strong and can change the character of the wetting transition from continuous to first order by creating artificial metastable wetting states. Common cutoff corrections such as the use of a force switching function, a shifted potential, or a shifted force do not avoid this. Such a qualitative difference urges caution and suggests that using truncated non-bonded potentials can induce unphysical behavior that cannot be fully accounted for by adjusting other interaction parameters.

Freezing of water is arguably one of the most common phase transitions on Earth and almost always happens heterogeneously. Despite its importance, we lack a fundamental understanding of what makes substrates efficient ice nucleators. Here we address this by computing the ice nucleation (IN) ability of numerous model hydroxylated substrates with diverse surface hydroxyl (OH) group arrangements. Overall, for the substrates considered, we find that neither the symmetry of the OH patterns nor the similarity between a substrate and ice correlate well with the IN ability. Instead, we find that the OH density and the substrate-water interaction strength are useful descriptors of a material's IN ability. This insight allows the rationalization of ice nucleation ability across a wide range of materials, and can aid the search and design of novel potent ice nucleators in the future., Comment: main text + SI

What makes a material a good ice nucleating agent? Despite the importance of heterogeneous ice nucleation to a variety of fields, from cloud science to microbiology, major gaps in our understanding of this ubiquitous process still prevent us from answering this question. In this work, we have examined the ability of generic crystalline substrates to promote ice nucleation as a function of the hydrophobicity and the morphology of the surface. Nucleation rates have been obtained by brute-force molecular dynamics simulations of coarse-grained water on top of different surfaces of a model fcc crystal, varying the water-surface interaction and the surface lattice parameter. It turns out that the lattice mismatch of the surface with respect to ice, customarily regarded as the most important requirement for a good ice nucleating agent, is at most desirable but not a requirement. On the other hand, the balance between the morphology of the surface and its hydrophobicity can significantly alter the ice nucleation rate and can also lead to the formation of up to three different faces of ice on the same substrate. We have pinpointed three circumstances where heterogeneous ice nucleation can be promoted by the crystalline surface: (i) the formation of a water overlayer that acts as an in-plane template; (ii) the emergence of a contact layer buckled in an ice-like manner; and (iii) nucleation on compact surfaces with very high interaction strength. We hope that this extensive systematic study will foster future experimental work aimed at testing the physiochemical understanding presented herein., Comment: Main + SI

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insight into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that in the last few decades have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state of the art computational methods, by reviewing simulations of e.g. ice nucleation or crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insight into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that in doing so the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that by improving (i.) existing interatomic potentials; and (ii.) currently available enhanced sampling methods, the community can move towards accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.