Coherent third-order spectroscopic probes of molecular chirality

Equivalence of quantum and classical third order response for weakly anharmonic coupled oscillators

Two-dimensional (2D) infrared (IR) spectra are commonly interpreted using a quantum diagrammatic expansion that describes the changes to the density matrix of quantum systems in response to light-matter interactions. Although classical response functions (based on Newtonian dynamics) have shown promise in computational 2D IR modeling studies, a simple diagrammatic description has so far been lacking. Recently, we introduced a diagrammatic representation for the 2D IR response functions of a single, weakly anharmonic oscillator and showed that the classical and quantum 2D IR response functions for this system are identical. Here, we extend this result to systems with an arbitrary number of bilinearly coupled, weakly anharmonic oscillators. As in the single-oscillator case, quantum and classical response functions are found to be identical in the weakly anharmonic limit or, in experimental terms, when the anharmonicity is small relative to the optical linewidth. The final form of the weakly anharmonic response function is surprisingly simple and offers potential computational advantages for application to large, multi-oscillator systems.

Optimized interferometric setup for chiral and achiral ultrafast IR spectroscopy

We report an actively stabilized interferometer-based set-up for the detection of vibrational circular dichroism (VCD) and optical rotatory dispersion (VORD) with femtosecond laser pulses. Our approach combines and improves elements of several previous measurement strategies, including signal amplification in a crossed polarizer configuration, precise control and modulation of polarization, phase stability, tight focusing, broad-band detection and spectral interferometry. Their importance for static and transient measurements is motivated by a signal analysis based on Jones matrices and response theory. Only depending on the pump-beam polarization, the set-up can selectively detect transient VCD and VORD or transient linear birefringence (LB) and linear dichroism (LD), which usually constitute the dominant artifacts in the chiral measurements. For illustration we present transient LB and LD data of an achiral Rhenium carbonyl complex, detected simultaneously by spectral interferometry, and we analyze residual background signals in the experimental configuration for transient chiral spectroscopy.

Two-dimensional x-ray correlation spectroscopy of remote core states

Nonlinear all-X-ray signals that involve large core-atom separation compared to the X-ray wavelengths may not be described by the dipole approximation since they contain additional phase factors. Expressions for the rotationally averaged 2D X-ray photon echo signals from randomly oriented systems that take this position-dependent phase into account for arbitrary ratio between the core separation and the resonant wavelength are presented. Application is made to the Se K-edge of a selenium dipeptide system.

Monitoring the folding of Trp-cage peptide by two-dimensional infrared (2DIR) spectroscopy

Protein folding is one of the most fundamental problems in modern molecular biology. Uncovering the detailed folding mechanism requires methods that can monitor the structures at high temporal and spatial resolution. Two-dimensional infrared (2DIR) spectroscopy is a new tool for studying protein structures and dynamics with high time resolution. Using atomistic molecular dynamics simulations, we illustrate the folding process of Trp-cage along the dominant pathway on the free energy landscape by analyzing nonchiral and chiral coherent 2DIR spectra along the pathway. Isotope labeling is used to reveal residue-specific information. We show that the high resolution structural sensitivity of 2DIR can differentiate the ensemble evolution of protein and thus provides a microscopic picture of the folding process.

Multidimensional attosecond resonant X-ray spectroscopy of molecules: lessons from the optical regime

New free-electron laser and high-harmonic generation X-ray light sources are capable of supplying pulses short and intense enough to perform resonant nonlinear time-resolved experiments in molecules. Valence-electron motions can be triggered impulsively by core excitations and monitored with high temporal and spatial resolution. We discuss possible experiments that employ attosecond X-ray pulses to probe the quantum coherence and correlations of valence electrons and holes, rather than the charge density alone, building on the analogy with existing studies of vibrational motions using femtosecond techniques in the visible regime.

Coherent control protocol for separating energy-transfer pathways in photosynthetic complexes by chiral multidimensional signals

Adaptive optimizations performed using a genetic algorithm are employed to construct optimal laser pulse configurations that separate spectroscopic features associated with the two main energy-transfer pathways in the third-order nonlinear optical response simulated for the Fenna-Matthews-Olson (FMO) photosynthetic complex from the green sulfur bacterium Chlorobium tepidum. Superpositions of chirality-induced tensor components in both collinear and noncollinear pulse configurations are analyzed. The optimal signals obtained by manipulating the ratios of various 2D spectral peaks reveal detailed information about the excitation dynamics.

Probing many-particle correlations in semiconductor quantum wells using double-quantum-coherence signals

Multidimensional analysis of coherent signals is commonly used in nuclear magnetic resonance to study correlations among spins. These techniques were recently extended to the femtosecond regime and applied to chemical, biological and semiconductor systems. In this work, we apply a two-dimensional correlation spectroscopy technique which employs double-quantum-coherence to investigate many-body effects in a semiconductor quantum well. The signal is detected along the direction k(1)+ k(2)- k(3), where k(1), k(2) and k(3) are the pulse wave vectors in chronological order. We show that this signal is particularly sensitive to many-body correlations which are missed by time-dependent Hartree-Fock approximation. The correlation energy of two-exciton can be probed with a very high resolution arising from a two-dimensional correlation spectrum, where two-exciton couplings spread the cross peaks along both axes of the 2D spectrum to create a characteristic highly resolved pattern. This level of detail is not available from conventional one-dimensional four-wave mixing or other two-dimensional correlation spectroscopy signals such as the photo echo (-k(1)+ k(2)+ k(3)).

Polarization shaping in the mid-IR and polarization-based balanced heterodyne detection with application to 2D IR spectroscopy

We demonstrate amplitude, phase and polarization shaping of femtosecond mid-IR pulses using a germanium acousto-optical modulator by independently shaping the frequency-dependent amplitudes and phases of two orthogonally polarized pulses which are then collinearly overlapped using a wire-grid polarizer. We use a feedback loop to set and stabilize the relative phase of the orthogonal pulses. We have also used a wire-grid polarizer to implement polarization-based balanced heterodyne detection for improved signal-to-noise of 2D IR spectra collected in a pump-probe geometry. Applications include coherent control of molecular vibrations and improvements in multidimensional IR spectroscopy.

Coherent multidimensional vibrational spectroscopy of biomolecules: concepts, simulations, and challenges

The response of complex molecules to sequences of femtosecond infrared pulses provides a unique window into their structure, dynamics, and fluctuating environments. Herein we survey the basic principles of modern two-dimensional infrared (2DIR) spectroscopy, which analogous to those of multidimensional NMR spectroscopy. The perturbative approach for computing the nonlinear optical response of coupled localized chromophores is introduced and applied to the amide backbone transitions of proteins, liquid water, membrane lipids, and amyloid fibrils. The signals are analyzed using classical molecular dynamics simulations combined with an effective fluctuating Hamiltonian for coupled localized anharmonic vibrations whose dependence on the local electrostatic environment is parameterized by an ab initio map. Several simulation methods, (cumulant expansion of Gaussian fluctuation, quasiparticle scattering, the stochastic Liouville equations, direct numerical propagation) are surveyed. Chirality-induced techniques which dramatically enhance the resolution are demonstrated. Signatures of conformational and hydrogen-bonding fluctuations, protein folding, and chemical-exchange processes are discussed.

Double-quantum resonances and exciton-scattering in coherent 2D spectroscopy of photosynthetic complexes

A simulation study demonstrates how the nonlinear optical response of the Fenna-Matthews-Olson photosynthetic light-harvesting complex may be explored by a sequence of laser pulses specifically designed to probe the correlated dynamics of double excitations. Cross peaks in the 2D correlation plots of the spectra reveal projections of the double-exciton wavefunctions onto a basis of direct products of single excitons. An alternative physical interpretation of these signals in terms of quasiparticle scattering is developed.

Selection rules and symmetry relations for four-wave mixing measurements of uniaxial assemblies

Uniaxial systems represent the next lowest symmetry below isotropic and are ubiquitous. The objective of the present work is to present a systematic foundation for interpreting polarization-dependent four-wave mixing measurements of oriented and aligned assemblies. Orientational averages connecting the molecular frame to the macroscopic frame in uniaxial assemblies were derived for several common molecular symmetry groups for coherent anti-Stokes Raman spectroscopy (CARS) measurements, coherent anti-Stokes two-photon spectroscopy (CATS) probing electronic transitions, resonant two-photon absorption (2PA), and traditional Raman measurements. First, the complete set of orientational averages connecting the molecular and macroscopic frames was compiled for the most general case of C1 molecular symmetry. Then, the orientational averages of a select few commonly occurring molecular symmetry groups (Cs, C2, C2v, and C3v) were explored in greater detail to illustrate the approach and to facilitate the interpretation of routine experimental measurements. One outcome of this analysis is the prediction of efficient electric dipole-allowed chiral-specific four-wave mixing in uniaxially oriented media.

Chirality-induced signals in coherent multidimensional spectroscopy of excitons

The nonlocal second- and third-order susceptibilities of an isotropic ensemble of aggregates are calculated by solving the nonlinear exciton equations which map the system into coupled anharmonic oscillators. Both electric and magnetic contributions are included using the minimal-coupling Hamiltonian. The various tensor components are evaluated to first order in the optical wave vector k. Additional structural information about the interchromophore distances, which is not accessible through zeroth-order contributions (the dipole approximation), is contained to the first order in k. New resonant second- and third-order signals predicted for chiral molecules provide multidimensional extensions of circular dichroism spectroscopy. Numerical simulations demonstrate the sensitivity of third-order signals to the secondary structural motiffs of peptides.

Coherent control of pump-probe signals of helical structures by adaptive pulse polarizations

The simplification of the pump-probe spectrum of excitons by pure-phase-polarization pulse shaping is investigated by a simulation study. The state of light is manipulated by varying the phases of two perpendicular polarization components of the pump, holding its total spectral and temporal intensity profiles fixed. Genetic and iterative Fourier transform algorithms are used to search for pulse phase functions that optimize the ratio of the signal at two frequencies. New features are extracted from the congested pump-probe spectrum of a helical pentamer by selecting a combination of Liouville space pathways. Tensor components which dominate the optimized spectra are identified.

Quadrupole contribution to the third-order optical activity spectroscopy

Time-resolved nonlinear optical activity measurement spectroscopy can be a useful tool for studying biomolecular and chemical reaction dynamics of chiral molecules. Only recently, the two-dimensional (2D) circularly polarized photon echo (CP-PE) spectroscopy of polypeptides and a photosynthetic light-harvesting complex were discussed, where the beam configuration was specifically controlled in such a way to eliminate the quadrupole contribution to the CP-PE signal. In this paper, we generalize the CP-PE spectroscopy by including the transition quadrupole contributions from peptide amide I vibrational transition and chlorophyll electronic transition. By using a density functional theory calculation method, the corresponding amide I vibrational and chlorophyll Q(y) electronic transition quadrupole tensor elements are determined. Amplitude of nonlinear optical transition pathway involving a quadrupole transition is found to be comparable to those of magnetic dipole terms for two different cases considered, i.e., dipeptides and photosynthetic antenna complex. However, due to the rotational averaging factors, the overall quadrupole contribution is an order of magnitude smaller than the magnetic dipole contribution. This suggests that the conventional 2D photon echo method and experimental scheme can be directly used to measure the 2D CP-PE signal from proteins and molecular complexes and that the 2D CP-PE signal is mainly dictated by the magnetic dipole contribution.

Simulation protocols for coherent femtosecond vibrational spectra of peptides

Two algorithms for simulating the response of peptides to sequences of IR pulses are developed and applied to N-methyl acetamide (NMA) and a 17 residue alpha-helical peptide (YKKKH17) in D(2)O. A fluctuating vibrational-exciton Hamiltonian for the amide I mode is constructed from molecular dynamics trajectories. Coupling with the environment is described using a density functional theory electrostatic map. The cumulant expansion of Gaussian fluctuation incorporates motional narrowing due to fast frequency fluctuations and is adequate for NMA and for isotopically labeled bands in large peptides. Real-space truncation of the scattering matrix of the nonlinear exciton equations significantly reduces the computational cost, making it particularly attractive for slow fluctuations in large globular proteins.

Intermittent search process and teleportation

The authors study an intermittent search process combining diffusion and "teleportation" phases in a d-dimensional spherical continuous system and in a regular lattice. The searcher alternates diffusive phases, during which targets can be discovered, and fast phases (teleportation) which randomly relocate the searcher, but do not allow for target detection. The authors show that this alternation can be favorable for minimizing the time of first discovery, and that this time can be optimized by a convenient choice of the mean waiting times of each motion phase. The optimal search strategy is explicitly derived in the continuous case and in the lattice case. Arguments are given to show that much more general intermittent motions do provide optimal search strategies in d dimensions. These results can be useful in the context of heterogeneous catalysis or in various biological examples of transport through membrane pores.

Two-Dimensional Circularly Polarized IR Photon Echo Spectroscopy of Polypeptides: Four-Wave-Mixing Optical Activity Measurement

A coherent two-dimensional (2D) optical spectroscopy utilizing circularly polarized (CP) beams, which was shown to be useful in studying molecular chirality in condensed phases, was theoretically proposed recently [Cho et al. J. Chem. Phys. 2003, 119, 7003]. A photon echo (PE) version of 2D optical activity spectroscopy is discussed in this paper. Considering various dipeptide and polypeptide systems, where the amide I local modes constitute the set of basis modes used to describe exciton and biexciton states as linear combinations of those basis modes, we present numerically simulated 2D circularly polarized IR PE spectra. It is shown that this novel spectroscopic method can provide additional information on the angles between the transition magnetic dipole and the transition electric dipole of two different vibrationally excited states, which are highly sensitive to the 3D structure and chirality of a given polypeptide. Also, a hierarchical relation of IR absorption, vibrational circular dichroism, 2D IR PE, and 2D CP-IR PE is discussed to show advantages of 2D optical activity spectroscopy in general.

Synchronized and configurable source of electrical pulses for x-ray pump-probe experiments

A method is described for the generation of software tunable patterns of nanosecond electrical pulses. The bipolar, high repetition rate (up to 250 MHz), fast rise time (<30 ps), square pulses are suitable for applications such as the excitation sequence in dynamic pump-probe experiments. Synchronization with the time structure of a synchrotron facility is possible as well as fine control of the relative delay in steps of 10 ps. The pulse generator described here is used to excite magnetic nanostructures with current pulses. Having an excitation system which can match the high repetition rate of a synchrotron allows for utilization of the full x-ray flux and is needed in experiments which require a large photon flux. The fast rise times allow for picosecond time resolution in pump-probe experiments. All pulse pattern parameters are configurable by software.

Manipulating multidimensional electronic spectra of excitons by polarization pulse shaping

A simulation study demonstrates how coherent control, combined with adaptive polarization pulse shaping and a genetic algorithm, may be used to simplify femtosecond coherent nonlinear optical signals of excitons. Cross peaks are amplified and resolved, and diagonal peaks are suppressed in the heterodyne-detected two-pulse echo signal from the Soret band of a porphyrin dimer coupled to a Brownian oscillator bath. Various optimization strategies involving the spectral, temporal, and polarization profiles of the second pulse are compared.

Appearance of intramolecular high-frequency vibrations in two-dimensional, time-integrated three-pulse photon echo data

An alternative experimental outline to measure homodyne detected three-pulse photon-echo data is presented. The novel experimental approach allowing for online monitoring and correction of experimental timing and stability is discussed in detail using the paradigm system of Nile blue in alcohol solution. It is shown that excellent signal-to-noise ratios together with high reproducibility of the data can be routinely achieved. We report in detail on the appearance of high-frequency intramolecular vibrations in the two-dimensional three-pulse photon-echo data and suggest that besides the conventionally discussed three-pulse photon-echo peak-shift the width of the integrated echo signal as a function of population time contains identical and easily accessible information on high-frequency intramolecular vibrations. A comparison of experimental data with theoretical modeling is performed showing that the observed echo-width oscillations are in line with predictions of the Brownian oscillator model.

Coherent control of cross-peaks in chirality-induced two-dimensional optical signals of excitons

Polarization tuning of the interference of chirality-induced tensor components is used to simulate the suppression of diagonal peaks and amplification of cross peaks in femtosecond two-dimensional photon echo signals of excitons in a chiral porphyrin dimer. Superpositions of various tensor components which generate the optimized signals are constructed using a genetic learning algorithm. Exciton couplings and bath correlations may be extracted from these highly resolved signals.

Water flow in poly(N-isopropylacrylamide) gels

Friction between a polymer network of poly(N-isopropylacrylamide) gels and solvent water was investigated. The gel was mechanically constrained in a glass capillary at gelation, and hydrostatic pressure was directly applied to the cross section of the cylinder. The temperature dependence of the flow velocity was extensively measured in the vicinity of the transition temperature for gels with different lengths, l(0), at gelation. As the temperature increased, the friction slightly decreased at the transition point and increased rapidly in the collapsed phase. Although the flow velocity depended on l(0), the friction in the vicinity of the transition point was well scaled by l(0) based on the Hagen-Poiseuille equation for the flux of water flow in a capillary. The results suggested that the assumption that the gel is a bundle of microcapillaries was applicable to the water flow through the hydrogel, which was largely deformed not only by the pressure applied to the solvent but also by the shrinking force caused by the temperature increment. Macroscopic deformation did not affect the friction between the three-dimensional polymer network and water.

Anisotropy of photofragment recoil as a function of dissociation lifetime, excitation frequency, rotational level, and rotational constant

Quantum mechanical calculations of photofragment angular distributions have been performed as a function of the frequency of excitation, the lifetime of the dissociative state, the rotational level, and the rotational constant. In the limit of high J values and white, incoherent excitation, the general results are found to agree exactly with both those of Mukamel and Jortner [J. Chem. Phys. 61, 5348 (1974)] and those of Jonah [J. Chem. Phys. 55, 1915 (1971)]. Example calculations describe how the anisotropy is dependent on the degree of broadening, the rotational constant, the initial rotational level, and the frequency of excitation. Applications are also made to interpret experimental results on the photodissociation of ClO via the 11-0, 10-0, and 6-0 bands of the A 2Pi3/2 -X 2Pi3/2 transition and on the photodissociation of O2 via the 0-0 band of the E 3Sigmau- -X 3Sigmag- transition.

Characterization of the solvation dynamics of an ionic liquid via molecular dynamics simulation

The solvation dynamics of ionic liquids have been the subject of intense experimental study but remain poorly understood. We present the results of molecular dynamics simulations of the solvation dynamics of the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate in response to photoexcitation of the fluorescent dye coumarin-153. We reproduce the time-resolved fluorescence Stokes shift using linear response theory, then use novel statistical techniques to analyze cation and anion contributions to the signal. We find that the solvation dynamics are dominated by collective ionic motion and characterize the time scale for various features of the collective response. Further, we use the Steele analysis [Mol. Phys. 61, 1031 (1987)] to characterize the contributions to the observed Stokes shift made by translational and rovibrational degrees of freedom. Our results indicate that in contrast to molecular liquids, the rovibrational response is trivial and the observed fluorescence response arises almost entirely from ionic translation. Our results resolve previously open questions in the literature about the nature of the rapid dynamics in room-temperature ionic liquids and offer insight into the physical principles governing ionic liquid behavior on longer time scales.

Homogeneous nucleation rate measurements of 1-propanol in helium: The effect of carrier gas pressure

Kinetics of homogeneous nucleation in supersaturated vapor of 1-propanol was studied using an upward thermal diffusion cloud chamber. Helium was used as a noncondensable carrier gas and the influence of its pressure on observed nucleation rates was investigated. The isothermal nucleation rates were determined by a photographic method that is independent on any nucleation theory. In this method, the trajectories of growing droplets are recorded using a charge coupled device camera and the distribution of local nucleation rates is determined by image analysis. The nucleation rate measurements of 1-propanol were carried out at four isotherms 260, 270, 280, and 290 K. In addition, the pressure dependence was investigated on the isotherms 290 K (50, 120, and 180 kPa) and 280 K (50 and 120 kPa). The isotherm 270 K was measured at 25 kPa and the isotherm 260 K at 20 kPa. The experiments confirm the earlier observations from several thermal diffusion chamber investigations that the homogeneous nucleation rate of 1-propanol tends to increase with decreasing total pressure in the chamber. In order to reduce the possibility that the observed phenomenon is an experimental artifact, connected with the generally used one-dimensional description of transfer processes in the chamber, a recently developed two-dimensional model of coupled heat, mass, and momentum transfer inside the chamber was used and results of both models were compared. It can be concluded that the implementation of the two-dimensional model does not explain the observed effect. Furthermore the obtained results were compared both to the predictions of the classical theory and to the results of other investigators using different experimental devices. Plotting the experimental data on the so-called Hale plot shows that our data seem to be consistent both internally and also with the data of others. Using the nucleation theorem the critical cluster sizes were obtained from the slopes of the individual isotherms and compared with the Kelvin prediction. The influence of total pressure on the observed isothermal nucleation rate was studied in another experiment, where not only temperature but also supersaturation was kept constant as the total pressure was changed. It was shown that the dependence of the nucleation rate on pressure gets stronger as pressure decreases.

Nonlinear optical spectroscopy of chiral molecules

We review nonlinear optical processes that are specific to chiral molecules in solution and on surfaces. In contrast to conventional natural optical activity phenomena, which depend linearly on the electric field strength of the optical field, we discuss how optical processes that are nonlinear (quadratic, cubic, and quartic) functions of the electromagnetic field strength may probe optically active centers and chiral vibrations. We show that nonlinear techniques open entirely new ways of exploring chirality in chemical and biological systems: The cubic processes give rise to nonlinear circular dichroism and nonlinear optical rotation and make it possible to observe dynamic chiral processes at ultrafast time scales. The quadratic second-harmonic and sum-frequency-generation phenomena and the quartic processes may arise entirely in the electric-dipole approximation and do not require the use of circularly polarized light to detect chirality. They provide surface selectivity and their observables can be relatively much larger than in linear optical activity. These processes also give rise to the generation of light at a new color, and in liquids this frequency conversion only occurs if the solution is optically active. We survey recent chiral nonlinear optical experiments and give examples of their application to problems of biophysical interest.