The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We investigate the evolution of triple, non-hierarchical, black hole (BH) systems making use of $2.9\times10^4$ 3-body simulations. Varying the mutual orbital inclination, the three BH masses and the inner and outer eccentricities, we show that retrograde, nearly planar configurations lead to a significant shrinkage of the inner binary. We find an universal trend of triple systems, that they tend to evolve toward prograde configurations, Moreover, we demonstrate that the orbital flip, driven by the torque exerted on the inner BH binary (BHB) by the outer BH, leads in general to tighter inner orbits. In some cases, the resulting BHB undergoes coalescence within a Hubble time, releasing gravitational waves (GWs). Frequently, the inner BHB merger occurs after a component swap between one of its components and the outer BH. The mass spectrum of the BHBs that underwent the component exchange differs significantly from the case in which the BHB merge without any swap. A large fraction of merging BHBs with initial separation $1$ AU enter the $10^{-3}-10^{-1}$ Hz frequency band with large eccentricities, thus representing potential LISA sources. Mergers originating from initially tighter BHB ($a\sim 0.01$ AU), instead, have a large probability to have eccentricities above 0.7 in the $1$ Hz band. We find that the mergers' mass distribution in this astrophysical channel maps the original BH binary spectrum. This might have interesting consequences in light of the growing population of BH mergers detected by LIGO.

In this paper we study the evolution of a primordial black hole binary (BHB) in a sample of over 1500 direct-summation $N-$body simulations of small-and intermediate-size isolated star clusters as proxies of galactic open clusters. The BHBs have masses in the range of the first LIGO/Virgo detections. Some of our models show a significant hardening of the BHB in a relatively short time. Some of them merge within the cluster, while ejected binaries, typically, have exceedingly long merger timescales. The perturbation of stars around BHB systems is key to induce their coalescence. The BHBs which merge in the cluster could be detected with a delay of a few years between space detectors, as future LISA, and ground-based ones, due to their relatively high eccentricity. Under our assumptions, we estimate a BHB merger rate of $R_{\rm mrg} \sim 2$ yr$^{-1}$ Gpc$^{-3}$. We see that in many cases the BHB triggers tidal disruption events which, however, are not linked to the GW emission. Open cluster-like systems are, hence, a promising environment for GWs from BHBs and tidal disruptions.

We investigate the evolution of intermediate-mass (IMBHs), stellar (BHs) and binary black holes (BHBs), deposited near a supermassive black hole (SMBH) by a population of massive star clusters. Stellar BHs rapidly segregate around the SMBH, driving the formation of extreme mass-ratio inspirals that coalesce at a rate $\Gamma= 0.02-0.2$ yr$^{-1}$ Gpc$^{-3}$ at redshift $z=0$. A few IMBHs orbiting the SMBH favour the formation of massive pairs that coalescence within a Hubble time, being the merger rate for this channel $\Gamma =0.03$ yr$^{-1}$ Gpc$^{-3}$. Recoiling kicks post-merger can eject the remnant from the galaxy centre, especially in dwarf galaxies. Our results suggest that this mechanism can lead to up to $10^5$ ejected SMBH within 1 Gpc. An IMBH co-existing with a few single and binary BHs in the same cluster can affect significantly their evolution, either driving binary disruption, yielding to intermediate-mass ratio inspirals (merger rate $\Gamma =9.5$ yr$^{-1}$ Gpc$^{-3}$), or boosting BHBs coalescence ($\Gamma =2-8$ yr$^{-1}$ Gpc$^{-3}$). In a few simulations, the SMBH boosts BHBs coalescence, leading this process to a merger rate $\Gamma =1$ yr$^{-1}$ Gpc$^{-3}$. We note that BHBs experiencing a merger in a galactic nucleus can be erroneously estimated $\sim 30\%$ heavier than it really is because of the Doppler shift of the wave frequency as caused by the rapid motion around the SMBH. All our simulations are carried out using an $N$-body code tailored to treat close encounters and post-Newtonian dynamics, that includes also the galaxy field and dynamical friction in the particles' equation of motion.

The Milky Way centre exhibits an intense flux in the gamma and X-ray bands, whose origin is partly ascribed to the possible presence of a large population of millisecond pulsars (MSPs) and cataclysmic variables (CVs), respectively. However, the number of sources required to generate such an excess is much larger than what is expected from in situ star formation and evolution, opening a series of questions about the formation history of the Galactic nucleus. In this paper we make use of direct $N$-body simulations to investigate whether these sources could have been brought to the Galactic centre by a population of star clusters that underwent orbital decay and formed the Galactic nuclear star cluster (NSC). Our results suggest that the gamma ray emission is compatible with a population of MSPs that were mass segregated in their parent clusters, while the X-ray emission is consistent with a population of CVs born via dynamical interactions in dense star clusters. Combining observations with our modelling, we explore how the observed $\gamma$ ray flux can be related to different NSC formation scenarios. Finally, we show that the high-energy emission coming from the galactic central regions can be used to detect black holes heavier than $10^5\Ms$ in nearby dwarf galaxies.

The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.

We investigated the evolution of a massive galactic nucleus hosting a super-massive black hole (SMBH) with mass $M_\mathrm{SMBH}=10^8 \mathrm{M}_\odot$ surrounded by a population of 42 heavy star clusters (GCs). Using direct $N$-body modelling, we show here that the assembly of an NSC through GCs orbital decay and merger is efficiently inhibited by the tidal forces exerted from the SMBH. The GCs mass loss induced by tidal forces causes a significant modification of their mass function, leading to a population of low-mass ($<10^4$) clusters. Nonetheless, the GCs debris accumulated around the SMBH give rise to well-defined kinematical and morphological properties, leading to the formation of a disk-like structure. Interestingly, the disk is similar to the one observed in the M31 galaxy nucleus, which has properties similar to our numerical model. The simulation produced a huge amount of data, which we used to investigate whether the GC debris deposited around the SMBH can enhance the rate of tidal disruption events (TDEs) in our galaxy inner density distribution. Our results suggest that the GCs disruption shapes the SMBH neighbourhoods leading to a TDE rate of $\sim 2 \times 10^{-4} $yr$^{-1}$, a value slightly larger than what expected in previous theoretical modelling of galaxies with similar density profiles and central SMBHs. The simulation presented here is the first of its kind, representing a massive galactic nucleus and its star cluster population on scales $\sim 100$ pc.

One of the leading scenarios for the formation of nuclear star clusters in galaxies is related to the orbital decay of globular clusters (GCs) and their subsequent merging, though alternative theories are currently debated. The availability of high-quality data for GCs structural and orbital parameters allow to test different nuclear star cluster formation scenarios. The Fornax dwarf spheroidal (dSph) galaxy is the heaviest satellite of the Milky Way and it is the only known dwarf spheroidal hosting 5 GCs, whereas there are no clear signatures for the presence of a central massive black hole. For this reason, it represents a suited place to study the orbital decay process in dwarf galaxies. In this paper we model the future evolution of the Fornax GCs simulating them and the host galaxy by means of direct $N$-body simulations. Our simulations take in account also the gravitational field generated by the Milky Way. We found that if the Fornax galaxy is embedded in a standard Cold Dark Matter Halo, the nuclear cluster formation would be significantly hampered by the high central galactic mass density. In this context, we discuss the possibility that infalling GCs drive the flattening of the galactic density profile, giving a possible alternative explanation to the so-called cusp/core problem. Moreover, we briefly discuss the link between GC infall process and the absence of massive black holes in the centre of dSphs.

The Fornax dwarf spheroidal galaxy is the most massive satellites of the Milky Way, claimed to be embedded in a huge dark matter halo, and the only among the Milky Way satellites hosting five globular clusters. Interestingly, their estimated masses, ages and positions seem hardly compatible with the presence of a significant dark matter component, as expected in the $\Lambda$ CDM scheme. Indeed, if Fornax would have a CDM halo with a standard density profile, all its globular clusters should have sunk to the galactic centre many Gyr ago due to dynamical friction. Due to this, some authors proposed that the most massive clusters may have formed out of Fornax and later tidally captured. In this paper we investigate the past evolution of the Fornax GC system by using both a recently developed, semi-analytical treatment of dynamical friction and direct $N$-body simulations of the orbital evolution of the globular clusters within Fornax and of Fornax galaxy around the Milky Way. Our results suggest that an "in-situ" origin for all the clusters is likely if their observed positions are close to their spatial ones and their orbits are almost circular. Moreover, the Milky Way seems to accelerate the GC decay reducing the decay time of $15\%$. Nevertheless, our results indicate that the GCs survival probability exceeds $50\%$, even in the case of cuspy density profiles. We conclude that more detailed data are required to shed light on the Fornax dark matter content, to distinguish between a cuspy or a cored profile.

Many observed globular clusters (GCs) seem to show a central overabundance of mass whose nature has not yet fully understood. Indeed, it is not clear whether it is due to a central intermediate mass black hole (IMBH) or to a massive stellar system (MSS) composed of mass segregated stars. In this contribution we present a semi-analytic approach to the problem complemented by 12 $N$-body simulations in which we followed the formation of MSSs in GCs with masses up to $3\times 10^5$ \Ms. Some implications for the formation of IMBHs and gravitational waves emission are discussed in perspective of a future work.

We present a series of $N$-body simulations representing the evolution of a galactic nucleus and its stellar content in a nearly one-to-one representation. The aim of this suite of simulations is to shed light on the interplay between nuclear clusters (NCs), super-massive black holes (SMBH) and the galactic nuclei in which they are contained. We modelled galaxies with masses from few times $10^8$ to $10^{11}$ M$_\odot$, hosting in their nucleus a number of globular clusters and, in some cases, a central SMBH.

Recent observations of the Virgo cluster and the Local Group suggested that some galaxies are flowing out from their parent cluster. This may be the signature that dark energy (DE) acts significantly also on small cosmological scales. By means of direct N-body simulations we performed several simulations, in which the effect of DE and gravity are taken into account, aiming to determine whether DE can produce an outflow of galaxies compatible with observations. Comparing the different simulations, our results suggest that the observed outflow of galaxies is likely due to the local effect of DE.

We performed a series of direct N-body simulations with the aim to follow the dynamical evolution of a galaxy cluster (GC) ($M_{clus}\simeq 10^{14} M_{\odot}$) in different environment. The results show the formation of heavy sub-structures in the cluster centre in consequence of multiple merging among the innermost galaxies. Moreover we investigate the dynamics of super-massive black holes (SMBHs) residing in the centre of galaxies that form the most massive sub-structure.

Recently, observational data and high precision mapping of the local velocity field of Local Group and Virgo cluster have revealed a linear velocity-distance relation of the outermost galaxies, properly referred to as Local Hubble Flow. By means of direct N-body method, we performed several simulations in which a galaxy cluster undergoes the action of the Dark Energy force and of the gravitational one induced by the gas. We reproduced the so-called Hubble diagrams, to highlight the outflow of the galaxies lying in the external region of the cluster. Our preliminary results suggest that the observed outflow of galaxies is likely due to the local effect of Dark Energy. Furthermore, the accuracy of the N-body method used, allows us to follow the merging process among some galaxies with the aim to reproduce the formation of a single compact object in the centre of the cluster.

We present results of a set of $N$-body simulations to model the future evolution of the 11 young massive clusters hosted in the central region of the dwarf starburst galaxy Henize 2-10, which contains at its center a massive black hole with a mass $M_{\rm BH} \simeq 2\times 10^6$ M$_\odot$. Nuclear star clusters are present in a great quantity of galaxies of mass similar to Henize 2-10. Our results \citep{ASCD15} show that the orbital decay and merging of the Henize 2-10 clusters will likely lead to the formation of a nuclear star cluster with mass $M_{\rm NSC} \simeq 4-6 \times 10^6$ M$_\odot$ and effective radius $r_{\rm NSC} \simeq 4.1$ pc. Additionally, we found that this mechanism can lead to the formation of disky structures with global properties similar to those of nuclear stellar disks, which reside in many "middle-weight" galaxies. This work confirms and enlarge recent results \citep{ASCDS15} that indicate how nuclear star clusters and super massive black holes are only partially correlated, since the formation process of nuclear star clusters is poorly affected by a black hole of the size of that in Henize 2-10. A new result is that nuclear star clusters and nuclear stellar disks may share the same formation path.

We investigate the interaction of a massive globular cluster (GC) with a super massive black hole (SMBH), located at the centre of its host galaxy, by means of direct N-body simulations. The results show that tidal distortions induced by the stellar background and the SMBH act on a time shorter than that of dynamical friction decay for a $10^6$ M$_\odot$ GC whenever the SMBH mass exceeds $\sim 10^8$ M$_\odot$ . This implies an almost complete dissolution of the infalling GC before it reaches the inner region ($\lesssim 5$ pc) of the parent galaxy. The generalization of this result to a larger sample of infalling GCs shows that such destructive process may prevent the formation and growth of a bright galactic nucleus. Another interesting, serendipitous, result we obtained is that the close interaction between the SMBH and the GC produces a "wave" of stars that escape from the cluster and, in a fraction, even from the whole galaxy.

During their evolution, star clusters undergo mass segregation, by which the orbits of the most massive stars shrink, while the lighter stars move outwards from the cluster centre. In this context, recent observations and dynamical modelling of several galactic and extra-galactic globular clusters (GCs) suggest that most of them show, close to their centre, an overabundance of mass whose nature is still matter of debate. For instance, many works show that orbitally segregated stars may collide with each other in a runaway fashion, leading to the formation of a very massive star or an intermediate mass black hole (IMBH) with a mass comparable to the observed mass excess. On the other hand, segregated stars can form a dense system if the IMBH formation fails. In this paper we study the early formation phase of a dense, massive sub-system (MSS) in several GCs models using a recently developed semi-analytical treatment of the mass segregation process. In order to investigate how the MSS properties depend on the host cluster properties, we varied initial mass function (IMF), total mass, spatial distribution and metallicity of our models. Our results show how the IMF contributes to determine the final mass of the MSS, while the metallicity and the spatial distribution play a minor role. The method presented in this paper allowed us to provide scaling relations that connect the MSS mass and the host cluster mass in agreement with the observed correlation. In order to follow the early formation stage of the MSSs and improve our statistical results, we performed several $N$-body simulations of stellar clusters with masses between $10^3$ and $2\times 10^5$ solar masses.

The central region of the galaxy Henize 2-10 has a central black hole (BH) with a mass of about $2\times 10^6$ M$_\odot$. While this black hole does not appear to coincide with any central stellar over density, it is surrounded by 11 young massive clusters with masses above $10^5$ M$_\odot$. The availability of high quality data on the structure of the galaxy and the age and mass of the clusters provides excellent initial conditions for studying the dynamical evolution of Henize 2-10's nucleus. Here we present a set of $N$-body simulations of the central clusters and black hole to understand whether and how they will merge to form a nuclear star cluster. Nuclear star clusters (NSCs) are present in a majority of galaxies with stellar mass similar to Henize 2-10. Despite the results depend on the choice of initial conditions, we find that a NSC with mass $M_{NSC}\simeq 4-6\times 10^6$ M$_\odot$ and effective radius $r_{NSC}\simeq 2.6-4.1$ pc will form within $0.2$ Gyr. This work is the first showing, in a realistic realization of the host galaxy and its star cluster system, that the formation of a bright nucleus is a process that can happen after the formation of a central massive BH leading to a composite NSC+BH central system. The merging process of the clusters does not affect significantly the kinematics of the BH, whose motion, after the globular cluster merger, is limited to a $\sim 1$ pc oscillation at less than $2$ kms$^{-1}$ speed.

Almost all galaxies along the Hubble sequence host a compact massive object (CMO) in their center. The CMO can be either a supermassive black hole (SMBH) or a very dense stellar cluster, also known as nuclear star cluster (NSC). Generally, heavier galaxies (mass $\gtrsim 10^{11}\mathrm{M}_{\odot}$ ) host a central SMBH while lighter show a central NSC. Intermediate mass hosts, instead, contain both a NSC and a SMBH. One possible formation mechanisms of a NSC relies on the dry-merger (migratory) scenario, in which globular clusters (GCs) decay toward the center of the host galaxy and merge. In this framework, the absence of NSCs in high-mass galaxies can be imputed to destruction of the infalling GCs by the intense tidal field of the central SMBH. In this work, we report preliminary results of N-body simulations performed using our high-resolution, direct, code HiGPUs, to investigate the effects of a central SMBH on a single GC orbiting around it. By varying either the mass of the SMBH and the mass of the host galaxy, we derived an upper limit to the mass of the central SMBH, and thus to the mass of the host, above which the formation of a NSC is suppressed.

Nuclear Star Clusters (NSCs) are often present in spiral galaxies as well as resolved Stellar Nuclei (SNi) in elliptical galaxies centres. Ever growing observational data indicate the existence of correlations between the properties of these very dense central star aggregates and those of host galaxies, which constitute a significant constraint for the validity of theoretical models of their origin and formation. In the framework of the well known 'migratory and merger' model for NSC and SN formation, in this paper we obtain, first, by a simple argument the expected scaling of the NSC/SN mass with both time and parent galaxy velocity dispersion in the case of dynamical friction as dominant effect on the globular cluster system evolution. This generalizes previous results by \cite{TrOsSp} and is in good agreement with available observational data showing a shallow correlation between NSC/SN mass and galactic bulge velocity dispersion. Moreover, we give statistical relevance to predictions of this formation model, obtaining a set of parameters to correlate with the galactic host parameters. We find that the correlations between the masses of NSCs in the migratory model and the global properties of the hosts reproduce quite well the observed correlations, supporting the validity of the migratory-merger model. In particular, one important result is the flattening or even decrease of the value of the NSC/SN mass obtained by the merger model as function of the galaxy mass for high values of the galactic mass, i.e. $\gtrsim 3\times 10^{11}$M$_\odot$, in agreement with some growing observational evidence.

In this paper we treat the problem of the dynamical friction decay of a massive object moving in an elliptical galaxy with a cuspidal inner distribution of the mass density. We present results obtained by both self-consistent, direct summation, N-body simulations, as well as by a new semi-analytical treatment of dynamical friction valid in such cuspy central regions of galaxies. A comparison of these results indicates that the proposed semi-analytical approximation is the only reliable in cuspy galactic central regions, where the standard Chandrasekhar's local approximation fails, and, also, gives estimates of decay times that are correct at 1% respect to those given by N-body simulations. The efficiency of dynamical friction in cuspy galaxies is found definitively higher than in core galaxies, especially on more radially elongated satellite orbits. As another relevant result, we find a proportionality of the dynamical friction decay time to the -0.67 power of the satellite mass, M, shallower than the standardly adopted 1/M dependence.

Dynamical friction is the process responsible for matter transport toward the inner regions of galaxies in form of massive objects, like intermediate mass black holes, globular clusters and small satellite galaxies. While very bright galaxies show an almost flat luminosity profile in the inner region, fainter ones have, usually, a peaked, cuspidal, profile toward the center. This makes unreliable, in these cases, the use of the classic Chandrasekhar's formula for dynamical friction in its local approximation. Using both N-body simulations and a semi analytical approach, we have obtained reliable results for the orbital decay of massive objects in cuspidal galaxies. A relevant result is that of a shallower dependence of dynamical friction braking on the satellite mass than in the usual Chandrasekhar's local expression, at least in a range of large satellite masses.

The central regions of galaxies show the presence of massive black holes and/or dense stellar systems. The question about their modes of formation is still under debate. A likely explanation of the formation of the central dense stellar systems in both spiral and elliptical galaxies is based on the orbital decay of massive globular clusters in the central region of galaxies due to kinetic energy dissipation by dynamical friction. Their merging leads to the formation of a nuclear star cluster, like that of the Milky Way, where a massive black hole (Sgr A*) is also present. Actually, high precision N-body simulations (Antonini, Capuzzo-Dolcetta et al. 2012, ApJ, 750, 111) show a good fit to the observational characteristics of the Milky Way nuclear cluster, giving further reliability to the cited `migratory' model for the formation of compact systems in the inner galaxy regions.

It is well known that a large fraction of galaxies have cuspy luminosity profiles in their central regions, at least within the observational resolution. In such cases, the often used, simplified, local approximation for the dynamical friction braking classical term fails when the massive satellite moves through the inner parts of the galaxy, although the scattering integral still converges for phase space distribution singularities that are not too sharp. Here we present preliminary results of our work aiming at finding better and more reliable results from the integration of motion of massive objects (globular clusters) in galaxies where the density diverges to the center in a power law form, with exponent greater than -2.