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.

Theoretical and experimental OD-stretch vibrational spectroscopy of heavy water

In view of the current situation in which the OD-stretch vibrational spectra have been scarcely computed with non-polarizable rigid D2O models, we investigate the IR and Raman spectra of D2O by using a newly-reported model TIP4P/2005-HW. From the comparison between the calculations and experimental data, we find the excellent performance of TIP4P/2005-HW for vibrational spectroscopy of D2O in the same manner as TIP4P/2005 for H2O, although one may still conveniently employ an alternative method that regards OH as putative OD to calculate the OD-stretch spectra with similar quality from TIP4P/2005 trajectories. We also demonstrate that the appropriate setting for the spectral simulation of D2O under the time-averaging approximation reflects the slower dynamics (i.e., slower motion of translation and rotation due to the heavier mass and stronger hydrogen bond) of D2O than H2O. Moreover, we show from the theoretical calculations that the established interpretation of the OH-stretch spectra of H2O is finely applicable to the OD-stretch of D2O.

Transferability of vibrational spectroscopic map from TIP4P to TIP4P-like water models

We computed the IR, Raman, and sum frequency generation spectra of water in the OH-stretch region by employing the quantum/classical mixed approach that consists of a vibrational spectroscopic map and molecular dynamics (MD) simulation. We carried out the MD simulation with the TIP4P, TIP4P/2005, and TIP4P/Ice models and applied the map designed for TIP4P by Skinner et al. to each MD trajectory. Although the map is not tuned for TIP4P-like models, TIP4P/2005 and TIP4P/Ice provide the best reproduction of the experimental vibrational spectra of liquid water and crystalline ice, respectively. This result demonstrates the transferability of the map from TIP4P to TIP4P/2005 and TIP4P/Ice, meaning that one can choose an appropriate TIP4P-like model to calculate the vibrational spectra of an aqueous system without rebuilding the map.

Progress in phase-sensitive sum frequency generation spectroscopy

Sum frequency generation (SFG) spectroscopy is a unique and powerful tool for investigating surfaces and interfaces at the molecular level. Phase-sensitive SFG (PS-SFG) is an upgraded technique that can overcome the inherent drawbacks of conventional SFG. Here we review several methods of PS-SFG developed and reported in 1990-2020. We introduce how and by which group each PS-SFG method was designed and built in terms of interferometer implementation for optical heterodyne detection, with one exception of a recent numerical method that does not rely on interferometry. We also discuss how PS-SFG solved some typical problems for aqueous interfaces that were once left open by conventional SFG. These problems and their solutions are good examples to demonstrate why PS-SFG is essential. In addition, we briefly note a few terminology issues related with PS-SFG to avoid confusion.

Structure of the electrical double layer at the ice-water interface

The surface of ice in contact with water contains sites that undergo deprotonation and protonation and can act as adsorption sites for aqueous ions. Therefore, an electrical double layer should form at this interface and existing models for describing the electrical double layer at metal oxide-water interfaces should be able to be modified to describe the surface charge, surface potential, and ionic occupancy at the ice-water interface. I used a surface complexation model along with literature measurements of the zeta potential of ice in brines of various strength and pH to constrain equilibrium constants. I then made predictions of ion site occupancy, surface charge density, and partitioning of counterions between the Stern and diffuse layers. The equilibrium constant for cation adsorption is more than 5 orders of magnitude larger than the other constants, indicating that this reaction dominates even at low salinity. Deprotonated OH sites are predicted to be slightly more abundant than dangling O sites, consistent with previous work. Surface charge densities are on the order of ±0.001 C/m² and are always negative at the moderate pH values of interest to atmospheric and geophysical applications (6-9). In this pH range, over 99% of the counterions are contained in the Stern layer. This suggests that diffuse layer polarization will not occur because the ionic concentrations in the diffuse layer are nearly identical to those in the bulk electrolyte and that electrical conduction and polarization in the Stern layer will be negligible due to reduced ion mobility.

A hydrogen-bonding structure in self-formed nanodroplets of water adsorbed on amorphous silica revealed via surface-selective vibrational spectroscopy

Water adsorption onto a material surface is known to change macroscopic surface properties such as wettability and friction coefficient.