document: it is thought that the vast majority of stars are formed in star clusters ( lada & lada  2003 ) . during the collapse and fragmentation of a giant molecular cloud into a star cluster 
   , only a modest percentage ( @xmath2  % ) of the gas is turned into stars ( e.g. lada & lada  2003 ) . 
   thus , during the initial phases of its lifetime , a star cluster will be made up of a combination of gas and stars . 
   however , at the onset of stellar winds and after the first supernovae explosions , enough energy is injected into the gas within the embedded cluster to remove the gas on timescales shorter than a crossing time ( e.g. hills  1980 ; lada et al . 
   1984 ; goodwin  1997a ) . 
   the resulting cluster , now devoid of gas , is far out of equilibrium , due to the rapid change in gravitational potential energy caused by the loss of a significant fraction of its mass . 
   while this process is fairly well understood theoretically ( e.g. hills  1980 ; mathieu  1983 ; goodwin  1997a , b ; boily & kroupa 2003a , b ) , its effects have received little consideration in observational studies of young massive star clusters . 
   in particular , many studies have recently attempted to constrain the initial stellar mass function ( imf ) in clusters by studying the internal dynamics of young clusters . by measuring the velocity dispersion and half - mass radius of a cluster , and assuming that the cluster is in virial equilibrium , an estimate of the dynamical mass can be made . by then comparing the ratio of dynamical mass to observed light of a cluster to simple stellar population models ( which require an input imf ) one can constrain the slope or lower / upper mass cuts of the imf required to reproduce the observations . studies 
   which have done such analyses have found discrepant results , with some reporting non - standard imfs ( e.g. smith & gallagher  2001 , mengel et al . 
   2002 ) and others reporting standard kroupa  ( 2002 ) or salpeter  ( 1955 ) type imfs ( e.g. maraston et al .  2004 ; 
   larsen & richtler  2004 ) .    however , bastian et al .  ( 2006 ) noted an age - dependence in how well clusters fit standard imfs , in the sense that all clusters @xmath1100  myr were well fit by kroupa or salpeter imfs , while the youngest clusters showed a significant scatter . 
   they suggest that this is due to the youngest ( tens of myr ) clusters being out of equilibrium , hence undercutting the underlying assumption of virial equilibrium needed for such studies .    in order to test this scenario , in the present work we shall look at the detailed luminosity profiles of three young massive clusters , namely m82-f , ngc  1569-a , & ngc  1705 - 1 , all of which reside in nearby starburst galaxies . m82-f and ngc  1705 - 1 
   have been reported to have non - standard stellar imfs ( smith & gallagher  2001 , mccrady et al . 
   2005 , sternberg  1998 ) . 
   here we provide evidence that they are likely not in dynamical equilibrium due to rapid gas loss , thus calling into question claims of a varying stellar imf . 
   ngc  1569-a appears to have a standard imf ( smith & gallagher  2001 ) based on dynamical measurements , however we show that this cluster is likely also out of equilibrium . throughout this work 
   we adopt ages of m82-f , ngc  1569-a , and ngc  1705 to be @xmath3  myr ( gallagher & smith  1999 ) , @xmath4  myr ( anders et al . 
   2004 ) and 1020  myr ( heckman & leitherer  1997 ) respectively . 
   studies of star clusters in the galaxy ( e.g.  lada & lada  2003 ) as well as extragalactic clusters ( bastian et al . 
   2005a , fall et al . 
   2005 ) have shown the existence of a large population of young ( @xmath5  10 - 20  myr ) short - lived clusters . 
   the relative numbers of young and old clusters can only be reconciled if many young clusters are destroyed in what has been dubbed `` infant - mortality '' . 
   it has been suggested that rapid gas expulsion from young cluster which leaves the cluster severely out of equilibrium would cause such an effect ( bastian et al . 
   we provide additional evidence for this hypothesis in the present work . 
   the paper is structured in the following way . in   [ data ] and 
   [ models ] we present the observations ( i.e. luminosity profiles ) and models of early cluster evolution , respectively . in   [ disc ] we compare the observed profiles with our @xmath0-body simulations and in   [ conclusions ] we discuss the implications with respect to the dynamical state and the longevity of young clusters . 
   for the present work , we concentrate on _ f555w _ ( v ) band observations of m82-f , ngc  1569-a , and ngc  1705 - 1 taken with the _ high - resolution channel _ ( hrc ) of the _ advanced camera for surveys _ ( acs ) on - board the _ hubble space telescope _ ( hst ) . 
   the acs - hrc has a plate scale of 0.027 arcseconds per pixel . 
   all observations were taken from the hst archive fully reduced by the standard automatic pipeline ( bias correction , flat - field , and dark subtracted ) and drizzled ( using the multidrizzle package - koekemoer et al . 
   2002 ) to correct for geometric distortions , remove cosmic rays , and mask bad pixels . 
   the observations of m82-f are presented in more detail in mccrady et al . 
   total exposures were 400s , 130s , and 140s for m82-f , ngc  1569-a , and ngc  1705 - 1 respectively .      due to the high signal - to - noise of the data , we were able to produce surface brightness profiles for each of the three clusters on a per - pixel basis . 
   the flux per pixel was background subtracted and transformed to surface brightness . 
   the inherent benefit of using this technique , rather than circular apertures , is that it does not assume that the cluster is circularly symmetric . 
   this is particularly important for m82-f , which is highly elliptical ( e.g. mccrady et al . 
   2005 ) . for m82-f we took a cut through the major axis of the cluster . 
   the results are shown in the top panel of fig . 
   [ fig : obs ] . 
   we note that a cut along the minor - axis of this cluster as well as using different filters ( u , b , and i - also from _ hst - acs / hrc _ imaging ) would not change the conclusions presented in   [ disc ]  &    [ conclusions ] .    for ngc  1569-a and ngc  1705 - 1 
   we were able to assume circular symmetry ( after checking the validity of this assumption ) and hence we binned the data as a function of radius from the centre . 
   the results for these clusters are shown in the centre and bottom panels of fig . 
   [ fig : obs ] , where the circular data points represent mean binning in flux and the triangles represent median binning . the standard deviation of the binned ( mean ) data points is shown . 
   we also note that our conclusions would remain unchanged (   [ disc ]  &    [ conclusions ] ) if we used the _ f814w _ 
   ( i ) _ hst - acs / hrc _ observations . 
   we did not correct the surface brightness profiles for the psf as the effects that we are interested in happen far from the centre of the clusters and therefore should not be influenced by the psf . in all panels of fig . 
   [ fig : obs ] we show the psf as a solid green line ( taken from an _ acs - hrc _ observation of a star in a non - crowded region ) . 
   the background of the area surrounding each cluster is shown by a horizontal dashed line .    in order to quantify our results , 
   we fit two analytical profiles to the observed lps . 
   the first is a king  ( 1962 ) function , which fits well the galactic globular clusters and is characterised by centrally concentrated profiles with distinct tidal cut - offs in their outer regions . 
   the second analytical profile used is an elson , fall , & freeman  ( eff - 1987 ) profile , which is also centrally concentrated with a non - truncated power - law envelope . 
   the eff profile has been shown to fit young clusters in the lmc better ( eff ) as well as young massive clusters in galaxies outside the local group ( e.g. larsen  2004 ; schweizer  2004 ) . 
   the best fitting king and eff profiles are shown as blue / dashed and red / solid lines respectively . 
   the fits were carried out on all points within 0.5  of the centre of the clusters , i.e. the point at which , from visual inspection , the profile deviates from a smoothly decreasing function .    as is evident in fig . 
   [ fig : obs ] all cluster profiles are well fit by both king and eff profiles in their inner regions . 
   _ however , 
   none of the clusters appear tidally truncated , in fact all three clusters display an excess of light at large radii with respect to the best fitting power - law profile_. the points of deviation from the best fitting eff profiles are marked with arrows . 
   this result will be further discussed in   [ disc ] . 
   due to the rather large distance of the galaxies as well as the non - uniform background around the clusters presented here , background subtraction is non - trivial . 
   however , we have checked the effect of selecting different regions surrounding the clusters and note that our conclusions remain unchanged . 
   we also note that in the lmc , where the background can be much more reliably determined , many clusters show excess light at large radii ( e.g.  eff ; elson  1991 ; & mackey & gilmore  2003 ) . 
   we model star clusters using @xmath0-body simulations . 
   star clusters are constructed as plummer ( 1911 ) spheres using the prescription of aarseth et al . 
   ( 1974 ) which require the plummer radius @xmath6 and total mass @xmath7 to be specified . 
   clusters initially contain 30000 equal - mass stars , or the gravitational softening . 
   this is because the dynamics we model are those of violent relaxation to a new potential and so 2-body encounters are unimportant . ] . 
   simulations were conducted on a grape-5a special purpose computer at the university of cardiff using a basic @xmath0-body integrator ( the speed of the grape hardware means that sophisticated codes are not required for a simple problem such as this ) 
   .    the expulsion of residual gas from star clusters has been modelled by several authors ( see in particular lada et al . 
   1984 ; goodwin 1997a , b ; geyer & burkert 2001 ; kroupa et al . 
   2001 ; boily & kroupa 2003a , b ) . 
   the typical method is to represent the gas as an external potential which is removed on a certain timescale . 
   gas removal is expected to be effectively instantaneous , i.e. to occur in less than a crossing time ( e.g. goodwin 1997a ; melioli & de gouveia dal pino 2006 ) . as such we require no gas potential , and can model the cluster as a system that is initially out of virial equilibrium ( equivalent to starting the simulations at the end of the gas expulsion ) . 
   the subsequent evolution is the violent relaxation ( lynden - bell 1967 ) of the cluster as it attempts to return to virial equilibrium . 
   we define an _ effective _ star formation efficiency @xmath8 which parameterises how far out of virial equilibrium the cluster is after gas expulsion . a cluster which initially contains @xmath9  % stars and @xmath9  % gas ( i.e. a @xmath9  % star formation efficiency ) which is initially in virial equilibrium will have a stellar velocity dispersion that is a factor of @xmath10 - more generally @xmath11 - too large to be virialised after the gas is ( instantaneously ) lost . 
   we define the efficiency as effective , as it assumes that the gas and stars are initially in virial equilibrium , which may not be true . 
   we choose as initial conditions , @xmath12  pc ( corresponding to a half mass radius of @xmath13  pc ) and @xmath14 or @xmath15 ( i.e. the total initial stellar plus gas mass was @xmath16 or @xmath15 ) as representative of young massive star clusters . in order to compare the simulations with our observations 
   we place the simulations at our assumed distance of m82 , namely 3.6  mpc ( assuming that it is at the same distance as m81 - freedman et al . 
   previous simulations have shown that for @xmath17 clusters are totally destroyed by gas expulsion , but for higher @xmath8 significant ( stellar ) mass loss occurs , but a bound core remains ( goodwin 1997a , b ; boily & kroupa 2003a , b ) . for @xmath18 and @xmath19 
   respectively , @xmath20 and @xmath21  % of the initial stellar mass is lost within @xmath22  myr . 
   we confirm those results . 
   the escaping stars are not lost instantaneously , however . 
   stars escape with a velocity of order of the initial velocity dispersion of the cluster , typically a few  km  s@xmath23 . 
   therefore , escaping stars will still be physically associated with the cluster for @xmath24 
   @xmath25  myr _ after _ gas expulsion . 
   these stars produce a ` tail ' in the surface brightness profile and produce the observed excess light at large radii . 
   we assume a constant mass - to - light ratio for the simulation and convert the projected mass density into a luminosity and hence surface brightness profile . 
   the normalisation of the surface brightness is arbitrary and scaled so that the central surface brightness is similar to that of the observed clusters . 
   two of the simulations are shown in fig . 
   [ fig : model ] . 
   the filled circles are the surface brightness of the simulated cluster , with the specific parameters ( total initial mass , @xmath8 , and time since gas expulsion ) of the simulations shown . 
   we follow the same fitting technique as with the observations , namely fitting king and eff profiles ( dashed blue and solid red lines respectively ) to the profile . as was seen in the observations , the simulations display excess light at large radii . 
   the detailed correspondence between the observations and simulations presented here lead us to conclude that m82-f , ngc  1569-a and ngc  1705 - 1 display the signature of rapid gas removal and hence are _ not in dynamical equilibrium_. in future works we will provide a large sample of luminosity profiles of young massive extragalactic star clusters , as well as a detailed set of models which can be used to constrain the star formation efficiency of the clusters . here 
   we simply note that models with a sfe between 40 - 50% best reproduce the observations . 
   similar surface brightness profiles with an excess of light at large radii are seen in young lmc clusters : see eff and elson  ( 1991 ) in which many clusters clearly show these unusual profiles , and also mackey & gilmore  ( 2003 ) - in particular for r136 . 
   these profiles are also well matched by our simulations . 
   mclaughlin & van der marel  ( 2005 ) have compiled a data base of structural parameters for young lmc / smc clusters and compare the m / l ratio from dynamical estimates to that predicted by simple stellar population models ( i.e. to check the dynamical state of the young clusters ) . however , the study was limited as the young clusters tend to be of relatively low - mass , making it difficult to measure accurate velocity dispersions . here 
   we simply note that the five clusters in their sample younger than 100  myr all show significant deviations in the m / l ratio , but also note that this may simply be due to stochastic measurement errors . 
   it appears likely that the excess light at large radii seen in many massive young star clusters is a signature of violent relaxation after gas expulsion . 
   this suggests that these clusters have effective star formation efficiencies of around @xmath25 
   @xmath9  % , such that they show a significant effect , but do not destroy themselves rapidly . 
   it should also be noted that the escaping stars are not just physically associated with a cluster in the surface brightness profiles . 
   measurements of the velocity dispersion of the cluster will also include the escaping stars . 
   this will result in an artificially high velocity dispersion that reflects the initial total stellar _ and _ gaseous mass . 
   thus , mass estimates based on the assumption of stellar virial equilibrium may be wrong by a factor of up to three for 1020  myr after gas expulsion as is shown in fig . 
   [ fig : virial ] for @xmath26 and  @xmath27  % clusters ( i.e.  at the ages of ngc  1569-a and ngc  1705 - 1 ) .    clusters with @xmath28 @xmath27  % rapidly readjust to their new potential and the virial mass estimates become fairly accurate @xmath24  @xmath21  myr after gas expulsion ( i.e. for a cluster age of @xmath21 
   @xmath29  myr ) . 
   however , for @xmath30  % , the virial mass is significantly greater than the actual mass for @xmath31  myr and clusters do not settle into virial equilibrium for @xmath32  myr indeed , between @xmath33 and @xmath25  myr after gas expulsion the virial mass estimate _ underestimates _ the total mass by up - to @xmath33  % as the cluster has over - expanded . 
   a few recent studies have reported non - kroupa  ( 2002 ) or non - salpeter  ( 1955 ) type initial stellar mass functions ( imf ) in young star clusters ( e.g. smith & gallagher  2001 ; mengel et al . 
   these results were based on comparing dynamical mass estimates ( found by measuring the velocity dispersion and half - mass radius of a cluster and assuming virial equilibrium ) and the light observed from the cluster with simple stellar population models ( which assume an input stellar imf ) . 
   other studies based on the same technique ( e.g.  larsen & ritchler  2004 ; maraston et al . 
   2004 ) have reported standard kroupa- or salpeter - type imfs . 
   recently , bastian et al .  ( 2006 ) noted a strong age dependence on how well young clusters fit ssp models with standard imfs , with all clusters older than @xmath34  myr being will fit by a kroupa imf . 
   based on this age dependence , they suggested that the youngest star clusters ( @xmath35  myr ) may not be in virial equilibrium . 
   the observations presented here strongly support this interpretation as m82-f and ngc  1705 - 1 both seem to have been strongly affected by rapid gas loss . while ngc  1569-a has been reported to have a salpeter - type imf ( smith & gallagher 2001 ) 
   , the excess light at large radii suggests that this cluster has also undergone a period of violent relaxation and stars lost during this are still associated with the cluster even though its velocity dispersion correctly measures its mass ..    it should be noted that the obvious signature of violent relaxation in the profile of m82-f suggests that it is at the lower end of its age estimate of @xmath36  myr ( gallagher & smith 1999 ) , as by @xmath25  @xmath9  myr the tail of stars becomes disassociated from the cluster . 
   another possibility is that m82-f has been tidally shocked and has had a significant amount of energy input into the cluster , thus mimicking the effects of gas expulsion . 
   whichever is the case , the tail of stars from m82-f - whatever its age - is a signature of violent relaxation and strongly suggests that it is out of virial equilibrium . 
   if a young star cluster has a low enough effective star formation efficiency ( @xmath37  % ) it can become completely unbound and dissolve over the course of a few tens of myr . 
   this mechanism has been invoked to explain the expanding ob associations in the galaxy ( hills  1980 ) . 
   recent studies of large extragalactic cluster populations in m  51 ( bastian et al . 
   2005a ) and ngc  4038/39 ( fall et al . 
   2005 ) have shown a large excess of young ( @xmath510  myr ) clusters relative to what would be expected for a continuous cluster formation history . 
   both of these studies suggest that the excess of extremely young clusters is due to a population of short - lived unbound clusters . 
   the rapid dissolution of these clusters has been dubbed `` infant mortality '' . 
   the observations and simulations presented here support such a scenario . if the star formation efficiency is less than 30% - no matter what the mass - the rapid removal of gas completely disrupts a cluster ( although see fellhauer & kroupa  2005 for a mechanism which can produce a bound cluster with @xmath38% ) . 
   even if @xmath8 is large enough to leave a bound cluster , the cluster may be out of equilibrium enough for external effects to completely dissolve it , such as the passage of giant molecular clouds ( gieles et al . 
   2006 ) or in the case of large cluster complexes , other young star clusters . 
   interestingly , gas expulsion often significantly lowers the _ stellar _ mass of the cluster even if a bound core remains ( see   [ models ] ) . 
   thus , relating the observed mass function of clusters to the birth mass function needs to account not only for infant mortality , but also for ` infant weight - loss ' in which a cluster could lose @xmath39  % of its initial _ stellar _ mass in @xmath40  myr . 
   the current simulations do not include either a stellar imf , nor the evolution of stars . 
   the inclusion of these effects do not significantly effect the results as the mass - loss due to stellar evolution is low compared to that due to gas expulsion ( see goodwin 1997a , b ) . 
   in particular , we do not expect the preferential loss of low - mass stars as these clusters are too young for equipartition to have occured , thus stars of all masses are expected to have similar velocities . 
   one caveat to this is the effect of primordial mass ( hence velocity ) segregation which may mean that the most massive stars are very unlikely to be lost as they have the lowest velocity dispersion . 
   we will consider such points in more detail in a future paper . 
   observations of the surface brightness profiles of the massive young clusters m82-f , ngc  1569-a , and ngc  1705 - 1 show a significant excess of light at large radii compared to king or eff profiles . 
   simulations of the effects of gas expulsion on massive young clusters produce exactly the same excess due to stars escaping during a period of violent relaxation . 
   gas expulsion can also cause virial mass estimates to be significantly wrong for several 10s  of  myr . 
   these signatures are also seen in many other young star clusters ( e.g. elson  1991 ; mackey & gilmore  2003 ) and suggest that gas expulsion is an important phase in the evolution of young clusters that can not be ignored . 
   in particular , this shows that claims of unusual imfs for young star clusters are probably in error as these clusters are _ not _ in virial equilibrium as is assumed . 
   in future work we will further explore the dynamical state of young clusters in order to constrain the star - formation efficiency within the clusters . 
   we would like to thank mark gieles and francois schweizer for interesting and useful discussions , as well as markus kissler - patig and linda smith for critical readings of earlier drafts of the manuscript . 
   the anonymous referee is thanked for useful suggestions and comments . 
   this paper is based on observations with the nasa / esa _ hubble space telescope _ which is operated by the association of universities for research in astronomy , inc . under nasa contract nas5 - 26555 . 
   spg is supported by a uk astrophysical fluids facility ( ukaff ) fellowship . 
   the grape-5a used for the simulations was purchased on pparc grant ppa / g / s/1998/00642 . 
   99 aarseth , s.j . , hnon , m. , & wielen , r.  1974 , a&a , 37 , 183 anders , p. , de grijs , r. , fritze - v . 
   alvensleben , u. , & bissantz , n.  2004 , mnras , 347 , 17 bastian , n. , gieles , m. , lamers , h.j.g.l.m . , 
   scheepmaker , r. a. , & de grijs , r.  2005a , a&a 431 , 905 bastian , n. , saglia , r.p . , goudfrooij , p. , kissler - patig , m. , maraston , c. , schweizer , f. , zoccali , m ..  2006 , a&a , 448 , 881 boily , c.m . & kroupa , p.  2003a , mnras , 338 , 665 boily , c.m . 
   & kroupa , p.  2003b , mnras , 338 , 673 elson , r.a.w . 
   1991 , apjs , 76 , 185 elson , r.a.w . , 
   fall , m.s . 
   , & freeman , k.c . 
   1987 , apj 323 , 54 ( eff ) fall , s.m . , 
   chandar , r. , whitmore , b.c . 
   2005 , apj , 631 , 133 fellhauer , m. & kroupa , p.  2005 , mnras , 359 , 223 freedman , w. , hughes , s.m . , 
   madore , b.f . 
   1994 , apj , 427 , 628 gallagher , j.s . , 
   iii & smith , l.j .  1999 , 
   mnras , 304 , 540 geyer , m.p . 
   & burkert , a.  2001 , mnras , 323 , 988 gieles , m. , portegies zwart , s. f. , sipior , m. , baumgardt , h. , lamers , h.j.g.l.m . , leenaarts , j. , 2006 , mnras in prep goodwin , s.p .  1997a , mnras , 284 , 785 goodwin , s.p . 
   1997b , mnras , 286 , 669 heckman , t.m . & 
   leitherer , c.  1997 , aj , 114 , 69 hills , j.g .  1980 , apj , 235 , 986 king , i.  1962 , aj 67 , 471 koekemoer , a. m. , fruchter , a. s. , hook , r. n. , & hack , w. 2002 , in the 2002 hst calibration workshop , ed . s. arribas , a. koekemoer , & b. whitmore ( baltimore : stsci ) , 339  2001 , mnras , 321 , 699 kroupa , p.  2002 , science , 295 , 82 lada , c.j . , 
   margulis , m. , & dearborn , d.  1984 , apj , 285 , 141 lada , c.j . & lada , e.a . 
   2003 , ara&a , 41 , 57 larsen , s.s . 
   2004 , a&a , 416 , 537 larsen , s.s . & richtler , t.  2004 , a&a , 427 , 495 lynden - bell , d.  1967 , mnras , 136 , 101 mackey , a.d . & gilmore , g.f . 
   2003 , mnras , 338 , 85 maraston , c. , bastian n. , saglia r. p. , kissler - patig m. , schweizer f. , & goudfrooij p.  2004 , a&a , 416 , 467 mathieu , r.d . 1983 , apj , 267 , l97 mccrady , n. , graham , j.r . , & vacca , w.d . 
   2005 , apj , 621 , 278 mclaughlin , d.e . & 
   van der marel , r.p . 
   2005 , apjs , 161 , 304 melioli , c. & de gouveia dal pino , e. m.  2006 , a&a , 445 , l23 mengel , s. , lehnert , m.d . , 
   thatte , n. , & genzel , r.  2002 , a&a , 383 , 137 meylan g.  1993 , in asp conf . 
   48 , the globular cluste - galaxy connection . 
   smith , j.p . 
   brodie , astron . 
   , san fransisco , p.  588 
   plummer , h.c . 
   1911 , mnras , 71 , 460 salpeter , e.e .  1955 , apj , 121 , 161 schweizer , f.  2004 in asp conf . 
   322 , `` the formation and evolution of massive star clusters '' , eds . 
   lamers , l.j . 
   smith , a. nota , p.  111 
   smith , l.j . , & gallagher , j.s . 
   2001 , mnras , 326 , 1027 sternberg , a.  1998 , apj , 506 , 721
summary: we present detailed luminosity profiles of the young massive clusters m82-f , ngc  1569-a , and ngc  1705 - 1 which show significant departures from equilibrium ( king and eff ) profiles . 
   we compare these profiles with those from @xmath0-body simulations of clusters which have undergone the rapid removal of a significant fraction of their mass due to gas expulsion . 
   we show that the observations and simulations agree very well with each other suggesting that these young clusters are undergoing violent relaxation and are also losing a significant fraction of their _ stellar _ mass . 
   that these clusters are not in equilibrium can explain the discrepant mass - to - light ratios observed in many young clusters with respect to simple stellar population models without resorting to non - standard initial stellar mass functions as claimed for m82-f and ngc  1705 - 1 . 
   we also discuss the effect of rapid gas removal on the complete disruption of a large fraction of young massive clusters ( `` infant mortality '' ) . finally we note that even bound clusters may lose @xmath1  50% of their initial _ stellar _ mass due to rapid gas loss ( `` infant weight - loss '' ) . 
   [ firstpage ]    galaxies : star clusters  stellar dynamics 
   methods : @xmath0-body simulations