Molecular simulations of heterogeneous ice nucleation. II. Peeling back the layers

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.

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

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.

Influences of 1-Propanol and Methanol Additives on Crystallization of Thin Amorphous Solid Water Films

In general, randomly oriented ice crystallites are formed by heating amorphous solid water (ASW) films at ∼160 K via homogeneous nucleation. Here, we demonstrate that monolayers of methanol and 1-propanol additives incorporated in the multilayer ASW film lead to heterogeneous nucleation at the substrate interface of Pt(111), as evidenced by the occurrence of epitaxial ice growth. The mobility of water in direct contact with the Pt(111) substrate is decreased relative to that in the bulk, but it can be increased via interactions with hydrophobic moieties of alcohols that are segregated to the interfacial region. As a result, heterogeneous nucleation occurs at ca. 160 K along with homogeneous nucleation in the film interior. However, the template effect is quenched when the alcohols are in direct contact with the substrate. The methanol adspecies deposited onto the ASW film surface induces heterogeneous nucleation at a temperature as low as 145 K, but the 1-propanol adspecies has no such an effect. Their different ability of heterogeneous nucleation at the free ASW film surface, as well as their uptake behaviors in the near surface region, is associated with the hydrophobic hydration of the alcohols resulting from different lengths of the aliphatic moiety.

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.

Local Ice-like Structure at the Liquid Water Surface

Experiments and computer simulations have established that liquid water's surfaces can deviate in important ways from familiar bulk behavior. Even in the simplest case of an air-water interface, distinctive layering, orientational biases, and hydrogen bond arrangements have been reported, but an overarching picture of their origins and relationships has been incomplete. Here we show that a broad set of such observations can be understood through an analogy with the basal face of crystalline ice. Using simulations, we demonstrate a number of structural similarities between water and ice surfaces, suggesting the presence of domains at the air-water interface with ice-like features that persist over 2-3 molecular diameters. Most prominent is a shared characteristic layering of molecular density and orientation perpendicular to the interface. Lateral correlations of hydrogen bond network geometry point to structural similarities in the parallel direction as well. Our results bolster and significantly extend previous conceptions of ice-like structure at the liquid's boundary and suggest that the much-discussed quasi-liquid layer on ice evolves subtly above the melting point into a quasi-ice layer at the surface of liquid water.

Patterning Configuration of Surface Hydrophilicity by Graphene Nanosheet towards the Inhibition of Ice Nucleation and Growth

Freezing of liquid water occurs in many natural phenomena and affects countless human activities. The freezing process mainly involves ice nucleation and continuous growth, which are determined by the energy and structure fluctuation in supercooled water. Herein, considering the surface hydrophilicity and crystal structure differences between metal and graphene, we proposed a kind of surface configuration design, which was realized by graphene nanosheets being alternately anchored on a metal substrate. Ice nucleation and growth were investigated by molecular dynamics simulations. The surface configuration could induce ice nucleation to occur preferentially on the metal substrate where the surface hydrophilicity was higher than the lateral graphene nanosheet. However, ice nucleation could be delayed to a certain extent under the hindering effect of the interfacial water layer formed by the high surface hydrophilicity of the metal substrate. Furthermore, the graphene nanosheets restricted lateral expansion of the ice nucleus at the clearance, leading to the formation of a curved surface of the ice nucleus as it grew. As a result, ice growth was suppressed effectively due to the Gibbs–Thomson effect, and the growth rate decreased by 71.08% compared to the pure metal surface. Meanwhile, boundary misorientation between ice crystals was an important issue, which also prejudiced the growth of the ice crystal. The present results reveal the microscopic details of ice nucleation and growth inhibition of the special surface configuration and provide guidelines for the rational design of an anti-icing surface.

Impact of surface nanostructure and wettability on interfacial ice physics

Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface-ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations.

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.

Nucleation and growth of water ice on oxide surfaces: the influence of a precursor to water dissociation

We have investigated how nucleation and growth processes of ice are influenced by interfacial molecular interactions on some oxide surfaces, such as rutile TiO2(110), TiO2(100), MgO(100), and Al2O3(0001), based on the diffraction patterns of electrons transmitted through ice crystallites under the experimental configuration of reflection high energy electron diffraction (RHEED). The cubic ice Ic grows on the TiO2(110) surface with the epitaxial relationship of (110)Ic//(110)TiO2 and [001]Ic//[11[combining macron]0]TiO2. The epitaxial ice growth tends to be disturbed on the TiO2(110) surface under the presence of oxygen vacancies and adatoms. The result is not simply ascribable to small misfit values between TiO2 and ice Ic lattices (∼2%) because ice grains are formed randomly on TiO2(100). No template effects are identified during ice nucleation on the pristine MgO(100) and Al2O3(0001) surfaces either. The water molecules are chemisorbed weakly on these surfaces as a precursor to dissociation via the acid-base interaction. Such anchored water species act as an inhibitor of epitaxial ice growth because the orientation flexibility of physisorbed water during nucleation is hampered at the interface by the preferential formation of hydrogen bonds.

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.

Hydrogen polarity of interfacial water regulates heterogeneous ice nucleation

Using all-atomic molecular dynamics (MD) simulations, we show that the structure of interfacial water (IW) induced by substrates characterizes the ability of a substrate to nucleate ice. We probe the shape and structure of ice nuclei and the corresponding supercooling temperatures to measure the ability of IW with various hydrogen polarities for ice nucleation, and find that the hydrogen polarization of IW even with the ice-like oxygen lattice increases the contact angle of the ice nucleus on IW, thus lifting the free energy barrier of heterogeneous ice nucleation. The results show that not only the oxygen lattice order but the hydrogen disorder of IW on substrates are required to effectively facilitate the freezing of top water.

Anomalous Stability of Two-Dimensional Ice Confined in Hydrophobic Nanopores

The freezing of water mostly proceeds via heterogeneous ice nucleation, a process in which an effective nucleation medium not only expedites ice crystallization but also may effectively direct the polymorph selection of ice. Here, we show that water confined within a hydrophobic slit nanopore exhibits a freezing behavior strongly distinguished from its bulk counterpart. Such a difference is reflected by a strong, non-monotonic pore-size dependence of freezing temperature but, more surprisingly, by an unexpected stacking ordering of crystallized two-dimensional ice containing just a few ice layers. In particular, confined trilayer ice is found to exclusively crystallize into a well-ordered, hexagonal stacking sequence despite the fact that nanopore exerts no explicit constraint on stacking order. The absence of cubic stacking sequence is found to be originated from the intrinsically lower thermodynamic stability of cubic ice over hexagonal ice at the interface, which contrasts sharply the nearly degenerated stability of bulk hexagonal and cubic ices. Detailed examination clearly reveals that the divergence is attributed to the inherent difference between the two ice polymorphs in their surface phonon modes, which is further found to generically occur at both hydrophobic and hydrophilic surfaces.

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.

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.

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.

On the time required to freeze water

By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avrami's expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude.

Formation of Methane Hydrate in the Presence of Natural and Synthetic Nanoparticles

Natural gas hydrates occur widely on the ocean-bed and in permafrost regions, and have potential as an untapped energy resource. Their formation and growth, however, poses major problems for the energy sector due to their tendency to block oil and gas pipelines, whereas their melting is viewed as a potential contributor to climate change. Although recent advances have been made in understanding bulk methane hydrate formation, the effect of impurity particles, which are always present under conditions relevant to industry and the environment, remains an open question. Here we present results from neutron scattering experiments and molecular dynamics simulations that show that the formation of methane hydrate is insensitive to the addition of a wide range of impurity particles. Our analysis shows that this is due to the different chemical natures of methane and water, with methane generally excluded from the volume surrounding the nanoparticles. This has important consequences for our understanding of the mechanism of hydrate nucleation and the design of new inhibitor molecules.

Is Water at the Graphite Interface Vapor-like or Ice-like?

Graphitic surfaces are the main component of soot, a major constituent of atmospheric aerosols. Experiments indicate that soots of different origins display a wide range of abilities to heterogeneously nucleate ice. The ability of pure graphite to nucleate ice in experiments, however, seems to be almost negligible. Nevertheless, molecular simulations with the monatomic water model mW with water-carbon interactions parameterized to reproduce the experimental contact angle of water on graphite predict that pure graphite nucleates ice. According to classical nucleation theory, the ability of a surface to nucleate ice is controlled by the binding free energy between ice immersed in liquid water and the surface. To establish whether the discrepancy in freezing efficiencies of graphite in mW simulations and experiments arises from the coarse resolution of the model or can be fixed by reparameterization, it is important to elucidate the contributions of the water-graphite, water-ice, and ice-water interfaces to the free energy, enthalpy, and entropy of binding for both water and the model. Here we use thermodynamic analysis and free energy calculations to determine these interfacial properties. We demonstrate that liquid water at the graphite interface is not ice-like or vapor-like: it has similar free energy, entropy, and enthalpy as water in the bulk. The thermodynamics of the water-graphite interface is well reproduced by the mW model. We find that the entropy of binding between graphite and ice is positive and dominated, in both experiments and simulations, by the favorable entropy of reducing the ice-water interface. Our analysis indicates that the discrepancy in freezing efficiencies of graphite in experiments and the simulations with mW arises from the inability of the model to simultaneously reproduce the contact angle of liquid water on graphite and the free energy of the ice-graphite interface. This transferability issue is intrinsic to the resolution of the model, and arises from its lack of rotational degrees of freedom.

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.

Balance between hydration enthalpy and entropy is important for ice binding surfaces in Antifreeze Proteins

Antifreeze Proteins (AFPs) inhibit the growth of an ice crystal by binding to it. The detailed binding mechanism is, however, still not fully understood. We investigated three AFPs using Molecular Dynamics simulations in combination with Grid Inhomogeneous Solvation Theory, exploring their hydration thermodynamics. The observed enthalpic and entropic differences between the ice-binding sites and the inactive surface reveal key properties essential for proteins in order to bind ice: While entropic contributions are similar for all sites, the enthalpic gain for all ice-binding sites is lower than for the rest of the protein surface. In contrast to most of the recently published studies, our analyses show that enthalpic interactions are as important as an ice-like pre-ordering. Based on these observations, we propose a new, thermodynamically more refined mechanism of the ice recognition process showing that the appropriate balance between entropy and enthalpy facilitates ice-binding of proteins. Especially, high enthalpic interactions between the protein surface and water can hinder the ice-binding activity.

Enhanced heterogeneous ice nucleation by special surface geometry

The freezing of water typically proceeds through impurity-mediated heterogeneous nucleation. Although non-planar geometry generically exists on the surfaces of ice nucleation centres, its role in nucleation remains poorly understood. Here we show that an atomically sharp, concave wedge can further promote ice nucleation with special wedge geometries. Our molecular analysis shows that significant enhancements of ice nucleation can emerge both when the geometry of a wedge matches the ice lattice and when such lattice match does not exist. In particular, a 45° wedge is found to greatly enhance ice nucleation by facilitating the formation of special topological defects that consequently catalyse the growth of regular ice. Our study not only highlights the active role of defects in nucleation but also suggests that the traditional concept of lattice match between a nucleation centre and crystalline lattice should be extended to include a broader match with metastable, non-crystalline structural motifs.

Understanding ice nucleation is important for the development of accurate cloud models. Here Bi et al. show that sharp wedges can enhance ice nucleation both when the wedge geometry matches the ice lattice and when such matching is absent, in which case nucleation is promoted by topological defects.

Ice Nucleation Efficiency of Hydroxylated Organic Surfaces Is Controlled by Their Structural Fluctuations and Mismatch to Ice

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.

Boreal pollen contain ice-nucleating as well as ice-binding 'antifreeze' polysaccharides

Ice nucleation and growth is an important and widespread environmental process. Accordingly, nature has developed means to either promote or inhibit ice crystal formation, for example ice-nucleating proteins in bacteria or ice-binding antifreeze proteins in polar fish. Recently, it was found that birch pollen release ice-nucleating macromolecules when suspended in water. Here we show that birch pollen washing water exhibits also ice-binding properties such as ice shaping and ice recrystallization inhibition, similar to antifreeze proteins. We present spectroscopic evidence that both the ice-nucleating as well as the ice-binding molecules are polysaccharides bearing carboxylate groups. The spectra suggest that both polysaccharides consist of very similar chemical moieties, but centrifugal filtration indicates differences in molecular size: ice nucleation occurs only in the supernatant of a 100 kDa filter, while ice shaping is strongly enhanced in the filtrate. This finding may suggest that the larger ice-nucleating polysaccharides consist of clusters of the smaller ice-binding polysaccharides, or that the latter are fragments of the ice-nucleating polysaccharides. Finally, similar polysaccharides released from pine and alder pollen also display both ice-nucleating as well as ice-binding ability, suggesting a common mechanism of interaction with ice among several boreal pollen with implications for atmospheric processes and antifreeze protection.

Heterogeneous Nucleation of an n-Alkane on Tetrahedrally Coordinated Crystals

Heterogeneous nucleation refers to the propensity for phase transformations to initiate preferentially on foreign surfaces, such as vessel walls, dust particles, or formulation additives. In crystallization, the form of the initial nucleus has ramifications for the crystallographic form, morphology, and properties of the resulting solid. Nevertheless, the discovery and design of nucleating agents remains a matter of trial and error because of the very small spatiotemporal scales over which the critical nucleus is formed and the extreme difficulty of examining such events empirically. Using molecular dynamics simulations, we demonstrate a method for the rapid screening of entire families of materials for activity as nucleating agents and for characterizing their mechanism of action. The method is applied to the crystallization of n-pentacontane, a model surrogate for polyethylene, on the family of tetrahedrally coordinated crystals, including diamond and silicon. A systematic variation of parameters in the interaction potential permits a comprehensive, physically based screening of nucleating agents in this class of materials, including both real and hypothetical candidates. The induction time for heterogeneous nucleation is shown to depend strongly on crystallographic registry between the nucleating agent and the critical nucleus, indicative of an epitaxial mechanism in this class of materials. Importantly, the severity of this registry requirement weakens with decreasing rigidity of the substrate and increasing strength of attraction to the surface of the nucleating agent. Employing this method, a high-throughput computational screening of nucleating agents becomes possible, facilitating the discovery of novel nucleating agents within a broad "materials genome" of possible additives.

Tuning Ice Nucleation with Supercharged Polypeptides

Supercharged unfolded polypeptides (SUPs) are exploited for controlling ice nucleation via tuning the nature of charge and charge density of SUPs. The results show that positively charged SUPs facilitate ice nucleation, while negatively charged ones suppress it. Moreover, the charge density of the SUP backbone is another parameter to control it.

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

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 insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past 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 such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, 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 toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Can Ice-Like Structures Form on Non-Ice-Like Substrates? The Example of the K-feldspar Microcline

Feldspar minerals are the most common rock formers in Earth's crust. As such they play an important role in subjects ranging from geology to climate science. An atomistic understanding of the feldspar structure and its interaction with water is therefore desirable, not least because feldspar has been shown to dominate ice nucleation by mineral dusts in Earth's atmosphere. The complexity of the ice/feldspar interface arising from the numerous chemical motifs expressed on the surface makes it a challenging system. Here we report a comprehensive study of this challenging system with ab initio density functional theory calculations. We show that the distribution of Al atoms, which is crucial for the dissolution kinetics of tectosilicate minerals, differs significantly between the bulk environment and on the surface. Furthermore, we demonstrate that water does not form ice-like overlayers in the contact layer on the most easily cleaved (001) surface of K-feldspar. We do, however, identify contact layer structures of water that induce ice-like ordering in the second overlayer. This suggests that even substrates without an apparent match with the ice structure may still act as excellent ice nucleating agents.