Pair instability (PI) and pulsational PI prevent the formation of black holes (BHs) with mass $\gtrsim{}60$ M$_\odot$ from single star evolution. Here, we investigate the possibility that BHs with mass in the PI gap form via stellar mergers and multiple stellar mergers, facilitated by dynamical encounters in young star clusters. We analyze $10^4$ simulations, run with the direct N-body code nbody6++gpu coupled with the population synthesis code MOBSE. We find that up to $\sim{}6$~% of all simulated BHs have mass in the PI gap, depending on progenitor's metallicity. This formation channel is strongly suppressed in metal-rich (Z = 0.02) star clusters, because of stellar winds. BHs with mass in the PI gap are initially single BHs but can efficiently acquire companions through dynamical exchanges. We find that $\sim{}$21%, 10% and 0.5% of all binary BHs have at least one component in the PI mass gap at metallicity Z = 0.0002, 0.002 and 0.02, respectively. Based on the evolution of the cosmic star formation rate and metallicity, and under the assumption that all stars form in young star clusters, we predict that $\sim{}5$~% of all binary BH mergers detectable by advanced LIGO and Virgo at their design sensitivity have at least one component in the PI mass gap.

Young star clusters are likely the most common birthplace of massive stars across cosmic time and influence the formation of compact binaries in several ways. Here, we simulate the formation of black hole -- neutron star binaries (BHNSs) in young star clusters, by means of the binary population synthesis code \texttt{MOBSE} interfaced with the $N$-body code \texttt{NBODY6++GPU}. BHNSs formed in young star clusters (dynamical BHNSs) are significantly more massive than BHNSs formed from isolated binaries (isolated BHNSs): $\sim{}40$~\% of the dynamical BHNS mergers have total mass $>15$ M$_\odot$, while only $\sim{}0.01$~\% of the isolated BHNS mergers have mass in excess of this value. Hence, our models strongly support a dynamical formation scenario for GW190814, given its total mass $\sim{}26$ M$_\odot$, if this event is a BHNS merger. All our dynamical BHNSs are ejected from their parent star cluster before they reach coalescence. Thus, a significant fraction of BHNS mergers occurring in the field might have originated in a young star cluster. The mass spectrum of BHNS mergers from gravitational-wave detections will provide a clue to differentiate between dynamical and isolated formation of BHNSs.

We investigate the properties of the host galaxies of compact binary mergers across cosmic time. To this end, we combine population synthesis simulations together with galaxy catalogues from the hydrodynamical cosmological simulation EAGLE to derive the properties of the host galaxies of binary neutron star (BNS), black hole-neutron star (BHNS) and binary black hole (BBH) mergers. Within this framework, we derive the host galaxy probability, i.e., the probability that a galaxy hosts a compact binary coalescence as a function of its stellar mass, star formation rate, $K_s$ magnitude and $B$ magnitude. This quantity is particularly important for low-latency searches of gravitational wave (GW) sources as it provides a way to rank galaxies lying inside the credible region in the sky of a given GW detection, hence reducing the number of viable host candidates. Furthermore, even if no electromagnetic counterpart is detected, the proposed ranking criterion can still be used to classify the galaxies contained in the error box. Our results show that massive galaxies (or equivalently galaxies with a high luminosity in $K_s$ band) have a higher probability of hosting BNS, BHNS, and BBH mergers. We provide the probabilities in a suitable format to be implemented in future low-latency searches.

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

Next generation ground-based gravitational-wave detectors will observe binary black hole (BBH) mergers up to redshift $\gtrsim{}10$, probing the evolution of compact binary (CB) mergers across cosmic time. Here, we present a new data-driven model to estimate the cosmic merger rate density (MRD) evolution of CBs, by coupling catalogs of CB mergers with observational constraints on the cosmic star formation rate density and on the metallicity evolution of the Universe. We adopt catalogs of CB mergers derived from recent $N-$body and population-synthesis simulations, to describe the MRD of CBs formed in young star clusters (hereafter, dynamical CBs) and in the field (hereafter, isolated CBs). The local MRD of dynamical BBHs is $\mathcal{R}_{\rm BBH}=67^{+29}_{-23}$ Gpc$^{-3}$ yr$^{-1}$, consistent with the 90% credible interval from the first and second observing run (O1 and O2) of the LIGO-Virgo collaboration, and with the local MRD of isolated BBHs ($\mathcal{R}_{\rm BBH}=49^{+79}_{-37}$ Gpc$^{-3}$ yr$^{-1}$). The local MRD of dynamical and isolated black hole-neutron star binaries (BHNSs) is $\mathcal{R}_{\rm BHNS}=38^{+32}_{-24}$ and $45^{+45}_{-32}$ Gpc$^{-3}$ yr$^{-1}$, respectively. Both values are consistent with the upper limit inferred from O1 and O2. Finally, the local MRD of dynamical binary neutron stars (BNSs, $\mathcal{R}_{\rm BNS}=150^{+56}_{-40}$ Gpc$^{-3}$ yr$^{-1}$) is a factor of two lower than the local MRD of BNSs formed in isolated binaries ($\mathcal{R}_{\rm BNS}=281^{+109}_{-75}$ Gpc$^{-3}$ yr$^{-1}$). The MRD for all CB classes grows with redshift, reaching its maximum at $z \in [1.5,2.5]$, and then decreases. This trend springs from the interplay between cosmic star formation rate, metallicity evolution and delay time of binary compact objects.

Characterizing the properties of the host galaxies of merging compact objects provides essential clues to interpret current and future gravitational-wave detections. Here, we investigate the stellar mass, star formation rate (SFR), metallicity and colours of the host galaxies of merging compact objects in the local Universe, by combining the results of MOBSE population-synthesis models together with galaxy catalogs from the EAGLE simulation. We predict that the stellar mass of the host galaxy is an excellent tracer of the merger rate per galaxy ${\rm n}_{\rm GW}$ of double neutron stars (DNSs), double black holes (DBHs) and black hole-neutron star binaries (BHNSs). We find a significant correlation also between ${\rm n}_{\rm GW}$ and SFR. As a consequence, ${\rm n}_{\rm GW}$ correlates also with the $r-$band luminosity and with the $g-r$ colour of the host galaxies. Interestingly, $\gtrsim{}60$ %, $\gtrsim{}64$ % and $\gtrsim{}73$ % of all the DNSs, BHNSs and DBHs merging in the local Universe lie in early-type galaxies, such as NGC 4993. We predict a local DNS merger rate density of $\sim{}238~{\rm Gpc}^{-3}~{\rm yr}~^{-1}$ and a DNS merger rate $\sim{}16-121$ Myr$^{-1}$ for Milky Way-like galaxies. Thus, our results are consistent with both the DNS merger rate inferred from GW170817 and the one inferred from Galactic DNSs.

The next generation ground-based gravitational wave interferometers will possibly observe mergers of binary black holes (BBHs) and binary neutron stars (BNSs) to redshift $z\gtrsim{}10$ and $z\gtrsim{}2$, respectively. Here, we characterize the properties of merging BBHs, BNSs and neutron star-black hole binaries across cosmic time, by means of population-synthesis simulations combined with the Illustris cosmological simulation. We find that the mass of merging compact objects does not depend (or depends very mildly) on the merger redshift. Even the mass distribution of black holes depends only mildly on redshift, because BBHs originating from metal-poor progenitors ($Z\leq{}4\times{}10^{-3}$) dominate the entire population of merging BBHs across cosmic time. For a common-envelope efficiency $\alpha{}\ge{}3$, the main difference between the mass distribution of BBHs merging in the last Gyr and that of BBHs merging more than 11 Gyr ago is that there is an excess of heavy merging black holes ($20-35$ M$_\odot$) in the last Gyr. This excess is explained by the longer delay time of massive BBHs.