The nucleation and growth of clathrate hydrates of a hydrophobic guest comparable to methane or carbon dioxide are studied by molecular dynamics simulations of two-phase systems. The crystallization proceeds in two steps: First, the guest molecules concentrate in “blobs”, amorphous clusters involving multiple guest molecules in water-mediated configurations. These blobs are in dynamic equilibrium with the dilute solution and give birth to the clathrate cages that eventually transform it into an amorphous clathrate nucleus. In a second step, the amorphous clathrate transforms into crystalline clathrate. At low temperatures, the system can be arrested in the metastable amorphous clathrate phase for times sufficiently long for it to appear as an intermediate in the crystallization of clathrates. The “blob mechanism” unveiled in this work synthesizes elements of the labile cluster and local structuring hypotheses of clathrate nucleation and bears strong analogies to the two-step mechanisms of crystallization of proteins and colloids.

Clathrate hydrates and ice I are the most abundant crystals of water. The study of their nucleation, growth, and decomposition using molecular simulations requires an accurate and efficient algorithm that distinguishes water molecules that belong to each of these crystals and the liquid phase. Existing algorithms identify ice or clathrates, but not both. This poses a challenge for cases in which ice and hydrate coexist, such as in the synthesis of clathrates from ice and the formation of ice from clathrates during self-preservation of methane hydrates. Here we present an efficient algorithm for the identification of clathrate hydrates, hexagonal ice, cubic ice, and liquid water in molecular simulations. CHILL+ uses the number of staggered and eclipsed water–water bonds to identify water molecules in cubic ice, hexagonal ice, and clathrate hydrate. CHILL+ is an extension of CHILL (Moore et al. Phys. Chem. Chem. Phys. 2010, 12, 4124–4134), which identifies hexagonal and cubic ice but not clathrates. In addition to the identification of hydrates, CHILL+ significantly improves the detection of hexagonal ice up to its melting point. We validate the use of CHILL+ for the identification of stacking faults in ice and the nucleation and growth of clathrate hydrates. To our knowledge, this is the first algorithm that allows for the simultaneous identification of ice and clathrate hydrates, and it does so in a way that is competitive with respect to existing methods used to identify any of these crystals.

The nucleation, growth, structure and melting of ice in 3 nm diameter hydrophilic nanopores are studied through molecular dynamics simulations with the mW water model. The melting temperature of water in the pore was T^pore _m = 223 K, 51 K lower than the melting point of bulk water in the model and in excellent agreement with experimental determinations for 3 nm silica pores. Liquid and ice coexist in equilibrium at the melting point and down to temperatures as low as 180 K. Liquid water is located at the interface of the pore wall, increasing from one monolayer at the freezing temperature, T^pore _f = 195 K, to two monolayers a few degrees below T^pore _m. Crystallization of ice in the pore occurs through homogeneous nucleation. At the freezing temperature, the critical nucleus contains ∼75 to 100 molecules, with a radius of gyration similar to the radius of the pore. The critical nuclei contain features of both cubic and hexagonal ice, although stacking of hexagonal and cubic layers is not defined until the nuclei reach ∼150 molecules. The structure of the confined ice is rich in stacking faults, in agreement with the interpretation of X-ray and neutron diffraction experiments. Though the presence of cubic layers is twice as prevalent as hexagonal ones, the crystals should not be considered defective Ic as sequences with more than three adjacent cubic (or hexagonal) layers are extremely rare in the confined ice.

Ice crystallized below 200 K has the diffraction pattern of a faulty cubic ice, and not of the most stable hexagonal ice polymorph. The origin and structure of this faulty cubic ice, presumed to form in the atmosphere, has long been a puzzle. Here we use large-scale molecular dynamics simulations with the mW water model to investigate the crystallization of water at 180 K and elucidate the development of cubic and hexagonal features in ice as it nucleates, grows and consolidates into crystallites with characteristic dimensions of a few nanometres. The simulations indicate that the ice crystallized at 180 K contains layers of cubic ice and hexagonal ice in a ratio of approximately 2 to 1. The stacks of hexagonal ice are very short, mostly one and two layers, and their frequency does not seem to follow a regular pattern. In spite of the high fraction of hexagonal layers, the diffraction pattern of the crystals is, as in the experiments, almost identical to that of cubic ice. Stacking of cubic and hexagonal layers is observed for ice nuclei with as little as 200 water molecules, but a preference for cubic ice is already well developed in ice nuclei one order of magnitude smaller: the critical ice nuclei at 180 K contain approximately ten water molecules in their core and are already rich in cubic ice. The energies of the cubic-rich and hexagonal-rich nuclei are indistinguishable, suggesting that the enrichment in cubic ice does not have a thermodynamic origin.

Heterogeneous nucleation of ice induced by organic materials is of fundamental importance for climate, biology, and industry. Among organic ice-nucleating surfaces, monolayers of long chain alcohols are particularly effective, while monolayers of fatty acids are significantly less so. As these monolayers expose to water hydroxyl groups with an order that resembles the one in the basal plane of ice, it was proposed that lattice matching between ice and the surface controls their ice-nucleating efficiency. Organic monolayers are soft materials and display significant fluctuations. It has been conjectured that these fluctuations assist in the nucleation of ice. Here we use molecular dynamic simulations and laboratory experiments to investigate the relationship between the structure and fluctuations of hydroxylated organic surfaces and the temperature at which they nucleate ice. We find that these surfaces order interfacial water to form domains with ice-like order that are the birthplace of ice. Both mismatch and fluctuations decrease the size of the preordered domains and monotonously decrease the ice freezing temperature. The simulations indicate that fluctuations depress the freezing efficiency of monolayers of alcohols or acids to half the value predicted from lattice mismatch alone. The model captures the experimental trend in freezing efficiencies as a function of chain length and predicts that alcohols have higher freezing efficiency than acids of the same chain length. These trends are mostly controlled by the modulation of the structural mismatch to ice. We use classical nucleation theory to show that the freezing efficiencies of the monolayers are directly related to their free energy of binding to ice. This study provides a general framework to relate the equilibrium thermodynamics of ice binding to a surface and the nonequilibrium ice freezing temperature and suggests that these could be predicted from the structure of interfacial water.