GW200105 and GW200115 are compatible with a scenario of primordial black hole binary coalescences

Two gravitational wave events, i.e. GW200105 and GW200115, were observed by the Advanced LIGO and Virgo detectors recently. In this work, we show that they can be explained by a scenario of primordial black hole binaries that are formed in the early Universe. The merger rate predicted by such a scenario could be consistent with the one estimated from LIGO and Virgo, even if primordial black holes constitute a fraction of cold dark matter. The required abundance of primordial black holes is compatible with the existing upper limits from microlensing, caustic crossing and cosmic microwave background observations.


Introduction
Based on a second half of the third observing run, the LIGO Scientific and Virgo Collaborations (LVC) [1] reported two gravitational wave events, namely GW200105 and GW200115, which are compatible with neutron star-black hole (NSBH) binaries. At 90% confidence level, the primary components were found to be black holes with masses of 8.9 +1.2 −1.5 M and 5.7 +1.8 −2.1 M , respectively, while the secondary ones 1.9 +0. 3 −0.2 M and 1.5 +0.7 −0.3 M , respectively. Compared with the maximal mass of neutron stars, both of the secondaries are compatible with neutron stars with a probability of ∼ 90%. When the events were assumed to be representative of the entire NSBH population, the merger rate densities of GW200105 and GW200115 were inferred to be 16 +38 −14 Gpc −3 year −1 and 36 +82 −30 Gpc −3 year −1 , respectively. When a broader distribution of component masses was assumed, they became 130 +112 −69 Gpc −3 year −1 . The above rates were found to be consistent with the model predictions of NSBH formation in isolated binaries or young star clusters. For either of the two events, however, LVC found no a e-mail: zhaozc@bnu.edu.cn (corresponding author) evidence of tides or tidal disruption meanwhile identified no electromagnetic counterparts. Therefore, there is no direct evidence that the secondaries are neutron stars.
In this work, we show that a scenario of primordial black hole (PBH) binary coalescences can explain the origin of GW200105 and GW200115. Gravitational collapse of enhanced overdensities could produce PBHs in the early Universe [2][3][4][5][6][7][8][9][10][11][12]. PBHs can form binaries through several formation channels [13,14]. Amongst these channels, the most efficient one was originally proposed in Refs. [14,15], and recently extensively studied in Refs. [16][17][18] by considering such as a generic mass distribution of PBHs and the tidal forces from all the neighboring PBHs and linear density perturbations. To determine if GW200105 and GW200115 are primordial, we have to check whether their merger rate densities can be accounted for with the corresponding PBH abundance allowed by the existing upper limits (see review in Ref. [19] and references therein). In the following, we will compute the merger rate density of PBHs in two specific models by following Ref. [20]. The first one assumes that the mass function is determined by GW200105 and GW200115. The second one assumes that all the black hole binary coalescences observed by LVC are primordial. The predicted rates will be compared with the estimated ones from LVC.
The rest of the paper is arranged as follows. In Sect. 2, we show the formula for the merger rate of PBH binaries. In Sect. 3, we briefly review the log-normal mass function of PBHs. Our results and conclusions are shown in Sects. 4 and 5, respectively.

Merger rate of PBH binaries
We consider the formation channel of PBH binaries in the early Universe [14]. It is known that this channel makes a dominant contribution to the PBH merger rate [17]. This means neglects of the binary formation mechanisms, such as dynamical captures and three-body interactions, in the late Universe [21]. The late-Universe formation mechanism was originally proposed in Ref. [14] and revisited recently in Ref. [15]. Two neighboring PBHs could form a binary due to torque from a third neighboring PBH, and merge with each other within the age of the Universe due to the energy loss via gravitational emissions. However, the mass distribution of PBHs was assumed to be monochromatic in the original literature. More recently, the merger rate of PBHs with an extended mass function was studied in Refs. [16][17][18]. To be specific, Ref. [16] took into account the tidal force from a PBH which is closest to the center of mass of PBH binary. Reference [17] assumed a flat mass function of PBHs within a narrow mass range. To generalize the above two works, Ref. [18] took into account the tidal forces from all neighboring PBHs and linear density perturbations, and meanwhile considered a generic mass function of PBHs. Therefore, we follow Ref. [18] to compute the merger rate of PBH binaries with a log-normal mass distribution.
For PBHs within mass intervals of (m i , m i + dm i ) and (m j , m j + dm j ), the merger rate per unit volume The comoving merger rate density is obtained as follows [18] where f = Ω pbh /Ω m is the total fraction of non-relativistic matter in PBHs, σ eq ≈ 0.005 is the variance of overdensities of the rest of dark matter on scales of order O(10 −2 − 10 5 )M at equality [17], min(x 1 , x 2 ) selects the minimal value between x 1 and x 2 , P(m) is the mass function of PBHs, t and τ denote the cosmic time and the present age of the universe, respectively. The abundance of PBHs in dark matter is Here, Ω pbh , Ω dm and Ω m denote the energy density fractions of PBHs, dark matter and non-relativistic matter, respectively, in the critical energy density of the Universe at present. In addition, we would not consider the redshift evolution of the merger rate density, since only the low-redshift events are detectable for LVC. Throughout this work, we fix all cosmological parameters to be the best-fit values from the Planck 2018 results [22].

Mass function of PBHs
The mass function of PBHs can be determined by the production mechanisms of PBHs in the early Universe. The most likely one is based on the gravitational collapse of overdensities in the radiation dominated epoch of the Universe [2][3][4][5][6][7][8]. Therefore, the mass function of PBHs depends on the properties of primordial curvature perturbations. It is usually parametrized to be a log-normal function as follows where m * and σ stand for the center mass and width, respectively [18,20,23]. In this work, we consider two different methods to determine values of m * and σ . First, we assume that all observed black hole binaries are primordial. The results have been shown to be m * = 19M and σ = 0.97 [20], which will be used in this work. There are also alternative choices of mass function [18], which would not significantly change our results of the abundance of PBHs. Second, we assume that only GW200105 and GW200115 are of primordial origin and they determine the mass function of PBHs.
Typically, we choose one half of the total mass to determine the value of m * , namely m * = 5.4M for GW200105 and m * = 3.6M for GW200115. Meanwhile, we set the width to be σ = m * /10M , implying σ = 0.54 for GW200105 and σ = 0.36 for GW200115. In such a case, slightly different choices would change our predictions on the merger rate density by a factor of O(1). Based on Eq. (1), they would not significantly alter our results of the abundance of PBHs [20].

Limits on PBH abundance
Our results are shown in Fig. 1. We depict the limits on f pbh estimated from the merger rates of GW200105 (red triangle) and GW200115 (green triangle), and the broad mass function (purple circle), respectively. The gray solid lines denote the error bars at 90% CL, due to uncertainties in the measurements of the merger rates [1]. The dotted curves stand for the existing upper limits at 90% CL from OGLE (blue) [24], EROS/MACHO (cyan) [25,26], ICARUS (purple) [27], and Planck (orange) [28], which are shown here for comparison. In the case of single PBH event, the mass function is determined by the considered event, i.e. GW200105 or GW200115. Given the observational uncertainty, we derive the allowed abundance of PBHs which matches the observed merger rate of the event [1]. We find that both of the events are well explained by such a scenario, when we take f pbh O(10 −2 ) that is compatible with the existing best constraint f pbh O(10 −1 ) from caustic crossing [27]. We also find that the inferred upper limits are close to the existing upper limit, implying that such a scenario may be further tested in the near future. Furthermore, it is interesting to find that the mass function with m * = 5.4M and σ = 0.54 can almost precisely explain GW200105 and GW200115 simultaneously. But this coincidence may be accidental.
In the case of the broad mass distribution, the mass function is determined by all the observed black hole binaries, since all of them are assumed to be primordial. This has been done by using maximum-likelihood analysis in an existing literature [20]. In our work, the merger rate is predicted via integration over m i ∈ [2.5, 40]M and m j ∈ [1, 3]M , following a method consistent with Ref. [1]. Given f pbh O(10 −2 ), we conclude that the merger rates predicted by this model are compatible with those of GW200105 and GW200115 reported by LVC [1], although the predicted constraint on f pbh is in tension with that inferred from the cosmic microwave background measured by Planck [28]. In addition, by combining the broad mass function from Ref. [20] with the two recently observed events, i.e. GW200105 and GW200115, we obtain a joint result of m * = 16.4M and σ = 1.0, which is slightly different with that of Ref. [20].

Conclusions
In this work, we have shown that the scenario of PBH binary coalescences can account for the observed merger rates of GW200105 and GW200115. We computed the merger rates in two different models. The first one assumed that only GW200105 and GW200115 have a primordial origin. The second one assumed that all the observed black hole binaries are primordial. Given f pbh O(10 −2 ), we found that both of the two models can explain the recently reported events. For the former, the required abundance of PBHs were found to be compatible with the existing upper limits from other observations. The inferred upper limits were shown to be close to the existing upper limits, implying that such scenarios can be further tested by observations of (SGWB) in the near future [29][30][31][32][33]. However, the latter was found to be in tension with the exiting observational limits. In addition, we found that the log-normal mass function of PBHs with m * = 5.4M and σ = 0.54 almost simultaneously explains the observed merger rates of GW200105 and GW200115, including the center values and error bars. In summary, depending on the abundance of PBHs, the scenario of PBH binaries formed through the early-Universe channel could well explain GW200105 and GW200115 events. of Sciences with Grant NO. XDPB15, and the science research grants from the China Manned Space Project with NO. CMS-CSST-2021-B01. Z.C.Z. is supported by the National Natural Science Foundation of China with Grant No. 12005016.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The data is deposited in the text.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .