Calculation of absorption spectra for light-harvesting systems using non-Markovian approaches as well as modified Redfield theory

Unraveling quantum coherences mediating primary charge transfer processes in photosystem II reaction center

Photosystem II (PSII) reaction center (RC) is a unique complex that is capable of efficiently separating electronic charges across the membrane. The primary energy- and charge-transfer (CT) processes occur on comparable ultrafast timescales, which makes it extremely challenging to understand the fundamental mechanism responsible for the near-unity quantum efficiency of the transfer. Here, we elucidate the role of quantum coherences in the ultrafast energy and CT in the PSII RC by performing two-dimensional (2D) electronic spectroscopy at the cryogenic temperature of 20 kelvin, which captures the distinct underlying quantum coherences. Specifically, we uncover the electronic and vibrational coherences along with their lifetimes during the primary ultrafast processes of energy and CT. We construct an excitonic model that provides evidence for coherent energy and CT at low temperature in the 2D electronic spectra. The principles could provide valuable guidelines for creating artificial photosystems with exploitation of system-bath coupling and control of coherences to optimize the photon conversion efficiency to specific functions.

The location of the low-energy states in Lhca1 favors excitation energy transfer to the core in the plant PSI-LHCI supercomplex

Lhca1 is one of the four pigment-protein complexes composing the outer antenna of plant Photosystem I-light-havesting I supercomplex (PSI-LHCI). It forms a functional dimer with Lhca4 but, differently from this complex, it does not contain 'red-forms,' i.e., pigments absorbing above 700 nm. Interestingly, the recent PSI-LHCI structures suggest that Lhca1 is the main point of delivering the energy harvested by the antenna to the core. To identify the excitation energy pathways in Lhca1, we developed a structure-based exciton model based on the simultaneous fit of the low-temperature absorption, linear dichroism, and fluorescence spectra of wild-type Lhca1 and two mutants, lacking chlorophylls contributing to the long-wavelength region of the absorption. The model enables us to define the locations of the lowest energy pigments in Lhca1 and estimate pathways and timescales of energy transfer within the complex and to the PSI core. We found that Lhca1 has a particular energy landscape with an unusual (compared to Lhca4, LHCII, and CP29) configuration of the low-energy states. Remarkably, these states are located near the core, facilitating direct energy transfer to it. Moreover, the low-energy states of Lhca1 are also coupled to the red-most state (red forms) of the neighboring Lhca4 antenna, providing a pathway for effective excitation energy transfer from Lhca4 to the core.

Recent progress in atomistic modeling of light-harvesting complexes: a mini review

In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.

Electronic Absorption Spectra from Off-Diagonal Quantum Master Equations

Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrodinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrodinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrodinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained in this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion (HEOM) method.

Accuracy of approximate methods for the calculation of absorption-type linear spectra with a complex system-bath coupling

The accuracy of approximate methods for calculating linear optical spectra depends on many variables. In this study, we fix most of these parameters to typical values found in photosynthetic light-harvesting complexes of plants and determine the accuracy of approximate spectra with respect to exact calculation as a function of the energy gap and interpigment coupling in a pigment dimer. We use a spectral density with the first eight intramolecular modes of chlorophyll a and include inhomogeneous disorder for the calculation of spectra. We compare the accuracy of absorption, linear dichroism, and circular dichroism spectra calculated using the Full Cumulant Expansion (FCE), coherent time-dependent Redfield (ctR), and time-independent Redfield and modified Redfield methods. As a reference we use spectra calculated with the Exact Stochastic Path Integral Evaluation method. We find the FCE method to be the most accurate for the calculation of all spectra. The ctR method performs well for the qualitative calculation of absorption and linear dichroism spectra when pigments are moderately coupled (∼15 cm-1), but ctR spectra may differ significantly from exact spectra when strong interpigment coupling (>100 cm-1) is present. The dependence of the quality of Redfield and modified Redfield spectra on molecular parameters is similar, and these methods almost always perform worse than ctR, especially when the interpigment coupling is strong or the excitonic energy gap is small (for a given coupling). The accuracy of approximate spectra is not affected by resonance with intramolecular modes for typical system-bath coupling and disorder values found in plant light-harvesting complexes.

Spectral densities and absorption spectra of the core antenna complex CP43 from photosystem II

Besides absorbing light, the core antenna complex CP43 of photosystem II is of great importance in transferring excitation energy from the antenna complexes to the reaction center. Excitation energies, spectral densities, and linear absorption spectra of the complex have been evaluated by a multiscale approach. In this scheme, quantum mechanics/molecular mechanics molecular dynamics simulations are performed employing the parameterized density functional tight binding (DFTB) while the time-dependent long-range-corrected DFTB scheme is applied for the excited state calculations. The obtained average spectral density of the CP43 complex shows a very good agreement with experimental results. Moreover, the excitonic Hamiltonian of the system along with the computed site-dependent spectral densities was used to determine the linear absorption. While a Redfield-like approximation has severe shortcomings in dealing with the CP43 complex due to quasi-degenerate states, the non-Markovian full second-order cumulant expansion formalism is able to overcome the drawbacks. Linear absorption spectra were obtained, which show a good agreement with the experimental counterparts at different temperatures. This study once more emphasizes that by combining diverse techniques from the areas of molecular dynamics simulations, quantum chemistry, and open quantum systems, it is possible to obtain first-principle results for photosynthetic complexes, which are in accord with experimental findings.

Quantum-Classical Approach for Calculations of Absorption and Fluorescence: Principles and Applications

Absorption and fluorescence spectroscopy techniques provide a wealth of information on molecular systems. The simulations of such experiments remain challenging, however, despite the efforts put into developing the underlying theory. An attractive method of simulating the behavior of molecular systems is provided by the quantum-classical theory─it enables one to keep track of the state of the bath explicitly, which is needed for accurate calculations of fluorescence spectra. Unfortunately, until now there have been relatively few works that apply quantum-classical methods for modeling spectroscopic data. In this work, we seek to provide a framework for the calculations of absorption and fluorescence lineshapes of molecular systems using the methods based on the quantum-classical Liouville equation, namely, the forward-backward trajectory solution (FBTS) and the non-Hamiltonian variant of the Poisson bracket mapping equation (PBME-nH). We perform calculations on a molecular dimer and the photosynthetic Fenna-Matthews-Olson complex. We find that in the case of absorption, the FBTS outperforms PBME-nH, consistently yielding highly accurate results. We next demonstrate that for fluorescence calculations, the method of choice is a hybrid approach, which we call PBME-nH-Jeff, that utilizes the effective coupling theory [Gelzinis, A.; J. Chem. Phys. 2020, 152, 051103] to estimate the excited state equilibrium density operator. Thus, we find that FBTS and PBME-nH-Jeff are excellent candidates for simulating, respectively, absorption and fluorescence spectra of real molecular systems.

Time-dependent atomistic simulations of the CP29 light-harvesting complex

Light harvesting as the first step in photosynthesis is of prime importance for life on earth. For a theoretical description of photochemical processes during light harvesting, spectral densities are key quantities. They serve as input functions for modeling the excitation energy transfer dynamics and spectroscopic properties. Herein, a recently developed procedure is applied to determine the spectral densities of the pigments in the minor antenna complex CP29 of photosystem II, which has recently gained attention because of its active role in non-photochemical quenching processes in higher plants. To this end, the density functional-based tight binding (DFTB) method has been employed to enable simulation of the ground state dynamics in a quantum-mechanics/molecular mechanics (QM/MM) scheme for each chlorophyll pigment. Subsequently, the time-dependent extension of the long-range corrected DFTB approach has been used to obtain the excitation energy fluctuations along the ground-state trajectories also in a QM/MM setting. From these results, the spectral densities have been determined and compared for different force fields and to spectral densities from other light-harvesting complexes. In addition, time-dependent and time-independent excitonic Hamiltonians of the system have been constructed and applied to the determination of absorption spectra as well as exciton dynamics.

Correlated vibration-solvent effects on the non-Condon exciton spectroscopy

Excitation energy transfer is crucially involved in a variety of systems. During the process, the non-Condon vibronic coupling and the surrounding solvent interaction may synergetically play important roles. In this work, we study the correlated vibration-solvent influences on the non-Condon exciton spectroscopy. Statistical analysis is elaborated for the overall vibration-plus-solvent environmental effects. Analytic solutions are derived for the linear absorption of monomer systems. General simulations are accurately carried out via the dissipaton-equation-of-motion approach. The resulted spectra in either the linear absorption or strong field regime clearly demonstrate the coherence enhancement due to the synergetic vibration-solvent correlation.

Effect of logarithmic perturbations in ohmic like spectral densities in dynamics of electronic excitation using variational polaron transformation approach

Electronic excitation energy transfer is a ubiquitous process that has generated prime research interest since its discovery. Recently developed variational polaron transformation-based second-order master equation is capable of interpolating between Förster and Redfield limits with exceptional accuracy. Forms of spectral density functions studied so far through the variational approach provide theoretical support for various experiments. Recently introduced ohmic like spectral density function that can account for logarithmic perturbations provides generality and exposition to a unique and practical set of environments. In this paper, we exploit the energy transfer dynamics of a two-level system attached to an ohmic like spectral density function with logarithmic perturbations using a variational polaron transformed master equation. Our results demonstrate that even for a relatively large bath coupling strength, quantum coherence effects can be increased by introducing logarithmic perturbations of the order of one and two in super-ohmic environments. Moreover, for particular values of the ohmicity parameter, the effect of logarithmic perturbations is observed to be insignificant for the overall dynamics. In regard to ohmic environments, as logarithmic perturbations increase, damping characteristics of the coherent transient dynamics also increase in general. It is also shown that, having logarithmic perturbations of the order of one in an ohmic environment can result in a less efficient energy transfer for relatively larger system bath coupling strengths.

Absorption and Circular Dichroism Spectra of Molecular Aggregates With the Full Cumulant Expansion

The exciton Hamiltonian of multichromophoric aggregates can be probed by spectroscopic techniques such as linear absorption and circular dichroism. To compare calculated Hamiltonians to experiments, a lineshape theory is needed, which takes into account the coupling of the excitons with inter- and intramolecular vibrations. This coupling is normally introduced in a perturbative way through the cumulant expansion formalism and further approximated by assuming a Markovian exciton dynamics, for example with the modified Redfield theory. Here, we present the implementation of the full cumulant expansion (FCE) formalism (J. Chem. Phys.142, 2015, 09410625747060) to efficiently compute absorption and circular dichroism spectra of molecular aggregates beyond the Markov approximation, without restrictions on the form of exciton–phonon coupling. By employing the LH2 system of purple bacteria as a challenging test case, we compare the FCE lineshapes with the Markovian lineshapes obtained with the modified Redfield theory, showing that the latter presents a less satisfying agreement with experiments. The FCE approach instead accurately describes the lineshapes, especially in the vibronic sideband of the B800 peak. We envision that the FCE approach will become a valuable tool for accurately comparing model exciton Hamiltonians with optical spectroscopy experiments.

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.

Effect of Time-Delayed Feedback on the Interaction of a Dimer System with its Environment

In this work, we report modeling of non-Markovian open quantum systems, consisting of an excitonic dimer that displays memory effect due to time delayed interaction with its environment. We, indeed investigate the effect of these time delays on quantum coherence and excitation dynamical behavior in the time domain generally considered for photosynthetic experiments (few hundred femtoseconds). In particular, we show that the coherence is maintained for periods proportional to time delays. Additionally, if delay is taken into account, coupling to the environment can be tuned to lower values, unlike in previous studies. This kind of intriguing effect can, therefore, when generalized to complete systems, permit more control on the experimental parameters, which may lead to more accurate description of the photosynthetic energy transfer functioning and subsequent applications in artificial photovoltaic research.

General Galilei Covariant Gaussian Maps

We characterize general non-Markovian Gaussian maps which are covariant under Galilean transformations. In particular, we consider translational and Galilean covariant maps and show that they reduce to the known Holevo result in the Markovian limit. We apply the results to discuss measures of macroscopicity based on classicalization maps, specifically addressing dissipation, Galilean covariance and non-Markovianity. We further suggest a possible generalization of the macroscopicity measure defined by Nimmrichter and Hornberger [Phys. Rev. Lett. 110, 16 (2013)PRLTAO0031-9007].

The role of interparticle interaction and environmental coupling in a two-particle open quantum system

The effects of bath coupling on an interacting two-particle quantum system are studied using tools from information theory. Shannon entropies of the one (reduced) and two-particle distribution functions in position, momentum and separable phase-space are examined. Results show that the presence of the bath leads to a delocalization of the distribution functions in position space, and a localization in momentum space. This can be interpreted as a loss of information in position space and a gain of information in momentum space. The entropy sum of the system, in the presence of a bath, is shown to be dependent on the strength of the interparticle potential and also on the strength of the coupling to the bath. The statistical correlation between the particles, and its dependence on the bath and interparticle potential, is examined using mutual information. A stronger repulsive potential between particles, in the presence of the bath, yields a smaller correlation between the particles positions, and a larger one between their momenta.

Electron transfer within a reaction path model calibrated by constrained DFT calculations: application to mixed-valence organic compounds

The quantum dynamics of electron transfer in mixed-valence organic compounds is investigated using a reaction path model calibrated by constrained density functional theory (cDFT). Constrained DFT is used to define diabatic states relevant for describing the electron transfer, to obtain equilibrium structures for each of these states and to estimate the electronic coupling between them. The harmonic analysis at the diabatic minima yields normal modes forming the dissipative bath coupled to the electronic states. In order to decrease the system-bath coupling, an effective one dimensional vibronic Hamiltonian is constructed by partitioning the modes into a linear reaction path which connects both equilibrium positions and a set of secondary vibrational modes, coupled to this reaction coordinate. Using this vibronic model Hamiltonian, dissipative quantum dynamics is carried out using Redfield theory, based on a spectral density which is determined from the cDFT results. In a first benchmark case, the model is applied to a series of mixed-valence organic compounds formed by two 1,4-dimethoxy-3-methylphenylene fragments linked by an increasing number of phenylene bridges. This allows us to examine the coherent electron transfer in extreme situations leading to a ground adiabatic state with or without a barrier and therefore to the trapping of the charge or to an easy delocalization.

Efficient tool to calculate two-dimensional optical spectra for photoactive molecular complexes

We combine the coherent modified Redfield theory (CMRT) with the equation of motion-phase matching approach (PMA) to calculate two-dimensional photon-echo spectra for photoactive molecular complexes with an intermediate strength of the coupling to their environment. Both techniques are highly efficient, yet they involve approximations at different levels. By explicitly comparing with the numerically exact quasiadiabatic path integral approach, we show for the Fenna-Matthews-Olson complex that the CMRT describes the decay rates in the population dynamics well, but final stationary populations and the oscillation frequencies differ slightly. In addition, we use the combined CMRT+PMA to calculate two-dimensional photon-echo spectra for a simple dimer model. We find excellent agreement with the exact path integral calculations at short waiting times where the dynamics is still coherent. For long waiting times, differences occur due to different final stationary states, specifically for strong system-bath coupling. For weak to intermediate system-bath couplings, which is most important for natural photosynthetic complexes, the combined CMRT+PMA gives reasonable results with acceptable computational efforts.

Path induced coherent energy transfer in light-harvesting complexes in purple bacteria

Features of path dependent energy transfer in a dual-ring light-harvesting (LH2) complexes (B850) system have been examined in detail systematically. The Frenkel-Dirac time dependent variational method with the Davydov D1 Ansatz is employed with detailed evolution of polaron dynamics in real space readily obtained. It is found that the phase of the transmission amplitude through the LH2 complexes plays an important role in constructing the coherent excitonic energy transfer. It is also found that the symmetry breaking caused by the dimerization of bacteriochlorophylls and coherence or correlation between two rings will be conducive in enhancing the exciton transfer efficiency.

Förster resonance energy transfer, absorption and emission spectra in multichromophoric systems. I. Full cumulant expansions and system-bath entanglement

We study the Förster resonant energy transfer rate, absorption and emission spectra in multichromophoric systems. The multichromophoric Förster theory (MCFT) is determined from an overlap integral of generalized matrices related to the donor's emission and acceptor's absorption spectra, which are obtained via a full 2nd-order cumulant expansion technique developed in this work. We calculate the spectra and MCFT rate for both localized and delocalized systems, and calibrate the analytical 2nd-order cumulant expansion with the exact stochastic path integral method. We present three essential findings: (i) The role of the initial entanglement between the donor and its bath is found to be crucial in both the emission spectrum and the MCFT rate. (ii) The absorption spectra obtained by the cumulant expansion method are nearly identical to the exact spectra for both localized and delocalized systems, even when the system-bath coupling is far from the perturbative regime. (iii) For the emission spectra, the cumulant expansion can give reliable results for localized systems, but fail to provide reliable spectra of the high-lying excited states of a delocalized system, when the system-bath coupling is large and the thermal energy is small. This paper also provides a simple golden-rule derivation of the MCFT, reviews existing methods, and motivates further developments in the subsequent papers.

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.

Modeling of fluorescence line-narrowed spectra in weakly coupled dimers in the presence of excitation energy transfer

This work describes simple analytical formulas to describe the fluorescence line-narrowed (FLN) spectra of weakly coupled chromophores in the presence of excitation energy transfer (EET). Modeling studies for dimer systems (assuming low fluence and weak coupling) show that the FLN spectra (including absorption and emission spectra) calculated for various dimers using our model are in good agreement with spectra calculated by: (i) the simple convolution method and (ii) the more rigorous treatment using the Redfield approach [T. Renger and R. A. Marcus, J. Chem. Phys. 116, 9997 (2002)]. The calculated FLN spectra in the presence of EET of all three approaches are very similar. We argue that our approach provides a simplified and computationally more efficient description of FLN spectra in the presence of EET. This method also has been applied to FLN spectra obtained for the CP47 antenna complex of Photosystem II reported by Neupane et al. [J. Am. Chem. Soc. 132, 4214 (2010)], which indicated the presence of uncorrelated EET between pigments contributing to the two lowest energy (overlapping) exciton states, each mostly localized on a single chromophore. Calculated and experimental FLN spectra for CP47 complex show very good qualitative agreement.

On the Controversial Nature of the 825 nm Exciton Band in the FMO Protein Complex

The nature of the low-energy 825 nm band of the Fenna-Matthews-Olson (FMO) protein complex from Chlorobaculum tepidum at 5 K is discussed. It is shown, using hole-burning (HB) spectroscopy and excitonic calculations, that the 825 nm absorption band of the FMO trimer cannot be explained by a single electronic transition or overlap of electronic transitions of noninteracting pigments. To explain the shape of emission and nonresonant HB spectra, downward uncorrelated excitation energy transfer (EET) between trimer subunits should be taken into account. Modeling studies reveal the presence of three sub-bands within the 825 nm band, in agreement with nonresonant HB and emission spectra. We argue that after light induced coherences vanish, uncorrelated EET between the lowest exciton levels of each monomer takes place. HB induced spectral shifts provide a new insight on the energy landscape of the FMO protein.

Photosynthetic light harvesting: excitons and coherence

Photosynthesis begins with light harvesting, where specialized pigment-protein complexes transform sunlight into electronic excitations delivered to reaction centres to initiate charge separation. There is evidence that quantum coherence between electronic excited states plays a role in energy transfer. In this review, we discuss how quantum coherence manifests in photosynthetic light harvesting and its implications. We begin by examining the concept of an exciton, an excited electronic state delocalized over several spatially separated molecules, which is the most widely available signature of quantum coherence in light harvesting. We then discuss recent results concerning the possibility that quantum coherence between electronically excited states of donors and acceptors may give rise to a quantum coherent evolution of excitations, modifying the traditional incoherent picture of energy transfer. Key to this (partially) coherent energy transfer appears to be the structure of the environment, in particular the participation of non-equilibrium vibrational modes. We discuss the open questions and controversies regarding quantum coherent energy transfer and how these can be addressed using new experimental techniques.

Normal mode analysis of the spectral density of the Fenna-Matthews-Olson light-harvesting protein: how the protein dissipates the excess energy of excitons

We report a method for the structure-based calculation of the spectral density of the pigment-protein coupling in light-harvesting complexes that combines normal-mode analysis with the charge density coupling (CDC) and transition charge from electrostatic potential (TrEsp) methods for the computation of site energies and excitonic couplings, respectively. The method is applied to the Fenna-Matthews-Olson (FMO) protein in order to investigate the influence of the different parts of the spectral density as well as correlations among these contributions on the energy transfer dynamics and on the temperature-dependent decay of coherences. The fluctuations and correlations in excitonic couplings as well as the correlations between coupling and site energy fluctuations are found to be 1 order of magnitude smaller in amplitude than the site energy fluctuations. Despite considerable amplitudes of that part of the spectral density which contains correlations in site energy fluctuations, the effect of these correlations on the exciton population dynamics and dephasing of coherences is negligible. The inhomogeneous charge distribution of the protein, which causes variations in local pigment-protein coupling constants of the normal modes, is responsible for this effect. It is seen thereby that the same building principle that is used by nature to create an excitation energy funnel in the FMO protein also allows for efficient dissipation of the excitons' excess energy.

Juxtaposing density matrix and classical path-based wave packet dynamics

In many physical, chemical, and biological systems energy and charge transfer processes are of utmost importance. To determine the influence of the environment on these transport processes, equilibrium molecular dynamics simulations become more and more popular. From these simulations, one usually determines the thermal fluctuations of certain energy gaps, which are then either used to perform ensemble-averaged wave packet simulations, also called Ehrenfest dynamics, or to employ a density matrix approach via spectral densities. These two approaches are analyzed through energy gap fluctuations that are generated to correspond to a predetermined spectral density. Subsequently, density matrix and wave packet simulations are compared through population dynamics and absorption spectra for different parameter regimes. Furthermore, a previously proposed approach to enforce the correct long-time behavior in the wave packet simulations is probed and an improvement is proposed.

Open Quantum Dynamics Calculations with the Hierarchy Equations of Motion on Parallel Computers

Calculating the evolution of an open quantum system, i.e., a system in contact with a thermal environment, has presented a theoretical and computational challenge for many years. With the advent of supercomputers containing large amounts of memory and many processors, the computational challenge posed by the previously intractable theoretical models can now be addressed. The hierarchy equations of motion present one such model and offer a powerful method that remained under-utilized so far due to its considerable computational expense. By exploiting concurrent processing on parallel computers the hierarchy equations of motion can be applied to biological-scale systems. Herein we introduce the quantum dynamics software PHI, that solves the hierarchical equations of motion. We describe the integrator employed by PHI and demonstrate PHI's scaling and efficiency running on large parallel computers by applying the software to the calculation of inter-complex excitation transfer between the light harvesting complexes 1 and 2 of purple photosynthetic bacteria, a 50 pigment system.