Hydrogen bonding and vibrational energy relaxation of interfacial water: A full DFT molecular dynamics simulation

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.

Unified picture of vibrational relaxation of OH stretch at the air/water interface

The elucidation of the energy dissipation process is crucial for understanding various phenomena occurring in nature. Yet, the vibrational relaxation and its timescale at the water interface, where the hydrogen-bonding network is truncated, are not well understood and are still under debate. In the present study, we focus on the OH stretch of interfacial water at the air/water interface and investigate its vibrational relaxation by femtosecond time-resolved, heterodyne-detected vibrational sum-frequency generation (TR-HD-VSFG) spectroscopy. The temporal change of the vibrationally excited hydrogen-bonded (HB) OH stretch band (ν=1→2 transition) is measured, enabling us to determine reliable vibrational relaxation (T1) time. The T1 times obtained with direct excitations of HB OH stretch are 0.2-0.4 ps, which are similar to the T1 time in bulk water and do not noticeably change with the excitation frequency. It suggests that vibrational relaxation of the interfacial HB OH proceeds predominantly with the intramolecular relaxation mechanism as in the case of bulk water. The delayed rise and following decay of the excited-state HB OH band are observed with excitation of free OH stretch, indicating conversion from excited free OH to excited HB OH (~0.9 ps) followed by relaxation to low-frequency vibrations (~0.3 ps). This study provides a complete set of the T1 time of the interfacial OH stretch and presents a unified picture of its vibrational relaxation at the air/water interface.

Ultrafast vibrational dynamics of the free OD at the air/water interface: Negligible isotopic dilution effect but large isotope substitution effect

Vibrational relaxation dynamics of the OH stretch of water at the air/water interface has been a subject of intensive research, facilitated by recent developments in ultrafast interface-selective nonlinear spectroscopy. However, a reliable determination of the vibrational relaxation dynamics in the OD stretch region at the air/D₂O interface has not been yet achieved. Here, we report a study of the vibrational relaxation of the free OD carried out by time-resolved heterodyne-detected vibrational sum frequency generation spectroscopy. The results obtained with the aid of singular value decomposition analysis indicate that the vibrational relaxation (T₁) time of the free OD at the air/D₂O interface and air/isotopically diluted water (HOD-H₂O) interfaces show no detectable isotopic dilution effect within the experimental error, as in the case of the free OH in the OH stretch region. Thus, it is concluded that the relaxation of the excited free OH/OD predominantly proceeds with their reorientation, negating a major contribution of the intramolecular energy transfer. It is also shown that the T₁ time of the free OD is substantially longer than that of the free OH, further supporting the reorientation relaxation mechanism. The large difference in the T₁ time between the free OD and the free OH (factor of ∼2) may indicate the nuclear quantum effect on the diffusive reorientation of the free OD/OH because this difference is significantly larger than the value expected for a classical rotational motion.

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.

Ab initio molecular dynamics study on energy relaxation path of hydrogen-bonded OH vibration in bulk water

The vibrational energy relaxation paths of hydrogen-bonded (H-bonded) OH excited in pure water and in isotopically diluted (deuterated) water are elucidated via non-equilibrium ab initio molecular dynamics (NE-AIMD) simulations. The present study extends the previous NE-AIMD simulation for the energy relaxation of an excited free OH vibration at an air/water interface [T. Ishiyama, J. Chem. Phys. 154, 104708 (2021)] to the energy relaxation of an excited H-bonded OH vibration in bulk water. The present simulation shows that the excited OH vibration in pure water dissipates its energy on a timescale of 0.1 ps, whereas that in deuterated water relaxes on a timescale of 0.7 ps, consistent with the experimental observations. To decompose these relaxation energies into the components due to intramolecular and intermolecular couplings, constraints are introduced on the vibrational modes except for the target path in the NE-AIMD simulation. In the case of pure water, 80% of the total relaxation is attributed to the pathway due to the resonant intermolecular OH⋯OH stretch coupling, and the remaining 17% and 3% are attributed to intramolecular couplings with the bend overtone and with the conjugate OH stretch, respectively. This result strongly supports a significant role for the Förster transfer mechanism of pure water due to the intermolecular dipole-dipole interactions. In the case of deuterated water, on the other hand, 36% of the total relaxation is due to the intermolecular stretch coupling, and all the remaining 64% arises from coupling with the intramolecular bend overtone.

Vibrational energy relaxation of interfacial OH on a water-covered α-Al2O3(0001) surface: a non-equilibrium ab initio molecular dynamics study

Vibrational relaxation of adsorbates is a sensitive tool to probe energy transfer at gas/solid and liquid/solid interfaces. The most direct way to study relaxation dynamics uses time-resolved spectroscopy. Here we report on a non-equilibrium ab initio molecular dynamics (NE-AIMD) methodology to model vibrational relaxation of OH vibrations on a hydroxylated, water-covered α-Al₂O₃(0001) surface. In our NE-AIMD approach, after exciting selected O-H bonds their coupling to surface phonons and to the water adlayer is analyzed in detail, by following both the energy flow in time, as well as the time-evolution of Vibrational Density of States (VDOS) curves. The latter are obtained from Time-dependent Correlation Functions (TCFs) and serve as prototypical, generic representatives of time-resolved vibrational spectra. As most important results, (i) we find a few-picosecond lifetime of the excited modes and (ii) identify both hydrogen-bonded aluminols and water molecules in the adsorbed water layer as main dissipative channels, while the direct coupling to Al₂O₃ surface phonons is of minor importance on the timescales of interest. Our NE-AIMD/TCF methodology is powerful for complex adsorbate systems, in principle even reacting ones, and opens a way towards time-resolved vibrational spectroscopy.

Energy relaxation path of excited free OH vibration at an air/water interface revealed by nonequilibrium ab initio molecular dynamics simulation

Nonequilibrium ab initio molecular dynamics (NE-AIMD) simulations are conducted at an air/water interface to elucidate the vibrational energy relaxation path of excited non-hydrogen-bonded (free) OH. A recent time-resolved vibrational sum frequency generation (TR-VSFG) spectroscopy experiment revealed that the relaxation time scales of free OH at the surface of pure water and isotopically diluted water are very similar to each other. In the present study, the dynamics of free OH excited at the surface of pure water and deuterated water are examined with an NE-AIMD simulation, which reproduces the experimentally observed features. The relaxation paths are examined by introducing constraints for the bonds and angles of water molecules relevant to specific vibrational modes in NE-AIMD simulations. In the case of free OH relaxation at the pure water surface, stretching vibrational coupling with the conjugate bond makes a significant contribution to the relaxation path. In the case of the isotopically diluted water surface, the bend (HOD)-stretching (OD) combination band couples with the free OH vibration, generating a relaxation rate similar to that in the pure water case. It is also found that the reorientation of the free OH bond contributes substantially to the relaxation of the free OH vibrational frequency component measured by TR-VSFG spectroscopy.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

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.

Intramolecular vibrational energy redistribution and the quantum ergodicity transition: a phase space perspective

The onset of facile intramolecular vibrational energy flow can be related to features in the connected network of anharmonic resonances in the classical phase space.

Vibration Spectral Dynamics of Weakly Coordinating Water Molecules near an Anion: FPMD Simulations of an Aqueous Solution of Tetrafluoroborate

The extent to which the ions affect the nearby water molecules will decide the structure-making or breaking nature of those ions in aqueous solutions. The effects of a weakly coordinating anion on the structure, dynamics, and vibrational properties of water molecules are not so significant as compared to an anion capable of making strong ion-water hydrogen bonds. The present work deals with the first-principles molecular dynamics study of an aqueous solution of such a weakly coordinating anion, tetrafluoroborate (BF₄^-), using dispersion-corrected DFT-based first-principles molecular dynamics (FPMD) simulations. Various structural, dynamical, and spectral properties, such as radial distribution functions (RDFs), rotational dynamics, vibrational density of states (VDOS), hydrogen bond as well as dangling OH autocorrelation functions, and residence dynamics, were calculated to investigate the effects of the anion on nearby water molecules. The process of spectral diffusion was assessed through a time series wavelet transformation of trajectories obtained from FPMD simulations. The first ion-water solvation shell extends up to 5.5 Å, containing around 20 water molecules. The lifetime of the ion-water hydrogen bond is found to be 1.19 ps, whereas the water-water hydrogen bond lifetime is found to be 1.13 ps. Inside the solvation shell, the persistence time of dangling OH chromophores and the average frequency of OH modes inside the solvation shell are found to be more compared to bulk. Three time scales are found for solvation shell OH modes from the frequency-frequency correlation function. A very short time scale is found for the intact ion-water interaction; the short time scale is for the ion-water hydrogen bond, and the long time scale is for escape dynamics of water molecules from the ion solvation shell. From the mean squared displacement, it is found that solvation water molecules diffuse slower than the bulk. However, solvation shell water molecules show faster relaxation from the analysis of rotational anisotropy. Within the longer time scale of spectral diffusion, this process (which is related to various dynamics of the molecules) is not yet complete, as compared to fast anisotropic decay. This fact is similar to the experimental finding of spectral diffusion and anisotropy time scales in the aqueous solution of borohydride anion. The calculated results are also compared with available experimental data wherever possible.

Ab initio Modeling of the Vibrational Sum-Frequency Generation Spectrum of Interfacial Water

Understanding the structural and dynamical features of interfacial water is of greatest interest in physics, chemistry, biology, and materials science. Vibrational sum-frequency generation (SFG) spectroscopy, which is sensitive to the molecular orientation and dynamics on the surfaces or at the interfaces, allows one to study a wide variety of interfacial systems. The structural and dynamical features of interfacial water at the air/water interface have been extensively investigated by SFG spectroscopy. However, the interpretations of the spectroscopic features have been under intense debate. Here, we report a simulated SFG spectrum of the air/water interface based on ab initio molecular dynamics simulations, which covers the OH stretching, bending, and libration modes of interfacial water. Quantitative agreement between our present simulations and the most recent experimental studies ensures that ab initio simulations predict unbiased structural features and electrical properties of interfacial systems. By utilizing the kinetic energy spectral density (KESD) analysis to decompose the simulated spectra, the spectroscopic features can then be assigned to specific hydrogen-bonding configurations of interfacial water molecules.

A Microscopic Interpretation of Pump-Probe Vibrational Spectroscopy Using Ab Initio Molecular Dynamics

What happens when extra vibrational energy is added to water? Using nonequilibrium molecular dynamics simulations, also including the full electronic structure, and novel descriptors, based on projected vibrational density of states, we are able to follow the flow of excess vibrational energy from the excited stretching and bending modes. We find that the energy relaxation, mostly mediated by a stretching-stretching coupling in the first solvation shell, is highly heterogeneous and strongly depends on the local environment, where a strong hydrogen bond network can transport energy with a time scale of 200 fs, whereas a weaker network can slow down the transport by a factor 2-3.

Structure from Dynamics: Vibrational Dynamics of Interfacial Water as a Probe of Aqueous Heterogeneity

The structural heterogeneity of water at various interfaces can be revealed by time-resolved sum-frequency generation spectroscopy. The vibrational dynamics of the O–H stretch vibration of interfacial water can reflect structural variations. Specifically, the vibrational lifetime is typically found to increase with increasing frequency of the O–H stretch vibration, which can report on the hydrogen-bonding heterogeneity of water. We compare and contrast vibrational dynamics of water in contact with various surfaces, including vapor, biomolecules, and solid interfaces. The results reveal that variations in the vibrational lifetime with vibrational frequency are very typical, and can frequently be accounted for by the bulk-like heterogeneous response of interfacial water. Specific interfaces exist, however, for which the behavior is less straightforward. These insights into the heterogeneity of interfacial water thus obtained contribute to a better understanding of complex phenomena taking place at aqueous interfaces, such as photocatalytic reactions and protein folding.