Nuclear quantum effects on the vibrational dynamics of liquid water

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.

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.

Ionic Dynamics and Vibrational Spectral Diffusion of a Protic Alkylammonium Ionic Salt through Intrinsic Cationic N-H Vibrational Probe from FPMD Simulations

We employed density functional theory (DFT)-based molecular dynamics simulations to explore the structure, dynamics, and spectral properties of the protic ionic entity trimethylammonium chloride (TMACl). Structural investigations include calculating the site-site radial distribution functions (RDFs), the distribution of constituent cations and anions in three-dimensional space, and combined distribution functions of the hydrogen-bonded pair RDF versus angle, revealing the structural characteristics of the ionic solvation and the intermolecular interactions within ions. Further, we determined the instantaneous vibrational stretching frequencies of the intrinsic N-H stretch probe modes by applying the time-series wavelet method. The associated ionic dynamics within the protic ionic compound were investigated by calculating the time-evolution of the fluctuating frequencies and the frequency-time correlation functions (FFCFs). The time scale related to the local structural relaxation process and the average hydrogen bond lifetime, ion cage dynamics, and mean squared displacement were investigated. The faster decay component of the FFCFs, depicting the intermolecular motion of intact hydrogen bonds in TMACl, is 0.07 ps for the Perdew-Burke-Ernzerhof (PBE)-based simulation and 0.06 ps for the PBE-D2 representation. The slower time scale of the longer picosecond decay time component of PBE and PBE-D2 representations are 3.13 and 2.87 ps, respectively. These picosecond time scales represent more significant fluctuations of the hydrogen-bonding partners in the ionic entity and hydrogen-bond jump events accompanied by large angular jumps. The longest picosecond time scales represent structural relaxation, including large angular jumps and ion-pair dynamics. Also, ion cage lifetimes correlate with the slowest time scale of the associated dynamics of vibrational spectral diffusion despite the type of DFT functional. This study benchmarks DFT treatments of the exchange-correlation functional with and without the van der Waals (vdW) dispersion correction scheme. The inclusion of vdW interactions to the PBE functional represents a less structured state of the ionic entity and faster dynamics of the molecular motions relative to the one predicted by the PBE system. All the results illustrate the necessity of accurately describing the Coulomb interactions, vdW dispersive interactive forces, and localized hydrogen bonds required to sustain the energetic balance in this ionic salt.

Structures of liquid and aqueous water isotopologues at ambient temperature from ab initio path integral simulations

The heavy hydrogen isotopes D and T are found in trace amounts in water, and when their concentration increases they can play an intricate role in modulating the physical properties of the liquid. We present an analysis of the microscopic structures of ambient light water (H₂O(l)), heavy water (D₂O(l)), T₂O(l), HDO(aq) and HTO(aq) studied by ab initio path integral molecular dynamics (PIMD). Unlike previous ab initio PIMD investigations of H₂O(l) and D₂O(l) [Chen et al., Phys. Rev. Lett., 2003, 91, 215503] [Machida et al., J. Chem. Phys., 2017, 148, 102324] we find that D₂O(l) is more structured than H₂O(l), as is predicted by the experiment. The agreement between the experiment and our simulation for H₂O(l) and D₂O(l) allows us to accurately predict the intra- and intermolecular structures of T₂O(l) HDO(aq) and HTO(aq). T₂O(l) is found to have a similar intermolecular structure to that of D₂O(l), while the intramolecular structure is more compact, giving rise to a smaller dipole moment than those of H₂O(l) and D₂O(l). For the mixed isotope species, HDO(aq) and HTO(aq), we find smaller dipole moments and fewer hydrogen bonds when compared with the pure species H₂O and D₂O. We can attribute this effect to the relative compactness of the mixed isotope species, which results in a lower dipole moment than that of the pure species.

Ab initio study of nuclear quantum effects on sub- and supercritical water

The structures of water in the ambient, subcritical, and supercritical conditions at various densities were studied systematically by ab initio path integral molecular dynamics simulations. It was found that the nuclear quantum effects (NQEs) have a significant impact on the structure of hydrogen bonds in close contact, not only in the ambient condition but also in the sub- and supercritical conditions. The NQEs on the structure beyond the hydrogen bond contact are important in ambient water, but not much for water in the sub- and supercritical conditions. The NQEs are furthermore important for determining the number of hydrogen bonds in the ambient conditions, and this role is, however, diminished in the sub- and supercritical conditions. The NQEs do, nevertheless, show their importance in determining the intramolecular structure of water and the close contact structures of the hydrogen bonds, even at sub- and supercritical conditions. Using the RPBE-D3 functional, the computed radial distribution functions for ambient water are in excellent agreement with experimental data, upgrading our previous results using the BLYP-D2 functional [Machida et al., J. Chem. Phys. 148, 102324 (2018)]. The computed radial distribution functions for water in the sub- and supercritical conditions were carefully compared with experiment. In particular, we found that the first peak in hydrogen pair distribution functions matches only when the NQEs are taken into account.

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.

Tumbling with a limp: local asymmetry in water's hydrogen bond network and its consequences

Ab initio molecular dynamics simulations of ambient liquid water and energy decomposition analysis have recently shown that water molecules exhibit significant asymmetry between the strengths of the two donor and/or the two acceptor interactions.

On the microscopic origin of the cryoprotective effect in lysine solutions

Lysine cryoprotective properties are due to the tight bonding of the first hydration Shell to the amino acid. However this effect is only possible for concentration up to 5.4 water molecules per lysine.

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.

Ab initio molecular dynamics simulation of vibrational energy redistribution of selective excitation of C-H stretching vibrations for solid nitromethane

Vibrational energy redistribution (VER) of energetic materials plays an important role in transferring the injected energy to the hot spots, but it is extremely challenging to understand the mechanism of VER from experimental or theoretical studies. Here, we combined nonequilibrium molecular dynamics with density functional theory to study the processes of VER for solid nitromethane after the selective excitation of the C-H stretching vibration. The VER processes are traced by monitoring the normal-mode kinetic energies of both excited and unexcited vibrations. To explore the underlying VER mechanism, we also analyzed the spectral energy density for the normal mode, obtained from the squared modulus of the short-time Fourier transition of their normal mode momentum. The results showed that the simulated VER progress was reproduced well compared with the previous 3D IR-Raman experiments of liquid nitromethane. Interestingly, the symmetric dependence of the coupling mechanism between the normal modes has been found.

Vibrational echo spectroscopy of aqueous sodium bromide solutions from first principles simulations

A theoretical study of the time-dependent vibrational echo spectroscopy of sodium bromide solutions in deuterated water at two different concentrations of 0.5 and 5.0 M and at temperatures of 300 and 350 K is presented using the method of ab initio molecular dynamics simulations. The instantaneous fluctuations in frequencies of local OD stretch modes are calculated using time-series analysis of the simulated trajectories. The third-order polarization and intensities of three pulse photon-echo are calculated from ab initio simulations. The timescales of vibrational spectral diffusion are determined from the frequency time correlation functions (FTCF) and short-time slope of three pulse photon echo (S3PE) calculated within the second-order cumulant and Condon approximations. It is found that under ambient conditions, the rate of vibrational spectral diffusion becomes slower with increase in ionic concentration. Decay of S3PE calculated for different systems give timescales, which are in close agreement with those of FTCF and also with the results of experimental time-dependent vibrational spectroscopic experiments. © 2019 Wiley Periodicals, Inc.

Temperature dependence of the ultrafast vibrational echo spectroscopy of OD modes in liquid water from first principles simulations

The temperature dependence of the vibrational spectral diffusion of OD modes in liquid water is investigated through calculations of vibrational echo spectral observables from first principles molecular dynamics.

On the Hydrogen Bond Strength and Vibrational Spectroscopy of Liquid Water

In the present work, we introduce two new metrics i.e. hydrogen-bond strength and charge-transfer between the donor/acceptor water molecules as a measure of hydrogen-bond rearrangement dynamics. Further, we also employ a simple model based on energy flux through the donor-acceptor water pairs to quantify the extent of the local hydrogen-bond network reorganization. Most importantly, we report a linear relationship between the OH stretch frequency and the charge and energy transfer through donor-acceptor water pairs. We demonstrate that the vibrational frequency fluctuations, which are used to determine third-order non-linear spectroscopic observables like the short-time slope of three pulse photon echo, can be used as an analog of the fluctuations in the hydrogen-bond strength and charge-transfer. The timescales obtained from our hydrogen-bond strength correlation and charge-transfer correlation decay are in excellent agreement with the computed frequency-time correlation function, as well as with recent vibrational echo experiments.

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.