The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

Binary black holes (BBHs) are thought to form in different environments, including the galactic field and (globular, nuclear, young and open) star clusters. Here, we propose a method to estimate the fingerprints of the main BBH formation channels associated with these different environments. We show that the metallicity distribution of galaxies in the local Universe along with the relative amount of mergers forming in the field or in star clusters determine the main properties of the BBH population. Our fiducial model predicts that the heaviest merger to date, GW170729, originated from a progenitor that underwent 2--3 merger events in a dense star cluster, possibly a galactic nucleus. The model predicts that at least one merger remnant out of 100 BBH mergers in the local Universe has mass $90 < M_{\rm rem}/ {\rm ~M}_\odot \leq{} 110$, and one in a thousand can reach a mass as large as $M_{\rm rem} \gtrsim 250$ M$_\odot$. Such massive black holes would bridge the gap between stellar-mass and intermediate-mass black holes. The relative number of low- and high-mass BBHs can help us unravelling the fingerprints of different formation channels. Based on the assumptions of our model, we expect that isolated binaries are the main channel of BBH merger formation if $\sim 70\%$ of the whole BBH population has remnants masses $<50$ M$_\odot$, whereas $\gtrsim{}6$\% of remnants with masses $>75$ M$_\odot$ point to a significant sub-population of dynamically formed BBH binaries.

The detection of gravitational waves emitted during a neutron star - black hole merger and the associated electromagnetic counterpart will provide a wealth of information about stellar evolution nuclear matter, and General Relativity. While the theoretical framework about neutron star - black hole binaries formed in isolation is well established, the picture is loosely constrained for those forming via dynamical interactions. Here, we use N-body simulations to show that mergers forming in globular and nuclear clusters could display distinctive marks compared to isolated mergers, namely larger masses, heavier black holes, and the tendency to have no associated electromagnetic counterpart. These features could represent a useful tool to interpreting forthcoming observations. In the Local Universe, gravitational waves emitted from dynamical mergers could be unravelled by detectors sensitive in the decihertz frequency band, while those occurring at the distance range of Andromeda and the Virgo Cluster could be accessible to lower-frequency detectors like LISA.

In this paper, we explore the mechanisms that regulate the formation and evolution of stellar black hole binaries (BHBs) around supermassive black holes (SMBHs). We show that dynamical interactions can efficiently drive "in-situ" BHB formation if the SMBH is surrounded by a massive nuclear cluster (NC), while orbitally segregated star clusters can replenish the BHB reservoir in SMBH-dominated nuclei. We discuss how the combined action of stellar hardening and mass segregation sculpts the BHB orbital properties. We use direct N-body simulations including post-Newtonian corrections up to 2.5 order to study the BHB-SMBH interplay, showing that the Kozai-Lidov mechanism plays a crucial role in shortening binaries lifetime. We find that the merging probability weakly depends on the SMBH mass in the $10^6-10^9{\rm ~M}_\odot$ mass range, leading to a merger rate $\Gamma \simeq 3-8$ yr$^{-1}$ Gpc$^{-3}$ at redshift zero. Nearly $40\%$ of the mergers have masses in the "BH mass gap", $50-140{\rm ~M}_\odot$, thus indicating that galactic nuclei are ideal places to form BHs in this mass range. We argue that gravitational wave (GW) sources with components mass $m_1>40{\rm ~M}_\odot$ and $m_2<30{\rm ~M}_\odot$ would represent a strong indicator of a galactic nuclei origin. The majority of these mergers could be multiband GW sources in the local Universe: nearly $40\%$ might be seen by LISA as eccentric sources and, a few years later, as circular sources by LIGO and the Einstein Telescope, making decihertz observatories like DECIGO unique instruments to bridge the observations during the binary inspiral.

We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$).

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

Our Galaxy and the nearby Andromeda Galaxy (M31) form a bound system, even though the relative velocity vector of M31 is currently not well constrained. Their orbital motion is highly dependent on the initial conditions, but all the reliable scenarios imply a first close approach in the next 3-5 Gyrs. In our study, we simulate this interaction via direct N-body integration, using the HiGPUs code. Our aim is to investigate the dependence of the time of the merger on the physical and dynamical properties of the system. Finally, we study the dynamical evolution of the two Supermassive Black Holes placed in the two galactic centers, with the future aim to achieve a proper resolution to follow their motion until they form a tight binary system.

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.

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.

According to our current knowledge about physical and dynamical properties of the Milky Way-M31 system, it seems likely that these two galaxies will collide and eventually merge in a time very sensitive to initial conditions. Using the HiGPUs code, we have performed several numerical simulations to study the dynamics of the system, trying to define the role of indeterminacy in the present day observed relative velocities of the two galaxies and the time of the merger. At the same time, we have followed the dynamics of the two massive black holes sitting in the galactic centers, to check (within the space and time resolution limits of our simulation) their relative motion upon the completion of the galaxies merger process.

In this paper we explore the interplay between intermediate-mass black holes (IMBHs) and their nursing globular clusters (GCs), taking advantage of over 2000 Monte Carlo GC models. We find that the average density of IMBHs sphere of influence can be uniquely connected to the host GCs luminosity and half-light radius via a fundamental plane. We propose a statistical approach to systematically identify potential Galactic GCs harbouring either an IMBH or a massive subsystem comprised of stellar BHs. Our models show that the IMBH is often bound to a stellar companion or a stellar BH, which can lead to tidal disruption events or to low-frequency gravitational waves. We show that GCs orbiting close to the Galactic Centre have a larger probability to witness IMBH formation during their early evolution. These low-orbit GCs can deliver several IMBHs into the galaxy innermost regions, with potential impact on both electromagnetic and GW emission. We discuss potential connections between IMBHs and SMBHs inhabiting galactic nuclei, exploring the possibility that in some cases they share similar formation pathways.

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.

A large fraction of galactic nuclei is expected to host supermassive black hole binaries (BHB), likely formed during the early phase of galaxies assembly and merging. In this paper, we use a large set of state-of-art numerical models to investigate the interplay between a BHB and a massive star cluster (GCs) driven toward the galactic centre by dynamical friction. Varying the BHB mass and mass ratio and the GC orbit, we show that the reciprocal feedback exerted between GCs and the BHB shapes their global properties. We show that, at GC-to-BHB mass ratios above 0.1, the GC affects notably the BHB orbital evolution, possibly boosting its coalescence. This effect is maximized if the GC moves on a retrograde orbit, and for a non-equal mass BHB. We show that the GC debris dispersed around the galactic nucleus can lead to the formation of a nuclear cluster, depending on the BHB tidal field, and that the distribution of compact remnants resulting from the GC disruption can carry information about the BHB orbital properties. We find that red giant stars delivered by the spiralling GC can be disrupted at a rate of $\simeq (0.7-7)\times 10^{-7}$ yr$^{-1}$ for BHB masses $\sim 10^7{\rm M}_\odot$, while tens to hundreds of stars can be possibly observed in the galactic halo as high-velocity stars, with velocities up to $\sim 2000$ km s$^{-1}$, depending on the BHB orbital properties.

In this paper we propose a novel technique to constrain the progenitor binary black hole (BBH) formation history using the remnant masses and spins of merged black holes (BHs). Exploring different models, we found that dynamically formed BBHs are distributed differently in the mass-spin plane than those that have formed in isolation. Stellar evolution recipes crucially affect the remnant mass distribution, suggesting that future efforts should be devoted to finding a common way of modelling the evolutionary phases of single and binary stars. Our simple approach has allowed us to place weak constraints on the origin of the presently observed population of merged BHs, although with high uncertainties. Our results show that the fingerprints of different BBH formation channels will emerge as soon as LIGO detects more than $\sim 10^2$ merger events. This work provides a way of thinking that can be easily used both by people working on isolated and on dynamical BBH formation and evolution.

There have been increasing theoretical speculations and observational indications that certain globular clusters (GCs) could contain a sizeable population of stellar mass black holes (BHs). In this paper, we shortlist at least 29 Galactic GCs that could be hosting a subsystem of BHs (BHS). In a companion paper, we analysed results from a wide array of GC models (simulated with the MOCCA code for cluster simulations) that retained few tens to several hundreds of BHs at 12 Gyr and showed that the properties of the BHS in those GCs correlate with the GC's observable properties. Building on those results, we use available observational properties of 140 Galactic GCs to identify 29 GCs that could potentially be harbouring up to a few hundreds of BHs. Utilizing observational properties and theoretical scaling relations, we estimate the density, size and mass of the BHS in these GCs. We also calculate the total number of BHs and the fraction of BHs contained in a binary system for our shortlisted Galactic GCs. Additionally, we mention other Galactic GCs that could also contain significant number of single BHs or BHs in binary systems.

We model the inspiral of globular clusters (GCs) towards a galactic nucleus harboring a supermassive black hole (SMBH), a leading scenario for the formation of nuclear star clusters. We consider the case of GCs containing either an intermediate-mass black hole (IMBH) or a population of stellar mass black holes (BHs), and study the formation of gravitational wave (GW) sources. We perform direct summation $N$-body simulations of the infall of GCs with different orbital eccentricities in the live background of a galaxy with either a shallow or steep density profile. We find that the GC acts as an efficient carrier for the IMBH, facilitating the formation of a bound pair. The hardening and evolution of the binary depends sensitively on the galaxy's density profile. If the host galaxy has a shallow profile the hardening is too slow to allow for coalescence within a Hubble time, unless the initial cluster orbit is highly eccentric. If the galaxy hosts a nuclear star cluster, the hardening leads to coalescence by emission of GWs within $3-4$ Gyr. In this case, we find a IMBH-SMBH merger rate of $\Gamma_{\rm IMBH-SMBH} = 2.8\times 10^{-3}$ yr$^{-1}$ Gpc$^{-3}$. If the GC hosts a population of stellar BHs, these are deposited close enough to the SMBH to form extreme-mass-ratio-inspirals with a merger rate of $\Gamma_{\rm EMRI} = 0.25$ yr$^{-1}$ Gpc$^{-3}$. Finally, the SMBH tidal field can boost the coalescence of stellar black hole binaries delivered from the infalling GCs. The merger rate for this merging channel is $\Gamma_{\rm BHB} = 0.4-4$ yr$^{-1}$ Gpc$^{-3}$.

In this paper, we discuss how globular clusters (GCs) structural and observational properties can be used to infer the presence of a black hole system (BHS) inhabiting their inner regions. We propose a novel way to identify the BHS size, defined as the GC radius containing a mass contributed equally from stars and stellar BHs. Using this definition, similar to the well-known concept of "influence radius", we found a "fundamental plane" connecting the BHS typical density with the GC central surface density profile, total luminosity and observational half-mass radius. Our approach allows us to define a unique way to connect the observational GCs parameters with their dark content. Comparing our results with observed Milky Way GCs, we found that many of them likely host, at the present time, as many as several hundreds of BHs. These BHS are characterized by a relatively low typical density, $\rho_\bhs \sim 10-10^5\Ms$ pc$^{-3}$ and composed of relatively massive BHs, with average masses in the range $m_\bhs = 14-22\Ms$. We also show that a similar approach can be used to find Milky Way GCs potentially hosting an intermediate-mass black hole.

We discuss how the centres of galaxy clusters evolve in time, showing the results of a series of direct N-body simulations. In particular, we followed the evolution of a galaxy cluster with a mass $M_{clus} \simeq 10^{14} $M$_{\odot}$ in different configurations. The dynamical evolution of the system leads in all the cases to the formation of dense and massive sub-structures in the cluster centre, that form in consequence of a series of collisions and merging among galaxies travelling in the cluster core. We investigate how the structural properties of the main merging product depends on the characteristics of those galaxies that contributed to its formation.

The application of high end computing to astrophysical problems, mainly in the galactic environment, is under development since many years at the Dep. of Physics of Sapienza Univ. of Roma. The main scientific topic is the physics of self gravitating systems, whose specific subtopics are: i) celestial mechanics and interplanetary probe transfers in the solar system; ii) dynamics of globular clusters and of globular cluster systems in their parent galaxies; iii) nuclear clusters formation and evolution; iv) massive black hole formation and evolution; v) young star cluster early evolution. In this poster we describe the software and hardware computational resources available in our group and how we are developing both software and hardware to reach the scientific aims above itemized.