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 dynamics and pattern formation of bacteria serve as paradigms in understanding many properties of multicellular interacting systems , such as collective behavior , self organization , evolution , and development @xcite . the characteristics of bacterial motility and aggregation depend on numerous biological and chemical parameters , which include light exposure , temperature , concentration of food or of other substances , and lead to the existence of many classes of cell motion . since a bacterium can not utilize electromagnetic or acoustic radiation to sense its environment , it must rely on physical contact and its motility may depend on the production , diffusion and degradation of a relevant set of sensed chemicals . these chemicals are called chemoattractants or chemorepellants if they attract or repel bacteria , respectively . examples of chemoattractants include food - sugars and amino acids - whereas antiobiotics , fatty acids or other noxious substances are chemorepellant . a prototype and well studied example of bacterial motility is the run and tumble mechanism used by _ escherichia coli _ in liquid environments . the specialized structures allowing for cellular motion are known as flagellae , spinning helical tails which extend from the cellular membrane into the surrounding environment : counter - clockwise rotation of the flagellar motor leads to a one - directional run whereas clockwise rotation causes tumbling along a random direction . during its run , e.coli _ cell periodically senses for chemosensitive substances , and , by comparing concentration levels in new and old environments , adjusts its motion and its likelihood to tumble in a new , randomly chosen , direction @xcite . if the cell is moving in the direction of increasing attractant , for instance , the probability of tumbling is lowered , so that the run time and ` mean free path ' length in this direction are increased . e.coli _ cell , thus , compares chemical concentrations at nearby points , so that its future motion is determined by a chemical gradient @xcite via the mechanism of _ chemotaxis may occur in other non - acqueous systems , such as in colonies of the bacterium _ myxococcus xanthus _ or of the eukaryote _ dictyostelium discoideum _ ( slime - mould ) which crawl on agar plates @xcite . the exact regulatory mechanism leading to the surface gliding of these rod - like and slow moving organisms is yet unknown , and responses to external stimuli appear to be more complex than in enteric bacteria such as @xmath0 @xcite . the cells undergo chemotaxis , but lack , or show very modest , responses to specific nutrient stimuli @xcite , suggesting that the chemotactic behavior is due mainly to self - generated signaling chemicals @xcite . the fact that cell density regulates typical reversal frequencies may be a consequence of one bacteria sensing the higher levels of signaling chemical in the presence of others @xcite . in particular , under starving conditions , the cells are driven by auto - released fibril trails to aggregate and form multicellular compounds called fruiting bodies . these compounds raise above the agar plate and eventually sporulate into new and more resistant organisms which are then released to the environment in search of new sources of nutrients . one major difference between these systems and enteric bacteria is that the chemoattractants embedded in or self - deposited on the agar plates are relatively immobile on the time scale of bacterial motility , whereas chemoattractant substances in acqueous environments have a finite diffusion constant . e.coli_ , the gliding cells can only sense what they come into immediate contact with , and since the chemosensitive substance diffuses slowly , large fibril concentration differences exist over the length of the cell . the gradient driving the chemotaxis mechanism is then determined between positions within a cell body length , as opposed to the large distances involved in enteric bacteria chemotaxis . another possible mechanism driving cellular motility is _ chemokinesis_. here , changes in the direction of motion or in the cell velocity are determined locally , without recourse to the determination of concentration gradients . chemokinesis may be a function of temperature , substrate adhesion , salt concentration or substances affecting the internal metabolism of the cell , slowing or enhancing its motion or turning frequency @xcite . the observed patterns of movement of _ m. xanthus _ and of slime - mould suggest that the motility of these systems may be affected by self - induced chemokinesis as well as self - induced chemotaxis . _ m. xanthus _ move by shooting out and retracting pili from the ends of their prolate bodies @xcite . in nutritionally abundant environments or under conditions of high population density , both cellular systems generate repulsive fibril trails , leading to cell dispersion and the colonization of new regions of the agar plate . the mechanism is hypothesized to be dictated by chemokinesis @xcite . previous theoretical studies into the collective dynamics of interacting motile bacteria , both of _ e. coli _ and _ m. xanthus _ type , include coarse - grained convective - diffusive transport type models @xcite and lattice - type simulations @xcite . these studies have considered the effects of interactions with concentration fields , modeling chemotaxis or chemokinesis or , as in ref.@xcite , both chemosensitive mechanisms . the continuum theories take into account turning rates , and changes in speed or direction of motion , however they generally ignore long time path history effects induced by the sensed chemicals . in these systems , the bias to the motion arises from an external , fixed stimulus and not a dynamical one . similar studies - with a fixed external trail - have been presented in different contexts , such as the aggregation of ants due to chemotaxis @xcite . other authors @xcite have presented systems in which the fibril trail that dictates the motion of the cells is determined over an extended area and is not sensed locally . the aim of this paper is to model and disentangle the effects of _ self - generated _ chemotaxis and chemokinesis , acting individually or concurrently and under diverse attractrant and repulsive conditions . in contrast to many of the previous studies , the chemical trail of our systems will emerge from the motion of the bacterial cells and evolve with the cellular motion itself . observations of individual organism paths may be more revealing of cellular dynamics than multicellular ensembles , and we restrict our attention to the behavior of a single bacterium . even if isolated , the dynamics of this single cell will be affected by its self - secreted fibril trail , and the resulting pattern will be that of a self - interacting random walker . we shall try to understand cellular motion by developing , and analytically solving , the equations of motion for an isolated _ m. xanthus _ cell under a mean - field approximation , and by means of simple monte - carlo simulations , incorporating chemotaxis and chemokinesis to various degrees . the underlying assumption of this paper is that the chemosensitive material acts as a regulator of cell motion affecting only the frequency of direction reversal , via chemokinesis or chemotaxis . the speed of the bacterium is assumed constant in either direction and is controlled by inherent cellular bio - energetics . experimental single cell trajectories show that an isolated cell travels along its axis , occasionally veering off , and for the sake of simplicity , we will only consider one - dimensional dynamics . in order to validate our models , cell particle tracking can be eventually performed in one - dimensional etchings or imprints into the agar plates . in this section we derive a general fokker - planck equation for the single one - dimensional random walker . as discussed in the introduction , we assume that the ` bacterium ' , or particle , travels with constant speed @xmath1 , and is subject to an _ ad hoc _ directional reverse mechanism that incorporates both chemotaxis and chemokinesis . to start , we ignore the effect of both phenomena , and we simply investigate the role of a constant reversal rate .     and length @xmath2 crawling in one dimension and depositing chemoattractant ( or repellent ) on the surface . the density of shading of the footprint is proportional to the total attractant deposited and arises from the hypothetical trajectory depicted by the arrows beneath . for the sake of simplicity , we have neglected the cell length @xmath2 in the depiction of the trajectory.,height=96 ]    the probability that a bacterium is centered at position @xmath3 at time @xmath4 with velocity @xmath5 will be denoted by @xmath6 and the probability of a directional switch in time @xmath7 is @xmath8 . if @xmath9 is independent of @xmath10 we can write :    @xmath11    in the above equations , we have included the possibility that @xmath9 depends on the type of reversal , since in general the reversal rates from one direction to the other are not the same . specifically : @xmath12    expanding eqns.([proa ] ) and ( [ prob ] ) and keeping only the @xmath13 terms , we find : @xmath14    the initial conditions are chosen so that at @xmath15 the particle is localized at @xmath16 with amplitudes @xmath17 satisfying the condition that @xmath18 : @xmath19    we obtain : @xmath20 \\ \nonumber & & \hspace{-1 cm } \left\ { \frac{a_+}{2 } [ \delta(x - v_0t)+\delta(x+v_0 t ) ] + \theta[v_0 t - |x| ] \right . \\ & & \nonumber \hspace{-1 cm } \left[\frac{a_+ q q_+}{4 v_0 } \sqrt{\frac{v_0 t+x}{v_0 t -x } } ~ i_1 \left(\frac{q_+ q}{2v_0 }   \sqrt{(v_0t)^2-x^2}\right )   \right . \\ & & \hspace{-1 cm } \nonumber + \left . \left . \frac{a_- \gamma_-}{2 v_0 } i_0 \left ( \frac{q_+ q}{2 v_0 } \sqrt{(v_0t)^2 -x^2}\right ) \right ]   \right\};\end{aligned}\ ] ]    @xmath21 \\ \nonumber & & \hspace{-1 cm } \left\ { \frac{a_-}{2 } [ \delta(x - v_0t)+\delta(x+v_0 t ) ] + \theta[v_0 t - |x| ] \right . \\ & & \nonumber \hspace{-1 cm } \left[\frac{a_- q q_+}{4 v_0 } \sqrt{\frac{v_0 t - x}{v_0 t + x } } ~ i_1 \left(\frac{q_+ q}{2v_0 }   \sqrt{(v_0t)^2-x^2}\right ) \right . \\ \nonumber & & \left . \hspace{-1 cm } + \frac{a_+ \gamma_+}{2 v_0 } i_0 \left ( \frac{q_+ q}{2 v_0 } \sqrt{(v_0t)^2 -x^2}\right ) \right ]   \right\}. \\ \nonumber\end{aligned}\ ] ]    here , the parameters @xmath22 , @xmath23 and @xmath24 are defined so that @xmath25 , @xmath26 and @xmath27 . the function @xmath28 is the heaviside function , @xmath29=1 if @xmath30 , and the @xmath31 functions are modified bessel functions . the @xmath32 function terms carry the probability of the random walker cell having performed no reversals , and the heaviside terms embody causality . the total particle probability distribution function is @xmath33 . under the symmetric condition @xmath34 , and @xmath35 ( or equivalently @xmath36 ) the solution reduces to that for the generalized smoluchowski diffusion equation determined by hemmer @xcite , which at large times is a spreading gaussian . the fourier transform of the distribution functions yields the moments of the random walker cell from which the cumulants can be extracted . we now consider the generation and sensing of chemoattactant or chemorepellant , and its effect on the reversal rates @xmath37 . _ m. xanthus _ cells deposit attractant chemical matter contained within a fibril slime under the area of cell contact . we assume that the release of the chemoattractant occurs at a constant rate and uniformly over the cell - substrate footprint . in general , for a cell of length @xmath2 , the fibril concentration @xmath38 obeys the following equation :    @xmath39 -r_d \phi(x , t),\end{aligned}\ ] ]    where the chemical is generated at a constant rate @xmath40 per unit area . the heaviside function @xmath28 in eq.([phi ] ) limits the increase in fibril concentration to the region strictly under the footprint of the bacterium centered at @xmath41 , whereas degradation and diffusion are included via the terms proportional to @xmath42 and @xmath43 . if the surface - deposited chemoattractant does not diffuse appreciably and @xmath44 is neglected , we find :    @xmath45 . \end{aligned}\ ] ]    note that @xmath38 depends on the entire history of the trajectory @xmath46 up to time @xmath47 . to couple the fibril concentration with cell dynamics , we must now include a relation between @xmath48 and @xmath38 . this relationship will depend on the particular mechanism of cell direction reversal and various choices will be discussed in the following subsections . it is important to note that , since we are assuming that fibril sensing is the main factor in determining direction reversals , we are also implicitly assuming that the sensing process , for example via a biochemical network within the bacterium , is much faster than the typical time for a cell to have moved appreciably . if we further take the molecular sensing apparatus to be distributed uniformly at the cell substrate footprint , the total chemoattractant @xmath49 sensed at time @xmath4 is just the local contribution integrated over the cellular length :    @xmath50    where @xmath51 is the width of the cell as depicted in fig.[1d ] . let us consider the case of chemokinesis first . here , the reversal rate depends only on the integrated fibril concentration . we will assume no intrinsic drift and set @xmath52 . it is also reasonable to expect the switching probability to saturate for large enough fibril concentration , when the cell can not respond to further increases in fibril levels . a plausible functional dependence for @xmath9 on the total sensed chemoattractant @xmath53 is of the type found in cooperative chemical binding , a michaelis - menten form with hill coefficient @xmath54 , so that for a cell at position @xmath55 , the reversal rate is :    @xmath56    where @xmath57 and @xmath54 , the inflection parameter and transition sharpness ( or hill coefficient ) , are peculiar to the system being considered . the effect of chemokinesis is measured by @xmath58 ; with the unbiased choice @xmath59 , @xmath58 is limited to @xmath60 . we must now solve eqns.([pro1 ] ) and ( [ pro2 ] ) with the above spatial and temporal dependence for @xmath61 , which in turn depends on the concentration @xmath62 via eqns.([phitotal ] ) and ( [ phint ] ) . the latter equations carry an explicit dependence on the past history of the walker . our approach will be to consider the @xmath63 contribution of eq.([hill ] ) as giving rise to eqns.([pro1 ] ) and ( [ pro2 ] ) with constant reversal rates , and find its associated green s function . the @xmath64 term in @xmath61 of eq.([hill ] ) then generates non - homogeneous differential equations , whose solution may be obtained via the known green s function . we will then be able to construct an iterative process to obtain @xmath65 as a function of all preceeding probabilities @xmath66 . to build the green s matrix function @xmath67 , we note that eqns.([pro1 ] ) and ( [ pro2 ] ) can be written as :    @xmath68    where @xmath69^t$ ] and    @xmath70    for the case of generic @xmath71 , independent of @xmath10 , the solution to this equation , with the initial conditions expressed by eqns.([initial ] ) is given by eqns.([pplus ] ) and ( [ pminus ] ) . in the case @xmath72 we refer to these solutions as @xmath73 and @xmath74 . the green s function for this problem stems from the modified matrix equation :    @xmath75    the fact that :    @xmath76    leads to the conclusion that the green s matrix function is given by :    @xmath77    this result can also be verified by direct evaluation of @xmath78 in fourier space . from eqns.(@xmath79 ) and ( @xmath80 ) it follows that reversal rates of the type @xmath81 yield a new matrix equation :    @xmath82    where :    @xmath83    the resulting equation for @xmath65 can be solved iteratively :    @xmath84    in the case of chemokinesis , the non - homogeneous term @xmath85 is contained in eq.([hill ] ) . in order to simplify the expression for @xmath62 , we assume unit cell width and take the limit of cellular length @xmath86 . we also neglect diffusion and decay of the fibril . under these assumptions , @xmath87 and @xmath88 . the problem is now a deterministic one , as exemplified by the delta - function in the expression for @xmath89 . in order to apply eq.([solution ] ) , which gives the statistical probability for one bacterium to lie at @xmath10 , we need to know its @xmath90 location at @xmath46 for all previous times . in other words , the complete history of the cell is needed to determine future motion and solutions can not be found .    nevertheless , we can approximate @xmath91 by an averaged density @xmath92 , utilizing a mean field theory approach . here , we average the probabilities @xmath93 for a single walker over many distinct independent realizations , so that the self interaction is to be expressed on the _ average _ of all replicas of the system , and does not depend on the history of individual walkers . we can now rewrite the equation for the fibril concentration as :    @xmath94   = r \ell \rho(x , t),\ ] ]    where we indicate the total bacterial density by @xmath95 . the above relationship signifies that , on average , fibril growth is proportional to the total density of bacteria at @xmath10 . integrating eq.([concentration ] ) and taking the initial condition @xmath96 , we obtain :    @xmath97    from which , using eq.([solution ] ) and the definition @xmath98 we can write :    @xmath99    here @xmath100 is a combination of the matrix elements of the green s function : @xmath101 and is even in @xmath3 . an equation similar to eq.([sdeltarho ] ) can be obtained for the total density @xmath95 , which depends on @xmath102 :    @xmath103    here , @xmath104 and is odd in @xmath3 , as required by the normalization of both @xmath95 and @xmath105 . it is useful to note that the above recursive equations may be used to evaluate the number density @xmath106 for a system of @xmath107 interacting bacteria as well . in this case , the distribution function @xmath108 and the evolution equations are the same as ( @xmath109 ) and ( @xmath110 ) provided we redefine @xmath111 . equations @xmath112 , ( @xmath109 ) and ( @xmath110 ) can now be solved numerically . once @xmath113 and @xmath114 are determined , @xmath115 can be evaluated and used to calculate the quantities of interest at subsequent times . the incorporation of chemotaxis into the determination of the reversal rate follows closely that of chemokinesis . let us consider chemotaxis as the sole reversal mechanism for a cell of length @xmath2 traveling at @xmath1 speed in the positive direction , and situated at position @xmath55 . the chemical gradient over the cellular length will drive the reversal rate . to specify the relation between gradient and reversal probabilities , we first introduce the probability terms @xmath116 :    @xmath117,\ ] ]    where @xmath118 is a constant that measures the strength of chemotaxis . for positive @xmath118 , the above expression translates into a low reversal probability whenever @xmath119 so that the cell is likely to keep moving speed at @xmath120 . this means we are modeling an attractant fibril trail . on the other hand , negative values of @xmath118 will represent a chemorepellant system . the same argument can be presented for a bacterium traveling at @xmath121 speed . in this case , the direction of motion tends to persist for negative gradients , hence the expression for @xmath122 , where positive @xmath118 values signify chemoattractant systems . in the limit of cellular length @xmath123 , the above expressions are written as :    @xmath124 % \\ % \label{htwo } % \mu_- ( x , t ) = e^ { + \sigma \phi_x ( x , t)},\end{aligned}\ ] ]    where the @xmath125 . we use this exponential function to capture the sensitive chemical signalling , such as that found in _ e. coli _ @xcite . another alternative for the function @xmath126 is to include a saturation in the form of @xmath127 in the exponent of ( [ hone ] ) . we assume the reversal probabilities @xmath128 at position @xmath10 to depend on @xmath116 through the following :    @xmath129,\ ] ]    so that for @xmath130 , in the absence of chemotaxis , the reversal probabilities are @xmath131 . we are now able to use the same formulation derived for the chemokinesis mechanism and apply it to the chemotactically driven case . again , we must resort to the mean field calculation of the densities @xmath95 and @xmath132 by means of the green s matrix function . the matrix @xmath133 appearing in the non - homogeneous term of eq.([nonhomo ] ) now has components @xmath134 given by :    @xmath135    and eq.([solution ] ) still holds . note that since @xmath136 , then @xmath137 . repeating the same type of calculations as for the chemokinesis case , for @xmath138 under chemotaxis we find :    @xmath139    here the previous history of the cell is contained in @xmath140 through the derivatives of the fibril concentration @xmath141 . the @xmath142 function is the same as defined in the case of chemokinesis . similarly , the density @xmath95 is :    @xmath143    and @xmath95 is properly normalized by the fact that @xmath144 ( defined in the preivous subsection ) is odd with respect to @xmath3 . if the two mechanisms of chemokinesis and chemotaxis are active , both @xmath58 and @xmath118 are non - zero , and @xmath128 must include both contributions . for @xmath145 , we write :    @xmath146 -1 } { p_k ( x , t ) [ \mu_\pm(x , t ) - 1 ] + 1 } \right\}.\ ] ]    here , @xmath147 is the contribution from chemokinesis in the form expressed by eq.([hill ] ) with @xmath148 , i.e. @xmath149 from the chemokinesis case . the contribution from chemotaxis is @xmath150 as defined by eq.([hone ] ) . the choice @xmath151 sets @xmath152 , that is , the rate is dictated only by the chemokinesis mechanism . conversely , @xmath153 sets @xmath154 and @xmath145 reduces to the purely chemotactic expressions of eq.([chemorev ] ) . for both @xmath155 the reversal rate is @xmath156 . the corresponding @xmath157 and @xmath158 are contained in eq.([gp ] ) and yield the following recursive relationship for @xmath132 :    @xmath159.\end{aligned}\ ] ]    the corresponding equation for @xmath95 is :    @xmath160.\end{aligned}\ ] ]    here , the @xmath161 and @xmath162 functions depend on chemotaxis and chemokinesis via :    @xmath163[\mu_+(p_k-1)-p_k ] } , \\ \nonumber \\ { \cal v}&= & \frac{\mu_+(1 - 2 p_k ) } { [ p_k(\mu_+-1)+1][\mu_+(p_k-1)-p_k]}.\end{aligned}\ ] ]    we have used @xmath164 and suppressed the @xmath10 dependence . note that in the absence of chemotaxis , for @xmath130 and @xmath165 , @xmath166 and @xmath167 , leading to the pure chemokinesis case . in the absence of chemokinesis , for @xmath168 , @xmath169 , as defined in the previous subsection , and @xmath170 . we shall discuss the numerical solution to the equations for @xmath95 and @xmath171 in the various cases of chemokinesis , chemotaxis or both , in the next section . in this section , by means of numerical integration , we solve equations ( [ srho ] ) , ( [ sdeltarho ] ) , ( [ srhochem ] ) , ( [ sdeltarhochem ] ) , ( [ sdall ] ) , and ( [ srall ] ) for @xmath95 and @xmath102 . under chemokinesis , eqns.([srho ] ) and ( [ sdeltarho ] ) are volterra equations of the second kind in the @xmath4 variable , for which standard methods can be applied @xcite . in particular , we discretize both the temporal and spatial axis according to a uniform mesh , of spacing @xmath172 and @xmath173 , and solve the two coupled equations iteratively . let us assume that @xmath174 and @xmath175 are known , where @xmath176 represent time steps up to @xmath177 . at the subsequent time step @xmath4 , @xmath132 is evaluated by calculating @xmath178 and by adding to it the sum over all previous time steps @xmath179 of the integrand . of these terms , the last one , at @xmath180 , contains the kernel function @xmath181 , which , when integrated over space , extracts a term proportional to the quantity of interest , @xmath132 . the other contributions of the sum , for @xmath182 , are spatial integrals , that are evaluated as sums after the discretization over the @xmath183 variable . these sums contain the known values of @xmath184 , as well as @xmath185 which depends on the fibril concentration @xmath141 . the latter is just the integrated @xmath186 as follows from eq.(@xmath187 ) . in order to evaluate @xmath141 , a time discretized summation over the @xmath188 values is carried out . the assumption here is that the bacterium senses the environment before laying down new fibril , so that the concentration @xmath141 is evaluated not by the trapeziodal rule , in which the extrema of the integration region are weighed each as 1/2 , but by weighing @xmath15 with 1 and ignoring @xmath180 . in the case of chemokinesis , once @xmath102 is determined , @xmath95 , which depends on @xmath189 , can be calculated in the same way .    in fig.[kin2 ] we plot the results for @xmath95 for a system under chemokinesis . the meshes are set at @xmath190 and the maximum time step , in these units is @xmath191 . cell velocity is fixed at @xmath192 and the parameter @xmath193 . for chemoattractant fibril , or @xmath194 , the distribution is narrower than the bare @xmath195 case , also plotted in the graph . negative values of @xmath58 represent chemorepulsion and broaden the distribution , eventually leading to a plateau - like region . the existence of the plateau is due to the fact that for large negative values of @xmath196 the reversal rate at saturation is almost zero , leading to ballistic motion of the cell in the region where the fibril has been deposited . once the cell reaches the edge of the plateau region of fibril concentration , the reversal rate increases dramatically and the cell tends to reverse its motion and travel backwards , until the opposite side of the plateau region is reached . this behavior accounts for the uniformity of the distribution function in the central region . time evolution also tends to broaden - and consequently lower - the plateau region , since the bacteria tend to build up fibril at the edges , thus increasing the range of motion .    for the case of chemotaxis , @xmath95 , is easier to evaluate , since it can be calculated without recourse to @xmath132 , as indicated by eq.(@xmath197 ) . the procedure is the same as that outlined for the case of chemokinesis . we take the gradient appearing in the expressions for @xmath198 to be evaluated over a distance @xmath199 , so that for @xmath200 and @xmath201 of the mesh we have :    @xmath202 \right].\ ] ]    results are plotted in fig.[chem1 ] . as discussed in the previous section , positive ( negative ) values of @xmath118 indicate chemoattraction ( chemorepulsion ) . it is interesting to note that for large negative values of @xmath118 two peaks develop ( only one is seen from fig.[chem1 ] , the other is its symmetrical image in the @xmath3 axis ) , with the bacteria traveling towards unexplored regions of space , void of fibril . when both mechanisms are present , chemotaxis and chemokinesis may enforce or compete against each other in their localizing or spreading effects . chemoattractant values of @xmath203 for instance , contrast the ballistic tendency of @xmath204 , and viceversa , chemorepulsive values of @xmath205 oppose the localizing effect of @xmath206 . in figures [ comp1 ] through [ comp5 ] , we examine the interplay of the two mechanisms for several values of @xmath118 and @xmath58 .    for both @xmath207 the distribution functions are localized about the central position , whereas for @xmath208 two peaks traveling in opposite directions develop . the sharpness of these peaks increases with increasing @xmath209 , as can be seen by comparing figures [ comp3 ] through [ comp5 ] . and the localization position increases with increasing @xmath210 as each figure also shows . competing effects arise when @xmath58 and @xmath118 are of opposite signs . for large negative values of @xmath118 , the effect of positive @xmath58 values is simply to shift the mean of the spreading peak towards the origin ( figures [ comp5 ] ) . for smaller negative values of @xmath118 increasing @xmath58 tends to localize the distribution functions at the origin , until the traveling peak is washed out and a plateau region , typical of large positive @xmath58 values forms ( fig.[comp3 ] ) . finally , in the regime of @xmath203 and @xmath211 , the localizing tendency due to chemoattraction is offset by the chemokinesis - plateau tendency . unusual distribution function patterns may arise as in the case of @xmath212 or @xmath213 and @xmath214 or @xmath215 in fig.[comp1 ] . to analyze the statistics of the one dimensional dynamics of cell motion under chemotaxis and chemokinesis , we also implemented a stochastic simulation on a 1d lattice . we discretize space and time and assume the random walker cell to have finite length @xmath2 . at time @xmath15 , the single cell is positioned at the origin of a one dimensional track and there is no fibril substance present . fibril is constantly secreted by the cell and the dynamics of the particle position obeys :    @xmath216    where @xmath217 is the position of the center of the cell at time @xmath218 and @xmath219 indicates the direction it was traveling prior to time @xmath4 . the stochastic random variable @xmath220 is defined as :    @xmath221 + 1 & \quad \mbox{with prob . } \,\ , 1-\gamma_t \end{array } \\ \nonumber\end{aligned}\ ] ]    where @xmath222 is the probability for the cell to reverse its direction , before the time step at @xmath218 . the initial direction of motion is chosen so that @xmath223 is @xmath224 with equal probability . we simultaneously model chemotaxis and chemokinesis through the dependence of @xmath222 on the self secreted fibril concentration @xmath38 , as in the continuum case . discretizing equations ( [ phi ] ) and ( [ phint ] ) we obtain :    @xmath225    where @xmath226 is the amount of chemoattractant that has built up at location @xmath3 and @xmath227 is the total amount of chemoattractant sensed by the cell . both @xmath226 and @xmath227 naturally depend on the past trajectory of the cell , which we keep track of as the simulation proceeds . we let @xmath228 represent the probability of direction reversal based on chemokinesis alone derived from equation ( [ hill ] ) . this probability is then stretched towards 0 or 1 depending on the local normalized chemical gradient in chemotaxis , through the quantity @xmath229 . analogous to eqs . ( [ hill ] ) and ( [ hopone ] ) in the continuum theory , the kinesis and chemotaxis effects are represented by :    @xmath230,\end{aligned}\ ] ]    where @xmath231 and @xmath58 and @xmath118 are weight parameters subject to the constraints :    @xmath232    the parameter @xmath118 scales the effect size of chemotaxis , with no chemotaxis corresponding to @xmath233 ; while @xmath58 scales the strength of chemokinesis , with no chemokinesis corresponding to @xmath153 . at each step of the monte - carlo iteration , the total reversal probability is calculated as :    @xmath234.\ ] ]    in fig.[compare ] we compare the results from the monte carlo simulations and the mean - field approximation of the previous section , for small times , for a typical case of chemokinesis and chemotaxis . we can now analyze results of the monte - carlo simulations for longer times than in the analytic case : @xmath235 . at these time scales , the general trends found in the analytic case persist for most combinations of @xmath236 . the simulation results are summarized in a phase - diagram . positive values of @xmath58 and @xmath118 both cause attractiveness in the system and we expect suppressed spreading of the distribution function . negative values of @xmath58 and @xmath118 on the other hand cause repulsion and enhanced spreading . to quantify these effects , we compared the distribution function @xmath237 with a normalized gaussian , corresponding to the parameters @xmath238 . more specifically , we compare the integral of the reference gaussian at the curvature inflexion points to the integral of the probability distribution function between the same limits . if the integral of the pdf is less than that of the corresponding gaussian , the dynamics is classified as being suppressed . this definition is shown in fig.[phase](top ) . we explore the entire phase space @xmath58 and @xmath118 and find the crossover from one regime to the other . the enhanced or suppressed diffusive character of the pdf curves may evolve from that of the early time regime , but it is found that after @xmath239 simulation time steps the classification is stable . the phase diagram is shown in fig.[phase](bottom ) . for @xmath240 , both kinesis and chemotaxis effects suppress the diffusivity of the trajectories . only for @xmath241 can the spreading dynamics be enhanced . due to our choice of the functional forms of the effects ( eqs . [ hill ] and [ hopone ] ) , the dynamics are much more sensitive to @xmath118 that @xmath58 , nonetheless , the effects of kinesis induces an asymmetry in @xmath58 in the phase diagram . the curve approaches the limits @xmath242 asymptotically . in this paper we presented a mean field model and monte - carlo simulations for bacterial dynamics under the mechanisms of chemokinesis and chemotaxis acting concurrently . the bacterial motion is that of a one - dimensional self - interacting random walker . the fibril trail that governs the dynamics of the system is itself a dynamical quantity that depends on the past motion of the cells . we model the two mechanisms in terms of two characterizing parameters , @xmath118 and @xmath58 , which represent , respectively , the degree of chemotaxis and chemokinesis . the mean - field results agree with monte - carlo simulations in the limit of short times . for long times we find a phase - diagram in ( @xmath118 - @xmath58 ) space that separates enhanced or suppressed diffusion regimes . in contrast to the short - time numerics displayed in figs.[kin2]-[comp5 ] , which show pdfs indicative of enhanced diffusion for @xmath211 , the long time phase diagram exhibits suppressed diffusion provided @xmath118 is negative and @xmath243 is not too large . therefore , in the long - time limit , kinesis modeled by the logistic saturation , eq.([hill ] ) , results in suppressed diffusion in the absence of chemorepellent effects ( @xmath241 ) . we have investigated the effects of increased cooperativity in the biochemical signaling through the @xmath54 parameter . as can be seen from fig.[phase ] increased @xmath54 values do not significantly affect the phase - diagram curve . the model exhibits variations in shapes of the particle pdf . the characteristic pdf shapes depend on whether long or short time dynamics are considered . one of the possible applications of this work is for isolated a - motility ( adventurous ) myxobacteria cells whose dynamics is driven by self - secreted slime and does not require the presence of neighboring cells in direct contact @xcite . our analysis is restricted to a one - dimensional system . for a single cell constrained to move along a 1d agar track , the relevant parameters ( @xmath118 , @xmath58 and @xmath54 ) of our model can be tuned to identify the correct mechanism of motion . moreover , the motion of freely moving a - type myxobacteria cells is _ locally _ one - dimensional . our model does not account for contact interactions between cells , rather , cells interact indirectly via trails of self - secreted slime . this simplification enables us to study the details of the effects of the history - dependent fibril concentration . direct cell - cell signalling interactions on a 2d aggregating colony are known to lead to the formation of complex structures such as propagating rippling waves and spiraling fruiting bodies @xcite . although the investigation of contact interactions between cells is beyond the scope of this paper , our approach does not preclude the formation of propagating waves and associated rippling patterns in systems where fibril - mediated cell - cell interactions are included . the authors are grateful for the assistance of g. lakatos in improving the implementation of our simulations and numerical computations . md and tc were supported by the national science foundation through grant dms-0206733 .
Abstract:

Answer: motivated by observations of the dynamics of _ myxococcus xanthus _ , we present a self - interacting random walk model that describes the competition between chemokinesis and chemotaxis . cells are constrained to move in one dimension , but release a chemical chemoattractant at a steady state . the bacteria senses the chemical that it produces . the probability of direction reversals is modeled as a function of both the absolute level of chemoattractant sensed directly under each cell as well as the gradient sensed across the length of the cell . if the chemical does not degrade or diffuse rapidly , the one dimensional trajectory depends on the entire past history of the trajectory . we derive the corresponding fokker - planck equations , use an iterative mean field approach that we solve numerically for short times , and perform extensive monte - carlo simulations of the model . cell positional distributions and the associated moments are computed in this feedback system . average drift and mean squared displacements are found . crossover behavior among different diffusion regimes are found . specialpagefalse oddhead_physical review e _ , * 68 * , 021925 , 2003 evenheadoddhead


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: viewing optimization problems that arise in computer science from the perspective of statistical mechanics has led to successful insights  @xcite . from this point of view , where optimization algorithms are intimately related to the energy or cost function , the features of the energy landscape are crucial in determining the success or otherwise of optimization . this connection is less evident in computer science since algorithms do not necessarily have such a direct connection with the cost function . one of the open questions of the field is to understand the extent to which the performance of algorithms can be determined by the character of the optimization problem and not the details of the algorithm itself  @xcite .    in this paper we hope to illuminate the issue by introducing a pair of models with the unusual feature that although they have the same static properties determined by the list of energy levels , the energy landscape is different since the notion of which states are close to each other is not the same . we perform dynamical simulations of relaxation after a quench to explore the energy landscape for each of the models following a similar strategy to that used to investigate other computer science problems such as satisfiability ( k - sat ) and graph colouring  @xcite . one of the best known optimization models is the traveling salesman problem ( tsp )  @xcite . states can be labeled by permutations since the paths in the tsp correspond to various orders of the cities visited . the combinatorial models we consider in this paper have states which are permutations of individual numbers rather than permutations of more general objects such as the set of city coordinates used in the tsp . in this case the natural optimization problem is sorting the numbers and we shall take that problem as sufficient motivation . the problem of sorting has been extensively treated in the computer science literature  @xcite and many efficient algorithms are known . most of this work is concerned with analyzing the time taken to perform the sort , though issues such a memory requirement and ability to use cache are also important .    for some optimization problems , such as tsp , the cost function is clearly the length of the path and this can be taken as the hamiltonian of the statistical mechanical model . the problem of sorting does not have an uniquely obvious cost function . we require that the lowest energy corresponds to the sorted state but have freedom to decide how to measure the degree other permutations are sorted . in this paper we consider two different energy functions that compute the degree of sortedness in different ways : one is similar to the tsp and measures the cost in terms of the length of a path , the other assumes knowledge of the sorted state and computes the distance from that state in a direct manner .    by investigating these two models we discover that the spectrum of energies is identical . this implies that the partition functions are also identical and all static properties will be the same . yet the hamiltonians are different , and the energy each hamiltonian assigns to a given state is not the same . we make some preliminary observations on the mapping between states with a given energy according to one hamiltonian and the same energy according to the other hamiltonian , but do not delve into mathematical details . of more concern to this presentation is the fact that although no physical distinction between the two models is visible at the static level , it is manifest in the dynamics . we investigate the local dynamics after a deep quench and show how it displays very different behavior for each model : in one case the ( sorted ) ground state is eventually found , in the other case it is not . the difference in behavior is identified as being due to rather different energy landscapes , with stable local minima in one hamiltonian but not the other .    from a physical point of view , study of the dynamical relaxation after a quench falls under the topic of phase ordering kinetics  @xcite . from the perspective of computer science , it is the natural way to study the efficiency of a local search algorithm in finding the optimum solution to a problem for statistical mechanical approaches to sorting , the states or configurations are the permutations of a set of numbers . in the following , we shall consider the set of @xmath0 integers @xmath1 . a more general model , in which the numbers are taken to have randomly chosen real values , is briefly described in appendix [ adisordered ] . most of the properties investigated in this paper hold for both models , but for the sake of clarity we shall work exclusively with the simpler integer model . the states correspond to all possible permutations @xmath2 of the set of @xmath0 integers . the permutation label can be regarded either passively as an ordered @xmath0-tuple or actively as the mapping needed to get to that @xmath0-tuple from the identity  @xcite . generally we shall employ the first interpretation , so @xmath3 signifies the @xmath4th element of the permutation , ( @xmath5 ) , but the second interpretation will be convenient later when we use it to write the permutation in terms of cycles . we imagine the @xmath3 as analogs of spins on a line , and sometimes we will refer to this line as a one dimensional spatial direction . the size of the state space grows as @xmath6 in comparison with typical ( ising ) spin models where the space grows as @xmath7 . it is well - known that models such as this , where the state space grows faster than exponentially , have difficulties of interpretation since scaling of the temperature or energy with @xmath0 is necessary to ensure that certain quantities are extensive . in this work , we do not investigate the phase structure , and never need to perform this scaling since we only consider zero ( or infinite ) temperature . we refer the reader to mzard and parisi  @xcite and to anderson and fu  @xcite for further discussion of this matter . we consider two energy or cost functions that introduce a measure of the distance of a permutation away from the identity permutation . the energy is lowest and vanishes when evaluated for the ordered or identity permutation . in writing these expressions it is convenient to introduce the analog of the kronecker delta overlap between individual spins as the one dimensional distance metric between the numbers @xmath3 : @xmath8    the first energy is familiar as the cost function for the tsp  @xcite , here evaluated for the one - dimensional case . the boundary conditions are slightly different from the usual tsp since we insist they be fixed rather than periodic . the connection with sorting is clearer with fixed boundary conditions and we need not consider the degeneracy of all states under cyclic shifts or inversions . however , it should be noted that the bulk part of the energy minimizes to either ascending or descending order , and it is only the boundary conditions that select ascending order . this energy , which henceforth we shall call the tsp energy ( @xmath9 ) , is :    @xmath10    the first form is written in the standard form for tsp and we have additionally assumed : @xmath11 , @xmath12 for any @xmath2 . in the second form the constant term is removed , by writing the end terms of the sum explicitly . instead of using the function ( [ eq:1dmetric ] ) , powers ( notably quadratic ) of the node position differences could be considered , but these energy functions do not appear to be natural in this context and do not obey the properties we will demonstrate below . configuration : the permutation 1364725 . in the top diagram , the tsp representation is shown . note the dependence on boundary conditions and that there is a continuous path from one wall to the other . the second diagram interprets the displacement energy as bipartite matching : each line could be regarded as a spring pulling the relevant point to its sorted location . the final diagram is a different representation of the displacement energy as an assignment problem using just one set of nodes . in this representation , the permutation splits up into several separate pieces corresponding to its disjoint cycles . ] the energy function for our second model , which we shall term the displacement energy due to its connection with the `` total displacement '' of r.w . floyd@xcite , is given as a sum of site distances from the sorted configuration @xmath13 . @xmath14    this form of cost function is familiar as an `` assignment problem ''  @xcite . the choice of distance measure relying on site differences is by no means unique . many other ways of defining the distance between configurations are possible . for example , in the physics literature , an overlap based on counting similar links is common  @xcite . another overlap ( that happens to yield a tractable model  @xcite ) is based on the matching problem and is determined by counting the number of sites that are in their correct relative positions . the different cost functions can be interpreted geometrically as illustrated in figure [ fig : lines ] . for the permutation , 1364725 with @xmath15 , used as the example in the figure , the energies may be computed as : @xmath16    these expressions and the diagrams in the figure make it clear that the energy may be computed as a sum over runs ( either ascending or descending ) in the case of the tsp energy , and as a sum over cycles in the case of the displacement energy . a run is the term used in the combinatorial literature for a subsequence of adjacent elements that are in sorted ( or antisorted ) order . for the tsp energy each contribution is the difference between the maximum and minimum value contained in the run . however , the contribution of each cycle to the displacement energy is not usually a single term ( for more complex permutations than that shown in the figure ) and it is necessary to look at runs within a cycle . the energy spectra of @xmath9 and @xmath17 are identical . that is , there is a one - one map between the full list of energies for all states computed with @xmath17 and the list computed with @xmath9 . of course there is a shuffling in the way the states are associated with energies .    . the multiplets and their energies are shown in the centre the left column shows the states according to the tsp energy and the right column shows them according to the displacement energy . states that are marked with a `` @xmath18 '' take different energies depending on the energy function used . ]    for example , figure [ fig : spectrum ] shows the spectrum of the @xmath19 model , with states labeled for each energy function . for many , but not all of these @xmath19 states , the energy is the same whether evaluated with the tsp or the displacement energy . the proportion of states with invariant energy decreases at larger @xmath0 . there are many degeneracies between energy levels : in the figure there are a 9-plet , a 7-plet , a 4-plet and a triplet besides the singlet ground state ( which is always unique ) that together make the @xmath20 states . in appendix [ amultordered ] we list the multiplicity structure for small values of @xmath0 . this information might be expected to lead to an expression for the partition function . however , we have been unable to obtain a general formula for these multiplicities and according to knuth@xcite , the generating function does not appear to have a simple form .      to demonstrate the relation between the two models formally , we present a mapping between the states @xmath21 , ( in fact a permutation in the state space ) that has the property that @xmath22 . the mapping relies on the representation of the state @xmath2 in terms of a permutation of the ordered state @xmath23 using the cyclic representation . this mapping is well - known and is for example described in knuth  @xcite . we write any permutation in terms of @xmath24 distinct cycles ( singleton cycles included explicitly ) . each cycle is labeled by a superscript @xmath25 and contains @xmath26 elements in the cycle as : @xmath27 we fix the freedom available in the way the cycle form is written by requiring that :    * the largest element appears at the start of each cycle : @xmath28 for all other @xmath29 in the cycle @xmath30 . * cycles appear in the order determined by their first elements : @xmath31 for @xmath32 . note that this is slightly different from the convention defined by knuth . the action of this permutation on the ordered state @xmath23 is to take the value @xmath33 to @xmath34 ( or cyclically ) . the resulting permutation is the original state @xmath2 of the map . the state @xmath35 obtained from the map is given simply by removing brackets from the definition of the element ( [ eq : cycles ] ) above and regarding the list of numbers as a permutation . for example with @xmath19 , the cyclic form @xmath36 takes @xmath37 to @xmath38 . thus @xmath2 = @xmath38 and @xmath35 = @xmath39 . in this case , both @xmath2 and @xmath35 lie in the same multiplet according to either the tsp or displacement energy . this is not the case for @xmath2 = @xmath40 and @xmath35 = @xmath41 that appear at the top of figure [ fig : spectrum ] and correspond to the cyclic form @xmath42 . nonetheless , @xmath43 = @xmath44 . the displacement energy in the state @xmath2 is given by a sum over cycles as : @xmath45    the tsp energy can be evaluated in the state @xmath35 to obtain three terms , one from the elements at the end of the original sum in ( [ eq : hstring ] ) , one from the contribution of each cycle , and a term corresponding to the intercycle contributions . @xmath46    the cycle sums appearing in ( [ eq : efield1 ] ) and the last term in ( [ eq : estring1 ] ) are identical except for one additional term in the first equation . the total difference between the energies is given by : @xmath47    now , by using the requirements on ordering of the @xmath33 s stated above , we find that the modulus signs may be removed in each sum and the total vanishes . this does not occur for other choices of distance function , but does continue to hold for the model based on real numbers rather than integers . the map is invertible due to the requirements listed above . note that the map must leave the ground state unchanged . however , the multiplet structure is not preserved , states that appear in a certain multiplet according to one energy function may appear in a different multiplet according to the other energy function . this can be seen in the case of @xmath19 as the states marked with a `` @xmath18 '' in figure [ fig : spectrum ] . several distinct versions of the map exist , based on different conventions for the way the cycle form is written . an approach to understanding the symmetries of the system would be to combine a map and the inverse of a different version . we do not consider this approach here since it takes us too far afield from the aim of the present paper . here we consider statistics of the energy ( either tsp or displacement ) spectrum : the average energy and its fluctuation . these are averages from a combinatorial point of view in which all states contribute equally ; thermodynamically they are effectively at infinite temperature . an estimate of the large @xmath0 behavior is most easily obtained by writing the tsp energy in terms of @xmath48 , the average distance between nodes : @xmath49    where the approximation consists in ignoring correlations between the distances between different pairs of nodes . in this case , the pairs can be imagined as two independent points thrown at random in the interval @xmath50 $ ] , so the average distance between them is @xmath51 .    a more formal computation based on the displacement energy proceeds as follows , where @xmath52 indicates a sum over all permutations . @xmath53 the order of summation is exchanged in the second line after which the sum over permutations becomes simple since there are @xmath54 permutations in which @xmath3 takes a given value @xmath29 . the resulting sum can be performed using standard techniques and a similar argument holds for the second moment . @xmath55    the width of the energy distribution , @xmath56 , scales as @xmath57 and therefore becomes relatively more peaked at large @xmath0 . although the energy spectra of the two models is identical and the partition functions are the same , the models still have distinct properties . this would be apparent by studying the response to some external field that couples to states in the same way in each model . the displacement energy itself could be regarded as an example of this kind of additional term in the cost function . however , in the context of this paper and the problem of sorting , the distinction is best studied by looking at the dynamical properties of the models . we only consider local dynamics : that is the basic moves are adjacent transpositions . of course there is no physical basis to these models requiring locality , and in real sorting algorithms non - locality of the elementary moves is essential to obtain efficient sorting . furthermore , we only consider a dynamics associated with the energy function - namely the monte carlo metropolis algorithm . this is natural from a theoretical physics perspective ( even though it does not correspond to any physical dynamics  @xcite ) , but algorithms that have much less clear relationship with the cost function are common in computer science . indeed we could imagine a reasonable algorithm that selects random sites and transposes with the neighbor if they are out of order . the usual approach physicists have used to investigate computer science problems is to study the dynamical relaxation under local search algorithms  @xcite . from the physical point of view , the approach consists in studying the phase ordering kinetics after a deep quench  @xcite .    . an average over 1000 different initial states is taken with error bars that are too small to be shown . ] the monte carlo moves are local adjacent transpositions and at zero temperature this effectively constitutes a randomized steepest descent algorithm . the metropolis algorithm we shall use selects a trial site at random and transposes with its ( right hand ) neighbor provided this move reduces the energy ( or does not change it ) . we perform numerical simulations starting from a random configuration which typically has energy very close to the average energy computed in section [ avenergy ] . figure [ fig : omontecarlo ] shows the results of numerical studies of dynamics according to this algorithm . the two curves in the figure correspond to dynamics based on each of the energy functions we have defined . starting from a randomly chosen initial state , the plot follows the evolution of the energy averaged over many choices of this initial state . this situation corresponds to a quench from very high temperature to zero temperature . the model based on the displacement energy evolves in the expected manner : the energy slowly reduces and eventually reaches the sorted ground state . the other model , based on the tsp energy starts with a similar initial energy then rapidly reduces to a plateau value at which level it continues indefinitely . no further decrease in energy is evident , even for much longer runs , and the energy never arrives at the sorted ground state . this behavior can be improved somewhat by using simulated annealing but it remains extremely slow and still tends to get stuck . to @xmath58 are shown with axes scaled to demonstrate data collapse . the energy axis is scaled by @xmath59 and is plotted logarithmically . the time axis is scaled by @xmath60 . an average over 1000 different initial states is taken with errors that are too small to see on the figure . ] according to figure [ fig : omontecarlo ] , the displacement energy appears to reduce at an exponential rate . the log - log plot ( not shown ) is only able to substantiate this for the first part of the decay . the quality of this initial exponential decay , and the form of the subsequent deviations from this form are shown in figure [ fig : displedynscaling ] .    within the initial exponential decay region it is possible to measure the @xmath0 dependence of the exponential timescale . at short times this characteristic time scales as @xmath61 ( fits using times up to order @xmath61 give the exponent with accuracy @xmath62 , and value tending to 2.0 as fewer points are taken ) , so the time axis of the figure has been scaled as @xmath63 in order to collapse the data in the initial region . the deviation from exponential form and absence of data collapse in the later part of the data is a finite size effect . with the help of small quantities of data for very large sizes up to @xmath64 , careful measurements of the time required to reach fixed values of @xmath65 can be made . the way these times scale with @xmath0 shows a consistent trend towards an exponent of 2.0 at larger @xmath0 .    , @xmath66 , and @xmath67 are shown with axes scaled to demonstrate data collapse . the energy axis is scaled by @xmath59 . the time axis is scaled by @xmath59 . an average over 1000 different initial states is taken with errorbars that are the size of the marks . ]    certainly for any size that can be simulated , the total time to arrive at zero energy is affected by the finite size effects and has a scaling exponent larger than 2 . not surprisingly , the resulting sort is rather slower than achieved with standard sorting algorithms that have best case behavior increasing as @xmath68 .    for the tsp energy , figure [ fig : stringedynscaling ] shows a scaled plot for different size systems . in this case the time is only scaled by a factor of @xmath59 and the energy axis is not logarithmic . the stability of the plateau energy is very clear in this figure . a detailed investigation finds that after the effect of finite size effects are removed , the per site energy on the plateau is @xmath69 . the fact that the evolution timescales of the tsp and displacement models are respectively @xmath0 and @xmath70 indicates another surprising distinction between the dynamics of the two models . , though the form is the same for all configurations on the plateau ) of the evolution of the tsp energy . size @xmath71 averaged over 1000 different initial states . ]      throughout this work we have implicitly considered the indices to form a line . in this section we consider correlations along this line and refer to such correlations as spatial . studying the evolution of these correlators may illuminate the dynamical properties of the model since dynamical lengthscales may appear . by analogy with the definition of the distance metric ( [ eq:1dmetric ] ) , we define the spatial correlator at distance @xmath72 as :    @xmath73    up to terms relating to boundary conditions , @xmath74 is nothing other than the tsp energy . a small value of @xmath75 signifies a strong correlation , and in the fully sorted state it grows linearly @xmath76 . from arguments similar to those used to derive the average energy , it can be shown that the correlator takes a constant value of @xmath77 when averaged over randomly chosen configurations ( more precisely , @xmath78 ) . in the figures to follow , we only show @xmath79 up to values of @xmath80 . this is because for larger @xmath72 , only a small number of pairs appear in the sum ( [ eq : cofr ] ) and the indices of these pairs are often near the boundaries , thus making this region excessively dependent on boundary effects . as the configuration evolves according to the tsp model , the spatial correlator remains similar to its initial constant value @xmath77 corresponding to the initial random configuration , but some structure develops at small @xmath72 . once the plateau is reached , there is no further change in the spatial correlator and figure [ fig : stringspatialc ] shows its final form . the non - trivial structure has a clear size , of less than 10 units , and neither the size scale nor the shape of the structure depend on @xmath0 . this scale indicates the distance over which ordering can take place according to the tsp dynamical process . the mechanism that blocks further growth of the order beyond this scale is discussed below . averaged over 1000 different initial states . in order of increasing slope the lines are for configurations taken every 10000 monte carlo steps after the start of the simulation . there is no change after 40000 steps . ]    in the case of the displacement energy , the correlator evolves as shown in figure [ fig : dispspatialc ] . here , no detailed structure ever appears , and the correlator is always a straight line with a gradient that smoothly evolves towards the steepest slope corresponding to the final sorted ground state . there is no characteristic lengthscale over which ordering takes place and then grows . it would rather appear that the system becomes organized on all scales simultaneously .      in order to investigate correlations between configurations at different times , we use the same overlap that was employed for the displacement energy . @xmath81    other possibilities are certainly possible . for example , an overlap based on matching counts the number of positions where the number has not changed . the graphs based on this choice do not convey any significantly different information from those given below . shown for different waiting times for the evolution according to tsp energy . size @xmath71 averaged over 1000 different initial states . waiting time is zero for the top curve and increases by 200 in each lower curve until the bottom curve which is the same for any waiting time greater than 1000 . this lowest curve indicates that aging does not occur . ] the simplest correlation to measure is against the completely sorted ground state . this is none other than the displacement energy itself , and was shown ( for dynamics based on this energy ) in figure [ fig : omontecarlo ] . another familiar comparison is against the initial state . for the dynamics according to the displacement energy , the correlation with the initial state disappears rapidly ( within about @xmath82 steps ) , and we do not show a figure in this case . for the dynamics according to the tsp energy , we show in figure [ fig : stringtwotimec ] the two - time correlators for different waiting time plotted against @xmath83 in the conventional manner . for @xmath84 the correlation against the initial state is included in this figure . we have drawn two - time correlators not because aging occurs in this model - it does not , but as a convenient way of demonstrating two features : that the configurations on the plateau retain some correlation with the original random state , and to show that evolution is not frozen on the plateau . with this aim , the waiting times shown in figure [ fig : stringtwotimec ] are quite short . the final value of the @xmath85 curve ( about @xmath86 ) is much less than that associated with correlation between random configurations ( @xmath87 in this case due to the factor of 2 in the definition ( [ eq : ctime ] ) ) , indicating that not all information in the original configuration has been lost by the time the plateau is reached . this agrees with the result of the spatial correlation that indicates that only local modification takes place . the lowest two - time curve holding for all @xmath88 ( for size @xmath71 ) corresponds to waiting times that have reached the plateau . this curve is not constantly zero as would be the case if there was no dynamics happening on the plateau . the configuration continues to evolve , though the energy does not change . however , the correlation is bounded : irrespective of how long after the waiting time , the correlator never rises above a certain value ( about @xmath89 ) . this limiting value is independent of @xmath0 . the reason for this behaviour is the limited size of the flat directions that are identified below .    . the four states in the centre have tsp energy @xmath90 while those on the periphery have @xmath91 . local transpositions allow three possible moves between states and these are all shown for the local minima quartet . ] the fact that a simple algorithm such as zero temperature monte carlo is unable to find the ground state is hardly a surprise in optimization models . there are many statistical mechanical examples of this state of affairs , and the effect is usually ascribed to the features of the free energy landscape . in the k - satisfiability problem the difficulty of finding the ground state via a local search procedure is due to the proliferation of states which trap the search into metastable phase . eventually , a large fluctuation provides a means of reaching the ground state  @xcite . however , in our tsp model the reason is more prosaic and the origin of the effect is due to the presence of stable local energy minima . in this respect , it is similar to the xor - sat model , that also suffers from such minima  @xcite . though the two energy functions have matching energy levels they have very different energy landscape characteristics . the displacement version has no stable local minima , but the tsp version does . moreover , these minima appear to proliferate as @xmath0 increases and are never avoided .    for small values of @xmath0 the local minima of the tsp model may be found explicitly . no such minima exist for @xmath92 . for @xmath19 there are four minima connected by monte carlo moves as shown in figure [ fig : localmin ] . note that these minima are precisely the states that appear with a `` @xmath18 '' in figure [ fig : spectrum ] , and indeed the energy assignments of figure [ fig : localmin ] are inverted for the displacement energy , so in that case there are no local minima . configuration : 5 9 11 15 7 4 10 12 13 2 . ]    for larger @xmath0 , the trapping configurations appear as shown in figure [ fig : plateauconfig ] . they consist of alternating ascending and descending runs separated by _ turning pairs_. for an upper turning pair , each element of the pair is greater than either of the elements neighbouring the pair , with a similar definition for a lower turning pair . these trapping configurations found on the plateau , are slightly more ordered than random configurations that consist of ascending and descending runs separated by ordinary turning points . a turning point is a requirement on a subsequence of length 3 , whereas a tuning pair is a requirement on a longer subsequence of length 4 . the dynamical monte carlo process performs this small scale ordering ( as observed in the spatial correlators ) to arrive at the plateau configurations . the trapping configurations are stable local minima : flat directions correspond to the moves that interchange the two elements of the turning pairs , and all other transpositions raise the energy . the number of states in the trap is therefore @xmath93 . since the number of turning pairs grows linearly with the size of the system , the size of the trap grows exponentially with @xmath0 . in appendix [ aplateau ] , an estimate for the number of turning points ( same as the number of turning pairs ) for plateau configurations is derived , so the exponential growth is approximately @xmath94 . of course , many different traps exist , each with this typical size . the picture of the trapping configurations in figure [ fig : plateauconfig ] makes the turning pairs appear like domain walls . this is a reasonable interpretation since the bulk part of the tsp energy has two different minima with ascending and descending order and it is these phases that are separated by the turning pairs . the boundary conditions raise the degeneracy of the sorted and antisorted phases , but this effect never comes into play here since the domain walls are frozen and do not move after the initial relaxation . the local minima provide a basis for attempting to understand the value of the plateau energy observed under the dynamics above . a rather naive argument given in appendix [ aplateau ] , based on characterizing the typical length of ascending or descending runs in the permutation defining the local minimum state gives the per site value of @xmath95 . this should be compared with the numerical value of @xmath69 found in section [ sec : timescales ] .    , with error bars indicating an average over 1000 random initial states . note the ledge near @xmath96 , and that only for @xmath96 does the algorithm succeed in finding the zero energy state . ]      given the rather different dynamic properties of the two models , it is natural to consider the behavior of models interpolating between them . the interpolating hamiltonian is : @xmath97 where @xmath98 is a parameter in ( 0,1 ) . the average energy ( in the combinatoric sense at infinite temperature ) is independent of @xmath98 . all other quantities , such as the energy of individual states and the degeneracy pattern are different from either of the models previously discussed . depending on the value of @xmath98 , the dynamics is expected to be more like one or other of the two lines shown in figure [ fig : omontecarlo ] . one might hope for a transition to take place for some intermediate value of @xmath98 , but in fact the dynamics reaches a plateau for any @xmath99 . figure [ fig : interpolate ] shows how the energy value of the plateau varies with @xmath98 . for small @xmath0 this curve is composed of a series of steps , ending in a horizontal ledge near @xmath96 . for larger @xmath0 the curve becomes smooth and the width of the ledge shrinks , eventually the curve becomes a sawtooth . the question arises as to why stable local minima appear even for an infinitesimal perturbation away from the displacement energy . consider the state corresponding to the completely reversed permutation ( with the central two elements transposed in the case of @xmath0 even ) . the displacement energy of this state is flat with respect to any of the local moves , though a series of moves will eventually lead to a lower energy state . on the other hand the same state is part of a local minimum quartet with respect to the tsp energy . an infinitesimal addition of the tsp component is therefore sufficient to make the state a local minimum , albeit with infinitesimally small barriers . although this argument correctly describes the reason that trapping can take place for @xmath98 so close to 1 , the identification of completely reversed states as being responsible is incorrect since their energy is considerably higher than the plateau . considering statistical mechanical models based on states that are permutations of numbers , we have demonstrated a relationship between two models with distinct energy functions . statically their partition functions are identical since the energy spectrum is the same . dynamically there are substantial differences since one model has local energy minima and the other does not . this strange situation , that models with identical static properties have distinct local dynamics , is not a paradox . in this case it is clearly due to the shuffling of states that modifies the energy landscape by rearranging the energies of states that are close to each other . the static analysis showed that the energy spectrum was the same by demonstrating a one - to - one map relating states with the same energy in each model . this map continues to hold for the more general model based on real numbers rather than integers . dynamically , we simulated a deep quench and showed that the energy decay proceeds in a completely different way in each of the two models . for the displacement model , no dynamical lengthscale appears as the system is reorganised . characteristic timescales do appear and vary as @xmath61 , though with strong finite size effects . for the tsp model , there was rapid ( timescale varying as @xmath0 ) decay to a plateau as a result of reorganisation over small spatial scales of size less than 10 units , that did not completely destroy the corelation with the initial state . evolution continued to occur on the plateau , but was limited in range . this was interpreted in terms of trapping configurations with turning pairs . flat movement within the trap was still possible between about @xmath94 local minima states . optimization problems of most interest , such as k - satisfiability , are much richer than the models presented here . nevertheless , in the context of the one parameter family of interpolating models , we have shown that our simple local search procedure exhibits a transition ( at the very edge of the interpolation region ) from being able to sort the numbers to becoming trapped by local minima . we hope that these models provide a simple arena for studying the general question of the role of the energy landscape in the performance of local search algorithms . instead of the integer model discussed in the text , this more general model is based on the set of @xmath0 real numbers @xmath100 , @xmath101 , with each @xmath100 in the range @xmath102 $ ] . without any loss of generality , we take the @xmath100 s to be ordered ( @xmath103 for all @xmath104 ) as this makes the identity permutation of the indices correspond to the lowest energy , or sorted state . the @xmath100 s should be regarded as quenched random variables , and for this reason the model might be regarded as a disordered model in contrast to the integer ordered model . the disorder however , is not of the independent variety that has been considered for both assignment and tsp problems using a replica approach in @xcite , but rather corresponds to euclidean distances in one dimension . here , the replica approach becomes intractable since the average over disorder couples sites .    for the disordered model , the definitions of the two energy functions remain exactly as given on the first lines of equations ( [ eq : hstring ] ) and ( [ eq : hfield ] ) , however the distance function is replaced by : @xmath105 irrespective of the choice of @xmath106 s , there are degeneracies in the spectrum , but not to the extent found in the ordered model . for example , for @xmath19 , the disordered model has energy levels with multiplicities ( 9,4,3,3,1,1,1,1,1 ) . in the ordered model these combine to give ( 9,7,4,3,1 ) . the mapping between the energy levels according to each model continues to hold , however there is a small difference in the moments . here the moments are defined after additional averaging over the @xmath106 s . the mean energy is the same as in the ordered model by design - indeed it was this criterion that selected the range of the @xmath100 s to be @xmath102 $ ] . the second moment is however slightly larger than quoted in ( [ eq : esdav ] ) . @xmath107    even for individual instances of the disordered model , similar dynamical behaviour to that described in the text for the integer model is observed . after averaging over realizations of the @xmath100 s , the similarity becomes even closer . the multiplicities are most easily obtained using the displacement energy @xmath17 this may be evaluated for every permutation and @xmath108 denotes the number of times that displacement @xmath109 appears . for small @xmath0 ( up to @xmath110 ) , it is straightforward to numerically enumerate these coefficients , and the first few are given in the table below . an estimate of the value of the average in the plateau is obtained by taking into account the effect of the monte carlo moves . we consider a subsequence of four points and look at the effect of a transposition on the central pair of the four , in all 24 possible orderings . the four points are labeled 1234 , but this only signifies their relative magnitude . table [ tab : mcn4 ] indicates whether the transposition causes a positive , negative or zero change in energy and lists the change in number of ( internal ) maxima or minima .    using this table , we find that moves that do not change the energy do not alter the number of runs . moreover , there are six configuration that both reduce the energy and number of runs . these are the pair 1324 , 4231 , and the quartet 2314 , 1423 , 3241 , 4132 . if we take their naive weights from the initial random configuration , then such transpositions for every pair of adjacent points leads to an average reduction in @xmath111 of @xmath112 . so the final value after one monte carlo sweep is @xmath113 . the argument for the average value at a maximum is extended from 3 ( given in the example above ) to 4 points , where it reproduces the 3 point case when all contributions are included . however , some of these contributions are removed by the monte carlo move ( we only include those configurations with @xmath114 in table [ tab : mcn4 ] ) resulting in @xmath115 . overall we obtain the estimate @xmath116 , to be compared with the numerical value of @xmath117 . the approximations in this approach are quite drastic . first we ignored correlations between the number of runs and the value of @xmath118 for that run in formula ( [ eq : erun ] ) . then we just looked at the effect of one monte carlo sweep , used naive weights and ignored any influence of one monte carlo move on another . it is therefore surprising that we obtain such a reasonable estimate . indeed , numerical studies show that the estimates of both @xmath111 and @xmath119 are incorrect by about 10% .
Abstract:

Answer: motivated by the problem of sorting , we introduce two simple combinatorial models with distinct hamiltonians yet identical spectra ( and hence partition function ) and show that the local dynamics of these models are very different . after a deep quench , one model slowly relaxes to the sorted state whereas the other model becomes blocked by the presence of stable local minima .