Impact of intermolecular vibrational coupling effects on the sum-frequency generation spectra of the water/air interface

Modeling Vibrational Sum Frequency Generation Spectra of Interfacial Water on a Gold Surface: The Role of the Fermi Resonance

Studying the hydrogen bonding structure of H2O at the metal-water interface is a highly complex yet fascinating endeavor. The intricate interactions and diverse orientations of water molecules on metal surfaces with varying potentials pose a significant challenge in elucidating the coupling between O-H stretching and H-O-H bending modes. In this study, we employed DFT-MD simulation to explore how the orientation of interfacial water molecules changes with the applied potential on the Au(111) surface. Based on the surface-specific velocity-velocity correlation function (ssVVCF) formula, we calculated vibrational sum frequency generation (VSFG) spectra for the O-H stretches. We found that three assigned peaks (∼3300, ∼3450, and 3650 cm-1) shifted toward lower frequencies as the potential moved toward more negative values. Our results align remarkably well with experimental Raman spectroscopy data. Notably, our VSFG analysis revealed a significant change in the VSFG spectra of the hydrogen-bonded O-H groups (∼3300 cm-1), switching from a negative to a positive sign with decreasing potential. This alteration suggests a substantial change in the orientation of these low-frequency O-H groups owing to their increased interactions with the Au surface. In contrast, the orientations of both the high-frequency O-H groups (∼3450 cm-1) and the dangling O-H groups (∼3650 cm-1) remained unaffected by the applied potentials. Furthermore, our analysis of the decomposed vibrational density of states (VDOS) for the H-O-H bending mode uncovered the coupling between the H-O-H bending and O-H stretching vibrations, known as the Fermi resonance. Our work suggests that the H-O-H bending vibration becomes restricted when water molecules transition from the ″one-H-down″ to the ″two-H-down″ conformation, leading to a redshift in the O-H stretching vibration through the Fermi resonance. By constructing the VSFG and decomposed VDOS spectra, we gained valuable insights into the structural changes that Raman spectra alone cannot fully interpret. Specifically, our analysis revealed the critical role of the Fermi resonance effect in shaping the spectroscopic signature of interfacial water molecules on the Au(111) surface.

Nuclear quantum effects on the vibrational dynamics of the water-air interface

We have applied path-integral molecular dynamics simulations to investigate the impact of nuclear quantum effects on the vibrational dynamics of water molecules at the water-air interface. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time series analysis, while the time scales of vibrational spectral diffusion are determined from frequency-time correlation functions and joint probability distributions. We find that the inclusion of nuclear quantum effects leads not only to a redshift in the vibrational frequency distribution by about 120 cm-1 for both the bulk and interfacial water molecules but also to an acceleration of the vibrational dynamics at the water-air interface by as much as 35%. In addition, a blueshift of about 45 cm-1 is seen in the vibrational frequency distribution of interfacial water molecules compared to that of the bulk. Furthermore, the dynamics of water molecules beyond the topmost molecular layer was found to be rather similar to that of bulk water.

Bulklike Vibrational Coupling of Surface Water Revealed by Sum-Frequency Generation Spectroscopy

Vibrational coupling between interfacial water molecules is important for energy dissipation after on-water chemistry, yet intensely debated. Here, we quantify the interfacial vibrational coupling strength through the linewidth of surface-specific vibrational spectra of the water's O─H (O─D) stretch region for neat H_{2}O/D_{2}O and their isotopic mixtures. The local-field-effect-corrected experimental SFG spectra reveal that the vibrational coupling between hydrogen-bonded interfacial water O─H groups is comparable to that in bulk water, despite the effective density reduction at the interface.

Fully First-Principles Surface Spectroscopy with Machine Learning

Our current understanding of the structure and dynamics of aqueous interfaces at the molecular level has grown substantially due to the continuous development of surface-specific spectroscopies, such as vibrational sum-frequency generation (VSFG). As in other vibrational spectroscopies, we must turn to atomistic simulations to extract all of the information encoded in the VSFG spectra. The high computational cost associated with existing methods means that they have limitations in representing systems with complex electronic structure or in achieving statistical convergence. In this work, we combine high-dimensional neural network interatomic potentials and symmetry-adapted Gaussian process regression to overcome these constraints. We show that it is possible to model VSFG signals with fully ab initio accuracy using machine learning and illustrate the versatility of our approach on the water/air interface. Our strategy allows us to identify the main sources of theoretical inaccuracy and establish a clear pathway toward the modeling of surface-sensitive spectroscopy of complex interfaces.

Characterizing Interfaces by Voronoi Tessellation

The chemistry of interfaces differs markedly from that of the bulk. Calculation of interfacial properties depends strongly on the definition of the interface, which can lead to ambiguous results that vary between studies. There is a need for a method that can explicitly define the interfaces and boundaries in molecular systems. Voronoi tessellation offers an attractive solution to this problem through its ability to determine neighbors among specified groups of atoms. Here we discuss three cases where Voronoi tessellation combined with modeling of vibrational sum frequency generation (SFG) spectroscopy yields relevant insights: the breakdown of the air-water interface into clear and intuitive molecular layers, the study of the hydration shell in biological systems, and the acceleration of difficult spectral calculations where intermolecular vibrational couplings dominate. The utility of Voronoi tessellation has broad applications that extend beyond any single type of spectroscopy or system.

Is Unified Understanding of Vibrational Coupling of Water Possible? Hyper-Raman Measurement and Machine Learning Spectra

The impact of the vibrational coupling of the OH stretch mode on the spectra differs significantly between IR and Raman spectra of water. Unified understanding of the vibrational couplings is not yet achieved. By using a different class of vibrational spectroscopy, hyper-Raman (HR) spectroscopy, together with machine-learning-assisted HR spectra calculation, we examine the impact of the vibrational couplings of water through the comparison of isotopically diluted H₂O and pure H₂O. We found that the isotopic dilution reduces the HR bandwidths, but the impact of the vibrational coupling is smaller than in the IR and parallel-polarized Raman. Machine learning HR spectra indicate that the intermolecular coupling plays a major role in broadening the bandwidth, while the intramolecular coupling is negligibly small, which is consistent with the IR and Raman spectra. Our result clearly demonstrates a limited impact of the intramolecular vibration, independent of the selection rules of vibrational spectroscopies.

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.

On the Fresnel factor correction of sum-frequency generation spectra of interfacial water

Insights into the microscopic structure of aqueous interfaces are essential for understanding the chemical and physical processes on the water surface, including chemical synthesis, atmospheric chemistry, and events in biomolecular systems. These aqueous interfaces have been probed by heterodyne-detected sum-frequency generation (HD-SFG) spectroscopy. To obtain the molecular response from the measured HD-SFG spectra, one needs to correct the measured ssp spectra for local electromagnetic field effects at the interface due to a spatially varying dielectric function. This so-called Fresnel factor correction can change the inferred response substantially, and different ways of performing this correction lead to different conclusions about the interfacial water response. Here, we compare the simulated and experimental spectra at the air/water interface. We use three previously developed models to compare the experiment with theory: an advanced approach taking into account the detailed inhomogeneous interfacial dielectric profile and the Lorentz and slab models to approximate the interfacial dielectric function. Using the advanced model, we obtain an excellent quantitative agreement between theory and experiment, in both spectral shape and amplitude. Remarkably, we find that for the Fresnel factor correction of the ssp spectra, the Lorentz model for the interfacial dielectric function is equally accurate in the hydrogen (H)-bonded region of the response, while the slab model underestimates this response significantly. The Lorentz model, thus, provides a straightforward method to obtain the molecular response from the measured spectra of aqueous interfaces in the H-bonded region.

Simulation of the Chiral Sum Frequency Generation Response of Supramolecular Structures Requires Vibrational Couplings

Chiral vibrational sum frequency generation (SFG) spectroscopy probes the structure of the solvation shell around chiral macromolecules. The dominant theoretical framework for understanding the origin of chiral SFG signals is based on the analysis of molecular symmetry, which assumes no interaction between molecules. However, water contains strong intermolecular interactions that significantly affect its properties. Here, the role of intermolecular vibrational coupling in the chiral SFG response of the O-H stretch of water surrounding an antiparallel β-sheet at the vacuum-water interface is investigated. Both intramolecular and intermolecular couplings between O-H groups are required to simulate the full lineshape of the chiral SFG signal. This dependence is also observed for a chiral water dimer, illustrating that this phenomenon is not specific to larger systems. We also find that a dimer of C3v molecules predicted to be chirally SFG-inactive by the symmetry-based theory can generate a chiral SFG signal when intermolecular couplings are considered, suggesting that even highly symmetric solvent molecules may produce chiral SFG signals when interacting with a chiral solute. The consideration of intermolecular couplings extends the prevailing theory of the chiral SFG response to structures larger than individual molecules and provides guidelines for future modeling.

Energy relaxation dynamics of hydrogen-bonded OH vibration conjugated with free OH bond at an air/water interface

Vibrational energy relaxation dynamics of the excited hydrogen-bonded (H-bonded) OH conjugated with free OH (OD) at an air/water (for both pure water and isotopically diluted water) interface are elucidated via non-equilibrium ab initio molecular dynamics (NE-AIMD) simulations. The calculated results are compared with those of the excited H-bonded OH in bulk liquid water reported previously. In the case of pure water, the relaxation timescale (vibrational lifetime) of the excited H-bonded OH at the interface is T₁ = 0.13 ps, which is slightly larger than that in the bulk (T₁ = 0.11 ps). Conversely, in the case of isotopically diluted water, the relaxation timescale of T₁ = 0.74 ps in the bulk decreases to T₁ = 0.26 ps at the interface, suggesting that the relaxation dynamics of the H-bonded OH are strongly dependent on the surrounding H-bond environments particularly for the isotopically diluted conditions. The relaxation paths and their rates are estimated by introducing certain constraints on the vibrational modes except for the target path in the NE-AIMD simulation to decompose the total energy relaxation rate into contributions to possible relaxation pathways. It is found that the main relaxation pathway in the case of pure water is due to intermolecular OH⋯OH vibrational coupling, which is similar to the relaxation in the bulk. In the case of isotopically diluted water, the main pathway is due to intramolecular stretch and bend couplings, which show more efficient relaxation than in the bulk because of strong H-bonding interactions specific to the air/water interface.

Molecular Fingerprints of Hydrophobicity at Aqueous Interfaces from Theory and Vibrational Spectroscopies

Hydrophobicity/hydrophilicity of aqueous interfaces at the molecular level results from a subtle balance in the water-water and water-surface interactions. This is characterized here via density functional theory-molecular dynamics (DFT-MD) coupled with vibrational sum frequency generation (SFG) and THz-IR absorption spectroscopies. We show that water at the interface with a series of weakly interacting materials is organized into a two-dimensional hydrogen-bonded network (2D-HB-network), which is also found above some macroscopically hydrophilic silica and alumina surfaces. These results are rationalized through a descriptor that measures the number of "vertical" and "horizontal" hydrogen bonds formed by interfacial water, quantifying the competition between water-surface and water-water interactions. The 2D-HB-network is directly revealed by THz-IR absorption spectroscopy, while the competition of water-water and water-surface interactions is quantified from SFG markers. The combination of SFG and THz-IR spectroscopies is thus found to be a compelling tool to characterize the finest details of molecular hydrophobicity at aqueous interfaces.

Hydrogen bond dynamics of interfacial water molecules revealed from two-dimensional vibrational sum-frequency generation spectroscopy

Vibrational sum-frequency generation (vSFG) spectroscopy allows the study of the structure and dynamics of interfacial systems. In the present work, we provide a simple recipe, based on a narrowband IR pump and broadband vSFG probe technique, to computationally obtain the two-dimensional vSFG spectrum of water molecules at the air-water interface. Using this technique, to study the time-dependent spectral evolution of hydrogen-bonded and free water molecules, we demonstrate that at the interface, the vibrational spectral dynamics of the free OH bond is faster than that of the bonded OH mode.

"On-The-Fly" Calculation of the Vibrational Sum-Frequency Generation Spectrum at the Air-Water Interface

In the present work, we provide an electronic structure based method for the "on-the-fly" determination of vibrational sum frequency generation (v-SFG) spectra. The predictive power of this scheme is demonstrated at the air-water interface. While the instantaneous fluctuations in dipole moment are obtained using the maximally localized Wannier functions, the fluctuations in polarizability are approximated to be proportional to the second moment of Wannier functions. The spectrum henceforth obtained captures the signatures of hydrogen bond stretching, bending, as well as low-frequency librational modes.

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.

Time-dependent vibrational sum-frequency generation spectroscopy of the air-water interface

Vibrational sum-frequency generation spectroscopy is a powerful method to study the microscopic structure and dynamics of interfacial systems. Here we demonstrate a simple computational approach to calculate the time-dependent, frequency-resolved vibrational sum-frequency generation spectrum (TD-vSFG) of the air-water interface. Using this approach, we show that at the air-water interface, the transition of water molecules with bonded OH modes to free OH modes occurs at a time scale of $$\sim$$ ~ 3 ps, whereas water molecules with free OH modes rapidly make a transition to a hydrogen-bonded state within $$\sim$$ ~ 2 ps. Furthermore, we also elucidate the origin of the observed differential dynamics based on the time-dependent evolution of water molecules in the different local solvent environments.

Nuclear Quantum Effect on the χ(2) Band Shape of Vibrational Sum Frequency Generation Spectra of Normal and Deuterated Water Surfaces

In the vibrational sum frequency generation (VSFG) spectrum of the air/water interface, there is an open question whether the imaginary spectrum of nonlinear susceptibility Im[χ⁽²⁾] has a positive band in the low-frequency tail of the OH stretching band. We previously elucidated by molecular dynamics (MD) analysis of the VSFG spectrum that the positive band could arise from particularly strong hydrogen-bonded water pairs at the water surface. This mechanism should be emphasized in OD stretching in comparison to OH, because OD forms stronger hydrogen bonds. Therefore, we calculated the Im[χ⁽²⁾] spectra of normal and deuterated water surfaces by MD simulation including the nuclear quantum effect and demonstrated that the low-frequency positive feature could arise in the tail of the OD stretching band. This positive feature is sensitive to the oxygen-oxygen distance of hydrogen-bonding pairs at the surface and hence the temperature and disappears with increasing temperature.