Sink or Swim: Ions and Organics at the Ice–Air Interface

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.

Shear-Enhanced Ion Rejection during Seawater Freezing

Ion rejection during seawater freezing is the basis for freeze desalination. A high ion rejection rate is desired for improving the performance of freeze desalination. In this work, we propose a method to enhance the ion rejection rate through external shear, which is demonstrated through molecular dynamics (MD) simulations and experiments. MD simulations show that the ion rejection rate increases with an increasing shear rate. This is attributed to the disruption of the hydration bonds between ions and water molecules in the hydration shell caused by the shear. Consequently, the mobility of ions is increased, and the energy barrier is reduced at the ice-water interface such that ions have a greater chance of diffusing into the aqueous solution, leading to an enhanced ion rejection rate. The MD results in this work are qualitatively confirmed by experiments and provide insights into the enhancement of the ion rejection rate through external parameters.

Charged Ice Particle Beams with Selected Narrow Mass and Kinetic Energy Distributions

Small ice particles play an important role in atmospheric and extraterrestrial chemistry. Circumplanetary ice particles that are encountered by space probes at hypervelocities play a critical role in the determination of surface and subsurface properties of their source bodies. Here we present an apparatus for the generation of low-intensity beams of single mass-selected charged ice particles under vacuum. They are produced via electrospray ionization of water at atmospheric pressure and undergo evaporative cooling when transferred to vacuum through an atmospheric vacuum interface. m/z selection is achieved through two subsequent quadrupole mass filters operated in the variable-frequency mode within a range of m/z values between 8 × 10⁴ and 3 × 10⁷. Velocity and charge of the selected particles are measured using a nondestructive single-pass image charge detector. From the known electrostatic acceleration potentials and settings of the quadrupoles the particle masses could be obtained and be accurately controlled. It has been shown that the droplets are frozen within the transit time of the apparatus such that ice particles are present after the quadrupole stages and finally detected. The demonstrated correspondence between particle mass and specific quadrupole potentials in this device allows preparation of beams of single particles with a repetition rate between 0.1 and 1 Hz with various diameter distributions from 50 to 1000 nm at 30-250 eV of kinetic energy per charge. This corresponds to velocities and particle masses quickly available between 600 m/s (80 nm) and 50 m/s (900 nm) and particle charge numbers (positive) between 10³ and 10⁴[e], depending upon size.

Effect of nitrogen molecules on the growth kinetics at the interface between a (111) plane of cubic ice and water

The molecular-scale growth kinetics of ice from water in the presence of air molecules are still poorly understood, despite their importance for understanding ice particle formation in nature. In this study, a molecular dynamics simulation is conducted to elucidate the molecular-scale growth kinetics at the interface between a (111) plane of cubic ice and water in the presence of N2 molecules. Two potential models of N2 molecules with and without atomic charges are examined. For both models, N2 molecules bind stably to the interface for a period of 1 ns or longer, and the stability of the binding is higher for the charged model than for the noncharged model. Free-energy surfaces of an N2 molecule along the interface and along an ideal (111) plane surface of cubic ice suggest that for both models, the position where an N2 molecule binds stably is different at the interface and on the ideal plane surface, and the stability of the binding is much higher for the interface than for the ideal plane surface. For both models, stacking-disordered ice grows at the interface, and the formation probability of a hexagonal ice layer in the stacking-disordered ice is higher for the charged model than for the uncharged model. The formation probability for the hexagonal ice layer in the stacking-disordered ice depends not only on the stability of binding but also on the positions where N2 molecules bind on the underlying ice, and the number of N2 molecules that bind stably to the underlying ice.

Direct Observation of Anthracene Clusters at Ice Surfaces

Heterogeneous processes can control atmospheric composition. Snow and ice present important, but poorly understood, reaction media that can greatly alter the composition of air in the cryosphere in polar and temperate regions. Atmospheric scientists struggle to reconcile model predictions with field observations in snow-covered regions due in part to experimental challenges associated with monitoring reactions at air-ice interfaces, and debate regarding reaction kinetics and mechanisms has persisted for over a decade. In this work, we use wavelength-resolved fluorescence microscopy to determine the distribution and chemical speciation of the pollutant anthracene at environmentally relevant frozen surfaces. Our results indicate that anthracene adsorbs to frozen surfaces in monomeric form, but that following lateral diffusion, molecules ultimately reside within brine channels at saltwater ice surfaces, and in micron-sized clusters at freshwater ice surfaces; emission profiles indicate extensive self-association. We also measure anthracene photodegradation kinetics in aqueous solution and artificial snow prepared from frozen freshwater and saltwater solutions. Our results suggest that anthracene─and likely other aromatic pollutants─undergo bimolecular photodegradation at the surface of freshwater ice and sea ice, but not at the surface of frozen organic matter. These results will improve predictions of pollutant fate and exposure risk in the cryosphere. The techniques used can be applied to numerous surfaces within and beyond the atmospheric sciences.

Molecular understanding of ion rejection in the freezing of aqueous solutions

In this work, we investigate the microscopic mechanism of ion rejection phenomena during the freezing of aqueous NaCl solutions through molecular dynamics simulations. It is found that the hydration energy for the ion-water interaction is stronger than that between ions and ice, which is the fundamental reason giving rise to the phenomenon of ion rejection. The probability of ions being rejected by ice is determined by the competition between the energy barrier at the ice-water interface and the thermal effect. The ion rejection rate increases with increasing temperature. Furthermore, it is found that the rejection rate of Na+ is higher than that of Cl- because of the relatively large hydration energy difference between Na+-water and Na+-ice interactions. The role of temperature in the applications of ion rejection in freeze desalination is also discussed.

Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation

Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.

Using Excimeric Fluorescence to Study How the Cooling Rate Determines the Behavior of Naphthalenes in Freeze-Concentrated Solutions: Vitrification and Crystallization

We utilized fluorescence spectroscopy to learn about the molecular arrangement of naphthalene (Np) and 1-methylnaphthalene (MeNp) in frozen aqueous solutions. The freezing induces pronounced compound aggregation in the freeze-concentrated solution (FCS) in between the ice grains. The fluorescence spectroscopy revealed prevalent formation of a vitrified solution and minor crystallization of aromatic compounds. The FCS is shown as a specific environment, differing significantly from not only the pure compounds but also the ice surfaces. The results indicate marked disparity between the behavior of the Np and the MeNp; the cooling rate has a major impact on the former but not on the latter. The spectrum of the Np solution frozen at a faster cooling rate (ca 20 K/min) exhibited a temperature-dependent spectral behavior, whereas the spectrum of the solution frozen at a slower rate (ca 2 K/min) did not alter before melting. We interpret the observation through considering the varied composition of the FCS: Fast freezing leads to a higher water content expressed by the plasticizing effect, allowing molecular rearrangement, while slow cooling produces a more concentrated and drier environment. The experiments were conceived as generalizable for environmentally relevant pollutants and human-made freezing.

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

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

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.

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

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

Emerging investigator series: spatial distribution of dissolved organic matter in ice and at air-ice interfaces

Dissolved organic matter (DOM) is a common solute in snow and ice at Earth's surface. Its effects on reaction kinetics in ice and at air-ice interfaces can be large, but are currently difficult to quantify. We used Raman microscopy to characterize the surface and bulk of frozen aqueous solutions containing humic acid, sodium dodecyl sulfate (SDS), and citric acid at a range of concentrations and temperatures. The surface-active species (humic acid and SDS) were distributed differently than citric acid. Humic acid and SDS are almost completely excluded to the air-ice interface during freezing, where they form a film that coats the surface nearly completely. A liquid layer that coats the majority of the surface was observed at all humic acid and SDS concentrations. Citric acid, which is smaller and less surface active, is excluded to liquid channels at the air-ice interface and within the ice bulk, as has previously been reported for ionic solutes such as sodium chloride. Incomplete surface wetting was observed at all citric acid concentrations and at all temperatures (up to -5 °C). Citric acid appears to be solvated in frozen samples, but SDS and humic acid do not. These results will improve our understanding of the effects of organic solutes on environmental and atmospheric chemistry within ice and at air-ice interfaces.

Photodegradation Rate Constants for Anthracene and Pyrene Are Similar in/on Ice and in Aqueous Solution

Snowpacks contain a variety of chemicals, including organic pollutants such as toxic polycyclic aromatic hydrocarbons (PAHs). While PAHs undergo photodegradation in snow and ice, the rates of these reactions remain in debate. Some studies report that photochemical reactions in snow proceed at rates similar to those expected in a supercooled aqueous solution, but other studies report faster reaction rates, particularly at the air-ice interface (i.e., the quasi-liquid layer, or QLL). In addition, one study reported a surprising nonlinear dependence on photon flux. Here we examine the photodegradation of two common PAHs, anthracene and pyrene, in/on ice and in solution. For a given PAH, rate constants are similar in aqueous solution, in internal liquid-like regions of ice, and at the air-ice interface. In addition, we find the expected linear relationship between reaction rate constant and photon flux. Our results indicate that rate constants for the photochemical loss of PAHs in, and on, snow and ice are very similar to those in aqueous solution, with no enhancement at the air-ice interface.

Pre-melting and the adsorption of formic acid at the air-ice interface at 253 K as seen by NEXAFS and XPS

Interactions between trace gases and ice are important in environmental chemistry and for Earth's climate. In particular, the adsorption of trace gases to ice surfaces at temperatures approaching the melting point has raised interest in the past, because of the prevailing pre-melting. Here, we present Near Edge X-ray Absorption Fine Structure (NEXAFS) spectroscopy data at ambient partial pressure of water to better define the onset temperature of pre-melting at the interfacial region of ice. Further, this study directly compares the interaction between an organic acid common in the atmosphere, formic acid, and that of an aliphatic carbon with ice at 253 K. It makes use of X-ray Photoelectron Spectroscopy (XPS) with its inherent narrow probing depth covering both the surface and near-surface bulk region when detecting electrons. We use the tender X-ray range for excitation to locate the organic species within the interfacial region with an extended probing depth compared to published XPS work. Electron kinetic energy dependent C1s photoemission data indicate that, at low coverage of a few 10¹⁴ molecules cm^-2, the presence of formic acid is restricted to the upper ice layers of the interfacial region. Increasing the dosage, formic acid penetrates 6-7 nm into the air-ice interface. The presence of the more hydrophobic aliphatic carbon is restricted to the upper ice monolayers. This direct comparison of an organic acid with an aliphatic compound confirms the emerging picture where solutes enter the interfacial region of ice at a depth related to their specific tendency to form solvation shells.

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

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

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

In this study, the solid-vapor equilibrium and the quasi liquid layer (QLL) of ice Ih exposing the basal and primary prismatic faces were explored by means of grand canonical molecular dynamics simulations with the monatomic mW potential. For this model, the solid-vapor equilibrium was found to follow the Clausius-Clapeyron relation in the range examined, from 250 to 270 K, with a Δ H_sub of 50 kJ/mol in excellent agreement with the experimental value. The phase diagram of the mW model was constructed for the low pressure region around the triple point. The analysis of the crystallization dynamics during condensation and evaporation revealed that, for the basal face, both processes are highly activated, and in particular cubic ice is formed during condensation, producing stacking-disordered ice. The basal and primary prismatic surfaces of ice Ih were investigated at different temperatures and at their corresponding equilibrium vapor pressures. Our results show that the region known as QLL can be interpreted as the outermost layers of the solid where a partial melting takes place. Solid islands in the nanometer length scale are surrounded by interconnected liquid areas, generating a bidimensional nanophase segregation that spans throughout the entire width of the outermost layer even at 250 K. Two approaches were adopted to quantify the QLL and discussed in light of their ability to reflect this nanophase segregation phenomena. Our results in the μVT ensemble were compared with NPT and NVT simulations for two system sizes. No significant differences were found between the results as a consequence of model system size or of the working ensemble. Nevertheless, certain advantages of performing μVT simulations in order to reproduce the experimental situation are highlighted. On the one hand, the QLL thickness measured out of equilibrium might be affected because of crystallization being slower than condensation. On the other, preliminary simulations of AFM indentation experiments show that the tip can induce capillary condensation over the ice surface, enlarging the apparent QLL.

Ice-Nucleating and Antifreeze Proteins Recognize Ice through a Diversity of Anchored Clathrate and Ice-like Motifs

Cold-adapted organisms produce antifreeze and ice-nucleating proteins to prevent and promote ice formation. The crystal structure of hyperactive bacterial antifreeze protein (AFP) MpAFP suggests that this protein binds ice through an anchored clathrate motif. It is not known whether other hyperactive AFPs and ice-nucleating proteins (INPs) use the same motif to recognize or nucleate ice. Here we use molecular simulations to elucidate the ice-binding motifs of hyperactive insect AFPs and a model INP of Pseudomonas syringae. We find that insect AFPs recognize ice through anchored clathrate motifs distinct from that of MpAFP. By performing simulations of ice nucleation by PsINP, we identify two distinct ice-binding sites on opposite sides of the β-helix. The ice-nucleating sequences identified in the simulations agree with those previously proposed for the closely related INP of Pseudomonas borealis based on the structure of the protein. The simulations indicate that these sites have comparable ice-nucleating efficiency, but distinct binding motifs, controlled by the amino acid sequence: one is an anchored clathrate and the other ice-like. We conclude that anchored clathrate and ice-like motifs can be equally effective for binding proteins to ice and promoting ice nucleation.