the formal expressions of equilibrium statistical mechanics strictly apply only to ergodic systems that are in thermodynamic equilibrium .
thus these expressions only strictly apply to systems which are at the global free energy minimum given the system hamiltonian and the macroscopic thermodynamic state variables ( number of particles , temperature and pressure or density ) . for such systems gibbsian
equilibrium statistical mechanics provides an exact prescription for how to calculate the various thermodynamic quantities @xcite .
however , these prescriptions are routinely applied to systems that are not in true thermodynamic equilibrium ( for example to metastable liquids @xcite , glasses @xcite , polymorphs @xcite and allotropes ) .
it is often observed empirically that within experimental uncertainties many expressions for thermodynamic quantities yield consistent results . in the present paper
we provide arguments for why many of the results of equilibrium statistical mechanics can be applied to such time independent nondissipative nonequilibrium systems .
we also point out some of the limits inherent in the application of the formulae of equilibrium statistical mechanics to such systems .    we choose to study the isothermal isobaric ensemble @xcite ( externally regulated pressure and temperature ) .
the methods and reasoning we use here can be directly transferred to other ensembles such as the canonical ( fixed volume and externally regulated temperature ) .
the gibbs free energy @xmath0 , which is the thermodynamic potential for the isothermal isobaric ensemble , is related to the partition function @xmath1 by the equation @xmath2 and the partition function is given by the integral@xmath3,\label{del0}\ ] ] where @xmath4 is the phase space vector describing the coordinates @xmath5 and momenta @xmath6 , of all the @xmath7 particles in the system , @xmath8 is the thermodynamic pressure , and @xmath9 where @xmath10 is boltzmann s constant and @xmath11 is the temperature .
the integration domain @xmath12 provides limits for both integrals and extends over all the available phase space @xmath13 .
this is @xmath14 for every component of the generalized momentum , @xmath15 for the volume @xmath16 , and over the volume for the cartesian coordinates of the particles . since the system hamiltonian @xmath17 is single valued ,
so too is the partition function and in turn the free energy .    if we require the distribution function of a single thermodynamic phase it is necessary that other phases do not contribute significantly to the partition function .
the full integration domain @xmath12 may include states that are characteristic of crystalline states or fluids states . in the thermodynamic limit
this does not cause problems because , as we shall see , the partition function will be _ completely _ dominated by those microscopic domains that have the lowest free energy . however
the application of these formulae to allotropes or metastable systems does present a problem .
the standard equilibrium statistical mechanical expressions for variables such as the enthalpy @xmath18 , the average volume @xmath19 and second order quantities such as the specific heat at constant pressure @xmath20 may all be computed from a knowledge of the partition function eq .
[ del0 ] or equivalently the thermodynamic potential eq .
[ g0 ] . if other phases of lower free energy exist this computation ( from eq .
[ del0 ] as written ) will strictly speaking be incorrect .
it is well known that the formulae for thermal properties such as entropy , free energy , temperature and specific heat do _ not _ hold for _ dissipative _ nonequilibrium systems outside the linear response regime @xcite . in this paper
we examine the question of whether they are correct for any _ nondissipative _ nonequilibrium systems such as allotropes , metastable systems or history dependent glasses .
we provide a statistical mechanical theory of time independent , nondissipative , nonequilibrium systems .
the theory is based on the fact that these systems are nonergodic and individual sample systems comprise ergodic domains that do not span all of phase space . we show that if these domains are robust with respect to small changes in thermodynamic state variables , a successful statistical mechanical treatment of these nonequilibrium systems can be given .
we provide direct evidence , from molecular dynamics simulations on a model glass former , that the resulting statistical mechanical formulae are satisfied within empirical errors .
finally we provide an independent test of the two key elements of our theory : boltzmann weights within the phase space domains and the robustness of those domains .
it happens that these two elements are the necessary and sufficient conditions for the application of the transient fluctuation relation to finite thermodynamic quenches ( in temperature or pressure ) for such systems @xcite . while the application of thermodynamics to a single time averaged system is usually straightforward the application to a ensemble , whose members may be locked in different phase space domains , can require modification to the standard formulae .    in the case of glasses
our treatment has some similarities with the energy landscape approach of stillinger and weber @xcite .
however there are significant differences ; our treatment makes no reference to the inherent structure and imposes no a priori knowledge of the inter - domain relative population levels .
the energy landscape approach has been extended to account for the phenomena of ageing or history dependence by the addition of a fictive parameter @xcite .
sciortino has convincingly shown that the addition of a single fictive parameter is inadequate to deal with glasses , which may have different properties at the same temperature and pressure if they are prepared by a different protocol ( different history dependence ) @xcite and poses the challenge to recover a thermodynamic description `` by decomposing the ageing system into a collection of substates '' .
the treatment we present here succeeds in doing just that by providing a rigorous development of equilibrium statistical mechanics and thermodynamics for ensembles of systems where the phase space breaks up into ensembles of domains whose inter - domain dynamics is nonergodic and whose inter - domain population levels may not be boltzmann weighted .
a dynamical system in equilibrium has the properties that it is nondissipative and that its macroscopic properties are time independent .
thus the n - particle phase space distribution function , @xmath21 , must be a time independent solution to the liouville equation @xcite,@xmath22 where @xmath23 is the phase space compression factor @xcite obtained by taking the divergence of the equations of motion ( see eq . [ phase - space - compression ] ) and @xmath24 is the extended phase space vector which consists of @xmath25 and may include additional dynamical variables such as the volume @xmath16 . since the system is assumed to be nondissipative both the ensemble average @xmath26 and the time average @xmath27 of the phase space compression factor ( which is directly proportional to the rate at which heat is exchanged with the fictitious thermostat ) are zero .
the time independent solution to eq .
[ liouville - eq0 ] depends on the details of the equations of motion .
equilibrium solutions to eq .
[ liouville - eq0 ] for the equations of motion , suitable for use in molecular dynamics simulations , are compatible with gibbsian equilibrium statistical mechanics @xcite .
microscopic expressions for mechanical properties like the pressure , the internal energy , the enthalpy and the volume can be derived without reference to gibbsian statistical mechanics and indeed can be proved to hold for nonequilibrium systems including nonequilibrium dissipative systems .
there are two ways in which the formulae derived from gibbsian equilibrium statistical mechanics can break down .
the most obvious way is that the relative weights of microstates may be non - boltzmann and the exponential factor @xmath28 $ ] , may be replaced by some other function ( either the exponential function itself may be modified as in tsallis statistics @xcite or the hamiltonian may be modified to some new function @xmath29 ) . in either circumstance
the standard expressions for the thermal quantities derived from equilibrium gibbsian statistical mechanics will not be valid .
this certainly happens in dissipative nonequilibrium systems where the distribution function is not a time independent solution to eq .
[ liouville - eq0 ] .
in deterministic nonequilibrium steady states the phase space may break down into ergodically separated domains ( each of which will be fractal and of lower dimension than the ostensible phase space dimension .
this is a consequence of dissipation . )
however for these steady states , the domains are always exquisitely sensitive to macroscopic thermodynamic parameters since they are strange fractal attractors @xcite .
often a deterministic _ nonequilibrium _ steady state approaches a unique fractal attractor . as time progresses the distribution function collapses ever closer to ( but never reaching ) the steady state attractor .
the second way that these expressions may fail is that the system may become nonergodic . in this case
three things happen .
a ) most obviously time averages no longer equal full ( domain @xmath12 ) ensemble averages .
b ) if we take an initial microstate the subsequent phase space trajectory will span some phase space domain @xmath30 where the initial phase is labeled @xmath31 . in this case for nondissipative nonequilibrium systems where the domains are robust ( i.e. small changes in thermodynamic state parameters , to leading order do not change the domain ) the standard equations of equilibrium statistical mechanics may continue to be valid but in a slightly modified form . we will examine this in some detail below .
c ) given robust domains the population densities between each domain may well depend on the history of the system .
macroscopic _ history can be expected to condition the ensemble s set of initial _ microstates _
@xmath32 from which the macroscopic material is formed .
this in turn can be expected to condition the set of nonergodic domains @xmath33 that characterize the ensemble . for a macroscopic sample spanning a single ergodic domain @xmath30 , the free energy @xmath34 then satisfies only a local extrema principle and thus looses much of its thermodynamic meaning .
we use the constant pressure nos - hoover equations of motion by combining the nos - hoover feedback mechanism with the so - called sllod or dolls equations of motion , which are equivalent for dilation .
it is known that these equations of motion do not produce artifacts in the systems linear response to an external field and that to leading order the effect on the dynamical correlation functions is at most @xmath35 , where @xmath7 is the number of particles @xcite .
the equations of motion are,@xmath36 where @xmath37 is the position , @xmath38 is the momentum and @xmath39 is the force on the @xmath40 particle , @xmath41 is the particle mass , @xmath42 is the barostat time constant , @xmath43 is the thermostat time constant , @xmath11 is the input temperature , @xmath8 is the input ( thermodynamic ) pressure and the instantaneous ( mechanical ) pressure is given by @xmath44 . because these equations of motion have additional dynamical variables the extended phase space vector is @xmath45 . in order to obtain the equilibrium distribution function
we first define the hamiltonian , in the absence of any external fields , dilation @xmath46 or thermostats @xmath47 , as @xmath48 , where @xmath49 is the total inter - particle potential energy .
to proceed further we identify the extended hamiltonian as @xmath50 and then obtain the phase space compression factor@xmath51 where the index @xmath52 sums over the components of the cartesian position and momentum vectors . using the heisenberg streaming representation ( rather than the more usual schrdinger representation eq .
[ liouville - eq0 ] ) of the liouville equation@xmath53=-\lambda\left({\mbox{\boldmath $ \gamma$}}^{\prime}(t)\right),\label{liouville - streaming}\ ] ] we can obtain the particular time independent solution for the distribution function@xmath54 where @xmath55 is the instantaneous enthalpy . the second exponential on the rhs of eq .
[ dist0 ] with @xmath56 and @xmath57 in the argument , which has no dependence on the input temperature @xmath11 or the input pressure @xmath8 , is statistically independent from the rest of the distribution function , which is the standard equilibrium isothermal isobaric distribution .
we can normalize eq .
[ dist0 ] by integrating over all space to obtain the thermodynamic equilibrium distribution function , @xmath58 where the standard isothermal isobaric distribution function is@xmath59 it should be emphasized that the derivation of eq .
[ dist0 ] says nothing about the existence or otherwise of any domains .
these must be considered when normalizing eq .
[ dist0 ] and thus eqs .
[ dist1 ] & [ dist - standard ] are only valid in thermodynamic equilibrium .
if we wish to use eq .
[ dist0 ] outside thermodynamic equilibrium we must consider domains .
we can also use the so called sllod equations of motion @xcite to apply strain rate controlled couette flow ( planar shear ) to our equations of motion .
the necessary modifications to the first two lines of eq .
[ eom ] result in@xmath60 where @xmath61 is the strain rate and the last three lines of eq .
[ eom ] remain unchanged .
as we have stated in the introduction , the full phase space domain includes phase points from many different thermodynamic phases ( gases , liquids and crystals ) . in the thermodynamic limit
this does not cause problems . to understand this suppose we can label microstates to be in either of two possible thermodynamic phases 1 or 2 bound by two phase space domains @xmath62 and @xmath63 . by assumption we are not presently considering the possibility of co - existence .
the system is assumed to be ergodic : atoms in one thermodynamic phase can in time , transform into the other phase .
assume the two thermodynamic phases have different free energies : @xmath64 is the gibbs free energy of the first phase and @xmath65 is that of the second phase . for a sufficiently large @xmath7 the free energy eq .
[ g0 ] is an extensive variable .
we may thus express the partition function as the sum of contributions from the two phases@xmath66 where the lower case @xmath67 on the second line is used to represent the intensive free energies which do not change with system size @xmath7 in the thermodynamic limit .
if @xmath68 is less than @xmath69 then in the thermodynamic limit , @xmath70 , the only significant contribution to the partition function @xmath1 will be due to the `` equilibrium '' phase namely phase 1 .
thus although the free energy defined in eq .
[ del0 ] , is given by an integral over all of phase space @xmath12 , in the thermodynamic limit this integral can be approximated to arbitrary precision , as an integral over the domain that includes the most stable phase .
suppose @xmath62 includes only crystalline phases and @xmath63 includes only amorphous phases and further suppose a particular crystalline phase has a lower free energy that any amorphous phase .
according to eq .
[ del0 ] , we should calculate the free energy by integrating over all crystalline and all amorphous phases . in practice in the thermodynamic limit
we can compute the free energy to arbitrary accuracy by integrating eq .
[ del0 ] , only over that part of phase space within which the thermodynamically stable state resides .    if we consider a nonergodic system that according to different preparative protocols can be formed in either phase 1 or phase 2 .
after preparation , because the system is nonergodic both phases are kinetically stable indefinitely . by restricting the phase space integrals for the free energy to those domains that contain the kinetically stable phase we can compute the free energy of that phase . however , although it may be possible to formally assign free energies to nonergodic systems , these free energies clearly fail to satisfy any global extremum principle .
as we will show these partition functions can be used formally to yield first and second order thermodynamic quantities by numerical differentiation .
the metastable domain is a subset of the thermodynamic equilibrium domain which contains all possible atom positions including ones belonging to the metastable phase .    within a single domain
the system is , by construction , ergodic .
thus for almost all microstates @xmath31@xmath71@xmath30 , ensemble averages , of some variable @xmath72 , @xmath73 , equal time averages @xmath74 , for phase space trajectories that start at time zero , @xmath75 microscopic expressions for mechanical variables may be used as a test of ergodicity in nondissipative systems which are out of equilibrium .
( note a nondissipative system does not on average exchange heat with any thermal reservoir with which it has been in contact for a long time . ) in the case of metastable fluids or allotropes we may introduce a single restricted domain and by construction the system remains ergodic within this domain .
a _ gedanken _ experiment can be used to justify the boltzmann weighting and the applicability of the zeroth law of thermodynamics for such systems .
consider a double well potential with an inner and outer potential well .
if the barrier between the inner and outer wells is much greater than @xmath76 , so that over the duration of observation ( which is much greater than any relaxation time in the ergodically restricted subsystem ) no particles cross the barrier , then the system considered as a double well system , will be , by construction , nonergodic . for systems composed of particles that are solely found in the inner potential well ,
our hypotheses are that the distribution of states in the inner well will be given by a boltzmann distribution taken over the inner domain only and that if such a system is in thermal contact with a body in true thermodynamic equilibrium , then the temperature of the ergodically restricted system must equal that of the system in true thermodynamic equilibrium .
we can justify these hypotheses by considering a fictitious system that only has the inner potential well and in which the potential function is positive infinity for all separations that are greater than the inner well ( this includes the position of the outer well ) . in accord with gibbsian statistical mechanics
the distribution of states is canonical over this ( single well ) potential .
furthermore the zeroth law of thermodynamics will apply to this single well system . now
if we dynamically generate the outer well , all the particles locked inside the inner well can not `` know '' that the outer well has been formed so their dynamics will be completely unchanged by the time dependent generation of the new outer well .
the generation of an inaccessible outer well will not alter the distribution of states in the inner well nor will it cause any flow of heat to the equilibrium heat bath surrounding the system .
this provides a compelling physical justification for our domain hypotheses over a single ergodic sub domain of phase space .    in order to recover many of the basic relationships of gibbsian statistical mechanics
it is also necessary that the system appears to be in dynamical equilibrium , i.e. @xmath77 , for some small @xmath78 , over the longest observation time @xmath79 .
we use the definition of the partition function eq .
[ del0 ] as before but now the domain @xmath30 in the integral is over a single contiguous hypervolume in the configuration space of the generalized position coordinate @xmath5 * * * * and volume @xmath16 .
the domain over the generalized momentum @xmath6 and multipliers @xmath56 and @xmath57 remains unchanged .
we then obtain the gibbs free energy by use of eqs .
[ g0 ] & [ del0 ] .
thus far all we have altered is our definition of the domain . in changing the definition of the domain we have opened a potential problem for gibbsian statistical mechanics .
if we change the temperature or the pressure of the system the domain may also change .
if the domain changes this may make a contribution to the derivatives of the partition function , eq .
[ del0 ] , and the direct connection with the standard outcomes of macroscopic thermodynamics will be lost .
thus the domains need to be robust with respect to changes in thermodynamic state variables .
there are three means by which a system could have robust domains .
the first and most obvious is that the domain does not change when the pressure or the temperature changes , @xmath80 , where @xmath81 is a thermodynamic state variable and @xmath82 is the other thermodynamic state variables .
when we lower the temperature only the inverse temperature @xmath83 in eq .
[ dist0 ] changes and when we change the pressure only the parameter @xmath8 changes .
if the domain s boundary is determined by a surface on which @xmath84 always has a very high value it will remain unchanged under infinitesimal changes in @xmath8 or @xmath83
. we will refer to a surface domain that does nt change with the state variables as completely robust .
the second way is that the distribution function is always identically zero on the domain boundary ,    @xmath85    where @xmath86 is the surface of the domain @xmath30 . because the domain is contiguous ( required for it to be ergodic ) it must have a single connected surface .
such a domain will be robust .
the third way the domain can be robust is less restrictive and allows for the possibility that the domain does change when the thermodynamic variables are changed substantially .
if @xmath87 is an infinitesimal change in a thermodynamic state variable then ,    @xmath88    where @xmath82 denotes the other thermodynamic state variables , we require that @xmath89 for first order thermodynamic property formulae to be correct @xmath90 for second order property formulae to be correct etc . obviously this third way will be satisfied in the first two cases as well .
later in the paper we will introduce an independent test of domain robustness .
however , if a system was not robust then we would expect that small changes in the state variables would change the macroscopic properties of the sample permanently - it would be as though the preparation history of the sample was continuing even for small changes in the state variable .
quite obviously if we produce huge changes in the state variables we will of course permanently change the properties of the system because we permanently deform the ergodic domain .
experience shows however that very many nondissipative nonequilibrium systems are quite robust with respect to small changes in state variables .
all that is required for fluctuation formulae for first second and third order thermodynamic quantities to be valid , is that the domains be unchanged , to first second or third order , by infinitesimal changes in the state variables .
obviously a robust domain is an ideal construct .
however on the typical time scale of interest , which is usually orders of magnitude less than the time scale on which the system will change to a new phase of lower free energy , this can be a very good approximation .
we are now able to recover most of the standard results of gibbsian equilibrium statistical mechanics .
for example we may calculate the enthalpy @xmath91 from the partition function @xmath1 , eq .
[ del0 ] , as@xmath92 here we are considering an ensemble of systems which occupy a single ergodic domain @xmath30 . since this domain is self ergodic the ensemble average is equal to the corresponding time average .
the term on the rhs of the second line is obviously the average value of the instantaneous enthalpy , @xmath93 , obtained by using the equilibrium distribution function , eq .
[ dist - standard ] or equivalently eq . [ dist1 ] with the integration limits restricted to the domain @xmath30 , where @xmath8 is the externally set thermodynamic pressure .
we can also obtain expressions for the average volume @xmath19 and the constant pressure specific heat @xmath20 by taking the appropriate derivatives of the partition function eq .
in other ensembles we can use the same procedure to find other variables e.g. the internal energy , the average pressure and the constant volume specific heat in the canonical ( _ n_,_v_,_t _ ) ensemble .
an important outcome is that this description remains compatible with macroscopic thermodynamics . here
the gibbs free energy is defined as,@xmath94 where @xmath95 is the internal energy and @xmath96 is the entropy . if we take the derivative of eq .
[ gibbs thermo ] with respect to one of the isobaric isothermal ensembles conjugate variables @xmath97 while keeping the others fixed we obtain@xmath98 we now write down the microscopic equilibrium equation for the gibbs entropy@xmath99 it is an easy matter to show that eqs .
[ g0 ] , [ del0 ] & [ gibbs entropy ] are consistent with the two derivatives given in eq .
[ gibbs conj deriv ]
. given our condition of a robust boundary we thus have a form of gibbsian statistical mechanics for metastable states which remains in agreement with macroscopic thermodynamics .
we now wish to consider an ensemble of systems which is prepared from an initial ergodic ( usually high temperature ) equilibrium ensemble .
there is some protocol @xmath100 , which involves a temperature quench or a sharp pressure increase etc , which breaks the ensemble into a set of sub - ensembles characterized by different macroscopic properties .
after the protocol @xmath100 , has been executed we allow all the ensemble members to relax to states which are macroscopically time independent - to within experimental tolerances .
we assume that the ensemble can be classified into a set of sub - ensembles @xmath101 whose macroscopic properties take on @xmath102 distinct sets of values .
for the longest observation times available a macroscopic system classified as an @xmath103 system is not observed to transform into a @xmath83 system , and vice versa .
the full ensemble of systems is thus non - ergodic .
however , in each individual sub - ensemble , say sub - ensemble @xmath103 , the constituent members are ergodic ( by construction ) .
thus we can partition the full phase space into a set of domains , @xmath33 .    from the arguments given above ( in section b ) , after the relaxation of initial transients
, we expect to observe a boltzmann distribution of states within an individual domain which is therefore independent of the quench protocol .
however the distribution between domains can not be expected to be boltzmann distributed and will instead be dependent on the quench protocol . within a given ensemble the proportion of ensemble members ultimately found in domain @xmath30
is given by a weight @xmath104 which is subject to the constraint @xmath105 we can calculate the full _ ensemble _ average of some macroscopic property @xmath72 as,@xmath106 since the full ensemble of states is nonergodic the phase space breaks up into disjoint domains which in themselves are ergodic .
thus each domain may be identified by any point in phase space , @xmath13 , that is a member of it so the subscript @xmath103 is a function of the phase vector @xmath107 allowing the following expression for the distribution function@xmath108 where @xmath109 if @xmath110 and @xmath111 otherwise , and @xmath112 the entropy is given by @xmath113 and using eq .
[ multi - d - dist ] we have the following expressions for the multidomain entropy,@xmath114=\sum_{\alpha=1}^{n_{d}}w_{\alpha}s_{\alpha}-k_{b}\sum_{\alpha=1}^{n_{d}}w_{\alpha}\ln(w_{\alpha}).\label{s - multi - d}\ ] ] the term @xmath115 is the inter - domain entropy , which is maximized by an even distribution of ensemble members over all domains , while @xmath86 is the intra - entropy of domain @xmath103 considered as a single _ n_-particle system .
if we substitute eq .
[ multi - d - dist - raw ] into eqs .
[ s - multi - d ] we find that ,    @xmath116    where @xmath117 .
combining eq .
[ gibbs thermo ] with eq .
[ s - multi - d - g ] we obtain the following expression for the gibbs free energy    @xmath118=\sum_{\alpha=1}^{n_{d}}w_{\alpha}g_{\alpha}-s_{d}t.\label{g - multi - d}\ ] ]    it is easy to verify that if we hold the local domain weights fixed and then vary the temperature or the pressure that eqs .
[ s - multi - d ] & [ g - multi - d ] are compatible with eqs .
[ gibbs conj deriv ] .
this means that if we have a fixed number of robust domains , whose population levels or weights are non - boltzmann distributed , eqs .
[ g - multi - d ] & [ s - multi - d ] provide a direct microscopic link to standard macroscopic thermodynamics . on the extremely long time scale
the weighting functions @xmath119 may vary and the system will tend towards the direction where the free energy eq . [ g - multi - d ] is reduced . without the inter - domain entropy term ,
@xmath120 , eq . [ g - multi - d ] would be minimized when the domain with the lowest free energy has all the ensemble members in it .
it turns out that eq .
[ g - multi - d ] is minimized when all the domain weights are boltzmann distributed , ie when@xmath121 here ( i.e. upon obeying eq .
[ weights boltzmann ] ) the entropy and free energy given by eqs .
[ s - multi - d ] & [ g - multi - d ] coincide with the standard equilibrium expressions so the free energy must be a minimum . to prove this we use eq .
[ g - multi - d ] and we remove the first weight @xmath122 , so that the constraint eq .
[ w - constraint ] is respected while the remaining weights are independent .
this means that the free energy can be written as @xmath123 .
the constrained partial derivatives are then @xmath124 using the fact that for boltzmann weights , eq . [ weights boltzmann],@xmath125\label{weights - boltz-2}\ ] ] where @xmath126 is the equilibrium free energy , we find that at equilibrium ,    @xmath127    it remains to prove that this is indeed a minimum . using the same approach for treating the constraint , we continue making the first weight a function of all the others , and obtain@xmath128 using the boltzmann weights it is easy to show that@xmath129\nonumber \\ \frac{\partial^{2}g}{\partial w_{1}\partial w_{\alpha } } & = & \frac{\partial^{2}g}{\partial w_{\gamma}\partial w_{1}}=0\nonumber \\ \frac{\partial^{2}g}{\partial w_{\gamma}\partial w_{\alpha } } & = & \delta_{\gamma\alpha}\frac{k_{b}t}{w_{\alpha}}=\delta_{\gamma\alpha}k_{b}texp[\beta(g_{\alpha}-g_{eq})].\label{hessian-2}\end{aligned}\ ] ] from these results it is easy to see that the hessian matrix @xmath130 is positive definite , thus we have proved the free energy to be a minimum for the case of equilibrium .
we can use a knowledge of the multiple domain thermodynamic potential , eq .
[ g - multi - d ] , to compute averages .
as an example we consider the average enthalpy again,@xmath131 eq . [ enthalpy - multi - d ] can easily be derived from eq .
[ gibbs thermo ] and is in essence the same as the first line of eq .
[ enthalpy ] .
it is straightforward to see that upon using eq .
[ g - multi - d ] to calculate the derivative in eq .
[ enthalpy - multi - d ] one obtains the average enthalpy as given by eq .
[ multi - d - average ] .
one can do the same for other quantities such as the specific heat etc .
we are now in a position to test various outcomes , using computer simulation , which may be derived by drawing on the previous material . if we start with an ensemble of systems which are initially in equilibrium at temperature @xmath132 and then at time @xmath133 we quench them by setting @xmath134 we can solve the liouville equation eq .
[ liouville - streaming ] to obtain@xmath135\label{quench dist}\end{aligned}\ ] ] where @xmath136 .
this nonequilibrium distribution function , valid for @xmath137 , explicitly requires the solution of the equations of motion and is of limited utility .
however it allows the identification of a formal condition to identify the amount of time , which must elapse after the quench , before eqs .
[ dist0 ] , [ dist - standard ] , [ multi - d - average ] or [ multi - d - dist ] can be applied to the ensemble .
that is the quantity @xmath55 must be statistically independent of @xmath138 .
thus we are interested in the correlation function @xmath139 where @xmath140 this function will equal 1 for a perfectly correlated system , -1 for a perfectly anticorrelated system and 0 for an uncorrelated system .
when we consider an ensemble of systems , occupying the various domains to different degrees , we see that eq .
[ correlation1 ] may not decay to zero given the trajectories are unable to leave their domains . if the transients , due to the quench , fully decay eq .
[ correlation1 ] will decay to zero and eqs .
[ dist0 ] & [ dist - standard ] will become valid for the ensemble .    if the correlation function , eq .
[ correlation1 ] decays to a plateau then we may have a situation where eqs .
[ multi - d - average ] and [ multi - d - dist ] are valid .
it may be that the material can still slowly age , due to processes , that occur on a time scale which is longer than the one we are monitoring . for a glass
we expect that correlation function eq .
[ correlation1 ] will not fully decay on a reasonable time scale .
if we give the system time to age , such that it appears to be time translationaly invariant , and then compute the following correlation function,@xmath141 where@xmath142 we may observe a full decay . if this occurs eqs .
[ multi - d - average ] and [ multi - d - dist ] will be valid . to compute this
correlation function @xmath143 trajectories are produced and for each of these the averages , @xmath144 , where @xmath72 is an arbitrary dynamical variable , appearing in eq .
[ correlation2 ] are approximately obtained by time averaging . when the system is ergodic and time translationaly invariant these two correlation functions , eqs .
[ correlation1 ] & [ correlation2 ] , will give the same result .
however for a nonergodic system @xmath145 will decay to zero on a reasonable time scale while @xmath146 will not .
rather @xmath146 may decay to some plateau on a reasonable time scale and then decay on a much slower time scale .
the preceding sections then rest on this clear separation of time scales in the correlation function eq .
[ correlation1 ] . for metastable fluids and allotropes
this separation is so extreme that we probably can not observe , even the early stages of , the later slow decay on any reasonable experimental time scale .
further for these systems there will be only a single domain and thus they appear ergodic . for glasses
some signs of the later decay can often be observed , however it is still very much slower than the initial decay . in the field of glassy dynamics the initial decay
is often called the @xmath83 decay and the slower long time decay is often called the @xmath103 decay . as the glass is further aged this second stage decay
is observed to slow down dramatically while the early decay does not change very much @xcite .    to allow the hypothesis of local phase space domains to be tested we will now introduce several relations whose derivation draws upon the equilibrium distribution function eq .
[ dist - standard ] .
we will also discuss the effect of the phase space domains on these relations .
first we introduce the configurational temperature@xmath147 where @xmath148 is a @xmath149 dimensional vector representing the inter - particle forces on each atom @xmath150 .
this relation is easily derived from eq .
[ dist - standard ] ( set @xmath72 to @xmath151 , integrate by parts and drop the boundary terms ) and will remain valid for nonergodic systems where eq .
[ weights boltzmann ] is not obeyed .
the relation involves the spatial derivative of the force , which is not zero at the cutoff radius for the potential we use in our simulations ( see below ) .
this along with finite size effects can result in a small disagreement between eq .
[ config temp ] and the kinetic temperature for our system when in equilibrium .
equilibrium fluctuation formulae may be easily derived from eq .
[ multi - d - average ] , see ref .
@xcite for some examples of how this is done .
we consider how the average enthalpy changes with the temperature at constant pressure,@xmath152\label{cp - regular}\ ] ] where @xmath153 is the constant pressure specific heat .
the calculation of the rhs of eq .
[ cp - regular ] using the ensemble average given by eq .
[ multi - d - average ] is subtle . if we simply calculate @xmath154 and @xmath155 with the use of eq .
[ multi - d - average ] and then plug the results into eq .
[ cp - regular ] we obtain , what we will refer to as , a _ single domain average _ which does not give us the correct change in the average enthalpy for a history dependent equilibrium ensemble .
this is because the average @xmath155 superimposes across the different domains while the quantity @xmath156 contains spurious cross terms .
if we derive the heat capacity by taking the second derivative of the gibbs free energy for the multiple domain distribution function , eqs .
[ g - multi - d ] & [ enthalpy - multi - d ] , with respect to the temperature , @xmath157 we obtain the following    @xmath158,\label{cp - concat}\ ] ]    where the quantity @xmath159 is obtained for each domain separately ( here @xmath160 represents an average taken where all ensemble members are in the @xmath161 domain ) .
we will refer to this as a _ multidomain average _ which is consistent with the nonergodic statistical mechanics and thermodynamics that we have introduced here .
it is obvious that both a single and multiple domain average will give the same result in the case of thermodynamic equilibrium and metastable equilibrium ( single domain ) .
the transition from the single domain average producing the correct result to an anomalous result is symptomatic of an ergodic to a history dependent nonergodic transition .
if we consider a large macroscopic system ( the super system ) to be made of @xmath162 independent subsystems the multidomain average remains self - consistent .
to see this we apply eq . [ cp - concat ] to fluctuations in the super system
and then we inquire how this relates to fluctuations in the subsystem .
the enthalpy of one instance of the super system will be given by @xmath163 . in principle
an ensemble of super systems can be prepared by applying the same history dependent macroscopic protocol to all members of this ensemble . due to the statistical independence of the subsystems , upon taking an ensemble average of super systems , we have @xmath164 for all @xmath165 . here
the average @xmath166 is taken over the ensemble of super systems and the @xmath161 subsystem in each super system is identified by its location .
it is then easy to show that the specific heat obtained from the ensemble average eq .
[ cp - concat ] of the super system is equivalent to that obtained from the subsystem due to the two quantities @xmath167 and @xmath168 possessing identical cross terms which cancel each other out ( as a result of the independence of the subsystems ) upon applying eq [ cp - concat ] .
if we ignore finite size effects , due to assuming the equivalence of ensembles , the constant volume specific heat is related to the constant pressure specific heat by the equation@xmath169 we may also obtain an expression for the constant volume specific heat @xmath170 by deriving equilibrium fluctuation formula for each of the derivatives appearing in eq .
[ cv ] in lieu of directly measuring them .
we may then obtain a single domain expression for @xmath170 which does not work for the history dependent glass and also a correctly weighted ensemble average ( multidomain average ) which does .
this is completely analogous to what has been shown in detail for @xmath20 .
as the equations are unwieldy , and their derivation ( given an understanding of the @xmath20 case ) is straightforward , we will not reproduce them here .
the application of the evans - searles transient fluctuation theorem @xcite to the systems treated in this paper provides a sharp test of the assumptions used to develop the theory given in this paper .
the theorem describes a time reversal symmetry satisfied by a generalized entropy production , namely the so - called dissipation function .
the precise mathematical definition of this function requires a knowledge of the dynamics and also of the initial distribution function .
the three necessary and sufficient conditions for the fluctuation theorem to be valid are that the initial distribution is known ( here we assume the distribution is boltzmann weighted over some initial domain of phase space ) , that the dynamics is time reversible ( all the equations of motion used here are time reversible ) and lastly that the system satisfies the condition known as ergodic consistency .
when applied to the systems studied here this requires that the phase space domains should be robust with respect to the sudden changes imposed on the system and that the number of inter - domain transitions remain negligible on the time over which the theory is applied .
if any one of these three conditions fails then the theorem can not be applied to the system and the corresponding fluctuation relation will not be satisfied@xcite .
we can use the fluctuation theorem to obtain relations for how the system responds upon suddenly changing the input temperature or pressure for a system , which is initially in equilibrium as specified by eq .
[ dist0 ] or [ dist1 ] .
firstly we consider a change in the pressure , while holding the temperature fixed , by changing the input variable @xmath8 in eq .
[ eom ] ( thermodynamic pressure ) to @xmath171 at time @xmath133 for a system initially in equilibrium with @xmath172 .
the probability density @xmath173 of observing a change in volume of @xmath174 relative to a change of equal magnitude but opposite sign is then given by@xmath175 to derive this expression we have had to assume that the intra - domain populations are boltzmann distributed according to eq .
[ dist0 ] .
ergodic consistency requires that for any initial phase space point @xmath176 that can be initially observed with nonzero probability , there is a nonzero probability of _ initially _ observing the time reversal map @xmath177of the end point @xmath178 , ( ie @xmath179 such that @xmath180@xmath181@xmath182 ) .
this condition obviously requires that the phase space domains remain robust and the number of inter - domain transitions remain negligible for at least a time @xmath183 , after the pressure ( or temperature ) quench .
if we sample all or our initial @xmath133 , @xmath172 configurations from the one trajectory which remains locked in a single domain even after the quench , we expect eq .
[ tft - p ] to be valid .
if we prepare an ensemble of initial configurations using the same protocol we still expect eq .
[ tft - p ] to remain valid even with different domain weightings @xmath184 , as defined in eqs .
[ multi - d - average ] & [ multi - d - dist ] , provided the domains are robust over the time @xmath183 appearing in eq .
[ tft - p ] .
note that suddenly reducing the pressure by a very large amount could result in a breakdown of the robustness condition .
[ tft - p ] may be partially summed to obtain what is referred to as the integrated fluctuation theorem @xmath185    for the case where we change the input temperature @xmath11 in eq .
[ eom ] while holding the input pressure @xmath8 fixed we obtain a relation for fluctuations in the extended instantaneous enthalpy @xmath186 .
we start with a system initially in equilibrium at temperature @xmath187 and we then subject it to a temperature quench by changing the input temperature in eq .
[ eom ] to @xmath188 , at time @xmath133 , while holding the input pressure fixed .
the probability density @xmath189 of observing a change in instantaneous enthalpy @xmath190 relative to a change of equal magnitude but opposite sign is then given by@xmath191 note that if we suddenly increase the temperature by a very large amount we could expect to violate the robustness or the negligible inter - domain transition condition . in common with eq .
[ tft - p ] we expect that this expression will be valid when all initial configurations are sampled from a single common domain and also when sampled from multiple arbitrarily populated domains under the assumption that the domains are robust and the number of transitions are negligible over time @xmath183 .
this equation may also be partially summed to obtain@xmath192
for our simulations we use a variation on the kob andersen glass former @xcite featuring a purely repulsive potential @xcite . the pairwise additive potential is@xmath193 & \;\;\;\forall\;\;\ ; r_{ij}<\sqrt[6]{2}\,\sigma_{\alpha\beta}\nonumber \\
u_{ij}(r_{ij})=0\;\;\;\;\;\ , &   & \;\;\;\forall\;\;\ ; r_{ij}>\sqrt[6]{2}\,\sigma_{\alpha\beta},\label{potential}\end{aligned}\ ] ] where the species identities of particles @xmath194 and @xmath195 , either @xmath196 or @xmath72 , are denoted by the subscripts @xmath103 and @xmath83 .
the energy parameters are set @xmath197 , @xmath198 and the particle interaction distances @xmath199 , @xmath200 . the energy unit is @xmath201 , the length unit is @xmath202 and the time unit is @xmath203 with both species having the same mass @xmath41 .
the composition is set at @xmath204 , the number of particles are @xmath205 , the pressure is set to @xmath206 and the temperature unit is @xmath207 .
the time constants are set at @xmath208 and @xmath209 .
note that the energy parameters are slightly different to the potential we used in ref .
the equations of motion were integrated using a fourth order runge - kutta method @xcite .
the time step used was @xmath210 and sometimes @xmath211 for very low temperatures .    from previous work on binary mixtures we know the basic reason why this system is vary reluctant to crystallize @xcite .
the chosen nonadditivity of the species _
interaction makes the mixture extremely miscible ; consider the present value of @xmath200 relative to the additive value of @xmath212 .
this effect dominates over the choice of the energy parameters . due to this extreme miscibility
the relatively large composition fluctuations necessary , about the average composition of @xmath213 , to form the crystal phases ( either the pure species _ a , _
@xmath214 , fcc crystal or the binary , @xmath215 , cscl crystal ) are strongly suppressed and crystallization is strongly frustrated .     with @xmath216.[fig-1 ] ]    we identify the nominal glass transition by calculating the diffusion coefficient as a function of temperature .
this is shown for both species in fig .
[ fig-1 ] on a logarithmic plot demonstrating how the diffusion coefficient approaches zero critically , @xmath217 where @xmath218 is the critical exponent , with a nominal glass transition temperature of @xmath216 .
it would be , perhaps , more customary to obtain a nominal glass transition temperature by analyzing the critical divergence of the viscosity .
given that the stokes einstein relation is strongly violated upon approaching the glass transition one might be concerned that this would give a very different result .
however the violation of the stokes einstein relation can largely be attributed to the exponent @xmath218 being different between the viscosity and the diffusion coefficient rather than the nominal glass transition temperature @xmath219 @xcite .
the correlation functions given in eqs .
[ correlation1 ] & [ correlation2 ] were calculated from ensembles of @xmath220 independent simulations at the two temperatures given in fig .
[ fig-2 ] ( @xmath221 and @xmath222 ) . in all cases
the systems were subject to an instantaneous quench , from an initial equilibrium at @xmath223 , by changing the value of the input temperature in eq .
the system was then run for a significant time , in the case of the glass ensemble @xmath224 , in an attempt to age it .
of course the longest time that can be accessed in a molecular dynamics simulation is rather short , and so the system is not very well aged , but we are still able to meaningfully treat it as a time invariant state . each of the @xmath220 independent simulations was interpreted as being stuck in its own domain @xmath30 and the correlation functions were calculated for each of these domains using time averaging .
the time averaging was approximately @xmath220 times longer than the longest time @xmath225 that the correlation functions were calculated out to .
obviously in the limit of an infinite number of independent simulations and the case where the domains are robust we will obtain the exact multidomain average given by eq .
[ multi - d - average ] .
we assume our limited ensemble of simulations is representative of this .
the data from each domain ( independent simulation ) was then used to obtain the correlation functions eqs.[correlation1 ] & [ correlation2 ] as seen in fig .
[ fig-2 ] .     &
[ correlation2 ] as a function of logarithmic time for the temperature @xmath221 .
the calculation of the function was started after the fluid had been given time to equilibrate .
the strong agreement between the two correlation functions is indicative of ergodicity .
+ b ) the instantaneous enthalpy correlation functions as defined by eqs .
[ correlation1 ] & [ correlation2 ] as a function of logarithmic time for the temperature @xmath222 .
the calculation of the correlation function was started at various times after the quench , all approximately at @xmath226 , in an attempt to age the system .
the difference between the correlation functions is indicative of nonergodicity .
notice that @xmath145 reaches a near full decay between @xmath227 , while @xmath146 reaches a non - decaying plateau .
[ fig-2 ] ]    at the higher temperature @xmath221 it can be seen that the two functions are equivalent demonstrating how the system is ergodic .
it can also be seen that the correlation function has decayed on a time scale of @xmath228 which is therefore ( by eq .
[ quench dist ] ) the time scale on which the ensemble becomes accurately represented by eq .
[ dist - standard ] with only one domain @xmath12 which does not necessarily extend over all phase space .
the oscillations , which can be seen in the correlation function , are due to both ringing in the feedback mechanisms of eq .
[ eom ] and the frequency dependent storage component of the bulk viscosity . the statistical uncertainty in the correlation function becomes larger than these oscillations somewhere between a time of @xmath227 .
if we constructed an experiment where the pressure was regulated by a piston and a spring , the correlation functions , eqs.[correlation1 ] & [ correlation2 ] , would depend on the details of the piston and spring parameters in a similar way to the simulations dependence on the details of the feed back mechanism .
when the system is quenched to the lower temperatures ( @xmath222 ) ergodicity is lost and we obtain a glass . the complete decay of the first correlation function , @xmath146 , eq . [ correlation1 ] may no longer occur because the individual trajectories remain stuck in local domains which have different average values for @xmath229 .
this is similar to the much studied density correlation function @xcite ( the intermediate scattering function ) which decays to a finite plateau for a glass or more generally a solid material . on the other hand
there is nothing to stop the second correlation function , @xmath145 , eq .
[ correlation2 ] , from fully decaying when the system is nonergodic .
if the second correlation function fully decays whilst the first is only able to decay to a plateau then we have a situation where eqs .
[ multi - d - dist ] & [ multi - d - average ] are valid as can be seen from eq .
[ quench dist ] . in fig .
[ fig-2 ] it can be seen that @xmath146 does indeed fail to decay while @xmath145 comes very close to fully decaying at the time where @xmath146 reaches the plateau .
the reason @xmath145 does nt fully decay here is due to the fact that the inter - domain transition rates , while small , are not exactly zero .
if the system had been aged more extensively this problem would be significantly reduced .
this effect is exacerbated by the time averaging , used to form the averages for each trajectory , being two orders of magnitude longer than the longest time the correlation function was calculated to .
the effect of the state slowly evolving due to finite inter - domain transition rates is too small to seriously compromise the modeling of the system as obeying eqs .
[ multi - d - average ] & [ multi - d - dist ] and thus we have obtained direct evidence for the validity of these equations . the height of the plateau for @xmath146 will depend on the history of the system , i.e. the protocol used to prepare the ensemble .
we move on to a comparison between the kinetic temperature and the configurational temperature , the results of which may be seen in fig . [
fig-3]a ) .
( controlled directly by the nos hoover thermostat ) and the configurational temperature given by eq . [ config temp ] as a function of the kinetic temperature .
the results for the system undergoing couette flow ( shear ) are for a constant strain rate of @xmath230 .
+ b ) the heat capacity calculated using the equilibrium fluctuation formula by the single domain averaging method eq .
[ cp - regular ] and the multidomain method eq .
[ cp - concat ] . also shown is data obtained by numerically differentiating the enthalpy by central difference . at the temperatures above the peak the three types of averages
give very similar results .
also shown is equivalent equilibrium fluctuation formula data for the constant volume specific heat eq .
[ fig-3 ] ]    the input temperature ranges from @xmath231 to @xmath232 . also shown are results for the system , undergoing constant planar shear , eq . [ sllod ] , with a strain rate of @xmath230 . at the higher temperatures we see a very small relative discrepancy between the two types of temperatures , which we attribute to the discontinuity in the first spatial derivative of the inter - particle force at the cutoff radius and to finite size effects .
these effects appear to diminish a little at lower temperatures . for temperatures above t = 1.5
the chosen strain rate has no significant effect on the configurational temperature indicating that our system is in the linear response domain @xcite .
at the lowest temperatures , well below the glass transition temperature , we observe good agreement between the configurational and input temperatures for the system without shear .
this provides further evidence of our assertion that the system obeys boltzmann statistics in the glass eqs .
[ multi - d - average ] & [ multi - d - dist ] .
however , at low temperatures , the system that is undergoing shear shows an increasing relative discrepancy between the two temperatures . at low temperatures
the system leaves the linear response domain @xcite demonstrating the fundamental difference between the nonequilibrium distribution of the history dependent glassy state and that of a strongly driven steady state .
if we drive the system hard enough , at any given temperature , we can always make a disagreement between the two types of temperature due to the steady state no longer being accurately represented by a boltzmann distribution i.e. due to a break down in local thermodynamic equilibrium .
when the system is not driven by an external field we have been unable to observe any difference in the two temperatures by deeply supercooling a glass forming mixture apart from the initial transient decay immediately following the quench , which falls off surprisingly rapidly .    in fig
[ fig-3]b ) results are presented for the heat capacities ( the specific heat multiplied by the volume ) at both constant pressure @xmath233 and constant temperature @xmath234 .
details of the protocol used to obtain this data is given in the end note @xcite .
the estimates from the multidomain average are compared with those from the single domain average .
the results from the multidomain averages exhibit the well - known peak , which is a signature of the onset of the glass transition , and has been observed directly by calorimetry in many experiments on real glass forming materials @xcite .
the temperature , where the peak is observed , depends on the history of the system .
no peak is observed for the single domain averages which continue to increase as the temperature is lowered . while the two methods for forming averages give the same results at temperatures above the peak , they diverge at temperatures below the peak .
it is the multidomain average that gives results consistent with the actual calorimetric behavior of the system .
this may be seen in the figure by comparing the data which has been computed by numerically differentiating the enthalpy using central difference . at the peak
neither the central difference ( due to rapid rate of change ) or the multidomain average data ( due to a lack of domain robustness ) are reliable and they show significant differences .
however below the peak they once again show quantitative agreement providing strong evidence that the domains are robust in this region .     and the solid line is@xmath235 and b ) a sudden temperature change , where the symbols are @xmath236 and the solid line is@xmath237 .
results from an ensemble of independent initial systems and in addition from a single initial trajectory ( with a time of @xmath238 being computed between each transient trajectory ) are shown for a total of @xmath239 pressure or temperature changes .
the serial results have been shifted up , for clarity , by adding 0.2 to the data [ fig-4 ] ]    if we substantially increase the duration the time averages ( for each domain ) are constructed on , the peak will be shifted to lower temperatures as previously shown @xcite .
this requires the time average to be constructed over some two orders of magnitude more time than the decay time for the correlation function in eq .
[ correlation1 ] .
this is necessary in order to obtain enough independent samples for a meaningful estimate of the variance , of the instantaneous enthalpy appearing in eq .
[ cp - regular ] . for a large macroscopic system
we would expect that the specific heat measured over the entire ensemble would differ very little to that measured from any one of its members .
we are now in a position to make an unambiguous interpretation of the peak in the specific heat .
the peak is observed at the temperature where the system leaves metastable equilibrium and enters a history dependent state that requires averages to be computed by eq .
[ multi - d - average ] rather than by direct use of eq .
[ dist - standard ] .
the calculation of both @xmath229 and @xmath240 will be different for each domain . if we use time averaging to calculate these variables on a time scale that falls within the plateau region for eq .
[ correlation1 ] , see fig .
[ fig-2]b , then the amount of time chosen to form the average is not critical .
the peak occurs because the various ensemble members have become locked in local domains on the time scale that we are able to access .
near the peak itself these domains are not expected to be robust .    the multidomain average , eq .
[ cp - concat ] , gives the heat capacity for a glass with robust domains . at the lowest temperatures the heat capacity reaches the beginning of a plateau , fig .
[ fig-3]b . for the constant volume heat capacity @xmath234 this plateau ( within uncertainties due to finite size effects ) has a value that is consistent with the dulong - petit law @xcite , as would be expected for an amorphous solid where the potential energy surface can be modeled as harmonic upon transformation to the orthogonal independent basis set .
this is exactly what we would expect from our local domain model at low temperatures .    testing the integrated transient fluctuation theorem ( itft ) for a sudden pressure change eq .
[ ift - p ] and a sudden temperature change eq .
[ ift - t ] provides further evidence that the boltzmann distribution may be used to accurately describe intra - domain statistics in the glassy state and also that in a properly aged glass , the domains are robust with respect to the pressure and temperature changes studied here , fig .
[ fig-4 ] .
these equations remain valid whether we subject an ensemble of simulations ( multidomain ) to a quench or we sample from a single trajectory ( single domain ) , which remains stuck in a single domain .
the accuracy with which these relations are satisfied is powerful independent evidence for the applicability of our assumptions to the systems studied here .
the fact that over the times shown in fig .
[ fig-4 ] , the itft does indeed yield correct results directly implies that , within experimental tolerance of the data , the phase space domains must be robust and the number of inter - domain transitions must be negligible . unlike the specific heat fluctuation formula this requires that the domains are robust to finite changes of the state rather than infinitesimal changes .
thus , given that we have aged the glass sufficiently that domains are robust the number of transitions are negligible over the longest time the fluctuation formulae are computed , we expect eqs .
[ ift - p ] & [ ift - t ] to be correct .
if we wished to apply the steady state fluctuation theorem matters become more difficult @xcite .
we have presented a rigorous development of statistical mechanics and thermodynamics for nonergodic systems where the macroscopic properties are sensibly time independent and the phase space for the ensemble is partitioned into robust domains .
using computer simulation we have carried out various tests on a glassy system and shown that apart from the immediate vicinity of the glass transition , the computed results are consistent with our theory .
while the intra - domain populations are individually boltzmann distributed , the inter - domain populations are not .
a correlation function whose decay to zero requires global boltzmann weighting , has been derived and it has been shown that it decays on a reasonable time scale for ergodic systems but not for nonergodic systems . a second correlation function which decays to zero if the intra - domain populations are boltzmann distributed but globally the inter - domain populations are not , has also been defined .
we have developed expressions for obtaining averages in a multiple domain ensemble and shown how single domain averages , which always give correct results in metastable equilibrium , can give spurious results in a history dependent nonergodic ensemble .
the statistical mechanics and thermodynamics developed here allow the derivation of expressions for multidomain ensemble averages which give the correct results for time nondissipative nonequilibrium ensembles .
the fundamental origin of the peak in the specific heat near the glass transition has been unambiguously shown to be a signature of a transition from metastable equilibrium to a nonergodic multi - domain ensemble .
we have shown that the transient fluctuation relations for temperature and pressure quenches provide independent tests of the fundamental hypotheses used in our theory : that intra - domain populations are individually boltzmann distributed , that except in the immediate vicinity of the glass transition the domains are robust with respect to small but finite variations in thermodynamic state variables , and that the inter - domain transition rates are negligible .