Sum-frequency generation: Polarization surface spectroscopy analysis of the vibrational surface modes on the basal face of ice Ih

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.

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

Phase-resolved sum-frequency generation measurements combined with molecular dynamics simulations are employed to study the effect of temperature on the molecular arrangement of water on the basal face of ice. The topmost monolayer, interrogated through its nonhydrogen-bonded, free O-H stretch peak, exhibits a maximum in surface H-bond density around 200 K. This maximum results from two competing effects: above 200 K, thermal fluctuations cause the breaking of H bonds; below 200 K, the formation of bulklike crystalline interfacial structures leads to H-bond breaking. Knowledge of the surface structure of ice is critical for understanding reactions occurring on ice surfaces and ice nucleation.

Observation and Identification of a New OH Stretch Vibrational Band at the Surface of Ice

We study the signatures of the OH stretch vibrations at the basal surface of ice using heterodyne-detected sum-frequency generation and molecular dynamics simulations. At 150 K, we observe seven distinct modes in the sum-frequency response, five of which have an analogue in the bulk, and two pure surface-specific modes at higher frequencies (∼3530 and ∼3700 cm^-1). The band at ∼3530 cm^-1 has not been reported previously. Using molecular dynamics simulations, we find that the ∼3530 cm^-1 band contains contributions from OH stretch vibrations of both fully coordinated interfacial water molecules and water molecules with two donor and one acceptor hydrogen bond.

Ice Surfaces

Ice is a fundamental solid with important environmental, biological, geological, and extraterrestrial impact. The stable form of ice at atmospheric pressure is hexagonal ice, Ih. Despite its prevalence, Ih remains an enigmatic solid, in part due to challenges in preparing samples for fundamental studies. Surfaces of ice present even greater challenges. Recently developed methods for preparation of large single-crystal samples make it possible to reproducibly prepare any chosen face to address numerous fundamental questions. This review describes preparation methods along with results that firmly establish the connection between the macroscopic structure (observed in snowflakes, microcrystallites, or etch pits) and the molecular-level configuration (detected with X-ray or electron scattering techniques). Selected results of probing interactions at the ice surface, including growth from the melt, surface vibrations, and characterization of the quasi-liquid layer, are discussed.

Theoretical and experimental examination of SFG polarization analysis at acetonitrile–water solution surfaces

Polarization analysis of SFG spectroscopy is thoroughly examined in collaboration of SFG measurements and MD simulations.

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

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

Insights into hydrogen bonding via ice interfaces and isolated water

Water in a confined environment has a combination of fewer available configurations and restricted mobility. Both affect the spectroscopic signature. In this work, the spectroscopic signature of water in confined environments is discussed in the context of competing models for condensed water: (1) as a system of intramolecular coupled molecules or (2) as a network with intermolecular dipole-dipole coupled O-H stretches. Two distinct environments are used: the confined asymmetric environment at the ice surface and the near-isolated environment of water in an infrared transparent matrix. Both the spectroscopy and the environment are described followed by a perspective discussion of implications for the two competing models. Despite being a small molecule, water is relatively complex; perhaps not surprisingly the results support a model that blends inter- and intramolecular coupling. The frequency, and therefore the hydrogen-bond strength, appears to be a function of donor-acceptor interaction and of longer-range dipole-dipole alignment in the hydrogen-bonded network. The O-H dipole direction depends on the local environment and reflects intramolecular O-H stretch coupling.

Best face forward: crystal-face competition at the ice-water interface

The ice-water interface plays an important role in determining the outcome of both biological and environmental processes. Under ambient pressure, the most stable form of ice is hexagonal ice (Ih). Experimentally probing the surface free energy between each of the major faces of Ih ice and the liquid is both experimentally and theoretically challenging. The basis for the challenge is the near-equality of the surface free energy for the major faces along with the tendency of water to supercool. As a result, morphology from crystallization initiated below 0 °C is kinetically controlled. The reported work circumvents supercooling consequences by providing a polycrystalline seed, followed by isothermal, equilibrium growth. Natural selection among seeded faces results in a single crystal. A record of the growth front is preserved in the frozen boule. Crystal orientation at the front is revealed by examining the boule cross section with two techniques: (1) viewing between crossed polarizers to locate the optical axis and (2) etching to distinguish the primary-prism face from the secondary-prism face. Results suggest that the most stable ice-water interface at 0 °C is the secondary-prism face, followed by the primary-prism face. The basal face that imparts the characteristic hexagonal shape to snowflakes is a distant third. The results contrast with those from freezing the vapor where the basal and primary-prism faces have comparable free energy followed by the secondary-prism face.

Hydrogen bonding in the prism face of ice I(h) via sum frequency vibrational spectroscopy

The prism face of single crystal ice I(h) has been studied using sum frequency vibrational spectroscopy focusing on identification of resonances in the hydrogen-bonded region. Several modes have been observed at about 3400 cm(-1); each mode is both polarization and orientation dependent. The polarization capabilities of sum frequency generation (SFG) are used in conjunction with the crystal orientation to characterize three vibrational modes. These modes are assigned to three-coordinated water molecules in the top-half bilayer having different bonding and orientation motifs.

Interpretation of IR and Raman line shapes for H2O and D2O ice Ih

Noticeable differences between the vibrational (IR and Raman) spectra of neat H(2)O and D(2)O ice Ih are observed experimentally. Here, we employ our theoretical mixed quantum/classical approach to investigate these differences. We find reasonable agreement between calculated and experimental line shapes at both high and low temperatures. From understanding the structure of ice Ih and its vibrational exciton Hamiltonian, we provide assignments of the IR and Raman spectral features for both H(2)O and D(2)O ice Ih. We find that in H(2)O ice these features are due to strong and weak intermolecular coupling, not to intramolecular coupling. The differences between H(2)O and D(2)O ice spectra are attributed to the significantly stronger intramolecular coupling in D(2)O ice. Our conclusion for both H(2)O and D(2)O ice is that the molecular symmetric and antisymmetric normal modes do not form a useful basis for understanding OH or OD stretch spectroscopy.

Origin of Vibrational Spectroscopic Response at Ice Surface

Since the basal plane surface of ice was first observed by sum frequency generation, an extraordinarily intense band for the hydrogen(H)-bonded OH stretching vibration has been a matter of debate. We elucidate the remarkable spectral feature of the ice surface by quantum mechanics/molecular mechanics calculations. The intense H-bonded band is originated mostly from the "bilayer-stitching" modes of a few surface bilayers, through significant intermolecular charge transfer. The mechanism of enhanced signal is sensitive to the order of the tetrahedral ice structure, as the charge transfer is coupled to the vibrational delocalization.

Hydrogen bonding in the hexagonal ice surface

A recently developed technique in sum frequency generation spectroscopy, polarization angle null (or PAN-SFG), is applied to two orientations of the prism face of hexagonal ice. It is found that the vibrational modes of the surface are similar in different faces. As in the basal face, the prism face of ice contains five dominant resonances: 3096, 3146, 3205, 3253, and 3386 cm(-1). On the basal face, the reddest resonance occurs at 3098 cm(-1); within the bandwidth, the same as the prism face. On both the prism and basal faces, this mode contains a significant quadrupole component and is assigned to the bilayer stitching hydrogen bonds. The bluest of the resonances, 3386 cm(-1), occurs slightly blue-shifted at 3393 cm(-1) in the basal face. The prism face has two orientations: one with the optic or c axis in the input plane (the plane formed by the surface normal and the interrogating beam propagation) and one with the c axis perpendicular to the input plane. The 3386 cm(-1) mode has significant intensity only with the c axis in the input plane. On the basis of these orientation characteristics, the 3386 cm(-1) mode is assigned to double-donor molecules in either the top half bilayer or in the lower half bilayer. On the basis of frequency considerations, it is assigned to double-donor molecules in the top half bilayer. These are water molecules containing a nonbonded lone pair. In addition to identification of the components of the broad hydrogen-bonded region, PAN-SFG measures the tangential vs longitudinal content of the vibrational modes. In accord with previous suggestions, the lower frequency modes are predominantly tangential, whereas the higher frequency modes are mainly longitudinal. On the prism face, the 3386 cm(-1) mode is entirely longitudinal.