Noticeable differences between the vibrational (IR and Raman) spectra of neat H(2)O and D(2)O ice Ih are observed experimentally. Here, we employ our theoretical mixed quantum/classical approach to investigate these differences. We find reasonable agreement between calculated and experimental line shapes at both high and low temperatures. From understanding the structure of ice Ih and its vibrational exciton Hamiltonian, we provide assignments of the IR and Raman spectral features for both H(2)O and D(2)O ice Ih. We find that in H(2)O ice these features are due to strong and weak intermolecular coupling, not to intramolecular coupling. The differences between H(2)O and D(2)O ice spectra are attributed to the significantly stronger intramolecular coupling in D(2)O ice. Our conclusion for both H(2)O and D(2)O ice is that the molecular symmetric and antisymmetric normal modes do not form a useful basis for understanding OH or OD stretch spectroscopy.

Conspectus Recent experimental advances in vibrational spectroscopy, such as ultrafast pulses and heterodyne detection, have made it possible to probe the structure and dynamics of bulk and interfacial water in unprecedented detail. We consider three aqueous interfaces: the water liquid/vapor interface, the interface between water and the surfactant headgroups of reverse micelles, and the interface between water and the lipid headgroups of aligned multi-bilayers. In the first case, sum-frequency spectroscopy is used to probe the interface, while in the second and third cases, the confined water pools are sufficiently small that techniques of bulk spectroscopy such as FTIR, pump-probe, 2DIR, etc. can be used to probe the interfacial water. In this review, we discuss our attempts to model these three systems and interpret the existing experiments. In particular, for the water liquid/vapor interface we find that three-body interactions are essential for reproducing the experimental sum-frequency spectrum, and presumably for the structure of the interface as well. The observed spectrum is interpreted as arising from overlapping and cancelling positive and negative contributions from molecules in different hydrogen-bonding environments. For the reverse micelles, our theoretical models confirm that the experimentally observed blue shift of the water OD stretch (for dilute HOD in H2O) arises from weaker hydrogen bonding to sulfonate oxygens. We interpret the observed slow-down in water rotational dynamics as arising from curvature-induced frustration. For the water confined between lipid bilayers, our theoretical models confirm that the experimentally observed red shift of the water OD stretch arises from stronger hydrogen bonding to phosphate oxygens. We develop a model for heterogeneous vibrational lifetime distributions, and implement the model to calculate isotropic and anisotropic pump-probe decays, and compare with experiment.

Infrared spectroscopy of the water OH stretch provides a sensitive probe of the local hydrogen-bonding structure and dynamics of water molecules. Previously, we have utilized a mixed quantum/classical model to calculate vibrational spectroscopic observables for bulk water, ice, the liquid/vapor interface, and small water clusters, as well as water interacting with ions and biological molecules. These studies rely on spectroscopic maps that relate the OH stretching frequency and transition dipole to the local environment around a water molecule. Our spectroscopic maps were parametrized based on water clusters taken from bulk water simulations; in this article, we test the robustness of these maps for water in nonbulk-liquid environments. We find that the frequency, transition dipole, and coupling maps work as well for the water surface, ice Ih, and the water hexamer as they do for liquid water. This suggests that these maps may be generally applied to study the vibrational spectroscopy of water in diverse, potentially heterogeneous environments.

The vibrational spectroscopy of hydration water in dilauroylphosphatidylcholine lipid multi-bilayers is investigated using molecular dynamics simulations and a mixed quantum/classical model for the OD stretch spectroscopy of dilute HDO in H(2)O. FTIR absorption spectra, and isotropic and anisotropic pump-probe decay curves have been measured experimentally as a function of the hydration level of the lipid multi-bilayer, and our goal is to make connection with these experiments. To this end, we use third-order response functions, which allow us to include non-Gaussian frequency fluctuations, non-Condon effects, molecular rotations, and a fluctuating vibrational lifetime, all of which we believe are important for this system. We calculate the response functions using existing transition frequency and dipole maps. From the experiments it appears that there are two distinct vibrational lifetimes corresponding to HDO molecules in different molecular environments. In order to obtain these lifetimes, we consider a simple two-population model for hydration water hydrogen bonds. Assuming a different lifetime for each population, we then calculate the isotropic pump-probe decay, fitting to experiment to obtain the two lifetimes for each hydration level. With these lifetimes in hand, we then calculate FTIR spectra and pump-probe anisotropy decay as a function of hydration. This approach, therefore, permits a consistent calculation of all observables within a unified computational scheme. Our theoretical results are all in qualitative agreement with experiment. The vibrational lifetime of lipid-associated OD groups is found to be systematically shorter than that of the water-associated population, and the lifetimes of each population increase with decreasing hydration, in agreement with previous analysis. Our theoretical FTIR absorption spectra successfully reproduce the experimentally observed red-shift with decreasing lipid hydration, and we confirm a previous interpretation that this shift results from the hydrogen bonding of water to the lipid phosphate group. From the pump-probe anisotropy decay, we confirm that the reorientational motions of water molecules slow significantly as hydration decreases, with water bound in the lipid carbonyl region undergoing the slowest rotations.

Using our newly developed explicit three-body (E3B) water model, we simulate the surface of liquid water. We find that the timescale for hydrogen-bond switching dynamics at the surface is about three times slower than that in the bulk. In contrast, with this model rotational dynamics are slightly faster at the surface than in the bulk. We consider vibrational two-dimensional (2D) sum-frequency generation (2DSFG) spectroscopy as a technique for observing hydrogen-bond rearrangement dynamics at the water surface. We calculate the nonlinear susceptibility for this spectroscopy for two different polarization conditions, and in each case we see the appearance of cross-peaks on the timescale of a few picoseconds, signaling hydrogen-bond rearrangement on this timescale. We thus conclude that this 2D spectroscopy will be an excellent experimental technique for observing slow hydrogen-bond switching dynamics at the water surface. Interfaces play important roles in many disciplines of science. The water liquid/vapor interface, for example, is of great interest in chemistry, biology, and earth science and is an important model system for water in a heterogeneous environment. Of particular interest is understanding the extent to which the structure and dynamics, and ultimately reactivity, of water at the interface differ from those in the bulk. For example, how does the distribution of hydrogen bonds differ between interfacial and bulk water? How anisotropic is the orientation of the water molecules at the interface? In terms of dynamics, how do the diffusion constant, rotational relaxation time, and hydrogen-bond rearrangement time vary as the interface is approached? One can also consider vibrational dynamics processes such as energy relaxation and transfer.One important technique for addressing these questions is computer simulation. Models used in these calculations for the water surface range from rigid, fixed-point-charge two-body models (1-3), to fluctuating charge or polarizable models (4, 5), to ab initio molecular dynamics calculations (6-10). Regarding static properties, for example, some effort has been expended toward understanding what fraction of H atoms in the surface layer are hydrogen bonded, and what fraction of molecules do not donate any hydrogen bonds (nondonors or "acceptor-only" molecules) (6, 9). In terms of dynamics, it is generally found that diffusion is faster at the interface than in the bulk (1, 4, 10), and rotational relaxation is also faster (3,6,7,10). On the other hand, two studies with fixed-charge two-body models show that hydrogenbond rearrangement is slower at the interface (2, 3), whereas one study with a fluctuating-charge model shows that hydrogen-bond rearrangement is faster (5). In this latter study the authors conclude that this is generally true for polarizable models.Because of its surface sensitivity, vibrational sum-frequency generation (SFG) spectroscopy (11, 12) has become one of the most powerful experimental techniques for the study of interfaces, including the one separa...