the theoretical description of slow dynamics is a crucial point to elucidate the nature of the glass transition in structural glass - forming liquids .
one of the commonly used approaches , mode - coupling theory ( mct ) , was from the very beginning designed for the supercooled simple liquids @xcite , i.e. for the non - disordered models ( as opposed to the models which contain quenched disorder naturally ) .
later it has been proven that mc - equations become exact for a number of spin - glass models @xcite as well as for the polymeric manifold in a random media @xcite ( i.e. for the models with quenched disorder ) provided that the number of variables components goes to the infinite .
actually applicability of the mc - equations has been substantially extended to the case when the time translation invariance and the fluctuation - dissipation theorem do not hold any more @xcite .
this striking similarity between the models with and without quenched disorder suggests that the effective disordered potential ( e.g. , in a supercooled liquid ) is in a sense `` self - induced '' and the difference between such a `` self - induced disorder '' and the quenched disorder might not be crucial @xcite .    in order to provide some insight into self - induced disorder we employed in ref.@xcite a feynman variational principle ( vp ) for a set of interacting particles .
indeed it was shown that the vp is capable to treat metastable states of the glass - forming system . the main point in ref .
@xcite was that the partition function representation in terms of functional integrals is twofold : ( a ) either as an integral over the local density , @xmath3 , or ( b ) over the conjugated to @xmath4 field @xmath5 .
it has been _ assumed _ that the _ component average _ free energy , @xmath6 , ( which is only meaningful in the supercooled regime ) is equal to the variational free energy @xmath7 .
there are at least four strong reasons in favor of this ( at first sight not obvious ) assumption .
* the variational free energy @xmath8 is a upper bound for the _ canonical _ free energy , i.e. , @xmath9 , as it should be since @xmath10 , where the complexity @xmath11 @xcite . * after implementation of vp the initial problem is reduced to a self - consistent random field ginzburg - landau model ( rfglm ) .
then , as was shown previously , the corresponding field @xmath12 must be upgraded to a _ replicated _ field @xmath13 , where @xmath14 ( with the final limit @xmath15 ) and the density field @xmath4 plays the role of an `` external '' field .
note that here the density @xmath16 is gaussian due to the use of the vp .
eventually the correlators of @xmath17 and @xmath18 - fields can be determined self - consistently . *
the resulting replicated partition function for rfglm has a typical form which may eventually lead to the replica symmetry breaking ( rsb ) , structural glass transition and the  self - induced  disorder . *
finally , in the case of the long - range interaction the partition function allows the expansion around the saddle point , or mean - field ( mf ) , solution .
it is possible to show then that the _ next to the mean - field approximation _ and vp merge and both become exact , i. e. @xmath19 and the glassy phase does not appear .
some evidence for this behavior was deduced from the results for the particles on a @xmath20- dimensional hypersphere @xcite at large dimensions , @xmath21 .
the aim of this paper is to face the full dynamical problem for a non - disorded model with a long range interaction . using the expansion around the saddle point solution we derive the full equation of motion for
the time dependent density - density correlator and show that a `` glassy '' solution does not exist .
conversly , if we add a term describing quenched disorder , by random distribution of the strength of the interaction potential , then the resulting equations of motion for two time density correlation and response functions fall in the same class as mc - equations which have been widely discussed @xcite .
this means that the `` self - induced '' disorder is not generic for the pure model with the long - range interaction and vice versa on addition of a quenched disorder the phase space becomes very rugged resulting in slow dynamical processes .
the paper is organized as follows . in the next section [ ii ]
we first introduce the theorotical model without quenched disorder .
its dynamics is discussed by using the functional integral technique .
the saddle point solution yields the mean field dynamcis .
expansions around the saddle point yield one loop corrections .
the legendre transformation provides the possibilty of the analysis of the full dynamic correlation matrix . in the section [ iii ]
quenched disorder is introduced by a `` random bond model '' and a gaussian disorder .
the corresponding generating functional ( gf ) is computed by the self - consistent hartree approximation , which results in a set of coupled langevin equations , that are solved in their asymptotic regimes .
more details on the calculations are laid out in the corresponding appendices .
we start from a simple model system which consists of interacting particles . to do so , let us consider a set of @xmath22 particles in @xmath23 - dimensional space interacting by a pair potential of the form @xmath24 this is a typical example of a long - range potential with a characteristic length @xmath25 and a coupling constant @xmath26 .
the choice of this potential is twofold .
it contains a cut off at a @xmath27 and allows thus to control the range of the interaction .
moreover at small scales @xmath28 it consists of a typical power law decay with long range character , if @xmath29 .
therefore the so chosen potential allows to keep control on range and nature of the interaction , which will become essential below . to ensure extensivity of the total interaction energy we require that the integral @xmath30 , i.e. it does not depend on the number of particles @xmath0 . as a result
we have @xmath31 .
the intermolecular potential ( [ 1 ] ) has the form of the generalized kac potential @xmath32 which has been used for the rigorous treatment of the van - der - waals theory @xcite . in order to provide conditions for the expansion around a saddle point , carried out later on ( see below )
, we should require that the length @xmath27 must be larger compared to the characteristic size of the system , which scales naturally as @xmath33 ) at @xmath34 . as a consequence we find the limits for the range parameter @xmath35 @xmath36 below we shall restrict our considerations to the case : @xmath37 , and the strength of the interaction @xmath38 ( pure repulsion ) without loss of generalization in the main statements we are going to predict .
then the fourier transformation of the potential ( [ 1 ] ) takes especially simple form @xmath39 which allows accurate analytic calculations . in the limit @xmath34
we have thus @xmath40 , but the relevant minimum wave vector is @xmath41 and thus @xmath42 can be actually neglected under the integration over the whole @xmath43 - space . as a result
we arrive formally at a _
one - component plasma model _ ( ocp )
@xcite where the electroneutrality is implicitly provided by a neutralizing background .      in the following we set up the relevant equations of motion for the model system . we restrict ourselves to the langevin dynamics , which can be comfortably formulated in terms of dynamic functionals which allows the systematic @xmath1 - expansion treatment .
the langevin dynamics of @xmath0 particles interacting via the potential ( [ 1 ] ) ( at @xmath37 and @xmath38 ) is described by the equation of motion @xmath44 where @xmath45 and @xmath46 are the mass and the friction coefficient respectively , @xmath47 and @xmath48 .
the random force in eq.([4 ] ) is gaussian with @xmath49 and the correlator @xmath50 where from now on we work in units where the boltzmann constant @xmath51 .
as was mentioned , it is more convenient to reformulate the langevin problem ( [ 4])-([5 ] ) by using the celebrated martin - siggia - rose generating functional ( gf ) method @xcite .
the method was first applied for the @xmath52 - model with the long - range interaction in @xcite and for the polymer melt dynamics in @xcite . despite the fact that the langevin equation ( [ 4 ] ) is of the second order it is possible to show that the jacobian which appear under transformation to the functional variables , still equal to one ( see appendix in @xcite ) . after using this technique for the problem ( [ 4 ] ) - ( [ 5 ] )
, gf takes the form @xmath53 + \int dt \sum_{p=1}^{n}\sum_{m=1}^{n } \frac{\mu}{n }    i{\hat r}_{j}^{(p)}(t ) \nabla_{j}^{(p)}v\left({\bf r}^{(p ) } - { \bf        r}^{(m)}\right)\right\ } \ ; , \label{6}\end{aligned}\ ] ] where the action of the free system @xmath54 = \int dt \left\ { t\gamma_{0 }    \left[i{\hat r}_{j}^{(p)}(t)\right]^{2 } + i{\hat      r}_{j}^{(p)}(t)\left[m_{0}\frac{\partial^2}{\partial t^2 }      { r}_{j}^{(p)}(t ) +   \gamma_{0}\frac{\partial}{\partial t }      { r}_{j}^{(p)}(t )   \right]\right\ } \ ; . \label{7}\end{aligned}\ ] ] in the following we are going to transform this functional to collective density variables . by using the transformations to the mass density @xmath55 and the longitudinal projection of the response field density @xmath56 for the gf one gets @xmath57 where the summation over the repeated greek indices is implied . in eq.([10 ] ) we have introduced 2-dimensional field @xmath58 where @xmath59 and @xmath60 .
the `` entropy '' of the free system is given as usual by @xmath61\nonumber\\ & \times&\delta\left[\pi({\bf r},t)-\sum_{p=1}^{n}\:i{\hat r}_{j}^{(p)}(t)\nabla_{j}\delta({\bf r}-{\bf    r}^{(p)}(t))\right ] \ ; , \label{12}\end{aligned}\ ] ] and @xmath62 is the 2@xmath63 2 - interaction matrix @xmath64 and @xmath65 is a source field .
an alternative valuable representation of gf can be obtained through the `` functional fourier transformation '' @xmath66 and its inversion @xmath67 the substitution of eq.([12 ] ) into eq.([14 ] ) leads to the explicit expression for the free - system  gf @xmath68\nonumber\\   & - & i\sum_{p=1}^{n}\int dt \psi \left({\bf r}^{(p)}\right )    + \left .
i \hat r_{j}^{(p)}(t)\nabla_{j}\phi ( { \bf      r})\right|_{{\bf r}= { \bf r}^{(p)}(t ) } \bigg\}\ ; , \label{16}\end{aligned}\ ] ] where @xmath69 and @xmath70 are components of the column - variable @xmath71 by making use ( [ 15 ] ) in ( [ 10 ] ) and after functional integration over @xmath72 one gets @xmath73\bigl[u^{-1}\bigl]_{\alpha\beta } ( 1,2)\bigl[i\psi_{\beta}(2)+ \chi_{\beta}(2)\bigl ] + \int d1\psi_{\alpha}(1)\lambda_{\alpha}(1)\bigg\}\ ; , \label{18}\end{aligned}\ ] ] where we have also add a source field @xmath74 conjugated to @xmath75 . as a result eqs.([10 ] ) and ( [ 18 ] ) provide two equivalent representations of gf . for the purpose of expansion around the saddle point we use representation ( [ 18 ] ) at @xmath76
which after the transformation @xmath77 , yields @xmath78\bigg\}\ ; , \label{19}\end{aligned}\ ] ] which is appropriate for a saddle point integration , since the particle number @xmath0 is large .
the action hereby is given as @xmath79&= & \frac{1}{2}\int dt \int d^3{\bf r}d^3{\bf r'}\psi_{\alpha}({\bf    r},t)\left[v^{-1}\right]_{\alpha\beta}({\bf r}- { \bf r ' } ; \kappa ) \psi_{\beta}({\bf r ' } , t ) -   \frac{1}{n}\log \int \prod_{p=1}^{n}d{\bf r}^{(p)}(t)d{\bf \hat r}^{(p)}(t)\nonumber\\ & \times&\exp \bigg\{\sum_{p=1}^{n}a_{0}[{\bf r}^{(p)},{\bf \hat r}^{(p ) } ] - i\sum_{p=1}^{n}\int dt r_{\alpha}^{(p)}(t)\left[\psi_{\alpha}\left({\bf r}^{(p)}(t)\right ) + i\chi_{\alpha}\left({\bf r}^{(p)}(t)\right)\right]\bigg\}\ ; , \label{20}\end{aligned}\ ] ] and the interaction matrix @xmath80 recall that the relation @xmath81 is necessary for the validity of the saddle point integration .
moreover we have defined the column - vector @xmath82 for convenience .
minimization of @xmath83 $ ] with respect to @xmath75 leads to the sp - equations for the mean fields @xmath84 @xmath85 where the average @xmath86 is calculated by using the _ cumulant _ gf @xmath87 -\nonumber\\ i\sum_{p=1}^{n}\int dt r_{\alpha}^{(p)}(t)\left[\overline{\psi_{\alpha}}\left({\bf r}^{(p)}(t)\right ) + i\chi_{\alpha}\left({\bf r}^{(p)}(t)\right)\right]\bigg\}\ ; .\label{24}\end{aligned}\ ] ] the correlation matrix in the random phase approximation ( rpa ) is defined in such a way @xmath88\ ; . \label{25}\end{aligned}\ ] ] after linearization of eq.([24 ] ) with respect to @xmath89 the 2@xmath632 - rpa correlation matrix is easily found to coincide with the well known form @xcite @xmath90^{-1}\right\}_{\alpha\beta}(1 , 2)\ ; , \label{26}\end{aligned}\ ] ] where @xmath91 is the interaction matrix ( [ 21 ] ) and @xmath92 is the correlation matrix for the _ free system _
@xmath93 has the form @xmath94 in eq.([27 ] ) @xmath95 and @xmath96 are response functions whereas @xmath97 stands for the correlation function .
the relation between them is given by the fluctuation dissipation theorem ( fdt ) which in @xmath98- representation has the form @xmath99 it is easy to check that in this case the fdt for the rpa type correlation matrix ( [ 26 ] ) also holds @xmath100 where @xmath101 is the inverse temperature .
the corresponding elements of the rpa - matrix ( [ 26 ] ) are of an especially simple form in the fourier - @xmath102 - representation , namely @xmath103\left[1+v(k)f_{01}({\bf k},\omega)\right ] } \label{rpares1}\\ s_{01}({\bf k},\omega)&=&\frac{f_{01}({\bf      k},\omega)}{1+v(k)f_{01}({\bf k},\omega ) } \label{rpares2}\\ s_{10}({\bf k},\omega)&=&\frac{f_{10}({\bf      k},\omega)}{1+v(k)f_{10}({\bf k},\omega)}\ ; .
\label{rpares3}\end{aligned}\ ] ] it turns out interesting to recover the wellknown form in the static limit , where we have @xmath104 $ ] and for the correlator @xmath105 one gets @xmath106 where we have used @xmath107 .
this expression is completely equivalent to the correlator for the ocp - model ( see eq.(10.1.7 ) in @xcite ) with the direct correlation function @xmath108 and the debye wavenumber @xmath109 .
now let us expand the action ( [ 20 ] ) around sp - solution ( [ 23 ] ) up to the second order with respect to the fluctuations @xmath110 .
after the functional integration we arrive at the following result for the gf @xmath111\ ; , \label{31}\end{aligned}\ ] ] where @xmath112 is the inverse matrix of the effective interactions @xcite @xmath113_{\alpha\beta}(1 , 2 ) +    \frac{1}{n}\left<\delta\rho_{\alpha}(1)\delta\rho_{\beta}(2)\right>_{{\rm sp}}\ ; .\label{32}\end{aligned}\ ] ] in eqs.([31 ] ) - ( [ 32 ] ) we deliberately keep external field @xmath65 nonzero because it to be used in the next subsection for the legendre transformation .
the functional legendre transformation is a general way to provide the dyson equation for the _ full correlation matrix _
in doing so the _ irreducible _ gf , @xmath115 , is defined by the identity @xmath116 by doing functional differentiation of ( [ 33 ] ) one gets @xmath117 and @xmath118_{\alpha\beta}(1 , 2 ) = \frac{\delta^{2}\gamma\left\{\left<\rho_{\alpha}(1)\right>\right\}}{\delta\left<\rho_{\alpha}(1)\right > \delta\left<\rho_{\beta}(2)\right > } \ ; .\label{35}\end{aligned}\ ] ] taking into account the result in eq.([31 ] ) we find the following result for gf @xmath119\ ; , \label{35'}\end{aligned}\ ] ] where @xmath120 in eq.([35 ] ) one should consider @xmath65 as a functional of @xmath121 given by eq.([34 ] ) .
double differentiation of eq.([35 ] ) leads to an equation of the dyson form @xmath118_{\alpha\beta}(1 , 2 ) = \left[s^{-1}\right]_{\alpha\beta}(1 , 2 ) - \sigma_{\alpha\beta}(1 , 2)\ ; , \label{37 } \end{aligned}\ ] ] where the rpa - correlation matrix , @xmath122 , is defined by eqs.([rpares1 ] ) - ( [ rpares3 ] ) and the `` self - energy '' functional @xmath123 has the form @xmath124 in eq.([38 ] ) the `` trace '' is taken over the variables 3,4 and indices @xmath125 .
the explicit differentiation in eq.([38 ] ) leads to the result @xmath126 where @xmath127 is a short hand notation of the matrix @xmath128 and we have took into account that @xmath129 at @xmath130 because the fluctuations are gaussian .
further calculation yields @xmath131 where @xmath132 and the full mean field @xmath133 the expression for @xmath134 can be easily found by differentiation of eq.([42 ] ) with respect to @xmath135 .
taking into account eqs.([23 ] ) , ( [ 34 ] ) and ( [ 35 ] ) at @xmath136 one gets @xmath137_{\alpha\beta}(1 , 2 ) - \mu \int d4 d3 v_{\alpha\omega}(1 , 4 ) s_{\omega\gamma}(4 , 3)r_{\gamma\beta}(3 , 2)\label{43}\end{aligned}\ ] ] or finally @xmath138^{-1}\right\}_{\alpha\gamma}(1 , 3)\left[\hat g^{-1}\right]_{\gamma\beta}(3 , 2)\ ; , \label{44 } \end{aligned}\ ] ] where the hatted variables stands for the corresponding 2@xmath632 matrices .
substitution eqs.([44 ] ) and ( [ 40 ] ) in eq.([39 ] ) yields @xmath139_{\gamma\alpha}(3 , 1)\left[g^{-1}\right]_{\delta\beta}(4 , 2)\;,\label{45}\end{aligned}\ ] ] where 2@xmath632 vertex - matrix has the form @xmath140^{-1}\right\}_{\delta\gamma}(4 , 3 ) s^{(4)}_{\gamma\delta\omega\chi}(3 , 4 , 5 , 6)\nonumber\\ & \times&\left\{\left[\hat 1 + \mu \hat v \hat s    \right]^{-1}\right\}_{\omega\alpha}(5 , 1)\left\{\left[\hat 1 + \mu \hat v \hat s    \right]^{-1}\right\}_{\chi\beta}(6 , 2)\ ; .\label{46}\end{aligned}\ ] ] in eq.([46 ] ) @xmath141 is the 4 - point ( response ) correlator matrix in rpa @xmath142 the explicit calculation of @xmath143 is implemented in appendix b. the vertex - matrix can be seen as a one - loop diagram ( see fig .
+    the higher loops contributions which include generally speaking 2 m - point correlators , @xmath144 , can be also considered , however the `` self - energy '' has still the same convolution structure : @xmath145 . here
the vertex - matrix @xmath146 is calculated in rpa only .
that is why these contributions basically do not change our results .
as a result in the @xmath147- representation the dyson equation ( [ 37 ] ) with the `` self - energy '' functional ( [ 45 ] ) reduces to a quadratic one @xmath148_{\gamma\delta}({\bf k } , \omega)g_{\delta\beta}({\bf k},\omega ) - g_{\alpha\beta}({\bf k } , \omega ) + k_{\alpha\beta}({\bf k } , \omega ) = 0\ ; .\label{48}\end{aligned}\ ] ] the coefficients of the eq.([48 ] ) trace the problem back to the _ free system dynamics _ which is embodied in the correlation matrices @xmath149 and @xmath150 .
it is not surprising then that a specification of the model for the free system dynamics is necessary , before going to the investigation of eq.([48 ] ) .      as we have mentioned the explicit solution of eq.([48 ] )
needs the specification of the free system dynamics .
two simple models are most amenable for the theoretical treatment   : free diffusion model ( fdm ) and the relaxation time approximation model ( rtam  ) @xcite . the latter provide more reasonable dynamical information also for short time intervals , @xmath151 , where the fdm completely failed ( e.g. , the sum rule does not hold ) .
it turns out that upon calculation of the trace in eq.([46 ] ) the integral is ultraviolet - divergent for fdm and only rtam leads to the finite result .
the matrix elements for rtam has the form @xmath152 where we introduced the diffusion coefficient @xmath153 , the characteristic time scale @xmath154 , and @xmath155 for the overall density . at @xmath156
we return to fdm . in the case of rtam
the solution of eq.([48 ] ) for the full correlation matrix reads @xmath157k_{01}({\bf k } ,      \omega)}}{2 \left[-i\omega\tau_{\rm c } -      \omega^{2}\tau_{0 } \tau_{\rm c } + \chi_{st}^{-1}(k)\right]}\label{49'}\\ g_{10}({\bf k},\omega ) & = & g_{01}(-{\bf k},-\omega ) \label{50 } \end{aligned}\ ] ] and @xmath158          k_{01}({\bf k},\omega)}}{- i\omega\tau_{\rm c } - \omega^2 \tau_{0 } \tau_{\rm c } + \chi_{st}^{-1}({\bf k } ) } \right|^2 - k_{00}({\bf k } ,    \omega)}{{\rm re } \left\{\sqrt{1 - 4\left [ - i \omega\tau_{\rm c } -\omega^2 \tau_{0 } \tau_{\rm c } + \chi_{st}^{-1}({\bf k } ) \right ] k_{01}({\bf k } ,    \omega)}\right\}}\ ; . \label{51}\end{aligned}\ ] ] the explicit calculation of the matrix @xmath159 ( see eq.([46 ] ) is given in the appendix c. the overall behavior of the correlation function @xmath160 according eq .
( [ 51 ] ) is shown in fig .
2 ( at @xmath161 ) .
it can be seen clearly that there no singularity appears at @xmath162 .
+    the singularity , however , might be responsible for a glass transition .
instead , the low frequency limit of @xmath160 slowly changed with control parameters ( which is not shown in fig.2 ) .
that means that for the non - disordered model with a general repulsive long ranged potential ( [ 1 ] ) the glass transition is not generic .
this very important conclusion suggest that for the model with a long range interaction the phase space is too smooth to show a glass transition . in order to obtain a glass like transition a competing interactions or a quenched disorder
should be added .
this leads to glassy dynamics , as we will show in the next section .
it is interesting to note , that for the generalized kac potential eq .
( [ kac ] ) , where @xmath163 and its fourier transformation are positive definite functions , mct - memory kernel vanishes at @xmath164@xcite .
the corresponding argumentation is relegated to the appendix d. the explanation for this result lies in the fact that the `` cage effect '' , which is a cornerstone of mct , is missing in the mf - limit .
the `` glass transition '' which has been studied in ref.@xcite for the particles interacted via the kac potential ( [ kac ] ) has a completely different nature . in ref.@xcite
the function @xmath163 has a step form so that its fourier transform @xmath165 is negative at some value of @xmath166 . as a result
the system becomes unstable and a nonuniform configuration where the particles are grouped into `` clumps '' shows up .
it was found that the slow dynamics of the mf - model is associated with these clumps and does not touch a single particle motion .
obviously it is different from the conventional glass transition @xcite .
in the previous sections we have shown in detail that in the absence of disorder the dynamical spectrum changed monotonically with a control parameter and no glassy dynamics can be seen .
the natural question which arises now is : how will the introduction of competing interactions and / or quenched disorder affect the dynamics of the system discussed above ? to provide an answer to this question we will use already existing models of heteropolymers and their disordered two body interaction @xcite .
the use of these models and techniques are here natural , since the behavior of heteropolymers is well discussed in the literature . in principle , two practical possibilities exist .    * the strength of the two- body interaction , @xmath167 , in eq.([1 ] ) is now a random function of all pairs of the interacted particles , @xmath168 ( `` random - bond model '' ) . *
each particle carry a single `` charge '' @xmath169 , so that @xmath170 and each @xmath169 is randomly distributed ( `` random sequence model '' ) .
it turns out to be sufficient for the purpose of this paper to restrict ourself only with the `` random - bond model '' , where @xmath168 does not depend from the choice of pairs and has a gaussian distribution @xmath171 the competing long range interactions frustrates the system of particles and the question is whether a glass transition exist or not .
normally frustration and frozen disorder is enough for the existence of glassy phases .
here the problem is more complicated , since the long range nature of the interaction may provide opposite effects .
the averaging over the quenched disorder in eq.([6 ] ) ( after the substitution @xmath172 can be carried out just in the same way as in ref.@xcite .
similarly , typical two - time dependent terms immediately appear .
they are also bilinear with respect to the forces of interaction @xmath173 . in order to rationalize these terms
it is convenient to introduce ( besides the mass density ( [ 8 ] ) and the response field density ( [ 9 ] ) ) the following collective variables @xmath174 after the introduction of the 4-dimensional column - fields @xmath175\ ; , \end{aligned}\ ] ] where @xmath176and 4@xmath634 - matrix @xmath177\label{55}\end{aligned}\ ] ] the whole expression for gf takes the form @xmath178 where the entropy is given by @xmath179\nonumber\\ & \times&\prod_{{\rm a}=1}^{4}\delta\left[q_{{\rm a}}(1 ;    2)-\sum_{p=1}^{n}\:p_{{\rm a}}^{(p)}(1 ; 2)\delta({\bf r}_{1}-{\bf    r}^{(p)}(t_{1}))\delta({\bf r}_{2}-{\bf    r}^{(p)}(t_{2 } ) ) \right]\ ;   .
\label{57}\end{aligned}\ ] ] we had used the column - operators @xmath180 and the external field , @xmath181 , conjugated to @xmath182 , has been introduced also .
the two- point collective fields ( [ 53 ] ) have a meaning of the dynamical `` overlaps '' .
it is a dynamical generalization of the parisi `` overlaps '' in a replica space @xcite .
for example @xmath183 quantify density - density and @xmath184 response - density overlaps respectively between two space - time points .
the `` entropy '' ( [ 57 ] ) corresponds to the volume in the dynamical phase space when not only fields @xmath72 but also overlaps @xmath185 are given . in a sense the `` entropy '' ( [ 57 ] ) is again the generalization of the entropy for the heteropolymer spanned in a replica space at the given set of `` overlaps '' @xcite .
let us introduce the functional @xmath186 by the functional fourier transformation @xmath187 after substitution in eq.([56 ] ) and integration over @xmath188 and @xmath189 one gets @xmath190_{\alpha\beta}(1,2)\psi_{\beta}(2)\nonumber\\ & - & \frac{n}{\chi_{0}^2}\int d1 d2 d3 d4 \phi_{\rm      a}(1;2)\left[\gamma^{-1}\right]_{{\rm ab}}(1 , 2 , 3 , 4 ) \phi_{\rm      b}(3;4)\bigg\}\ ; , \label{60}\end{aligned}\ ] ] where @xmath191\nonumber\\ - i\sum_{p=1}^{n}\int dt r_{\alpha}^{(p)}(t)\left[\psi_{\alpha}\left({\bf r}^{(p)}(t)\right ) +    i\chi_{\alpha}\left({\bf r}^{(p)}(t)\right)\right]\nonumber\\ - i\sum_{p=1}^{n}\int dt dt ' p^{(p)}_{a}(t;t')\left[\phi_{a}\left[{\bf      r}^{(p)}(t);{\bf r}^{(p)}(t')\right ] + ih_{a}\left[{\bf      r}^{(p)}(t);{\bf r}^{(p)}(t')\right]\right]\bigg\}\ ; . \label{61}\end{aligned}\ ] ] in order to ensure the extensivity of the whole effective action in eq.([60 ] ) we put the variance @xmath192 ( so that the variance of the whole strength factor in eq.([3 ] ) scaled as @xmath193 akin to ref.@xcite ) .
this enable to represent gf in a similar to eq.([19 ] ) form @xmath194\bigg\}\ ; , \label{62}\end{aligned}\ ] ] where @xmath195= \frac{1}{2\xi_{0}}\int dt \int d1 d2\psi_{\alpha}(1)\left[v^{-1}\right]_{\alpha\beta}(1 , 2)\psi_{\beta}(2)\nonumber\\     - \frac{1}{\chi_{0}^2}\int d1 d2 d3 d4 \phi_{\rm      a}(1;2)\left[\gamma^{-1}\right]_{{\rm ab}}(1 , 2 , 3 , 4 ) \phi_{\rm      b}(3;4 ) -\frac{1}{n}\log
\int \prod_{p=1}^{n}d{\bf r}^{(p)}(t)d{\bf \hat r}^{(p)}(t)\nonumber\\ \times\exp \bigg\{\sum_{p=1}^{n}a_{0}[{\bf r}^{(p)},{\bf \hat r}^{(p ) } ] - i\sum_{p=1}^{n}\int dt r_{\alpha}^{(p)}(t)\left[\psi_{\alpha}\left({\bf r}^{(p)}(t)\right ) +    i\chi_{\alpha}\left({\bf r}^{(p)}(t)\right)\right ] \nonumber\\ -i\sum_{p=1}^{n}\int dt dt ' p^{(p)}_{a}(t;t')\left[\phi_{a}\left[{\bf      r}^{(p)}(t);{\bf r}^{(p)}(t')\right ] + ih_{a}\left[{\bf      r}^{(p)}(t);{\bf r}^{(p)}(t')\right]\right ] \bigg\}\ ; .\label{63}\end{aligned}\ ] ] the resulting sp - equation reads @xmath196 where the average @xmath197 is calculated with the gf @xmath198\nonumber\\    & - & i\sum_{p=1}^{n}\int dt r_{\alpha}^{(p)}(t)\left[\overline{\psi_{\alpha}}\left({\bf r}^{(p)}(t)\right ) +    i\chi_{\alpha}\left({\bf r}^{(p)}(t)\right)\right]\nonumber\\   & -&i\sum_{p=1}^{n}\int dt dt ' p^{(p)}_{a}(t;t')\left[\overline{\phi_{a}}\left[{\bf      r}^{(p)}(t);{\bf r}^{(p)}(t')\right ] + ih_{a}\left[{\bf      r}^{(p)}(t);{\bf r}^{(p)}(t')\right]\right]\bigg\}\ ; . \label{65}\end{aligned}\ ] ] thereby we are left with the gf of a free system which experiences the external mean - fields @xmath199 and @xmath200 .      in order to calculate gf given by eq.([65 ] ) we will use the self - consistent hartree approximation ( scha ) . for this approximation
we replace the real action by an appropriate gaussian one in such a way that all terms which include more then two fields @xmath201 or / and @xmath202 are written in all possible ways as products of pairs of @xmath201 or @xmath202 coupled to self - consistent averages of the remaining fields .
the analogy between scha and sp - approximation at @xmath34 for the special case when the non - quadratic terms in the action are only the functions of the mean - squared displacement @xmath203 ^ 2\right>/n$ ] has been proven in ref.@xcite . in our case the action in eq.([65 ] ) has a more general form . in the appendix
e we show that the scha and the _ next to the saddle point approximation _ ( nspa ) merge and both become exact , if the gf with an arbitrary action can be treated by a steepest descent approach at @xmath204 .
let us make the fourier transformation of the mean - fields @xmath205 and insert it in eq.([65 ] ) .
then for eq.([65 ] ) we use the hartree - type action ( see eq.([140 ] ) . by doing
so we put for simplicity the expectation value @xmath206 .
it is easy to assure oneselves also that the `` response - response overlap '' @xmath207 ( similar to @xmath208 in ref.@xcite ) . in the curse of the derivation
we have used sp - equation ( [ sp ] ) and defined the correlator ( or the incoherent scattering function ) @xmath209 as well as the response functions @xmath210 where @xmath211 stands for the averaging with the hartree - type of action .
after collection of all terms the final result ( at @xmath212 ) reads then @xmath213 + \int d t d t ' i\hat{\bf r}^{(p)}(t ) { \bf r}^{(p)}(t)\lambda(t , t')\nonumber\\ & - & \int d t d t ' i\hat{\bf r}^{(p)}(t ) { \bf r}^{(p)}(t')\lambda(t , t ' ) +   \int d t d t ' i\hat{\bf r}^{(p)}(t)i\hat{\bf r}^{(p)}(t')\eta(t , t ' )   \bigg\}\ ; , \label{71}\end{aligned}\ ] ] where @xmath214 and @xmath215 ^ 2 \ ; .\label{73}\end{aligned}\ ] ] in eqs.([71 ] ) - ( [ 73 ] ) we have restricted ourselves to the homogeneous case @xmath216 for the correlation and response function .
the equation of motion for the one particle correlator @xmath217 and the corresponding response function @xmath218 ( which actually does not depend from the particle index @xmath2 ) can be derived from eq.([71 ] ) by using the standard techniques @xcite .
the resulting equations read @xmath219{\cal p}(t , t')\nonumber -   \int\limits_{-\infty}^{t}d\tau    \lambda ( t , \tau){\cal p}(t , t ' ) \nonumber\\ + \int\limits_{-\infty}^{t'}d\tau   \ : \eta ( t , \tau){\cal g}(t ' , \tau ) = -2 t \gamma_{0}{\cal g}(t ' , t)\label{76}\end{aligned}\ ] ] and @xmath219{\cal g}(t , t')\nonumber - \int\limits_{-\infty}^{t}d\tau    \lambda ( t , \tau){\cal g}(\tau , t ' ) \nonumber\\ = - \delta(t - t')\ ; . \label{77}\end{aligned}\ ] ] eqs.([76 ] ) - ( [ 77 ] ) should be supplemented with the initial conditions @xmath220 and @xmath221 . by making use this condition , equipartition @xmath222 , causality @xmath223 at @xmath224 as well as the condition @xmath225 one find from eq.([76 ] ) the following equation @xmath226{\cal p}(t , t )   - \int\limits_{-\infty}^{t}d\tau   \lambda ( t , \tau){\cal p}(\tau , t ' ) + \int\limits_{-\infty}^{t}d\tau   \ : \eta ( t , \tau){\cal g}(t , \tau ) = 2 t .\label{76'}\end{aligned}\ ] ]    the set of eqs.([76 ] ) - ( [ 77 ] ) have the same structure as the dyson eq.([37 ] ) .
after the matrices inversions and going to the time domain eqs.([37 ] ) ( in time - translational invariant case ) take the form @xmath227g_{01}(t , t ' ) & - & \int\limits_{-\infty}^{t}d\tau \sigma_{10}(t - \tau ) g_{01}(\tau - t ' ) = \delta(t - t')\\\label{78 } \left[m_{0}\frac{\partial^{2}}{\partial t^{2 } } + \gamma_{0}\frac{\partial}{\partial t } + \mu(0)\right]g_{00}(t , t ' ) & - & \int\limits_{-\infty}^{t}d\tau \sigma_{10}(t - \tau ) g_{00}(\tau - t ' ) \nonumber\\ & -&\int d\tau\sigma_{11}(t - \tau ) g_{10}(\tau - t ' ) = 2t\gamma_{0 } g_{10}(t   - t')\ ; , \label{79}\end{aligned}\ ] ] where @xmath228 and rpa - fourier spectrum @xmath229 eqs.([76 ] ) - ( [ 77 ] ) are turned to the dyson equations ( [ 78 ] ) - ( [ 79 ] ) provided that @xmath230 we can show @xcite that the relation @xmath231 holds , provided that the fdt is satisfied for @xmath232 .
we then have in addition @xmath233 bearing eqs.([81 ] ) in mind the eq.([82 ] ) takes in our case the form ( at @xmath234 ) @xmath235 the validity of the relationship ( [ 84 ] ) can be checked by replacing ( [ 72 ] ) and ( [ 73 ] ) in ( [ 84 ] ) .
the general eqs.([76 ] ) - ( [ 77 ] ) are equivalent , _ mutatis mutandis _ , to the corresponding equations for the @xmath2 - spin system or a particle in the random potential at the large dimension @xcite
. the most important features of these equations are the _ glassy dynamical behavior _ and the universal _ aging regime_. at low temperatures the system tries to minimize the energy and each particle ( with a number @xmath2 ) tends to surround itself with other particles which assure the strength parameter @xmath236 . on the other hand
the long - range interaction tries to support other pairs @xmath237 corresponding to @xmath238 . as a result
the system becomes `` frustrated '' and many local free energy minima appear .    in the spirit of ref.@xcite when @xmath239 we have to discriminate between different cases : ( i ) the _ asymptotic regime _
when @xmath240 and ( ii ) the _ aging regime _
when @xmath241 .
the aging regime is much more complicated because the time - translational invariance and fdt are violated .
this regime has been extensively investigated both theoretically @xcite and by the computer simulation @xcite . in the following we
restrict ourselves only to the asymptotic regime , for the sake of clearness and simplicity , and since the main features will be already visible .
this asymptotic regime is characterized by the large time scales , i.e. , @xmath242 but keeping the difference @xmath243 finite . under these circumstances
we can define @xmath244 then the equation for the displacement @xmath245 $ ] , response function @xmath246 and the static correlator @xmath247 takes correspondingly the forms @xmath248 { \cal d}_{\rm as}(\tau ) - \int\limits_{0}^{\tau}d\tau ' \lambda_{\rm as } ( \tau - \tau'){\cal d}_{\rm as}(\tau ' ) - \nonumber\\ -   \int\limits_{0}^{\infty}d\tau'\left [ \lambda_{\rm as } ( \tau + \tau ' ) - \lambda_{\rm as}(\tau')\right]{\cal d}_{\rm as}(\tau ' ) - 2 \int\limits_{0}^{\infty}d\tau'\left [ \eta_{\rm as } ( \tau + \tau ' ) - \eta_{\rm as}(\tau')\right]{\cal g}_{\rm as}(\tau')= 2 t \label{86}\end{aligned}\ ] ] @xmath248 { \cal g}_{\rm as}(\tau ) - \int\limits_{0}^{\tau}d\tau ' \lambda_{\rm as } ( \tau - \tau'){\cal g}_{\rm as}(\tau ' ) = 0\label{87}\end{aligned}\ ] ] @xmath249 , \label{88 } \end{aligned}\ ] ] where @xmath250 @xmath251 however , it is also convenient to define the `` anomaly '' @xmath252 @xcite .
the eqs.([86 ] ) - ( [ 88 ] ) has been analyzed first in the context of polymeric manifold in the random media @xcite and the random - phase sine - gordon model @xcite .
the peculiarity of our model is defined by its memory functions @xmath253 and @xmath254 .
for example , let as give an explicit expression for @xmath254 .
the gaussian form of the correlator , @xmath255 , leads from eq.([73])to the result @xmath256 usually it is assumed that at the high temperature fdt holds , i.e. @xmath257 and @xmath258 in this case eqs.([86 ] ) and ( [ 87 ] ) merge and take a simple form @xmath248 { \cal d}_{\rm as}(\tau ) - \int\limits_{0}^{\tau}d\tau ' \lambda_{\rm as } ( \tau - \tau'){\cal d}_{\rm as}(\tau ' ) = 2 t \ ; .\label{94}\end{aligned}\ ] ] it turns out @xcite that the solution which satisfies the fdt is only stable above a critical temperature @xmath259 . for the stability analysis
it is convenient to represent eq.([94 ] ) in the form @xmath260 { \cal d}_{\rm as}(\tau ) -\int\limits_{0}^{\tau}d\tau'\left[\eta_{\rm as } ( \tau - \tau ' ) - \eta_{\rm as}(\tau)\right]\frac{\partial}{\partial\tau'}{\cal d}_{\rm as}(\tau ' ) = 2t\ ; , \label{95}\end{aligned}\ ] ] where @xmath261 for @xmath262 the stability condition which comes out of eq.([95 ] ) reads @xmath263{\cal    d}_{\rm as}(\tau ) \le 2 t \quad.\label{97}\end{aligned}\ ] ] then the stationary value of the displacement @xmath264 reads    @xmath265    by taking into account eqs.([92 ] ) and ( [ 93 ] ) the stability condition becomes @xmath266 for @xmath267 , where @xmath268 q - \left(\frac{\chi_{0}}{t}\right)^2 \frac{\sqrt{\pi}}{12 } \sqrt{q } + 1 \ ; .\label{100}\end{aligned}\ ] ] the critical values @xmath269 and @xmath259 at which the condition ( [ 99 ] ) first becomes violated is defined by equations @xmath270 consequently , eqs .
( [ 102 ] ) have the simple solution @xmath271    the fig.3 shows the behavior of @xmath272 in the vicinity of the critical point .
it can be seen that the minimum , @xmath273 , at which @xmath274 appears continuously , i.e. the instability of the fdt solution shows up as a @xmath275 order phase transition .
this is analogous to the dynamics of polymeric manifolds in a medium with the long range correlation in disorder @xcite . in particular , if `` anomaly '' @xmath276 then @xmath277 and @xmath278 , so in this case the fdt solution is unstable for any finite temperature .
+    let us consider the dynamics at the temperature slightly above the critical point : @xmath279 , where @xmath280 . for large @xmath281
the decomposition @xmath282 is possible , where @xmath283 .
the substitution of this decomposition into eq.([95 ] ) and the expansion up to the second order with respect to @xmath284 yields @xmath285 ^ 2 + \frac{1}{2}\int\limits_{0}^{\tau}d\tau'\left[f(\tau - \tau ' ) -    f(\tau)\right]\frac{\partial}{\partial\tau'}f(\tau ' ) = 0 \ ; .\label{105}\end{aligned}\ ] ] following ref.@xcite let us make the laplace transformation @xmath286 \equiv { \tilde f}(z)$ ] and introduce the scaling functions , @xmath287 or @xmath288 , in a such way @xmath289 where @xmath290 and @xmath291 . if @xmath292 and @xmath293 then one can write eq.([105 ] ) in the form @xmath294 ( see eq.(2.68b ) from ref.@xcite ) .    in the critical regime @xmath295 ( or @xmath296 )
the solution of eq.([108 ] ) has a form @xmath297 . in this limit the first term in eq.([108 ] )
is dropped out and the exponent is defined by the equation @xmath298 the solution of eq .
( [ 109 ] ) gives @xmath299 . in the opposite limit @xmath300 ( or @xmath301 ) the last term in eq.([108 ] ) can be neglected . in this case
the solution has the form @xmath302 , where @xmath303 . as a result
the overall scaling reads @xmath304 where @xmath305 is some constant .    at @xmath306
the fdt is violated for the large time separation @xmath281 and the aging regime is arising .
it should be mentioned that the asymptotic regime can not be decoupled from the aging one @xcite .
in actual fact , the `` anomaly '' @xmath307 in the asymptotic eqs .
( [ 86 ] ) - ( [ 88 ] ) strictly speaking can be calculated only from the aging regime . because of the distinct aim in this paper we are not going to discuss the aging regime here expecting to return to it in a later publication .
in the present paper we have considered the dynamics of two models with the long range repulsive interaction .
the interaction potential was designed in a way to enable the saddle point treatment , as well as a fluctuation expansion .    for the pure model we have derived eq.([48 ] ) for the full correlation matrix @xmath308 in the one - loop approximation , which has an the explicit solution ( see eqs.([49 ] ) - ( [ 51 ] )  ) .
this solution has a `` boring behavior '' at @xmath162 which manifests the absence of the glass dynamics .
the physical background of this stems from the fact that the potential is much too soft and the `` cage effect '' is completely missing .
this conclusion is in accordance with the interacting particles statistical thermodynamics analysis , which was given in ref.@xcite .
it was shown there that for the infinite range interaction potential , which allows a well defined saddle point treatment , the glassy phase is simply suppressed .
on the other hand , the same model but with a randomly distributed strength of interaction ( the `` random - bond model '' ) leads to the continuous glass transition .
this type of transition is also the case for the polymeric manifolds in the disordered medium with long range correlation @xcite as well as for the @xmath2 - spin interaction spin - glass model at the large external field @xcite
. it would be also interesting to investigate the more realistic `` random sequence model '' when each particle carry a random `` charge '' .
qualitatively the same glassy behavior have been found in the pure spin models with the deterministic but very rapidly oscillating coupling between variables @xcite .
it was assumed that the effective quenched disorder is `` self - induced''@xcite .
this means that because of the slow dynamics some degrees of freedom freeze and play the role of the effectively quenched disorder .    as a conclusion ,
the glass transition in the pure systems of the interacting particles , where the disorder is actually `` self - induced '' , goes beyond the mean - field level @xcite .
this appears to hard to implement in the present context , because it implies the consideration of the short range interaction potential as well as activated processes .
the authors have benefited from discussions with m. fuchs , a. latz , g. migliorini and greatly indebted to the deutsche forschungsgemeinschaft ( dfg ) , the sonderforschungsbereich sfb 262 for financial support of the work .
in the full analogy with eq.([25 ] ) the expression for @xmath309 reads @xmath310 \label{113}\end{aligned}\ ] ] the expansion of the @xmath311 up to the @xmath312 order with respect to the mean field @xmath313 can be easily obtained from eq.([24 ] ) @xmath314\nonumber\\ & + & \frac{1}{2 ! } \int d2 d3\left<\delta\rho_{\alpha}(1)\delta\rho_{\beta}(2)\delta\rho_{\gamma}(3 ) \right>_{0 } \left[\chi_{\beta}(2 )   -    i\bar{\psi_{\beta}}(2)\right ] \left[\chi_{\gamma}(3 )   -    i\bar{\psi_{\gamma}}(3)\right]\nonumber\\ & + & \frac{1}{3!}\int d2 d3 d4\left<\delta\rho_{\alpha}(1)\delta\rho_{\beta}(2)\delta\rho_{\gamma}(3 ) \delta\rho_{\delta}(4 ) \right>_{0 } \nonumber\\ & \times&\left[\chi_{\beta}(2 )   -    i\bar{\psi_{\beta}}(2)\right ] \left[\chi_{\gamma}(3 )   -    i\bar{\psi_{\gamma}}(3)\right]\left[\chi_{\delta}(4 )   -    i\bar{\psi_{\delta}}(4)\right]\quad . \label{114 } \end{aligned}\ ] ] by using the sp - equation ( [ 23 ] ) and after threefold differentiation with respect to @xmath65 ( see eq.([113 ] ) ) we find @xmath315^{-1 } \right\}_{{\bar\alpha}\alpha}({\bar 1 } , 1)\left\{\left[{\hat 1 } + \mu { \hat v } { \hat        f}^{(2)}\right]^{-1 } \right\}_{{\bar\beta}\beta}({\bar 2 } , 2)\nonumber\\ & \times&\left\{\left[{\hat 1 } + \mu { \hat v } { \hat        f}^{(2)}\right]^{-1 } \right\}_{{\bar\gamma}\gamma}({\bar 3 } , 3)\left\{\left[{\hat 1 } + \mu { \hat v } { \hat        f}^{(2)}\right]^{-1 } \right\}_{{\bar\delta}\delta}({\bar 4 } , 4)\quad , \label{115}\end{aligned}\ ] ] where the 4-point free system correlation matrix @xmath316 in eq.([115 ] ) we imply the summation ( integration ) over the barred indices ( barred space - time variables ) . when deriving eq.([115 ] )
we have also kept in mind that @xmath317 and @xmath318 , etc .
the fact that the matrix @xmath319 is symmetrical with respect to simultaneous permutations of greek indices and space - time arguments as well as eq.([26 ] ) have been used .
it is easy to show that @xmath319 is factorized @xmath320 on the other side it is instructive to check that even in this case @xmath309 can not be factorized .
the substitution of the eg.([115 ] ) into eq.([46 ] ) after the straightforward algebra yields @xmath322^{-1 } - \left [ \mu^{-1 } { \hat v}^{-1 } +      { \hat f}\right]^{-1}\right\}_{{\bar\beta}{\bar\alpha}}({\bar 2 } , { \bar 1})\nonumber\\ & \times&\left\{f_{{\bar\alpha}{\bar\beta}}({\bar 1},{\bar      2})f_{{\bar\gamma}{\bar\delta}}({\bar 3},{\bar 4 } ) + f_{{\bar\alpha}{\bar\gamma}}({\bar 1},{\bar      3})f_{{\bar\beta}{\bar\delta}}({\bar 2},{\bar 4 } ) + f_{{\bar\alpha}{\bar\delta}}({\bar 1},{\bar      4})f_{{\bar\beta}{\bar\gamma}}({\bar 2},{\bar      3})\right\}\nonumber\\ & \times&\left\{\left[{\bar 1 } + 2\mu { \bar v } { \bar        f}\right]^{-1}\right\}_{{\bar\gamma}\alpha}({\bar 3 } , 1)\left\{\left[{\bar 1 } + 2\mu { \hat v } { \hat        f}\right]^{-1}\right\}_{{\bar\delta}\beta}({\bar 4 } , 2)\quad,\end{aligned}\ ] ] where as before for the repeated barred indices ( variables ) the summation ( integration ) is implied . for the time - space - translational invariant case the respective fourier transformation leads to the result : @xmath323^{-1 }
- \left[\mu^{-1}{\hat v}^{-1 } +      { \hat f}\right]^{-1}\right\}_{{\bar\beta}{\bar\alpha}}(-{\bf k } , -\omega )   f_{{\bar\alpha}{\bar\gamma}}(-{\bf     k},-\omega)f_{{\bar\beta}{\bar\delta}}({\bf k},\omega)\nonumber\\ & + & \left\{\left[\left(2 \mu \right)^{-1}{\hat v}^{-1 } +      { \hat f}\right]^{-1 } - \left[\mu^{-1}{\hat v}^{-1 } +      { \hat f}\right]^{-1}\right\}_{{\bar\beta}{\bar\alpha}}({\bf k } , \omega )   f_{{\bar\alpha}{\bar\delta}}({\bf     k},\omega)f_{{\bar\beta}{\bar\gamma}}(-{\bf k},-\omega)\bigr\}\nonumber\\ & \times&\left\{\left[{\bar 1 } + 2\mu { \hat v } { \hat        f}\right]^{-1}\right\}_{{\bar\gamma}\alpha}(-{\bf k } , -\omega)\left\{\left[{\bar 1 } + 2\mu { \bar v } { \bar        f}\right]^{-1}\right\}_{{\bar\delta}\beta}({\bf k } , \omega)\quad,\end{aligned}\ ] ] where the trace @xmath324^{-1 } -    \left[{\hat 1 } + \mu^{-1}{\hat f}^{-1}{\hat v}^{-1}\right]^{-1 } \right\}_{{\bar\alpha}{\bar\alpha}}({\bf q } , \omega)\quad.\end{aligned}\ ] ] with the correlation matrix @xmath325 given by eqs.([49 ] ) , by doing integration over @xmath326 one can check that the trace @xmath327 .
this gives finally @xmath328 where @xmath329^{-1}\right\}_{\alpha{\bar\gamma}}({\bf k } , \omega)f_{{\bar\gamma}{\bar\beta}}({\bf k},\omega)\nonumber\\ & \times&\left\{\left[{\hat 1 } + \left (        2 \mu \right)^{-1}{\hat f}^{-1}{\hat v}^{-1}\right]^{-1 } -    \left[{\hat 1 } + \mu^{-1}{\hat f}^{-1}{\hat v}^{-1}\right]^{-1 } \right\}_{{\bar\beta}{\bar\delta}}({\bf q } , \omega)\nonumber\\ & \times&\left\{\left[{\bar 1 } + 2\mu { \hat v } { \hat        f}\right]^{-1}\right\}_{\alpha{\bar\gamma}}({\bf k } , \omega)\quad .\end{aligned}\ ] ]
in this case the direct correlation function @xmath330 and its fourier transformation takes the scaling form @xmath331 let us insert this expression to the mct - memory kernel ( see eq.(3.32 ) in @xcite ) .
it is reasonable then to rescale the integration variables in the memory kernel , @xmath332 as well as to put for the external wave vector @xmath333 , where @xmath334 is some reference wave vector .
the last scaling means that in the mf - limit an experiment probes a very small wave vector : @xmath335 .
the resulting scaling of the memory kernel , @xmath336 , reads @xmath337\right\}^2}{q_{0}^2}\ :
\widetilde{s}({\bf k } , t ) \:\widetilde{s}({\bf p } , t)\ ; , \label{kernel}\end{aligned}\ ] ] where we have took into account the scaling form of the correlator : @xmath338 .
thus we finally arrive at the scaling @xmath339 and the glass transition dies out .
let us prove that scha becomes exact for gf given by ( [ 65 ] ) in the limit @xmath34 .
we will consider even a more general gf @xmath340 + \sum_{p=1}^{n}\int   d1 x_{\alpha}^{(p)}(1)\chi_{\alpha}(1)\bigg\ } \;,\label{117}\end{aligned}\ ] ] where we have used the shorthand notations @xmath341 and `` 1 '' embraces cartesian indices as well as time variable : @xmath342 . in eq.([117 ] ) @xmath343 $ ] is an arbitrary functional of @xmath344 .    instead of the exact action functional in eq.([117 ] )
we consider now the trial one which has a gaussian form @xmath345 & = & \sum_{p=1}^{n}\bigg\ { \frac{1}{2 } \int d1 d2 \:x_{\alpha}^{(p)}(1)a_{\alpha\beta}(1 , 2)x_{\beta}^{(p)}(2)\nonumber\\ & -&\int d1 d2\:x_{\alpha}^{(p)}(1){\gamma}_{\alpha\beta}(1 , 2)x_{\beta}^{(p)}(2 ) - \int d1 l_{\alpha}(1)x_{\alpha}^{(p)}(1)\bigg\ }   \ ; .\label{119}\end{aligned}\ ] ] let us look for the `` best '' coefficients @xmath346 and @xmath347 in a sense that the exact `` free energy '' @xmath348 = - \log z\left\{\chi_{\alpha}\right\}$ ] tends to the trial one @xmath349 = - \log \int\prod d x_{\alpha}^{(p)}\exp \ { - s [ x_{\alpha}^{(p)}]\}$ ] , i. e. @xmath350 \longrightarrow f_{0}\left[\chi_{\alpha}\right]\label{120}\end{aligned}\ ] ] and both becomes exact at @xmath34 .
we can show that the property ( [ 119 ] ) is satisfied by @xmath351 and @xmath352 which are obtained by extremization of the functional @xmath353\bigg\ } \nonumber\\ & + & \sum_{p=1}^{n}\bigg\{\int d1 d2 \gamma_{\alpha\beta}(1 , 2 ) \left < x_{\alpha}^{(p)}(1)x_{\beta}^{(p)}(2)\right>_{\rm s}\nonumber\\   & + & \int d1
\left[l_{\alpha}(1 ) - \chi_{\alpha } \right ] \left < x_{\alpha}^{(p)}(1)\right>_{\rm s } - \left < w\left[x_{\alpha}^{(p)}(1)\right]\right>_{\rm s } \bigg\}\ ; , \label{121}\end{aligned}\ ] ] where we use the notations @xmath354\bigr\}}{\int\prod_{p=1}^{n}\prod_{\alpha=0 ,    1}dx_{\alpha}^{(p)}(1)\:\exp\bigl\ { - s\left [    x_{\alpha}^{(p)}\right]\bigr\ } } \;.\label{122}\end{aligned}\ ] ] the extremization conditions reads @xmath355 the variations in eqs.([122 ] ) can be done directly . during the calculation
the generalized wick s theorem @xcite should be also taken into account .
namely , because the averaging ( [ 121 ] ) is simply the gaussian integral the wick s theorem yields @xmath356\right>_{\rm s }   & = & \left < x_{\alpha}^{(p)}(1)\right>_{\rm s}\left < w\left[x_{\alpha}^{(p)}\right]\right>_{\rm s } \nonumber\\ & + & \int d2 \:\left<\delta x_{\alpha}^{(p)}(1)\delta x_{\beta}^{(p)}(2 ) \right>_{\rm s}\:\left<\frac{\delta}{\delta x_{\beta}^{(p)}(2 ) }   \ : w\left[x_{\alpha}^{(p)}\right]\right>_{\rm s}\ ; , \label{124}\end{aligned}\ ] ] where @xmath357 .
after the straightforward calculation we find @xmath358\right>_{\rm s}\label{125}\end{aligned}\ ] ] and @xmath359\:\left <    x_{\alpha}^{(p)}(2)\right>_{\rm s}\ ; .\label{126}\end{aligned}\ ] ] then equations for the two moments take the form @xmath360_{\alpha\beta]}(1 , 2 ) \left[\left<\frac{\delta}{\delta          x_{\beta}^{(p)}(2 ) } w\left[x_{\alpha}^{(p)}\right]\right>_{\rm        s } + \chi_{\beta}(2)\right]\label{127}\end{aligned}\ ] ] and @xmath361^{-1}\right\}_{\alpha\beta}(1 , 2 ) \ ; , \label{128}\end{aligned}\ ] ] where @xmath362 and @xmath363 stands for the corresponding @xmath364 - matrices .    on the other side
the saddle point ( sp ) treatment of eq.([117 ] ) at @xmath365 yields @xmath366 and @xmath367^{-1}\right\}_{\alpha\beta}(1 , 2 ) \ ; , \label{130}\end{aligned}\ ] ] where @xmath368 and @xmath369 stands for the field in sp .
in order to show the analogy between eqs.([127])-([128 ] ) and eqs.([129])-([130 ] ) let us make the functional fourier transformation @xmath370\right\ } = \int d x_{\alpha}^{(p)}(1)\exp\left\{w\left[x_{\alpha}^{(p)}(1)\right ] - i\int d1\:x_{\alpha}^{(p)}(1)y_{\alpha}^{(p)}(1)\right\}\label{132}\end{aligned}\ ] ] and its inversion @xmath371\right\ } = \int d y_{\alpha}^{(p)}(1)\exp\left\{k\left[y_{\alpha}^{(p)}(1)\right ] + i\int d1\:x_{\alpha}^{(p)}(1)y_{\alpha}^{(p)}(1)\right\}\ ; .\label{133}\end{aligned}\ ] ] then eqs.([129 ] ) - ( [ 130 ] ) can be written as @xmath372_{\alpha\beta}(1 , 2)\left[i\left < y_{\beta}^{(p)}(2)\right>_{\rm sp } + \chi_{\beta}(2)\right]\label{134}\end{aligned}\ ] ] and @xmath373^{-1}\right\}_{\alpha\beta}(1 , 2)\ ; , \label{135}\end{aligned}\ ] ] where the correlation matrix @xmath374 and @xmath375 + i\int d1\:\bar{x}_{\alpha}^{(p)}(1)y_{\alpha}^{(p)}(1)\right\}}{\int d y_{\alpha}^{(p ) }    \exp\left\{k\left[y_{\alpha}^{(p)}\right ] + i\int d1\:\bar{x}_{\alpha}^{(p)}(1)y_{\alpha}^{(p)}(1)\right\}}\ ; .\label{137}\end{aligned}\ ] ] at @xmath376 by making use eq.([133 ] ) one can immediately see that @xmath377\right>_{\rm         s } \longrightarrow
i\left < y_{\alpha}^{(p)}(1)\right>_{\rm sp}\label{138}\end{aligned}\ ] ] and @xmath378\right>_{\rm s } \longrightarrow   - \left<\delta   y_{\alpha}^{(p)}(1)\delta   y_{\beta}^{(p)}(2 ) \right>_{\rm sp}\ ; , \label{139}\end{aligned}\ ] ] and scha exactly corresponds to nspa . for the case which was treated in sec.iii c @xmath379 and the hartree - type action ( [ 119 ] ) cast the form @xmath345 & = & \sum_{p=1}^{n}\bigg\ { \frac{1}{2 } \int d1 d2 \:x_{\alpha}^{(p)}(1)a_{\alpha\beta}(1 , 2)x_{\beta}^{(p)}(2)\nonumber\\ & -&\frac{1}{2}\int d1 d2\:\left<\frac{\delta^2}{\delta x_{\alpha}^{(p)}(1 ) \delta x_{\beta}^{(p)}(2 ) }   \ : w\left[x_{\alpha}^{(p)}\right]\right>_{\rm s }   x_{\alpha}^{(p)}(1 )   x_{\beta}^{(p)}(2 ) \bigg\ }   \ ; .\label{140}\end{aligned}\ ] ]    99 w.gtze in , _ liquids , freezing and glass transition _ , ed . by j.p.hansen , d.levesque and j.zinn-justin , amsterdam , north - holland , 1991 v.n .
prigodin , j. phys .
* c 4 * , 785 ( 1992 ) t.r .
kirkpatrick and d. thirumalai , phys . rev . * b36 * , 5388 ( 1987 ) t.r .
kirkpatrick and d. thirumalai , phys .
* b37 * , 5342 ( 1988 ) t.r .
kirkpatrick and d. thirumalai , trans .
theor . and
stat . phys .
* 24*,927 ( 1995 ) a. crisanti , h. horner and h .- j .
sommers , z. phys . *
92 * , 257 ( 1993 ) j. - p. bouchaud , l.f .
cugliondolo , j. kurchan and m. mzard , physica * a226 * , 243 ( 1996 ) j. - p. bouchaud , l.f .
cugliondolo , j. kurchan and m. mzard , in , _ spin glasses and random fields _ , ed . by a.p .
young , world scientific , 1997 h. kinzelbach and h. horner j. phys.i ( france ) * 3 * , 1329 ( 1993 ) h. kinzelbach and h. horner j. phys.i ( france ) * 3 * , 1901 ( 1993 ) s. franz and m. mzard , physica * a210 * , 48 ( 1993 ) v.g .
rostiashvili and t.a .
vilgis , europhys . lett . * 49 * , 162 ( 2000 ) r.g .
palmer , adv .
phys . * 31 * , 669 ( 1982 ) l.f .
cugliondolo , j. kurchan , g. parisi and f. retort , phys .
rev . lett.*74 * 1012 ( 1995 ) j.l . lebowitz and o. penrose , j. math . phys.*7 * , 98 ( 1996);o .
penrose and j.l .
lebowitz , j. stat .
phys.*3 * , 211 ( 1971 ) j.- p. hansen and i.r .
mcdonald , _ theory of simple liquids _ academic press , london , 1986 c. dedominicis and l. pelityphys . rev . * b18*,353 ( 1978 ) ; r. bauch , h.k .
janssen and h.wagner .
z. phys.*b24 * , 113 ( 1976 ) v.g .
rostiashvili and r. schilling z. phys.*b102 * , 117 ( 1997 ) v.g .
rostiashvili , m. rehkopf and t.a .
vilgis , eur .
phys . j.*7 * , 233 ( 1998 ) v.g .
rostiashvili , m. rehkopf and t.a .
vilgis , j .
110 * , 639 ( 1999 ) v.g .
rostiashvili , physica*a231 * , 575 ( 1996 ) j. zinn - justin , _ quantum field theory and critical phenomena _
clarendon press , oxford , 1989 d. forster , _ hydrodynamic fluctuations , brocen symmetry and correlation function _ , benjamin , 1975 a. weyersberg and t.a .
vilgis , phys .
e , * 48 * , 277 , 1993 , t.a .
vilgis , a. weyersberg , m.g .
brereton phys .
e * 49 * , 3031 , 1994 m. fuchs , private communication a.i .
melcuk , r.a .
ramos , h. gould , w. klein and r.d .
mountain phys .
rev . lett.*75 * , 2522 ( 1995 ) ; w. klein , h. gould , r.a .
ramos , i. clejan and a.i .
melcuk physica * a205 * 738 ( 1994 ) t. garel , h. orland and e. pitardin , in _ spin glasses and random fields _
young , world scientific , 1997 t. garel , h. orland , europhys . lett .
* 7 * , 307 ( 1988 ) e.i .
shakhnovich and a.m. gutin , j. phys.*a22 * , 1647 ( 1989 ) c. d. sfatos , a.m .
gutin and e.i .
shakhnovich phys .
* e48*,465 ( 1993 ) h. sompolinsky and a. zippelius , phys .
* b25*,6860 ( 1982 ) h.g .
schuster , z .
phys.*b45*,99 ( 1981 ) m. mzard , g. parisi and m.a .
virazoro , _ spin glass and beyond _ , singapore , world scientific , 1987 l.f.cugliondolo and j. kurchan phys .
rev . lett.*71 * , 173 ( 1993 ) l.f.cugliondolo and p. le doussal phys . rev.*e53
* , 1525 ( 1996 ) g. parisi phys .
lett . * 79 * , 3660 ( 1997 ) w. kob and j .- l
. barrat phys .
lett . * 78 * , 4581 ( 1997 ) ; physica * a263 * , 234 ( 1999 ) ; cond - mat/9905248 d. cule and y. shapir phys .
rev.*e53 * , 1553 ( 1996 ) j.p .
bouchaud and m. mzard j. phys .
i ( france ) * 4 * 1109 ( 1994 ) g. migliorini and f. ritort j. phys .
* a27 * 7669 ( 1994 ) w. krauth and m. mzard z. phys.*b97 * 127 ( 1995 ) s. franz and j. hertz phys .
74 * , 2114 ( 1995 )