glassy systems like spin - glasses and structural glasses below the glass transition temperatures are out of equilibrium even on the macroscopic time scale .
thus , the slow dynamics of glassy systems has been a subject of continuous interest in the past years @xcite .
experimentally , the anomalous dynamical behaviors are characterized by slow relaxation with long time tail and aging phenomena .    among them ,
aging phenomena are the most striking dynamical behaviors as follows : two - time quantities like the correlation functions @xmath2 explicitly depend on the time elapsed after the quench @xmath3 ( the waiting time ) .
if the waiting time @xmath3 is of microscopic time scale , the phenomena are merely transient on the way of relaxation to equilibrium .
however , dependence on the waiting time continues even when @xmath3 is so large that one - time quantities like the magnetization are asymptotically close to time - independent values @xcite .
since these phenomena mean that the dynamics is not stationary , aging is a sign showing that these systems are out of equilibrium even in macroscopic time scale , i.e. several days or weeks .
these aging phenomena appear also in mean - field models of spin - glasses and do not disappear even in the infinite waiting time limit @xcite .
in addition , the phase space of these models decomposes into a large number of areas separated with infinitely high free energy barriers .
thus , glassy systems never reach true equilibrium and hence they are beyond the scope of thermodynamics and equilibrium statistical mechanics . hence , in order to describe thermodynamic properties of glassy systems , out - of - equilibrium thermodynamics based on a dynamical description without the assumption of ergodicity is necessary .    in investigation of such anomalous dynamical behavior of glassy systems
, it was found that aging phenomena have some universal properties .
theoretical analyses have suggested that there are two time regimes characterized by different dynamical properties of the two - time correlation function @xmath2 and the associated linear response function @xmath4 @xcite . in the first time regime where the time difference @xmath5 is short compared to @xmath3 , the dynamics looks stationary and the usual fluctuation - dissipation theorem ( fdt ) holds . on the other hand , in the second time regime where the time difference is comparable to @xmath3 , aging phenomena occur ,
i.e. @xmath2 depends on @xmath3 apart from dependence on @xmath5 .
in addition , it is known that the usual fdt between the correlation @xmath2 and the response function @xmath4 should be modified in a well - defined way which involves the rescaling of the temperature @xcite .
the modification was found to be valid not only for mean - field models but also for other glassy systems : spin - glass models with finite - range interactions @xcite , real spin - glasses @xcite , structural glasses @xcite and a model of phase separation @xcite .
in addition , it is known for some glassy systems that the correlation function @xmath2 obeys the scaling law that it depends on @xmath6 and @xmath3 only through the value of @xmath7 , where @xmath8 is a system - dependent increasing function of time .
aging of the correlation function and the modification of fdt imply that properties of the work done by modulating an external field in an isothermal process are completely different from properties predicted by traditional thermodynamics .
in addition , the existence of the two time scales implies that thermodynamic properties must strongly depend on experimental time scales .
hence , the experimental time scale dependence of the thermodynamic properties of glassy systems should be investigated to construct out - of - equilibrium thermodynamics for glassy systems .    in order to describe our results precisely , we summarize thermodynamics for an isothermal process .
thermodynamics tells that when one quasistatically changes the external field the work needed for the change is independent of a path of changing the external field .
in addition , the quasistatic work is equal to the change of the helmholtz free energy .
when the process is not quasistatic , the work is larger than the quasistatic work .
this fact is called the principle of the minimum work and is derived from the second law of thermodynamics .
we show that the properties described above do not hold in systems with aging .
more precisely , there are two experimental time regimes characterized by different thermal properties .
the first regime is a time - domain where the waiting time is much longer than the time lapse of the process .
we call the time lapse the measurement time . in this regime ,
the principle of the minimum work holds even though a system is out of equilibrium .
more precisely , when the process is not infinitely slow , the work needed for the process is larger than the work for an infinitely slow process .
in addition , value of the work for the infinitely slow process depends only on the initial state and the final state and hence it can play a role of a free energy .
the second regime is the experimental time regime where the length of the measurement time are comparable to that of the waiting time . in this regime , for the single - correlation - scale systems such as @xmath0-spin spherical spin - glasses the work done in an infinitely slow operation depends on the path of changing the field even when the waiting time is infinite .
this property form a striking contrast to the consequence of traditional thermodynamics described above . on the other hand , for the multi - correlation - scale systems such as sherrington - kirkpatrick model ,
the work done in an infinitely slow operation is independent of the path .    in sec .
[ section2 ] , we describe an isothermal process considered in this paper and introduce the two time regimes which characterize the experimental time scales and play a significant role in this paper . in sec .
[ section3 ] , we see that in the first time regime , usual thermodynamic relations hold even though the system is out of equilibrium .
the only difference is that the value of the work in an infinitely slow operation @xmath9 is different from that of the change of the free energy calculated from equilibrium statistical mechanics . in sec .
[ section4 ] , we present general discussion on properties of the work @xmath9 in an infinitely slow operation in the second time regime and derive conditions when @xmath9 depends on the path of changing the external field . in sec .
[ section5 ] , by using the results obtained in the previous section , we show that for the single - correlation - scale systems the work in an infinitely slow operation depends on the path of changing the external field as a consequence of aging .
possibility of observation of this path - dependence is also discussed .
in addition , we show that for the multi - correlation - scale systems @xmath9 is not path - dependent .
our results are summarized in sec .
[ section6 ] , where implications of our results on thermodynamics of glassy systems and experimental protocols to observe quasiequilibrium properties are discussed .
we describe an isothermal process to consider thermodynamic properties when aging occurs . a simple way to observe aging phenomena is through the following field cooling process ( fig .
[ protocol ] ) ; the temperature of the heat bath is decreased in a small field @xmath10 to a sub - critical temperature at time 0 . after a waiting time @xmath11 the field begins to change according to a given time - dependence @xmath12 .
this change of the field continues for a period of @xmath13 .
the initial and final values of @xmath14 are @xmath15 and @xmath16 .
we refer to @xmath13 , the time lapse of the change of the field , as the measurement time . the waiting time
@xmath11 and the measurement time @xmath13 characterize the experimental time scales .
when the field is so weak that the response is linear , the work @xmath17 done on the sample during the process is given in terms of the response function @xmath4 by @xmath18 where the `` magnetization '' @xmath19 is given by @xmath20 since we are interested in long waiting time behavior , the contribution of the first term of eq .
( [ m ] ) is ignored . for feasibility of showing long time behavior
, we rewrite eq .
( [ m ] ) by integration by part to express @xmath19 in terms of a susceptibility @xmath21 instead of the response function ; @xmath22 where the susceptibility is defined as @xmath23    in order to discuss the dependence of the work on the waiting time and the measurement time , we rewrite the expression of the work by the transformations : @xmath24 , @xmath25 .
we assume that @xmath26 and @xmath27 are finite , in order to exclude the unrealistic cases where the speed and the acceleration of changing the field is infinity .
for example , the path such as @xmath28 is excluded , since @xmath26 and @xmath27 are infinite at @xmath29 .
thus , the work is reduced to @xmath30 which implies that the dependence of the work on the experimental time scales is determined by that of the susceptibility . since the first term of eq .
( [ w_s ] ) becomes a constant @xmath31 in the long waiting time limit which we are interested in , we will analyze properties of the second term of eq .
( [ w_s ] ) , i.e. the work done in the zero - field cooling process ( @xmath32 ) , in the rest of this chapter for simplicity ; @xmath33    in order to discuss behavior of the susceptibility when aging occurs , we recapitulate the long - time behavior of the response function referred to in the previous section .
it is known that there are two time regimes characterized by different behavior of the correlation function @xmath2 and the fdt @xcite ;    * at long times @xmath6 and @xmath3 such that @xmath34 , the correlation function is the function of only the time difference @xmath5 , i.e. the time - translational invariance ( tti ) holds .
in addition , although the sample is out of equilibrium , the usual fdt holds @xcite as @xmath35 where @xmath36 is defined as @xmath37 .
since the properties of two - time quantities in this time regime are the same as that in equilibrium , this time regime is called the quasiequilibrium regime . * whereas , at long and well - separated times such that @xmath38 , aging occurs , i.e. the two - time correlation function @xmath2 depends on @xmath3 even in the long time limit @xmath39 .
this implies no tti and this time regime is called the aging regime .
in addition , it is known that in the aging regime the fdt is modified as @xmath40}{k_b t}\frac{\partial c(t , t')}{\partial t ' } , \label{generalized_fdt}\ ] ] where the fdt violation factor @xmath41 is a function which depends on @xmath42 only through the dependence of the correlation function @xmath2 @xcite .
thus , a system can not be considered to be in a quasiequilibrium state , since the usual fdt is strongly violated .
when only one correlation scale exists apart from the quasiequilibrium regime , it is known that the fdt violation factor is a constant .    in order to clarify the meaning of the time region of the aging regime ( @xmath38 ) , we give an explicit expression of aging of the correlation function as @xmath43 \equiv   \hat{c}(\mu ) , \label{explicit_expression_aging}\ ] ] where @xmath44 is a system - dependent function which characterizes the aging regime and @xmath45 is the time difference @xmath5 .
the correlation function @xmath46 depends on the waiting time @xmath3 through the value of @xmath47 .
occurrence of aging means that the limiting function @xmath46 takes a non - trivial value such that @xmath48 . here
, @xmath49 is the dynamical e - a order parameter .
thus , the time region of the aging regime is the time region where @xmath50 is finite . here
, we assume that @xmath47 and @xmath51 are positive .
this definition of aging is illustrated in terms of contour plot of the correlation function on @xmath52 plain ( fig .
[ contour ] ( a ) ) .
the plot when aging occurs is completely different from that when aging does not occur ( fig .
[ contour ] ( b ) ) .
the contour lines give the system dependent function @xmath44 which characterizes the aging regime , since the contour lines are given by @xmath53 in the contour plot .
we show two examples of the function @xmath44 and the aging regime for systems like spherical spin - glasses and real spin - glasses . in these systems
, it is known that the correlation function behaves as a function of @xmath7 , where the scaling function @xmath8 is a system - dependent function @xcite .
the function @xmath44 is given in terms of @xmath8 as @xmath54-t$ ] , where @xmath55 is a constant larger than unity .
the inverse function of @xmath8 exists since @xmath8 is a monotonically increasing function .
for example ,    1 .   when @xmath56 , @xmath57 thus , @xmath58 is a function of @xmath47 and hence the aging regime is the time regime where @xmath59 is finite .
when @xmath60 ( 0<\alpha < 1)$ ] , @xmath61 thus , @xmath62 is a function of @xmath47 and hence the aging regime is the time regime when @xmath63 is finite .
from these results , it is shown by the definition of the susceptibility eq .
( [ chi ] ) that the behavior of the susceptibility depends on the time regimes according to the dependence of the behavior of the response function . in the quasiequilibrium regime ,
the susceptibility depends only on the time difference @xmath64 and is given by the correlation function as @xmath65 \label{chi_f}\ ] ] where @xmath66 . whereas , in the aging regime , time translational invariance ( tti ) does not hold and the susceptibility is given by the correlation function as @xmath67 where @xmath41 is the fdt violation factor defined by eq .
( [ generalized_fdt ] ) .
it implies that the susceptibility is a function of the correlation function .    since the difference of the two arguments of the susceptibility in eq .
( [ w_s_2 ] ) is @xmath68 , there are also two time regimes characterized by different dependence of the work on the experimental time scales @xmath13 and @xmath69 .
* when the waiting time is much longer than the measurement time ( @xmath70 ) , the susceptibility which appears in eq .
( [ w_s_2 ] ) obeys eq .
( [ chi_f ] ) since @xmath71 holds .
hence , from eqs .
( [ w_s_2 ] ) and ( [ chi_f ] ) the work in this regime is written as @xmath72 \label{w_f}.\end{aligned}\ ] ] in the derivation , @xmath73 for any path @xmath74 is used .
* on the other hand , when both the waiting time and the measurement time are long ( @xmath75 ) , the susceptibility obeys eq .
( [ chi_a ] ) since @xmath76 holds @xcite . in this regime
, it is seen from eq .
( [ chi_a ] ) that the work is a functional of the correlation function which shows aging .
the ambiguous relation @xmath77 is described explicitly in sec .
[ section4 ] with the function @xmath44 defined in eq .
( [ explicit_expression_aging ] ) .
we refer to the former case as the quasiequilibrium regime and the latter case as the aging regime without any confusion with the time regimes characterized by the behavior of the correlation and the response function .
in this section , we investigate the properties of the work done on the sample when the measurement time @xmath13 is short compared to the long waiting time @xmath11 ( the quasiequilibrium regime ) .      first , we derive the work done in an slow operation such that the measurement time @xmath13 is so long that the correlation function relax to a time - independent value .
this process is formulated by taking the infinite measurement time limit @xmath78 .
glassy systems are out of equilibrium even in such an infinitely slow process .
hence , we call a process where the external field is changed infinitely slowly a `` slow process '' instead of a `` quasistatic process '' .    since for the slow process
the time lapse of the process is infinity , the work done in the slow process in the quasiequilibrium regime is given by taking the limit @xmath78 after taking the infinite waiting time limit of eq .
( [ w_f ] ) .
since the integrand , the correlation function , is finite , the order of limit and integration can be changed .
thus , we can write the work @xmath9 for the slow process as @xmath79.\end{aligned}\ ] ]    here , in order to give the expression of the work @xmath9 , we introduce the dynamical edwards - anderson ( e - a ) order parameter defined as @xmath80 using eq .
( [ integration ] ) and the above definition , the work @xmath81 for the slow process is given by @xmath82 \label{w_qs_f } \label{quasistatic_work_f}.\ ] ] therefore , the work for the slow process in the quasiequilibrium regime is the difference of a state function , since the right hand side of eq .
( [ w_qs_f ] ) is independent of the path of changing the field @xmath74 and dependent only on the thermodynamic variables ( @xmath83 and @xmath84 ) and constants intrinsic to the system ( @xmath85 and @xmath49 ) .
it is important to note that this property is derived by using only the fdt .    on the other hand , for usual systems apart from glassy systems
, one may expect that equilibrium is usually achieved when @xmath69 is long ( @xmath86 ) .
so , thermodynamics can be applied and the work for the slow process is equal to the change of the helmholtz free energy , @xmath87 , calculated by equilibrium statistical mechanics .
however , we show below that this naive expectation fails for glassy systems . more precisely , the work for the slow process , eq .
( [ quasistatic_work_f ] ) , is different from the change of the helmholtz free energy calculated by equilibrium statistical mechanics .
it is because the glassy systems are out of equilibrium even when the waiting time is infinite .    assuming that @xmath88 is a physical quantity coupled to the external field @xmath14 linearly , we see from statistical mechanics with the assumption of ergodicity that the isothermal susceptibility @xmath89 is given by @xmath90 where @xmath91 denotes average over the gibbs - boltzmann distribution and @xmath92 denotes disorder - average . thus , the free energy difference @xmath87 due to change @xmath93 of the external field is @xmath94 where @xmath95 is the usual e - a order parameter .
since ergodicity is broken for glassy systems , the phase space decomposes into many pure states . assuming that @xmath96 denotes the thermal average in local equilibrium in the pure state @xmath97 and @xmath98 denotes the free energy at the pure state @xmath97
, we see that @xmath99 where @xmath100 is the probability that the system is found at the pure state @xmath97 in true equilibrium and is given by @xmath101 .    in order to compare @xmath9 , eq .
( [ w_qs_f ] ) , with @xmath87 , eq .
( [ free_energy_difference ] ) , we rewrite eq .
( [ w_qs_f ] ) in terms of pure states . since the correlation function is given by @xmath102 and the local equilibrium in pure states
is achieved in the long time limit , @xmath103 where @xmath104 is the probability that the system is found in pure state @xmath97 at time @xmath105 . since the system is out of equilibrium at time @xmath105 , @xmath106 .
thus , @xmath107 is not equal to @xmath85 . in addition , from eqs .
( [ ea_parameter_pure_state ] ) and ( [ q_pure_state ] ) , @xmath108 is not equal to @xmath49 .
consequently , we conclude that the work @xmath9 for the slow process does not coincide with the free energy difference calculated by statistical mechanics with assumption of ergodicity .
we discuss the properties of the work when the measurement time is finite . from eqs .
( [ w_f ] ) and ( [ quasistatic_work_f ] ) , the difference between the work when the process is not slow and the work for the slow process is given by @xmath109-q \ } \label{devi_work}.\ ] ] @xmath110>q$ ] when @xmath13 is finite , since the quasiequilibrium regime is considered .
hence , the difference from the work for the slow process is positive for any path of changing the field , when the measurement time is finite .
this implies the principle of minimum work ; @xmath111 where the equality holds only when the measurement time is infinite , i.e. the slow - process limit .
our result is very similar to the consequence of thermodynamics which tells the work for non - quasistatic process is larger than the change of the helmholtz free energy .
however , our result is different from that and beyond the scope of thermodynamics since in our discussion the initial state and the final state of the process are out of equilibrium and the value of the work for the slow process is different from the value of change in the helmholtz free energy derived by equilibrium statistical mechanics .
the long measurement time behavior of @xmath112 is analyzed in this subsection .
we show that the behavior of @xmath112 when the measurement time @xmath13 is long but finite is determined by the long time behavior of the correlation function .
we describe below the results for four types of behavior of the correlation function which include almost all types of relaxation of the correlation , e.g. , the exponential , the power law @xcite , the logarithmic @xcite and the stretched exponential relaxation @xcite .
the derivations are given in appendix  [ appendix ] .    1 .
when the correlation function behaves as @xmath113=0 $ ] , the difference is given by @xmath114 where @xmath115 \label{k}\ ] ] and @xmath81 denotes the work for the slow process .
this case includes the power law relaxation such that @xmath116 when @xmath117 and the stretched exponential relaxation ( @xmath118 ) as well as the exponential relaxation @xcite .
when the correlation function behaves as @xmath113=c$ ] , the difference is given by @xmath119 these two results show that @xmath120 is proportional to the inverse of the measurement time only when @xmath113=0 $ ] . in these two cases ,
the difference from the work for the slow process takes the minimum value when the field increases linearly as @xmath121 , since we assume that @xmath26 and @xmath27 is finite .
when the correlation function obeys the power law relaxation , such that @xmath116 ( @xmath122 ) , the difference is given by @xmath123 in this case , the difference from @xmath9 obeys the power law whose exponent is equal to the exponent of the correlation function .
when the correlation function obeys the logarithmic relaxation @xmath124 , the difference also obeys the logarithmic relaxation as @xmath125 \label{log}.\ ] ] in this case , the difference does not depend on the path @xmath74 .
experimentally , these results tell how long the measurement should take and the suitable path of changing the external field in order to determine the value of the work for the slow process , i.e. the difference of a state function in quasiequilibrium regime .
in this section , we discuss the properties of the work in the isothermal slow process when both the waiting time and the measurement time are long , i.e. in the aging regime , by using the modified fdt and aging of the correlation function introduced in eqs .
( [ generalized_fdt ] ) and ( [ explicit_expression_aging ] ) .
condition when the work for the slow process depends on the path changing the external field is obtained . by using the results obtained in this section ,
the path - dependence of the work for particular systems is discussed in the next section , sec .
[ section5 ] .      as shown in sec .
[ section2 ] , in the aging regime the work is a functional of the correlation function which shows aging . the slow process of the aging regime is given by taking the long measurement time limit @xmath78 with holding the relation @xmath77 , which guarantees that the slow - process limit is taken within the aging regime .
the slow - process limit has always this meaning in this section .
the order of the slow - process limit and the integration can be changed since the susceptibility @xmath126 is finite . from eq .
( [ w_s_2 ] ) the work @xmath81 for the slow process in the aging regime is given by @xmath127 \label{qs_w_a}. \end{aligned}\ ] ] thus , the dependence of @xmath128 on @xmath129 and @xmath130 determines the dependence of the work for the slow process on the path of changing the field @xmath74 ; @xmath9 is independent of @xmath74 if @xmath131 does not depend on @xmath129 and @xmath130 .
on the other hand , @xmath9 is a functional of @xmath74 if @xmath131 depends on @xmath129 and @xmath130 .
we derive conditions when @xmath9 depends on @xmath74 in the rest of this section .    in order to analyze @xmath128 , we clarify the meaning of @xmath77 or the aging regime by using @xmath51 in eq .
( [ explicit_expression_aging ] ) . by substitution @xmath132 in eq .
( [ explicit_expression_aging ] ) , we see that the following equation holds in the region except for the point @xmath133 which does not contribute to the value of the work by itself ; @xmath134 } } }   c(s_1 \triangle t+t_w , s_2\triangle t+t_w ) = \hat{c}(\mu ) .
\label{2}\ ] ] this equation implies that the aging regime is the experimental time regime where @xmath135 $ ] is finite and the slow - process limit of the correlation function in the integrand in eq .
( [ qs_w_a ] ) is a function of @xmath136 $ ] .
we discuss below the aging regimes and dependence of @xmath131 on @xmath129 and @xmath130 in four cases of the different long - time behaviors of @xmath44 which exhaust all possibilities ( see fig .
[ g ] ) .    1 .
when @xmath137 , @xmath50 is finite only in the region @xmath138 . in the region of integration of eq .
( [ qs_w_a ] ) except for the point @xmath139 , @xmath140\rightarrow\infty \label{aging_infty}\ ] ] as @xmath141 since @xmath142 .
an example of such @xmath51 is @xmath143 . when @xmath143 and @xmath144 , @xmath145=\infty$ ] as shown in eq .
( [ aging_infty ] ) .
since the point @xmath139 does not contribute to the value of the work by itself , the experimental time regime where @xmath146 $ ] is finite does not exist in this case .
if there is no other aging regime except for that characterized by such @xmath51 exists , only the quasiequilibrium regime contributes to the value of the work and hence the work for the slow process is independent of @xmath74 and is given by eq .
( [ quasistatic_work_f ] ) .
2 .   when @xmath147 , @xmath148 } =   \frac{s_2 \triangle t+t_w}{(s_1-s_2)\triangle t } =   \frac{s_2+t_w/\triangle t}{s_1-s_2}.\ ] ] thus , from eq .
( [ 2 ] ) , @xmath149 is a function of @xmath150 ; @xmath151 } } } c(s_1\triangle t+t_w , s_2\triangle t+t_w )   \nonumber \\   & = & \hat{c}(\frac{s_2+\mu'}{s_1-s_2 } ) .
\label{t}\end{aligned}\ ] ] it implies that the aging regime is the region where @xmath152 is finite . in this aging regime the slow - process limit of the correlation function in the integrand of eq .
( [ qs_w_a ] ) depends on @xmath129 and @xmath130 if @xmath153 in eq .
( [ 2 ] ) is not a constant .
it means that the work for the slow process , @xmath9 , depends on the path of change of the external field , @xmath74 .
3 .   when @xmath154 and @xmath155 exists , @xmath156}= \frac{s_2 \triangle t+t_w}{g(\triangle t)}\frac{g(\triangle t ) } { g[(s_1-s_2)\triangle t ] } \nonumber \\
\rightarrow \frac{t_w}{g(\triangle t)}g_r(s_1-s_2 ) , \label{aging_regime}\end{aligned}\ ] ] as @xmath78 .
an example of such @xmath51 is @xmath157 . when @xmath157 , eq .
( [ aging_regime ] ) holds ; @xmath158}\rightarrow   \frac{t_w}{g(\triangle t)}\frac{1}{(s_1-s_2)^{\alpha}}\ ] ] as @xmath141 . from eqs .
( [ 2 ] ) and ( [ aging_regime ] ) , @xmath159 } } } c(s_1\triangle t+t_w , s_2\triangle t+t_w )   \nonumber \\   & = & \hat{c}[\mu'g_r(s_1-s_2 ) ] .
\label{st}\end{aligned}\ ] ] it implies that the aging regime is the region where @xmath160 is finite . in this aging regime the slow - process limit of the correlation function in the integrand of eq .
( [ qs_w_a ] ) is a function of @xmath129 and @xmath130 if @xmath153 in eq .
( [ 2 ] ) is not a constant .
it means that the work @xmath9 depends on @xmath74 .
4 .   when @xmath154 and @xmath161 , we prove below that the aging regime for experimental time scales does _ not _ exist .
we assume that the aging regime exists .
it means that a function @xmath162 exists such that @xmath163 $ ] is finite . then , from the condition @xmath164 , @xmath165 } \nonumber
\\ & = & \lim_{{\triangle t \rightarrow \infty }   \atop{\ ; \mu=\frac{t_w}{p(\triangle t ) } } } \frac{t_w}{g[(s_1-s_2)\triangle t ] } \nonumber \\ & = & \lim_{\triangle t \rightarrow \infty } \frac{\mu p(\triangle t)}{g[(s_1-s_2)\triangle t]}.\end{aligned}\ ] ] assuming that @xmath166 , we see from the condition @xmath167 that @xmath165}/ \lim_{{\triangle t \rightarrow \infty }   \atop{\ ; \mu=\frac{t_w}{p(\triangle t ) } } } \frac{s'_2\triangle t+t_w}{g[(s'_1-s'_2)\triangle t ] }   \nonumber \\   & = &   \lim_{\triangle t \rightarrow \infty } \frac{\mu p(\triangle t)}{g[(s_1-s_2)\triangle t]}\frac{g[(s'_1-s'_2 ) \triangle t]}{\mu p(\triangle t ) } \nonumber \\    & = &    \lim_{\triangle t \rightarrow \infty}\frac{g[(s'_1-s'_2)\triangle t ] } { g[(s_1-s_2)\triangle t ] } \nonumber \\    & = & \infty .
\label{mujun}\end{aligned}\ ] ] since eq .
( [ mujun ] ) contradicts the assumption that @xmath163 $ ] is finite , the aging regime for experimental time scales does _ not _ exist in this case . hence , if only this aging regime characterized by such @xmath44 exists , @xmath168 is equal to @xmath105 or @xmath49 .
it implies that the work is independent of @xmath74 .
an example of such @xmath44 is @xmath169 .
consequently , we conclude that _ if there are aging regimes characterized by @xmath51 of case 2 . or case 3 . and the function @xmath46 is not a constant and the fdt violation factor @xmath41 is not equal to @xmath105 , the work for the slow process depends on the path of changing the external field @xmath74 . in this case , contrary to the fundamental assumption of thermodynamics , the work for the slow process is not the difference of a state function . _
it has been explicitly checked on several disordered models that in the long - time limit two situations with different dynamical behavior seem to exist @xcite :    * there are systems with only one correlation scale apart from the quasiequilibrium regime , which we call `` single - correlation - scale systems '' . for these systems ,
the correlation function scales as @xmath170 $ ] and the fdt violation factor @xmath41 is a constant .
equilibrium states of these systems are solved by a one step replica symmetry breaking ansatz .
examples are the @xmath0-spin spherical spin - glasses @xcite and a lennard - jones binary mixture which is a model of structural glasses in a glassy state @xcite .
the real spin - glasses such as agmn seem to belong to this class since it is known that the correlation function is scaled with single scaling function @xmath8 .
* there are systems such as sherrington - kirkpatrick ( s - k ) model which have an infinite number of correlation - scales apart from the quasiequilibrium regime @xcite , which we call `` multi - correlation - scale systems '' .
for these systems , ultrametricity in time holds for any correlation such that @xmath171 ;
@xmath172 \ \mbox{when } \ t_1>t_2>t_3 \
\mbox{and } \         t_3\rightarrow \infty$ ] .
the fdt violation factor @xmath41 is a nontrivial function of @xmath173 .
equilibrium properties of these systems are solved by a full replica - symmetry breaking ansatz .    in this section ,
we analyze the path - dependence of the work for the slow process in the aging regime for these two classes of systems by using the general results obtained in the previous section , sec .
[ section4 ] .
it is important to note that only these two classes of systems seem to exist @xcite .      for the single - correlation - scale systems , scaling of the correlation function
@xmath170 $ ] holds .
hence , the work for the slow process in the aging regime is given as @xmath174\}. \label{w_qs_a_2}\end{aligned}\ ] ] three explicit choices of the scaling function @xmath8 have been proposed so far ; @xmath175 $ ] and @xmath176 $ ] @xcite .
we analyze below the properties of the work for the slow process for each three cases .    1 .
when @xmath56 , the correlation function decays in power - law as found in the trap model @xcite and the @xmath0-spin spherical spin - glasses @xcite .
the slow - process limit of the correlation function in eq .
( [ w_qs_a_2 ] ) depends on @xmath129 and @xmath130 as @xmath177 this is an example of eq .
( [ t ] ) in previous section .
2 .   when @xmath178         ( \alpha < 1)$ ] @xcite , the slow - process limit of the correlation function in eq .
( [ w_qs_a_2 ] ) depends on @xmath129 and @xmath130 through the scaling function which behaves as @xmath179.\ ] ] this is an example of eq .
( [ st ] ) .
3 .   when @xmath180 ( a > 1)$ ] @xcite , the slow - process limit of the correlation function which appears in eq .
( [ w_qs_a_2 ] ) depends on @xmath129 and @xmath130 through the scaling function which behaves as @xmath181.\ ] ] this is another example of eq .
( [ st ] )    to my knowledge , only three types of the scaling function @xmath8 discussed above are known so far for single - correlation - scale systems . for each of three cases the slow - process limit of the correlation function in the integrand of eq .
( [ qs_w_a ] ) is a function of @xmath129 and @xmath130 .
it means for the single - correlation - scale systems that the work for the slow process in the aging regime depends on the path of changing the external field @xmath74 and contrary to the assumption of thermodynamics the work for the slow process is not a difference of a state function .    in order to see more explicit form of the work for the slow process in the aging regime
, we introduce the modified fdt eq .
( [ generalized_fdt ] ) with constant violation factor @xmath41 which has been verified to be valid for the single - correlation - scale systems @xcite .    from eq .
( [ generalized_fdt ] ) , the susceptibility is given by @xmath182 where @xmath183 .
hence , from eq .
( [ qs_w_a ] ) , the work for the slow process is given by @xmath184\ } \label{qs_w_a_2}.\end{aligned}\ ] ]    since @xmath185 when @xmath186 , the second term on the right hand side is non - zero .
hence , if @xmath187 is finite the work for the slow process is not the difference of a state function since the work depends on the path of changing the field @xmath74 . in addition , the value of the work is different from the value of the work for the slow process in the quasiequilibrium regime : the first term on the right hand side of eq .
( [ qs_w_a_2 ] ) .
in addition , it is also shown that since the correlation function is positive and smaller than the dynamical e - a order parameter in this regime , there are bounds for the value of the quasistatic work as @xmath188 eq .
( [ bound ] ) shows that the work @xmath9 in the aging regime is independent of the path of changing the field and coincides with the work @xmath9 in the quasiequilibrium regime when the effective temperature @xmath187 is infinite , which holds for the @xmath189-spin spherical spin - glass model @xcite and a model of phase separation @xcite .      in the previous subsection
, we saw that the work for the slow process depends on the path of changing the external field when @xmath75 .
we show here that the path - dependence of the work for the slow process is strong enough to be observed in experiments by evaluating the path - dependence of the @xmath190-spin spherical spin - glass model .
first , we describe the several pieces of information needed to evaluate eq .
( [ qs_w_a_2 ] ) .
it is known that the scaling of the correlation function , such that @xmath191 and the modified fdt with the effective temperature , i.e. eq .  ( [ generalized_fdt ] ) , hold for the model @xcite . when @xmath192 @xcite , @xmath193 has been obtained .
the dynamical e - a parameter @xmath49 and the effective temperature are given by @xmath194 and @xmath195 when @xmath192 .
in computation of these values , we uses the formulae : @xmath196    in order to evaluate the magnitude of the path - dependence , we compute the relative difference between the works for two paths ( @xmath121 and @xmath197 ) which is defined as the difference between the works for the two paths divided by the average of the two works .
the results obtained from eq .
( [ qs_w_a_2 ] ) are as follows ; when the ratio of the measurement time to the waiting time is @xmath198 , the relative difference is equal to @xmath199 , i.e. about @xmath200 ; when the ratio is @xmath201 , the relative difference is equal to @xmath202 , i.e. about @xmath203 .
hence , we can say that the path - dependence of the work in the slow process can be observed for the @xmath190-spin spherical spin - glass model .
the relative difference is plotted against the ratio of the waiting time to the measurement time in fig .
[ figure1 ] .
it tells that contrary to intuition the path - dependence becomes larger , as one increases the measurement time with keeping the waiting time fixed .      for the multi - correlation - scale systems such as s - k model
, it is known that the ultrametricity in time holds for all correlations such that @xmath171 in the long - time limit @xcite ; @xmath204 \ \ \mbox{when } \
t_1 > t_2 > t_3 .
\label{ultrametric_relation_2}\ ] ]    by using this ultrametric relation eq .
( [ ultrametric_relation_2 ] ) , we prove below that the correlation function in the aging regime is a constant in the two cases ( case 2 and case 3 in sec .
[ section4 ] ) where the work for the slow process can depend on the path of changing the external field @xmath74 .
thus , for these two cases , the work for the slow process is independent of the path @xmath74 .
it has been proven in sec .
[ section4 ] for the other two cases that the work for the slow process can not depend on @xmath74 .
thus , the work for the slow process for the multi - correlation - scale systems is independent of the path of changing the field in all the experimental time regimes .    1 .   when @xmath147 , by substitution @xmath205 in eq .
( [ ultrametric_relation_2 ] ) and assuming that @xmath206 , we see that @xmath207 .
\label{ult2}\ ] ] hence , @xmath208 where @xmath209 .
thus , from eq .
( [ ult2 ] ) , @xmath210.\ ] ] since it implies that @xmath211 @xmath212 is a constant and hence @xmath213 $ ] is a constant .
2 .   when @xmath154 and @xmath155 exists , the ultrametric relation eq .
( [ ultrametric_relation_2 ] ) is written as @xmath214 \label{ult3}\end{aligned}\ ] ] where @xmath97 is a constant such that @xmath215 .
the second argument of @xmath216 in eq .
( [ ult3 ] ) is rewritten as @xmath217}}c(\tau+t',t ' ) \nonumber \\ & = & \lim_{{\tau\rightarrow\infty}\atop{\ ; \mu g_r(a)=t'/g(\tau)}}c(\tau+t',t ' ) \nonumber \\ & \equiv & \bar{c}[\mu g_r(a ) ] , \label{cat}\end{aligned}\ ] ] where @xmath218 .
similarly , the first argument of @xmath216 in eq .
( [ ult3 ] ) is rewritten as @xmath219 \nonumber \\ & = & \bar{c}[\mu g_r(1-a)]\end{aligned}\ ] ] where the last equality is shown from eq .
( [ cat ] ) .
consequently , @xmath220,\bar{c}[\mu g_r(a)]\}.\ ] ] since @xmath221 , @xmath222 is a constant .
we proved above only that the correlation function is a constant within a correlation - scale since ultrametricity in time is valid within a correlation - scale .
on the other hand , even for the multi - correlation - scale systems , the value of the correlation function becomes @xmath105 when the time difference is much larger than the waiting time .
so , value of the correlation function must depend on an infinite number of correlation - scales .
hence , there is possibility that there is an infinite number of experimental time regimes characterized by different values of the work for the slow process .
summarizing our results , for systems with aging we have considered the experimental time scale dependence of the work exerted by modulating the external field in an isothermal process .
we have shown that there are two experimental time regimes characterized by different behavior of the work . in the quasiequilibrium regime @xmath223 ,
the work for the slow process is independent of the path of changing the external field and the value of the work for the slow process is different from the change of the helmholtz free energy obtained by statistical mechanics with the assumption of ergodicity . when the process is not infinitely slow , the difference between the work for non - slow process and the work for the slow process is positive for any path of changing the filed , i.e. the principle of minimum work holds .
in addition , we derived the dependence of the difference on the measurement time from the long - time behaviors of the correlation function .    in the aging regime @xmath224 ,
whose precise meaning is given in sec .
[ section4 ] , contrary to a fundamental assumption of thermodynamics , the work exerted in the slow process on the single - correlation - scale systems for which the fdt violation factor @xmath41 is not equal to @xmath105 depends on the path of changing the external field .
it was examined that for the @xmath190-spin spherical spin - glass model the magnitude of the path - dependence is strong enough to be observed in experiments .
the conditions when the work for the slow process depends on the path is also derived generally . on the other hand , for the multi - correlation - scale systems , the work for the slow process is independent of the path of changing the field .
however , the work for the slow process still may depend on history of a system since there is possibility that there is an infinite number of experimental time regimes characterized by different values of the work for the slow process . in order to attack on this problem , information on time scale of the aging regime of multi - correlation - scale systems beyond time
reparametrization invariance is necessary @xcite .
it is important to note that these results have a much broader range of validity , because they were derived without any assumption on microscopic properties of a system .
we consider the implications of our results on thermodynamics for systems with aging . in the quasiequilibrium regime , our results are the same as the consequences of thermodynamics , where the work for the slow process is independent of the path of changing the external field and the work for a non - slow process is larger than that for the slow process .
only discrepancy is difference between the value of the work for the slow process and that of the free energy obtained by statistical mechanics with the assumption of ergodicity .
however , this discrepancy does not mean that the framework of usual thermodynamics is invalid since the free energy can be defined by the work for the slow process .
this suggests the usual framework of thermodynamics is valid in the quasiequilibrium regime even though the system is out of equilibrium .
we considered only the properties of the work in the linear response regime .
thus , in order to prove validity of usual thermodynamics in the quasiequilibrium regime , the relation between the heat and the entropy has to be considered and extension of our analysis outside the linear response regime has to be performed . for these purposes
, it would be interesting to analyze langevin dynamics of spin - glass models with stochastic energetics @xcite .
since the work for the slow process in the quasiequilibrium regime is independent of history of a system in linear response regime , it gives a fundamental state function for thermodynamics in the quasiequilibrium regime and is worth being measured in experiments .
we discussed an experiment to measure it .
our results in sec .
[ section3 ] for the measurement time dependence of the work shows how long measurement time is needed to reach the expected accuracy .
however , since the waiting time can be long but should be finite in reality , the measurement time has to be finite in order to keep the experimental time scale in the quasiequilibrium regime .
otherwise the value of the work does not come close to that of the work for the slow process in the quasiequilibrium regime since the experimental time scale is in the aging regime .
hence , one can not expect that the value of the measured work coincides with that of the work for the slow process with arbitrary accuracy .
in other words , the accuracy of measurement of the state function is bounded by the length of the waiting time .    on the other hand , in the aging regime ,
our results imply that in contrast to the usual framework of thermodynamics it is essential to consider the experimental time scales such as the waiting time and the measurement time . furthermore , history of a system as a state variable is necessary to describe thermodynamic properties of glassy materials .
although constructing thermodynamics for the aging regime is a tough work , it is a fruitful and challenging problem to be investigated since some universal properties in the aging regime are known @xcite .
recently , a framework of thermodynamics for the aging regime when the temperature changes in time is proposed @xcite .
the framework avoids treating the problem of path - dependence by fixing a path to change the temperature .
it is not clear whether thermodynamics for a single path to change the temperature is useful .
its validity and usefulness should be checked .
i wish to acknowledge valuable discussions with takashi odagaki , akira yoshimori , hiizu nakanishi , kazuo kitahara , osamu narikiyo , miki matsuo , ken sekimoto , shin - ichi sasa , yoshihisa miyamoto , hajime takayama , koji hukushima and hitoshi inoue .
in this appendix , we derive the long measurement time behavior of the deviation from the work for the slow process for four types of relaxation of the correlation function described in sec .
[ difference ] .          at first
, we consider the deviation when @xmath225 $ ] is equal to @xmath105 . here , in order to prove that the deviation is proportional to the inverse of the measurement time @xmath13 , we show that the proportional constant is finite .
the proportional constant is given as @xmath230.\end{aligned}\ ] ] since @xmath231 exists in the case , @xmath227 is finite .
hence , the order of the limit @xmath78 and the integration can be changed .
thus , the proportional constant is given by @xmath232 where @xmath233 is defined by eq .
( [ k ] ) .
consequently , it is proved that the proportional constant is finite and the deviation is proportional to @xmath234 .
next , we consider the deviation when @xmath225 $ ] is equal to a finite value @xmath226 . in order to prove that the deviation is proportional to @xmath235 ,
we show that the proportional constant is finite . dividing the region of integration of eq .
( [ app : dev1 ] ) into five regions , the proportional constant is given by @xmath236 \label{app : dev2}.\end{aligned}\ ] ]    we discuss the behavior of the function @xmath237 in each region .
when @xmath238 , from the definition of @xmath227 , eq .
( [ k ] ) , it is shown that @xmath239.\ ] ] thus , @xmath237 is finite even when @xmath78 . on the other hand ,
when @xmath240 , we express the function @xmath237 by using a function defined as @xmath241 : @xmath242   \nonumber \\   & = & \int_0 ^ 1 d\tau [ c(\tau)-q]+\int_1^{s \triangle t } d\tau   [ \frac{c}{\tau}+f(\tau ) ]   \nonumber \\   & = & \int_0 ^ 1 d\tau [ c(\tau)-q]+c \log s+c
\log \triangle t+ \int_1^{s \triangle t } d\tau f(\tau ) .
\end{aligned}\ ] ] since @xmath243 , the fourth term on the right hand side is finite even when @xmath13 is infinite .
thus , the behavior of the function is given as @xmath244 when @xmath240 . by evaluating the each term of eq .
( [ app : dev2 ] ) with using the above behavior of the function @xmath237 , it is straightforward to show that @xmath245 the case where @xmath225 $ ] does not exist ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~      in order to find out the leading term in the behavior for the long measurement time , we divide the region of integration of eq .
( [ devi_work ] ) into three regions as @xmath248 where we assume that the parameter @xmath97 satisfies @xmath215 . since all the integrands ( the time derivative of the field and the correlation function ) are finite
, it can be shown that the first two terms on the right hand side are bounded as @xmath249 and @xmath250 .\ ] ] in addition , we notice that since @xmath251 in the region of integration of the third term , @xmath252 is large when the measurement time @xmath13 is large .
thus , the third term is determined by the long - time behavior of the correlation function .    at first
, we treat the case where the correlation function behaves at a long time as @xmath253 .
we divide the region of integration of the third term of eq .
( [ division ] ) into three regions as @xmath254 since all the integrands are finite , it is shown that the last two terms are bounded by @xmath255 and @xmath256    on the other hand , by using the long - time behavior of the correlation function , the first term is given as @xmath257 thus , we find that this term is @xmath258 $ ] . assuming the free parameter @xmath97 is larger than the exponent @xmath259
, one sees that the leading term is eq .
( [ power1 ] ) by comparing the orders of all the terms in eqs .
( [ division ] ) and ( [ division2 ] ) .
hence , the long measurement time behavior of the deviation is given by eq .
( [ power ] ) .
next , we consider the case where the long - time behavior of the correlation function is given by @xmath260 . again , we start by evaluating the third term on the right hand side of eq .
( [ division ] ) . by using the long - time behavior of the correlation function ,
the term is given by @xmath261 since in the region of integration the following inequality holds : @xmath262 by neglecting numbers of @xmath263 the third term of the right hand side of eq .
( [ division ] ) is reduced to @xmath264 furthermore , by neglecting numbers of @xmath265 $ ] , the term is reduced to @xmath266    consequently , one sees that the leading term is the above one by comparing the orders of all the terms in eqs .
( [ division ] ) and ( [ division2 ] ) .
thus , the long measurement time behavior of the deviation is given by eq .
( [ log ] ) .    finally , in order to determine the order of the neglected term we determine the value of the free parameter @xmath97 , since the neglected term contain the free parameter @xmath97 .
since in terms of @xmath97 the order of the neglected terms are given as @xmath263 and @xmath267 $ ] , we can determine the value of @xmath97 so that the value of the neglected terms are minimized .
thus , assuming that the neglected term is given by @xmath268 where @xmath259 and @xmath269 are constants , the value of @xmath97 is determined by @xmath270   = 0.\ ] ] since the solution is @xmath271 , one sees that the order of neglected term is @xmath272 $ ] as shown in eq .
( [ log ] ) .
the present address is research institute for applied mechanics , kyushu university , kasuga 816 - 8580 , japan .
electronic address : mituhiro@riam.kyushu-u.ac.jp for recent reviews , see e.g. e. vincent , j. hammann , m. ocio , j .-
bouchaud , and l.f .
cugliandolo , in _ proceedings of sitges conference on glassy systems _ , edited by e. rubi ( springer , berlin , 1996 ) ; j .-
bouchaud , l. cugliandolo , j. kurchan , and m. mzard , in _ spin glasses and random fields _ , edited by a.p .
young ( world scientific , singapore , 1997 ) .
these articles are found also at e - print cond - mat/9607224 and cond - mat/9702070 .
l. lundgren , p. svedlindh , p. nordblad , and o. beckman , phys .
* 51 * , 911 ( 1983 ) ; m. alba , j. hammann , m. ocio , and ph .
refregier , j. appl . phys . * 61 * , 3983 ( 1987 ) . l.f .
cugliandolo and j. kurchan , phys .
* 71 * , 173 ( 1993 ) .
cugliandolo and j. kurchan , j. phys .
a * 27 * , 5749 ( 1994 ) .
cugliandolo , j. kurchan , and l. peliti , phys .
e * 55 * , 3898 ( 1997 ) .
s. franz and h. rieger , j. stat .
phys . * 79 * , 749 ( 1995 ) ; e. marinari , g. parisi , f. ricci - tersenghi , and j.j .
ruiz - lorenzo , j. phys .
a * 31 * , 2611 ( 1998 ) .
cugliandolo , d.r .
grempel , j. kurchan , and e. vincent , europhys . lett . * 48 * , 699 ( 1999 ) .
g. parisi , phys .
* 79 * , 3660 ( 1997 ) .
barrat and w. kob , europhys . lett .
* 46 * , 637 ( 1999 ) .
grigera and n.e .
israeloff , phys .
lett . * 83 * , 5038 ( 1999 ) .
l. berthier , j .- l
. barrat , and j. kurchan , eur .
j. b * 11 * , 635 ( 1999 ) .
the reason why the usual fdt holds in the quasiequilibrium regime though a system is out - of - equilibrium is that the change of the free energy stops at the long time limit , since the fdt violation is bounded by the time - derivative of the free energy @xcite .
cugliandolo , d.s .
dean , and j. kurchan , phys .
lett . * 79 * , 2168 ( 1997 ) . when @xmath273 , @xmath274 does not hold even if @xmath275 since the time difference @xmath252 is zero .
however , the case can be ignored , since the susceptibility is finite there and the point @xmath273 does not contribute to the value of the integral by itself .
ogielski , phys .
b * 32 * , 7384 ( 1985 ) ; a. ito , h. aruga , e. torikai , m. kikuchi , y. syono , and h. takei , phys . rev
. lett . * 57 * , 483 ( 1986 ) .
m. ocio , h. bouchiat , and p. monod , j. mag . mag .
* 54 * , 11 ( 1986 ) .
ngai , comments solid state phys . *
9 * , 127 ( 1979 ) ; * 9 * , 141 ( 1980 ) ; a.t .
ogielski , phys .
b * 32 * , 7384 ( 1985 ) ; m.h . cohen and g.s .
grest , phys .
b * 20 * , 1077 ( 1979 ) .
it has been shown from the massless langevin equation without the assumption of a linear response that @xmath276 @xcite . since the correlation function derived from the massless langevin equation
obeys the exponential relaxation , our result agrees with that of ref .
k. sekimoto and s. sasa , j. phys .
jpn . * 66 * , 3326 ( 1997 ) .
bouchaud , j. phys .
i ( france ) * 2 * , 1705 ( 1992 ) .
this form of the scaling function was proposed to account for the results of experiments in polymer glasses @xcite , then used in analyses of aging effects in decay of the thermoremanent magnetization @xcite , and found in the exact solution of the asymmetric spherical sk model @xcite .
struik , _ physical aging in amorphous polymers and other materials _ ( houston , 1978 ) .
m. ocio , m. alba , and j. hammann , j. physique lett .
( france ) * 46 * , l-1101 ( 1985 ) ; m. alba , m. ocio , and j. hammann , europhys . lett . * 2 * , 45 ( 1986 ) . s. franz and m. mzard , europhys . lett . * 26 * , 209 ( 1994 ) : physica a * 209 * , 1 ( 1994 ) ; l.f .
cugliandolo , j. kurchan , p.le doussal , and l. peliti , phys .
lett . * 78 * , 350 ( 1997 ) .
this form of the scaling function is suggested by the numerical data from the analysis of a point particle in a random potential with infinite dimension @xcite and is used to scale data for the thermoremanent magnetization and the out - of - phase susceptibility .
cugliandolo and p.le doussal , phys .
e * 53 * , 1525 ( 1996 ) ; l.f .
cugliandolo , j. kurchan , and p.le doussal , phys .
lett . * 76 * , 2390 ( 1996 ) .
cugliandolo and d.s .
dean , j. phys .
a * 28 * , 4213 ( 1995 ) .
the energy is scaled with a constant @xmath277 which determines the variance of the coupling constant between spins as @xmath278 the equations of motion of @xmath2 and @xmath4 of the multi - correlation - scale systems in the long waiting time limit have an infinite set of invariance called `` time reparametrization invariance '' @xcite ; if we perform an arbitrary reparametrization of time @xmath279 and we redefine the two - time quantities as @xmath280 \ \
\hat{g}(\hat{t},\hat{t}')=g[\xi(t),\xi(t')].\ ] ] then , the redefined functions @xmath281 and @xmath282 satisfy the same equations of motion . it is important to note that this invariance is a consequence of having neglected the time derivatives in taking the long waiting time limit .
the full equations of motion have no such invariance and the solution is unique .
k. sekimoto , j. phys .
66 * , 1234 ( 1997 ) .
m. nieuwenhuizen , phys .
e * 61 * , 267 ( 2000 ) .


due to the observations of cosmic microwave background ( cmb ) radiation  @xcite , large scale structure ( lss )  @xcite , and ia supernovae  @xcite we known that the universe is actually in accelerated expansion . a direct way to incorporate this accelerated expansion of universe into einstein
general relativity theory  @xcite is to include the cosmological constant @xmath6 .
the inclusion of this constant has the same effect of a non - zero vacuum energy density , @xmath7 , in such way that the pressure @xmath8 and the density @xmath9 are related by a state equation with the form @xmath10 .    on the other hand ,
the accelerated expansion of universe can also be described by adding to the universe content a homogeneous fluid with energy density @xmath11 , the so called _ dark energy _
this fluid has positive energy density @xmath12 , and negative pressure @xmath13 , in such way that @xmath14 .
this pressure is then responsible for the accelerated universe expansion .
it is a remarkable fact that @xmath15 of universe content must be in the dark energy form .
also , dark energy would be uniformly distributed in all space , and so , its density is constant in all points and times .
this is in contrast with the time evolution of baryonic matter density @xmath16 , that diminishes due to the expansion of the universe .
it is called _ the cosmic coincidence _ to the fact that today , the baryonic matter density , @xmath16 has approximately the same value of dark energy density @xmath12  @xcite . to compare different cosmological models , it is common to define the density parameter @xmath17 that it is the ratio between the density @xmath18 with the critical density @xmath19 of the friedmann universe , @xmath20 . in this sense
the cosmic coincidence implies in the equality @xmath21 .
this coincidence can be viewed as an indicative of the existence of some dynamical effect that relates both scales of baryonic matter and of the dark energy density .
nevertheless , among all the models in the literature that are devoted to explain the accelerated expansion of universe , there is a class of dynamical models that obtain the desired negative pressure due to the inclusion of a scalar field that is the responsible for the variations in the expected value of vacuum energy . as the neutrino squared mass difference , @xmath22 ev@xmath23 , where @xmath24 and @xmath25 are the masses of third and first neutrino mass state , is of same order of dark energy scale , it is straightforward to think a model in which the scalar field couple to neutrinos and hence , the total energy of the fluid can vary slowly as the neutrino density decreases  @xcite . in the _ mass varying neutrino _ models
the inclusion of a scalar field allow the coupling of neutrino and dark energy densities due to the non - standard neutrino - scalar coupling . the main consequence of this coupling to the neutrino physics is the fact that now the neutrino masses depends on the medium density .
this field could couples neutrinos to the baryonic matter and also to the neutrino background  @xcite .
the consequences in cosmology of mavan s had been studied in last years  @xcite and references therein .    for the neutrino phenomenology the consequences of mavan s were studied for the solar and atmospheric neutrinos  @xcite .
we will use the data from super - kamiokande ( sk for now ) experiment  @xcite . in the experiment
sk , down - going neutrinos , the ones that are produce in the atmosphere immediately over super - kamiokande detector , reaches sk with the the cosine of zenith angle @xmath26 and travels approximately 20 km in the atmosphere . on the other hand ,
up - going neutrinos are produced in the opposite side of the planet and crosses all the earth before reaches sk in the the cosine of zenith angle @xmath27 direction .
the so - called _ up - down _ asymmetry of neutrino events in sk is recognized as the first experimental corroboration of neutrino flavor oscillations .
we explicit here that atmospheric neutrinos that arrives at sk from different directions travels different distances and crosses regions of very distinct densities ( see next section for details ) .
this fact makes the angular distribution of events in sk a good place to looking for dependence of medium density effects in the propagation of neutrinos . in fact , the sk experiment reported no improvement of data fit due to the inclusion of diagonal mavan s mechanism in the propagation of atmospheric neutrinos  @xcite . for our knowledge
there is no analysis made for non - diagonal mavan s mechanism for atmospheric neutrinos .
we also apply the mavan s model to describe the beam muon neutrinos and anti - neutrinos at minos experiment  @xcite , where a muon neutrino or muon anti - neutrino beam with energy between few hundred of mev and few gev .
the neutrino or anti - neutrino beam travels a few hundred of kilometers inside earth s crust before reach the far detector . by the comparison between the number of neutrino events in near and far detector minos collaboration
has measure precisely the standard neutrino parameters @xmath28 and @xmath29 .
see for example  @xcite and references therein .
in mavan s context , the main difference between minos beam and sk is that in the former neutrinos cross only the upper crust of earth , that can be described by a constant matter density and the latter cross a different non - constant densities . in this way
the mavan s effective oscillation for minos , with a constant density and with sk with a variable density allow us to test the essential characteristic of mavan mechanism , the density dependence of the neutrino mass differences .
the paper is organized as follows . in the section  ( [ sec : model ] ) we describe the framework to mavans that we adopt and in sec .
( [ oscilograms ] ) we show the changes in neutrino oscillations due to mavans .
then in sec .
( [ sec : muons ] ) we compare our results for the angular distribution of events in sk for neutrino oscillations without and with mavans .
also in section  ( [ sec : minos ] ) we show how mavan s mechanism affects the allowed range of oscillations parameters for minos experiment and in section  ( [ sec : chi2 ] ) we perform the @xmath30 analysis of sk atmospheric neutrino combined with minos beam data , and the constrains to mavan s model that we obtain .
conclusions are in section   ( [ sec : conclusion ] ) .
in the neutrino mass - mixing formalism , the time evolution of neutrino flavor eigenstates is given in terms of neutrino mixing , in which one has to describe the flavor eigenstate and hence , the time evolution of atmospheric neutrinos , in the two neutrino flavor approximation , is given by the evolution equation @xmath31\left(\begin{array}{c } \nu_{\mu}\\ \nu_{\tau } \end{array}\right)~ , \label{sch}\end{aligned}\ ] ] and the hamiltonian assumes the form , in the mass basis , @xmath32~ , \label{mmatrix}\end{aligned}\ ] ] here @xmath33 is the neutrino energy , @xmath34 is the square difference of mass eigenstates , and @xmath35 is the mixing matrix @xmath36 where @xmath37 .
the muon neutrino survival probability is @xmath38 where @xmath39 is the squared mass difference , l is the distance and @xmath33 is the neutrino energy .
the formalism that we adopt to include mavans in the neutrino propagation has the standard model of particles plus a light scalar field @xmath40 with mass @xmath41 that couples with neutrinos @xmath42 , i=1,2,3 and with fermion fields f = e , n , p . in such model , the modification due to mavan
s is the introduction of a fermion density dependent term in each one of matrix elements in the neutrino mass matrix , eq .
( [ mmatrix ] ) , and so , including mavan s , the flavor neutrino evolution is described by the generalization of eq .
( [ sch ] ) .
explicitly we have @xmath43= h^{{\rm flavor}}_{{\rm mavan } } \left[\begin{array}{c } \nu_{\mu}\\ \nu_{\tau } \end{array}\right ] ,   \label{eqevo - drg}\end{aligned}\ ] ] where the modified hamiltonian in the mavan s framework is @xmath44 where the mixing matrix u is defined in eq .
( [ u ] ) and the mavan hamiltonian has the form @xmath45~. \label{mevo - drg}\end{aligned}\ ] ] here @xmath46 and @xmath47 are the mavan s parameters and @xmath48 is the function of the earth matter density @xmath18 that neutrinos feels while cross the earth . when @xmath49 we recover the standard neutrino evolution given in eq .
( [ sch ] ) .
we can classify the behavior of mavan mechanism given in eq .
( [ mevo - drg ] ) in two types : ( a ) when @xmath50 and @xmath51 and ( b ) when @xmath52 and @xmath53 . in the former case
the neutrino survival probability , for constant density , is given by @xmath54 where the effective mass scale is given by @xmath55 $ ] .
in this case ,
when @xmath51 the amplitude of oscillations , @xmath56 is the same as in the standard neutrino oscillations in eq .
( [ standardoscillation ] ) , and the phase of the oscillations , proportional to @xmath57 , have now a matter density dependence .
this was the case mostly worked in the literature  @xcite .
the latter case , the hamiltonian in the mass basis is non - diagonal and for our knowledge it was _ not explored _ in the literature for atmospheric neutrino phenomenology .
the probability , for constant density @xmath18 , can be written as @xmath58 where the amplitude , @xmath59 and the phase of oscillations , @xmath60 are different from the usual standard oscillation scenario .
the mavan mass difference is given by @xmath61 where the mavan mass difference depends on the medium density and the explicit expression for the amplitude is @xmath62 ^ 2 \label{eq : sineff}\ ] ] that also have a dependence on the medium density .
the angle @xmath63 is defined as @xmath64 ^ 2 } { 1   + \left [ 2\alpha^{2}_{32}g(\rho)\right]^{2 } } \label{eq : neta}\ ] ] in mavan s case for @xmath65 induce that the mass basis it is not diagonal and the parameter @xmath63 is the mixing angle that diagonalizes the mass basis as shown in reference  @xcite .
this made that the non - diagonal mavan s would induces neutrino oscillations even if the vacuum mixing angle was zero , @xmath66 .    in the case of standard oscillations , see eq .
( [ standardoscillation ] ) we have the symmetry @xmath67 , but in mavan case with @xmath65 we broken this symmetry and the eq .
( [ eq : sineff ] ) have different results for @xmath68 compared when @xmath69 . for vanishing mavan parameters , @xmath70 , we have @xmath71 , and we recover the standard neutrino oscillation . for very large values of mavan parameters we have @xmath72 ,
that implies that @xmath73 .
if we have large mixing angles @xmath74 we will have suppression of the amplitude , but for smaller mixing angles will have a enhancement of the amplitude .
most of the previous analysis of mavan works for the first two generations  @xcite and then theirs constrains are related to parameters of the first generation @xmath75 and @xmath76 similar to the parameters defined in eq .
( [ mevo - drg ] ) . from the reference  @xcite that provide a upper bound for the elements of matrix @xmath77  ev at 90 % c.l .
the only case that works in the mavan scenario for the second and third families are the references  @xcite . in this
works they assume the non - zero parameters @xmath78 and the medium matter dependence is given by @xmath79 , where @xmath18 is the matter density , @xmath80 and _ n _ is a free parameter .
it was ruled out diagonal mavan to be the dominant oscillation scenario for atmospheric neutrinos and they constraint the _
n _ parameter to be in the range @xmath81 at 90 % c.l .
@xcite .
now we are going to compare the oscillation probabilities without mavan s and with mavan s to understand the changes in the oscillation probabilities .
most of previous analysis use a model for the mavan s mass @xmath82 given as a function of parameter with dimensions of energy such as @xmath83 , with the index i denoting the diagonal mass eigenstates i=1,2 and the non - diagonal mass eigenstates i=3 ; with different functions @xmath48 as function of the matter density @xmath18 , and a parameter @xmath84 with dimension of energy .
we decide to adopt a dimensionless parameters in this work , the @xmath85 as described in eq .
( [ mevo - drg ] ) , but we can related our parameters @xmath85 with the previous analysis by the replacing of @xmath86 , where @xmath87 is the relevant mass difference of the analysis .    in this work
we decide to study the case when for non - diagonal mavan parameter , @xmath52 and @xmath53 , vide eq .
( [ mevo - drg ] ) . as
an example e will use in this work the mavan density dependence as used in ref .
@xcite , given by @xmath88 in this work we will use the matter profile of earth that we take from  @xcite . where @xmath18 is the matter density that neutrino crosses and
@xmath89 g/@xmath90 is the matter density at earth core .
this choice is motivate to generates a soft and well behaved function even with the abrupt variations of earth density profile , other choices made the computation numerical unstable as reported in ref .
also it have a finite limit for @xmath91 and @xmath92 always . in this work
we will use the matter profile of earth that we take from  @xcite .
we solve numerically compute the muon neutrino survival probability , @xmath93 for the general case of a neutrino crossing the earth from different chords .
we can related the traveled distance by neutrino _
l _ to the zenith angle @xmath94 by @xmath95 , where @xmath96 km is earth radius and _ h _ is the point of atmosphere were neutrinos are produced , approximately @xmath97 km .
for @xmath98 we have the maximum ( mininum ) distance .
we shown in fig .
( [ fig : pmumu - cz ] ) the muon neutrino probability as a function of cosine of zenith angle , @xmath99 for fixed values of the amplitude of the mixing angle @xmath100 , for the mass difference @xmath101 ev@xmath23 and for a fixed energy @xmath102  gev and several values of @xmath47 : @xmath103 we have respectively the mavan probabilities as dashed curves and without mavan as solid curves in fig .
( [ prob - cz-04 ] , ( [ prob - cz-01]),([prob - cz-02 ] ) and ( [ prob - cz-03 ] ) . for values of @xmath104 , the neutrino traveled in vacuum only and the oscillation probability with and without mavan are the same . for @xmath105 ,
the neutrino cross inside the earth and the mavan effect begin to pop up . to see more clearly the effect of mavan parameters we can compare the fig .
( [ prob - cz-04 ] ) that have @xmath106 , with fig .
( [ prob - cz-03 ] ) that have @xmath107 .
when can see when we increase the mavan parameter @xmath47 , for maximal mixing , the mavan amplitude get smaller and the oscillation phase increases giving more fast wiggles that appear in the probability .
the increase of wiggles made the maximums and minimums move to smaller values of @xmath108 .
the concept of oscillograms is interesting tool to understand the complete behavior of neutrino probability in some model for neutrino oscillation .
we plot in fig .
( [ fig : pmumu_alpha0.1 ] ) the oscillograms of muon neutrino survival probability , denoted by @xmath109 , as function of neutrino energy @xmath110 and the cosine of zenith angle , @xmath111 . in the fig  ( [ oscillogram-1-mavan ] ) we show for the standard neutrino case and in the others plots for increasing bigger values of mavan parameter .
for the cosine of zenith angle , @xmath112 , when the neutrino did not cross the earth , we have zero mavan effect and for @xmath113 , the muon survival probability is modified due mavan mechanism . for comparison
we show in the upper - left panel , in the upper - right panel , lower - left panel lower - right panel of n fig .
( [ fig : pmumu_alpha0.1 ] ) for mavan parameters @xmath47 respectively equal to 0.0 , 0.1 , 0.5 and 1.0 .
an enhancement of @xmath114 also implies that the positions of maximums and minimums is dislocated to higher values of neutrino energy . as an example of this
let we look to the first minimum in @xmath109 for @xmath115 . in the former three panels of fig  ( [ oscillogram-1-mavan ] )
we see that the first minimum ocurs for @xmath116  gev . in the lower - left panel , where @xmath107 this minimum in @xmath109 had its intensity
reduced ( due to the increase of @xmath117 but also we see that the minimum was dislocated to @xmath118  gev . at @xmath116
gev now we see the first maximum of oscillation that for s.o .
occur at @xmath119  gev .
the same kind of displacement is found for all the maximums and minimums in the lower - left panel of fig  ( [ oscillogram-1-mavan ] ) when compared with the cases in which @xmath120 .
atmospheric neutrinos , composed muon neutrinos and electron neutrinos , are produced all around the earth atmosphere and traveled to super - kamiokande from all directions .
once in the detector they interact producing muons and electrons and super - kamiokande measure the zenith angle dependence of these muons and electrons .
the rate for these events can be computed as @xmath121\theta[e^{{\rm max}}_{\mu}-e_{\mu}(\cos\theta_{z},e_{\nu},\cos\theta_{\mu } ) ]   ~ , \label{nmu}\end{aligned}\ ] ] where @xmath122 is the number of targets in sk , @xmath123 is the livetime , @xmath33 is the neutrino energy , @xmath124 is the cosine of zenith angle ( @xmath94 ) of the neutrino , @xmath125 is the azimuth angle of the neutrino , @xmath126 is the cosine of zenith angle of the muon , @xmath127 is the muon neutrino flux , @xmath128 is the muon neutrino survival probability , @xmath129 is the differential charged current muon - neutrino cross - section . for the integration boundaries , @xmath130 and @xmath131 the initial and final neutrino energies ; @xmath132 and @xmath133 the bin of zenith angle distribution of super - kamiokande experiment , in this we use 10 equal bins of muon zenith angle between @xmath134 .
we compute this expression for the 2 different types of events in sk experiment : the so - called _ sub - gev _ data set and _ mult - gev _ data set .
they correspond respectively to the intervals of @xmath135  gev , @xmath136  gev ) and @xmath137  gev , @xmath138  gev ) .
we use the kinematics of reaction , defined by @xmath139 and @xmath33 variables to set - up the allowed range of muon energy , given by the function @xmath140 and we constrain to be in the _ sub - gev _ and _ mult - gev _ energy range . notice that the zenith angle of leptons it is a function of zenith angle of neutrino , the scattering angle and the energy of neutrino and this produce a stronger averaging effect on the original neutrino direction .    for the electron - like zenith distribution of events in sk we can write @xmath141\theta[e^{\rm max}_{e}-e_{e}(\cos\theta_{z},e_{\nu},\cos\theta_{e } ) ]   ~ , \label{ne}\end{aligned}\ ] ] which it is very similar to the muon events , with the exception that the electron neutrino did not oscillate , due our initial assumption only have oscillations between muon and tau neutrinos .
the muon and electron neutrino fluxes are taken from  @xcite .
the differential cross section follows  @xcite , whe ere we divide the cross section in three parts , first the quasi - elastic neutrino scattering with finite mass corrections  @xcite ; second the one pion contribution and third the dis contribution  @xcite .
due the super - kamiokande did not discriminate between particles and anti - particles , we sum over neutrino and anti - neutrino types .
the muon and electron rate for atmospheric neutrinos have the uncertainties from the prediction of the atmospheric neutrino fluxes , @xmath127 and @xmath142 , computation that can be @xmath143 and the relative error @xmath144 of 5%  @xcite .
because of this error in the absolute normalization we will use as the physical observable is the zenith distribution of number of events and the absolute value of our prediction will be scaled with the experimental data .
also the smallness of the relative error of muon and electron neutrino fluxes compared the error in the absolute number implies a stronger correlation between the fluxes of muon and electron neutrinos , that we should take into account .
we use this formalism to describe two energy ranges of sk data , the so - called _ sub - gev _ and _ mult - gev _ data set for electrons and muons , eq .  ( [ nmu ] ) and eq .  ( [ ne ] ) .
we compute the muon neutrino survival probability , given in eq .
( [ standardoscillation ] ) , with the oscillation parameters given by @xmath101  ev@xmath23 and the amplitude @xmath145 shown as the black curve in fig .
( [ fig : muonsub ] ) and fig .
( [ fig : muonmulti ] ) respectively for sub - gev and multi - gev sample . for these parameters we have no oscillation for @xmath146 and average out oscillation for @xmath147 .
our results match the theoretical curves for the number of events for no - oscillation and standard oscillation of sk experiment . for our computation with mavan s probability
, we use the numerical solution of eq .
( [ eqevo - drg ] ) using the matter profile given in eq .
( [ dmevo1 ] ) and we compute the rate for electron and muon ( given by eq .
( [ nmu ] ) and eq .
( [ ne ] ) ) , the result is the dashed curve in fig .
( [ fig : muonsub ] ) and fig .
( [ fig : muonmulti ] ) .
the effect of non - zero mavan , using the same oscillation parameters and the mavan parameter @xmath148 it to distort the muon distribution making a small oscillation for cosine of lepton zenith @xmath149 , coming from neutrino zenith angle @xmath150 , due the averaging effect mentioned after the eq .
( [ nmu ] ) ; and suppressing the averaging of neutrino oscillations for @xmath147 , both behavior are disfavored by the sk data and from this we expect to have a constrain in the mavan parameter from this data .    in atmospheric neutrinos
, neutrino can came from different directions and they probe different medium density making the mass difference @xmath151 to be different in each point of the travel . on the other hand , for minos , the neutrino travel only cross the upper crust , and we can assume that density along this short chord is constant , we will assume @xmath152  g/@xmath90 from the prem model .
we can use the results of section  ( [ sec : frame ] ) of mavan mechanism for constant medium density where the muon survival tprobability is given by eq .
( [ case2mavan ] ) with the mavan mass difference given by eq .
( [ eq : effdm ] ) and the amplitude by eq .
( [ eq : sineff ] ) .
the interesting is that the functional form of mavan oscillation probability is exactly like the standard oscillation probability , with the replacement @xmath153 and @xmath154 .
are indicated in the figure .
also we assume @xmath155 and @xmath156 ev@xmath23 .
points refers to minos data from  @xcite . ]
until now we present a phenomenological framework for off - diagonal mavan s mechanism and show it changes the @xmath157 atmospheric neutrino oscillation pattern in the sk experiment . however , the sk detector is not the only one sensitive to the neutrino oscillations in the @xmath157 sector and we can use the minos data  @xcite to constrain the non - diagonal mavan parameter . we show in fig .
( [ minoscte ] ) the comparison between the standard neutrino oscillation and the mavan s probability for the setup of minos experiment ( using the l = l@xmath158 and the density @xmath152  g/@xmath90 ) .
the oscillations parameters are fixed as @xmath159  ev@xmath1 and maximal mixing @xmath160 .
when we increase the mavan parameter we have for maximal mixing that the mavan amplitude get suppressed @xmath161 , as you can see that the minimum , in fig .
( [ minoscte ] ) , is less deeper for higher values of non - diagonal mavan parameter @xmath47 . at same that the oscillation phase , @xmath162 is bigger for higher values of @xmath47 , that it implies that the minimum should be for higher values of energy .
we also show for comparison in fig .
( [ minoscte ] ) the ratio of experimental number of events over the theoretical prediction without oscillation as data points to emphasize that when we increase the value of non - diagonal mavan parameter we became further apart from the experimental data .
so far we have shown examples of mavan effects for maximal vacuum mixing , @xmath160 , e.g. in fig .
( [ fig : pmumu - cz],[fig : pmumu_alpha0.1],[fig : neve],[minoscte ] ) .
a subtle effect can happen for the constant medium density case for larger values of mavan parameter , @xmath47 , in this case we have that the mavan amplitude @xmath163 as commented in section  ( [ sec : model ] ) . for small vacuum mixing angles
@xmath164 this implies that the amplitude is enlarged compared without mavan and for larger vacuum mixing angles @xmath164 the situation is the opposite , and you have suppression of oscillation . to give a example we display in fig .
( [ sineff - ii ] ) , the vacuum mixing angle @xmath165 as a function of mavan parameter @xmath47 that gives a fixed value of mavan amplitude @xmath166 ( in another words the inverse function of eq .
( [ eq : sineff ] ) ) .
we can see that it have two solutions : one , the black curve , with small vacuum mixing angle @xmath165 and another solution given by the dashed red curve with large vacuum mixing angle .
this implies that we can have a small mixing angle and large @xmath167 parameter or with can begin a large mixing angle and small parameter @xmath167 and both give a effective large mavan amplitude . for @xmath168 ,
both curves coincide and we have full suppression of oscillation amplitude .
we will see the consequences of subtle effect later in our analysis in section  ( [ sec : chi2 ] ) .
parameter that can made a fixed valur for the mavan amplitude @xmath166 ]
here we give the details of data analysis .
first we begin with the analysis of minos experiment and later to the super - kamiokande analysis for the standard vacuum oscillation scenario and for the mavan scenario .     c.l . for oscillations parameters @xmath170 without ( s.o ) , dotted blue line , and with mavan s , dashed red line for the minos analysis . ]
the minos experiment made a likelihood analysis of theirs data for _ neutrinos and anti - neutrinos _ for standard neutrino oscillation and produces a table with the values of @xmath171 , where @xmath172 , where @xmath173 is the likelihood value and @xmath174 is the likelihood for the best fit .
this table is publicly available in ref .
we can translated the likelihood language into @xmath175 language using @xmath176 .
for the mavan analysis , we can use the property that the mavan probability have the same functional form of standard oscillation , as discussed in section  ( [ sec : minos ] ) , any function of probability also should have similar behavior . therefore the @xmath177 function given in ref .
@xcite as a function of @xmath178 and @xmath179 , should be equal to mavan @xmath180 as a function of the @xmath181 and @xmath162 parameters .
numerically we have @xmath182 where for minos experiment we can use the expression for mavan amplitude @xmath181 and mass difference @xmath162 respectively eq.([eq : sineff ] ) and eq.([eq : effdm ] ) that give the mavan parameters as functions of the vacuum oscillation parameters @xmath178 and @xmath179 .
we use this procedure to get the allowed region for minos only , for the standard oscillation scenario and for the mavan case .
the results is shown in fig .
( [ minosonly - regi ] ) where the dotted blue curve correspond to standard neutrino oscillations region of @xmath183 of c.l .. as we plot the region of allowed region in the standard scenario as an function of @xmath184 ( and not in function of @xmath185 ) then two degenerated minimums ( denoted by up triangles ) do appear .
when we include mavans and minimize with respect to @xmath47 we get the allowed region is given by the red dashed curve which has the best fit for non - zero @xmath47 , and for vacuum mixing angles smaller then maximal @xmath186 and mass differences @xmath187  ev@xmath1 .
in the mavan s allowed region the mixing parameters can have larger values of @xmath188 and smaller values of @xmath28 that are not allowed in standard oscillation scenario .
smaller values of @xmath189  ev@xmath1 are not allowed in the standard scenario due it implies smaller oscillation effect that is in contradiction with the minos data , but in mavan mechanism , the effective mass difference @xmath162 can be larger then vacuum oscillation mass difference compensating the smaller value of @xmath179 .
also the value of @xmath190 in the mavan solution , if we are working in standard scenario will implies a smaller oscillation amplitude that also it is not compatible with minos data , but also we can allow values @xmath184 far for maximal for non - zero @xmath47 .
now we will work with the analysis of sk data , where we will use the muon and electron data for sub - gev and multi - gev samples .
the sample of atmospheric neutrino data is specially interesting for the mavan oscillation effect because it is composed by events produced by neutrinos traveling in vacuum and in matter from the use the muon and electron data for sub - gev and multi - gev sample as discussed in section iii .    to settle the basics of our analysis
we should be aware that the predicted atmospheric neutrino flux have a uncertainty of 30 % in the normalization and also a stronger correlation between the fluxes , the relative error between the muon and electron flux is around 5%  @xcite .
as said before , due this normalization error we will test the shape of muon and electron distribution and not the absolute number of events .
we made the following way , we are going to analysis the shape of atmospheric neutrino data and not include the comparison of the absolute value of data . to do this we will change our theoretical prediction of oscillation @xmath191 , where @xmath192 is the normalization parameter with a error of @xmath193 , with @xmath194 for each flavor . also we should use a correlation between the electron and muon number of events due the correlation of the neutrino fluxes .
we made two analysis :    1 .
for the standard oscillation scenario to test our ability to reproduce the results of sk analysis for the oscillation parameters @xmath179 and @xmath195 .
we reproduce the main characteristics of with @xmath196  ev@xmath1 and near maximal mixing ; 2 .
the mavan scenario for non - diagonal parameter @xmath167 , with @xmath197 ,    where our goal is reproduce the angular distribution predicted by sk experiment using our computation made independently of the sk experiment and from this to do a reliable analysis of mavan phenomena .
our choice of @xmath175 function is @xmath198 where @xmath199 is the number of events in the bin @xmath200 measured by sk , @xmath201 is our prediction for the number of events for the bin @xmath200 that depend on the oscillation model used ; the non - diagonal matrix @xmath202 such the diagonal entries have error of @xmath203 as said before , and the non - diagonal entries fixed by the correlation between the muon and electron flux  @xcite .
the sum is over the 40 bins : @xmath204 , counting 2 flavors and 2 samples : sub - gev or multi - gev , bins .
we add the second term in eq .
( [ xq2-new ] ) to introduce a penalty factor when @xmath205 assume values to far from @xmath206 with respect to the error in normalization ,  @xmath207 .    .
summary of our @xmath30 analysis for combined data from super - kamiokande and minos experiment .
the first line corresponds to pure standard neutrino oscillations(s.o . ) and the second line corresponds to mavan scenario . in each case
are show the best fit ( b.f . )
values of the parameters .
[ cols="^,^,^,^,^,^",options="header " , ]     the minos experiment test the mavan scenario for a constant density that implies that the mavan parameters are fixed , and combining with sk experiment that test the mavan for variable density made specially adequate the main hypothesis of mavan idea , to have a mass difference and the mixing amplitude that depends on the local density .
to achieve this we combine the two analysis , from minos experiment and from sk experiment , we will use the @xmath30 test for both experiments super - kamiokande and minos , @xmath208 with this @xmath175 with three oscillation parameters , first we analyze the standard oscillation scenario and second the mavan scenario for this combined analysis of minos and sk .
w show the best fit values for both scenarios in tab .
( [ tab : tab01 ] ) . from this information
we can see that there is a milder improvement of the mavan solution over the standard neutrino oscillation for the combining fit . for the mavan analysis of the combination of minos and sk
, we have found that the best fit is for a non - zero value of mavan parameter ( @xmath209 and the mixing parameters @xmath210  ev@xmath23 and mixing angle @xmath211 and the best fit parameters for standard oscillation are very similar @xmath212  ev@xmath23 and @xmath213 ( see tab .
( [ tab : tab01 ] ) .     for the cases with ( dashed red curve ) and without mavan s effect ( solid black line )
is compared with the result from minos collaboration ( dotted blue curve ) .
we also show the values of best fit points respectively in each analysis . ]
we shown in the _ central panel _ of fig .
( [ fig : deltachi - dmsin-2 ] ) the allowed region of @xmath214 and @xmath188 parameters for the standard oscillation scenario ( s.o . for now ) and for the mavan scenario for the following cases    1 .
first for analysis of the _ standard oscillations ( s.o . )
_ in minos experiment alone shown the dotted blue curve , 2 .
second , for the _ standard oscillations ( s.o . ) analysis _ for the _ combination _ of minos and sk experiments shown by the black curve . we can compare with the minos only result the improvement on the determination of range of @xmath188 for the combination .
the recent values from global fits for these parameters  @xcite are included in our @xmath215 allowed region for @xmath28 and @xmath216 , 3 .
third , for the _ mavan scenario _
analysis for the _ combination _ of minos and sk experiments shown by dashed red curve ,    all plots shown the 90 % c.l . allowed region . comparing the mavan allowed region for minos experiment only shown in fig .
( [ minosonly - regi ] ) and for the combination of minos and sk data ( shown in central panel of fig .
( [ fig : deltachi - dmsin-2 ] ) ) we see that combining the two experiments we constrain more the allowed region of parameters .
another good tool to understand our solution is the plot of the projection of @xmath217 function with respect the one of three oscillation parameters , @xmath218 and @xmath47 .
when we show the projection , e. g. for example @xmath219 we have minimized over the other parameters , @xmath188 and @xmath47 and so on .
we begin with the the _ top panel _ of fig .
( [ fig : deltachi - dmsin-2 ] ) we shown the plot of the projection of @xmath220 for the standard oscillation ( black curve ) and mavan ( dashed red curve ) together with the 66 % , 90 % and 95% c.l .
respectively @xmath221 .
we can see that the combination of minos and sk suppress the high values of @xmath222 that are present in the minos only analysis , they appear only at high c.l .
( not shown ) when we have @xmath223 .
these high values are only possible in minos analysis for @xmath224 , however this large value will induce stronger changes in oscillation probabilities for the sk experiment due larger density that the neutrino feels when cross the earth , typically we have @xmath225 .
we can notice in fig .
( [ fig : pmumu - cz ] ) that for @xmath224 we expect to see a stronger suppression of the muon neutrino oscillation that will conflict the sk data . in the _ right panel _ of fig .
( [ fig : deltachi - dmsin-2 ] ) we shown the plot of the projection of @xmath226 for also both standard oscillation and mavan scenario respectively by black and dashed red curve .
we can see a slightly increase of allowed region for @xmath179 when we have a non - zero mavan parameter .
for @xmath227(dotted - dashed blue line ) , @xmath183(dashed red line ) , and @xmath228(solid green line ) c.l .. black square is the best fit point .
the auxiliary upper plot is the same of fig .
( [ fig : deltachi - dmsin-2 ] ) .
the left auxiliary plot refers to @xmath229 as function of @xmath47 and minimized with respect to the other mixing parameters , @xmath230 and @xmath179 .
]    we show in central panel of fig .
( [ fig : deltachi_aisin ] ) the allowed region of parameters @xmath231 and minimized with respect to @xmath28 at @xmath227 , 90 % and 95 % c.l .
you can see from this plot that there is a correlation between the @xmath184 and @xmath167 for range of values allowed , for the highest @xmath47 value we have the widest range for @xmath184 .
this is left - over of behavior , discussed in section ( [ sec : minos ] ) , that we can have a smaller vacuum mixing amplitude for higher @xmath47 or a large vacuum mixing amplitude for smaller @xmath47 .
in the right panel we show the plot of the projection of @xmath232 where we can see the standard oscillation scenario ( when @xmath51 ) is compatible at @xmath233 and the more interesting information from this plot is that there is no allowed region for @xmath234 at 90 % c.l . .
the solutions for higher @xmath235 present in the minos analysis and only appear now for the combined analysis at @xmath236% c.l .
for the same reasons discussed in the previous paragraph .    at end , to completeness , in central panel of fig .
( [ fig : deltachi_dmai ] ) we show the allowed region of parameters as function of @xmath237 . in this central panel , we can see that the correlation between the values of @xmath238 and @xmath47 are much milder then between the @xmath239 as shown in fig .
( [ fig : deltachi_aisin ] ) .
the reason for this is the change in @xmath162 due @xmath47 it is not so stronger compared with the change in @xmath240 .     for @xmath227(dotted - dashed blue line ) , @xmath183(dashed red line ) , and @xmath228(solid green line ) c.l ..
black square is the best fit point .
the auxiliary upper(right ) plot refers to @xmath229 as function of @xmath28(@xmath47 ) and minimized with respect to other parameters . ]
in this work we test the possibility of neutrinos have their mass dependent on the medium density .
the so - called mavan s model includes this new feature through a new neutrino interaction mediated by a scalar field .
we investigate the consequences of such model in the phenomenology of atmospheric neutrino data from super - kamiokande experiment and for the data from minos experiment using for the first time to test the existence of non - diagonal mavan parameter .    using the fact the mavan for the constant matter density , as it is the case for the minos experiment , have the same functional form of standard oscillation scenario we use the minos analysis for the standard neutrino oscillation and we extend for the mavan scenario . in this analysis larger values of non - diagonal mavan parameter , @xmath241 , are allowed and we have significant changes in the allowed region for oscillation parameters , @xmath214 and @xmath242 .
we compute by ourselves the event rate for super - kamiokande experiment , involves the correct description of the _ sub - gev _ and _ multi - gev _ neutrino energy samples and we also use the analysis of the minos experiment to made a combined fit of these two experiments . in minos experiment the neutrino crosses only a small chord of the earth and then the neutrino feels a constant matter density and for other side the super - kamiokande have neutrinos coming from different directions and therefore feels different medium densities .
then by combining the two data , we can test the essence of mavan hypothesis , the mass dependence on the medium density .
we have found that the best fit values for @xmath0  ev@xmath1 , @xmath2 and mavan parameter @xmath3 and the best fit values for the standard oscillation give similar values(see table  ( [ tab : tab01 ] ) ) .
the allowed region for the oscillation parameters , @xmath214 and @xmath242 is shown in fig  ( [ fig : deltachi - dmsin-2 ] ) that shows that the allowed region is very stable with the addition of mavan scenario .
although the best fit is for non - zero mavan parameter we have all values of mavan parameter @xmath243 are ruled out at 90% c.l that allow us to conclude that the non - diagonal mavan coupling should be give smaller contribuition to neutrino oscillation .
a.  g.  riess _ et al .
_ [ supernova search team collaboration ] , astron .
j.   * 116 * , 1009 ( 1998 ) [ astro - ph/9805201 ] . s.  perlmutter _ et al .
_ [ supernova cosmology project collaboration ] , astrophys .
j.   * 517 * , 565 ( 1999 ) [ astro - ph/9812133 ] .
p.  astier _ et al .
_ [ snls collaboration ] , astron .
astrophys .
* 447 * , 31 ( 2006 ) [ astro - ph/0510447 ] .
p.  gu , x.  wang and x.  zhang , phys .
d * 68 * , 087301 ( 2003 ) [ hep - ph/0307148 ] .
x.  -j .
bi , p.  -h .
gu , x.  -l .
wang and x.  -m .
zhang , phys .
d * 69 * , 113007 ( 2004 ) [ hep - ph/0311022 ] .
n.  afshordi , m.  zaldarriaga and k.  kohri , phys .
d * 72 * , 065024 ( 2005 ) [ astro - ph/0506663 ] .
m.  honda , r.  takahashi and m.  tanimoto , jhep * 0601 * , 042 ( 2006 ) [ hep - ph/0510018 ] .
v.  barger , p.  huber and d.  marfatia , phys .
* 95 * , 211802 ( 2005 ) [ hep - ph/0502196 ] .
m.  cirelli , m.  c.  gonzalez - garcia and c.  pena - garay , nucl .
b * 719 * , 219 ( 2005 ) [ hep - ph/0503028 ] .
m.  c.  gonzalez - garcia , p.  c.  de holanda and r.  zukanovich funchal , phys .
d * 73 * , 033008 ( 2006 ) [ hep - ph/0511093 ] .
p.  c.  de holanda , jcap * 0907 * , 024 ( 2009 ) [ arxiv:0811.0567 [ hep - ph ] ] .
f.  rossi - torres , m.  m.  guzzo , p.  c.  de holanda and o.  l.  g.  peres , phys .
d * 84 * , 053010 ( 2011 ) [ arxiv:1012.5304 [ hep - ph ] ] .
m. f. carneiro and p. c. de holanda , advances in high energy physics , * article i d 293425 * , ( 2013 ) ; available at http://dx.doi.org/10.1155/2013/293425 doi:10.1155/2013/293425 .
k.  k.  shiraishi , _ super - kamiokande atmospheric neutrino analysis of matter - dependent neutrino oscillation models _ , phd .
thesis , at http://www - sk.icrr.u - tokyo.ac.jp / doc / sk / pub/. k.  abe _ et al .
_ [ super - kamiokande collaboration ] , phys .  rev .
d * 77 * , 052001 ( 2008 ) [ arxiv:0801.0776 [ hep - ex ] ] .
p.  adamson _ et al .
_ [ minos collaboration ] , phys .
rev .  lett .
* 110 * , 251801 ( 2013 ) [ arxiv:1304.6335 [ hep - ex ] ] .
m.  c.  gonzalez - garcia , h.  nunokawa , o.  l.  g.  peres , t.  stanev and j.  w.  f.  valle , phys .
d * 58 * , 033004 ( 1998 ) [ hep - ph/9801368 ] .
d. r. gratieri , phd thesis , appendix e , http://webbif.ifi.unicamp.br/tesesonline/teses/if1564.pdf a.  d.  dziewonski and d.  l.  anderson , physics of the earth and planetary interiors * 25 * , 297 ( 1981 ) .
m.  honda , t.  kajita , k.  kasahara , s.  midorikawa and t.  sanuki , phys .
d * 75 * , 043006 ( 2007 ) [ astro - ph/0611418 ] .