Environment-Assisted Quantum Coherence in Photosynthetic

Oct 30, 2017 - Historically, most of the studies of EET considered either the Markovian limit or quantum master equation approach where perturbative t...
11 downloads 15 Views 1MB Size
Subscriber access provided by READING UNIV

Letter

Environment Assisted Quantum Coherence in Photosynthetic Complex Rajesh Dutta, and Biman Bagchi J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b02480 • Publication Date (Web): 30 Oct 2017 Downloaded from http://pubs.acs.org on November 1, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Environment Assisted Quantum Coherence in Photosynthetic Complex Rajesh Dutta1 and Biman Bagchi1,* 1

SSCU, Indian Institute of Science, Bangalore 560012, India. *Email: [email protected]

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 22

Abstract Recent experiments [Nature, 2007, 446, 782-786] revealed the existence of surprisingly long lived quantum coherence in noisy biological environment of photosynthetic Fenna-MatthewsOlson (FMO) complex. Such coherence can clearly play 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 [Science, 2009, 323, 369-373]. Even though these systems are strongly correlated, most theoretical studies invoke Markovian approximation where the temporal correlation of bath fluctuations is neglected. We use an elegant non-perturbative 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.

TOC Graphic

KEYWORDS: FMO complex, non-Markovian bath, Excitation energy transfer, Fluctuation assisted coherence, Coherent and Incoherent transfer

2 ACS Paragon Plus Environment

Page 3 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Quantum transport of charge and energy can be quite distinct from their classical analogues. In the classical limit, we can imagine a relation between diffusion and fluctuating forces in terms of Einstein’s picture. The situation is quite different in the quantum case, as exemplified by profound effects of disorder in Anderson localization.1 Existence of coherence between different sites and/or energy levels plays an important role in quantum transport. Presence of disorder can destroy this coherence. A large number of studies have been devoted to include the effects of static disorder on quantum transport or diffusion (as in Anderson localization problem or in phonon localization).2 In contrast, relatively fewer studies have considered effects of dynamical disorder in transport properties.3 The situation has undergone a see-saw change in recent years, triggered by recent results in excitation energy transfer (EET) in photosynthetic complex4-9 and conjugate polymers.10-13 The main reason behind the attention is the discovery of coherent transport in the noisy biological environment. The important application of EET is that the efficiency of energy transfer from chromophores to reaction center in photosynthetic system is near about unity.14 If it can be mimicked in artificial solar cells, one can design more efficient solar cells. Historically, most of the studies of EET considered either the Markovian limit or quantum master equation approach where perturbative truncation is necessary. However, these approaches could be problematic in the cases of both photosynthetic complex and conjugate polymers where the EET timescales and the environment bath correlation time scales are comparable.5-7 In a series of pioneering studies using 2D Fourier transform electronic spectroscopy, Fleming and co-workers4 observed long lived coherent EET in photosynthetic complexes. Collini and Scholes10,11 have observed coherent intra-chain EET in conjugated polymer at room temperature. Later Ishizaki and Fleming studied theoretically explore the role of environment in energy

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 22

transport through the consideration of reorganization for FMO15 complex. In this approach infinite set of equations are truncated after finite numbers. Another well-known approach is based on polaron transfer technique developed by Jang and co-workers16 to investigate the energy transfer dynamics. However, the second-order perturbative truncation with respect to the renormalized electron-phonon coupling is very crude approximation in case photosynthetic systems. Few years ago, Aspuru-Guzik and co-workers17,18 have described environment assisted quantum transport in real photosynthetic complex and model systems. Silbey and co-workers19-23 have computed the efficiency of exciton transfer in case of dimer model system and population relaxation for FMO complex. However, use of Bourret’s integrals is limited to fast fluctuations or short bath correlation time. Despite all the studies, however, there is yet no study on the effects of temporal correlation of bath in a non-Markovian environment on transport in an extended system. In this study we explore coherent and incoherent transport of excitation energy as a function of bath correlation time and the fluctuation strength. In our earlier studies we observed quantum diffusion, line shape and occupation probability in presence of dynamical disorder for several model systems.24,25 Recently Moix, Khasin and Cao26 have concluded that although the static disorder leads to lack of diffusion and localization, a small amount of dynamic disorder is enough to continue the exciton transport. In the present study we consider real photosynthetic FMO complex such that we can investigate the localization and delocalization excitation energy in presence of both heterogeneity of site energies and dynamic disorders. Inclusion of memory effects in the combined system-bath dynamics is highly non-trivial and is rarely attempted. Usual treatments invoke a Markovian approximation where bath dynamics is delta correlated in time. There is however, an alternative powerful approach to

4 ACS Paragon Plus Environment

Page 5 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

include memory effects introduced by Kubo which is known as quantum stochastic Liouville equation (QSLE).27 In this approach, one needs to specify the form of bath evolution operator as described below. In the present study, we consider the total Hamiltonian in the interaction representation of the fluctuating bath Hamiltonian, so that the interaction V is time dependent and can be modeled by a stochastic function with known statistical properties. In this interaction representation, the Hamiltonian can be represented as follows,

H (t) = HS + V (t)

(1)

We use the Hamiltonian based on the theories of Haken-Strobl-Reineker28,29 and Silbey30 where the system (exciton) Hamiltonian (Hs) is defined as

HS =

∑E

k

k

k

k +

∑J

kl

k

l

.

(2)

k ,l k ≠l

Here, Ek is the energy of an electronic exciton localized at site k in absence of any coupling and

Jkl is the time independent off-diagonal interaction or inter-site coupling between excitations at site k and l. We assume J kl = J kl* = Jlk = J . Thus in the interaction representation, the time dependent coupling potential can be written as,

V (t ) = Vd (t )∑ k k + Vod (t )∑ k l k

k ,l k ≠l

.

(3)

Vint is the interaction Hamiltonian between the system (exciton) and the environment. Here, Vd ( t ) and Vod ( t ) denote diagonal (local) and off-diagonal (non-local) parts of the fluctuating Hamiltonian Vint . We assume V(t) to be described either by a Poisson stochastic process or a 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 22

Gaussian Markov process.31 We use Kubo’s QSLE (see Supporting Information (SI)) and the equation is provided as follows ∂σ i = − [ H (t ), σ ] + Γ vσ . ∂t h

(4)

where σ is the reduced density matrix (averaged over system variables) defined using joint probability distribution, Γ v is the stochastic diffusion operator (see SI) and V is the same random energy variable. The QSLE was used to explain exciton transport in one dimensional lattice considering correlated and uncorrelated bath for continuum model.32 The non-perturbative approach can be used to study both strong coupling non-Markovian limit and weak coupling Markovian limit. In weak coupling and Markovian limit EET ensues with the equilibrium phonon states and classical rate expression is appropriate to explain the EET dynamics. The behavior can be used to explain FRET.33 We consider two systems cyclic hetero trimer and heptamer; and use system parameters from FMO Hamiltonian (SI). All the calculations are performed considering uncorrelated diagonal fluctuations because of negligibly small off-diagonal fluctuations in photosynthetic systems. Here uncorrelated diagonal fluctuations refer to the situation where all the fluctuating diagonal elements in the Hamiltonian are independent of each other. The coupled equations of motion (EOM) for cyclic trimer is given as follows x

∂σ jmn ∂t

x

i  = −  ∑ Ek k k  σ jmn h  k =1  3

 3  1 x i i −  ∑ J kl k l  σ jmnp − Vd ∑ (δ j +1, j ′ + δ j −1, j′ ) ( 1 1 ) σ j ′mn   h  k =1,l =1 h j′= 0   k ≠l 

1 i − Vd ∑ (δ m +1,m′ + δ m −1,m′ ) ( 2 2 h m′ = 0

)

x

i h

1

σ jm′n − Vd ∑ (δ n +1,n′ + δ n −1,n′ ) ( 3 3 ) σ jmn′ x

n′ = 0

− jbdσ jmnp − mbdσ jmnp − nbdσ jmnp .

(5) 6

ACS Paragon Plus Environment

Page 7 of 22

where, Vd is the strength of the fluctuation and bd corresponds to the rate of fluctuation. We assume uncorrelated nearest neighbor diagonal fluctuations have same strength and rate. We have numerically solved the coupled EOM using Runge-Kutta fourth order method. Population of each site can be represented as follows, Pn ( t ) = n σ

N

∏ ai

n

(6)

i =1

where, n is the site number, N is the total number of site and a1 , a2 , a3 ,....aN all the elements are zero that denotes the equilibrium states of Poisson bath. 1

0.6

0.8

0.4 0.2 0 0

0.6 0.4 0.2

0.5

1 1.5 Time (ps)

2

0 0

2.5

0.5

1 1.5 Time (ps)

2

2.5

0.3

0.5

(c)

||

0.4

||

(d) 0.25

||

Off-diagonal

||

0.3

||

0.2 0.1 0 0

Site1 Site2 Site3 Site1 Site2 Site3

(b) Population

0.8 Population

1

Site1 Site2 Site3 Site1 Site2 Site3

(a)

Off-diagonal

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

|| ||

0.2

||

0.15 0.1 0.05

0.5

1 1.5 Time (ps)

2

2.5

0 0

0.5

1 1.5 Time (ps)

2

2.5 −1

Figure 1.(a)-(b) depict population dynamics in trimer system. (a) Solid line :Vd= 50 cm-1 and bd = −1

−1

10 fs., dashed line :Vd= 50 cm-1 and bd = 250 fs. (b) Solid line :Vd= 350 cm-1 and bd = 10 fs., dashed 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

−1

line : Vd= 350 cm-1 and bd = 250 fs. Figure (c)-(d) indicate absolute value of coherence(c) Solid line −1

−1

:Vd= 50 cm-1 and bd = 10 fs., dashed line Vd= 50 cm-1 and bd = 250 fs. (d) Solid line :Vd= 350 cm-1 −1

−1

and bd = 10 fs, dashed line : Vd= 350 cm-1 and bd = 250 fs. 1

1 Site1 Site2 Site3 Site4 Site5 Site6 Site7

Population

0.8 0.6

0.8

0.4

0.6 0.4 0.2

0.2 0 0

0.5

1 1.5 Time (ps)

0.5

2

0 0

2.5

Off-diagonal

|| ||

0.2

2

||

0.1

2.5

|| ||

0.4

||

0.3

1 1.5 Time (ps)

(d)

||

0.4

0.5

0.5

||

(c) Off-diagonal

Site1 Site2 Site3 Site4 Site5 Site6 Site7

(b) Population

(a)

|| ||

0.3

||

0.2

||

0.1

0 0

0.5

1 1.5 Time (ps)

2

0 0

2.5

0.5

1 1.5 Time (ps)

2

2.5

Figure 2. .(a)-(b) depict population dynamics of FMO complex (cyclic seven sites system). (a) Vd= 50 −1 −1 cm-1 and bd = 10 fs. (b) Vd= 50 cm-1 and bd = 250 fs. Figure (c)-(d) signify the absolute value of −1

coherences of FMO complex (cyclic seven sites system), with (c): Vd= 50 cm-1 and bd = 10 fs. −1

(d):Vd= 50 cm-1 and bd = 250 fs..

1

1 Site1 Site2 Site3 Site4 Site5 Site6 Site7

0.8 0.6 0.4

0.8

0.2 0 0

Site1 Site2 Site3 Site4 Site5 Site6 Site7

(b) Population

(a) Population

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 22

0.6 0.4 0.2

0.5

1 1.5 Time (ps)

2

2.5

0 0

ACS Paragon Plus Environment

0.5

1 1.5 Time (ps)

2

2.5 8

Page 9 of 22

0.25

0.3

||

(c)

Off-diagonal

|| ||

0.15

||

0.1

||

0.05 0 0

||

(d)

0.25

||

0.2 Off-diagonal

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

||

0.2

||

0.15

||

|| ||

0.1 0.05

0.2

0.4 0.6 Time (ps)

0.8

1

0 0

0.2

0.4 0.6 Time (ps)

0.8

1

Figure 3.(a)-(b) depict population dynamics of FMO complex (cyclic seven sites system). (a)Vd= 350 −1 −1 cm-1 and bd = 10 fs. (b) Vd= 350 cm-1 and bd = 250 fs. Figure (c)-(d) indicates absolute value of −1

coherences of FMO complex (cyclic seven sites system), with (c): Vd= 350 cm-1 and bd = 10 fs. −1

(d):Vd= 350 cm-1 and bd = 250 fs..

In the absence of bath fluctuations, coherence is trivial because in this case it occurs due to the presence of the constant inter-site coupling J. However, in the presence of bath fluctuations several possibilities arise, one must consider both the role of bath fluctuation strength (Vd) and −1

bath correlation time ( bd ). In the absence of any off-diagonal coupling J, transition from nonMarkovian to Markovian limit can be obtained by increasing bd alone, at fixed Vd which is the real picture in case of vibrational energy relaxation.31 However, the situation is complex in the presence of a significant J, as shown below. According to recent studies mentioned earlier, EET in photosynthetic complex as well as in conjugated polymers occurs in the intermediate regime where EET and bath correlation time scales are almost similar. These studies indicate that EET in photosynthetic or conjugated polymer occurs in non-Markovian limit. In this study strong coupling denotes dominance of time independent off-diagonal coupling J. For the fixed value of J in case of FMO complex, value of −1 Vd= 50 cm-1 and bd = 250 fs. falls in the strong coupling limit. However, in this case weak −1 coupling refers Vd= 350 cm-1 and bd = 10 fs. In this study we consider intermediate limits as

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 22

−1 −1 Vd= 50 cm-1 and bd = 10 fs. as well as Vd= 350 cm-1 and bd = 250 fs. Both strong coupling and

intermediate regime correspond to non-Markovian limit. Markovian or weak coupling limit can −1 be obtained when Vd is large and bd is small. In this limit

Vd2 act as an effective bath relaxation bd

Vd2 rate. If is much larger than the energy transfer time or time the scale associated with J, bd energy transfer can be incoherent. The Figures 1-3 show population and absolute value of coherences for trimer as well as FMO complex. Initially site 2 gains most of the population by taking the advantage of quantum coherence. In case of strong coupling limit, population shows oscillatory dynamics i.e. EET is completely coherent up to 3 ps. However, in intermediate coupling regime duration of oscillation decreases and we observe transition from coherent to incoherent transport of energy at 1 ps. In Fig. 1a and 1b we plot population for small as well as large bath correlation time, for two fixed values of fluctuation strength. Even when the bath correlation time (inverse of bd) is small, increase in the fluctuation strength (Vd) may not be suffice to reduce or remove the effect of the electronic coupling J. Consequently, coherence and oscillations may survive for a long time. One of the main parameters that determines the transport mechanism is

Vd2 and not only bd bd

alone, but also J. A major outcome of the present work is the understanding of the precise role played by the static (time averaged) off-diagonal coupling J in a non-Markovian description. This aspect is often not realized in a Markovian description that does not consider Vd and bd separately.

10 ACS Paragon Plus Environment

Page 11 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Vd2 If the ratio ( ) is much greater than J, and at the same time both bd and Vd are larger bd than J, then the motion can become incoherent. However, if bd alone becomes large at a fixed Vd, the ratio

Vd2 may become smaller than J. In this limit, oscillation could last for a long time. bd

Similarly one can explain the results for FMO complex as well. In the intermediate coupling regime when both fluctuation strength and bath correlation time are large, coherence could be destroyed by high fluctuation strength. In this case, one can induce a crossover at intermediate times from coherence to incoherence by changing the bath correlation time. One can also explain similar behavior where fluctuation strength and bath correlation time scale are small or fluctuation is very fast. The asymmetry between the behavior of both Vd and

bd−1 as well as the

Vd2 ratio helps in coherent transport and consequently one can term this as an bd

environment or fluctuation assisted coherent EET. We also plot the absolute value of coherences in equilibrium bath states in all the limits. All the non-local coherences actually denote the non-local quantum correlations and create new pathways of EET and interference between pathways as well. The non-local terms have less contribution when environment effect is strong. Coherences in excited bath states have negligibly small contributions in the FMO complex that vary from 10-2 to 10-6 with increasing bath correlation time scales. In case of FMO complex all the site energies and coupling elements are different. Consequently there is a possibility of competition between localization due to site energy heterogeneities and delocalization due to dynamic disorders (bath fluctuations). In this study we observe slow non11 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 22

oscillatory EET dynamics which essentially denote the dominance of the effect of heterogeneity in site energies. Thus the only possibility of energy transfer is the incoherent or hopping mechanism in the long time limit. We consider the temperature corrected quantum stochastic Liouville equation derived by Tanimura and Kubo.34,35 Temperature independent quantum stochastic Liouville equation ignores the reaction of the system to the bath. The high temperature classical bath approximation leads to equal population of all sites in long time and the system does not attain Boltzmann distribution. The main difference between temperature corrected and the temperature independent stochastic Liouville equation is that, in addition to the stochastic interaction term it contains dissipation term which is missing in the stochastic model and helps to capture the effect of bath on system dynamics. These two terms are related through the fluctuation–dissipation theorem leads the system to attain Boltzmann distribution at finite temperature and at long time limit. Temperature corrected stochastic Liouville equation for Gaussian modulation case can be written as follows ∂σ ( V,t ) ∂t

 i ∂  ∂  iβ b  ∂  o =  − H (t ) x − b V + − V +  V σ ( V,t ) ∂V  ∂V  2  ∂V    h

(7)

where, as before, H (t ) = H ex + V (t ) and V is the random variable for the Gaussian modulation. First two terms on the right hand side of Eq. (7) is essentially same as the QSLE equation. Here,

β=

1 o , O f = Of + fO . The term associated with β is the reaction of the system to the bath k BT

and is the temperature correction term that helps the system to reach thermal equilibrium. The set of coupled EOM in the hierarchy is now given as follows ∂σ m β bh o  i i i = − H exx σ m − V xσ m +1 −  V x − V  σ m −1 − mbσ m ∂t h h h 2 

(8) 12

ACS Paragon Plus Environment

Page 13 of 22

It is well known that the truncation after second bath state (i.e. considering only 0 and 1 bath states) in QSLE equation provides the exact expression for Poisson bath or the two state jump model. In the present case, where we consider fluctuation as a Gaussian stochastic process, we include up to fourth bath state from the infinite number of bath states. (see SI, Eq.S13 for coupled EOM). 1

0.6

0.8

0.4

0.4

1 1.5 Time (ps)

2

0.8 0.6

0.5

1 1.5 Time (ps)

1 Site1 (Temp. Dep.) Site2 (Temp. Dep.) Site3 (Temp. Dep.) Site1 (Temp. Indep.) Site2 (Temp. Indep.) Site3 (Temp. Indep.)

(c)

0 0

2.5

2

2.5 Site1 (Temp. Dep.) Site2 (Temp. Dep.) Site3 (Temp. Dep.) Site1 (Temp. Indep.) Site2 (Temp. Indep.) Site3 (Temp. Indep.)

(d) 0.8 Population

0.5

1

0.4

0.6 0.4 0.2

0.2 0 0

0.6

0.2

0.2 0 0

Site1 (Temp. Dep.) Site2 (Temp. Dep.) Site3 (Temp. Dep.) Site1 (Temp. Indep.) Site2 (Temp. Indep.) Site3 (Temp. Indep.)

(b) Population

0.8 Population

1

Site1 (Temp. Dep.) Site2 (Temp. Dep.) Site3 (Temp. Dep.) Site1 (Temp. Indep.) Site2 (Temp. Indep.) Site3 (Temp. Indep.)

(a)

Population

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

0.5

1 1.5 Time (ps)

2

2.5

0 0

0.5

1 1.5 Time (ps)

2

2.5 −1

Figure 4.(a)-(d) corresponds to population dynamics of trimer. (a): Vd= 50 cm-1 and bd = 10 fs. (b): −1

−1

−1

Vd= 50 cm-1 and bd = 250 fs. (c): Vd= 350 cm-1 and bd = 10 fs. (d): Vd= 350 cm-1 and bd = 250 fs. Solid line corresponds to temperature independent study whereas dotted line corresponds to temperature independent study. Effect of temperature is investigated at T = 300 K.

Figure 4 shows the effects of temperature in excitation transfer dynamics. Temperature promotes incoherent transport by destruction of coherence. In strong coupling limit when environment effect is negligible (Figure 4a), EET occurs via coherent mechanism and there is very small difference between temperature dependent and independent theories. However, in

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 22

long time limit after the transition from coherent to incoherent transport, the difference increases. Similarly when environment effect is very strong there is negligibly small effect of temperature on dynamics (Figure 4c) as environment contribution is quite larger than temperature. In the intermediate coupling regime when there is transition from coherent to incoherent transport the effect of temperature is prominent (Figures 4b and 4d). In this study, we use both temperature dependent and independent QSLE to investigate the EET dynamics of trimer and FMO complex. We consider fluctuation as a Poisson stochastic process for temperature independent studies. However, for temperature dependent studies we consider Gaussian bath. The main results of the present study are summarized as follows (1) The condition of coherence is that, the constant off-diagonal coupling J has to be larger than Vd. That is, to get incoherence, J must be overwhelmed. Thus, in extended systems with multiple chromophores, J may be significant, and coherence is natural. In nonMarkovian or strong coupling limit energy transport is fully coherent up to long time (~3 ps.). However, in the intermediate coupling regime we observe transition from coherent to incoherent excitation transport at 1 ps.. The asymmetry between the behavior Vd2 of both Vd and b as well as the ratio makes the study of EET dynamics non-trivial. bd −1 d

−1

In the Markovian and weak coupling limit, Vd and bd are combined together in the form Vd2 , the ratio dictates the nature EET dynamics and energy transfer occurs through bd incoherent or hopping mechanism. Excitation transfer in the intermediate regime can be termed as fluctuation or environment assisted coherent transport. In the absence of J (or,

14 ACS Paragon Plus Environment

Page 15 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

very small J) where one would expect incoherence, one can still get coherence in the −1 non-Markovian limit where Vd is large, bd is also large and baths are temporally

correlated. This is an important result of this work. (2) Absolute values of coherence in different limits indicate the non-local quantum correlation between different sites that lead to create new pathways for EET and consequently the EET dynamics is facilitated. (3) Presence of site energy heterogeneities and dynamic disorders creates competition between localization and delocalization of EET.

At intermediate time regime

heterogeneities in the site energies localizes the excitation energy whereas dynamic disorders help in delocalization of the same. At long time limit slow dynamics indicate the localization of excitation energy on each site. (4) Effects of temperature are investigated by the comparison between results obtained from temperature dependent and independent theories. In the strong coupling and in the nonMarkovian limit, the difference between both theories can be observed in the long time limit. As the temperature independent theory provides equal population in long time limit, difference between two theories increases in this limit. In the weak coupling or Markovian limit when environment effect is strong, both theories predict more or less same dynamics. However, in the intermediate regime when transition from coherent to incoherent transport occurs, temperature effect is more prominent. Experimentally it was observed that for FMO complex coherence can survive upto 1 ps. In our study Figures 2(b) and 3(b) show that coherence in the intermediate limit can survive upto 1 ps. which is similar to the experimental time scale of coherent EET in photosynthetic complex.

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 22

Sibey and co-workers studied efficiency of conversion and dependence on the trapping rate using extended Haken-Strobl model and they observed that efficiency is less than that obtained from full quantum calculation. They also concluded that EET efficiency in FMO complex is dominated by Förster hopping rate. To obtain greater efficiency, just coherent energy transfer is not enough. Presence of incoherent mechanism in long time limit after initial coherent dynamics helps in increasing the efficiency. To study the scenario in detail, one needs to consider both the effects of temperature and vibration in the dynamics in the non-Markovian limit.

ACKNOWLEDGEMENTS BB thanks Department of Science and Technology (DST, India) and Sir J. C. Bose fellowship for providing partial financial support. RD thanks Dr. Sarmistha Sarkar, Mr. Sayantan Mondal and Mr. Saumyak Mukherjee for carefully reading the manuscript.

Supporting Information Available Detailed description about temperature independent quantum stochastic Liouville equation and temperature dependent coupled equations of motion is provided.

REFERENCES 1. Anderson P. W. Absence of diffusion in certain random lattices, Phys. Rev. 1958, 109, 1492-1505. 2. Logan D. E.; Wolynes P. G. Anderson Localization in Topologically Disordered Systems, Phys. Rev. B 1985, 31, 2437-2450. 3. Zwanzig R., Rate Processes with Dynamical Disorder, Acc. Chem. Res. 1990, 23, 148152. 4. Engel G. S.; Calhoun T. R.; Read E. L.; Ahn T.-K.; Mančal T.; Cheng Y.-C.; Blankenship R. E.; Fleming G. R. Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems, Nature 2007, 446, 782-786.

16 ACS Paragon Plus Environment

Page 17 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

5. Ishizaki A.; Fleming G. R. Theoretical examination of quantum coherence in a photosynthetic system at physiological temperature, Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 17255–17260. 6. Ishizaki A.; Fleming G. R. Unified treatment of quantum coherent and incoherent hopping dynamics in electronic energy transfer: Reduced hierarchy equation approach, J. Chem. Phys. 2009, 130, 234111. 7. Ishizaki A.; Fleming G. R. On the adequacy of the Redfield equation and related approaches to the study of quantum dynamics in electronic energy transfer, J. Chem. Phys. 2009, 130, 234110. 8. Fassioli F.; Dinshaw R.; Arpin P. C.; Scholes G. D. Photosynthetic light harvesting: excitons and coherence, J. R. Soc. Interface 2014, 11, 20130901. 9. Scholes G. D. Quantum-Coherent Electronic Energy Transfer: Did Nature Think of It First? J. Phys. Chem. Lett., 2010, 1, 2–8. 10. Collini E.; Scholes G. D. Coherent Intrachain Energy Migration in a Conjugated Polymer at Room Temperature, Science 2009, 323, 369-373. 11. Collini E.; Scholes G. D. Electronic and Vibrational Coherences in Resonance Energy Transfer along MEH-PPV Chains at Room Temperature, J. Phys. Chem. A, 2009, 113, 4223–4241. 12. Sterpone F.; Martinazzo R.; Panda A. N; Burghardt I. Coherent Excitation Transfer Driven by Torsional Dynamics: a Model Hamiltonian for PPV Type Systems, Z. Phys. Chem. 2011, 225, 541–551. 13. Tozer O. R.; Barford W. Intrachain exciton dynamics in conjugated polymer chains in solution, J. Chem. Phys. 2015, 143, 084102. 14. Blankenship R. E. Molecular Mechanisms of Photosynthesis, Blackwell Science, Oxford/Malden, 2002. 15. Fenna R. E.; Matthews B. W. Chlorophyll arrangement in a bacteriochlorophyll protein from Chlorobium limicola, Nature 1975, 258, 573-577. 16. Jang S.; Cheng Y. C.; Reichman D. R.; Eaves J. D. Theory of coherent resonance energy transfer, J. Chem. Phys. 2008, 129, 101104. 17. Mohseni M.; Rebentrost P.; Lloyd S.; Aspuru-Guzik A. Environment-assisted quantum walks in photosynthetic energy transfer, J. Chem. Phys. 2008, 129, 174106. 18. Rebentrost P.; Mohseni M.; Kassal I.; Lloyd S.; Aspuru-Guzik A. Environment-assisted quantum transport, New J. Phys. 2009, 11, 033003. 19. Cao J.; Silbey R. J. Optimization of Exciton Trapping in Energy Transfer Processes, J. Phys. Chem. A 2009, 113, 13825-13838. 20. Wu J.; Liu F.; Shen Y.; Cao J.; Silbey R. J. Efficient energy transfer in light-harvesting systems, I: optimal temperature, reorganization energy and spatial–temporal correlations, New J. Phys. 2010, 12, 105012.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 22

21. Wu J.; Liu F.; Ma J.; Cao J; Silbey R. J. Efficient energy transfer in light-harvesting systems: Quantum-classical comparison, flux network, and robustness analysis, J. Chem. Phys. 2012, 137, 174111. 22. Chen X.; Silbey R. J. Effect of correlation of local fluctuations on exciton coherence, J. Chem. Phys. 2010, 132, 204503. 23. Chen X.; Silbey R. J. Excitation Energy Transfer in a Non-Markovian Dynamical Disordered Environment: Localization, Narrowing, and Transfer Efficiency, J. Phys. Chem. B 2011, 115, 5499–5509. 24. Dutta R.; Bagchi B. Effects of dynamic disorder on exciton migration: Quantum diffusion, coherences, and energy transfer, J. Chem. Phys. 2016, 145, 164907. 25. Dutta R.; Bagchi K.; Bagchi B. Role of quantum coherence in shaping the line shape of an exciton interacting with a spatially and temporally correlated bath, J. Chem. Phys. 2017, 146, 194902. 26. Moix J. M.; Khasin M.; Cao J. Coherent quantum transport in disordered systems: I. The influence of dephasing on the transport properties and absorption spectra on onedimensional systems, New J. Phys.2013, 15, 085010. 27. Kubo R. Stochastic Liouville Equations, J. Math. Phys. 1963, 4, 174-183. 28. Haken H.; Strobl G. An Exactly Solvable Model for Coherent and Incoherent Exciton Motion, Z. Phys. 1973, 262, 135-148. 29. Haken H.; Reineker P. The Coupled Coherent and Incoherent Motion of Excitons and its Influence on the Line Shape of Optical Absorption, Z. Phys. 1972, 249, 253-268. 30. Silbey R. Electronic Energy Transfer in Molecular Crystals, Annu. Rev. Phys. Chem. 1976, 27, 203-223. 31. Bagchi B.; Oxtoby D. W. Effect of excitation on non-Markovian vibrational energy relaxation, J. Phys. Chem., 1982, 86, 2197–2205. 32. Bagchi B.; Oxtoby D. W. Theory of non‐Markovian exciton transport in a one dimensional lattice, J. Chem. Phys.1983, 79, 6211-6221. 33. Bagchi B. Molecular Relaxation in Liquids, Oxford University press, New York, 2012. 34. Tanimra Y.; Kubo R. Time Evolution of a Quantum System in Contact with a Nearly Gaussian-Markoffian Noise Bath, J. Phys. Soc. Jpn. 1989, 58, 101-114. 35. Tanimura Y. Stochastic Liouville, Langevin, Fokker–Planck, and Master Equation Approaches to Quantum Dissipative Systems, J. Phys. Soc. Jpn. 2006, 75, 082001.

18 ACS Paragon Plus Environment

Page 19 of 22

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 20 of 22

Page 21 of 22

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 22 of 22