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

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.

Vibrational dynamics and spectroscopy of water at porous g-C3N4 and C2N surfaces

Porous graphitic materials containing nitrogen are promising catalysts for photo(electro)chemical reactions, notably water splitting, but can also serve as "molecular sieves". Nitrogen increases the hydrophilicity of the graphite parent material, among other effects. A deeper understanding of how water interacts with C- and N-containing layered materials, if and which differences exist between materials with different N content and pore size, and what the role of water dynamics is - a prerequsite for catalysis and sieving - is largely absent, however. Vibrational spectroscopy can answer some of these questions. In this work, the vibrational dynamics and spectroscopy of deuterated water molecules (D2O) mimicking dense water layers at room temperature on the surfaces of two different C/N-based materials with different N content and pore size, namely graphitic C3N4 (g-C3N4) and C2N, are studied using ab initio molecular dynamics (AIMD). In particular, time-dependent vibrational sum frequency generation (TD-vSFG) spectra of the OD modes and also time-averaged vSFG spectra and OD frequency distributions are computed. This allows us to distinguish "free" (dangling) OD bonds from OD bonds that are bound in a H-bonded water network or at the surface - with subtle differences between the two surfaces and also to a pure water/air interface. It is found that the temporal decay of OD modes is very similar on both surfaces with a correlation time near 4 ps. In contrast, TD-vSFG spectra reveal that the interconversion time from "bonded" to "free" OD bonds is about 8 ps for water on C2N and thus twice as long as for g-C3N4, demonstrating a propensity of the former material to stabilize bonded OD bonds.

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.

Full Assignment of Ab-Initio Raman Spectra at Finite Temperatures Using Wannier Polarizabilities: Application to Cyclohexane Molecule in Gas Phase

We demonstrate how to fully ascribe Raman peaks simulated using ab initio molecular dynamics to specific vibrations in the structure at finite temperatures by means of Wannier functions. Here, we adopt our newly introduced method for the simulation of the Raman spectra in which the total polarizability of the system is expressed as a sum over Wannier polarizabilities. The assignment is then based on the calculation of partial Raman activities arising from self- and/or cross-correlations between different types of Wannier functions in the system. Different types of Wannier functions can be distinguished based on their spatial spread. To demonstrate the predictive power of this approach, we applied it to the case of a cyclohexane molecule in the gas phase and were able to fully assign the simulated Raman peaks.

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.

Assessing the Molecular Specificity and Orientation Sensitivity of Infrared, Raman, and Vibrational Sum-Frequency Spectra

Linear programming was used to assess the ability of polarized infrared absorption, Raman scattering, and visible–infrared sum-frequency generation to correctly identify the composition of a mixture of molecules adsorbed onto a surface in four scenarios. The first two scenarios consisted of a distribution of species where the polarity of the orientation distribution is known, both with and without consideration of an arbitrary scaling factor between candidate spectra and the observed spectra of the mixture. The final two scenarios have repeated the tests, but assuming that the polarity of the orientation is unknown, so the symmetry-breaking attributes of the second-order nonlinear technique are required. The results indicate that polarized Raman spectra are more sensitive to orientation and molecular identity than the other techniques. However, further analysis reveals that this sensitivity is not due to the high-order angle dependence of Raman, but is instead attributed to the number of unique projections that can be measured in a polarized Raman experiment.