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: discretized models of @xmath2dimensional euclidean quantum gravity , known as simplicial gravity or dynamical triangulations , have been extensively studied in the past decade . the two dimensional model was proposed quite some time ago in ref . @xcite , and generalized to three and four dimensions in ref . the model is defined by the grand - canonical partition function @xmath4 where @xmath5 and @xmath6 are the discrete cosmological and inverse newton s constants , @xmath7 is the number of @xmath2dimensional simplices , and @xmath8 is the canonical partition function @xmath9 the micro - canonical partition function , @xmath10 , is defined as @xmath11 where the sum is over an ensemble of combinatorial triangulations @xmath12 with a spherical topology and with @xmath13 vertices and @xmath7 @xmath2simplices . for @xmath14 these two numbers are independent . the symmetry factor @xmath15 of a triangulation @xmath16 is given by the number of equivalent labeling of the vertices . including this factor is equivalent to summing over all _ labeled _ triangulations @xmath17 .    as a statistical system simplicial gravity displays a wealth of intriguing features . in @xmath14 it exhibits a geometric phase transition separating a strong - coupling _ crumpled _ phase , where the internal geometry of the triangulations collapses , from a weak - coupling _ elongated _ phase characterized by triangulations consisting of `` bubbles '' glued together into a tree - like structure or branched polymers @xcite . both in three @xcite and four dimensions @xcite this transition is discontinuous .    in the absence of a continuous phase transition with a divergent correlation length any sensible continuum limit where a theory of gravity might be defined appears to be ruled out this has led to the studies of modified models of simplicial gravity , e.g.  by changing the relative weight of the triangulations with a modified measure @xcite , or by coupling matter fields to the theory @xcite . suitable modified the phase structure of the model changes ; the polymerization of the geometry in the weak - coupling limit  the elongated phase  is suppressed and a new _ crinkled _ phase appears @xcite . the crinkled phase is separated from the other two phases by lines of either a soft continuous phase transitions , or possible a cross - over ( a divergent specific heat is not observed ) . qualitatively the same phase structure , shown schematically in figure  1 , is observed both in three and four dimensions . both in the elongated and in the crinkled phase the canonical partition function is seen to have the asymptotic behavior @xmath18 which defines the string susceptibility exponent @xmath19 . in the crumpled phase , on the other hand , the sub - leading behavior is exponential , corresponding to @xmath20 . many of the properties of the model eq . ( [ part.grand ] ) , in the weak - coupling regime @xmath21 , are reflected by triangulations close to the upper kinematic bound @xmath22 . this is evident both from monte carlo ( mc ) simulations and from the analysis of a strong - coupling expansion ( sce ) of the partition function . for combinatorial triangulations , where each simplex is uniquely defined by a set of @xmath23 distinct vertices , the upper kinematic bound is @xmath24 . although this bound is only saturated in the thermodynamic limit , @xmath25 , for each finite volume there exist triangulations with a well - defined maximal number of vertexes @xmath26 . the ensemble of these triangulations , the _ maximal _ ensemble @xmath27 , defines the micro - canonical partition function @xmath28    in this paper we investigate the properties of the maximal ensemble , mostly restricted to three dimensions , using both mc simulations and a sce , i.e.  a direct enumeration of the micro - canonical partition function eq . ( [ part.max ] ) . this investigation is separated into two parts . in the first part we study the pure ( un - modified ) ensemble and show that @xmath29 effectively separates into _ three _ distinct series ( or two series in @xmath30 ) as in general there are several different volumes @xmath7 corresponding to a given @xmath31 . one of those series is an ensemble of _ stacked spheres _ , a particular simple class of triangulations which can be enumerated exactly . this we do in section  3 . the other two series can be considered as `` almost '' stacked spheres , i.e.  with one , or two , _ defects _ respectively . heuristic arguments , supported by numerical evidence , show that although the leading asymptotic behavior of the three series is identical , the effect of inserting a defect increases the exponent @xmath19 by @xmath32 ( or @xmath33 ) compared to stacked spheres . this is discussed in section  4 .    in the second part of the paper , section  5 , we study the model eq . ( [ part.max ] ) with a modified measure , @xmath34 here @xmath35 is the _ order _ of a vertex @xmath36 , i.e.  the number of simplices incident to that vertex . as the measure is modified , varying @xmath37 , we observe a soft geometric phase transition from the elongated phase to a crinkled phase for @xmath38 , in agreement with the phase diagram figure  1 . moreover , a scaling analysis of the transition is consistent with a third order transition ( @xmath39 ) .    in the crinkled phase the fractal structure is dominated by a gas of sub - singular vertices , i.e.  vertices whose order grows sub - linearly with the volume . the intrinsic fractal structure is characterized by a set of critical exponents : ( i ) a negative string susceptibility exponent @xmath19 , which decreases with @xmath37 . ( ii ) a hausdorff dimension which is either @xmath40 , or @xmath41 ( as @xmath42 ) , depending on whether it is measured on the triangulation itself or its dual graph . ( iii ) a spectral dimension , measured on the dual graph , which increases smoothly from a branched polymer value @xmath43 , at @xmath44 , to @xmath45 at @xmath46 . we have collected some of the more technical issues in two appendixes : in appendix  a we list the sce of the three different series , and in appendix  b we describe some complications in extracting @xmath19 from measurements of minbu distributions on maximal triangulations . a maximal triangulation is a triangulation which for given volume @xmath7 has the maximal allowed number of vertices @xmath13 . for a combinatorial triangulation this implies : @xmath47 here @xmath48 denotes the floor function , i.e.  the biggest integer not greater than @xmath49 . from this definition it follows that in general more then one volume will correspond to a given maximal vertex number . in @xmath50 this defines three different _ series _ , labeled by the @xmath51pairs : @xmath52 the first series @xmath53 has the smallest volume for a given number of vertices and we will call it the _ minimal _ series . likewise , in @xmath30 there are in principle four different series . however , as only combinatorial four - triangulations of _ even _ volume are allowed this reduces to only two : @xmath54    in the space of all triangulations , the different series are related through a set of ( topology preserving ) geometric changes , the @xmath55moves . in a @xmath55move , where @xmath56 , a ( @xmath57)sub - simplex is replaced by its  dual  ( @xmath58)sub - simplex , provided no manifold constraint is violated . these moves are used in the mc simulations and are ergodic for @xmath59 , i.e.  any two triangulations are related through a finite sequence of moves @xcite . = 5.0 in    in three dimensions there are two sets of moves : move @xmath60 inserts a vertex into a tetrahedra and its inverse , move @xmath61 , deletes a vertex of order four . in move @xmath62 a triangle separating two adjacent tetrahedra is replaced by a link connecting the opposite vertices . this is shown in figure  [ fig.move3d ] . these moves change the number of different @xmath2simplices in the following way : @xmath63 note that only the first set of moves changes both the volume and the number of vertices , and when applied to a maximal triangulation we stay within the particular series , i.e.  the different series are internally connected _ via _ this move . the different series are connected , for fixed @xmath13 , by the @xmath62 and the @xmath64moves . however , as it is not possible to apply move @xmath64 to the series @xmath53 ( as triangulations in this series have the minimal volume for a given @xmath13 ) , the only connection which the maximal ensemble has with the rest of the triangulation space is _ via _ move @xmath62 applied to the series @xmath65 . the relations between the different series , _ via _ the @xmath55moves , are shown in figure  3 . [ r][r]@xmath66 [ c][c]@xmath66   maximal ensemble , @xmath53 , @xmath67 and @xmath65 , _ via _ the @xmath55moves.,title="fig:",width=377 ]    [ moves ]      a @xmath2dimensional _ stacked sphere_. however , in @xmath68 this does not define a maximal ensemble as @xmath69 for every ( spherical ) triangulation . ] is a triangulation constructed by successivelly gluing together smallest @xmath2spheres made out of @xmath70 , @xmath2dimensional regular simplexes@xcite . gluing is defined as cutting away one @xmath2simplex in each triangulation and joining the two triangulations by the resulting boundary . gluing a @xmath2sphere corresponds to applying the ( @xmath71)move .    from the definitions of stacked spheres and the maximal ensemble it follows that stacked spheres belong to the minimal series @xmath53 . the inverse statement , although less trivial , is also true . this was proven , in three and four dimensions , in two theorems due to walkup @xcite :    for any combinatorial triangulation of a 3sphere the inequality @xmath72    holds with equality if and only if it is a stacked sphere .    for any combinatorial triangulation of a 4sphere the inequality @xmath73    holds with equality if and only if it is a stacked sphere .    from those theorems follows a non - trivial statement : the move @xmath74 , and it inverse @xmath75 , is ergodic when applied to the ensemble of stacked spheres . in contrast , for the non - minimal series this statement is not true . there exist triangulations , @xmath76 , @xmath77 , that can not be reduced to the minimal configuration of the corresponding series by a repeated application of the move @xmath78 . this we have verified numerically . as stacked spheres have a simple tree - like structure , it is possible to calculate the partition function restricted to this ensemble explicitly . in this section we present algorithms for doing so , both for _ labeled _ triangulations ( which enter the mc simulations ) and for the number of _ distinct _ or _ unlabeled _ triangulations . we also show how to generalize the algorithm for counting labeled triangulations to include a modified measure eq . ( [ part.measure ] ) . all the results in this section are based on ref . @xcite which deals with the enumeration of _ simplicial clusters_. a ( @xmath70)dimensional simplicial cluster is a simplicial complex obtained by gluing together @xmath79dimensional regular simplices along their @xmath2dimensional faces . comparing this to the definition of stacked spheres one sees that a @xmath2dimensional stacked sphere is the boundary of a @xmath79dimensional simplicial cluster where the outer faces of the simplicial cluster correspond to the simplices of the stacked sphere . so number of @xmath80dimensional clusters with @xmath81 @xmath79simplices corresponds to the number of stacked spheres with @xmath82 vertices , i.e.  @xmath83 simplices . the basic `` building block '' is the expression for the number of @xmath79dimensional simplicial clusters build out of @xmath84simplices rooted at a marked outer face . denoting this number by @xmath85 we have the following recursive relation  @xcite , @xmath86 introducing the counting series , @xmath87 from lagrange s inversion formula ( see for example @xcite page 147 ) it follows that @xmath88    to calculate the number of labeled simplicial clusters we note that once the outer face has been labeled , the remaining @xmath89 vertices can be labeled in @xmath90 distinct ways . the four labels on the marked face can be chosen in one of @xmath91 ways . including the @xmath92 factor from eq . ( [ part.micro ] ) , a normalization factor @xmath93 ( as used in ref . @xcite ) , and dividing by @xmath7 to adjust for the fact that we count rooted configurations , we arrive at the micro - canonical partition function for stacked spheres , @xmath94 the first 20 terms in the above series are listed in tables  2 and 3 in appendix  a for @xmath95 and 4 . this corresponds to the enumeration of all stacked spheres of volume @xmath96 and @xmath97 , respectively .    to investigate the leading behavior of the above series , we perform an asymptotic expansion of eq . ( [ wd ] ) in powers of @xmath98 . for @xmath95 and 4 this gives : @xmath99 in both cases , the leading asymptotic behavior is of the form eq . ( [ def.gamma ] ) with @xmath100 . moreover , for both series the corrections to the leading order are analytic in @xmath98 , i.e.  there are no confluent singularities . this behavior is essential for applying the particular variant of the ratio method we use to analyze the series in next section . in contrast , the enumeration of the unlabeled triangulations  the number of distinct graphs  is much harder as the symmetry factor of each triangulation has to be canceled . this is the main subject of ref . @xcite which provides an algorithm to construct the generating function @xmath101 where @xmath102 is the number of distinct ( @xmath2 + 1)dimensional simplicial clusters build out of @xmath103 simplices . we list the functions @xmath104 and @xmath105 in appendix  a together with the first 20 terms ( tables  2 and 3 ) . it is interesting to note that although the series of unlabeled triangulations has the same asymptotic behavior as the labeled one , i.e.  the same @xmath106 and @xmath19 , the sub - leading corrections are not analytic in this case . this is evident from an analysis of the singularities of the functions @xmath104 and @xmath105 . as a consequence , for unlabeled triangulations it is is impossible to extract the leading behavior using a series extrapolation method such as the ratio method ( as is demonstrated in figure  [ fig.zeta_c ] ) . it is possible to extend the above algorithm for enumerating labeled stacked spheres to accommodate the measure term @xmath107 . the additional complication is to keep track of the order of each of the vertices on the marked face . here we present a recursive algorithm to do this . however , for simplicity we consider only @xmath50 stacked spheres , i.e.  @xmath30 simplicial clusters . we denote by @xmath108 the sum over all triangulations rooted at a marked outer face , where @xmath109 is the order of the vertices at the root . each triangulation is weighted by a factor @xmath110 where the product is over all vertices _ not _ belonging to the marked face . thus the number of rooted simplicial clusters becomes @xmath111    & \hspace{30pt }    \times\bigl ( a_4+b_4+c_4+d_4 - 4 \bigr ) ^{\beta }     \\[4pt ] \nonumber    & \hspace{-55pt}\times\delta_{a_1 + b_1 + c_1 -2 , q_1}\ ;    \delta_{d_1 + a_2 + b_2 -2 , q_2 } \ ;    \delta_{c_2 + d_2 + a_3 -2 , q_3}\ ;    \delta_{b_3 + c_3 + d_3 -2 , q_4}\;,\end{aligned}\ ] ] with the normalization @xmath112 we have to sum over all possible vertex orders of the marked face , @xmath113 to obtain the micro - canonical partition function @xmath114 where @xmath115 is the weight of the smallest triangulation . although a solution of eq . ( [ enrecmes ] ) is not known in a closed form , it is possible to calculate the series @xmath116 recursively . a simple program written in fortran 90 using quadruple precision takes about 7 hours of cpu time on a 600mhz alpha station to compute @xmath116 for @xmath117 @xmath118 . each additional term in the series requires four times the amount of computation as the previous one . compared to a method of directly identifying distinct triangulations in mc simulations , discussed in next section , using the algorithm presented above it is possible to calculate more terms in the sce . however , method presented above is limited to the enumeration of stacked spheres , i.e.  to a sub - class of the maximal ensemble , and moreover the series has to be calculated separately for each value of @xmath37 . we will start our exploration of the maximal ensemble with the pure ( unmodified ) case @xmath44 . as stated , the partition function eq . ( [ part.max ] ) effectively separates into three distinct series , @xmath53 , @xmath67 and @xmath65 , corresponding to regions of the triangulation - space only loosely connected _ via _ the ( @xmath119)moves . in this section we would like to explore in more details the difference between those series , in particular their asymptotic behavior . the following heuristic argument gives some idea of what asymptotic behavior we can expect . the series @xmath67 and @xmath65 , which are close to but do not saturate the upper bound in walkup s theorem eq  ( [ walk3 ] ) , can be considered as stacked spheres with one , or two , defects , respectively . a defect is defined as an one application of the move ( @xmath120 ) , i.e.  replacing a triangle by a link . for a sufficiently large triangulation one naively expects that the number of ways such a defect can be inserted should grow linearly with the volume ( inserting defects into two distinct non - symmetric stacked spheres does not lead to identical triangulations ) . this implies that although the leading exponential behavior of the three series might be the same , the sub - leading behavior is should be multiplied by @xmath121 ( or @xmath122 ) , hence @xmath123 ( or @xmath124 ) , with @xmath19 defined in eq . ( [ def.gamma ] ) .    to investigate the asymptotic behavior , and to verify the statement above , we have performed several numerical experiments . in figure  [ fig.free ] we show the micro - canonical partition function for the three series measured numerically for @xmath125 . we have divided out the leading behavior of the partition functions , i.e.  we plot @xmath126 where @xmath127 is some convenient normalization and we assume the exact values @xmath128 the behavior in figure  [ fig.free ] is indeed compatible with this assumption . if the series had different @xmath129 the curves should approach a straight line with a non - zero slope . likewise , a different value of @xmath130 would manifest itself in a very slow convergence . neither behavior is observed , all three curves in figure  [ fig.free ] converge rapidly on a horizontal line . it should be noted , however , that the two non - minimal series appear to have larger finite - size corrections .      a more precise determination of the critical parameters , @xmath129 and @xmath130 , comes from the sce of the micro - canonical partition function eq . ( [ part.max ] ) . in section  3 we presented a algorithm for calculating the sce restricted to the minimal series @xmath53 . a corresponding recursive algorithm for constructing the two non - minimal series is , however , not available . in fact , it is not clear that such algorithm exist . unlike @xmath131 , the non - minimal series can not be constructed from some ( different ) minimal configuration by a successive application of move ( @xmath132 ) . we thus resort to an alternative method to enumerate the two non - minimal series . this method , based on a direct identification of all distinct triangulations in mc simulations , has previously been applied in four dimensions to identify all distinct triangulations of volume @xmath133 @xcite . the two main ingrediences of the method are : ( i ) a sufficiently complicated hash function @xmath134 that uniquely identifies combinatorially distinct triangulations . ( ii ) a calculation of the corresponding symmetry factors , @xmath15 , by an explicit permutation of the vertex labels . the sce s for @xmath135 and @xmath136 are shown in table  4 in appendix  a. as the two non - minimal series grow much faster than the minimal one , this method is limited to volume @xmath137 . this corresponds to 10 and 8 terms respectively . to calculate the next term would require the identification , and storing , of @xmath138 distinct triangulations  a formidable task with the computer resources available .    to analyze the sce and to extract the asymptotic behavior we apply an appropriate series extrapolation method . for these particular series the most powerful method is a variant of the ratio method , a neville - aitken extrapolation , in which successive correction terms to the partition function are eliminated by a suitable sequence extrapolation technique @xcite . for the minimal series this method yields impressively accurate estimates of the critical parameters . applied to the first 20 terms in table  2 we get @xmath139 compared to the exact values @xmath140 and @xmath141 the quoted errors indicate how the estimate change as the last term in the series is excluded .    for the non - minimal series the finite - size corrections are larger and , in addition , we have less terms available . the estimates of the critical parameters are correspondingly less accurate , especially for @xmath65 : @xmath142 these values are though consistent with the conjecture that @xmath106 is the same for all the three series , but that @xmath143 .    the same behavior is observed for the two series of the maximal ensemble in @xmath30 . the first 20 terms in table  3 yield equally impressive agreement with the exact values : @xmath144{3125/256 } \approx 0.625503\ldots$ ] and @xmath140 . and the first 7 terms in the non - minimal series @xmath67 , given in ref . @xcite , give @xmath145 and @xmath146 .    some final comments regarding the non - minimal series . although an exact enumeration algorithm is not available it is possible to obtain a lower bound on the series @xmath67 and @xmath65 by assuming that all non - minimal triangulations can be constructed by move @xmath60 from the appropriate minimal configuration . the resulting sce has , for both series this assumes that the two defects are connected , thus excludes triangulations with two randomly distributed defects . this results in a lower estimate @xmath147 , instead of @xmath148 which we expect for @xmath65 . ] , the same @xmath106 and @xmath147 . the different asymptotic behavior for the three series is very pronounced when one considers the distribution of minimal neck baby universes ( minbu s ) measured on maximal triangulations . a minbu is a part of a triangulation connected to the rest _ via _ a minimal neck  on a three dimensional combinatorial triangulation a minimal neck consists of four triangles glued along their edges to form a closed surface ( a tetrahedra ) . by counting in how many ways a `` baby '' of size @xmath149 can be connected to a `` mother '' of size @xmath150 , the distribution of minbu s can be expressed in terms of the canonical partition functions of the corresponding volumes @xcite .    for a maximal triangulation the situation is more complicated as the baby , the mother , and the `` universe '' can belong to any one of the different series @xmath151 . it is not even obvious that if a maximal triangulation is split along a minimal neck the two parts belong to the maximal ensemble at all . this is though easy to verify and @xmath152 , the result only belongs to the maximal ensemble provided @xmath153 . ] . in general the distribution of baby universes measured on maximal triangulations becomes @xmath154    & \sim & b^{\gamma(k_b)-2 } \ ; ( n_3-b+2)^{\gamma(k_m)-2}\;.\end{aligned}\ ] ] the indices @xmath155 , @xmath156 , and @xmath157 refer to the baby , the mother , and the universe , the factors @xmath149 and @xmath158 count the possible locations of the neck , and @xmath159 is a _ contact _ term due to interactions on the neck . the contact term is discussed further in appendix  b ; for pure simplicial gravity this is just a trivial constant . in the last step we have used the expected asymptotic behavior of the series eq . ( [ def.gamma ] ) .    in principle , eq . ( [ babymax ] ) defines nine different distributions , however the additional restriction , @xmath160 , reduces the number of possible combinations to six . two of those , @xmath161 and @xmath162 , both corresponding to @xmath163 , are plotted in figure  [ fig.baby3d ] . both distributions are measured on triangulations with @xmath164 , but on different volumes , @xmath165 , and 502 . the former distribution , @xmath161 , which is composed solely of the minimal series , has remarkably small finite - size effects . a fit to eq . ( [ babymax ] ) , including minbu s @xmath166 , gives @xmath167 . in contrast , for an acceptable fit to the distribution @xmath162 a lower cut - off , @xmath168 , had to be imposed . the result , @xmath169 , is consistent with our previous estimate . the remaining four distributions in eq . ( [ babymax ] ) correspond to a baby and a mother belonging to different series . however , the mixture of different @xmath19 s and large finite - size effects makes a reliable estimate of @xmath130 very difficult from these distributions . we now turn to the maximal ensemble modified by varying the measure eq . ( [ part.measure ] ) for @xmath170 . for most part we restricted our investigation to the ensemble of stacked spheres where we expect less finite - size effects . this provides the additional benefit that the mc simulations are much simplified as only move @xmath60 , and its inverse , are required for the updating process to be ergodic . however , we also did extensive simulations of the full model ( including finite values of @xmath6 ) to verify that this did not introduce any bias in the measurement process . typically we collected about 10 to 50 thousand independent measurements at each value of @xmath37 and on different volumes @xmath171 to 8000 .      as the fluctuations in the vertex orders are enhanced by modifying the measure , i.e.  by decreasing @xmath37 , the internal geometry changes . the polymerization of the triangulations is suppressed and the model enters a new phase akin to the crinkled phase observed in four dimensions @xcite . the nature of the internal geometry in the crinkled phase is discussed in section  5.2 , here we present what evidence we have for the existence of a continuous phase transition at a critical value @xmath172 .    in figure [ fig.trans]_a _ we show the fluctuations in the energy density , @xmath173 , as we vary @xmath37 . this is for @xmath174 to 8000 and we define @xmath175 . we observe on all volumes a peak in @xmath176 consistent with the existence of a phase transitions at @xmath177 . the peak value appears to saturate as the volume is increased ; this makes it difficult to distinguish between a soft continuous phase transition and a cross - over behavior . the scaling of the peak value is , however , consistent with a third order phase transition . this is shown in figure  [ fig.cvmax ] where we fit the maximal value to the scaling behavior : @xmath178 . the fit yields @xmath179 , with a @xmath180 . assuming hyper - scaling , @xmath181 , is valid this implies a specific heat exponent @xmath182 . additional evidence for a transition is provided by other geometric observables such as the maximal vertex order @xmath183 whose value jumps at the transition . as discussed in the introduction , and is quantified in the section  5.2 , the crinkled phase is characterized by a gas of sub - singular vertices , i.e.  vertices whose local volume grows like @xmath184 , @xmath185 . as a consequence both the fluctuations in the maximal vertex order , @xmath186 , and its energy derivative , @xmath187 , diverge at the transition . an example of the former is shown in figure  [ fig.trans]_b_. in fact , the susceptibility @xmath188 appears to diverge in the whole crinkled phase in contrast to its behavior in the elongated phase . note that some of the data in figure  [ fig.trans ] are measured on non - minimal triangulations ( e.g.  @xmath189 and 4000 ) . this shows that the phase transition is a generic behavior of the maximal ensemble , not just of the minimal series . the same is true for the nature of the crinkled phase explored in the next section .      to investigate the nature of the internal geometry in the crinkled phase , and to compare it to the elongated phase , we have measured several aspects of the fractal structure : the distribution of vertex orders @xmath190 , the string susceptibility exponent @xmath19 , and both the fractal and spectral dimensions @xmath191 and @xmath192 . the shape of the vertex order distribution @xmath190 is a good indicator for the nature of the internal geometry in the different phases . in the elongated phase the vertex orders are more or less equdistributed and the maximal vertex order grows logarithmically with the system size , @xmath193 . in the crinkled phase , on the other hand , vertices of large order are favored and the maximal vertex order grows sub - linearly , @xmath194 . estimates of @xmath195 , from finite - size scaling including volume @xmath171 to 8000 , indicate that @xmath196 for sufficiently large negative @xmath37 . examples of the vertex order distribution in the two phases are shown in figure  [ fig.vord ] for @xmath197 and @xmath46 . in the crinkled phase the distribution @xmath190 developes a `` bump '' in the tail , corresponding to a condensation of vertices of large order . the bump persists as the volume is increased , i.e.   it captures a finite fraction of the volume as @xmath198 . this suggests that the crinkled phase can be characterized as a _ gas _ of sub - singular vertices . when the ensemble of triangulations in eq . ( [ part.micro ] ) includes degenerate triangulations @xcite . however , in @xmath30 the bump is disconnected from the rest of the distribution  this justifies the nomenclature _ singular structure_. ]      to measure the string susceptibility exponent for various values of @xmath37 , and restricted to the minimal series , we have used several methods :    _ the pseudo - critical cosmological constant _ , @xmath199 + assuming the asymptotic behavior eq . ( [ def.gamma ] ) , a saddle - point approximation of the grand - canonical partition function gives : @xmath200    by measuring @xmath199 with high accuracy on volumes @xmath201 to 2400 , and for @xmath202 , -2 , -3 , and -4 , we determined @xmath19 by a fit to eq . ( [ eq.mu_pseudo ] ) .    _ the minbu distribution _ , @xmath203 + we measured the minbu distribution on volume @xmath204 and for the same values of @xmath37 . however , in contrast to the elongated phase , in the crinkled phase the contact term @xmath159 has a non - trivial dependence on @xmath149 and has to be eliminated prior to extracting @xmath19 from the ( corrected ) distribution @xmath205 by a fit to eq . ( [ babymax ] ) .    _ the strong - coupling expansion _ i + using the enumeration methods for a non - zero @xmath37 , presented in section  3.2 , we calculated the first 17 terms of the sce of stacked spheres for few values of @xmath37 . from this sce we extracted @xmath19 using the ratio method .    _ the strong - coupling expansion _ ii + by explicitly identifying all distinct triangulations of volume @xmath206 , and calculating the corresponding measure term , we calculate the first 12 terms in the sce for an arbitrary value of @xmath37 .    [ tab.gamma ]    .estimates of the exponent @xmath19 , for @xmath202 , @xmath207 , @xmath208 , and @xmath209 ; ( a ) from a fit to @xmath199 , eq .  ( [ eq.mu_pseudo ] ) , ( with [ @xmath210 ) ; ( b ) from the minbu distribution , eq . ( [ babymax ] ) , both including and excluding the contact term , @xmath161 and @xmath205 ; and ( c ) from the first 17 terms in the sce . [ cols="^ , < , < , < , < , < " , ] as discussed in section  5.2 , for minbu s on maximal triangulations , the interactions at neck  the _ contact _ term  play a ( unusually ) big and non - trivial role in the crinkled phase . it is essential to take these interactions into account when estimating @xmath19 from a fit to the measured minbu distributions . in this appendix we analyze the behavior of the contact term and explain why it plays such an important role for this particular ensemble .    in the nomenclature of section  3 we can write the minbu distribution , eq . ( [ babymax ] ) , for @xmath197 and measured on a stacked sphere of size @xmath103 , as @xmath211 however , for @xmath212 the situation is different . we must take into account the different possible distributions of vertex orders , @xmath213 and @xmath214 , on the marked simplex ( the neck ) , both on the baby and on the mother . the distribution is then @xmath215 defining @xmath216 we can rewrite eq . ( [ bab ] ) as @xmath217     & \equiv   & { \cal b}_n(b;\beta)\cdot { \cal c}_n(b;\beta)\;,\end{aligned}\ ] ] where we define @xmath218 as the _ contact _ term ..    the following argument gives some idea of the relevance of the contact term in the different phases . if we assume that the distribution @xmath219 factorizes , @xmath220 then the distribution @xmath221 is related to the vertex order distribution @xmath190 , i.e.  the probability of finding a vertex with given order @xmath222 , through @xmath223 in the elongated phase the distribution @xmath221 decays exponentially , whereas in the crinkled phase vertices of large order are favored and the distribution has a power - law decay ( figure  [ fig.vord ] ) . an exponential decay of @xmath221 effectively introduces a cut - off in the sum over vertex orders in eq . ( [ ct ] ) . hence , for sufficiently large triangulations the contact term becomes independent of the volume . in contrast , in the crinkled phase there is no such cut - off and we can expect a strong volume dependence . this is demonstrated by numerical measurements of the contact term , shown in figure  [ fig.cont ] . in the elongated phase , @xmath224 , the contact term @xmath225 is more or less constant for @xmath226 , whereas for @xmath46 the contact term @xmath227 depends strongly on @xmath149 . although this effect can be expected for any ensemble modified by the measure eq . ( [ part.measure ] ) , the particular construction of maximal triangulations , by a repeated application of move @xmath228 , enhances the role of the contact term . for every time move @xmath228 is applied one minimal neck is generated , hence for a maximal triangulation the number of minbu s is fixed . thus the minbu distributions , measured on @xmath53 , is more constrained than measured on a general ensemble where triangulations with fewer minbu s are favored as @xmath229 . how does all this affect the measurements of @xmath230 ? the standard method for extracting @xmath19 from a minbu distribution is to impose a lower cut - off @xmath231 on the minbu size included in the fit to eq . ( [ babymax ] ) . the cut - off is then increased until a stable estimate of @xmath232 is obtained @xcite . in the elongated phase , and in general for most models of simplicial gravity , this methods works fine . however , for stacked spheres this method fails in the crinkled phase . unless the contact term is eliminated , i.e.divided out of the minbu distribution in the measurement process , a stable plateau in @xmath232 is not observed . an example of this is shown in figure  [ fig.cont ] for @xmath46 .
Abstract:

Answer: in the weak - coupling limit , @xmath0 , the partition function of simplicial quantum gravity is dominated by an ensemble of triangulations with the ratio @xmath1 close to the upper kinematic limit . for a combinatorial triangulation of the @xmath2sphere this limit is @xmath3 . defining an ensemble of _ maximal _ triangulations , i.e.  triangulations that have the maximal possible number of vertices for a given volume , we investigate the properties of this ensemble in three dimensions using both monte carlo simulations and a strong - coupling expansion of the partition function , both for pure simplicial gravity and a with a suitable modified measure . for the latter we observe a continuous phase transition to a _ crinkled _ phase and we investigate the fractal properties of this phase . bi / th-98/38 + december 1998    _ g.  thorleifsson , p.  bialas and b.  petersson _ +    facultt fr physik , universitt bielefeld d-33615 , bielefeld , germany +


Instruction: Write an abstract for the following acadamic paper.

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

Paper: in a contemporary context , the newtonian dynamics of a large number of self - gravitating particles describes the dynamics of putative dark matter , and to first order the formation and evolution of galaxies and clusters . dark matter halos identified in cosmological simulations have similar ( and nearly ) universal profiles , independent of the cosmological parameters or the masses of the halos or the epoch of their identification . these profiles have essentially the same form as those fitting early type galaxies ( e.g. , navarro et . al . 1997 ; merritt et al . 2005 ) . yet there is no satisfactory theory predicting these preferred choices of dynamical equilibria . steady states are characterized by the absence of coherent collective modes , which decay due to differences between motions on neighboring particle trajectories and interactions between the coherent modes and individual particle motions . these effects represent generic mechanisms driving gravitational systems toward collisionaless equilibria , ending the initial phase of violent relaxation , with its singnificant potential fluctuations and accompanying changes in particle energies . the associated processes are often referred to as phase mixing and landau damping respectively ( binney & tremaine 2008 , hereafter referred to as bt ) and should be most efficient when nearby trajectories rapidly diverge and spread in phase space . however , attempts at relating these trajectory stability properties to the choice of @xmath0-body equilibrium  or extracting virtually any interesting information at all from them  are frustrated by the following phenomenon .    as @xmath1 , particles are essentially affected by the mean field potential corresponding to the smoothed out mass distribution ( indeed , it can be rigorously proven that the mean field dynamics should describe asymptotically the time evolution of @xmath0-body system trajectories as @xmath1 ; braunn & hepp 1977 ) . the @xmath0-body problem is then transformed to @xmath0 1-body problems in this self consistent field , and steady state systems with sufficient symmetry are near integrable  e.g. , spherical systems are completely integrable when in a steady state , supporting only regular orbits , which conserve as many isolating independent integrals of motion as degrees of freedom     the phase space distances between trajectories moving on regular orbits generally diverges linearly in time  ( bt ) . thus the distance between trajectories of systems with only regular orbits should also diverge linearly in time . this should imply that , for an @xmath0-body system with smoothed out potentials supporting only regular orbits , any exponential divergence timescale ( indicating irregular trajectories and chaotic dynamics ) should tend to infinity as @xmath0 is increased . nevertheless , paradoxically , it has long been realized that point particle @xmath0-body systems invariably reveal exponential divergences on a timescale that does not increase with @xmath0 , even for systems with smoothed background potentials that only support regular orbits ( miller 1964 ; goodman , heggie & hut 1993 ; kandrup & sideris 2001 ; hemsendorf & merritt 2002 ) . this miller instability however saturates at progressively smaller scales as @xmath0 increases , and trajectories moving within a distribution of fixed point particles have been found to behave more and more like smooth potential counterparts as the number of background particles increased ( valluri & merritt 2000 ; kandrup & sideris 2001 ; hut & heggie 2002 ) . this suggests that the miller instability is a direct consequence of representation in terms of point particles and that it is of little relevance if one is concerned with the gross structure and macroscopic evolution of large-@xmath0 gravitational systems . furthermore , when the singularity in the potential is removed ( e.g. when @xmath2 ) , the exponential divergence timescale between neighboring trajectories does increase with @xmath0 ( goodman , heggie & hut 1993 ) and correlate with the evolution of dynamical variables like angular momentum ( el - zant 2002 ) . it also worth noting that rigorous derivation of the approach toward the mean field limit in large-@xmath0 gravitational systems requires a smoothed two body potential ( more precisely , @xmath3 is required to have continuous and bounded first and second derivatives ; braun & hepp 1977 ) . this again suggests that , if progress is to be made in understanding the the evolution of macroscopic quantities ( e.g. , density and velocity dispersion profiles ) in terms of the dynamical stability of trajectories , one has to move beyond the anomalous case of discrete point particles . the purpose of the present study is to show that , when the singularity in the two - body potential is removed , physically interesting information , pertaining to the choice of equilibrium , can be extracted from the stability properties of @xmath0-body trajectories . this will be done by using a geometric method for characterizing stability ( section  2 ) , where the incompleteness of the pure point particle description of trajectory stability is clearest . when the singularity in the two body potential is smoothed the miller instability is removed . what transpire instead are relations between macroscopic quantities characterizing equilibrium configurations ( as described in section  3 ) . in section  4 , we propose an interpretation of the approach to equilibrium and associated damping of coherent modes in terms of the stability properties of a system s trajectories ; in section  5 this interpretation is tested and discussed . an explanation of the origin of equilibrium ( spatial and phase space ) density profiles of dark matter halos and elliptical galaxies follows from the context thus set . the divergence of nearby trajectories can be also be deduced from geometric representations of motion and stability in gravitational systems ( e.g. , gurzadyan & savvidy 1986 ; kandrup 1990 ; cerruti - sola & pettini 1995 ; el - zant 1997 ; el - zant 1998 ) . though the geometric description is not unique ( e.g. , casetti et . al . 2000 ) , one well known formulation has the advantage of involving a positive definite metric defined directly on the lagrangian configuration manifold . for free particles , this manifold is simply the euclidean space spanned by the @xmath4 particle coordinates . in general , a system will move on a curved subspace , with velocity vectors defining its tangent space ( the phase space being the cotangent bundle ; arnold 1989 ) . a locally euclidean structure is generated on this space by the quadratic ( in the coordinate derivatives ) form of the kinetic energy ( arnold 1989 ) . for a system conserving its total energy @xmath5 , the associated conformally flat metric can be expressed in terms of the cartesian coordinates and the system kinetic energy @xmath6 as ( e.g. , lanczos 1986 ) @xmath7 the curvatures associated with this metric determine the stability of its trajectories . at the simplest intuitive level , it is clear why curvature should be related to stability . for example , the distance between neighboring geodesics ( great circles ) on a sphere oscillates sinusoidally as one moves from pole to pole , while geodesics on a negatively curved saddle are hyperbolic and diverge away from each other . the case of a flat plane corresponds to marginally stability , with distance between straight lines increasing linearly along their lengths .    in the present case , the variation of the normal distance between neighboring geodesics ( and therefore dynamical trajectories ) defined by the above metric is determined by the two dimensional sectional curvatures @xmath8 in @xmath9-dimensional configuration space , where @xmath10 is the system velocity vector and @xmath11 define directions normal to it . and it can be rigorously shown that when all the @xmath8 are negative , trajectories are exponentially unstable  ( anosov 1967 ; arnold 1989 ) . isotropic @xmath0-body systems with singular two - body potentials are predicted to tend toward this state as @xmath0 is increased and direct two - body collisions are ignored  ( gurzadyan & savvidy 1986 ; kandrup 1990 ) . perplexing trends , similar to those obtained from direct calculations ( i.e. exponential divergence timescales that do not increase with @xmath0 ) , are then recovered when the full gravitational field is taken into account  ( kandrup 1990 ) . this is a direct consequence of representing the matter distribution in terms of discrete point particles . in this paper we wish to move beyond this anomalous case .    we will be interested in how spherical @xmath0-body systems with isotropic velocities damp out fluctuations , we therefore seek an averaged measure of their response to random perturbations . averaging the two dimensional curvatures over all possible directions @xmath12 and trajectories ( geodesics ) at a given point in configuration space one obtains the the ricci scalar ; which , for @xmath13 , can be expressed as @xmath14 where @xmath15 and @xmath16 are the gradient and laplacian , taken with respect to coordinates @xmath17 , and @xmath18 are their particle masses ( which will be assumed to be equal ) .    the complete equations for trajectory stability will involve the two dimensional curvatures , and are usually written in terms of the geodesic distance @xmath19 ( e.g. , kandrup 1990 for a clear derivation ) . a correspondingly transformned equation measuring average stability  involving the scalar curvature and describing stability in the time domain  is given by  ( e.g. , cerruti - sola & pettini 1995 ) @xmath20 here @xmath21 measures the mean divergence of trajectories and @xmath22 near dynamical equilibrium , @xmath23 and the dominant term on the right hand side will be the curvature term . for singular newtonian potentials , eq . ( [ scalar ] ) predicts negative @xmath24 ( since the first term vanishes ) , and exponential instability is present even in the case of spherical equilibria . however , because of the singularities present when any two particle trajectories cross , one is here dealing with an incomplete manifold ; it is not clear if straightforward application of the geometric approach remains appropriate in this case  ( abraham & marsden 1978 ) . the miller instability thus figures again as a direct consequence of discrete representation in terms of point particles . as it will now be show however , the predictions are qualitatively different , and much more interesting , when one moves beyond this particular case . when the singularity in the potential is softened , and as @xmath1 , the first term in eq . ( [ scalar ] ) represents ( through the poisson equation ) a smoothed out density field integrated over the particle distribution ( with the softening length acting as smoothing parameter ) . in dynamical equilibrium , the two terms on the right hand side of eq . ( [ scalar ] ) are then comparable , reflecting a kinetic pressure versus gravity balance . to see this , assume @xmath25 . in this case , one finds that @xmath26 where @xmath27 and @xmath28 refer to the density and acceleration at the position of particle @xmath29 . this reflects a balance between a pressure term ( a density multiplied by the velocities squared ) , and a gravitational binding term ; it is therefore expected to be at least approximately satisfied for systems in dynamical equilibrium . in fact , by dividing by @xmath0 and defining @xmath30 , @xmath31 , @xmath32 as the means square speed of particles and @xmath33 , we can rewrite the above equation as @xmath34 assuming that , in terms of a system characteristic length scale @xmath35 , we can write @xmath36 and @xmath37 , then eq . ( [ eq : equ ] ) simply implies that @xmath38 and @xmath39 , which is characteristic of systems in dynamical equilibrium . the results obtained thus far suggest that the two terms of eq  ( [ scalar ] ) should be of the same order in the case of systems in dynamical equilibrium . it is therefore reasonable to expect that such configurations will correspond to small values of @xmath40 . nevertheless , this does not imply that @xmath24 will exactly be zero for all such systems . in particular , as shown in appendix  a , @xmath24 will not vanish for systems dominated by nearly constant density isothermal cores , which are the natural product of the process of classic violent relaxation gone to completion ; ; lynden - bell 1967 ; nakamura 2000 ) . however , it has been known for sometime now that such distributions do not match those of elliptical galaxies or dark matter halos . instead , of all possible collisionless near - spherical equilibria , these structures exhibit a narrow range of density profiles , characterized by central cusps that can be reasonably described by single power laws and extended halos with steeper power laws . it is therefore natural in our context to ask whether the latter correspond to @xmath41 ; and whether , in general , this condition gives us any insight into choice of gravitational equilibriia resulting from violent relaxation . we tackle the second question first , leaving the first one to the next section . in this regard we recall that @xmath24 describes the stability characteristics of trajectories ( through eqs .  3 and 4 ) . in dynamical equilibrium trajectories are stable if @xmath24 is positive and diverge exponentially if @xmath24 is negative . before equilibrium is reached however the density and potential of a system started far from equilibrium will fluctuate . these fluctuations can lead to trajectory instability even if the average @xmath24 is positive . the associated divergence and phase space mixing should result in the damping of coherent modes ; and large scale fluctuations can only last up to and untill till these collective motions are damped . this means that the system should settle in a collisionless steady state when it finds a configuration where its response to fluctuations is least coherent ; and this should correspond to dynamical equilibria where the trajectory divergence and phase mixing in response to fluctuations are maximal . at the same time , it is required that , in the absence of fluctuations , trajectories of a spherical system should not be exponentially unstable in the infinite-@xmath0 limit , so as to recover an integrable system . this implies that @xmath24 can not be negative when the system is in exact equilibrium . ( [ eq : stab ] ) suggests that configurations that satisfy these requirements have @xmath41 at equilibrium . this is easy to understand qualitatively : if average @xmath24 is positive and large , fluctuations ( in @xmath24 and @xmath42 ) have a smaller effect , and solutions of ( [ eq : stab ] ) tend to those of a harmonic oscillator ; conversely , if the average @xmath24 is near zero , @xmath43 is dominated by fluctuations , which can lead to significant instability , especially that in that case fluctuations in @xmath43 can lead to negative values . a system starting far from equilibrium , with initial conditions that allow for violent relaxation , should then end up in such marginally stable states . this intuitive argument is investigated more detail in appendix  b for the case of periodic oscillations and stochastic fluctuations . in the next section numerical simulations will be presented that suggest that systems starting further from equilibrium do indeed end up with @xmath24 closer to zero . but before closing this section it is important to note that though the above interpretation stands in seeming contradiction with that proposed by kandrup and collaborators , whereby it was suggested that trajectory divergence causes relaxation rather than subduing it , the clash here is only partial . the aforementioned authors reached their conclusions while studying trajectories moving through either stationary non - integrable potentials ( e.g. kandrup & siopis 2003 ) , or in particular spherical potentials with time dependent mass or length scales ( but which were otherwise rigid and fixed ; e.g kandrup et . al . 2003 ) . in the latter case they found that the time dependence can be a significant source of trajectory divergence , which may drive phase space mixing associated with violent relaxation . up to this point there is no contradiction with the interpretation outlined above ; fluctuations in the density and potential will indeed lead to trajectory instability and phase space mixing . however , when one takes into account the evolution of the self consistent field , it will be realized that the trajectory divergence and phase space mixing should also lead to the decay of the very collective modes that are causing the time dependence ( and associated trajectory instability , mixing and relaxation ) in the first place . when these are efficiently damped , no further evolution can occur ( in the collisionless limit ) . this is in line with our contention that the final product of violent relaxation do not correspond to classic most probable configurations but to ( quite different ) ones , which are most efficient in damping out collective modes via trajectory instability and associated phase space mixing . we now examine the above interpretation by means of direct @xmath0-body simulations . since we are interested in estimating the densities and accelerations in the continuum limit , we use a technique  ( hernquist & ostriker 1992 ) , which expands the density and potential in smooth functional series . in this study , we are interested in strictly spherical configurations and the results described below only make use of the radial expansions in this scheme , which is carried up to order @xmath44 ( however it has been verified that the results are not significantly altered if azimuthal expansion terms are included , as noted in the discussion below ) . systems are sampled using a @xmath45 particles ( except for one case with 500 000 as also described below ) , and started from homogeneous spatial initial conditions inside a unit sphere of unit mass . this configuration has a natural timescale @xmath46 , which will be our time unit . we have run configurations starting with isotropic velocities that are either constant or following various decreasing or increasing functions of radius . the results shown here correspond to the latter case , for the following reason . unless the initial velocities are vanishingly small , uniform initial velocities will imply that the initial pseudo - phase space density @xmath47 , does not diverge anywhere . on the other hand , collisionless evolution implies that @xmath48 generally decreases along particle trajectories  ( bt ) , so that the maximal final @xmath48 can not exceed the initial one . one of our aims is to explain the dynamical origin of cosmological halos , which have ( a centrally divergent ) @xmath49  ( taylor & navarro 2001 ) , which in turn requires an intial phase space density that is also divergent ( so that this initial density can , in principle , be mapped via collisionless dynamics into that characterizing dark halos . from among various functional forms satisfying this condition that were tried , results will be shown for initial @xmath50 , noting that the trends reported seem generic . we vary the initial virial ratio @xmath51 from @xmath52 , corresponding to near equilibrium , to @xmath53 , corresponding to a system started far from equilibrium . given the sampling errors for finite @xmath0 , @xmath24 will never tend to zero exactly in practice ; and the absolute values of the densities and accelerations ( and local dynamical times ) strongly depend on the central concentrations of the final configurations . a relative measure is thus required for comparing these . we use the normalized quantity @xmath54 which measures the departure from equality in eq . the results are in fig .  1 , which shows that the equilibrium @xmath55 is indeed always of order unity or less , and that the further from equilibrium the initial state , the closer the evolved system is to a configuration that minimizes @xmath56 , in line with the general interpretation presented in the previous section . the initial @xmath55 is zero when @xmath57 ; systems with initial @xmath58 start with @xmath59 and then explore state space until they achieve configurations where ( according to eq . [ eq : stab])oscillations are efficiently damped . we note that a similar trend was found when systems were started with constant initial phase space density @xmath48 , except that generally larger final values of @xmath55 are obtained . the origin of the density profiles of bound structures identified in the context of cosmological simulations have been intensely investigated ; the mechanism for their advent was sought for example in terms of the merging histories of smaller structures in the way of forming larger halos ( e.g. , syer & white 1998 ; dekel et . al . 2003 ) . however , the same sort of structures tend to materialize even in simulations where a cutoff has been applied to the power spectrum and the collapse is monolithic ( e.g. , huss et . ; moore et . 1999 ; wang & white 2009 ) . these profiles thus seem to arise nearly independently of the condition of their formation ; merging may conserve their form ( e.g. , boylan - kolchin & ma 2004 ; el - zant 2008 ) , but is not necessary to its genesis . it is thus natural to inquire if a simple explanation may be sought in terms of the stability issues we have addressing in this paper . we know from the previous subsection that systems started far from equilibrium to undergo simple monolythic collapse have near - zero final average curvatures ; they are therefore marginally stable according to the interpretation outlined in section  4 . if , as it is suggsted , violent relaxation should generically terminate in such states , and if violent relaxation is predominatly responsible for the structure of dark halos , then those halos should correspond to configurations with near - zero curvature . this would be the case if they corresponded to the final configurations of our simple collapse simulations  despite the lack of cosmological context or associated merging processes in these simulations . indeed , one finds , as shown in fig .  2 , that the equilibrium density and phase space density profiles of systems started far from equilibrium and undergoing simple monolithic collapse , starting from cold initial conditions , are remarkably similar to those of cosmological halos ( and also shallow cusp elliptical galaxies ; merritt et al . 2005 ) , despite the lack of cosmological context . systems started from virial equilibrium on the other hand tended to end a up with large central cores ( as do systems started with initial phase space density profiles that disallow the final configurations shown ) . this simulation was repeated with several values of radial and azimuthal expansion orders in the herquist - ostriker scheme , including one simulation with @xmath60 particles and 20 radial expansion terms and 12 azimuthal ones , without significant change in the results . there are some difference between the configuration shown and idealized halo profiles . our outer density profile is steeper than the @xmath61 form of cosmological halos ; however finite configurations can not have @xmath61 outer profiles . ( cosmological halos are not isolated equilibrium structures as their outer regions are subject to mass infall , which reduces the steeness of the outer profile ; e.g. , gott 1975 ) . the functional form of our empirical spatial density fit is also slightly different from standard nfw ( where @xmath62 ) , and there seems to be some flattening toward the center of the @xmath63 cusp . yet , similar departures from ideal nfw often also apply to cosmological halos ( navarro et . al . 2004 ) . our argument for the production of the universal profiles  when the initial phase space density allows for their advent  is that the associated configurations are most efficient in damping out large scale coherent oscillations . we now explicitly verify that among the configurations we have simulated ( which always have positive or near zero average curvature once they reach a quasi - steady state ) , those with smaller @xmath55 do efficiently damp out fluctuations . this is already apparent in fig .  1 , where fluctuations are quickly damped for such systems . this may be intuitively expected , for these systems lack homogeneous harmonic cores with nearly constant orbital frequencies and coherent response , and similar trends were indeed observed previously ( e.g. , david & theuns 1989 ) . indeed ( 3 ) allows in principle for the existence of long lived coherent modes when the average @xmath24 is large ( as discussed in appendix  b ) . further verification is obtained by perturbing the final states and measuring variations in the density distributions as systems , so perturbed , are evolved . the results of one such experiment , where the spatial coordinates were decreased by @xmath64 , are shown in fig . , the system that was perturbed from a state with smaller @xmath56 is seen to quickly settle back to an equilibrium that is closer to the unperturbed one , and to gradually close in further on it ( as the outer regions , with large crossing times , evolve .    from the above , one can conjure that the observed robustness of the universal profiles when perturbed may also be a reason behind their conservation under processes involving significant potential fluctuations  such as mergers , accretions and fly bys  that are important in the context of cosmological simulations . finally we note that , despite the formal similarity between profiles of elliptical galaxies and dark halos , the above interpretation strictly speaking applies only to the advent of dark halo configurations . the case of early type galaxies is much more complicated . the currently favored theory for the formation of early type galaxies involves a two - phase process : a highly dissipative initial stage accompanied by intense star formation , followed by a succession of dissipationless mergers ( e.g. , oser et . al . 2010 ) . the simple dissipationless collapse simulations presented here are evidently not of direct relevance to this rather involved two - phase model . however , though the initial phase is highly dissipative , it sets the initial conditions for largely dissipationless dynamics , presumably dominated by significant fluctuations characterizing an evolving time dependent system ( this should be especially the case if star formation is sufficiently efficient and the associated burst sufficiently intense ; indeed , oser et . al . suggest this phase bears an uncanny resemblance to claasic monolythic collapse models ) . if the general interpretation ( of section  4 ) is correct , marginally stable systems should in principle be generic products of such a situation . furthermore , as noted above , resulting profiles should be conserved in the second phase ( of dissapationless merging ) . this seems plausible ; however it is evidently beyond the scope of this paper to explicitly test the interpretation presented here in the complex context of the two - phase theory of the formation of early type galaxies . our main concern in this paper has been with the way mean field fluctuations are damped as a result of the mixing that occurs due to the instability of phase space trajectories . historically , however , it has been ( controversially ) argued that when trajectory instability stems not from large scale fluctuations but from discreteness noise , a new mode of relaxation , different from standard two body relaxation and with a correspondingly different timescale , arises ( gurzadyan & savvidy 1986 ) . as noted in the introduction and in section  2 , when the full self consistent gravitational field is taken into account , and with a singular potential , the prediction is that system trajectories will always be unstable on the timescale of a dynamical time or less ( the miller instability ) . gurzadyan & savvidy on the other hand employed a holtsmark distribution ( e.g. , chandrasekhar & von neumann 1942 ) , where the force function is not self consistent but is dominated by nearest neighbor interactions ( and thus varies as  @xmath65 ) . they thence deduced that the instability timescale scales varies as @xmath66 , and identified this with a new collective relaxation time . three things are to be noted in our present context . first , assuming that that fluctuations due to discreteness noise ( in the curvature and associated quantities , as density and force ) scale as @xmath67 , then eqs .  ( [ eq : stab ] ) and ( a3 ) predict that , when the curvature is much larger than discreteness noise fluctuations , unstable trajectories would diverge on a timescale @xmath68 ; while @xmath69 when the reverse is true . thus , for statistically significant @xmath0 and unless the curvature is extremely closely fine tuned to zero , the divergence timescale for nearby trajectories will scale more as the standard two body relaxation time than the gurzadyan & savvidy timescale . secondly , even with an @xmath70 scaling , the source of the instability in our case are fluctuations around a near - zero curvature , calculated taking into account the full self consistent field , and not a small negative curvature arising from closest neighbor forces , as in the calculation of gurzadyan & savvidy . the third point concerns interpretation . relaxation on the two body timescale , as usually understood in stellar dynamics , is accompanied by structural changes in the collisionless equilibrium configuration of the system , due to diffusive random walks undergone by the integrals of motion associated with the mean field potential ; e.g. , energy relaxation leads to core contraction and collapse . whether divergence between neighboring trajectories directly leads to such effects on the corresponding exponential timescale is not clear , since such divergence will include changes in phases between particle trajectories as well vartiations in the integrals of motion associated with the corresponding trajectories in the mean field potential . much further detailed investigation would be needed to establish a precise relation ( if any ) with relaxation as normally understood in stellar dynamics .    on the other hand , it seems clear that the efficiency of phase space mixing due to exponential divergence should directly correlate with the decay of coherent collective modes , as even simple phase mixing ( without any exponential instability or changes in the integrals of motion ) is effective in damping of these modes . thus the interpretation of the effect of trajectory divergence in the context of the present paper appears to rest on firmer ground . an analogous interpretation in the case of discreteness noise is that the timescale of decoherence of collective modes should depend on the number of particles in the system ( as observed , e.g. , by david & theuns 1989 ) . it is hoped that the present exploratory study has shown that the methods employed , which have hitherto been of little use , are potentially powerful tools in understanding the structure of gravitational systems ; and that it would stimulate further work on these fundamental  yet lately little researched  problems pertaining to the evolution of newtonian gravitational systems . according to the classic idea of violent relaxation , a collisionless system starting far from dynamical equilibrium reaches such a state via collective oscillations . if the associated mean field fluctuations persist until the bulk of the system is described by a most probable phase space distribution , the resulting distribution functions are of maxwellian form , which are quite different from the distribution functions of dark matter halos for example ( e.g. , widrow 2000 ) . the associated density profiles are dominated by isothermal cores ( lynden - bell 1967 ) , which are not observed in elliptical galaxies or in simulated cosmological halos . in contrast , in the scenario proposed here , violent relaxation should persist only until a system , starting far from equilibrium , finds a configuration where collective oscillations are efficiently damped ; and such configurations support dynamical trajectories that are marginally stable . this condition necessarily leads to equilibria that are different from isothermal ( except for exact singular isothermal pheres ) , and which are remarkably similar to those observed and modeled , provided that the initial conditions do not impose phase space constraints that disallow evolution to such states . these systems do not ring and are more robust ; in the sense that , when perturbed , they rapidly return near their initial states , as opposed to configurations dominated by harmonic cores , where small perturbations give rise to significant structural changes accompanied by long - lived collective modes . there is an intuive explanation of this result ; as opposed to configurations dominated by nearly constant density harmonic cores , these systems display relatively large orbital frequency ranges , and are thus more efficient in damping collective modes . it is also possible that this is a reason behind the conservation of the universal profiles through major mergers and other processes involving significant potential fluctuations that are important in a cosmological context , though this contention has not been investigated explicitly here . i would like to thank the referee for comments that lead to significant improvement in substance and presentation , replacing @xmath28 and @xmath27 in eq . ( [ eq : press ] ) by their continuum limit counterparts , and integrating over the spherical mass distribution , one gets @xmath71 where @xmath72 is the radial density . defining the local acceleration @xmath73 ( with @xmath74 the mass enclosed within radius @xmath75 ) , one can write @xmath76 in dynamical equilibrium one can use the jeans equation for an isotropic system , @xmath77 ( e.g. , bt ) , to get @xmath78 where @xmath79 is the one - dimensional velocity dispersion . there is one easily identifiable solution of eq . ( [ curcrude ] ) ; it corresponds to @xmath80 and @xmath81 . this is a singular isothermal sphere ; however its gravitational force , and both terms in eq . ( [ curcrude ] ) , diverge as @xmath82 . well behaved systems  including those with inner power law density profiles shallower than @xmath83 , as is the case of simulated dark matter halos  have @xmath84 as @xmath82 and @xmath85 ( e.g. , when an outer density profile is integrating by parts , noting this , and that the kinetic energy is related to the average velocity dispersion by @xmath86 , and @xmath87 , one gets @xmath88 if effective pressure and temperature functions can be assigned , then @xmath89 and @xmath90 , and @xmath91 save for the integrations over the mass distribution , this form is that of an ideal gas equation of state . nevertheless , because of the factor @xmath92 , there are no associated isothermal states .    thus nearly isothermal systems ( with non - diverging @xmath93 ) can not be marginally stable in the sense of @xmath41 . therefore they should not result from violent relaxation according to the scenario proposed in this paper , even though thay are predicted via entropy maximization in the context of classic violent relaxation theory . ( [ eq : equi ] ) also shows that this is the case with any configuration dominated by a nearly constant density core . this is in line weith the numerical results of section  5 . eq .  ( [ eq : stab ] ) can be written in terms of the average @xmath43 and time varying term @xmath94 as @xmath95 in general , phase space averages are not well defined for gravitational systems ( which have a non - compact phase space ) . nevertheless , every exact ( @xmath1 ) collisionless equilibrium has a characteristic @xmath24 and associated @xmath96 . a system oscillating about such an equilibrium will have an average @xmath43 close to this value ( as the second and third terms in eq . 4 are zero in exact equilibrium , their contributions to average @xmath43 being second order in the fluctuations ) , with rms fluctuations @xmath97 . if the time variations are periodic this is a hill equation , with its well known instability strip increasing in measure as @xmath98 increases relative to @xmath99 , even for positive @xmath99 ( e.g. magnus & winkler 2004 ) . solutions will always be unstable if the curvature can reach negative values . on the other hand , for positive @xmath99 that is large relative to the time varying term , the instability regions are vanishingly small , since n that case , solutions eq . ( [ eq : av ] ) tend to those of a harmonic oscillator . these results suggest that for a given level of periodic fluctuations instability is more likely to occur ( and it occurs with smaller divergence timescale ) , as the average curvature @xmath99 decreases .    the existence of stable solutions for large @xmath99 may account for the apparently locked long lived oscillatory modes observed in systems starting near virial equilibrium ( fig .  1 and fig  3 ) . as pointed out in the discussion of fig .  3 , there is also an intuitive explanation of these oscillations : the associated configurations exhibit dominant harmonic cores , where orbital frequencies change little , facilitating the persistence of coherent modes . this result is in line with those of the previous appendix , showing that systems dominated by nearly constant density cores are unlikely to correspond to near zero curvature .    if the fluctuations in @xmath43 can be considered as stochastic ( more precisely as normally distributed random variables ) , one can use a standard formulation for solving stochastic equations , developed by van kampen ( 1976 ) , to show that solutions are always unstable , even for positive @xmath99 , with the timescale of exponential divergence between nearby trajectories decreasing as the rms fluctuations @xmath98 increase relative to @xmath99 . to see this , assume that the fluctuations @xmath98 have characteristic timescale @xmath100 . this gives rise to the following estimate of the exponential divergence timescale of the quantity @xmath21 ( e.g. , eqs . 79 of  casetti et . al . 2000 ) @xmath101 with @xmath102 . the natural timescale for fluctuations in a gravitational system is the dynamical ( or crossing ) time , so @xmath103 . and since @xmath43 comes in units of inverse time squared , we can write @xmath104 and @xmath105 , and express the exponentiation time associated with trajectory divergence in terms of the natural timescale as @xmath106 this is proportional to @xmath107 for @xmath108 and tends to @xmath109 as @xmath110 . for a given fluctuation level @xmath111 , the decoherence of trajectories due to local divergence will be more efficient , i.e. its exponential time will be smaller , for smaller @xmath112 , and therefore @xmath99 and average @xmath24 . valluri , m. , & merritt , d. 2000 , in the chaotic universe , proceedings of the second icra network workshop , advanced series in astrophysics and cosmology , vol.10 , ed . v. g. gurzadyan & r. ruffini ( world scientific ) 229
Abstract:

Answer: the origin of equilibrium gravitational configurations is sought in terms of the stability of their trajectories , as described by the curvature of their lagrangian configuration manifold of particle positions  a context in which subtle spurious effects originating from the singularity in the two body potential become particularly clear . we focus on the case of spherical systems , which support only regular orbits in the collisionless limit , despite the persistence of local exponential instability of @xmath0-body trajectories in the anomalous case of discrete point particle representation even as @xmath1 . when the singularity in the potential is removed , this apparent contradiction disappears . in the absence of fluctuations , equilibrium configurations generally correspond to positive scalar curvature , and thus support stable trajectories . a null scalar curvature is associated with an effective , averaged , equation of state describing dynamically relaxed equilibria with marginally stable trajectories . the associated configurations are quite similar to those of observed elliptical galaxies and simulated cosmological halos , and are necessarily different from the systems dominated by isothermal cores , expected from entropy maximization in the context of the standard theory of violent relaxation . it is suggested that this is the case because a system starting far from equilibrium does not reach a most probable state via violent relaxation , but that this process comes to an end as the system finds and ( settles in ) a configuration where it can most efficiently wash out perturbations . we explicitly test this interpretation by means of direct simulations .