Open Access
Issue
A&A
Volume 666, October 2022
Article Number A194
Number of page(s) 5
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/202244257
Published online 26 October 2022

© G.-P. Li 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1. Introduction

Since the first gravitational-wave (GW) event, GW150914 (Abbott et al. 2016), detected by Advanced LIGO (The LIGO Scientific Collaboration 2015) and Advanced Virgo (Acernese et al. 2015), more than 80 stellar-mass binary black hole (BBH) mergers have been reported by the LIGO-Virgo-KAGRA (LVK) Collaboration (Abbott et al. 2019, 2021; The LIGO Scientific Collaboration 2021a,b). The two main formation channels of BBH mergers have been extensively discussed, and include isolated binary evolution and dynamical formation (Mapelli 2021; Mandel & Farmer 2022). Dynamical formation channels occur only in dense environments such as nuclear star clusters (NSCs) or globular clusters (GCs), although a body of work has shown that dynamical formation channels produce mergers in the field of the host galaxies (i.e., Michaely & Perets 2019, 2020; Raveh et al. 2022; Michaely & Naoz 2022). In the above classical classification, the dynamical formation channel from wide systems in the field is ignored.

Some special GW events (e.g., Abbott et al. 2019, 2020a,b) with significant asymmetric masses, massive components, or high spins are promising candidates for hierarchical BBH mergers – with at least one of the black holes (BHs) being the remnant of a previous merger (Gerosa & Fishbach 2021). In particular, GW190521 (Abbott et al. 2020b) is regarded as a promising candidate to date (Anagnostou et al. 2020; Dall’Amico et al. 2021; Kimball et al. 2021; Baibhav et al. 2021; Fragione et al. 2022). Hierarchical mergers can occur via dynamical formation channels (Miller & Colbert 2004; Gerosa & Berti 2019; Yang et al. 2019; Dall’Amico 2021; Li 2022) or triple-body (and higher) configurations (Hamers & Safarzadeh 2020; Liu & Lai 2021). A primary condition for the occurrence of hierarchical mergers is that the remnants of first-generation BBH mergers are retained by the host star cluster (the binary also needs to survive few-body interactions). Therefore, the escape speed of the host star cluster needs to be larger than the typical recoil velocity (or “kick”) of remnant BHs.

Energy, angular momentum, and linear momentum are carried by GW radiation, in which the loss of linear momentum causes a kick imparted to a remnant BH. Measuring the kick velocity is significantly important in the inference of the hierarchical merger rate. The kick can be extracted from GW signals but with great uncertainty (Calderón Bustillo et al. 2018; Varma et al. 2020). The kick velocities inferred from the Gravitational Wave Transient Catalog (GWTC) events lie in a wide range: ∼50−2000 km s−1 (Mahapatra et al. 2021; Varma et al. 2022). In comparison, the escape speed of a typical GC is ∼2−100 km s−1 and that of an NSC is ∼10−600 km s−1 (Antonini & Rasio 2016).

Recently, Fragione & Loeb (2021) computed the distribution of kick velocities of the GWTC-2 events (Abbott et al. 2021) as a function of the spins of two BHs in a binary and found that remnant BHs with dimensionless spins as small as 0.1 would be retained by NSCs with escape speeds of higher than 100 km s−1. Doctor et al. (2021) obtained the kick velocity distributions of the GWTC-2 events using a numerical relativity (NR) surrogate waveform model (Varma et al. 2019a,b) and found that GCs and NSCs may retain up to ∼3% and ∼46% of their remnant BHs, respectively. Mahapatra et al. (2021) estimated the kick velocities of the GWTC-2 events using NR fitting formulas (Campanelli et al. 2007) and found that ∼50% of remnant BHs can be retained by NSCs with escape speeds of 200 km s−1. These results will independently constrain the rates of hierarchical mergers.

In this paper, we present a method to estimate the rates of hierarchical mergers in star clusters and give a simple formula in Sect. 2. We demonstrate that the rate of hierarchical merger is approximately proportional to the retention probability. In Sect. 3 we apply the formula to estimate hierarchical merger rates in GCs and NSCs. Finally, in Sect. 4 we discuss the comparison of our results with the simulations, and we summarise our main results in Sect. 5.

2. Method

2.1. Formula derivation

We assume that there are no hierarchical mergers currently taking place, that the BBH merger rate is average within a Hubble time, and that single stars burst at t = 0. We can approximate the merger rate in generic star clusters to

R cl = N cl N 1 g f BH 2 τ Hub , $$ \begin{aligned} R_{\rm cl} = \frac{N_{\rm cl}N_{\rm 1g}f_{\rm BH}}{2\tau _{\rm Hub}}, \end{aligned} $$(1)

where Ncl is the number density of star clusters, N1g (i.e., NBH) is the total number of first-generation (1g) BHs formed from the collapse of stars in the cluster, and fBH represents the fraction of BHs that form BH–BH binaries, harden, and merge within a Hubble time τHub. For example, if the cluster has N1g = 100 and fBH = 1, there are 50 mergers. We note that fBH in practice is not constant, which depends on mass ratio, total mass, and spins of BBHs. We then relax the above assumption and consider hierarchical mergers, rewriting Eq. (1) as

R cl N cl τ Hub ( 1 4 N 1 g f BH + 1 4 N 2 g f BH ) . $$ \begin{aligned} R_{\rm cl} \approx \frac{N_{\rm cl}}{\tau _{\rm Hub}} (\frac{1}{4}N_{\rm 1g}f_{\rm BH}+\frac{1}{4}N_{\rm 2g} f_{\rm BH}). \end{aligned} $$(2)

Hierarchical mergers will occur in the 2g population with N2g = N1g − N1gfBH + NremPret, with Nrem = RclτHub/Ncl being the total BBH merger number and Pret the retention probability of the remnant BHs in the cluster. We expect that the relationship between the 1g and 2g populations is that N2g ∼ N1g, because the number of BHs in the cluster that do not participate in mergers is much greater than the number that do, and some of the remnant BHs can be retained by the host cluster if the escape speed of the host is larger than their kick velocities. This illustrates that Eq. (2) can be derived naturally and correctly. For N2g (below), we use the exact value with N1g − N1gfBH + NremPret, and not the approximate value (∼N1g). Moreover, the mergers of a higher generation than 2g + 2g are extremely rare, and we therefore neglect them (Gerosa & Berti 2019; Rodriguez et al. 2019).

The 2g population (N2g) includes two parts: one is an initial BH population (N1g(1 − fBH)) and the other is a remnant BH population (NremPret); here, we limit our discussion to the 2g population without reference to the 1g population. Therefore, the merger from the 2g population could be a first-generation (1g + 1g) merger or a hierarchy (1g + 2g or 2g + 2g merger; we note that the merger from the 1g population must be a 1g + 1g merger). A key piece of information is whether any of the BHs come from the remnant BH population, because the condition for a hierarchical merger to occur is that at least one of the BHs comes from the remnant of a previous merger (i.e., the BH in the remnant BH population). Here, we introduce a hierarchical probability factor ϵ to represent the probability that a merger produced from the 2g population is a hierarchy, and correspondingly 1 − ϵ stands for the probability that this merger is a 1g + 1g merger. For ϵ, it can also be understood as the ratio of the number of 1g + 1g mergers from the 2g population to the total number of mergers, including 1g + 1g, 1g + 2g, and 2g + 2g mergers in the 2g population. If all the BHs come from the initial BH population, then the number of all combinations of 1g + 1g mergers from random pairings is C N 1 g ( 1 f B H ) 2 $ {C}^2_{N_{\rm 1g}\left(1-f_{\rm BH}\right)} $, where C n m = n ( n 1 ) ( n 2 ) ( n m + 1 ) m ! $ {C}^m_n = \frac{n(n-1)(n-2) \cdots (n-m+1)}{m!} $ is the number of all such combinations for taking m different elements from n (n ≥ m) different elements at a time regardless of their order. The number of all combinations of 1g + 1g, 1g + 2g, and 2g + 2g mergers in the 2g population is C N 2 g 2 $ {C}^2_{N_{\rm 2g}} $. Hence, ϵ has the following form:

ϵ = 1 C N 1 g ( 1 f BH ) 2 C N 2 g 2 = 1 N 1 g ( 1 f BH ) [ N 1 g ( 1 f BH ) 1 ] N 2 g ( N 2 g 1 ) 1 [ N 1 g ( 1 f BH ) N 2 g ] 2 , $$ \begin{aligned} \epsilon&= 1-\frac{{C}^2_{N_{\rm 1g}\left(1-f_{\rm BH}\right)}}{{C}^2_{N_{\rm 2g}}}\nonumber \\&= 1-\frac{N_{\rm 1g}\left(1-f_{\rm BH}\right) \left[N_{\rm 1g}\left(1-f_{\rm BH}\right)-1\right]}{N_{\rm 2g}(N_{\rm 2g}-1)} \approx 1-\left[\frac{N_{\rm 1g}\left(1-f_{\rm BH}\right)}{N_{\rm 2g}}\right]^2, \end{aligned} $$(3)

because N1g(1−fBH) ≫ 1.

The number of mergers from the 2g population (N2g) is 1 4 N 2 g f BH $ \frac{1}{4} N_{\mathrm{2g}}f_{\mathrm{BH}} $. Thus, the number of hierarchical mergers should be ϵ multiplied by 1 4 N 2 g f BH $ \frac{1}{4} N_{\mathrm{2g}}f_{\mathrm{BH}} $ in the cluster. In addition, the total number of BBH mergers that come from both the 1g population and the 2g population is Nrem. Consequently, the fraction of hierarchical mergers is obtained by dividing these two items:

f hier = ϵ N 2 g f BH 4 N rem . $$ \begin{aligned} f_{\rm hier} = \frac{\epsilon N_{\rm 2g}f_{\rm BH}}{4N_{\rm rem}}. \end{aligned} $$(4)

Therefore, the rate of hierarchical merger should be

R hier = R cl f hier = 1 2 R cl f BH P ret [ 1 f BH P ret 4 f BH ( 4 2 P ret ) ] . $$ \begin{aligned} R_{\rm hier}&= R_{\rm cl}f_{\rm hier}\nonumber \\&= \frac{1}{2} R_{\rm cl}f_{\rm BH}P_{\rm ret} \left[1-\frac{f_{\rm BH}P_{\rm ret}}{4-f_{\rm BH}\left(4-2P_{\rm ret}\right)}\right]. \end{aligned} $$(5)

Equation (5) can be approximated to (this is verified in Fig. 1)

R hier 1 2 R cl f BH P ret . $$ \begin{aligned} R_{\rm hier} \sim \frac{1}{2} R_{\rm cl}f_{\rm BH}P_{\rm ret}. \end{aligned} $$(6)

thumbnail Fig. 1.

Fractions (left) and rates (right) of hierarchical mergers as a function of the BBH merger rates (top) and retention probabilities (bottom) in GCs (red) and NSCs (black). The derived values of the fractions and rates are obtained according to Eqs. (4) and (5) (star symbols), respectively. The approximate values were calculated using Eqs. (7) and (6) (dashed line). The fractions and rates are plotted as a function of BBH merger rates (retention probabilities) with fixed retention probabilities of 4% and 40% (BBH merger rates fixed to 10 Gpc−3 yr−1 and 1 Gpc−3 yr−1) for GCs and NSCs, respectively. Ncl ∼ 1 Mpc−3 and ∼0.01 Mpc−3 and NBH ∼ 2000 and ∼2 × 104 are adopted for GCs and NSCs, respectively.

Equation (6) suggests that the rate of hierarchical merger is approximately proportional to the retention probability. Interestingly, the form of Eq. (6) is essentially the same as that of Eq. (17) of Kimball et al. (2020) who based their calculations on a comparison with simulations (Rodriguez et al. 2019). Then, combining Eqs. (5) and (6), we can rewrite Eq. (4) into an approximate but simple form as

f hier 1 2 f BH P ret . $$ \begin{aligned} f_{\rm hier} \sim \frac{1}{2}f_{\rm BH}P_{\rm ret}. \end{aligned} $$(7)

We further consider the branching ratio ( Γ 1 g + 2 g 2 g + 2 g $ \Gamma_\mathrm{\frac{\mathrm{1g + 2g}}{\mathrm{2g + 2g}}} $) of 1g + 2g to 2g + 2g merger rates. In the same way as Eq. (3), we can obtain hierarchical probability factors ϵ 1 g + 2 g = 1 C N 1 g ( 1 f BH ) 2 C N 2 g 2 C N rem P ret 2 C N 2 g 2 $ \epsilon_{\mathrm{1g+2g}} = 1-\frac{{C}^2_{N_{\mathrm{1g}}\left(1-f_{\mathrm{BH}}\right)}}{{C}^2_{N_{\mathrm{2g}}}}- \frac{{C}^2_{N_{\mathrm{rem}}P_{\mathrm{ret}}}}{{C}^2_{N_{\mathrm{2g}}}} $ for 1g + 2g mergers and ϵ 2 g + 2 g = C N rem P ret 2 C N 2 g 2 $ \epsilon_{\mathrm{2g+2g}} =\frac{{C}^2_{N_{\mathrm{rem}}P_{\mathrm{ret}}}}{{C}^2_{N_{\mathrm{2g}}}} $ for 2g + 2g mergers. Therefore, we expect

Γ 1 g + 2 g 2 g + 2 g = 1 C N 1 g ( 1 f BH ) 2 C N 2 g 2 C N rem P ret 2 C N 2 g 2 C N rem P ret 2 C N 2 g 2 4 ( 1 / f BH 1 ) P ret . $$ \begin{aligned} \Gamma _{\rm \frac{1g + 2g}{2g + 2g}}&= \frac{1- \frac{{C}^2_{N_{\rm 1g}\left(1-f_{\rm BH}\right)}}{{C}^2_{N_{\rm 2g}}}- \frac{\mathrm{C}^2_{N_{\rm rem}P_{\rm ret}}}{{C}^2_{N_{\rm 2g}}}}{\frac{{C}^2_{N_{\rm rem}P_{\rm ret}}}{{C}^2_{N_{\rm 2g}}}}\nonumber \\&\sim \frac{4(1/f_{\rm BH}-1)}{P_{\rm ret}}. \end{aligned} $$(8)

If we assume that the fraction of BHs fBH ≲ 20% (Antonini 2014; McKernan et al. 2018) and the retention probability Pret ≲ 80% (Mahapatra et al. 2021, also see Sect. 2.3), the branching ratio of hierarchical mergers is ≳20. This indicates that more than ∼95% of hierarchical mergers in clusters are 1g + 1g mergers rather than 2g + 2g mergers.

It is noteworthy that (1) some of the 2g binaries we calculate as successful mergers may not merge within a Hubble time because they must go through 1g mergers first; (2) we add remnant BHs to the 2g population at once, but the process should take place during the whole Hubble time; and (3) the densities of star clusters and the numbers of BHs in them are varying over cosmic time, and it is therefore unrealistic to expect them to be well-approximated as constants over a Hubble time, although they may be constant over shorter evolutionary windows. For these reasons, the simplicity of our assumptions should be taken into account when interpreting our results.

2.2. Density of star clusters and BH numbers

In general, roughly one star per 1000 M will be massive enough to form a BH, which is simply a result of the initial mass function of stars and the lowest birth mass needed to form a BH (∼20 M; O’Leary et al. 2009). Therefore, the number of BHs – which is going to be a strong function of radius in star clusters (especially in NSCs) – will scale with the mass of the cluster, and there is no unique answer because star clusters have different masses (O’Leary et al. 2009; Tsang 2013). We simply assume that there are ∼2000 (∼2 × 104) BHs in a GC (NSC) with a space density of around 1 (0.01) per Mpc3 (Baldry et al. 2012; Antonini 2014; Antonini & Rasio 2016). In comparison, there are ∼1 − 4 × 104 stellar mass BHs in the central parsec of the Milky Way (Miralda-Escudé & Gould 2000; Antonini 2014; Generozov et al. 2018).

2.3. Retention probability of remnant BHs

The retention probabilities of remnant BHs of ∼3% and 46% were obtained in GCs and NSCs respectively, from Doctor et al. (2021) who used a NR surrogate waveform model (Varma et al. 2019a,b). Following the calculations of Mahapatra et al. (2021), we find that clusters with escape speeds of 50 km s−1 and 250 km s−1 can retain ∼1%−12% and ∼14%−70% of the GWTC-2 events, respectively. Here, we adopt the retention probabilities of 4% and 40% in GCs and NSCs, respectively, as fiducial values and consider a wide range for comparison below.

3. Rate of hierarchical mergers

We assume that the BBH merger rate is uniform over the ∼10 Gyr history of the Universe, allowing us to obtain the fractions and rates of hierarchical mergers according to Eqs. (4) and (5) (star symbols) respectively, which we show in Fig. 1. We also use simplified Eqs. (6) and (7) (dashed line) to verify that they are correct.

From Fig. 1 (a1), we see that the fraction of hierarchical mergers in GCs is two orders of magnitude smaller than that in NSCs if they have the same BBH merger rates. If we consider that the typical BBH merger rates are ∼10 Gpc−3 yr−1 and ∼1 Gpc−3 yr−1 in GCs and NSCs (Mandel & Broekgaarden 2022) respectively, then the fraction is ∼2% in NSCs, while this value is ∼0.2% in GCs. We find that GCs and NSCs have hierarchical merger rates of ∼0.02 Gpc−3 yr−1 (see Fig. 1 (a2)).

Generally, the escape speeds in NSCs are larger than those in GCs, which causes the retention probability of remnant BHs in NSCs to be about one order of magnitude larger than that in GCs. This therefore suggests that hierarchical mergers are rarer in GCs. However, from Fig. 1 (a2) we find that the rates of hierarchical mergers in GCs and NSCs are the same, considering that the BBH merger rate in GCs is ten times larger than that in NSCs (Mandel & Broekgaarden 2022). We note that the similar rates of hierarchical mergers of GCs and NSCs are due to the fact that the retention fraction is taken here to be 4% and 40% (see Sect. 2.3), respectively, and under the assumption that mergers have no dependency on physical parameters such as masses and mass ratio (see Sect. 2). Further, we calculate the ratio of hierarchical mergers in NSCs to those in GCs using Eq. (6):

Γ NSC GC ( hier ) = f BH , NSC f BH , GC R NSC R GC P ret , NSC P ret , GC , $$ \begin{aligned} \Gamma _{\rm \frac{NSC}{GC}}^\mathrm{(hier)}&= \frac{f_{\rm BH,NSC}}{f_{\rm BH,GC}}\frac{R_{\rm NSC}}{R_{\rm GC}}\frac{P_{\rm ret,NSC}}{P_{\rm ret,GC}}, \end{aligned} $$(9)

where fBH, NSC (fBH, GC) is the fraction of BHs that form BH–BH binaries, harden, and merge in NSCs (GCs), and Pret, NSC and Pret, GC are the retention probabilities for NSCs and GCs, respectively. We set that RNSC/RGC ∼ 0.1 (Mandel & Broekgaarden 2022) and Pret, NSC/Pret, GC ∼ 10 (Doctor et al. 2021; Mahapatra et al. 2021). Therefore, we obtain a ratio of hierarchical merger rate in NSCs to that in GCs of ∼1 if fBH, NSC ≈ fBH, GC. Even if the ratio (fBH, NSC/fBH, GC) is not equal to 1, its fluctuation would not be significantly large. In consequence, the hierarchical merger rate in GCs is comparable to that in NSCs.

We relax the retention probability for GCs up to ∼15% in Fig. 1 (a3) and see that the fraction of hierarchical mergers in GCs can reach as high as that in NSCs. The corresponding hierarchical merger rate in GCs can even exceed that in NSCs (see Fig. 1 (a4)). We also adopt wide ranges in BH numbers and consider different BBH merger rates and retention probabilities to estimate fractions, rates, and branching ratios. The obtained results are listed in Table 1. We see that the hierarchical merger rate is affected by significant uncertainty, with a wide range from 𝒪(10−5) to 𝒪(1) Gpc−3 yr−1. Overall, the rate of hierarchical mergers is very small. If we assume that all the BBH mergers detected by LIGO-Virgo came from GCs or NSCs, then at most one event might instead be a hierarchical merger. We note that a wide range of BBH merger rates is adopted in Table 1. Therefore, it is normal that merger rates may be overestimated or underestimated compared with the results from the LVK Collaboration, who found the BBH merger rate to be between 17.9 Gpc−3 yr−1 and 44 Gpc−3 yr−1 at a fiducial redshift (z = 0.2; The LIGO Scientific Collaboration 2021c).

Table 1.

Fractions, rates, and branching ratios of hierarchical mergers for different cases.

4. Discussion

Rodriguez et al. (2019) used a grid of 96 dynamical models of dense star clusters and a cosmological model of cluster formation to explore the production of hierarchical BBH mergers. These latter authors found that if all initial BHs that formed from the collapse of stars are nonspinning (i.e., BH dimensionless spin parameter χ = 0), then more than 10% of mergers are hierarchical mergers. Following the method of Baibhav et al. (2021) to pair binaries, we estimate the retention probability of 1g + 1g mergers by calculating kick velocities of merging binaries using NR fitting formulas (Campanelli et al. 2007). We find that more than 80% of remnant BHs would be retained by clusters with escape speeds of ∼200 km s−1 if all 1g BHs are nonspinning. This is relatively high compared to existing observations and calculations (Fragione & Loeb 2021; Mahapatra et al. 2021; Kimball et al. 2021), therefore the contribution of hierarchical mergers from GCs to the BBH merger rate might drop from ∼10% to ≲1% for moderate χ ∼ 0.5 (Rodriguez et al. 2019). This is consistent with our result that the fractions of hierarchical mergers fhier are all less than 10%; see Table 1 (also see Fig. 1 (a1) and (a3)). The simplified model of Gerosa & Berti (2019) also shows a relatively small fraction of hierarchical mergers in clusters if the BH spin is as much as 0.2. Moreover, Mapelli et al. (2021b) found that hierarchical BBHs in NSCs account for ∼10−2 − 0.2 Gpc−3 yr−1, which broadly agrees with our estimates (see Fig. 1 (a4)).

The 2g BHs may preferentially form binaries in dynamical interactions due to more massive components. As a result, the fraction fBH is different in the 1g population and the 2g population in Eq. (2). Furthermore, there are some initial binaries in the 1g population, which would make fBH smaller in 1g + 1g mergers and larger in higher generation mergers. To justify the use of the same fBH for the 1g and 2g populations, we set the fraction in the 2g population to n multiplied by fBH in the 1g population. The factor n is larger than 1 but limited to 1/fBH. Then, Eq. (2) is

R cl N cl τ Hub ( 1 4 N 1 g f BH + 1 4 N 2 g n f BH ) . $$ \begin{aligned} R_{\rm cl} \approx \frac{N_{\rm cl}}{\tau _{\rm Hub}} (\frac{1}{4}N_{\rm 1g}f_{\rm BH}+\frac{1}{4}N_{\rm 2g} nf_{\rm BH}). \end{aligned} $$(10)

We obtain the rate of hierarchical mergers using Eqs. (3) and (4):

R hier ( n ) = 1 2 n R cl f BH P ret [ 1 f BH P ret 4 f BH ( 4 2 P ret ) ] . $$ \begin{aligned} R^{(n)}_{\rm hier} = \frac{1}{2} nR_{\rm cl}f_{\rm BH}P_{\rm ret} \left[1-\frac{f_{\rm BH}P_{\rm ret}}{4-f_{\rm BH}\left(4-2P_{\rm ret}\right)}\right]. \end{aligned} $$(11)

We can define a ratio Γ n 1 $ \Gamma_{\frac{n}{1}} $, which is equal to the rate from Eq. (11) divided by that from Eq. (5):

Γ n 1 = R hier ( n ) R hier = n . $$ \begin{aligned} \Gamma _{\frac{n}{1}} = \frac{R^{(n)}_{\rm hier}}{R_{\rm hier}}=n. \end{aligned} $$(12)

We expect n to be between 1 and 2 in general due to the fact that the binary needs to survive few-body interactions, and therefore the difference between considering the change of the fraction and not considering the change of the fraction results in a change of less than a factor two for the rates of hierarchical mergers. This effect is relatively small and the results from Eq. (5) are reliable.

5. Conclusion

The kick can be extracted from the GW signals. This information can then be used to obtain the retention probability of remnant BHs in different types of star clusters assuming all GWTC events occurred in clusters.

In this study, we attempted to estimate the rates of hierarchical merger from the retention probability of remnant BHs. A simple formula is given by Eq. (6) that suggests the rate of hierarchical merger is proportional to the retention probability. We find that ∼2% of BBH mergers in NSCs may be hierarchical mergers, while this percentage in GCs is a few tenths of a percent. However, this does not mean the rate of hierarchical merger in GCs is smaller than that in NSCs. We find that hierarchical merger rates in GCs and NSCs are about the same with ∼𝒪(10−2) Gpc−3 yr−1 because the BBH merger rate in GCs is actually higher than that in NSCs.

Similar conclusions on hierarchical mergers (and rates) have already been obtained by various authors (Rodriguez et al. 2019; Kimball et al. 2020; Arca Sedda et al. 2020; Fragione & Silk 2020; Mapelli et al. 2021a,b) using even more sophisticated techniques (i.e., less approximate). Nevertheless, simple analytic arguments can provide useful insights because our results are not based on simulations. Moreover, we also point out that although the efficiency of hierarchical mergers in GCs is lower than that in NSCs, this does not mean that the hierarchical merger rate in GCs is lower than that in NSCs (see Sect. 3). This provides the warning that if an event detected by LIGO-Virgo is determined as a hierarchical merger, then the probability that it originates from a GC is approximately equal to the probability that it originates from an NSC.

Understanding the distribution of BHs at the centers of galaxies is crucial for making predictions about the expected merger rates. We expect that with third-generation GW detectors in operation (Punturo et al. 2010a,b; Abbott et al. 2017), the increasing data on GW events will help us to constrain the retention probability and merger rate more precisely.

Acknowledgments

G.-P. Li is grateful to Da-Bin Lin, Ilya Mandel, Xiao-Yan Li, Jia-Xin Cao, Cheng Gao, and Hao-Qiang Zhang for useful discussions. This work is supported by the National Natural Science Foundation of China (Grant No. 11773007) and the Guangxi Science Foundation (Grant No. 2018GXNSFFA281010).

References

  1. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, Phys. Rev. Lett., 116, 061102 [Google Scholar]
  2. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, Class. Quant. Grav., 34, 044001 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2019, Phys. Rev. X, 9, 031040 [Google Scholar]
  4. Abbott, R., Abbott, T. D., Abraham, S., et al. 2020a, Phys. Rev. D, 102, 043015 [NASA ADS] [CrossRef] [Google Scholar]
  5. Abbott, R., Abbott, T. D., Abraham, S., et al. 2020b, Phys. Rev. Lett., 125, 101102 [Google Scholar]
  6. Abbott, R., Abbott, T. D., Abraham, S., et al. 2021, Phys. Rev. X, 11, 021053 [Google Scholar]
  7. Acernese, F., Agathos, M., Agatsuma, K., et al. 2015, Class. Quant. Grav., 32, 024001 [Google Scholar]
  8. Anagnostou, O., Trenti, M., & Melatos, A. 2020, ArXiv e-prints [arXiv:2010.06161] [Google Scholar]
  9. Antonini, F. 2014, ApJ, 794, 106 [NASA ADS] [CrossRef] [Google Scholar]
  10. Antonini, F., & Rasio, F. A. 2016, ApJ, 831, 187 [Google Scholar]
  11. Arca Sedda, M., Mapelli, M., Spera, M., Benacquista, M., & Giacobbo, N. 2020, ApJ, 894, 133 [NASA ADS] [CrossRef] [Google Scholar]
  12. Baibhav, V., Berti, E., Gerosa, D., Mould, M., & Wong, K. W. K. 2021, Phys. Rev. D, 104, 084002 [NASA ADS] [CrossRef] [Google Scholar]
  13. Baldry, I. K., Driver, S. P., Loveday, J., et al. 2012, MNRAS, 421, 621 [NASA ADS] [Google Scholar]
  14. Calderón Bustillo, J., Clark, J. A., Laguna, P., & Shoemaker, D. 2018, Phys. Rev. Lett., 121, 191102 [CrossRef] [Google Scholar]
  15. Campanelli, M., Lousto, C., Zlochower, Y., & Merritt, D. 2007, ApJ, 659, L5 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dall’Amico, M. 2021, ArXiv e-prints [arXiv:2112.02020] [Google Scholar]
  17. Dall’Amico, M., Mapelli, M., Di Carlo, U. N., et al. 2021, MNRAS, 508, 3045 [CrossRef] [Google Scholar]
  18. Doctor, Z., Farr, B., & Holz, D. E. 2021, ApJ, 914, L18 [NASA ADS] [CrossRef] [Google Scholar]
  19. Fragione, G., & Loeb, A. 2021, MNRAS, 502, 3879 [NASA ADS] [CrossRef] [Google Scholar]
  20. Fragione, G., & Silk, J. 2020, MNRAS, 498, 4591 [NASA ADS] [CrossRef] [Google Scholar]
  21. Fragione, G., Kocsis, B., Rasio, F. A., & Silk, J. 2022, ApJ, 927, 231 [NASA ADS] [CrossRef] [Google Scholar]
  22. Generozov, A., Stone, N. C., Metzger, B. D., & Ostriker, J. P. 2018, MNRAS, 478, 4030 [Google Scholar]
  23. Gerosa, D., & Berti, E. 2019, Phys. Rev. D, 100, 041301 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gerosa, D., & Fishbach, M. 2021, Nat. Astron., 5, 749 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hamers, A. S., & Safarzadeh, M. 2020, ApJ, 898, 99 [NASA ADS] [CrossRef] [Google Scholar]
  26. Kimball, C., Talbot, C., Berry, C. P. L., et al. 2020, ApJ, 900, 177 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kimball, C., Talbot, C., Berry, C. P. L., et al. 2021, ApJ, 915, L35 [CrossRef] [Google Scholar]
  28. Li, G.-P. 2022, Phys. Rev. D, 105, 063006 [NASA ADS] [CrossRef] [Google Scholar]
  29. Liu, B., & Lai, D. 2021, MNRAS, 502, 2049 [NASA ADS] [CrossRef] [Google Scholar]
  30. Mahapatra, P., Gupta, A., Favata, M., Arun, K. G., & Sathyaprakash, B. S. 2021, ApJ, 918, L31 [NASA ADS] [CrossRef] [Google Scholar]
  31. Mandel, I., & Broekgaarden, F. S. 2022, Liv. Rev. Relativ., 25, 1 [NASA ADS] [CrossRef] [Google Scholar]
  32. Mandel, I., & Farmer, A. 2022, Phys. Rep., 955, 1 [NASA ADS] [CrossRef] [Google Scholar]
  33. Mapelli, M. 2021, Handbook of Gravitational Wave Astronomy, 4 [Google Scholar]
  34. Mapelli, M., Dall’Amico, M., Bouffanais, Y., et al. 2021a, MNRAS, 505, 339 [CrossRef] [Google Scholar]
  35. Mapelli, M., Santoliquido, F., Bouffanais, Y., et al. 2021b, Symmetry, 13, 1678 [NASA ADS] [CrossRef] [Google Scholar]
  36. McKernan, B., Ford, K. E. S., Bellovary, J., et al. 2018, ApJ, 866, 66 [Google Scholar]
  37. Michaely, E., & Naoz, S. 2022, ApJ, 936, 2 [NASA ADS] [CrossRef] [Google Scholar]
  38. Michaely, E., & Perets, H. B. 2019, ApJ, 887, L36 [NASA ADS] [CrossRef] [Google Scholar]
  39. Michaely, E., & Perets, H. B. 2020, MNRAS, 498, 4924 [NASA ADS] [CrossRef] [Google Scholar]
  40. Miller, M. C., & Colbert, E. J. M. 2004, Int. J. Mod. Phys. D, 13, 1 [NASA ADS] [CrossRef] [Google Scholar]
  41. Miralda-Escudé, J., & Gould, A. 2000, ApJ, 545, 847 [CrossRef] [Google Scholar]
  42. O’Leary, R. M., Kocsis, B., & Loeb, A. 2009, MNRAS, 395, 2127 [CrossRef] [Google Scholar]
  43. Punturo, M., Abernathy, M., Acernese, F., et al. 2010a, Class. Quant. Grav., 27, 194002 [Google Scholar]
  44. Punturo, M., Abernathy, M., Acernese, F., et al. 2010b, Class. Quant. Grav., 27, 084007 [NASA ADS] [CrossRef] [Google Scholar]
  45. Raveh, Y., Michaely, E., & Perets, H. B. 2022, MNRAS, 514, 4246 [NASA ADS] [CrossRef] [Google Scholar]
  46. Rodriguez, C. L., Zevin, M., Amaro-Seoane, P., et al. 2019, Phys. Rev. D, 100, 043027 [Google Scholar]
  47. The LIGO Scientific Collaboration (Aasi, J., et al.) 2015, Class. Quant. Grav., 32, 074001 [Google Scholar]
  48. The LIGO Scientific Collaboration, the Virgo Collaboration (Abbott, R., et al.) 2021a, ArXiv e-prints [arXiv:2108.01045] [Google Scholar]
  49. The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, et al. 2021b, ArXiv e-prints [arXiv:2111.03606] [Google Scholar]
  50. The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, et al. 2021c, ArXiv e-prints [arXiv:2111.03634] [Google Scholar]
  51. Tsang, D. 2013, ApJ, 777, 103 [NASA ADS] [CrossRef] [Google Scholar]
  52. Varma, V., Field, S. E., Scheel, M. A., et al. 2019a, Phys. Rev. Res., 1, 033015 [NASA ADS] [CrossRef] [Google Scholar]
  53. Varma, V., Gerosa, D., Stein, L. C., Hébert, F. M. C., & Zhang, H. 2019b, Phys. Rev. Lett., 122, 011101 [NASA ADS] [CrossRef] [Google Scholar]
  54. Varma, V., Isi, M., & Biscoveanu, S. 2020, Phys. Rev. Lett., 124, 101104 [NASA ADS] [CrossRef] [Google Scholar]
  55. Varma, V., Biscoveanu, S., Islam, T., et al. 2022, Phys. Rev. Lett., 128, 191102 [NASA ADS] [CrossRef] [Google Scholar]
  56. Yang, Y., Bartos, I., Gayathri, V., et al. 2019, Phys. Rev. Lett., 123, 181101 [CrossRef] [Google Scholar]

All Tables

Table 1.

Fractions, rates, and branching ratios of hierarchical mergers for different cases.

All Figures

thumbnail Fig. 1.

Fractions (left) and rates (right) of hierarchical mergers as a function of the BBH merger rates (top) and retention probabilities (bottom) in GCs (red) and NSCs (black). The derived values of the fractions and rates are obtained according to Eqs. (4) and (5) (star symbols), respectively. The approximate values were calculated using Eqs. (7) and (6) (dashed line). The fractions and rates are plotted as a function of BBH merger rates (retention probabilities) with fixed retention probabilities of 4% and 40% (BBH merger rates fixed to 10 Gpc−3 yr−1 and 1 Gpc−3 yr−1) for GCs and NSCs, respectively. Ncl ∼ 1 Mpc−3 and ∼0.01 Mpc−3 and NBH ∼ 2000 and ∼2 × 104 are adopted for GCs and NSCs, respectively.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.