Mapping of exciton–exciton annihilation in a molecular dimer via fifth-order femtosecond two-dimensional spectroscopy

Exciton annihilation in molecular aggregates suppressed through qu antum interference

Exciton-exciton annihilation (EEA), an important loss channel in optoelectronic devices and photosynthetic complexes, has conventionally been assumed to be an incoherent, diffusion-limited process. Here we challenge this assumption by experimentally demonstrating the ability to control EEA in molecular aggregates using the quantum phase relationships of excitons. We employed time-resolved photoluminescence microscopy to independently determine exciton diffusion constants and annihilation rates in two substituted perylene diimide aggregates featuring contrasting excitonic phase envelopes. Low-temperature EEA rates were found to differ by more than two orders of magnitude for the two compounds, despite comparable diffusion constants. Simulated rates based on a microscopic theory, in excellent agreement with experiments, rationalize this EEA behaviour based on quantum interference arising from the presence or absence of spatial phase oscillations of delocalized excitons. These results offer an approach for designing molecular materials using quantum interference where low annihilation can coexist with high exciton concentrations and mobilities.

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.

Fluorescence-Detected Two-Quantum Photon Echoes via Cogwheel Phase Cycling

Two-dimensional (2D) electronic spectroscopy can separate homogeneous and inhomogeneous broadening. While established methods usually probe a photon-echo signal, i.e., a third-order response, to access the homogeneous line width of singly excited states, the homogeneous line width of doubly excited states remained spectroscopically inaccessible. Here we demonstrate the acquisition of two-quantum (2Q) photon echoes using fluorescence-detected 2D spectroscopy. In these eighth-order signals, 2Q coherences are rephased with themselves, leading to line-narrowed 2Q-2Q 2D spectra. By using cogwheel phase cycling, adapted from nuclear magnetic resonance spectroscopy, we isolate the 2Q-2Q 2D spectra of a squaraine dimer and a squaraine polymer and verify the same selectivity of cogwheel phase cycling compared to traditional "nested" phase cycling. The observed difference, between the two systems, in the homogeneous line width of the biexciton can be rationalized as a signature of the interplay of exciton-exciton annihilation, exciton diffusion, and exciton delocalization.

Phase stable, shot-to-shot measurement of third- and fifth-order two-quantum correlation spectra using a pulse shaper in the pump-probe geometry

We demonstrate the first phase stable measurement of a third-order 2Q spectrum using a pulse shaper in the pump-probe geometry. This measurement was achieved by permuting the time-ordering of the pump pulses, thus rearranging the signal pathways that are emitted in the probe direction. The third-order 2Q spectrum is self-heterodyned by the probe pulse. Using this method, one can interconvert between a 1Q experiment and a 2Q experiment by simply reprogramming a pulse shaper or delay stage. We also measure a fifth-order absorptive 2Q spectrum in the pump-probe geometry, which contains similar information as a third-order experiment but does not suffer from dispersive line shapes. To do so, we introduce methods to minimize saturation-induced artifacts of the pulse shaper, improving fifth-order signals. These techniques add new capabilities for 2D spectrometers that use pulse shapers in the pump-probe beam geometry.

Third-order pump-probe spectroscopy applied to molecular dimers: characterization of relaxation dynamics and exciton–exciton annihilation

Exciton–exciton annihilation in a dimer, described within the basis of localizes monomer states.

Intrachain and Interchain Exciton-Exciton Annihilation in Donor-Acceptor Copolymers

PM6 is a donor-acceptor copolymer widely used in high-efficiency organic solar cells (OSC). Here, we report the observation of exciton-exciton annihilation (EEA) in both PM6 copolymer solution and film. We find that the EEA mechanisms are distinct in these two systems. In the PM6 solution, EEA occurs via the long-range dipole-dipole interaction between intrachain excitons, while in the PM6 film it takes place between the interchain polaron pairs (spatially indirect excitons) via the short-range Coulomb interaction. The quantum yield of free polaron formation in the film is also quantified, which indicates that EEA leads to extra interchain polaron-pair dissociation, accounting for the enhanced free polaron formation.

Automatic Feynman diagram generation for nonlinear optical spectroscopies and application to fifth-order spectroscopy with pulse overlaps

Perturbative nonlinear optical spectroscopies are powerful methods to understand the dynamics of excitonic and other condensed phase systems. Feynman diagrams have long provided the essential tool to understand and interpret experimental spectra and to organize the calculation of spectra for model systems. When optical pulses are strictly time ordered, only a small number of diagrams contribute, but in many experiments, pulse-overlap effects are important for interpreting results. When pulses overlap, the number of contributing diagrams can increase rapidly, especially with higher order spectroscopies, and human error is especially likely when attempting to write down all the diagrams. We present an automated Diagram Generator (DG) that generates all the Feynman diagrams needed to calculate any nth-order spectroscopic signal. We characterize all perturbative nonlinear spectroscopies by their associated phase-discrimination condition as well as the time intervals where pulse amplitudes are nonzero. Although the DG can be used to automate impulsive calculations, its greatest strength lies in automating finite pulse calculations where pulse overlaps are important. We consider third-order transient absorption spectroscopy and fifth-order exciton-exciton interaction 2D (EEI2D) spectroscopy, which are described by six or seven diagrams in the impulsive limit, respectively, but 16 or 240 diagrams, respectively, when pulses overlap. The DG allows users to automatically include all relevant diagrams at a relatively low computational cost, since the extra diagrams are only generated for the inter-pulse delays where they are relevant. For EEI2D spectroscopy, we show the important effects of including the overlap diagrams.

Efficient numerical method for predicting nonlinear optical spectroscopies of open systems

Nonlinear optical spectroscopies are powerful tools for probing quantum dynamics in molecular and nanoscale systems. While intuition about ultrafast spectroscopies is often built by considering impulsive optical pulses, actual experiments have finite-duration pulses, which can be important for interpreting and predicting experimental results. We present a new freely available open source method for spectroscopic modeling, called Ultrafast Ultrafast (UF²) spectroscopy, which enables computationally efficient and convenient prediction of nonlinear spectra, such as treatment of arbitrary finite duration pulse shapes. UF² is a Fourier-based method that requires diagonalization of the Liouvillian propagator of the system density matrix. We also present a Runge-Kutta-Euler (RKE) direct propagation method. We include open system dynamics in the secular Redfield, full Redfield, and Lindblad formalisms with Markovian baths. For non-Markovian systems, the degrees of freedom corresponding to memory effects are brought into the system and treated nonperturbatively. We analyze the computational complexity of the algorithms and demonstrate numerically that, including the cost of diagonalizing the propagator, UF² is 20-200 times faster than the direct propagation method for secular Redfield models with arbitrary Hilbert space dimension; it is similarly faster for full Redfield models at least up to system dimensions where the propagator requires more than 20 GB to store; and for Lindblad models, it is faster up to Hilbert space dimension near 100 with speedups for small systems by factors of over 500. UF² and RKE are part of a larger open source Ultrafast Software Suite, which includes tools for automatic generation and calculation of Feynman diagrams.

Fifth-order two-quantum absorptive two-dimensional electronic spectroscopy of CdSe quantum dots

Two-quantum variants of two-dimensional electronic spectroscopy (2DES) have previously been used to characterize multi-exciton interactions in molecules and semiconductor nanostructures though many implementations are limited by phasing procedures or non-resonant signals. We implement 2DES using phase-cycling to simultaneously measure one-quantum and two-quantum spectra in colloidal CdSe quantum dots. In the pump-probe geometry, fully absorptive spectra are automatically acquired by measuring the sum of the rephasing and nonrephasing signals. Fifth-order two-quantum spectroscopy allows for direct access to multi-exciton states that may be obscured in excited state absorption signals due to population relaxation or third-order two-quantum spectra due to the non-resonant response.

Exciton-exciton annihilation in a molecular trimer: Wave packet dynamics and 2D spectroscopy

We theoretically study the exciton-exciton annihilation (EEA) in a molecular trimer MMM. The system is treated within a model of electronic states, and the coupling to a bath is incorporated using the quantum jump method. Two situations of initial excitation are compared. In the first one, a two-photon process populates configurations M^*M^*M and MM^*M^* so that two excitons reside on neighboring monomers M. Then, EEA can immediately proceed. In contrast, if the trimer initially is in the local configuration M^*MM^*, exciton diffusion must occur before the annihilation process can take place. For the trimer, this excitonic motion takes place on a very short time scale. In both cases, wave packets are prepared which show a different quantum dynamics where the latter depends on the couplings and decay rates. It is documented how fifth-order coherent two-dimensional spectroscopy can be used to directly map the EEA as a function of time.

Signatures of exciton dynamics and interaction in coherently and fluorescence-detected four- and six-wave-mixing two-dimensional electronic spectroscopy

Two-dimensional electronic spectroscopy (2DES) can be realized in increasing nonlinear orders of interaction with the electric field, bringing new information about single- and multi-particle properties and dynamics. Furthermore, signals can be detected both coherently (C-2DES) and by fluorescence (F-2DES), with fundamental and practical differences. We directly compare the simultaneous measurements of four- and six-wave mixing C-2DES and F-2DES on an excitonic heterodimer of squaraine molecules. Spectral features are described in increasing orders of nonlinearity by an explicit excitonic model. We demonstrate that the four-wave-mixing spectra are sensitive to one-exciton energies, their delocalization and dynamics, while the six-wave-mixing spectra include information on bi-exciton and higher excited states including the state energies, electronic coupling, and exciton-exciton annihilation. We focus on the possibility to extract the dynamics arising from exciton-exciton interaction directly from the six-wave-mixing spectra. To this end, in analogy to previously demonstrated fifth-order coherently detected exciton-exciton-interaction 2DES (EEI2D spectroscopy), we introduce a sixth-order fluorescence-detected EEI2D spectroscopy variant.

Effects of high pulse intensity and chirp in two-dimensional electronic spectroscopy of an atomic vapor

The effects of high pulse intensity and chirp on two-dimensional electronic spectroscopy signals are experimentally investigated in the highly non-perturbative regime using atomic rubidium vapor as clean model system. Data analysis is performed based on higher-order Feynman diagrams and non-perturbative numerical simulations of the system response. It is shown that higher-order contributions may lead to a fundamental change of the static appearance and beating-maps of the 2D spectra and that chirped pulses enhance or suppress distinct higher-order pathways. We further give an estimate of the threshold intensity beyond which the high-intensity effects become visible for the system under consideration.

A wave packet picture of exciton-exciton annihilation: Molecular dimer dynamics

The usual view of exciton-exciton annihilation (EEA) processes in molecular aggregates is based on locally excited states of the monomer units. However, the corresponding localized configurations can only be assumed if the system is in a coherent superposition of eigenstates, i.e., a wave packet. We study a molecular dimer and focus on the characterization of EEA by a wave packet motion induced in the system by ultrashort pulse excitation. Here, coherences that appear are destroyed by dissipation processes. We discuss the influence of interband and intraband relaxation on the dynamics. The states that participate in the annihilation process are directly accessible by fifth-order optical two-dimensional spectroscopy. Such spectra are calculated, and spectral features are related to the annihilation process.

From wavelike to sub-diffusive motion: exciton dynamics and interaction in squaraine copolymers of varying length

Exciton transport and exciton-exciton interactions in molecular aggregates and polymers are of great importance in natural photosynthesis, organic electronics, and related areas of research. Both the experimental observation and theoretical description of these processes across time and length scales, including the transition from the initial wavelike motion to the following long-range exciton transport, are highly challenging. Therefore, while exciton dynamics at small scales are often treated explicitly, long-range exciton transport is typically described phenomenologically by normal diffusion. In this work, we study the transition from wavelike to diffusive motion of interacting exciton pairs in squaraine copolymers of varying length. To this end we use a combination of the recently introduced exciton-exciton-interaction two-dimensional (EEI2D) electronic spectroscopy and microscopic theoretical modelling. As we show by comparison with the model, the experimentally observed kinetics include three phases, wavelike motion dominated by immediate exciton-exciton annihilation (10-100 fs), sub-diffusive behavior (0.1-10 ps), and excitation relaxation (0.01-1 ns). We demonstrate that the key quantity for the transition from wavelike to diffusive dynamics is the exciton delocalization length relative to the length of the polymer: while in short polymers wavelike motion of rapidly annihilating excitons dominates, in long polymers the excitons become locally trapped and exhibit sub-diffusive behavior. Our findings indicate that exciton transport through conjugated systems emerging from the excitonic structure is generally not governed by normal diffusion. Instead, to characterize the material transport properties, the diffusion presence and character should be determined.

Interplay between structural hierarchy and exciton diffusion in artificial light harvesting

Unraveling the nature of energy transport in multi-chromophoric photosynthetic complexes is essential to extract valuable design blueprints for light-harvesting applications. Long-range exciton transport in such systems is facilitated by a combination of delocalized excitation wavefunctions (excitons) and exciton diffusion. The unambiguous identification of the exciton transport is intrinsically challenging due to the system's sheer complexity. Here we address this challenge by employing a spectroscopic lab-on-a-chip approach: ultrafast coherent two-dimensional spectroscopy and microfluidics working in tandem with theoretical modeling. We show that at low excitation fluences, the outer layer acts as an exciton antenna supplying excitons to the inner tube, while under high excitation fluences the former converts its functionality into an exciton annihilator which depletes the exciton population prior to any exciton transfer. Our findings shed light on the excitonic trajectories across different sub-units of a multi-layered artificial light-harvesting complex and underpin their great potential for directional excitation energy transport.