No Arabic abstract
Hierarchical mergers are one of the distinctive signatures of binary black hole (BBH) formation through dynamical evolution. Here, we present a fast semi-analytic approach to simulate hierarchical mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs). Hierarchical mergers are more common in NSCs than they are in both GCs and YSCs, because of the different escape velocity. The mass distribution of hierarchical BBHs strongly depends on the properties of first-generation BBHs, such as their progenitors metallicity. In our fiducial model, we form black holes (BHs) with masses up to $sim{}10^3$ M$_odot$ in NSCs and up to $sim{}10^2$ M$_odot$ in both GCs and YSCs. When escape velocities in excess of 100 km~s$^{-1}$ are considered, BHs with mass $>10^3$ M$_odot$ are allowed to form in NSCs. Hierarchical mergers lead to the formation of BHs in the pair instability mass gap and intermediate-mass BHs, but only in metal-poor environments. The local BBH merger rate in our models ranges from $sim{}10$ to $sim{} 60$ Gpc$^{-3}$ yr$^{-1}$; hierarchical BBHs in NSCs account for $sim{}10^{-2}- 0.2$ Gpc$^{-3}$ yr$^{-1}$, with a strong upper limit of $sim{}10$ Gpc$^{-3}$ yr$^{-1}$. When comparing our models with the second gravitational-wave transient catalog, we find that multiple formation channels are favored to reproduce the observed BBH population.
The growing population of binary black holes (BBHs) observed by gravitational wave detectors is a potential Rosetta stone for understanding their formation channels. Here, we use an upgraded version of our semi-analytic codes FASTCLUSTER and COSMO$mathcal{R}$ATE to investigate the cosmic evolution of four different BBH populations: isolated BBHs and dynamically formed BBHs in nuclear star clusters (NSCs), globular clusters (GCs), and young star clusters (YSCs). We find that the merger rate density of BBHs in GCs and NSCs is barely affected by stellar metallicity ($Z$), while the rate of isolated BBHs changes wildly with $Z$. BBHs in YSCs behave in an intermediate way between isolated and GC/NSC BBHs. The local merger rate density of Nth-generation black holes (BHs), obtained by summing up hierarchical mergers in GCs, NSCs and YSCs, ranges from $sim{1}$ to $sim{4}$ Gpc$^{-3}$ yr$^{-1}$ and is mostly sensitive to the spin parameter. We find that the mass function of primary BHs evolves with redshift in GCs and NSCs, becoming more top-heavy at higher $z$. In contrast, the primary BH mass function almost does not change with redshift in YSCs and in the field. This signature of the BH mass function has relevant implications for Einstein Telescope and Cosmic Explorer. Finally, our analysis suggests that multiple channels contribute to the BBH population of the second gravitational-wave transient catalogue.
We investigate properties of black hole (BH) binaries formed in globular clusters via dynamical processes, using direct N-body simulations. We pay attention to effects of BH mass function on the total mass and mass ratio distributions of BH binaries ejected from clusters. Firstly, we consider BH populations with two different masses in order to learn basic differences from models with single-mass BHs only. Secondly, we consider continuous BH mass functions adapted from recent studies on massive star evolution in a low metallicity environment, where globular clusters are formed. In this work, we consider only binaries that are formed by three-body processes and ignore stellar evolution and primordial binaries for simplicity. Our results imply that most BH binary mergers take place after they get ejected from the cluster. Also, mass ratios of dynamically formed binaries should be close to one or likely to be less than 2:1. Since the binary formation efficiency is larger for higher-mass BHs, it is likely that a BH mass function sampled by gravitational-wave observations would be weighed toward higher masses than the mass function of single BHs for a dynamically formed population. Applying conservative assumptions regarding globular cluster populations such as small BH mass fraction and no primordial binaries, the merger rate of BH binaries originated from globular clusters is estimated to be at least 6.5 per yr per Gpc^3. Actual rate can be up to more than several times of our conservative estimate.
We derive the probability for a newly formed binary black hole (BBH) to undergo an eccentric gravitational wave (GW) merger during binary-single interactions inside a stellar cluster. By integrating over the hardening interactions such a BBH must undergo before ejection, we find that the observable rate of BBH mergers with eccentricity $>0.1$ at $10 rm{Hz}$ relative to the rate of circular mergers can be as high as $sim 5%$ for a typical globular cluster (GC). This further suggests that BBH mergers forming through GW captures in binary-single interactions, eccentric or not, are likely to constitute $sim 10%$ of the total BBH merger rate from GCs. Such GW capture mergers can only be probed with an $N$-body code that includes General Relativistic corrections, which explains why recent Newtonian cluster studies not have been able to resolve this population. Finally, we show that the relative rate of eccentric BBH mergers depends on the compactness of their host cluster, suggesting that an observed eccentricity distribution can be used to probe the origin of BBH mergers.
Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($chi_mathrm{eff}$), the precession spin parameter ($chi_mathrm{p}$), and the chirp mass ($m_mathrm{chirp}$) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with ${bar chi}_mathrm{typ}equiv overline{(chi_mathrm{eff}^2+chi_mathrm{p}^2)^{1/2}}$, up to roughly the maximum $m_mathrm{chirp}$ among first-generation (1g) BHs, and (ii) a plateau at ${bar chi}_mathrm{typ}sim 0.6$ at higher $m_mathrm{chirp}$. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from ${bar chi}_mathrm{typ}$ as a function of $m_mathrm{chirp}$. The GW data observed in LIGO/Virgo O1--O3a prefers an increase in ${bar chi}_mathrm{typ}$ at low $m_mathrm{chirp}$, which is consistent with the growth of the BH spin magnitude by hierarchical mergers, at $sim 2 sigma$ confidence. A Bayesian analysis suggests that 1g BHs have the maximum mass of $sim 15$--$30,M_odot$ if the majority of mergers are of high-generation BHs (not among 1g-1g BHs), which is consistent with mergers in active galactic nucleus disks and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favored to have non-zero spin magnitudes of $sim 0.3$. We also forecast that signatures for hierarchical mergers in the ${bar chi}_mathrm{typ}$ distribution can be confidently recovered once the number of GW events increases to $gtrsim O(100)$.
Black hole-main sequence star (BH-MS) binaries are one of the targets of the future data releases of the astrometric satellite {it Gaia}. They are supposed to be formed in two main sites: a galactic field and star clusters. However, previous work has never predicted the number of BH-MS binaries originating in the latter site. In this paper, we estimate the number of BH-MS binaries formed in open clusters and detectable with {it Gaia} based on the results of {it N}-body simulations. By considering interstellar extinction in the Milky Way (MW) and observational constraints, we predict $sim 10$ BH-MS binaries are observable. We also find that chemical abundance patterns of companion MSs will help us to identify the origin of the binaries as star clusters. Such MSs are not polluted by outflows of the BH progenitors, such as stellar winds and supernova ejecta. Chemical anomalies might be a good test to confirm the origin of binaries with relatively less massive MSs ($lesssim 5M_{odot}$), orbital periods ($sim 1.5;$year) and higher eccentricities ($e gtrsim 0.1$).