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: it is generally assumed that diffusion in macroscopic systems is a result of chaos on a microscopic scale , following einstein s 1905 explanation for brownian motion , that is , the movement of a colloidal particle due to thermal fluctuations in the surrounding fluid . we describe chaos , diffusion and related properties in detail in sec . [ s : fund ] . we define microscopic chaos ( usually shortened to just `` chaos '' ) in terms of unpredictability as quantified by positive lyapunov exponents or a positive kolmogorov - sinai entropy ; later we will find that properties based on periodic orbits also provide a useful description of chaotic aspects of the microscopic dynamics . note that we restrict the use of the term `` chaotic '' to its usual meaning , that is a positive kolmogorov - sinai entropy ; while the properties of periodic orbits are important , we do not classify systems as chaotic or nonchaotic based on these properties . we call a system diffusive when the mean square displacement of a particle is proportional to the time , or a distribution of particles satisfies the diffusion equation .    a recent experiment of gaspard et al  @xcite , described below , purports to show that the diffusion of a brownian particle is due to microscopic chaos . while we believe that brownian motion ( including both the brownian particle and the solvent ) is most likely chaotic , our simulations using nonchaotic models , in a brief comment  @xcite and in greater detail here lead to the same results as found in the experiment , so that no experimental proof of microscopic chaos has been given in  @xcite . here we explore the question of what kind of experimental measurements or data analysis might be required to identify microscopic chaotic dynamics , and the connection between microscopic chaos and diffusion . we consider generalizations of models of diffusion due to lorentz and ehrenfest ( figs . [ f : fix]-[f : lor ] ) , where a single point particle undergoes elastic collisions with a fixed arrangement of circular or square scatterers in two dimensions , respectively . collisions with the circular scatterers of the ( chaotic ) lorentz gas lead to exponential separation of nearby trajectories , that is , a positive lyapunov exponent , while collisions with the square scatterers of the ( nonchaotic ) ehrenfest model lead to at most linear separation of nearby trajectories , and the lyapunov exponents are all zero . actually , we consider models with three types of scatterers , all of which are non - overlapping : circles as in the lorentz gas , squares oriented such that their diagonals are parallel to the coordinate axes and with only four particle velocity directions as in the ehrenfest wind - tree model , and squares of arbitrary orientations with arbitrary particle velocity directions as a generalization of the ehrenfest model . ( 280,280)(-50,0 ) ( 40,10 )    ( 280,280)(-50,0 ) ( 40,10 )    ( 280,280)(-50,0 ) ( 40,10 )    for each of these three models , we consider the following three cases : an infinite number of randomly placed fixed scatterers as in the original models of lorentz and ehrenfest , a small number of fixed scatterers in an elementary cell subject to periodic boundary conditions , and an arrangement of a finite number of fixed scatterers enclosed by absorbing boundaries . note that a model with periodic boundary conditions can be viewed as an infinite periodic system . for the purposes of studying properties such as the mean free time , it is better to view the model as a finite system with periodic boundaries , but to compute the diffusion coefficient in terms of the growth of the mean square displacement with time , it is necessary to view the model as an infinite periodic system . a link between chaos and diffusion involves a fundamental question of statistical mechanics , since one has to make a connection between the microscopic and macroscopic behavior of large systems . the models we consider here are very special from a physical point of view , but very attractive from a mathematical point of view , because there is only one moving particle , rather than a large number as in , for example , brownian motion . these can be considered bona fide statistical mechanical models if the large number of scatterers are included , as long as their lack of motion is irrelevant to the questions we ask . we believe that our models incorporate the essential features needed for a discussion of microscopic chaos and diffusion . we now describe the experiment and analysis of gaspard et al . in more detail . the position of a brownian particle in a fluid was measured at regular intervals of @xmath0 . the experimental time series data was then interpreted using standard techniques of chaotic time series analysis , suggesting a positive lower bound on the kolmogorov - sinai entropy , hence microscopic chaos . the method they used , due to cohen and procaccia  @xcite was adapted from an approach pioneered by grassberger and procaccia  @xcite . in this method one analyzes the distribution of recurrences , _ i.e. _ instances where the system approximately retraces part of its previous trajectory in phase space for a certain length of time , leading to a determination of the kolmogorov - sinai entropy . there is detailed mathematical justification for the method ( under assumptions such as the length of the data set being sufficient to approximate infinite time limits , see sec . [ s : gp ] ) , but the idea is very simple : recurrences give useful information about the predictability of the system . if long recurrences occur very often , it is easy to predict the future of the system from previous instances similar to the recent part of the trajectory , so the system has a small or zero kolmogorov - sinai entropy , whereas a rapid decay of the frequency with length of the recurrences indicates a high degree of microscopic chaos . thus , the cohen - procaccia method used by gaspard can be used to calculate the kolmogorov - sinai entropy from a time series in principle , but certain mathematical conditions and limits apply , sometimes restricting its applicability in practice . gaspard et al . concluded from their analysis that the kolmogorov - sinai entropy of the system containing the brownian particle was positive , and hence that the microscopic dynamics was chaotic . subsequently , we showed  @xcite that the same approach applied to an `` identical '' time series generated by a numerical simulation of the nonchaotic ( infinite , fixed orientation ) ehrenfest model yielded virtually identical results . this is surprising , since the collisions with the flat sides of square scatterers do not lead to positive lyapunov exponents , as described above . in ref . @xcite we attributed the discrepancy to the physical issue of time scales  the time interval between measurements ( @xmath0 ) was vastly greater than the typical collision times of the brownian particle with the solvent particles in the fluid ( approximately @xmath1 ) . while this is certainly an experimental problem , it leaves open the question of whether in principle a similar analysis with a much higher resolution could ever prove from experimental data that the microscopic dynamics is chaotic or not . one aim of the current work is to shed light on this , and related questions . we perform the same cohen - procaccia analysis on all the models mentioned above , finding that the same results are obtained even for a model with a periodic array of squares . this rules out the possibility that the apparently positive value of the kolmogorov - sinai entropy is due to a loss of information associated with the randomness of the many scatterers . we discuss the relevant time scales in our models , finding a single microscopic time scale , and we can then qualitatively explain the results of the cohen - procaccia algorithm when the time step is shorter or longer than this time . we conclude that due to the rarity of recurrences in a diffusive system , the determination of chaoticity in such a system requires vastly longer time series than are practical for experimental or computational work , even when the measurements are taken at microscopic time and distance scales .    at that point we return to the question of how microscopic dynamics manifests itself in a time series , as well as in the macroscopic behavior . we develop a new method which we call `` almost periodic recurrences '' that selects a particular class of recurrences different from those used in the cohen - procaccia method . our approach is designed to distinguish systems based on their periodic orbit properties , which are related to , but not equivalent to the usual definition of chaos as positive kolmogorov - sinai entropy . we find that , in contrast to the cohen - procaccia method , this method can distinguish between the chaotic lorentz model and the nonchaotic ehrenfest model and thus reveals at least one way in which microscopic chaos is manifest , although on macroscopic scales both models exhibit diffusion . the basic idea of our method is simple : almost periodic recurrences are related to periodic orbits which have very different properties in chaotic and nonchaotic systems , so a search for periodic orbits using recurrences may distinguish between the two classes of systems . the periodic orbits of chaotic systems ( such as the lorentz gas ) are exponentially unstable , so many repetitions close to a periodic orbit are unlikely . in contrast , the periodic orbits of nonchaotic systems ( for example our wind - tree models ) can be power law unstable , allowing many repetitions . the markedly different properties of periodic orbits in the lorentz and ehrenfest models leads to the following important observations about diffusion in finite geometries . solutions of the diffusion equation with absorbing boundaries exhibit exponential decay , corresponding to the probability of a particle remaining in the system for a given time . the randomly oriented wind - tree model has periodic orbits with power law escape , so that the escape from the whole system can not be exponential at long times . the fixed orientation wind - tree model has no periodic orbits at all , so the particle can not remain in the system longer than a fixed maximum time . both of these are examples of situations in which the infinite model satisfies the diffusion equation by having a gaussian distribution with the mean square displacement proportional to the time , but the corresponding finite model does not satisfy the diffusion equation . the outline of this paper is as follows : first we introduce our models ( sec .  [ s : mod ] ) , and discuss in detail their relevant time scales and chaotic and diffusive properties ( sec . [ s : fund ] ) . then we discuss the cohen - procaccia method and its inability to distinguish chaotic from nonchaotic dynamics in our case ( sec . [ s : gp ] ) . finally we introduce our alternative `` almost periodic recurrence '' method and discuss the properties of the periodic orbits ( sec . [ s : isr ] ) , and their consequences for finite systems ( sec . [ s : esc ] ) . we conclude with a discussion of our results and some open questions ( sec.[s : conc ] ) . the ehrenfest wind - tree model describes a point `` wind '' particle moving in straight lines in the plane punctuated by elastic collisions with fixed square scatterers , the `` trees '' . in the original model ( fig . [ f : fix ] ) , the trees are oriented with their diagonals along the @xmath2 and @xmath3 axes , and the wind particle moves with a fixed velocity in only four possible directions , along these axes . the trees are located at random positions , with a given number density , and not overlapping . there is an overlapping version we will not consider ; it exhibits anomalous ( sub)-diffusion and non - gaussian distributions  @xcite . secondly , we consider the case of randomly oriented scatterers ( fig . [ f : ran ] ) , which allows all possible wind particle directions . in both the fixed and randomly oriented cases all lyapunov exponents are zero , because a collision from a flat scatterer causes only a linear separation of nearby trajectories ; it does not cause the exponential separation associated with convex curved scatterers . thirdly , we consider the ( two dimensional ) lorentz gas ( fig . [ f : lor ] ) , that is , a model where the scatterers are circular , and we choose the area and number density of the scatterers to correspond to the wind - tree models . this model is known to have normal diffusive properties  @xcite and is used here as a comparison to examine the effects of positive versus zero lyapunov exponents .      for the numerical simulations , we take the velocity of the wind particle to be @xmath4 . for the wind - tree models we take the side length @xmath5 equal to @xmath6 . this is done for computational convenience ; space is divided into unit squares ( of side length 1 ) aligned with the coordinate axes , so that at most one scatterer can be contained in each unit square . for the lorentz gas the circular scatterers have a radius @xmath7 to give them the same area @xmath8 as the squares . in all cases the total area considered has a number of unit squares , @xmath9 , up to @xmath10 in both @xmath2 and @xmath3 directions , and periodic boundary conditions are used . to simulate an infinite system , a large value of @xmath9 ( up to 3500 ) is used ; for the maximum time of @xmath11 time units the particle never travels far enough to sample the periodic boundary conditions . to simulate a periodic system , a small size such as @xmath12 is used ; we always consider the case of finite horizon , so that the time between collisions is bounded . to simulate a system with absorbing boundaries , an intermediate size such as @xmath13 is used ; the absorbing boundaries are one distance unit inside the edge of the system to avoid the possibility of the particle colliding with a scatterer as well as with its periodic image . we use a number of scatterers @xmath14 equal to @xmath15 . this means that the number density @xmath16 is @xmath17 in all cases and the proportion of the total area covered by the scatterers is @xmath18 . thus we have an intermediate density of scatterers , which is most convenient from the point of view of the simulations and also the time scales . for a significantly lower density , an infinite horizon is much more likely in all but the largest systems . for a significantly higher density ( limited by a maximum of @xmath19 ) the metropolis algorithm ( see below ) used to position the scatterers is very much slower . also , both for low and high density systems , there is an additional time scale ( see section  [ s : time ] ) , which complicates the analysis . we now describe the method used to position the fixed scatterers . each configuration of random positions and orientations is obtained from a version of the metropolis algorithm , used in  @xcite . the square scatterers of the wind - tree models are placed initially at points belonging to a square lattice defined such that they do not overlap , that is , with their centers at integer coordinates @xmath20 such that @xmath21 is even , and their orientation as in fig . [ f : fix ] . not every lattice site is then occupied , depending on the number density . each scatterer is shifted and rotated a small random amount in turn , typically a few tenths of a distance unit and of order ten degrees respectively . if the shift causes an overlap , the move is rejected and the previous configuration is used . the procedure is repeated with different shifts and rotations . a large fixed number of shifts are attempted , sufficient for reasonable measures of the correlation between scatterers to have long converged to their final values .    if the circles of the lorentz gas are placed on the same lattice as the squares they overlap slightly ; at the chosen density of @xmath22 there is enough freedom to allow them to find non - overlapping locations , so that the final configuration is non - overlapping , as we could check . at higher densities circles may not find non - overlapping locations , and a close packed triangular lattice should then be used instead . we now define the notation we will use to distinguish between the various models . there are a very large number of possible parameters that could appear in general , referring to the shape of the scatterers ( squares , circles , etc . ) , different densities and different shaped boundary conditions to name a few , but for conciseness we limit the notation to include only those models we use here , in particular keeping the areal density always at @xmath22 . using square brackets [ , ] to denote different cases , our notation is [ f , r , l][@xmath23,[p , a]@xmath9 ] . the first symbol [ f , r , l ] denotes the type and orientation of the scatterers , f for fixed oriented squares , r for randomly oriented squares and l for circles ( the lorentz gas ) . the next symbol [ @xmath23,p , a ] gives information about the boundary conditions , either @xmath23 , that is , no boundaries , or p for periodic boundary conditions or a for absorbing boundaries . the symbols p and a are followed by the size of the system , @xmath9 , which is an even integer . thus the original ehrenfest model is denoted f@xmath23 , while a periodic lorentz gas with @xmath25 ( containing @xmath26 scatterers at the density used in our simulations ) is denoted lp8 ( see fig . [ f : lor ] ) . the latter model is generated by random shifts of the scatterers as above , ensuring that periodic images of the scatterers do not overlap . when no confusion can arise , we will sometimes refer to classes of models with a simplified notation , for example we denote all randomly oriented wind - tree models as r , and all lorentz models with absorbing boundaries as la . this section collects a number of fundamental properties and results that form a basis for sections  [ s : gp ] ,  [ s : isr ] and  [ s : esc ] . first we consider the important issue of _ time scales _ , which puts the time series analysis methods of sections  [ s : gp ] and  [ s : isr ] into a physical context . then we consider _ chaotic properties _ , giving known results and conjectures for our models ; all of the following sections use this . after this we study _ diffusive properties _ , obtaining some new results . these are of direct importance to sec . [ s : esc ] ; the time series analysis methods are applied only to our diffusive models , . ] but the diffusion is somewhat incidental as the methods can be applied to arbitrary time series . at the end of this section we relate the chaotic and diffusive properties of our models . one of our difficulties of the gaspard et al experiment  @xcite is that the interval between their measurements ( 1/60 s ) is so much greater than the relevant time scale ( @xmath1 ) of the dynamics , that is the time interval determined by collisions between the brownian particle and the other particles in the fluid . it is thus important to clarify the issue of what time scales are relevant in our models , so that the simulation results can be put in context . there are three relevant time scales here , the time taken for the particle to traverse the length of a scatterer ( defined to be of order one time unit ; see sec . [ s : comp ] ) , the mean free time between collisions @xmath27 , and the time at which the particle begins to notice the finite size @xmath9 of the system ( when it is not infinite ) . it turns out ( see sec .  [ s : diff ] ) that all the random and periodic models except fp have a well defined diffusion coefficient @xmath28 , so in these cases the time taken to reach the boundary is of order @xmath29 . the mean free time can be calculated exactly ( see the appendix ) and is given by @xmath30 , @xmath31 and @xmath32 for the f , r , and l models respectively , independent of the size of the system ( ignoring the case of absorbing boundaries ) . we have checked these values numerically , and the results agree within their uncertainties of about @xmath33 . it is clear that all the time scales in these models are of order one time unit , except those defined by the boundary . this means that a time step of one time unit should be sufficient to observe effects due to the microscopic chaos in the time series analysis methods of sections  [ s : gp ] and  [ s : isr ] . we began with the notion of microscopic chaos as the presence of a positive lyapunov exponent arising from collisions with strictly convex scatterers , or as a positive kolmogorov - sinai entropy quantifying the lack of predictability of a chaotic system . we now want to make these ideas more precise and clarify what is known and what is conjectured about our models .    from a mathematical point of view , there are a large number of dynamical properties that a system may exhibit  @xcite . those of interest to us are    1 . we need ergodicity for the time series analysis methods of secs .  [ s : gp ] and  [ s : isr ] ; that is , a single long trajectory is supposed to sample the dynamics of the whole phase space . we either know or conjecture that ergodicity holds in all of our models ( see below ) . 2 .   decay of velocity correlations . we note that eqs . ( [ e : gk],[e : gkb ] ) below give the diffusion and burnett coefficients as integrals over velocity autocorrelation functions : the diffusion coefficient involves two - time correlations , while the burnett coefficient involves four - time correlations . for these coefficients to be defined the appropriate correlations must decay sufficiently quickly . numerical evidence for the existence or nonexistence of these coefficients for our models is given in sec . [ s : diff ] . 3 .   a positive lyapunov exponent , corresponding to exponential separation of initially close trajectories , follows from the shape of the scatterers . the lorentz gas has a positive lyapunov exponent due to its convex scatterers , while the wind - tree models have all zero lyapunov exponents due to their ( piecewise ) flat scatterers . 4 .   a positive kolmogorov - sinai entropy , which is our definition of chaos , and which the cohen - procaccia method ( see sec . [ s : gp ] ) is designed to compute . we either know or conjecture that the ks entropy is equal to the sum of the positive lyapunov exponents in our models ( pesin s theorem , see below ) . actually there is at most one positive lyapunov exponent in these systems . periodic orbit properties sensitive to chaoticity are described in sec . [ s : isr ] and used in secs . [ s : isr ] and  [ s : esc ] . we now briefly discuss what is known about our models regarding these dynamical properties . the periodic lorentz gas has exponential decay of ( all two time ) correlations  @xcite , which implies ergodicity and the convergence of the integral in  ( [ e : gk ] ) below . for  ( [ e : gkb ] ) see ref . the periodic wind - tree models are generically ergodic  @xcite . pesin s theorem holds for all of these [ f , r , l]p models , see refs . @xcite for the lorentz gas and ref . @xcite for the wind - tree models . the infinite models are much more difficult to treat rigorously because the phase space is not compact . in order to make reasonable conjectures about the above properties , we use the solution of the diffusion equation  ( [ e : gauss ] ) below , which is verified numerically in sec . [ s : diff ] for the models [ f , r , l]@xmath23 , to argue that the number of different scatterers hit by the particle in time @xmath34 is proportional to @xmath35 , as follows . in two dimensions , the probability density of the particle being at any point decays as @xmath36 . this means that in a time @xmath34 , the particle comes close ( say , within a radius @xmath37 ) to the point a number of times proportional to the integral of @xmath36 , that is @xmath38 . if the point corresponds to the location of a scatterer , we could say that the particle collides with each scatterer a number of times proportional to @xmath38 . in this case the total number of different scatterers hit by the particle in time @xmath34 might be expected to grow proportional to @xmath35 , which is what we find numerically in fig . [ f : stat ] . ( 300,320)(-70,0 ) ( 0,20 ) ( 155,10)@xmath34 ( -15,165)@xmath39 ( 40,285)l@xmath23 ( 40,275)r@xmath23 ( 40,265)f@xmath23    these observations have the following consequences for the chaotic properties of the infinite models . because the particle returns to each scatterer an infinite number of times , we expect that it passes arbitrarily close to every point in phase space , that is , the system is ergodic . the ks entropy gives the amount of information per unit time required to predict the trajectory knowing its infinite past . the unpredictability has two sources , the instability associated with the collisions in the lorentz gas ( that is , the positive lyapunov exponent ) , and the random positions of the scatterers with which the particle has not yet collided in all the models . the number of different scatterers hit by the particle grows as @xmath35 , so the rate decreases to zero as @xmath40 . we conclude from this argument that the random positions of the scatterers do not contribute to the ks entropy , although results at finite time might suggest otherwise , given that @xmath40 decays so slowly . in other words , pesin s theorem is satisfied for these infinite models also . ) below ) in two dimensions . in three dimensions , the density decays as @xmath41 which is integrable at infinity , so that the particle collides with each scatterer only a finite number of times and thus does not sample the whole phase space . in addition , there is a contribution to the ks entropy from the random positions of the scatterers , thus violating pesin s theorem . there is no contradiction here ; the proofs of pesin s theorem all demand that the phase space be compact . ] we have discussed a number of dynamical properties with regard to our models . the chaotic lorentz gas satisfies more of these properties than the nonchaotic wind - tree models . chaos is defined here by positive ks entropy and is clearly relevant fro diffusion , yet we show in subsequent sections that the periodic orbit properties also appear to be relevant to diffusion . the conclusion discusses an example of a model which is chaotic as it has positive ks entropy , but has periodic orbit properties and diffusive properties more similar to the nonchaotic wind - tree models . although we have a substantial number of relevant numerical results below , a full understanding of the dynamical properties ( ks entropy , periodic orbit properties or perhaps others ) most closely related to diffusion still eludes us . now we discuss diffusive properties . we begin with the diffusion equation , and then relate it to the mean square displacement of a particle . the diffusion equation @xmath42 for a probability density @xmath43 with a constant diffusion coefficient @xmath28 is linear . thus its general solution is an integral over green s functions given by the solution for an initial dirac delta distribution @xmath44 , that is , @xmath45\ ] ] is the conditional probability density that a particle initially at the point @xmath46 will be at the point @xmath47 a time @xmath34 later ; this is therefore a 2-time probability distribution ( or density ) function . the spatial dimension @xmath48 in our case . the diffusion equation is known to hold for the probability density of a particle undergoing deterministic dynamics in a number of systems including the lorentz gas  @xcite . in the diffusion equation a macroscopic approximation is used , implying that @xmath43 is an average over space and time scales large compared to the characteristic microscopic lengths and times of the dynamics . the mean square displacement of the particle in the @xmath2-direction after a macroscopic time @xmath34 is obtained from ( [ e : gauss ] ) , @xmath49 where @xmath50 . this is the einstein relation for diffusion . note that @xmath51 is now an explicit function of time . similarly , the fourth and sixth cumulants are @xmath52 respectively . it is also possible to obtain dimensionless forms of the above expressions dividing them by the appropriate power of @xmath53 . in this form the fourth cumulant is usually called the kurtosis , and is a common measure of how close a probability density is to a gaussian distribution . straightforward manipulations applied to the mean square displacement lead to the green - kubo expression for the diffusion coefficient ( if it exists ) @xmath54 where @xmath55 is one ( say the @xmath2- ) component of the velocity and the subscript denotes the time . this shows that the existence of a diffusion coefficient is also related to the rate of decay of the velocity autocorrelation function . the integral diverges if the decay rate of the integrand is @xmath56 or slower , unless the integrand is oscillatory . an oscillatory integrand can not be ruled out , but is not typical for the lorentz gas as observed in numerical simulations . the diffusion equation is the lowest order macroscopic description of the diffusion process . the next approximation , in the direction of microscopic space and time scales , _ ie . _ when the probability density @xmath57 varies sufficiently rapidly , the right hand side of eq . ( [ e : diff ] ) may be augmented by terms such as @xmath58 , where @xmath59 is called the linear super - burnett coefficient  @xcite . like the diffusion coefficient , it obeys an einstein relation @xmath60 this is consistent with the zero kurtosis found previously  ( [ e : kurt ] ) : to find the kurtosis we divided by @xmath61 which is proportional to @xmath62 if the diffusion coefficient exists , while the burnett coefficient in  ( [ e : burn ] ) gives a smaller term ( for large @xmath34 ) , proportional to @xmath34 . the burnett coefficient can also be written in the form of a green - kubo relation in terms of integrals over four - time velocity correlation functions , @xmath63 where @xmath64 as above and @xmath65 etc . a slow decay of the relevant correlation functions leads more easily to a divergence in this case than in eq . ( [ e : gk ] )  @xcite . this means that it is possible to have a well defined diffusion coefficient , but a divergent burnett coefficient , which occurs if the fourth cumulant increases faster than @xmath34 but slower than @xmath62 . we estimate the cumulants defined in eqs . ( [ e : msd]-[e : six ] ) numerically by choosing a large number ( typically @xmath66 ) of random initial conditions for the particle not inside a scatterer , for a single configuration of scatterers . we report only the cumulants in the @xmath2-direction , although we have checked that the @xmath3 distribution is similar . see figs . [ f : msd]-[f : six ] and table  [ t : diff ] . ( 500,530)(-100,0 ) ( -100,30 ) ( 170,65 ) ( 155,20)@xmath34 ( -115,265)@xmath53 ( 0,442)@xmath34 ( 0,426)@xmath67 ( 0,410)@xmath62 ( 0,394)l@xmath23 ( 0,378)lp4 ( 0,362)r@xmath23 ( 0,346)rp4 ( 0,330)f@xmath23 ( 0,314)fp4    ( 300,320)(-70,0 ) ( 0,20 ) ( 165,10)@xmath34 ( -25,165)@xmath68 ( 60,261)@xmath34 ( 60,251)@xmath62 ( 60,241)@xmath69 ( 60,232)l@xmath23 ( 60,222)lp4 ( 60,213)r@xmath23 ( 60,203)rp4 ( 60,194)f@xmath23 ( 60,184)fp4    ( 300,320)(-70,0 ) ( 0,20 ) ( 165,10)@xmath34 ( -30,165)@xmath70 ( 60,261)@xmath34 ( 60,251)@xmath62 ( 60,241)@xmath71 ( 60,231)@xmath72 ( 60,222)l@xmath23 ( 60,212)lp4 ( 60,203)r@xmath23 ( 60,193)rp4 ( 60,184)f@xmath23 ( 60,174)fp4    .diffusive properties of our models , in terms of the diffusion coefficient @xmath28 , and the burnett coefficient @xmath59 if they exist . the existence of these coefficients is related to the rate of decay of the velocity autocorrelation functions of  ( [ e : gk],[e : gkb ] ) . note that by this criterion , the periodic lorentz gas lp4 is more `` chaotic '' than the infinite lorentz model l@xmath23 , while the periodic fixed orientation wind - tree model fp4 is less `` chaotic '' than the infinite model f@xmath23 . only typical values of the coefficients for the periodic models are given , since they depend on the position of the scatterers . [ t : diff ] [ cols="^,^,^,^",options="header " , ]     for five models , [ f , r , l]@xmath23 , [ r , l]p4 , the mean square displacement is found to be proportional to the time , and the fourth and sixth cumulants increase slower than @xmath62 and @xmath71 respectively at large times , indicating that the distribution approaches a gaussian . the exception to this rule is fp4 , for which the mean square displacement grows faster than the time , approximately as @xmath67 and the kurtosis is also nonzero . this superdiffusive behavior is not particularly stable in the sense that the mean square displacement for this model has larger fluctuations than the other models , and the exponent of approximately @xmath73 varies unpredictably between 1 and 2 with the exact positions of the scatterers and the size of the cell ( for example for fp6 or fp8 ) . we now look at how the chaotic and diffusive properties of our models are related in the light of our results , specifically of : ( a ) normal diffusion without chaos , ( b ) superdiffusion without chaos , ( c ) a possible relation between the diffusion coefficient and the periodic orbits , and ( d ) the burnett coefficient of the periodic lorentz gas . in no way does the data presented in figs . [ f : msd]-[f : six ] distinguish between chaotic and nonchaotic models . [ [ normal - diffusion ] ] normal diffusion + + + + + + + + + + + + + + + +    our results show that it is possible to have a well defined diffusion coefficient and a gaussian distribution function without positive lyapunov exponents , that is , in the f@xmath23 , r@xmath23 and rp4 models . the rp4 model is particularly interesting in this regard , since the two most obvious sources of unpredictability : dispersive collisions and collisions with new randomly placed scatterers , are absent . recall that we argued ( in section  [ s : chaos ] ) that these two processes might together determine the ks entropy , that is , the rate of loss of information about the particle s position and velocity . the rp4 model is thus the clearest example of a model in which there is normal diffusion , and zero ks entropy . an example with the same properties , using a periodic rhombus model has been studied in ref . @xcite , where a mean square displacement proportional to the time was observed , but the higher moments were not studied .    [ [ superdiffusion ] ] superdiffusion + + + + + + + + + + + + + +    the fp4 model , in contrast , has a mean square displacement growing faster than linearly with the time , and a non - gaussian distribution . this can be related to a slow decay of the velocity autocorrelation function ( see above ) , but a physical explanation of why this occurs in the fp4 model but not in the rp4 model is lacking .    [ [ diffusion - and - periodic - orbits ] ] diffusion and periodic orbits + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    it is perhaps somewhat hazardous to deduce more from this data than the existence of normal diffusion ; however we would nevertheless like to point out a possible connection between the value of the diffusion coefficient ( given by the vertical displacement of the curves in the inset of fig . [ f : msd ] ) in the infinite models and the properties of periodic orbits discussed in section  [ s : isr ] . we make no statement about the diffusion coefficient of the periodic models because this coefficient depends on the position of the scatterers . the idea is that the presence of periodic orbits in a model with randomly positioned scatterers might be expected to lower the diffusion coefficient , since the periodic orbits may keep the particle trapped in roughly the same place ; this takes time without contributing to the mean square displacement . thus we observe that the fixed orientation wind - tree model f@xmath23 has no periodic orbits and has the largest diffusion coefficient ; the random lorentz gas l@xmath23 has periodic orbits , but they are exponentially unstable , and has the next largest diffusion coefficient ; the random orientation wind - tree model r@xmath23 has periodic orbits with power law instability that can inhibit diffusion for a long time , and the diffusion coefficient is the smallest of the three .    [ [ the - burnett - coefficient ] ] the burnett coefficient + + + + + + + + + + + + + + + + + + + + + + +    the fourth cumulant grows linearly with time for the lp4 model , indicating a finite burnett coefficient , as expected , given its exponential decay of correlations  @xcite . the other models , f@xmath23 , r@xmath23 , rp4 and l@xmath23 have a divergent burnett coefficient , which is not surprising since their decay of correlations is most likely as a power law . thus , from the more subtle properties of the two time ( displacement ) correlation function , we can deduce some properties relating to the rate of decay of higher order correlations , even special four - time ( velocity ) correlations , but not whether there is a positive lyapunov exponent , since l@xmath23 behaves similarly to r@xmath23 from this point of view . even in the intermediate region between the ballistic behavior at small times and the ( super)-diffusive behavior at large times in fig . [ f : msd ] , that is , where the distribution functions are measured at microscopic time scales , there is no clear distinction between the chaotic and nonchaotic models . we have now discussed at some length the connection between chaotic and diffusive properties of our models . while this has yielded valuable insight into how little chaoticity is needed for diffusion to occur , it has not provided us with an understanding of how chaotic and nonchaotic diffusion differ , let alone a method to distinguish chaotic and nonchaotic diffusive systems experimentally . the reason , given in ref . @xcite is that two time correlations ( for example our cumulants ) are insufficient to characterize chaos ; multi - time correlations ( or their underlying probability distributions ) are required . the remaining sections of this paper all use multi - time distributions of one form or another to investigate our models . we begin with a more detailed study of the grassberger - procaccia method than in ref . @xcite , where we found it to incorrectly classify the f@xmath23 model as chaotic . grassberger and procaccia  @xcite gave the earliest practical methods for computing chaotic properties such as entropies and dimensions from experimental or computer generated time series . these are still the most popular methods in use , sometimes with minor variations . in ref . @xcite the conclusion of microscopic chaos using data from a brownian motion experiment was based on a slightly different method of cohen and procaccia  @xcite . both methods are clearly reviewed in ref . @xcite . we consider these methods as applied to the calculation of the ks entropy @xmath74 , the non - zero value of which characterizes chaos . we recall the discussion of section  [ s : chaos ] where we concluded that the wind - tree models all have zero ks entropy , but the lorentz gas has a ks entropy equal to its positive lyapunov exponent .    the original grassberger - procaccia method  @xcite computes a slightly different dynamical entropy @xmath75 ( defined using the square of the probability , @xmath76 rather than the conventional @xmath77 ) satisfying the inequality @xmath78 . thus a positive estimate @xmath79 implies that @xmath80 and hence chaoticity . we follow ref . @xcite in using an adaptation due to cohen and procaccia  @xcite that allows @xmath74 to be estimated directly . we replace the experimental time series of the positions of the brownian particle in ref . @xcite by a numerical time series of any of our infinite or periodic diffusive models , containing the @xmath2 and @xmath3 positions of the particle at @xmath11 times uniformly spaced by a separation @xmath81 . unlike ref . @xcite , by removing data points . the new time step was denoted @xmath82 . we do not need to do this in our numerical work , and we do not use the notation @xmath82 for time steps to avoid confusion with the mean free time . ] our values of @xmath81 , equal to 1 and 0.01 time units are close to the dynamical time scale determined by the mean free time @xmath83 ( see sec .  [ s : time ] ) .    a single point on our numerical trajectory is denoted @xmath84 , where @xmath85 is an integer ranging from 0 to @xmath86 . for the fixed oriented square ( f ) models the trajectory contains all the phase space information except the velocity which can take only four values in this case . for the randomly oriented square ( r ) and the lorentz ( l ) models the velocity has one continuous degree of freedom . the brownian motion experiment uses only one component ( say , @xmath2 ) of the particle position , ignoring the huge number of other degrees of freedom contained in the fluid system . from a mathematical point of view , as long as the degrees of freedom are all coupled ( that is , all the degrees of freedom depend on each other directly or indirectly ) , this type of method usually converges to the correct entropy ( or dimension ) using as few as one measured variable , however the efficiency may then well be such that an unreasonably large number of data points is required for a reliable estimate of these quantities . a number ( typically @xmath87 to @xmath88 ) of ( in our case non - overlapping ) uniformly spaced subsequences of the long trajectory are denoted `` reference segments '' ( `` reference trajectories '' in ref . @xcite ) , having lengths @xmath89 ranging from 1 to 100 . these reference segments are compared with all subsequences of the trajectory having the same length . the comparison is made according to the position relative to the initial point of each subsequence , with a euclidean metric in real space and a maximum metric over the trajectory . thus , if @xmath85 and @xmath90 give the positions in the full trajectory of the initial points of the @xmath91 reference segments and @xmath92 total segments respectively , the distance between the segments beginning at @xmath85 and @xmath90 is defined as @xmath93 where @xmath94 etc . we are now in a position to count the number of recurrences , in what is called `` pattern probability '' in  @xcite . for a given tolerance or spatial resolution @xmath37 we can define the probability of a certain pattern defined by the reference sequence beginning at @xmath85 as the proportion of all the subsequences of the full trajectory that are within a distance ( as defined above ) @xmath37 of the reference sequence : @xmath95\ ] ] and from this the `` pattern entropy '' @xmath96 which gives the information theoretic entropy of the dynamics with respect to the pattern ( of length @xmath89 and resolution @xmath37 ) probability distribution . the pattern entropy thus contains information about the @xmath89-time ( displacement ) distribution functions , as suggested by the discussion at the end of the previous section . when the pattern entropy increases linearly with time in the limit of large @xmath89 , it is possible to define an entropy per unit time , @xmath97\ ] ] in practice @xmath89 is finite , and @xmath74 is always bounded above by @xmath98 ( since it is an average of @xmath99 , where @xmath57 bounded below by @xmath100 ) and tends to saturate towards this value at long times . however , we can often find a large enough linear region in a plot of @xmath74 as a function of @xmath89 ( for fixed @xmath37 and @xmath81 ) to compute @xmath101 . in ref . @xcite @xmath101 was plotted as a function of @xmath37 with different curves corresponding to different @xmath81 . from dimensional arguments it follows that any diffusive process has a scale invariance leading to a dependence @xmath102 , which is in fact observed in such plots .    assuming that there is sufficient data that all the limits ( of large @xmath14 , @xmath91 and @xmath89 ) can be approximated reliably , the ks entropy @xmath74 is equal to the maximum value of @xmath101 as @xmath37 and @xmath81 are varied . the measured value @xmath101 will be less than the true ks entropy @xmath74 if the values of @xmath37 and or @xmath81 are so large that not all of the information contained in the dynamics is represented in the time series . this effect does not explain why a nonchaotic system may appear to be chaotic , that is , why the measured value of @xmath101 is greater than @xmath103 , as in ref . @xcite : this must be due to a problem with the above mentioned limits , and we defer a discussion of this point until after we have presented our results .    our results for each of the six models [ f , r , l][@xmath23,p4 ] are illustrated in fig . [ f : gpl ] with @xmath104 and in fig . [ f : gps ] with @xmath105 . the case @xmath104 samples the motion of the particle as it begins to diffuse among the scatterers . the models ( both chaotic and nonchaotic ) are indistinguishable , with the linear growth of the pattern entropy suggesting positive ks entropy , except that the superdiffusive model fp4 is showing signs that the ks entropy is zero , as a nonchaotic model should .    ( 330,530)(-70,0 ) ( 0,360 ) ( 170,360 ) ( 0,190 ) ( 170,190 ) ( 0,20 ) ( 170,20 ) ( 85,10)@xmath89 ( -65,95)@xmath106 ( -10,480)f ( -10,310)r ( -10,140)l ( 85,515)@xmath23 ( 245,515)p4    ( 330,530)(-70,0 ) ( 0,360 ) ( 170,360 ) ( 0,190 ) ( 170,190 ) ( 0,20 ) ( 170,20 ) ( 85,10)@xmath89 ( -65,95)@xmath106 ( -10,480)f ( -10,310)r ( -10,140)l ( 85,515)@xmath23 ( 245,515)p4    for @xmath105 ( fig . [ f : gps ] ) , where mostly ballistic behavior occurs , the graphs of pattern entropy vs. number of time steps is no longer linear . the r and l models look much the same , while the f models have an unusual feature when the pattern entropy is equal to @xmath107 which is due to the four available wind directions , that is , only a fourth of the trajectories at a given @xmath108 point are close in phase space after a certain time . the pattern entropy is becoming linear at larger times in all models , suggesting positive ks entropy . we can draw the following conclusions about the use of grassberger - procaccia type methods for this class of problems : ( a ) measurements need to be made at microscopic time scales , ( b ) even at microscopic time scales , the method does not seem to distinguish between chaotic and nonchaotic dynamics , and ( c ) the reason seems to be that unfeasibly long time series are needed in order to do that . [ [ macroscopic - measurements ] ] macroscopic measurements + + + + + + + + + + + + + + + + + + + + + + + +    it is clear that in an analysis as explained above , using large distance and time scales a diffusive system shows itself as almost independent of its microscopic dynamics ( whether chaotic , non - chaotic or stochastic ) , although we can not exclude the possibility that an unreasonably long trajectory might still give some evidence of weak correlations that survive for macroscopic times . in the brownian motion experiment , for example , there are of the order of @xmath109 collisions of the brownian particle with the surrounding solvent between each measurement of its position , so any correlations arising from , say , nonchaotic microscopic dynamics would be effectively averaged out between measurements . this is essential from the point of view of designing future experiments , but it is not the full story as it does not explain our results , which are obtained using short times . the brownian motion experiment could not investigate the motion at short time scales , but it is important to know whether data from an improved experiment might determine the microscopic chaoticity in principle . [ [ microscopic - measurements ] ] microscopic measurements + + + + + + + + + + + + + + + + + + + + + + + +    our results do not suggest that the grassberger - procaccia type methods can distinguish between chaotic and nonchaotic diffusive models even at short times because the nonchaotic r models consistently give the same results as the chaotic l models . the only distinctions we have been able to make are between diffusion and superdiffusion for @xmath104 , and between the continuous velocity space of the r and l models and the four velocity directions of the f models . both of these distinctions can be made from a time series without using the intensive analysis of the grassberger - procaccia method . [ [ how - long - a - time - series - is - required ] ] how long a time series is required ? + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    our wind - tree models appear unpredictable ( in the sense of positive ks entropy ) because it takes the system some time to `` realize '' that it is nonchaotic , in other words , quite a long time series length @xmath14 is required . for example , the most pronounced nonchaotic diffusive model we have is rp4 . once the particle `` knows '' the positions and orientations of all the scatterers , the system is completely predictable , and hence a zero ks entropy is manifest . however it is necessary for the particle to enumerate the large number of ways of achieving this combinations . the eight collisions could occur in any order , giving an extra factor of @xmath110 . thus our estimate is @xmath111 . note that collisions may not be able to occur in any order ( thus lowering the estimate ) , but we have ignored situations where the particle collides with three or four sides of a scatterer , or with the same sides more than once ( thus greatly increasing the estimate ) . the trajectory length would need to be substantially larger than the true number of ways the particle can determine the scatterer configuration , in order to obtain reasonable statistics , so that nonchaotic recurrences dominate . ] with sufficient statistics that nonchaotic recurrences dominate the gp calculation , for which our trajectory length of @xmath112 is apparently insufficient . for the infinite models such as r@xmath23 the motion of the particle is never completely predictable , but the proportion of collisions that the particle encounters new scatterers decreases as @xmath40 according to the discussion of sec . [ s : chaos ] . the measurement of such a slow decrease is far beyond the capabilities of any feasible implementation of the grassberger - procaccia method .    since the grassberger - procaccia approach does not seem to distinguish our chaotic and nonchaotic diffusive models , we turn to an alternative time series analysis method , with a view towards suggesting possible methods for analyzing future experiments that are able to sample the dynamics on microscopic time scales . such a method can not possibly enumerate all possible trajectory segments , as discussed above , but there are certain types of trajectory segments , namely those that are almost periodic , that stand out as predictable in a nonchaotic system . the method of the next section takes advantage of this property , and in effect distinguishes chaotic and nonchaotic systems by singling out this subset of all possible recurrences . we have seen in sec .  [ s : gp ] that the variant of the grassberger - procaccia method used by gaspard et al can not distinguish between nonchaotic and chaotic models , for the reason that the time series ( either experimental or numerical ) is not long enough to provide a reasonable approximation to the infinite time limit implied by the method . we can attempt to circumvent this difficulty by looking for specific types of recurrences , as opposed to taking an average over all types of recurring trajectory segments . in particular we will show that a promising candidate for a useful specific recurrence is given by one that corresponds to an orbit that is almost periodic ( again , within a spatial tolerance @xmath37 ) . we can identify these almost periodic orbits as those that repeat ( within a distance @xmath37 ) their previous motion at equally spaced intervals of time , hence the name `` almost periodic recurrences '' ( apr ) . if we compare the apr approach with the grassberger - procaccia ( gp ) method , there are a number of similarities and differences . both methods are designed to deduce dynamical properties from time series data . gp aims to make quantitative statements about the ks entropy , whereas apr as yet only provides a qualitative statement about chaoticity as expressed in periodic orbit properties , which are not equivalent to the usual definition of chaos as positive ks entropy . both methods make use of the probability of recurrences . gp then averages over many probabilities , while apr singles out a few especially significant recurrences . this means , for example that in an intermittent system , where the dynamics switches irregularly from chaotic to regular behavior , gp measures only the average , chaotic , dynamics , while apr is sensitive to the lack of chaoticity in the regular regions of phase space . in other words , an intermittent system can have a positive lyapunov exponent , and also periodic orbits with nonchaotic properties . the apr method characterizes such a system as nonchaotic , which is an appropriate designation with regard to diffusive properties as we observe in sec . [ s : esc ] below , but which is inappropriate from the usual point of view , which is that chaos is defined as positive ks entropy . an example of such an intermittent system is given in the final discussion . we now discuss the properties of periodic orbits in chaotic and nonchaotic systems ( particularly our models ) , before describing the apr method in more detail and presenting the results . periodic orbits are of great importance in chaotic systems . although typical trajectories ( as defined by the liouville or natural measures ) are not periodic , there are many systems for which they are approximated by periodic orbits when the periodic orbits are dense in the phase space or on an attractor . this allows many properties to be computed from expansions involving periodic orbits  @xcite . the properties of the periodic orbits are not directly related to the other dynamical properties discussed in sec . [ s : chaos ] , yet they do tend to differ between chaotic and nonchaotic systems , and so they constitute additional dynamical properties related to chaos . the properties of the periodic orbits in the lorentz gas are of most interest to us , and in some sense typical for chaotic systems , so we focus on these . in this and the following section we examine two relevant properties of periodic orbits , and sketch proofs of these properties .    [ [ existence ] ] existence + + + + + + + + +    there is an easy constructive proof of the existence of periodic orbits in the lorentz gas : choose two disks arbitrarily ; the line joining their centers gives a periodic orbit unless there is a disk interposing ; consider one of the original disks and the interposing disk , and repeat until no disk interposes . [ [ stability ] ] stability + + + + + + + + +    the periodic orbits are exponentially unstable , that is , almost all trajectories beginning a small distance @xmath37 from the periodic orbit at time @xmath113 are at a distance @xmath114 at time @xmath34 , where @xmath115 is the maximum ( here the only positive ) lyapunov exponent , which depends on the periodic orbit @xmath116 . this exponential instability makes it very unlikely for a typical trajectory to remain near a given periodic orbit for long times . we study here the properties of periodic orbits of our wind - tree models , which are typical of nonchaotic systems , without attempting to classify all possible nonchaotic periodic behavior . for a complementary study of the periodic orbits in polygonal billiards , see ref . @xcite . not all wind - tree models have periodic orbits . the argument given above for the lorentz gas fails because an orbit connecting two square scatterers is periodic only if the scatterers have the same orientation ( of zero probability in the randomly oriented model ) and the particle moves perpendicular to the surface of the scatterers ( excluded by definition in the fixed orientation model ) . in general , the existence of periodic orbits depends on the locations and orientations of the scatterers . we will find that using generic locations of scatterers periodic orbits exist in the r@xmath23 model , but not in the f@xmath23 model . the periodic models ( such as fp4 and rp4 ) are harder to treat mathematically , and are not discussed here . [ [ f - model - nonexistence ] ] f model : nonexistence + + + + + + + + + + + + + + + + + + + + +    the absence of periodic orbits for generic configurations of the f@xmath23 model follows directly from arguments in a proof of aarnes  @xcite . in order for the particle to return to its starting point , an expression linear in the positions of the scatterers must vanish . the coefficients of the @xmath108 position of a scatterer in this expression are integers depending on which faces the particle hits , and are zero only if the particle hits opposite faces an equal number of times . a simple _ reductio ad absurdum _ argument shows then that there is a scatterer ( perhaps more than one ) , the furthest to the top and right ( largest value of @xmath117 ) , which has collisions on only one of its faces ( the lower left ) , and hence the coefficients corresponding to the position of this scatterer are nonzero . if the scatterers are randomly placed , no linear combination of the positions can vanish , so no periodic orbit can exist . of course , there are special configurations of the scatterers ( such as at the corners of a square aligned with the coordinate axes ) that allow periodic orbits . since there are therefore generically no periodic orbits in the f@xmath23 model , the remainder of this section refers to the r@xmath23 model . [ [ r - model - existence ] ] r model : existence + + + + + + + + + + + + + + + + + +    the existence of periodic orbits for generic configurations of the r@xmath23 model follows from the observation that any set of configurations of a finite number of ( _ eg . _ 3 ) scatterers with nonzero ( lesbesgue ) measure includes configurations appearing somewhere in a generic infinite configuration . a period 3 orbit can always be found in an acute triangular billiard by minimizing the total path length , see fig . [ f : tri ] . if three scatterers outline an acute triangle and are positioned so that their faces contain the path of minimal length ( clearly of nonzero measure ) , a period 3 orbit exists . ( 280,280)(-50,0 ) ( 0,0 )    [ [ r - model - stability ] ] r model : stability + + + + + + + + + + + + + + + + + +    finally we investigate the stability of periodic orbits in the r@xmath23 model . the combination of the linear dynamics of the free particle with length preserving reflections leads to a linear separation of nearby trajectories . that is , an initial separation of @xmath37 in the direction of the velocity leads to a separation @xmath118 in the position after time @xmath34 . another way of saying this is that the number of uniformly distributed trajectories remaining within a distance @xmath37 of the periodic orbit is @xmath119 for large @xmath34 . thus a particle in the r@xmath23 model is quite likely to spend a long time near a periodic orbit , in contrast to in the lorentz gas . now we return to our method of almost periodic recurrences for distinguishing chaotic properties of diffusive systems from a time series . this consists of counting the number of almost periodic sequences in the time series , thus hoping to exploit the different properties of periodic orbits in chaotic and nonchaotic systems as discussed in the previous two sections .    as in the grassberger - procaccia method we begin with a time series containing @xmath11 positions of the particle @xmath84 spaced at a time interval @xmath104 , see sec . [ s : time ] . analogous to eq . ( [ e : dist ] ) we define a distance between two segments of length @xmath120 ( the period of an almost periodic orbit ) that begin at points @xmath85 and @xmath90 on the trajectory , @xmath121    using eq . ( [ e : dt ] ) we compute the number of initial points @xmath85 for which the orbit repeats within a tolerance @xmath37 , @xmath122\ ] ] where we use @xmath123 for the results presented in fig . [ f : isr ] . ( 330,470)(-50,0 ) ( 0,320 ) ( 150,320 ) ( 0,170 ) ( 150,170 ) ( 0,20 ) ( 150,20 ) ( 85,10)@xmath120 ( -15,95)@xmath124 ( -10,430)f ( -10,280)r ( -10,130)l ( 85,460)@xmath23 ( 235,460)p4    it is clear that the chaotic lorentz models , with an exponential decay of @xmath124 with @xmath120 can be easily distinguished from the nonchaotic f and r models which mostly have a slower and noisier decay , so we have achieved our objective of finding a time series analysis method that can distinguish between time series of chaotic and nonchaotic diffusive systems . it remains for us to explain the difference , and the plots obtained using the properties of the periodic orbits in these systems . for that , we note that there are two types of ( almost ) periodic orbits that might appear in the expression for @xmath124 : an orbit of length @xmath120 that has been ( almost ) repeated , and a shorter orbit of length @xmath125 that has appeared @xmath126 times , so that the total time is equal to @xmath127 . [ [ the - l - models ] ] the l models + + + + + + + + + + + +    in the lorentz gas , the probability of remaining near a periodic orbit for a time @xmath120 is of order @xmath128 ( see sec .  [ s : cpo ] ) , where @xmath115 is the maximum ( here the only ) positive lyapunov exponent , and depends on the periodic orbit @xmath116 . this is why both the infinite and periodic lorentz gases give an exponential form in fig . [ f : isr ] . because @xmath124 can be due to long orbits or repeats of short orbits , the exponential form depends on both ( 1 ) the exponential escape from short orbits with many repetitions and ( 2 ) the exponential instability of orbits with length , that is , that @xmath115 does not approach zero for long orbits . could become arbitrarily small . however , the probability of a long free path is also exponentially small , so the result is still exponential . ]    [ [ the - r - models ] ] the r models + + + + + + + + + + + +    in the randomly oriented wind - tree models , the probability of remaining near a periodic orbit for a time @xmath120 is of order @xmath129 ( see sec .  [ s : npo ] ) . this form is not immediately apparent in fig . [ f : isr ] because orbits ( particularly those surviving for long times ) can be counted more than once due to contributions from an interval of different initial points @xmath85 . both the r@xmath23 model and the rp4 model continue beyond @xmath130 due to increasingly rare long orbits and many repeats of shorter orbits . for example , the almost periodic orbit giving rise to the conspicuous point with period @xmath131 in the rp4 model is repeated many times , also contributing to multiples of the period , double @xmath132114,115 ; triple @xmath132172,173 ; and so on up to @xmath1321033,1034 . it is clear that the dynamics is completely different to that of a chaotic system . [ [ the - f - models ] ] the f models + + + + + + + + + + + +    in the fixed oriented wind - tree models , there are no periodic orbits ( see sec .  [ s : npo ] ) . this means that any contribution to @xmath124 is due to orbits that are close to periodic . the infinite model f@xmath23 contains many arrangements of a few scatterers that come arbitrarily close to generating a true periodic orbit , so the plot is similar to that of the r models . the periodic model fp4 , in contrast , has only a few short orbits that are at all close to periodic , and none that can be repeated to give further contributions at higher values of @xmath120 , or these would be observed in the figure . what we have presented here barely scratches the surface of the almost periodic recurrences method and its variants . it would be quite easy to search for 3 or more occurrences of an ( almost ) periodic orbit , or to distinguish between a single repeat of a long orbit and many repeats of a short orbit . we have restricted the presentation here to a demonstration of the method as a means to distinguish time series of chaotic and nonchaotic diffusive systems , but we hope that it will be useful in more general contexts . the case of absorbing boundary conditions provides a rich context for illuminating the differences between chaotic and nonchaotic diffusion . this is because the diffusion equation with absorbing boundary conditions leads naturally to an exponential escape process , while the periodic orbits of the randomly oriented wind - tree model require a power law escape process ( see sec . [ s : npo ] ) , thus violating the diffusion equation . we consider here the probability that a single randomly placed particle will remain in an open system for a given time , or equivalently , the number of identical noninteracting particles remaining in the system after the same time . for chaotic dynamics the escape rate formalism of gaspard and nicolis  @xcite uses the escape rate ( defined below ) to connect the diffusion coefficient to dynamical properties such as the ks entropy and the positive lyapunov exponents . for nonchaotic systems the escape process is qualitatively different depending on the properties of the periodic orbits as we have noted above , and the escape rate formalism must be generalized , if it is to make any sense at all . we now describe the escape rate formalism as it applies to chaotic systems , discuss the results for our models , and then attempt to generalize the formalism to include nonchaotic systems .    a macroscopic description of escape on a square ( for simplicity ; other geometries are analogous ) is given by the diffusion equation ( [ e : diff ] ) for the particle density @xmath43 , together with the boundary condition @xmath133 along the lines @xmath134 , @xmath135 , @xmath136 and @xmath137 . the general solution is @xmath138 with the decay rate @xmath139 at long times , the solution is dominated by the slowest decaying mode , corresponding to the escape rate @xmath140 note that the escape of particles is exponential in  ( [ e : diffsol ] ) . the escape rate formalism equates this macroscopic escape rate @xmath141 with the ( exponential ) escape rate of a ( microscopically ) chaotic system , which is related to its ks entropy @xmath74 and the sum of its positive lyapunov exponents @xmath142 ( here at most a single positive lyapunov exponent ) by  @xcite @xmath143 the `` escape rate formula '' which generalizes pesin s formula ( sec .  [ s : chaos ] ) to open systems . not only can the escape rate of chaotic systems be calculated from periodic orbit theory  @xcite , but , as we will show below , there is an intimate connection between periodic orbits and the escape process in nonchaotic systems . we now present our numerical results , from which we learn how the above theory for chaotic systems is modified in the nonchaotic case . as described in sec . [ s : comp ] , the arrangement of scatterers used is the same as in the periodic case with @xmath13 ; the absorbing boundaries are one distance unit from the edge of the periodic cell , leading to an open system of size @xmath144 . as before we use fixed oriented squares ( f ) , randomly oriented squares ( r ) or circles ( l ) . we denote absorbing boundary conditions by a , so the full notation for these models is [ f , r , l]a18 . we place @xmath145 particles uniformly ( without overlapping the scatterers ) in the square of size @xmath144 and compute the number of particles remaining in the system as a function of time , see figs . [ f : esc ] ,  [ f : esc2 ] . ( 300,320)(-50,0 ) ( 0,20 ) ( 165,10)@xmath34 ( -15,165)@xmath146    ( 300,320)(-50,0 ) ( 0,20 ) ( 165,10)@xmath34 ( -15,165)@xmath146    the chaotic lorentz gas exhibits exponential decay as predicted by the diffusion equation for all times , but the nonchaotic wind - tree models deviate from exponential decay at late times . the results are well described by the empirical expressions @xmath147 for the l , r and f models respectively over the range of times considered . the initial exponential decay may be compared with eqs . ( [ e : diffsol]-[e : er ] ) to give an effective diffusion coefficient . we find @xmath148 , @xmath149 and @xmath150 , which are certainly consistent with the results of the mean square displacement , tab . [ t : diff ] in sec . [ s : diff ] , given that eqs . ( [ e : diffsol]-[e : er ] ) depend on the macroscopic diffusion equations , and so require a large system @xmath151 limit . the coefficient @xmath152 is obtained by matching eq . ( [ e : diffsol ] ) to the uniform probability density of the initial conditions ; this comes to @xmath153 with @xmath154 as the initial number of particles in the system . from this close agreement , we conclude that the early escape is well described by the diffusion equation with the same diffusion coefficient that appears in the mean square displacement of sec . [ s : diff ] ; there is no trace of non - diffusive behavior at short times . [ [ the - r - model ] ] the r model + + + + + + + + + + +    we now discuss the long time behavior of the nonchaotic models . the randomly oriented wind - tree model has a @xmath36 decay in  ( [ e : r ] ) due to its periodic orbits ; particles that are initially close to a periodic orbit lead to this power law as described in sec . [ s : npo ] . the coefficient gives an estimate of the density of periodic orbits ; @xmath155 is quite small compared to the number of particles , @xmath145 , so periodic orbits are relatively rare in a sense that is difficult to define precisely . in terms of the escape rate formalism , a power law decay corresponds to an escape rate @xmath156 which trivially satisfies  ( [ e : erf ] ) , but yields no information about the diffusion coefficient @xmath28 , that is , ( [ e : er ] ) is not satisfied . [ [ the - f - model ] ] the f model + + + + + + + + + + +    the fixed oriented wind - tree model shows a complete escape of all the particles in a finite time , as might be expected from the lack of periodic orbits . the dramatic transition from the initial exponential decay to a much slower ( at first ) linear decay is somewhat surprising . we interpret this as follows : while there are no exactly periodic orbits in the f model , it is possible for the particle to remain in an orbit that is nearly periodic for some time , allowing the particle to remain in the system longer than the exponential decay would predict . if there is a bundle of trajectories that survive for just 800 time units , then this would lead to a linear law because particles are initially evenly distributed along the bundle of trajectories . if there was , in addition , a large bundle of trajectories that survive for , say 500 time units , there would be a kink in the plot , with a sharp decrease in gradient at @xmath157 . the observation that the curve is close to linear ( see fig . [ f : esc2 ] ) thus implies that the lengths of these long living trajectories are strongly peaked around 800 time units . one such long lived orbit is depicted in fig . [ f : traj ] ; its length is seen to be due to two almost periodic orbits . the complete escape of the particles corresponds to an escape rate @xmath158 , which flagrantly violates both  ( [ e : er ] ) and  ( [ e : erf ] ) and yields no information about the diffusion coefficient @xmath28 . ( 350,350)(-50,0 ) ( 0,0 ) ( 55,150)a ( 65,270)b    we conclude from these observations that , although the wind - tree models satisfy the diffusion equation when there are no boundaries , that is , the two time distribution function is of the correct gaussian form , such a macroscopic approximation is not valid when there are absorbing boundary conditions . this means we can not determine the diffusion coefficient in the same way as for the lorentz gas , eqs . ( [ e : diffsol]-[e : er ] ) , that is , from @xmath159 because the infinite time limit leads to zero or infinity as explained above . here @xmath160 is the probability of a particles remaining in a system of size @xmath9 for time @xmath34 ; this is obtained by observing many independent particles , and taking the limit of an infinite number of particles . we can , however , attempt to calculate the diffusion coefficient in the finite wind - tree models by taking the @xmath161 limit at the same time as the @xmath151 limit . this is because the diffusion equation is a good approximation even in the open case as long as the time is not too long , and the decay remains exponential . for example , assuming that the `` density of periodic orbits '' in the ra model is independent of @xmath9 ( which seems to be the case numerically ) , the transition time from exponential to power law decay , as determined by comparing the two terms in eq . ( [ e : r ] ) using eq . ( [ e : er ] ) is of order @xmath162 . this should be compared with the time scale for which the slowest decaying mode of the diffusion equation dominates , @xmath163 , see  ( [ e : gmn ] ) . thus , for a given @xmath9 , there is a narrow range in time @xmath164 in which the decay is exponential , and in which a limit can be taken to obtain @xmath28 from the escape process . for example we can combine the @xmath9 and @xmath34 limits by setting @xmath165 and @xmath166 , @xmath167 we expect this relation to hold for all the models [ l , r , f]a . we emphasize that eq . ( [ e : narrow ] ) , if it is mathematically valid , still does not provide a practical method for computing the diffusion coefficient beyond the estimates we gave following eq . ( [ e : f ] ) . this is because the range of time scales @xmath164 grows so slowly with @xmath9 that extremely large system sizes are required for more precise estimates of @xmath28 . we can also turn the argument around , and suggest that the escape at long times of finite systems may be a good experimental technique for detecting ( or ruling out at some level ) nonchaotic microscopic dynamics . in any case , the range of validity of the diffusion equation as a macroscopic description of nonchaotic systems is restricted . we have explored the connections between microscopic chaos and diffusion . on a superficial level , the details of the microscopic dynamics seem to have little effect on the diffusion process : gaussian diffusion is observed in our rp4 model , which has scatterers with flat sides ( hence no exponential separation of nearby trajectories ) and is periodic ( hence there is no source of randomness from a disordered environment ) . in addition , the grassberger - procaccia method of time series analysis can not distinguish our chaotic and nonchaotic diffusive models because it would require impractically long data sets . on a deeper level , however , the subtle differences between chaotic and nonchaotic diffusion are quite apparent if you know where to look : our time series analysis method based on periodic orbits has no trouble distinguishing between the chaotic and nonchaotic models , and the long time behavior of escape from an open system is also determined by the properties of the chaotic versus nonchaotic periodic orbits .    in the light of our results , we can make a few concrete suggestions for an experimental determination of the chaotic or nonchaotic properties in a diffusive system . firstly , as we remarked in ref . @xcite , it is necessary to make measurements on the distance and time scales relevant to the microscopic dynamics . in our models , this question is simplified by the fact that all microscopic time scales are of order one in our units , sec . [ s : time ] . once this is achieved , one could use an approach based on the method of almost periodic recurrences , which searches for periodic behavior ; alternatively the late time decay of particles from an open system can also elicit the chaotic or nonchaotic nature of the dynamics . this raises again the question of exactly what dynamical properties would be measured ; there are systems with the positive ks entropy of the lorentz gas as well as the power law unstable periodic orbits of the wind - tree models . we will now make some remarks about such systems , although we make no attempt to enumerate all degrees and classes of chaotic dynamics . our example in this discussion is a model containing both circular and randomly oriented square scatterers . the circular scatterers lead , as in the lorentz gas , to a positive lyapunov exponent , while the square scatterers lead , as in the randomly oriented wind - tree model , to power law unstable periodic orbits , using the same argument as in sec . [ s : npo ] . we have not simulated such a model numerically ; however the arguments used in the previous sections may be applied , leading to the following predictions :    sec . [ s : conn ] : there will be a positive diffusion coefficient with a gaussian distribution function , as in both the l and r models . [ s : gp ] : the grassberger - procaccia method will describe the dynamics as chaotic , due to a positive lyapunov exponent and hence positive ks entropy , as in the l models .    sec . [ s : isr ] : the apr method will describe the dynamics as nonchaotic , due to the power law unstable periodic orbits , as in the r models .    sec .  [ s : esc ] : such a system with absorbing boundaries will exhibit power law escape at long times , again due to the periodic orbits , as in the r model . the different conclusions reached by the gp and apr methods exemplify the fact that these methods are based on different dynamical properties ; for diffusion in an open system , it is clear that the mixed model is most similar to the r model , and hence that the apr method distinguishes the relevant dynamical property in this case . we conclude with some open questions and possibilities for further work . in terms of _ microscopic chaos _ , we have studied only a few models . we note that there are examples of lack of chaoticity we have not considered , such as the coexistence of chaotic and nonchaotic regions found in kam theory . in terms of _ diffusion _ , we note that in the light of sec . [ s : esc ] the diffusion equation is not a complete macroscopic description of nonchaotic diffusion . in terms of _ time series analysis _ , our apr method , while sufficient for our purposes here , could be much more developed and applied . with regard to _ periodic orbit theory _ , the methods developed for exploiting the importance of periodic orbits to compute properties of chaotic systems  @xcite do not apply to nonchaotic systems , and yet we have seen that periodic orbits also play an important role in nonchaotic systems . we note that our infinite models , while sharing some of the properties of corresponding finite billiards , appear to fall outside the domain of current mathematics . for example , a recent study of infinite billiard systems  @xcite is mostly restricted to cases with finite areas that are finitely connected ; our models obey neither of these conditions . it would be interesting to see if it is possible to rigorously determine the status of our models , particularly their ergodic properties , rate of decay of correlations and kolmogorov - sinai entropy within the framework of some mathematical theory of infinite nonchaotic systems . more generally , a mathematical understanding of statistical mechanics including the thermodynamic limit naturally leads to the study of infinite systems . finally , we remark that all the models we have considered  with the exception of one ( fp4 )  are diffusive , limiting our investigation as to the presence of microscopic chaos or not to such systems . in addition , the precise role played by microscopic chaos  as represented by the lyapunov exponents  and `` macroscopic chaos '' , as embodied by the randomly placed scatterers for the existence of a diffusion process and the value of the diffusion coefficient , remains open . a similar but more complicated situation obtains when diffusion of momentum ( viscosity ) or energy ( heat conduction ) and other transport processes are considered . first of all we want to acknowledge that h. van beijeren suggested the wind - tree model as an important alternative to the lorentz gas for the study of microscopic chaoticity and diffusion , which led to the work of ref . @xcite and this paper . we thank him also , along with l. bunimovich , j. r. dorfman , p. gaspard and i. procaccia for stimulating and helpful discussions . this work was supported by the engineering research program of the office of basic energy sciences of the us department of energy under contract # de - fg02 - 88-er13847 . we give here a calculation of the mean free time for our models quoted in sec . [ s : time ] . our derivation is due to chernov  @xcite , extended to the fixed oriented model , and with a technical caveat for the infinite models . refer to sec . [ s : comp ] for the definitions of @xmath168 , @xmath9 , @xmath14 and @xmath169 . the mean free time @xmath83 is equal to the mean free path , since the velocity is one . the mean free time is known exactly for billiard systems@xcite , which include the [ r , l]p models discussed here . we make a minor extension to allow the fp model ( which differs because we do not want to allow all velocity directions ) . we expect that the formula for the mean free time would still hold in the infinite models on physical grounds , but we can not justify this mathematically . briefly , the argument in ref . @xcite observes that the total volume of phase space @xmath170 can be computed in two different ways .    one expression for the volume of phase space is as an integral over the area of the billiard , giving @xmath171 for the original wind - tree model with fixed orientations , and @xmath172 for random orientations or the lorentz gas . here , @xmath170 is the volume of phase space , including both position and velocity , while @xmath173 is the area accessible to the point particle . for our infinite billiards we think of a large but finite periodic system of length @xmath9 , and formally take the limit @xmath151 at the end of the calculation . since the area @xmath173 is equal to @xmath174 , we have @xmath175 here as @xmath22 . the prefactor @xmath176 comes from the four possible wind directions , and the @xmath177 from integrating over all possible directions . putting these expressions together we have @xmath178 for fixed oriented squares and @xmath179 for randomly oriented squares or circles . an alternative expression for the volume of phase space is an integral over the boundary of the billiard , where the contribution from each point is given by the free path / time @xmath82 . in this way we find @xmath180 for fixed oriented squares and @xmath181 for randomly oriented squares or circles . here @xmath57 is the length of the boundary , that is , the total perimeter of all the scatterers , @xmath83 is the _ mean _ free path / time , and the numerical prefactor is the integral over the component of velocity perpendicular to the boundary . for the wind - tree models we have @xmath182 where @xmath14 is the number of scatterers , so @xmath183 ; for the lorentz gas , we have @xmath184 .    comparing the expressions for @xmath170 in both calculations , we find @xmath185 ( using our value of @xmath22 ) for the case of fixed oriented squares , @xmath186 for the case of random oriented squares , and @xmath187 for the lorentz gas . p. gaspard , m. e. briggs , m. k. francis , j. v. sengers , r. w. gammon , j. r. dorfman , and r. v. calabrese , nature * 394 * , 865 ( 1998 ) . c. p. dettmann , e. g. d. cohen and h. van beijeren , nature * 401 * , 875 ( 1999 ) . a. cohen and i. procaccia , phys . a * 31 * , 1872 ( 1985 ) . p. grassberger and i. procaccia , phys . rev . a * 28 * , 2591 ( 1983 ) . e. h. hauge and e. g. d. cohen , j. math . phys . * 10 * , 397 ( 1969 ) . h. van beijeren , rev . mod . phys . * 54 * , 195 ( 1982 ) . p. gaspard , _ chaos , scattering and statistical mechanics _ ( cambridge university press , cambridge , 1998 ) . w. w. wood and f. lado , j. comp . * 7 * , 528 ( 1971 ) . eckmann and d. ruelle , rev . mod . phys . * 57 * , 617 ( 1985 ) . young , ann . math . * 147 * , 585 ( 1998 ) . n. i. chernov and c. p. dettmann , physica a * 279 * , 37 ( 2000 ) . e. gutkin , j. stat . phys . * 83 * , 7 ( 1996 ) . m. wojtkowski , commun . phys . * 105 * , 391 ( 1986 ) . s. lepri , l. rondoni , and g. benettin , chao - dyn/9909004 . gaspard , adv . * xcix * , 369 ( 1997 ) . cvitanovi et al . d. biswas , pramana - j . phys . * 48 * , 487 ( 1997 ) . j. f. aarnes `` an ergodic problem for r - polygons '' dept . of mathematics , technical university of trondheim , norway ( unpublished,1974 ) . p. gaspard and g. nicolis , phys . lett . * 65 * , 1693 ( 1990 ) . p. gaspard and j. r. dorfman , phys . e * 52 * , 3525 ( 1995 ) . s. troubetzkoy , nonlinearity * 12 * , 293 ( 1999 ) . n. chernov , j. stat . phys . * 88 * , 1 ( 1997 ) .
Abstract:

Answer: we investigate the connections between microscopic chaos , defined on a dynamical level and arising from collisions between molecules , and diffusion , characterized by a mean square displacement proportional to the time . we use a number of models involving a single particle moving in two dimensions and colliding with fixed scatterers . we find that a number of microscopically nonchaotic models exhibit diffusion , and that the standard methods of chaotic time series analysis are ill suited to the problem of distinguishing between chaotic and nonchaotic microscopic dynamics . however , we show that periodic orbits play an important role in our models , in that their different properties in our chaotic and nonchaotic models can be used to distinguish them at the level of time series analysis , and in systems with absorbing boundaries . our findings are relevant to experiments aimed at verifying the existence of chaoticity and related dynamical properties on a microscopic level in diffusive systems .