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.

Ten binary black-hole mergers have already been detected during the first two observing runs of advanced LIGO and Virgo, and many more are expected to be observed in the near future. This opens the possibility for gravitational-wave astronomy to better constrain the properties of black hole binaries, not only as single sources, but as a whole astrophysical population. In this paper, we address the problem of using gravitational-wave measurements to estimate the proportion of merging black holes produced either via isolated binaries or binaries evolving in young star clusters. To this end, we use a Bayesian hierarchical modeling approach applied to catalogs of merging binary black holes generated using state-of-the-art population synthesis and N-body codes. In particular, we show that, although current advanced LIGO/Virgo observations only mildly constrain the mixing fraction $f \in [0,1]$ between the two formation channels, we expect to narrow down the fractional errors on $f$ to $10-20\%$ after a few hundreds of detections.

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.

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 existence of massive stellar black hole binaries (MBHBs), with primary black hole (BH) masses $\ge 31 \, M_\odot$, was proven by the detection of the gravitational wave (GW) event GW150914 during the first LIGO/Virgo observing run (O1), and successively confirmed by seven additional GW signals discovered in the O1 and O2 data. By adopting the galaxy formation model \texttt{GAMESH} coupled with binary population synthesis (BPS) calculations, here we investigate the origin of these MBHBs by selecting simulated binaries compatible in mass and coalescence redshifts. We find that their cosmic birth rates peak in the redshift range $6.5 \leq z \leq 10$, regardless of the adopted BPS. These MBHBs are then old systems forming in low-metallicity ($Z \sim [0.01-0.1] \, Z_{\odot}$), low-stellar-mass galaxies, before the end of cosmic reionization, i.e. significantly beyond the peak of cosmic star formation. GW signals generated by coalescing MBHBs open up new possibilities to probe the nature of stellar populations in remote galaxies, at present too faint to be detected by available electromagnetic facilities.

Natal kicks are matter of debate and significantly affect the merger rate density of compact objects. Here, we present a new simple formalism for natal kicks of neutron stars (NSs) and black holes (BHs). We describe the magnitude of the kick as $v_{\rm kick}\propto{}f_{\rm H05}\,{}\,{}m_{\rm ej}\,{}\,{}m_{\rm rem}^{-1}$, where $f_{\rm H05}$ is a normalization factor, drawn from a Maxwellian distribution with one-dimensional root-mean-square velocity $\sigma{}=265$~km~s$^{-1}$, $m_{\rm ej}$ is the mass of the supernova (SN) ejecta and $m_{\rm rem}$ is the mass of the compact object. This formalism matches the proper motions of young Galactic pulsars and can naturally account for the differences between core-collapse SNe of single stars, electron-capture SNe and ultra-stripped SNe occurring in interacting binaries. Finally, we use our new kick formalism to estimate the local merger rate density of binary NSs ($R_{\rm BNS}$), BH--NS binaries ($R_{\rm BHNS}$) and binary BHs ($R_{\rm BBH}$), based on the cosmic star formation rate density and metallicity evolution. In our fiducial model, we find $R_{\rm BNS}\sim{}600$~Gpc$^{-3}$~yr$^{-1}$, $R_{\rm BHNS}\sim{}10$~Gpc$^{-3}$~yr$^{-1}$ and $R_{\rm BBH}\sim{}50$~Gpc$^{-3}$~yr$^{-1}$, fairly consistent with the numbers inferred from the LIGO-Virgo collaboration.

We investigate the impact of stellar rotation on the formation of black holes (BHs), by means of our population-synthesis code SEVN. Rotation affects the mass function of BHs in several ways. In massive metal-poor stars, fast rotation reduces the minimum zero-age main sequence (ZAMS) mass for a star to undergo pair instability and pulsational pair instability. Moreover, stellar winds are enhanced by rotation, peeling-off the entire hydrogen envelope. As a consequence of these two effects, the maximum BH mass we expect from the collapse a rotating metal-poor star is only $\sim{}45$ M$_\odot$, while the maximum mass of a BH born from a non-rotating star is $\sim{}60$ M$_\odot$. Furthermore, stellar rotation reduces the minimum ZAMS mass for a star to collapse into a BH from $\sim{}18-25$ M$_\odot$ to $\sim{}13-18$ M$_\odot$. Finally, we have investigated the impact of different core-collapse supernova (CCSN) prescriptions on our results. While the threshold value of compactness for direct collapse and the fallback efficiency strongly affect the minimum ZAMS mass for a star to collapse into a BH, the fraction of hydrogen envelope that can be accreted onto the final BH is the most important ingredient to determine the maximum BH mass. Our results confirm that the interplay between stellar rotation, CCSNe and pair instability plays a major role in shaping the BH mass spectrum.

Observations of gravitational waves and their electromagnetic counterparts may soon uncover the existence of coalescing compact binary systems formed by a stellar-mass black hole and a neutron star. These mergers result in a remnant black hole, possibly surrounded by an accretion disk. The mass and spin of the remnant black hole depend on the properties of the coalescing binary. We construct a map from the binary components to the remnant black hole using a sample of numerical-relativity simulations of different mass ratios $q$, (anti-)aligned dimensionless spins of the black hole $a_{\rm BH}$, and several neutron star equations of state. Given the binary total mass, the mass and spin of the remnant black hole can therefore be determined from the three parameters $(q,a_{\rm BH},\Lambda)$, where $\Lambda$ is the tidal deformability of the neutron star. Our models also incorporate the binary black hole and test-mass limit cases and we discuss a simple extension for generic black hole spins. We combine the remnant characterization with recent population synthesis simulations for various metallicities of the progenitor stars that generated the binary system. We predict that black-hole-neutron-star mergers produce a population of remnant black holes with masses distributed around $7M_\odot$ and $9M_\odot$. For isotropic spin distributions, nonmassive accretion disks are favoured: no bright electromagnetic counterparts are expected in such mergers.

We investigate the properties of the host galaxies of compact binary mergers across cosmic time, by means of population-synthesis simulations combined with galaxy catalogues from the EAGLE suite. We analyze the merger rate per galaxy of binary neutron stars (BNSs), black hole--neutron star binaries (BHNSs) and binary black holes (BBHs) from redshift zero up to redshift six. The binary merger rate per galaxy strongly correlates with the stellar mass of the host galaxy at any redshift considered here. This correlation is significantly steeper for BNSs than for both BHNSs and BBHs. Moreover, we find that the merger rate per galaxy depends also on host galaxy's star formation rate and metallicity. We derive a robust fitting formula that relates the merger rate per galaxy with galaxy's star formation rate, stellar mass and metallicity at different redshifts. The typical masses of the host galaxies increase significantly as redshift decreases, as a consequence of the interplay between delay time distribution of compact binaries and cosmic assembly of galaxies. Finally, we study the evolution of the merger rate density with redshift. At low redshift ($z\leq{}0.1$) early-type galaxies give a larger contribution to the merger rate density than late-type galaxies. This trend reverts at $z\ge{}1$.

Using simulations performed with the population synthesis code MOBSE, we compute the merger rate densities and detection rates of compact binary mergers formed in isolation for second- and third-generation gravitational-wave detectors. We estimate how rates are affected by uncertainties on key stellar-physics parameters, namely common envelope evolution and natal kicks. We estimate how future upgrades will increase the size of the available catalog of merger events, and we discuss features of the merger rate density that will become accessible with third-generation detectors.

We explore the host galaxies of compact-object binaries (black hole--black hole binaries, BHBs; neutron star--black hole binaries, NSBHs; double-neutron stars; DNSs) across cosmic time, by means of population-synthesis simulations combined with the Illustris cosmological simulation. At high redshift ($z\gtrsim{}4$) the host galaxies of BHBs, NSBHs and DNSs are very similar and are predominantly low-mass galaxies (stellar mass $M<10^{11}$ M$_\odot$). If $z\gtrsim{}4$ most compact objects form and merge in the same galaxy, with a short delay time. At low redshift ($z\leq{}2$), the host galaxy populations of DNSs differ significantly from the host galaxies of both BHBs and NSBHs. DNSs merging at low redshift tend to form and merge in the same galaxy, with relatively short delay time. The stellar mass of DNS hosts peaks around $\sim{}10^{10}-10^{11}$ M$_\odot$. In contrast, BHBs and NSBHs merging at low redshift tend to form in rather small galaxies at high redshift and then to merge in larger galaxies with long delay times. This difference between DNSs and black hole binaries is a consequence of their profoundly different metallicity dependence.

Six gravitational wave events have been reported by the LIGO-Virgo collaboration (LVC), five of them associated with black hole binary (BHB) mergers and one with a double neutron star (DNS) merger, while the coalescence of a black hole-neutron star (BHNS) binary is still missing. We investigate the progenitors of double compact object binaries with our population-synthesis code MOBSE. MOBSE includes advanced prescriptions for mass loss by stellar winds (depending on metallicity and on the Eddington ratio) and a formalism for core-collapse, electron-capture and (pulsational) pair instability supernovae. We investigate the impact of progenitor's metallicity, of the common-envelope parameter $\alpha{}$ and of the natal kicks on the properties of DNSs, BHNSs and BHBs. We find that neutron-star (NS) masses in DNSs span from 1.1 to 2.0 M$_\odot$, with a preference for light NSs, while NSs in merging BHNSs have mostly large masses ($1.3-2.0$ M$_\odot$). BHs in merging BHNSs are preferentially low mass ($5-15$ M$_\odot$). BH masses in merging BHBs strongly depend on the progenitor's metallicity and span from $\sim{}5$ to $\sim{}45$ M$_\odot$. The local merger rate density of both BHNSs and BHBs derived from our simulations is consistent with the values reported by the LVC in all our simulations. In contrast, the local merger rate density of DNSs matches the value inferred from the LVC only if low natal kicks are assumed. This result adds another piece to the intricate puzzle of natal kicks and DNS formation.

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.

Studying the formation and evolution of black hole binaries (BHBs) is essential for the interpretation of current and forthcoming gravitational wave (GW) detections. We investigate the statistics of BHBs that form from isolated binaries, by means of a new version of the SEVN population-synthesis code. SEVN integrates stellar evolution by interpolation over a grid of stellar evolution tracks. We upgraded SEVN to include binary stellar evolution processes and we used it to evolve a sample of $1.5\times{}10^8$ binary systems, with metallicity in the range $\left[10^{-4};4\times 10^{-2}\right]$. From our simulations, we find that the mass distribution of black holes (BHs) in double compact-object binaries is remarkably similar to the one obtained considering only single stellar evolution. The maximum BH mass we obtain is $\sim 30$, $45$ and $55\, \mathrm{M}_\odot$ at metallicity $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. A few massive single BHs may also form ($\lesssim 0.1\%$ of the total number of BHs), with mass up to $\sim 65$, $90$ and $145\, \mathrm{M}_\odot$ at $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. These BHs fall in the mass gap predicted from pair-instability supernovae. We also show that the most massive BHBs are unlikely to merge within a Hubble time. In our simulations, merging BHs like GW151226 and GW170608, form at all metallicities, the high-mass systems (like GW150914, GW170814 and GW170104) originate from metal poor ($Z\lesssim{}6\times 10^{-3}$) progenitors, whereas GW170729-like systems are hard to form, even at $Z = 10^{-4}$. The BHB merger rate in the local Universe obtained from our simulations is $\sim 90 \mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$, consistent with the rate inferred from LIGO-Virgo data.

Searching for distinctive signatures, which characterize different formation channels of binary black holes (BBHs), is a crucial step towards the interpretation of current and future gravitational wave detections. Here, we investigate the demography of merging BBHs in young star clusters (SCs), which are the nursery of massive stars. We performed $4\times{} 10^3$ N-body simulations of SCs with metallicity $Z=0.002$, initial binary fraction $0.4$ and fractal initial conditions, to mimic the clumpiness of star forming regions. Our simulations include a novel population-synthesis approach based on the code MOBSE. We find that SC dynamics does not affect the merger rate significantly, but leaves a strong fingerprint on the properties of merging BBHs. More than 50 % of merging BBHs in young SCs form by dynamical exchanges in the first few Myr. Dynamically formed merging BBHs are significantly heavier than merging BBHs in isolated binaries: merging BBHs with total mass up to $\sim{}120$ M$_\odot$ form in young SCs, while the maximum total mass of merging BBHs in isolated binaries with the same metallicity is only $\sim{}70$ M$_\odot$. Merging BBHs born via dynamical exchanges tend to have smaller mass ratios than BBHs in isolated binaries. Furthermore, SC dynamics speeds up the merger: the delay time between star formation and coalescence is significantly shorter in young SCs. In our simulations, massive systems such as GW170729 form only via dynamical exchanges. Finally $\sim{}2$ % of merging BBHs in young SCs have mass in the pair-instability mass gap ($\sim{}60-120$ M$_\odot$). This represents a unique fingerprint of merging BBHs in SCs.

Natal kicks are one of the most debated issues about double neutron star (DNS) formation. Several observational and theoretical results suggest that some DNSs have formed with low natal kicks ($\lesssim{}50$ km s$^{-1}$), which might be attributed to electron-capture supernovae (ECSNe). We investigate the impact of ECSNe on the formation of DNSs by means of population synthesis simulations. In particular, we assume a Maxwellian velocity distribution for the natal kick induced by ECSNe with one dimensional root-mean-square $\sigma_{\rm ECSN} = 0,7,15,26,265$ km s$^{-1}$. The total number of DNSs scales inversely with $\sigma_{\rm ECSN}$ and the number of DNS mergers is higher for relatively low kicks. This effect is particularly strong if we assume low efficiency of common-envelope ejection (described by the parameter $\alpha=1$), while it is only mild for high efficiency of common-envelope ejection ($\alpha{}=5$). In most simulations, more than 50 per cent of the progenitors of merging DNSs undergo at least one ECSN and the ECSN is almost always the first SN occurring in the binary system. Finally, we have considered the extreme case in which all neutron stars receive a low natal kick ($\lesssim{}50$~km~s$^{-1}$). In this case, the number of DNSs increases by a factor of ten and the percentage of merging DNSs which went through an ECSN is significantly suppressed ($<40$ per cent).

The high mass X-ray binaries (HMXBs) provide an exciting framework to investigate the evolution of massive stars and the processes behind binary evolution. HMXBs have shown to be good tracers of recent star formation in galaxies and might be important feedback sources at early stages of the Universe. Furthermore, HMXBs are likely the progenitors of gravitational wave sources (BH--BH or BH--NS binaries that may merge producing gravitational waves). In this work, we investigate the nature and properties of HMXB population in star-forming galaxies. We combine the results from the population synthesis model MOBSE (Giacobbo et al. 2018) together with galaxy catalogs from EAGLE simulation (Schaye et al. 2015). Therefore, this method describes the HMXBs within their host galaxies in a self-consistent way. We compute the X-ray luminosity function (XLF) of HMXBs in star-forming galaxies, showing that this methodology matches the main features of the observed XLF.

We investigate the host galaxies of compact objects merging in the local Universe, by combining the results of binary population-synthesis simulations with the Illustris cosmological box. Double neutron stars (DNSs) merging in the local Universe tend to form in massive galaxies (with stellar mass $>10^{9}$ M$_\odot$) and to merge in the same galaxy where they formed, with a short delay time between the formation of the progenitor stars and the DNS merger. In contrast, double black holes (DBHs) and black hole $-$ neutron star binaries (BHNSs) form preferentially in small galaxies (with stellar mass $<10^{10}$ M$_\odot$) and merge either in small or in larger galaxies, with a long delay time. This result is an effect of metallicity: merging DBHs and BHNSs form preferentially from metal-poor progenitors ($Z\leq{}0.1$ Z$_\odot$), which are more common in high-redshift galaxies and in local dwarf galaxies, whereas merging DNSs are only mildly sensitive to progenitor's metallicity and thus are more abundant in massive galaxies nowadays. The mass range of DNS hosts we predict in this work is consistent with the mass range of short gamma-ray burst hosts.

The recent detection of gravitational waves has proven the existence of massive stellar black hole binaries (BBHs), but the formation channels of BBHs are still an open question. Here, we investigate the demography of BBHs by using our new population-synthesis code MOBSE. MOBSE is an updated version of the widely used binary population-synthesis code BSE (Hurley et al. 2000, Hurley et al. 2002) and includes the key ingredients to determine the fate of massive stars: up-to-date stellar wind prescriptions and supernova models. With MOBSE, we form BBHs with total mass up to $\sim{}120$ M$_\odot$ at low metallicity, but only systems with total mass up to $\sim{}80$ M$_\odot$ merge in less than a Hubble time. Our results show that only massive metal-poor stars ($Z\lesssim 0.002$) can be the progenitors of gravitational wave events like GW150914. Moreover, we predict that merging BBHs form much more efficiently from metal-poor than from metal-rich stars.

Six gravitational wave detections have been reported so far, providing crucial insights on the merger rate of double compact objects. We investigate the cosmic merger rate of double neutron stars (DNSs), neutron star-black hole binaries (NSBHs) and black hole binaries (BHBs) by means of population-synthesis simulations coupled with the Illustris cosmological simulation. We have performed six different simulations, considering different assumptions for the efficiency of common envelope (CE) ejection and exploring two distributions for the supernova (SN) kicks. The current BHB merger rate derived from our simulations spans from $\sim{}150$ to $\sim{}240$ Gpc$^{-3}$ yr$^{-1}$ and is only mildly dependent on CE efficiency. In contrast, the current merger rates of DNSs (ranging from $\sim{}20$ to $\sim{}600$ Gpc$^{-3}$ yr$^{-1}$) and NSBHs (ranging from $\sim{}10$ to $\sim{}100$ Gpc$^{-3}$ yr$^{-1}$) strongly depend on the assumptions on CE and natal kicks. The merger rate of DNSs is consistent with the one inferred from the detection of GW170817 only if a high efficiency of CE ejection and low SN kicks (drawn from a Maxwellian distribution with one dimensional root mean square $\sigma{}=15$ km s$^{-1}$) are assumed.

The first four gravitational wave events detected by LIGO were all interpreted as merging black hole binaries (BHBs), opening a new perspective on the study of such systems. Here we use our new population-synthesis code MOBSE, an upgraded version of BSE (Hurley et al. 2002), to investigate the demography of merging BHBs. MOBSE includes metallicity-dependent prescriptions for mass loss of massive hot stars. It also accounts for the impact of the electron-scattering Eddington factor on mass loss. We perform >10^8 simulations of isolated massive binaries, with 12 different metallicities, to study the impact of mass loss, core-collapse supernovae and common envelope on merging BHBs. Accounting for the dependence of stellar winds on the Eddington factor leads to the formation of black holes (BHs) with mass up to 65 Msun at metallicity Z~0.0002. However, most BHs in merging BHBs have masses <40 Msun. We find merging BHBs with mass ratios in the 0.1 - 1.0 range, even if mass ratios >0.6 are more likely. We predict that systems like GW150914, GW170814 and GW170104 can form only from progenitors with metallicity Z<=0.006, Z<=0.008 and Z<=0.012, respectively. Most merging BHBs have gone through a common envelope phase, but up to ~17 per cent merging BHBs at low metallicity did not undergo any common envelope phase. We find a much higher number of mergers from metal-poor progenitors than from metal-rich ones: the number of BHB mergers per unit mass is ~10^-4 Msun^-1 at low metallicity (Z = 0.0002 - 0.002) and drops to ~10^-7 Msun^-1 at high metallicity (Z ~ 0.02).

The cosmic merger rate density of black hole binaries (BHBs) can give us an essential clue to constraining the formation channels of BHBs, in light of current and forthcoming gravitational wave detections. Following a Monte Carlo approach, we couple new population-synthesis models of BHBs with the Illustris cosmological simulation, to study the cosmic history of BHB mergers. We explore six population-synthesis models, varying the prescriptions for supernovae, common envelope, and natal kicks. In most considered models, the cosmic BHB merger rate follows the same trend as the cosmic star formation rate. The normalization of the cosmic BHB merger rate strongly depends on the treatment of common envelope and on the distribution of natal kicks. We find that most BHBs merging within LIGO's instrumental horizon come from relatively metal-poor progenitors (<0.2 Zsun). The total masses of merging BHBs span a large range of values, from ~6 to ~82 Msun. In our fiducial model, merging BHBs consistent with GW150914, GW151226 and GW170104 represent ~6, 3, and 12 per cent of all BHBs merging within the LIGO horizon, respectively. The heavy systems, like GW150914, come from metal-poor progenitors (<0.15 Zsun). Most GW150914-like systems merging in the local Universe appear to have formed at high redshift, with a long delay time. In contrast, GW151226-like systems form and merge all the way through the cosmic history, from progenitors with a broad range of metallicities. Future detections will be crucial to put constraints on common envelope, on natal kicks, and on the BHB mass function.

The mass spectrum of stellar black holes (BHs) is highly uncertain. Theoretical models of BH formation strongly depend on the efficiency of stellar winds of the progenitor star and on the supernova (SN) explosion mechanism. We discuss the BH mass spectrum we obtain using SEVN, a new public population-synthesis code that includes up-to-date stellar-wind prescriptions and several SN explosion models. Our models indicate a sub-solar metallicity environment for the progenitors of the gravitational wave source GW150914. We show that our models predict substantially larger BH masses (up to ~100 Msun) than other population synthesis codes, at low metallicity. In this proceeding, we also discuss the impact of pair-instability SNe on our previously published models.