Formation and dissociation of intra-intermolecular hydrogen-bonded solute-solvent complexes: chemical exchange two-dimensional infrared vibrational echo spectroscopy

2-Methoxyphenol (2MP) solutes form weak complexes with toluene solvent molecules. The complexes are unusual in that the 2MP hydroxyl has an intramolecular hydrogen bond and simultaneously forms an intermolecular hydrogen bond with toluene and other aromatic solvents. In the equilibrated solute-solvent solution, there exists approximately the same concentration of 2MP-toluene complex and free 2MP. The very fast formation and dissociation (chemical exchange) of this type of three-centered hydrogen bond complex were observed in real time under thermal equilibrium conditions with two-dimensional (2D) infrared vibrational echo spectroscopy. Chemical exchange is manifested in the 2D spectrum by the growth of off-diagonal peaks. Both the formation and dissociation can be characterized in terms of the dissociation time constant, which was determined to be 3 ps for the 2MP-toluene complex. The intra-intermolecular hydrogen bond formation is influenced by subtle details of the molecular structure. Although 2MP forms a complex with toluene, it is demonstrated that 2-ethoxyphenol (2EP) does not form complexes to any significant extent. Density functional calculations at the B3LYP/6-31+G(d,p) level suggest that steric effects caused by the extra methyl group in 2EP are responsible for the difference.     

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.     

Frequency-frequency correlation functions and apodization in two-dimensional infrared vibrational echo spectroscopy: a new approach

Ultrafast two-dimensional infrared (2D-IR) vibrational echo spectroscopy can probe structural dynamics under thermal equilibrium conditions on time scales ranging from femtoseconds to approximately 100 ps and longer. One of the important uses of 2D-IR spectroscopy is to monitor the dynamical evolution of a molecular system by reporting the time dependent frequency fluctuations of an ensemble of vibrational probes. The vibrational frequency-frequency correlation function (FFCF) is the connection between the experimental observables and the microscopic molecular dynamics and is thus the central object of interest in studying dynamics with 2D-IR vibrational echo spectroscopy. A new observable is presented that greatly simplifies the extraction of the FFCF from experimental data. The observable is the inverse of the center line slope (CLS) of the 2D spectrum. The CLS is the inverse of the slope of the line that connects the maxima of the peaks of a series of cuts through the 2D spectrum that are parallel to the frequency axis associated with the first electric field-matter interaction. The CLS varies from a maximum of 1 to 0 as spectral diffusion proceeds. It is shown analytically to second order in time that the CLS is the T(w) (time between pulses 2 and 3) dependent part of the FFCF. The procedure to extract the FFCF from the CLS is described, and it is shown that the T(w) independent homogeneous contribution to the FFCF can also be recovered to yield the full FFCF. The method is demonstrated by extracting FFCFs from families of calculated 2D-IR spectra and the linear absorption spectra produced from known FFCFs. Sources and magnitudes of errors in the procedure are quantified, and it is shown that in most circumstances, they are negligible. It is also demonstrated that the CLS is essentially unaffected by Fourier filtering methods (apodization), which can significantly increase the efficiency of data acquisition and spectral resolution, when the apodization is applied along the axis used for obtaining the CLS and is symmetrical about tau=0. The CLS is also unchanged by finite pulse durations that broaden 2D spectra.     

Taking apart the two-dimensional infrared vibrational echo spectra: more information and elimination of distortions

Ultrafast two-dimensional infrared (2D-IR) vibrational echo spectroscopy can probe the fast structural evolution of molecular systems under thermal equilibrium conditions. Structural dynamics are tracked by observing the time evolution of the 2D-IR spectrum, which is caused by frequency fluctuations of vibrational mode(s) excited during the experiment. However, there are a variety of effects that can produce line shape distortions and prevent the correct determination of the frequency-frequency correlation function (FFCF), which describes the frequency fluctuations and connects the experimental observables to a molecular level depiction of dynamics. In addition, it can be useful to analyze different parts of the 2D spectrum to determine if dynamics are different for subensembles of molecules that have different initial absorption frequencies in the inhomogeneously broadened absorption line. Here, an important extension to a theoretical method for extraction of the FFCF from 2D-IR spectra is described. The experimental observable is the center line slope (CLSomega(m)) of the 2D-IR spectrum. The CLSomega(m) is obtained by taking slices through the 2D spectrum parallel to the detection frequency axis (omega(m)). Each slice is a spectrum. The slope of the line connecting the frequencies of the maxima of the sliced spectra is the CLSomega(m). The change in slope of the CLSomega(m) as a function of time is directly related to the FFCF and can be used to obtain the complete FFCF. CLSomega(m) is immune to line shape distortions caused by destructive interference between bands arising from vibrational echo emission, from the 0-1 vibrational transition (positive), and from the 1-2 vibrational transition (negative) in the 2D-IR spectrum. The immunity to the destructive interference enables the CLSomega(m) method to compare different parts of the bands as well as comparing the 0-1 and 1-2 bands. Also, line shape distortions caused by solvent background absorption and finite pulse durations do not affect the determination of the FFCF with the CLSomega(m) method. The CLSomega(m) can also provide information on the cross correlation between frequency fluctuations of the 0-1 and 1-2 vibrational transitions.     

Solute-solvent complex kinetics and thermodynamics probed by 2D-IR vibrational echo chemical exchange spectroscopy

The formation and dissociation kinetics of a series of triethylsilanol/solvent weakly hydrogen bonding complexes with enthalpies of formation ranging from -1.4 to -3.3 kcal/mol are measured with ultrafast two-dimensional infrared (2D IR) chemical exchange spectroscopy in liquid solutions at room temperature. The correlation between the complex enthalpies of formation and dissociation rate constants can be expressed with an equation similar to the Arrhenius equation. The experimental results are in accord with previous observations on eight phenol/solvent complexes with enthalpies of formation from -0.6 to -2.5 kcal/mol. It was found that the inverse of the solute-solvent complex dissociation rate constant is linearly related to exp(-DeltaH0/RT) where DeltaH0 is the complex enthalpy of formation. It is shown here, that the triethylsilanol-solvent complexes obey the same relationship with the identical proportionality constant, that is, all 13 points, five silanol complexes and eight phenol complexes, fall on the same line. In addition, features of 2D IR chemical exchange spectra at long reaction times (spectral diffusion complete) are explicated using the triethylsilanol systems. It is shown that the off-diagonal chemical exchange peaks have shapes that are a combination (outer product) of the absorption line shapes of the species that give rise to the diagonal peaks.     

Accidental vibrational degeneracy in vibrational excited states observed with ultrafast two-dimensional IR vibrational echo spectroscopy

The coupling between the OD stretch v=2 level and benzene-ring modes in 2-methoxyphenol-OD (hydroxyl H replaced by D) is observed with ultrafast two-dimensional (2D) IR vibrational echo spectroscopy. Because of this coupling, the 1-2 transition peak in the 2D spectrum is split into a doublet with peaks of approximately equal amplitudes. Several molecules and solvents were used to study this phenomenon. Near-IR (NIR) spectroscopy measurements and density-functional theory calculations (B3LYP6-31+G(d,p) level) were also applied. Experimental results and calculations show that the OD stretch 1-2 transition is coupled to a combination band related to the benzene-ring motions. A simple quantum-mechanical model indicates that the combination band has a frequency of 5172 and 5176.5 cm(-1) in CCl4 and hexane, respectively. The transition between this combination band and the ground state is too weak to detect by NIR. The transition between this band and the OD stretch first excited state is also so weak that most of the intensity of the doublet comes from the oscillator strength produced by coupling to the OD stretch. The model gives the coupling strengths as 6.5 and 7 cm(-1) in CCl4 and hexane, respectively.     

Site-specific vibrational dynamics of the CD3zeta membrane peptide using heterodyned two-dimensional infrared photon echo spectroscopy

Heterodyned two-dimensional infrared (2D IR) spectroscopy has been used to study the amide I vibrational dynamics of a 27-residue peptide in lipid vesicles that encompasses the transmembrane domain of the T-cell receptor CD3zeta. Using 1-(13)C[Double Bond](18)O isotope labeling, the amide I mode of the 49-Leucine residue was spectroscopically isolated and the homogeneous and inhomogeneous linewidths of this mode were measured by fitting the 2D IR spectrum collected with a photon echo pulse sequence. The pure dephasing and inhomogeneous linewidths are 2 and 32 cm(-1), respectively. The population relaxation time of the amide I band was measured with a transient grating, and it contributes 9 cm(-1) to the linewidth. Comparison of the 49-Leucine amide I mode and the amide I band of the entire CD3zeta peptide reveals that the vibrational dynamics are not uniform along the length of the peptide. Possible origins for the large amount of inhomogeneity present at the 49-Leucine site are discussed.     

Calculation of two-dimensional infrared spectra of ultrafast chemical exchange with numerical Langevin simulations

We combine numerical Langevin simulations with numerical integration of the Schrodinger equation to calculate two-dimensional infrared spectra of ultrafast chemical exchange. This provides a tool to model and interpret such spectra of molecules undergoing chemical processes, such as isomerization and solvent exchange reactions. Two-dimensional infrared spectroscopy has already been used to extract reaction rates for ultrafast chemical reactions. We demonstrate that these spectra are not only sensitive to the rates, but also to the finite duration of the exchange. This is emphasised by comparing with the popular Kubo two-state jump models, which do not account for finite exchange times.     

Ion-pairing dynamics of Li+ and SCN- in dimethylformamide solution: chemical exchange two-dimensional infrared spectroscopy

Ultrafast two-dimensional infrared (2DIR) spectroscopy has been proven to be an exceptionally useful method to study chemical exchange processes between different vibrational chromophores under thermal equilibria. Here, we present experimental results on the thermal equilibrium ion pairing dynamics of Li(+) and SCN(-) ions in N,N-dimethylformamide. Li(+) and SCN(-) ions can form a contact ion pair (CIP). Varying the relative concentration of Li(+) in solution, we could control the equilibrium CIP and free SCN(-) concentrations. Since the CN stretch frequency of Li-SCN CIP is blue-shifted by about 16 cm(-1) from that of free SCN(-) ion, the CN stretch IR spectrum is a doublet. The temperature-dependent IR absorption spectra reveal that the CIP formation is an endothermic (0.57 kJ∕mol) process and the CIP state has larger entropy by 3.12 J∕(K?mol) than the free ion states. Since the two ionic configurations are spectrally distinguishable, this salt solution is ideally suited for nonlinear IR spectroscopic investigations to study ion pair association and dissociation dynamics. Using polarization-controlled IR pump-probe methods, we first measured the lifetimes and orientational relaxation times of these two forms of ionic configurations. The vibrational population relaxation times of both the free ion and CIP are about 32 ps. However, the orientational relaxation time of the CIP, which is ～47 ps, is significantly longer than that of the free SCN(-), which is ～7.7 ps. This clearly indicates that the effective moment of inertia of the CIP is much larger than that of the free SCN(-). Then, using chemical exchange 2DIR spectroscopy and analyzing the diagonal peak and cross-peak amplitude changes with increasing the waiting time, we determined the contact ion pair association and dissociation time constants that are found to be 165 and 190 ps, respectively. The results presented and discussed in this paper are believed to be important, not only because the ion-pairing dynamics is one of the most fundamental physical chemistry problems but also because such molecular ion-ion interactions are of critical importance in understanding Hofmeister effects on protein stability.     

Ultrafast vibrational dynamics of water at a charged interface revealed by two-dimensional heterodyne-detected vibrational sum frequency generation

Two-dimensional heterodyne-detected vibrational sum frequency generation (2D HD-VSFG) spectroscopy is performed for an aqueous interface for the first time. The 2D HD-VSFG spectra in the OH stretch region are obtained from a positively charged surfactant∕water interface with isotopically diluted water (HOD∕D(2)O) to reveal the femtosecond vibrational dynamics of water at the charged interface. The 2D HD-VSFG spectrum is diagonally elongated immediately after photoexcitation, clearly demonstrating inhomogeneity in the interfacial water. This elongation almost disappears at 300 fs owing to the spectral diffusion. Interestingly, the 2D HD-VSFG spectrum at the 0 fs shows an oppositely asymmetric shape to the corresponding 2D IR spectrum in bulk water: The bandwidth of the bleach signal gets narrower when the pump wavenumber becomes higher. This suggests that the dynamics and mechanism of the hydrogen bond rearrangement at the charged interface are significantly different from those in bulk water.     

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.     

Ultrafast infrared spectroscopy of riboflavin: dynamics, electronic structure, and vibrational mode analysis

Femtosecond time-resolved infrared spectroscopy was used to study the vibrational response of riboflavin in DMSO to photoexcitation at 387 nm. Vibrational cooling in the excited electronic state is observed and characterized by a time constant of 4.0 +/- 0.1 ps. Its characteristic pattern of negative and positive IR difference signals allows the identification and determination of excited-state vibrational frequencies of riboflavin in the spectral region between 1100 and 1740 cm (-1). Density functional theory (B3LYP), Hartree-Fock (HF) and configuration interaction singles (CIS) methods were employed to calculate the vibrational spectra of the electronic ground state and the first singlet excited pipi* state as well as respective electronic energies, structural parameters, electronic dipole moments and intrinsic force constants. The harmonic frequencies of the S 1 excited state calculated by the CIS method are in satisfactory agreement with the observed band positions. There is a clear correspondence between computed ground- and excited-state vibrations. Major changes upon photoexcitation include the loss of the double bond between the C4a and N5 atoms, reflected in a downshift of related vibrations in the spectral region from 1450 to 1720 cm (-1). Furthermore, the vibrational analysis reveals intra- and intermolecular hydrogen bonding of the riboflavin chromophore.     

Retrieval of spectral and dynamic properties from two-dimensional infrared pump-probe experiments

We have developed a fitting algorithm able to extract spectral and dynamic properties of a three level oscillator from a two-dimensional infrared spectrum (2D-IR) detected in time resolved nonlinear experiments. Such properties go from the frequencies of the ground-to-first and first-to-second vibrational transitions (and hence anharmonicity) to the frequency-fluctuation correlation function. This last is represented through a general expression that allows one to approach the various strategies of modeling proposed in the literature. The model is based on the Kubo picture of stochastic fluctuations of the transition frequency as a result of perturbations by a fluctuating surrounding. To account for the line-shape broadening due to pump pulse spectral width in double-resonance measurements, we supply the fitting algorithm with the option to perform the convolution of the spectral signal with a Lorentzian function in the pump-frequency dimension. The algorithm is tested here on 2D-IR pump-probe spectra of a Gly-Ala dipeptide recorded at various pump-probe delay times. Speedup benchmarks have been performed on a small Beowulf cluster. The program is written in FORTRAN language for both serial and parallel architectures and is available free of charge to the interested reader.     

Probing intermolecular couplings in liquid water with two-dimensional infrared photon echo spectroscopy

Two-dimensional infrared photon echo and pump probe studies of the OH stretch vibration provide a sensitive probe of the correlations and couplings in the hydrogen bond network of liquid water. The nonlinear response is simulated using numerical integration of the Schrodinger equation with a Hamiltonian constructed to explicitly treat intermolecular coupling and nonadiabatic effects in the highly disordered singly and doubly excited vibrational exciton manifolds. The simulated two-dimensional spectra are in close agreement with our recent experimental results. The high sensitivity of the OH stretch vibration to the bath dynamics is found to arise from intramolecular mixing between states in the two-dimensional anharmonic OH stretch potential. Surprisingly small intermolecular couplings reproduce the experimentally observed intermolecular energy transfer times.     

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.     

Signatures of beta-peptide unfolding in two-dimensional vibrational echo spectroscopy: a simulation study.

An ensemble of exciton Hamiltonians for the amide-I band of the folded and unfolded states of a helical beta-heptapeptide is generated using a molecular dynamics (MD) simulation. The correlated fluctuations of its parameters and their signatures in two-dimensional (2D) vibrational echo spectroscopy are computed. This technique uses infrared pulse sequences to provide ultrafast snapshots of molecular structural fluctuations, in analogy with multidimensional NMR. The present study demonstrates that, by combining a method of calculating the vibrational Hamiltonian from MD snapshots and the nonlinear exciton equations (NEE), it may be possible to simulate realistic multidimensional IR spectra of chemically and biologically interesting systems.     

Ultrafast intermolecular dynamics of liquid water: a theoretical study on two-dimensional infrared spectroscopy

Physical and chemical properties of liquid water are dominated by hydrogen bond structure and dynamics. Recent studies on nonlinear vibrational spectroscopy of intramolecular motion provided new insight into ultrafast hydrogen bond dynamics. However, our understanding of intermolecular dynamics of water is still limited. We theoretically investigated the intermolecular dynamics of liquid water in terms of two-dimensional infrared (2D IR) spectroscopy. The 2D IR spectrum of intermolecular frequency region (<1000 cm(-1)) is calculated by using the equilibrium and nonequilibrium hybrid molecular dynamics method. We find the ultrafast loss of the correlation of the libration motion with the time scale of approximately 110 fs. It is also found that the energy relaxation from the libration motion to the low frequency motion takes place with the time scale of about 180 fs. We analyze the effect of the hindered translation motion on these ultrafast dynamics. It is shown that both the frequency modulation of libration motion and the energy relaxation from the libration to the low frequency motion significantly slow down in the absence of the hindered translation motion. The present result reveals that the anharmonic coupling between the hindered translation and libration motions is essential for the ultrafast relaxation dynamics in liquid water.     

Theoretical calculations of infrared absorption, vibrational circular dichroism, and two-dimensional vibrational spectra of acetylproline in liquids water and chloroform

Infrared absorption, vibrational circular dichroism, and two-dimensional infrared pump-probe and photon echo spectra of acetylproline solutions are theoretically calculated and directly compared with experiments. In order to quantitatively determine interpeptide interaction-induced amide I mode frequency shifts, high-level quantum chemistry calculations were performed. The solvatochromic amide I mode frequency shift and fluctuation were taken into account by carrying out molecular dynamics simulations of acetylproline dissolved in liquids water and chloroform and by using the extrapolation method developed recently. We first studied correlation time scales of the two amide I vibrational frequency fluctuations, cross correlation between the two fluctuating local mode frequencies, ensemble averaged conformations of the acetylproline molecule in liquids water and chloroform. The corresponding conformations of the acetylproline in liquids water and chloroform are close to the ideal 3(10) helix and the C(7) structure, respectively. A few methods proposed to determine the angle between the two transition dipoles associated with the amide I vibrations were tested and their limitations are discussed.     

Water dynamics in small reverse micelles in two solvents: two-dimensional infrared vibrational echoes with two-dimensional background subtraction

Water dynamics as reflected by the spectral diffusion of the water hydroxyl stretch were measured in w(0) = 2 (1.7 nm diameter) Aerosol-OT (AOT)/water reverse micelles in carbon tetrachloride and in isooctane solvents using ultrafast 2D IR vibrational echo spectroscopy. Orientational relaxation and population relaxation are observed for w(0) = 2, 4, and 7.5 in both solvents using IR pump-probe measurements. It is found that the pump-probe observables are sensitive to w(0), but not to the solvent. However, initial analysis of the vibrational echo data from the water nanopool in the reverse micelles in the isooctane solvent seems to yield different dynamics than the CCl(4) system in spite of the fact that the spectra, vibrational lifetimes, and orientational relaxation are the same in the two systems. It is found that there are beat patterns in the interferograms with isooctane as the solvent. The beats are observed from a signal generated by the AOT/isooctane system even when there is no water in the system. A beat subtraction data processing procedure does a reasonable job of removing the distortions in the isooctane data, showing that the reverse micelle dynamics are the same within experimental error regardless of whether isooctane or carbon tetrachloride is used as the organic phase. Two time scales are observed in the vibrational echo data, ~1 and ~10 ps. The slower component contains a significant amount of the total inhomogeneous broadening. Physical arguments indicate that there is a much slower component of spectral diffusion that is too slow to observe within the experimental window, which is limited by the OD stretch vibrational lifetime.     

DNA vibrational coupling revealed with two-dimensional infrared spectroscopy: insight into why vibrational spectroscopy is sensitive to DNA structure

Two-dimensional infrared (2D IR) spectroscopy was used to study the carbonyl vibrational modes of guanine and cytosine bases in A- and B-form DNA. Located between 1600 and 1700 cm(-1), these modes are often used to monitor DNA secondary structure with traditional infrared spectroscopies such as FTIR, but traditional spectroscopies lack the necessary observables to unravel the coupling mechanisms that make these modes sensitive to secondary structure. By using 2D IR spectroscopy and electronic structure calculations on d(G(5)C(5)) and d(GC)(8) model nucleic acids, we find that hydrogen-bonded guanine/cytosine base pairs are primarily electrostatically coupled and that the coupling between these modes can be modeled with a transition dipole density approach. In comparison, electrostatics is insufficient to model stacked bases because of cooperative charge-sharing effects, but the coupling can be accurately calculated using a finite difference method. We find that the coupling is very strong for both hydrogen-bonded and stacked base geometries, creating vibrational modes that extend both across the base pairs and along the lengths of the helices. Our results provide a physical basis for understanding how strong coupling gives rise to the empirically established relationship between infrared spectroscopy and DNA/RNA secondary structure.