How ice grows from premelting films and water droplets

Freezing of water and melting of ice: theoretical modeling at the nanoscale

Freezing of water and melting of ice at the nanoscale play critical roles in science and technology fields, including aviation systems, infrastructures, and other broad spectrum of technologies. To cope with the icing challenge, nanoscale anti-icing surface technology has been developed. The freezing and melting temperatures can be tailored by manipulating the size (the radius of water or ice); however, it lacks systemic research. In this work, the size effect on the melting temperature of ice nanocrystals was first established, which considered the variation of bond energy and equivalent heat energy from the perspective of the force-heat equivalence energy density principle. Based on the heterogeneous nucleation mode and by further considering the size and temperature effects on the interface energy involved solid-liquid energy and liquid-vapor energy as well as the above developed melting temperature model, another model is established to accurately predict the freezing temperature of water nanodroplets. The parameters required by the two models established in this paper have a clear physical meaning and establish the quantitative relationships among freezing temperature, melting temperature, surface stress, interface energy, and other thermodynamic parameters. The agreement between model prediction and experimental simulation data confirms the validity and universality of the established models. The higher prediction accuracy of this work compared to the other theoretical models, due to the more detailed consideration and the reference point, captures the errors introduced by the experiment or simulation. This study contributes to a deeper understanding of the underlying mechanism of freezing of water and melting of ice nanocrystals and provides theoretical guidance for the design of cryopreservation systems and anti-icing systems for aviation.

Surface premelting of ice far below the triple point

Premelting of ice, a quasi-liquid layer (QLL) at the surface below the melting temperature, was first postulated by Michael Faraday 160 y ago. Since then, it has been extensively studied theoretically and experimentally through many techniques. Existing work has been performed predominantly on hexagonal ice, at conditions close to the triple point. Whether the same phenomenon can persist at much lower pressure and temperature, where stacking disordered ice sublimates directly into water vapor, remains unclear. Herein, we report direct observations of surface premelting on ice nanocrystals below the sublimation temperature using transmission electron microscopy (TEM). Similar to what has been reported on hexagonal ice, a QLL is found at the solid-vapor interface. It preferentially decorates certain facets, and its thickness increases as the phase transition temperature is approached. In situ TEM reveals strong diffusion of the QLL, while electron energy loss spectroscopy confirms its amorphous nature. More significantly, the premelting observed in this work is thought to be related to the metastable low-density ultraviscous water, instead of ambient liquid water as in the case of hexagonal ice. This opens a route to understand premelting and grassy liquid state, far away from the normal water triple point.

Time domain self-bending photonic hook beam based on freezing water droplet

Tunable optical devices are of great interest as they offer adjustability to their functions. Temporal optics is a fast-evolving field, which may be useful both for revolutionizing basic research of time-dependent phenomena and for developing full optical devices. With increasing focus on ecological compatibility, bio-friendly alternatives are a key subject matter. Water in its various forms can open up new physical phenomena and unique applications in photonics and modern electronics. Water droplets freezing on cold surfaces are ubiquitous in nature. We propose and demonstrate the effectual generation of time domain self-bending photonic hook (time-PH) beams by using mesoscale freezing water droplet. The PH light bends near the shadow surface of the droplet into large curvature and angles superior to a conventional Airy beam. The key properties of the time-PH (length, curvature, beam waist) can be modified flexibly by changing the positions and curvature of the water-ice interface inside the droplet. Due to the modifying internal structure of freezing water droplets in real time, we showcase the dynamical curvature and trajectory control of the time-PH beams. Compared with the traditional methods, our phase-change- based materials (water and ice) of the mesoscale droplet have advantages of easy fabrication, natural materials, compact structure and low cost. Such PHs may have applications in many fields, including temporal optics and optical switching, microscopy, sensors, materials processing, nonlinear optics, biomedicine, and so on.

Deicing of Sessile Droplets Using Surface Acoustic Waves

Deicing has significant relevance in various applications such as transportation, energy production, and telecommunication. The use of surface acoustic waves (SAWs) is an attractive option for deicing as it offers several advantages such as localized heating, in situ control, low power, and system integration for highly efficient deicing. Here, we report an understanding of the dynamics of deicing of microlitre volume water droplets (1 to 30 μL) exposed to low power (0.3 W) SAW actuation using an interdigitated electrode on a piezoelectric (LiNbO₃) substrate. We study the time variation of the volume of liquid water from the onset of SAW actuation to complete deicing, which takes 2.5 to 35 s depending on the droplet volume. The deicing phenomenon is attributed to acoustothermal heating which is found to be greatly influenced by the loss of ice adhesion with the substrate and the acoustic streaming within the liquid water. Acoustothermal heating inside the droplet is characterized by the temperature distribution inside the droplet using infrared thermography, and acoustic streaming is observed using dye-based optical microscopy. A rapid enhancement in deicing is observed upon the detachment of ice from the substrate and the onset of acoustic streaming, marked by a sudden increase in the liquid water volume, droplet temperature, and heat transfer coefficient. The deicing time is found to increase linearly with droplet volume as observed from experiments and further verified using a theoretical model. Our study provides an improved understanding of the recently introduced SAW-based deicing technique that may open up the avenue for a suitable alternative to standard deicing protocols.

Ice friction at the nanoscale

The origin of ice slipperiness has been a matter of great controversy for more than a century, but an atomistic understanding of ice friction is still lacking. Here, we perform computer simulations of an atomically smooth substrate sliding on ice. In a large temperature range between 230 and 266 K, hydrophobic sliders exhibit a premelting layer similar to that found at the ice/air interface. On the contrary, hydrophilic sliders show larger premelting and a strong increase of the first adsorption layer. The nonequilibrium simulations show that premelting films of barely one-nanometer thickness are sufficient to provide a lubricating quasi-liquid layer with rheological properties similar to bulk undercooled water. Upon shearing, the films display a pattern consistent with lubricating Couette flow, but the boundary conditions at the wall vary strongly with the substrate's interactions. Hydrophobic walls exhibit large slip, while hydrophilic walls obey stick boundary conditions with small negative slip. By compressing ice above atmospheric pressure, the lubricating layer grows continuously, and the rheological properties approach bulk-like behavior. Below 260 K, the equilibrium premelting films decrease significantly. However, a very large slip persists on the hydrophobic walls, while the increased friction on hydrophilic walls is sufficient to melt ice and create a lubrication layer in a few nanoseconds. Our results show that the atomic-scale frictional behavior of ice is a combination of spontaneous premelting, pressure melting, and frictional heating.

Ice friction at the nanoscale

The origin of ice slipperiness has been a matter of great controversy for more than a century, but an atomistic understanding of ice friction is still lacking. Here, we perform computer simulations of an atomically smooth substrate sliding on ice. In a large temperature range between 230 and 266 K, hydrophobic sliders exhibit a premelting layer similar to that found at the ice/air interface. On the contrary, hydrophilic sliders show larger premelting and a strong increase of the first adsorption layer. The nonequilibrium simulations show that premelting films of barely one-nanometer thickness are sufficient to provide a lubricating quasi-liquid layer with rheological properties similar to bulk undercooled water. Upon shearing, the films display a pattern consistent with lubricating Couette flow, but the boundary conditions at the wall vary strongly with the substrate's interactions. Hydrophobic walls exhibit large slip, while hydrophilic walls obey stick boundary conditions with small negative slip. By compressing ice above atmospheric pressure, the lubricating layer grows continuously, and the rheological properties approach bulk-like behavior. Below 260 K, the equilibrium premelting films decrease significantly. However, a very large slip persists on the hydrophobic walls, while the increased friction on hydrophilic walls is sufficient to melt ice and create a lubrication layer in a few nanoseconds. Our results show that the atomic-scale frictional behavior of ice is a combination of spontaneous premelting, pressure melting, and frictional heating.

Stress accumulation by confined ice in a temperature gradient

When materials freeze, they often undergo damage due to ice growth. Although this damage is commonly ascribed to the volumetric expansion of water upon freezing, it is usually driven by the flow of water toward growing ice crystals that feeds their growth. The freezing of this additional water can cause a large buildup of stress. Here, we demonstrate a technique for characterizing this stress buildup with unprecedented spatial resolution. We create a stable ice-water interface in a controlled temperature gradient and measure the deformation of the confining boundary. Analysis of the deformation field reveals stresses applied to the boundary with [Formula: see text](micrometers) spatial resolution. Globally, stresses increase steadily over time as liquid water is transported to more deeply undercooled regions. Locally, stresses increase until ice growth is stalled by the confining stresses. Importantly, we find a strong localization of stresses, which significantly increases the likelihood of damage caused by the presence of ice, even in apparently benign freezing situations. Ultimately, the limiting stress that the ice exerts is proportional to the local undercooling, in accordance with the Clapeyron equation, which describes the equilibrium between a stressed solid and its melt. Our results are closely connected to the condensation pressure during liquid-liquid phase separation and the crystallization pressure for growing crystals. Thus, they are highly relevant in fields ranging from cryopreservation and frost heave to food science, rock weathering, and art conservation.

Soft wetting with (a)symmetric Shuttleworth effect

The wetting of soft polymer substrates brings in multiple complexities when compared with the wetting on rigid substrates. The contact angle of the liquid is no longer governed by Young’s Law, but is affected by the substrate’s bulk and surface deformations. On top of that, elastic interfaces exhibit a surface energy that depends on how much they are stretched—a feature known as the Shuttleworth effect (or as surface-elasticity). Here, we present two models through which we explore the wetting of drops in the presence of a strong Shuttleworth effect. The first model is macroscopic in character and consistently accounts for large deformations via a neo-Hookean elasticity. The second model is based on a mesoscopic description of wetting, using a reduced description of the substrate’s elasticity. While the second model is more empirical in terms of the elasticity, it enables a gradient dynamics formulation for soft wetting dynamics. We provide a detailed comparison between the equilibrium states predicted by the two models, from which we deduce robust features of soft wetting in the presence of a strong Shuttleworth effect. Specifically, we show that the (a)symmetry of the Shuttleworth effect between the ‘dry’ and ‘wet’ states governs horizontal deformations in the substrate. Our results are discussed in the light of recent experiments on the wettability of stretched substrates.

Homogeneous ice nucleation rates for mW and TIP4P/ICE models through Lattice Mold calculations

Water freezing is the most common liquid-to-crystal phase transition on Earth, however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240K. We employ two different water models, mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic non-polarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih, Ic and a stacking mixture of ice Ih/Ic; reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner which contributes to benchmarking the freezing behaviour of two popular water models.

Intermolecular forces at ice and water interfaces: premelting, surface freezing and regelation

Using Lifshitz theory we assess the role of van der Waals forces at interfaces of ice and water. The results are combined with measured structural forces from computer simulations to develop a quantitative model of the surface free energy of premelting films. This input is employed within the framework of wetting theory and allows us to predict qualitatively the behavior of quasi-liquid layer thickness as a function of ambient conditions. Our results emphasize the significance of vapor pressure. The ice vapor interface is shown to exhibit only incomplete premelting, but the situation can shift to a state of complete surface melting above water saturation. The results obtained serve also to assess the role of subsurface freezing at the water-vapor interface, and we show that intermolecular forces favor subsurface ice nucleation only in conditions of water undersaturation. We show ice regelation at ambient pressure may be explained as a process of capillary freezing, without the need to invoke the action of bulk pressure melting. Our results for van der Waals forces are exploited in order to gauge dispersion interactions in empirical point charge models of 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.

Control of phase ordering and elastic properties in phase field crystals through three-point direct correlation

Effects of three-point direct correlation on properties of the phase field crystal (PFC) modeling are examined for the control of various ordered and disordered phases and their coexistence in both three-dimensional and two-dimensional systems. Such effects are manifested via the corresponding gradient nonlinearity in the PFC free-energy functional that is derived from classical density functional theory. Their significant impacts on the stability regimes of ordered phases, phase diagrams, and elastic properties of the system, as compared to those of the original PFC model, are revealed through systematic analyses and simulations. The nontrivial contribution from three-point direct correlation leads to the variation of the critical point of order-disorder transition to which all the phase boundaries in the temperature-density phase diagram converge. It also enables the variation and control of system elastic constants over a substantial range as needed in modeling different types of materials with the same crystalline structure but different elastic properties. The capability of this PFC approach in modeling both solid and soft matter systems is further demonstrated through the effect of three-point correlation on controlling the vapor-liquid-solid coexistence and transitions for body-centered cubic phase and on achieving the liquid-stripe or liquid-lamellar phase coexistence. All these provide a valuable and efficient method for the study of structural ordering and evolution in various types of material systems.

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.

Gradient-dynamics model for liquid drops on elastic substrates

The wetting of soft elastic substrates exhibits many features that have no counterpart on rigid surfaces. Modelling the detailed elastocapillary interactions is challenging, and has so far been limited to single contact lines or single drops. Here we propose a reduced long-wave model that captures the main qualitative features of statics and dynamics of soft wetting, but which can be applied to ensembles of droplets. The model has the form of a gradient dynamics on an underlying free energy that reflects capillarity, wettability and compressional elasticity. With the model we first recover the double transition in the equilibrium contact angles that occurs when increasing substrate softness from ideally rigid towards very soft (i.e., liquid). Second, the spreading of single drops of partially and completely wetting liquids is considered showing that known dependencies of the dynamic contact angle on contact line velocity are well reproduced. Finally, we go beyond the single droplet picture and consider the coarsening for a two-drop system as well as for a large ensemble of drops. It is shown that the dominant coarsening mode changes with substrate softness in a nontrivial way.

Crystal growth of bcc titanium from the melt and interfacial properties: A molecular dynamics simulation study

The crystal growth kinetics and interfacial properties of titanium (Ti) are studied using molecular dynamics computer simulation. The interactions between the Ti atoms are modeled via an embedded atom method potential. First, the free solidification method (FSM) is used to determine the melting temperature T_m at zero pressure where the transition from liquid to body-centered cubic crystal occurs. From the simulations with the FSM, the kinetic growth coefficients are also determined for different orientations of the crystal, analyzing how the coupling to the thermostat affects the estimates of the growth coefficients. At T_m, anisotropic interfacial stiffnesses and free energies as well as kinetic growth coefficients are determined from capillary wave fluctuations. The so-obtained growth coefficients from equilibrium fluctuations and without the coupling of the system to a thermostat agree well with those extracted from the FSM calculations.