Femtosecond pump-probe spectroscopy of electron-transfer systems: a nonperturbative approach

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.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

Modeling and correction of distorted two-dimensional Fourier transform spectra from pixelated pulse shaping devices

Two-dimensional Fourier transform spectra of a three level model system are simulated using a non-perturbative density matrix formalism. The electric field distortions resultant from using pixelated pulse shaping devices to produce phase-locked pulse pairs are modeled and the effects on the recovered spectra are examined. To minimize spectral distortions, a temporal filtering scheme is employed which eliminates contributions from spurious sample polarizations.

Efficient Simulation of Three-Pulse Photon-Echo Signals with Application to the Determination of Electronic Coupling in a Bacterial Photosynthetic Reaction Center

A time-nonlocal quantum master equation coupled with a perturbative scheme to evaluate the third-order polarization in the phase-matching direction k(s) = -k(1) + k(2) + k(3) is used to efficiently simulate three-pulse photon-echo signals. The present method is capable of describing photon-echo peak shifts including pulse overlap and bath memory effects. In addition, the method treats the non-Markovian evolution of the density matrix and the third-order polarization in a consistent manner, thus is expected to be useful in systems with rapid and complex dynamics. We apply the theoretical method to describe one- and two-color three-pulse photon-echo peak shift experiments performed on a bacterial photosynthetic reaction center and demonstrate that, by properly incorporating the pulse overlap effects, the method can be used to describe simultaneously all peak shift experiments and determine the electronic coupling between the localized Q(y) excitations on the bacteriopheophytin (BPhy) and accessory bateriochlorophyll (BChl) in the reaction center. A value of J = 250 cm(-1) is found for the coupling between BPhy and BChl.

Site-selective photoinduced electron transfer from alcoholic solvents to the chromophore facilitated by hydrogen bonding: a new fluorescence quenching mechanism

Solute-solvent intermolecular photoinduced electron transfer (ET) reaction was proposed to account for the drastic fluorescence quenching behaviors of oxazine 750 (OX750) chromophore in protic alcoholic solvents. According to our theoretical calculations for the hydrogen-bonded OX750-(alcohol)(n) complexes using the time-dependent density functional theory (TDDFT) method, we demonstrated that the ET reaction takes place from the alcoholic solvents to the chromophore and the intermolecular ET passing through the site-specific intermolecular hydrogen bonds exhibits an unambiguous site selectivity. In our motivated experiments of femtosecond time-resolved stimulated emission pumping fluorescence depletion spectroscopy (FS TR SEP FD), it could be noted that the ultrafast ET reaction takes place as fast as 200 fs. This ultrafast intermolecular photoinduced ET is much faster than the diffusive solvation process, and even significantly faster than the intramolecular vibrational redistribution (IVR) process of the OX750 chromophore. Therefore, the ultrafast intermolecular ET should be coupled with the hydrogen-bonding dynamics occurring in the sub-picosecond time domain. We theoretically demonstrated for the first time that the selected hydrogen bonds are transiently strengthened in the excited states for facilitating the ultrafast solute-solvent intermolecular ET reaction.

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.

Two-dimensional optical three-pulse photon echo spectroscopy. I. Nonperturbative approach to the calculation of spectra

The nonperturbative approach to the calculation of nonlinear optical spectra of Seidner et al. [J. Chem. Phys. 103, 3998 (1995)] is extended to describe four-wave mixing experiments. The system-field interaction is treated nonperturbatively in the semiclassical dipole approximation, enabling a calculation of third order nonlinear spectroscopic signals directly from molecular dynamics and an efficient modeling of multilevel systems exhibiting relaxation and transfer phenomena. The method, coupled with the treatment of dynamics within the Bloch model, is illustrated by calculations of the two-dimensional three-pulse photon echo spectra of a simple model system-a two-electronic-level molecule. The nonperturbative calculations reproduce well-known results obtained by perturbative methods. Technical limitations of the nonperturbative approach in dealing with a dynamic inhomogeneity are discussed, and possible solutions are suggested. An application of the approach to an excitonically coupled dimer system with emphasis on the manifestation of complex exciton dynamics in two-dimensional optical spectra is presented in paper II Pisliakov et al. [J. Chem. Phys. 124, 234505 (2006), following paper].

Three-pulse photon echo of an excitonic dimer modeled via Redfield theory

In this article the third-order response of an excitonically coupled dimer is studied. The three-pulse photon echo signals were calculated by extracting polarization components from the total polarization in the corresponding phase-matched directions. The total nonlinear response was obtained by numeric propagation of the density matrix, with the exciton-vibrational coupling modeled via Redfield relaxation theory. The full two-dimensional three-pulse photon echo signals and the peak shift were analyzed in terms of the density-matrix dynamics of coherence dephasing and population relaxation. The location of the two-exciton state was found to be essential for proper modeling of the three-pulse photon echo. In particular, an oscillation in the three-pulse photon echo peak shift is found if the two-exciton state is displaced. The oscillations can be related to the dynamics of the one-exciton coherences.

Mixed quantum-classical Redfield master equation

Redfield master equation is derived from mixed quantum-classical Liouville equation using product initial conditions. Simple two-level system example is given and comparison with Fermi golden rule is made.

Efficient method for the calculation of time- and frequency-resolved four-wave mixing signals and its application to photon-echo spectroscopy

An efficient method has been developed for the calculation of third-order time- and frequency-resolved optical signals. To obtain the general four-wave mixing signal, seven auxiliary density matrices have to be propagated in time. For the special cases of two-pulse photon-echo and transient-grating signals, two or three density matrices, respectively, are required. The method is limited to weak laser fields (it is thus valid within the third-order perturbation theory) but allows for any pulse durations and automatically accounts for pulse-overlap effects. To illustrate the method, we present the explicit derivation of the three-pulse photon-echo signal. Any other third-order optical signal can be calculated in the same manner. As an example, two- and three-pulse photon-echo and transient-grating signals for a weakly damped displaced harmonic oscillator have been calculated.

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.

Exciton exciton annihilation dynamics in chromophore complexes. II. Intensity dependent transient absorption of the LH2 antenna system

Using the multiexciton density matrix theory of excitation energy transfer in chromophore complexes developed in a foregoing paper [J. Chem. Phys. 118, 746 (2003)], the computation of ultrafast transient absorption spectra is presented. Beside static disorder and standard mechanisms of excitation energy dissipation the theory incorporates exciton exciton annihilation (EEA) processes. To elucidate signatures of EEA in intensity dependent transient absorption data the approach is applied to the B850 ring of the LH2 found in rhodobacter sphaeroides. As main indications for two-exciton population and resulting EEA we found (i) a weakening of the dominant single-exciton bleaching structure in the transient absorption, and (ii) an intermediate suppression of long-wavelength and short-wavelength shoulders around the bleaching structure. The suppression is caused by stimulated emission from the two-exciton to the one-exciton state and the return of the shoulders follows from a depletion of two-exciton population according to EEA. The EEA-signature survives as a short-wavelength shoulder in the transient absorption if orientational and energetic disorder are taken into account. Therefore, the observation of the EEA-signatures should be possible when doing frequency resolved transient absorption experiments with a sufficiently strongly varying pump-pulse intensity.