We present the results of the analysis of deep photometric data of 32 Galactic globular clusters. We analysed 69 parallel field images observed with the Wide Field Channel of the Advanced Camera for Surveys of the Hubble Space Telescope which complemented the already available photometry from the globular cluster treasury project covering the central regions of these clusters. This unprecedented data set has been used to calculate the relative fraction of stars at different masses (i.e. the present-day mass function) in these clusters by comparing the observed distribution of stars along the cluster main sequence and across the analysed field of view with the prediction of multimass dynamical models. For a subsample of 31 clusters, we were able to obtain also the half-mass radii, mass-to-light ratios and the mass fraction of dark remnants using available radial velocity information. We found that the majority of globular clusters have single power law mass functions $F(m) \propto m^\alpha$ with slopes $\alpha>-1$ in the mass range $0.2<m/\text{M}_{\odot}<0.8$. By exploring the correlations between the structural/dynamical and orbital parameters, we confirm the tight anticorrelation between the mass function slopes and the half-mass relaxation times already reported in previous works, and possible second-order dependence on the cluster metallicity. This might indicate the relative importance of both initial conditions and evolutionary effects on the present-day shape of the mass function.

Observed mass-to-light ratios (M/L) of metal-rich globular clusters (GCs) disagree with theoretical predictions. This discrepancy is of fundamental importance since stellar population models provide the stellar masses that underpin most of extragalactic astronomy, near and far. We have derived radial velocities for 1,622 stars located in the centres of 59 Milky Way GCs - twelve of which have no previous kinematic information - using integral-field unit data from the WAGGS project. Using N-body models, we then determine dynamical masses and M/L ratios for the studied clusters. Our sample includes NGC 6528 and NGC 6553, which extend the metallicity range of GCs with measured M/L up to [Fe/H] ~ -0.1 dex. We find that metal-rich clusters have M/L more than two times lower than what is predicted by simple stellar population models. This confirms that the discrepant M/L-[Fe/H] relation remains a serious concern. We explore how our findings relate to previous observations, and the potential causes for the divergence, which we conclude is most likely due to dynamical effects.

The internal dynamics of multiple stellar populations in Globular Clusters (GCs) provides unique constraints on the physical processes responsible for their formation. Specifically, the present-day kinematics of cluster stars, such as rotation and velocity-dispersion, could be related to the initial configuration of the system. In recent work, we provided the first study of the kinematics of different stellar populations in NGC$\,$0104 over a large field of view in the plane of the sky, exploiting Gaia Data Release 2 (DR2) proper motions combined with multi-band ground-based photometry. In this paper, we combine Gaia DR2 proper motions with Very Large Telescope radial velocities to investigate the kinematics along the line of sight and in the plane of the sky of multiple populations in seven GCs, namely NGC$\,$0104, NGC$\,$0288, NGC$\,$5904, NGC$\,$6121, NGC$\,$6254, NGC$\,$6752 and NGC$\,$6838. Among the analyzed clusters only NGC$\,$0104 and NGC$\,$5904 show significant rotation. Separating our sample into two groups of first- and second-population stars (1P and 2P) we find that overall these two populations exhibit a similar rotation pattern in NGC$\,$0104. However, some hints of different rotation are observed in the external regions of this cluster. Interestingly, 1P and 2P stars in NGC$\,$5904 show different rotation curves, with distinct phases and such difference is significant at the $\sim$2.5-$\sigma$ level. The analysis of the velocity-dispersion profiles of multiple populations confirms that 2P stars of NGC$\,$0104 show stronger anisotropy than the 1P.

The 'chromosome map' diagram (ChM) proved a successful tool to identify and characterize multiple populations (MPs) in 59 Galactic Globular Clusters (GCs). Here, we construct ChMs for 11 GCs of both Magellanic Clouds (MCs) and with different ages to compare MPs in Galactic and extra-Galactic environments, and explore whether this phenomenon is universal through 'place' and 'time'. MPs are detected in five clusters. The fractions of 1G stars, ranging from about 50% to more than 80%, are significantly higher than those observed in Galactic GCs with similar present-day masses. By considering both Galactic and MC clusters, the fraction of 1G stars exhibits: (i) a strong anti-correlation with the present-day mass, and (ii) with the present-day mass of 2G stars; (iii) a mild anti-correlation with 1G present-day mass. All Galactic clusters without MPs have initial masses smaller than ~1.5 10^5 solar masses but a mass threshold governing the occurrence of MPs seems challenged by massive simple-population MC GCs; (iv) Milky Way clusters with large perigalactic distances typically host larger fractions of 1G stars, but the difference disappears when we use initial cluster masses. These facts are consistent with a scenario where the stars lost by GCs mostly belong to the 1G. By exploiting recent work based on Gaia, half of the known Type II GCs appear clustered in a distinct region of the integral of motions space, thus suggesting a common progenitor galaxy. Except for these Type II GCs,we do not find any significant difference in the MPs between clusters associated with different progenitors.

We collected radial velocities of more than 50.000 individual stars in 156 Galactic globular clusters (GGC) and matched them with HST photometry and Gaia DR2 proper motions. This allowed us to derive the GGC's mean proper motions and space velocities. By fitting a large set of N-body simulations to their velocity dispersion and surface density profiles, combined with new measurements of their internal radially dependent mass functions, we have determined their present-day masses and structural parameters, and for 144 GGCs their internal kinematics. We also derive the initial cluster masses by calculating the cluster orbits backwards in time applying suitable recipes to account for mass-loss and dynamical friction. The new fundamental parameters of GGCs are publicly available via an online database, which will regularly be updated.

We compare the results of a large grid of N-body simulations with the surface brightness and velocity dispersion profiles of the globular clusters $\omega$ Cen and NGC 6624. Our models include clusters with varying stellar-mass black hole retention fractions and varying masses of a central intermediate-mass black hole (IMBH). We find that an $\sim 45,000$ M$_\odot$ IMBH, whose presence has been suggested based on the measured velocity dispersion profile of $\omega$ Cen, predicts the existence of about 20 fast-moving, $m>0.5$ M$_\odot$ main-sequence stars with a (1D) velocity $v>60$ km/sec in the central 20 arcsec of $\omega$ Cen. However no such star is present in the HST/ACS proper motion catalogue of Bellini et al. (2017), strongly ruling out the presence of a massive IMBH in the core of $\omega$ Cen. Instead, we find that all available data can be fitted by a model that contains 4.6% of the mass of $\omega$ Cen in a centrally concentrated cluster of stellar-mass black holes. We show that this mass fraction in stellar-mass BHs is compatible with the predictions of stellar evolution models of massive stars. We also compare our grid of $N$-body simulations with NGC 6624, a cluster recently claimed to harbor a 20,000 M$_\odot$ black hole based on timing observations of millisecond pulsars. However, we find that models with $M_{IMBH}>1,000$ M$_\odot$ IMBHs are incompatible with the observed velocity dispersion and surface brightness profile of NGC 6624,ruling out the presence of a massive IMBH in this cluster. Models without an IMBH provide again an excellent fit to NGC 6624.

In this paper, we analyze stellar proper motions in the core of the globular cluster 47 Tucanae to explore the possibility of an intermediate-mass black hole (IMBH) influence on the stellar dynamics. Our use of short-wavelength photometry affords us an exceedingly clear view of stellar motions into the very center of the crowded core, yielding proper motions for $>$50,000 stars in the central 2'. We model the velocity dispersion profile of the cluster using an isotropic Jeans model. The density distribution is taken as a central IMBH point mass added to a combination of King templates. We individually model the general low-mass cluster objects (main sequence/giant stars), as well as the concentrated populations of heavy binary systems and dark stellar remnants. Using unbinned likelihood model fitting, we find that the inclusion of the concentrated populations in our model plays a crucial role in fitting for an IMBH mass. The concentrated binaries and stellar-mass black holes (BHs) produce a sufficient velocity dispersion signal in the core so as to make an IMBH unnecessary to fit the observations. We additionally determine that a stellar-mass BH retention fraction of $\gtrsim 8.5\%$ becomes incompatible with our observed velocities in the core.

We derived the three-dimensional velocities of individual stars in a sample of 62 Galactic globular clusters using proper motions from the second data release of the Gaia mission together with the most comprehensive set of line-of-sight velocities with the aim of investigating the rotation pattern of these stellar systems. We detect the unambiguous signal of rotation in 15 clusters at amplitudes which are well above the level of random and systematic errors. For these clusters, we derived the position and inclination angle of the rotation axis with respect to the line of sight and the overall contribution of rotation to the total kinetic energy budget. The rotation strengths are weakly correlated with the half-mass radius, the relaxation time and anticorrelated with the destruction rate, while no significant alignment of the rotation axes with the orbital poles has been observed. This evidence points toward a primordial origin of the systemic rotation in these stellar systems.

We present new deep imaging of the central regions of the remote globular cluster NGC 2419, obtained with the F343N and F336W filters of HST/WFC3. The new data are combined with archival imaging to constrain nitrogen and helium abundance variations within the cluster. We find a clearly bimodal distribution of the nitrogen-sensitive F336W-F343N colours of red giants, from which we estimate that about 55% of the giants belong to a population with about normal (field-like) nitrogen abundances (P1), while the remaining 45% belong to a nitrogen-rich population (P2). On average, the P2 stars are more He-rich than the P1 stars, with an estimated mean difference of Delta Y = 0.05, but the P2 stars exhibit a significant spread in He content and some may reach Delta Y = 0.13. A smaller He spread may be present also for the P1 stars. Additionally, stars with spectroscopically determined low [Mg/Fe] ratios ([Mg/Fe]<0) are generally associated with P2. We find the P2 stars to be slightly more centrally concentrated in NGC 2419 with a projected half-number radius of about 10% less than for the P1 stars, but the difference is not highly significant (p=0.05). We find evidence of rotation for the P1 stars, whereas the results are inconclusive for the P2 stars, which are consistent with no rotation as well as the same average rotation found for the P1 stars. Because of the long relaxation time scale of NGC 2419, the radial trends and kinematic properties of the populations are expected to be relatively unaffected by dynamical evolution. Hence, they provide constraints on formation scenarios for multiple populations, which must account not only for the presence of He spreads within sub-populations identified via CNO variations, but also for the relatively modest differences in the spatial distributions and kinematics of the populations.

We analyse high spatial resolution near infra-red (NIR) imaging of NGC6384, a Milky Way like galaxy, using ARGOS commissioning data at the Large Binocular Telescope (LBT). ARGOS provides a stable PSF$_{\rm FWHM}\!=\!0.2"\!-\!0.3"$ AO correction of the ground layer across the LUCI2 NIR camera $4'\!\times4'$ field by using six laser guide stars (three per telescope) and a natural guide star for tip-tilt sensing and guiding. Enabled by this high spatial resolution we analyse the structure of the nuclear star cluster (NSC) and the central kiloparsec of NGC6384. We find via 2D modelling that the NSC ($r_{\rm eff}\!\simeq\!10$pc) is surrounded by a small ($r_{\rm eff}\!\simeq\!100$pc) and a larger Sersi\'c ($r_{\rm eff}\!\simeq\!400$pc), all embedded within the NGC\,6384 large-scale boxy/X-shaped bulge and disk. This proof-of-concept study shows that with the high spatial resolution achieved by ground-layer AO we can push such analysis to distances previously only accessible from space. SED-fitting to the NIR and optical HST photometry allowed to leverage the age-metallicity-extinction degeneracies and derive the effective NSC properties of an young to old population mass ratio of $8\%$ with ${\cal M}_{\rm\star,old}\!\simeq\!3.5\times10^7M_\odot$, Age$_{\rm old,\ young}\!=\!10.9\pm1.3$Gyr and 226\,Myr $\pm62\%$, metallicity [M/H]$=\!-0.11\pm0.16$ and $0.33\pm39\%$dex, and $E(B\!-\!V)\!=\!0.63$ and 1.44mag.

The recent discovery of super-massive black holes (SMBHs) in the centers of high-mass ultra compact dwarf galaxies (UCDs) suggests that at least some UCDs are the stripped nuclear star clusters of lower mass galaxies. Tracing these former nuclei of stripped galaxies provides a unique way to track the assembly history of a galaxy or galaxy cluster. In this paper we present a new method to estimate how many UCDs host an SMBH in their center and thus are stripped galaxy nuclei. We revisit the dynamical mass measurements that suggest many UCDs have more mass than expected from stellar population estimates, which recent observations have shown is due to the presence of an SMBH. We revise the stellar population mass estimates using a new empirical relation between the mass-to-light ratio (M/L) and metallicity, and use this to predict which UCDs are most likely to host an SMBH. This enables us to calculate the fraction of UCDs that host SMBHs across their entire luminosity range for the first time. We then apply the SMBH occupation fraction to the observed luminosity function of UCDs and estimate that in the Fornax and Virgo cluster alone there should be $69^{+32}_{-25}$ stripped nuclei with SMBHs. This analysis shows that stripped nuclei with SMBHs are almost as common in clusters as present-day galaxy nuclei. We estimate the local SMBH number density in stripped nuclei to $3-8\times10^{-3}Mpc^{-3}$, which represents a significant fraction (10-40\%) of the SMBH density in the local Universe. These SMBHs hidden in stripped nuclei will increase expected event rates for tidal disruption events and SMBH-SMBH and SMBH-BH mergers. The existence of numerous stripped nuclei with SMBHs are a direct consequence of hierarchical galaxy formation, but until now their impact on the SMBH density had not been quantified.

As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution-function (DF)-based models, Jeans models, and a grid of N-body models) by applying them to mock observations from a star-by-star N-body simulation of the GC M 4 by Heggie. The mocks mimic existing and anticipated data for GCs: surface brightness or number density profiles, local stellar mass functions, line-of-sight velocities, and Hubble Space Telescope- and Gaia-like proper motions. We discuss the successes and limitations of the methods. We find that multimass DF-based models, Jeans, and N-body models provide more accurate mass profiles compared to single-mass DF-based models. We highlight complications in fitting the kinematics in the outskirts due to energetically unbound stars associated with the cluster ("potential escapers", not captured by truncated DF models nor by N-body models of clusters in isolation), which can be avoided with DF-based models including potential escapers, or with Jeans models. We discuss ways to account for mass segregation. For example, three-component DF-based models with freedom in their mass function are a simple alternative to avoid the biases of single-mass models (which systematically underestimate the total mass, half-mass radius, and central density), while more realistic multimass DF-based models with freedom in the remnant content represent a promising avenue to infer the total mass and the mass function of remnants.

We have derived the mean proper motions and space velocities of 154 Galactic globular clusters and the velocity dispersion profiles of 141 globular clusters based on a combination of Gaia DR2 proper motions with ground-based line-of-sight velocities. Combining the velocity dispersion profiles derived here with new measurements of the internal mass functions allows us to model the internal kinematics of 144 clusters, more than 90% of the currently known Galactic globular cluster population. We also derive the initial cluster masses by calculating the cluster orbits backwards in time applying suitable recipes to account for mass loss and dynamical friction. We find a correlation between the stellar mass function of a globular cluster and the amount of mass lost from the cluster, pointing to dynamical evolution as one of the mechanisms shaping the mass function of stars in clusters. The mass functions also show strong evidence that globular clusters started with a bottom-light initial mass function. Our simulations show that the currently surviving globular cluster population has lost about 80% of its mass since the time of formation. If globular clusters started from a log-normal mass function, we estimate that the Milky Way contained about 500 globular clusters initially, with a combined mass of about $2.5 \cdot 10^8$ $M_\odot$. For a power-law initial mass function, the initial mass in globular clusters could have been a factor of three higher.

The large columns of dusty gas enshrouding and fuelling star-formation in young, massive stellar clusters may render such systems optically thick to radiation well into the infrared. This raises the prospect that both "direct" radiation pressure produced by absorption of photons leaving stellar surfaces and "indirect" radiation pressure from photons absorbed and then re-emitted by dust grains may be important sources of feedback in such systems. Here we evaluate this possibility by deriving the conditions under which a spheroidal, self-gravitating, mixed gas-star cloud can avoid catastrophic disruption by the combined effects of direct and indirect radiation pressure. We show that radiation pressure sets a maximum star cluster formation efficiency of $\epsilon_{\rm max} \sim 0.9$ at a (very large) gas surface density of $\sim 10^5 M_\odot$ pc$^{-2} (Z_\odot/Z) \simeq 20$ g cm$^{-2} (Z_\odot/Z)$, but that gas clouds above this limit undergo significant radiation-driven expansion during star formation, leading to a maximum stellar surface density very near this value for all star clusters. Data on the central surface mass density of compact stellar systems, while sparse and partly confused by dynamical effects, are broadly consistent with the existence of a metallicity-dependent upper-limit comparable to this value. Our results imply that this limit may preclude the formation of the progenitors of intermediate-mass black holes for systems with $Z \gtrsim 0.2 Z_\odot$.

We have determined masses, stellar mass functions and structural parameters of 112 Milky Way globular clusters by fitting a large set of N-body simulations to their velocity dispersion and surface density profiles. The velocity dispersion profiles were calculated based on a combination of more than 15,000 high-precision radial velocities which we derived from archival ESO/VLT and Keck spectra together with ~20,000 published radial velocities from the literature. Our fits also include the stellar mass functions of the globular clusters, which are available for 47 clusters in our sample, allowing us to self-consistently take the effects of mass segregation and ongoing cluster dissolution into account. We confirm the strong correlation between the global mass functions of globular clusters and their relaxation times recently found by Sollima & Baumgardt (2017). We also find a correlation of the escape velocity from the centre of a globular cluster and the fraction of first generation stars (FG) in the cluster recently derived for 57 globular clusters by Milone et al. (2017), but no correlation between the FG star fraction and the global mass function of a globular cluster. This could indicate that the ability of a globular cluster to keep the wind ejecta from the polluting star(s) is the crucial parameter determining the presence and fraction of second generation stars and not its later dynamical mass loss.

A unique signature for the presence of massive black holes in very dense stellar regions is occasional giant-amplitude outbursts of multiwavelength radiation from tidal disruption and subsequent accretion of stars that make a close approach to the black holes. Previous strong tidal disruption event (TDE) candidates were all associated with the centers of largely isolated galaxies. Here we report the discovery of a luminous X-ray outburst from a massive star cluster at a projected distance of 12.5 kpc from the center of a large lenticular galaxy. The luminosity peaked at ~10^{43} erg/s and decayed systematically over 10 years, approximately following a trend that supports the identification of the event as a TDE. The X-ray spectra were all very soft, with emission confined to be <3.0 keV, and could be described with a standard thermal disk. The disk cooled significantly as the luminosity decreased, a key thermal-state signature often observed in accreting stellar-mass black holes. This thermal-state signature, coupled with very high luminosities, ultrasoft X-ray spectra and the characteristic power-law evolution of the light curve, provides strong evidence that the source contains an intermediate-mass black hole (IMBH) with a mass of a few ten thousand solar mass. This event demonstrates that one of the most effective means to detect IMBHs is through X-ray flares from TDEs in star clusters.

Recent research has been constraining the retention fraction of black holes (BHs) in globular clusters by comparing the degree of mass segregation with $N$-body simulations. They are consistent with an upper limit of the retention fraction being $50\,\%$ or less. In this work, we focus on direct simulations of the dynamics of BHs in star clusters. We aim to constrain the effective distribution of natal kicks that BHs receive during supernova (SN) explosions and to estimate the BH retention fraction. We used the collisional $N$-body code nbody6 to measure the retention fraction of BHs for a given set of parameters, which are: the initial mass of a star cluster, the initial half-mass radius, and $\sigma_\mathrm{BH}$, which sets the effective Maxwellian BH velocity kick distribution. We compare these direct $N$-body models with our analytic estimates and newest observational constraints. The numerical simulations show that for the one-dimensional (1D) velocity kick dispersion $\sigma_\mathrm{BH} < 50\,\mathrm{km\,s^{-1}}$, clusters with radii of 2 pc and that are initially more massive than $5 \times 10^3\,M_\odot$ retain more than $20\,\%$ of BHs within their half-mass radii. Our simple analytic model yields a number of retained BHs that is in good agreement with the $N$-body models. Furthermore, the analytic estimates show that ultra-compact dwarf galaxies (UCDs) should have retained more than $80\,\%$ of their BHs for $\sigma_\mathrm{BH} \leq 190\,\mathrm{km\,s^{-1}}$. Although our models do not contain primordial binaries, in the most compact clusters with $10^3$ stars, we have found evidence of delayed SN explosions producing a surplus of BHs compared to the IMF due to dynamically formed binary stars. These cases do not occur in the more populous or expanded clusters.

The origin of ultracompact dwarfs (UCDs), a class of compact stellar systems discovered two decades ago, still remains a matter of debate. Recent discoveries of central supermassive black holes in UCDs likely inherited from their massive progenitor galaxies provide support for the tidal stripping hypothesis. At the same time, on statistical grounds, some massive UCDs might be representatives of the high luminosity tail of the globular cluster luminosity function. Here we present a detection of a $3.3^{+1.4}_{-1.2}\times10^6\,M_{\odot}$ black hole ($1\sigma$ uncertainty) in the centre of the UCD3 galaxy in the Fornax cluster, that corresponds to 4 per cent of its stellar mass. We performed isotropic Jeans dynamical modelling of UCD3 using internal kinematics derived from adaptive optics assisted observations with the SINFONI spectrograph and seeing limited data collected with the FLAMES spectrograph at the ESO VLT. We rule out the zero black hole mass at the $3\sigma$ confidence level when adopting a mass-to-light ratio inferred from stellar populations. This is the fourth supermassive black hole found in a UCD and the first one in the Fornax cluster. Similarly to other known UCDs that harbour black holes, UCD3 hosts metal rich stars enhanced in $\alpha$-elements that supports the tidal stripping of a massive progenitor as its likely formation scenario. We estimate that up to 80 per cent of luminous UCDs in galaxy clusters host central black holes. This fraction should be lower for UCDs in groups, because their progenitors are more likely to be dwarf galaxies, which do not tend to host central black holes.

We examine the internal properties of the most massive ultracompact dwarf galaxy (UCD), M59-UCD3, by combining adaptive optics assisted near-IR integral field spectroscopy from Gemini/NIFS, and Hubble Space Telescope (HST) imaging. We use the multi-band HST imaging to create a mass model that suggests and accounts for the presence of multiple stellar populations and structural components. We combine these mass models with kinematics measurements from Gemini/NIFS to find a best-fit stellar mass-to-light ratio ($M/L$) and black hole (BH) mass using Jeans Anisotropic Models (JAM), axisymmetric Schwarzschild models, and triaxial Schwarzschild models. The best fit parameters in the JAM and axisymmetric Schwarzschild models have black holes between 2.5 and 5.9 million solar masses. The triaxial Schwarzschild models point toward a similar BH mass, but show a minimum $\chi^2$ at a BH mass of $\sim 0$. Models with a BH in all three techniques provide better fits to the central $V_{rms}$ profiles, and thus we estimate the BH mass to be $4.2^{+2.1}_{-1.7} \times 10^{6}$ M$_\odot$ (estimated 1$\sigma$ uncertainties). We also present deep radio imaging of M59-UCD3 and two other UCDs in Virgo with dynamical BH mass measurements, and compare these to X-ray measurements to check for consistency with the fundamental plane of BH accretion. We detect faint radio emission in M59cO, but find only upper limits for M60-UCD1 and M59-UCD3 despite X-ray detections in both these sources. The BH mass and nuclear light profile of M59-UCD3 suggests it is the tidally stripped remnant of a $\sim$10$^{9-10}$ M$_\odot$ galaxy.

The recent discovery of massive black holes (BHs) in the centers of high-mass ultra compact dwarf galaxies (UCDs) suggests that at least some are the stripped nuclear star clusters of dwarf galaxies. We present the first study that investigates whether such massive BHs, and therefore stripped nuclei, also exist in low-mass ($M<10^{7}M_{\odot}$) UCDs. We constrain the BH masses of two UCDs located in Centaurus A (UCD320 and UCD330) using Jeans modeling of the resolved stellar kinematics from adaptive optics VLT/SINFONI data. No massive BHs are found in either UCD. We find a $3\,\sigma$ upper limit on the central BH mass in UCD\,330 of $M_{\bullet}<1.0\times10^{5}M_{\odot}$, which corresponds to 1.7\% of the total mass. This excludes a high mass fraction BH and would only allow a low-mass BHs similar to those claimed to be detected in Local Group GCs. For UCD320, poorer data quality results in a less constraining $3\,\sigma$ upper limit of $M_{\bullet}<1\times10^{6}M_{\odot}$, which is equal to 37.7\% of the total mass. The dynamical $M/L$ of UCD320 and UCD330 are not inflated compared to predictions from stellar population models. The non-detection of BHs in these low-mass UCDs is consistent with the idea that elevated dynamical $M/L$s do indicate the presence of a substantial BH. Despite not detecting massive BHs, these systems could still be stripped nuclei. The strong rotation ($v/\sigma$ of 0.3 to 0.4) in both UCDs and the two-component light profile in UCD330 support the idea that these UCDs may be stripped nuclei of low-mass galaxies where the BH occupation fraction is not yet known.

The distribution of stars around a massive black hole (MBH) has been addressed in stellar dynamics for the last four decades by a number of authors. Because of its proximity, the centre of the Milky Way is the only observational test case where the stellar distribution can be accurately tested. Past observational work indicated that the brightest giants in the Galactic Centre (GC) may show a density deficit around the central black hole, not a cusp-like distribution, while we theoretically expect the presence of a stellar cusp. We here present a solution to this long-standing problem. We performed direct-summation $N-$body simulations of star clusters around massive black holes and compared the results of our simulations with new observational data of the GC's nuclear cluster. We find that after a Hubble time, the distribution of bright stars as well as the diffuse light follow power-law distributions in projection with slopes of $\Gamma \approx 0.3$ in our simulations. This is in excellent agreement with what is seen in star counts and in the distribution of the diffuse stellar light extracted from adaptive-optics (AO) assisted near-infrared observations of the GC. Our simulations also confirm that there exists a missing giant star population within a projected radius of a few arcsec around Sgr A*. Such a depletion of giant stars in the innermost 0.1 pc could be explained by a previously present gaseous disc and collisions, which means that a stellar cusp would also be present at the innermost radii, but in the form of degenerate compact cores.

This is the second of three papers that search for the predicted stellar cusp around the Milky Way's central black hole, Sagittarius A*, with new data and methods. We aim to infer the distribution of the faintest stellar population currently accessible through observations around Sagittarius A*. We use adaptive optics assisted high angular resolution images obtained with the NACO instrument at the ESO VLT. Through optimised PSF fitting we remove the light from all detected stars above a given magnitude limit. Subsequently we analyse the remaining, diffuse light density. The analysed diffuse light arises from sub-giant and main-sequence stars with KS ~ 19 - 20 with masses of 1 - 2 Msol . These stars can be old enough to be dynamically relaxed. The observed power-law profile and its slope are consistent with the existence of a relaxed stellar cusp around the Milky Way's central black hole. We find that a Nuker law provides an adequate description of the nuclear cluster's intrinsic shape (assuming spherical symmetry). The 3D power-law slope near Sgr A* is \gamma = 1.23 +- 0.05. At a distance of 0.01 pc from the black hole, we estimate a stellar mass density of 2.3 +- 0.3 x 10^7 Msol pc^-3 and a total enclosed stellar mass of 180 +- 20 Msol. These estimates assume a constant mass-to-light ratio and do not take stellar remnants into account. The fact that no cusp is observed for bright (Ks 16) giant stars at projected distances of roughly 0.1-0.3 pc implies that some mechanism has altered their appearance or distribution.

(abridged) In this paper we revisit the problem of inferring the innermost structure of the Milky Way's nuclear star cluster via star counts, to clarify whether it displays a core or a cusp around the central black hole. Through image stacking and improved PSF fitting we push the completeness limit about one magnitude deeper than in previous, comparable work. Contrary to previous work, we analyse the stellar density in well-defined magnitude ranges in order to be able to constrain stellar masses and ages. The RC and brighter giant stars display a core-like surface density profile within a projected radius R<0.3 pc of the central black hole, in agreement with previous studies, but show a cusp-like surface density distribution at larger R. The surface density of the fainter stars can be described well by a single power-law at R<2 pc. The cusp-like profile of the faint stars persists even if we take into account the possible contamination of stars in this brightness range by young pre-main sequence stars. The data are inconsistent with a core-profile for the faint stars.Finally, we show that a 3D Nuker law provides a very good description of the cluster structure. We conclude that the observed stellar density at the Galactic Centre, as it can be inferred with current instruments, is consistent with the existence of a stellar cusp around the Milky Way's central black hole, Sgr A*. This cusp is well developed inside the influence radius of about 3 pc of Sgr A* and can be described by a single three-dimensional power-law with an exponent gamma=1.23+-0.05. The apparent lack of RC stars and brighter giants at projected distances of R < 0.3 pc (R<8") of the massive black hole may indicate that some mechanism has altered their distribution or intrinsic luminosity.

In this paper we compare the mass function slopes of Galactic globular clusters recently determined by Sollima & Baumgardt (2017) with a set of dedicated N-body simulations of star clusters containing between 65,000 to 200,000 stars. We study clusters starting with a range of initial mass functions (IMFs), black hole retention fractions and orbital parameters in the parent galaxy. We find that the present-day mass functions of globular clusters agree well with those expected for star clusters starting with Kroupa or Chabrier IMFs, and are incompatible with clusters starting with single power-law mass functions for the low-mass stars. The amount of mass segregation seen in the globular clusters studied by Sollima & Baumgardt (2017) can be fully explained by two-body relaxation driven mass segregation from initially unsegregated star clusters. Based on the present-day global mass functions, we expect that a typical globular cluster in our sample has lost about 75% of its mass since formation, while the most evolved clusters have already lost more than 90% of their initial mass and should dissolve within the next 1 to 2 Gyr. Most clusters studied by Sollima & Baumgardt also show a large difference between their central and global MF slopes, implying that the majority of Galactic globular clusters is either near or already past core collapse. The strong mass segregation seen in most clusters also implies that only a small fraction of all black holes formed in globular clusters still reside in them.

We have derived the global mass functions of a sample of 35 Galactic globular clusters by comparing deep Hubble Space Telescope photometry with suitable multimass dynamical models. For a subset of 29 clusters with available radial velocity information we were also able to determine dynamical parameters, mass-to-light ratios and the mass fraction of dark remnants. The derived global mass functions are well described by single power-laws in the mass range $0.2 < m/M_\odot < 0.8$ with mass function slopes $\alpha>-1$. Less evolved clusters show deviations from a single-power law, indicating that the original shape of their mass distribution was not a power-law. We find a tight anticorrelation between the present-day mass function slopes and the half-mass relaxation times, which can be understood if clusters started from the same universal IMF and internal dynamical evolution is the main driver in shaping the present-day mass functions. Alternatively, IMF differences correlated with the present-day half-mass relaxation time are needed to explain the observed correlation. The large range of mass function slopes seen for our clusters implies that most globular clusters are dynamically highly evolved, a fact that seems difficult to reconcile with standard estimates for the dynamical evolution of clusters. The mass function slopes also correlate with the dark remnant fractions indicating a preferential retention of massive remnants in clusters subject to high mass-loss rates.

We constrain the possible presence of a central black hole (BH) in the center of the Large Magellanic Cloud (LMC). This requires spectroscopic measurements over an area of order a square degree, due to the poorly known position of the kinematic center. Such measurements are now possible with the impressive field of view of the Multi Unit Spectroscopic Explorer (MUSE) on the ESO Very Large Telescope. We used the Calcium Triplet (~850nm) spectral lines in many short-exposure MUSE pointings to create a two-dimensional integrated-light line-of-sight velocity map from the ~$10^8$ individual spectra, taking care to identify and remove Galactic foreground populations. The data reveal a clear velocity gradient at an unprecedented spatial resolution of 1 arcmin$^{2}$. We fit kinematic models to arrive at a $3\sigma$ upper-mass-limit of $10^{7.1}$ M$_{Sun}$ for any central BH - consistent with the known scaling relations for supermassive black holes and their host systems. This adds to the growing body of knowledge on the presence of BHs in low-mass and dwarf galaxies, and their scaling relations with host-galaxy properties, which can shed light on theories of BH growth and host system interaction.

Our Galaxy hosts the annihilation of a few $\times 10^{43}$ low-energy positrons every second. Radioactive isotopes capable of supplying such positrons are synthesised in stars, stellar remnants, and supernovae. For decades, however, there has been no positive identification of a main stellar positron source leading to suggestions that many positrons originate from exotic sources like the Galaxy's central super-massive black hole or dark matter annihilation. %, but such sources would not explain the recently-detected positron signal from the extended Galactic disk. Here we show that a single type of transient source, deriving from stellar populations of age 3-6 Gyr and yielding ~0.03 $M_\odot$ of the positron emitter $^{44}$Ti, can simultaneously explain the strength and morphology of the Galactic positron annihilation signal and the solar system abundance of the $^{44}$Ti decay product $^{44}$Ca. This transient is likely the merger of two low-mass white dwarfs, observed in external galaxies as the sub-luminous, thermonuclear supernova known as SN1991bg-like.

Intermediate mass black holes play a critical role in understanding the evolutionary connection between stellar mass and super-massive black holes. However, to date the existence of these species of black holes remains ambiguous and their formation process is therefore unknown. It has been long suspected that black holes with masses $10^{2}-10^{4}M_{\odot}$ should form and reside in dense stellar systems. Therefore, dedicated observational campaigns have targeted globular cluster for many decades searching for signatures of these elusive objects. All candidates found in these targeted searches appear radio dim and do not have the X-ray to radio flux ratio predicted by the fundamental plane for accreting black holes. Based on the lack of an electromagnetic counterpart upper limits of $2060 M_{\odot}$ and $470 M_{\odot}$ have been placed on the mass of a putative black hole in 47 Tucanae (NGC 104) from radio and X-ray observations respectively. Here we show there is evidence for a central black hole in 47 Tuc with a mass of M$_{\bullet}\sim2300 M_{\odot}$$_{-850}^{+1500}$ when the dynamical state of the globular cluster is probed with pulsars. The existence of an intermediate mass black hole in the centre of one of the densest clusters with no detectable electromagnetic counterpart suggests that the black hole is not accreting at a sufficient rate and therefore contrary to expectations is gas starved. This intermediate mass black hole might be a member of electromagnetically invisible population of black holes that are the elusive seeds leading to the formation of supermassive black holes in galaxies.

We present the detection of supermassive black holes (BHs) in two Virgo ultracompact dwarf galaxies (UCDs), VUCD3 and M59cO. We use adaptive optics assisted data from the Gemini/NIFS instrument to derive radial velocity dispersion profiles for both objects. Mass models for the two UCDs are created using multi-band Hubble Space Telescope (HST) imaging, including the modeling of mild color gradients seen in both objects. We then find a best-fit stellar mass-to-light ratio ($M/L$) and BH mass by combining the kinematic data and the deprojected stellar mass profile using Jeans Anisotropic Models (JAM). Assuming axisymmetric isotropic Jeans models, we detect BHs in both objects with masses of $4.4^{+2.5}_{-3.0} \times 10^6$ $M_{\odot}$ in VUCD3 and $5.8^{+2.5}_{-2.8} \times 10^6$ $M_{\odot}$ in M59cO (3$\sigma$ uncertainties). The BH mass is degenerate with the anisotropy parameter, $\beta_z$; for the data to be consistent with no BH requires $\beta_z = 0.4$ and $\beta_z = 0.6$ for VUCD3 and M59cO, respectively. Comparing these values with nuclear star clusters shows that while it is possible that these UCDs are highly radially anisotropic, it seems unlikely. These detections constitute the second and third UCDs known to host supermassive BHs. They both have a high fraction of their total mass in their BH; $\sim$13% for VUCD3 and $\sim$18% for M59cO. They also have low best-fit stellar $M/L$s, supporting the proposed scenario that most massive UCDs host high mass fraction BHs. The properties of the BHs and UCDs are consistent with both objects being the tidally stripped remnants of $\sim$10$^9$ M$_\odot$ galaxies.

Palomar 4 is a low-density globular cluster with a current mass $\approx30000 M_{\odot}$ in the outer halo of the Milky Way with a two-body relaxation time of the order of a Hubble time. Yet, it is strongly mass segregated and contains a stellar mass function depleted of low-mass stars. Pal 4 was either born this way or it is a result of extraordinary dynamical evolution. Since two-body relaxation cannot explain these signatures alone, enhanced mass loss through tidal shocking may have had a strong influence on Pal 4. Here, we compute a grid of direct N-body simulations to model Pal 4 on various eccentric orbits within the Milky Way potential to find likely initial conditions that reproduce its observed mass, half-light radius, stellar MF-slope and line-of-sight velocity dispersion. We find that Pal 4 is most likely orbiting on an eccentric orbit with an eccentricity of $e\approx 0.9$ and pericentric distance of $R_p\approx5$ kpc. In this scenario, the required 3D half-mass radius at birth is similar to the average sizes of typical GCs ($R_h\approx4-5$ pc), while its birth mass is about $M_0\approx10^5 M_{\odot}$. We also find a high degree of primordial mass segregation among the cluster stars, which seems to be necessary in every scenario we considered. Thus, using the tidal effect to constrain the perigalactic distance of the orbit of Pal 4, we predict that the proper motion of Pal 4 should be in the range $-0.52\leq\mu_\delta\leq-0.38$ mas\,yr$^{-1}$ and $-0.30\leq\mu_{\alpha\cos\delta}\leq-0.15$ mas\,yr$^{-1}$.

We obtained precise line-of-sight radial velocities of 23 member stars of the remote halo globular cluster Palomar 4 (Pal 4) using the High Resolution Echelle Spectrograph (HIRES) at the Keck I telescope. We also measured the mass function of the cluster down to a limiting magnitude of V~28 mag using archival HST/WFPC2 imaging. We derived the cluster's surface brightness profile based on the WFPC2 data and on broad-band imaging with the Low-Resolution Imaging Spectrometer (LRIS) at the Keck II telescope. We find a mean cluster velocity of 72.55+/-0.22 km/s and a velocity dispersion of 0.87+/-0.18 km/s. The global mass function of the cluster, in the mass range 0.55<=M<=0.85 M_solar, is shallower than a Kroupa mass function and the cluster is significantly depleted in low-mass stars in its center compared to its outskirts. Since the relaxation time of Pal 4 is of the order of a Hubble time, this points to primordial mass segregation in this cluster. Extrapolating the measured mass function towards lower-mass stars and including the contribution of compact remnants, we derive a total cluster mass of 29800 M_solar. For this mass, the measured velocity dispersion is consistent with the expectations of Newtonian dynamics and below the prediction of Modified Newtonian Dynamics (MOND). Pal 4 adds to the growing body of evidence that the dynamics of star clusters in the outer Galactic halo can hardly be explained by MOND.

We have determined the masses and mass-to-light ratios of 50 Galactic globular clusters by comparing their velocity dispersion and surface brightness profiles against a large grid of 900 N-body simulations of star clusters of varying initial concentration, size and central black hole mass fraction. Our models follow the evolution of the clusters under the combined effects of stellar evolution and two-body relaxation allowing us to take the effects of mass segregation and energy equipartition between stars self-consistently into account. For a subset of 16 well observed clusters we also derive their kinematic distances. We find an average mass-to-light ratio of Galactic globular clusters of $<M/L_V>=1.98 \pm 0.03$, which agrees very well with the expected M/L ratio if the initial mass function of the clusters was a standard Kroupa or Chabrier mass function. We do not find evidence for a decrease of the average mass-to-light ratio with metallicity. The surface brightness and velocity dispersion profiles of most globular clusters are incompatible with the presence of intermediate-mass black holes (IMBHs) with more than a few thousand $M_\odot$ in them. The only clear exception is $\omega$ Cen, where the velocity dispersion profile provides strong evidence for the presence of a $\sim$40,000 $M_\odot$ IMBH in the centre of the cluster.

We present MUSE observations of the debated ultra faint stellar system Crater. We spectroscopically confirm 26 member stars of this system via radial velocity measurements. We derive the systematic instrumental velocity uncertainty of MUSE spectra to be 2.27$\rm \,km\,s^{-1}$. This new dataset increases the confirmed member stars of Crater by a factor of 3. One out of three bright blue stars and a fainter blue star just above the main-sequence-turn-off are also found to be likely members of the system. The observations reveal that Crater has a systemic radial velocity of $v_{\rm sys}=148.18^{\rm +1.08}_{\rm -1.15}\rm \,km\,s^{-1}$, whereas the most likely velocity dispersion of this system is $\sigma_{\rm v}=2.04^{\rm +2.19}_{\rm -1.06} \rm \,km\,s^{-1}$. The total dynamical mass of the system, assuming dynamical equilibrium is then $M_{\rm tot}=1.50^{+4.9}_{-1.2}\cdot 10^{\rm 5}M_{\odot}$ implying a mass-to-light ratio of M/L$_{\rm V}$=8.52$^{+28.0}_{-6.5}\, M_{\odot}/L_{\odot}$, which is consistent with a purely baryonic stellar population within its errors and no significant evidence for the presence dark matter was found. We also find evidence for a velocity gradient in the radial velocity distribution. We conclude that our findings strongly support that Crater is a faint intermediate-age outer halo globular cluster and not a dwarf galaxy.

We compare the predictions of a semi-analytic model for ultra-compact dwarf galaxy (UCD) formation by tidal stripping to the observed properties of globular clusters (GCs) and UCDs in the Fornax and Virgo clusters. For Fornax we find the predicted number of stripped nuclei agrees very well with the excess number of GCs$+$UCDs above the GC luminosity function. GCs$+$UCDs with masses $>10^{7.3}$ M$_\odot$ are consistent with being entirely formed by tidal stripping. Stripped nuclei can also account for Virgo UCDs with masses $>10^{7.3}$ M$_\odot$ where numbers are complete by mass. For both Fornax and Virgo, the predicted velocity dispersions and radial distributions of stripped nuclei are consistent with that of UCDs within $\sim$50-100 kpc but disagree at larger distances where dispersions are too high and radial distributions too extended. Stripped nuclei are predicted to have radially biased anisotropies at all radii, agreeing with Virgo UCDs at clustercentric distances larger than 50 kpc. However, ongoing disruption is not included in our model which would cause orbits to become tangentially biased at small radii. We find the predicted metallicities and central black hole masses of stripped nuclei agree well with the metallicities and implied black hole masses of UCDs for masses $>10^{6.5}$ M$_\odot$. The predicted black hole masses also agree well with that of M60-UCD1, the first UCD with a confirmed central black hole. These results suggest that observed GC$+$UCD populations are a combination of genuine GCs and stripped nuclei, with the contribution of stripped nuclei increasing toward the high-mass end.

Many of the scenarios proposed to explain the origin of chemically peculiar stars in globular clusters (GCs) require significant mass-loss ($\ge95\%$) to explain the observed fraction of such stars. In the GCs of the Fornax dwarf galaxy significant mass-loss could be a problem. Larsen et al. (2012) showed that there is a large ratio of GCs to metal-poor field stars in Fornax and about $20-25\%$ of all the stars with ${\rm [Fe/H]}<-2$ belong to the four metal-poor GCs. This imposes an upper limit of $\sim80\%$ mass-loss that could have happened in Fornax GCs. In this paper, we propose a solution to this problem by suggesting that stars can leave the Fornax galaxy. We use a series of $N$-body simulations, to determine the limit of mass-loss from Fornax as a function of the initial orbital radii of GCs and the speed with which stars leave Fornax GCs. We consider a set of cored and cuspy density profiles for Fornax. Our results show that with a cuspy model for Fornax, the fraction of stars which leave the galaxy, can be as high as $\sim90\%$, when the initial orbital radii of GCs are $R=2-3\,\rm{kpc}$ and the initial speed of stars is $v>20\,\rm{km\,s^{-1}}$. We show that such large velocities can be achieved by gas expulsion induced mass-loss but not stellar evolution induced mass-loss. Our results imply that one cannot interpret the metallicity distribution of Fornax field stars as evidence against significant mass-loss in Fornax GCs, if mass-loss is due to gas expulsion.

We have obtained membership probabilities of stars within a field of radius $\sim3^\circ$ around the centre of the open cluster Alpha Persei using proper motions and photometry from the PPMXL and WISE catalogues. We have identified 810 possible stellar members of Alpha Persei. We derived the global and radial present-day mass function (MF) of the cluster and found that they are well matched by two-stage power-law relations with different slopes at different radii. The global MF of Alpha Persei shows a turnover at $m=0.62\,\mathrm{M}_{\odot}$ with low and high-mass slopes of $\alpha_\mathrm{low}=0.50\pm0.09$ ($0.1<m/\mathrm{M}_{\odot}<0.62$) and $\alpha_\mathrm{high}=2.32\pm0.14$ ($0.62\leq m/\mathrm{M}_{\odot}<4.68$) respectively. The high-mass slope of the cluster increases from $2.01$ inside $1\hbox{$.\!\!^\circ$}10$ to $2.63$ outside $2\hbox{$.\!\!^\circ$}2$, whereas the mean stellar mass decreases from $0.95$ to $0.57\,\mathrm{M}_{\odot}$ in the same regions, signifying clear evidence of mass segregation in the cluster. From an examination of the high-quality colour-magnitude data of the cluster and performing a series of Monte Carlo simulations we obtained a binary fraction of $f_{\rm bin}=34\pm12$ percent for stars with $0.70<m/\mathrm{M}_{\odot}<4.68$. This is significantly larger than the observed binary fraction, indicating that this open cluster contains a large population of unresolved binaries. Finally, we corrected the mass-function slopes for the effect of unresolved binaries and found low- and high-mass slopes of $\alpha_\mathrm{low}=0.89\pm0.11$ and $\alpha_\mathrm{high}=2.37\pm0.09$ and a total cluster mass of $352\,\mathrm{M}_{\odot}$.

Outer-halo globular clusters show large half-light radii and flat stellar mass functions, depleted in low-mass stars. Using N-body simulations of globular clusters on eccentric orbits within a Milky Way-like potential, we show how a cluster's half-mass radius and its mass function develop over time. The slope of the central mass function flattens proportionally to the amount of mass a cluster has lost, and the half-mass radius grows to a size proportional to the average strength of the tidal field. The main driver of these processes is mass segregation of dark remnants. We conclude that the extended, depleted clusters observed in the Milky Way must have had small half-mass radii in the past, and that they expanded due to the weak tidal field they spend most of their lifetime in. Moreover, their mass functions must have been steeper in the past but flattened significantly as a cause of mass segregation and tidal mass loss.

We present DEIMOS multi-object spectroscopy (MOS) of 22 star-forming dwarf galaxies located in four gas-rich groups, including six newly-discovered dwarfs. Two of the galaxies are strong tidal dwarf galaxy (TDG) candidates based on our luminosity-metallicity relation definition. We model the rotation curves of these galaxies. Our sample shows low mass-to-light ratios (M/L=0.73$\pm0.39M_\odot/L_\odot$) as expected for young, star-forming dwarfs. One of the galaxies in our sample has an apparently strongly-falling rotation curve, reaching zero rotational velocity outside the turnover radius of $r_{turn}=1.2r_e$. This may be 1) a polar ring galaxy, with a tilted bar within a face-on disk; 2) a kinematic warp. These scenarios are indistinguishable with our current data due to limitations of slit alignment inherent to MOS-mode observations. We consider whether TDGs can be detected based on their tidal radius, beyond which tidal stripping removes kinematic tracers such as H$\alpha$ emission. When the tidal radius is less than about twice the turnover radius, the expected falling rotation curve cannot be reliably measured. This is problematic for as much as half of our sample, and indeed more generally, galaxies in groups like these. Further to this, the H$\alpha$ light that remains must be sufficiently bright to be detected; this is only the case for three (14%) galaxies in our sample. We conclude that the falling rotation curves expected of tidal dwarf galaxies are intrinsically difficult to detect.

We perform a series of direct $N$-body calculations to investigate the effect of residual gas expulsion from the gas-embedded progenitors of present-day globular clusters (GCs) on the stellar mass function (MF). Our models start either tidally filling or underfilling, and either with or without primordial mass segregation. We cover 100 Myr of the evolution of modeled clusters and show that the expulsion of residual gas from initially mass-segregated clusters leads to a significantly shallower slope of the stellar MF in the low- ($m\leq 0.50 M_\odot$) and intermediate-mass ($\simeq 0.50-0.85 M_\odot$) regime. Therefore, the imprint of residual gas expulsion and primordial mass segregation might be visible in the present-day MF. We find that the strength of the external tidal field, as an essential parameter, influences the degree of flattening, such that a primordially mass-segregated tidally-filling cluster with $r_h/r_t\geq 0.1$ shows a strongly depleted MF in the intermediate stellar mass range. Therefore, the shape of the present-day stellar MF in this mass range probes the birth place of clusters in the Galactic environment. We furthermore find that this flattening agrees with the observed correlation between the concentration of a cluster and its MF slope, as found by de Marchi et al.. We show that if the expansion through the residual gas expulsion in primordial mass segregated clusters is the reason for this correlation then GCs most probably formed in strongly fluctuating local tidal fields in the early proto-Milky Way potential, supporting the recent conclusion by Marks \& Kroupa.

We aim to unveil the most massive central cluster black holes in the universe. We present a new search strategy which is based on a black hole mass gain sensitive 'calorimeter' and which links the innermost stellar density profile of a galaxy to the adiabatic growth of its central SMBH. In a first step we convert observationally inferred feedback powers into SMBH growth rates by using reasonable energy conversion efficiency parameters, $\epsilon$. In the main part of this paper we use these black hole growth rates, sorted in logarithmically increasing steps encompassing our whole parameter space, to conduct $N$-Body computations of brightest cluster galaxies with the newly developed MUESLI software. For the initial setup of galaxies we use core-Sersic models in order to account for SMBH scouring. We find that adiabatically driven core re-growth is significant at the highest accretion rates. As a result, the most massive black holes should be located in BCGs with less pronounced cores when compared to the predictions of empirical scaling relations which are usually calibrated in less extreme environments. For efficiency parameters $\epsilon<0.1$, BCGs in the most massive, relaxed and X-ray luminous galaxy clusters might even develop steeply rising density cusps. Finally, we discuss several promising candidates for follow up investigations, among them the nuclear black hole in the Phoenix cluster. Based on our results, it might have a mass of the order of $10^{11} M_\odot$.

We have carried out a large grid of N-body simulations in order to investigate if mass-loss as a result of primordial gas expulsion can be responsible for the large fraction of second generation stars in globular clusters (GCs) with multiple stellar populations (MSPs). Our clusters start with two stellar populations in which $10\%$ of all stars are second generation stars. We simulate clusters with different initial masses, different ratios of the half-mass radius of first to second generation stars, different primordial gas fractions and Galactic tidal fields with varying strength. We then let our clusters undergo primordial gas-loss and obtain their final properties such as mass, half-mass radius and the fraction of second generation stars. Using our N-body grid we then perform a Monte Carlo analysis to constrain the initial masses, radii and required gas expulsion time-scales of GCs with MSPs. Our results can explain the present-day properties of GCs only if (1) a substantial amount of gas was present in the clusters after the formation of second generation stars and (2) gas expulsion time-scales were extremely short ($\lesssim 10^5$ yr). Such short gas expulsion time-scales are in agreement with recent predictions that dark remnants have ejected the primordial gas from globular clusters, and pose a potential problem for the AGB scenario. In addition, our results predict a strong anti-correlation between the number ratio of second-generation stars in GCs and the present-day mass of GCs. So far, the observational data show only a significantly weaker anti-correlation, if any at all.

Recently, two independent groups found very different results when measuring the central velocity dispersion of the galactic globular cluster NGC 6388 with different methods. While L\"utzgendorf et al. (2011) found a rising profile and a high central velocity dispersion (23.3 km/s), measurements obtained by Lanzoni et al. (2013) showed a value 40% lower. The value of the central velocity dispersion has a serious impact on the mass and possible presence of an intermediate-mass black hole at the center of NGC 6388. We use a photometric catalog of NGC 6388 to create a simulated SINFONI and ARGUS dataset. The construction of the IFU data cube is done with different observing conditions reproducing the conditions reported for the original observations as closely as possible. In addition, we produce an N-body realization of a 10^6 M_SUN stellar cluster with the same photometric properties as NGC 6388 to account for unresolved stars. We find that the individual radial velocities, i.e. the measurements from the simulated SINFONI data, are systematically biased towards lower velocity dispersions. The reason is that due to the wings in the point spread function the velocities get biased towards the mean cluster velocity. This study shows that even with AO supported observations, individual radial velocities in crowded fields are likely to be biased. The ARGUS observations do not show this kind of bias but were found to have larger uncertainties than previously obtained. We find a bias towards higher velocity dispersions in the ARGUS pointing when fixing the extreme velocities of the three brightest stars but find those variations are within the determined uncertainties. We rerun Jeans models and fit the kinematic profile with the new uncertainties. This yields a BH mass of M_BH = (2.8 +- 0.4) x 10^4 M_SUN and M/L ratio M/L = (1.6 +- 0.1) M_SUN/L_SUN, consistent with our previous results.

The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and under different levels of tidal stress to derive mass, mass function and degree of anisotropy. We find that masses can be under/over-estimated up to 50% depending on the degree of relaxation reached by the cluster, the available range of observed masses and distances of radial velocity measures from the cluster center and the strength of the tidal field. The mass function slope appears to be better constrainable and less sensitive to model inadequacies unless strongly dynamically evolved clusters and a non-optimal location of the measured luminosity function are considered. The degree and the characteristics of the anisotropy developed in the N-body simulations are not adequately reproduced by popular analytic models and can be detected only if accurate proper motions are available. We show how to reduce the uncertainties in the mass, mass-function and anisotropy estimation and provide predictions for the improvements expected when Gaia proper motions will be available in the near future.

We examine the radial distributions of stellar populations in the globular cluster (GC) M15, using HST/WFC3 photometry of red giants in the nitrogen-sensitive F343N-F555W color. Surprisingly, we find that giants with "primordial" composition (i.e., N abundances similar to those in field stars) are the most centrally concentrated within the WFC3 field. We then combine our WFC3 data with SDSS u, g photometry and find that the trend reverses for radii >1' (3 pc) where the ratio of primordial to N-enhanced giants increases outwards, as already found by Lardo et al. The ratio of primordial to enriched stars thus has a U-shaped dependency on radius with a minimum near the half-light radius. N-body simulations show that mass segregation might produce a trend resembling the observed one, but only if the N-enhanced giants are ~0.25 Mo less massive than the primordial giants, which requires extreme He enhancement (Y~0.40). However, such a large difference in Y is incompatible with the negligible optical color differences between primordial and enriched giants which suggest Delta Y < 0.03 and thus a difference in turn-off mass of Delta M < 0.04 Mo between the different populations. The radial trends in M15 are thus unlikely to be of dynamical origin and presumably reflect initial conditions, a result that challenges all current GC formation scenarios. We note that population gradients in the central regions of GCs remain poorly investigated and may show a more diverse behavior than hitherto thought.

The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. Our group investigated the presence of intermediate-mass black holes for a sample of 10 Galactic globular clusters. We measured the inner kinematic profiles with integral-field spectroscopy and determined masses or upper limits of central black holes in each cluster. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M_bh - sigma) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). We ran N-body simulations of globular clusters containing IMBHs in a tidal field and studied their effects on mass-loss rates and remnant fractions and showed that an IMBH in the center prevents core collapse and ejects massive objects more rapidly. These simulations were further used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighbouring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole.

Ultracompact dwarf galaxies (UCDs) are among the densest stellar systems in the universe. These systems have masses up to 200 million solar masses, but half light radii of just 3-50 parsecs. Dynamical mass estimates show that many UCDs are more massive than expected from their luminosity. It remains unclear whether these high dynamical mass estimates are due to the presence of supermassive black holes or result from a non-standard stellar initial mass function that causes the average stellar mass to be higher than expected. Here we present the detection of a supermassive black hole in a massive UCD. Adaptive optics kinematic data of M60-UCD1 show a central velocity dispersion peak above 100 km/s and modest rotation. Dynamical modeling of these data reveals the presence of a supermassive black hole with mass of 21 million solar masses. This is 15% of the object's total mass. The high black hole mass and mass fraction suggest that M60-UCD1 is the stripped nucleus of a galaxy. Our analysis also shows that M60-UCD1's stellar mass is consistent with its luminosity, implying many other UCDs may also host supermassive black holes. This suggests a substantial population of previously unnoticed supermassive black holes.

We use the Millennium II cosmological simulation combined with the semi-analytic galaxy formation model of Guo et al. (2011) to predict the contribution of galactic nuclei formed by the tidal stripping of nucleated dwarf galaxies to globular cluster (GC) and ultra-compact dwarf galaxy (UCD) populations of galaxies. We follow the merger trees of galaxies in clusters back in time and determine the absolute number and stellar masses of disrupted galaxies. We assume that at all times nuclei have a distribution in nucleus-to-galaxy mass and nucleation fraction of galaxies similar to that observed in the present day universe. Our results show stripped nuclei follow a mass function $N(M) \sim M^{-1.5}$ in the mass range $10^6 < M/M_\odot < 10^8$, significantly flatter than found for globular clusters. The contribution of stripped nuclei will therefore be most important among high-mass GCs and UCDs. For the Milky Way we predict between 1 and 3 star clusters more massive than $10^5 M_\odot$ come from tidally disrupted dwarf galaxies, with the most massive cluster formed having a typical mass of a few times $10^6 M_\odot$, like omega Centauri. For a galaxy cluster with a mass $7 \times 10^{13} M_\odot$, similar to Fornax, we predict $\sim$19 UCDs more massive than $2\times10^6 M_\odot$ and $\sim$9 UCDs more massive than $10^7 M_\odot$ within a projected distance of 300 kpc come from tidally stripped dwarf galaxies. The observed number of UCDs are $\sim$200 and 23, respectively. We conclude that most UCDs in galaxy clusters are probably simply the high mass end of the GC mass function.

We present a new version of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). While previous versions of EMACSS reproduced clusters of single-mass stars, this version models clusters with an evolving stellar content. Stellar evolution dominates early evolution, and leads to: (1) reduction of the mean mass of stars due to the mass loss of high-mass stars; (2) expansion of the half-mass radius; (3) for (nearly) Roche Volume filling clusters, the induced escape of stars. Once sufficient relaxation has occurred (~ 10 relaxation times-scales), clusters reach a second, 'balanced' state whereby the core releases energy as required by the cluster as a whole. In this state: (1) stars escape due to tidal effects faster than before balanced evolution; (2) the half-mass radius expands or contracts depending on the Roche volume filling factor; and (3) the mean mass of stars increases due to the preferential ejection of low-mass stars. We compare the EMACSS results of several cluster properties against N-body simulations of clusters spanning a range of initial number of stars, mass, half-mass radius, and tidal environments, and show that our prescription accurately predicts cluster evolution for this database. Finally, we consider applications for EMACSS, such as studies of galactic globular cluster populations in cosmological simulations.

We use direct $N$-body calculations to study the evolution of the unusually extended outer halo globular cluster Palomar 4 (Pal~4) over its entire lifetime in order to reproduce its observed mass, half-light radius, velocity dispersion and mass function slope at different radii. We find that models evolving on circular orbits, and starting from a non-mass segregated, canonical initial mass function (IMF) can reproduce neither Pal 4's overall mass function slope nor the observed amount of mass segregation. Including either primordial mass segregation or initially flattened IMFs does not reproduce the observed amount of mass segregation and mass function flattening simultaneously. Unresolved binaries cannot reconcile this discrepancy either. We find that only models with both a flattened IMF and primordial segregation are able to fit the observations. The initial (i.e. after gas expulsion) mass and half-mass radius of Pal~4 in this case are about 57000 M${\odot}$ and 10 pc, respectively. This configuration is more extended than most globular clusters we observe, showing that the conditions under which Pal~4 formed must have been significantly different from that of the majority of globular clusters. We discuss possible scenarios for such an unusual configuration of Pal~4 in its early years.

We present the adaptable MUESLI code for investigating dynamics and erosion processes of globular clusters (GCs) in galaxies. MUESLI follows the orbits of individual clusters and applies internal and external dissolution processes to them. Orbit integration is based on the self-consistent field method in combination with a time-transformed leapfrog scheme, allowing us to handle velocity-dependent forces like triaxial dynamical friction. In a first application, the erosion of globular cluster systems (GCSs) in elliptical galaxies is investigated. Observations show that massive ellipticals have rich, radially extended GCSs, while some compact dwarf ellipticals contain no GCs at all. For several representative examples, spanning the full mass scale of observed elliptical galaxies, we quantify the influence of radial anisotropy, galactic density profiles, SMBHs, and dynamical friction on the GC erosion rate. We find that GC number density profiles are centrally flattened in less than a Hubble time, naturally explaining observed cored GC distributions. The erosion rate depends primarily on a galaxy's mass, half-mass radius and radial anisotropy. The fraction of eroded GCs is nearly 100% in compact, M 32 like galaxies and lowest in extended and massive galaxies. Finally, we uncover the existence of a violent tidal disruption dominated phase which is important for the rapid build-up of halo stars.

We introduce version two of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). The first version (Alexander & Gieles) assumed that cluster evolution is balanced for the majority of the life-cycle, meaning that the rate of energy generation in the core of the cluster equals the diffusion rate of energy by two-body relaxation, which makes the code suitable for modelling clusters in weak tidal fields. In this new version we extend the model to include an unbalanced phase of evolution to describe the pre-collapse evolution and the accompanying escape rate such that clusters in strong tidal fields can also be modelled. We also add a prescription for the evolution of the core radius and density and a related cluster concentration parameter. The model simultaneously solves a series of first-order ordinary differential equations for the rate of change of the core radius, half-mass radius and the number of member stars N. About two thousand integration steps in time are required to solve for the entire evolution of a star cluster and this number is approximately independent of N. We compare the model to the variation of these parameters following from a series of direct N-body calculations of single-mass clusters and find good agreement in the evolution of all parameters. Relevant time-scales, such as the total lifetimes and core collapse times, are reproduced with an accuracy of about 10% for clusters with various initial half-mass radii (relative to their Jacobi radii) and a range of different initial N up to N = 65536. We intend to extend this framework to include more realistic initial conditions, such as a stellar mass spectrum and mass loss from stars. The EMACSS code can be used in star cluster population studies and in models that consider the co-evolution of (globular) star clusters and large scale structures.

Ibata et al. (2013) recently reported the existence of a vast thin plane of dwarf galaxies (VTPD) orbiting around Andromeda. We investigate whether such a configuration can be reproduced within the standard cosmological framework and search for similar planes of co-rotating satellite galaxies around Andromeda-like host haloes in data from the Millennium II simulation combined with a semi-analytic galaxy formation model. We apply a baryonic mass cut of $2.8\times 10^4 \text{M}_{\text{sun}}$ for the satellite haloes and restrict the data to a PAndAS like field. If we include the so-called orphan galaxies in our analysis, we find that planes with a rms lower than the VTPD are common in Millennium II. This is partially due to the strongly radially concentrated distribution of orphan galaxies. Excluding part of the orphan galaxies brings the radial distributions of Millennium II satellites into better agreement with the satellite distribution of Andromeda while still producing a significant fraction of planes with a lower rms than the VTPD. We also find haloes in Millennium II with an equal or higher number of co-rotating satellites than the VTPD. This demonstrates that the VTPD is not in conflict with the standard cosmological framework, although a definite answer of this question might require higher resolution cosmological simulations that do not have to consider orphan galaxies. Our results finally show that satellite planes in Millennium II are not stable structures, hence the VTPD might only be a statistical fluctuation of an underlying more spherical galaxy distribution.

Intermediate-mass black holes (IMBHs) may provide the missing link to understanding the growth of supermassive black holes in the early Universe. Some formation scenarios predict that IMBHs could have formed by runaway collisions in globular clusters (GCs). Understanding the influence of an IMBH in the center of a GC on its environment might provide indirect detection methods. Our goal is to test the effects of different initial compositions of GCs on their evolution in a tidal field. We pin down the crucial observables that indicate the presence of an IMBH at the center of the cluster. In addition to central IMBHs, we also consider the effects of different stellar-mass black hole retention and stellar binary fractions. We performed a set of 22 N-body simulations and varied particle numbers, IMBH masses, stellar-mass black-hole retention fractions, and stellar binary fractions. These models are all run in an external tidal field to study the effect of black holes on the cluster mass loss, mass function, and life times. Finally, we compared our results with observational data. We found that a central massive black hole increases the escape rate of high-mass stars from a star cluster, implying that the relative depletion of the mass function at the low-mass end proceeds less rapidly. Furthermore, we found a similar behavior for a cluster hosting a high number of stellar-mass black holes instead of one massive central IMBH. The presence of an IMBH also weakly affects the fraction of the cluster mass that is constituted by stellar remnants, as does the presence of primordial binaries. We compared our simulations with observational data from the literature and found good agreement between our models and observed mass functions and structural parameters of GCs. We exploited this agreement to identify GCs that could potentially host IMBHs.

CONTEXT: The dynamical mass-to-light (M/L) ratios of massive ultra-compact dwarf galaxies (UCDs) are about 50% higher than predicted by stellar population models. AIMS: Here we investigate the possibility that these elevated M/L ratios are caused by a central black hole (BH), heating up the internal motion of stars. We focus on a sample of ~50 extragalactic UCDs for which velocity dispersions and structural parameters have been measured. METHODS: Using up-to-date distance moduli and a consistent treatment of aperture and seeing effects, we calculate the ratio Psi=(M/L)_{dyn}/(M/L)_{pop} between the dynamical and the stellar population M/L of UCDs. For all UCDs with Psi>1 we estimate the mass of a hypothetical central BH needed to reproduce the observed integrated velocity dispersion. RESULTS: Massive UCDs (M>10^7 M_*) have an average Psi = 1.7 +-0.2, implying notable amounts of dark mass in them. We find that, on average, central BH masses of 10-15% of the UCD mass can explain these elevated dynamical M/L ratios. The implied BH masses in UCDs range from several 10^5 M_* to several 10^7 M_*. In the M_BH-Luminosity plane, UCDs are offset by about two orders of magnitude in luminosity from the relation derived for galaxies. Our findings can be interpreted such that massive UCDs originate from progenitor galaxies with masses around 10^9 M_*, and that those progenitors have SMBH occupation fractions of 60-100%. The suggested UCD progenitor masses agree with predictions from the tidal stripping scenario. Lower-mass UCDs (M<10^7 M_*) exhibit a bimodal distribution in Psi, suggestive of a coexistence of massive globular clusters and tidally stripped galaxies in this mass regime. CONCLUSIONS: Central BHs as relict tracers of tidally stripped progenitor galaxies are a plausible explanation for the elevated dynamical M/L ratios of UCDs.

We have determined possible cluster members of the nearby open cluster Praesepe (M44) based on J and K photometry and proper motions from the PPMXL catalogue and z photometry from the Sloan Digital Sky Survey (SDSS). In total we identified 893 possible cluster members down to a magnitude of J = 15.5 mag, corresponding to a mass of about 0.15 Msun for an assumed cluster distance modulus of (m-M)_0 = 6.30 mag (d ~ 182 pc), within a radius of 3.5{\deg} around the cluster centre. We derive a new cluster centre for Praesepe ({\alpha}_centre = 8h 39m 37s, {\delta}_centre = 19{\deg} 35' 02"). We also derive a total cluster mass of about 630 Msun and a 2D half-number and half-mass radius of 4.25 pc and 3.90 pc respectively. The global mass function (MF) of the cluster members shows evidence for a turnover around m = 0.65 Msun. While more massive stars can be fit by a power-law {\xi}(m) ~ m^(-{\alpha}) with slope {\alpha} = 2.88 +/- 0.22, stars less massive than m = 0.65 Msun are best fitted with {\alpha} = 0.85 +/- 0.10. In agreement with its large dynamical age, we find that Praesepe is strongly mass segregated and that the mass function slope for high mass stars steepens from a value of {\alpha} = 2.32 +/- 0.24 inside the half-mass radius to {\alpha} = 4.90 +/- 0.51 outside the half-mass radius. We finally identify a significant population of binaries and triples in the colour-magnitude diagram of Praesepe. Assuming non-random pairing of the binary components, a binary fraction of about 35% for primaries in the mass range 0.6 < m/Msun < 2.20 is required to explain the observed number of binaries in the colour-magnitude diagram (CMD).

Ultra Compact Dwarf Galaxies (UCDs) and dwarf galaxy nuclei have many common properties, such as internal velocity dispersions and colour-magnitude trends, suggesting tidally stripped dwarf galaxies as a possible UCD origin. However, UCDs typically have sizes more than twice as large as nuclei at the same luminosity. We use a GPU-enabled version of the particle-mesh code \textsc{superbox} to study the possibility of turning nucleated dwarf galaxies into UCDs by tidally stripping them in a Virgo-like galaxy cluster. We find that motion in spherical potentials, where close passages happen many times, leads to the formation of compact ($r_h \lesssim 20$ pc) star clusters/UCDs. In contrast, orbital motion where close passages happen only once or twice leads to the formation of extended objects which are large enough to account for the full range of observed UCD sizes. For such motion, we find that dwarf galaxies need close pericentre passages with distances less than 10 kpc to undergo strong enough stripping so that UCD formation is possible. As tidal stripping produces objects with similar properties to UCDs, and our estimates suggest dwarf galaxies have been destroyed in sufficient numbers to explain the observed number of UCDs in M87, we consider tidal stripping to be a likely origin of UCDs. However, comparison with cosmological simulations is needed to determine if the number and spatial distribution of UCDs formed by tidal stripping matches the observations of UCDs in galaxy clusters.

The evolution of the global stellar mass function (MF) of star clusters is studied based on a large set of N-body simulations of clusters with a range of initial masses, initial concentrations, in circular or elliptical orbits in different tidal environments. Models with and without initial mass segregation are included. The depletion of low mass stars in initially Roche-volume (tidal) filling clusters starts typically on a time scale of the order of the core collapse time. In clusters that are initially underfilling their Roche-volume it takes longer because the clusters have to expand to their tidal radii before dynamical mass loss becomes important. We introduce the concept of the differential mass function (DMF), which describes the changes with respect to the initial mass function (IMF). We show that the evolution of the DMF can be described by a set of very simple analytic expressions that are valid for a wide range of initial cluster parameters and for different IMFs. The agreement between this description and the models is very good, except for initially Roche-volume underfilling clusters that are severely mass segregated.

For galaxies hosting supermassive black holes (SMBHs), it has been observed that the mass of the central black hole (M_BH) tightly correlates with the effective or central velocity dispersion (sigma) of the host galaxy. The origin of this M_BH - sigma scaling relation is assumed to lie in the merging history of the galaxies but many open questions about its origin and the behavior in different mass ranges still need to be addressed. The goal of this work is to study the black-hole scaling relations for low black-hole masses, where the regime of intermediate-mass black holes (IMBHs) in globular clusters (GCs) is entered. We collect all existing reports of dynamical black-hole measurements in globular clusters, providing black-hole masses or upper limits for 14 candidates. We plot the black-hole masses versus different cluster parameters including total mass, velocity dispersion, concentration and half-mass radius. We search for trends and test the correlations in order to quantify their significance using a set of different statistical approaches. For correlations showing a large significance we perform a linear fit, accounting for uncertainties and upper limits. We find a clear correlation between the mass of the IMBH and the velocity dispersion of the globular cluster. As expected, the total mass of the globular cluster then also correlates with the mass of the IMBH. While the slope of the M_BH - sigma correlation differs strongly from the one observed for SMBHs, the other scaling relations M_BH - M_TOT, and M_BH - L are similar to the correlations in galaxies. Significant correlations of black-hole mass with other cluster properties were not found in the present sample.

Intermediate-mass black holes (IMBHs, 10^2-10^5 M_sun) fill the gap between stellar-mass black holes and supermassive black holes (SMBHs). Simulations have shown that IMBHs may form in dense star clusters, and therefore may still be present in these smaller stellar systems. We investigate the Galactic globular cluster NGC 5286 for indications of a central IMBH using spectroscopic data from VLT/FLAMES, velocity measurements from the Rutgers Fabry Perot at CTIO, and photometric data from HST. We run analytic spherical and axisymmetric Jeans models with different central black-hole masses, anisotropy, mass-to-light ratio, and inclination. Further, we compare the data to a grid of N-body simulations without tidal field. Additionally, we use one N-body simulation to check the results of the spherical Jeans models for the total cluster mass. Both the Jeans models and the N-body simulations favor the presence of a central black hole in NGC 5286 and our detection is at the 1- to 1.5-sigma level. From the spherical Jeans models we obtain a best fit with black-hole mass M_BH=(1.5+-1.0)x10^3 M_sun. The error is the 68% confidence limit from Monte Carlo simulations. Axisymmetric models give a consistent result. The best fitting N-body model is found with a black hole of 0.9% of the total cluster mass (4.38+-0.18)x10^5 M_sun, which results in an IMBH mass of M_BH=(3.9+-2.0)x10^3 M_sun. Jeans models give lower values for the total cluster mass. Our test of the Jeans models with N-body simulation data shows that this discrepancy has two reasons: The influence of a radially varying M/L profile, and underestimation of the velocity dispersion as the measurements are limited to bright stars. We conclude that detection of IMBHs in Galactic globular clusters remains a challenging task unless their mass fractions are above those found for SMBHs in nearby galaxies. [abridged]

For a sample of nine Galactic globular clusters we measured the inner kinematic profiles with integral-field spectroscopy that we combined with existing outer kinematic measurements and HST luminosity profiles. With this information we are able to detect the crucial rise in the velocity-dispersion profile which indicates the presence of a central black hole. In addition, N-body simulations compared to our data will give us a deeper insight in the properties of clusters with black holes and stronger selection criteria for further studies. For the first time, we obtain a homogeneous sample of globular cluster integral- field spectroscopy which allows a direct comparison between clusters with and without an intermediate-mass black hole.

The Large Magellanic Cloud is one of the nearest galaxies to us and is one of only few galaxies where the star formation history can be determined from studying resolved stellar populations. We have compiled a new catalogue of ages, luminosities and masses of LMC star clusters and used it to determine the age distribution and dissolution rate of LMC star clusters. We find that the frequency of massive clusters with masses M>5000 Msun is almost constant between 10 and 200 Myr, showing that the influence of residual gas expulsion is limited to the first 10 Myr of cluster evolution or clusters less massive than 5000 Msun. Comparing the cluster frequency in that interval with the absolute star formation rate, we find that about 15% of all stars in the LMC were formed in long-lived star clusters that survive for more than 10 Myr. We also find that the mass function of LMC clusters younger than 1 Gyr can be fitted by a power-law mass function with slope \alpha=-2.3, while older clusters follow a significantly shallower slope and interpret this is a sign of the ongoing dissolution of low-mass clusters. Our data shows that for ages older than 200 Myr, about 90% of all clusters are lost per dex of lifetime. The implied cluster dissolution rate is significantly faster than that based on analytic estimates and N-body simulations. Our cluster age data finally shows evidence for a burst in cluster formation about 1 Gyr ago, but little evidence for bursts at other ages.

The formation of supermassive black holes at high redshift still remains a puzzle to astronomers. Their growth becomes reasonable only when starting from a massive seed black hole with mass of the order of 10^2 - 10^5 M_SUN. Intermediate-mass black holes (IMBHs) are therefore an important field of research. Especially the possibility of finding them in the centers of globular clusters has recently drawn attention. The search for IMBHs in the centers of globular clusters could therefore shed light on the process of black-hole formation and cluster evolution. We are investigating six galactic globular clusters for the presence of an IMBH at their centers. Based on their kinematic and photometric properties, we selected the globular clusters NGC 1851, NGC 1904 (M79), NGC 5694, NGC 5824, NGC 6093 (M80) and NGC 6266 (M62). We use integral field spectroscopy in order to obtain the central velocity-dispersion profile of each cluster. We compute the cluster photometric center and the surface brightness profile using HST data. After combining these datasets we compare them to analytic Jeans models. We use varying M/L_V profiles for clusters with enough data points in order to reproduce their kinematic profiles in an optimal way. Finally, we vary the mass of the central black hole and test whether the cluster is better fitted with or without an IMBH. We present the statistical significance, including upper limits, of the black-hole mass for each cluster. NGC 1904 and NGC 6266 provide the highest significance for a black hole. Jeans models in combination with a M/L_V profile obtained from N-body simulations (in the case of NGC 6266) predict a central black hole of M_BH = (3 +- 1) x 10^3 M_SUN for NGC 1904 and M_BH = (2 +- 1) x 10^3 M_SUN for NGC 6266. Furthermore, we discuss the possible influence of dark remnants and mass segregation at the center of the cluster on the detection of an IMBH.

We investigate the dynamical evolution of the Orion Nebula Cluster (ONC) by means of direct N-body integrations. A large fraction of residual gas was probably expelled when the ONC formed, so we assume that the ONC was much more compact when it formed compared to its current size, in agreement with the embedded cluster radius-mass relation from Marks & Kroupa (2012). Hence, we assume that few-body relaxation played an important role during the initial phase of evolution of the ONC. In particular, three body interactions among OB stars likely led to their ejection from the cluster and, at the same time, to the formation of a massive object via runaway physical stellar collisions. The resulting depletion of the high mass end of the stellar mass function in the cluster is one of the important points where our models fit the observational data. We speculate that the runaway-mass star may have collapsed directly into a massive black hole (Mbh > 100Msun). Such a dark object could explain the large velocity dispersion of the four Trapezium stars observed in the ONC core. We further show that the putative massive black hole is likely to be a member of a binary system with appr. 70 per cent probability. In such a case, it could be detected either due to short periods of enhanced accretion of stellar winds from the secondary star during pericentre passages, or through a measurement of the motion of the secondary whose velocity would exceed 10 km/s along the whole orbit.

We highlight the impact of cluster-mass-dependent evolutionary rates upon the evolution of the cluster mass function during violent relaxation, that is, while clusters dynamically respond to the expulsion of their residual star-forming gas. Mass-dependent evolutionary rates arise when the mean volume density of cluster-forming regions is mass-dependent. In that case, even if the initial conditions are such that the cluster mass function at the end of violent relaxation has the same shape as the embedded-cluster mass function (i.e. infant weight-loss is mass-independent), the shape of the cluster mass function does change transiently {\it during} violent relaxation. In contrast, for cluster-forming regions of constant mean volume density, the cluster mass function shape is preserved all through violent relaxation since all clusters then evolve at the same mass-independent rate. On the scale of individual clusters, we model the evolution of the ratio between the dynamical mass and luminous mass of a cluster after gas expulsion. Specifically, we map the radial dependence of the time-scale for a star cluster to return to equilibrium. We stress that fields-of-view a few pc in size only, typical of compact clusters with rapid evolutionary rates, are likely to reveal cluster regions which have returned to equilibrium even if the cluster experienced a major gas expulsion episode a few Myr earlier. We provide models with the aperture and time expressed in units of the initial half-mass radius and initial crossing-time, respectively, so that our results can be applied to clusters with initial densities, sizes, and apertures different from ours.

We report the detection of five high-velocity stars in the core of the globular cluster NGC 2808. The stars lie on the the red giant branch and show total velocities between 40 and 45 km/s. For a core velocity dispersion sigma_c = 13.4 km/s, this corresponds to up to 3.4 sigma_c. These velocities are close to the estimated escape velocity (~ 50 km/s) and suggest an ejection from the core. Two of these stars have been confirmed in our recent integral field spectroscopy data and we will discuss them in more detail here. These two red giants are located at a projected distance of ~ 0.3 pc from the center. According to their positions on the color magnitude diagram, both stars are cluster members. We investigate several possible origins for the high velocities of the stars and conceivable ejection mechanisms. Since the velocities are close to the escape velocity, it is not obvious whether the stars are bound or unbound to the cluster. We therefore consider both cases in our analysis. We perform numerical simulations of three-body dynamical encounters between binaries and single stars and compare the resulting velocity distributions of escapers with the velocities of our stars. We compare the predictions for a single dynamical encounter with a compact object with those of a sequence of two-body encounters due to relaxation. If the stars are unbound, the encounter must have taken place recently, when the stars were already in the giant phase. After including binary fractions and black-hole retention fractions, projection effects, and detection probabilities from Monte-Carlo simulations, we estimate the expected numbers of detections for all the different scenarios. Based on these numbers, we conclude that the most likely scenario is that the stars are bound and were accelerated by a single encounter between a binary of main-sequence stars and a ~ 10 M_sun black hole.

Globular clusters are an excellent laboratory for stellar population and dynamical research. Recent studies have shown that these stellar systems are not as simple as previously assumed. With multiple stellar populations as well as outer rotation and mass segregation they turn out to exhibit high complexity. This includes intermediate-mass black holes which are proposed to sit at the centers of some massive globular clusters. Today's high angular resolution ground based spectrographs allow velocity-dispersion measurements at a spatial resolution comparable to the radius of influence for plausible IMBH masses, and to detect changes in the inner velocity-dispersion profile. Together with high quality photometric data from HST, it is possible to constrain black-hole masses by their kinematic signatures. We determine the central velocity-dispersion profile of the globular cluster NGC 2808 using VLT/FLAMES spectroscopy. In combination with HST/ACS data our goal is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and constrain the cluster mass to light ratio as well as its total mass. We derive a velocity-dispersion profile from integral field spectroscopy in the center and Fabry Perot data for larger radii. High resolution HST data are used to obtain the surface brightness profile. Together, these data sets are compared to dynamical models with varying parameters such as mass to light ratio profiles and black-hole masses. Using analytical Jeans models in combination with variable M/L profiles from N-body simulations we find that the best fit model is a no black hole solution. After applying various Monte Carlo simulations to estimate the uncertainties, we derive an upper limit of the back hole mass of M_BH < 1 x 10^4 M_SUN (with 95 % confidence limits) and a global mass-to-light ratio of M/L_V = (2.1 +- 0.2) M_SUN/L_SUN.

Most recent progress in understanding the dynamical evolution of star clusters relies on direct N-body simulations. Owing to the computational demands, and the desire to model more complex and more massive star clusters, hardware calculational accelerators, such as GRAPE special-purpose hardware or, more recently, GPUs (i.e. graphics cards), are generally utilised. In addition, simulations can be accelerated by adjusting parameters determining the calculation accuracy (i.e. changing the internal simulation time step used for each star). We extend our previous thorough comparison (Anders et al. 2009) of basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY6. Here we focus on differences arising from using different hardware accelerations (including the increasingly popular graphic card accelerations/GPUs) and different calculation accuracy settings. We use the large number of star cluster models (for a fixed stellar mass function, without stellar/binary evolution, primordial binaries, external tidal fields etc) already used in the previous paper, evolve them with STARLAB/KIRA (and NBODY6, where required), analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we apply the bootstrap algorithm for functional dependencies developed in our previous study. In general we find very high comparability of the simulation results, independent of the used computer hardware (including the hardware accelerators) and the used N-body code. For the tested accuracy settings we find that for reduced accuracy (i.e. time step at least a factor 2.5 larger than the standard setting) most simulation results deviate significantly from the results using standard settings. The remaining deviations are comprehensible and explicable.

Supermassive black holes (SMBHs) are fundamental keys to understand the formation and evolution of their host galaxies. However, the formation and growth of SMBHs are not yet well understood. One of the proposed formation scenarios is the growth of SMBHs from seed intermediate-mass black holes (IMBHs, 10^2 to 10^5 M_{\odot}) formed in star clusters. In this context, and also with respect to the low mass end of the M-sigma relation for galaxies, globular clusters are in a mass range that make them ideal systems to look for IMBHs. Among Galactic star clusters, the massive cluster $\omega$ Centauri is a special target due to its central high velocity dispersion and also its multiple stellar populations. We study the central structure and dynamics of the star cluster $\omega$ Centauri to examine whether an IMBH is necessary to explain the observed velocity dispersion and surface brightness profiles. We perform direct N-body simulations to follow the dynamical evolution of $\omega$ Centauri. The simulations are compared to the most recent data-sets in order to explain the present-day conditions of the cluster and to constrain the initial conditions leading to the observed profiles. We find that starting from isotropic spherical multi-mass King models and within our canonical assumptions, a model with a central IMBH mass of 2% of the cluster stellar mass, i.e. a 5x10^4 M_{\odot} IMBH, provides a satisfactory fit to both the observed shallow cusp in surface brightness and the continuous rise towards the center of the radial velocity dispersion profile. In our isotropic spherical models, the predicted proper motion dispersion for the best-fit model is the same as the radial velocity dispersion one. (abridged)

This paper presents a new scheme to treat escaping stars in the orbit-averaged Fokker-Planck models of globular star clusters in a galactic tidal field. The existence of a large number of potential escapers, which have energies above the escape energy but are still within the tidal radius, is taken into account in the models. The models allow potential escapers to experience gravitational scatterings before they leave clusters and thus some of them may lose enough energy to be bound again. It is shown that the mass evolution of the Fokker-Planck models are in good agreement with that of N-body models including the full tidal-force field. The mass-loss time does not simply scale with the relaxation time due to the existence of potential escapers; it increases with the number of stars more slowly than the relaxation time, though it tends to be proportional to the relaxation time in the limit of a weak tidal field. The Fokker-Planck models include two parameters, the coefficient gamma in the Coulomb logarithm ln(gamma N) and the coefficient nu_e controlling the efficiency of the mass loss. The values of these parameters are determined by comparing the Fokker-Planck models with the N-body models. It is found that the parameter set (gamma, nu_e)=(0.11, 7) works well for both single-mass and multi-mass clusters, but that the parameter set (gamma, nu_e)=(0.02, 40) is another possible choice for multi-mass clusters.

Surface photometry is a necessary tool to establish the dynamical state of stars clusters. We produce realistic HST-like images from N-body models of star clusters with and without central intermediate-mass black holes (IMBHs) in order to measure their surface brightness profiles. The models contain ~600,000 individual stars, black holes of various masses between 0% to 2% of the total mass, and are evolved for a Hubble time. We measure surface brightness and star count profiles for every constructed image in order to test the effect of intermediate mass black holes on the central logarithmic slope, the core radius, and the half-light radius. We use these quantities to test diagnostic tools for the presence of central black holes using photometry. We find that the the only models that show central shallow cusps with logarithmic slopes between -0.1 and -0.4 are those containing central black holes. Thus, the central logarithmic slope seems to be a good way to choose clusters suspect of containing intermediate-mass black holes. Clusters with steep central cusps can definitely be ruled out to host an IMBH. The measured r_c/r_h ratio has similar values for clusters that have not undergone core-collapse, and those containing a central black hole. We notice that observed Galactic globular clusters have a larger span of values for central slope and r_c/r_h than our modeled clusters, and suggest possible reasons that could account for this and contribute to improve future models.

By analysing models of the young massive cluster R136 in 30 Doradus, set-up using the herewith introduced and publicly made available code McLuster, we investigate and compare different methods for detecting and quantifying mass segregation and substructure in non-seeing limited N-body data. For this purpose we generate star cluster models with different degrees of mass segregation and fractal substructure and analyse them. We quantify mass segregation by measuring, from the projected 2d model data, the mass function slope in radial annuli, by looking for colour gradients in radial colour profiles, by measuring Allison's Lambda parameter, and by determining the local stellar surface density around each star. We find that these methods for quantifying mass segregation often produce ambiguous results. Most reliable for detecting mass segregation is the mass function slope method, whereas the colour gradient method is the least practical in an R136-like configuration. The other two methods are more sensitive to low degrees of mass segregation but are computationally much more demanding. We also discuss the effect of binaries on these measures. Moreover, we quantify substructure by looking at the projected radial stellar density profile, by comparing projected azimuthal stellar density profiles, and by determining Cartwright & Whitworth's Q parameter. We find that only high degrees of substructure affect the projected radial density profile, whereas the projected azimuthal density profile is very sensitive to substructure. The Q parameter is also sensitive to substructure but its absolute value shows a dependence on the radial density gradient of the cluster and is strongly influenced by binaries. (abridged)

The disruption rate of stars by supermassive black holes (SMBHs) is calculated numerically with a modified version of Aarseth's NBODY6 code. The initial stellar distribution around the SMBH follows a S\'{e}rsic n=4 profile representing bulges and early type galaxies. In order to infer relaxation driven effects and to increase the statistical significance, a very large set of N-body integrations with different particle numbers N, ranging from 10^{3} to 0.5 \cdot 10^{6} particles, is performed. Three different black hole capture radii are taken into account, enabling us to scale these results to a broad range of astrophysical systems with relaxation times shorter than one Hubble time, i.e. for SMBHs up to M_bh \approx 10^{7} M_sun. The computed number of disrupted stars are driven by diffusion in angular momentum space into the loss cone of the black hole and the rate scales with the total number of particles as dN/dt \propto N^{b}, where b is as large as 0.83. This is significantly steeper than the expected scaling dN/dt \propto ln(N) derived from simplest energy relaxation arguments. Only a relatively modest dependence of the tidal disruption rate on the mass of the SMBH is found and we discuss our results in the context of the M_bh/sigma relation. The number of disrupted stars contribute a significant part to the mass growth of black holes in the lower mass range as long as a significant part of the stellar mass becomes swallowed by the SMBH. This also bears direct consequences for the search and existence of IMBHs in globular clusters. For SMBHs similar to the galactic center black hole SgrA*, a tidal disruption rate of 55 \pm 27 events per Myr is deduced. Finally relaxation driven stellar feeding can not account for the masses of massive black holes M_bh \geq 10^{7} M_sun. (abridged)

Intermediate-mass black holes (IMBHs) are of interest in a wide range of astrophysical fields. In particular, the possibility of finding them at the centers of globular clusters has recently drawn attention. IMBHs became detectable since the quality of observational data sets, particularly those obtained with HST and with high resolution ground based spectrographs, advanced to the point where it is possible to measure velocity dispersions at a spatial resolution comparable to the size of the gravitational sphere of influence for plausible IMBH masses. We present results from ground based VLT/FLAMES spectroscopy in combination with HST data for the globular cluster NGC 6388. The aim of this work is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and to compare the results with the expected value predicted by the $M_{\bullet} - \sigma$ scaling relation. The spectroscopic data, containing integral field unit measurements, provide kinematic signatures in the center of the cluster while the photometric data give information of the stellar density. Together, these data sets are compared to dynamical models and present evidence of an additional compact dark mass at the center: a black hole. Using analytical Jeans models in combination with various Monte Carlo simulations to estimate the errors, we derive (with 68% confidence limits) a best fit black-hole mass of $ (17 \pm 9) \times 10^3 M_{\odot}$ and a global mass-to-light ratio of $M/L_V = (1.6 \pm 0.3) \ M_{\odot}/L_{\odot}$.

We present the internal kinematics of UCD3, the brightest known ultra-compact dwarf galaxy (UCD) in the Fornax cluster, making this the first UCD with spatially resolved spectroscopy. Our study is based on seeing-limited observations obtained with the ARGUS Integral Field Unit of the VLT/FLAMES spectrograph under excellent seeing conditions (0.5 - 0.67 arcsec FWHM). The velocity field of UCD3 shows the signature of weak rotation, comparable to that found in massive globular clusters. Its velocity dispersion profile is fully consistent with an isotropic velocity distribution and the assumption that mass follows the light distribution obtained from Hubble Space Telescope imaging. In particular, there is no evidence for the presence of an extended dark matter halo contributing a significant (>~33 per cent within R < 200 pc) mass fraction, nor for a central black hole more massive than ~5 per cent of the UCD's mass. While this result does not exclude a galaxian origin for UCD3, we conclude that its internal kinematics are fully consistent with it being a massive star cluster.

We determine the mean velocity dispersion of six Galactic outer halo globular clusters, AM 1, Eridanus, Pal 3, Pal 4, Pal 15, and Arp 2 in the weak acceleration regime to test classical vs. modified Newtonian dynamics (MOND). Owing to the non-linearity of MOND's Poisson equation, beyond tidal effects, the internal dynamics of clusters is affected by the external field in which they are immersed. For the studied clusters, particle accelerations are much lower than the critical acceleration a_0 of MOND, but the motion of stars is neither dominated by internal accelerations (a_i >> a_e) nor external accelerations (a_e >> a_i). We use the N-body code N-MODY in our analysis, which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo, and Nipoti (2006) to derive the line-of-sight velocity dispersion by adding the external field effect. We show that Newtonian dynamics predicts a low-velocity dispersion for each cluster, while in modified Newtonian dynamics the velocity dispersion is much higher. We calculate the minimum number of measured stars necessary to distinguish between Newtonian gravity and MOND with the Kolmogorov-Smirnov test. We also show that for most clusters it is necessary to measure the velocities of between 30 to 80 stars to distinguish between both cases. Therefore the observational measurement of the line-of-sight velocity dispersion of these clusters will provide a test for MOND.

We use direct N-body simulations of gas embedded star clusters to study the importance of stellar collisions for the formation and mass accretion history of high-mass stars. Our clusters start in virial equilibrium as a mix of gas and proto-stars. Proto-stars then accrete matter using different mass accretion rates and the amount of gas is reduced in the same way as the mass of stars increases. During the simulations we check for stellar collisions and we investigate the role of these collisions for the build-up of high-mass stars and the formation of runaway stars. We find that a significant number of collisions only occur in clusters with initial half-mass radii r_h < 0.1 pc. After emerging from their parental gas clouds, such clusters end up too compact compared to observed young, massive open clusters. In addition, collisions lead mainly to the formation of a single runaway star instead of the formation of many high mass stars with a broad mass spectrum. We therefore conclude that massive stars form mainly by gas accretion, with stellar collisions only playing a minor role if any at all. Collisions of stars in the pre-main sequence phase might however contribute to the formation of the most massive stars in the densest star clusters and possibly to the formation of intermediate-mass black holes with masses up to a few 100 Msun.

We present the first ever direct $N$-body computations of an old Milky Way globular cluster over its entire life time on a star-by-star basis. Using recent GPU hardware at Bonn University, we have performed a comprehensive set of $N$-body calculations to model the distant outer halo globular cluster Palomar 14 (Pal 14). By varying the initial conditions we aim at finding an initial $N$-body model which reproduces the observational data best in terms of its basic parameters, i.e. half-light radius, mass and velocity dispersion. We furthermore focus on reproducing the stellar mass function slope of Pal 14 which was found to be significantly shallower than in most globular clusters. While some of our models can reproduce Pal 14's basic parameters reasonably well, we find that dynamical mass segregation alone cannot explain the mass function slope of Pal 14 when starting from the canonical Kroupa initial mass function (IMF). In order to seek for an explanation for this discrepancy, we compute additional initial models with varying degrees of primordial mass segregation as well as with a flattened IMF. The necessary degree of primordial mass segregation turns out to be very high. This modelling has shown that the initial conditions of Pal 14 after gas expulsion must have been a half-mass radius of about 20 pc, a mass of about 50000 M$_{\odot}$, and possibly some mass segregation or an already established non-canonical IMF depleted in low-mass stars. Such conditions might be obtained by a violent early gas-expulsion phase from an embedded cluster born with mass segregation. Only at large Galactocentric radii are clusters likely to survive as bound entities the destructive gas-expulsion process we seem to have uncovered for Pal 14. In addition we compute a model with a 5% primordial binary fraction to test if such a population has an effect on the cluster's evolution.

The Galactic globular cluster omega Centauri is a prime candidate for hosting an intermediate mass black hole. Recent measurements lead to contradictory conclusions on this issue. We use VLT-FLAMES to obtain new integrated spectra for the central region of omega Centauri. We combine these data with existing measurements of the radial velocity dispersion profile taking into account a new derived center from kinematics and two different centers from the literature. The data support previous measurements performed for a smaller field of view and show a discrepancy with the results from a large proper motion data set. We see a rise in the radial velocity dispersion in the central region to 22.8+-1.2 km/s, which provides a strong sign for a central black hole. Isotropic dynamical models for omega Centauri imply black hole masses ranging from 3.0 to 5.2x10^4 solar masses depending on the center. The best-fitted mass is 4.7+-1.0x10^4 solar masses.

Most globular clusters have half-mass radii of a few pc with no apparent correlation with their masses. This is different from elliptical galaxies, for which the Faber-Jackson relation suggests a strong positive correlation between mass and radius. Objects that are somewhat in between globular clusters and low-mass galaxies, such as ultra-compact dwarf galaxies, have a mass-radius relation consistent with the extension of the relation for bright ellipticals. Here we show that at an age of 10 Gyr a break in the mass-radius relation at ~10^6 Msun is established because objects below this mass, i.e. globular clusters, have undergone expansion driven by stellar evolution and hard binaries. From numerical simulations we find that the combined energy production of these two effects in the core comes into balance with the flux of energy that is conducted across the half-mass radius by relaxation. An important property of this `balanced' evolution is that the cluster half-mass radius is independent of its initial value and is a function of the number of bound stars and the age only. It is therefore not possible to infer the initial mass-radius relation of globular clusters and we can only conclude that the present day properties are consistent with the hypothesis that all hot stellar systems formed with the same mass-radius relation and that globular clusters have moved away from this relation because of a Hubble time of stellar and dynamical evolution.

We describe the interplay between stellar evolution and dynamical mass loss of evolving star clusters, based on the principles of stellar evolution and cluster dynamics and on a grid of N-body simulations of cluster models. The cluster models have different initial masses, different orbits, including elliptical ones, and different initial density profiles. We use two sets of cluster models: initially Roche-lobe filling and Roche-lobe underfilling. We identify four distinct mass loss effects: (1) mass loss by stellar evolution, (2) loss of stars induced by stellar evolution and (3) relaxation-driven mass loss before and (4) after core collapse. Both the evolution-induced loss of stars and the relaxation-driven mass loss need time to build up. This is described by a delay-function of a few crossing times for Roche-lobe filling clusters and a few half mass relaxation times for Roche-lobe underfilling clusters. The relaxation-driven mass loss can be described by a simple power law dependence of the mass dM/dt =-M^{1-gamma}/t0, (with M in Msun) where t0 depends on the orbit and environment of the cluster. Gamma is 0.65 for clusters with a King-parameter W0=5 and 0.80 for more concentrated clusters with W0=7. For initially Roche-lobe underfilling clusters the dissolution is described by the same gamma=0.80. The values of the constant t0 are described by simple formulae that depend on the orbit of the cluster. The mass loss rate increases by about a factor two at core collapse and the mass dependence of the relaxation-driven mass loss changes to gamma=0.70 after core collapse. We also present a simple recipe for predicting the mass evolution of individual star clusters with various metallicities and in different environments, with an accuracy of a few percent in most cases. This can be used to predict the mass evolution of cluster systems.

Based on our recent work on tidal tails of star clusters (Kuepper et al. 2009) we investigate star clusters of a few 10^4 Msun by means of velocity dispersion profiles and surface density profiles. We use a comprehensive set of $N$-body computations of star clusters on various orbits within a realistic tidal field to study the evolution of these profiles with time, and ongoing cluster dissolution From the velocity dispersion profiles we find that the population of potential escapers, i.e. energetically unbound stars inside the Jacobi radius, dominates clusters at radii above about 50% of the Jacobi radius. Beyond 70% of the Jacobi radius nearly all stars are energetically unbound. The velocity dispersion therefore significantly deviates from the predictions of simple equilibrium models in this regime. We furthermore argue that for this reason this part of a cluster cannot be used to detect a dark matter halo or deviations from Newtonian gravity. By fitting templates to the about 10^4 computed surface density profiles we estimate the accuracy which can be achieved in reconstructing the Jacobi radius of a cluster in this way. We find that the template of King (1962) works well for extended clusters on nearly circular orbits, but shows significant flaws in the case of eccentric cluster orbits. This we fix by extending this template with 3 more free parameters. Our template can reconstruct the tidal radius over all fitted ranges with an accuracy of about 10%, and is especially useful in the case of cluster data with a wide radial coverage and for clusters showing significant extra-tidal stellar populations. No other template that we have tried can yield comparable results over this range of cluster conditions. All templates fail to reconstruct tidal parameters of concentrated clusters, however. (abridged)

We investigate the early evolution of two distinct populations of low-mass stars in globular clusters under the influence of primordial gas expulsion driven by supernovae to study if this process can increase the fraction of second generation stars at the level required by observations. We analyse N-body models that take into account the effect of primordial gas expulsion. We divide the stars into two populations which mimic the chemical and dynamical properties of stars in globular clusters so that second generation stars start with a more centrally concentrated distribution. The main effect of gas expulsion is to eject preferentially first generation stars while second generation stars remain bound to the cluster. In the most favourable cases second generation stars can account for 60% of the bound stars we see today. We also find that at the end of the gas expulsion phase, the radial distribution of the two populations is still different, so that long-term evolution will further increase the fraction of second generation stars. The large fraction of chemically anomalous stars is readily explainable as a second generation of stars formed out of the slow winds of rapidly rotating massive stars if globular clusters suffer explosive residual gas expulsion for a star formation efficiency of about 0.33.

It has been shown recently that the dynamical V-band mass-to-light ratios of compact stellar systems with masses from 10^6 to 10^8 Solar masses are not consistent with the predictions from simple stellar population (SSP) models. Top-heavy stellar initial mass functions (IMFs) in these so-called ultra compact dwarf galaxies (UCDs) offer an attractive explanation for this finding, the stellar remnants and retained stellar envelopes providing the unseen mass. We therefore construct a model which quantifies by how much the IMFs of UCDs would have to deviate in the intermediate-mass and high-mass range from the canonical IMF in order to account for the enhanced M/L_V ratio of the UCDs. The deduced high-mass IMF in the UCDs depends on the age of the UCDs and the number of faint products of stellar evolution retained by them. Assuming that the IMF in the UCDs is a three-part power-law equal to the canonical IMF in the low-mass range and taking 20% as a plausible choice for the fraction of the remnants of high-mass stars retained by UCDs, the model suggests the exponent of the high-mass IMF to be approximately 1.6 if the UCDs are 13 Gyr old (i.e. almost as old as the Universe) or approximately 1.0 if the UCDs are 7 Gyr old, in contrast to 2.3 for the Salpeter-Massey IMF. If the IMF was as top-heavy as suggested here, the stability of the UCDs might have been threatened by heavy mass loss induced by the radiation and evolution of massive stars. The central densities of UCDs must have been in the range 10^6 to 10^7 Solar masses per cubic parsec when they formed with star formation rates of 10 to 100 Solar masses per year.

We investigate the color-magnitude relation for globular clusters (GCs) -- the so-called "blue tilt" -- detected in the ACS Fornax Cluster Survey and using the combined sample of GCs from the ACS Fornax and Virgo Cluster Surveys. We find a tilt of gamma_z=d(g-z)/dz=-0.0257 +- 0.0050 for the full GC sample of the Fornax Cluster Survey (~5800 GCs). This is slightly shallower than the value gamma_z=-0.0459 +- 0.0048 found for the Virgo Cluster Survey GC sample (~11100 GCs). The slope for the merged Fornax and Virgo datasets (~16900 GCs) is gamma_z=-0.0293 +- 0.0085, corresponding to a mass-metallicity relation of Z ~ M^0.43. We find that the blue tilt sets in at GC masses in excess of M ~ 2*10^5 M_sun. The tilt is stronger for GCs belonging to high-mass galaxies (M_* > 5 * 10^10 M_sun) than for those in low-mass galaxies (M_* < 5 * 10^10 M_sun). It is also more pronounced for GCs with smaller galactocentric distances. Our findings suggest a range of mass-metallicity relations Z_GC ~ M_GC^(0.3-0.7) which vary as a function of host galaxy mass/luminosity. We compare our observations to a recent model of star cluster self-enrichment with generally favorable results. We suggest that, within the context of this model, the proto-cluster clouds out of which the GCs formed may have had density profiles slightly steeper than isothermal and/or star formation efficiencies somewhat below 0.3. We caution, however, that the significantly different appearance of the CMDs defined by the GC systems associated with galaxies of similar mass and morphological type pose a challenge to any single mechanism that seeks to explain the blue tilt. We therefore suggest that the merger/accretion histories of individual galaxies have played a non-negligible role determining the distribution of GCs in the CMDs of individual GC systems.

We analyse the potential migration of massive planets forming far away from an inner planetary system. For this, we follow the dynamical evolution of the orbital elements of a massive planet undergoing a dissipative process with a gas disc centred around the central sun. We use a new method for post-Newtonian, high-precision integration of planetary systems containing a central sun by splitting the forces on a particle between a dominant central force and additional perturbations. In this treatment, which allows us to integrate with a very high-accuracy close encounters, all gravitational forces are integrated directly, without resorting to any simplifying approach. After traversing the disc a number of times, the planet is finally trapped into the disc with a non-negligible eccentricity

(abridged) Here we discuss the question whether the extreme circumstances in the centre of the Milky Way may be the reason for a significant variation of the IMF. By means of stellar evolution models using different codes we show that the observed luminosity in the central parsec is too high to be explained by a long-standing top-heavy IMF, considering the limited amount of mass inferred from stellar kinematics in this region. In contrast, continuous star formation over the Galaxy's lifetime following a canonical IMF results in a mass-to-light ratio and a total mass of stellar black holes (SBHs) consistent with the observations. Furthermore, these SBHs migrate towards the centre due to dynamical friction, turning the cusp of visible stars into a core as observed in the Galactic Centre. For the first time here we explain the luminosity and dynamical mass of the central cluster and both the presence and extent of the observed core, since the number of SBHs expected from a canonical IMF is just enough to make up for the missing luminous mass. We conclude that the Galactic Centre is consistent with the canonical IMF and do not suggest a systematic variation as a result of the region's properties such as high density, metallicity, strong tidal field etc.

We study the dynamics of stellar-mass black holes (BH) in star clusters with particular attention to the formation of BH-BH binaries, which are interesting as sources of gravitational waves (GW). We examine the properties of these BH-BH binaries through direct N-body simulations of star clusters using the GPU-enabled NBODY6 code. We perform simulations of N <= 10^5 Plummer clusters of low-mass stars with an initial population of BHs. Additionally, we do several calculations of star clusters confined within a reflective boundary mimicking only the core of a massive cluster. We find that stellar-mass BHs with masses ~ 10 solar mass segregate rapidly into the cluster core and form a sub-cluster of BHs within typically 0.2 - 0.5 pc radius, which is dense enough to form BH-BH binaries through 3-body encounters. While most BH binaries are ejected from the cluster by recoils received during super-elastic encounters with the single BHs, few of them harden sufficiently so that they can merge via GW emission within the cluster. We find that for clusters with $N \ga 5\times 10^4$, typically 1 - 2 BH-BH mergers occur within them during the first ~ 4 Gyr of evolution. Also for each of these clusters, there are a few escaping BH binaries that can merge within a Hubble time, most of the merger times being within a few Gyr. These results indicate that intermediate-age massive clusters constitute the most important class of candidates for producing dynamical BH-BH mergers. Old globular clusters cannot contribute significantly to the present-day BH-BH merger rate since most of the mergers from them would have occurred earlier. In contrast, young massive clusters are too young to produce significant number of BH-BH mergers. Our results imply significant BH-BH merger detection rates for the proposed "Advanced LIGO" GW detector. (Abridged)

Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) 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 disk 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 sub-structures in the haloes of galaxies.

We investigate the ratio between the half-mass radii r_h of Galactic globular clusters and their Jacobi radii r_J given by the potential of the Milky Way and show that clusters with galactocentric distances R_{GC}>8 kpc fall into two distinct groups: one group of compact, tidally-underfilling clusters with r_h/r_J<0.05 and another group of tidally filling clusters which have 0.1 < r_h/r_J<0.3. We find no correlation between the membership of a particular cluster to one of these groups and its membership in the old or younger halo population. Based on the relaxation times and orbits of the clusters, we argue that compact clusters and most clusters in the inner Milky Way were born compact with half-mass radii r_h < 1 pc. Some of the tidally-filling clusters might have formed compact as well, but the majority likely formed with large half-mass radii. Galactic globular clusters therefore show a similar dichotomy as was recently found for globular clusters in dwarf galaxies and for young star clusters in the Milky Way. It seems likely that some of the tidally-filling clusters are evolving along the main sequence line of clusters recently discovered by Kuepper et al. (2008) and are in the process of dissolution.

Observations of the Galactic Centre show evidence of one or two disc-like structures of very young stars orbiting the central super-massive black hole within a distance of a few 0.1 pc. A number of analyses have been carried out to investigate the dynamical behaviour and consequences of these discs, including disc thickness and eccentricity growth as well as mutual interaction and warping. However, most of these studies have neglected the influence of the stellar cusp surrounding the black hole, which is believed to be 1-2 orders of magnitude more massive than the disc(s). By means of N-body integrations using our bhint code, we study the impact of stellar cusps of different compositions. We find that although the presence of a cusp does have an important effect on the evolution of an otherwise isolated flat disc, its influence on the evolution of disc thickness and warping is rather mild in a two-disc configuration. However, we show that the creation of highly eccentric orbits strongly depends on the graininess of the cusp (i.e. the mean and maximum stellar masses): While Chang (2009) recently found that full cycles of Kozai resonance are prevented by the presence of an analytic cusp, we show that relaxation processes play an important role in such highly dense regions and support short-term resonances. We thus find that young disc stars on initially circular orbits can achieve high eccentricities by resonant effects also in the presence of a cusp of stellar remnants, yielding a mechanism to create S-stars and hyper-velocity stars. Furthermore, we discuss the underlying initial mass function (IMF) of the young stellar discs and find no definite evidence for a non-canonical IMF.

Evolutionary synthesis models are the prime method to construct models of stellar populations, and to derive physical parameters from observations. One of the assumptions for such models so far has been the time-independence of the stellar mass function. However, dynamical simulations of star clusters in tidal fields have shown the mass function to change due to the preferential removal of low-mass stars from clusters. Here we combine the results from dynamical simulations of star clusters in tidal fields with our evolutionary synthesis code GALEV to extend the models by a new dimension: the total cluster disruption time. We reanalyse the mass function evolution found in N-body simulations of star clusters in tidal fields, parametrise it as a function of age and total cluster disruption time and use this parametrisation to compute GALEV models as a function of age, metallicity and the total cluster disruption time. We study the impact of cluster dissolution on the colour (generally, they become redder) and magnitude (they become fainter) evolution of star clusters, their mass-to-light ratios (off by a factor of ~2 -- 4 from standard predictions), and quantify the effect on the cluster age determination from integrated photometry (in most cases, clusters appear to be older than they are, between 20 and 200%). By comparing our model results with observed M/L ratios for old compact objects in the mass range 10^4.5 -- 10^8 Msun, we find a strong discrepancy for objects more massive than 10^7 Msun (higher M/L). This could be either caused by differences in the underlying stellar mass function or be an indication for the presence of dark matter in these objects. Less massive objects are well represented by the models. The models for a range of total cluster disruption times are available online. (shortened)

Precise radial velocity measurements from HIRES on the Keck I telescope are presented for 40 stars in the outer halo globular cluster NGC 2419. These data are used to probe the cluster's stellar mass function and search for the presence of dark matter in this cluster. NGC 2419 is one of the best Galactic globular clusters for such a study due to its long relaxation time (T_{r0} ~ 10^{10} yr) and large Galactocentric distance (R_{GC} ~ 90 kpc) -- properties that make significant evolutionary changes in the low-mass end of the cluster mass function unlikely. We find a mean cluster velocity of <v_r>=-20.3 +- 0.7 km/sec and an internal velocity dispersion of \sigma = 4.14 +- 0.48 km/sec, leading to a total mass of (9.0 +- 2.2) * 10^5 Msun and a global mass-to-light ratio of M/L_V = 2.05 +- 0.50 in solar units. This mass-to-light ratio is in good agreement with what one would expect for a pure stellar system following a standard mass function at the metallicity of NGC 2419. In addition, the mass-to-light ratio does not appear to rise towards the outer parts of the cluster. Our measurements therefore rule out the presence of a dark matter halo with mass larger than ~10^7 Msun inside the central 500 pc, which is lower than what is found for the central dark matter densities of dSph galaxies. We also discuss the relevance of our measurements for alternative gravitational theories such as MOND, and for possible formation scenarios of ultra-compact dwarf galaxies.

We use the distant outer halo globular cluster Palomar 14 as a test case for classical vs. modified Newtonian dynamics (MOND). Previous theoretical calculations have shown that the line-of-sight velocity dispersion predicted by these theories can differ by up to a factor of three for such sparse, remote clusters like Pal 14. We determine the line-of-sight velocity dispersion of Palomar 14 by measuring radial velocities of 17 red giant cluster members obtained using the Very Large Telescope (VLT) and Keck telescope. The systemic velocity of Palomar 14 is (72.28+-0.12)km/s. The derived velocity dispersion of (0.38+-0.12)km/s of the 16 definite member stars is in agreement with the theoretical prediction for the classical Newtonian case according to Baumgardt et al. (2005). In order to exclude the possibility that a peculiar mass function might have influenced our measurements, we derived the cluster's main sequence mass function down to 0.53Msolar using archival images obtained with the Hubble Space Telescope. We found a mass function slope of 1.27+-0.44, which is, compared to the canonical mass function, a significantly shallower slope. The derived lower limit on the cluster's mass is higher than the theoretically predicted mass in case of MOND. Our data are consistent with a central density of 0.1 Msolar pc^-3. We need no dark matter in Palomar 14. If the cluster is on a circular orbit, our spectroscopic and photometric results argue against MOND, unless this cluster experienced significant mass loss.

(Abridged) Using luminosities and structural parameters of globular clusters (GCs) in the nuclear regions (nGCs) of low-mass dwarf galaxies from HST/ACS imaging we derive the present-day escape velocities (v_esc) of stellar ejecta to reach the cluster tidal radius and compare them with those of Galactic GCs with extended (hot) horizontal branches (EHBs-GCs). For EHB-GCs, we find a correlation between the present-day v_esc and their metallicity as well as (V-I)_0 colour. The similar v_esc, (V-I)_0 distribution of nGCs and EHB-GCs implies that nGCs could also have complex stellar populations. The v_esc-[Fe/H] relation could reflect the known relation of increasing stellar wind velocity with metallicity, which in turn could explain why more metal-poor clusters typically show more peculiarities in their stellar population than more metal-rich clusters of the same mass do. Thus the cluster v_esc can be used as parameter to describe the degree of self-enrichment. The nGCs populate the same Mv vs. rh region as EHB-GCs, although they do not reach the sizes of the largest EHB-GCs like wCen and NGC 2419. We argue that during accretion the rh of an nGC could increase due to significant mass loss in the cluster vicinity and the resulting drop in the external potential in the core once the dwarf galaxy dissolves. Our results support the scenario in which Galactic EHB-GCs have originated in the centres of pre-Galactic building blocks or dwarf galaxies that were later accreted by the Milky Way.

N-body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N-body simulations. They provide insight in the structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N-body codes STARLAB/KIRA and NBODY4 are widely used for a range of applications, there is no thorough comparison study yet. Here we thoroughly compare basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY4. We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields etc), evolve them in parallel with STARLAB/KIRA and NBODY4, analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison we develop a bootstrap algorithm for functional dependencies. We find an overall excellent agreement between the codes, both for the clusters' structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies. Our results reassure the comparability and the possibility to combine results from these two major N-body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. (abridged)

We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in MOND, using the N-body code N-MODY, which is a particle-mesh based code with a numerical MOND potential solver developed by Ciotti, Londrillo and Nipoti (2006). We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of $10^4 M_\odot$ to $10^9 M_\odot$ and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime ($a_i, a_e \ll a_0$), where the motion of stars is either dominated by internal accelerations ($a_i \gg a_e$) or constant external accelerations ($a_e \gg a_i$). In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime ($a_i \sim a_e \sim a_0$). This allows for a much improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.

Observations of the Galactic Centre show evidence of disc-like structures of very young stars orbiting the central super-massive black hole within a distance of a few 0.1 pc. While it is widely accepted that about half of the stars form a relatively flat disc rotating clockwise on the sky, there is a substantial ongoing debate on whether there is a second, counter-clockwise disc of stars. By means of N-body simulations using our bhint code, we show that two highly inclined stellar discs with the observed properties cannot be recognised as two flat circular discs after 5 Myr of mutual interaction. Instead, our calculations predict a significant warping of the two discs, which we show to be apparent among the structures observed in the Galactic Centre. While the high eccentricities of the observed counter-clockwise orbits suggest an eccentric origin of this system, we show the eccentricity distribution in the inner part of the more massive clockwise disc to be perfectly consistent with an initially circular disc in which stellar eccentricities increase due to both non-resonant and resonant relaxation. We conclude that the relevant question to ask is therefore not whether there are two discs of young stars, but whether there were two such discs to begin with.

Stars in globular clusters (GCs) exhibit a peculiar chemical pattern with strong abundance variations in light elements along with a constant abundance in heavy elements. These abundance anomalies can be explained by a primordial pollution due to a first generation of fast rotating massive stars which released slow winds into the ISM from which a second generation of chemically anomalous stars can be formed. In particular the observed ratio of anomalous and standard stars in clusters can be used to constrain the dynamical evolution of GCs as around 95% of the standard stars need to be lost by the clusters. We show that both residual gas expulsion during the cluster formation and long term evolution are needed to achieve this ratio.

A significant fraction of stars in globular clusters (about 70%-85%) exhibit peculiar chemical patterns with strong abundance variations in light elements along with constant abundances in heavy elements. These abundance anomalies can be created in the H-burning core of a first generation of fast rotating massive stars and the corresponding elements are convoyed to the stellar surface thanks to rotational induced mixing. If the rotation of the stars is fast enough this matter is ejected at low velocity through a mechanical wind at the equator. It then pollutes the ISM from which a second generation of chemically anomalous stars can be formed. The proportion of anomalous to normal star observed today depends on at least two quantities : (1) the number of polluter stars; (2) the dynamical history of the cluster which may lose during its lifetime first and second generation stars in different proportions. Here we estimate these proportions based on dynamical models for globular clusters. When internal dynamical evolution and dissolution due to tidal forces are accounted for, starting from an initial fraction of anomalous stars of 10% produces a present day fraction of about 25%, still too small with respect to the observed 70-85%. In case gas expulsion by supernovae is accounted for, much higher fraction is expected to be produced. In this paper we also address the question of the evolution of the second generation stars that are He-rich, and deduce consequences for the age determination of globular clusters.

We investigate the long-term dynamical evolution of two distinct stellar populations of low-mass stars in globular clusters in order to study whether the energy equipartition process can explain the high number of stars harbouring abundance anomalies seen in globular clusters. We analyse N-body models by artificially dividing the low-mass stars (m<0.9 Msun) into two populations: a small number of stars (second generation) consistent with an invariant IMF and with low specific energies initially concentrated towards the cluster-centre mimic stars with abundance anomalies. These stars form from the slow winds of fast-rotating massive stars. The main part of low-mass (first generation) stars has the pristine composition of the cluster. We study in detail how the two populations evolve under the influence of two-body elaxation and the tidal forces due to the host galaxy.Stars with low specific energy initially concentrated toward the cluster centre need about two relaxation times to achieve a complete homogenisation throughout the cluster. For realistic globular clusters, the number ratio between the two populations increases only by a factor 2.5 due to the preferential evaporation of the population of outlying first generation stars. We also find that the loss of information on the stellar orbital angular momentum occurs on the same timescale as spatial homogenisation.

Various studies have established that the dynamical M/L ratios of ultra-compact dwarf galaxies (UCDs) tend to be at the limit or beyond the range explicable by standard stellar populations with canonical IMF. We discuss how IMF variations may account for these high M/L ratios and how observational approaches may in the future allow to discriminate between those possibilities. We also briefly discuss the possibility of dark matter in UCDs.

Ultra-compact dwarf galaxies (UCDs) are stellar systems with masses of around 10^7 to 10^8 Msun and half mass radii of 10-100 pc. They have some properties in common with massive globular clusters, however dynamical mass estimates have shown that UCDs have mass-to-light ratios which are on average about twice as large than those of globular clusters at comparable metallicity, and tend to be larger than what one would expect for old stellar systems with standard mass functions. One possible explanation for elevated high mass-to-light ratios in UCDs is the existence of a substantial amount of dark matter, which could have ended up in UCDs if they are the remnant nuclei of tidally stripped dwarf galaxies. Tidal stripping of dwarf galaxies has also been suggested has the origin of several massive globular clusters like Omega Cen, in which case globular clusters could have also formed with substantial amounts of dark matter. In this paper, we present collisional N-body simulations which study the co-evolution of a system composed out of stars and dark matter. We find that the dark matter gets removed from the central regions of such systems due to dynamical friction and mass segregation of stars. The friction timescale is significantly shorter than a Hubble time for typical globular clusters, while most UCDs have friction times much longer than a Hubble time. Therefore, a significant dark matter fraction remains within the half-mass radius of present-day UCDs, making dark matter a viable explanation for the elevated M/L ratios of UCDs. If at least some globular clusters formed in a way similar to UCDs, we predict a substantial amount of dark matter in their outer parts.

We study the escape rate, dN/dt, from clusters with different radii in a tidal field using analytical predictions and direct N-body simulations. We find that dN/dt depends on the ratio R=r_h/r_j, where r_h is the half-mass radius and r_j the radius of the zero-velocity surface. For R>0.05, the "tidal regime", there is almost no dependence of dN/dt on R. To first order this is because the fraction of escapers per relaxation time, t_rh, scales approximately as R^1.5, which cancels out the r_h^1.5 term in t_rh. For R<0.05, the "isolated regime", dN/dt scales as R^-1.5. Clusters that start with their initial R, Ri, in the tidal regime dissolve completely in this regime and their t_dis is insensitive to the initial r_h. We predicts that clusters that start with Ri<0.05 always expand to the tidal regime before final dissolution. Their t_dis has a shallower dependence on Ri than what would be expected when t_dis is a constant times t_rh. For realistic values of Ri, the lifetime varies by less than a factor of 1.5 due to changes in Ri. This implies that the "survival" diagram for globular clusters should allow for more small clusters to survive. We note that with our result it is impossible to explain the universal peaked mass function of globular cluster systems by dynamical evolution from a power-law initial mass function, since the peak will be at lower masses in the outer parts of galaxies. Our results finally show that in the tidal regime t_dis scales as N^0.65/w, with w the angular frequency of the cluster in the host galaxy. [ABRIDGED]

Over the last 15 years, around a hundred very young stars have been observed in the central parsec of our Galaxy. While the presence of young stars forming one or two stellar disks at approx. 0.1 pc from the supermassive black hole (SMBH) can be understood through star formation in accretion disks, the origin of the S stars observed a factor of 10 closer to the SMBH has remained a major puzzle. Here we show the S stars to be a natural consequence of dynamical interaction of two stellar disks at larger radii. Due to precession and Kozai interaction, individual stars achieve extremely high eccentricities at random orientation. Stellar binaries on such eccentric orbits are disrupted due to close passages near the SMBH, leaving behind a single S star on a much tighter orbit. The remaining star may be ejected from the vicinity of the SMBH, thus simultaneously providing an explanation for the observed hypervelocity stars in the Milky Way halo.

We have obtained high-resolution spectra of 23 ultra-compact dwarf galaxies (UCDs) in the Fornax cluster with -10.4>M_V>-13.5 mag (10^6<M/M_*<10^8), using FLAMES/Giraffe at the VLT. This is the largest homogeneous data set of UCD internal dynamics assembled to date. We derive dynamical M/L ratios for 15 UCDs covered by HST imaging. In the M_V-sigma plane, UCDs with M_V<-12 mag are consistent with the extrapolated Faber-Jackson relation for luminous ellipticals, while fainter UCDs are closer to the extrapolated globular cluster (GC) relation. At a given metallicity, Fornax UCDs have on average 30-40% lower M/L ratios than Virgo UCDs, suggesting possible differences in age or dark matter content between Fornax and Virgo UCDs. For our sample of Fornax UCDs we find no significant correlation between M/L ratio and mass. We combine our data with available M/L ratio measurements of compact stellar systems with 10^4<M/M_*<10^8, and normalise all M/L estimates to solar metallicity. We find that UCDs (M > 2*10^6 M_*) have M/L ratios twice as large as GCs (M < 2*10^6 M_*). We show that stellar population models tend to under-predict dynamical M/L ratios of UCDs and over-predict those of GCs. Considering the scaling relations of stellar spheroids, UCDs align well along the 'Fundamental Manifold', constituting the small-scale end of the galaxy sequence. The alignment for UCDs is especially clear for r_e >~ 7 pc, which corresponds to dynamical relaxation times that exceed a Hubble time. In contrast, GCs exhibit a broader scatter and do not appear to align along the manifold. We argue that UCDs are the smallest dynamically un-relaxed stellar systems, with M > 2*10^6 M_* and 7<r_e<100 pc. Future studies should aim at explaining the elevated M/L ratios of UCDs and the environmental dependence of their properties.

A novel way of looking at the evolution of star clusters is presented. With a dynamical temperature, given by the mean kinetic energy of the cluster stars, and a dynamical luminosity, which is defined as the kinetic energy of the stars leaving the cluster in analogy to the energy of photons emitted by a star, the dissolution of star clusters is studied using a new dynamical temperature-luminosity diagram for star clusters. The investigation contains a parameter-space study of open clusters of up to N = 32768 single-mass stars with different initial density distributions, half-mass radii, tidal conditions and binary fractions. The clusters show a strong correlation between dynamical temperature and dynamical luminosity and most of the investigated cluster families share a common sequence in such a dynamical temperature-luminosity diagram. Deviations from this sequence are analyzed and discussed. After core collapse, the position of a cluster within this diagram can be defined by three parameters: the mass, the tidal conditions and the binary fraction. Due to core collapse all initial conditions are lost and the remaining stars adjust to the given tidal conditions. Binaries as internal energy sources influence this adjustment. A further finding concerns the Lagrange radii of star clusters: Throughout the investigated parameter space nearly all clusters show a constant half-mass radius for the time after core collapse until dissolution. Furthermore, the ratio of half-mass radius to tidal radius evolves onto a common sequence which only depends on the mass left in the cluster.

We have studied the dissolution of initially mass segregated and unsegregated star clusters due to two-body relaxation in external tidal fields, using Aarseth's collisional N-body code NBODY4 on GRAPE6 special-purpose computers. When extrapolating results of initially not mass segregated models to globular clusters, we obtain a correlation between the time until destruction and the slope of the mass function, in the sense that globular clusters which are closer to dissolution are more strongly depleted in low-mass stars. This correlation fits observed mass functions of most globular clusters. The mass functions of several globular clusters are however more strongly depleted in low-mass stars than suggested by these models. Such strongly depleted mass functions can be explained if globular clusters started initially mass segregated. Primordial mass segregation also explains the correlation between the slope of the stellar mass function and the cluster concentration which was recently discovered by De Marchi et al. (2007). In this case, it is possible that all globular clusters started with a mass function similar to that seen in young open clusters in the present-day universe, at least for stars below m=0.8 Msun. This argues for a near universality of the mass function for different star formation environments and metallicities in the range -2 < [Fe/H] < 0. We finally describe a novel algorithm which can initialise stationary mass segregated clusters with arbitrary density profile and amount of mass segregation.

We have investigated the effect of initial mass segregation on the runaway merging of stars. The evolution of multi-mass, dense star clusters was followed by means of direct N-body simulations of up to 131.072 stars. All clusters started from King models with dimensionless central potentials of 3.0 <= W_0 <= 9.0. Initial mass segregation was realized by varying the minimum mass of a certain fraction of stars whose either (1) distances were closest to the cluster center or (2) total energies were lowest. The second case is more favorable to promote the runaway merging of stars by creating a high-mass core of massive, low-energy stars. Initial mass segregation could decrease the central relaxation time and thus help the formation of a high-mass core. However, we found that initial mass segregation does not help the runaway stellar merger to happen if the overall mass density profile is kept constant. This is due to the fact that the collision rate of stars is not increased due to initial mass segregation. Our simulations show that initial mass segregation is not sufficient to allow runaway merging of stars to occur in clusters with central densities typical for star clusters in the Milky Way.

We have developed a new method for post-Newtonian, high-precision integration of stellar systems containing a super-massive black hole (SMBH), splitting the forces on a particle between a dominant central force and perturbations. We used this method to perform fully collisional N-body simulations of inspiralling intermediate-mass black holes (IMBHs) in the centre of the Milky Way. We considered stellar cusps of different power-law indices and analysed the effects of IMBHs of different masses, all starting from circular orbits at an initial distance of 0.1 pc. Our simulations show how IMBHs deplete the central cusp of stars, leaving behind a flatter cusp with slope consistent with what has recently been observed. If an additional IMBH spirals into such a flat cusp, it can take 50 Myr or longer to merge with the central SMBH, thus allowing for direct observation in the near future. The final merger of the two black holes involves gravitational wave radiation which may be observable with planned gravitational wave detectors. Furthermore, our simulations reveal detailed properties of the hypervelocity stars (HVSs) created, and how generations of HVSs can be used to trace IMBHs in the Galactic Centre. We find that significant rotation of HVSs (which would be evidence for an IMBH) can only be expected among very fast stars (v > 1000 km/s). Also, the probability of creating a hypervelocity binary star is found to be very small.

Recently de Marchi, Paresce & Pulone (2007) studied a sample of twenty globular clusters and found that all clusters with high concentrations have steep stellar mass-functions while clusters with low concentration have comparatively shallow mass-functions. No globular clusters were found with a flat mass-function and high concentration. This seems curious since more concentrated star clusters are believed to be dynamically more evolved and should have lost more low-mass stars via evaporation, which would result in a shallower mass-function in the low-mass part. We show that this effect can be explained by residual-gas expulsion from initially mass-segregated star clusters, and is enhanced further through unresolved binaries. If gas expulsion is the correct mechanism to produce the observed trend, then observation of these parameters would allow to constrain cluster starting conditions such as star formation efficiency and the time-scale of gas expulsion.

We explore how the expulsion of gas from star-cluster forming cloud-cores due to supernova explosions affects the shape of the initial cluster mass function, that is, the mass function of star clusters when effects of gas expulsion are over. We demonstrate that if the radii of cluster-forming gas cores are roughly constant over the core mass range, as supported by observations, then more massive cores undergo slower gas expulsion. Therefore, for a given star formation efficiency, more massive cores retain a larger fraction of stars after gas expulsion. The initial cluster mass function may thus differ from the core mass function substantially, with the final shape depending on the star formation efficiency. A mass-independent star formation efficiency of about 20 per cent turns a power-law core mass function into a bell-shaped initial cluster mass function, while mass-independent efficiencies of order 40 per cent preserve the shape of the core mass function.

Mass segregation stands as one of the most robust features of the dynamical evolution of self-gravitating star clusters. In this paper we formulate parametrised models of mass segregated star clusters in virial equilibrium. To this purpose we introduce mean inter-particle potentials for statistically described unsegregated systems and suggest a single-parameter generalisation of its form which gives a mass-segregated state. We describe an algorithm for construction of appropriate star cluster models. Their stability over several crossing-times is verified by following the evolution by means of direct N-body integration.

We present new results on the evolution of the mass function of the globular cluster system of the Milky Way, taking the effect of residual gas expulsion into account. We assume that gas embedded star clusters start with a power-law mass function with slope \beta=2. The dissolution of the clusters is then studied under the combined influence of residual gas expulsion driven by energy feedback from massive stars, stellar mass-loss, two-body relaxation and an external tidal field. The influence of residual gas expulsion is studied by applying results from a large grid of N-body simulations computed by Baumgardt & Kroupa (2007). In our model, star clusters with masses less than 10^5 M_sun lose their residual gas on timescales much shorter than their crossing time and residual gas expulsion is the main dissolution mechanism for star clusters, destroying about 95% of all clusters within a few 10s of Myr. We find that in this case the final mass function of globular clusters is established mainly by the gas expulsion and therefore nearly independent of the strength of the external tidal field, and that a power-law mass function for the gas embedded star clusters is turned into a present-day log-normal one. Another consequence of residual gas expulsion and the associated strong infant mortality of star clusters is that the Galactic halo stars come from dissolved star clusters. Since field halo stars would come mainly from low-mass, short-lived clusters, our model provides an explanation for the observed abundance variations of light elements among globular cluster stars and the absence of such variations among the halo field stars.

The dissolution time (t_dis) of clusters in a tidal field does not scale with the ``classical'' expression for the relaxation time. First, the scaling with N, and hence cluster mass, is shallower due to the finite escape time of stars. Secondly, the cluster half-mass radius is of little importance. This is due to a balance between the relative tidal field strength and internal relaxation, which have an opposite effect on t_dis, but of similar magnitude. When external perturbations, such as encounters with giant molecular clouds (GMC) are important, t_dis for an individual cluster depends strongly on radius. The mean dissolution time for a population of clusters, however, scales in the same way with mass as for the tidal field, due to the weak dependence of radius on mass. The environmental parameters that determine t_dis are the tidal field strength and the density of molecular gas. We compare the empirically derived t_dis of clusters in six galaxies to theoretical predictions and argue that encounters with GMCs are the dominant destruction mechanism. Finally, we discuss a number of pitfalls in the derivations of t_dis from observations, such as incompleteness, with the cluster system of the SMC as particular example.

Dynamical mass estimates of ultra-compact dwarfs galaxies and massive globular clusters in the Fornax and Virgo clusters and around the giant elliptical Cen A have revealed some surprising results: 1) above about 10^6 M_sun the mass-to-light (M/L) ratio increases with the objects' mass; 2) some UCDs/massive GCs show high M/L values (4 to 6) that are not compatible with standard stellar population models; and 3) in the luminosity-velocity dispersion diagram, UCDs deviate from the well defined relation of `normal' GCs, being more in line with the Faber-Jackson relation of early-type galaxies. In this contribution, we present the observational evidences for high mass-to-light ratios of UCDs and discuss possible explanations for them.

Two aspects of our recent N-body studies of star clusters are presented: (1) What impact does mass segregation and selective mass loss have on integrated photometry? (2) How well compare results from N-body simulations using NBODY4 and STARLAB/KIRA?

An exciting recent finding regarding scaling relations among globular clusters is the so-called 'blue tilt': clusters of the blue sub-population follow a trend of redder colour with increasing luminosity. In this paper we evaluate to which extent field star capture over a Hubble time may explain the 'blue tilt'. We perform collisional N-body simulations to quantify the amount of field star capture occuring over a Hubble time to star clusters with 10^3 to 10^6 stars. In the simulations we follow the orbits of field stars passing through a star cluster and calculate the energy change that the field stars experience due to gravitational interaction with cluster stars during one passage through the cluster. The capture condition is that their total energy after the passage is smaller than the gravitational potential at the cluster's tidal radius. By folding this with the fly-by rates of field stars with an assumed space density as in the solar neighbourhood and a range of velocity dispersions, we derive estimates on the mass fraction of captured field stars as a function of environment. We find that integrated over a Hubble time, the ratio between captured field stars and total number of clusters stars is very low (<~ 10^(-4)), even for the smallest considered field star velocity dispersion sigma=15 km/s. This holds for star clusters in the mass range of both open clusters and globular clusters. We furthermore show that tidal friction has a negligible effect on the energy distribution of field stars after interaction with the cluster. We conclude that field star capture is not a probable mechanism for creating the colour-magnitude trend of metal-poor globular clusters.

We have carried out a large set of N-body simulations studying the effect of residual-gas expulsion on the survival rate and final properties of star clusters. We have varied the star formation efficiency, gas expulsion timescale and strength of the external tidal field, obtaining a three-dimensional grid of models which can be used to predict the evolution of individual star clusters or whole star cluster systems by interpolating between our runs. The complete data of these simulations is made available on the Internet. Our simulations show that cluster sizes, bound mass fraction and velocity profile are strongly influenced by the details of the gas expulsion. Although star clusters can survive star formation efficiencies as low as 10% if the tidal field is weak and the gas is removed only slowly, our simulations indicate that most star clusters are destroyed or suffer dramatic loss of stars during the gas removal phase. Surviving clusters have typically expanded by a factor 3 or 4 due to gas removal, implying that star clusters formed more concentrated than as we see them today. Maximum expansion factors seen in our runs are around 10. If gas is removed on timescales smaller than the initial crossing time, star clusters acquire strongly radially anisotropic velocity dispersions outside their half-mass radii. Observed velocity profiles of star clusters can therefore be used as a constraint on the physics of cluster formation.

We investigate the evolution of collisionally merged stars with mass of ~100 MSun which might be formed in dense star clusters. We assumed that massive stars with several tens Msun collide typically after ~1Myr of the formation of the cluster and performed hydrodynamical simulations of several collision events. Our simulations show that after the collisions, merged stars have extended envelopes and their radii are larger than those in the thermal equilibrium states and that their interiors are He-rich because of the stellar evolution of the progenitor stars. We also found that if the mass-ratio of merging stars is far from unity, the interior of the merger product is not well mixed and the elemental abundance is not homogeneous. We then followed the evolution of these collision products by a one dimensional stellar evolution code. After an initial contraction on the Kelvin-Helmholtz (thermal adjustment) timescale (~1000-10000 yr), the evolution of the merged stars traces that of single homogeneous stars with corresponding masses and abundances, while the initial contraction phase shows variations which depend on the mass ratio of the merged stars. We infer that, once runaway collisions have set in, subsequent collisions of the merged stars take place before mass loss by stellar winds becomes significant. Hence, stellar mass loss does not inhibit the formation of massive stars with mass of ~1000Msun.

(Simplified) We present a theoretical framework which establishes how the core radius of a star cluster varies with the mass of an assumed central black hole. Our result is that the ratio of core to half-mass radius varies as the 3/4 power of the ratio of the black hole to cluster mass, when the system is well relaxed. The theory compares favourably with a number of simulations of this problem, which extend to black hole masses of order 10% of the cluster mass. Though strictly limited as yet to clusters with stars of equal mass, our conclusion strengthens the view that clusters with large core radii are the most promising candidates in which to find a massive black hole.

Observations indicate that present-day star formation in the Milky Way disk takes place in stellar ensembles or clusters rather than in isolation. Bound, long lived stellar groups are known as open clusters. They gradually lose stars and in their final evolutionary stages they are severely disrupted leaving an open cluster remnant made of a few stars. In this paper, we study in detail the stellar content and kinematics of the poorly populated star cluster NGC1901. This object appears projected against the Large Magellanic Cloud. The aim of the present work is to derive the current evolutionary status, binary fraction, age and mass of this stellar group. These are fundamental quantities to compare with those from N-body models in order to study the most general topic of star cluster evolution and dissolution.The analysis is performed using wide-field photometry in the UBVI pass-band, proper motions from the UCAC.2 catalog, and 3 epochs of high resolution spectroscopy, as well as results from extensive N-body calculations.The star group NGC1901 is found to be an ensemble of solar metallicity stars, 400+/-100 Myr old, with a core radius of 0.23 pc, a tidal radius of 1.0 pc, and located at 400+/-50 pc from the Sun. Out of 13 confirmed members, only 5 single stars have been found. Its estimated present-day binary fraction is at least 62%. The calculated heliocentric space motion of the cluster is not compatible with possible membership in the Hyades stream.Our results show that NGC1901 is a clear prototype of open cluster remnant characterized by a large value of the binary fraction and a significant depletion of low-mass stars. In the light of numerical simulations, this is compatible with NGC1901 being what remains of a larger system initially made of 500-750 stars.

Aims. We determine masses and mass-to-light ratios of five ultra-compact dwarf galaxies and one dwarf elliptical nucleus in the Fornax cluster from high resolution spectroscopy. Methods. Velocity dispersions were derived from selected wavelength regions using a direct-fitting method. To estimate the masses of the UCDs a new modelling program has been developed that allows a choice of different representations of the surface brightness profile (i.e. Nuker, Sersic or King laws) and corrects the observed velocity dispersions for observational parameters (i.e. seeing, slit size). Results. The observed velocity dispersions range between 22 and 30 km/s. The resulting masses are between 1.8 and 9.5x10^7M_sun. These, as well as the central and global projected velocity dispersions, were derived from the generalized King model which turned out to give the most stable results. The masses of two UCDs, that are best fitted by a two-component profile, were derived from a combined King+Sersic model. The mass-to-light ratios of the Fornax UCDs range between 3 and 5 (M/L_V)_sun. Within 1-2 half-mass radii dark matter is not dominating UCDs. Conclusions. We show that the mass-to-light ratios of UCDs in Fornax are consistent with those expected for pure stellar populations. Thus UCDs seem to be the result of cluster formation processes within galaxies rather than being compact dark matter dominated substructures themselves. Whether UCDs gained their mass in super-star cluster complexes of mergers or in nuclear star cluster formation processes remains an open question. It appears, however, clear that star clusters more massive than about 5\times10^6M_sun exhibit a more complex formation history than the less massive `ordinary' globular clusters.

This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discussion at the meeting. Our meeting, which had close to 300 registered participants, was broadly aimed at the large community of astronomers around the world working on the formation and evolution of compact objects and interacting binary systems in dense star clusters, such as globular clusters and galactic nuclei. Great advances have occurred in this field during the past few years, including the introduction of fundamentally new theoretical paradigms for the formation and evolution of compact objects in binaries as well as countless new discoveries by astronomers that have challenged many accepted models. Some of the highlights include: a nearly complete census of all the millisecond pulsars in 47 Tuc; first detections of many new radio pulsars in other clusters, particularly Terzan 5; detailed studies of X-ray binary populations and their luminosity functions in many galaxies and extragalactic globular clusters; increasing evidence for intermediate-mass black holes in clusters and greatly improved theoretical understanding of their possible formation processes.

We have performed N-body simulations of the formation of hyper-velocity stars (HVS) in the centre of the Milky Way due to inspiralling intermediate-mass black holes (IMBHs). We considered IMBHs of different masses, all starting from circular orbits at an initial distance of 0.1 pc. We find that the IMBHs sink to the centre of the Galaxy due to dynamical friction, where they deplete the central cusp of stars. Some of these stars become HVS and are ejected with velocities sufficiently high to escape the Galaxy. Since the HVS carry with them information about their origin, in particular in the moment of ejection, the velocity distribution and the direction in which they escape the Galaxy, detecting a population of HVS will provide insight in the ejection processes and could therefore provide indirect evidence for the existence of IMBHs. Our simulations show that HVS are generated in short bursts which last only a few Myrs until the IMBH is swallowed by the supermassive black hole (SMBH). HVS are ejected almost isotropically, which makes IMBH induced ejections hard to distinguish from ejections due to encounters of stellar binaries with a SMBH. After the HVS have reached the galactic halo, their escape velocities correlate with the distance from the Galactic centre in the sense that the fastest HVS can be found furthest away from the centre. The velocity distribution of HVS generated by inspiralling IMBHs is also nearly independent of the mass of the IMBH and can be quite distinct from one generated by binary encounters. Finally, our simulations show that the presence of an IMBH in the Galactic centre changes the stellar density distribution inside r<0.02 pc into a core profile, which takes at least 100 Myrs to replenish.

We investigate encounters between giant molecular clouds (GMCs) and star clusters. We propose a single expression for the energy gain of a cluster due to an encounter with a GMC, valid for all encounter distances and GMC properties. This relation is verified with N-body simulations of cluster-GMC encounters and excellent agreement is found. The fractional mass loss from the cluster is 0.25 times the fractional energy gain. This is because 75% of the injected energy goes to the velocities of escaping stars, that are higher than the escape velocity. We derive an expression for the cluster disruption time (t_dis) based on the mass loss from the simulations, taking into account the effect of gravitational focusing by the GMC. The disruption time depends on the cluster mass (M_c) and half-mass radius (r_h) as t_dis=2.0 S (M_c/10^4 M_sun)(3.75 pc/r_h)^3 Gyr, with S=1 for the solar neighbourhood and inversely proportional to the GMC density. The observed shallow relation between cluster radius and mass gives t_dis a power-law dependence on the mass with index 0.7, similar to that found from observations and from simulations of clusters dissolving in tidal fields (0.62). The constant of 2.0 Gyr is about a factor of 3.5 shorter than found from earlier simulations of clusters dissolving under the combined effect of galactic tidal field and stellar evolution. It is somewhat higher than the observationally determined value of 1.3 Gyr. It suggests, however, that the combined effect of tidal field and encounters with GMCs can explain the lack of old open clusters in the solar neighbourhood. GMC encounters can also explain the (very) short disruption time that was observed for star clusters in the central region of M51, since there rho_n is an order of magnitude higher than in the solar neighbourhood.

Globular clusters are among the oldest objects in galaxies, and understanding the details of their formation and evolution can bring valuable insight into the early history of galaxies. This review summarises the current knowledge about the dissolution of star clusters and discusses the implications of star cluster dissolution for the evolution of the mass function of star cluster systems in galaxies.

We present the methods and preparatory work for our study of the collisional runaway scenario to form a very massive star (VMS, M > 400 Msun) at the centre of a young, compact stellar cluster. In the first phase of the process, a very dense central core of massive stars (M =~ 30-120 Msun) forms through mass segregation and gravothermal collapse. This leads to a collisional stage, likely to result in the formation of a VMS (itself a possible progenitor for an intermediate-mass black hole) through a runaway sequence of mergers between the massive stars. In this paper we present the runaway scenario in a general astrophysical context. We then explain the numerical method used to investigate it. Our approach is based on a Monte Carlo code to simulate the stellar dynamics of spherical star clusters using a very large number of particles (a few 1e5 to several 1e6). Finally, we report on test computations carried out to ensure that our implementation of the important physics is sound. In a second paper, we present results from more than 100 cluster simulations realized to determine the conditions leading to the collisional formation of a VMS and the characteristics of the runaway sequences.

We simulate the inner 100pc of the Milky-Way Galaxy to study the formation and evolution of the population of star clusters and intermediate mass black holes. For this study we perform extensive direct N-body simulations of the star clusters which reside in the bulge, and of the inner few tenth of parsecs of the super massive black hole in the Galactic center. In our N-body simulations the dynamical friction of the star cluster in the tidal field of the bulge are taken into account via (semi)analytic soluations. The N-body calculations are used to calibrate a (semi)analytic model of the formation and evolution of the bulge. We find that about 10% of the clusters born within 100pc of the Galactic center undergo core collapse during their inward migration and form intermediate-mass black holes (IMBHs) via runaway stellar merging. After the clusters dissolve, these IMBHs continue their inward drift, carrying a few of the most massive stars with them. We predict that region within about 10 parsec of the SMBH is populated by about 50IMBHs of some 1000Msun. Several of these are expected to be accompanied still by some of the most massive stars from the star cluster. We also find that within a few milliparsec of the SMBH there is a steady population of several IMBHs. This population drives the merger rate between IMBHs and the SMBH at a rate of about one per 10Myr, sufficient to build the accumulate majority of mass of the SMBH. Mergers of IMBHs with SMBHs throughout the universe are detectable by LISA, at a rate of about two per week.

Recent X-ray observations and theoretical modelling have made it plausible that some ultraluminous X-ray sources (ULX) are powered by intermediate-mass black holes (IMBHs). N-body simulations have also shown that runaway merging of stars in dense star clusters is a way to form IMBHs. In the present paper we have performed N-body simulations of young clusters such as MGG-11 of M82 in which IMBHs form through runaway merging. We took into account the effect of tidal heating of stars by the IMBH to study the tidal capture and disruption of stars by IMBHs. Our results show that the IMBHs have a high chance of capturing stars through tidal heating within a few core relaxation times and we find that 1/3 of all runs contain a ULX within the age limits of MGG-11, a result consistent with the fact that a ULX is found in this galaxy. Our results strengthen the case for some ULX being powered by intermediate-mass black holes.

We present a simple analytical description of the disruption of star clusters in a tidal field, which agrees excellently with detailed N-body simulations. The analytic expression can be used to predict the mass and age histograms of surviving clusters for any cluster initial mass function and any cluster formation history. The method is applied to open clusters in the solar neighbourhood, based on the new cluster sample of Kharchenko et al. From a comparison between the observed and predicted age distributions in the age range between 10 Myr to 3 Gyr we find the following results: (1) The disruption time of a 10^4 M_sun cluster in the solar neighbourhood is about 1.3+/-0.5 Gyr. This is a factor 5 shorter than derived from N-body simulations of clusters in the tidal field of the galaxy. (2) The present starformation rate in bound clusters within 600 pc from the Sun is 5.9+/-0.8 * 10^2 M_sun / Myr, which corresponds to a surface star formation rate in bound clusters of 5.2+/-0.7 10^(-10) M_sun/yr/pc^2. (3) The age distribution of open clusters shows a bump between 0.26 and 0.6 Gyr when the cluster formation rate was 2.5 times higher than before and after. (4) The present star formation rate in bound clusters is half as small as that derived from the study of embedded clusters. The difference suggests that half of the clusters in the solar neighbourhood become unbound within 10 Myr. (5) The most massive clusters within 600 pc had an initial mass of 3*10^4 M_sun. This is in agreement with the statistically expected value based on a cluster initial mass function with a slope of -2, even if the physical upper mass limit is as high as 10^6 M_sun.

We propose to determine the stellar velocity dispersions of globular clusters in the outer halo of the Milky Way in order to decide whether the dynamics of the universe on large scales is governed by dark matter or modified Newtonian dynamics (MOND). We show that for a number of galactic globular clusters, both the internal and the external accelerations are significantly below the critical acceleration parameter $a_0$ of MOND. This leads to velocity dispersions in case of MOND which exceed their Newtonian counterparts by up to a factor of 3, providing a stringent test for MOND. Alternatively, in case high velocity dispersions are found, these would provide the first evidence that globular clusters are dark matter dominated.

Close encounters and physical collisions between stars in young dense clusters can result in new channels for stellar evolution, and may lead to the formation of very massive stars and black holes via runaway merging. We present some details of this process, using the results of N-body simulations and simple analytical estimates to place limits on the cluster parameters for which it expected to occur. For small clusters, the mass of the runaway is effectively limited by the total number of high-mass stars in the system. For larger clusters, the runaway mass is determined by the fraction of stars that can mass-segregate to the cluster core while still on the main sequence. In typical cases, the result is in the range commonly cited for intermediate-mass black holes. This mechanism may therefore have important implications for the formation of massive black holes and black-hole binaries in dense cluster cores.

It has been assumed that intermediate-mass black holes (IMBHs) in globular clusters can only reside in the most centrally concentrated clusters, with a so-called `core-collapsed' density profile. While this would be a natural guess, it is in fact wrong. We have followed the evolution of star clusters containing IMBHs with masses between 125 \le M_{BH} \le 1000 M_{\odot} through detailed N-body simulations, and we find that a cluster with an IMBH, in projection, appears to have a relatively large `core' with surface brightness only slightly rising toward the center. This makes it highly unlikely that any of `core-collapsed' clusters will harbor an IMBH. On the contrary, the places to look for an IMBH are those clusters that can be fitted well by medium-concentration King models. The velocity dispersion of the visible stars in a globular cluster with an IMBH is nearly constant well inside the apparent core radius. For a cluster of mass M_C containing an IMBH of mass M_{BH}, the influence of the IMBH becomes significant only at a fraction 2.5 M_{BH}/M_C of the half-mass radius, deep within the core, where it will affect only a small number of stars. In conclusion, observational detection of an IMBH may be possible, but will be challenging.

In this paper we report results of collisional N-body simulations of the dynamical evolution of equal-mass star clusters containing a massive central black hole. Each cluster is composed of between 5,000 to 180,000 stars together with a central black hole which contains between 0.2% to 10% of the total cluster mass. We find that for large enough black hole masses, the central density follows a power-law distribution with slope \rho \sim r^{-1.75} inside the radius of influence of the black hole, in agreement with predictions from earlier Fokker Planck and Monte Carlo models. The tidal disruption rate of stars is within a factor of two of that derived in previous studies. It seems impossible to grow an intermediate-mass black hole (IMBH) from a M \le 100 Msun progenitor in a globular cluster by the tidal disruption of stars, although M = 10^3 Msun IMBHs can double their mass within a Hubble time in dense globular clusters. The same is true for the supermassive black hole at the centre of the Milky Way. Black holes in star clusters will feed mainly on stars tightly bound to them and the re-population of these stars causes the clusters to expand, reversing core-collapse without the need for dynamically active binaries. Close encounters of stars in the central cusp also lead to an increased mass loss rate in the form of high-velocity stars escaping from the cluster. A companion paper will extend these results to the multi-mass case.

We have followed the evolution of multi-mass star clusters containing massive central black holes through collisional N-body simulations done on GRAPE6. Each cluster is composed of between 16,384 to 131,072 stars together with a black hole with an initial mass of M_BH=1000 Msun. We follow the evolution of the clusters under the combined influence of two-body relaxation, stellar mass-loss and tidal disruption of stars. The (3D) mass density profile follows a power-law distribution \rho \sim r^{-\alpha} with slope \alpha=1.55. This leads to a constant density profile of bright stars in projection, which makes it highly unlikely that core collapse clusters contain intermediate-mass black holes (IMBHs). Instead globular clusters containing IMBHs can be fitted with standard King profiles. The disruption rate of stars is too small to form an IMBH out of a M_BH \approx 50 Msun progenitor black hole, unless a cluster starts with a central density significantly higher than what is seen in globular clusters. Kinematical studies can reveal 1000 Msun IMBHs in the closest clusters. IMBHs in globular clusters are only weak X-ray sources since the tidal disruption rate of stars is low and the star closest to the IMBH is normally another black hole. For globular clusters, dynamical evolution can push compact stars near the IMBH to distances small enough that they become detectable through gravitational radiation. If 10% of all globular clusters contain IMBHs, extragalactic globular clusters could be one of the major sources for {\it LISA}. (abridged)

Observations and N-body simulations both support a simple relation for the disruption time of a cluster as a function of its mass of the form: t_dis = t_4 * (M/10^4 Msun)^gamma. The scaling factor t_4 seems to depend strongly on the environment. Predictions and observations show that gamma ~ 0.64 +/- 0.06. Assuming that t_dis ~ M^0.64 is caused by evaporation and shocking implies a relation between the radius and the mass of a cluster of the form: r_h ~ M^0.07, which has been observed in a few galaxies. The suggested relation for the disruption time implies that the lower mass end of the cluster initial mass function will be disrupted faster than the higher mass end, which is needed to evolve a young power law shaped mass function into the log-normal mass function of old (globular) clusters.

We have followed the evolution of multi-mass star clusters containing massive central black holes by N-body simulations on the GRAPE6 computers of Tokyo University. We find a strong cluster expansion and significant structural changes of the clusters. Star clusters with IMBHs have power-law density profiles $rho ~ r^{-alpha}$ with slopes $alpha=1.55$ inside the influence sphere of the central black hole. This leads to a constant density profile of bright stars in projection, which rules out the presence of intermediate mass black holes in core collapse clusters. If the star clusters are surrounded by a tidal field, a central IMBH speeds up the destruction of the cluster until a remnant of a few hundred stars remains, which stays bound to the IMBH for a long time. We also discuss the efficiency of different detection mechanisms for finding IMBHs in star clusters.

We discuss the structural change and degree of mass segregation of young dense star clusters within about 100pc of the Galactic center. In our calculations, which are performed with GRAPE-6, the equations of motion of all stars and binaries are calculated accurately but the external potential of the Galaxy is solved (semi)analytically. The simulations are preformed to model the Arches star cluster. We find that star clusters with are less strongly perturbed by the tidal field and dynamical friction are much stronger affected by mass segregation; resulting in a significant pile-up of massive stars in the cluster center. At an age of about 3.5Myr more than 90 per cent of the stars more massive than ~10Msun are concentrated within the half-mass radius of the surviving cluster. Star clusters which are strongly perturbed by the tidal field of the parent Galaxy are much less affected by mass segregation.

A luminous X-ray source is associated with a cluster (MGG-11) of young stars \~200pc from the center of the starburst galaxy M82. The properties of the X-ray source are best explained by a black hole with a mass of at least 350Msun, which is intermediate between stellar-mass and supermassive black holes. A nearby but somewhat more massive star cluster (MGG-9) shows no evidence of such an intermediate mass black hole, raising the issue of just what physical characteristics of the clusters can account for this difference. Here we report numerical simulations of the evolution and the motions of stars within the clusters, where stars are allowed to mergers with each other. We find that for MGG-11 dynamical friction leads to the massive stars sinking rapidly to the center of the cluster to participate in a runaway collision, thereby producing a star of 800-3000Msun, which ultimately collapses to an black hole of intermediate mass. No such runaway occurs in the cluster MGG-9 because the larger cluster radius leads to a mass-segregation timescale a factor of five longer than for MGG-11.

We know from observations that globular clusters are very efficient catalysts in forming unusual short-period binary systems or their offspring, such as low-mass X-ray binaries (LMXBs; neutron stars accreting matter from low-mass stellar companions), cataclysmic variables (CVs; white dwarfs accreting matter from stellar companions), and millisecond pulsars (MSPs; rotating neutron stars with spin periods of a few ms). Although there has been little direct evidence, the overabundance of these objects in globular clusters has been attributed by numerous authors to the high densities in the cores, which leads to an increase in the formation rate of exotic binary systems through close stellar encounters. Many such close binary systems emit X-radiation at low luminosities (L_x < 10^{34} erg/s) and are being found in large numbers through observations with the Chandra X-ray Observatory. Here we present conclusive observational evidence for a link between the number of close binaries observed in X-rays in a globular cluster and the stellar encounter rate of the cluster. We also make an estimate of the total number of LMXBs in globular clusters in our Galaxy.

We present a comparison between the observational data on the kinematical structure of G1 in M31, obtained with the Hubble WFPC2 and STIS instruments, and the results of dynamical simulations carried out using the special-purpose computer GRAPE-6. We have obtained good fits for models starting with single cluster King-model initial conditions and even better fits when starting our simulations with a dynamically constructed merger remnant of two star clusters. In the latter case, the results from our simulations are in excellent agreement with the observed profiles of luminosity, velocity dispersion, rotation, and ellipticity. We obtain a mass-to-light ratio of $M/L = 4.0 \pm 0.4$ and a total cluster mass of $M=(8\pm 1)\times 10^6 M_\odot$. Given that our dynamical model can fit all available observational data very well, there seems to be no need to invoke the presence of an intermediate-mass black hole in the center of G1.

This paper considers the phenomenon of deep core collapse in collisional stellar systems, with stars of equal mass. The collapse takes place on some multiple, $\xi ^{-1}$, of the central relaxation time, and produces a density profile in which $\rho \propto r^{-\alpha}$, where $\alpha$ is a constant. The parameters $\alpha$ and $\xi$ have usually been determined from simplified models, such as gas and Fokker-Planck models, often with the simplification of isotropy. Here we determine the parameters directly from N-body simulations carried out using the newly completed GRAPE-6.

We present a detailed comparison between the latest observational data on the kinematical structure of the core of M15, obtained with the Hubble STIS and WFPC2 instruments, and the results of dynamical simulations carried out using the special-purpose GRAPE-6 computer. The observations imply the presence of a significant amount of dark matter in the cluster core. In our dynamical simulations, neutron stars and/or massive white dwarfs concentrate to the center through mass segregation, resulting in a sharp increase in $M/L$ toward the center. While consistent with the presence of a central black hole, the Hubble data can also be explained by this central concentration of stellar-mass compact objects. The latter interpretation is more conservative, since such remnants result naturally from stellar evolution, although runaway merging leading to the formation of a black hole may also occur for some range of initial conditions. We conclude that no central massive object is required to explain the observational data, although we cannot conclusively exclude such an object at the level of $\sim500-1000$ solar masses. Our findings are unchanged when we reduce the assumed neutron-star retention fraction in our simulations from 100% to 0%.

We report results of a large set of N-body calculations aimed to study the evolution of multi-mass star clusters in external tidal fields. Our clusters start with the same initial mass-functions, but varying particle numbers, orbital types and density profiles. Our main focus is to study how the stellar mass-function and other cluster parameters change under the combined influence of stellar evolution, two-body relaxation and the external tidal field. We find that the lifetimes of star clusters moving on similar orbits scale as T sim T_RH^x, where T_RH is the relaxation time, and the exponent x depends on the initial concentration of the cluster and is around x approx 0.75. From the results for the lifetimes, we predict that between 53% to 67% of all galactic globular clusters will be destroyed within the next Hubble time. Low-mass stars are preferentially lost and the depletion is strong enough to turn initially increasing mass-functions into mass-functions which decrease towards the low-mass end. The details of this depletion are insensitive to the starting condition of the cluster. The preferential depletion of low-mass stars from star clusters leads to a decrease of their mass-to-light ratios except for a short period close to final dissolution. The fraction of compact remnants is increasing throughout the evolution and they are more strongly concentrated towards the cluster cores than main-sequence stars. For a sample of galactic globular clusters with well observed parameters, we find a correlation between the observed slope of the mass-function and the lifetimes predicted by us. It seems possible that galactic globular clusters started with a mass-function similar to what one observes for the average mass-function of the galactic disc and bulge. (Abridged)

We report results of N-body simulations of isolated star clusters, performed up to the point where the clusters are nearly completely dissolved. Our main focus is on the post-collapse evolution of these clusters. We find that after core collapse, isolated clusters evolve along nearly a single sequence of models whose properties are independent of the initial density profile and particle number. Due to the slower expansion of high-N clusters, relaxation times become almost independent of the particle number after several core collapse times, at least for the particle range of our study. As a result, the dissolution times of isolated clusters exhibit a surprisingly weak dependence on N. We find that most stars escape due to encounters between single stars inside the half-mass radius of the cluster. Encounters with binaries take place mostly in the cluster core and account for roughly 15% of all escapers. Encounters between single stars at intermediate radii are also responsible for the build up of a radial anisotropic velocity distribution in the halo. For clusters undergoing core oscillations, escape due to binary stars is efficient only when the cluster center is in a contracted phase. Our simulations show that it takes about 10^5 N-body time units until the global anisotropy reaches its maximum value. The anisotropy increases with particle number and it seems conceivable that isolated star clusters become vulnerable to radial orbit instabilities for large enough N. However, no indication for the onset of such instabilities was seen in our runs.

Interacting galaxies like the famous Antennae (NGC 4038/4039) or Stephan's Quintet (HCG 92) show considerable star forming activity in their tidal arms. High resolution images (e.g. from HST-observations) indicate that these regions consist of up to hundreds of massive stellar clusters or tidal dwarf galaxies (TDG). In this paper we want to investigate the future fate of these clusters of massive star clusters (in this work called super-clusters). We simulate compact super-clusters in the tidal field of a host-galaxy and investigate the influence of orbital and internal parameters on the rate and timescale of the merging process. We show that it is possible that such configurations merge and build a dwarf galaxy, which could be an important mechanism of how long-lived dwarf satellite galaxies form. A detailed study of the merger object will appear in a follow-up paper.

We report on UBVRI photometry of two $5^{\prime} \times 5^{\prime}$ fields in the region of the young open cluster NGC 3324. One of our fields covers the core region, while the other is closer to the tidal radius of the cluster. Our study provides the first CCD photometry of NGC 3324. We find that the cluster is very young and probably contains several pre Main Sequence (MS) stars. 25 members are identified on the basis of their position in the (U-B) vs (B-V) diagram. We investigate the relation of the red super-giant HD 92207 with NGC 3324, suggesting that it probably does not belong to the cluster. Our second field is close to Loden 165, a possible cluster of stars that has never been studied so far. We show that this object is a probable open cluster, much older than NGC 3324 and much closer to the Sun.

We study the evolution of the galactic globular cluster system to determine its initial mass-function and the fraction of halo stars that could have come from disrupted globular clusters. We study the cluster evolution under the influence of two destruction mechanisms: Two-body relaxation and dynamical friction. New results of N-body simulations are used for the lifetimes of clusters dissolving under the influence of two-body relaxation. Two different mass-functions are studied: A gaussian initial distribution in log(M_C) with mean mass and scatter similar to what one observes for the galactic globular cluster system, and a power-law distribution resembling what is seen for clusters in merging and interacting galaxies. We find that in the inner parts of the galaxy, both distributions evolve in such a way that they are consistent with the observations. In the outer parts, a gaussian initial distribution gives the better fit. This might change however if elliptic orbits are considered, or there are undiscovered low-mass clusters in the galactic halo. In any case, only a small fraction of the stellar halo of the Milky Way originated in globular clusters.

We report results of collisional N-body simulations aimed to study the N-dependance 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 which gain the energy necessary for escape are scattered to lower energies before they can leave the cluster. Since 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.

Mean proper motions and parallaxes of 205 open clusters were determined from their member stars found in the Hipparcos Catalogue. 360 clusters were searched for possible members, excluding nearby clusters with distances D < 200 pc. Members were selected using ground based information (photometry, radial velocity, proper motion, distance from the cluster centre) and information provided by Hipparcos (proper motion, parallax). Altogether 630 certain and 100 possible members were found. A comparison of the Hipparcos parallaxes with photometric distances of open clusters shows good agreement. The Hipparcos data confirm or reject the membership of several Cepheids in the studied clusters.

The statistical test described by Wielen et al.(1994) is used to derive new zero-points of ground-based Cepheid period-luminosity (PL) and period-luminosity-colour (PLC) relations. Eleven relations are compared with the Hipparcos data. Our results argue for a typical increase of the adopted distance scale by about 8% +- 8%. Our zero-point for the PL relation of Caldwell & Laney (1991) is in agreement with that of Feast & Catchpole (1997).