soon after the experimental discovery  @xcite of carbon nanotubes , it was suggested  @xcite that there might also exist carbon _ nanotori _  , carbon nanotubes in which the ends of the tube are identified and sewn together . within a few years ,
experimental evidence for such structures emerged  @xcite , and since then they have attracted a great deal of experimental  @xcite and theoretical  @xcite attention .
there are a variety of reasons for this intense interest .
for example , certain species of carbon nanotori exhibit unusual magnetic properties  @xcite , including persistent magnetic moments at nearly zero flux  @xcite and colossal paramagnetic moments  @xcite .
these objects can also display a diverse variety of electric properties : some nanotori are inherently metallic , while others are semiconducting and still others are insulators .
as we shall discuss , the most general nanotorus can be parametrized by four integers @xmath0 . together and in various combinations , these describe the radius of the underlying nanotube , the chirality of the underlying nanotube , the length of the underlying nanotube , and the relative twist  @xcite that might occur upon sewing opposite ends of the tube together to form the torus .
the important point , however , is that nanotori with different @xmath0 are fundamentally physically distinct : they have different sizes , different shapes , and different twisted `` honeycomb '' patterns of carbon atoms laid out on their surfaces .
of course , there are certain symmetries of the underlying graphene sheet which lead to trivial equivalences amongst these nanotori .
for example , @xmath1 rotations of the underlying graphene sheet will produce identical nanotori .
such nanotori are therefore not physically distinct .    in this paper , however , we shall demonstrate that there are additional symmetries that relate _ physically distinct _
nanotori to each other , forcing such nanotori to exhibit identical energy spectra and electrical properties .
these so - called `` modular '' symmetries therefore transcend the traditional hexagonal lattice symmetries of the graphene sheet , and arise solely in the process of compactifying the graphene sheet in order to form the nanotorus . in some sense
, the appearance of these modular symmetries is entirely expected , for they result directly from the geometry of the compactification . however , to the best of our knowledge , the significance and consequences of these symmetries have not been fully appreciated thus far in the nanotorus literature .
as we shall demonstrate , these symmetries have a number of profound effects .
first , we shall see that use of these symmetries allows us to partition the set of carbon nanotori into distinct equivalence classes as far as their spectral properties are concerned , and leads to a dramatic reduction in the numbers of _ spectrally distinct _
carbon nanotori as compared with the numbers of _ physically distinct _ nanotori .
moreover , as we shall show , many of the standard rules of thumb advanced in the literature in order to describe the conductivity properties of these nanotori actually fail to respect these symmetries
. such rules of thumb are therefore incomplete as descriptions of the physics of these nanotori , and must be replaced by statements which respect the full symmetry structure of the compactified graphene sheet .
finally , we demonstrate that these symmetries even extend to situations in which the carbon nanotori enclose different types of magnetic flux .
they thus should have applicability for many of the fascinating magnetic properties of carbon nanotori , including the possibility of persistent currents .
in order to explain the origins of these spectral symmetries , we begin with a brief review which will also serve to highlight our notation and conventions .      3.0 truein    for the uncompactified graphene sheet , all electron wavevectors @xmath9 are allowed .
however , this situation changes when we `` roll up '' one dimension of the graphene sheet to produce a single - walled carbon nanotube , or equivalently when we identify any two carbon atoms on the graphene sheet whose positions differ by an arbitrary lattice vector @xmath14 with @xmath15 . imposing the bloch condition on the electron wavefunctions @xmath16 in addition to the new periodicity condition
@xmath17 then restricts @xmath18 to the set of wavevectors satisfying the condition @xmath19 where @xmath20 and @xmath21 .
these allowed values of @xmath22 therefore form parallel lines in the @xmath23 plane , and the locations at which these lines intersect the bragg planes in fig .
[ fig : energycontours ] determine whether the corresponding @xmath24 nanotube is metallic , semiconducting , or insulating .
nanotubes for which @xmath25 ( , @xmath26 ) or @xmath27 ( , @xmath28 ) are dubbed `` zigzag '' or `` armchair '' respectively ; nanotubes with other values of @xmath24 are generally referred to as `` chiral '' . in general , the vector @xmath29 perpendicular to @xmath30 is the tube axis .    2.15 truein    we now consider imposing _ two _ non - parallel identifications on the graphene sheet , as illustrated in fig .
[ fig : torusparam ] . in general
, we consider two arbitrary lattice identification vectors @xmath31 and @xmath32 ; without loss of generality we shall assume that @xmath33 .
we shall also assume without loss of generality that the quantity @xmath34 is positive , or equivalently that the relative angle @xmath35 from @xmath30 to @xmath36 lies in the range @xmath37 , as illustrated in fig .
[ fig : torusparam ] . such identifications together result in a so - called @xmath0 nanotorus which may be viewed as a @xmath24 nanotube in which opposite ends are joined to each other , potentially with a relative twist angle @xmath35 .
this description is especially appropriate if @xmath38 .
note that @xmath39 describes the number of hexagonal cells which tile the surface of the resulting donut , and is equal to precisely half the number of carbon atoms in the nanotorus . imposing the double periodicity condition
@xmath40 on the electron wavefunction @xmath16 in conjunction with the bloch condition then leads to the constraint in eq .
( [ eq : tubek ] ) as well as the additional constraint @xmath41 where @xmath42 . solving eqs .
( [ eq : tubek ] ) and ( [ eq : tubek2 ] ) simultaneously then yields @xmath43 \nonumber\\    k_y & = & \frac{2\pi}{l_1l_2\sin\theta }      \left[-\ell_1 l_2\cos\delta + \ell_2 l_1\cos\beta\right ]   \label{eq : torik}\end{aligned}\ ] ] or equivalently @xmath44\nonumber\\    k_y & = &   { 2\pi \over \sqrt{3 } n_{\rm hex } } \left [ -\ell_1 ( 2p + q ) + \ell_2 ( 2m+n ) \right]~. \label{eq : torik2}\end{aligned}\ ] ] thus , we see that the allowed electron wavevectors now form a two - dimensional grid of points in the @xmath23 plane .
note that there are always exactly @xmath39 allowed wavevectors lying inside the fundamental region bounded by the bragg planes in fig .
[ fig : energycontours ] , where points lying on a single bragg plane itself are counted as half and where points lying at the intersections of two bragg planes are counted as a third .    whether such nanotori are metallic , semiconducting , or insulating depends on whether these points hit or come particularly close to the intersections of the bragg planes in fig .
[ fig : energycontours ] .
it is found that nanotori in which both @xmath45 and @xmath46 are multiples of three are metallic , with solutions to eqs .
( [ eq : torik ] ) and ( [ eq : torik2 ] ) that lie precisely on the fermi surface . for practical applications
, it proves useful to focus on only those nanotori with @xmath47 , as this condition allows one to bend the nanotube into a nanotorus without excessive strain or deformation of the underlying graphene sheet ( the effects of which would otherwise distort the dispersion relation in fig .
[ fig : energycontours ] , and hence the band structure of the torus ) . within this limit , it is then conventional to regard nanotori with @xmath48  ( mod  3 ) but @xmath49  ( mod  3 ) as semiconducting , since the allowed wavevectors come `` close '' to the fermi surface in this limit .
all other nanotori are then considered insulating .
all of the above results are completely standard , and are well known in the carbon - nanotorus literature .
in particular , tori with different values of @xmath0 are physically distinct : they have entirely different arrangements of hexagons tiling their surfaces , with different values of the physical radii @xmath50 , chiral angle @xmath51 , and twist angle @xmath35
. we shall take this to be our definition of `` physically distinct '' .
there are , of course , certain trivial identifications which relate different nanotori to each other : for example , the @xmath0 torus and the @xmath52 torus are actually identical , since they correspond to @xmath53 pairs which are related to each other by a uniform @xmath1 rotation .
such trivial identifications reflect the underlying hexagonal lattice symmetries of the graphene sheet from which these nanotori are constructed , and result in identical carbon nanotori with identical patterns of carbon atoms on their surfaces .
such nanotori are therefore not physically distinct .    by contrast , an important question is whether there exist physically _ distinct _ nanotori ( , tori which are _ not _ related by symmetries of the hexagonal lattice ) which nevertheless yield identical grids of allowed wavevectors @xmath23 . if so , we will have found cases of physically distinct tori with different values of @xmath54 which are nevertheless spectrally identical . in other words ,
such tori will have identical energy spectra and electrical conducting properties .    at first sight
, it might appear that no such spectral equivalences exist for carbon nanotori .
after all , the results in eq .
( [ eq : torik ] ) for carbon nanotori may initially appear to be nothing more than a two - dimensional generalization of the results in eq .
( [ eq : tubek ] ) for carbon nanotubes , and no such spectral equivalences exist for carbon nanotubes .
however , it is important to realize that this is not true for carbon nanotori . in particular
, it turns out that the set of allowed values of @xmath23 in eq .
( [ eq : torik2 ] ) is actually invariant under two additional symmetry transformations which we shall denote @xmath55 and @xmath56 : @xmath57 under these transformations , it is straightforward to demonstrate that @xmath34 is invariant , while the physical parameters @xmath58 , @xmath59 , @xmath35 , and @xmath51 transform in the following manner : @xmath60 it is important to stress that the individual equations in eqs .
( [ eq : torik ] ) and ( [ eq : torik2 ] ) are not invariant under @xmath55 or @xmath56 ; rather , what is invariant is the _ set _ of values of @xmath23 to which these equations lead .
moreover , since these sets of solutions for @xmath23 are invariant under @xmath55 and @xmath56 individually , they are also invariant under any sequence of @xmath55 and @xmath56 transformations . for example ,
under the @xmath61 transformation we find @xmath62 , and this too is a symmetry of the solutions to eq .
( [ eq : torik2 ] ) .
it then follows from these observations that any two tori whose defining vectors @xmath30 and @xmath36 differ through @xmath55 and @xmath56 transformations share the same electronic spectra but are intrinsically different from each other  , that they are isospectral but physically distinct .
that they are isospectral follows from the fact that the same set of @xmath23 solutions are selected in each case .
by contrast , that they are physically distinct follows from the fact that the fundamental identifications between carbon atoms on the graphene sheet are altered by @xmath55 and @xmath56 transformations in a manner that transcends trivial hexagonal lattice symmetries .
this is perhaps easiest to see in the case of the @xmath56 transformation , which corresponds to the action @xmath63 .
this changes not only @xmath59 but @xmath35 , and thus produces a new torus which has a greater `` twist '' when the ends of the nanotube are joined
.    it may be less obvious that the @xmath55 transformation also connects physically distinct tori .
indeed , this transformation corresponds to the action @xmath64 , and at first glance it might appear that this is merely a trivial relabeling of the two independent periodicities , along with a reflection ( sign flip ) .
however , we must remember that the @xmath55 transformation also changes the associated @xmath35 and @xmath51 angles in non - trivial ways .
alternatively , we can also appreciate the non - trivial nature of @xmath55 by considering the combination t   s t^-1 s^-1  =   t s t [ tprimedef ] which corresponds to the action @xmath65 .
[ the second equality in eq .
( [ tprimedef ] ) follows as a result of the identity @xmath66 , or equivalently @xmath67 . ]
note that @xmath68 is in some sense `` dual '' to @xmath56 : each performs a full twist around a different cycle of the torus .
in other words , as illustrated in fig .
[ tori ] , while @xmath56 corresponds to cutting through one side of the torus and re - attaching the two edges along with a twist , @xmath68 corresponds to cutting the torus along the other cycle and then twisting along the cut in an analogous manner before attaching the newly adjacent bonds .
thus , @xmath68 results in a physically distinct torus , just as @xmath56 does , and furthermore these resulting tori are not related to each other through hexagonal lattice symmetries .
it then follows from eq .
( [ tprimedef ] ) that @xmath55 also can not correspond to a lattice symmetry  , @xmath55 must connect physically distinct tori as well .
indeed , @xmath55 and @xmath68 are interchangeable in the sense that the two generators of the modular group can be considered to be either @xmath69 or @xmath70
.    3.5 truein    note that if we define the complex quantity @xmath71 , then the @xmath55 and @xmath56 transformations correspond to @xmath72 and @xmath73 respectively .
together , these transformations generate the so - called `` modular group '' , which is one of the primordial symmetries associated with toroidal compactifications . in general , any transformation which takes the form @xmath74 where @xmath75 and @xmath76 is a modular transformation .
what we have shown , then , is that any two carbon nanotori whose defining parameters are related by a modular transformation are spectrally identical even though they are physically distinct . in other words , any transformation which can be generated through repeated actions of the @xmath55 and @xmath56 generators leads to a spectral equivalence between physically distinct tori .    it is important to emphasize that the appearance of these modular symmetries is not a total surprise .
modular symmetries often arise in the presence of toroidal compactifications , and are a direct consequence of the underlying geometry of the carbon nanotorus .
indeed , the zone - folding of a macroscopic lattice necessarily introduces periodic equivalences of the sort discussed here , and in the case of two dimensions with non - parallel identifications , such periodicities generically lead to modular symmetries . in other words ,
the underlying toroidal graph is not altered by modular symmetries , and thus the adjacencies of the graph ( and the isospectrality which results ) are preserved intact . however , when applied to the case of carbon nanotori , these symmetries are _ physical _ , in the sense that the nanotori they relate are isospectral but physically distinct .
as we shall see , this is particularly significant , allowing these symmetries to have profound effects on the physics of these nanotorus systems and in particular on the range of behaviors they may exhibit .    in a sense
, these modular symmetries can be viewed as a generalization of a similar set of symmetries that emerge for mbius graphs . to see this ,
let us imagine taking a long strip of hexagons and gluing opposite ends of this strip together with an arbitrary number of half - twists . at first glance
, one might suspect that the resulting spectrum should depend on the actual number of half - twists .
however , it turns out that the resulting spectrum depends only on whether this number is even or odd  @xcite .
this mbius case is , of course , a relatively trivial example of this phenomenon , since it deals with only a single possible kind of twist corresponding to the single periodicity that defines the mbius strip .
indeed , for the case of carbon nanotori , our point is that there exists a whole _
group _ of transformations  the modular group  which is built upon the actions of _ two _ non - commuting underlying generators ( @xmath55 and @xmath56 ) which can be combined and performed in a variety of inequivalent ways .
this leads to `` sewing '' and
`` gluing '' configurations on the torus which are much more complex than on the mbius strip , and which require greater mathematical machinery to study  and this in turn leads to distinctions between the notions of `` physically distinct '' and `` spectrally distinct '' which are even more pronounced than in the mbius case .
however , even the relatively simple mbius case illustrates one of our basic themes : there can be lots of physically distinct twisted hexagon strips , yet very few of them are spectrally distinct .
it is the preservation of the underlying graph in each case which makes this possible .    as a result of these modular symmetries
, we can henceforth establish an unambiguous convention for uniquely describing the spectral properties of a given carbon nanotorus : we calculate its complex parameter @xmath77 as defined above , and then use a sequence of @xmath55 and @xmath56 modular transformations as needed in order to bring @xmath77 into a special region of the complex @xmath77-plane known as the `` fundamental domain '' @xmath78 defined by  :  -0  ,   ||1where @xmath79 and @xmath80 .
note that the condition @xmath81 already ensures that @xmath82 .
we can then use the underlying symmetries of the hexagonal graphene lattice in order to bring the resulting angle @xmath51 into the range @xmath83 . note that this convention is tantamount to describing the spectrum of a given carbon nanotorus in terms of that spectrally equivalent nanotorus which is as close to being rectangular as possible . following this procedure
therefore provides an unambiguous test of whether any two physically distinct carbon nanotori in fact have the same spectral properties .
the existence of these modular symmetries and the spectral equivalences they induce has profound implications for the physics of carbon nanotori . indeed , as we shall see
, the use of these modular symmetries will sharpen our ability to classify nanotori on the basis of those physical properties ( such as their metallicities , _ etc .
_ ) which ultimately stem from their band structures . furthermore
, use of these modular symmetries can even provide _ statistical _ insights into the properties of randomly produced carbon nanotori  insights which may ultimately have practical consequences for the manufacture of such objects in environments in which their resulting properties often can not be controlled or engineered in advance .    as discussed above ,
the first and most direct implication of these modular symmetries is that there exist physically distinct nanotori  often with markedly different radii and chiral angles  which nevertheless possess identical energy spectra and conductivity properties . as an example of this phenomenon ,
let us consider the @xmath84 , @xmath85 , and @xmath86 carbon nanotori .
clearly , we see that these nanotori are physically distinct and are characterized by different sets of cycle lengths @xmath50 , different chiral angles @xmath51 , and different twist angles @xmath35 . as a result , these three tori have entirely different patterns of carbon atoms tiling their surfaces . yet , the latter two nanotori are in fact related to each another by @xmath55 and @xmath56 transformations , and consequently they each yield the same spectrum of allowed @xmath22 vectors shown in fig .  [
fig : energyspecthreetori ] . similarly the first nanotorus is related to the second two through not only @xmath55 and @xmath56 transformations , but also trivial 60@xmath87 rotations of the underlying graphene sheet .
thus , these tori all share identical spectra , and according to the convention specified above , we can describe this spectrum uniquely as having @xmath88 , @xmath89 , and @xmath90 .    3.0 truein    second
, as a corollary to this observation , we also learn that the usual notions of `` zigzag '' and `` armchair ''  notions which are critical for describing the chirality and metallicity of carbon nanotubes  no longer hold as special indicators of the spectral properties of corresponding carbon nanotori .
in other words , nanotori built from zigzag or armchair nanotubes might have spectral properties which are identical to those of nanotori built from non - zigzag or non - armchair nanotubes ; moreover , these spectral properties might or might not correspond to the special zigzag or armchair angles @xmath91 . as an example , the @xmath92 nanotorus is built from a zigzag nanotube while the @xmath93 nanotorus is built from an armchair nanotube and the @xmath94 and @xmath95 nanotori are built from chiral nanotubes with different chiralities . yet all four nanotori have identical spectral properties which are the same as those of a chiral nanotorus with @xmath96 , @xmath97 , and @xmath98 ( or @xmath99 ) .
indeed , using the mathematical results of ref .
@xcite , we can show that any carbon nanotorus is spectrally equivalent to one built from a zigzag carbon nanotube .
conversely , a carbon nanotorus can have a spectrum which exhibits a zigzag or armchair property that is lacking in the original nanotube from which it is constructed .
as an example , the @xmath100 nanotorus is built from a zigzag nanotube and the @xmath101 nanotorus is built from a chiral nanotube . yet both nanotori have the spectral properties of an armchair nanotorus , with @xmath96 , @xmath102 , and @xmath103 .
likewise , the @xmath104 nanotorus is built from an armchair nanotube and the @xmath105 nanotorus is built from a chiral nanotube . yet both nanotori have the spectral properties of a zigzag nanotorus , with @xmath96 , @xmath106 , and @xmath107 .
third , it turns out that not every possible spectral signature @xmath108 can be realized , even in principle .
instead , these three quantities experience internal constraints and correlations which ultimately reflect the fixed hexagonal lattice structure of the underlying graphene sheet and which require that any allowed spectral signature @xmath108 originate from four real integers @xmath0
. these correlations amongst @xmath108 can take a variety of different forms .
for example , for odd values of @xmath39 , it turns out that there exist no self - consistent spectral signatures with purely imaginary values of @xmath77 .
similarly , for @xmath109 , there are solutions with @xmath110 , but no other inverse integers ; likewise , all of these except for @xmath111 have @xmath26 .
in fact , these internal constraints amongst the four real parameters embodied in @xmath112 are sufficiently strong that they effectively eliminate one real degree of freedom from within this parametrization .
in other words , although four real degrees of freedom are required in order to parametrize a given physically distinct nanotorus , only three real degrees of freedom are required in order to uniquely describe its spectral properties  @xcite .
it is easy to understand why such constraints arise .
ordinarily , as a question of topology , tori can exist with all shapes and volumes , for there are literally an infinite number of ways in which we can construct a torus by rolling up an unmarked sheet of paper .
however , in the case of carbon nanotori , our original `` sheet of paper '' is not unmarked : it is actually a graphene sheet of carbon atoms which has its own hexagonal lattice structure .
the existence of such a lattice structure has a number of critical consequences : it forces our toroidal defining vectors @xmath30 to @xmath36 to be _ lattice vectors _ ; it restricts the resulting possible combinations of @xmath39 and @xmath77 to those values consistent with the periodicity of the lattice ; and it breaks the rotational symmetry of our original uncompactified two - dimensional sheet and necessitates the introduction of a new measurable parameter , the angle @xmath51 defined in fig .  [ fig : torusparam ] , which describes the orientation of the torus relative to the underlying lattice .
conversely , the presence of the hexagonal lattice gives a clear meaning to the twist angle @xmath35 .
this quantity would have had no meaning when rolling up an unmarked sheet of paper .
fourth , it is clear that the existence of spectral equivalences between different nanotori implies that the number of _ spectrally _ distinct nanotori in any set will necessarily be smaller than the number of _ physically _ distinct nanotori in that set .
however , it turns out that the magnitude of this truncation can easily become quite staggering . as an example , let us consider the set of physically distinct @xmath0 nanotori with fixed @xmath113 which can be formed from integers in the range @xmath114 for some cutoff @xmath115 . for concreteness , we shall take @xmath116 . in order to count only physically distinct nanotori , we shall require that @xmath24 be chosen such that @xmath117 .
we shall also require , as a rough measure of their physical consistency in three - dimensional space , that each such nanotorus have an outer circumference @xmath59 which is at least triple its inner cross - sectional circumference @xmath58 ; note that other similar mathematical conditions may alternatively be imposed , but the qualitative results to follow are essentially unchanged .
given these constraints , we then find through direct enumeration that there are exactly @xmath118 physically distinct carbon nanotori which have @xmath88 . yet , as a result of these spectral equivalences , it turns out that these @xmath118 physically distinct carbon nanotori give rise to only @xmath119 distinct energy spectra !
these distinct energy spectra are listed in table  [ nhex18table ] , along with a representative carbon nanotorus in each class .
.for @xmath88 and @xmath116 , there are @xmath118 physically distinct carbon nanotori @xmath0 with @xmath120 . however , these exhibit only 14 spectrally distinct energy spectra and band structures , and only four of these correspond to metals .
these 14 spectrally distinct values of @xmath121 and @xmath51 are listed above , along with a sample @xmath0 nanotorus in each class . [ cols="^,^,^,^,^",options="header " , ]     this is clearly a major truncation .
even with the relatively small value @xmath88 and relatively small cutoff @xmath116 , each spectral signature listed in table  [ nhex18table ] is shared by literally hundreds or thousands of physically distinct carbon nanotori .
fifth , it is also worth noting that this truncation does not treat metallic and non - metallic nanotori equally . in general , a given @xmath0 carbon nanotorus will be metallic if at least one of its allowed wavevectors @xmath9 lies on the fermi surface , or in this case on one of the six points at which two bragg planes intersect .
as is well known , this occurs only when the differences @xmath45 and @xmath46 are each a multiple of three .
this implies that in any large set of physically distinct carbon nanotori , approximately one - ninth of the nanotori should be metallic .
however , we see from table  [ nhex18table ] that four out of the fourteen possible distinct spectral signatures are metallic .
this is almost triple what would have been expected , implying that the fraction of spectrally distinct nanotori which are metallic is nearly triple the fraction of physically distinct nanotori which are metallic . in other words ,
the truncation from the four - parameter space of physically distinct nanotori to the three - parameter space of spectrally distinct nanotori is remarkably sensitive to the metallicity of the nanotori in question .
one might argue that both of these effects  the huge truncation in the number of distinct spectral signatures and the relative abundance of those which are metallic  merely reflect the fact that we restricted our integers @xmath0 to lie within a fixed range bounded by @xmath122 , or that we took a relatively small value of @xmath39 . however , it is easy to demonstrate that the second of these effects is relatively insensitive to the choice of @xmath115 , and that the first of these effects only becomes even more dramatic as @xmath115 is increased .
for example , if we restrict our attention to carbon nanotori with @xmath96 , we find that the number of physically distinct nanotori and the number of spectrally distinct nanotori both rise as a function of @xmath115 . however , we see from fig .  [ panels ] ( left panel ) that the number of physically distinct nanotori with fixed @xmath39 grows as @xmath123 , as expected , while the number of spectrally distinct nanotori quickly hits a plateau ( right panel ) which remains flat for an increasingly long interval in @xmath115 before a new , hitherto - unrealizable spectral signature becomes possible and a new plateau develops .
[ this growth in the number of physically distinct nanotori as a function of @xmath115 can easily be deduced from the observation that the number of quadruplets of integers @xmath0 grows as @xmath124 , while the corresponding number of attainable values of @xmath125 grows as @xmath123 . ] as a result , the number of physically distinct nanotori quickly outpaces the number of spectrally distinct nanotori .
moreover , we see from fig .
[ panels ] that the number of distinct metallic spectral signatures remains roughly one third ( and not one ninth ) the total number of distinct spectral signatures .
we may also consider how these results vary with the choice of @xmath39 .
this is shown in fig .
[ nhexplot ] where , as functions of @xmath39 and for @xmath116 , we have plotted the number of spectrally distinct carbon nanotori as well as the number of those spectrally distinct carbon nanotori which are metallic .
it is clear from these results that the numbers of spectrally distinct nanotori remain relatively small , even though they rise with @xmath39 , as expected .
3.4 truein 0.2truein 3.4 truein    -0.1 truein    3.25 truein    it is also easy to understand the jagged , oscillating nature of the results in fig .  [ nhexplot ] .
nanotori must have @xmath126 in order to be metallic , while they must have @xmath127 if they are rectangular , with @xmath128 .
thus , one difference between tori with @xmath129 and those with @xmath130 is the existence of additional non - rectangular carbon nanotori with @xmath131 .
however , it must also be borne in mind that a random selection of four integers @xmath0 is @xmath132 times more likely to result in an even value for @xmath133 than an odd one . as we see from fig .
[ nhexplot ] , the combined effect from these two features is fairly significant .    finally , another important implication of these spectral equivalences between physically distinct carbon nanotori concerns the traditional rules of thumb which allow us to determine whether a given carbon nanotorus is metallic , semiconducting , or insulating . throughout the existing literature on this topic
, one finds what we shall call the `` rule of three '' : a given @xmath0 carbon nanotorus with @xmath134 will be metallic if both @xmath45 and @xmath46 are multiples of three , semiconducting if @xmath45 is multiple of three while @xmath46 is not , and insulating in all other cases . it is , of course , easy to verify that this characterization of a metallic nanotorus is modular invariant . specifically , if a given nanotorus has @xmath0 parameters which are metallic according to the rule of three , then all other nanotori to which it is spectrally equivalent will also have parameters which are metallic according to the rule of three .
in other words , if @xmath45 and @xmath46 are both multiples of three , modular transformations of these parameters will not disturb this property .
by contrast , the rule - of - three definition of semiconductors is _ not _ modular invariant  even when we preserve the condition that @xmath135 . as a graphic illustration of this point , consider the @xmath136 torus .
note that indeed @xmath135 for this torus . according to the standard rule of three
, such a torus can be identified as a semi - conductor because @xmath45 is a multiple of three . by contrast , let us now consider the @xmath137 nanotorus .
this torus also clearly has @xmath135 .
however , because @xmath45 is not a multiple of three , we would expect this nanotorus to be an insulator .
indeed , the rule of three tells us to expect this even though our first impression might be that the second nanotorus has larger radii in both directions , and therefore might have energy levels which are more closely spaced .
however , even though these tori have very different physical parameters , it turns out that they are related through modular transformations and therefore have identical energy spectra .
they therefore also have identical metallicity properties .
this provides a graphic illustration that as a mathematical statement , the standard `` rule of three '' fails to characterize the conductivity properties of such nanotori because it is inconsistent with the modular transformations which reflect the additional symmetries of the compactified graphene sheet . in other words ,
the standard definition for a semiconducting carbon nanotorus fails to be modular invariant , and thus can not be complete as a description of the underlying conductivity properties of the nanotorus .
this is yet another consequence of the fact that modular invariance in this context is an actual physical symmetry relating the spectra of physically distinct nanostructures .
we close this section with an important comment . throughout this section
, our goal has been to illustrate various mathematical ramifications of the modular symmetries which govern the spectra of different carbon nanotori .
the specific examples we have provided throughout this section were therefore chosen for their mathematical simplicity as opposed to their phenomenological practicality .
for example , we restricted our attention in this section to nanotori with relatively small values of @xmath39 , while realistic carbon nantoroi can be expected to have @xmath138
. likewise , realistic carbon nanotori will generally be quite long and thin , with quantities such as @xmath139 exceeding @xmath58 by an order of magnitude or more .
however , all of the conclusions we have drawn in this section continue to hold even when more realistic tori are considered . as an example , let us restrict our attention to carbon nanotori with @xmath140 for which @xmath141 .
the latter condition guarantees that our nanotorus remains relatively long and thin even if there are multiple windings of the hexagonal carbon lattice around the tube axis of the torus .
taking @xmath142 , we find there are @xmath143 physically distinct nanotori which satisfy these conditions , but only @xmath144 of these are spectrally distinct .
furthermore , of these @xmath144 spectral equivalence classes , @xmath119 correspond to metals .
we see , then , that even for realistic nanotori , our modular symmetries continue to lead to large classes of spectrally equivalent nanotori and a relative overabundance of classes which are metallic . indeed , in the limit @xmath145 , it is straightforward to show that the number of physically distinct carbon nanotori in each spectral equivalence class also grows to infinity .
with an eye towards potential implications of modular symmetries for the magnetic phenomena associated with twisted carbon nanotori  @xcite , we now consider the introduction of magnetic fluxes . for full generality
, we consider the possibility of _ two _ distinct fluxes : one which travels all the way around ( and through ) the length of the nanotube which forms the torus ; and another , namely the usual aharonov - bohm flux , which pierces the plane of the nanotorus , coming up through the donut hole .
we shall denote these fluxes @xmath146 and @xmath147 respectively , as their associated vector potentials @xmath148 and @xmath149 lie parallel to the vectors @xmath30 and @xmath29 in fig .
[ fig : torusparam ] , respectively .
as is typical for such systems , we then find that we can incorporate the effects of these fluxes by keeping our previous band - structure energy function @xmath150 in eq .
( [ eq : bandstightbinding ] ) unchanged , and simply modifying our constraint equations for @xmath23 in eqs .
( [ eq : tubek ] ) and ( [ eq : torik ] ) so that the integers @xmath151 are shifted according to _
i   _ i + _ i_0   where @xmath152 is the flux quantum and where _ 2   _ t + _ 1 _ 1  .
note that it is the possible existence of a non - trivial twist angle @xmath35 which is responsible for the distinction between @xmath153 and @xmath147 .
we then find that the resulting system continues to exhibit a spectral equivalence under the modular transformations in eqs .
( [ eq : stonindices ] ) and ( [ eq : stonvars ] ) as long as we allow @xmath154 to remain invariant under the @xmath56 transformation and to mix with each other under the @xmath55 transformation : s :   . thus
, a carbon nanotorus parametrized by @xmath155 will have the same spectral properties as one parametrized by @xmath156 , where these two sets of parameters are related through the modular transformations discussed above .
in this paper , we have highlighted and investigated the implications of a geometric symmetry  modular invariance  which emerges upon the compactification of a graphene sheet to form a carbon nanotorus .
although not traditionally considered in the carbon - nanotorus literature , modular invariance plays a critical role in describing the spectral properties of these nanotori and leads to spectral equivalences between physically distinct nanotori .
as we have shown , this has profound implications for the classification of carbon nanotori , indicating that large numbers of seemingly unrelated nanotori are in fact completely identical in terms of their spectral properties . along the way
, we also showed that the traditional `` rule of three '' classification rubric is incomplete , as it is based on quantities which are do not respect these spectral equivalences .
we also found that the fraction of spectrally distinct carbon nanotori which are metals is approximately three times greater than would naively be expected on the basis of standard results in the literature .
finally , we also showed that these spectral equivalences can easily be extended to cases in which non - trivial magnetic fluxes are present .
the existence of these spectral symmetries also provides a deeper theoretical underpinning to certain results which already exist in the literature .
for example , it is well known that many carbon nanotori exhibit persistent currents in the presence of a non - zero magnetic flux @xmath157 . as functions of @xmath158 ,
these currents typically follow complicated `` sawtooth '' patterns which have a natural periodicity under shifts @xmath159 .
it has also separately been observed ( see , , ref .
@xcite ) that any such sawtooth pattern is preserved but shifted horizontally upon the introduction of a nanotorus twist in which @xmath160 , where @xmath161 is chosen such that @xmath162 is also a lattice vector . indeed ,
when @xmath163 , the magnitude of this horizontal shift exactly matches the periodicity of the sawtooth pattern and the net result is unchanged .    remarkably , this coincidence is now easy to understand from the point of view of modular transformations : when @xmath163 , the mapping from @xmath164 is nothing but the @xmath56 modular transformation , and as we have shown , modular transformations preserve the electrical properties of the torus , including its persistent currents .
we thus see that the periodicity of the sawtooth pattern for persistent currents under shifts @xmath159  a periodicity which can be understood on elementary grounds having nothing to do with modular transformations  can now also be interpreted as a spectral equivalence under the @xmath56 modular transformation .
moreover , we now see that this is merely the tip of the iceberg , and that these sorts of spectral equivalences actually have a richer structure and context that not only corresponds to the entire modular group but also transcends the specific example of persistent currents .
as we have discussed above , the traditional `` rule of three '' is not formulated in a modular - invariant way and is therefore incomplete as a description of the electronic properties of carbon nanotori .
however , the full implications of these modular symmetries are significantly broader than just the rule of three : no theoretical result concerning the electronic properties of carbon nanotori can be correct unless it respects these modular symmetries . in other words ,
no calculation of any electronic property of a given carbon nanotorus in terms of its fundamental defining parameters @xmath0 can be correct unless it yields a result which is invariant under the modular transformations in eq .
( [ eq : stonindices ] ) . in this
respect , modular invariance functions for carbon nanotori in much the same way as gauge invariance functions for electromagnetic systems : no theoretical result can be correct unless it can be phrased in terms of quantities which are invariant under the symmetry in question . of course , at a mathematical level , both modular invariance and gauge invariance rest on relatively simple algebraic identities .
however , they both provide powerful organizing principles , and have significant physical manifestations and implications for the symmetry structures of the systems in which they appear .
needless to say , several additional comments are in order .
first , it should be noted that in this paper we have focused on what might called the `` ideal '' nanotorus . in particular , we have not accounted for the fact that the actual physical construction of such a torus in three - dimensional space requires that we introduce both an intrinsic and extrinsic curvature onto our otherwise flat graphene sheet . in this sense , the construction of a carbon nanotorus is different from that of a carbon nanotube ( in which only the extrinsic curvature is non - vanishing ) .
the introduction of intrinsic curvature requires that we subject our underlying graphene sheet to considerable strain , deforming not only the positions of carbon atoms but also their relative spacings .
these effects have been addressed by a number of authors , using a variety of different techniques  @xcite .
that said , these spectral equivalences should continue to hold , even in the presence of such deformations .
there are several reasons for this .
first , in the limit @xmath165 , all effects due to these deformations will be suppressed .
however , this is precisely the limit in which nanotori can be constructed from purely hexagonal graphene sheets without the introduction of curvature - inducing pentagonal or heptagonal carbon rings .
second , it can be shown that even when such strain is present , maximum toroidal stability occurs when this strain is uniformly distributed along the nanotorus  @xcite , and this is precisely the situation in which the techniques of ref .
@xcite can be used in order to mathematically rewrite the effects of such strain as arising due to a fictitious magnetic flux .
as we have seen , the spectral equivalences we have found continue to exist even when such fluxes are introduced .
but most importantly , while _ any _ deformations of the underlying graphene sheet can be expected to have effects on the corresponding band - structure energy function @xmath8 shown in fig .
[ fig : energycontours ] , such deformations will not disturb the symmetries inherent in the constraint equations for @xmath166 derived in eqs .
( [ eq : tubek ] ) and ( [ eq : tubek2 ] ) .
indeed , these equations reflect nothing more than the effects of toroidal compactification , and their structure leads directly to the modular symmetries inherent in eqs .
( [ eq : torik ] ) and ( [ eq : torik2 ] ) . thus , since these spectral equivalences ultimately stem from these symmetry properties , we are assured that any two tori related by modular transformations will sample the same set of wavevectors @xmath9 .
such tori will therefore continue to be spectrally identical regardless of the function @xmath8 , provided the deformations to @xmath8 for the two tori are themselves identical .
moreover , as we have argued above , even in cases where these deformations are not identical , they can at most differ by terms which are suppressed by factors of the presumably large nanotori radii .
such differences can therefore be safely neglected .
a similar conclusion also holds for thermal effects .
it might seem , at first glance , that thermal effects could also destroy the spectral equivalences , since they too can have a dramatic effect on the band structure of the underlying graphene sheet  @xcite . however , as noted above , these spectral equivalences are a consequence of the geometric symmetries that arise upon _ compactifying _ this sheet ; they are largely independent of the symmetries of the sheet itself . indeed ,
these modular symmetries exist for _ any _ choices of identification vectors @xmath30 and @xmath36 , even if those vectors are altered by other effects . therefore , as long as the temperature is sufficiently low that the electron coherence length exceeds both the inner and outer toroidal circumferences , the spectral equivalences we have been discussing should remain intact up to terms suppressed by the large nanotori radii .    these comments notwithstanding
, it still remains true that not all such nanotori are equally likely to appear in nature .
in particular , those nanotori whose constructions implicitly involve large numbers of twists are likely to be rather difficult to construct or stablize , as the carbon atoms that constitute the underlying graphene sheet are likely to experience significant strain , leading to major deformations of the underlying carbon bond lengths and angles away from their ideal values .
this is especially relevant , given that our statistical discussions in sect .
iv implicitly assume that each of the relevant @xmath0 nanotori can appear with equal probability .    at a mathematical level
, it is difficult to draw a firm boundary between those carbon nanotori which are realizable in nature as _ bona - fide _ molecules and those which are not .
however , in this paper we have restricted ourselves to statistical examinations of only those nanotori which already obey certain critical constraints .
for example , we have limited ourselves to nanotori in which @xmath58 and @xmath59 are both significantly bigger than the carbon - carbon bond length @xmath7 .
thus each twist experienced for the whole nanotorus has only a minor effect on the bond lengths and angles corresponding to each individual carbon atom .
in addition , we have further limited ourselves to nanotori which also satisfy @xmath167 .
this ensures that each such nanotorus can be realized in three - dimensional space , without the unphysical self - overlapping that would arise if this condition were not met .
together , these conditions help to ensure that the strains induced by these modular transformations are not too severe .
needless to say , there are many implications of these spectral equivalences which we have not yet explored .
it would be interesting , for example , to consider the implications of these symmetries for the existence or absence of persistent currents as well as the existence or absence of colossal magnetic moments .
this work is currently in progress
. it would also be interesting to understand these modular transformations in terms of fictitious fluxes , using analogues of the techniques presented in ref .
@xcite .
first , as we have seen , modular invariance is the symmetry which underlies most of the results we have presented in this paper . at a mathematical level , modular invariance is normally just a relabelling symmetry in the sense that two sets of torus parameters related by a modular transformation normally correspond to the same physical torus .
this is certainly the case in string theory , and in most other situations in theoretical high - energy physics in which modular transformations have played a significant role ( see , , ref .
@xcite ) .
however , the case of carbon nanotori is quite different . here , modular transformations relate parameters corresponding to carbon nanotori which are physically distinct .
as we have discussed , this is because we are not merely rolling up an unmarked sheet of paper when we subject it to two non - parallel identifications ; we are rolling up a graphene sheet which already has a hexagonal carbon lattice imprinted on it .
viewed from this perspective , it is therefore somewhat remarkable that modular transformations continue to play a role , indicating when two distinct nanotori will have the same spectral properties . indeed , in the case of carbon nanotori , we see that modular invariance is thus promoted from a mere relabeling symmetry to something far deeper : _ modular invariance becomes an outright physical symmetry between physically distinct entities_. we know of no other physical situation in which modular invariance plays such a role .
second , it is also exciting at a practical level that physically distinct carbon nanotori can have identical energy spectra and electrical properties . since these nanotori are physically distinct , some are likely to be far more complicated to construct in the laboratory than others .
nevertheless , these spectral equivalences suggest that it may not be necessary to fabricate a very complex nanotorus ( or generate sizable magnetic fluxes ) in order to obtain a desired spectral property ; there are likely to be far simpler nanotori and/or fluxes which can perform the same function .
this could have significant implications for the production and use of such nano - materials .
this work was supported in part by the department of energy under grant  de - fg02 - 04er-41298 .
we are happy to thank p.  mceuen , d.  ralph , c.  stafford , and h.  tye for discussions , and we wish to acknowledge the hospitality of the kavli institute for theoretical physics ( kitp ) in santa barbara , california , where this paper was completed .
the opinions and conclusions expressed here are those of the authors , and do not represent either the department of energy or the national science foundation .
m.  terrones , w.  k.  hsu , j.  p.  hare , h.  w.  kroto , h.  terrones , and d.  r.  m.  walton , phil .
a * 354 * , 2025 ( 1996 ) ; h.  r.  she , r.  martel , and p.  avouris , phys .
* 84 * , 4441 ( 2000 ) ; m.  sano , a.  kamino , j.  okamura , and s.  shinkai , science  * 293 * , 1299 ( 2001 ) ; h.  watanabe , c.  manabe , t.  shigematsu , and m.  shimizu , appl .  phys .  lett .  * 78 * , 2928 ( 2001 ) .
m.  f.  lin , r.  b.  chen , and f.  l.  shyu , solid state commun .
107 , 227 ( 1998 ) ; a.  a.  odintsov , w.  smit , and h.  yoshioka , europhys . lett . * 45 * ( 5 ) , 598 ( 1999 ) ; k.  sasaki , phys .
b * 65 * , 155429 ( 2002 ) ; s.  latil , s.  roche , and a.  rubio , phys .  rev .
b * 67 * , 165420 ( 2003 ) ; z.  zhang , z.  yang , x.  wang , j.  yuan , h.  zhang , m.  qiu , and j.  peng , j.  phys .  condens .
matter * 17 * , 4111 ( 2005 ) ; c.
feng and k.  m.  liew , carbon * 47 * , 1664 ( 2009 ) .


various approaches have been used to understand ageing experiments @xcite @xcite @xcite : droplet picture @xcite , mean field model @xcite @xcite , and trap models @xcite @xcite @xcite @xcite .
ageing is by definition the property that the two points correlation function @xmath1 depends explicitly on the waiting time @xmath2 .
@xmath3 is a macroscopic observable , for instance the magnetization in the zero field cooled experiment @xcite or the thermoremanent magnetization experiment @xcite . on the theoretical level
, dynamics can be implemented in several ways .
one possibility is langevin dynamics which was used for instance in @xcite where mean field ageing was first obtained , or in @xcite where ageing in the presence of flat directions was studied .
another possibility is glauber dynamics @xcite . within the framework of glauber dynamics and
given the probability distribution @xmath4 , for the system to be in any microcanonical state @xmath5 at time @xmath6 , the probability distribution at time @xmath7 is given by @xmath8 where @xmath9 is the glauber matrix . by `` microcanonical states '' , we mean a state such as the system is in a given microcanonical configuration with probability one .
this article is devoted to the study of the glauber dynamics of various models , phenomenological or microscopic , in connection with ageing .
more precisely , for a finite size system of size @xmath10 , we calculate from the glauber matrix a time @xmath0 such as the finite size system equilibrates on a time of the order of @xmath0 , starting from out  of  equilibrium initial conditions at time @xmath6 . in terms of ageing , the correlation function ( [ eqcorr ] ) depends explicitely on the waiting time @xmath2 provided the system is out  of  equilibrium at time @xmath2 , namely @xmath2 is smaller than @xmath0 .
in other words , ageing for a finite size system is interupted after @xmath0 [ see @xcite where the denomination `` interrupted ageing '' was proposed ] . in order to understand ageing in the thermodynamic limit ( @xmath11 )
, the variations of @xmath0 with the system size have to be carefully studied .
if @xmath0 goes to a constant in the thermodynamic limit , ageing is interrupted , whereas if @xmath0 diverges ageing is not interrupted in the thermodynamic limit : in the thermodynamic limit , the system reaches equilibrium on infinite time scales .
we also discuss the possibility of metastability without ageing on experimental time scales , as it is the case for instance in diamond .    in general
, @xmath0 depends on the eigenvalues of the glauber matrix , but also on the overlaps between microcanonical states , and the transformation matrix elements between the eigenstates of the glauber matrix and the microcanonical states .
however , in some simple cases , with only self  overlaping states , the time @xmath0 can be expressed only in terms of the density of eigenvalues of the glauber matrix .
the form of the glauber matrix will be given in each case in the core of the paper .
whatever the system , the eigenvalues of @xmath9 are negative .
the spectrum is bounded above by zero ( zero is always an eigenvalue , corresponding to the equilibrium boltzmann distribution ) .
the eigenvalues of the glauber matrix have the dimension of the inverse of a time , so that they can be viewed as the inverse relaxation times of the system . since we are interested in long time behavior , a special attention will be paid to the smallest eigenvalues .
this article is organized as follows . in section [ gm ] , we present a detailed study of the average density of small eigenvalues of the random energy model ( rem ) like trap model proposed and studied in @xcite and inspired by the rem model @xcite .    in the case of the rem
 like trap model , two regimes are predicted and checked numerically for the density of small eigenvalues in the glass phase : an arrhenius  like regime and a non arrhenius regime for the smallest eigenvalues . by arrhenius  like regime , we mean a regime in which the density of eigenvalues is the same as the density of eigenvalues calculated from the arrhenius law , assuming a correspondence between the eigenvalues of the glauber matrix and the inverse relaxation times of the system .
the eigenvalue density in the arrhenius regime are given in section [ arrh - def ] .
we also establish a connection between the behavior of the density of small eigenvalues of rem  like trap models and the spectral dimension of the associated diffusion process .    next ,
in section [ general ] , we give a general relation for the existence of ageing from a glauber dynamics point of view .
this section is the main part of the present article . when applied to the rem  like trap model , our criterium is equivalent to the fact noticed in @xcite that the average trapping time diverges below @xmath12 .
the rest of the paper is devoted to examine the case of other models : brownian motion which was shown in @xcite to be an ageing phenomenon , a model introduced by barrat and mzard @xcite ( bm model ) , the one dimensional ising chain ( in this case , our results are consistent with the existence of interrupted ageing ) , and the sherrington  kirkpatrick model ( sk model ) . in all these cases ,
we find consistent results between the existence of ageing calculated from the behavior of the two times autocorrelation functions and our criterium based on the glauber matrix .
we end  up with some final remarks .
trap models can be defined without a priori reference to a microscopic model
. however , some trap models , such as the rem  like trap model @xcite are related to the parisi solution of spin glasses : even though these models are toy models , some of them retain some feature of the dynamics of microscopic models .
the rem  like trap model consists of a trap model in infinite dimension , that is , the system can `` hop '' from one trap to any other trap , according to the rules of the dynamics .
other toplogies will be analyzed below , in reduced dimension .
the rem  like trap model corresponds to a single replica symmetry breaking step .
generalizations of the rem  like trap model including multi replica symmetry breaking steps were studied in @xcite and their relevance to explain experimental data was shown .
the models of this section consist of @xmath10 traps of depth @xmath13 , with the @xmath14 chosen among the distribution @xmath15 .
the reason why this particular energy trap distribution is chosen is that there exists a finite transition temperature between a high and a low temperature glassy phase @xcite .
the dynamics rules are given by @xmath16 for @xmath17 if the traps @xmath18 and @xmath19 are connected and zero otherwise .
the diagonal terms are defined by @xmath20 which enforces probability conservation .
equation ( [ propcons ] ) is a general constraint for any glauber dynamics .
another constraint is the detailed balance condition @xmath21 it is clear that the aforementioned rules for the dynamics of the rem  like trap model satisfy detailed balance .
however , as we shall see below , other types of dynamics also satisfy detailed balance and probability conservation . as a consequence of detailed balance , the eigenvectors of the glauber matrix are not physical states ( in the sense that they are not probability distributions ) except for the boltzmann distribution ( the stationary eigenvector corresponding to the @xmath22 eigenvalue ) .
notice that if @xmath9 is a glauber matrix @xmath23 , with @xmath24 a positive number , is also a glauber matrix .
multiplying @xmath9 by @xmath24 amounts to changing the unit time .
the proportionality constant is imposed by physical considerations .
for the models considered in this section , this constant is chosen such as the diagonal elements of the glauber matrix do not scale with the number @xmath10 of traps in the large @xmath10 limit .
for instance , in the infinite dimensional case , we choose @xmath25 .
this dynamics will be studied with different phase space topologies : infinite dimension , where all the traps are connected ( section [ infinite ] ) , one dimensional topology , where a given trap is connected only to its two nearest neighbors ( section [ disor ] ) , and also one dimensional topology with `` ordered '' traps , in the sense that the trap energies have been ordered @xmath26 , and the trap at site @xmath18 has an energy @xmath14 .
a useful property that we will use throughout this paper is that glauber matrices can be mapped onto symmetric matrices with the same spectrum .
let us call @xmath27 the diagonal matrix @xmath28 .
then , because of detailed balance , the matrix @xmath29 is symmetric and has the same spectrum as @xmath9 . using a quantum mechanical analogy ,
we denote by @xmath30 the microcanonical state corresponding the system in trap @xmath18 with a probability unity .
the analogy to quantum mechanics is useful , even though there is no hermitian structure .
in particular , any linear combination of microcanonical states is not a probability distribution .
we propose ourselves to carry out a detailed study of the density of eigenvalues of the rem  like trap model and related models in reduced dimension . since the arrhenius law is introduced `` by hand '' in the dynamics , it is interesting to question in which regime the density of eigenvalues is given by the arrhenius law , assuming a correspondence @xmath31 between the eigenvalues of the glauber matrix and relaxation times
. we will naturally call _ arrhenius regime _ such a regime . as we will see below ,
the density of eigenvalues is not always in the arrhenius regime .    in this arrhenius regime , and by definition , the off - diagonal terms are neglectable as far as the density of eigenvalues is concerned .
the diagonal coefficients of the glauber matrix are proportional to @xmath32 , so that the distribution of eigenvalues is @xmath33 this distribution is the same as the one expected from an arrhenius law with life times @xmath34 .
if one assumes an arrhenius picture , the distribution of relaxation times is @xmath35 which is equivalent to ( [ aha ] ) .
as noticed in @xcite , the average trapping time diverges if @xmath36 , leading to ageing . distribution of trapping times with no first moment were already studied in @xcite , and shown to lead to anomalous diffusion , namely to a non trivial spectral dimension , even above the upper critical dimension .
we will compare our results for the rem  like trap model with the results of @xcite derived of the context of anomalous diffusion on site
disordered hypercubic lattices .
the density of eigenvalues is thus described by the arrhenius law provided the off - diagonal terms of the glauber matrix are a small perturbation , so that the density of eigenvalues of the glauber matrix is given by the diagonal terms only .
we now examine several cases with different phase space topologies and determine in the different cases on which condition the density of eigenvalues is determined by the arrhenius picture .
we also compare our analytical approach to numerical diagonalizations of the symmetrized matrix @xmath37 .      in this section ,
the system can hop from one trap to any other trap ( infinite dimension ) .
the diagonal coefficients are @xmath38 and the off  diagonal coefficients of the symmetrized glauber matrix are @xmath39      we treat the off
diagonal coefficients @xmath40 ( @xmath17 ) as a perturbation .
this perturbative treatment is consistent since we are looking for a regime in which the off  diagonal coefficients
do not modify the density of eigenvalues with respect to the purely diagonal case .
the second order correction to the eigenvalue @xmath41 is @xmath42 the first order corrections being identically zero since the perturbation consists only of off  diagonal elements .
we assume that the density of trap energies in the vicinity of @xmath14 is large enough for the discrete summation to be replaced by an integral .
if this condition is not satisfied , the second order corrections @xmath43 may have large fluctuations because the divergence for @xmath44 is no more cancelled .
the condition that the unperturbed eigenvalue density is large in the vicinity of @xmath14 simply reads @xmath45 .
we assume in the rest of this section that this condition is fulfilled , and we can safely replace the discrete sum by an integral : @xmath46 the singularities for @xmath47 and @xmath48 cancel each other , and we can approximate @xmath43 as @xmath49 with @xmath50 .
we have used the following approximation : if @xmath51 is small , then @xmath52 and we use @xmath53 as a cut  off .
clearly ( [ app ] ) is only an approximation , but we are only interested in the dominant behavior of @xmath54 as a function of @xmath14 .
now , we want to compare the approximation ( [ app ] ) of @xmath43 to the typical level spacing between @xmath41 and @xmath55 , where we assume that the traps are ordered such as @xmath56 .
to to so , we first calculate the distribution of @xmath57 and we obtain @xmath58 , with @xmath59 notice here that @xmath60 diverges as @xmath61 goes to zero if @xmath36 .
this accumulation of small relaxation times is related to the existence a glass transition for this model for @xmath62 .
let us now calculate the level spacing statistics @xmath63 , which is the density probability that the spacing between the absolute value of the eigenvalue @xmath9 and the absolute value of the next eigenvalue is @xmath64 .
because of the statistical independence of the trap depths , @xmath63 is readily obtained as @xmath65 in this expression , @xmath66 is the number of levels in the interval @xmath67 $ ] ( with @xmath68 ) .
the other terms correspond to the condition that there is no eigenvalue in the interval @xmath69 $ ] . assuming that @xmath61 is small compared to unity , and that @xmath10 is large , we get @xmath70 the off - diagonal terms are irrelevant to the density of eigenvalues provided @xmath71 is much smaller than unity .
this condition reads @xmath72 , with @xmath73 defined in ( [ deltai ] ) . since we are looking for a necessary condition valid
if @xmath36 , we neglect the temperature  dependent prefactor and conclude that the density of states is in the arrhenius regime only for traps such that @xmath74 .
notice that the traps such as @xmath75 are deeper than the typical deepest trap ( of the order of @xmath76 ) .
the arrhenius picture is thus valid for a _ typical _ trap realization , with a lowest trap not deeper than @xmath76 .
however , non typical disorder realizations may be generated .
for these realizations , the lowest eigenvalues are not described by an arrhenius picture .
the existence or absence of the arrhenius regime even for a very large number of traps is thus sample dependent .
this non self averaging behavior is reminiscent on the sample dependence of the quantity @xmath77 with @xmath78 the probability to find the system in trap @xmath18 @xcite .
the density of small eigenvalues of the rem  like trap model is pictured on figure [ fig1 ] . for practical purposes , we do not calculate the density @xmath79 since this quantity is too small for small eigenvalues ( see ( [ aha ] ) ) .
we rather calculate @xmath80 , where @xmath81 is the boltzmann weight associated to the eigenvalue @xmath82 , which is calculated as the average of the boltzmann weight operator @xmath83 over the eigenstate associated to the eigenvalue @xmath82 . in the arrhenius picture , we have @xmath84 since @xmath36 , @xmath85 increases exponentially as a function of @xmath86 when @xmath87 .
we see on figure [ fig1 ] the existence of the arrhenius regime , in good agreement with the condition @xmath74 .
we can now fit the density of eigenvalues in the arrhenius regime to the form @xmath88 and see if we recover @xmath89 , as predicted from the arrhenius law .
the fit to the arrhenius law is shown on figure [ fig2 ] .
the variations of the exponent @xmath90 as a function of the inverse temperature are plotted on figure [ fig3 ] , and we conclude to a good agreement with the arrhenius law .
random walks on a finite @xmath91dimensional lattice with a broad distribution of waiting times have been analyzed in the past ( see @xcite and references therein ) . as we shall see
, we can relate our results concerning the existence of a arrhenius  like regime for the rem  like trap model to the results of @xcite , where the long time behavior of random walk on such lattices was investigated by means of the real space renormalization method introduced in @xcite .
the rem  like trap model falls in the universality class described by machta since the ingredients are the same in the two cases , rather than the universality class of random walks in a random field force @xcite .
in order to make contact with the work of machta , we first calculate the probability @xmath92 of being at the origin : @xmath93 where @xmath94 is the probability to find the system in the trap @xmath18 at time @xmath7 , starting from the trap @xmath18 at time @xmath95 , and @xmath96 is the set of eigenvalues of the glauber matrix .
we have shown that , in the infinite dimensional case , the typical small eigenvalue density ( relevant to the long time behavior of @xmath97 ) has the arrhenius form @xmath98 .
a simple scaling analyzis then shows that @xmath99 in the long time limit .
this result was derived for a graph such as any site is connected with any other site , that is we started from the begining from an infinite dimensional situation . in order to extrapolate this result to a finite dimension @xmath91 larger than the upper critical dimension 2 ( as shown in @xcite ) we first notice that , within the framework of the rem  like trap model , the microscopic time @xmath100 scales like @xmath101 , with @xmath102 the lattice spacing , whereas , in the context of random walks , the correct scaling is @xmath103 [ an illustration of the validity of such a scaling will be given in section [ brownian - motion ] ] .
this implies that @xmath104 , which gives the correct prescription to go from the rem  like trap model in infinite dimension to the finite @xmath91 dimensional random walk model : one should replace the time @xmath7 in ( [ p0rem ] ) by @xmath105 .
this leads to @xmath106 for @xmath91 larger than the upper critical dimension @xmath107 .
this is exactly what was found in @xcite using a real space renormalization group approach , with a non trivial spectral dimension @xmath108 .
we aim to calculate the structure of the eigenstates of the rem  like trap model .
to do so , we decompose a given microcanonical state @xmath30 into the eigenvectors @xmath109 , and calculate the absolute value of the coefficients of this decomposition .
the result is plotted on figure [ fig4 ] .
we see that a microcanonical state @xmath30 with an eigenvalue @xmath110 is mainly a linear combination of eigenstates with eigenvalues close to @xmath110 .
this observation can be understood in the framework of perturbation theory as follows .
the first order correction to the eigenvectors are @xmath111 so that , at first order @xmath112 at the first order of perturbation theory , the pure states @xmath30 is thus mainly a linear combination of eigenvectors with an eigenvalue close to the eigenvalue of @xmath113 .    starting from a pure state @xmath30 at time @xmath6 , we can calculate the evolution of the system : @xmath114 \nonumber , \end{aligned}\ ] ] where we have used ( [ per1 ] ) and ( [ per2 ] ) .
as expected from the arrhenius law , the relaxation time is of the order of @xmath115 , which diverges in the limit @xmath116 .
we will compare below the structure of the eigenvectors of the rem  like trap model to the structure of eigenvectors of another model ( bm model ) .
assuming periodic boundary conditions , the diagonal terms of the glauber matrix are given by @xmath117 the perturbation induced by the off - diagonal terms does not affect significatively the density of eigenvalues provided the perturbation @xmath118 is small compared to the level spacing between @xmath41 and @xmath55 .
since @xmath119 is uncorrelated from @xmath14 , the matrix element @xmath120 must be averaged over @xmath119 , so that @xmath40 is of the order of @xmath121 , that is @xmath122 the distribution of eigenvalues of the unperturbed system is @xmath123 from what we deduce the level spacing statistics in the long time ( small @xmath61 ) approximation @xmath124 the off - diagonal terms do not affect the density of eigenvalues provided @xmath125 which leads to the condition @xmath126 .
the average number of traps generated among a sequence of @xmath10 independent traps chosen among the distribution @xmath127 and fulfilling the aforementioned condition is @xmath128 which is vanishing in the large @xmath10 limit . we thus conclude that the arrhenius regime with this topology does not exist , that is off - diagonal terms will always modify significatively the density of eigenvalues with respect to the arrhenius case .
as we will see below , the arrhenius regime is replaced by a renormalized arrhenius regime , where the density of eigenvalues is given by ( [ lawarr ] ) , but with @xmath90 lower than @xmath129 .
again , we generate numerically the density of small eigenvalues @xmath85 .
the fit is shown on figure [ fig2 ] , and the variations of the exponent in the arrhenius regime are plotted on figure [ fig3 ] .
we conclude to the existence of a renormalized arrhenius regime , in the sense that the density of eigenvalues is given by ( [ lawarr ] ) but with an exponent @xmath90 smaller than @xmath129 , indicating a slowing down of the dynamics due to the one dimensional topology compared to the infinite dimensional case : in one dimension , to go from one trap to a deeper trap , the system has to jump to intermediate traps .
section [ relation1 ] was devoted to the connection between the spectrum of the glauber matrix in infinite dimension and the long time behavior of random walks on finite @xmath91 dimensional trapped hypercubic lattices above the upper critical dimension .
we now analyze the analogy in one dimension .
we have shown numerically that the density of eigenvalues behaves like @xmath130 for small @xmath82 , where @xmath131 is different from @xmath9 . in terms of random walks , this result becomes , with the notations of section [ relation1 ] @xmath132 namely the spectral dimension is @xmath133 .
we now use the results of @xcite where @xmath134 was calculated : @xmath135 from what we conclude that the coefficient @xmath90 plotted on figure [ fig3 ] is @xmath136 this analytic form is consistent with numerical calculations of the density of eigenvalues .
we now assume that @xmath10 energy traps were generated among the distribution @xmath127 , and that these energies have been ordered in such a way that @xmath56 .
a site @xmath66 with a trap of energy @xmath137 is connected to the sites @xmath138 ( with an energy @xmath139 ) and @xmath140 ( with an energy @xmath141 , except for the sites @xmath142 and @xmath10 which have no neighbor .
intuitively , such a trap ordering should make the dynamics faster with respect to the case analyzed in section [ disor ] since the system can decrease its energy continuously during the relaxation .
we have @xmath143 since the density of level spacings @xmath63 ( [ p(s ) ] ) is large for small values of @xmath61 , we expect that @xmath55 is close to @xmath41 at least if @xmath144 is small , so that we can use the following approximation : @xmath145 .
the off - diagonal terms will not modify the density of eigenvalues with respect to the arrhenius case provided @xmath146 that is @xmath75 .
the number of traps fulfilling this condition in a distribution of @xmath10 traps is on average equal to unity , so that only a small number of eigenvalues contribute to the arrhenius regime .
nevertheless , it is still possible to make a statistics over a large ensemble of @xmath10 traps distribution and to test numerically if the arrhenius regime can be observed .
we generated numerically the density of small eigenvalues @xmath85 .
as shown on figure [ fig2 ] , we still can fit @xmath147 to the form @xmath148 , and the exponent @xmath90 is close to @xmath129 , as in the arrhenius regime .
( see figure [ fig3 ] ) .
we conclude that ordering the traps in one dimension restores the arrhenius regime for the lowest eigenvalues .
again , the existence of the arrhenius regime is sample dependent .
we are now going to give a general relation between ageing and the glauber matrix .
this relation is a generalization of a criterium given in @xcite . in the context of the rem
like trap model , the average trapping time @xmath149 diverges if @xmath150 . on the other hand , ageing exists below @xmath12 only , in the limit where the microscopic time scale @xmath100 goes to zero .
we are now going to argue that a similar criterium can be found for any statistical mechanical model , provided the master equation is linear .
this is not always the case .
for instance , the master equation of the backgammon model in a mean field dynamics at zero temperature is not linear @xcite .
let us now start our general argument .
the idea is to start with a system in a given microcanonical state @xmath30 and write the glauber dynamics equation for the evolution of this pure state .
first , we call @xmath151 the transformation matrix @xmath152 where @xmath30 is a microcanonical state and @xmath153 are the eigenstates of the glauber matrix . by microcanonical state , we mean that if the system is in the state @xmath30 , the system is with probability one in the microcanonical state @xmath5 . in order to avoid confusion ,
all the labels related to the eigenvectors will be denoted by greek letters whereas the labels for pure states will be denoted by latin letters .
the evolution of the pure state ( [ pure2 ] ) is readily obtained as @xmath154 in other words , the probability to be in microcanonical state @xmath155 at time @xmath7 starting from the microcanonical state @xmath5 is @xmath156 we now consider the time
dependent correlation function @xmath157 where @xmath158 is the number of microcanonical states .
we use @xmath158 instead of @xmath10 to avoid confusion between the number of microstates @xmath158 and the number of sites @xmath10 in a spin model . in a trap model , @xmath159 , whereas in an ising spin system , @xmath160 . in ( [ correlation ] ) , we have carried out an average over all the initial microcanonical states , that is the system is first quenched from an infinite temperature , and the correlation ( [ correlation ] ) is averaged over the realizations of the infinite temperature states . in ( [ correlation ] ) , @xmath161 is the overlap between the pure states @xmath5 and @xmath155 , which has to be defined within each model . as we will see below in section
[ brownian - motion ] , even within a given dynamics , the choice of the overlap depends on what type of observables are considered in the two times correlation functions .
we now want to calculate from ( [ correlation ] ) a typical time by integrating ( [ correlation ] ) .
however , the boltzmann eigenstate leads to a trivial divergence of this integral .
we thus need a prescription to eliminate the equilibrium distribution in ( [ correlation ] ) .
we first discuss this prescription in the absence of metastability and postpone for later discussion the case of metastability .
we define the typical time @xmath162 by @xmath163 where the summation over the eigenstates excludes the boltzmann distribution .
the overline denotes an average over the disorder realizations , if necessary .
the significance of @xmath0 is clear : for a finite size system with @xmath10 sites , the out  of  equilibrium dynamics occurs before @xmath0 and , for times larger than @xmath0 , the finite size system has equilibriated . in other words
, @xmath0 is the cross  over time before which the finite size system exhibits ageing .
notice that in this argument , the system is a finite size system .
one should thus carefully study finite size effects .
two cases may occur : first , @xmath0 tends to a constant in the thermodynamic limit . in this case , ageing is interrupted @xcite in the thermodynamic limit .
as we shall see , this is the case for the ising chain at a finite temperature . on the other hand ,
if @xmath0 diverges in the thermodynamic limit , ageing is not interrupted . in this case , the exitence of a divergent relaxation time with respect to times @xmath95 and @xmath7 implies the existence of a divergent relaxation time with respect to any finite waiting time @xmath2 and @xmath164 , and thus the existence of ageing in the thermodynamic limit .
the relation between the existence / absence , the nature of ageing in the thermodynamic limit and the behavior of @xmath0 was derived on heuristic basis .
we could not find a rigourous derivation using a general argument .
however , in what follows , we investigate the behavior of several models ( rem  like trap model , random walks , bm model , one dimensional ising chain and sk model ) and we find in all these cases that our criterium based on the behavior of @xmath0 is correct .
as far as non interrupted ageing in the thermodynamic limit is concerned , we will see below that this is the case for the bm and sk models .
this is also the case for the rem  like trap model .
we have shown in section [ infinite ] that for all the typical realizations of the trap configurations ( typical is defined in section [ infinite ] ) , the density of eigenvalues is given by the arrhenius law , which means that , in the case of the rem  like trap model , our criterium for the existence of ageing is equivalent to the criterium of bouchaud @xcite that we have recalled at the beginning of this section , provided the overlap between traps is taken to be @xmath165 , where @xmath166 is the edwards  anderson order parameter ( within the one replica symmetry breaking step picture ) .
notice that the existence of a non arrhenius regime generated by non typical trap configurations does not modify our conclusion regarding the divergence of @xmath167 since the non arrhenius regime only occurs for traps such as @xmath168 if @xmath11 .
the correlation @xmath169 is @xmath170 which is non integrable if @xmath171 .
the typical time @xmath172 only depends on the density of eigenvalues @xmath173 for a system of @xmath10 traps , and is independent on the transition matrix , which motivates the study of the density of eigenvalues carried out in section [ gm ] . in general , in any model with only self  overlapping pure states , the existence of ageing is coded only in the behavior of the density of small eigenvalues : there is ageing in the thermodynamic limit if @xmath174 tends to a constant if @xmath82 goes to zero ( an exemple will be provided in section [ bm - model ] ) or if @xmath174 diverges if @xmath82 goes to zero ( this is the case for the rem  like trap model already studied in section [ gm ] ) .
let us now discuss the behavior of the correlation function ( [ correlation ] ) in the presence of metastability .
we have in mind the case of diamond for instance , which is metastable but in which case there is no ageing .
we do not have a glauber matrix description of diamond , but take this case as a generic situation for metastability and draw qualitative conclusions . in this case
we expect that for a finite size system , ( [ correlation ] ) will first decay to a finite value over a first time scale @xmath175 .
after @xmath175 , the system reaches its metastable configuration ( for instance diamond ) and thermalizes in a portion of its phase space .
we thus expect that , for times larger than @xmath175 , the correlation ( [ correlation ] ) reaches a plateau . on larger time scales , the correlation ( [ correlation ] ) should decay to its equilibrium value with a cross  over time @xmath176 , corresponding for instance to the transition from diamond to graphite .
the time scale @xmath176 should increase with the number of sites , so that , with a macroscopic number of sites , no transition from diamond to graphite is observable on experimental time scales . on the other hand ,
@xmath175 may remain small even in the thermodynamic limit , signaling the absence of ageing in the graphite phase .
in such a situation where the cross  over time scale @xmath176 is not observable , the only relevant time scale is @xmath175 . in order to calculate @xmath175 one should not use the prescription ( [ prescription ] ) but rather substract to the correlation ( [ correlation ] ) its plateau value for @xmath177 .
brownian motion was shown in @xcite to be an ageing phenomenon .
this model is simple and thus provides a simple test of the criterium of section [ general ] .
brownian motion in an euclidian space of dimension @xmath91 has the following glauber dynamics : @xmath178 and @xmath179 if @xmath18 and @xmath19 are nearest neighbors on a hypercubic lattice .
@xmath100 is the microscopic hopping time . in this case
, the glauber matrix is symmetric and , because of translational invariance , its eigenstates are plane waves @xmath180 where @xmath181 ( we have assumed cyclic boundary conditions with a lattice spacing @xmath102 ) .
the matrix elements of the glauber matrix on the plane wave basis are @xmath182 where @xmath90 runs over the @xmath91 directions and @xmath183 is the component of the wave vector @xmath184 on the @xmath90 direction . the long time behavior of the random walk is readily obtained by expanding the cosine up to the second order , and the probability to find the brownian particle at time @xmath2 on site @xmath185 is , in the long time limit : @xmath186 where the particle is at the origin at time @xmath6 ( out of equilibrium state ) . as noticed first in @xcite , the two times correlation functions of the position are proportional to the waiting time : @xmath187 in the large @xmath2 limit . in order to make contact with @xcite , one should adopt the scaling @xmath188 in the thermodynamic / long time limits , which allows to recover the langevin dynamics results and thus the existence of ageing in any dimension for the position  position correlation function .    in order to make contact between the correlation ( [ c(1 ) ] ) and the general form of the correlation ( [ correlation ] ) , we should choose the overlap @xmath189 . with this choice of the overlap , ( [ correlation ] ) and ( [ c(1 ) ] )
are both related to the same physical quantity : the mean square displacement . using this specific form for the overlap
, one can calculate @xmath190 , integrate with respect to @xmath7 and show that @xmath191 , in agreement with the existence of ageing in this system for the position correlations .
one could also be interested in the following correlation : @xmath192 in the large @xmath7 , @xmath2 limit , this correlation is found to be @xmath193 in the thernodynamic limit , @xmath194 , @xmath188 , so that @xmath195 is vanishing in the thermodynamic limit , signaling the absence of ageing for the correlation ( [ c(2 ) ] ) .
the overlap associated to the correlation ( [ c(2 ) ] ) is the local one : @xmath196 .
we can thus calculate @xmath197 using the integration ( [ taurem ] ) over the density of eigenvalues . in order to calculate the behavior of @xmath0
, we consider a finite system of volume @xmath198 with a lattice spacing @xmath102 .
the number @xmath10 of sites is @xmath199 . by thermodynamic limit , we mean that @xmath200 is fixed , @xmath194 and @xmath188 . using ( [ taurem ] )
, we have then @xmath201 where the summation is restricted over the wave vectors @xmath184 such as @xmath202 .
clearly , the dominant contribution to ( [ tau(2 ) ] ) comes from the smallest eigenvalues only .
it is thus legitimate to expand the eigenvalues up to the quadratic order in @xmath184 : @xmath203 one then obtains @xmath204 if @xmath205 , @xmath206 since @xmath188 , @xmath207 in the thermodynamic limit in two dimensions .
now if @xmath208 , @xmath209 since @xmath188 , this expression goes to zero in the thermodynamic limit whatever the dimension .
we have thus shown that , in the case of random walks , and in the case of the correlations ( [ c(1 ) ] ) and ( [ c(2 ) ] ) , the existence / absence of ageing in the thermodynamic limit is consistent with the general argument of section [ general ] .
we assume that the traps have been ordered such as @xmath56 . unlike the rem  like trap model
, we do not make any asumption about the trap distribution .
the finite temperature dynamics is @xmath210 if @xmath17 and satisfies detailed balance .
the zero temperature form of the glauber matrix of this model is very simple since @xmath211 , @xmath212 if @xmath213 and @xmath214 if @xmath215 at zero temperature
. the spectrum is @xmath216 .
the eigenvector @xmath217 for the eigenvalue @xmath22 is @xmath218 .
it is easy to show that the eigenvector @xmath219 associated to the eigenvalue @xmath220 is @xmath221 we can now express the pure states in terms of eigenvectors @xmath222 this means that , at zero temperature , a pure state @xmath223 [ @xmath66 is the @xmath66-ieth traps , with @xmath26 ] overlaps only with eigenstates associated to eigenvalues @xmath225 , and this overlap does not depend on the eigenstate .
compared to the result of section [ eigi ] , the structure of the decomposition of a pure state into eigenstates of the glauber dynamics is very different from the rem  like trap model .
if one starts from a pure state @xmath226 at time @xmath6 , the state of the system at time @xmath7 is given by @xmath227 we recover the fact that this model can reach its ground state in a time of order @xmath10 at zero temperature , whereas in the limit of zero temperatures , and for the rem  like trap model , the time spent in a trap @xmath18 is proportional to @xmath228 , which goes to infinity if @xmath229 and @xmath10 fixed .
the two times autocorrelation functions were calculated in @xcite by means of a laplace transform formalism .
the result is @xmath230 if @xmath231 .
we rederive this result in the appendix using a treatment different from @xcite , and based only on the diagonalization of the glauber matrix .
the bm dynamics thus exhibits ageing .
provided only self overlaps are non zero ( @xmath232 ) , the correlation function introduced in section [ general ] is @xmath233 , provided @xmath234 .
the typical time @xmath235 below which ageing occurs is thus logarithmically divergent . on the other hand , this can be also shown directly from the spectrum : @xmath236 which diverges logarithmically .
as far as the finite temperature behavior of the two times autocorrelation function is concerned , this model exhibits interrupted ageing @xcite since at any finite temperature the correlation length is finite , of the order of @xmath237 , ageing does not exist after a time independent of the system size .
it is thus interesting to apply to this case the general argument of section [ general ] .
we first begin with the definition of the glauber dynamics .
we use the dynamics introduced in @xcite .
the master equation of this single spin flip dynamics is @xmath238 as a consequence of detailed balance , @xmath239 notice that , unlike the case of the rem  like trap or bm models , there is no need to multiply the matrix elements by @xmath240 .
it is shown in @xcite that in the infinite temperature limit , the eigenvalues of the glauber matrix are negative integers between @xmath241 and @xmath95 , with degeneracies given by the binomial coefficients .
the largest relaxation time in the infinite temperature limit is equal to unity whatever the system size , and there is thus no need to rescale the coefficients of the glauber matrix .      in the case of a spin model ,
the overlap @xmath242 is given by @xmath243 since ageing is interrupted , we should observe a saturation of the time @xmath244 as @xmath10 increases .
moreover , @xmath244 should increase as the temperature decreases . in practise
, we work with an open chain and we take advantage of the existence of the global @xmath245 symmetry ( invariance of the glauber dynamics under the @xmath246 transformation in a zero magnetic field ) , and the reflection symmetry of the chain .
the glauber matrix is diagonalized in each of these four symmetry sectors .
the results are plotted on figure [ fig5 ] .
the difference @xmath247 decreases as @xmath10 increases , in agreement with the existence of interrupted ageing .
this model was proposed in 1975 as an `` exactly solvable '' model of spin glasses @xcite . for reviews on this model and its developments ,
we refer the reader to the reviews @xcite .
the hamiltonian of the sk model is @xmath248 where the sum is carried out over all the pairs of sites ( infinite dimensional model ) and the quenched random interactions are distributed according to a gaussian distribution @xmath249 this model is a spin glass if @xmath250 and a paramagnet otherwise .
the two times correlation function were shown in @xcite to exhibit ageing in the glassy phase .
we are going to apply our general argument of section [ general ] to the case of the rem  like trap model .
to do so , we calculate numerically @xmath251 for a large number of disorder realizations .
the result is plotted on figure [ fig6 ] , as well as the differences @xmath252 for @xmath253 sites .
we observe that @xmath254 increases at low temperatures and decreases at high temperature .
the crossing point on the insert of figure [ fig6 ] should be identified to the glass temperature .
this temperature is not equal to unity since the system sizes are very small .
notice that even though the sizes are quite small , a large number of disorder configurations is required to have a significant statistics ( of the order of 500000 independent disorder configurations ) .
the computer time required was equivalent to 1000 hours of cpu time on a single processor .
our results are thus consistent with the existence of ageing in the low temperature phase , and no ageing in the high temperature phase .
we have thus carried out a detailed analysis of ageing from a glauber dynamics point of view , with an emphasis on finite size effects .
we have shown how to extract information from the glauber matrix about ageing , and how to analyze finite size effects . in some simple cases with only self overlaps , the behavior of the time after which there is no ageing only involves the knowledge of the average density of eigenvalues @xmath173 when @xmath87 .
this is the case for the rem  like trap model , in which case we have carried out a careful analysis of the density of small eigenvalues .
our analysis was confirmed by numerical diagonalizations of glauber matrices of small clusters . in this case
, we were able to related the typical behavior of the density of eigenvalues to the spectral dimension of the corresponding diffusion process .    :
r.m . acknowledges interesting discussions with a. barrat , j.p .
bouchaud , s. franz and n. wingreen .
the authors acknowledge b. douot for his interest to this work and his encouragements .
r.m . also acknowledges the hospitality of nec research institute at princeton where part of this work was done .
we are going to calculate the two times autocorrelation function @xmath255 using a method based only on the exact diagonalization of the glauber matrix at a zero temperature . starting from an infinite temperature state at time @xmath6 ,
@xmath256 with @xmath30 a pure state , the natural definition of the two times autocorrelation function within this model is @xmath257 where @xmath258 is the probability to find the system in the trap number @xmath18 at time @xmath2 and @xmath259 is the escape rate from trap @xmath18 , that is the diagonal element of the glauber matrix @xmath260 . this definition of the waiting time dependent autocorrelation function is similar to the one used in @xcite . using the expression ( [ eig ] ) of the eigenvectors and the decomposition
( [ dec1 ] ) and ( [ dec2 ] ) of the pure states into eigenvectors , we can calculate in a straightforward fashion the evolution of the system , starting from the state ( [ initial ] ) : @xmath261 |e_k \rangle .\ ] ] the probabilities @xmath258 are easily recovered : @xmath262 . after straightforward calculations , the dominant behavior of @xmath255 is obtained , in the limit @xmath263 : @xmath264 which is equivalent to the result of @xcite .
notice that if @xmath2 is of the order of the number of traps , the system has reached equilibrium : the dynamics is stationary ( the system is stuck in the lowest trap ) and there is no ageing .
m. alba , m. ocio , j. hamman , europhys .
lett . , 45 ( 1986 ) ; j. phys .
( 1985 ) l-1101 ; m. alba , j. hamman , m. ocio , p. refregier , j. appl .
phys . , 3683 ( 1987 ) ; e. vincent , j. hamman , m. ocio , _ recent progress in random magnets _ , d.h .
( world scientific , singapore , 1992 ) .
m. mzard , g. parisi and m. virasoro , _ spin glass theory and beyond _ , world scientific ed .
( 1993 ) ; k.h .
fisher and j.a .
hertz , _ spin glasses _ , cambridge studies in magnetism ( 1991 ) ; v.s .
dotsenko , m.v .
feigelman and l.b .
ioffe , _ spin glasses and related problem _ , physics reviews , hardwood academic publishers ( 1990 ) .
density of eigenvalues in the arrhenius regime for the rem  like trap model and fit to the exponential form ( @xmath265,@xmath266 traps ) .
the insert shows the same quantity , but with also the smallest eigenvalues .
for these small eigenvalues , the arrhenius picture is not valid .
density of small eigenvalues of the models of section [ gm ] , for @xmath267 .
( a ) : infinite dimensional case ( rem  like trap model ) ; ( b ) : one dimensional case , with `` disordered '' traps ; ( c ) : one dimensional case with `` ordered '' traps . the lines in cases ( b ) and ( c ) are parallel , and in good agreement with the arrhenius law , whereas we observe a significant deviation to the arrhenius law in ( b ) case .
exponent @xmath90 deduced from the fit ( [ lawarr ] ) of @xmath85 .
the circles correspond to the infinite dimensional case ( rem  like trap model ) , the squares to the one dimensional `` disordered '' case , and the diamonds to the `` ordered '' one dimensional case .
the solid line corresponds the the arrhenius regime ( in other words , to site  disordered random walks above the upper critical dimension 2 ) , and the dashed line corresponds to the one dimensional renormalized arrhenius regime ( in other words , to site  disordered random walks in one dimension ) .
structure of the eigenstates of the rem  like trap model for 20 traps .
the label @xmath18 is an eigenvalue label , with eigenvalues in growing order and @xmath19 labels the traps ( the depth increases if @xmath19 increases ) .
( a ) represents @xmath268 for the @xmath269 pure states .
( b ) represents the depth of the traps .
( c ) represents the 20 eigenvalues between @xmath270 and @xmath95 for 4 different temperatures ( circles : @xmath271 , squares : @xmath267 , diamonds : @xmath272 , triangles : @xmath273 ) .