The origin of multiple stellar populations in Globular Clusters (GCs) is one of the greatest mysteries of modern stellar astrophysics. N-body simulations suggest that the present-day dynamics of GC stars can constrain the events that occurred at high redshift and led to the formation of multiple populations. Here, we combine multi-band photometry from the Hubble Space Telescope (HST) and ground-based facilities with HST and Gaia Data Release 2 proper motions to investigate the spatial distributions and the motions in the plane of the sky of multiple populations in the type II GCs NGC 5139 ($\omega\,$Centauri) and NGC 6656 (M 22). We first analyzed stellar populations with different metallicities. Fe-poor and Fe-rich stars in M 22 share similar spatial distributions and rotation patterns and exhibit similar isotropic motions. Similarly, the two main populations with different iron abundance in $\omega\,$Centauri share similar ellipticities and rotation patterns. When analyzing different radial regions, we find that the rotation amplitude decreases from the center towards the external regions. Fe-poor and Fe-rich stars of $\omega\,$Centauri are radially anisotropic in the central region and show similar degrees of anisotropy. We also investigate the stellar populations with different light-element abundances and find that their N-rich stars exhibit higher ellipticity than N-poor stars. In $\omega\,$Centauri Centauri both stellar groups are radially anisotropic. Interestingly, N-rich, Fe-rich stars exhibit different rotation patterns than N-poor stars with similar metallicities. The stellar populations with different nitrogen of M 22 exhibit similar rotation patterns and isotropic motions. We discuss these findings in the context of the formation of multiple populations.

We present a detailed study of the composition of 20 M giants in the Galactic center with 15 of them confirmed to be in the Nuclear Star Cluster. As a control sample we have also observed 7 M giants in the Milky Way Disk with similar stellar parameters. All 27 stars are observed using the NIRSPEC spectograph on the KECK II telescope in the K-band at a resolving power of R=23,000. We report the first silicon abundance trends versus [Fe/H] for stars in the Galactic center. While finding a disk/bulge like trend at subsolar metallicities, we find that [Si/Fe] is enhanced at supersolar metallicities. We speculate on possible enrichment scenarios to explain such a trend. However, the sample size is modest and the result needs to be confirmed by additional measurements of silicon and other \textalpha-elements. We also derive a new distribution of [Fe/H] and find the most metal rich stars at [Fe/H]=+0.5 dex, confirming our earlier conclusions that the Galactic center hosts no stars with extreme chemical composition.

The Sagittarius dwarf spheroidal galaxy (Sgr dSph) is in an advanced stage of disruption but still hosts its nuclear star cluster (NSC), M54, at its center. In this paper, we present a detailed kinematic characterization of the three stellar populations present in M54: young metal-rich (YMR); intermediate-age metal-rich (IMR); and old metal-poor (OMP), based on the spectra of $\sim6500$ individual M54 member stars extracted from a large MUSE/VLT dataset. We find that the OMP population is slightly flattened with a low amount of rotation ($\sim0.8$ km s$^{-1}$) and with a velocity dispersion that follows a Plummer profile. The YMR population displays a high amount of rotation ($\sim5$ km s$^{-1}$) and a high degree of flattening, with a lower and flat velocity dispersion profile. The IMR population shows a high but flat velocity dispersion profile, with some degree of rotation ($\sim2$ km s$^{-1}$). We complement our MUSE data with information from \textit{Gaia DR2} and confirm that the stars from the OMP and YMR populations are comoving in 3D space, suggesting that they are dynamically bound. While dynamical evolutionary effects (e.g. energy equipartition) are able to explain the differences in velocity dispersion between the stellar populations, the strong differences in rotation indicate different formation paths for the populations, as supported by an $N$-body simulation tailored to emulate the YMR-OMP system. This study provides additional evidence for the M54 formation scenario proposed in our previous work, where this NSC formed via GC accretion (OMP) and in situ formation from gas accretion in a rotationally supported disc (YMR).

We present multi-band photometry covering $\sim$ 5\deg $\times$ 5\deg across $\omega$ Cen collected with the Dark Energy Camera, combined to Hubble Space Telescope and Wide Field Imager data for the central regions. The unprecedented photometric accuracy and field coverage allowed us to confirm the different spatial distribution of blue and red main-sequence stars, and of red-giant branch (RGB) stars with different metallicities. The ratio of the number of blue to red main-sequence stars shows that the blue main-sequence sub-population has a more extended spatial distribution compared to the red main-sequence one, and the frequency of blue main-sequence stars increases at a distance of $\sim$ 20 arcmin from $\omega$ Cen center. Similarly, the more metal-rich RGB stars show a more extended spatial distribution compared to the more metal-poor ones in the outskirts of the cluster. Moreover, the centers of the distributions of metal-rich and metal-poor RGB stars are shifted in different directions with respect to the geometrical center of $\omega$ Cen. We constructed stellar density profiles for the blue and red main-sequence stars; they confirm that the blue main-sequence sub-population has a more extended spatial distribution compared to the red main-sequence one in the outskirts of $\omega$ Cen, as found based on the star number ratio. We also computed the ellipticity profile of $\omega$ Cen, which has a maximum value of 0.16 at a distance of $\sim$ 8 arcmin from the center, and a minimum of 0.05 at $\sim$ 30 arcmin; the average ellipticity is $\sim0.10$. The circumstantial evidence presented in this work suggests a merging scenario for the formation of the peculiar stellar system $\omega$ Cen.

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, could be related to the initial configuration of the system. In recent work, we provided the first study of the kinematics of different stellar populations in NGC$\,$0104 over a large field of view in the plane of the sky, exploiting Gaia Data Release 2 (DR2) proper motions combined with multi-band ground-based photometry. In this paper, we combine Gaia DR2 proper motions with Very Large Telescope radial velocities to investigate the kinematics along the line of sight and in the plane of the sky of multiple populations in seven GCs, namely NGC$\,$0104, NGC$\,$0288, NGC$\,$5904, NGC$\,$6121, NGC$\,$6254, NGC$\,$6752 and NGC$\,$6838. Among the analyzed clusters only NGC$\,$0104 and NGC$\,$5904 show significant rotation. Separating our sample into two groups of first- and second-population stars (1P and 2P) we find that overall these two populations exhibit a similar rotation pattern in NGC$\,$0104. However, some hints of different rotation are observed in the external regions of this cluster. Interestingly, 1P and 2P stars in NGC$\,$5904 show different rotation curves, with distinct phases and such difference is significant at the $\sim$2.5-$\sigma$ level. The analysis of the velocity-dispersion profiles of multiple populations confirms that 2P stars of NGC$\,$0104 show stronger anisotropy than the 1P.

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 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 NGC0104 (47Tucanae) 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 NGC0288, NGC5904 (M5), NGC6121 (M4), NGC6752, NGC6838 (M71) and further explore NGC0104. Among the analyzed clusters only NGC0104 and NGC5904 show significant rotation on the plane of the sky. Interestingly, multiple stellar populations in NGC5904 exhibit different rotation curves.

We propose a new formation channel for intermediate mass black hole (IMBH) binaries via globular cluster collisions in the Galactic disc. Using numerical simulations, we show that the IMBHs form a tight binary that enters the gravitational waves (GWs) emission dominated regime driven by stellar interactions, and ultimately merge in $\lesssim 0.5$ Gyr. These events are clearly audible to LISA and can be associated with electromagnetic emission during the last evolutionary stages. During their orbital evolution, the IMBHs produce runaway stars comparable with GAIA and LAMOST observations.

Globular clusters (GCs), the oldest stellar systems observed in the Milky Way, have for long been considered single stellar populations. As such, they provided an ideal laboratory to understand stellar dynamics and primordial star formation processes. However, during the last two decades, observations unveiled their real, complex nature. Beside their pristine stars, GCs host one or more helium enriched and possibly younger stellar populations whose formation mechanism is still unknown. Even more puzzling is the existence of GCs showing star by star iron spreads. Using detailed N-body simulations we explore the hypothesis that these anomalies in metallicity could be the result of mutual stripping and mergers between a primordial population of disc GCs. In the first paper of this series we proved, both with analytical arguments and short-term N-body simulations, that disc GCs have larger fly-by and close encounter rates with respect to halo clusters. These interactions lead to mass exchange and even mergers that form new GCs, possibly showing metallicity spreads. Here, by means of long-term direct N-body simulations, we provide predictions on the dynamical properties of GCs that underwent these processes. The comparison of our predictions with available and future observational data could provide insights on the origin of GCs and on the Milky Way build-up history as a whole.

We present the results of a self-consistent $N$-body simulation following the evolution of a primordial population of thick disc globular clusters (GCs). We study how the internal properties of such clusters evolve under the action of mutual interactions, while they orbit a Milky Way-like galaxy. For the first time, through analytical and numerical considerations, we find that physical encounters between disc GCs are a crucial factor that contributed to the shape of the current properties of the Galactic GC system. Close passages or motion on similar orbits may indeed have a significant impact on the internal structure of clusters, producing multiple gravitationally bound sub-populations through the exchange of mass and even mergers. Our model produces two major mergers and a few small mass exchanges between pairs of GCs. Two of our GCs accrete stars from two companions, ending up with three internal sub-populations. We propose these early interactions and mergers between thick disc GCs with slightly different initial chemical compositions as a possible explanation for the presence of the spreads in metallicity observed in some of the massive Milky Ways GCs.

The series of events, which occurred at high redshift and originated multiple stellar populations in Globular Clusters (GCs) are still poorly understood. Theoretical work suggests that the present-day dynamics of stars in nearby GCs, including the rotation and velocity dispersion, may retain important clues on the formation of multiple populations. So far, the dynamics of multiple populations have been investigated either from radial velocities of a relatively-small sample of stars, or from relative proper motions of stars in the small field of view provided by the Hubble Space Telescope. In this context, Gaia provides the unique opportunity to investigate the dynamics of thousands GC stars over a wide field of view. For the first time, we combine Gaia DR2 proper motions and multi-band photometry to study the internal motions of the two main stellar populations of 47 Tucanae in a wide field of view. We confirm that this cluster exhibits high rotation on the plane of the sky and find that both stellar generations share similar rotation patters. Second-generation stars show stronger anisotropies and smaller tangential-velocity dispersion than the first generation, while there is no significant difference between their radial-velocity dispersion profiles. We discuss the impact of these results in the context of the formation scenarios for multiple stellar populations in GCs.

The Galactic center of the Milky-Way harbors a massive black hole (BH) orbited by a diverse population of young and old stars. A significant fraction of the youngest stars ($\sim4-7$ Myr) reside in a thin stellar disk with puzzling properties; the disk appears to be warped, shows asymmetries, and contains one or more clumpy structures (e.g. IRS 13). Models explaining the clumping invoked the existence of an intermediate-mass BH of $10^{3}-10^{4}$ M$_{\odot}$, but no kinematic evidence for such a BH has been found. Here we use extended $N$-body simulations and hybrid self-consistent field method models to show that naturally formed residual temporal asphericity of the hosting nuclear star cluster gives rise to torques on the disk, which lead to changes in its orientation over time, and to recurrent formation and dissolution of single spiral arm ($m=1$ modes) structures. The changing orientation leads to a flapping-like behavior of the disk and to the formation of a warped disk structure. The spiral arms may explain the over-densities in the disk (clumping) and its observed asymmetry, without invoking the existence of an intermediate-mass BH. The spiral arms are also important for the overall disk evolution and can be used to constrain the structure and composition of the nuclear stellar cluster.

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.

[abridged] In the typical giant-impact scenario for the Moon formation most of the Moon's material originates from the impactor. Any Earth-impactor composition difference should, therefore, correspond to a comparable Earth-Moon composition difference. Analysis of Moon rocks shows a close Earth-Moon composition similarity, posing a challenge for the giant-impact scenario, given that impactors were thought to significantly differ in composition from the planets they impact. Here we use a large set of 140 simulations to show that the composition of impactors could be very similar to that of the planets they impact; in $4.9\%$-$18.2\%$ ($1.9\%$-$6.7\%$) of the cases the resulting composition of the Moon is consistent with the observations of $\Delta^{17}O<15$ ($\Delta^{17}O<6$ ppm). These findings suggest that the Earth-Moon composition similarity could be resolved as to arise from the primordial Earth-impactor composition similarity. Note that although we find the likelihood for the suggested competing model of very high mass-ratio impacts (producing significant Earth-impactor composition mixing) is comparable ($<6.7\%$), this scenario also requires additional fine-tuned requirements of a very fast spinning Earth. Using the same simulations we also explore the composition of giant-impact formed Mars-moons as well as Vesta-like asteroids. We find that the Mars-moon composition difference should be large, but smaller than expected if the moons are captured asteroids. Finally, we find that the left-over planetesimals ('asteroids') in our simulations are frequently scattered far away from their initial positions, thus potentially explaining the mismatch between the current position and composition of the Vesta asteroid.

Nuclear Star Clusters (NSCs) are commonly observed in the centres of most galactic nuclei, including our own Milky Way. While their study can reveal important information about the build-up of the innermost regions of galaxies, the physical processes that regulate their formation are still poorly understood. NSCs might have been formed through gas infall and subsequent in situ star formation, and/or through the infall and merging of multiple star clusters into the centre of the galaxy. Here, we investigate the viability of the latter, by studying direct N-body simulations of inspiralling clusters to the centre of a Milky-Way-like nuclear bulge that hosts a massive black hole. We find that the NSC that forms through this process can show both morphological and kinematical properties that make it comparable with observations of the Milky Way NSC, including significant rotation- a fact that has been so far attributed mainly to gas infall. We explore its kinematic evolution, to see if and how the merger history can imprint fossil records on its dynamical structure. Moreover, we study the effect of stellar foreground contamination in the line-of-sight kinematics of the NSC. Our study shows that no fine tuning of the orientation of the infalling globular clusters is necessary to result in a rotating NSC. We suggest that cluster-inspiral is a viable mechanism for the formation of rotating NSCs.

Globular Clusters (GCs) and Nuclear Star Clusters (NSCs) are typically composed by several stellar populations, characterized by different chemical compositions. Different populations show different ages in NSCs but not necessarily in GCs. The youngest populations in NSCs appear to reside in disk-like structures, as observed in our Galaxy and in M31. Gas infall followed by formation of second generation (SG) stars in GCs may similarly form disk-like structures in the clusters nuclei. Here we explore this possibility and follow the long term evolution of stellar disks embedded in GCs, and study their effects on the evolution of the clusters. We study disks with different masses by means of detailed N-body simulations and explore their morphological and kinematic signatures on the GC structures. We find that as a second generation disk relaxes, the old, first generation, stellar population flattens and becomes more radially anisotropic, making the GC structure become more elliptical. The second generation stellar population is characterized by a lower velocity dispersion, and a higher rotational velocity, compared with the primordial older population. The strength of these kinematic signatures depends both on the relaxation time of the system and on the fractional mass of the second generation disk. We therefore conclude that SG populations formed in flattened configurations will give rise to two systematic trends: (1) Positive correlation between GC ellipticity and fraction of SG population (2) Positive correlation between GC relaxation time and ellipticity. Thereby GC ellipticities and rotation could be related to the formation of SG stars and their initial configuration.

Most of the properties of the Earth-Moon system can be explained by a collision between a planetary embryo and the growing Earth late in the accretion process. Simulations show that most of the material that eventually aggregates to form the Moon originates from the impactor. However, analysis of the terrestrial and lunar isotopic composition show them to be highly similar. In contrast, the compositions of other solar system bodies are significantly different than the Earth and Moon. This poses a major challenge to the giant impact scenario since the Moon-forming impactor is then thought to also have differed in composition from the proto-Earth. Here we track the feeding zones of growing planets in a suite of simulations of planetary accretion, in order to measure the composition of Moon-forming impactors. We find that different planets formed in the same simulation have distinct compositions, but the compositions of giant impactors are systematically more similar to the planets they impact. A significant fraction of planet-impactor pairs have virtually identical compositions. Thus, the similarity in composition between the Earth and Moon could be a natural consequence of a late giant impact.

Nuclear star clusters (NSCs) are dense stellar clusters observed in galactic nuclei, typically hosting a central massive black hole. Here we study the possible formation and evolution of NSCs through the inspiral of multiple star clusters hosting intermediate mass black holes (IMBHs). Using an N-body code we examine the dynamics of the IMBHs and their effects on the NSC. We find that IMBHs inspiral to the core of the newly formed NSC and segregate there. Although the IMBHs scatter each other and the stars, none of them is ejected from the NSC. The IMBHs are excited to high eccentricities and their radial density profile develops a steep power-law cusp. The stars also develop a power-law cusp (instead of the central core that forms in their absence), but with a shallower slope. The relaxation rate of the NSC is accelerated due to the presence of IMBHs, which act as massive-perturbers. This in turn fills the loss-cone and boosts the tidal disruption rate of stars both by the MBH and the IMBHs to a value excluded by rate estimates based on current observations. Rate estimates of tidal disruptions can therefore provide a cumulative constraint on the existence of IMBHs in NSCs.

Nuclear stellar cluster (NSCs) are known to exist around massive black holes (MBHs) in galactic nuclei. They are thought to have formed through in-situ star formation following gas inflow to the nucleus of the galaxy and/or through the infall of multiple stellar clusters. Here we study the latter, and explore the composite structure of the NSC, and its relation to the various stellar populations originating from its progenitor infalling clusters. We use N-body simulations of clusters infall, and show that this scenario may produce observational signatures in the form of age segregation: the distribution of the stellar properties (e.g. stellar age and/or metallicity) in the NSCs reflect the infall history of the different clusters. The stellar populations of clusters infalling at different times (dynamical ages), are differentially segregated in the NSC, and are not fully mixed even after few Gyrs of evolution. Moreover, the radial properties of stellar populations in the progenitor cluster are mapped to their radial distribution in the final NSC, potentially leading to efficient mass segregation in NSCs, even those where relaxation times are longer than a Hubble time. Finally, the overall structures of the stellar populations present non-spherical configurations and show significant cluster to cluster population differences.

Globular clusters (GCs) and many nuclear clusters (NCs) show evidence for hosting multiple generations of stellar populations. Younger stellar populations in NCs appear to reside in disk like structures, including the nuclear cluster in our own Galactic center as well as in M31. Kinematic studies of the anomalous globular cluster Omega Centauri, thought to possibly be a former dwarf galaxy (or a galactic nucleus), show evidence for its hosting of a central, kinematically cold disk component. These observations suggest that formation of second (or multiple) generation stars may occur in flattened disk like structures. Here we use detailed N-body simulations to explore the possible evolution of such stellar disks, embedded in globular clusters. We follow the long term evolution of a disk like structure similar to that observed in Omega Centauri and study its properties. We find that a stellar-disk like origin for second generation stellar populations can leave behind significant kinematic signatures in properties of the clusters, including an anisotropic distribution, and lower velocity dispersions, which can be used to constrain the origin of second generations stars and their dynamical evolution.

In this paper we report a study (see Mastrobuono-Battisti et al., 2012) about the formation and characteristics of the tidal tails around Palomar 5 along its orbit in the Milky Way potential, by means of direct N-body simulations and simplified numerical models. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested.

By means of direct N-body simulations and simplified numerical models, we study the formation and characteristics of the tidal tails around Palomar 5, along its orbit in the Milky Way potential. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested. Our models can reproduce the density contrast between the clumps and the surrounding tails found in the observed streams, without including any lumpiness in the dark halo, suggesting new upper limits on its granularity.

Abridged: In one widely discussed model for the formation of nuclear star clusters (NSCs), massive globular clusters spiral into the center of a galaxy and merge to form the nucleus. It is now known that at least some NSCs coexist with supermassive black holes (SBHs); this is the case, for instance, in the Milky Way (MW). In this paper, we investigate how the presence of a SMBH at the center of the MW impacts the merger hypothesis for the formation of its NSC. Starting from a model consisting of a low-density nuclear stellar disk and the SMBH, we use N-body simulations to follow the successive inspiral and merger of globular clusters. The clusters are started on circular orbits of radius 20 pc, and their initial masses and radii are set up in such a way as to be consistent with the galactic tidal field at that radius. The total accumulated mass by ~10 clusters is about 1.5x10^7 Solar masses. Each cluster is disrupted by the SMBH at a distance of roughly one parsec. The density profile that results after the final inspiral event is characterized by a core of roughly this radius, and an envelope with density that falls off as 1/r^2. These properties are similar to those of the MW NSC, with the exception of the core size, which in the MW is a little smaller. But by continuing the evolution of the model after the final inspiral event, we find that the core shrinks substantially via gravitational encounters in a time (when scaled to the MW) of 10 Gyr as the stellar distribution evolves toward a Bahcall-Wolf cusp. We also show that the luminosity function of the MW NSC is consistent with the hypothesis that a large fraction of the mass comes from (~10Gyr) old stars, brought in by globular clusters. We conclude that a model in which a large fraction of the mass of the MW NSC arose from infalling globular clusters is consistent with existing observational constraints.

Nuclear Star Clusters are observed at the center of many galaxies. In particular in the center of the Milky Way the Nuclear Star Cluster coexists with a cen- tral supermassive black hole. The origin of these clusters is still unknown; a possible formation mechanism is the decay of massive globular clusters driven inward to the galactic center by dynamical friction and their subsequent merging. By investigating this scenario by means of sophisticated N-body simulations we found that this process could lead to a final product which actually shows many of the observed features of the Milky Way Nuclear Star Cluster.

In the center of the Milky Way, as well as in many other galaxies, a compact star cluster around a very massive black hole is observed. One of the possible explanations for the formation of such Nuclear Star Clusters is based on the 'merging' of globular clusters in the inner galactic potential well. By mean of sophisticated N-body simulations, we checked the validity of this hypothesis and found that it may actually has been the one leading to the formation of the Milky Way Nuclear Star Cluster.

We present and discuss the characteristics and performances, both in term of computational speed and precision, of a numerical code which numerically integrates the equation of motions of N 'particles' interacting via Newtonian gravitation and move in an external galactic smooth field. The force evaluation on every particle is done by mean of direct summation of the contribution of all the other system's particle, avoiding truncation error. The time integration is done with second-order and sixth-order symplectic schemes. The code, NBSymple, has been parallelized twice, by mean of the Computer Unified Device Architecture to make the all-pair force evaluation as fast as possible on high-performance Graphic Processing Units NVIDIA TESLA C 1060, while the O(N) computations are distributed on various CPUs by mean of OpenMP Application Program. The code works both in single precision floating point arithmetics or in double precision. The use of single precision allows the use at best of the GPU performances but, of course, limits the precision of simulation in some critical situations. We find a good compromise in using a software reconstruction of double precision for those variables that are most critical for the overall precision of the code. The code is available on the web site astrowww.phys.uniroma1.it/dolcetta/nbsymple.html

In this paper we analyze data of 8 elliptical galaxies in order to study the difference between their globular cluster systems (GCSs) radial distributions and those of the galactic stellar component. In all the galaxies studied here the globular cluster system density profile is significantly flatter toward the galactic centre than that of stars. If this difference is interpreted as a depauperation of the initial GC population, the estimated number of missing globular clusters is significant, ranging from 21% to 71% of their initial population abundance in the eight galaxies examined. The corresponding mass lost to the central galactic region is 7x10^7-1.85x10^9 solar masses. All this mass carried toward central galactic regions have likely had an important feedback on the innermost galactic region, including its violent transient activity (AGN) and local massive black hole formation and growth.