## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Direct evidence of quantum transport in photosynthetic light-harvesting complexes

Edited by* Peter M. Rentzepis, University of California, Irvine, CA, and approved October 17, 2011 (received for review April 04, 2011)

## Abstract

The photosynthetic light-harvesting apparatus moves energy from absorbed photons to the reaction center with remarkable quantum efficiency. Recently, long-lived quantum coherence has been proposed to influence efficiency and robustness of photosynthetic energy transfer in light-harvesting antennae. The quantum aspect of these dynamics has generated great interest both because of the possibility for efficient long-range energy transfer and because biology is typically considered to operate entirely in the classical regime. Yet, experiments to date show only that coherence persists long enough that it can influence dynamics, but they have not directly shown that coherence does influence energy transfer. Here, we provide experimental evidence that interaction between the bacteriochlorophyll chromophores and the protein environment surrounding them not only prolongs quantum coherence, but also spawns reversible, oscillatory energy transfer among excited states. Using two-dimensional electronic spectroscopy, we observe oscillatory excited-state populations demonstrating that quantum transport of energy occurs in biological systems. The observed population oscillation suggests that these light-harvesting antennae trade energy reversibly between the protein and the chromophores. Resolving design principles evident in this biological antenna could provide inspiration for new solar energy applications.

Photosynthetic organisms employ light-harvesting antennae to capture and transport solar energy to the reaction center where charge separation occurs. This energy transport process proceeds through a complex network of coupled chromophores embedded in protein matrices of light-harvesting antenna complexes. Because of static Coulombic dipole couplings, the excitation typically delocalizes among two or more chromophores—these delocalized excited states are known as “excitons.” Though the excitonic states only delocalize across the chromophores, the protein bath is necessary for enabling energy transport by allowing the system to dissipate energy.

The precise mechanism of dissipation and whether the protein helps to steer the transport remain interesting and open questions regarding optimal design of energy transport in disordered systems. In most electronic systems, coherences among states dephase far faster than the states themselves can relax thereby precluding contributions of coherence to relaxation processes. We define transport in such systems as “classical.” Microscopically, classical transport arises from small, independent fluctuations within the protein that enable relaxation of excitonic populations through resonance energy transfer (1⇓⇓–4). This incoherent mechanism gives rise to exponential relaxation dynamics and ignores coherent dynamics.

Recent studies on photosynthetic complexes reveal that quantum coherence persists on the same timescale as population transfer—long enough to impact transport dynamics (5⇓⇓⇓–9). This experimental data implies that quantum coherences should not be ignored in photosynthetic antenna systems. Long-lived quantum coherence alone, however, is not sufficient to produce high quantum efficiency. For quantum coherences to affect efficiency or biological function, the coherences must drive dynamics and change the probability of finding the excitation in a given state. That is, coherences must couple to populations. We denote this regime as “quantum transport.”

In this study, we revisit two-dimensional electronic spectra of photosynthetic systems published previously to seek evidence of coupling between quantum coherences and excited-state populations that provide direct signatures of quantum transport (9). These data show long-lived quantum coherence in the Fenna–Matthews–Olson (FMO) complex of the green sulfur bacteria, *Chlorobaculum tepidum*. This phenomenon has been the subject of many studies, both theoretical and experimental (10⇓–12). Such long-lived coherence could potentially lead to quantum transport.

Recent theoretical work has identified oscillatory probability distributions as the signature of quantum transport (13). Conclusively identifying the origin of an oscillatory signal as arising from either population or coherence dynamics, however, presents an experimental challenge. To help discern population oscillation signals from those of quantum coherence, we start by examining transport dynamics using Redfield master equation. We show that the population oscillations appear in our data and they can be linked to quantum transport through a unique phase relationship. We conclude by discussing the implications of quantum transport both for our understanding of the microscopic protein environment and our phenomenological understanding of energy transfer.

To derive the origin of the oscillating probabilities, we represent the excitonic states of the FMO complex statistically as a quantum ensemble using the reduced density operator formalism (14). We expand the density operator, *ρ*, in the excitonic basis, , which are time-independent eigenstates constructed from linear combinations of the coupled excited states of the bacteriochlorophyll *a* chromophores (the sites). This excitonic basis corresponds to the eigenbasis of the time-averaged system Hamiltonian, and we will restrict our discussion throughout this paper to this basis. In this representation, the reduced density matrix is [1]The diagonal elements of the density matrix, , represent ensemble populations. The off-diagonal elements, *ρ*_{ij} for *i* ≠ *j*, represent coherences and are complex numbers carrying phase information. Experimentally, we probe the ensemble average of the populations or coherences over bath realizations.

Quantum dynamics within photosynthetic complexes are typically modeled with second-order perturbation theory using Lindblad or Redfield master equations. Although each method has different strengths and weaknesses (and both reproduce the dynamics discussed here), we choose to use the Redfield master equation for this work because it conveniently separates unitary time evolution from dissipative dynamics and clearly illustrates the origin of quantum transport due to coupling between populations and coherences. The Redfield master equation can be expressed as [2]where is the system Hamiltonian, and is the relaxation superoperator. In essence, all the dynamics important for function are contained within the relaxation superoperator. Treating the interaction perturbatively, an expansion of the Redfield equation for the system of excitonic states yields [3]The Redfield equation shows that temporal evolution of a particular element of the density matrix depends on the energy of the system eigenstates, *ε*_{n}, as well as on relaxation within the system. The relaxation dynamics are captured within the relaxation superoperator, ; elements *κ*_{ij,kl} determine how element *ρ*_{kl} will contribute to the time evolution of the element *ρ*_{ij} of the density matrix. The physical origin of the relaxation superoperator arises from fluctuations from the time-averaged Hamiltonian that instantaneously affect the eigenstates by effectively mixing them. The relaxation superoperator elements can be divided into three different classes: population transfer terms, *ρ*_{kk} → *ρ*_{ii} (*κ*_{ii,kk}), representing classical transport; coherence transfer terms, *ρ*_{kl} → *ρ*_{ij} for *k* ≠ *l* and *i* ≠ *j* (*κ*_{ij,kl}); and coupling between populations and coherences, *ρ*_{kl} → *ρ*_{ii} for *k* ≠ *l* (*κ*_{ii,kl}). The *κ*_{ii,kl} terms represent quantum transport by our definition.

Simply interpreted, the Redfield equation represents a first-order differential rate law that describes the change in each element in the density matrix. Two different processes control the rate of change of *ρ*_{ij}: unitary dynamics, , and nonunitary contributions (relaxation) due to coupling to other elements of the density matrix, . Unitary dynamics is responsible for the sinusoidal phase oscillation of coherences. As a result, we observe quantum coherence as an oscillation in peak amplitude, or quantum beat, with a fixed frequency proportional to the energy gap. For a population term, however, the rate of change depends only on relaxation dynamics because the unitary oscillation vanishes when the energy gap *ε*_{i} - *ε*_{j} is zero. For a population term, the summation notation expands to [4]

From Eq. **4**, we can see that population dynamics is governed by two distinct mechanisms: The first term represents transfer to and from other populations via resonance energy transport, and the second term causes oscillations driven by coherences. Population transfer via incoherent, resonance energy transport leads to exponential growth or decay depending on the direction of transport. The final term in Eq. **4** shows that populations will take on the oscillatory character of coherences.

Fig. 1 schematically depicts the three types of transfer possible among elements of the density matrix. In classical resonance energy transfer (blue), the dynamics is mostly confined to the electronic states of the chromophores. That is, energy is transferred monotonically from one state to the next, and the population dynamics is entirely decoupled from the coherence dynamics. Next, coherence transfer (green) can accompany the population transfer allowing relaxation from one superposition state to another. In quantum transport (red), coherences and populations couple, allowing coherences to drive transfer between populations. This process yields oscillatory population dynamics.

## Results

To search for evidence of quantum transport in photosynthetic systems, we revisit two-dimensional electronic spectroscopic data on the FMO complex published previously (9). (A complete discussion of methods and materials can be found in ref. 9.) Phase-sensitive two-dimensional electronic spectroscopy permits observation and characterization of population and coherence dynamics (15). The experimental and theoretical details for this spectroscopy have been described elsewhere (16⇓⇓–19). In essence, the data from 2D spectroscopy can be interpreted as a correlation map (*ω*_{τ},*ω*_{t}) of the “input” energy, *ℏω*_{τ}, and the “output” energy, *ℏω*_{t}. A series of 2D spectra as a function of waiting time reveals time-dependent evolution of a system after excitation. In rephasing pathways, the photon echo eliminates inhomogeneous broadening in the antidiagonal direction, which further improves spectral resolution. The diagonal peak (*ω*_{τ} = *ω*_{t}) contains information on the population states and their relaxation dynamics. The existence of a cross-peak (*ω*_{τ} ≠ *ω*_{t}) indicates a coupling and transfer pathway between population states that cannot be observed directly using linear absorption spectroscopy. The spectral positions of the cross-peaks help to identify couplings and coherences between excitonic states. Theoretical studies predict quantum beating due to unitary evolution of coherences will appear in cross-peaks in rephasing pathways but not in diagonal peaks (20).

We start our investigation by first isolating the coherence beating signal due to unitary evolution using the same procedure discussed in detail in prior works (9, 21). This coherence beating signal is used only for comparison to determine whether coherences couple to populations. A representative absorptive rephasing 2D spectrum of the FMO complex is shown in Fig. 2 (*Inset*). We consider only the absorptive (real) portion of the rephasing signal, which eliminates contributions from interfering response pathways and improves spectral resolution. On the diagonal, excitonic features blend together, but the lowest energy exciton 1 peak is most well resolved. Downhill energy transfer between two coupled excitonic states gives rise to cross-peaks below the diagonal. Because of peak broadening and the small energy difference between excitons 2 and 3, the cross-peaks between excitons 1 and 2 and between excitons 1 and 3 overlap with each other, apparently forming one well-resolved cross-peak below the main diagonal.

We observe two pronounced quantum beating frequencies within the amplitude oscillations of the most well-resolved cross-peak (green dot) (21). The coherence beating signal (green trace) in Fig. 2 was obtained by subtracting two exponentials corresponding to population transfer dynamics. A Fourier transform of this signal reveals approximately twice as great a contribution from the lower frequency (*ω*_{12}) component as from the high frequency component (*ω*_{13}).

Next, we isolate the population signal. In rephasing 2D spectra, the perturbative pathways used to probe the populations all involve interactions with light of the same color. Therefore, signals associated with populations will appear on the main diagonal (*ω*_{τ} = *ω*_{t}). We hypothesize that coupling between populations and coherences will be strongest for coherence and populations involving a common state. Therefore, we expect the coherence signal isolated above to drive population oscillations on the main diagonal features either directly above or directly to the left of the observed coherence signal. The width of our laser spectrum precludes precise measurements of the lowest energy exciton population. We therefore extrapolate vertically from the cross-peak signal to the main diagonal to extract a population signal. As we will show later, the observed signal appears throughout the lower part of the diagonal feature and is not strongly dependent on this choice.

To isolate the oscillatory component of the population signal from the exponential relaxation dynamics, we again remove a biexponential decay. The data prior to *T* = 80 fs have been excluded from our analysis to avoid pulse overlap effects. The resulting oscillatory population signal (red) is plotted in Fig. 2 above the quantum coherence beating signal (green) described previously. We observe pronounced oscillation in the population with similar frequencies to those in the coherence signal but phase-shifted from the coherence signal. As a control, we have checked that this result does not change qualitatively if the initial coherence time points are discarded providing discrimination against pulse overlap in the indirect, *τ*, domain as well.

We examine this phase shift further by quantifying the frequency and phase of the population signal relative to the quantum coherence beating. Using the procedure described by Hayes et al. (21), we first fit the cross-peak quantum beating signal with a linear combination of two exponentially damped sine functions, [5]where *A*_{12(13)}, *ω*_{12(13)}, *ϕ*_{12(13)}, and *γ*_{12(13)} represent the magnitude, frequency, phase, and lifetime associated with the coherence signal from excitons 1 and 2 (1 and 3) (all regression parameters are given in Supporting Information). Next, we model the functional form of the population oscillation signal. For this regression, we hold *A*_{12(13)}, *ω*_{12(13)}, *ϕ*_{12(13)}, and *γ*_{12(13)} constant based on the above fit to the coherence signal and fit only three parameters: a phase shift, ΔΦ, and two weighted contributions, *κ*_{22,12} and *κ*_{22,13}, representing the relaxation superoperator elements that link coherences to population, [6]A comparison of the fits (solid) and raw data (dashed) is shown in Fig. 2. The resultant phase shift is approximately 90° (ΔΦ = 85 ± 9°). The 90° phase shift of the population oscillation signal allows us to exclude the possibility that this signal is quantum beating due to unitary evolution of the coherences. This phase shift also effectively discriminates against signals arising from the wings of the coherence feature; such a signal would be in phase with the coherence signal. Similarly, coherence beating due to an excited-state absorption pathway would be 180° out of phase.

This 90° phase shift is unusual in two-dimensional spectroscopy. We designed a filter to better visualize the extent of this signal. We employ a z-transform filter because it maintains phase information while providing better filtering than a simple Fourier transform (22). Applying a z-transform filter across the waiting time, *T*, axis of the dataset, we isolate long-lived beating signals (dephasing rate, Γ, less than 30 cm^{-1}) with a beat frequency, *ω*_{T}, between 155 and 163 cm^{-1}, corresponding to the difference frequency between excitons 1 and 2 (23). This beating signal appears in our data only at the position of the expected exciton 1-2 cross-peak and on the diagonal directly above it (Fig. 3*A*). The phase of the signal (Fig. 3*B*) shows the characteristic phase shift of approximately 90° throughout the spectral feature. A similar relationship exists for the 200 cm^{-1} beat frequency corresponding to the difference frequency between excitons 1 and 3 (see Supporting Information). This filter confirms that the signal is specific to the exciton and population in question and does not arise from vibration, laser fluctuations, or other noise sources.

## Discussion

Oscillations occur in both the coherence beating signal and in the population signal shown in Fig. 2. The existence and position of the coherence beating signal agree well with a theoretical study by Cheng and Fleming (20). However, electronic quantum coherence cannot explain beating on the diagonal in the rephasing pathway. Quantum transport simultaneously explains the oscillation observed on the diagonal signal and the observed phase relationship.

To understand this 90° phase shift, we again examine the relationship between coherences and populations created by quantum transport terms derived in Eq. **4**. The quantum transport term demonstrates that the time derivative of population will include an oscillatory component due to coupling to coherences. Just as the integral of cosine is sine, the Redfield equation predicts that population oscillations should have a 90° phase shift from the coherence terms that drive them. No other mechanism has been suggested to create a signal phase-shifted by 90°. We therefore interpret the observed population oscillation as evidence of quantum transport and verify coupling between populations and coherences.

Beyond the oscillations observed here, quantum transport has many ramifications that help explain other observations in photosynthetic complexes. For example, although the current analysis focuses on how coherences drive population oscillations, it is interesting to consider the opposite process through which populations couple to coherences allowing them to borrow lifetime from populations. This process directly explains observations of long-lived quantum coherence. In a more microscopic analogy, the population oscillations can be thought of as Rabi-like oscillations that indicate strong mixing between the excitonic system and the bath (as compared to a true Rabi oscillation arising from strong coupling to an optical field). Following this analogy to Rabi oscillations, we can speculate that the coupling between the system and bath is strong enough that fluctuations in the protein bath effectively entangle bath modes with the system in a nontrivial way. This entangling allows the excitonic states of the FMO to trade energy with the protein bath on a femtosecond timescale.

The observed coupled dynamics of populations and coherences demonstrates a unique and more active role of the protein matrix in energy transport beyond the coherence protection mechanisms proposed by Lee et al., which invoked weakly coupled correlated baths (6). Furthermore, an oscillation in probability distribution in the nondegenerate excitonic states suggests that the total energy of the “system” may not be conserved. An oscillation in total energy of the chromophores does not impugn energy conservation, but again suggests that the system is strongly coupled to the bath such that a varying amount of the excitation is delocalized onto the bath modes as a function of time. We hypothesize that some energy is reversibly transferred between the chromophores and these bath modes throughout the transport process. To fully capture this time-dependent mixing of system and bath, we must consider the protein bath (or some portion of it) as an integral part of the energy transport system.

In summary, the oscillatory population dynamics that we observe suggests protein-assisted coupling between populations and coherences. As population feeds into coherences, coherence lifetime increases to match nonequilibrium population lifetime. In addition, transfer of coherence into population drives the population dynamics and imparts an oscillatory behavior to the energy transfer process. This observation provides direct evidence of quantum transport through the FMO complex. Because the excited states of the FMO complex are nondegenerate in energy, oscillating populations imply oscillations of the total excitonic energy. This dynamics requires that the system become mixed with the protein bath, leading to reversible exchange of energy between the two. Such a time-dependent mixing of system and bath demands a rethinking of the nature of electronic excitations and their dynamics in complex solvation environments.

## Materials and Methods

### Sample Preparation.

Trimers of the FMO from *Chlorobaculum tepidum* (24) in a buffer of 800 mM tris/HCl (pH 8.0) were mixed 35∶65 vol/vol in glycerol with 0.1% lauryldimethylamine 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.

### 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 elsewhere (18, 25).

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 (18).

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.

## Acknowledgments

We acknowledge the experimental efforts of Kelly Fransted and Dugan Hayes in collecting this data. We also acknowledge Jianzhong Wen and Robert E. Blankenship for providing the FMO sample. G.S.E. and G.P. acknowledge support from the Air Force Office of Scientific Research (FA9550-09-1-0117), the Defense Advanced Research Projects Agency (Grants N66001-10-1-4060 and N66001-10-1-4022), and the Searle Foundation for support of this work. S.M. gratefully acknowledges the support of National Science Foundation Grant CHE-1058791, Defense Advanced Research Projects Agency Grant BAA-10-40 QUBE, and the Chemical Sciences, Geosciences and Biosciences Division, Office of Basic Energy Sciences, Office of Science, US Department of Energy.

## Footnotes

↵

^{1}G.P and D.V.V. contributed equally to this work.- ↵
^{2}To whom correspondence should be addressed. E-mail: gsengel{at}uchicago.edu.

Author contributions: G.P., D.V.V., D.A., S.M., and G.S.E. designed research; G.P. D.V.V., D.A., S.M., and G.S.E. performed research; G.P., D.V.V., D.A., J.R.C., N.H.C.L., S.M., and G.S.E. analyzed data; and G.P., D.V.V., D.A., S.M., and G.S.E. wrote the paper.

The authors declare no conflict of interest.

*This Direct Submission article had a prearranged editor.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1105234108/-/DCSupplemental.

## References

- ↵
- ↵
- Redfield AG

- ↵
- ↵
- ↵
- ↵
- Lee H,
- Cheng YC,
- Fleming GR

- ↵
- ↵
- ↵
- Panitchayangkoon G,
- et al.

- ↵
- ↵
- ↵
- ↵
- ↵
- Mukamel S

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵

## Citation Manager Formats

## Sign up for Article Alerts

## Jump to section

## You May Also be Interested in

### More Articles of This Classification

### Physical Sciences

### Chemistry

### Related Content

- No related articles found.

### Cited by...

- Elucidation of near-resonance vibronic coherence lifetimes by nonadiabatic electronic-vibrational state character mixing
- Light harvesting in phototrophic bacteria: structure and function
- Engineering Coherence Among Excited States in Synthetic Heterodimer Systems
- Photosynthetic reaction center as a quantum heat engine
- Electronic resonance with anticorrelated pigment vibrations drives photosynthetic energy transfer outside the adiabatic framework