Ice formation on kaolinite: Insights from molecular dynamics simulations

Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation

Heterogeneous nucleation is the main path to ice formation on Earth. The ice nucleating ability of a certain substrate is mainly determined by both molecular interactions and the structural mismatch between the ice and the substrate lattices. We focus on the latter factor using molecular simulations of the mW model. Quantifying the effect of structural mismatch alone is challenging due to its coupling with molecular interactions. To disentangle both the factors, we use a substrate composed of water molecules in such a way that any variation on the nucleation temperature can be exclusively ascribed to the structural mismatch. We find that a 1% increase in structural mismatch leads to a decrease of ∼4 K in the nucleation temperature. We also analyze the effect of orientation of the substrate with respect to the liquid. The three main ice orientations (basal, primary prism, and secondary prism) have a similar ice nucleating ability. We finally assess the effect of lattice flexibility by comparing substrates where molecules are immobile to others where a certain freedom to fluctuate around the lattice positions is allowed. Interestingly, we find that the latter type of substrate is more efficient in nucleating ice because it can adapt its structure to that of ice.

The molecular scale mechanism of deposition ice nucleation on silver iodide

Heterogeneous ice nucleation is a ubiquitous process in the natural and built environment. Deposition ice nucleation, i.e. heterogeneous ice nucleation that - according to the traditional view - occurs in a subsaturated water vapor environment and in the absence of supercooled water on the solid, ice-forming surface, is among the most important ice formation processes in high-altitude cirrus and mixed-phase clouds. Despite its importance, very little is known about the mechanism of deposition ice nucleation at the microscopic level. This study puts forward an adsorption-based mechanism for deposition ice nucleation through results from a combination of atomistic simulations, experiments and theoretical modelling. One of the most potent laboratory surrogates of ice nucleating particles, silver iodide, is used as a substrate for the simulations. We find that water initially adsorbs in clusters which merge and grow over time to form layers of supercooled water. Ice nucleation on silver iodide requires at minimum the adsorption of 4 molecular layers of water. Guided by the simulations we propose the following fundamental freezing steps: (1) Water molecules adsorb on the surface, forming nanodroplets. (2) The supercooled water nanodroplets merge into a continuous multilayer when they grow to about 3 molecular layers thick. (3) The layer continues to grow until the critical thickness for freezing is reached. (4) The critical ice cluster continues to grow.

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

Molecular simulation of the confined crystallization of ice in cement nanopore

Freezing of water under nanoconfinement exhibits physical peculiarities with respect to the bulk water. However, experimental observations are extremely challenging at this scale, which limits our understanding of the effect of confinement on water properties upon freezing. In this study, we use molecular dynamic simulations to investigate how confinement affects the kinetics of growth of ice and the thermodynamic equilibrium of ice-liquid coexistence. TIP4P/Ice water model and CSH-FF model were applied to simulate ice crystallization in a confined cement system at temperatures down to 220 K. We adapted an interface detection algorithm and reparameterized the CHILL/CHILL+ algorithm to capture ice growth. The confinement leads to a shift of the maximum growth rate of ice to a higher temperature than for bulk water. Both the confinement and surface impurities contribute to slowing down the ice growth. For the ice-liquid coexistence at equilibrium, we derive a formulation of Thomson's equation adapted to statistical physics quantities accessible by molecular simulation, and we show that this adapted equation predicts accurately the melting line of bulk and confined ice Ih as a function of pressure. The confinement decreases systematically the melting temperature of ice of about 5 K compared with bulk ice Ih. A premelted water film about 1 nm thick is observed between the solid wall and ice, and its thickness is found to decrease continuously as temperature is lowered. We note that the surface impurities are key to the formation of the premelted water nanofilm when the temperature is lower than 250 K.

A fully quantum-mechanical treatment for kaolinite

Neural network potentials for kaolinite minerals have been fitted to data extracted from density functional theory calculations that were performed using the revPBE + D3 and revPBE + vdW functionals. These potentials have then been used to calculate the static and dynamic properties of the mineral. We show that revPBE + vdW is better at reproducing the static properties. However, revPBE + D3 does a better job of reproducing the experimental IR spectrum. We also consider what happens to these properties when a fully quantum treatment of the nuclei is employed. We find that nuclear quantum effects (NQEs) do not make a substantial difference to the static properties. However, when NQEs are included, the dynamic properties of the material change substantially.

RSeeds: Rigid Seeding Method for Studying Heterogeneous Crystal Nucleation

Heterogeneous nucleation is the dominant form of liquid-to-solid transition in nature. Although molecular simulations are most uniquely suited to studying nucleation, the waiting time to observe even a single nucleation event can easily exceed the current computational capabilities. Therefore, there exists an imminent need for methods that enable computationally fast and feasible studies of heterogeneous nucleation. Seeding is a technique that has proven to be successful at dramatically expanding the range of computationally accessible nucleation rates in simulation studies of homogeneous crystal nucleation. In this article, we introduce a new seeding method for heterogeneous nucleation called Rigid Seeding (RSeeds). Crystalline seeds are treated as pseudorigid bodies and simulated on a surface with metastable liquid above its melting temperature. This allows the seeds to adapt to the surface and identify favorable seed-surface configurations, which is necessary for reliable predictions of crystal polymorphs that form and the corresponding heterogeneous nucleation rates. We demonstrate and validate RSeeds for heterogeneous ice nucleation on a flexible self-assembled monolayer surface, a mineral surface based on kaolinite, and two model surfaces. RSeeds predicts the correct ice polymorph, exposed crystal plane, and rotation on the surface. RSeeds is semiquantitative and can be used to estimate the critical nucleus size and nucleation rate when combined with classical nucleation theory. We demonstrate that RSeeds can be used to evaluate nucleation rates spanning many orders of magnitude.

Accurate prediction of ice nucleation from room temperature water

Crystal nucleation is one of the most fundamental processes in the physical sciences and almost always occurs heterogeneously with the aid of a nucleating substrate. No example of nucleation is more ubiquitous and impactful than the formation of ice, vital to fields as diverse as geology, biology, aeronautics, and climate science. However, despite considerable effort, we still cannot predict a priori the efficacy of a nucleating agent. Here we utilize deep learning methods to accurately predict nucleation ability from images of room temperature liquid water-generated from molecular dynamics simulations-on a broad range of substrates. The resulting model, named IcePic, can rapidly and accurately infer nucleation ability, eliminating the requirement for either notoriously expensive simulations or direct experimental measurement. In an online poll, IcePic was found to significantly outperform humans in predicting the ice nucleating efficacy of materials. By analyzing the typical errors made by humans, as well as the application of reverse interpretation methods, physical insights into the role the water contact layer plays in ice nucleation have been obtained. Moving forward, we suggest that IcePic can be used as an easy, cheap, and rapid way to discern the nucleation ability of substrates, also with potential for learning other properties related to interfacial water.

The role of structural order in heterogeneous ice nucleation

The freezing of water into ice is a key process that is still not fully understood. It generally requires an impurity of some description to initiate the heterogeneous nucleation of the ice crystals. The molecular structure, as well as the extent of structural order within the impurity in question, both play an essential role in determining its effectiveness. However, disentangling these two contributions is a challenge for both experiments and simulations. In this work, we have systematically investigated the ice-nucleating ability of the very same compound, cholesterol, from the crystalline (and thus ordered) form to disordered self-assembled monolayers. Leveraging a combination of experiments and simulations, we identify a “sweet spot” in terms of the surface coverage of the monolayers, whereby cholesterol maximises its ability to nucleate ice (which remains inferior to that of crystalline cholesterol) by enhancing the structural order of the interfacial water molecules. These findings have practical implications for the rational design of synthetic ice-nucleating agents.

The freezing of water into ice is still not fully understood. Here, we investigate the role of structural disorder within the biologically relevant impurities that facilitate this fundamental phase transition.

Can molecular simulations reliably compare homogeneous and heterogeneous ice nucleation?

In principle, the answer to the posed titular question is undoubtedly "yes." But in practice, requisite reference data for homogeneous systems have been obtained with a treatment of intermolecular interactions that is different from that typically employed for heterogeneous systems. In this article, we assess the impact of the choice of truncation scheme when comparing water in homogeneous and inhomogeneous environments. Specifically, we use explicit free energy calculations and a simple mean field analysis to demonstrate that using the "cut-and-shift" version of the Lennard-Jones potential (common to most simple point charge models of water) results in a systematic increase in the melting temperature of ice I_h. In addition, by drawing an analogy between a change in cutoff and a change in pressure, we use existing literature data for homogeneous ice nucleation at negative pressures to suggest that enhancements due to heterogeneous nucleation may have been overestimated by several orders of magnitude.

Lipid bilayers as potential ice nucleating agents

Cellular damage is a key issue in the context of cryopreservation. Much of this damage is believed to be caused by extracellular ice formation at temperatures well above the homogeneous freezing point of pure water. Hence the question: what initiates ice nucleation during cryopreservation? In this paper, we assess whether cellular membranes could be responsible for facilitating the ice nucleation process, and what characteristics would make them good or bad ice nucleating agents. By means of molecular dynamics simulations, we investigate a number of phospholipids and lipopolysaccharide bilayers at the interface with supercooled liquid water. While these systems certainly appear to act as ice nucleating agents, it is likely that other impurities might also play a role in initiating extracellular ice nucleation. Furthermore, we elucidate the factors which affect a bilayer's ability to act as an ice nucleating agent; these are complex, with specific reference to both chemical and structural factors. These findings represent a first attempt to pinpoint the origin of extracellular ice nucleation, with important implications for the cryopreservation process.

Volcanic ash ice nucleation activity is variably reduced by aging in water and sulfuric acid: the effects of leaching, dissolution, and precipitation

Volcanic ash nucleates ice when immersed in supercooled water droplets, giving it the potential to influence weather and climate from local to global scales. This ice nucleation activity (INA) is likely derived from a subset of the crystalline mineral phases in the ash. The INA of other mineral-based dusts can change when exposed to various gaseous and aqueous chemical species, many of which also interact with volcanic ash in the eruption plume and atmosphere. However, the effects of aqueous chemical aging on the INA of volcanic ash have not been explored. We show that the INA of two mineralogically distinct ash samples from Fuego and Astroni volcanoes is variably reduced following immersion in water or aqueous sulfuric acid for minutes to days. Aging in water decreases the INA of both ash samples by up to two orders of magnitude, possibly due to a reduction in surface crystallinity and cation availability accompanying leaching. Aging in sulfuric acid leads to minimal loss of INA for Fuego ash, which is proposed to reflect a quasi-equilibrium between leaching that removes ice-active sites and dissolution that reveals or creates new sites on the pyroxene phases present. Conversely, exposure to sulfuric acid reduces the INA of Astroni ash by one to two orders of magnitude, potentially through selective dissolution of ice-active sites associated with surface microtextures on some K-feldspar phases. Analysis of dissolved element concentrations in the aged ash leachates shows supersaturation of certain mineral species which could have precipitated and altered the INA of the ash. These results highlight the key role that leaching, dissolution, and precipitation likely play in the aqueous aging of volcanic ash with respect to its INA. Finally, we discuss the implications for understanding the nature and reactivity of ice-active sites on volcanic ash and its role in influencing cloud properties in the atmosphere.

Effect of Graphene on Ice Polymorph

Recently, ice with stacking disorder structure, consisting of random sequences of cubic ice (Ic) and hexagonal ice (Ih) layers, was reported to be more stable than pure Ih/Ic. Due to a much lower free energy barrier of heterogeneous nucleation, in practice, the freezing process of water is controlled by heterogeneous nucleation triggered by an external medium. Therefore, we carry out molecular dynamic simulations to explore how ice polymorphism depends on the lattice structure of the crystalline substrates on which the ice is grown, focusing on the primary source of atmospheric aerosols, carbon materials. It turns out that, during the nucleation stage, the polymorph of ice nuclei is strongly affected by graphene substrates. For ice nucleation on graphene, we find Ih is the dominant polymorph. This can be attributed to structural similarities between graphene and basal face of Ih. Our results also suggest that the substrate only affects the polymorph of ice close to the graphene surface, with the preference for Ih diminishing as the ice layer grows.

Effect of interfacial dipole on heterogeneous ice nucleation

In this work, we performed molecular dynamics simulations of ice nucleation on a rigid surface model of cubic zinc blende structure with different surface dipole strength and orientation. Our results show that, although substrates are excellently lattice-matched to cubic ice, ice nucleation merely occurred as the interfacial water molecules (IWs) show identical or similar orientations to that of water molecules in cubic ice. Free energy landscapes revealed that, as substrates have non-suitable dipole strength/orientation, there exist large free energy barriers for rotating dipole IWs to the right orientation to trigger ice formation. This study stresses that, beyond the traditional view of lattice match and the similarity of lattice length between the substrate and new-formed crystal, the similarity between molecular orientations of interfacial component and component in the specific new-formed crystalline face is also critical for promoting ice nucleation.

Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation

Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.

Routes to cubic ice through heterogeneous nucleation

The freezing of water into ice is one of the most important processes in the physical sciences. However, it is still not understood at the molecular level. In particular, the crystallization of cubic ice ([Formula: see text])-rather than the traditional hexagonal polytype ([Formula: see text])-has become an increasingly debated topic. Although evidence for [Formula: see text] is thought to date back almost 400 y, it is only in the last year that pure [Formula: see text] has been made in the laboratory, and these processes involved high-pressure ice phases. Since this demonstrates that pure [Formula: see text] can form, the question naturally arises if [Formula: see text] can be made from liquid water. With this in mind, we have performed a high-throughput computational screening study involving molecular dynamics simulations of nucleation on over 1,100 model substrates. From these simulations, we find that 1) many different substrates can promote the formation of pristine [Formula: see text]; 2) [Formula: see text] can be selectively nucleated for even the mildest supercooling; 3) the water contact layer's resemblance to a face of ice is the key factor determining the polytype selectivity and nucleation temperature, independent of which polytype is promoted; and 4) substrate lattice match to ice is not indicative of the polytype obtained. Through this study, we have deepened understanding of the interplay of heterogeneous nucleation and ice I polytypism and suggest routes to [Formula: see text] More broadly, the substrate design methodology presented here combined with the insight gained can be used to understand and control polymorphism and stacking disorder in materials in general.

Macroscopic surface charges from microscopic simulations

Attaining accurate average structural properties in a molecular simulation should be considered a prerequisite if one aims to elicit meaningful insights into a system's behavior. For charged surfaces in contact with an electrolyte solution, an obvious example is the density profile of ions along the direction normal to the surface. Here, we demonstrate that, in the slab geometry typically used in simulations, imposing an electric displacement field D determines the integrated surface charge density of adsorbed ions at charged interfaces. This allows us to obtain macroscopic surface charge densities irrespective of the slab thickness used in our simulations. We also show that the commonly used Yeh-Berkowitz method and the "mirrored slab" geometry both impose vanishing integrated surface charge densities. We present results both for relatively simple rocksalt (1 1 1) interfaces and the more complex case of kaolinite's basal faces in contact with an aqueous electrolyte solution.

Predicting heterogeneous ice nucleation with a data-driven approach

Water in nature predominantly freezes with the help of foreign materials through a process known as heterogeneous ice nucleation. Although this effect was exploited more than seven decades ago in Vonnegut's pioneering cloud seeding experiments, it remains unclear what makes a material a good ice former. Here, we show through a machine learning analysis of nucleation simulations on a database of diverse model substrates that a set of physical descriptors for heterogeneous ice nucleation can be identified. Our results reveal that, beyond Vonnegut's connection with the lattice match to ice, three new microscopic factors help to predict the ice nucleating ability. These are: local ordering induced in liquid water, density reduction of liquid water near the surface and corrugation of the adsorption energy landscape felt by water. With this we take a step towards quantitative understanding of heterogeneous ice nucleation and the in silico design of materials to control ice formation.

Effects of Inorganic Ions on Ice Nucleation by the Al Surface of Kaolinite Immersed in Water

Molecular dynamics simulations are employed to investigate the influence of inorganic salts on ice nucleation by the Al surface of kaolinite, terminated with hydroxyl groups. Seven salt solutions (LiI(Cl), NaI(Cl), KI(Cl), and NH₄I) are considered. Simulations were performed at 300 K to obtain equilibrium surface-ion and surface-water density profiles. These simulations show no specific ion adsorption at the kaolinite surface. There are weak surface-ion correlations, with cations preferring to be closer to the surface than the anions. At a supercooling of 26 K (taking account of freezing point depression), 1 M salt solutions slowed ice nucleation by a factor of 2-3 compared with pure water and significantly reduced the rate of ice growth after nucleation. All salt solutions had similar influences on ice nucleation, and no specific ion effects were identified. Ice nucleation simulations for 1 M NaI(Cl), KI(Cl), and LiI solutions were performed for a range of temperatures. In all cases, the supercooling required for ice nucleation was larger by ∼1-6 K, after accounting for freezing point depression, than that required for pure water. For 1 M LiI solution an earlier laboratory study using kaolin as ice nucleating particles (INP) reported that the supercooling required for ice nucleation was ∼11 K smaller than that required for pure water. Our simulation results are not consistent with this finding. In this paper, we report new laboratory results for 1 M LiI solution employing kaolinite as INP. In our experiments ice nucleation in the LiI solution required the same supercooling as pure water, which is more consistent with our simulations.

Insights into the Emerging Networks of Voids in Simulated Supercooled Water

The structural evolution of supercooled liquid water as we approach the glass transition temperature continues to be an active area of research. Here, we use molecular dynamics simulations of TIP4P/ice water to study the changes in the connected regions of empty space within the liquid, which we investigate using the Voronoi-voids network. We observe two important features: supercooling enhances the fraction of nonspherical voids and different sizes of voids tend to cluster forming a percolating network. By examining order parameters such as the local structure index (LSI), tetrahedrality and topological defects, we show that water molecules near large void clusters tend to be slightly more tetrahedral than those near small voids, with a lower population of under- and overcoordinated defects. We show further that the distribution of closed rings of water molecules around small and large void clusters maintain a balance between 6 and 7 membered rings. Our results highlight the changes of the dual voids and water network as a structural hallmark of supercooling and provide insights into the molecular origins of cooperative effects underlying density fluctuations on the subnanometer and nanometer length scale. In addition, the percolation of the voids and the hydrogen bond network around the voids may serve as useful order parameters to investigate density fluctuations in supercooled water.

Seeding approach to nucleation in the NVT ensemble: The case of bubble cavitation in overstretched Lennard Jones fluids

Simulations are widely used to study nucleation in first order phase transitions due to the fact that they have access to the relevant length and time scales. However, simulations face the problem that nucleation is an activated process. Therefore, rare event simulation techniques are needed to promote the formation of the critical nucleus. The Seeding method, where the simulations are started with the nucleus already formed, has proven quite useful in efficiently providing estimates of the nucleation rate for a wide range of orders of magnitude. So far, Seeding has been employed in the NPT ensemble, where the nucleus either grows or redissolves. Thus, several trajectories have to be run in order to find the thermodynamic conditions that make the seeded nucleus critical. Moreover, the nucleus lifetime is short and the statistics for obtaining its properties is consequently poor. To deal with these shortcomings we extend the Seeding method to the NVT ensemble. We focus on the problem of bubble nucleation in a metastable Lennard Jones fluid. We show that, in the NVT ensemble, it is possible to equilibrate and stabilise critical bubbles for a long time. The nucleation rate inferred from NVT-Seeding is fully consistent with that coming from NPT-Seeding. The former is quite suitable to obtain the nucleation rate along isotherms, whereas the latter is preferable if the dependence of the rate with temperature at constant pressure is required. Care should be taken with finite size effects when using NVT-Seeding. Further work is required to extend NVT seeding to other sorts of phase transitions.

Heterogeneous ice nucleation correlates with bulk-like interfacial water

Establishing a direct correlation between interfacial water and heterogeneous ice nucleation (HIN) is essential for understanding the mechanism of ice nucleation. Here, we study the HIN efficiency on polyvinyl alcohol (PVA) surfaces with different densities of hydroxyl groups. We find that the HIN efficiency increases with the decreasing hydroxyl group density. By explicitly considering that interfacial water molecules of PVA films consist of "tightly bound water," "bound water," and "bulk-like water," we reveal that bulk-like water can be correlated directly to the HIN efficiency of surfaces. As the density of hydroxyl groups decreases, bulk-like water molecules can rearrange themselves with a reduced energy barrier into ice due to the diminishing constraint by the hydroxyl groups on the PVA surface. Our study not only provides a new strategy for experimentally controlling the HIN efficiency but also gives another perspective in understanding the mechanism of ice nucleation.

The Effect of Crystallinity and Crystal Structure on the Immersion Freezing of Alumina

Determining the factors that constitute an efficient ice nucleus is an ongoing area of research in the atmospheric community. In particular, surface characteristics such as functional groups and surface defects impact the ice nucleation efficiency. Crystal structure has been proposed to be a possible factor that can dictate ice nucleation activity through the templating of water molecules on the surface of the aerosol particle. If the crystal structure of the surface matches that of the crystal structure of ice, it has been shown to increase ice nucleation activity. In this study, alumina was chosen as a model system because crystal structure and crystallinity can be tuned, and the effect on immersion freezing was explored. The nine alumina samples include polymorphs of AlOOH, Al(OH)₃, and Al₂O₃, which have a range of crystal structures and crystallinities. The samples were characterized with X-ray diffraction (XRD), scanning electron microscopy (SEM), transmission electron microscopy (TEM), energy dispersive spectroscopy (EDS), and Brunauer-Emmett-Teller (BET) analysis. From the immersion freezing experiments, corundum [α-Al₂O₃] was shown to have the highest ice nucleation activity likely because of its high lattice match and high degree of crystallinity. Crystal structure alone did not show a strong correlation with ice nucleation activity, but a combination of a hexagonal crystal structure and a highly crystalline surface was seen to nucleate ice at warmer temperatures than the other alumina samples. This study provides experimental results in the study of ice nucleation of a range of alumina samples, which have possible implications for alumina-based mineral dust particles. Our findings suggest that crystallinity and crystal structure are important to consider when evaluating the ice nucleation efficiency of aerosol particles in laboratory and modeling studies.

Stabilization of AgI's polar surfaces by the aqueous environment, and its implications for ice formation

AgI is a potent inorganic ice nucleating particle, a feature often attributed to the lattice match between its {0001} surfaces and ice. Dissolved ions are found to be essential to the stability of these polar surfaces, and crucial to ice formation.

Does Local Structure Bias How a Crystal Nucleus Evolves?

The broad scientific and technological importance of crystallization has led to significant research probing and rationalizing crystal nucleation processes. Previous work has generally neglected the possibility of the molecular-level dynamics of individual crystal nuclei coupling to local structures. However, recent experimental work has conjectured that this can occur. Therefore, to address a deficiency in scientific understanding of crystallization, we have probed the nucleation of prototypical single and multicomponent crystals (specifically, ice and mixed gas hydrates). We establish that local structures can bias the evolution of nascent crystal phases on a nanosecond time scale by, for example, promoting the appearance or disappearance of specific crystal motifs and thus reveal a new facet of crystallization behavior. Moreover, we demonstrate structural biases are likely present during crystallization processes beyond ice and gas hydrate formation. Structurally biased dynamics are a lens for understanding existing computational and experimental results while pointing to future opportunities.

Unravelling the origins of ice nucleation on organic crystals

Organic molecules such as steroids or amino acids form crystals that can facilitate the formation of ice - arguably the most important phase transition on earth. However, the origin of the ice nucleating ability of organic crystals is still largely unknown. Here, we combine experiments and simulations to unravel the microscopic details of ice formation on cholesterol, a prototypical organic crystal widely used in cryopreservation. We find that cholesterol - which is also a substantial component of cell membranes - is an ice nucleating agent more potent than many inorganic substrates, including the mineral feldspar (one of the most active ice nucleating materials in the atmosphere). Scanning electron microscopy measurements reveal a variety of morphological features on the surfaces of cholesterol crystals: this suggests that the topography of the surface is key to the broad range of ice nucleating activity observed (from -4 to -20 °C). In addition, we show via molecular simulations that cholesterol crystals aid the formation of ice nuclei in a unconventional fashion. Rather than providing a template for a flat ice-like contact layer (as found in the case of many inorganic substrates), the flexibility of the cholesterol surface and its low density of hydrophilic functional groups leads to the formation of molecular cages involving both water molecules and terminal hydroxyl groups of the cholesterol surface. These cages are made of 6- and, surprisingly, 5-membered hydrogen bonded rings of water and hydroxyl groups that favour the nucleation of hexagonal as well as cubic ice (a rare occurrence). We argue that the phenomenal ice nucleating activity of steroids such as cholesterol (and potentially of many other organic crystals) is due to (i) the ability of flexible hydrophilic surfaces to form unconventional ice-templating structures and (ii) the different nucleation sites offered by the diverse topography of the crystalline surfaces. These findings clarify how exactly organic crystals promote the formation of ice, thus paving the way toward deeper understanding of ice formation in soft and biological matter - with obvious reverberations on atmospheric science and cryobiology.

Heterogeneous seeded molecular dynamics as a tool to probe the ice nucleating ability of crystalline surfaces

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 scales (ns) and short length scales (nm) involved. On the other hand, conventional molecular simulations struggle to cope with the relatively long time scale 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 (HSEED) approach 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.

Heterogeneous Ice Nucleation: Interplay of Surface Properties and Their Impact on Water Orientations

Ice is ubiquitous in nature, and heterogeneous ice nucleation is the most common pathway of ice formation. How surface properties affect the propensity to observe ice nucleation on that surface remains an open question. We present results of molecular dynamics studies of heterogeneous ice nucleation on model surfaces. The models surfaces considered emulate the chemistry of kaolinite, an abundant component of mineral dust. We investigate the interplay of surface lattice and hydrogen bonding properties in affecting ice nucleation. We find that lattice matching and hydrogen bonding are necessary but not sufficient conditions for observing ice nucleation at these surfaces. We correlate this behavior to the orientations sampled by the metastable supercooled water in contact with the surfaces. We find that ice is observed in cases where water molecules not only sample orientations favorable for bilayer formation but also do not sample unfavorable orientations. This distribution depends on both surface-water and water-water interactions and can change with subtle modifications to the surface properties. Our results provide insights into the diverse behavior of ice nucleation observed at different surfaces and highlight the complexity in elucidating heterogeneous ice nucleation.

Simulations of water nano-confined between corrugated planes

Water confined to nanoscale widths in two dimensions between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water to the imposition of a periodicity in the confinement by means of computer simulations, both using empirical potentials and from first-principles. For that we propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice and other values that would correspond to more realistic substrates. For the former, the phase diagram shows an overall rise of the melting temperature. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude, the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials (TIP4P/2005), both agree remarkably well for the main conclusions of the study.

The effects of ice on methane hydrate nucleation: a microcanonical molecular dynamics study

The NVE simulations realize the ice shrinking when methane hydrate nucleates both heterogeneously and homogeneously.