Detection of Electronic and Vibrational Coherence Effects in Electron-Transfer Systems by Femtosecond Time-Resolved Fluorescence Spectroscopy: Theoretical Aspects

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.

Modeling Excited-State Proton Transfer Using the Lindblad Equation: Quantification of Time-Resolved Spectroscopy with Mechanistic Insights

The quantum dynamics of excited-state intramolecular proton transfer (ESIPT) is studied using a multilevel vibronic Hamiltonian and the Lindblad master equation. We simulate time-resolved fluorescence spectroscopy of 2-(2'-hydroxyphenyl) benzothiazole (HBT) and 10-hydroxybenzo[h]quinoline (HBQ), which suggests that the underlying mechanism behind the initial ultrafast rise and decay in the spectra is electronic state population that evolves simultaneously with proton wave packet dynamics. The results predict that the initial rise and decay signals at different wavelengths vary significantly with system properties in terms of their shape, the time, and the intensity of the maximum. These findings provide clues for data interpretation, mechanism validation, and control of the dynamics, and the model serves as an attempt toward clarifying ESIPT by direct comparison to time-resolved spectroscopy.

The fate of molecular excited states: modeling donor-acceptor dyes

We investigate the relaxation of a coherently excited molecule in the Redfield approximation. The molecular model, parametrized to describe donor-acceptor dyes that represent a large family of molecules of interest...

Interplay of vibrational wavepackets during an ultrafast electron transfer reaction

Electron transfer reactions facilitate energy transduction and photoredox processes in biology and chemistry. Recent findings show that molecular vibrations can enable the dramatic acceleration of some electron transfer reactions, and control it by suppressing and enhancing reaction paths. Here, we report ultrafast spectroscopy experiments and quantum dynamics simulations that resolve how quantum vibrations participate in an electron transfer reaction. We observe ballistic electron transfer (~30 fs) along a reaction coordinate comprising high-frequency promoting vibrations. Along another vibrational coordinate, the system becomes impulsively out of equilibrium as a result of the electron transfer reaction. This leads to the generation (by the electron transfer reaction, not the laser pulse) of a new vibrational coherence along this second reaction coordinate in a mode associated with the reaction product. These results resolve a complex reaction trajectory composed of multiple vibrational coordinates that, like a sequence of ratchets, progressively diminish the recurrence of the reactant state.

Monitoring of Nonadiabatic Effects in Individual Chromophores by Femtosecond Double-Pump Single-Molecule Spectroscopy: A Model Study

We explore, by theoretical modeling and computer simulations, how nonadiabatic couplings of excited electronic states of a polyatomic chromophore manifest themselves in single-molecule signals on femtosecond timescales. The chromophore is modeled as a system with three electronic states (the ground state and two non-adiabatically coupled excited states) and a Condon-active vibrational mode which, in turn, is coupled to a harmonic oscillator heat bath. For this system, we simulate double-pump single-molecule signals with fluorescence detection for different system-field interaction strengths, from the weak-coupling regime to the strong-coupling regime. While the signals are determined by the coherence of the electronic density matrix in the weak-coupling regime, they are determined by the populations of the electronic density matrix in the strong-coupling regime. As a consequence, the signals in the strong coupling regime allow the monitoring of nonadiabatic electronic population dynamics and are robust with respect to temporal inhomogeneity of the optical gap, while signals in the weak-coupling regime are sensitive to fluctuations of the optical gap and do not contain information on the electronic population dynamics.

Dissipative dynamics at conical intersections: simulations with the hierarchy equations of motion method

The effect of a dissipative environment on the ultrafast nonadiabatic dynamics at conical intersections is analyzed for a two-state two-mode model chosen to represent the S₂(ππ*)–S₁(nπ*) conical intersection in pyrazine (the system) which is bilinearly coupled to infinitely many harmonic oscillators in thermal equilibrium (the bath). The system–bath coupling is modeled by the Drude spectral function. The equation of motion for the reduced density matrix of the system is solved numerically exactly with the hierarchy equation of motion method using graphics-processor-unit (GPU) technology. The simulations are valid for arbitrary strength of the system–bath coupling and arbitrary bath memory relaxation time. The present computational studies overcome the limitations of weak system–bath coupling and short memory relaxation time inherent in previous simulations based on multi-level Redfield theory [A. Kühl and W. Domcke, J. Chem. Phys. 2002, 116, 263]. Time evolutions of electronic state populations and time-dependent reduced probability densities of the coupling and tuning modes of the conical intersection have been obtained. It is found that even weak coupling to the bath effectively suppresses the irregular fluctuations of the electronic populations of the isolated two-mode conical intersection. While the population of the upper adiabatic electronic state (S₂) is very efficiently quenched by the system–bath coupling, the population of the diabatic ππ* electronic state exhibits long-lived oscillations driven by coherent motion of the tuning mode. Counterintuitively, the coupling to the bath can lead to an enhanced lifetime of the coherence of the tuning mode as a result of effective damping of the highly excited coupling mode, which reduces the strong mode–mode coupling inherent to the conical intersection. The present results extend previous studies of the dissipative dynamics at conical intersections to the nonperturbative regime of system–bath coupling. They pave the way for future first-principles simulations of femtosecond time-resolved four-wave-mixing spectra of chromophores in condensed phases which are nonperturbative in the system dynamics, the system–bath coupling as well as the field-matter coupling.

Ultrafast photochemistry in liquids

Ultrafast photochemical processes can occur in parallel with the relaxation of the optically populated excited state toward equilibrium. The latter involves both intra- and intermolecular modes, namely vibrational and solvent coordinates, and takes place on timescales ranging from a few tens of femtoseconds to up to hundreds of picoseconds, depending on the system. As a consequence, the reaction dynamics can substantially differ from those usually measured with slower photoinduced processes occurring from equilibrated excited states. For example, the decay of the excited-state population may become strongly nonexponential and depend on the excitation wavelength, contrary to the Kasha and Vavilov rules. In this article, we first give a brief account of our current understanding of vibrational and solvent relaxation processes. We then present an overview of important classes of ultrafast photochemical reactions, namely electron and proton transfer as well as isomerization, and illustrate with several examples how nonequilibrium effects can affect their dynamics.

Two-dimensional optical three-pulse photon echo spectroscopy. II. Signatures of coherent electronic motion and exciton population transfer in dimer two-dimensional spectra

Using the nonperturbative approach to the calculation of nonlinear optical spectra developed in a foregoing paper [Mancal et al., J. Chem. Phys. 124, 234504 (2006), preceding paper], calculations of two-dimensional electronic spectra of an excitonically coupled dimer model system are presented. The dissipative exciton transfer dynamics is treated within the Redfield theory and energetic disorder within the molecular ensemble is taken into account. The manner in which the two-dimensional spectra reveal electronic couplings in the aggregate system and the evolution of the spectra in time is studied in detail. Changes in the intensity and shape of the peaks in the two-dimensional relaxation spectra are related to the coherent and dissipative dynamics of the system. It is shown that coherent electronic motion, an electronic analog of a vibrational wave packet, can manifest itself in two-dimensional optical spectra of molecular aggregate systems as a periodic modulation of both the diagonal and off-diagonal peaks.

Quantum modeling of transient infrared spectra reflecting photoinduced electron-transfer dynamics

A theoretical description of transient vibrational spectra following the impulsive optical excitation of a molecular system is presented. The approach combines the nonsecular evaluation of the Redfield equations to describe the dissipative dynamics of the system with an efficient implementation of the doorway-window formalism to calculate optical pump/infrared probe (vis/IR) spectra. Both parts of the calculation scale with N2, thus facilitating the treatment of systems with a dimension up to 10(4). The formulation is applied to a simple model of photoinduced electron transfer, which takes into account two coupled electronic states and a single anharmonic vibrational mode. Despite its simplicity, the model is found to exhibit quite complex electronic and vibrational relaxation dynamics, which in turn give rise to rather complex time- and frequency-resolved vis/IR spectra. Interestingly, the calculated IR spectra of the electron-transfer system predict the appearance of novel vibronically induced sidebands, which may even dominate the spectrum at early times.

Detecting electronic coherence in excited-state electron transfer in fluorinated benzenes

Photoinduced electron transfer (ET) in heavily fluorinated benzenes in solution has been studied with broadband transient absorption spectroscopy. Spectrally resolved kinetics exhibit oscillations with amplitude up to 70% of the signal. The oscillation frequency is specific for each probe; in addition, for pentafluorobenzene it markedly depends on solvent, being 86 cm(-1) in hexane and 94 cm(-1) in acetonitrile. We argue that the observed behavior is not related to vibrational coherences, but originate from coherent electronic motion between an optically excited and an ET state.

Electronic decoherence time for non-Born-Oppenheimer trajectories

An expression is obtained for the electronic decoherence time of the reduced density electronic matrix in mixed quantum-classical molecular-dynamics simulations. The result is obtained by assuming that decoherence is dominated by the time dependence of the overlap of minimum-uncertainty packets and then maximizing the rate with respect to the parameters of the wave packets. The expression for the decay time involves quantities readily available in non-Born-Oppenheimer molecular-dynamics simulations, and it is shown to have a reasonable form when compared with two other formulas for the decay time that have been previously proposed.

Exciton Analysis in 2D Electronic Spectroscopy

A theoretical description of femtosecond two-dimensional electronic spectroscopy of multichromophoric systems is presented. Applying the stationary phase approximation to the calculation of photon echo spectra and taking into account exciton relaxation processes, we obtain an analytic expression for numerical simulations of time- and frequency-resolved 2D photon echo signals. The delocalization of one-exciton states, spatial overlaps between the probability densities of different excitonic states, and their influences on both one- and two-dimensional electronic spectra are studied. The nature of the off-diagonal cross-peaks and the time evolution of both diagonal and off-diagonal peak amplitudes are discussed in detail by comparing experimentally measured and theoretically simulated 2D spectra of the natural Fenna-Matthews-Olson (FMO) photosynthetic light-harvesting complex. We find that there are two noncascading exciton energy relaxation pathways.

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.