Grand Canonical Investigation of the Quasi Liquid Layer of Ice: Is It Liquid?

Exploring surface properties and premelting in crystals

Crystal surfaces play a pivotal role in governing various significant processes, such as adsorption, nucleation, wetting, friction, and wear. A fundamental property that influences these processes is the surface free energy, γ. We have directly calculated γ(T) for low-index faces of Lennard-Jones (LJ), germanium, and silicon crystals along their sublimation lines using the computational cleavage technique. Our calculations agree well with experimental values for Si(111) and Ge(111), highlighting the accuracy of the method and models used. For LJ crystals, we identified a premelting onset at Tpm = 0.75Tm, marked by a sharp increase in atom mobility within the second outermost surface layer. Notably, Tpm closely aligned with the endpoint of the LJ melting line at negative pressures, Tend = 0.76Tm. We hypothesize that the emergence and coexistence of a liquid film atop the LJ crystal at Tpm < T < Tm correspond to the metastable melting line under negative pressures experienced by stretched crystal surfaces. Furthermore, our study of thin LJ crystal slabs reveals that premelting-induced failure leads to recrystallization below the homogeneous freezing limit, offering a promising avenue to explore crystal nucleation and growth at extremely deep supercoolings. Finally, no evidence of premelting was detected in the model crystals of Ge and Si, which is consistent with the experimental observations. Overall, our findings offer valuable insights into crystal surface phenomena at the atomic scale.

In-layer inhomogeneity of molecular dynamics in quasi-liquid layers of ice

Quasi-liquid layers (QLLs) are present on the surface of ice and play a significant role in its distinctive chemical and physical properties. These layers exhibit considerable heterogeneity across different scales ranging from nanometers to millimeters. Although the formation of partially ice-like structures has been proposed, the molecular-level understanding of this heterogeneity remains unclear. Here, we examined the heterogeneity of molecular dynamics on QLLs based on molecular dynamics simulations and machine learning analysis of the simulation data. We demonstrated that the molecular dynamics of QLLs do not comprise a mixture of solid- and liquid water molecules. Rather, molecules having similar behaviors form dynamical domains that are associated with the dynamical heterogeneity of supercooled water. Nonetheless, molecules in the domains frequently switch their dynamical state. Furthermore, while there is no observable characteristic domain size, the long-range ordering strongly depends on the temperature and crystal face. Instead of a mixture of static solid- and liquid-like regions, our results indicate the presence of heterogeneous molecular dynamics in QLLs, which offers molecular-level insights into the surface properties of ice.

Hierarchical Structuring of Black Silicon Wafers by Ion-Flow-Stimulated Roughening Transition: Fundamentals and Applications for Photovoltaics

Ion-flow-stimulated roughening transition is a phenomenon that may prove useful in the hierarchical structuring of nanostructures. In this work, we have investigated theoretically and experimentally the surface texturing of single-crystal and multi-crystalline silicon wafers irradiated using ion-beam flows. In contrast to previous studies, ions had relatively low energies, whereas flow densities were high enough to induce a quasi-liquid state in the upper silicon layers. The resulting surface modifications reduced the wafer light reflectance to values characteristic of black silicon, widely used in solar energetics. Features of nanostructures on different faces of silicon single crystals were studied numerically based on the mesoscopic Monte Carlo model. We established that the formation of nano-pyramids, ridges, and twisting dune-like structures is due to the stimulated roughening transition effect. The aforementioned variety of modified surface morphologies arises due to the fact that the effects of stimulated surface diffusion of atoms and re-deposition of free atoms on the wafer surface from the near-surface region are manifested to different degrees on different Si faces. It is these two factors that determine the selection of the allowable "trajectories" (evolution paths) of the thermodynamic system along which its Helmholtz free energy, F, decreases, concomitant with an increase in the surface area of the wafer and the corresponding changes in its internal energy, U (dU>0), and entropy, S (dS>0), so that dF=dU - TdS<0, where T is the absolute temperature. The basic theoretical concepts developed were confirmed in experimental studies, the results of which showed that our method could produce, abundantly, black silicon wafers in an environmentally friendly manner compared to traditional chemical etching.

Controlled Adhesion of Ice-Toward Ultraclean 2D Materials

The scalable 2D device fabrication and integration demand either the large-area synthesis or the post-synthesis transfer of 2D layers. While the direct synthesis of 2D materials on most targeted surfaces remains challenging, the transfer approach from the growth substrate onto the targeted surfaces offers an alternative pathway for applications and integrations. However, the current transfer techniques for 2D materials predominantly involve polymers and organic solvents, which are liable to contaminate or deform the ultrasensitive atomic layers. Here, novel ice-aided transfer and ice-stamp transfer methods are developed, in which water (ice) is the only medium in the entire process. In practice, the adhesion between various 2D materials and ice can be well controlled by temperature. Through such controlled adhesion of ice, it is shown that the new transfer methods can yield ultrahigh quality and exceptional cleanliness in transferred 2D flakes and continuous 2D films, and are applicable for a wide range of substrates. Furthermore, beyond transfer, ice can also be used for cleaning the surfaces of 2D materials at higher temperatures. These novel techniques can enable unprecedented ultraclean 2D materials surfaces and performances, and will contribute to the upcoming technological revolutions associated with 2D materials.

Local Ice-like Structure at the Liquid Water Surface

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

Premelting layer during ice growth: role of clusters

The premelting layer plays an important role in ice growth, but there is a significant gap in our knowledge between the atomistic premelting surface structure and the macroscopic growth mechanism. In this work, using large-scale molecular dynamics simulation, we reveal the existence of clusters on the premelting surface, as an intermediate feature bridging the gap. We show the spontaneous formation and evolution of clusters, and they form a stable distribution determined by the growth rate. We demonstrate how this stable distribution is related to the growth mode of ice, connected by the growth of clusters. We come to a bilayer-by-bilayer growth mode at simulation-reachable high growth rates, but another mechanism, namely "cluster stacking", is speculated to exist at lower growth rates. This work builds a connection between the microscopic structure of the premelting layer and the macroscopic growth of ice, making a step forward toward the full understanding of premelting and ice growth.

Investigating the quasi-liquid layer on ice surfaces: a comparison of order parameters

Ice surfaces are characterized by pre-melted quasi-liquid layers (QLLs), which mediate both crystal growth processes and interactions with external agents. Understanding QLLs at the molecular level is necessary to unravel the mechanisms of ice crystal formation. Computational studies of the QLLs heavily rely on the accuracy of the methods employed for identifying the local molecular environment and arrangements, discriminating between solid-like and liquid-like water molecules. Here we compare the results obtained using different order parameters to characterize the QLLs on hexagonal ice (Ih) and cubic ice (Ic) model surfaces investigated with molecular dynamics (MD) simulations in a range of temperatures. For the classification task, in addition to the traditional Steinhardt order parameters in different flavours, we select an entropy fingerprint and a deep learning neural network approach (DeepIce), which are conceptually different methodologies. We find that all the analysis methods give qualitatively similar trends for the behaviours of the QLLs on ice surfaces with temperature, with some subtle differences in the classification sensitivity limited to the solid-liquid interface. The thickness of QLLs on the ice surface increases gradually as the temperature increases. The trends of the QLL size and of the values of the order parameters as a function of temperature for the different facets may be linked to surface growth rates which, in turn, affect crystal morphologies at lower vapour pressure. The choice of the order parameter can be therefore informed by computational convenience except in cases where a very accurate determination of the liquid-solid interface is important.

Scratch-Healing Behavior of Ice by Local Sublimation and Condensation

We show that the surface of ice is scratch healing: micrometer-deep scratches in the ice surface spontaneously disappear by thermal relaxation on the time scale of roughly an hour. Following the dynamics and comparing it to different mass transfer mechanisms, we find that sublimation from and condensation onto the ice surface is the dominant scratch-healing mechanism. The scratch-healing kinetics shows a strong temperature dependence, following an Arrhenius behavior with an activation energy of ΔE = 58.6 ± 4.6 kJ/mol, agreeing with the proposed sublimation mechanism and at odds with surface diffusion or fluid flow or evaporation-condensation from a quasi-liquid layer.

How ice grows from premelting films and water droplets

Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes. At low saturation, freezing proceeds by terrace spreading, but the motion of the buried solid is conveyed through the liquid to the outer liquid-vapor interface. At higher saturations water droplets condense, a large crater forms below, and freezing proceeds undetectably beneath the droplet. Our approach is a general framework that naturally models freezing close to three phase coexistence and provides a first principle theory of ice growth and melting which may prove useful in the geosciences.

Relating Structure and Ice Nucleation of Mixed Surfactant Systems Relevant to Sea Spray Aerosol

Ice nucleating particles (INPs) influence weather and climate by their effect on cloud phase state. Fatty alcohols present within aerosol particles confer a potentially important source of ice nucleation activity to sea spray aerosol produced in oceanic regions. However, their interactions with other aerosol components and the influence on freezing were previously largely unknown. Here, we report quantitative measurements of fatty alcohols in model sea spray aerosol and examine the relationships between the composition and structure of the surfactants and subphase in the context of these measurements. Deposited mixtures of surfactants retain the ability to nucleate ice, even in fatty acid-dominant compositions. Strong refreezing effects are also observed, where previously frozen water-surfactant samples nucleate more efficiently. Structural sources of refreezing behavior are identified as either kinetically trapped film states or three-dimensional (3D) solid surfactant particles. Salt effects are especially important for surfactant INPs, where high salt concentrations suppress freezing. A simple water uptake model suggests that surfactant-containing aerosol requires either very low salt content or kinetic trapping as solid particles to act as INPs in the atmosphere. These types of INPs could be identified through comparison of different INP instrument responses.

Ice-assisted synthesis of functional nanomaterials: the use of quasi-liquid layers as nanoreactors and reaction accelerators

The discovery of peculiar quasi-liquid layers on ice surfaces marks a major breakthrough in ice-related sciences, as the facile tuning of the reactions and morphologies of substances in contact with these layers make ice-assisted chemistry a low-cost, environmentally benign, and ubiquitous methodology for the synthesis of nanomaterials with improved functionality. Ice-templated synthesis of porous materials offers the appealing features of rapid self-organization and remarkable property changes arising from confinement effects and affords materials that have found a diverse range of applications such as batteries, supercapacitors, and gas separation. Moreover, much attention has been drawn to the acceleration of chemical reactions and transformations on the ice surface due to the freeze concentration effect, fast self-diffusion of surface water, and modulated surface potential energy. Some of these results are related to the accumulation of inorganic contaminants in glaciers and the blockage of natural gas pipelines. As an emerging theme in nanomaterial design, the dimension-controlled synthesis of hybrid materials with unprecedentedly enhanced properties on ice surfaces has attracted much interest. However, a deep understanding of quasi-liquid layer characteristics (and hence, the development of cutting-edge analytical technologies with high surface sensitivity) is required to achieve the current goal of ice-assisted chemistry, namely the preparation of tailor-made materials with the desired properties.

Rounded Layering Transitions on the Surface of Ice

Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers. However, disjoining pressure curves, which closely follow expectations from a renormalized mean field liquid state theory, show that there are no layering phase transitions in the thermodynamic sense along the sublimation line. Instead, we find that transitions at mean field level are rounded due to capillary wave fluctuations. We see signatures that true first order layering transitions could arise at low temperatures, for pressures between the metastable line of water-vapor coexistence and the sublimation line. The extrapolation of the disjoining pressure curve above water-vapor saturation displays a true first order phase transition from a thin to a thick film consistent with experimental observations.

Premelting of ice adsorbed on a rock surface

Considering ice-premelting on a quartz rock surface (i.e. silica) we calculate the Lifshitz excess pressures in a four layer system with rock–ice–water–air.

Structure and water attachment rates of ice in the atmosphere: role of nitrogen

In this work we perform computer simulations of the ice surface in order to elucidate the role of nitrogen in the crystal growth rates and crystal habits of snow in the atmosphere. In pure water vapor at temperatures typical of ice crystal formation in cirrus clouds, we find that basal and primary prismatic facets exhibit a layer of premelted ice, with thickness in the subnanometer range. For partial pressures of 1 bar, well above the expected values in the troposphere, we find that only small amounts of nitrogen are adsorbed. The adsorption takes place onto the premelted surface, and hardly any nitrogen dissolves within the premelting film. The premelting film thickness does not change either. We quantify the resulting change of the ice/vapor surface tension to be in the hundredth of mN m^-1 and find that the structure of the pristine ice surface is not changed in a significant manner. We perform a trajectory analysis of colliding water molecules, and find that the attachment rates from direct ballistic collision are very close to unity irrespective of the nitrogen pressure. Nitrogen is however at sufficient density to deflect a fraction of trajectories with smaller distance than the mean free path. Our results show explicitly that the reported differences in growth rates measured in pure water vapor and a controlled nitrogen atmosphere are not related to a significant disruption of the ice surface due to nitrogen adsorption. On the contrary, we show clearly from our trajectory analysis that nitrogen slows down the crystal growth rates due to collisions between water molecules with bulk nitrogen gas. This clarifies the long standing controversy of the role of inert gases on crystal growth rates and demonstrates their influence is solely related to the diffusion limited flow of water vapor across the gas phase.

The Surface of Ice under Equilibrium and Nonequilibrium Conditions

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.

Fluctuation spectroscopy of surface melting of ice with and without impurities

Water is ubiquitous, and the surface properties of ice have been studied for some time, due to their importance. A liquidlike layer (LLL) is known to exist on ice, below the melting point. We use surface thermal fluctuation spectroscopy to study the LLL, including its thickness, for pure ice, and for ice with impurities. We find that the properties of the LLL are experimentally those of liquid water, with thickness much smaller than previous results. We also find that impurities cause the LLL to be thicker, and quite inhomogeneous, with properties depending on the dopant.

How Do Ice Crystals Grow inside Quasiliquid Layers?

A microscopic understanding of crystal-melt interfaces, inseparably involved in the dynamics of crystallization, is a long-standing challenge in condensed matter physics. Here, using an advanced optical microscopy, we directly visualize growing interfaces between ice basal faces and quasiliquid layers (QLLs) during ice crystal growth. This system serves as a model for studying the molecular incorporation process of the crystal growth from a supercooled melt (the so-called melt growth), often hidden by inevitable latent heat diffusion and/or the extremely high crystal growth rate. We reveal that the growth of basal faces inside QLLs proceeds layer by layer via two-dimensional nucleation of monomolecular islands. We also find that the lateral growth rate of the islands is well described by the Wilson-Frenkel law, taking into account the slowing down of the dynamics of water molecules interfaced with ice. These results clearly indicate that, after averaging surface molecular fluctuations, the layer by layer stacking still survives even at the topmost layer on basal faces, which supports the kink-step-terrace picture even for the melt growth.

Why Is It So Difficult to Identify the Onset of Ice Premelting?

Premelting of ice at temperatures below 0 °C is of fundamental importance for environmental processes. Various experimental techniques have been used to investigate the temperature at which liquid-like water first appears at the ice-vapor interface, reporting onset temperatures from -160 to -2 °C. The signals that identify liquid-like order at the ice-vapor interface in these studies, however, do not show a sharp initiation with temperature. That is at odds with the expected first-order nature of surface phase transitions, and consistent with recent large-scale molecular simulations that show the first premelted layer to be sparse and to develop continuously over a wide range of temperatures. Here we perform a thermodynamic analysis to elucidate the origin of the continuous formation of the first layer of liquid at the ice-vapor interface. We conclude that a negative value of the line tension of the ice-liquid-vapor three-phase contact line is responsible for the continuous character of the transition and the sparse nature of the liquid-like domains in the incomplete first layer.