Proceedings of the National Academy of Scienceswww.pnas.org
- > Current Issue
- > vol. 107 no. 29
- >
Long-lived quantum coherence in photosynthetic complexes at physiological temperature
- Gitt Panitchayangkoona,
- Dugan Hayesa,
- Kelly A. Fransteda,
- Justin R. Carama,
- Elad Harela,
- Jianzhong Wenb,
- Robert E. Blankenshipb, and
- Gregory S. Engela,1
- Communicated by Graham R. Fleming, University of California, Berkeley, CA, May 3, 2010 (received for review January 26, 2010)
Abstract
Photosynthetic antenna complexes capture and concentrate solar radiation by transferring the excitation to the reaction center that stores energy from the photon in chemical bonds. This process occurs with near-perfect quantum efficiency. Recent experiments at cryogenic temperatures have revealed that coherent energy transfer—a wave-like transfer mechanism—occurs in many photosynthetic pigment-protein complexes. Using the Fenna–Matthews–Olson antenna complex (FMO) as a model system, theoretical studies incorporating both incoherent and coherent transfer as well as thermal dephasing predict that environmentally assisted quantum transfer efficiency peaks near physiological temperature; these studies also show that this mechanism simultaneously improves the robustness of the energy transfer process. This theory requires long-lived quantum coherence at room temperature, which never has been observed in FMO. Here we present evidence that quantum coherence survives in FMO at physiological temperature for at least 300 fs, long enough to impact biological energy transport. These data prove that the wave-like energy transfer process discovered at 77 K is directly relevant to biological function. Microscopically, we attribute this long coherence lifetime to correlated motions within the protein matrix encapsulating the chromophores, and we find that the degree of protection afforded by the protein appears constant between 77 K and 277 K. The protein shapes the energy landscape and mediates an efficient energy transfer despite thermal fluctuations.
Energy transfer through photosynthetic pigment-protein complexes operates with exceptionally high quantum efficiency (1). Recent studies have demonstrated that energy moves through antennae using not only a classical hopping mechanism but also a manifestly quantum mechanical wave-like mechanism at cryogenic temperatures (2–5). Theoretical studies of this process within the Fenna–Matthews–Olson antenna complex (FMO) show that this quantum transport mechanism requires a balance between unitary (oscillatory) and dissipative (dephasing) dynamics; further, this balance appears to be optimized near room temperature and contributes to the robustness of the process (6–9). This theory demands that quantum coherence persist long enough to affect transport, but quantum beating has never been observed in FMO at physiological temperature.
The FMO pigment-protein complex from Chlorobium tepidum serves as a model system for photosynthetic energy transfer processes (2, 10–13). This complex conducts energy from the larger light-harvesting chlorosome to the reaction center in green sulfur bacteria (14, 15). Each noninteracting FMO monomer contains seven coupled bacteriochlorophyll-a chromophores arranged asymmetrically, yielding seven nondegenerate, delocalized molecular excited states called excitons (11, 16). The small number of distinct states makes this particular complex attractive for theoretical studies of transport dynamics. As shown by Ishizaki and Fleming (13), the arrangement of the chromophores in FMO results in a downhill, rugged energetic landscape with two distinct routes through which an excitation can travel to reach the lowest energy state. While classical trajectories can navigate such funnel-like landscapes, the wave-like motion through the complex improves efficiency by avoiding kinetic traps. In higher plants, this mechanism likely becomes more important because the landscape is more rugged without a downhill arrangement (17).
Recent investigations of photosynthetic systems at 77 K have found evidence of coherent energy transfer in many antenna complexes and even in the reaction center of purple bacteria (2–4). This wave-like energy transfer mechanism, however, can contribute to the near-perfect quantum efficiency of photosynthesis only if coherences survive in these systems during energy transfer at physiological temperatures. As temperature increases, thermally excited vibrational modes of the protein bath drive larger energetic fluctuations, thereby accelerating decoherence (14, 18). Although this dephasing seems unfavorable, Mohseni et al. (6) and Plenio and Huelga (7) have independently shown that the delicate interplay between quantum coherence and dephasing can create fast and unidirectional transfer pathways within these complexes, resulting in highly efficient electronic energy transfer (8, 9, 19). This scheme exploits quantum coherence to overcome an energy barrier, but subsequent dephasing processes trap the excitation at the target site. Optimal transport therefore requires both dephasing and coherent energy transfer.
The initial excitation or transfer event necessarily creates quantum coherence because both the dipole and site operators do not commute with the system Hamiltonian. For a system of two excitons described by Ψ(t) = c1ϕ1 + c2ϕ2, the time evolution of the density matrix is given by 
[1]The first two terms represent populations in the excitonic basis, whereas the latter two describe coherences. The phase factors in the coherence terms are responsible for quantum beating, which appears as a periodic modulation of population in the site basis and peak amplitude. The frequency of this beating corresponds to the energy difference between the two excitons giving rise to that particular quantum coherence. Traditionally this phenomenon is ignored in transport dynamics because fast electronic dephasing generally destroys quantum coherence before it can impact the transport process. For example, at cryogenic temperature, coherences between ground and excited states in FMO dephase in approximately 70 fs. In contrast, coherences among excited states have been shown to persist beyond 660 fs—long enough to improve transport efficiency (2). Such a coherence can persist only if the electronic spectral motion among chromophores is strongly correlated, as demonstrated in a conjugated polymer system by Collini and Scholes (20). Biologically, these correlations arise because the protein environment forces transition energies of chromophores to fluctuate together due to spatial uniformity of the dielectric bath. The resultant long-lived quantum coherences evolve in time, creating periodic oscillations in both spectral signals and wave-packet position. This quantum beat provides the signature of quantum coherence.
[1]The first two terms represent populations in the excitonic basis, whereas the latter two describe coherences. The phase factors in the coherence terms are responsible for quantum beating, which appears as a periodic modulation of population in the site basis and peak amplitude. The frequency of this beating corresponds to the energy difference between the two excitons giving rise to that particular quantum coherence. Traditionally this phenomenon is ignored in transport dynamics because fast electronic dephasing generally destroys quantum coherence before it can impact the transport process. For example, at cryogenic temperature, coherences between ground and excited states in FMO dephase in approximately 70 fs. In contrast, coherences among excited states have been shown to persist beyond 660 fs—long enough to improve transport efficiency (2). Such a coherence can persist only if the electronic spectral motion among chromophores is strongly correlated, as demonstrated in a conjugated polymer system by Collini and Scholes (20). Biologically, these correlations arise because the protein environment forces transition energies of chromophores to fluctuate together due to spatial uniformity of the dielectric bath. The resultant long-lived quantum coherences evolve in time, creating periodic oscillations in both spectral signals and wave-packet position. This quantum beat provides the signature of quantum coherence.
Experimentally, we probe quantum coherences in FMO using two-dimensional Fourier transform electronic spectroscopy to directly observe electronic couplings and quantum coherences as a function of time (21–24). The experimental method and theory have been described in detail elsewhere (23). In short, three laser pulses interact in the weak field limit with the sample to produce a third-order polarization (25). The first pulse creates a superposition of ground and excited states, and the phase of this coherence evolves for a time τ (coherence time) before a second pulse creates either a population or a superposition of excited states for a time T (waiting time). During this time, both population and coherence transfer occur, and coherences evolve phase. A third pulse then generates a second ground-excited state coherence. Finally, the signal pulse appears after a time t (rephasing time) and is heterodyne-detected with a local oscillator pulse in a unique phase-matched direction. A Fourier transform of the signal along the coherence and rephasing time dimensions at a fixed waiting time gives a 2D spectrum that correlates the “input” coherence frequency to an “output” rephasing frequency. Cross-peaks (off-diagonal peaks) appear in this spectrum due to coupling between states, providing evidence of energy transfer.
Two types of response pathways—rephasing and nonrephasing—contribute to the overall signal in 2D spectroscopy. We probe these two pathways individually by switching the order of the first two pulses, which controls the direction in which the phase of the coherence evolves during the coherence time. In the rephasing (nonrephasing) pathway, phase evolution proceeds in the opposite (same) direction during the coherence and rephasing times, resulting in a photon echo (free induction decay) signal. When analyzed independently, the rephasing and nonrephasing signals provide complementary information on coherences because quantum beats during the waiting time appear in different positions in the spectra (24, 26). For rephasing pathways, beating appears in the cross-peaks, whereas for nonrephasing pathways it appears in the peaks along the diagonal. The sinusoidal beating pattern observed in a cross-peak consists of a single frequency corresponding to the transition energy gap between the two excited states that give rise to this peak. The quantum beating observed in a diagonal peak, however, is the sum of several different frequencies arising from coherences with all other states, which complicates attempts to measure individual quantum coherence lifetimes. We therefore choose to focus on the off-diagonal peaks to enable more accurate measurement of the lifetime of an individual quantum coherence.
Results
Representative 2D spectra in Fig. 1 A–D show the real part of the third-order nonlinear response taken at four different temperatures (77 K, 125 K, 150 K, and 277 K) at a waiting time T = 400 fs. The lowest energy peak is well resolved at low temperatures, but more vibrational modes of the protein are occupied at 277 K, resulting in faster dephasing between the ground and excited states. According to the time-bandwidth product, spectral resolution is inversely related to signal lifetime, so the short duration of the signal along τ and t results in loss of spectral resolution in the respective frequency dimensions. Furthermore, the spectrum acquired from the flowing sample also suffers from rearrangement of the hydration shell around the protein within the first 70 fs of waiting time, causing a rapid growth in the antidiagonal linewidths. In contrast, the low-temperature spectra maintain resolution and narrow peak shapes even after 1 ps because the translational motions of the solvent are frozen.
Fig. 1.
Temperature dependence data. Representative two-dimensional electronic spectra of FMO are shown at the waiting time T = 400 fs and 77 K (A), 125 K (B), 150 K (C), and 277 K (D). The data are shown with an arcsinh color scale to highlight small features in both negative and positive portions of the real third-order nonlinear response. Peaks broaden at higher temperature due to faster dephasing between ground and excited states, preventing resolution of the lowest excited state. The quantum beat signals are extracted at the spectral position (white circle) corresponding to the location of 1–3 cross-peak predicted by a theoretical study (27). The beating signals (E) demonstrate agreement in phase and beating frequency among all four temperatures while showing shorter quantum beat lifetimes at higher temperatures.
Fig. 1E shows the absolute value of the amplitude of the highlighted cross-peak as a function of waiting time for each temperature. This feature was chosen because it is well separated from the congested diagonal peaks and the amplitude therefore beats with only a single frequency. The spectral coordinates of the cross-peak were chosen according to the energies of excitons 1 and 3 taken from the Hamiltonian calculated by Adolphs and Renger (27). According to the projection-slice theorem, we can separate the third-order response into real (absorptive) and imaginary (dispersive) portions by fitting the phase of the signal to the pump-probe data taken at the same waiting time as described by Brixner et al. (23). Separating the data into real and imaginary parts will improve spectral resolution by eliminating wide dispersive features. However, scatter from the frozen samples creates an interferometric signal that cannot be separated from the pump-probe signal; this interference introduces errors in the real portion of the phased 2D datasets that could be mistaken for quantum beating. Therefore, we present the peak amplitude in absolute value to eliminate possible phase errors. Examining the cross-peak amplitudes at all four temperatures, we see clear evidence of quantum beating. The beating signals demonstrate excellent agreement in both the phase and frequency across all temperatures, indicating that the same phenomenon discovered at 77 K extends to at least 277 K. At 277 K, the coherence envelope fits well to a single exponential decay. The observed dephasing rate shows strong temperature dependence with coherence observable only to about 300 fs at 277 K, corresponding to a dephasing rate approximately four times faster than the rate at 77 K.
Discussion
Extracting the beating signal from the 277-K spectrum presents an additional challenge because of the broad lineshapes. Even using only the absorptive portion of the response to improve resolution, we cannot resolve individual cross-peaks in the spectra as illustrated in Fig. 2A. With such broad lineshapes, beating from neighboring peaks confounds attempts to isolate a single quantum beating signal. We therefore isolate rephasing and nonrephasing signals separately both on and off the main diagonal to demonstrate that the beating occurs at the cross-peak only in the rephasing portion of the signal as predicted by theory (24).
Fig. 2.
Room temperature quantum beating. Integrated peak amplitudes taken from the absorptive part of the 2D signal at 277 K are plotted as a function of waiting time. The rephasing signal at the indicated cross-peak (red line) shows multiple periods of quantum beating as it decays, while the nonrephasing signal (green line) shows no beating. As a further comparison, the rephasing signal at the indicated diagonal position (blue line) shows only smooth population decay. The magnitude of the diagonal signal was scaled by a factor of 0.2 to facilitate comparison with the off-diagonal signal. The inset shows the 2D spectrum of the absorptive rephasing signal taken at T = 140 fs. Below the spectra, Feynman diagrams representing the relevant Liouville pathways illustrate why beating arises only in the rephasing cross-peak signal.
The Feynman diagrams in Fig. 2 show the Liouville-space pathways giving rise to the cross-peak signal in both rephasing and nonrephasing data indicated by the red and green arrows, respectively. While both rephasing (red) and nonrephasing (green) pathways produce a signal at the same spectral position, these two pathways yield significantly different information during the waiting time T. In the rephasing pathway, we observe quantum beating on the cross-peak due to the phase factor exp(-i(E3 - E1)t/ℏ) in Green’s function. The expected frequency of the beating is 
. In contrast, Green’s function of the cross-peak during waiting time in the nonrephasing pathway equals 1, which indicates that it does not beat.
. In contrast, Green’s function of the cross-peak during waiting time in the nonrephasing pathway equals 1, which indicates that it does not beat.
To verify that we have extracted actual quantum beating, we first contrast the rephasing signal on the main diagonal to the off-diagonal signal. The amplitude of the rephasing diagonal peak (blue line) shows a smooth decay over the waiting time T, whereas the rephasing cross-peak amplitude (red line) shows multiple periods of quantum beating as it decays. Second, we compare the beating in the rephasing signal to the nonrephasing signal at the cross-peak (green line). The nonrephasing signal in the region of the cross-peak demonstrates similar population dynamics to the rephasing signal but shows no quantum beating. These data confirm that the beating observed at room temperature arises from an electronic coherence, despite the loss of spectral resolution.
After elimination of population transfer dynamics from the beating signal by subtraction of a multiexponential fit, we extract the beating signals from the coherence between excitons 1 and 3 as shown in Fig. 3. Each signal is fit with the product of a single sine function and a single exponential decay from which the dephasing rate can be estimated. All regressions start from 80 fs to avoid pulse overlap effects and contributions from solvent reorganization at room temperature. The beating frequency found for each fit was approximately 160 cm-1; in comparison, we would predict the frequency to be 198 cm-1, corresponding to the energy difference between excitons 1 and 3 given in the Hamiltonian referenced earlier. The dephasing rate taken from the exponential part of the fitting function is plotted as a function of temperature in the inset of Fig. 3. We observe a linear relationship between temperature and dephasing rate; the slope obtained from the linear regression is 0.52 ± 0.07 cm-1/K (SD). The data agree with prior work by Aspuru-Guzik and co-workers (8) who use an uncorrelated local ohmic spectral density with a frequency cutoff to model the FMO protein. This work predicts the dephasing rate as a function of temperature to be 
, where ER is a reorganization energy (35 cm-1) and ωc is the cutoff frequency (150 cm-1). Such an uncorrelated bath model can accurately represent the observed excited state dynamics. However, correlated bath motions are required to simultaneously fit dephasing dynamics of both zero-quantum coherences and single-quantum coherences. This correlation is the mechanism by which the protein enables long-lived quantum coherence and coherent energy transfer.
, where ER is a reorganization energy (35 cm-1) and ωc is the cutoff frequency (150 cm-1). Such an uncorrelated bath model can accurately represent the observed excited state dynamics. However, correlated bath motions are required to simultaneously fit dephasing dynamics of both zero-quantum coherences and single-quantum coherences. This correlation is the mechanism by which the protein enables long-lived quantum coherence and coherent energy transfer. 
Fig. 3.
Temperature dependence of coherence dephasing. Integrated cross-peak amplitudes are taken from the absolute value of the combined 2D signal (rephasing + nonrephasing) after removal of exponential population decay at 77 K, 125 K, and 150 K (colored solid lines). The amplitude at 277 K is taken from the absorptive portion of the rephasing signal (dashed red line). The beating signals are normalized to their respective maxima and fit to the product of a sine function and an exponential decay (solid black lines). The beating frequency is given for each temperature. The dephasing rate taken from the exponential part of the fit is plotted as a function of temperature along with standard errors in the inset. The statistically weighted linear fit of these points (dashed black line) has a slope of 0.52 ± 0.07 cm-1/K (SD).
Furthermore, Aspuru-Guzik and co-workers (8) found that at room temperature the dephasing frequency is on the order of the excitonic energy gaps and couplings. This temperature corresponds to optimal transfer efficiency because dephasing traps an excitation at an energetic minimum, but does not collapse a coherence before it has a chance to complete at least a period of quantum beating and overcome the initial energy barrier. Intuitively, for a complete sampling of all possible pathways, all quantum coherences should have a chance to complete at least one oscillatory cycle. While coherence prevents loss to local minima, fast dephasing prevents loss to exciton recombination and nonradiative decay. Thus, optimal efficiency is attained when dephasing completely destroys coherences shortly after a full period of quantum beating. Our experimental results show that not only does quantum coherence persist at physiological temperature long enough to impact overall transport dynamics, but it lasts for nearly two full cycles, corresponding to the predictions for optimal transport efficiency (8). Therefore, these data support the hypothesis that environmentally assisted quantum transport in FMO is relevant for biological function.
Protection at cryogenic temperature by correlated spectral motion within a protein was demonstrated in 2007 by Lee et al. (3). In theory, the origin of this correlation could be as simple as a local site fluctuation in a site that contributes to two excitons. Despite their spatial proximity, the correlation expected between this pair of excitons, however, is only 2.6% because no single site contributes significantly to both excitons 1 and 3 based on the Hamiltonian from Adolphs and Renger (27). Yet, this coherence shows strong, measurable beating in the data above. We therefore speculate that this phenomenon does not derive from a simple entanglement in the site basis but rather from the complex bath environment. Further, the time scale of these fluctuations must be fast relative to our measurement time or we would observe oscillation of the locations of spectral features. No change in spectral location is observed. Motions of large structural elements of the protein provide the necessary spatial characteristics but are far too slow. Future work will explore modeling microscopic fluctuations of the dielectric environment within the protein.
These data offer insight into the microscopic requirements for coherent transport in disordered molecular assemblies although the spectroscopy cannot directly probe local excitonic environments. The agreement in both phase and frequency of the quantum beating as a function of temperature indicates that the protection of the zero-quantum coherences in the one-exciton manifold does not require specific microscopic vibrational motions. The soft, anharmonic nature of the low-energy vibrational modes of the protein should change the vibrational environment significantly as these modes become thermally populated. If the precise nature of the modes were critical, the dephasing would increase faster than linearly, which is not observed. Thus, this effect likely owes to the polymeric nature of the protein and electronic contributions to the complex dielectric environment within the protein bath. It is reasonable to expect that synthetic systems can emulate such an environment.
Conclusion
These data prove the same quantum beating signals observed at 77 K persist to physiological temperature and show agreement in both phase and frequency, indicating that the experiment is following the same quantum coherence at all temperatures. We observe a 130-fs e-folding lifetime for this excited state coherence at 277 K and observe quantum coherence lasting beyond 300 fs, showing that evolution has had the opportunity to exploit the theorized environmentally assisted quantum transport (EnAQT) mechanism for biological function. Beating may persist even longer, but we cannot separate actual dephasing from deleterious interference from other spectral features. Thus, this measurement represents a lower bound. Microscopically, this beating survives because the energies of the excited states involved fluctuate such that the energy gap remains largely constant. We hypothesize the fluctuations of the protein dielectric provide the long-range spatial correlations on the appropriate time scale to protect quantum coherence and enable coherence transfer. As suggested by Cheng et al. (26), this wave-like transfer mechanism results in oscillatory population dynamics, allowing particular sites to have momentary populations higher than their respective equilibrium populations. When coupled to an energy trap such as the reaction center, this process can greatly increase the quantum yield of a system relative to the classical limit.
Materials and Methods
Sample Preparation.
FMO from C. tepidum (28) in a buffer of 800 mM tris/HCl (pH 8.0) was mixed 35∶65 vol/vol in glycerol with 0.1% laurydimethylamine oxide detergent and loaded into a 200-μm fused quartz cell (Starna). The optical density at 809 nm was 0.32. A cryostat (Oxford Instruments) was used to cool the sample to 77 K, 125 K, and 150 K. The room temperature sample was prepared by mixing the sample 35∶65 vol/vol in 800 mM tris/HCl buffer (pH 8.0) and pumped (Masterflex, Cole-Palmer) through a 200-μm quartz flow cell (Starna). The sample reservoir was kept in a water bath at 4 °C to cool the sample and prevent thermal degradation. The optical density of the room temperature sample at 809 nm was 0.35.
Data Acquisition.
A self mode-locking Ti:sapphire oscillator (Coherent, Micra) was used to seed a regenerative amplifier (Coherent, Legend Elite), which produced a 5.0-kHz pulse train of 38-fs pulses centered at 806 nm with a spectral bandwidth of 35 nm. The 10-Hz stability of the laser power during data acquisition ranged from 0.08% to 0.19%. Two pairs of phase-locked beams were generated by a diffractive optic (Holoeye). All beams were incident on identical optics except for one-degree fused silica wedges (Almaz Optics) on delay stages (Aerotech), which determined the time delays. Each delay was calibrated using spectral interferometry as described earlier (23, 29).
Neutral density filters with total optical density 3.1 at 809 nm attenuated the local oscillator beam. The total power incident on the sample was 4.8 nJ (1.6 nJ/pulse), which was focused to a spot size less than 70 μm. A 0.3-m spectrometer (Andor Shamrock) frequency-resolved the emitted signal and local oscillator beam, which were captured on a 1,600 × 5 pixel region of a back-illuminated, thermoelectrically cooled CCD (Andor Newton). Scatter subtraction, Fourier windowing, and transformation to frequency-frequency space were performed as described previously (23).
Two-dimensional data were collected at waiting times (T) in 20-fs increments for all temperatures. At each waiting time, the coherence time was scanned from -500 to 500 fs in steps of 4 fs. A 2D spectrum at T = 0 fs was taken every 200 fs to monitor sample integrity. No degradation was observed after 100 continuous 2D data acquisitions at and below 150 K. Pump-probe data were used to phase only data from the 277-K flowing sample; scatter from cryogenic samples prevented accurate phasing of those datasets.
Acknowledgments
This work was supported by Defense Advanced Research Projects Agency Grant HR0011-09-1-0051 and Air Force Office of Scientific Research Grant FA9550-09-1-0117, as well as funding from the Dreyfus Foundation and the Searle Foundation. Funding for E.H. was provided by National Science Foundation Grant DMR-0844115 and the Institute for Complex Adaptive Matter Branches Cost-Sharing Fund. Funding for J.W. and R.E.B. was provided by Grant DEFG02-07ER15846 from the Photosynthetic Systems program of the Basic Energy Sciences division of Department of Energy.
Footnotes
- 1To whom correspondence should be addressed at: Gordon Center for Integrative Science, 929 East 57th Street, GCIS E119, Chicago, IL 60637. E-mail: gsengel@uchicago.edu.
- Author contributions: G.P., K.A.F., and G.S.E. designed research; G.P., D.H., K.A.F., and J.R.C. performed research; J.W. and R.E.B. contributed reagents; G.P., D.H., J.R.C., E.H., and G.S.E. analyzed data; and G.P., D.H., and G.S.E. wrote the paper.
- The authors declare no conflict of interest.
References
Facebook
Twitter- Google+
CiteULike
Connotea
Delicious
Digg
Mendeley
HighWire Press-hosted articles citing this article
- This Issue
July 20, 2010
vol. 107 no. 29
Masthead (PDF)Table of Contents
PNAS July 20, 2010 vol. 107 no. 29 12766-12770
Access
- AbstractFree
- Full Text (HTML)
- Full Text (PDF)
- Figures Only
Article Tools
Share
- Facebook
- Twitter
- Google+
- CiteULike
- Connotea
- Delicious
- Digg
- Mendeley
No comments:
Post a Comment