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: transverse collective oscillations have been recently excited in the radial plane of an elongated cylindrically symmetric condensate @xcite . they are quasi-2d excitations and indeed , at very low temperatures , the measured frequencies of the transverse monopole and quadrupole modes are found to be in very good agreement with the theoretical predictions obtained by linearizing the gross - pitaevskii equation for a trapped two - dimensional ( 2d ) system . moreover , it has been experimentally found that the transverse breathing mode of an elongated condensate exhibits unique features . first , the excitation frequency is close to twice the radial trapping frequency , it is almost independent of the strength of the two - body interaction and the number of particles , and it is nearly independent of temperature . and second , the transverse breathing mode has a very small damping rate compared to other modes . the measured damping rate of the transverse quadrupole mode is approximately one order of magnitude larger than the one of the transverse breathing mode . this is in contrast to the 3d case , where both monopole and quadrupole modes have similar decay rates . these striking properties experimentally found for the transverse breathing mode are in agreement with the unique features predicted for the breathing mode in 2d isotropic condensates @xcite : its frequency is a universal quantity and it is an undamped mode . one can conclude that these features are characteristic of the reduced dimensionality thus , providing a demonstration of the two - dimensional nature of the transverse oscillations of an elongated condensate . therefore , a highly anisotropic cigar - shaped trap can be considered as a first approximation as an infinite cylinder along the @xmath0-direction .    in the present paper we consider a cylindrical condensate and calculate at zero temperature the transverse spectrum of excitations , by solving the linearized gross - pitaevskii equation in 2d , and the modes with non - zero momentum @xmath1 along the longitudinal axis , as well . the main purpose is to investigate the coupling between quasi-2d ( transverse ) collective oscillations and thermal excitations in the collisionless regime . we study the landau damping as a decay mechanism of transverse collective modes as a function of temperature , in which a transverse collective oscillation is annihilated in a collision with a thermal excitation to give rise to another excitation . this paper is organized as follows . in sec . ii we calculate the transverse spectrum of excitations of a cylindrical condensate by solving the linearized gross - pitaevskii equation in 2d within the bogoliubov theory . we obtain also the modes with non - zero @xmath1 . in sec . iii we briefly recall the perturbation theory for the interaction between collective modes of a condensate and thermal excitations ( developed in ref . we apply this theory to calculate the decay rate due to landau damping of a transverse oscillation . in sec . iv we discuss the main results for the transverse quadrupole mode . the situation with the singular transverse breathing mode is also discussed . at low temperature , the elementary excitations of a trapped weakly interacting degenerate bose gas are described by the time - dependent gross - pitaevskii ( gp ) equation for the order parameter : @xmath2 where @xmath3 is the interaction coupling constant , fixed by the @xmath4-wave scattering length @xmath5 , @xmath6 is the atomic mass and the order parameter is normalized to the number of atoms in the condensate @xmath7 . the confining potential is usually cylindrically symmetric and it is given by @xmath8 , with @xmath9 . the ratio between axial ( @xmath10 ) and radial ( @xmath11 ) trapping frequencies defines the anisotropy parameter of the trap @xmath12 . experimentally , in ref . @xcite , the transverse modes have been excited in the radial plane of a highly anisotropic cigar - shaped trap , with @xmath13 . thus , as a reasonable approximation , we consider an idealized cylindrical trap that is uniform in the @xmath0-direction ( @xmath14 ) and has isotropic trapping potential in the radial plane @xmath15 the harmonic trap frequency @xmath16 provides a typical length scale for the system , @xmath17 . we assume a cylindrical condensate with large longitudinal size @xmath18 compared to the radial size @xmath19 , and with a number of atoms per unit length @xmath20 .    at low temperatures , the dynamics of the condensate is described by the linearized time - dependent gp equation . assuming small oscillations of the order parameter around the ground - state value @xmath21 \ , ,   \label{fluctuation}\ ] ] where @xmath22 is the chemical potential , the normal modes of the condensate can be found by seeking fluctuations such as @xmath23 the functions @xmath24 and @xmath25 are the  particle " and  hole " components characterizing the bogoliubov transformations , and @xmath26 is the frequency of the corresponding excited state . the transverse excitations of a highly elongated condensate can be obtained by neglecting oscillations along the longitudinal direction and thus , looking for fluctuations in the transverse plane of the form @xcite @xmath27 where we have used cylindrical coordinates ( @xmath28 ) . transverse excitations are quasi-2d modes . to calculate landau damping we must also consider thermal excitations having a quasi - continuous spectrum at @xmath29 with the finite longitudinal wave vector @xmath1 along @xmath0 direction . to study these modes we look for solutions that correspond to propagating waves of the form @xcite @xmath30    inserting ( [ fluctuation ] ) , ( [ deltapsi ] ) and ( [ deltapsilong ] ) in eq . ( [ tdgp ] ) and retaining terms up to first order in @xmath24 and @xmath25 , it follows three equations . the first one is the nonlinear equation for the order parameter of the ground state , @xmath31 \psi _ { 0}(r_{\perp } ) = \mu \psi _ { 0}(r_{\perp } ) \ , ,   \label{groundstate}\ ] ] where @xmath32 , and @xmath33 ; while @xmath34 and @xmath35 obey the following coupled equations  @xcite : @xmath36 u_{k}(r_{\perp } )   \nonumber \\ & & + g\psi _ { 0}^{2}v_{k}(r_{\perp } )   \label{coupled1 } \\ -\hbar \omega v_{k}(r_{\perp } ) & = & \left [ h_{0}+\frac{\hbar ^{2}}{2m}\left ( \frac{m^{2}}{r_{\perp } ^{2}}+k^{2}\right ) -\mu + 2g\psi _ { 0}^{2}\right ] v_{k}(r_{\perp } )   \nonumber \\ & & + g\psi _ { 0}^{2}u_{k}(r_{\perp } ) \;.   \label{coupled2}\end{aligned}\ ] ] we solve the linearized gp equations eqs . ( [ groundstate]-[coupled2 ] ) at zero temperature , to obtain the ground state wave function @xmath37 , the eigenfrequencies @xmath38 and hence the energies @xmath39 of the excitations , as well as the corresponding functions ( @xmath40 ) .    in a cylindrical trap the modes of the condensate are labeled by @xmath41 , where @xmath42 is the number of nodes in the radial solution , @xmath43 is the @xmath0-projection of the orbital angular momentum , and @xmath1 is the longitudinal wave vector . the quantum numbers @xmath42 and @xmath43 are discrete ( @xmath44 , and @xmath45 ) whereas @xmath1 for highly elongated condensates ( @xmath46 ) is a real number , and thus , it is a continuous index . moreover , it can be seen from eqs . ( [ coupled1 ] ) and ( [ coupled2 ] ) , that levels with @xmath47 or @xmath48 are exactly degenerated . the excitation spectrum of a cylindrical condensate presents different branches labeled by ( @xmath49 ) that are continuous and increasing functions of @xmath50 @xcite . for a given mode , i.e. , fixed @xmath42 and @xmath43 , the frequency of excitation is @xmath38 . we have numerically checked from eqs . ( [ coupled1 ] ) and ( [ coupled2 ] ) that the state with ( @xmath51 ) and @xmath52 corresponds to the transverse breathing mode with excitation frequency equal to @xmath53 , and that it is independent of the strength of the two - body interaction potential and the number of particles . of course , according to the universal nature of this mode noted above , this general result is also valid in the thomas - fermi ( tf ) limit , where the dispersion relation can be found analytically@xcite . however , the functions @xmath24 and @xmath25 for this mode do not coincide in general with their thomas - fermi limit .    the state with ( @xmath54 ) and @xmath52 is a transverse quadrupole mode whose frequency depends on the two - body interaction , ranging from the tf value @xmath55 , valid for large condensates ( @xmath56 or equivalently @xmath57 ) , to the noninteracting value @xmath53 . therefore , the transverse quadrupole has to be obtained numerically by solving eqs . ( [ coupled1 ] ) and ( [ coupled2 ] ) in each particular value of the dimensionless parameter @xmath58 . in the present work , we use the numerical solutions to calculate the landau damping of collective modes . it is well known that the dipole mode ( @xmath59 ) is unaffected by two - body interactions due to the translational invariance of the interatomic force which can not affect the motion of the center of mass @xcite and its excitation frequency is @xmath60 . it is worth stressing that the independence of the interaction of the transverse monopole mode , is a unique property of 2d systems , related to the presence of a hidden symmetry of the problem described by the two - dimensional lorentz group so(2,1 ) @xcite . let us consider a collective mode with frequency @xmath61 , and the excited states @xmath62 available by thermal activation , with energies @xmath63 and @xmath64 , respectively . suppose that this collective mode has been excited and , therefore , the condensate oscillates with the corresponding frequency @xmath61 . due to interaction effects , the thermal cloud of excitations can either absorb or emit quanta of this mode producing a damping of the collective oscillation . we want to study the decay process in which a quantum of oscillation @xmath65 is annihilated ( created ) and the @xmath66-th excitation is transformed into the @xmath67-th one ( or viceversa ) . the energy is conserved during the transition process , therefore @xmath68 where we assume that @xmath69 . this mechanism is known as landau damping . another possible decay mechanism , also due to the coupling between collective and thermal excitations , is the beliaev damping @xcite , which is based on the decay of an elementary excitation into a pair of excitations . let us define the dissipation rate @xmath70 through the following relation between the energy of the system @xmath71 and its dissipation @xmath72 : @xmath73 assuming that the damping is small , one can use perturbation theory to calculate the probabilities for the transition between the @xmath66-th excitation and the @xmath74-th one , both available by thermal activation , yielding the following expression for the landau damping rate @xcite : @xmath75 where the transition frequencies @xmath76 are positive , and the delta function ensures the energy conservation during the transition process . we have assumed that the thermal cloud is at thermodynamic equilibrium and the states @xmath62 are thermally occupied with the usual bose factor @xmath77^{-1}$ ] . the `` damping strength '' has the dimensions of a frequency and is given by @xmath78 the matrix element that couples the low - energy collective mode ( @xmath79 ) with the higher energy single - particle excitations ( for which we use the indices @xmath62 ) is @xcite @xmath80 . \label{matrixel}\end{aligned}\ ] ]    in this work we calculate the quantities @xmath81 by using the numerical solutions @xmath24 and @xmath25 of eqs . ( [ groundstate]-[coupled2 ] ) into the integrals ( [ matrixel ] ) , avoiding the use of further approximations in the spectrum of excitations . in cylindrically symmetric traps , @xmath82 and from ( [ deltapsilong ] ) the longitudinal excitations are @xmath83 and the energies @xmath84 are not degenerate . we are interested in the decay of transverse collective excitations ( [ deltapsitrans ] ) @xmath85 , due to the coupling with thermally excited states ( @xmath62 ) . the functions @xmath86 and @xmath87 do not depend on @xmath0 , and from eq . ( [ matrixel ] ) it is straightforward to see that the matrix element @xmath88 couples only energy levels ( @xmath62 ) with the same quantum number @xmath89 , and @xmath90 . to calculate the damping rate @xmath91 of the collective oscillation we have to sum up over all couples of excited states ( @xmath62 ) , or equivalently , over all pairs of sets ( @xmath92 ) , that satisfy the above conditions . since @xmath1 is a continuous quantum number , the sum over @xmath1 becomes and integral and eq . ( [ damping2 ] ) yields @xmath93 where the sum is restricted to pairs that satisfy @xmath94 . the frequency of a mode is a continuous function of @xmath1 , therefore the transition frequency is @xmath95 . then , in eq . ( [ gamma2 ] ) the integration over @xmath1 yields @xmath96 ^{-1 }   \label{gamma3}\ ] ] where @xmath97 is the wave vector in which the conservation of energy is verified , i.e. , @xmath98 , and we have used the known property of delta functions @xmath99 we want to calculate the characteristic decay rates due to landau damping of the transverse low - lying collective modes of an elongated bose - einstein condensate as a function of @xmath100 . we consider the collective excitations in the collisionless regime , which is achieved at low enough temperatures .    for a fixed number of trapped atoms , the number of atoms in the condensate depends on temperature . at zero temperature the quantum depletion is negligible and all the atoms can be assumed to be in the condensate @xcite . at finite temperature the condensate atoms coexist with the thermal bath . however , at low enough temperatures the excitation spectrum can be safely calculated by neglecting the coupling between the condensate and thermal atoms @xcite . it means that the excitation energies at a given @xmath100 can be obtained within bogoliubov theory at @xmath101 normalizing the number of condensate atoms to the corresponding condensate fraction at that temperature and assuming that they are thermally occupied with the bose factor .    in order to present numerical results we choose a gas of @xmath102rb atoms ( scattering length @xmath103 cm ) confined in an elongated cylindrical trap with frequencies @xmath104 and @xmath105 hz , that corresponds to an oscillator length @xmath106  cm . the number of condensate atoms per unit length is taken to be @xmath107 . these conditions are close to the experimental ones of ref . @xcite . we have solved the linearized gp equations to obtain an exact description of the ground state @xmath108 and the normal modes of the condensate within bogoliubov theory without using the thomas - fermi or hartree - fock approximations . we have obtained the following numerical results : @xmath109 , and the excitation frequencies of the transverse monopole @xmath110 , and of the transverse quadrupole mode @xmath111 . we have calculated the branches of the excitation spectrum of the cylindrical condensate , labeled by ( @xmath49 ) as a function of @xmath1 .    to calculate the damping rate of a collective mode @xmath112 we have to obtain , first , all pairs of levels ( @xmath62 ) of the excitation spectrum that satisfy both the energy conservation of the transition process @xmath113 , and the selection rules given by the particular collective mode under study . and then , calculate the corresponding matrix elements . we have restricted our calculations to levels with energy @xmath114 . the contribution of higher excited levels can be neglected since their occupation becomes negligible in the range of temperatures we have considered . it is worth recalling that the dipole mode is undamped @xcite since it is not affected by two - body interactions . we have checked that within our formalism the landau damping of the transverse dipole mode is zero , as expected . first of all , we have calculated the landau damping of the transverse quadrupole mode ( @xmath52 , @xmath115 and @xmath116 ) as a function of @xmath100 , @xmath117 . for such a mode , the functions @xmath118 and @xmath87 depend on the radial coordinate @xmath119 and have also angular dependence ( @xmath120 ) . therefore , from eq . ( [ matrixel ] ) it follows the selection rules corresponding to transverse quadrupole - like transitions : @xmath121 and @xmath122 . only the energy levels ( @xmath62 ) that satisfy these conditions can be coupled by the interaction . the integration of the radial part has to be done numerically .    fixed the number of condensate atoms , we have calculated the landau damping for the transverse quadrupole mode as a function of temperature . at zero temperature we have obtained the excitation spectrum within bogoliubov theory . we have found that there are @xmath123 pairs of levels that verify the transverse quadrupole selection rules and the energy condition ( [ selection ] ) at finite wave vectors , i.e. , at @xmath124 . we have calculated the matrix element @xmath88 for each allowed transition , and then , the corresponding damping strength @xmath125 at a given @xmath100 . summing up the damping strengths for all transitions , we have obtained the damping rate for the transverse quadrupole oscillation , due to the coupling with thermal excitations .    in figure [ fig1 ] we plot the damping rate versus @xmath126 for the transverse quadrupole mode for a fixed number of condensate atoms per unit length in a cylindrical trap . as expected , landau damping increases with temperature since the number of excitations available at thermal equilibrium is larger when @xmath100 increases . one can distinguish two different regimes , one at very low @xmath100 ( @xmath127 ) and the other at higher @xmath128 where damping is linear with temperature . however , this linear behaviour appears at relatively small temperature ( @xmath129 ) in comparison to the homogeneous system @xcite where this regime occurs at @xmath130 . for the highly elongated condensate of ref . @xcite , it has been measured at @xmath131 a quality factor @xcite for the transverse quadrupole mode @xmath132 which gives a damping rate @xmath133 . this value is in agreement with our calculations of the landau damping for this mode ( see fig .  1 ) . a more quantitative comparison is not possible at present due to the absence of detailed measurements of the damping of this mode . however , it is reasonable to conclude that the transverse quadrupole mode in elongated condensates decays via landau damping . it is worth noting that the order of magnitude of the damping rate of the transverse quadrupole mode is the same as the one previously estimated for a uniform gas @xcite , for spherical traps @xcite and for anisotropic traps @xcite as well . the landau damping for other modes can be calculated analogously . however , the physical situation for the decay of the singular transverse monopole mode turns out to be much more complicated . formally the calculations can be committed in the same way . the quantum numbers for this mode are @xmath52 , @xmath134 and @xmath135 . the functions @xmath136 and @xmath87 depend only on the radial coordinate @xmath137 , and the matrix element @xmath88 couples only energy levels ( @xmath138 ) with the same quantum numbers @xmath43 and @xmath1 . that is , the selection rules corresponding to transverse monopole - like transitions are @xmath122 ( as we have already discussed in the previous section ) and @xmath139 .    from the spectrum of elementary excitations , we have found that there are only @xmath140 pairs of levels that verify the transverse monopole selection rules and the energy condition ( [ selection ] ) at @xmath141 . the damping strengths associated to the allowed transitions are smaller than the @xmath125 values associated to quadrupole - like transitions . as a result from our calculations the damping occurs quite small . at @xmath142 we obtain a landau damping @xmath143 that is one order of magnitude smaller than the experimental decay @xcite @xmath144 measured in @xcite . thus , landau damping can not explain the experimentally observed decay of the transverse breathing mode .    actually , there are reasons to believe that even this small value of landau damping obtained in the present calculations , could be larger than the true one . the point is that in our calculations , based on perturbation theory , the thermal cloud is assumed to be in the state of thermodynamic equilibrium and thus one can neglect its motion . this is correct for the case of quadrupole oscillations , where the frequencies of the collective motion of the condensate and the cloud are different . however , due to the unique peculiar nature of the 2d breathing mode , the breathing oscillation of the thermal cloud has the same frequency @xmath145 as the condensate . then , if the oscillation is excited by the deformation of the trap , as it takes place in the present experiment @xcite , the cloud will oscillate in phase with the condensate what will further decrease the landau damping @xcite . therefore , this quasi-2d monopole mode in highly elongated condensates may decay via another damping mechanism , different from the landau damping .    after finishing this work it has appeared a preprint by jackson and zaremba @xcite in which the effect of the collisions between elementary excitations has been taken also into account . this effect is usually small in comparison with the landau damping . nevertheless , for this peculiar mode the collisions seems to be the dominant damping mechanism , yielding a decay rate that it is in agreement with the experimentally measured value @xcite . note , however , that the experiments are produced at quite high amplitude of oscillations and non - linear effects can play an important role . for example , as it was noted in @xcite , the transverse breathing mode exhibits a  parametric  instability due to the decay into two or more excitations with non - zero momentum along the axis . in ref . @xcite it has been considered the process in which the transverse breathing oscillation decays into two excitations propagating along @xmath0 with momentum @xmath1 and @xmath146 . this damping mechanism is active also at zero temperature . but in this case the damping rate depends on the amplitude of the breathing oscillation . we have investigated the decay of low - lying transverse oscillations of a large cylindrical condensate . first of all , we have calculated the transverse normal modes , as well as the excitations with finite @xmath1 , of the condensate by solving the linearized time - dependent gross - pitaevskii equation . and then , within the formalism of ref . @xcite , we have calculated numerically the matrix elements associated with the transitions between excited states allowed by the selection rules of the transverse collective modes . within a first - order perturbation theory and assuming the thermal cloud to be in thermodynamic equilibrium , we have studied the landau damping of transverse collective modes due to the coupling with thermal excited levels as a function of temperature . for the damping rate of the transverse quadrupole mode , we have found a good agreement with the experimentally measured value of ref .  @xcite . whereas our result for the transverse breathing mode is one order of magnitude smaller than the experimental decay . one can conclude that in a highly elongated cylindrically symmetric condensate , transverse quadrupole oscillations decay via landau damping mechanism , but the transverse breathing mode which has an anomalously small measured damping rate , may decay via other damping mechanisms . we thank j.  dalibard , g.  shlyapnikov and e.  zaremba for useful discussions . m.  g. thanks the `` ramon y cajal program '' of the mcyt ( spain ) for financial support .                                        note that the  perturbative  calculations of the landau damping can correspond to a different type of experiment , in which the oscillations of the condensate are excited by the action of a thin laser beam on the condensate , but not on the thermal cloud . in this case the  perturbative  calculations will describe the initial stage of the process , when the cloud is still at rest . as it was presented in the lecture by e. zaremba @xcite , the landau damping corresponding to this type of experiment was calculated using the method developed in ref . @xcite in an elongated condensate with the same experimental parameters of ref . however , they have obtained larger values for the landau damping of the transverse breathing mode in comparison with our results . this discrepancy could be related to the semiclassical nature of their method @xcite , which does not take into account the discrete nature of the energy spectrum , and also to the fact that they consider a cigar - shaped condensate instead of an infinite cylinder as in the present calculations . the finite value of the longitudinal frequency @xmath147 could be important for this transverse breathing mode . therefore , this is still an interesting open problem .
Abstract:

Answer: we study the interaction between low - lying transverse collective oscillations and thermal excitations of an elongated bose - einstein condensate by means of perturbation theory . we consider a cylindrically trapped condensate and calculate the transverse elementary excitations at zero temperature by solving the linearized gross - pitaevskii equations in two dimensions . we use them to calculate the matrix elements between thermal excited states coupled with the quasi-2d collective modes . the landau damping of transverse collective modes is investigated as a function of temperature . at low temperatures , the damping rate due to the landau decay mechanism is in agreement with the experimental data for the decay of the transverse quadrupole mode , but it is too small to explain the slow experimental decay of the transverse breathing mode . the reason for this discrepancy is discussed .


Instruction: Write an abstract for the following acadamic paper.

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

Paper: the most trivial disorder originates in media due to random static impurities ( see , e.g. @xcite ) . however , another but very interesting possibilities are glasses , which have a huge number of ( meta)stable degenerated states . originally glassy system has a homogeneous hamiltonian with no intrinsic random parameter , and disorder occurs in it dynamically . recently @xcite we have demonstrated that a popular frenkel - kontorova model @xcite presents an example of glassy system , which has a lot of static states , known as _ configurational excitations _ of the classical ground state , with energy splitting extremely ( exponentially ) small . as it was shown in @xcite , this model has a nontrivial quantum dynamics , the quantum phase transition : if quantum parameter exceeds some critical value , the `` pinned '' glassy phase turns into `` sliding '' phonon gas .    the frenkel - kontorova model ( fk ) @xcite is widely used @xcite in the solid state physics to get insight on generic properties of noncommesurate systems . its ground state , which is rather quasiperiodical @xcite than periodical , attracts also an attention @xcite as some interplay @xcite between order and disorder @xcite . this model describes a chain of atoms / particles interacting with elastic forces , placed in periodic potential , which period differs from a mean interparticle distance . the ground state ( gs ) of this model is defined as a static , equilibrium configuration of the chain , that corresponds to the _ absolute _ minimum of the chain potential energy . the ground state is unique and has some special order of particles in the chain , that was discovered by aubry @xcite more than twenty years ago . the positions of atoms in the chain are described by an area preserving map , which is well known in the field of dynamical chaos as the chirikov standard map @xcite . the ratio of the mean interparticle distance to a period of the external potential in the fk model determines the rotation number of the invariant curves of the map , while the amplitude of the periodic potential gives the value of the dimensionless parameter @xmath0 . for @xmath1 the kam curves are smooth and the spectrum of longwave phonon excitations in the chain is characterized by a linear dispersion law starting from zero frequency . in this regime the chain can freely slide along the external field ( the `` sliding '' phase ) . on the contrary , for @xmath2 the kam curves are destroyed and replaced by an invariant cantor set , which is called cantorus . in this regime the phonon spectrum has a gap , and the chain is pinned ( `` pinned '' phase ) . later , on the example of ising spin model to which the fk model can be approximately reduced @xcite it has been shown @xcite that the gs has some well defined hierarchical structure , which particular detailes are determined by number properties of the ratio of the mean interparticle distance to the period of the external field . recently our numerical study @xcite of the original fk model in the pinned phase has shown , that the gs has indeed an hierarchical structure , but in some important detailes different from predicted in @xcite .    in short , we put our attention to the striking fact , that in the pinned phase of the fk chain there are some particles , which positions exponentially close to bottoms of wells of the external potential : corresponding external force acting to such particle is extremely close to zero . obviously , in static equilibrium each of these particles can be considered as some dummy `` glue '' that only couples two adjacent parts of the chain , which ends has almost identical ( exponentially close ) tension forces . another important observation is that small deviations of glue particles from their well bottoms are _ groupped _ into well defined hierarchically ordered scales . now , if one cut the chain into fragments via glue particles , which belong to level with least deviation from well bottoms , one gets several fragments of _ two _ sizes , or two species of some `` bricks '' . the most simple picture quoted here assumes that this rotational number approximates the golden mean value . ] . then one may repeat the procedure , cutting `` bricks '' via `` glue '' particles that belong to next scale of deviations and getting two new species of smaller bricks , and so far . at the last step one gets two species of smallest possible bricks with _ no _ glue particle inside : smaller bricks @xmath3 which consist of 2 particles inside a single well , and larger brick @xmath4 which consist of 4 particles ( two pairs in two adjacent wells ) . in this way one gets an hierarchically ordered set of brick species @xmath5 with very simple composition rules @xcite :    @xmath6    where symbol @xmath7 denotes of an insertion of glue particle , which `` glues '' two adjacent bricks the difference of tension forces at boundaries of bricks @xmath8 and @xmath9 is exponentially small and decrease rapidly with a number of the hierarcy level . in principle , these rules are sufficient to construct a gs for a fk chain of any length , if the rotation number parameter of the chain approximates the mean golden value @xmath10 . besides the gs there exists `` configurational excitation states '' ( ces ) , presented by static equilibrium configurations corresponding to _ local _ ( rather than absolute ) minima of the chain potential , with energy very close to gs . within the picture just outlined above ces correspond to different permutations of bricks @xcite , therefore the number of them can be combinatorically huge . at any accessible small temperature their contributions can dominate over the contribution of gs .    in this paper we address to phonon excitations of the chain , small vibrations around static gs and ces configurations of the chain . these excitations are relevant as for heat transport properties @xcite of the chain , as for some quantum effects @xcite , especially in the quasiclassical limit . as in the previous paper , we concentarte on the case of pinned phase of the chain , which corresponds to a nonzero phonon gap . we start with analysis of the structure of the phonon frequency spectrum in the gs . it is well known , that this spectrum is splitted into bands @xcite but , to our knowledge , up to now there is no clear explanation for origin of its splitting into particular bands and subbands . we show , that this splitting is a direct consequence of particular spatial structure of the chain in its gs . we have found , that this structure is also hierarchically ordered , with definite resemblance and distinctions with respect to a spatial structure of underlying gs . localisation properties of phonons in incommensurate one - dimentional chains are intesively studied in recent works @xcite , with strong indications @xcite , that phonon modes in the gs of fk chain are not localized , and even at edges of frequency bands they are rather _ prelocalized _ , than localized . we study also phonon properties as in gs , as in ces of the chain . our results confirm , that in the gs phonon modes are only prelocalized . however , the situation appears to be quite different for ces . even for ces , which energy splitting ( in natural problem scale ) @xmath11 there are phonon modes , which are localized exponentially . moreover , for ces with higher splitting @xmath12 we see , that there are entire bands of exponentially localized phonon modes . the hamiltonian of the fk model is    @xmath13    the first term in the hamiltonian is a kinetic energy , where we put masses of particles @xmath14 , the second term describes interparticle interaction with elasticity coefficient put to unity , while the third term corresponds to particle interaction with external periodical field with coupling constant @xmath15 . all @xmath16 particles are distributed over @xmath17 period / wells of the external potential , which period , without any loss of generality is taken equal to @xmath18 . the ratio @xmath19 gives @xcite the rotational number of corresponding standard map @xcite . we assume periodical boundary conditions : @xmath20 , @xmath21 , where @xmath22 is the length of the chain . in our subsequent analysis we take ( as some _ typical _ example of fk chain ) the chain with @xmath23 as an approximation of the golden mean value @xmath24 , and parameter @xmath25 well above the critical value @xmath26@xcite . the technique to obtain gs and ces is described in our previous work @xcite . let us stress , that validity of the gs can be checked either by direct analysis @xcite of its spatial structure , or proved by monotonous behavior @xcite of its hall function . note , that the gs is almost degenerate with enormous number of ces : there are hundreds of states with @xmath27 , while the number of ecs with @xmath11 exceeds @xmath28 ! effective hamiltonian for phonon modes can be obtained by expansion up to second order terms of the original fk chain hamiltonian ( [ hfk ] ) around a chosen static configuration of the chain : @xmath29 where @xmath30 are small deviations of particles from their equilibrium static positions @xmath31 , and @xmath32 are corresponding particle momenta . the elasticity matrix @xmath33 solving the eigenvalue problem @xmath34 numerically , we get both the spectrum of phonon frequencies and corresponding vectors of phonon modes . it is well known for a long time @xcite , that frequency spectrum of phonons in the frenkel - kontorova chain is splitted into several bands . main features of the spectrum for a chain in the gs look as very universal : ( i ) the number of main bands are independent of the parameter @xmath15 as well as of a length of the chain ; ( ii ) it was also noticed in @xcite ( without any explanation , either ) that boundaries of the bands correspond to fibonacci numbers ; ( iii ) replacement in interparticle interactions elastic forces by the lenard - jones potential does not change qualitatively a pattern of splitting phonon spectrum into band @xcite . now we can explain these and more features as direct consequence of particular spatial structure @xcite of the gs , which main relevant detailes are summarized in introduction . the key point of our explanation is that the gs consists of large number of almost identical elements . for clarity sake , we consider a particular case of the fk chain with @xmath23 , but all our arguments can be easily generalizied for a chain of different size . applying the picture of cs spatial structure outlined in introduction to our particular case , we can see , that our cs can be presented as a composition of two bricks of 3d level of hierarcy : @xmath35 with subsequent expansion , in accordance with composition rules ( [ cr1]),([cr2 ] ) , into bricks of 2nd level : @xmath36 which , in turn , are expanded into bricks of 1st level : @xmath37 that , at last , are expanded into basic bricks of zero level : @xmath38 in this way it is easy to calculate , that the chain can be cut into pieces either as @xmath39 or @xmath40 or @xmath41 or @xmath42 , see ( [ xl3 ] ) . mutual static disturbance of bricks is small and decreases exponentially with a level of hierarcy , therefore we have at any level of hierarcy a sequence in some order of two species of almost identical structures ( bricks ) . now let us make some important remark . in our study in @xcite of _ static _ structure of the chain , the glue particles play a some passive role only , which results in their specification as some dummy `` glue '' . however , in _ dynamical _ problem of motion particles in the chain , that we now address to , the role of glue particles becomes more important , since they have masses . actually , in this case one should consider as repeating structures _ elementary cells _ , which we get if attach to each brick ( e. g. at the left side ) one glue particle . at the basic level of hierarcy one have 55 elementary cells of 3 particles and 89 elementary cells of 5 particles , let them be called as cells @xmath43 and @xmath44 , respectively . each elementary cell , being isolated ( and , e.g. periodically closed ) has its own eigenfrequencies : three for a cell @xmath45 and five for a cell @xmath46 . then , if these cells in the chain be really isolated , then , see fig.[frgs ] we get three bands each of 55 degenerated states ( solid lines ) which belongs to cells @xmath45 , and five bands each of 89 degenerated states ( dotted lines ) . due to cells interactions all degenerated frequencies are really splitted into finite width bands , see fig.[frgs ] . the closer inspection of frequency spectrum shows , that seven largest breaks cut the exact spectrum in bands @xmath47 , @xmath48 , @xmath49 , @xmath50 , @xmath51 , @xmath52 , @xmath53 , @xmath54 so the band widths have the same order 89 , 55 , 89 , 89 , 55 , 89 , 55 , 89 , exactly as the sequence of solid and dotted horizontal lines in fig.[frgs ] ! now , let us remind , that the sequence of cells @xmath45 and @xmath46 in our gs is not random but they belong to 34 cells of the next hierarcy level , see ( [ xs1 ] ) : 13 cells @xmath55 and 21 cells @xmath56 . again , if this cells would be decoupled , we could see bands , that consist of 13 and 21 degenerated states . indeed , in the middle part of the fig.[frgs ] , where the central band is shown with greater resolution , we see a sequence of bands , that contain 21,13 , 21 , ( 21 + 13 ) , 21 , 13 , 21 number of states , which is similar to that we have at the main level of hierarcy ! however , at higher resolution ( the right box in fig.[frgs ] ) the spectrum in central band becomes structureless . this is not surprising , since the `` cells '' we have introduced are not weakly coupled objects with respect to phonon modes . in fact , what cells provides , is that some chain fragments have fixed periods ; these fragments form locally their band structure . if bands of different fragments overlap , their levels are collectivised into one common band . the central band is the case , however in this case band levels are less sensitive to small variations introduced by extra regularities of next levels of hierarcy of spatial chain structure . however , a quite new interesting phenomenon occurs , if one consider band , which belongs to a cell @xmath46 but does _ not _ overlap with any band of the cell @xmath45 . phonons with frequencies inside this band will be damped along the cells of the kind @xmath45 , these cells will play a role of some `` potential barriers '' , which decouple cells of the kind @xmath46 each from other . our chain contains 89 cells @xmath46 and 55 cells @xmath45 , which is like to 89 particles separated by 55 barriers , or distributed among 55 wells . now we have got a new effective frenkel - kontorova chain , where cells @xmath46 play role some particles with potential barriers @xmath45 among them . for a clarity sake , let us take more graphical notations for these effective `` barrier '' @xmath57 , and `` particle '' @xmath58 . then , seconary bricks , which occur in this effective fk chain , are @xmath59 and @xmath60 , and corresponding elementary cells can be obtained adding to bricks one `` glue '' particle at the left : @xmath61 , @xmath62 , having again three and five eigenfrequencies , respectively . our effective chain has 8 cells @xmath63 and 13 cells @xmath64 . now , if in the frequency spectrum of new _ effective _ fk chain we take a band , which belongs to the cell @xmath64 but not to the cell @xmath63 , it must contain 13 states . we can repeat the whole procedure ones more ; then we come to next level effective frenkel - kontorova chain with 13 effective particles @xmath64 distributed among 8 potential barriers @xmath63 , and to next generation cells , one of the kind @xmath65 and two of the kind @xmath66 . numerical data presented in fig.[frgs4 ] confirm our picture in all the detailes . here we consider the fourth band of spectra , which is well resolved from other bands , see , the left box of fig.[frgs4 ] . it contains 89 states , the number of cells @xmath46 . in the middle box we take again the fourth band which in turn contains 13 states , the number of cells @xmath64 . at last , in the right box we show by open circles 3 states of the cell @xmath65 and by closed ones @xmath67 states of two cells @xmath66 . note , that this new kind of hierarcy is quite different from that we found in @xcite with respect to spatial structure , since the transformation rules between levels of hierarcy are more complicated . in conclusion of this section it should be stressed , that all the universal features of the global band structure mentioned at the beginning are goverened by nearest order in the chain . in contrast , the fine structure depends crucially on the far order .. the latter is destroyed in the ces , therefore in ces the fine structure is washed out , and frequency spectra become smoother , see fig.[frgs_ex ] . in particular , we see also in fig.[frgs_ex ] how two lowest bands are merged into common one . the ordering of bricks , which persists in the gs , is gradually destroyed with configurational excitation of the chain . the lowest excitations are equivalent destruction of largest bricks due to permutations of bricks of preceeding level of hierarcy : the smaller bricks are destroyed , the higher is the excitation energy @xcite . the number of almost degenerate ces , originating from different bricks permutation is combunatorically large , therefore some arbitrary chosen ces has rather _ random _ sequence of bricks , and is similar to disordered media , e.g. spin glass system .    from the very beginning it is clear , that perfect randomization of bricks order may result in anderson localization of phonon modes . a new interesting point is , that fk model does _ not _ require any external disorder : randomization of bricks order in fk chain occurs dynamically . another interesting feature , which characterizes a degree of chaotization of bricks in ces , is that all the examples of exponential localization are obtained from _ single _ arbitrary taken ces , without averaging over any encemble of nearest ces . a typical quantity used traditionally in studies of localization phenomena ( see , in particular @xcite ) is the participation ratio ( pr ) , defined as @xmath68 where @xmath69 is a normalized ( @xmath70 ) vector of the phonon eigenstate . its value correspond to a chain fraction occupied by the localized state , but whether the state is localized exponentially can be unclear . typical feature of exponentially localized state is that components outside the center of localization are exponentially small . meanwhile , this components do not contribute in pr at all .    to our opinion a better characteristics to indicate an _ exponential _ localization can be the generalized mean geometrical value ( mgv ) , defined as@xmath71 where , a normalization factor @xmath16 provides that for extended states @xmath72 be order of unity . the reason in favor to @xmath72 is that it essentially better probes the exponentially small components ( tails ) of the phonon eigenstate . in particular , for a typical exponentially localized state @xmath73 , @xmath74 is a localization length , and value of the mgv : @xmath75 , i.e. becomes _ exponentially _ small . moreover , one can get from the value of mgv an estimate of the localization length as @xmath76 note an important difference between the estimate of the localization length as inverse participation ratio and our estimate ( [ lloc ] ) . the former estimates a size of the domain where the eigenstate is localized , and is sensitive to particular short range dynamics for formation of the given state . on the contrary , the latter is related to the rate of the exponential falloff at the tails of the eigenstates , and characterizes properties of the disordered media . now let us turn to localization properties of phonons in the fk chain . as in the previous section , we concentrate on the numerical study of the chain with the rotation number @xmath77 , and parameter @xmath25 . in fig.[fgs_pr - mgv ] we present our comparison of pr and mgv for gs . one can see , that both quantities @xmath78 and @xmath72 look very similar . note , that @xmath72 for gs is well distinct from zero ( dotted line ) . this means , that all phonon states are _ not _ exponentially localized ; the smallness of @xmath78 and @xmath72 for some phonon modes means only that this modes are prelocalized only , which agree with earlier studies @xcite .    in order to see , how the exponential localization manifest itself in pr and mgv , let us address to ces . we start with a typical ces , which is still in energy very close to gs : @xmath79 . despite to , that energy splitting is small , this ces belong to 5th band in the energy structure of ces @xcite , and it is one of @xmath28 cess with @xmath80 . in fig.[eloc ] we plot @xmath78 and @xmath72 for the same ces . it is seen that the behaviour of pr remains in main details the same , while in the mgv plot there are points , where the curve touches a _ zero _ line . this means , that @xmath72 at these points is exponentially small , i.e. there is an _ exponential _ localization of these modes . note , that the pr plot at the same points has no clear indications that these modes are localized exponentially . some typical examples of nonlocalized mode @xmath81 and localized one @xmath82 for this ces is presented in fig.[exmods ] . we see , that the latter is perfectly localized exponentially , with a localization length @xmath83 , which characterizes a correlation length of disorder . for comparison we show also the same mode in the gs in a `` prelocalized '' state . in fact , ces with @xmath84 corresponds an early localization of phonon modes : only small fraction of the is exponentially localized , as seen from fig.[eloc ] . to get insight , what modes are localized at the first turn , in fig.[wfr ] we plot mgv as a function of phonon frequency . we see , that localized modes are located at edges of the frequency spectrum bands . now , how this picture localization looks for ces with higher splitting from gs ? we expect @xcite , that at the splitting @xmath84 the largest robust elementary cells are @xmath85 and @xmath86 , that take part in mutual permutations only , while larger cells are destroyed . in the range of splitting @xmath87 the cell @xmath86 can dissociate @xcite , that increase a number of smaller cells @xmath85 and decrease a corelation length of the disorder . next , at splitting @xmath88 the cell @xmath85 can dissociate too , and permutations of cells @xmath89 come into play , that decrease a correlation length of the disorder even more . in apparent agreement with our expectations decrease of the disorder scale results in a total exponential localization of substantial fraction of phonon modes , especially in high frequency regoin , as seen from fig.[exw ] .    in fig.[loc2ex ] we plot fo the same two ces our estimate for localization length ( [ lloc ] ) . two dotted lines show levels @xmath90 and @xmath85 , which correspond to expected sizes of largest robust elementary cells , which survive at energy splitting @xmath91 and @xmath92 , respectively . we see , that minimal localization length follows the size of maximal robust structure in the chain . since the chain is not homogeneous , localization properties are not homogeneous too . as seen from figs . [ exw],[loc2ex ] , the localization is maximal at the high frequency part of phonon spectrum , while low frequency part seems nonlocalized . in fact , at least longwave modes of phonons show a clear tendency to be localized too , see fig.[lwstates ] . in this paper we have studied properties of phonon modes in the frenkel - kontorova chain , taken in the regime of pinned phase . spatial `` brick '' structure of the ground state of the chain found in @xcite , appears to be very useful for understanding the fine structure of the chain phonon spectrum . actually , similar analysis can also be performed for electronic spectrum , studied in recent work @xcite .    however , it is obvious , that in a real physical scale the ground state is highly degenerated with a huge number of the static configurational excitations states ( ces ) of the chain . actually this means that the true ground state of the chain is practically inaccessible . ces has properties quite different from that of the ground state : their spatial structure is rather chaotic @xcite . as a result , they can cause the anderson - like ( exponential ) localization of phonon modes similar to that seen in disordered media . this means , that results of previous studies of phonon@xcite and electron properties@xcite should be revized or extended to more realistic states of the chain . configurational excitation states ( ces ) can provide a possibility to study a gradual transition from the order to disorder . it is important , that the number of ces grows with energy splitting from the gs very fast . even at very small splitting the number of ecs is huge . it is curious that each particular ces has intrinsic chaos , which reminds the situation in classical spin glass @xcite , but in contrast to the latter , the chaos in ecs arises dynamically , without any external noise . the quantization of fk model in small @xmath93 limit is in essence the quantization of its phonon modes @xcite . therefore localization of phonon modes means localization of quantum states . note also , that phonon quantization problem is very close to electron quantization problem @xcite , where transition from localization to delocalized state is interesting as a insulator - metal transition . this work was supported in part by ec rtn network contract hprn - ct-2000 - 0156 . one of us ( o.v.z . ) thanks cariplo fundation , infn and rfbr grant no.04 - 02 - 16570 for financial support . support from the pa infn `` quantum transport and classical chaos '' is gratefully acknowledged . ziman , _ models of disorder _ ( cambridge university press , cambridge,1979 ) . zhirov , g.  casati and d.l.shepelyansky , phys . e * 65 * , 026220 ( 2002 ) . y.i . frenkel and t.k . kotorova , zh.eksp.teor.fiz.*8 * , 1340 ( 1938 ) . zhirov , g.  casati and d.l . shepelyansky , phys . e * 67 * , 056209 ( 2003 ) . f. nabarro , _ theory of crystall dislocations _ ( clarendon , oxford , 1967 ) . pokrovsky and a.l . talapov , _ theory of incommesurate crystals _ , soviet scientifical reviews supplement series physics , vol.1 ( harwood , london , 1984 ) . ying , phys . rev . b*3 * , 4160 ( 1971 ) . l.  pietronio , w.r . schneider and s.  strssler , phys . b * 24 * , 2187(1981 ) . s.  aubry , in _ solitons and condensed matter physics _ , ed . by bishop  a.r . and schneider  t. , ( springer , n.y . , 1978 ) . s.  aubry , j. phys . ( france ) , * 44 * , 147 ( 1983 ) . l.m .  floria and j.j . mazo , adv . * 45 * , 505 ( 1996 ) . m.  weiser and f.j . elmer , phys . b*53 * , 7539 ( 1996 ) ; z.  phys . b**104 * * , 55 ( 1997 ) . braun , t.  dauxois , m.v . paliy , and m.  peyrard , phys . * 78 * , 1295 ( 1997 ) ; phys . e , * 55 * , 3598 ( 1997 ) . l.  consoli , h.j.f . knops , and a.  fasolino , phys . * 85 * , 302 ( 2000 ) . s.  aubry , physica  d , * 7 * , 240 ( 1983 ) ; s.  aubry and p.y . le  daeron , physica  d * 8 * , 381 ( 1983 ) . s.  aubry and g.  abramovici , physica  d * 43 * , 199 ( 1990 ) . e.  burkov , b.e.c . koltenbach , and l.w . bruch , phys . b * 53 * , 14179 ( 1996 ) . b.  hu , b.  li , and h.  zhao , phys . e * 61 * , 3828 ( 2000 ) . ya . g.  sinai,_introduction to ergodic theory _ , ( princeton university press , 1977 ) . chirikov , phys . rep . * 52 * , 263 ( 1979 ) . beyeler , l.  pietronero and s.  strssler , phys . b * 22 * , 2988(1980 ) . f.  vallet , r.  schilling and s.  aubry , europhys . , * 2 * , 815 ( 1986 ) ; j. phys . c * 21 * , 67 ( 1988 ) . p.  tong , b.  li , and b.  hu , phys . b * 59 * , 8639 ( 1999 ) . b.  hu , b.  li , and h.  zhao , europhys . * 53 * 342 ( 2001 ) . ketoja and i.i . satija , physica  d * 104 * , 239 ( 1997 ) . m.  peyrard and s.  aubry , j.  phys . c * 16 * , 1593 ( 1983 ) . greene , j. math . * 20 * , 6 ( 1979 ) ; + r.s . mackay , physica  d * 7 * , 283 ( 1983 ) . m.  mzard , g.  parizi , and m.a . virasoro,_spin glass theory and beyond _ , ( wold sci . , singapore , 1997 ) . p.  tong , b.  li and b.  hu , phys . . lett . * 88 * , 046804 ( 2002 ) .
Abstract:

Answer: we study numerically phonon modes of the classical one - dimensional frenkel - kontorova chain , in the regime of pinned phase characterized by the phonon gap and devil s staircase , as well as by a large number of states ( configurational excitations ) , which energy splitting from the ground state is exponentially small . we demonstrate , these states behave like disorder media : their phonon modes are _ exponentially _ localized , in contrast to the phonon modes in the ground state , where phonons are _ prelocalized _ only . we demonstrate also , the phonon frequency spectrum of the ground state has an hierarchical structure , a direct manifestation of hierarchical spatial structure , found for the ground state of the fk chain in our recent work .


Instruction: Write an abstract for the following acadamic paper.

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

Paper: the _ swift _ mission , with its ability to localize gamma - ray bursts ( grbs ) in real time , has revolutionized our understanding of these phenomena in many ways . one of the most interesting discoveries is that the light curve of the x - ray afterglow displays a large diversity of behaviors ( nousek et al . 2006 ) , rather than being a relatively featureless power - law .    the x - ray afterglow sets in as a rapidly fading source at the end of the prompt emission ( tagliaferri et al . this early phase is understood as the radiation of the prompt phase reaching the observer from off - axis angles @xmath1 ( kumar & panaitescu 2000 ; kumar et al . 2006 ; lazzati & begelman 2006 ; zhang et al . 2007 ) . the steep decay phase is usually followed by a flat component , whose origin is still highly debated ( granot & kumar 2006 ; uhm & beloborodov 2007 ; genet et al . 2007 ; ghisellini et al . finally , at random times between a few hundred seconds up to several tens of thousands of seconds after the onset , the x - ray afterglow displays sudden rebrightenings , known as x - ray flares ( burrows et al . 2005 ; falcone et al . 2006 , 2007 ; chincarini et al . 2007 , hereafter c07 ) . x - ray flares could be due to a variety of causes , either related to external shock activity or to the inner engine itself . it can be shown that any mechanism related to the external shock would produce flares with a characteristically long timescale ( lazzati et al . 2002 ; lazzati & perna 2007 ) . most of the observed flares have fast rise and decay times ( c07 ; kocevski , butler & bloom 2007 ) and must therefore be related to activity of the central engine at times comparable to those at which the flare is observed . for this reason , they are of great importance for our understanding of the mechanism that powers the grb outflows . they are potentially unique laboratories to investigate the properties of relativistic outflows from compact objects over a broad range of luminosities . as we will see in the following , the isotropic equivalent luminosity of the outflow ranges from @xmath2  erg / s during the prompt phase to @xmath3  erg / s during flares at late times .    we study a sample of grbs with x - ray flares that have been observed by the swift x - ray telescope ( xrt ) for a sufficiently long time and for which redshift information is available . we compute the average energy output from the inner engine as a function of the time elapsed since the grb explosion , and we compare it to several mechanisms for energy extraction from a magnetar or a black hole . this letter is organized as follows : in   2 we discuss the sample and the procedure used to derive the cumulative light curve , in   3 we discuss different mechanisms that could power the flares , in   4 we discuss the role of the progenitor star in shaping the late time outflow , and in   5 we discuss our results , their limitations and their implications . we selected from the tables in falcone et al . ( 2007 ; hereafter f07 ) the sub - sample of bursts with x - ray flares for which redshift information is available . the bursts and some of their key properties are listed in tab .  [ tab : sample ] . the sample consists of 10 grbs for a total of 24 flares . all the bursts have xrt observations up to comoving time @xmath4  s.    the flare isotropic equivalent energy was computed from the fluence @xmath5 in tab . the average flare light curve between the times @xmath6 and @xmath7 was computed as : @xmath8 where @xmath9 is the average luminosity of the @xmath10 flare during its active time , @xmath11 is the total number of grbs considered , @xmath12 is the total number of flares and @xmath13 is the fraction of time in the interval @xmath14 during which the @xmath15 flare is active . this prescription is equivalent to the assumption that the flares have a square shape . even though this is a very poor approximation to the real shape of the flares , its effect is negligible when many flares are averaged to compute the light curve . we tested both a gaussian shape ( c07 ) and a triangular shape , and we obtained consistent results within the errors . the computation of the errors in the average light curve is nontrivial , since the main contribution is not the uncertainty in each fluence measurement , but rather the uncertainty in the number of flares that are observed on average during a given time interval . for this reason we assume that the uncertainty of each fluence is as large as the measurement itself and we compute the errors on the average light curve as : @xmath16^{1/2}.\ ] ] this prescription ensures that if the average flux at a certain time is dominated by a single bright flare , the error is large . small errors occur only in time intervals with a large number of flares . .properties of the bursts selected for the analysis [ tab : sample ] [ cols="<,^,^,^",options="header " , ]     the resulting light curve is shown in fig . [ fig : z ] , where @xmath17 is shown versus time ( @xmath18 ) . the solid dark points with error bars show the results assuming a band model spectrum for the flares ( band et al . 1993 ; f07 ) , while the white lozenges show the results of the power - law spectral model ( f07 ) . we model the average light - curve as a power - law , excluding the first point ( at @xmath19  s ) since it is still contaminated by the prompt emission in many cases . we find that the average light curve is very well described by a power - law with index ( @xmath20 errors ) : @xmath21 the fit has @xmath22d.o.f.@xmath23 , a small value that is not surprising given the very conservative assumptions on the uncertainties . the formal error of 0.16 is , as a consequence , very conservative . in order to check the result , we computed the average light curve also as the derivative of the cumulative energy produced as a function of time since the grb onset , finding analogous results . averaging many bursts together allowed us to increase the signal - to - noise ratio of the flare light curve . we now examine whether the resulting slope can be recovered in individual grbs that display a large number of flares . when we compute the flare light curve for an individual event , the knowledge of the redshift is not necessary and we can select the grbs from a larger pool of events . there are three grbs in the f07 catalog with six or more flares . one is grb  050904 , which is in our sample and has 7 flares . another is grb  051117a , which has 7 flares according to f07 . however , the seven flares of grb  051117a overlap one another and therefore we discard this event from our analysis . finally , grb  050803 has six flares . figure  [ fig:050803 ] shows the flare light curve , in observed time and flux , for these two bursts . power - law fits were performed for these cases , yielding @xmath24 and @xmath25 for grb  050803 and grb  050904 , respectively . even though the uncertainties are larger , as expected given the smaller statistics , the individual cases confirm that flares produce more luminosity at early times than at late times , their light curve is consistent with a power - law , and the index of the power - law is consistent with the cumulative curve of fig . [ fig : z ] .    grb  050904 has seven flares and is also part of our main sample . since it has so many flares , we checked that the results for the average flare energy are not dependent on the presence of grb  050904 in the sample . reproducing fig . [ fig : z ] without the seven flares of grb  050904 yields similar results , with a slope @xmath26 , in agreement with the result of eq . [ eq : slope ] . even though collimated outflows originating from compact objects accreting matter from a disk are ubiquitous in the universe ( e.g. , ferrari 1998 ; blandford 2001 ) , the actual mechanism that produces the jet and determines its properties ( luminosity , entropy , opening angle , structure , and magnetization ) remains elusive . in addition , we see jets in other types of objects , such as pulsars , where mechanisms other than accretion could play a role in the jet production ( but see blackman & perna 2004 ) . in the case of grbs , several mechanisms have been proposed to produce the relativistic jet : conversion of internal energy into bulk motion with hydrodynamic collimation ( cavallo & rees 1978 ; lazzati & begelman 2005 ) , energy deposition from neutrinos , energy released from a rapidly spinning , newly born magnetar ( usov 1992 ) , and magnetic collimation and acceleration ( vlahakis & knigl 2001 ) . it has proven so far extremely challenging to prune some of these possibilities , since the grb radiation is produced far away from the place where the outflow is accelerated , and we have not been able to connect the properties of the radiation ( light curves and spectra ) to the ( magneto)-hydrodynamical properties of the plasma producing it . let us first consider a system made by a newly formed stellar mass black hole accreting matter at a high rate from a disk . during the prompt grb phase , such a disk is so hot and dense that neutrino losses provide an effective cooling mechanism ( popham , woosley & fryer 1999 ; narayan , piran & kumar 2001 ; chen & beloborodov 2007 ) . in this phase , the jet may be powered by neutrino annihilation , even though magnetic effects could play a dominant role . in the latter case , the outflow may either originate on the disk surface , as material escapes along low inclination magnetic field lines ( e.g. , levinson 2006 ) , or by blandford - znajek mechanisms ( e.g. , mckinney & gammie 2004 ) . 1d steady - state calculations of chen and beloborodov ( 2007 ) showed that , for a rotating bh , the switch from a neutrino - cooled to an advective disk takes place at an accretion rate of @xmath27s , depending on the viscosity prescription used . in all cases , the switch - off of the neutrino cooling is expected to be associated to a change in the grb outflow characteristics , since the magnetic field that can be anchored to a dense , high pressure , thin disk is expected to be stronger than the one anchored to a lower - density , lower - pressure , thick inflow .    figure  [ fig : mdot ] shows the average flare luminosities in units of the accretion rate in solar masses per second , where we assumed that the grb is beamed into 1 per cent of the sky , that the efficiency of converting the accretion rate into jet luminosity is 0.1 per cent , and that these two parameters are constant . the two parameters are rather standard ( e.g. , chen & beloborodov 2007 ; podsiadlowski et al . 2004 ) . changing them by a factor up to ten would not affect the main conclusion since the prompt and flare luminosities would scale in the same way . the average luminosity of the prompt phase of the 10 grbs is also shown in the figure with a thick solid line . such accretion rates are compared to the one at which neutrino cooling is switching on / off ( chen & beloborodov 2007 ) . this transition lies suggestively at the boundary between the prompt phase and the flaring phase . we argue that this could be an indication of the fact that the prompt phase of grbs , the most luminous one , is powered by accretion onto a black hole in the form of a geometrically thin  neutrino cooled  disk . as the accretion rate drops under the critical value , the fast accretion is switched off and the prompt phase ends . the disk swells under the effect of internal pressure and a new accretion geometry , with a thick disk , sets in . this gives rise to the flares that we observe during the afterglow phase .    in the case of a magnetically driven jet , the change in the outflow dynamics would be brought about by the fact that a thick disk has a lower pressure and therefore can anchor a magnetic field with lower intensity . it would also change the pitch angle of low - inclination magnetic field lines . in the case of a neutrino powered outflow , the switch off of the neutrino cooling would shut off completely the outflow . it seems therefore that the observations of late - time , low - luminosity flares provide evidence against the powering of grb outflows exclusively by neutrino annihilation .    the material necessary to provide the accretion at late stages can be provided , in long - duration gamma - ray bursts , by the fallback of material that did not reach the escape velocity in the stellar explosion ( chevalier 1989 ; macfadyen , woosley & heger 2001 ; zhang & woosley 2008 ) . alternatively , especially in the case of short grbs , the natural evolution of an accretion disk can supply the accretion rate necessary to explain the observations . consider an accretion disk that forms in a short timescale ( comparable to or smaller than the duration of the grb prompt phase ) and is left to evolve without any sizable mass input thereafter . the accretion rate depends on the assumptions made on the nature of viscosity . the case of a thin disk has been studied thoroughly . frank , king & raine ( 2002 ) report the exact solution for the case of a thin disk with constant viscosity . they show that , after several tens of viscous time - scales , the accretion rate onto the central object approaches the asymptotic form @xmath28 . cannizzo , lee & goodman ( 1990 ) performed numerical simulations for the same initial conditions using the prescriptions of @xmath29-viscosity ( shakura & sunyaev 1973 ) . they find that the accretion rate at late times scales with time as @xmath30 , independently of the value of @xmath31 assumed . the case of a thick disk needs to be studied in detail ( lazzati & begelman in preparation ) . the comparison of the observed flare luminosity to the theoretical rates of late - time accretion implies that there is an almost linear relation between accretion and the luminosity of the outflow . in our favored scenario , where the outflow is powered by magnetic processes , this relation would be caused by the decay of the magnetic field as the disk becomes less massive and dense . the reason why we observe such a linear relation is not entirely clear and deserves further investigation .    an alternative to the accretion disk  bh system is a rapidly spinning magnetar that powers the outflow as a consequence of spin - down ( thompson , chang & quataert 2004 ; bucciantini et al . this system can provide late - time energy either through neutrino emission or through dipole radiation . the neutrino emission decays exponentially and can not provide the power required ( thompson et al . 2004 ) . in the simple vacuum dipole scenario the decay of the late - time energy deposition is too steep ( @xmath32 ) . however , the interaction of the field with the stellar material can create shallower slopes consistent with the observations . in this scenario , the fact that the neutrino cooling switch coincides with the transition from the prompt to the flaring phase is purely coincidental . in addition , it is not clear how the continuous luminosity produced by the spin - down of the magnetar can be converted into a flaring source . even though the average flare curve of fig . [ fig : z ] is featureless , we ought to keep in mind that flares are episodic events . the fact that an accretion disk or a spinning magnetar can provide energy at late time does not ensure that the energy will be released intermittently . it is not surprising to observe variability of relativistic sources on timescales comparable to the dynamical time of the system . however , in the case of flares , the variability timescale is many orders of magnitude longer than the dynamical one for a solar mass black hole , and it grows approximately linearly with time ( c07 ) . it is more likely that this variability is associated to a viscous time - scale rather than to the dynamical time scale . we must therefore seek either a mechanism that can release energy episodically from the inner engine , or a way to transform a continuous output of energy into a highly variable one as it propagates from the engine to the radiation zone . the propagation of the jet through the cold stellar progenitor material provides in principle a way by which a continuous outflow can be converted into a succession of individual fireballs . morsony et al . ( 2007 ) showed that even if a continuous jet is injected into the core of a massive star , the ensuing light curve is highly variable . consider an engine releasing an outflow with continuous but decreasing luminosity . inside the star the jet is in pressure equilibrium with a high pressure cocoon of material ( lazzati & begelman 2005 ) . as the jet luminosity decreases , the cocoon becomes over - pressured and squeezes the jet . at the same time , the cocoon pressure decreases due to the fact that the cocoon material is escaping from the surface of the star . if , at any point in time , the cocoon pressure overcomes the stagnation pressure of the jet , the jet would be choked and the flow of energy interrupted . the stagnation pressure of a relativistic outflow with lorentz factor @xmath33 and cross section @xmath34 is given by @xmath35 , where @xmath36 is the luminosity of the jet . following lazzati & begelman ( 2005 ) , the cocoon pressure can be written as : @xmath37 where @xmath38 is the average jet luminosity before breakout , @xmath39 is the average stellar progenitor density and @xmath40 its radius , and @xmath41 is the breakout time . a condition for the stagnation of the jet can be obtained by comparing the stagnation pressure with eq . [ eq : pcocoon ] . it can be seen that the cocoon pressure can not reach the jet stagnation pressure for any reasonable parameter set . we conclude therefore that the instability giving rise to the flaring behavior of the late time activity of the inner engine has to be intrinsic to the jet release process or to its transition from the non - relativistic to the relativistic stage , and can not be brought about by the propagation of the relativistic outflow in the star . we have computed the average energy released in the form of x - ray flares overlaid on the power - law decay of the afterglow of long - duration gamma - ray bursts . a sample of 10 long duration grbs from the catalog of f07 with redshift measurements was used for the analysis . we conclude that , on average , the late - time energy release approximately follows the power - law scaling @xmath42 . our analysis is possibly affected by several biases . first , the definition of a flare is fraught with uncertainty , as the comparison of the c07 and the f07 catalogs easily reveals . we here adopt the definition of f07 and refer to that paper for a discussion of their selection criteria . second , the measurement of the slope of the power - law depends on our capability to select flares . low brightness flares could be lost in the early phases when the afterglow can easily outshine them . on the other hand , very short duration flares could be missed at late times when the observations are not continuous in time . it appears that the first bias is most serious ( since no very short duration flare was ever detected in the late phases ) and that the slope is possibly underestimated . we do not believe this should affect our conclusions , since almost all the grbs we considered have an early flare and so the contribution of shallow flares would be minimal . one important effect could , however , create a systematic overestimation of the slope . the late time flares may be less beamed than the early ones , and therefore their isotropic equivalent luminosities would appear smaller than the one of early flares due , in part , to geometric end not intrinsic effects . according to simulations , however , most of the opening angle evolution takes place at very early times ( morsony et al . 2007 ) , as long as the injection opening angle does not evolve . the energy to power the flares may come from accretion or from the spinning down of a magnetar . in the first case , we were able to estimate the accretion rate required to power the prompt emission as well as the flaring phase . interestingly , during the prompt phase the accretion rate is so large that neutrino emission cools the flow and accretion takes place in the form of a thin disk . at accretion rates of about 0.001 to 0.01 solar masses per year , the neutrino luminosity drops and the accretion disk becomes thick . we find that the switch - off of the neutrino cooling , and the transition from a thin to a thick disk , take place between the prompt and the flaring phase . indeed , the prompt emission spikes and the late - time flares exhibit some differences . while the prompt spikes do not show any evolution in their duration or peak luminosity , the flares become longer and shallower with time . this property naturally arises in a disk that fragments and accretes the various blobs of material on their viscous timescales ( perna , armitage & zhang 2006 ) . in addition , flares are very episodic and the engine is  off " most of the time at late stages , while during the prompt phase the engine is  on " most of the time . we propose that the differences are due to the switch from a thin disk to a thick disk configuration . a more detailed analysis of the thick disk accretion dynamics is on - going ( lazzati & begelman in preparation ) . an alternative scenario is that of a spinning - down pulsar . in this case , the fact that the neutrino switch coincides with the transition from the prompt to the flare phase would be a pure coincidence . more observations and a better understanding of both engine models are needed before definitive conclusions can be drawn . although we have focused on long - duration grbs in this paper , from a theoretical point of view short - duration grbs may also have a flaring phase , as long as they are powered by accretion from a disk or by the spin - down of a magnetar . the late - time accretion provided by a disk that is draining into the bh is sufficient to power the flares independently of the presence or not of fall - back material from the progenitor star . we thank andrei beloborodov , brian metzger , eliot quataert and todd thompson for useful discussions . this work was supported by nsf grants ast-0307502 and ast-0507571 , nasa atp grant nng06gi06 g , and swift gi program nnx06ab69 g . 99 band d. , et al . , 1993 , apj , 413 , 281 blackman , e. g. , perna , r. 2004 , apj , 601 , l71 blandford r.  d. , 2001 , progress of theoretical physics supplement , 143 , 182 bucciantini n. , quataert e. , arons j. , metzger b.  d. , thompson t.  a. , 2008 , mnras , 383 , l25 burrows d.  n. , et al . , 2005 , sci , 309 , 1833 cannizzo j.  k. , lee h.  m. , goodman j. , 1990 , apj , 351 , 38 cavallo g. , rees m.  j. , 1978 , mnras , 183 , 359 chen w .- x . , beloborodov a.  m. , 2007 , apj , 657 , 383 chevalier r.  a. , 1989 , apj , 346 , 847 chincarini g. , et al . , 2007 , apj , 671 , 1903 ( c07 ) falcone a.  d. , et al . , 2006 , apj , 641 , 1010 falcone a.  d. , et al . , 2007 , apj , 671 , 1921 ( f07 ) ferrari a. , 1998 , ara&a , 36 , 539 frank j. , king a. , raine d.  j. , 2002 , `` accretion power in astrophysics '' , cambridge university press , cambridge , uk genet f. , daigne f. , mochkovitch r. , 2007 , mnras , 381 , 732 ghisellini g. , ghirlanda g. , nava l. , firmani c. , 2007 , apj , 658 , l75 granot j. , kumar p. , 2006 , mnras , 366 , l13 kocevski d. , butler n. , bloom j.  s. , 2007 , apj , 667 , 1024 kumar p. , panaitescu a. , 2000 , apj , 541 , l51 kumar p. , mcmahon e. , barthelmy s.  d. , burrows d. , gehrels n. , goad m. , nousek j. , tagliaferri g. , 2006 , mnras , 367 , l52 lazzati d. , rossi e. , covino s. , ghisellini g. , malesani d. , 2002 , a&a , 396 , l51 lazzati d. , begelman m.  c. , 2005 , apj , 629 , 903 lazzati d. , begelman m.  c. , 2006 , apj , 641 , 972 lazzati d. , perna r. , 2007 , mnras , 375 , l46 levinson a. , 2006 , apj , 648 , 510 macfadyen a.  i. , woosley s.  e. , heger a. , 2001 , apj , 550 , 410 mckinney j.  c. , gammie c.  f. , 2004 , apj , 611 , 977 morsony b.  j. , lazzati d. , begelman m.  c. , 2007 , apj , 665 , 569 nousek j.  a. , et al . , 2006 , apj , 642 , 389 perna , r. , armitage , p. j. , zhang , b , 2006 , apj , 636 , l29 podsiadlowski p. , mazzali p.  a. , nomoto k. , lazzati d. , cappellaro e. , 2004 , apj , 607 , l17 popham r. , woosley s.  e. , fryer c. , 1999 , apj , 518 , 356 shakura n.  i. , syunyaev r.  a. , 1973 , a&a , 24 , 337 tagliaferri g. , et al . , 2005 , nature , 436 , 985 thompson t.  a. , chang p. , quataert e. , 2004 , apj , 611 , 380 uhm z.  l. , beloborodov a.  m. , 2007 , apj , 665 , l93 usov v.  v. , 1992 , nature , 357 , 472 vlahakis n. , knigl a. , 2001 , apj , 563 , l129 zhang b .- b . , liang e .- w . , zhang b. , 2007 , apj , 666 , 1002 zhang w. , woosley s. e. , heger a. , 2008 , apj in press ( arxiv : astro - ph/0701083v2 )
Abstract:

Answer: we compute the average luminosity of x - ray flares as a function of time , for a sample of 10 long - duration gamma - ray burst afterglows . the mean luminosity , averaged over a timescale longer than the duration of the individual flares , declines as a power - law in time with index @xmath0 . we elaborate on the properties of the central engine that can produce such a decline . assuming that the engine is an accreting compact object , and for a standard conversion factor between accretion rate and jet luminosity , the switch between a neutrino - cooled thin disk and a non - cooled thick disk takes place at the transition from the prompt to the flaring phase . we discuss the implications of this coincidence under different scenarios for the powering of the grb outflow . we also show that the interaction of the outflow with the envelope of the progenitor star can not produce flares out of a continuous relativistic flow , and conclude that it is the dynamics of the disk or the jet - launching mechanism that generates an intrinsically unsteady outflow on timescales much longer than the dynamical timescale of the system . this is consistent with the fact that x - ray flares are observed in short - duration grbs as well as in long - duration ones . gamma - ray : bursts