Quasi-liquid Layer Research Articles

Experimental and theoretical evidence for bilayer-by-bilayer surface melting of crystalline ice

On the surface of water ice, a quasi-liquid layer (QLL) has been extensively reported at temperatures below its bulk melting point at 273 K. Approaching the bulk melting temperature from below, the thickness of the QLL is known to increase. To elucidate the precise temperature variation of the QLL, and its nature, we investigate the surface melting of hexagonal ice by combining noncontact, surface-specific vibrational sum frequency generation (SFG) spectroscopy and spectra calculated from molecular dynamics simulations. Using SFG, we probe the outermost water layers of distinct single crystalline ice faces at different temperatures. For the basal face, a stepwise, sudden weakening of the hydrogen-bonded structure of the outermost water layers occurs at 257 K. The spectral calculations from the molecular dynamics simulations reproduce the experimental findings; this allows us to interpret our experimental findings in terms of a stepwise change from one to two molten bilayers at the transition temperature.

Thermodynamic origin of surface melting on ice crystals

Since the pioneering prediction of surface melting by Michael Faraday, it has been widely accepted that thin water layers, called quasi-liquid layers (QLLs), homogeneously and completely wet ice surfaces. Contrary to this conventional wisdom, here we both theoretically and experimentally demonstrate that QLLs have more than two wetting states and that there is a first-order wetting transition between them. Furthermore, we find that QLLs are born not only under supersaturated conditions, as recently reported, but also at undersaturation, but QLLs are absent at equilibrium. This means that QLLs are a metastable transient state formed through vapor growth and sublimation of ice, casting a serious doubt on the conventional understanding presupposing the spontaneous formation of QLLs in ice-vapor equilibrium. We propose a simple but general physical model that consistently explains these aspects of surface melting and QLLs. Our model shows that a unique interfacial potential solely controls both the wetting and thermodynamic behavior of QLLs.

Premelting-Induced Smoothening of the Ice-Vapor Interface.

We perform computer simulations of the quasiliquid layer of ice formed at the ice-vapor interface close to the ice Ih-liquid-vapor triple point of water. Our study shows that the two distinct surfaces bounding the film behave at small wavelengths as atomically rough and independent ice-water and water-vapor interfaces. For long wavelengths, however, the two surfaces couple, large scale parallel fluctuations are inhibited, and the ice-vapor interface becomes smooth. Our results could help explain the complex morphology of ice crystallites.

Characteristics of quasi-liquid metal layer in the cumulative jet and barrier contact zone

According to the known experimental data on explosive metal processing, as well as the theoretical and experimental studies on the sintered metal powders, we estimated the minimum pressure and temperature values at which this layer and its thickness are formed. The study used a hypothetical model of forming a quasi-liquid metal layer between the leader of the cumulative jet, and a barrier at the time of its destruction.

Ice Recrystallization in a Solution of a Cryoprotector and Its Inhibition by a Protein: Synchrotron X-Ray Diffraction Study

Ice formation and recrystallization is a key phenomenon in freezing and freeze-drying of pharmaceuticals and biopharmaceuticals. In this investigation, high-resolution synchrotron X-ray diffraction is used to quantify the extent of disorder of ice crystals in binary aqueous solutions of a cryoprotectant (sorbitol) and a protein, bovine serum albumin. Ice crystals in more dilute (10 wt%) solutions have lower level of microstrain and larger crystal domain size than these in more concentrated (40 wt%) solutions. Warming the sorbitol–water mixtures from 100 to 228 K resulted in partial ice melting, with simultaneous reduction in the microstrain and increase in crystallite size, that is, recrystallization. In contrast to sorbitol solutions, ice crystals in the BSA solutions preserved both the microstrain and smaller crystallite size on partial melting, demonstrating that BSA inhibits ice recrystallization. The results are consistent with BSA partitioning into quasi-liquid layer on ice crystals but not with a direct protein–ice interaction and protein sorption on ice surface. The study shows for the first time that a common (i.e., not-antifreeze) protein can have a major impact on ice recrystallization and also presents synchrotron X-ray diffraction as a unique tool for quantification of crystallinity and disorder in frozen aqueous systems.

Direct Visualization of Quasi-Liquid Layers on Ice Crystal Surfaces Induced by Hydrogen Chloride Gas

Surface melting of ice crystals forms quasi-liquid layers (QLLs) on ice surfaces, and affects a wide variety of natural phenomena. Since QLLs enhance various chemical reactions in ice clouds, the formation of QLLs by atmospheric gases has been studied intensively. However, such studies were performed using spectroscopy techniques, which have low spatial resolution. Here we show the first direct visualization of QLLs on ice basal faces in the presence of hydrogen chloride (HCl) gas (model atmospheric gas) by advanced optical microscopy, which can visualize individual 0.37 nm-thick elementary steps on ice crystal surfaces. We found that the HCl gas induced the appearances of QLLs with a droplet shape in the temperature range from −15.0 to −1.5 °C, where no QLL appears in the absence of HCl gas. This result indicates that HCl gas adsorbed on ice crystal surfaces probably changed the surface structure of ice crystals and then induced the subsequent melting of ice surfaces. We also observed the movement, shape...

Two types of quasi-liquid layers on ice crystals are formed kinetically

Surfaces of ice are covered with thin liquid water layers, called quasi-liquid layers (QLLs), even below their melting point (0 °C), which govern a wide variety of phenomena in nature. We recently found that two types of QLL phases appear that exhibit different morphologies (droplets and thin layers) [Sazaki G. et al. (2012) Proc Natl Acad Sci USA 109(4):1052-1055]. However, revealing the thermodynamic stabilities of QLLs remains a longstanding elusive problem. Here we show that both types of QLLs are metastable phases that appear only if the water vapor pressure is higher than a certain critical supersaturation. We directly visualized the QLLs on ice crystal surfaces by advanced optical microscopy, which can detect 0.37-nm-thick elementary steps on ice crystal surfaces. At a certain fixed temperature, as the water vapor pressure decreased, thin-layer QLLs first disappeared, and then droplet QLLs vanished next, although elementary steps of ice crystals were still growing. These results clearly demonstrate that both types of QLLs are kinetically formed, not by the melting of ice surfaces, but by the deposition of supersaturated water vapor on ice surfaces. To our knowledge, this is the first experimental evidence that supersaturation of water vapor plays a crucially important role in the formation of QLLs.

Insitu Determination of Surface Tension-to-Shear Viscosity Ratio for Quasiliquid Layers on Ice Crystal Surfaces.

We have experimentally determined the surface tension-to-shear viscosity ratio (the so-called characteristic velocity) of quasiliquid layers (QLLs) on ice crystal surfaces from their wetting dynamics. Using an advanced optical microscope, whose resolution reaches the molecular level in the height direction, we directly observed the coalescent process of QLLs and followed the relaxation modes of their contact lines. The relaxation dynamics is known to be governed by the characteristic velocity, which allows us to access the physical properties of QLLs in a noninvasive way. Here we quantitatively demonstrate that QLLs, when completely wetting ices, have a thickness of 9±3 nm and an approximately 200 times lower characteristic velocity than bulk water, whereas QLLs, when partially wetting ices, have a velocity that is 20 times lower than the bulk. This indicates that ice crystal surfaces significantly affect the physical properties of QLLs localized near the surfaces at a nanometer scale.

Molecular Dynamics Simulation of Ice Indentation by Model Atomic Force Microscopy Tips

We have performed extensive molecular dynamics simulations of nanoindentation of an ice slab with model atomic force microscopy (AFM) tips. We found the presence of a quasi-liquid layer between the tip and the ice for all explored indentation depths. For the smallest tip studied (R = 0.55 nm), the force versus indentation depth curves present peaks related to the melting of distinct monolayers of ice, and we were able to calculate the work (free energy) associated with it. For a larger tip (R = 1.80 nm) having a size not commensurate with the average monolayer thickness, we did not find a clear structure in force curves. This work can help guide the interpretation of experimental AFM indentation of ice and other crystalline solids. More specifically, it provides guidelines for tip sizes where layer-by-layer melting can be achieved and for the order of magnitude of forces that need to be detected.

Revealing Surface Waters on an Antifreeze Protein by Fusion Protein Crystallography Combined with Molecular Dynamic Simulations.

Antifreeze proteins (AFPs) adsorb to ice through an extensive, flat, relatively hydrophobic surface. It has been suggested that this ice-binding site (IBS) organizes surface waters into an ice-like clathrate arrangement that matches and fuses to the quasi-liquid layer on the ice surface. On cooling, these waters join the ice lattice and freeze the AFP to its ligand. Evidence for the generality of this binding mechanism is limited because AFPs tend to crystallize with their IBS as a preferred protein-protein contact surface, which displaces some bound waters. Type III AFP is a 7 kDa globular protein with an IBS made up two adjacent surfaces. In the crystal structure of the most active isoform (QAE1), the part of the IBS that docks to the primary prism plane of ice is partially exposed to solvent and has clathrate waters present that match this plane of ice. The adjacent IBS, which matches the pyramidal plane of ice, is involved in protein-protein crystal contacts with few surface waters. Here we have changed the protein-protein contacts in the ice-binding region by crystallizing a fusion of QAE1 to maltose-binding protein. In this 1.9 Å structure, the IBS that fits the pyramidal plane of ice is exposed to solvent. By combining crystallography data with MD simulations, the surface waters on both sides of the IBS were revealed and match well with the target ice planes. The waters on the pyramidal plane IBS were loosely constrained, which might explain why other isoforms of type III AFP that lack the prism plane IBS are less active than QAE1. The AFP fusion crystallization method can potentially be used to force the exposure to solvent of the IBS on other AFPs to reveal the locations of key surface waters.

Structural evolution and wear of ice surface during rubber–ice contact

This study investigates different possible contact mechanisms of rubber–ice sliding friction using in-situ and ex-situ ice microscopy at test temperatures below melting point. Several frictional processes were identified based on the evolution of ice surface morphology. The nature of rubber–ice contact is significantly different for softer and stiffer rubber compounds. For softer compounds the ice microstructure evolves through frost rounding/removal, surface scratching by hard filler plowing, and shearing of quasi-liquid layer by capillary drag. Cold recrystallization and local melting/refreezing of ice was not found to occur owing to homogenous pressure distribution. For stiffer compounds, the pressure distribution is uneven and contact is more localized, resulting in enhanced local deformation features such as subgrain boundaries, visible etch pits and large-scale melting/refreezing.

Melting of crystalline Si nanoparticle investigated by simulation

In the present work, we use molecular dynamics (MD) simulations to investigate melting of the crystalline Si nanoparticle. Atoms in the nanoparticle interact with each other via the Stillinger-Weber potential. Two heating rates are used. We find that melting of the nanoparticle occurs via propagation of quasi-liquid layer from the surface into the core of the nanoparticle until this layer reaches the critical thickness. We find heating rate affects on mechanism of melting of Si nanoparticle, i.e. coexistence of the two melting mechanisms (homogeneous and heterogeneous ones) occurs if low heating rate is used and it is unlike that proposed in the past. Size affects on melting of Si nanoparticle are found and discussed. In addition, we find that the global bond order parameters Q l can be used to detect melting of Si system unlike some calculations presented in the past.

Relationship between densification behavior and stabilization of quasi-liquid grain boundary layers in CuO-doped 0.7CaTiO3–0.3NdAlO3 microwave ceramics

In a combination with aberration-corrected TEM and water quench treatment, CuO–xCu2O grain boundary quasi-liquid layers in multicomponent 0.7CaTiO3–0.3NdAlO3 perovskite microwave ceramics doped by 0.5wt% CuO are observed. Their structure and chemistry are identified. The relationship between their stabilization and densification behavior during the early sintering stage is clarified by constant-heating-rate/isothermal sintering experiments. It is revealed that the densification of CuO-doped 0.7CaTiO3–0.3NdAlO3 ceramics is dominated by viscous flow mechanism. Additionally, the appearance of quasi-liquid layers during the early sintering stage promotes mass transport by short-circuit diffusion. These results can be used as a guide to control the microstructure and properties of microwave ceramics.

Dynamics of ozone and nitrogen oxides at Summit, Greenland. II. Simulating snowpack chemistry during a spring high ozone event with a 1-D process-scale model

Observed depth profiles of nitric oxide (NO), nitrogen dioxide (NO2), and ozone (O3) in snowpack interstitial air at Summit, Greenland were best replicated by a 1-D process-scale model, which included (1) geometrical representation of snow grains as spheres, (2) aqueous-phase chemistry confined to a quasi-liquid layer (QLL) on the surface of snow grains, and (3) initialization of the species concentrations in the QLL through equilibrium partitioning with mixing ratios in snowpack interstitial air. A comprehensive suite of measurements in and above snowpack during a high O3 event facilitated analysis of the relationship between the chemistry of snowpack and the overlying atmosphere. The model successfully reproduced 2 maxima (i.e., a peak near the surface of the snowpack at solar noon and a larger peak occurring in the evening that extended down from 0.5 to 2 m) in the diurnal profile of NO2 within snowpack interstitial air. The maximum production rate of NO2 by photolysis of nitrate (NO3−) was approximately 108 molec cm−3 s−1, which explained daily observations of maxima in NO2 mixing ratios near solar noon. Mixing ratios of NO2 in snowpack interstitial air were greatest in the deepest layers of the snowpack at night and were attributed to thermal decomposition of peroxynitric acid, which produced up to 106 molec NO2 cm−3 s−1. Highest levels of NO in snowpack interstitial air were confined to upper layers of the snowpack and observed profiles were consistent with photolysis of NO2. Production of nitrogen oxides (NOx) from NO3− photolysis was estimated to be two orders of magnitude larger than NO production and supports the hypothesis that NO3− photolysis is the primary source of NOx within sunlit snowpack in the Arctic. Aqueous-phase oxidation of formic acid by O3 resulted in a maximum consumption rate of ∼106–107 molec cm−3 s−1 and was the primary removal mechanism for O3.

Prism and Other High-Index Faces of Ice Crystals Exhibit Two Types of Quasi-Liquid Layers

Surface melting of ice significantly governs a wide variety of phenomena in nature. We recently succeeded in directly observing the surface melting on ice basal faces (Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 1052–1055) by our advanced optical microscopy, which can detect 0.37-nm-thick elementary steps on basal faces. However, the direct observation of surface melting on the other prism and high-index faces remains an experimental challenge. To fully obtain a comprehensive understanding, we need to examine the surface melting on these faces. Here we show the appearance of two types of quasi-liquid layers (QLLs) on prism and high-index faces just below the melting point. We quickly raised the temperature of ice crystals and then observed prism and high-index faces during a roughening process. We found that with increasing temperature, round liquid-like droplets (α-QLLs) first appeared at temperatures higher than −1.4 to −0.5 °C, and then thin liquid-like layers (β-QLLs) also appeared at temperatures higher...

Surface-enhanced nitrate photolysis on ice.

Heterogeneous nitrate photolysis is the trigger for many chemical processes occurring in the polar boundary layer and is widely believed to occur in a quasi-liquid layer (QLL) at the surface of ice. The dipole-forbidden character of the electronic transition relevant to boundary layer atmospheric chemistry and the small photolysis/photoproduct yields in ice (and in water) may confer a significant enhancement and interfacial specificity to this important photochemical reaction at the surface of ice. Using amorphous solid water films at cryogenic temperatures as models for the disordered interstitial air-ice interface within the snowpack suppresses the diffusive uptake kinetics, thereby prolonging the residence time of nitrate anions at the surface of ice. This approach allows their slow heterogeneous photolysis kinetics to be studied, providing the first direct evidence that nitrates adsorbed onto the first molecular layer at the surface of ice are photolyzed more effectively than those dissolved within the bulk. Vibrational spectroscopy allows the ∼3-fold enhancement in photolysis rates to be correlated with the nitrates' distorted intramolecular geometry, thereby hinting at the role played by the greater chemical heterogeneity in their solvation environment at the surface of ice than that in the bulk. A simple 1D kinetic model suggests (1) that a 3(6)-fold enhancement in photolysis rate for nitrates adsorbed onto the ice surface could increase the photochemical NO2 emissions from a 5(8) nm thick photochemically active interfacial layer by 30(60)%, and (2) that 25(40)% of the NO2 photochemical emissions to the snowpack interstitial air are released from the topmost molecularly thin surface layer on ice. These findings may provide a new paradigm for heterogeneous (photo)chemistry at temperatures below those required for a QLL to form at the ice surface.

Computational studies of atmospherically-relevant chemical reactions in water clusters and on liquid water and ice surfaces.

CONSPECTUS: Reactions on water and ice surfaces and in other aqueous media are ubiquitous in the atmosphere, but the microscopic mechanisms of most of these processes are as yet unknown. This Account examines recent progress in atomistic simulations of such reactions and the insights provided into mechanisms and interpretation of experiments. Illustrative examples are discussed. The main computational approaches employed are classical trajectory simulations using interaction potentials derived from quantum chemical methods. This comprises both ab initio molecular dynamics (AIMD) and semiempirical molecular dynamics (SEMD), the latter referring to semiempirical quantum chemical methods. Presented examples are as follows: (i) Reaction of the (NO(+))(NO3(-)) ion pair with a water cluster to produce the atmospherically important HONO and HNO3. The simulations show that a cluster with four water molecules describes the reaction. This provides a hydrogen-bonding network supporting the transition state. The reaction is triggered by thermal structural fluctuations, and ultrafast changes in atomic partial charges play a key role. This is an example where a reaction in a small cluster can provide a model for a corresponding bulk process. The results support the proposed mechanism for production of HONO by hydrolysis of NO2 (N2O4). (ii) The reactions of gaseous HCl with N2O4 and N2O5 on liquid water surfaces. Ionization of HCl at the water/air interface is followed by nucleophilic attack of Cl(-) on N2O4 or N2O5. Both reactions proceed by an SN2 mechanism. The products are ClNO and ClNO2, precursors of atmospheric atomic chlorine. Because this mechanism cannot result from a cluster too small for HCl ionization, an extended water film model was simulated. The results explain ClNO formation experiments. Predicted ClNO2 formation is less efficient. (iii) Ionization of acids at ice surfaces. No ionization is found on ideal crystalline surfaces, but the process is efficient on isolated defects where it involves formation of H3O(+)-acid anion contact ion pairs. This behavior is found in simulations of a model of the ice quasi-liquid layer corresponding to large defect concentrations in crystalline ice. The results are in accord with experiments. (iv) Ionization of acids on wet quartz. A monolayer of water on hydroxylated silica is ordered even at room temperature, but the surface lattice constant differs significantly from that of crystalline ice. The ionization processes of HCl and H2SO4 are of high yield and occur in a few picoseconds. The results are in accord with experimental spectroscopy. (v) Photochemical reactions on water and ice. These simulations require excited state quantum chemical methods. The electronic absorption spectrum of methyl hydroperoxide adsorbed on a large ice cluster is strongly blue-shifted relative to the isolated molecule. The measured and calculated adsorption band low-frequency tails are in agreement. A simple model of photodynamics assumes prompt electronic relaxation of the excited peroxide due to the ice surface. SEMD simulations support this, with the important finding that the photochemistry takes place mainly on the ground state. In conclusion, dynamics simulations using quantum chemical potentials are a useful tool in atmospheric chemistry of water media, capable of comparison with experiment.

Mechanism of anisotropic surface self-diffusivity at the prismatic ice-vapor interface.

Predictive theoretical models for mesoscopic roughening of ice require improved understanding of attachment kinetics occurring at the ice-vapor interface. Here, we use classical molecular dynamics to explore the generality and mechanics of a transition from anisotropic to isotropic self-diffusivity on exposed prismatic surfaces. We find that self-diffusion parallel to the crystallographic a-axis is favored over the c-axis at sub-melt temperatures below about -35 °C, for three different representations of the water-water intermolecular potential. In the low-temperature anisotropic regime, diffusion results from interstitial admolecules encountering entropically distinct barriers to diffusion in the two in-plane directions. At higher temperatures, isotropic self-diffusion occurring deeper within the quasi-liquid layer becomes the dominant mechanism, owing to its larger energy of activation.

Fluctuating Local Recrystallization of Quasi-Liquid Layer of Sub-Micrometer-Scale Ice: A Molecular Dynamics Study

Molecular dynamic simulation of a faceted ice-Ih crystal with the largest dimension of 0.06 µm in a vacuum is performed by employing the TIP4P intermolecular potential at temperatures Tm − 23, Tm − 13, and Tm − 1 K, where Tm is the melting point of the TIP4P bulk ice. We thereby observe at all the temperatures that the quasi-liquid layers (QLLs) formed on the basal (0001) surfaces are bumpy and that the liquid bumps repeatedly form and break at various places in a random manner. At Tm − 1 K, a liquid sheet appears under the liquid bumps. At an intermediate temperature for the bilayer-by-bilayer surface melting, the molecules under the local areas of either thin or thick QLLs have respectively the tendency to melt or recrystallize.

Water accommodation on ice and organic surfaces: insights from environmental molecular beam experiments.

Water uptake on aerosol and cloud particles in the atmosphere modifies their chemistry and microphysics with important implications for climate on Earth. Here, we apply an environmental molecular beam (EMB) method to characterize water accommodation on ice and organic surfaces. The adsorption of surface-active compounds including short-chain alcohols, nitric acid, and acetic acid significantly affects accommodation of D2O on ice. n-Hexanol and n-butanol adlayers reduce water uptake by facilitating rapid desorption and function as inefficient barriers for accommodation as well as desorption of water, while the effect of adsorbed methanol is small. Water accommodation is close to unity on nitric-acid- and acetic-acid-covered ice, and accommodation is significantly more efficient than that on the bare ice surface. Water uptake is inefficient on solid alcohols and acetic acid but strongly enhanced on liquid phases including a quasi-liquid layer on solid n-butanol. The EMB method provides unique information on accommodation and rapid kinetics on volatile surfaces, and these studies suggest that adsorbed organic and acidic compounds need to be taken into account when describing water at environmental interfaces.