photosynthesis , arguably the most important photo - induced process on earth , converts the energy of light into its chemically / biologically useful form .
it is often argued that this conversion is remarkably efficient .
however , it has to be distinguished between the _ quantum _ efficiency , reaching almost unity@xcite , and light - to - chemical energy efficiency , which is significantly lower , in the order of several percent@xcite .
this is given by an evolutionary pressure on the development of a robust photosynthetic machinery optimized for survival rather than energy conversion efficiency .
energy relaxation processes are integral part of the photosynthetic function as they enable the energy transfer to proceed unidirectionally @xcite , and provide protection of the photosynthetic apparatus against harmful over - excitation .
the photosynthetic machinery of plants has developed a complex hierarchy of self - regulatory mechanisms to avoid excess excitation or ( when unavoidable ) to dissipate it into heat @xcite .
starting from processes controlled on the macroscopic level ( e. g. orientation of leaves ) , over spontaneous microscopic ( cellular ) events such as chloroplast movements , to truly nano- and sub - nanoscopic mechanisms such as reorganization of antenna complexes and direct regulation of energy transfer on the level of small groups of interacting chromophores , plants actively react to changing illumination conditions .
the sub - nanoscopic processes , which are the focus of the present study , operate in response to the increase of @xmath0 across the thylakoid membrane .
such an increase is an indicator of high illumination .
nowadays it is generally accepted that carotenoids are involved in these energy dissipation processes .
the precise molecular mechanism is , however , still subject of discussion@xcite .
it is likely that different mechanisms evolved in different classes of organisms and/or that several mechanisms operate at once .
most of our current knowledge about the early processes of photosynthesis was obtained by ultrafast spectroscopy . while conventional bulk spectroscopies are extremely useful in following ultrafast photo - induced events in photosynthesis , whenever structural inhomogeneity of the sample is involved , the information obtained from these spectroscopies becomes obscured by an inevitable ensemble averaging .
some well established spectroscopic methods , such as hole - burning @xcite , and some modern multidimensional methods , such as coherent two - dimensional spectroscopy @xcite , provide certain access to the homogeneous properties .
information on the behaviour of individual molecules has been , however , available only since the advent of single - molecule spectroscopy ( sms ) @xcite .
the system studied in this work , the light harvesting complex ii ( lhcii ) of higher plants , is the plants major light - harvesting antenna containing almost half of all the chlorophyll in the chloroplast .
correspondingly , most of the light absorption and subsequent energy transfer processes in plants and algae occur in this complex .
the lhcii antenna occurs naturally in a trimeric form and its main function is to deliver excitation energy to the nearby photosystem ii ( psii ) .
given the major role of lhcii in light - harvesting and energy transfer , it is not surprising that it is also implicated in participating in regulated energy dissipation , the so - called non - photochemical quenching ( npq ) @xcite .
the crystal structure of the complex @xcite enables us efficient theoretical modeling of the complex s spectroscopic properties using its chromophores ( chlorophylls ( chl ) and carotenoids ) as the model units .
the model parameters such as electronic coupling are greatly constrained by the known mutual orientation and distances .
existence of a large body of previous measurements together with theoretical attempts to fit this whole body of data to a single model @xcite gives a great confidence in extending the modeling towards single molecular experiments .    in recent years , single - molecule spectroscopy ( sms ) experiments on several photosynthetic antenna complexes including the lhcii were performed .
fluorescence spectral peak distributions @xcite , spectral diffusion @xcite , fluorescence intensity distributions @xcite and time traces @xcite were obtained from these measurements . in many cases also
fluorescence intermittency ( blinking ) was observed , and it was conjectured that in the case of lhcii the process behind the fluorescence intermittency plays a role in npq @xcite .
as far as the theory of the sms of lhcii is concerned , the ensemble - averaged spectra and also the peak distributions of lhcii trimers can be successfully explained by the disordered excitonic model @xcite . in general , the excitonic model was successfully applied in the past on fs to ps timescale , and it represents an indispensable tool in analyzing ultrafast spectroscopic experiments on molecular aggregates and in particular on photosynthetic antennae@xcite .
its application to longer time scales of seconds and minutes is conditioned by the assumption of a certain separation of time - scales . over the course of the excitation - emission cycle ( nanoseconds ) ,
individual chromophores of a complex are assumed to be found in fixed spatial arrangements and experiencing an environment described by a fixed set of parameters .
the emission spectrum of a complex in such a fixed spatial arrangement is predicted by the excitonic model which gives the population distribution of the excited state manifold and the probability of emission at the corresponding transition frequencies .
the spectrum of exciting light matters only to the extent to which the excited states reached at a given excitation wavelength are connected to the final state by some relaxation pathways .
once the pathways are available , the final state is given irrespective of the initial state after absorption of light . because only the final distribution of excited state population matters , the changes ( switches ) of the spatial arrangement or environmental conditions occurring on the sub - nanosecond time scale are only observed as sudden changes ( with respect to the nanosecond fluorescence time scale ) of the fluorescence spectra . despite the fact that the experiment we analyze in this work does not provide more insight into the sub - nanosecond dynamics of individual complexes than previous works
, we nevertheless formulate the theory consistently in such a way that it enables the description of such processes .
this is done in order to highlight the existence of a less studied type of processes on an intermediate time scale and to stress the need to search for experimental techniques which can cover the range of timescales from femtoseconds to nanoseconds in single molecular spectroscopy .
we thus provide theoretical techniques to treat these future experiments .
we have recently reported a progress towards measuring processes on the order of hundreds of femtoseconds in single light - harvesting complexes @xcite .
the paper is organized as follows : in section [ sec : excited - state - population ] we first discuss the frenkel exciton model as a basis for the formulation of equations of motion for the populations of the excited states of chromophore aggregates with strong interchromophore couplings .
then , in section [ sub : excited - state - dynamics ] , we introduce equations of motion for excitonic populations valid over a broad ( from ps to min ) timescale range .
we discuss their generality and the range of validity .
it is argued that these equations provide an ideal means for the description of the sms experiments . in section [ sec : lhcii ] the results of section [ sub : excited - state - dynamics ] are applied to lhcii photosynthetic antenna complexes .
all results are compared to the experiment .
it is shown that our equations give correct fluorescence spectrum and peak statistics , i. e. appropriate steady - state population , under typical sms experiment condition .
then we model the intensity traces , while the switching behaviour is included by incorporating one particular previously proposed npq mechanism , namely energy transfer to lutein lut 1 .
the switching between on and off states is controlled by a 2-level model , where the switching causes a change of the chl a612 - lut 1 coupling .
it is shown that using realistic parameters we are able to reproduce the experimental results .
the details of our energy relaxation theory and of the stochastic model of switching are given in supporting information ( si ) .
in the usual sms experiments , time - resolved ( on the times scale of 10 ms to s ) fluorescence of the studied molecules is observed .
if the excited state life time of the studied chromophores is sufficiently long ( nanoseconds in the case of chlorophylls studied here ) , the expected fluorescence spectra can be calculated from the steady state , quasi - equilibriated populations of the excited electronic eigenstates of the molecular system , assuming canonical equilibrium .
depending on the strengths of the chromophore - chromophore resonance interaction and the system - bath couplings , the eigenstates can be approximated either by the excited eigenstates of the electronic subsystem or the excited state of the individual chromophores forming the aggregate . in the present paper , we treat both cases within one formalism provided by the framework of the frenkel exciton model .
the frenkel exciton model provides an excellent tool for the treatment of pigment - protein complexes on femto- and pico - second time scales @xcite .
the basic notion of the frenkel exciton model is the one of localized excited states of the chromophores .
these states have negligible wavefunction overlap with neighboring excited states , and they can therefore be assumed orthogonal to each other , forming a suitable basis for the aggregate hilbert space . in the treatment of linear absorption and fluorescence experiments ,
only singly excited collective states of the molecular aggregate have to be included .
the localized excited states can thus be denoted as @xmath1where @xmath2 and @xmath3 are the electronic ground and excited states of the chromophore denoted by index @xmath4 , respectively .
the states @xmath5 form a complete hilbert subspace for the case that exactly one molecule of the aggregate is excited .
the molecular system hamiltonian , @xmath6 , is however rarely diagonal in the basis of the states @xmath5 and resonance couplings @xmath7 between excited electronic states @xmath5 and @xmath8 occur in all situations interesting for light - harvesting . in the absence of a protein environment ,
the light would resonantly excite eigenstates of the system hamiltonian . in all realistic cases ,
the interaction of the system with its environment co - determines the excitation frequencies , and prescribes thus a `` preferred basis '' of electronic states in which it is the most advantageous to formulate the energy transfer theory .    in the case of the resonance interaction @xmath7 exceeding the typical reorganization energy @xmath9 of the protein bath of the chromophores @xmath4 and @xmath10 , ( @xmath11 , we include @xmath7 explicitly into the system hamiltonian .
the hamiltonian is diagonalized to obtain electronic eigenstates , and the effect of the protein bath is included via perturbation theory yielding a redfield type relaxation tensor . in the opposite case
@xmath12 we neglect the resonance coupling in the system hamiltonian , and its effects are included perturbatively yielding frster type energy transfer rates between localized excited states @xcite .
let us first focus on deriving a closed set of equations for the excited state populations .
extensive work was done in the last years on developing methods to accurately describe the system dynamics following an ultrafast excitation by external light . the traditional redfield and
frster approaches were superseded by more accurate ( exact in some cases ) methods , such as heom @xcite , tedopa @xcite and other methods @xcite .
these methods bring unprecedented accuracy at an increased numerical cost .
it seems , however , that for analyzing many of the recent experiments , it is still possible to rely on the traditional tools , as they capture the physics of the problem ( and often even the quantitative aspects of the problem ) very well @xcite .
commonly , equations of motion for some relevant degrees of freedom ( dof ) , electronic states in our case , are derived by reducing the liouville - von neumann equation for the total density matrix to an equation for the so - called reduced density matrix ( rdm ) .
these equations describe time evolution of a molecular system for a fixed configuration of the protein environment ( assuming fast fluctuations around this fixed configuration ) , and they are therefore suitable for the description of ultrafast laser experiments .
however , on the timescale relevant in sms ( up to tens of seconds ) , usage of these equations is actually not appropriate .
first of all , some transient effect at the short times affect even the long time properties of the system and the steady state population dynamics , and second , slow changes in the protein environment can entirely change the energy landscape including the case that one has to change the theoretical limit ( localized states , delocalized states ) in which one works .
the latter case is especially difficult to treat and requires to go beyond the traditional master equation approaches which we apply here . in this work we therefore concentrate on the extension of the validity of the master equation approach towards long times under the assumption that the dependence of the hamiltonian of the system on time is negligible within one absorption emission time scale ( nanoseconds ) .
as for the transient effects at short time , when dealing with fast dephasing of optical coherences , short time transient effects are responsible for the absorption lineshape . this aspect of the transient effects will be taken into account in full .
spontaneous emission of photons by chlorophylls occurs with a nanosecond life - time .
another class of transient effects , namely dynamic electronic coherence due to excitation by light , is therefore also unimportant and its treatment can be avoided .
we therefore derive approximative equations for the populations dynamics only , with the validity in the range from picoseconds to tens of seconds .
initially , our hamiltonian consists of the system , bath and system - bath interaction terms , and we formally assume some total density matrix @xmath13 which follows the liouville  von neumann equation .
there are several methods how to arrive to a master equation ( in a convolution - less form ) rdm @xmath14 @xcite . in general ,
such equation reads :    @xmath15-{\cal r}(t , t_{0})\rho(t),\label{eq : equation_of_motion}\ ] ]    @xmath16    the total hamiltonian @xmath17 consists of the bath renormalized system hamiltonian @xmath18 and the system - light interaction , which is ( in the dipole approximation ) given by the dipole moment operator @xmath19 and electric field @xmath20 .
the bath is completely eliminated in the reduced description , and its effects are represented by the relaxation tensor @xmath21 .
it is important to note that in eq .
( [ eq : equation_of_motion ] ) , the relaxation tensor @xmath21 explicitly depends on some initial time @xmath22 , in which the initial condition for the propagation is set , or more precisely , the relaxation tensor depends on the quantity @xmath23 .
an exact relaxation tensor would also depend on @xmath20 .
in general , the choice of @xmath22 is arbitrary , and we should choose it such that the temporal profile of the excitation field is non - zero for @xmath24 only .
however , in most practical theoretical approaches , the evolution of the bath due to excitation of the system is reflected in the relaxation tensor @xmath21 by the bath correlation function , which decays rather quickly . by choosing @xmath22 sufficiently small
, we could always make such a relaxation tensor constant .
this is obviously an artifact of the approximations used . for a smooth envelope of the external field changing on the same or slower timescale than the correlation function of the bath
, there is no good choice of @xmath22 .
the relaxation tensor always becomes constant too early .
luckily , as we will see below , we do not need to account for the external field @xmath20 beyond first order perturbation theory , and in this regime , there is a natural choice of time @xmath22 which enables us to correctly account for the transient time dependence of the relaxation tensor even for a steady state externally driven by light .
it is important to note , that the relaxation tensor @xmath21 is completely abstract up to now .
it can represent some exact relaxation superoperator , or a result of perturbation theory with respect to some parameters , such as the redfield or frster tensors .
let us now formulate the equations of motion for the excited state evolution in the linear regime of the system s interaction with the electric field .
the validity of the linear regime has been discussed e.g. in refs .
@xcite , and it is the same as the validity of the third order response theory for non - linear laser spectroscopy .
we will also use the secular approximation ( equations of motion for populations @xmath25 and coherences @xmath26 , @xmath27 are assumed independent ) , although this we do only to simplify the resulting equations .
secular approximation could be avoided if one so wishes , at the cost of treating the full dm .
we write explicitly the elements of the rdm of eq .
( [ eq : equation_of_motion ] ) in which we keep the double time dependence of the hamiltonian and the relaxation tensor .
we get the following set of coupled equations : @xmath28 @xmath29 here , the @xmath30 are transfer rates between populations ( from @xmath10 to @xmath4 ) , @xmath31 is the population relaxation rate of the state @xmath4 , including the rate of radiative depopulation @xmath32 ( we assume it is constant ) , and @xmath33 is the optical coherence dephasing rate .    the equation of motion for the optical coherences , eq .
( [ eq : coherences eq ] ) , can be solved by introducing the interaction picture using the evolution operator element @xmath34 , and integrating eq .
( [ eq : coherences eq ] ) in the interaction picture . by returning back to the schrdinger picture
, we obtain the actual field - induced and field driven optical coherence element in the form : @xmath35e(t - t^{\prime}).\label{eq : coherence_optical}\ ] ] we can see that eq .
( [ eq : coherence_optical ] ) is actually not a solution of eq .
( [ eq : coherences eq ] ) . in solving eq .
( [ eq : coherences eq ] ) numerically , the dephasing rate would quickly become constant for @xmath36 and the transient effects for @xmath37 would be lost completely .
however , in eq .
( [ eq : coherence_optical ] ) these effects are properly taken into account .
the discrepancy between eq .
( [ eq : coherences eq ] ) and eq .
( [ eq : coherence_optical ] ) is due to a different treatment of the bath . in a response function approach ,
of which the linearization of the full equations of motion with respect to the field is a variant , it is in general possible to account for the bath in a more consistent way than in master equations . in ordinary master equations ,
the bath is correctly described at @xmath22 , and the description of its subsequent evolution after @xmath22 is extremely limited .
an example of such limitation is discussed e.g. in ref .
@xcite .    in most biological energy transferring systems , the pure dephasing is much faster ( hundreds of fs ) than the population relaxation ( units and tens of ps ) .
we can therefore assume that the populations remain constant during the integration in eq .
( [ eq : coherence_optical ] ) .
moreover , it is reasonable to set the upper limit of the integration in eq .
( [ eq : coherence_optical ] ) to infinity by sending @xmath38 .
we can expect that the simultaneous action of pure dephasing and external driving by a field with a slowly varying envelop creates a steady state optical coherence .
the dependence on @xmath22 in the upper limit of the integral describes a transient evolution of the optical coherence after switching on the interaction with the field .
now that the short time time - dependent nature of the dephasing rates is taken into account correctly , we can set @xmath38 and write :    @xmath39.\label{eq : coherence ft}\ ] ]    it is important to note that now @xmath40 is a global time which can run through the whole minutes long sms experiment .
the evolution operator changes on an ultrafast time scale , but this timescale is scanned in the integration over the variable @xmath41 .
the properties of the hamiltonian @xmath42 change on a very slow time scale ( with respect to optical dephasing ) , and so does the evolution operator element @xmath43 .
unlike the eq .
( [ eq : coherences eq ] ) , which is valid for @xmath23 small , eq . ( [ eq : coherence ft ] ) is valid for all times .
the purpose of deriving eq .
( [ eq : coherence ft ] ) was to insert it eventually into eq .
( [ eq : populations eq ] ) . also here
we face the problem of transient time dependence of the rates .
however , because populations change much more slowly , these effects are not as important as in the case of coherences .
they can be , however , treated rigorously , even including bath memory effects between evolutions by eq .
( [ eq : coherences eq ] ) and eq .
( [ eq : populations eq ] ) , as we have shown elsewhere @xcite . in this work , we will assume the energy transfer and relaxation rates not to depend on the difference @xmath23 , although they may depend weakly on the time @xmath40 , i. e. @xmath44 and @xmath45 .
we are ready now to insert eq .
( [ eq : coherence ft ] ) for optical coherences into the equation for populations , eq .
( [ eq : populations eq ] ) .
we obtain    @xmath46    here , we defined @xmath47 , and we used the fact that @xmath48 and even @xmath34 are ordinary c - numbers .    in eq .
( [ eq : driven_pop_evol ] ) , the populations are driven by a second order field term .
we have treated the field classically so far .
if we did that quantum mechanically , we would now have to trace over the field dof in order to obtain reduced equations of motion for the electronic state populations only .
the term @xmath49 would thus be replaced by @xmath50 which can be interpreted as a quantum mechanical expectation value .
the latter expression has the form of two - time correlation function of the electric fields of the light and its fourier transform is the power spectrum of the light @xcite ,    @xmath51    now , inserting eq .
( [ eq : power_spectrum ] ) into eq .
( [ eq : driven_pop_evol ] ) , switching the order of integrations , and using the definition of absorption lineshape of the i - th excitonic transitions ( see e.g. @xcite )    @xmath52    we arrive at    @xmath53\int d\omega w(\omega)\chi_{i}(\omega;t).\label{eq : popevolution}\ ] ]    our effort has yielded a closed set of equations for excitonic populations only .
the population changes are given by the transfer rates between electronic levels , population quenching and a source terms expressed as an overlap of the excitonic spectra with the light spectrum .
all quantities are in principle dependent on time , most importantly the excitonic absorption spectrum and all rates can weakly depend on time to simulate slow changes of the protein and chromophore configurations . also the light spectrum can be considered time - dependent by generalizing the wiener - khintchine theorem for the instaneous power spectrum @xmath54@xcite .
the changes can be faster than the time resolution of the sms experiment , but they have to be slower than the dephasing , or even energy transfer dynamics .
we should also note here that if we wish to stay strictly in the linear regime , the saturation term @xmath55 in eq .
( [ eq : popevolution ] ) should be very close to one .
in fact , we can set it equal to one as a reasonable approximation with the same validity as the linear approximation for fields .        according to crystalographic studies , lhcii complex consists of three monomeric units , each containing 14 chlorophylls and for carotenoids : two luteins , neoxanthin and a carotenoid of the xanthophyll cycle@xcite .
the experiments described in this section were performed on monomeric lhcii complexes . in accord with the experiment , we focus on one such a monomeric unit in our model .
we treat the lhcii monomer as a strongly coupled systems of chromophores , weakly coupled to the bath and weakly interacting with light . because only absorption and fluorescence are measured , we do not attempt to fit the site energies in our simulations from scratch , as the fitting of this limited set of experiments would not be unique . instead , we take the pigment transition energies from ref .
@xcite , where both lhcii trimers and monomers where treated .
the coupling energies between the pigments were calculated in the dipole - dipole approximation , and the dipole orientations were taken from the crystal structure using an effective dipole strengths of 3.4 d for chl b and 4.0 d for chl a. the bath is described by means of a spectral density obtained by fluorescence line narrowing experiment ( fln ) ( see ref .
excitonic absorption and fluorescence lineshapes are calculated by means of the second order cumulant expansion ( see ref .
@xcite ) , and the population transfer rates are calculated by the multilevel redfield theory ( see ref .
@xcite ) . for comparison
, we also calculated the rates by modified redfield theory , ref.@xcite , and we concluded that the results remain essentially the same .
the population relaxation rates of chlorophylls were taken to be 3 ns in accord with the experiment @xcite .
the equations of motion , eq .
( [ eq : popevolution ] ) , allow us to use light with any spectral composition . in the experiment described in this paper , we use spectrally narrow ( laser ) illumination at 630 nm . for details on the calculations see si .    in figs .
[ fig : excitonic - model]a and [ fig : excitonic - model]b we present the calculated bulk absorption and fluorescence spectra of lhcii monomers at 5 c compared to experimental values taken from @xcite .
we note here that the bulk spectra of lhcii in monomeric and trimeric form are practically identical , see fig .
s1 in si for comparison . in calculations ,
the spectra were averaged over a gaussian disorder of site energies with full width at half maximum ( fwhm ) of 110 @xmath56 .
although the blue chl b shoulder is not perfectly reproduced , the agreement between measured and theoretical absorption is good in the region of our excitation , and the fluorescence ( fl ) spectrum shows a good agreement in general . we therefore conclude that our excitonic model captures correctly the features of the studied system that are the most relevant in the present study .    the fl spectrum in fig .
[ fig : excitonic - model]b is dominated by the lowest four excitons , which are the most populated ones .
these excitons are formed by strongly coupled pigments chl 610 - 611 - 612 , chl 602 - 603 and chl 613 - 14 ( see ref .
@xcite for nomenclature ) .
this is in agreement with previous modeling results for the trimeric lhcii @xcite .
bulk ( a ) absorption and ( b ) fluorescence spectrum of the qy band of lhcii monomers .
the points are experimental values taken from ref .
@xcite , the lines are calculated by our exciton model .
the coloured lines are individual excitonic contributions , the black line is the overal spectrum . ]
the calculated bulk spectra seem to be in a good agreement with the experiment .
our model also reproduces successfully the statistics from the sms experiments . in fig .
[ fig : fl - peak - distribution]a we present the fl spectral peak distribution compared to the experiment .
we can see that the calculated distribution is a little broader than the experimental one , but the agreement is again reasonably good .
the larger distribution width of the calculated spectra can be explained by a relatively long integration time in the experiment ( 1 s ) , during which the system samples several individual realizations of the disorder . in extreme cases the measured values
get averaged towards the mean value . as a result
the measured distribution is narrower . from our comparisons between experiments and calculations
we conclude that the exciton model with gaussian disorder not only reproduces the averaged absorption spectra and equilibriated populations of excitons ( resulting in characteristic fl spectra ) , but also the individual realizations provided by this model are in a good agreement with the experiment .    experimental ( red ) and calculated ( green ) fluorescence peak distribution .
( a ) peak position histogram , ( b ) fl peak position and intensity plot .
theoretical points are calculated as individual realizations of energetic disorder . ]
the next feature of the lhcii single molecular spectra that we need to address is the significant amount of blinking , i. e. reversible switching to the off state . since the measured dwell times in the off state are in the range of seconds , which is significantly longer than the lifetime of any long - living species such as triplet states @xcite
, the off states must correspond to states with efficient excitation energy dissipation .
correspondingly , our model has to be extended by including some fluorescence quenching mechanism .
one of the possible mechanism of fl quenching in lhcii , proposed by ruban and coworkers @xcite , is an excitation energy transfer from the lowest chl a exciton states to a lutein molecule , lut 1 ( lut620 ) , see ref .
@xcite for nomenclature ) .
the lut 1 molecule resides in the vicinity of the so - called terminal emitter group of chlorophylls , composed of chls a610 , 611 and 612 , and it is supposed to be coupled mainly to chl a612 @xcite .
the s1 state of the lut molecule is optically forbidden , and it has a short ( 10 ps ) lifetime due to a decay through a non - radiative channel @xcite .
the transition energy from the s0@xmath57s1 of lut is in the vicinity of the transition to the chl @xmath58 state . due to its short excited state life - time , lut could in principle act as an excitation ( and fluorescence ) quencher .
let us test this mechanism within our model to see if it can account for the observed blinking .
the important parameters of the lutein in context of our model are its s1 state site energy relative to its groundstate and the coupling to chlorophylls , in particular to chl a612
.     role of lut 1 parameters .
dependence of the fluorescence yield on ( a ) lut 1 s1 energy and ( b ) lut 1 - chl a612 coupling .
the energy dependence is calculated with redfield ( blue ) and modified redfield ( red ) theory for comparison .
the dependence on the coupling strength depicts calculated points ( redfield theory ) fitted with exponential dependence .
already a realistically small coupling around @xmath59 leads to significant fl quenching . ]    in fig . [
fig : the - lut-1]a we present the dependence of the relative fl quantum yield on the lut energy for fixed value of the lut - chl coupling of 14 @xmath56 .
the energy dependence agrees well with the one obtained by ruban @xcite .
the quenching is only efficient when the lut energy is below one of the red chlorophylls ( around 15100 @xmath56 ) and the plateau enables lut to act as an efficient quencher even in disordered systems .    in fig .
[ fig : the - lut-1]b we show the dependency of the fl quantum yield on the lut - chl coupling for fixed lut energy 14500 @xmath56 . due to large reorganization energy ,
14500 @xmath56 corresponds to the zero - phonon line at 13900 @xmath56 and thus agrees with experimental observations @xcite . from the coupling dependence of fluorescence in fig .
[ fig : the - lut-1]b we can conclude that weak coupling is sufficient for lut 1 to act as a fluorescence quencher . very importantly
, even small changes in the lut - chl coupling can result in a big difference in the fluorescence intensity .
based on this analysis we decided to use lutein s1 site energy of 14500 @xmath56 in our simulations of blinking .
we define the quenched state by the value of 12 @xmath56 for coupling of the lut 1 to chl 612 and the unquenched state by the zero coupling .
our model allows any type of time dependency of the chl - lut coupling to be used , and it could in principle accommodate input from structure based md studies and quantum chemical treatment of the ( dexter type ) coupling of the chl @xmath58 and the lut s1 states .
however , a much rougher phenomenological model of the chl - lut coupling changes enables a better discussion of the feasibility of the proposed quenching mechanism than the parameter free _ ab initio _ calculations suffering from uncertainties in the structural information .
next we proceed to model the blinking statistics .
as mentioned in the introduction , the blinking statistics alone can be well described by a two - level model proposed by valkunas et al . in ref .
@xcite . by random fluctuations
, the protein samples its potential energy surface ( pes ) performing thus a random walk ( rw ) .
the model of ref .
@xcite assumes that there are two stable conformations of the protein corresponding to two minima of the protein pes .
these are approximated by two harmonic potentials .
the protein undergoes a rw in this potential , and at every step it has a certain probability to switch from its current pes to the other pes . in our treatment
we use a discrete rw description , which enables us to follow individual trajectories of the proteins . for the details of the approach taken in this study and the differences from the original model by valkunas et al . , see si and refs .
@xcite .
to connect the two pes model to our particular lut quenching model , we assume that the change of protein conformation somehow changes the lut - chl coupling .
the lut s1 state does not have a dipole moment , and the resonance coupling similar to those between allowed states does not occur here .
the two different protein conformations responsible for the quenched / unquenched states would then result in two slightly different orientations / positions of the pigments , leading to different strengths of the coupling .
this mechanism is in accord with recent quantum - chemical study by duffy et al .
@xcite , where small configurational changes were found to lead to substantial changes in chl - car couplings .
the switching between the pes is controlled by the rw model with diffusion parameters adjusted to fit the experimental dwell time distributions .
the comparison between the calculated and experimentally determined dwell - times is presented in figs .
[ fig : dwell - times]a and [ fig : dwell - times]b , for the on and off times , respectively .
the agreement is again fairly good letting us believe that our phenomenological model captures the most important features of the protein dynamics affecting the blinking behaviour .     experimental ( red ) and calculated ( green ) probability distribution of dwell times in the ( a ) on and ( b ) off state , logarithmic scale .
while the off state distribution follows a power law , the on state distribution has an exponential tail at longer times . ]      finally , we can connect the two models described above and simulate the blinking behaviour . to this end
we continuously model the fluorescence of the lhcii complexes , and output the intensity ( and spectrum ) every 10 ms corresponding to the experimental integration bins .
simultaneously , we let the protein do the rw on its pes , and we adjust the chl / lut coupling , when the protein switches between pes . to obtain more realistic traces , either the site energies or couplings
can be slightly varied after every jump , resulting in different energy levels .
such a procedure , however , does not lead to qualitatively different conclusions and it can be in principle omitted . in the calculations presented here , we used gaussian disorder with the fwhm of 0.3 @xmath56 for couplings and 1 @xmath56 for energies . for every realization of the disorder
a 60 s trace is modelled .
this is repeated for 200 realizations of the disorder , reflecting the experimental conditions .
the resulting statistics are presented in fig .
[ fig : intensity - traces - statistics ] .
total dwell times in fig .
[ fig : intensity - traces - statistics]a represent the overall amount of time spent in a given intensity level . from this
we can see the two - level character of the blinking and simultaneously also the presence of some intermediate levels , which result from particular disorder realizations .
the dwell times are similar for the on and off states .
[ fig : intensity - traces - statistics]b shows how often are the intensity levels visited per fixed amount of time .
the modelled access frequency distribution is naturally very symmetrical , a direct result of the fact that in the model the complex switches only between the on and off state .
the number of on / off states visited per minute therefore has to be the same .
the experimentally analyzed intensity traces contain also jumps between levels within the on / off states in an amount which can be , to some extent , modified by adjusting the sensitivity of the level analysis .
the presence of these intra - state jumps results in higher switching frequency and wider distribution in intensities , causing a moderate discrepancy between experiment and simulation . in order to include this kind of switching into the model
, dynamic sampling of the disorder would have to be incorporated .
work in this direction will be presented elsewhere .    for the reasons stated above , the level access
frequency distribution is not well suited for comparison of the model with the experiment .
a more appropriate measure of the blinking would be the intensity - intensity correlation function defined as @xmath60 .
this quantity is well - known from single - molecule measurements , where it is often used to characterize the blinking behaviour @xcite . in fig .
[ fig : intensity - traces - statistics]c we present @xmath61 obtained from 50 measured long enough traces , compared with the model .
we can see that the agreement between experiment and theory is good , indicating that our model gives reasonable switching between the intensity levels .
the shape of the correlation function is given by the dwell time statistics , see also fig .
[ fig : dwell - times ] .
the initial fast drop implies the abundance of short blinking events .
this results from the mechanism of the protein switching between its potential surfaces , where the short succesive blinking events are caused by the dynamics in the vicinity of the pes intersection .
experimental ( red ) and calcuated ( green ) intensity blinking statistics .
( a ) the percentage of time the complexes dwelled on respective intensity levels .
the two - state structure of the low - intensity off states and higher intensity on states is apparent .
( b ) how often per minute the complexes accessed the respective intensity levels .
the experimental frequency is higher due to switching within the on / off states .
( c ) the intensity correlation function . ]
in section [ sub : excited - state - dynamics ] of the results we derived approximate equations of motion for populations evolution on the timescale of the sms experiment . we have shown that the weak illumination regime , in which the sms experiment is performed , allows for an effective source term description of the light - matter interaction in which short time transient effects arising in the photoinduced evolution of molecular systems can be consistently accommodated side by side with the slow evolution of the protein bath observed in sms experiments .
we demonstrated the validity and scope of application of our equations by simulating our single molecule experiment on lhcii monomers .
based on the recent research in elucidating the npq mechanism in lhcii and the connection between fluorescence quenching and energy dissipation we implemented energy transfer to lut 1 as a blinking mechanism .
based on our calculations we were able to confirm the findings of duffy et al .
@xcite and chmeliov et al.@xcite that , within a reasonable range of parameters , lut 1 can indeed act as an efficient quencher . since our model extends the previous treatment by including realistic excitation conditions and population transfer rates , it is remarkable , how similar our fluorescence quenching dependence on the lut 1 energy is to the one in ref .
moreover , we were able to confirm that lut 1 acts as an efficient quencher also under am1.5 illumination ( data not shown since the dependence is very similar ) . at the same time we can see that the amount of quenching is very sensitive to the change of coupling of the lut to the chlorophylls . since the coupling itself is very sensitive to the distance and orientation between the pigments , it provides a possible link to the protein conformation change working as a switching mechanism as proposed in ref .
indeed , when using the 2-level switching model to control the change of coupling , we are able to reproduce the experimentally obtained blinking statistics .
although far from being exclusive in any way , our argument strongly supports the notion of the protein acting as a conformational switch regulating the amount of quenching in the system .
the agreement between the theory and experiment also serves as a good demonstration of the scope of our equations .
they provide a description for controlling the energy transfer in the system by modulating the parameters of the excitonic model .
we are aware of the remaining phenomenological nature of our connection of the 2-level switching model with the excitonic model .
further improvements in the direction of introducing more parameters with a particular physical meaning , for example the relation to the actual pes shape , are needed and will be subject to further study . also recent experimental observations indicate the presence of more relevant timescales in the intensity traces suggesting the inadequacy of a 2-level model with one reaction coordinate . finally ,
although some connection between the fluorescence blinking and npq was already shown by krger et al .
@xcite , their exact relation is yet to be elucidated .
the equations derived here are a suitable tool for these future investigations .
[ [ section ] ]
trimeric lhcii complexes of spinach were isolated from freshly prepared thylakoid membranes as described earlier @xcite .
monomeric complexes were obtained by incubating lhcii trimers with 1% ( w / v ) octyl glucoside and 10 @xmath62 phospholipase a2 ( sigma ) @xcite .
subsequent fast protein liquid chromatography ( fplc ) ensured a homogeneous sample preparation .
the ensemble fluorescence absorption and emission spectra were measured on a lambda40 spectro - photometer ( perkin - elmer ) and a fluorolog tau-3 ( jobin yvon ) , respectively . for sms experiments ,
the sample was diluted down to a concentration of ~ 10 pm in a measuring buffer ( 20 mm hepes , ph 8 and 0.03% ( w / v ) n - dodecyl @xmath63-d - maltoside ) and then immobilized on a pll ( poly - l - lysine , sigma ) coated cover glass .
addition of an oxygen scavenging mix of 750 @xmath64 glucose oxidase , 100 @xmath65 catalase and 7.5 @xmath66 glucose ( all sigma ) to the closed sample chamber inhibited the formation of highly reactive singlet oxygen and improved the photostability of complexes .
a confocal microscope was used to measure the fluorescence of single complexes at 5 @xmath67 as described previously @xcite .
the sample was excited at 630 nm utilizing a ti : sapphire laser ( coherent mira 900f ) with a pulse width of 200 fs and a repetition rate of 76 mhz coupled to a tunable optical parametric oscillator ( coherent mira opo ) .
near - circular polarized light was obtained by utilizing a berek polarization compensator ( 5540 new focus ) . a fluorescence beam splitter ( 70:30 , thorlabs ) allowed us to simultaneously measure the fluorescence spectrum _ via _ a ccd camera ( spec10:100br , roper scientific ) with an integration time of one second and the wavelength integrated fluorescence intensity _ via _ an avalanche diode ( spcm - aqr 16 , perkin elmer ) with a binning time of 10 milliseconds .
the fluorescence of one complex was analyzed for either one minute or until it photo - bleached and a set of 200 complexes served as the basis for statistical analysis .
the fluorescence peak distribution was obtained by fitting of a skewed gaussian to the fluorescence spectrum as shown in kruger et al . @xcite and
the blinking analysis was performed equivalently to the algorithm described elsewhere @xcite .
the equations of motion , eq .
( [ eq : popevolution ] ) , have quasi constant coefficients , and they can therefore be written in the form    @xmath68    where @xmath69 is a matrix of relaxation and population transfer rates , and @xmath70 are the source terms .
( [ eq : matrix_eq-1 ] ) can be solved analytically :    @xmath71    this expression enables us to find the populations at any time without the need to solve for all the previous times .
if we aim at steady state , @xmath22 can be send to @xmath72 and the populations then depend only parametrically on time @xmath40 .
the weak dependence of @xmath73 and @xmath69 on @xmath40 makes it possible to explain changes in the populations of the emitting states of a molecular system due to slow changes of the protein environment and the structure of the molecular system .
10 url # 1`#1`urlprefix[2]#2 [ 2][]#2    _ _ ( , ) .    & .
_ _ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .    , & _ _ ( , , ) .    , & .
_ _ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .
_ et  al . _ _ _ * * , ( ) .
_ et  al . _ . _
_ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .    , , , &
_ _ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .
et  al . _ .
_ _ * * ( ) .    & .
_ _ * * , ( ) .    , , , & .
_ _ * * , ( ) .    , , , & .
_ _ * * , ( ) .
_ _ * * , ( ) .    & .
_ _ * * , ( ) .
_ et  al . _ . _ _ ( ) .
_ et  al . _ . _
_ * * , ( ) .    , , & .
_ _ * * , ( ) .    , , & _ _ * * , ( ) .    , & ( ) .    , , & _ _ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .    , , & _ _ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .    , & _ _ * * , ( ) .    , , , & .
_ _ ( ) .    & _ _ ( , ) .    , & _ _ ( , ) .
& _ _ * * , ( ) .
_ _ * * , ( ) .    & .
_ _ ( ) .    , , & .
_ _ * * , ( ) .    &
_ _ * * , ( ) .    & _ _ * * , ( ) .    , , , & .
_ _ * * , ( ) .    , , & .
_ _ * * , ( ) .
_ et  al . _ . _ _ * * , ( ) .    _
_ ( , , ) .    & .
_ _ * * , ( ) .    & .
_ _ * * , ( ) .    & _ _ * * , ( ) .    , & .
_ _ * * , ( ) .    & .
_ _ * * , ( ) .    & .
_ _ * * , ( ) .
_ _ ( , , ) .    .
_ _ * * , ( ) .    , , & .
_ _ * * , ( ) .    & .
_ _ * * , ( ) .    , , & .
_ _ * * , ( ) .
_ et  al .
_ _ _ * * , ( ) .    , & .
_ _ * * , ( ) .
_ _ * * , ( ) .    , , , &
_ _ * * , ( )
.    _ et  al . _ .
_ _ * * , ( ) .    , , & ( eds . ) _ _ ( , ) .    , , , & .
_ _ * * , ( ) .
_ et  al . _
_ _ * * , ( ) .
_ _ ( , ) .
p.m. , j.m.g . and r.v.g .
were supported by the vu university and by an advanced investigator grant from the european research council ( no . 267333 , photprot ) to r.v.g .
was also supported by the nederlandse organisatie voor wetenschappelijk onderzoek , council of chemical sciences ( nwo - cw ) via a top - grant ( 700.58.305 ) , and by the eu fp7 project papets ( ga 323901 ) .
r.v.g . gratefully acknowledges his academy professor grant from the netherlands royal academy of sciences ( knaw ) .
p.m. and t.m .
received financial support from the czech science foundation ( gacr ) , grant no .
14 - 25752s .
t. m. acknowledges the kind support of the nederlandse organisatie voor wetenschappelijk onderzoek ( nwo ) visitor grant nr . 040.11.423 during his stay at the vu amsterdam .
the interaction of the pigments , i. e. the system , with the vibrations , i. e. its environment , is treated by second order perturbation theory . there
the bath is completely described by its correlation function or , equivalently , in the frequency domain , by the spectral density of the bath vibrations .
we use the spectral density from ref .
it is constructed from one overdamped oscillator representing the slow protein vibrations and 48 high - frequency intrapigment modes ( for more context for the spectral density and excitonic model see ref .
@xcite ) :      the ( temperature dependent ) correlation function given by this spectral density is assumed to be uncorrelated between individual sites and differs only between chl a and chl b , while the difference is in the coupling strength @xmath75 :      the time - dependent correlation function is then obtained by fourier transform of ( [ eq : correl.func . ] ) .
we use the same values of the spectral density parameters as in ref .
the difference between the vertical , franck - condon transition of the pigments , which are called site energies in this text , and their 0 - 0 transitions is given by the reorganization energy due to the interaction with the bath :      because the pigments are strongly coupled , the preferred basis of calculations is the excitonic basis in which the system hamiltonian is diagonal with eigenvalues , exciton energies , @xmath78 .
all quantities , including correlation functions @xmath79 , reorganization energies @xmath9 and transition dipole moments @xmath19 have to be transformed into the excitonic basis :                                        the population relaxation rates of chlorophylls are a result of energy transfer and radiative and non - radiative decay , i.e. @xmath92 , where @xmath93 , @xmath94 is inverse lifetime @xmath95 of site @xmath96 excited state , which is taken to be 3 ns for all @xmath96 .      to check for sample degradation a control bulk measurement of absorption and fluorescence was performed . in fig .
[ fig : bulk - spectra - measured ] the measured spectra are given together with the spectra of lhcii trimers taken from @xcite .
their excellent agreement confirms absence of degradation in course of the experiment and also justifies the usage of the trimer spectra for the bulk spectra modelling in the main text .        in the original description by valkunas @xcite and , in more detail
, in @xcite , the protein diffusive motion was described as a continuous - time random walk ( ctrw ) on a two - dimensional potential energy surface . here
we somewhat simplify this description in the following way .
the two generalized coordinates represent fast and slow degrees of freedom and thus , when following the slow dynamics , the fast fluctuations can be adiabatically eliminated .
it can be shown that when we consider the potential dependence on the fast coordinate the same in both the on and off states for simplicity , i.e. setting @xmath97 in @xcite , the fast coordinate can be completely removed .
this leaves us with effectively one - dimensional problem .
@xmath104 is a rate of falling from the @xmath105 potential , the first exponential term reflects the energy gap law and the min term ensures detailed balance condition .
the coefficient @xmath106 can be treated as a constant and @xmath107 is a characteristic frequency of the protein environment vibrations responsible for the tunneling .    the protein diffusion under this conditions can be described either by the ctrw or by a discrete random walk ( rw ) .
the former approach was employed by valkunas et al.@xcite .
however , we believe that for our purpose it is better to solve this problem as a discrete rw for two reasons . first ,
the coupled smoluchovski equations for the ctrw on the two potentials can be decoupled only for conditional probabilities , i.e. assuming that the system was in the opposite state in the previous interval , and , in the same time , employing the same , equilibrium initial condition for each dwell time .
when continuously modelling the trajectory of a single protein , we do not have to include the resetting after switching and also the conditioning will be inherent , as the system is observed being in a particular state . and second , when we want to simulate the individual intensity time traces , it is more natural to really follow the trajectories of the individual proteins on their pes .
if we want to follow the time trace of every single molecule , we should follow its particular trajectory .
the blinking statistics will then be recovered by averaging over a large number of molecules , exactly as in the experiment . to the purpose of following trajectories of individual proteins
, we need to describe its discrete rw ( drw ) in the potential .
we will denote probability of going right ( left ) as @xmath108 ( @xmath109 ) . in the symmetrical rw we have @xmath110 . if the protein is at a coordinate @xmath111 , in the next step it will move with probability @xmath108 to @xmath112 and with probability @xmath109 to @xmath113 , where @xmath114 is the length of the step . inspired by classical approach by van kampen @xcite , we augment the position dependent probabilities in the presence of the potential @xmath115 as      we note that defined probabilities defined in eq .
( [ eq : pprobandqprob ] ) reflect the detailed balance condition @xmath117 . considering a small step @xmath114
, we can use taylor expansion in @xmath111 @xmath118 , obtaining                where @xmath125 is the time duration of the step . during @xmath125
the protein walks either to the left or right with the respective probability and with the probability @xmath126 switches to the other surface .    similarly to @xcite the position near the potential intersection
can be chosen as an initial condition . however , as we do not include resetting after switching in our approach , this determines only the starting point of each trajectory and is therefore not of significant importance . in table
[ tab : rw - model - paramaters ] we present the parameters of our rw model compared to the ones used by valkunas et al @xcite . considering the differences - lhcii monomers vs trimers , discrete vs continuous rw ,
no resetting vs resetting - the agreement is satisfactory .