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...

Understanding the structure of water near cell membranes is crucial for characterizing water-mediated events such as molecular transport. To obtain structural information of water near a membrane, it is useful to have a surface-selective technique that can probe only interfacial water molecules. One such technique is vibrational sum-frequency generation (VSFG) spectroscopy. As model systems for studying membrane headgroup/water interactions, in this paper we consider lipid and surfactant monolayers on water. We adopt a theoretical approach combining molecular dynamics simulations and phase-sensitive VSFG to investigate water structure near these interfaces. Our simulated spectra are in qualitative agreement with experiments and reveal orientational ordering of interfacial water molecules near cationic, anionic, and zwitterionic interfaces. OH bonds of water molecules point toward an anionic interface leading to a positive VSFG peak, whereas the water hydrogen atoms point away from a cationic interface leading to a negative VSFG peak. Coexistence of these two interfacial water species is observed near interfaces between water and mixtures of cationic and anionic lipids, as indicated by the presence of both negative and positive peaks in their VSFG spectra. In the case of a zwitterionic interface, OH orientation is toward the interface on the average, resulting in a positive VSFG peak.

Self-assembled monolayers of 1, 4-phenylene diisocyanide (PDI) were formed on Au and Pt-group transition metals and examined by surface-enhanced Raman spectroscopy under controlled applied potential. On all of the metals examined, PDI adsorbs in an edge-on manner, with one NC group bound to the surface and the other pointing away from the surface. The N-C stretching frequency (nu(NC)) suggests that depending on the metal, PDI adsorbs on different binding sites: terminal sites on Au, both terminal and bridging on Rh and Pt, and predominantly 3-fold hollow sites for Pd. This binding site preference can be understood in terms of the difference in d-band center energy and d-orbital filling among the metals. The applied potential affects the N-C bonding differently as inferred from the potential dependence of nu(NC). On Au, Rh, and Pd, the nu(NC) increases linearly with the applied potential, yielding a Stark tuning slope, dnu(NC)/dE, of 25, 12, and 10 cm(-1)/V, respectively. On Pt, the nu(NC) is nearly independent of the applied potential. On all of the metals studied, the frequencies of benzene ring vibration modes are not dependent on the applied potential, consistent with the edge-on orientation in which the ring does not directly interact with the surface. Several ring vibrations are, however, sensitive to the nature of metal substrate due to different binding sites involved. The ability of the free NC group to function as an anchoring point is demonstrated by the attachment of gold nanoparticles on PDI-covered Au and Pd. The study provides useful NC-metal bonding information for isocyanide-based molecular electronic developments.

Application of the Herman-Kluk semiclassical propagator to the calculation of spectroscopic response functions for anharmonic oscillators has demonstrated the quantitative accuracy of these approximate dynamics. In this approach, spectroscopic response functions are expressed as multiple phase-space integrals over pairs of classical trajectories and their associated stability matrices. Here we analyze the Herman-Kluk semiclassical approximation to a linear response function and determine the origin of the capacity of this method to reproduce quantum effects in a response function from classical dynamical information. Our analysis identifies those classical trajectories that contribute most significantly to the response function on different time scales. This finding motivates a procedure for computing the linear response function in which the interference between pairs of classical trajectories is treated approximately, resulting in an integral over a single average trajectory, as in a purely classical calculation.