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.

Once the age and metallicity are fixed, the colour distribution of horizontal branch stars in a globular cluster depends on few parameters: the helium abundance of the population and the mass lost during the pre-HB stages. These parameters are usually derived from the HB itself, hence they are degenerate. Breaking this degeneracy and understanding their role is a tricky and challenging problem that no study has solved yet. Combining the information obtained from the chromosome maps and the analysis of multi-band photometry with state of the art stellar evolution models, we can obtain a solid estimate of Y for the various stellar populations in a GC. We will then have, for the first time, the possibility to break the parameters’ degeneracy on the HB, understand the role of the mass loss, and lay the foundation to build another piece of the multiple populations mosaic.

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 collected radial velocities of more than 50.000 individual stars in 156 Galactic globular clusters (GGC) and matched them with HST photometry and Gaia DR2 proper motions. This allowed us to derive the GGC’s mean proper motions and space velocities. By fitting a large set of N-body simulations to their velocity dispersion and surface density profiles, combined with new measurements of their internal radially dependent mass functions, we have determined their present-day masses and structural parameters, and for 144 GGCs their internal kinematics. We also derive the initial cluster masses by calculating the cluster orbits backwards in time applying suitable recipes to account for mass-loss and dynamical friction. The new fundamental parameters of GGCs are publicly available via an online database, which will regularly be updated.

Using direct N-body simulations of self-gravitating systems we study the dependence of dynamical chaos on the system size N. We find that the N-body chaos quantified in terms of the largest Lyapunov exponent Λmax decreases with N. The values of its inverse (the so-called Lyapunov time tλ) are found to be smaller than the two-body collisional relaxation time but larger than the typical violent relaxation time, thus suggesting the existence of another collective time scale connected to many-body chaos.

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.

Supermassive black holes are found in most galactic nuclei. A large fraction of these nuclei also contain a nuclear stellar cluster surrounding the black hole. Here we consider the idea that the nuclear stellar cluster formed first and that the supermassive black hole grew later. In particular we consider the merger of three stellar clusters to form a nuclear stellar cluster, where some of these clusters contain a single intermediate-mass black hole (IMBH). In the cases where multiple clusters contain IMBHs, we discuss whether the black holes are likely to merge and whether such mergers are likely to result in the ejection of the merged black hole from the nuclear stellar cluster. In some cases, no supermassive black hole will form as any merger product is not retained. This is a natural pathway to explain those galactic nuclei that contain a nuclear stellar cluster but apparently lack a supermassive black hole; M33 being a nearby example. Alternatively, if an IMBH merger product is retained within the nuclear stellar cluster, it may subsequently grow, e.g. via the tidal disruption of stars, to form a supermassive black hole.

We investigate, for the first time, the formation and evolution of the tidal tail released from a young Pleiades-like star cluster due to expulsion of primordial gas in a realistic gravitational field of the Galaxy. The tidal tails (as well as clusters) are integrated by nbody6 from their embedded phase for more than 300 Myr. We vary the star formation efficiency (SFE) from 33% to 100% and the timescales of gas expulsion as free parameters, and provide predictions for the morphology and kinematics of the evolved tail for each of the models. The resulting tail properties are intended for comparison with Gaia measurements, where an inverse analysis of our findings might constrain some of the poorly understood conditions and processes in embedded star clusters during the gas phase and gas expulsion.

We present preliminary results of the wide-field photometric study of the isolated elliptical galaxy NGC 1172, and its globular cluster system. Our data was obtained with the GMOS camera mounted on the Gemini South telescope, in the g′, r′, i′ and z′ bands. The aim of this work is to further our understanding of the evolution of NGC 1172, and to look for possible explanations for its unusual high specific frequency.

We present an analysis of the globular cluster system (GCS) of the galaxy NGC 3613, an intrinsically bright elliptical galaxy (MV = −21.5) in a low density environment (it is the central galaxy of a group of a dozen galaxies). Based on Gemini/GMOS photometry of NGC 3613 we obtained the following properties for this GCS. A ‘blue tilt’ is detected in the colour-magnitude diagram. The colour distribution is bimodal, presenting the two classical globular cluster (GC) sub-populations. The spatial and azimuthal projected distributions show that red sub-population correlates with the stellar component of the host galaxy.

We investigate the dynamics of supermassive black holes (SMBHs) in galactic cores by means of a semi-analytic model based on the Langevin equation, including dynamical friction and stochastic noise accounting for the gravitational interactions with stars. The model is validated against direct N-body simulations of intermediate-mass black holes in stellar clusters where a realistic number of particles is accessible. For the galactic case, we find that the SMBH experiences a Brownian-like motion with a typical displacement from the geometric center of the Galaxy of a few parsecs, for system parameters compatible with M87.

The study of ages, helium mass fraction (Y) and chemical composition of globular clusters in dwarf galaxies is important for understanding the physical conditions at the main evolutionary stages of the host galaxies and for constraining the build-up histories of large galaxies. We present the analysis of integrated-light spectra of 8 extragalactic and 20 Galactic globular clusters (GCs) using our population synthesis method. We calculate synthetic spectra of GCs according to the defined stellar mass functions using model atmospheres and stellar parameters ([Fe/H], Teff, and logg) set by theoretical isochrones. The main advantage of our method is the ability to determine not only chemical composition but also the age and mean Y in a cluster by modelling and analysis of Balmer absorption lines. The knowledge of Y and anomalies of light elements in star clusters is one of the key points for understanding the phenomenon of multiple stellar populations.

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.

Over the last decade, much of the key questions in Galactic Archaeology have been asnwered by studying the Milky Way’s globular cluster (GC) system. Following on this, it has been shown that a substantial fraction of the Milky Way’s stellar halo field arises from GC dissolution. In this work, we make use of the latest data release fromn the APOGEE survey to study GC dissolution ratios in different spatial regions of the Galaxy. Our results will allow us to constrain many astrophysical questions, such as: the origin of N-Rich stars, the mass contribution from GCs to the stellar halo of the Galaxy, the origin of the Galactic GC system and the mass assembly of the Milky Way.

We present the results obtained from Near-UV observations of the cluster NGC 5466 taken with UVIT onboard AstroSat to study the radial distribution of Blue Straggler Stars (BSSs), covering the cluster region up to a radius of 14'. Our study confirms that the BSSs are more centrally concentrated than Horizontal Branch (HB) stars in the cluster. We do not find a statistically significant rising trend in the radial distribution of BSSs for the outer regions, as most of the previously catalogued BSSs that are located in the cluster outskirts were found to be non-members, quasars or galaxies. This study highlights the importance of UV imaging combined with membership information to probe the radial distribution of BSSs.

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.

For the first time, we report the identification of NUV bright red clump (RC) stars and the extension of RC stars over two magnitudes both in color and magnitude axis in NUV vs (NUV – optical) color magnitude diagram. We find that the extension of RC is not due to photometric uncertainties. We suggest that the extension could be an effect of field star contamination. We also suggest that if it is an intrinsic property of the cluster then age and/or metallicity spread within the cluster could be the possible reasons for extended RC.

With the hundreds of merging binary black holes (BHs) expected to be detected by LIGO, LISA, and other upcoming instruments, the modelling of astrophysical channels that lead to the formation of compact BH binaries has become of crucial importance. BHs of any size can form bound systems in every astrophysical environment, from the field to galactic nuclei. If a binary is too wide, it needs a catalysis process to harden and merge, as in the case a third objects orbiting the BH binary on a distant orbit. In this case, Kozai-Lidov cycles can pump up the binary eccentricity, thus driving it to a merger thanks to efficient energy dissipation at the pericenter. Some remarkable scenarios where the Kozai-Lidov mechanism operates are in triple and quadruple systems of stellar BHs, and in intermediate-mass BH-stellar BH binaries in orbit around a central supermassive BH in galactic nuclei.

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.