Excess Hydrogen Bond at the Ice-Vapor Interface around 200 K

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.

Deciphering the Properties of Nanoconfined Aqueous Solutions by Vibrational Sum Frequency Generation Spectroscopy

When confined between walls at nanometer distances, water exhibits surprisingly different properties with reference to bare interfacial water. Based on computer simulations, we demonstrate how vibrational sum frequency generation (VSFG) spectroscopy can be used-even with very mild symmetry breaking-to discriminate multilayer water in wide slit pores from both bilayer and monolayer water confined within molecularly narrow pores. Applying the technique, the VSFG lineshapes of monolayer, bilayer, and multilayer water are found to differ in characteristic ways, which is explained by their distinct density stratifications giving rise to different H-bonding patterns in the respective solvation layers.

Evolution of Chemical Bonding and Crystalline Swelling-Shrinkage of Montmorillonite upon Temperature Changes Probed by in Situ Fourier Transform Infrared Spectroscopy and X-ray Diffraction

Clay minerals are distributed in Earth's crust and troposphere and in Martian crust where temperature varies. Understanding the changes of chemical bonding and crystalline swelling-shrinkage of montmorillonite (Mnt) upon temperature changes is fundamental for studying its surface reactivity and interaction in specific surroundings. However, such an issue remains poorly understood. Here, in situ high- and low-temperature Fourier transform infrared (HT- and LT-FTIR) spectroscopy and X-ray diffraction (HT- and LT-XRD) were performed to study the evolution of chemical bonding and crystalline swelling-shrinkage of sodium-montmorillonite (NaMnt) upon temperature changes. The FTIR results show that the hydroxyl content in NaMnt decreased when the temperature increased from 20 to 700 °C, while it is independent of temperature from 0 to -150 °C. The formation of hydroxyls at the "broken" layer edges of NaMnt is related to adsorbed water molecules on the surfaces, and its content increased when the particle size of the NaMnt decreased. The water molecules in the interlayer space of NaMnt could bond to the tetrahedral sheet of NaMnt through Si₂O-H₂O bonds. HT- and LT-XRD results reveal that all of those water molecules in NaMnt were removed after heating to 100 °C in the heating chamber. The NaMnt was transformed from a state of monolayer interlayer water molecules at 20 °C to a dehydrated state at 100 °C, and then to a dehydroxylated state at 700 °C. Accordingly, the basal spacings of NaMnt were changed from 1.27 to 0.97 nm and then to 0.96 nm, respectively. When NaMnt was cooled from 20 to -268 °C, a crystalline swelling of NaMnt occurred with an increase of 0.03 nm of basal spacing. This work demonstrates that high/low temperature has a remarkable effect on the hydroxyls and the water molecules in NaMnt, which in turn affects its swelling-shrinkage performance. These findings provide some in-depth understanding of the changes of chemical bonding and crystalline swelling-shrinkage of montmorillonite upon temperature changes and the reasons behind these, which might be helpful for the design of engineering Mnt in high-/low-temperature applications.

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.

Accurate molecular orientation at interfaces determined by multimode polarization-dependent heterodyne-detected sum-frequency generation spectroscopy via multidimensional orientational distribution function

Many essential processes occur at soft interfaces, from chemical reactions on aqueous aerosols in the atmosphere to biochemical recognition and binding at the surface of cell membranes. The spatial arrangement of molecules specifically at these interfaces is crucial for many of such processes. The accurate determination of the interfacial molecular orientation has been challenging due to the low number of molecules at interfaces and the ambiguity of their orientational distribution. Here, we combine phase- and polarization-resolved sum-frequency generation (SFG) spectroscopy to obtain the molecular orientation at the interface. We extend an exponentially decaying orientational distribution to multiple dimensions, which, in conjunction with multiple SFG datasets obtained from the different vibrational modes, allows us to determine the molecular orientation. We apply this new approach to formic acid molecules at the air–water interface. The inferred orientation of formic acid agrees very well with ab initio molecular dynamics data. The phase-resolved SFG multimode analysis scheme using the multidimensional orientational distribution thus provides a universal approach for obtaining the interfacial molecular orientation.

Structure of the Water Molecule Layer between Ice and Amorphous/Crystalline Surfaces Based on Molecular Dynamics Simulations

The structure of the water layer between the ice interface and the hydroxylated amorphous/crystalline silica surfaces was investigated using molecular dynamics simulations. The results indicate that the density profile in the direction perpendicular to the surface has two density peaks in the water layer at the ice-silica interface, which are affected by the silanol group density on the wall and the degree of supercooling in the system. In the two density peaks, the one facing the ice interface side has the same structure as the ice crystal, while the other density peak facing the silica surface has an icelike structure. In the solidification process, the ice and icelike structures in the layer progress more on the amorphous silica surface where the density of the silanol groups is low. The relationship between the ice crystallization and the thickness of the layer has been studied in detail; the lower the temperature, the more the ice crystallization progresses and the thinner the layer becomes.

Intrinsic reconstruction of ice-I surfaces

Understanding the precise atomic structure of ice surfaces is critical for revealing the mechanisms of physical and chemical phenomena at the surfaces, such as ice growth, melting, and chemical reactions. Nevertheless, no conclusive structure has been established. In this study, noncontact atomic force microscopy was used to address the characterization of the atomic structures of ice Ih(0001) and Ic(111) surfaces. The topmost hydrogen atoms are arranged with a short-range (2 × 2) order, independent of the ice thickness and growth substrates used. The electrostatic repulsion between non-hydrogen-bonded water molecules at the surface causes a reduction in the number of the topmost hydrogen atoms together with a distortion of the ideal honeycomb arrangement of water molecules, leading to a short-range-ordered surface reconstruction.

Effects of Temperature on Cs+(H2O)20 Clathrate Structure

Clusters consisting of 20 water molecules and a single cesium ion are especially stable due to their clathrate structure that is composed exclusively of three-coordinate water molecules. Clathrate stability was investigated using infrared photodissociation (IRPD) spectroscopy in the free-OH stretching region (∼3600-3800 cm^-1) at ion cell temperatures between 135 and 355 K. At 275 K and colder, IRPD spectra of Cs⁺(H₂O)₂₀ have just one acceptor-acceptor-donor band. At higher temperatures, a higher-energy acceptor-donor band emerges and grows in intensity. Non-clathrate Na⁺(H₂O)₂₀ structures contain both of these bands, which do not change significantly in intensity over the temperature range. These results indicate a rapid onset in the conversion from clathrate to non-clathrate structures with temperature and suggest that some clathrate population remains even at the highest temperatures investigated. These results provide new insights into the role of entropy in clathrate stability.

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.

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.

Tailoring electric field standing waves in reflection-absorption infrared spectroscopy to enhance absorbance from adsorbates on ice surfaces

The spectroscopic detection of molecules adsorbed onto ice surfaces at coverages similar to those encountered under typical environmental conditions requires high surface selectivity and sensitivity that few techniques can afford. An experimental methodology allowing a significant enhancement in the absorbance from adsorbed molecules is demonstrated herein. It exploits Electric Field Standing Wave (EFSW) effects intrinsic to grazing incidence Reflection-Absorption Infrared (RAIR) spectroscopy, where film thickness dependent optical interferences occur between the multiple reflections of the IR beam at the film-vacuum and the substrate-film interfaces. In this case study, CH₄ is used as a probe molecule and is deposited on a 20 ML coverage dense amorphous solid water film adsorbed onto solid Ar underlayers of various thicknesses. We observe that, at thicknesses where destructive interferences coincide with the absorption features from the CH stretching and HCH bending vibrational modes of methane, their intensity increases by a factor ranging from 10 to 25. Simulations of the RAIR spectra of the composite stratified films using a classical optics model reproduce the Ar underlayer coverage dependent enhancements of the absorbance features from CH₄ adsorbed onto the ice surface. They also reveal that the enhancements occur when the square modulus of the total electric field at the film's surface reaches its minimum value. Exploiting the EFSW effect allows the limit of detection to be reduced to a coverage of (0.2 ± 0.2) ML CH₄, which opens up interesting perspectives for spectroscopic studies of heterogeneous atmospheric chemistry at coverages that are more representative of those found in the natural environment.

Sum frequency generation, calculation of absolute intensities, comparison with experiments, and two-field relaxation-based derivation

The experimental sum frequency generation (SFG) spectrum is the response to an infrared pulse and a visible pulse and is a highly surface-sensitive technique. We treat the surface dangling OH bonds at the air/water interface and focus on the absolute SFG intensities for the resonant terms, a focus that permits insight into the consequences of some approximations. For the polarization combinations, the calculated linewidths for the water interface dangling OH SFG band at 3,700 [Formula: see text] are, as usual, too large, because of the customary neglect of motional narrowing. The integrated spectrum is used to circumvent this problem and justified here using a Kubo-like formalism and theoretical integrated band intensities rather than peak intensities. Only relative SFG intensities are usually reported. The absolute integrated SFG intensities for three polarization combinations for sum frequency, visible, and infrared beams are computed. We use molecular dynamics and the dipole and the polarizability matrix elements obtained from infrared and Raman studies of [Formula: see text]O vapor. The theoretical expressions for two of the absolute susceptibilities contain only a single term and agree with experiment to about a factor of 1.3, with no adjustable parameters. The Fresnel factors are included in that comparison. One of the susceptibilities contains instead four positive and negative terms and agrees less well. The expression for the SFG correlation function is normally derived from a statistical mechanical formulation using a time-evolving density matrix. We show how a derivation based on a two-field relaxation leads to the same final result.

Hydrogen order at the surface of ice Ih revealed by vibrational spectroscopy

A combination of heterodyne-detected sum frequency generation spectroscopy and theoretical modeling elucidates that the surface of ice I_h at 100 K has hydrogen order with the OH group pointing upward to the air (“H-up” orientation).

Capturing intrinsic site-dependent spectral signatures and lifetimes of isolated OH oscillators in extended water networks

The extremely broad infrared spectrum of water in the OH stretching region is a manifestation of how profoundly a water molecule is distorted when embedded in its extended hydrogen-bonding network. Many effects contribute to this breadth in solution at room temperature, which raises the question as to what the spectrum of a single OH oscillator would be in the absence of thermal fluctuations and coupling to nearby OH groups. We report the intrinsic spectral responses of isolated OH oscillators embedded in two cold (~20 K), hydrogen-bonded water cages adopted by the Cs⁺·(HDO)(D₂O)₁₉ and D₃O⁺·(HDO)(D₂O)₁₉ clusters. Most OH oscillators yield single, isolated features that occur with linewidths that increase approximately linearly with their redshifts. Oscillators near 3,400 cm^-1, however, occur with a second feature, which indicates that OH stretch excitation of these molecules drives low-frequency, phonon-type motions of the cage. The excited state lifetimes inferred from the broadening are considered in the context of fluctuations in the local electric fields that are available even at low temperature.

Accessing the Accuracy of Density Functional Theory through Structure and Dynamics of the Water-Air Interface

Density functional theory-based molecular dynamics simulations are increasingly being used for simulating aqueous interfaces. Nonetheless, the choice of the appropriate density functional, critically affecting the outcome of the simulation, has remained arbitrary. Here, we assess the performance of various exchange-correlation (XC) functionals, based on the metrics relevant to sum-frequency generation spectroscopy. The structure and dynamics of water at the water-air interface are governed by heterogeneous intermolecular interactions, thereby providing a critical benchmark for XC functionals. We find that the XC functionals constrained by exact functional conditions (revPBE and revPBE0) with the dispersion correction show excellent performance. The poor performance of the empirically optimized density functional (M06-L) indicates the importance of satisfying the exact functional condition. Understanding the performance of different XC functionals can aid in resolving the controversial interpretation of the interfacial water structure and direct the design of novel, improved XC functionals better suited to describing the heterogeneous interactions in condensed phases.

Structure and Dynamics of Water at the Water-Air Interface Using First-Principles Molecular Dynamics Simulations. II. NonLocal vs Empirical van der Waals Corrections

van der Waals (vdW) correction schemes have been recognized to be essential for an accurate description of liquid water in first-principles molecular dynamics simulation. The description of the structure and dynamics of water is governed by the type of the vdW corrections. So far, two vdW correction schemes have been often used: empirical vdW corrections and nonlocal vdW corrections. In this paper, we assess the influence of the empirical vs nonlocal vdW correction schemes on the structure and dynamics of water at the water–air interface. Since the structure of water at the water–air interface is established by a delicate balance of hydrogen bond formation and breaking, the simulation at the water–air interface provides a unique platform to testify as to the heterogeneous interaction of water. We used the metrics [Ohto et al. J. Chem. Theory Comput., 2019, 15, 595−60230468702] which are directly connected with the sum-frequency generation spectroscopic measurement. We find that the overall performance of nonlocal vdW methods is either similar or worse compared to the empirical vdW methods. We also investigated the performance of the optB88-DRSLL functional, which showed slightly less accuracy than the revPBE-D3 method. We conclude that the revPBE-D3 method shows the best performance for describing the interfacial water.

Deconstructing water's diffuse OH stretching vibrational spectrum with cold clusters

The diffuse vibrational envelope displayed by water precludes direct observation of how different hydrogen-bond topologies dictate the spectral response of individual hydroxy group (OH) oscillators. Using cold, isotopically labeled cluster ions, we report the spectral signatures of a single, intact water (H₂O) molecule embedded at various sites in the clathrate-like cage structure adopted by the Cs⁺·(D₂O)₂₀ ion. These patterns reveal the site-dependent correlation between the frequencies of the two OH groups on the same water molecule and establish that the bound OH companion of the free OH group exclusively accounts for bands in the lower-energy region of the spectrum. The observed multiplet structures reveal the homogeneous linewidths of the fundamentals and quantify the anharmonic contributions arising from coupling to both the intramolecular bending and intermolecular soft modes.

The Surface of Ice under Equilibrium and Nonequilibrium Conditions

The ice premelt, often called the quasi-liquid layer (QLL), is key for the lubrication of ice, gas uptake by ice, and growth of aerosols. Despite its apparent importance, in-depth understanding of the ice premelt from the microscopic to the macroscopic scale has not been gained. By reviewing data obtained using molecular dynamics (MD) simulations, sum-frequency generation (SFG) spectroscopy, and laser confocal differential interference contrast microscopy (LCM-DIM), we provide a unified view of the experimentally observed variation in quasi-liquid (QL) states. In particular, we disentangle three distinct types of QL states of disordered layers, QL-droplet, and QL-film and discuss their nature.

The topmost ice layer is energetically unstable, as the topmost interfacial H₂O molecules lose a hydrogen bonding partner, generating a disordered layer at the ice–air interface. This disordered layer is homogeneously distributed over the ice surface. The nature of the disordered layer changes over a wide temperature range from −90 °C to the bulk melting point. Combined MD simulations and SFG measurements reveal that the topmost ice surface starts to be disordered around −90 °C through a process that the topmost water molecules with three hydrogen bonds convert to a doubly hydrogen-bonded species. When the temperature is further increased, the second layer starts to become disordered at around −16 °C. This disordering occurs not in a gradual manner, but in a bilayer-by-bilayer manner.

When the temperature reaches −2 °C, more complicated structures, QL-droplet and QL-film, emerge on the top of the ice surface. These QL-droplets and QL-films are inhomogeneously distributed, in contrast to the disordered layer. We show that these QL-droplet and QL-film emerge only under supersaturated/undersaturated vapor pressure conditions, as partial and pseudopartial wetting states, respectively. Experiments with precisely controlled pressure show that, near the water vapor pressure at the vapor-ice equilibrium condition, no QL-droplet and QL-film can be observed, implying that the QL-droplet and QL-film emerge exclusively under nonequilibrium conditions, as opposed to the disordered layers formed under equilibrium conditions.

These findings are connected with many phenomena related to the ice surface. For example, we explain how the disordering of the topmost ice surface governs the slipperiness of the ice surface, allowing for ice skating. Further focus is on the gas uptake mechanism on the ice surface. Finally, we note the unresolved questions and future challenges regarding the ice premelt.

Fluctuation spectroscopy of surface melting of ice with and without impurities

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

Enhanced structural disorder at a nanocrystalline ice surface

Enhanced structural disorder at the surface of nanocrystalline ice is studied by heterodyne-detected sum-frequency generation spectroscopy.

Molecular Insight into the Slipperiness of Ice

Measurements of the friction coefficient of steel-on-ice over a large temperature range reveal very high friction at low temperatures (-100 °C) and a steep decrease in the friction coefficient with increasing temperature. Very low friction is only found over the limited temperature range typical for ice skating. The strong decrease in the friction coefficient with increasing temperature exhibits Arrhenius behavior with an activation energy of E_a ≈ 11.5 kJ mol^-1. Remarkably, molecular dynamics simulations of the ice-air interface reveal a very similar activation energy for the mobility of surface molecules. Weakly hydrogen-bonded surface molecules diffuse over the surface in a rolling motion, their number and mobility increasing with increasing temperature. This correlation between macroscopic friction and microscopic molecular mobility indicates that slippery ice arises from the high mobility of its surface molecules, making the ice surface smooth and the shearing of the weakly bonded surface molecules easy.

Water at surfaces with tunable surface chemistries

Aqueous interfaces are ubiquitous in natural environments, spanning atmospheric, geological, oceanographic, and biological systems, as well as in technical applications, such as fuel cells and membrane filtration. Where liquid water terminates at a surface, an interfacial region is formed, which exhibits distinct properties from the bulk aqueous phase. The unique properties of water are governed by the hydrogen-bonded network. The chemical and physical properties of the surface dictate the boundary conditions of the bulk hydrogen-bonded network and thus the interfacial properties of the water and any molecules in that region. Understanding the properties of interfacial water requires systematically characterizing the structure and dynamics of interfacial water as a function of the surface chemistry. In this review, we focus on the use of experimental surface-specific spectroscopic methods to understand the properties of interfacial water as a function of surface chemistry. Investigations of the air-water interface, as well as efforts in tuning the properties of the air-water interface by adding solutes or surfactants, are briefly discussed. Buried aqueous interfaces can be accessed with careful selection of spectroscopic technique and sample configuration, further expanding the range of chemical environments that can be probed, including solid inorganic materials, polymers, and water immiscible liquids. Solid substrates can be finely tuned by functionalization with self-assembled monolayers, polymers, or biomolecules. These variables provide a platform for systematically tuning the chemical nature of the interface and examining the resulting water structure. Finally, time-resolved methods to probe the dynamics of interfacial water are briefly summarized before discussing the current status and future directions in studying the structure and dynamics of interfacial water.