Instruction: Write an abstract for the following acadamic paper.

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

Paper: in the last decade , studies about the structure and dynamics of complex networks have blossomed , thanks in particular to the versatility of the network representation , which has turned out to be adequate for systems as diverse as the internet or social networks . a large body of knowledge about the empirical description of networked systems has thus been accumulated , together with a wealth of modeling techniques ; a good level of understanding of how dynamical processes taking place on networks depend on their structure has been as well reached @xcite . many network studies have been concerned with systems of interest in several scientific areas a priori remote from physics ( social sciences , biology , computer science , epidemiology , ... ) , and they have also reached more traditional fields of statistical physics , such as the study of glassy systems , as we now describe . the many puzzles raised by the glass transition , and in particular the slow dynamics displayed by glassy systems at low temperatures , have been the subject of a large interest in the past decades @xcite . one of the approaches which has led to promising insights consists in the description of the dynamics of a glassy system inside its configuration space . the energy landscape of a glassy system is typically rugged , made of many local minima ( metastable states ) , whose huge number makes it difficult to reach equilibrium . in this framework , the energy landscape is seen as a set of basins of attractions of local minima ( `` traps '' ) , and the system evolves through a succession of harmonic vibrations inside traps and jumps between minima @xcite . this picture has stimulated the definition and study of various simplified models of dynamical evolution between traps , in order to reproduce the phenomenology of glassy dynamics @xcite . on the other hand , several studies have focused on obtaining a better understanding of the structure of these local minima . a way to attain this goal is to perform numerical simulations of small systems , at a fixed temperature , and quenching them at regular time intervals in order to make them reach the nearest local minimum . information is then gathered on the various local minima , and on the sizes of their basins of attraction . various studies have investigated , among other issues , the detailed structure of the potential energy landscape , the substructure of minima , and the properties of energy barriers between minima @xcite . several works have also used the information on the energy landscape to study a master equation for the time evolution of the probability to be in each minimum . the considered systems range from clusters of lennard - jones atoms to proteins or heteropolymers @xcite . an interesting property of the modeling of the energy landscape in terms of a set of traps linked by energy barriers , lies in the possibility to define and study its network representation within the context of network theory . in this representation , each local minimum is associated to a node , and a link is drawn between two nodes whenever it is possible for the system to jump between the basins of attraction of the corresponding minima . the links can then be defined as weighted and directed , as jumps between minima are not equiprobable , and may be easier in one direction than in the other . networks of local minima of the energy landscape have thus been built and studied . these networks have been found to exhibit a small - world character @xcite . the number of links of each node ( its degree ) turns out to be strongly heterogeneous , possibly with scale - free degree distributions , which have been linked to scale - free distributions of the areas of the basins of attraction @xcite . complex network analysis tools have also been used to investigate the structure of energy landscapes of various systems of interest , such as lennard - jones atoms , proteins , or spin glasses , among others @xcite . the energy of a minimum and its degree ( i.e. , the number of other minima which can be reached from this minimum ) have been shown to be correlated , as well as the barriers to overcome to escape from a minimum . in particular , a logarithmic dependence of the energy of a minimum on its degree has been exhibited , as well as energy barriers increasing as a ( small ) power of the degree of a node @xcite . no systematic study of these issues has however been performed , and most investigations have been limited to relatively small systems because of computational limitations . most importantly , the investigations cited above have focused on the topology of the network of minima , conceived as a tool to characterize the energy landscape . the structure of a network has however a deep impact on the properties of the dynamical processes which take place on it @xcite . it seems thus adequate to put to use the tools and techniques developed for the analysis of dynamical processes on networks to achieve a better understanding of how the energy landscape structure , represented as a network , affects the system performing a random walk in it , and how the onset of glassy dynamics can be described in this way in a general framework . in a previous paper @xcite , we have made a first step to fill this gap by focusing on the trap model put forward in ref . @xcite . in this paper , we generalize our approach to more involved transition rates between energy minima . we show how the heterogeneous mean - field ( hmf ) theory @xcite can be used in this context to highlight the connection between the topological properties of the network of minima and the dynamical exploration of these minima . we show in particular that the relationship between energy and degree of the minima is a crucial ingredient for the existence of a transition and the subsequent glassy phenomenology . our results shed light on the empirically found relationship between the energy of a local minimum and its degree , and we hope that they will stimulate more systematic investigations on this issue .    we have organized our paper as follows : in section  [ sec : model ] we define our model of energy landscape dynamics as a random walk on a complex network . different physical transition rates are proposed , and the corresponding numerical implementation is discussed . in section  [ sec : formalism ] we present a theoretical analysis based on the heterogeneous mean - field approximation for dynamical processes on complex networks . this formalism is applied in section [ sec : general - formalism ] , where general analytical approximate expression are presented for the main quantities characterizing the glassy transition and dynamics . these expressions are applied to the different physical transition rates considered in section  [ sec : appl - phys - trans ] , where checks against numerical simulations are also presented . in section  [ sec : barriers ] we discuss the relation between energy basins and energy barriers . finally , in section  [ sec : conclusions ] we present our conclusions . we consider a network of @xmath0 nodes , in which each vertex @xmath1 corresponds to a minimum in the energy landscape , and a link is drawn between two minima @xmath1 and @xmath2 if the system can jump directly from @xmath1 to @xmath2 . to each node @xmath1 is associated the energy @xmath3 of the corresponding minimum ( energies are defined from a reference level , in such a way that @xmath4 for all @xmath1 ) . moreover , an energy gap @xmath5 is associated to the edge between vertices @xmath1 and @xmath2 , as depicted in fig .  [ landscape ] : @xmath5 is a symmetric function , such that the energy barrier that must be overcome to jump from vertex @xmath1 to vertex @xmath2 can be written as @xmath6 and , analogously , @xmath7 . obviously , we will have in general @xmath8 . are measured positive downwards . energy gaps @xmath5 are defined positive upwards.,width=340 ]    the system under investigation is pictured as a walker exploring the network through a biased random walk . the rate ( probability per unit time ) @xmath9 to go from vertex @xmath1 to vertex @xmath2 depends a priori on the energy at vertices @xmath1 and @xmath2 and/or the energy barrier between @xmath1 and @xmath2 that must be overcome . the random walk model is defined in discrete time @xmath10 as follows :    * at time @xmath10 , the walker is in vertex @xmath1 . * it chooses at random a neighbor of @xmath1 , namely @xmath2 . * with a probability @xmath9 , that depends on the energy @xmath11 and/or on the energy barrier @xmath12 , the walker hops to vertex @xmath2 . * time is updated @xmath13 . the relationship between the probabilities @xmath9 and the energy and energy barriers can be of different forms . in usual unbiased random walks , @xmath9 is a constant independent of both @xmath1 and @xmath2 @xcite . as a first step to introduce a dependence on the nodes , a possible approach is that in which the energy barriers depend only on the local minima themselves , i.e. we consider @xmath14 . for example , in the bouchaud trap model considered in ref . @xcite , the probability to exit from a trap is just an arrhenius law depending only on the departing trap s depth , namely @xmath15 where @xmath16 is the inverse temperature and @xmath17 is a constant that determines a global timescale . other possible definitions include the metropolis one @xmath18 and the glauber rate @xmath19 we note that the rates considered in the bouchaud trap model are quite different from the metropolis or glauber rates . indeed , while the former depends only on the depth of the originating trap , the latter depend also on the energy of the arriving vertex . this translates in the fact that , in the limit of zero temperature , the dynamics of the bouchaud trap model is frozen for any @xmath11 , while metropolis and glauber dynamics still allow jumps to lower energy minima @xcite . within an even more realistic representation of glassy dynamics , one can also contemplate the case @xmath20 , allowing for the transition rates to depend explicitly on the energy barriers between adjacent minima . as a paradigm of this choice , we propose a rate of the arrhenius form @xmath21 which acts a straightforward generalization of the local transition rate ( [ eq : rate_traps ] ) . the case of rate ( [ eq : rate_traps ] ) ( local trapping ) was studied in a previous publication @xcite . in the following , we will consider in turn non - local rates ( [ eq : rate_metropolis ] ) , ( [ eq : rate_glauber ] ) and ( [ eq : rate_barriers ] ) and discuss the fundamental differences due to the introduction of energy barriers in the model . we emphasize that our model differs both from usual unbiased random walks , as the local energy determines the transition rates , and from mean - field trap models in which jumps between any pair of energy minima are a priori possible : here , the system can jump only between neighboring nodes . the dynamical evolution depends therefore both on the network topology and on the energies associated to the nodes .      to implement numerically the random walk , it is convenient to resort to the techniques developed for general diffusion processes on complex weighted networks @xcite . the main advantage of this method is to avoid rejection steps , thus improving dramatically the computational efficiency @xcite . therefore , at each simulation step the random walker sitting at node @xmath1 selects a neighbor @xmath2 with probability @xmath22 , where the sum in the normalizing factor is extended to all of @xmath1 s neighbors . as the walker hops on node @xmath2 , the physical time is incremented by an interval @xmath23 drawn from the exponential distribution @xmath24 , where @xmath25 is the inverse of the average escape rate out of node @xmath1 . in this way the simulation time is disentangled from the physical time and the latter has no impact on the simulation efficiency . no matter how much physical time a walker spends in a node , from the simulation time point of view it is always just a time step .    the network substrates on which we will focus are scale - free networks with a degree distribution of the form @xmath26 and @xmath27 . we will generate them with the uncorrelated configuration model ( ucm ) @xcite , that allows us to tune the degree distribution to the desired form and prevents the formation of degree - degree correlations . networks are therefore generated as follows : a number of stubs ( or semi - links ) extracted from the desired final degree distribution is assigned to each node . stubs are then randomly paired to form links between nodes , with the prescription that multiple links as well as self - loops must be avoided . a minimal degree @xmath28 is fixed a priori . to avoid spurious effects due to the possible presence of tree - like structures @xcite it is convenient to adopt @xmath29 . we will choose @xmath30 in all of our simulations . so far the algorithm coincides with the one of the configuration model @xcite , but the ucm introduces moreover a cutoff to the degree distribution , @xmath31 , which avoids the formation of degree correlations by limiting the size of the hubs @xcite . in order to gain analytical understanding of the role of the different transition rates in the corresponding glassy dynamics , we apply a standard heterogeneous mean - field ( hmf ) formalism @xcite . the basic tenet of hmf is the assumption that all the dynamical properties of a vertex depend only on its degree . vertices are thus grouped into classes according to their degree , and vertices with the same degree are treated as equivalent . this approximation is consistent with previous findings that have uncovered the correlations between the energy of a local minimum and the degree of the corresponding node in the network @xcite . we therefore make the assumption that there exists a relationship @xmath32 where the function @xmath33 is a characteristic of the model . this also means that the distributions of energies @xmath34 of the system s landscape , and the degree distribution @xmath35 of the corresponding network are linked through @xmath36 . in the same spirit , we make the further assumption that the energy gap between minima @xmath1 and @xmath2 depends only on the degrees of @xmath1 and @xmath2 , i.e. , that it can be written as @xmath37 , where @xmath38 is a symmetric function of @xmath39 and @xmath40 .    under the hmf approximation the dynamics will thus focus on the transitions between different degree classes . the rate to go from a vertex @xmath39 to a vertex @xmath40 can be written as @xmath41 the function @xmath42 , defined as the conditional probability that a vertex of degree @xmath39 is connected with another vertex of degree @xmath40 @xcite , takes into account the topological features of the network , by gauging the probability of selecting a vertex @xmath40 as neighbor of @xmath39 . the function @xmath43 measures the rate of jumping from a vertex of degree @xmath39 to a vertex of degree @xmath40 ( given that they are connected by an edge ) , and depends on @xmath39 and @xmath40 through the rates @xmath44 and the functions @xmath36 and @xmath45 . obviously , the rate @xmath43 is not in general a symmetric function of @xmath39 and @xmath40 . it is worth noting that , apart from a normalization , equation  ( [ eq : weighted_propagator ] ) is simply the so - called weighted propagator describing the probability that a node in class @xmath39 interacts with a node in class @xmath40 @xcite . we also note that the rates @xmath43 depend on the inverse temperature @xmath46 through the microscopic rates @xmath9 . in this section , we apply the hmf theory to compute different quantities relevant for the characterization of the dynamics of a random walk in a complex energy landscape represented in terms of a network of minima . the description of a random walk dynamics starts from the occupation probability @xmath47 , defined as the probability for the walker to be in _ any _ node of degree @xmath39 at a time @xmath48 . its time evolution can be easily represented in terms of a master equation of the form @xmath49 upon describing the state at time @xmath48 with the row vector @xmath50 , where @xmath51 is the cutoff or largest degree in the network , equation ( [ rate_scalar ] ) can be rewritten in vector form as @xmath52 where the matrix @xmath53 , with elements @xmath54 is a generalization of a laplacian matrix to the case of directed weighted graphs . the matrix elements satisfy @xmath55 which ensures conservation of probability and states that the columns of @xmath53 are not linearly independent . the real part of every eigenvalue of @xmath53 is non - negative @xcite . as a consequence , all solutions of eq . ( [ rate_vectorial ] ) , which can be formally written as @xmath56 are stable according to lyapunov criteria . in particular , since @xmath57 , @xmath53 always has the eigenvalue @xmath58 , which corresponds to a constant solution of the problem . at this point one can proceed in close analogy with discrete - time regular markov chains @xcite . by making the assumption that the matrix @xmath59 is non - negative and irreducible ( indeed it is for every choice of @xmath60 in the following ) , we can prove that the @xmath58 eigenvalue of @xmath53 has algebraic multiplicity @xmath61 . hence , the stationary solution of eq . ( [ rate_vectorial ] ) is unique .      in order to calculate the steady solution @xmath62 in the limit @xmath63 , one can impose @xmath64 . this leads to the condition @xmath65 so that we are left with the task of finding the left nullspace of @xmath53 . ( [ homogeneous_v ] ) is a homogeneous system of algebraic linear equations . it admits non - trivial solutions since @xmath66 . in our case , the solution to ( [ homogeneous_v ] ) can be easily found by imposing the detailed balance condition . namely , writing eq . ( [ homogeneous_v ] ) as @xmath67=0 ,      \label{eq:1}\ ] ] we can obtain a solution by imposing that the terms inside the summation in eq . ( [ eq:1 ] ) cancel individually , that is @xmath68 substituting the form of @xmath69 , we obtain @xmath70 where in the last step we have used the degree detailed balance condition @xmath71 which simply expresses that the number of edges from a node of degree @xmath39 to a node of degree @xmath40 is equal to the number of edges from a node of degree @xmath40 to a node of degree @xmath39 @xcite . from eq . ( [ eq:2 ] ) , we see that its right - hand - side must be expressible as a simple ratio of a function of @xmath39 over a function of @xmath40 . a general way to obtain this is to impose a coarse - grained rate @xmath60 taking the general form @xmath72 in other words , we assume that the rate @xmath60 can be written as the product of a function of @xmath40 , a function of @xmath39 , and a symmetric function @xmath73 ( where @xmath39 and @xmath40 need not be separable ) . we will see later that all the rates @xmath74 defined in sec . [ sec : definition ] ( traps , glauber , metropolis , and energy barriers ) can be written in such a form . the stationary solution is then given by @xmath75 where @xmath76 is a normalizing constant determined by the condition @xmath77 . such a solution is unique , as proven above . interestingly , the symmetric function @xmath78 does not enter the steady solution , although it will play a role in affecting the transient behavior , as we will see in the next sections .      the steady state solution found above for the occupation probability is defined if and only if the normalization constant @xmath79 is finite . when this condition is met , the random walker reaches an equilibrium state with a distribution @xmath80 . on the other hand , whenever such condition is not met , the random walker is unable to reach a steady state , i.e. the steady solution to the rate equation does not correspond to any physical steady state in equilibrium @xmath81 . we identify this region of the phase space with the glass phase for our random walker @xcite .    the functions @xmath82 and @xmath83 depend on the temperature , on the precise dynamics chosen ( traps , glauber , metropolis ) , and encode the relationship @xmath36 between energy and degree of the minima . the degree distribution moreover enters explicitly the expression @xmath84 . as the various parameters of the model are changed , it is thus a priori possible to go from one phase in which @xmath84 is finite to one in which @xmath84 diverges . in a physical system in particular , the control parameter is usually the temperature , while the topology of the network of minima and the function @xmath36 are given . it is then clear from eq . ( [ eq:4 ] ) that the presence or absence of a finite glass transition temperature @xmath85 , such that @xmath84 becomes infinite for @xmath86 , depends on the interplay between the topology of the landscape network ( as determined by @xmath35 ) and the functions @xmath82 and @xmath83 . interestingly , at this mean - field level , the existence of a transition does not depend on the network degree correlations , since the conditional probability @xmath42 do not enter eq . ( [ eq:4 ] ) . let us consider for instance a network of minima with a heavy - tailed degree distribution such as @xmath26 . a transition between a finite and an infinite @xmath84 can be observed if and only if @xmath87 behaves at large @xmath39 in the form @xmath88 where the exponent @xmath89 depends on the temperature , and can take values smaller or larger than @xmath90 depending on the temperature . another example is given by a stretched exponential form for @xmath35 , @xmath91 , in which case a transition is observed if and only if @xmath87 is of the form @xmath92 , with @xmath93 a function of the temperature ( the transition is then given by @xmath94 ) .      in any finite system , unless the product function @xmath87 exhibits some sort of singularity , the normalization constant eq . ( [ eq:4 ] ) is finite and the steady state distribution @xmath95 exists , the occupation probability @xmath47 converging to it after an equilibration time , i.e. @xmath96 the corresponding thermalization of the occupation probability occurs in a way depending on the function @xmath36 . shallow energy minima are indeed explored first , while deep traps ( large @xmath97 ) are visited at larger times @xcite . if @xmath36 is a growing function of @xmath39 , as indeed found empirically @xcite , small degree nodes correspond to shallow minima , and deeper minima are associated to larger nodes . the evolution of @xmath47 takes then place in a hierarchical fashion : the small degree region equilibrates first , and progressive equilibration of larger degree regions takes place at larger times . in this respect , we obtain a strong difference between the biased random walk that the glassy system experiences and usual diffusion processes corresponding to unbiased random walks , which visit first large degree vertices and then cascade down towards small degree nodes @xcite , in the present case we observe an `` inverse cascade process '' from small vertices to hubs . we have found in @xcite that , in the case of a random walk among traps , this hierarchical thermalization is summarized in a scaling form for @xmath47 , which can be written as @xmath98 where @xmath99 represents the maximum degree of the vertices equilibrated up to time @xmath48 , and @xmath100 interpolates between @xmath101 at small @xmath102 and the short time form of @xmath103 which is proportional to @xmath104 . we will see in the next section that a similar scaling is obeyed for other transition rates . in general , for the glassy dynamics , the functional form of @xmath99 can moreover be obtained through the following argument : the total time @xmath48 can be written as the sum of the trapping times spent in the vertices that have been visited since the beginning of the dynamics . trapping times increase with the depth of the minimum , hence with the degree ( we are still considering the case of an increasing function @xmath33 ) , and , in the glassy phase , the consequence is that the sum of trapping times is dominated by the vertex with the largest degree visited up to that point , namely @xmath105 . moreover , the average trapping time @xmath106 at a given vertex @xmath39 can be estimated as the inverse of the average escape rate @xmath107 from that vertex : @xmath108 we can therefore estimate @xmath99 , the typical degree up to which nodes are `` equilibrated '' at time @xmath48 , by approximating @xmath109 , and solving the equation @xmath110 to obtain @xmath105 as a function of @xmath48 . note that the result depends here on the function @xmath111 and not only on @xmath82 , @xmath83 , @xmath36 and @xmath35 . the properties of the system can be further quantified by measuring the average time @xmath112 required by the random walker to escape from the vertex it occupies at time @xmath48 @xcite . for small waiting times @xmath48 , @xmath113 increases as a result of the transient equilibration of @xmath114 . for large @xmath48 , such that @xmath115 is close enough to the equilibrium @xmath116 , @xmath113 can be calculated instead as the average @xmath117 \tau_k\ ] ] where @xmath118 is the inverse of the equilibrium escape rate , cf eq . ( [ eq:7 ] ) , yielding @xmath119 most interestingly , the explicit form of the average escape time @xmath113 depends explicitely on the symmetric function @xmath111 as well as on the network degree correlations , as expressed by the conditional probability @xmath42 . let us go back to the issue of the existence of a glass transition in the model . we first recall the phenomenology of the fully connected trap model , with transition rates @xmath120 for any @xmath1 and @xmath2 , where the energies @xmath11 are random numbers extracted from a distribution @xmath34 @xcite . as all traps are connected with each other , all traps are equiprobable after a jump , so that the probability for the system to be in a trap of depth @xmath97 is simply @xmath34 , and the average rest time spent in a trap is @xmath121 . a transition between a high temperature phase and a glassy one is thus obtained if and only if , when @xmath46 increases , @xmath122 is finite at small @xmath46 and diverges at a finite @xmath85 . such a phenomenology is obtained if and only if @xmath34 is of the form @xmath123 at large @xmath97 ( else the transition temperature is either @xmath58 or @xmath124 ) , and the transition temperature is then @xmath125 @xcite .    in the present case of a nework of minima , the average rest time that the walker spends in a minimum is @xmath126 where the symbol @xmath127 refers to the average performed over the measure @xmath128 , which represents the probability that the walker is in any vertex of degree @xmath39 after a hop . note that we disregard here the _ physical _ time , which is the sum of times spent in the various minima , but consider only the number of hops between minima . in the case of the traps model , @xmath129 is simply given by the probability to be in a node of degree @xmath39 after a hop in a random walk , i.e. by @xmath130 @xcite , since the transition rates do not depend on the arrival node . in a non - local trapping model instead , we need to write a master equation of the form @xmath131 where the matrix @xmath132 is now stochastic and the derivative is intended with respect to the number of hops . in the long time limit , we impose @xmath133 and calculate @xmath128 as we did for @xmath95 , imposing the detailed balance condition , and obtaining @xmath134   r_k \ , \ ] ] where @xmath135 is a normalization factor , given by @xmath136    we finally obtain for the average @xmath137 : @xmath138 where @xmath84 is the normalization factor of @xmath139 defined in eq . ( [ eq:4 ] ) . as for the average escape time , the average rest time @xmath137 thus depends on all the parameters of the system , including the network s degree correlations and the symmetric function @xmath140 . in this section , we apply the general hmf results obtained in the previous section to physical transition probabilities between local minima given by the trap model , glauber , metropolis and barrier - mediated rates . we will focus for definiteness on scale - free networks characterized by a power - law degree distribution @xmath141 with @xmath142 , which turns out to be the interval reported in the literature @xcite . let us first consider the explicit form of the transition rates in each case , to show that they can be cast in the form given by eq . ( [ eq:3 ] ) . in the case of the trap model , the rate to jump from a vertex @xmath39 to a vertex @xmath40 is simply @xmath143 , where we recall that @xmath33 gives the depth of a node of degree @xmath39 : it depends only on the degree of the starting node , and not on the node reached after the jump . we can therefore use @xmath144 the glauber rate can be written as @xmath145 leading to @xmath146 the metropolis transition , on its turn , reads @xmath147.\ ] ] since , for positive @xmath89 , @xmath148 , we can choose @xmath149 finally , in the presence of energy barriers , the transition rate reads @xmath150 where @xmath151 is a symmetric function of its arguments , so that we can use @xmath152      interestingly , for all the transition rates considered above , the product of the functions @xmath153 and @xmath83 , which controls the existence of the steady state solution of the occupation probability , takes the form @xmath154 the normalization constant @xmath84 defined in eq . ( [ eq:4 ] ) can thus be written as @xmath155 for a power - law degree distribution @xmath26 , a finite glass transition temperature is then obtained if and only if @xmath36 is of the form @xmath156 which is precisely what has been found , in conjunction with a scale - free degree distribution , in ref . @xmath84 is then indeed given by a sum of terms of the form @xmath157 , which converges if and only if @xmath158 in other words , a transition between a high temperature phase in which @xmath95 exists and a low temperature glassy phase is obtained at the critical temperature  @xcite @xmath159    quite noticeably , the existence of a transition at a finite temperature , as well as the value of this temperature , does not depend on the form of the transition rates between the local minima , but only on the existence of a particular interplay between the topology of the network of minima and the relationship between energy and degree in this network , as determined by the function @xmath36 . we emphasize that this result is also independent of the network degree correlations @xmath42 , as already noted in the previous section . let us focus on the case of a scale - free network of minima , with @xmath160 and @xmath161 . for any of the rates discussed above , the steady state measure , when it exists , is given by @xmath162 where @xmath163 is the riemann @xmath163 function . a plot of @xmath164 as a function of @xmath39 and @xmath165 is given in fig . [ q_k ] , while data from simulations are reported in fig . [ q_k_sim ] for the evolution of @xmath115 under glauber dynamics . above the transition ( @xmath166 ) , low-@xmath39 states ( i.e. , shallow minima ) are more probable . as the temperature decreases , @xmath95 becomes less and less peaked at low values of @xmath39 , and large-@xmath39 states , which correspond to lower energies , become more and more probable . for the random walker to be in any node of degree @xmath39 . for @xmath167 the system undergoes a transition to a glassy state.,width=340 ]    in any finite system , the sum defining @xmath84 , is finite at any temperature as the degree distribution has a cut - off at a finite @xmath51 : @xmath168 for instance , for @xmath169 , and with @xmath170 , the sum @xmath171 is analytic in @xmath172 for any finite @xmath51 . here @xmath173 is the harmonic number of order @xmath174 , which tends to @xmath175 for @xmath176 . the probability @xmath95 is thus well defined for every @xmath165 and for any finite system . in particular , performing a continuous degree approximation in eq . ( [ eq:19 ] ) , we can obtain an estimate of the network size dependence of @xmath84 as @xmath177 for @xmath178 , @xmath84 tends to a constant as the network size ( and thus @xmath51 ) increases . on the other hand , for @xmath179 , @xmath84 diverges as @xmath180 , i.e. , as @xmath181 in uncorrelated scale - free networks , which obey @xmath182 .     for glauber dynamics . the distribution measured after a small waiting time @xmath48 is determined by a usual unbiased random walk behavior , i.e @xmath183 , while at larger times it relaxes to the equilibrium @xmath95 . the relaxation towards equilibrium starts from small degree nodes . data refer to ucm networks with @xmath184 and @xmath185 for @xmath186 ( @xmath187 ) : with these parameter values , for small times @xmath188 , while at large enough times @xmath189.,width=340 ]      at low temperatures , even for a finite system , the evolution of @xmath190 towards @xmath95 is slow , as displayed in fig . [ q_k_sim ] , and an ageing regime takes place , in which the function @xmath115 obeys the scaling form @xmath191 where the characteristic degree @xmath105 can be estimated from eq . ( [ eq:8 ] ) . in order to simplify its computation , we will consider an uncorrelated network of minima , such that @xmath192 , and we will work in the continuous degree approximation , using the normalized form @xmath193 , where @xmath28 is the minimum degree present in the network .    in the case of the glauber dynamics , the escape rate can be expressed , within the above approximations , as @xmath194 ,   \label{eq:21}\end{aligned}\ ] ] where we have used the relation @xmath195 and where @xmath196 $ ] is the gauss hypergeometric function . using the asymptotic expansion for @xmath197 @xcite , we obtain that the leading behavior for large @xmath39 yields @xmath198 which leads to @xmath199 from here , using the relation @xmath109 , we obtain @xmath200     at different temperatures ( glauber dynamics ) . data refer to ucm networks with @xmath184 and @xmath201 , so that @xmath202 ( where we have taken @xmath187 ) . the top panel presents data for @xmath203 while both middle and bottom panels concern the low temperature case @xmath204 . accordingly , for the top panel we use @xmath205 for the rescaling , while for the central and the bottom ones it holds @xmath206 ( see eq . ( [ eq:24 ] ) ) . the curves corresponding to different @xmath48 collapse well under this rescaling . note that we use rather small values of @xmath48 , because the equilibration time , defined by @xmath207 , is @xmath208 for @xmath204 and @xmath209 for @xmath210 . each curve is obtained by averaging over @xmath211 simulation runs.,width=340 ]    in fig . [ collapsep ] we check the validity of eqs . ( [ eq:25 ] ) and ( [ eq:24 ] ) by performing a data collapse analysis for different values of @xmath48 . the curves obtained for different @xmath48 collapse indeed as predicted .    in the case of the metropolis transition rates , a similar analysis yields @xmath212 leading to the same asymptotic behavior as in eq . ( [ eq:22 ] ) , and therefore to the same scaling picture as for the glauber rate .    as pointed out for the case of the traps model @xcite , however , in finite systems the scaling relations in eq . ( [ eq:24 ] ) hold only as far as @xmath213 , i.e. it exists an equilibration time @xmath214 , obtained by inverting eq . ( [ eq:24 ] ) , above which the system has completely relaxed and eq . ( [ eq:25 ] ) is no more valid . finally , it is worth stressing that , while for large temperatures the scaling exponent relating @xmath105 and @xmath48 depends on the temperature , in the low temperature phase it becomes independent , being just proportional to the transition temperature . we note that this saturation of the exponent at @xmath215 is very different from the phenomenology obtained in the trap model @xcite , for which @xmath216 . an immediate consequence is that the equilibration time strongly depends on @xmath46 , as @xmath217 , for a system described by traps , but is given by @xmath218 for any @xmath204 for glauber and metropolis rates . contrarily to results for the glass transition temperature @xmath219 and the steady state , the glassy dynamics for barrier - mediated rates does not yield the same results as for glauber and metropolis rates , since @xmath107 does depend on the symmetric function @xmath220 . in particular , we need here to choose a functional form for @xmath45 . we propose to use @xmath221 which will be justified in section [ sec : barriers ] . in this case we obtain @xmath222 in order to keep an interesting phenomenology , the constant @xmath223 can not be chosen arbitrarily . if @xmath223 were independent of the system size , the escape rate would be dominated by the exponential at all temperatures . this behavior would reflect the fact that in this case the rate is suppressed in transitions involving nodes of large @xmath39 , eventually generating unphysically large barriers at large @xmath51 . to prevent the system from building up infinite barriers , one can impose @xmath224 such that the maximum barrier is always comparable to the lowest energy minimum and neither term dominates the other . as a consequence , we take @xmath225 where @xmath226 is now constant and size independent . contrarily to the previous cases , @xmath99 is hard to determine , as no explicit inversion of eqs . ( [ eq:8 ] ) and ( [ eq:29 ] ) can be provided for the range of parameters of interest in our study . a numerical evaluation of @xmath99 , is reported in figure [ f : kw_barriers ] . the maximum degree of equilibrated vertices @xmath105 has an initial power law increase in time , which is reminiscent of the local trapping model . however , as larger degree nodes are equilibrated , exponential barriers come into play and the hierarchical thermalization becomes logarithmic in time .     for barrier - mediated dynamics . curves are obtained from the numerical inversion of eq . ( [ eq:8 ] ) . values of parameters are : @xmath201 , @xmath227 , @xmath228 , @xmath229 , @xmath230 . , width=340 ]      , so that @xmath231 ( @xmath187 ) . the scaling forms of eq . ( [ eq : escape ] ) produce a collapse of the curves concerning different systems sizes in the three regimes of high , intermediate and low temperature ( top , center and bottom panels , respectively ) . the slightly worse collapse obtained for @xmath232 may be due to logarithmic corrections as @xmath46 is close to both @xmath85 and @xmath233 . each point is averaged over @xmath234 simulation runs ( @xmath235 runs on each of @xmath235 network realizations ) . , width=340 ]    the average escape time , defined as the average time required by the system to escape from the vertex it occupies , can be computed in the long time limit from eq . ( [ eq : t_esc ] ) , as a function of the average trapping time @xmath106 in vertices of degree @xmath39 . from the asymptotic expansions of @xmath106 in eq . ( [ eq:23 ] ) , valid for glauber and metropolis dynamics , evaluation of eq . ( [ eq : t_esc ] ) allows us to observe that , whenever a finite @xmath219 exists , @xmath236 diverges at @xmath237 in an infinite system , as was already observed in the case of local trapping @xcite . it is noticeable that the same divergence temperature is obtained , as eq . ( [ eq : t_esc ] ) a priori involves the network s degree correlations and the function @xmath140 . within the continuous degree approximation , the divergence of the escape time with the system size follows the scaling laws @xmath238 as noted in the previous paragraph for the equilibration time , we note that the scaling for @xmath239 differs from the form @xmath240 encountered for local trapping @xcite . figure [ f : escape ] reports simulation data that confirm the validity of eq . ( [ eq : escape ] ) . as the temperature is lowered , the initial transient becomes longer , but for large enough times @xmath48 the asymptotic behavior predicted in eq . ( [ eq : escape ] ) is reached , as made clear from the collapse of curves concerning different system sizes .    in the case of barrier - mediated dynamics , @xmath106 depends on the symmetric function @xmath151 , as expressed in eq . ( [ eq:29 ] ) , namely @xmath241 proceeding as above we obtain @xmath242e_0 } & \beta_c/(2+\epsilon ) <        \beta < \beta_c \\ \mbox{const . }   & \beta<\beta_c/(2+\epsilon )        \\      \end{array }    \right .   , \end{aligned}\ ] ] where we recall that @xmath226 does not depend on the system size .      , as a function of the inverse temperature , for different system sizes . data are obtained by numerical computation of eq . ( [ tau_zi ] ) ( with @xmath187 ) . note the exponential increase with @xmath46 for @xmath243 , which saturates for @xmath244 as predicted by eq . ( [ e : rest]).,width=340 ]    , as a function of the degree distribution cut - off and for different temperatures . data are obtained by numerical computation of eq . ( [ tau_zi ] ) . @xmath245 grows as a power - law of @xmath51 , with an exponent that grows as @xmath46 increases ( going from bottom to top in the figure ) . the thick - gray line corresponds to @xmath246 , while the thick - black line corresponds to @xmath247 . for larger values of @xmath46 , the power - law behavior corresponds to the predicted @xmath248 , which no more depends on @xmath46 . the @xmath249 curve is reported as a dashed line for reference . values of @xmath46 represented here are comprised between @xmath250 and @xmath251.,width=340 ]    the hmf expression for the asymptotic average rest time , defined as the average time spent by the system in a minimum , is given by eq . ( [ tau_zi ] ) , namely @xmath252 , where the quantities @xmath84 and @xmath135 , for uncorrelated scale - free networks and a degree - energy relation @xmath253 , take the form , in the continuous degree approximation , @xmath254    let us first recall the case of the local trap model . both @xmath255 and @xmath111 are then constants , so that @xmath256 . thus , the average rest time behaves as @xmath84 : it is finite for @xmath210 , and diverges with the system size as @xmath257 for @xmath204 , signaling the emergence of the glassy regime at low temperatures .    in the case of glauber and metropolis dynamics ( both leading to the same results ) , the situation is more involved , since the product @xmath258 entering @xmath135 is not constant . in fact , @xmath135 diverges with @xmath51 for @xmath259 , that is , at a lower temperature given by @xmath260 . the interplay of these two temperatures determines the behavior of the system for finite sizes within the low temperature phase . in particular , for the glauber dynamics with @xmath261 and @xmath262 $ ] , we have @xmath263 upon considering lower values of @xmath46 , @xmath137 first encounters the divergence of @xmath84 at @xmath85 , which is then partially regularized by the divergence of @xmath135 at @xmath233 . from these results , we obtain the emergence of three scaling regimes for the behavior of @xmath137 as a function of the system size : @xmath264     for different transition rates , for a random walker on ucm networks with @xmath265 ( @xmath187 ) and @xmath184 . discrete points : simulation results ; continuous lines : theoretical predictions from @xmath266 , based on simulation parameters . the glauber and metropolis transition rates induce two changes of behavior at @xmath85 and at @xmath233 , the first being a steep increase of the average rest time @xmath137 and the second a smoothing / saturation of this increase . no saturation of @xmath137 is instead observed when barriers are present . the agreement with theoretical predictions is remarkable , thus corroborating the validity of the hmf assumptions . moderate deviations are found only in the case of barriers , where exponential growth is expected to add greater fluctuations . in the inset , the difference between the two behaviors is more evident thanks to a different scale of the plot . note that , since in the simulations @xmath187 , the high temperature limit of the rest time is different for the glauber and metropolis dynamics , being @xmath267 and @xmath268 respectively . for the case of barriers we have chosen @xmath269 . each point is obtained by averaging the rest times corresponding to the first @xmath270 hops of the random walker in each of @xmath271 network realizations . , width=340 ]    the direct numerical computation of eq . ( [ tau_zi ] ) is shown in figs . [ 2temperatures ] and [ kcscaling ] , showing the validity of this analysis . in particular , the exponential increase of @xmath272 with @xmath46 in the intermediate temperature range @xmath273 is clearly apparent in fig . [ 2temperatures ] , and fig . [ kcscaling ] confirms that , for @xmath274 , the exponent in the scaling law for the system size @xmath51 does not depend on the temperature . while the temperature @xmath85 signals the onset of the low temperature phase with glassy dynamics for all considered transition rates , for glauber / metropolis dynamics the low temperature phase can be further divided into two regions that correspond to different behaviors of the timescales with the system size . ( @xmath187 , @xmath231 ) . data for different system sizes collapse well when rescaled according to the theoretical values of eq . ( [ e : rest ] ) . while the agreement is excellent both for high and low temperature ( top and bottom panels , respectively ) , logarithmic corrections are probably present for the regime of intermediate temperatures @xmath275 ( central panel ) . in each simulation run the rest interval starting before @xmath48 and ending after @xmath48 is considered , and each point in the figure is averaged over @xmath234 simulation runs ( @xmath235 runs on each of @xmath235 network instances ) . , width=340 ]    figures [ simul_2temperatures ] and [ f : rest ] moreover show the result of numerical simulations of random walkers on scale - free networks for glauber and metropolis dynamics as well as in the case of barriers , globally confirming the above discussed picture . dynamics in the presence of barriers do not yield the same phenomenology as glauber and metropolis rate . in this case , we have @xmath276 and @xmath277 . selecting @xmath278 , as in section  [ sec : glassy - dynamics-1 ] , we are led to @xmath279 ^ 2.\ ] ] as for the escape time @xmath113 , upon choosing @xmath45 size independent , the rest time @xmath122 will be diverging exponentially with @xmath51 at every temperature . by introducing the size dependence as in eq . ( [ eq : size_depen ] ) , instead , one can see that the @xmath135 integral converges to a constant for large @xmath51 so that one is left with @xmath280 we therefore obtain the same picture as in the case of traps , with an exponential increase of @xmath122 as @xmath46 increases , as confirmed by numerical simulations in fig . [ simul_2temperatures ] . inspired by analogies with systems governed by the arrhenius law , we have introduced a transition rate that takes into account the energy barriers between states . within the heterogeneous mean - field approximation , in which all variables depend only on the degree of the vertices , and choosing @xmath281 , the transition rate we have considered becomes @xmath282 where @xmath151 is a symmetric function of the degrees of the two nodes . this model represents in essence an extension of the local trap model , where non - locality enters only in the form of symmetric energy gaps @xmath5 ( see fig . [ landscape ] ) . the steady state has exactly the same form as the ones discussed so far , which incidentally is the same as for the local trap model . as shown in the previous sections , the presence of barriers affect transient relaxation phenomena , but not the steady state . a different question is whether one can be more specific about the realistic functional form of the coarse - grained function @xmath151 . in the previous section we have already introduced a definition of @xmath151 . here we provide the rationale behind that choice . numerical simulations of the energy - landscape network of lennard - jones clusters show that the average barrier to escape from state @xmath39 follow the power law @xmath283 , with @xmath284 @xcite . in our model , such average can be computed as @xmath285.\ ] ] for simplicity we focus on uncorrelated networks , as simulations indeed show weak degree correlations . under this assumption , the first term of the sum on the right - hand side of eq . ( [ avg_barrier ] ) will contribute as a logarithm of @xmath39 and the power - law behavior of @xmath286 is possible whenever @xmath287 , which leads to consider the form proposed in previous sections , @xmath288 where @xmath289 has the dimensions of an energy ( a discussion about the possible values of @xmath223 is given in section [ sec : glassy - dynamics-1 ] ) . more complicated functional forms can also be proposed , for example accounting for barriers of different signs , as long as they retain the same power - law behavior of eq . ( [ sigma_linear ] ) in the large @xmath39 limit . it is interesting to notice that @xmath283 implies that the average escape rate @xmath290 has the form of a stretched exponential @xmath291 , if we neglect the logarithmic correction . in this paper , we have presented a simple mathematical framework for the description of the dynamics of glassy systems in terms of a random walk in a complex energy landscape . we have shown how to incorporate into this picture the network representation of this landscape , put forward and studied by several authors @xcite , in order to go beyond simple mean - field models of random walks between traps that are all connected to each other . while our previous work had focused on the case of a landscape consisting of traps connected by a network @xcite , we have here generalized our study to more involved and realistic transition rates between minima , including glauber or metropolis rates , and the possibility of energy barriers between minima . we have shown how the interplay between the topology of the network of minima and the relationship between the energy and the degree of a minimum may determine a rich phenomenology , with the existence of two phases and of glassy dynamics at low temperature . interestingly , the existence of these phases , and the transition temperature , do not depend on the network s degree correlations nor on the precise form of the transition rates , but other more detailed properties do . in the case of glauber and metropols dynamics , the low temperature phase can be further divided into two regions with different scaling properties of the average trapping time as a function of the temperature . overall , our results rationalize and link the empirical findings about correlations between the energy of the minima and their degree , and should stimulate further investigations on this issue . our work has also interesting applications in terms of diffusion phenomena on complex networks , and shows that non trivial transition rates can lead to a very interesting phenomenology . usual random walks lead to a higher probability for the random walker to be in a large degree node ( @xmath292 ) , with respect to the random choice of a node ( @xmath293 ) ; here , the models we have studied can lead to various stationary probabilities , such as for instance a uniform coverage which does not depend anymore on the degree . interestingly , the biased random walks among traps that we have studied can even display a phase transition phenomenon , as either a temperature parameter or the network s properties are changed , with the possible presence of a glassy phase with slow dynamics . p. m. , r.p .- s . , and a. baronchelli acknowledge financial support from the spanish mec , under project fis2010 - 21781-c02 - 01 , and the junta de andaluca , under project no . p09-fqm4682 .. r.p .- s . acknowledges additional support through icrea academia , funded by the generalitat de catalunya . a. baronchelli acknowledges support of spanish mci through the juan de la cierva program funded by the european social fund .
Abstract:

Answer: we present a simple mathematical framework for the description of the dynamics of glassy systems in terms of a random walk in a complex energy landscape pictured as a network of minima . we show how to use the tools developed for the study of dynamical processes on complex networks , in order to go beyond mean - field models that consider that all minima are connected to each other . we consider several possibilities for the transition rates between minima , and show that in all cases the existence of a glassy phase depends on a delicate interplay between the network s topology and the relationship between energy and degree of a minimum . interestingly , the network s degree correlations and the details of the transition rates do not play any role in the existence ( nor in the value ) of the transition temperature , but have an impact only on more involved properties . for glauber or metropolis rates in particular , we find that the low - temperature phase can be further divided into two regions with different scaling properties of the average trapping time . overall , our results rationalize and link the empirical findings about correlations between the energy of the minima and their degree , and should stimulate further investigations on this issue .