Ice plays crucially important roles in various phenomena because of its abundance on Earth. However, revealing the dynamic behavior of quasi-liquid layers (QLLs), which governs the surface properties of ice crystals at temperatures near the melting point, remains an experimental challenge. Here we show that two types of QLL phases appear that exhibit different morphologies and dynamics. We directly visualized the two types of QLLs on ice crystal surfaces by advanced optical microscopy, which can visualize the individual 0.37-nm-thick elementary steps on ice crystal surfaces. We found that they had different stabilities and different interactions with ice crystal surfaces. The two immiscible QLL phases appeared heterogeneously, moved around, and coalesced dynamically on ice crystal surfaces. This picture of surface melting is quite different from the conventional picture in which one QLL phase appears uniformly on ice crystal surfaces.

Due to the abundance of ice on earth, the phase transition of ice plays crucially important roles in various phenomena in nature. Hence, the molecular-level understanding of ice crystal surfaces holds the key to unlocking the secrets of a number of fields. In this study we demonstrate, by laser confocal microscopy combined with differential interference contrast microscopy, that elementary steps (the growing ends of ubiquitous molecular layers with the minimum height) of ice crystals and their dynamic behavior can be visualized directly at air-ice interfaces. We observed the appearance and lateral growth of two-dimensional islands on ice crystal surfaces. When the steps of neighboring two-dimensional islands coalesced, the contrast of the steps always disappeared completely. We were able to discount the occurrence of steps too small to detect directly because we never observed the associated phenomena that would indicate their presence. In addition, classical two-dimensional nucleation theory does not support the appearance of multilayered two-dimensional islands. Hence, we concluded that two-dimensional islands with elementary height (0.37 and 0.39 nm on basal and prism faces, respectively) were visualized by our optical microscopy. On basal and prism faces, we also observed the spiral growth steps generated by screw dislocations. The distance between adjacent spiral steps on a prism face was about 1/20 of that on a basal face. Hence, the step ledge energy of a prism face was 1/20 of that on a basal face, in accord with the known lower-temperature roughening transition of the prism face.

Significance: Phase transitions of ice are a major source of a diverse set of natural phenomena on Earth. In particular, quasi-liquid layers (QLLs) resulting from surface melting are recognized to be key players involved in various natural phenomena spanning from making snowballs to electrification of thunderclouds. With the aid of in situ observations with our advanced optical microscopy combined with two-beam interferometry, we elucidate a thermodynamic origin of the formation of QLLs and their unique wetting behavior (pseudo-partial wetting and wetting transitions) on ice surfaces. We show that QLLs are a metastable transient state formed through vapor growth and sublimation of ice that are absent at equilibrium.

The ice premelt, often called the quasi-liquid layer (QLL), is key for the lubrication of ice, gas uptake by ice, and growth of aerosols. Despite its apparent importance, in-depth understanding of the ice premelt from the microscopic to the macroscopic scale has not been gained. By reviewing data obtained using molecular dynamics (MD) simulations, sum-frequency generation (SFG) spectroscopy, and laser confocal differential interference contrast microscopy (LCM-DIM), we provide a unified view of the experimentally observed variation in quasi-liquid (QL) states. In particular, we disentangle three distinct types of QL states of disordered layers, QL-droplet, and QL-film and discuss their nature. The topmost ice layer is energetically unstable, as the topmost interfacial H₂O molecules lose a hydrogen bonding partner, generating a disordered layer at the ice–air interface. This disordered layer is homogeneously distributed over the ice surface. The nature of the disordered layer changes over a wide temperature range from −90 °C to the bulk melting point. Combined MD simulations and SFG measurements reveal that the topmost ice surface starts to be disordered around −90 °C through a process that the topmost water molecules with three hydrogen bonds convert to a doubly hydrogen-bonded species. When the temperature is further increased, the second layer starts to become disordered at around −16 °C. This disordering occurs not in a gradual manner, but in a bilayer-by-bilayer manner. When the temperature reaches −2 °C, more complicated structures, QL-droplet and QL-film, emerge on the top of the ice surface. These QL-droplets and QL-films are inhomogeneously distributed, in contrast to the disordered layer. We show that these QL-droplet and QL-film emerge only under supersaturated/undersaturated vapor pressure conditions, as partial and pseudopartial wetting states, respectively. Experiments with precisely controlled pressure show that, near the water vapor pressure at the vapor-ice equilibrium condition, no QL-droplet and QL-film can be observed, implying that the QL-droplet and QL-film emerge exclusively under nonequilibrium conditions, as opposed to the disordered layers formed under equilibrium conditions. These findings are connected with many phenomena related to the ice surface. For example, we explain how the disordering of the topmost ice surface governs the slipperiness of the ice surface, allowing for ice skating. Further focus is on the gas uptake mechanism on the ice surface. Finally, we note the unresolved questions and future challenges regarding the ice premelt.