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 central parsec of our galaxy harbors different groups of young stars orbiting the central smbh . these include one or two stellar disks at @xmath0pc from sgr a@xmath1 made up of 56myr old stars @xcite , and a number of stars on highly eccentric orbits even closer , the so - called s stars @xcite . with orbital periods as small as 15 years , the s stars provide powerful constraints on the mass and size of the central dark object @xcite . while the origin of the stellar disk(s ) observed in the galactic center can be understood @xcite , the s stars show a puzzling combination of interesting features @xcite : they are very young ( @xmath2myr ) , with less than 0.01pc a factor of 10 closer to the smbh than the closest young population known ( the stellar disks ) , and move around sgr a@xmath1 on randomly oriented eccentric orbits .    a number of mechanisms to create the s stars have been proposed , but so far none of them was able to account for all their properties ( see also review by * ? ? ? * ) : it has been suggested that the s stars are not young but rejuvenated stars or exotic objects , but this is in conflict with their spectral properties , showing they are ordinary main - sequence b - type stars @xcite . in situ formation is excluded by the inexplicably high pressure required to overcome the strong tidal forces , and will also lead to a stellar disk @xcite , which can not account for the random orientation and high eccentricity of the orbits within the stars lifetime @xcite . the latter argument also argues against the suggested scenario of an infalling cluster of young stars @xcite , which would moreover require a very high mass @xcite , or an intermediate - mass black hole in its center for it to survive the tidal force of sgr a@xmath1 @xcite . exchange capture with compact remnants @xcite only works with a large cluster of black holes and fails to explain the youth and high initial eccentricities . the lense - thirring effect may account for the innermost star @xcite , s2 , being dragged out of the disk plane , but can not explain the orbital orientation of the longer - period s stars , nor the proximity of the orbit itself . interaction of stars in a central cluster with a fossil accretion disk can lead to high eccentricities and orbital decay @xcite if the disk mass is sufficiently high , but fails to explain the observed bias toward young b - type stars . furthermore , such an accretion disk is expected to fragment and form stars and should thus have been observed .    tidal disruption of binaries falling into the galactic center on highly eccentric orbits has been shown to be an effective mechanism to create s stars , but until now required very massive star clusters passing close to sgr a@xmath1 @xcite or the presence of a large enough population of compact massive perturbers in combination with many young stars on eccentric orbits @xcite .    in the following we show how the dynamical interaction of two stellar disks at larger radii , as they are observed in the galactic center , naturally leads to the creation of s stars . section [ sec : model ] describes the setup of our numerical simulations . in   [ sec : high - ecc ] , we present the outcome of these simulations , showing how the disk stars can eventually gain very high eccentricities . section [ sec : disrupt ] explains how binary stars on such eccentric orbits around the smbh get disrupted , leaving behind an s star and possibly a hypervelocity star . in  [ sec : discuss ] we discuss our results , showing how the mechanism suggested here explains all the observed properties of the s stars . an extensive analysis of young massive stars in the innermost parsec of our galaxy shows that almost all stars reside in two @xmath3myr old rotating disks  except for the s stars in the central arcsecond @xcite . while the clockwise disk is very distinct , there is some debate over the disk structure of the counterclockwise system @xcite . however , @xcite have shown that a single stellar disk can not explain the high inclinations and eccentricities observed . using direct @xmath4-body calculations as described below , we find that the distribution of eccentricities and inclinations of the observed young stars orbits is consistent with an initial configuration of two relatively flat disks , since the less massive disk is destroyed over time due to dynamical interaction ( see   [ sec : high - ecc ] ) .    to investigate the stellar dynamical evolution of two interacting stellar disks , we have computed a number of @xmath4-body integrations . the parameters of our model are chosen so as to best fit the observations by @xcite and are as follows : our model consists of a @xmath5 smbh and two flat circular stellar disks , mutually inclined at an angle of @xmath6 , with a surface density profile that scales with distance as @xmath7 . the disks have well - defined initial extents of @xmath8pc and @xmath9pc , respectively . we assume that the lower mass limit of the stars observed by @xcite is @xmath10 . for our models , we assume a disk age of 5myr , and hence an upper initial mass limit of @xmath11 for the observed disk stars @xcite . using the initial mass function @xmath12 as derived from the @xmath13-band luminosity function by @xcite with mass limits 1 and 120@xmath14 , this leads to a total initial stellar mass of @xmath15 to explain the 73 observed o stars . allowing for 30% binaries and another 10% for stars disrupted by the smbh or unobserved , and using a mass ratio of 2:1 for the disks , we have initial stellar masses of @xmath16 and @xmath17 , respectively . as a consequence , our disks initially consist of 535 and 267 multimass stars ( and binaries ) , respectively . the above parameters have proven to result in a distribution of eccentricities and semimajor axes as well as an inclination between the two disks after 5 myr which best match the observations . we used our bhint code @xcite for the integration , which has been developed specifically to calculate the dynamics of stars orbiting a smbh , and includes post - newtonian treatment up to order 2.5 to account for the effects of general relativity . furthermore , we have added the sse package by @xcite to account for the effects of stellar mass loss . the stars and binaries are modeled as point masses , and all binaries are assumed to be equal - mass . tidal disruption is not considered ; however , we assume that stars on very tight orbits ( @xmath18au ) are swallowed by the smbh . we switch on the post - newtonian terms for the central motion around the smbh in our bhint integrator @xcite for eccentricities @xmath19 , and switch it off when @xmath20 falls below 0.8 . a fully post - newtonian calculation shows that in the regime considered in this letter , relativistic effects are negligible for orbits with lower eccentricities within the considered timescale . we finally note that the detailed choice of parameters seems not to be too important . a number of models with variations in total mass , mass function , and disk age lead to comparable results . due to mutual torques , the two stellar disks described above precess about each other , with a frequency proportional to the other disk s mass : the precession frequency of a star orbiting a smbh of mass @xmath21 at radius @xmath22 at an inclination @xmath23 relative to a narrow disk of mass @xmath24 at a radius @xmath25 can be approximated as @xcite @xmath26    as can be seen , the precession frequency depends also on the distance of a star from the central smbh @xcite , and since both disks have a finite extent , stars at different central distances precess with different frequencies , thus warping the originally flat disks . in particular , the innermost stars quickly get inclined with respect to their disk of origin , and so start precessing about both disks . this leads to further thickening and deformation of the disks as well as large inclinations . figure [ fig : discevol1 ] depicts the evolution of the stellar disks shape and thickness . being exposed to the larger torque , the less massive disk is warped much faster , explaining why the observed thickness of the less massive disk in the center of our galaxy is substantially higher , making its present existence controversial @xcite . apart from precession , all stars are subject to the kozai mechanism @xcite : stars on highly inclined orbits relative to an axisymmetric perturbation periodically obtain small inclinations , and simultaneously gain high eccentricities . at their pericenter , stars on orbits with eccentricity @xmath27 pass the central smbh within @xmath28au , less than the s stars pericenter distances . the effects of general relativity causes eccentric orbits to precess , which may damp the kozai mechanism and prevent further growth in eccentricity @xcite . however , our simulations show that in the presence of very massive , nearby disks , relativistic effects do not prevent eccentricities as high as @xmath29 . figure [ fig : eccent ] shows the eccentricity evolution of some stars in our simulation , and the total fraction of stars that have reached eccentric orbits within a given time . a significant fraction of the stars obtain very high eccentricities . in contrast , a single stellar disk around a smbh evolves only very slowly : in a corresponding integration , no star obtained an orbital eccentricity above 0.8 within a comparable time span . this is mainly due to the lack of an external torque , and the lack of large inclinations required for the kozai mechanism to become effective . as most stars form in binaries @xcite , it is natural to assume that a significant fraction of stars in the stellar disks are followed by a companion . the orbits of these binaries around the smbh are subject to the same kozai mechanism as single stars , but a very close pericenter passage may break up the binary and leave one star on a tightly bound orbit around the smbh .    in a set of three - body integrations of a stellar binary orbiting a smbh , we have investigated the distribution of binary disruption remnants . as can be seen in figure [ fig : eccent ] , the eccentricity evolution due to the kozai effect is smooth . we calculated a set of models including a smbh , an intermediate - mass black hole ( imbh ) to drive the kozai , and a stellar binary on an initially circular orbit at high inclination with respect to the imbh s orbit in order to test the impact of a kozai oscillation on the binary s inner orbital parameters . figure [ fig : binkozai ] shows the time evolution of a 0.1au binary s both outer and inner eccentricity and semimajor axis . it can be seen that until binary disruption , only the outer eccentricity is changing . in particular , the binary does not widen . only during the last few pericenter passages does the inner eccentricity grow to a value of @xmath30 , but this does not significantly change the resulting semimajor axis of the bound member , nor the velocity of the ejected star . computations with wider binaries ( 1au ) yield similar results . we thus conclude that any tight binary undergoing kozai resonance reaches its point of disruption basically unperturbed . we have therefore computed a set of models where the pericenter distance of a binary orbiting a smbh is equal to the disruption distance ( fig . [ fig : disrupt ] ) . we find that the disruption remnants of binaries with separations larger than a few au orbit the central smbh on similar trajectories as the initial binary , and so can still be observed as disk stars . however , the disruption of a tight binary ( with separation of 1au or less ) forces one star onto a very tight and eccentric s - star - type orbit , while the other star gains a significant amount of energy and may even be ejected from the galaxy , thus becoming a progenitor of the hypervelocity stars observed in the galactic halo ( e.g. , * ? ? ? * ; * ? ? ? for example , disruption of binaries with a separation of 0.2au leaves one star on an orbit with semimajor axis consistent with those of the observed s stars , while the other star is ejected from the galaxy with a velocity of up to 1000kms@xmath31 in more than half of our integrations . binary disruption has been shown to be an effective mechanism to create the s stars @xcite . one way to create young close - passage binaries is by massive perturbers like star clusters or molecular clouds at distances of a few pc from the galactic center @xcite . our simulations show that these binaries are also a natural consequence of the interaction of the two stellar disks : using a binary fraction of 30% , 18 of the b - type binaries achieve pericenter distances below their tidal radius , assuming a binary separation of 0.1au . given that not much is known about the binary fraction in the galactic center , this is fully consistent with the observed number of 15 s stars . the mechanism described above naturally explains the reported mysteries of the s stars : the disruption of stellar binaries as tight as 0.11au can leave a single star orbiting the central smbh a factor of 10 closer than the original binary . this is just the distance ratio of s stars and stellar disks observed in the center of our galaxy . the orientation of the disks rapidly changes due to precession . furthermore , the change of orientation of a binary s orbit is a natural consequence of the kozai effect . once a binary is disrupted and leaves a single star on a much closer orbit , this orbit is decoupled from the originating stellar disk s motion , as both the precession frequency and the kozai oscillation period are functions of central distance , explaining the widely varying orbital orientations observed . the eccentricities created for the tightly bound stars after tidal disruption are as high as 0.99 . resonant relaxation @xcite and kozai interactions between the s stars are expected to turn these extremely high eccentricities into an isotropic distribution on a short timescale , just as is observed for the s stars .    while all of the s stars observed so far are main - sequence b - type stars , observations of the stellar disks in the galactic center mostly revealed ob giants and wolf - rayet stars . but even for the unusually flat mass function suggested @xcite , one would expect a factor of 5 more b - type main - sequence stars in the disks than o stars . most of these stars have not been observed yet , mainly due to crowding of the observed region , limited spatial resolution , and strong nebular emission @xcite . our computation assuming equal - mass binaries predicts only one s star above 20@xmath14 per 18 b - type s stars . the lack of observed o - type s stars is therefore not statistically significant . the number of expected massive s stars would drop further for random pairing of binary component masses , since disruption of a high mass ratio binary will only moderately change the orbit of the massive member . furthermore , the massive o stars might not form in as tight binaries as are required to produce s stars . hence , it is legitimate to assume that both s stars and disk stars formed in the same environment , thus providing a solution to the `` paradox of youth '' .    in this letter , we studied the evolution of two isolated stellar disks orbiting a smbh . in a forthcoming publication , we will analyze the effects of other contributions of a perturbing potential , such as a stellar cusp ( e.g. , * ? ? ? * ) , a massive circumnuclear disk @xcite , an intermediate - mass black hole and/or a star cluster ( probably irs 13 ; * ? ? ? preliminary results indicate , e.g. , that the presence of a cusp of stellar black holes ( e.g. , * ? ? ? * ) does not prevent the dynamical formation of s stars as described in this letter . we thank ladislav ubr for helpful discussions . this work was supported by the german research foundation ( dfg ) through the priority program 1177 `` witnesses of cosmic history : formation and evolution of black holes , galaxies and their environment . ''                  , f. , genzel , r. , alexander , t. , abuter , r. , paumard , t. , ott , t. , gilbert , a. , gillessen , s. , horrobin , m. , trippe , s. , bonnet , h. , dumas , c. , hubin , n. , kaufer , a. , kissler - patig , m. , monnet , g. , strbele , s. , szeifert , t. , eckart , a. , schdel , r. , & zucker , s. 2005 , , 628 , 246        , r. , schdel , r. , ott , t. , eisenhauer , f. , hofmann , r. , lehnert , m. , eckart , a. , alexander , t. , sternberg , a. , lenzen , r. , clnet , y. , lacombe , f. , rouan , d. , renzini , a. , & tacconi - garman , l.  e. 2003 , , 594 , 812                                        , t. , genzel , r. , martins , f. , nayakshin , s. , beloborodov , a.  m. , levin , y. , trippe , s. , eisenhauer , f. , ott , t. , gillessen , s. , abuter , r. , cuadra , j. , alexander , t. , & sternberg , a. 2006 , , 643 , 1011
Abstract:

Answer: over the last 15 years , around a hundred very young stars have been observed in the central parsec of our galaxy . while the presence of young stars forming one or two stellar disks at @xmath0pc from the supermassive black hole ( smbh ) can be understood through star formation in accretion disks , the origin of the s stars observed a factor of 10 closer to the smbh has remained a major puzzle . here we show the s stars to be a natural consequence of dynamical interaction of two stellar disks at larger radii . due to precession and kozai interaction , individual stars achieve extremely high eccentricities at random orientation . stellar binaries on such eccentric orbits are disrupted due to close passages near the smbh , leaving behind a single s star on a much tighter orbit . the remaining star may be ejected from the vicinity of the smbh , thus simultaneously providing an explanation for the observed hypervelocity stars in the milky way halo .


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 heisenberg model@xcite expresses the energy dependence of localized electron spins by means of the exchange interaction . this assumption creates a starting point for an analysis on magnetic states in insulators . there are a few well founded results for this model . for instance , one - dimensional systems , the bethe ansatz and bosonization solutions provide a set of exact solutions of the model@xcite . rigorous results at zero temperature have been obtained for two and three hypercubic lattices@xcite . by numerical approach , physical properties are exact results for finite systems and have been conjectured as an exact result for the properties of infinite systems@xcite . more generallly , the study of the heisenberg model is made using methods of approximated solutions . however , these approaches certainly contain partial information that can differ from the exact solution . another natural possibility is to introduce a new model that makes the analytical analysis easier or that can evidently be solved exactly . it is important , however , that the model remains appropriate for describing magnetic systems . this picture can be emphasized considering the ising model@xcite which is the most important simplification of the heisenberg model . the crucial assumption is that the ising hamiltonian is explicitly diagonal , but nontrivial from the point of view of statistical mechanics . the ising model has played a key role in the development of the theory of interacting systems . it has a phase transition at finite temperature that can be worked out with mathematical rigor and yields an exact solution in two dimensions . although it has only been considered a theoretical model for many years , experiments have brougt forth evidence of its applicability for real magnetic materials ( see @xcite and references therein ) . it should be noted that another virtue of the model is its correspondence with systems such as lattice gases , binary alloys and neural networks among othershuang .    in conductors , the study of magnetism is more complex because the electrons of magnetic states are often not localized . the basic model to describe the interaction between electrons having translational degrees of freedom is the hubbard model@xcite .    suggested in 1963 , the hubbard model also much like the heisenberg model remains with many of its basic features in debate . an important problem to be answered is to establish which limits within the parameters of the hubbard model that favor ferromagnetism when neither the electronic interaction nor the electronic itinerance favors this situation . many believe that ferromagnetissm can only be observed when there are more than one electronic orbital per site . in the simple hubbard model when only one orbital per site is present the induced antiferromagnetism due to the electronic itinerant term disfavors ferromagnetism . however , in the extreme limits of electronic repulsion and in special lattice cases , we can find ferromagnetism in this model@xcite . indeed , the small number of accurate solutions ( one - dimensional systems , ground state theorems and numerical approach)@xcite and the difficulty of studying the hubbard model indicate that the understanding of this problem remains vague . for example , we do not have an accurate solution for the transition metals .    knowing the difficulties in dealing with the hubbard model and the importance of studying the effect of the electronic itinerance term , i will shortly introduce a model that remains with the essential features of the hubbard model , but simplifies it using a classical situation for the hopping . despite the simplification , this model must be capable of describing the magnetism of materials with itinerant electrons . this paper is organized as follows . in sec . [ sec : model ] the model is described . the one - dimensional ring considering nearest - neighbor hopping is solved in sec . [ sec : one - dimensional ] . [ sec : mean - field ] presents the mean - field solution for the model . [ sec : conclusion ] is devoted to conclusions . the starting point for the present model is the hubbard model , which has the following hamiltonian @xmath0where @xmath1 , @xmath2 are the creation and annihilation operators for electrons of spin @xmath3 at site @xmath4 . the density of electrons is denoted by @xmath5 where @xmath6 is the number of lattice sites , and @xmath7 is the number operator for electrons of spin @xmath3 at site @xmath4 . the first term of the hamiltonian corresponds to the hopping of the electrons between sites @xmath4 and @xmath8 , and @xmath9 is the hopping integral representing the overlap of electron wave functions . usually , the overlapping is assumed only over nearest neighbour sites . the second term represents the on - site coulomb interaction @xmath10 between the electrons . the crucial ingredient for the construction of the model is to consider the energy associated to the electron itinerance as being equal to the total number of possible electron jumps . hence , it is easy to see that the hamiltonian can be defined as    @xmath11    in this case , the hopping term commutes with the coulombic term . the states of the number operator are eigenstates of the hamiltonian . as in the ising model , the hamiltonian is explicitly diagonal and the model is much easier to study compared to the hubbard model . the hamiltonian of eq . ( @xmath12 ) can be seen as an approximation of the extended hubbard model , constructed by the inclusion of a third term in the simple hubbard model ( eq . @xmath13 ) which describes the off - site coulomb repulsion between electrons@xcite,@xcite,@xcite . the off - site interaction has been appropriate for compounds that exibit near - neighbors interaction compared to the bandwidth , in particular quasi one - dimensional organic conductors@xcite . the hamiltonian of eq . ( @xmath12 ) is the zero - band - width limit of the extended hubbard model assuming off - site interaction only between electrons of equal spins . this assumption can be justified  considering that the electron repulsion of different spins has its main contribution on the same site , while electron repulsion of equal spins contributes only to different sites by the pauli exclusion principle . hubbard showed that the zero - band - width limit of the one - dimensional extended hubbard model can present an interpretation of the optical spectra of many tetracyanoquinodimethane salts . indeed , some features of finite bandwidth may well survive in the present model if the interaction effects are of dominant importance@xcite . these considerations indicate to be reasonable the application of the model represented by the hamiltonian of eq . ( 2 ) to the quasi one - dimensional organic conductors as well as to other materials with the same electronic structure . since we wish a model capable of describing the magnetic properties of materials with itinerant electrons , it should be important to observe that relevant symmetry transformations of the hubbard model are preserved , such as the u(1 ) charge symmetry and the particle - hole transformation in the case of bipartite lattice . notice that this has an important consequence : the condition for the chemical potential @xmath14 for the half - filled band , independently of the temperature , is identical to that of hubbard model . it is also crucial to observe the ground state configurations of the new model . here , only general features will be presented . detailed analyses depend on specific works using this model . it is easy to see that at half - filled band the two nel antiferromagnetic configurations are ground states of the new hamiltonian for @xmath15 in non frustrated lattices . when the number of electrons is smaller than half of the number of lattice sites the ground state becomes highly degenerate and the ferromagnetic configuration has also the lower energy . the model favors the antiferromagnetic order , but for low or high electron densities there is a tendency towards the coexistence of ferromagnetism with antiferromagnetism . as will be observed in the one - dimensional case , for the half - filled band the antiferromagnetic ordered phase exists only at zero temperature . with finite temperature the ordered phase does not exist .    in the square lattice , for the half - filled band case , the ground state is antiferromagnetic , but at finite temperature(@xmath16 ) the solution has not yet been found . for this same lattice , the hubbard model is not solvable and even its ground state is not known . using the monte - carlo approach , hirsch@xcite predicts antiferromagnetic ordered ground state only at @xmath17 . i believe that both the new model and the hubbard model yield a similar result to the relation between the ising and heisenberg models . with finite temperature , both the new and the ising models have an ordered phase while the hubbard and the heisenberg models do not@xcite . one sees that in this case the hamiltonian ( 2 ) can be written as @xmath18 , where @xmath19let us take @xmath20 over nearest neighbor sites , otherwise @xmath21 . the thermodynamic properties of the model can be obtained using the transfer matrix method . the grand - canonical partition function , at temperature @xmath16 , is written as @xmath22 , where @xmath23 are eigenvalues of the transfer matrix    @xmath24    where rows and columns are labelled by @xmath25 and @xmath26 , respectively . the matrix @xmath27 is 4x4 and will be numerically diagonalized . the largest eigenvalues @xmath28of @xmath27 are found . considering the thermodynamic limit ( @xmath29 ) , the  free energy is @xmath30 , and the thermodynamic properties follow from properly differentiating the free energy . for the half - filled band the analytical expression for @xmath28 is    @xmath31    where @xmath32 $ ] and @xmath33 $ ] . figure 1 presents the specific heat @xmath34 as a function of temperature for typical values of @xmath35 for the half - filled band . for @xmath36 it has a peak . whose value increases with @xmath37 until the maximum @xmath38 for @xmath39 and decreases for larger values of @xmath37 . for @xmath40 the peak splits in two , which reflects a rearrangement of the fermionic structure in the system@xcite . the picture presented is analogous to those ones of the one - dimensional half - filled band hubbard@xcite and falikov - kimball models@xcite . however , the critical values of @xmath37 which causes this splitting are different : @xmath41 for the falikov - kimball model , @xmath42 for the hubbard model and @xmath43 for the present model . for @xmath44 the second peak ( high - temperature peak ) gets further from the first peak as @xmath45 and then the first peak yields the specific heat with the maximum value of @xmath46 which is exactly the specific heat of the ising model@xcite . figure 2 shows the temperature dependence of the magnetic susceptibility @xmath47 for the several choices of @xmath37 in the half - filled case . the @xmath47 has a maximum and vanishes at zero temperature . for large @xmath37 the magnetic susceptibility is well described by the ising model as was evidenced by the specific heat in fig . finally , it is obtained without further difficulties no phase transition at any finite @xmath16 and no long - range order for all @xmath48 .     of infinite chain versus temperature for typical values of @xmath37 : i ) @xmath41 ; ii ) @xmath49 ; iii ) @xmath50 ; iv ) @xmath51 ; v ) @xmath52.,width=302 ]     : i ) @xmath41 ; ii ) @xmath49 ; iii ) @xmath53.,width=302 ] having obtained solutions for the one - dimensional case one presents below the mean - field solution . the self - consistent equations for translationally invariant systems are given by    @xmath54    and    @xmath55    where @xmath56 and @xmath57 are the mean values of the electron number ( @xmath58 ) and local magnetization ( @xmath59 ) , respectively , on a lattice where each site has @xmath60 nearest neighbours . in the half - filled band @xmath61 , which implies that @xmath14 independently of the temperature and it is easy to see that @xmath62 ^{-1}$ ] . the dependence of the magnetization on @xmath37 is considered . the critical temperature @xmath63 is ( assuming @xmath64 , @xmath65 )    @xmath66    the phase diagram is shown in fig . 3 . below @xmath63 there is a spontaneous magnetization . the inset in fig . 3 shows the temperature dependence of the magnetization . the magnetization increases with decreasing temperature and attains the saturate value at zero temperature . mean - field solutions suggest that the  magnetic order phase increases if @xmath37 increases in agreement with the hubbard model predictions@xcite . values.,width=302 ]       notice that for the @xmath67 limit the model is reduced to a lattice gas of two types of atoms . this means that the present model can also be interpreted as a generalization of the ising model which is exactly the @xmath44 case . in summary , a classical lattice model was introduced to describe the magnetism of itinerant electron systems . it is much easier to analyze than the hubbard model . i have illustrated some features of its magnetic properties . these general results clearly show that the presented model can be viewed as a classical hubbard model in the same way as the ising model can be considered with relation to the quantum heisenberg model . naturally , it is not able to capture the quantum frustration of the hubbard model much like the ising model can not capture the quantum frustration of the heisenberg model . however , i believe that it is a promising help to qualitatively clarify our understanding on the nature of magnetism . in addition , the proposed model is a generalization of the ising model . the possibility of generalizing the onsager solution for the square lattice case would be useful in future phase transition research . i thank m. e. de souza and a. reikdal for reading the manuscript . this work was financially supported by cnpq and fap - se . k. harada , n. kawashima , and m. troyer . lett . * 90 * , 117203 ( 2003 ) ; e. manousakis , rev . phys . * 63 * , 1 ( 1991 ) ; t. kennedy , e. h. lieb , and s. shastry , j. stat . phys . * 53 * , 1019 ( 1988 ) ; k. kubo and t. kishi , phys . rev . lett . * 61 * , 2585 ( 1988 ) ; f. j. dyson , e. h. lieb , and b. simon , j. stat . phys . * 18 * , 335 ( 1978 ) . j. c. bonner and m. e. fisher , phys . rev . * 135 * , a640 ( 1964 ) ; k. lefmann , c. rischel , eur . j. b * 21 * , 313 ( 2001 ) ; a. a. cuccoli , t. roscilde , v. tognetti , r. vaia , and p. verrucchi , phys . b * 67 * , 104414 ( 2003 ) .              e. h. lieb and f. y. wu . lett . * 20 * , 1445 ( 1968 ) ; r. strack and d. vollhardt , phys . . lett . * 72 * , 2425 ( 1994 ) ; h. tasaki , phys . . lett . * 69 * , 1608 ( 1992 ) ; s. a. trugman , phys . b * 42 * , 6612 ( 1990 ) ; e. h. lieb , phys . rev . lett . * 62 * , 1201 ( 1989 ) . gu , s .- s . deng , y. -q . li , and h. -q . lin , phys . . lett . * 93 * , 086402 ( 2004 ) ; y. z. zang , phys . 92 * , 246404 ( 2004 ) ; s. onari , r. arita , k. kuroki , h. aoki , phys . b * 70 * , 094523 ( 2004 ) .
Abstract:

Answer: a new lattice model of interacting electrons is presented . it can be viewed as a classical hubbard model in which the energy associated to electron itinerance is proportional to the total number of possible electron jumps . symmetry properties of the hubbard model are preserved . in the half - filled band with strong interaction the model becomes the ising model . the main features of the magnetic behavior of the model in the one - dimensional and mean - field cases are studied . classical spin model , lattice model , magnetism


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: supersymmetry ( susy ) is one of the most fascinating models beyond the standard model ( sm ) , which can simultaneously solve the naturalness problem , explain the cosmic dark matter and achieve the gauge coupling unification . searching for the susy particles is an important task for lhc especially after the discovery of higgs . with the recent 13  fb@xmath3 dataset of lhc , a gluino and a stop lighter than 1.8 tev and 900 gev have been excluded out respectively if the mass splitting between the stop and the lsp is large enough , and the mass of an electroweakino has been excluded to 0.4 - 1 tev depending on the decay mode . all the results push susy scale much higher than before , which further challenges our understanding on the naturalness problem .    among those tens of super - particles , the top partner - stop , plays the most important role to understand the naturalness in susy models . in the framework of natural supersymmetry , a light stop and higgsino are usually predicted  @xcite . lots of works have been done on searching for such light stops and higgsinos  @xcite . recently atlas reported some excesses in the stop search with @xmath4 channel  @xcite . seven signal regions are defined and three of them observed more than @xmath2 excesses , among which the most significant one could reach @xmath0 . although no significant excess is reported by cms in the same channel  @xcite , they observe mild excesses in several signal regions . even in the @xmath5 lepton channel searches @xcite , atlas and cms both show small excesses with a lower significance . if they are true signals it is likely due to new particles .    in this paper , we study the interpretation of the excesses as light stop pair production in the minimal supersymmetric standard model  ( mssm ) . although atlas and cms have done the similar analyses based on the simplified models , the results might be changed in a concrete model of mssm . we show that the simplified models where stop decays into a higgsino or a bino are not favored by other stop searches while extended models could provide a better fit to the excess . this paper is organized as follows . we first briefly overview the atlas @xmath4 search and the other constraints at the lhc in section  [ sec : stop ] . in section  [ sec : model ] we interpret the excess by a light stop pair production in the mssm . in section  [ sec : darkmatter ] , we discuss the dark matter constraints , such as the relic abundance and the direct detection constraints . we draw our conclusions in section [ sec : conclusion ] . the atlas @xmath4 search mainly aims at stop pair production followed by stop decay into top and neutralino , where one of the resulting top decay leptonically . although more details are found in  @xcite we summarize some of the selection cuts and the results in table  [ tab : atlas ] . seven signal regions are defined for the searches . among them , the signal regions of dm_low , bc2x_diag and sr1 observe @xmath0 , @xmath7 , @xmath8 excesses , respectively . interestingly , the sr1 signal region was already defined before and an excess around @xmath9 had been reported even with the 3.2  fb@xmath3 data  @xcite . [ cols="^,^,^,^,^,^,^,^",options="header " , ]      there are other stop searches based on 0 lepton , 1 lepton , and 2 leptons both from atlas and cms , which constrain the stop parameter space . we have found the 1 lepton search constraints from cms  @xcite gives more or less similar to that of the atlas , therefore , considering the signal regions with no excess in atlas 1 lepton search is enough . the two lepton search channels do not set strong constraints  @xcite , and we did not consider them in this paper . the cms 0 lepton search channels set slightly stronger bound than those of atlas , therefore , we consider the following two hadronic stop search channels at cms to set the 95% cl exclusion . the cms boosted top search requires top tagging and is especially sensitive to the stop with a large mass splitting with lsp , where boosted tops are expected in the final state . the cms hadronic stop search aims at more conventional topologies from stop decay . in the following , the atlas 1 lepton search is denoted as atlas_1l , the cms hadronic stop search  @xcite is denoted as cms_hadronic and cms boosted top search  @xcite is denoted as cms_boosted . in this paper we interpret the excess from a light stop pair production and their decays . we assume the lightest stop dominantly consists of the right handed stop  ( @xmath10 ) . the @xmath11 usually is accompanied by @xmath12 which is close in mass , and it is constrained strongly with a bound up to about 1 tev  @xcite . we first consider a simple 1-step decay from the right - handed stops . depending on the dominant component of the lsp , the constraints are modified as the decay modes change . we scan the parameter space to figure out the corresponding @xmath2 favored region in dm_low . our simulation for atlas_1l and cms_hadronic searches are based on madgraph + pythia  @xcite , and the generated signal events are passed to delphes3  @xcite for a detector simulation . for recasting cms_boosted search we simply compare the cross sections rescaled with @xmath13 with the upper bound of the cross section given in the cms paper assuming the best sensitivity is from two boosted top tagged events . we first consider the simplest model : the lightest supersymmetric particle  ( lsp ) is bino and the next - lightest supersymmetric particle  ( nlsp ) is @xmath10 . as shown in fig  [ decay_bino ] , the only decay channel is @xmath14 in this model , which is frequently assumed for many analyses in atlas and cms . in fig  [ bino ] we show the @xmath2 favored region by the excess in dm_low and the 95% c.l . exclusion contours on the @xmath15 plane from atlas_1l , cms_hadronic , and cms_boosted searches . the strongest limit is from cms_boosted search because the tops in the final state tend to be boosted . we find all the @xmath16 favored region are excluded out by cms_boosted and almost all the @xmath2-region is excluded .      to ease the tension from cms_boosted stop search , we may consider a model where the stop has other decay channels than the direct decay into top and lsp . as an example , we consider another simplified model where lsp is higgsino . for the higgsino lsp case , the higgsinos ( @xmath17 ) and charged higgsino ( @xmath18 ) is naturally degenerate . as the typical mass difference is @xmath19 gev , we can not observe the decay products from those particles essentially and all particles behaves like lsp in terms of the collider signature at the lhc . the branching ratios are @xmath20 when the phase space suppression due to the top mass is negligible . thus , the contribution of the two boosted tops is reduced by a factor @xmath21 . since cms_boosted stop search relays on two boosted tops in the final states , the constraint becomes much weaker . on the other hand , the @xmath22 signal can originate from events with one stop decay into a top and the other into a bottom , thus , there are no strong reduction factor . in fig  [ higgsino ] we show the signal preferred region and several exclusion contours . in the plot we fixed @xmath23 . note the branching ratios of stop decay do not change a lot when we vary the @xmath24 because higgsino - stop coupling only comes from the yukawa . the gray region is limited by direct chargino search from lep@xcite . the cms_boosted stop search becomes significantly weaker , and there appears a large region not excluded by either atlas_1l and cms_boosted but in the @xmath2 favored region . however , the cms_hadronic constraints are still strong enough to exclude the whole @xmath2 signal favored region . it is because reducing @xmath25 also reduce the number of 1 lepton signals , while the conventional 0 lepton signals are not reduced and it results in a similar sensitivity to the bino lsp case .     favored region and the excluded region from the bino lsp model . , scaledwidth=40.0% ]     favored region and the excluded region from the higgsino lsp model . , scaledwidth=40.0% ] although by reducing @xmath26 the tension between cms_boosted and atlas_1l searches could be eased , it also reduces the signal events and makes the conventional hadronic stop search relatively more effective . to avoid this situation , keeping more signal events while reducing top branching ratio is necessary , therefore , it is preferable to find a way to make the @xmath27 also contribute to enhance the lepton signals . we consider here a model where nlsp is higgsino and lsp is bino . since the chargino will decay into @xmath28 plus neutralino , one lepton could come from the @xmath28 decay through a two step decay . we set the mass difference between stop and higgsino 150  gev to forbid @xmath29 decay to make the discussion simpler . we have checked that opening @xmath30 mode also gives similar final results .        in fig [ fig : higgsinobino ] , we show our simulation results recasting the atlas_1l , cms_hadronic , and cms_boosted searches on the @xmath15 plane assuming @xmath31  gev . large parameter region in the higgsino + bino lsp model satisfying @xmath32  gev with @xmath33  gev or @xmath34  gev@xmath35  gev with @xmath36  gev is found consistent within @xmath2 to all constraints although @xmath16 favored region is still excluded by both atlas_1l and cms_hadronic . favored region and the excluded region in higgsino+bino lsp model . we fix @xmath37 gev and @xmath38 . left side of the lux line has been excluded by the dark matter direct detection experiments.,scaledwidth=40.0% ] atlas_1l also provides the @xmath39 and @xmath40 distributions in dm_low . we have selected three benchmark points , which are indicated with crosses in fig  [ fig : higgsinobino ] , and show the expected distributions . the background distributions we just take from the atlas plots . the benchmark points with stop - bino mass 650 - 350 gev , 750 - 300 gev , 800 - 200 gev predict the signal events in dm_low to be 9.4 ( @xmath41 ) , 9.8 ( @xmath42 ) and 8.3 ( @xmath43 ) , respectively . the numbers in the parentheses indicate the statistical deviations assuming the corresponding signal injection , to be compared with @xmath44 of no signal assumption . although all the benchmark points are consistent within @xmath45 based on the total number of events in dm_low the predicted distributions are different . we find compressed spectrum is slightly preferred as the overflow bin does nt contain much events . in the future these distributions would be important to distinguish between models .     and @xmath40 distribution for the three benchmark points in the consistent region of the higgsino + bino lsp model . the numbers in legend are the stop and the lsp masses . the right most bin contains overflow events.,title="fig:",scaledwidth=40.0% ]   and @xmath40 distribution for the three benchmark points in the consistent region of the higgsino + bino lsp model . the numbers in legend are the stop and the lsp masses . the right most bin contains overflow events.,title="fig:",scaledwidth=40.0% ]    although we only consider the higgsino - bino case ( we denote the case with higgsino nlsp and bino lsp as higgsino - bino , etc . in the following ) , we would like give some comments on other possibilities . since the coupling between stop  ( @xmath10 ) and wino is suppressed by the neutralino mixing , the wino - bino ( wino - higgsino ) cases are essentially reduced to the bino ( higgsino ) lsp case . with the large l - r mixing in stop sector we can tune the relative branching ratios by the stop mixing angle and the collider signature of the wino - bino case could be similar to the higgsino - bino model . three remaining possibilities are bino - higgsino , bino - wino , and higgsino - wino cases . for the bino - higgsino case , the stop dominantly decays into a higgsino , as stop couples higgsinos thorough the top yukawa coupling , which is much stronger than bino through the gauge coupling , and due to the phase space suppression , therefore , it is again reduced to the higgsino lsp case . the higgsino - wino case is similar to the higgsino - bino case but the higgsino decay pattern is more complex , and branching ratio into @xmath28 would be suppressed . the bino - wino case might be an interesting possibility but constrained by the multi - lepton searches . on the other hand , direct detection constraint must be significantly weak for this case whereas indirect searches might be dangerous . also notice that the long lived charged particle searches at 8 tev  lhc already set the lower bound of the wino mass around 270 gev  @xcite , and the allowed parameter space would be highly constrained . in the previous section , we study the case where bino is the lsp . it might overclose the universe because the pair annihilation cross section of bino like dark matter is rather small . an interesting possibility is to resort to the slepton co - annihilation . if a slepton is nearly degenerate with the bino lsp @xmath46  gev  @xcite it would help to reduce the relic abundance to the observed one . one may wonder inserting slepton in the spectrum will affect the decay mode of the chargino . as the coupling between slepton and higgsino is proportional to the corresponding lepton yukawa coupling only a light stau with a large @xmath24 could affect the branching ratios . note that the chargino decaying into neutralino plus @xmath28 is dominated by the longitudinal mode of @xmath28 ( from @xmath47-@xmath48-@xmath49 interaction ) , and not suppressed by the higgsino - bino mixing . the stau decaying channel only becomes sizable when the @xmath50 is enhanced by very large @xmath24 . taking @xmath38 as an example , we find in our @xmath2 favored region , the stau decay mode at most has a branching ratio of @xmath51 . if we reduce the @xmath52 , this decay branching ratio only takes 3% , so our fig  [ fig : higgsinobino ] is essentially unchanged . moreover , the @xmath53 emitted from the @xmath54 decay is not detectable at the collider due to the small mass difference .    the dark matter direct search also constrains the higgsino - bino model due to the mixing between higgsino and bino . in fig [ fig : higgsinobino ] we draw the sensitivity of the latest lux results  @xcite the parameter space on the left side of the line has been excluded . in the plot we fix @xmath37 gev and @xmath38 , sign@xmath55 , and the higgs mass is also tuned to be 125 gev . we find one of our benchmark points has been excluded out and one just lives the boundary of exclusion line . therefore the future dark matter direct search will be important to test this scenario if the signal is confirmed at the lhc . note that in the case of sign@xmath56 , there exists a blind spot region for dark matter direct search when @xmath57 is satisfied@xcite . we may also assume r - parity is violated and bino could decay into sm particles with the life time long enough in a collider time scale . in such a case , the constraints discussed in this section are not applicable . in this paper , we try to interpret the @xmath0 excess reported in atlas @xmath1 channel by a light stop pair production in the mssm . after considering the consistency with other stop searches , we find : ( 1 ) simple models where stop decays into a higgsino or a bino are not favored by the other search channels ( 2 ) an extension of them , higgsino - bino lsp , can explain the data at @xmath2 level without conflicting with the other search channels . it is remarkable that the susy spectrum in the surviving stop decay scenario for the possible excess , that is a light stop and a light higgsino , is nothing but what we expect in a natural susy . * acknowledgements * the authors wish to thank dr . lei wu for valuable discussions and for the technical helps . cch , mt , mmn and tty are supported by world premier international research center initiative ( wpi initiative ) , mext , japan . mmn is also supported by grant - in - aid for scientific research nos . jp16h06492 , jp16h03991 and jp26287039 . mt is supported by grant - in - aid for scientific research numbers jp16h03991 and jp16h02176 . tty is supported by jsps grants - in - aid for scientific research nos . jp26287039 , jp26104009 and jp16h02176 . r.  kitano and y.  nomura , phys . d * 73 * ( 2006 ) 095004 ; l.  j.  hall , d.  pinner and j.  t.  ruderman , jhep * 1204 * ( 2012 ) 131 ; m.  papucci , j.  t.  ruderman and a.  weiler , jhep * 1209 * ( 2012 ) 035 ; c.  brust _ et al , _ jhep * 1203 * ( 2012 ) 103 ; h.  abe , t.  kobayashi and y.  omura , phys . d * 76 * ( 2007 ) 015002 ; h.  baer _ et al , _ phys . d * 89 * ( 2014 ) no.11 , 115019 ; z.  kang , j.  li and t.  li , jhep * 1211 * ( 2012 ) 024 ; c.  han , f.  wang and j.  m.  yang , jhep * 1311 * ( 2013 ) 197 . m.  drees , m.  hanussek and j.  s.  kim , phys .  rev . d * 86 * ( 2012 ) 035024 ; j.  cao _ et al , _ jhep * 1211 * ( 2012 ) 039 ; c.  han _ et al , _ jhep * 1310 * ( 2013 ) 216 ; c.  han _ et al , _ jhep * 1402 * ( 2014 ) 049 ; z.  han , g.  d.  kribs , a.  martin and a.  menon , phys .  rev . d * 89 * ( 2014 ) no.7 , 075007 ; c.  han , arxiv:1409.7000 [ hep - ph ] ; h.  baer , a.  mustafayev and x.  tata , phys .  rev . d * 90 * ( 2014 ) no.11 , 115007 ; c.  han _ et al , _ jhep * 1504 * ( 2015 ) 132 ; k.  i.  hikasa _ et al , _ phys . d * 93 * ( 2016 ) no.3 , 035003 ; a.  kobakhidze _ et al , _ phys . b * 755 * ( 2016 ) 76 ; c.  han _ et al , _ arxiv:1609.02361 [ hep - ph ] ; t.  plehn _ et al , _ jhep * 1010 * , 078 ( 2010 ) ; t.  plehn , m.  spannowsky and m.  takeuchi , jhep * 1105 * , 135 ( 2011 ) ; t.  plehn , m.  spannowsky and m.  takeuchi , jhep * 1208 * , 091 ( 2012 ) ; m.  r.  buckley , t.  plehn and m.  takeuchi , jhep * 1308 * , 086 ( 2013 ) ; d.  gonalves , k.  sakurai and m.  takeuchi , arxiv:1604.03938 [ hep - ph ] . cms collaboration [ cms collaboration ] , cms - pas - sus-16 - 029 . cms collaboration [ cms collaboration ] , cms - pas - sus-16 - 030 . the atlas collaboration [ atlas collaboration ] , atlas - conf-2016 - 077 . m.  aaboud _ et al . _ [ atlas collaboration ] , arxiv:1606.03903 [ hep - ex ] . j.  alwall _ et al . _ , jhep * 1407 * ( 2014 ) 079 doi:10.1007/jhep07(2014)079 [ arxiv:1405.0301 [ hep - ph ] ] . t.  sjostrand , s.  mrenna and p.  z.  skands , jhep * 0605 * ( 2006 ) 026 doi:10.1088/1126 - 6708/2006/05/026 [ hep - ph/0603175 ] . g.  abbiendi _ et al . _ [ opal collaboration ] , eur . j.  c * 29 * ( 2003 ) 479 ; j.  abdallah _ et al . _ [ delphi collaboration ] , eur . j.  c * 31 * ( 2003 ) 421 ; g.  abbiendi _ et al . _ [ opal collaboration ] , eur .  phys . j.  c * 35 * ( 2004 ) 1 . g.  aad _ et al . _ [ atlas collaboration ] , phys . d * 88 * ( 2013 ) no.11 , 112006 doi:10.1103/physrevd.88.112006 [ arxiv:1310.3675 [ hep - ex ] ] . m.  citron , j.  ellis , f.  luo , j.  marrouche , k.  a.  olive and k.  j.  de vries , phys . d * 87 * , no . 3 , 036012 ( 2013 ) doi:10.1103/physrevd.87.036012 [ arxiv:1212.2886 [ hep - ph ] ] . d.  s.  akerib _ et al . _ [ lux collaboration ] , phys . * 116 * ( 2016 ) no.16 , 161302 doi:10.1103/physrevlett.116.161302 [ arxiv:1602.03489 [ hep - ex ] ] . c.  cheung , l.  j.  hall , d.  pinner and j.  t.  ruderman , jhep * 1305 * ( 2013 ) 100 doi:10.1007/jhep05(2013)100 [ arxiv:1211.4873 [ hep - ph ] ] . p.  huang and c.  e.  m.  wagner , phys . d * 90 * ( 2014 ) no.1 , 015018 doi:10.1103/physrevd.90.015018 [ arxiv:1404.0392 [ hep - ph ] ] .
Abstract:

Answer: recently atlas reported a @xmath0 excess in the stop search with @xmath1 channel . we try to interpret the signal by a light stop pair production in the mssm . we find : ( 1 ) simple models where stop decays into a higgsino or a bino are not favored . ( 2 ) an extension of them can explain the data at @xmath2 level without conflicting with the other search channels . a surviving possibility includes a light stop and a light higgsino , which is expected in a natural susy scenario .


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 random phase approximation ( rpa ) for electron correlation was originally developed in the context of the many - body perturbation treatment of the uniform electron gas @xcite and later reformulated in the framework of density functional theory ( dft ) via the adiabatic - connection fluctuation - dissipation theorem @xcite . in the last years , rpa correlation witnessed an increasing interest in quantum chemistry @xcite , due to its ability to describe with good accuracy a rather large variety of reaction and interaction energies . in this context rpa is generally employed as a post kohn - sham ( ks ) approach where the exact exchange ( exx ) and the correlation energies are evaluated using the eigenvalues and eigenfunctions from a ks calculation with conventional exchange - correlation ( xc ) functionals . exx ks orbitals @xcite have been also employed via the optimized effective potential ( oep ) method , while for a fully variational approach self - consistent exx - rpa orbitals are required @xcite . in all cases practical applications of the method are often hindered by the the very slow convergence of the rpa - correlation energy with the basis set dimension @xcite , which is related to the electron cusp problem @xcite , in analogy with wavefunction methods . the cusp problem is strongly reduced for range - separated rpa @xcite , whereas for full - range rpa approaches extrapolation to the complete basis set ( cbs ) must be employed . basis set extrapolation to the cbs limit has been longley studied for post hartree - fock correlated methods , with the aim of reducing as much as possible the basis set truncation error without resorting to the brute force employment of extremely large basis sets . numerical work @xcite on various correlation treatments of the helium atom has shown that the basis set incompleteness error is approximately proportional to the inverse third power of the maximum angular momentum , provided that all angular momenta are radially saturated . however , the ideal conditions of fully saturated angular momenta in an atomic system are rather distant from those encountered in real applications on molecules with finite basis sets . nevertheless , the correlation energy has been empirically observed to converge regularly enough to allow the development of various accurate extrapolation schemes which have been assessed in several studies in literature @xcite . on the other hand , no such study has been performed for rpa correlation energies .    in this paper we aim at filling this gap and shed light on the usefulness and limits of different extrapolation schemes for rpa correlation energies . in this respect the goal of the present work is to : ( 1 ) create a set of benchmark cbs rpa absolute correlation energies to be used as reference for future assessments ; ( 2 ) validate some extrapolation schemes , usually applied in the context of coupled cluster theory , for the calculation of cbs rpa absolute energies , including semiempirical techniques to approach the extrapolation from relatively small basis sets ; ( 3 ) test the convergence behavior of different basis sets and the reliability of the extrapolations schemes for atomization energies , which exploit a significant error cancellation effect . absolute correlation energies will be considered , as in many corresponding studies of correlated methods , because these converge regularly with the basis set and provide an ideal quantity for the study of extrapolation schemes . in this paper we focus on light elements because they provide a simpler convergence of the correlation energy and can be more easily investigated with large basis sets . rpa calculations were performed using the eigenvalues and the orbitals from dft calculations based on the perdew - burke - ernzerhof ( pbe ) exchange - correlation functional @xcite . in all calculations dunning s correlation - consistent basis sets @xcite ( cc - pv@xmath0z , @xmath1 ) augmented with core and core - valence basis functions @xcite were employed . these are hereafter denoted as v@xmath0zcv with @xmath0=4 , 5 , 6 , and 7 . augmented diffuse basis functions were not considered in the present work since they were shown to slow down rpa basis set convergence without bringing substantial benefits in most cases @xcite . also , no basis set superposition error correction was considered for atomization energies . for molecules experimental geometries were considered @xcite . all calculations were performed with the turbomole program package @xcite using the implementation described in ref . @xcite .    [ cols="<,>,>,>,>,>,>,>,>,>",options="header " , ]     the inspection of tab . [ tab4 ] reveals that several approaches are capable to yield very good results , with mean absolute errors close to 1 mha and mean relative errors smaller than 1% . in particular , all the calculations using raw basis - set results ( except v4zcv ) perform remarkably well , showing a significantly better performance with respect to the case of absolute energies , thanks to the cancellation of systematic errors in the present case . interestingly little or no benefit comes instead from the use of extrapolation techniques . in fact , 67-extrapolation results are in line with v7zcv results , while 56-extrapolation atomization energies are even worst that v5zcv ones and of similar quality as the v4zcv results . this result shall be rationalized in terms of an inevitable small but random computational noise that affects the extrapolation data ( especially if small basis sets are used as a base for the extrapolation ) . the small inaccuracies of extrapolated data can be in fact almost irrelevant for absolute energies ( below 1 mha ) , but can easily sum to several mha in the case of atomization energies ( especially for many atoms molecules ) . finally , we remark the good performance of the semiempirical extrapolation method , which yields atomization energies in good agreement with the reference ones , with a mean absolute relative error of only 0.75% and very small maximum absolute deviation and standard deviation .    the results of tab . [ tab4 ] show that the quality of each extrapolation procedure to yield accurate atomization energies can not be directly inferred from the results obtained for absolute energies , because atomization energies are energy differences and thus a complex error propagation can occur . as we showed this effect is more important for low - level approaches and is somehow amplified for extrapolated results , because of the presence of unsystematic errors due to the numerical procedure . it is thus finally important to note that this issue applies also for our reference atomization energies , that were obtained from the best estimated rpa absolute energies of tab . [ tab1 ] , although in this case we may expect the effect to be rather limited due to the high quality of the extrapolation . nevertheless , we have to remark that possibly errors of few mha are plausible for some systems , so that in general the set can not be fully considered as an accurate benchmark set of atomization energies , although it could be safely considered as reference for a qualitative discussion of the trends obtained from different extrapolation methods . the slow basis set convergence of the rpa correlation energy is a major problem for the development and the application of the method in quantum chemistry and solid - state physics . in fact , a careful assessment of the rpa methodology and related methods ( e.g. , rpa+ @xcite , rpa+se @xcite ) , requires the existence of a set of benchmark absolute and relative rpa - correlation energies as well as the availability of accurate extrapolation techniques . practical applications , on the other hand , would greatly benefit from efficient extrapolation schemes with well calibrated levels of confidence . surprisingly however only little effort has been dedicated in literature to study the cbs limit of rpa correlation energies and assess the effectiveness of different extrapolation procedures in this context . in this paper we aimed at pushing this work one step further and considered the problem in a more systematic way , constructing a benchmark set of reference rpa - correlation energies and studying the behavior of different extrapolation schemes in the rpa framework . our work can be summarized in the following conclusions :    * the extremely slow convergence of the rpa correlation energy with the basis set dimension prevents the computation of highly accurate reference energies . already at the v7zcv level in fact the energies can be incorrect up to 10 mha . nevertheless , we showed that a careful extrapolation procedure allows to construct an accurate benchmark set with a well defined estimation of the uncertainty . in the present work we estimate our reference data to have an accuracy close to 1 mha for most of the systems . * we analyzed several two - point extrapolation formulas against our benchmark set of absolute rpa correlation energies . it turns out that good results can only be achieved if global parameters are optimized in eqs . ( [ e5 ] ) , ( [ e6 ] ) , and ( [ e7 ] ) , as to effectively reproduce high - order terms in the theoretical asymptotic energy expansion ( @xmath2 ) . when this procedure is followed all the considered formulas provide equivalent results . moreover , with a careful optimization of the global parameters very good results can also be achieved from 56-extrapolation procedures . * rpa correlation energies evaluated at the v4zcv level appear to be of too low quality to serve as a basis in an extrapolation procedure with globally optimized parameters . this fact indicates that in the case of 45-extrapolations high - order terms beyond the @xmath3 one , play a crucial role and are thus dominant in the asymptotic expansion . therefore , rpa energies obtained at the quadruple - zeta level of theory can not be considered well converged with respect to the one - particle expansion in basis set and shall not be used in applications . nevertheless , the introduction of atomic - based parameters within a semiempirical extrapolation scheme can strongly reduce the problem and provide extrapolated rpa absolute energies of good quality ( within 2 mha from the reference ) . the semiempirical extrapolation method can thus be considered as a good tool for the cbs extrapolation in rpa applications on large systems . * when energy differences are considered , an unpredictable error propagation occurs and accurate cbs results are even more hard to achieve than the absolute correlation energies . however , because an error cancellation occurs for the systematic errors ( e.g. systematic underestimation of the energy by incomplete basis sets ) we finally found that most approaches yield results that agree within 2 mha . because no highly accurate reference data exist to assess such small differences , no clear preference can be expressed for any of the considered approaches and 2 mha ( @xmath4 kcal / mol , @xmath5 mev ) must be considered the level of confidence of state - of - the - art rpa calculations on atomization energies . therefore , caution must be always employed when small energy difference are considered . we mention finally that , following the encouraging results of the present study , further investigations will be needed to assess in detail the basis set - dependence of rpa energies and the accuracy of extrapolation schemes for energy differences ( e.g. atomization energies ) , which are a main quantity in many practical computational studies . in fact , when energy differences are considered , an unpredictable error propagation may occur and accurate cbs results are harder to achieve than the absolute correlation energies , despite the former can be directly derived as an algebraic sum of the latter . from our calculations it can be easily seen indeed that all the approaches considered in this work yield very similar atomization energies for the molecules of the test set , with differences mostly below 2 - 3 mha . however , no clear trend can be identified between the different approaches , so that no highly accurate reference data can be established to assess such small differences . this traces back to the fact that systematic errors ( e.g. systematic underestimation of the energy by incomplete basis sets ) may be expected to cancel out in this case , but extrapolation procedures may introduce in general unsystematic errors which worsen the results with respect to the raw data obtained from the corresponding basis sets . we thank turbomole gmbh for providing us with the turbomole program package , and m. margarito for technical support . this work was funded by the erc starting grant fp7 project dedom , grant agreement no . 207441 . 99 bohm d , pines d ( 1953 ) phys rev 92:609 gell - mann m , brueckner ka ( 1957 ) phys rev 106:364 langreth dc , perdew jp ( 1977 ) phys rev b 15:2884 furche f ( 2001 ) phys rev b 64:195120 fuchs m , gonze x ( 2002 ) phys rev b 65:235109 furche f , van voorhis t ( 2005 ) j chem phys 122:164106 dobson jf , wang j , dinte bp , mclennan k , le hm ( 2005 ) int j quantum chem 101:579 scuseria ge , henderson tm , sorensen dc ( 2008 ) j chem phys 129:231101 grneis a , marsman m , harl j , schimka l , kresse g ( 2009 ) j chem phys 131:154115 janesko bg , scuseria ge ( 2009 ) j chem phys 131:154106 janesko bg , henderson tm , scuseria ge ( 2009 ) j chem phys 130:081105 toulouse j , gerber ic , jansen g , savin a , ngyn jg ( 2009 ) phys rev lett 102:096404 toulouse j , zhu w , ngyn jg , savin a ( 2010 ) phys rev a 82:032502 ruzsinszky a , perdew jp , csonka gi ( 2010 ) j chem theor comput 6:127 . nguyen hv , galli g ( 2010 ) j chem phys 132:044109 paier j , janesko bg , henderson tm , scuseria ge , grneis a , kresse g ( 2010 ) j chem phys 132:094103 eshuis h , yarkony j , furche f ( 2010 ) j chem phys 132:234114 jansen g , liu rf , ngyn jg ( 2010 ) j chem phys 133:154106 ren x , tkatchenko a , rinke p , scheffler m ( 2011 ) phys rev lett 106:153003 eshuis h , furche f ( 2011 ) j phys chem lett 2:983 ruzsinszky a , perdew jp , csonka gi ( 2011 ) j chem phys 134:114110 ngyn jg , liu rf , toulouse j , jansen g ( 2011 ) j chem theory comput 7:3116 heelmann a ( 2011 ) j chem phys 134:204107 eshuis h , bates je , furche f ( 2012 ) theor chem acc 131:1084 paier j , ren x , rinke p , scuseria ge , grneis a , kresse g , scheffler m ( 2012 ) new j phys 14:043002 ren x , rinke p , blum v , wieferink j , tkatchenko a , sanfilippo a , reuter k , scheffler m ( 2012 ) new j phys 14:053020 ren x , rinke p , joas c , scheffler m ( 2012 ) j mater sci , doi : 10.1007/s10853 - 012 - 6570 - 4 heelmann a , grling a ( 2010 ) mol phys 108:359 heelmann a , grling a ( 2011 ) phys rev lett 106:093001 heelmann a , grling a ( 2011 ) mol phys 109:2473 godby rw , schlter m , sham lj ( 1986 ) phys rev lett 56:2415 kotani t ( 1998 ) j phys cond matt 10:9241 niquet yn , fuchs m , gonze x ( 2003 ) phys rev a 68:032507 grning m , marini a , rubio a ( 2006 ) j chem phys 124:154108 hellgren m , von barth u ( 2007 ) phys rev b 76:075107 hellgren m , von barth u ( 2010 ) j chem phys 132:044101 hellgren m , rohr dr , gross eku ( 2012 ) j chem phys 136:034106 verma p , bartlett rj ( 2012 ) j chem phys 136:044105    harl j , kresse g ( 2008 ) phys rev b 77:045136 eshuis h , furche f ( 2012 ) j chem phys 136:084105 janesko bg , henderson tm , scuseria ge ( 2009 ) j chem phys 131:034110 zhu w , toulouse j , savin a , ngyn jg ( 2010 ) j chem phys 132:244108    schwartz c ( 1962 ) phys rev 126:1015 schwartz c ( 1963 ) methods comput phys 2:241 carroll dp , silverstone hj , metzger rm ( 1979 ) j chem phys 71:4142 schmidt hm , von hirschhausen h ( 1983 ) phys rev a 28:3179 hill rn ( 1985 ) j chem phys 83:1173 kutzelnigg w , morgan iii jd ( 1992 ) j chem phys 96:4484 schmidt hm , linderberg j ( 1994 ) phys rev a 49:4404 mller t , basis set accuracy and calibration in quantum chemistry , in computational nanoscience : do it yourself ! , grotendorst j , blugel s , marx d ( eds . ) , john von neumann institute for computing , julich , nic series , vol . 31 , isbn 3 - 00 - 017350 - 1 , pp . 19 - 43 , 2006 . feller d , peterson ka , grant hill j ( 2011 ) j chem phys 135:044102 barkowies d ( 2007 ) j chem phys 127:084105 barkowies d ( 2007 ) j chem phys 127:164109 varandas ajc ( 2007 ) j chem phys 126:244105 perdew jp , burke k , ernzerhof m ( 1996 ) phys rev lett 77:3865 dunning jr th ( 1989 ) j chem phys 90:1007 woon de , dunning jr . th ( 1994 ) j chem phys 100:2975 feller d , peterson ka ( 1999 ) j chem phys 110:8384 feller d , sordo ja ( 2000 ) j chem phys 113:485 woon de , dunning jr . th ( 1995 ) j chem phys 103:4572 peterson ka , dunning jr . th ( 2002 ) j chem phys 117:10548 grant hill j , mazumder s , peterson ka ( 2010 ) j chem phys 132:054108 boyd drj ( 1955 ) j chem phys 23:922 venkateswarlu p , gordy w ( 1955 ) j chem phys 23:1200 redington rl , olson wb , cross pc ( 1962 ) j chem phys 36:1311 herzberg g , electronic spectra and electronic structure of polyatomic molecules , van nostrand , new york , 1966 . sverdlov lm , kovner ma , krainov ep , `` vibrational spectra of polyatomic molecules '' wiley , new york , 1974 . tsuboi m , overend j ( 1975 ) j mol spec 52:256 hirota e ( 1979 ) j mol spect 77:213 hoy ar , bunker pr ( 1979 ) j mol spect 74:1 huber kp , herzberg g , molecular spectra and molecular structure . constants of diatomic molecules , van nostrand reinhold co. , 1979 . herbst e , messer jk , delucia fc ( 1984 ) j mol spect 108:42 junttila ml , lafferty wj , burkholder jb ( 1994 ) j mol spect 164:583 gurvich lv , veyts iv , alcock cb , thermodynamic properties of individual substances , fouth edition , hemisphere pub . co. , new york , 1989 . kuchitsu k ( editor ) `` structure of free polyatomic molecules - basic data '' springer , berlin , 1998 . nist diatomic spectral database ( www.physics.nist.gov/physrefdata/molspec/diatomic/index.html ) turbomole v6.2 2010 , a development of university of karlsruhe and forschungszentrum karlsruhe gmbh , 1989 - 2007 , turbomole gmbh , since 2007 ; available from http://www.turbomole.com . kutzelnigg w ( 1985 ) theor chem acta 68:445 helgaker t , klopper w , koch h , noga j ( 1997 ) j chem phys 106:9639 halkier a , helgaker t , jrgensen p , klopper w , koch h , olsen j , wilson ak ( 1998 ) chem phys lett 286:243 martin jml ( 1996 ) chem phys lett 259:669 truhlar dg ( 1998 ) chem phys lett 294:45 fast pl , sanchez ml , truhlar dg ( 1999 ) j chem phys 111:2921 pitok m , riley ke , neogrdy p , hobza p ( 2008 ) chemphyschem 9:1636
Abstract:

Answer: we construct a reference benchmark set for atomic and molecular random - phase - approximation ( rpa ) correlation energies in a density functional theory ( dft ) framework at the complete basis set limit . this set is used to evaluate the accuracy of some popular extrapolation schemes for rpa all - electron molecular calculations . the results indicate that for absolute energies accurate results , clearly outperforming raw data , are achievable with two - point extrapolation schemes based on quintuple- and sextuple - zeta basis sets . moreover , we show that results in good agreement with the benchmark can also be also obtained by using a semiempirical extrapolation procedure based on quadruple- and quintuple - zeta basis sets . finally , we analyze the performance of different extrapolation schemes for atomization energies .


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 existence of self - organized critical dynamics in complex systems has traditionally been demonstrated through numerical simulation of certain classes of cellular automata referred to as sandpile models @xcite . non - linear , spatio - temporal dynamics is always essential for the emergence of soc behavior , but the details of this dynamics for a specific natural system is often poorly understood and/or not accessible to observation . in many cases the information available is in the form of time - series of spatially averaged data like stock - price indices , geomagnetic indices , or global temperature data . for scientists who deal with such data a natural question to ask is : are there specific signatures of soc dynamics that can be detected from such data ?    in this letter we shall report some results which provide a partial answer to such a question . some important statistical features of the toppling activity are common to most weakly driven sandpile models described in the literature , and these are used to formulate a stochastic model for the toppling activity signal . a benchmark case against which our results are tested , is a numerical study of the bak - tang - wiesenfeld ( btw ) sandpile @xcite . a crucial step in our work is a re - scaling of the dynamical variables which allows a natural passage to the thermodynamic ( continuum ) limit . we demonstrate that this leads to new results concerning soc scaling laws . we find that the probability density function ( pdf ) for the toppling activity is a stretched exponential or close to the bramwell - holdsworth - pinton distribution @xcite , depending on whether the sandpile is so slowly driven that avalanches are well separated , or it is driven so hard that several avalanches run simultaneously . the pdf for avalanche durations is unique in the thermodynamic limit , but is not a power law , unless we redefine the meaning of an avalanche to be the activity burst between successive times for which the activity rises above a positive threshold . implementing such a threshold yields an exponent for the avalanche duration pdf of 1.63 , in agreement with @xcite , but in contradiction to @xcite . it also gives power - law quiet - time statistics as in @xcite and thus refuting the claim in @xcite that soc implies power - law distributed avalanche durations , but poisson - distributed quiet times . the sandpile models considered in this short paper deal with a @xmath0-dimensional lattice of @xmath1 sites each of which are occupied by a certain integer number of quanta which we conveniently can think of as sand grains . the dynamics on the lattice is given by a toppling rule which implies that if the number of grains on a site exceeds a prescribed threshold , the grains on that site are distributed to its nearest neighbors . if the occupation number of some of these neighbors exceed the toppling threshold these sites will topple in the next time step , and the dynamics continues as an avalanche until all sites are stable . the details of this toppling rule can vary , but a useful theory for a broad class of natural phenomena should not be very sensitive to such detail .    in natural systems the soc dynamics is usually driven by some weak random external forcing . in sandpile models this can be modeled by dropping of sand grains at randomly selected sites at widely separated times . in numerical algorithms this is often done by dropping sand grains only at those times when no avalanche is running . this ensures that the drive does not interfere with the avalanching process . usually it will then only take a few time steps from one avalanche has stopped until a new starts , so for a large system the quiet times between avalanches will appear insignificant compared to their durations . a more physical drive would be to drop sand also during avalanches . if the dropping rate is slower than the typical duration of a system - size avalanche the drive would still not interfere with the avalanche dynamics , but the quiet times would depend on the statistical distribution of dropping times , which is typically a poisson distribution . in many natural systems , however , avalanching occurs all the time , corresponding to a higher driving rate . in such cases , and also because there will always be noise in time - series data , we can not identify the start and termination of an avalanche from a zero condition of our observable . in practice we have to define avalanches as _ bursts _ in the time series identified by a threshold on the signal @xcite . in a sandpile simulation such bursts are correlated and therefore the quiet times between the bursts are power - law distributed even if the dropping of sand grains is chosen to be a poisson process . hence if focus is on modeling features that can be detected in observational data we shall think of avalanches as activity bursts starting and terminating at a non - zero threshold value . moreover , one of the main results of this work is that power - law shape of the pdf for avalanche duration is true _ only _ if one defines avalanches in this way . in the btw sandpile . b ) : the increments @xmath2 of the trace in ( a ) , showing that @xmath3 is large when @xmath4 is large . c ) : conditional pdfs of @xmath5 for @xmath6 respectively . d ) the conditional mean and variance of @xmath7 versus @xmath8.,width=321 ]    we shall assume that the lattice has linear extent @xmath9 with @xmath1 sites , so the thermodynamic limit @xmath10 can be thought of as a continuum limit . the sandpile evolves in discrete time steps labeled by @xmath11 , and the number of sites whose occupation number exceeds the toppling threshold at time @xmath12 is called the toppling activity @xmath13 . the toppling increment is @xmath14 . let us define two active sites as dynamically connected if they have at least one common nearest neighbor , and define a connected cluster as a collection of active sites which are linked trough such connections . from numerical simulations of sandpiles we observe that such clusters never consist of more than a few elements and that the instantaneous number of clusters @xmath15 increases in proportion to @xmath8 . this implies that at each time @xmath12 we can label the clusters by @xmath16 , where @xmath17 is a constant depending on the specific toppling rule and the dimension @xmath18 of the sandpile . we can then decompose the increment @xmath19 into a sum of local increment contributions @xmath20 produced by each of the clusters , i.e. @xmath21 we think of the local increment contributions as random variables which take values in a finite sample space . indeed , if each cluster @xmath22 only consists of a single overcritical site , then @xmath23 takes values in the set @xmath24 .    as a first step to a stochastic model we make a mean - field assumption @xcite , which impiles that @xmath20 and @xmath25 are statistically independent for @xmath26 . then the central limit theorem states that in the limit @xmath27 , @xmath28 the conditional probability density @xmath29 $ ] of an increment @xmath19 , given @xmath13 , is gaussian with variance @xmath30 , where @xmath31 -(e[\xi_{n , i}|x_n])^2)$ ] . this has been verified numerically in the two - dimensional btw - model as shown in fig . the figure demonstrates the need to introduce a conditional probability : the conditional variance of the increments is proportional to @xmath8 and the conditional mean is not zero . in fact , numerical simulations show that the the conditional mean increment , @xmath32 $ ] , is positive for small @xmath8 , reflecting the natural tendency for the activity to grow when it is small . on the other hand the mean increment decays exponentially to zero for moderate @xmath8 , and becomes negative when @xmath8 is comparable to the activity of a system - size avalanche , reflecting the limiting influence of the finite system size . these effects will be incorporated as a drift - term correction to the model , but for now we consider for simplicity of argument a gaussian process with non - stationary increments and no drift term : @xmath33 where @xmath34 is a stationary gaussian stochastic process with unit variance . from the numerical sandpile data ( see fig . [ fig1 ] ) we observe that the _ normalized _ toppling process @xmath35 has the characteristics of a fractional brownian walk with hurst exponent @xmath36 on time scales shorter than the characteristic growth time for a system - size avalanche , consistent with a power spectrum which scales like @xmath37 . thus we model the normalized increment process as @xmath38 , where @xmath39 is a fractional brownian walk with hurst exponent @xmath40 . for the transition to the thermodynamic limit , where time will become a continuous variable , we can think about @xmath39 as the result of a discrete sampling of the ( continuous - time ) fractional brownian motion ( fbm ) @xmath41 . this process has the property @xmath42 . we now have a stochastic difference equation @xmath43 numerical simulations show that @xmath44 , is not valid for very large activity . the effect of imperfect scaling with increasing @xmath45 can be built into eq . [ sde ] through an @xmath45-dependent drift term . however , the distributions of duration and size of sub - system size avalanches ( defined by a threshold @xmath46 ) is not sensitive to this feature of the bwt model . we have given a detailed treatment of this problem in @xcite . ] where @xmath47 can be interpreted as a fractal dimension of the set of active sites imbedded in the @xmath18-dimensional lattice space . this property is used to re - scale @xmath13 such that it has a well - defined limit as @xmath48 . we also have to re - scale the time variable by letting @xmath49 , where @xmath50 . the value of @xmath51 will become apparent if we define the normalized activity variable @xmath52 , such that the corresponding increment becomes @xmath53 where @xmath54 . a well - defined thermodynamic limit @xmath27 requires @xmath55 , for which eq . ( [ stocheq1 ] ) , by introduction of the limit function @xmath56 , reduces to the stochastic differential equation @xmath57 where we have heuristically added a drift term @xmath58 to account for the non - zero mean of the conditional increment . we take @xmath59 to be an exponentially decaying function based on the numerical results from the sandpile . in the @xmath60-dimensional btw model we find that @xmath61 and hence @xmath62 . this defines re - scaled coordinates @xmath63 and @xmath64 . from simulations of the @xmath60-dimensional btw sandpile for @xmath65 . also shown is @xmath66 found from simulations of eq . ( [ sde ] ) , and a stretched exponential fit ( dashed curve , vertically shifted for visibility ) . all pdfs are scaled to unit variance.,width=321 ]    a time series @xmath67 , representing a succession of avalanches with zero quiet times , can be constructed numerically from the discrete - time version of eq . ( [ sde ] ) by integrating the equation using realizations of the fractional gaussian noise process @xmath68 . at those times when @xmath69 drops below zero we consider the avalanche as terminated , and a new , independent realization of @xmath68 is generated and used to produce the next avalanche . from long , stationary time - series generated from the stochastic model and from the sandpile model this way , we can construct pdfs @xmath66 which turn out to give almost identical results for the two models ( see fig .  [ fig2 ] ) . the shape of this pdf is universal in the thermodynamic limit : a stretched exponential @xmath70 with @xmath71 . a different pdf appears if the time - series are constructed by launching the avalanches at random times ( poisson - distributed ) with characteristic time between launches shorter than the growth time of a system size avalanche . in this case several avalanches may run simultaneously , and @xmath66 from both models are close to the bramwell - holdsworth - pinton distribution , which was claimed to be valid for the toppling - activity in the btw - model in @xcite . consider a solution of eq . ( [ sde ] ) with initial condition @xmath72 , and let @xmath73 be the evolution of the density distribution in @xmath74-space of an ensemble of realizations of the stochastic process @xmath69 all launched at activity @xmath75 at time @xmath76 . every realization @xmath69 will sooner or later terminate at a finite time @xmath77 for which @xmath78 and @xmath79 , and then we remove it from the ensemble . @xmath73 contains information about all commonly considered avalanche characteristics . for example , it is easily found from from eq . ( [ sde ] ) that , on time scales shorter than the growth time of a system - size avalanche , @xmath69 is a self - similar process with non - stationary increments and self - similarity exponent @xmath80 @xcite . hence the variance of @xmath69 with respect to @xmath73 will scale as @xmath81 . that this relation holds for the 2-dimensional btw model can easily be verified through numerical simulation ( fig . [ fig3](a ) ) .    we can also compute the survival probability @xmath82 , which is the probability that a realization of an avalanche has not terminated at the time @xmath83 . this function is related to the pdf for avalanche durations by @xmath84 , so that @xmath85 is a power law if and only if @xmath86 is a power law . [ fig3](b ) shows the function @xmath86 for numerical simulations of the btw sandpile in the re - scaled coordinates @xmath87 and @xmath88 , demonstrating that the pdf for avalanche durations does _ not _ represent a power law . the power - law form @xmath89 proposed in @xcite can only be obtained as a tangent to the log - log plot of @xmath86 at a given duration time @xmath83 , and the slope of this tangent depends crucially on the duration time @xmath83 for which this tangent is drawn . the situation changes if we let all avalanches terminate when @xmath74 drops below a small threshold @xmath46 as proposed in @xcite . in this case avalanche durations are the return times to the line @xmath90 , and by changing coordinates to @xmath91 we see that this corresponds to the return times to the time axis of the process given by the stochastic differential equation @xmath92 . for small avalanches where @xmath93 we can approximate this expression with @xmath94 , i.e. can approximate @xmath95 by a fractional brownian motion with hurst exponent @xmath40 . using the result of ding and yang @xcite on the return times of a fractional brownian motion we get @xmath96 . numerical simulations of the btw model verify this result : the survival function @xmath86 becomes a power law on time scales shorter than a system - size avalanche ( see fig . [ fig4](a ) ) , and the slope of the graph in a log - log plot is approximately @xmath97 , which corresponds to a scaling of the pdf for duration times on the form @xmath98 . the result is also reproduced by simulations of eq . ( [ sde ] ) with an exponentially decaying drift term . [ fig4](b ) shows the log - log plot of the pdf for duration times in the stochastic differential equation and a line with slope @xmath99 , demonstrating that the avalanche statistics in the btw sandpiles is captured by the stochastic differential equation .     with respect to the pdf @xmath73 . the variance grows like @xmath100 , with @xmath101 for times less than the duration of a system size avalanche . b ) double - logarithmic plots of the survival function @xmath86 in the re - scaled coordinates @xmath87 and @xmath88 , demonstrating that the pdf of avalanche durations is not a power - law . the dotted line has slope @xmath102.,width=321 ]     and @xmath88 for @xmath103 when the durations are defined by putting a small threshold @xmath104 on the toppling activity . the function shows power - law behavior with exponent @xmath97 for avalanches smaller than system size . b ) the pdf for duration times from simulations of eq . ( [ sde ] ) when avalanches are defined in the same way as for the sandpiles . the dotted line has slope @xmath99.,width=321 ]    from the scaling @xmath105 we can deduce an exponent for the pdf of avalanche size as well . on the time scales where the toppling activity can be approximated by a fractional brownian motion @xmath41 , the signal disperses with time as @xmath106 , the size of an avalanche of duration @xmath83 scales like @xmath107 . assuming that the pdf for avalanche size is on the form @xmath108 , the relation @xmath109 yields @xmath110 , so @xmath111 with @xmath112 we obtain @xmath113 . we also remark that if we omit the drift term and let @xmath114 and @xmath115 we obtain the so - called mean - field theory of sandpiles . in this case the stochastic differential equation has a corresponding fokker - planck equation @xmath116 if we solve this equation on the interval @xmath117 with absorbing boundary conditions in @xmath118 we can obtain an analytical expression for @xmath73 , and from some straightforward algebra we find for large @xmath83 that @xmath119 @xcite . since @xmath115 we can not approximate the toppling activity by a brownian motion on any scale and thus @xmath69 diverges like @xmath120 , where @xmath80 . by replacing @xmath40 with @xmath80 in eq . ( [ rel ] ) we get @xmath121 , in agreement with previous mean field approaches @xcite . we point out that the validity of eq .  [ sde ] is not restricted to the btw model . for instance , the equation has been verified for the zhang model @xcite , though with a different hurst exponent @xmath40 . time series of global quantities derived from numerical simulation of different sandpile and turbulent fluid systems can be shown to be adequately described by eq .   4 , where @xmath40 and the specific form of the drift term depend on the system at hand @xcite .
Abstract:

Answer: a stochastic theory for the toppling activity in sandpile models is developed , based on a simple mean - field assumption about the toppling process . the theory describes the process as an anti - persistent gaussian walk , where the diffusion coefficient is proportional to the activity . it is formulated as a generalization of the it stochastic differential equation with an anti - persistent fractional gaussian noise source . an essential element of the theory is re - scaling to obtain a proper thermodynamic limit , and it captures all temporal features of the toppling process obtained by numerical simulation of the bak - tang - wiesenfeld sandpile in this limit .