Antonio Sollima

We present the result of a sample of B-stars in the Large Magellanic Cloud young double stellar cluster NGC 1850 A and NGC 1850 B, observed with the integral-field spectrograph at the Very Large Telescope, the Multi Unit Spectroscopic Explorer. We compare the observed equivalent widths (EWs) of four He lines (4922 Å 5015 Å 6678 Å and 7065 Å) with the ones determined from synthetic spectra computed with different He mass fraction (Y=0.25, 0.27, 0.30 and 0.35) with the code SYNSPEC, that takes into account the non-LTE effect. From this comparison, we determined the He mass fraction of the B stars, finding a not homogeneous distribution. The stars can be divided in three groups, He-weak (Y< 0.24) and the He-normal (0.24≤Y≤ 0.26) belonging to the MS of NGC 1850 A, and the He-rich stars 0.33≤Y≤0.38) situated in the MS associated to NGC 1850 B. We have analyzed the stellar rotation as possible responsible of the anomalous features of the He lines in the He-rich stars. We provide a simple analysis of the differences between the observed EWs and the ones obtained from the theoretical models with different rotation velocity (Vsin i = 0 and 250 Km/s). The resolution of the MUSE spectra do not allow to get a conclusive result, however our analysis support the He-enhanced hypothesis.

We developed self-consistent dynamical models of stellar systems in the framework of quasi-linear modified Newtonian dynamics (QUMOND). The models are constructed from the anisotropic distribution function of Gunn and Griffin, combined with the modified Poisson equation defining this gravitation theory and take into account the external field effect. We have used these models, and their Newtonian analogues, to fit the projected density and the velocity dispersion profiles of a sample of 18 Galactic globular clusters, using the most updated datasets of radial velocities and Gaia proper motions. We have thus obtained, for each cluster, estimates of the dynamical mass-to-light ratio (M/L) for each theory of gravity. The selected clusters have accurate proper motions and a well sampled mass function down to the very low mass regime. This allows us to constrain the degree of anisotropy and to provide, from comparison with stellar evolution isochrones, a dynamics-independent estimate of the minimum mass-to-light ratio (M/L)_min. Comparing the best-fitting dynamical M/L with (M/L)_min, we find that for none of the analyzed clusters the two gravity theories are significantly incompatible with the observational data, although for one of them (NGC 5024) the dynamical M/L predicted by QUMOND lies at 2.8σ below (M/L)_min. Though the proposed approach suffers from some limitations (in particular the lack of a treatment of mass segregation), the obtained results suggest that the kinematics of globular clusters in a relatively weak external field can be a powerful tool to prove alternative theories of gravitation.

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) ∼ m^α with break masses at 0.4 M_⊙ 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 log-normal mass function with a characteristic mass log M_C = −0.36 and width σ_C = 0.28 for low-mass stars and a power-law mass function for stars with m > 1 M_⊙ 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.

We present the results of a spectroscopic survey of the outskirts of 4 globular clusters—NGC 1261, NGC 4590, NGC 1904, and NGC 1851— covering targets within 1 degree from the cluster centres, with 2dF/AAOmega on the Anglo-Australian Telescope (AAT) and FLAMES on the Very Large Telescope (VLT). We extracted chemo-dynamical information for individual stars, from which we estimated the velocity dispersion profile and the rotation of each cluster. The observations are compared to direct N-body simulations and appropriate limepy/spes models for each cluster to interpret the results. In NGC 1851, the detected internal rotation agrees with existing literature, and NGC 1261 shows some rotation signal beyond the truncation radius, likely coming from the escaped stars. We find that the dispersion profiles for both the observations and the simulations for NGC 1261, NGC 1851, and NGC 1904 do not decrease as the limepy/spes models predict beyond the truncation radius, where the N-body simulations show that escaped stars dominate; the dispersion profile of NGC 4590 follows the predictions of the limepy/spes models, though the data do not effectively extend beyond the truncation radius. The increasing/flat dispersion profiles in the outskirts of NGC 1261, NGC 1851 and NGC 1904, are reproduced by the simulations. Hence, the increasing/flat dispersion profiles of the clusters in question can be explained by the tidal interaction with the Galaxy without introducing dark matter.

We determine the mass functions (MFs) and the dynamical parameters of 15 nearby open clusters (OCs), using the unprecedented data set of the Gaia Early Data Release 3. We select the members of each cluster by combining the photometric (colour and magnitude) and astrometric (parallax and proper motions) parameters of stars, minimizing the contamination from Galactic field interlopers. By comparing the observed distribution of stars along the cluster main sequence with the best-fitting synthetic population, we find the present-day MF and the binary fraction of the OCs, along with their dynamical parameters like mass, half-mass radius, and half-mass relaxation time. We found that the global present-day MF of OCs are consistent with a single power-law function, F(m) ∝ m^α, with slopes −3 < α < −0.6 including both subsolar, 0.2 < m/M_⊙ < 1, and supersolar mass regimes. A significant correlation between the MF slope and the ratio of age to half-mass relaxation time is evidenced, similarly to the same conclusion already observed among Galactic globular clusters. However, OCs evolve along different tracks in comparison with the globular clusters, possibly indicating primordial differences in their initial mass function (IMF). The comparison with Monte Carlo simulations suggests that all the analysed OCs could be born with an IMF with slope α_IMF < −2.3. We also show that the less evolved OCs have an MF consistent with that of the solar neighbourhood, indicating a possible connection between the dissolution of OCs and the formation of the Galactic disc.

Astronomy & Astrophysics (A&A) is an international journal which publishes papers on all aspects of astronomy and astrophysics

We present the result of a survey of Monte Carlo simulations of globular clusters hosting two generations of stars including a large (⁠|$f_{b}=50{{\ \rm per\ cent}}$|⁠) fraction of primordial binaries in both populations. The dynamical evolution of the two stellar populations is followed for a Hubble time taking into account the effect of the tidal field, two-body relaxation, stellar evolution and three/four-bodies interactions. The fraction of surviving binaries, once accounted for the observational bias and uncertainties, is compared with the available radial velocity time-series performed in real globular clusters, and it is used to constrain the initial spatial concentration of the second generation. The fraction of second generation binaries appears to depend only on the ratio between the total cluster mass and the initial size of the second generation which determines the average velocity dispersion across the extent of this stellar population. In spite of the various uncertainties, we find that the observed fraction can be obtained only assuming a strong initial concentration of the second generation (r_{h, S} ∼ 0.1 (M/10⁶ M_⊙) pc). The evolution of the first generation binary fraction is more sensitive to the tidal field strength (with a non negligible effect of the cluster orbital eccentricity) since the tidal field has a direct impact on the first generation structural properties.

We present observations of the stellar kinematics of the centre of the core collapsed globular cluster M15 obtained with the MUSE integral field spectrograph on the Very Large Telescope operating in narrow field mode. Thanks to the use of adaptive optics, we obtain a spatial resolution of 0.1 arcsec and are able to reliably measure the radial velocities of 864 stars within 8 arcsec of the centre of M15, thus providing the largest sample of radial velocities ever obtained for the innermost regions of this system. Combined with previous observations of M15 using MUSE in wide field mode and literature data, we find that the central kinematics of M15 are complex with the rotation axis of the core of M15 offset from the rotation axis of the bulk of the cluster. While this complexity has been suggested by previous work, we confirm it at higher significance and in more detail.

I simulate the evolution of a stellar system hosting two stellar populations whose initial set-up is defined according to the two main scenarios proposed for the origin of multiple populations in Galactic globular clusters: (i) formation of a second generation from a cooling flow of pristine+polluted gas and (ii) accretion of polluted gas on to the proto-stellar discs of a fraction of low-mass stars. For this purpose, Monte Carlo simulations containing from 105 up to 3 × 106 particles have been run including the effect of stellar evolution, binary interactions, external tidal field, and a detailed modelling of the proto-stellar disc structure. The early accretion of gas on to proto-stellar discs is unable to produce discrete populations and to alter the chemical composition of a significant ($\gt 10{{\ \rm per\ cent}}$) fraction of stars unless a disc lifetime larger (tdisc ∼ 20 Myr) than that predicted by models is assumed. Moreover, in this scenario the mixing time-scale of the two populations is too short to reproduce the observed segregation of the chemically enriched population. On the other hand, simulations run within the cooling flow scenario can evolve after a Hubble time into stellar systems with a first-to-second population mass ratio similar to that observed in globular clusters, provided that an initial filling-factor rh/rJ > 0.15 is adopted. However, in the weak tidal field regime a radial segregation of the second population stronger than what observed in Milky Way globular clusters at large Galactocentric distances is predicted. This discrepancy disappears in simulations following eccentric orbits in a realistic axisymmetric potential.

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σ 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 (NGC 288, NGC 2298, NGC 5139, NGC 6341, and NGC 7099) spectacular tidal tails extend up to the border of the analysed field of view at 5 deg 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.