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.

Present-day structural and kinematical properties of multiple populations (MPs) can provide useful information about the physical mechanisms driving the formation and early evolution of globular clusters (GCs). As part of a large project aimed at characterizing the kinematics of MPs, here we present a detailed multi-epoch analysis of the low-mass GC NGC6362. We find that MPs in this system show significant differences in their line-of-sight velocity dispersion profiles. This result is totally unexpected given the dynamical age and fraction of mass lost by NGC6362. We also find that the binary fraction is remarkably larger in the first (FP) than in the second population (SP). We show that such a difference can efficiently inflate the velocity dispersion of FP at intermediate/large cluster-centric distances with respect to SP. Indeed, our results nicely match the predictions of state-of-the art N-body simulations of the co-evolution of MPs in GCs including the effect of binaries.

N-body simulation is the necessary tool to investigate the evolution of star clusters. It is important to develop a time integration method which guarantees the appropriate accuracy and calculation cost.

Here we present a new simple method for long term simulation of stars around a massive black hole in stellar systems. Usually the time integration orbits of stars revolving a massive black hole requires much simulation time. We introduce a time transformation which is a kind of ”inverse KS regularization” of time. Using our method, the integration of the long term evolution near a black hole (BH) becomes easier, especially applied to relatively large star clusters and the Galactic Center.