Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical ...abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars on to the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of N-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other.
Dalessandro et al. observed a similar distribution for blue straggler stars and main-sequence turn-off stars in the Galactic globular cluster NGC 6101, and interpreted this feature as an indication ...that this cluster is not mass-segregated. Using direct N-body simulations, we find that a significant amount of mass segregation is expected for a cluster with the mass, radius and age of NGC 6101. Therefore, the absence of mass segregation cannot be explained by the argument that the cluster is not yet dynamically evolved. By varying the retention fraction of stellar-mass black holes, we show that segregation is not observable in clusters with a high black hole retention fraction (>50 per cent after supernova kicks and >50 per cent after dynamical evolution). Yet all model clusters have the same amount of mass segregation in terms of the decline of the mean mass of stars and remnants with distance to the centre. We also discuss how kinematics can be used to further constrain the presence of a stellar-mass black hole population and distinguish it from the effect of an intermediate-mass black hole. Our results imply that the kick velocities of black holes are lower than those of neutron stars. The large retention fraction during its dynamical evolution can be explained if NGC 6101 formed with a large initial radius in a Milky Way satellite.
We use radial velocities from spectra of giants obtained with the WIYN telescope, coupled with existing chemical abundance measurements of Na and O for the same stars, to probe the presence of ...kinematic differences among the multiple populations of the globular cluster (GC) M13. To characterize the kinematics of various chemical subsamples, we introduce a method using Bayesian inference along with a Markov chain Monte Carlo algorithm to fit a six-parameter kinematic model (including rotation) to these subsamples. We find that the so-called extreme population (Na-enhanced and extremely O-depleted) exhibits faster rotation around the centre of the cluster than the other cluster stars, in particular, when compared with the dominant 'intermediate' population (moderately Na-enhanced and O-depleted). The most likely difference between the rotational amplitude of this extreme population and that of the intermediate population is found to be ~4 km s super( -1), with a 98.4 per cent probability that the rotational amplitude of the extreme population is larger than that of the intermediate population. We argue that the observed difference in rotational amplitudes, obtained when splitting subsamples according to their chemistry, is not a product of the long-term dynamical evolution of the cluster, but more likely a surviving feature imprinted early in the formation history of this GC and its multiple populations. We also find an agreement (within uncertainties) in the inferred position angle of the rotation axis of the different subpopulations considered. We discuss the constraints that these results may place on various formation scenarios.
The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analog of large star-formation events in the distant universe. We determined the recent formation history and the initial ...mass function (IMF) of massive stars in 30 Doradus on the basis of spectroscopic observations of 247 stars more massive than 15 solar masses (Formula: see text). The main episode of massive star formation began about 8 million years (My) ago, and the star-formation rate seems to have declined in the last 1 My. The IMF is densely sampled up to 200 Formula: see text and contains 32 ± 12% more stars above 30 Formula: see text than predicted by a standard Salpeter IMF. In the mass range of 15 to 200 Formula: see text, the IMF power-law exponent is Formula: see text, shallower than the Salpeter value of 2.35.
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/overestimated up to 50 per cent 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 centre 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.
ABSTRACT
We report on the detection of a black hole (NGC 1850 BH1) in the ∼100-Myr-old massive cluster NGC 1850 in the Large Magellanic Cloud. It is in a binary system with a main-sequence turn-off ...star (4.9 ± 0.4 M⊙), which is starting to fill its Roche lobe and is becoming distorted. Using 17 epochs of Very Large Telescope/Multi-Unit Spectroscopic Explorer observations, we detected radial velocity variations exceeding 300 km s−1 associated with the target star, linked to the ellipsoidal variations measured by the fourth phase of the Optical Gravitational Lensing Experiment in the optical bands. Under the assumption of a semidetached system, the simultaneous modelling of radial velocity and light curves constrains the orbital inclination of the binary to 38° ± 2°, resulting in a true mass of the unseen companion of $11.1_{-2.4}^{+2.1}\,{\rm M}_{\odot }$. This represents the first direct dynamical detection of a black hole in a young massive cluster, opening up the possibility of studying the initial mass function and the early dynamical evolution of such compact objects in high-density environments.
Abstract
Using data from Gaia DR2, we study the radial number density profiles of the Galactic globular cluster sample. Proper motions are used for accurate membership selection, especially crucial ...in the cluster outskirts. Due to the severe crowding in the centres, the Gaia data are supplemented by literature data from HST and surface brightness measurements, where available. This results in 81 clusters with a complete density profile covering the full tidal radius (and beyond) for each cluster. We model the density profiles using a set of single-mass models ranging from King and Wilson models to generalized lowered isothermal limepy models and the recently introduced spes models, which allow for the inclusion of potential escapers. We find that both King and Wilson models are too simple to fully reproduce the density profiles, with King (Wilson) models on average underestimating (overestimating) the radial extent of the clusters. The truncation radii derived from the limepy models are similar to estimates for the Jacobi radii based on the cluster masses and their orbits. We show clear correlations between structural and environmental parameters, as a function of Galactocentric radius and integrated luminosity. Notably, the recovered fraction of potential escapers correlates with cluster pericentre radius, luminosity, and cluster concentration. The ratio of half mass over Jacobi radius also correlates with both truncation parameter and PE fraction, showing the effect of Roche lobe filling.
ABSTRACT
We have determined stellar mass functions of 120 Milky Way globular clusters and massive Large Magellanic Cloud/Small Magellanic Cloud 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 and 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) ∼ mα with break masses at 0.4 and 1.0 M⊙ and mass function slopes of αLow = −0.3 for stars with masses m < 0.4 M⊙, αHigh = −2.30 for stars with m > 1.0 M⊙, and αMed = −1.65 for intermediate-mass stars. Alternatively, a lognormal mass function with a characteristic mass log MC = −0.36 and width σC = 0.28 for low-mass stars and a power-law mass function for stars with m > 1 M⊙ also fit our data. We do not find a significant environmental dependence of the initial mass function on 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.
ABSTRACT
The distribution of stars and stellar remnants (white dwarfs, neutron stars, and black holes) within globular clusters holds clues about their formation and long-term evolution, with ...important implications for their initial mass function (IMF) and the formation of black hole mergers. In this work, we present best-fitting multimass models for 37 Milky Way globular clusters, which were inferred from various data sets, including proper motions from Gaia EDR3 and HST, line-of-sight velocities from ground-based spectroscopy and deep stellar mass functions from HST. We use metallicity-dependent stellar evolution recipes to obtain present-day mass functions of stars and remnants from the IMF. By dynamically probing the present-day mass function of all objects in a cluster, including the mass distribution of remnants, these models allow us to explore in detail the stellar (initial) mass functions of a large sample of Milky Way GCs. We show that, while the low-mass mass function slopes are strongly dependent on the dynamical age of the clusters, the high-mass slope (α3; m > 1 M⊙) is not, indicating that the mass function in this regime has generally been less affected by dynamical mass loss. Examination of this high-mass mass function slope suggests an IMF in this mass regime consistent with a Salpeter IMF is required to reproduce the observations. This high-mass IMF is incompatible with a top-heavy IMF, as has been proposed recently. Finally, based on multimass model fits to our sample of Milky Way GCs, no significant correlation is found between the high-mass IMF slope and cluster metallicity.