Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: photosynthesis is the natural mechanism for the capture and storage of energy from sunlight by living organisms . excitation energy is absorbed by pigments in the photosynthetic antennae and subsequently transferred to a reaction center where an electron - transfer event initiates the process of biochemical energy conversion . in certain bacterial systems and higher plants light harvesting efficiency is indeed above 99% @xcite . although this phenomenon has been studied for decades damjanovi97,ritz01,grondelle04,cho05 , a full description of the underlying mechanism leading to this remarkably high efficiency is yet not available . it has already been demonstrated experimentally that the excitation energy transfer within chromophoric arrays of photosynthetic complexes could involve quantum coherence under certain physical conditions . in particular , this phenomenon has been observed via electronic spectroscopy of delocalized exciton states of light - harvesting complexes @xcite and the fenna - matthews - olson ( fmo ) protein complex @xcite . the energy transfer mechanism in multichromophoric arrays can often be described by a semiclassical frster theory which involves incoherent hopping of the excitations between energy levels forster65,scholes03,maybook . in this method , the coulomb interaction among different sites is treated perturbatively to calculate the probability of exciton hopping . the more general approach for including coherent effects is given by redfield theory which provides a microscopic description of excitation dynamics via a master equation in a reduced space of excitons in the weak phonon coupling and born - markov approximation @xcite . an equivalent approach to redfield theory for calculating the diffusion constant of excitons was proposed by silbey and grover @xcite . alternative methods to explore coherent and incoherent exciton transfer were also introduced using a stochastic model ( haken and strobl @xcite ) , and a generalized master equation formalism ( kenkre and knox @xcite ) .    in order to study the nonlinear spectroscopy of molecular aggregates , zhang _ et al . _ @xcite introduced a modified redfield equation for statically disordered exciton systems which treats the diagonal elements of exciton - bath coupling in a non - perturbative fashion . this approach was later used to model energy transfer dynamics in light - harvesting complexes of higher plants @xcite . yang and fleming provide a comprehensive comparison of frster , standard redfield , and modified redfield theories in ref . a generalized theory for multichromophoric frster resonance energy transfer which includes coherence effects within donors and acceptors , while considering donor - acceptor interactions according to the standard frster model was also proposed by jang , newton , and silbey jang04,cheng06,jang07 . in another study , the effects of geometry and trapping on energy transfer were examined in simple chromophoric arrays within the haken - strobl model @xcite . recently , direct evidence of quantum coherence in the dynamics of energy transfer has been observed experimentally in the fmo complex @xcite and also in the reaction center of purple bacteria @xcite . these previous studies led us to further explore and characterize quantum interference , decoherence effects , and their interplay within the dynamics of photosynthetic complexes as potential mechanisms for the enhancement of the energy transfer efficiency . here , we develop a quantum walk approach , based on a quantum trajectory picture in the born - markov and secular approximations , as a natural framework for incorporating quantum dynamical effects in energy transfer , as opposed to a classical random walk picture that can effectively describe the excitation hopping in the frster model . the concept of quantum walks originated by feynman works in connection with diffusion in quantum dynamics , in particular to model the dynamics of a quantum particle on a lattice @xcite , and also path integral formalism for discretizing the dirac equation @xcite . continuous - time quantum walks were also used by klafter and silbey to find hopping time distribution functions in exciton dynamics klaftersilbey80 . the formal discrete and continuous - time approaches to quantum walks  were developed later , e.g. , see refs . aharonov93,kempe03 , including some investigations of model decoherence effects @xcite . purely _ unitary _ continuous - time approaches to quantum walks were employed in the context of quantum information science , where they yield potential exponential speedups over classical algorithms @xcite . the quantum walks are of particular interest as potential computational tools @xcite , and applications to quantum cellular automata @xcite , quantum optical systems @xcite , and coherent excitation transport @xcite .    in this work , we develop a theoretical framework for studying the role of quantum coherence in energy transfer dynamics in molecular systems within the born - markov approximation in the lindblad formalism . our approach is essentially equivalent to a redfield theory with the secular approximation . however , our approach naturally leads to quantum trajectory picture in a fixed - excitation reduced hilbert space that can be described by the concept of directed quantum walks in liouville space . quantum walks in actual physical systems differ from idealized models of quantum walks in several significant ways . first , hamiltonians of physical systems typically possess energy mismatches between sites that lead to anderson localization @xcite . second , actual quantum walks are subject to relatively high levels of environment - induced noise and decoherence . the key result of this paper is that the interplay between the coherent dynamics of the system and the incoherent action of the environment can lead to significantly greater transport efficiency than coherent dynamics on its own . we introduce the concepts of energy transfer efficiency ( ete ) and its susceptibility and robustness and explore the dynamical effects of coherent evolution and environmental effects at various temperatures from a microscopic hamiltonian formalism . for the fmo protein , we show that a grover - type quantum search @xcite can not explain the high ete of this complex . however , we demonstrate that a directed quantum walk approach can be used for studying the energy transfer efficiency as a function of temperature , reorganization energy , trapping rate , and quantum jumps from sites to sites . moreover , we explore similar dependencies for the susceptibilities of ete with respect to basic processes contributing to the fmo dynamics including the free hamiltonian , the phonon bath  jumps , dephasing in the energy basis , transfer to the acceptor , and exciton decay . we demonstrate that the efficiency increases from 70% for a purely unitary quantum walk to 99% in the presence of environment - assisted quantum jumps . this article is organized as follows . in sec . [ lindblad ] ,  we develop a lindblad master equation in the site basis for studying energy transfer of multichromophoric channel systems in the born - markov approximation . in sec . [ qw ] , we introduce a quantum walks formalism in liouville space to describe the energy transfer pathways . the definition of ete is presented in sec . [ ete ] . in sec . [ fmo ] , we apply our theoretical approach for studying the dynamics of fmo complex . some concluding remarks are given in sec . [ conclusion ] . the fenna - matthews - olson protein acts as an energy transfer channel in the biological process of photosynthesis connecting the base plate of the antenna complex to the reaction center of green sulfur bacteria . this type of functional role of an interacting multichromophoric system can be formalized by the hamiltonian for an consisting of @xmath0 donors , @xmath1 channel chromophores , and @xmath2 acceptors as : @xmath3the @xmath4 and @xmath5 are the creation and annihilation operators for an electron - hole pair ( exciton ) at chromophore @xmath6 and @xmath7 are the site energies ( not including the bchl a transition frequency @xmath8@xmath9 ) and @xmath10 . the @xmath11 are coulomb couplings of the transition densities of the chromophores , often taken to be of the frster dipole - dipole form , @xmath12 with @xmath13 the distance between site @xmath6 and @xmath14 and @xmath15 the transition dipole moment of chromophore @xmath6 @xcite . note that in systems where chromophores are closely packed ( e.g. , the fmo complex of green sulfur bacteria @xcite ) or the site energies are ( almost ) resonant ( e.g. , the lh1 ring of purple bacteria @xcite ) , @xmath7 can be of the same order of magnitude as @xmath16 such cases require a non - perturbative treatment of the coupling . in this work , the @xmath17 include all intra - donor / channel / acceptor couplings and inter - chromophoric couplings for donor - channel and channel - acceptor . here , we ignore the @xmath11 for inter donor - acceptor coupling due to large spatial separation , as for example chlorosomes and reaction center in the green sulfur bacteria @xcite . thus , excitation transfer from donor to acceptor always occurs via the channel . we also assume that for donor - channel and channel - acceptor coherent couplings are weak , i.e. , @xmath18 for couplings into and out of the channel . this implies that the energy transfer to and from the channel can be described by semi - classical frster theory .    in this work , we study the role of the fmo protein acting as an energy transfer channel in green sulfur bacteria . thus , we focus on the dynamics of a chromophoric channel of @xmath1 sites with denoting the free hamiltonian in the reduced channel hilbert space as @xmath19 the hamiltonian @xmath20 is formally equivalent to the hamiltonian of eq . ( generalfreehamiltonian - eq ) with @xmath21 we consider only the zero and single excitation manifolds given by the states @xmath22 and @xmath23 . we denote the eigenbasis of the hamiltonian @xmath24 as exciton basis @xmath25 , where @xmath26 . here , the effect of the donor is modeled by a static initialization of the channel . to account for the channel - acceptor coupling we introduce an effective non - hermitian channel - acceptor hamiltonian , @xmath27 , which can be obtained by a projector - operator method analogous to ref . @xcite , with the acceptor transfer rates @xmath28 where @xmath29 denotes the density of states for the acceptor .    in general , the multichromophoric channel is subject to a thermal phonon bath and a radiation field . the interaction hamiltonian can be written as @xmath30 with the phonon coupling @xmath31where @xmath32 is an operator acting in the bath hilbert space . the exciton__-__phonon and exciton - thermal photon interaction @xmath33 describes the coupling of the bath operators @xmath34 to the transition dipole moment of each chromophore , i.e.,@xmath35the phonon terms @xmath36 induce relaxation and dephasing without changing the number of excitations . in other words , the state of the multichromophoric system remains in a fixed excitation manifold under the evolution generated by @xmath37 . the @xmath33 hamiltonian leads to transitions between exciton manifolds . generally , one can also consider diagonal static disorder in the hamiltonian as : @xmath38 this disorder can be generated for example by variation in the structure of the protein environment in the time scales which are usually much slower than excitation transfer time scale of @xmath39 ps @xcite . the dynamics of the system to second order in the system - bath coupling , can be described by the lindblad master equation in the born - markov and secular approximations as @xcite :    @xmath40+l_{p}(\rho ( t))+l_{r}(\rho ( t ) ) . \label{master eq}\ ] ]    where @xmath41 are the lamb shifts due to phonon and photon - bath coupling . the respective lindblad superoperators @xmath42 and @xmath43 are given by  ( @xmath44 ) :    @xmath45 ,   \nonumber\end{aligned}\ ] ]    for @xmath37 the corresponding lindblad generators are @xmath46 , where the summation is over all transitions with frequency @xmath47 in the single - excitation manifold and @xmath48 denotes the exciton with frequency @xmath49 . the secular approximation is valid when the relevant time scale of the intrinsic evolution of the system , @xmath50 , is much faster than the relaxation time scale . for example , an energy difference of 200cm@xmath9 between excitons , e.g. , in the fmo complex , translates to a time scale of 0.16 ps which is much smaller than the typical 1 ps time scale of energy relaxation in the single - excitation manifold . the rates @xmath52 are given by the fourier transform of the bath correlation function as @xmath53 where we assume that off - diagonal fluctuations are small compared to diagonal fluctuations @xcite . the rate can be further simplified to the site - independent expression @xmath54 @xmath55 $ ] where @xmath56 $ ] is the bosonic distribution function at temperature @xmath57 . here , we assume an ohmic spectral density with  @xmath58 for @xmath59 and @xmath60  elsewhere , with cutoff @xmath61 and reorganization energy @xmath62 @xcite .    in the case of the hamiltonian @xmath33 we obtain the lindblad generators @xmath63 where @xmath64 is the molecular transition frequency , separated into @xmath65@xmath9 for the @xmath66 band of bacteriochlorophyll @xmath67 and excitonic energies @xmath68 of the order of 300@xmath69 @xcite . the respective rate is again assumed to be diagonal , @xmath70 @xmath71 and site independent , @xmath72 for ohmic @xmath73 and super - ohmic spectral densities one is able to approximate @xmath74 this approximation yields a simplified lindblad superoperator similar to eq . ( [ lindblad superoperator ] ) without the summation over frequencies and with the generators @xmath75 finally , we choose  the rates @xmath76 such that it leads to a 1ns exciton life - time , as experimentally measured for chromophoric complexes e.g. , in ref . the lamb - shifts @xmath77 are explicitly given by @xmath78 ( for @xmath79 ) where @xmath80 is the imaginary part of the half - sided fourier - transform of the bath - correlator . the lamb shift usually contributes only marginally to the dynamics of the system , e.g. in the fmo complex @xcite . the master equation eq . ( [ master eq ] ) and its lindblad superoperators eq . ( [ lindblad superoperator ] ) contain a significant degree of complexity , reflecting the non - trivial form of the fmo complex and its diverse sources of environmental interaction . to deal with this complexity one needs to separate the contributions of the different physical processes to the dynamics . to this end , in the next section , we explicitly construct a quantum trajectory master equation to study the effects of free hamiltonian , damping and relaxation . specifically , we demonstrate that the energy transfer in multi - chromophoric systems of the chromophoric channel can be considered as a generalized ( directed ) continuous - time quantum walk in the single - excitation manifold interrupted by jumps to the zero - excitation manifold . in general , the frster theory for energy transfer leads to a classical random walk description of the transport in photosynthetic units sener02,sener04 . the equation of motion for the classical probabilities of an excitation being at site @xmath67 , @xmath81 is given by @xmath82the @xmath83 denotes the markov transition matrix elements which describes the classical frster rates between site @xmath67 and @xmath84 . in  closely - packed chromophoric arrays , however , one has to consider not only populations of states but also coherence between states ; consequently , the equation of motion is given by a master equation for the density matrix such as eq . ( master eq ) . for the purposes of evaluating the dynamics of this master equation , and for investigating the interplay between coherence and decoherence in the resulting quantum walk , it is convenient to re - express the master equation for a single excitation in terms of a quantum trajectory picture of open quantum system dynamics , as in @xcite : @xmath85^{\star }   \label{quantumtrajectory - eq } \\ & & + \sum_{m , m^{\prime } , n , n^{\prime } } ^{n_{c}}\gamma ^{p}(m , m^{\prime } , n , n^{\prime } ) w_{m , m^{\prime } } \rho ( t)w_{n , n^{\prime } } ^{\dagger } \nonumber \\ & & + \sum_{m}^{n_{c}}\gamma _ { m}^{r}r_{m}\rho ( t)r_{m}^{\dagger } , \nonumber\end{aligned}\]]where the @xmath86 generate jumps in the single - exciton manifold and the @xmath87 generate jumps between exciton manifolds . the jump rates are given by@xmath88  and @xmath89 . we have also defined @xmath90^{\star } $ ] as a generalized commutation relation for any two operators @xmath91 and @xmath92 as @xmath93^{\star } = ab - b^{\dagger } a^{\dagger } $ ] . in order to describe absorption of the exciton at the acceptor site and the damping due to the phonon bath , we introduced an effective anti - hermitian hamiltonian as : @xmath94 , with @xmath95    now taking the trace of the overall master equation leads to the probability density that no jump to the zero manifold occurs between time @xmath96 and @xmath97 as : @xmath98dt+\sum_{m , m^{\prime } , n , n^{\prime } } ^{n_{c}}\gamma ^{p}(m , m^{\prime } , n , n^{\prime } ) \mathrm{tr}[w_{m , m^{\prime } } \rho _ { eff}(t)w_{n , n^{\prime } } ^{\dagger } ] dt,$ ] and the probability density of a jump event becomes :  @xmath99dt.$ ] in the case of a no - jump trajectory one has the master equation @xmath100^{\star } + \sum_{m , m^{\prime } , n , n^{\prime } } ^{n_{c}}\gamma ^{p}(m , m^{\prime } , n , n^{\prime } ) w_{m , m^{\prime } } \rho _ { eff}(t)w_{n , n^{\prime } } ^{\dagger } , $ ] which can be considered as a directed quantum walk on the one - exciton manifold described by the density operator @xmath101 . * *    we note that the contribution due to the free hamiltonian , i.e. , the equation @xmath102 $ ] represents a continuous - time unitary quantum walk in the single - excitation manifold . the off - diagonal elements of the free hamiltonian terms represent quantum coherent hopping with amplitudes @xmath103 between the localized sites with energy @xmath104 the analogy of continuous - time unitary quantum walks to classical random walks was addressed within the context of quantum algorithms in ref . @xcite and further explored in childs02,kempe03 . the damping contribution to the dynamics within @xmath105^{\star } $ ] leads to site - dependent relaxation for both diagonal and off - diagonal elements of the density operator and the terms @xmath106 induce quantum jumps in the single - excitation manifold . the similarity between quantum walks and classical random walks , eq . ( classical walk eq ) , can be readily emphasized in liouville space . the quantum trajectory picture is re - expressed as a matrix equation by representing the density operator in a vectorial form . thus , quantum walks in multi - chromophoric complexes are described by the equation@xmath107where @xmath108 this equation manifests a quantum walk in the @xmath109-dimensional liouville space , with the transition ( super- ) matrix , @xmath110the time variation of the vector @xmath111 in liouville space represents the dynamics of population elements @xmath112 as well as coherence elements @xmath113 which are the signature of quantum dynamics . the real part of @xmath114 is responsible for the directionality in the quantum walk which enhances the excitation energy transfer for the fmo complex as we will demonstrate later . this is due to the * *  * * _ effective _ interplay * *  * * between the free hamiltonian and phonon - bath coupling that generates quantum jumps . next we study the properties of the quantum walk picture in terms of the energy transfer efficiency as an universal measure . the energy transfer efficiency of the channel is defined as the integrated probability of the excitation successfully leaving the channel to the acceptor . this definition of the channel efficiency has the advantage of being independent of the detailed dynamics within the acceptor , charge separation , and the energy storage via a chemical reaction . more formally we define :    @xmath115    similar definitions ,  using _ integrated _ success probabilities , were used in context of energy transfer from donor to acceptors ritz01,castro07 , and quantum random walks @xcite . note that the efficiency in eq . ( [ ete eq ] ) has the upper limit of @xmath116 yet in most practical cases there is always a natural cutoff , @xmath117 due to finite excitation life - time , e.g. , for chromophoric complexes and gaas quantum dots @xmath118ns sener02,nazir05 . moreover , the relevant dynamical time scale for the excitation to be transferred to the acceptor is about @xmath119 a more quantitative measure for the transfer time through the quantum channel to the acceptor is given by,@xmath120    in this work we mainly focus on the energy transfer efficiency . in order to study its optimality and robustness , we define a decomposition @xmath121 of the superoperator by @xmath122 @xmath123the different @xmath124 can be chosen in any desired way in order to isolate the effect of different parts of the master equation . for example , one @xmath121 could represent the effect of the system hamiltonian on its own , while another could represent the effect of the phonon bath . each term @xmath121 is associated with a scalar quantity @xmath125 , which represents the overall strength of the of the @xmath126th term in the dynamics . for example , the @xmath125 associated with the hamiltonian gives the overall energy scale of the hamiltonian , while the @xmath125 associated with the photon bath gives the strength of the coupling to that bath . it is natural to investigate the energy transfer efficiency in terms of the susceptibilities @xmath127 and the hessian @xmath128 using the scaling parameters @xmath129 @xmath130 in the neighborhood of @xmath131 @xmath132 for any chromophoric complex the gradient of the efficiency defined by the set of @xmath133 can be used as a measure for local optimality of its performance with respect to independent parameters @xmath134 moreover , such efficiency gradient can be utilized for engineering photovoltaic materials with optimal and robust energy transfer in their respective parameter space . one can verify that the susceptibilities satisfy the identity @xmath135 where the sum consists of a complete decomposition of the superoperator . this identity is the consequence of energy conservation for all times @xmath96 : in the single - excitation manifold , the energy that is transferred by the physical processes @xmath136 out of the chromophoric channel is absorbed by the reaction center in a given time interval . additionally , the hessian @xmath137 is a measure for the second - order robustness of natural or engineered chromophoric systems . note that the robustness also satisfies a similar conservation property , i.e. , @xmath138 next , we study the quantum walk approach in the context the ete in the fmo complex . the structure of the fenna - matthews - olson ( fmo ) complex of green sulfur bacteria was revealed by x - ray crystallography @xcite , as the first pigment - protein complex structure to ever be determined in this method , and since then it has been extensively studied @xcite . the fmo complex consist of a trimer , formed by three identical monomers , each constituting of seven bacteriochlorophyll molecules ( bchl a ) supported by a rigid protein backbone . the fmo complex essentially acts as a molecular wire , transferring excitation energy from the chlorosomes , which are the main light - harvesting antennae of green sulfur bacteria , to the membrane - embedded type i reaction center . in the recent study of the _ chlorobium tepidum _ fmo complex by engel _ et al . _ , @xcite , using two - dimensional electronic spectroscopy , direct evidence of long - lived coherence in the form of quantum beatings was demonstrated at 77k . the presence of quantum coherence prompted speculations about the presence of quantum computation in fmo . indeed , it was argued that fmo acts as a dedicated computational device @xcite , since excitons are able to explore many states simultaneously and select the correct answer , which here is the lowest energy excitonic state . this operation was claimed to be _ analogous _ to grover s algorithm , which is known to provide a quadratic speed - up over its classical counterparts for searching elements of an unstructured database @xcite . here , we argue that a purely unitary grover - type search algorithm can not explain the efficiency of the exciton transfer in fmo complex . however , we employ our directed quantum walk approach to study the quantum effects in the dynamics of the fmo complex . we use the free hamiltonian of the fmo complex as given in ref . the site - energy differences and inter - site couplings lead to exciton energies with separations of around @xmath139cm@xmath9 . the highest energy states are exciton 6 and 7 which are mainly delocalized over sites 5/6 and 1/2 , respectively . the lowest exciton state 1 involves the site 3 . for an overview of the structure of the fmo complex see fig . [ fig3 ] ( a ) . the initial states for the simulation are taken to be sites 1 and/or 6 which are close to the chlorosome antenna ( donor ) @xcite . transfer of the excitation from the fmo channel to the acceptor occurs via site 3 with the rate @xmath140 , which is a free parameter in our simulations and , if not otherwise stated , taken to be @xmath141ps@xmath142 the exciton lifetime is assumed to be @xmath143ns @xcite . the bath spectrum is taken to be as described above with the reorganization energy @xmath144cm@xmath9 and cutoff @xmath145@xmath146 inferred from fig .  2 of ref . @xcite . the purely unitary evolution generated by the seven - site hamiltonian is not performing a grover - type search . in general , for a unitary evolution to be qualified as a quantum search a certain set of conditions have to realized : ( i ) @xmath147 , where @xmath148 is an equal superposition of the basis states @xmath149 of the hilbert space including the solution the search , the target state @xmath150 . @xmath151 is the pure initial state of the system at some arbitrary time @xmath152 and @xmath153 . for a standard grover search algorithm @xmath154 is equal to zero . this condition assures that the initial input state of the search is unbiased with respect to the target state . ( ii ) @xmath155 where @xmath156 is the smallest time in which the free evolution of the initial state , @xmath157 , has significant overlap with the target state , and @xmath158 . ( iii ) @xmath159 where @xmath160 is either a constant or a @xmath161 in summary , a hamiltonian generating the grover algorithm should map an equal coherent superposition of all possible ( database ) states into a desired target state , within a time polynomial in the size of database and with a probability of close to one . for the fmo complex , we have investigated a variety of different reasonable initial states ( e.g. , an equal superposition of all bchl states , or a localized excitation on bchls 1 and 6 ) and target states ( e.g. , the state localized in bchl 3 ) . in neither case did the coherent hamiltonian dynamics result in a significant overlap to the target state . we found that the overlap oscillates in time but never exceeds the value of 0.4 . we conclude that the _ unitary _ grover search algorithm can not explain the efficiency of the exciton transfer in fmo complex . however , in principle , certain non - unitary generalizations of quantum search algorithms could still be developed to be relevant in this context ; especially for the case of non - unitary oracles interacting with a non - markovian and/or spatially correlated environment . we investigate other non - trivial quantum dynamical effects in the presence of phonon - bath fluctuations and exciton recombination and trapping . the results are shown in figs .  1 , 2 , and 3 . in fig . [ fig1 ] ( a ) , we illustrate the functional dependence of the ete , eq . ( [ ete eq ] ) , on temperature for two initial states localized at site 1 or 6 , respectively . the overall dependence is less than 1% for reasonable temperatures . this can be explained by the relatively small size of the fmo and the approximately three orders of magnitude separation of lifetime ( @xmath162 ) and acceptor transfer ( @xmath163 ) timescales . in order to see this , note that at zero temperature there are only quantum jumps originated from spontaneous emission of energy into the phonon bath , leading to relaxation down the energy funnel . this phenomenon in itself leads to a high efficiency of transport , due to the presence of irreversible trapping on a time scale much faster than the lifetime of the excitation . at higher temperatures quantum jumps due to stimulated emission and absorption enter the dynamics . both processes have the same rates and a temperature dependence which is determined by the bosonic distribution function @xmath164 . in the fmo protein , stimulated emission of excitonic energy helps the transport , since the target state has the lowest energy . the effect of absorption is twofold : it facilitates the overlap with the trapping site when there is an energy barrier in the transport path . at the same time , absorption processes can lead to transfer away from low - energy target sites ( detrapping ) . these effects explain the difference in the temperature dependence for the two initial states localized either at site 1 or at site 6 . the ete increases slightly with temperature for the initial state being at site 1 ( blue line in fig . [ fig1 ] ( a ) ) . this site has a large overlap with exciton 3 and the spatial pathways to site 3 involve energetically higher excitons 4,5,6 , and 7 , see fig .  [ fig3 ] ( a ) . detrapping explains the slight decrease of ete as a function of temperature for initial state 6 ( red line in fig . [ fig1 ] ( a ) ) . this site has a high overlap with exciton 5 and 6 and energetically funnels down to exciton 1 ( site 3 ) . for larger photosynthetic complexes and/or in the absence of low - energy trapping sites , the temperature dependence of the efficiency is expected to be more significant . first , temperature - independent spontaneous emission will play a less prominent role . second , the energy transfer time will become more comparable to the exciton lifetime . this translates to a higher temperature dependence for the overall efficiency since any variation in the transport becomes more pronounced . mes .    the functional dependence of the ete on the reorganization energy at room temperature is shown in fig . [ fig1 ] ( b ) . the reorganization energy can be understood as a linear scaling of the phonon bath . thus , at @xmath165 we have the ( pure ) quantum walk limit leading to efficiencies of 70% and 85% respectively for both pathways . the difference of 15% between the two different initial states can be explained by the higher localization of initial state 6 due larger energy mismatch of site 6 and target site 3**. * * an increase in reorganization energy results in an efficiency increase up to about 99% , demonstrating the effect of the environment - assisted quantum walk . as mentioned before the actual value of the irreversible transfer rate @xmath166 from the fmo complex to the reaction center ( acceptor ) is not well known * * , * * since there exists insufficient chrystallographic data for the combined fmo / rc structure @xcite . in fig . [ fig1 ] ( c ) , we explore the dependence of the ete on this unknown parameter at room temperature**. * * the ete increases monotonically from zero to almost one within a range of five orders of magnitude . one obtains the largest increase in efficiency when the acceptor transfer rate is @xmath39ns@xmath9 thus surpasses the lifetime rate @xmath167 in the limit of large transfer rates to the acceptor , @xmath168 the lifetime of the excitation does not significantly reduce the efficiency . on the other hand , in the limit of small transfer rates to the acceptor , @xmath169 most of the excitation dissipates into the environment before being transferred to the acceptor . figure [ fig1 ] ( d ) , shows the transfer time , eq . ( [ transfertime eq ] ) , of the excitation initially at site 1 or 6 to the acceptor as a function of reorganization energy . in the fully quantum limit it takes the excitation more than 50ps to arrive at the acceptor . this value dramatically improves for higher reorganization energies . at the value of @xmath144cm@xmath9 one finds a transfer time of around 4ps@xmath170 which was reported based on different considerations in @xcite . figure [ fig2 ] illustrates the susceptibilities @xmath171 for the basic processes in the fmo dynamics including the hamiltonian , the phonon bath coupling , transfer to the acceptor , and the loss of the excitation . in fig . [ fig2 ] ( a ) , the susceptibilities for those processes are shown as a function of temperature . the susceptibility of the ete to the phonon bath shows several interesting features . the system is rather susceptible to perturbations of the phonon bath coupling at zero temperature . in this limit , the directionality of the quantum walk is maximized , since , due to spontaneous emission , the excitation can only move down in energy towards the target site 3 while energy absorption from the phonon bath is suppressed . increasing temperature leads to increased stimulated emission _ and _ absorption , thus the system is less susceptible to perturbations of the phonon bath since their effect on emission and absorption processes is the same . concomitantly , the efficiency becomes more susceptible to the transfer process to the reaction center at higher temperature . this is again readily explained by the increased phonon - bath absorption rate : in the presence of detrapping processes it becomes more important to capture the excitation immediately once it arrives at site 3 . we note that the efficiency is not susceptible to a small variation of the free hamiltonian due to inherent irreversibility of the fully unitary quantum dynamics . the susceptibility with respect to dephasing in the energy basis has also insignificant variation for all the temperatures considered here due to the dominant role of the quantum jumps to the overall ete , and thus is not presented in fig . [ fig2 ] . the dependence of the susceptibilities of the ete as a function of the reorganization energy is shown in fig . [ fig2 ] ( b ) . the susceptibility of the ete to the phonon bath peaks at around 1cm@xmath9 for the initial state 1 and around 0.01cm@xmath9 for the initial state 6@xmath172 note that , for the ohmic spectral density considered here , the reorganization energy scales linearly the strength of the coupling to the phonon bath , therefore this curve is the derivative of the ete in fig . [ fig1 ] ( b ) , where the steepest ascent occurs around 1cm@xmath9 or 0.01cm@xmath9 respectively . the susceptibility of the ete on the transfer rate to the acceptor becomes smaller for higher reorganization energy . larger reorganization energy leads to faster thermalization . in this regime the system is more resilient to perturbations of the transfer process to the acceptor and the dissipation to the environment . the dependence of the susceptibility of the ete on the acceptor transfer rate is illustrated in fig . [ fig2 ] ( c ) . we observe a maximum in the susceptibility at around @xmath173ns@xmath9 , which is the regime where the time scales of transfer to the acceptor and the lifetime are comparable , @xmath174 , compare to fig . [ fig1 ] ( c ) . the susceptibilities of the energy transfer efficiency on inter - site quantum jumps is shown in fig . [ fig3 ] ( b ) . more formally , we look at the effect of perturbations of site to site jump terms and the corresponding damping , @xmath175 in the quantum walk master eq . ( transitionsupermatrix ) . we obtain a schematic picture of the most susceptible pathways in the fmo complex . the susceptibilities are correlated with the incoherent transport pathways . the susceptibilities are large when relaxation or absorption of energy due to the phonon bath is important for a particular transport pathway . this is the case e.g.  when there is a relatively large off - resonance between different sites and thus coherent coupling does not lead to significant overlap between the sites , for example for sites 3/7 . on the other hand , when coherent coupling is similar to the site - energy splitting , incoherent processes are less significant , leading to small susceptibilities such as for sites 1/2 and 5/6 . a positive ( negative ) susceptibility of a process @xmath176 indicates an increasing ( decreasing ) ete . for the fmo complex , jumps towards ( away from ) site 3 have positive ( negative ) susceptibility , a signature of the directionality in the irreversible dynamics . thus , the quantum walk formalism together with a rather straightforward measure of the energy transfer susceptibility provides a valuable method to investigate spatial exciton transfer pathways in multichromophoric arrays . we have developed a general theoretical framework within the lindblad formalism for studying the role of quantum effects in energy transfer dynamics of arbitrary chromophoric arrays interacting with a thermal bath from a microscopic hamiltonian . we have shown that a quantum walk approach , which has been widely used in quantum information science , provides an appropriate mathematical framework for studying energy transfer . we have generalized the concept of continuous - time quantum walks to non - unitary and temperature - dependent dynamics in liouville space . this approach can also be used to generally study decoherence effects in quantum walks in arbitrary geometries . the energy transfer efficiency was used as a universal measure to study the transfer properties of the environment - assisted quantum walk . we we have applied our method to explore the energy transfer efficiency and its susceptibilities for the fenna - matthews - olson protein complex as a function of temperature , reorganization energy , trapping rates , and quantum jumps . in particular the energy transfer susceptibilities were studied with respect to the free hamiltonian , the phonon bath , dephasing , trapping rate , and the exciton loss . this approach provides valuable insight into the dynamical role of various generators of the master equation and into spatial exciton transfer pathways in multichromophoric arrays . we have shown that the overall environment coupling strength leads to a substantial enhancement of the ete of about 25% for the fmo complex . thus , the overall energy transfer efficiency of 99% can be explained with the open nature of the multichromophoric dynamics involving an effective interplay between free hamiltonian and fluctuations in the protein and solvent . we applied the general formalism presented in this work to reveal the contributions of underlying physical mechanism to quantum transport @xcite . we also quantified the phenomena of environment - assisted quantum transport due to an effective interplay of quantum dynamical coherence and a pure - dephasing noise model @xcite . using pure dephasing model , others observed similar effects @xcite . motivated by recent observations of non - local effects in the structure and the dynamics of the purple bacteria reaction center @xcite and the fmo protein @xcite , and also based on recent advances in the understanding and control of non - markovian open quantum systems @xcite , generalizations of our scheme to include environments with temporal and spatial correlations are currently underway . these results could lead to new ways for engineering optimal state transfer in quantum spin networks @xcite which interact with realistic environments . since we employ the lindblad equation , our quantum walk formalism can be implemented numerically using the monte carlo wavefunction ( mcwf )  approach . the mcwf method was originally developed for dissipative processes in quantum optics @xcite . the main advantage of mcwf  is the fact that one only needs to simulate the wave function rather than the density operator @xcite . our approach can potentially be used to enhance energy transfer efficiency via engineering quantum interference effects . for certain binary tree chromophoric arrays , e.g. , dendrimers , * *  * * quantum walks could lead to an exponential speed - up over classical walks of excitations @xcite . in general , the combined biology and quantum information inspired approach of this study could provide new insight for engineering artificial photosystems , such as quantum dots and dendrimers @xcite to achieve optimal energy transport by exploiting their environmental effects . we would like to acknowledge useful discussions with j. biamonte , g.r . fleming , i. kassal , a. najmaie , a.t . rezakhani , and l. vogt . we thank the faculty of arts and sciences of harvard university , the army research office ( project w911nf-07 - 1 - 0304 ) and harvard s initiative for quantum science and engineering for funding .
Abstract:

Answer: energy transfer within photosynthetic systems can display quantum effects such as delocalized excitonic transport . recently , direct evidence of long - lived coherence has been experimentally demonstrated for the dynamics of the fenna - matthews - olson ( fmo ) protein complex [ engel _ et al . _ , nature 446 , 782 ( 2007 ) ] . however , the relevance of quantum dynamical processes to the exciton transfer efficiency is to a large extent unknown . here , we develop a theoretical framework for studying the role of quantum interference effects in energy transfer dynamics of molecular arrays interacting with a thermal bath within the lindblad formalism . to this end , we generalize continuous - time quantum walks to non - unitary and temperature - dependent dynamics in liouville space derived from a microscopic hamiltonian . different physical effects of coherence and decoherence processes are explored via a universal measure for the energy transfer efficiency and its susceptibility . in particular , we demonstrate that for the fmo complex an effective interplay between free hamiltonian and thermal fluctuations in the environment leads to a substantial increase in energy transfer efficiency from about 70% to 99% .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: strongly correlated electrons and heavy fermions display a number of phenomena not understandable within simple fermi - liquid theory of metals . these unusual effects include dynamical formation of magnetic local moments and magnetic orders , mixed - valence fluctuations , kondo and heavy - fermion behavior , non - fermi - liquid phases and metal - semiconductor transitions . at present there are two principal approaches for treating electron correlations . ab - initio investigation within the density - functional theory and the local - density approximation takes explicitly into account the coulomb repulsion to equilibrium states within elementary cells ( muffin - tin sphere ) spread over a few interatomic distances . short - range effects of the electron correlations in the ground state are quite well captured . much less reliable are the results of the local - density approximations for long - range effects and low - energy excitations above the ground state . this failure can be attributed to the underestimation of the pauli exclusion principle and long - range fluctuations spreading over thousands of lattice constants . the other approach to the treatment of the electron correlations is based on tight - binding model calculations . the microscopic structure of the electron lattice gas is significantly simplified by taking into account only a few atomic ( wannier ) orbitals . on the other hand , large - scale effects and low - energy excitations are well described . presently the two approaches are more or less complementary . a realistic description of the electron correlations has to join the two existing approaches to enable a sufficiently accurate treatment of both short as well as long distances . however , to succeed in a more reliable description of the electron correlations on small and large scales , it is mandatory to pick up and maintain the relevant aspects of each of the construction . the former supply us with realistic static properties of the ground - state . the latter offers qualitatively accurate estimates of dynamical fluctuations and collective responses to low - energy excitations of strongly correlated systems . it appears that most of the physics of strongly correlated electrons is covered by the archetypal anderson impurity and hubbard lattice models with only a local electron interaction . the weak - coupling regime is governed by a hartree - fock mean field and fermi - liquid theories . extended systems at low temperatures are pauli paramagnets with smeared out local magnetic moments . for bipartite lattices antiferromagnetic long - range order sets in at half filling at arbitrarily small interaction . in the strong - coupling regime the hubbard model at half filling maps onto a heisenberg antiferromagnet with pronounced local magnetic moments and the curie - weiss law for the staggered magnetic susceptibility . the spectral structure is dominated by separated lower and upper hubbard bands and the strongly correlated system seems insulating even in the paramagnetic phase . the most difficult task of the model calculations is to address the transition region between the weak and strong coupling limits . we present in this paper a theoretical approach based on many - body perturbation theory with sophisticated renormalizations of one as well as two - particle green functions ( propagators ) . in particular we show how the complicated algebra of two - particle multiple scatterings can be simplified to a manageable form . the presented single - band theory can be extended to the multi - band case and can hence be used for a more realistic electron - structure calculations . we have at our disposal relatively reliable techniques to describe the two extreme limits of weak and strong couplings of the hubbard model . we can use perturbation expansions in each of the limits . in the former case it is the expansion in the interaction strength while in the latter in the kinetic energy ( hopping matrix ) . the problem is that the two expansions do not match each other and break down before the transition regime is reached . the effective coulomb repulsion becomes comparable with the kinetic energy in the transition region and there is no apparent small parameter at intermediate coupling . dynamical fluctuations control the low - temperature physics of interacting electrons at intermediate coupling and neither weak - coupling nor atomic - like perturbation theories are adequate . weak and strong - coupling perturbation expansions with the bare interaction can no longer reflect the relevant physics with creation and annihilation of long - living electron - hole pairs . the electron - hole pairs start to condense and lead to metal - insulator and magnetic transitions generating poles in appropriate two - particle green functions .    the failure of the existing analytic - numerical methods to describe properly the transition between weak and strong coupling lies in an inadequate treatment of the two - particle criticality . while the one - particle functions are usually renormalized due to a self - consistency in the self - energy , the coupling constant and two - particle vertices remain unrenormalized . however , the two - particle criticality demands a theory with a dynamical renormalization of the relevant two - particle green functions .    to succeed in a quantitative description of this complex situation , it is necessary to reformulate perturbation theory in such a way that the bare interaction be systematically replaced with fully renormalized two - particle functions determined from multiple two - particle scattering processes . it means that we have to use a perturbation theory at the level of two - particle green functions including dynamical vertex renormalizations . the most advanced approximation for two - particle functions are the so - called parquet diagrams . parquet diagrams were introduced to describe interaction of mesons more effectively @xcite . since then a number of attempts have been made to utilize the nontrivial renormalization scheme of the parquet algebra in condensed matter . kondo effect @xcite , x - ray edge problem @xcite , formation of the local magnetic moment @xcite are among the most well known applications . inability to solve the parquet equations effectively has impeded broader application of the method . parquet diagrams represent a systematic way of a summation and renormalization of feynman graphs for two - particle green functions . instead of concentrating on the one - particle irreducible diagrams and the dyson equation , the parquet approach sums two - particle diagrams contributing to vertex functions for which bethe - salpeter equations are constructed . the resulting algebra is much more complicated than that of the one - particle approximations . there are namely three topologically inequivalent definitions of a two - particle irreducibility . it may be defined according to cutting an electron - hole or electron - electron pair propagation , or according to cutting a polarization bubble shielding the electron - electron interaction . each possibility defines a two - particle channel of multiple pair scatterings characterized by a different binding of independent variables in the vertex functions . a general two - particle quantity will be denoted in this paper as in fig .  1 . each two - particle function carries three independent four - momenta and two spin indices . we generally denote the fermionic four - momenta @xmath0 and the transferred bosonic ones as @xmath1 , where @xmath2 and @xmath3 are the respective matsubara frequencies at temperature @xmath4 . it is convenient to introduce a matrix notation in the spin indices to distinguish different two - particle channels in the parquet diagrams . we define a @xmath5 matrix for the generic two - particle function @xmath6 : @xmath7 we speak about singlet and triplet contributions to a two - particle function if the spins of the involved fermions are antiparallel or parallel , respectively . we define three matrix multiplication schemes for the two - particle quantities    @xmath8_{\sigma\sigma'}(k , k';q)&= &    \frac 1{\beta{\cal n}}\sum_{q '' } x_{\sigma\sigma'}(k , k';q'')y_{\sigma      \sigma'}(k+q'',k'+q'';q - q''),\\   \label{eq : conv - ee }    \left[\widehat{x}\circ\widehat{y}\right]_{\sigma\sigma'}(k , k';q)&=&\frac    1{\beta{\cal n}}\sum_{q '' } x_{\sigma\sigma'}(k , k'+q'';q - q'')y_{\sigma      \sigma'}(k+q - q'',k';q''),\\    \label{eq : conv - u }    \left[\widehat{x}\star\widehat{y}\right]_{\sigma\sigma'}(k , k';q)&=&\frac    1{\beta{\cal n}}\sum_{\sigma''k''}x_{\sigma\sigma''}(k , k'';q)y_{\sigma ''      \sigma'}(k'',k';q ) \end{aligned}\ ] ]    representing summations over the intermediate states in the three inequivalent two - particle channels , electron - hole , ( @xmath9 ) , electron - electron , ( @xmath10 ) , and interaction , ( @xmath11 ) , respectively . we see that the variables of the two - particle functions are convoluted differently in inequivalent channels . note that only the interaction channel mixes the singlet and triplet contributions . we decompose the full two - particle green function into a sum of always reducible and irreducible projections onto each inequivalent channel @xmath12 where @xmath13 refers to a two - particle channel . the parquet diagrams can at best be represented graphically . having in mind the above introduced notation and a general convention that double primed variables are summed over , we can write    @xmath14    equations ( 4 ) are the bethe - salpeter equations and represent one set of the full parquet approximation . to complete it we must add relations connecting the two - particle reducible , @xmath15 , and irreducible , @xmath16 , functions . for this purpose we introduced two - particle bubbles for each channel    @xmath17 \label{eq : g2v }    g^{(2)u}_{\sigma\sigma'}(k , k';q)&=&g_{\sigma'}(k')g_{\sigma ' }   ( k'+q)\ .\end{aligned}\ ] ]    the functions @xmath16 are then defined as @xmath18g^{(2)\alpha}_{\sigma\sigma ' }    ( k , k';q )    \end{aligned}\ ] ] where @xmath19 is a sum of all feynman diagrams simultaneously irreducible in each two - particle channel . the sum contains only higher - order diagrams where three and more particles are multiply interconnected . in the usual treatment with only two - particle multiple scatterings this irreducible part is neglected . we hence put @xmath20 . parquet equations ( 4 ) define the two - particle functions as a functional of the renormalized one - particle propagators containing the self - energy . the self - energy is a functional of the full two - particle function @xmath21 from the parquet equations . hence to complete the parquet approximation we have to add formulas determining the dependence of the one - particle propagator on the self - energy ( dyson equation ) and the self - energy on the full two - particle function . the consistency between one and two - particle quantities demands @xmath22    ^{-1}\ , \\[4pt ]   \label{eq : sigma }    \sigma_\sigma(k)&=&-\frac u{\beta^2{\cal n}^2}\sum_{\sigma'}\sum_{k'q }    g_{\sigma}(k+q){\cal k}_{\sigma \sigma'}(k+q , k';-q)g_{\sigma'}(k ' )    g_{\sigma'}(k'+q )    \end{aligned}\ ] ] where @xmath23 is the particle density , @xmath24 , and @xmath25 is the chemical potential and @xmath26 an external magnetic field . in our notation measures corrections to the static hartree approximation contained in the particle densities @xmath23 determined from a sum rule @xmath27 . ] equations ( [ eq : k - def])-([eq : sigma ] ) constitute the parquet approximation . in this approximation only fully renormalized one and two - particle functions appear . it means that not only all one - particle propagators are fully renormalized with the self - energy insertions but also the bare electron - electron interaction in the two - particle scatterings is replaced by corresponding irreducible two - particle ( vertex ) functions . what is missing in the parquet diagrams are three - particle and higher - order cumulant functions . they are , however , not expected to be relevant for the static coulomb interaction . they contribute via dipole and higher - order multi - pole interactions and are weak in the hubbard model . parquet diagrams differ significantly from simpler , single - channel approximations using multiple two - particle scatterings such as @xmath28 , @xmath29 or @xmath30-matrix ( tma ) . the difference manifest itself in the effective interaction in the two - particle scattering processes . the bare interaction in the two - particle scattering events remains unrenormalized in the single - channel approximations . in the parquet equations it is replaced by renormalized vertex functions representing a dynamically screened interparticle potential . a qualitative difference between a solution to the full parquet equations and to the single - channel approximations becomes evident at the two - particle criticality . the vertex functions from one or more channels get critical and sharply peaked around the fermi energy . the exchange between the scattered particles becomes strongly energy dependent and can not be approximated by a static effective interaction . moreover , the singularity at the critical point in the parquet approximation must be integrable , while the single - channel approximations usually display logarithmic divergences in the integrated quantities ( e.  g. self - energy ) . on the other hand the nonlinear integral equations represented by the parquet diagrams are extremely complicated with an intricate analytic structure of the solution . each two - particle function entering the parquet equations contains three independent ( complex ) variables . neither of them can simply be neglected , since the variables are mutually interconnected due to the scatterings in different inequivalent channels . consequently , a nonperturbative solution to the parquet equations has not yet been found . that is also why the parquet diagrams are relatively little used in condensed matter although being in existence more than thirty years . the only hope that the parquet approximation may become useful for quantitative calculations of the effects of the electron correlations is that one succeeds in simplifying the full parquet algebra to a manageable form . to this purpose we have to decide where we want the parquet approximation to be reliable . as discussed in the introduction it is the intermediate coupling with critical two - particle functions where we need to improve upon easy - minded approximations . we hence demand asymptotic reliability of the parquet diagrams at the two - particle critical points and in particular at the metal - insulator transition . we retain only the diagrams making the effective mass of the electrons divergent at the metal - insulator transition . it was recently argued that at least at half filling only two of the three inequivalent channels are relevant in the critical region of the metal - insulator transition @xcite . only multiple scatterings in the electron - hole and interaction channels contribute to divergent diagrams at the critical point . this simplification , called dipole approximation , however , does not reduce the number of the variables in the two - particle functions . it only simplifies the algebra of the parquet approximation . to lower the number of the relevant variables we utilize the fact that the singularities in the two - particle functions in the parquet approximation must be integrable . it means that when integrating over the variable in which the singularity arises we obtain a finite result . our approximation consists in neglecting all finite ( nondiverging ) contributions at the critical point and keeping only the divergent ones . this must be done at the level of the diverging two - particle functions . there is always one relevant bosonic variable in each two - particle channel . it is the variable remaining untouched in the scattering processes . it can easily be seen from multiplication rule ( 2 ) that it is @xmath31 , @xmath32 , and @xmath33 for the electron - hole , electron - electron , and interaction channels , respectively . divergences in the reducible functions appear in the respective conserved variables . for the metal - insulator transition only @xmath31 and @xmath33 are relevant .    in the first step we reduce the number of the relevant variables in the two - particle functions to two . this is achieved when we neglect mixing of the reducible functions from different channels . we replace the full two - particle function @xmath21 on the right - hand side of ( 4 ) by @xmath34 , i.  e. by the sum of the completely irreducible vertex and the reducible function in the @xmath35-channel . such a replacement does not change the ( universal ) critical behavior of the resultant reducible function @xmath36 . the bare interaction is a constant and consequently the solution does not depend on the outgoing variables . further on , the triplet functions can be eliminated from this approximation , since @xmath37 gets irrelevant . if we denote @xmath38 and @xmath39 , the parquet equations reduce to @xmath40    \label{eq : parq2b }     i^h_{\sigma,-\sigma}(k;q)&=&u+\frac 1{\beta^2{\cal n}^2}\sum_{q',q '' }    i^v_{\sigma,-\sigma}(k;q')g_{-\sigma}(k+q ' ) g_{-\sigma}(k+q'+q )    \nonumber\\   & & \hspace*{-10pt } \times i^v_{\sigma,-\sigma}(k+q';q''-q ' ) g_\sigma   ( k+q '' ) g_\sigma ( k+q''+q)i^h_{\sigma,-\sigma}(k+q'';q)\ . \end{aligned}\ ] ] the bosonic variable @xmath33 has different meaning in the functions @xmath41 and @xmath42 . it is the horizontal , vertical transfer momentum for the former and the latter , respectively .    because of breaking of the symmetry between the incoming and outgoing variables in the two - particle functions we must be careful in deriving the equation for the self - energy . to be consistent we use a symmetrized formula with the bare interaction at the incoming and outgoing momentum . we obtain @xmath43\ .\end{aligned}\ ] ]    the reduced two - channel parquet approximation with two - variable two - particle functions has the critical behavior of the full parquet approximation . the universal quantities derived from the divergent functions of the parquet equations ( 4 ) are reproduced in equations ( [ eq : parq2a])-([eq : sigma2 ] ) . only nonuniversal quantities such as the critical interaction do not come out in this simplification as in the full parquet approximation . we achieved a simplification of the full parquet algebra without loosing leading critical behavior of the solution , but the resulting equations have nonlinear convolutive character . we have to find a nonperturbative solution to these equations when addressing the metal - insulator transition . this is , however , still beyond the reach of our analytic - numerical tools . we must find a further reduction of complexity of the parquet equations to make them useful for a quantitative analysis . it is the bosonic variable @xmath33 in the two - particle functions @xmath44 that is relevant at the critical point and in which the singularity arises . the fermionic variable @xmath45 in both functions labels the eigenvalues of the integral kernel in equations ( [ eq : parq2a ] ) , ( [ eq : parq2b ] ) . the maximal eigenvalue governs the critical behavior . if we are interested only in the critical asymptotic behavior we can approximate the actual maximal eigenvalue by a value at the lowest - lying fermionic four - momentum . by such an approximation we neglect mixing of the fermionic four - momenta in the equations and replace @xmath45-dependence in the two - particle functions by a constant , their value at @xmath46 . this approach corresponds to a low - energy expansion used by hamann @xcite by assessing the strong - coupling limit of the renormalized rpa of suhl in the single - impurity problem . the suggested simplification is possible only at zero temperature , where the difference between the fermionic and bosonic matsubara frequencies vanishes . we hence put @xmath47 , which is the most interesting case for the metal - insulator transition . using the above ansatz , neglecting @xmath45-dependence in @xmath44 , and putting @xmath46 in the one - electron propagators we reduce the parquet equations to a manageable form . we introduce new functions @xmath48 with the aid of which we obtain for the two - particle functions @xmath49    \label{eq : i^h } i^h_{\sigma,-\sigma}(q)&= & \frac u{1-\mathcal{k}_{-\sigma}(q )    \mathcal{k}_{\sigma}(q ) } \ .\end{aligned}\ ] ] the above equations have the complexity of the single - channel approximations , but contain the full vertex renormalization in the two - particle scatterings . it is always the renormalized vertex function that enters the denominator of equations ( [ eq : i^v ] ) , ( [ eq : i^h ] ) via the functions @xmath50 and @xmath21 . the single channel approximations replace the two - particle vertex functions @xmath44 in ( [ eq : gam ] ) , ( [ eq : k ] ) by the bare interaction @xmath11 . it is easy to derive the corresponding equation for the self - energy . it reads @xmath51\end{aligned}\ ] ] where we introduced a single bubble @xmath52    equations ( [ eq : gam])-([eq : x - def ] ) replace the original parquet approximation . they fully determine the generating two- and one - particle functions . all thermodynamic and spectral quantities can be calculated from them . since we used a number of approximate steps in the derivation of the reduced parquet algebra , any new quantity must first be formulated within the full parquet approximation . the simplifications are then introduced via the generating one- and two - particle functions . in this way a consistent theory is built up , since the parquet approximation can be viewed upon as a conserving approximation derived from a generating functional @xcite . an exact solution as well as the full parquet approximation fulfill a number of identities and symmetries . one of the most important features of the exact solution is the electron - hole symmetry . it means , that the physics must not depend on whether we use the particle or the hole picture . the electron - hole symmetry is obeyed by the complete parquet approximation . it is not clear whether we have not lost some of the required symmetries in course of simplifying the parquet algebra . we lost the symmetry between the incoming and outgoing variables in our approximate equations . further on we also lost the electron - hole symmetry in only one particle of the pair in the two - particle scatterings . the former was lost during the reduction to two - variable two - particle functions and the latter due to the suppression of the electron - electron channel . these losts are not so serious if we are in the critical region of the metal - insulator transition . the one - sided electron - hole transformation changes the sign of the effective interaction and hence breaks the bound electron - hole pair . the full electron - hole symmetry in all variables must hold for any approximate theory , otherwise we were unable to guarantee the fermi liquid properties at weak coupling . the electron - hole transformation may be defined in the parameter space of the hubbard model as @xmath53 where the bar denotes a hole representation . it is easy to check by inspection that a solution to the reduced parquet equations ( [ eq : gam])-([eq : x - def ] ) fulfill the following symmetry relations @xmath54   \label{eq : eh - sym2 } \gamma_{\sigma}(q)&=&\ \bar{\gamma}_{\sigma}(-q)\ ,   \hspace{10pt}\mathcal{k}_{\sigma}(q)=\ \bar{\cal k}_{\sigma}(-q ) \ . \end{aligned}\ ] ]    this electron - hole symmetry together with the analyticity of the one- and two - electron functions is used to prove fermi - liquid properties of solutions to ( [ eq : gam])-([eq : x - def ] ) for sufficiently weak interaction , i.  e. before poles in the vertex functions @xmath55 appear . up to know we have worked with a very general formulation of the parquet approximation applicable in any spatial dimension . however , recent studies indicate that most of the effects of strong electron correlations can qualitatively be understood already within a dynamical mean - field theory or impurity models @xcite . it is sufficient in many situations to reduce the effects of strong correlations to a dynamical mean - field where the self - energy and vertex functions are local , but frequency - dependent . such an approach significantly reduces the numerical complexity but keeps much of the interesting physics of strongly correlated electrons . we apply the simplified parquet approximation to a single impurity or a dynamical mean field model . the fluctuations in space do not contribute to the dynamics and the four - momenta from the general formulation of the parquet approximation collapse to frequencies . this simplification enables one to perform explicitly analytic continuation from the matsubara to the real frequencies using contour integrals . realizing that the reduced equations hold for @xmath47 and introducing @xmath56 for the bosonic , fermionic complex frequencies , respectively , we can write for the two - particle functions @xmath57,\ ] ] @xmath58   & & \hspace*{80pt}\left . + \frac{g_\sigma(\omega-\zeta)}{1-{\cal k}_\sigma(\omega-\zeta){\cal         k}_{-\sigma}(\omega-\zeta ) } \mbox{im}g_{-\sigma}(\omega_+)\right ] \end{aligned}\ ] ] and analogously for the one - particle self - energy @xmath59\right.\nonumber\\    & & \hspace*{-50pt}\left . + x_{-\sigma,-\sigma } ( \omega - z)\left(i^h_{-\sigma ,      \sigma}(\omega_+ ) -u\right)\mbox{im}g_\sigma(\omega_+ )    + g_{-\sigma}(\omega+z)\mbox{im}\left[x_{\sigma,-\sigma}(\omega_+ )       i^v_{\sigma,-\sigma}(\omega_+)\right ] \right.\nonumber \\     & & \left . + x_{\sigma,-\sigma}(\omega - z)i^v_{\sigma,-\sigma}(\omega - z )      \mbox{im}g_{-\sigma}(\omega_+)\right\ } \ \end{aligned}\ ] ] where @xmath60 . the local propagator is for the mean field @xmath61 and for the impurity model @xmath62^{-1 } \end{aligned}\ ] ] where @xmath63 for the mean - field and @xmath64 for the impurity model . in the former @xmath65 is the density of states and in the latter , the parameter @xmath66 is the energy of the impurity orbital , @xmath67 hybridization to the conduction electrons , and @xmath68 the one - electron local green function of the conduction electrons . note that there is no difference at the two - particle level between the impurity and mean - field equations . we use the local formulation of the simplified parquet approximation to demonstrate that intermediate coupling is dominated by a two - particle criticality . we confine ourselves only to half filling and the spin symmetric case ( paramagnet ) where the singularity indicates a metal - insulator transition . we further assume a bipartite lattice , i.  e. @xmath69 . a solution to equations ( [ eq : k - cont])-([eq : sigma1 ] ) has at weak coupling the required fermi - liquid properties . the two - particle functions @xmath70 are real at the fermi energy . iterations starting with @xmath71 converge to a unique solution . such a procedure becomes numerically unstable if we increase the interaction beyond @xmath72 defined as @xmath73^{-1 } , \quad   x_0(0)=-\int_{-\infty}^0    \frac{d\omega } \pi\ \mbox{im } \left[g^{(0)}(\omega_+)^2\right ] \end{aligned}\ ] ] where @xmath74 is a noninteracting one - particle propagator . the value @xmath72 corresponds to the critical interaction determined from the @xmath28 pole in the electron - hole susceptibility . explicitly , @xmath75 for a mean - field ( @xmath76 ) model with a lorentzian dos ( infinite bandwidth ) and @xmath77 for a bethe @xmath76 lattice ( finite bandwidth ) . this interaction need not necessarily represent a critical point . it only indicates that we can not start iterations from an unperturbed solution and that weak - coupling perturbation theories get numerically unstable . we can say that the weak - coupling regime goes over at @xmath72 to an intermediate one . the transition is , however , smooth for effective impurity models .    at intermediate coupling the functions @xmath42 and @xmath41 approach a pole at the fermi energy @xmath78 . it means that @xmath79 as well as @xmath80 approach zero . a new small ( dimensionless ) scale @xmath81 arises . this scale can not be derived from the one - electron model parameters . it is a consequence of a two - particle criticality . in the single - impurity case it is related to the kondo scale . the scale vanishes at the metal - insulator transition if it is of second order . the vanishing scale @xmath82 causes nonanalyticities in the two - particle functions @xmath83 and @xmath84 , and in the one - electron self - energy @xmath85 . since the singularity in the vertex functions @xmath86 must be integrable , the functions @xmath87 remain bounded ( continuous ) . only their derivatives show divergences at the fermi energy at the critical point . it is straightforward to show that leading - order divergences in ( [ eq : k - cont])-([eq : sigma1 ] ) are @xmath88 ^ 2}{1 + \gamma(0)}\ , \hspace{10pt}\mbox{im}\      \gamma'(0_+)\doteq -\frac u\pi\ \frac{\left [        \mbox{im}\ g(0_+)\right]^2}{1-\mathcal{k}(0)^2}\ , \\[4pt ]    \label{eq : sigma - sing }     \mbox{re}\ \sigma'(0)&\doteq&-\frac u\pi ux(0)\ \mbox{im}\        g(0_+)\left[\frac 1{1         + \gamma(0)}+\frac 1{1-\mathcal{k}(0)^2}\right ] \end{aligned}\ ] ] where the prime denotes the derivative w.r.t . the frequency variable . we utilized the symmetry of the half - filled case , @xmath89 in the derivation of ( [ eq : k - gam - sing ] ) , ( [ eq : sigma - sing ] ) . the mott - hubbard metal - insulator transition is defined by divergence of the effective mass @xcite . since @xmath90 we have proved that the singularity in the two - particle vertex functions @xmath55 causes divergence of the effective electron mass and hence the mott - hubbard metal - insulator transition . the way the effective electron mass diverges is essentially determined by vanishing of the kondo scale @xmath82 . it means that the kondo scale sets an interval around the fermi energy within which the landau quasiparticle picture holds not only for the one - electron but also for the two - electron functions . the narrow quasiparticle kondo - suhl peak is then a consequence of sharply peaked two - particle vertex functions . no effective static interaction can describe the kondo strong - coupling asymptotics and the mott - hubbard metal - insulator transition with ( [ eq : k - gam - sing ] ) , ( [ eq : sigma - sing ] ) reliably . the simplified parquet approximation describes qualitatively well the weak - coupling fermi - liquid regime . it also correctly binds divergence of the effective mass with the kondo scale arising from a singularity in two - particle functions . the derived approximate theory does well on the metallic side of the metal - insulator transition . whether the approximation can correctly capture the physics up to the transition point must be decided from the dependence of the kondo scale @xmath82 onto the interaction strength @xmath11 or better on the parameter @xmath91 . the correct behavior must reproduce the exponential law @xmath92 for the single - impurity model . we need a deeper analysis of equations ( [ eq : k - cont])-([eq : sigma1 ] ) to decide whether this really is the case . it will be done in a separate paper . next we have to address the strong - coupling regime for @xmath93 . fermi - liquid theory no longer holds and it is unclear how to iterate towards a stable solution . the problem is to choose the value of the self - energy at the fermi level . it vanishes in the fermi - liquid , metallic phase but diverges in the atomic , insulating solution . what value the self - energy acquires in a self - consistent insulating solution to ( [ eq : k - cont])-([eq : sigma1 ] ) needs to be investigated in more details . one can , however , easily check that the simplified parquet approximation allows for a consistent @xmath94 limit . after an appropriate scaling of the variables one recovers an atomic solution at zero temperature . dynamical fluctuations then enter as perturbations in @xmath95 , where @xmath96 is an effective bandwidth of the conduction electrons . a perturbation expansion in @xmath95 can be a good starting point for iterations to the full solution of equations ( [ eq : k - cont])-([eq : sigma1 ] ) in the strong - coupling , insulating phase . we demonstrated in this paper that a successful description of intermediate and strong coupling regimes in correlated electron systems demands sophisticated approximations with a dynamical renormalization of two - particle vertex functions . a transition regime between the weak and strong coupling limits is governed by the formation and condensation of long - living electron - hole pairs resulting in poles in two - particle green functions . these dynamical effects can not properly be rendered by theories with only effective static interactions . the well - known single - channel two - particle approximations such as rpa , tma , or gwa are inadequate at intermediate coupling , since the interaction in the multiple scatterings remains unrenormalized . only the full parquet approximation with all two - particle irreducible diagrams can be expected to capture the relevant physics at the two - particle criticality of intermediate coupling . the parquet diagrams lead to a qualitatively different critical behavior than that obtained from the simple - minded single - channel approximations . the only singularities that can appear in the parquet approximation are integrable ones , i.  e. convolutions of singular and regular functions must be finite . this feature distinguishes the parquet theory from the other existing approximations in correlated electron systems . the parquet algebra is too intricate to allow for a nonperturbative solution . only an appropriate reduction in the complexity of the parquet approximation can lead to a successful application of the method . we proposed in this paper a scheme how to reduce the nonlinear convolutive integral parquet equations to a system of quasi - algebraic equations not more complicated than the single - channel approximations . the only difference is that the effective interaction in the scattering events in one channel is dynamically renormalized by scatterings from other channels . the reduction of the parquet algebra is enabled by taking into account only leading - order divergences at the metal - insulator transition and a low - energy expansion around the maximal eigenvalue of the parquet equations . the resulting equations allow only for integrable singularities and lead to sharply peaked two - particle vertex functions in the vicinity of the critical metal - insulator transition . they correctly bind divergence of the effective one - electron mass with the kondo scale stemming from a two - particle singularity . what still remains to do is to show how well the simplified parquet approximation reproduces the strong - coupling asymptotics of the kondo scale in the single - impurity problem known from an exact solution . only then we can be sure we have a reliable interpolation between the weak and strong - coupling regimes of hubbard - like models . this work was supported in part by the grant no . 202/98/1290 of the grant agency of the czech republic . 99 v.  v.  sudakov , dokl . nauk sssr * 111 * , 338 ( 1956 ) , [ english transl . soviet phys . - doklady * 1 * , 662 ( 1957 ) ] . + a.  a.  abrikosov , physics * 2 * , 5 ( 1964 ) . + b.  roulet , j.  gavoret , and p. nozires , phys . rev . * 178 * , 1078 ( 1969 ) . + r.  a. weiner , phys . rev . lett . * 24 * , 1071 ( 1970 ) . + v.  jani , j. phys . : matter * 10 * , 2915 ( 1998 ) . + d.  r.  hamann , phys . rev . * 186 * , 549 ( 1969 ) . + v.  jani , preprint _ cond - mat/9806118 _ + a.  georges , g.  kotliar , w.  krauth , and m.  rozenberg , rev . phys . * 68 * , 13 ( 1996 ) . + w.  f.  brinkman and t.  m.  rice , phys . rev . b*2 * , 4302 ( 1970 ) .
Abstract:

Answer: single - band hubbard model at criticality of the metal - insulator transition is studied using approximations derived from parquet theory . it is argued that only the electron - hole and interaction two - particle channels in the parquet algebra are relevant . a scheme is proposed how to reduce the parquet equations to a manageable form the complexity of which is comparable with single - channel approximations such as the renormalized rpa . the newly derived approximation , however , contains dynamical vertex corrections , remains self - consistent at the two - particle level and allows only for integrable singularities . a qualitatively new approach for studying two - particle singularities at zero temperature is obtained .    * international workshop on electron correlations and material properties , crete greece , june 28 - july 3 , 1998 * + _ to appear in the conference proceedings , eds . n kioussis and a. gonis , plenum press _    * strongly correlated electrons : + dynamical vertex renormalization * +    = vaclav jani + institute of physics , academy of sciences of the czech republic , + na slovance 2 , cz-180 40 praha , czech republic +