Hidden Isolated OH at the Charged Hydrophobic Interface Revealed by Two-Dimensional Heterodyne-Detected VSFG Spectroscopy

Water around hydrophobic groups mediates hydrophobic interactions that play key roles in many chemical and biological processes. Thus, the molecular-level elucidation of the properties of water in the vicinity of hydrophobic groups is important. We report on the structure and dynamics of water at two oppositely charged hydrophobic ion/water interfaces, that is, the tetraphenylborate-ion (TPB⁻)/water and tetraphenylarsonium-ion (TPA⁺)/water interfaces, which are clarified by two-dimensional heterodyne-detected vibrational sum-frequency generation (2D HD-VSFG) spectroscopy. The obtained 2D HD-VSFG spectra of the anionic TPB⁻ interface reveal the existence of distinct π-hydrogen bonded OH groups in addition to the usual hydrogen-bonded OH groups, which are hidden in the steady-state spectrum. In contrast, 2D HD-VSFG spectra of the cationic TPA⁺ interface only show the presence of usual hydrogen-bonded OH groups. The present study demonstrates that the sign of the interfacial charge governs the structure and dynamics of water molecules that face the hydrophobic region.     

Ultrafast vibrational energy relaxation of the water bridge

We report the energy relaxation of the OH stretch vibration of HDO molecules contained in an HDO:D(2)O water bridge using femtosecond mid-infrared pump-probe spectroscopy. We found that the vibrational lifetime is shorter (~630 ± 50 fs) than for HDO molecules in bulk HDO:D(2)O (~740 ± 40 fs). In contrast, the thermalization dynamics following the vibrational relaxation are much slower (~1.5 ± 0.4 ps) than in bulk HDO:D(2)O (~250 ± 90 fs). These differences in energy relaxation dynamics strongly indicate that the water bridge and bulk water differ on a molecular scale.     

Femtosecond Hydrogen Bond Dynamics of Bulk‐like and Bound Water at Positively and Negatively Charged Lipid Interfaces Revealed by 2D HD‐VSFG Spectroscopy

Interfacial water in the vicinity of lipids plays an important role in many biological processes, such as drug delivery, ion transportation, and lipid fusion. Hence, molecular‐level elucidation of the properties of water at lipid interfaces is of the utmost importance. We report the two‐dimensional heterodyne‐detected vibrational sum frequency generation (2D HD‐VSFG) study of the OH stretch of HOD at charged lipid interfaces, which shows that the hydrogen bond dynamics of interfacial water differ drastically, depending on the lipids. The data indicate that the spectral diffusion of the OH stretch at a positively charged lipid interface is dominated by the ultrafast (<～100 fs) component, followed by the minor sub‐picosecond slow dynamics, while the dynamics at a negatively charged lipid interface exhibit sub‐picosecond dynamics almost exclusively, implying that fast hydrogen bond fluctuation is prohibited. These results reveal that the ultrafast hydrogen bond dynamics at the positively charged lipid–water interface are attributable to the bulk‐like property of interfacial water, whereas the slow dynamics at the negatively charged lipid interface are due to bound water, which is hydrogen‐bonded to the hydrophilic head group.     

Structural inhomogeneity of interfacial water at lipid monolayers revealed by surface-specific vibrational pump-probe spectroscopy

We report vibrational lifetime measurements of the OH stretch vibration of interfacial water in contact with lipid monolayers, using time-resolved vibrational sum frequency (VSF) spectroscopy. The dynamics of water in contact with four different lipids are reported and are characterized by vibrational relaxation rates measured at 3200, 3300, 3400, and 3500 cm(-1). We observe that the water molecules with an OH frequency ranging from 3300 to 3500 cm(-1) all show vibrational relaxation with a time constant of T(1) = 180 ± 35 fs, similar to what is found for bulk water. Water molecules with OH groups near 3200 cm(-1) show distinctly faster relaxation dynamics, with T(1) < 80 fs. We successfully model the data by describing the interfacial water containing two distinct subensembles in which spectral diffusion is, respectively, rapid (3300-3500 cm(-1)) and absent (3200 cm(-1)). We discuss the potential biological implications of the presence of the strongly hydrogen-bonded, rapidly relaxing water molecules at 3200 cm(-1) that are decoupled from the bulk water system.     

Ultrafast vibrational dynamics of interfacial water

We report investigations of the vibrational dynamics of water molecules at the water–air and at the water–lipid interface. Following vibrational excitation with an intense femtosecond infrared pulse resonant with the O–H stretch vibration of water, we follow the subsequent relaxation processes using the surface-specific spectroscopic technique of sum frequency generation. This allows us to selectively follow the vibrational relaxation of the approximately one monolayer of water molecules at the interface. Although the surface vibrational spectra of water at the interface with air and lipids are very similar, we find dramatic variations in both the rates and mechanisms of vibrational relaxation. For water at the water–air interface, very rapid exchange of vibrational energy occurs with water molecules in the bulk, and this intermolecular energy transfer process dominates the response. For membrane-bound water at the lipid interface, intermolecular energy transfer is suppressed, and intramolecular relaxation dominates. The difference in relaxation mechanism can be understood from differences in the local environments experienced by the interfacial water molecules in the two different systems.     

Spectrally- and time-resolved vibrational surface spectroscopy: ultrafast hydrogen-bonding dynamics at D2O/CaF2 interface

Time- and frequency-domain three-wave mixing spectroscopy (IR+visible sum frequency generation) is developed as the lowest-order nonlinear technique that is both surface selective and capable of measuring spectral evolution of vibrational coherences. Using 70 fs infrared and 40 fs visible pulses, we observe ultrafast spectral dynamics of the OD stretch of D2O at the CaF2 surface. Spectral shifts indicative of the hydrogen-bond network rearrangement occur on the 100 fs time scale, within the observation time window determined by the vibrational dephasing. By tuning the IR pulse wavelength to the blue or red side of the OD-stretch transition, we selectively monitor the dynamics of different subensembles in the distribution of the H-bond structures. The blue-side excitation (weaker H-bonding structures) shows monotonic decay and nu(OD) frequency shift to the red on a 100 fs time scale, which is better described by a Gaussian than an exponential frequency correlation function. In contrast, the red-side excitation (stronger H-bonding structures) results in a blue spectral shift and a recursion in the signal at 125+/-10 fs, indicating the presence of an underdamped intermolecular mode of interfacial water.     

Theoretical studies of solute vibrational energy relaxation at liquid interfaces

Recent advances in the theoretical understanding of solute vibrational energy relaxation at liquid interfaces and surfaces are described. Non-equilibrium molecular dynamics simulations of the relaxation of an initially excited solute molecule are combined with equilibrium force autocorrelation calculations to gain insight into the factors that influence the vibrational relaxation rate. Diatomic and triatomic nonpolar, polar, and ionic solute molecules adsorbed at the liquid/vapor interface of several liquids as well as at the water/CCl(4) liquid/liquid interface are considered. In general, the vibrational relaxation rate is significantly slower (a factor of 3 to 4) at the liquid/vapor and liquid/liquid interface than in the bulk due to the reduced density, which gives rise to a reduced contribution of the repulsive solvent-solute forces on the vibrational mode. The surface effects on the ionic solutes are much smaller (50% or less slower relaxation relative to the bulk). This is due to the fact that ionic solutes at the interface are able to keep part of their solvation shell to a degree that depends on their size. Thus, a significant portion of the repulsive forces is maintained. A high degree of correlation is found between the peak height of the solvent-solute radial distribution function and the vibrational relaxation rate. The relaxation rate at the liquid/liquid interface strongly depends on the location of the solute across the interface and correlates with the change in the density and polarity profile of the interface.     

Nonperturbative vibrational energy relaxation effects on vibrational line shapes

A general formulation of nonperturbative quantum dynamics of solutes in a condensed phase is proposed to calculate linear and nonlinear vibrational line shapes. In the weak solute-solvent interaction limit, the temporal absorption profile can be approximately factorized into the population relaxation profile from the off-diagonal coupling and the pure-dephasing profile from the diagonal coupling. The strength of dissipation and the anharmonicity-induced dephasing rate are derived in Appendix A. The vibrational energy relaxation (VER) rate is negligible for slow solvent fluctuations, yet it does not justify the Markovian treatment of off-diagonal contributions to vibrational line shapes. Non-Markovian VER effects are manifested as asymmetric envelops in the temporal absorption profile, or equivalently as side bands in the frequency domain absorption spectrum. The side bands are solvent-induced multiple-photon effects which are absent in the Markovian VER treatment. Exact path integral calculations yield non-Lorentzian central peaks in absorption spectrum resulting from couplings between population relaxations of different vibrational states. These predictions cannot be reproduced by the perturbative or the Markovian approximations. For anharmonic potentials, the absorption spectrum shows asymmetric central peaks and the asymmetry increases with anharmonicity. At large anharmonicities, all the approximation schemes break down and a full nonperturbative path integral calculation that explicitly accounts for the exact VER effects is needed. A numerical analysis of the O-H stretch of HOD in D(2)O solvent reveals that the non-Markovian VER effects generate a small recurrence of the echo peak shift around 200 fs, which cannot be reproduced with a Markovian VER rate. In general, the nonperturbative and non-Markovian VER contributions have a stronger effect on nonlinear vibrational line shapes than on linear absorption.     

A first principles simulation study of fluctuations of hydrogen bonds and vibrational frequencies of water at liquid-vapor interface

We present a theoretical study of the structure and dynamics of water-vapor interface by means of ab initio molecular dynamics simulations. The inhomogeneous density, hydrogen bond and orientational profiles, voids and vibrational frequency distributions are investigated. We have also studied various dynamical properties of the interface such as diffusion, orientational relaxation, hydrogen bond dynamics and vibrational frequency fluctuations. The diffusion and orientational relaxation of water molecules are found to be faster at the interface which can be correlated with the voids present in the system. The hydrogen bond dynamics, however, is found to be slightly slower at the interface than that in bulk water. The correlations of hydrogen bond relaxation with the dynamics of vibrational frequency fluctuations are also discussed.     

Microscopic inhomogeneity and ultrafast orientational motion in an organic photovoltaic bulk heterojunction thin film studied with 2D IR vibrational spectroscopy

Two-dimensional infrared vibrational spectroscopy is used to examine conformational inhomogeneity and ultrafast orientational motion within local environments of an organic photovoltaic bulk heterojunction thin film. The bulk heterojunction material consists of a mixture of the electron donor poly[2-methoxy-5-(2'-ethylhexyloxy)-1,4-(1-cyanovinylene)phenylene] (CN-MEH-PPV) and the electron acceptor [6,6]-phenyl-C(61)-butyric acid methyl ester (PCBM). PCBM species reside in a distribution of environments within large domains of the molecules that cause their C=O stretch modes to be inhomogeneously broadened. The molecular inhomogeneity also results in frequency dependent vibrational relaxation dynamics. The butyric acid methyl ester group of PCBM undergoes ultrafast wobbling-in-the-cone orientational motion on the 110 fs time scale within a cone semiangle of 29 degrees . The vibrational dynamics are sensitive metrics of molecular order in the material and have implications for charge mobility and degradation phenomena in organic photovoltaic devices. This report represents the first study of organic photovoltaic materials using ultrafast two-dimensional infrared vibrational spectroscopy.     

Three distinct water structures at a zwitterionic lipid/water interface revealed by heterodyne-detected vibrational sum frequency generation

Lipid/water interfaces and associated interfacial water are vital for various biochemical reactions, but the molecular-level understanding of their property is very limited. We investigated the water structure at a zwitterionic lipid, phosphatidylcholine, monolayer/water interface using heterodyne-detected vibrational sum frequency generation spectroscopy. Isotopically diluted water was utilized in the experiments to minimize the effect of intra/intermolecular couplings. It was found that the OH stretch band in the Imχ((2)) spectrum of the phosphatidylcholine/water interface exhibits a characteristic double-peaked feature. To interpret this peculiar spectrum of the zwitterionic lipid/water interface, Imχ((2)) spectra of a zwitterionic surfactant/water interface and mixed lipid/water interfaces were measured. The Imχ((2)) spectrum of the zwitterionic surfactant/water interface clearly shows both positive and negative bands in the OH stretch region, revealing that multiple water structures exist at the interface. At the mixed lipid/water interfaces, while gradually varying the fraction of the anionic and cationic lipids, we observed a drastic change in the Imχ((2)) spectra in which spectral features similar to those of the anionic, zwitterionic, and cationic lipid/water interfaces appeared successively. These observations demonstrate that, when the positive and negative charges coexist at the interface, the H-down-oriented water structure and H-up-oriented water structure appear in the vicinity of the respective charged sites. In addition, it was found that a positive Imχ((2)) appears around 3600 cm(-1) for all the monolayer interfaces examined, indicating weakly interacting water species existing in the hydrophobic region of the monolayer at the interface. On the basis of these results, we concluded that the characteristic Imχ((2)) spectrum of the zwitterionic lipid/water interface arises from three different types of water existing at the interface: (1) the water associated with the negatively charged phosphate, which is strongly H-bonded and has a net H-up orientation, (2) the water around the positively charged choline, which forms weaker H-bonds and has a net H-down orientation, and (3) the water weakly interacting with the hydrophobic region of the lipid, which has a net H-up orientation.     

Dynamics of nanoscopic water: vibrational echo and infrared pump-probe studies of reverse micelles

The dynamics of water in nanoscopic pools 1.7-4.0 nm in diameter in AOT reverse micelles were studied with ultrafast infrared spectrally resolved stimulated vibrational echo and pump-probe spectroscopies. The experiments were conducted on the OD hydroxyl stretch of low-concentration HOD in the H2O, providing a direct examination of the hydrogen-bond network dynamics. Pump-probe experiments show that the vibrational lifetime of the OD stretch mode increases as the size of the reverse micelle decreases. These experiments are also sensitive to hydrogen-bond dissociation and reformation dynamics, which are observed to change with reverse micelle size. Spectrally resolved vibrational echo data were obtained at several frequencies. The vibrational echo data are compared to data taken on bulk water and on a 6 M NaCl solution, which is used to examine the role of ionic strength on the water dynamics in reverse micelles. Two types of vibrational echo measurements are presented: the vibrational echo decays and the vibrational echo peak shifts. As the water nanopool size decreases, the vibrational echo decays become slower. Even the largest nanopool (4 nm, approximately 1000 water molecules) has dynamics that are substantially slower than bulk water. It is demonstrated that the slow dynamics in the reverse micelle water nanopools are a result of confinement rather than ionic strength. The data are fit using time-dependent diagrammatic perturbation theory to obtain the frequency-frequency correlation function (FFCF) for each reverse micelle. The results are compared to the FFCF of water and show that the largest differences are in the slowest time scale dynamics. In bulk water, the slowest time scale dynamics are caused by hydrogen-bond network equilibration, i.e., the making and breaking of hydrogen bonds. For the smallest nanopools, the longest time scale component of the water dynamics is approximately 10 times longer than the dynamics in bulk water. The vibrational echo data for the smallest reverse micelle displays a dependence on the detection wavelength, which may indicate that multiple ensembles of water molecules are being observed.     

Vibrational and orientational dynamics of water in aqueous hydroxide solutions

We report the vibrational and orientational dynamics of water molecules in isotopically diluted NaOH and NaOD solutions using polarization-resolved femtosecond vibrational spectroscopy and terahertz time-domain dielectric relaxation measurements. We observe a speed-up of the vibrational relaxation of the O-D stretching vibration of HDO molecules outside the first hydration shell of OH(-) from 1.7 ± 0.2 ps for neat water to 1.0 ± 0.2 ps for a solution of 5 M NaOH in HDO:H(2)O. For the O-H vibration of HDO molecules outside the first hydration shell of OD(-), we observe a similar speed-up from 750 ± 50 fs to 600 ± 50 fs for a solution of 6 M NaOD in HDO:D(2)O. The acceleration of the decay is assigned to fluctuations in the energy levels of the HDO molecules due to charge transfer events and charge fluctuations. The reorientation dynamics of water molecules outside the first hydration shell are observed to show the same time constant of 2.5 ± 0.2 ps as in bulk liquid water, indicating that there is no long range effect of the hydroxide ion on the hydrogen-bond structure of liquid water. The terahertz dielectric relaxation experiments show that the transfer of the hydroxide ion through liquid water involves the simultaneous motion of ~7 surrounding water molecules, considerably less than previously reported for the proton.     

Dynamics of proteins encapsulated in silica sol-gel glasses studied with IR vibrational echo spectroscopy

Spectrally resolved infrared stimulated vibrational echo spectroscopy is used to measure the fast dynamics of heme-bound CO in carbonmonoxy-myoglobin (MbCO) and -hemoglobin (HbCO) embedded in silica sol-gel glasses. On the time scale of approximately 100 fs to several picoseconds, the vibrational dephasing of the heme-bound CO is measurably slower for both MbCO and HbCO relative to that of aqueous protein solutions. The fast structural dynamics of MbCO, as sensed by the heme-bound CO, are influenced more by the sol-gel environment than those of HbCO. Longer time scale structural dynamics (tens of picoseconds), as measured by the extent of spectral diffusion, are the same for both proteins encapsulated in sol-gel glasses compared to that in aqueous solutions. A comparison of the sol-gel experimental results to viscosity-dependent vibrational echo data taken on various mixtures of water and fructose shows that the sol-gel-encapsulated MbCO exhibits dynamics that are the equivalent of the protein in a solution that is nearly 20 times more viscous than bulk water. In contrast, the HbCO dephasing in the sol-gel reflects only a 2-fold increase in viscosity. Attempts to alter the encapsulating pore size by varying the molar ratio of silane precursor to water (R value) used to prepare the sol-gel glasses were found to have no effect on the fast or steady-state spectroscopic results. The vibrational echo data are discussed in the context of solvent confinement and protein-pore wall interactions to provide insights into the influence of a confined environment on the fast structural dynamics experienced by a biomolecule.     

Multidimensional infrared spectroscopy of water. I. Vibrational dynamics in two-dimensional IR line shapes

In this and the following paper, we describe the ultrafast structural fluctuations and rearrangements of the hydrogen bonding network of water using two-dimensional (2D) infrared spectroscopy. 2D IR spectra covering all the relevant time scales of molecular dynamics of the hydrogen bonding network of water were studied for the OH stretching absorption of HOD in D2O. Time-dependent evolution of the 2D IR line shape serves as a spectroscopic observable that tracks how different hydrogen bonding environments interconvert while changes in spectral intensity result from vibrational relaxation and molecular reorientation of the OH dipole. For waiting times up to the vibrational lifetime of 700 fs, changes in the 2D line shape reflect the spectral evolution of OH oscillators induced by hydrogen bond dynamics. These dynamics, characterized through a set of 2D line shape analysis metrics, show a rapid 60 fs decay, an underdamped oscillation on a 130 fs time scale induced by hydrogen bond stretching, and a long time decay constant of 1.4 ps. 2D surfaces for waiting times larger than 700 fs are dominated by the effects of vibrational relaxation and the thermalization of this excess energy by the solvent bath. Our modeling based on fluctuations with Gaussian statistics is able to reproduce the changes in dispersed pump-probe and 2D IR spectra induced by these relaxation processes, but misses the asymmetry resulting from frequency-dependent spectral diffusion. The dynamical origin of this asymmetry is discussed in the companion paper.     

Water dynamics--the effects of ions and nanoconfinement

Hydrogen bond dynamics of water in highly concentrated NaBr salt solutions and reverse micelles are studied using ultrafast 2D-IR vibrational echo spectroscopy and polarization-selective IR pump-probe experiments performed on the OD hydroxyl stretch of dilute HOD in H(2)O. The vibrational echo experiments measure spectral diffusion, and the pump-probe experiments measure orientational relaxation. Both experimental observables are directly related to the structural dynamics of water's hydrogen bond network. The measurements performed on NaBr solutions as a function of concentration show that the hydrogen bond dynamics slow as the NaBr concentration increases. The most pronounced change is in the longest time scale dynamics which are related to the global rearrangement of the hydrogen bond structure. Complete hydrogen bond network randomization slows by a factor of approximately 3 in approximately 6 M NaBr solution compared to that in bulk water. The hydrogen bond dynamics of water in nanoscopically confined environments are studied by encapsulating water molecules in ionic head group (AOT) and nonionic head group (Igepal CO 520) reverse micelles. Water dynamics in the nanopools of AOT reverse micelles are studied as a function of size by observing orientational relaxation. Orientational relaxation dynamics deviate significantly from bulk water when the size of the reverse micelles is smaller than several nm and become nonexponential and slower as the size of the reverse micelles decreases. In the smallest reverse micelles, orientational relaxation (hydrogen bond structural randomization) is almost 20 times slower than that in bulk water. To determine if the changes in dynamics from bulk water are caused by the influence of the ionic head groups of AOT or the nanoconfinement, the water dynamics in 4 nm nanopools in AOT reverse micelles (ionic) and Igepal reverse micelles (nonionic) are compared. It is found that the water orientational relaxation in the 4 nm diameter nanopools of the two types of reverse micelles is almost identical, which indicates that confinement by an interface to form a nanoscopic water pool is a primary factor governing the dynamics of nanoscopic water rather than the presence of charged groups at the interface.     

Probing dynamics of complex molecular systems with ultrafast 2D IR vibrational echo spectroscopy

Ultrafast 2D IR vibrational echo spectroscopy is described and a number of experimental examples are given. Details of the experimental method including the pulse sequence, heterodyne detection, and determination of the absorptive component of the 2D spectrum are outlined. As an initial example, the 2D spectrum of the stretching mode of CO bound to the protein myoglobin (MbCO) is presented. The time dependence of the 2D spectrum of MbCO, which is caused by protein structural evolution, is presented and its relationship to the frequency-frequency correlation function is described and used to make protein structural assignments based on comparisons to molecular dynamics simulations. The 2D vibrational echo experiments on the protein horseradish peroxidase are presented. The time dependence of the 2D spectra of the enzyme in the free form and with a substrate bound at the active site are compared and used to examine the influence of substrate binding on the protein's structural dynamics. The application of 2D vibrational echo spectroscopy to the study of chemical exchange under thermal equilibrium conditions is described. 2D vibrational echo chemical exchange spectroscopy is applied to the study of formation and dissociation of organic solute-solvent complexes and to the isomerization around a carbon-carbon single bond of an ethane derivative.