# Global coherence of quantum evolutions based on decoherent histories: theory and application to photosynthetic quantum energy transport

###### Abstract

Assessing the role of interference in natural and artificial quantum dyanamical processes is a crucial task in quantum information theory. To this aim, an appopriate formalism is provided by the decoherent histories framework. While this approach has been deeply explored from different theoretical perspectives, it still lacks of a comprehensive set of tools able to concisely quantify the amount of coherence developed by a given dynamics. In this paper we introduce and test different measures of the (average) coherence present in dissipative (Markovian) quantum evolutions, at various time scales and for different levels of environmentally induced decoherence. In order to show the effectiveness of the introduced tools, we apply them to a paradigmatic quantum process where the role of coherence is being hotly debated: exciton transport in photosynthetic complexes. To spot out the essential features that may determine the performance of the transport we focus on a relevant trimeric subunit of the FMO complex and we use a simplified (Haken-Strobl) model for the system-bath interaction. Our analysis illustrates how the high efficiency of environmentally assisted transport can be traced back to a quantum recoil avoiding effect on the exciton dynamics, that preserves and sustains the benefits of the initial fast quantum delocalization of the exciton over the network. Indeed, for intermediate levels of decoherence, the bath is seen to selectively kill the negative interference between different exciton pathways, while retaining the initial positive one. The concepts and tools here developed show how the decoherent histories approach can be used to quantify the relation between coherence and efficiency in quantum dynamical processes.

###### pacs:

03.65.-w,03.65.Yz,05.60.Gg,42.50.Lc,03.67.Ac## I Introduction

Coherence is ultimately the most distinctive feature of quantum systems. Finding a proper measure of the coherence present at different times scales in a quantum dynamical system is the first essential step for assessing the role of quantum interference in natural and artificial processes. This is of particular relevance for those quantum evolutions in which information or energy are transformed and transferred in order to achieve a given task with high efficiency. In this context the relevant questions are: how much and what kind of coherence is created vs destroyed by the dynamical evolution? How does coherence determine/enhance the performance of the given process? These are in general difficult questions and to be answered they require an appropriate and sufficiently comprehensive framework. A general and fundamental formalism to describe quantum interference is provided by the decoherent histories (DH) approach to quantum mechanics. DH have mainly found applications to foundational issues of quantum mechanics such as the formulation of a consistent framework to describe closed quantum systems, the emergence of classical mechanics from a quantum substrate, the solution of quantum paradoxes, decoherence theory, quantum probabilities Griffiths (); Griffiths-Book (); Omnes (); GMannHartle (); GellMannHartle (); ZurekHist (); Hartle_NegProb2008 (); GMann-Hartle_2014-1 (). However, DH can also be a systematic tool for quantifying interference in quantum processes, and discussing its relevance therein. Indeed, DH provide a precise mathematical formalization of interference by means of the the so called decoherence matrix . The latter is built on the elementary notion of histories and allows one to describe the quantum features vs the classical ones in tems of interference between histories, or pathways if one resorts to the mental picture of the double slit experiments. It is however difficult to quantify in a compact and meaningful way the content of and its implications for the dynamics of specific systems. Our first main goal is therefore to define and test appropriate measures allowing for the investigation of how interference can determine the performance of a given quantum information processing task. Starting from and by its sub-blocks we define different functionals. In particular, we introduce a global measure of coherence able to describe the coherence content of a general quantum evolution at its various time scales; an average (over different time-scales) measure of coherence; and average mesure of interference between histories leading to a specific output.

While the tools we introduce are of general interest and application, in order to test them we apply them to a specific but relevant instance of quantum dynamics taking place in photosynthetic membranes of bacteria and plants: quantum energy transport. Here the basic common mechanism is the following: a quantum excitation is first captured by the system and then migrates through a network of sites (chromophores) towards a target site, e.g., a reaction center, where the energy is transformed and used to trigger further chemical reactions. There is now an emerging consensus that efficient transport in natural and biologically-inspired artificial light-harvesting systems builds on a finely tuned balance of quantum coherence and decoherence caused by environmental noise Ishizaki-PCCP (); Ishizaki-PNAS (); Belcher-Virus (); Kevin-Bardeen_Original (); GiordaJCP (), a phenomenon known as environment-assisted quantum transport (ENAQT). This paradigm has emerged with clarity in recent years, as modern spectroscopic techniques first suggested that exciton transport within photosynthetic complexes might be coherent over appreciable timescales Engel (). Indeed, a growing number of experiments has provided solid evidence that coherent dynamics occurs even at room temperature for unusually long timescales (of the order of ) Collini (); Panit (). Efforts to describe these systems have led to general models of ENAQT Kevin-Bardeen_Original (); Mohseni (); Mohseni2 (); Plenio (); Cao (), depicting the complex interplay of three key factors: coherent motion, i.e., quantum delocalization of the excitation over different sites, environmental decoherence, and localization caused by a disordered energy landscape. So far, the presence of coherence in light-harvesting systems has been qualitatively associated to the observation of distinctive ‘quantum features’. Originally, coherence was identified with ‘quantum wavelike’ behavior as reflected by quantum beats in the dynamics of chromophore populations within a photosynthetic complex. Later works, employing quantum-information concepts and techniques, have switched attention towards quantum correlations between chromophores, in particular quantum entanglementCaruso (); Whaley (); GiordaLHCII (); Ishizaki (). Besides being open to criticism (see, e.g., Miller (); Tiersch ()), these approaches do not provide direct quantitative measures of coherence in the presence of noise. Therefore, in what follows, we shall apply the novel tools based on DH to a simple yet fundamental model of quantum energy transfer. We will focus on a relevant trimeric subunit of the Fenna-Matthews-Olson (FMO) complex, the first pigment-protein complex to be structurally characterized Blankenship (). The trimer is virtually the simplest paradigmatic model retaining the basic charcteristics of a disordered transfer network and it can also be conceived as an essential building block of larger networks. For simplicity, we will use the well-known Haken-Strobl model Haken () to describe the interplay between Hamiltonian and dephasing dynamics. While the model is an oversimplified description of the actual dynamics taking place in real systems, it allows to spot out the essential features that may determine the high efficiency of the transport. We shall initially focus on a new coherence measure , based on the decoherence matrix, and characterize its behavior verifying that it can consistently identify the bases and timescales over which quantum coherent phenomena are present during the evolution of the system. We shall then show how the average coherence exhibited on those time scales can be connected with the delocalization process. A more detailed analysis will be aimed at distinguishing between constructive and destructive interference affecting the histories ending at the site where the excitation exits the photosynthetic structure. By using the decoherence functional, we will show that the beneficial role of dephasing for the transport efficiency lies in a selective suppression of destructive interference, a fact that has been systematically suggested in the literature, but never expressed within a general and comprehensive framework that allows the quantitative evaluation of coherence and its effects.

The application of the introduced tools and methods based on DH to a simple yet paradigmatic system shows how one can properly quantify the coherence content of a complex quantum dynamics and elucidate the role of coherence in determining the overall efficiency of the process Pleniocoherence ().

The paper is structured as follows. In Section II we review the basic decoherent histories formalism. In Section III we define the measure of coherence and describe its meaning and properties. In Section IV we first introduce the used model for describing the energy transport in the selected trimeric complex. We then discuss the coherence properties of the excitonic transport: By means of the appropriate measures based on the decoherent histories formalism we identify the essential features that may determine the high efficiency of the transport. In Section V we briefly discuss how to extend our results to the whole FMO complex. In Section VI we summarize our results an draw our conclusions.

## Ii Decoherent histories

The formalism of decoherent (or consistent) histories was developed in slightly different flavors by Griffiths Griffiths (); Griffiths-Book (), Gell-Mann GMannHartle (); GMann-Hartle_2014-1 (), HartleHartle_NegProb2008 () and Omnès Omnes (). DH provide a consistent formulation of quantum mechanics where probabilities of measurement outcomes are replaced by probabilities of histories. In this formulation, external measurement apparatuses are not needed, and then one does not need to postulate a “classical domain" of observers. As a consequence, quantum mechanics becomes a theory that allows the calculation of probabilities of sequences of events within any closed system, including the whole universe, without the necessity of invoking postulates about the role of measurement. In this framework, the “classical domain" can be seen to emerge as the description of the system becomes more and more coarse-grained.

The idea of ‘histories’ stems from Feynman’s ‘sum-over-histories’ formulation of quantum mechanics. As is known, any amplitude between an initial and a final state can be expressed as a sum over paths, or histories: Upon inserting the identity decomposition at differerent times we get

where we use the Heisenberg notation . Thus the total amplitude is decomposed as a sum of amplitudes, each one corresponding to a different history identified by a sequence of projectors .

The decoherent histories formalism assumes that histories are the fundamental objects of quantum theory and gives a prescription to attribute probabilities to (sets of) histories. A history is defined as a sequence of projectors at times . Probabilites can be assigned within exhaustive sets of exclusive histories, i.e., sets of histories where subscripts label different alternatives at times . Histories are exhaustive and exclusive in the sense that the projectors at each time satisfy relations of orthogonality , and completeness, . In other words, the projectors define a projective measurement. Within a specified set, any history can be identified with the sequence of alternatives realized at times .

Different alternative histories can be grouped together with a procedure called coarse-graining. Starting from histories and we can define a new, coarse-grained history by summing projectors for all times such that and differ:

for all . By iterating this procedure, one can obtain
more and more coarse-grained histories. A special type of coarse-graining is the temporal coarse-graining:
we group together histories such that
such that at some time we have .
Then the coarse-grained history
contains only one-projector (equal to the identity) at time ,
that can be neglected and hence removed from the string of projectors
defining the history. On the other hand temporal fine-graining
can be implemented for example by allowing different alternatives
at a times . In particular, one can
create new sets of histories
from a given one by adding different alternatives
at time ; the sets are fine grained versions
of the sets .

Once we specify the initial state and the (unitary) time
evolution , we can assign any history a weight

where we use the Heisenberg notation . When the initial state is pure, and the final projectors are one dimensional, , this formula takes the simple form of a squared amplitude

(1) |

Weights cannot be interpreted as true probabilities, in general. Indeed, due to quantum interference between histories, the do not behave as classical probabilities. Indeed, consider two exclusive histories and the relative coarse-grained history by: . If the were real probabilities, we would expect . Instead, what we find is

Due to the non-classical term , representing quantum interference between the histories and , the classical probability-sum-rule is violated. The matrix

(2) |

is called decoherence functional or decoherence matrix. The decoherence matrix can be thought of as a “density matrix over histories”: Its diagonal elements are the weights of histories and its off-diagonal elements are interferences between pairs of histories. The decoherence matrix has the following properties: i) it is Hermitian ii) it is semipositive definite iii) it is trace one iv) it is block-diagonal in the last index, . Weights of coarse-grained histories can be obtained by summing matrix entries in an block of the decoherence matrix corresponding to the original fine-grained histories. For istance, the weight of history is obtained by summing entries of a block of the decoherence matrix:

A necessary and sufficient condition to guarantee that the probability sum rule apply within a set of histories is

This condition is termed as weak decoherenceGellMannHartle (); the necessary and sufficient condition that is typically satisfied GellMannHartle () and that we will adopt in the following is the stronger one termed as medium decoherence

(3) |

Medium decoherence implies weak decoherence. Any exhaustive and set
of exclusive histories satisfying medium decoherence is called a decoherent
set. The fundamental rule of DH approach is that probabilities can
be assigned within a decoherent set, each history being assigned a
probabability equal to its weight. If medium decoherence holds, the
diagonal elements of the decoherence matrix can be identified as real
probabilities for histories and we can write .

Due to property iv), if we perform a temporal coarse-graining over
all times except the last, we obtain ‘histories’ with only one projection,
at the final time . These histories automatically
satisfy medium decoherence:

where is the probability that the system is in at time . Due to interference, the probability of being in at time is not simply the sum of probabilities of all alternative paths leading to , i.e, of all alternative histories with final projection . In formulas,

The probability and the global interference of histories ending in can be thus expressed as

(4) |

with . Destructive interference will happen when , constructive interference when .

The decoherent histories formalism is consistent with and encompasses the model of environmentally induced decoherenceZurekHist (). Given a factorization of the Hilbert space into a subsystem of interest and the rest (environment), , the events of a history take the form where and are projectors onto Hilbert subspaces of and respectively. Histories for alone can be obtained upon considering appropriate coarse-grainings over the degrees of freedom of the environment, such that the events are where is the identity over . Upon introducing the time-evolution propagator as , we can rewrite the decoherence matrix as:

(5) | ||||

If the initial state is factorized , then the reduced density matrix evolves according to where is the (non-unitary) reduced propagator defined by

(6) |

If the evolution of the system and environment is Markovian, we can write . As proved by Zurek ZurekHist (), under the assumption of Markovianity we can rewrite the decoherence matrix in terms of reduced quantities alone, i.e., quantities pertaining to the system only:

(7) |

That is, the model of environmentally induced decoherence can be obtained by applying the decoherent histories formalism to system and environment together, and by coarse-graining over the degrees of freedom of the environment.

## Iii The coherence measure

The DH approach provides the most fundamental framework in which the transition from the quantum to the classical realm can be expressed. Indeed, it is based on the most basic feature characterizing the quantum world: interference and the resulting coherence of the dynamical evolution. Despite being a well developed field of study, the DH history approach lacks for a proper global measure of the coherence produced by the dynamics at the different time scales. We therefore introduce a measure that quantifies the global amount of coherence within a set of histories. Assume projectors for all times are taken in a fixed basis , . Assume further that histories are composed by taking equally spaced times between consecutive projections i.e., . (in other words, histories correspond to projections applied in the same basis and repeated at regular times). For such a set of histories, consider the decoherence matrix

where . Take the von Neumann entropy of the decoherence matrix,

(8) |

Due to coherence between histories, differs from the ‘classical-like’ Shannon entropy of history weights

(9) |

where are the diagonal elements of , i.e., the weights. The difference between the two quantities is wider if off-diagonal elements of the decoherence matrix are bigger, i.e., if the set of histories is more coherent. Let us define:

(10) |

We argue that is suitable to be used as a general
measure of coherence within the set of histories defined by .
Indeed, we can readily prove the following properties:

i) . is obvious. To prove , let
us define a matrix
where off-diagonal entries are set to zero. Since

we obtain that the numerator of (10) can be expressed as a quantum relative entropy:

where is the relative entropy between and .

ii) iff ,
i.e., vanishes if medium decoherence holds for
the set of histories, since the two quantities
and coincide in this case.

Thus is in essence a (statistical) distance between the decoherence matrix and the corresponding diagonal matrix , renormalized so that its value lies between and . The greater are the off-diagonal elements of , the greater the distance. The meaning of can be easily understood if we use the linear entropy, a lower bound to the logarithmic version:

In this case, we obtain a ‘linear entropy’ proxy of as:

which is a simplified version that, by avoiding the diagonalization of , helps containing the numerical complexity.

The measure introduced is well grounded on physical considerations. In the following we will apply it to a simple system in order to check its consistency, and later use it to characterize the coherence properties of the evolutions induced by various regimes of interaction with the environment. First, one has to check whether the measure properly takes into account the action of the bath. In particular, if the bath is characterized by a decoherence time , it is known (ZurekHist ()) that on time scales the decoherence matrix becomes diagonal: The probability of a history at time can be fully determined by its probability at time , since no interference can occur between different histories. Indeed, the action of the bath is to create a decoherent set of histories that are defined by a proper projection basis: the pointer basis (ZurekHist ()). Therefore, the fine-graining procedure obtained by constructing a set of histories via the addition of a new complete set of projections in the same basis at time to the set , should leave the coherence functional invariant, i.e., . If instead the same fine-graining procedure should lead to .

Before passing to analyze a specific system, we want to focus on the complexity of the evaluation of and . The dimension of the decoherent matrix grows with the dimension of the basis and the number of time instants that define each history as . This exponential growth in principle limits the application of the DH approach to small systems. However, as for the system considered in this paper the computational effort is contained due to the small number of subsystems (chromophores) and the small dimension of the Hilbert space which is limited to the single-exciton manifold. As we shall see, by limiting the choice of to a reasonable number, the analysis can be fruitfully carried even on a laptop.

## Iv Trimer

We now start to analyze decoherent histories in simple models of energy transfer comprising a small number of chromophores (sites). Neglecting higher excitations, each site can be in its ground or excited state. We work in the single-excitation manifold, and define the site basis as

i.e., state represents the exciton localized at site . On-site energies and couplings are represented by a Hamiltonian that is responsible for the unitary part of the dynamics. Interaction with the environment is implemented by the Haken-Strobl model, that has been extensively used in models of ENAQT Mohseni (); Plenio (); GiordaLHCII (); Cao (). The effect of the environment is represented by a Markovian dephasing in the site basis, expressed by Lindblad terms in the evolution, as follows:

(11) |

where are projectors onto the site basis, and are the (local) dephasing rates. Furthermore, site can be incoherently coupled to an exciton sink, represented by a Linblad term

where and is the trapping
rate. Contrary to other works, we neglect exciton recombination, as
it acts on much longer timescales than dephasing and trapping.

The global evolution is Markovian and can be represented by means
of the Liouville equation

(12) |

that can be simply solved by exponentiation,

(13) |

In the notation above, the propagator has the form . The efficiency of the transport can be evaluated as the leak of the population of the exit site towards the sink:

(14) |

The overall efficiency of the process is obtained by letting

While its Markovianity limits the faithful description of decoherence processes actually taking place in real photosynthetic systems, the model retains the basic and commonly accepted aspects of decoherence, that acts in the site basis: albeit in a complex non-Markovian way, the protein enviroment measures the system locally (i.e., on each site), thus destroying the coherence in the site basis and creating it in the exciton basis. Note that the formalism can also be applied to a ‘dressed’ or polaronic basis where we include strong interactions between chromophores and vibrational modes. That is, to apply the DH method, one only needs a model in which an exciton hops between sites, dressed or undressed. The model is therefore suitable to readily implement the decoherent histories paradigm and to spot the main basic features we are interested in and that are at the basis of the success of ENAQT.

The FMO unit has chromophores and a complex energy and coupling landscape with no symmetries. Energies and couplings (i.e., the Hamiltonian ) can be obtained by different techniques: They can be extracted by means of 2D spectroscopy as in Cho () or computed through ab initio calculations as in Renger (), with similar but not exactly equal results. This very complex struxture makes FMO far from ideal as a first example to study. We thus prefer to start by working with a much simpler, yet fully relevant subsystem: the trimeric unit composed by the sites and of the FMO complex in the notation of Renger (); Cho ()). The first chromophore is the site in which the energy transfer begins, while the third chromophore is the site from which the excitation leaves the complex. The Hamiltonian of the trimeric subunit is Renger ()

(15) |

The eigenenergies of the system ar given by
which yields the eigenperiods
.
Due its structure, the trimer is a chain composed by a pair of chromophores
(), degenerate in energy and forming a strongly coupled dimer,
and a third chromophore moderately coupled with the second one only.
Since in the following we suppose that the exciton starts from site
, we expect a prominent role of the dimer in the dynamics, at least
in the first tens of femtoseconds.

In order to show how the DH analysis can be implemented, in the following we are going to consider histories in the site and the energy bases, with projections at times . We first use the coherence function introduced above (10) to evaluate the global coherence of the exciton transport process. In order to test the behavior of for different values of dephasing, in Fig 1 we first plot as a function of the time interval between projections for two values of the dephasing rate: , corresponding to the full quantum regime (Fig.1) for the site basis (a); corresponding to an intermediate value of dephasing (Fig. 1) for the site basis (b) and the energy basis (c).

Before entering the discussion of the various regimes, we note that as a function of the number of projections all curves display the expected behavior: The increase (decrease) of the number of projections corresponds to a temporal fine-graining (coarse-graining) of the evolution; therefore, an increase (decrease) of should imply an increase (decrease) of the amount of coherence between histories. As shown in Fig. (1) the function correctly reproduces the fine (coarse) graining feature: the qualitative behavior of as a function of is not affected by the choice of , while an increase of corresponds, at fixed , to an increase of . We will therefore use in the following the value that allows for a neat description of the phenomena and for a reasonable computational time.

As for the behavior at fixed , we have that in the full quantum regime (), the system obviously displays coherence in the site basis only since

and the decoherence matrix in the energy basis is diagonal and independent on and . This simply means that in the full quantum regime histories in the exciton basis are fully decohered, since the system is not able to create coherence among excitons. Still in the full quantum regime, in the site basis, the coherence oscillates as the exciton, starting at site , goes back an forth along the trimer, and the evolution builds up coherence in this basis, see Fig.1 (a). In this regime, the trimer can be approximately seen as a dimer composed by the first two chromophores, and the exciton performs Rabi oscillations with a period given by ; oscillates with half the period: for the exciton is migrated mostly on site and has a minumum – which is different from zero since the exciton is partly delocalized on site , and the system therefore exhibits a non vanishing coherence.

For intermediate values of , Fig.1 (b), the coherence in site basis as measured by correctly drops down at ZurekHist (). The dephasing has a strong and obvious effect on the coherence between pathways: coherence in this basis is a monotonically decreasing function of . This is well highlighted by the global coherence function , whose maximal values are reduced by a factor of with respect to those corresponding to full quantum regime. After a time the histories are fully decohered. Indeed, due to the specific model of decoherence (11), which amounts to projective measurements on at each site with a rate , the system kills the coherence in the site basis, which in turn corresponds to the stable pointer basis for this model ZurekHist (), i.e., the basis in which the density matrix is forced to be diagonal by the specific decoherence model. On the other hand, and for the same reason, the dynamics starts to build up coherence in the exciton basis , see Fig. 1(c). However, this coherence is later destroyed - on a time scale of approximately - since the stationary state of the model is the identity. This effect is even more evident if one compares the behavior of in the exciton basis for different values of , as shown in Fig. 1(d): grows with and it lasts over longer time scales. This feature is coherent with the expectations: the equilibrium state for high is the identity. Due to the projections implemented by the environment in the site basis, the system is forced to create coherence in the exciton basis. When is very high a quantum Zeno effect takes in, the dynamics is blocked, and the time required to reach the equilibrium, and to destroy coherences in all bases, consequently grows.

This first analysis therefore shows that is indeed a good candidate for assessing the global coherence properties of quantum evolutions. For a fixed number of projections , can be interpreted as a measure of the global coherence exhibited by the dynamics over the time scale .

We now analyze in detail the specific features of quantum transport for the trimer. The dynamics starts at site and evolves by delocalizing the exciton on the other chromophores. In order to study this process, we first use a measure of delocalization introduced inGiordaLHCII () for the study of LHCII complex dynamics:

(16) |

that is simply the Shannon entropy of , the populations of the three chromophores. This measure allows one to follow how much the exciton gets delocalized over the trimer with time and in different dephasing situations: , i.e., is zero when the exciton is localized on a chromophore and it takes its maximal value when the population of the three sites are equal. In Fig. 2 we plot both and the population of site for different values of . Due to the presence of interference, in the mainly quantum regime (), the exiton first delocalizes mainly over the dimer and partly on the third site: The first maximum corresponds to when the system builds up a (close to uniform) coherent superposition between sites and , while a non negligible part of the exciton is found in site ; indeed , the last value corresponding to i.e., to a uniform superposition over the sites and only. As the dynamics of the sytems extends to later times we see that and have an oscillatory behavior, whose main period is , and which approximately corresponds to Rabi oscillations between site and , although the initial state fully localized in site cannot be rebuilt due to the presence of site . As for the transport, we see that in this regime the system cannot take advantage of the initial fast and high delocalization: the exciton bounces back and forth over the trimer. In the intermediate regime , due, as we will later see, to the selective suppression of interference processes, the initial speed up in delocalization is sustained by the dynamical evolution, and the transfer rate to site is correspondingly increased. For very high values of decoherence () the role of initial interference is suppressed and the initial speed-up disappears: the environment measures the system in site basis at high rates and the delocalization process is highly reduced.

The optimal delocalization occurs in correspondence of and it can be interpolated with a double exponential function

with . The first time scale describes the initial fast quantum delocalization process described above, while the second time scale the slower subsequent delocalization and the reaching of the equilibrium situation, .

We now pass to systematically analyze the behavior of the coherence of the evolution with respect to the strength of the interaction with the environment and its relevance for the energy transport process. As a first step we plot both and for different values of , Fig. (3). The plots show that the coherence function exhibits the required behavior: For small , oscillates with period , following the Rabi oscillations of the dimer. The minima occur at , showing that the exciton is “partially” localized on site or , and partially delocalized on site . As grows, the system becomes unable to create coherence on large time scales; the decay of is mirrored by a the reduction of the amplitude in the oscillations of .

We now focus on the relevant time scales for the initial fast delocalization process highlighted by our previous analysis, which are of the order of tens to hundreds of femtoseconds. We therefore introduce the following average measure of global coherence of the evolution

is the average of the coherence exhibited by the dynamics of the system at the time scales . In Fig. 4 we show for the trimer (15) in the site basis for different values of . We first focus on the behavior of for values of dephasing in the range . In this range, for small time scales to the average global coherence is approximately constant and equals the value attained in the full quantum regime i.e., . For larger time scales () rapidly decreases with . This analysis shows that the behavior of matches the expectations: the higher the smaller the time scales over which decoherence takes place, the lower the global coherence of the dynamics. Along with the functional is therefore in general a good candidate for the evaluation of the global coherence of open quantum systems evolution. As for the transport dynamics, we focus on the timescale identified with the analysis of for optimal dephasing; for and we see that the system indeed retains most of the average coherence of the purely quantum regime up to the optimal values of decoherence ( in the figure), losing it afterwards; this is a clear indication that this phenomenon is at the basis of the the fast initial delocalization process. Over longer time scales, the relevance of coherence is highly suppressed.

We now deepen our analysis about the relevance of the coherence of the evolution for the energy tranfer efficiency. To this aim we focus on the basic feature that distingushes the classical and the quantum regime: interference. In particular we focus on the sub-block of the decoherence matrix pertaining to the third chromophore, which describes the set of histories in site basis ending at site . Due to interference the probability of occupation of the site at time can be written in terms of the the histories ending at site , see (4) .

In Fig. 5 we show for different values of dephasing. One has different regimes: for , the set of histories in site basis is fully decohered; , the histories do not interfere with each other and i.e, the probability is simply the sum of diagonal elements of . In the mainly quantum regime , : after the initial positive peak the histories interfere with each other, globally the interference is mostly negative and therefore . For intermediate values of decoherence the interference has a positive peak and then reduces to zero. While the first initial fingersnap of positive interference that takes place in the first is common for all curves corresponding to small and intermediate values of , the main effect of the bath is displayed after this initial period of time: the decoherence gradually suppresses interference, both the positive and the negative one; however, for intermediate values of the effect is stronger as for the negative part of the interference patterns. The environment thus implements what can be called a quantum recoil avoiding effect: it prevents the part of the exciton that – thanks to constructive interference – has delocalized on site to flow back to the the other sites.

In order to evaluate a possible advantage provided by the initial speed-up in the delocalization process and by the interference phenomena showed above one has to take into account another relevant time scale of the transport process: the trapping time. Indeed, if the system is to take advantage of the fast delocalization due to the coherent behavior, the exit of the exciton should take place on time scales of the order of the delocalization process. The theoretical and experimental evidences show that this is the case: the trapping time for the FMO complex is estimated in the literature to be of the order of i.e., the exit of the exiton starts soon after the fast delocalization due to quantum coherence has taken place. The role of the interference between paths, in particular those leading to site , can therefore be appreciated by numerically evaluating

(17) |

i.e., the average over the trapping time scale of of the total (), negative () and positive () average interference between the histories ending in site , with . In particular, in Fig. 6(a) the different kinds of interference are plotted for histories terminating at site : on average, the negative interference highly reduces the total interference for small values of decoherence strength; when , vanishes, the average total interference equals the positive one , and it is maximal for values of comparable to those that maximize (). In Fig. 6(b), we compare the behavior of for all sites. The results again suggest that decoherence acts on the interference provided by the quantum engine in order to favor the flow of the exciton towards the exit chromophore: the average positive interference between histories ending at sites and grows in modulus with and attains a maximum for intermediate values of decoherence; while the average negative interference between histories ending at site 1 decreases and attains a minimum for intermediate values of . The combined effect of decoherence and interference thus helps depopulating site 1 and populating site 2 and 3.

We can now tackle one of the most relevant aspects of our discussion: the net effect of the above described phenomena on the overall efficiency of the transport. The latter can be fully appreciated by evaluating the efficiency of the process (14) and by recognizing that, in the decoherent histories language, it can be expressed as:

where and . This split allows one to appriciate the role of interference for the efficiency. In Fig. (7) is plotted for different values of dephasing. In agreement with what discussed above, we have three regimes: for very small values of the overall efficiency is poor; this is due to the presence of high negative interference that in average prevents the exciton to migrate to the exit site. For large values of the interference processes are completely washed out and the system cannot take advantage of the fast quantum delocalization. For intermediate (optimal) values of only the negative interference has been washed out: is positive, it acts on short time scales, and it provides on average an enhancement of the global efficiency.

These results, within the limits of the simple model of decoherence taken into account, undoubtedly show for the first time that the so called ENAQT phenomenon can well and properly be understood both qualitatively and quantitatively within the decoherence histories approach, i.e., in terms of very the basic concepts of coherence and interference between histories. The often recalled “convergence” of time scales or “Goldilocks” effect (Goldilocks ()) in biological quantum transport systems seems therefore to be well rooted in the processes discussed above: if decoherence is too small the system shows both positive and negative interference (see Fig. 2), the delocalization has an ocillatory behavior, and the exciton bounces back and forth along the network thus preventing its efficient extraction. If instead decoherence is very high one has that the complete washing out of intereference and coherence implies the delocalization process to be very slow, no matter how fast the trapping mechanism try to suck the exciton out of the system. In order to take advantage of the effects of quantum coherent dynamics: the bath must act on the typical time scales of quantum evolution in order to implement the quantum recoil avoiding process; the extraction of the exciton from the complex, characterized by , must then start soon after the initial fast delocalization has taken place. Should the extraction take place on longer time scales, the benefits of the fast initial delocalization would be spoiled: Waiting long enough, the system would eventually reach together with equilibrium a decent delocalization even for moderately high values of , but in this case the transfer would be obviously much slower.

## V Fmo

The above arguments can be easily applied to the whole FMO complex. Fig. 8 and 9 show the application of the decoherent histories method to excitonic transport in FMO. The main features of the behavior of and are maintained although obvious differences can be found since the dynamics in now determined by the interplay of different eigenperiods and interference paths are more complex. In particular, Fig.8(b), one can observe a revival of positive interference for small values of , that does enhance the efficiency for , Fig. 9(a); but this is not sufficient to compensate the initial and subsequent negative interference, thus impeding the reach of optimal values of . In general, compared to the trimer and as suggested by Fig. 9(a) the maximum average positive coherence on short time scales is attained for smaller values of . The overall picture is not significantly affected if one decides to start the dynamics from site instead of site , as it often is reported in the literature.

## Vi Conclusions

The decoherent histories approach provides a general theory to study
the distinctive feature exhibited by quantum systems: coherence. However,
despite its generality and foundational character, in order to measure
the effects of coherence and decoherence the DH approach needs to
be complemented with a quantitative way to condense the information
contained in the basic object of the theory, i.e., the decoherence
matrix . In this paper we introduce a set of tools that
allow one to assess the (global) coherence properties of quantum (Markovian)
evolution and that can be used to relate the coherence content of
a general quantum dynamical process to the relevant figure of merits
of the given problem. We first define the coherence functional
, that can be interpreted as a measure of the global
coherence exhibited by the dynamics in the basis over the time
scale . While this measure is completely general,
one can further introduce other relevant tools tailored to the specific
system and type of system-environment interaction at hand. We thus
focus on a simple yet paradigmatic model of environmentally assisted
energy transfer where coherence effects have been shown to play a
significant role in determining the efficiency of the process: a trimeric
subunit of the Fenna-Matthews-Olson photosynthetic complex. Based
on and we define: a measure
able to characterize the average coherence
exhibited by the dynamics of the system over the time scales
for a fixed value of the dephasing ; a measure of the
average interference occurring between the histories
ending at a given “site” .

Within the specific model, we
first thoroughly assess the consistency of the behavior of
in the various regimes. We then show how the introduced tools allow
to study the intricate connections between the efficiency of the transport
process and the coherence properties of the dynamics. In particular
we show that the delocalization of the exciton over the chromophoric
subunit is strongly affected by the amount of (average) coherence
allowed by the interaction with the bath in the first tens to hundreds
of femtoseconds. If the system-bath interaction is too strong, coherence
is suppressed alongside the interference between different histories,
in particular those ending at the site where the excitation leaves
the complex. If the interaction is too weak the system exhibits high
values of coherence even on long time scales, but it also exhibits
negative interference between pathways ending at the exit site, a
manifestation of the fact that the exciton bounces back and forth
over the network thus preventing its efficient extraction. In the
intermediate regime i.e., when the different time scales of the system
(quantum oscillations, decoherence and trapping rate) converge, the
system shows high values of coherence on those time scales. The action
of the bath has a quantum recoil avoiding effect on the dynamics
of the excitation: the benefits of the fast initial quantum
delocalization of the exciton over the network are preserved and sustained
in time by the dynamics; in terms of pathways leading to the exit
site, the action is to selectively kill the negative interference
between pathways, while retaining the initial positive one. These
effects can be explicitly connected to the overall efficiency of the
environment-assisted quantum transport: the gain in efficiency
for intermediate (optimal) values of decoherence can thus be traced
back to the basic concepts of coherence and interference between pathways
as expressed in the decoherent histories language.

While the specific
decoherence model used (Haken-Strobl) is an oversimplified description
of the actual dynamics taking place in real systems, we believe
that our analysis allows to spot out the essential features that may
determine the high efficiency of the transport even in more complex
system-environment scenarios.

In conclusion, the tools introduced
in this paper allow to thoroughly assess the coherence properties
of quantum evolutions and can be applied to a large variety of quantum
systems, the only limits being the restriction to Markovian dynamics
and the computational efforts required for high dimensional systems.
However, the extension to non-Markovian realms is indeed possible Allegraprep (),
and the use of parallel computing may allow the treatement of reasonably
large systems.

###### Acknowledgements.

P. Giorda and M. Allegra would like to thank Dr. Giorgio “Giorgione” Villosio for his friendship, his support, his always reinvigorating optimism, and his warm hospitality at the Institute for Women and Religion - Turin, where this paper was completed (cogitato, mus pusillus quam sit sapiens bestia, aetatem qui non cubili uni umquam committit suam, quin, si unum obsideatur, aliud iam perfugium elegerit).P. Giorda would like to thank Prof. A. Montorsi, Prof. M.G.A. Paris and Prof. M. Genovese for their kind help.

S. Lloyd would like to thank M. Gell-Mann for helpful discussions.

## References

- (1) R. B. Griffiths, J. Stat. Phys. 36, 219 (1984).
- (2) R. B. Griffiths, Consistent quantum theory, Cambridge University Press (2003).
- (3) M. Gell-Mann and J. B. Hartle, Phys.Rev. D 47, 3345 (1993).
- (4) M. Gell-Mann and J. B. Hartle, in Complexity, entropy, and the physics of information (Addison-Wesley, Reading, Massachusetts, 1990); in Proc. of the 25th International Conference on High Energy Physics, Singapore (World Scientific, Singapore, 1990).
- (5) M. Gell-Mann and J. B. Hartle. Physical Review A 89, 052125 (2014).
- (6) J. B. Hartle, Phys. Rev. A 78, 012108 (2008).
- (7) W. H. Zurek, J. P. Paz, Phys. Rev. D 48, 2728 (1993).
- (8) R. Omnès, J. Stat. Phys. 53, 893, 933, 957 (1988), 57, 357 (1989); Rev. Mod. Phys. 64, 339 (1992);
- (9) A. Ishizaki, T. R Calhoun,G. S. Schlau-Cohen, and G. R. Fleming, Phys. Chem. Chem. Phys. 12, 7319 (2010).
- (10) A. Ishizaki, and G. R. Fleming, Proc. Natl. Acad. Sci. USA 106, 17255 (2009).
- (11) H. Park, N. Heldman, P. Rebentrost, L. Abbondanza, A. Iagatti, A. Alessi, B. Patrizi, M. Salvalaggio, L. Bussotti, M. Mohseni, F. Caruso, H. C. Johnsen, R. Fusco, P. Foggi, P. F. Scudo, S. Lloyd & A. M. Belcher, Nature Materials (2015).
- (12) L. Banchi, G. Costagliola, A. Ishizaki, and P. Giorda, J. Chem. Phys. 138, 184107 (2013).
- (13) K. M. Gaab and C. J. Bardeen, J. Chem. Phys. 121, 7813 (2004).
- (14) R.E. Blankenship, Molecular Mechanisms of Photosynthesis, Blackwell Science Ltd., London (2002).
- (15) G.S. Engel, T. R. Calhoun, E. L. Read, T. K. Ahn, T. Mancal, Y. C. Cheng, R. E. Blankenship and G. R. Fleming, Nature 446, 782 (2007).
- (16) E. Collini, C. Y.Wong, K. E.Wilk, P.M. Curmi, P. Brumer, and G. D. Scholes, Nature 463, 644 (2010).
- (17) G. Panitchayangkoon, D. Hayes, K. A. Fransted, J. R. Caram, E. Harel, J. Wen, R. E. Blankenship, G. S. Engel, Proc. Nat. Acad. Sci. 107, 12766 (2010).
- (18) M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru-Guzik, J. Chem. Phys. 129, 174106 (2008).
- (19) P. Rebentrost, M. Mohseni, I. Kassal, S. Lloyd, and A. Aspuru- Guzik, New J. Phys. 11, 033003 (2009).
- (20) M.B. Plenio and S.F. Huelga, New J. Phys. 10, 113019 (2008).
- (21) J.M. Moix, M. Khasin, J. Cao, New Journal of Physics 15, 085010 (2013)
- (22) S. Lloyd, M. Mohseni, A. Shabani, H. Rabitz, arXiv:1111.4982 (2011).
- (23) H. Haken, G. Strobl, in The Triplet State, Proceedings or the International Symposium, Am. Univ. Beirut, Lebanon (1967), A.B. Zahlan, ed., Cambridge University Press, Cambridge (1967); H. Haken, P. Reineker, Z. Phys. 249, 253 (1972).
- (24) M. Cho, H. M. Vaswani, T. Brixner, J. Stenger, and Gr. R. Fleming, J. Phys. Chem. B 109, 10542-10556 (2005).
- (25) J. Adolphs and T. Renger, Biophysical Journal 91, 2778-2797 (2006).
- (26) F. Caruso, A.W. Chin, A. Datta, S.F. Huelga, M.B. Plenio, Phys. Rev. A 81, 062346 (2010).
- (27) M. Sarovar, A. Ishizaki, G.R. Fleming, K.B. Whaley, Nat. Phys. 6, 462 (2010).
- (28) P. Giorda, S. Garnerone, P. Zanardi, S. Lloyd, arXiv:1106.1986 (2011).
- (29) A. Ishizaki and G. R. Fleming, New J. Phys. 12, 055004 (2010).
- (30) M. Tiersch, S. Popescu and H. J. Briegel, Phil. Trans. R. Soc. A, 370, 1972, 3771-3786 (2012).
- (31) W. H. Miller, J. Chem. Phys. 136, 210901 (2012).
- (32) for a definion and discussion of coherence in quantum states, see T. Baumgratz, M. Cramer, and M. B. Plenio Phys. Rev. Lett. 113, 140401 (2014) and references therein.
- (33) M. Allegra, P. Giorda and S. Lloyd, in preparation.