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: frictional interfaces are abundant in biological ( e.g. adherent cells and cell locomotion ) , engineering ( e.g. micro - electro - mechanical devices ) and geophysical ( e.g. earthquake faults ) systems around us , and are of fundamental and practical importance . consequently , understanding the dynamics of dry frictional interfaces has been the focus of intense scientific activity in the last few decades @xcite . it has been established that under steady - state sliding conditions , the frictional resistance features a non - trivial velocity dependence , and that this dependence has dramatic consequences on the dynamic response of frictional interfaces @xcite . specifically , it has been shown that for a broad range of materials friction is velocity - weakening  that is , the steady frictional resistance is a decreasing function of the sliding velocity  at least in the regime of low velocities , up to a few hundreds of microns per second . this feature favors various instabilities and stick - slip motion @xcite . a very recent compilation of a large set of experimental data for a broad range of materials , however , has revealed that for higher slip velocities , friction generically becomes velocity - strengthening over some range of slip velocities @xcite . the existence of velocity - strengthening behavior might have significant effects on various aspects of frictional dynamics . in particular , recent laboratory experiments on fault - zone materials have documented slow slip interfacial events  an intensely debated issue  and have linked it to a crossover in the frictional response , from velocity - weakening to velocity - strengthening friction , with increasing slip velocity @xcite . while the possible implications of the existence of velocity - strengthening friction have been rather sporadically discussed in the literature @xcite , to the best of our knowledge a comprehensive and systematic theoretical exploration of these important issues is currently missing .    as a first step in closing this gap , we study here the effect of velocity - strengthening friction on spatiotemporal interfacial dynamics , energy dissipation and radiation , and the global interfacial strength , with a special focus on the nucleation , propagation and arrest of rupture fronts . we explore three variants of a realistic rate - and - state friction law , one which is purely velocity - weakening , one which crosses over at higher velocities to logarithmic velocity - strengthening friction , and one which crosses over to linear velocity - strengthening friction . we show that the presence or absence of velocity - strengthening friction at relatively high slip velocities can significantly affect the global interfacial resistance ( strength ) and the energy released during frictional instabilities ( `` event magnitude '' ) , even under quasi - static loading conditions . different forms of velocity - strengthening friction , in our case logarithmic and linear , give rise to events of similar magnitude , yet with dramatically different dissipation and radiation rates . the difference stems from the broad range of the underlying rupture propagation velocities , where stronger velocity - strengthening friction promotes slower rupture , possibly orders of magnitude slower than elastic wave - speeds . this result is directly related to the recent experimental observations of @xcite . all in all , our results show that velocity - strengthening friction should be properly quantified and incorporated into friction theory as it appears to affect many basic properties of spatially extended frictional interfaces . the ideas to be presented below have been originally influenced by the works in @xcite and then further developed in @xcite . the rate - and - state friction model we study has been introduced recently in @xcite , and is reviewed here briefly . we start by considering a multi - contact interface and write the ratio @xmath0 of the real contact area ( the area of all contact asperities ) to the nominal one , in terms of a state parameter @xmath1 ( of time dimensions ) as @xmath2 \ ,   \label{eq : a}\ ] ] where @xmath3 is the normal ( compressive ) stress at the interface , @xmath4 is the material hardness , @xmath5 is a dimensionless material parameter of order @xmath6 , and @xmath7 is a short time cutoff @xcite . @xmath1 is usually interpreted as the interface s effective age , and its evolution is given by @xmath8 where @xmath9 is the local interfacial slip velocity and @xmath10 is a lengthscale related to the contact asperities geometry . @xmath11 , with an extremely small @xmath12nm / s , is a regularization function that plays no important role , and is actually omitted in all of the analytic results that follow . as we focus here on unidirectional motion , we do not distinguish between @xmath9 and @xmath13 . the frictional stress @xmath14 is written as a sum of an elastic contribution , @xmath15 , and a viscous contribution @xmath16 , @xmath17 the viscous contribution takes the form @xmath18 , where at least at low velocities , the rheological part @xmath19 corresponds to a stress - biased thermally - activated process @xcite @xmath20 here , @xmath21 is boltzmann s constant , @xmath22 is the absolute temperature , @xmath23 is an activation volume and @xmath24 is a velocity scale , related to a microscopic attempt rate . a higher velocity variant of eq . ( [ eq : tauvis ] ) will be discussed below . the elastic stress follows the evolution equation @xmath25 where @xmath26 is the interfacial shear modulus and @xmath27 is the effective height of the interface . equations - describe the first variant of the friction model we study below . we begin by describing its behavior under steady sliding at a velocity @xmath28 . the steady solution of eq . is @xmath29 , from which it follows that the contact area is a logarithmically decreasing function of @xmath9 @xcite . the fixed point of eq . reads @xmath30 , and hence the overall frictional resistance is given by @xmath31 where a higher order logarithmic term was omitted and the following definitions were used @xmath32 in the low velocity regime , i.e. @xmath33 , @xmath34 is a logarithmic function of @xmath9 , with @xmath35 . therefore , if @xmath36 ( which is quite generically the case ) , friction is logarithmic velocity - weakening .    physically , friction is velocity - weakening because the real contact area is a decreasing function of the sliding velocity , and its velocity dependence is stronger than the rheological dependence of @xmath16 . however , as discussed at length in @xcite , when @xmath37 the contact area saturates , and friction becomes logarithmically velocity - strengthening . we term this model the _ logarithmic velocity - strengthening _ ( ls ) friction model . the resulting steady - state friction curve is shown in fig . [ fig : ss ] .     vs. the slip velocity @xmath9 for the three model variants ( different colors , also marked by labels ) . note that all of the curves coincide at low velocities and that the driving velocity @xmath28 is marked . ]    in case the contact area continues to decrease indefinitely with increasing @xmath9 , friction remains velocity - weakening for arbitrarily high velocities . this is formally achieved by removing the `` 1 '' in the argument of the logarithm in eq . , that is , replacing eq . by @xmath38 \ . \label{eq : a_wrong}\ ] ] although this is somewhat unphysical , this choice was widely used in the literature @xcite , and we term it the _ pure velocity - weakening _ ( pw ) friction model . the resulting steady - state friction curve is also shown in fig . [ fig : ss ] . a third variant of the model is obtained by modifying the rheological function @xmath19 , cf . as discussed extensively in @xcite , and to some extent in @xcite , the simple picture of a single barrier , linearly biased , thermally - activated process is expected to break down when asperity - level stresses become sufficiently large . when this happens , a different dissipation mechanism is expected to dominate friction . while at the moment there is no general quantitative theory for the velocity dependence of friction in this regime , it is not expected to be logarithmic , but rather to exhibit a significantly stronger dependence on the slip velocity . in this work , we consider a simple model in which the logarithmic dependence crosses over continuously ( but not smoothly ) to a linear viscous rheology . explicitly , we replace eq . by @xmath39 & v > v_c\\                   \end{cases}\ , \label{eq : tau_lts}\ ] ] where @xmath40 is a dimensionless parameter . we term this model the _ stronger - than - logarithmic _ ( stl ) velocity - strengthening model . the resulting steady - state friction curve is shown in fig . [ fig : ss ] . we stress that all three variants coincide in the low velocity regime , where they feature logarithmic velocity - weakening friction . at higher slip velocities , the ls variant , which is described by eqs . - , features a crossover to logarithmic velocity - strengthening friction . the pw variant does not feature any strengthening at all ( i.e. it remains velocity - weakening ) , and is obtained from the ls model by using eq . instead of . the stl variant features linear velocity - strengthening friction , and is obtained from the ls model by using eq . instead of eq . .    in order to investigate the implications of the different constitutive laws on frictional dynamics , we need to consider a spatially - extended interface under inhomogeneous sliding conditions . to this end , we consider a long elastic block of height @xmath41 ( in the @xmath42-direction ) and length @xmath43 ( in the @xmath44-direction ) , in frictional contact ( at @xmath45 ) with a rigid substrate ( i.e. no deformation of the substrate is considered ) , see fig . [ fig : sketch ] . the trailing edge of the elastic block ( at @xmath46 ) is moved at a constant velocity @xmath28 in the positive @xmath44-direction , while the leading edge ( at @xmath47 ) is stress - free . the block is driven quasi - statically with @xmath48 m/s , which is representative of typical laboratory experiments @xcite and generically belongs to the steady - state velocity - weakening friction branch ( cf . [ fig : ss ] ) . the upper edge of the elastic block ( at @xmath49 ) experiences a constant normal stress @xmath3 , @xmath50 , but no shear stress , i.e. @xmath51 .     ( @xmath41 is the block s height ) and driven by a velocity @xmath28 at its trailing edge ( @xmath46 ) . the leading edge is at @xmath47 . the shear stress at the interface , @xmath52 , equals to the frictional stress @xmath53 . ] we focus on plane - strain deformation conditions and furthermore assume that @xmath41 is smaller than the smallest lengthscale @xmath54 characterizing the spatial variation of various fields in the @xmath44-direction . under the stated conditions , the momentum balance equation @xmath55 where @xmath56 and @xmath57 ( @xmath58 ) are the components of the displacement vector and cauchy s stress tensor , respectively , and @xmath59 is the mass density , reduces to ( see @xcite for derivation ) @xmath60 where the plane - strain hooke s law was used . here @xmath62 ( where @xmath63 is the shear modulus of the bulk and @xmath64 is poisson s ratio ) and the shear stress at @xmath45 simply equals the frictional stress , @xmath65 . note also that @xmath66 . corrections to eqs . ( [ eq : eom1])-([eq : eom2 ] ) appear only to order @xmath67 , a situation reminiscent of the shallow water approximation in fluid mechanics . finally , note that the lateral force required to maintain the velocity boundary condition at the trailing edge , @xmath68 , reads @xmath69 and the traction - free boundary condition at the leading edge implies @xmath70 . equation , with the stated boundary conditions and with @xmath53 corresponding to one of the three friction laws described above , has been solved numerically using a commercial differential equations solver . the model parameters for polymethyl - methacrylate(pmma ) , a polymeric glass that is widely used in laboratory experiments @xcite , were extracted from a large set of experimental data . the parameters are listed in table [ tab : params ] , and the procedure for obtaining them is described in @xcite . the initial conditions are @xmath71 , @xmath72 , @xmath73 , and @xmath74s , the latter is typical of laboratory scale experiments . the results presented here are largely insensitive to the choice of the initial value of @xmath1 . vs. time for the three models ( color code as in fig . [ fig : ss ] ) . it is seen that all of the curves coincide for short times , and then begin to diverge . the ls and stl models maintain the same `` envelope '' , while the pw model features more pronounced stress drops , larger inter - event time and a lower overall resistance . ( inset ) the same data as in the main panel , but this time @xmath75 is plotted vs. time . the red lines are linear fits to the values of @xmath75 at the rupture arrest times @xmath76 ( i.e. @xmath75 right after the force drops ) , cf . the prediction in eq . . ]    we begin by studying the macroscopic response of the system . figure [ fig : fd ] shows the total frictional force exerted by the loading machine as a function of time , @xmath77 . it is seen that the friction force increases gradually until it experiences an abrupt drop , followed by repeated cycles of gradual increases and abrupt drops , typical of frictional systems @xcite . the drops in the friction force , which appear as vertical lines in this figure , occur when sliding becomes unstable , and involve nucleation and propagation of rupture fronts , as will be discussed below .    before the first drop , the friction force corresponding to the three variants is identical , as can be expected because the dynamics in this regime are slow and governed by the loading velocity @xmath28 . in this range of velocities , the three variants coincide and consequently the first drop occurs almost exactly at the same point in time for all of the variants , suggesting that the instability mechanism is insensitive to the high velocity behavior ( as predicted in @xcite ) . however , since the instabilities are accompanied by much larger velocities , the high velocity behavior of the friction law becomes important . figure [ fig : fd ] demonstrates that while the ls and stl models give rise to almost identical force profiles , the pw model results in significantly larger force drops , and a lower overall interfacial resistance . this suggests , and will be further substantiated in what follows , that while the total energy dissipated during these drops is similar in the ls and stl models , the energy dissipated in the pw model is significantly larger . other features of the global friction curves shown in fig . [ fig : fd ] , such as the lower envelope of @xmath77 ( corresponding to the values of @xmath77 after each drop ) , will be discussed and explained theoretically below . is the spatial location of the front tip , cf . [ fig : before_after ] . @xmath78 is the time where the front starts to propagate , cf . [ fig : fd ] . note also the vast difference in timescales between the panels . the wave speed in this system is @xmath79 m / s . ]    in order to understand the origin of these differences , one must examine the complex spatiotemporal dynamics that give rise to the `` force drop events '' , which are described at length in @xcite . as stated above , the instabilities result in the nucleation , propagation and arrest of rupture fronts , a scenario reported by many experimental , numerical and analytical works @xcite . most of the remainder of this paper will be focused on the first event , which is marked in fig . [ fig : fd ] by @xmath78 . the rationale for focussing on the first event ( rather than some later event ) is that it ensures that the state of the interface is the same for all three model variants at the onset of instability ( with no history effects ) , cleanly isolating the effects of the existence and form of the velocity - strengthening branch . having said that , we note that it is clear from fig . [ fig : fd ] that the differences between the three variants persist to _ any _ event . furthermore , multiple - event properties will be explicitly discussed in relation to eqs . - and the inset of fig . [ fig : fd ] .    figure [ fig : fronts ] shows the propagation and arrest of rupture fronts during the first event . first , we note the vast difference in the timescales involved : while rupture fronts in the ls and pw models arrest after a few 10@xmath80s , in the stl model they last for a few ms . it is observed , however , that while the penetration depth of the front into the interface in the ls and stl models is comparable , for the pw model it is an order of magnitude larger . furthermore , the rupture propagation velocity in the ls model is an order of magnitude smaller than in the pw model ( the latter is of the order of the elastic wave - speed ) , and the propagation velocity in the stl model is yet two orders of magnitude smaller .    both the ls and stl models give rise to rupture fronts that are much slower than the elastic wave - speed . these remarkably low rupture propagation velocities , three orders of magnitude slower than the elastic wave - speed in the stl model , might be related to the important , and rather intensely debated , issue of slow rupture @xcite . our calculations suggest that the emergence of slow rupture might be directly related to the existence and form of velocity - strengthening friction . this is in accord with recent laboratory experiments on fault - zone materials , which documented slow slip events together with a clear crossover from velocity - weakening to velocity - strengthening friction with increasing slip velocity @xcite .    . it is seen that the stress left at the tail of the rupture fronts , @xmath81 , is roughly homogeneous in space , and that it is much lower in the pw model than in the ls and stl models . the location of the fronts after the event is marked by @xmath82 . the deeper penetration of the pw model , also shown in fig . [ fig : fronts ] , is clearly visible . ]    a lot can be learned from the state of the interface after the rupture front has passed . in fig . [ fig : before_after ] we plot the spatial distribution of the ( normalized ) friction stress just before the first rupture event and immediately after it for the three variant models . in both of these states , the higher slip rates associated with the rupture fronts are not present ( before the event they have not yet been generated and after the event they have died off ) , and the mechanical state is quasi - static . in line with the previous results , prior to the inception of the first event the stress profiles in the three models essentially coincide . when the fronts propagate and eventually arrest , they leave behind them a residual stress profile , which is much smaller in the pw model compared to the ls and stl models . this residual stress is approximately homogeneous in space and is lower than the stress prior to the event . the elastic energy release during this stress relaxation process is the driving force to frictional dissipation . the approximate spatial homogeneity of @xmath14 left behind _ any _ rupture front when it arrests , allows us to estimate the loading force @xmath83 at the discrete arrest times @xmath84 ( that is , there is @xmath76 corresponding to each rupture event ) . for that aim , we neglect the contribution to the integral in the region @xmath85 , where @xmath86 is the location of the peak of @xmath14 slightly after a rupture front arrested ( cf . fig . [ fig : before_after ] ) , and then assume that @xmath87 can be replaced by a constant residual stress @xmath81 , obtaining @xmath88    to calculate @xmath86 , we note that at the arrest times @xmath76 eq . takes the form @xmath89 ( i.e. in the range @xmath90 and neglecting inertia ) . with the approximate boundary conditions @xmath91 , this equation can be readily solved as @xmath92 this can be substituted in eq . to give @xmath93 where @xmath94 was used ( which is , of course , valid at any time , not only at the discrete arrest times @xmath95 ) .    the prediction in eq . , i.e. @xmath96 , is tested in the inset of fig . [ fig : fd ] for all three models over many events ( i.e. this is a multiple - event property , not only a property of the first event , which was the focus of the discussion up to now ) . the analytic prediction is observed to be in favorable agreement with the numerical data for all three models , where the prefactor ( slope ) in the relation @xmath96 is the same for the ls and stl models , but is significantly smaller for the pw model . these results show that @xmath81 is the same for every rupture event and lend direct support to the assumption that spatial variations of the residual stress left behind _ any _ rupture front can be neglected , consistent with the explicit stress profiles shown in fig . [ fig : before_after ] ( for the first event in the three different models ) . the latter observation allows us to extract @xmath81 , the only unknown quantity in eq . ( all other quantities are known parameters , which are the same for all three models ) , yielding @xmath97 for the ls and stl models and @xmath98 for the pw model . the fact that the models that feature a nonmonotonic velocity dependence , i.e. the ls and stl models , give rise to an essentially identical residual stress @xmath81 is intimately related to the value of the steady state stress at the minimum of the friction curve ( cf . [ fig : ss ] ) , which is the same for both . equation then shows that the fact that the pw model produces a lower overall frictional resistance ( and deeper force drops ) compared to the ls and stl models is intimately related to the fact that the residual stress left behind the rupture fronts in the pw model is significantly lower than that of the ls and stl models . furthermore , eq . suggests an explanation for why the penetration depth , i.e. @xmath86 , is significantly larger in the pw model than in the other two models . the `` static friction coefficient '' @xmath99 is ordinarily defined as the tangential force , normalized by the normal force , needed to initiate global motion of the block . this force also corresponds to the peak of the loading curve . from this perspective , all of the spatiotemporal dynamics discussed up to now are precursory @xcite , as they precede global motion which sets in only when a rupture front reaches the leading edge of the block ( i.e. when @xmath100 ) . hence , we can estimate @xmath99 , which quantifies the global frictional resistance , as @xmath101 where eq . was used . this shows that the `` static '' frictional resistance of the interface , measured at slow loading velocities ( here @xmath48 m/s ) , is influenced by dynamic processes at much higher slip rates and furthermore that the existence of velocity - strengthening friction behavior strongly affects @xmath99 through @xmath81 @xcite . the results discussed above highlight two important points . first , an effectively constant residual stress @xmath81 is left behind rupture fronts in all of the models studied here . this property emerges spontaneously , unlike conventional slip - weakening models in which it is assumed a priori ( see , for example , @xcite and the discussion in @xcite ) . the value of @xmath81 depends on the existence of velocity - strengthening friction , which in turn has significant implications on the strength of the interface , as evident from fig . [ fig : fd ] and eq . . note also that the constancy of the residual stress @xmath81 implies that the mechanical fields associated with frictional shear cracks in 2d are well described by the classical theory of fracture @xcite . second , once @xmath81 is known , the arrest of rupture fronts is determined by global equilibrium conditions @xcite , rather than by dynamic considerations ( cf . as energy dissipation is at the heart of frictional phenomena , it will be interesting and instructive to consider the energy budget in the system . as a starting point , we briefly remind the reader that the linear momentum conservation law of eq . ( [ eq : momentum ] ) can be transformed into a continuity equation for the energy density ( using hooke s law and integration by parts ) . the result reads @xmath102 the first term is the rate of variation of the energy density ( both kinetic and elastic ) , and the second term is the divergence of the energy flux vector . their sum vanishes when energy is conserved . following the same procedure , one can derive the energy continuity equation for our model by combining eqs . and , obtaining @xmath103 where we defined @xmath104 here @xmath105 is the kinetic energy density , @xmath106 is the ( bulk ) linear elastic strain energy density , @xmath107 is the interfacial elastic energy density and @xmath108 is the energy flux . the interfacial energy density , @xmath107 , is dissipated during sliding due to the rupture of asperities , resulting in a dissipation rate @xmath109 @xcite , in addition to the standard dissipation rate @xmath110 . equation has the same structure as eq . , except for the non - vanishing dissipation rate @xmath111 , which exists because frictional dynamics are dissipative , and the existence of an interfacial elastic contribution @xmath107 ( both in the stored energy and in the dissipation power @xmath109 ) . velocity strengthening & absent & logarithmic & linear + rupture propagation speed & @xmath112 m / s & @xmath113 m / s & @xmath114 m / s + event s duration @xmath115 & @xmath116s & @xmath117s & @xmath118s + total dissipated energy . ] & 1.4 j / m & 0.26 j / m & 0.26 j / m + maximal dissipation rate . ] & 52 kw / m & 16 kw / m & 0.19 kw / m + total radiated energy . ] & 0.27 j / m & 4.5 mj / m & 0.54 @xmath80j / m + penetration length @xmath86 & 11.3 cm & 5.52 cm & 5.46 cm + * * *    the quantities defined in eq . are densities that exhibit complex spatiotemporal behaviors during frictional instabilities ( which result in rupture events ) . in order to gain some insight into these complex energy - exchange processes , it will be useful to consider the corresponding space - integrated quantities @xmath119 and @xmath120 . the interplay between these various quantities during frictional instabilities ( `` events '' ) , shown for all three models in fig . [ fig : partition ] , is an essential feature of interfacial dynamics . our goal is to quantify generic energy - exchange processes during frictional instabilities @xcite and in particular to understand the differences between the three models in this respect . as the dynamics during frictional instabilities are much faster than typical loading rates , we expect them to be exclusively driven by the already stored elastic energy . that is , we expect the rate of change of the sum of bulk and interfacial elastic energies , @xmath121 , to be negative during an event . figure [ fig : partition ] clearly demonstrates this , and that @xmath122 is negligible compared to @xmath123 ( hence we neglect the former compared to the latter in what follows ) . the time integral of @xmath123 over the event duration is the total energy released , which is a natural measure of the magnitude of the event ( other measures exist as well ) . the elastic energy released is either being dissipated directly or is being first transformed into kinetic energy ( `` radiation '' ) . eventually , the kinetic energy is also dissipated . this generic picture is demonstrated in fig . [ fig : partition ] for all three models . in particular , it is observed that the dissipation contributions @xmath124 and @xmath125 are comparable , where the former is typically larger than the latter . kinetic energy generation ( `` radiation '' ) , @xmath126 , is observed in the first part of the event . in the second part of the event @xmath127 , when the kinetic energy decays and is being dissipated .    while this generic qualitative picture is similar in all three models , there are large quantitative differences that we wish to discuss now . the main characteristics of the first rupture event in the ls , stl and pw models are summarized in table [ tab : dissipation ] . as we already know from fig . [ fig : fronts ] , the events are mediated by rupture fronts of vastly different velocities in the three models ( @xmath128 m/s in the pw model , @xmath129 m/s in the ls model and @xmath1301 m/s in the stl model ) . the event duration is about 40% larger in the pw model as compared to the ls model , both in the few @xmath131s range , while it is two orders of magnitude larger in the stl model ( @xmath130ms ) . despite the vast differences in the rupture propagation velocity and event duration , the total dissipated energy ( which equals the amount of elastic energy released during the event ) in the ls and stl models is essentially identical . this is in line with fig . [ fig : fd ] , which shows that the two models feature nearly identical stress drops and frictional resistance , and with fig . [ fig : before_after ] and the inset of fig . [ fig : fd ] , which show that the residual stress @xmath81 in the two models is essentially identical . this result clearly demonstrates that depending on the form of the velocity - strengthening friction branch ( e.g. logarithmic vs. linear ) one can observe events of the same magnitude ( i.e. integrated dissipation / energy release ) accompanied by very different dissipation rates ( see table [ tab : dissipation ] ) . this result might be related to geophysical observations indicating that slow rupture does not necessarily imply smaller integrated slip and energy release ( * ? ? ? * especially figure 5 ) . the total dissipation in the pw model is about @xmath132 times larger than the total dissipation in the ls and stl models , consistent with the much larger stress drops and the significantly reduced interfacial resistance observed in fig . [ fig : fd ] . moreover , the amount of kinetic energy generated during the event is much larger in the pw model as compared to the other two models , and is about 19% of the total energy released ( though eventually it is also dissipated ) . in systems of larger heights @xmath41 , this radiated kinetic energy will decay on longer timescales , allowing it to interact with remote boundaries . the kinetic energy generated in the stl model is negligibly small , while in the ls it makes about 1.7% of the released energy ( a similar value was reported in @xcite , although direct comparison is precarious ) . all in all , these results provide strong evidence that the existence and form of velocity - strengthening friction has significant implications on frictional dynamics and strength . in conclusion , by studying the spatiotemporal dynamics in three variants of a realistic rate - and - state friction model under quasi - static side - loading conditions , we showed that the existence and form of velocity - strengthening friction may significantly affect various aspects of the frictional response of interfaces . these include the propagation velocity of coherent fronts that mediate interfacial rupture events , the emergence of slow rupture , the elastic energy released during events ( i.e. their magnitude ) , the dissipation and radiation rates , and the global frictional resistance ( strength ) . the clear connection between the existence of velocity - strengthening friction and slow rupture appears to be directly related to the recent experimental results of @xcite . it is also shown that events of similar magnitude ( and hence stress drops ) can be accompanied by substantially different dissipation and kinetic energy radiation rates .    our theoretical results , together with extensive experimental evidence @xcite , highlight the need to quantitatively characterize the velocity - strengthening frictional response of interfaces , both experimentally and theoretically , and to systematically incorporate it into friction theory . since frictional instabilities spontaneously lead to accelerated slip that probes relatively high - velocity properties of frictional interfaces , the latter  which include velocity - strengthening friction  affect the frictional response even under quasi - static loading conditions . this understanding may offer new ways to interpret existing observations in a broad range of frictional systems and to develop predictive theories of the dynamics of spatially extended frictional interfaces .
Abstract:

Answer: frictional interfaces are abundant in natural and man - made systems and their dynamics still pose challenges of fundamental and technological importance . a recent extensive compilation of multiple - source experimental data has revealed that velocity - strengthening friction , where the steady - state frictional resistance increases with sliding velocity over some range , is a generic feature of such interfaces . moreover , velocity - strengthening friction has very recently been linked to slow laboratory earthquakes and stick - slip motion . here we elucidate the importance of velocity - strengthening friction by theoretically studying three variants of a realistic rate - and - state friction model . all variants feature identical logarithmic velocity - weakening friction at small sliding velocities , but differ in their higher velocity behaviors . by quantifying energy partition ( e.g. radiation and dissipation ) , the selection of interfacial rupture fronts and rupture arrest , we show that the presence or absence of velocity - strengthening friction can significantly affect the global interfacial resistance and the total energy released during frictional instabilities ( `` event magnitude '' ) . furthermore , we show that different forms of velocity - strengthening friction ( e.g. logarithmic vs. linear ) may result in events of similar magnitude , yet with dramatically different dissipation and radiation rates . this happens because the events are mediated by interfacial rupture fronts with vastly different propagation velocities , where stronger velocity - strengthening friction promotes slower rupture . these theoretical results may have significant implications on our understanding of frictional dynamics .


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: overwhelming observational evidence @xcite implies the existence of dark matter ( dm ) , whose nature is still unknown . searching for dm is now a major research theme in both particle physics and astrophysics . the former mainly focusses on direct detection by experiments sensitive to the interaction of dm particles with normal matter , while the latter uses indirect methods , including measurements of the end products of dm annihilation / decay ( e.g. @xmath1-ray , @xmath2 pairs ) or mapping of gravitationally lensed structures . a number of particle physics scenarios have been proposed for the nature of dm . among these scenarios , the most compelling is that of the weakly interacting massive particle ( wimp ) that acts as cold dark matter and effectively explains the origin of the large - scale structure of the universe . theories generally hypothesize wimps that interact with normal matter no more strongly than the weak force and therefore direct observation is challenging . if dm is a form of a thermal relic particle that was once in thermal equilibrium in the very early universe , then it may light up galactic substructures by self - pair annihilations . in addition , dm decay products may also be a source that is comparable to annihilation . the absolute @xmath1-ray intensity depends on the wimp annihilation cross section or decay rate , particle type , particle mass , and astrophysical distributions , all of which are poorly known . the relative annihilation / decay luminosity is density - dependent , which may offer insight into the nature of dm , e.g. the decay signal is proportional to density while the annihilation luminosity is proportional to density squared . additional correction factors may be applied to the standard velocity - weighted thermal cross section @xmath3 . one example is the sommerfeld boost  a nonperturbative increase in the annihilation cross - section at low velocities that is the result of a generic attractive force between the incident dm particles . baryons could also play an important role in determining the profiles of dm halos . adiabatic contraction induced by the infall of baryonic matter could steepen the dm density profile substantially , e.g. to @xmath4 from a navarro - frenk - white ( nfw ) @xcite profile @xmath5 @xcite . another possible correction to the annihilation cross - section that has not been previously discussed in the galactic context is that of @xmath6-wave suppressed annihilations , as might be important for spin - dependent interactions . diffuse galactic @xmath1-ray emission ( dge ) is believed to be produced by interactions between cosmic rays and the interstellar medium ( ism ) . the launch of the fermi gamma - ray satellite , with its large area telescope ( lat ) , enabled a detailed study of cosmic - ray origin and propagation , and of the interstellar medium . in a recent study @xcite , the fermi - lat team published an analysis of the measurements of the diffuse @xmath1-ray emission from the first 21 months of the fermi mission . they compared the data with models generated by the galprop code , and showed that these emission templates underpredict the fermi - lat data at energies of a few gev in the inner galaxy . this can possibly be explained by undetected point - source populations and variations of the cosmic ray spectra . @xcite have performed an analysis of the dge in the milky way halo region to search for a signal from dark matter annihilations or decays . they considered a large region covering the central part of the galactic halo while masking out the galactic plane . in such a region , the dm signal would have a large s / n ratio and would not depend on detailed assumptions about the center profile , e.g. the assumptions of a nfw profile or a cored profile would only differ by a factor of @xmath7 . this paper provided conservative limits on the dm cross - section assuming that all the @xmath1-ray signal comes from dm in this region . they also provided more stringent limits based on modeling the foreground @xmath1-ray signal with the galprop code . their results impact the possible mass range over which dm is produced thermally in the early universe and challenge the dm annihilation interpretation of pamela / fermi - lat cosmic ray anomalies . however , several points could be improved : 1 ) although using a nfw profile for modeling pure dm and a cored profile for modeling the baryonic matter effect covers extreme cases , it is still worth trying a more realistic profile directly from simulations ; 2 ) in addition to the mass and cross - section of the dm particle , the different annihilation / decay schemes may also play an important role , e.g. subhalo effects would become considerably strong in the region of interest where tidal disruption plays a role and especially when considering possible sommerfeld enhancements ; and 3 ) instead of detecting dm signal from the region that masks out the galactic center region , it may be of interest to search for a possible signal from the center or anticenter regions . therefore , a good template for the galactic center or the full sky should provide more constraints on the dm physics ; 4 ) the spatially dependent velocity dispersion of the dm due to substructure plays an important role both in considering sommerfeld or @xmath8wave enhanced annihilation signals , and simulations can provide this information . the contribution of dm annihilation / decay processes to the fermi - lat @xmath1-ray signal in the galactic center ( or its expected signal from subhalos in the galactic halo ) is difficult to quantify without knowing the galactic distribution of dark matter and the dark matter particle properties leading to annihilations or decays . to compare the observations with predictions of dm emission , we must rely on numerical simulations that follow the formation of cosmic structure in the highly non - linear regime . a few high - resolution simulations of milky way - sized halos have been completed over the last few years : the aquarius project @xcite , the via lactea series @xcite , and the ghalo run @xcite . the ability of the fermi - lat satellite to detect dm annihilation signals from galactic subhalos has been previously studied by @xcite and @xcite using the high - resolution via lactea ii ( vlii ) n - body simulation . to date , there has been no clear signal of dark matter annihilations or decays in the outer galactic halo . there are numerous possible explanations for non - detection in the outer galactic halo combined with the observed dge in the galactic center , including simulated dm - only structure being inconsistent with the actual galactic structure , the predicted dm annihilation / decay luminosity from particle physics being overly optimistic , or some misunderstanding of the observational data . motivated by the range of possible explanations for the excess of dge data to the model , we present a study of various possible galactic dm contributions to the observed dge emission . we construct several maps of gamma - ray emission from different dm annihilation / decay scenarios in order to illustrate how each scenario may offer additional templates to fit the fermi data . we additionally compared dm annihilation signals in each scenario to the dge model from the galprop code of @xcite in different regions on the sky . we organize the rest of this paper as follows : in @xmath9[sec : darkmatter ] , we give a description of different dm annihilation / decay schemes . our data processing method and @xmath1-ray generating algorithm are described in @xmath9[sec : simulation ] . we show how we add our correction to the normal annihilations / decays in @xmath9[sec : corrections ] . the all - sky maps results and comparisons to fermi data for different sky regions are presented in  [ sec : regions ] . conclusions and discussion are presented in  [ sec : discussion ] . gamma - rays are one of the final products from dm self - annihilations or decays . in the case of annihilations , the @xmath1-ray flux ( photons @xmath10 sr @xmath11s@xmath11 ) in a solid angle d@xmath12 of a given line - of - sight @xmath13 can be written as @xmath14 where @xmath3 is the velocity - weighted thermal cross section , which is usually treated as a constant , and @xmath15 is the density of dm at position @xmath16 , and @xmath17 is the distance from the dm particle to the observer . in the case of decay , the @xmath1-ray intensity is then proportional to @xmath18 , @xmath19 where @xmath20 is the decay lifetime . eqn-([eqn : dm_ann ] ) and eqn-([eqn : dm_decay ] ) are usually written as a combination of a particle physics factor @xmath21 and an astrophysical factor @xmath22 , while we combined them as @xmath23 to indicate what has been considered in our @xmath1-ray producing code . since @xmath24 and @xmath20 are unknown , the normalization is arbitrary . for a given astrophysical dm distribution , the @xmath23 factors given here determine the relative intensity of the predicted flux . calculating the absolute value of @xmath23 requires the absolute value of the particle physics parameters . as shown in several papers @xcite , @xmath24 is less than @xmath25 @xmath26s@xmath11 , and particle mass @xmath27 varies from 100 gev to 10 tev for the case of wimps . for our purpose , only the relative intensity is needed . the @xmath1-ray maps in this work are produced from the simulation data ( see  [ sec : simulation ] ) by using one of these formulae in addition to possible cross - section correction factors . we consider in total six different @xmath1-ray map production scenarios : ( a ) pure annihilation without any correction factors , assuming @xmath28 is a constant ; ( b ) annihilation with a sommerfeld enhancement of @xmath29 included in the annihilation cross section ( @xmath30 is the relative velocity of the annihilating particles ) ; ( c ) annihilation with a sommerfeld enhancement of @xmath31 , as expected near resonance ; ( d ) annihilation proportional to @xmath32 consistent with @xmath33wave channel suppression ; ( e ) adiabatic contraction to the dark matter density profiles , and ( f ) pure decay , taken to be linear in dark matter density . the correction cases are described in @xmath9[sec : corrections ] . the via lactea series includes some of the highest resolution collisionless simulations of the assembly of a milky way - sized halo . vli , the first of the series , contains @xmath34 dm particles , covering the virial volume and surroundings of a host halo of @xmath35 . the host halo and subhalo properties are presented in @xcite , @xcite , and @xcite . for this work , we used the second generation of the simulation series , the vlii run @xcite , which has a slightly higher resolution . it employs about 1 billion particles each of mass @xmath36 to simulate a milky - way - sized host halo and its substructures . we extract the roughly @xmath37 particles of the @xmath38 snapshot within @xmath39 ( @xmath40 kpc , the radius where the density is 200 times larger than the critical density @xmath41 ) of the host halo . about @xmath42 individual subhalos have been identified in this region . the halos found within @xmath39 show considerable self - similarity . an nfw - like profile fitting to the density could be generalized by , @xmath43^{\frac{\beta-\alpha}{\alpha}}},\ ] ] @xcite , fixes @xmath44 , @xmath45 and gives the best - fitting parameters @xmath46 , @xmath47 kpc and @xmath48pc@xmath49 . the host halo is not spherically symmetric but is ellipsoidally shaped . the detailed host halo properties can be found in @xcite .    to produce the @xmath1-ray maps , we follow @xcite and locate a fiducial ` observer ' 8 kpc from the host halo center along the intermediate principle axis of the ellipsoid . given the position of an observer , the density field of each particle can be represented by a dirac-@xmath50 function . eqn-([eqn : dm_ann ] ) and ( [ eqn : dm_decay ] ) can be rewritten as the summation of flux over each particle , i.e. @xmath51 . for annihilation , @xmath52 ; and decay , @xmath53 , where @xmath54 is the distance from the particle to the observer . to avoid shot noise , we smoothed the particles with a sph kernel so that each particle is represented by a sphere of radius @xmath55 instead of a point . the flux @xmath56 , correspondingly , spreads out as a gaussian on the sky . we computed @xmath57 , @xmath55 and @xmath58 ( the velocity dispersion ) by finding the volume encompassing the nearest 32 particles ( using the code smooth ) . the velocity dispersion @xmath58 is read out simultaneously from this process . our final result for the omnidirectional @xmath1-ray signal has been calculated using healpix @xcite . we use a nside=512 to model the angular resolution of fermi - lat ( roughly 0.2@xmath59 ) . this is an extremely compute - intensive task , it takes more than 8 hours on a regular cpu to produce a single map with the 400 m particles . in order to be able to create the maps involving the different annihilation scenarios and different viewpoints , we have developed a novel gpu - based algorithm@xcite .    in this method , we projected @xmath60 onto two separate tangential planes of the two celestial hemispheres using stereographic projection , then remapped these to the final healpix projections . we have realized that the problem can be reduced to rendering a projected density profile of a particle , weighted by different factors depending on the physics of the annihilation . as these computations became quite similar to those used in visualization , we were able to implement them in the shader language of opengl ( open graphics library ) on a high - end nvidia gpu for maximum performance . as a result , we are able to render the projected profiles at a rate of more than 10 m particles per second and hence build a map in 24 seconds . this thousand - fold speedup was an essential factor , that enabled us to carry out the large number of numerical experiments needed to complete this paper . even with the high resolution of vlii , only a portion of the hierarchal structure of dm clumps is resolvable . in principle , the wimp dm substructure has mass power all the way down to the mass of the dm kinetic or thermal decoupling scales , of order @xmath61 these unresolved structures will boost the final brightness of the @xmath1-ray emission according to the @xmath62 dependence ( while decay is not affected due to the @xmath18 dependence ) . detailed studies of the boost factor @xmath63 can be found in @xcite . although as pointed out by springel et al @xcite and fornasa et al . @xcite the unresolved subhalo component provides a important source of emission , springel et al also note that this unresolved smooth component is quite flat ( i.e. the contrast of the brightest and faintest points is only about 1.5 ) . in practice , this component is unlikely to be distinguishable from other emission components . to include the emission of unresolved subhalos , we suggest that use of an all - sky uniform component would be an alternative solution . on the other hand , the boost factor is also highly model - dependent ( e.g. kamionkoski et al . @xcite ) . in terms of simplicity and clarity , we choose to not include the boost factor calculation in this work , since we only care about the possible shapes of the profiles .     and the dotted dash line is the nfw@xmath65 profile describing the adiabatic contraction with @xmath66 and @xmath67 . the normalization is chosen such that the intersection point of the contracted and non - contracted profiles are at a density of 1 . see text for details.,scaledwidth=49.0% ] the pamela @xcite and ams-02 results @xcite show a rise in positron fraction at high energy and a hardening of the spectral index , while no antiproton excesses were found . if this signal were due to annihilation , the standard thermal cross - section @xmath24 would be too small to simultaneously fit @xmath2 and avoid antiproton excesses . to gain a larger cross - section , a force carrier @xmath68 @xcite is proposed for mediating the dm particle interaction . this mediator could be a standard model particle or an unknown boson responsible for dark sector forces . assuming the dm particle is a majorana or dirac fermion , such a process is denoted as @xmath69 . the sommerfeld enhancement @xcite can be therefore calculated . for the resonance case ( when @xmath70 , where @xmath71 is the coupling strength , and @xmath72 is an integer ) , this model gives an enhancement of @xmath73 to the cross - section . for the non - resonance case , @xmath74 . for @xmath75 , @xmath76 . the sommerfeld cross - section must saturate in a viable model : we take this saturation to be @xmath77 km / s , above the resolution limit of our simulations . from our point of view , this amounts to a renormalisation that we ignore : it is the profile shape and clumpiness that we care about . the cross - section for dm particle pair annihilation changes significantly over the simulation range caused by velocity changes in different environments , i.e. by a factor @xmath78 or more @xcite . in such a case , the subhalos , whose velocity dispersion is less than that of the central halo , light up . for the @xmath31 case , the substructures almost dominate , as shown in fig-[fig : allskymap](c ) .    to calculate an accurate sommerfeld enhancement value at a given location , one needs the whole phase - space distribution of each particle . this is not possible for our case using the vlii simulation . we therefore follow @xcite and assume a boltzmann velocity distribution of the relative velocity @xmath79 : @xmath80.\ ] ] the sommerfeld enhancement is then given by @xmath81 . if @xmath82 has no saturation , then @xmath83 is then @xmath84 for the @xmath29 case and @xmath85 for the @xmath31 case . if one takes into consideration the saturation velocity @xmath86 , below which @xmath82 no longer increases , we find @xmath83 presents similar properties . we therefore approximate @xmath83 completely by @xmath82 without losing much accuracy . in our calculation , we calculated @xmath87 as @xmath88 or @xmath85 when @xmath89 , and as @xmath90 when @xmath91 . here @xmath92 is the velocity dispersion mentioned in the last section . fig-[fig : allskymap ] ( b ) and ( c ) shows the sommerfeld enhancement case for the @xmath1-ray map . to avoid the cross - section blowing up at very low velocity , we applied saturation velocities equal to 1 km / s for the @xmath29 case and 5 km / s for the @xmath31 case . as indicated in the sommerfeld enhancement case , the relative velocity of particles varies significantly after thermal freeze out . if the cross - section is velocity - dependent , it will either suppress high velocity or low velocity annihilation . constraints on annihilation processes from cmb and @xmath1-ray observations could be considerably weakened if the @xmath33wave channel were suppressed . @xmath8wave annihilations are orders of magnitudes larger than @xmath33wave annihilations at recombination and are an interesting case to be considered . for example , neutralino annihilation could be dominated by the @xmath8wave process @xcite . in contrast to the sommerfeld case , we consider an annihilation case whose @xmath33wave process is suppressed . for the annihilation rate , we phenomenologically set @xmath93 , where @xmath94 for @xmath33wave(@xmath8wave ) annihilation @xcite . fig-[fig : allskymap ] ( d ) shows a map from the @xmath8wave annihilation with the @xmath33wave channel suppressed . more discussion will be given in  [ sec : discussion ] . when structure is forming , baryonic matter dissipates its thermal and kinetic energy and falls to the center of the dm halo . the condensation of gas and stars in the inner regions of dm halos will adiabatically contract the dm density distribution and lead to a denser profile in the center . this effect is implied from both theoretical and observational studies @xcite . adiabatic contraction is generally studied using the standard contraction model ( sac ) as introduced by @xcite . however , hydrodynamic simulations show that the contraction is weaker and a modified adiabatic contraction model was introduced in @xcite . we consider the baryonic matter contraction effect only in the host halo , where there is enough baryonic matter to cause the contraction . to simplify the calculation , we follow @xcite and modify the host halo density profile to be a nfw@xmath95 profile with @xmath96 , @xmath97 , @xmath98 and @xmath99kpc . to conserve the the total mass of the host halo within @xmath39 , we use @xmath100 . we convolve a contraction factor @xmath101 with the density of each particle in the simulation . the mass of each particle is changed accordingly . the contracted density profile and the original profile are shown in fig-[fig : densityprofile ] . we normalized the profile such that the intersection point of the original and the contracted profile is at a density of @xmath102 . we include the nfw and nfw@xmath95 profile in fig-[fig : densityprofile ] as points of comparison . fig-[fig : allskymap ] ( e ) shows an annihilation map of the contracted host halo plus uncontracted subhalos . the all - sky maps of the six different annihilation / decay scenarios are shown in fig-[fig : allskymap ] . the normalization is chosen such that the host halo has flux unity . all the maps are smoothed by a gaussian kernel with fwhm 0.5@xmath59 . they show significant differences in both the center halo profile and the relative subhalo flux . compared to the pure annihilation case of fig-[fig : allskymap]a , the @xmath29 and @xmath31 sommerfeld enhancement cases ( fig-[fig : allskymap]b ) have considerably more substructures visible . the @xmath31 resonant sommerfeld enhancement case has the largest contribution from substructures to the annihilation flux of the six scenarios . for the @xmath8wave annihilation case ( fig-[fig : allskymap]d ) , the adiabatic contraction case ( fig-[fig : allskymap]e ) and pure decay case ( fig-[fig : allskymap]f ) , substructures are less visible than in the pure annihilation scenario . the smooth component also differs for the six scenarios in terms of both the central maximum and the profile slope , as shown in fig-[fig : fluxprof ] . the adiabatic contraction case has a steep central profile such that the normalization of the total flux overwhelms the relatively small contribution from substructure . it looks more spherical in the all - sky map because of the assumption of spherical nfw@xmath95 contraction ( fig-[fig : allskymap]e ) . both sommefeld enhancement cases have a larger central flux than the pure annihilation case due to the distribution of the velocity dispersion . with a @xmath32 dependence , the @xmath8wave annihilation case therefore has a flatter central flux profile . the pure decay case has the shallowest all - sky map and radial profile because of the linear @xmath18 dependence . we compare our results to the fermi - lat dge model . @xcite presented a measurement of the first 21 months of the fermi - lat mission and compared a grid of models of dge emission produced by the galprop code . these models incorporate the observed astrophysical distribution of cosmic - ray sources , interstellar gas and radiation fields . they compare the predicted model intensities and spectra with the observations in various sky regions . the models are consistent with the fermi data in the anti - center regions but underpredict the data in the inner regions of the galaxy at energies around a few gev . their conclusions concerning the discrepancy mainly focused on undetected point sources and cosmic ray spectral index variations throughout the galaxy . although it is generally accepted that there is an  excess " @xmath1-ray signal most probably from the galactic center , it is still worth considering signals from off - center regions .    in @xcite , various sky regions are taken into consideration for the comparison of models and observations . we take the same sky regions to plot the angular profiles , namely , @xmath103 , @xmath104 , @xmath105 and @xmath106 . we also show the dge model in each subplot . this dge model is based on the galprop model @xmath107 . here we do not make comparisons to models with other choices of parameters . however , since all the models are required to fit the observational data , they show very similar features . in order to get the total gamma ray flux we integrated over energy , from 50mev to 800gev , in the galprop model . the results are shown in fig-[fig : cutplots ] . the upper left panel shows the longitudinal profile of the region of the sky with the central @xmath108 removed , while the upper right panel shows the longitudinal profile of the central @xmath109 of the galaxy . the lower left panel shows the latitude profiles of the region with @xmath105 . the lower right panel shows the latitude profiles of the region with @xmath110 removed . two apparent large subhalos that are visible in the all - sky maps ( fig-[fig : allskymap ] ) would make the latitude profiles off - central in this panel . to better illustrate the comparison between the fermi - lat model and the profile from our model , we remove the influence from the 2 sub halos by separating the smooth main halo component and the subhalo component . we select a random position sitting 8kpc away from the gc to get an all - sky map of the subhalo component such that there is no subhalo too bright or too close to the observer . after that , we combine the subhalo component and main halo component to get a properly centered profile . the variations in the longitudinal and latitude profiles of the annihilation cases are the result of the contribution from the subhalos to the annihilation signal . the contribution from the subhalos is most important for the two sommerfeld - corrected annihilation cases . the longitudinal profile of @xmath31 sommerfeld - corrected case in the upper left panel is completely dominated by emission from subhalos . note that these profiles , if used to interpret the fermi data , must be combined with the main dge model , since they are in any case a small contribution to the full @xmath1-ray emission . these comparisons provide us with intuition on how these different profiles would work if used to fit the real data . for example , in the upper panels , the dge signal shows a relatively flat profile compared to the longitudinal profiles of the five annihilation cases . the decay profile is the most degenerate with the dge signal and would provide a relatively larger weight when fit to data . on the contrary , in the lower two panels , due to the disk component of the milky way , the dge signal has a steep central profile that most closely resembles the latitude profile of the pure annihilation case . it may lead to a larger weight to the pure annihilation contribution . the subhalos contribution can not be directly used as a part of the profile model since their positions are arbitrary compared to the real case . as shown in the figures , this contribution is dominated by fluctuations and has a flat envelope . they provide an additional flat smooth component if unresolved ; and a measure of the statistical uncertainties , if resolved .     correction . _ thick dash - dotted _ : annihilation with sommerfeld enhancement @xmath31 correction . _ dotted _ : annihilation with @xmath32 correction . _ thick dashed _ : annihilation with adiabatic contraction . _ thin dashed _ : pure decay . the fluxes are normalized such that the host halo has flux unity . see text for details . , scaledwidth=49.0% ]    we take the inner galactic region into special consideration since it might be a promising region for detecting dm signals . other works propose alternative astrophysical models for the @xmath1-ray excess at the galactic center , such as unresolved millisecond pulsars @xcite . @xcite strongly disfavor this explanation for both the spatial distribution and spectrum of observed millisecond pulsars . assuming that the @xmath1-ray excess mainly comes from dm decays / annihilations , we show here the predicted flux map ( fig-[fig : center55 ] ) and angular profile ( fig-[fig : fluxprof ] right panel ) of the inner @xmath111 region of the simulated maps for the different @xmath1-ray generating cases . the angular profiles in fig-[fig : fluxprof ] right panel are normalized to the total flux of the host halo , as used in the all - sky maps . the flux maps in fig-[fig : center55 ] are normalized such that the center pixel has flux unity and the color scale is linear . the sommerfeld - enhanced ( fig-[fig : center55]b and fig-[fig : center55]c ) cases and adiabatic contraction case ( fig-[fig : center55]e ) have a more contracted profile than the pure annihilation case ( fig-[fig : center55]a ) . in contrast , the @xmath8wave annihilation ( fig-[fig : center55]d ) and pure decay case ( fig-[fig : center55]c ) have flatter profiles illustrated by the greater spread in @xmath1-ray flux over the @xmath111 region . the six cases of the @xmath1-ray flux maps of fig-[fig : center55 ] show considerable variations and complex asymmetrical contours in the central region . this asymmetry would be absent from simple parameterized model of the the dm annihilation signal from a nfw / cored profile considered previously @xcite . note that the simulation used in this paper has a convergence radius @xmath112 pc which corresponds to approximately @xmath113 in the center region . therefore , our method would underestimate the real flux at the innermost region . to estimate this underestimation , we simply assume the real converged density profile is equation-[eqn : nfw ] and the simulation profile is a cored profile @xmath114\}$ ] with @xmath115kpc , @xmath116pc@xmath49@xcite . it turns out for annihilation , the nfw flux is 2 times larger of the cored profile occured at @xmath117 ; for decay , it occurred at @xmath118 . when using the templates to fit the real data , one should pay special attention the central few pixels .    . the legend is the same as fig . [ fig : fluxfunc ] . see text for details . , scaledwidth=49.0% ] the morphologies and positions of each subhalo differ completely in different realizations of simulations , and may not directly apply to fitting templates to the observational data , whereas their statistical properties can be well studied using the maps produced by our tools . one useful statistical property is the angular power spectrum ( aps ) of anisotropies . as pointed out by @xcite , aps could be used along with the energy spectrum to decompose the emission components . @xcite did actual comparisons using simulation - generated aps @xcite with fermi - lat measurements @xcite , yielding reasonable preliminary results that constrain the dm annihilation cross section . our results which consider six different physical scenarios present an extension of these results . we leave a more detailed study of anisotropies ( which must have the boost factor carefully calculated ) to a future project . in addition to the diffuse emission , we are able to use our maps to study the point sources from fermi - lat . for annihilations , when we add in the sommerfeld enhancement , the number density of detected small halos increases . for the @xmath31 case , the substructure almost dominates the radiation field other than in the galactic center . fig-[fig : fluxfunc ] shows the flux functions of sub - halos for different cases . the flux of the subhalos are calculated by summing over the flux of all pixels inside the halo s angular radius . all the curves converge at @xmath119 since the total number of sub - halos in the simulation is fixed at @xmath120 . the curve for the sommerfeld enhancement @xmath31 is much steeper than the other curves , which indicates that the subhalos are brighter than in other models . we further plot the detectable point sources ( with angle extending less than @xmath121 @xcite ) in fig-[fig : numfunc ] . compared to @xcite , these curves are shallower than the identified and unidentified sources in the second fermi - lat source catalog for the inner region , but are relatively flat at larger latitudes . all of these features imply that observable dm signals could come from off - center subhalos .    in summary , we are not aiming to set new constraints on dm annihilations or decays , but rather to provide possible new templates ( the smooth components ) to fit the dge data with dark matter models and illustrate the possible range of profiles for plausible models . compared to a simple assumption that the dm halo has a nfw or an isothermal profile , our results based on the simulation data are more realistic ( e.g. fig-[fig : densityprofile ] ) . the effects of subhalos could be also be statistically included in the templates , e.g. aps or used to provide as estimates of the uncertainties . the galactic central region might be a promising region for detecting a dm signal as recently pointed out by @xcite . the various dm maps we present here show significant observational differences when the dm annihilation / decay physics are changed . the morphology of the center region is also asymmetric and complex , and simple analytical templates may fail to match the real signal . there are three key points that we wish to emphasize here . firstly , compared to the fermi dge result , the annihilation signal , with or without sommerfeld enhancement , has distinctly different central profiles . the shape of the fermi excess profile in the fermi bubble region has recently been used as a means of discriminating dark matter from millisecond pulsar models of the central excess @xcite . our philosophy here is similar but we extend this approach to probing the nature of dm annihilations . we predict different profiles for the alternative dark matter annihilation models . the differences are especially notable for decaying dark matter , for which the flux profile is very flat . secondly , it is clear that some of our models , especially the sommerfeld cases , provide more bright substructures than are allowed by the observations in terms of unresolved and unidentified point sources . we leave it to others to make detailed comparisons with the data , but our predicted cumulative counts of subhalos should provide a basic comparison . the @xmath8wave annihilation case is a promising model since it removes many of these bright substructures . the model including adiabatic contraction has a similar effect of de - emphasizing the substructure by increasing the contrast between host halo core and the ( unadiabatically contracted ) subhalos , although the quantitative role of adiabatic contraction of the dm in halos remains controversial . finally we note that we have not included any spectral signatures , which differ among the various adopted annihilation channels and particle masses : for example the sommerfeld models require high values for the particle mass ( @xmath122 tev ) and hence would be best constrained through observational data at relatively high photon energies . we thank michael kuhlen for his help on developing our algorithm . we acknowledge support from national science foundation grant oia-1124403 , and oci-1040114 . the research of js has also been supported at iap by erc project 267117 ( dark ) hosted by universit pierre et marie curie - paris 6 . support for p.m. was provided by the nsf through grants oia-1124453 , and by nasa through grant nnx12af87 g . the research of a.s . and of rfgw was also supported by grant 109285 from the gordon and betty moore foundation .
Abstract:

Answer: observations of diffuse galactic gamma ray emission ( dge ) by the _ fermi _ large area telescope ( lat ) allow a detailed study of cosmic rays and the interstellar medium . however , diffuse emission models of the inner galaxy underpredict the fermi - lat data at energies above a few gev and hint at possible non - astrophysical sources including dark matter ( dm ) annihilations or decays . we present a study of the possible emission components from dm using the high - resolution via lactea ii n - body simulation of a milky way - sized dm halo . we generate full - sky maps of dm annihilation and decay signals that include modeling of the adiabatic contraction of the host density profile , sommerfeld - enhanced dm annihilations , @xmath0-wave annihilations , and decaying dm . we compare our results with the dge models produced by the fermi - lat team over different sky regions , including the galactic center , high galactic latitudes , and the galactic anti - center . this work provides possible smooth component templates of dm to fit the observational data . the subhalo contributions can be considered to provide statistically meaningful templates , and demonstrate how spatial profiles are significantly modified according to different annihilation / decay scenarios . we argue that a subhalo - based approach can help constrain the dm physics .