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: recent progress in the description of integrable partial difference equations is to a great part due to the consistency approach @xcite , in particular in the form of 3 dimensional consistency - around - a - cube ( cac ) . one of the highlights of this approach is the immediate existence of a lax pair and bcklund transforms ( bt ) , which can be directly constructed from the `` side - equations '' of the cube @xcite . one can say , in effect , that the set of side equations yields the bt and the lax pair . one situation where side equations have been used effectively is in constructing soliton solutions to the lattice equations @xcite . originally it was assumed that on all faces of the cube the equations were the same in form , depending on the relevant corner variables ( one component at each corner ) and spectral parameters . recently , 3d consistent sets have appeared , with different equations on the faces @xcite . in this more general context it is of interest to take a closer look at the lax pairs , bcklund transforms and consistency , and investigate what they are good for . it should be noted that in the context of partial differential equations the existence of trivial lax pairs is well known @xcite and similar examples have also been noted for some discrete equations ( see @xcite , chapt . we will first recall how consistency around the cube , existence of lax pairs and bcklund transformations are intimately related for lattice maps on the square lattice given by multi - affine relations ( sections [ sec-3d - reminder ] and [ sec - lax - bt ] ) . these considerations apply to the elementary cells , and are _ local_. in section [ sec - examples ] , we describe specific examples having the cac property . we detail in particular the already known explicit forms of the equations and lax pairs which we use in the rest of this paper . in particular , we find that in some cases the zero curvature condition ( zcc ) yields two different equations that can be used to define rational evolution in the lattice . we then address the _ global _ problem of defining equations over the whole lattice , with the guideline given by the 3 dimensional structure coming from the cac construction ( section [ secalge ] ) , and check integrability with the calculation of algebraic entropy . in section [ sec-2x2 ] , we push the use of @xmath0-matrix lax pairs to its limits , by constructing discrete systems over a larger sub - lattice of the original lattice . this brings to the light interesting ( integrable ) structures related to a generalized form of the functional yang - baxter equations  @xcite . the starting point is a regular 2d square lattice , with vertices labeled by integers @xmath1 . functions @xmath2 are associated to the vertices , and they are subject to a constraint at all elementary cells . this constraint is expressed by an equation @xmath3 , assumed to be multi - affine in the four vertex variables . it should depend on all 4 vertex variables , and it should not factorize . it may also depend on some parameters . sometimes the parameters can be associated to specific directions of the lattice , in which case they appear as `` spectral parameters '' . to ease the notation , one usually denotes the running value @xmath2 by @xmath4 , and for neighboring values one only indicates the shifts : @xmath5 , or in 3d setting , @xmath6 etc .    for multidimensional consistency one needs to build a cube on top of a square and give equations on all six faces of the cube , see figure [ fcac ] , the bottom equation being the original one . ( 3482,2813)(0,-10 ) ( 450,1883 ) ( -150,1883)(0,0)[lb]@xmath7 ( 1275,2708 ) ( 800,2708)(0,0)[lb]@xmath8 ( 3075,2708 ) ( 3375,2633)(0,0)[lb]@xmath9 ( 2250,83 ) ( 2550,8)(0,0)[lb]@xmath10 ( 1275,908 ) ( 925,908)(0,0)[lb]@xmath11 ( 2250,1883 ) ( 1750,2008)(0,0)[lb]@xmath12 ( 450,83 ) ( -00,8)(0,0)[lb]@xmath13 ( 3075,908 ) ( 3300,833)(0,0)[lb]@xmath14 ( 1275,2708)(3075,2708 ) ( 1275,2708)(450,1883 ) ( 450,1883)(450,83 ) ( 3075,2708)(2250,1883 ) ( 450,1883)(2250,1883 ) ( 3075,2633)(3075,908 ) ( 1275,908)(450,83 ) ( 1275,908)(3075,908 ) ( 2250,1883)(2250,83 ) ( 450,83)(2250,83 ) ( 3075,908)(2250,83 ) ( 1275,2633)(1275,908 )    supposing @xmath15 are given , one can compute @xmath10 using the bottom equation , @xmath7 using the left side equation , and @xmath9 using the back side equation . we may then calculate the final value @xmath12 in three different ways , using the front- , right- and top - equations , respectively . the equations are said to be consistent - around - a - cube ( cac ) if the three ways yield the same value for @xmath12 . note that it is also possible to check cac with other initial values that allow full evolution , for example , @xmath16 . the lax / bt approaches differ from cac in that only the side equations are used as input , while the top and bottom equations are supposed to be derived from them . the lax pair for the bottom equation is constructed from the side equations by isolating the variables @xmath18 , and writing the equations as :    [ side - eqs ] @xmath19    the coefficient functions @xmath20 will be affine linear in their arguments . we now introduce the homogeneous coordinates @xmath21 for @xmath22 and its shifts by @xmath23 this amounts to considering @xmath22 as belonging to the projective space @xmath24 . we denote by @xmath25 the pair @xmath26 which is defined up to a global factor . equations can then be written as @xmath27 where @xmath28 and similarly for the bar - quantities . since we are working in @xmath29 all equalities are projective equalities , and as a reminder of this we use the symbol @xmath17 to indicate that two matrices are equivalent if their entries are proportional , in other words , @xmath30 and @xmath31 belong to @xmath32 . the lax matrices provide parallel transport of @xmath25 along the bonds of the lattice . the zero curvature condition means that the parallel transport along any closed path on the lattice is trivial . it is necessary and sufficient to ensure that taking @xmath25 from position @xmath33 to @xmath34 via the two routes @xmath35 and @xmath36 gives the same result : @xmath37 this matrix relation yields three scalar equations . these equations are written in terms of the variables @xmath38 , and in order to satisfy the zcc we must impose a constraint on these variables . in standard cases , this constraint will just be the bottom equation . indeed , in integrable cases the three equations in have a common factor . it may also happen that the zcc is satisfied automatically , if the side equations are simple enough , or sometimes the common factor may factorize , yielding multiple choices . we will examine specific examples of this phenomenon below . a bcklund transformation ( bt ) is generated by a set of equations @xmath39 on two sets of functions @xmath40 and @xmath41 , such that eliminating one set of functions say @xmath40 ( resp . @xmath41 ) one gets an equations @xmath42 ( resp . @xmath43 ) on the other set ( for details , see proposition 4.1 in @xcite . )    for a 3d consistent set of equations the side equations provide a bt between the top and bottom equations . in order to derive the bottom equation from the side equations one proceeds as follows :    1 . solve @xmath7 from the left equation . 2 .   solve @xmath9 from the back equation . 3 .   solve @xmath12 from the front equation . 4 .   after this the right equation can be written as @xmath44 , where @xmath45 are polynomials in @xmath38 . the greatest common divisor ( gcd ) of @xmath46 ( or one of its factors if it factorizes ) will be the bottom equation . one can find the following explicit forms for the @xmath46 : @xmath47    that the lax and the bt approaches yield the same equations can be seen as follows : let us denote @xmath48 , then the matrix elements of @xmath49 are related to @xmath46 as follows ( here subscripts indicate the matrix element ) : @xmath50 thus @xmath51 iff @xmath52 . of course one can equally well use the side equations to solve for the variables @xmath53 from the left- , back- , and front - equations , after which the right - equation will be a polynomial in @xmath54 with coefficient depending on @xmath18 , with their gcd yielding the top equation . as a first example we consider the case where all the side equations are linear , i.e. , the left side equation is : @xmath55 and we have the same equation with suitable subscript changes on the other vertical sides . one then finds that neither the lax nor the bt approach yields anything , the conditions are satisfied without reference to the bottom equation . what about cac ? it will involve both the bottom and the top equations , so it should give some conditions . indeed , if one tries cac with and a completely general multi - affine bottom and top equations , related by a shift ( but with same parameters ) one finds a consistent set with the bottom equation @xmath56 this is a combination of the lattice modified kdv equation ( lmkdv ) ( aka q1 in the abs list @xcite ) , and the linear equation . a similar analysis can be done starting from the linearizable side equation @xmath57 the lax matrix for this system is diagonal and the zcc is automatically satisfied . from cac analysis we find that is compatible with the six parameter family of homogeneous bottom / top equations of degree @xmath58 @xmath59    we will examine the integrability of these equations in section ( [ secalge - lin ] ) . [ [ h1 ] ] h1 ~~    the lattice potential kdv ( lpkdv ) , which describes the permutability property of continuous kdv , is a paradigm of integrable lattice equations ( aka h1 in the abs list @xcite ) . for this system everything works well . [ we will return to this model in a new context in section  [ h1_2x2 ] . ] the model is given by @xmath60 after imposing this equation on the sides of the cac cube , with suitable variable changes , i.e. ,    [ side - h1 ] @xmath61    one easily finds    @xmath62\cdot { \rm h1},\\ \mathcal r_0&=&[x_1x_2(x_{12}-x)+ ( \beta-\gamma)x_1 + ( \alpha-\gamma)x_2]\cdot { \rm h1},\end{aligned}\ ] ]    with h1 , as given in , as the gcd . similarly by working on the bottom variables @xmath38 one obtains the 3-shifted h1 as the top equation . the lax matrices in this case are @xmath63 and one easily finds that @xmath64      this is an asymmetric deformation of @xcite    [ side - h1e ] @xmath67 the bt or lax computations give as gcd the bottom equation @xmath68 and similarly for the top equation we get @xmath69    if we draw a line connecting the corners appearing in the deformation term as a product , the lines form a tetrahedron inside the cube , see figure [ fcach1 ] . ( 3482,2513)(0,200 ) ( 450,2183 ) ( -150,2183)(0,0)[lb]@xmath7    ( 1275,2708 ) ( 800,2708)(0,0)[lb]@xmath8 ( 3075,2708 ) ( 3375,2633)(0,0)[lb]@xmath9 ( 2250,383 ) ( 2550,308)(0,0)[lb]@xmath10 ( 1275,908 ) ( 825,908)(0,0)[lb]@xmath11 ( 2250,2183 ) ( 1650,2258)(0,0)[lb]@xmath12    ( 450,383 ) ( -00,308)(0,0)[lb]@xmath13    ( 3075,908 ) ( 3300,833)(0,0)[lb]@xmath14    ( 1275,2708)(3075,2708 ) ( 1275,2708)(450,2183 )    ( 450,2183)(450,383 ) ( 450,2183)(2250,2183 ) ( 450,383)(2250,383 )    ( 3075,2708)(2250,2183 ) ( 3075,2633)(3075,908 ) ( 1275,908)(450,383 ) ( 1275,908)(3075,908 ) ( 2250,2183)(2250,383 ) ( 3075,908)(2250,383 ) ( 1275,2633)(1275,908 ) ( 1275,908)(450,2183)(2150,383)(2975,2708)(2925,2708)(1125,908)(2875,2708)(250,2183)(1025,908)(2000,383)(50,2183)(1850,383 )    the special thing about these equations is that the parallel sides are _ not _ obtained by shifts , but the shift must be accompanied by a @xmath70 rotation , as described in figure [ fcach1 ] . the fact that the parallel sides are not identical implies that the lax matrices associated to parallel bonds come in two different forms . letting @xmath71 and @xmath72 we have @xmath73 since we have two different parallel sides we can have two different bottom equations , one derived from @xmath74 as above , the other possible choice being @xmath75 the totality of equations on this cube , obtained from the previous one with @xmath76 etc actually corresponds to an inversion of the cube by @xmath77 which yields    [ side - h1e2 ] @xmath78    these are also consistent . we thus have two kinds of cubes , one defined with , and one with . these cubes can be glued together : note that the 2-shift of is equal to . in fact under the central reflection the diagonal lines in figure [ fcach1 ] would change to the other diagonals . we will return later to the question of gluing the cubes to fill the space . this model was proposed in @xcite ( see sec . 3.1 , case @xmath79 ) . in that paper the model was given in a cube with flipped coordinates @xmath80 , @xmath81 . in the coordinates of figure [ fcac ] the side equations of this model are given by    [ side - h1w ] @xmath82    [ [ bt - derivation - of - bottom - and - top - equations . ] ] bt derivation of bottom and top - equations . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    now computing the values of @xmath83 from left- , back- and right - equations , respectively , gives for the front equation an expression that does not contain @xmath22 at all . this expression factorizes into _ two _ factors and thus we could have two different bottom - equations :    [ h1w - bot ] @xmath84    similarly , working with the values at the bottom square we get two candidates for the top equations :    [ h1w - top ] @xmath85    notice that in down2 and top2 , there is an explicit dependence on @xmath86 even though these are 2d equations . in fact , one should now consider @xmath86 as a global parameter , although it was associated with the third dimension in the above derivation .    performing the usual cac computations with the given sides and different top and bottom equations reveals that the set of equations is consistent in two cases : with the pair ( bottom1,top1 ) , or with ( bottom2,top2 ) , which was given in @xcite . in other words , given the side equations there are two consistent ways of completing the cube . continuing further , we can also interpret the cube to provide a bt between the left and right equations . indeed , if we do the above bt construction on the corners of the left side equation the result is identically zero with ( bottom1,top1 ) , while ( bottom2,top2 ) produces the right - equation .    [ [ the - lax - matrices . ] ] the lax matrices . + + + + + + + + + + + + + + + + +    the standard procedure gives @xmath87 since the matrices are upper triangular the zcc implies @xmath88 remarkably enough , the above sum factorize as @xmath89.\ ] ] note that we can write also in the form @xmath90 where @xmath91 is a shift in @xmath92 and @xmath93 a shift in @xmath94 . this of course is in the form of a conservation law . so far we have only considered a single cube and its cac / lax / bt . but as the name indicates , lattice equations should be defined over the whole lattice . this brings further complications , for example with one cube we could freely do different mbius transformations in each corner of a cube , but when the cube is part of lattice such seemingly innocuous actions will affect neighboring cubes as well and can destroy the lattice structure . the rule is simple : cover the two dimensional lattice with consistent cubes , with the condition that adjacent vertical faces coincide exactly , that is to say their four corner values satisfy the same equation . this is expected to produce integrable lattice equations . we will follow this guideline for various models mentioned before , and systematically check integrability of the lattice equations so obtained , by calculating their algebraic entropy  @xcite : the vanishing of the entropy is a yes / no test which gives a clear cut separation between integrable and non integrable cases . we may briefly recall how to calculate the entropy . the local equation determines an evolution , starting from initial conditions given for example on a diagonal staircase ( lattice points of coordinates @xmath95 with @xmath96 or @xmath97 ) . the solution is then calculated on diagonals moving away from the diagonal of initial conditions , explicitly in terms of these initial conditions . the algebraic entropy is defined as the rate of growth of the degrees on these diagonals . exponential growth is generic , and polynomial growth is characteristic of integrability , while linear growth is associated with linearizable equations  @xcite .    the exact shape of the diagonal line on which the initial values are given is not important . if one modifies this shape locally , the sequence of degrees will change , but not its asymptotic rate of growth . this should be kept in mind for some of the models studied below . for example one could very well change the initial diagonal ( steps of height 1 and width 1 ) to a diagonal with bigger steps , for example with height 2 and width 2 , as this will not affect the calculation of the entropy . we can check the integrability of the two quad equations given in section ( [ linear - sides ] ) . the first one , i.e. , equation ( [ q1-lin ] ) , leads to the sequence of degrees : @xmath98 that is to say @xmath99 . this quadratic growth confirms integrability .    * remark * : since it is integrable , it should also be consistent with a nontrivial set of side - equations , such that one can produce it via lax / bt computations , but we will leave this open .    for the second equation of section ( [ linear - sides ] ) , that is to say the general homogeneous relation ( [ quad - lin ] ) of degree @xmath58 , the result is different . we get the following sequence of degrees : @xmath100 which may be fitted with the generating function @xmath101 indicating a non vanishing entropy @xmath102 , showing non - integrability . thus even cac is not sharp in this case . it was actually shown in  @xcite that the simple additional condition @xmath103 , renders ( [ lin - side2 ] ) integrable . we have already noted that in the h1@xmath104 model there are two different cubes related by inversion , and that these cubes can be glued together in a unique way , providing a black - white lattice . this problem has been discussed in detail in @xcite .    since the gluing process is completely fixed and each cube has the cac property it is expected that the composite lattice is integrable . we have calculated the sequence of degrees @xmath105 for the evolution defined by this model , with initial conditions on a diagonal as prescribed above . the outcome is the sequence @xmath106 fitted by the generating function @xmath107 the above sequence has quadratic growth . thus the entropy is vanishing and the model is integrable .      in the flipped h1@xmath104 case ( sec . [ fh1e ] ) , the side equations are the same , allowing simple gluing together of the cubes , but , as we saw in section ( [ fh1e ] ) , the side - equations are somewhat weak and allow two different compatible pairs of bottom / top equations . we may then cover the two dimensional lattice with consistent cubes , assigning either equation bottom1 ( top1 ) , which we call _ white _ , or equation bottom2(top2 ) which we call _ black _ to each elementary cell . this can be done in an arbitrary way if one just insists on having a compatible 3d structure of cubes over the 2d lattice . it is then natural to ask which of the configurations obtained in this way are integrable . let us consider periodic distributions . the lattice is divided into rectangular groups of cells of width @xmath108 and height @xmath109 . within such a rectangle , a fixed assignment is made , and the pattern is repeated periodically in both directions . ( a pattern with @xmath110 and @xmath111 gives a uni - colored assignment . ) consider for example @xmath112 . there are a priori @xmath113 possible patterns of that size , but only 3 inequivalent ones which can not be reduced to configurations having smaller periods ( see figure [ fpatt1 ] ) . the naming convention is to list the colors starting from the lower left corner onwards , denoting bottom1/top1 with 0 , alias white , bottom2/top2 with 1 , alias black . the equivalence of patterns comes from the fact that we have to look at the lattice globally . it is easy to see , for example , that in the case @xmath112 we have the equivalences @xmath114 \simeq [ 0010 ] \simeq [ 1000 ] \simeq [ 0001]$ ] , @xmath115 \simeq [ 1101 ] \simeq [ 0111 ] \simeq [ 1110 ] $ ] , and @xmath116 \simeq [ 0101 ] $ ] ( checkerboard lattice ) . moreover @xmath117 $ ] and @xmath118 $ ] have periods @xmath34 , @xmath119 $ ] and @xmath120 $ ] , @xmath121 $ ] , @xmath122 $ ] have periods @xmath123 and @xmath124 . patterns.[fpatt1],height=75 ]    * claim : * _ some of the distributions are integrable , and some are not . _ although the pattern is 3d consistent , the lax pair is weak and can not precisely fix the bottom and top equations . [ [ times-1-patterns - unicolor - distributions . ] ] @xmath125 patterns ( unicolor distributions ) . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    both unicolor distributions @xmath126 have vanishing entropy . the purely white one is linear . the purely black one is non - trivially integrable , showing quadratic growth of the sequence of degrees @xmath127    [ [ times-2-patterns . ] ] @xmath128 patterns . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    both @xmath128 patterns @xmath129 or @xmath130 , that is to say alternating black and white stripes , are integrable , with quadratic growth of the degrees . [ [ times-2-patterns.-1 ] ] @xmath0 patterns . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    for @xmath112 we have different results for the different patterns in figure [ fpatt1 ] . * both @xmath120 $ ] and @xmath114 $ ] are integrable , with quadratic growth of the degrees . * the calculation of the degrees for @xmath131 $ ] yields the sequence @xmath132 this sequence if fitted by the rational generating function @xmath133 and gives a non vanishing entropy @xmath134 with @xmath135 the largest root of @xmath136 , approximately @xmath137 .    * caveat : * when computing sequences of degrees , one should in principle consider iterations of the whole pattern , but that tends to make the calculations heavier . for the sequence ( [ 1011_seq ] ) , this would mean considering only the subsequence formed by odd terms , leading to a growth given by the maximal root @xmath138 of @xmath139 . of course @xmath140 .    [ [ times-3-patterns . ] ] @xmath141 patterns . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we have examined all the period @xmath142 patterns . the various nonequivalent patterns are depicted in figure [ fpatt23 ] . patterns.[fpatt23],height=264 ]    the computations show that what matters is not just the proportion of black and white cells , but the actual conformation of the pattern . for example the period @xmath142 patterns @xmath143 $ ] and @xmath144 $ ] have an equal number of black and white cells . the first one is integrable ( quadratic growth of the degrees ) while the latter is not , as may be seen from the sequence : @xmath145 this sequence has exponential growth , but is not long enough to determine an exact value of the entropy . the approximate value is @xmath146 .    out of the @xmath147 period @xmath142 nonequivalent patterns , we have one linear case ( all white @xmath148 $ ] ) , eight integrable cases ( @xmath149 $ ] , @xmath150 $ ] , @xmath151 $ ] , @xmath152 $ ] , @xmath153 $ ] , @xmath143 $ ] , @xmath154 $ ] , @xmath155 $ ] ) , and four non integrable ones ( @xmath156 $ ] , @xmath144 $ ] , @xmath157 $ ] , @xmath158 $ ] ) . the following pictures show the aspect of two integrable cases and two non integrable ones . patterns.,height=120 ]     patterns.,height=120 ]    * remark : * the entropy calculations for these patterns can be made equally well with the relations \ { white = bottom1 , black = bottom2 } or with \ { white = top1 , black = top2}. both would give the same results .    from the above results , one may already conclude that random distributions are expected to be non - integrable . we will next push the lax concept and zcc to a @xmath159 sublattice described in figure [ fl2x2 ] . ( such sub - lattices have been discussed previously , e.g. , in @xcite . ) to determine the evolution we need 5 initial values , marked with black disks and expect to get values for the vertices at the open circles . since the lax matrices belong to @xmath160 the zero curvature conditions can provide at most three equations , and thus if everything works well the evolution is determined . it is also clear that this will not work for bigger sub - lattices as we would then need to provide more than 3 values . ( 3482,2513)(0,-10 ) ( 0,2000 ) ( 2000,0 ) ( 1000,0 ) ( 0,1000 ) ( 0,0 ) ( 1000,2000 ) ( 2000,1000 ) ( 2000,2000 )    ( 0,2000)(0,0 ) ( 1000,2000)(1000,0 ) ( 0,2000)(2000,2000 ) ( 2000,2000)(2000,0 ) ( 0,0)(2000,0 ) ( 0,1000)(2000,1000 )    ( 100,40)(0,0)[lb]@xmath161 ( 1050,40)(0,0)[lb]@xmath162 ( -50,2090)(0,0)[lb]@xmath163 ( 1000,2090)(0,0)[lb]@xmath164    ( -1000,450)(0,0)[lb]@xmath165 ( 2050,450)(0,0)[lb]@xmath166    ( -1000,1450)(0,0)[lb]@xmath167 ( 2050,1450)(0,0)[lb]@xmath168    the zero curvature condition for the @xmath159 sublattice is given by @xmath169      in the flipped h1@xmath104 case the lax matrices were given in eq . now using this on a @xmath170 sub - lattice we get just one condition , namely @xmath171 from this one can in principle solve @xmath172 in terms of the other variables . however , this equation does not determine the values for @xmath173 or @xmath174 and therefore these lax matrices fail to give the evolution . [ [ h1_2x2 ] ] h1 ~~    for this basic model the lax matrices were given in eq .  . condition leads to equations that have _ two rational _ solutions : the regular one    [ h1reg ] @xmath175    and an exotic solution    [ h1ex ] @xmath176 } { ( p - r)(x - x_{22})-(q - r)(x - x_{11})-(x - x_{11})(x - x_{22})(x_1-x_2)},\end{aligned}\ ] ]    the regular solution could also be obtained using the evolution on the original lattice , first solving for @xmath10 . as a consequence @xmath173 depends only on @xmath177 and @xmath174 only on @xmath178 . the exotic solution is different , as @xmath173 and @xmath174 both depend on _ all _ initial values . furthermore , it depends on @xmath179 and @xmath180 , and not solely on @xmath181 , as is the case for the regular solution . we may view the variables @xmath182 as associated to the vertices and @xmath183 as associated to the bonds of the @xmath170 sublattice .    in the algebraic entropy analysis the vertex variables are linear and for the bond variables we find the sequence of degrees @xmath184 this sequence can be fitted with the generating function @xmath185 the sequence has quadratic growth , signaling integrability .    what it the nature of the exotic solution ? since the vertex variables have independent linear evolution we can solve the equation with @xmath186 . when this is substituted into the bond equations they give a non - autonomous generalization of a yang - baxter map : using the coarse grained indexing @xmath187 , @xmath188 , i.e. , @xmath189 we have    [ yb - fv ] @xmath190 the solution @xmath191 constant is not allowed and if either @xmath192 or @xmath193 is constant , @xmath194 collapses to @xmath195 . in the generic case , denoting @xmath196 we get @xmath197}{x - y - f(n)+g(m)}.\ ] ] after the further translation @xmath198 , where @xmath91 is a solution of @xmath199 we finally get @xmath200    which is a non - autonomous version of the adler map @xcite ( aka @xmath201 in the classification @xcite ) . the situation can be described by the following diagram : @xmath202    [ [ h3_2x2 ] ] h3 ~~    the phenomenon described in the previous section is not generic . indeed , applying the same coarse - graining to an arbitrary integrable quad - equation will lead to a system having only one rational solution ( the regular one coming from the original lattice ) . we have , however , found more examples where an exotic rational solution exists . here is one , provided by the lattice modified kdv ( lmkdv ) ( aka h3@xmath203 ) . in that case the defining relations of the exotic solution are :    @xmath204    in the algebraic entropy analysis the sequence of degrees for the vertex variables has linear growth as expected , while the sequence for the bonds is the same as for h1 ( see above ) . now the equation on the vertex variables can be solved with @xmath205 and if we introduce @xmath206 we obtain the bond equations in the form @xmath207 using the previously introduced notation . with the further scaling @xmath208 where @xmath91 solves @xmath209 equation reduces to @xmath210 where @xmath211,\quad \beta(n)=\lambda q^2/[g(n)(g(n)-q^2-r^2)].\ ] ] this is nothing but a non - autonomous version of @xmath212 in the classification of @xcite . the diagram presented for h1 works also for h3@xmath213 . we have found that this phenomenon occurs also for the lattice modified kdv and for the lattice schwarzian kdv . we have examined in some detail the properties of the models defined in this way and we will present these results elsewhere @xcite . we have discussed the strength of the lax pairs ( or the zero curvature condition ) and bt through various examples . we have found several cases where the zcc does not uniquely determine the evolution but allows two possibilities . this happens , e.g. , in the flipped h1@xmath104 model . if one then builds an infinite lattice by arbitrarily choosing for each cell one of the two allowed relations , the result is sometimes integrable and sometimes shows nonzero entropy . if the zcc is pushed to a @xmath170 sublattice we get more examples where it is ambiguous and yields both the regular solution as well as an exotic one . the latter can not be generated by some equation in the sublattice , because in the exotic solution the variables @xmath173 and @xmath174 depend on both @xmath214 and @xmath215 which is not possible using the equations on the elementary squares . the equations we have obtained in this way can be interpreted as having vertex and edge variables : the variables with even number of indices live at the vertices while the ones with an odd number of indices live on the edges . the edge variables evolve as in a non - autonomous functional yang - baxter equation . for more details .        99 fw nijhoff , a ramani , b grammaticos and y ohta , _ on discrete painlev equations associated with the lattice kdv systems and the painlev vi equation _ appl . math . * 106 * ( 2001 ) 261 . nijhoff and a. walker , _ the discrete and continuous painlev vi hierarchy and the garnier systems _ , glasgow math . j. * 43a * , 109 ( 2001 ) .      v adler , a bobenko and yu suris , _ classification of integrable equations on quad - graphs . the consistency approach _ * 233 * ( 2003 ) 513 . nijhoff , _ lax pair for the adler ( lattice krichever  novikov ) system _ , phys . lett . * a297 * 49 - 58 ( 2002 ) . j. atkinson . _ bcklund transformations for integrable lattice equations _ j. phys . a : math . theor . * 41 * 135202 ( 2008 ) . p d xenitidis and v g papageorgiou . _ symmetries and integrability of discrete equations defined on a black  white lattice _ , a : math . theor . * 42 * ( 2009 ) 454025 ( 13pp . ) . j. hietarinta and c. viallet . _ searching for integrable lattice maps using factorization _ , j. phys . * 40 * ( 2007 ) 1262912643 v.e . adler , _ recuttings of polygons _ , funct . app . , * 27 * , 141143 ( 1993 ) .
Abstract:

Answer: we consider various 2d lattice equations and their integrability , from the point of view of 3d consistency , lax pairs and bcklund transformations . we show that these concepts , which are associated with integrability , are not strictly equivalent . in the course of our analysis , we introduce a number of black and white lattice models , as well as variants of the functional yang - baxter equation .


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: since the discovery of bose - einstein condensation in magnetically trapped bose gases , the study of the low - energy collective excitations has attracted a big interest both from experimental and theoretical point of view . mean field theory has proven to be a good framework to study static , dynamic and thermodynamic properties of these trapped gases . in particular , it provides predictions of the frequencies of collective excitations that very well agree with the observed ones . recently @xcite the energy shifts and damping rates of these low - lying collective excitations have been measured as a function of temperature . however , these phenomena have not yet been completely understood theoretically .    in this paper we study the influence of thermal excitations on collective oscillations of the condensate in the collisionless regime . previous papers on this subject has been devoted mainly to calculation of the landau damping by means of perturbation theory . in refs.@xcite only the uniform system has been considered , whereas in refs.@xcite landau damping in trapped bose gases has also been studied but using the semiclassical approximation for thermal excitations and the hydrodynamic approximation for collective oscillations . an important point of refs.@xcite is that the authors discuss the possible chaotic behavior of the excitations in an anisotropic trap . the frequency shift has also been studied for a trapped condensate in the collisionless regime in ref.@xcite .    in the present work we study the interaction between collective and thermal excitations using the gross - pitaevskii equation and perturbation theory . we consider spherically symmetric traps , since in this case the spectrum of excitations is easily calculated , avoiding the use of further approximations . even though the case of anisotropic traps can be significantly different in the final results , a detailed investigation of spherical traps is instructive . we explore , in particular , the properties of monopole oscillations by studying the temperature - induced resonances that occur in the condensate when an external perturbation of the trapping frequency is applied and , also , the landau damping associated with the interaction with thermally excited states . this paper is organized as follows . in section ii we introduce the general equations that describe the elementary excitations of the condensate within the bogoliubov theory @xcite . in sec . iii we recall the perturbation theory for a trapped bose - condensed gas in order to study the interaction between elementary excitations . in sec . iv we introduce the linear response function formalism and calculate the response function of the condensate when a small perturbation of the trapping frequency is applied . we derive analytic equations for the response function at zero temperature and treat perturbatively the contribution of the elementary excitations , which is related to landau damping . in sec . v we discuss the main results . we consider a weakly interacting bose - condensed gas confined in an external potential @xmath0 at @xmath1 . the elementary excitations of a degenerate bose gas are associated with the fluctuations of the condensate . at low temperature they are described by the time dependent gross - pitaevskii ( gp ) equation for the order parameter @xcite : @xmath2 where @xmath3 is the number of atoms in the condensate . at zero temperature it coincides with the total number of atoms @xmath4 , except for a very small difference @xmath5 due to the quantum depletion of the condensate . the coupling constant @xmath6 is proportional to the @xmath7-wave scattering length @xmath8 through @xmath9 . in the present work we will discuss the case of positive scattering length , as for @xmath10rb atoms . the trap is included through @xmath0 , which is chosen here in the form of an isotropic harmonic potential : @xmath11 . the harmonic trap provides a typical length scale for the system , @xmath12 . so far experimental traps have axial symmetry , with different radial and axial frequencies , but experiments with spherical traps are also feasible @xcite . the choice here of a spherical trap has two different reasons . first , it greatly reduces the numerical effort and will allow us to study the interaction of oscillations with elementary excitations without any further approximations . second , the energy spectrum of the excitations in such a trap is well resolved yielding to the appearance of well - separated resonances . in anisotropic traps , conversely , the spectrum of excitations is much denser .    the normal modes of the condensate can be found by linearizing equation ( [ tdgp ] ) , i.e. , looking for solutions of the form @xmath13 \label{linearized}\ ] ] where @xmath14 is the chemical potential and functions @xmath15 and @xmath16 are the  particle " and  hole " components characterizing the bogoliubov transformations . after inserting in eq . ( [ tdgp ] ) and retaining terms up to first order in @xmath15 and @xmath16 , one finds three equations . the first one is the nonlinear equation for the order parameter of the ground state , @xmath17 where @xmath18 ; while @xmath19 and @xmath20 obey the following coupled equations  @xcite : @xmath21   u ( { \bf r } ) + g   \psi_0 ^ 2 v ( { \bf r } ) \label{coupled1 } \\ - \hbar \omega v({\bf r } ) & = & [ h_0 - \mu + 2 g \psi_0 ^ 2 ]   v ( { \bf r } ) + g   \psi_0 ^ 2 u ( { \bf r } ) \ ; . \label{coupled2}\end{aligned}\ ] ] numerical solutions of these equations have been found by different authors @xcite . in the present work , we use them to calculate the response function of the condensate under an external perturbation and the landau damping of collective modes . when the adimensional parameter @xmath22 is large , the time - dependent gp equation reduces to the hydrodynamic equations @xcite : @xmath23 @xmath24 where @xmath25 is the particle density and the velocity field is @xmath26 . the static solution of equations ( [ hydron])-([hydrov ] ) gives the thomas - fermi ground state density , which in the spherical symmetric trap reads @xmath27 \label{tf}\ ] ] in the region where @xmath28 , and @xmath29 elsewhere . the chemical potential @xmath14 is fixed by the normalization of the density to the number of particles @xmath30 in the condensate . the density profile ( [ tf ] ) has the form of an inverted parabola , which vanishes at the classical turning point @xmath31 defined by the condition @xmath32 . for a spherical trap , this implies @xmath33 it has been shown @xcite that the hydrodynamic equations ( [ hydron ] ) and ( [ hydrov ] ) correctly reproduce the low - lying normal modes of the trapped gas in the linear regime when @xmath22 is large ( see however ref.@xcite ) . let us briefly recall the perturbation theory for the interaction between collective modes of a condensate and thermal excitations as it was developed in ref.@xcite . suppose that a certain mode of the condensate has been excited and , therefore , it oscillates with the corresponding frequency @xmath34 . we assume that this oscillation is classical , i.e. the number of quanta of oscillation ( @xmath35 ) is very large . then , the energy of the system associated with the occurrence of this classical oscillation can be calculated as @xmath36 with @xmath37 . due to interaction effects , the thermal bath can either absorb or emit quanta of this mode producing a damping of the collective oscillation . the energy loss can be written as @xmath38 where @xmath39 and @xmath40 are the probabilities of absorption and emission of one quantum @xmath41 , respectively . the interaction between excitations is small , so one can use perturbation theory to calculate the probabilities for the transition between a @xmath42-th excitation and a @xmath43-th one , available by thermal activation @xmath44 let @xmath45 and @xmath46 be the corresponding energies and assume @xmath47 . since energy is conserved during the transition process , one has @xmath48 . the interaction term in second quantization is given by @xmath49 in the framework of bogoliubov theory , the field operator @xmath50 can be written as the sum of the condensate wave function @xmath51 , which is the order parameter at equilibrium , and its fluctuations @xmath52 , where @xmath53 [ see eq . ( [ linearized ] ) ] . the fluctuations can be expressed in terms of the annihilation ( @xmath54 ) and creation ( @xmath55 ) operators of the elementary excitations of the system : @xmath56 \ , , \label{uv}\ ] ] where the functions @xmath15 and @xmath16 are properly normalized solutions of equations ( [ coupled1])-([coupled2 ] ) . in the sum ( [ uv ] ) one can select a low energy collective mode , for which we use the notation @xmath57 , and investigate its interaction with higher energy single - particle excitations , for which we use the indices @xmath58 as in ( [ prob ] ) . these latter excitations are assumed to be thermally excited . inserting expression ( [ uv ] ) into eq . ( [ vint ] ) one rewrites the interaction term @xmath59 in terms of the annihilation and creation operators . since we want to study the decay process in which a quantum of oscillation @xmath41 is annihilated ( created ) and the @xmath42-th excitation is transformed into the @xmath43-th one ( or viceversa ) , we will keep only terms linear in @xmath60 and in the product @xmath61 . and the energy conservation during the transition process will be ensured by the delta function @xmath62 . this mechanism is known as landau damping @xcite .    assuming that at equilibrium the states @xmath58 are thermally occupied with the usual bose factor @xmath63^{-1}$ ] , the rate of energy loss can be calculated as @xcite @xmath64 where @xmath65 . \label{matrixel}\end{aligned}\ ] ]    let us define the dissipation rate @xmath66 through the following relation between the energy of the system @xmath67 and its dissipation @xmath68 : @xmath69 using expression ( [ de / dt ] ) @xmath66 can be calculated as @xmath70 where the transition frequencies @xmath71 are positive . the  damping strength " @xmath72 has the dimensions of a frequency . in this work we calculate the quantities @xmath73 by using the numerical solutions @xmath15 and @xmath16 of eqs . ( [ linearized]-[coupled2 ] ) into the integrals ( [ matrixel ] ) . the results will be discussed in section v. the results of the previous section can be used also to study the effect that an external perturbation of the trap has on the collective excitations of the condensate . let us assume the trapping frequency in the form [ @xmath74 , where @xmath75 is a time - dependent modulation . assuming that the perturbation is small , one can use the response function formalism to describe the fluctuations of the system . let us briefly recall the basic formalism @xcite . the behaviour of a system under an external perturbation can be described by studying the fluctuations that may generate the external interaction to a certain physical quantity of the system . an external perturbation acting on the system is described by a new term in the hamiltonian of the type @xmath76 where @xmath77 is the quantum operator of the physical quantity that may fluctuate , and @xmath78 is the  perturbing force " . the mean value @xmath79 is zero in the equilibrium state , in absence of perturbation , and is not zero when it is present . for a periodic perturbation @xmath80 , the relation between @xmath79 and @xmath81 is @xmath82 where @xmath83 is the response function also called generalised susceptibility .    in general @xmath83 is a complex function . it can be seen that the imaginary part of the susceptibility determines the absorption of energy @xmath84 of the external force @xmath85 by the system through the following relation : @xmath86\ , |f|^2 \ , , \label{dissipation}\ ] ] and that the real and imaginary parts of @xmath83 satisfy the kramers - krnig relation @xmath87=\frac{2}{\pi } \,p \int_{0}^{\infty } \frac{im[\alpha(\xi)]}{\xi^2-\omega^2}\,\xi\ , d\xi \ , , \label{kramerskronig}\ ] ] where @xmath88 means the principal value of the integral . the time - dependent external drive @xmath89 induces oscillations of the condensate density @xmath90 with frequency @xmath91 ; @xmath92 . expanding the energy due to the confining potential , @xmath93 , with respect to @xmath89 and @xmath90 one obtains the `` mixed '' term , corresponding to the hamiltonian ( [ v ] ) : @xmath94 comparing it with eq . ( [ v ] ) one can identify the perturbing force and the corresponding coordinate as @xmath95 note that the first order term @xmath96 can be omitted because it gives an additive shift in the hamiltonian which does not contribute to the equations of motion of the system . once we have identified @xmath85 and @xmath97 , we can calculate the response function of the condensate @xmath83 . according to the definition one has @xmath98 let us present the response function in the form @xmath99 , where @xmath100 corresponds to the response function of the condensate at @xmath1 , i.e. , calculated without elementary excitations , and @xmath101 is the contribution of the excitations . at low temperatures it can be assumed that @xmath102 and then @xmath103 can be treated as a perturbation . we proceed as follows . first , we use the hydrodynamic approximation to obtain the response function at @xmath1 . then , within a perturbation theory , we introduce the contribution of the elementary excitations at finite @xmath104 to obtain @xmath101 .      for a spherically symmetric breathing mode @xcite , one can easily prove that the hydrodynamic equations of motion ( [ hydron ] ) and ( [ hydrov ] ) admit analytic solutions of the form @xcite @xmath105 these equations are restricted to the region where @xmath106 . notice that they include the ground state solution ( [ tf ] ) in the thomas - fermi limit . this is recovered by putting @xmath107 , @xmath108 , and @xmath109 . inserting eqs . ( [ n(r , t ) ] ) into the hydrodynamic equations , one obtains two coupled differential equations for the time dependent coefficients @xmath110 and @xmath111 , while at any time @xmath112 is fixed by the normalization of the density to the total number of atoms . the form ( [ n(r , t ) ] ) for the density and velocity distributions is equivalent to a scaling transformation of the order parameter . that is , at each time , the parabolic shape of the density is preserved , while the classical radius @xmath31 , where the density ( [ n(r , t ) ] ) vanishes , scales in time as @xcite @xmath113 where the unperturbed radius @xmath114 is given by eq . ( [ tfr ] ) . the relation between the scaling parameter @xmath115 and the coefficient @xmath110 is @xmath116 . inserting it into eq . ( [ n(r , t ) ] ) we obtain @xmath117 \ , . \label{n(r , t)2}\ ] ] and the hydrodynamic equations then yield @xmath118 and @xmath119 ^ 2 \,b -\frac{\omega_{\rm ho}^2}{b^4}=0\ , . \label{b(t)}\ ] ] the second and third terms of ( [ b(t ) ] ) give the effect of the external trap and of the interatomic forces , respectively . from ( [ r(t ) ] ) and ( [ b(t ) ] ) it follows that at equilibrium @xmath120 and @xmath121 . for a small driving strength @xmath89 , one can assume that the radius of the cloud is perturbed around its equilibrium value , so @xmath122 where @xmath123 it means that @xmath124 is the fractional amplitude of oscillations of the radius and , therefore , it is a measurable quantity .    in the small amplitude limit , one can linearize eq . ( [ b(t ) ] ) with respect to @xmath125 and @xmath124 yielding the following equation @xmath126 the solution is @xmath127 where @xmath128 corresponds to the frequency of the normal mode of monopole in the hydrodynamic limit @xcite . keeping only the lowest order in the small perturbation @xmath124 , eq . ( [ n(r , t)2 ] ) yields @xmath129\delta b\ , , \label{n(r , t)3}\ ] ] and using the thomas - fermi radius at equilibrium ( [ tfr ] ) , it follows that the density fluctuation is given by @xmath130     \delta b(t ) \ , . \label{dn}\end{aligned}\ ] ] we can calculate now @xmath97 using ( [ f ] ) and ( [ dn ] ) , finding @xmath131 where @xmath132 . then , from eq . ( [ db ] ) , one gets @xmath133 at @xmath1 there are no thermally excited states and , hence , @xmath134 . by comparing the definition ( [ alphadef ] ) with ( [ dx3 ] ) one has @xmath135 this is the response function at zero temperature without including any dissipation . therefore @xmath100 is real , i.e. , the induced oscillations at @xmath1 are undamped . the energy of oscillation can be calculated as twice the mean kinetic energy associated to the mode , @xmath136 . for a monopole mode in an isotropic trap , the calculation @xcite gives @xmath137 , where @xmath138 is the amplitude of the oscillation of the cloud ( [ fract ] ) . using eqs . ( [ dx2 ] ) and ( [ alphadef ] ) at @xmath1 , it follows that @xmath139      now we want to calculate the contribution of the thermally excited states to the response function . we study the low temperature regime , where @xmath140 and the energy of oscillation ( [ eosc ] ) can be estimated using @xmath141 instead of @xmath142 . the effect of @xmath103 will be introduced within a perturbation theory . we have already seen that the thermal excitations can either absorb or emit quanta of oscillation @xmath143 and thus they will dissipate energy . the contribution of the elementary excitations to the susceptibility will be a complex function , @xmath144+{\rm i } im[\alpha_1]$ ] , whose imaginary part is related to the absorption of energy @xmath84 of the external perturbation . however , in a stationary solution which is the case under consideration , the absorption @xmath84 must be compensated by the energy dissipation ( [ dampdef ] ) due to the interaction with the elementary excitations . therefore , @xmath145 let us rewrite the definition of the damping rate ( [ dampdef ] ) by using ( [ damping2 ] ) and ( [ dampingik ] ) with a generic oscillation frequency @xmath91 : @xmath146 inserting eq . ( [ eosc ] ) and defining @xmath147 $ ] one obtains the energy dissipation @xmath148 let us recall that the energy dissipation according to eqs . ( [ dissipation ] ) and ( [ q ] ) can be calculated also from the imaginary part of the response function @xmath99 . since @xmath100 is real , eq . ( [ dissipation ] ) becomes @xmath149\ , |f|^2=-\dot{e } \ , . \label{dissipation2}\ ] ] comparing eqs . ( [ edot2 ] ) and ( [ dissipation2 ] ) one can calculate the imaginary part of @xmath101 as @xmath150= 4 \,\frac{15 \mu n}{7 } \sum_{ik }     \gamma_{ik } \,\delta(\omega_{ik}-\omega ) \ , |\beta(\omega)|^2\ , . \label{ima1}\ ] ] and using the kramers - krnig relation ( [ kramerskronig ] ) one finds the real part @xmath151=\frac{8}{\pi } \sum_{ik }      \frac{\omega_{ik } \gamma_{ik}}{\omega_{ik}^2-\omega^2 }      \frac{15 \mu n}{7 } |\beta(\omega_{ik})|^2 \ , . \label{rea1}\ ] ] now we have all the ingredients to calculate the response function of a spherically symmetric trapped condensate when the monopole mode is excited and a small perturbation of the trapping frequency @xmath75 is applied . it can be calculated within first order perturbation as @xmath152 + { \rm i } im[\alpha_1(\omega)]$ ] , by using eqs . ( [ alpha0 ] ) , ( [ rea1 ] ) and ( [ ima1 ] ) , respectively . it is worth stressing that the real part of the susceptibility diverges at @xmath153 ( resonance of the condensate at @xmath1 ) but also at @xmath154 which are the frequencies of the thermal excited modes that due to the interaction are coupled with the monopole . actually , the resonances of the condensate can be found by measuring the fractional amplitude of oscillations of the cloud radius @xmath124 at different perturbing frequencies . this measurable quantity can be easily related to the response function @xmath83 from eqs . ( [ dx2 ] ) and ( [ alphadef ] ) @xmath155    note that the perturbation theory we have used is valid when @xmath156 . this condition becomes very restrictive at @xmath157 near @xmath158 . however it is not difficult to improve the approximation in this region by taking benefit of the analogy between the response function and the green function @xmath159 . it is well known that the green function obeys the dyson equation @xcite which relates the perturbed quantity ( @xmath159 ) and the unperturbed one ( @xmath160 ) through the inversed functions ( @xmath161 and @xmath162 ) in such a way that a perturbation theory for @xmath161 has a more wide applicability that for @xmath159 . analogously , we will find a relation between the inverse response functions , perturbed ( @xmath163 ) and unperturbed ( @xmath164 ) . one has @xmath165 and formally with the same accuracy @xmath166 now the applicability of ( [ alpha-1b ] ) is restricted only by the condition that the second term is small with compare to @xmath167 . it is worth noting that according to the equation ( [ alpha-1b ] ) the poles of @xmath168 related to the resonances are shifted with compare to frequencies @xmath169 and are given by the equation : @xmath170 . however these shifts are very small . in order to present numerical results we choose a gas of @xmath10rb atoms ( scattering length @xmath171 cm ) . for the spherical trap we fix the frequency @xmath172 hz , which is the geometric average of the axial and radial frequencies of ref . @xcite , and corresponds to the oscillator length @xmath173 cm . we solve the linearized gross - pitaevskii equations eqs . ( [ linearized]-[coupled2 ] ) at zero temperature , to obtain the ground state wave function @xmath51 and the spectrum of excited states @xmath174 as well as the corresponding functions @xmath175 . in spherically symmetric traps the eigenfunctions are labeled by @xmath176 , where @xmath177 is the number of nodes in the radial solution , @xmath178 is the orbital angular momentum and @xmath179 its projection . the eigenfunctions are @xmath180 , the energies @xmath181 are @xmath182 degenerate and the occupation of the thermally excited states is fixed by the bose factor .    for a fixed number of trapped atoms , @xmath4 , the number of atoms in the condensate , @xmath30 , depends on temperature @xmath104 . at zero temperature all the atoms are in the condensate , except a negligible quantum depletion @xcite . at finite temperature the condensate atoms coexists with the thermal bath . in the thermodynamic limit @xcite the @xmath104-dependence of the condensate fraction is @xmath183 $ ] . we consider the collective excitations in the collisionless regime . this regime is achieved at low enough temperature . the excitation spectrum at low temperature can be safely calculated by neglecting the coupling between the condensate and thermal atoms @xcite . it means that the excitation energies at a given @xmath104 can be obtained within bogoliubov theory at @xmath1 normalizing the number of condensate atoms to @xmath184 . we investigate the monopole mode ( @xmath185 and @xmath186 ) . the functions @xmath187 and @xmath188 do not present angular dependence , and from eq . ( [ matrixel ] ) it is straightforward to see that the matrix element @xmath189 couples only those energy levels ( @xmath58 ) with the same quantum numbers @xmath178 and @xmath179 . that is , the selection rules corresponding to the monopole - like transition are @xmath190 and @xmath191 . it is obvious , also , that different pairs of levels with the same quantum numbers @xmath177 and @xmath178 but different @xmath179 give the same contribution . therefore , only the integration of the radial part has to be done numerically .    fixed @xmath30 and at a given temperature , we calculate the damping strengths ( [ dampingik ] ) for the transitions @xmath154 coupled with the monopole . in figure 1 we show the values of @xmath73 ( in units of the frequency of the monopole @xmath153 ) for @xmath192 @xmath10rb atoms at @xmath193 . the arrow points to the frequency of the breathing mode @xmath194 , and the chemical potential is @xmath195 [ these values are numerical results of the linearized gross - pitaevskii equations ( [ linearized])-([coupled2 ] ) for @xmath192 rubidium atoms ] . the position of the bars correspond to the allowed transition frequencies @xmath169 ( in units of @xmath196 ) whereas their height defines the numerical value of @xmath73 @xcite . one can see that there are two different types of allowed transitions @xmath154 . the damping strength associated to most of them is very small . conversely , there are a few transitions which give relatively large values of @xmath73 . the latter correspond to transitions between the lowest levels ( @xmath197 ) for different values of @xmath178 ( @xmath198 ) . the main reason for these `` strong transitions '' is that the temperature occupation factor for these low - lying levels is large . moreover the calculation shows that the matrix elements are also enhanced compared to other transitions . this is due to the fact that the radial wave functions involved in the integration have either one ( @xmath199 ) or no node ( @xmath200 ) , differently from the oscillating character of the radial wave functions associated to higher levels @xcite . the contribution of the other transitions is like a small `` background '' which is difficult to resolve in the scale of the figure . a close - up view of the damping strengths of the transition frequencies around the monopole is displayed in the inset of fig .  1 in order to show the dense background . it is worth stressing that such a distinction between `` background '' and `` strong '' transitions depends on the number of condensed atoms in the system and , of course , on temperature . when the number of atoms in the condensate increases , the number of excited states available by thermal excitations also increases , leading to a denser and less resoluble background . in figure 2 we present the same as in fig . 1 but for @xmath201 atoms of rubidium at @xmath193 , where here @xmath202 . in this case , one can see that the difference between the `` strong '' and `` weak '' transitions is not so impressive as in a bigger condensate since all damping strengths can be appreciated in the same scale . we can conclude that at large @xmath30 we have actually two different phenomena . the strong transitions create temperature induced resonances which can be observed in direct experiments . the background transitions give rise to landau damping of the collective oscillations ( see subsection b ) .      using the transition frequencies @xmath154 and the corresponding damping strength @xmath73 , we have calculated the response function @xmath83 . at zero temperature , the response function @xmath100 given by eq . ( [ alpha0 ] ) , gives a resonance at the monopole frequency @xmath203 evaluated in the hydrodynamic regime . due to interaction , thermal excited modes are coupled with the monopole . it means that when one excites the breathing mode of the condensate , the elementary excitations can give rise to other resonances at @xmath154 , which are the frequencies where @xmath204 $ ] diverges [ see eq . ( [ rea1 ] ) ] . we will now discuss the conditions for the observation of these effects in actual experiments . in particular , we calculate the contribution of these resonances to the response function and estimate the associated strengths . let us study the resonances at @xmath205 for @xmath206 atoms of @xmath10rb . the behavior of the damping coefficients @xmath73 is analogous to the one for @xmath207 condensate atoms ( see fig .  1 ) but in this case the difference between `` strong '' resonances and small background is even bigger : the dense background is not more resoluble in the scale of the strong resonances . there are five resonances that stand out the others , and that we label as @xmath208 and @xmath209 the corresponding damping strength ( see table 1 for numerical values ) .    for perturbing frequencies close to the monopole @xmath210 , the monopole susceptibility , eq . ( [ alpha0 ] ) , can be approximated to @xmath211 where @xmath212 . analogously , @xmath213 near each resonance @xmath214 can be presented in the form @xmath215 . the ratio @xmath216 is a measure of the relative intensity between temperature - induced and monopole resonance . table 1 displays the numerical values of the relative intensity for each temperature - induced resonance @xmath208 with respect to the monopole one , for @xmath206 atoms in the condensate at @xmath205 . the relative strenght of the response function ( @xmath216 ) at @xmath208 depends not only on the damping coefficient @xmath209 but also on @xmath217 . it means that one mode @xmath208 will be easier to excite , i.e. , the strength of the response will be bigger , when it is close to the frequency of the monopole . note also that the resonance strength increases with temperature through @xmath209 .    from table 1 one can see that the biggest resonance occurs at @xmath218 which is resoluble from the monopole frequency @xmath219 and has a large enough relative strength to be observed . it means that , tuning the perturbation frequency @xmath91 to this value , a fluctuation of the fractional amplitude of oscillations can be observed .    in figure 3 we have plotted the frequency dependence of the real part of response function @xmath220 calculated according to equation ( [ alpha-1b ] ) for @xmath206 . the response function is given in arbitrary units , and frequency is in units of @xmath196 . the dashed line shows the monopole resonance at @xmath153 , whereas the other divergences of @xmath83 correspond to the temperature - induced resonances at @xmath208 . from this figure one can see that the thermal induced resonances are quite distinct one from the other and from the monopole one . therefore , temperature - induced resonances could be observed in experiments with good enough frequency resolution and good accuracy in the measurement of the radius fluctuations . we would like to stress that the phenomenon we have discussed is related to quite delicate features of interaction between elementary excitations , and therefore , its observation would give rich information about properties of bose - einstein condensed gases at finite temperature .      from fig . 1 one can see that the weak background transitions @xmath154 have , generally speaking , very small frequency separation . to estimate this distance quantitatively let us renominate the resonances by an index @xmath42 in the order of increasing value of @xmath91 . then , one can define the average distance between resonances @xmath221 according to : @xmath222 in a small interval around the collective oscillation @xmath223 , we sum up all the transition frequencies allowed by the monopole selection rules and find the following values for the average distance between two consecutive transition frequencies : @xmath224 and @xmath225 for @xmath206 , @xmath207 and @xmath226 , respectively . it is hopeless , of course , to try to resolve these resonances . actually , there are reasons to believe that these resonances are smoothed and overlapped . first of all , because a real trap can not be exactly isotropic . it means that levels with different @xmath179 have slightly different energies , only levels with @xmath227 are exactly degenerated . therefore , each energy level with a given @xmath178 will be splitted on @xmath228 closer sublevels making more dense the energy spectrum . furthermore , all excitations at finite temperature have associated a finite life time . excitations with @xmath229 which are the ones that mainly contribute in the `` background '' transitions , have the shortest life time . this can be accounted for phenomenologically by assuming that these levels have a finite lorentzian width @xmath230 . that is , instead of delta functions in the equation for the damping rate ( [ damping2 ] ) , we will consider a lorentzian distribution centered at @xmath154 with a fixed width @xmath231 : @xmath232 $ ] . in this case the damping rate becomes a smooth function of @xmath233 and its value when @xmath234 defines the landau damping of the monopole oscillations . at conditions @xmath235 the damping rate will have only a weak dependence on the exact value of @xmath230 . in figure 4 we plot the the dimensionless damping rate @xmath236 as a function of the lorentzian width @xmath230 ( in units of @xmath196 ) for @xmath237 at different temperatures . the summation in ( [ damping2 ] ) has been done over all resonances excluding of course the  strong resonances " presented in figures 1 , 2 and 3 . one can see that the @xmath230-dependence is weak indeed in the interval @xmath238 and @xmath239 can be reliable extrapolated from this interval to the value @xmath240 . we take as landau damping this extrapolated value of @xmath239 . one can estimate the accuracy of this extrapolation procedure to be of the order of 10 % according to the change of @xmath239 over this interval . in figure 5 we plot the damping rate versus @xmath241 for @xmath206 and @xmath207 atoms in the condensate . as expected , landau damping increases with temperature since the number of excitations available at thermal equilibrium is larger when @xmath104 increases . one can distinguish two different regimes in fig . 5 one at very low @xmath104 ( @xmath242 ) and the other at higher @xmath104 . the behaviour of the damping rate becomes linear at relatively small temperature ( @xmath243 ) in comparison to the homogeneous system @xcite where this regime occurs at @xmath244 . moreover , the damping rate increases for larger number of condensed atoms because the density of states available to the system also increases . it is interesting to note that the order of magnitude of the damping rate is the same as the one previously estimated for an uniform gas @xcite and for anisotropic traps @xcite . we have considered the monopole oscillation of a bose - condensed dilute atomic gas in an isotropic trap . first of all , we have calculated the normal modes of the condensate by solving the time - dependent gross - pitaevskii equation within bogoliubov theory @xcite and then we have used the formalism developed in ref.@xcite to calculate the matrix elements associated with the transitions between excited states allowed by the monopole selection rules . within a first order perturbation theory we have studied the landau damping of collective modes due to the coupling with thermal excited levels . we have developed the response function formalism to study the fluctuations of the system due to an external perturbation . the contribution of the elementary excitations has been introduced also perturbatively as in the calculation of the damping strength , and we have derived analytic equations for the response function at zero temperature and at low temperature regime . we have seen that when the condensate oscillates with the monopole mode and a small perturbation to the trap frequency is applied , one can excite new resonances at the transition frequencies . these thermal - induced resonances are coupled with the monopole due to interaction effects . one can not exclude _ a priori _ the possibility to observe such resonances also in anisotropic traps . this problem deserves further investigation . observation of these resonances would give important and unique information about the interaction between elementary excitations in bose - einstein condensed gases . m. edwards , p.  a.  ruprecht , k.  burnett , r.  j.  dodd , and c.  w.  clark , phys . rev . lett . * 77 * , 1671 ( 1996 ) ; p. a. ruprecht , mark edwards , k. burnett , and charles w. clark , phys . a * 54 * , 4178 ( 1996 ) ; m. edwards , r.  j.  dodd , c.  w.  clark , and k. burnett , j. res . * 101 * , 553 ( 1996 ) . when @xmath245 is large one can use the hydrodynamic approximation for the functions @xmath15 and @xmath16 of the low - energy excitations [ see l.p . pitaevskii , _ recent progress in many - body theories _ , ed . d. nielson and r. bishop , world scientific ( singapore , 1998 ) , p.3 ] . however , we have found that the matrix elements @xmath189 , defined in section iii , are quite sensitive to the accuracy of these functions . hydrodynamics can give values for @xmath73 an order of magnitude smaller than the values calculated numerically by using the bogoliubov functions ( [ coupled1],[coupled2 ] ) . on the contrary , it is completely safe to use hydrodynamics to calculate the response function at zero temperature @xmath246 ( see section iv ) , since the difference between hydrodynamics and exact values of the monopole frequency @xmath153 is very small . beliaev decay of an elementary excitation into a pair of excitations [ s.t . beliaev , soviet phys . jetp * 34 * ( 1958 ) 323 ] is not active for the lowest energy modes in the case of trapping potential because of the discretization of levels . we have neglected all the transitions that involve the lowest dipole mode ( @xmath247 , @xmath248 ) because in an external harmonic potential this mode , corresponding to the oscillation of the center of mass , is unaffected by the interatomic forces and then the transition probability due to interaction effects ( [ prob ] ) must be zero @xcite . however , the perturbation formalism we have presented in section iii does not take into account this physical consideration automatically and therefore , we have omitted by hand all transitions with the dipole mode .      . damping coefficients @xmath249 ( in units of @xmath196 ) of the  strong resonances " @xmath208 ( in units of @xmath196 ) and relative intensities @xmath250 between temperature - induced the monopole resonance , for @xmath206 condensate atoms of @xmath10rb in a spherical trap with @xmath251 cm at @xmath193 . [ cols="^,^,^ " , ]
Abstract:

Answer: interaction between collective monopole oscillations of a trapped bose - einstein condensate and thermal excitations is investigated by means of perturbation theory . we assume spherical symmetry to calculate the matrix elements by solving the linearized gross - pitaevskii equations . we use them to study the resonances of the condensate induced by temperature when an external perturbation of the trapping frequency is applied and to calculate the landau damping of the oscillations .