Heterogeneous nucleation of ice on carbon surfaces

Ice-Nucleating and Antifreeze Proteins Recognize Ice through a Diversity of Anchored Clathrate and Ice-like Motifs

Cold-adapted organisms produce antifreeze and ice-nucleating proteins to prevent and promote ice formation. The crystal structure of hyperactive bacterial antifreeze protein (AFP) MpAFP suggests that this protein binds ice through an anchored clathrate motif. It is not known whether other hyperactive AFPs and ice-nucleating proteins (INPs) use the same motif to recognize or nucleate ice. Here we use molecular simulations to elucidate the ice-binding motifs of hyperactive insect AFPs and a model INP of Pseudomonas syringae. We find that insect AFPs recognize ice through anchored clathrate motifs distinct from that of MpAFP. By performing simulations of ice nucleation by PsINP, we identify two distinct ice-binding sites on opposite sides of the β-helix. The ice-nucleating sequences identified in the simulations agree with those previously proposed for the closely related INP of Pseudomonas borealis based on the structure of the protein. The simulations indicate that these sites have comparable ice-nucleating efficiency, but distinct binding motifs, controlled by the amino acid sequence: one is an anchored clathrate and the other ice-like. We conclude that anchored clathrate and ice-like motifs can be equally effective for binding proteins to ice and promoting ice nucleation.

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.

Impact of sequential surface-modification of graphene oxide on ice nucleation

Base-washed graphene-oxide which has been sequentially-modified by thiol-epoxy chemistry, results in materials with ice-nucleation activity. The role of hydro-philic/phobic grafts and polymers was evaluated with the most potent functioning at just 0.25 wt%. These 2-D hybrid materials may find use in cryopreservation and fundamental studies on ice formation.

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.

Control of ice nucleation: freezing and antifreeze strategies

Water freezing remains a perennial topic of great relevance to many important aspects of our lives; from the climate to human society and from economics to medicine, frozen water profoundly influences our living environment and life activities.

Pre-critical fluctuations and what they disclose about heterogeneous crystal nucleation

Heterogeneous crystal nucleation is ubiquitous in nature and at the heart of many industrial applications. At the molecular scale, however, major gaps in understanding this phenomenon persist. Here we investigate through molecular dynamics simulations how the formation of precritical crystalline clusters is connected to the kinetics of nucleation. Considering heterogeneous water freezing as a prototypical scenario of practical relevance, we find that precritical fluctuations connote which crystalline polymorph will form. The emergence of metastable phases can thus be promoted by templating crystal faces characteristic of specific polymorphs. As a consequence, heterogeneous classical nucleation theory cannot describe our simulation results, because the different substrates lead to the formation of different ice polytypes. We discuss how the issue of polymorphism needs to be incorporated into analysis and comparison of heterogeneous and homogeneous nucleation. Our results will help to interpret and analyze the growing number of experiments and simulations dealing with crystal polymorph selection.

Promotion of Homogeneous Ice Nucleation by Soluble Molecules

Atmospheric aerosols nucleate ice in clouds, strongly impacting precipitation and climate. The prevailing consensus is that ice nucleation is promoted heterogeneously by the surface of ice nucleating particles in the aerosols. However, recent experiments indicate that water-soluble molecules, such as polysaccharides of pollen and poly(vinyl alcohol) (PVA), increase the ice freezing temperature. This poses the question of how do flexible soluble molecules promote the formation of water crystals, as they do not expose a well-defined surface to ice. Here we use molecular simulations to demonstrate that PVA promotes ice nucleation through a homogeneous mechanism: PVA increases the nucleation rate by destabilizing water in the solution. This work demonstrates a novel paradigm for understanding ice nucleation by soluble molecules and provides a new handle to design additives that promote crystallization.

CH4 Hydrate Formation between Silica and Graphite Surfaces: Insights from Microsecond Molecular Dynamics Simulations

Microsecond simulations have been performed to investigate CH₄ hydrate formation from gas/water two-phase systems between silica and graphite surfaces, respectively. The hydrophilic silica and hydrophobic graphite surfaces exhibit substantially different effects on CH₄ hydrate formation. The graphite surface adsorbs CH₄ molecules to form a nanobubble with a flat or negative curvature, resulting in a low aqueous CH₄ concentration, and hydrate nucleation does not occur during 2.5 μs simulation. Moreover, an ordered interfacial water bilayer forms between the nanobubble and graphite surface thus preventing their direct contact. In contrast, the hydroxylated-silica surface prefers to be hydrated by water, with a cylindrical nanobubble formed in the solution, leading to a high aqueous CH₄ concentration and hydrate nucleation in the bulk region; during hydrate growth, the nanobubble is gradually covered by hydrate solid and separated from the water phase, hence slowing growth. The silanol groups on the silica surface can form strong hydrogen bonds with water, and hydrate cages need to match the arrangements of silanols to form more hydrogen bonds. At the end of the simulation, the hydrate solid is separated from the silica surface by liquid water, with only several cages forming hydrogen bonds with the silica surface, mainly due to the low CH₄ aqueous concentrations near the surface. To further explore hydrate formation between graphite surfaces, CH₄/water homogeneous solution systems are also simulated. CH₄ molecules in the solution are adsorbed onto graphite and hydrate nucleation occurs in the bulk region. During hydrate growth, the adsorbed CH₄ molecules are gradually converted into hydrate solid. It is found that the hydrate-like ordering of interfacial water induced by graphite promotes the contact between hydrate solid and graphite. We reveal that the ability of silanol groups on silica to form strong hydrogen bonds to stabilize incipient hydrate solid, as well as the ability of graphite to adsorb CH₄ molecules and induce hydrate-like ordering of the interfacial water, are the key factors to affect CH₄ hydrate formation between silica and graphite surfaces.

Heterogeneous Freezing of Carbon Nanotubes: A Model System for Pore Condensation and Freezing in the Atmosphere

Heterogeneous ice nucleation is an important mechanism for cloud formation in the upper troposphere. Recently, pores on atmospheric particles have been proposed to play a significant role in ice nucleation. To understand how ice nucleation occurs in idealized pores, we characterized the immersion freezing activity of various sizes of carbon nanotubes. Carbon nanotubes are used both as a model for pores and proxy for soot particles. We determined that carbon nanotubes with inner diameters between 2 and 3 nm exhibit the highest ice nucleation activity. Implications for the freezing behavior of porous materials and nucleation on soot particles will be discussed.

Reaction Coordinate for Ice Crystallization on a Soft Surface

The control of assembly and crystallization of molecules is becoming increasingly important in chemistry, engineering, and materials sciences. Crystallization is also central to understand natural processes that include the formation of atmospheric ice and biomineralization. Organic surfaces, biomolecules, and even liquid/vapor interfaces can promote the nucleation of crystals. These soft surfaces present significant structural fluctuations, which have been shown to strongly impact the rate of crystallization. This raises the question of whether degrees of freedom of soft surfaces play a role in the reaction coordinate for crystal nucleation. Here we use molecular simulations to investigate the mechanism of ice nucleation promoted by an alcohol monolayer. Our analysis indicates that while the flexibility of the surface strongly depresses its ice nucleation ability, it does not play a role in the coordinate that controls the transformation from liquid to ice. We find that the variable that drives the transformation is the size of the crystalline cluster, the same as that for the homogeneous crystallization. We argue that this is a general result that arises from the separation of time scales between surface fluctuations and the crossing of the transition state barrier for crystallization.

Inhibition of Heterogeneous Ice Nucleation by Bioinspired Coatings of Polyampholytes

Control of heterogeneous ice nucleation (HIN) on foreign surfaces is of great importance for anti-ice-nucleation material design. In this work, we studied the HIN behaviors on various ion-modified poly(butylene succinate) (PBS) surfaces via chain-extension reaction. Inspired by antifreeze proteins (AFPs), the PBS-based polyampholytes, containing both negative and positive charge groups on a single chain, show excellent performance of ice nucleation inhibition and freezing delay. Unlike the extremely high price and low availability of AFPs, these PBS-based polyampholytes can be commercially synthesized under mild reaction conditions. Through water freezing tests on a wide range of substrates at different temperatures, these PBS-based polyampholytes have shown application value of tuning ice nucleation via a simple spin-coating method.

Interfacial Water Arrangement in the Ice-Bound State of an Antifreeze Protein: A Molecular Dynamics Simulation Study

Molecular dynamics (MD) simulations have been carried out to study the heterogeneous ice nucleation on modeled peptide surfaces. Simulations show that large peptide surfaces made by TxT (threonine-x-threonine) motifs with the arrangements of threonine (Thr) residues identical to the periodic arrangements of waters on either the basal or prism plane of ice are capable of ice nucleation. Nucleated ice plane is the (0001) basal plane of hexagonal ice (Ih) or (111) plane of cubic ice (Ic). However, due to predefined simulation cell dimensions, the ice growth is only observed on the surface where the Thr residues are arranged like the water arrangement on the basal plane of ice Ih. The γ-methyl and γ-hydroxyl groups of Thr residue are necessary for such ice formation. From this ice nucleation and growth simulation, the interfacial water arrangement in the ice-bound state of Tenebrio molitor antifreeze protein (TmAFP) has been determined. The interfacial water arrangement in the ice-bound state of TmAFP is characterized by five-membered hydrogen bonded rings, where each of the hydroxyl groups of the Thr residues on the ice-binding surface (IBS) of the protein is a ring member. It is found that the water arrangement at the protein-ice interface is distorted from that in bulk ice. Our analysis further reveals that the hydroxyl groups of Thr residues on the IBS of TmAFP form maximum three hydrogen bonds each with the waters in the bound state and methyl groups of Thr residues occupy wider spaces than the normal grooves on the (111) plane of ice Ic. Methyl groups are also located above and along the 3-fold rotational axes of the chair-formed hexagonal hydrogen bonded water rings on the (111) plane.

Overview of Ice Nucleating Particles

Ice particle formation in tropospheric clouds significantly changes cloud radiative and microphysical properties. Ice nucleation in the troposphere via homogeneous freezing occurs at temperatures lower than −38°C and relative humidity with respect to ice above 140%. In the absence of these conditions, ice formation can proceed via heterogeneous nucleation aided by aerosol particles known as ice nucleating particles (INPs). In this chapter, new developments in identifying the heterogeneous freezing mechanisms, atmospheric relevance, uncertainties, and unknowns about INPs are described. The change in conventional wisdom regarding the requirements of INPs as new studies discover physical and chemical properties of these particles is explained. INP sources and known reasons for their ice nucleating properties are presented. The need for more studies to systematically identify particle properties that facilitate ice nucleation is highlighted. The atmospheric relevance of long-range transport, aerosol aging, and coating studies (in the laboratory) of INPs are also presented. Possible mechanisms for processes that change the ice nucleating potential of INPs and the corresponding challenges in understanding and applying these in models are discussed. How primary ice nucleation affects total ice crystal number concentrations in clouds and the discrepancy between INP concentrations and ice crystal number concentrations are presented. Finally, limitations of parameterizing INPs and of models in representing known and unknown processes related to heterogeneous ice nucleation processes are discussed.

Insight on Structure of Water and Ice Next to Graphene Using Surface-Sensitive Spectroscopy

The water/graphene interface has received considerable attention in the past decade due to its relevance in various potential applications including energy storage, sensing, desalination, and catalysis. Most of our knowledge about the interfacial water structure next to graphene stems from simulations, which use experimentally measured water contact angles (WCAs) on graphene (or graphite) to estimate the water-graphene interaction strength. However, the existence of a wide spectrum of reported WCAs on supported graphene and graphitic surfaces makes it difficult to interpret the water-graphene interactions. Here, we have used surface-sensitive infrared-visible sum frequency generation (SFG) spectroscopy to probe the interfacial water structure next to graphene supported on a sapphire substrate. In addition, the ice nucleation properties of graphene have been explored by performing in situ freezing experiments as graphitic surfaces are considered good ice nucleators. For graphene supported on sapphire, we observed a strong SFG peak associated with highly coordinated, ordered water next to graphene. Similar ordering was not detected next to bare sapphire, implying that the observed ordering of water molecules in the former case is a consequence of the presence of graphene. Our analysis indicates that graphene behaves like a hydrophobic (or negatively charged) surface, leading to enhanced ordering of water molecules. Although liquid water orders next to graphene, the ice formed is proton disordered. This research sheds light on water-graphene interactions relevant in optimizing the performance of graphene in various applications.

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.

Parameterization of a coarse-grained model with short-ranged interactions for modeling fuel cell membranes with controlled water uptake

The coarse-grained model FF_pvap reproduces the experimental activity coefficient of water in tetramethylammonium chloride solutions over a wide range of concentrations, with a hundred-fold gain in computing efficiency with respect to atomistic models.

Active sites in heterogeneous ice nucleation—the example of K-rich feldspars

Ice formation on aerosol particles is a process of crucial importance to Earth’s climate and the environmental sciences, but it is not understood at the molecular level. This is partly because the nature of active sites, local surface features where ice growth commences, is still unclear. Here we report direct electron-microscopic observations of deposition growth of aligned ice crystals on feldspar, an atmospherically important component of mineral dust. Our molecular-scale computer simulations indicate that this alignment arises from the preferential nucleation of prismatic crystal planes of ice on high-energy (100) surface planes of feldspar. The microscopic patches of (100) surface, exposed at surface defects such as steps, cracks, and cavities, are thought to be responsible for the high ice nucleation efficacy of potassium (K)–feldspar particles.

Foreign Particle Promoted Crystalline Nucleation for Growing High-Quality Ultrathin Rubrene Films

Introducing foreign particles or agents as nucleator is an efficient way to promote crystallization in the crystal growth field, with the advantage to speed up the crystallizing rate and control the growth process. However, in the field of organic crystalline film growth, where the crystallization and morphology modulation are of significant importance in optoelectronics, this method has rarely been utilized. Particularly, some potential high-performance materials such as rubrene face the problem of crystallization during film formation. Here a strategy is reported to promote the crystallization of rubrene films in the initial stage assisted by foreign particles. Highly ordered thin film from the sub-monolayer stage can be achieved. Efficient charge transport and high mobility up to 2.95 cm(2) V(-1) s(-1) are achieved on thus ultrathin crystalline films. Such a method enables the well controlling of the film growth from the very early stage and produces uniform crystalline films with good reproducibility, thus highly promising to yield desired optoelectrical properties and applications.

Effect of Surface Energy on Freezing Temperature of Water

Previous studies have found that superhydrophobic surfaces are effective in delaying freezing of water droplets. However, the freezing process of water droplets on superhydrophobic surfaces depends on factors such as droplet size, surface area, roughness, and cooling rate. The role of surface energy, independent of any other parameters, in delaying freezing of water is not understood. Here, we have used infrared-visible sum frequency generation spectroscopy (SFG) to study the freezing of water next to solid substrates with water contact angles varying from 5° to 110°. We find that the freezing temperature of water decreases with increasing surface hydrophobicity only when the sample volume is small (∼10 μL). For a larger volume of water (∼300 μL), the freezing temperature is independent of surface energy. For water next to the surfaces with contact angle ≥54°, we observe a strong SFG peak associated with highly coordinated water. This research sheds new light on understanding the key factors in designing new anti-icing coatings.

What Determines the Ice Polymorph in Clouds?

Ice crystals in the atmosphere nucleate from supercooled liquid water and grow by vapor uptake. The structure of the ice polymorph grown has strong impact on the morphology and light scattering of the ice crystals, modulates the amount of water vapor in ice clouds, and can impact the molecular uptake and reactivity of atmospheric aerosols. Experiments and molecular simulations indicate that ice nucleated and grown from deeply supercooled liquid water is metastable stacking disordered ice. The ice polymorph grown from vapor has not yet been determined. Here we use large-scale molecular simulations to determine the structure of ice that grows as a result of uptake of water vapor in the temperature range relevant to cirrus and mixed-phase clouds, elucidate the molecular mechanism of the formation of ice at the vapor interface, and compute the free energy difference between cubic and hexagonal ice interfaces with vapor. We find that vapor deposition results in growth of stacking disordered ice only under conditions of extreme supersaturation, for which a nonequilibrium liquid layer completely wets the surface of ice. Such extreme conditions have been used to produce stacking disordered frost ice in experiments and may be plausible in the summer polar mesosphere. Growth of ice from vapor at moderate supersaturations in the temperature range relevant to cirrus and mixed-phase clouds, from 200 to 260 K, produces exclusively the stable hexagonal ice polymorph. Cubic ice is disfavored with respect to hexagonal ice not only by a small penalty in the bulk free energy (3.6 ± 1.5 J mol(-1) at 260 K) but also by a large free energy penalty at the ice-vapor interface (89.7 ± 12.8 J mol(-1) at 260 K). The latter originates in higher vibrational entropy of the hexagonal-terminated ice-vapor interface. We predict that the free energy penalty against the cubic ice interface should decrease strongly with temperature, resulting in some degree of stacking disorder in ice grown from vapor in the tropical tropopause layer, and in polar stratospheric and noctilucent clouds. Our findings support and explain the evolution of the morphology of ice crystals from hexagonal to trigonal symmetry with decreasing temperature, as reported by experiments and in situ measurements in clouds. We conclude that selective growth of the elusive cubic ice polymorph by manipulation of the interfacial properties can likely be achieved at the ice-liquid interface but not at the ice-vapor interface.

Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite

Most ice in nature forms because of impurities which boost the exceedingly low nucleation rate of pure supercooled water. However, the microscopic details of ice nucleation on these substances remain largely unknown. Here, we have unraveled the molecular mechanism and the kinetics of ice formation on kaolinite, a clay mineral playing a key role in climate science. We find that the formation of ice at strong supercooling in the presence of this clay is about 20 orders of magnitude faster than homogeneous freezing. The critical nucleus is substantially smaller than that found for homogeneous nucleation and, in contrast to the predictions of classical nucleation theory (CNT), it has a strong two-dimensional character. Nonetheless, we show that CNT describes correctly the formation of ice at this complex interface. Kaolinite also promotes the exclusive nucleation of hexagonal ice, as opposed to homogeneous freezing where a mixture of cubic and hexagonal polytypes is observed.

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.

Tuning ice nucleation with counterions on polyelectrolyte brush surfaces

Heterogeneous ice nucleation (HIN) on ionic surfaces is ubiquitous in a wide range of atmospheric aerosols and at biological interfaces. Despite its great importance in cirrus cloud formation and cryopreservation of cells, organs, and tissues, it remains unclear whether the ion-specific effect on ice nucleation exists. Benefiting from the fact that ions at the polyelectrolyte brush (PB)/water interface can be reversibly exchanged, we report the effect of ions on HIN on the PB surface, and we discover that the distinct efficiency of ions in tuning HIN follows the Hofmeister series. Moreover, a large HIN temperature window of up to 7.8°C is demonstrated. By establishing a correlation between the fraction of ice-like water molecules and the kinetics of structural transformation from liquid- to ice-like water molecules at the PB/water interface with different counterions, we show that our molecular dynamics simulation analysis is consistent with the experimental observation of the ion-specific effect on HIN.

Ice-nucleating bacteria control the order and dynamics of interfacial water

Ice-nucleating organisms play important roles in the environment. With their ability to induce ice formation at temperatures just below the ice melting point, bacteria such as Pseudomonas syringae attack plants through frost damage using specialized ice-nucleating proteins. Besides the impact on agriculture and microbial ecology, airborne P. syringae can affect atmospheric glaciation processes, with consequences for cloud evolution, precipitation, and climate. Biogenic ice nucleation is also relevant for artificial snow production and for biomimetic materials for controlled interfacial freezing. We use interface-specific sum frequency generation (SFG) spectroscopy to show that hydrogen bonding at the water-bacteria contact imposes structural ordering on the adjacent water network. Experimental SFG data and molecular dynamics simulations demonstrate that ice-active sites within P. syringae feature unique hydrophilic-hydrophobic patterns to enhance ice nucleation. The freezing transition is further facilitated by the highly effective removal of latent heat from the nucleation site, as apparent from time-resolved SFG spectroscopy.

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.

Free energy contributions and structural characterization of stacking disordered ices

Crystallization of ice from deeply supercooled water and amorphous ices - a process of fundamental importance in the atmosphere, interstellar space, and cryobiology - results in stacking disordered ices with a wide range of metastabilities with respect to hexagonal ice. The structural origin of this high variability, however, has not yet been elucidated. Here we use molecular dynamics simulations with the mW water model to characterize the structure of ice freshly grown from supercooled water at temperatures from 210 to 270 K, the thermodynamics of stacking faults, line defects, and interfaces, and to elucidate the interplay between kinetics and thermodynamics in determining the structure of ice. In agreement with experiments, the ice grown in the simulations is stacking disordered with a random distribution of cubic and hexagonal layers, and a cubicity that decreases with growth temperature. The former implies that the cubicity of ice is determined by processes at the ice/liquid interface, without memory of the structure of buried ice layers. The latter indicates that the probability of building a cubic layer at the interface decreases upon approaching the melting point of ice, which we attribute to a more efficient structural equilibration of ice at the liquid interface as the driving force for growth wanes. The free energy cost for creating a pair of cubic layers in ice is 8.0 J mol(-1) in experiments, and 9.7 ± 1.9 J mol(-1) for the mW water model. This not only validates the simulations, but also indicates that dispersion in cubicity is not sufficient to explain the large energetic variability of stacking disordered ices. We compute the free energy cost of stacking disorder, line defects, and interfaces in ice and conclude that a characterization of the density of these defects is required to predict the degree of metastability and vapor pressure of atmospheric ices.

Ice nucleation on nanotextured surfaces: the influence of surface fraction, pillar height and wetting states

Ice nucleation and growth on nanostructured surfaces.