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

Molecular Fingerprints of Ice Surfaces in Sum Frequency Generation Spectra: A First-Principles Machine Learning Study

Understanding the molecular-level structure and dynamics of ice surfaces is crucial for deciphering several chemical, physical, and atmospheric processes. Vibrational sum-frequency generation (SFG) spectroscopy is the most prominent tool for probing the molecular-level structure of the air-ice interface as it is a surface-specific technique, but the molecular interpretation of SFG spectra is challenging. This study utilizes a machine-learning potential, along with dipole and polarizability models trained on ab initio data, to calculate the SFG spectrum of the air-ice interface. At temperatures below ice surface premelting, our simulations support the presence of a proton-ordered arrangement at the Ice I h surface, similar to that seen in Ice XI. Additionally, our simulations provide insight into the assignment of SFG peaks to specific molecular configurations where possible and assess the contribution of subsurface layers to the overall SFG spectrum. These insights enhance our understanding and interpretation of vibrational studies of environmental chemistry at the ice surface.

Ice Sintering by Sublimation and Condensation

The sintering behavior of ice has been the subject of controversy for more than 160 years. Various factors have led to confusion about the mechanisms behind mass transport during sintering; erroneously derived growth rate exponents, experimental challenges in achieving equilibrium conditions, and incorrect comparisons between ice sintering and snow densification have all played a role. Here we demonstrate that sintering of ice under equilibrium conditions proceeds primarily through sublimation and condensation. Mass transfer occurs through the vapor phase, driven by increased volatility at the formed neck due to its high curvature. Our findings on the sintering of ice spheres are consistent with the healing of micrometer-sized scratches in ice under similar conditions.

Simulating atmospheric freezing of single aqueous droplets to ice in a cryogenically cooled ultrasonic levitator

Atmospheric freezing of water droplets suspended in air followed by cloud formation and precipitation represent fundamental steps of the terrestrial water cycle. These aqueous droplets exhibit distinct freezing mechanisms and thermodynamic requirements compared to bulk water often forming metastable supercooled water at subzero temperatures on the Celsius scale (<273 K) prior to crystallization. Here, we report on a real-time spectroscopic investigation combined with simultaneous visualizations of single aqueous droplet freezing events inside a cryogenically cooled ultrasonic levitation chamber with the ultimate goal of probing the molecular structure evolution and stages of ice formation. The observed droplet freezing follows a pseudoheterogeneous ice nucleation mechanism mimicking the process that occurs for atmospherically supercooled water droplets at the air-water interface. This proof-of-concept experimental setup allows future crystallization studies of homo- and heterogeneously doped aqueous droplets under simulated atmospheric environments-also in the presence of reactive trace gases, thus untangling dynamic molecular interactions and chemical reactions, which are of fundamental interest to low-temperature atmospheric chemistry delineating with ice nucleation mechanisms.

Unveiling the face-dependent ice growth kinetics: Insights from molecular dynamics on the basal and prism surfaces

Ice nucleation and growth are critical in many fields, including atmospheric science, cryobiology, and aviation. However, understanding the detailed mechanisms of ice crystal growth remains challenging. In this work, crystallization at the ice/quasi-liquid layer (QLL) interface of the basal and primary prism (prism1) surfaces of hexagonal ice (Ih) was investigated using molecular dynamics simulations across a wide range of temperatures for the TIP4P/Ice model, with comparisons to the mW coarse-grained model. Together with elucidating the temperature-dependent mechanisms of crystallization, face-specific growth rates were systematically estimated. While the prism surface generally exhibits faster growth rates than the basal surface, a temperature-dependent crossover in growth rates between the basal and prism surfaces is observed in TIP4P/Ice simulations, which correlates with crossovers in QLL thickness and properties and with the well-known column to platelets transition in ice-crystal habits at low vapor pressure. This observation helps decode the complex dependence between crystal morphology and temperature in ice crystals.

Crafting Hollow Spheres via Bulk Ice Melting with ppb-Level Gas Sensing Performance

Ice melting, a common yet complex phenomenon, remains incompletely understood. While theoretical studies suggest that preexisting defects in ice generate "off-lattice" water molecules, triggering bulk ice melting, direct experimental evidence of their form has been lacking as the transparent and transient nature of ice poses significant challenges for observation with current techniques. Here, we introduce an ice-melting-induced lyophilization (IMIL) technique that employs graphene-based nanoprobes to replicate and track liquid evolution within melting bulk ice. Our experimental data and theoretical calculations indicate that "off-lattice" water molecules form spherical droplets that enlarge and coalesce as the melting progresses. Notably, the IMIL technique represents a novel nanotechnology for crafting high-quality hollow spheres by leveraging naturally occurring droplets as templates, offering advantages in simplicity, environmental friendliness, scalability, and size adjustability over traditional methods. Additionally, platinum-loaded graphene-based hollow spheres fabricated via the IMIL technique demonstrate ultrasensitive formaldehyde detection with a 5 parts per billion detection limit, rapid response and recovery times (∼4.9 s), and room-temperature operation without auxiliary technology, outperforming WHO standards and current detection methods. These findings highlight the potential of the IMIL technique for creating versatile hollow spheres for diverse applications.

Molecular dynamics simulation of film water thickness and properties at different interfaces in partially saturated frozen soil systems

The film water, with an exceptional capacity to maintain a premelting, liquid-like state even under subzero conditions, provides a potential dynamic conduit for the movement of water in frozen soils. However, the distinctive structural and dynamic characteristics of film water have not been comprehensively elucidated. In this study, molecular dynamics (MD) simulations were conducted to examine the freezing of a system containing ice, water, silica, and gas. The simulations revealed that as the temperature approaches the melting point, the air-water interface tends to possess a thicker layer of unfrozen water, characterized by a higher diffusion coefficient and lower viscosity. In contrast, the film water near the silica-water interface tends to be thinner and remains relatively unaffected by temperature, with only one twentieth of the diffusion coefficient and nearly 20 times the viscosity observed at the air-water interface. These distinct characteristics resulted from the varying interactions between water molecules and their immediate surroundings. Consequently, the film water in proximity of the silica can be assumed to be relatively immobile compared to that of air-water interface. These findings have implications for the study of unsaturated frozen soil systems, in particular, the importance of considering the film water at the air-water interface in the modeling framework.

Dissecting the Molecular Structure of the Air/Ice Interface from Quantum Simulations of the Sum-Frequency Generation Spectrum

Ice interfaces are pivotal in mediating key chemical and physical processes such as heterogeneous chemical reactions in the environment, ice nucleation, and cloud microphysics. At the ice surface, water molecules form a quasi-liquid layer (QLL) with properties distinct from those of the bulk. Despite numerous experimental and theoretical studies, a molecular-level understanding of the QLL has remained elusive. In this work, we use state-of-the-art quantum dynamics simulations with a realistic data-driven many-body potential to dissect the vibrational sum-frequency generation (vSFG) spectrum of the air/ice interface in terms of contributions arising from individual molecular layers, orientations, and hydrogen-bonding topologies that determine the QLL properties. The agreement between experimental and simulated spectra provides a realistic molecular picture of the evolution of the QLL as a function of the temperature without the need for empirical adjustments. The emergence of specific features in the experimental vSFG spectrum suggests that surface restructuring may occur at lower temperatures. This work not only underscores the critical role of many-body interactions and nuclear quantum effects in understanding ice surfaces but also provides a definitive molecular-level picture of the QLL, which plays a central role in multiphase and heterogeneous processes of relevance to a range of fields, including atmospheric chemistry, cryopreservation, and materials science.

Freeze-Induced Protein Assembly of α-Synuclein into Stable Microspheres to Fabricate Light-Induced Cargo Release Systems

Stable hollow-type microspheres (MSs) have been fabricated using α-synuclein (αS), an amyloidogenic protein, via freeze-induced protein self-assembly. This assembly process involves three steps: rapid freezing to form spherical protein condensates from αS oligomers, frozen annealing to form a crust on the condensate and freeze-drying to create an interior lumen via the three-dimensional (3D) coffee-stain effect. The crust produced during the frozen-annealing step is a β-sheet-mediated protein structure that is presumed to be created at the quasi-liquid layer of the protein-ice interface and thus contributes to the stability of MSs in aqueous solutions at room temperature without any additional surface stabilization. MSs transform into amyloid fibril condensates when heated to 70 °C, and the drug is loaded via centrifugal membrane filtration. Additionally, the MSs were shielded with an iron-alginate layer embedded with gold nanoparticles (AuNPs) to prevent premature leakage and to control drug release. This takes advantage of the photothermal effect of AuNPs, resulting in combined cytotoxicity between the drug and heat. Therefore, drug-loaded MSs comprising αS and AuNPs can be suggested as light-controllable drug delivery systems that exhibit chemical and physical anticancer therapeutic effects.

Capillary-Enhanced Biomimetic Adhesion on Icy Surfaces for High-Performance Antislip Shoe-Soles

The World Health Organization (WHO) reports 684,000 deaths/year due to slips and falls (SFs), with ∼38 million people requiring medical attention per annum. In particular, SFs on ice surfaces account for 45% of all SF incidents, costing over $100 billion globally in healthcare, intensive care, and insurance expenses. Current antislip solutions focus on hydrophobicity to repel interfacial fluids, aiming to maintain solid-to-solid contact. However, these solutions often wear out quickly, clog, or become ineffective. Wet ice is particularly challenging due to its nanometer-thick quasi-liquid layer (QLL), which makes it extremely slippery. Inspired by the capillary suction adhesion observed in gecko footpads and the slip resistance of frog toepads on wet surfaces, we developed an innovative approach to regulate ice adhesion and deadhesion. The solution presented in this work mimics this mechanism by employing textured microcavities into silicone rubber (SR)/zirconia (ZrO2) closely mirroring the properties of gecko and frog toepads. Given the dynamics of walking, the surface exhibited hydrophilicity-induced capillary suction of the QLL, facilitating their rapid frost to achieve greater mechanical interlocking. The developed textures displayed capillary suction within 1.5 ms, resulting in a maximum friction coefficient of 3.46 on wet ice. This breakthrough outcome provides a robust, durable solution to significantly reduce SFs on ice surfaces, saving lives and livelihoods.

Investigating Shear Stress of Ice Accumulated on Surfaces with Various Roughnesses: Effects of a Quasi-Water Layer

The investigation of the anti-icing/deicing is essential because the icing phenomenon deteriorates the natural environment and various projects. By conducting molecular dynamics simulation, this work analyzes the effect of the quasi-water layer on the ice shear stress over smooth and rough surfaces, along with the underlying physics of the quasi-water layer. The results indicate that the thickness of the quasi-water layer monotonically increases with temperature, resulting in a monotonic decrease in the ice shear stress on the smooth surface. Due to the joint effects of the smooth surface wettability and the quasi-water layer, the ice shear stress increases and then decreases to almost a constant value when the surface changes from a hydrophobic to a hydrophilic one. For rough surfaces with stripe nanostructures, when the width of the bump for one case equals the depression for the other case, the variations of shear stress with height for these two cases are almost the same. The rough surface is effective in reducing the ice shear stress compared to the smooth surface due to the thickening of the quasi-water layer. Each molecule in the quasi-water layer and its four nearest neighboring molecules gradually form a tetrahedral ice-like structure along the direction away from the surface. The radial distribution function also shows that the quasi-water layer resembles the liquid water rather than the ice structure. These findings shed light on developing anti-icing and deicing techniques.

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.

Organic Nucleation: Water Rearrangement Reveals the Pathway of Ibuprofen

The organic nucleation of the pharmaceutical ibuprofen is investigated, as triggered by the protonation of ibuprofen sodium salt at elevated pH. The growth and aggregation of nanoscale solution species by Analytical Ultracentrifugation and Molecular Dynamics (MD) simulations is tracked. Both approaches reveal solvated molecules, oligomers, and prenucleation clusters, their size as well as their hydration at different reaction stages. By combining surface-specific vibrational spectroscopy and MD simulations, water interacting with ibuprofen at the air-water interface during nucleation is probed. The results show the structure of water changes upon ibuprofen protonation in response to the charge neutralization. Remarkably, the water structure continues to evolve despite the saturation of protonated ibuprofen at the hydrophobic interface. This further water rearrangement is associated with the formation of larger aggregates of ibuprofen molecules at a late prenucleation stage. The nucleation of ibuprofen involves ibuprofen protonation and their hydrophobic assembly. The results highlight that these processes are accompanied by substantial water reorganization. The critical role of water is possibly relevant for organic nucleation in aqueous environments in general.

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.

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.

Precisely Controlled Polymerization of Two-Dimensional Conducting Polymers in Quasi-Liquid Layer Enables Ultrahigh Sensing Performance

Gas sensors based on conducting polymers offer great potential for high-performance room temperature applications due to their cost-effectiveness, high-sensitivity, and operational advantage. However, their current performance is limited by the deficiency of control in conventional polymerization methods, leading to poor crystallinity and inconsistent material properties. Here, the quasi-liquid layer (QLL) on the ice surface acts as a self-regulating nano-reactor for precise control of thermodynamics and kinetics in the polymerization, resulting in a 7.62 nm thick two-dimensional (2D) polyaniline (PANI) film matching the QLL thickness. The ultra-thin film optimizes the exposure of active sites, enhancing the detection of analyte gases at low concentrations. It is validated by fabricating a chemiresistive gas sensor with the 2D PANI film, demonstrating stable room-temperature detection of ammonia down to 10 ppt in ambient air with an impressive 10% response. This achievement represents the highest sensitivity among sensors of this kind while maintaining excellent selectivity and repeatability. Moreover, the QLL-controlled polymerization strategy offers an alternative route for precise control of the polymerization process for conducting polymers, enabling the creation of advanced materials with enhanced properties.

First-principles spectroscopy of aqueous interfaces using machine-learned electronic and quantum nuclear effects

Vibrational spectroscopy is a powerful approach to visualising interfacial phenomena. However, extracting structural and dynamical information from vibrational spectra is a challenge that requires first-principles simulations, including non-Condon and quantum nuclear effects. We address this challenge by developing a machine-learning enhanced first-principles framework to speed up predictive modelling of infrared, Raman, and sum-frequency generation spectra. Our approach uses machine learning potentials that encode quantum nuclear effects to generate quantum trajectories using simple molecular dynamics efficiently. In addition, we reformulate bulk and interfacial selection rules to express them unambiguously in terms of the derivatives of polarisation and polarisabilities of the whole system and predict these derivatives efficiently using fully-differentiable machine learning models of dielectric response tensors. We demonstrate our framework's performance by predicting the IR, Raman, and sum-frequency generation spectra of liquid water, ice and the water-air interface by achieving near quantitative agreement with experiments at nearly the same computational efficiency as pure classical methods. Finally, to aid the experimental discovery of new phases of nanoconfined water, we predict the temperature-dependent vibrational spectra of monolayer water across the solid-hexatic-liquid phases transition.

Spiers Memorial Lecture: Water at interfaces

In this article we discuss current issues in the context of the four chosen subtopics for the meeting: dynamics and nano-rheology of interfacial water, electrified/charged aqueous interfaces, ice interfaces, and soft matter/water interfaces. We emphasize current advances in both theory and experiment, as well as important practical manifestations and areas of unresolved controversy.

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.

Ice-confined synthesis of highly ionized 3D-quasilayered polyamide nanofiltration membranes

Existing polyamide (PA) membrane synthesis protocols are underpinned by controlling diffusion-dominant liquid-phase reactions that yield subpar spatial architectures and ionization behavior. We report an ice-confined interfacial polymerization strategy to enable the effective kinetic control of the interfacial reaction and thermodynamic manipulation of the hexagonal polytype (Ih) ice phase containing monomers to rationally synthesize a three-dimensional quasilayered PA membrane for nanofiltration. Experiments and molecular simulations confirmed the underlying membrane formation mechanism. Our ice-confined PA nanofiltration membrane features high-density ionized structure and exceptional transport channels, realizing superior water permeance and excellent ion selectivity.

Pursuing colloidal diamonds

The endeavor to selectively fabricate a cubic diamond is challenging due to the formation of competing phases such as its hexagonal polymorph or others possessing similar free energy. The necessity to achieve this is of paramount importance since the cubic diamond is the only polymorph exhibiting a complete photonic bandgap, making it a promising candidate in view of photonic applications. Herein, we demonstrate that due to the presence of an external field and delicate manipulation of its strength we can attain selectivity in the formation of a cubic diamond in a one-component system comprised of designer tetrahedral patchy particles. The driving force of such a phenomenon is the structure of the first adlayer which is commensurate with the (110) face of the cubic diamond. Moreover, after a successful nucleation event, once the external field is turned off, the structure remains stable, paving an avenue for further post-synthetic treatment.

Mass Accommodation of Water on Ice at Environmentally Relevant Temperatures: Insights from Fast Scanning Calorimetry

Using a conceptually simple, quasi-adiabatic, fast scanning calorimetry technique, we have investigated the sublimation kinetics of ice films with thicknesses ranging from 14 to 400 nm at environmentally relevant temperatures, between 223 and 268 K. The technique enables accurate determination of ice sublimation rates into vacuum under the conditions of free molecular flow during rapid yet quasistatic heating. The measured sublimation fluxes yield the vapor pressure of the ice samples, which is indistinguishable from that derived from experiments under near-equilibrium conditions. Thus, in agreement with the microscopic reversibility principle, we conclude that the mass accommodation coefficient of water by ice is unity and temperature-independent in the temperature range of the studies. We discuss these findings in the context of current computational and theoretical research into the chemistry and physics of aqueous interfaces.

Quasi-water layer sandwiched between hexagonal ice and wall and its influences on the ice tensile stress

The presence of a quasi-water/premelting layer at the interface between wall and ice when the temperature (T) is below the melting point was extensively observed in experiments. In this work, molecular dynamics simulations are performed to analyze the underlying physics of the quasi-water layer and the effects of the layer on the ice tensile stress. The results indicate that each molecule and its four nearest neighbours in the quasi-water layer representing an equilibrium structure gradually form a tetrahedral ice-like arrangement from an unstructured liquid-like structure along the direction away from the wall. The average density of the quasi-water layer is equal to or higher than the bulk density of water at T ≥ 240 K or T ≤ 240 K respectively, and reaches 1.155 g cm^-3 at T = 210 K, suggesting a structural correlation with the high-density liquid phase of water. Depending on the temperature and wall wettability, the thickness of the quasi-water layer (H_q) ranges from ∼2 Å to ∼25 Å. For prescribed hydrophilic walls, H_q increases monotonically with temperature, and is almost proportional to(T_m - T)^-1/3, where T_m is the melting temperature of ice. H_q keeps an almost constant value (2 Å) as the temperature increases and rises sharply after passing a threshold temperature of T ≈ 250 K. In the joint effects of the wall wettability and quasi-water layer's thickness, the ice tensile stress decreasing monotonically at a larger temperature shows an upward trend and then reduces to almost a constant value as the wall changes from a hydrophobic to a hydrophilic one. The results reveal the potential development of anti-icing/de-icing techniques by heating the wall or modifying its surface to increase H_q.

Effect of sodium chloride adsorption on the surface premelting of ice

We characterise the structural properties of the quasi-liquid layer (QLL) at two low-index ice surfaces in the presence of sodium chloride (Na⁺/Cl^-) ions by molecular dynamics simulations. We find that the presence of a high surface density of Na⁺/Cl^- pairs changes the surface melting behaviour from step-wise to gradual melting. The ions lead to an overall increase of the thickness and the disorder of the QLL, and to a low-temperature roughening transition of the air-ice interface. The local molecular structure of the QLL is similar to that of liquid water, and the differences between the basal and primary prismatic surface are attenuated by the presence of Na⁺/Cl^- pairs. These changes modify the crystal growth rates of different facets and the solvation environment at the surface of sea-water ice with a potential impact on light scattering and environmental chemical reactions.

Solvents and Stabilization in Ionic Liquid Films

We report the interfacial structures and chemical environments of ionic liquid films as a function of dilution with molecular solvents and over a range of film thicknesses (a few micrometers). Data from spectroscopic ellipsometry and infrared spectroscopy measurements show differences between films comprised of neat ionic liquids, as well as films comprised of ionic liquids diluted with two molecular solvents (water and acetonitrile). While the water-diluted IL films follow thickness trends predicted by the Landau-Levich model, neat IL and IL/MeCN films deviate significantly from predicted behaviors. Specifically, these film thicknesses are far greater than the predicted values, suggesting enhanced intermolecular interactions or other non-Newtonian behaviors not captured by the theory. We correlate film thicknesses with trends in the infrared intensity profiles across film thicknesses and IL-solvent dilution conditions and interpret the changes from expected behaviors as varying amounts of the film volume existing in isotropic (bulk) vs anisotropic (interfacial) states. The hydrogen bonding network of water-diluted ionic liquids is implicated in the agreement of this system with the Landau-Levich model's thickness predictions.

Crystal Nucleation and Growth of Inorganic Ionic Materials from Aqueous Solution: Selected Recent Developments, and Implications

In this review article, selected, latest theoretical, and experimental developments in the field of nucleation and crystal growth of inorganic materials from aqueous solution are highlighted, with a focus on literature after 2015 and on non-classical pathways. A key point is to emphasize the so far underappreciated role of water and solvent entropy in crystallization at all stages from solution speciation through to the final crystal. While drawing on examples from current inorganic materials where non-classical behavior has been proposed, the potential of these approaches to be adapted to a wide-range of systems is also discussed, while considering the broader implications of the current re-assessment of pathways for crystallization. Various techniques that are suitable for the exploration of crystallization pathways in aqueous solution, from nucleation to crystal growth are summarized, and a flow chart for the assignment of specific theories based on experimental observations is proposed.

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.

Insights into Ionic Liquids: From Z-Bonds to Quasi-Liquids

Ionic liquids (ILs) hold great promise in the fields of green chemistry, environmental science, and sustainable technology due to their unique properties, such as a tailorable structure, the various types available, and their environmentally friendly features. On the basis of multiscale simulations and experimental characterizations, two unique features of ILs are as follows: (1) strong coupling interactions between the electrostatic forces and hydrogen bonds, namely in the Z-bond, and (2) the unique semiordered structure and properties of ultrathin films, specifically regarding the quasi-liquid. In accordance with the aforementioned theoretical findings, many cutting-edge applications have been proposed: for example, CO₂ capture and conversion, biomass conversion and utilization, and energy storage materials. Although substantial progress has been made recently in the field of ILs, considerable challenges remain in understanding the nature of and devising applications for ILs, especially in terms of e.g. in situ/real-time observation and highly precise multiscale simulations of the Z-bond and quasi-liquid. In this Perspective, we review recent developments and challenges for the IL research community and provide insights into the nature and function of ILs, which will facilitate future applications.

Photo-generated hydroxyl radicals contribute to the formation of halogen radicals leading to ozone depletion on and within polar stratospheric clouds surface

Polar stratospheric clouds (PSCs), of which the surface is a dynamic liquid water layer and might consist of aqueous HNO₃ and H₂O₂, is a well-known key meteorological condition contributing to the ozone hole in the polar stratosphere. PSCs has been considered to provide abundant surface for the heterogeneous reactions causing the formation of the Cl₂ and HOCl, which are further photolyzed into Cl and ClO radicals leading to the ozone destruction. Here we demonstrated that the sunlight drives the massive and stable production of OH radicals in aqueous HNO₃ and its main photo-induced byproduct HNO₂. We also found that the photo-generated OH radicals in aqueous HNO₃, HNO₂ and H₂O₂ have the remarkable capability to react with the dissolved HCl, Cl^- and Br^- to form halogen radicals. In addition, we observed that the H₂O₂ can react with dissolved HCl and Br^- in darkness to form and release Cl₂ and Br₂ gases, which could further be photolyzed into reactive halogen radicals whenever sunlight is available. All these findings suggest that, except for the well-known heterogeneous reactions, photochemical reactions involving the aqueous HNO₃ and H₂O₂ on and within PSCs surface might constitute another important halogen activation pathway for ozone destruction. This study may shed deeper insights into the mechanism of halogen radicals resulting in ozone depletion in polar stratosphere.

UV Resonance Raman Spectroscopy of weakly hydrogen-bonded water in the liquid phase and on ice and snow surfaces

The hydrogen bond network has a major role in determining the physical and chemical properties of water both in the solid and in the liquid state. In the bulk liquid...

Interior Melting of Rapidly Heated Gold Nanoparticles

Normal melting invariably starts from surfaces or interfaces due to the weaker bonding constraints in these regions. However, we show that melting can be initiated from the interior of gold nanoparticles with high heating rates. We find that melting starts from the surface with the formation of a premelting layer, as usual, but that the premelting layer does not extend to the interior under certain conditions. Instead, liquid nucleation occurs in the core of the nanoparticle. This unexpected interior melting is connected to the slower melting kinetics, which is related to heat transfer near the premelted surface. The required conditions for interior melting are a suitable size of the nanoparticle and a sufficiently fast heating rate. The present results point to a novel melting regime in nanoparticles. We note that the time scales are now accessible using ultrafast tools such as X-ray lasers that can probe dynamical structure changes, suggesting opportunities for experiments.

Low Density Interior in Supercooled Aqueous Nanodroplets Expels Ions to the Subsurface

The interaction between water and ions within droplets plays a key role in the chemical reactivity of atmospheric and man-made aerosols. Here we report direct computational evidence that in supercooled aqueous nanodroplets a lower density core of tetrahedrally coordinated water expels the cosmotropic ions to the denser and more disordered subsurface. In contrast, at room temperature, depending on the nature of the ion, the radial distribution in the droplet core is nearly uniform or elevated toward the center. We analyze the spatial distribution of a single ion in terms of a reference electrostatic model. The energy of the system in the analytical model is expressed as the sum of the electrostatic and surface energy of a deformable droplet. The model predicts that the ion is subject to a harmonic potential centered at the droplet's center of mass. We name this effect "electrostatic confinement". The model's predictions are consistent with the simulation findings for a single ion at room temperature but not at supercooling. We anticipate this study to be the starting point for investigating the structure of supercooled (electro)sprayed droplets that are used to preserve the conformations of macromolecules originating from the bulk solution.

Edge premelting of two-dimensional ices

The surface of a three-dimensional ice crystal naturally has a quasi-liquid layer (QLL) at temperatures below its bulk melting point, due to a phenomenon called surface premelting. Here, we show that the edges of a two-dimensional (2D) bilayer hexagonal ice adsorbed on solid surfaces undergo premelting as well, resulting in the formation of quasi-liquid bands (QLBs) at the edges. Our extensive molecular dynamics simulations show that the QLB exhibits structure and dynamics indistinguishable from the bilayer liquid phase, acting as a lower-dimensional analog of the QLL on the bulk ice. We further find that at low temperatures, the width of the QLBs at armchair-type edges of the 2D ice is almost identical to that at zigzag-type edges but becomes far greater than the latter at temperatures near the melting point. The chirality-dependent edge premelting of 2D ices should add an important new ingredient to the heterogeneity of premelting.

Ion Pairing Mediates Molecular Organization Across Liquid/Liquid Interfaces

Liquid/liquid interfaces play a central role in scientific fields ranging from nanomaterial synthesis and soft matter electronics to nuclear waste remediation and chemical separations. This diversity of functions arises from an interface's ability to respond to changing conditions in its neighboring bulk phases. Understanding what drives this interfacial flexibility can provide novel avenues for designing new functional interfaces. However, limiting this progress is an inadequate understanding of the subtle intermolecular and interphase interactions taking place at the molecular level. Here, we use surface-specific vibrational sum frequency generation spectroscopy combined with atomistic molecular dynamics simulations to investigate the self-assembly and structure of model ionic oligomers consisting of an oligodimethylsiloxane (ODMS) tail covalently attached to a positively charged methyl imidazolium (MIM⁺) head group at buried oil/aqueous interfaces. We show how the presence of seemingly innocuous salts can impart dramatic changes to the ODMS tail conformations in the oil phase via specific ion effects and ion-pairing interactions taking place in the aqueous phase. These specific ion interactions are shown to drive enhanced amphiphile adsorption, induce morphological changes, and disrupt emergent hydrogen-bonding structures at the interface. Tuning these interactions allows for independent control over the oligomer structure in the oil phase versus interfacial population changes and represents key mechanistic insight that is needed to control chemical reactions at liquid/liquid interfaces.

Water Mobility in the Interfacial Liquid Layer of Ice/Clay Nanocomposites

At solid/ice interfaces, a premelting layer is formed at temperatures below the melting point of bulk water. However, the structural and dynamic properties within the premelting layer have been a topic of intense debate. Herein, we determined the translational diffusion coefficient Dt of water in ice/clay nanocomposites serving as model systems for permafrost by quasi-elastic neutron scattering. Below the bulk melting point, a rapid decrease of Dt is found for charged hydrophilic vermiculite, uncharged hydrophilic kaolin, and more hydrophobic talc, reaching plateau values below -4 °C. At this temperature, Dt in the premelting layer is reduced up to a factor of two compared to supercooled bulk water. Adjacent to charged vermiculite the lowest water mobility was observed, followed by kaolin and the more hydrophobic talc. Results are explained by the intermolecular water interactions with different clay surfaces and interfacial segregation of the low-density liquid water (LDL) component.

Ice-Inspired Superlubricated Electrospun Nanofibrous Membrane for Preventing Tissue Adhesion

Inspired by the superlubricated surface (SLS) of ice, which consists of an ultrathin and contiguous layer of surface-bound water, we built a SLS on the polycaprolactone (PCL)/poly(2-methacryloxyethylphosphorylcholine) (PMPC) composite nanofibrous membrane via electrospinning under controlled relative humidity (RH). The zwitterionic PMPC on the nanofiber provided a surface layer of bound water, thus generating a hydration lubrication surface. Prepared under 20% RH, electrospun PCL/PMPC nanofibers reached a minimum coefficient of friction (COF) of about 0.12 when the weight ratio of PMPC to PCL was 0.1. At a higher RH, a SLS with an ultralow COF of less than 0.05 was formed on the composite nanofibers. The high stability of the SLS hydration layer on the engineered nanofibrous membrane effectively inhibited fibroblast adhesion and markedly reduced tissue adhesion during tendon repair in vivo. This work demonstrates the great potential of this ice-inspired SLS approach in tissue adhesion-prevention applications.

Photodecay of guaiacol is faster in ice, and even more rapid on ice, than in aqueous solution

Snowpacks contain a wide variety of inorganic and organic compounds, including some that absorb sunlight and undergo direct photoreactions. How the rates of these reactions in, and on, ice compare to rates in water is unclear: some studies report similar rates, while others find faster rates in/on ice. Further complicating our understanding, there is conflicting evidence whether chemicals react more quickly at the air-ice interface compared to in liquid-like regions (LLRs) within the ice. To address these questions, we measured the photodegradation rate of guaiacol (2-methoxyphenol) in various sample types, including in solution, in ice, and at the air-ice interface of nature-identical snow. Compared to aqueous solution, we find modest rate constant enhancements (increases of 3- to 6-fold) in ice LLRs, and much larger enhancements (of 17- to 77-fold) at the air-ice interface of nature-identical snow. Our computational modeling suggests the absorption spectrum for guaiacol red-shifts and increases on ice surfaces, leading to more light absorption, but these changes explain only a small portion (roughly 2 to 9%) of the observed rate constant enhancements in/on ice. This indicates that increases in the quantum yield are primarily responsible for the increased photoreactivity of guaiacol on ice; relative to solution, our results suggest that the quantum yield is larger by a factor of roughly 3-6 in liquid-like regions and 12-40 at the air-ice interface.

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.

Interfacial Vibrational Dynamics of Ice Ih and Liquid Water

Insights into energy flow dynamics at ice surfaces are essential for understanding chemical dynamics relevant to atmospheric and geographical sciences. Here, employing ultrafast surface-specific spectroscopy, we report the interfacial vibrational dynamics of ice I_h. A comparison to liquid water surfaces reveals accelerated vibrational energy relaxation and dissipation at the ice surface for hydrogen-bonded OH groups. In contrast, free-OH groups sticking into the vapor phase exhibit substantially slower vibrational dynamics on ice. The acceleration and deceleration of vibrational dynamics of these different OH groups at the ice surface are attributed to enhanced intermolecular coupling and reduced rotational mobility, respectively. Our results highlight the unique properties of free-OH groups on ice, putatively linked to the high catalytic activities of ice surfaces.

Surface phase transitions and crystal habits of ice in the atmosphere

With climate modeling predicting a raise of at least 2°C by year 2100, the fate of ice has become a serious concern, but we still do not understand how ice grows (or melts). In the atmosphere, crystal growth rates of basal and prism facets exhibit an enigmatic temperature dependence and crossover up to three times in a range between 0° and -40°. Here, we use large-scale computer simulations to characterize the ice surface and identify a sequence of previously unidentified phase transitions on the main facets of ice crystallites. Unexpectedly, we find that as temperature is increased, the crystal surface transforms from a disordered phase with proliferation of steps to a smooth phase with small step density. This causes the anomalous increase of step free energies and provides the long sought explanation for the enigmatic crossover of snow crystal growth rates found in the atmosphere.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Direct Experimental Evidence for Markedly Enhanced Surface Proton Activity Inherent to Water Ice

Autoionization and subsequent proton transfer processes determine the proton activity inherent to water molecular systems. In this study, we provide direct experimental evidence that the proton activity is markedly enhanced at the surface of crystalline ice, on the basis of the simultaneous observation of H/D exchange of water molecules at the surface and in the interior of well-defined double-layer ice films composed of H₂O and D₂O. Thermal desorption mass spectrometry showed clear signatures derived from the surface H/D exchange equilibrium, whereas infrared absorption spectroscopy indicated no appreciable H/D exchange progress in the interior. Detailed kinetic analyses revealed that the rate of H/D exchange at the surface is at least 3 orders of magnitude higher than in the interior. This drastic enhancement of the proton activity suggests an extremely high concentration of surface-hydrated protons in comparison with those in the bulk. Our results also highlight the impact of the local hydrogen-bond structure on the autoionization of water molecules.

Multiscale coupled Maxwell's equations and polarizable molecular dynamics simulation based on charge response kernel model

A computational scheme of coupled Maxwell's equations and polarizable molecular dynamics simulation has been developed based on a multi-scale model to describe the coupled dynamics of light electromagnetic waves and molecules in crystalline solids, where the charge response kernel model is employed to incorporate electronic polarization of the molecules. The method is applicable to electronically non-resonant light-matter interaction systems that involve atomic motions in spectroscopy and photonics. Since the scheme simultaneously traces the light propagation in a medium on a macroscopic scale and the microscopic molecular motion under the light electric field, this enables us to treat the experimental setup and mimic its measurement process. As the first applications, we demonstrate three numerical examples of basic spectroscopies of an ice crystalline solid: simulations of reflection and transmission of visible light, infrared absorption measurement, and stimulated Raman scattering measurement. These examples show the detailed behaviors of the interacting light fields and molecules in the spectroscopic processes.

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.