Non-classicality of the molecular vibrations assisting exciton energy transfer at room temperature.

Advancing the debate on quantum effects in light-initiated reactions in biology requires clear identification of non-classical features that these processes can exhibit and utilize. Here we show that in prototype dimers present in a variety of photosynthetic antennae, efficient vibration-assisted energy transfer in the sub-picosecond timescale and at room temperature can manifest and benefit from non-classical fluctuations of collective pigment motions. Non-classicality of initially thermalized vibrations is induced via coherent exciton-vibration interactions and is unambiguously indicated by negativities in the phase-space quasi-probability distribution of the effective collective mode coupled to the electronic dynamics. These quantum effects can be prompted upon incoherent input of excitation. Our results therefore suggest that investigation of the non-classical properties of vibrational motions assisting excitation and charge transport, photoreception and chemical sensing processes could be a touchstone for revealing a role for non-trivial quantum phenomena in biology.

The experimental demonstration of oscillatory electronic dynamics in light-harvesting complexes12345 has triggered wide-spread interest in uncovering quantum phenomena that may have an impact on the function of the molecular components of living organisms. In general, however, oscillatory patterns in dynamics are not sufficient argument to rule out classical descriptions of the same behaviour. Indeed, recent work has discussed how classical coherence models can predict electronic coherence beating67. Therefore, an important challenge for the growing field of quantum effects in biomolecules is to clearly identify which quantum features with no classical counterpart may manifest in these systems and how they may influence the process of interest.

The question of non-classicality of the dynamics of electronic excitations in light-harvesting systems has been addressed by investigating Leggett–Garg inequalities8. This work concludes that, under Markovian evolution, temporal correlations of the observables of individual pigments should violate classical bounds, and in consequence certain classical theories are unsuitable to describe electronic dynamics. Other works have investigated the quantumness of the electronic degrees of freedom910111213. Despite these efforts, it is still far from understood which non-classical phenomena are directly correlated with efficient energy distribution in a prototype light-harvesting system.

What is clear is that exciton energy transport depends not only on the topology of electronic couplings among pigments but is critically determined by exciton–phonon interactions: molecular motions14 and environmental fluctuations141516 drive efficient transport processes in light-harvesting antennae. In fact, it is well known that exciton–phonon interactions in these complexes have a rich structure as a function of energy and generally include coupling to both continuous and discrete modes associated to low-energy solvated protein fluctuations and underdamped intramolecular vibrations, respectively14. Moreover, evidence is mounting that the interaction between excitons and underdamped vibrations whose energies commensurate exciton splittings may be at the heart of the coherence beating probed in two-dimensional (2D) photon echo spectroscopy17181920212223. Although some insights into the importance of such resonances can be gained from Förster theory24, the wider implications for optimal spatio-temporal distribution of energy192526, for modulation of exciton coherences19202122 and for collective pigment motion dynamics23 have just recently started to be clarified.

The current state of the debate then suggests that a conceptual advance in understanding non-trivial quantum phenomena assisting excitation transport could emerge precisely from investigating non-classical features of the molecular motions and phonon environments that play such a key role. Techniques able to manipulate vibrational states2728 and probe their quantum properties2930 may indeed provide the experimental platform to address this issue. Here we investigate this question in a prototype dimer ubiquitous in light-harvesting antennae of cyanobacteria31, cryptophyte algae3233 and higher plants343536 and show that commensurate energies of exciton splitting and underdamped high-energy vibrations allows exciton–vibration dynamics to induce and harness non-classical fluctuations of collective pigment motions for efficient energy transfer. Negative values of the Mandel Q-parameter37 indicating sub-Poissonian phonon occupation fluctuations and, correspondingly, negative regions in a regularized quasi-probability P distribution in phase space3839 unambiguously preclude any classical description of such fluctuations or its correlations with transport. Our results show a potential functional relevance of non-classicality of molecular fluctuations for exciton transport and therefore provide a framework to investigate similar non-trivial quantum phenomena in the large variety of biomolecular transport4041, photoreception42 and chemical sensing processes434445 that are known (or hypothesized) to be assisted by unequilibrated vibrational motion.

Results

Characterizing non-classicality

The field of quantum optics has developed a solid framework to quantify the quantum properties of bosonic fields46. It therefore provides excellent conceptual and quantitative tools to investigate non-classicality of the harmonic vibrational degrees of freedom of interest in this work. From the perspective of quantum optics, quantum behaviour with no classical counterpart—that is, non-classicality—arises if the state of the system of interest cannot be expressed as a statistical mixture of coherent states defining a valid probability measure47. This then leads to non-positive values of a phase–space quasi-probability distribution such as the Glauber–Sudarshan P(α)-function47.

where χ(ξ) is the characteristic function of the bosonic quantum state. However, highly singular behaviour of P(α) can make its characterization challenging both theoretically and experimentally. To overcome this, verification of the non-classicality of a quantum state can be performed by constructing a regularized distribution P(α) as the Fourier transform of a filtered quantum characteristic function χ(ξ)3839 as explained in the Methods section. Negativities in this regularized distribution are necessary and sufficient condition of quantum behaviour with no classical analogue38 and offer a significant advantage over other distributions, such as the Wigner distribution, which can be positive for quantum states that are truly non-classical.

As an alternative to phase–space distributions, signatures of non-classicality can be observed in the fluctuations of the bosonic field. For a single mode, negative values of the Mandel’s Q–parameter37 can be a signature of non-classical behaviour. It characterizes the departure of the occupation number distribution P(n) from Poissonian statistics through the inequality

where and denote the first and second moments of the bosonic number operator , respectively. Vanishing Q indicates Poissonian number statistics where the mean of equals its variance as it is characteristic of classical wave-like behaviour—that is, a coherent state of light. For a chaotic thermal state, one finds that Q=>0 indicating that particles are ‘bunched’. A Fock state is characterized by Q<0 indicating that particle occupation is restricted to a particular level. Inequalities involving occupation probabilities of nearest number states can similarly witness non-classical occupation fluctuations48.

In what follows we use the above framework to investigate the non-classical behaviour of the vibrational motions that drive excitation dynamics in prototype dimers present in a variety of antennae proteins of photosynthetic systems.

Prototype dimers and collective motion

We consider a prototype dimer where each chromophore has an excited electronic state of energy ε strongly coupled to a quantized vibrational mode of frequency ωvib much larger than the thermal energy scale KBT and described by the bare Hamiltonians Hel=∑εσ+σ− and Hvib=ωvib(b1†b1+b2†b2), respectively. Inter-chromophore coupling is generated by dipole–dipole interactions of the form Hd−d=V(σ1+σ2−+σ2+σ1−). The electronic excited states interact with their local vibrational environments with strength g, linearly displacing the corresponding mode coordinate, Hel−vib=g∑σ+σ−(b†+b). In the above b†(b) creates (annihilates) a phonon of the vibrational mode of chromophore i, while σ± creates or destroys an electronic excitation at site i. The eigenstates |X› and |Y› of Hel+Hd−d denote exciton states with energy splitting given by ΔE=((Δε)2+4V2)1/2 and Δε=ε1−ε2. Transformation into collective mode coordinates shows that centre of mass mode b(†)cm=(b1(†)+b2(†))/ decouples from the electronic degrees of freedom and that only the mode corresponding to the relative displacement mode with bosonic operators

couples to the excitonic system. It is the non-classical properties of this collective mode that we investigate. In collective coordinates, the effective exciton–vibration Hamiltonian then reads

with σ=σ2+σ2−−σ1+σ1− and σ=(σ1+σ2−+σ2+σ1−). Tiwary et al.23 have recently pointed out that 2D spectroscopy can probe the involvement of these anticorrelated, relative displacement motions in electronic dynamics. From now on and for simplicity, we denominate this relative displacement mode as the collective mode.

We are interested in dimers that satisfy ΔE~ωvib>g>V where the effects of underdamped high-energy vibrational motions are expected to be most important193136. Several natural light-harvesting antennae include pairs of chromophores that clearly fall in this regime. Two important examples of such dimers are illustrated in Fig. 1a,b and correspond to the central PEB (phycoeritrobilin)50–PEB50 dimer in the cryptophyte antennae PE545 (ref. 33) and a Chl-Chl pair in the light-harvesting complex II (LHCII) complex of higher plants36; both corresponding to systems that have exhibited coherence beating in 2D spectroscopy3449. Importantly, in each case, the dimer considered contributes to an important energy transfer pathway towards exit sites333, suggesting that the the phenomena we discuss will have an effect in the performance of the whole complex. Moreover, synthetic versions of such prototype dimers could be available50. Most remarkably, LHCII is likely the most abundant light-harvesting complex on Earth35, while cryptophyte antennae such as PE545 are ecologically important as they support photosynthesis under extreme low-light conditions5152. From this perspective, the dimers of interest are exceptionally relevant biomolecular prototypes. Spectroscopy studies indicate that these dimers are subject to a structured exciton–phonon interaction as considered in our model. For the PEB50 dimer, the intramolecular mode of interest has frequency around 1,111 cm−1 (ref. 33), which compares with the frequency of the breathing mode of the tetrapyrrole53 (Carles Curutchet, personal communication). In the case of the Chl pair, it has been shown that a mode around 750 cm−1 is coupled to the electronic dynamics36 and this energy is close to the frequency of in-plane deformations of the pyrrole54. Furthermore, vibrational dephasing in chromophores55 and in other systems such as photoreceptors56 is known to be of the order of picoseconds. Some aspects of the influence of non-equilibrium vibrational motion in these specific dimers have been considered before1957; however, none of these studies have addressed the question of interest: can vibration-assisted transport exploit quantum phenomena that have no classical analogue?

Figure 1

Prototype dimers.

(a,b) Cryptophyte antennae phycoerythrin 545 (PE545) and LHCII present in higher plants have pairs of pigments whose electronic and vibrational parameters fall in the regime of our vibration-assisted transport model. (a) Representation of the pigments and protein environment of a PE545 complex of Rhodomonas CS24 (Protein Data Bank ID code 1XG0, ref. 32). The central PEB dimer pigments PEB50 and PEB50 are highlighted in red and green, respectively. For this PEB50 dimer, there is an uncertainty in the value of the energy gap3233. We take parameters from refs 32, 33 such that Δε=1,042 cm−1 and V=92 cm−1, so ΔE=1,058.2 cm−1 being quasi-resonant with an intramolecular mode of frequency ωvib=1,111 cm−1. The strength of linear coupling to this mode is g=ωvib(0.0578)1/2=267.1 cm−1. (b) Representation of the LHCII antennae of Spinacia oleracea (Protein Data Bank ID code 1RWT, ref. 34). Several pairs of close Chl-Chl (red-green) chlorophylls satisfy the conditions of our model. In particular, we consider the Chl-Chl pair for which Δε=661 cm−1 and V=−47.1 cm−1, resulting in ΔE=667.7 cm−1 (ref. 66). An intramolecular vibrational mode of frequency ωvib=742.0 cm−1 is close to this energy gap and each chromophore couples to this mode with strength g=ωvib(0.03942)1/2=147.3 cm−1 as obtained from (ref. 40). Scale bar: 1 nm.

Non-classicality via coherent exciton–vibration dynamics

We first consider the quantum coherent dynamics associated to Hex−vib to illustrate how non-classical behaviour of the collective motion emerges out of an initial thermal phonon distribution and an excitonic state with no initial superpositions: , which in the basis of exciton–vibration states of the form |X,n› (see Fig. 2a), becomes . Here n denotes the phonon occupation number of the relative displacement mode coupled to exciton dynamics (see Eq. 4), while Pth(n) denotes the thermal occupation of such level. The observables of interest are the population of the lowest excitonic state ρ(t)=∑〈Y,n|ρ(t)|Y,n›, the absolute value of the coherence ρ(t)=〈X,0|ρ(t)|Y,0›, which denotes the inter-exciton coherence in the ground state of the collective vibrational mode, and the non-classicality given by negative values Q(t) and corresponding negativities in the regularized quasi-probability distribution P(α). Hamiltonian evolution generates coherent transitions from states |X,n› to |Y,n+1› (see Fig. 2a) with a rate f that depends on the exciton delocalization (|V|/Δε), the coupling to the mode g and the phonon occupation n—that is, . Since ωvib≫KBT the ground state of the collective mode is largely populated, such that the Hamiltonian evolution of the initial state is dominated by the evolution of the state |X,0›. This implies that the energetically close exciton–vibration state |Y,1› becomes coherently populated at a rate , leading to the oscillatory pattern observed in the probability of occupation ρ(t) as illustrated in Fig. 2b. The low-frequency oscillations of the dynamics of ρ(t) cannot be assigned to the exciton or the vibrational degrees of freedom alone as expected from quantum coherent evolution of the exciton-plus-effective mode system. For instance, if the mode occupation is restricted to at most n=1, the period of the amplitude of ρ(t) is given approximately by the inverse of

Figure 2

Exciton-collective mode states and free exciton dynamics.

(a) The energy levels of the exciton-collective mode states used to describe energy transfer . The red arrow denotes population transfer from to . (b) Quantum coherent dynamics of the PEB50 dimer in PE545 illustrating population of the lowest exciton ρ(t) (thick blue curve) and the inter-exciton coherence in the ground-state of the collective mode |ρ(t)| (thin red curve).

with α2=2g2+ and (2g24V2/α2ΔE2)<<1. Coherent exciton population transfer is accompanied by beating of the inter-exciton coherence |ρ(t)|, with the main amplitude modulated by the same low-frequency oscillations of ρ(t) and a superimposed fast oscillatory component of frequency close to ωvib (see Fig. 2b). This fast-driving component arises from local oscillatory displacements: when V≃0 the time evolution of each local mode is determined by the displacement operator with amplitude (ref. 58). As the state |Y,1› is coherently populated, the collective quantized mode is driven out of equilibrium towards a non-classical state in which selective occupation of the first vibrational level takes place, thereby modulating occupation of higher levels. This manifests itself in sub-Poissonian phonon statistics as indicated by negative values of Q(t) shown in Fig. 3a. Similar phenomena have been described in the context of electron transport in a nanoelectromechanical system5960. Moreover, Fig. 3b shows that at times when Q(t) is negative—that is, t=0.2 ps—the regularized quasi-probability distribution P(α) at this time exhibits negativities, thereby ruling out any classical description of the same phenomena. Interestingly, the non-classical properties of the collective vibrational motion resemble non-classicality of bosonic thermal states (completely incoherent states) that are excited by a single quanta3961. Importantly, such non-classical behaviour of the vibrational motion arises only when the electronic interaction between pigments is finite. For comparison, Fig. 3c shows that if V=0, an electronic excitation drives the local underdamped vibration towards a thermal displaced state with super-Poissonian statistics (Q(t)>0), which has an associated positive probability distribution in phase space as illustrated in Fig. 3d. In short, non-classicality of the collective mode quasi-resonant with the excitonic transition arises through the transient formation of exciton–vibration states.

Figure 3

Non-classicality of collective and local vibrational modes.

(a) Mandel Q-parameter of the relative displacement mode of the considered PEB50 dimer from PE545 when a biological electronic coupling is considered. Shaded regions denote times of non-classicality. (b) The associated regularized quasi-probability distribution P(α) at t=0.2 ps. Shaded regions denote areas of negative probability. (c) Mandel Q-parameter of the intramolecular high-energy vibration when dipole coupling is zero and (d) the associated regularized quasi-probability distributions P(α) at t=0.2 ps.

Energy and coherence dynamics under thermal relaxation

We now investigate the dynamics of the exciton–vibration dimer when each local electronic excitation interacts additionally with a low-energy thermal bath described by a continuous distribution of harmonic modes. The strength of this interaction is described by a Drude spectral density with associated reorganization energy λ and cutoff frequency Ωc2g|V|/Δ∈ population transfer to ρ(t) is incoherent. For the PEB50 dimer λ~110 cm−1 that place this dimer in this incoherent regime where ρ(t) has a slow but continuous exponential rise reflecting the fact that thermal fluctuations inducing transitions from |X,n› to |Y,n› now have a large contribution to exciton transport. However, even in this regime, transfer to ρ(t) is always more efficient with the quasi-resonant mode than in the situations where only thermal-bath-induced transitions are considered (see dashed lines in Fig. 4a,d,g). The underlying reason is that before vibrational relaxation takes place (around t=1 ps), the system is transiently evolving towards a thermal configuration of exciton-collective mode states. Hence, in both coherent and incoherent population transfer regimes, transfer to the lowest exciton state involves a transient, selective population of the first vibrational level of the collective mode. The transition from coherent to incoherent exciton population dynamics is then marked by the onset of energy dissipation of the exciton–vibration system as shown in Fig. 4c,f,i, where E(t)=Tr{Hex−vibρ(t)} has been depicted for different values of λ. While exciton population growth is non-exponential, energy dissipation into the thermal bath is transiently prevented as indicated by periods of positive slope of E(t) as happens in Fig. 4c,f. Quantification of the energy that is transiently ‘extracted’ from the low-energy thermal bath can provide an interesting physical interpretation of the advantages of non-exponential exciton transfer in the framework of non-equilibrium thermodynamics62.

Figure 4

Energy and coherence evolution under thermal relaxation.

The dynamics of ρ(t) with (thick blue curve) and without (dashed curve) coupling to vibration (top row), |ρ(t)| (middle row) and energy of the exciton vibration system E(t)=Tr{Hex−vibρ(t)} (bottom row) for the exciton–vibration parameters of the PEB50 dimer and three interaction strengths to the low-energy thermal bath: (a–c): λ=6 cm−1, (d–f) λ=35 cm−1 and (g–i): λ=110 cm−1. Initial energy E(0), displayed as dashed line and times where E(t)>E(0) are shaded.

For completeness, we present in Figure 4b,e,h how the beating patterns of the coherence ρ(t) reveal the structured nature of the exciton–phonon interaction and witnesses whether there is coherent exciton–vibration evolution as it has been pointed out by recent studies2122. The frequency components of such oscillatory exciton coherences vary depending on the coupling to the thermal bath. In the coherent regime, as ρ(t) follows exciton populations, the main amplitude is modulated by the same relevant energy difference between exciton–vibration states (see Fig. 4b,e). This behaviour is relevant for the parameter regime of the Chl dimer (see Supplementary Fig. S1). In contrast, for the PEB50 dimer, the short-time oscillations of ρ(t) (between t=0 and t=0.1 ps) arise from purely electronic correlations because of bath-induced renormalization of the electronic Hamiltonian63. This exciton coherence retains the superimposed driving at a frequency ωvib and is accompanied by non-classicality as it will be described shortly, indicating that vibrational motion is still out of thermal equilibrium. The dynamical features presented in Fig. 4g,h agree with previous findings based on a perturbative approach19 and with the timescales of the exciton coherence beating reported for cryptophyte algae418.

Non-classicality under thermal relaxation

Interaction with the thermal environment would eventually lead to the emergence of classicality in longer timescales. However, in the picosecond timescale of interest, the collective mode exhibits periods of non-classicality across a wide range of thermal bath couplings λ as indicated by sub-Possonian fluctuations with Q(t)<0 in Fig. 5a–c and the corresponding negativities in the distributions P(α) shown in Fig. 5d–f. This survival of non-classicality is concomitant with a slow decay of the exciton–vibration coherence ρ(t) (not shown). Non-classical behaviour of collective fluctuations are then expected for the parameters of both the PEB50 dimer for which λ=110 cm−1 and the Chl dimer for which λ=37 cm−1. The non-classical fluctuations predicted by Q(t) also agree with those witnessed by a parameter quantifying correlations between nearest-neighbours’ occupations48, which we present in Supplementary Fig. S2. As expected, the maximum non-classicality indicated by the most negative value of Q(t) decreases for larger reorganization energies. Nonetheless, the time average of these non-classicalities is not a monotonic function of λ. For moderate values of λ, the collective mode spends longer periods in states with non-classical fluctuations—that is, periods for which Q(t)<0 as seen in Fig. 5b, thereby stabilizing non-classicality at a particular level. This sub-picosecond stabilization of non-classicality is expected in the regime of the Chl dimers as illustrated in Supplementary Fig. S1.

Figure 5

Non-classicality of collective motions under thermal relaxation.

(a–c) Dynamics of the Mandel Q-parameter for λ=6, 35, and 110 cm−1, respectively. Shaded regions denote times of non-classicality. (d–f) Regularized distribution P(α) at t=0.2 ps for each corresponding value of λ. Shaded regions denote areas of negative probability.

Functional role of non-classicality

Non-classical fluctuations of collective motions correlate to exciton population transfer. In order to demonstrate this, we investigate quantitative relations between non-classicality and exciton energy transport by considering relevant integrated averages in the timescale of the Hamiltonian evolution of the exciton–vibration system denoted by τ. For the parameters of the PEB50, this timescale is about half a picosecond and is comparable to the timescale in which excitation energy would be distributed away to other chromophores or to a trapping state. The time-integrated averages over τ are defined as: 〈F[ρ(t)]›=(1/τ)∫0dt F[ρ(t)], where F[ρ(t)] corresponds to the exciton population ρ(t) and the non-classicality of the underdamped collective mode through periods of sub-Poissonian statistics Q(t)Θ[−Q(t)] as functions of the coupling to the bath λ. As shown in Fig. 6, the average exciton population and non-classicality follow a similar non-monotonic trend as a function of the coupling to the thermal bath, indicating a direct quantitative relation between efficient energy transfer in the timescale τ and the degree of non-classicality. The appearance of a maximal point in the average non-classicality as a function of the system-bath coupling indicates that the average quantum response of the collective anticorrelated motion to the impulsive electronic excitation is optimal for a small amount of thermal noise.

Figure 6

Correlations between non-classicality and exciton populations.

Time-integrated averages of exciton population ρ(t) and non-classicality as quantified by Q(t)Θ[−Q(t)] (blue and green, respectively) as functions of coupling to the thermal background by fixing environment cutoff frequency Ωc=100 cm−1 and varying reorganization energy λ with exciton–vibration parameters corresponding to the PEB50 dimer.

It is worth emphasizing that the above functional role of non-classicality holds for vibration-assisted transport where the high-energy intramolecular modes considered are quasi-resonant with the excitonic energy splitting. When vibrational motions are significantly detuned with the bare exciton transition ΔE, transport is dominated by the thermal background and no selective population of the state |Y,1› takes place; hence, periods of sub-Poissonian fluctuations vanish. To illustrate the difference with the off-resonance case, Fig. 7 shows the same time-integrated averages as in Fig. 6, with the electronic parameters of the PEB50 dimer but now considering an intramolecular vibration of frequency ωvib=1,520 cm−1 significantly detuned from ΔE. In this case, thermally activated transport (see dashed line in Fig. 7) and vibration-assisted transport (solid line in Fig. 7) are practically indistinguishable and the average Q(t) is positive with a value near zero as expected for a thermal distribution of a high-energy harmonic mode.

Figure 7

Effects of off-resonance vibrations.

Time-integrated averages as in Fig. 6 with electronic parameters of the PEB50 dimer but now considering an intramolecular vibration which is ωvib=1,520 cm−1—that is, significantly detuned from the energy gap ΔE=1,058.2 cm−1 and similar coupling strength g=ωvib(0.0265)1/2≈247 cm−1. Parameters obtained from ref. 33. The dashed curve shows time average of ρ(t) when no intramolecular vibration is considered and transport results only because of the thermal background.

The degree of purity of the initial exciton state is also important in enabling and harnessing non-classical fluctuations of the collective mode. One therefore should expect that statistical mixtures of excitons with finite purity can still trigger such non-classical response. To illustrate this, we now consider mixed initial states of the form ρ(t0)=ρex⊗ρthvib where with 1/2≤r≤1. The associated linear entropy quantifying the mixedness of the initial exciton states is given by SL=2r(r−1). The time-averaged non-classicality (Fig. 8a) and average population transfer (Fig. 8b) follow similar decreasing, yet non-zero, monotonic trends for mixed states. These results suggest that non-classical vibrational motion can be prompted and exploited even under incoherent conditions creating statistical mixture of excitons64. The trends presented in Figs 6, 8 constitute theoretical evidence of direct quantitative correlations between non-classicality and exciton population transport in a relevant sub-picosecond timescale and therefore illustrate a functional role for quantum phenomena with no classical counterpart in a prototype light-harvesting system. Our results remain valid when the picosecond-dephasing rate of the vibrational motion is included in the dynamics as can be seen in Supplementary Fig. S3.

Figure 8

Non-classicality under incoherent exciton input.

The time-integrated averages of exciton population ρ(t)−ρ(0) and non-classicality as quantified by Q(t)Θ[−Q(t)] as functions of the linear entropy SL=2r(r−1) for 0

Discussion

Light-harvesting complexes are fundamental components of the photosynthetic machinery on Earth. While there has been an extraordinary advance in the techniques to probe these systems at ultrafast timescales, it is conceptually unclear what truly quantum features with no classical counterpart these systems may exhibit and exploit for efficient energy transport65. Besides depending on long-range electronic interactions, excitation energy distribution is fundamentally affected by the molecular motions and harmonic (or in cases anharmonic) environments that modulate the electronic dynamics. We therefore put forward the idea that, precisely, investigation of non-classical phenomena associated to such molecular motions can pave the way towards understanding which non-trivial quantum phenomena can have an impact on efficient energy distribution and trapping. Within this line of thought, we investigate non-classical behaviour in a prototype dimer that conveys fundamental physical principles of vibration-assisted transport in a variety of light-harvesting antennae. We demonstrate that quasi-resonances between excitonic transitions and underdamped high-energy intramolecular vibrations can trigger and harness non-trivial quantum behaviour of collective pigment motions that are initially in a thermal, fully incoherent, state. In this scenario, correlations between non-classical fluctuations of collective pigment motions and efficient population of the target exciton state are found. Negative values both of the Mandel Q-parameter and of the quasi-probability distribution P(α) of the collective motion assure that no classical distribution can describe the same behaviour. These non-trivial quantum phenomena are predicted for a variety of initial excitation conditions including statistical mixtures of excitons indicating that such non-classicality can be activated even under incoherent input of photoexcitations. Transient coherent ultrafast phonon spectroscopy29, which is sensitive to low-phonon populations of high-energy vibrations30, may provide an interesting experimental approach to probe the phenomena we describe.

The prototype exciton–vibration dimers here investigated are representative of interband-like transfer pathways present in the majority of light-harvesting complexes. For instance, in the LHCII complex, the considered dimer contributes to the fastest component of the Chl→Chl interband transfer pathway that directs excitation energy towards exit sites66. The demonstration that non-classicality is concomitant with efficient vibration-assisted transfer in these dimers therefore suggests that non-classical phenomena will have a contribution to the efficiency of the whole complex. A rigorous quantitative estimation of such contribution requires both a careful extension of our formalism to quantify these features in a multi-modal system (as likely other non-equilibrated vibrational motion will be involved) and a careful weighting of the vibration-assisted processes in the overall spatio-temporal distribution of energy.

The framework we propose can also be applied to gain insights into the non-classical response of vibrational motion in a variety of transport4041 and sensing processes in biomolecules4243. Of particular interest are charge transfer in reaction centres40 and isomerization of photoreceptors42 where specific intramolecular vibrational motions are known to be driven out of thermal equilibrium during the light-initiated electronic dynamics. It will also be interesting to use this framework to understand possible non-trivial quantum behaviour of molecular motions in chemical sensors43 that are conjectured to operate through weak electronic interactions to sense molecular vibrations of the order of a thousand wavenumbers4445.

We have also illustrated how in our prototype dimer with biologically relevant parameters, exciton–vibration dynamics can lead to non-exponential excitonic energy distribution whereby dissipation into a low-energy thermal bath can be transiently prevented. From this view, coherent vibrational motions that do not relax quickly and whose fluctuations cannot be described classically may be seen as an internal quantum mechanism controlling energy distribution and storage. Further insights into the advantage of these non-trivial quantum behaviour may therefore be gained in a thermodynamic framework6267.

In conclusion, we have provided theoretical evidence that vibration-assisted exciton transport in prototype dimers, representative of interband-like transitions in a variety of photosythetic light-harvesting antennae, can exploit non-trivial quantum phenomena that cannot be reproduced by any classical counterpart, namely, non-classical fluctuations of collective pigment motions. Given that a variety of transport4041 and sensing phenomena42 in biomolecules are known to involve non-equilibrium vibrational motion, our findings have broad implications for the field of quantum effects in biology as they suggest that investigating the non-classical nature of molecular fluctuations harnessed in these processes could be the key to reveal a role for truly non-trivial quantum features.

Methods

Open quantum system dynamics

To accurately account for the effects of the low-energy thermal bath, we have adopted a hierarchical expansion of the exciton–vibration dynamics686970. We split the high-energy mode from the bath of harmonic oscillators and explicitly treat the quantum interactions between electronic excitations and these modes of frequency ωvib by including it within the definition of the system, Hex−vib=Hex⊗11vib+11ex⊗Hvib+Hex−vib. The electronic operators then couple to the remaining vibrational modes HI=∑g(σ+σ−⊗11vib)(b†+b). This approach allows the effects of the low-energy thermal bath on the exciton–vibration dynamics to be accurately accounted for. Truncating the quantized mode at Fock level M=6 adequately describes both reduced dynamics of the collective quantized mode and the electronic dynamics of the prototype dimer at room temperature. A spectral density of the form J(ω)=2λΩcω/(Ωc2+ω2) is assumed, where λ and Ωc are the reorganization energy and cutoff frequency, respectively. Supplementary Note 1 furnishes further details of the hierarchical expansion of exciton–vibration dynamics. Converged dynamics are obtained by terminating the hierarchical expansion at level N=11 and including just the K=0 Matsubara term. No additional Matsubara terms are necessary as Ωc

Regularized quasi-probability distributions

The quasi-probability distribution we calculate is a regularized version of the P-representation P(α)=(1/π2)∫d2ξ e*χ(ξ)Ω(ξ), where the quantum characteristic function χ(ξ)=Tr{erdrdρvib} of the reduced state the vibration ρvib is reconstructed in the truncated Fock basis . Ω(ξ) is a non-classicality filter38 taken as the triangular function with w=5, which fulfils the necessary condition of a filter such that P(ξ) detects non-classicality. More details are given in Supplementary Note 2.

Author contributions

E.J.O. performed the calculations. A.O.-C. conceived, designed and supervised the research. Both authors analysed and discussed the results and co-wrote the manuscript.

Additional information

How to cite this article: O’Reilly, E. J. and Olaya-Castro, A. Non-classicality of the molecular vibrations assisting exciton energy transfer at room temperature. Nat. Commun. 5:3012 doi: 10.1038/ncomms4012 (2014).