Origin of Vibrational Spectroscopic Response at Ice Surface

Effect of counterions on the structure and dynamics of water near a negatively charged surfactant: a theoretical vibrational sum frequency generation study

Charged aqueous interfaces are of paramount importance in electrochemical, biological and environmental sciences. The properties of aqueous interfaces with ionic surfactants can be influenced by the presence of counterions. Earlier experiments involving vibrational sum frequency generation (VSFG) spectroscopy of aqueous interfaces with negatively charged sodium dodecyl sulfate (Na+DS- or SDS) surfactants revealed that the hydrogen bonding strength of the interfacial water molecules follows a certain order when salts of monovalent and divalent cations are added. It is known that cations do not directly participate in hydrogen bonding with water molecules, rather they only influence the hydrogen bonded network through their electrostatic fields. In the current work, we have simulated the aqueous interfacial systems of sodium dodecyl sulfate in the presence of chloride salts of mono and divalent countercations. The electronic polarization effects on the ions are considered at a mean-field level within the electronic continuum correction model. Our calculations of the VSFG spectra show a blue shift in the presence of added countercations whose origin is traced to different relative contributions of water molecules from the solvation shells of the surfactant headgroups and the remaining water molecules in the presence of countercations. Furthermore, the cations shield the electric fields of the surfactant headgroups, which in turn influences the contributions of water molecules to the total VSFG spectrum. This shielding effect becomes more significant when divalent countercations are present. The dynamics of water molecules is found to be slower at the interface in comparison to the bulk. The interfacial depth dependence of various dynamical quantities shows that the interface is structurally and dynamically more heterogeneous at the microscopic level.

Theoretical Study of the Two-Dimensional Vibrational Sum Frequency Generation Spectroscopy of the Air-Water Interface at Varying Temperature and Its Connections to the Interfacial Structure and Dynamics

We performed a theoretical study of the temperature variation of two-dimensional vibrational sum frequency generation (2D-VSFG) spectra of the OH stretch modes at air-water interfaces in the mid-IR region. The calculations are performed at four different temperatures from 250 to 325 K by using a combination of techniques involving response function formalism of nonlinear spectroscopy, electronic structure calculations, and molecular dynamics simulations. Also, the calculations are performed for isotopically dilute solutions so that the intra- and intermolecular coupling between the vibrational modes of interest can be ignored. We have established the connections of temperature variation of various frequency- and time-dependent features of the calculated spectra to the changes in the underlying structure and dynamics of the interfaces. The results reveal that interfacial water is dynamically more heterogeneous than bulk water, with three dominant dynamical processes exhibiting their corresponding time-dependent features in the 2D-VSFG spectrum. These are the spectral diffusion of hydrogen-bonded OH groups at the interface, conversion of an initially hydrogen-bonded OH group to a dangling OH which is a stable state for surface water, unlike the bulk water, and the third one, which involves the conversion of an initially free or dangling OH group to its hydrogen-bonded state at the interface. The temporal appearance of the cross peaks corresponding to interconversion of the hydrogen-bonded state to the dangling state or vice versa of an interfacial OH group is found to take place at a slower rate than the dynamics of spectral diffusion of hydrogen-bonded molecules at the interface, which, in turn, is slower than the corresponding spectral diffusion of bulk water molecules. The temperature variation of these dynamic processes can be linked to the decay of appropriate hydrogen-bond and non-hydrogen-bond time correlation functions of interfacial water molecules for the different air-water systems studied in this work.

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.

Thickness dependent homogeneous crystallization of ultrathin amorphous solid water films

The crystallization mechanism and kinetics of amorphous materials are of paramount importance not only in basic science but also in the application field because they are closely related to their thermal stability. In the case of amorphous nanomaterials, thermal stability distinctively different from that of bulk materials often emerges. Despite intensive studies in the past, a thorough understanding of the stability at the molecular level has not been reached particularly on how crystallization processes depend on size and are influenced by their surface and interface. In this article, we report the film-size-dependent crystallization of thermally relaxed nonporous ASW ultrathin films on a Pt(111) surface as a benchmark system of amorphous molecular films. The crystallization processes at the surface and interior of the ASW ultrathin films are monitored simultaneously with thermal desorption and infrared reflection absorption, respectively, as a function of the film thickness. Here, we demonstrate that the crystallization is initiated solely by "homogeneous nucleation" irrespective of the film thickness while the crystallization rate remarkably depends on the thickness; the rate of 5-layer (∼1.5 nm) ASW films is one order of magnitude higher than that of 20-layer (∼6 nm) films. Moreover, we found a clear correlation between the film-thickness-dependent crystallization kinetics and microscopic structural disorder associated with the broad distribution of hydrogen-bond lengths between water molecules.

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).

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.

Supersolidity of undercoordinated and hydrating water

Electrostatic polarization or molecular undercoordination endows the supersolidity by shortening and stiffening the H–O bond and lengthening and softening the O:H nonbond, deepening the O 1s energy level, and prolonging the photoelectron and phonon lifetime. The supersolid phase is less dense, viscoelastic, mechanically and thermally more stable, which offsets boundaries of structural phases and critical temperatures for phase transition of the coordination-resolved core–shell structured ice such as the ‘no man's land’ supercooling and superheating.

Proton Order toward the Surface of Ice Ih Revealed by Heterodyne-Detected Sum Frequency Generation Spectroscopy

Using heterodyne-detected sum frequency generation (HD-SFG) spectroscopy, we investigated surface proton order at the basal, primary prism, and secondary prism faces of single-crystalline ice I_h at ca. 130 K. The complex phase of the obtained spectra clearly indicates that second-order nonlinear polarization from which the HD-SFG signal arises is generated exclusively at the surfaces. This suggests surface proton ordering along the normal, whereas the bulk remains proton-disordered, as is well known for ice I_h. A strong positive peak observed in the HD-SFG spectra enables us to determine the "direction" of the surface proton order as "H-up", that is, the hydrogen atom of the OH group pointing away from the bulk, irrespective of the ice faces. Reliable HD-SFG measurements carried out in the present study have greatly advanced our understanding of surface structure of ice I_h.

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.

Transient Phase of Ice Observed by Sum Frequency Generation at the Water/Mineral Interface During Freezing

We observed a transient noncentrosymmetric phase of ice at water/mineral interfaces during freezing, which enhanced the intensity of the IR-visible sum frequency generation intensity by up to 20-fold. The lifetime of the transient phase was several minutes. Since the most stable form of ice, hexagonal and cubic ice, are centrosymmetric, our study suggests the transient existence of stacking-disordered ice during the freezing process at water/mineral interfaces. Stacking-disordered ice, which has only been observed in bulk ice at temperatures lower than -20 °C, is a random mixture of layers of hexagonal ice and cubic ice. However, the transient phase at the ice/mineral interface was observed at temperatures as high as -1 °C. It suggests that the mineral surface may play a role in promoting and stabilizing the formation of stacking-disordered ice at the interface.

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.

First-Principles Framework to Compute Sum-Frequency Generation Vibrational Spectra of Semiconductors and Insulators

We present a first-principles framework to compute sum-frequency generation (SFG) vibrational spectra of semiconductors and insulators. The method is based on density functional theory and the use of maximally localized Wannier functions to compute the response to electric fields, and it includes the effect of electric field gradients at surfaces. In addition, it includes quadrupole contributions to SFG spectra, thus enabling the verification of the dipole approximation, whose validity determines the surface specificity of SFG spectroscopy. We compute the SFG spectra of ice I_{h} basal surfaces and identify which spectra components are affected by bulk contributions. Our results are in good agreement with experiments at low temperature.

The surface roughness, but not the water molecular orientation varies with temperature at the water-air interface

We examine the temperature dependence of the interfacial molecular structure at the water-air interface by combining experimental and simulated sum-frequency generation (SFG) spectroscopy. The experimental SFG spectra of the OH-stretching mode show a decrease in the amplitude at ∼3300 cm(-1) with increasing temperature, while the 3700 cm(-1) 'free OH' SFG feature is insensitive to temperature changes. The simulated spectra are in excellent agreement with experiment. A comparison between interfacial SFG spectra and bulk infrared/Raman spectra reveals that the variation of the SFG signal due to the temperature change is not caused by a temperature-dependent OH bond orientation of the interfacial water molecules, but can be fully accounted for by the temperature dependence of the optical response of water. These results indicate that while the thickness of the interfacial region varies with temperature, the molecular organization of interfacial water at the water-air interface is surprisingly insensitive to temperature changes.

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.

A direct evidence of vibrationally delocalized response at ice surface

Surface-specific vibrational spectroscopic responses at isotope diluted ice and amorphous ice are investigated by molecular dynamics (MD) simulations combined with quantum mechanics/molecular mechanics calculations. The intense response specific to the ordinary crystal ice surface is predicted to be significantly suppressed in the isotopically diluted and amorphous ices, demonstrating the vibrational delocalization at the ordinary ice surface. The collective vibration at the ice surface is also analyzed with varying temperature by the MD simulation.

Size, separation, structural order, and mass density of molecules packing in water and ice

The structural symmetry and molecular separation in water and ice remain uncertain. We present herewith a solution to unifying the density, the structure order and symmetry, the size (H-O length dH), and the separation (d(OO) = d(L) + d(H) or the O:H length d(L)) of molecules packing in water and ice in terms of statistic mean. This solution reconciles: i) the d(L) and the d(H) symmetrization of the O:H-O bond in compressed ice, ii) the d(OO) relaxation of cooling water and ice and, iii) the d(OO) expansion of a dimer and between molecules at water surface. With any one of the d(OO), the density ρ(g·cm⁻³), the d(L), and the d(H), as a known input, one can resolve the rest quantities using this solution that is probing conditions or methods independent. We clarified that: i) liquid water prefers statistically the mono-phase of tetrahedrally-coordinated structure with fluctuation, ii) the low-density phase (supersolid phase as it is strongly polarized with even lower density) exists only in regions consisting molecules with fewer than four neighbors and, iii) repulsion between electron pairs on adjacent oxygen atoms dictates the cooperative relaxation of the segmented O:H-O bond, which is responsible for the performance of water and ice.

Density, Elasticity, and Stability Anomalies of Water Molecules with Fewer than Four Neighbors

Goldschmidt-Pauling contraction of the H-O polar-covalent bond elongates and polarizes the other noncovalent part of the hydrogen bond (O:H-O), that is, the O:H van der Waals bond, significantly, through the Coulomb repulsion between the electron pairs of adjacent oxygen (O-O). This process enlarges and stiffens those H2O molecules having fewer than four neighbors such as molecular clusters, hydration shells, and the surface skins of water and ice. The shortening of the H-O bond raises the local density of bonding electrons, which in turn polarizes the lone pairs of electrons on oxygen. The stiffening of the shortened H-O bond increases the magnitude of the O1s binding energy shift, causes the blue shift of the H-O phonon frequencies, and elevates the melting point of molecular clusters and ultrathin films of water, which gives rise to their elastic, hydrophobic, highly-polarized, ice-like, and low-density behavior at room temperature.

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.

Quantum Calculations of Intramolecular IR Spectra of Ice Models Using Ab Initio Potential and Dipole Moment Surfaces

We report the IR spectra of two forms of ice in the monomer bend and OH-stretching regions, using recently developed ab initio potential and dipole moment surfaces for arbitrarily many water monomers. Coupling and anharmonicity of the intramolecular vibrational modes are taken into account using coupled three-mode variational calculations, within the local-monomer model. Spectra for the surface and core regions of these ice models are presented. The calculated spectra for the core region, with no adjustments, are in good agreement with experiment for the intramolecular OH-stretch and bend regions. Our analysis also shows a significant contribution from the overtone of the monomer bend to the OH-stretch region of the spectra.