Source author record

H. Baumgardt

H. Baumgardt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2020arXiv

New Insight into the stellar mass function of Galactic globular clusters

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^α$ with slopes $α>-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.

preprint2020arXiv

The WAGGS project -- III. Discrepant mass-to-light ratios of Galactic globular clusters at high metallicity

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.

preprint2019arXiv

Multiple populations in globular clusters and their parent galaxies

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.

preprint2019arXiv

Three-Component Kinematics of Multiple Stellar Populations in Globular Clusters with Gaia and VLT

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-$σ$ level. The analysis of the velocity-dispersion profiles of multiple populations confirms that 2P stars of NGC$\,$0104 show stronger anisotropy than the 1P.

preprint2015arXiv

Biases in the determination of dynamical parameters of star clusters: today and in the Gaia era

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.

preprint2013arXiv

Indication for an intermediate-mass black hole in the globular cluster NGC 5286 from kinematics

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]

preprint2012arXiv

How well do STARLAB and NBODY compare? II: Hardware and accuracy

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.

preprint2011arXiv

A Dynamical N-body Model for the Central Region of $ω$ Centauri

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 $ω$ 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 $ω$ 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 $ω$ 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)

preprint2011arXiv

Kinematic signature of an intermediate-mass black hole in the globular cluster NGC 6388

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} - σ$ 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}$.

preprint2011arXiv

Mass segregation and fractal substructure in young massive clusters: (I) the McLuster code and method calibration

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)

preprint2011arXiv

Spatially Resolved Kinematics of an Ultra-Compact Dwarf Galaxy

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.

preprint2011arXiv

The role of stellar collisions for the formation of massive stars

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.

preprint2011arXiv

Tidal disruption rate of stars by supermassive black holes obtained by direct N-body simulations

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é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)

preprint2010arXiv

A top-heavy stellar initial mass function in starbursts as an explanation for the high mass-to-light ratios of ultra compact dwarf galaxies

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.

preprint2010arXiv

Evolution of two stellar populations in globular clusters II. Effects of primordial gas expulsion

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.

preprint2010arXiv

Peculiarities in Velocity Dispersion and Surface Density Profiles of Star Clusters

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)

preprint2010arXiv

The ACS Fornax Cluster Survey. IX. The Color-Magnitude Relation of Globular Cluster Systems

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.

preprint2009arXiv

How well do STARLAB and NBODY4 compare? I: Simple models

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)

preprint2009arXiv

The photometric evolution of dissolving star clusters: II. Realistic models. Colours and M/L ratios

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)

preprint2009arXiv

The velocity dispersion and mass-to-light ratio of the remote halo globular cluster NGC 2419

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 σ= 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.

preprint2009arXiv

Tidal Tails of Star Clusters

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.

preprint2008arXiv

The ages of Galactic globular clusters in the context of self-enrichment

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.

preprint2006arXiv

Neutron Stars and Black Holes in Star 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.