Ultrafast two-dimensional infrared vibrational echo chemical exchange experiments and theory

Ultrafast two-dimensional (2D) infrared vibrational echo experiments and theory are used to examine chemical exchange between solute-solvent complexes and the free solute for the solute phenol and three solvent complex partners, p-xylene, benzene, and bromobenzene, in mixed solvents of the partner and CCl4. The experiments measure the time evolution of the 2D spectra of the hydroxyl (OD) stretching mode of the phenol. The time-dependent 2D spectra are analyzed using time-dependent diagrammatic perturbation theory with a model that includes the chemical exchange (formation and dissociation of the complexes), spectral diffusion of both the complex and the free phenol, orientational relaxation of the complexes and free phenol, and the vibrational lifetimes. The detailed calculations are able to reproduce the experimental results and demonstrate that a method employed previously that used a kinetic model for the volumes of the peaks is adequate to extract the exchange kinetics. The current analysis also yields the spectral diffusion (time evolution of the dynamic line widths) and shows that the spectral diffusion is significantly different for phenol complexes and free phenol.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Stochastic simulation of chemical exchange in two dimensional infrared spectroscopy

The stochastic Liouville equations are employed to investigate the combined signatures of chemical exchange (two-state jump) and spectral diffusion (coupling to an overdamped Brownian oscillator) in the coherent response of an anharmonic vibration to three femtosecond infrared pulses. Simulations reproduce the main features recently observed in the OD stretch of phenol in benzene.     

Relating linear vibrational spectroscopy to condensed-phase hydrogen-bonded structures: Liquid-to-supercritical water

The pressure and temperature-dependent linear absorption spectrum of partially deuterated water HOD dissolved in heavy water D(2)O was measured in the OH-stretching spectral region. The temperature was varied in the interval of 298 K相似文献   

Interpreting DNA vibrational circular dichroism spectra using a coupling model from two-dimensional infrared spectroscopy

Two-dimensional infrared spectroscopy was recently used to measure the vibrational couplings between carbonyl bonds located on DNA nucleobases (Krummel, A. T.; Mukherjee, P.; Zanni, M. T. J. Phys. Chem. B 2003, 107, 9165 and Krummel, A. T.; Zanni, M. T. J. Phys. Chem. B 2006, 110, 13991). Here, we extend the coupling model derived from these 2D IR experiments to simulate the vibrational absorption and vibrational circular dichroism (VCD) spectra of three double-stranded DNA oligomers: poly(dG)-poly(dC), poly(dG-dC), and dGGCC. Using this model, we determine that the VCD spectrum of A-form poly(dG)-poly(dC) is dominated by interactions between stacked bases, whereas the coupling between base pairs and stacked bases carries equal importance in the VCD spectrum of B-form poly(dG-dC). We also simulate the absorption and VCD spectra of dGGCC, which is a combination of A- and B-form configurations. These simulations give insight into the structural interpretation of VCD and absorption spectroscopies that have long been used to monitor DNA secondary structure and kinetics.     

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.     

Dynamics of water probed with vibrational echo correlation spectroscopy

Vibrational echo correlation spectroscopy experiments on the OD stretch of dilute HOD in H(2)O are used to probe the structural dynamics of water. A method is demonstrated for combining correlation spectra taken with different infrared pulse bandwidths (pulse durations), making it possible to use data collected from many experiments in which the laser pulse properties are not identical. Accurate measurements of the OD stretch anharmonicity (162 cm(-1)) are presented and used in the data analysis. In addition, the recent accurate determination of the OD vibrational lifetime (1.45 ps) and the time scale for the production of vibrational relaxation induced broken hydrogen bond "photoproducts" ( approximately 2 ps) aid in the data analysis. The data are analyzed using time dependent diagrammatic perturbation theory to obtain the frequency time correlation function (FTCF). The results are an improved FTCF compared to that obtained previously with vibrational echo correlation spectroscopy. The experimental data and the experimentally determined FTCF are compared to calculations that employ a polarizable water model (SPC-FQ) to calculate the FTCF. The SPC-FQ derived FTCF is much closer to the experimental results than previously tested nonpolarizable water models which are also presented for comparison.     

Water penetration into protein secondary structure revealed by hydrogen-deuterium exchange two-dimensional infrared spectroscopy

Two-dimensional infrared spectroscopy in conjunction with hydrogen-deuterium exchange experiments provides detailed information about solvent penetration into protein structure. Correlating the secondary-structure sensitivity of the amide I vibration and the solvent-exposure sensitivity of amide II provides a direct probe of solvent-inaccessible residues of proteins embedded in the hydrophobic core or those involved in strong hydrogen bonds in secondary structures. Distinct spectral signatures of the cross-peak region arising from the coupling of the amide I and II modes imply a significant degree of structural stability of hydrogen-bonded contacts in alpha-helices and beta-sheets in a series of proteins. Ubiquitin, an alpha/beta-protein, exhibits strong alpha-helical signatures and lacks those of the beta-sheet in the cross-peak region, demonstrating that ubiquitin's beta-sheet exchanges protons with the surrounding solvent and is conformationally unstable.     

Solid-state vibrational spectroscopy: Part 11. An infrared and raman vibrational spectroscopic study of metal(II) halide aniline complexes

Infrared (4000?200 cm^?1) and Raman (3500?50 cm^?1) spectra are reported for metal(II) halide aniline complexes of the following stoichiometries: (MX₂an₂) (M  Co, Ni or Hg, X  Cl; M  Mn, X  Cl or Br; M  Zn or Cd, X  Cl, Br or I); (MX₂an₃) (M  Mn, X  Cl or Br; M  Ni, X  Cl); (CdCl₂an) and an assignment is proposed for all the observed bands. Low-temperature (83 K) IR spectra are also reported and it is noted that whilst the aniline ring and CH mode values are virtually insensitive to temperature, the NH₂ rocking and metal-ligand stretching mode values increase with decreasing temperature, whilst the NH₂ stretching mode values decrease with decreasing temperature.     

Hydrogen bond lifetimes and energetics for solute/solvent complexes studied with 2D-IR vibrational echo spectroscopy

Weak pi hydrogen-bonded solute/solvent complexes are studied with ultrafast two-dimensional infrared (2D-IR) vibrational echo chemical exchange spectroscopy, temperature-dependent IR absorption spectroscopy, and density functional theory calculations. Eight solute/solvent complexes composed of a number of phenol derivatives and various benzene derivatives are investigated. The complexes are formed between the phenol derivative (solute) in a mixed solvent of the benzene derivative and CCl4. The time dependence of the 2D-IR vibrational echo spectra of the phenol hydroxyl stretch is used to directly determine the dissociation and formation rates of the hydrogen-bonded complexes. The dissociation rates of the weak hydrogen bonds are found to be strongly correlated with their formation enthalpies. The correlation can be described with an equation similar to the Arrhenius equation. The results are discussed in terms of transition state theory.     

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.