The Global Mass Functions of 35 Galactic globular clusters: II. Clues on the Initial Mass Function and Black Hole Retention Fraction

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.Comment: 8 pages, 6 figures, MNRAS, 472, 74

The effect of tides on the Fornax dwarf spheroidal galaxy

Estimates of the mass distribution and dark-matter (DM) content of dwarf spheroidal galaxies (dSphs) are usually derived under the assumption that the effect of the tidal field of the host galaxy is negligible over the radial extent probed by kinematic data-sets. We assess the implications of this assumption in the specific case of the Fornax dSph by means of N-body simulations of a satellite orbiting around the Milky Way. We consider observationally-motivated orbits and we tailor the initial distributions of the satellite's stars and DM to match, at the end of the simulations, the observed structure and kinematics of Fornax. In all our simulations the present-day observable properties of Fornax are not significantly influenced by tidal effects. The DM component is altered by the interaction with the Galactic field (up to 20% of the DM mass within 1.6 kpc is lost), but the structure and kinematics of the stellar component are only mildly affected even in the more eccentric orbit (more than 99% of the stellar particles remain bound to the dwarf). In the simulations that successfully reproduce Fornax's observables, the dark-to-luminous mass ratio within 1.6 kpc is in the range 5-6, and up to 16-18 if measured within 3 kpc.Comment: 19 pages, 16 figures. Accepted for publication in MNRA

Deep multi-telescope photometry of NGC 5466. II. The radial behaviour of the mass function slope

We use a combination of data acquired with the Advanced Camera for Survey (ACS) on board the Hubble Space Telescope and the Large Binocular Camera (LBC-blue) mounted on the Large Binocular Telescope, to sample the main sequence stars of the globular cluster NGC~5466 in the mass range $0.3<M/M_\odot<0.8$. We derive the cluster's Luminosity Function in several radial regions, from the center of the cluster out to the tidal radius. After corrections for incompleteness and field-contamination, this has been compared to theoretical Luminosity Functions, obtained by multiplying a simple power law Mass Function in the form dN/dm$\propto m^{\alpha}$ by the derivative of the mass-luminosity relationship of the best-fit isochrone. We find that $\alpha$ varies from -0.6 in the core region to -1.9 in the outer region. This fact allows us to observationally prove that the stars in NGC 5466 have experienced the effects of mass segregation. We compare the radial variation of $\alpha$ from the center out to 5 core radii (r$_c$) in NGC 5466 and the globular cluster M10, finding that the gradient of $\alpha$ in the first 5r$_c$ is more than a factor of 2 shallower in NGC 5466 than in M10, in line with the differences in the clusters' relaxation timescales. NGC 5466 is dynamically younger than M10, with two-body relaxation processes only recently starting to shape the distribution of main sequence stars. This result fully agrees with the conclusion obtained in our previous works on the radial distribution of Blue Straggler Stars, further confirming that this can be used as an efficient clock to measure the dynamical age of stellar systems.Comment: Accepted for publications on Ap

Terzan 8: a Sagittarius-flavoured globular cluster

Massive globular clusters (GCs) contain at least two generations of stars with slightly different ages and clearly distinct light elements abundances. The Na-O anticorrelation is the best studied chemical signature of multiple stellar generations. Low-mass clusters appear instead to be usually chemically homogeneous. We are investigating low-mass GCs to understand what is the lower mass limit where multiple populations can form, mainly using the Na and O abundance distribution. We used VLT/FLAMES spectra of giants in the low-mass, metal-poor GC Terzan 8, belonging to the Sagittarius dwarf galaxy, to determine abundances of Fe, O, Na, alpha-, Fe-peak, and neutron-capture elements in six stars observed with UVES and 14 observed with GIRAFFE. The average metallicity is [Fe/H]=-2.27+/-0.03 (rms=0.08), based on the six high-resolution UVES spectra. Only one star, observed with GIRAFFE, shows an enhanced abundance of Na and we tentatively assign it to the second generation. In this cluster, at variance with what happens in more massive GCs, the second generation seems to represent at most a small minority fraction. We discuss the implications of our findings, comparing Terzan 8 with the other Sgr dSph GCs, to GCs and field stars in the Large Magellanic Cloud, Fornax, and in other dwarfs galaxies.Comment: 15 pages, 12 figures, 10 tables; accepted for publication on Astronomy and Astrophysic

Evolution of Binary Stars in Multiple-Population Globular Clusters - II. Compact Binaries

We present the results of a survey of N-body simulations aimed at exploring the evolution of compact binaries in multiple-population globular clusters.We show that as a consequence of the initial differences in the structural properties of the first-generation (FG) and the second-generation (SG) populations and the effects of dynamical processes on binary stars, the SG binary fraction decreases more rapidly than that of the FG population. The difference between the FG and SG binary fraction is qualitatively similar to but quantitatively smaller than that found for wider binaries in our previous investigations.The evolution of the radial variation of the binary fraction is driven by the interplay between binary segregation, ionization and ejection. Ionization and ejection counteract in part the effects of mass segregation but for compact binaries the effects of segregation dominate and the inner binary fraction increases during the cluster evolution. We explore the variation of the difference between the FG and the SG binary fraction with the distance from the cluster centre and its dependence on the binary binding energy and cluster structural parameters. The difference between the binary fraction in the FG and the SG populations found in our simulations is consistent with the results of observational studies finding a smaller binary fraction in the SG population.Comment: 9 pages, 12 figures. Accepted for publication in MNRA

The eye of Gaia on globular clusters structure: tidal tails

I analyse the projected density distribution of member stars over a wide area surrounding 18 Galactic globular clusters using the photometric and astrometric information provided by the second data release of the Gaia mission. A 5D mixture modelling technique has been employed to optimally isolate the signal of the cluster stellar population from the contamination of the Galactic field, taking advantage of its different distribution in the space formed by colours, magnitudes, parallaxes and proper motions. In 7 clusters I detect collimated overdensities at a >3 sigma level above the background density extending well beyond the cluster tidal radius, consistent with the distortion expected as a result of the tidal interaction with the Milky Way potential. In five of these clusters (NGC288, NGC2298, NGC5139, NGC6341 and NGC7099) spectacular tidal tails extend up to the border of the analysed field of view at 5 degrees from the centre. At large distances from the cluster centre, the orientation of the detected overdensities appears to be systematically aligned with the cluster orbital path, in agreement with the predictions of N-body simulations. The fraction of stars contained in the tidal tails of these clusters is also used to determine the first observational estimate of their present-day destruction rates

Evolution of Binary Stars in Multiple-Population Globular Clusters

The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and survival of their binary stars. Our simulations show that initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars can leave a fingerprint in the current properties of the binary population. SG binaries are disrupted more efficiently than those of the FG population resulting in a global SG binary fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We have also studied the dependence of the binary properties on the distance from the cluster centre. Although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries. The extent of the differences between the surviving FG and SG binary binding energy distribution also varies radially within the cluster and is larger in the cluster inner regions.Comment: 10 pages, 12 figures. Accepted for publication in MNRA

Evidence for a bottom-light initial mass function in massive star clusters

We have determined stellar mass functions of 120 Milky Way globular clusters and massive LMC/SMC star clusters based on a comparison of archival Hubble Space Telescope photometry with a large grid of direct N-body simulations. We find a strong correlation of the global mass function slopes of star clusters with both their internal relaxation times as well as their lifetimes. Once dynamical effects are being accounted for, the mass functions of most star clusters are compatible with an initial mass function described by a broken power-law distribution $N(m) \sim m^\alpha$ with break masses at 0.4 M$_\odot$ and 1.0 M$_\odot$ and mass function slopes of $\alpha_{Low}=-0.3$ for stars with masses $m1.0$ M$_\odot$ and $\alpha_{Med}=-1.65$ for intermediate-mass stars. Alternatively, a log-normal mass function with a characteristic mass $\log M_C = -0.36$ and width $\sigma_C=0.28$ for low-mass stars and a power-law mass function for stars with $m>1$ M$_\odot$ also fits our data. We do not find a significant environmental dependency of the initial mass function with either cluster mass, density, global velocity dispersion or metallicity. Our results lead to a larger fraction of high-mass stars in globular clusters compared to canonical Kroupa/Chabrier mass functions, increasing the efficiency of self-enrichment in clusters and helping to alleviate the mass budget problem of multiple stellar populations in globular clusters. By comparing our results with direct N-body simulations we finally find that only simulations in which most black holes are ejected by natal birth kicks correctly reproduce the observed correlations.Comment: 20 pages, 9 figures, accepted for publication in MNRAS, data is available from http://people.smp.uq.edu.au/HolgerBaumgardt/globular/index.htm

Chemical characterization of the globular cluster NGC 5634 associated to the Sagittarius dwarf spheroidal galaxy

As part of our on-going project on the homogeneous chemical characterization of multiple stellar populations in globular clusters (GCs), we studied NGC 5634, associated to the Sagittarius dwarf spheroidal galaxy, using high-resolution spectroscopy of red giant stars collected with FLAMES@VLT. We present here the radial velocity distribution of the 45 observed stars, 43 of which are member, the detailed chemical abundance of 22 species for the seven stars observed with UVES-FLAMES, and the abundance of six elements for stars observed with GIRAFFE. On our homogeneous UVES metallicity scale we derived a low metallicity [Fe/H]=-1.867 +/-0.019 +/-0.065 dex (+/-statistical +/-systematic error) with sigma=0.050 dex (7 stars). We found the normal anti-correlations between light elements (Na and O, Mg and Al), signature of multiple populations typical of massive and old GCs. We confirm the associations of NGC 5634 to the Sgr dSph, from which the cluster was lost a few Gyr ago, on the basis of its velocity and position and the abundance ratios of alpha and neutron capture elements.Comment: 16 pages, 10 figures, 11 tables; accepted for publication on Astronomy and Astrophysic