I present the first results of multiple stellar populations in globular clusters from James Webb Space Telescope (JWST) data. We obtained combined NIRSpec and NIRCam data of the GC 47 Tucanae to investigate the properties of the different stellar populations in the low-mass regime, down to ∼0.5 M⊙. Our analysis of both the photometric data and the spectra suggests that the multiple stellar populations among M dwarfs share similar properties of those extensively studied on bright red giant branch stars. The comparison between stellar populations’ properties at different stellar mass, e.g. at the red giants and M dwarfs’ phases, can potentially provide strong constraints to the formation scenarios of the still eluding phenomenon of the multiple populations in star clusters.

I present the first evidence of multiple populations in the globular cluster (GCs) 47 Tucanae based on images collected with the near-infrared camera (NIRCam) on board the James Webb Space Telescope (JWST). While NIRCam photometry is poorly sensitive to multiple populations among stars brighter than the main-sequence (MS) knee, the M-dwarfs more-massive than ∼0.1 define a wide color range due to multiple populations. The star-to-star color differences are mostly due to the different amounts of water vapor (hence oxygen) that affect the spectra of M-dwarfs. The chromosome map unveils an extended first population (1P) composed of M-dwarfs with different metallicities and three main groups of second-population (2P) stars that are depleted in oxygen with respect to the 1P. I present the discovery of an MS of very-low-mass stars and tentatively associated it with a sequence composed of O-rich stars alone.

We explore the relationship between globular cluster total number, NGC, and central black hole mass, M•, in spiral galaxies. Including cosmic scatter, log M• ∝ (1.64 ± 0.24) log NGC. Whereas in ellipticals the correlation is linear [log M• ∝ (1.02 ± 0.10) log NGC], and hence could be due to statistical convergence through mergers, this mechanism cannot explain the much steeper correlation in spirals. Additionally, we derive total stellar galaxy mass, M*, from its two-slope correlation with NGC (Hudson et al. 2014). In the M• versus M* parameter space, with M* derived from NGC, M• ∝ (1.48 ± 0.18) log M* for ellipticals, and M• ∝ (1.21 ± 0.16) log M* for spirals. The observed agreement between ellipticals and spirals may imply that black holes and galaxies co-evolve through “calm” accretion, AGN feedback and other secular processes.

We have used Hubble Space Telescope and ground-based photometry to determine total V-band magnitudes and mass-to-light ratios of more than 150 Galactic globular clusters. We do this by summing up the magnitudes of their individual member stars, using colour-magnitude information, Gaia DR2 proper motions, and radial velocities to distinguish cluster stars from background stars. Our new magnitudes confirm literature estimates for bright clusters with $V<8$, but can deviate by up to two magnitudes from literature values for fainter clusters. They lead to absolute mass-to-light ratios that are confined to the narrow range $1.4<M/L_V<2.5$, significantly smaller than what was found before. We also find a correlation between a cluster’s $M/L_V$ value and its age, in agreement with theoretical predictions. The $M/L_V$ ratios of globular clusters are also in good agreement with those predicted by stellar isochrones, arguing against a significant amount of dark matter inside globular clusters. We finally find that, in agreement with what has been seen in M 31, the magnitude distribution of outer halo globular clusters has a tail towards faint clusters that is absent in the inner parts of the Milky Way.

To look at propagating winds from evolved stars into the interstellar medium is to look at how they are sustained. To understand their origins, we must look to the circumstances that create them in the first instance. In this article, I examine the physical conditions under which pulsation-enhanced, dust-driven winds are first generated. These initial conditions can help constrain the late evolutionary stages of these stars and provide insight into the mechanisms that cause the mass loss itself.

The opportunities offered by JWST and the ELT for the detection and study of forming/just formed globular clusters at high redshifts are illustrated, also alluding at the unique insight we may get on the very early stages of galaxy formation.

Recent observations of globular clusters encourage to revise some aspects of the traditional paradigm, in which they were considered to be isotropic in velocity space and non-rotating. However, the theory of collisionless spheroids with some kinematic richness has seldom been studied. We present here a further step in this direction, owing to new results regarding the linear stability of rotating Plummer spheres, with varying anisotropy in velocity space and total amount of angular momentum. We extend the well-known radial orbit instability to rotating systems, and discover a new regime of instability in fast rotating, tangentially anisotropic systems.

The blooming era of precision astrometry for Galactic studies truly brings the rich internal dynamics of globular clusters to the centre stage. But several aspects of our current understanding of fundamental collisional stellar dynamics cannot match such new-generation data and the theoretical ambitions they trigger. This rapidly evolving context offers the stimulus to address a number of old and new questions concerning the phase space properties of this class of stellar systems.

The Asymptotic Giant Branch (AGB) scenario for the formation of multiple populations has been quantitatively studied in the course of the last twenty years, examining the aspects both of nucleosynthesis and of the dynamics of formation of new stars in a cooling flow at the center of the first generation cluster, and of the following N–body evolution. The large complexity of these studies finds many validations in the properties of multiple populations. Here I shortly summarize recent accomplishments in the study of the evolution of massive AGBs and super–AGBs including the explanation of anomalous high lithium abundances in ‘extreme’ second generation stars in ω Cen and NGC 2808.

Recent observations of globular clusters imposed major revisions to the previous paradigm, in which they were considered to be isotropic in velocity space and non-rotating. However, the theory of collisionless spheroids with some kinematic richness has seldom been studied. We present here a first step in this direction, owing to new results regarding the linear stability of rotating Plummer spheres, with varying anisotropy in velocity space and total amount of angular momentum. We extend the well-known radial orbit instability to rotating systems, and discover a new regime of instability in fast rotating, tangentially anisotropic systems.

Stars in globular clusters lose mass through slow stellar winds that are retained by the stellar cluster and should contribute to build up a non-negligible intracluster medium over time. However, all the observations so far found only a negligible amount of gas in GCs. We propose here to test different mechanisms such as ram-pressure stripping by the motion of the GC in the Galactic halo medium and the inclusion of ionising sources to explain the lack of gas in GCs. We use full 3D hydrodynamical simulations taking into account stellar winds, ionising radiation, radiative heating and radiative pressure. We find that the combined effect of ram-pressure and ionisation are able to explain the negligible amount of gas observed in the core of intermediate-mass and massive GCs.

Over a hundred millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic N-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses ≍ 2 × 105M⊙ can produce up to 10 – 20 MSPs, while a very massive GC model with mass ≍ 106M⊙ can produce close to 100. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations.

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, seems to be related to the initial configuration of the system. In recent work (Milone et al. 2018), we analyzed for the first time the kinematics of the different stellar populations in NGC 0104 (47 Tucanae) over a large field of view, exploiting the Gaia Data Release 2 proper motions combined with multi-band ground-based photometry. In this paper, based on the work by Cordoni et al. (2019), we extend this analysis to six GCs, namely NGC 0288, NGC 5904 (M 5), NGC 6121 (M 4), NGC 6752, NGC 6838 (M 71) and further explore NGC 0104. Among the analyzed clusters only NGC 0104 and NGC 5904 show significant rotation on the plane of the sky. Interestingly, multiple stellar populations in NGC 5904 exhibit different rotation curves.

Setting the formation of globular clusters (GCs) within a cosmological context and characterising the properties of proto-GCs at high redshift is currently a major challenge. In this work, we address that challenge by exploring a suit of high-resolution cosmological simulations from the First Billion Years (FiBY) project z at ⩾6 to investigate theoretical scenarios concerning the formation of old, low-mass stellar systems with a particular focus on GCs. Two distinct groups of objects are identified in the simulations. The first group of objects, with a high baryon fraction, we associate with proto-GCs. The second group, that exhibit a high stellar fraction, could be forming ultra-faint dwarf galaxies (UFDs). The objects with high baryon fraction are promising proto-GC candidates because they have little to no dark matter (DM), have number densities consistent with predictions from the literature, are very compact and have a high stellar density. We fit and also assess the redshift-zero globular system mass - halo mass relation and find it provides a reasonable fit to our proto-GC objects, indicating that this relation is likely set at formation.

The Galactic globular cluster system went and is still going through dynamical processes that require to be explored in detail. Here we illustrate how primordial massive globular clusters born in the Milky Way’s disc evolved by stripping material from each other or even merging very early during their lives. These processes might explain the puzzling presence of star-by-star spreads in iron content observed in massive globular clusters and should be taken into account when studying globular cluster stellar populations. In this context, we show how the direct comparison between the predictions provided by our direct N-body simulations and observations can shed light on the origin and chemo-dynamical evolution of globular clusters.

We investigate the dissolution process of star clusters embedded in an external tidal field and harboring a subsystem of stellar-mass black hole. For this purpose we analyzed the MOCCA models of real star clusters contained in the Mocca Survey Database I. We showed that the presence of a stellar-mass black hole subsystem in tidally filling star cluster can lead to abrupt cluster dissolution connected with the loss of cluster dynamical equilibrium. Such cluster dissolution can be regarded as a third type of cluster dissolution mechanism. We additionally argue that such a mechanism should also work for tidally under-filling clusters with a top-heavy initial mass function.

I present the results of a survey of the kinematics of a large sample of Galactic globular clusters performed thanks to the synergy between the 2nd Gaia data release and the most extensive collection of radial velocities. This unprecedented dataset of 3D velocities of thousand of stars in 62 globular clusters has been used to investigate the rotation patterns of these stellar systems providing insight into the impact of two-body relaxation and tides on the formation and evolution of their rotation.

We review the implications of the Gaia Data Release 2 catalogue for studying the dynamics of Milky Way globular clusters, focusing on two separate topics.

The first one is the analysis of the full 6-dimensional phase-space distribution of the entire population of Milky Way globular clusters: their mean proper motions (PM) can be measured with an exquisite precision (down to 0.05 mas yr−1, including systematic errors). Using these data, and a suitable ansatz for the steady-state distribution function (DF) of the cluster population, we then determine simultaneously the best-fit parameters of this DF and the total Milky Way potential. We also discuss possible correlated structures in the space of integrals of motion.

The second topic addresses the internal dynamics of a few dozen of the closest and richest globular clusters, again using the Gaia PM to measure the velocity dispersion and internal rotation, with a proper treatment of spatially correlated systematic errors. Clear rotation signatures are detected in 10 clusters, and a few more show weaker signatures at a level ∼0.05 mas yr−1. PM dispersion profiles can be reliably measured down to 0.1 mas yr−1, and agree well with the line-of-sight velocity dispersion profiles from the literature.

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.

The recent measurements of internal variations of helium in Galactic and extragalactic Globular Clusters (GCs) set binding constraints to the models of formation of Multiple Populations (MPs) in GCs, and gave rise, at the same time, to crucial questions related with the influence of the environment on MP formation as well as with the role played by GCs in the early galactic formation. We present the most recent estimates of helium enrichment in the main populations of a large sample of Galactic and extragalactic GCs.