in the last decade , studies about the structure and dynamics of complex networks have blossomed , thanks in particular to the versatility of the network representation , which has turned out to be adequate for systems as diverse as the internet or social networks .
a large body of knowledge about the empirical description of networked systems has thus been accumulated , together with a wealth of modeling techniques ; a good level of understanding of how dynamical processes taking place on networks depend on their structure has been as well reached @xcite .
many network studies have been concerned with systems of interest in several scientific areas a priori remote from physics ( social sciences , biology , computer science , epidemiology , ... ) , and they have also reached more traditional fields of statistical physics , such as the study of glassy systems , as we now describe .
the many puzzles raised by the glass transition , and in particular the slow dynamics displayed by glassy systems at low temperatures , have been the subject of a large interest in the past decades @xcite .
one of the approaches which has led to promising insights consists in the description of the dynamics of a glassy system inside its configuration space .
the energy landscape of a glassy system is typically rugged , made of many local minima ( metastable states ) , whose huge number makes it difficult to reach equilibrium . in this framework ,
the energy landscape is seen as a set of basins of attractions of local minima ( `` traps '' ) , and the system evolves through a succession of harmonic vibrations inside traps and jumps between minima @xcite .
this picture has stimulated the definition and study of various simplified models of dynamical evolution between traps , in order to reproduce the phenomenology of glassy dynamics @xcite .
on the other hand , several studies have focused on obtaining a better understanding of the structure of these local minima .
a way to attain this goal is to perform numerical simulations of small systems , at a fixed temperature , and quenching them at regular time intervals in order to make them reach the nearest local minimum .
information is then gathered on the various local minima , and on the sizes of their basins of attraction .
various studies have investigated , among other issues , the detailed structure of the potential energy landscape , the substructure of minima , and the properties of energy barriers between minima @xcite .
several works have also used the information on the energy landscape to study a master equation for the time evolution of the probability to be in each minimum .
the considered systems range from clusters of lennard - jones atoms to proteins or heteropolymers @xcite .
an interesting property of the modeling of the energy landscape in terms of a set of traps linked by energy barriers , lies in the possibility to define and study its network representation within the context of network theory . in this representation
, each local minimum is associated to a node , and a link is drawn between two nodes whenever it is possible for the system to jump between the basins of attraction of the corresponding minima .
the links can then be defined as weighted and directed , as jumps between minima are not equiprobable , and may be easier in one direction than in the other .
networks of local minima of the energy landscape have thus been built and studied .
these networks have been found to exhibit a small - world character @xcite .
the number of links of each node ( its degree ) turns out to be strongly heterogeneous , possibly with scale - free degree distributions , which have been linked to scale - free distributions of the areas of the basins of attraction @xcite .
complex network analysis tools have also been used to investigate the structure of energy landscapes of various systems of interest , such as lennard - jones atoms , proteins , or spin glasses , among others @xcite .
the energy of a minimum and its degree ( i.e. , the number of other minima which can be reached from this minimum ) have been shown to be correlated , as well as the barriers to overcome to escape from a minimum .
in particular , a logarithmic dependence of the energy of a minimum on its degree has been exhibited , as well as energy barriers increasing as a ( small ) power of the degree of a node @xcite .
no systematic study of these issues has however been performed , and most investigations have been limited to relatively small systems because of computational limitations .
most importantly , the investigations cited above have focused on the topology of the network of minima , conceived as a tool to characterize the energy landscape .
the structure of a network has however a deep impact on the properties of the dynamical processes which take place on it @xcite .
it seems thus adequate to put to use the tools and techniques developed for the analysis of dynamical processes on networks to achieve a better understanding of how the energy landscape structure , represented as a network , affects the system performing a random walk in it , and how the onset of glassy dynamics can be described in this way in a general framework .
in a previous paper @xcite , we have made a first step to fill this gap by focusing on the trap model put forward in ref .
@xcite . in this paper
, we generalize our approach to more involved transition rates between energy minima .
we show how the heterogeneous mean - field ( hmf ) theory @xcite can be used in this context to highlight the connection between the topological properties of the network of minima and the dynamical exploration of these minima .
we show in particular that the relationship between energy and degree of the minima is a crucial ingredient for the existence of a transition and the subsequent glassy phenomenology .
our results shed light on the empirically found relationship between the energy of a local minimum and its degree , and we hope that they will stimulate more systematic investigations on this issue .    we have organized our paper as follows : in section  [ sec : model ] we define our model of energy landscape dynamics as a random walk on a complex network .
different physical transition rates are proposed , and the corresponding numerical implementation is discussed . in section  [ sec : formalism ] we present a theoretical analysis based on the heterogeneous mean - field approximation for dynamical processes on complex networks .
this formalism is applied in section [ sec : general - formalism ] , where general analytical approximate expression are presented for the main quantities characterizing the glassy transition and dynamics .
these expressions are applied to the different physical transition rates considered in section  [ sec : appl - phys - trans ] , where checks against numerical simulations are also presented . in section  [ sec : barriers ] we discuss the relation between energy basins and energy barriers .
finally , in section  [ sec : conclusions ] we present our conclusions .
we consider a network of @xmath0 nodes , in which each vertex @xmath1 corresponds to a minimum in the energy landscape , and a link is drawn between two minima @xmath1 and @xmath2 if the system can jump directly from @xmath1 to @xmath2 . to each node @xmath1
is associated the energy @xmath3 of the corresponding minimum ( energies are defined from a reference level , in such a way that @xmath4 for all @xmath1 ) . moreover , an energy gap @xmath5 is associated to the edge between vertices @xmath1 and @xmath2 , as depicted in fig .  [ landscape ] : @xmath5 is a symmetric function , such that the energy barrier that must be overcome to jump from vertex @xmath1 to vertex @xmath2 can be written as @xmath6 and , analogously , @xmath7 . obviously , we will have in general @xmath8 .
are measured positive downwards .
energy gaps @xmath5 are defined positive upwards.,width=340 ]    the system under investigation is pictured as a walker exploring the network through a biased random walk .
the rate ( probability per unit time ) @xmath9 to go from vertex @xmath1 to vertex @xmath2 depends a priori on the energy at vertices @xmath1 and
@xmath2 and/or the energy barrier between @xmath1 and @xmath2 that must be overcome . the random walk model is defined in discrete time @xmath10 as follows :    * at time @xmath10 , the walker is in vertex @xmath1 .
* it chooses at random a neighbor of @xmath1 , namely @xmath2 . * with a probability @xmath9 , that depends on the energy @xmath11 and/or on the energy barrier @xmath12 , the walker hops to vertex @xmath2 . * time is updated @xmath13 .
the relationship between the probabilities @xmath9 and the energy and energy barriers can be of different forms . in usual unbiased random walks , @xmath9 is a constant independent of both @xmath1 and @xmath2 @xcite . as a first step to introduce a dependence on the nodes ,
a possible approach is that in which the energy barriers depend only on the local minima themselves , i.e. we consider @xmath14 .
for example , in the bouchaud trap model considered in ref .
@xcite , the probability to exit from a trap is just an arrhenius law depending only on the departing trap s depth , namely @xmath15 where @xmath16 is the inverse temperature and @xmath17 is a constant that determines a global timescale .
other possible definitions include the metropolis one @xmath18 and the glauber rate @xmath19 we note that the rates considered in the bouchaud trap model are quite different from the metropolis or glauber rates . indeed , while the former depends only on the depth of the originating trap , the latter depend also on the energy of the arriving vertex .
this translates in the fact that , in the limit of zero temperature , the dynamics of the bouchaud trap model is frozen for any @xmath11 , while metropolis and glauber dynamics still allow jumps to lower energy minima @xcite . within an even more realistic representation of glassy dynamics , one can also contemplate the case @xmath20 , allowing for the transition rates to depend explicitly on the energy barriers between adjacent minima . as a paradigm of this choice , we propose a rate of the arrhenius form @xmath21 which acts a straightforward generalization of the local transition rate ( [ eq : rate_traps ] ) .
the case of rate ( [ eq : rate_traps ] ) ( local trapping ) was studied in a previous publication @xcite . in the following
, we will consider in turn non - local rates ( [ eq : rate_metropolis ] ) , ( [ eq : rate_glauber ] ) and ( [ eq : rate_barriers ] ) and discuss the fundamental differences due to the introduction of energy barriers in the model .
we emphasize that our model differs both from usual unbiased random walks , as the local energy determines the transition rates , and from mean - field trap models in which jumps between any pair of energy minima are a priori possible : here , the system can jump only between neighboring nodes .
the dynamical evolution depends therefore both on the network topology and on the energies associated to the nodes .      to implement numerically the random walk , it is convenient to resort to the techniques developed for general diffusion processes on complex weighted networks @xcite .
the main advantage of this method is to avoid rejection steps , thus improving dramatically the computational efficiency @xcite .
therefore , at each simulation step the random walker sitting at node @xmath1 selects a neighbor @xmath2 with probability @xmath22 , where the sum in the normalizing factor is extended to all of @xmath1 s neighbors . as the walker hops on node @xmath2 , the physical time is incremented by an interval @xmath23 drawn from the exponential distribution @xmath24 , where @xmath25 is the inverse of the average escape rate out of node @xmath1 . in this way
the simulation time is disentangled from the physical time and the latter has no impact on the simulation efficiency .
no matter how much physical time a walker spends in a node , from the simulation time point of view it is always just a time step .    the network substrates on which we will focus are scale - free networks with a degree distribution of the form @xmath26 and @xmath27
. we will generate them with the uncorrelated configuration model ( ucm ) @xcite , that allows us to tune the degree distribution to the desired form and prevents the formation of degree - degree correlations .
networks are therefore generated as follows : a number of stubs ( or semi - links ) extracted from the desired final degree distribution is assigned to each node .
stubs are then randomly paired to form links between nodes , with the prescription that multiple links as well as self - loops must be avoided . a minimal degree @xmath28 is fixed a priori . to avoid spurious effects due to the possible presence of tree - like structures
@xcite it is convenient to adopt @xmath29 .
we will choose @xmath30 in all of our simulations .
so far the algorithm coincides with the one of the configuration model @xcite , but the ucm introduces moreover a cutoff to the degree distribution , @xmath31 , which avoids the formation of degree correlations by limiting the size of the hubs @xcite .
in order to gain analytical understanding of the role of the different transition rates in the corresponding glassy dynamics , we apply a standard heterogeneous mean - field ( hmf ) formalism @xcite .
the basic tenet of hmf is the assumption that all the dynamical properties of a vertex depend only on its degree .
vertices are thus grouped into classes according to their degree , and vertices with the same degree are treated as equivalent .
this approximation is consistent with previous findings that have uncovered the correlations between the energy of a local minimum and the degree of the corresponding node in the network @xcite .
we therefore make the assumption that there exists a relationship @xmath32 where the function @xmath33 is a characteristic of the model .
this also means that the distributions of energies @xmath34 of the system s landscape , and the degree distribution @xmath35 of the corresponding network are linked through @xmath36 . in the same spirit
, we make the further assumption that the energy gap between minima @xmath1 and @xmath2 depends only on the degrees of @xmath1 and @xmath2 , i.e. , that it can be written as @xmath37 , where @xmath38 is a symmetric function of @xmath39 and @xmath40 .    under the hmf approximation
the dynamics will thus focus on the transitions between different degree classes .
the rate to go from a vertex @xmath39 to a vertex @xmath40 can be written as @xmath41 the function @xmath42 , defined as the conditional probability that a vertex of degree @xmath39 is connected with another vertex of degree @xmath40 @xcite , takes into account the topological features of the network , by gauging the probability of selecting a vertex @xmath40 as neighbor of @xmath39 .
the function @xmath43 measures the rate of jumping from a vertex of degree @xmath39 to a vertex of degree @xmath40 ( given that they are connected by an edge ) , and depends on @xmath39 and @xmath40 through the rates @xmath44 and the functions @xmath36 and @xmath45 . obviously , the rate @xmath43 is not in general a symmetric function of @xmath39 and @xmath40 .
it is worth noting that , apart from a normalization , equation  ( [ eq : weighted_propagator ] ) is simply the so - called weighted propagator describing the probability that a node in class @xmath39 interacts with a node in class @xmath40 @xcite .
we also note that the rates @xmath43 depend on the inverse temperature @xmath46 through the microscopic rates @xmath9 .
in this section , we apply the hmf theory to compute different quantities relevant for the characterization of the dynamics of a random walk in a complex energy landscape represented in terms of a network of minima .
the description of a random walk dynamics starts from the occupation probability @xmath47 , defined as the probability for the walker to be in _ any _ node of degree @xmath39 at a time @xmath48 .
its time evolution can be easily represented in terms of a master equation of the form @xmath49 upon describing the state at time @xmath48 with the row vector @xmath50 , where @xmath51 is the cutoff or largest degree in the network , equation ( [ rate_scalar ] ) can be rewritten in vector form as @xmath52 where the matrix @xmath53 , with elements @xmath54 is a generalization of a laplacian matrix to the case of directed weighted graphs .
the matrix elements satisfy @xmath55 which ensures conservation of probability and states that the columns of @xmath53 are not linearly independent .
the real part of every eigenvalue of @xmath53 is non - negative @xcite . as a consequence , all solutions of eq .
( [ rate_vectorial ] ) , which can be formally written as @xmath56 are stable according to lyapunov criteria .
in particular , since @xmath57 , @xmath53 always has the eigenvalue @xmath58 , which corresponds to a constant solution of the problem . at this point
one can proceed in close analogy with discrete - time regular markov chains @xcite . by making the assumption that the matrix @xmath59 is non - negative and irreducible ( indeed it is for every choice of @xmath60 in the following )
, we can prove that the @xmath58 eigenvalue of @xmath53 has algebraic multiplicity @xmath61 .
hence , the stationary solution of eq .
( [ rate_vectorial ] ) is unique .      in order to calculate the steady solution @xmath62 in the limit @xmath63
, one can impose @xmath64 .
this leads to the condition @xmath65 so that we are left with the task of finding the left nullspace of @xmath53 .
( [ homogeneous_v ] ) is a homogeneous system of algebraic linear equations .
it admits non - trivial solutions since @xmath66 . in our case ,
the solution to ( [ homogeneous_v ] ) can be easily found by imposing the detailed balance condition .
namely , writing eq .
( [ homogeneous_v ] ) as @xmath67=0 ,      \label{eq:1}\ ] ] we can obtain a solution by imposing that the terms inside the summation in eq .
( [ eq:1 ] ) cancel individually , that is @xmath68 substituting the form of @xmath69 , we obtain @xmath70 where in the last step we have used the degree detailed balance condition @xmath71 which simply expresses that the number of edges from a node of degree @xmath39 to a node of degree @xmath40 is equal to the number of edges from a node of degree @xmath40 to a node of degree @xmath39 @xcite . from eq .
( [ eq:2 ] ) , we see that its right - hand - side must be expressible as a simple ratio of a function of @xmath39 over a function of @xmath40 . a general way to obtain this is to impose a coarse - grained rate @xmath60 taking the general form @xmath72 in other words , we assume that the rate @xmath60 can be written as the product of a function of @xmath40 , a function of @xmath39 , and a symmetric function @xmath73 ( where @xmath39 and @xmath40 need not be separable ) .
we will see later that all the rates @xmath74 defined in sec .
[ sec : definition ] ( traps , glauber , metropolis , and energy barriers ) can be written in such a form .
the stationary solution is then given by @xmath75 where @xmath76 is a normalizing constant determined by the condition @xmath77 .
such a solution is unique , as proven above .
interestingly , the symmetric function @xmath78 does not enter the steady solution , although it will play a role in affecting the transient behavior , as we will see in the next sections .      the steady state solution found above for the occupation probability
is defined if and only if the normalization constant @xmath79 is finite .
when this condition is met , the random walker reaches an equilibrium state with a distribution @xmath80 . on the other hand ,
whenever such condition is not met , the random walker is unable to reach a steady state , i.e. the steady solution to the rate equation does not correspond to any physical steady state in equilibrium @xmath81 .
we identify this region of the phase space with the glass phase for our random walker @xcite .    the functions @xmath82 and @xmath83
depend on the temperature , on the precise dynamics chosen ( traps , glauber , metropolis ) , and encode the relationship @xmath36 between energy and degree of the minima .
the degree distribution moreover enters explicitly the expression @xmath84 .
as the various parameters of the model are changed , it is thus a priori possible to go from one phase in which @xmath84 is finite to one in which @xmath84 diverges . in a physical system in particular , the control parameter is usually the temperature , while the topology of the network of minima and the function @xmath36 are given .
it is then clear from eq .
( [ eq:4 ] ) that the presence or absence of a finite glass transition temperature @xmath85 , such that @xmath84 becomes infinite for @xmath86 , depends on the interplay between the topology of the landscape network ( as determined by @xmath35 ) and the functions @xmath82 and @xmath83 .
interestingly , at this mean - field level , the existence of a transition does not depend on the network degree correlations , since the conditional probability @xmath42 do not enter eq .
( [ eq:4 ] ) .
let us consider for instance a network of minima with a heavy - tailed degree distribution such as @xmath26 .
a transition between a finite and an infinite @xmath84 can be observed if and only if @xmath87 behaves at large @xmath39 in the form @xmath88 where the exponent @xmath89 depends on the temperature , and can take values smaller or larger than @xmath90 depending on the temperature .
another example is given by a stretched exponential form for @xmath35 , @xmath91 , in which case a transition is observed if and only if @xmath87 is of the form @xmath92 , with @xmath93 a function of the temperature ( the transition is then given by @xmath94 ) .      in any finite system , unless the product function @xmath87 exhibits some sort of singularity , the normalization constant eq .
( [ eq:4 ] ) is finite and the steady state distribution @xmath95 exists , the occupation probability @xmath47 converging to it after an equilibration time , i.e. @xmath96 the corresponding thermalization of the occupation probability occurs in a way depending on the function @xmath36 .
shallow energy minima are indeed explored first , while deep traps ( large @xmath97 ) are visited at larger times @xcite .
if @xmath36 is a growing function of @xmath39 , as indeed found empirically @xcite , small degree nodes correspond to shallow minima , and deeper minima are associated to larger nodes .
the evolution of @xmath47 takes then place in a hierarchical fashion : the small degree region equilibrates first , and progressive equilibration of larger degree regions takes place at larger times . in this respect , we obtain a strong difference between the biased random walk that the glassy system experiences and usual diffusion processes corresponding to unbiased random walks , which visit first large degree vertices and then cascade down towards small degree nodes @xcite , in the present case we observe an `` inverse cascade process '' from small vertices to hubs .
we have found in @xcite that , in the case of a random walk among traps , this hierarchical thermalization is summarized in a scaling form for @xmath47 , which can be written as @xmath98 where @xmath99 represents the maximum degree of the vertices equilibrated up to time @xmath48 , and @xmath100 interpolates between @xmath101 at small @xmath102 and the short time form of @xmath103 which is proportional to @xmath104 .
we will see in the next section that a similar scaling is obeyed for other transition rates .
in general , for the glassy dynamics , the functional form of @xmath99 can moreover be obtained through the following argument : the total time @xmath48 can be written as the sum of the trapping times spent in the vertices that have been visited since the beginning of the dynamics .
trapping times increase with the depth of the minimum , hence with the degree ( we are still considering the case of an increasing function @xmath33 ) , and , in the glassy phase , the consequence is that the sum of trapping times is dominated by the vertex with the largest degree visited up to that point , namely @xmath105 .
moreover , the average trapping time @xmath106 at a given vertex @xmath39 can be estimated as the inverse of the average escape rate @xmath107 from that vertex : @xmath108 we can therefore estimate @xmath99 , the typical degree up to which nodes are `` equilibrated '' at time @xmath48 , by approximating @xmath109 , and solving the equation @xmath110 to obtain @xmath105 as a function of @xmath48 . note that the result depends here on the function @xmath111 and not only on @xmath82 , @xmath83 , @xmath36 and @xmath35 .
the properties of the system can be further quantified by measuring the average time @xmath112 required by the random walker to escape from the vertex it occupies at time @xmath48 @xcite . for small waiting times
@xmath48 , @xmath113 increases as a result of the transient equilibration of @xmath114 . for large @xmath48 , such that @xmath115 is close enough to the equilibrium @xmath116
, @xmath113 can be calculated instead as the average @xmath117 \tau_k\ ] ] where @xmath118 is the inverse of the equilibrium escape rate , cf eq .
( [ eq:7 ] ) , yielding @xmath119 most interestingly , the explicit form of the average escape time @xmath113 depends explicitely on the symmetric function @xmath111 as well as on the network degree correlations , as expressed by the conditional probability @xmath42 .
let us go back to the issue of the existence of a glass transition in the model .
we first recall the phenomenology of the fully connected trap model , with transition rates @xmath120 for any @xmath1 and @xmath2 , where the energies @xmath11 are random numbers extracted from a distribution @xmath34 @xcite .
as all traps are connected with each other , all traps are equiprobable after a jump , so that the probability for the system to be in a trap of depth @xmath97 is simply @xmath34 , and the average rest time spent in a trap is @xmath121 .
a transition between a high temperature phase and a glassy one is thus obtained if and only if , when @xmath46 increases , @xmath122 is finite at small @xmath46 and diverges at a finite @xmath85 . such a phenomenology is obtained if and only if @xmath34 is of the form @xmath123 at large @xmath97 ( else the transition temperature is either @xmath58 or @xmath124 ) , and the transition temperature is then @xmath125 @xcite .    in the present case of a nework of minima ,
the average rest time that the walker spends in a minimum is @xmath126 where the symbol @xmath127 refers to the average performed over the measure @xmath128 , which represents the probability that the walker is in any vertex of degree @xmath39 after a hop .
note that we disregard here the _ physical _ time , which is the sum of times spent in the various minima , but consider only the number of hops between minima . in the case of the traps model , @xmath129
is simply given by the probability to be in a node of degree @xmath39 after a hop in a random walk , i.e. by @xmath130 @xcite , since the transition rates do not depend on the arrival node . in a non - local trapping model instead
, we need to write a master equation of the form @xmath131 where the matrix @xmath132 is now stochastic and the derivative is intended with respect to the number of hops . in the long time limit ,
we impose @xmath133 and calculate @xmath128 as we did for @xmath95 , imposing the detailed balance condition , and obtaining @xmath134   r_k \ , \ ] ] where @xmath135 is a normalization factor , given by @xmath136    we finally obtain for the average @xmath137 : @xmath138 where @xmath84 is the normalization factor of @xmath139 defined in eq .
( [ eq:4 ] ) .
as for the average escape time , the average rest time @xmath137 thus depends on all the parameters of the system , including the network s degree correlations and the symmetric function @xmath140 .
in this section , we apply the general hmf results obtained in the previous section to physical transition probabilities between local minima given by the trap model , glauber , metropolis and barrier - mediated rates .
we will focus for definiteness on scale - free networks characterized by a power - law degree distribution @xmath141 with @xmath142 , which turns out to be the interval reported in the literature @xcite .
let us first consider the explicit form of the transition rates in each case , to show that they can be cast in the form given by eq .
( [ eq:3 ] ) . in the case of the trap model ,
the rate to jump from a vertex @xmath39 to a vertex @xmath40 is simply @xmath143 , where we recall that @xmath33 gives the depth of a node of degree @xmath39 : it depends only on the degree of the starting node , and not on the node reached after the jump .
we can therefore use @xmath144 the glauber rate can be written as @xmath145 leading to @xmath146 the metropolis transition , on its turn , reads @xmath147.\ ] ] since , for positive @xmath89 , @xmath148 , we can choose @xmath149 finally , in the presence of energy barriers , the transition rate reads @xmath150 where @xmath151 is a symmetric function of its arguments , so that we can use @xmath152      interestingly , for all the transition rates considered above , the product of the functions @xmath153 and @xmath83 , which controls the existence of the steady state solution of the occupation probability , takes the form @xmath154 the normalization constant @xmath84 defined in eq .
( [ eq:4 ] ) can thus be written as @xmath155 for a power - law degree distribution @xmath26 , a finite glass transition temperature is then obtained if and only if @xmath36 is of the form @xmath156 which is precisely what has been found , in conjunction with a scale - free degree distribution , in ref .
@xmath84 is then indeed given by a sum of terms of the form @xmath157 , which converges if and only if @xmath158 in other words , a transition between a high temperature phase in which @xmath95 exists and a low temperature glassy phase is obtained at the critical temperature  @xcite @xmath159    quite noticeably , the existence of a transition at a finite temperature , as well as the value of this temperature , does not depend on the form of the transition rates between the local minima , but only on the existence of a particular interplay between the topology of the network of minima and the relationship between energy and degree in this network , as determined by the function @xmath36 .
we emphasize that this result is also independent of the network degree correlations @xmath42 , as already noted in the previous section .
let us focus on the case of a scale - free network of minima , with @xmath160 and @xmath161 .
for any of the rates discussed above , the steady state measure , when it exists , is given by @xmath162 where @xmath163 is the riemann @xmath163 function .
a plot of @xmath164 as a function of @xmath39 and @xmath165 is given in fig .
[ q_k ] , while data from simulations are reported in fig .
[ q_k_sim ] for the evolution of @xmath115 under glauber dynamics . above the transition ( @xmath166 ) ,
low-@xmath39 states ( i.e. , shallow minima ) are more probable .
as the temperature decreases , @xmath95 becomes less and less peaked at low values of @xmath39 , and large-@xmath39 states , which correspond to lower energies , become more and more probable .
for the random walker to be in any node of degree @xmath39 . for @xmath167
the system undergoes a transition to a glassy state.,width=340 ]    in any finite system , the sum defining @xmath84 , is finite at any temperature as the degree distribution has a cut - off at a finite @xmath51 : @xmath168 for instance , for @xmath169 , and with @xmath170 , the sum @xmath171 is analytic in @xmath172 for any finite @xmath51 . here
@xmath173 is the harmonic number of order @xmath174 , which tends to @xmath175 for @xmath176 .
the probability @xmath95 is thus well defined for every @xmath165 and for any finite system . in particular , performing a continuous degree approximation in eq .
( [ eq:19 ] ) , we can obtain an estimate of the network size dependence of @xmath84 as @xmath177 for @xmath178 , @xmath84 tends to a constant as the network size ( and thus @xmath51 ) increases . on the other hand , for @xmath179 , @xmath84 diverges as @xmath180 , i.e. , as @xmath181 in uncorrelated scale - free networks , which obey @xmath182 .     for glauber dynamics . the distribution measured after a small waiting time
@xmath48 is determined by a usual unbiased random walk behavior , i.e @xmath183 , while at larger times it relaxes to the equilibrium @xmath95 . the relaxation towards equilibrium starts from small degree nodes .
data refer to ucm networks with @xmath184 and @xmath185 for @xmath186 ( @xmath187 ) : with these parameter values , for small times @xmath188 , while at large enough times @xmath189.,width=340 ]      at low temperatures , even for a finite system , the evolution of @xmath190 towards @xmath95 is slow , as displayed in fig .
[ q_k_sim ] , and an ageing regime takes place , in which the function @xmath115 obeys the scaling form @xmath191 where the characteristic degree @xmath105 can be estimated from eq .
( [ eq:8 ] ) . in order to simplify its computation
, we will consider an uncorrelated network of minima , such that @xmath192 , and we will work in the continuous degree approximation , using the normalized form @xmath193 , where @xmath28 is the minimum degree present in the network .    in the case of the glauber dynamics , the escape rate can be expressed , within the above approximations , as @xmath194 ,   \label{eq:21}\end{aligned}\ ] ] where we have used the relation @xmath195 and where @xmath196 $ ] is the gauss hypergeometric function .
using the asymptotic expansion for @xmath197 @xcite , we obtain that the leading behavior for large @xmath39 yields @xmath198 which leads to @xmath199 from here , using the relation @xmath109 , we obtain @xmath200     at different temperatures ( glauber dynamics ) .
data refer to ucm networks with @xmath184 and @xmath201 , so that @xmath202 ( where we have taken @xmath187 ) .
the top panel presents data for @xmath203 while both middle and bottom panels concern the low temperature case @xmath204 .
accordingly , for the top panel we use @xmath205 for the rescaling , while for the central and the bottom ones it holds @xmath206 ( see eq .
( [ eq:24 ] ) ) .
the curves corresponding to different @xmath48 collapse well under this rescaling .
note that we use rather small values of @xmath48 , because the equilibration time , defined by @xmath207 , is @xmath208 for @xmath204 and @xmath209 for @xmath210 .
each curve is obtained by averaging over @xmath211 simulation runs.,width=340 ]    in fig .
[ collapsep ] we check the validity of eqs .
( [ eq:25 ] ) and ( [ eq:24 ] ) by performing a data collapse analysis for different values of @xmath48 .
the curves obtained for different @xmath48 collapse indeed as predicted .    in the case of the metropolis transition rates , a similar analysis yields @xmath212 leading to the same asymptotic behavior as in eq .
( [ eq:22 ] ) , and therefore to the same scaling picture as for the glauber rate .    as pointed out for the case of the traps model @xcite , however , in finite systems the scaling relations in eq .
( [ eq:24 ] ) hold only as far as @xmath213 , i.e. it exists an equilibration time @xmath214 , obtained by inverting eq .
( [ eq:24 ] ) , above which the system has completely relaxed and eq .
( [ eq:25 ] ) is no more valid .
finally , it is worth stressing that , while for large temperatures the scaling exponent relating @xmath105 and @xmath48 depends on the temperature , in the low temperature phase it becomes independent , being just proportional to the transition temperature .
we note that this saturation of the exponent at @xmath215 is very different from the phenomenology obtained in the trap model @xcite , for which @xmath216 .
an immediate consequence is that the equilibration time strongly depends on @xmath46 , as @xmath217 , for a system described by traps , but is given by @xmath218 for any @xmath204 for glauber and metropolis rates .
contrarily to results for the glass transition temperature @xmath219 and the steady state , the glassy dynamics for barrier - mediated rates does not yield the same results as for glauber and metropolis rates , since @xmath107 does depend on the symmetric function @xmath220 . in particular , we need here to choose a functional form for @xmath45 .
we propose to use @xmath221 which will be justified in section [ sec : barriers ] . in this case
we obtain @xmath222 in order to keep an interesting phenomenology , the constant @xmath223 can not be chosen arbitrarily . if @xmath223 were independent of the system size , the escape rate would be dominated by the exponential at all temperatures .
this behavior would reflect the fact that in this case the rate is suppressed in transitions involving nodes of large @xmath39 , eventually generating unphysically large barriers at large @xmath51 . to prevent the system from building up infinite barriers
, one can impose @xmath224 such that the maximum barrier is always comparable to the lowest energy minimum and neither term dominates the other . as a consequence , we take @xmath225 where @xmath226 is now constant and size independent .
contrarily to the previous cases , @xmath99 is hard to determine , as no explicit inversion of eqs .
( [ eq:8 ] ) and ( [ eq:29 ] ) can be provided for the range of parameters of interest in our study .
a numerical evaluation of @xmath99 , is reported in figure [ f : kw_barriers ] .
the maximum degree of equilibrated vertices @xmath105 has an initial power law increase in time , which is reminiscent of the local trapping model .
however , as larger degree nodes are equilibrated , exponential barriers come into play and the hierarchical thermalization becomes logarithmic in time .     for barrier - mediated dynamics .
curves are obtained from the numerical inversion of eq .
( [ eq:8 ] ) .
values of parameters are : @xmath201 , @xmath227 , @xmath228 , @xmath229 , @xmath230 . ,
width=340 ]      , so that @xmath231 ( @xmath187 ) .
the scaling forms of eq .
( [ eq : escape ] ) produce a collapse of the curves concerning different systems sizes in the three regimes of high , intermediate and low temperature ( top , center and bottom panels , respectively ) .
the slightly worse collapse obtained for @xmath232 may be due to logarithmic corrections as @xmath46 is close to both @xmath85 and @xmath233 . each point
is averaged over @xmath234 simulation runs ( @xmath235 runs on each of @xmath235 network realizations ) .
, width=340 ]    the average escape time , defined as the average time required by the system to escape from the vertex it occupies , can be computed in the long time limit from eq .
( [ eq : t_esc ] ) , as a function of the average trapping time @xmath106 in vertices of degree @xmath39 . from the asymptotic expansions of @xmath106 in eq .
( [ eq:23 ] ) , valid for glauber and metropolis dynamics , evaluation of eq .
( [ eq : t_esc ] ) allows us to observe that , whenever a finite @xmath219 exists , @xmath236 diverges at @xmath237 in an infinite system , as was already observed in the case of local trapping @xcite .
it is noticeable that the same divergence temperature is obtained , as eq .
( [ eq : t_esc ] ) a priori involves the network s degree correlations and the function @xmath140 . within the continuous degree approximation
, the divergence of the escape time with the system size follows the scaling laws @xmath238 as noted in the previous paragraph for the equilibration time , we note that the scaling for @xmath239 differs from the form @xmath240 encountered for local trapping @xcite .
figure [ f : escape ] reports simulation data that confirm the validity of eq .
( [ eq : escape ] ) .
as the temperature is lowered , the initial transient becomes longer , but for large enough times @xmath48 the asymptotic behavior predicted in eq .
( [ eq : escape ] ) is reached , as made clear from the collapse of curves concerning different system sizes .    in the case of barrier - mediated dynamics , @xmath106 depends on the symmetric function @xmath151 , as expressed in eq .
( [ eq:29 ] ) , namely @xmath241 proceeding as above we obtain @xmath242e_0 } & \beta_c/(2+\epsilon ) <        \beta < \beta_c \\ \mbox{const . }   & \beta<\beta_c/(2+\epsilon )        \\      \end{array }    \right .   , \end{aligned}\ ] ] where we recall that @xmath226 does not depend on the system size .      , as a function of the inverse temperature , for different system sizes .
data are obtained by numerical computation of eq .
( [ tau_zi ] ) ( with @xmath187 ) .
note the exponential increase with @xmath46 for @xmath243 , which saturates for @xmath244 as predicted by eq .
( [ e : rest]).,width=340 ]    , as a function of the degree distribution cut - off and for different temperatures .
data are obtained by numerical computation of eq .
( [ tau_zi ] ) .
@xmath245 grows as a power - law of @xmath51 , with an exponent that grows as @xmath46 increases ( going from bottom to top in the figure ) .
the thick - gray line corresponds to @xmath246 , while the thick - black line corresponds to @xmath247 . for larger values of @xmath46 ,
the power - law behavior corresponds to the predicted @xmath248 , which no more depends on @xmath46 .
the @xmath249 curve is reported as a dashed line for reference .
values of @xmath46 represented here are comprised between @xmath250 and @xmath251.,width=340 ]    the hmf expression for the asymptotic average rest time , defined as the average time spent by the system in a minimum , is given by eq .
( [ tau_zi ] ) , namely @xmath252 , where the quantities @xmath84 and @xmath135 , for uncorrelated scale - free networks and a degree - energy relation @xmath253 , take the form , in the continuous degree approximation , @xmath254    let us first recall the case of the local trap model . both @xmath255 and @xmath111
are then constants , so that @xmath256 .
thus , the average rest time behaves as @xmath84 : it is finite for @xmath210 , and diverges with the system size as @xmath257 for @xmath204 , signaling the emergence of the glassy regime at low temperatures .    in the case of glauber and metropolis dynamics ( both leading to the same results ) ,
the situation is more involved , since the product @xmath258 entering @xmath135 is not constant .
in fact , @xmath135 diverges with @xmath51 for @xmath259 , that is , at a lower temperature given by @xmath260 .
the interplay of these two temperatures determines the behavior of the system for finite sizes within the low temperature phase .
in particular , for the glauber dynamics with @xmath261 and @xmath262 $ ] , we have @xmath263 upon considering lower values of @xmath46 , @xmath137 first encounters the divergence of @xmath84 at @xmath85 , which is then partially regularized by the divergence of @xmath135 at @xmath233 . from these results ,
we obtain the emergence of three scaling regimes for the behavior of @xmath137 as a function of the system size : @xmath264     for different transition rates , for a random walker on ucm networks with @xmath265 ( @xmath187 ) and @xmath184 .
discrete points : simulation results ; continuous lines : theoretical predictions from @xmath266 , based on simulation parameters .
the glauber and metropolis transition rates induce two changes of behavior at @xmath85 and at @xmath233 , the first being a steep increase of the average rest time @xmath137 and the second a smoothing / saturation of this increase .
no saturation of @xmath137 is instead observed when barriers are present .
the agreement with theoretical predictions is remarkable , thus corroborating the validity of the hmf assumptions .
moderate deviations are found only in the case of barriers , where exponential growth is expected to add greater fluctuations . in the inset ,
the difference between the two behaviors is more evident thanks to a different scale of the plot .
note that , since in the simulations @xmath187 , the high temperature limit of the rest time is different for the glauber and metropolis dynamics , being @xmath267 and @xmath268 respectively . for the case of barriers
we have chosen @xmath269 .
each point is obtained by averaging the rest times corresponding to the first @xmath270 hops of the random walker in each of @xmath271 network realizations .
, width=340 ]    the direct numerical computation of eq .
( [ tau_zi ] ) is shown in figs .
[ 2temperatures ] and [ kcscaling ] , showing the validity of this analysis .
in particular , the exponential increase of @xmath272 with @xmath46 in the intermediate temperature range @xmath273 is clearly apparent in fig . [ 2temperatures ] , and fig .
[ kcscaling ] confirms that , for @xmath274 , the exponent in the scaling law for the system size @xmath51 does not depend on the temperature .
while the temperature @xmath85 signals the onset of the low temperature phase with glassy dynamics for all considered transition rates , for glauber / metropolis dynamics the low temperature phase can be further divided into two regions that correspond to different behaviors of the timescales with the system size .
( @xmath187 , @xmath231 ) .
data for different system sizes collapse
well when rescaled according to the theoretical values of eq .
( [ e : rest ] ) . while the agreement is excellent both for high and low temperature ( top and bottom panels , respectively ) , logarithmic corrections are probably present for the regime of intermediate temperatures @xmath275 ( central panel ) . in each simulation
run the rest interval starting before @xmath48 and ending after @xmath48 is considered , and each point in the figure is averaged over @xmath234 simulation runs ( @xmath235 runs on each of @xmath235 network instances ) .
, width=340 ]    figures [ simul_2temperatures ] and [ f : rest ] moreover show the result of numerical simulations of random walkers on scale - free networks for glauber and metropolis dynamics as well as in the case of barriers , globally confirming the above discussed picture .
dynamics in the presence of barriers do not yield the same phenomenology as glauber and metropolis rate . in this case
, we have @xmath276 and @xmath277 . selecting @xmath278 , as in section  [ sec : glassy - dynamics-1 ] ,
we are led to @xmath279 ^ 2.\ ] ] as for the escape time @xmath113 , upon choosing @xmath45 size independent , the rest time @xmath122 will be diverging exponentially with @xmath51 at every temperature . by introducing the size dependence as in eq .
( [ eq : size_depen ] ) , instead , one can see that the @xmath135 integral converges to a constant for large @xmath51 so that one is left with @xmath280 we therefore obtain the same picture as in the case of traps , with an exponential increase of @xmath122 as @xmath46 increases , as confirmed by numerical simulations in fig .
[ simul_2temperatures ] .
inspired by analogies with systems governed by the arrhenius law , we have introduced a transition rate that takes into account the energy barriers between states . within the heterogeneous mean - field approximation , in which all variables depend only on the degree of the vertices , and choosing @xmath281 ,
the transition rate we have considered becomes @xmath282 where @xmath151 is a symmetric function of the degrees of the two nodes .
this model represents in essence an extension of the local trap model , where non - locality enters only in the form of symmetric energy gaps @xmath5 ( see fig .
[ landscape ] ) .
the steady state has exactly the same form as the ones discussed so far , which incidentally is the same as for the local trap model .
as shown in the previous sections , the presence of barriers affect transient relaxation phenomena , but not the steady state .
a different question is whether one can be more specific about the realistic functional form of the coarse - grained function @xmath151 . in the previous section
we have already introduced a definition of @xmath151 .
here we provide the rationale behind that choice .
numerical simulations of the energy - landscape network of lennard - jones clusters show that the average barrier to escape from state @xmath39 follow the power law @xmath283 , with @xmath284 @xcite . in our model ,
such average can be computed as @xmath285.\ ] ] for simplicity we focus on uncorrelated networks , as simulations indeed show weak degree correlations . under this assumption , the first term of the sum on the right - hand side of eq .
( [ avg_barrier ] ) will contribute as a logarithm of @xmath39 and the power - law behavior of @xmath286 is possible whenever @xmath287 , which leads to consider the form proposed in previous sections , @xmath288 where @xmath289 has the dimensions of an energy ( a discussion about the possible values of @xmath223 is given in section [ sec : glassy - dynamics-1 ] ) .
more complicated functional forms can also be proposed , for example accounting for barriers of different signs , as long as they retain the same power - law behavior of eq .
( [ sigma_linear ] ) in the large @xmath39 limit .
it is interesting to notice that @xmath283 implies that the average escape rate @xmath290 has the form of a stretched exponential @xmath291 , if we neglect the logarithmic correction .
in this paper , we have presented a simple mathematical framework for the description of the dynamics of glassy systems in terms of a random walk in a complex energy landscape .
we have shown how to incorporate into this picture the network representation of this landscape , put forward and studied by several authors @xcite , in order to go beyond simple mean - field models of random walks between traps that are all connected to each other .
while our previous work had focused on the case of a landscape consisting of traps connected by a network @xcite , we have here generalized our study to more involved and realistic transition rates between minima , including glauber or metropolis rates , and the possibility of energy barriers between minima .
we have shown how the interplay between the topology of the network of minima and the relationship between the energy and the degree of a minimum may determine a rich phenomenology , with the existence of two phases and of glassy dynamics at low temperature .
interestingly , the existence of these phases , and the transition temperature , do not depend on the network s degree correlations nor on the precise form of the transition rates , but other more detailed properties do . in the case of glauber and metropols dynamics
, the low temperature phase can be further divided into two regions with different scaling properties of the average trapping time as a function of the temperature .
overall , our results rationalize and link the empirical findings about correlations between the energy of the minima and their degree , and should stimulate further investigations on this issue .
our work has also interesting applications in terms of diffusion phenomena on complex networks , and shows that non trivial transition rates can lead to a very interesting phenomenology .
usual random walks lead to a higher probability for the random walker to be in a large degree node ( @xmath292 ) , with respect to the random choice of a node ( @xmath293 ) ; here , the models we have studied can lead to various stationary probabilities , such as for instance a uniform coverage which does not depend anymore on the degree .
interestingly , the biased random walks among traps that we have studied can even display a phase transition phenomenon , as either a temperature parameter or the network s properties are changed , with the possible presence of a glassy phase with slow dynamics .
p. m. , r.p .- s . , and
a. baronchelli acknowledge financial support from the spanish mec , under project fis2010 - 21781-c02 - 01 , and the junta de andaluca , under project no .
p09-fqm4682 .. r.p .- s
. acknowledges additional support through icrea academia , funded by the generalitat de catalunya .
a. baronchelli acknowledges support of spanish mci through the juan de la cierva program funded by the european social fund .


low mass x - ray binaries ( lmxbs ) are systems in which a low mass ( @xmath1 m@xmath2 ) star transfers matter to either a low magnetic field ( typically @xmath3 g ) neutron star or a black hole . in both cases
plasma flows down to a few stellar radii and produces observable properties in x - rays ( spectra and timing ) . over the last several years comprehensive studies of these properties
have provided crucial information on fundamental properties of compact objects , e.g. , the equation of state of neutron stars ( see van der klis 2000 for a review ) , or general relativistic effects near black hole horizons ( see tanaka & lewin 1995 , esin et al .
however , there exist several sources whose x - ray properties indicate they may also be lmxbs but whose x - ray emission was so faint that no instruments were sensitive enough to study them until the launch of _
this paper describes our _ chandra  _ observations of eight such objects listed as candidate lmxbs in @xcite ; the accurate locations we derived with _
chandra  _ will enable follow up observations in the infrared / optical .
below we briefly summarize prior observations for each of these objects .
table  [ tab : fluxsum ] provides a convenient comparison of the observed fluxes discussed below with our _ chandra  _ observations .
[ [ u-170840 ] ] 4u 170840 + + + + + + + + + +    is a persistent x - ray source that was observed by most of the early x - ray instruments , e.g. _ uhuru _ , _ oso-7 _ , _ ariel v _ , and _ heao-1 _
* and references therein ) .
rosat  _ source , 1rxs j17141224.8405034 , at r.a .
@xmath4 , decl .
= 405034.5 ( j2000 , 1-@xmath5 error radius = 10 ) , detected in the _ rosat  _ all - sky survey , lies in the error box for 4u 170840 @xcite .
since 1996 , 4u 170840 has been a persistent source in the _ rossi x - ray timing explorer ( rxte ) _ all - sky monitor ( asm ) .
its 2 - 10 kev flux ergs @xmath6 s@xmath7  ( 2 - 10 kev ) . ] slowly declined from about @xmath8 ergs @xmath6 s@xmath7  in 1996 - 1997 to a minimum of about @xmath9 ergs @xmath6 s@xmath7  in mid-1998 , followed by a slow rise to about @xmath10 ergs @xmath6 s@xmath7  in mid-2001 , and a second slow decline that continued in 2002 - 2003 . at the time of our _ chandra  _ observations , the 2 - 10 kev flux level in the _ rxte  _ asm was about @xmath11 ergs @xmath6 s@xmath7 .
@xcite reported the discovery of x - ray bursts in a _
bepposax  _ observation of 4u 1708408 in 1999 august . during this observation , the spectrum of the persistent flux (
before and after the bursts ) was well fitted with an absorbed power - law with a column density @xmath12 @xmath6 , a photon index of @xmath13 , and a gaussian iron line .
the unabsorbed 2 - 10 kev flux was @xmath14 ergs @xmath6 s@xmath7 .
in addition , @xcite also reported results for several _
rxte  _ observations in 1997 and 2000 along with an additional _ bepposax  _ observation in 2001 during which no bursts were seen ; however , steady persistent emission was observed . the energy spectra in these observations required more complicated models , an absorbed power law with photon index @xmath15 for an assumed @xmath16 @xmath6 , and a blackbody component with @xmath17 kev .
the unabsorbed 2 - 10 kev fluxes were @xmath18 ergs @xmath6 s@xmath7  and @xmath19 ergs @xmath6 s@xmath7  in the 2000 june 18 _ rxte  _ and 2001 august _
bepposax  _ observations , respectively .
[ [ s-1711339 ] ] 2s 1711339 + + + + + + + + + + +    is an x - ray burster originally detected with _
ariel v_@xcite and more precisely located with _
a radio source was detected in the _
sas-3  _ error circle at r.a .
= @xmath20 , decl .
= 340036(b1950 , 90% confidence error radius 2.2 ) in 1978 july @xcite .
rosat  _ source , rxs j171419.334023 @xcite , also lies within the _
sas-3  _ error circle , but just outside the radio source error circle .
more recently , @xcite reported on observations of 2s 1711339 with _ bepposax _ , the _
rxte  _ asm , and _
chandra_. the _
rxte  _ asm detected an x - ray outburst of 2s 1711339 from 1998 july until 1999 may at a 2 - 10 kev flux level of @xmath21 ergs @xmath6 s@xmath7 . during that time period , _
bepposax  _ detected 10 short x - ray bursts , with a persistent flux level ( before and after the bursts ) of @xmath22 ergs @xmath6 s@xmath7  ( 2 - 28 kev ) .
the persistent emission was fitted with a cutoff power - law with photon index 0.7 , high energy cutoff 2.8 kev , and an assumed @xmath23 of @xmath24 @xmath6 , derived from _
bepposax  _ narrow field instrument ( nfi ) observations on 2000 february 29 , when the 2 - 6 kev flux was @xmath25 ergs @xmath6 s@xmath7 . on 2000 march 22
, one burst was detected with _ bepposax _ , but the 3-@xmath5 upper limit on the persistent flux ( before and after the burst ) was @xmath26 ergs @xmath6 s@xmath7  ( 2 - 28 kev ) .
chandra  _ observation , also described in this paper , on 2000 june 9 , yielded a location of r.a .
= @xmath27 , decl .
= 340247  ( j2000 , 90% confidence error radius 1 ) , which was consistent with the wide field camera @xcite , _ ariel v _
@xcite , and _
rosat  _ @xcite positions .
chandra  _ position , however , lies 0.2  outside the 90% confidence error circle of the radio source .
[ [ slx-1735269 ] ] slx 1735269 + + + + + + + + + + + +    is an x - ray burster first reported in 1985 from spacelab-2 x - ray telescope ( xrt ) observations @xcite .
it was present in the _ einstein  _ slew survey 5 years earlier @xcite .
rosat  _ source , 1rxs j173817.0265940 , lies within the _ einstein  _ error circle @xcite .
_ granat  _ sigma observations showed it to be a persistent hard x - ray source @xcite .
type i x - ray bursts were discovered with the _ bepposax  _ wide field camera ( * ? ? ?
_ asca  _ observations showed the 0.6 - 10 kev spectrum to be well fitted with a power - law with photon index 2.15 and @xmath28 @xmath6 @xcite , which was consistent with the values derived from _
rosat  _ and _ granat  _ art - p observations @xcite . using _
rxte  _ observations from 1997 february - may , @xcite showed that the power spectrum was characterized by a strong band limited noise component which was approximately flat below a break frequency of 0.1 - 2.3 hz .
above the break frequency the power spectrum declines as a power law of index 0.9 . at the highest count rates ,
a broad bump was observed around 0.9 hz .
the strength of the noise ( 2 - 60 kev , integrated over 0.01 - 100 hz ) was @xmath29% and @xmath30% rms in the high and low count rate data , respectively .
fits to energy spectra yielded power law indices of 2.2 and 2.4 ( with an assumed @xmath23 of @xmath31 @xmath6 ) and 3 - 25 kev fluxes of @xmath32 ergs @xmath6 s@xmath7  and @xmath33 ergs @xmath6 s@xmath7 , respectively .
power spectra , from extended _ rxte  _ observations in 1997 october @xcite at a 1 - 20 kev flux level of @xmath34 ergs @xmath6 s@xmath7 , had a noise strength of 27.6% ( 2 - 40 kev , integrated over 0.005 - 300 hz ) and a break frequency of 0.08 - 0.15 hz that appeared to increase with intensity .
@xcite fitted this power spectrum with a four - lorentzian model and found a remarkable resemblance to observations of the low - luminosity x - ray burster 1e 17243045 .
slx 1735269 has been detected as a persistent source with the _ rxte  _ asm since 1996 at an average 2 - 10 kev flux level of about @xmath35 ergs @xmath6 s@xmath7 .
[ [ grs-1736297 ] ] grs 1736297 + + + + + + + + + + + +    was discovered with the _ granat  _ art - p telescope in september 1990 and was subsequently observed about a month later @xcite at a similar flux level .
its spectrum was characterized by a power law with photon index of 1.8 and a 3 - 12 kev flux of @xmath36 ergs @xmath6 s@xmath7 .
@xcite identified rx j1739.52942 with grs 1736297 , since it lies well within the 90% confidence 90  radius _ granat  _ error circle @xcite .
further , @xcite identified a be star in the _ rosat  _ error circle , suggesting that , although grs 1736297 has been previously classified as a candidate lmxb @xcite , it may instead be a transient be / x - ray binary system .
[ [ ks-1739304 ] ] ks 1739304 + + + + + + + + + + +    is a transient x - ray source discovered with _ mir_/kvant in 1989 at a 2 - 30 kev flux level of about 9 mcrab and located to within a 1.6  error circle @xcite .
subsequent observations with _
granat  _ art - p in autumn 1990 did not detect the source , with a 3-@xmath5 upper limit of 2 mcrab ( 8 - 20 kev ) @xcite .
very little is known about this object .
[ [ slx-1746331 ] ] slx 1746331 + + + + + + + + + + + +    was discovered with the spacelab-2 xrt in 1985 august @xcite .
it had a very soft spectrum , best fitted with a thermal bremsstrahlung with @xmath37 kev and a 2 - 10 kev flux of @xmath38 ergs @xmath6 s@xmath7 .
based on the very soft spectrum , @xcite suggested slx 1746331 as a potential black hole candidate .
the position of slx 1746331 fell between fields in the _
einstein  _ galactic plane survey of @xcite and it was not detected in the _ exosat  _ galactic plane survey @xcite . during the _ rosat  _ survey observation in 1990 september , 1rxs j174948.4331215 , which lies in the error circle of slx 1746331 , was detected in outburst , but it was not detected in subsequent _
rosat  _ high resolution imager ( hri ) observations in 1994 october @xcite .
hence the _ exosat  _ and _ rosat  _ non - detections suggest this object is a transient system .
[ [ e-1746.73224 ] ] 1e 1746.73224 + + + + + + + + + + + + + +    was discovered during the _ einstein  _ galactic plane survey @xcite , at a 2 - 10 kev flux level of @xmath39 ergs @xmath6 s@xmath7 .
spacelab-2 xrt observations in 1985 august gave an upper limit on the 2 - 10 kev flux of @xmath40 ergs @xmath6 s@xmath7  @xcite . during the _ rosat  _ all - sky survey , 1rxs j175003.8 - 322622 was detected in the _ einstein  _ error circle of 1e 1746.73224 @xcite .
very little is known about this object .
[ [ u-181212 ] ] 4u 181212 + + + + + + + + + +    was first detected with _
uhuru_@xcite and was later observed by several satellites , including _
oso-7 _ , _ ariel v _ , _
heao-1 _ , and _
exosat  _ with a typical 2 - 10 kev flux of @xmath41 ergs @xmath6 s@xmath7   ( * ? ? ? * and references therein ) .
4u 181212 was detected in the _ rosat  _ all - sky survey as a bright source denoted 1rxs j181506.1120545 @xcite .
x - ray bursts were first detected from this source in 1982 with hakucho @xcite and later with _
bepposax  _ @xcite .
4u 181212 is classified as an atoll source and shows band limited noise with a break frequency of 0.095 hz and a qpo at 0.85 hz @xcite .
4u 181212 has been detected as a persistent source with the _ rxte  _ asm since 1996 at an average 2 - 10 kev flux level of about @xmath42 ergs @xmath6 s@xmath7 .
recent observations in 2000 april with _ rxte  _ and _ bepposax  _ show 4u 181212 in a hard state with a hard tail extending above 100 kev .
the power spectrum was characterized by a @xmath43 hz qpo and three broad noise components , extending above @xmath44 hz @xcite .
each of the eight objects was observed using the advanced ccd imaging spectrometer ( acis ) on the _
chandra  _ x - ray observatory .
table  [ tab : obs ] lists the details of the observations . for some sources ,
data were collected in two different observing modes , timed exposure ( te ) mode and continuous clocking ( cc ) mode .
data obtained in te mode allow for two - dimensional imaging .
accurate spectroscopy of bright targets , however , is limited due to pulse pile - up .
timing studies are limited by the total number of photons collected and the 3.24 s time resolution . in cc mode , the amount of pile - up is negligible owing to the 2.85 ms time resolution , allowing for accurate spectroscopy and timing . unfortunately only one of the sources , slx 1735269 , was detected in cc mode .
the best known location for each object was positioned on the nominal target position of acis - s3 , a back illuminated ccd on the spectroscopic array ( acis - s ) with good charge transfer efficiency and spectral resolution , with the exception of the cc mode observation of slx 1735269 which fell on acis - s2 , a front illuminated ccd on the spectroscopic array ( acis - s ) .
standard processing was performed by the _
chandra  _ data center .
the data were filtered to exclude events with grades 1 , 5 , 7 , hot pixels , bad columns , and events on ccd node boundaries .
analysis in this paper was done using _
chandra  _ interactive analysis of observations ( ciao ) version 2.2.1 and _
chandra  _ calibration database ( caldb ) version 2.10 . for piled - up sources , the tool _
acis_detect_afterglow _ in the standard processing pipeline can reject valid source photons in addition to afterglow photons resulting from residual charge from cosmic ray events in the ccd pixels .
we examined the spatial distribution of afterglow events for each observation using methods described in the ciao documentation . for te mode observations of 2s 1711339 and 1e 1746.73224
we found that @xmath45% of events consistent with the point source had been rejected ; hence we reprocessed the data to retain the afterglow flagged events .
for all other observations , @xmath46% of events were rejected , so the standard processing was retained .      using te mode data
, we extracted source locations for each of the detected sources . these locations are listed in table  [ tab : loc ] .
sources were located by one of two methods : ( 1 ) 2s 1711 - 339 and 1e 1746.73224 were located using the ciao tool _
wavdetect_. the latter source was observed twice in te mode and locations from both observations were consistent .
( 2 ) 4u 170840 , slx 1735269 , and 4u 181212 were brighter and suffered from considerable pile - up of photons in the image core , resulting in a source that looked ring shaped with a hole in the center .
the ciao detect tools were inadequate to provide a good location of the centroid , therefore we modeled the data using a gaussian function multiplied by a hyperbolic tangent in radius , scaled to approach zero at r = 0.0 @xcite .
the uncertainty of these locations is limited by systematic effects to a circle with @xmath47 radius @xcite .
no other known sources were found in the images , hence we were unable to use astrometry to further improve the locations .
three objects , grs 1736297 , slx 1746331 , and ks 1739304 went undetected in all observations . to derive upper limits for these sources , we used the ciao tool _
dmextract _ to extract count rates for a 6  radius source region and a background annulus ( inner radius = 6 , outer radius = 30 ) centered on the best known position .
the 99% confidence upper limits were count rates of @xmath48 cts s@xmath7 , @xmath49 cts s@xmath7 , and @xmath50 cts s@xmath7 ( total counts = 6.0 cts , 5.8 cts , and 4.6 cts ) for grs 1736297 , slx 1746331 , and ks 1739304 , respectively . using webpimms , we estimated an upper limit 1 - 10 kev absorbed flux of @xmath51 ergs @xmath6 s@xmath7   for grs 1736297 , assuming @xmath52 @xmath6 , derived using colden , a web tool based on @xcite , and a power law photon index = 1.8 @xcite . similarly , for slx 1746331 , we estimated an upper limit 1 - 10 kev absorbed flux of @xmath53 ergs @xmath6 s@xmath7 , assuming @xmath54 @xmath6 , using colden , and a thermal bremsstrahlung spectrum with kt = 1.5 kev @xcite . for ks 1739304 , we estimated an upper limit 1 - 10 kev absorbed flux of @xmath55 ergs @xmath6 s@xmath7 , assuming @xmath56 @xmath6 , using colden , and a power law photon index = 1.5 . a fourth object , 2s 1711339 , that was detected in a te mode observation in 2000 , went undetected in the 2002 cc mode observation . to derive an upper limit for it , we again used dmextract , with a 6  radius source region and a background annulus ( inner radius = 6 , outer radius = 30 ) centered on the _
chandra  _ position ( shifted to correspond to the 1-dimensional cc mode image ) from the 2000 observation , to obtain 99% confidence upper limits of @xmath57 cts s@xmath7 or a total 30.3 counts for the observation . using webpimms with our spectral parameters from the 2000 observation table  [ tab : spec
] , we estimate a 1 - 10 kev upper limit absorbed flux of @xmath58 ergs @xmath6 s@xmath7 , a factor of more than 100 fainter than in the 2000 observation . in the 5 may 2002 observation of ks 1739304 , an object was present at r.a
. = @xmath59 , decl.= @xmath60 , approximately @xmath61  from the expected location of ks 1739304 .
we measured a total of approximately 45 counts from this object in the 745 s observation .
given its large separation from ks 1739304 , and the fact that it is well outside the 1.6  radius _ mir_/kvant error circle @xcite , the two objects are most likely unrelated .
energy spectra were difficult to study for most of these sources because in all but one of the observations where a source was detected , we had only imaging ( te ) mode data and pile - up was a problem .
for these sources we used two approaches to extract and fit the energy spectra : ( 1 ) direct extraction from the te mode data , fitting the data using the pile - up model and ( 2 ) extraction of trailed image spectra . to directly extract spectra and associated response files for each source , we used a source region consisting of a 6  radius circle and a background region consisting of an annulus with inner and outer radii of 6  and 30 , respectively , in the ciao tool _ psextract_. these spectra were then fitted using xspec 11.2 . for all sources , we fitted data for energies from 0.3 - 10.0 kev using an absorbed power law model with pile - up ( phabs*powerlaw*pile - up ) . for most sources we found that the psf fraction treated for pile - up ( parameter 5 ) in the pile - up model , needed to be reduced from its default value of 0.95 to 0.85 to allow the grade morphing parameter ( parameter 4 ) to not be pegged at 1.0 and to allow the model to fit the data at higher energies .
for some sources , the pile - up model had no effect on the spectral fits , even though the source was obviously piled - up ; hence we used trailed images to characterize those spectra .
results of our fits with the pile - up model are listed in table  [ tab : spec ] .
figure  [ fig : pile - up ] shows examples of spectra before and after the pile - up model is applied .
trailed images were extracted using the ciao tool _
acisreadcorr_. during each row transfer , the detector is exposed to the source for 40@xmath62s , but the photons are recorded in the wrong y position , resulting in a streak along the x - axis .
we calculated the total transfer streak exposure time using the following expression @xmath63 where @xmath64s , the time to transfer charge from one row to the next ; 1024 is the number of pixels in a row ; @xmath65 is the subarray fraction , @xmath66 is the number of pixels excluded near the source ; @xmath67 is the total exposure time in seconds ; and @xmath68 s , the time for one readout of the subarray .
spectra were extracted using the ciao tool _ psextract _ with the source region defined as two boxes along the y - axis , each 6 pixels wide in the x - direction .
the two boxes included the entire y - axis , except a region @xmath69 pixels from the y - coordinate of the piled - up source .
the exposure and backscal keywords in the spectra were corrected for using the results of the ciao tool _
acisreadcorr _ and the calculation above .
the resulting exposure times and total counts measured from the source in each trailed image are listed in table  [ tab : trspec ] .    for five of the observations ,
trailed image spectra could be extracted .
we grouped the source spectra into bins containing 15 counts .
then these spectra were fitted in xspec 11.2 with an absorbed power - law model ( phabs*powerlaw ) . for two of the sources , 2s 1711339 and 1e 1746.73224
, very few counts remained in the trailed image spectra .
the results of fitting trailed image spectra are listed in table  [ tab : spec ] and examples are shown in figure  [ fig : trailed ] .
for one of the sources , slx 1735269 , energy spectra were extracted from cc mode data taken on 2000 april 4 . a source region consisting of a 6   radius circle and a background region consisting of an annulus with inner and outer radii of 6  and 30 , respectively , was used in the ciao tool _
psextract _ to extract source and background spectra for the cc mode data and to create associated detector response files .
we then grouped the data in the source spectrum so that each bin contained 25 counts .
this observation fell on the acis - s2 chip , a front illuminated chip , so the low - energy response was reduced . in xspec 11.2 @xcite , we fitted data for energies of 1.0 - 1.8 kev and 2.2 - 9.0 kev with an absorbed power - law model ( phabs*powerlaw ) .
fit results are listed in table  [ tab : spec ] and are shown in figure  [ fig : ccspec ] .
the cc mode observations of slx 1735269 resulted in the largest number of source counts of any of our sources .
hence we used these observations to investigate if more complicated spectral models are warranted .
for example , @xcite reported an iron line with energy @xmath70 kev , fixed width 0.9 kev , and equivalent width 344 ev in addition to an absorbed power - law in observations of 4u 1708408 . to search for evidence of a similar line in slx 1735269
, we included a gaussian with a fixed centroid at 6.5 kev and fixed width of 0.9 kev .
our fit resulted in a 3-sigma upper limit on the equivalent width of 975 ev , indicating that we had insufficient statistics to detect such a line .
next , since a soft thermal component has been observed in several lmxbs , we fitted a blackbody in addition to the absorbed power - law .
the fit gave @xmath71 with 253 degrees of freedom .
an f - test comparing this fit to the absorbed power law had a value of 6.02 and a probability of @xmath72 , suggesting that a softer component may be present .
the resulting fit parameters were : @xmath73 @xmath6 , blackbody temperature @xmath74 kev , blackbody normalization @xmath75 , power law index @xmath76 , and power law normalization @xmath77 photons kev@xmath7 @xmath6 s@xmath7 at 1 kev .
since this observation with the largest number of source counts resulted in only a suggestion of a softer component , we concluded that we had insufficient statistics to distinguish spectral models in the rest of our observations , which had significantly fewer source counts .
rxte  _ observations of slx 1735269 with 3 - 25 kev fluxes @xmath78 ergs @xmath6 s@xmath7  show band limited noise levels with fractional amplitudes of 24 - 28% and break frequencies of 0.1 - 0.2 hz @xcite and a qpo near 0.9 hz with a fractional rms of @xmath79% and fwhm = 0.3 - 0.8 hz .
@xcite also reported a very different power spectrum from  low count rate " 1997 _ rxte  _ observations of slx 1735269 with 3 - 25 kev fluxes @xmath80 ergs @xmath6 s@xmath7 , that had a break frequency of 2.3 hz and a fractional rms of 17% .
@xcite fitted _ bepposax  _ power spectra for 4u 181212 with a 4 lorentzian model after @xcite .
the low - frequency lorentzian , fitting the low - frequency end of the band limited noise , had a fractional amplitude of 9.8% rms and a break frequency of 0.15 hz . a low - frequency qpo was also detected with _
bepposax  _ at 0.73 hz with a fractional amplitude of 3.9% rms .
we computed the sensitivity of our chandra observations to features in the power spectrum using @xmath81 where @xmath82 is the fractional root - mean - squared ( rms ) of the signal ; @xmath83 is the significance of the expected detection ; @xmath84 is the total observed count rate ; @xmath85 is the source count rate ; @xmath86 is the full - width half - max ( fwhm ) of the signal , which is approximately equal to twice the break frequency obtained with a broken power law model for band limited noise ( see belloni , psaltis , & van der klis 2002 ) ; and @xmath87 is the exposure time . for te - mode observations of slx 1735269 , 4u 1708408 , 2s 1711339 , 4u 181212 , and
the 2000 aug 30 observation of 1e 1746.73224 , @xmath84 = 0.2 - 0.7 cts s@xmath7 , @xmath85 = 0.2 - 0.7 cts s@xmath7 , and @xmath87 = 950 - 6800 s , resulting in 3-@xmath5 sensitivities of 30 - 50% rms for features with @xmath88 hz .
the 2002 july 16 observation of 1e 1746.73224 was longer , with @xmath87 = 8500 s , and had slightly higher count rates of @xmath89 cts s@xmath7 and @xmath90 cts s@xmath7 , resulting in a 3-@xmath5 sensitivity of 23% rms .
our _ chandra  _ te mode observations were not sensitive enough to detect broadband features at previously detected levels .
however , the cc mode observation of slx 1735269 , with @xmath91 cts s@xmath7 , @xmath92 cts s@xmath7 , and @xmath93 s , had a much better 3-@xmath5 sensitivity of 9.3% rms for @xmath94 hz .
hence band limited noise at levels observed in the brighter state of slx 1735269 @xcite should have been detectable in the cc mode observations ; however , the observed qpo was not detectable .
we estimated a 3-@xmath5 sensitivity of 18% rms for @xmath95 hz , corresponding to the fainter state observed by @xcite , for our cc mode observation of slx 1735269 ; hence we were unable to confirm this state with our _ chandra  _ observations .
for the cc mode observation of slx 1735269 , we created a light curve using ciao tools that account for the fact that event times recorded in cc mode data were the times the events were read , not the times when the charge was deposited on the detector @xcite .
we generated a power spectrum for the energy range 0.3 - 10 kev from the 2.85 ms cc mode light curve .
no obvious variability was present in either the light curve or the power spectrum . integrating from 0.01 - 1 hz
, we estimated a noise strength of @xmath96% rms for this observation ; since this noise level was near our detection threshold , we can not easily confirm that this was entirely source related .
however , the lack of strong band limited noise at levels @xmath97% rms , combined with our 1 - 10 kev flux measurement of @xmath98 ergs @xmath6 s@xmath7  suggests that slx 1735269 was most likely in a state similar to the  lowest count rate " state reported in @xcite .
we obtained images of slx 1735269 , 1e 1746.73224 , and 4u 181212 using the united kingdom infrared telescope ( ukirt ) to search for near - infrared counterparts . for slx 1735269 ,
the seeing was approximately 0.7 , the airmass was 1.46 , the upper limit on the presence of a star in the _
chandra  _ error circle was @xmath99 .
the astrometric solution was performed using five stars from the 2mass in j catalog .
the fit was good with an rms of 0.26 .
the astrometric accuracy of the 2mass catalogue is @xmath100 . for 1e 1746.73224
the seeing was approximately 0.6 , the airmass was 1.64 , and the upper limit for a source in the _
chandra  _ error circle was @xmath101 .
the astrometric solution was performed using eight stars from the 2mass in j catalog .
the fit was good with an rms of 0.23 .
there was no 2mass image available for 4u 181212 ; hence , the astrometric solution of the ukirt image of the field containing 4u 181212 ( seeing @xmath102 , airmass 1.18 ) was obtained by identifying four usno - a1.0 stars in an optical digital sky survey ( dss ) image with a near infrared star . however , with an rms of 0.8  the astrometric solution was poor .
a star appears to be present near the error circle , but its j band magnitude is difficult to assess since we have no j band magnitudes for other stars in the field for comparison .
chandra  _ we observed eight faint , little studied likely lmxbs .
of the eight systems , we detected five in at least one observation : 4u 1708408 , 2s 1711339 , slx 1735269 , 1e 1746.73224 , and 4u 181212 .
chandra  _ observations reported in this paper resulted in precise locations for all five of these objects , allowing for future deep optical and infrared observations .
ukirt images of slx 1735269 , 1e 1746.73224 , and 4u 181212 did not reveal counterparts in the j band .
energy spectra for these five systems were generally consistent with hard power laws ( with photon index 1.5 - 2 ) typical for faint lmxbs ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) .
absorbed fluxes ( 1 - 10 kev ) for the detected sources ranged from @xmath103 ergs @xmath6 s@xmath7 .
the relatively narrow range of power law indices and observed fluxes suggest that all five of the detected systems may be in a similar state .
energy spectra and fluxes for 4u 1708408 , slx 1735269 , and 4u
181212 were consistent with previous measurements in similar energy bands , supporting earlier findings that these systems are persistent .
we observed slx 1735269 at a 1 - 10 kev flux level of @xmath104 ergs @xmath6 s@xmath7  in both observations .
this combined with the lack of detection of strong band limited noise in the cc mode observation suggests that slx 1735269 was in a similar state to the  low count rate " state reported in @xcite . the energy spectrum and flux measured using the pile - up model with _
chandra  _ on 2000 june 9 for 2s 1711339 were consistent with that measured with the _ bepposax  _ nfi on 2000 feb 29 and with _
bepposax  _ wfc upper limits on 2000 march 22 , suggesting that we were observing either an extended tail of the bright 1998 - 1999 outburst observed with _
bepposax  _ and _ rxte  _ or a faint outburst .
an estimate of the 2 - 6 kev flux reported in @xcite from the 2000 june 9 _ chandra  _ observation was a factor of @xmath105 fainter than the flux we measured with _
chandra  _ for the same observation .
we believe that @xcite did not properly account for the effects of pulse pile - up . in our second _
chandra  _ observation of 2s 1711339 , in 2002 march , it had faded below detectability in cc mode .
1e 1746.73224 , the least studied of our five detected sources , and the only one of the five not to have previously shown x - ray bursts , was observed with _
chandra  _ at absorbed 1 - 10 kev fluxes of @xmath106 ergs @xmath6 s@xmath7  and @xmath107 ergs @xmath6 s@xmath7  in 2000 and 2002 , respectively . in the _ rosat  _ all - sky survey @xcite , 1e 1746.73224 was measured at a count rate of 0.08 counts s@xmath7 in the position sensitive proportional counter ( pspc ) .
this corresponds to an absorbed 1 - 10 kev flux of @xmath108 ergs @xmath6 s@xmath7 , if the spectral shape we measured with _ chandra  _ is assumed .
these fluxes are a factor of @xmath96 higher than that observed in the _ einstein  _ galactic plane survey @xcite .
the similarity of the _
chandra  _ and _ rosat  _ fluxes suggests this source is persistent , and the _
einstein  _ measurement suggests that it , like most persistent sources , is also variable .
the three undetected systems , grs 1736297 , ks 1739304 , and slx 1746331 , had shown previous transient behavior .
these non - detections were not due to positional errors .
upper limits on their fluxes were @xmath109 ergs @xmath6 s@xmath7 , two to three orders of magnitude fainter than our faintest detection .
if these sources were located at the galactic center , at a distance of 8 kpc , then these fluxes would correspond to luminosities of @xmath110 ergs s@xmath7 .
the united kingdom infrared telescope is operated by the joint astronomy centre on behalf of the u.k .
particle physics and astronomy research council , and some of the data reported here were obtained as part of the ukirt service programme .
ck and sp are party supported by sao grants go0 - 1054a and go2 - 3046b .
whgl is grateful for support from nasa .
aldcroft , t.l . ,
karovska , m. , cresitello - ditmar , m.l . ,
cameron , r.a .
, & markevitch , m.l .
2000 , proc .
spie , 4012 , 650 arnaud , k.a .
1996 , astronomical data analysis software and systems v , eds . g. jacoby & j. barnes , asp conf .
series , 101 , 17 barret , d. et al .
2000 , apj 533 , 329 barret , d. , mcclintock , j.e . , &
grindlay , j.e .
1996 , apj , 473 , 963 barret , d. , olive , j.f . , &
oosterbroek , t. 2003 , in press , astro - ph/0301130 bazzano , a. et al .
1997 , iau circ . , 6668 , 2 belloni , t. , psaltis , d. , & van der klis , m. 2002 , apj , 572 , 392 carpenter , g.f . , eyles , c.j . ,
skinner , g. k. , wilson , a. m. , & willmore , a.p .
1977 , mnras , 179 , 297 cherepashchuk , a.m. et al .
1994 , a&a , 289 , 419 cocchi , m. et al .
2000 , , 357 , 527 cornelisse , r. et al .
2002 , , 392 , 885 david , p. et al .
1997 , , 322 , 229 dickey , j.m & lockman , f.j .
1990 , ara&a , 28 , 215 elvis , m. et al .
1992 , apjs , 80 , 257 esin , a.a .
2001 , apj , 555 , 483 forman , w. , jones , c. , tananbaum , h. 1976 , apj , 206 , l29 forman , w. et al .
1978 , apjs , 38 , 357 goldwurm , a. et al .
1996 , , 310 , 857 grebenev , s.a . ,
pavilinsky , m.n . , & sunyaev , r.a .
1996 , in the 2nd integral workshop : the transparent universe , esa sp-382 ( noordwijk : esa ) , 183 greenhill , j.g .
thomas , r.m . & duldig , m.l .
1979 , nature , 279,620 hertz , p. & grindlay , j.e . 1984 , apj , 278 , 137 hulleman , f. , tennant , a.f . ,
van kerkwijk , m.h . ,
kulkarni , s.r .
, kouveliotou , c. & patel , s.k .
2001 , apj , 563 , l49 migliari , s. et al .
2003 , mnras , in press , astro - ph/0303136 motch , c. et al .
1998 , a&as , 132 , 341 murakami , t. et al . 1983 , pasj , 35 , 531 pavilinsky , m.n . ,
grebenev , s.a . , & sunyaev , r.a .
1992 , sval , 18 , 88 pavilinsky , m.n . , grebenev , s.a . , & sunyaev , r.a .
1994 , apj , 425 , 110 skinner , g. k. et al .
1987 , nature , 330 , 544 skinner , g.k . ,
foster , a.j . ,
willmore , a.p . , & eyles , c.j .
1990 , mnras , 243 , 72 sunyaev , r. et al .
1991 , adv .
space res .
, 11(8 ) , 177 tanaka , y. & lewin , w.h.g 1995 , in x - ray binaries , eds .
lewin , j. van paradijs , & e.p.j .
van den heuvel , ( new york : cambridge university press ) , 126 van der klis , m. 1994 , apjs , 92 , 511 van der klis , m. 2000 , ann .
astrophys . , 38 , 717 van paradijs , j. 1995 , in x - ray binaries , eds .
lewin , j. van paradijs , & e.p.j .
van den heuvel , ( new york : cambridge university press ) , 536 van paradijs , j. & van der klis , m. 1994 , , 281 , l17 voges , w. et al .
1999 , , 349 , 389 warwick , r.s . ,
norton , a.j . ,
turner , m.j.l . ,
watson , m.g .
& willingdale , r. 1988 , mnras , 232 , 551 wijnands , r. & van der klis , m. 1999a , apj , 514 , 939 wijnands , r. & van der klis , m. 1999b , , 345 , l35 zavlin , v.e . , pavlov , g.g , sanwal , d. , & trmper , j. 2000 , apj , 540 , l25    lcccc 4u 170840 & 1996 - 1997 & _ rxte  _ asm & 2 - 10 kev & @xmath8 + & mid-1998 & _ rxte  _ asm & 2 - 10 kev & @xmath9 + & 1999 august & _ bepposax  _ & 2 - 10 kev & @xmath14 + & 2000 may 15 & _ chandra  _ & 1 - 10 kev & @xmath111 + & 2000 jun 18 & _ rxte  _ pca & 2 - 10 kev & @xmath18 + & mid-2001 & _ rxte  _ asm & 2 - 10 kev & @xmath10 + & 2001 august & _ bepposax  _ & 2 - 10 kev & @xmath112 + 2s 1711339 & 1998 july - 1999 may & _ rxte  _ asm & 2 - 10 kev & @xmath113 + & 1998 july - 1999 may & _ bepposax  _ wfc & 2 - 28 kev & @xmath22 + & 2000 february 29 & _ bepposax  _ nfi & 2 - 6 kev & @xmath25 + & 2000 march 22 & _ bepposax  _ wfc & 2 - 28 kev & @xmath114 + & 2000 june 9 & _ chandra  _ & 1 - 10 kev & @xmath115 + & 2002 march 12 & _ chandra  _ & 1 - 10 kev & @xmath116 + slx 1735269 & 1996 - 2003 & _ rxte  _ asm & 2 - 10 kev & @xmath35 + & 1997 february - may & _ rxte  _ pca & 3 - 25 kev & @xmath117 + & 1997 october & _ rxte  _ pca & 1 - 20 kev & @xmath118 + & 2000 april 4 & _ chandra  _ & 1 - 10 kev & @xmath98 + & 2000 may 23 & _ chandra  _ & 1 - 10 kev & @xmath119 + grs 1736297 & 1990 september - october & _ granat  _
art - p & 3 - 12 kev & @xmath120 + & 2000 may 31 & _ chandra  _ & 1 - 10 kev & @xmath121 + ks 1739304 & 1989 & _ mir_/kvant & 2 - 30 kev & @xmath104 + & 1990 & _ granat  _ art - p & 8 - 20 kev & @xmath122 + & 2002 may 5 & _ chandra  _ & 1 - 10 kev & @xmath116 + slx 1746331 & 1985 august & spacelab-2 xrt & 2 - 10 kev & @xmath123 + & 2000 june 9 & _ chandra  _ & 1 - 10 kev & @xmath124 + 1e 1746.73224 & 1978 - 1981 & _
einstein  _ & 2 - 10 kev & @xmath39 + & 1985 august & spacelab-2 xrt & 2 - 10 kev & @xmath125 + & 1990 - 1991 & _ rosat  _ pspc & 1 - 10 kev & @xmath108 + & 2000 august 30 & _
chandra  _ & 1 - 10 kev & @xmath106 + & 2002 july 15 - 16 & _ chandra  _ & 1 - 10 kev & @xmath107 + 4u 181212 & & _ oso-7 _ , _ ariel v _ , _
heao-1 _ , _
exosat  _ & 2 - 10 kev & @xmath126 + & 1996 - 2003 & _ rxte  _ asm & 2 - 10 kev & @xmath42 + & 2000 june 14 & _ chandra  _ & 1 - 10 kev & @xmath127 + [ tab : fluxsum ]    lclcccc 4u 170840 & 661 & 2000 may 15 & te & 0.841 s & 1.2 & y + 2s 1711339 & 662 & 2000 jun 9 & te & 0.841 s & 1.0 & y + & 2695 & 2002 mar 12 & cc & 2.85 ms & 4.7 &
n + slx 1735269 & 664 & 2000 apr 4 & cc & 2.85 ms & 1.4 & y + & 663 & 2000 may 23 & te & 0.841 s & 1.7 & y + grs 1736297 & 665 & 2000 may 31 & te & 0.841 s & 1.3 & n + & 666 & 2000 jul 7 & cc & 2.85 ms & 12.2 &
n + ks 1739304 & 2696 & 2002 apr 27 & cc & 2.85 ms & 3.4 & n + & 2697 & 2002 may 5 & te & 3.24 s & 0.8 & n + slx 1746331 & 667 & 2000 jun 9 & te & 0.841 s & 1.7 & n + & 668 & 2000 jul 24 & cc & 2.85 ms & 4.1 & n + 1e 1746.73224 & 669 & 2000 aug 30 & te & 3.24 s & 6.8 & y + & 2698 & 2002 jul 15 - 16 & te & 0.841 s & 8.5 & y + 4u 181212 & 670 & 2000 jun 14 & te & 3.24 s & 1.0 & y + [ tab : obs ]    lll 4u 170840 & @xmath128 & @xmath129 + 2s 1711339 & @xmath130 & @xmath131 + slx 1735269 & @xmath132 & @xmath133 + 1e 1746.73224 & @xmath134 & @xmath135 + & @xmath136 & @xmath137 + 4u 181212 & @xmath138 & @xmath139 + [ tab : loc ]    lllllllc 4u 170840 & 2000 may 15 & tr & @xmath140 & @xmath141 & & @xmath142 & 87 + 2s 1711339 & 2000 jun 9 & te & @xmath143 & @xmath144 & @xmath145 & @xmath146 & 4.4 + & & tr & @xmath147 & @xmath148 & & @xmath149 & 11 + & 2002 mar 12 & cc & 1.4 & 1.9 & & & @xmath150 + slx 1735269 & 2000 apr 4 & cc & @xmath151 & @xmath152 & & @xmath153 & 19 + & 2000 may 23 & tr & @xmath154 & @xmath155 & & @xmath156 & 21 + grs 1736297 & 2000 may 31 & te & 1.1 & 1.8 & & & @xmath157 + ks 1739304 & 2002 may 5 & te & 1.4 & 1.5 & & & @xmath150 + slx 1746331 & 2000 jun 9 & te & 0.4 & 1.5 & & & @xmath158 + 1e 1746.73224 & 2000 aug 30 & te & @xmath159 & @xmath160 & @xmath161 & @xmath162 & 2.1 + & 2002 jul 15 - 16 & te & @xmath163 & @xmath164 & @xmath165 & @xmath166 & 3.2 + & & tr & @xmath167 & @xmath168 & & @xmath169 & 3.3 + 4u 181212 & 2000 jun 14 & tr & @xmath170 & @xmath171 & & @xmath172 & 44 + [ tab : spec ]    llcc 4u 170840 & 2000 may 15 & 11.2 & 412 + 2s 1711339 & 2000 jun 9 & 8.84 & 86 +
slx 1735269 & 2000 may 23 & 16.0 & 228 + 1e 1746.73224 & 2002 jul 15 - 16 & 79.4 & 117 + 4u 181212 & 2000 jun 14 & 11.9 & 312 + [ tab : trspec ]