Effects of dynamic disorder on exciton migration: Quantum diffusion, coherences, and energy transfer

Memory effects in the efficiency control of energy transfer under incoherent light excitation in noisy environments

Fluctuations in energy gap and coupling constants between chromophores can play an important role in absorption and energy transfer across a collection of two-level systems. In photosynthesis, light-induced quantum coherence can affect the efficiency of energy transfer to the designated "trap" state. Theoretically, the interplay between fluctuations and coherence has been studied often, employing either a Markovian or a perturbative approximation. In this study, we depart from these approaches to incorporate memory effects by using Kubo's quantum stochastic Liouville equation. We introduce the effects of decay of the created excitation (to the ground state) on the desired propagation and trapping that provides a direction of flow of the excitation. In the presence of light-induced pumping, we establish a relation between the efficiency, the mean survival time, and the correlation decay time of the bath-induced fluctuations. A decrease in the steady-state coherence during the transition from the non-Markovian regime to the Markovian limit results in a decrease in efficiency. As in the well-known Haken-Strobl model, the ratio of the square of fluctuation strength to the rate plays a critical role in determining the mechanism of energy transfer and in shaping the characteristics of the efficiency profile. We recover a connection between the transfer flux and the imaginary part of coherences in both equilibrium and excited bath states, in both correlated and uncorrelated bath models. We uncover a non-monotonic dependence of efficiency on site energy heterogeneity for both correlated and uncorrelated bath models.

Quantum dynamical study of inter-chain exciton transport in a regioregular P3HT model system at finite temperature: HJ vs. H-aggregate models

We report on quantum dynamical simulations of inter-chain exciton transport in a model of regioregular poly(3-hexylthiophene), rr-P3HT, at finite temperature, using the Multi-Layer Multi-Configuration Time-Dependent Hartree (ML-MCTDH) method for a system of up to 63 electronic states and 180 vibrational modes. A Frenkel Hamiltonian of HJ aggregate type is used, along with a reduced H-aggregate representation; electron-phonon coupling includes local high-frequency modes as well as anharmonic intermolecular modes. The latter are operative in mediating inter-chain transport, by a mechanism of transient localization type. Strikingly, this mechanism is found to be of quantum coherent character and involves non-adiabatic effects. Using periodic boundary conditions, a normal diffusion regime is identified from the exciton mean-squared displacement, apart from early-time transients. Diffusion coefficients are found to be of the order of 3 x 10-3 cm2/s, showing a non-monotonous increase with temperature.

Role of Spatially Correlated Fluctuations in Photosynthetic Excitation Energy Transfer with an Equilibrium and a Nonequilibrium Initial Bath

The transfer of excitation energy in photosynthetic light-harvesting complexes has inspired growing interest for its scientific and engineering significance. Recent experimental findings have suggested that spatially correlated environmental fluctuations may account for the existence of long-lived quantum coherent energy transfer observed even at physiological temperature. In this paper, we investigate the effects of spatial correlations on the excitation energy transfer dynamics by including a nonequilibrium initial bath in a simulated donor-acceptor model. The initial bath state, which is assumed to be either equilibrium or nonequilibrium, is expanded in powers of coupling strength within the polaron formalism of a quantum master equation. The spatial correlations of bath fluctuations strongly influence the decay of coherence in the dynamics. The role of a nonequilibrium initial bath is also influenced by spatial correlations and becomes the most conspicuous for certain degrees of spatial correlations from which we propose a picture that the spatial correlations of bath fluctuations open up new energy transfer pathways, playing a role of protecting coherence. Besides, we apply the polaron master equation approach to study the dynamics in a two-site subsystem of the FMO complex and provide a practical example that shows the versatility of this approach.

Excitation Energy Transfer Efficiency in Fluctuating Environments: Role of Quantum Coherence in the Presence of Memory Effects

Several recent studies have interrogated the role of quantum coherence in affecting the transfer efficiency of an optical excitation to the designated "trap" state where the energy can be used for subsequent reactions, as in photosynthesis. However, these studies invoke a Markovian approximation for the time correlation function describing the environment-induced stochastic fluctuations. Here, we employ Kubo's quantum stochastic Liouville equation (QSLE) to include memory effects. We extend the existing QSLE scheme to introduce decay of a newly created excitation due to radiative and nonradiative channels and also by desired trapping toward the targeted chromophore. We show that the theoretical formalism based on the QSLE correctly reproduces the rate equation description in the Markovian limit, with the rate constants determined by an appropriate quantum limiting procedure. We find that under certain conditions, the efficiency of excitation transfer to the trap gains from the combined presence of quantum coherence and temporally correlated stochastic fluctuations. We work out different limiting situations in order to discover and quantify the optimum conditions for the energy transfer to the trapped state. We find that maximum energy transfer efficiency is achieved in the intermediate limit between coherent and incoherent transport.

An exact solution in the theory of fluorescence resonance energy transfer with vibrational relaxation

The elegant expression of Förster that predicts the well-known 1/R⁶ distance (R) dependence of the rate of energy transfer, although widely used, was derived using several approximations. Notable among them is the neglect of the vibrational relaxation in the reactant (donor) and product (acceptor) manifolds. Vibrational relaxation can play an important role when the energy transfer rate is faster than the vibrational relaxation rate. Under such conditions, donor to acceptor energy transfer can occur from the excited vibrational states. This phenomenon is not captured by the usual formulation based on the overlap of donor emission and acceptor absorption spectra. Here, we develop a Green's function-based generalized formalism and obtain an exact solution for the excited state population relaxation and the rate of energy transfer in the presence of vibrational relaxation. We find that the application of the well-known Förster's expression might lead to overestimation of R.

Coherent Speedup of Excitation Energy Transfer in PC645

The fundamental principle governing the robust excitation energy transfer (EET) in the light-harvesting complexes remained elusive. Recent spectroscopic observations of the EET attract attention toward the role of the quantum effects in the biological systems. In this study, along with the time evolution, the distance between the quantum states of the pigments is also considered for a detailed illustration of the effect of the quantum coherence on the speed of the excitation transfer. The comparative analysis of the coherently delocalized EET with the incoherent discrete hopping reveals that the coherent superposition speedup the excitation transfer process, for instance, with the initially localized excitation, coherence provides a ∼4-fold enhancement to the excitation transfer rate.

Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature

We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.

Polaritons and excitons: Hamiltonian design for enhanced coherence

The primary questions motivating this report are: Are there ways to increase coherence and delocalization of excitation among many molecules at moderate electronic coupling strength? Coherent delocalization of excitation in disordered molecular systems is studied using numerical calculations. The results are relevant to molecular excitons, polaritons, and make connections to classical phase oscillator synchronization. In particular, it is hypothesized that it is not only the magnitude of electronic coupling relative to the standard deviation of energetic disorder that decides the limits of coherence, but that the structure of the Hamiltonian-connections between sites (or molecules) made by electronic coupling-is a significant design parameter. Inspired by synchronization phenomena in analogous systems of phase oscillators, some properties of graphs that define the structure of different Hamiltonian matrices are explored. The report focuses on eigenvalues and ensemble density matrices of various structured, random matrices. Some reasons for the special delocalization properties and robustness of polaritons in the single-excitation subspace (the star graph) are discussed. The key result of this report is that, for some classes of Hamiltonian matrix structure, coherent delocalization is not easily defeated by energy disorder, even when the electronic coupling is small compared to disorder.

First-Principles Quantum and Quantum-Classical Simulations of Exciton Diffusion in Semiconducting Polymer Chains at Finite Temperature

We report on first-principles quantum-dynamical and quantum-classical simulations of photoinduced exciton dynamics in oligothiophene chain segments, representative of intrachain exciton migration in the poly(3-hexylthiophene) (P3HT) polymer. Following up on our recent study (Binder R.; Burghardt, I. Faraday Discuss. 2020, 221, 406), multilayer multiconfiguration time-dependent Hartree calculations for a short oligothiophene segment comprising 20 monomer units (OT-20) are carried out to obtain full quantum-dynamical simulations at finite temperature. These are employed to benchmark mean-field Ehrenfest calculations, which are shown to give qualitatively correct results for the present system. Periodic boundary conditions turn out to significantly improve earlier estimates of diffusion coefficients. Using the Ehrenfest approach, a series of calculations are subsequently carried out for larger lattices (OT-40 to OT-80), leading to estimates for temperature-dependent mean-squared displacements, which are found to exhibit a near-linear dependence as a function of time. The resulting diffusion coefficient estimates are an increasing function of temperature, whose detailed functional form depends on the degree of static disorder. With a realistic static disorder parameter (σ_s ≃ 0.06 eV), the diffusion coefficients decrease from D ∼ 1 × 10^-2 cm² s^-1 to D ∼ 1 × 10^-3 cm² s^-1, in qualitative agreement with experimental data for P3HT. The dynamical scenario obtained from our simulations shows that exciton migration in P3HT-type chains is a largely adiabatic process throughout the temperature regime we investigated (i.e., T = 50-300 K). The resulting picture of exciton migration is a coherent, but not bandlike, motion of an exciton-polaron driven by fluctuations induced by low-frequency modes. This process acquires partial hopping character if static disorder becomes prominent and Anderson localization sets in.

First-principles description of intra-chain exciton migration in an oligo(para-phenylene vinylene) chain. I. Generalized Frenkel-Holstein Hamiltonian

A generalized Frenkel-Holstein Hamiltonian is constructed to describe exciton migration in oligo(para-phenylene vinylene) chains, based on excited state electronic structure data for an oligomer comprising 20 monomer units (OPV-20). Time-dependent density functional theory calculations using the ωB97XD hybrid functional are employed in conjunction with a transition density analysis to study the low-lying singlet excitations and demonstrate that these can be characterized to a good approximation as a Frenkel exciton manifold. Based on these findings, we employ the analytic mapping procedure of Binder et al. [J. Chem. Phys. 141, 014101 (2014)] to translate one-dimensional (1D) and two-dimensional (2D) potential energy surface (PES) scans to a fully anharmonic, generalized Frenkel-Holstein (FH) Hamiltonian. A 1D PES scan is carried out for intra-ring quinoid distortion modes, while 2D PES scans are performed for the anharmonically coupled inter-monomer torsional and vinylene bridge bond length alternation modes. The kinetic energy is constructed in curvilinear coordinates by an exact numerical procedure, using the TNUM Fortran code. As a result, a fully molecular-based, generalized FH Hamiltonian is obtained, which is subsequently employed for quantum exciton dynamics simulations, as shown in Paper II [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)].

Importance of Intermonomer Couplings of the FMO Complex in Coherently Initiated Dynamics

In this work, the hierarchical equation of motion approach was employed to study the excitation transfer efficiency (ETE) from the Fenna-Matthews-Olson (FMO) complex to the reaction center core complex (RCC). Further, the realistic initial pure state (created by the interaction of femtosecond laser pulse with the isolated FMO complex) was used in the presence of an inhomogeneous protein environment. We observed the same amount of excitation transfer from the FMO to the RCC, in agreement with the experiment using femtosecond absorption spectroscopy. Simulation analysis shows weak intermonomer couplings are strongly affecting the dynamics of ETE, which implies that all the monomers of the FMO complex need to be considered as connected with the RCC pigment for a better estimate of the ETE.

First-principles quantum simulations of exciton diffusion on a minimal oligothiophene chain at finite temperature

High-dimensional multiconfigurational quantum dynamics simulations are carried out at finite temperature to simulate exciton diffusion on an oligothiophene chain, representative of a segment of the poly(3-hexylthiophene) (P3HT) polymer. The ab initio parametrized site-based Hamiltonian of Binder et al. [Phys. Rev. Lett., 2018, 120, 227401] is employed to model a 20-site system, including intra-ring and inter-ring high-frequency modes as well as torsional modes which undergo thermal fluctuations induced by an explicit harmonic oscillator bath. The system-bath dynamics is treated within the setting of a stochastic mean-field Schrödinger equation. For the 20-site excitonic system, a total of 20 Frenkel states and 248 modes are propagated using the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method. The resulting dynamics can be interpreted in terms of the coherent motion of an exciton-polaron quasi-particle stochastically driven by torsional fluctuations. This dynamics yields a near-linear mean squared displacement (MSD) as a function of time, from which a diffusion coefficient can be deduced which increases with temperature, up to 5.7 × 10-3 cm2 s-1 at T = 300 K.

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.

Delocalization and Quantum Entanglement in Physical Systems

Quantum coherence and entanglement in an extended interacting system where energy levels are nondegenerate and coupled to a dissipative environment is a common occurrence in nature, like in photosynthetic reaction systems and conjugated polymers. The temperature dependence of quantum coherence in a trimer complex (first three subunits of the Fenna-Matthews-Olson complex) is studied using a temperature-dependent quantum stochastic Liouville equation. In the non-Markovian limit, the lowering of temperature induces long-lasting quantum coherence that, in turn, leads to delocalization, whose length grows. The entanglement and coherence length determine the nature of the dynamic localization.

Role of quantum coherence in shaping the line shape of an exciton interacting with a spatially and temporally correlated bath

Kubo's fluctuation theory of line shape forms the backbone of our understanding of optical and vibrational line shapes, through such concepts as static heterogeneity and motional narrowing. However, the theory does not properly address the effects of quantum coherences on optical line shape, especially in extended systems where a large number of eigenstates are present. In this work, we study the line shape of an exciton in a one-dimensional lattice consisting of regularly placed and equally separated optical two level systems. We consider both linear array and cyclic ring systems of different sizes. Detailed analytical calculations of line shape have been carried out by using Kubo's stochastic Liouville equation (SLE). We make use of the observation that in the site representation, the Hamiltonian of our system with constant off-diagonal coupling J is a tridiagonal Toeplitz matrix (TDTM) whose eigenvalues and eigenfunctions are known analytically. This identification is particularly useful for long chains where the eigenvalues of TDTM help understanding crossover between static and fast modulation limits. We summarize the new results as follows. (i) In the slow modulation limit when the bath correlation time is large, the effects of spatial correlation are not negligible. Here the line shape is broadened and the number of peaks increases beyond the ones obtained from TDTM (constant off-diagonal coupling element J and no fluctuation). (ii) However, in the fast modulation limit when the bath correlation time is small, the spatial correlation is less important. In this limit, the line shape shows motional narrowing with peaks at the values predicted by TDTM (constant J and no fluctuation). (iii) Importantly, we find that the line shape can capture that quantum coherence affects in the two limits differently. (iv) In addition to linear chains of two level systems, we also consider a cyclic tetramer. The cyclic polymers can be designed for experimental verification. (v) We also build a connection between line shape and population transfer dynamics. In the fast modulation limit, both the line shape and the population relaxation, for both correlated and uncorrelated bath, show similar behavior. However, in slow modulation limit, they show profoundly different behavior. (vi) This study explains the unique role of the rate of fluctuation (inverse of the bath correlation time) in the sustenance and propagation of coherence. We also examine the effects of off-diagonal fluctuation in spectral line shape. Finally, we use Tanimura-Kubo formalism to derive a set of coupled equations to include temperature effects (partly neglected in the SLE employed here) and effects of vibrational mode in energy transfer dynamics.

Environment-Assisted Quantum Coherence in Photosynthetic Complex

Recent experiments [ Engel et al. Nature, 2007, 446, 782-786 ] revealed the existence of surprisingly long-lived quantum coherence in the noisy biological environment of the photosynthetic Fenna-Matthews-Olson (FMO) complex. Such coherence can clearly play an important role in facilitating efficient energy transfer. The occurrence of quantum coherence in quantum transport is also implicated in excitation transport processes in conjugated polymers [ Collini et al. Science, 2009, 323, 369-373 ]. Even though these systems are strongly correlated, most theoretical studies invoke a Markovian approximation where the temporal correlation of bath fluctuations is neglected. We use an elegant nonperturbative method based on Kubo's quantum stochastic Liouville equation (QSLE) to study the effects of correlated non-Markovian bath fluctuations in several different limits and find the interesting result that fluctuations not only destroy coherence but under appropriate conditions can also facilitate it. We show that temperature has the most pronounced effect in the intermediate coupling limit where it can promote transition from coherent to incoherent transfer.