A simple model for the calculation of nonlinear optical response functions and femtosecond time-resolved spectra

Two-dimensional fluorescence excitation spectroscopy: A novel technique for monitoring excited-state photophysics of molecular species with high time and frequency resolution

We propose a novel UV/Vis femtosecond spectroscopic technique, two-dimensional fluorescence-excitation (2D-FLEX) spectroscopy, which combines spectral resolution during the excitation process with exclusive monitoring of the excited-state system dynamics at high time and frequency resolution. We discuss the experimental feasibility and realizability of 2D-FLEX, develop the necessary theoretical framework, and demonstrate the high information content of this technique by simulating the 2D-FLEX spectra of a model four-level system and the Fenna-Matthews-Olson antenna complex. We show that the evolution of 2D-FLEX spectra with population time directly monitors energy transfer dynamics and can thus yield direct qualitative insight into the investigated system. This makes 2D-FLEX a highly efficient instrument for real-time monitoring of photophysical processes in polyatomic molecules and molecular aggregates.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

RhoDyn: A ρ-TD-RASCI Framework to Study Ultrafast Electron Dynamics in Molecules

This article presents the program module RhoDyn as part of the OpenMOLCAS project intended to study ultrafast electron dynamics within the density-matrix-based time-dependent restricted active space configuration interaction framework (ρ-TD-RASCI). The formalism allows for the treatment of spin-orbit coupling effects, accounts for nuclear vibrations in the form of a vibrational heat bath, and naturally incorporates (auto)ionization effects. Apart from describing the theory behind and the program workflow, the paper also contains examples of its application to the simulations of the linear L_2,3 absorption spectra of a titanium complex, high harmonic generation in the hydrogen molecule, ultrafast charge migration in benzene and iodoacetylene, and spin-flip dynamics in the core excited states of iron complexes.

Simulating Fluorescence-Detected Two-Dimensional Electronic Spectroscopy of Multichromophoric Systems

We present a theory for modeling fluorescence-detected two-dimensional electronic spectroscopy of multichromophoric systems. The theory is tested by comparison of the predicted spectra of the light-harvesting complex LH2 with experimental data. A qualitative explanation of the strong cross-peaks as compared to conventional two-dimensional electronic spectra is given. The strong cross-peaks are attributed to the clean ground-state signal that is revealed when the annihilation of exciton pairs created on the same LH2 complex cancels oppositely signed signals from the doubly excited state. This annihilation process occurs much faster than the nonradiative relaxation. Furthermore, the line shape difference is attributed to slow dynamics, exciton delocalization within the bands, and intraband exciton-exciton annihilation. This is in line with existing theories presented for model systems. We further propose the use of time-resolved fluorescence-detected two-dimensional spectroscopy to study state-resolved exciton-exciton annihilation.

Mapping the Evolution of Spatial Exciton Coherence through Time-Resolved Fluorescence

Quantum coherence is expected to have a positive effect on the transfer efficiency of excitation energy through photosynthetic aggregates and conjugated polymers, but its significance to the functioning of these molecular assemblies remains largely unknown. We propose a new experimental means to monitor the coherence between distant molecular sites on a time scale relevant to energy transfer. Through numerical calculations, we demonstrate that the range of such spatial coherence continually scales as the 0-0 to 0-1 vibronic peak ratio in time-resolved fluorescence spectroscopy. As such, this observable allows one to monitor the coherent evolution of an excited state, displaying the large coherence length following optical excitation, and the subsequent dephasing over time.

Exact quantum master equation for a molecular aggregate coupled to a harmonic bath

We consider a molecular aggregate consisting of N identical monomers. Each monomer comprises two electronic levels and a single harmonic mode. The monomers interact with each other via dipole-dipole forces. The monomer vibrational modes are bilinearly coupled to a bath of harmonic oscillators. This is a prototypical model for the description of coherent exciton transport, from quantum dots to photosynthetic antennae. We derive an exact quantum master equation for such systems. Computationally, the master equation may be useful for the testing of various approximations employed in theories of quantum transport. Physically, it offers a plausible explanation of the origins of long-lived coherent optical responses of molecular aggregates in dissipative environments.

Efficient calculation of time- and frequency-resolved four-wave-mixing signals

"Four-wave-mixing" is the generic name for a family of nonlinear electronic and vibrational spectroscopies. These techniques are widely used to explore dissipation, dephasing, solvation, and interstate coupling mechanisms in various material systems. Four-wave-mixing spectroscopy needs a firm theoretical support, because it delivers information on material systems indirectly, through certain transients, which are measured as functions of carrier frequencies, durations, and relative time delays of the laser pulses. The observed transients are uniquely determined by the three-pulse-induced third-order polarization. There exist two conceptually different approaches to the calculation of the nonlinear polarization. In the standard perturbative approach to nonlinear spectroscopy, the third-order polarization is expressed in terms of the nonlinear response functions. As the material systems become more complex, the evaluation of the response functions becomes cumbersome and the calculation of the signals necessitates a number of approximations. Herein, we review alternative methods for the calculation of four-wave-mixing signals, in which the relevant laser pulses are incorporated into the system Hamiltonian and the driven system dynamics is simulated numerically exactly. The emphasis is on the recently developed equation-of-motion phase-matching approach (EOM-PMA), which allows us to calculate the three-pulse-induced third-order polarization in any phase-matching direction by performing three (with the rotating wave approximation) or seven (without the rotating wave approximation) independent propagations of the density matrix. The EOM-PMA is limited to weak laser fields (its domain of validity is equivalent to the approach based on the third-order response functions) but allows for arbitrary pulse durations and automatically accounts for pulse-overlap effects. As an illustration, we apply the EOM-PMA to the calculation of optical three-pulse photon-echo two-dimensional (2D) signals for a generic model system, which represents a characteristic photophysical dynamics of large molecules or chromophores in condensed phases. The EOM-PMA is easy to implement and can straightforwardly be incorporated into any computational scheme, which provides the time-dependent density matrix or wave function of the material system of interest. In particular, EOM-PMA-based computer codes can efficiently be implemented on parallel computers. The EOM-PMA facilitates considerably the computation of four-wave-mixing signals and 2D spectra, in both vibrational and electronic spectroscopy. The EOM-PMA can be extended to higher order optical responses, e.g., heterodyned 3D IR, transient 2D IR, and other six-wave-mixing techniques.

A nonperturbative calculation of nonlinear spectroscopic signals in liquid solution

Nonlinear spectroscopic signals in liquid solution were calculated without treating the field-matter interaction in a perturbative manner. The calculation is based on the assumption that the intermolecular degrees of freedom can be treated classically, while the time evolution of the electronic state is treated quantum mechanically. The calculated overall electronic polarization is then resolved into its directional components via the method of Seidner et al. [J. Chem. Phys. 103, 3998 (1995)]. It is shown that the time dependence of the directional components is independent of laser intensity in the impulsive pulse regime, which allows for flexibility in choosing the procedure for calculating optical response functions. The utility and robustness of the nonperturbative procedure is demonstrated in the case of a two-state chromophore solvated in a monoatomic liquid, by calculating nonlinear time-domain signals in the strong-field, weak-field, impulsive, and nonimpulsive regimes.

Mixed quantum-classical description of spectroscopy of dissipative systems

Mixed quantum-classical statistical mechanics is employed to calculate dipole moment correlation function and linear absorption spectra. A quantum two-level subsystem interacting with quantum vibrations (primary oscillators) which in turn are coupled to a classical bath composed of infinite set of harmonic oscillators is used as a dissipative system. Starting with mixed quantum-classical Liouville equation for the evaluation of the mixed quantum-classical dipole moment correlation function and using coherent states and the inverse of Baker-Campbell-Hausdorf formula to evaluate the trace over the primary oscillators, whereby, a closed analytical expression for the electronic dipole moment correlation function is obtained. Illustrations of several absorption spectra at different temperatures are provided. An approximate optical four-point correlation is obtained in the high temperature limit. A strategy for deriving an exact optical four-point correlation is suggested.

Transient Phenomena in Time- and Frequency-Gated Spontaneous Emission

The effect of overlapping pump and gate pulses on time- and frequency-gated spontaneous emission spectra is explored for a model of material dynamics that accounts for strong nonadiabatic and electron-vibrational coupling effects, vibrational relaxation, and optical dephasing, thus representing characteristic features of photoinduced excited-state dynamics in large molecules in the gas phase or in condensed phases. The behaviors of the sequential, coherent, and doorway-window contributions to the spontaneous emission spectrum are studied separately. The interrelation between the sequential and coherent contributions is demonstrated to be sensitive to the carrier frequencies of the pump and gate pulses and also to the optical dephasing rate, opening the possibility of an experimental determination of the latter. The coherent contribution is shown to dominate the spectrum at specific emission frequencies.

Time- and frequency-resolved fluorescence spectra of nonadiabatic dissipative systems: What photons can tell us

The monitoring of the excited-state dynamics by time- and frequency-resolved spontaneous emission spectroscopy has been studied in detail for a model exhibiting an excited-state curve crossing. The model represents characteristic aspects of the photoinduced ultrafast dynamics in large molecules in the gas or condensed phases and accounts for strong nonadiabatic and electron-vibrational coupling effects, as well as for vibrational relaxation and optical dephasing. A comprehensive overview of the dependence of spontaneous emission spectra on the characteristics of the excitation and detection processes (such as carrier frequencies, pump/gate pulse durations, as well as optical dephasing) is presented. A systematic comparison of ideal spectra, which provide simultaneously perfect time and frequency resolution and thus contain maximal information on the system dynamics, with actually measurable time- and frequency-gated spectra has been carried out. The calculations of real time- and frequency-gated spectra demonstrate that complementary information on the excited-state dynamics can be extracted when the duration of the gate pulse is varied.