Analysis of cross peaks in two-dimensional electronic photon-echo spectroscopy for simple models with vibrations and dissipation

Theoretical model of femtosecond coherence spectroscopy of vibronic excitons in molecular aggregates

When used as pump pulses in transient absorption spectroscopy measurements, femtosecond laser pulses can produce oscillatory signals known as quantum beats. The quantum beats arise from coherent superpositions of the states of the sample and are best studied in the Fourier domain using Femtosecond Coherence Spectroscopy (FCS), which consists of one-dimensional amplitude and phase plots of a specified oscillation frequency as a function of the detection frequency. Prior works have shown ubiquitous amplitude nodes and π phase shifts in FCS from excited-state vibrational wavepackets in monomer samples. However, the FCS arising from vibronic-exciton states in molecular aggregates have not been studied theoretically. Here, we use a model of vibronic-exciton states in molecular dimers based on displaced harmonic oscillators to simulate FCS for dimers in two important cases. Simulations reveal distinct spectral signatures of excited-state vibronic-exciton coherences in molecular dimers that may be used to distinguish them from monomer vibrational coherences. A salient result is that, for certain relative orientations of the transition dipoles, the key resonance condition between the electronic coupling and the frequency of the vibrational mode may yield strong enhancement of the quantum-beat amplitude and, perhaps, also cause a significant decrease of the oscillation frequency to a value far lower than the vibrational frequency. Future studies using these results will lead to new insights into the excited-state coherences generated in photosynthetic pigment-protein complexes.

Plasmon mediated coherent population oscillations in molecular aggregates

The strong coherent coupling of quantum emitters to vacuum fluctuations of the light field offers opportunities for manipulating the optical and transport properties of nanomaterials, with potential applications ranging from ultrasensitive all-optical switching to creating polariton condensates. Often, ubiquitous decoherence processes at ambient conditions limit these couplings to such short time scales that the quantum dynamics of the interacting system remains elusive. Prominent examples are strongly coupled exciton-plasmon systems, which, so far, have mostly been investigated by linear optical spectroscopy. Here, we use ultrafast two-dimensional electronic spectroscopy to probe the quantum dynamics of J-aggregate excitons collectively coupled to the spatially structured plasmonic fields of a gold nanoslit array. We observe rich coherent Rabi oscillation dynamics reflecting a plasmon-driven coherent exciton population transfer over mesoscopic distances at room temperature. This opens up new opportunities to manipulate the coherent transport of matter excitations by coupling to vacuum fields.

Phonon Signatures in Photon Correlations

We show that the second-order, two-time correlation functions for phonons and photons emitted from a vibronic molecule in a thermal bath result in bunching and antibunching (a purely quantum effect), respectively. Signatures relating to phonon exchange with the environment are revealed in photon-photon correlations. We demonstrate that cross-correlation functions have a strong dependence on the order of detection giving insight into how phonon dynamics influences the emission of light. This work offers new opportunities to investigate quantum effects in condensed-phase molecular systems.

Vibrational response functions for multidimensional electronic spectroscopy in nonadiabatic models

The interplay of nuclear and electronic dynamics characterizes the multidimensional electronic spectra of various molecular and solid-state systems. Theoretically, the observable effect of such interplay can be accounted for by response functions. Here, we report analytical expressions for the response functions corresponding to a class of model systems. These are characterized by coupling between the diabatic electronic states and the vibrational degrees of freedom, resulting in linear displacements of the corresponding harmonic oscillators, and by nonadiabatic couplings between pairs of diabatic states. In order to derive the linear response functions, we first perform the Dyson expansion of the relevant propagators with respect to the nonadiabatic component of the Hamiltonian, then derive and expand with respect to the displacements the propagators at given interaction times, and finally provide analytical expressions for the time integrals that lead to the different contributions to the linear response function. The approach is then applied to the derivation of third-order response functions describing different physical processes: ground state bleaching, stimulated emission, excited state absorption, and double quantum coherence. Comparisons between the results obtained up to sixth order in the Dyson expansion and independent numerical calculation of the response functions provide evidence of the series convergence in a few representative cases.

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.

Vibrational response functions for multidimensional electronic spectroscopy in the adiabatic regime: A coherent-state approach

Multi-dimensional spectroscopy represents a particularly insightful tool for investigating the interplay of nuclear and electronic dynamics, which plays an important role in a number of photophysical processes and photochemical reactions. Here, we present a coherent state representation of the vibronic dynamics and of the resulting response functions for the widely used linearly displaced harmonic oscillator model. Analytical expressions are initially derived for the case of third-order response functions in an N-level system, with ground state initialization of the oscillator (zero-temperature limit). The results are then generalized to the case of Mth order response functions, with arbitrary M. The formal derivation is translated into a simple recipe, whereby the explicit analytical expressions of the response functions can be derived directly from the Feynman diagrams. We further generalize to the whole set of initial coherent states, which form an overcomplete basis. This allows one, in principle, to derive the dependence of the response functions on arbitrary initial states of the vibrational modes and is here applied to the case of thermal states. Finally, a non-Hermitian Hamiltonian approach is used to include in the above expressions the effect of vibrational relaxation.

Broadband visible two-dimensional spectroscopy of molecular dyes

Two-dimensional Fourier transform spectroscopy is a promising technique to study ultrafast molecular dynamics. Similar to transient absorption spectroscopy, a more complete picture of the dynamics requires broadband laser pulses to observe transient changes over a large enough bandwidth, exceeding the inhomogeneous width of electronic transitions, as well as the separation between the electronic or vibronic transitions of interest. Here, we present visible broadband 2D spectra of a series of dye molecules and report vibrational coherences with frequencies up to ∼1400 cm^-1 that were obtained after improvements to our existing two-dimensional Fourier transform setup [Al Haddad et al., Opt. Lett. 40, 312-315 (2015)]. The experiment uses white light from a hollow core fiber, allowing us to acquire 2D spectra with a bandwidth of 200 nm, in a range between 500 and 800 nm, and with a temporal resolution of 10-15 fs. 2D spectra of nile blue, rhodamine 800, terylene diimide, and pinacyanol iodide show vibronic spectral features with at least one vibrational mode and reveal information about structural motion via coherent oscillations of the 2D signals during the population time. For the case of pinacyanol iodide, these observations are complemented by its Raman spectrum, as well as the calculated Raman activity at the ground- and excited-state geometry.

A quantum Langevin equation approach for two-dimensional electronic spectra of coupled vibrational and electronic dynamics

We present an efficient method to simulate two-dimensional (2D) electronic spectra of condensed-phase systems with an emphasis on treating quantum nuclear wave packet dynamics explicitly. To this end, we combine a quantum Langevin equation (QLE) approach for dissipation and a perturbative scheme to calculate three-pulse photon-echo polarizations based on wave packet dynamics under the influence of external fields. The proposed dynamical approach provides a consistent description of nuclear quantum dynamics, pulse-overlap effects, and vibrational relaxation, enabling simulations of 2D electronic spectra with explicit and non-perturbative treatment of coupled electronic-nuclear dynamics. We apply the method to simulate 2D electronic spectra of a displaced-oscillator model in the condensed phase and discuss the spectral and temporal evolutions of 2D signals. Our results show that the proposed QLE approach is capable of describing vibrational relaxation, decoherence, and vibrational coherence transfer, as well as their manifestations in spectroscopic signals. Furthermore, vibrational quantum beats specific for excited-state vs ground-state nuclear wave packet dynamics can also be identified. We anticipate that this method will provide a useful tool to conduct theoretical studies of 2D spectroscopy for strong vibronically coupled systems and to elucidate intricate vibronic couplings in complex molecular systems.

Intermolecular conical intersections in molecular aggregates

Conical intersections (CoIns) of multidimensional potential energy surfaces are ubiquitous in nature and control pathways and yields of many photo-initiated intramolecular processes. Such topologies can be potentially involved in the energy transport in aggregated molecules or polymers but are yet to be uncovered. Here, using ultrafast two-dimensional electronic spectroscopy (2DES), we reveal the existence of intermolecular CoIns in molecular aggregates relevant for photovoltaics. Ultrafast, sub-10-fs 2DES tracks the coherent motion of a vibrational wave packet on an optically bright state and its abrupt transition into a dark state via a CoIn after only 40 fs. Non-adiabatic dynamics simulations identify an intermolecular CoIn as the source of these unusual dynamics. Our results indicate that intermolecular CoIns may effectively steer energy pathways in functional nanostructures for optoelectronics.

Spectral Filtering as a Tool for Two-Dimensional Spectroscopy: A Theoretical Model

Two-dimensional optical spectroscopy is a powerful technique for the probing of coherent quantum superpositions. Recently, the finite width of the laser spectrum has been employed to selectively tune experiments for the study of particular coherences. This involves the exclusion of certain transition frequencies, which results in the elimination of specific Liouville pathways. The rigorous analysis of such experiments requires the use of ever more sophisticated theoretical models for the optical spectroscopy of electronic and vibronic systems. Here we develop a nonimpulsive and non-Markovian model, which combines an explicit definition of the laser spectrum, via the equation of motion-phase matching approach (EOM-PMA), with the hierarchical equations of motion (HEOM). This theoretical framework is capable of simulating the 2D spectroscopy of vibronic systems with low frequency modes, coupled to environments of intermediate and slower time scales. In order to demonstrate the spectral filtering of vibronic coherences, we examine the elimination of lower energy peaks from the 2D spectra of a zinc porphyrin monomer upon blue-shifting the laser spectrum. The filtering of Liouville pathways is revealed through the disappearance of peaks from the amplitude spectra for a coupled vibrational mode.

Vibronic coupling in organic semiconductors for photovoltaics

Ultrafast two-dimensional electronic spectroscopy reveals vibronically-assisted coherent charge transport and separation in organic materials and opens up new perspectives for artificial light-to-current conversion.

Linear and third- and fifth-order nonlinear spectroscopies of a charge transfer system coupled to an underdamped vibration

We study hole, electron, and exciton transports in a charge transfer system in the presence of underdamped vibrational motion. We analyze the signature of these processes in the linear and third-, and fifth-order nonlinear electronic spectra. Calculations are performed with a numerically exact hierarchical equations of motion method for an underdamped Brownian oscillator spectral density. We find that combining electron, hole, and exciton transfers can lead to non-trivial spectra with more structure than with excitonic coupling alone. Traces taken during the waiting time of a two-dimensional (2D) spectrum are dominated by vibrational motion and do not reflect the electron, hole, and exciton dynamics directly. We find that the fifth-order nonlinear response is particularly sensitive to the charge transfer process. While third-order 2D spectroscopy detects the correlation between two coherences, fifth-order 2D spectroscopy (2D population spectroscopy) is here designed to detect correlations between the excited states during two different time periods.

Two-dimensional optical spectroscopy of homo- and heterodimers

We theoretically study the two-dimensional (2D) spectroscopy of molecular dimers.

Mutations to R. sphaeroides Reaction Center Perturb Energy Levels and Vibronic Coupling but Not Observed Energy Transfer Rates

The bacterial reaction center is capable of both efficiently collecting and quickly transferring energy within the complex; therefore, the reaction center serves as a convenient model for both energy transfer and charge separation. To spectroscopically probe the interactions between the electronic excited states on the chromophores and their intricate relationship with vibrational motions in their environment, we examine coherences between the excited states. Here, we investigate this question by introducing a series of point mutations within 12 Å of the special pair of bacteriochlorophylls in the Rhodobacter sphaeroides reaction center. Using two-dimensional spectroscopy, we find that the time scales of energy transfer dynamics remain unperturbed by these mutations. However, within these spectra, we detect changes in the mixed vibrational-electronic coherences in these reaction centers. Our results indicate that resonance between bacteriochlorophyll vibrational modes and excitonic energy gaps promote electronic coherences and support current vibronic models of photosynthetic energy transfer.

Simulation of femtosecond two-dimensional electronic spectra of conical intersections

We have simulated femtosecond two-dimensional (2D) electronic spectra for an excited-state conical intersection using the wave-function version of the equation-of-motion phase-matching approach. We show that 2D spectra at fixed values of the waiting time provide information on the structure of the vibronic eigenstates of the conical intersection, while the evolution of the spectra with the waiting time reveals predominantly ground-state wave-packet dynamics. The results show that 2D spectra of conical intersection systems differ significantly from those obtained for chromophores with well separated excited-state potential-energy surfaces. The spectral signatures which can be attributed to conical intersections are discussed.

Quantum diffusion wave-function approach to two-dimensional vibronic spectroscopy

We apply the quantum diffusion wavefunction approach to calculate vibronic two-dimensional (2D) spectra. As an example, we use a system consisting of two electronic states with harmonic oscillator potentials which are coupled to a bath and interact with three time-delayed laser pulses. The first- and second-order perturbative wave functions which enter into the expression for the third-order polarization are determined for a sufficient number of stochastic runs. The wave-packet approach, besides being an alternative technique to calculate the spectra, offers an intuitive insight into the dissipation dynamics and its relation to the 2D vibronic spectra.

Polaron dynamics in two-dimensional photon-echo spectroscopy of molecular rings

We have developed a new approach to the computation of third-order spectroscopic signals of molecular rings, by incorporating the Davydov soliton theory into the nonlinear response function formalism. The Davydov D1 and D Ansätze have been employed to treat the interactions between the excitons and the primary phonons, allowing for a full description of arbitrary exciton-phonon coupling strengths. As an illustration, we have simulated a series of optical 2D spectra for two models of molecular rings.

Two-dimensional vibronic spectroscopy of coherent wave-packet motion

We theoretically study two-dimensional (2D) spectroscopic signals obtained from femtosecond pulse interactions with diatomic molecules. The vibrational wave-packet dynamics is monitored in the signals. During the motion in anharmonic potentials the wave packets exhibit vibrational revivals and fractional revivals which are associated with particular quantum phases. The time-dependent phase changes are identified by inspection of the complex-valued 2D spectra. We use the Na(2) molecule as a numerical example and discuss various pulse sequences which yield information about vibrational level structure and phase relationships in different electronic states.

Interference Effects in Vibronic 2D-Spectra of Diatomic Molecules

We theoretically study interference effects in two-dimensional (2D) vibronic spectra which arise from two electronically excited states taking part in the multi-photon process initiated by femtosecond laser pulses. Therefore, a model is employed which mimiques the situation encountered in many halogen and interhalogen molecules. There, upon excitation from the ground state, an excited bound state and a dissociative state exist which are close in energy. We demonstrate that the different pathways to final states which enter into the third-order polarization result in pronounced interference patterns in the 2D-spectra.

High frequency vibrational modulations in two-dimensional electronic spectra and their resemblance to electronic coherence signatures

In this work we analyze how nuclear coherences modulate diagonal and off-diagonal peaks in two-dimensional electronic spectroscopy. 2D electronic spectra of pinacyanol chloride are measured with 8 fs pulses, which allows coherent excitation of the 1300 cm(-1) vibrational mode. The 2D spectrum reveals both diagonal and off-diagonal peaks related to the vibrational mode. On early time scales, up to 30 fs, coherent dynamics give rise to oscillations in the amplitudes, positions, and shapes of the peaks in the 2D spectrum. We find an anticorrelation between the amplitude and the diagonal width of the two diagonal peaks. The measured data are reproduced with a model incorporating a high frequency mode coupled to an electronic two-level-system. Our results show that these anticorrelated oscillations occur for vibrational wavepackets and not exclusively for electronic coherences as has been assumed previously.