in the last few years , great progress has been made towards the understanding of the out of equilibrium dynamics of _ classical _ systems with slow dynamics and aging effects .
experiments in several glassy materials such as polymer glasses,@xcite spin - glasses,@xcite orientational glasses,@xcite simple liquids like glycerol@xcite and gels , @xcite have shown that all these systems share an extremely slow dynamics at low temperatures .
although some aspects of the dynamic evolution as the precise scaling laws or the dependence on cooling rates vary from system to system , all the examples mentioned above are characterized by the existence of a non - equilibrium low temperature phase with _
aging effects_. this means that measurements show a strong dependence upon the time elapsed since its preparation by a thermal quench or annealing from a high temperature phase .
needless to say , the standard equilibrium approach is not suited to describe such a phenomenology .    in order to describe the spin - glass case ,
a number of theoretical methods have been exploited , ranging from scaling arguments based on coarsening ideas;@xcite phenomenological models founded on particular assumptions on the structure of phase space;@xcite analytical solutions to models with infinite - range interactions;@xcite and numerical simulations of finite and infinite dimensional models.@xcite all of these approaches have succeeded in capturing some of the effects observed experimentally .
importantly enough , the solution of simplified models with infinite range interactions has been fruitful for several reasons :
( i ) it has shown that indeed a large family of such models do capture aging effects that resemble , at least qualitatively , the experimental observations .
( ii ) it has provided a general framework to analyze the non - equilibrium dynamics with new predictions , like scaling laws for the time - delayed correlation function and asymptotic violations of the fluctuation - dissipation theorem , that are not restricted to the infinite dimensional case . some of these have been confirmed by numerical studies of more realistic models@xcite and will be probed experimentally soon.@xcite ( iii ) it has allowed to extend the _
connection@xcite between mean - field theories of disordered systems and the mode - coupling theory of super - cooled liquids@xcite to the low temperature glassy phase.@xcite this is most welcome since many of the similarities observed in the behavior of real glasses of very different origin can now be understood theoretically .
( iv ) problems of a more general nature such as the dynamics of manifolds in a random quenched environment , relevant to the description of e.g. the pinning of the flux lattice by disorder in dirty high @xmath1 superconductors,@xcite can be studied in a similar way.@xcite the dependence on the internal dimension of the manifold can be included in the formalism and the analysis yields , among other results , dynamic scaling laws for the displacements with aging effects .
the studies we have mentioned so far are mostly concerned with the analysis of _ classical _ systems .
this treatment is in general justified since , in many cases , the critical temperature at which the transition to the  glassy phase " occurs is sufficiently high as to make quantum effects irrelevant . nevertheless , in many cases of practical interest
, the critical temperature can be lowered by tuning an external parameter .
experimental realizations of this situation are given by insulating magnetic materials@xcite such as liho@xmath2y@xmath3f@xmath4 and randomly mixed hydrogen - bonded ferro - antiferroelectric crystals of the rb@xmath3 ( nh@xmath4)@xmath5h@xmath6po@xmath4 type ( known as radp).@xcite    theoretical studies addressing the effects of quantum fluctuations in disordered media have been considered over the last twenty years .
a variety of techniques including replica theory,@xcite renormalization group@xcite and monte carlo simulations@xcite have been employed in order to understand the _
static _ properties and critical behavior of quantum systems such as quantum rotors@xcite , transverse field ising models in finite dimensions@xcite and mean - field systems like the sherrington - kirkpatrick model in a transverse field.@xcite    from the study of _ classical _ glassy systems described in the first paragraph , one knows that at low temperatures their dynamics becomes so slow that the systems are not able to reach thermal equilibrium at any experimentally accessible time .
the same occurs regarding numerical simulations of large systems .
it is then plausible that even at very low temperatures , say @xmath7 , with only quantum fluctuations driving the dynamics , the materials and/or models will need an unaffordable long time to reach equilibrium and that the relevant evolution will be characterized by non - equilibrium effects .
indeed , it is well - known that certain glasses in the limit of very low temperatures have aging effects@xcite and hence the approaches based on the assumption of equilibration have a restricted domain of application . instead
, one should necessarily start from a quantum _ dynamic _ description in order to obtain sensible information about the systems .
to the best of our knowledge this problem has not been theoretically tackled , from a microscopic point of view , yet .
the question then arises as to how do quantum fluctuations modify the _ real time _
dynamics of glassy systems . in this paper
, we treat the case of a disordered model in contact with a thermal quantum bath .
our first aim is to present a general framework that allows us to study the real time dynamics of a system with quantum _ and _ thermal fluctuations starting from an arbitrary initial condition . in the case of classical models
this is most easily done within the langevin approach .
one generally starts by postulating that the dynamics of the microscopic variables is given by a stochastic langevin equation .
the noise and friction terms account for the bath - system interaction . a path - integral formulation of the generating functional , known as the martin - siggia - rose formalism@xcite ( msr ) , yields an elegant and simple formulation which is particularly well suited for problems with disorder . however , the generalization of this procedure to the quantum case is far from being straightforward .
the langevin equation does not seem to be a good starting point since , being phenomenological , it can not be quantized in an obvious way.@xcite in other words , it is not clear how to make the langevin approach compatible with the intrinsic microscopic dynamics of a quantum system .
the approach we use here results from the combination of different techniques which have been intensively discussed in the literature and successfully applied to a wide range of phenomena .
essentially , it consists of the schwinger and keldysh@xcite closed - time path - integral formalism applied to a system in contact with a thermal bath , combined with the integration of the bath variables leading to the feynman and vernon influence functional.@xcite a subsequent integration over the disorder gives the quantum analogue of the msr generating functional for disordered systems .
though the formalism allows us to consider general cooling procedures or / and various temperature variations , we shall not study involved thermal histories here but concentrate instead on a quench from a configuration associated to the disordered phase .    for the sake of definiteness we focus on a system
whose classical counterpart , the so - called `` @xmath0-spin spherical spin - glass model '' , has been studied in great detail .
this model is described by a potential energy @xmath8 = - \sum_{i_1 < \dots
< i_p}^n   j_{i_i\dots i_p } \ , \sigma_{i_1 } \dots \sigma_{i_p}\ ] ] with @xmath9 quenched independent gaussian variables and @xmath10 a parameter .
let us attempt a short review of the behavior of the _ classical _ @xmath0-spin models .    in its ising version
it was first introduced by derrida in relation with the random energy model.@xcite several studies of its static properties , in particular using replicas , followed .
@xcite kirkpatrick and thirumalai@xcite studied the equilibrium dynamics ,  la sompolinsky,@xcite and remarked that the dynamic equations , in the high temperature phase , are those appearing as @xmath11 models in the mode - coupling theory of super - cooled liquids developed to describe the glass transition.@xcite    the spherical version , where a time - dependent lagrange multiplier is introduced to enforce the spherical constraint @xmath12 , was first studied by jones , kosterlitz and thouless@xcite for @xmath13 .
its statics is fully solvable either by using replicas  a replica symmetric ansatz is exact  or through a direct calculation .
its dynamics@xcite is also fully solvable and it renders explicit the fact that the model is formally related@xcite to the @xmath14 model in the large @xmath15 limit and dimension @xmath16 . in the classical limit , it is then a toy model for ferromagnetic domain - growth . in the quantum case
one feels tempted to relate it to coarsening in quantum systems .
the analysis of the spherical model for @xmath17 is due to crisanti and sommers,@xcite who also studied the equilibrium dynamics  la sompolinsky in collaboration with horner.@xcite in the classical limit and at the static level , this model is the prototype spin - glass solved by a parisi one - step replica symmetric ansatz .
still at the classical level and in the large @xmath15 limit , one derives exact evolution equations that exhibit a first order  discontinuous  dynamic transition at @xmath18 . above @xmath18 , as in the ising case , one recovers the mode - coupling equation for model @xmath11.@xcite below @xmath18 it is well - known that the dynamic approach used in ref.[@xcite ] is incorrect.@xcite the dynamics in the low - temperature phase was solved in ref .
[ ] showing that the dynamics of mean - field spin - glasses strikingly resembles that of real systems .
in the spin - glass phase the evolution is intrinsically out of equilibrium and the systems age in similar way to the one observed in experiments .    in order to define the model at the quantum level
, we need to specify the dynamics of the variables @xmath19 . for simplicity
, we consider a kinetic term of the form @xmath20 where the dot denotes a time - derivative .
calling @xmath21 the lagrange multiplier enforcing the spherical constraint , the model we here consider is then described by the action @xmath22 =   \int dt \ ,   \left [ \frac{m}{2 } { \dot { \bbox \sigma}}^2 -   \frac{z}{2 } \left ( { \bbox \sigma}^2 -n \right )   + \sum_{i_1 < \dots < i_p}^n   j_{i_i\dots i_p } \ , \sigma_{i_1 } \dots \sigma_{i_p }   \right ]   \ ; .\ ] ]    our main goals are to detect a dynamic transition from a disordered phase to an ordered glassy phase and to characterize the dynamics in these phases . with this purpose ,
we pay particular attention to the departure of correlation and response functions from stationarity and , we search for one of the hallmarks of a glassy phase , namely , violations of the fluctuation - dissipation theorem .
we also intend to discuss our results in view of previous and possible experimental tests .
the paper is organized as follows . in section [ ctp ]
we discuss the closed - time path - integral formalism of keldysh and schwinger and the coupling of the quantum system to a quantum heat bath .
we obtain a generating functional with an effective action that takes into account the effect of the coupling to the bath . in three appendices
we give details of the keldysh schwinger formalism ( appendix [ appkeldysh ] ) , we recall the quantum fdt ( appendix [ appfdt ] ) and we discuss the classical limit ( appendix [ classicallimit ] ) .
section [ meanfieldequations ] is devoted to the derivation of the dynamic equations through a saddle - point approximation of the generating functional .
the method is similar to the one used in the classical case when manipulating the martin - siggia - rose generating functional .
the general scenario for the real - time dynamics of this system is discussed in section [ scenario ] .
it is inspired by the behavior of the classical counterpart but some generalizations are necessary . in section [ equilibriumdynamics ] the equations and their solution in the paramagnetic phase
are studied with special emphasis set on the analysis of the out of phase susceptibility close to the transition line . in section [ nonequilibrium ]
the dynamic equations in the glassy phase are studied .
the non - equilibrium dynamics is shown to survive quantum fluctuations in a finite region of the phase diagram .
further details of the approach to solve the equations are given in appendix [ integrals ] .
the modification of fdt needed to construct an ansatz to solve these equations and its relation to effective temperatures are discussed in section [ generalfdt ] . finally , in section [ conclusions ] we present a summary of our results , our conclusions and some ideas as to how to continue investigating along the lines here proposed
. a short account of some of these results has appeared in ref .
the closed - time path - integral ( ctp ) formalism of keldysh and schwinger@xcite for treating the non - equilibrium dynamics of a quantum system has been reviewed several times in the literature .
@xcite in order to set the notation we briefly present it in section  [ generalframework ] .
we then consider the situation in which the system under study is coupled to a quantum thermal bath and derive the effective action following feynman and vernon @xcite ( section  [ couplingtoaheatbath ] ) . in section  [ ctpfordisorderedsystems ]
we finally address the case of a disordered system and we obtain the effective action from which the dynamic equations are derived in section  [ meanfieldequations ] .
the basic ingredient of this formulation is the `` in - in '' generating functional from which , after derivation with respect to the external sources @xmath23 , @xmath24 , the different green functions can be obtained .
let us consider a model with degrees of freedom described by a field @xmath25 in the heisenberg picture . denoting @xmath26 the density matrix at the initial time @xmath27 ( which from now on we set to zero ) , the generating functional
is defined as @xmath28   & = & \mbox{tr } \left[t^ * \exp\left(- \frac{i}{\hbar }   \int_{0}^\infty dt \ ,   { \bbox \xi}^-(t ) { \bbox \phi}(t)\right )   t \exp\left (     \frac{i}{\hbar } \int_{0}^\infty dt \ ,   { \bbox \xi}^+(t ) { \bbox \phi}(t ) \right )    \hat\rho(0 )   \right ] \
;   \label{generating}\end{aligned}\ ] ] where @xmath29 and @xmath30 are the @xmath15-vector external sources and a scalar product between @xmath15 vectors is assumed .
the symbols @xmath31 and @xmath32 are the time and anti - time ordering operators : @xmath33 defining ensemble averages in the usual way @xmath34 one obtains @xmath35 which is directly related to the feynman causal propagator . in the same way ,
@xmath36 { \ln } { \cal z } \right|_{{\bbox \xi}=0 } \ ; \end{aligned}\ ] ] where @xmath37 is the _ symmetrized correlation function _ ( @xmath38 ) .    the response function @xmath39 is defined as the variation of the averaged field @xmath40 with respect to the `` strength '' @xmath41 of a perturbation that modifies the potential energy according to @xmath42 : @xmath43 in linear - response theory , it can be expressed in terms of the averaged commutator : @xmath44   \rangle   \\ & = & \left .
\frac{\hbar}{i }    \left [ \frac{\delta^2}{\delta \xi_i^+(t ) \delta \xi_j^+(t ' ) } + \frac{\delta^2}{\delta \xi_j^+(t ) \delta \xi_i^-(t ' ) } \right ] { \ln } { \cal z } \right|_{{\bbox \xi}=0 } \ ; .\end{aligned}\ ] ]    one of the main advantages of this formalism is that the generating functional ( [ generating ] ) admits a path - integral representation which is a slight modification of the usual `` in - out '' feynman path integral .
this will allow us to make formal manipulations similar to the ones performed in the study of the dynamics of classical systems via the martin - siggia - rose@xcite ( msr ) formalism . in a short - hand notation , the generating functional reads @xmath28   & = & \int   { \cal d}{\bbox \phi}^+ { \cal d}{\bbox \phi}^- \exp\left [ \frac{i}{\hbar } \ ,
\left ( s[{\bbox \phi^+ } ] -s [ { \bbox \phi^- } ] + \int dt \ , { \bbox \xi}^+(t ) { \bbox \phi}^+(t ) - \int dt \ , { \bbox \xi}^-(t ) { \bbox \phi}^-(t ) \right ) \right ] \nonumber\\ & &   \;\;\;\;\;\;\;\;\;\ ;   \times   \langle { \bbox \phi}^+   | \hat
\rho(0 ) |   { \bbox \phi}^-   \rangle \label{aaab}\end{aligned}\ ] ] where @xmath45 $ ] is the action of the system and @xmath46 stands for the matrix element of the density matrix which has support only at @xmath47 ( see appendix [ appkeldysh ] for a proof of this equation ) .
hereafter , we omit the limits of the time - integrals that , unless otherwise stated , go from @xmath48 to @xmath49 .    in the following
it will be useful to write the correlation and response in terms of the fields @xmath50 , @xmath51 : @xmath52    notice that as a result of having two external sources ( one related to the evolution with time ordering and the other with anti - time ordering ) a doubling of degree of freedom is necessary when writing the path - integral . as we explicitly show in appendix [ classicallimit ] , the doubling of degrees of freedom characteristic of this description is intimately linked to the introduction of lagrange multipliers in the msr description of the classical model .
the time - integration in eq .
( [ aaab ] ) can be interpreted as being closed , going forward from @xmath53 to @xmath54 and then backwards from @xmath54 to @xmath53 .
this motivates the name of the method .      up to now
, the discussion has been completely general and it applies to systems with an arbitrary number of degrees of freedom .
we now consider the coupling of the system of interest to a thermal quantum bath , assumed to be in equilibrium
. a convenient way to deal with the interaction between the system and the bath in the path - integral formalism is due to feynman and vernon .
@xcite in this approach one starts by describing the system plus bath by @xmath55 with @xmath56 denoting the variables of the system and @xmath57 , @xmath58 , denoting the variables of the bath .
the total action is @xmath59= s_s[{\bbox \sigma } ] + s_b[{\bbox v}^a ] +   s_{sb}[{\bbox{\sigma , v}^a}]\ ] ] where @xmath60 is the action characterizing the system ( and eventually depending upon disorder ) , @xmath61 is the action for the bath , and @xmath62 contains the system - bath interaction terms . for definiteness
we assume that our system interacts linearly with the bath which we model as a set of independent harmonic oscillators , @xmath63 = s _ { s}[{\bbox \sigma } ] + \frac12   \sum_{a=1}^{n_b } \ ,
m_a \ ,   \left ( \dot { \bbox v}_a^2 - \omega_a^2 \ , { \bbox v}_a^2 \right )   - \sum_{a=1}^{n_b } c_a { \bbox v}_a { \bbox \sigma } \ ; .\ ] ] @xmath64 are their masses , @xmath65 their frequencies and @xmath66 the coupling constants between the oscillators and the system .
of course , more general couplings between bath and system are possible but this is the simplest choice that allows us to go a bit farther analytically . if we further assume that the initial density matrix @xmath67 is factorizable as @xmath68 where @xmath69 is the boltzman distribution for the bath variables at equilibrium at a temperature @xmath70 , one can integrate out the bath variables , obtaining in this way an effective thermal action @xmath71 $ ] , that enters the _ feynman - vernon influence functional _ : @xmath72 & = &   - \int dt \int dt ' \ ; \left ( { \bbox \sigma}^+(t ) -   { \bbox \sigma}^-(t ) \right ) \ , \eta(t - t ' ) \ ,    \left ( { \bbox \sigma}^+(t ' ) + { \bbox \sigma}^-(t ' ) \right )   \nonumber\\ & &   + i\int dt \int dt ' \ ; \left ( { \bbox \sigma}^+(t ) -   { \bbox \sigma}^-(t ) \right ) \ , \nu(t - t ' ) \ ,    \left ( { \bbox \sigma}^+(t ' ) - { \bbox \sigma}^-(t ' ) \right )   \ ; . \label{sthermal}\end{aligned}\ ] ] the noise and dissipative kernels @xmath73 and @xmath74 are given by @xcite @xmath75 in these equations , @xmath76 is the spectral density of the bath : @xmath77 the system - bath interaction manifests via the appearance of two quadratic and non - local terms ( [ sthermal ] ) in the effective action .
as it can be most easily seen in the classical limit presented in appendix [ classicallimit ] they can be related to the noise and dissipative terms of the associated langevin formulation of the classical dynamics .
in addition , it is easy to see that the kernels @xmath73 and @xmath74 are related to the correlation and response functions of the bath .
in fact @xmath78 is the response of the bath while @xmath79 is the correlation of the bath .
they are related in a way dictated by the fluctuation - dissipation theorem ( fdt ) : @xmath80 ( see appendix [ appfdt ] for the derivation and properties of the quantum fdt . )
this equation is the quantum counterpart of the fluctuation - dissipation relation of classical relaxational dynamics , that fixes the dissipative coefficient in terms of the strength of the noise correlation of the langevin equation .
it expresses the fact that the bath is assumed to be in equilibrium at the initial time @xmath48 and that it stays in equilibrium at all subsequent times @xmath81 .
it is important to notice that the assumption of thermal equilibrium of the bath _ does not imply _ that the system  that is in contact with it  is in equilibrium at the initial time nor that it will reach equilibrium at a given equilibration time @xmath82 .
it is known that `` mean - field '' classical disordered models with relaxational dynamics
the @xmath0 spin - glass model is an example  do not reach equilibrium at any time if the limit of @xmath15 going to infinity is taken at the outset of the calculation.@xcite it is one of the main aims of this paper to show that this result carries through to the quantum case . as a consequence of the out of equilibrium behavior , we shall find that the fdt _ does not hold _ between response and correlation of the system in a sub - critical region of the phase diagram .
the integrals in the @xmath74 and @xmath73 definitions are , in general , ill - defined ( divergent ) and might need to be regularized .
a way to regularize them@xcite is to observe that one expects @xmath76 to decrease to zero for large @xmath83 and introduce an explicit cut - off : @xmath84 by modifying @xmath76 in this way , we regularize the kernels @xmath74 and @xmath73 at the same time preserving the fdt relation for the bath variables ; in other words , we do not break the assumption of equilibrium of the bath .
different environments are characterized by different behaviors of @xmath76 such that for @xmath85 : @xmath86 a typical example is @xcite @xmath87 where @xmath88 is a constant that plays the rle of a friction coefficient .    throughout this paper
, we concentrate in the ohmic case .
the kernel @xmath74 is @xmath89 in the limit where the cut - off tends to infinity it becomes @xmath90 of particular importance is the zero temperature limit of the kernel @xmath73 : @xmath91 one sees that the quantum fluctuations yield a non - trivial kernel even in the absence of thermal fluctuations . in the zero temperature limit
the characteristic time goes to infinity and the long - time decay of @xmath73 becomes power - law : @xmath92 in the classical limit , for any temperature , @xmath73 becomes @xmath93 note that for a finite cut - off @xmath94 , the classical limit of the kernel @xmath73 also has a rather slow , power - law decay at large time differences .
if , next , the cut - off is sent to infinity one recovers a delta function @xmath95 characteristic of e.g. a white thermal noise .
the classical limit of these terms is further discussed in appendix [ classicallimit ] where we show that the dynamic equations of the corresponding classical model correspond to a friction coefficient @xmath96 .
the coupling to the quantum thermal bath introduces three time - scales into the problem .
first , the inverse cut - off @xmath97 characterizes the memory of the bath .
second , the coupling strength @xmath98 yields the inverse relaxational characteristic time .
third , @xmath99 accounts for the relative importance of quantum to thermal effects .      the ctp formalism results to be well suited for problems with disorder.@xcite let us suppose that the system under consideration is described by the action , @xmath100 = s_{\sc o}[\bbox{\sigma } ] -   \int dt \
, v[\bbox{\sigma } , j]\ ] ] where @xmath101 is the disordered independent part of the action and @xmath102 $ ] is the potential energy that depends upon random couplings collectively denoted by @xmath103 . as usual in disordered systems , we concentrate on quantities that are averaged over the disorder distribution @xmath104 $ ] .
we wish to consider the real - time dynamics of the system starting from a _ random initial condition _ at time @xmath48 when it is set in contact with the environment ( that is itself in equilibrium at a constant temperature @xmath31 ) .
given that the initial condition is chosen to be random , it is not correlated with the disorder .
therefore , the density operator @xmath67 does not depend upon disorder and the generating functional without sources @xmath105   = \mbox{tr } \ , [ \ , \hat\rho(0 ) \ , ] \ ] ] _ is also independent of disorder_.    this property is equivalent to the independence of the classical martin - siggia - rose ( msr ) generating functional @xcite without sources upon disorder.@xcite as in the classical case , it allows us to write dynamic equations for random initial conditions without having to compute the average over disorder of @xmath106 $ ] and hence without resorting to the use of replicas .
we are then interested in the averaged generating functional : @xmath107 } = \int d j p[j ] \ ;   { \cal z } [ \bbox{\xi}^+ , \bbox{\xi}^- , j ] \ ; , \ ] ] from which any averaged operator can be computed as @xmath108 } { \delta { \bbox   \xi}^+(t ) }   } \right| _ { { \bbox \xi } = 0 } = \left . \frac{1}{{\cal z}[0,0,j ] } \ ; \frac{\delta \overline{{\cal z}[{\bbox \xi}^+,{\bbox \xi}^-,j ] } } { \delta { \bbox   \xi}^+(t ) }   \right| _ { { \bbox \xi } = 0 } \ ; , \end{aligned}\ ] ] with all sources set to zero and @xmath109 .
here and in what follows , the overline represents an average over the disorder .    in many cases of interest
, the integration over the disorder can be performed explicitly , @xmath110 \ ; \exp \left[-\frac{i}{h }   \int dt   \left ( \ , v[{\bbox \sigma}^+ , j ] - v[\bbox{\sigma } ^- , j ] \ ,
\right ) \right]= \exp \left[-
\frac{i}{h } v_{\sc d } [ \bbox { \sigma}^+ , \bbox{\sigma}^- ] \right ] \;\ ] ] and it introduces a non - linear interaction term that is , in general , non - local both in time and in the `` space '' indices @xmath111 .
a rather general way of introducing randomness is to consider a model with a gaussian random potential energy term @xmath112 $ ] depending upon the variables @xmath56 and correlated as @xmath113 v[{\bbox \sigma } ' , j ] }   =   - \frac{n}{2 } { \cal v } \left ( \frac { ( { \bbox \sigma } - { \bbox \sigma}')^2}{n } \right )   \label{randompotential } \ ; .\ ] ] the variables @xmath19 may represent spins leading to a spin - glass model or a position in an @xmath15-dimensional space becoming the problem of the dynamics and/or diffusion of a particle in a random potential .
since we are interested in following the real time dynamics of such systems , @xmath19 will be a time - dependent @xmath15-vector .    for the sake of definiteness
, we shall solve a particular realization of the random potential that in the classical case defines the so - called @xmath0 spin - glass : @xcite @xmath114 =   - \sum_{i_1 < \dots
< i_p}^n   j_{i_i\dots i_p } \ ,
\sigma_{i_1 } \dots \sigma_{i_p }   \;.\ ] ] @xmath0 is a parameter , @xmath115 .
the interaction strengths are quenched independent random variables with a gaussian distribution @xmath116   =   \sqrt { \frac{n^{p-1}}{{\tilde j}^2\pi p ! } } \ ; \exp \left ( -   \frac{n^{p-1}}{{\tilde j}^2 p ! } \sum_{i_1
\neq \dots \neq
i_p }   ( j_{i_i\dots i_p})^2   \right )   \;\;\;\;\;\;\;\ ; \rightarrow \;\;\;\;\;\;\;\ ; \overline {   ( j_{i_i\dots i_p})^2 } = \frac{{\tilde j}^2 p!}{2 n^{p-1 } } \ ; .\ ] ] this corresponds to a random potential correlation @xmath117 and it leads to an effective , non - local in time , interaction @xmath118 & = &   \frac{{\tilde j}^2 ni}{4 \hbar } \int dt \int dt '   \left [ \ ; \left ( \frac1n { \bbox \sigma}^+(t ) { \bbox \sigma}^+(t ' ) \right)^p +   \left (   \frac1n   { \bbox \sigma}^-(t ) { \bbox \sigma}^-(t ' ) \right)^p \right . \nonumber\\   & & \left .   \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\ ; -   \left ( \frac1n { \bbox \sigma}^+(t ) { \bbox \sigma}^-(t ' ) \right)^p -   \left (   \frac1n   { \bbox \sigma}^-(t ) { \bbox \sigma}^+(t ' ) \right)^p \;\ ; \right ] \ ; . \label{actionv}\end{aligned}\ ] ]    in its spherical version,@xcite @xmath119 are continuous dynamic variables @xmath120 , @xmath121 , that satisfy the global spherical constraint @xmath122 at each instant .
we enforce this constraint by introducing a time - dependent lagrange multiplier @xmath123 .
the disordered - independent part of the action is @xmath124 =   \int dt \ ,   \left [ \frac{m}{2 } { \dot { \bbox \sigma}}^2 -   \frac{z}{2 } \left ( { \bbox \sigma}^2 -n \right )   \right ]   \label{actiono } \ ; .\ ] ]
as explained in the introduction , for simplicity we have chosen a kinetic term that dictates a microscopic dynamics with second derivatives in time .
other choices , for example leading to first time - derivatives , are also possible .
the derivation of the dynamic equations that we explain in section  [ meanfieldequations ] applies with only minor modifications .
we now have all the elements to derive the equations which describe the dynamic evolution of a disordered quantum system . according to the discussion in section [ ctp ] ,
the system is described by the generating functional @xmath125 } = \int   { \cal d}{\bbox \sigma}^+ { \cal d}{\bbox \sigma}^- \exp\left [ \frac{i}{\hbar } \ , \left ( s_{\sc eff}[{\bbox \sigma^+ } ,   { \bbox \sigma^- } ] + \int dt \ , { \bbox \xi}^+(t ) { \bbox \sigma^+}(t ) - \int dt \ , { \bbox \xi}^-(t ) { \bbox \sigma^-(t ) } \right ) \right]\ ] ] with @xmath126 = s_{\sc o}[{\bbox \sigma}^+ ] - s_{\sc o}[{\bbox \sigma}^- ] + s_{\sc t}[{\bbox \sigma}^+ , { \bbox \sigma}^- ]   - v_{\sc d}[{\bbox \sigma}^+ , { \bbox \sigma}^- ] \ ; . \label{action1}\ ] ]
we remind that @xmath101 is the disorder independent part of the action given by eq .
( [ actiono ] ) , @xmath127 accounts for the the system - bath interaction , eq .
( [ sthermal ] ) , and @xmath128 is the effective potential which arises as a result of the integration over the disorder and is given by eq .
( [ actionv ] ) .
the quantum dynamic equations of motion follow from similar steps to those usually used @xcite to obtain the classical equations of motion in the martin - siggia - rose formalism .
we introduce macroscopic order parameters and derive , through a saddle - point point approximation of the ks generating functional ( that becomes exact when @xmath129 ) , the dynamic equations of motion . since we are interested in considering the dynamics of the system for _ large but finite _ times with respect to @xmath15 , out of equilibrium effects are expected .
we shall therefore make no assumption about the time - dependence of the order - parameters .
the quadratic terms in the action can be condensed into one term by introducing the operator @xmath130 in such a way that @xmath126 =   -\frac12   \int dt   \int dt ' { \bbox \sigma}^\alpha(t ) \ , \mbox{op}^{\alpha\beta}(t , t ' ) { \bbox \sigma}^\beta(t ' )   - v_{\sc d}[{\bbox \sigma}^+ , { \bbox \sigma}^-]\ ] ] where greek indices label @xmath131 and the sum convention is assumed .    introducing the identity @xmath132 into the generating functional ,
the full action can be rewritten as @xmath133 & = &   - \frac{1}{2}\int dt   \int dt ' \ ;   { \bbox \sigma^\alpha}(t ) \ , \left ( \ , \mbox{op}^{\alpha\beta}(t , t ' )   +   \lambda^{\alpha\beta}(t , t ' )
\right ) { \bbox \sigma}^\beta(t ' )    \nonumber\\ & &   + \frac{n}{2 } \int dt \int dt ' \lambda^{\alpha\beta}(t , t ' ) \ ,   q^{\alpha\beta}(t , t ' ) + \frac{n}{2 }   \int dt \left ( z^+(t ) - z^-(t )
\right )   \nonumber\\ & &   + \frac{i { \tilde j}^2 n}{4 \hbar } \int dt \int dt '   \left [ \ , \left(q^{++}(t , t ' ) \right)^p +   \left ( q^{--}(t , t ' ) \right)^p \right .
\nonumber\\ & &   \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\ ; - \left .
\left ( q^{+-}(t , t ' ) \right)^p -   \left ( q^{-+}(t , t ' ) \right)^p \ , \right ] \ ; .\end{aligned}\ ] ] the stationary - point values of the order parameters @xmath134 with @xmath135 are related to the `` physical '' correlations and responses defined in eqs .
( [ defcorr ] ) and ( [ defresp ] ) as follows @xmath136 with @xmath137 it is easy to check that these functions satisfy the identity for all pairs of times @xmath138 . in the classical case
this is a requirement to be imposed upon the saddle - point two - point functions in order to preserve causality . in the quantum case
, this condition emerges as a property of the green functions . ]
@xmath139    the functional integration over @xmath140 and @xmath141 is now quadratic and can be performed .
this amounts to replace the quadratic term in @xmath142 by @xmath143      at this stage , all terms in the action depend upon the `` macroscopic '' quantities @xmath144 and @xmath145 and are proportional to @xmath15 .
since it is easier to write the equations in matrix notation , we encode @xmath146 and @xmath147 in two matrices @xmath148 and we define @xmath149(t , t ' )   & \equiv & \left ( \begin{array}{rl } ( q^{++}(t , t'))^{p-1 }   \ ; &   \ ; -(q^{+-}(t , t'))^{p-1 } \\
-(q^{-+}(t , t'))^{p-1 } \ ; &   \ ; ( q^{--}(t , t'))^{p-1 } \end{array } \right ) \label{saddlef } \ ; .\end{aligned}\ ] ] we denote with a cross the standard operational product @xmath150 where a sum over @xmath151 is assumed .
the saddle - point with respect to @xmath152 yields @xmath153 the matrix and time - operator inverse of @xmath154 is denoted @xmath155 .
the saddle - point equation with respect to @xmath134 yields @xmath156(t , t ' ) \ ; . \label{saddleq}\ ] ] equations  ( [ saddlelambda ] ) and ( [ saddleq ] ) imply , in a compact matrix and time - operator notation , @xmath157 \otimes { \cal q } \ ; , \label{saddleq}\ ] ] where @xmath158 is the identity : @xmath159 .
the saddle - point with respect to @xmath145 yields @xmath160 and these equations lead , as expected , to the spherical constraint .    in the limit @xmath129
one could also proceed as in refs .
[ ] and write the full action @xmath161 in terms of a _ single _ variable .
this is at the expense of modifying the thermal kernel and the interaction term in a self - consistent way , through the introduction of terms arising from the non - linear interactions ( the vertex @xmath162 and the self - energy @xmath163 , respectively ) .
this procedure is not of particular usefulness for the analysis of the model we treat in this paper since the single variable effective action is gaussian .
it does however become useful for dealing with quantum models whose single - spin effective action has higher order interaction terms .
an example is the quantum sherrington - kirkpatrick model .
one could then envisage to derive an effective quantum langevin equation for the single variable@xcite and study this equation with an adequate numerical algorithm as the one developed by eissfeller and opper.@xcite      the dynamic equations for the auto - correlation and response follow from the set of equations ( [ saddlef])-([saddleq ] ) and the definitions of the dynamic order parameters given in eqs .
( [ q++])-([q ] ) .
more precisely , the equation of motion for the response function follows from the subtraction of the @xmath164 and @xmath165 components of eq .
( [ saddleq ] ) : @xmath166 \ ,   r(t'',t ' ) \ ; , \label{eqr1}\end{aligned}\ ] ] and the equation of motion for the auto - correlation function follows from the addition of the @xmath165 and @xmath167 components of eq .
( [ saddleq ] ) : @xmath168 \label{eqc1 } \ ; .\end{aligned}\ ] ] written in this way , eq .
( [ eqc1 ] ) is complex .
its imaginary part yields @xmath169 moreover , since the response is causal , products of advanced @xmath170 and retarded @xmath171 responses vanish identically for all @xmath172 : @xmath173 and one can show that for any integer @xmath174 and any constants @xmath175 , @xmath176 @xmath177 ^k= \left [ c(t , t ' ) + c_1 r(t , t ' )   \right ] ^k + \left [ c(t , t ' ) + c_2 r(t',t ) \right ] ^k - c^k(t , t ' ) \ ; .\ ] ] using this property one has @xmath178^{p-1}\ ] ] and @xmath179 = \nonumber\\ & & 2 c(t'',t ' )   \mbox{im } \left [ c(t ,
t'')-\frac{i \hbar}{2 } r(t , t '' ) \right]^{p-1 } -\hbar r(t',t '' )   \mbox{re}\left [ c(t , t'')-\frac{i\hbar}{2 } ( r(t , t'')+ r(t'',t ) ) \right]^{p-1 } \ ; . \nonumber\\\end{aligned}\ ] ] we can identify the self - energy @xmath163 and the vertex @xmath162 as @xmath180^{p-1 } \ ; , \label{sigmatilde } \\
d(t , t')- 2 \hbar \nu(t - t ' ) & \equiv & \tilde{d}(t , t ' ) \equiv \frac{p \tilde{j}^2}{2 } \mbox{re}\left [ c(t , t')-\frac{i\hbar}{2 } ( r(t , t')+ r(t',t ) ) \right]^{p-1 } \ ; . \label{dtilde}\end{aligned}\ ] ] for @xmath181 , they can be encoded in a single complex equation : @xmath182 note that the total self - energy @xmath183 and vertex @xmath184 are real and have two contributions of different origin : one arises from the interaction of the system and the bath ( @xmath74 and @xmath73 ) and one is caused by the non - linearities stemming from the average over disorder ( that we called @xmath163 and @xmath162 in eqs .
( [ sigmatilde ] ) and ( [ dtilde ] ) ) . if @xmath13 , @xmath163 and @xmath162 are identical to the classical ones . instead , if @xmath185 , the non - linear terms acquire an explicit dependence upon @xmath186 .
the dynamic equations can then be written in a compact form @xmath187 that we call later the @xmath188-eq . and @xmath189-eq
. , respectively .
it is important to realize that the self - energy @xmath190 is proportional to the response function @xmath191 , which in turns implies @xmath192 this means that the upper limit of integration in eqs .
( [ schwingerr ] ) and ( [ schwingerc ] ) is @xmath193 , which renders the equations explicitly causal .
there are no more independent equations for @xmath188 and @xmath189 .
the other two equations that can be obtained from eq .
( [ saddleq ] ) are the equation for @xmath194 , that is equivalent to the @xmath188-eq above , and one equation that identically cancels by virtue of the identity ( [ identityqo ] ) .    in their integrated form as schwinger - dyson equations the dynamic equations read : @xmath195 with the propagator given by @xmath196    real and imaginary parts of eqs .
( [ z+ ] ) and ( [ z- ] ) combined with eq .
( [ saddlelambda ] ) imply the equal - times conditions @xmath197 in addition , from eq .
( [ eqr1 ] ) one obtains that the first derivative of the response function is discontinuous at equal times : @xmath198 while from eq .
( [ eqc1 ] ) one obtains that the correlation is continuous : @xmath199    the equation for @xmath200 can be derived from the schwinger - dyson equation ( [ intschwingerc ] ) , by imposing the spherical constraint through the evaluation at @xmath201 .
multiplying operationally by @xmath202 one obtains @xmath203   \nonumber\\ & &   + m \int_0^t dt '' \int_0^t   dt ' '' \ ,     ( \partial_t r(t , t '' ) ) \,d(t'',t''')\ , ( \partial_t r(t , t ' '' ) ) \ ; .\end{aligned}\ ] ] the last term is a consequence of having a kinetic term with second derivatives . had we chosen a first derivative term in the kinetic energy the last term would have not appeared .
it can be easily identified with minus the second - derivative of the correlation at equal times by taking the limit @xmath204 in eq .
( [ schwingerc ] ) .
thus @xmath205   - m \left .
\frac{\partial^2 } { \partial t^2 } c(t , t ' )   \right|_{t'\to t^- } \ ; .
\label{zeq}\ ] ]    in conclusion , eqs .
( [ schwingerr ] ) , ( [ schwingerc ] ) and ( [ zeq ] ) are the complete set of equations that determines the dynamics of the system . in the following
we make a time reparametrization @xmath206 that implies @xmath207 and @xmath208 the system of units we use is @xmath209 & = & [ \hat r ] = [ \hat \hbar ] = 1 \ ; , \nonumber\\ \left[\hat m \right ]   & = &   [ \hat \beta ] = [ 1/\hat{\tilde j } ]   = [ 1/\hat{\lambda } ] = [ \hat{t } ]   \ ; .\end{aligned}\ ] ] hereafter we write the equation in reparametrized variables and drop all hats .
as already mentioned , the set of integro - differential equations ( [ schwingerr ] ) , ( [ schwingerc ] ) and ( [ zeq ] ) is causal  as in the classical case  and one can attempt a numerical solution by a simple algorithm .
the correlation @xmath210 and response @xmath191 are _ two time quantities _ , that is , they depend on @xmath193 ( which physically corresponds to the time of observation ) and @xmath211 ( which corresponds to the age of the system ) . in a normal situation
, one expects that after an equilibration time @xmath82 ( in general model dependent ) all two - time functions will depend only on the time difference @xmath212 .
in other words , all two - time functions become time - translation invariant ( tti ) . under these circumstances ,
to which we refer as equilibrium evolution , response and correlation are not independent quantities but are linked via the fluctuation - dissipation theorem ( fdt ) .
nevertheless , it is also known that certain systems never reach ( at least for experimentally accessible time - scales ) an equilibrium dynamics . for arbitrary large @xmath211 and in a certain region of the parameter space , although one - time quantities , such as the energy - density , tend to an asymptotic limit , two - time quantities , such as the correlation and response , depend both on @xmath193 and @xmath211 in a non - trivial way .
this is indeed the typical situation for a glassy system.@xcite    the question we would like to explore is whether eqs .
( [ schwingerr ] ) , ( [ schwingerc ] ) and ( [ zeq ] ) encode a non - equilibrium evolution as described above .
we know that indeed this is the case for classical disordered systems . as shown in ref .
[ ] , this kind of equations describe , at least partially , the non equilibrium dynamics of glassy systems , notably _ aging _ effects of spin - glasses
. the salient feature of these solutions is that , below a critical temperature @xmath213 , and for large @xmath211 , there exist at least two time scales : for @xmath212 small with respect to a characteristic time @xmath214 , that depends upon the age of the system , the dynamics is similar to an equilibrium evolution ( stationary regime ) but for @xmath212 large with respect to @xmath214 the dynamics becomes extremely slow with aging effects . in this regime ,
correlation and response are related in a manner that violates fdt .    in this paper
we show that this also happens in the quantum problem in a region of parameter space . in what follows we concentrate on the case @xmath17 . as in the classical case
, the @xmath13 model is connected to @xmath14 models and has a different behavior . in fig .
[ phasediagram ] we present a schematic phase diagram @xmath215 for @xmath216 .
as we show below , the dynamics in the disordered (  paramagnetic ) and ordered (  spin - glass ) phases are characterized by
equilibrium and `` out - of - equilibrium '' effects , respectively .
the dynamics for the quantum system in the glassy phase has a similar pattern to the classical one , in the sense that there are at least two time - sectors in which the evolution is of a different nature .
inspired by the classical problem,@xcite we propose that also in the quantum case the weak - ergodicity breaking @xcite and the weak long - term memory @xcite properties hold .
these two assumptions are based in part on the insight coming from the numerical solution of the full equations .
the _ weak - ergodicity breaking scenario _ states that , for @xmath181 , the correlation function decays in such a way that @xmath217 in these equations @xmath218 is the edwards - anderson order parameter that characterizes the spin - glass phase .
this means that for large @xmath181 such that @xmath212 is small with respect to @xmath214 , the correlation function first decays from 1 to @xmath218 in a time - translational invariant manner .
@xmath214 is a growing function of @xmath211 whose precise form depends on the model .
the correlation goes further below @xmath218 to eventually reach @xmath219 in a manner that depends both upon @xmath193 and @xmath211 ( the aging effect ) .
this behavior suggests the presence of at least two time - sectors in which the dynamics is stationary and non - stationary , respectively .
we then write the correlation as the sum of a stationary and an aging contribution : @xmath220 the matching conditions at equal times between @xmath221 and @xmath222 are @xmath223 with @xmath224 together with eq .
( [ limitqq ] ) they ensure that in the two - time sector in which @xmath221 decays from @xmath225 to @xmath219 , @xmath222 is just a constant @xmath218 .
instead , in the time - sector in which @xmath222 decays from @xmath218 to @xmath219 , @xmath221 vanishes identically .
the name weak - ergodicity breaking @xcite reflects the fact that for short time - differences the system behaves as if it would be trapped in some region of phase space of `` size '' @xmath218  suggesting ergodicity breaking .
however , it is always able to escape this region in a time - scale @xmath214 that depends upon its age @xmath211 .
hence , trapping is gradual and ergodicity breaking is _
weak_.    we recall that in the classical and purely relaxational case ( langevin dynamics ) the correlation functions are monotonic with respect to both times @xmath193 and @xmath211 . however , in the presence of inertial terms there can be oscillations and the decay can be non - monotonic .
this will depend upon the relative value of the mass @xmath226 with respect to the other parameters in the problem . for the parameters we study in this paper , the oscillations appear only in the stationary regime , the aging dynamics having a monotonic decay towards zero .
this is relevant since it will allow us to use the general properties of monotonic correlation functions proven in ref .
[ ] to find the two - time scaling of @xmath227 .
regarding the response function , we propose @xmath228 with @xmath229 when a system is in equilibrium , the response is simply related to the correlation via fdt . here
we need to extend this relation to a non - equilibrium quantum situation . as in the classical case
we assume that the dynamics in the stationary regime is like an equilibrium dynamics and satisfies fdt .
we then assume that @xmath230 and @xmath231 are related by @xmath232 ( see appendix [ appfdt ] ) .
in contrast , in the time - sector in which the dynamics is non - stationary and manifestly non - equilibrium , fdt will be modified . in section [ generalfdt ]
we suggest an extension of fdt to quantum non - equilibrium systems with slow dynamics , that we later check in the particular model under study .    the _ weak long - term memory scenario _ states@xcite that the system keeps a weak memory of what happened in the past .
more precisely , the response function tends to zero when times get far apart , @xmath233 and its integral over a _
finite _ time - interval also vanishes @xmath234 but , its integral over an interval that grows with time gives a finite contribution : @xmath235 note that the property ( [ shortr ] ) implies @xmath236    an important quantity that reflects the decay of the response function is the two - time susceptibility defined as the integral of the response over the time - interval @xmath237 $ ] : @xmath238 ( @xmath239 is the strength of the applied perturbation , see eq .
( [ respdef ] ) ) . in the classical case
, a useful description of the departure from fdt is obtained by plotting @xmath240 vs @xmath241 for several @xmath242 using @xmath243 as a parameter . since the classical fdt implies @xmath244 , violations of the theorem manifest as departures from a straight line plot . in the quantum case ,
the relation between @xmath188 and @xmath189 is more involved and _ a priori _ one does not expect these plots to be useful . as we show in section [ generalfdt ] , it turns out that even for a quantum system these parametric plots are relevant .
in the paramagnetic phase ( see fig .  [ phasediagram ] ) one expects that , after a short non - equilibrium transient , a time - translational invariant ( tti ) solution establishes . we show that a tti ansatz implies that fdt holds between @xmath188 and c and between @xmath163 and @xmath162 . by solving numerically the full set of equations ( [ schwingerr ] ) , ( [ schwingerc ] ) and ( [ zeq ] )
we show the correctness of this ansatz and later discuss some properties of the solutions .    after the short initial transient , one - time quantities such as @xmath123
should reach a limit @xmath245 .
we evaluate it as @xmath246 - m \left .
\partial^2_{t } c(t - t ' )   \right|_{t'\to t^{-}\to\infty } \ ; .\ ] ] the fourier - transform of eqs .
( [ schwingerr ] ) and ( [ schwingerc ] ) yields _ the quantum mode - coupling equations _ : @xmath247 thus , @xmath248    if we now assume that @xmath163 and @xmath162 are related by fdt , it is easy to check that eqs .
( [ 111 ] ) , ( [ 222 ] ) and ( [ imromega ] ) imply that @xmath188 and @xmath189 are also related by fdt . the proof is completed by showing that fdt between @xmath188 and @xmath189 implies fdt between @xmath163 and @xmath162 .
this is easier to do in the time - domain , by verifying that eq .
( [ fdttemporal ] ) holds between @xmath163 and @xmath162 if it does between @xmath188 and @xmath189 .
a solution of this type exists only in a restricted region of the phase diagram . in order to demonstrate the existence of a dynamic transition
, we solve the _ full _ dynamic equations numerically and find that at a critical line , both the tti and fdt assumptions break down .
for the purpose of solving eqs .
( [ schwingerr ] ) , ( [ schwingerc ] ) and ( [ zeq ] ) numerically , we proceed as in the purely relaxational classical case.@xcite notice though that the presence of non local kernels @xmath74 and @xmath73 , that resemble the delta - function and appear convoluted with the correlation and response , renders the numerical solution harder .
the larger the cut - off @xmath94 , the smaller the iteration step @xmath249 we need to compute these integrals with a good precision .
we found an acceptable data collapse for @xmath250 if @xmath251 .
note that the values of @xmath249 typically used in the purely relaxational case are at least one order of magnitude larger .
this is the reason why the time - window we explore here is rather narrow though it suffices to show the trend of the solution .    even after rescaling time and all parameters in the problem , there are still many free parameters left : @xmath252
we here concentrate on the dependence upon only a small fraction of this set .
we fix hereafter @xmath216 , @xmath253 , @xmath254 and @xmath251 and study the behavior upon the remaining parameters @xmath255 .
we discuss the data for the correlation , response and integrated response for three points in the paramagnetic phase :    * @xmath256 .
this corresponds to the classical problem with colored noise and inertia .
the classical critical temperature for the purely relaxational case is @xmath257 . in fig .
[ fight2 ] we show the correlation @xmath258 vs. the time - difference @xmath259 for @xmath260 . in fig .
[ fight3 ] we plot the response @xmath261 vs @xmath259 for the same waiting times . in both cases
the decay is very fast , typical of upper critical dynamics , and tti quickly sets in ( cfr . the curves for @xmath262 and @xmath263 that are indistinguishable in both plots ) .
in the inset , we plot the integrated response @xmath264 vs. @xmath258 for @xmath265 using @xmath259 as a parameter .
the @xmath266 vs. @xmath189 plot very soon becomes a straight line of slope @xmath267 showing that the classical fdt is satisfied .
* @xmath268 .
this corresponds to a strong quantum regime without thermal noise .
the data for correlations and response in figs .
[ fighhbar3 - 1 ] and [ fighhbar3 ] are qualitatively similar to the ones for the classical problem at high temperatures in the sense that an equilibrium regime is quickly attained .
the decay is fast and tti holds . in the inset of fig .
[ fighhbar3 ] we display the @xmath266 vs. @xmath189 plot . in this strong quantum case at @xmath7 , the curve severely deviates from the classical expectation of a straight line of infinite slope .
the quantum fluctuation - dissipation theorem is verified in fig .
[ rfoul5hbar6t0 ] by comparing the response function obtained numerically with the convolution of the kernel @xmath269 with the numerical correlation function , see eq .
( [ fdtfourier ] ) .
these two curves should be equal in equilibrium , as imposed by fdt , and we here check that this is indeed the case in the paramagnetic phase . *
this is the quantum problem with both thermal and quantum fluctuations in the paramagnetic phase . in fig .
[ mixedhb1t2 ] we check fdt for these parameters and observe , in the @xmath266 vs. @xmath189 curve displayed in the inset , that the fdt becomes the classical .
we continue the study of the paramagnetic phase by analyzing the variation with @xmath31 and @xmath186 , close to the transition line , of the out of phase susceptibility @xmath271 this is the quantity that is most commonly measured experimentally . for instance , in the context of the study of quantum phase transitions , the dynamics of the randomly diluted , dipolar coupled , ising magnet liho@xmath2y@xmath3f@xmath4 was studied by measuring,@xcite among other quantities , @xmath272 .
we follow two routes to the transition :    * @xmath273 , @xmath274 .
one approaches the classical model by decreasing the strength of quantum fluctuation at fixed temperature @xmath275 . in fig .
[ chit13hbar ] we plot @xmath272 vs. @xmath83 . in the curves for small @xmath186
we observe a local minimum that is associated to the appearance of a plateau in the correlation function , at the value @xmath218 , when approaching @xmath18 .
this is typical of `` discontinuous transitions '' as the one in the classical @xmath0-spin model@xcite ( and in super - cooled - liquids@xcite ) .
we here note that this feature survives the presence of small quantum fluctuations and progressively disappears when @xmath186 increases . for @xmath276 ,
the minimum has been substituted by a single peak .
* @xmath277 , @xmath278 .
one approaches the critical line @xmath279 from above by modifying the strength of the quantum fluctuations while keeping the temperature fixed @xmath280 .
figure [ chit01hbar ] displays @xmath272 vs. @xmath83 for this series of @xmath186 .
( we can not get closer to the transition line because the dynamics gets so slow that the stationary limit can not be reached . )
the curves have a distinct peak that moves towards decreasing frequencies for decreasing @xmath186 .
the spectral width of the bell - shaped curves also increases when approaching the critical line .
the low frequency part of the curves is importantly suppressed with increasing @xmath186 .
all these features show that quantum fluctuations radically affect the relaxation . +
contrary to fig .
[ chit13hbar ] , there is no local minimum in any of the curves .
this is in accordance with a decreasing and eventually vanishing plateau value @xmath218 when the parameters approach the quantum critical point .
we shall further discuss the behavior of the model close to the transition line in section [ criticalline ] .
these results demonstrate the existence of a cross - over in the disordered phase , separating a region where the scaling laws are controlled by the quantum critical point @xmath281 from another region with different scaling laws controlled by the classical critical point @xmath282 .    before concluding our analysis of the paramagnetic phase ,
let us compare our results with the situation encountered for the system@xcite liho@xmath2y@xmath3f@xmath4 , which is believed to be an experimental realization of the transverse ising model .
naturally , we do not expect to find the same behavior in both models as we know that in the classical limit they are different:@xcite the classical transition being first order for the @xmath0-spin model and second order for the ising case .
our analysis indicates that the dynamic transition of the @xmath0-spin model remains first order along the critical line but becomes second order at the quantum critical point .
this is the mirror situation to what has been observed in the dipolar glass . in ref .
[ ] , it has been suggested that the quantum spin - glass transition may be first order .
figure [ chit13hbar ] displays the effect of quantum fluctuations close to the classical critical point and looks rather similar to fig .  2 in ref .  [ ] which instead shows the dynamics close to the quantum critical point .
figure [ chit01hbar ] looks similar to fig .  4 in ref .
[ ] , where the dynamics close to the classical critical point is presented .
this comparison indicates that this model belongs to a different class from the system liho@xmath2y@xmath3f@xmath4 and probably from the transverse field ising model.@xcite
in this section we focus on the non - equilibrium dynamics observed in the glassy phase . in order to justify the scenario discussed in section [ scenario ] , we first present the results from the numerical solution of the full equations .
next , we analyze these equations in the stationary and aging regimes and obtain a self - consistent relation that determines the edwards - anderson parameter .      the numerical solution shows a very different behavior below the critical line @xmath283 . in order to illustrate the slow dynamics , the breakdown of tti and the weak - memory of the system we choose the following parameters :    * @xmath284
this is the classical , glassy model studied in ref .
[ ] but with the addition of inertia .
the effect of inertia amounts to the presence of oscillations around @xmath218 but it does not change the qualitative behavior of the model .
for example , as we shown analytically below , the value of @xmath218 is not modified by inertia if @xmath226 is small .
* @xmath285 , @xmath286 .
this corresponds to pure quantum fluctuations .
figure  [ ct0hb1 ] shows the correlation @xmath258 vs. @xmath259 for the waiting times @xmath287 .
the curves have a stationary and an aging regime separated by the edwards - anderson parameter @xmath288 , see eq .
( [ limitqq ] ) .
the correlations decay to zero at far apart times , as in eq .
( [ limit0 ] ) .
this demonstrates the weak ergodicity breaking scenario .
+ in fig .
[ rt0hb1 ] we plot the response @xmath261 vs. @xmath259 for the same waiting times .
one observes a small time - difference regime , @xmath289 , where tti rapidly establishes and a large time - difference regime , @xmath290 , where the response shows a decaying queue with a @xmath242 dependence .
the integration of this queue over a time - interval scaling with @xmath242 will tell us about the violations of fdt in the quantum case ( see fig .
[ ml5hb01t0 ] below ) .
it also shows that the system has a weak long - term memory as described in section  [ scenario ] . *
@xmath291 , @xmath286 .
the qualitative behavior is similar to that obtained for @xmath7 and @xmath292 apart from the fact that thermal fluctuations slightly decrease the value of the edwards - anderson parameter and accelerate significantly the decay in the aging regime , as demonstrated in fig .
[ ct05hb1 ] .
the lagrange multiplier @xmath123 is a one - time quantity , i.e. it depends only upon the total time @xmath193 . in this formalism
we assume that all one - time quantities reach a well - defined limit asymptotically . in appendix [ integrals ]
we describe how we obtain the following equation for the asymptotic limit of @xmath123 : @xmath294 - m \left .
\partial^2_\tau c_{\sc st}(\tau)\right|_{\tau\to 0 } \ ; .\end{aligned}\ ] ] let us discuss each term in this expression .
we called @xmath295 the aging contribution : @xmath296 \ ; . \end{aligned}\ ] ] the kernels @xmath73 and @xmath74 , that are related to the bath , do not contribute to @xmath295 .
this is so because we have assumed that for @xmath297 the kernels have already decayed to zero in such a way that they do not contribute to these integrals .
more precisely , we are neglecting terms of the form @xmath298 where @xmath299 is either @xmath73 or @xmath74 and @xmath300 is either @xmath222 or @xmath301 .
the second and third terms come from the constant ( non - zero ) limit of the first decay of the correlation @xmath302 and equivalently of the vertex @xmath303 in the model under study @xmath304 if we use @xmath305 , a property of the weak long - term memory scenario .      if @xmath306 are such that @xmath307 , the discussion in section [ scenario ] implies @xmath308 and @xmath309 .
the schwinger - dyson equation for @xmath188 in this time sector reads @xmath310 and it keeps the same form as in the high - temperature regime , apart from the fact that the constant @xmath245 has contributions from the aging regime .
the schwinger - dyson equation for @xmath189 reads @xmath311 \ ; . \label{cst}\end{aligned}\ ] ] one can now fourier - transform both equations @xmath312 the formal solution to the equation for @xmath221 is @xmath313 the first term on the right - hand - side has an imaginary and a real part . the imaginary part vanishes identically since , due to fdt , both @xmath314 and @xmath315 are proportional to @xmath316 which is zero at @xmath317 .
concerning the real part of this first term , as we have assumed that @xmath231 goes to zero for @xmath318 , we need to impose the self - consistent condition @xmath319 this is the condition that fixes the edwards - anderson parameter
. we shall find it again in the next section as the matching condition between the stationary and aging regimes .
the final equation for @xmath320 is @xmath321 one can check that these calculations are consistent with the results from @xmath245 . actually , the integrals in @xmath21-eq . involving the stationary parts
can be evaluated with the help of the equations for @xmath322 and @xmath221 , eqs .
( [ rstat ] ) and ( [ cstat ] ) , and yield once again eq .
( [ cond2 ] ) .    similarly to the high - temperature case
one can now show that fdt for @xmath323 and @xmath324 implies fdt for @xmath322 and @xmath221 .
the remainder of the proof , i.e. to show that fdt between @xmath322 and @xmath221 implies fdt between @xmath323 and @xmath324 depends only upon the from of @xmath323 and @xmath324 as functions of @xmath322 and @xmath221 and is not modified from the one discussed in section  [ equilibriumdynamics ] .      if we now choose the times @xmath138 to be well - separated so as to have @xmath325 and @xmath326 , the weak - ergodicity breaking and weak long - term memory hypotheses allow us to throw the second time derivatives on the left - hand - side .
we assume that their contribution is much weaker than the one of each of the integral terms on the right - hand - side .
this is an assumption that we have to verify at the end of the calculation , once the solution for @xmath222 and @xmath301 is known .
it corresponds to the over - damped limit : @xmath327    using the approximation described in appendix [ integrals ] , the @xmath188-eq in the aging regime becomes @xmath328 and we call it the @xmath301-eq .
similarly , the @xmath189-eq becomes @xmath329 and we call it the @xmath222-eq . in
all integrals we approximated @xmath330 and @xmath331 since @xmath332 and @xmath333 are very rapidly decreasing functions .
the edwards - anderson parameter @xmath218 is determined self - consistently from the matching of @xmath334 .
taking the limit @xmath204 in the @xmath301-eq and @xmath222-eq one obtains @xmath335    the first equation admits the solution @xmath336 since @xmath337 is proportional to @xmath338  see eq .
( [ sigmatilde ] ) .
this corresponds to the high - temperature solution where there is no aging regime .
here we concentrate on the other possibility , that is to say when @xmath339 in actual fact , we can further approximate this equation by using one of the assumptions already used all over this section : that the response becomes smaller and smaller as time passes  though its integral over an infinite interval gives a finite contribution .
if we neglect all terms that are proportional to @xmath338 with respect to terms that are proportional to @xmath218 , only the first term in the power expansions of @xmath163 and @xmath162 survive : @xmath340 that in our model become @xmath341 in accord with the large @xmath259 limit of the stationary values ( [ dqstat ] ) .
equations ( [ hola ] ) and ( [ cq ] ) become @xmath342 the second equation is the same as the one arising from the end of the stationary regime , eq .
( [ cond2 ] ) .
we now study these equations with the help of eqs .
( [ rstat ] ) and ( [ cstat ] ) for the stationary functions . using @xmath343 one
has @xmath344 we remind that the factor @xmath345 can be written in terms of the stationary correlation function using fdt ; therefore this is a closed equation for the correlation that gives the edwards - anderson parameter @xmath218 . in the case of the @xmath0-spin model it reads @xmath346 in the classical case , the integral can be readily computed and the final equation for @xmath218 is @xmath347 that coincides with the result for the purely relaxational dynamics.@xcite for @xmath17 fixed , @xmath218 is a function of temperature ; it equals one at @xmath7 and tends to @xmath348 at the dynamic critical temperature @xmath213 .
if @xmath349 , the edwards - anderson parameter @xmath218 depends upon @xmath186 and @xmath31 , @xmath350 .
it decreases when thermal and/or quantum fluctuations increase .
one observes @xmath351 . at the critical line @xmath352 , @xmath353
as long as @xmath354 . at the _ quantum critical point _ the nature of the transition changes and @xmath355 .
we shall explain this result after introducing the modification of the quantum fdt needed to complete the solution .
figure  [ plotq ] shows the tendency of @xmath356 for @xmath357 .
one of the important characteristics of the dynamic evolution in equilibrium is the fluctuation - dissipation theorem theorem ( fdt ) that establishes a model independent relationship ( in the sense that only involves @xmath31 and @xmath186 but no other parameter in the model ) between correlation and response functions . as recalled in appendix [ appfdt ] , this result can be easily obtained , for the quantum and classical cases , under the hypothesis of equilibrium . in general , such a relation is not expected to hold in a non equilibrium situation and we can then say that fdt is  violated " ( of course what is violated is not the theorem but the hypothesis under which the result is derived ) . it is nevertheless surprising that in the glassy phase of certain classical models , the violation of fdt @xcite takes such a simple form that it is possible to formulate a generalized version of this relation through the introduction of an effective temperature @xmath358 .
@xcite in this way the fdt takes the standard equilibrium form but with a model dependent effective temperature which is different from the bath temperature @xmath31 ( in general @xmath358 is higher or equal than @xmath31 and it can take different values in different time - sectors ) .
we would like to explore in this section how these ideas carry over to the quantum case .
the standard quantum fluctuation - dissipation theorem states @xmath359 with @xmath360 guided by the classical case , we introduce the effective temperature@xcite @xmath361 and we propose , as a generalization of fdt to the out of equilibrium case , the following relation @xmath362 with @xmath363 and @xmath364 a function of both times @xmath193 and @xmath211 . evidently , eqs .
( [ xfdt1 ] ) and ( [ xfdt2 ] ) reduce to eqs .
( [ fdt1 ] ) and ( [ fdt2 ] ) when @xmath365 if the evolution is tti and the factor @xmath366 is set to one .
let us assume , for reasons that are intimately related to the time - separation discussed in section [ scenario ] , that @xmath364 is characterized by @xmath367 for times in the stationary regime @xmath368 this implies : @xmath369 with @xmath370 if we suppose that ( even at @xmath7 ) @xmath371 is finite in the glassy phase , it is then clear that in the aging sector the fourier integral is dominated by @xmath372 . after using @xmath373
we obtain @xmath374    for a standard quantum system in equilibrium at finite temperature , there is a characteristic time , that depends upon temperature and the strength of quantum fluctuations , beyond which the evolution becomes classical .
concerning fdt , this means that the classical relation establishes with @xmath375 .
this is clear , for instance , for the model under consideration at finite temperature in the paramagnetic phase ( see the inset in fig .
[ fight3 ] ) .
if our assumptions are correct , the situation in the glassy phase is different . for small times
compared with @xmath376 , now a _
waiting - time dependent _ characteristic time , correlation and response are related by the standard ( quantum ) fdt . instead , for time differences larger than @xmath376 ,
they are related in the same way as in the classical case _ but _ with a @xmath31 and @xmath186 dependent effective temperature that is different from the one of the environment .    in general , the proposal given in eq .
( [ fdt6 ] ) will allow us to reduce the @xmath222-eq .
and @xmath301-eq equations to single equation for @xmath222 .
it is simple to check that , at least for the model under consideration , this equation has the same structure as the one for the classical problem .
all dependence on @xmath186 enters through @xmath218 , @xmath377 and @xmath322 .
hence , once the explicit form of @xmath377 is known , the aging equations are solved by the same ansatz as in the classical case : @xmath378 note that the general arguments introduced in ref .
[ ] to characterize the behavior of _ monotonic _ correlation functions apply to the correlations in the aging regime and justify the solution ( [ agingc])-([agingr ] ) . even if the inertia might make the correlation oscillate , it does only in the stationary regime since it oscillates around @xmath218 , and it becomes monotonic in the aging    in eq .
( [ xtwoscales ] ) , we assumed that @xmath364 can take only two values in view of the model we study .
however , in the classical case , a more general situation arises in models with multiple time - scales , for which the factor measuring the violation of fdt is a non - constant function of the correlation in the aging regime.@xcite we expect the same structure to carry through to quantum problems whose classical counterparts are of the multi - scale type .
we check , both analytically and numerically , the viability of our assumptions for the @xmath0-spin model . the factor @xmath377 measuring the violation of fdt is determined by eqs .
( [ firstq ] ) and ( [ secondqq ] ) @xmath380 using eq .
( [ fdt6 ] ) and the equivalent relation between @xmath381 and @xmath382 , we obtain @xmath383 and @xmath384 in the classical limit @xmath385 and the result in ref .
[ ] is recovered .
note that both in the classical and quantum case , @xmath386 if @xmath13 . since the case @xmath13 is formally connected to ferromagnetic domain growth ( in the mean - field approximation )
there is still no - memory in the quantum domain growth .    the classical critical point @xmath387 is characterized by@xcite @xmath388 following the critical line towards the quantum critical point @xmath389 , one obtains @xmath390 and a finite zero - frequency stationary susceptibility @xmath391 . at the quantum critical point
, @xmath218 vanishes and the zero - frequency stationary susceptibility diverges : @xmath392 ( it is interesting to notice that the case @xmath13 is clearly different . ) on the other hand , if approaching the quantum critical point from below along the axis @xmath393 @xmath394 then @xmath395    [ ml5hb01t0 ]    we next check numerically the relation between response and correlation in the glassy phase , @xmath396 corresponding to pure quantum fluctuations .
first , we use fourier analysis to check the quantum fdt for short time differences . in fig .
[ fdthbarfourier ] we compare two different ways of obtaining @xmath397 for total time @xmath193 fixed and equal to @xmath398 .
the thin curve is the result from the direct numerical solution of the dynamic equations and shows a fast evolution for short time differences and a tail that gives rise to the weak memory of the system .
the bold curve is obtained from eqs .
( [ fdt3 ] ) and ( [ fdt4 ] ) using @xmath399 with @xmath400 obtained numerically .
we have estimated the edwards - anderson parameter as @xmath401 .
one observes that for small time - differences the two curves coincide showing that fdt holds in this time - sector .
second , we check the generalized fdt in the aging regime . for this purpose
, it is convenient to use the @xmath266 vs. @xmath189 plots defined in section [ scenario ] . in fig .
[ chic1 ] we observe , as predicted by eqs .
( [ fdt3 ] ) , ( [ fdt4 ] ) and ( [ fdt6 ] ) , two behaviors according to the relative value of @xmath189 and @xmath218 .
the portion of the graph that is significant for the aging regime corresponds to @xmath402 . as claimed in eq .
( [ fdt6 ] ) the @xmath266 vs @xmath189 curve approaches a time - independent asymptotic limit that is a straight line of slope @xmath403 . in this case , @xmath404 and this shows that even a @xmath7 , @xmath405 .
the effective temperature depends on @xmath31 and @xmath186 in a non trivial way . in figure [ chic2 ]
we show the @xmath266 vs. @xmath189 plot for @xmath186 = 1 and @xmath7 . in this case @xmath406 .    in fig .
[ teffmixed ] we work at fixed @xmath407 and study the variation of @xmath358 with t. from the @xmath266 vs @xmath189 plots we determine how @xmath358 varies with @xmath31 by measuring the slopes of the different curves and by using eq .
( [ xfinal ] ) .
the six curves above the dotted straight line correspond to temperatures below @xmath1 . from top to bottom @xmath408 .
one notices that @xmath409 in all these curves .
the curve below the dotted straight line corresponds to @xmath410 .
it has a slope @xmath411 ( indicated by the dots ) and @xmath412 . in the inset
we plot the temperature dependence of the effective temperature .
@xmath358 is linear in @xmath31 when the bath temperature is above @xmath1 but severely deviates from the linear behavior for temperatures below the critical line .
two curious features are that , in the glassy phase , @xmath358 increases when either @xmath31 or @xmath186 decrease .
the surprising increase of @xmath358 as @xmath31 decreases is also obtained in the classical limit.@xcite
in this paper we presented a formalism that allows us to analytically study the behavior of quantum glassy models within the typical situation encountered in experiments , i.e. out of equilibrium .
faced to the question as to how relevant quantum fluctuations are for the asymptotic real - time dynamics , one would have naively answered that they are irrelevant since quantum mechanics is not expected to play a relevant role in the large time sector , meaning large @xmath242 and for all @xmath413 .
we show in this paper that this would have been a wrong conclusion and that the dynamic behavior is much richer :    quantum effects , together with temperature , determine _ where _ in parameter space a glassy phase exists .
we show that sufficiently strong quantum fluctuations destroy the glassy phase at arbitrary low temperatures driving the system towards a paramagnetic phase .
we predict then the existence of a dynamic quantum phase transition .
this has no analogue in the classical model where glassiness is always present at @xmath7 .
quantum fluctuations also dictate _ when _ , in the glassy phase , non - equilibrium effects manifest .
for any waiting - time @xmath242 after preparation , there exists a characteristic time @xmath376 that determines the end of a stationary regime and the entrance into the aging regime .
this means that for time differences @xmath414 the dynamics is stationary while for @xmath415 two - point functions have explicit waiting - time dependences .
the characteristic time @xmath376 depends on temperature , the strength of quantum fluctuations and the waiting - time @xmath242 .
the dynamics in the paramagnetic phase satisfies the fluctuation - dissipation theorem . in the glassy phase
, we found that the quantum fdt is satisfied in the stationary regime and that a generalization is needed in the aging regime .
this generalization defines a two - time function @xmath364 that is consistent with the definition of an effective temperature @xmath416 and that we later characterized , based on the numerical evidence for this model .
when @xmath358 is different from zero , the modification of the quantum fdt takes a similar form as the classical one but with a @xmath31 and @xmath186 dependent effective temperature .
the characteristic time @xmath376 is reminiscent of a decoherence time in the sense that it separates the stationary regime in which the influence of quantum effects is very strong and explicit from the aging regime in which the influence of quantum effects is implicit , appearing through , e.g. , the @xmath186 dependence of @xmath358 .
one of the virtues of the approach we have followed is that it is close in spirit to the dynamic approach to classical systems .
this encourages us to attempt to extend many of the recent advances in the understanding of classical glassy dynamics to the quantum problem .
we mention below some of the questions that merit further attention .
the model here studied has a modification of the quantum fdt given by a piecewise function @xmath417 if @xmath418 ( @xmath419 ) and a constant @xmath420 if @xmath421 ( @xmath422 ) .
it is known that in classical systems a more general situation can arise .
another family of models exists with the factor measuring the fdt violation given by @xmath417 if @xmath418 and a non - constant function of the auto - correlation @xmath423 if @xmath421 .
we expect the arguments put forward in section [ generalfdt ] to apply to this kind of system with minor modifications .
it would be of particular interest to see how this structure appears in a concrete model as , for instance , the problem of a particle in a long - range correlated random potential .    in the classical case ,
an intriguing connection between the factor that measures the fdt violation@xcite in the non - equilibrium approach @xmath424 and the parisi function @xmath425 , @xmath426 of the replica analysis in equilibrium , has been noticed and checked in several occasions .
the structure of both functions is the same and , on top of this , their explicit values coincide if an argument of `` marginality '' is used to determine @xmath425 in the equilibrium case . in ref .
[ ] , several quantum models related , but not equal , to the one considered here were studied with a gaussian variational method and replica theory in the imaginary time matsubara representation . in order to obtain a physical value for the conductivity of the model , giamarchi and le doussal used the marginality condition to fix the replica breaking point parameter @xmath427 , @xmath428 .
a careful comparison of the dynamic @xmath377 with the static @xmath429 is now in order .
nevertheless , this can not be readily done by comparing the result in ref .
[ ] with ours since the way in which both systems are coupled to the environment are not necessarily equivalent .
thence , if a formal relation between fdt - violations and the replica @xmath425 exists at the quantum level is an interesting question that deserves further study .
the last remark we wish to make about the fdt violations is that in the classical case it has been related to the production of entropy during evolution.@xcite it would be desirable to see whether such a connection also holds at the quantum level .
concerning the formalism here developed , it can be extended in several directions .
on the one hand , one would like to study models with a spatial structure such as the one associated to the motion of a _ manifold _ of internal dimension @xmath430 in a random quenched potential .
this quantum problem has many interesting applications and , in the case of an embedding space with a large dimensionality can be treated , for example , within a dynamic gaussian variational approach.@xcite on the other hand , it would be extremely interesting to study the opposite limiting case of low dimensions and analyze , for example , the real - time dynamics of the ising chain in a transverse field .
in addition , one could investigate the effect induced on the dynamics by different environments like considering sub or super ohmic baths , non - linear couplings and even baths of a completely different nature as those generated by spin variables.@xcite    the martin - siggia rose generating functional for classical stochastic processes takes a particularly simple form if the time - space is enlarged to a superspace by introducing a couple of grassmann coordinates , and the variables of the problem are encoded in a superfield.@xcite this can also be done for systems with disorder.@xcite the dynamic equations are then written in terms of a single super - correlator that encodes the auto - correlation and response .
the equilibrium theorems , viz .
causality , invariance under time - translations and the fluctuation - dissipation theorem , follow from the symmetries of the msr effective action written in this way.@xcite the supersymmetric formulation of the classical dynamics has been helpful to find a solution of mean - field like models with an internal dimension like the problem of the manifold in a random media but , most importantly , it has highlighted the formal connection between the out of equilibrium formalism and the replica approach to the statics .
it would be very interesting to identify a supersymmetric structure in the quantum problem , probably at the level of the keldysh - schwinger generating functional .
this could provide a means to establish a relation between the quantum dynamic formalism here presented and the quantum replica approach .
we believe that an analytical study of aging effects along the lines presented in this paper will be relevant in the understanding of realistic glassy systems .
indeed , even at very low temperatures , where quantum effects are important , a glass is a non - equilibrium system that evolves over many orders of magnitude in time.@xcite aging effects , demonstrating this non - equilibrium evolution , have been observed with , for example , burning hole experiments in organic glasses ( see ref .  [ ] and references therein ) . in the context of magnetic disordered systems with explicit quenched disorder like liho@xmath2y@xmath3f@xmath4
, a great effort has been devoted to the study of the quantum phase transition and the dynamics close to the quantum critical point .
it would be interesting to perform dynamic measurements in the glassy phase of these compounds to search for aging effects and check if a scenario as the one here discussed with , for example , a two step relaxation of correlation functions and fdt violations is also present .
another area of research where the relevance of quantum effects in the very low temperature dynamics is currently under study are systems of magnetic nano - particles such as @xmath431fe@xmath6o@xmath432 particles ( maghemite ) embedded in a silica matrix .
the search for a non - stationary evolution in samples at small concentrations ( @xmath433 ) has given a negative result .
@xcite however , in experiments performed at higher temperatures ( where quantum effects are not expected to play an important role ) with samples with a larger concentration ( @xmath434 ) much of the phenomenology of spin - glasses has been recovered.@xcite it is then licit to ask whether there is a region in the temperature - concentration parameter - space where the systems exhibit aging driven by quantum fluctuations . in the area of magnetic nano - particles , a notion of an effective temperature greater than the temperature of the environment was advocated in ref . [ ] and it would be interesting to explore the connection with effective temperature defined in this paper .
we have chosen to comment on three kind of materials where effects like the ones obtained in this article might be observed .
surely enough , many other glassy systems with or without explicit disorder might exhibit a similar phenomenology .    *
acknowledgments *    we want to especially thank j. kurchan for very useful discussions and constant encouragement during the preparation of this work .
we also wish to thank t. evans , a. georges , t. giamarchi , d. grempel , p. le doussal , r. revers , m. rozenberg and e. vincent for useful discussions . at the early stages of this work
we have also benefited from discussions with l. bettencourt and h. de vega .
the authors want to thank the international center for theoretical physics , trieste , italy for hospitality during the `` summer college in condensed matter on statistical physics of frustrated systems '' .
l.f.c . wishes to thank the institute for theoretical physics , university of california at santa barbara for hospitality during the work - shop `` jamming and rheology '' partially supported by the grant no phy94 - 07194 .
g.l is financially supported by the ec grant n erbfmbict 961226 .
[ appkeldysh ]    the generating functional ( [ generating ] ) is more conveniently expressed with the help of path - integrals : @xmath28   & = & \int d{\bbox \phi }      \ ;   \langle { \bbox \phi }   | t^ * \exp\left ( -\frac{i}{\hbar }   \int_{0}^\infty dt \ , { \bbox \xi}^-(t ) { \bbox \phi}(t ) \right )   t \exp\left ( \frac{i}{\hbar }   \int_{0}^\infty dt \ , { \bbox \xi}^+(t ) { \bbox \phi}(t ) \right ) |     \hat \rho(0 ) | { \bbox \phi } \rangle \nonumber\\ & = & \int d{\bbox \phi}d{\bbox \phi}_1 d{\bbox \phi}_2      \ ;   \langle { \bbox \phi }   |   t^ * \exp\left ( -\frac{i}{\hbar } \int_{0}^\infty dt \ , { \bbox \xi}^-(t ) { \bbox \phi}(t ) \right )    \nonumber\\ & &   \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\ ; \langle { \bbox \phi}_1 |   t \exp\left ( \frac{i}{\hbar } \int_{0}^\infty dt \ , { \bbox \xi}^+(t ) { \bbox \phi}(t ) \right )      \langle { \bbox \phi}_2   | \hat \rho(0 ) | { \bbox \phi } \rangle \ ; .\end{aligned}\ ] ] the two matrix elements are now expressed as usual feynman path - integrals @xmath28    & = &    \int d{\bbox \phi }     d{\bbox \phi}_1     d{\bbox \phi}_2       \int_{{\bbox \phi}^+(0)={\bbox \phi}_2}^{{\bbox \phi}^+(\infty)=    { \bbox \phi}_1 }    { \cal d}{\bbox \phi}^+    \int_{{\bbox \phi}^-(0)={\bbox \phi}}^{{\bbox \phi}^-(\infty)=   { \bbox \phi}_1 }   { \cal d}{\bbox \phi}^-    \nonumber\\ & & \exp\left [ \frac{i}{\hbar } \ , \left ( s[{\bbox \phi^+ ] - s [ \bbox \phi^- } ] + \int dt \ ,   { \bbox \xi}^+ { \bbox \phi}^+ - \int dt \ , { \bbox \xi}^- { \bbox \phi}^- \right ) \right ]   \;\ ; \langle { \bbox \phi}_2   | \hat \rho(0 ) |   { \bbox \phi }   \rangle \ ; .\end{aligned}\ ] ] the first integrals @xmath435 , etc .
are standard while we denote @xmath436 the path - integrals .
the doubling of degrees of freedom @xmath437 is a consequence of having two matrix elements , one for each source @xmath438 .
formally this expression is rewritten as @xmath28    & = &    \int { \cal d}{\bbox \phi}^+    \int { \cal d}{\bbox \phi}^-    \exp\left [ \frac{i}{\hbar } \ , \left ( s[{\bbox \phi^+ ] - s [ \bbox \phi^- } ] + \int dt \ ,   { \bbox \xi}^+ { \bbox \phi}^+ - \int dt \ , { \bbox \xi}^- { \bbox \phi}^- \right ) \right ]   \nonumber\\ & &   \;\;\;\;\;\;\;\;\;\;\ ; \times   \langle { \bbox \phi}^+ |
\hat \rho(0 ) |   { \bbox \phi}^- \rangle \ ; .\end{aligned}\ ] ]
[ appfdt ]    in this appendix we recall the quantum fluctuation - dissipation theorem . proofs and descriptions of this theorem can be found in several textbooks.@xcite we express it in the time - domain in a form that we use in section [ equilibriumdynamics ] to show that the dynamic equations are compatible with a tti - fdt ansatz .
we also write it in a mixed time - fourier notation that gives us insight as to how extend it to the case of glassy non - equilibrium dynamics .
if at time @xmath211 the system is characterized by a density functional @xmath439 , the two - time correlation functions read @xmath440 \label{corr}\ ] ] where the time - dependent operators , in the heisenberg representation , are defined as @xmath441 the trace is defined in the usual way , @xmath442 \equiv \sum_\alpha \langle \psi_\alpha | \bullet |   \psi_\alpha \rangle$ ] , with @xmath443 an orthonormal basis in fock space .
the normalization is given by @xmath444 $ ] .
since operators do not commute , the quantum two - time auto - correlation functions are not symmetric in times @xmath445 one thus defines the symmetrized and antisymmetrized correlation functions : @xmath446}(t , t ' ) & = &   \frac12 \ ,   \langle a(t ) b(t ' ) - b(t ' ) a(t ) \rangle \ ; , \end{aligned}\ ] ] respectively .
in linear response theory @xmath447 ( see eq .
( [ respdef ] ) for the definition ) and the correlation @xmath448}(t , t')$ ] are related by the kubo formula @xmath449 \ ,
\rangle = \frac{2 i}{\hbar } \ ; \theta(t - t ' ) \ ,   c_{[a , b]}(t , t ' ) \label{fdtasym } \ ; .\ ] ]    if the system has reached equilibrium with a heat - bath at temperature @xmath31 at time @xmath211 , the density functional @xmath439 is just the boltzmann factor @xmath450 .
it is then immediate to show that , in equilibrium , time - translation invariance holds @xmath451 in addition , eq .
( [ corr ] ) with @xmath452 imply the kms condition @xmath453 using now the kms properties and assuming , for definiteness , that @xmath454 it is easy to verify the following equation @xmath455 where @xmath456 .
this is a way to express fdt through an analytic continuation to complex times that we use in section  [ equilibriumdynamics ] to show that a tti and fdt solution holds in the paramagnetic phase .    in terms of the fourier
transformed @xmath457 defined by @xmath458 the kms relations read @xmath459 and lead to the following relation between fourier transforms of the correlation functions : @xmath460}(\omega ) & = & \left(1 - \exp(-\beta \hbar \omega ) \right ) \
,   c_{ab}(\omega )   \ ; , \nonumber\\ 2 c_{\{a , b \}}(\omega ) & = & \left(1 + \exp(-\beta \hbar \omega ) \right ) \ ,   c_{ab}(\omega ) \ ; ,   \nonumber\\ c_{[a , b]}(\omega ) & = & \tanh\left ( \frac{\beta \hbar \omega}{2 } \right ) \ ,   c_{\{a , b \}}(\omega ) \ ; .\end{aligned}\ ] ] back in eq .
( [ fdtasym ] ) this implies fdt @xmath461 using @xmath462 one has @xmath463 from which we obtain the real and imaginary relations @xmath464 if @xmath465 , @xmath466 and eq .
( [ fdtfourier ] ) becomes the classical fdt : @xmath467 with @xmath468 .
in the zero temperature limit , and for @xmath469 , fdt reads @xmath470 a finite integrated response @xmath471 implies @xmath472 . in the glassy model
we study in the paper , below @xmath283 , the response and correlation turn out to be the sum of two contributions .
one is a stationary part that satisfies fdt .
since one expects it to yield a finite contribution to the susceptibility @xmath264 , one must have : @xmath473 thus , the stationary part of the correlation function oscillates around zero or , in other words , the full correlation oscillates around @xmath218 at zero temperature .
this is observed in the numerical solution to the dynamic equations .
[ classicallimit ]    the classical limit can be checked at every stage of the calculation .
we shall here show how one can recover ( i ) the martin - siggia - rose ( msr ) action from the keldysh - shwinger ( ks ) action ( [ action1 ] ) ; ( ii ) the equations of motion , for the classical correlation and response from the equations of motion for their quantum partners .
this appendix is intended for readers already familiar with the msr formalism .
we refer to ref .
[ ] for a detailed description of the method .
let us identify , in the limit @xmath474 , the quantum linear combinations on the left with the classical variables on the right @xmath475    in the classical limit @xmath476 , and when the cut - off tends to infinity , the kernels @xmath74 and @xmath73 become @xmath477 with these identifications , we establish a connection between the ks and the msr actions .
the friction and thermal terms become @xmath478 by calling @xmath96 we recover two terms of the msr action .
the kinetic and constraint terms reduce to @xmath479 by calling @xmath480 and @xmath481 the last two terms impose the spherical constraint and are usually included in the path - integral measure , while the second term is the one left in the msr action : @xmath482    finally , the terms that depend upon disorder become @xmath483 - v[{\bbox \sigma^- } , j ]    =   - \int_0^\infty dt \ ; i\hat{\bbox s } \ , \frac{\delta v[{\bbox s } , j ] } { \delta { \bbox s } }   \ ; .\ ] ]    putting everything together , in the large friction limit @xmath484 such that one neglects the second derivative terms in the kinetic part of the action , one recovers the usual msr action for the langevin dynamics of the classical model in a thermal bath at temperature @xmath31 : @xmath485 } { \delta { \bbox s}(t ) }   i \hat{\bbox s}(t )    \right ] \ ; .\ ] ]    in order to study the classical limit of the dynamic equations ( [ schwingerr ] ) and ( [ schwingerc ] ) we recognize the classical limit of the correlation and response as follows : @xmath486 in a similar way , it is now easy to verify that the classical dynamic equation for the @xmath0-spin model follow from the @xmath487 limit of eqs .
( [ schwingerr ] ) and ( [ schwingerc ] ) .
[ integrals ]    integrals over time windows spanning the interval @xmath488 $ ] appear in the dynamic equations ( [ schwingerr ] ) and ( [ schwingerc ] ) . in the large times limit we approximate them in the way we here describe .
integrals of the form @xmath489 appear , for example , in @xmath123 .
the idea is to separate the integration time - interval as @xmath490 = [ 0,\delta ] \;\ ; u \;\ ; [ \delta , t^- ] \;\ ;   u \;\ ; [ t^-,t ] \ ; .\ ] ] if @xmath249 is chosen to be a finite time , all functions can be approximated by @xmath491 that vanishes when @xmath492 . since the integration interval is finite , this term can be neglected . in the second interval
the functions vary in the aging regime and in the third interval they vary in the stationary regime .
thus @xmath493 \ ; . \nonumber\\\end{aligned}\ ] ] we assume that this separation is sharp and that we can neglect the corrections associated to mixing of the three regimes . in the third term we replaced @xmath299 and @xmath300 in terms of @xmath494 , @xmath495 .
since in all cases either @xmath299 or @xmath300 is proportional to the response , we can neglect the first term using the property ( [ finitetimes ] ) of the weak long - term memory scenario . we can then replace the lower limit of the second integral by @xmath219 and the upper limit of the second integral by @xmath193 .
in addition , assuming that @xmath300 is proportional to the response , @xmath496 this yields @xmath497      another type of integrals is : is @xmath498 in particular , if @xmath499 , @xmath500 and @xmath501 , this integrals yields the susceptibility .
let us assume that @xmath193 and @xmath211 are far apart ; we start as for the first type of integrals , by dividing the time interval in three subintervals @xmath502 = [ t',t'^+ ] \;\ ; u \;\ ; [ t'^+,t^- ] \;\ ; u \;\;[t^-,t]\ ] ] and by approximating the functions by their functional form inside each of the intervals : @xmath503 in the first term @xmath504 can be replaced by @xmath505 .
the same applies to @xmath299 in the last term .
all functions vary fast in the stationary regime but very slowly in the aging regime .
the next assumption is that functions in the aging regime that are convoluted with functions in the stationary regime , can be considered to be constant and taken out of the integral .
that is to say @xmath506 where we used @xmath507 and @xmath508 .
typically , @xmath509 and @xmath510 are proportional to the response function . using fdt the integrals in the first and third term
can then be performed in the classical limit or they can be expressed as functions of the correlation in the quantum case .
the second term instead depends exclusively on the aging dynamic sector .
e - mail address : leticia@physique.ens.fr e - mail address : lozano@ipno.in2p3.fr unit propre du cnrs , associe   lecole normale suprieure et   luniversit de paris sud .
unit de recherche des universits paris 11 et 6 associe au cnrs .
l. lundgren , p. svedlindh , p. nordblad and o. beckmann , phys .
lett . * 51 * , 911 ( 1983 ) .
e. vincent , j. hammann , m. ocio , j.p .
bouchaud and l. f. cugliandolo , _ slow dynamics and aging _ in sitges conference on glassy systems , m. rub ed .
( springer - verlag , 1997 ) , cond - mat/9607224 .
f. alberici , p. doussineau and a. levelut , j. phys .
i france * 7 * , 329 ( 1997 ) ; europhys . lett . * 39 * 329 ( 1997 ) .
f. alberici - kious , j - p bouchaud , l. f. cugliandolo , p. doussineau and a. levelut , _ aging in _
ta@xmath432 : _ a domain growth interpretation _ , cond - mat/9805208
.        a. j. bray and m. a. moore , j. phys .
* c17 * , l463 ( 1984 ) and in _ heidelberg colloquium on glassy dynamics _ , lecture notes in physics * 275 * , ed .
j. l. van hemmen and i. morgenstern ( springer - verlag , berlin ) .
g.j . koper and h.j .
hilhorst , j. phys .
( france ) * 49 * , 429 ( 1988 ) .
d.s . fisher and d.a .
huse , phys .
lett . * 56 * , 1601 ( 1986 ) ; phys . rev . *
b38 * , 373 ( 1988 ) .
v. dotsenko , j. phys .
* c18 * , 6023 ( 1985 ) . v. dotsenko , m. feigelman and l. b. ioffe , _ spin - glasses and related problems _ , soviet scientific reviews * 15 * ( harwood 1990 ) . m. v. feigelman and v. vinokur , j. phys .
i ( france ) * 49 * , 1731 ( 1988 ) .
g. parisi , phys .
lett . * 79 * , 3660 ( 1997 ) .
w. kob and j - l barrat , phys .
* 78 * , 4581 ( 1997 ) .
a. barrat , phys . rev . *
e57 * , 3629 ( 1998 ) .
h. yoshino , _ aging dynamics of an elastic string diffusing in a disordered media _ , cond - mat/9802283 .
j - l barrat and w. kob , _ fluctuation dissipation ratio in an aging lennard - jones glass _ , cond - mat/9806027 .
t. d. kirkpatrick and d. thirumalai , phys .
* b36 * , 5388 ( 1987 ) . t. r. kirkpatrick and p. wolynes , phys . rev . * a35 * , 3072 ( 1987 ) ; phys
* b36 * , 8552 ( 1987 ) .
t. r. kirkpatrick , d. thirumalai , p. g. wolynes ; phys . rev . *
a40 * , 1045 ( 1989 ) .
j - p bouchaud , l. f. cugliandolo , j. kurchan and m. mzard , _ out of equilibrium dynamics in spin - glasses and other glassy systems _ , in  spin - glasses and random fields , a. p. young ed .
( world scientific , singapore , 1997 ) , cond - mat/9702070 .
w. wu , b. ellmann , t. f. rosenbaum , g. aeppli and d. h. reich , phys .
* 67 * , 2076 ( 1991 ) .
w. wu , d. bitko , t. f. rosenbaum and g. aeppli , phys .
lett . * 71 * , 1919 ( 1993 ) . t. f. rosenbaum , j. phys . *
c8 * , 9759 ( 1996 ) .
e. courtens , j. phys.lett .
( paris ) * 43 * , l199 ( 1982 ) ; phys .
lett . * 52 * , 69 ( 1984 ) .
e. matsushita , t. matsubara , prog .
. phys . * 71 * , 235 ( 1984 ) .
r. pirc , b. tadic and r. blinc , z. phys . *
b61 * , 69 ( 1985 ) ; phys . rev .
* b36 * , 8607 ( 1987 ) .
p. shukla and s. singh , phys . rev . * b23 * , 4661 ( 1981 ) .
a. j. bray and m. a. moore , * 13 * l665 ( 1985 ) .
k. d. usadel , solid state commun . * 58 * , 629 ( 1986 ) .
t. yamamoto and h. ishii , j. phys .
* c20 * , 6053 ( 1987 ) . y. y. godschmidt and p. y. lai , phys .
64 * , 2467 ( 1990 ) .
y. y. godschmidt , phys . rev . *
b41 * , 4858 . v dobrosavljevic and d thirumalai , j. phys . *
a23 * , l767 ( 1990 ) . l. de cesare , k. lubierska - walasek , i. rabuffo and k. walasek , physica * a214 * , 499 ( 1995 ) ; j. phys . * a29 * , 1605 ( 1996 ) . t. k. kope , phys . rev . * b54 * , 3367 ( 1996 ) .
h. rieger and a. p. young , _ quantum spin - glasses _ ( springer - verlag , berlin , 1996 ) , cond - mat/9607005 .
r. n. bhatt , _ quantum spin - glasses _ , in  spin - glasses and random fields , a. p. young ed .
( world scientific , singapore , 1997 ) .
c. p. martin , e. siggia and h. a. rose , phys .
* a8 * , 423 ( 1973 ) . h. k. janssen , z. phys . *
b23 * , 377 ( 1976 ) ; _ dynamics of critical phenomena and related topics _ , lecture notes in physics * 104 * , c. p. enz ed . , ( springer - verlag , berlin , 1979 ) .
the relevance of the quantum langevin equation has been extensively treated in the literature .
some references are g. w. ford , m. kac and m. mazur , j. math . phys .
* 6 * , 504 ( 1965 ) .
g. w. ford and m. kac , j. stat .
phys * 46 * , 803 ( 1987 ) .
t. jonsson , j. mattsson , c. djurberg , f. a. kahn , p. nordblad and p. svedlindh , phys .
lett . * 75 * , 4138 ( 1995 ) .
t. jonsson , j. mattsson , p. nordblad , p. svedlindh , j. magn .
magn . mat . *
168 * , 269 ( 1997 )