star clusters are vital laboratories for testing theories of stellar dynamics and stellar evolution  ( e.g. * ? ? ?
the globular clusters orbiting the milky way constitute a sample of about @xmath0 evolved stellar systems with ages @xmath1 years .
these provide useful constraints on the long term evolution of clusters but information about the paths which led to these end states can only be inferred via comparisons with the end points of numerical cluster simulations . in contrast , the system of rich star clusters in the large magellanic cloud ( lmc ) includes clusters of comparable mass to galactic globular clusters but which span the entire age range from less than @xmath2 years ( e.g. r 136 ) to about @xmath3 years ( e.g. ngc 2019 ) .
this allows us to observe massive star clusters at various epochs in their evolution and makes it possible to identify evolutionary trends , under the assumption that the external processes affecting cluster evolution have not changed significantly over the lifetime of the lmc . with this aim in mind ,
hst project go7307  ( e.g. * ? ? ?
* ) obtained deep images of four pairs of lmc clusters . within each selected pair
the clusters are approximately coeval , have similar masses and metallicities and lie at comparable radii from the lmc centre but differ significantly in other properties such as core radius .
the pairs range in age from about @xmath4 years up to @xmath3 years .
the youngest clusters ( ngc 1805 and ngc 1818 ) are old enough to have expelled all their residual gas and so are pure systems .    fig .
[ fig : lmccores ] shows a plot of the observed core radii of star clusters in the lmc as a function of their ages , with the clusters in the go7307 hst sample identified by name .
there is a clear trend that the spread in core radius is an increasing function of cluster age .
this trend was first identified by  @xcite from ground - based images and was recently confirmed by  @xcite using hst archive data .
@xcite calculated core radii for @xmath5 lmc clusters in a uniform manner with the aim of obtaining an unbiased data set .
they discuss in some detail the uncertainties in the data and conclude that selection effects are not responsible for the observed trend .
they suggest , as  @xcite had earlier argued , that the increased spread in core radius is an evolutionary phenomenon : all clusters formed with relatively small core radii ( @xmath6 pc ) and subsequently some clusters experienced core expansion while others did not .
@xcite discuss the possibility that inter - cluster variations in the slope of the initial mass function ( imf ) could have produced the observed spread in core radius .
they point out that mass loss from the increased numbers of massive stars found in clusters with imfs flatter than the standard salpeter imf could lead to core expansion of the observed magnitude on a timescale of about @xmath7 years . as  @xcite notes , however , the imf slope required to explain the largest core radii would lead to disruption of the clusters on a timescale of a few @xmath8 years . more directly , based on a detailed observational analysis of the clusters in the go7307 hst sample ,  @xcite
have shown that the slopes of their luminosity functions are similar to each other within the observational uncertainties , strongly suggesting that their imfs were very similar down to stellar masses of about @xmath9m@xmath10 .
the absence of inter - cluster population variations implies that stellar evolutionary mass loss can not be invoked to explain the observations .
as a star cluster evolves , encounters lead to a redistribution of energy among the stars .
heavier stars and binaries acquire lower velocities , on average , and sink towards the cluster core .
this process of mass segregation could enhance the heating effect of a population of primordial binaries because the binaries would tend to fall to the centre where encounters occur more frequently
. observations of young clusters  ( e.g. * ? ? ?
* ) and numerical simulations of star formation  ( e.g. * ? ? ?
* ) show that so - called primordial mass segregation can also result from the cluster formation process .
@xcite have shown that there is strong evidence for mass segregation in all the young and intermediate - age clusters in the go7307 sample .
further , the amount of mass segregation appears to be similar in all clusters and can not therefore account for differences in the cluster core radii .    the aim of this paper is to investigate , by means of simulations , whether dynamical processes could give rise to the observed core radius  age relation . in this context , the pair of clusters ngc 1868 and ngc 1831 are of particular interest .
these clusters have very similar ages and metallicities and their masses differ by at most a factor of two  ( see table  1 of * ? ? ?
however , as fig .
[ fig : lmccores ] shows , the core radius of ngc 1831 is almost three times larger than that of ngc 1868 .
given that the cluster ages are only around @xmath11myr any process which has affected the core radii must be capable of producing significant evolution even on this relatively short time scale . in this paper , we consider the dynamical effects of ( 1 ) variations in the external tidal field experienced by the clusters and ( 2 ) differences in the initial numbers of dynamically active binary stars in the clusters .
given that the crossing times at the half - mass radii of these clusters are about @xmath12 yrs and the times to orbit the lmc are about the same order of magnitude , it is reasonable to suppose that either of these processes could have begun to produce observable effects on cluster structure on the relevant time scale .    during the past decade ,
numerical simulations of star clusters have become increasingly realistic owing to dramatic increases in computational speeds and significant improvements in programming algorithms  ( e.g. * ? ? ?
the code used in this paper incorporates the stellar evolution models of  @xcite .
this makes it possible to analyse the output from the simulations in an identical manner to that used for observational data , thus allowing meaningful confrontation of models with observations .
this is particularly important when discussing cluster properties such as core radius because the relation between the theoretical and observational definitions of this quantity depend on the conversion between the mass function and luminosity function of the cluster stars . for a real cluster
this conversion is non - trivial because it depends on , among other quantities , the cluster metallicity and is further complicated by the possible presence of mass segregation  ( see e.g. * ? ? ?
* ) . in this paper ,
all comparisons between models and observations are made in terms of directly observable properties .    in the interests of clarity , we include here a note on the definition of the core radius .
as  @xcite discuss , theoreticians , numerical simulators and observers assign different meanings to the term `` core radius '' , reflecting the quantity of greatest relevance to their purpose . in theoretical work
, the core radius @xmath13 refers to the natural length scale of the system under study , for example the scale radius of the widely used  @xcite profile .
the operational definition of the core radius commonly used in discussions of @xmath14 simulations is that of the density radius @xmath15 .
this is the density - weighted average of the distance of each star from the density centre  @xcite , where the local density at each star is estimated from the mass within a sphere containing the star s six nearest neighbours and the density centre is defined as the density - weighted mean position of the stars .
observationally , since the available information is either in the form of cluster surface brightness profiles or luminosities of individual stars instead of mass densities and stellar masses , the core radius @xmath16 corresponds to the radius at which the projected surface brightness falls to half its central value . in the case of a cluster of equal - mass stars in equilibrium
, there is generally a well - defined relation between @xmath13 and @xmath16 .
for example , for a  @xcite model it is straightforward to show that @xmath17 .
it is important to note , however , that during the evolution of the cluster this relation is not static and that there is no general relation between @xmath15 and either @xmath13 or @xmath16 . in this paper
, we use the term `` core radius '' exclusively to refer to the observational quantity @xmath16 . in section
[ sec : tides_discussion ] we will compare the evolution of @xmath15 and @xmath16 , in order to emphasize the care which must be taken in comparing observations with simulations .     of lmc star clusters as a function of age .
the clusters in the go7307 hst sample are identified by name , and matched pairs are joined by dotted lines .
( following fig .  1 of @xcite
, data taken from @xcite).,width=317 ]    in section  [ sec : nbody ] we describe briefly the code used in the simulations , as well as the pipeline used to calculate observational properties of the simulated clusters .
section  [ sec : tides ] considers the effects of time variation of an external tidal field on the evolution of the cluster core radius while section  [ sec : binfrac ] deals with clusters moving in a steady tidal field but with different initial binary fractions . in section
[ sec : bintide ] we investigate the combined effect of both mechanisms by studying binary - rich clusters in a time - varying tidal field and the results of these simulations are discussed in the context of the lmc clusters in section  [ sec : bintidedisc ] . finally in section  [ sec : conclusions ] we summarise our conclusions .
.summary of runs performed .
@xmath19 , @xmath20 and @xmath21 are the initial number of stars , total stellar mass and primordial fraction of hard binaries , respectively .
the clusters were all simulated in the tidal field of a point mass representing the lmc , and each simulation ran until time @xmath22 . [ cols="<,^,^,^,^,^,^",options="header " , ]     simulations are a very valuable tool in the study of star clusters because they allow many different processes to be studied with a minimum of simplifying assumptions .
the simulations discussed below were performed using the nbody4 code  ( e.g. * ? ? ?
* ) running on the harp-3 and grape-6 special purpose computers  @xcite at the institute of astronomy , cambridge and on a grape-6 board at the american museum of natural history .
the nbody4 code uses the 4th order hermite integration scheme to follow the orbits of the cluster stars and utilises state - of - the - art regularisation schemes  @xcite for the accurate and efficient treatment of the internal motion of binary and higher - order subsystems .
it also incorporates a full set of stellar evolution models , in the form of parameterized functions , which are used to follow the evolution of single stars and binaries  @xcite .
of particular importance to the present study is the treatment of mass loss .
any gas which is lost from an evolving star and is not accreted on to a binary companion is assumed to leave the cluster instantaneously , thereby reducing the depth of the cluster potential well .
a realistic treatment of mass loss is therefore essential to allow the study of the evolution of such cluster properties as the core radius .    within the code
the equations of motion are integrated in scaled units such that @xmath23 and the total mass and virial radius of the cluster are also set to unity  @xcite . for a star cluster in virial equilibrium
the initial energy in these units is @xmath24 and the crossing time is @xmath25 .
given the total mass of the cluster in solar masses and an appropriately chosen length scale ( which determines the conversion from length units to physical units ) it is straightforward to obtain the conversion factors for time and velocity from units to myr and @xmath26 , respectively .    in order to facilitate direct comparison of the results of our simulations with the observed data , we determine the core radius for each model cluster using the same processing pipeline as that used by  @xcite . in summary ,
the stellar data from the model clusters ( in the form of position and luminosity information for each member star ) are binned radially to obtain a projected surface brightness profile .
the three - parameter model  @xcite @xmath27 is then fitted to the cluster profile . here
, @xmath28 is the central surface brightness , @xmath29 determines the fall - off of the surface brightness profile at large radii and @xmath16 is related to the length scale @xmath30 via @xmath31 .
full details of the fitting procedure used and discussion of the associated uncertainties are presented in  @xcite . in order to reduce the numerical noise due to the small number of stars in our simulated clusters
we observe each cluster along three perpendicular lines of sight and obtain the median core radius from these three estimates .
this helps to reduce the impact of individual bright stars which might otherwise skew the core radius estimation .
further , each run is performed four times using different random realisations of the initial cluster model and the results are averaged in order to ensure the robustness of any observed trends .
initially , all the model clusters are generated according to a  @xcite density law @xmath32 following the prescription of  @xcite .
here @xmath33 is the plummer length scale and @xmath34 is the central mass density .
observations of the young clusters in the lmc show that they do not appear to be tidally limited  @xcite . in this case , the plummer profile is a reasonable model to use for the initial density profile given that the fall - off in its outer regions lies between that of a tidally truncated king model and the observed profiles of the young lmc clusters .
the initial distribution of stellar masses is drawn from the imf of  @xcite which is a three - part power law fit to stellar data in the solar neighbourhood .
we take the lower and upper mass limits to be @xmath35m@xmath10 and @xmath36m@xmath10 , respectively .
we assume solar metallicity for the stars in our clusters .
this is an appropriate choice for the younger clusters  for example ,  @xcite have shown that ngc 1805 and ngc 1818 , whose ages are about @xmath37 and @xmath38myr respectively , are well matched by isochrones of approximately solar metallicity . since the pairs of clusters in the go7307 hst sample
were chosen to have similar metallicity , we know a priori that the observed range of cluster core radii can not be due to metallicity variations .
we have therefore chosen to perform all our simulations at solar metallicity , and postpone discussion of the effects of metallicity variations on star cluster evolution to future work .
table  [ tab : runs ] summarises the simulations which we have performed .
the clusters in these runs initially contain @xmath39 or @xmath40 stars and have a mean mass of about @xmath41 m@xmath10 .
clearly , the lmc clusters are more massive and contain greater numbers of stars than our model clusters and therefore direct comparison of the results of our simulations with the observed data must be approached with caution .
however , the effects of increasing @xmath18 in our simulations can be estimated by considering the various time scales in the problems we study .
detailed discussion of these scaling issues is presented in sections  [ sec : tides_discussion ] and  [ sec : bintidedisc ] where we argue that none of our conclusions is strongly dependent on the small size of our simulated clusters .
we first consider the effects of the external tidal field of the lmc on the evolution of the cluster core radii .
theoretical arguments and numerical simulations strongly suggest that tides have played a significant role in the evolution of star clusters and cluster systems in galaxies such as the milky way  ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
in addition , observations of tidal tails around at least twenty milky way globular clusters  ( e.g. * ? ? ? * ; * ? ? ?
* ) provide compelling evidence of the ongoing disruption of star clusters by the milky way .
the gravitational field of the lmc ( mass @xmath42m@xmath10 ) is much weaker than that of a galaxy like the milky way ( mass @xmath43m@xmath10 ) .
however , as  @xcite show ( their fig .  15 ) , there is a deficit of low - mass clusters older than @xmath3 years .
if we assume that the cluster population originally contained a spectrum of cluster masses similar to that of the present - day young clusters , their absence at late times strongly suggests that tides have been at work on the cluster system in exactly the manner expected from the models of  @xcite .
given that the go7307 hst sample of 8 clusters was selected on the basis that they lie at similar distances from the lmc centre but display a wide range of core radii , the steady tidal field experienced by a cluster on a roughly circular orbit can not be the origin of the observed trends .
further , given the intrinsic weakness of the lmc field at a radius of @xmath44kpc , where our clusters lie , one might not expect the tidal field to play a significant role in cluster evolution .
we therefore consider clusters on non - circular orbits .
such clusters experience tidal shocks during their orbits and we investigate whether these impart sufficient energy to the cluster to produce differences in the evolution of @xmath16 .
@xcite discusses the possible role of disk shocking in the internal evolution of the intermediate - age lmc clusters , but questions whether any of the clusters is on an appropriate disk - crossing orbit for this to be an important effect . in this section
we consider instead the tidal shocking produced by a point mass galaxy .    in order to investigate the role of the lmc tidal field in the evolution of its cluster system
, it is necessary to incorporate an appropriate external gravitational field into the simulations .
we model the lmc potential as that of a point mass of @xmath45m@xmath10 .
this is a rather drastic oversimplification of the actual situation .
we note , however , that this model exaggerates the effect of the tidal field , as any extended distribution of mass will have a more slowly varying gravitational potential .
further , in the radius range @xmath46  kpc , in which our model clusters orbit , the gradient of the potential is within a factor of two of that of the detailed lmc mass model of  @xcite .
thus our simulations constitute a plausible upper limit to the possible effects of tides .    in the simulations ,
the equations of motion are integrated in an accelerating but non - rotating reference frame centred on the centre of mass of the cluster and whose axes are aligned with fixed spatial directions .
this choice of frame simplifies the equations of motion for a cluster on a non - circular orbit despite having less symmetry than the more usual rotating frame used for steady tidal fields whose axes always point towards the attracting centre  @xcite .
implementation of the time varying tidal field requires some care because mass loss from evolving stars or the escape of stars from the cluster causes the centre of mass to acquire a non - zero velocity relative to the frame of integration .
it is therefore necessary to re - centre the frame of reference to coincide with the rest frame of the centre of mass following each such event .
all internal and external forces must then be reinitialized to ensure consistency with the re - centred frame during subsequent integrations .
we have performed simulations for clusters on elliptical orbits with perigalactic radii @xmath47kpc and apogalactic radii @xmath48kpc , probing the region in which most of our sample clusters lie but largely avoiding the complications of the inner regions of the lmc which may contain a stellar bar  ( e.g. * ? ? ?
we also modelled clusters on circular orbits at a radius of @xmath49kpc which corresponds to the time - averaged radius of the elliptical orbits .
this allows us to isolate the effect of the time variation of the tidal field because both sets of clusters experience the same mean tidal field .
we calculate the initial perigalactic tidal radius @xmath50 in the approximation that the cluster is on a circular orbit at a radius @xmath51 . in this case
, @xmath50 is determined from the cluster mass @xmath52 , the mass of the point mass galaxy @xmath53 and @xmath51 via  ( e.g. * ? ? ?
* ) @xmath54 the apocentric tidal radius @xmath55 is determined in a similar way from @xmath56 .
the unit of length , which determines the conversion between length units and physical units , is then chosen such that the limiting radius of all clusters is initially equal to @xmath50 . however , at the start of the simulations the clusters are placed at the apocentre of their orbits about the lmc which means that they do not fill their tidal radii . if , instead , the clusters initially extended to @xmath55 and were not truncated at @xmath50 then the removal of stars outside @xmath50 during the first few perigalactic passages could mask any heating effects at radii interior to @xmath50 due to the time variation of the tidal field .
in contrast to a king model , the plummer density law ( [ eqn : plummer ] ) formally extends to infinity and so the outer radius of a cluster is not well - defined .
the length unit we use yields a pericentric tidal radius close to the ( random ) position of the outermost star in our initial clusters .
the internal stellar velocities are scaled to ensure virial equilibrium at the start of the simulation .
the tidal radius of a cluster on an elliptical orbit does not have a unique value .
it varies continuously both as the cluster moves along its orbit and when the cluster loses mass through stellar evolution or the escape of member stars .
however , for clusters whose internal relaxation times are longer than their orbital periods , tidal stripping tends to truncate the star cluster at @xmath57  @xcite .
these two time scales are of similar magnitude for our simulated clusters so we expect two - body relaxation processes to lead to some expansion of the cluster beyond @xmath50 .
we take account of this by defining escaping stars to be those stars outside twice the apogalactic tidal radius @xmath55 .
stars whose orbits carry them outside this radius are generally lost to the cluster  ( e.g. * ? ? ?
stars are thus free to populate the region between @xmath50 and @xmath58 .
the choice of length unit determines the physical density of the cluster and hence the physical rate at which two - body relaxation processes affect the evolution of the cluster .
our aim is to investigate the influence of a time - varying external tidal field so it is important to ensure that all internal evolutionary processes proceed in an identical manner for each of our simulated clusters .
we therefore use the same length unit for the clusters on circular orbits as we chose for those on elliptical orbits .
this means that the clusters on circular orbits also do not initially fill their tidal radii and so necessarily allows expansion of these clusters .
however , since the circular orbits are placed at the time - averaged radius of the elliptical orbits this expansion should on average be the same for both sets of clusters .
any differences in the rates of expansion are due to the time variation in the tidal field .    the effect of steady tidal fields on star cluster evolution has been the subject of papers by a number of authors  ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
however , until recently the problem of clusters moving in time - varying tidal fields had received relatively little attention from simulations .
@xcite uses the results of some simulations of single - mass star clusters orbiting a point mass to determine the appropriate choice of escape radius for simulations of mass loss .
more recently ,  @xcite presented an extensive survey of the dynamical evolution of star clusters in external tidal fields .
our investigation is complementary to theirs because we are looking at the evolution of the observable global structure of the cluster while  @xcite focus their attention mostly on the evolution of the stellar mass function within the cluster .
the results of all our simulations are presented in terms of quantities which can be compared directly with those calculated from the observed data , thereby enabling us to discuss how different dynamical processes impact on observations .
we follow the evolution of the clusters for about @xmath59 ( where @xmath60 is the half - mass relaxation time ) this corresponds to @xmath61gyr in most cases . by this time , the clusters on elliptical orbits have lost more than @xmath36 per cent of their stars making estimates of properties such as the core radius increasingly noisy . at regular intervals during the evolution ,
the data from the clusters are recorded and all relevant properties of the clusters are calculated .
( upper curves ) and @xmath62myr ( lower curves ) .
there is clear evidence for mass segregation developing in the cluster.,width=317 ]    the solid curve in fig .
[ fig : circecc ] presents the evolution of the core radius for clusters on circular orbits about the lmc at a distance of about @xmath49kpc from the centre .
the broken curve shows the evolution of @xmath16 for clusters on elliptical orbits with the same time - averaged galactocentric radius .
it is apparent from this figure that the core radius evolution in both cases is indistinguishable , especially in view of the noise associated with the output of small-@xmath18 simulations . during the first @xmath63yrs
, the core radius is seen to decrease slightly from its initial value .
this is because the calculation of @xmath16 is dominated by the light of the most massive stars which rapidly become more centrally concentrated due to the effects of mass segregation . as fig .
[ fig : massseg ] shows , after @xmath64myr the half - mass radius of the more massive stars is already @xmath65 per cent smaller than that of the less massive stars .
after its initial contraction phase , the core radius begins to increase .
this occurs because the clusters do not fill the tidal radius corresponding to their galactocentric positions at the start of the simulations .
a combination of mass loss from the rapid evolution of massive stars and two - body relaxation can therefore drive the overall expansion of the cluster seen in fig .  [
fig : massseg ] .
the top panel of fig .
[ fig : velevol ] shows that the cluster s half - mass radius @xmath66 is increasing with time ( see fourth panel of fig .
[ fig : internal_props ] for long - term @xmath66 behaviour ) . during this period its @xmath15 (
not shown ) is roughly constant at about @xmath67pc .
( solid curve ) and core radius @xmath16 ( broken curve ) ; middle : plot of the virial ratio @xmath68 ( whose equilibrium value is @xmath69 ) ; bottom : velocity anisotropy parameter @xmath70 of stars within the inner @xmath37 per cent lagrangian radius ( broken curve ) and between the @xmath71 per cent and @xmath72 per cent radii ( solid curve ) .
see text for discussion.,width=317 ]    the expansion of @xmath66 is driven by a combination of two processes which are illustrated in the lower panels of fig .
[ fig : velevol ] .
the middle panel shows the evolution of the virial ratio @xmath68 .
initially the cluster is in equilibrium with @xmath73 .
after about @xmath74myr , the cluster is temporarily displaced from equilibrium ( @xmath75 ) by mass loss from the most massive stars .
this is accompanied by the increase in the half - mass radius of the cluster shown in the top panel of fig .
[ fig : velevol ] .
it is interesting to note that the @xmath16 of the cluster does not change as a result of this mass loss .
stellar evolution can not account for the sustained expansion of the cluster , however .
the cluster loses about @xmath76 per cent of its mass as a result of the above mass loss events . using the virial theorem it is straightforward to show that the gravitational radius @xmath77 of a stellar system with initial radius @xmath78 and mass @xmath52 which loses an amount of mass @xmath79 is given by  @xcite @xmath80 a @xmath76 per cent mass loss should lead to an increase of less than @xmath37 per cent in the size of the cluster which is insufficient to account for the cluster expansion seen in figs .
[ fig : circecc ] and  [ fig : massseg ] .
the bottom panel of fig .
[ fig : velevol ] shows the evolution of the velocity anisotropy parameter @xmath81  @xcite inside the @xmath37 per cent lagrangian radius of the cluster ( the radius containing @xmath37 per cent of the total mass of the cluster ) and in the outer parts ( between the @xmath71 per cent and @xmath72 per cent lagrangian radii ) .
initially , the velocity distribution throughout the cluster is isotropic ( @xmath82 ) .
while the central parts of the cluster retain this isotropy , in the outer parts the distribution becomes increasingly radially anisotropic ( @xmath83 ) .
this indicates that the continued growth of @xmath66 is being fueled by encounters in the inner parts of the cluster which produce a population of stars on extended radial orbits .
the velocity distribution in the outer parts of isolated clusters simulated by  @xcite exhibited similar behaviour during post core - collapse evolution ; those authors similarly concluded that the growth of the outer lagrangian radii at late times was due to two - body processes in the inner regions .
given the similarity between the core radius evolution of clusters on circular and elliptical orbits , we now proceed to investigate the extent to which the clusters on elliptical orbits have been affected by the variation in the tidal field . fig .
[ fig : internal_props ] presents the evolution of certain key internal properties of two typical clusters , one on a circular orbit and one on an elliptical orbit .
the top panel shows the galactocentric radius @xmath84 along the elliptical orbit to emphasize the clear correlation between perigalactic passages and changes in the internal properties of the cluster on the elliptical orbit .
the second panel shows that , as expected , during a perigalactic passage the rms velocity @xmath85 of the stars in the cluster increases sharply .
this velocity kick ejects the most energetic stars from the cluster . as a result of this mass loss
, @xmath85 settles to a value slightly below that of the cluster on the circular orbit once it returns to virial equilibrium .
the third and fourth panels of fig .
[ fig : internal_props ] show how @xmath60 and @xmath66 evolve during the cluster orbit .
the most noticeable effect is that the relaxation time for the cluster on the elliptical orbit is significantly reduced after the second perigalactic passage , so that the cluster evolves more rapidly following this encounter .
this is due both to the decrease in the total number of stars in the cluster and to the increase in the mean stellar mass @xmath86 ( shown in the bottom panel of fig .
[ fig : internal_props ] ) .
the increase in @xmath86 occurs because , on average , low - mass stars acquire larger velocities during pericentre passages and are hence more likely to escape from the cluster .
the value of @xmath66 is seen to increase during the perigalactic passage but rapidly settles back to a value slightly below the circular - orbit value .
thus while the tidal shocks experienced by the clusters on elliptical orbits do not affect the evolution of @xmath16 , they are nonetheless strong enough to modify other important properties of the cluster . as fig .
[ fig : internal_props ] shows , the impact of the shocks on the cluster properties increases with time  this is due to the decreasing number of bound stars in the cluster .     with time for a cluster on an elliptical orbit ( broken curve ) .
the solid line shows the position of the circular orbit of equal time - averaged radius .
( lower four panels ) evolution of rms velocity @xmath85 , half - mass relaxation time @xmath60 , half - mass radius @xmath66 and mean stellar mass @xmath86 .
curves for two cases are shown : ( 1 ) a cluster on a circular orbit ( solid curves ) and ( 2 ) a cluster on an elliptical orbit ( broken curves ) .
note that the half - mass radius is calculated using all stars within a radius of @xmath58 of the cluster centre . in the simulations all such stars are classified as non - escapers ,
although at any given time a fraction of them are formally unbound.,width=309 ]      before drawing our final conclusions based on these simulations , we now investigate whether the lack of impact on @xmath87 is due to the weakness of the shocks we have simulated or whether the value of @xmath16 is intrinsically robust to tidal shocking .
the strength of a tidal shock is generally presented in terms of the dimensionless parameter @xmath88 given by  ( e.g. * ? ? ?
* ) @xmath89 where @xmath90 is an estimate of the destruction time of the cluster owing to energy input from the gravitational shocks .
the time scale @xmath90 is given by @xmath91 where @xmath92 is the energy of the cluster evaluated at the half - mass radius and @xmath93 is the rate at which stars at the half - mass radius acquire energy due to the shock . in a simulation containing a finite number of stars , evaluation of this expression
requires some care . averaging the energy changes of stars near the half - mass radius
does not produce robust results because individual stars may gain and lose energy during the perigalactic passage through other means than interaction with the external tidal field .
instead we calculate the total energy binding the cluster as a function of time and estimate the energy input to the cluster during perigalactic passages .
[ fig : energy ] compares the time evolution of the binding energy @xmath94 of a cluster on a circular orbit with that of a cluster on an elliptical orbit . for a cluster in an external potential @xmath94
is defined as  ( e.g. * ? ? ?
* ) @xmath95 where @xmath96 and @xmath97 are the total kinetic and potential energy of the cluster , respectively , and @xmath98 is the energy in the tidal field
. the binding energy must be negative for the cluster to remain bound .
the effects of the external field are visible as dramatic changes in @xmath94 at about @xmath99 and @xmath100myr .
note that the sudden , temporary change in @xmath94 for the circular run at about @xmath101myr is due to a collision between two binaries that results in the the break - up of one of the binaries and the subsequent ejection of its components .
the ratio of the energy produced by this event to the binding energy of the cluster emphasises the relatively fragile nature of the clusters we are simulating .    .
curves for two cases are shown : ( 1 ) a cluster on a circular orbit ( solid curve ) and ( 2 ) a cluster on an elliptical orbit ( broken curve ) .
energies are given in standard units ( @xmath102).,width=317 ]     for clusters on circular ( solid curve ) and elliptical ( broken curve ) orbits in units of the initial mass @xmath20 .
enhanced mass loss during the period of revirialisation which follows each perigalactic passage is seen in the elliptical orbit case.,width=317 ]    during perigalactic passages , some stars acquire velocities above the escape velocity of the cluster and this leads to an enhanced mass - loss rate ( see fig .
[ fig : masstime ] ) .
however , the finite time required for such stars to escape means that @xmath94 can temporarily become positive because it is calculated by summation over all stars within twice the tidal radius of the cluster . after the cluster has passed pericentre , stars continue to leave across the tidal boundary and @xmath94 returns to a negative value indicating that the remaining cluster is still gravitationally bound . following each shock
, the cluster undergoes a revirialisation process during which the energy acquired from the tidal interaction is redistributed among the cluster stars .
this process results in further loss of mass from the cluster because , on average , the revirialisation process gives velocity kicks to stars which are three times larger than those arising from the initial shock  ( see e.g. * ? ? ? * ) .    , @xmath65 , @xmath36 , @xmath103 and @xmath72 per cent of the cluster mass .
the heating effect from the external tidal field is evident , even at the radius containing only @xmath65 per cent of the cluster mass.,width=317 ]    we take the energy input by the shock to be the difference between the energy of the cluster after this revirialisation process is complete and the energy of the cluster on the circular orbit at this time ( about @xmath104myr ) .
the time between the shocks is @xmath105myr , which leads to @xmath106myr . at this time @xmath107 for the cluster
is about @xmath108myr , so the shock strength is about @xmath109 . comparing this value with fig .
13 of  @xcite and noting that the initial concentrations @xmath110 of our clusters are about @xmath111 ( where we assume that @xmath112 is the initial outer limit of the clusters ) , we find that our clusters lie approximately on the borderline of the region where shocks are important . we conclude that the shocks our clusters are experiencing would significantly affect their evolution on long time scales . indeed , as fig .
[ fig : massradcomp ] illustrates , even the @xmath65 per cent lagrangian radius of the cluster is affected by the tidal shocks , indicating energy input into the cluster over a wide range of radii . nevertheless , on a timescale of up to @xmath113 yrs , the effects of tidal heating are insufficient to modify the observed @xmath16 of the clusters significantly .
( broken curve ) with the observational core radius @xmath16 ( solid curve ) for a single cluster .
top : un - smoothed results from the simulations .
middle : curves smoothed by averaging over five time bins .
bottom : ratio of @xmath16 to @xmath15 calculated using the smoothed curves from the middle panel .
see text for discussion.,width=317 ]    the top panel of fig .
[ fig : velevol ] illustrates an important point which is worth emphasizing in the context of comparisons between simulated data and observations , namely that @xmath16 and @xmath114 are not , in general , related in a simple way . in particular circumstances , for example at the start of our simulations where the model clusters are in strict dynamical equilibrium , it is easy to show that , for a plummer model composed of equal mass stars , the two quantities are related by @xmath115 .
however , fig .
[ fig : velevol ] shows that this relation quickly breaks down as the cluster evolves away from its initial conditions .
thus , in the absence of further information , observational estimation of @xmath87 must not be assumed to imply knowledge of the mass distribution in the cluster .
the same is true of the density radius @xmath15 defined in section  [ sec : introduction ] .
[ fig : corecomp ] shows the evolution of @xmath16 ( solid curve ) and @xmath15 ( broken curve ) for one of the clusters on an elliptical orbit . as the top panel of the figure shows , in a small-@xmath18 simulation , both @xmath16 and @xmath15
are intrinsically noisy quantities : this is emphasised by the middle panel where we have smoothed the curves from the top panel using a window function of width five time bins .
the bottom panel of fig .
[ fig : corecomp ] , which shows the ratio of @xmath16 to @xmath15 for the smoothed data , demonstrates that the relation between the two quantities is not a simple proportionality and that the evolutionary trends in the two quantities can be quite different .
the complicated relation between light and mass in star clusters is the result of a variety of processes such as stellar evolution and mass segregation which modify the mass and light distributions in different ways .
given the absence of simple relations between the different core radii , great care must be exercised when we interpret observational data through comparison with numerical simulations or theoretical models .
the simulations we have performed contain relatively low numbers of stars and are smaller and of lower total mass than the actual lmc clusters .
however , the mean densities of our clusters are very similar to those of the lmc clusters and it is this which determines the extent to which an external tidal field can influence the internal dynamics of a cluster . on this basis , we can be confident that conclusions based on our simulations are relevant to the clusters in the lmc .
our simulations have shown that , on a timescale of @xmath116gyr , which corresponds to several @xmath60 for our model clusters , tidal forces are not able to produce a spread in core radius similar to that observed among the lmc clusters .
there are two obvious avenues for building on the simulations we have performed .
first , simulations of clusters containing a larger numbers of stars @xmath18 could be performed .
this would allow us to achieve cluster masses comparable to those of the real lmc clusters .
however , increasing @xmath18 in a cluster while holding the physical density of the cluster constant leads to an increase in @xmath60 .
the expansion seen in @xmath16 for our simulated clusters is driven primarily by two - body processes ( which proceed on a time scale @xmath60 ) rather than by tidal heating and it will therefore be slowed by increasing @xmath18 .
further , the disruptive effect of pericentre passages depends on the ability of a cluster to re - fill the region outside its pericentric tidal radius @xmath50 during the time between close passages .
as @xmath60 increases this re - filling is reduced .
the cluster will therefore have fewer stars outside @xmath50 and hence will be less affected by shocks than small-@xmath18 clusters ( see fig .  2 of  @xcite ) . a cluster on an elliptical orbit and which
, as in our simulations , initially filled its pericentric tidal radius would remain truncated at this radius throughout its evolution because there would not be sufficient time between perigalactic passages for the region outside @xmath50 to be re - filled  @xcite .
thus shocks should have less impact on the @xmath16 of larger-@xmath18 clusters and the observed expansions of clusters on both circular and elliptical orbits should remain similar .    the second possible improvement would be to use a more realistic model of the lmc potential , for example one based upon the recent mass models of  @xcite . however , the actual lmc mass distribution is more extended than the keplerian potential we have been using in our simulations and therefore the gradients in its gravitational field are smaller .
consequently , the associated strength of its tidal effects on clusters moving on elliptical orbits is reduced .
thus , we conclude that neither of the above modifications to our simulations would tend to increase the effect of the tidal field .
the fact that the time - varying tidal field in our simulations did not result in changes in the evolution of @xmath16 in the extreme case of a low - mass cluster in orbit about a point - mass galaxy leads us to conclude that tides have not yet had time to play a significant role in the internal evolution of the intermediate - age lmc clusters .
in this section we present results for models moving in a steady tidal field but with different numbers of primordial binaries , in order to isolate the impact of large binary fractions on the evolution of the cluster properties .
the clusters discussed here move on the same orbits as the circular orbit runs of section  [ sec : tides ] .
based on an analysis of the observed binary sequence in the young cluster ngc 1818 ,  @xcite have shown that the binary fraction @xmath117 in the centre of this cluster is about @xmath118 . here @xmath119,where @xmath120 and @xmath121 are the numbers of single stars and binary systems , respectively . to date , the binary fractions of the older clusters in our sample have not been determined .
thus , at present the binary fractions of these clusters are essentially free parameters .
heggie s
@xcite law implies that only hard binaries are dynamically important in the evolution of a star cluster .
a binary is said to be hard when the magnitude of its binding energy @xmath122 exceeds the kinetic energy of a typical cluster star  @xcite . in the case of a binary with component masses @xmath123 and @xmath124 at a separation @xmath30 moving in a cluster with mean mass @xmath86 and velocity dispersion @xmath125 the hardness condition is @xmath126 ( we note that the nbody4 code uses a dispersion of @xmath127 in this equation to determine the hardness limit for binaries , in order to take account of the larger velocity dispersion in the core of a centrally concentrated cluster . ) in our simulations the boundary between hard and soft binaries occurs at a separation of about @xmath128au for binaries containing stars of the mean stellar mass in the cluster .
hard binaries act as heat sources because they become more tightly bound as a result of interactions with single stars and thereby impart kinetic energy to the stars they encounter .
@xcite found that the density radii of clusters containing primordial binaries experienced shallower contraction during core collapse and subsequently tended to be larger than those of clusters containing only single stars .
this result , based on simulations of single - mass clusters with binary fractions of about @xmath129 , motivates our investigation of clusters containing large numbers of primordial binaries .
hereinafter , we will use @xmath117 to refer to the fraction of _ hard _ binaries in the cluster .    as an extreme case
, we consider first clusters with initial @xmath130 and with the hard binaries randomly distributed throughout the cluster . in each binary
the mass of the primary star is drawn from the imf of  @xcite and the mass of its companion is obtained by randomly choosing a mass ratio for the binary from a uniform distribution constrained by the single star mass limits  ( for a more detailed discussion , see * ? ? ?
the binary orbital parameters are generated from the distribution of  @xcite .
the clusters are otherwise generated exactly as in the single star runs of section  [ sec : tides ] .
we note that the inclusion of a large number of hard primordial binaries dramatically increases the computation time required for each simulation .     per cent ,
@xmath131 per cent and @xmath36 per cent hard primordial binaries ( short - long dashed , long - dashed and solid curves , respectively ) .
all clusters are on circular orbits about the lmc .
see text for discussion.,width=317 ]     ( top panel ) and the energy binding the cluster @xmath94 ( middle panel ) .
two cases are shown ( 1 ) a cluster without primordial binaries ( broken curves ) and ( 2 ) a cluster containing @xmath36 per cent primordial binaries ( solid curves ) .
the total mass is normalised by its initial value @xmath20 and @xmath94 is given in standard units .
the bottom panel shows the evolution of the internal binding energy of the primordial binaries @xmath132 ( in standard units ) in the run containing primordial binaries .
the decrease in @xmath132 is the result of interaction - induced binary hardening .
note that @xmath133 includes the internal energy of escaped binaries .
see text for discussion.,width=317 ]     per cent primordial binaries ( solid curves ) and with single stars only ( broken curves ) . from bottom to top , the curves show the radii containing @xmath37 , @xmath65 , @xmath36 , @xmath103 and @xmath72 per cent of the cluster mass .
the heating effect of the primordial binaries is visible at all radii.,width=317 ]    fig .
[ fig : sinbin ] compares the evolution of @xmath16 for clusters containing @xmath134 single stars and no primordial binaries with that for clusters containing a significant initial binary population .
an increased binary fraction produces a noticeable effect on the evolution of @xmath16 .
after @xmath116gyr , the core radii of clusters containing @xmath36 per cent primordial binaries are about a factor of two larger than those of clusters with no binaries .
the results presented in fig .  [ fig : sinbincomp ] give some insight into the mechanism by which a large binary fraction can drive an expansion of @xmath16 .
the top panel shows that the rate at which mass is lost from a cluster with initial @xmath135 is significantly greater than in the case of a cluster with no initial binary population .
after @xmath116gyr this cluster has lost almost @xmath131 per cent more mass than the cluster containing only single stars .
the evolution of the energy binding the cluster @xmath94 is shown in the second panel of fig .
[ fig : sinbincomp ] . in a cluster with a large initial binary fraction
, @xmath94 rapidly attains a value significantly higher than that of a cluster of single stars .
the enhanced mass - loss rate due to the presence of the binaries decreases the depth of the cluster s potential - well ; the magnitude of the binding energy is further reduced by the heating of the remaining bound stars through weak encounters . as fig .  [
fig : sinbinpercent ] shows , the energy input due to binary heating is distributed throughout the cluster .
all the lagrangian radii of the cluster rapidly attain larger values when the cluster contains primordial binaries .
we note that , as was discussed in section  [ sec : tides_results ] , the growth of the clusters containing only single stars is due to the fact that the clusters do not fill their tidal radii at the start of the simulations .
we use identical initial conditions for the clusters containing binaries  the excess heating observed in those clusters is due to the presence of the hard binary population .
the third panel of fig .
[ fig : sinbincomp ] shows the internal binding energy of the primordial hard binaries @xmath132 as a function of time . in calculating @xmath132
we include the binding energy of those hard binaries which have escaped from the cluster because we are interested in the rate at which encounters lead to the hardening of binaries .
this plot demonstrates that , as expected , the hard binaries become harder with time and so impart energy to the cluster .
it is well known that energy generation by hard binaries halts core collapse  ( e.g. * ? ? ?
however , fig .
[ fig : sinbincomp ] demonstrates that even in the very early stages of a cluster s evolution , the heating effect of a large primordial hard binary population can have a significant impact on the energy budget of the cluster .
we note that the occasional large changes in @xmath132 are the result of stellar evolutionary processes and therefore do not contribute to the direct heating of the cluster .
for example , at @xmath136myr the components of one binary system merge to form a single star : prior to the merger the binary is hardened by tidal interactions between its components .
the massive star thus produced explodes as a supernova soon after and the velocity kick received by the remnant as a result of the explosion is reflected in the spike in @xmath137 .
the prominence of this event is due to the relatively small number of stars in the simulation .
[ fig : sinbin ] shows that the @xmath16 evolution of clusters with @xmath138 is very similar to that of clusters with @xmath135 initially .
this suggests that the heating effect of primordial binaries saturates above a certain value of the binary fraction @xmath139 . for larger values of @xmath117 ,
the @xmath16 evolution is unchanged , although the mass - loss rate is greater for increased @xmath117 . from fig .
[ fig : sinbin ] we see that @xmath139 for the clusters we have simulated appears to lie in the range @xmath140 because an initial binary fraction of @xmath37 per cent has a noticeably smaller impact on @xmath16 than one of @xmath131 per cent .
we are currently investigating this effect in more detail , including its dependence on the number of stars in the simulation : the results will be presented elsewhere . here
we note only that the saturation may arise from the fact that at early times only a relatively small number of binaries take part in the heating process . in our simulations only about three per cent of the binaries are dynamically hardened during the first @xmath141myr .
we find that this fraction appears to be independent of the value of @xmath117 . whether a particular binary experiences an encounter which leads to hardening depends on the masses and separation of the components as well as the local stellar density in the neighbourhood of the binary . for a fixed binary mass ,
however , the hardest binaries have the smallest separations and hence are least likely to interact with other stars .
this is reflected in our simulations : for a given binary mass , all the binaries which undergo early hardening are initially among the more weakly bound pairs .          in our final set of simulations we examine clusters containing large numbers of primordial binaries and moving on non - circular orbits about the lmc .
we have considered the extreme case of @xmath135 but , as was shown in the previous section , these results also apply to somewhat smaller values of @xmath117 . as fig .  [
fig : eccbincorecomp ] shows , the core radius of a cluster on an elliptical orbit and containing primordial binaries exhibits stronger growth than that of a cluster of single stars on a circular orbit : after about @xmath142myr the core radius of the former cluster is almost a factor of two larger .
the increase in core radius is solely due to the presence of binaries . as in the case of clusters containing only single stars ( section  [ sec : tides ] )
, a cluster with the same primordial binary fraction but on a circular orbit exhibits almost identical core radius evolution .
however , the increased mass - loss rate which results from binary heating leaves the cluster on the elliptical orbit more vulnerable to tidal disruption . after only @xmath142myr
, the cluster has lost almost @xmath143 per cent of its mass and the core radius determination is consequently very uncertain .
by @xmath116gyr it is virtually disrupted and retains only about @xmath37 per cent of its initial mass while the circular orbit cluster still contains about @xmath144 per cent of its initial mass .
our simulations demonstrate that the presence of a large population of hard binaries can produce significant expansion of the core radius of a cluster . in agreement with many previous studies  ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) we have thus shown that binaries constitute a dynamically significant population in star clusters .
in particular , we have found that primordial binaries can rapidly affect observable cluster properties such as @xmath16 . in the context of the lmc clusters , however , it is clear that the presence of binaries alone can not account for the the observed trends in @xmath16 and in particular for the factor of three difference in core radius between ngc 1831 and ngc 1868 . even in the extreme case in which one cluster has no primordial binaries at all
while the other cluster has a hard binary fraction of @xmath145 per cent , the difference in @xmath16 over the lifetime of the intermediate - age clusters is at most a factor of two .
the core radius of a cluster containing @xmath37 per cent binaries expands by only a factor of @xmath146 over the same time scale .
further , the similarity of the observed present - day luminosity functions of the clusters in the go7307 sample  @xcite argues against large inter - cluster variations in the imf and hence makes significant differences in binary fraction unlikely .
the initial radial distribution of hard binaries has an impact on their heating effect on the cluster stars  ( e.g. * ? ? ?
* ) and may also have a bearing on the saturation of the heating effect which we observed in section  [ sec : binfrac ] .
a binary population which was more centrally concentrated than the single stars as a result of either primordial or dynamical mass segregation would be expected to produce more heating because of the increased interaction rates in the inner parts of the cluster . in our future simulations
we will consider alternatives to the uniform distribution used in the present runs . however , based on our current simulations it is possible to estimate the significance of centrally concentrating the binary population .
the clusters with @xmath36 per cent binary fractions contain approximately the same number of binaries within @xmath147 as a cluster with a @xmath37 per cent overall binary fraction in which all the binaries are located within @xmath147 .
given that none of our simulations exhibited @xmath147 growth of a magnitude sufficient to match the observations , a core binary fraction in excess of @xmath36 per cent would be required to produce greater @xmath147 growth .
the observed similarity of the degree of mass segregation in the intermediate - age clusters  @xcite render it unlikely that the radial distribution of binaries could vary so greatly between clusters .    as in the case of the simulations in section  [ sec : tides ] , it is natural to ask whether our conclusions are also valid for clusters containing larger numbers of stars than those in our current simulations .
the simulation of large-@xmath18 clusters containing significant numbers of primordial binaries has only recently become feasible and we are currently investigating the impact of binaries on the evolution of such clusters ; the results of this study will be presented elsewhere .
the present simulations show that the disruption time for a small-@xmath18 cluster in a time - varying external tidal field is significantly reduced if it contains a large fraction of primordial binaries .
it is not trivial to extrapolate these results to larger-@xmath18 clusters .
however , binary heating of the cluster is essentially a two - body process and therefore proceeds on the relaxation time scale @xmath60 . given that a larger-@xmath18 cluster with the same mean physical density has a larger @xmath60 , the physical rate at which the cluster is heated by its binary population will be reduced .
thus , while larger-@xmath18 clusters would obviously survive for longer than those in our present simulations , their physical rate of @xmath16 expansion would be reduced .
large variations in primordial binary fraction are therefore unlikely to be the explanation of the lmc observations .
for the lmc cluster sample . clusters with @xmath16 larger ( smaller ) than @xmath148 pc are shown as open ( filled ) circles .
data are taken from @xcite.,width=317 ]    the core radii of star clusters in the lmc display a strong trend with age , namely that the older clusters exhibit a much greater spread in core radius than their younger counterparts .
in particular , the @xmath16 of the intermediate - age clusters ngc 1868 and ngc 1831 differ by about a factor of three . it has been suggested  @xcite that the trend represents a physical evolution of the clusters and that the clusters all formed with similar , small core radii which were modified by subsequent dynamical evolution .
the possibility that this might be caused by variations in the cluster imf has been ruled out based on hst photometry of a sample of eight clusters  @xcite . in this paper
we have used simulations of star clusters initially containing up to @xmath134 stars to investigate some possible alternative explanations .
we have considered two possible dynamical effects which could produce significant @xmath16 evolution : time - varying external tidal fields and a varying population of primordial binaries .    throughout this paper
we have analysed our simulations using the same procedures as were applied by  @xcite to the observed data .
we emphasised the importance of such an approach by comparing the evolution of the observational core radius to that of the half - mass radius and the density radius .
as has been pointed out by a number of authors  ( e.g. * ? ? ?
* ) , in general there are no simple relations between these quantities
. thus observations of simulated data sets are essential to ensure that comparisons between observed and simulated data are made in a consistent manner .
we considered clusters on circular and elliptical orbits of equal time - averaged radius about a point mass representing the lmc .
the clusters are scaled so that their limiting radii correspond to the pericentric tidal radius of the clusters on elliptical orbits .
however , at the start of each simulation the clusters lie well within the tidal radius which corresponds to their initial galactocentric position ( apocentre for the clusters on elliptical orbits and the circular radius for the clusters on circular orbits ) .
initially , the clusters expand due to stellar mass loss but the competing effect of mass segregation means that @xmath16 decreases slightly . at later times ,
two - body relaxation processes drive a further expansion of the clusters because the clusters are not tidally limited .
this expansion is mirrored by a growth in @xmath87 .
however , the rate of @xmath16 expansion is found to be identical irrespective of the orbit on which the cluster is moving . we conclude that differences in the cluster orbits can not be responsible for the observed differences in @xmath16 among the lmc clusters .    given that the binary fractions of all but the youngest clusters are currently unconstrained by observations , we compared the evolution of the @xmath16 of a cluster containing no primordial binaries with that of a cluster with a large initial hard binary fraction .
the presence of significant numbers of binaries leads to increased growth of @xmath16 compared with the evolution of clusters containing only single stars .
binary fractions of @xmath36 per cent and @xmath131 per cent produce similar degrees of @xmath16 expansion , while a @xmath37 per cent binary fraction results in noticeably less growth . even in the most extreme case of a cluster with a binary fraction of @xmath36 per cent on an elliptical orbit
, however , @xmath16 increases by at most a factor of two . not only is this smaller than the increase required to explain the lmc observations but the required spread in primordial binary fraction is already difficult to reconcile with the observed similarity of the cluster stellar luminosity functions  @xcite .
our simulations have demonstrated that neither large variations in primordial binary fraction nor differences in cluster orbits about the lmc can account for the @xmath16age relation seen in the lmc clusters .
the origin of the spread in @xmath16 among the intermediate - age clusters therefore remains an open issue .
we now return briefly to the observed data on the lmc cluster system as a whole to see whether there are any further clues to the origin of the trend .
[ fig : rhom ] presents evidence for a systematic difference in the central density of the large-@xmath16 clusters relative to the rest of the cluster population : these clusters have systematically lower central densities for a given cluster mass .
this could be indicative of variations in star formation efficiency ( sfe ) between the clusters : low sfe leads to rapid cluster expansion and consequently a drop in the central density  ( e.g. * ? ? ?
simulations of cluster expansion due to residual gas expulsion from clusters with low sfe  @xcite do not appear to show sustained core expansion of the magnitude seen in the intermediate - age lmc clusters .
however , the core radii quoted in  @xcite are those of the mass distribution in the simulated clusters which , as was discussed in section  [ sec : tides ] , may not accurately reflect the changes in the light distribution .
we are currently investigating the long - term evolution of clusters which experience different degrees of primordial gas loss owing to sfes in the range @xmath149 , required to explain the range of core radii seen in the youngest lmc clusters  @xcite . we are also simulating the effects of primordial binaries on larger-@xmath18 clusters .
of particular interest is the possibility that the heating effects of a large binary population could enable a cluster to retain a large initial @xmath16 produced by gas expulsion .
miw acknowledges financial support from pparc .
jrh acknowledges financial support from a kalbfleisch fellowship and thanks the institute of astronomy ( cambridge , uk ) for supporting a visit during this work .
adm would like to acknowledge the support of a trinity college ers grant and a british government ors award .
cat thanks churchill college for a fellowship .
we are very grateful to sverre aarseth for assistance with the nbody4 code and for numerous useful discussions .
we thank pete bunclark at the institute of astronomy , cambridge , for handling the installation and setup of the grape-6 hardware .
we also thank hongsheng zhao for helpful discussions and richard de grijs and douglas heggie for useful discussions and comments on an earlier version of this paper .