Hydrogen polarity of interfacial water regulates heterogeneous ice nucleation

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.

Nanoporous ices: an emerging class in the water/ice family

The history of scientific research on diverse ice structures dates back to more than a century. To date, 20 three-dimensional crystalline ice phases (ice I-ice XX) have been identified in the laboratory, among which ice XVI and ice XVII belong to a class of low-density nanoporous ices. Nanoporous ices can also be viewed as a special class of porous materials or water ice, as they possess a relatively high fraction of nano-cavities and/or nano-channels built into the hydrogen-bonded water framework. As such, like the prototypical class of porous materials (e.g., MOFs and COFs), nanoporous ices can be named as water oxygen-vertex frameworks (WOFs). Because of their large surface-to-volume ratio, WOFs may be potential media for gas storage, gas purification and separation. They may be applied to the biomedical field owing to their excellent biocompatibility. The field of porous ices is still emerging, as many porous ice structures that are predicted to be stable by computer simulations require future experimental confirmation. For future theoretical/computational studies, as the machine-learning method becomes an increasingly popular research tool in the material science and chemical science fields, more reliable porous ice structures and phase diagrams will be predicted with the development of more accurate machine-learning force fields.

Quantitative Structure-Activity Relationship Studies on Alkane Chemistry Tuning Ice Nucleation

Understanding how surface chemistry influences ice nucleation is essential for both forecasting icing phenomena and designing surfaces with desired ice-control abilities. Although alkylating is one of the most common and simplest ways for surface chemical modification, the effect of alkane chemistry on ice nucleation remains ambiguous as a result of the usually accompanying interferences of substrate morphology or heat transfer. Here, we decouple the effect of alkane chemistry on ice nucleation by investigating the ice nucleation behaviors on alkane self-assembled monolayers (SAMs) with atomic-level roughness and (sub)nanoscale thickness. Our results indicate that the introduction of alkane chemistry leads to decreased ice nucleation activities, i.e., increased anti-icing abilities, and the longer alkyl chain endows the SAM surface with the more inert ability to promote ice nucleation. The alkyl-chain-length-dependent ice nucleation activities are found to be correlated with the surface polarity. This work sheds light on a long-standing question of how alkane chemistry influences ice nucleation and offers a useful strategy for tuning ice nucleation.

Chemical Nature of Heterogeneous Electrofreezing of Supercooled Water Revealed on Polar (Pyroelectric) Surfaces

The ability to control the icing temperature of supercooled water (SCW) is of supreme importance in subfields of pure and applied sciences. The ice freezing of SCW can be influenced heterogeneously by electric effects, a process known as electrofreezing. This effect was first discovered during the 19th century; however, its mechanism is still under debate. In this Account we demonstrate, by capitalizing on the properties of polar crystals, that heterogeneous electrofreezing of SCW is a chemical process influenced by an electric field and specific ions. Polar crystals possess a net dipole moment. In addition, they are pyroelectric, displaying short-lived surface charges at their hemihedral faces at the two poles of the crystals as a result of temperature changes. Accordingly, during cooling or heating, an electric field is created, which is negated by the attraction of compensating charges from the environment. This process had an impact in the following experiments. The icing temperatures of SCW within crevices of polar crystals are higher in comparison to icing temperatures within crevices of nonpolar analogs. The role played by the electric effect was extricated from other effects by the performance of icing experiments on the surfaces of pyroelectric quasi-amorphous SrTiO₃. During those studies it was found that on positively charged surfaces the icing temperature of SCW is elevated, whereas on negatively charged surfaces it is reduced. Following investigations discovered that the icing temperature of SCW is impacted by an ionic current created within a hydrated layer on top of hydrophilic faces residing parallel to the polar axes of the crystals. In the absence of such current on analogous hydrophobic surfaces, the pyroelectric effect does not influence the icing temperature of SCW. Those results implied that electrofreezing of SCW is a process influenced by specific compensating ions attracted by the pyroelectric field from the aqueous solution. When freezing experiments are performed in an open atmosphere, bicarbonate and hydronium ions, created by the dissolution of atmospheric CO₂ in water, influence the icing temperature. The bicarbonate ions, when attracted by positively charged pyroelectric surfaces, elevate the icing temperature, whereas their counterparts, hydronium ions, when attracted by the negatively charged surfaces reduce the icing temperature. Molecular dynamic simulations suggested that bicarbonate ions, concentrated within the near positively charged interfacial layer, self-assemble with water molecules to create stabilized slightly distorted “ice-like” hexagonal assemblies which mimic the hexagons of the crystals of ice. This occurs by replacing, within those ice-like hexagons, two hydrogen bonds of water by C–O bonds of the HCO₃^– ion. On the basis of these simulations, it was predicted and experimentally confirmed that other trigonal planar ions such as NO₃^–, guanidinium⁺, and the quasi-hexagonal biguanidinium⁺ ion elevate the icing temperature. These ions were coined as “ice makers”. Other ions including hydronium, Cl^–, and SO₄^–2 interfere with the formation of ice-like assemblies and operate as “ice breakers”. The higher icing temperatures induced within the crevices of the hydrophobic polar crystals in comparison to the nonpolar analogs can be attributed to the proton ordering of the water molecules. In contrast, the icing temperatures on related hydrophilic surfaces are influenced both by compensating charges and by proton ordering.

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.

Use of Ion Exchange To Regulate the Heterogeneous Ice Nucleation Efficiency of Mica

Heterogeneous ice nucleation (HIN) triggered by mineral surfaces typically exposed to various ions can have a significant impact on the regional atmosphere and climate. However, the dependence of HIN on the nature of the mineral surface ions is still largely unexplored due to the complexity of mineral surfaces. Because K⁺ on the atomically flat (001) surface of mica can be readily replaced by different cations through ion exchange, muscovite mica was selected; its simple nature provides a very straightforward system that can serve as the model for investigating the effects of mineral surface ions on HIN. Our experiments show that the surface (001) of H⁺-exchanged mica displays markedly higher HIN efficiencies than that of Na-/K-mica. Vibrational sum-frequency generation spectroscopy reveals that H-mica induces substantially less orientation ordering than Na-/K-mica within the contact water layer at the interface. Molecular dynamics simulations suggest that the HIN efficiency of mica depends on the positional arrangement and orientation of the interfacial water. The formation of the hexagonal ice Ih basal-type structure in the first water layer atop the mica surface facilitates HIN, which is determined by the size of the protruding ions atop the mica surface and by the surface adsorption energy. The orientational distribution is optimal for HIN when 25% of the water molecules in the first water layer atop the mica surface have one OH group pointing up and 25% have one OH group pointing down, which, in turn, is determined by the surface charge distribution.