differential analysis is based on studying the time evolution of the distance between trajectories emanating from different initial conditions .
a dynamical system is called contractive if any two trajectories converge to one other at an exponential rate .
this implies many desirable properties including convergence to a unique attractor ( if it exists ) , and entrainment to periodic excitations  @xcite .
contraction theory proved to be a powerful tool for analyzing nonlinear dynamical systems , with applications in control theory  @xcite , observer design  @xcite , synchronization of coupled oscillators  @xcite , and more .
recent extensions include : the notion of partial contraction  @xcite , analyzing networks of interacting agents using contraction theory  @xcite , a lyapunov - like characterization of incremental stability  @xcite , and a lasalle - type principle for contractive systems  @xcite .
a contractive system with added diffusion terms or random noise still satisfies certain asymptotic properties  @xcite . in this respect , contraction is a robust property .    in this paper , we introduce three forms of generalized contractive systems  ( gcss ) .
these are motivated by requiring to take place only after arbitrarily small transients in time and/or amplitude .
indeed , contraction is usually used to prove _ asymptotic _ properties , and thus allowing ( arbitrarily small ) transients seems reasonable . in some cases as we change the parameters in a contractive system it becomes a  gcs just before it looses contractivity . in this respect ,
a gcs is the analogue of marginal stability in lyapunov stability theory .
we provide several sufficient conditions for a system to be a  gcs .
these conditions are checkable , and we demonstrate their usefulness using examples of systems that are _ not _ contractive with respect to any norm , yet are gcss .
we begin with a brief review of some ideas from contraction theory . for more details , including the historic development of contraction theory , and the relation to other notions ,
see e.g.  @xcite .
consider the time - varying system [ eq : fdyn ] = f(t , x ) , with the state @xmath0 evolving on a convex set  @xmath1 .
we assume that  @xmath2 is differentiable with respect to  @xmath0 , and that both  @xmath2 and  @xmath3 are continuous in  @xmath4 .
let  @xmath5 denote the solution of   at time  @xmath6 with  @xmath7 ( for the sake of simplicity , we assume from here on that @xmath5 exists and is unique for all  @xmath8 and all  @xmath9 ) .
we say that   is _ contractive _ on  @xmath10 with respect to a norm  @xmath11 if there exists  @xmath12 such that [ eq : contdef ]    @xmath13 and all  @xmath14 .
in other words , any two trajectories contract to one another at an exponential rate .
this implies in particular that the initial condition is `` quickly forgotten '' .
recall that a vector norm  @xmath15 induces a matrix measure @xmath16 defined by @xmath17 where  @xmath18 is the matrix norm induced by  @xmath19 .
a standard approach for proving   is based on bounding some matrix measure of the jacobian  @xmath20 . indeed ,
it is well - known ( see , e.g.  @xcite ) that if there exist a vector norm  @xmath19 and  @xmath12 such that the induced matrix measure  @xmath21 satisfies [ eq : jtc ] ( j(t , x ) ) -c , for all @xmath13 and all  @xmath22 then   holds .
( this is in fact a particular case of using a lyapunov - finsler function to prove contraction  @xcite . )
it is well - known  ( * ? ? ?
3 ) that the matrix measure induced by the  @xmath23 vector norm is [ eq : muqdef ] _
1(a)=\{c_1(a ) ,  , c_n(a ) } , where [ eq : ccstac ] c_j(a):=a_jj+ _ |a_ij| , i.e. , the sum of the entries in column  @xmath24 of  @xmath25 , with non diagonal elements replaced by their absolute values .
the matrix measure induced by the  @xmath26 norm is [ eq : mat_meas_inf ] _ ( a)=\{d_1(a ) ,  ,
d_n(a ) } , where [ eq : dstac ] d_j(a):=a_jj+ _ |a_ji| , i.e. , the sum of the entries in row  @xmath24 of  @xmath25 , with non diagonal elements replaced by their absolute values .
often it is useful to work with scaled norms .
let  @xmath27 be some vector norm , and let  @xmath28 denote its induced matrix measure .
if  @xmath29 is an invertible matrix , and @xmath30 is the vector norm defined by @xmath31 then the induced matrix measure is @xmath32    one important implication of contraction is _ entrainment _ to a periodic excitation . recall that  @xmath33 is called  _
@xmath34-periodic _ if @xmath35 for all @xmath36 and all  @xmath22 .
it is well - known that if   is contractive and  @xmath37 is  @xmath34-periodic then for any  @xmath38 there exists a unique periodic solution @xmath39 of  , of period  @xmath34 , and every trajectory converges to  @xmath40 .
entrainment is important in various applications ranging from biological systems  @xcite to the stability of a power grid  @xcite .
note that for the particular case where  @xmath37 is time - invariant , this implies that if  @xmath10 contains an equilibrium point  @xmath41 then it is unique and all trajectories converge to  @xmath41 .
the remainder of this paper is organized as follows . .
section  [ sec : main ] details sufficient conditions for their existence , and describes their implications .
the proofs of all the results are detailed in section  [ sec : proofs ] .
we begin by defining three generalizations of  .    [
def : qcont ] the time - varying system   is said to be :    * _ contractive after a small overshoot and short transient _  ( sost ) on  @xmath10 w.r.t . a norm  @xmath11 if for each  @xmath42 and each  @xmath43 there exists  @xmath44 such that @xmath45 for all @xmath46 and all @xmath14 . * _ contractive after a small overshoot _ ( so ) on  @xmath10 w.r.t . a norm  @xmath11 if for each  @xmath42 there exists  @xmath47 such that @xmath48 for all @xmath46 and all @xmath14 . * _ contractive after a short transient _ ( st ) on  @xmath10 w.r.t . a norm  @xmath11 if for each  @xmath43 there exists  @xmath49 such that @xmath50 for all @xmath46 and all @xmath14 .
the definition of sost is motivated by requiring contraction at an exponential rate , but only after an ( arbitrarily small ) time  @xmath51 , and with an ( arbitrarily small ) overshoot  @xmath52 . however , as we will see below when the convergence rate  @xmath53 may depend on  @xmath54 a somewhat richer behavior may occur .
the definition of  so is similar to that of  sost , yet now the convergence rate  @xmath55 depends only on  @xmath54 , and there is no time transient  @xmath51 ( i.e. ,  @xmath56 ) . in other words ,
so is a uniform ( in  @xmath51 ) version of  sost .
the third definition ,  st , allows the contraction to `` kick in '' only after a time transient of length  @xmath51 .
it is clear that every contractive system is  sost , so , and  st .
thus , all these notions are generalizations of contraction . also , both so and st imply  sost and , as we will see below , under a mild technical condition on   and sost are equivalent .
figure  [ fig : graphbn2 ] on p.   summarizes the relations between these gcss ( as well as other notions defined below ) .
the motivation for these definitions stems from the fact that important applications of contraction are in proving _
asymptotic _ properties .
for example , proving that an equilibrium point is globally attracting or that the state - variables entrain to a periodic excitation .
these properties describe what happens as  @xmath57 , and so it seems natural to generalize contraction in a way that allows initial transients in time and/or amplitude .
the next simple example demonstrates a system that does not satisfy  , but is a gcs .
[ exa : scalarsys ] consider the _
scalar _ time - varying system [ eq : scals ] x(t)=-(t)x(t ) , with the state @xmath0 evolving on  @xmath58 $ ] , and  @xmath59 is a class  k function ( i.e.  @xmath40 is continuous and strictly increasing , with  @xmath60 ) .
it is straightforward to show that this system does not satisfy   w.r.t . _
any _ norm ( note that the jacobian  @xmath61 satisfies  @xmath62 ) , yet it is  st , with @xmath63 , for any given @xmath43 .
the next section presents our main results .
the proofs are placed in section  [ sec : proofs ] .
the next three subsections study the three forms of  gcss defined above .      just like contraction
,  sost implies entrainment to a periodic excitation . to show this ,
assume that the vector field  @xmath37 in   is  @xmath34 periodic .
pick  @xmath64 .
define  @xmath65 by  @xmath66 .
in other words ,  @xmath67 maps  @xmath68 to the solution of   at time  @xmath69 for the initial condition  @xmath70 . then  @xmath67 is continuous and maps the convex and compact set  @xmath10 to itself , so by the brouwer fixed point theorem ( see , e.g.  ( * ? ? ?
* ch .  6 ) ) there exists  @xmath71 such that  @xmath72 , i.e.  @xmath73 .
this implies that   admits a periodic solution  @xmath74 with period  @xmath34 . assuming that the system is also  sost , pick  @xmath75 .
then there exists  @xmath44 such that @xmath76 for all  @xmath77 and all  @xmath78 .
taking  @xmath79 implies that every solution converges to  @xmath80 .
in particular , there can not be two distinct periodic solutions .
thus , we proved the following .
[ prop : st_entrain ] suppose that the time - varying system  , with state  @xmath0 evolving on a compact and convex state - space  @xmath81 , is sost , and that the vector field  @xmath37 is  @xmath34-periodic .
then for any  @xmath82 it admits a unique periodic solution  @xmath83 with period  @xmath34 , and  @xmath84 converges to  @xmath80 for any  @xmath85 .
since both  so and  st imply  sost , proposition  [ prop : st_entrain ] holds for all three forms of  gcss .
[ def : nc ] system   is said to be _ nested contractive _  ( nc ) on  @xmath10 with respect to a norm  @xmath19 if there exist convex sets  @xmath86 , and norms  @xmath87 , where  @xmath88 $ ] , such that the following conditions hold .
a.   @xmath89 } \omega_\zeta=\omega$ ] , and [ eq : setsinc ] _ _ 1 _ _ 2 , _ 1 _ 2 .
b.   for every  @xmath90 there exists  @xmath91 $ ] , with  @xmath92 as  @xmath93 , such that for every  @xmath77 and every  @xmath94 [ eq : enter ] x(t , t_1,a ) _ , tt_1 + , and   is contractive on  @xmath95 with respect to  @xmath96 . c.   [ item : cc ] the norms @xmath96 converge to  @xmath97 as  @xmath98 , i.e. , for every  @xmath99 there exists  @xmath100 , with  @xmath101 as  @xmath102 , such that @xmath103    eq .
means that after an arbitrarily short time every trajectory enters and remains in a subset  @xmath95 of the state space on which we have contraction with respect to  @xmath96 .
[ thm : qcon ] if the system   is nc w.r.t .
the norm  @xmath19 then it is sost w.r.t .
the norm  @xmath19 .
the next example demonstrates theorem  [ thm : qcon ] .
it also shows that as we change the parameters in a contractive system , it may become a  gcs when it hits the `` verge '' of contraction .
this is reminiscent of an asymptotically stable system that becomes marginally stable as it looses stability .
[ exa : bio_smith ] consider the system @xmath104 where  @xmath105 , and @xmath106 as explained in  ( * ? ? ?
4 ) this may model a simple biochemical feedback control circuit for protein synthesis in the cell .
the  @xmath107s represent concentrations of various macro - molecules in the cell and therefore must be non - negative .
it is straightforward to verify that  @xmath108 implies that  @xmath109 for all  @xmath36 .
let  @xmath110 , and for  @xmath111 let @xmath112 we show in section  [ sec : proofs ] that if [ eq : km1 ] k-1 < then   is contractive on  @xmath113 w.r.t .
the scaled norm  @xmath114 for all  @xmath111 sufficiently small .
if @xmath115 then   does not satisfy  , w.r.t . any norm , on  @xmath113 , yet it is  sost on  @xmath113 w.r.t . the norm  @xmath116 .
note that for all  @xmath117 , [ eq : gderi ] g(x_n)= = g(0 ) .
thus   implies that the system satisfies   if and only if the `` total dissipation '' @xmath40 is strictly larger than  @xmath118 .
using the fact that  @xmath119 for all  @xmath120 it is straightforward to show that the set @xmath121\times[0 , ( \alpha_1 \alpha_2)^{-1 } ] \times\dots\times [ 0 , \alpha^{-1 } ]   ) \ ] ] is an invariant set of the dynamics for all  @xmath122 .
thus ,  , with  @xmath123 , admits a unique equilibrium point  @xmath124 and @xmath125 this property also follows from a more general result  ( * ? ? ?
4.2.1 ) that is proved using the theory of irreducible cooperative dynamical systems .
yet the contraction approach leads to new insights .
for example , it implies that the distance between trajectories can only decrease , and can also be used to prove entrainment to suitable generalizations of   that include periodically - varying inputs .
[ exa : phos : relay ] consider the system @xmath126 where  @xmath127 , and  @xmath128 . in the context of phosphorelay
@xcite ,  @xmath129 is the strength at time  @xmath130 of the stimulus activating the hk , @xmath131 is the concentration of the phosphorylated form of the protein at the @xmath132th layer at time  @xmath130 , and @xmath133 denotes the total protein concentration at that layer . note that @xmath134 is the flow of the phosphate group to an external receptor molecule .    in the particular case where  @xmath135 for all  @xmath132   becomes the _ ribosome flow model _  ( rfm ) @xcite .
this is the mean - field approximation of an important model from non - equilibrium statistical physics called the _ totally asymmetric simple exclusion process _  ( tasep ) @xcite . in the rfm ,
@xmath136 $ ] is the normalized occupancy at site  @xmath132 , where
@xmath137 [ @xmath138 means that site  @xmath132 is completely free [ full ] , and  @xmath139 is the capacity of the link that connects site  @xmath132 to site  @xmath140 .
this has been used to model mrna translation , where every site corresponds to a group of codons on the mrna strand ,  @xmath131 is the normalized occupancy of ribosomes at site  @xmath132 at time  @xmath130 , @xmath129 is the initiation rate at time  @xmath130 , and @xmath139 is the elongation rate from site  @xmath132 to site  @xmath140 .    our original motivation for generalizing   was to prove entrainment in the rfm  @xcite . for more results on the  rfm ,
see  @xcite .
assume that there exists  @xmath141 such that  @xmath142 for all  @xmath143 .
let  @xmath144\times\dots\times[0,p_n]$ ] denote the state - space of  . then
, as shown in section  [ sec : proofs ] ,   does not satisfy  , w.r.t .
any norm , on  @xmath10 , yet it is  sost on  @xmath10 w.r.t .
the  @xmath23 norm .    considering theorem  [ thm : qcon ] in the special case where all the sets  @xmath95 in definition  [ def : nc ] are equal to  @xmath10 yields the following result .
[ coro : new ] suppose that is contractive on  @xmath10 w.r.t . a set of norms  @xmath96 ,  @xmath88 $ ] , and that condition  ( [ item : cc ] ) in definition  [ def : nc ] holds . then   is sost on  @xmath10 w.r.t .
@xmath19 .
corollary  [ coro : new ] may be useful in cases where some matrix measure of the jacobian  @xmath20 of turns out to be non positive on  @xmath10 , but not strictly negative , suggesting that the system is `` on the verge '' of satisfying  .
the next result demonstrates this for the time - invariant system [ eq : time_in_var_sys ] = f(x ) , and the particular case of the matrix measure  @xmath145 induced by the  @xmath23 norm . recall that this is given by   with the  @xmath146s defined in  .
[ prop : new_meas_zero ] consider the jacobian  @xmath147 of the time - invariant system  .
suppose that  @xmath10 is compact and that the set  @xmath148 can be divided into two non - empty disjoint sets  @xmath149 and  @xmath150 such that the following properties hold for all  @xmath22 :    1 .   for any  @xmath151 , @xmath152 ; [ item : s0 ] 2 .   for any  @xmath153 , @xmath154 ; [ item : sminus ] 3
for any  @xmath155 there exists an index  @xmath156 such that  @xmath157 .
[ item : rec ]    then   is sost on  @xmath10 w.r.t .
the  @xmath23 norm .
the proof of proposition  [ prop : new_meas_zero ] is based on the following idea . by compactness of  @xmath10
, there exists  @xmath158 such that [ eq : comp_ass ] c_j(j(x))<- , j s_- x. the conditions stated in the proposition imply that there exists a diagonal matrix  @xmath159 such that  @xmath160 for all  @xmath161 .
furthermore , there exists such a  @xmath159 with diagonal entries _ arbitrarily close _ to  @xmath162 , so  @xmath163 for all  @xmath153 .
thus ,  @xmath164 . now
corollary  [ coro : new ] implies sost .
[ exa : new_from_russo ] consider the system : @xmath165 where  @xmath166 , and  @xmath167 $ ] .
this is a basic model for a transcriptional module that is ubiquitous in both biology and synthetic biology ( see , e.g. ,  @xcite ) . here
@xmath168 is the concentration at time  @xmath130 of a transcriptional factor  @xmath169 that regulates a downstream transcriptional module by binding to a promoter with concentration  @xmath170 yielding a protein - promoter complex  @xmath171 with concentration  @xmath172 .
the binding reaction is reversible with binding and dissociation rates  @xmath173 and @xmath174 , respectively .
the linear degradation rate of  @xmath169 is  @xmath175 , and as the promoter is not subject to decay , its total concentration , @xmath176 , is conserved , so  @xmath177 .
the jacobian of   is  @xmath178 , and all the properties in prop .
[ prop : new_meas_zero ] hold with  @xmath179 and  @xmath180 .
indeed ,  @xmath181 for all  @xmath182 .
thus , is sost on  @xmath10 w.r.t . the  @xmath23 norm .
[ exa : sec_new_from_russo ] a more general example studied in  @xcite is where the transcription factor regulates several independent downstream transcriptional modules .
this leads to the following model : @xmath183 where  @xmath184 is the number of regulated modules .
the state - space is  @xmath185\times\dots\times [ 0,e_{t , n}]$ ] .
the jacobian of   is @xmath186 and all the properties in prop .
[ prop : new_meas_zero ] hold with  @xmath179 and  @xmath187 .
thus , this system is sost on  @xmath10 w.r.t . the  @xmath23 norm .
arguing as in the proof of proposition  [ prop : new_meas_zero ] for the matrix measure  @xmath188 induced by the  @xmath26 norm ( see  ) yields the following result .
[ prop : new_meas_zero_infty ] consider the jacobian  @xmath147 of the time - invariant system  .
suppose that  @xmath10 is compact and that the set  @xmath148 can be divided into two non - empty disjoint sets  @xmath149 and  @xmath150 such that the following properties hold for all  @xmath189 :    1 .
@xmath190 for all  @xmath191 ; [ item : s0_inf ] 2 .
@xmath192 for all  @xmath193 ; [ item : sminus_inf ] 3 .   for any  @xmath191
there exists an index  @xmath194 such that @xmath195 .
[ item : rec_inf ]    then   is sost on  @xmath10 w.r.t .
the  @xmath26 norm .
a natural question is under what conditions  so and  sost are equivalent . to address this issue
, we introduce the following definition .
[ eq : defntr ] we say that   is _ weakly expansive _ ( we ) if for each  @xmath158 there exists  @xmath196 such that for all  @xmath14 and all  @xmath64 [ eq : sep ] | x(t , t_0,a)-x(t , t_0,b )     [ prop : sepimp ] suppose that is we .
then is sost if and only if it is so .
[ rem : global_lip ] suppose that  @xmath37 in   is lipschitz globally in  @xmath10 uniformly in  @xmath130 , i.e. there exists  @xmath197 such that @xmath198 then by gronwall s lemma ( see , e.g.  ( * ? ? ?
* appendix  c ) ) @xmath199 for all @xmath8 , and this implies that   holds for  @xmath200 .
in particular , if  @xmath10 is compact and  @xmath37 is periodic in  @xmath130 then  we holds under rather weak continuity arguments on  @xmath37 .      for _ time - invariant _ systems whose
state evolves on a convex and compact set it is possible to give a simple sufficient condition for  st .
let  @xmath201 [ @xmath202 denote the interior [ boundary ] of a set  @xmath203 .
we require the following definitions .
we say that   is _ non expansive _
( ne ) w.r.t . a norm  @xmath19 if for all  @xmath204 and all  @xmath205 [ eq : exp ] |x(s_2,s_1,a)-x(s_2,s_1,b)|    with  @xmath206 replaced by  @xmath207 .
[ def : ic ] the time - invariant system   with the state @xmath0 evolving on a compact and convex set  @xmath208 , is said to be _ interior contractive _  ( ic ) w.r.t .
a norm  @xmath209 if the following properties hold :    a.   for every  @xmath210 , [ eq : cond_a_enu ] x(t , x_0 ) , t>0 ; b.   for every  @xmath211 , [ eq : mucom_inv ] ( j(x))<0 , where  @xmath212 is the matrix measure induced by  @xmath19 .    in other words ,
the matrix measure is negative in the interior of  @xmath10 , and the boundary of  @xmath10 is `` repelling '' . note that these conditions do not necessarily imply that the system satisfies   on  @xmath10 , as it is possible that  @xmath213 for some  @xmath214 . yet ,   does imply that   is  ne on  @xmath10 .
[ thm : time_invar ] if the system   is ic w.r.t .
a norm  @xmath19 then it is  st w.r.t .
@xmath19 .
the proof of this result is based on showing that ic implies that for each @xmath215 there exists  @xmath216 such that @xmath217 and then using this to conclude that for any  @xmath218 all the trajectories of the system are contained in a convex and compact set  @xmath219 . in this set
the system is contractive with rate  @xmath220 .
the next example , that is a variation of a system studied in  @xcite , demonstrates this reasoning .
[ exa : new_from_russo_inf ] consider a transcriptional factor  @xmath169 that regulates a downstream transcriptional module by irreversibly binding , at a rate  @xmath221 , to a promoter  @xmath222 yielding a protein - promoter complex  @xmath171 .
the promoter is not subject to decay , so its total concentration , denoted by  @xmath223 , is conserved .
assume also that  @xmath169 is obtained from an inactive form  @xmath224 , for example through a phosphorylation reaction that is catalyzed by a kinase with abundance  @xmath225 satisfying  @xmath226 for all  @xmath227 .
the sum of the concentrations of @xmath224 , @xmath169 , and  @xmath171 is constant , denoted by  @xmath228 , with  @xmath229 . letting  @xmath230 denote the concentrations of  @xmath231 at time
@xmath130 yields the model @xmath232 with the state evolving on  @xmath233\times[0,e_t]$ ] .
let  @xmath234 , and consider the matrix measure  @xmath235 .
a calculation yields @xmath236 so  @xmath237 , and @xmath238 letting  @xmath239 $ ] , we conclude that  @xmath240 for all  @xmath241 . for any  @xmath242 , @xmath243 , and arguing as in the proof of theorem  [ thm : time_invar ] ( see section  [ sec : proofs ] ) , we conclude that for any  @xmath43 there exists  @xmath244 such that @xmath245 in other words , after time  @xmath51 all the trajectories are contained in the closed and convex set  @xmath246\times[0,e_t]$ ] .
letting  @xmath247 yields  @xmath248 and @xmath249 so   is  st w.r.t .
@xmath250 .
as noted above , the introduction of gcss is motivated by the idea that contraction is used to prove asymptotic results , so allowing initial transients should increase the class of systems that can be analyzed while still allowing to prove asymptotic results .
the next result demonstrates this .
[ coro : attract ] if   is  ic with respect to some norm then it admits a unique equilibrium point  @xmath251 , and  @xmath252 for all  @xmath77 .
the proof of corollary  [ coro : attract ] , given in the appendix , is based on theorem  [ thm : time_invar ] .
consider the _ variational system _ ( see , e.g. ,  @xcite ) associated with  : @xmath253 an alternative proof of corollary  [ coro : attract ] is possible , using the lyapunov - finsler function  @xmath254 , where  @xmath209 is the vector norm corresponding to the matrix measure  @xmath255 in  , and the lasalle invariance principle described in  @xcite .
since  ic implies  st and this implies  sost , it follows from proposition  [ prop : st_entrain ] that ic implies entrainment to  @xmath34-periodic vector fields .
the next example demonstrates this .
[ exa : st_ent_inf ] consider again the system in example  [ exa : new_from_russo_inf ] , and assume that the kinase abundance  @xmath225 is a strictly positive and periodic function of time with period  @xmath34 . since we already showed that this system is  st
, it admits a unique periodic solution  @xmath80 , of period  @xmath34 , and any trajectory of the system converges to  @xmath80 .
figure  [ fig : ent_inf ] depicts the solution of   for  @xmath256 , @xmath257 , @xmath258 , @xmath259 ,  @xmath260 , and initial condition  @xmath261 .
it may be seen that both state - variables converge to a periodic solution with period  @xmath262 .
( in particular ,  @xmath263 converges to the constant function  @xmath264 that is of course periodic with period  @xmath34 . )
( solid line ) and  @xmath265 ( dashed line ) of the system in example  [ exa : st_ent_inf ] as a function of  @xmath130 . , height=264 ]    contraction can be characterized using a lyapunov - finsler function  @xcite .
the next result describes a similar characterization for  st . for simplicity
, we state this for the time - invariant system  .
[ prop : equ_lyap ] the following two conditions are equivalent .
the time - invariant system   is st w.r.t . a norm  @xmath19 .
[ imp_clf : a ] b.   [ imp_clf : b ] for any  @xmath43 there exists  @xmath49 such that for any  @xmath14 and any  @xmath266 on the line connecting  @xmath68 and  @xmath267 the solution of   with  @xmath268 and  @xmath269 satisfies [ eq : delta_cond ] |x(t+)|(-t ) |x(0)|,t0 .    note that   implies that the function  @xmath270 is a _
generalized _ lyapunov - finsler function in the following sense . for any  @xmath43
there exists  @xmath49 such that along solutions of the variational system : @xmath271    in the next section , we describe several more related notions and explore the relations between them .
it is straightforward to show that each of the three generalizations of contraction in definition  [ def : qcont ] implies that   is  ne .
one may perhaps expect that any of the three generalizations of contraction in definition  [ def : qcont ] also implies  wc . indeed
,  st does imply wc , because @xmath272 for all  @xmath273 if st holds ( simply apply the definition with  @xmath274 , @xmath275 , and  @xmath276 in  ) .
however , the next example shows that  so does not imply wc .
[ exa : eps ] consider the scalar system [ eq : shift ] =    -2x , & 0|x| < 1/2 , + - , & |x|1 ,    with  @xmath0 evolving on  @xmath58 $ ] .
clearly , this system is not wc .
however , it is not difficult to show that it satisfies the definition of  so with  @xmath277 .
the next result presents two conditions that are equivalent to  sost .
[ lem : eqdef ] the following conditions are equivalent .    1 .
[ item1 ] system   is sost on  @xmath10 w.r.t .
some vector norm @xmath278 .
[ item2 ] for each  @xmath279 there exists  @xmath280 such that @xmath281 for all  @xmath46 and all  @xmath14 .
[ item3 ] for each  @xmath42 and each  @xmath43 there exists  @xmath282 such that [ eq : qcontnew ] |x(t , t_1,a)-x(t , t_1,b)|_v ( 1 + ) ( - ( t -t_1 ) _ 1 ) |a - b|_v , for all  @xmath283 and all  @xmath14 .
[ fig : graphbn2 ] summarizes the relations between the various contraction notions .
_ proof of theorem  [ thm : qcon ] .
_ fix arbitrary @xmath111 and @xmath38 .
the function @xmath284 $ ] is as in the statement of the theorem . for each @xmath90 ,
let @xmath285 be a contraction constant on  @xmath95 , where we write @xmath286 here and in what follows .
pick  @xmath14 and  @xmath43 . by  ,
@xmath287 for all  @xmath288 , so @xmath289 in particular , @xmath290 from the convergence property of norms in the theorem statement , there exist  @xmath291 such that @xmath292 and @xmath293 , @xmath294 as @xmath295 . combining this with   yields @xmath296 note
that taking  @xmath297 yields @xmath298 now for  @xmath288 let  @xmath299 .
then @xmath300 where the last inequality follows from  .
since @xmath293 , @xmath294 as @xmath295 , @xmath301 for @xmath43 small enough .
summarizing , there exists  @xmath302 such that for all  @xmath303 $ ] @xmath304 for all  @xmath305 and all  @xmath306 .
now pick  @xmath307 . for any  @xmath143 ,
let  @xmath308 .
then @xmath309 and this completes the proof .
@xmath310        , so @xmath311 thus , @xmath312 suppose that  @xmath313 .
then for all  @xmath117 , @xmath314 combining this with   implies that there exists a sufficiently small  @xmath111 such that @xmath315 for all  @xmath316 , so the system is contractive on  @xmath317 w.r.t .
@xmath114 .
we now use theorem  [ thm : qcon ] to prove that   is sost .
since  @xmath323 and  @xmath324 , @xmath325 for  @xmath326 $ ] , let @xmath327 it is straightforward to verify that   satisfies condition  ( br ) in  ( * ? ? ?
* lemma  1 ) , and this implies that for every  @xmath43 there exists  @xmath328 such that @xmath329 for all  @xmath218 .
then @xmath330 we already showed that this implies that there exists a  @xmath99 and a norm  @xmath331 such that   is contractive on  @xmath332 w.r.t . this norm .
summarizing , all the conditions in theorem  [ thm : qcon ] hold , and we conclude that is  sost on  @xmath113 w.r.t .
@xmath116 .
@xmath310    _ analysis of the system in example  [ exa : phos : relay ] .
_ for @xmath85 , let  @xmath333 denote the solution of   at time  @xmath143 for the initial condition  @xmath334 .
pick  @xmath43 .
satisfies condition  ( br ) in ( * ? ? ?
* lemma  1 ) , and this implies that there exists  @xmath335 such that for all  @xmath85 , all  @xmath336 , and all  @xmath337 @xmath338 furthermore , if we define  @xmath339 , @xmath340 , then the  @xmath341 system also satisfies condition  ( br ) in ( * ? ? ?
* lemma  1 ) , and this implies that there exists  @xmath342 such that for all  @xmath85 , all  @xmath336 , and all  @xmath343 @xmath344 we conclude that after an arbitrarily short time  @xmath345 every state - variable  @xmath131 , @xmath346 , is separated from  @xmath347 and from  @xmath133 .
this means the following . for
@xmath348 $ ] , let @xmath349 note that  @xmath350 , and that  @xmath351 is a strict subcube of  @xmath10 for all  @xmath326 $ ] .
then for any  @xmath94 , and any  @xmath43 there exists  @xmath352 , with  @xmath353 as  @xmath354 , such that [ eq : zinomeps ] x(t , t_1,a ) _ , tt_1+a .
the jacobian of   satisfies  @xmath355 , where @xmath356 note that  @xmath357 is metzler , tridiagonal , and has zero sum columns for all  @xmath358 .
note also that for any  @xmath359 every entry  @xmath360 on the sub- and super - diagonal of  @xmath361 satisfies  @xmath362 , with  @xmath363 .
note also that there exist  @xmath364 such that  @xmath365 is singular ( e.g. , when  @xmath366 and  @xmath367 the second column of  @xmath20 is all zeros ) , and this implies that the system does not satisfy   on  @xmath10 w.r.t .
any norm .    by  ( * ? ? ? * theorem  4 ) , for any  @xmath368 $ ] there exists @xmath369 , and a diagonal matrix  @xmath370 , with  @xmath371 , such that   is contractive on  @xmath95 w.r.t .
the the scaled  @xmath23 norm defined by  @xmath372 .
furthermore , we can choose  @xmath54 such that  @xmath373 as  @xmath98 , and  @xmath374 as  @xmath375 . summarizing , all the conditions in definition  [ def : nc ] hold , so   is nc on  @xmath10 and applying theorem  [ thm : qcon ] concludes the analysis .
@xmath310    _ proof of proposition  [ prop : new_meas_zero ] . _ without loss of generality , assume that  @xmath376 , with  @xmath377 , so that  @xmath378 .
fix  @xmath379 .
let  @xmath380 with the  @xmath381s defined as follows .
for every  @xmath155 ,  @xmath382 and  @xmath383 .
all the other  @xmath381s are one .
let  @xmath384
. then  @xmath385 .
we now calculate  @xmath386 .
fix  @xmath387 .
then  @xmath388 , so @xmath389 where the inequality follows from the fact that  @xmath390 for all  @xmath391 , and for the specific value  @xmath392 we have  @xmath393 and  @xmath394 .
we conclude that for every  @xmath191 , @xmath395 .
it follows from property  [ item : sminus ] ) in the statement of proposition  [ prop : new_meas_zero ] and the compactness of  @xmath10 that there exists  @xmath158 such that  @xmath396 for all  @xmath153 and all  @xmath397 , so for  @xmath111 sufficiently small we have  @xmath398 for all  @xmath153 and all  @xmath397 .
we conclude that for all  @xmath111 sufficiently small ,  @xmath399 , i.e. the system is contractive w.r.t .
@xmath400 . clearly ,  @xmath401 as  @xmath402 , and applying corollary  [ coro : new ] completes the proof .
_ proof of proposition  [ prop : sepimp ] .
_ suppose that   is sost w.r.t .
some norm  @xmath403 .
pick  @xmath111 .
since the system is  we , there exists  @xmath404 such that @xmath405 for all @xmath406 $ ] .
letting  @xmath407 yields [ eq : nnsep ] | x(t , t_0,a)-x(t , t_0,b ) |_v ( 1 + ) ( - ( t -t_0 ) _ 2 )          _ proof of lemma  [ lem : strict ] . _
pick @xmath43 and  @xmath415 . since  @xmath10 is an invariant set ,
@xmath416 is also an invariant set ( see , e.g. ,  ( * ? ? ?
* lemma iii.6 ) ) , so   implies that @xmath417 for all  @xmath418 . since  @xmath419 is compact , @xmath420 .
thus , there exists a neighborhood  @xmath421 of @xmath422 , such that @xmath423 for all @xmath424 .
cover  @xmath10 by such  @xmath421 sets . by compactness of  @xmath10 , we can pick a finite subcover .
pick smallest  @xmath41 in this subcover , and denote this by  @xmath425 . then  @xmath426 and we have that @xmath427 for all @xmath9 .
now , pick  @xmath428 .
let @xmath429 .
then @xmath430 and this completes the proof of lemma  [ lem : strict ] .
@xmath310    we can now prove theorem  [ thm : time_invar ] .
we recall some definitions from the theory of convex sets .
let  @xmath431 denote the closed ball of radius  @xmath432 around  @xmath0 ( in the euclidean norm ) .
let  @xmath433 be a compact and convex set with  @xmath434 .
let  @xmath435 denote the _ inradius _ of  @xmath391 , i.e. the radius of the largest ball contained in  @xmath436 . for  @xmath437 $ ] the _ inner parallel set of  @xmath433 at distance
@xmath438 _ is @xmath439 note that  @xmath440 is a compact and convex set ; in fact ,  @xmath440 is the intersection of all the translated support hyperplanes of  @xmath433 , with each hyperplane translated `` inwards '' through a distance  @xmath438 ( see  ( * ? ? ?
* section  17 ) ) .
assume , without loss of generality , that  @xmath441 .
pick  @xmath43 .
let  @xmath442 .
by lemma  [ lem : strict ] , @xmath443 and  @xmath444 for all  @xmath445 .
let  @xmath446
. then  @xmath447 .
pick  @xmath448 .
we claim that  @xmath449 . to show this ,
assume that there exists  @xmath450 such that  @xmath451 .
then there is a point  @xmath452 on the line connecting  @xmath453 and  @xmath454 such that  @xmath455 .
therefore , @xmath456 and this is a contradiction as  @xmath457 .
we conclude that  @xmath458 .
let  @xmath459 .
then   implies that  @xmath460 .
thus , the system is contractive on  @xmath461 , and for all  @xmath14 and all  @xmath227 @xmath462 where  @xmath97 is the vector norm corresponding to the matrix measure  @xmath463 .
this establishes  st , and thus completes the proof of theorem  [ thm : time_invar ] .
@xmath310    _ proof of corollary  [ coro : attract ] .
_ since  @xmath10 is convex , compact , and invariant , it includes an equilibrium point  @xmath41 of  .
clearly ,  @xmath464 . by theorem  [ thm : time_invar ] , the system is  st . pick  @xmath465 and  @xmath43 , and let  @xmath49 . applying   with  @xmath466 yields @xmath467 for all  @xmath468 . taking
@xmath469 completes the proof .
@xmath310    another possible proof of corollary  [ coro : attract ] is based on defining  @xmath470 by  @xmath471 .
then for any  @xmath85 ,  @xmath472 is nondecreasing , and the lasalle invariance principle tells us that  @xmath473 converges to an invariant subset of the set  @xmath474 , for some  @xmath475 . if  @xmath476 then we are done .
otherwise , pick  @xmath341 in the omega limit set of the trajectory .
then  @xmath477 , so   implies that  @xmath478 is strictly decreasing .
this contradiction completes the proof .
_ proof of proposition  [ prop : equ_lyap ] .
_ pick  @xmath204 .
let  @xmath479\to \omega$ ] be the line  @xmath480 .
note that since  @xmath10 is convex ,  @xmath481 for all  @xmath482 $ ] .
let @xmath483 this measures the sensitivity of the solution at time  @xmath130 to a change in the initial condition along the line  @xmath80 .
note that  @xmath484 , and @xmath485 comparing this to   implies that  @xmath486 is equal to the second component ,  @xmath487 , of the solution of the variational system   with initial condition @xmath488    suppose that the time - invariant system   is st .
pick  @xmath43 .
let  @xmath49 .
then for any  @xmath489 and any  @xmath490 $ ] , @xmath491 dividing both sides of this inequality by  @xmath54 and taking  @xmath492 implies that [ eq : condwexp ] |w(t + , r ) | ( - t ) |b - a| , so @xmath493 this proves the implication   @xmath494 . to prove the converse implication ,
assume that   holds
. then   holds and thus @xmath495 so the system is  st .     _
proof of lemma  [ lem : eqdef ] . _
if   is sost then   holds for the particular case  @xmath496 in definition  [ def : qcont ] . to prove the converse implication , assume that   holds .
pick  @xmath497 .
let [ eq : deftaumin ] : = \ { , } , and let  @xmath49 . pick  @xmath498 , and let  @xmath499 .
then @xmath500 where the last inequality follows from  .
thus , @xmath501 and recalling that  @xmath497 were arbitrary , we conclude that condition  2 ) in lemma  [ lem : eqdef ] implies sost .    to prove that condition  3 )
is equivalent to sost , suppose that   holds . then for any  @xmath502 , @xmath503 so we have sost .
conversely , suppose that   is sost .
pick any  @xmath75 .
then there exists  @xmath504 such that for any  @xmath505 @xmath506 thus , for any  @xmath507 @xmath508 taking  @xmath12 sufficiently small such that  @xmath509 implies that   holds for  @xmath510 .
this completes the proof that   is equivalent to sost .
@xmath310                z.  aminzare and e.  d. sontag , `` contraction methods for nonlinear systems : a brief introduction and some open problems , '' in _ proc .
53rd ieee conf . on decision and control _
, los angeles , ca , 2014 , pp . 38353847 .
s.  bonnabel , a.  astolfi , and r.  sepulchre , `` contraction and observer design on cones , '' in _ proc .
50th ieee conf . on decision and control and european control conference _ , orlando , florida , 2011 , pp . 71477151 .
h.  l. smith , _ monotone dynamical systems : an introduction to the theory of competitive and cooperative systems _ , ser .
mathematical surveys and monographs.1em plus 0.5em minus 0.4emprovidence , ri : amer .
soc . , 1995 , vol .
41 .
y.  zarai , m.  margaliot , and t.  tuller , `` explicit expression for the steady state translation rate in the infinite - dimensional homogeneous ribosome flow model , '' _ ieee / acm trans .
computational biology and bioinformatics _
10 , pp . 13221328 , 2013 .    g.  poker , y.  zarai , m.  margaliot , and t.  tuller , `` maximizing protein translation rate in the nonhomogeneous ribosome flow model : a convex optimization approach , '' _ j. royal society interface _ , vol .  11 , no . 100 , p. 20140713


according to greek mythology the goddess of wisdom , pallas athene , emerged clad in full armor after hephaestus split open zeus s head . in much the same way the local group sprang forth suddenly , and almost complete , in chapter vi of the realm of the nebulae ( hubble 1936 , pp .
124 - 151 ) .
hubble describes the local group as `` a typical small group of nebulae which is isolated in the general field '' .
he assigned ( in order of decreasing luminosity ) m31 , the galaxy , m33 , the large magellanic cloud , the small magellanic cloud , m32 , ngc205 , ngc 6822 , ngc 185 , ic 1613 and ngc 147 to the local group , and regarded ic 10 as a possible member . in the 2/3 century since hubble s work , the number of known local group members has increased from 12 to 36 ( see table 1 ) by the addition of almost two dozen low - luminosity galaxies .
recent detailed discussions of individual local group galaxies are given in mateo ( 1998 ) , grebel ( 2000 ) and van den bergh ( 2000a ) .
hubble ( 1936 , p. 128 ) pointed out that investigations of the local group were important for two reasons : [ 1 ] `` [ t]he members have been studied individually , as the nearest and most accessible examples of their particular types , in order to determine the[ir ] internal structures and stellar contents '' . in the second place [ 2 ] , `` the [ g]roup may be examined as a sample collection of nebulae , from which criteria can be derived for further exploration '' .
small galaxy groups , like the local group , are quite common . from inspection of the prints of the _ palomar sky survey _
van den bergh ( 2002a ) has estimated that 16% of nearby galaxies are located in such small groups .
hubble ( 1936 , p. 128 ) emphasized that `` the groups [ such as the local group ] are aggregations drawn from the general field , and are not additional colonies superposed on the field '' . from its observed radial velocity dispersion of 61 @xmath0 8 km @xmath1
the local group is found to have a virial mass of ( 2.3 @xmath0 0.6 ) x @xmath2 m@xmath3 ( courteau & van den bergh 1999 ) . the zero - velocity surface of the local group has a radius of 1.18 @xmath0 0.15 mpc , but @xmath480% of the local group members are actually situated within 0.4 mpc of the barycenter of the local group , which is located between m31 and the galaxy . from its virial mass , and the integrated luminosity of group members of 4.2 x 10 @xmath5 l@xmath3 , the dynamical mass - to - light ratio of the local group
is found to be @xmath6 = 44 @xmath0 12 in solar units .
such a high m / l value is an order of magnitude larger than the mass - to - light ratio in the solar neighborhood of the galaxy .
this supports the conclusion by kahn & woltjer ( 1959 ) that the mass of the local group is dominated by invisible matter .
the local group is situated in the outer reaches of the virgo supercluster .
the nearest neighbor of the local group is the small antlia group ( van den bergh 1999a ) .
this tiny cluster is located at a distance of only 1.7 mpc from the barycenter of the local group , i.e. well beyond the zero - velocity surface of the group .
the antlia group has a mean radial velocity of + 114 @xmath0 12 km s@xmath7 .
the number of galaxies brighter than m@xmath8 = -11.0 in the local group is 22 , compared to only four such objects in the antlia group .
however , because the antlia group contains no supergiant galaxies like m31 and the milky way system , its integrated luminosity is @xmath4150 times smaller than that of the local group .
since the local group is a relatively small cluster , this result suggests that clusters of galaxies have a range in luminosities ( and masses ? ) that extends over at least four orders of magnitude .
the centaurus group ( van den bergh 2000b ) , at a distance of @xmath43.9 mpc , which contains m 83 ( ngc 5236 ) and centaurus a ( ngc 5128 ) , is the nearest massive cluster . if one assumes ngc 5128 and ngc 5236 to be members of a single cluster one obtains a total virial mass of 1.4 x 10@xmath9 m@xmath3 and a zero - velocity radius of 2.3 mpc for the centaurus group . in other words
the zero - velocity surfaces of the centaurus group and the local group would almost touch each other . however , karachentsev et al .
( 2002 ) conclude that the galaxies surrounding ngc5128 and m83 , respectively , actually form dynamically distinct clusterings .
if that is indeed the case then the total mass of these clusters is reduced to only 3 x 10@xmath10 m@xmath3 , and the radius of the zero velocity surface around cen a is less than 1.3 mpc .
a second massive nearby cluster contains ic342 and the highly obscured elliptical maffei 1 .
mccall ( 1989 ) concluded that `` it is likely that ic342 and maffei 1 had a significant impact on the past dynamical evolution of the major members of the local group '' .
more recently fingerhut et al . (
2003ab ) have , however , found that the galactic absorption in front of maffei 1 is lower than was previously believed . as a result the distances to maffei 1 , and its companions , are too large for these objects to have had significant dynamical inetractions with individual local group members since the big bang .
if the peculiar velocities of galaxies are induced by gravitational interactions , then one might have expect massive field galaxies to have a lower velocity dispersion than dwarfs .
data on the nearby field ( whiting 2003 ) do not appear to support this expectation .
whiting finds the mean radial velocity dispersion among field galaxies within 10 mpc to be 113 km s@xmath7 .
a much lower dispersion of @xmath430 km s@xmath7 for the local hubble flow has , however , been found by karachentsev et al .
it is noted in passing that the former value appears to be significantly larger than the radial velocity dispersion of 61 @xmath0 8 km @xmath1 that courteau & van den bergh ( 1999 ) found within the local group itself .
to first approximation the local group is a binary system with massive clumps of galaxies centered on m31 and on the galaxy .
van den bergh ( 2000 , p. 290 ) estimates a mass m(a ) = ( 1.15 - 1.5 ) x 10@xmath10 m@xmath3 for the andromeda subgroup of the local group , compared to m(g ) = ( 0.46 - 1.25 ) x 10@xmath10 m@xmath3 for the galactic subgroup .
more recently sakamoto , chiba & beers ( 2003 ) have given somewhat higher mass estimates .
if leo i is included they find m(g ) = ( 1.5 - 3.0 ) x @xmath2 m@xmath3 , compared to m(g ) = ( 1.1 - 2.2 ) x 10@xmath10 m@xmath3 if it is assumed that leo i is not a member of the galactic subgroup .
recent proper motion observations by piatek et al .
( 2002 ) suggest that the fornax dwarf spheroidal galaxy may not , as had previously been thought , be a distant satellite of the galaxy .
instead the data appear to indicate that fornax is a free - floating member of the local group that is presently near perigalacticon . within each of the two main local group subclusters
there are additional subclumps such as the m31 + m32 + ngc205 triplet , the ngc147 + ngc185 binary , and the lmc + smc binary . for the entire local group
one finds , at the 99% confidence level , that low - luminosity early - type dsph galaxies are more concentrated in subclumps than are late - type dir galaxies . in other words
most dir galaxies appear to be free - floating members of the local group , whereas the majority ( but not all ) dsph galaxies seem to be directly associated with either m31 or the galaxy .
it is not yet clear if the mean ages of the stellar populations in dsph galaxies are themselves a function of location .
van den bergh ( 1995 ) has tentatively suggested that there is some evidence to suggest that star formation in dsph galaxies in the dense regions close to m31 and the galaxy might typically have started earlier than star formation in remote dsph galaxies .
the data that are shown in table 2 show a strong correlation between the morphological types of faint galaxies with @xmath11 ( i.e. objects fainter than m32 ) and their environment .
this dependence was first noticed by einasto et al .
almost all sph + dsph galaxies are seen to be associated with the two dense subclusters within the local group , whereas most ir galaxies appear to be more - or - less isolated group members .
it seems quite possible possible ( cf .
skillman et al .
2003 ) that the faint dir and dsph galaxies have similar progenitors , and that the differences between them that are observed now are due to environmental factors that favored gas loss from those dwarfs that occurred in dense environments , i.e. near giant galaxies .    within the local group the m31 , m32
, m33 subgroup has a total luminosity l@xmath8 = 3.0 x 10@xmath5 l@xmath3 , which is significantly greater than that of the subgroup centered on the milky way system which has l @xmath8 = 1.1 x @xmath12 l@xmath3 .
these luminosities of the m31 and galactic subgroups account for 71% and 24% of the total local group luminosity , respectively .
it should , however , be noted that some uncertainty in the luminosity ratio of m31 , to that of the galaxy , is introduced by the fact that both of these systems are viewed edge - on . as a result
the internal absorption corrections ( which may be quite large ) are uncertain . nevertheless , the notion that m31 is more massive than the galaxy receives some support from the observation that m31 appears to have 450 @xmath0 100 globular clusters , compared to only 180 @xmath0 20 such clusters associated with the galaxy ( barmby et al .
2000 ) . furthermore ( freeman 1999 ) the bulge mass of m31 is 3.6 x 10@xmath5 m@xmath3 , which is almost twice as large as the 2 x 10@xmath5 m@xmath3 mass of the galactic bulge . on the basis of these results
one might have expected the total mass of the m31 subgroup of the local group to be two or three times larger than that of the galaxy subgroup .
surprisingly this does not appear to be the case evans & wilkinson ( 2000 ) .
using radial velocity observations evans et al .
( 2000 ) conclude that `` there is no dynamical evidence for the widely held belief that m31 is more massive - it may even be less massive '' .
more recently gottesman et al .
( 2002 ) have also concluded from dynamical arguments that the mass of m31 `` is unlikely to be as great as that of our own milky way '' .
these authors even make the heretical suggestion that m31 might not have a massive halo at all !
either the well - known perfidity of small - number statistics has mislead us about the relative masses of m31 and the galaxy , or the mass - to - light ratio in the milky way system is much higher than that of the andromeda galaxy .
if the latter conclusion is correct then one would have to accept the existence of significant galaxy - to - galaxy variations in the ratio of visible to dark matter among giant spirals .
it has been known for many years ( van den bergh 1969 ) that the halo of m31 contains some globular clusters that are much more metal rich than those that occur in the galactic halo .
perhaps the best know example of such an object is the luminous globular mayall ii .
furthermore the color - magnitude diagrams for individual m31 halo stars ( mould & kristian 1986 , pritchet & van den bergh 1988 , durrell , harris & pritchet 1994 ) all show that ( 1 ) halo stars have a wide range in metallicity , and ( 2 ) the mean metallicity of stars in the halo of m31 is surprisingly high .
the mean values of [ fe / h ] for m31 halo stars that were obtained by mould & kristian , by pritchet & van den bergh and by durrell , harris & pritchet are @xmath13>$ ] = -0.6 , -1.0 and -0.6 , respectively .
[ it is noted in passing that the halo of m31 does contain a metal - poor component which includes clusters such as mayall iv , some rr lyrae variables , and non - variable horizontal - brach stars ( sarajedini & van duyne 2001 ) .
the observation that the stars in the halo of m31 appear , on average , to be so much more metal - rich than than those in the galactic halo ( durrell , harris & pritchet 1994 ) suggests that these two giant spiral galaxies had quite different evolutionary histories .
the higher metallicity of m31 halo stars indicates that the building blocks of the andromeda halo had much higher masses than those of the galactic halo .
simulations of murali et al .
( 2002 ) show that a significant fraction of the mass that was originally in such merging ancestral galaxies will end up as intergalactic debris , and presumably also in an extended halo of the final merged object .
an independent check on the metallicities of m31 giants is provided by the recent spectroscopic observations undertaken by reitzel & guhathakurta ( 2002 ) .
these authors find that their spectra of stars in the halo of m31 have a mean metallicity @xmath13>~$]= -1.3 .
this value is significantly lower that that derived from photometric observations of stars in the halo of the andromeda galaxy .
a possible explanation for this difference is that insufficient correction was made for the fact that old metal - rich red giants are fainter ( and therefore more difficult to observe ) than are old metal - poor giants .
non homogeneity of the andromeda halo might also have contributed to the observed difference in the mean metallicities derived from photometry and from spectroscopy .
for example , reitzel & guhathakurta find four stars of solar metallicity in the halo of m31 .
they suggest that these objects might represent metal - rich debris from an accretion event . other evidence for such accretion events is provided by ibata et al .
( 2001 ) and ferguson et al .
recently yanny et al .
( 2003 ) have also found possible evidence for such a tidal stream in the galaxy .
this stream is located beyond the plane of the milky way at a distance of @xmath420 kpc from the galactic center .
an interesting clue regarding the origin of the difference between the andromeda and the galactic halos is provided by the observations of pritchet & van den bergh ( 1994 ) ] [ see figure 1 ] which show that m31 has an r@xmath14 profile out to @xmath4 20 kpc from its nucleus .
such a structure is likely to have resulted from violent relaxation .
this suggests that the overall morphology of the andromeda galaxy was determined by violent relaxation resulting from the early merger of two ( or more ) massive metal - rich ancestral objects . on the other hand
the main body of the milky way system may have been assembled  la eggen , lynden - bell & sandage ( 1964 ) ] , with its metal - poor halo forming  la searle & zinn ( 1978 ) via late infall and capture of `` small bits and pieces '' .
however , it appears that these fragments differed in a rather fundamental way from those that produced dwarf spheroidal galaxies .
shetrone , ct , & sargent ( 2001 ) find that dsph galaxies are iron - rich and have 0.02 @xmath15 [ @xmath16 /fe ] @xmath15 0.13 , compared to typical galactic values of [ @xmath16/fe ] @xmath4 0.28 dex over the same range of metallicities .
this shows that the bulk of the galactic halo stars can not have formed in dwarf spheroidal galaxies that subsequently disintegrated .
in particular fulbright ( 2002 ) finds that less than 10% of local metal - poor stars with @xmath17 < -1.2 $ ] have alpha - to - iron abundance ratios similar to those found in dsph galaxies . more generally tolstoy et al .
( 2003 ) conclude that the observed element abundance patterns make it difficult to form a significant proportion of the stars observed in our galaxy in small galaxies that subsequently merged to form the disk , bulge , and inner halo of the milky way .
bekki , harris & harris ( 2003 ) have studied the distribution of stars of various metallicities after the merger of two spirals .
however , their model is not likely to be applicable to the early merger of the ancestral objects of m31 .
the reason for this is that extended disks would be destroyed by frequent tidal interactions at large look - back times .
_ hubble space telescope _
images show that galaxies with obvious disks mostly have @xmath18 , whereas the vast majority of the objects seen at @xmath19 appear to have either compact or chaotic morphologies ( van den bergh 2002b ) .
so the bekki et al .
model , in which extended disks merge , is probably inappropriate for galaxies at @xmath19 , i.e. for mergers that took place more than 9 gyr ago . in their pioneering study of the metallicities of individual stars in galactic halos
mould & kristian ( 1986 ) also observed stars in the halo of the late - type spiral m33 .
their color - magnitude diagram suggested that the halo of m33 was very metal - poor and had @xmath13>~= -2.2 $ ] . this value is more than an order of magnitude lower than that for stars in the halo of the andromeda galaxy .    the discussion given above
may be summarised by saying that m31 may have formed from the early merger of the two or three most massive galaxies in the core of the andromeda subgroup of the local group .
the less - massive andromeda companions , such as m32 and ngc205 , may represent objects in the core of the andromeda subgroup which had such low masses that they managed to survive individually .
due to differences in evolutionary history the halo of m31 contains relatively metal - rich globular clusters , whereas the galactic halo does not . another difference between the
m31 and galactic globular cluster systems has been noted by rich et al .
( 2002 ) who find that m31 does not appear to contain globular clusters with extremely blue horizontal branches such as m92 in the galactic halo .
an additional example of major differences between globular cluster systems is provided by m33 and the lmc .
these two late - type disk galaxies have comparable luminosities ( m33  m@xmath8 = -18.9 , lmc   m@xmath8 = -18.5 ) but radically different globular cluster systems .
surprisingly the lmc globulars , which are both very old and quite metal - poor , appear to have disk kinematics ( schommer et al .
1993 ) . on the other hand the metal - poor globular clusters associated with m33 ( schommer et al .
1991 ) seem to have halo - like kinematics .
sarajedini et al .
( 1998 found that eight out of 10 globular clusters in their m33 halo sample had significantly redder horizontal branches than galactic globulars of similar metallicity .
this difference might be interpreted as a second parameter effect . alternatively , and perhaps
more plausibly , the m33 globulars may be a few gyr younger than their galactic counterparts .
on the latter interpretation the m33 halo globular clusters exhibit an unexpectedly large age dispersion of @xmath43 - 5 gyr .
it is presently a mystery why the m33 halo globular clusters would have formed a few gyr later than typical galactic halo and lmc disk globulars .
some light might eventually be shed on these questions by observations of the radial velocities of rr lyrae stars in the lmc , and perhaps also in the near future , of rr lyrae stars in m33 .
the main differences between the globular cluster systems of m33 ( sc ) and the galaxy ( sbc ) could perhaps be understood ( van den bergh 2002b ) by assuming that late - type galaxies take significantly longer to arrive at their final morphology than do spirals of earlier morphological types .
however , the great age of the lmc cluster system appears to conflict with this simple explanation .
the observations of mould & kristian ( 1986 ) ] appear to show that the field stars in the halo of m33 are extremely metal - poor and have @xmath13>~= -2.2 $ ] .
it would be important to confirm this result by new photometry and to compare this value with the mean metallicity of stars in the outer halo of the large magelanic cloud .
not unexpectedly the majority of local group dwarf galaxies are surrounded by small families of metal - poor globular clusters .
however , it is puzzling that the sagittarius dwarf has one globular cluster companion ( terzan 7 , [ fe / h ] = -0.36 ) that is quite metal - rich .
how could such a relatively high metallicity have been built up within a dwarf galaxy ?
the only other galactic halo ( r@xmath20 10 kpc ) globular clusters that are known to have metallicities higher than [ fe / h ] = -1.0 are pal 1 and pal 12 .
the latter object is itself suspected of also being associated with the disintegrating sagittarius dwarf ( irwin 1999 ) .
this speculation is supported by the observations of martnez - delgado et al .
( 2002 ) who have found that pal 12 is possibly imbedded in tidal debris of the sagittarius dwarf .
rosenberg et al . ( 1998ab ) have also shown that the cluster pal 1 is significantly younger than most other galactic globular clusters .
it is presently not clear which kind of evolutionary scenarios would allow halo clusters like pal 1 and pal 12 to attain such relatively high metallicities . from its present luminosity and morphological type
the dwarf elliptical galaxy m32 would have been expected to be embedded in a swarm of 10 - 20 globular clusters .
it is therefore puzzling that not a single globular cluster appears to be associated with this galaxy .
perhaps some of the innermost m32 clusters were dragged into its compact luminous nucleus by dynamical friction .
also loosely bound outer globulars , that were originally associated with m32 , might have been detached by tidal interactions with the main body of m31 .
such detached m32 clusters would remain in the halo of m31 and might be recognized by being unusually compact .
it would be very worthwhile to undertake a systematic search for such m32 clusters with small r@xmath21 values in the halo of m31 .    for the vast majority of galaxies the specific globular cluster frequency
s ( harris & van den bergh 1981 ) is less than 10 .
however , the fornax dwarf , which has 5 globulars associated with it , has @xmath22 .
recent work by kleyna et al .
( 2003 ) appears to indicate that the ursa minor dwarf may have an even higher s value .
if a dynamically cold clustering of stars that these authors find in umi is a globular cluster ( or a disintegrated cluster ) then s @xmath4400 for the umi system .
taken at face value this result suggests , that the fraction of the light of dwarf spheroidals , that is in the form of globular clusters , may be much higher in dwarf spheroidal galaxies than it is in more luminour ( massive ) systems .
if the milky way system had collapsed in the fashion advocated by eggen , lynden - bell & sandage ( 1962 ) , then one would have expected the stars and globular clusters in the galactic halo to exhibit a radial metallicity gradient . on the other hand the halo of the milky way system
would not be expected to have such a metallicity gradient if , as envisioned by searle & zinn ( 1978 ) , it had formed by the accretion of many `` bits and pieces '' . using the 1999 version of the globular cluster catalog of harris ( 1996 ) , van den bergh ( 2003 ) , found a possible hint for the existence of such a radial metallicity gradient among galactic halo ( r@xmath2310 kpc ) globular clusters . however , the reality of this a gradient is not supported by the more recent data contained in the 2003 version of harris s catalog . on the other hand the clusters in the main body of the galaxy ,
i.e. those with r@xmath2410 kpc , do appear to show a radial abundance gradient .
globulars at r@xmath25 kpc are , on average , found to be more metal - rich than those having 4.0 @xmath15 r@xmath26 kpc .
a kolmogorov - smirnov test shows that there is only a 4% probability that the metallicities of these inner and outer globular clusters samples were drawn from the same parent population . taken at face value the existence of a galactic metallicity gradient between 4 kpc and 10
kpc favors the suggestion that the els model provides an adequate description of the formation of the main body of the galactic halo , whereas the sz model predictions agree with the observed lack of a metallicity gradient in the region with r@xmath2310 kpc .
it is a curious ( and unexplained ) fact that the distribution of flattening values of globular clusters differs significantly from galaxy to galaxy .
both open and globular clusters in the lmc are , for example , on average more flattened than their galactic counterparts .
furthermore ( van den bergh 1983 ) the luminous clusters in the large cloud are typically more flattened than are the less luminous ones .
finally it is of interest to note that the most luminous globular in many local group galaxies also seem to be among the most flattened .
examples are the globular mayall ii ( @xmath27 = 0.22 ) in m31 , @xmath28 centauri ( @xmath27 = 0.19 ) in the galaxy , and ngc1835 ( @xmath27 = 0.21 ) in the lmc .
the presently known members of the local group exhibit a flat luminosity function with slope @xmath16 = -1.1 @xmath0 0.1 ( pritchet & van den bergh 1999 ) .
this value is significantly lower than the slope @xmath16 = -1.8 that is predicted by the cold cold dark matter theory ( klypin et al .
the low frequency of faint local group dwarfs has been confirmed in a recent hunt for new faint local group members by whiting , hau & irwin ( 2002 ) .
this observed deficiency of faint lg members suggests that many of the progenitors of low mass galaxies were destroyed before they had a chance to form significant numbers of stars .
alternatively it might be assumed that the missing faint galaxies can be identified with compact high - velocity clouds .
however , maloney & putman ( 2003 ) have recently shown that such objects , if they were located at distances of @xmath41 mpc , would be largely ionized .
these authors therefore conclude that the compact high - velocity clouds are not at cosmologically significant distances , but that they are instead associated with the galactic halo .
all attempts to search for evidence of star formation in compact high - velocity clouds ( e.g. simon & blitz 2002 ) have so far remained unsuccessful .
it is therefore the deficiency of faint local group members is real .
figure 2 appears to show ( van den bergh 2000a , p. 281 ) that the luminosity function of dsph galaxies in the local group is steeper than that for dir galaxies .
this conclusion should , however , be regarded as provisional because a kolmogorov - smirnov test shows that the difference between the luminosity distributions of local group dsph and dir galaxies is only significant at the 75% confidence level .
if the luminosity function of dsph galaxies is , indeed , steeper than that for dir galaxies , then future discoveries are most likely to turn up very faint dsph ( rather than dir ) members of the local group .
it would clearly be very important to undertake sky surveys in two ( or more ) colors to search for the signatures of the color - magnitude diagrams of extremely faint ( and so far undiscovered ) resolved dwarf members of the local group .
three distinct explanations might be invoked to account for the apparent excess of faint dsph galaxies among presently known local group members : ( 1 ) perhaps gas was more likely to escape from faint ( low - mass ) galaxies than from more massive objects . as a result
low - mass galaxies would most often end up as gas free dsph galaxies .
alternatively ( 2 ) the mass spectrum with which galaxies form might depend on environmental density in such a way that high density regions ( i.e. the neighborhood of m31 and the galaxy ) form a larger fraction of low mass objects ( the majority of which end up as dsph galaxies ) .
the latter assumption would be consistent with the work of trentham et al .
( 2001 ) and trentham & tully ( 2003 ) who found that the galaxian luminosity function of the dense virgo cluster is much steeper than that for less dense clusters such as the ursa major cluster and the local group .
on the other hand the view that dense regions produce galaxies with steep galaxian luminosity functions appears to conflict with the result of sabatini et al .
( 2003 ) , who find that the virgo cluster luminosity function seems to be steeper in the low density outer regions of the virgo cluster than it is in the high density core of this cluster .
this observation might , however , be accounted for by assuming that tides preferentially destroy fragile dwarfs in the cores of dense clusters .
finally , ( 3 ) and perhaps most plausibly , the gas in the progenitors of the missing dwarfs might have been photoevaporated during reionization .
it is difficult to tease out information on the morphological evolution of galaxies from the distribution of stars of various ages .
the central bulges of giant spirals , such as m31 and the galaxy , are dominated by old stars .
this supports the notion that these objects were built up inside out , with their bulges forming first and the disk possibly being accreted at a later time
. it would be very important to establish how old the first ( presumably quite metal - poor ) generation of galactic disk stars is .
this problem is made more intractable because such thin disk stars have to be disentangled from stars that are physically associated with , and embedded within , an older thick disk population .
in fact , tidal interactions might pump energy into ( and hence thicken ) an initially thin disk of very metal - poor stars .
_ hubble space telescope _ observations of galaxies at at large look - back times suggest ( van den bergh 2002b ) that most disk star formation occurs at @xmath291.5 , i.e. during the last 9 gyr .
one reason for the paucity of disks at larger redshifts is , presumably , that such extended structures would be destroyed by tidal forces during the frequent encounters between galaxies at high redshifts . from two slightly metal - poor stars with disk kinematics liu & chaboyer ( 2000 ) find a thin disk age of 9.7 @xmath0 0.6 gyr .
such an age is consistent with the ages of spiral disks that are inferred from the fact that the hst images of distant galaxies start to show obvious disks at z @xmath41.5 .
since bars are generally assumed to have formed from global instabilities in disks one would not expect to see barred galaxies at @xmath19 .
if bars can not form from from initially chaotic protodisks then bar formation might be delayed to even smaller look - back times .
this suspicion appears to be confirmed by observations ( van den bergh et al .
1996 , 2002 ) which seem to suggest that the frequency of barred galaxies declines precipitously beyond redhifts of z @xmath40.7 .
if this conclusion is correct then one would expect the bar of the lmc to be younger than 6 gyr .
this conclusion is consistent with ( but not proved by ) the observation ( smecker - hane et al .
2002 ) that bursts of star formation occurred in the bar of the large cloud 4 - 6 gyr and 1 - 2 gyr ago .
almost all local group galaxies are found to contain some very old stars like rr lyrae variables .
this shows that these galaxies started to form stars quite early in the history of the universe , i.e. more than @xmath410 gyr ago .
the best candidate for a `` young '' local group member is the dwarf irregular leo a. however , dolphin et al .
( 2002 ) have discovered a few rr lyrae variables , which have ages @xmath30 gyr , in this galaxy .
furthermore , recent observations by schulte - ladbeck et al . ( 2002 ) show that this object also contains some metal - poor red horizontal branch stars .
this clearly demonstrates that leo a is not a young galaxy .
m32 and many of the local group dsph galaxies have not experienced recent star formation . on the other hand the gas rich spiral and irregular group members are still forming stars at a significant rate .
qualitative data on the past rate of star formation in such galaxies can be obtained from both their integrated colors and from the intensity ratios of various spectral lines . however , the hope that the age distribution of star clusters in galaxies could provide more detailed information on the past rate of star formation has been shattered by the work of larsen & richtler ( 1999 , 2000 ) which appears to show that the rate of cluster formation varies as a rather high power of the rate of star formation .
in other words there is not a one - to - one correspondence between the rate of cluster formation and the general rate of star formation . in the local group
this phenomenon is beautifully illustrated by the difference between the quiescent dwarf irregular ic1613 , which contains few star clusters of any kind ( i.e. baade 1963 , p. 231 ; van den bergh 1979 ) , and the large magellanic cloud that is presently forming both stars and clusters quite actively .
it seems likely that the present specific frequency of globular clusters in galaxies was mainly determined by their peak rates of star formation , with elevated peak rates resulting in high present specific cluster frequencies .
only fragmentary information is so far available on the luminosity evolution of individual local group members .
few star clusters in the lmc have ages between the 3.2 gyr age of the oldest open clusters and the @xmath413 gyr ( rich , shara & zurek 2001 ) age of the lmc globular clusters .
this probably means that the large cloud experienced a quiescent period that extended for @xmath4 10 gyr . during this `` dark age ''
no violent bursts of star formation ( which could have triggered the formation of star clusters ) occurred .
however , it is quite likely that a trickle of star formation ( such as that which presently occurs in ic1613 ) continued during the dark ages between 3.2 gyr and 13 gyr ago .
this speculation is supported by the data of da costa ( 2002 ) which seem to show that the metallicity in the lmc increased between the beginning and the end of the dark age , i.e. between the termination of globular cluster formation @xmath413 gyr ago , and the beginning of open cluster formation 3.2 gyr ago . a possibly more complicated scenario is hinted at by the work of smecker - hane et al .
( 2002 ) who conclude that star formation in the bar of the lmc was episodic , while the rate of star formation remained more or less constant within the disk of the large cloud . however , an important caveat is that the rate of star formation in lmc the disk was so low that the data do not provide strong constraints on the lmc disk star formation history .
the carina dsph galaxy seems to have experienced a major burst of star formation 7 gyr ago ( hurley - keller , mateo & nemeic 1998 ) ] .
however , at maximum this object probably only reached m@xmath8 @xmath4 -16 making it too faint to have become what babul & ferguson ( 1996 ) have called a `` boojum '' .
there has been a long controversy among astronomers regarding the nature ( or even the existence of ) a fundamental difference between open clusters and globular clusters .
the present consensus is that all clusters populations initially formed with an power law mass spectrum , and that globular clusters are simply the oldest and most massive population component that was best able to withstand the erosion caused by the destruction of lower mass clusters via evaporation , encounters with massive interstellar clouds and disk / bulge shocks .
however , a different scenario has been proposed by van den bergh ( 2001 ) .
he suggested that there have , in fact , been two ( perhaps quite distinct ) epochs of cluster formation . during the first of these globular clusters might have formed as halo gas was being compressed by shocks that were driven inwards by ionization fronts generated during cosmic reionization at @xmath31 .
such effects would presumably be greatest in the halos of small protogalaxies that are relatively easy to ionize . a second generation of massive clusters might have formed by the compression ( and subsequent collapse ) of giant molecular clouds that was triggered by heating of the interstellar medium induced by collisions between gas rich protogalaxies .
a similar view has recently been expressed by schweizer ( 2003 ) who also argues that the first generation of globular clusters formed nearly simultaneously from pristine molecular clouds that were heated and shocked by the strong pressure increase that accompanied cosmological reionization .
schweizer argues that this hypothesis might also account for the similarity of metal - poor globular clusters in all types of galaxies and environments .
both van den bergh ( 2001 ) and schweizer ( 2003 ) argue that second generation globular clusters were mainly formed during subsequent collisions and mergers between galaxies .
from dynamical arguments kahn & woltjer ( 1959 ) ] first showed that the local group can only be stable if it contains a significant amount of invisible matter . using radial velocity observations of local group members courteau & van den bergh ( 1999 )
have estimated that the local group has a total mass of ( 2.3 @xmath0 0.6 ) x 10@xmath10 m@xmath3 , from which the mass - to - light ratio ( in solar units ) is found to be m / l @xmath8 = 44 @xmath0 14 .
this high value shows that the total mass of the local group exceeds that of the visible parts of its constituent galaxies by an order of magnitude . in their
1959 paper kahn & woltjer suggested that this `` missing mass '' in the local group might be in the form of hot ( 5 x 10@xmath32 degrees ) low density ( 1 x 10@xmath33 protons @xmath34 ) gas , which would be difficult to detect observationally .
hui & haiman ( 2003 ) have recently shown that the thermal history of such gas has probably been quite complex . in recent years
the notion that hot gas is responsible for the missing mass in the local group has been overshadowed by the idea that this missing mass is actually in the form of cold dark matter ( blumenthal et al .
however , recent _ far ultraviolet explorer _
satellite observations ( nicastro et al . 2003 ) of the absorption lines of o vi ( which have radial velocities of only a few hundred km s @xmath7 ) suggest that hot gas may provide a non - negligible contribution to the missing mass in the local group .
alternatively sternburg 2003 the hot gas clouds observed by nicastro et al .
might just be reprocessed metal - enriched gas , that was ejected from the neighborhood of the galactic plane in supernova driven fountains .
if the latter suggestion is correct , then one might expect that the hot clouds would be relatively metal -rich .
on the other hand they would be expected to be metal - poor if they are composed of hot primordial gas .
the final fate of the local group has been discussed by forbes et al .
( 2000 ) who conclude that dynamic friction will eventually result in the merger of m31 and the galaxy .
this merged object will resemble an elliptical galaxy with m@xmath8 @xmath4 -21 that contains @xmath4 700 globular clusters . in a universe that continues to expand for ever ( bennett et al .
2003 , spergel et al .
2003 ) this object will , in the distant future , be the only remaining visible object in the universe .
* both the high metallicity of the m31 halo , and the r@xmath14 luminosity profile of the andromeda galaxy , suggest that this object might have formed from the early merger , and subsequent violent relaxation , of two ( or more ) relatively massive metal - rich ancestral objects . *
the main body of the galaxy may have formed in the manner suggested by eggen , lynden - bell & sandage ( 1962 ) , whereas its halo is more likely to have been assembled by accretion of `` bits and pieces '' in the manner first suggested by searle & zinn ( 1978 ) . *
it is profoundly puzzling that the old metel - poor globular clusters in the lmc appear to have been formed in an early disk , whereas the globulars associated with m33 seem to have originated in a slightly younger halo . *
it is speculated that the oldest generation of globular clusters in the universe might have formed as halo gas was compressed and heated in shocks that were driven inwards by ionization fronts generated during cosmic reionization . on the other hand second generation globular clusters formed as a result of the heating of molecular clouds during collisions between gas - rich galaxies .
it is emphasized that the history of star formation in galaxies is often very different from the history of cluster formation . *
it is presently not understood how globular clusters like terzan 7 ( which is associated with the sagittarius dwarf ) were able to attain a relatively high metallicity . neither do we know why the globular clusters associated with some galaxies are so much more flattened than are those in others .
* it is suggested that the specific globular frequency of galaxies was mainly determined by the peak rate of star formation during evolution . *
all local group galaxies appear to contain a very old population component , i.e. all nearby galaxies started to form stars just after they were formed .
in other words there are no truly young galaxies in the local group . *
the local group has a mass of ( 2.3 @xmath0 0.6 ) x 10@xmath10 m@xmath3 , a luminosity l@xmath8 = 4.2 x 10@xmath5 l@xmath3 , and a zero - velocity radius of 1.18 @xmath0 0.16 mpc .
most of the mass and luminosity of the local group appears to be concentrated in two subgroups that are centered on m31 and the galaxy , respectively .
there is presently a lively controversy about which of these two subgroups is the most massive .
if the galactic subgroup turns out to be more massive than the m31 group , then the ratio of dark to visible matter must differ significantly from group to group . * it is not yet clear if hot low density gas provides a significant contribution to the total mass of local group galaxies .
rich , r. m. , corsi , c. e. , ballazzini , m. , frederici , l. , cacciari , c. & fusi pecci , f. 2002 in extragalactic star clusters (= iau symposium no.207 ) , eds .
d. geisler , e. k. grebel , and d. minniti ( san francisco : asp ) , p.140                                                                                m31 & ngc 224 & sb i - ii & 760 & -21.2 + milky way & galaxy & sbc i - ii : & 8 & -20.9 + m33 & ngc 598 & sc ii - iii & 795 & -18.9 + lmc & ... & ir iii - iv & 50 & -18.5 + smc & ... & ir iv / iv - v & 59 & -17.1 + m32 & ngc 221 & e2 & 760 & -16.5 + ngc 205 & ... & sph & 760 & -16.4 + ic 10 & ... & ir iv : & 660 & -16.3 + ngc 6822 & ... & ir iv - v & 500 & -16.0 + ngc 185 & ... & sph & 660 & -15.6 + ic 1613 & ... & ir v & 725 & -15.3 + ngc 147 & ... & sph & 660 & -15.1 + wlm & ddo 221 & ir iv - v & 925 & -14.4 + sagittarius & ... & dsph(t ) & 24 & -13.3 + fornax & ... & dsph & 138 & -13.1 + pegasus & ddo 216 & ir v & 760 & -12.3 + leo a & ddo69 & ir v & 800 & -12.2 + sagdig & ... & ir v & 1180 & -12.0 + leo i & regulus & dsph & 250 & -11.9 + and i & ... & dsph & 810 & -11.8 + and ii & ... & dsph & 700 & -11.8 + aquarius & ddo 210 & v & 1025 & -11.3 + pegasus ii & and vi & dsph & 815 & -10.5 + and v & ... & dsph & 810 & -10.2 + and iii & ... & dsph & 760 & -10.2 + cetus & ... & dsph & 775 & -10.2 + leo ii & ... & dsph & 210 & -10.1 + pisces & lgs 3 & dir / dsph & 620 &  -9.8 + phoenix & ... & dir / dsph & 395 &  -9.8 + sculptor & ... & dsph & 87 &  -9.8 + tucana & ... & dsph & 895 &  -9.6 + cassiopeia & and vii & dsph & 690 &  -9.5 + sextans & ... & dsph & 86 &  -9.5 + carina & ... & dsph & 100 &  -9.4 + draco & ... & dsph & 79 &  -8.6 + ursa minor & ... & dsph & 63 &  -8.5 +