Theoretical examination of quantum coherence in a photosynthetic system at physiological temperature

Akihito Ishizaki and Graham R. Fleming
PNAS October 13, 2009 106 (41) 17255-17260; https://doi.org/10.1073/pnas.0908989106

  1. Contributed by Graham R. Fleming, August 7, 2009 (received for review June 4, 2009)

Abstract

The observation of long-lived electronic coherence in a photosynthetic pigment–protein complex, the Fenna–Matthews–Olson (FMO) complex, is suggestive that quantum coherence might play a significant role in achieving the remarkable efficiency of photosynthetic electronic energy transfer (EET), although the data were acquired at cryogenic temperature [Engel GS, et al. (2007) Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems. Nature 446:782–786]. In this paper, the spatial and temporal dynamics of EET through the FMO complex at physiological temperature are investigated theoretically. The numerical results reveal that quantum wave-like motion persists for several hundred femtoseconds even at physiological temperature, and suggest that the FMO complex may work as a rectifier for unidirectional energy flow from the peripheral light-harvesting antenna to the reaction center complex by taking advantage of quantum coherence and the energy landscape of pigments tuned by the protein scaffold. A potential role of quantum coherence is to overcome local energetic traps and aid efficient trapping of electronic energy by the pigments facing the reaction center complex.

Photosynthesis starts with the absorption of a photon of sunlight by one of the light-harvesting pigments, followed by transfer of the energy to the reaction center, where the primary electron transfer reactions convert the solar energy into an electrochemical gradient. (1) The transfer of this excitation energy towards the reaction center (RC) occurs with a near unity quantum yield.

The Fenna–Matthews–Olson (FMO) pigment–protein complex (2–5) is found in low light-adapted green sulfur bacteria. Under physiological conditions, this complex is situated between the so-called base-plate protein of the large peripheral chlorosome antenna and the RC complex, and is tasked with transporting sunlight energy harvested in the chlorosome to the RC pigments. The complex is a trimer made of identical subunits, each of which contains seven bacteriochlorophyll (BChl) molecules. By virtue of its relatively small size, it represents an important model in photosynthetic excitation energy transfer (EET), and has been extensively studied experimentally and theoretically. In the late 90s, Savikhin et al. (6) observed quantum beating in the FMO complex using fluorescence anisotropy technique although it had been generally assumed that electronic coherence decay so rapidly that it does not affect the EET. (7) However, it was difficult to unambiguously assign this experimental observation to an electronic origin. Recently, Engel et al. (8) investigated the FMO complex isolated from Chlorobaculum tepidum by means of 2D electronic spectroscopy, (9–11) and succeeded in observing long-lasting quantum beating providing direct evidence for long-lived electronic coherence. (12, 13) The observed coherence lasts for time scales similar to the EET timescales, implying that electronic excitations move coherently through the FMO protein rather than by incoherent hopping motion. Observation of long-lasting quantum coherence is not unique to the FMO complex. Lee et al. (14) revealed coherent dynamics in the RC of a purple bacteria, Rhodobacter sphaeroides, by applying two-color electronic coherent photon echo technique. These observations have led to the suggestion that quantum coherence might play significant roles in achieving remarkable efficiency of photosynthetic EET. Although the observations of long-lasting quantum coherence provide valuable insights into the inner working of photosynthetic complexes, the measurements were performed outside the physiological range of temperatures. Generally, quantum coherence at physiological temperatures is fragile compared to that at cryogenic temperatures because amplitude of environmental fluctuations increases with increasing temperature. Hence, the robustness and roles of quantum coherence under physiological conditions are to a large extent unknown. In order to explore these questions, theoretical investigations with models of appropriate sophistication are required in addition to further experimental studies (15–18). Recently, Aspuru-Guzik and coworkers (19–21), and Plenio and Huelga (22) discussed effects of quantum coherence on enhancement of photosynthetic EET efficiency from the perspective of a quantum walk.

Photosynthetic EET processes are usually analyzed in two perturbative limits. When the electronic coupling between pigments is small in comparison to the electron-environment coupling, which can be specified by magnitude of reorganization energy, the electronic coupling can be treated perturbatively. This treatment yields Förster theory (23), which describes incoherent hopping between pigments. In the opposite limit, when the reorganization energies are small, it is possible to treat the electron-environment coupling perturbatively to obtain a quantum master equation. The most commonly used theory from this approach in the literature of photosynthetic EET is the conventional Redfield equation (24), which is one of the few viable paths to explore quantum coherence. However, in a typical photosynthetic EET system the reorganization energies are not small in comparison to the electronic coupling. In the FMO complex, the electronic coupling strengths span a wide range, 1 ≈ 100 cm−1 while the suggested reorganization energies span a similar range (11, 25–27). Hence, the Redfield equation approach might lead to erroneous insights or incorrect conclusions regarding the quantum coherence and its interplay with the protein environment (17, 18). Recently, the present authors developed a quantum dynamics equation for EET (18) by taking into account effects of the environmental reorganization dynamics on the pigments in a nonperturbative fashion based on the hierarchical expansion technique (28). The developed equation can describe EET processes irrespective of the electron-environment coupling strength, and reduces to the conventional Redfield and Förster theories in their respective limits of validity. In the regime of coherent wave-like motion, moreover, the equation predicts several times longer lifetimes of the electronic coherence between pigments than does the Redfield equation.

In this article, we apply the aforementioned quantum dynamic equation in order to investigate the robustness and roles of quantum coherence in the C. tepidum FMO complex under physiological conditions. Here, special attention is paid to the EET dynamics in the site basis in order to obtain insights into the spatial and temporal evolution of electronic excitation energy through the FMO complex. Calculations are performed on the basis of the electron-environment coupling parameters adopted for simultaneous fitting of linear and 2D rephasing, nonrephasing and polarization-dependent electronic spectra (27). The numerical results reveal that electronic coherence persists for several hundred femtoseconds even at physiological temperature, T = 300 K. Based on the numerical results, the physical roles of quantum coherence under physiological conditions are discussed.

Theory

To describe EET dynamics in a photosynthetic complex containing N pigments, each pigment is modeled by a two-level system describing its S0 → S1 transition (the Qy transition of a BChl). The total Hamiltonian of the EET dynamics consists of three parts, Embedded Image The first part Ĥe is the Hamiltonian of electronic states of the pigments expressed as Embedded Image where |j〉 represents the state where only the jth pigment is in its S1 state and all others are in their S0 states. ɛj is the so-called site energy of the jth pigment defined as the optical transition energy at the equilibrium configuration of environmental phonons associated with the S0 state. Jjk denotes the electronic coupling between the jth and kth pigments. The second term Ĥph is the Hamiltonian of the environmental phonons, where the ξth mode dynamics is described by the dimensionless coordinate, q̂ξ. The last Hamiltonian Ĥel-ph is responsible for fluctuations in the site energies by the phonon dynamics, and is expressed as Embedded Image where we have defined V̂j ≡ |j〉〈j| and ûj ≡ −∑ξcjξq^ξ with cjξ being the coupling constant between the jth pigment and ξth phonon mode. The fluctuations in site energies of different pigments may be correlated or uncorrelated, depending on whether the site energies couple to the same or to different phonon modes of the protein environment. At present, however, there is no direct way to access accurate information on the correlated fluctuations in the FMO complex; hence, we assume uncorrelated fluctuations in this work.

An adequate description of the EET dynamics is given by the reduced density operator ρ^(t) ≡ Trph{ρ^tot(t)}, where ρ^tot(t) denotes the density operator for the total system. For the reduction we suppose that the total system at the initial time, t = 0, is in the factorized product state of the form, ρ^tot(0) ∝ ρ^(0) exp(−βH^ph) with β ≡ 1/kBT. Generally, this initial condition is unphysical since it neglects an inherent correlation between a system and its environment (28, 29). In electronic excitation processes, however, this initial condition is appropriate because it corresponds to excitation in accordance to the vertical Franck–Condon transition (18).

The coupling of an electronic transition of the jth pigment to the environmental phonons can be specified by the spectral density of the electron-phonon coupling constants, Jj(ω). This includes information contents on dynamics of the environmental phonons. The equilibrium phonon state in the S0 state is generally a nonequilibrium state in the S1 state because of the electron-phonon coupling. After the excitation, hence, the coupling relaxes the phonons toward the actual equilibrium with dissipating the reorganization energy. The reorganization dynamics, in principle, can be measured by fluorescence Stokes shift experiment (30), where the observable quantity is the phonon relaxation function, Embedded Image Note that Γj(0) ≡ 2λj is the Stokes shift magnitude expressed as two times the reorganization energy λj. Simultaneously, the coupling induces fluctuations in the site energies of pigments. Information on the timescale of the processes can be obtained by means of three-pulse photon-echo peak shift experiment (30), where the observable quantity is the phonon symmetrized correlation function, Sj(t). If the environmental phonons can be described classically, the relaxation function and the symmetrized correlation function satisfy the classical fluctuation-dissipation relation, Γj(t) = βSj(t) ; hence, these two functions contain the same information on the phonon dynamics, whose characteristic timescale is given by ref. 31Embedded Image Generally, the relaxation function has a complicated form involving several components (30). In this article, we model the relaxation function by an exponential form, Γj(t) = 2λjexp(−γjt), in order to focus on the timescale of the phonon relaxation dynamics. For this modeling, the timescale of the phonon relaxation is simply τc = γj−1, and the spectral density is expressed as the overdamped Brownian oscillator model, Jj(ω) = 2λjγjω/(ω2 + γj2). Although this spectral density has been successfully used for analyses of experimental results (27, 32, 33), it may produce qualitatively different phonon sidebands from the experimental results in zero temperature limit (7, 34, 35). However, this limitation is not the major concern of the present study.

When the high-temperature condition characterized by βℏγj < 1 is assumed for the overdamped Brownian oscillator model, the following hierarchically coupled equations of motion for the reduced density operator can be derived in a nonperturbative manner (18): Embedded Image for sets of nonnegative integers, n ≡ (n1,n2,…,nN). nj± differs from n only by changing the specified nj to nj ± 1, i.e., nj± ≡ (n1,…,nj ± 1,…,nN). In Eq. 6, only the element σ^(0,t) is identical to the reduced density operator ρ^(t), while the others {σ^(n ≠ 0,t)} are auxiliary operators to take into account fluctuations in site energies and dissipation of reorganization energies. The Liouvillian corresponding to the electronic Hamiltonian H^e is denoted by L^e, while the relaxation operators are given by Embedded Image Embedded Image where we have introduced the hyper-operator notations, O^ ×f^ ≡ O^f^ − f^O^ and O^ ○ f^ ≡ O^f^ + f^O^, for any operator O^ and operand f^. Note that Eq. 6 is a set of equations of motion in the operator form, and hence the numerical results are independent of the representation in which the equations are integrated. The hierarchically coupled equations Eq. 6 continue to infinity, which is impossible to treat computationally. In order to terminate Eq. 6 at a finite stage, we replace Eq. 6 by Embedded Image for the integers n = (n1,n2,…,nN) satisfying Embedded Image where ωe is a characteristic frequency for L^e (18). Thus, the required number of the operators {σ^(n,t)} is evaluated as ∑k=0N(N−1k+N−1)=(N+N)!/(N!N!) . Note that low-temperature correction terms explained in Appendix should be included into Eq. 6 when the high-temperature condition βℏγj < 1 is not satisfied.

As demonstrated in ref. 18, the reduced hierarchy equation Eq. 6 can describe quantum coherent wave-like motion, incoherent hopping, and an elusive intermediate EET regime in a unified manner, and reduces to the conventional Redfield and Förster theories in their respective limits of validity. Recently, Jang et al. (15) developed another quantum dynamic equation to interpolate between the Redfield and Förster limits by employing the small polaron transformation (36). The small polaron approach is based on the second-order perturbative truncation with respect to the renormalized electron-phonon coupling, whereas the present hierarchy approach is derived in a non perturbative fashion. It will be interesting to compare the two approaches for future works.

Results and Discussion

In this section, we present and discuss numerical results regarding EET dynamics in the FMO complex of C. tepidum. The complex is a trimer made of identical subunits, each containing seven BChls. Because the strongest electronic coupling between two BChls in different FMO monomeric subunit is about an order of magnitude smaller than the local reorganization energies, the coherence between them is rapidly destroyed by the environmental disturbance (26). Therefore, we assume that the intersubunit coupling is vanishingly small and we consider the EET dynamics within one subunit. To simulate the EET dynamics, we use the Hamiltonian for the trimeric structure of the FMO complex given in ref. 26 (supporting information (SI) Table S1.) We use the usual numbering of the BChls, which was originally chosen by Fenna and Matthews (2) (see Fig. 1).

Fig. 1.
Fig. 1.

Seven BChl molecules belonging to the monomeric subunit of the FMO complex. The complex is oriented with BChl 1 and 6 toward the baseplate protein whereas BChl 3 and 4 define the target region in contact with the reaction center complex. The spiral strands are α-helices that are part of protein environment.

Although the spectral density Jj(ω) is a crucial factor to describe EET dynamics, no direct and detailed information on its form is available for the FMO complex at present. Hence, several empirical forms have been employed under the assumption that the spectral densities for the different BChls are equivalent (11, 25–27, 33), and then the phonon relaxation time τc estimated by Eqs. 4 and 5 ranges from 35 fs (26)* to 166 fs (25) in the literature. Recently, Read et al. (27) conducted 2D electronic spectroscopic experiments to visualize excitonic structure in the FMO complex of Prosthecochloris aestuarii, and they performed simultaneous fitting of the linear and 2D rephasing, nonrephasing, and polarization-dependent spectra by employing the overdamped Brownian oscillator model. To obtain excellent agreement between the experimental data and numerical fitting, they adopted λj = 35 cm−1 and τc = γj−1 = 50 fs as the values of reorganization energy and relaxation time of the phonons, respectively. Therefore, we also employ these values with the assumption that the phonon spectral densities for the individual BChls are equivalent. For numerical integration of Eq. 6, the depth of hierarchy we employed here is at most N = 12. For all calculations, the accuracy were checked by changing the values of N to make sure that the numerical results are converged.

EET Dynamics and Temperature Dependence

Following the previous proposal based on theoretical calculations (26), the FMO complex has been assumed to be oriented with BChl 1 and 6 toward the baseplate protein whereas BChls 3 and 4 define the target region in contact with the RC complex. Recently, this orientation was verified experimentally by Wen et al. (5). Accordingly, we adopt BChls 1 or 6 as the initial excited pigment for numerical calculations.

Fig. 2 presents the EET dynamics at cryogenic temperature, T = 77 K. These results clearly show that the energy flow in the FMO complex occurs primarily through two EET pathways, which connect spatially proximate and excitonically coupled BChls as demonstrated by Brixner et al. with 2D electronic spectroscopy (11, 25): Embedded Image Embedded Image

Fig. 2.
Fig. 2.

Time evolution of the population of each BChl in the FMO complex. Calculations were done for cryogenic temperature, T = 77 K. The reorganization energy and the phonon relaxation time are set to be λj = 35 cm−1 and τc = γj−1 = 50 fs, respectively.

The double arrow (⇌) in Eq. 11a indicates that excitation energy equilibrates between BChls 3 and 4 after BChl 3 is populated. In Fig. 2A, the long lifetime of BChl 1 is due to the relatively high site energy of BChl 2. Hence, populations of BChls 1 and 3 are dominant at time t = 1 ps. In Fig. 2B, on the other hand, those of BChls 3 and 4 are dominant at time t = 1 ps owing to very fast transport of electronic excitation energy. These behaviors are consistent with the transient pump-probe experiment by Vulto et al. (7, 37) because BChl 1 participates in exciton 3 and BChls 3 and 4 construct exciton 1, where the numbers give the rank starting with the smallest exciton energy. Furthermore, in Fig. 2A and B, quantum coherent wave-like motions are visible up to 700 fs. This time scale is consistent with the observation by Engel et al. (8) with 2D electronic spectroscopy.

Fig. 3 gives the EET dynamics at physiological temperature, T = 300 K. The dynamics at T = 300 K is also dominated by the two EET pathways given in Eq. 11, as described by Adolphs and Renger with the conventional and modified Redfield theories. (26, 28) More noteworthy is that coherent wave-like motions can be clearly observed up to 350 fs even at the physiological temperature. Note that the phonon relaxation time used here, τc = 50 fs, was obtained as a numerical fitting parameter for the experimental data; it was not measured directly. For numerical fitting of other 2D spectra of the FMO complex from C. tepidum (11), an Ohmic spectral function with an exponential cutoff was employed, i.e., Jj(ω) ∝ λjωexp(−ω/ωc) with λj = 35 cm−1 and ωc = 50 cm−1 (25), which yields the phonon relaxation time of τc = 166 fs. Therefore, attention should be also paid to possibilities of slower phonon relaxation times. Most recently, the present authors have demonstrated that slower phonon relaxation times are preferable in terms of longer lifetime of electronic coherence among pigments (18). Indeed, calculations employing the slower value, γj−1 = 166 fs, predict longer-lasting (up to 550 fs) wave-like motions at T = 300 K as shown in Fig. 4, which can not be reproduced with conventional Redfield theory because the dynamics are now in the strong non-Markovian regime. In Eq. 6, slower fluctuations preserve longer-lived quantum coherent oscillation, whereas the Markov approximation in the Redfield theory causes infinitely fast fluctuation and thus collapses the coherence quickly. On the other hand, in the fastest relaxation case in the literature, τc = 35 fs (35), numerical results with Eq. 6 also show the wave-like behavior lasting for 350fs at T = 300 K, virtually the same as in the case of τc = 50 fs. (Fig. S1) Therefore, we expect that long-lived electronic coherence among BChls is present in the FMO complex even at physiological temperatures, irrespective of the phonon relaxation times.

Fig. 3.
Fig. 3.

Time evolution of the population of each BChl in the FMO complex. Calculations were done for physiological temperature, T = 300 K. The other parameters are the same as in Fig. 2.

Fig. 4.
Fig. 4.

Time evolution of the population of each BChl in the FMO complex. Calculations were done for the same parameters as in Fig. 3, except for the phonon relaxation times τc = γj−1.

Roles of Quantum Coherence

The observation of long-lasting and robust quantum coherence prompts the speculation that quantum effects may play a significant role in achieving the remarkable efficiency of photosynthetic EET. Engel et al. (8) proposed that the FMO complex performs a quantum search algorithm that is more efficient than a classical random walk suggested by the hopping mechanism. Quantum coherence enables the excitation to rapidly and reversibly sample multiple pathways to search for BChl 3 that connects to the RC.

Unidirectional energy flow facilitates the achievement of a near-unity quantum yield. For the unidirectionality, once excitation migrates to the linker pigments between the FMO complex and the RC, it has to be trapped and directed to the RC (27). Here, we explore the idea that the FMO complex may work as a type of “rectifier” for the unidirectional flow from the chlorosome antenna to the RC by taking advantage of quantum coherence and the site energy landscape tuned by the protein scaffold. The basic idea is as follows: If the EET were dominated only by diffusive incoherent hopping, trapping in subsidiary energetic minima would be inevitable. However, quantum delocalization can allow avoidance of the traps to aid the subsequent trapping of excitation by the linker pigments, BChls 3 and 4.

More insight can be gained if attention is given to the site energy of the BChl in the baseplate. It has been reported that the Qy absorption of the baseplate BChl in C. tepidum is located around 800 nm (12,500 cm−1) (39, 40). This lies between the site energies of BChl 1 (12,410 cm−1) and BChl 6 (12,630 cm−1) in ref. 26, and implies that the FMO complex has two windows and their respective onward pathways given in Eq. 11 in order to capture and convey electronic excitation energy from the baseplate without loss. Fig. 5 presents the energy landscapes along the two primary pathways.

Fig. 5.
Fig. 5.

The energy landscapes along the two primary transfer pathways in the FMO complex: baseplate → BChls 1 → 2 → 3 → 4 (A) and baseplate → BChls 6 → 5, 7, 4 → 3 (B). The relatively strong couplings between BChls are depicted by solid lines.

In terms of an “energy funnel,” it is advantageous that the site energy of BChl 1 is low in comparison to that of the baseplate BChl. Nevertheless, the site energy of BChl 2 is higher than BChl 1, and hence BChl 1 is a local energetic minimum. It is not clear whether there exists a physiological significance to the barrier created by BChl 2. Here, we notice that the energy gap between BChls 1 and 3 is 200 cm−1, which corresponds to the magnitude of thermal energy kBT at physiological temperature T = 300 K. If the BChls in the pathway shown in Fig. 5A were arranged in such a way that their site energies monotonically decrease, the relatively flat energy landscape compared to kBT would allow facile backward transfer of excitation away from BChl 3 at physiological temperatures, and the equilibrated population of BChl 3 would decrease. However, the FMO energy landscape is set up such that the site energy of BChl 2 is higher by ≈300 cm−1 than that of BChl 3; hence, the back transfer is suppressed and thus BChl 3 is able to have larger amount of population of excitation in equilibrium. (Fig. S2) If the energy transfer from BChl 1 to BChl 2 were a slow incoherent hopping process, the excitation energy would be trapped on BChl 1. However, the wave-like motion induced by the quantum superposition between BChls 1 and 2 enables such trapping to be avoided. Indeed, Figs. 2A and 3A clearly show that the detrapping of excitation on BChl 1 involves quantum coherent wave-like motion. (Figs. S3 and S4)

The site energy of BChl 6 is higher than the energy of the baseplate BChl. If the deexcitation of the BChl 6 were slow, the excitation could return to the baseplate easily according to the detailed balance condition between BChl 6 and a baseplate BChl. In order to avoid this backward transfer, BChl 6 needs to convey its excitation energy to other BChls as quickly as possible. This is accomplished by having BChl 6 be strongly coupled to both of BChls 5 and 7, which in turn are strongly coupled to BChl 4. Figs. 2, 3, and 4 demonstrate that the initial rises of populations in BChls 5, 7, and 4 are based on quantum wave-like motion irrespective of temperature, and then the excitation of BChl 6 is delocalized coherently over BChls 4, 5, 6, and 7 on the 100-fs time scale. After the ultrafast delocalization, the excitation begins to explore the lowest energy site, BChl 3. This enables the energy flow though the pathway shown in Fig. 5B to be unidirectional and highly efficient.

Here, we have not taken into account the effect of static disorder because the suggested standard deviation of the Gaussian distribution of the site energies is ≈20 cm−1 (the corresponding fwhm is ≈50 cm−1) (25, 27), which is smaller than the site energy gaps in the FMO complex, and therefore the static disorder of the individual sites has little influence on the qualitative nature of the energy landscape. However, even if the static disorder is relatively large, quantum coherence is still expected to help excitation to overcome local energetic traps created by the static disorder as discussed by Renger and Marcus (41) for the variable-range hopping mechanism of electron transfer through DNA.

Concluding Remarks

In this article, we have investigated theoretically the spatial and temporal dynamics of the EET through the FMO complex of C. tepidum in order to address the robustness and role of the quantum coherence under physiological conditions. For the description of the dynamics, we applied the hierarchically coupled equations approach, which reduces to the conventional Redfield and Förster theories in their respective limits of validity. The numerical calculations were performed based on the electron-phonon coupling parameters adopted for the fitting of 2D electronic spectroscopic data, where the reorganization energy and the phonon relaxation time are 35 cm−1 and τc = 50 fs, respectively. The results revealed the presence of quantum coherent wave-like motion persisting for 350 fs even at physiological temperature T = 300 K. Taking into account the possibility of slower phonon relaxation, we also carried out numerical calculations for τc = 166 fs. For this case, the dynamic equation used here yields wave-like motion persisting for 550 fs at T = 300 K, which cannot be reproduced by the conventional Markovian Redfield equation because the dynamics are in the strong non-Markovian regime. Recently, Collini and Scholes (42, 43) showed the presence of a quantum coherent EET process in the conjugated polymer MEH-PPV by means of 2D electronic spectroscopy. Their 2D spectra show that long-lived electronic coherences persist for at least 250 fs after photoexcitation of MEH-PPV in solution at 293 K. Although MEH-PPV is not a photosynthetic complex, their experimental observation of long-lived coherence at room temperature is corroborative for our theoretical predictions. Moreover, as one of potential roles of quantum coherence in the FMO complex, we suggested that quantum coherence allows the FMO complex to work as a rectifier for unidirectional energy flow from the chlorosome antenna to the RC complex. It is reasonable to suppose that quantum coherence can overcome local energetic traps to aid the subsequent trapping of electronic energy by the BChl molecules in contact with the RC complex.

In the present calculations, it has been assumed that the phonon spectral densities for the individual BChls are equivalent and fluctuations in different site energies are uncorrelated. This is because there exists no accurate and detailed information on the coupling between BChls and phonons of the FMO protein. Nevertheless, our results predict the presence of quantum coherence lasting for several hundred femtoseconds at physiological temperature. In the RC complex of a purple bacterium, Groot et al. (44) reported that the phonon relaxation time and electron-phonon coupling strength of one of the BChls are both about 1.5 times different from those of a record BChl in the same complex. Additionally, Lee et al. (14) showed that fluctuations in the bacteriopheophytin and accessory BChl excited energies are strongly correlated enabling long-lasting electronic coherence between them. It has also been reported that correlated fluctuations may enhance the contribution of quantum coherence to EET efficiency of the FMO complex (20). Hence, it is also likely that the couplings between the FMO protein and the embedded BChls are heterogeneous and may be optimized for individual BChl sites in order to preserve longer-lived coherence or to achieve higher efficiency in comparison to the present results. A more detailed understanding of quantum coherent effects in photosynthetic systems demands comprehensive investigations combining the quantum dynamic theories with quantum chemical calculations and molecular dynamics simulations, in addition to further experimental studies. Investigations from the perspective of quantum information science should also provide insights into the issue.

Appendix: Low-Temperature Correction Terms

For the overdamped Brownian oscillator model of the phonon spectral density, the phonon symmetrized correlation function can be expressed as Embedded Image where the first term is the classical correlation function and the second term involving the bosonic Matsubara frequency, νk ≡ 2πk/βℏ (k ≥ 1), represents the quantum corrections. Under the high-temperature condition characterized by βℏγj < 1, the second term is vanishingly small and then the symmetrized correlation function can be described classically. For low temperature βℏγj > 1, however, the first term is not sufficient. The quantum dynamic equation based only on the first term does not ensure detailed balance at low temperature; hence, some diagonal elements of the density matrices may be negative.

In this article, we have included only the lowest-order quantum correction involving only νk=1 to calculate the EET dynamics at T = 77K, because the breakdown of the detailed balance condition is minor for the parameters employed here. For the cases of γj−1 ≥ 50 fs and T ≥ 77 K, the value of ν1 is sufficiently large in comparison with γj ; hence, ν1exp(−ν1t) can be replaced by Dirac's δ function. As the result, the symmetrized correlation function containing the lowest-order quantum correlation is given by Embedded Image Following the derivation in ref. 18 with Eq. A2, we can obtain a modified equation applicable to the case of T = 77 K. The resultant expression can be obtained by the following replacements in Eq. 6: Embedded Image Embedded Image

Acknowledgments

We thank Yuan-Chung Cheng, Tae Kyu Ahn, Thomas Renger, and Alán Aspuru-Guzik for insightful discussions and critical reading of the manuscript. This work was supported by the Director, Office of Science, Office of Basic Energy Sciences, of the U.S. Department of Energy under Contract. DE-AC02-05CH11231 and by the Chemical Sciences, Geosciences and Biosciences Division, Office of Basic Energy Sciences, U.S. Department of Energy under contract DE-AC03-76SF000098. A.I. was supported by Postdoctoral Fellowship for Research Abroad by the Japan Society for the Promotion of Science.

Footnotes

  • 1To whom correspondence should be addressed. E-mail: grfleming{at}lbl.gov.

  • Author contributions: G.F. and A.I. designed research; A.I. performed research; and G.F. and A.I. wrote the paper.

  • The authors declare no conflict of interest.

  • See Commentary on page 17247.

  • This article contains supporting information online at www.pnas.org/cgi/content/full/0908989106/DCSupplemental.

  • ↵* The relaxation function and the symmetrized correlation function at 300 K produced by the spectral density in ref. 26 show nonoscillatory negative values from 80 fs to 1 ps, whose physical origin is not yet clear. Hence, it is impossible to estimate the relaxation time τc by employing Eq. 5. Here, we have estimated the relaxation time only from the initial decay (<80 fs) of the relaxation function.

References

  1. ↵
    1. Blankenship RE
    (2002) Molecular Mechanisms of Photosynthesis (World Scientific, London).
  2. ↵
    1. Fenna RE,
    2. Matthews BW
    (1975) Chlorophyll arrangement in a bacteriochlorophyll protein from Chlorobium limicola. Nature 258:573–577.
    OpenUrlCrossRef
  3. ↵
    1. Li Y-F,
    2. Zhou W,
    3. Blankenship RE,
    4. Allen JP
    (1997) Crystal structure of the bacteriochlorophyll a protein from Chlorobium tepidum. J Mol Biol 271:456–471.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Camara-Artigas A,
    2. Blankenship RE,
    3. Allen JP
    (2003) The structure of the FMO protein from Chlorobium tepidumat 2.2Å resolution. Photosynth Res 75:49–55.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Wen J,
    2. Zhang H,
    3. Gross ML,
    4. Blankenship RE
    (2009) Membrane orientation of the FMO antenna protein from Chlorobaculum tepidumas determined by mass spectrometry-based footprinting. Proc Natl Acad Sci USA 106:6134–6139.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Savikhin S,
    2. Buck DR,
    3. Struve WS
    (1997) Oscillating anisotropies in a bacteriochlorophyll protein: Evidence for quantum beating between exciton levels. Chem Phys 223:303–312.
    OpenUrlCrossRef
  7. ↵
    1. Vulto SIE,
    2. et al.
    (1999) Excited state dynamics in FMO antenna complexes from photosynthetic green sulfur bacteria: A kinetic model. J Phys Chem B 103:8153–8161.
    OpenUrl
  8. ↵
    1. Engel GS,
    2. et al.
    (2007) Evidence for wavelike energy transfer through quantum coherence in photosynthetic systems. Nature 446:782–786.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Jonas DM
    (2003) Two-dimensional femtosecond spectroscopy. Annu Rev Phys Chem 54:425–463.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Brixner T,
    2. Mancăl T,
    3. Stiopkin IV,
    4. Fleming GR
    (2004) Phase-stabilized two-dimensional electronic spectroscopy. J Chem Phys 121:4221–4236.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Brixner T,
    2. et al.
    (2005) Two-dimensional spectroscopy of electronic couplings in photosynthesis. Nature 434:625–628.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Pisliakov AV,
    2. Mancăl T,
    3. Fleming GR
    (2006) Two-dimensional optical three-pulse photon echo spectroscopy. II. Signatures of coherent electronic motion and exciton population transfer in dimer two-dimensional spectra. J Chem Phys 124:234505.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Cheng Y-C,
    2. Fleming GR
    (2008) Coherence quantum beats in two-dimensional electronic spectroscopy. J Phys Chem A 112:4254–4260.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Lee H,
    2. Cheng Y-C,
    3. Fleming GR
    (2007) Coherence dynamics in photosynthesis: Protein protection of excitonic coherence. Science 316:1462–1465.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Jang S,
    2. Cheng Y-C,
    3. Reichman DR,
    4. Eaves JD
    (2008) Theory of coherent resonance energy transfer. J Chem Phys 129:101104.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Cheng Y-C,
    2. Fleming GR
    (2009) Dynamics of light harvesting in photosynthesis. Annu Rev Phys Chem 60:241–262.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Ishizaki A,
    2. Fleming GR
    (2009) On the adequacy of the Redfield equation and related approaches to the study of quantum dynamics in electronic energy transfer. J Chem Phys 130:234110.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Ishizaki A,
    2. Fleming GR
    (2009) Unified treatment of quantum coherent and incoherent hopping dynamics in electronic energy transfer: Reduced hierarchy equations approach. J Chem Phys 130:234111.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Mohseni M,
    2. Rebentrost P,
    3. Lloyd S,
    4. Aspuru-Guzik A
    (2008) Environment-assisted quantum walks in photosynthetic energy transfer. J Chem Phys 129:174106.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Rebentrost P,
    2. Mohseni M,
    3. Aspuru-Guzik A
    (2009) Role of quantum coherence in chromophoric energy transport. J Phys Chem B 113:9942–9947.
    OpenUrlPubMed
  21. ↵
    1. Rebentrost P,
    2. Mohseni M,
    3. Kassal I,
    4. Lloyd S,
    5. Aspuru-Guzik A
    (2009) Environment-assisted quantum transport. New J Phys 11:033003.
    OpenUrlCrossRef
  22. ↵
    1. Plenio MB,
    2. Huelga SF
    (2008) Dephasing-assisted transport: quantum networks and biomolecules. New J Phys 10:113019.
    OpenUrlCrossRef
  23. ↵
    1. Förster T
    (1948) Zwischenmolekulare energiewanderung und fluoreszenz. Ann Phys 2:55–75.
    OpenUrl
  24. ↵
    1. Redfield AG
    (1957) On the theory of relaxation processes. IBM J Res Dev 1:19–31.
    OpenUrl
  25. ↵
    1. Cho M,
    2. Vaswani HM,
    3. Brixner T,
    4. Stenger J,
    5. Fleming GR
    (2005) Exciton analysis in 2D electronic spectroscopy. J Phys Chem B 109:10542–10556.
    OpenUrlPubMed
  26. ↵
    1. Adolphs J,
    2. Renger T
    (2006) How proteins trigger excitation energy transfer in the FMO complex of green sulfur bacteria. Biophys J 91:2778–2797.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Read EL,
    2. et al.
    (2008) Visualization of excitonic structure in the Fenna-Matthews-Olson photosynthetic complex by polarization-dependent two-dimensional electronic spectroscopy. Biophys J 95:847–856.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Tanimura Y
    (2006) Stochastic Lioville, Langevin, Fokker-Planck, and master equation approaches to quantum dissipative systems. J Phys Soc Jpn 75:082001.
    OpenUrlCrossRef
  29. ↵
    1. Grabert H,
    2. Schramm P,
    3. Ingold G-L
    (1988) Quantum Brownian motion: The functional integral approach. Phys Rep 168:115–207.
    OpenUrlCrossRef
  30. ↵
    1. Fleming GR,
    2. Cho M
    (1996) Chromophore-solvent dynamics. Annu Rev Phys Chem 47:109–134.
    OpenUrlCrossRef
  31. ↵
    1. Kubo R
    (1963) Stochastic Liouville equations. J Math Phys 4:174–183.
    OpenUrlCrossRef
  32. ↵
    1. Zigmantas D,
    2. et al.
    (2006) Two-dimensional electronic spectroscopy of the B800-B820 light-harvesting complex. Proc Natl Acad Sci USA 103:12672–12677.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Read EL,
    2. et al.
    (2007) Cross-peak-specific two-dimensional electronic spectroscopy. Proc Natl Acad Sci USA 104:14203–14208.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Wendling M,
    2. et al.
    (2000) Electron-Vibrational coupling in the Fenna-Matthews-Olson complex of Prosthecochloris aestuariidetermined by temperature-dependent absorption and fluorescence line-narrowing measurements. J Phys Chem B 104:5825–5831.
    OpenUrl
  35. ↵
    1. Renger T,
    2. Marcus RA
    (2002) On the relation of protein dynamics and exciton relaxation in pigment-protein complexes: An estimation of the spectral density and a theory for the calculation of optical spectra. J Chem Phys 116:9997–10019.
    OpenUrlCrossRef
  36. ↵
    1. Rackovsky S,
    2. Silbey R
    (1973) Electronic-energy transfer in impure solids. I. Two molecules embedded in a lattice. Mol Phys 25:61–72.
    OpenUrlCrossRef
  37. ↵
    1. Vulto SIE,
    2. et al.
    (1998) Excited-state structure and dynamics in FMO antenna complexes from photosynthetic green sulfur bacteria. J Phys Chem B 102:10630–10635.
    OpenUrl
    1. Müh F,
    2. et al.
    (2007) α-Helices direct excitation energy flow in the Fenna-Matthews-Olson protein. Proc Natl Acad Sci USA 104:16862–16867.
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Francke C,
    2. Amesz J
    (1997) Isolation and pigment composition of the antenna system of four species of green sulfur bacteria. Photosynth Res 52:137–146.
    OpenUrlCrossRef
  40. ↵
    1. Frigaard N-U,
    2. et al.
    (2005) Isolation and characterization of carotenosomes from a bacteriochlorophyll c-less mutant of Chlorobium tepidum. Photosynth Res 86:101–111.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Renger T,
    2. Marcus RA
    (2003) Variable-range hopping electron transfer through disordered bridge states: Application to DNA. J Phys Chem A 107:8404–8419.
    OpenUrlCrossRef
  42. ↵
    1. Collini E,
    2. Scholes GD
    (2009) Coherent intrachain energy migration in a conjugated polymer at room temperature. Science 323:369–373.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Collini E,
    2. Scholes GD
    (2009) Electronic and vibrational coherences in resonance energy transfer along MEH-PPV chains at room temperature. J Phys Chem A 113:4223–4241.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Groot ML,
    2. Yu JY,
    3. Agarwal R,
    4. Norris JR,
    5. Fleming GR
    (1998) Three-pulse photon echo measurements on the accessory pigments in the reaction center of Rhodobacter sphaeroides. J Phys Chem B 102:5923–5931.
    OpenUrl
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Theoretical examination of quantum coherence in a photosynthetic system at physiological temperature
Akihito Ishizaki, Graham R. Fleming
Proceedings of the National Academy of Sciences Oct 2009, 106 (41) 17255-17260; DOI: 10.1073/pnas.0908989106
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

More Articles of This Classification

Physical Sciences

Show more

Chemistry

Show more

Biological Sciences

Show more

Biophysics and Computational Biology

Show more

Related Content

Cited by...

Similar Articles

You May Also be Interested in

PNAS Profile with NAS member and anthropologist Michael Tomasello
PNAS Profile
PNAS Profile with NAS member and anthropologist Michael Tomasello
Early monumental burial sites
Researchers report an early monumental burial site near Lake Turkana in Kenya that may have served as a stable landmark for mobile herders in a changing physical environment and as a social anchor point to foster communal identity and interaction among mobile herders.
Moon. Image courtesy of Pixabay/Ponciano.
Evidence of surface water ice on the moon
A study reports evidence of water ice on the moon’s surface, discerned via a signature in the near-infrared reflectance spectra that suggests the ice was formed by slow condensation due to impact or water migration through the lunar exosphere.
Image courtesy of Pixabay/Ponciano.
Proceedings of the National Academy of Sciences: 115 (38)
Current Issue

Submit

Jump to section