Electronic coherence effects in photosynthetic light harvesting

Photosynthetic light harvesting is a paradigmatic example for quantum effects in biology. In this work, we review studies on quantum coherence effects in the LH2 antenna complex from purple bacteria to demonstrate how quantum mechanical rules play important roles in the speedup of excitation energy transfer, the stabilization of electronic excitations, and the robustness of light harvesting in photosynthesis. Subsequently, we present our recent theoretical studies on exciton dynamical localization and excitonic coherence generation in photosynthetic systems. We apply a variational-polaron approach to investigate decoherence of exciton states induced by dynamical fluctuations due to system-environment interactions. The results indicate that the dynamical localization of photoexcitations in photosynthetic complexes is significant and imperative for a complete understanding of coherence and excitation dynamics in photosynthesis. Moreover, we use a simple model to investigate quantum coherence effects in intercomplex excitation energy transfer in natural photosynthesis, with a focus on the likelihoods of generating excitonic coherences during the process. Our model simulations reveal that excitonic coherence between acceptor exciton states and transient nonlocal quantum correlation between distant pairs of chromophores can be generated through intercomplex energy transfer. Finally, we discuss the implications of these theoretical works and important open questions that remain to be answered.     

Population and coherence dynamics in light harvesting complex II (LH2)

The electronic excitation population and coherence dynamics in the chromophores of the photosynthetic light harvesting complex 2 (LH2) B850 ring from purple bacteria (Rhodopseudomonas acidophila) have been studied theoretically at both physiological and cryogenic temperatures. Similar to the well-studied Fenna-Matthews-Olson (FMO) protein, oscillations of the excitation population and coherence in the site basis are observed in LH2 by using a scaled hierarchical equation of motion approach. However, this oscillation time (300 fs) is much shorter compared to the FMO protein (650 fs) at cryogenic temperature. Both environment and high temperature are found to enhance the propagation speed of the exciton wave packet yet they shorten the coherence time and suppress the oscillation amplitude of coherence and the population. Our calculations show that a long-lived coherence between chromophore electronic excited states can exist in such a noisy biological environment.     

Inhomogeneous dephasing masks coherence lifetimes in ensemble measurements

An open question at the forefront of modern physical sciences is what role, if any, quantum effects may play in biological sensing and energy transport mechanisms. One area of such research concerns the possibility of coherent energy transport in photosynthetic systems. Spectroscopic evidence of long-lived quantum coherence in photosynthetic light-harvesting pigment protein complexes (PPCs), along with theoretical modeling of PPCs, has indicated that coherent energy transport might boost efficiency of energy transport in photosynthesis. Accurate assessment of coherence lifetimes is crucial for modeling the extent to which quantum effects participate in this energy transfer, because such quantum effects can only contribute to mechanisms proceeding on timescales over which the coherences persist. While spectroscopy is a useful way to measure coherence lifetimes, inhomogeneity in the transition energies across the measured ensemble may lead to underestimation of coherence lifetimes from spectroscopic experiments. Theoretical models of antenna complexes generally model a single system, and direct comparison of single system models to ensemble averaged experimental data may lead to systematic underestimation of coherence lifetimes, distorting much of the current discussion. In this study, we use simulations of the Fenna-Matthews-Olson complex to model single complexes as well as averaged ensembles to demonstrate and roughly quantify the effect of averaging over an inhomogeneous ensemble on measured coherence lifetimes. We choose to model the Fenna-Matthews-Olson complex because that system has been a focus for much of the recent discussion of quantum effects in biology, and use an early version of the well known environment-assisted quantum transport model to facilitate straightforward comparison between the current model and past work. Although ensemble inhomogeneity is known to lead to shorter lifetimes of observed oscillations (simply inhomogeneous spectral broadening in the time domain), this important fact has been left out of recent discussions of spectroscopic measurements of energy transport in photosynthesis. In general, these discussions have compared single-system theoretical models to whole-ensemble laboratory measurements without addressing the effect of inhomogeneous dephasing. Our work addresses this distinction between single system and ensemble averaged observations, and shows that the ensemble averaging inherent in many experiments leads to an underestimation of coherence lifetimes in individual systems.     

Femtosecond pump-probe spectroscopy of I2 in a dense rare gas environment: a mixed quantum/classical study of vibrational decoherence

The process of decoherence of vibrational states of I2 in a dense helium environment is studied theoretically using the mixed quantum/classical method based on the Bohmian formulation of quantum mechanics [E. Gindensperger, C. Meier, and J. A. Beswick, J. Chem. Phys. 113, 9369 (2000)]. Specifically, the revival of vibrational wave packets is a quantum phenomena which depends sensitively on the coherence between the vibrational states excited by an ultrafast laser pulse. Its detection by a pump-probe setup as a function of rare gas pressure forms a very accurate way of detecting vibrational dephasing. Vibrational revivals of I2 in high pressure rare gas environments have been observed experimentally, and the very good agreement with the simulated spectra confirms that the method can accurately describe decoherence processes of quantum systems in interaction with an environment.     

Quantitative investigations of quantum coherence for a light-harvesting protein at conditions simulating photosynthesis

Recent measurements using two-dimensional electronic spectroscopy (2D ES) have shown that the initial dynamic response of photosynthetic proteins can involve quantum coherence. We show how electronic coherence can be differentiated from vibrational coherence in 2D ES. On that basis we conclude that both electronic and vibrational coherences are observed in the phycobiliprotein light-harvesting complex PC645 from Chroomonas sp. CCMP270 at ambient temperature. These light-harvesting antenna proteins of the cryptophyte algae are suspended in the lumen, where the pH drops significantly under sustained illumination by sunlight. Here we measured 2D ES of PC645 at increasing levels of acidity to determine if the change in pH affects the quantum coherence; quantitative analysis reveals that the dynamics are insensitive to the pH change.     

Decoherence in an anharmonic oscillator coupled to a thermal environment: a semiclassical forward-backward approach

The decoherence of an anharmonic oscillator in a thermal harmonic bath is examined via a semiclassical approach. A computational strategy is presented and exploited to calculate the time dependence of the purity and the decay of individual matrix elements in the energy representation for a variety of initial states. The time dependence of the decoherence is found to depend on the temperature of the bath, the coupling strength, the initial state of the oscillator, and the choice of quantity measuring the decoherence. Recurrences in the purity and in the off-diagonal matrix elements are observed, as well as the collapse of these matrix elements to the diagonal, providing evidence for the retention of quantum coherence for time scales longer than that indicated by the purity. The results are used to analyze the utility of the Caldeira-Leggett and Redfield models of decoherence and to assess the dependence of dephasing rates on the degree of structure in phase space. In several cases we find that the dephasing dynamics can be described as an initial Zeno-effect regime, followed by a Caldeira-Leggett region, followed by recurrences.     

Multichromophoric Förster resonance energy transfer from b800 to b850 in the light harvesting complex 2: evidence for subtle energetic optimization by purple bacteria

This work provides a detailed account of the application of our multichromophoric F?rster resonance energy transfer (MC-FRET) theory (Phys. Rev. Lett. 2004, 92, 218301) for the calculation of the energy transfer rate from the B800 unit to the B850 unit in the light harvesting complex 2 (LH2) of purple bacteria. The model Hamiltonian consists of the B800 unit represented by a single bacteriochlorophyll (BChl), the B850 unit represented by its entire set of BChls, the electronic coupling between the two units, and the bath terms representing all environmental degrees of freedom. The model parameters are determined, independent of the rate calculation, from the literature data and by a fitting to an ensemble line shape. Comparing our theoretical rate and a low-temperature experimental rate, we estimate the magnitude of the BChl-Qy transition dipole to be in the range of 6.5-7.5 D, assuming that the optical dielectric constant of the medium is in the range of 1.5-2. We examine how the bias of the average excitation energy of the B800-BChl relative to that of the B850-BChl affects the energy transfer time by calculating the transfer rates based on both our MC-FRET theory and the original FRET theory, varying the value of the bias. Within our model, we find that the value of bias 260 cm-1, which we determine from the fitting to an ensemble line shape, is very close to the value at which the ratio between MC-FRET and FRET rates is a maximum. This provides evidence that the bacterial system utilizes the quantum mechanical coherence among the multiple chromophores within the B850 in a constructive way so as to achieve efficient energy transfer from B800 to B850.     

Signatures of correlated excitonic dynamics in two-dimensional spectroscopy of the Fenna-Matthew-Olson photosynthetic complex

Long-lived excitonic coherence in photosynthetic proteins has become an exciting area of research because it may provide design principles for enhancing the efficiency of energy transfer in a broad range of materials. In this publication, we provide new evidence that long-lived excitonic coherence in the Fenna-Mathew-Olson pigment-protein (FMO) complex is consistent with the assumption of cross correlation in the site basis, indicating that each site shares bath fluctuations. We analyze the structure and character of the beating crosspeak between the two lowest energy excitons in two-dimensional (2D) electronic spectra of the FMO Complex. To isolate this dynamic signature, we use the two-dimensional linear prediction Z-transform as a platform for filtering coherent beating signatures within 2D spectra. By separating signals into components in frequency and decay rate representations, we are able to improve resolution and isolate specific coherences. This strategy permits analysis of the shape, position, character, and phase of these features. Simulations of the crosspeak between excitons 1 and 2 in FMO under different regimes of cross correlation verify that statistically independent site fluctuations do not account for the elongation and persistence of the dynamic crosspeak. To reproduce the experimental results, we invoke near complete correlation in the fluctuations experienced by the sites associated with excitons 1 and 2. This model contradicts ab initio quantum mechanic∕molecular mechanics simulations that observe no correlation between the energies of individual sites. This contradiction suggests that a new physical model for long-lived coherence may be necessary. The data presented here details experimental results that must be reproduced for a physical model of quantum coherence in photosynthetic energy transfer.     

Quantum coherence in photosynthesis

Quantum coherence improves the quantum efficiency of excitonic energy transport within the Fenna-Matthews-Olson photosynthetic complex from the green sulphur bacterium, Chlorobium tepidum. Experimental evidence from third-order nonlinear spectroscopies provides clear evidence of quantum coherence among excited states persisting for picoseconds despite rapid (<100fs) dephasing of quantum coherence between ground and excited states. This protection of quantum coherence can arise from multiple mechanisms, but the net effect is the same: the energetic landscape is course-grained thereby improving efficiency by effectively smoothing the rugged energetic landscape while simultaneously eliminating trap states. The protein bath enables the unusual observed dynamics and illustrates some simple design principles that provide direction to synthetic efforts to mimic the effect. This communication provides an overview of experimental and theoretical notions for those interested in exploiting design principles of photosynthetic energy transfer in synthetic systems.     

Quantum nonlocality in the excitation energy transfer in the Fenna–Matthews–Olson complex

The Fenna–Matthews–Olson (FMO) complex—a pigment protein complex involved in photosynthesis in green sulfur bacteria—is remarkably efficient in transferring excitation energy from light harvesting antenna molecules to a reaction center. Recent experimental and theoretical studies suggest that quantum coherence and entanglement may play a role in this excitation energy transfer (EET). We examine whether bipartite quantum nonlocality, a property that expresses a stronger‐than‐entanglement form of correlation, exists between different pairs of chromophores in the FMO complex when modeling the EET by the hierarchically coupled equations of motion method. We compare the results for nonlocality with the amount of bipartite entanglement in the system. In particular, we analyze in what way these correlation properties are affected by different initial conditions. It is found that bipartite nonlocality only exists when the initial conditions are chosen in an unphysiological manner and probably is absent when considering the EET in the FMO complex in its natural habitat. It is also seen that nonlocality and entanglement behave quite differently in this system. In particular, for localized initial states, nonlocality only exists on a very short time scale and then drops to zero in an abrupt manner. As already known from previous studies, quantum entanglement between chromophore pairs, on the other hand, is oscillating and exponentially decaying and follow thereby a pattern more similar to the chromophore population dynamics. The abrupt disappearance of nonlocality in the presence of nonvanishing entanglement is a phenomenon we call nonlocality sudden death; a striking manifestation of the difference between these two types of correlations in quantum systems.     

Time-resolved coherent anti-Stokes Raman-scattering measurements of I2 in solid Kr: vibrational dephasing on the ground electronic state at 2.6-32 K

Time-resolved coherent anti-Stokes Raman-scattering (CARS) measurements are carried out for iodine (I2) in solid krypton matrices. The dependence of vibrational dephasing time on temperature and vibrational quantum number v is studied. The v dependence is approximately quadratic, while the temperature dependence of both vibrational dephasing and spectral shift, although weak, fits the exponential form characteristic of dephasing by pseudolocal phonons. The analysis of the data indicates that the frequency of the pseudolocal phonons is approximately 30 cm(-1). The longest dephasing times are observed for v = 2 being approximately 300 ps and limited by inhomogeneous broadening. An increase in the dephasing rate of v = 2 as the temperature is lowered to T = 2.6 K is taken as a clear indication of lattice-strain-induced inhomogeneity of the ensemble coherence.     

Semiclassical Liouville method for the simulation of electronic transitions: single ensemble formulation

In this paper, we describe a single ensemble implementation of the semiclassical Liouville method for simulating quantum processes using classical trajectories. In this approach, one ensemble of trajectories supports the evolution of all semiclassical density matrix elements, rather than employing a distinct ensemble for each. The ensemble evolves classically under a single reference Hamiltonian, which is chosen based on physical grounds; for electronic relaxation of an initially excited state, the initially populated upper surface Hamiltonian is the natural choice. Classical trajectories evolving on the reference potential then represent the time-dependent upper state population density and also the electronic coherence and the ground state density created by electronic transition. The error made in the classical motion of the trajectories for these latter distributions is compensated for by incorporating the difference between the correct and reference Liouville propagators into the calculation of the coefficients of the individual trajectories. This approach gives very accurate results for a number of model problems and cases describing ultrafast electronic relaxation dynamics.     

Theory and Simulation of the Environmental Effects on FMO Electronic Transitions

Long-lived quantum coherence has been experimentally observed in the Fenna-Matthews-Olson (FMO) light-harvesting complex. It is much debated which role thermal effects play and if the observed low-temperature behavior arises also at physiological temperature. To contribute to this debate we use molecular dynamics simulations to study the coupling between the protein environment and the vertical excitation energies of individual bacteriochlorophyll molecules in the FMO complex of the green sulphur bacterium Chlorobaculum tepidum. The so-called spectral densities, which account for the environmental influence on the excited state dynamics, are determined from temporal autocorrelation functions of the energy gaps between ground and first excited states of the individual pigments. Although the overall shape of the spectral density is found to be rather similar for all pigments, variations in their magnitude can be seen. Differences between the spectral densities for the pigments of the FMO monomer and FMO trimer are also presented.     

Electronic coherence dynamics in trans-polyacetylene oligomers

Electronic coherence dynamics in trans-polyacetylene oligomers are considered by explicitly computing the time dependent molecular polarization from the coupled dynamics of electronic and vibrational degrees of freedom in a mean-field mixed quantum-classical approximation. The oligomers are described by the Su-Schrieffer-Heeger Hamiltonian and the effect of decoherence is incorporated by propagating an ensemble of quantum-classical trajectories with initial conditions obtained by sampling the Wigner distribution of the nuclear degrees of freedom. The electronic coherence of superpositions between the ground and excited and between pairs of excited states is examined for chains of different length, and the dynamics is discussed in terms of the nuclear overlap function that appears in the off-diagonal elements of the electronic reduced density matrix. For long oligomers the loss of coherence occurs in tens of femtoseconds. This time scale is determined by the decay of population into other electronic states through vibronic interactions, and is relatively insensitive to the type and class of superposition considered. By contrast, for smaller oligomers the decoherence time scale depends strongly on the initially selected superposition, with superpositions that can decay as fast as 50 fs and as slow as 250 fs. The long-lived superpositions are such that little population is transferred to other electronic states and for which the vibronic dynamics is relatively harmonic.     

Simulation of vibrational dephasing of I(2) in solid Kr using the semiclassical Liouville method

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.     

Multipartite quantum entanglement evolution in photosynthetic complexes

We investigate the evolution of entanglement in the Fenna-Matthew-Olson (FMO) complex based on simulations using the scaled hierarchical equations of motion approach. We examine the role of entanglement in the FMO complex by direct computation of the convex roof. We use monogamy to give a lower bound for entanglement and obtain an upper bound from the evaluation of the convex roof. Examination of bipartite measures for all possible bipartitions provides a complete picture of the multipartite entanglement. Our results support the hypothesis that entanglement is maximum primary along the two distinct electronic energy transfer pathways. In addition, we note that the structure of multipartite entanglement is quite simple, suggesting that there are constraints on the mixed state entanglement beyond those due to monogamy.