Instruction: Write an abstract for the following acadamic paper.

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

Paper: the emergence of slow modes of relaxations is a major signature of the glassy behavior . @xcite the first consequence of large equilibration time is that the properties of a glass evolve on experimental timescales , the system ages . despite the huge amount of work that has been devoted to glassy systems , a general understanding of slow relaxation processes and of their relation to macroscopic rheology is still lacking . for example , experimental and numerical data converge to a stretched exponential shape of relaxation phenomena , the so - called kohlrausch - williams - watts ( kww ) relaxation ,  @xcite while mode coupling theory predicts power law behavior . @xcite when it comes to extend these results to non - linear rheology , the situation is even more difficult from the theoretical point of view , since the system is driven far from equilibrium . at an applied strain rate , power law viscosities or the emergence of a yield stress , are common features of complex fluids . @xcite md simulations of lennard - jones systems  @xcite and recent experiments on colloidal glasses  @xcite are also consistent with such properties . recent theoretical approaches to non - linear rheology  @xcite have called on ideas originating in mode coupling theory , or trap models of spin glasses . @xcite these approaches , however , fail to provide a description of non - linear rheology that is consistent with experimental observations ;  @xcite moreover , these models , based only on a phase - space picture of glassiness , do not lead to a clear identification of real - space physical mechanism that are responsible for an observed macroscopic behavior . some other attempts have been directed towards phenomenological approaches ,  @xcite but those are still unsatisfactory , and do not offer a microscopic picture , even heuristic , of the undergoing physical mechanisms . the current work lies at the convergence point of several recent theoretical ideas in various fields . firstly , it is based on the idea by liu , nagel and coworkers , that some unique mechanism lies behind the profound similarities displayed by structural glassy materials , be they colloidal glasses , granular materials , foams ,  @xcite drawing on this idea , the fundemental assumption underlying the current work is that macroscopic properties of dense materials results from major _ dynamical _ properties of structural rearrangement of the contact network . another important element of the current work comes from the study of density relaxation in granular materials . @xcite under vertical tapping , granular materials present a time - logarithmic relaxation of the occupied volume ; this is in strong analogy with free - volume relaxation , commonly observed in glasses , and of central importance for aging . @xcite however , in the case of tapped granular materials , density relaxation can not result from thermal activation since there is no thermal bath : the logarithmic density relaxation has been explained by introducing an equation of motion for free - volume , which characterizes the internal state of the granular material .    finally , this work uses a theory introduced recently by falk and langer to account for elasto - plastic transition in amorphous solids . @xcite the so - called shear transformation zone ( stz ) theory , originates from ideas by argon , spaepen and others to describe creep in metallic alloys ;  @xcite it has been primarily developed in the framework of fracture mechanics , and provides a general scheme for jamming as resulting from a structuration of the material by creation of arrangements oriented along principal directions of an applied stress . stz theory calls to free - volume arguments , but in fact , free - volume remains a parameter and enters constants of the theory . these ideas lead me to propose a general approach to the rheology of dense materials , leading to a simple set of constitutive equations . @xcite the purpose of the this article is to elaborate the theoretical arguments underlying these equations a to detail some of their consequences . two types of internal state variables are used : free - volume , which is an isotropic property of a molecular packing , and populations of arrangements , which are related to the distortion of the contact network between molecules ( hence introduce anisotropy ) . those variables evolves dynamically by free - volume activated processes : two types of rearrangements are associated to compaction and shear motion . free - volume enters activation factors that control both types of elementary rearrangements , and it is also a dynamical quantity . kww relaxation has been shown to result from the existence of a distribution of timescales ;  @xcite and it seems that this is often thought to be a necessary requirement . the first contribution of the current work is to show that it is not necessary , and that kww relaxation also result from a more simple mechanism : the interplay between free - volume logarithmic relaxation and free - volume activated shear deformation . the second contribution of this work is to establish a relation between kww relaxation , time - logarithmic density relaxation , and power law viscosity of a strongly sheared glassy material . a single parameter of the theory determines these three types of phenomena . power law viscosity will be shown to result from the existence of shear induced dilatancy . i propose one of the simplest set of equations that can account for those three properties is a single framework . the third contribution of this work is to present a theoretical framework which accounts for two coexisting mechanisms for jamming . in mode - coupling theory , jamming results from an entropy crisis : at low temperature , the system supports any applied stress because no motion of molecules is allowed . resistance to shear is therefore an intrinsic property of a material , determined by the thermodynamic parameters of the system , and results from a phase transition . in stz theory , jamming results from a structuration of the contact network : jamming requires some amount of plastic deformation to induce the creation of anisotropic structures that can support a static stress . those two mechanisms are shown here to be incorporated in a unique theory . the difference between hard and soft glassy materials will be identified with the cases when one or the other process dominates . the interplay between those two mechanisms lead to a novel interpretation of brittleness and ductility , and provide age - dependent dynamical yield criteria . the main assumption made in this work is that the out - of - equilibrium dynamics of a glass can be modeled in the framework of a mean - field approximation . by mean - field , i mean that spatially averaged quantities account for the local state of the material . constitutive equations are proposed that relate the average shear stress , to the average shear strain . the internal state of the material is determined by free - volume and by densities of arrangements ( to be defined further ) . the forces in the material are supposed to be characterized by a stress tensor , which is written @xmath0 where @xmath1 is the pressure , and @xmath2 the shear stress . the deformation tensor is defined accordingly , and composed of shear , and of isotropic deformations , directly related to dilatancy . the average free - volume per molecule is denoted @xmath3 , and is related to the total volume @xmath4 by , @xmath5 where @xmath6 is the volume per molecule of the material in random close packing , and @xmath7 is the number of molecules .    under strong shear , soft glassy materials , like clays , foams , or granular materials , can achieve large deformations ; the evolution of the system is more similar to the flow of a liquid than to the deformation of a solid . in this work i present a theory of deformation that is expected to hold for large deformations of amorphous solids or soft glassy materials . since there is a long - standing matter of debate in the mechanical community about linear versus non - linear deformations ,  @xcite it seems necessary to clarify this issue . i show here that the notion of inherent states  @xcite offers a simple physical interpretation of small and large deformations .    at low temperature , a supercooled liquid or a glassy material evolves in a complex energy landscape which can be partitioned in domains of influence of local minima . for almost all geometrical configuration of the molecules , a unique inherent state is defined as the local minimum to which the systems relaxes if suddenly quenched to zero temperature ,  @xcite and in the absence of forcing . when a material is forced , a constant stress or strain pulls the system through the energy landscape . each geometrical configuration reached lies in the domain associated with an instantaneous inherent state . the time series of instantaneous inherent states provides a coarse - grained description of the deformation of the material , and is taken here as the definition of irreversible , plastic , deformation . the instantaneous inherent state is a natural reference state for a given molecular configuration , and deviation from the instantaneous inherent state defines an _ inherent ( shear ) deformation _ starting from a known initial configuration , the total shear deformation , @xmath9 , results both from inherent deformation , and from the _ flow _ of inherent states in the phase space . by definition , the deformation @xmath10 is measured between two inherent states : it is not a state variable , since it is defined only from the knowledge of a configuration designed as `` initial '' . @xmath8 , however , is a state variable , since it is well defined for any instantaneous configuration of the system . moreover , @xmath8 is expected to be small , since it measures the deformation between two nearby configurations in phase space ; elasticity results from a harmonic approximation for every local minimum : it permits to write the stress as given by a hooke law , @xmath11 . the equation that governs @xmath2 can finally be written : @xmath12 the rate of plastic deformation @xmath13 will come out of a statistical analysis of hopping motion in the phase - space . it will be written as a function of the state variables of the system , @xmath2 , @xmath3 , an others , to define the constitutive law of the material . this work relies on the idea that , in dense materials , macroscopic deformation results from free - volume activated rearrangements . this is expected because temperature is so low that thermal activation becomes irrelevant , and activated processes are controlled by entropic fluctuations of the free - volume . a rearrangements occurs if sufficient space exists in the neighborhood of a given point ; the transition probabilities are thus determined by the size distribution of voids . in order to work at a mean - field level , the volume distribution is characterized by the single scalar value @xmath3 ; in order to write activation factors , it is necessary to provide an ansatz for the distribution of voids in the material .    in the original works by cohen and turnbull on free - volume theory , it was argued that the voids in a material are given by a poisson distribution . @xcite in their picture , the motion of a molecule in a glassy material is permitted only if it can move into a hole of size , say @xmath14 . the probability to find a hole larger than @xmath14 , is proportional to @xmath15 : this leads to estimate a diffusion constant @xmath16 , which is a basis of free - volume approaches . important differences , however , exist between various uses and interpretations of free - volume theory . falk and langer use free - volume activation factors in the definition of stz theory . @xcite their argument borrowed from original ideas by edwards and coworkers for granular materials :  @xcite the only extensive variable is supposed to be the volume @xmath4 ; the number of states available to the system is roughly proportional to @xmath17 , where @xmath18 is an arbitrary constant with dimension of a volume ; the entropy is defined as , @xmath19 and the intensive quantity , analogous to temperature , is @xmath20 : @xmath21 free - volume thus enters activation factors of the form @xmath15 .    in stz theory , this discussion turns out to be somewhat formal since free - volume is taken constant : @xmath3 is absorbed in the phenomenological constants of the resulting equations . it also leads to complicated expressions for transition rates as functions of the applied stress ( because @xmath14 is supposed to be a function of @xmath2 ) . however , free - volume should vary when the material dilates or contracts . this idea has been considered recently by several authors in the analysis of density relaxation of granular materials subjected to vertical tapping . @xcite in these works , free - volume is understood as a purely dynamical quantity , and relaxes by the motion of single grains in holes created by volume fluctuations . a poisson distribution of volume fluctuations is assumed at all times : this leads to an equation of motion of the form , @xmath22 \quad . \label{eqn : vf:0}\ ] ] the activation factor @xmath23 $ ] is the probability that a hole of size larger than @xmath14 is present at a given point in the material ; @xmath24 is determined by the tapping frequency . the current work borrows from those different lines of thought . drawing on the works initiated on granular compaction , free - volume is a dynamical quantity ; voids are expected to redistribute fast enough in the material ; fast means , faster than free - volume activated collective rearrangements themselves . this is justified because redistribution of free - volume requires very tiny displacements of the molecules , and voids , considered as particles , are expected to diffuse faster than molecules themselves : around a void , there is by definition an excess of free - volume . there is an important difference , however , from the picture used in  @xcite where compaction results from the motion of a single grain in a hole . here , borrowing on the argument of stz theory , compaction results from elementary rearrangements , involving several molecules at a mesoscopic scale . probabilities of transitions are estimated by entropic arguments . due to the redistribution of voids , the number of states available to the system is supposed to be proportional to @xmath17 at all times : the fast redistribution of voids allows the system to realize a quasi - equilibrium thermodynamic ensemble , determined by an intensive variable , @xmath3 . the existence of an intensive variable , associated to the fluctuations of molecular configurations , is supported by the observation of an effective fluctuation - dissipation theorem in sheared fluids . @xcite for reasons of clarity , and to be specific , i reserve the word temperature for the intensive quantity associated with the fluctuations of energy , or for the specific kinetic energy;therefore , i will not comply with to the recent use of `` effective temperature '' for something that is neither thermodynamic nor granular temperature . nevertheless , the observation of an effective fluctuation - dissipation theorem indicates that such an intensive quantity exists , without identifying what this quantity actually is . here , it is identified to @xmath3 and is associated with entropic fluctuations of molecular configurations at a mesoscopic scale . the important novelty introduced in the current work it that @xmath3 enters activation factors and is also a dynamical quantity . this is in weak analogy with granular temperature , which is a temperature , and evolves dynamically . but granular temperature couples with properties of the material only as far as it determines a collision frequency ,  @xcite while free - volume enters exponential prefactors in transformation rates . another important element of the current theory is that the dynamics of free - volume results both from density relaxation by activated rearrangements  which exists in the absence of forcing  and from dilatancy induced by the macroscopic shear flow  which occurs when the system is driven out - of - equilibrium . the non - linear coupling with the mean - flow will be addressed later ; i will first continue the description of rearrangement processes . the typical elementary compaction process , corresponds to the system going through a saddle point in phase space . the rate of such transformations is controlled by free - volume fluctuations and this saddle is characterized by a typical activation volume @xmath14 , which depends on the shape of the molecules , and on the details of microscopic interactions . the activated hopping of the system through this type of saddle permits density relaxation , and leads to an equation of motion for @xmath3 of the form  ( [ eqn : vf:0 ] ) . other types of elementary rearrangements occur in the material , in particular leading to elementary shear motion . before presenting a more elaborate scheme for macroscopic shear deformation , let me use a very minimal argument : shear strain is supposed to be proportional to shear stress : @xmath25 ; where @xmath26 is proportional to a free - volume activation factor . the current approach relies on the remark that the relative motion of molecules that permit compaction and shear are different ; different types of collective rearrangements , should require different activation volumes . in phase space , shear and compaction correspond to two types of saddle points , and those saddle points are entropic barriers of different heights . the activation volume for an elementary shear is denoted @xmath27 , and _ a priori _ differs from @xmath14 . this leads to @xmath28\,\sigma \quad . \label{eqn : dotgamma:0}\ ] ]    one major purpose of this work is to study the consequences of introducing complexity , not through a non - trivial distribution of timescales ,  @xcite but through the existence of two types of entropy barriers of unequal heights . the set of equations  ( [ eqn : sigma:0 ] ) ,  ( [ eqn : vf:0 ] ) , and  ( [ eqn : dotgamma:0 ] ) will be shown to account for logarithmic relaxation of @xmath3 and for a kww relaxation of the stress after a strain increment . [ [ shear - transformation - zones . ] ] shear transformation zones . + + + + + + + + + + + + + + + + + + + + + + + + + + +    equation  ( [ eqn : dotgamma:0 ] ) accounts for a free - volume dependent viscosity : the material described is a liquid ( in fact a non - newtonian liquid ) . it is however , a very rough description of shear motion . the theoretical description of shear deformation that i present is directly inspired by stz theory . @xcite the first assumption on which stz theory is based is that elementary shear deformation does not occur just anywhere in the material , but in zones , where a rearrangement is permitted by the local configuration of molecules . a shear transformation zone ( stz ) is thus defined as a locus within a material , where an elementary shear is possible . another essential assumption is that stz s are two - state systems . this can be understood as follows : an elementary arrangement corresponds to the opening and closing of some contacts between molecules ; once such a rearrangement has occurred somewhere in the material , these molecules can not shear further in the same direction , but they can shear backwards . a symmetry is thus introduced by shearing motion , and ( at least ) two types of arrangements can be identified that are transformed into one another , as pictured here :    = 0.0011    ( 100,60)(0,0 ) ( 50,40)(0,0)@xmath29 ( 50,0)(0,0)@xmath30 ( -50,0 )    note that an stz is not supposed to contain only four molecules , but it is a mesoscopic object that accounts for the structure of several neighboring molecules ( say of order 10 ) . to simplify the analysis , only one pair of types of arrangements is considered , aligned along the principal directions of the stress tensor . macroscopic motion results from the statistics of elementary rearrangements , and reads : @xmath31 where @xmath32 denote the number density of arrangements , and where @xmath33 denote the rate of transition @xmath34 . those rates result from activation processes determined by free - volume , but also by force fluctuations . a stz shears if there is sufficient free - volume at its location , and if the local bias of the force network triggers the shear in the appropriate direction . in this work , force and volume fluctuations are supposed to be uncorrelated ; the probabilities associated with those two types of fluctuations factorize : @xmath35 . from the preceding discussion , @xmath36 . in order to estimate the factors @xmath37 associated with the fluctuations of the force network , let me first remark that @xmath38 is a positive , increasing function of @xmath2 . moreover , an elementary shear is expected to be triggered by _ large forces_. the distributions of static forces in glassy systems have been shown to share strong similarities with the force distribution in a granular material ,  @xcite in which case , large forces are distributed exponentially .  @xcite this justifies an exponential dependency of @xmath37 on @xmath2 : @xmath39\,\exp\left[{\pm{\sigma\over\bar\mu}}\right ] \quad,\ ] ] with @xmath40 , the update frequency of microscopic processes . note however , that the exponential form for the activation factors @xmath37 plays no role in the following , since these factors will be linearized . the parameter @xmath41 measures the typical stress that must be overcome to trigger a rearrangement . in hard - sphere systems , it is expected to be proportional to the pressure @xmath1 . in general , @xmath41 is essentially the energy required to break a bond between two molecules , times the density of bonds . here , @xmath41 will be taken constant , either because @xmath1 fixed , or because it depends on some interaction potential . at some point , @xmath41 will be set to unity , which fixes the unit of forces .    [ [ dynamics - of - arrangements . ] ] dynamics of arrangements . + + + + + + + + + + + + + + + + + + + + + + + + +    following  @xcite , the equation of motion for the populations @xmath32 is written , @xmath42 the first two terms of the rhs account for the transitions between both types of stz . the last term accounts for the renewal of molecular configuration by the macroscopic flow : stz constantly appear and disappear during plastic deformation . the renewal of configurations is supposed to be proportional to the work of plastic deformation , @xmath43 , which is thus a common factor of creation and annhilitation terms . [ [ shear - induced - dilatancy . ] ] shear - induced dilatancy . + + + + + + + + + + + + + + + + + + + + + + + +    in this work , the dynamics of rearrangements is coupled to free - volume dynamics . part of this coupling results from the factor @xmath44 that controls shear transformation , but the macroscopic flow is also expected to influence the relaxation of free - volume . free - volume dynamics is derived in analogy with the dynamics of the populations @xmath32 : plastic deformations corresponds to the evolution of the system in the phase - space along trajectories that do not necessarily permit a minimization of free - volume . in real space , plastic deformation requires the constant renewal of local arrangements , and new arrangements are local configurations of large entropy ; by definition , their density is not optimized . therefore , the renewal of configurations by the shear flow leads to dilatancy . in experiments on suspensions , for example , the system is fluidized by applying a strong shear which rejuvenates the glassy structure . @xcite    in order to quantify the average dilatancy induced by the shear flow , i assume that a fraction of the work of plastic deformation , @xmath45 in used for enthalpy production , @xmath46 . the equation for @xmath3 is : @xmath47 + { { \cal a}_v\over p}\,\sigma\,\dot\gamma \label{eqn : vf:1}\ ] ] this equation is an essential element of the current theory . the non - linear coupling with the macroscopic flow shows up when the system is strongly driven out - of - equilibrium . the competition between self - relaxation and shear induced dilatancy will be shown to lead to a power law viscosity .      from  ( [ eqn : stz : npm ] ) , the equation of motion for the total density , @xmath48 , of arrangements reads , @xmath49 the creation and destruction terms account for an entropic mixing of molecular configurations . the previous equation shows that this mixing leads to a convergence of @xmath50 to an asymptotic value , @xmath51 ; therefore , @xmath52 defines a state of maximal entropy of molecular configurations . for this reason , i do not expect @xmath50 to be different from @xmath53 unless in very special circumstances , _ e.g. _ if some level of crystallization exists in the material . in particular , an amorphous solid is usually produced by a rapid quench of a liquid at low temperature : the system comes from a state of high entropy , therefore , in the resulting glass , @xmath50 is expected to achieve its asymptotic value . in this work , i assume that @xmath52 at all times . the dynamics of local arrangements is thus determined only by the bias @xmath54 between the populations of stz . equations  ( [ eqn : stz:0 ] ) and  ( [ eqn : stz : npm ] ) are written in a more suitable form by introducing the variable @xmath55 and parameters , @xmath56 , @xmath57 , and @xmath58 . it is also convenient to introduce @xmath59 , @xmath60 and @xmath61 , and to rescale free - volume as , @xmath62    the complete set of constitutive equations is finally : @xmath63\ , \left(\sinh\left({\sigma\over\bar\mu}\right)-\delta\,\cosh\left({\sigma\over\bar\mu}\right)\right)\\ \label{eqn : stzdil:2 } \dot\delta&=&{\ddpl\over\epsilon_0}\,\left(1- \mu_0\,\sigma\,\delta\right)\\ \label{eqn : stzdil:3 } \dot\chi & = & -e_1 \exp\left[-{\kappa\over\chi}\right ] + \alpha\,\sigma\,\ddpl \quad.\end{aligned}\ ] ] those equations involve two state variables , @xmath64 and @xmath20 which play very different roles . @xmath64 evolves only when the material is sheared , therefore accounts for some sort of structural memory induced by deformation imposed on the material . on the contrary , @xmath20 evolves spontaneously , and relaxes towards 0 . since this relaxation is slow , @xmath20 also encodes some sort of memory , but a very different type of memory than @xmath64 , a memory which constantly evolves as the system ages . note that the asymptotic value for the relaxation of @xmath20 ( which here is 0 ) may , in fact , be temperature dependent ; this question is eluded in the current work : it seems a very reasonable assumption for materials like suspensions if not for glasses below the glass transition , because the logarithmic slow - down of density relaxation around @xmath65 prevents it to reach asymptotic values far below ( @xmath66 ) . the parameters involved in the theory should depend on pressure and temperature , and this dependency results from particularities of a given material ; for example the update frequencies @xmath24 and @xmath40 are expected to be proportional to @xmath67 ( where @xmath68 is the temperature ) in a dense hard - sphere material . @xcite such questions , however , are not the focus of this study ; in this article , the values of the parameters are constant and can be understood as a given thermodynamical state of the system . equations  ( [ eqn : stzdil:1]-[eqn : stzdil:3 ] ) present non - linear couplings between shear motion and density relaxation . before studying the complete set of equations in various tests , i start with two simple cases , when one or the other process dominates . i show here that jamming transitions are captured by equations  ( [ eqn : stzdil:1]-[eqn : stzdil:2 ] ) , and free - volume relaxation by equation  ( [ eqn : stzdil:3 ] ) .      in the absence of shear stress , @xmath69 and @xmath70 ; the system  ( [ eqn : stzdil:1]-[eqn : stzdil:3 ] ) reduces to the single equation : @xmath71 \quad,\ ] ] which accounts for density relaxation .    at short time , this process is dominated by the initial value , @xmath72 : @xmath73 \quad;\ ] ] this defines a time scale for short - time free - volume relaxation : @xmath74 \quad.\ ] ]    the long time behavior is estimated from the integral expression : @xmath75 \d \chi \simeq   \kappa\,\exp\left[{\kappa\over \chi(t)}\right]\ ] ] at long time , the integral is dominated by the small values of @xmath20 , _ i.e. _ by the ultimate value reached , @xmath76 . therefore , @xmath77 free - volume relaxes logarithmically with time : in the absence of any forcing , the system ages .      in order to isolate the stz mechanism for jamming , i now assume that free - volume is constant , at some value @xmath20 ; the evolution of the system is governed by equations  ( [ eqn : stzdil:1 ] ) and  ( [ eqn : stzdil:2 ] ) only , and a constant stress @xmath2 is applied . these equations admit two stationary solutions : a jammed state for which , @xmath78 , and @xmath79 a steady deformation regime , for which , @xmath80 and , @xmath81\ , \left(\sinh\left({\sigma\over\bar\mu}\right)-{1\over\mu_0 \sigma}\,\cosh\left({\sigma\over\bar\mu}\right)\right ) \quad.\ ] ] the latter solution is unstable at small stresses , in which case , the material jams . for large stresses , the jammed state loose stability , and steady deformation results . the two solutions exchange stability at the yield stress @xmath82 , which is the solution of : @xmath83      in the following , i will consider several limiting cases of the general constitutive equations  ( [ eqn : stzdil:1]-[eqn : stzdil:3 ] ) . in particular , the exponential dependency of @xmath84 as a function of @xmath2 will not be studied , although this non - linearity might be important in some cases . @xcite activation factors involving @xmath2 will thus be linearized ; this is valid under the assumption that @xmath85 ; @xmath41 is then incorporated into constants by taking it to unity . it leads to : @xmath86\,\left(\sigma-\delta\right)\\ \label{eqn : stzdil : lin:2 } \dot\delta&=&{\ddpl\over\epsilon_0}\,\left(1- \mu_0\,\sigma\,\delta\right)\\ \label{eqn : stzdil : lin:3 } \dot \chi & = & -e_1 \exp\left[-{\kappa/\chi}\right]+\alpha\,\sigma\ddpl \quad.\end{aligned}\ ] ]    i will also consider an _ isotropic ( liquid ) limit _ , where the structural anisotropy , measured by @xmath64 plays no role . this is the case either at timescales when the dynamics of @xmath64 is too slow , @xmath87 ( and @xmath88 ) or because the typical stress @xmath89 vanishes . in this case , the rheology is determined by a single state variable : @xmath90\,\sigma\\ \label{eqn : stzdil : iso:2 } \dot \chi & = & -e_1 \exp\left[-{\kappa/\chi}\right]+\alpha\,\sigma\ddpl \quad.\end{aligned}\ ] ] the sets of equations  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) and  ( [ eqn : stzdil : iso:1]-[eqn : stzdil : iso:2 ] ) will be helpful to emphasize the roles of @xmath64 versus @xmath20 . many aspects of glassiness will be shown to be captured by equations  ( [ eqn : stzdil : iso:1 ] ) and  ( [ eqn : stzdil : iso:2 ] ) . jamming , however , is not captured at this level of simplification , and requires the existence of the variable @xmath64 . a word about the preparation of the sample : at time @xmath91 , the system is quenched from a highly mixed state , in which , @xmath88 and @xmath20 is large . the sample then ages during a waiting time , @xmath92 without any forcing . after this relaxation , some test is performed on the system . two main experimental procedures consist in forcing the stress , or forcing the deformation . if a stress , @xmath2 , is applied , the plastic deformation @xmath93 is determined by constitutive equations , and the total deformation is , @xmath94 . if the deformation , @xmath95 , is forced , the evolution of @xmath2 is given by equation  ( [ eqn : sigma:0 ] ) , and couples to constitutive equations .    if a small perturbation is applied to the sample , the response of the system results from first order terms in the linearized constitutive equations  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) , or  ( [ eqn : stzdil : iso:1]-[eqn : stzdil : iso:2 ] ) . in relaxation processes , however , the linear approximation is not always valid at long times ,  @xcite because , the activation factors can become small , possibly smaller than non - linear terms . in the current theory , all non - linear terms involve the work of plastic deformation @xmath96 . therefore , the linear analysis is valid so long as this work is small compared to relaxation processes . stress relaxation is monitored after a small strain step . the step can be understood as a ramp in the deformation , which occurs on a very small time interval around @xmath92 . the strain @xmath97 determines the initial value of the stress , @xmath98 . during the relaxation , @xmath99 ; the subsequent evolution of the stress is governed by @xmath100 , and accounts for the plastic adaptation of the material to the imposed strain . the relaxation modulus , @xmath101 is the correlator associated with rearrangement processes : it measures the correlation between the initial state of the material , and its state at any following time , @xmath102 . assuming that the diffusion of molecules is dominated by collective rearrangements , stress relaxation can then be compared with the decay of self - intermediate scattering functions , which measure the spatial decorrelation resulting from the diffusion of particles in the dense medium . @xcite the long - time behavior of correlation functions is one of the major properties of glassy materials . @xcite ( see also recent reviews by debenedetti and stillinger ,  @xcite or angell _ et al_.  @xcite ) .      in the isotropic limit , the constitutive equations are  ( [ eqn : stzdil : iso:1 ] ) and  ( [ eqn : stzdil : iso:2 ] ) . stress relaxation is governed by the linearized equations : @xmath103\sigma\\ \dot\chi & = & -e_1 \exp\left[-{\kappa\over\chi}\right ] \quad.\end{aligned}\ ] ] to determine the relaxation modulus , @xmath104 , the only relevant initial condition is the initial value of @xmath20 , when the strain is applied ; this value is also determined by the initial free - volume @xmath72 at time @xmath91 and by the waiting time @xmath92 . moreover , either @xmath105 or @xmath106 can be incorporated in the time scale , therefore , only two effective parameters determine all the possible outcomes of those equations : @xmath65 , and @xmath107 . while @xmath65 measures the relative heights of entropy barriers , @xmath107 measures the relative update frequencies of shear and compaction processes .      at long time , the relaxation of @xmath20 is given by  ( [ eqn : vf : relax ] ) , and stress relaxation verifies : @xmath108\sigma \simeq -\mu\,e_0\ , \left({e_1\,t\over\kappa}\right)^{-1/\kappa}\sigma \quad.\ ] ]    the response modulus @xmath109 is , @xmath110\ ] ] with an exponent @xmath111 which is directly related to the ratio @xmath65 between the heights of entropy barriers : @xmath112    the validity of linear approximation is checked by considering the time - dependency of the work of plastic deformations @xmath96 ; this term should be compared with the first term is free - volume relaxation , @xmath113\propto { \kappa\over t } \quad.\ ] ] the stress relaxation , leads to the following evolution for the work of plastic deformations : @xmath114 \quad.\ ] ] if @xmath115 , or @xmath116 , the relaxation of @xmath96 is dominated by the exponential factor ; if @xmath117 , the stress goes to a constant , but @xmath118 decays faster than @xmath119 . in all cases , the non - linear term decays faster than the linear term , which justifies the linear approximation at all times .    for @xmath116 , _ i.e. _ @xmath115 , the stress undergoes kww relaxation . this property results directly from the existence of two types of saddle points of unequal heights , when hopping motion is controlled by a logarithmically relaxing free - volume . measured values of exponent @xmath111 in glassy materials , ranging from @xmath120 to @xmath121  @xcite are consistent with , @xmath122 $ ] . the ratio @xmath65 is related to the shape of molecules , and to the details of their interactions . there is no _ a priori _ reason , from a phase space point of view , to assume a particular value for @xmath65 . in this work i will systematically study both cases in order to draw a complete picture of the phenomenology described by the proposed constitutive equations .    if @xmath117 , then @xmath123 , and the previous calculation shows that the stress assumes a non - vanishing asymptotic value after relaxation . in this case , @xmath124 , and the stress relaxation towards its ultimate constant value can be further estimated as a power law : @xmath125\\ & \simeq & \exp\left[-|a|\,t_w^{-|\beta|}\right]\;\left(1+|a|\,t^{-|\beta|}\right ) \quad.\end{aligned}\ ] ] interpreting @xmath126 as a correlation function , it also means that the material never completely decorrelates from its initial state . this non - ergodicity is an indication of the glass transition as understood in mode - coupling theory  @xcite ; it can be characterized by the asymptotic value @xmath127 , the edwards - anderson parameter : @xmath128 \quad,\ ] ] which increases with the age of the sample . such a dependency of the edwards - anderson parameter of the age of the sample has been observed experimentally by bonn and coworkers on a colloidal glass at very low concentrations . @xcite      at short time , if free - volume can be assumed large or constant , the exponential term @xmath129 $ ] does not vary very much , and stress undergoes an exponential relaxation : @xmath130\;\left(t - t_w\right)\right ] \quad.\ ] ] this relaxation occurs on a time scale of order , @xmath131 \quad,\ ] ] which should be compared with the time scale of free - volume relaxation , @xmath132 . if @xmath133 , free - volume relaxes faster than stress , and the long time , free - volume controlled , stress relaxation can be observed . if @xmath134 , the stress relaxes exponentially before free - volume achieves sufficiently low values to slow down stress relaxation . the crossover between exponential and kww relaxation is determined by , @xmath135 or , @xmath136 \sim   { \chi_0\over e_1}\;\exp\left[{\kappa\over\chi_0}\right ] \quad.\ ] ] for a large initial free - volume , this reduces to , @xmath137 kww relaxation can be observed for a larger value of @xmath105 , or for sufficiently small initial value of the free - volume , @xmath138 , ( that is for sufficiently large waiting time @xmath92 ) .    in order to look at the short and long time response of the system , equations  ( [ eqn : sigma:0 ] ) ,  ( [ eqn : stzdil : iso:1 ] ) and  ( [ eqn : stzdil : iso:2 ] ) are integrated numerically with a runge - kutta algorithm with adaptive timestep . the system is prepared in a dilute state , @xmath139 , which is allowed to relax during time @xmath92 , in the absence of any forcing ( @xmath140)=0 $ ] ) . at time @xmath92 , a strain @xmath141 is applied and produces an elastic response , @xmath142 . the dynamics of free - volume and stress are then monitored . [ [ crossover - stress . ] ] crossover stress . + + + + + + + + + + + + + + + + +    in a first series of tests , @xmath92 is set to 0 , and the relaxation spectra are displayed for different values of the parameter @xmath105 . parameters @xmath143 , @xmath144 and @xmath145 are set to 1 ; @xmath146 and @xmath147 . the results are displayed figure  [ fig : relax1 ] . for the first value @xmath148 , only exponential relaxation is observed because all the stress has been released before @xmath20 entered the logarithmic relaxation . the value @xmath149 corresponds to the crossover @xmath150 : kww relaxation is observed at later stages of the relaxation , with a non - vanishing values of @xmath126 . for a given @xmath72 , the stress release during the exponential part of the relaxation is proportional to @xmath106 , hence inversely proportional to @xmath105 . the plateau in the relaxation spectrum can be observed only for sufficiently large values of @xmath151 so that @xmath126 is still observable by the time @xmath20 enters the logarithmic relaxation . these relaxation curves are very similar to those observed in md simulations ,  @xcite or in experiments ;  @xcite the general appearance of those curves depends sensibly on the parameters considered . from a practical point of view , it means that it depends on the time window that is available , numerically or experimentally . = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath152 ( 5,50)(0,0)@xmath126 ( 90,2)(0,0)@xmath153    [ [ aging ] ] aging + + + + +    the crossover between exponential relaxation can also be studied by varying the initial free - volume . this amounts to varying the waiting time @xmath92 for a given @xmath154 . in order to observe a plateau in the relaxation spectrum , parameter @xmath105 is set to 20 , and parameters @xmath143 , @xmath144 and @xmath145 are set to 1 ; @xmath146 . in order to obtain an exponential relaxation for short waiting times , and to observe a wide variety of behavior , the initial value of @xmath20 is set to 100 . the results are displayed figure  [ fig : relax2 ] . with this choice of parameter , the curves resemble experimental data . @xcite    such curves strongly depend on the available experimental window , and on the width of the crossover region which is unobservable in a log - lin plot for smaller values of @xmath151 . to emphasize this point , another set of relaxation spectra is presented figure  [ fig : relax ] for values of the parameters : @xmath155 , @xmath156 , and @xmath146 . the initial condition is @xmath157 and waiting times are : @xmath158 . the relaxation spectra are presented in log - lin , log - log , and log - log(log ) plot ; the crossover between exponential and kww relaxation is observable around values of the relaxation modulus of order @xmath159 , which is featureless in a log - lin plot . = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath152 ( 5,50)(0,0)@xmath126 ( 90,2)(0,0)@xmath160    = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath126 ( 5,50)(0,0)@xmath161 ( 90,2)(0,0)@xmath160    [ [ ergodic - to - non - ergodic - transition . ] ] ergodic to non - ergodic transition . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    i now consider the dependency of relaxation spectra on the parameter @xmath65 . this parameter is of primary importance in the current approach and is directly related to the existence of entropy barriers of unequal heights in the phase space . the parameters are taken to be , @xmath155 , @xmath162 , the initial condition is @xmath157 , and the waiting is set to @xmath163 . this choice is primarily motivated to obtain easily observable features on the relaxation spectrum . the results are displayed figure  [ fig : relax3 ] . when @xmath65 goes to 1 , the system presents a transition to non - ergodic behavior . for @xmath117 , the asymptotic value of the stress is no longer vanishing . at the point @xmath164 , the long time relaxation of @xmath2 verifies : @xmath108\sigma \simeq -\mu\,{e_0\over e_1\,t}\sigma \quad,\ ] ] whence an asymptotic power law relaxation : @xmath165    = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath152 ( 5,50)(0,0)@xmath126 ( 90,2)(0,0)@xmath160      kww relaxation is also achieved with the complete set of equations  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) . at time @xmath91 , @xmath166 , and no constraint is applied until time @xmath92 : @xmath140 ) = 0 $ ] . during the ramp , the deformation of the material is essentially elastic , no rearrangement occurs , hence , @xmath167 , while the imposed strain determines the initial value the stress : @xmath142 .    during the ensuing relaxation , @xmath100 and , for a small deformation , @xmath168 . it results that , at all times , the quantity @xmath169 is constant : @xmath170 here , @xmath171 , but this expression may vary , depending on the preparation of a sample . to treat initial values of @xmath2 and @xmath64 in the greatest generality , the initial value for the difference @xmath172 is denoted @xmath173 ( and here , @xmath174 ) .    at long times , the dynamics of @xmath64 and @xmath2 verify , @xmath175\left(\sigma-\delta\right)\\ & \simeq & -\left(\mu+{1\over\epsilon_0}\right )   \,e_0\,\left({e_1\,t\over\kappa}\right)^{-1/\kappa}\left(\sigma-\delta\right)\end{aligned}\ ] ] whence , @xmath176\ ] ] with , @xmath177 the evolution in time of the quantity @xmath172 is governed by the same equation as the dynamics of @xmath2 in the isotropic limit . therefore , all the previous discussion presented in the isotropic limit pertains to the general case , up to a constant term , which is fixed by the conservation of the quantity @xmath169 . the dynamics of @xmath64 and @xmath2 follow .    in particular , the evolution of @xmath64 reads , @xmath178\over 1+\mu\epsilon_0 } \quad,\ ] ] and the relaxation modulus verifies , @xmath179\over 1+\mu\epsilon_0 } \quad.\ ] ] the isotropic expression for @xmath180 is recovered by taking @xmath181 . note that the relaxation of the stress could lead to a negative value if @xmath182 , which depends on the preparation of the sample .    for @xmath116 ( _ i.e. _ @xmath115 ) , the stress undergoes a kww relaxation towards a non - vanishing value , which depends on the properties of the material and on the preparation of the sample . for @xmath123 power law relaxation is observed towards a non - vanishing value of the stress . in the current work , two origins are therefore identified for the existence of a non - vanishing stress at long times : the creation of mechanical anisotropy of the contact network ( jamming ) , or the existence of a ergodic to non - ergodic transition for the hopping motion in phase space ( freezing ) . recent observations by bonn _ et al _ indicate that for a given material , jamming could not be identified only to a loss of ergodicity , and that some amount of jamming could coexist with liquid - like behavior . @xcite such features are allowed in the current theory and show the importance of the stz mechanism for soft glassy materials . the validity of the linear approximation is checked by considering the long times relaxation of the rate of plastic deformation : @xmath183\ ] ] and the work of plastic deformation decays like @xmath13 since @xmath2 goes to a constant ; the non - linear term in equation  ( [ eqn : stzdil : lin:3 ] ) is negligible at all times . when a constant strain rate @xmath184 is applied , the material is driven strongly out - of - equilibrium ; the response of the system involves non - linear terms . i study here the stationary states of a sheared material , in which , @xmath185 . in the first sections , i make a very primitive stability analysis of those stationary states , by assuming that @xmath185 is fixed at all times , and that the stress is given by the constitutive relation  ( [ eqn : stzdil : lin:1 ] ) or  ( [ eqn : stzdil : iso:1 ] ) ; this allows to decouple equation  ( [ eqn : sigma:0 ] ) . the existence of an unstability in the complete set of equations accounts for the emergence of stick - slip motion in boundary lubrication . @xcite this question will be briefly discussed in section  [ sec : hopf ] , but raises so many issues that it will be tackled thoroughly in another publication ; for the sake of completeness , the calculation of the hopf criterion is given in appendix  [ app : hopf ] .      in the limit @xmath186 ( no yield stress ) , the relation between shear rate and shear stress reads ( from  ( [ eqn : stzdil : iso:1 ] ) ) : @xmath187 \quad.\ ] ] eliminating @xmath2 between equations  ( [ eqn : stzdil : iso:1 ] ) and  ( [ eqn : stzdil : iso:2 ] ) , yields : @xmath188 + \alpha { \dot\epsilon^2\over e_0}\ , \exp\left[{1\over\chi}\right ] \label{eqn : chidyn}\ ] ] which is positive iff , @xmath189 in the stationary flow , the relation between shear strain and volume reads , @xmath190 \quad.\ ] ] for future use , i introduce the critical rate of deformation , @xmath191    = 0.005    ( 110,60)(8,-8 ) ( 13,40)(0,0)@xmath192 ( 10,48)(0,0)@xmath184 ( 85,-13)(0,0)@xmath20 ( 10,-13 )    the bifurcation diagram for @xmath20 is presented figure  [ fig : vf ] . if @xmath193 , then @xmath194 at all times : free - volume diverges , the system blows up . if @xmath195 , the equation for @xmath20 admits one solution : @xmath196 the existence of a stationary value for the free - volume indicates that for any finite rate of deformation , aging stops . in fact , this is not completely exact if a vanishingly small shear rate was applied , such that the associated steady free - volume is so small , that it is not accessible at reasonable timescales . as a consequence , for very small @xmath13 , the apparent viscosity saturates to a constant value this is a common feature of complex liquids ( see _ e.g. _ @xcite ) .    in the steady state the value of the shear stress is , @xmath197 with , @xmath198 the system behaves like a power law fluid ; the exponent @xmath199 of the power law viscosity is related to the exponent @xmath111 of kww relaxation by , @xmath200 those two exponents are directly related to the slope @xmath65 of @xmath152 as a function of @xmath201 . the assumptions underlying the current approach do not permit to evaluate the activation volumes associated with density relaxation and shear motion . however , the single parameter @xmath65 determines two exponents associated with the linear and non - linear response of a strained material . the existence of such a relation constitutes a important result of the current theory , although the form given to this relation depends on non - generic features : how non - linear dilatancy is introduced in the equation of motion for free - volume ; how the update frequencies ( here taken constant ) are written . in a hard - sphere material , for example , collision frequencies depend on free - volume as @xmath152 , and the relation between @xmath2 and @xmath184 is expected to be modified by logarithmic corrections . molecular dynamics simulations of sheared fluid are consistent with power law rheologies for @xmath199 equal to 0.1  @xcite or 0.2  @xcite recent measurement of power law rheology have also been obtained on a colloidal glass , leading to values of @xmath199 ranging from 0.1 to 0.35 . @xcite a value of @xmath111 of order @xmath202 corresponds to @xmath203 , which seems in reasonable agreement with the currently available data .    the theory of structural rearrangement described here , is expected to present some degree of universality , and hold , in a given form , for wide classes of materials . what this model shows is that there is an intimate relation between time - logarithmic density relaxation , kww relaxation , and power law rheology of a dense material , and that this relation may depend only on generic features for a given class of materials . therefore the existence of a relation between @xmath111 and @xmath199 could be tested by considering families of materials . the experimental or numerical evaluation of such a relation appears as an important mean of characterization of the coupling between two major modes of deformation . the constitutive equations are now  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) . in the steady state , @xmath204 plugging expression  ( [ eqn : stzdil : lin:1 ] ) for the plastic deformation in equation  ( [ eqn : stzdil : lin:3 ] ) leads to the following relation between @xmath20 and @xmath2 in the steady state : @xmath205 = e_0e_1(^2 - 1_0 ) . @xmath205 using  ( [ eqn : stzdil : lin:3 ] ) again , @xmath184 can then be written as a function of @xmath2 : @xmath206 for large @xmath2 or @xmath184 , the system behaves as a power law fluid , with the same exponent @xmath199 as in the isotropic limit . however , @xmath184 vanishes for a non - zero value of @xmath2 , which determines the yield stress : @xmath207    the following expression can also be obtained for @xmath184 as a function of @xmath20 : @xmath208   \over \sqrt{\alpha}\ ; \sqrt{e_0\alpha\,\exp\left[{\kappa/\chi}\right ]               + e_1\mu_0\,\exp\left[{1/\chi}\right ] } } \quad,\ ] ] which is an increasing function of @xmath20 , and which has a maximum for @xmath209 : @xmath210 for @xmath211 , like in the isotropic case , the shear deformation leads to a constant increase of free - volume , and the system blows up . when a shear rate is suddenly imposed on a material , the transient dynamics often lead to a first rise of the stress before a dynamical yield stress is reached at which plastic deformations begin ; once plastic deformation occur , the stress relaxes towards its strain - rate dependent asymptotic value . @xcite these features are observed in sheared bulk solids  @xcite but also in sheared thin films ,  @xcite granular materials  @xcite or in solid friction . @xcite another generic property of transient regimes is that the dynamical yield stress depends on the age of the material : during aging , a material strengthens . these phenomena are captured by the proposed equations . strain softening and strengthening result essentially from the equation of motion for free - volume and are captured by the isotropic equations .    figure  [ fig : bump ] present some examples of transient dynamics obtained by numerical integration of equations  ( [ eqn : sigma:0 ] ) ,  ( [ eqn : stzdil : iso:1]-[eqn : stzdil : iso:2 ] ) . parameters are set to @xmath212 . the initial value for the free - volume ( at time @xmath91 ) is set to @xmath213 , and the system relaxes without forcing during a time @xmath92 . at time @xmath92 , a constant strain rate is applied . to prevent the blow - up of the system , @xmath184 must be smaller than @xmath192 , which equals @xmath214 with the chosen parameters ; for this reason , the strain rate is set to @xmath215 . two different values of the parameter @xmath65 have been used to explore the two important cases when @xmath65 is larger or smaller than 1 . the stress is displayed as a function of the total time @xmath102 since preparation of the sample . = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath152 ( 5,50)(0,0)@xmath2 ( 90,2)(0,0)@xmath102    a dynamical yield stress clearly emerges for a sufficiently aged material , and increases with time . the yielding of the material is accompanied by a sudden rise of the free - volume allowing transformation rates to display non - vanishing values . these transients result from a rather complex interplay between the dynamics of free - volume and shear transformation . the complete study of this process may lead to many developments that would divert this work from the main questions that i would like to address , and i prefer to leave the discussion here . the stability analysis of equations  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) and  ( [ eqn : sigma:0 ] ) shows the existence of a hopf bifurcation which is responsible for oscillatory response to a given strain rate . this leads to some understanding of stick - slip motion in sheared lubricants  @xcite or in granular materials . @xcite in this work , i mention this instability for the sake of completeness , but it raises many questions , and will require a more thorough consideration . the calculation of the hopf criterion of bifurcation is given in appendix  [ app : hopf ] . in the isotropic limit , the steady plastic deformation is stable iff @xmath216 , with : @xmath217 ^ 2 \quad.\ ] ] for @xmath115 , @xmath218 : there is no hopf bifurcation . for @xmath117 , @xmath219 on the interval @xmath220 $ ] , and the steady state becomes unstable iff @xmath221 . this criterion is represented figure  [ fig : hopf : iso ] . for driving strain rates above @xmath192 , the system looses stability , because there is no stationary solution for @xmath20 . below @xmath192 , the hopf criterion separates a domain of stability , above @xmath222 where the system reaches steady plastic deformation , and and domain below this curve where stick - slip motion is achieved . = 0.005    ( 110,60)(8,-8 ) ( 6,40)(0,0)@xmath222 ( 68,-13)(0,0)@xmath192 ( 85,-13)(0,0)@xmath184 ( 10,-13 )    the introduction of the variable @xmath64 in the set of equations  ( [ eqn : sigma:0 ] ) ,  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) slightly modifies this picture ( see appendix  [ app : hopf ] ) . two limit curves @xmath223 emerge , which delimit the unstable domain in the plane @xmath224 . these curves are displayed figure  [ fig : hopf : lin ] along with the previous curves @xmath222 of the isotropic limit . the unstable domain is a subspace of what it is in the isotropic limit , and the system present a creeping zone , close to @xmath99 , when steady sliding is stable for very small values of the strain rate . = 0.005    ( 90,90)(0,0 ) ( 10,5 ) ( 5,75)(0,0)@xmath144 ( 70,2)(0,0)@xmath184    some examples of stick - slip motion are displayed figure  [ fig : stickslip ] , with the set of parameters used in the phase diagram of figure  [ fig : hopf : lin ] : @xmath225 , and @xmath226 in the isotropic limit ( solid lines ) and @xmath227 for the complete set of equations ( dashed lines ) . the fast bursts of plastic deformation are accompanied by sudden jumps in the free - volume . @xmath64 does not play a major role for the definition of the dynamical yield stress , as seen from the almost identical values of the maximum stress ; @xmath64 , however , strongly influences smaller values of the stress as it shifts upwards the unstable fixed point ; this also modifies the periodicity of stick - slip motion . = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath152 ( 5,50)(0,0)@xmath2 ( 90,2)(0,0)@xmath102 in a creep test , a material is tested by forcing a constant stress . since @xmath2 is constant , the elastic deformation is @xmath228 at all times , and the plastic deformation is determined by the constitutive equations .        a small constant stress @xmath2 is applied at time @xmath229 . the plastic deformation of the material is determined by equations  ( [ eqn : stzdil : iso:1]-[eqn : stzdil : iso:2 ] ) , for a fixed @xmath2 . at large times , the rate of plastic deformation reads @xmath230 this lead to a power law variation of the compliance , @xmath231 : @xmath232 if @xmath115 , or @xmath116 , the compliance diverges at long times ; if @xmath117 , or @xmath123 , the compliance saturates at long times : @xmath233 the system jams . this jamming is purely entropic , and does not result from a structural anisotropy of the force network . note that the asymptotic value of the compliance decays with the increasing age of the sample . the validity of the linear approximation must be checked by considering the evolution of the work of plastic deformation : @xmath234 the linear approximation is justified so long as this term is negligible , when compared to the first term in free - volume relaxation  ( [ eqn : stzdil : iso:2 ] ) , which decays as @xmath235 . when @xmath117 , the work @xmath45 decays faster than @xmath235 , and therefore , the linear approximation is valid at all times ; when @xmath115 , however , this approximation fails when : @xmath236 or , @xmath237 the time @xmath238 is a typical time after which the linear approximation fails ; this time diverges when @xmath239 . it will appear clearer in the following that , when @xmath115 , free - volume presents a non - vanishing stationary value for any small applied stress . this asymptotic value is determined by the balance between density relaxation and shear induced dilatancy . the failure of the linear approximation , for @xmath115 , is due to the existence of a non - zero stationary value of @xmath20 for any applied stress .      in the general case , the constitutive equations are  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) . the constant applied stress deforms the material elastically , and at time @xmath240 , no plastic rearrangement has yet occurred . note that the initial value of the state variable @xmath64 depends on the preparation of the sample . in the experiments by cloitre _ et al _ , for example , the sample is prepared by applying a strong shear stress , above the yield stress during a rather long time . @xcite in this case , the initial value of @xmath64 is expected to be non - vanishing , and this leads to their observation of strain recovery . for the sake of generality , i will consider that the initial value of @xmath64 is non - vanishing : @xmath241 . in this section , i will assume that this value of @xmath64 is small enough so that the linear approximation is still valid . when plastic rearrangements occur in the material , they induce a structuration of the material , and the variable @xmath64 starts to evolve as @xmath168 . the plastic rate of deformation evolves as , @xmath242 and @xmath64 as , @xmath243    from this equation , the relaxation of the variable @xmath172 is easily obtained : @xmath244 \quad,\ ] ] with @xmath245 when @xmath115 , @xmath172 undergoes a kww relaxation towards 0 : the material relaxes towards a jammed state in which the stress is supported by structural anisotropy . when @xmath117 , @xmath172 saturates to a non - vanishing value : at long times , a fraction of the stress is supported by the bias between arrangements , measured by @xmath64 , but another fraction is supported by the loss of ergodicity resulting from the freezing of elementary shear processes .    from equation @xmath168 , the compliance reads , @xmath246 ) \quad.\ ] ] the deformation relaxes to a constant value , which depends on the preparation of the sample ; the relaxation is either a kww or a power law relaxation depending on @xmath65 . note that the power law variation of the compliance that was found in the isotropic limit is recovered when @xmath247 ; for finite @xmath248 , this power law behavior shows up at short times , when @xmath64 is still far from saturation .      in this section i start by considering the complete set of equations  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) , and in particular , i recall how the dynamics of @xmath64 lead to jamming . then i show how this leads to an effective equation for free - volume which is valid in the isotropic limit or in the general case . a constant stress @xmath2 is applied . the dynamics of @xmath64 is governed by the following equation : @xmath249\ , \left(\sigma-\delta\right)\;\left(1-\mu_0\,\sigma\,\delta\right ) \quad . \label{eqn : deltadyn}\ ] ] since the normalized free - volume @xmath20 enters this equation only through a positive factor , the time - dependent value of @xmath20 does not change the stability analysis performed at any fixed free - volume . the bifurcation diagram for @xmath64 is drawn figure  [ fig : delta ] . = 0.005    ( 110,60)(8,-8 ) ( 12,40)(0,0)@xmath64 ( 85,-13)(0,0)@xmath2 ( 10,-13 )    the dynamics of the variable @xmath64 leads to jamming for any stress smaller than the yield stress  @xcite @xmath250 in which case , @xmath251 and @xmath99 . above the yield stress , the steady shear flow becomes a stable solution , in which case , @xmath252 and @xmath253\,\left(\sigma-{1\over\mu_0\,\sigma}\right ) \quad . \label{eqn : eps : creep}\ ] ]    if @xmath254 , the systems jams , and the second term in the equation of motion for the free - volume vanishes . free - volume continues to relax logarithmically in time while the system does not deform . when @xmath255 , the system presents a non - zero plastic flow , and the evolution of free - volume is governed by the following asymptotic equation : @xmath188 + \alpha { e_0}\,\exp\left[-{1\over\chi}\right]\,\left(\sigma^2-{1\over\mu_0}\right ) \quad . \label{eqn : chi : creep}\ ] ] at the level of this equation , the jamming resulting from stz s shows up only through the term @xmath89 in equation  ( [ eqn : chi : creep ] ) , which shifts the square of the applied stress . the stability analysis pertains essentially to the isotropic limit . from this point on , i will consider the response to a constant stress as governed by equations  ( [ eqn : eps : creep ] ) and ( [ eqn : chi : creep ] ) ( which amounts to neglect the influence of @xmath64 during transients ) . i now study the free - volume dynamics which result from equation  ( [ eqn : chi : creep ] ) . the material dilates iff @xmath256 , with : @xmath257}\ ] ] and contracts otherwise . the resulting bifurcation diagram for @xmath20 is drawn figure  [ fig : chi:1 ] and  [ fig : chi:2 ] for @xmath115 and @xmath117 respectively .      for @xmath115 , free - volume dynamics admit a stationary solution for any @xmath258 $ ] , with , @xmath259 for values of the applied stress in the interval @xmath260 $ ] , the material undergoes steady plastic deformation : it is ductile . the shear deformation reads , @xmath261 and the strain rate behaves as a power of the stress for large stresses . when the applied stress varies in the interval , @xmath258 $ ] , the strain rate can display a whole interval of values @xmath262 \quad.\ ] ] for larger values of @xmath2 , above @xmath263 , free - volume diverges , leading to the break - up of the material . = 0.005    ( 110,60)(8,-8 ) ( 10,48)(0,0)@xmath2 ( 13,40)(0,0)@xmath263 ( 13,-1)(0,0)@xmath264 ( 85,-13)(0,0)@xmath20 ( 10,-13 )    the resulting dynamics for the compliance and for the free - volume are presented figure  [ fig : compl : soft ] for different values of the driving stress @xmath2 . the first part of the dynamics is determined by the initial value of free - volume and the corresponding timescale . as easily seen on the bifurcation diagram  ( [ fig : chi:1 ] ) , for small stresses , an initially large free - volume decays , and slows down the plastic deformation . this results in a plateau for the dynamics of the compliance ; a larger times , free - volume converges towards a constant value , and the dynamics of @xmath265 resumes a steady increase . these curves are to compare with measurements of particle diffusion : the particle diffusion in a dense material is expected to be dominated by collective processes ; the diffusive constant is then @xmath266 $ ] , and the average motion of a particle in a glassy medium is essentially governed by the same equation as @xmath265 . the features displayed of figure  [ fig : compl : soft ] are amazingly similar to the curves found in several recent experiments . @xcite    = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath152 ( 5,50)(0,0)@xmath267 ( 90,2)(0,0)@xmath153    if free - volume diverges , just before break - up , the current equations indicate that the strain rate should display a behavior of the form , @xmath268 however , this expression must be taken with caution , since the strongly out - of - equilibrium behavior of a material during break - up , when @xmath20 is large , is beyond the scope of the current theory . in particular , assumption has been made that the distribution of voids is close to an equilibrium state . during break - up , the distribution of free - volume is expected to display strong heterogeneities leading to the nucleation of fractures .      for @xmath117 , and for any @xmath269 , @xmath270 : free - volume relaxes to 0 . in this case , the applied stress does not stop aging , shear deformation is vanishingly small : the material creeps . at long time , the linear term dominates free - volume relaxation : the calculation of linear response that has been presented in section  [ sec : lin : iso ] is valid , and indicates that compliance saturates : the system jams . in general , part of this jamming is supported by the structuration of the material , which lowers the value of the stress that enters equation  ( [ eqn : chi : creep ] ) ; the remaining stress is supported by entropic freezing of the rearrangement dynamics . the sign of free - volume relaxation can be inverted only for larger stresses , above @xmath271 ; it this case , the response depends sensitively on the value of @xmath2 and on the state of the material at the time @xmath92 when stress is applied : a larger stress is required to break an older material . the stationary state is unstable : when @xmath2 is large enough to invert the sign of free - volume dynamics , the material breaks up . in this case , there are two types of response : creep for low stress , and break - up for large values of the stress . the material is brittle . = 0.005    ( 110,60)(8,-8 ) ( 10,40)(0,0)@xmath2 ( 13,20)(0,0)@xmath263 ( 85,-13)(0,0)@xmath20 ( 10,-13 )    the resulting dynamics for free - volume and for the compliance are displayed figure  [ fig : compl : hard ] . the saturation of the compliance is clearly seen for small stresses while the divergence occurs at a value at @xmath272 . the shape of these curves changes drastically in a very small interval of @xmath2 around @xmath271 because of long transient to escape from the unstable fixed point . these features are commonly observed in creep tests . @xcite    = 0.005    ( 100,110)(5,0 ) ( 10,5 ) ( 5,95)(0,0)@xmath152 ( 5,50)(0,0)@xmath265 ( 90,2)(0,0)@xmath102      materials with @xmath115 are liable to display ductile behavior , but can break ; those with @xmath117 either creep or break - up , and could be classified as unconditionally brittle . the study of stz equations one a one - dimensional elasto - plastic decohesion has shown the importance of dynamics in ductile to brittle transition . @xcite the current analysis slightly contradicts those results , and indicates that two major types of materials can be identified , which display different qualitative behavior : this is in agreement with the common idea that some materials are ductile , some other brittle . in fact , for hard materials ( @xmath117 ) , the existence of a strong entropic crisis precludes dynamics of rearrangements to come into play : intrinsic properties of the material completely dominate the dynamical process of deformation . the existence of a dynamically driven transition between ductile and brittle behavior is still expected for soft materials ( @xmath115 ) for which the dynamics or arrangements can fully play its role . the interplay between dynamics and intrinsic glassy properties of material , as displayed in the current approach , seems a promising perspective . as far as fracture mechanics is concerned , this picture indicates that during the breaking of a hard glass , the motion of molecules in a material should be essentially localized around the surface of a crack , where there is some free - volume , whereas in the bulk , plastic deformation can occur only in places where a very large stress in concentrated . this is apparently related to a theoretical approach of fracture mechanics , recently introduced by brener and spatschek ,  @xcite where the dynamics of the fracture is controlled by diffusion of molecules along the surface . for soft glassy materials , on the contrary , the propagation of a fracture is expected to induce more easily observable plastic deformation in the bulk of the material . the current work relies on the introduction of dynamical equations for the intensive quantity associated with entropic fluctuations of molecular configurations . in this work , this quantity is identified to free - volume , but other physical pictures could be given . for example , free - volume is related to the average number of contacts per molecule ; in a foam , the deformation of the bubbles can be accounted for a effective free - volume that is determined by the packing of effective overlapping particles as those used in the simulations by durian ,  @xcite even though the relative proportions of a fluid mixture is kept constant ; in a colloidal suspension , the existence of electrostatic screening contributes to an effective size of the particle that depends on the packing structure . in fact , the assumption made in this work is rather general , and many situation can be expected to lead to a dynamically driven effective free - volume ; this may explain the ubiquity of glassy behavior . the theoretical approach to glassy materials can be summarized as : how ergodicity breaks ? or what ergodicity breaking is sufficient to account for slow modes of relaxation . the introduction of a whole distribution of timescales is a very general way to account for non - ergodicity . here i assume that a much weaker ergodicity breaking is sufficient , and this is directly supported by the observation that fluctuation - dissipation theorem holds in a weak ( but , not so weak ) form . this indicates that the underlying distributions of state variables should not ( in fact ) be too weird , and that an effective width accounts for its essential features . however , since the system is out - of - equilibrium , this width must , at least , become a dynamical quantity . the equation proposed appear as a very minimal first order closure of the dynamics of an underlying distribution of voids . the set of curves displayed in figure  [ fig : relax3 ] is very reminiscent of the temperature dependency of relaxation spectra . @xcite obviously , the parameters of the theory should depend on temperature and other thermodynamical parameters , but here i purposely do not intend to address this question . there are ( at least ) two reasons for that : firstly , because glassy behavior is observed in many systems , in particular in suspensions or in granular materials , where thermodynamical temperature is not a relevant parameter . secondly , because it seems to me more interesting to extract all possible conclusions from the simple equations proposed without trying to incorporate them into a more specific setup . i do think that an effective dynamically driven free - volume should emerge from microscopic dynamics under minimal assumptions . therefore , it seems to me important to try to get some fundamental understanding of the consequences of such assumptions for a fixed set of parameters , without complicating the discussion . the dependency of the parameters of this theory on thermodynamic quantities is definitely an important issue to address in the future , but it first requires to have some intuition about the outcome of those equations . this work is intended to provide a grid of possible behavior displayed by simple equations , and i hope it can be helpful for further developments .    finally , one major interest of the current work is , in my opinion , that is gives a unified picture where different mechanisms for jamming coexist . this coexistence of jamming with non - linear rheology has been observed experimentally , but could not be captured by earlier model of glassy rheology . moreover , it opens the way towards a description of glassy materials where thermodynamics ( that determine the parameters of constitutive equations ) and dynamics coexists . this work has primarily benefited from the support and encouragements of jean carlson and jim langer who received me among their people in the busy heights of ucsb ; i acknowledge the opportunity they gave me to pursue my ideas , and their continuous interest for my work . i have enjoyed useful discussions with them , and also with others : ralph archuleta , pascal favreau , anthony foglia , delphine gourdon , jacob israelachvili , daniel lavalle , jean - christophe nave . i also would like to thank daniel bonn for giving me early versions of his papers , which were helpful in the writing of this manuscript . this work was supported by the w. m. keck foundation , and the nsf grant no . dmr-9813752 , and epri / dod through the program on interactive complex networks . the system  ( [ eqn : sigma:0 ] ) ,  ( [ eqn : stzdil : iso:1 ] ) ,  ( [ eqn : stzdil : iso:2 ] ) reduces to two equations for @xmath2 and @xmath20 , @xmath274\sigma\\ \dot\chi & = & -e_1 \exp\left[-{\kappa\over\chi}\right ] + \alpha\,e_0\,\exp\left[-{1\over\chi}\right]\,\sigma^2\end{aligned}\ ] ] and the jacobian of this dynamical system reads : @xmath275 ^ 2\cr 2\,e_0\,\alpha\,\sigma\,\phi &   \left(-\kappa\,e_1\phi^\kappa + e_0\,\alpha \,{\sigma } ^2\,\phi\right)\log[\phi]^2\cr } \right ) \quad.\ ] ]    in the steady state , the variable @xmath2 verifies , @xmath276 and @xmath277 takes on the value , @xmath278 for all @xmath279 , @xmath280 is a strictly increasing function of @xmath184 . in terms of @xmath281 , the eigenvalues @xmath282 of the jacobian verify : @xmath283 ^ 2 \right ) & + & \\ { \mu\,e_0\,e_1\,(\kappa+1)}\;\phi_{\rm s}^{\kappa+1}\,\ln\left[\phi_{\rm s}\right]^2 & = & 0 \quad,\end{aligned}\ ] ] and the hopf bifurcation is determined by a critical value of @xmath144 as a function of @xmath280 : @xmath284 ^ 2\ ] ] or , @xmath217 ^ 2 \quad.\ ] ]      in the general case , the system  ( [ eqn : sigma:0 ] ) ,  ( [ eqn : stzdil : lin:1]-[eqn : stzdil : lin:3 ] ) involves one additional variable ( @xmath64 ) , and the bifurcation analysis is somewhat more complicated . in the steady state , the relation between strain rate and free - volume can be rewritten as , @xmath285 since this expression is an strictly increasing function of @xmath277 , it is possible to carry out the analysis in terms of @xmath286 as an effective parameter . when the strain rate varies between @xmath287 and @xmath192 , @xmath277 varies on the interval @xmath288 $ ] . the values above @xmath214 are meaningless because there is no solution to the equations of motion .    the hopf criterion can be written in the form , @xmath289 with : @xmath290 the quantity @xmath291 has the opposite sign of the real part of the complex eigenvalues of the jacobian : the steady state is stable , iff @xmath292 .    for @xmath115 , and for any @xmath293 $ ] , @xmath294 is positive : the two values of @xmath144 , solutions of @xmath295 have the same sign . this sign can change only at a point where @xmath294 vanishes , that is @xmath296 or @xmath214 . the sign shared by these solutions on the interval @xmath293 $ ] can then be evaluated by taking the limit value of the coefficient @xmath297 close to 0 or 1 : it is positive . therefore , the equation @xmath295 has no real positive solution for @xmath115 : there is no hopf instability . it can be checked that there is no other instability , therefore , for @xmath115 , the system is stable . i am now considering the case when @xmath117 . to get a first idea of the phase diagram , in the parameter space @xmath298 , the quantity @xmath291 can be evaluated on the line @xmath299 : @xmath300 . on this line , the positiveness of @xmath291 is then equivalent to , @xmath301 which requires that : @xmath302 \quad.\ ] ]    for @xmath303 the coefficient @xmath294 is negative , hence , equation @xmath295 has a single positive solution . on the interval @xmath304 $ ] this equation admits two solutions which are positive because @xmath297 is negative close to @xmath296 : the solution @xmath305 which vanishes at the point @xmath306 and the solution @xmath307 which exists on the whole interval @xmath288 $ ] . the unstability occurs for @xmath308 $ ] . at small shear rates , or small @xmath277 , the criterion @xmath291 is dominated by : @xmath309 whence , @xmath310 the solution @xmath307 can then be identified to : @xmath311 close to @xmath296 . the same expression was found in the isotropic limit , but here , the relation between @xmath277 and the shear rate , @xmath184 is different . the solution @xmath305 can be expanded as , @xmath312
Abstract:

Answer: typical properties of glassy materials are shown to be captured by a mean - field free - volume theory . relaxation processes are supposed to be free - volume activated , and different entropy barriers are associated with density relaxation and shear motion . free - volume time logarithmic relaxation , kohlrausch - williams - watts , and power law viscosity result from the non - linear dynamics of spatially averaged quantities . the exponents associated with these phenomena are related to a single parameter of the theory . the theory also accounts for coexistence of jamming transitions and non - linear rheology .