Pathways for the formation of ice polymorphs from water predicted by a metadynamics method

Anisotropy in spinodal-like dynamics of unknown water at ice V-water interface

Experimentally demonstrating the existence of waters with local structures unlike that of common water is critical for understanding both the origin of the mysterious properties of water and liquid polymorphism in single component liquids. At the interfaces between water and ices Ih, III, and VI grown/melted under pressure, we previously discovered low- and high-density unknown waters, that are immiscible with the surrounding water. Here, we show, by in-situ optical microscopy, that an unknown water appears at the ice V-water interface via spinodal-like dynamics. The dewetting dynamics of the unknown water indicate that its characteristic velocity is ~ 90 m/s. The time evolution of the characteristic length of the spinodal-like undulation suggests that the dynamics may be described by a common model for spinodal decomposition of an immiscible liquid mixture. Spinodal-like dewetting dynamics of the unknown water transiently showed anisotropy, implying the property of a liquid crystal.

Collective Variables for Crystallization Simulations-from Early Developments to Recent Advances

Crystallization is an important physicochemical process which has relevance in material science, biology, and the environment. Decades of experimental and theoretical efforts have been made to understand this fundamental symmetry-breaking transition. While experiments provide equilibrium structures and shapes of crystals, they are limited to unraveling how molecules aggregate to form crystal nuclei that subsequently transform into bulk crystals. Computer simulations, mainly molecular dynamics (MD), can provide such microscopic details during the early stage of a crystallization event. Crystallization is a rare event that takes place in time scales much longer than a typical equilibrium MD simulation can sample. This inadequate sampling of the MD method can be easily circumvented by the use of enhanced sampling (ES) simulations. In most of the ES methods, the fluctuations of a system's slow degrees of freedom, called collective variables (CVs), are enhanced by applying a bias potential. This transforms the system from one state to the other within a short time scale. The most crucial part of such CV-based ES methods is to find suitable CVs, which often needs intuition and several trial-and-error optimization steps. Over the years, a plethora of CVs has been developed and applied in the study of crystallization. In this review, we provide a brief overview of CVs that have been developed and used in ES simulations to study crystallization from melt or solution. These CVs can be categorized mainly into four types: (i) spherical particle-based, (ii) molecular template-based, (iii) physical property-based, and (iv) CVs obtained from dimensionality reduction techniques. We present the context-based evolution of CVs, discuss the current challenges, and propose future directions to further develop effective CVs for the study of crystallization of complex systems.

Effect of nitrogen molecules on the growth kinetics at the interface between a (111) plane of cubic ice and water

The molecular-scale growth kinetics of ice from water in the presence of air molecules are still poorly understood, despite their importance for understanding ice particle formation in nature. In this study, a molecular dynamics simulation is conducted to elucidate the molecular-scale growth kinetics at the interface between a (111) plane of cubic ice and water in the presence of N2 molecules. Two potential models of N2 molecules with and without atomic charges are examined. For both models, N2 molecules bind stably to the interface for a period of 1 ns or longer, and the stability of the binding is higher for the charged model than for the noncharged model. Free-energy surfaces of an N2 molecule along the interface and along an ideal (111) plane surface of cubic ice suggest that for both models, the position where an N2 molecule binds stably is different at the interface and on the ideal plane surface, and the stability of the binding is much higher for the interface than for the ideal plane surface. For both models, stacking-disordered ice grows at the interface, and the formation probability of a hexagonal ice layer in the stacking-disordered ice is higher for the charged model than for the uncharged model. The formation probability for the hexagonal ice layer in the stacking-disordered ice depends not only on the stability of binding but also on the positions where N2 molecules bind on the underlying ice, and the number of N2 molecules that bind stably to the underlying ice.

Stable Binding Conformations of Polymaleic and Polyacrylic Acids at a Calcite Surface in the Presence of Countercations: A Metadynamics Study

Elucidating the stable binding conformations of additives at the surface of CaCO3 crystals is essential to biomineralization, scale inhibition, and materials technology. However, accomplishing this by experimental means is rather difficult. In this study, molecular dynamics simulations based on a metadynamics approach were conducted to elucidate the stable binding conformations of a deprotonated polymaleic acid (PMA) additive and two deprotonated poly(acrylic acid) (PAA) additives with different polymerization degrees in the presence of various countercations at a hydrated calcite (104) surface. The simulated free-energy surfaces suggested the existence of several slightly different stable binding conformations for each additive. The appearance of these distinct binding conformations is speculated to originate from different balances of interactions between the additive, the calcite surface, and the countercations. The binding conformations and binding stabilities at the calcite surface were affected by the countercations, with Ca2+ ions producing a more pronounced effect than Na+ ions. Furthermore, the simulation results suggested that the binding stability at the calcite surface was higher for the PMA additive than for the PAA additives, and the PAA additive with a polymerization degree of 10 displayed a binding stability that was similar to or lower than that of the PAA additive with a polymerization degree of 5. The present simulation method provides a new strategy for analyzing the binding conformations of complex additives at material surfaces, developing additives that stably bind to these surfaces, and designing additives to control crystal growth.

High-Density Liquid Water at a Water-Ice Interface

Because ice surfaces catalyze various key chemical reactions impacting nature and human life, the structure and dynamics of interfacial layers between water vapor and ice have been extensively debated with attention to the quasi-liquid layer. Other interfaces between liquid water and ice remain relatively underexplored, despite their importance and abundance on the Earth and icy extraterrestrial bodies. By in situ optical microscopy, we found that a high-density liquid layer, distinguishable from bulk water, formed at the interface between water and high-pressure ice III or VI, when they were grown or melted in a sapphire anvil cell. The liquid layer showed a bicontinuous pattern, indicating that immiscible waters with distinct structures were separated on the interfaces in a similar manner to liquid-liquid phase separation through spinodal decomposition. Our observations not only provide a novel opportunity to explore ice surfaces but also give insight into the two kinds of structured water.