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.

Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called "upper-mass gap" would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through pre-collapse stellar collisions in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as $\sim20\%$ of all BH progenitors undergo one or more collisions prior to stellar collapse and up to $\sim1\%$ of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the "multiple-generation" merger channel described in other analyses. Finally, we describe how stellar collisions in clusters may provide a unique pathway to pair-instability supernovae and briefly discuss the expected rate of these events and other electromagnetic transients.

More and more observations indicate that young star clusters could retain imprints of their formation process. In particular, the degree of substructuring and rotation are possibly the direct result of the collapse of the parent molecular cloud from which these systems form. Such properties can, in principle, be washed-out, but they are also expected to have an impact on the relaxation of these systems. We ran and analyzed a set of ten hydrodynamical simulations of the formation of embedded star clusters through the collapse of turbulent massive molecular clouds. We systematically studied the fractality of our star clusters, showing that they are all extremely substructured (fractal dimension $D=1.0-1.8$). We also found that fractality is slowly reduced, with time, on small scales, while it persists on large scales on longer timescales. Signatures of rotation are found in different simulations at every time of the evolution, even for slightly supervirial substructures, proving that the parent molecular gas transfers part of its angular momentum to the new stellar systems.

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.

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.

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.

Recent interpretation of the color$-$magnitude diagrams of the Milky Way (MW) bulge has suggested that the observed double red-clump feature can be a natural consequence of He-enhanced stellar populations in the MW bulge. This implies that globular clusters (GCs), where the He-enhanced second-generation (SG) stars can be efficiently created, are the most likely candidate contributors of He-rich stars to the MW bulge. We extend this idea to the Galactic inner halo and investigate the fraction of the SG stars as a function of the Galactocentric distance. We use bluer blue-horizontal branch (bBHB) stars, which are assumed to be originated from He-rich SG populations, as proxies of SG stars, and find that the fraction of bBHB stars increases with decreasing Galactocentric distance. Simulations of the GC evolution in the MW tidal field qualitatively support the observed trend of bBHB enhancement in the inner halo. In these simulations, the increasing tidal force with decreasing Galactocentric distance leads to stripping of stars not only from the outskirts but also from the central regions of GCs, where SG stars are more abundant. We discuss the implication and prospect of our findings concerning the formation history of the bulge and inner halo of the MW.

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.