Evolution of multimass globular clusters in the Galactic tidal field with the effects of velocity anisotropy

We study the evolution of globular clusters with mass spectra under the influence of the steady Galactic tidal field, including the effects of velocity anisotropy. Similarly to single-mass models, velocity anisotropy develops as the cluster evolves, but the degree of anisotropy is much smaller than in isolated clusters. Except for very early epochs of the cluster evolution, the velocity distributions of nearly all mass components become tangentially anisotropic at the outer parts. We examine how the mass function (MF) changes in time. Specifically, we find that the power-law index of the MF decreases monotonically with the total mass of the cluster, in agreement with previous findings based on isotropic models or N -body studies. This is also consistent with the behaviour of the observed slopes of MFs for a limited number of clusters. We attempt to compare our results with multimass King models, although it is almost impossible to fit the entire density profiles for all mass components. When the MF is fixed, the central densities of individual components show significant differences between Fokker–Planck and King models. We obtain 'best-fitting' multimass King models, for which the central density of individual components as well as the total density distribution agrees with the Fokker–Planck models by adjusting the MF. The MFs obtained in this way closely resemble the MF within the half-mass radius of the Fokker–Planck result. Also, we find that the local MFs predicted by Fokker–Planck calculations vary more rapidly with radius than best-fitting multimass King models. The projected velocity profiles for anisotropic models show significant flattening toward the tidal radius compared with the isotropic model. This is caused by the fact that the tangential velocity dispersion becomes dominant at the outer parts. Such a behaviour of velocity profile appears to be consistent with the observed profiles of the collapsed cluster M15.     

Tidal effects of the galaxy on a globular cluster

Tidal effects of a Schmidt (1965) model galaxy on a typical globular cluster moving in an orbit along the axis of symmetry of the spheroid is studied under the impulsive approximation. Ann-body simulation is made for comparison. Results show that in both cases, a high concentration cluster gains enough energy to be totally dissolved within a distance of 2 kpc. A significant mass loss occurs as the cluster approaches the Roche distance.     

Time variation of ellipticity of globular clusters in the Large Magellanic Cloud

Dynamical evolution of globular clusters in the Large Magellanic Cloud (LMC) is investigated by means of N-body simulations; particular attention is paid to time evolution in the ellipticitical figure of globular clusters. The simulations were started with a binary globular cluster. It merged into a single cluster with ellipticity of about 0.3. The simulations were continued until the cluster became rounder due to the effects of two body relaxation and of tidal field of LMC. It is found that the outward angular momentum transport due to the gravothermal contraction makes the inner region rounder; the ellipticity at about the initial half-mass radius (r _h) decreases with the e-folding time of 20 relaxation times. On the other hand, the outer region becomes rounder due to the stripping of stars by the tidal field; the ellipticity at about 3r _h decreases with the e-folding time of 80 crossing times therein, though the time scale depends on the direction of the tidal field relative to the spin of the cluster. These two effects are comparable at about the half-mass radius. Taking account of such theoretical results we reanalyzed observed data for the ellipticity at about the half-mass radius of LMC clusters. We estimated the relaxation time and crossing time for each of the observed clusters, from which we calculated the effective time of getting round of the cluster. We plotted the observed ellipticity of the clusters against their non-dimensional age — i.e., the age normalized by the effective time. We found that observed ellipticity distribution is consistent with our picture.     

Resolving the timing problem of the globular clusters orbiting the Fornax dwarf galaxy

We re-investigate the old problem of the survival of the five globular clusters (GCs) orbiting the Fornax dwarf galaxy in both standard and modified Newtonian dynamics (MOND). For the first time in the history of the topic, we use accurate mass models for the Fornax dwarf, obtained through Jeans modelling of the recently published line-of-sight (LOS) velocity dispersion data, and we are also not resigned to circular orbits for the GCs. Previously conceived problems stem from fixing the starting distances of the globulars to be less than half the tidal radius. We relax this constraint since there is absolutely no evidence for it and show that the dark matter (DM) paradigm, with either cusped or cored DM profiles, has no trouble sustaining the orbits of the two least massive GCs for a Hubble time almost regardless of their initial distance from Fornax. The three most massive globulars can remain in orbit as long as their starting distances are marginally outside the tidal radius. The outlook for MOND is also not nearly as bleak as previously reported. Although dynamical friction (DF) inside the tidal radius is far stronger in MOND, outside DF is negligible due to the absence of stars. This allows highly radial orbits to survive, but more importantly circular orbits at distances more than 85 per cent of Fornax's tidal radius to survive indefinitely. The probability of the GCs being on circular orbits at this distance compared with their current projected distances is discussed and shown to be plausible. Finally, if we ignore the presence of the most massive globular (giving it a large LOS distance), we demonstrate that the remaining four globulars can survive within the tidal radius for the Hubble time with perfectly sensible orbits.     

Effects of external tidal field on the evolution of the outer regions of multi-mass star clusters

We present N -body simulations (including an initial mass function) of globular clusters in the Galaxy in order to study effects of the tidal field systematically on the properties of the outer parts of globular clusters. Using nbody6 , which correctly takes into account the two-body relaxation, we investigate the development of tidal tails of globular clusters in the Galactic tidal field. For simplicity, we have employed only the spherical components (bulge and halo) of the Galaxy, and ignored the effects of stellar evolution which could have been important in the very early phase of the cluster evolution. The total number of stars in our simulations is about 20 000, which is much smaller than the realistic number of stars. All simulations had been done for several orbital periods in order to understand the development of the tidal tails. In our scaled-down models, the relaxation time is sufficiently short to show the mass segregation effect, but we did not go far enough to see the core collapse, and the fraction of stars lost from the cluster at the end of the simulations is only ∼10 per cent. The radial distribution of extra-tidal stars can be described by a power law with a slope around −3 in surface density. The directions of tidal tails are determined by the orbits and locations of the clusters. We find that the length of tidal tails increases towards the apogalacticon and decreases towards the perigalacticon. This is an anti-correlation with the strength of the tidal field, caused by the fact that the time-scale for the stars to respond to the potential is similar to the orbital time-scale of the cluster. The escape of stars in the tidal tails towards the pericentre could be another reason for the decrease of the length of tidal tails. We find that the rotational angular velocity of tidally induced clusters shows quite different behaviour from that of initially rotating clusters.     

Test computations on the dynamical evolution of star clusters

Test calculations have been carried out on the evolution of star clusters using the fluiddynamical method devised by Larson (1970a). Large systems of stars have been considered with specific concern with globular clusters. With reference to the analogous standard model by Larson, the influence of varying in turn the various free parameters (cluster mass, star mass, tidal radius, mass concentration of the initial model) has been studied for the results. Furthermore, the partial release of some simplifying assumptions with regard to the relaxation time and distribution of the target stars has been considered. The change of the structural properties is discussed, and the variation of the evolutionary time scale is outlined. An indicative agreement of the results obtained here with structural properties of globular clusters as deduced from previous theoretical models is pointed out.     

Monte Carlo simulations of star clusters – III. A million-body star cluster

A revision of Stodółkiewicz's Monte Carlo code is used to simulate the evolution of million-body star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. All models consist of 1 000 000 stars. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods. During the evolution, the initial mass function (IMF) changes significantly. The local mass function around the half-mass radius closely resembles the actual global mass function. At the late stages of evolution, the mass of the evolved stars inside the core can be as high as 97 per cent of the total mass in this region. For the whole system, the evolved stars can compose up to 75 per cent of the total mass. The evolution of cluster anisotropy strongly depends on initial cluster concentration, IMF and the strength of the tidal field. The results presented are the first step in the direction of simulating the evolution of real globular clusters by means of the Monte Carlo method.     

Spectroscopy of globular cluster candidates in the Sculptor group galaxies NGC 253 and 55

We have obtained spectra for 103 published globular cluster candidates in the Sculptor group galaxies NGC 253 and 55. On the basis of radial velocities and digitized plate images, 14 globular clusters are identified in NGC 253 and one probable globular cluster is identified in NGC 55. The majority of the objects in the sample appear to be background galaxies. We have obtained and analysed COSMOS plate scans of NGC 253 and 55 and use these along with the spectroscopically identified clusters to define new samples of globular cluster candidates in the two galaxies which should have reduced contamination.     

Globular cluster systems in low-luminosity early-type galaxies near the Fornax cluster centre

We present a photometric study of the globular cluster systems (GCSs) of the Fornax cluster galaxies NGC 1374, NGC 1379 and NGC 1387. The data consist of images from the wide-field MOSAIC imager of the Cerro Tololo Inter-American Observatory (CTIO) 4-m telescope, obtained with Washington C and Kron–Cousins R filters. The images cover a field of  36 × 36 arcmin²  , corresponding to  200 × 200 kpc²  at the Fornax distance. Two of the galaxies, NGC 1374 and NGC 1379, are low-luminosity ellipticals while NGC 1387 is a low-luminosity lenticular. Their cluster systems are still embedded in the cluster system of NGC 1399. Therefore, the use of a large field is crucial and some differences to previous work can be explained by this. The colour distributions of all GCSs are bimodal. NGC 1387 presents a particularly distinct separation between red and blue clusters and an overproportionally large population of red clusters. The radial distribution is different for blue and red clusters, red clusters being more concentrated towards the respective galaxies. The different colour and radial distributions point to the existence of two globular cluster subpopulations in these galaxies. Specific frequencies are in the range   S_N = 1.4–2.4  , smaller than the typical values for elliptical galaxies. These galaxies might have suffered tidal stripping of blue globular clusters by NGC 1399.     

The mass and anisotropy profiles of galaxy clusters from the projected phase-space density: testing the method on simulated data

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase-space density, which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters, we analyse mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard Λ cold dark matter model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 per cent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case, it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.     

Scaling of N-body calculations

We report results of collisional N -body simulations aimed at studying the N dependence of the dynamical evolution of star clusters. Our clusters consist of equal-mass stars and are in virial equilibrium. Clusters moving in external tidal fields and clusters limited by a cut-off radius are simulated. Our main focus is to study the dependence of the lifetimes of the clusters on the number of cluster stars and the chosen escape condition.
We find that star clusters in external tidal fields exhibit a scaling problem in the sense that their lifetimes do not scale with the relaxation time. Isolated clusters show a similar problem if stars are removed only after their distance to the cluster centre exceeds a certain cut-off radius. If stars are removed immediately after their energy exceeds the energy necessary for escape, the scaling problem disappears.
We show that some stars that gain the energy necessary for escape are scattered to lower energies before they can leave the cluster. As the efficiency of this process decreases with increasing particle number, it causes the lifetimes not to scale with the relaxation time. Analytic formulae are derived for the scaling of the lifetimes in the different cases.     

The formation of a bound star cluster: from the Orion nebula cluster to the Pleiades

Direct N -body calculations are presented of the formation of Galactic clusters using GasEx , which is a variant of the code Nbody6 . The calculations focus on the possible evolution of the Orion nebula cluster (ONC) by assuming that the embedded OB stars explosively drove out 2/3 of its mass in the form of gas about 0.4 Myr ago. A bound cluster forms readily and survives for 150 Myr despite additional mass loss from the large number of massive stars, and the Galactic tidal field. This is the very first time that cluster formation is obtained under such realistic conditions. The cluster contains about 1/3 of the initial 10⁴ stars, and resembles the Pleiades cluster to a remarkable degree, implying that an ONC-like cluster may have been a precursor of the Pleiades. This scenario predicts the present expansion velocity of the ONC, which will be measurable by upcoming astrometric space missions. These missions should also detect the original Pleiades members as an associated expanding young Galactic-field subpopulation. The results arrived at here suggest that Galactic clusters form as the nuclei of expanding OB associations.
The results have wide implications, also for the formation of globular clusters and the Galactic-field and halo stellar populations. In view of this, the distribution of binary orbital periods and the mass function within and outside the model ONC and Pleiades is quantified, finding consistency with observational constraints. Advanced mass segregation is evident in one of the ONC models. The calculations show that the primordial binary population of both clusters could have been much the same as is observed in the Taurus–Auriga star-forming region. The computations also demonstrate that the binary proportion of brown dwarfs is depleted significantly for all periods, whereas massive stars attain a high binary fraction.     

The distance to Galactic globular clusters through RR Lyrae pulsational properties

By adopting the same approach outlined by De Santis & Cassisi, we evaluate the absolute bolometric magnitude of the zero-age horizontal branch (ZAHB) at the level of the RR Lyrae variable instability strip in selected Galactic globular clusters. This allows us to estimate the ZAHB absolute visual magnitude for these clusters and to investigate its dependence on the cluster metallicity. The derived M _V (ZAHB)–[Fe/H] relation, corrected in order to account for the luminosity difference between the ZAHB and the mean RR Lyrae magnitude, has been compared with some of the most recent empirical determinations in this field, such as the one provided by Baade–Wesselink analyses, RR Lyrae periods, Hipparcos data for field variables and main-sequence fitting based on Hipparcos parallaxes for field subdwarfs. As a result, our relation provides a clear support to the 'long' distance scale. We discuss also another method for measuring the distance to Galactic globular clusters. This method is quite similar to the one adopted for estimating the absolute bolometric magnitude of the ZAHB but it relies only on the pulsational properties of the Lyrae variables in each cluster. The reliability and accuracy of this method have been tested by applying it to a sample of globular clusters for which, owing to the morphology of their horizontal branch (HB), the use of the commonly adopted ZAHB fitting is a risky procedure. We notice that the two approaches for deriving the cluster distance modulus provide consistent results when applied to globular clusters, the RR Lyrae instability strip is well populated. As the adopted method relies on theoretical predictions on both the fundamental pulsational equation and the allowed mass range for fundamental pulsators, we give an estimate of the error affecting present results, owing to systematic uncertainties in the adopted theoretical framework.     

Warped discs and the directional stability of jets in active galactic nuclei

The radial distribution of globular clusters in galaxies is always less peaked to the centre than that of the halo stars. Extending previous work to a sample of Hubble Space Telescope globular cluster systems in ellipticals, we evaluate the number of clusters potentially lost to the galactic centre as the integrals of the difference between the observed globular cluster system distribution and the underlying halo light profile. In the sample of galaxies examined it is found that the initial populations of globular clusters may have been ∼30 per cent to 50 per cent richer than now. If these 'missing' globular clusters have decayed and have been partly destroyed in the very central galactic zones, they have carried there a significant quantity of mass that, plausibly, contributed to the formation and feeding of a massive object therein. It is relevant to note that the observed correlation between the core radius of the globular cluster system and the parent galaxy luminosity can be interpreted as a result of evolution.     

Tidal tails of star clusters

Based on recent findings of a formation mechanism of substructure in tidal tails by Küpper et al., we investigate a more comprehensive set of N -body models of star clusters on orbits about a Milky Way like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for the overdensities to build up.
The distance of the overdensities along the tidal tail from the cluster centre depends for circular orbits only on the mass of the cluster and the strength of the tidal field, and therefore decreases monotonically with time, while for eccentric orbits the orbital motion influences the distance, causing a periodic compression and stretching of the tails and making the distance oscillate with time. We provide an approximation for estimating the distance of the overdensities in this case.
We describe an additional type of overdensity which arises in extended tidal tails of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or a disc shock is not the direct origin of an overdensity within a tidal tail. Escape due to such tidal perturbations does not take place immediately after the perturbation but is rather delayed and spread over the orbit of the cluster. All observable overdensities are therefore of the mentioned two types. In particular, we note that substructured tidal tails do not imply the existence of dark matter substructures in the haloes of galaxies.     

Quantitative analysis of clumps in the tidal tails of star clusters

Tidal tails of star clusters are not homogeneous but show well-defined clumps in observations as well as in numerical simulations. Recently, an epicyclic theory for the formation of these clumps was presented. A quantitative analysis was still missing. We present a quantitative derivation of the angular momentum and energy distribution of escaping stars from a star cluster in the tidal field of the Milky Way and derive the connection to the position and width of the clumps. For the numerical realization we use star-by-star N -body simulations. We find a very good agreement of theory and models. We show that the radial offset of the tidal arms scales with the tidal radius, which is a function of cluster mass and the rotation curve at the cluster orbit. The mean radial offset is 2.77 times the tidal radius in the outer disc. Near the Galactic Centre the circumstances are more complicated, but to lowest order the theory still applies. We have also measured the Jacobi energy distribution of bound stars and showed that there is a large fraction of stars (about 35 per cent) above the critical Jacobi energy at all times, which can potentially leave the cluster. This is a hint that the mass loss is dominated by a self-regulating process of increasing Jacobi energy due to the weakening of the potential well of the star cluster, which is induced by the mass loss itself.     

Density profiles and luminosity functions of the globular clusters M3 and M92

Tautenburg Schmidt plates with the globular clusters M3 and M92 were measured on the APM facility in Cambridge, mainly for astrometric purposes. The data were used to derive the radial density profiles and the luminosity functions of these clusters. There is a good agreement of our results with formerly ones derived by other authors. In both clusters effects of mass segregation are observed.     

Stellar kinematics in the centers of globular clusters

The high central stellar densities in globular clusters provide a unique environment to study the fundamental dynamical process of two‐body relaxation. This process is the main driver of the dynamical evolution in the center of a globular cluster and has a profound effect on the structure of the cluster and on its stellar environment. We have obtained stellar absorption line spectra with STIS to measure the radial velocities of individual stars in the crowded center of the globular cluster M15. These data increase the number of stars with known radial velocities within the central arcsec by a factor of about three and significantly improve the constraints on the mass distribution in M15. The data provide the most detailed look of the central structure of any globular cluster and show that there is a compact dark central mass component. Similar studies using ground based facilities can be efficiently performed by employing Integral Field Units. We have started a project to better constrain the central mass density in the globular cluster M3 using the GMOS‐IFU on Gemini North. The data will also allow us to better understand the central rotation which is neither explained nor predicted by any globular cluster model. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mass segregation in globular clusters

Photoelectric aperture-photometry ofω Cen inU, B, V, R andI bands has established that the cluster is bluer between 2 arcmin and 4 arcmin from the centre, than it is elsewhere. The difference inB - I colour between the centre and this blue zone is ≃ 0.45 mag. The core radius is found to be dependant on the wavelength band chosen for observation, the smallest core radius being for theI band. Equidensitometry ofω Cen inB, V and infrared bands shows a wavelength dependence with the cluster being nearly spherical in the infrared band. It shows a maximum ellipticity around 3 arcmin from the cluster centre. The blue contribution in this zone comes from both a diffuse background of unresolved stars and an increase in the relative abundance of horizontal branch (HB) stars. The similarity between the diffuse background and the HB stars has been demonstrated. A photographic subtraction technique is used to study the distribution of HB stars in the cluster. Results of equidensitometry of the cluster 47 Tuc, obtained in the present study, are compared with the earlier results of photoelectric photometry. Here too an increase in ellipticity is associated with an increase in the blueness of the cluster. All globular clusters studied so far for ellipticity show a similarity in the dependence of the ellipticity on the distance from the centre. The ellipticity has small values near the centre and in the outer regions, with the maximum value in between. We suggest that the red stars in globular clusters have a nearly spherical distribution. The blue stars form a bulge around the core with a more elliptical distribution and a different orientation. A similarity between the ellipticity aspects of both the globular clusters and rotation in the nucleus of M 31 is pointed out; the M 31 nucleus may thus show a bluer colour and smaller UV excess around the region where the rotation curve shows a peak. An erratum to this article is available at .     

球状星团的天体测量(英)

球状星团是银河系中最年老的天体之一，是储存着银河系早期演化珍贵信息的“化石”．球状星团的天体测量，主要包括球状星因天区内恒星相对自行的测定，并由这些相对自行数据采用适当方法定出星团的绝对自行，或者直接测定绝对自行．利用这些自行数据，或者进一步与测光和视向速度数据结合，可以开展与球状星团的距离、运动、动力学状况、质量、年龄、演化等等以及银河系的结构和演化等有关的一系列重要的研究．在本文中对本世纪70年代中期以来在球状星团相对自行测定和成员概率估计、内部运动检测、绝对自行测定和空间运动研究这三方面取得的成果和进展以及存在的问题作了评述．