Electronic excitation dynamics in multichromophoric systems described via a polaron-representation master equation

Capturing non-Markovian polaron dressing with the master equation formalism

Understanding the dynamics of open quantum systems in strong coupling and non-Markovian regimes remains a formidable theoretical challenge. One popular and well-established method of approximation in these circumstances is provided by the polaron master equation (PME). In this work, we re-evaluate and extend the validity of the PME to capture the impact of non-Markovian polaron dressing, induced by non-equilibrium open system dynamics. By comparing with numerically exact techniques, we confirm that while the standard PME successfully predicts the dynamics of system observables that commute with the polaron transformation (e.g., populations in the Pauli z-basis), it can struggle to fully capture those that do not (e.g., coherences). This limitation stems from the mixing of system and environment degrees of freedom inherent to the polaron transformation, which affects the accuracy of calculated expectation values within the polaron frame. Employing the Nakajima-Zwanzig projection operator formalism, we introduce correction terms that provide an accurate description of observables that do not commute with the transformation. We demonstrate the significance of the correction terms in two cases, the canonical spin-boson model and a dissipative time-dependent Landau-Zener protocol, where they are shown to impact the system dynamics on both short and long timescales.

Controlling Vibronic Coupling in Chlorophyll Proteins: The Effects of Excitonic Delocalization and Vibrational Localization

Vibrational-electronic (vibronic) coupling plays a critical role in excitation energy transfer in molecular aggregates and pigment-protein complexes (PPCs). But the interplay between excitonic delocalization and vibronic interactions is complex, often leaving even qualitative questions as to what conceptual framework (e.g., Redfield versus Förster theory) should be used to interpret experimental results. To shed light on this issue, we report here on the interplay between excitonic delocalization and vibronic coupling in site-directed mutants of the water-soluble chlorophyll protein (WSCP), as reflected in 77 K fluorescence spectra. Experimentally, we find that in PPCs where excitonic delocalization is disrupted (either by mutagenesis or heterodimer formation), the relative intensity of the vibrational sideband (VSB) in fluorescence spectra is suppressed by up to 37% compared to that of the native protein. Numerical simulations reveal that this effect results from the localization of high-frequency vibrations in the coupled system; while excitonic delocalization suppresses the purely electronic transition due to H-aggregate-like dipole-dipole interference, high-frequency vibrations are unaffected, leading to a relative enhancement of the VSB. By comparing VSB intensities of PPCs both in the presence and absence of excitonic delocalization, we extract a set of "local" Huang-Rhys (HR) factors for Chl a in WSCP. More generally, our results suggest a significant role for geometric effects in controlling energy-transfer rates (which depend sensitively on absorption/fluorescence line shapes) in molecular aggregates and PPCs.

Quantum Work Statistics at Strong Reservoir Coupling

Determining the statistics of work done on a quantum system while strongly coupled to a reservoir is a formidable task, requiring the calculation of the full eigenspectrum of the combined system and reservoir. Here, we show that this issue can be circumvented by using a polaron transformation that maps the system into a new frame where weak-coupling theory can be applied. Crucially, this polaron approach reproduces the Jarzynski fluctuation theorem, thus ensuring consistency with the laws of stochastic thermodynamics. We apply our formalism to a system driven across the Landau-Zener transition, where we identify clear signatures in the work distribution arising from a non-negligible coupling to the environment. Our results provide a new method for studying the stochastic thermodynamics of driven quantum systems beyond Markovian, weak-coupling regimes.

Tunability in 3-Site Electronic Excitation Transfer Dynamics: Insights into the Role of Perturbative Coupling

Electronic excitation transfer dynamics in photosynthetic systems, including the Fenna-Matthews-Olson complex, are often modeled using the interaction picture of three two-level systems, also known as the 3-site system. Among the two possible configurations, uphill and downhill, a recent publication reported an intriguing correlation between population dynamics and the intersite coupling. Specifically, the uphill configuration has been shown to have a pronounced dependence on perturbations in the intersite coupling, whereas the downhill configuration displays negligible dependence. In this study, we consider a generic 3-site configuration and model site interactions through the Markovian master equation. Through this approach, we derive succinct analytical expressions for the population dynamics between the sites, shedding light on the differences in behavior between the two configurations. Using these analytical expressions, we demonstrate the range of tunability in population dynamics achievable with minimal changes in intersite coupling, and we validate these findings through simulation results. These insights into the population dynamics of a 3-site system are expected to play a crucial role in facilitating the design of efficient energy-transfer pathways in molecular systems.

Vibronic effects on the quantum tunnelling of magnetisation in Kramers single-molecule magnets

Single-molecule magnets are among the most promising platforms for achieving molecular-scale data storage and processing. Their magnetisation dynamics are determined by the interplay between electronic and vibrational degrees of freedom, which can couple coherently, leading to complex vibronic dynamics. Building on an ab initio description of the electronic and vibrational Hamiltonians, we formulate a non-perturbative vibronic model of the low-energy magnetic degrees of freedom in monometallic single-molecule magnets. Describing their low-temperature magnetism in terms of magnetic polarons, we are able to quantify the vibronic contribution to the quantum tunnelling of the magnetisation, a process that is commonly assumed to be independent of spin-phonon coupling. We find that the formation of magnetic polarons lowers the tunnelling probability in both amorphous and crystalline systems by stabilising the low-lying spin states. This work, thus, shows that spin-phonon coupling subtly influences magnetic relaxation in single-molecule magnets even at extremely low temperatures where no vibrational excitations are present.

Resonant vibrations produce quantum bridge over high-energy states in heterogeneous antenna

Photosynthetic light-harvesting complexes usually contain several pools of molecules with a big difference in transition energies, for example, chlorophylls a and b in plant antennas. Some pathways of the excitation energy transfer may include pigments from the low-energy pool separated by a site occupied by a high-energy molecule. We demonstrate that such pathways may be functional if high-frequency intramolecular vibrations fall in resonance with the energy gap between the neighboring molecules belonging to different pools. In this case, a vibration-assisted mixing of the excited states can produce delocalized vibronic states playing a role of 'quantum bridge' that facilitates a passage over high-energy barrier. We perform calculations of the excitation dynamics in the model three-state system with the parameters emerging from our previous studies of real antennas. Simulation of the dynamics in an explicit electron-vibrational basis demonstrates that the rate of transfer between the two chlorophylls a through the chlorophyll b intermediate is increased by a factor of 1.7-2 in the presence of resonant vibration. A possible influence of energetic disorder and other (non-resonant) vibrations on this effect is discussed.

Mechanism of Delocalization-Enhanced Exciton Transport in Disordered Organic Semiconductors

Large exciton diffusion lengths generally improve the performance of organic semiconductor devices, because they enable energy to be transported farther during the exciton lifetime. However, the physics of exciton motion in disordered organic materials is not fully understood, and modeling the transport of quantum-mechanically delocalized excitons in disordered organic semiconductors is a computational challenge. Here, we describe delocalized kinetic Monte Carlo (dKMC), the first model of three-dimensional exciton transport in organic semiconductors that includes delocalization, disorder, and polaron formation. We find that delocalization can dramatically increase exciton transport; for example, delocalization across less than two molecules in each direction can increase the exciton diffusion coefficient by over an order of magnitude. The mechanism for the enhancement is 2-fold: delocalization enables excitons to hop both more frequently and further in each hop. We also quantify the effect of transient delocalization (short-lived periods where excitons become highly delocalized) and show that it depends strongly upon the disorder and transition dipole moments.

Compact and complete description of non-Markovian dynamics

Generalized master equations provide a theoretically rigorous framework to capture the dynamics of processes ranging from energy harvesting in plants and photovoltaic devices to qubit decoherence in quantum technologies and even protein folding. At their center is the concept of memory. The explicit time-nonlocal description of memory is both protracted and elaborate. When physical intuition is at a premium, one would desire a more compact, yet complete, description. Here, we demonstrate how and when the time-convolutionless formalism constitutes such a description. In particular, by focusing on the dissipative dynamics of the spin-boson and Frenkel exciton models, we show how to: easily construct the time-local generator from reference reduced dynamics, elucidate the dependence of its existence on the system parameters and the choice of reduced observables, identify the physical origin of its apparent divergences, and offer analysis tools to diagnose their severity and circumvent their deleterious effects. We demonstrate that, when applicable, the time-local approach requires as little information as the more commonly used time-nonlocal scheme, with the important advantages of providing a more compact description, greater algorithmic simplicity, and physical interpretability. We conclude by introducing the discrete-time analog and a straightforward protocol to employ it in cases where the reference dynamics have limited resolution. The insights we present here offer the potential for extending the reach of dynamical methods, reducing both their cost and conceptual complexity.

Even a little delocalization produces large kinetic enhancements of charge-separation efficiency in organic photovoltaics

In organic photovoltaics, charges can separate efficiently even if their Coulomb attraction is an order of magnitude greater than the available thermal energy. Delocalization has been suggested to explain this fact, because it could increase the initial separation of charges in the charge-transfer (CT) state, reducing their attraction. However, understanding the mechanism requires a kinetic model of delocalized charge separation, which has proven difficult because it involves tracking the correlated quantum-mechanical motion of the electron and the hole in large simulation boxes required for disordered materials. Here, we report the first three-dimensional simulations of charge-separation dynamics in the presence of disorder, delocalization, and polaron formation, finding that even slight delocalization, across less than two molecules, can substantially enhance the charge-separation efficiency, even starting with thermalized CT states. Delocalization does not enhance efficiency by reducing the Coulomb attraction; instead, the enhancement is a kinetic effect produced by the increased overlap of electronic states.

On simulating the dynamics of electronic populations and coherences via quantum master equations based on treating off-diagonal electronic coupling terms as a small perturbation

Quantum master equations provide a general framework for describing the dynamics of electronic observables within a complex molecular system. One particular family of such equations is based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. In this paper, we show how different choices of projection operators, as well as whether one starts out with the time-convolution or the time-convolutionless forms of the generalized quantum master equation, give rise to four different types of such off-diagonal quantum master equations (OD-QMEs), namely, time-convolution and time-convolutionless versions of a Pauli-type OD-QME for only the electronic populations and an OD-QME for the full electronic density matrix (including both electronic populations and coherences). The fact that those OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. To address this, we also extend a procedure for extracting Schrödinger picture electronic coherences from interaction picture populations recently introduced by Trushechkin in the context of time-convolutionless Pauli-type OD-QME to the other three types of OD-QMEs. The performance of the aforementioned four types of OD-QMEs is explored in the context of the Garg-Onuchic-Ambegaokar benchmark model for charge transfer in the condensed phase across a relatively wide parameter range. The results show that time-convolution OD-QMEs can be significantly more accurate than their time-convolutionless counterparts, particularly in the case of Pauli-type OD-QMEs, and that rather accurate Schrödinger picture coherences can be obtained from interaction picture electronic inputs.

Strong Exciton-Vibrational Coupling in Molecular Assemblies. Dynamics Using the Polaron Transformation in HEOM Space

In Frenkel exciton dynamics of aggregated molecules, the polaron transformation (PT) technique leads to decoupling of diagonal elements in the subspace of excited electronic states from vibrations. In this article we describe for the first time how PT becomes applicable in the framework of the "Hierarchical Equations of Motion" (HEOM) approach for treatment of open quantum systems. We extend the concept of formulating operators in HEOM space by deriving hierarchical equations of PT which lead to a shift in the excited state potential energy surface to compensate its displacement. While the assumption of thermal equilibration of the vibrational oscillators, introduced by PT, results in a stationary state in a monomer, in a dimer under the same assumption nonequilibrium dynamics appears because of the interplay of the transfer process and vibrational equilibration. Both vertical transitions generating a vibrationally hot state and initially equilibrated vibrational oscillators evolve toward the same stationary asymptotic state associated with polaron formation. The effect of PT on the dynamics of this process depends on initial excitation and basis representation of the electronic system. The developed approach facilitates a generic formulation of quantum master equations involving perturbative treatment of polaron dynamics.

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.

Analytical derivation of equilibrium state for open quantum system

Calculation of the equilibrium state of an open quantum system interacting with a bath remains a challenge to this day, mostly due to a huge number of bath degrees of freedom. Here, we present an analytical expression for the reduced density operator in terms of an effective Hamiltonian for a high temperature case. Comparing with numerically exact results, we show that our theory is accurate for slow baths and up to intermediate system-bath coupling strengths. Our results demonstrate that the equilibrium state does not depend on the shape of spectral density in the slow bath regime. The key quantity in our theory is the effective coupling between the states, which depends exponentially on the ratio of the reorganization energy to temperature and, thus, has opposite temperature dependence than could be expected from the small polaron transformation.

Improved population operators for multi-state nonadiabatic dynamics with the mixed quantum-classical mapping approach

Application to the 7-state Frenkel-exciton Hamiltonian for the Fenna–Matthews–Olson complex shows that using a different representation of the electronic population operators can drastically improve the accuracy of the quasiclassical mapping approach without increasing the computational effort.

Linear and nonlinear spectroscopy from quantum master equations

We investigate the accuracy of the second-order time-convolutionless (TCL2) quantum master equation for the calculation of linear and nonlinear spectroscopies of multichromophore systems. We show that even for systems with non-adiabatic coupling, the TCL2 master equation predicts linear absorption spectra that are accurate over an extremely broad range of parameters and well beyond what would be expected based on the perturbative nature of the approach; non-equilibrium population dynamics calculated with TCL2 for identical parameters are significantly less accurate. For third-order (two-dimensional) spectroscopy, the importance of population dynamics and the violation of the so-called quantum regression theorem degrade the accuracy of TCL2 dynamics. To correct these failures, we combine the TCL2 approach with a classical ensemble sampling of slow microscopic bath degrees of freedom, leading to an efficient hybrid quantum-classical scheme that displays excellent accuracy over a wide range of parameters. In the spectroscopic setting, the success of such a hybrid scheme can be understood through its separate treatment of homogeneous and inhomogeneous broadening. Importantly, the presented approach has the computational scaling of TCL2, with the modest addition of an embarrassingly parallel prefactor associated with ensemble sampling. The presented approach can be understood as a generalized inhomogeneous cumulant expansion technique, capable of treating multilevel systems with non-adiabatic dynamics.

Applicability of transfer tensor method for open quantum system dynamics

Accurate simulations of open quantum system dynamics is a long standing issue in the field of chemical physics. Exact methods exist, but are costly, while perturbative methods are limited in their applicability. Recently a new black-box type method, called transfer tensor method (TTM), was proposed [J. Cerrillo and J. Cao, Phys. Rev. Lett. 112, 110401 (2014)]. It allows one to accurately simulate long time dynamics with a numerical cost of solving a time-convolution master equation, provided many initial system evolution trajectories are obtained from some exact method beforehand. The possible time-savings thus strongly depend on the ratio of total versus initial evolution lengths. In this work, we investigate the parameter regimes where an application of TTM would be most beneficial in terms of computational time. We identify several promising parameter regimes. Although some of them correspond to cases when perturbative theories could be expected to perform well, we find that the accuracy of such approaches depends on system parameters in a more complex way than it is commonly thought. We propose that the TTM should be applied whenever system evolution is expected to be long and accuracy of perturbative methods cannot be ensured or in cases when the system under consideration does not correspond to any single perturbative regime.

Quantum design of photosynthesis for bio-inspired solar-energy conversion

Photosynthesis is the natural process that converts solar photons into energy-rich products that are needed to drive the biochemistry of life. Two ultrafast processes form the basis of photosynthesis: excitation energy transfer and charge separation. Under optimal conditions, every photon that is absorbed is used by the photosynthetic organism. Fundamental quantum mechanics phenomena, including delocalization, underlie the speed, efficiency and directionality of the charge-separation process. At least four design principles are active in natural photosynthesis, and these can be applied practically to stimulate the development of bio-inspired, human-made energy conversion systems.

Spatio-Temporal Steering for Testing Nonclassical Correlations in Quantum Networks

We introduce the concept of spatio-temporal steering (STS), which reduces, in special cases, to Einstein-Podolsky-Rosen steering and the recently-introduced temporal steering. We describe two measures of this effect referred to as the STS weight and robustness. We suggest that these STS measures enable a new way to assess nonclassical correlations in an open quantum network, such as quantum transport through nano-structures or excitation transfer in a complex biological system. As one of our examples, we apply STS to check nonclassical correlations among sites in a photosynthetic pigment-protein complex in the Fenna-Matthews-Olson model.

Quantum - coherent dynamics in photosynthetic charge separation revealed by wavelet analysis

Experimental/theoretical evidence for sustained vibration-assisted electronic (vibronic) coherence in the Photosystem II Reaction Center (PSII RC) indicates that photosynthetic solar-energy conversion might be optimized through the interplay of electronic and vibrational quantum dynamics. This evidence has been obtained by investigating the primary charge separation process in the PSII RC by two-dimensional electronic spectroscopy (2DES) and Redfield modeling of the experimental data. However, while conventional Fourier transform analysis of the 2DES data allows oscillatory signatures of vibronic coherence to be identified in the frequency domain in the form of static 2D frequency maps, the real-time evolution of the coherences is lost. Here we apply for the first time wavelet analysis to the PSII RC 2DES data to obtain time-resolved 2D frequency maps. These maps allow us to demonstrate that (i) coherence between the excitons initiating the two different charge separation pathways is active for more than 500 fs, and (ii) coherence between exciton and charge-transfer states, the reactant and product of the charge separation reaction, respectively; is active for at least 1 ps. These findings imply that the PSII RC employs coherence (i) to sample competing electron transfer pathways, and ii) to perform directed, ultrafast and efficient electron transfer.

Origin of long-lived quantum coherence and excitation dynamics in pigment-protein complexes

We explore the mechanism for the long-lived quantum coherence by considering the discrete phonon modes: these vibrational modes effectively weaken the exciton-environment interaction, due to the new composite (polaron) formed by excitons and vibrons. This subsequently demonstrates the role of vibrational coherence which greatly contributes to long-lived feature of the excitonic coherence that has been observed in femtosecond experiments. The estimation of the timescale of coherence elongated by vibrational modes is given in an analytical manner. To test the validity of our theory, we study the pigment-protein complex in detail by exploring the energy transfer and coherence dynamics. The ground-state vibrational coherence generated by incoherent radiations is shown to be long-survived and is demonstrated to be significant in promoting the excitation energy transfer. This is attributed to the nonequilibriumness of the system caused by the detailed-balance-breaking, which funnels the downhill migration of excitons.

How exciton-vibrational coherences control charge separation in the photosystem II reaction center

In photosynthesis absorbed sun light produces collective excitations (excitons) that form a coherent superposition of electronic and vibrational states of the individual pigments. Two-dimensional (2D) electronic spectroscopy allows a visualization of how these coherences are involved in the primary processes of energy and charge transfer. Based on quantitative modeling we identify the exciton-vibrational coherences observed in 2D photon echo of the photosystem II reaction center (PSII-RC). We find that the vibrations resonant with the exciton splittings can modify the delocalization of the exciton states and produce additional states, thus promoting directed energy transfer and allowing a switch between the two charge separation pathways. We conclude that the coincidence of the frequencies of the most intense vibrations with the splittings within the manifold of exciton and charge-transfer states in the PSII-RC is not occurring by chance, but reflects a fundamental principle of how energy conversion in photosynthesis was optimized.

The Role of Resonant Vibrations in Electronic Energy Transfer

Nuclear vibrations play a prominent role in the spectroscopy and dynamics of electronic systems. As recent experimental and theoretical studies suggest, this may be even more so when vibrational frequencies are resonant with transitions between the electronic states. Herein, a vibronic multilevel Redfield model is reported for excitonically coupled electronic two-level systems with a few explicitly included vibrational modes and interacting with a phonon bath. With numerical simulations the effects of the quantized vibrations on the dynamics of energy transfer and coherence in a model dimer are illustrated. The resonance between the vibrational frequency and energy gap between the sites leads to a large delocalization of vibronic states, which then results in faster energy transfer and longer-lived mixed coherences.

Modelling exciton-phonon interactions in optically driven quantum dots

We provide a self-contained review of master equation approaches to modelling phonon effects in optically driven self-assembled quantum dots. Coupling of the (quasi) two-level excitonic system to phonons leads to dissipation and dephasing, the rates of which depend on the excitation conditions, intrinsic properties of the QD sample, and its temperature. We describe several techniques, which include weak-coupling master equations that are perturbative in the exciton-phonon coupling, as well as those based on the polaron transformation that can remain valid for strong phonon interactions. We additionally consider the role of phonons in altering the optical emission characteristics of quantum dot devices, outlining how we must modify standard quantum optics treatments to account for the presence of the solid-state environment.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

Exact vs. asymptotic spectral densities in the Garg-Onuchic-Ambegaokar charge transfer model and its effect on Fermi's golden rule rate constants

The Garg-Onuchic-Ambegaokar model [J. Chem. Phys. 83, 4491 (1985)] has been used extensively for benchmarking methods aimed at calculating charge transfer rates. Within this model, the donor and acceptor diabats are described as shifted parabolas along a single primary mode, which is bilinearly coupled to a harmonic bath consisting of secondary modes, characterized by an Ohmic spectral density with exponential cutoff. Rate calculations for this model are often performed in the normal mode representation, with the corresponding effective spectral density given by an asymptotic expression derived at the limit where the Ohmic bath cutoff frequency is much larger than the primary mode frequency. We compare Fermi's golden rule rate constants obtained with the asymptotic and exact effective spectral densities. We find significant deviations between rate constants obtained from the asymptotic spectral density and those obtained from the exact one in the deep inverted region. Within the range of primary mode frequencies commonly employed, we find that the discrepancies increase with decreasing temperature and with decreasing primary mode frequency.

Energy transfer efficiency in the chromophore network strongly coupled to a vibrational mode

Using methods from condensed matter and statistical physics, we examine the transport of excitons through the photosynthetic complex from a receiving antenna to a reaction center. Writing the equations of motion for the exciton creation-annihilation operators, we are able to describe the exciton dynamics, even in the regime when the reorganization energy is of the order of the intrasystem couplings. We determine the exciton transfer efficiency in the presence of a quenching field and protein environment. While the majority of the protein vibrational modes are treated as a heat bath, we address the situation when specific modes are strongly coupled to excitons and examine the effects of these modes on the energy transfer efficiency in the steady-state regime. Using the structural parameters of the Fenna-Matthews-Olson complex, we find that, for vibrational frequencies below 16 meV, the exciton transfer is drastically suppressed. We attribute this effect to the formation of a "mixed exciton-vibrational mode" where the exciton is transferred back and forth between the two pigments with the absorption or emission of vibrational quanta, instead of proceeding to the reaction center. The same effect suppresses the quantum beating at the vibrational frequency of 25 meV. We also show that the efficiency of the energy transfer can be enhanced when the vibrational mode strongly couples to the third pigment only, instead of coupling to the entire system.

On the accuracy of coherent modified Redfield theory in simulating excitation energy transfer dynamics

In this study, we investigate the accuracy of a recently developed coherent modified Redfield theory (CMRT) in simulating excitation energy transfer (EET) dynamics. The CMRT is a secular non-Markovian quantum master equation that is derived by extending the modified Redfield theory to treat coherence dynamics in molecular excitonic systems. Herein, we systematically survey the applicability of the CMRT in a large EET parameter space through the comparisons of the CMRT EET dynamics in a dimer system with the numerically exact results. The results confirm that the CMRT exhibits a broad applicable range and allow us to locate the specific parameter regimes where CMRT fails to provide adequate results. Moreover, we propose an accuracy criterion based on the magnitude of second-order perturbation to characterize the applicability of CMRT and show that the criterion summarizes all the benchmark results and the physics described by CMRT. Finally, we employ the accuracy criterion to quantitatively compare the performance of CMRT to that of a small polaron quantum master equation approach. The comparison demonstrates the complementary nature of these two methods, and as a result, the combination of the two methods provides accurate simulations of EET dynamics for the full parameter space investigated in this study. Our results not only delicately evaluate the applicability of the CMRT but also reveal new physical insights for factors controlling the dynamics of EET that should be useful for developing more accurate and efficient methods for simulations of EET dynamics in molecular aggregate systems.

Dissipation enhanced vibrational sensing in an olfactory molecular switch

Motivated by a proposed olfactory mechanism based on a vibrationally activated molecular switch, we study electron transport within a donor-acceptor pair that is coupled to a vibrational mode and embedded in a surrounding environment. We derive a polaron master equation with which we study the dynamics of both the electronic and vibrational degrees of freedom beyond previously employed semiclassical (Marcus-Jortner) rate analyses. We show (i) that in the absence of explicit dissipation of the vibrational mode, the semiclassical approach is generally unable to capture the dynamics predicted by our master equation due to both its assumption of one-way (exponential) electron transfer from donor to acceptor and its neglect of the spectral details of the environment; (ii) that by additionally allowing strong dissipation to act on the odorant vibrational mode, we can recover exponential electron transfer, though typically at a rate that differs from that given by the Marcus-Jortner expression; (iii) that the ability of the molecular switch to discriminate between the presence and absence of the odorant, and its sensitivity to the odorant vibrational frequency, is enhanced significantly in this strong dissipation regime, when compared to the case without mode dissipation; and (iv) that details of the environment absent from previous Marcus-Jortner analyses can also dramatically alter the sensitivity of the molecular switch, in particular, allowing its frequency resolution to be improved. Our results thus demonstrate the constructive role dissipation can play in facilitating sensitive and selective operation in molecular switch devices, as well as the inadequacy of semiclassical rate equations in analysing such behaviour over a wide range of parameters.

Optimal Energy Transfer in Light-Harvesting Systems

Photosynthesis is one of the most essential biological processes in which specialized pigment-protein complexes absorb solar photons, and with a remarkably high efficiency, guide the photo-induced excitation energy toward the reaction center to subsequently trigger its conversion to chemical energy. In this work, we review the principles of optimal energy transfer in various natural and artificial light harvesting systems. We begin by presenting the guiding principles for optimizing the energy transfer efficiency in systems connected to dissipative environments, with particular attention paid to the potential role of quantum coherence in light harvesting systems. We will comment briefly on photo-protective mechanisms in natural systems that ensure optimal functionality under varying ambient conditions. For completeness, we will also present an overview of the charge separation and electron transfer pathways in reaction centers. Finally, recent theoretical and experimental progress on excitation energy transfer, charge separation, and charge transport in artificial light harvesting systems is delineated, with organic solar cells taken as prime examples.

Coherent quantum transport in disordered systems: A unified polaron treatment of hopping and band-like transport

Quantum transport in disordered systems is studied using a polaron-based master equation. The polaron approach is capable of bridging the results from the coherent band-like transport regime governed by the Redfield equation to incoherent hopping transport in the classical regime. A non-monotonic dependence of the diffusion coefficient is observed both as a function of temperature and system-phonon coupling strength. In the band-like transport regime, the diffusion coefficient is shown to be linearly proportional to the system-phonon coupling strength and vanishes at zero coupling due to Anderson localization. In the opposite classical hopping regime, we correctly recover the dynamics described by the Fermi's Golden Rule and establish that the scaling of the diffusion coefficient depends on the phonon bath relaxation time. In both the hopping and band-like transport regimes, it is demonstrated that at low temperature, the zero-point fluctuations of the bath lead to non-zero transport rates and hence a finite diffusion constant. Application to rubrene and other organic semiconductor materials shows a good agreement with experimental mobility data.

Quantum mechanics of excitation transport in photosynthetic complexes: a key issues review

For a long time microscopic physical descriptions of biological processes have been based on quantum mechanical concepts and tools, and routinely employed by chemical physicists and quantum chemists. However, the last ten years have witnessed new developments on these studies from a different perspective, rooted in the framework of quantum information theory. The process that more, than others, has been subject of intense research is the transfer of excitation energy in photosynthetic light-harvesting complexes, a consequence of the unexpected experimental discovery of oscillating signals in such highly noisy systems. The fundamental interdisciplinary nature of this research makes it extremely fascinating, but can also constitute an obstacle to its advance. Here in this review our objective is to provide an essential summary of the progress made in the theoretical description of excitation energy dynamics in photosynthetic systems from a quantum mechanical perspective, with the goal of unifying the language employed by the different communities. This is initially realized through a stepwise presentation of the fundamental building blocks used to model excitation transfer, including protein dynamics and the theory of open quantum system. Afterwards, we shall review how these models have evolved as a consequence of experimental discoveries; this will lead us to present the numerical techniques that have been introduced to quantitatively describe photo-absorbed energy dynamics. Finally, we shall discuss which mechanisms have been proposed to explain the unusual coherent nature of excitation transport and what insights have been gathered so far on the potential functional role of such quantum features.

Reduced quantum dynamics with arbitrary bath spectral densities: hierarchical equations of motion based on several different bath decomposition schemes

We investigated applications of the hierarchical equation of motion (HEOM) method to perform high order perturbation calculations of reduced quantum dynamics for a harmonic bath with arbitrary spectral densities. Three different schemes are used to decompose the bath spectral density into analytical forms that are suitable to the HEOM treatment: (1) The multiple Lorentzian mode model that can be obtained by numerically fitting the model spectral density. (2) The combined Debye and oscillatory Debye modes model that can be constructed by fitting the corresponding classical bath correlation function. (3) A new method that uses undamped harmonic oscillator modes explicitly in the HEOM formalism. Methods to extract system-bath correlations were investigated for the above bath decomposition schemes. We also show that HEOM in the undamped harmonic oscillator modes can give detailed information on the partial Wigner transform of the total density operator. Theoretical analysis and numerical simulations of the spin-Boson dynamics and the absorption line shape of molecular dimers show that the HEOM formalism for high order perturbations can serve as an important tool in studying the quantum dissipative dynamics in the intermediate coupling regime.

Quantifying non-Markovianity for a chromophore-qubit pair in a super-Ohmic bath

An approach based on a non-Markovian time-convolutionless polaron master equation is used to probe the quantum dynamics of a chromophore-qubit in a super-Ohmic bath. Utilizing a measure of non-Markovianity based on dynamical fixed points, we study the effects of the environmental temperature and the coupling strength on the non-Markovian behavior of the chromophore in a super-Ohmic bath. It is found that an increase in the temperature results in a reduction in the backflow information from the environment to the chromophore, and therefore, a suppression of non-Markovianity. In the weak coupling regime, increasing the coupling strength will enhance the non-Markovianity, while the effect is reversed in the strong coupling regime.