Upon cooling, water freezes to ice. This familiar phase transition occurs widely in nature, yet unlike the freezing of simple liquids, it has never been successfully simulated on a computer. The difficulty lies with the fact that hydrogen bonding between individual water molecules yields a disordered three-dimensional hydrogen-bond network whose rugged and complex global potential energy surface permits a large number of possible network configurations. As a result, it is very challenging to reproduce the freezing of 'real' water into a solid with a unique crystalline structure. For systems with a limited number of possible disordered hydrogen-bond network structures, such as confined water, it is relatively easy to locate a pathway from a liquid state to a crystalline structure. For pure and spatially unconfined water, however, molecular dynamics simulations of freezing are severely hampered by the large number of possible network configurations that exist. Here we present a molecular dynamics trajectory that captures the molecular processes involved in the freezing of pure water. We find that ice nucleation occurs once a sufficient number of relatively long-lived hydrogen bonds develop spontaneously at the same location to form a fairly compact initial nucleus. The initial nucleus then slowly changes shape and size until it reaches a stage that allows rapid expansion, resulting in crystallization of the entire system.

The adsorption of guest and kinetic inhibitor molecules on the surface of methane hydrate is investigated by using molecular dynamics simulations. We calculate the free energy profile for transferring a solute molecule from bulk water to the hydrate surface for various molecules. Spherical solutes with a diameter of ∼0.5 nm are significantly stabilized at the hydrate surface, whereas smaller and larger solutes exhibit lower adsorption affinity than the solutes of intermediate size. The range of the attractive force is subnanoscale, implying that this force has no effect on the macroscopic mass transfer of guest molecules in crystal growth processes of gas hydrates. We also examine the adsorption mechanism of a kinetic hydrate inhibitor. It is found that a monomer of the kinetic hydrate inhibitor is strongly adsorbed on the hydrate surface. However, the hydrogen bonding between the amide group of the inhibitor and water molecules on the hydrate surface, which was believed to be the driving force for the adsorption, makes no contribution to the adsorption affinity. The preferential adsorption of both the kinetic inhibitor and the spherical molecules to the surface is mainly due to the entropic stabilization arising from the presence of cavities at the hydrate surface. The dependence of surface affinity on the size of adsorbed molecules is also explained by this mechanism.

We report a molecular dynamics simulation demonstrating a fast spontaneous liquid-liquid phase separation of water and a subsequent slow crystallization to ice. It is found that supercooled water separates rapidly into low- and high-density domains so as to reduce the surface energy in the rectangular simulation cell at certain thermodynamic states. The liquid-liquid phase separation, which is about two orders of magnitude faster than the crystallization, suggests a possibility to observe this phenomenon experimentally.

GenIce is an efficient and user‐friendly tool to generate hydrogen‐disordered ice structures. It makes ice and clathrate hydrate structures in various file formats. More than 100 kinds of structures are preset. Users can install their own crystal structures, guest molecules, and file formats as plugins. The algorithm certifies that the generated structures are completely randomized hydrogen‐disordered networks obeying the ice rule with zero net polarization. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Recent reports indicate that the crystallization of clathrate hydrates occurs in multiple steps that involve amorphous intermediates and metastable clathrate crystals. The elucidation of the reaction coordinate for clathrate crystallization requires the use of order parameters able to identify the reactants, products, and intermediates in the crystallization pathway. Nevertheless, existing order parameters cannot distinguish between amorphous and crystalline clathrates or between different clathrate crystals. In this work, we present the first set of order parameters that discern between the sI and sII clathrate crystals, the amorphous clathrates, the blob of solvent-separated guests and the liquid solution. These order parameters can be used to monitor the advance of the crystallization and for the efficient implementation of methods to sample the rare clathrate nucleation events in molecular simulations. We illustrate the use of these order parameters in the analysis of the growth and the dissolution of clathrate crystals and the spontaneous nucleation and growth of clathrates under conditions of high supercooling.