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: [ [ stress - fluctuation - formalism . ] ] stress fluctuation formalism . + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    among the fundamental properties of any solid material are its elastic constants @xcite . they describe the macroscopic response to weak external deformations and encode information about the potential landscape of the system . in their seminal work @xcite squire , holt and hoover derived an expression for the isothermal elastic constants extending the classical born theory @xcite to canonical ensembles of imposed particle number @xmath11 , volume @xmath12 and _ finite _ ( mean ) temperature @xmath5 . they found a correction term @xmath13 to the born expression for the elastic constants which involves the mean - square fluctuations of the stress @xcite @xmath14 with @xmath15 being the inverse temperature , @xmath16 boltzmann s constant , @xmath17 the instantaneous value of a component of the stress tensor @xmath18 and greek letters denoting the spatial coordinates in @xmath19 dimension . perhaps due to the fact that the demonstration of the  stress fluctuation formalism " in ref . @xcite assumed explicitly a well - defined reference position and a displacement field for the particles , it has not been appreciated that this approach is consistent ( at least if the initial pressure is properly taken into account ) with the well - known pressure fluctuation formula for the compression modulus @xmath20 of isotropic fluids @xcite which has been derived several decades earlier , presumably for the first time in the late 1940s by rowlinson @xcite . in any case , the stress fluctuation formalism provides a convenient route to calculating the elastic properties in computer simulations which has been widely used in the past @xcite . it has been generalized to systems with nonzero initial stress @xcite , hard - sphere interactions @xcite and arbitrary continuous potentials @xcite , or to the calculation of local mechanical properties @xcite . [ [ low - temperature - limit - and - non - affine - displacements . ] ] low - temperature limit and non - affine displacements . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in particular from ref . @xcite it has become clear that the stress fluctuation correction to the leading born term does not necessarily vanish in the zero - temperature limit . this is due to the fact that when a system is subjected to a homogeneous deformation , the ensuing particle displacements need not follow the macroscopic strain affinely @xcite . the stress fluctuation term quantifies the extent of these non - affine displacements , whereas the born term reflects the affine part of the particle displacements . how important the non - affine motions are , depends on the system under consideration @xcite . while the elastic properties of crystals with one atom per unit cell are given by the born term only , stress fluctuations are significant for crystals with more complex unit cells @xcite . they become particularly pronounced for polymer like soft materials @xcite and amorphous solids @xcite . this is revealed by recent simulation studies of various glass formers @xcite , like lennard - jones ( lj ) mixtures @xcite , polymer glasses @xcite or silica melts @xcite . many of these simulation studies concentrate on the mechanical properties deep in the glassy state , exploring for instance correlations between the non - affine displacement field and vibrational anomalies of the glass (  boson peak " ) @xcite , the mechanical heterogeneity at the nanoscale @xcite or the onset of molecular plasticity in the regime where the macroscopic deformation is still elastic @xcite . [ [ shear - stress - and - shear - modulus . ] ] shear stress and shear modulus . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    surprisingly , the temperature dependence of the elastic constants of glassy materials does not appear to have been investigated extensively using the stress fluctuation formalism @xcite . following the pioneering work by barrat _ @xcite , we present in this work such a study where we focus on the behavior of the shear modulus @xmath0 of two isotropic colloidal glass - formers in @xmath1 and @xmath2 dimensions sampled , respectively , by means of molecular dynamics ( md ) and monte carlo ( mc ) simulations @xcite . since we are interested in isotropic systems one can avoid the inconvenient tensorial notation of the full stress fluctuation formalism as shown in fig . [ fig_sketch ] . focusing on the shear in the @xmath21-plane we use @xmath22 for the mean shear stress and @xmath23 for its instantaneous value . as an important contribution to the shear modulus we shall monitor the shear stress fluctuations @xmath24 as a function of temperature @xmath5 . quite generally , the shear modulus @xmath0 may then be computed according to the stress fluctuation formula @xcite @xmath25 where @xmath26 stands for the so - called  affine shear elasticity " corresponding to an assumed _ affine _ response to an external shear strain @xmath3 . for pairwise interaction potentials it can be shown that @xmath27 i.e. @xmath26 comprises the so - called  born - lam coefficient " @xmath28 , a simple moment of the first and the second derivative of the pair potential with respect to the distance of the interacting particles @xcite , and the excess contribution @xmath29 to the total pressure @xmath6 . these notations will be properly defined in sec . [ theo_stressfluctu ] below where a simple derivation will be given . since eq . ( [ eq_gmodtt ] ) is not the only operational definition of the shear modulus discussed in this work , an index has been added reminding that the shear stress @xmath30 is used for the determination of @xmath0 . [ [ physical - motivation - and - systems - of - interest . ] ] physical motivation and systems of interest . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    if an infinitesimal shear stress @xmath4 is applied at time @xmath31 , an isotropic solid exhibits , quite generally , an elastic response quantified in strength by a finite static ( @xmath32-independent ) shear modulus ( @xmath33 ) as shown by the top curve in fig .  [ fig_sketch](b ) , whereas a liquid has a vanishing modulus ( @xmath34 ) and flows on long time scales as shown by the bottom curve . the shear modulus @xmath0 can thus serve as an  order parameter " distinguishing liquid and solid states @xcite . upon melting , the shear modulus of crystalline solids is known to vanish discontinuously with @xmath5 . a natural question which arises is that of the behavior of @xmath35 for amorphous solids and glasses in the vicinity of the glass transition temperature @xmath9 . interestingly , qualitative different theoretical predictions have been put forward by mode - coupling theory ( mct ) @xcite and some versions of replica theory @xcite predicting a _ discontinuous jump _ at the glass transition whereas other versions of replica theory @xcite suggest a _ continuous _ transition in the static limit of large sampling times @xmath32 . ( see ref . @xcite for a topical recent discussion of the replica approach . ) the numerical characterization of @xmath35 for colloidal glass - former is thus of high current interest . since the stress fluctuation formalism can readily be added to the standard simulation codes of colloidal glass - formers @xcite at imposed particle number @xmath11 , box volume @xmath12 , box shape ( @xmath36 ) and temperature @xmath5 , called @xmath37-ensembles , this suggests the use of eq . ( [ eq_gmodtt ] ) to tackle this issue . note that the latter stress fluctuation relation should also be useful for wide range of related systems , beyond the scope of the present work , including the glass transition of polymeric liquids @xcite , colloidal gels @xcite and self - assembled networks @xcite as observed experimentally for hyperbranched polymer chains with sticky end - groups @xcite or bridged networks of telechelic polymers in water - oil emulsions @xcite . [ [ validity - through - a - solid - liquid - transition . ] ] validity through a solid - liquid transition . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    this begs the delicate question of whether the general stress fluctuation formalism and especially eq . ( [ eq_gmodtt ] ) for the shear modulus only holds for elastic solids with a well - defined displacement field , as suggested in the early literature @xcite , or if the approach may actually be also applied _ through _ a solid - liquid transition up to the liquid state @xcite with @xmath38 and thus @xmath39 as it should for a liquid . by reworking the theoretical arguments leading to eq . ( [ eq_gmodtt ] ) it will be seen that this is in principal the case , at least if standard thermodynamics can be assumed to apply . the latter point is indeed not self - evident for the colloidal glasses we are interested in where _ ( i ) _ some degrees of freedom are frozen on the time scale probed and _ ( ii ) _ the measured response to an imposed infinitesimal shear @xmath4 may depend on the sampling time @xmath32 , especially for temperatures @xmath5 around @xmath9 as shown in fig . [ fig_sketch](b ) . to be on the safe side , it thus appears to be appropriate to test numerically the applicability of the stress fluctuation formula , eq . ( [ eq_gmodtt ] ) , using the corresponding strain fluctuation relations which are more fundamental on experimental grounds . [ [ strain - fluctuation - relations . ] ] strain fluctuation relations . + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    panel ( a ) of fig . [ fig_sketch ] shows a schematic setup for probing experimentally the shear modulus @xmath0 of isotropic systems confined between two rigid walls . we assume that the bottom wall is fixed and the particle number @xmath11 , the ( mean ) normal pressure @xmath6 and the ( mean ) temperature @xmath5 are kept constant . in linear response either a shear strain @xmath40 or a ( mean ) shear stress @xmath4 may be imposed . the thin solid line indicates the original unstrained reference system with @xmath36 and @xmath41 , the dashed line the deformed state . from both the mechanical and the thermodynamical point of view , the shear modulus is defined as @xmath42 where the derivative may also be taken at @xmath41 . in practice , one obtains @xmath0 by fitting a linear slope to the @xmath43-data at vanishing strain . obviously , the noise level is normally large for small strains . from the numerical point of view it would thus be more appropriate to work at imposed mean stress @xmath41 , to let the strain freely fluctuate and to sample the _ instantaneous _ values @xmath44 . the shear modulus @xmath0 is then simply obtained by linear regression @xmath45 where the index indicates that both the instantaneous values @xmath46 and @xmath30 are used . ( note that imposing the strain @xmath3 also fixes the instantaneous strain @xmath46 which thus can not fluctuate . ) the associated dimensionless correlation coefficient @xmath47 allows to determine the quality of the fit . we do not know whether eq . ( [ eq_gmodgt ] ) and eq . ( [ eq_cmodgt ] ) have actually been used in a real experiment , however since no difficulty arises to work in a computer simulation at constant @xmath4 and to record @xmath44 , we use this linear regression as our key operational definition of @xmath0 which is used to judge the correctness of the stress fluctuation relation , eq .  ( [ eq_gmodtt ] ) , computed at constant @xmath3 . to avoid additional physics at the walls , we use of course periodic boundary conditions @xcite ( the central image being sketched ) as shown by the panels on the right - hand side of fig . [ fig_sketch ] . as discussed in sec . [ theo_legendre ] , it is possible to determine @xmath48 at constant volume ( @xmath49-ensemble ) as shown by panel ( c ) or constant normal pressure ( @xmath50-ensemble ) as shown by panel ( e ) . the latter ensemble has the advantage that one can easily impose in addition the same normal pressure @xmath6 for all temperatures @xmath5 as it is justified for experimental reasons . if thermodynamics holds , it is known ( sec . [ theo_strainfluctu ] ) that @xmath51 should hold @xcite . using eq . ( [ eq_gmodgt ] ) this implies a second operational definition @xmath52 which has the advantage that only the instantaneous shear strain @xmath46 has to be recorded . the main technical point put forward in this paper is that all operational definitions yield similar results @xmath53 even for temperatures where a strong dependence on the sampling time @xmath32 is observed . in other words , eq . ( [ eq_gmods ] ) states that this time dependence applies similarly ( albeit perhaps not exactly ) to the various moments computed for the three operational definitions of @xmath0 in different ensembles . [ [ outline . ] ] outline . + + + + + + + +    the paper is organized as follows : several theoretical issues are regrouped in sec . [ sec_theo ] . we begin the discussion in sec . [ theo_strainfluctu ] by reminding a few useful general thermodynamic relations . a simple derivation of the stress fluctuation formula , eq . ( [ eq_gmodtt ] ) , for the shear modulus @xmath54 in constant-@xmath3 ensembles is given in sec . [ theo_stressfluctu ] . ( the analogous stress fluctuation formula for the compression modulus @xmath20 in constant-@xmath12 ensembles @xcite is rederived in appendix  [ app_kstressfluctu ] . ) shear stress fluctuations @xmath55 in ensembles with different boundary conditions are discussed in sec . [ theo_legendre ] . for constant-@xmath4 ensembles @xmath56 will be shown to reduce to the affine shear elasticity @xmath26 , i.e. the functional @xmath54 must vanish for all temperatures . as pointed out recently @xcite , the truncation of the interaction potentials used in computational studies implies an impulsive correction for the affine shear elasticity @xmath26 . this is summarized in sec . [ theo_trunc ] . the numerical models used in this study are presented in sec . [ sec_algo ] . we turn then in sec . [ sec_res ] to our main numerical results . starting with the high-@xmath5 liquid limit , sec . [ res_thigh ] shows that all our operational definitions correctly yield a vanishing shear modulus . moving to the opposite low-@xmath5 limit we present in sec . [ res_dynmat ] the elastic response for a topologically fixed network of harmonic springs constructed from the  dynamical matrix " of our 2d model glass system @xcite . as shown in sec . [ res_tlow ] , the same behavior is found qualitatively for low-@xmath5 glassy bead systems . the temperature dependence of stress fluctuations and shear moduli in different ensembles is then discussed in sec . [ res_tvari ] . as stated above , one central goal of this work is to determine the behavior of the shear modulus around @xmath9 . we conclude the paper in sec . [ sec_conc ] . thermodynamic relations and numerical findings related to the compression modulus @xmath20 are regrouped in appendix  [ app_k ] . we assume in the following that standard thermodynamics and thermostatistics @xcite can be applied and remind a few useful relations . let us denote an extensive variable by @xmath57 and its instantaneous value by @xmath58 , the conjugated intensive variable by @xmath59 and its instantaneous value by @xmath60 @xcite . with @xmath61 being the free energy at temperature @xmath5 and imposed @xmath57 , the intensive variable @xmath59 and the associated modulus @xmath62 are given by @xmath63 alternatively , one may consider the legendre transformation @xmath64 for which @xmath65 in our case , the extensive variable is @xmath66 , the intensive variable @xmath67 @xcite and the modulus @xmath68 @xcite . we emphasize that there is an important difference between extensive and intensive variables : if the extensive variable @xmath57 is fixed , @xmath58 can not fluctuate . if , on the other hand , the intensive variable @xmath69 is imposed , not only @xmath58 but also @xmath60 may fluctuate @xcite . in our case this means , that @xmath70 for the @xmath37-ensemble , while the shear stress fluctuations @xmath55 defined above do clearly not vanish in the @xmath49-ensemble shown in fig . [ fig_sketch](d ) . assuming @xmath59 to be imposed , it is well known that @xcite @xmath71 with @xmath15 being the inverse temperature . for our variables @xmath72 and @xmath73 , eq . ( [ eq_dxdi ] ) implies eq . ( [ eq_dgamdtau ] ) from which it is seen that eq . ( [ eq_gmodgt ] ) is consistent with eq . ( [ eq_gmodgg ] ) . the latter relation is also obtained directly using eq . ( [ eq_dxdx ] ) . as discussed in textbooks @xcite a simple average @xmath74 of some observable @xmath75 does not depend of the chosen ensemble , at least not if the system is large enough ( @xmath76 ) . a correlation function @xmath77 of two observables @xmath75 and @xmath78 may differ , however , depending on whether @xmath57 or @xmath59 are imposed . as demonstrated by lebowitz , percus and verlet in 1967 @xcite one verifies that @xmath79 if the extensive variable and at least one of the observables are identical , the left - hand side of eq . ( [ eq_dadb ] ) vanishes : if @xmath80 and @xmath81 , this yields eq .  ( [ eq_dxdi ] ) . if @xmath82 , this implies eq . ( [ eq_dxdx ] ) . more importantly , for @xmath83 one obtains using eq . ( [ eq_ftom ] ) @xmath84 if it is possible to probe the intensive variable fluctuations in both ensembles , this does in principle allow to determine the modulus @xmath62 . ( the latter equation might actually be taken as the definition of @xmath62 in cases where the thermodynamic reasoning becomes dubious . ) we emphasize that the modulus @xmath62 is an intrinsic property of the system and does _ not _ depend on the ensemble , albeit it is determined using ensemble depending correlation functions . for a thermodynamically stable system , @xmath85 , eq . ( [ eq_didi ] ) implies that @xmath86 and that both fluctuations must become similar in the limit where the modulus @xmath62 becomes small . for the shear stress fluctuations @xmath55 , eq . ( [ eq_didi ] ) corresponds to the transformation @xmath87 and one thus expects @xmath55 to be larger for @xmath49-systems than for @xmath37-systems , while in the liquid limit where @xmath88 the imposed boundary conditions should not matter . we shall verify eq . ( [ eq_muftransform ] ) numerically below . a glance at the stress fluctuation formula , eq .  ( [ eq_gmodtt ] ) , suggests the identity @xmath89 assuming eq . ( [ eq_dgamdtau ] ) and eq . ( [ eq_muamuftau ] ) , it follows for the correlation coefficient @xmath90 defined above that @xmath91 if the measured @xmath44 were indeed perfectly correlated , i.e. @xmath92 , this implies @xmath93 and thus @xmath94 . in fact , for all situations studied here we always have @xmath95 i.e. the affine shear elasticity sets only an upper bound to the shear modulus and a theory which only contains the affine strain response must _ overpredict _ @xmath0 . we show now directly that eq . ( [ eq_gmodtt ] ) and , hence , eqs . ( [ eq_muamuftau],[eq_cmodgt_mua],[eq_gbound ] ) indeed hold and this under quite general conditions .      [ [ introduction . ] ] introduction . + + + + + + + + + + + + +    as shown in fig . [ fig_sketch](a ) , we demonstrate eq . ( [ eq_gmodtt ] ) from the free energy change associated with an imposed arbitrarily small pure shear strain @xmath3 in the @xmath21-plane assuming a constant particle number @xmath11 , a constant volume @xmath12 and a constant temperature @xmath5 ( @xmath37-ensemble ) . it is supposed that the total system hamiltonian may be written as the sum of a kinetic energy and a potential energy . since for a plain shear strain at constant volume the ideal free energy contribution does not change , i.e. is irrelevant for @xmath4 and @xmath0 , we may focus on the excess free energy contribution @xmath96 due to the conservative interaction energy of the particles . the ( excess ) partition function @xmath97 of the unperturbed system at @xmath36 is the boltzmann - weighted sum over all states @xmath98 of the system which are _ accessible _ within the measurement time @xmath32 . the argument @xmath99 is a short - hand for the unstrained reference . following the derivation of the compression modulus @xmath20 by rowlinson @xcite the partition function @xmath100 of the sheared system is supposed to be the sum over the _ same _ states @xmath98 , but with a different _ metric _ @xcite corresponding to the macroscopic strain which changes the total interaction energy @xmath101 of state @xmath98 and , hence , the weight of the sheared configuration for the averages computed . this is the central hypothesis made . interestingly , it is not necessary to specify explicitly the states of the unperturbed or perturbed system , e.g. , it is irrelevant whether the particles are distinguishable or not or whether they have a well - defined reference position for defining a displacement field . [ [ shear - stress - and - modulus - for - general - potential . ] ] shear stress and modulus for general potential . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    assuming eq . ( [ eq_fexgamma ] ) to hold we compute now the mean shear stress @xmath4 and the shear modulus @xmath0 for a general interaction potential @xmath101 . we note for later convenience that @xmath102 for the derivatives of the free energy and @xmath103 for the derivatives of the excess partition function where a prime denotes the derivative of a function @xmath104 with respect to its argument @xmath105 . using eq . ( [ eq_dzdgam ] ) and taking finally the limit @xmath106 one verifies for the shear stress that @xmath107 _ defining _ the instantaneous shear stress @xcite . note that the average taken is defined as @xmath108 using the weights of the unperturbed system . the shear stress thus measures the average change of the total interaction energy @xmath101 taken at @xmath36 . the shear modulus @xmath0 is obtained using in addition eq . ( [ eq_ddlogzddgam ] ) and eq . ( [ eq_ddzddgam ] ) and taking finally the @xmath106 limit . confirming thus the stress fluctuation formula eq . ( [ eq_gmodtt ] ) stated in the introduction , this yields @xmath109 with @xmath26 being the  affine shear elasticity " @xmath110 already mentioned above and @xmath55 as defined in eq . ( [ eq_muf ] ) . the comparison of eq . ( [ eq_gstressfluctu ] ) and eq . ( [ eq_muftransform ] ) confirms eq . ( [ eq_muamuftau ] ) . [ [ comment . ] ] comment . + + + + + + + +    the affine shear elasticity @xmath26 corresponds to the change ( second derivative ) of the total energy which would be obtained if one actually strains _ affinely _ in a computer simulation a given state @xmath98 without allowing the particles to relax their position . as shown for athermal ( @xmath111 ) amorphous bodies @xcite , the positions of the particles of such a strained configuration will of course in general change slightly to minimize the interaction energy relaxing thus the elastic moduli . this is also of relevance for thermalized solids where the non - affine displacements of the particles are driven by the minimization of the free energy . it is for this reason that the shear - stress fluctuation term @xmath112 must occur in eq . ( [ eq_gstressfluctu ] ) correcting the _ overprediction _ of the shear modulus @xmath0 by the affine shear elasticity @xmath26 , eq . ( [ eq_gbound ] ) . this point has been overlooked in the early literature @xcite and only appreciated much later @xcite as discussed in barrat s review @xcite . interestingly , as has been shown by lutsko @xcite , @xmath55 and other similarly defined stress fluctuations become temperature independent in the harmonic ground state approximation for @xmath111 . probing the stress fluctuations in a low - temperature simulation allows thus to determine the elastic moduli of athermal solids . please note that up to now we have not used explicitly the coordinate transformations ( metric change ) associated with an affine shear strain . this is needed to obtain operational definitions for the instantaneous shear stress @xmath30 ( and thus for @xmath4 and @xmath55 ) and the affine shear elasticity . [ [ coordinate - transformation . ] ] coordinate transformation . + + + + + + + + + + + + + + + + + + + + + + + + + +    as shown in fig . [ fig_sketch](a ) , we assume a pure shear strain which transforms the @xmath105-coordinate of a particle position or the distance between two particles as @xmath113 leaving all other coordinates unchanged @xcite . the squared distance @xmath114 between two particles thus transforms as @xmath115 where we need to keep the @xmath116-term for calculating the shear modulus below . we note for later reference that this implies @xmath117 for , respectively , the first and the second derivative with respect to @xmath3 . let us now consider an arbitrary function @xmath118 . using eq . ( [ eq_strain_rone ] ) it is readily seen that its first derivative with respect to @xmath3 may be written as @xmath119 in the second step we have introduced the components @xmath120 and @xmath121 of the normalized distance vector between both particles . more generally , we denote by @xmath122 the @xmath123-component of a normalized vector with @xmath124 stating only the small-@xmath3 limit for the second derivative of @xmath118 with respect to @xmath3 we note finally that @xmath125 where we have dropped the argument @xmath99 for the distance @xmath126 of the unperturbed reference system .    [ [ pair - interaction - potentials . ] ] pair interaction potentials . + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we assume now that the interaction energy @xmath101 of a configuration @xmath98 is given by a pair interaction potential @xmath127 between the particles @xcite with @xmath128 being the distance between two particles , i.e. @xmath129 where the index @xmath130 labels the interaction between the particles @xmath131 and @xmath132 with @xmath133 . using the general result eq . ( [ eq_tauhat_us ] ) for the shear stress it follows using eq . ( [ eq_strain_fronea ] ) that @xmath134 where we have dropped the argument @xmath99 in the @xmath106 limit . we have thus rederived for the shear stress ( @xmath135 ) the well - known kirkwood expression for the general excess stress tensor @xcite @xmath136 the affine shear elasticity @xmath26 is obtained using eq . ( [ eq_muadef ] ) and eq . ( [ eq_strain_frtwo ] ) . this yields @xmath137 with the first contribution @xmath28 being the so - called  born - lam coefficient " @xcite @xmath138 it corresponds to the first term in eq . ( [ eq_strain_frtwo ] ) . the second contribution @xmath139 stands for the normal excess stress in the @xmath140-direction , i.e. the @xmath141 component of the excess stress tensor . [ [ isotropic - systems . ] ] isotropic systems . + + + + + + + + + + + + + + + + + +    in this work we focus on isotropic materials under isotropic external loads . the normal ( excess ) stresses @xmath142 for all directions @xmath123 are thus identical , i.e. @xmath143 with @xmath144 being the excess pressure , @xmath6 the total pressure , @xmath145 the ideal gas pressure and @xmath19 the spatial dimension . this implies that the affine shear elasticity is given by @xmath146 as already stated in the introduction , eq . ( [ eq_muamubpex ] ) . we note that for isotropic @xmath19-dimensional systems it can be shown @xcite that the born - lam coefficient can be simplified as @xmath147 the @xmath19-dependent prefactor stems from the assumed isotropy of the system and the mathematical formula @xcite @xmath148 ( @xmath149 being the kronecker symbol @xcite ) for the components of a unit vector in @xmath19 dimensions pointing into arbitrary directions . by comparing the excess pressure @xmath29 , eq . ( [ eq_sigaapex ] ) , with the underlined second contribution to the born coefficient in eq . ( [ eq_mub ] ) this implies that @xmath150 it is thus inconsistent to neglect the explicit excess pressure in eq . ( [ eq_muamubpex ] ) while keeping the second term of the born - lam coefficient @xmath28 . this approximation is only justified if the excess pressure @xmath29 is negligible ( and not the total pressure @xmath6 ) . [ [ different - ensembles . ] ] different ensembles . + + + + + + + + + + + + + + + + + + + +    being simple averages neither @xmath29 , @xmath26 or @xmath28 does depend for sufficiently large systems on the chosen ensemble , i.e. irrespective of whether @xmath3 or @xmath4 is imposed one expects to obtain the same values . as already reminded in sec . [ theo_strainfluctu ] , this is different for fluctuations in general @xcite . for instance , it is obviously pointless to use the shear - strain fluctuation formulae @xmath48 or @xmath151 in a constant-@xmath3 ensemble . a more interesting result is predicted if @xmath152 is computed in the  wrong " constant-@xmath4 ensemble . since @xmath153 , eq . ( [ eq_muftransform ] ) , one actually expects to observe @xmath154 for all temperatures @xmath5 if our thermodynamic reasoning holds through the glass transition up to the liquid state . we will test numerically this non - trivial prediction in sec . [ sec_res ] . up to now we have only considered for simplicity of the presentation @xmath37- and @xmath49-ensembles at constant volume @xmath12 . since on general grounds pure deviatoric and dilatational ( volumetric ) strains are decoupled ( both strains commute ) , one expects the observable @xmath54 also to be applicable in the @xmath155-ensemble shown in fig . [ fig_sketch](d ) and the observables @xmath48 and @xmath151 to be applicable in the @xmath50-ensemble as sketched in fig . [ fig_sketch](e ) . in the latter ensemble one also expects to find @xmath156 and @xmath153 . this will also be tested below . [ [ direct - derivation - of - ensuremathmu_mathrmf - for - imposed - tau . ] ] direct derivation of @xmath55 for imposed @xmath4 . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    we demonstrate now directly eq . ( [ eq_muamuftau ] ) for the shear stress fluctuations using the general mathematical identity @xmath157 with @xmath105 being an _ unconstrained _ coordinate , @xmath158 some property , @xmath159 a  force " with respect to some  energy " @xmath160 and the average being boltzmann weighted , i.e. @xmath161 . ( it is assumed in eq . ( [ eq_intbyparts ] ) that @xmath160 decays sufficiently fast at the boundaries . ) following work by zwanzig @xcite we express first the instantaneous shear stress @xmath30 given above by the kirkwood formula , eq . ( [ eq_tauhatkirkwood ] ) , by the alternative virial representation @xcite @xmath162 the second term , sometimes called the  inner virial " , stands for a sum over all particles @xmath131 with @xmath163 being the @xmath140-coordinate of the particle and @xmath164 the @xmath105-coordinate of the force acting on the particle . this contribution indeed vanishes on average as can be seen using eq . ( [ eq_intbyparts ] ) , @xmath165 which shows that @xmath166 as it must . the stress fluctuations may thus be expressed as @xmath167 where we have used the integration by part , eq . ( [ eq_tauhatfluctu ] ) . this step requires that all particle positions are unconstrained and independent ( generalized ) coordinates . note that this is possible at imposed @xmath4 , but can not hold at fixed @xmath3 . assuming then pair interactions and using eq . ( [ eq_tauhatkirkwood ] ) one can reformulate eq . ( [ eq_stressfluctuvirial ] ) within a few lines . without further approximation this yields @xmath168 where the sum runs now over all interactions @xmath130 between particles @xmath133 . @xmath169 and @xmath170 refer to components of the distance vector of the interaction @xmath130 and @xmath171 to the @xmath105-component of the central force @xmath172 between a particle pair at a distance @xmath173 . note that in eq . ( [ eq_stressfluctdecoupled ] ) the stress fluctuations stemming from different interactions are decoupled , i.e. @xmath174 characterizes the self- or two - point contributions of directly interacting beads . since @xmath175 and , hence , @xmath176 this is identical to the affine shear elasticity @xmath26 derived at the end of sec . [ theo_stressfluctu ] . we have thus confirmed eq . ( [ eq_muamuftau ] ) . interestingly , the argument can be turned around :    1 . the stress fluctuation contribution term @xmath177 is obtained by the simple derivation given in this paragraph . 2 .   using the general legendre transformation , eq . ( [ eq_muftransform ] ) , this implies @xmath178 for the @xmath37- and the @xmath155-ensemble .    as already stressed , our thermodynamic derivation of eq . ( [ eq_gmodtt ] ) is rather general , not using in particular a well - defined displacement field for the particle positions . [ [ truncation . ] ] truncation . + + + + + + + + + + +    it is common practice in computational condensed matter physics @xcite to truncate a pair interaction potential @xmath127 , with @xmath126 being the distance between two particles @xmath131 and @xmath132 , at a conveniently chosen cutoff @xmath179 . this allows to reduce the number of interactions to be computed and energy or force calculations become @xmath180-processes . however , the truncation also introduces technical difficulties , e.g. , instabilities in the numerical solution of differential equations as investigated especially for the md method @xcite . without restricting much in practice the generality of our results , we assume below that    * the pair potential @xmath127 is short - ranged , i.e. that it decays within a few particle diameters , * it scales as @xmath181 with the reduced dimensionless distance @xmath182 where @xmath183 characterizes the length scale of the interaction @xmath130 and * the same reduced cutoff @xmath184 is set for all interactions @xmath130 .    for instance , for monodisperse particles with constant diameter @xmath185 , as for the standard lennard - jones ( lj ) potential @xcite , @xmath186 the scaling variable becomes @xmath187 and the reduced cutoff @xmath188 . the effect of introducing @xmath189 is to replace @xmath190 by the truncated potential @xmath191 with @xmath192 being the heaviside function @xcite . even if eq . ( [ eq_utrunc ] ) is taken _ by definition _ as the new system hamiltonian , it is well known that impulsive corrections at the cutoff have to be taken into account in general for the excess pressure @xmath29 and other moments of the first derivatives of the potential @xcite . this is seen by considering the derivative of the truncated potential @xmath193    [ [ shifting - of - truncated - potential . ] ] shifting of truncated potential . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    these corrections with respect to first derivatives can be avoided by considering a properly shifted potential @xcite @xmath194 since @xmath195 . with this choice _ no _ impulsive corrections arise for moments of the instantaneous shear stress @xmath30 and of other components of the excess stress tensor @xmath196 , eq . ( [ eq_sigab ] ) . specifically , if the potential is shifted , all impulsive corrections are avoided for the shear stress fluctuations @xmath55 , eq . ( [ eq_muf ] ) .    [ [ correction - to - the - born - coefficient . ] ] correction to the born coefficient . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    as shown in ref . @xcite , the standard shifting of a truncated potential is , however , insufficient in general for properties involving second ( and higher ) derivatives of the potential . this is particulary the case for the born coefficient @xmath28 , eq . ( [ eq_mub ] ) , which is required to compute the affine shear elasticity @xmath26 , eq . ( [ eq_muamubpex ] ) . the second derivative of the truncated and shifted potential being @xmath197 this implies that the born coefficient reads @xmath198 with @xmath199 being the bare born coefficient and @xmath200 the impulsive correction which needs to be taken into account . the latter correction may be readily computed numerically from the configuration ensemble using @xcite @xmath201 being a weighted radial pair distribution function which is related to the standard radial pair distribution function @xmath202 @xcite . [ [ mixtures - and - polydisperse - systems . ] ] mixtures and polydisperse systems . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    below we shall consider model systems for mixtures and polydisperse systems where @xmath183 may differ for each interaction @xmath130 . for such mixed potentials @xmath190 , @xmath203 and @xmath204 and their derivatives take in principal an explicit index @xmath130 , i.e. one should write @xmath205 , @xmath206 , @xmath207 and so on . ( this is often not stated explicitly to keep a concise notation . ) for example one might wish to consider the generalization of the monodisperse lj potential , eq . ( [ eq_lj ] ) , to a mixture or polydisperse system with @xmath208 where @xmath209 and @xmath183 are fixed for each interaction @xmath130 . in practice , each particle @xmath131 may be characterized by an energy scale @xmath210 and a  diameter " @xmath211 . the interaction parameters @xmath212 and @xmath213 are then given in terms of specified functions of these properties @xcite . the extensively studied kob - andersen ( ka ) model for binary mixtures of beads of type a and b @xcite , is a particular case of eq . ( [ eq_ljgeneralized ] ) with fixed interaction ranges @xmath214 , @xmath215 and @xmath216 and energy parameters @xmath217 , @xmath218 and @xmath219 characterizing , respectively , aa- , bb- and ab - contacts . the expression eq . ( [ eq_histolame ] ) remains valid for such explicitly @xmath130-dependent potentials . [ [ shear - modulus - and - compression - modulus . ] ] shear modulus and compression modulus . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    since @xmath220 , eq .  ( [ eq_histolame ] ) implies in turn @xmath221 with @xmath222 being the uncorrected bare shear modulus . as we shall see in sec . [ res_thigh ] , the correction @xmath200 is of importance for the _ precise _ determination of @xmath54 close to the glass transition where the shear modulus must vanish . an impulsive correction has also to be taken into account for the compression modulus @xmath223 computed from the normal pressure fluctuations in a constant volume ensemble as discussed in appendix  [ app_k ] . using the symmetry of isotropic systems , eq . ( [ eq_etabmub ] ) , one verifies @xcite @xmath224 with @xmath225 being the uncorrected compression modulus . [ [ introduction.-1 ] ] introduction . + + + + + + + + + + + + +    the computational data presented in this work have been obtained using two extremely well studied models of colloidal glass - formers which are described in detail elsewhere @xcite . we present first both models and the molecular dynamics ( md ) and monte carlo ( mc ) @xcite methods used to compute them in the different ensembles ( @xmath37 , @xmath49 , @xmath155  and @xmath50 ) compared in this study . we comment then on the quench protocol and locate the glass transition temperature @xmath9 for both models via dilatometry , as shown in fig . [ fig_quench_rho ] . this is needed as a prerequisite to our study of the elastic behavior in sec . [ sec_res ] . as sketched in fig . [ fig_sketch ] , periodic boundary conditions are applied in all spatial directions with @xmath226 for the linear dimensions of the @xmath19-dimensional simulation box . time scales are measured in units of the lj time @xmath227 @xcite for our md simulations ( @xmath228 being the particle mass , @xmath185 the reference length scale and @xmath229 the lj energy scale ) and in units of monte carlo steps ( mcs ) for our mc simulations . lj units are used throughout this work ( @xmath230 ) and boltzmann s constant @xmath16 is also set to unity . [ [ kob - andersen - model . ] ] kob - andersen model . + + + + + + + + + + + + + + + + + + +    the kob - andersen ( ka ) model @xcite for binary mixtures of lj beads in @xmath1 dimensions has been investigated by means of md simulation @xcite taking advantage of the public domain lammps implementation @xcite . we use @xmath231 beads per simulation box and molar fractions @xmath232 and @xmath233 for both types of beads @xmath234 and @xmath235 . following ref . @xcite we set @xmath236 , @xmath237 and @xmath238 for the interaction range and @xmath239 , @xmath240 and @xmath241 for the lj energy scales . only data for the usual ( reduced ) cutoff @xmath242 are presented . for the @xmath155-ensemble we use the nos - hoover barostat provided by the lammps code (  fix npt command " ) @xcite which is used to impose a constant pressure @xmath243 for all temperatures . starting with an equilibrated configuration at @xmath244 well above the glass transition temperature @xmath245 ( as confirmed in the inset of fig . [ fig_quench_rho ] ) , the system is slowly quenched in the @xmath155-ensemble . the imposed mean temperature @xmath246 varies linearly as @xmath247 with @xmath248 being the constant quench rate . after a tempering time @xmath249 we compute over a sampling time @xmath250 various properties for the isobaric system at fixed temperature such as the moduli @xmath54 represented by the filled spheres in fig . [ fig_gt ] . fixing only then the volume ( @xmath37-ensemble ) the system is again tempered ( @xmath249 ) and various properties as the ones recorded in table  [ tab_ka ] are computed using a sampling time @xmath250 . note that in the @xmath37-ensemble the equations of motion are integrated using a velocity verlet algorithm @xcite with a time step @xmath251 and systems are kept at the imposed temperature @xmath5 using a langevin thermostat @xcite with friction coefficient @xmath252 . by redoing for several temperatures around and below @xmath9 the tempering and the production runs we have checked that the presented results remain unchanged ( within numerical accuracy ) and that ageing effects are absolutely negligible . averages are taken over four independent configurations . as one may expect from the discussion in sec . [ theo_legendre ] , very similar results have been found for both ensembles for simple averages , for instance for the affine shear elasticity @xmath26 , and for the shear modulus @xmath54 , eq . ( [ eq_gmodtt ] ) . unfortunately , we are yet unable to present data for the ka model data at imposed shear stress @xmath4 ( @xmath49- , @xmath50-ensembles ) . .some properties for the ka model at normal pressure @xmath243 and shear stress @xmath41 : the mean density @xmath253 , the mean energy per volume @xmath254 , the impulsive correction term @xmath200 ( sec . [ theo_trunc ] ) , the ( corrected ) born - lam coefficient @xmath28 , the affine shear elasticity @xmath255 , the shear modulus @xmath54 obtained in the @xmath37-ensemble , the hypervirial @xmath256 and the affine dilatation elasticity @xmath257 discussed in appendix  [ app_k ] and the compression modulus @xmath258 obtained for the @xmath37-ensemble . the affine elasticities @xmath26 and @xmath259 are the natural scales for , respectively , the shear modulus @xmath0 and the compression modulus @xmath20 . [ tab_ka ] [ cols="^,^,^,^,^,^,^,^,^,^,^,^",options="header " , ]     [ [ local - mc - moves - for - the - plj - model . ] ] local mc moves for the plj model . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    for all ensembles considered here local mc moves are used ( albeit not exclusively ) where a particle is chosen randomly and a displacement @xmath260 , uniformly distributed over a disk of radius @xmath261 , from the current position @xmath262 of the particle is attempted . the corresponding energy change @xmath263 is calculated and the move is accepted using the standard metropolis criterion @xcite . the maxium displacement distance @xmath261 is chosen such that the acceptance rate @xmath234 remains reasonable , i.e. @xmath264 @xcite . as may be seen from the inset in fig . [ fig_quench_e_plj ] where the acceptance rate @xmath265 is shown as a function of temperature @xmath5 , we find that temperatures below @xmath266 are best sampled using @xmath267 , whereas @xmath268 is a good choice for the interesting temperature regime between @xmath269 and @xmath270 around the glass transition temperature @xmath271 . ( see the main panel of fig . [ fig_quench_rho ] for the determination of @xmath9 . ) various data sets are presented in sec . [ sec_res ] for a fixed value @xmath272 which allows a reasonable comparison of the sampling time dependence for temperatures around the glass transition . this value is not necessarily the best choice for tempering the system and for computing static equilibrium properties at the given temperature , especially below @xmath273 .    [ [ collective - plane - wave - mc - moves . ] ] collective plane - wave mc moves . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    local mc jumps in the @xmath37-ensemble become obviously inefficient for relaxing large scale structural properties @xcite for the lowest temperatures computed for the plj model ( @xmath274 ) as may be seen , e.g. , from the finite , albeit very small shear stresses @xmath4 which happen to occur naturally in some of the quenched configurations ( especially if the quench rate @xmath275 is too large ) and can only with difficulty be relaxed using local jumps _ alone _ @xcite . we have thus crosschecked and improved the values obtained using only local mc moves by adding global mc move attempts for all particles with longitudinal and ( more importantly ) transverse plane waves commensurate with the simulation box . the random amplitudes and phases of the waves are assumed to be uniformly distributed . note that the maximum amplitude @xmath276 of each wave is chosen inversely proportional to the length @xmath277 of the wavevector in agreement with continuum theory . this allows to keep the same acceptance rate for each global mode if @xmath278 . as one expects , the plane - wave displacement attempts becomes inappropriate for large @xmath279 where continuum elasticity breaks down @xcite and , hence , the acceptance rate @xmath280 gets too large to be efficient @xcite . note that the maximum amplitude @xmath276 has to be chosen much smaller for the longitudinal waves as for the transverse ones , i.e. since the systems are essentially incompressible , the transverse modes naturally are more effective for quenching and equilibrating the systems . these global plane - wave moves have been added for the tempering of all presented configurations below @xmath270 with @xmath281 and for production runs over sampling times @xmath282 for static properties for @xmath283 . details concerning the collective plane - wave mc moves must be given elsewhere . [ [ mc - sampling - of - different - ensembles . ] ] mc sampling of different ensembles . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in this paper we shall be concerned more with the consequences of affine displacements associated with collective mc moves changing the overall box volume @xmath12 and shear strain @xmath3 . using again a metropolis criterion for accepting a suggested change of @xmath12 or @xmath3 , these pure dilatational and pure shear strain fluctuations are used to impose a mean normal pressure @xmath284 and a mean shear stress @xmath41 . for a shear strain altering mc move one first chooses randomly a strain fluctuation @xmath285 with @xmath286 . assuming an imposed shear stress @xmath4 the suggested @xmath285 is accepted if @xmath287 with @xmath288 being a uniformly distributed random variable with @xmath289 . the energy change @xmath263 associated with the affine strain may be calculated using eq . ( [ eq_strain_r ] ) . since @xmath290 is supposed in this study , the last term in the exponential drops out . ( as noted in sec . [ theo_stressfluctu ] , the volume remains unchanged in a pure shear strain and , hence , no translational entropy contribution appears . ) the maximum shear strain step @xmath291 is fixed such that at the given temperature the acceptance rate remains reasonable @xcite . as described in more detail in ref . @xcite , one similarly chooses a relative volume change @xmath292 with @xmath12 being the current volume @xcite . the volume altering move is accepted at an imposed pressure @xmath6 if @xmath293 where @xmath263 may be computed using eq . ( [ eq_strain_r_epsilon ] ) . the logarithmic contribution corresponds to the change of the translational entropy . using eq . ( [ eq_timposed ] ) after each mcs for local mc moves ( and global plane - wave mc moves if these are added ) realizes the @xmath49-ensemble shown in fig . [ fig_sketch](c ) , using eq . ( [ eq_pimposed ] ) the @xmath155-ensemble shown in fig . [ fig_sketch](d ) and using in turn both strain fluctuations the @xmath50-ensemble in fig . [ fig_sketch](e ) . we remind that @xmath37- and @xmath155-ensembles are used to determine @xmath54 and @xmath49- and @xmath50-ensembles to obtain @xmath48 and @xmath151 . as shown in fig . [ fig_quench_rho ] , we quench the configurations starting from @xmath294 in the @xmath155- and the @xmath50-ensemble using a constant quench rate @xmath275 . note that the smallest quench rate @xmath295 for one configuration in the @xmath155-shown did require alone a run over six months using one core of an intel xeon e5410 processor . interestingly , similar results are obtained with the @xmath50-ensemble only using a rate @xmath296 . ( due to the additional computation for the shear strain fluctuations this is , however , only a factor five faster . ) the configurations created by the @xmath155-quenches are used ( after tempering ) for the sampling in the @xmath37- and @xmath155-ensembles , the configurations obtained by the @xmath50-quenches for the sampling in the @xmath49- and @xmath50-ensembles . only two independent configurations have been sampled following the full protocol for each temperature and each of the four ensembles . instead of increasing further the number of configurations ( which will be done in the future ) we have focused in the present preliminary study on long tempering times and production runs over @xmath282 mcs which have been redone several times for a few temperatures ( summing up to total runs of @xmath297 mcs for @xmath269 and @xmath298 ) to check for equilibration problems and to verify that ageing effects can be ignored . we firmly believe that the breaking of the time translational invariance is not a relevant issue for the strong sampling time dependence reported below for the shear moduli .    [ [ algo_sample ] ] glass transition temperature . + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the inset of fig . [ fig_quench_rho ] presents for the ka model how the glass transition temperature @xmath9 may be obtained using standard dilatometry where the rescaled density @xmath299 is traced as a function of temperature @xmath5 as in ref . @xcite for a similar polymer glass problem . this reveals a linear low-@xmath5 ( solid line ) and a linear high-@xmath5 regime . the intercept of both lines confirms the well - known value @xmath300 from the literature @xcite . the main panel presents the unscaled number density @xmath253 for plj systems at constant normal pressure @xmath284 as a function of temperature @xmath5 for different quench rates @xmath275 in mcs . only local mc moves with a fixed maximum particle displacement @xmath301 for all temperatures @xmath5 are used for the examples given . the crosses refer to a quench in the @xmath50-ensemble where the box shape is allowed to fluctuate at @xmath41 . all other data refer to the @xmath155-ensemble with @xmath36 . the solid line and the dashed line are linear fits to the best density estimates indicated by the stars obtained for , respectively , low and high temperatures . by matching both lines one determines @xmath302 . a similar value is given if @xmath299 is plotted as a function of @xmath5 . the main panel of fig . [ fig_quench_e_plj ] shows the interaction energy per particle @xmath303 of the plj model as a function of temperature for the same quench protocols as in fig . [ fig_quench_rho ] . the solid line and the dashed line are linear fits to the energy for , respectively , low and high temperatures . this confirms @xmath302 for the plj model @xcite . we begin our discussion by focusing on the high temperature liquid limit where all reasonable operational definitions of the shear modulus @xmath0 must vanish . the results presented in fig . [ fig_gt_thigh ] have been obtained for the plj model at mean pressure @xmath284 , mean shear stress @xmath304 and mean temperature @xmath305 , i.e. far above the glass transition temperature @xmath302 . only local mc monomer displacements with a maximum jump displacement @xmath272 are reported here . instantaneous properties relevant for the moments are written down every @xmath306 mcs and averaged using standard gliding averages @xcite , i.e. we compute mean values and fluctuations for a given time interval @xmath307 $ ] and average over all possible intervals of length @xmath32 . the horizontal axis indicates the latter interval length @xmath32 in units of mcs . the vertical axis is made dimensionless by rescaling the moduli with the ( corrected ) affine shear elasticity @xmath308 from table  [ tab_plj ] . ( being a simple average , @xmath26 is found to be identical for all computed ensembles . ) in agreement with eq . ( [ eq_gbound ] ) , the ratio @xmath309 for all operational definitions drops rapidly below unity . the dashed line indicates the asymptotic power - law slope @xmath310 @xcite . this can be understood by noting that at imposed @xmath4 the shear strain freely diffuses in the liquid limit , i.e. @xmath311 , which implies that @xmath312 . for both @xmath49- and @xmath50-ensembles we compare our key definition @xmath48 , eq . ( [ eq_gmodgt ] ) , to the observable @xmath151 , eq . ( [ eq_gmodgg ] ) . in both cases we confirm @xmath313 for larger times as suggested by the thermodynamic relation eq . ( [ eq_dgamdtau ] ) which has been directly checked @xcite . we have also computed the dimensionless correlation coefficient @xmath314 , eq . ( [ eq_cmodgt ] ) , which characterize the correlation of the measured instantaneous shear strains @xmath46 and shear stresses @xmath30 . as indicated for the @xmath50-ensemble ( stars ) , @xmath90 vanishes rapidly with @xmath32 , i.e. @xmath46 and @xmath30 are decorrelated . the dash - dotted line indicates the power - law exponent @xmath315 expected from eq . ( [ eq_cmodgt_mua ] ) and @xmath316 . note that by plotting @xmath314 _ vs. _ @xmath317 , eq . ( [ eq_cmodgt_mua ] ) has been directly verified for both @xmath49- and @xmath50-ensembles ( not shown ) .    the stress fluctuation formula , eq . ( [ eq_gmodtt ] ) , is represented in fig . [ fig_gt_thigh ] by squares and diamonds for , respectively , @xmath37- and @xmath155-ensembles . the bare shear modulus @xmath222 _ without _ the impulsive correction @xmath200 for the born coefficient @xmath28 discussed in sec . [ theo_trunc ] is given by small filled symbols . the impulsive correction @xmath318 ( table  [ tab_plj ] ) computed independently from the weighted radial pair correlation function , eq .  ( [ eq_histolame ] ) , is shown by a horizontal line . the corrected moduli @xmath319 vanish indeed as expected . see ref . @xcite for a systematic numerical characterization of @xmath200 as a function of the reduced cutoff @xmath189 . as may also be seen there , a truncation correction is also necessary for the ka model . the @xmath200 for different temperatures may be found in table  [ tab_ka ] for the ka model and in table  [ tab_plj ] for the plj model . we assume from now on that the impulsive corrections are properly taken into account without stating this technical point explicitly . [ [ introduction.-2 ] ] introduction . + + + + + + + + + + + + +    we turn now to the opposite low-@xmath5 limit focusing specifically on the plj model at @xmath320 . obviously , such deeply quenched colloidal glasses must behave as amorphous solids with a constant shear modulus for large times ( albeit not infinite times ) as sketched in fig . [ fig_sketch](b ) by the top curve . before presenting in sec . [ res_tlow ] the elastic properties of these glasses , we discuss first conceptually simpler substitute systems formed by _ permanent _ spring networks . we do this for illustration purposes since questions related to ergodicity and ageing are by construction irrelevant and the thermodynamical relations reminded in sec . [ sec_theo ] should hold rigorously .    [ [ network - hamiltonian . ] ] network hamiltonian . + + + + + + + + + + + + + + + + + + + +    assuming ideal harmonic springs the interaction energy reads @xmath321 with @xmath322 being the spring constant and @xmath323 the reference length of spring @xmath130 . the sum runs over all springs @xmath130 between topologically connected vertices @xmath131 and @xmath132 of the network . we assume here a strongly and homogeneously connected network where every vertex @xmath131 is in contact with many neighbors @xmath132 . such a network may be constructed from a plj configuration at @xmath320 keeping the particle positions for the positions of the vertices and replacing each lj interaction @xmath130 by a spring of spring constant @xmath322 and reference length @xmath323 . the simplest choice we have investigated is to set @xmath324 and @xmath325 which corresponds to a well - defined ground state of energy @xmath326 and pressure @xmath327 at @xmath328 ( not shown ) . we discuss here instead a slightly more realistic choice of @xmath322 and @xmath323 which corresponds to the _ same _  dynamical matrix " as the plj configuration at @xmath329 and @xmath284 , i.e. the same second derivative of the total interaction potential with respect to the particle positions @xcite . this choice imposes the setting @xmath330 with @xmath331 being the interaction potential for reduced distances @xmath332 . ( impulsive corrections at @xmath333 are neglected here . ) [ [ sampling . ] ] sampling . + + + + + + + + +    taking advantage of the fact that all interacting vertices are known , these networks can be computed using global mc moves with longitudinal and transverse planar waves as described in sec .  [ sec_algo ] . no local mc jumps for individual vertices have been added . note that although some @xmath322 and @xmath323 may even be negative , this does not affect the global or local stability of the network . presumably due to the fact that our spring constants are rather large , we have not observed either any buckling instability if volume fluctuations are allowed . for networks at the same pressure @xmath6 , shear stress @xmath4 and temperature @xmath5 as the original plj configuration , we obtain , not surprisingly , the same affine shear elasticity @xmath334 as for the reference . the same applies for other simple averages . [ [ shear - modulus . ] ] shear modulus . + + + + + + + + + + + + + +    figure  [ fig_gt_dynmat ] presents various measurements of the shear modulus @xmath0 for different ensembles as a function of sampling time @xmath32 . as in sec . [ res_thigh ] , we use gliding averages and the vertical axis has been rescaled by the affine shear elasticity @xmath26 . as can be seen , we obtain in the long - time limit @xmath335 ( bold line ) irrespective of whether the shear modulus is determined from @xmath48 or @xmath151 in the @xmath49- or @xmath50-ensembles or using @xmath54 in the @xmath37- or @xmath155-ensembles at fixed shear strain @xmath3 . note that even in these cases the shear modulus decreases first with time albeit the network is perfectly at thermodynamic equilibrium and no ageing occurs by construction . as shown by the dash - dotted line , @xmath336 for short times , i.e. in this limit the response is affine . as indicated by the stars , we have also computed the correlation coefficient @xmath314 as a function of time . from @xmath92 for small @xmath32 , it decreases somewhat becoming @xmath337 for large times in agreement with eq . ( [ eq_cmodgt_mua ] ) . if @xmath54 is computed in the  wrong " @xmath49- or @xmath50-ensembles , it is seen to vanish rapidly with time as shown by the filled squares . this decay is of course expected from eq . ( [ eq_muamuftau ] ) as discussed in sec . [ theo_stressfluctu ] and sec . [ theo_legendre ] . [ [ shear - stress - fluctuations . ] ] shear - stress fluctuations . + + + + + + + + + + + + + + + + + + + + + + + + + +    to emphasize the latter point we show in fig . [ fig_muft_dynmat ] the stress fluctuations @xmath338 for different ensembles . this is also done to have a closer look at the time dependence of various contributions to @xmath339 . as before , the vertical axis has been rescaled with @xmath340 . we also indicate the rescaled values of @xmath341 ( thin top line ) and @xmath342 ( dash - dotted line ) computed by a gliding average over time intervals @xmath307 $ ] . being simple averages both properties do not depend on the ensemble probed ( not shown ) and become virtually immediately @xmath32-independent as can be seen from the indicated horizontal lines . this is quite different for the different shear stress fluctuations @xmath55 presented which increase monotonously from @xmath343 for short times to their plateau value for long times . in agreement with eq . ( [ eq_muamuftau ] ) we find that @xmath344 for @xmath49- and @xmath50-ensembles and as predicted by eq . ( [ eq_muftransform ] ) we confirm @xmath345 for @xmath37- and @xmath155-ensembles . the asymptotic limit is reached after about @xmath346 mcs in the first case for imposed @xmath4 , but only after @xmath347 for imposed @xmath3 . this is due to the fact that @xmath174 corresponds to the fluctuations of the self - contribution of individual springs as discussed in sec . [ theo_legendre ] , while @xmath348 also contains contributions from stress correlations of different springs . it should be rewarding to analyze in the future finite system - size effects and the time dependence of three - dimensional ( 3d ) systems using the dynamical matrix associated to the ka model in the @xmath329 limit . we return now to the colloidal glasses where the topology of the interactions is , of course , _ not _ permanently fixed ( at least not at finite temperatures ) and the shear stresses may thus fluctuate more strongly . as an example , we present in fig . [ fig_gt_tlow ] the shear modulus @xmath349 as a function of sampling time @xmath32 for plj beads at mean pressure @xmath284 , mean shear stress @xmath350 and mean temperature @xmath320 , i.e. the same conditions as in the previous subsection . all simple averages , such as the affine shear elasticity @xmath351 , are the same for all ensembles albeit different quench protocols have been used . for all examples given we use local mc jumps with @xmath301 which allows to keep a small , but still reasonable acceptance rate @xmath352 . we have also performed runs with additional global mc moves using longitudinal and transverse plane waves . this yields similar data which is shifted horizontally to the left ( not shown ) . as one would expect , the shear moduli presented in fig . [ fig_gt_tlow ] are similar to the ones shown in fig . [ fig_gt_dynmat ] for the permanent spring network . the filled triangles show @xmath54 computed in the  wrong " @xmath50-ensemble . as expected from eq . ( [ eq_muftransform ] ) these values vanish rapidly with time . if on the other hand @xmath48 , @xmath151 and @xmath54 are computed in their natural ensembles , all these measures of @xmath0 are similar and approach , as one expects for such a low temperature , the same finite asymptotic value @xmath353 indicated by the bold horizontal line . the stars indicate the squared correlation coefficient @xmath354 which is seen to collapse perfectly on the rescaled modulus @xmath317 which shows that eq . ( [ eq_cmodgt_mua ] ) even holds for small times before the thermodynamic plateau has been reached . similar reduced shear moduli @xmath355 have been also obtained for the ka model in the low-@xmath5 limit as may be seen from table  [ tab_ka ] . these results confirm that about half of the affine shear elasticity is released by the shear stress fluctuations in agreement with refs . interestingly , this value is slightly smaller as the corresponding value for the permanent dynamical matrix model presented above . since some , albeit very small , rearrangement of the interactions must occur for the plj particles at finite temperature , this is a reasonable finding . the rescaled shear stress fluctuations @xmath356 for the same ensembles are presented in fig . [ fig_muft_tlow ] . while the two - point correlations @xmath357 ( thin horizontal line ) and @xmath358 ( dash - dotted horizontal line ) do not dependent on time , one observes again that @xmath356 increases monotonously from zero to its long - time asymptotic plateau value . in agreement with eq . ( [ eq_muamuftau ] ) we find @xmath344 for the @xmath49  and the @xmath50  ensembles , while @xmath345 ( dashed horizontal line ) for the @xmath37  and the @xmath155  ensembles . ( similar behavior has been obtained for the ka model in the @xmath37-ensemble for temperatures @xmath274 . ) as may be better seen for the data set obtained for @xmath359 using the @xmath37-ensemble ( small filled squares ) where @xmath360 , the stress fluctuations @xmath361 do not rigorously become constant , but increase extremely slowly on the logarithmic time scales presented . there is , hence , even at low temperatures always some sampling time dependence if the interactions are not permanently fixed as for the topologically fixed networks discussed above . however , this effect becomes only sizeable above @xmath269 for the ka model and above @xmath266 for the plj model . we shall now attempt to characterize this temperature dependence . [ [ stress - fluctuations . ] ] stress fluctuations . + + + + + + + + + + + + + + + + + + + +    we turn now to the characterization of the temperature dependence of the shear modulus @xmath0 and various related properties . we focus first on the best values obtained for the longest simulation runs available . figure  [ fig_muamuf_t ] presents the stress fluctuations obtained for the ka model using the @xmath37-ensemble ( open spheres ) and for the plj model using @xmath37- and the @xmath50-ensembles . the affine shear elasticity @xmath26 obtained for each system is represented by small filled symbols . as shown for both models by a bold solid line for the low-@xmath5 and by a dashed line for the high-@xmath5 regime , @xmath8 decays roughly linearly with @xmath5 . interestingly , the slopes match precisely at @xmath362 for both models . as may be seen from the large triangles for the @xmath50-ensemble , we confirm for all temperatures that @xmath363 as expected from eq . ( [ eq_muamuftau ] ) . for systems without box shape fluctuations , it is seen that @xmath364 becomes _ non - monotonous _ with a clear maximum at @xmath365 . in the liquid regime for @xmath366 where the boundary conditions do not matter , we obtain @xmath367 , i.e. all shear stress fluctuations decay with temperature . below @xmath9 the boundary constraint ( @xmath36 ) becomes more and more relevant with decreasing @xmath5 which reduces the shear stress fluctuations . according to eq . ( [ eq_muftransform ] ) , @xmath364 must thus decrease with increasing @xmath0 as the system is further cooled down . that @xmath364 must necessarily be non - monotonous with a maximum at @xmath9 is one of the central results of the presented work . [ [ shear - modulus - as - a - function - of - temperature . ] ] shear modulus as a function of temperature . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the temperature dependence of the ( unscaled ) shear moduli @xmath35 is presented in fig . [ fig_gt ] . data obtained by md simulation of the ka model ( @xmath368 , @xmath300 ) are indicated for the @xmath37-ensemble ( open spheres ) and the @xmath155-ensemble ( filled spheres ) . all other results refer to the best values obtained by mc simulations of the plj model ( @xmath284 , @xmath302 ) using both local and global moves and different ensembles as indicated . for the plj model the values @xmath54 obtained in the @xmath37  ensemble ( squares ) are seen to be essentially identical for all @xmath5 to the moduli @xmath48 ( crosses ) and @xmath151 ( triangles ) obtained in ensembles with strain fluctuations . as shown by the small filled triangles , @xmath54 vanishes in the @xmath50  ensemble for all @xmath5 in agreement with eq . ( [ eq_muamuftau ] ) . ( we remind that the impulsive correction @xmath200 has to be taken into account if @xmath54 is used . ) decreasing the temperature further below @xmath9 the shear moduli are seen to increase rapidly for both models . as shown by the solid and the dashed lines both models are well fitted by a cusp - singularity @xmath369 with fit constants @xmath370 for the ka model and @xmath371 for the plj model where we have set in both cases @xmath372 . confirming the md simulations by barrat _ @xcite , this suggests that the transition is very sharp , albeit _ continuous _ in qualitative agreement with the predictions from ref .  @xcite . if we fit the data close to @xmath9 with an additional off - set @xmath373 , a slightly _ negative _ value is obtained which is not compatible with mct @xcite . however , admittedly both the number of data points close to @xmath9 and their precision ( due to the small number of independent configurations ) are yet not sufficient to rule out completely a positive , albeit ( presumably ) small off - set .    a slightly different representation of the data is given in fig . [ fig_gt_rescal ] where the reduced shear modulus @xmath374 is plotted as a function of the reduced temperature @xmath375 for both models . using the independently determined glass transition temperature @xmath9 and affine shear elasticity @xmath8 as scales to make the axes dimensionless , the rescaled data are found to collapse ! this is of course a remarkable and rather unexpected result considering that two different models in two different dimensions have been compared . the bold line indicates the cusp - singularity @xmath376 with a prefactor which is compatible to the ones used in eq . ( [ eq_cuspcont ] ) considering the typical values of @xmath8 in both models . whether this striking collapse is just due to some lucky coincidence or makes manifest a more general universal scaling is impossible to decide at present . ( please note that the corresponding data for the compression modulus shown in fig . [ fig_kt ] does not scale . ) no rescaling of the vertical axis is necessary for the squared correlation coefficient @xmath377 shown in the inset for two different ensembles of the plj model . the same continuous cusp - like decay ( bold line ) is found as for the modulus @xmath0 . this is again consistent with the thermodynamic reasoning put forward at the end of sec . [ theo_strainfluctu ] , eq . ( [ eq_cmodgt_mua ] ) . [ [ sampling - time - dependence . ] ] sampling time dependence . + + + + + + + + + + + + + + + + + + + + + + + + +    a word of caution must be placed here . the results presented in figs .  [ fig_muamuf_t]-[fig_gt_rescal ] correspond to the longest simulation runs we have at present been able to perform . obviously , this does _ not _ mean that they correspond to the limit for asymptotically long sampling times , if the latter limit exists , or at least an intermediate plateau of respectable period of time . focusing on the plj model in the @xmath50-ensemble we attempt in fig . [ fig_gt_t ] to give a tentative characterization of the sampling time dependence . we plot @xmath378 as a function of @xmath32 for different temperatures as indicated . to make the data comparable all simulations for @xmath379 have been performed with local mc jumps of same maximum distance @xmath272 . as we have already seen in fig . [ fig_gt_thigh ] , the shear modulus decays inversely with time in the liquid regime at high temperatures . this limit is indicated by the thin line . with decreasing @xmath5 the ( unscaled ) shear modulus increases and a shoulder develops . note that even for @xmath298 , i.e. slightly below the glass transition temperature @xmath380 estimated via dilatometry ( fig . [ fig_quench_rho ] ) , @xmath378 decays strongly with time . a more or less @xmath32-independent shoulder ( on the logarithmic scales used ) can only be seen below @xmath273 . we emphasize that we present here a sampling time effect expressing the fact that the configuration space is more and more explored with increasing @xmath32 and _ not _ a time correlation function @xcite . please note also that this time dependence has nothing to do with equilibration problems or ageing effects . time translational invariance is perfectly obeyed in our simulations as we have checked by rerunning the sampling _ after _ additional tempering over at least @xmath347 mcs for all temperatures . similar @xmath32-dependencies for different temperatures @xmath5 have also been observed for @xmath378 and @xmath339 for the plj model and for @xmath339 for the ka model ( not shown ) .    a different representation of the data for both models is shown in fig .  [ fig_gt_t ] where the reduced shear modulus @xmath381 is plotted against the reduced temperature @xmath382 for different sampling times @xmath32 as indicated . similar behavior is found for both models . the shear modulus decreases with @xmath32 and this the more the larger the temperature @xmath5 . for smaller temperatures the different data sets appear to approach more readily as one expects from the @xmath32-independent shoulder for the plj model in fig . [ fig_gt_t ] . as a consequence the glass transition is seen to become sharper with increasing @xmath32 . the data approach the continuous cusp - singularity indicated by the bold line . however , it is not clear from the current data if the data converge to an at least _ intermediately _ stable behavior . longer simulation runs are clearly necessary for both models to clarify this issue . similar data have also been obtained for the plj model from the independent determination of @xmath378 and @xmath383 in the @xmath49-ensemble and for @xmath339 in the @xmath37- and @xmath155-ensemble . the stress fluctuation formalism is a powerful method for computing elastic moduli for canonical ensembles using simulation boxes of constant volume and shape @xcite . focusing on the shear modulus @xmath0 of two well - known glass - forming colloidal model systems in @xmath1 and @xmath2 dimensions , which we have sampled by means of md and mc simulations , we have addressed the general question of whether the stress fluctuation method may be used through a solid - liquid transition where a well - defined displacement field ceases to be defined . deliberately setting the experimentally motivated linear regression eq . ( [ eq_gmodgt ] ) as the fundamental definition , the shear modulus @xmath0 has been computed from the strain and stress fluctuations ( assuming a fixed finite measurement time window @xmath384 ) in ensembles where either a shear strain @xmath36 ( @xmath37- and @xmath155-ensembles ) or a ( mean ) shear stress @xmath41 ( @xmath49- and @xmath50-ensembles ) have been imposed ( fig . [ fig_sketch ] ) . working at constant mean normal pressure we have computed and compared the temperature dependence for various simple averages and fluctuations contributing to the shear modulus in different ensembles . for the ka model only the currently available results for @xmath37- and @xmath155-ensembles have been presented .    as has been stressed in sec . [ theo_stressfluctu ] and sec . [ theo_legendre ] , the stress fluctuation representation @xmath54 at fixed shear strain does not rely _ in principle _ on a solid - like reference state for a displacement field of tagged particles if thermodynamics can be assumed to hold . it thus holds formally through the glass - transition temperature @xmath9 up to the liquid regime ( albeit with a trivial value @xmath34 ) as does the better known relation for the compression modulus @xmath20 ( appendix  [ app_k ] ) @xcite . as emphasized in sec . [ theo_trunc ] , impulsive corrections due to the truncation of the pair potentials have to be properly taken into account , especially at high temperatures ( fig .  [ fig_gt_thigh ] ) , for the precise determination of the affine shear elasticity @xmath255 @xcite . confirming eq . ( [ eq_gmods ] ) , it has been shown for the plj model that @xmath54 at @xmath36 has the _ same _ long - time behavior as the conceptually more direct observables @xmath48  and @xmath151   obtained using the strain fluctuations at constant mean shear stress @xmath41 . the standard thermodynamic relations comparing different ensembles appear thus to hold also _ in practice _ ( at least within our numerical precision ) through the glass transition for all @xmath5 and this albeit    1 . some degrees of freedom get quenched at low temperatures on the time window @xmath384 computationally accessible , 2 . all operational definitions of @xmath0 are transient in the sense that they vanish for @xmath385 ( figs . [ fig_gt_t ] and [ fig_gt_t ] ) and 3 . it is not self - evident that in the latter limit and for large temperatures @xmath3 and @xmath4 can still be treated as a pair of conjugated thermostatistical variables . the latter point compactly epitomized by the thermodynamic relation eq . ( [ eq_dgamdtau ] ) has been shown to hold , however , with remarkable precision up to very high temperatures for the plj model .    as predicted by general legendre transformation , @xmath54 vanishes for all temperatures @xmath5 , if @xmath4 rather than @xmath3 is imposed ( figs . [ fig_gt_dynmat ] ,  [ fig_gt_tlow ] and [ fig_gt ] ) . the shear - stress fluctuations @xmath174 are thus given by the affine response under an external load @xmath26 ( fig . [ fig_muamuf_t ] ) , i.e. @xmath174 reduces to a simple two - point pair correlation function if pair potentials are considered as discussed in sec . [ theo_legendre ] . the same holds above @xmath9 for constant @xmath3 , since the boundary conditions are irrelevant for the liquid state . this symmetry with respect to the boundary conditions ( @xmath386 ) is broken at the glass transition for a large , but finite @xmath384 : at constant @xmath3 the stress fluctuations reveal a strong non - monotonous behavior with a clear maximum at @xmath9 ( fig . [ fig_muamuf_t ] ) . the shear modulus is the order parameter characterizing this symmetry breaking . alternatively , as we have discussed in sec . [ theo_legendre ] , @xmath387 may be seen as an order parameter comparing the ratio of the non - affine to the affine shear responses . since @xmath388 for @xmath389 , this implies that the stress fluctuations contain higher - order correlations _ reducing _ the free energy of the system . the increase of @xmath0 below @xmath9 is reasonably fitted for both models by a _ continuous _ cusp singularity , eq . ( [ eq_cuspcont ] ) , in qualitative agreement with recent replica calculations @xcite . a jump discontinuity , as suggested by mode - coupling theory @xcite and another replica theory @xcite , is not compatible with the currently available data shown in fig . [ fig_gt ] and fig . [ fig_gt_rescal ] . however , as shown in fig . [ fig_gt_t ] and fig . [ fig_gt_t ] , our data depend strongly on sampling time @xmath32 and the computation of larger @xmath32 could lead to a sharper transition . thus a jump discontinuity can not be ruled out completely . ( we are not aware of any other investigation of the sampling time effect for the shear modulus . ) at present we believe , however , that it would be surprising if our data could be reconciled with a discontinuity at @xmath9 . it seems more likely that below @xmath9 , where the choice of the boundary conditions does matter as shown , the definition of the glassy shear modulus used in refs . @xcite do not correspond to the key operational definition @xmath48 of the shear modulus considered by us . in any case , our work shows that any theory and numerical scheme used to determine the shear modulus using a stress or strain fluctuation relation ( in the low - wavevector limit or at a small finite @xmath279 ) must cleary specify and take into account whether the shear stress or the shear strain are macroscopically imposed . otherwise , such an approach is void . [ [ beyond - the - presented - work . ] ] beyond the presented work . + + + + + + + + + + + + + + + + + + + + + + + + + +    future work should clearly focus on the more detailed description of the sampling time dependence . by extrapolating appropriately for the large-@xmath32 asymptotic behavior , this may allow to settle the theoretical debate . as emphasized above , one shortcoming of the present work is that the time - scale used in our mc simulations was slightly arbitrary , especially if additional non - local particle moves are used in the low-@xmath5 limit and if box volume and shape changing strains are included . the comparison of @xmath32-dependent properties for different temperatures @xmath5 becomes thus delicate . our mc study for the 2d soft particles should thus in any case be recomputed with langevin thermostat md dynamics as used for the data of the ka model presented . the @xmath32-dependence of the latter model , which is a better glass former than the 2d plj model @xcite , has still to be worked out for the key operational definition @xmath48 in the @xmath49- or @xmath50-ensembles @xcite . we plan also in the near future to reconsider more carefully our previous investigation of the glass transition of polymer melts @xcite taking properly into account the impulsive truncation corrections and comparing the shear moduli @xmath390 ( @xmath50-ensemble ) and @xmath54 ( @xmath50-ensemble ) around @xmath9 . in addition , it should be rewarding to compare our shear moduli with the values obtained from the displacement of particles following the procedure chosen in ref . this procedure relies on the idea that the particles fluctuate around a well - defined reference position . in the low-@xmath5 limit where the interaction network can be replaced by the dynamical matrix , this should yield the same results . however , for larger temperatures where the harmonic approximation must break down , this approach is questionable ( both for the analysis of experimental and computational data ) and should be compared to the moduli obtained directly from the strain and stress fluctuations of the overall simulation box .    [ [ a - simple - scalar - model . ] ] a simple scalar model . + + + + + + + + + + + + + + + + + + + + + +    we are currently working on an extremely simplified scalar model for permanent and transient networks in @xmath2 dimensions which may help to clarify the theoretical debate @xcite . as sketched in the panel ( a ) of fig . [ fig_band ] , in this model beads ( filled spheres ) are glued on rigid rods which are only allowed to move horizontally . in its most simple implementation this is done by means of global mc moves using transverse plane waves with a wavevector pointing in the @xmath140-direction . the glued beads are connected by ideal harmonic springs where the spring constants @xmath322 and reference lengths @xmath323 are chosen randomly according to fixed , narrow or broad distributions . the stiffness of the contacts between rods may differ ( even including negative @xmath322 and @xmath323 as for the dynamical matrix studied in sec . [ res_dynmat ] ) . the shear modulus may be computed either from the shear strain fluctuations in the @xmath49- or the shear stress fluctuations in the @xmath37-ensemble . the shear stresses can be relaxed by either allowing the particles to move along the rods using longitudinal plane waves in the @xmath105-direction and/or by breaking and reconnection of springs . the latter step is currently done by applying a chemical potential for the springs . physically , the study of such transient networks is motivated by systems of hyperbranched polymer chains with sticky end - groups @xcite or telechelic polymers in water - oil emulsions @xcite . the associated relaxation frequency @xmath391 for a given spring @xmath130 is assumed to be set by an activation barrier @xmath392 which may differ for each spring . depending on the distribution of this activation barrier the shear modulus @xmath349 thus decays more or less rapidly with sampling time @xmath32 as sketched in panel ( b ) of fig . [ fig_band ] . the shear modulus becomes constant , if the random percolating network is permanent ( @xmath393 , @xmath394 ) or the relaxation negligible ( top solid line ) . the system behaves as a liquid ( bottom line ) for large ( mean ) @xmath395 , while a shoulder is seen in the intermediate frequency range . while this is confirmed by preliminary simulations , we do at present not understand the role of the quenched noise injected in the model which is seen to strongly influence the shear stress fluctuations @xmath55 ( especially if fluffy layers happen to appear ) . also the understanding of finite system - size effects is of course crucial for such a low - dimensional model . ( only systems containing @xmath396 particles have been considered so far . ) as sketched in panel ( c ) the key question is to describe accurately for asymptotically long simulation runs and asymptotically large box sizes the shape of @xmath35 around the glass transition temperature @xmath397 where the shear modulus vanishes and @xmath348 should again have a maximum . thanks the cnrs and the irtg soft matter for supporting her sabbatical stay in strasbourg , p.p . the rgion alsace and the irtg soft matter and f.w . the daad for funding . we are indebted to a.  blumen ( freiburg ) and h.  meyer , o.  benzerara and j. farago ( all ics , strasbourg ) for helpful discussions . we have focused in the main part of this paper on the shear modulus @xmath0 of our glass - forming model systems and the comparison of the three observables @xmath48 , @xmath151 and @xmath54 . being isotropic the complete macroscopic elastic description of our systems only requires the numerical determination of _ one _ additional elastic modulus , the isothermal compression modulus @xcite , @xmath398 characterizing the elastic response with respect to a volumetric ( dilatational ) strain fluctuation as shown in fig . [ fig_sketch](d ) . as discussed in the standard textbooks @xcite , the extensive variable is here the volume @xmath399 and the conjugated intensive variable the negative pressure @xmath400 where @xmath401 and @xmath402 stand for the instantaneous values of the volume and the pressure . with these notations eq . ( [ eq_pkv ] ) is obviously consistent with eq . ( [ eq_ftom ] ) or eq . ( [ eq_htom ] ) . in analogy with @xmath48 for the shear modulus @xmath0 , the compression modulus @xmath20 may be determined in an @xmath155- or an @xmath50-ensemble using the linear regression relation @xmath403 where the index @xmath404 indicates that both the measured instantaneous volume @xmath401 and pressure @xmath402 have been used . the corresponding correlation coefficient is @xmath405 since from eq . ( [ eq_dxdi ] ) we have @xmath406 if @xmath6 is imposed , this implies that @xmath407 should be identical to the better known relation @xcite @xmath408 which corresponds to the operational definition @xmath151 discussed above . the compression modulus @xmath20 may of course also be determined using the pressure fluctuations in systems of imposed volume @xmath12 as sketched in fig .  [ fig_sketch](c ) . as ( to our knowledge ) first derived by rowlinson @xcite , the corresponding stress fluctuation formula reads @xmath409 being the fluctuation of the total normal pressure ( corresponding to @xmath55 ) and @xmath259 the  affine dilatational elasticity " ( corresponding to @xmath26 ) . we shall give below a proper operational definition for @xmath259 . here we only note that due to the general legendre transformation eq . ( [ eq_didi ] ) for intensive variables one knows immediately that @xmath410 this implies by comparison with eq . ( [ eq_kmodpp ] ) that @xmath411 in analogy to eq . ( [ eq_muamuftau ] ) for the shear , i.e. @xmath259 can be directly be determined from the pressure fluctuations in @xmath155- and @xmath50-ensembles . as before for the correlation coefficient @xmath90 , eq . ( [ eq_cmodgt_mua ] ) , it follows from eq . ( [ eq_etaaetafp ] ) that @xmath412 since @xmath413 , the affine dilatational elasticity sets an upper bound to the compression modulus @xmath20 . we reformulate here briefly the derivation of the stress fluctuation relation @xmath258 for the compression modulus @xmath20 of an isotropic system at imposed volume @xmath12 by rowlinson @xcite . with @xmath414 being the volume of the unperturbed simulation box we assume a small relative volume change @xmath415 with all coordinates equally deformed as , e.g. @xmath416 for the @xmath105-coordinate in @xmath19-dimensions . ( the argument @xmath99 denotes again the reference system . ) using a similar rescaling trick as for the shear modulus @xmath0 in sec . [ theo_stressfluctu ] , the interaction energy @xmath417 of a configuration @xmath98 of the strained system is expressed in terms of the coordinates ( state ) of the unperturbed system and the explicit metric parameter @xmath229 . due to the volume change we have now also to take into account the ideal gas ( kinetic ) contributions . we derive first the excess contribution @xmath29 to the total pressure @xmath418 and the excess contribution @xmath419 to the total compression modulus @xmath420 valid for an arbitrary conservative potential . the general relations eqs . ( [ eq_dlogzdgam]-[eq_ddzddgam ] ) stated above for an imposed shear strain @xmath3 still hold with the relative volume change @xmath229 replacing @xmath3 . using eq . ( [ eq_dzdgam ] ) it follows for the excess pressure that @xmath421 _ defining _ the instantaneous excess pressure @xcite . in the second step we have taken the limit @xmath422 and have dropped the argument @xmath99 . the average is again evaluated using the weights of the unperturbed system , eq . ( [ eq_average ] ) . using eq . ( [ eq_ddlogzddgam ] ) and eq . ( [ eq_ddzddgam ] ) one obtains for the compression modulus @xmath423 which is again understood to be taken at @xmath424 . the first term @xmath425 in eq . ( [ eq_kexstressfluctu ] ) corresponds to the change of the system energy assuming an affine strain transformation for all particle positions . it gives the excess contribution @xmath426 to the total affine dilatational elasticity @xmath427 . the second term corresponds to the excess contribution @xmath428 to the total normal pressure fluctuation @xmath429 mentioned in sec . [ app_kstrainfluctu ] . it corrects the overprediction of the affine strain contribution @xmath426 @xcite . we remind that for an ideal gas the isothermal compression modulus @xmath430 is given by the ideal gas pressure @xmath145 and that the ideal gas pressure fluctuation @xmath431 becomes @xmath432 at constant volume @xcite . interestingly , using again the general legendre transform eq . ( [ eq_didi ] ) one sees that @xmath433 for the ideal gas pressure fluctuations at imposed pressure .    due the assumed additivity of the total system hamiltonian , the stress fluctuation formula for the total compression modulus may now be written @xmath434 while this result is useful for constant-@xmath12 ensembles , especially if mc simulations are used , it is yet not in a form exactly equivalent to the general legendre transform eq . ( [ eq_kmodpp ] ) . since the ideal and excess pressure fluctuations are decoupled , @xmath435 , this implies @xmath429 . one may thus rewrite eq . ( [ eq_kconstantv ] ) as @xmath436 with @xmath437 where we have used eq . ( [ eq_etafidp ] ) in the second step . assuming eq . ( [ eq_etaaexid ] ) assures that @xmath259 becomes equivalent to @xmath438 . this has the nice feature that not only the functional @xmath258 , eq . ( [ eq_kmodpp ] ) , should give the correct compression modulus @xmath20 in the @xmath37  and the @xmath49  ensembles , but also that @xmath258 must vanish properly if the  wrong " @xmath155  and @xmath155  ensembles are used , just as @xmath54 vanishes for @xmath49- or @xmath50-ensembles .      up to now we have not taken advantage of the fact that the conservative interaction potential is assumed to be a pair potential , eq . ( [ eq_upairinteraction ] ) . this will be used now to express @xmath29 and @xmath426 in terms of simple ( ensemble independent ) two - point correlation functions . equation  ( [ eq_strain_x_epsilon ] ) implies that the squared distance @xmath114 between two particles transforms as @xmath439 it follows that @xmath440 where we have taken finally for both derivatives the limit @xmath422 . for the first two derivatives of a general function @xmath441 with respect to @xmath229 this implies @xmath442 for @xmath422 and dropping the argument @xmath99 on the right - hand sides . according to the general result eq . ( [ eq_peshat_us ] ) for the normal pressure and using eq . ( [ eq_strain_fronebeps ] ) one confirms for pair interactions the expected virial relation @xcite @xmath443 for the instantaneous excess pressure where the interactions between pairs of particles @xmath133 are again labeled by an index @xmath130 . using eq . ( [ eq_strain_rtwo_epsilon ] ) the excess contribution to the affine dilatational elasticity becomes @xmath444 is sometimes called the  hypervirial " contribution to the compression modulus @xcite . it corresponds to the first term indicated in eq . ( [ eq_strain_frtwobeps ] ) . by comparing eq . ( [ eq_etab ] ) with eq . ( [ eq_mub ] ) it is readily seen that for isotropic systems @xmath445 i.e. @xmath256 and @xmath28 and , hence , the affine dilatational elasticity @xmath259 and affine shear elasticity @xmath26 are closely related . ( this justifies to call @xmath256 a  born coefficient " . ) please note that eq . ( [ eq_kmodpp ] ) or eq . ( [ eq_kconstantv ] ) together with eq . ( [ eq_etaaex ] ) are perfectly consistent with the relations stated in the literature @xcite . we emphasize finally that , as discussed for the affine shear elasticity @xmath26 in sec . [ theo_stressfluctu ] , it is inconsistent to neglect the explicit excess pressure contribution to @xmath259 in eq . ( [ eq_etaaex ] ) but to keep the underlined term of the born coefficient @xmath256 , eq . ( [ eq_etab ] ) @xcite . the affine dilatational elasticity @xmath446 obtained for both models is presented in the inset of fig . [ fig_kt ] as a function of the reduced temperature @xmath447 . ( the values may also be found in table  [ tab_ka ] and table  [ tab_plj ] . ) we remind that for determining the born coefficient @xmath256 impulsive corrections need to be considered as discussed at the end of sec . [ theo_trunc ] . as shown by the small filled triangles , @xmath259 compares well with the total pressure fluctuation @xmath448 in an @xmath50-ensemble . the observed data collapse is expected from eq . ( [ eq_etaaetafp ] ) and eq . ( [ eq_etafidp ] ) . as indicated by bold and dashed lines for , respectively , the low- and high-@xmath5 limit @xmath259 decreases essentially linearly with temperature . the coefficients of the lines ( not shown ) are consistent with eq . ( [ eq_etabmub ] ) relating @xmath259 and @xmath26 and the fits given in fig . [ fig_muamuf_t ] for @xmath26 . we emphasize that the affine dilatational elasticity @xmath259 ( and thus @xmath26 ) is not only larger for the ka model but , more importantly that , is also increases more rapidly with decreasing @xmath5 than the plj model . the main panel of fig . [ fig_kt ] shows the best long - time values of the compression modulus @xmath20 for both models . the vertical axis is rescaled by @xmath259 . the rescaled stress fluctuation formula , eq . ( [ eq_kconstantv ] ) , for the @xmath37-ensemble is given for both models . for the plj model we also indicate @xmath407 ( crosses ) and @xmath449 ( large triangles ) for the @xmath50-ensemble . no rescaling with @xmath259 is necessary for the squared correlation coefficient @xmath450 ( stars ) also presented . a perfect collapse of all plj data is found confirming the equivalence of all operational definitions used . for both models we observe a strong increase of @xmath451 around @xmath452 with decreasing temperature . interestingly , the rescaled data of both models is similar but not identical . the mentioned increase around @xmath452 is clearly more pronounced for the ka model . note that if the unscaled @xmath20 is directly plotted against @xmath5 , this effect becomes even stronger ( not shown ) consistently with the scaling of @xmath259 presented in the inset . we note finally that for very low temperatures we observe for both models that the ratio @xmath453 approaches unity , i.e. the compression modulus is essentially dominated by the affine response . we remind that this is different for the reduced shear modulus for both models which is observed to approach @xmath454 in the low-@xmath5 limit ( fig . [ fig_gt_rescal ] ) , i.e. the shear stress fluctuations @xmath364 are more relevant than the pressure fluctuations @xmath455 in this limit . for the plj model we present in fig . [ fig_kt_thigh ] the sampling time dependence of various operation definitions of @xmath20 for the liquid limit at temperature @xmath456 . the horizontal axis indicates the sampling time @xmath32 in units of mcs , the vertical axis is made dimensionless using @xmath259 as natural scale . as shown for @xmath407 and @xmath449 obtained in the @xmath50-ensemble and for @xmath258 obtained using eq . ( [ eq_kconstantv ] ) in the @xmath37-ensemble , all operational definitions applied in their natural ensemble converge _ rapidly _ within @xmath457 mcs to the same plateau . ( this is a general finding for all temperatures . ) the plateau value for @xmath456 is indicated by the bold solid line . note that @xmath451 is smaller unity as expected from the discussion at the end of sec . [ app_kstrainfluctu ] . the small filled symbols refer to @xmath258 computed in the  wrong " @xmath50-ensemble . the small diamonds have been computed using directly eq . ( [ eq_kconstantv ] ) . as indicated by the dashed line , this yields @xmath458 for large times . only if we compute @xmath436 using eq . ( [ eq_etaaexid ] ) , it does vanish properly as shown by the small triangles . this shows numerically that eq . ( [ eq_kconstantv ] ) , albeit perfectly fine for constant-@xmath12 ensembles , is just not the correct formulation allowing to make manifest the general legendre transformation eq . ( [ eq_etaftransform ] ) . this is , however , achieved using eq . ( [ eq_etaaexid ] ) . this finding holds rigorously for all temperatures as can be seen from the scaling of the small triangles in the inset of fig . [ fig_kt ] .
Abstract:

Answer: the shear modulus @xmath0 of two glass - forming colloidal model systems in @xmath1 and @xmath2 dimensions is investigated by means of , respectively , molecular dynamics and monte carlo simulations . comparing ensembles where either the shear strain @xmath3 or the conjugated ( mean ) shear stress @xmath4 are imposed , we compute @xmath0 from the respective stress and strain fluctuations as a function of temperature @xmath5 while keeping a constant normal pressure @xmath6 . the choice of the ensemble is seen to be highly relevant for the shear stress fluctuations @xmath7 which at constant @xmath4 decay monotonously with @xmath5 following the affine shear elasticity @xmath8 , i.e. a simple two - point correlation function . at variance , non - monotonous behavior with a maximum at the glass transition temperature @xmath9 is demonstrated for @xmath7 at constant @xmath3 . the increase of @xmath0 below @xmath9 is reasonably fitted for both models by a continuous cusp singularity , @xmath10 , in qualitative agreement with some recent replica calculations . it is argued , however , that longer sampling times may lead to a sharper transition . the additive jump discontinuity predicted by mode - coupling theory and other replica calculations thus can not ultimately be ruled out .