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: in recent years , the annihilation - diffusion problem has generated significant interest , both theoretical and experimental . the annihilation - diffusion problem usually corresponds to that of the kinetics of particle density decay in the annihilation reaction @xmath1 ( one - species annihilation ) and @xmath2 ( two - species annihilation ) . in the latter case , in most physical systems that one is interested in , in addition to the thermal diffusion and kinetic annihilation , the particles a and b are charged and interact via a power - law long - range interaction ( lri ) , usually of coulomb type , although here we will consider a more general lri . such physically important interaction can clearly strongly influence the annihilation dynamics by introducing an additional time scale in the annihilation process , and can lead to a new mechanism for slow dynamics . the two - species annihilation reaction with lri can be studied in a variety of experiments . a thermal quench of a freely - suspended liquid crystal film from smectic - a to smectic - c phase @xcite is one experimental system where this annihilation process has been studied in great detail . in such experiments , immediately after quench , the singularities of the smectic director ( 2d vector ) field appear as positive and negative vortices interacting ( due to elastic forces ) via a logarithmic potential . as time elapses , vortices of opposite sign slowly annihilate , exhibiting complex dynamics that is clearly and strongly influenced by temperature , initial particle distribution , and lri that leads to an attraction between annihilating partners . similar annihilation problems also appear in turbulent flow , superconductivity , spinodal decomposition , and many other condensed matter systems . the annihilation process is also relevant to coarsening of topological defects produced by symmetry breaking field in particle physics models , after an early temperature quench due to the fast initial universe expansion , a process that is thought in part to determine the large scale structure of today s universe.@xcite    it is well known that in classical chemical kinetics , density decay for both one - species and two - species annihilation is described by the kinetic rate - equation ( see , e.g. , @xcite ): @xmath3 with the large-@xmath4 asymptotics given by : @xmath5 where @xmath6 is a reaction constant . equation ( [ eq : kinrate ] ) completely neglects all spatial fluctuations and correlations of the particle density . although the dynamics described by equations  ( [ eq : kinrate ] ) - ( [ eq : mfsoln ] ) is correct above an upper critical space dimension @xmath7 , for @xmath8 , both diffusion and fluctuations play an important role , substantially slowing down the density decay . it was indeed shown  @xcite , that for the one - species annihilation without lri , the @xmath9 . that is , for systems of dimensionality less than @xmath10 , the density decay is given by : @xmath11 where @xmath12 is diffusion constant and @xmath13 is the initial particle density . equation  ( [ eq : difsolnaa ] ) can be understood either by invoking a single length scaling arguments  @xcite or by using a more elaborate smoluchowski approach  @xcite . for the two - species case , it was shown ( refs.@xcite , @xcite ) that the upper critical dimension @xmath14 , and that for @xmath15 , the large-@xmath4 asymptotics is : @xmath16 the decay law  ( [ eq : difsolnab ] ) was confirmed in several numerical simulations of one- , two- , and three - dimensional systems  @xcite . it is important to emphasize that in order to observe such power - law decay , it is necessary initially to have an equal number of positive and negative charges , distributed at random . if initial numbers of positive and negative charges are different , one should observe an exponential density decay to the non - vacuum equilibrium ( see , e.g. , sokolov et al . if the system is well - stirred , i.e. , long - wavelength fluctuations of charge density are suppressed , the decay law would also be different . we will not consider such cases in this paper . to better study the role of correlations in reactions without lri , description in terms of secondary quantization operators of creation and annihilation was proposed . @xcite further developing this approach , peliti  @xcite proposed an `` exact '' ( valid to all orders in perturbation expansion ) renormalization - group theory for the one - species annihilation , and lee and cardy  @xcite suggested a similar approach to the two - species annihilation , both based on the rigorous master equation converted into a field - theoretic formulation . in the latter case , it was shown that the decay law  ( [ eq : difsolnab ] ) is an asymptotic one for large @xmath4 and small @xmath13 . it was also rigorously proven that the assumptions used in the derivation of equation  ( [ eq : difsolnab ] ) are correct only for @xmath17 , while for @xmath18 , a more elaborate renormalization group procedure which includes an effective noise with nontrivial correlations should be carried out . the addition of the lri further complicates the picture . equations of motion become strongly coupled via the nonlocal interaction , and only approximate solutions of these equations can be found . up to date , only coulombic  @xcite systems in two dimensions ( @xmath19 , the case corresponding to particles interacting via a logarithmic potential on a two - dimensional substrate or film ) were studied numerically  @xcite , and results appear to be inconclusive . in these simulations , particle density exhibited a power - law decay @xmath20 with the exponent @xmath21 varying from @xmath22  @xcite to @xmath23  @xcite . the authors proposed a scaling theory  @xcite suggesting that for a coulombic ( @xmath24 ) two - dimensional diffusion - annihilation system , annihilation exponent is equal to @xmath25 - close to but less than the mean - field exponent @xmath26 . oshanin et al . @xcite suggested that the exponent should be exactly 1 ; ispolatov and krapivsky  @xcite came to the same conclusion using an `` inpenetrable domain '' scaling theory . in their theoretical and computational study of the defect annihilation in two - dimensional xy - model , yurke et al . @xcite argued that annihilation exponent is 1 , with logarithmic correction due to the logarithmic dependence of the mobility on the defect size . thus , the problem of annihilation behavior for coulombic system in two dimensions has not been completely resolved , although it seems rather plausible that the final asymptotics is governed by a classical exponent 1 , with possible logarithmic correction . it seems to be certain , however , that annihilation in three - dimensional coulombic system has a mean - field type final asymptotics with exponent 1 . finally , we should point out that until recently , there has been no discussion of lri other than coulombic . recently . several scaling theories have been suggested to analyze arbitrary power - law interaction ( more short - ranged than coulombic )  @xcite.such problems may arise , e.g. , in describing interactions between vacancies and interstitials in a two- or three - dimensional crystal .    in this paper , we propose a self - consistent theory based on coarse - grained hydrodynamic equations of motion for particle number density and charge density fields . this theory allows us to systematically calculate the dependence of the density decay law on initial conditions and to investigate the role of lri . it can be shown@xcite that the self - consistent approximation that we employ here is equivalent to a resummation of an infinite class of feynman graphs which take into account the conserved charge density fluctuations but ignore the less important nonconserved number density fluctuations . this approximation can be further systematically improved by the use of perturbation theory and the renormalization group analysis that will be a subject of future publication.@xcite in the limit of weak long - range interactions , this approach agrees well with the known theoretical results ( see refs . @xcite ) for the two - species annihilation @xmath27 , thereby further clarifying the underlying assumptions that led to these results . the paper is organized as follows : in section  [ sec : eqns ] we define all variables and present the equations of motion on which our analysis and results will be based . in section  [ sec : solns ] the self - consistent approximation is described and the asymptotic solutions are obtained and analyzed . finally , in section  [ sec : concl ] we analyze the resulting phase diagram and discuss it in the context of the previously obtained results for various diffusion - annihilation systems , which are selected points on our phase diagram . let us consider a system consisting of two kinds of particles , a and b , with a having a positive charge @xmath28 and b having a negative charge @xmath29 . we label their , time- and position - dependent concentrations as @xmath30 and @xmath31 , respectively , and impose the condition that @xmath32 . the equation of motion for the densities is based on the generalized law of mass conservation , violated by the annihilation process , @xmath33 where the mass current is given by @xmath34 and @xmath35 is the electrostatic long - range potential at a point * r * at time @xmath4 due to local charge fluctuations away from neutrality ; @xmath36 is the power exponent of the long - range force , @xmath37 is particle mobility taken to be a constant , and @xmath38 is particle charge . equations  ( [ eq : evol ] ) - ( [ eq : poten ] ) should be solved in conjunction with initial conditions for @xmath39 and @xmath40 . since here we are interested in statistical averages , rather than in a dynamic solution for a given system with specific initial conditions , we will focus on the density correlation functions , with averages over random @xmath41 initial conditions . equations  ( [ eq : evol ] ) - ( [ eq : poten ] ) represent the coarse - grained continuum limit for the `` real '' equations of motion ; all variables in these equations are averaged over `` elementary volume '' @xmath42 . thus , only long - wavelength modes are actually described by these equations , and , therefore , only intermediate - time and large - time regimes can be analyzed . it is more natural and convenient to describe the system in terms of the particle number and charge densities . we denote the former one as @xmath43 , and the latter one as @xmath44 and relate them to densities @xmath45 and @xmath46 as follows : @xmath47 @xmath48 rewriting equations  ( [ eq : evol ] ) - ( [ eq : poten ] ) using densities @xmath49 and @xmath50 , we obtain : @xmath51 where @xmath52 .    in the absence of lri , equations  ( [ eq : rhoevol ] ) - ( [ eq : fevol ] ) are those analyzed in ref . however , the presence of the additional long - range interactions , makes the analysis of their asymptotic solutions rather nontrivial and , as we will show below , leads to new dynamic regimes . it is important to note that equations  ( [ eq : rhoevol ] ) - ( [ eq : fevol ] ) do not contain noise terms in their right - hand sides . it has been rigorously shown that such noise terms represent important correlations and in some cases may even become predominant in determining the asymptotic decay rate . in equations of motion describing a near - equilibrium dynamics , powerful fluctuation - dissipation theorem determines the form of noise correlations . in contrast , in systems far from equilibrium , such as a system of annihilating particles , it can be shown@xcite that the effective hydrodynamic equations of motion derived from the fundametal master equations contain noise terms with very nontrivial correlations , of the form which could not be easily guessed a priori . lee and cardy @xcite proved that in a two - species reaction without lri , such noise leads only to the renormalization of the reaction rate @xmath6 , but not to the change of the scaling exponents , provided that space dimensionality @xmath53 . we have shown , in a similar fashion @xcite , that for systems with lri , the noise has no effect on the asymptotic dynamics for @xmath53 , if the renormalization of both @xmath6 and @xmath54 is implied . in addition , since equations  ( [ eq : rhoevol ] ) - ( [ eq : fevol ] ) provide a coarse - grained description on the length scale larger than the interparticle spacing @xmath55 , the kinetic coefficients , e.g. @xmath6 , are effective coefficients that incorporate finite renormalization due to the correlations on short length scales .    in order to simplify further analysis , we divide each of the equations  ( [ eq : rhoevol ] ) -([eq : fevol ] ) by @xmath13 and transform everything to dimensionless variables as follows :    / _ 0 , f f / _ 0 , d d ( _ 0)^2/d ,     + @xmath56    it is important to notice that equation  ( [ eq : fevol ] ) is linear with respect to @xmath49 , while equation  ( [ eq : rhoevol ] ) is quadratic with respect to @xmath49 ( this points to the system s invariance with respect to the simultaneous charge sign reversal for all particles ) . in the next section , we will describe the self - consistent approximation and its solutions . we now make an important approximation in order to further simplify the analytical treatment of equations  ( [ eq : rhoevol])-([eq : fevol ] ) , namely , we choose to ignore the fluctuations of the particle density @xmath50 and concentrate on the fluctuations of the conserved charge density @xmath49 . this assumption is somewhat similar in spirit to the approach of glotzer and coniglio  @xcite for the problem of spinodal decomposition , and to the spherical approximation for the ising model in the limit of @xmath57 . unlike the `` classical '' mean - field approach , however , the proposed approximation does take into account exactly the charge density fluctuations , and is expected , therefore , to describe at least some of the features of the fluctuation - dominated kinetics . the justification of the proposed assumption lies in a simple observation that , while average particle density at any time is nonzero , so that @xmath58 is finite and likely to be small , average charge density is always zero , and therefore in comparison the charge fluctuations @xmath59 are large . thus , we expect the former fluctuations to be less important that the latter , and that we can approximate the particle number density by its average ( time - dependent ) value in the equations of motion without losing their important features . in a sense , our approximation is a generalization of an argument used by toussaint and wilczek  @xcite , in which they based their scaling decay law on a suggestion that @xmath60 . it seems clear that the approximation of ignoring the number density fluctuations must breakdown at least below some upper - critical dimension @xmath7 , since asymptotically @xmath50 vanishes . in this case our approximation will be valid for @xmath61 for all times and in systems below @xmath7 it will be a good approximation up to a crossover time beyond which the asymptotics will be modified by the number density fluctuations . systematically taking into account these additional fluctuations will be a subject of future work.@xcite    taking into account the above approximation , we rewrite the equation  ( [ eq : fevol ] ) in fourier representation , taking @xmath62 as a spatially - independent but time - dependent function : @xmath63 where @xmath64 .    the equation  ( [ eq : rhoevol ] ) in this self - consistent approximation is rewritten as : @xmath65 where @xmath66 denotes averaging over initial conditions . these equations of motion have to be supplemented with initial conditions . it is well known that the initial density distribution plays an important role in determining the scaling decay law . although the self - consistent approximation employed here is well suited for a comprehensive study of the influence of initial conditions on the dynamics , here we limit our study to a single type of an initial condition . throughout this paper we will focus on the dynamics initiated with a random gaussian particle distribution , completely characterized by @xmath67 which constraints the system to charge neutrality at all times and for simplicity we take the variance @xmath68 ( @xmath69 in physical units ) . for most charged experimental systems a more relevant initial conditions incorporate the suppression of long - wavelength fluctuations in charge density , which can be modeled ( within e.g. the debye - hckel approximation ) by @xmath70 in equation  ( [ eq : fkinit ] ) . the diffusion - only ( do ) case ( @xmath71 ) and coulombic case ( @xmath72 ) are the simplest systems with a relatively clear , yet interesting asymptotic behavior . all the intermediate interactions ( arbitrary @xmath36 and @xmath73 ) lead to a more complicated scaling behavior , with several regimes and crossovers . we will devote a subsection to each of these three cases . there are two ways of approaching the limit of `` no long - range interactions '' : by decreasing the force constant @xmath54 to 0 or by increasing the power exponent @xmath36 to infinity ( interaction with an effectively vanishing range ) . obviously , these limits should give the same answer . for simplicity , we will set @xmath74 and show that our self - consistent approximation yields the well - known results.@xcite-@xcite @xmath75 for @xmath71 , the kinetic equation for @xmath49 reduces to a simple diffusion equation , with the solution @xmath76 substituting this solution  ( [ eq : twfk ] ) for @xmath77 into equation   ( [ eq : rhonewevol ] ) , and taking into account the initial condition  ( [ eq : f1kinit ] ) - ( [ eq : fkinit ] ) , we obtain : @xmath78 an the exact solution expressible in terms of confluent hypergeometric functions is possible.@xcite it can also be easily shown that equation  ( [ eq : twsoln ] ) describes the asymptotic solution of equation  ( [ eq : twequation ] ) . this is expected , since , as we argued above , the approximations made by toussaint and wilczek  @xcite are very similar to our self - consistent model . the case of @xmath79 will be analyzed in subsection  [ intermediate systems ] , where it will be shown that for all @xmath80 , the decay law is asymptotically the same as for do systems .      in coulombic systems , the long - range interaction is the strongest possible that one can achieve without making the system thermodynamically unstable ( systems with interactions stronger than coulombic have infinite pressure and chemical potential even if their total charge is zero ) . because of this , one would expect the particle density decay for coulombic systems to be very close or equal to the mean - field law @xmath81 . @xcite as we describe below , the self - consistent approximation predicts the decay exponent @xmath26 consistent with this expectation and with some simulations reported in the literature.@xcite however , it can be shown@xcite that correction to this mean - field like decay law can arise from the number density fluctuations and noise for @xmath82 , both of which have been neglected in the self - consistent theory presented here . equation for the evolution of charge density @xmath49 for the coulombic systems in the self - consistent approximation can be exactly solved to yield : @xmath83\;.      \label{eq : coulfksoln}\ ] ] using this solution  ( [ eq : coulfksoln ] ) and the initial condition  ( [ eq : f1kinit ] ) - ( [ eq : fkinit ] ) , @xmath62 can be easily shown to satisfy the following differential equation : @xmath84 in order to find the asymptotic solutions , we introduce a new variable : @xmath85 equation  ( [ eq : coulrhoevol ] ) then transforms to : @xmath86 let us find the `` critical '' dimension @xmath7 , above which the mean - field behavior is manifested . the mean - field solution for @xmath87 is given by @xmath88 where @xmath89 corresponds to subdominant constant terms and terms decreasing with time . by counting powers of @xmath4 in the right - hand side and the left - hand side of equation  ( [ eq : thetaevol ] ) , we obtain :    power(lhs ) = @xmath90 ;    power(rhs ) = @xmath91 . + obviously , for the mean - field solution to be valid asymptotically , the power(lhs ) should be larger than the power(rhs ) , which happens for systems with dimensionality larger than the critical dimension , @xmath92 if @xmath93 , the asymptotic kinetics is determined by the slower process , which is the annihilation , with possibly interaction renormalized @xmath94 ( implicitely assumed here ) , and the coulomb interaction and diffusion are asymptotically irrelevant . for @xmath95 , the coulomb interaction has an interesting effect of continuously lowering the upper - critical dimension from @xmath96 ( for @xmath97 , in which case for @xmath15 the diffusion dominates giving @xmath98 ) down to @xmath7 given above .    in order to analyze the kinetics when space dimensionality @xmath73 is below @xmath7 , we employ the `` steady - state '' approximation , which suggests that at long times the time derivative in the lhs of equation  ( [ eq : coulrhoevol ] ) is the smallest of the three terms . in this case , the equation of motion can be written as : @xmath99 an exact solution of this equation is : @xmath100\ ; , \label{eq : theta2soln}\ ] ] and @xmath101 it can be easily shown that for large @xmath4 , the asymptotic solution for the particle density decay ( @xmath102 ) is : @xmath103 predicting the asymptotic decay exponent @xmath26 for coulomb systems , as in the mean - field regime , although with @xmath54- rather than @xmath104-determined amplitude . this large-@xmath4 limit is achieved when : @xmath105^{{4}/(4-d)}\ ; , \label{eq : coullarget}\ ] ] and it can be easily seen that in the limit of @xmath106 ( vanishing interactions ) the transition time @xmath107 to this region becomes infinite , i.e. this time is never reached . if diffusion is faster than the deterministic coulomb interaction - driven relaxation , i.e. , @xmath108 , then for times less than @xmath109 , the annihilation is governed by the intermediate asymptotics : @xmath110 so the particle density is described by the toussaint - wilczek solution up to the crossover time @xmath107 : @xmath111 this intermediate asymptotics , which exists only for @xmath112 , reflects early times diffusion - dominated decay , with the slower deterministic coulomb interaction driven classical @xmath113 decay appearing only at times later than @xmath107 . in contrast , for @xmath114 or if @xmath115 , there is no extended intermediate regime and one should see a quick transition to a classical decay law . although within the self - consistent approximation the asymptotic @xmath116 decay is not affected by the choice of reasonable initial conditions , the intermediate diffusion - dominated decay is certainly affected by our choice of random gaussian uncorrelated initial conditions given in equation  ( [ eq : fkinit ] ) . for instance if the screened debye - hckel initial conditions are used with @xmath117 , then this @xmath118 will appear as a multiplicative kernel under the @xmath119 integral in equation  ( [ eq : coulrhoevol ] ) . it will then modify the intermediate decay exponent from @xmath120 to @xmath121 ( and , for @xmath53 , eliminating this intermediate region altogether ) , without modifying the asymptotic decay of equation  ( [ eq : coulrholarget ] ) .    to sum up , we find that within the self - consistent approximation , coulombic systems @xmath122 asymptotically exhibit the @xmath113 density decay , consistent with several scaling arguments and simulations.@xcite      let us now consider the general case of long - range interactions with a power - law @xmath123 that is of shorter range ( weaker ) than the coulomb interaction considered in the previous section . equations  ( [ eq : fkevol ] ) - ( [ eq : rhonewevol ] ) can be solved to yield : @xmath124\;,\label{eq : intfksoln}\\      \frac{d \rho}{d t } + { \cal k}\rho^{2 } & = &       { \cal k}\int \frac{d^{d } k}{(2 \pi)^{d } }       \exp [ -2dk^{2}t - 2qk^{2-\sigma } \int_{0}^{t } \rho(\tau ) d \tau]\ ; , \label{eq : intrhoevol}\end{aligned}\ ] ] where @xmath125 . equation  ( [ eq : intrhoevol ] ) is significantly more complicated than its analogs for either coulombic or non - interacting cases . nevertheless , it is possible to find its power - law asymptotic solutions . using an asymptotic analysis analogous to that described in the previous subsection , we find several kinetic regimes depending on the values of @xmath73 and @xmath36 . these regimes depend crucially on the charge density relaxation mechanism , i.e. , whether the lri or diffusion determines the relaxation rate of @xmath126 at late times . in order to analyze the asymptotic behavior of the system , we again introduce the integrated density @xmath127 as defined in equation ( [ eq : theta ] ) . we also assume a power - law for the density and , for @xmath112 , neglect the term @xmath128 . in this case , the equations of motion are : @xmath129}\;,\label{eq : ims1}\\      \rho ( t ) & = & \frac{d \theta}{d t}\;.\end{aligned}\ ] ] depending on @xmath130 and @xmath73 , either the first or second term in the exponential in equation ( [ eq : ims1 ] ) dominates for large @xmath4 , corresponding to either diffusive or superdiffusive relaxation . we first assume that diffusive relaxation is prevalent and determine the conditions when it is true . in the case of diffusive relaxation mechanism and at large @xmath4 , equation  ( [ eq : ims1 ] ) can be simplified to yield : @xmath131}\ ; , \label{eq : ims2}\ ] ] where @xmath132 is a dimensionless constant absorbing integration over angular variables , and the integral becomes time - independent for large @xmath4 if the second term in the exponential vanishes with time . if we assume the power - law dependence for @xmath50 and @xmath87 : @xmath133 then it follows from equation  ( [ eq : ims2 ] ) that for the predominantly diffusive system @xmath134 , as expected . thus , in order for this solution to be self - consistent , we must require that : @xmath135 and , from the definition of @xmath130 , we determine the region where the relaxation and density decay are diffusion - limited : @xmath136 this region is marked fd ( fluctuation - dominated ) in figure 1 . in it , the lri are irrelevant at large @xmath4 , although they may influence the density decay kinetics for intermediate @xmath4 . the asymptotic decay law in the fd region is given by : @xmath137 in the region referred to as the ir ( intermediate region ) and which lies below fd in figure 1 ( @xmath138 ) the lri are strong and dominate the diffusion at large @xmath4 . to investigate the asymptotics of the decay in this region , we rewrite equation  ( [ eq : ims1 ] ) in the following way : @xmath139}. \label{eq : ims4}\ ] ] using @xmath140 . by solving equation  ( [ eq : ims4 ] ) approximately we find two asymptotics in this region : @xmath141 valid at intermediate times , and @xmath142 where @xmath143 for asymptotically large times . the crossover time @xmath144 from the diffusion - dominated decay to the lri - dominated decay is : @xmath145 thus , in this region ( marked ir - intermediate region - in figure 1 ) the lri accelerates the relaxation of the initial density fluctuations and thereby speed up the annihilation . if @xmath146 , the diffusive relaxation and the ( @xmath120)-law may be observed for the intermediate @xmath4 before transition to the superdiffusive relaxation and faster decay takes place for @xmath147 . if @xmath115 ( region mf - mean - field - in figure 1 ) , spatial fluctuations become irrelevant and the classical kinetic - rate equation becomes asymptotically correct , so the decay law in this region is given by : @xmath148 as previously discussed . in the presented work we derive approximate kinetic equations for the annihilation - diffusion process with long - range forces . to analyze the asymptotic decay law for systems with @xmath53 , we proposed a new self - consistent method of calculating the average particle density as function of time . since the total particle density is a nonconserved order parameter with positive average at all times , we argued that its fluctuations are less important in determining dynamics of annihilation than that of a conserved order parameter - the charge density . this approximation self - consistently decouples two kinetic equations and make it possible to find the asymptotic solutions .      for coulombic systems in more than two dimensions , our model yields the mean - field exponent @xmath26 , yet the role of segregation ( i.e. , charge density fluctuations ) is important and can not be simply left out . ispolatov and krapivsky  @xcite proposed the `` unpenetrable domain '' scaling concept , which in coulombic case results in decay exponent @xmath150 independently of space dimensionality . both their model and our self - consistent approximation neglect possible fluctuation modes due to the spatial variation of particle density , as well as noise , which can lead to some slowing down of the reaction kinetics , as indicated by simulations . elucidation of such modes and their role should require detailed account of noise and possibly use of renormalization group analysis when @xmath151 , since it is the critical dimension for the annihilation - diffusion problem . since neglecting of noise in this problem appears to be justified for long - time asymptotics for systems with @xmath53,@xcite it is possible that the exponent @xmath26 for @xmath152 coulombic systems is exact , even though there is no experimental evidence to support this conclusion . the analysis of the self - consistent equations of motion  ( [ eq : fkevol ] ) - ( [ eq : rhonewevol ] ) suggests that there is a new region in the ( n , d ) phase diagram ( the region labeled ir in figure 1 ) , in which the large - t asymptotics of the density decay is determined by long - range forces . the annihilation initally depletes positively charged region of negative particles and vice versa , and then the decay rate is determined by the speed of particle drift from such regions . again , our large - t asymptotics here agrees with the `` unpenetrable domain '' theory of ispolatov and krapivsky , although their model predicts different boundaries of the ir region in the ( n , d ) phase space ( @xmath153 is the upper boundary of the ir region in the self - consistent model , and @xmath154 is the upper boundary of the ir region in the `` unpenetrable domain model '' ; the lower boundary in both theories is the coulombic line @xmath24 ) . the self - consistent theory also predicts a crossover from diffusion - dominated decay to the lri - dominated decay at large times for the systems in this region .    because the self - consistent model is a semi - mean - field approximation ( it completely neglects particle density fluctuations and only takes into account the concerved charge density fluctuations ) , it should be considered only as a first step . a systematic perturbative analysis of equations  ( [ eq : rhoevol ] ) - ( [ eq : fevol ] ) around our self - consistent solution is needed to assess the role of left out number density fluctuations and noise.@xcite    the proposed self - consistent model , its somewhat uncontrolled approximations notwithstanding , represents an important tool in qualitative analysis of dynamic processes in two - component systems with one conserved and one non - conserved variable . it predicts new annihilation behavior ( ir regime ) and new crossovers between diffusion - driven and lri - driven decay regions , reproduces all known results for the annihilation problem in special limits , and can be used to systematically study the role of initial conditions in such processes .                                            by coulombic , we mean such a lri that its potential obeys poisson s equation in d - dimensional system . this transforms to a relationship between system dimensionality d and lr force power exponent n : @xmath24 . k. lindenberg , p. argyrakis and r. kopelman , the hierarchies of nonclassical regimes for diffusion - limited binary reactions . in : fluctuations and order : the new sysnthesis , springer - verlag , 1996 , and references therein . phase diagram of the annihilation - diffusion reaction with long - range forces . fr - forbidden region ( below coulombic line ) . ir - intermediate region - large-@xmath4 asymptotics is determined by lri . fd - fluctuation - dominated region - large-@xmath4 asymptotics is determined by diffusion and initial fluctuations . mf - mean - field region - large-@xmath4 asymptotics is determined by the kinetic - rate equation .
Abstract:

Answer: we introduce coarse - grained hydrodynamic equations of motion for diffusion - annihilation system with a power - law long - range interaction . by taking into account fluctuations of the conserved order parameter - charge density - we derive an analytically solvable approximation for the nonconserved order parameter - total particle density . asymptotic solutions are obtained for the case of random gaussian initial conditions and for system dimensionality @xmath0 . large - t , intermediate - t and small - t asymptotics were calculated and compared with existing scaling theories , exact results and simulation data .


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: the fermion determinant is a highly non - local object as can e.g.  be seen from the hopping expansion for lattice regularized fermions ( see @xcite for a basic introduction ) . the hopping expansion expresses the fermion determinant as the exponential of a sum over all possible closed loops and the external field variables along a loop are collected as factors for this loop . the exponential function can then be expanded and the result is a representation of the fermion determinant in terms of loops . the loop - sum in the exponent , however , contains loops of arbitrary length , which can also iterate parts or all of their contour arbitrarily often . thus , the exponent contains arbitrarily high powers of the external fields . on the other hand we know that , at least on a finite lattice , the fermion determinant is a finite polynomial in the external fields , and thus infinitely many contributions present in the exponent have to cancel each other when expanding the exponential . the result has to be a relatively simple loop representation for the fermion determinant .    for the case of staggered fermions several papers can be found in the literature @xcite-@xcite where polymer representations for the fermions are obtained . for the case of wilson fermions relatively few is known due to the more involved spinor structure of the fermions . here we concentrate on 2-d wilson fermions . an instance where the above discussed cancellation of contributions was brought under control for a model with wilson fermions is salmhofer s mapping of the strongly coupled lattice schwinger model to a self avoiding loop model @xcite . in a first step the gauge fields at strong ( = infinite ) coupling were integrated out and the remaining grassmann integral was then represented as a sum over loops . subsequently scharnhorst studied the two - flavor lattice schwinger model with this method @xcite and extended the techniques to find a two - color loop model for the 2-d lattice thirring model @xcite . in all of these cases , however the external field was either integrated out in the strong coupling limit @xcite or was nt present at all @xcite . on the other hand , our arguments given in the first paragraph show that the cancellation of higher winding loops is a universal phenomenon and that it should be possible to find a simple loop representation also for the fermion determinant in an external field .    in this article we present a simple formula for the fermion determinant of 2-d wilson lattice fermions in a scalar background field ( eq .  ( [ equivalence ] ) ) . the proof is based on a result for the hopping expansion for a generalized 8-vertex model where the vertices are coupled to a background field @xcite . with a proper choice of the vertex weights the square of the partition function for this model represents the 2-d fermion determinant in a scalar background field . the resulting expression reduces the hopping expansion to a finite sum ( on a finite lattice ) of loops . in addition it is possible to explicitly integrate out the external field . by doing so with a gaussian action for the scalar we generate a loop representation of the 2-d gross - neveu model . such loop representations for lattice field theories allow for a considerably more accurate numerical treatment as has e.g.  been demonstrated for the strongly coupled schwinger model @xcite . the basic idea for the proof of our loop representation is to identify the hopping expansion of the wilson fermion determinant with the hopping expansion of a generalized 8-vertex model @xcite . here we briefly rederive the hopping expansion for wilson fermions in a form suitable for comparison with the generalized 8-vertex model ( see e.g.  @xcite for a more detailed discussion of the hopping expansion ) . we study a lattice model of 2-d fermions which in the continuum corresponds to the action @xmath1 \ ; \psi(x ) \ ; , \label{continuum}\ ] ] where @xmath2 is a scalar external field . the lattice fermion determinant is expressed as a path integral @xmath3 \ ; = \ ; \int d\psi d\overline{\psi } \ ;   \exp\left(-\sum_{x , y \in \lambda } \overline{\psi}(x )   m[\theta](x , y ) \psi(y ) \right ) \ ; , \label{partfunct}\ ] ] where the kernel for the lattice fermion action ( wilson fermions ) which regularizes the continuum action ( [ continuum ] ) is given by @xmath4(x , y ) \ ; = \ ; \big[2 + m + \theta(x)\big ] \ ; \delta_{x , y }   \ ; - \ ; \sum_{\mu = \pm 1}^{\pm 2 } \gamma_\mu \delta_{x+\mu , y } \ ; , \ ] ] and we defined @xmath5 \ ; \ ; \ ; \ ; , \ ; \ ; \ ; \ ; \mu = 1,2 \ ; .\ ] ] here @xmath6 are pauli matrices . the sum in the exponent of ( [ partfunct ] ) runs over the whole lattice @xmath7 , which for simplicity we assume to be a finite rectangular piece of zz@xmath8 ( the generalization to e.g.  a torus is straightforward ) . the boundary conditions are open , i.e.  hopping terms that would lead to the outside of our lattice are omitted . we define @xmath9 and assume that @xmath2 is such that @xmath10 for all lattice points @xmath11 . this is a purely technical assumption due to the particular techniques we use for computing the determinant . the final result will be a finite polynomial in the @xmath2 and the above restriction is irrelevant then . we now can write @xmath3 \ ; = \ ; \prod_{x \in \lambda } h(x)^2   \det\big(1 - r[\theta ] \big ) \ ; = \ ; \prod_{x \in \lambda } h(x)^2 \exp \left ( - \sum_{n=1}^\infty \frac{1}{n } \mbox{tr } \ ;   r[\theta]^n \right ) \ ; . \label{hopexp}\ ] ] in the last step the hopping expansion was performed , i.e. the determinant was expressed using the well known trace - logarithm formula and the logarithm was expanded in a power series . the hopping matrix @xmath12 $ ] is defined as @xmath13(x , y ) \ ; = \ ; \sum_{\mu = \pm 1}^{\pm 2 } \ ; \gamma_{\mu } \ ; \frac{1}{h(x ) } \ ; \delta_{x+\mu , y } \ ; . \label{hoppingmatrix}\ ] ] the series in the exponent of ( [ hopexp ] ) converges for @xmath14|| < 1 $ ] , which can be enforced by choosing large enough @xmath15 , i.e.  suitable @xmath2 . again this is only a technical restriction and can be abandoned in the final result . due to the kronecker delta in ( [ hoppingmatrix ] ) , the contributions to tr@xmath12^n$ ] are supported on closed loops on the lattice , and since closed loops are of even length , the contributions for odd @xmath16 vanish . for we obtain @xmath18^{2k } \ ; = \ ;   \sum_{x \in \lambda } \ ; \sum_{l \in { \cal l}^{(2k)}_x } \ ; \prod_{y \in p(l ) } \frac{1}{h(y ) } \ ; \ ; \mbox{tr } \ ; \prod_{\mu \in l } \gamma_\mu \ ; . \label{hoptrace}\ ] ] here @xmath19 is the set of all closed , connected loops of length @xmath20 and base point @xmath11 . by @xmath21 we denote the set of all sites visited by the loop @xmath22 . note that a factor @xmath23 is produced whenever @xmath22 runs through @xmath11 which can be arbitrary often for long enough loops . the last term in ( [ hoptrace ] ) is the trace of the ordered product of the hopping generators @xmath24 as they appear along the loop @xmath22 . we remark , that @xmath25 , which implies that whenever a loop turns around at a site and runs back along its last link this contribution vanishes . thus all these _ back - tracking _ loops can be excluded from @xmath26 .    evaluating the trace over the matrices @xmath24 for a given loop is the remaining problem in the hopping expansion . it first has been solved in @xcite by realizing that the pauli matrices give rise to a representation of discrete rotations on the lattice . alternatively one can decompose the loop using four basic steps and compute the trace in an inductive procedure along the lines of @xcite . the result is @xmath27 by @xmath28 we denote the number of self - intersections of the loop @xmath22 and @xmath29 gives its number of corners . the result is independent of the orientation of the loop . inserting ( [ hoptrace ] ) and ( [ gtrace ] ) in ( [ hopexp ] ) we obtain @xmath3 \ ; =   \ ; \prod_{x \in \lambda } h(x)^2 \exp \left (   \sum_{k=1}^\infty \frac{1}{2k }   \sum_{y \in \lambda } \ ; \sum_{l \in { \cal l}^{(2k)}_y }   ( -1)^{s(l ) } \big(\frac{1}{\sqrt{2}}\big)^{c(l ) }    \prod_{z \in p(l ) } \frac{1}{h(z ) } \right ) . \label{prelim}\ ] ] finally we further simplify this expression by removing the explicit summation over the base points @xmath11 . a loop of length @xmath20 without complete iteration of its contour allows for @xmath20 different choices of a base point thus cancelling the factor @xmath30 in ( [ prelim ] ) . a loop which iterates its whole contour @xmath31 times allows only for @xmath32 different base points and a factor @xmath33 remains . the final expression is @xmath3 \ ; =   \ ; \prod_{x \in \lambda } h(x)^2 \ ; \exp \left ( 2 \sum_{l \in { \cal l } }   \frac{(-1)^{s(l)}}{i(l ) } \big(\frac{1}{\sqrt{2}}\big)^{c(l ) } \ ;   \prod_{y \in p(l ) } \frac{1}{h(y ) } \right ) \ ; , \label{finalhop}\ ] ] where @xmath34 is the set of all closed , connected , non back - tracking loops of arbitrary length . each loop is included in @xmath34 with only one of its two possible orientations and we collect an overall factor of 2 in the exponent . as already outlined , the next step is to compare the result ( [ finalhop ] ) to the hopping expansion for a generalized 8-vertex model studied in detail in @xcite . in this generalized model , the vertices are coupled to a locally varying external field @xmath35 . before we proceed let s first discuss this generalized model and its relation to the standard 8-vertex model . the standard 8-vertex model @xcite can be viewed as a model of 8 quadratic tiles ( vertices ) and each of them is assigned a weight @xmath36 ( compare fig .  [ tiles ] ) . a _ tiling _ of our lattice @xmath7 ( the same rectangular piece of @xmath37 as before ) is a covering of @xmath7 with the tiles such that on each site of @xmath7 we place one of our tiles with the centers of the tiles sitting on the sites . the set @xmath38 of _ admissible tilings _ is given by those arrangements of tiles where the black lines on the tiles never have an open end . the partition function of the standard 8-vertex model is the sum over all admissible tilings @xmath39 and the boltzmann weight for a particular tiling @xmath40 is given by the product of the weights @xmath41 for all tiles used in this tiling @xmath40 .    in our generalization of the model we now couple the vertices to an external field @xmath35 located on the sites @xmath11 of @xmath7 . the partition function of the generalized model is given by @xmath42 \ ; = \ ;   \sum_{t \in { \cal t } } \ ; \prod_{i = 1}^8 { w_i}^{n_i(t ) } \ ! \prod_{x \in p(t ) } \varphi(x ) \ ; . \label{8vertex}\ ] ] here @xmath43 denotes the set of all sites occupied by the tiling @xmath40 . when a site @xmath11 is occupied by tile nr .  2 , this site is counted twice giving a factor @xmath44 . in case @xmath11 is occupied by tile nr .  1 , @xmath45 and the factor is 1 . for all other tiles @xmath11 is counted once and the factor is @xmath35 . the generalized model ( [ 8vertex ] ) also allows for a hopping expansion which is derived in @xcite . the central step is to rewrite the partition function as an integral over grassmann variables along the lines of @xcite . the action is a quadratic form in the grassmann variables and the partition function gives rise to a pfaffian . the pfaffian however can be expanded similar to the expansion in section 2 and it is furthermore possible to explicitly evaluate the traces over the hopping matrix . the resulting expression reads @xcite @xmath42 \ ; = \ ;   ( -1)^{|\lambda| } w_1^{|\lambda| } \ ;   \exp \left ( \sum_{l \in { \cal l } } \frac{(-1)^{s(l)}}{i(l ) }   \big ( \frac{1}{w_1 } \big)^{|l| } \prod_{i=3}^8 w_i^{n_i(l ) } \ ; \prod_{x \in p(l ) } \varphi(x )    \right ) . \label{8vertexhop}\ ] ] as for the wilson fermions the sum runs over the set @xmath34 of all closed , non back - tracking loops of arbitrary length . each loop @xmath22 is included in @xmath34 with only one of its two possible orientations . by @xmath46 we denote the length of the loop , @xmath47 is the number of iterations of its complete contour and @xmath28 is the number of self - intersections . @xmath21 again denotes the set of sites visited by the loop @xmath22 with @xmath11 being included in @xmath21 whenever the loop @xmath22 runs through @xmath11 . by @xmath48 we denote the size of the lattice . in the following we will work on lattices with an even number of sites and the overall sign factor is 1 .    the exponents @xmath49 in ( [ 8vertexhop ] ) give the numbers for the abundance of the line elements as they are depicted in fig . [ tiles ] . e.g.  when the loop @xmath22 changes from heading east to heading north at a site it picks up a factor of @xmath50 and similarly for the other tiles @xmath51 . note that the loops @xmath52 occur as an ordered set of instructions for the directions the loop takes as it hops from one site to the next . the element corresponding to the tile 2 with weight @xmath53 ( compare fig .  [ tiles ] ) is not needed to describe the loop @xmath22 . equation ( [ 8vertexhop ] ) does not contain @xmath53 explicitly at all . the weight @xmath53 is related to the other weights through the free fermion condition @xcite @xmath54 the free fermion condition is a sufficient condition for finding an explicit solution to the standard 8-vertex model without external fields . the above mentioned grassmann representation automatically enforces the free fermion condition as is discussed in more detail in @xmath55 .    by comparing ( [ finalhop ] ) with ( [ 8vertexhop ] ) we find that the sums over the loops in the exponent become identical when setting @xmath56 using the free fermion condition ( [ freecond ] ) we find @xmath57 . this implies that when representing wilson fermions using eq . ( [ 8vertex ] ) , the sum over the tilings can be replaced by a sum over the set @xmath58 of closed , self - avoiding loops , i.e.  loops that are not allowed to self intersect or touch each other . the loops can have several disconnected but closed pieces and each piece is included in @xmath58 with only one of its two possible orientations . it is important to note , that the exponents in ( [ finalhop ] ) and ( [ 8vertexhop ] ) differ by an overall factor of 2 . this is due to the fact that the action for the wilson fermions is a bilinear form giving rise to a determinant while the grassmann action for ( [ 8vertex ] ) is a quadratic form giving rise to a pfaffian when integrating out the grassmann variables . since the kernel of the latter action is anti - symmetric the pfaffian , however , is given by the square root of a determinant causing the difference by a factor of 2 in the exponent .    putting things together we obtain the final formula @xmath3 \ ; = \ ; \prod_{x \in \lambda } [ 2 + m + \theta(x)]^2 \ ; \left ( \sum_{l \in { \cal l}_{sa } } \big ( \frac{1}{\sqrt{2 } } \big)^{c(l ) } \prod_{y \in p(l ) } \frac{1}{2 + m + \theta(y ) } \right)^2 \ ; , \label{equivalence}\ ] ] where as above @xmath29 is the number of corners of @xmath22 . the loops in @xmath58 are self - avoiding and thus for a given loop configuration @xmath22 each site of the lattice is occupied only once . thus in the sum each inverse field @xmath59^{-1}$ ] can only occur linearly . when taking the square of this sum we can only produce terms which are at most quadratic in the inverse field . the overall factor , however , still can cancel the quadratic terms and the final result is a finite polynomial in the fields @xmath2 . the above imposed technical restrictions on the range of @xmath2 can now be lifted . + several remarks on the result ( [ equivalence ] ) are in order : when setting all external fields to zero we find that free wilson fermions are equivalent to the square of the partition function of the self avoiding loop model @xcite with bending rigidity @xmath60 and bond weight @xmath61^{-1}$ ] . thus for a trivial background field we reproduce the result for free fermions obtained by scharnhorst @xcite with a different method . ( [ equivalence ] ) is a direct generalization of the trivial case to fermions in a scalar background field .    from an algebraic point of view eq . ( [ equivalence ] ) is exactly the expression we had in mind when discussing the general algebraic structure of the determinant . in two dimensions the fermion determinant has to be a polynomial which can at most be quadratic in the external field . the feature that the terms in this polynomial are organized according to closed loops is inherited from the hopping expansion . actually , for single contributions to ( [ equivalence ] ) it is even possible to trace their emergence from an expansion of the exponential in ( [ finalhop ] ) . when doing so , one finds that the intersection factor @xmath62 provides the mechanism which ensures the cancellation of loops with multiply occupied links . ( [ 8vertexhop ] ) establishes that this cancellation mechanism is also independent of the corner weights . finally , it is obvious from ( [ equivalence ] ) that in the loop representation it is straightforward to integrate out the scalar fields . this can be done with different actions for @xmath63 . in the next section we discuss the case of a simple gaussian which will produce the gross - neveu model @xcite . we remark , that a loop representation for the gross - neveu model with staggered fermions has been analyzed in @xcite and a numerical study of the model with conventional methods ( introduction of the auxialiary field ) is given in @xcite . in order to give an application of our eq . ( [ equivalence ] ) we now integrate the scalar field with a gaussian measure @xmath64 \ ; = \ ;   \int \prod_{x \in \lambda}\frac{d \theta(x)}{\sqrt{2 \pi g } } \exp\left(- \frac{1}{2 g } \theta(x)^2 \right ) \ ; .\ ] ] when integrating @xmath65)^n$ ] with this measure we generate a lattice version of the @xmath0-component gross - neveu model @xcite with ( continuum ) action @xmath66 \psi(x ) \ ; - \ ; g^2 \big [ \overline{\psi}(x ) \psi(x ) \big]^2 \big ) \ ; , \label{grossneveu}\ ] ] where @xmath67 now have @xmath0 flavor components . using ( [ equivalence ] ) the lattice partition function for the gross - neveu model reads @xmath68 \prod_{x \in \lambda }   [ 2\!+\!m\!+\!\theta(x)]^{2n } \ ; \left ( \sum_{l \in { \cal l}_{sa}}\!\big ( \frac{1}{\sqrt{2 } } \big)^{c(l ) } \!\prod_{y \in p(l ) } \frac{1}{2\!+\!m\!+\!\theta(y ) } \right)^{2n }   \nonumber \\ & = & \int d\mu[\theta ] \sum_{l_1 , ... , l_{2n } } \big ( \frac{1}{\sqrt{2 } } \big)^{c(l_1 ) + ... + c(l_{2n } ) } \prod_{x \in \lambda } [ 2\!+\!m\!+\!\theta(x)]^{2n -   o_x(l_1 , ... , l_{2n } ) } \nonumber \\   & = & \sum_{l_1 , ... , l_{2n } } \big ( \frac{1}{\sqrt{2 } } \big)^{c(l_1 ) + ... + c(l_{2n } ) } \prod_{x \in \lambda } g[2n - o_x(l_1 , ... , l_n ) ] \ ; . \label{loopgn}\end{aligned}\ ] ] the function @xmath69 counts how many of the independent loops @xmath70 occupy the site @xmath11 for a given loop configuration . each loop can either leave the site empty or occupy it once ( the loops are self - avoiding ) and thus @xmath71 has values between 0 and @xmath72 . the function @xmath73 , j = 0 , 1 , ... , 2n$ ] is given by @xmath74 \ ; = \ ; \int\!\frac{d\theta}{\sqrt{2\pi g } } e^{- \frac{1}{2 g } \theta^2 } [ 2\!+\!m\!+\!\theta ] ^j \ ; = \ ;   \sum_{k=0}^{[j/2 ] } { j \choose 2k } ( 2k-1 ) ! ! \ ; g^k \ ; ( 2+m)^{j-2k}.\ ] ] equation ( [ loopgn ] ) establishes that the @xmath0-component gross - neveu model is equivalent to a model of @xmath72 independent self - avoiding loops . the partition function is a sum over all loop configurations with a weight corresponding to the total number of corners times a simple function of the occupation number for each site . the partition function can easily be reformulated as a @xmath75-vertex model . the new tiles are obtained by using @xmath72 different colors and with each color one of the line patterns of fig . [ tiles ] is drawn onto our new tile , giving a total of @xmath75 different tiles . the weight for the tiles is given by a product of the weights @xmath41 for each color ( 1 for @xmath76 , 0 for @xmath53 and @xmath60 for @xmath77 ) times the occupation factor @xmath73 $ ] . + let s discuss the case of @xmath78 in more detail . here the tiles show lines in two colors , say red and blue . the weight factors are the products of the @xmath41 for the blue and red lines multiplied by 1 when the tile is empty , by a factor of @xmath79 when there is only one color and a factor of @xmath80 when both colors are used on the tile . when setting @xmath81 the weights can be factorized into two terms corresponding to the two colors and we recover the result for the free case already discussed above . a second choice of parameters also leads to a particularly simple model . when setting @xmath82 , we find that all tiles which show only one color vanish . thus for a non - vanishing contribution every site has to either be empty or is visited by both , a blue and a red loop . for our simple rectangular lattice this implies , that the red and blue loop configurations have to sit on top of each other and integrating out the scalar field at @xmath83 has produced a dirac delta on the space of loops . the resulting model is again a self avoiding loop model , now with bending rigidity @xmath84 and bond weight @xmath85 . this is the same model which was shown @xcite to be equivalent to the strong coupling schwinger model . it has been analyzed with monte - carlo methods in @xmath86 where the existence of a phase transition at @xmath87 was established . for a study of this model , using computer algebra on small lattices see @xcite . in the form of the above discussed two - color ( 64-vertex ) model , the @xmath88 gross - neveu model is now relatively simple to analyze numerically in the whole @xmath89-plane . in particular , the representation as a vertex model may allow for the use of efficient cluster algorithms @xcite . the motivation for this article was to find simple representations for the wilson lattice fermion determinant in an external field . counting the powers of the external field in the determinant one finds that a large set of contributions appearing in the standard hopping expansion has to cancel . the determinant can only be a finite polynomial in the external variables . finding , however , a reasonably simple and useful expression for the determinant is a hard problem .    in this article we succeeded in finding such a simple representation for the case of 2-d wilson fermions in a scalar background field . the determinant can be written as the product of two self - avoiding loop models coupled to the external field . to prove this result , grassmann techniques and hopping expansion for a generalized 8-vertex model were used . in the obtained loop representation it is straightforward to integrate out the scalar fields , and the application to the case of the gross - neveu model was discussed in more detail .    when one instance of a simple expression for a wilson type fermion determinant in an external field can be found it is natural to ask whether other and more realistic models also allow for a simple representation of the fermion determinant . the next candidate are 2-d fermions interacting with an abelian vector field . the vector field is coupled to the fermions through its gauge transporters which are supported on links instead of sites . this does not pose a fundamental problem since the generalized 8-vertex model @xcite can also be formulated for external field living on links . all terms of the hopping expansion can still be evaluated explicitly . a certain difficulty is , however , given by the fact that the link variables have to be complex conjugated when hopping in backward direction . thus the loops in the hopping expansion give contributions complex conjugated to each other when the orientation is reversed . we believe , however , that when properly organizing the loops this problem can be overcome . a hint in this direction is scharnhorst s proof for the existence of a vertex model for the 2-d thirring model @xcite . if a loop representation for the determinant in an external vector field exists then it should be possible to obtain scharnhorst s result by integrating out the vector field , similar to the loop representation of the gross - neveu model which was obtained by integrating out the scalar field . the case of non - abelian gauge fields poses the additional difficulty , that due to the traces over the gauge field matrices the product of two loops can not be written as a new loop . nevertheless , the power counting argument discussed in the introduction still suggests , that representations simpler than the standard hopping expansion should exist . similar in 4 dimensions where at least the problem of computing the trace of the 4-d @xmath90-matrices is solved @xcite . +   + _ acknowledgement : _ the author thanks christian lang , klaus scharnhorst and uwe - jens wiese for discussions and remarks on the literature .      i. montvay and g. mnster , _ quantum fields on a lattice _ , cambridge university press , cambridge 1994 . p. rossi and u. wolff , nucl . b248 ( 1984 ) 105 . m. karowski , r. schrader and h.j . thun , commun . math . phys . 97 ( 1985 ) 5 . i. montvay , phys . lett . b227 ( 1989 ) 260 . i. montvay , _ probabilistic methods in quantum field theory _ , p.h . damgaard et al . plenum press , new york 1990 . m. salmhofer , nucl . b362 ( 1991 ) 641 . k. scharnhorst , nucl . b479 ( 1996 ) 727 . k. scharnhorst , nucl . b503 ( 1997 ) 479 . c. gattringer , _ a formula for the hopping expansion of 8-vertex models coupled to an external field _ , report cond - mat/9811139 . h. gausterer , c.b . lang and m. salmhofer , nucl . b388 ( 1992 ) 275 ; h. gausterer and c.b . lang , nucl . b455 ( 1995 ) 785 . stamatescu , phys . d25 ( 1981 ) 1130 . gattringer , s. jaimungal and g.w . semenoff , _ loops , surfaces and grassmann representation for two- and three - dimensional ising models _ , report hep - th/9801098 , to appear in int . j. of mod . a. ; + c. gattringer , _ hopping expansion as a tool for handling dual variables in lattice models _ , report hep - lat/9809176 , parallel talk given at ` lattice-98 ' , to appear in the proceedings . c. fan and f.y . wu , phys . 179 ( 1969 ) 179 . c. fan and f.y . wu , phys . rev . b2 ( 1970 ) 723 . r. baxter , _ exactly solved models in statistical mechanics _ , academic press , london 1982 . s. samuel , j. math . 21 ( 1980 ) 2806 . e. fradkin , m. srednicki and l. susskind , phys . d21 ( 1980 ) 2885 . c. itzykson , nucl . b210 ( 1982 ) 448 . priezzhev , sov . jetp 47 ( 1978 ) 619 . f. rys and w. helfrich , j. phys . a15 ( 1982 ) 599 . gross and a. neveu , phys . d10 ( 1974 ) 3235 . y. cohen , s. elitzur and e. rabinovici , nucl . phys . b220 ( 1983 ) 102 . f. karsch , e. meggiolaro and l. turko , phys . d51 ( 1995 ) 6417 . evertz , g.  lana and m.  marcu , phys . 70 ( 1993 ) 875 .
Abstract:

Answer: we show that the fermion determinant for 2-d wilson lattice fermions coupled to an external scalar field is equivalent to self avoiding loops interacting with the external field . in an application of the resulting formula we integrate the scalar field with a gaussian action to generate the @xmath0-component gross - neveu model . the loop representation for this model is discussed . mit - ctp-2800 + hep - lat/9811014    * christof gattringer *    massachusetts institute of technology + center for theoretical physics + 77 massachusetts avenue + cambridge ma 02129 , usa    pacs : 11.15.ha + key words : lattice field theory , fermion determinant


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: collisional @xmath0-body simulations , in which the equations of motion of @xmath0 particles integrated numerically , have been extensively used in studies of dense star clusters , such as globular cluster , open cluster , and clusters with black holes , and also in studies of planetary formation . one feature of the collisional @xmath0-body simulations is need for relatively high accuracy in the force calculations , because the total number of timestep is very large to simulate relatively long simulation span , such as the relaxation timescale . another feature is a wide difference in orbital timescale of particles since two particles can approach arbitrary close . the individual timestep algorithm , first developed by aarseth ( @xcite ) , has been a powerful tool that handles the collisional @xmath0-body system , whose basic idea is to assign different times and timesteps to particles in the system . grape(gravity pipe)(@xcite ) is a special purpose hardware that can accelerate the individual timestep algorithm . grape hardware has specialized pipelines for the gravitational - force calculation , which is the most expensive part of the collisional @xmath0-body simulations . among the individual timestep algorithm , the hermite scheme ( @xcite ) can efficiently use the grape hardware , in which the block individual timestep algorithm ( @xcite ) and the 4th - other hermite integration are used . grape-6 ( @xcite ) is a massive - parallel hardware for the collisional @xmath0-body simulations using the hermite scheme . it consists of 1728 pipeline chips and has a peak speed of around 64 teraflops . although direct summation algorithm was used for the force calculations on the grape-6 system , whether it is the really best solution or not remains a question . the barnes - hut tree algorithm ( @xcite ) is one of algorithms that reduces the calculation cost by replacing forces from distant particles by those from a virtual particle at their center of mass . mcmillan and aarseth ( @xcite ) have demonstrated that it is possible to implement a combination of the barnes - hut tree algorithm and the individual timestep algorithm that runs efficiently on a single - processor computers . however , on the grape-6 system , the combination of the tree algorithm and the individual timestep algorithm was not possible , because its memory addressing scheme was limited only to sequential access to a full set of particle data , and there is not enough memory size for particle data . we successfully implemented the combination of the tree algorithm and the individual timestep algorithm on grape-9 system . grape-9 is a newly - developed system that uses fpga(field programmable gate array ) device and the force and predictor pipelines the same as grape-6 chip are integrated in the device . the grape-9 system also has an indirect memory addressing unit and a relatively large - sized particle memory , implemented by widely - used dram device . interaction lists for the tree algorithm can be stored in the grape-9 system , and the force pipelines can calculate only the interactions necessary . by our implementation , the interaction calculations are significantly reduced from the direct summation in the hermite scheme . the plan of this paper is as follows . in section 2 we describe implementation of the tree algorithm on the hermite scheme using grape-9 . in section 3 , we present the performance and accuracy of our implementation . section 4 is for discussion . in this section , we describe how the interactions are calculated using the tree algorithm in our implementation . in an ideal way , the interaction list should be created at every block timestep using predicted particle data , but it is not practical . in our implementation , the tree structure and the interaction lists are created only at intervals of @xmath3 , and the same interaction lists are used during @xmath3 . therefore , the interval @xmath3 becomes a cause of error in the interactions calculation since the tree structure is deformed as time advances . the interval @xmath3 has to be small enough not to affect simulation results and its performance . in our implementation , the maximum size of timestep is set to be @xmath3 for simplicity .    with the original hermite scheme , the previous grape system ( grape-6 system ) performs the integration of one step in the following way :    * as the initialization procedure , the host computer sends all data of all particles to the memory on grape . * the host computer selects particles to be integrated at the present system time . * repeat 4 - 6 for all particles selected . * the host computer predicts the position and velocity of the particle , and sends them to grape . * grape calculates the force from all other particles , and then returns the results to the host computer . * the host computer integrates the orbits of the particles and determines the new timestep . the updated particle data are sent to the memory on grape . * the host computer updates the present system time and go back to step 2 .    in our new implementation with the tree algorithm , the grape system ( grape-9 system ) performs the integration of one step in the following way ( the bold item number shows the step that changes from the original algorithm ) :    * at intervals of @xmath3 ( and at initial ) , the host computer makes tree data . the procedure includes construction of a tree structure , identification of groups of particles for which the same interaction list is used by traversing the tree structure , and creation of interaction lists for the groups . the host computer sends all particles data , the interaction lists for all groups , and tree node data listed up in the interaction lists , to the memory on grape . the tree node data are stored in the memory as ( pseudo-)particles that have positions , velocities , accelerations , and their time derivatives . * the host computer selects particles to be integrated at the present system time . * repeat 4 - 6 for all particles selected . * the host computer predicts the position and velocity of the particle , and sends them to grape . * the host computer sends the index number of the interaction list for the particle to grape . * grape calculates the force from particles in the interaction list , and then returns the results to the host computer . * the host computer integrates the orbits of the particles and determines the new timestep . the updated particle data are sent to the memory on grape . * the host computer updates the present system time and go back to step 2 . the differences from the original algorithm are in three steps : in step 1 , at intervals of @xmath3 , the tree structure and the interaction list are created and sent to grape . in step 4a , the index number of the interaction list is sent to grape . in step 5 , the force are calculated from particles in the interaction list , instead of from all particles .    in order to use efficiently the grape hardware , we use the modified tree algorithm ( as already described in step 1 ) , which was developed by barnes ( @xcite ) and implemented on the grape hardware by makino ( @xcite ) . with this algorithm , tree traversal is performed for a group of neighboring particles and an interaction list is created for the group . the maximum number of particles in the group , @xmath4 , is set to be optimal at which the total computing time is the minimum . as we increase @xmath5 , the interaction calculation on grape increases since interactions between particles in a group are calculated directly and the interaction list becomes longer . on the other hand , as we decrease @xmath5 , an efficiency of usage of the grape hardware becomes lower , since the number of particles in the same group at each block step becomes smaller on average . for the present system , @xmath6 , is close to optimal . note that , with such @xmath5 , interactions with a rather large number of the neighboring particles , about @xmath7 ( for @xmath8 ) , are directly calculated . the part of the tree algorithm in our code is almost the same as that used in the previous studies ( @xcite,@xcite ) . the simulation program is written using the grape-6 compatible api library and two additional functions for step 1 and steps 4a-5 , respectively . we implemented these algorithm on the grape-9 ( model 5000 ) system . the grape-9 system consists of 8 - 16 grape-9 cards , connected to the host computer via a pci express switch device ( plx pex8696 ) . the grape-9 card is a pci express extension card on which one fpga(field programmable gate ) device and one ddr2 sdram ( so - dimm module ) memory are mounted . in the fpga device , the force and predictor pipelines , almost identical to the grape-6 chip , and an indirect memory addressing unit are integrated , which is illustrated in figure [ fig : g9 ] . the interaction lists for the tree algorithm is stored in the indirect memory addressing unit of the fpga device , actually in on - chip memory of the fpga device , and all particles data and tree node data are stored in the memory unit , which consists of the ddr2 sdram memory . according to the interaction list , the indirect memory addressing unit outputs an address entry for the memory unit , and the force pipelines calculate only the interactions necessary . we use altera cyclone v 5cgxc9 for the fpga device . the wealthy amount of the on - chip memory in this device is one of reasons that enables our implementation . for the present implementation , we use a configuration in which grape-6-compatible 14 force pipelines and one predictor pipeline are integrated , which operates at 98mhz . other details on the grape-9 system will be discussed elsewhere .    as for the indirect memory addressing unit , we use the particle index unit same as in grape-5 ( @xcite ) , which was designed for the cell - index method ( @xcite ) to handle short - range forces in a periodic boundary condition . figure [ fig : indirect ] shows a block diagram for the indirect memory addressing unit . it consists of the cell - index memory and two counters : the cell counter and the particle index counter . in the cell - index memory , sets of start address and count number are stored . according to the output of the cell - index memory , the particle index counter generates entries to the memory unit . the cell counter indicates address entry for the cell - index memory . actually , in step 4a , the host computer sends the start address and count number of the cell counter for the group of the interaction list . because size of the on - chip memory of the fpga device is limited ( about several mbits ) , we store the interaction list for the tree algorithm in such form , instead of full sets of indices . the entry size of the cell - index memory is 98304 for the present fpga device . in order to reduce total length of the interaction lists in the cell - index memory , we rearrange all particles in the peano - hilbert order and store the tree node data for each group in a consecutive location in the memory unit . since we typically use 1 gb ddr2 sdram so - dimm ( 8 gb at maximum ) , the memory unit can store 10 millions particles at the maximum for each card .    when we perform calculations using multiple grape-9 cards , we use two parallelization methods in combination : ( 1 ) multiple cards calculate the force on the same set of ( @xmath9-)particles , but from different set of ( @xmath10-)particles . ( called @xmath10-parallel ) ( 2 ) multiple cards calculate the force on different set of ( @xmath9-)particles whose interaction list ( group ) are different . ( called @xmath9-parallel ) in this section , we discuss performance and accuracy for our implementation . as the benchmark runs , we integrated the plummer model with equal - mass particles . we use the standard unit ( @xcite ) in which @xmath11 . here @xmath12 and @xmath13 are the total mass and energy of the system , and @xmath14 is the gravitational constant . the timestep criterion is that of aarseth ( @xcite ) with @xmath15 . for the softening parameter , we used an @xmath0-dependent softening , @xmath16 . we set @xmath17 for all runs except for runs of @xmath18 , @xmath0=512k and 1 m ( @xmath19 and @xmath20 ) , and the interval @xmath21 . we used the grape-9 system that consists of 8 grape-9 cards and whose peak speed is 630 gflops . here we count operations for the gravitational force and its time derivatives as 57 floating - point operations . host computer has intel core i7 - 3820 ( 4core , 3.6ghz ) cpu . communications between the host computer and each grape-9 card is pci express gen1 4lane ( 1gb / s peak for each direction ) . in order to use 8 cards simultaneously , we used a parallelization method whose degree is 4 for @xmath10-parallel and 2 for @xmath9-parallel . figure [ fig : nt ] shows the calculation time , @xmath22 , to integrate the system for one time unit as a function of the number of particles , @xmath0 , for @xmath23 and @xmath24 , where @xmath25 is a opening parameter for the tree algorithm . for comparison , we also plot the calculation time for the original hermite scheme . we measured the calculation time , @xmath22 , from simulation time @xmath26 to @xmath27 ( and multiplied it by four ) to avoid the complication due to the startup procedure . figure [ fig : nsp ] shows another plot with an equivalent - performance , @xmath28 , defined by @xmath29 where @xmath30 is the total number of individual timestep to integrate one time unit , and @xmath31 is that for the hermite scheme . the equivalent - performance means the performance in the case that we perform the same simulation within the same time using the hermite scheme . for the ratio , @xmath32 , of @xmath331 m , we used instead those of @xmath33512k for each @xmath25 . the ratio @xmath32 itself is close to unity , for example , @xmath34 for @xmath33512k , @xmath35 , which means even if we use the tree algorithm total number of individual timestep does not increase so much . we can see that about a factor ( @xmath36k ) of speedup ( for @xmath8 ) against the original hermite scheme is achieved . figure [ fig : err256k ] shows errors in the total energy as a function of time up to simulation time @xmath37 for @xmath38k . the calculation of the potential energy is obtained with the direct summation on grape , not with the tree algorithm . we can see that the errors in our implementation of the tree algorithm increase linearly as time advances . in figure [ fig : err ] , the errors in the total energy at time @xmath37 are summarized . from both figures , we can see that the errors for @xmath39 are comparable to those for the original hermite scheme . [ tab : timing ]    0    [ tab : tree ]    we discuss breakdown of the calculations using a simple performance model . the calculation time per one particle step is expressed as @xmath40 in table [ tab : timing ] , the terms of the right - hand side measured in actual runs are listed . the first term of the right - hand side , @xmath41 , is the time to calculate the force and its time derivative for one particle on grape-9 , expressed as @xmath42 where @xmath43 is average numbers of the interaction list for the tree algorithm . in the case of one grape-9 card , @xmath44 ( s ) . the factor @xmath45 expresses a decrease in performance when the number of particles that calculate interactions simultaneously ( at step 5 ) is less than @xmath46 . here , @xmath46 and @xmath47 are the maximum and average number of the particles that calculate interactions simultaneously , respectively . the number @xmath48 becomes much smaller for the tree algorithm than in the original hermite scheme , because , even in the same block step , the particles belong to several different groups . the number @xmath49 in actual runs , which does not depends on @xmath0 nor @xmath25 so much . for the present system , @xmath50 , since it has 14 real force pipelines and each real pipeline serves as 4 virtual multiple pipelines ( @xcite ) . the second term , @xmath51 , is the time for the tree data processing spent in the host computer ( step 1 ) , which are listed in table [ tab : timing ] for @xmath52 . the time @xmath51 is proportional to @xmath53 and does not depend on @xmath0 so much . the third term , @xmath54 , expresses the time to transfer data between the host computer and grape , which include data conversion . since about 200 byte data transfer are required per one particle step , the sustained transfer speed is about 200mb / s . the fourth term , @xmath55 , is the time for the host computer to perform computations to integrate one particle other than @xmath51 .    as for the breakdown , at first , we note that , in the first term @xmath56 , the decrease in performance due to small @xmath57 is rather large and the sustained performance decreases to about half of its peak performance . this is partly because @xmath58 for the present system is not small enough . the system which has @xmath46 less than 30 is desirable for our implementation . second , the largest term is @xmath54 among three terms other than @xmath56 , and the fraction to @xmath56 is not small compared in the original hermite scheme because the number of interaction @xmath43 is not large , of course . third , the time @xmath51 is small enough compared to other terms , in the case of @xmath52 . we successfully implemented the hierarchical tree algorithm on the individual timestep algorithm ( the hermite scheme ) for collisional @xmath0-body simulations on the grape-9 system . the present grape-9 system has the indirect memory addressing unit and the memory unit large enough to store all particles data and also tree nodes data . in our implementation , the interaction calculations are significantly reduced , compared to direct @xmath1 summation in the original hermite scheme .    in comparison to other methods that also reduce calculation amounts for the individual timestep algorithm successfully , our implementation has one advantageous feature that interactions from particles at an intermediate range are evaluated in more accurate way . the neighbor scheme ( @xcite , @xcite ) is an example of such methods . in the scheme , the force on a particle is divided into two components , the neighbor force and the regular force and calculations amount are reduced by evaluating the regular force less frequently . p@xmath59 t ( particle - particle particle - tree , @xcite , @xcite ) is another example . in p@xmath59 t , the force on a particle is split into short - range and long - range contribution . the short - range force are evaluated with the hermite scheme and the long - range force are evaluated with the tree algorithm and leapfrog integrator . it is reported that less accurate evaluation for intermediate range force might influence the angular momentum evolution ( see @xcite ) .    at present , the grape-9 system is probably good solution for the implementation of the hierarchical tree algorithm on the individual timestep algorithm . further improvement with the next generation fpga device would provide more powerful computing systems . shipment within one year of a new fpga device ( altera arria 10 ) that has more than 4 times number of logic elements , 3 times operation speed , 8 times data transfer speed ( pcie gen3 8lane ) , 4 times size of the on - chip memory and 10 times memory bandwidth ( ddr3/ddr4 sdram ) , compare to the current fpga device ( altera cyclone v ) , has been announced . new system using such fpga device would be able to provide about 10 times of performance with keeping smaller @xmath60 , which is another required ingredient for our implementation . porting of our implementation on other accelerators , such as gpgpu device , is presumably feasible and in preparation . typically , very large number of parallel operations must be executed on such accelerator . since , for our implementation , the number of the interaction calculations that can be executed in parallel becomes smaller , some ingenuities would be necessary for an efficient use of the accelerator . aarseth ,  s.  j. 1963 , mnras , 126 , 223 aarseth ,  s.  j. 1999 , celest . mech . , 73 , 127 ahmad ,  a. & cohen ,  l. 1973 , j. comput . phys , 12 , 389 barnes ,  j.  e. 1990 , j. comput . phys , 87 , 161 barnes ,  j.  e. , & hut ,  p. 1986 , nature , 326 , 446 fukushige ,  t. , makino ,  j. , & kawai ,  a. 2005 , pasj , 57 , 1009 heggie ,  d.  c. , & mathieu ,  r.  d. 1986 , in the use of supercomputer in stellar dynamics , ed . p.hut & s.mcmillan ( new york : springer ) , 233 iwasawa ,  m. , portegies zwart ,  s. , & makino ,  j. 2015 , computional astrophysics and cosmology , 2 , 6 kawai ,  a. , fukushige ,  t. , makino ,  j. , & taiji ,  m. 2000 , pasj , 52 , 659 makino ,  j. 1991a , pasj , 43 , 621 makino ,  j. & aarseth ,  s.  j. 1992 , pasj , 44 , 141 makino ,  j. , fukushige ,  t. , koga ,  m. , & namura ,  k. 2003 , pasj , 55 , 1163 makino ,  j. , taiji ,  m. , ebisuzaki ,  t. , & sugimoto ,  d. 1997 , apj , 480 , 432 mcmillan ,  s.  l.  w. 1986 , in the use of supercomputer in stellar dynamics , ed . p.hut & s.mcmillan ( new york : springer ) , 156 mcmillan ,  s.  l.  w. & aarseth ,  s.  j. 1993 , apj , 414 , 200 nitadori ,  k. & aarseth ,  s.  j. 2012 , mnras , 424 , 545 oshino ,  s. , makino ,  j. , & funato ,  y. 2011 , pasj , 63 , 881 quentrec ,  b. & brot ,  c. 1975 , j. comput . phys , 13 , 430 sugimoto ,  d. , chikada ,  y. , makino ,  j. , ito ,  t. , ebisuzaki ,  t. , & umemura ,  m. 1990 , nature , 345 , 33 yoshikawa ,  k. & fukushige ,  t. 2005 , pasj , 57 , 849            , to integrate the system for one time unit as a function of the number of particles , @xmath0 . the triangle , square , and pentagon indicate the results with opening angles for the tree algorithm @xmath23 and @xmath24 , respectively . the star indicates the result for the hermite scheme . , width=340 ]    , defined in the text , as a function of the number of particles , @xmath0 . the triangle , square , and pentagon indicate those with opening angles for the tree algorithm @xmath61 and @xmath24 , respectively . the star indicates that for the hermite scheme . the thin dashed line indicates the peak performance for the hermite scheme . , width=340 ]       as a function of the number of particles , @xmath0 . the triangle , square , and pentagon indicate those with opening angles for the tree algorithm @xmath61 and @xmath24 , respectively . the star indicates that for the hermite scheme . , width=340 ]
Abstract:

Answer: we present an implementation of the hierarchical tree algorithm on the individual timestep algorithm ( the hermite scheme ) for collisional @xmath0-body simulations , running on grape-9 system , a special - purpose hardware accelerator for gravitational many - body simulations . such combination of the tree algorithm and the individual timestep algorithm was not easy on the previous grape system mainly because its memory addressing scheme was limited only to sequential access to a full set of particle data . the present grape-9 system has an indirect memory addressing unit and a particle memory large enough to store all particles data and also tree nodes data . the indirect memory addressing unit stores interaction lists for the tree algorithm , which is constructed on host computer , and , according to the interaction lists , force pipelines calculate only the interactions necessary . in our implementation , the interaction calculations are significantly reduced compared to direct @xmath1 summation in the original hermite scheme . for example , we can archive about a factor 30 of speedup ( equivalent to about 17 teraflops ) against the hermite scheme for a simulation of @xmath2 system , using hardware of a peak speed of 0.6 teraflops for the hermite scheme .