Higher-order optical correlation spectroscopy in liquids

Modeling, calculating, and analyzing multidimensional vibrational spectroscopies

Spectral line shapes in a condensed phase contain information from various dynamic processes that modulate the transition energy, such as microscopic dynamics, inter- and intramolecular couplings, and solvent dynamics. Because nonlinear response functions are sensitive to the complex dynamics of chemical processes, multidimensional vibrational spectroscopies can separate these processes. In multidimensional vibrational spectroscopy, the nonlinear response functions of a molecular dipole or polarizability are measured using ultrashort pulses to monitor inter- and intramolecular vibrational motions. Because a complex profile of such signals depends on the many dynamic and structural aspects of a molecular system, researchers would like to have a theoretical understanding of these phenomena. In this Account, we explore and describe the roles of different physical phenomena that arise from the peculiarities of the system-bath coupling in multidimensional spectra. We also present simple analytical expressions for a weakly coupled multimode Brownian system, which we use to analyze the results obtained by the experiments and simulations. To calculate the nonlinear optical response, researchers commonly use a particular form of a system Hamiltonian fit to the experimental results. The optical responses of molecular vibrational motions have been studied in either an oscillator model or a vibration energy state model. In principle, both models should give the same results as long as the energy states are chosen to be the eigenstates of the oscillator model. The energy state model can provide a simple description of nonlinear optical processes because the diagrammatic Liouville space theory that developed in the electronically resonant spectroscopies can easily handle three or four energy states involved in high-frequency vibrations. However, the energy state model breaks down if we include the thermal excitation and relaxation processes in the dynamics to put the system in a thermal equilibrium state. The roles of these excitation and relaxation processes are different and complicated compared with those in the resonant spectroscopy. Observing the effects of such thermal processes is more intuitive with the oscillator model because the bath modes, which cause the fluctuation and dissipation processes, are also described in the coordinate space. This coordinate space system-bath approach complements a realistic full molecular dynamics simulation approach. By comparing the calculated 2D spectra from the coordinate space model and the energy state model, we can examine the role of thermal processes and anharmonic mode-mode couplings in the energy state model. For this purpose, we employed the Brownian oscillator model with the nonlinear system-bath interaction. Using the hierarchy formalism, we could precisely calculate multidimensional spectra for a single and multimode anharmonic system for inter- and intramolecular vibrational modes.

Femtosecond Stimulated Raman Spectroscopy

Femtosecond stimulated Raman spectroscopy (FSRS) is a new ultrafast spectroscopic technique that provides vibrational structural information with high temporal (50-fs) and spectral (10-cm(1)) resolution. As a result of these unique capabilities, FSRS studies of chemical and biochemical reaction dynamics are expected to grow rapidly, giving previously unattainable insight into the structural dynamics of reactively evolving systems with atomic spatial and femtosecond temporal resolution. This review discusses the experimental and theoretical concepts behind FSRS, with an emphasis on the origins of its unique temporal and spectral capabilities. We illustrate these capabilities with vibrational studies of ultrafast electronic dynamics, as well as the direct structural observation of nonstationary vibrational wave-packet motion in small molecules and in complex biochemical reaction dynamics.

Hybrid molecular dynamics-quantum mechanics simulations of solute spectral properties in the condensed phase: Evaluation of simulation parameters

We have explored the impact of a number of basic simulation parameters on the results of a recently developed hybrid molecular dynamics-quantum mechanics (MD-QM) method (Mercer et al., J Phys Chem B 1999, 103, 7720). The method utilizes MD simulations to explore the ground-state configuration space of the system and QM evaluation of those structures to yield the time-dependent electronic transition energy, which is transformed into the optical line-broadening function using the second-order cumulant expansion. Both linear and nonlinear optical spectra can then be generated for comparison to experiment. The dependence of the resulting spectra on the length of the MD trajectory, the QM sampling rate, and the QM model chemistry have all been examined. In particular, for the system of oxazine-4 in methanol studied here, at least 20 ps of MD trajectory are needed for qualitative convergence of linear spectral properties, and >100 ps is needed for quantitative convergence. Surprisingly, little difference is found between the 3-21G and 6-31G(d) basis sets, and the CIS and TD-B3LYP methods yield remarkably similar spectra. The semiempirical INDO/s method yields the most accurate results, reproducing the experimental Stokes shift to within 5% and the FWHM to within 20%. Nonlinear 3-pulse photon echo peak shift (3PEPS) decays have also been simulated. Decays are generally poorly reproduced, though the initial peak shift which depends on the overall coupling of motions to the solute transition energy is within 15% of experiment for all model chemistries other than those using the STO-3G basis.

Fifth-order contributions to ultrafast spectrally resolved vibrational echoes: heme-CO proteins

The fifth order contributions to the signals of ultrafast infrared spectrally resolved stimulated vibrational echoes at high intensities have been investigated in carbonmonoxy heme proteins. High intensities are often required to obtain good data. Intensity dependent measurements are presented on hemoglobin-CO (Hb-CO) and a mutant of myoglobin, H64V-CO. The spectrally resolved vibrational echoes demonstrate that fifth order effects arise at both the 1-0 and the 2-1 emission frequencies of the stretching mode of the CO chromophore bound at the active site of heme proteins. Unlike one-dimensional experiments, in which the signal is integrated over all emission frequencies, spectrally resolving the signal shows that the fifth order contributions have a much more pronounced influence on the 2-1 transition than on the 1-0 transition. By spectrally isolating the 1-0 transition, the influence of fifth order contributions to vibrational echo data can be substantially reduced. Analysis of fifth order Feynman diagrams that contribute in the vibrational echo phase-matched direction demonstrates the reason for the greater influence of fifth order processes on the 1-2 transition, and that the fifth order contributions are heterodyne amplified by the third order signal. Finally, it is shown that the anharmonic oscillations in vibrational echo data of Hb-CO that previous work had attributed strictly to fifth order effects arise even without fifth order contributions.

Multidimensional infrared spectroscopy of the N-H bond motions in formamide

The heterodyned two-dimensional (2D) IR spectra and equilibrium dynamics of the N-H stretching motion of DCONHD in deuterated formamide, DCOND(2), were studied with 80 fs pulses at 3 microm. The time evolution of the heterodyned 2D IR spectra, pump-probe spectra, and photon echo peak shift demonstrate that interstate dynamics is occurring by relaxation of the original N-H excitation. The N-H vibrational frequency correlation function can be expressed as a sum of three exponentials with correlation times 0.24 ps, 0.8 ps, and 11 ps. The intermediate component is attributed to motions of the N-Hcdots, three dots, centeredO unit involving only slight angular variations of the N-H bond. The slow component is attributed to the structure breaking and making. The anisotropy decay confirmed that the significant angular N-H bond motion occurs on the 11 ps time scale. The fast component, which is the least well determined, might correspond to the modulation of the H-bond distance without angular motion. The correlation coefficient between the pumped and relaxed state distributions was +0.51, implying that the excited state phase memory is only slightly diminished by the relaxation of the N-H excitation. The relaxed modes are concluded to be local to the driven N-H mode.

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.

Vibrational dynamics of N–H, C–D, and CO modes in formamide

By means of heterodyned two-dimensional IR photon echo experiments on liquid formamide and isotopomers the vibrational frequency dynamics of the N-H stretch mode, the C-D mode, and the C=O mode were obtained. In each case the vibrational frequency correlation function is fitted to three exponentials representing ultrafast (few femtoseconds), intermediate (hundreds of femtoseconds), and slow (many picoseconds) correlation times. In the case of N-H there is a significant underdamped contribution to the correlation decay that was not seen in previous experiments and is attributed to hydrogen-bond librational modes. This underdamped motion is not seen in the C-D or C=O correlation functions. The motions probed by the C-D bond are generally faster than those seen by N-H and C=O, indicating that the environment of C-D interchanges more rapidly, consistent with a weaker C-D...O=C bond. The correlation decays of N-H and C=O are similar, consistent with both being involved in strong H bonding.

Ultrafast molecular dynamics of liquid aromatic molecules and the mixtures with CCl4

The ultrafast molecular dynamics of liquid aromatic molecules, benzene, toluene, ethylbenzene, cumene, and 1,3-diphenylpropane, and the mixtures with CCl(4) have been investigated by means of femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy. The picosecond Kerr transients of benzene, toluene, ethylbenzene, and cumene and the mixtures with CCl(4) show a biexponential feature. 1,3-Diphenylpropane and the mixtures with CCl(4) show triexponential picosecond Kerr transients. The slow relaxation time constants of the aromatic molecules and the mixtures with CCl(4) are qualitatively described by the Stoke-Einstein-Debye hydrodynamic model. The ultrafast dynamics have been discussed based on the Kerr spectra in the frequency range of 0-800 cm(-1) obtained by the Fourier transform analysis of the Kerr transients. The line shapes of the low-frequency intermolecular spectra located at 0-180 cm(-1) frequency range have been analyzed by two Brownian oscillators ( approximately 11 cm(-1) and approximately 45 cm(-1) peaks) and an antisymmetric Gaussian function ( approximately 65 cm(-1) peak). The spectrum shape of 1,3-diphenylpropane is quite different from the spectrum shapes of the other aromatic molecules for the low magnitude of the low-frequency mode of 1,3-diphenylpropane and/or an intramolecular vibration. Although the concentration dependences of the low- and intermediate-frequency intermolecular modes (Brownian oscillators) do not show a significant trend, the width of high-frequency intermolecular mode (antisymmetric Gaussian) becomes narrower with the higher CCl(4) concentration for all the aromatics mixtures with CCl(4). The result indicates that the inhomogeneity of the intermolecular vibrational mode in aromatics/CCl(4) mixtures is decreasing with the lower concentration of aromatics. The intramolecular vibrational modes of the aromatic molecules observed in the Kerr spectra are also shown with the calculation results based on the density functional theory.

Two-dimensional femtosecond spectroscopy

The simplest two-dimensional (2D) spectra show how excitation with one (variable) frequency affects the spectrum at all other frequencies, thus revealing the molecular connections between transitions. Femtosecond 2D Fourier transform (2D FT) spectra are more flexible and share some of the remarkable properties of their conceptual parent, 2D FT nuclear magnetic resonance. When 2D FT spectra are experimentally separated into real absorptive and imaginary refractive parts, the time resolution and frequency resolution can both reach the uncertainty limit set for each resonance by the sample itself. Coherent four-level contributions to the signal provide new molecular phase information, such as relative signs of transition dipoles. The nonlinear response can be picked apart by selecting a single coherence pathway (e.g., specifying the relative signs of energy level difference frequencies during different time intervals as in the photon echo). Because molecules are frozen on the femtosecond timescale, femtosecond 2D FT experiments can separate a distribution of instantaneous molecular environments and intramolecular geometries as inhomogeneous broadening. This review provides an introduction to two-dimensional Fourier transform experiments exploiting second- and third-order vibrational and electronic nonlinearities.