The dynamics of water are examined using ultrafast IR stimulated vibrational echo correlation spectroscopy. The OD hydroxyl stretch of HOD in H2O is probed with 45-fs pulses that have sufficient bandwidth (>400 cm-1) to span the entire broad spectrum. High-quality 2D correlation spectra are obtained having the correct phase relations across the broad hydroxyl band. The correlation spectra are found to evolve on multiple time scales. The time evolution of the vibrational echo correlation spectrum reflects the structural evolution of the hydrogen bond networks. The extended vibrational lifetime of the OD hydroxyl stretch of HOD in H2O facilitates the measurement of hydrogen bond dynamics for longer times than possible in previous studies of the OH stretch. Molecular dynamics simulations/electronic structure calculations are used to obtain the time correlation functions (TCF) for two water models, TIP4P and SPC/E. The TCFs are inputs to full time-dependent diagrammatic perturbation theory calculations, which yield theoretical correlation spectra. Quantitative comparison with the data demonstrates that the two water models somewhat overemphasize the fast fluctuations in water and do not contain a slow enough component to account for the slowest fluctuations. Fits to the data using a phenomenological triexponential TCF yield a slowest component of ∼2 ps, and TIP4P and SPC/E have slowest components of <1 ps. The TCF obtained from the water models and the triexponential TCF reproduce the linear absorption line shape equally well, but all miss to some extent the asymmetric “wing” on the low-energy side of the line. Therefore, the time dependence of the vibrational echo correlation spectra provides a good test for the TCF, but the absorption spectrum does not.

No abstract

IR and Raman (parallel- and perpendicular-polarized) spectra in the OH stretch region for liquid water were measured some years ago, but their interpretation is still controversial. In part, this is because theoretical calculation of such spectra for a neat liquid presents a formidable challenge due to the coupling between vibrational chromophores and the effects of motional narrowing. Recently we proposed an electronic structure/molecular dynamics method for calculating spectra of dilute HOD in liquid D(2)O, which relied on ab initio calculations on clusters to provide a map from nuclear coordinates of the molecules in the liquid to OH stretch frequencies, transition dipoles, and polarizabilities. Here we extend this approach to the calculation of couplings between chromophores. From the trajectories of the fluctuating local-mode frequencies, transition moments, and couplings, we use our recently developed time-averaging approximation to calculate the line shapes. Our results are in good agreement with experiment for the IR and Raman line shapes, and capture the significant differences among them. Our analysis shows that while the coupling between chromophores is relatively modest, it nevertheless produces delocalization of the vibrational eigenstates over up to 12 chromophores, which has a profound effect on the spectroscopy. In particular, our results demonstrate that the peak in the parallel-polarized Raman spectrum at about 3250 wavenumbers is collective in nature.