Growth and Structure of the First Layers of Ice on Ru(0001) and Pt(111)

Phase Transition of Cubic Ice to Hexagonal Ice during Growth and Decomposition

Ice, one of the most enigmatic materials on Earth, exhibits diverse polymorphism, with research mainly focusing on the most commonly observed phases: hexagonal ice (Ih), cubic ice (Ic), and stacking-disordered ice (Isd). While their formation or structural changes are crucial for advancements in cloud science, climate modeling, and cryogenic technology, the molecular mechanisms driving these phenomena remain unexplored. Herein, utilizing cryogenic transmission electron microscopy, we investigate the formation of ice at two different temperatures, demonstrating a size-dependent phase shift from Ic to Isd. Furthermore, a relatively metastable cubic phase in Isd transitions to a hexagonal phase under electron beam radiation. This transition, facilitated by crystal defects, contrasts with perfect crystalline Ic, which maintains its original phase, emphasizing the importance of defects in polymorphic phase transitions. Our findings provide novel insights on phase control during the ice growth processes and polymorphic phase transitions from the cubic-to-hexagonal phases.

Icing and Adhesion Behaviors on Surfaces with Varied Lattice Constants

Investigating droplet wetting and icing behavior is crucial for comprehending the principles of surface icing and the design of anti-icing surfaces. In this study, we present the evidence from molecular dynamics (MD) simulations that reveal a hitherto unreported behavior of droplet wetting and icing adhesion on surfaces with lattice constants from 2.7 to 4.5 Å. Here, we observe that the contact angles (CA) of droplets on a face-centered cubic (FCC) lattice surface consistently correlate positively with the lattice constant. Further examination of droplet behavior on an idealized crystal surface reveals that hydrophilic surfaces (e.g., CA = 85°) inhibit freezing more effectively than hydrophobic surfaces (e.g., CA = 97°). This finding contradicts the conventional explanation that hydrophobic surfaces reduce heterogeneous nucleation, thereby delaying icing. This study introduces a mechanistic explanation for the promotion of water icing by hydrophobic surfaces and offers a novel design concept for the development of anti-ice surfaces in future applications.

Two-dimensional bilayer ice in coexistence with three-dimensional ice without confinement

Icing plays an important role in various physical-chemical process. Although the formation of two-dimensional ice requires nanoscale confinement, two-dimensional bilayer ice in coexistence with three-dimensional ice without confinement remains poorly understood. Here, a critical value of a surface energy parameter is identified to characterize the liquid-solid interface interaction, above which two-dimensional and three-dimensional coexisting ice can surprisingly form on the surface. The two-dimensional ice growth mechanisms could be revealed by capturing the growth and merged of the metastable edge structures. The phase diagram about temperature and pressure vs energy parameters is predicted to distinguish liquid water, two-dimensional ice and three-dimensional ice. Furthermore, the deicing characteristics of coexisting ice demonstrate that the ice adhesion strength is linearly related to the ratio of ice-surface interaction energy to ice temperature. In addition, for gas-solid phase transition, the phase diagram about temperature and energy parameters is predicted to distinguish gas, liquid water, two-dimensional ice and three-dimensional ice. This work gives a perspective for studying the singular structure and dynamics of ice in nanoscale and provides a guide for future experimental realization of the coexisting ice.

Two-dimensional ice-like water adlayers on a mica surface with and without a graphene coating under ambient conditions

Water tends to wet all hydrophilic surfaces under ambient conditions, and the first water adlayers on solids are important for a broad range of physicochemical phenomena and technological processes, including corrosion, wetting, lubrication, anti-icing, catalysis, and electrochemistry. Unfortunately, challenges in characterizing the first water adlayer in the laboratory have hampered molecular-level understanding of the contact water structure. Herein, we present the first ab initio molecular dynamics simulation evidence of a previously unreported ice-like adlayer structure (named as Ice-AL-II) on a prototype mica surface under ambient conditions. Calculation showed that the newly identified Ice-AL-II structure is more stable than the widely recognized ice-adlayer structure on mica surfaces (named as Ice-AL-I). Ice-AL-II exhibited a face-centered corner-cut tetragon (or a face-centered irregular pentagon) pattern of a hydrogen-bonded network. The center of the corner-cut tetragon was occupied by either a K+ cation or a water molecule with two H atoms pinned by the mica (100) via double hydrogen bonds. Our simulation also suggested that bilayer Ice-AL-II favors AA stacking rather than AB stacking. Interestingly, when a graphene sheet was coated on top of the ice-like adlayer, the stability of Ice-AL-II was further enhanced. In contrast, due to its strongly puckered structure, the Ice-AL-I structure could be crushed into a near-Ice-AL-II structure by the graphene coating. Ice-AL-II is thus proposed as a promising candidate for the ice-like structure on a mica surface detected by scanning polarization force microscopy and by atomic force microscopy between a graphene coating and a mica surface.

Two rhombic ice phases from aqueous salt solutions under graphene confinement

Water exhibits rich ice phases depending upon its respective formation conditions, and in particular, the two-dimensional ice with nonhexagonal symmetry adsorbed on solids relates to the exceptional arrangement of water molecules. Despite extensive reporting of two-dimensional ice on various solid surfaces, the geometry and thermodynamics of ice formation from an aqueous salt solution are still unknown. In this Letter, we show the formation of single- and two-phase mixed two-dimensional rhombic ice from aqueous salt solutions with different concentrations under strong compressed confinement of graphene at ambient temperature by using classical molecular dynamics simulations and first-principles calculations. The two rhombic ice phases exhibit identical geometry and thermodynamic properties, but different projections of the oxygen atoms against solid surface symmetry, where they relate to the stable and metastable arrangements of water molecules confined between two graphene layers. A single-phase rhombic ice would grow from the confined saturated aqueous solutions since the previously stable rhombic molecular arrangement becomes an unstable high-energy state by introducing salt ions nearby. Our result reveals different rhombic ice phases growing from pure water and aqueous solutions, highlighting the deciding role of salt ions in the ice formation process due to their common presence in liquids.

Imaging surface structure and premelting of ice Ih with atomic resolution

Ice surfaces are closely relevant to many physical and chemical properties, such as melting, freezing, friction, gas uptake and atmospheric reaction1-8. Despite extensive experimental and theoretical investigations9-17, the exact atomic structures of ice interfaces remain elusive owing to the vulnerable hydrogen-bonding network and the complicated premelting process. Here we realize atomic-resolution imaging of the basal (0001) surface structure of hexagonal water ice (ice Ih) by using qPlus-based cryogenic atomic force microscopy with a carbon monoxide-functionalized tip. We find that the crystalline ice-Ih surface consists of mixed Ih- and cubic (Ic)-stacking nanodomains, forming 19 × 19 periodic superstructures. Density functional theory reveals that this reconstructed surface is stabilized over the ideal ice surface mainly by minimizing the electrostatic repulsion between dangling OH bonds. Moreover, we observe that the ice surface gradually becomes disordered with increasing temperature (above 120 Kelvin), indicating the onset of the premelting process. The surface premelting occurs from the defective boundaries between the Ih and Ic domains and can be promoted by the formation of a planar local structure. These results put an end to the longstanding debate on ice surface structures and shed light on the molecular origin of ice premelting, which may lead to a paradigm shift in the understanding of ice physics and chemistry.

Atomic force microscopy observation of surface morphologies and measurements of local contact potential differences of amorphous solid water samples deposited at 15 and 100 K

We use an ultra-high vacuum cryogenic atomic force microscope to investigate the surface morphology of amorphous solid water (ASW) prepared by oblique deposition of water vapor onto Si(111)7 × 7 substrates at temperatures of 15 and 100 K. Height-height correlation function analysis of topographic images suggests that ASW at 15 K has a columnar structure and that the typical diameter of the column is 5-10 nm. At 100 K, the typical diameter is 10-30 nm, although columnar features are less prominent. The surface roughness (i.e., deviation of the height) is greater at 15 K than at 100 K, indicating that the surface at 100 K exhibits a relatively flat morphology. This result implies that transient diffusion of deposited water molecules affects the surface morphology at 100 K. In addition, measurements of the local contact potential difference between the tip and the ASW surface suggest that the magnitude of the negative surface potential at the microscopic scale, which is attributed to spontaneous polarisation, cannot simply be scaled by the thickness of ASW as predicted in previous experiments with Kelvin probes.

The structure-giving role of Rb+ ions for water-ice nanoislands supported on Cu(111)

We characterize the effect of rubidium ions on water-ice nanoislands in terms of area, fractal dimension, and apparent height by low-temperature scanning tunneling microscopy. Water nanoislands on the pristine Cu(111) surface are compared to those at similar coverage on a Rb+ pre-covered Cu(111) surface to reveal the structure-giving effect of Rb+. The presence of Rb+ induces changes in the island shape, and hence, the water network, without affecting the nanoisland volume. The broad area distribution shifts to larger values while the height decreases from three bilayers to one or two bilayers. The nanoislands on the Rb+ pre-covered surface are also more compact, reflected in a shift in the fractal dimension distribution. We relate the changes to a weakening of the hydrogen-bond network by Rb+.

The effect of surface hydrophobicity and hydrophilicity on ion-ion interactions at water-solid interfaces

Condensation and arrangement of ions at water-solid interfaces are of great importance in the formation of electrical double layers (EDL) and the transport of ions under a confined geometry. So far, the microscopic understanding of interfacial ion configurations is still far from complete, especially when the local ion concentration is high and ion-ion interactions become prominent. In this study, we directly visualized alkali metal cations within the hydrogen-bonding network of water on graphite and Cu(111)-supported graphene surfaces, using qPlus-based noncontact atomic force microscopy (NC-AFM). We found that the codeposition of the alkali cations and water molecules on the hydrophobic graphite surface leads to the formation of an ion-doped bilayer hexagonal ice (BHI) structure, where the ions are repelled from each other and scattered in a disordered distribution. In contrast, the hydrated alkali cations aggregate in one dimension on the more hydrophilic graphene/Cu(111) surface, forming a nematic state with a long-range order. Such a nematic state arises from the delicate interplay between water-ion and water-water interactions under surface confinement. These results reveal the high sensitivity of ion-ion interactions and ionic ordering to the surface hydrophobicity and hydrophilicity.

Structure Discovery in Atomic Force Microscopy Imaging of Ice

The interaction of water with surfaces is crucially important in a wide range of natural and technological settings. In particular, at low temperatures, unveiling the atomistic structure of adsorbed water clusters would provide valuable data for understanding the ice nucleation process. Using high-resolution atomic force microscopy (AFM) and scanning tunneling microscopy, several studies have demonstrated the presence of water pentamers, hexamers, and heptamers (and of their combinations) on a variety of metallic surfaces, as well as the initial stages of 2D ice growth on an insulating surface. However, in all of these cases, the observed structures were completely flat, providing a relatively straightforward path to interpretation. Here, we present high-resolution AFM measurements of several water clusters on Au(111) and Cu(111), whose understanding presents significant challenges due to both their highly 3D configuration and their large size. For each of them, we use a combination of machine learning, atomistic modeling with neural network potentials, and statistical sampling to propose an underlying atomic structure, finally comparing its AFM simulated images to the experimental ones. These results provide insights into the early phases of ice formation, which is a ubiquitous phenomenon ranging from biology to astrophysics.

Evidence of Spontaneous Formation of Two-Dimensional Amorphous Clathrates on Superhydrophilic Surfaces

Clathrate hydrates reserved in the seabed are often dispersed in the pores of coarse-grained sediments; hence, their formation typically occurs under nanoconfinement. Herein, we show the first molecular dynamics (MD) simulation evidence of the spontaneous formation of two-dimensional (2D) clathrate hydrates on crystal surfaces without conventional nanoconfinement. The kinetic process of 2D clathrate formation is illustrated via simulated single-molecule deposition. 2D amorphous patterns are observed on various superhydrophilic face-centered cubic surfaces. Notably, the formation of 2D amorphous clathrate can occur over a wide range of temperatures, even at room temperature. The strong water-surface interaction, the characteristic properties of guest-gas molecules, and the underlying surface structure dictate the formation of 2D amorphous clathrates. Semiquantitative phase diagrams of 2D clathrates are constructed where representative patterns of 2D clathrates for characteristic gas molecules on prototypical Pd(111) and Pt(111) surfaces are confirmed by independent MD simulations. A tunable pattern of 2D amorphous clathrates is demonstrated by changing the lattice strain of the underlying substrate. Moreover, ab initio MD simulations confirm the stability of 2D amorphous clathrate. The underlining physical mechanism for 2D clathrate formation on superhydrophilic surfaces is elucidated, which offers deeper insight into the crucial role of water-surface interaction.

Identification of a common ice nucleus on hydrophilic and hydrophobic close-packed metal surfaces

Establishing a general model of heterogeneous ice nucleation has long been challenging because of the surface water structures found on different substrates. Identifying common water clusters, regardless of the underlying substrate, is one of the key steps toward solving this problem. Here, we demonstrate the presence of a common water cluster found on both hydrophilic Pt(111) and hydrophobic Cu(111) surfaces using scanning tunneling microscopy and non-contact atomic force microscopy. Water molecules self-assemble into a structure with a central flat-lying hexagon and three fused pentagonal rings, forming a cluster consisting of 15 individual water molecules. This cluster serves as a critical nucleus during ice nucleation on both surfaces: ice growth beyond this cluster bifurcates to form two-dimensional (three-dimensional) layers on hydrophilic (hydrophobic) surfaces. Our results reveal the inherent similarity and distinction at the initial stage of ice growth on hydrophilic and hydrophobic close-packed metal surfaces; thus, these observations provide initial evidence toward a general model for water-substrate interaction.

Wetting of a Stepped Platinum (211) Surface

Steps stabilize water adsorption on metal surfaces, providing favorable binding sites for water during wetting or ice nucleation, but there is limited understanding of the local water arrangements formed on such surfaces. Here we describe the structural evolution of water on the stepped Pt(211) surface using thermal desorption, low-energy electron diffraction, and scanning tunneling microscopy to probe the water structure. At low coverage water forms linear structures comprising zigzag chains along the steps that are decorated by H-bonded rings every one or two units along the terrace. Simple 2-coordinate H-bonded chains are not observed, indicating the Pt step binds too weakly to compensate entirely for a low water H-bond coordination number. As the coverage increases, water chains assemble into a disordered (2 × 1) structure, likely made up of the same narrow water chains along the steps with little or no H-bonding between adjacent structures. The chain structure disappears as water adsorption saturates the surface to form an incommensurate, disordered network of water rings of different size. Although the steps on Pt(211) clearly stabilize water adsorption and direct growth, the surface does not support the simple 1D chains previously proposed or an ordered 2D network such as seen on other surfaces. We discuss reasons for this and the factors that determine the behavior of the first water layer on stepped metal surfaces.

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.

Water adsorption and dynamics on graphene and other 2D materials: Computational and experimental advances

The interaction of water and surfaces, at molecular level, is of critical importance for understanding processes such as corrosion, friction, catalysis and mass transport. The significant literature on interactions with single crystal metal surfaces should not obscure unknowns in the unique behaviour of ice and the complex relationships between adsorption, diffusion and long-range inter-molecular interactions. Even less is known about the atomic-scale behaviour of water on novel, non-metallic interfaces, in particular on graphene and other 2D materials. In this manuscript, we review recent progress in the characterisation of water adsorption on 2D materials, with a focus on the nano-material graphene and graphitic nanostructures; materials which are of paramount importance for separation technologies, electrochemistry and catalysis, to name a few. The adsorption of water on graphene has also become one of the benchmark systems for modern computational methods, in particular dispersion-corrected density functional theory (DFT). We then review recent experimental and theoretical advances in studying the single-molecular motion of water at surfaces, with a special emphasis on scattering approaches as they allow an unparalleled window of observation to water surface motion, including diffusion, vibration and self-assembly.

Robustness of Bilayer Hexagonal Ice against Surface Symmetry and Corrugation

Two-dimensional (2D) bilayer hexagonal ice (BHI) is regarded as the first intrinsic 2D ice crystal. However, the robustness of such a structure or its derivatives against surface symmetry and corrugation is still unclear. Here, we report the formation of 2D BHI on gold surfaces with 1D corrugation, using noncontact atomic force microscopy. The hexagonal arrangement of the first wetting layer was visualized on the Au(110)-1×2 surface. Upon depositing more water molecules, the first layer would rearrange and shrink, resulting in the formation of buckled BHI. Such a buckled BHI is hydrophobic despite the appearance of dangling OH, due to the strong interlayer bonding. Furthermore, the BHI is also stable on the Au(100)-5×28 surface. This work reveals the unexpected generality of the BHI on corrugated surfaces with nonhexagonal symmetry, thus shedding new light on the microscopic understandings of the low-dimensional ice formation on solid surfaces or under confinement.

Submolecular Insights into Interfacial Water by Hydrogen-Sensitive Scanning Probe Microscopy

ConspectusWater-solid interfaces have attracted extensive attention because of their crucial roles in a wide range of chemical and physical processes, such as ice nucleation and growth, dissolution, corrosion, heterogeneous catalysis, and electrochemistry. To understand these processes, enormous efforts have been made to obtain a molecular-level understanding of the structure and dynamics of water on various solid surfaces. By the use of scanning probe microscopy (SPM), many remarkable structures of H-bonding networks have been directly visualized, significantly advancing our understanding of the delicate competition between water-water and water-solid interactions. Moreover, the detailed dynamics of water molecules, such as diffusion, clustering, dissociation, and intermolecular and intramolecular proton transfer, have been investigated in a well-controlled manner by tip manipulation. However, resolving the submolecular structure of surface water has remained a great challenge for a long time because of the small size and light mass of protons. Discerning the position of hydrogen in water is not only crucial for the accurate determination of the structure of H-bonding networks but also indispensable in probing the proton transfer dynamics and the quantum nature of protons.In this Account, we focus on the recent advances in the H-sensitive SPM technique and its applications in probing the structures, dynamics, and nuclear quantum effects (NQEs) of surface water and ion hydrates at the submolecular level. First, we introduce the development of high-resolution scanning tunneling microscopy/spectroscopy (STM/S) and qPlus-based atomic force microscopy (qPlus-AFM), which allow access to the degrees of freedom of protons in both real and energy space. qPlus-AFM even allows imaging of interfacial water in a weakly perturbative manner by measuring the high-order electrostatic force between the CO-terminated tip and the polar water molecule, which enables the subtle difference of OH directionality to be discerned. Next we showcase the applications of H-sensitive STM/AFM in addressing several key issues related to water-solid interfaces. The surface wetting behavior and the H-bonding structure of low-dimensional ice on various hydrophilic and hydrophobic solid surfaces are characterized at the atomic scale. Then we discuss the quantitative assessment of NQEs of surface water, including proton tunneling and quantum delocalization. Moreover, the weakly perturbative and H-sensitive SPM technique can be also extended to investigations of water-ion interactions on solid surfaces, revealing the effect of hydration structure on the interfacial ion transport. Finally, we provide an outlook on the further directions and challenges for SPM studies of water-solid interfaces.

Advances in Atomic Force Microscopy: Imaging of Two- and Three-Dimensional Interfacial Water

Interfacial water is closely related to many core scientific and technological issues, covering a broad range of fields, such as material science, geochemistry, electrochemistry and biology. The understanding of the structure and dynamics of interfacial water is the basis of dealing with a series of issues in science and technology. In recent years, atomic force microscopy (AFM) with ultrahigh resolution has become a very powerful option for the understanding of the complex structural and dynamic properties of interfacial water on solid surfaces. In this perspective, we provide an overview of the application of AFM in the study of two dimensional (2D) or three dimensional (3D) interfacial water, and present the prospect and challenges of the AFM-related techniques in experiments and simulations, in order to gain a better understanding of the physicochemical properties of interfacial water.

Motion of water monomers reveals a kinetic barrier to ice nucleation on graphene

The interfacial behaviour of water remains a central question to fields as diverse as protein folding, friction and ice formation. While the properties of water at interfaces differ from those in the bulk, major gaps in our knowledge limit our understanding at the molecular level. Information concerning the microscopic motion of water comes mostly from computation and, on an atomic scale, is largely unexplored by experiment. Here, we provide a detailed insight into the behaviour of water monomers on a graphene surface. The motion displays remarkably strong signatures of cooperative behaviour due to repulsive forces between the monomers, enhancing the monomer lifetime ( ≈ 3 s at 125 K) in a free-gas phase that precedes the nucleation of ice islands and, in turn, provides the opportunity for our experiments to be performed. Our results give a molecular perspective on a kinetic barrier to ice nucleation, providing routes to understand and control the processes involved in ice formation.

Metalated Graphyne-Based Networks as Two-Dimensional Materials: Crystallization, Topological Defects, Delocalized Electronic States, and Site-Specific Doping

Graphyne-based two-dimensional (2D) carbon allotropes feature extraordinary physical properties; however, their synthesis as crystalline single-layered materials has remained challenging. We report on the fabrication of large-area organometallic Ag-bis-acetylide networks and their structural and electronic properties on Ag(111) using low-temperature scanning tunneling microscopy combined with density functional theory (DFT) calculations. The metalated graphyne-based networks are robust at room temperature and assembled in a bottom-up approach via surface-assisted dehalogenative homocoupling of terminal alkynyl bromides. Large-area networks of several hundred nanometers with topological defects at domain boundaries are obtained due to the Ag-acetylide bonds' reversible nature. The thermodynamically controlled growth mechanism is explained through the direct observation of intermediates, which differ on Ag(111) and Au(111). Scanning tunneling spectroscopy resolved unoccupied states delocalized across the network. The energy of these states can be shifted locally by the attachment of a different number of Br atoms within the network. DFT revealed that free-standing metal-bis-acetylide networks are semimetals with a linear band dispersion around several high-symmetry points, which suggest the presence of Weyl points. These results demonstrate that the organometallic Ag-bis-acetylide networks feature the typical 2D material properties, which make them of great interest for fundamental studies and electronic materials in devices.

Water Dissociation and Hydroxyl Formation on Ni(110)

Nickel is an active catalyst for hydrogenation and re-forming reactions, with the reactions showing a strong dependence on the surface exposed. Here, we describe the mixed hydroxyl-water phases formed during water dissociation on Ni(110) using scanning tunneling microscopy and low-current low-energy electron diffraction. Water dissociation starts between 150 and 180 K as the H-bond structure evolves from linear one-dimensional (1D) chains of intact water into a two-dimensional (2D) network containing short rows of face-sharing hexagonal rings. As further water desorbs, the hexagonal rows adopt a local (2 × 3) arrangement, forming small, disordered domains separated by strain relief features. Decomposition of this phase occurs near 220 K to form linear 1D structures consisting of flat, zigzag water chains, with each water stabilized by donating one H to hydroxyl to form a branched chain structure. The OH-H₂O chains repel each other, with the saturation layer ordering into a (2 0, 1 4) structure that decomposes to OH near 245 K as further water desorbs. The structure of the mixed OH/H₂O phases is discussed and contrasted with those found on the related Cu(110) surface, with the differences attributed to strain in the 2D H-bond network caused by the short Ni lattice spacing and strong bond to OH/H₂O.

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.

A folded ice monolayer

A highly stable ice monolayer with folded structural motifs is predicted by means of a novel tiling method augmented with ab initio calculations. This ice monolayer has every two neighboring water hexamers connected by a water square yet folded into two distinct planes, and is thus coined as a folded ice model. It is in the ground state in a range of water densities from 0.08 to 0.12 Å-2, with a stronger energy preference at a lower water density. Its stability shown by ab initio molecular dynamics simulations can sustain up to a temperature of 100 K. The tiling method also enables the prediction of a family of considerably stable ice monolayers with a variety of puckered structures. These results enrich our knowledge of low-dimensional water structures and pave a way to explore more exotic ice nanostructures under confinements.

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.

Computational Prediction of Novel Ice Phases: A Perspective

Although computational prediction of new ice phases is a niche field in water science, the scientific subject itself is representative of two important areas in physical chemistry, namely, statistical thermodynamics and molecular simulations. The prediction of a variety of novel ice phases has also attracted general public interest since the 1980s. In particular, the prediction of low-dimensional ice phases has gained momentum since the confirmation of a number of low-dimensional "computer ice" phases in the laboratory over the past decade. In this Perspective, the research advancements in computational prediction of novel ice phases over the past few years are reviewed. Particular attention is placed on new ice phases whose physical properties or dimensional structures are distinctly different from conventional bulk ices. Specific topics include the (i) formation of superionic ices, (ii) electrofreezing of water under high pressure and in a high external electric field, (iii) prediction of low-density porous ice at strongly negative pressure, (iv) ab initio computational study of two-dimensional (2D) ice under nanoscale confinement, and (v) 2D ices formed on a solid surface near ambient temperature without nanoscale confinement. Clearly, the formation of most of these novel ice phases demands certain extreme conditions. Ongoing challenges and new opportunities for predicting new ice phases from either classical molecular dynamics simulation or high-level ab initio computation are discussed.

A Helical Monolayer Ice

Understanding water wetting layers on solid surfaces is essential for many natural and industrial processes. Here we find a helical ice monolayer with every six water molecules helically arranged along the normal of the basal plane by performing an intensive structural search based on ab initio calculations. The helical ice is more stable than all previous models of monolayer and bilayer ices in a wide range of water densities both in vacuum and on weakly interacting substrates due to a stronger network of hydrogen bonds enabled by the helical geometry. More compelling is the fact that this model adequately explains a recent experimental ice monolayer grown on graphite in terms of the lattice parameter, water density, and Moiré pattern. The helical character in the new ice model echoes previously reported helical motifs in one-dimensional ice structures and suggests an unexpected capability of hydrogen bonds in driving the surface reconstruction of ice structures.

Formation of Linear Water Chains on Ni(110)

Materials that bind strongly to water structure the contact layer, modifying its chemical and physical properties in a manner that depends on the symmetry and reactivity of the surface. Although detailed models have been developed for several inert surfaces, much less is known about reactive surfaces, particularly those with a symmetry different from that of ice. Here we investigate water adsorption on a rectangular surface, Ni(110), an active re-forming catalyst that interacts strongly with water. Instead of forming a network of H-bonded cyclic rings, water forms flat 1D water chains, leaving half the Ni atoms exposed. Second layer water also follows the surface symmetry, forming chains of alternating pentamer and heptamer rings in preference to an extended 2D structure. This behavior is different from that found on other surfaces studied previously and is driven by the short lattice spacing of the solid and the strength of the Ni-water bond.

Orientational ordering in heteroepitaxial water ice on metal surfaces

Sum frequency generation spectroscopy uncovers the orientational ordering in crystalline ice films of water grown on Pt(111) and Rh(111).

Reflection high energy electron diffraction (RHEED) study of ice nucleation and growth on Ni(111): influences of adspecies and electron irradiation

How interfacial molecular interactions influence nucleation and growth processes of water ice is explored using pristine, oxygenated, and CO-adsorbed Ni(111) substrates based on RHEED, together with the effects of high-energy electron irradiation on the crystallization kinetics. A monolayer of amorphous solid water deposited onto the pristine Ni(111) substrate crystallizes into ice Ic at ca. 150 K, whereas ice Ih (Ic) is formed preferentially during water vapor deposition at 135 K (125 K). The ice nucleation tends to be hampered on the oxygenated Ni(111) surface because of the hydrogen bond formation with chemisorbed oxygen, leading to the growth of randomly-oriented ice Ic crystallites via spontaneous nucleation. The amorphization and recrystallization of initially crystalline ices are observed during prolonged RHEED measurements at 20 and 70 K, respectively, signifying that high-energy electron irradiation has both thermal and non-thermal effects on the water phase transition. The epitaxial growth (non-epitaxial growth) of ice occurs during electron irradiation of amorphous solid water formed on the pristine and oxygenated Ni(111) substrates (CO-adsorbed Ni(111) substrate) even at 100 K (120 K) because nucleation and growth are initiated at the substrate interface (in the ASW film interior).

Direct observation of 2-dimensional ices on different surfaces near room temperature without confinement

Water-solid interfaces play important roles in a wide range of fields, including atmospheric science, geochemistry, electrochemistry, and food science. Herein, we report simulation evidence of 2-dimensional (2D) ice formation on various surfaces and the dependence of the 2D crystalline structure on the hydrophobicity and morphology of the underlying surface. Contrary to the prevailing view that nanoscale confinement is necessary for the 2D liquid-to-bilayer ice transition, we find that the liquid-to-bilayer hexagonal ice (BHI) transition can occur either on a model smooth surface or on model fcc-crystal surfaces with indices of (100), (110), and (111) near room temperature. We identify a critical parameter that characterizes the water-surface interaction, above which the BHI can form on the surface. This critical parameter increases as the temperature increases. Even at temperatures above the freezing temperature of bulk ice (Ih ), we find that BHI can also form on a superhydrophilic surface due to the strong water-surface interaction. The tendency toward the formation of BHI without confinement reflects a proper water-surface interaction that can compensate for the entropy loss during the freezing transition. Furthermore, phase diagrams of 2D ice formation are described on the plane of the adsorption energy versus the fcc lattice constant (Eads-afcc), where 4 monolayer square-like ices are also identified on the fcc model surfaces with distinct water-surface interactions.

Strongly Hole-Doped and Highly Decoupled Graphene on Platinum by Water Intercalation

Scanning tunneling microscopy and spectroscopy experiments under ultrahigh vacuum and low-temperature conditions have been performed on water-intercalated graphene on Pt(111). We find that the confined water layer, with a thickness around 0.35 nm, induces a strong hole doping in graphene, i.e., the Dirac point locates at round 0.64 eV above the Fermi level. This can be explained by the presence of a single "puckered bilayer" of ice-Ih, which has not been experimentally found on bare Pt(111), being confined in between graphene and Pt(111) surface. Moreover, the water intercalation makes graphene highly decoupled from the substrate, allowing us to reveal the intrinsic graphene phonons and double Rydberg series of even and odd symmetry image-potential states. Our work not only demonstrates that the electronic properties of graphene can be tuned by the confined water layer between graphene and the substrate, but also provides a generally applicable method to study the intrinsic properties of graphene as well as of other supported two-dimensional materials.

Strain Relief during Ice Growth on a Hexagonal Template

Heterogeneous ice nucleation at solid surfaces impacts many areas of science, from environmental processes, such as precipitation, to microbiological systems and food processing, but the microscopic mechanisms underpinning nucleation remain unclear. Discussion of ice growth has often focused around the role of the surface in templating the structure of water, forcing the first layer to adopt the registry of the underlying substrate rather than that of ice. To grow a thick ice film, water in the first few ice layers must accommodate this strain, but understanding how this occurs requires detailed molecular-scale information that is lacking. Here we combine scanning tunneling microscopy, low-energy electron diffraction, and work-function measurements with electronic structure calculations to investigate the initial stages of ice growth on a Pt alloy surface, having a lattice spacing 6% larger than ice. Although the first layer of water forms a strictly commensurate hexagonal network, this behavior does not extend to the second layer. Instead, water forms a 2D structure containing extended defect rows made from face-sharing pentamer and octamer rings. The defect rows allow the majority of second-layer water to remain commensurate with the solid surface while compensating lateral strain by increasing the water density close to that of an ice surface. The observation of octamer-pentamer rows in ice films formed on several surfaces suggests that the octamer-pentamer defect motif acts as a flexible strain relief mechanism in thin ice films, providing a mechanism that is not available during the growth of strained films in other materials, such as semiconductors.

Influence of the support on stabilizing local defects in strained monolayer oxide films

Two-dimensional materials with a honeycomb lattice, such as graphene and hexagonal boron nitride, often contain local defects in which the hexagonal elements are replaced by four-, five-, seven-, and eight-membered rings. An example is the Stone-Wales (S-W) defect, where a bond rotation causes four hexagons to be transformed into a cluster of two pentagons and two heptagons. A further series of similar defects incorporating divacancies results in larger structures of non-hexagonal elements. In this paper, we use scanning tunneling microscopy (STM) and density functional theory (DFT) modeling to investigate the structure and energetics of S-W and divacancy defects in a honeycomb (2 × 2) Ti2O3 monolayer grown on an Au(111) substrate. The epitaxial rumpled Ti2O3 monolayer is pseudomorphic and in a state of elastic compression. As a consequence, divacancy defects, which induce tension in freestanding films, relieve the compression in the epitaxial Ti2O3 monolayer and therefore have significantly lower energies when compared with their freestanding counterparts. We find that at the divacancy defect sites there is a local reduction of the charge transfer between the film and the substrate, the rumpling is reduced, and the film has an increased separation from the substrate. Our results demonstrate the capacity of the substrate to significantly influence the energetics, and hence favor vacancy-type defects, in compressively strained 2D materials. This approach could be applied more broadly, for example to tensile monolayers, where vacancy-type defects would be rare and interstitial-type defects might be favored.

Crystallization kinetics of thin water films on Pt(111): effects of oxygen and carbon-monoxide adspecies

This paper describes nucleation, epitaxial growth, and wettability of water on Pt(111) and how they are influenced by oxygen and carbon-monoxide adspecies, based on reflection high energy electron diffraction (RHEED), time-of-flight secondary ion mass spectrometry (TOF-SIMS), and temperature-programmed desorption (TPD). Amorphous solid water deposited onto the pristine Pt(111) substrate crystallizes into ice Ih together with a 2D layer at 150 K, whereas ice Ic (stacking disordered ice or a mixture of ice Ic and Ih) is formed preferentially onto oxygenated Pt(111) (CO-adsorbed Pt(111)) at 155-160 K (150 K). The ice nucleation and epitaxial growth tend to be hampered on the oxygenated Pt(111) surface via hydrogen bond formation with chemisorbed oxygen. The CO-adsorbed Pt(111) surface is hydrophobic, as evidenced by the fact that water forms a complex with CO during evaporation of crystallites at 160-165 K. A disordered 2D layer remains on pristine Pt(111) up to 175 K, whereas an ordered 2D layer exhibiting the (√3 ×√3)R30° structure formed on oxygenated Pt(111) up to 200 K.

Enhanced structural disorder at a nanocrystalline ice surface

Enhanced structural disorder at the surface of nanocrystalline ice is studied by heterodyne-detected sum-frequency generation spectroscopy.

Hexagonal Monolayer Ice without Shared Edges

When adsorbed on solids, water molecules are usually arranged into a honeycomb hydrogen-bond network. Here we report the discovery of a novel monolayer ice built exclusively from water hexamers but without shared edges, distinct from all conventional ice phases. Water grown on graphite crystalizes into a robust monolayer ice after annealing, attaining an exceedingly high density of 0.134  Å^{-2}. Unlike chemisorbed ice on metal surfaces, the ice monolayer can translate and rotate on graphite terraces and grow across steps, confirming its two-dimensional nature. First-principles calculations identify the monolayer ice structure as a robust self-assembly of closely packed water hexamers without edge sharing, whose stability is maintained by maximizing the number of intralayer hydrogen bonds on inert surfaces.

Ice Nucleation on a Corrugated Surface

Heterogeneous ice nucleation is a key process in many environmental and technical fields and is of particular importance in modeling atmospheric behavior and the Earth's climate. Despite an improved understanding of how water binds at solid surfaces, no clear picture has emerged to describe how 3D ice grows from the first water layer, nor what makes a particular surface efficient at nucleating bulk ice. This study reports how water at a corrugated, hydrophilic/hydrophobic surface restructures from a complex 2D network, optimized to match the solid surface, to grow into a continuous ice film. Unlike the water networks formed on plane surfaces, the corrugated Cu(511) surface stabilizes a buckled hexagonal wetting layer containing both hydrogen acceptor and donor sites. First layer water is able to relax into an "icelike" arrangement as further water is deposited, creating an array of donor and acceptor sites with the correct spacing and corrugation to stabilize second layer ice and allow continued commensurate multilayer ice growth. Comparison to previous studies of flat surfaces indicates nanoscale corrugation strongly favors ice nucleation, implying surface corrugation will be an important aspect of the surface morphology on other natural or engineered surfaces.

Nanoscale Structure and Dynamics of Water on Pt and Cu Surfaces from MD Simulations

The interaction of liquid water with Pt(111) is investigated with classical molecular dynamics (MD) simulations, where the forces are determined using the third-generation charge optimized many-body (COMB3) interatomic potential. In cases of sub-monolayer water coverage, the parameterized empirical potential predicts experimentally observed and energetically favorable √37 and √39 reconstructed water structures with "575757" di-interstitial defects. At both sub-monolayer and multilayer water coverages, the structure of the first wetting layer of liquid water on Pt(111) exhibits a characteristic distribution where the molecules form two distinct buckled layers as a result of the interplay between water-metal adsorption and water-water hydrogen bonds. The dynamic spreading rate of water nanodroplets on large Pt surfaces (>200 nm²) characterized by molecular kinetic spreading theory is an order of magnitude slower than the molecular kinetic rate of the same droplet on close-packed Cu surfaces due to variation in molecular distributions at the water-metal interface. These nanoscale MD simulation predictions using the COMB3 interatomic potential demonstrate the capability of capturing both many-body interactions between H₂O and Pt or Cu and hydrogen bonding in liquid water.

Identifying Few-Molecule Water Clusters with High Precision on Au(111) Surface

Revealing the nature of a hydrogen-bond network in water structures is one of the imperative objectives of science. With the use of a low-temperature scanning tunneling microscope, water clusters on a Au(111) surface were directly imaged with molecular resolution by a functionalized tip. The internal structures of the water clusters as well as the geometry variations with the increase of size were identified. In contrast to a buckled water hexamer predicted by previous theoretical calculations, our results present deterministic evidence for a flat configuration of water hexamers on Au(111), corroborated by density functional theory calculations with properly implemented van der Waals corrections. The consistency between the experimental observations and improved theoretical calculations not only renders the internal structures of absorbed water clusters unambiguously, but also directly manifests the crucial role of van der Waals interactions in constructing water-solid interfaces.

Preparation-Dependent Orientation of Crystalline Ice Islands on Ag(111)

We observe the transformation of fractal ice islands grown at 96 K to compact ones annealed at 118 K and compare those to compact islands grown directly at 118 K. The low-temperature grown islands form a four bilayer high wetting layer. The annealing causes a crystallization and reshaping of the islands and a substantial increase in height and roughness in particular at higher coverage. Moreover, it leads to a dewetting of the ice film. The islands grown at the higher temperature show qualitative similarities to the annealed ones at smaller nucleation density. However, their orientation with respect to the surface differs by 30° as compared to the annealed islands.