Studies of the Galactic Centre suggest that in-situ star formation may have given rise to the observed stellar population near the central supermassive black hole (SMBH). Direct evidence for a recent starburst is provided by the currently observed young stellar disc (2-7 Myr) in the central 0.5 pc of the Galaxy. This result suggests that star formation in galactic nuclei may occur close to the SMBH and produce initially flattened stellar discs. Here we explore the possible build-up and evolution of nuclear stellar clusters near SMBHs through in-situ star formation producing stellar discs similar to those observed in the Galactic Centre and other nuclei. We make use of N-body simulations to model the evolution of multiple young stellar discs and explore the potential observable signatures imprinted by such processes. Each of the five simulated discs is evolved for 100 Myr before the next one is introduced in the system. We find that populations born at different epochs show different morphologies and kinematics. Older and presumably more metal poor populations are more relaxed and extended, while younger populations show a larger amount of rotation and flattening. We conclude that star formation in central discs can reproduce the observed properties of multiple stellar populations in galactic nuclei differing in age, metallicity and kinematic properties.

The origin of the Nuclear Star Cluster in the centre of our Galaxy is still unknown. One possibility is that it formed after the disruption of stellar clusters that spiralled into the Galactic Centre due to dynamical friction. We trace the formation of the Nuclear Star Cluster around the central black hole, using state-of-the-art N-body simulations, and follow the dynamics of the neutron stars born in the clusters. We then estimate the number of Millisecond Pulsars (MSPs) that are released in the Nuclear Star Cluster, during its formation. The assembly and tidal dismemberment of globular clusters lead to a population of MSPs distributed over a radius of about 20 pc, with a peak near 3 pc. No clustering is found on the sub-parsec scale. We simulate the detectability of this population with future radio telescopes like the MeerKAT radio telescope and SKA1, and find that about of order ten MSPs can be observed over this large volume, with a paucity of MSPs within the central parsec. This helps discriminating this scenario from the in-situ formation model for the Nuclear Star Cluster that would predict an over abundance of MSPs closer to the black hole. We then discuss the potential contribution of our MSP population to the gamma-ray excess at the Galactic Centre.

The internal dynamics of globular clusters (GCs) is strongly affected by two-body interactions that bring the systems to a state of partial energy equipartition. Using a set of Monte Carlo clusters simulations, we investigate the role of the onset of energy equipartition in shaping the mass-to-light ratio (M/L) in GCs. Our simulations show that the M/L profiles cannot be considered constant and their specific shape strongly depends on the dynamical age of the clusters. Dynamically younger clusters display a central peak up to M/L $\simeq25$ $M_\odot/L_\odot$ caused by the retention of dark remnants; this peak flattens out for dynamically older clusters. Moreover, we find that also the global values of M/L correlate with the dynamical state of a cluster quantified as either the number of relaxation times a system has experienced $n_{rel}$ or the equipartition parameter $m_{eq}$: clusters closer to full equipartition (higher $n_{rel}$ or lower $m_{eq}$) display a lower M/L. We show that the decrease of M/L is primarily driven by the dynamical ejection of dark remnants, rather than by the escape of low-mass stars. The predictions of our models are in good agreement with observations of GCs in the Milky Way and M31, indicating that differences in relaxation state alone can explain variations of M/L up to a factor of $\simeq3$. Our characterization of the M/L as a function of relaxation state is of primary relevance for the application and interpretation of dynamical models.

Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment we see direct confirmation of the merging of two stellar remnant black-holes in a dynamically-formed binary, a gravitational wave source. We describe how the formation, evolution and ultimate ejection/destruction of binary systems containing black holes impacts the evolution of the cluster core. Also, through comparison with previous models of lower density, we show that the period distribution of black hole binaries formed through dynamical interactions in this high-density model favours the production of gravitational wave sources. We confirm that the number of black holes remaining in a star cluster at late times and the characteristics of the binary black hole population depend on the nature of the star cluster, critically on the number density of stars and by extension the relaxation timescale.

While tens or hundreds of stellar-remnant black holes are expected to form in globular star clusters, it is still unclear how many of those will be retained upon formation, and how many will be ejected through subsequent dynamical interactions. No such black holes have been found in any Milky Way globular cluster until the recent discovery of stellar-mass black holes in the globular cluster M22 (NGC 6656) with now an estimated population of 5-100 black holes. We present a direct N-body model of a star cluster of the same absolute and dynamical age as M22. Imposing an initial retention fraction of approx. 10% for black holes, 16 stellar-remnant black holes are retained at a cluster age of 12 Gyr, in agreement with the estimate for M22. Of those 16 BHs, two are in a binary system with a main sequence star each while also one pure black hole binary is present. We argue that multiple black holes can be present in any Milky Way cluster with an extended core radius, such as M22 or the model presented here.

Size differences of approx. 20% between red (metal-rich) and blue (metal-poor) sub-populations of globular clusters have been observed, generating an ongoing debate as to weather these originate from projection effects or the difference in metallicity. We present direct N-body simulations of metal-rich and metal-poor stellar populations evolved to study the effects of metallicity on cluster evolution. The models start with N = 100000 stars and include primordial binaries. We also take metallicity dependent stellar evolution and an external tidal field into account. We find no significant difference for the half-mass radii of those models, indicating that the clusters are structurally similar. However, utilizing observational tools to fit half-light (or effective) radii confirms that metallicity effects related to stellar evolution combined with dynamical effects such as mass segregation produce an apparent size difference of 17% on average. The metallicity effect on the overall cluster luminosity also leads to higher mass-to-light ratios for metal-rich clusters.

Direct N-body simulations of star clusters in a realistic Milky Way-like potential are carried out using the code NBODY6. Based on these simulations a new relationship between scale size and galactocentric distance is derived: the scale size of star clusters is proportional to the hyperbolic tangent of the galactocentric distance. The half-mass radius of star clusters increases systematically with galactocentric distance but levels off when star clusters orbit the galaxy beyond ~40 kpc. These simulations show that the half-mass radius of individual star clusters varies significantly as they evolve over a Hubble time, more so for clusters with shorter relaxation times, and remains constant through several relaxation times only in certain situations when expansion driven by the internal dynamics of the star cluster and the influence of the host galaxy tidal field balance each other. Indeed, the radius of a star cluster evolving within the inner 20 kpc of a realistic galactic gravitational potential is severely truncated by tidal interactions and does not remain constant over a Hubble time. Furthermore, the half-mass radius of star clusters measured with present day observations bears no memory of the original cluster size. Stellar evolution and tidal stripping are the two competing physical mechanisms that determine the present day size of globular clusters. These simulations also show that extended star clusters can form at large galactocentric distances while remaining fully bound to the host galaxy. There is thus no need to invoke accretion from an external galaxy to explain the presence of extended clusters at large galactocentric distances in a Milky Way-type galaxy.