Nuclear Star Cluster Research Articles

Constraining nuclear star cluster formation using MUSE-AO observations of the early-type galaxy FCC 47

Context. Nuclear star clusters (NSCs) are found in at least 70% of all galaxies, but their formation path is still unclear. In the most common scenarios, NSCs form in-situ from the galaxy’s central gas reservoir, through the merging of globular clusters (GCs), or through a combination of both.Aims. As the scenarios pose different expectations for angular momentum and stellar population properties of the NSC in comparison to the host galaxy and the GC system, it is necessary to characterise the stellar light, NSC, and GCs simultaneously. The large NSC (reff = 66 pc) and rich GC system of the early-type Fornax cluster galaxy FCC 47 (NGC 1336) render this galaxy an ideal laboratory to constrain NSC formation.Methods. Using Multi Unit Spectroscopic Explorer science verification data assisted by adaptive optics, we obtained maps for the stellar kinematics and stellar-population properties of FCC 47. We extracted the spectra of the central NSC and determined line-of-sight velocities of 24 GCs and metallicities of five.Results. The galaxy shows the following kinematically decoupled components (KDCs): a disk and a NSC. Our orbit-based dynamical Schwarzschild model revealed that the NSC is a distinct kinematic feature and it constitutes the peak of metallicity and old ages in FCC 47. The main body consists of two counter-rotating populations and is dominated by a more metal-poor population. The GC system is bimodal with a dominant metal-poor population and the total GC system mass is ∼17% of the NSC mass (∼7 × 108 M⊙).Conclusions. The rotation, high metallicity, and high mass of the NSC cannot be explained by GC-inspiral alone. It most likely requires additional, quickly quenched, in-situ formation. The presence of two KDCs likely are evidence of a major merger that has significantly altered the structure of FCC 47, indicating the important role of galaxy mergers in forming the complex kinematics in the galaxy-NSC system.

Asymmetric Dark Matter Imprint on Low-mass Main-sequence Stars in the Milky Way Nuclear Star Cluster

Abstract In this work, we study the impact of asymmetric dark matter (ADM) on low-mass main-sequence stars in the Milky Way’s nuclear star cluster, where the dark matter (DM) density is expected to be orders of magnitude above what is found near the Sun ( ). Using a modified stellar evolution code and considering a DM particle (m χ = 4 GeV) with a spin-dependent interaction cross section close to the limits allowed by direct detection, we found that the interactions of ADM with baryons in the star’s core can have two separate effects on the evolution of these stars: a decrease in the hydrogen burning rate, extending the duration of the main-sequence of stars with M ∼ 1M ⊙ by a few Gyr; the suppression of the onset of convection in the core of stars with M ≲ 1.5M ⊙ and consequent quench of supply for the nuclear reactions. If we consider ρ DM > 103 GeV cm−3 (corresponding to the inner 5 pc of the Milky Way), stars lighter than the Sun will have a main-sequence life span comparable to the current age of the universe. Stars heavier than two solar masses are not sensitive to the DM particles considered here.

The inner two degrees of the Milky Way

Context. Although there have been numerous studies of chemical abundances in the Galactic bulge, the central two degrees have been relatively unexplored due to the heavy and variable interstellar extinction, extreme stellar crowding, and the presence of complex foreground disk stellar populations. Aims. In this paper we discuss the metallicity distribution function, vertical and radial gradients, and chemical abundances of α-elements in the inner two degrees of the Milky Way, as obtained by recent IR spectroscopic surveys. Methods. We used a compilation of recent measurements of metallicities and α-element abundances derived from medium-high resolution spectroscopy. We compare these metallicities with low-resolution studies. Results. Defining “metal-rich” as stars with [Fe/H] > 0, and “metal-poor” as stars with [Fe/H] < 0, we find compelling evidence for a higher fraction (∼80%) of metal-rich stars in the Galactic Center (GC) compared to the values (50–60%) measured in the low latitude fields within the innermost 600 pc. The high fraction of metal-rich stars in the GC region implies a very high mean metallicity of +0.2 dex, while in the inner 600 pc of the bulge the mean metallicity is rather homogenous around the solar value. A vertical metallicity gradient of −0.27 dex kpc−1 in the inner 600 pc is only measured if the GC is included, otherwise the distribution is about flat and consistent with no vertical gradient. Conclusions. In addition to its high stellar density, the Galactic center/nuclear star cluster is also extreme in hosting high stellar abundances, compared to the surrounding inner bulge stellar populations; this has implications for formation scenarios and strengthens the case for the nuclear star cluster being a distinct stellar system.

The Next Generation Virgo Cluster Survey. XXIII. Fundamentals of Nuclear Star Clusters over Seven Decades in Galaxy Mass

Abstract Using deep, high-resolution optical imaging from the Next Generation Virgo Cluster Survey, we study the properties of nuclear star clusters (NSCs) in a sample of nearly 400 quiescent galaxies in the core of Virgo with stellar masses 105 ≲ / ≲ 1012. The nucleation fraction reaches a peak value f n ≈ 90% for ≈ 109 galaxies and declines for both higher and lower masses, but nuclei populate galaxies as small as ≈ 5 × 105 . Comparison with literature data for nearby groups and clusters shows that at the low-mass end nucleation is more frequent in denser environments. The NSC mass function peaks at M NSC ≈ 7 × 105 , a factor 3–4 times larger than the turnover mass for globular clusters (GCs). We find a nonlinear relation between the stellar masses of NSCs and those of their host galaxies, with a mean nucleus-to-galaxy mass ratio that drops to M NSC/M ⋆ ≈ 3.6 × 10−3 for ≈ 5 × 109 galaxies. Nuclei in both more and less massive galaxies are much more prominent: at the low-mass end, where nuclei are nearly 50% as massive as their hosts. We measure an intrinsic scatter in NSC masses at a fixed galaxy stellar mass of 0.4 dex, which we interpret as evidence that the process of NSC growth is significantly stochastic. At low galaxy masses we find a close connection between NSCs and GC systems, including very similar occupation distributions and comparable total masses. We discuss these results in the context of current dissipative and dissipationless models of NSC formation.

On the Orbital Decay of Globular Clusters in NGC 1052-DF2: Testing a Baryon-only Mass Model

Abstract The dark matter content of the ultra-diffuse galaxy NGC 1052-DF2, as inferred from globular cluster (GC) and stellar kinematics, carries a considerable amount of uncertainty, with current constraints also allowing for the complete absence of dark matter. We test the viability of such a scenario by examining whether in a “baryon-only” mass model the observed GC population experiences rapid orbital decay due to dynamical friction. Using a suite of 50 multi-GC N-body simulations that match observational constraints on both the stellar component of NGC 1052-DF2 and its GC population but differ in the initial line-of-sight positions and the tangential velocities of the GCs, we show that there is a substantial amount of realization-to-realization variance in the evolution of the GCs. Nevertheless, over ∼10 Gyr, some of the GCs experience significant orbital evolution. Others evolve less. A combination of reduced dynamical friction in the galaxy core and GC–GC scattering keeps the GCs afloat, preventing them from sinking all the way to the galaxy center. While the current phase-space coordinates of the GCs are not unlikely for a baryon-only mass model, the GC system does evolve over time. Therefore, if NGC 1052-DF2 has no dark matter, some of its GCs must have formed farther out, and the GC system must have been somewhat more extended in the past. The presence of a low-mass cuspy halo, while allowed by the kinematics, seems improbable, as significantly shorter inspiral timescales in the central region would quickly lead to the formation of a nuclear star cluster.

Black hole growth through hierarchical black hole mergers in dense star clusters: implications for gravitational wave detections

In a star cluster with a sufficiently large escape velocity, black holes (BHs) that are produced by BH mergers can be retained, dynamically form new BH binaries, and merge again. This process can repeat several times and lead to significant mass growth. In this paper, we calculate the mass of the largest BH that can be formed through repeated mergers of stellar seed BHs and determine how its value depends on the physical properties of the host cluster. We adopt an analytical model in which the energy generated by the black hole binaries in the cluster core is assumed to be regulated by the process of two-body relaxation in the bulk of the system. This principle is used to compute the hardening rate of the binaries and to relate this to the time-dependent global properties of the parent cluster. We demonstrate that in clusters with initial escape velocity $\gtrsim 300\rm km\ s^{-1}$ in the core and density $\gtrsim 10^5\ M_\odot\rm pc^{-3}$, repeated mergers lead to the formation of BHs in the mass range $100-10^5 \,M_\odot$, populating any upper mass gap created by pair-instability supernovae. This result is independent of cluster metallicity and the initial BH spin distribution. We show that about $10\%$ of the present-day nuclear star clusters meet these extreme conditions, and estimate that BH binary mergers with total mass $\gtrsim 100\,M _\odot$ should be produced in these systems at a maximum rate $\approx 0.05 \,\rm Gpc^{-3} yr^{-1}$, corresponding to one detectable event every few years with Advanced LIGO/VIRGO at design sensitivity. The contribution of globular clusters is likely to be negligible instead because the first BH merger remnant escapes following the relativistic kick. A possible connection of our results to the formation of massive BH seeds in galaxy nuclei and globular clusters is discussed.

The global stability of M33: still a puzzle

The inner disc of the local group galaxy M33 appears to be in settled rotational balance, and near IR images reveal a mild, large-scale, two-arm spiral pattern with no strong bar. We have constructed N-body models that match all the extensive observational data on the kinematics and surface density of stars and gas in the inner part of M33. We find that currently favoured models are unstable to the formation of a strong bar of semi-major axis 2 < a_B < 3 kpc within 1 Gyr, which changes the dynamical properties of the models to become inconsistent with the current, apparently well-settled, state. The formation of a bar is unaffected by how the gas component is modelled, by increasing the mass of the nuclear star cluster, or by making the dark matter halo counter-rotate, but it can be prevented by either reducing the mass-to-light ratio of the stars to Upsilon_V ~ 0.6 or Upsilon_K ~ 0.23 in solar units or by increasing the random motions of the stars. Also a shorter and weaker bar results when the halo is rigid and unresponsive. However, all three near-stable models support multi-arm spirals, and not the observed large-scale bi-symmetric spiral. A two-arm spiral pattern could perhaps be tidally induced, but such a model would require an implausibly low mass disc to avoid a bar and there is no visible culprit. Thus the survival of the current state of this exceptionally well-studied galaxy is not yet understood. We also suspect that many other unbarred galaxies present a similar puzzle.

The properties of nuclear star clusters and their host galaxies

AbstractNuclear star clusters are found at the centers of most galaxies. They are the densest stellar systems in the Universe, and thus have unique and interesting stellar dynamics. We review how common nuclear star clusters are in galaxies of different masses and types, and then discuss the typical properties of NSCs. We close by discussing the formation of NSCs, and how a picture is emerging of different formation mechanisms being dominant in lower and higher mass galaxies.

Using MUSE-AO observations to constrain the formation of the large nuclear star cluster in FCC 47

AbstractNuclear star clusters (NSCs) are found in at least 70% of all galaxies, but their formation path is still unclear. In the most common scenarios, NSCs form in-situ from the galaxy’s central gas reservoir, through merging of globular clusters (GCs), or through a combination of the two. As the scenarios pose different expectations for angular momentum and stellar population properties of the NSC in comparison to the host galaxy and the GC system, it is necessary to characterise the stellar light, NSC, and GCs simultaneously. Wide-field observations with modern integral field units such as the Multi Unit Spectroscopic Explorer (MUSE) allow to perform such studies. However, at large distances, NSCs usually are not resolved in MUSE observations. The particularly large NSC (Reff ∼ 66 pc) of the early-type galaxy FCC 47 at distance of ∼20 Mpc is an exception and is therefore an ideal laboratory to constrain NSC formation of external galaxies.

Dark Matter: An Efficient Catalyst for Intermediate-mass-ratio-inspiral Events

Abstract Gravitational waves (GWs) can be produced if a stellar compact object, such as a black hole (BH) or neutron star, inspirals into an intermediate-massive black hole of . Such a system may be produced in the center of a globular cluster or a nuclear star cluster and is known as an intermediate- or extreme-mass-ratio inspiral (IMRI or EMRI). Motivated by the recent suggestions that dark matter minispikes could form around IMBHs, we study the effect of dynamical friction against DM on the merger rate of IMRIs/EMRIs. We find that the merger timescale of IMBHs with BHs and NSs would be shortened by two to three orders of magnitude. As a result, the event rate of IMRIs/EMRIs are enhanced by orders of magnitude relative to that in the case of no DM minispikes. In the most extreme case where IMBHs are small and the DM minispikes have a steep density profile, all the BH in GCs and NSCs might be exhausted so that the mergers with NSs would dominate the current IMRIs/EMRIs. Our results suggest that the mass function of the IMBHs below 104 M ⊙ would bear imprints of the distribution of DM minispikes because these low-mass IMBHs can grow efficiently in the presence of DM minispikes by merging with BHs and NSs. Future space-based GW detectors, like LISA, Taiji, and Tianqin, can measure the IMRI/EMRI rate and hence constrain the distribution of DM around IMBHs.

Constraining the stellar mass function from the deficiency of tidal disruption flares in the nuclei of massive galaxies

ABSTRACT The rate of tidal disruption flares (TDFs) per mass of the disrupting black hole encodes information on the present-day mass function (PDMF) of stars in the clusters surrounding super massive black holes. We explore how the shape of the TDF rate with black hole mass can constrain the PDMF, with only weak dependence on black hole spin. We show that existing data can marginally constrain the minimum and maximum masses of stars in the cluster, and the high-mass end of the PDMF slope, as well as the overall TDF rate. With $\mathcal {O}(100)$ TDFs expected to be identified with the Zwicky Transient Facility, the overall rate can be highly constrained, but still with only marginal constraints on the PDMF. However, if ${\lesssim } 10 {{\ \rm per\ cent}}$ of the TDFs expected to be found by LSST over a decade ($\mathcal {O}(10^3)$ TDFs) are identified, then precise and accurate estimates can be made for the minimum stellar mass (within a factor of 2) and the average slope of the high-mass PDMF (to within $\mathcal {O}(10{{\ \rm per\ cent}})$) in nuclear star clusters. This technique could be adapted in the future to probe, in addition to the PDMF, the local black hole mass function and possibly the massive black hole binary population.

The stellar remnants of high-redshift nuclear starburst discs: a potential origin for nuclear star clusters?

Nuclear starburst discs (NSDs) are very compact star-forming regions in the centers of galaxies that have been studied as a possible origin for the absorbing gas around a central active galactic nucleus. NSDs may be most relevant at $z\sim 1$ when obscured accretion onto supermassive black holes (SMBHs) is common. This paper describes the characteristics of the stellar remnants of NSDs at $z=0.01$, taking into account the evolution from $z=1$. Using a stellar synthesis model, the colours, masses, and luminosities of the stellar remnants are computed for a suite of 192 two-dimensional NSD models. These properties are compared to observations of local nuclear star clusters (NSCs), and a good match is found between the predicted and observed properties. Dynamical effects will likely cause the final remnant to be a rotating, nearly spherical distribution. In addition, $\approx 20$% of the NSD remnants have half-light radii <~ 10 pc, consistent with NSCs hosted in both late-type and early-type galaxies, and all the remnants follow similar size-luminosity relationships as observed in nearby NSCs. NSDs require the presence of a central SMBH and the most massive and compact stellar remnants are associated with the most massive SMBHs, although stellar clusters with a variety of sizes can be produced by all considered SMBH masses. Overall, NSDs at $z\sim 1$ appear to be a promising origin for the $\gg 1$ Gyr NSC population in early- and late-type galaxies with large SMBHs.

Nuclear Star Clusters in Coma confirmation of an unusually high nucleation fraction

AbstractNuclear star clusters (NSCs) are stellar systems similar in size to globular clusters (GCs) but extremely dense, comparable only to some GCs and ultra-compact dwarfs. They are present in galaxies with a wide range of masses, morphologies and gas content. There are several formation scenarios proposed for the formation of such objects, such as the merger of GCs or extreme star formation caused by the inflow of gas. Recent studies show that the presence of an NSC is related to galaxy stellar mass. Moreover, it has been suggested that NSCs are more often found in high density environments. In our work, we use deep imaging of the core regions of the Coma cluster down to an absolute magnitude of –8.2 and found that in this environment the nucleation fraction is higher than in the Virgo and Fornax clusters. We find nucleated galaxies in Coma as faint as –11.2 mag.

The Galactic Center: Improved Relative Astrometry for Velocities, Accelerations, and Orbits near the Supermassive Black Hole

Abstract We present improved relative astrometry for stars within the central half parsec of our Galactic Center (GC) based on data obtained with the 10 m W. M. Keck Observatory from 1995 to 2017. The new methods used to improve the astrometric precision and accuracy include correcting for local astrometric distortions, applying a magnitude-dependent additive error, and more carefully removing instances of stellar confusion. Additionally, we adopt jackknife methods to calculate velocity and acceleration uncertainties. The resulting median proper motion uncertainty is 0.05 mas yr−1 for our complete sample of 1184 stars in the central 10″ (0.4 pc). We have detected 24 accelerating sources, 2.6 times more than the number of previously published accelerating sources, which extend out to 4″ (0.16 pc) from the black hole. Based on S0-2's orbit, our new astrometric analysis has reduced the systematic error of the supermassive black hole (SMBH) by a factor of 2. The linear drift in our astrometric reference frame is also reduced in the north–south direction by a factor of 4. We also find the first potential astrometric binary candidate S0-27 in the GC. These astrometric improvements provide a foundation for future studies of the origin and dynamics of the young stars around the SMBH, the structure and dynamics of the old nuclear star cluster, the SMBH’s properties derived from orbits, and tests of general relativity in a strong gravitational field.

The Impact of Stripped Nuclei on the Supermassive Black Hole Number Density in the Local Universe

Abstract The recent discovery of supermassive black holes (SMBHs) in high mass ultra-compact dwarf galaxies (UCDs) suggests that at least some UCDs are the nuclear star clusters of stripped galaxies. In this paper we present a new method to estimate how many UCDs host an SMBH and thus are stripped galaxy nuclei. We revisit the dynamical mass measurements that suggest many UCDs have more mass than expected from stellar population estimates, which observations have shown is due to the presence of an SMBH. We revise the stellar population mass estimates using a new empirical relation between the mass-to-light ratio (M/L) and metallicity to predict which UCDs most likely host an SMBH. We calculate the fraction of UCDs that host SMBHs across their entire luminosity range for the first time. We then apply the SMBH occupation fraction to the observed luminosity function of UCDs and estimate that in the Fornax and Virgo clusters alone there should be stripped nuclei with SMBHs. This analysis shows that stripped nuclei are almost as common in clusters as present-day galaxy nuclei. We estimate the SMBH number density caused by stripped nuclei to be (2–8) × 10−3 Mpc−3, which represents a significant fraction (8%–32%) of the SMBH density in the local universe. These SMBHs hidden in stripped nuclei increase expected event rates for tidal disruption events and SMBH–SMBH and SMBH–BH mergers. The existence of numerous stripped nuclei with SMBHs are a direct consequence of hierarchical galaxy formation, but until now their impact on the SMBH density had not been quantified.

The Milky Way like galaxy NGC 6384 and its nuclear star cluster at high NIR spatial resolution using LBT/ARGOS commissioning data

We analyse high spatial resolution near infra-red (NIR) imaging of NGC6384, a Milky Way like galaxy, using ARGOS commissioning data at the Large Binocular Telescope (LBT). ARGOS provides a stable PSF$_{\rm FWHM}\!=\!0.2"\!-\!0.3"$ AO correction of the ground layer across the LUCI2 NIR camera $4'\!\times4'$ field by using six laser guide stars (three per telescope) and a natural guide star for tip-tilt sensing and guiding. Enabled by this high spatial resolution we analyse the structure of the nuclear star cluster (NSC) and the central kiloparsec of NGC6384. We find via 2D modelling that the NSC ($r_{\rm eff}\!\simeq\!10$pc) is surrounded by a small ($r_{\rm eff}\!\simeq\!100$pc) and a larger Sersi\'c ($r_{\rm eff}\!\simeq\!400$pc), all embedded within the NGC\,6384 large-scale boxy/X-shaped bulge and disk. This proof-of-concept study shows that with the high spatial resolution achieved by ground-layer AO we can push such analysis to distances previously only accessible from space. SED-fitting to the NIR and optical HST photometry allowed to leverage the age-metallicity-extinction degeneracies and derive the effective NSC properties of an young to old population mass ratio of $8\%$ with ${\cal M}_{\rm\star,old}\!\simeq\!3.5\times10^7M_\odot$, Age$_{\rm old,\ young}\!=\!10.9\pm1.3$Gyr and 226\,Myr $\pm62\%$, metallicity [M/H]$=\!-0.11\pm0.16$ and $0.33\pm39\%$dex, and $E(B\!-\!V)\!=\!0.63$ and 1.44mag.

Evolution of supermassive black hole binaries and tidal disruption rates in non-spherical galactic nuclei

Binary supermassive black holes (SMBH) are expected to form naturally during galaxy mergers. After the dynamical friction phase, when the two SMBHs become gravitationally bound to each other, and a brief stage of initial rapid hardening, the orbit gradually continues to shrink due to three-body interactions with stars that enter the loss cone of the binary. Using the stellar-dynamical Monte Carlo code Raga, we explore the co-evolution of the binary SMBH and the nuclear star cluster in this slow stage, for various combinations of parameters (geometry of the star cluster, primary/secondary SMBH mass, initial eccentricity, inclusion of stellar captures/tidal disruptions). We compare the rates of stellar captures in galactic nuclei containing a binary to those of galaxies with a single SMBH. At early times, the rates are substantially higher in the case of a binary SMBH, but subsequently they drop to lower levels. Only in triaxial systems both the binary hardening rates and the capture rates remain sufficiently high during the entire evolution. We find that the hardening rate is not influenced by star captures, nor does it depend on eccentricity; however, it is higher when the difference between the black hole masses is greater. We confirm that the eccentricity of the binary tends to grow, which may significantly shorten the coalescence time due to earlier onset of gravitational-wave emission. We also explore the properties of the orbits entering the loss cone, and demonstrate that it remains partially full throughout the evolution in the triaxial case, but significantly depleted in the axisymmetric case. Finally, we study the distribution of ejected hypervelocity stars and the corresponding mass deficits in the central parts of the galaxies hosting a binary, and argue that the missing mass is difficult to quantify observationally.

An Adaptive Optics Survey of Stellar Variability at the Galactic Center

Abstract We present an ≈11.5 yr adaptive optics (AO) study of stellar variability and search for eclipsing binaries in the central ∼0.4 pc (∼10″) of the Milky Way nuclear star cluster. We measure the photometry of 563 stars using the Keck II NIRC2 imager (K′-band, λ 0 = 2.124 μm). We achieve a photometric uncertainty floor of Δm K′ ∼ 0.03 (≈3%), comparable to the highest precision achieved in other AO studies. Approximately half of our sample (50% ± 2%) shows variability: 52% ± 5% of known early-type young stars and 43% ± 4% of known late-type giants are variable. These variability fractions are higher than those of other young, massive star populations or late-type giants in globular clusters, and can be largely explained by two factors. First, our experiment time baseline is sensitive to long-term intrinsic stellar variability. Second, the proper motion of stars behind spatial inhomogeneities in the foreground extinction screen can lead to variability. We recover the two known Galactic center eclipsing binary systems: IRS 16SW and S4-258 (E60). We constrain the Galactic center eclipsing binary fraction of known early-type stars to be at least 2.4% ± 1.7%. We find no evidence of an eclipsing binary among the young S-stars nor among the young stellar disk members. These results are consistent with the local OB eclipsing binary fraction. We identify a new periodic variable, S2-36, with a 39.43 days period. Further observations are necessary to determine the nature of this source.

How nucleation and luminosity shape faint dwarf galaxies

ABSTRACT We study the intrinsic shapes of a sample of over 400 quiescent galaxies in the cores of the Virgo and Fornax clusters with luminosities 106 ≤ Lg/L⊙ ≤ 108. Similar to satellites of the Local Group and Centaurus A, these faint, low surface brightness cluster galaxies are best described as a family of thick (〈C/A〉 &amp;gt; 0.5), oblate-triaxial spheroids. However, the large sample size allows us to show that the flattening of their stellar distribution depends both on luminosity and on the presence of a nuclear star cluster. Nucleated satellites are thicker at all luminosities compared to their non-nucleated counterparts, and fainter galaxies are systematically thicker as well, regardless of nucleation. Once nucleation is accounted for, we find no evidence that the environment the satellites live in plays a relevant role in setting their 3D structure. We interpret both the presence of stellar nuclei and the associated thicker shapes as the result of preferential early and rapid formation, effectively making these faint nucleated galaxies the first generation of cluster satellites.

DirectN-body simulation of the Galactic centre

We study the dynamics and evolution of the Milky Way nuclear star cluster performing a high resolution direct one-million-body simulation. Focusing on the interactions between such stellar system and the central supermassive black hole, we find that different stellar components adapt their overall distribution differently. After 5 Gyr, stellar mass black holes are characterized by a spatial distribution with power-slope $-1.75$, fully consistent with the prediction of Bahcall-Wolf pioneering work. Using the vast amount of data available, we infer the rate for tidal disruption events, being $4 \times 10^{-6}$ per yr, and estimate the number of objects that emit gravitational waves during the phases preceding the accretion onto the super-massive black hole, $\sim 270$ per Gyr. We show that some of these sources could form extreme mass-ratio inspirals. We follow the evolution of binary stars population, showing that the initial binary fraction of $5\%$ drops down to $2.5\%$ inside the inner parsec. Also, we explored the possible formation of binary systems containing a compact object, discussing the implications for millisecond pulsars formation and the development of Ia Supernovae.