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: a number of astronomical observations , such as type ia supernovae ( snia ) @xcite , the cosmic microwave background ( cmb ) anisotropy @xcite , and the large scale structure ( lss ) @xcite , have shown that the expansion of our universe is accelerating . the acceleration of the universe strongly indicates the existence of a mysterious exotic matter , named dark energy , which has negative pressure and takes the largest proportion of the total density in the universe . the ultimate fate of our universe is to be determined by the dark energy , yet we know little about its properties . furthermore , while both the dark matter and the dark energy is commonly thought to interact weakly with ordinary matter , it is not immediately clear whether there is some interaction between the two , or if the two are in some way linked .    in this paper we consider the dark matter - dark energy interaction in the framework of the holographic dark energy ( hde ) model @xcite . the basic ideas of the hde model are based on the holographic principle @xcite , which was inspired by the bekenstein entropy bound of black holes@xcite . in this model , the dark energy is related the quantum zero - point energy density @xmath0 , the total energy of the whole system with size @xmath1 should not exceed the mass of a black hole of the same size , thus we have @xmath2 . the largest infrared cut - off @xmath1 is chosen by saturating the inequality so that the hde density is @xmath3where @xmath4 is a numerical constant , and @xmath5 is the reduced planck mass . if we take @xmath1 as the size of the current universe , for instance the hubble scale @xmath6 , then the dark energy density estimated with eq.([holo - density ] ) will be close to the observed value . however , hsu @xcite pointed out that this yields a wrong equation of state for dark energy . li @xcite subsequently proposed that the future event horizon might be chosen as the ir cut - off @xmath1 , and viable hde models can then be constructed . as demonstrated in ref . @xcite , the hde scenario may help to answer the _ fine tuning problem _ and the _ coincidence problem _ of dark energy . various aspects of the hde model and its extensions have been investigated subsequently . for example , it has been generalized to the case of non - flat universe @xcite , and it was studied from the perspective of the anthropic principle @xcite . the holographic dark energy was also considered for driving the inflation , and its resulting power spectrum index calculated @xcite . ref . @xcite defined a new future event horizon as the infrared cut - off @xmath1 and found that the coincidence problem could be alleviated to some extent . @xcite investigated the possible form of the cut - off scale @xmath1 for the solution of coincidence problem . inspired by the hde model , ref . @xcite proposed the ricci scalar curvature radius as @xmath1 . the model was reconstructed with a scalar field @xcite . observational constraints have also been studied @xcite . the growth of density perturbations in the holographic dark energy models has been studied in @xcite . interactions between the de and dm components were considered for the hde model @xcite . ref.@xcite suggested that for the @xmath7 case , with interactions between the two dark components , the dark energy could acquire an equation of state in agreement with observations , thus partly solving the no - acceleration problem noted in ref.@xcite for interaction - free models . ref.@xcite generalized the discussion in ref.@xcite to the case of a future event horizon . in that paper the proposed form of interaction is @xmath8 , which could also lead to a negative dark energy equation of state . however , ref.@xcite pointed out that the interaction terms such as @xmath9 , @xmath10 , @xmath11 all lead to unphysical results , because the dark matter density could be driven to a negative value in these models .    here we propose a form of interaction which we considered _ natural _ and _ physically plausible_. in general , when two types of matter interact with each other , the interaction should be proportional to the product of the powers of the number density of both . this is most easily understood in the following way : if one component does not exist , then there certainly would be no interaction at all ; with its amount increasing , the interaction rate should increase . thus , at least in the diffuse limit , the interaction rate should be proportional to some powers of the density ( the relation could become more complicated when the densities become high ) . for instance , this is the case for ordinary two - body chemical reactions , namely @xmath12i.e . the reaction rate is proportional to the product of number densities of the two particle species . however , besides the interaction in the usual sense ( `` scattering '' ) , one may also consider a more general case where the dark energy is converted to dark matter or vice versa by `` annihilation '' or `` decay '' . in such a case , the `` interaction term '' or more precisely the energy exchange term would depend only on one of the densities and not on the other .    in general , we consider the interaction term to be proportional to some powers of the density , @xmath13@xmath14 in which @xmath15 and @xmath0 are the energy densities of the matter and hde respectively , and @xmath16 is a parameter with mass dimension @xmath17}^{-1}$ ] . the power index @xmath18 is assumed to satisfy @xmath19 where if one of them is 0 , we will have the case of `` decay '' ( if the other is 1 ) or `` annihilation '' ( if the other is 2 ) . if @xmath20 , the interaction suggests that dark energy is converted into dark matter , while @xmath21 suggests the inverse process . this form of interaction was considered in the study of coupled quintessence model @xcite . however , not much research on this form of interaction was carried out to understand how such a coupling term would affect the evolution of the universe .    in the next section , we present the evolution equations for the interacting holographic dark energy model ( details of our derivation are given in the appendix ) . in sec . 3 , we study the dynamical behavior of the model . in sec.4 we explore the parameter space and find the best - fit and allowed region of @xmath18 . finally , concluding remarks are given in sec . now let us consider a spatially flat friedmann - robertson - walker ( frw ) universe filled with matter component @xmath15 , hde component @xmath0 and the radiation @xmath22 . for simplicity we do not distinguish dark matter and ordinary matter below but incorporate both in the matter component @xmath23 . it is certainly possible or even likely that dark energy only couples to dark matter . however , in that case the baryons simply evolve as @xmath24 where @xmath25 is the scale factor , the baryon component would then become dynamically irrelevant as the universe expands , so it only affects the details of the model but not its qualitative behavior , hence we will neglect this effect in the present paper . the friedmann equation reads @xmath26where @xmath22 is the energy density of the radiation which satisfies @xmath27 . thus , eqs . ( [ holo - density ] ) , ( [ dark matter ] ) , ( [ dark energy ] ) , ( [ friedman1 ] ) are the dynamic equations which uniquely determine the evolution of @xmath28 @xmath29 @xmath30 @xmath31 and @xmath32 . in eq.([holo - density ] ) , @xmath1 could be hubble horizon @xmath6 , particle horizon and future event horizon , i.e. @xmath33or @xmath34we define the following dimensionless quantities : @xmath35,@xmath36 , @xmath37 , where @xmath38 is the current critical density , and @xmath39 , and introduce @xmath40 the dimensionless @xmath41 parameter determines the characteristic strength of the interaction .    the equations of motion of the hde model for the particle ( event ) horizon cases are @xmath42,\notag \\ \frac{de(z)}{dz } & = & \frac{1}{2(1+z)}\left[3e(z)+\frac{(\tilde{\rho}_{r}-\tilde{\rho}_{\lambda } ) } { e(z ) } \pm \frac{2 \tilde{\rho}_{\lambda } ^{\frac{3}{2}}}{c\cdot e^{2}(z ) } -\frac{\lambda ( e(z)^{2}-\tilde{\rho}_{r}-\tilde{\rho}_{\lambda } ) ^{\alpha } \tilde{\rho}_{\lambda } ^{\beta } } { e(z)^{2}}\right ] ,   \notag \\ w_{\lambda } ( z ) & = & -\frac{1}{3}\pm \frac{2}{3}\frac{\tilde{\rho}_{\lambda } ^{\frac{1}{2}}}{c\cdot e(x ) } -\frac{\lambda } { 3}\frac{1}{e(x)}\left[e(x)^{2}-\tilde{\rho}_{\lambda } -\tilde{\rho}_{r}\right]^{\alpha } \tilde{\rho}_{\lambda } ^{\beta -1 } ,   \label{event case}\end{aligned}\]]in which the upper(lower ) sign represents the case of particle(future event ) horizon ( c.f . the appendix for details of derivation ) . it is easy to verify that if @xmath1 in eq.([holo - density ] ) equals the particle horizon @xmath43 or hubble horizon @xmath6 , the model can not fit the current observational data , because current data prefer the dark energy with equation of state near @xmath44 , while in the above two cases the equation of state @xmath45 will be always greater than @xmath46 ( see also the analysis in ref.@xcite ) . therefore , in the following we will focus on the interacting model with @xmath47 ( see eq.([event case ] ) ) . in this section we study the dynamical behavior in models with interaction between de and dm , which may give rise to some interesting behaviors . first we consider the case where @xmath18 are fixed to be 1 , then consider the more general cases .    as a typical example we consider the model with @xmath48 , which happens to be the @xmath49 model which best - fits the current data ( for details of the fitting procedure see next section ) . first we plot the evolution of densities of different components in fig.rho . while one might think that the coupling between de and dm might help resolve the coincidence problem , we find that for these particular parameter values the evolution of different densities is rather similar to that of the standard @xmath50 model , i.e. the dm and de densities could still differ by orders of magnitude , so it does not seem to help much in this case . rather , as in other hde models , the fine tuning problem could be alleviated : the density of the dark energy is associated with the scale of horizon rather than the planck scale , due to the holographic principle . we have also investigated cases with larger coupling constant @xmath41 . we find that indeed , for the strong coupling case ( @xmath51 ) , the dark matter and dark energy could be made to evolve with comparable density for an extended range of redshifts , thus in some sense solving the `` coincidence problem '' . this is perhaps also true for any strongly coupled dark matter - dark energy model . however , for strong coupling the fit to observation is not good .        to investigate how the dynamics of the components are affected by coupling with the other component , one may compare the expansion term @xmath52 ( @xmath53 ) with the interaction term @xmath54 . if the interaction term is greater than or comparable with the expansion term of both of the components , the two components are strongly coupled and evolve together . as the universe expands , typically the density of the component would drop as @xmath55 , the interaction term often drops much more rapidly than the expansion term . when the expansion term of a component is much greater than the interaction term , that component evolves almost freely , and we may say the two components are decoupled . the coupling and decoupling effects between the two components could be asymmetric , as each has its own expansion term . it is possible that for one component the expansion term is greater , so it is decoupled from the other component , but for the other component , its expansion term is still smaller than or comparable with the interaction , thus it is still coupled to the other component .         in fig . [ compbest ] , we compare the expansion term of the dark energy and the dark matter with the interaction term as a function of redshift . in the left panel , we plot the magnitude of the dark energy expansion term @xmath56 , and in the right panel , we plot the magnitude of the dark matter expansion term @xmath57 , as red solid lines . the interaction term @xmath58 is plotted in both panels for comparison with the expansion terms . in the early universe , the dark matter takes a greater proportion of the total density , and for dark energy the interaction term is comparable with the expansion term , so the dark energy is moderately coupled to the dark matter . later , as the density of dark matter is diluted by expansion , the interaction term becomes negligible compared with the expansion term , and the dark energy is decoupled from dark matter . in the right panel of fig . [ compbest ] , we plot the dark matter expansion term @xmath59 and interaction term @xmath60 as a function of redshift @xmath61 . we see that throughout the cosmic history , for the dark matter the expansion term is always greater than the interaction term , so the dark matter is at most weakly coupled to dark energy . however , at lower redshifts , the difference becomes larger , and the dark matter become further decoupled from dark energy . such a decoupling history is quite typical for many interactions . however , there could also be some quite unexpected and interesting behavior in the presence of dark energy . in fig . [ couple ] , we make plots similar to fig . [ compbest ] but for the model with @xmath62 , and this model is chosen to illustrate such an unusual behavior . in this model , because @xmath63 , the dark energy density actually increases at later time . we see from the left panel of fig . [ couple ] , since initially the dark matter density is greater , that the interaction term is comparable to the expansion term for dark energy . however , later on the dark energy becomes dominant , and the interaction term falls below the expansion term of dark energy , so we may say that dark energy becomes decoupled . however , more interestingly , from the right panel of the figure , we see that initially the interaction term of the dark matter is smaller than the expansion term , so the dark matter is only weakly coupled to dark energy , but at later time , the interaction term is comparable with the expansion term of dark matter , so the dm becomes coupled to de . this is the reverse of the familiar decoupling process , which we shall call _ incoupling_. this strange and interesting behavior could only occur when a component is coupled to dark energy , whose density does not decrease during the cosmic expansion . as the de keeps being converted into dm , the dm evolution is greatly affected by de . this model does not yield a good fit to the observational data , but is only used to illustrate an interesting dynamical behavior in the presence of @xmath63 dark energy .    in fig.[hubble ] , we plot the evolution of @xmath64 for three choices of @xmath4 parameter , i.e. the best fit @xmath65 , and @xmath66 , while keeping @xmath67 . the hubble parameter for the @xmath68cdm model is also plotted for comparison . we also plot the equation of state of the dark energy for these models in fig . [ eos1 ] . on the left panel , curves for different @xmath4 but identical @xmath41 are plotted ; on the right panel , curves for identical @xmath4 but different @xmath41 are plotted . in all of these models , the phantom divide @xmath69 is either crossed , or will be crossed in the near future . the dark energy density would keep increasing , the interaction with dm would not slow it down , as the dark energy has become decoupled to the dm . eventually a `` big rip '' would occur . the current value of the equation of state @xmath70 for the best - fit model @xmath71 is around @xmath72 , which is within the @xmath73 cl region of @xmath70 in the wmap 5-year resultswmap5 , suggesting that this model is consistent with the current observations .        the value of @xmath41 in these models does not significantly affect the fate of the universe , but it does greatly affect the past evolution . this result is easy to understand ; as the dark energy dynamically decouples from dark matter , the interaction become negligible , regardless the precise value of coupling constant . the eventual fate of the universe in this interacting holographic dark energy model is qualitatively similar to those without interactions , though the details of the evolution would be different . thus , as in the case of non - interacting holographic models , for the fate of the universe , @xmath74 leads to phantom - like behavior and big rip@xcite , while @xmath75 leads to quintessence behavior . this is very different from the models discussed in ref .  @xcite . for further discussions on the fate of the universe and big rip singularity in our model , see appendix b.    we now consider the more general case of @xmath76 . we plot the equation of state in fig . [ eos2 ] , varying @xmath77 in the left panel while keeping @xmath78 , and varying @xmath79 in the right panel while keeping @xmath80 . the other parameters are kept as before . for @xmath81 , the equation of state for dark energy is slightly more negative . for @xmath82 , at higher redshift @xmath83 increases rapidly , and the coupling term would become very large , the whole system would be dominated by the evolution of matter . however , such a model would not really fit the low redshift observation , so it would not be relevant to our purpose . for large @xmath79 , we can also see that the equation of state of dark energy @xmath84 at high @xmath61 , but the variation is not as rapid as in the case of @xmath85 , and this is because the dark energy component does not increase as much as the dark matter component . still , this makes a poor fit to the observational data . anyway , for the case @xmath86 , eq.([eeos ] ) is reduced to @xmath87 substituting this equation into eq.([newrhol ] ) and making use of eq.(fredman - new1 ) , we obtain the following approximate equation for @xmath88 @xmath89 for which the solution will simply be @xmath90 where @xmath91 is a constant . this scaling behavior of dark energy density ( if @xmath86 ) is very similar to that of the chaplygin gas model @xcite , which at the high redshift also has the scaling behavior as @xmath92 while at the low redshift has the constant magnitude . however , in the present case it does not provide a good fit to the observations . now we fit the model parameters with the current observational data . we use three kinds of observations ; the first one is type ia supernovae ( sne ia ) which serves as the standard candle ; the second is the baryon acoustic oscillation ( bao ) in large scale structure of galaxy distribution , the third is the cosmic microwave background anisotropies , the latter two provide standard rulers at different redshifts . for the flat geometry we are considering , @xmath93 , where @xmath94 are the luminosity distance , angular diameter distance and comoving coordinate distance , respectively , and @xmath95 is given by @xmath96 which could be calculated for each model by solving eq . ( [ event case ] ) or eqs . ( [ motion3])-([motion4  ] ) .     versus @xmath4 for @xmath49.,width=326 ]    for the sne ia , we use the data set published by supernova cosmology project ( scp ) team recently @xcite . this data set contains 307 selected sne ia which includes several current widely used sne ia data sets , such as those collected by the hubble space telescope ( hst ) @xcite , the supernova legacy survey ( snls ) @xcite and the equation of state : supernovae trace cosmic expansion ( essence ) surveys@xcite . using the same analysis procedure and improved selection approach , all of the sub data sets are analyzed to obtain a consistent and high - quality  union  sne ia data set , which provides tighter and more reliable constraints for our models . for the bao data , we use the data given in percival _ et al_. @xcite and wmap team @xcite . we use the quantity @xmath97 which is determined by the sloan digital sky survey ( sdss ) ( @xmath98 ) and the two degree field galaxy redshift survey ( 2dfgrs ) ( @xmath99 ) in our constraints , where @xmath100 is the comoving sound horizon size at the decoupling epoch , and @xmath101 is the effective distance defined by eisenstein _ et al_. @xcite . since this quantity is just a ratio of two `` standard rulers '' , it could give a reliable measurement of the expansion history of the universe . for the cmb , we use the position of the first peak of the cosmic microwave background ( cmb ) angular power spectrum @xmath102 . the details of this technique were described in our earlier paper ref . @xcite . we use the markov chain monte carlo ( mcmc ) method to constrain the parameters with the combined data set of sn ia+bao+cmb ( for details of our implementation of this technique , see ref . @xcite ) . in the present model , the relevant parameters are the matter density @xmath103 , hubble constant @xmath104 , and the holographic dark energy parameter @xmath105 . we marginalize over @xmath104 and @xmath103 . for @xmath106 , we limit their range to @xmath107 . as discussed in the last section , greater values of @xmath106 would not yield a good fit to observation . for the fixed @xmath49 case , as can be seen in fig . [ lambda0 ] , the best - fit value and error on @xmath105 are @xmath108 the behavior of the best fit model has been discussed as the example given in the previous section . the models within the @xmath109 contour also exhibit similar behavior . now we let @xmath18 be free parameters , the pdf of @xmath110 and @xmath111 are shown in fig . [ alpha ] . the best - fit values for @xmath18 are @xmath112 and @xmath113 . if one demand @xmath18 to be integers for physical reasons , then @xmath114 is still the best fit . alternatively , a non - integer power may rise if the interaction is complicated . the best fit for @xmath4 and @xmath41 are @xmath115 and @xmath116 , as shown in fig . [ lambda ] , it is not too different from the @xmath49 case . the holographic dark energy ( hde ) model is a dynamical model of dark energy , motivated by the holographic principle , in which the event horizon of the space - time plays an important role in the dynamics . the microphysical basis of the model is not well understood , but the dynamics of cosmic expansion could be studied . observationally , the hde model is distinguishable from the @xmath68cdm , as its equation of state varies with time . the model could also be distinguished from other dynamical dark energy models . for example , in the hde model the dark energy behaves as the chaplygin gas at high redshift if @xmath117 ( see the discussion at the end of sec.3 ) , but for the low redshift region , the equation of state of the chaplypin gas approaches a constant , while for this model it continues to evolve . in properly constructed hde models such as the ones discussed in this paper , the dark energy component is subdominant at high redshifts compared with matter or radiation , so it would not significantly affect the cmb or large scale structure power spectrum , though at low redshift it could affect the growth of structure in detail . we plan to investigate the details of the structure growth in subsequent studies .    in the framework of the holographic dark energy model , we proposed a form of interaction between dark energy and dark matter which we considered to be natural and physically plausible . that is , the interaction rate is proportional to the product of powers of the dark energy and dark matter densities @xmath118 . for this type of interaction , we find that for the weak coupling case ( the dimensionless coupling constant @xmath119 , see eq . ( [ eq : lambda ] ) ) , the ultimate fate of the universe is similar to the case of non - interacting holographic dark energy , namely for @xmath74 the equation of state of the dark energy will cross @xmath69 , and behave like a phantom in the future , and the universe ends in a big rip . we also discovered the interesting phenomenon of _ incoupling _ , which occurs only when the dark matter is coupled to dark energy . this is the opposite of the familiar process of _ decoupling _ , the dark matter becomes increasingly coupled to the dark energy as universe expands . we then use the currently available observational data to constrain the parameters of this model . the data set we used include 307 selected sne ia published by the scp team @xcite , baryon acoustic oscillation from sdss and 2dfgrs @xcite , as well as the position of the first peak of the cmb angular power spectrum @xcite . we find that the best - fit form of interaction is @xmath120 . if one requires @xmath18 to be integers for reasons of fundamental physics , then the most favored interaction rate is proportional to the product of the density of dark matter and dark energy . our mcmc chain computation was performed on the supercomputing center of the chinese academy of sciences and the shanghai supercomputing center . this work is supported by the chinese academy of sciences under grant kjcx3-syw - n2 , by the ministry of science and technology national basic science program ( project 973 ) under grant no.2007cb815401 , and by the national science foundation of china under the distinguished young scholar grant 10525314 , and the key project grants no.10503010 . we give the detailed derivation of the dynamical equations of the interacting holographic dark energy model here . using the dimensionless notations given in sec.2 , we have @xmath121@xmath122@xmath123@xmath124@xmath125so @xmath126@xmath127 . \label{fredman - new2}\]]by substituting eq.([rho - derivative ] ) into eq . ( [ newrhol ] ) , we have @xmath128if @xmath1 is taken to be the particle horizon ( eq.([particle horizon ] ) ) or future event horizon ( eq.([event horizon ] ) ) , the derivative of @xmath1 is @xmath129where the upper ( lower ) sign represents particle horizon ( future event horizon ) . by substituting eq . ( [ l derivative ] ) and ( [ rhol ] ) into eq . ( [ equation of state ] ) , we have the equation of state in the particle horizon ( upper sign ) and future event horizon ( lower sign ) cases @xmath130 substituting eq . ( [ eos for ph , eh ] ) into ( [ fredman - new2 ] ) , we have @xmath131taking the derivative of eq . ( [ fredman - new1 ] ) with respect to @xmath132 , we have @xmath133substituting eq . ( [ new rhom ] ) into eq . ( [ intermediate2 ] ) , we obtain @xmath134subtracting equations ( [ intermediate3 ] ) from ( [ intermediate3 ] ) , we obtain the following simple dynamical equation : @xmath135 ,   \label{motion3}\]]eliminate @xmath136 in eq . ( [ intermediate1 ] ) using eq.(fredman - new1 ) , @xmath137and @xmath138to make eq . ( [ motion3 ] ) and ( [ motion4 ] ) easy to solve numerically , we change the variable to @xmath61 , @xmath139,\text { } \tilde{\rho}_{\lambda } ( 0)=1-\omega _ { m0 } ,   \label{motion3'}\]]@xmath140 ,   \notag \\ e(0 ) & = & 1 ,   \label{motion4''}\end{aligned}\]]where @xmath141and the upper ( lower ) sign in eq . ( [ motion3 ] ) and ( [ motion4  ] ) corresponds to the model with particle ( future event ) horizon . thus , eq . ( [ motion3 ] ) and ( [ motion4  ] ) are the complete system of equations . eqs.([eeos])-([motion4  ] ) reproduce eqs . ( [ event case ] ) . in the limit @xmath142 and @xmath143 the equations reduce to the non - interacting holographic dark energy model @xcite . the evolution of event horizon size @xmath1 and the corresponding dark energy density are shown in fig . [ diverge ] . as @xmath144 , the size of the event horizon shrinks to zero , and the dark energy density diverges , i.e. we are approaching a big rip singularity . there have been extensive discussions on the future singularity of the dark energy models , particularly for the case where the future universe is dominated by the phantom dark energy . @xcite investigated the appearance of the future singularity , and showed that the singularities could be classified into four types . the big rip is however , the only relevant singularity in the holographic dark energy model @xcite . ref @xcite showed that the coupling between dark matter and dark energy may help to solve the coincidence problem , but not the big rip ( type ) singularity issue . in particular , @xcite discussed the singularity problem in the holographic dark energy model , and proved that :    1 .   if holographic dark energy @xmath145 there exists a de sitter solution to the dynamic equation system of the hde . but if @xmath146 the de sitter solution does not exist . 2 .   for the de sitter solution ( @xmath147 ) , one should consider the back reaction of the quantum effects near the singularity . the contribution of the conformal anomaly as a back reaction near singularity could indeed moderate or even prevent the future singularity ( eqs.(53)-(61 ) in @xcite ) .    in our interacting hde model , the best - fit parameter @xmath148 , and for such model the fate of the universe is a big rip in the future . in the following , we investigate whether this conclusion is affected by the back reaction of quantum effects . since the future of the universe is dominated by dark energy , we will neglect all of the terms which contain @xmath15 and @xmath149 , as they will become subdominant when compared with the dark energy density @xmath150 the dynamical equations then simplify to @xmath151@xmath152@xmath153where we only consider the case for future event horizon . from eqs . ( rho1_sim)-([l_sim ] ) , the asymptotical behavior of @xmath64 , @xmath88 and @xmath154 can be obtained easily : @xmath155@xmath156@xmath157where @xmath158 and @xmath159 is the fractional hde density taken at the present day . from eqs . ( [ h_sim ] ) , the derivative terms are @xmath160@xmath161@xmath162 the conformal anomaly @xmath163 ( back reaction term ) has the following form @xcite : @xmath164thus in the hde model @xmath165where @xmath166 -3b^{\prime \prime } \left [ \frac{1}{c}\left ( \frac{1}{c}+1\right ) \left ( \frac{1}{c}-1\right ) \right ]    \notag \\ & & -b\left ( \frac{1}{c}-1\right ) \left [ 2\left ( \frac{1}{c}-1\right ) ^{2}+7\left ( \frac{1}{c}-1\right ) + 4\right ] . \label{fbc}\end{aligned}\]]with eqs.([ta ] ) and ( [ fbc ] ) , the energy density and pressure of the conformal anomaly @xmath167 are @xcite , @xmath168under this back reaction , the expansion rate is @xmath169 . \label{hubble2}\end{aligned}\]]since the energy density of back reaction is positive ( @xmath170 , the function @xmath171 therefore , @xmath172 increase monotonically in the future if @xmath173 thus , we prove that for the @xmath74 case , the back reaction effect from quantum corrections can not prohibit the occurrence of the big rip . 99 a.  g.  riess _ et al . _ [ supernova search team collaboration ] , astron . j.  * 116 * , 1009 ( 1998 ) [ astro - ph/9805201 ] ; s.  perlmutter _ et al . _ [ supernova cosmology project collaboration ] , astrophys . j.  * 517 * , 565 ( 1999 ) [ astro - ph/9812133 ] . m.  tegmark _ et al . _ [ sdss collaboration ] , phys . rev.d * 69 * , 103501 ( 2004 ) [ astro - ph/0310723 ] ; k. abazajian _ et al . _ [ sdss collaboration ] , astron . j.  * 128 * , 502 ( 2004 ) [ astro - ph/0403325 ] ; k.  abazajian _ et al . _ [ sdss collaboration ] , astron . j.  * 129 * , 1755 ( 2005 ) [ astro - ph/0410239 ] . a.  g.  cohen , d.  b.  kaplan and a.  e.  nelson , phys . lett .  * 82 * , 4971 ( 1999 ) [ hep - th/9803132 ] . p.  horava and d.  minic , phys . * 85 * , 1610 ( 2000 ) [ hep - th/0001145 ] ; s.  d.  thomas , phys . lett .  * 89 * , 081301 ( 2002 ) . s.  d.  h.  hsu , phys . b * 594 * , 13 ( 2004 ) [ hep - th/0403052 ] . m.  li , phys . b * 603 * , 1 ( 2004 ) [ hep - th/0403127 ] . g.  t hooft , gr - qc/9310026 ; l.  susskind , j.  math . phys .  * 36 * , 6377 ( 1995 ) [ hep - th/9409089 ] . j. d. bekenstein , phys . d * 7 * , 2333 ( 1973 ) ; j. d. bekenstein , phys . d * 9 * , 3292 ( 1974 ) ; j. d. bekenstein , phys . d * 23 , * 287 ( 1981 ) ; j. d. bekenstein , phys . d * 49 * , 1912 ( 1994 ) ; s. w. hawking , commun . * 43 , * 199 ( 1975 ) ; s. w. hawking , phys . d * 13 * , 191 ( 1976 ) . n. bilic , g. b. tupper , r. d. viollier , astro - ph/0207423 ; n. bilic , g. b. tupper , r. d. viollier , phys . b * 535 * ( 2002 ) 17 - 21 , [ astro - ph/0111325 ] ; m. c. bento , o. bertolami , a. a. sen , phys . d * 66 * ( 2002 ) 043507 , [ gr - qc/0202064 ] .
Abstract:

Answer: we consider the interaction between dark matter and dark energy in the framework of holographic dark energy , and propose a natural and physically plausible form of interaction , in which the interacting term is proportional to the product of the powers of the dark matter and dark energy densities . we investigate the cosmic evolution in such models . the impact of the coupling on the dark matter and dark energy components may be asymmetric . while the dark energy decouples from the dark matter at late time , just as other components of the cosmic fluid become decoupled as the universe expands , interestingly , the dark matter may actually become coupled to the dark energy at late time . we shall call such a phenomenon as _ incoupling_. we use the latest type ia supernovae data from the scp team , baryon acoustics oscillation data from sdss and 2df surveys , and the position of the first peak of the cmb angular power spectrum to constrain the model . we find that the interaction term which is proportional to the first power product of the dark energy and dark matter densities gives an excellent fit to the current data .


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: over the last fifteen years there has been a surge of studies in anomalous diffusion , characterized by the deviation of the mean squared displacement ( msd ) @xmath5 of a stochastic process with the probability density function @xmath1 to find the particle at position @xmath6 at time @xmath7 , from the linear time dependence @xmath8 of ordinary brownian motion @xcite . anomalous diffusion is usually characterized in terms of the power - law form @xmath9 with the anomalous diffusion coefficient @xmath10 of physical dimension @xmath11 and the anomalous diffusion exponent @xmath12 . depending on the value of @xmath12 we distinguish subdiffusion ( @xmath13 ) and superdiffusion ( @xmath14 ) @xcite .    since the milestone discoveries of superdiffusion in turbulence already in 1926 @xcite and of subdiffusion in amorphous semiconductors @xcite the recent vast increase of interest in anomalous diffusion is due to its discovery in numerous microscopic systems , in particular in biological contexts . the cytoplasm of biological cells is heavily crowded with various obstacles , including proteins , nucleic acids , ribosomes , the cytoskeleton , as well as internal membranes compartmentalizing the cell @xcite . diffusion of natural and artificial tracers in this complex environment is often subdiffusive . similar situations are encountered in cell membranes @xcite . the experimental evidence for subdiffusion in the crowded cytoplasm of living cells ranges from the motion of small labeled proteins @xcite over mrna molecules and chromosomal loci @xcite , lipid and insulin granules @xcite , virus particles @xcite , as well as chromosomal telomeres @xcite and cajal bodies @xcite inside the nucleus . subdiffusion was observed for membrane resident proteins experimentally @xcite and for membrane lipid molecules in computer simulations @xcite . in controlled _ in vitro _ experiments with artificial crowders , anomalous diffusion was consistently observed @xcite . on larger scales anomalous diffusion was observed , for instance , for the motion of bacteria in a biofilm @xcite . the observed subdiffusion was ascribed to various physical mechanisms @xcite . apart from the apparent transient anomalous diffusion caused by a crossover from free normal diffusion to the plateau value of the msd @xcite , typically , three main families of anomalous diffusion processes are considered : ( i ) diffusion in a fractal environment where dead ends and bottlenecks slow down the motion on all scales @xcite ; ( ii ) motion in a viscoelastic environment , in which the effective anomalous motion of the tracer particle in the correlated many - body environment shows long - ranged antipersistent motion @xcite . the latter process is associated with fractional brownian motion ( fbm ) and generalized langevin equation motion with a power - law memory form of the friction kernel @xcite . ( iii ) and continuous time random walk ( ctrw ) models , in which the moving particle is successively trapped by binding events to the environment or caging effects for waiting times @xmath15 distributed like a power - law @xmath16 with @xmath13 @xcite . all three mechanisms lead to the power - law msd ( [ msd ] ) and they were indeed identified as processes generating the motion of different tracers in different cellular environments @xcite or in colloidal systems _ in vitro _ @xcite . these three anomalous diffusion mechanisms characterized by anomalous diffusion with a constant in time generalized diffusivity were identified in experiments involving fairly large endogenous as well as artificial tracers . recent experiments on eukaryotic cells @xcite using considerably smaller tracer proteins sampling over much larger subvolumes of the cell indicated a systematic variation of the cytoplasmic diffusivity with the separation from the cell nucleus . these spatial diffusivity gradients are due partly to the non - uniform distribution of crowders in the cytoplasm . this distribution can non - trivially affect the diffusion of tracers of different sizes in the cell cytoplasm . _ in vitro _ , fast gradients of the diffusivity can be realized , for instance , via a local variation of the temperature in thermophoresis experiments @xcite , as sketched in fig . [ fig - model ] . the diffusion of brownian particles in explicit solvents with temperature gradients was recently studied by multi - particle collision dynamics @xcite . fluctuation - dissipation relations and spurious drift effects in systems with spatially varying friction coefficient were recently treated theoretically in ref . @xcite . on larger scales , the diffusion of water molecules monitored by diffusive magnetic resonance imaging in the brain white matter was demonstrated to be heterogeneous and strongly anisotropic @xcite . the anisotropy is due to the presence of some spatially - oriented structures in the tissue and obstacles which give rise to a tensorial character of the apparent diffusion coefficient . the existence of a population splitting into two pools of water molecules with slow and fast diffusivities was shown for brain white matter @xcite . finally , spatial heterogeneities are also abundant in the completely different context of anomalous diffusion in subsurface hydrology @xcite .    . the diffusivity of the random walker is coupled to the spatial variation of the temperature of the environment , as symbolized by the speed of the runner.,width=302 ]    in the present study we examine the effects of the strength of the diffusivity gradient as defined by the scaling exponent @xmath4 and of the initial particle position @xmath17 in heterogeneous diffusion processes ( hdps ) with space - dependent diffusion coefficient @xmath18 . we pay particular attention to a phenomenon , that recently received considerable attention for its immediate relevance to the surging field of single particle tracking experiments in microscopic systems , namely , the so - called _ weak ergodicity breaking_. this is the distinct disparity between physical observables depending on whether they are evaluated in the conventional ensemble sense or , from measured time series @xmath19 of the particle position , in terms of time averages @xcite . we find that the hdp process gives rise to weakly non - ergodic behavior with a pronounced amplitude scatter of the time averaged msd of individual trajectories . we obtain details of the distribution of this scatter as well as the frequently used ergodicity breaking and non - gaussianity parameters . our analysis uncovers remarkable similarities of the hdp process with those of ctrw motion . in particular , we study the behavior of the hdp at the critical value of the scaling exponent of the diffusivity . the paper is organized as follows . we introduce the hdp model with @xmath6-dependent diffusivity in sec . [ sec - model ] . the main results for the evolution of the msd , the time averaged msd , the probability density function , the ergodicity breaking and non - gaussianity parameters in the whole range of the model parameters are then presented in sec . [ sec - main - results ] . we discuss our results and point out the directions for future research in sec . [ sec - discussion ] . hdps are defined in terms of the multiplicative yet markovian langevin equation @xcite @xmath20 where @xmath18 is the position - dependent diffusion coefficient and @xmath21 represents white gaussian noise . in what follows we concentrate on the power - law form @xmath22 for the diffusivity , where the amplitude @xmath23 has dimension @xmath24 . logarithmic and exponential forms for @xmath18 were considered in ref . the offset @xmath25 in eq . ( [ eq - dx ] ) avoids either divergencies of @xmath18 ( @xmath26 ) or stalling of the particle ( @xmath27 ) around @xmath28 in the simulations . in the following calculations we use the scaling form @xmath29 . we interpret the langevin equation ( [ lang ] ) in the stratonovich sense @xcite . the msd following from the stochastic equation ( [ lang ] ) with diffusivity ( [ eq - dx ] ) takes on the power - law form @xcite @xmath30 where we introduced the scaling exponent @xmath31 which denotes superdiffusion for @xmath32 and subdiffusion for @xmath26 . for @xmath27 the diffusivity grows away from the origin , leading to a progressive acceleration of the particle as it ventures into more distant regions from the origin , and vice versa for @xmath26 . in the special case @xmath3 , the theory developed in ref . @xcite breaks down as the scaling exponent ( [ pexp ] ) diverges . in refs . @xcite an exponential growth of the msd was found , see the discussion below . for even larger values of @xmath4 the anomalous diffusion exponent @xmath33 becomes negative , i.e. , we observe a strong localization , see the discussion below . the probability density function ( pdf ) of the hdp is given by the stretched or compressed gaussian @xcite @xmath34 for positive @xmath4 it exhibits a cusp , while for negative @xmath4 the pdf features a dip to zero at the origin .    in single particle tracking experiments one measures the time series @xmath19 of the particle position for a time span @xmath35 . it is usually evaluated in terms of the time averaged msd @xcite @xmath36 ^ 2dt , \label{eq - tamsd}\ ] ] where @xmath37 is the lag time . brownian motion is ergodic and for sufficiently long measurement times @xmath35 the equality @xmath38 holds @xcite . anomalous diffusion ( [ msd ] ) described by fbm or fractional langevin equation motion is asymptotically ergodic @xcite but may feature transiently non - ergodic behavior @xcite in confinement as well as transient aging @xcite , the explicit dependence on the time span elapsing between initial preparation of the system and start of the recording of the position , @xmath19 . ctrw processes with diverging time scales of the waiting time distribution @xmath39 show weakly non - ergodic behavior for all @xmath37 . namely , despite the scaling ( [ msd ] ) of the msd , @xmath40 scales _ linearly _ with @xmath37 . more precisely , if we average over sufficiently many trajectories , the quantity @xmath41 for ctrws scales like @xmath42 @xcite . this linear dependence on @xmath37 is preserved for aging ctrws @xcite . apart from this linear @xmath37-dependence we also observe the dependence on the process time @xmath35 , a signature of aging @xcite . interestingly , the hdp with diffusivity ( [ eq - dx ] ) displays weak ergodicity breaking of the form @xcite @xmath43 this behavior is analogous to that of scale - free ctrw motion , despite of the fact that the increment correlation function of hdps are ( anti)persistent in analogy to the ergodic fbm @xcite . we note that also other processes such as aging and correlated ctrws @xcite as well as scaled brownian motion with time - dependent diffusivity @xcite exhibit the duality between ( [ msd ] ) and a linear @xmath37-dependence of @xmath40 . ( [ eq - tamsd - scaling ] ) shows the @xmath44-scaling as function of the process time @xmath35 . for subdiffusion with @xmath45 , that is , the effective diffusivity of the process decays over time , as the particle ventures into low - diffusivity areas . conversely , for @xmath46 the diffusivity increases over time , the particle discovers areas with increasingly higher @xmath18 .    in the current paper we perform a detailed analysis of the msd and the time averaged msd in the entire range of @xmath4 , including the critical value @xmath3 . we are particularly interested in the extent of the weakly non - ergodic behavior , especially the fluctuations of the time averaged msd around the mean value @xmath47 characterized by the ergodicity breaking parameter . moreover we analyze the non - gaussianity of the process . in our analysis we study effects of the scaling exponent @xmath4 of @xmath18 as well as the initial position of the particle . the latter is known to affect the time - scales at which anomalous diffusion becomes significant @xcite . we combine analytical and numerical approaches . the simulations scheme for hdps was introduced in ref . in this section we start with the analysis of the msd and the time averaged msd with its amplitude fluctuations , the latter being quantified by the corresponding scatter distribution . we then analyze the ergodicity breaking and non - gaussianity parameters . finally , we study the pdf of the function that quantifies the degree of particle dispersion .      fig . [ fig - msd - tamsd - alfa ] shows the results from computer simulations for the msd @xmath48 and the time averaged msd @xmath40 from individual realizations along with the average @xmath49 taken over all single time traces . the values for the scaling exponent @xmath4 of the diffusivity ( [ eq - dx ] ) studied in fig . [ fig - msd - tamsd - alfa ] cover both the sub- and superdiffusive domains and include , in particular , the critical value @xmath3 where the scaling exponent @xmath33 of the msd ( [ msd ] ) diverges . for each @xmath4 we show @xmath50 single trajectories . in all cases , apart from the critical point , the scaling of the msd @xmath51 and both individual ( @xmath40 ) and mean ( @xmath47 ) time averaged msds agree well with the expected analytical behavior shown by the dashed lines : the scaling exponent of the msd varies consistently with @xmath4 , while that of the time averaged msd remains unity throughout . the initial discrepancy between theory and msds is due to the choice for the initial position @xmath52 , whose influence relaxes on a time scale depending on both @xmath17 and @xmath4 . the deviations of individual time traces @xmath40 at long lag times from the predicted behavior is due to unavoidable , bad statistics when the lag time gets close to the overall length @xmath35 of the time series . irreproducibility of time averages of physical observables such as the msd is an intrinsic property of weakly non - ergodic processes @xcite . for ctrw processes with scale - free waiting time distribution individual traces @xmath19 contain different , few dominating waiting time events that cause the amplitude scatter between different trajectories or , in other words , fluctuations of the apparent effective diffusion constant . this phenomenon occurs no matter how long the measurement time @xmath35 is taken . for hdps the scatter is due to the difference in the extent of excursions to regions of significantly different diffusivity .    the amplitude scatter between individual realizations @xmath40 in our hdp shown in fig . [ fig - msd - tamsd - alfa ] varies significantly with the value of the scaling exponent @xmath4 : away from the brownian value @xmath53 in both subdiffusive ( @xmath26 ) and superdiffusive ( @xmath27 ) cases the fluctuations of @xmath40 become more pronounced when the modulus of @xmath4 increases . at the critical point the fluctuations of @xmath40 appear to diverge , while beyond this critical point , the fluctuations decrease again . moreover , a population splitting in a faster ( steeper slope of @xmath54 ) and slower ( shallower slope ) fraction of trajectories appears @xcite , see especially the panel for @xmath55 . in terms of the dimensionless variable @xmath56 the amplitude scatter distribution @xmath57 reflects the randomness of individual time averages of the msd . for a sub- and a superdiffusive @xmath4 it was analyzed in ref . as shown here for a whole spectrum of @xmath4 values , there is a clear trend towards extreme fluctuations at the critical point @xmath3 , but even for considerably smaller values such as @xmath58 the fluctuations around the mean @xmath59 are significant . the width of the fluctuations of @xmath40 in each panel varies only marginally with the lag time @xmath37 , apart from the behavior at @xmath60 . as studied in ref . @xcite , the relative amplitude scatter distribution @xmath57 can be fitted with the generalized gamma - distribution . in particular , it tends to zero at @xmath61 , in contrast to subdiffusive ctrw processes , for which @xmath62 is always positive , indicating completely stalled trajectories @xcite . we note that when the exponent @xmath4 approaches the critical value @xmath3 , the number of steps @xmath63 necessary to approach the theoretically predicted asymptote ( [ eq - msd - scaling ] ) increases significantly . thus , as seen from fig . [ fig - msd - tamsd - alfa ] , for @xmath52 only few simulation steps , @xmath64 are needed for negative @xmath4 with larger modulus . it increases to @xmath65 steps for @xmath66 , @xmath67 for @xmath68 , and already @xmath69 for @xmath58 . once the anomalous diffusion exponent @xmath33 becomes negative , that is , for values of @xmath4 larger than the critical value @xmath3 , the msd ( [ eq - msd - scaling ] ) becomes a decreasing function of time . this agrees with our simulations if @xmath17 is chosen sufficiently large to enable the relaxation to the theoretical asymptote . for instance , for the extreme value @xmath70 and values of the initial position of @xmath71 and above the msd indeed follows the theoretical prediction @xmath72 ( not shown ) . in this region @xmath73 the diffusivity grows fast away from the origin the decreasing msd corresponds to the localization of particles in regions of slow diffusivity , see the detailed discussion in sec . [ sec - spread ] below .    in the limit @xmath74 , due to the huge spread one or few extremely large amplitudes @xmath40 of the time averaged msd may substantially affect the mean @xmath47 . the msd in this limit follows an exponential growth @xmath75 , as indicated by the dashed curve in the panel for @xmath3 in fig . [ fig - msd - tamsd - alfa ] . such a fast increase of the msd is consistent with the divergence of the scaling exponent @xmath33 as function of @xmath4 in eq . ( [ eq - msd - scaling ] ) as well as with the exponential msd growth predicted for a parabolically space - varying diffusivity in refs . this property can be straightforwardly inferred from the diffusion equation , @xmath76.\ ] ] multiplying both sides with @xmath77 and integrating over @xmath6 one arrives at @xcite @xmath78.\ ] ] this also rationalizes the observation that the scatter of @xmath40 is maximal for @xmath3 as the particles perform extremely far - reaching excursions relative to other @xmath4-values .      as introduced in refs . @xcite the fluctuations of the time averaged msd @xmath40 can be quantified by their variance , the ergodicity breaking parameter @xmath79 when @xmath80 , it means that the process is perfectly reproducible and all time averages over sufficiently long trajectories give the same value . this case corresponds to the sharp scatter distribution @xmath81 @xcite . for the canonical brownian motion the ergodicity breaking parameter reaches the zero value in the form @xmath82 at finite ratio @xmath83 . weakly non - ergodic processes have a positive value of @xmath84 . an alternative , weaker condition for ergodicity is when the msd ( [ msd ] ) and the time averaged msd ( [ eatamsd ] ) coincide . to measure the relative deviations from ergodicity , the @xmath85 parameter was introduced @xcite @xmath86 for ergodic dynamics its value is unity .     on the lag time @xmath37 . the brownian asymptote @xmath87 is shown as the dashed curve in panel a ) . the asymptotes ( [ eq - eb2 ] ) are shown in panel b ) as dashed lines of the corresponding color . the initial position of the particle was @xmath52 and @xmath88 traces were used for each value of @xmath4 , all other parameters are the same as those used in fig . [ fig - msd - tamsd - alfa ] . the symbols in both panels correspond to the same parameters.,width=302 ]    , eq .  ( [ pexp ] ) , of @xmath89 in panels a ) and the magnification in panel b ) as well as @xmath90 in panel c ) . the theoretical predictions from ref . @xcite correspond to the dashed black curve in panels a ) and b ) . the ergodicity breaking parameter of the subdiffusive ctrw given by eq .  ( [ eq - eb - ctrw ] ) is shown by the dot - dashed red curve in panel b ) . to compute each point in the graphs for @xmath50 traces of @xmath91 steps takes some two hours on a standard 3 ghz core workstation . different sets of the shown points correspond to varying particle initial conditions @xmath17 indicated by different colors . note the alternative scale for the power exponent @xmath92 on the top axes . parameters are the same as in fig . [ fig - msd - tamsd - alfa].,width=264 ]    to calculate the @xmath84 parameter analytically is not always an easy task . to see this , note that the msd follows from the change of the stochastic variable to the standard wiener process in the form @xmath93^{-1/2}$ ] for the time averaged msd @xmath40 the calculation is already more complicated , as it involves the two - point position correlation function . the latter is expressed via fox @xmath94-functions for the hdp @xcite . the analytical derivation of the ergodicity breaking parameter @xmath84 , however , requires fourth - order moments , whose calculation is a formidable task . so far only approximate methods are known for hdps @xcite . for subdiffusive ctrw processes it is possible to obtain the @xmath84 parameter more easily from the conjectured and numerically proven equivalence @xmath95 of @xmath96 with the ratio of the number of steps @xmath97 in an individual realization and the average @xmath98 @xcite . as for the limiting distribution for @xmath99 is known @xcite , this allows straightforward calculation of @xmath57 and its moments @xcite . for the multiplicative process studied here such a scheme does not work . we also note that another parameter involving the fourth moment of the particle displacement is the non - gaussianity parameter @xmath100 @xcite discussed below . due to the lack of an analytical theory , a major reason for the current simulations study is to explore the behavior of these two parameters in the whole range of the model parameters .    figs . [ fig - eb - lag ] and [ fig - eb - p ] display the dependence of the ergodicity breaking parameters on the lag time @xmath37 , the initial position @xmath17 , as well as the scaling exponent @xmath33 . we find that for standard brownian motion with @xmath101 the ergodicity breaking parameter eb follows the known asymptote @xmath102 from above . concurrently @xmath103 as expected for ergodic motion , except for very small @xmath37 values because of the initial relaxation of the influence of the initial position @xmath17 . as we depart from the value @xmath53 of brownian motion , the magnitude of @xmath84 grows and its functional dependence on the lag time becomes less pronounced , see fig . [ fig - eb - lag]a . approaching the critical point @xmath104 , due to the huge fluctuations of @xmath40 and the ensuing possibility of extreme events the ergodicity breaking parameter also explicitly depends on the number of traces @xmath105 used for the averaging and reaches values of @xmath106 and higher . on both sides of the critical point , corresponding to large positive or negative values of @xmath33 , the values of the parameter @xmath107 approach one another , as seen in fig . [ fig - eb - p]a , while @xmath108 exhibits a jump , fig .  [ fig - eb - p]b . the ratio of the different msds given by @xmath85 also depends on the power exponent @xmath4 and the initial value @xmath17 . as anticipated already from fig . [ fig - msd - tamsd - alfa ] , for the chosen initial condition @xmath52 the magnitude of @xmath85 decreases as @xmath4 gets progressively negative , while @xmath85 grows as @xmath4 increases from 0 to 2 , reaching very large values @xmath109 at @xmath3 . this is illustrated in fig . [ fig - eb - lag]b which is also consistent with the scaling @xmath110 with the lag time at different @xmath111 values . because of the initial condition @xmath17 , similar to our statements for @xmath47 this asymptote is approached later when @xmath74 .    in fig . [ fig - eb - p ] we also analyze the effects of the initial position @xmath17 and show the dependence of the ergodicity breaking parameters for long traces or short lag times , i.e. , when @xmath112 which is practically equivalent to @xmath113 , versus the anomalous diffusion exponent @xmath33 and the scaling exponent @xmath4 . we observe that in the region @xmath46 the values of @xmath107 obtained from simulations is in good agreement with the analytical estimate @xmath113 from ref . @xcite , represented by the dashed curve in fig .  [ fig - eb - p]a and b. these approximations are based on the asymptotic scaling of @xmath18 and do not consider the effect of the initial positions @xmath17 , assumed to be relaxed in the relevant @xmath114 limit . to see their actual impact on the dynamics we study the ergodicity breaking parameters numerically . we observe that the value of @xmath107 grows from the small pre - asymptotic brownian value at @xmath115 to progressively larger values at larger modulus of @xmath33 , as @xmath4 approaches the critical value @xmath3 from below and above . the functional dependence of @xmath116 obtained from our simulations follows quite well the predictions from the approximate calculation of @xmath84 @xcite , particularly in the range @xmath117 , see fig . [ fig - eb - p]b . the deviations for even more pronounced variation of @xmath18 as @xmath4 tends to 2 are likely due to insufficient statistics . as we show in fig . [ fig - eb - p]a by the circles and triangles , the value @xmath118 reveals measurable deviations for @xmath119 as compared to @xmath50 traces used for the averaging , while for less extreme @xmath4 values the two sets yield nearly identical results . this is a further indication towards a divergence of the fluctuations at the critical point . namely , the chance to find an even larger amplitude of @xmath40 in a given realization grows with the number of simulated trajectories .    as already anticipated in ref . @xcite the hdp process considered here is in some aspects reminiscent of subdiffusive ctrw processes with scale - free , power - law waiting time distributions with regard to the scaling of the msd and the time averaged msd , compare also sec . [ sec - discussion ] . apparently , the ergodicity breaking properties of these two processes also reveal similar features : the fluctuations of @xmath40 increase when the anomalous diffusion exponent @xmath33 becomes successively smaller than the brownian value 1 . indeed , as shown in fig . [ fig - eb - p]b for the relevant region @xmath45 there exists even a close agreement of the quantitative behavior of ctrw and hdp . we show the ergodicity breaking parameter obtained for the ctrw process @xcite @xmath120 where @xmath13 is the exponent in the pdf @xmath16 of the waiting times @xmath15 @xcite . for the ctrw process the exponent @xmath12 also occurs in the msd ( [ msd ] ) . we note that in contrast to hpds in superdiffusive ctrw processes the nonergodic behavior is merely ultraweak , effecting a different prefactor in the time averaged msd compared to the msd @xcite .    for @xmath121 we find that @xmath122 for initial positions @xmath17 close to the origin , in agreement with the results presented in ref . @xcite . for @xmath73 or @xmath123 the analytical model of ref . @xcite no longer applies but the simulations yield another region of growth for the value of @xmath107 . note that @xmath107 also depends on the trace length @xmath35 ( not shown ) , as discussed for two - dimensional hdps @xcite . as we show in fig . [ fig - eb - p]c , the magnitude of @xmath108 strongly varies with @xmath17 . in the region @xmath124 the ratio of the msds grows as @xmath125 indicated by the dashed asymptote in fig . [ fig - eb - p]c , while for large negative @xmath33 we observe @xmath126 . we finally note that @xmath108 reveals a much weaker dependence on the number @xmath105 of simulated traces , see the circles and triangles in fig . [ fig - eb - p]c . the reason is that @xmath85 involves only the second moment of the time averaged msd , while the @xmath84 is defined in terms of the fourth order moment , which is more sensitive to large variations .       as function of the lag time , plotted for the parameters of fig . [ fig - eb - lag ] with @xmath50 trajectories for each @xmath4 value . the approximately zero - valued trace @xmath127 for brownian motion ( @xmath53 ) is only shown partially . color coding is the same as in fig . [ fig - eb - lag ] . * b ) * the value of @xmath128 evaluated for varying power exponent @xmath4 and initial position @xmath17 . the color scheme is the same as in fig . [ fig - eb - p].,width=317 ]    the non - gaussianity parameter @xmath100 is a sensitive experimental indicator that often enables one to distinguish the type of particular diffusion processes observed in single - particle tracking experiments @xcite . it is related to the stationarity of increments of the diffusion process and involves the fourth moment along the time - averaged trajectory . for the diffusion process in an embedding space of dimension @xmath129 it is defined via the experimentally relevant time averaged quantities as @xcite @xmath130 where , in analogy to eq . ( [ eq - tamsd ] ) , the fourth moment is defined via @xmath131 ^ 4~dt.\ ] ] for brownian and fractional brownian motion , both gaussian processes , one consistently finds @xmath132 . for diffusion processes revealing a transient anomalous diffusion behavior , this parameter deviates substantially from zero , see , e.g. , the discussion in ref . @xcite . we systematically examine the behavior of @xmath127 in fig . [ fig - non - gauss]a . corroborating the results for the ergodicity breaking parameter , for brownian motion we obtain approximately zero values . as the power @xmath4 in @xmath133 deviates from zero , the non - gaussianity parameter reveals a rich behavior as a function of the lag time @xmath37 and the initial position @xmath17 . more specifically , in the region @xmath26 the non - gaussianity parameter progressively grows and reaches considerably large values for large @xmath2 , see the curve for @xmath134 in fig . [ fig - non - gauss]a . the value @xmath127 systematically decreases with the lag time @xmath37 along the trace for negative @xmath4 . in the region @xmath135 the value of @xmath100 also grows with @xmath4 but , in contrast to the case @xmath26 , the function @xmath127 stays rather constant with @xmath37 . as we approach the critical @xmath3 the non - gaussianity parameter reaches very high values . in the region @xmath73 the value of @xmath127 decreases again . these features of the functional behavior of @xmath136 correspond with the properties of @xmath137 , compare the curves in fig . [ fig - non - gauss]a and [ fig - eb - lag]a . note that similar to @xmath84 the statistics required to obtain a smooth curve for @xmath127 strongly depends on @xmath4 : the result from @xmath50 traces in fig . [ fig - non - gauss]a acquires pronouncedly higher fluctuations when we approach the critical value @xmath3 . we find that variations in the initial positions @xmath17 only have a marginal effect on @xmath138 . in fig . [ fig - non - gauss]b we illustrate the range of @xmath100 values for exponents in the range @xmath139 . even for the wide range from @xmath71 to @xmath140 the value of @xmath128 does not reveal any appreciable variation , in stark contrast to the strong @xmath17-dependence of @xmath141 . this indicates that for the hdp the non - gaussianity parameter is more robust than @xmath84 with respect to the choice of the initial condition . in fig . [ fig - non - gauss]b the dashed line indicates the proportionality to @xmath142 , which nicely matches the measured shape of @xmath138 . we finally analyze the time evolution of the pdf of the particle position in fig . [ fig - spread ] for three different values of the scaling exponent @xmath4 , distinguishing superdiffusion ( top row for @xmath68 , i.e. , ballistic motion with @xmath143 ) , normal diffusion ( middle row for @xmath53 and @xmath115 ) , and subdiffusion ( bottom row for @xmath144 , i.e. , @xmath145 ) . in each row the leftmost panel represents a very early evolution of the pdf close to the initial condition , while in the rightmost panel the pdf in some of the cases has almost reached the diffusion limit , in which the analytical asymptote ( [ pdf - hdp ] ) is valid . we observe that for hdps with positive @xmath4 , i.e. , when the diffusivity @xmath18 grows with the modulus of the position , the pdf for particles with an initial position far away from the origin , an asymmetric shape of the pdf is effected . namely , they progressively move towards regions of small diffusivity and accumulate there . this focusing due to the quenched nature of the diffusivity erases any memory of the initial condition and the common asymptote ( [ pdf - hdp ] ) of stretched gaussian shape ( @xmath135 ) is approached . in the opposite case with negative @xmath4 the focusing of particles in lower diffusivity regions applies again , albeit at longer times than shown here . this time , however , instead of the cusp at the origin for positive @xmath4 , the pdf drops down to zero at the origin and acquires the bimodal shape predicted by ( [ pdf - hdp ] ) , a compressed gaussian . for the brownian case with vanishing @xmath4 , no focusing takes place . in the final panel the distribution is already so broad that the individual , shifted pdfs appear on top of each other . similarly to the evolution of the msd shown in fig . [ fig - msd - tamsd - alfa ] , the relaxation time required for the system of diffusing particles to approach this long - time limit , @xmath1 depends on the power @xmath4 of the diffusivity in eq . ( [ eq - dx ] ) . for more negative @xmath4 values the traces are not yet relaxed to the theoretical long - time shape even after @xmath91 steps . the particles starting at larger @xmath17 remain trapped for the whole length of the simulation . conversely , when the diffusivity increases fast away from the origin ( moderate positive @xmath4 values ) , the equilibration to the long - time pdf is relatively fast , see the panel for @xmath68 in fig . [ fig - spread ] . for hpd processes whose diffusivity @xmath18 varies in power - law form with the distance @xmath6 from the origin and which give rise to anomalous diffusion we analyzed in detail the weakly non - ergodic behavior for varying power exponents @xmath4 and initial positions @xmath17 of the particles . in particular , we examined the functional dependencies and magnitudes of the variation of the ensemble and time averaged characteristics of such hdps with these parameters . the fluctuations of the amplitude of the time averaged msd of individual time traces were shown to systematically grow with increasing departure of the anomalous diffusion exponent @xmath33 from the brownian value @xmath115 , corresponding to the scaling exponent @xmath53 in the power - law form for @xmath18 . the fluctuations increase dramatically when the scaling exponent approaches its critical value @xmath3 . at that critical point we corroborated the turnover from the power - law scaling in time of the msd to an exponential growth . beyond the critical point the fluctuations of the time averaged msd become smaller again . concurrently we observe a distinct population splitting into a faster and slower fraction of time averaged msds . we paid particular emphasis on several parameters used to classify the departure from ergodic behavior , namely the two ergodicity breaking parameters @xmath84 and @xmath85 as well as the non - gaussianity parameter @xmath100 . while @xmath146 is simply defined as the ratio of the mean time averaged msd versus the msd , both @xmath84 and @xmath100 are based on the fourth order moments of the particle position and are known analytically only from approximate theories . a detailed numerical analysis of their properties was therefore used to obtain more concrete information on their behavior for different values of the scaling exponent @xmath4 and the initial particle position @xmath17 in the heterogeneous environment @xcite . within the analyzed parameter range we find that the behavior of @xmath147 is indeed in agreement with the heuristic theoretical analysis from ref . the parameter portraits for the ergodicity breaking and non - gaussianity parameters obtained from our numerical analysis will be useful for actual data evaluation . we also explored the detailed behavior of the hdp dynamics at the critical value @xmath3 and its vicinity , in particular with respect to the dramatic values reached by @xmath84 reflecting the dramatic fluctuations of the time averaged msd . a systematic numerical analysis of the particle pdf for varying initial conditions @xmath17 sheds additional light on hdps with different exponent @xmath4 , in particular , the visualization of the particle focusing in low diffusivity regions .    in hdps the non - ergodic behavior arises due to the heterogeneity of the environment . physically , this represents a space - dependent mobility or temperature . it is not a property of the particle , and each time the particle revisits a given point @xmath6 in space it has the same diffusivity @xmath18 . in a random walk sense , this scenario could also be translated into a local dependence of the waiting time for a jump event . in that interpretation the hdp corresponds to a motion in a quenched energy landscape @xcite , albeit a deterministic ( in contrast to random ) one . as such , the process is intrinsically different from renewal ctrw processes , which correspond to the motion in an annealed environment @xcite . nevertheless we observed that subdiffusive hdps share a number of features with subdiffusive ctrws with scale - free , power - law waiting time distribution . these features include the scaling laws for the ensemble and time averaged msds , and the form of the ergodicity breaking parameter @xmath85 ( in the relevant range @xmath45 ) . from the sole analysis of the msds @xmath51 and @xmath40 as well as the ergodicity breaking parameters , a significant distinction between ctrw and hdp is therefore difficult . yet there exist some crucial differences between the hdp and ctrw processes , that can be measured . thus , in hdps with @xmath45 the distribution @xmath57 of the relative amplitude @xmath40 of individual realizations decays to zero at @xmath61 , in contrast to the ctrw s finite fraction of immobile particles reflected in the finite value @xmath148 . indeed the scatter distribution @xmath57 was previously advocated as a good diagnosis tool for different stochastic processes @xcite , complementing other stochastic analysis methods @xcite . hdps are physically meaningful alternatives to the standard stochastic models for anomalous diffusion processes , namely , the ctrw process with long - tailed waiting time distribution , fbm and fractional langevin equation motion based on gaussian yet long - ranged in time correlations , diffusion in fractal environments , or their combinations @xcite . hpds are weakly non - ergodic , sharing some yet not all properties of ctrw processes . it is therefore important to explore the properties of hdps in more detail in the future as well as to develop more sophisticated tools to interpret measured time series and with confidence identify the underlying stochastic mechanism . in particular , we will study the properties of confined and aged hdps as well as spatially varying diffusivities with a random component , and the coupling of hdps with active motion . i. bronstein , y. israel , e. kepten , s. mai , y. shav - tal , e. barkai , and y. garini , phys . lett . * 103 * , 018102 ( 2009 ) ; k. burnecki , e. kepten , j. janczura , i. bronshtein , y. garini , and a. weron , biophys . j. * 103 * , 1839 ( 2012 ) . g. r. kneller , k. baczynski , and m. pasenkiewicz - gierula , j. chem . phys . * 135 * , 141105 ( 2011 ) ; j .- jeon , h. martinez - seara monne , m. javanainen , and r. metzler , phys . * 109 * , 188103 ( 2012 ) ; m. javanainen , h. hammaren , l. monticelly , j .- h . jeon , m. s. miettinen , h. martinez - seara monne , r. metzler , and i. vattulainen , faraday disc . * 161 * , 397 ( 2013 ) . s. c. kou , ann . * 2 * , 501 ( 2008 ) ; e. lutz , phys . e * 64 * , 051106 ( 2001 ) ; r. kupferman , j. stat . phys . * 114 * , 291 ( 2004 ) ; w. min , g. luo , b. j. cherayil , s. c. kou , and x. s. xie , phys . lett . * 94 * , 198302 ( 2005 ) .      q. xu , l. feng , r. sha , n. c. seeman , and p. m. chaikin , phys . 106 * , 228102 ( 2011 ) ; i. y. wong , m. l. gardel , d. r. reichman , e. r. weeks , m. t. valentine , a. r. bausch , and d. a. weitz , _ ibid . _ * 92 * , 178101 ( 2004 ) . bouchaud , j. phys . i ( paris ) * 2 * , 1705 ( 1992 ) ; g. bel and e. barkai , phys . lett . * 94 * , 240602 ( 2005 ) ; a. rebenshtok and e. barkai , _ ibid . _ * 99 * , 210601 ( 2007 ) ; m. a. lomholt , i. m. zaid , and r. metzler , phys . lett . * 98 * , 200603 ( 2007 ) . v. tejedor and r. metzler , j. phys . a * 43 * , 082002 ( 2010 ) ; m. magdziarz , r. metzler , w. szczotka , and p. zebrowski , phys . rev . e * 85 * , 051103 ( 2012 ) . see also the alternative approach in a. v. chechkin , m. hofmann , and i. m. sokolov , phys . e * 80 * , 031112 ( 2009 ) ; j. h. p. schulz , a. v. chechkin , and r. metzler , j. phys . a. * 46 * , 475001 ( 2013 ) . g. zumofen and j. klafter , physica d * 69 * , 436 ( 1993 ) ; a. godec and r. metzler , phys . e * 88 * , 012116 ( 2013 ) ; d. froemberg and e. barkai , _ ibid . _ * 87 * , 030104(r ) ( 2013 ) ; _ ibid . _ * 88 * , 024101 ( 2013 ) ; t. akimoto , phys . * 108 * , 164101 ( 2012 ) . m. magdziarz , a. weron , k. burnecki , and j. klafter , phys . . lett . * 103 * , 180602 ( 2009 ) ; m. magdziarz and j. klafter , phys . e * 82 * , 011129 ( 2010 ) ; y. meroz , i. eliazar , and j. klafter , j. phys . a * 42 * , 434012 ( 2009 ) ; y. meroz , i. m. sokolov , and j. klafter , phys . rev . lett . * 110 * , 090601 ( 2013 ) ; s. condamin , v. tejedor , r. voituriez , o. bnichou , and j. klafter , proc . usa * 105 * , 5675 ( 2008 ) ; v. tejedor , o. bnichou , r. voituriez , r. jungmann , f. simmel , c. selhuber - unkel , l. oddershede , and r. metzler , biophys . j. * 98 * , 1364 ( 2010 ) ; m. bauer , r. valiullin , g. radons , and j. kaerger , j. chem . phys . * 135 * , 144118 ( 2011 ) ; j .- h . jeon , e. barkai , and r. metzler , j. chem . phys . * 139 * , 121916 ( 2013 ) .
Abstract:

Answer: we study the stochastic behavior of heterogeneous diffusion processes with the power - law dependence @xmath0 of the generalized diffusion coefficient encompassing sub- and superdiffusive anomalous diffusion . based on statistical measures such as the amplitude scatter of the time averaged mean squared displacement of individual realizations , the ergodicity breaking and non - gaussianity parameters , as well as the probability density function @xmath1 we analyze the weakly non - ergodic character of the heterogeneous diffusion process and , particularly , the degree of irreproducibility of individual realization . as we show , the fluctuations between individual realizations increase with growing modulus @xmath2 of the scaling exponent . the fluctuations appear to diverge when the critical value @xmath3 is approached , while for even larger @xmath4 the fluctuations decrease , again . at criticality , the power - law behavior of the mean squared displacement changes to an exponentially fast growth , and the fluctuations of the time averaged mean squared displacement do not seem to converge for increasing number of realizations . from a systematic comparison we observe some striking similarities of the heterogeneous diffusion process with the familiar subdiffusive continuous time random walk process with power - law waiting time distribution and diverging characteristic waiting time .