Open Access
Issue
A&A
Volume 667, November 2022
Article Number A126
Number of page(s) 15
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202039418
Published online 16 November 2022

© P. Drozda et al. 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

Detailed analyses of X-ray binary observations have found a paucity of compact remnants, neutron stars (NSs) and black holes (BHs), in a mass range from ∼2 − 5 M (Bailyn et al. 1998; Özel et al. 2010). These observations contradicted the prevailing theory at the time that argued for a continuous distribution of compact remnants across this range (Fryer & Kalogera 2001). One explanation for this mass gap is that the current systems were the result of an observational bias (Farr et al. 2011; Kreidberg et al. 2012). Certainly, there has been a set of observed systems with measurements within the mass gap region, for example: the recent observation of an NS with a mass of 2 . 14 0.09 + 0.10 M $ 2.14^{+0.10}_{-0.09}\,{M}_{\odot} $ (Cromartie et al. 2020), the low-mass X-ray binary with a black hole mass of 3 . 3 0.7 + 2.8 M $ 3.3^{+2.8}_{-0.7}\,{M}_{\odot} $ (Thompson et al. 2019), and the ambiguous smaller component in the merger event GW190814 with mass m 2 = 2 . 59 0.09 + 0.08 M $ m_2=2.59^{+0.08}_{-0.09}\,{M}_{\odot} $ (Abbott et al. 2020). It was argued that the existence, or nonexistence, of compact objects in this gap region could be used to constrain the properties of the core-collapse supernova engine (Belczynski et al. 2012b). Alternatively, compact objects in the mass gap could be created by mergers of lighter compact objects, and thus probe globular cluster dynamics. This would allow such compact objects to acquire companions, which leads to second generation merges (Fragione & Loeb 2019; Fragione et al. 2020).

So far, no BH-NS binaries have been detected through electromagnetic (EM) observations (Liu et al. 2014; Bhattacharya et al. 2019). Although observations of X-ray binaries suggest that such systems should exist, the formation rate of merging BH-NS binaries is not well constrained through EM observations (Belczynski et al. 2011, 2012a, 2013; Beldycki & Belczynski 2016).

The LIGO/Virgo gravitational wave detectors are ideally suited to studying the properties of compact remnants and answering questions about the existence of a mass gap. The primary gravitational wave (GW) signals for this ground-based detector consortium arise from the merger of double compact objects: BH-BH, BH-NS, and NS-NS. The first two LIGO/Virgo science runs (O1 and O2) detected 11 merging compact binaries: 10 BH-BH and 1 NS-NS (Abbott et al. 2019a). Although these observations brought a wealth of information about BHs and NSs, no BH-NS mergers were observed and no compact object was reported to have mass in the mass gap region (Abbott et al. 2019b).

This picture is changing as data from the third (O3) LIGO/Virgo science run are published. Table 1 shows the most confident candidates associated with the mass gap from the O3 science run. Four events were classified as BH-NS mergers (Abbott et al. 2021; The LIGO Scientific Collaboration 2021), and one as either a BH-BH or a BH-NS merger because the lighter compact object lies in the first mass gap (FMG) and its nature is currently unknown (Abbott et al. 2020).

Table 1.

LIGO/Virgo BH-NS merger candidates.

We note that LIGO/Virgo’s prompt classification scheme adopts specific choices about how to characterize compact objects that differ from those typically assumed in the literature: NSs are objects with mass M < 3 M, FMG objects have mass 3 < M < 5 M, and BHs have mass M > 5 M. This classification means that any type of actual double compact object (NS-NS, BH-NS, or BH-BH) can have one or two FMG objects as long as the true mass limit discriminating between NSs and BHs is anywhere between ∼2 − 5 M.

GW (gravitational wave) ground based observations can identify the FMG objects, but they cannot easily distinguish between low-mass BHs or high-mass NSs (e.g., Hinderer et al. 2019). EM observations can possibly help to identify the nature of a given FMG object. For example, the expectation is that BH-BH mergers do not produce EM emission, and no EM counterpart has been found so far for such a merger (e.g., Greiner et al. 2016). For BH-NS and NS-NS mergers, various EM counterparts are expected across the EM spectrum. The most notable are short gamma-ray bursts (Paczyński 1991; Lee & Kluzniak 1995; Ruffert et al. 1997) and kilonovae (Li & Paczyński 1998). Although there are expected differences between the composition of the ejecta, depending on the size and nature (NS vs. BH) of the compact remnant, the differences are currently difficult to disentangle from the many model uncertainties (e.g., Metzger 2017; Miller et al. 2019a; Korobkin et al. 2021). However, by better understanding BH-NS properties, there is hope to better differentiate between BH-NS and NS-NS systems.

In this study we focus on BH-NS mergers that form in classical (common envelope (CE)) isolated binary evolution (e.g., Belczynski et al. 2016a). We estimate the typical physical properties of these mergers that LIGO/Virgo could detect. Using population synthesis models that include updated estimates of BH spins and masses (Belczynski et al. 2020), we estimate the fraction of BH-NS mergers that host FMG objects. We also estimate the amount of mass ejection during BH-NS mergers, estimating the nature and detectability of any associated kilonova emission.

2. Calculations

2.1. Evolutionary calculations

In this paper, we evolve populations of massive star binaries in order to estimate the population of BH-NS binaries, focusing on those systems that merge in the Hubble time. All models we present were obtained using the StarTrack population synthesis code (Belczynski et al. 2002, 2008a, 2020) under the isolated classical (CE) binary evolution scenario. We calculate a series of models with different prescriptions for stellar evolution and binary input physics. For each model, we use a range of metallicity Z from 0.0001 up to 0.03: 32 sub-models with the same input physics but with different metallicity, each calculated for 2 × 107 massive binaries. In this section, we review the assumed initial conditions and the basic prescriptions used for binary evolution.

We used initial binary parameters from Sana et al. (2012) as modified by de Mink & Belczynski (2015), and stellar winds were adopted from Vink et al. (2001) and Belczynski et al. (2010). We used a 50% binary fraction and we assumed a maximum NS mass of 2.5 M. We adopted a solar metallicity of Z = 0.014.

We adopt three component broken power-law initial mass functions (IMFs) from Kroupa (2001) for the primary (more massive) component in each binary: MZAMS, A ∝ M−1.3 for 0.08 ≤ M < 0.5 M, ∝M−2.2 for 0.5 ≤ M < 1.0 M, and ∝Mα3 for 1.0 ≤ M < 150 M, where we adopt α3 = 2.3. The initial mass of the secondary (MZAMS, B) binary component is taken from the flat mass ratio distribution in a range q0 = [0.08/MZAMS, A, 1], where q0 = MZAMS, B/MZAMS, A. The lower range limit is chosen in such a way in order to provide second component mass above the hydrogen burning limit (0.08 M).

In our CE calculations, we use an energy balance approach from Webbink (1984), with updates on binding envelope energy from Dominik et al. (2012). We do not take into account systems in which a CE with a Hertzsprung gap donor occurred, as we conservatively assume these merge during the CE phase. The core-envelope structure of these stars is not well known (Belczynski et al. 2007; Ivanova et al. 2013) and the survival of these systems during CE is highly uncertain. Additionally, such systems may evolve through thermal timescale Roche-lobe overflow (RLOF) rather than through CE (Pavlovskii et al. 2017). We adopted an accretion rate of 5% of the Bondi-Hoyle rate (MacLeod et al. 2017) onto the NS and BH during the CE phase.

During stable RLOF with an NS and BH accretor, we calculate binary evolution and component masses using prescriptions from Mondal et al. (2020). In all the other cases, we assume a nonconservative mass transfer, with 50% of the transferred mass being lost from a binary with high specific angular momentum jloss = 1 (Podsiadlowski et al. 1992). The remaining 50% of the transferred mass is attached to a companion star.

Following Fryer et al. (2012), we use a neutrino supported convective supernova engine to deduce the masses of NSs and BHs in our simulations. We allow for a different development time of such an engine. In the rapid supernova model, for intermediate-mass stars, the engine develops very quickly (∼100 ms) and is followed by a supernova explosion and the ejection of stellar outer layers, and light NSs form. For massive stars, the engine is not able to overcome the weight of in-falling outer stellar layers and stars collapse to form rather massive BHs. This naturally creates a mass gap between NSs and BHs (Belczynski et al. 2012b). In the delayed supernova model, the engine develops after a relatively long time (∼500 − 1000 ms) after core collapse. This allows for significant accretion (∼1 − 2 M) onto light proto-NSs, and the formation of heavy NSs and light BHs producing a continuous mass spectrum of compact objects.

The FMG is typically defined as the range: ∼2 − 5 M. To reflect the recently discovered (in radio) NSs with masses of about 2.1 M (see Cromartie et al. 2020; Zhang & Li 2019), we adopt a slightly narrower range as a definition of the FMG: 2.1 − 5 M in our calculations.

For massive stars, we take into account pair-instability pulsation supernova mass loss. We adopt the weak mass loss that allows for the formation of BHs up to ∼55 M (Belczynski et al. 2020). We also allow for the most massive stars to be totally disrupted (no BH remnant) by pair-instability supernovae. During NS and BH formation, we assume 10% and 1% neutrino mass loss, respectively.

Compact objects during formation may receive natal kicks (Hobbs et al. 2005). We use a one dimensional Maxwellian natal kick magnitude distribution with σ = 265 km s−1. This generates a three dimensional average speed of ∼420 km s−1. The orientation of the natal kick is random. For some models, we allow kicks to be decreased in magnitude through fallback of matter during compact object formation. The natal kick decreases with decreasing ejecta mass (increasing fallback mass):

V kick = ( 1 f fb ) V , $$ \begin{aligned} V_{\rm kick}= (1-f_{\rm fb}) V, \end{aligned} $$(1)

where V is the kick magnitude drawn from the Maxwellian distribution with σ = 265 km s−1, and ffb is a fraction of matter that falls back onto a compact object after the supernova explosion. The fallback fraction fb is obtained from formulae by Fryer et al. (2012). It should be noted that this fallback of matter during core collapse is always used in the calculation of compact object mass, but is used only for some models in the calculation of natal kicks. In other words, when fallback is applied, we assume asymmetric mass ejection natal kicks, and when fallback is not applied, we assume asymmetric neutrino emission driven natal kicks. For a more detailed discussion about natal kicks, see Sect. 6 of Belczynski et al. (2016b).

The natal BH spin magnitude (aspin = (cJ)/(GM2), where J and M are BH angular momentum and mass, respectively) is estimated using one of three different prescriptions of angular momentum transport in stellar interiors, depending on the simulation performed. The most efficient angular momentum transport is adopted from Fuller & Ma (2019), Fuller et al. (2019) and we assign aspin = 0.01 for each BH natal spin in this scenario. The next scenario employs the standard Tayler-Spruit magnetic dynamo and efficient angular momentum transport from Spruit (1999) as adopted in the MESA evolutionary code; this model results in aspin ≈ 0.05 − 0.15 depending on progenitor star mass and metallicity. Finally, we also adopt an inefficient angular momentum transport driven by meridional currents (Ekström et al. 2012) as used in the Geneva evolutionary code; this results in aspin ≈ 0.9 for low-mass BHs (≲20 M) and aspin ≈ 0.2 for high-mass BHs (≳20 M). Full details of the three models are given in Belczynski et al. (2020).

The natal BH spin magnitude may be affected by tidal interactions between the stars in progenitor binaries of BH-NS mergers. Therefore, we allow for tidal spin-up. If a Wolf-Rayet (WR) star forms in a sufficiently close binary, it is subject to strong tidal interactions, which affect its rotation in comparison to single stellar evolution models (see Hotokezaka & Piran 2017, Qin et al. 2018). The resultant spin of a BH formed from such a star will be different from the spin of a BH that was formed either in isolation or in a wide binary. We adopt natal BH spin magnitudes formed from tidally affected WR stars from Belczynski et al. (2020):

a spin = e 0.1 ( P orb / P 0 1 ) 1.1 + 0.125 , $$ \begin{aligned} a_{\rm spin}=e^{-0.1(P_{\rm orb}/P_0 -1)^{1.1}}+0.125, \end{aligned} $$(2)

where Porb[s] is the orbital period and P0 = 4000 s. We apply this formula for Porb = 0.1 − 1.3 days. We assume that wider systems (with Porb > 1.3 days) do not experience significant tidal interactions. BHs originating from binaries with Porb < 0.1 days are assigned a spin magnitude of aspin = 1.0.

Table 2 contains list of our models.

Table 2.

Calculated models, Supernova engine: R/D–rapid/delayed.

2.2. BH-NS merger mass ejection and kilonovae

When an NS is disrupted in a BH-NS merger, the radioactive ejecta can emit a burst of electromagnetic radiation known as a kilonova. We do not consider other potential electromagnetic signals in this study (e.g., jets producing short gamma-ray bursts considered by Postnov et al. 2019).

During the merger there are two basic scenarios: (i) the NS is disrupted outside the BH event horizon (potential EM signal), or (ii) the NS is disrupted inside the BH event horizon (no EM signal). What happens depends on the BH spin, the NS equation of state (EoS), the mass ratio of two compact objects and the orientation of the NS orbit with respect to the BH spin plane. For each BH-NS binary system that we form in our population synthesis calculations, we consider all these parameters. We determine the ejecta mass from Kawaguchi et al. (2016):

M ej M NS , b = max { a 1 Q n 1 ( 1 2 C NS ) C NS 1 a 2 Q n 2 r ISCO ( a spin · cos ( i tilt ) ) + a 3 ( 1 M NS M NS , b ) + a 4 , 0 } , $$ \begin{aligned} \dfrac{M_{\rm ej}}{M_{\mathrm{NS},b}}&= \max \left\{ a_1 Q^{n_1} (1-2C_{\rm NS})C_{\rm NS}^{-1}\right.\nonumber \\&\,\,\, \left.- a_2 Q^{n_2} \tilde{r}_{\rm ISCO} \left(a_{\rm spin} \cdot \cos (i_{\rm tilt}) \right) + a_3 \left(1-\dfrac{M_{\rm NS}}{M_{\mathrm{NS},b}}\right) +a_4 ,0 \right\} , \end{aligned} $$(3)

where Q = MBH/MNS (MBH is BH mass and MNS is NS mass), itilt is the angle between BH spin and orbital angular momentum, and CNS = GMNS/(c2RNS) is the NS compactness parameter that depends on the EoS. MNS and RNS are neutron star mass and radius, respectively.

We use the MPA1 NS EoS, which is consistent with observational constraints on neutron star tidal deformation from the first NS-NS merger observation (see Fig. 1), as well as neutron star interior composition explorer (NICER) observations (Miller et al. 2019b). We employ (and this figure shows) mass-radius estimates for several EoS s drawn from 1. The maximum NS mass for the MPA1 equation of state is ∼2.46 M. Rather than treat objects with masses between 2.46 M and 2.5 M as black holes (≲1% of NSs in BH-NS mergers in our models), we assume a constant mass-radius relation for objects in this mass range and we treat them as neutron stars. A 2.46 M NS has 11.3 km radius, and then all NSs above this mass (up to 2.5 M) have 11.3 km radius.

thumbnail Fig. 1.

Comparison of the NS mass-radius relation of four different equations of state (EoS s) with the LIGO/Virgo constraints from NS-NS merger GW170817 (90% credible limits shown by green and blue boxes: Abbott et al. 2018a). We employ the MPA1 EoS, which reaches ∼2.5 M (our adopted maximum NS mass) and agrees with the LIGO/Virgo constraints. Other relations are WFF1 (Wiringa et al. 1988), APR4 (Akmal et al. 1998), and DD2 (Alvarez-Castillo et al. 2016).

The parameters a1 − a4, n1, and n2 are discussed in Kawaguchi et al. (2016). MNS, b is the total baryonic mass of an NS. We use the relation between MNS and MNS, b for nonrotating NS from Gao et al. (2020) and Gupta et al. (2020), lowered by 1.8% to account for the relative error of NS baryonic mass estimation (a lower MNS, b provides more realistic ejecta masses). r ISCO $ \tilde{r}_{\mathrm{ISCO}} $ is the innermost stable circular orbit radius normalized by MBH (e.g., from Foucart et al. 2018):

r ISCO = r ISCO G M BH c 2 = 3 + Z 2 sign ( a spin ) ( 3 Z 1 ) ( 3 + Z 1 + 2 Z 2 ) , $$ \begin{aligned} \tilde{r}_{\rm ISCO}=\dfrac{r_{\rm ISCO}}{GM_{\rm BH} c^{-2}}=3+Z_2-\mathrm{sign}(a_{\rm spin})\sqrt{(3-Z_1)(3+Z_1+2Z_2)}, \end{aligned} $$(4)

where Z1 and Z2 are functions of BH spin, and Z 1 =1+ (1 a spin 2 ) 1/3 ( (1+ a spin ) 1/3 + (1+ a spin ) 1/3 ) $ Z_1=1+(1-a_{\rm spin}^2)^{1/3}((1+a_{\rm spin})^{1/3} +(1+a_{\rm spin})^{1/3}) $ and Z 2 = 3 a spin 2 + Z 1 2 $ Z_2=\sqrt{3a_{\mathrm{spin}}^2+Z_1^2} $. However, we use the projection of BH spin on the binary orbital momentum direction to calculate r ISCO $ \tilde{r}_{\mathrm{ISCO}} $ ( r ISCO ( a spin · cos ( i tilt ) ) $ \tilde{r}_{\mathrm{ISCO}}(a_{\mathrm{spin}} \cdot \cos(i_{\mathrm{tilt}})) $).

Kilonova emission remains an active area of research and the peak luminosity depends on a wide range of factors concerning the properties of the ejecta, such as their composition, their morphology, and the viewing angle (see Wollaeger et al. 2019, 2018; Zhu et al. 2018; Fontes et al. 2020; Korobkin et al. 2021). For this paper, we use the following parameterized formula for the o-band (1260 − 1360 nm) peak luminosity:

L peak = f × 10 40 ( M ej 0.01 M ) 0.2 ( v 0.1 c ) 0.2 erg s 1 , $$ \begin{aligned} L_{\rm peak}=f \times 10^{40} \left( \dfrac{M_{\rm ej}}{0.01\,{M}_{\odot }}\right)^{0.2} \left( \dfrac{{ v}}{0.1 c}\right)^{0.2} \mathrm {erg\ s^{-1}}, \end{aligned} $$(5)

where f is varied to match the wide range of current results: 0.3 ≤ f ≤ 4.5. We assume that f has uniform distribution in this range. For ejecta velocity, we use the formula from Kawaguchi et al. (2016):

v / c = 0.01533 Q + 0.1907 . $$ \begin{aligned} { v}/c=0.01533 Q + 0.1907. \end{aligned} $$(6)

We calculate apparent flux from the kilonova for each merger event with mass ejection. We assume that mass ejection is isotropic. Therefore, the flux can be calculated as:

F = L peak / ( 1 + z ) 2 4 π D L 2 , $$ \begin{aligned} F=\dfrac{L_{\rm peak}/(1+z)^2}{4 \pi D_L^2}, \end{aligned} $$(7)

where z is redshift and DL is the luminosity distance of a given merger.

This flux corresponds to the apparent o-band magnitude on Earth. The apparent kilonova magnitude can be calculated using Pogson’s equation:

M kn = m Vega 2.5 log 10 ( F F Vega ) , $$ \begin{aligned} M_{\rm kn}=m_{\rm Vega}-2.5 \log _{10} \left( \frac{F}{F_{\rm Vega}} \right), \end{aligned} $$(8)

where mVega = 0.026 mag is the apparent Vega magnitude as seen from Earth, FVega = 2.18072 × 10−8 W m−2 is flux from Vega, and F is flux from the kilonova.

2.3. Cosmology

We assume a flat Universe with ΩM = 0.3, ΩΛ = 0.7, Ωk = 0 (Ade et al. 2016), and Hubble constant H0 = 68.7 km s−1 Mpc−1. This gives the age of Universe t0 ∼ 13.7 Gyr (Planck Collaboration VI 2020; Hinshaw et al. 2013).

We assume that binary stars (Population I and II) are formed in the redshift range z = 0 − 15. The star formation rate density history is adopted from Madau & Fragos (2017):

sfrd ( z ) = 0.015 ( z + 1 ) 2.6 1 + ( ( z + 1 ) / 3.2 ) 6.2 M Mpc 3 yr 1 . $$ \begin{aligned} \mathrm{sfrd}(z)=0.015 \frac{(z+1)^{2.6}}{1+((z+1)/3.2)^{6.2}} \,{M}_{\odot }\, \mathrm{Mpc}^{-3}\ \mathrm{yr}^{-1}. \end{aligned} $$(9)

We adopt an average cosmic metallicity evolution from Madau & Fragos (2017) and we adopt a Gaussian distribution (with 0.5 dex sigma) of metallicity around the average at each redshift (Belczynski et al. 2020). After evolution, some binaries form BH-NS systems and some of these systems merge. We record their merger redshifts (and corresponding distances). We also calculate an intrinsic BH-NS merger rate density and an intrinsic merger rate as a function of redshift for all of our models, following the method described in detail by Belczynski et al. (2016b).

2.4. Result presentation

The results of our simulations are presented in the context of intrinsic properties of binary systems, as well as the properties of the subset of systems likely to be detected by LIGO/Virgo. We mostly focus on low redshift systems (z < 1) because BH-NS mergers will not be detected by LIGO/Virgo above redshift z = 1 even at full design sensitivity. We also present results based on the current sensitivity limits of the LIGO/Virgo O3 science run.

We present the intrinsic fraction of BH-NS mergers that have compact object or objects in the FMG, which we estimate from:

η 1 = R d fmg , z < 1 R d z < 1 , $$ \begin{aligned} \eta _1= \frac{Rd_{\mathrm{fmg},z < 1}}{Rd_{z < 1}}, \end{aligned} $$(10)

where Rdfmg, z < 1 (Gpc−3 yr−1) is the intrinsic merger rate density for BH-NS mergers that occur for redshifts z < 1 and that have two, one, or no objects within the FMG (assumed 2.1 − 5 M range, see Sect. 2.1). If there is only one compact object within the FMG, we identify whether it is an NS or a BH (using our 2.5 M maximum mass limit for an NS). Rdz < 1 (Gpc−3 yr−1) is the intrinsic merger rate density for all z < 1 BH-NS mergers.

To present results for mass ejection in mergers of BH-NS binaries, we use a fraction of intrinsic merger rate density. The intrinsic fraction of z < 1 BH-NS mergers with some mass ejection is found from:

η 2 = R d Mej , z < 1 R d z < 1 , $$ \begin{aligned} \eta _2= \frac{Rd_{\mathrm{Mej},z < 1}}{Rd_{z < 1}}, \end{aligned} $$(11)

where RdMej, z < 1 (Gpc−3 yr−1) is the intrinsic merger rate density for BH-NS mergers that occur at redshifts z < 1 and has an ejecta mass Mej > 0.001 M or Mej > 0.01 M, and Rdz < 1 (Gpc−3 yr−1) is the intrinsic merger rate density for all z < 1 BH-NS mergers.

Using simple mass scaling, we assess the detectability of BH-NS mergers in gravitational waves by LIGO/Virgo. We assume that BH-NS mergers are only detectable if the signal-to-noise ratio (S/R) in one LIGO/Virgo detector is larger than 8, and we estimate a fiducial S/R from:

S / R = 8.0 120 Mpc D L ( M chirp 1.2 M ) 5 / 6 , $$ \begin{aligned} S/R=8.0\dfrac{120\,\mathrm {Mpc}}{D_L} \left(\frac{M_{\rm chirp}}{1.2\,{M}_{\odot }}\right)^{5/6}, \end{aligned} $$(12)

where Mchirp = (m1m2)3/5(m1 + m2)−1/5 is the chirp mass of the BH-NS system, and m1 and m2 are BH-NS components masses. In this formula, we assume that the current O3 LIGO/Virgo detection range for an NS-NS merger with a typical chirp mass of 1.2 M is 120 Mpc (Abbott et al. 2018b), and we average over the source orientation and sky location. This gives us a detection rate of BH-NS mergers by LIGO/Virgo: Rligo yr−1. The fraction of BH-NS mergers that have compact objects in the FMG that are detectable by LIGO/Virgo is calculated from:

η 3 = R fmg , ligo R ligo , $$ \begin{aligned} \eta _3= \frac{R_{\rm fmg,ligo}}{R_{\rm ligo}}, \end{aligned} $$(13)

where Rfmg, ligo (yr−1) is the LIGO/Virgo detection rate of BH-NS mergers that have two, one or no objects within the FMG.

We also calculate the fraction of LIGO/Virgo detectable BH-NS mergers (S/R > 8) that show some mass ejection (Mej) during the merger process as:

η 4 = R Mej , ligo R ligo , $$ \begin{aligned} \eta _4= \frac{R_{\rm Mej,ligo}}{R_{\rm ligo}}, \end{aligned} $$(14)

where Rligo (yr−1) is the LIGO/Virgo detection rate of BH-NS mergers, and RMej, ligo (yr−1) is the detection rate of BH-NS mergers with Mej larger than a specified value (i.e., Mej > 0.001 M).

BH-NS mergers that are detectable as kilonovae are defined as follows: first, they need to be detectable by LIGO/Virgo (S/R > 8) and then they need to have enough mass ejecta escaping with appropriate speed during the merger to produce emission bright enough to be detected by a given telescope (see Eq. (5)). For each BH-NS merger with mass ejection (and kilonova), we use its apparent magnitude (Mkn, see Eq. (8)) and compare it with the threshold sensitivity of three different telescopes that are employed in the search for kilonovae.

We consider a range of telescopes, from small instruments that can cover a large fraction of the sky, such as ATLAS (2 telescopes, each with a diameter of 0.5 m) with a maximum reach of Mmax = 19.5 mag (exposure time 0.5 min), to the medium-size Canada-France-Hawaii Telescope (CFHT; 1 telescope with a diameter of 3.6 m) with 24.1 mag limit (exposure time 60 min), and the large Subaru telescope (1 telescope with a diameter of 8.2 m) with a limit of 26.0mag (exposure time 4 − 8 min). These exposure times are adopted from the official cameras’ specifications as characteristic times at certain bands.

ATLAS’s detection limit is given directly for o-band, while the CFHT and Subaru limits are adopted for J-band (WIRCam camera limit, band centered at 1220 nm) and Ic filter (Suprime-Cam camera limit, see Miyazaki et al. 2002), filter centered at 806 nm), respectively. If a kilonova is brighter than the threshold for a given telescope:

M kn < M max $$ \begin{aligned} M_{\rm kn} < M_{\rm max} \end{aligned} $$(15)

then we call it a detectable EM/kilonova counterpart to a given LIGO/Virgo BH-NS merger signal. Then we calculate the fraction of detectable BH-NS mergers as kilonova that are brighter than chosen telescope brightness threshold (see Eq. (15)) in all LIGO/Virgo detectable BH-NS mergers from:

η 5 = R kn , ligo R ligo , $$ \begin{aligned} \eta _5= \frac{R_{\rm kn,ligo}}{R_{\rm ligo}}, \end{aligned} $$(16)

where Rkn, ligo (yr−1) is the LIGO/Virgo detection rate of BH-NS mergers with detectable kilonovas, and Rligo (yr−1) is the LIGO/Virgo detection rate of BH-NS mergers.

It should be noted that the above scheme does not account for any localization issues that may arise during the search for kilonovae associated with LIGO/Virgo sources (e.g., Nissanke et al. 2013; Gomez et al. 2019). This approximation serves only as a guide to inform us whether any given kilonova is bright enough to be detected (with some typical exposure time) by a given telescope, if this telescope was pointed right at the kilonova. Additionally, we use peak brightness, so this is a very optimistic approximation of kilonova detectability.

3. First mass gap objects in BH-NS mergers

First we study the BH-NS merger compact object mass distribution (z < 1) for several evolutionary models (M230, M233, M280, and M283) using both rapid (Fig. 2) and delayed (Fig. 3) supernova engine prescriptions for remnant masses. As expected, the FMG is clearly visible in models that employ a rapid supernova engine (e.g., M230 and M233). Therefore, for these models, we do not expect BH-NS mergers with compact objects of comparable mass. However, for models that employ a delayed supernova engine (e.g., M280 and M283), there is no mass gap between NSs and BHs, and we predict some BH-NS mergers with comparable mass compact objects. This is important in context of mass ejection in BH-NS mergers (see Sect. 6).

thumbnail Fig. 2.

Mass distribution of compact objects in BH-NS mergers in the low-redshift Universe (z < 1) for models that employ a rapid supernova engine (M230, M233). There is a very clear mass gap between NSs and BHs. It should be noted that our data become noisy for masses above ≳20 M. However, apparent decline in merger rate density for these higher masses is real and expected in our models.

In Figs. 2 and 3 we also clearly see the effect of natal kicks on the merger rate density of BH-NS systems. For models with decreased natal kicks (e.g., M230 and M280) the rates are higher by about an order of magnitude than for models with high natal kicks (e.g., M233 and M283). This is a consequence of BH-NS progenitor binaries being disrupted more easily in models in which we employed higher natal kicks. One can also see some noisy features that appear at high mass ends. Those are the results of low statistics.

thumbnail Fig. 3.

Mass distribution of compact objects in BH-NS mergers in the low-redshift Universe (z < 1) for models that employ a delayed supernova engine (M280, M283). There is no mass gap between NSs and BHs. It should be noted that our data become noisy for masses above ≳20 M, although apparent decline in merger rate density for these higher masses is real and expected in our models.

In Table 3 we list fractions of BH-NS mergers that have both compact objects within the FMG (FMGboth), only a BH in the gap (FMGBH), only an NS in the gap (FMGNS), or no components in the gap (FMGnone). This table contains BH-NS mergers that take place for redshifts with z < 1 and the fractions are calculated using merger rate densities (see Eq. (10)). For models that the employ rapid supernova engine, there are no BH-NS mergers with objects in the FMG (assumed 2.1 − 5 M range). For models that employ the delayed supernova engine remnant prescription, about half or a little more than half (0.576 − 0.708) of BH-NS mergers do not have any mass gap objects. With these models, those that employ high natal kicks independent of compact object mass (M283, M383, and M483) show a significant fraction of BH-NS mergers with both compact objects within the mass gap (0.197 − 0.212). Since all compact objects have an equal probability of getting a high natal kick and disrupting the progenitor binary, the mass function of these compact objects falls steeply off with mass (approximately following the initial mass function of the stars); see the black (NS) and green (BH) lines in Fig. 3. That results in a significant fraction of BH-NS mergers with both compact objects within the FMG. In contrast, for our model in which natal kicks decrease with compact object mass (M280), this fraction is very small (0.036) because many BH-NS systems with heavy (∼10 M) BHs outside the FMG are not disrupted by the kick; see the blue line in Fig. 3. Depending on the model, significant fractions of BH-NS mergers (0.211 − 0.256) have one compact object (whether it is an NS or a BH) within the FMG. Although we differentiate between NSs and BHs in the mass gap in Table 3, we note that this division is arbitrary as we simply assume that any compact object above 2.5 M is a BH.

Table 3.

Intrinsic fraction of (z < 1) BH-NS mergers η1 in which both (FMGboth) components, only a BH (FMGBH) component, only an NS (FMGNS) component, or no (FMGnone) components are within the FMG.

In Table 4 we list fractions of BH-NS mergers with compact objects within the FMG, but only for systems that are detectable by LIGO/Virgo (S/R > 8; see Eq. (13)). The detectable LIGO/Virgo population is predicted to have a significant fraction of mergers (0.199 − 0.334) with at least one compact object within the mass gap for delayed supernova engine models (M280, M283, M383, and M483). A small but noticeable fraction (0.025 − 0.186) of BH-NS mergers are found to have both compact objects within the gap. Compared with the intrinsic population (see Table 3), the LIGO/Virgo detectable population has a smaller fraction of BHs and a larger fraction of NSs in the gap. This simply reflects the fact that LIGO/Virgo can detect heavier objects from larger distances and that light BHs within the gap tend to be less represented in the LIGO/Virgo detectable population, while NSs need to be heavy to be within the gap. Thus, NSs stand out more in this population.

Table 4.

Fraction of LIGO/Virgo (S/R > 8) BH-NS systems η3 where both (FMGboth) components, only a BH (FMGBH) component, only an NS (FMGNS) component, or no (FMGnone) components are within the FMG.

The fraction of FMG+FMG mergers among the entire LIGO/Virgo detectable BH-NS merger population is significantly smaller for the low BH natal kick model (M280: ∼0.04) than for high kick models (M283, M383, and M483: ∼0.2; see Table 4). FMG+FMG systems tend to survive high natal kicks more often than other BH-NS merger progenitors because of their specific evolutionary history that typically leads to small separations before supernova explosions (see Sect. 4.3). Therefore, models with high natal kicks contain higher fractions of FMG+FMG systems, which are mergers that are typically accompanied by mass ejection, (see Sect. 5) and thus show kilonova emission (see Sect. 6). However, these models have small BH-NS merger rates (see Sect. 7).

We find that BH-NS mergers may still be detectable at distances as far as z ∼ 0.1 with LIGO/Virgo during O3 (S/R > 8). Due to GW selection effects, the most distant mergers would be the most massive ones, for example: 2.2 M NS + 30.8 M BH.

4. Mass ratio of BH-NS systems

In Fig. 4 we show the mass ratio distribution (mass of an NS to a BH: q = MNS/MBH) for models with a rapid supernova engine (M230 and M233) and models with delayed supernova engine (M280 and M283). This figure includes the population of BH-NS mergers at redshifts z < 1 and shows the intrinsic mass ratio expressed by the merger rate density in bins of Δq = 0.01.

thumbnail Fig. 4.

Mass ratio distribution for BH-NS mergers in the low-redshift Universe (z < 1). For each BH-NS system, we calculated the mass ratio q = MNS/MBH ∈ [0, 1]. Values indicate merger rate densities estimated from the systems within each q bin, that is to say dRate dt−1dV−1dq−1Δq, with Δq = 0.01. For transparency purposes, we use simple merger rate density as the vertical axis label. Most mergers have rather small mass ratios q ∼ 0.1 − 0.2, but there are also extreme mass ratio systems such as 56 M BH + 1.3 M NS (q = 0.023). For models with the delayed supernova engine (M280 and M283), there is a secondary peak for high mass ratios q > 0.4, which consists mostly of double compact objects with both components in the FMG (see Sect. 4).

As seen immediately from Fig. 4, the majority of BH-NS mergers have small mass ratios q ∼ 0.1 − 0.2 independent of model assumptions. This comes from the fact that NSs have typical masses of 1 − 2 M and the majority of BHs found in BH-NS mergers in our models have a typical mass of ∼10 M (see Figs. 2 and 3). Since the rapid supernova engine does not populate the mass gap, when adopting this model, we thus predict no mergers with a mass ratio larger than q ∼ 0.4. However, for models with a delayed supernova engine, the mass ratio distribution shows a secondary peak for q ∼ 0.5 − 0.8. These are BH-NS mergers with both compact objects within the FMG, as such systems must have q ≳ 0.4.

In Fig. 5 we show cumulative mass ratio distribution for LIGO/Virgo detectable (S/R > 8) BH-NS mergers. In Table 5 we list fractions of LIGO/Virgo detectable BH-NS mergers (S/R > 8) that have a mass ratio smaller than 1/50, 1/30, 1/20, and 1/10 for all our models. Additional models that include BH-NS mergers with input physics assumptions that are different to those used in this study may be found online2 under the tab “Download/2020: Double Compact Objects/Belczynski et al. (2020)”.

thumbnail Fig. 5.

Cumulative distribution of LIGO/Virgo detectable (S/R > 8) BH-NS mergers in terms of mass ratio (see also Table 5).

Table 5.

Fraction of LIGO/Virgo detectable (S/R > 8) BH-NS mergers with extremely small ratios and extremely small mass ratios.

Our results show that it is possible to create extreme mass ratio systems q < 0.05. For example, the most extreme systems detectable by LIGO/Virgo (S/R > 8) in our models are: M230 q = 0.021 (MBH = 55.0 M, MNS = 1.2 M), M233 q = 0.023 (MBH = 49.0 M, MNS = 1.1 M), M280 q = 0.025 (MBH = 52.5 M, MNS = 1.3 M), M283 q = 0.029 (MBH = 29.9 M, MNS = 0.9 M), M383 q = 0.022 (MBH = 59.5 M, MNS = 1.3 M), and M483 q = 0.021 (MBH = 60.1 M, MNS = 1.3 M).

In the next subsection we discuss in detail the formation of one extreme mass ratio system (38.9 M + 1.3 M; Sect. 4.1), as infrequent as it may be. We follow with a detailed description of the formation of a BH-NS merger with a typical mass ratio (14.7 M + 1.8 M; Sect. 4.2), and finish with a description of a BH-NS merger with comparable mass components, both within the FMG (3.5 M + 2.1 M; Sect. 4.3).

4.1. Extreme mass ratio BH-NS: 38.9 + 1.3 M merger

In this study we use the terms primary and secondary component for binary components with higher and lower zero age main sequence (ZAMS) mass, respectively. Our calculations show that BH-NS mergers (independent of their mass ratio) follow a similar evolutionary history. In Fig. 6 we show the major stages of evolution of a massive binary (primary M1 = 91.7 M, secondary M2 = 9.5 M) that ultimately produces an extreme mass ratio BH-NS merger (MBH = 38.9 M, MNS = 1.3 M: q = 0.033). Below we describe the major evolutionary stages.

thumbnail Fig. 6.

Evolutionary route leading to the formation of the 38.9 M + 1.3 M (q = 0.033) BH-NS merger for metallicity Z = 0.0001 (model M280). The evolution from ZAMS to the formation of the BH-NS system takes ∼30.4 Myr, and then it takes ∼2.29 Gyr for the system to in-spiral (angular momentum loss due to the emission of gravitational waves), leading to the final merger of the two compact objects. (CHeB: core helium burning star, He star: naked helium star, HG: Hertzsprung gap star). No tidal spin-up is expected due to large orbital separations when binary components are naked helium stars: a ∼ 25 R (Porb = 2.1d; star A) and a ∼ 200 R (Porb = 60d; star B; see Sect. 2.1).

The evolution starts with a very massive primary (MA = 91.7 M) and an intermediate mass secondary (MB = 9.5 M) on a very wide orbit (semimajor axis a = 7400 R) with a small eccentricity (e = 0.09) Next, the primary becomes a core helium burning star while the secondary is still on the main sequence. Then, a CE phase is initiated by expansion of the primary; the primary becomes a naked helium star after H-rich envelope ejection. Afterwards, the primary experiences a pair-instability pulsation supernova mass loss and immediately collapses directly to a BH. The secondary becomes a Hertzsprung Gap (HG) star and, 14 kyr later, it initiates a stable Roche lobe overflow mass transfer that ends when the secondary is in the core helium burning stage, making a blue loop on the H-R diagram (radius decreases). Then, the Roche-lobe overflow restarts while the secondary is still in the core helium burning stage and it is moving redwards on the H-R diagram loop (radius increases). The Roche lobe overflow ends when the secondary loses most of its H-rich envelope. Afterwards, the secondary becomes a naked helium star and explodes in a Type Ib (core-collapse) supernova, forming an NS.

A highly eccentric BH-NS binary with an extreme mass ratio is formed. After more than 2 billion yr of tightening the orbit, the merger of the BH with the NS leads to a burst of gravitational waves, but there is no mass ejection during the merger (no kilonova).

There is nothing specifically unusual in the formation of extreme mass ratio BH-NS mergers. Similar evolutionary routes are responsible for nonextreme BH-NS mass ratio systems (see Klencki et al. 2018). BH-NS progenitors experience two major binary interactions (CE and RLOF) and are subject to one strong supernova explosion forming an NS. The natal kick (243 km s−1 in the case of the example shown in Fig. 6) produced during the supernova can make the BH-NS system highly eccentric. This reduces the time to final merger, allowing such systems to merge within Hubble time and to be detected by LIGO/Virgo. However, for extreme mass ratio systems there is no mass ejection during the merger and therefore these mergers are not expected to produce EM (e.g., kilonova or short GRB) signals (see Sect. 6).

4.2. Typical mass ratio BH-NS: 14.7 + 1.8 M merger

There are several variations to the formation of BH-NS mergers in respect to the evolution presented in Fig. 6. Here we pick one such variation to show the formation of a typical mass ratio (q = 0.12) BH-NS merger in our simulations.

The formation of a 14.7 + 1.8 M BH-NS merger at metallicity Z = 0.0005 (model M230) involves the evolutionary stages, which we describe below.

The evolution starts with a massive primary (MA = 35.8 M) and an intermediate mass secondary (MB = 20.6 M) on a wide (a = 6641 R) and eccentric orbit (e = 0.56). At some point, the primary becomes an asymptotic giant branch (AGB) star, while the secondary is still on the main sequence. As a result, the CE phase is initiated by the primary, which loses about half of its mass. Then, the primary undergoes a direct collapse (only a 0.15 M mass loss in neutrinos, no supernova explosion) to a BH and we assume no natal kick. Next, the secondary becomes a core helium burning star, initiates a second CE episode and becomes a naked helium star. Later the secondary explodes in a Type Ib/c (core-collapse) supernova, forming an NS with a high natal kick (229 km s−1). Subsequently, an eccentric (e = 0.90) BH-NS binary with a moderate mass ratio forms after 10.5 Myr of evolution since ZAMS. After 2.7 Gyr, the merger of the BH with the NS leads to a burst of gravitational waves but without mass ejection or kilonova (Fig. 7).

thumbnail Fig. 7.

Evolution of two massive stars that leads to the formation of the moderate mass ratio (q = 0.12) BH-NS (14.7 + 1.8 M) merger (model M230). The top panel shows the mass evolution of the primary and secondary stars, while the bottom panel shows the changes in the semimajor axis and eccentricity. Evolution from ZAMS to the BH-NS system formation takes 10.5 Myr, and then it takes another 2.7 Gyr for the two compact objects to merge and produce gravitational wave radiation that is potentially detectable by LIGO/Virgo. No mass ejection (no kilonova or short GRB) is expected in this case.

This typical q = 0.12 system was formed from a ZAMS binary with comparable mass components: qZAMS = 0.58, and it evolved through two CE phases. Both components evolve along similar tracks. They form remnants shortly after the CE they initiated. No tidal spin-up before the BH formation is expected: a = 557 R and Porb = 257d. Both stars lose a significant fraction (≳50%) of their mass in CE events, and the secondary component ejects some extra mass from the binary during a Type Ib/c supernova that forms a heavy NS. A high natal kick at the NS formation leads to the formation of a very eccentric (e = 0.9) BH-NS binary that can merge within the Hubble time, despite its rather large separation (a ∼ 40 R).

4.3. Comparable mass BH-NS: 3.5 + 2.1 M merger

Comparable mass BH-NS systems, by construction, have a high-mass NS and a low-mass BH, and typically one or both of the components have mass within the FMG. Therefore, such systems appear only within delayed supernova models. Below, we present the BH-NS merger with both components within the FMG (3.5 M BH +2.1 M NS; q = 0.6) that has formed at metallicity Z = 0.0015 in model M280. Firstly, the ZAMS binary is formed with almost equal-mass components: 18.80 M and 18.66 M at a moderately wide (a = 1100 R) and eccentric (e = 0.36) orbit. Later, when both components become core helium burning stars at a similar time, stable RLOF starts from the primary to the secondary. As a result, the primary component becomes a naked helium star after losing its H-rich envelope. Then, the CE is initiated by the secondary component, which is still a core helium burning star. After the CE, both components are naked helium stars. However, this state is not permanent: another stable RLOF takes place from the primary (evolved naked helium star) to the secondary component (∼2.0 M mass loss). As a result the primary component explodes as a stripped supernova (∼2.6 M mass loss) with a 228 km s−1 natal kick forming an NS, and the binary becomes eccentric (e = 0.24). The secondary component supernova explosion forms a BH and a 217 km s−1 natal kick reduces system eccentricity (e = 0.06). Then, a close (a = 3.6 R) BH-NS binary forms on an almost circular orbit after ∼12 Myr evolution since ZAMS. After another 0.6 Gyr, the BH-NS merger occurs.

In this case, the NS forms ∼0.2 Myr before the BH forms from the primary star due to early mass ratio reversal between the primary and the secondary (first stable RLOF). Figure 8 shows the evolution of the binary progenitor of this system in more detail. The small orbital separation prior to the BH formation (a = 2.3 R, Porb = 0.11d) causes the WR star progenitor of the BH to be spun up due to tidal interactions. The BH spin becomes very large: aspin = 0.983 in comparison to aspin = 0.092 expected from single stellar evolution without tidal interactions. Due to the high NS mass and its compactness, no mass ejection is expected even for the high spin and low mass of a BH.

thumbnail Fig. 8.

Evolution of binary system that forms a 3.5 + 2.1  M BH-NS merger (mass ratio of q = 0.6). We note that both compact objects are within the FMG and we expect that this merger will not lead to any mass ejection due to a high NS mass.

5. Mass ejection in BH-NS mergers

Some BH-NS mergers may produce an electromagnetic counterpart. For this to take place, the NS must disrupt outside the BH event horizon. In this section we discuss the amount of mass (if any) that is ejected during the BH-NS merger process.

Table 6 shows the fraction of BH-NS mergers in our intrinsic population (z < 1; Eq. (11)) with any mass ejecta (Mej > 0.001 M) and with significant mass ejecta (Mej > 0.01 M) for all our models. The corresponding distribution and reverse cumulative3 distributions are shown in Figs. 9 and 10, respectively.

thumbnail Fig. 9.

Intrinsic distribution of ejecta mass in BH-NS mergers with z < 1. The structures at Mej ≳ 0.03 M are shot noise.

thumbnail Fig. 10.

Intrinsic reverse cumulative distribution of ejecta mass (> Mej) in BH-NS mergers with z < 1.

Table 6.

Intrinsic fraction of BH-NS mergers (z < 1) η2 with any mass ejecta (Mej > 0.001 M) and with significant mass ejecta (Mej > 0.01 M).

The lowest ejecta mass found for BH-NS mergers produced by our simulations and calculated with equation 3 is Mej = 0.0002 M, but we use Mej = 0.001 M as the threshold for marking systems with any mass ejection. For our models that eject mass, the number of systems with Mej < 0.001 M is negligible. The largest ejecta mass is Mej ∼ 0.07 M, which is consistent with the results from Rosswog (2005), Kyutoku et al. (2015), Barbieri et al. (2021). The ejecta mass estimated for the first confirmed NS-NS merger, GW170817: ∼0.0002 − 0.03 M (Côté et al. 2017; Abbott et al. 2017) is in the range of our estimates.

For high BH spins (model M383; Geneva inefficient angular momentum transport), we obtain the maximum fraction of BH-NS mergers with any ejecta among our models: η2 = 0.087 (see Eq. (11)). This result is more pessimistic than earlier population synthesis predictions that showed the fraction of BH-NS mergers for rapidly spinning BHs with mass ejecta is as high as 0.4 (Belczynski et al. 2008b). This difference comes from the updated numerical simulations of ejecta mass, better constraints on the EoS, and an improved understanding of input stellar and binary physics. However, it should be noted that for low BH spins, our results are consistent with previous predictions that showed a 0.01 fraction of BH-NS mergers with mass ejecta (Belczynski et al. 2008b). For comparison, our reference model with low BH spins (model M280: efficient MESA angular momentum transport) shows that the fraction of BH-NS mergers with any mass ejecta is comparable: η2 = 0.006 (see Table 6).

The rapid supernova engine models (M230 and M233) produce BH-NS mergers with low-mass NSs (≲2 M) and show the larger fraction of systems with mass ejecta than the corresponding delayed supernova engine models (M280 and M283). In the delayed models, we also allow for the formation of heavy NSs (> 2 M), which are more compact than lower-mass NSs (see Fig. 1) and thus are harder to disrupt.

BH spin plays a significant role in setting the size of the event horizon and this regulates whether NS disruption can produce any mass ejecta (Foucart et al. 2018; Zappa et al. 2019; Sedda 2020). For rapidly rotating BHs (small event horizons), mass ejection is found in a small fraction of BH-NS mergers (M383) and the fraction decreases for slowly spinning BHs (M230, M233, M280, and M283). There are still some BH-NS mergers that generate mass ejection for the (almost) nonspinning BH model (M483).

Table 7 shows only those BH-NS mergers with any mass ejecta (Mej > 0.001 M; see Table 6) and we subdivide these to show the contribution of mergers in which none, one (BH or NS), or both components are within the FMG. These are, therefore, intrinsic fractions for the sub-population of BH-NS mergers (with any mass ejection) within z < 1. Models with a rapid supernova engine (M230 and M233) have both merger components always outside the FMG. However, for our reference model (M280) and two other models with a delayed supernova engine (M283 and M483) the fraction of BH-NS mergers with mass ejecta having both components in the FMG is very high ( ∼ 0.9).

Table 7.

Intrinsic BH-NS merger population with any mass ejecta (Mej > 0.001 M, z < 1).

Because binaries with both objects in the FMG correspond to more favorable mass ratios for tidal disruption outside the BH horizon, the fraction of binaries in the FMG that produce mass ejecta is high for delayed supernova engine models.

In the most extreme case of our reference model (M280), we find that the fraction of BH-NS mergers with any mass ejecta that have both components within the FMG is 0.959. These BH-NS systems have evolved through spin-up of a WR star and BHs have a very large spin that allows for easy mass ejection (aspin > 0.9). The lowest mass BHs (i.e., in the FMG) are naturally selected for systems with mass ejecta as they have the smallest event horizon. The heaviest NSs (i.e., in the FMG) are also favored since in order to spin up the BH progenitor (a WR star), either an NS forms first (so it must be heavy; mass ratio reversal), or the NS progenitor is a WR star at the time when the BH progenitor is a WR star (so both stars evolve at almost similar timescales and must be similar in mass). Additionally, this model allows for fallback-moderated natal kicks, so heavy NSs receive smaller kicks and they have an increased chance of surviving in a binary after a supernova explosion than the lighter NSs.

For the model with high initial BH spins (M383), the fraction of mergers with mass ejecta and both compact objects in the FMG is smaller (0.365), as specific processes (i.e., tidal spin up) are not the only way to produce mass ejecta.

In Table 8 we list fractions of BH-NS mergers with compact objects in and out of the FMG, but only for systems that have any mass ejection (Mej > 0.001 M) and are detectable by LIGO/Virgo (S/R > 8; see Eq. (12)). These fractions are calculated using merger rates (see Eq. (14)). This table shows how likely it is for a BH and an NS to be within the FMG for a LIGO/Virgo BH-NS merger that may have detectable kilonova (see Sect. 6).

Table 8.

LIGO/Virgo detectable BH-NS mergers (S/R > 8) with any mass ejecta Mej > 0.001 M.

The results are qualitatively similar to those we found for the intrinsic population of BH-NS mergers presented above in Table 7. The small to moderate quantitative differences arise from the fact that our adopted criteria for the LIGO/Virgo detectability set the horizon redshift for (light) BH-NS mergers with mass ejection to only zhor ∼ 0.07 (Table 8). For comparison, the horizon redshift for the overall BH-NS merger LIGO/Virgo detectable population is zhor ∼ 0.1. In contrast, the redshift cut that we employed in Table 7 is z = 1.

6. Electromagnetic counterpart

BH-NS mergers in which the NS is disrupted outside the BH event horizon produce mass ejecta and accretion disks that may produce EM counterparts. In this section we discuss our results on kilonova emission associated with BH-NS mergers with mass ejection. For LIGO/Virgo detectable events (S/R > 8), we calculate kilonova luminosities for each BH-NS merger with some mass ejecta (Mej > 0.001 M) following the scheme outlined in Sect. 2.2. For each merger event we employ the merger distance from our synthetic Universe (see Sect. 2.3).

Figure 11 shows the distribution of brightness in o-band (1260 − 1360 nm) of BH-NS mergers with kilonova emission for all our models. The distribution of kilonova brightness has a very large range (Mkn ∼ 15 − 28 mag), but dim kilonovae dominate (Mkn ∼ 24 mag). Since such events can be detected only at small redshifts (see Sect. 5), we do not take into account k-correction. The large brightness range comes from the combination of distance range for BH-NS mergers that are detectable by LIGO/Virgo during O3 (z ≲ 0.07; luminosity distance: ≲320 Mpc) and the fact that ejecta mass is found in a wide range of values (Mej ∼ 0.0002 − 0.07 M). We find that most kilonovae are dim because most events are found at distances (average redshift: zave ∼ 0.04; average luminosity distance of ∼180 Mpc) close to the LIGO/Virgo O3 detection horizon for BH-NS mergers (z ∼ 0.07; where most of the searched volume is located). Additionally, the ejecta mass found in our simulations is rather small; typically Mej < 0.03 M (see Fig. 9).

thumbnail Fig. 11.

Kilonova brightness distribution for LIGO/Virgo detectable (S/R > 8) BH-NS mergers.

Kilonova brightness varies with different model assumptions. In Table 9, for all our models, we show the fraction of BH-NS mergers with detectable kilonova among LIGO/Virgo detectable events (S/R > 8; see Eq. (16)). We identify each simulated kilonova as detectable based on the detection thresholds of some typical optical instruments used in the search for kilonovae. We show results for small instruments such as ATLAS (maximum reach of Mmax = 19.5 mag), the medium-size Canada-France-Hawaii Telescope (Mmax = 24.1 mag) and the large Subaru telescope (Mmax = 26.0 mag).

Table 9.

Fraction of BH-NS mergers with detectable kilonova among LIGO/Virgo detectable (S/R > 8) events for assumed exposure times.

We find that the fraction of potentially detectable kilonovae associated with BH-NS mergers is negligible for small instruments: in 10 000 LIGO/Virgo detected BH-NS mergers, only ≲1 could be accompanied by a kilonova. For medium size instruments, the fraction of BH-NS mergers with associated detectable kilonovae varies from small (≲4 kilonova detections per 100 LIGO/Virgo detections) to negligible (2 in 1000). For large telescopes, a slightly larger fraction of kilonovae can be associated with LIGO/Virgo candidates: up to ∼7 kilonovae per 100 LIGO/Virgo BH-NS detections. Much smaller fractions are possible; for example, this fraction drops to 4 kilonovae per 1000 LIGO/Virgo BH-NS detections for our reference model (M280) for an 8-m class telescope.

7. BH-NS merger rates

In Table 10 we list the BH-NS merger rate density (Rd) obtained from our models at redshifts z = 0,  0.5, and 1. Merger rates increase with redshift (at least until redshift z = 2) correlated with an increasing star formation rate. Merger rates for models with high NS/BH natal kicks (M233, M283, M383, and M483) are lower by more than 1 order of magnitude than for models with natal kicks moderated by fallback (M230 and M280).

Table 10.

Intrinsic BH-NS merger rate density Rd (Gpc−3 yr−1).

Comparing between models with fixed physics but different supernova engines, our merger rates are a factor of ∼2 larger in model M230 (rapid supernova engine) compared to model M280 (delayed supernova engine). Both models employ fallback decreased natal kicks. As shown in Fig. 12, BHs begin forming just above MZAMS ∼ 20 M for the rapid supernova engine. In particular, BHs that form from the lowest mass stars (MZAMS ∼ 20 − 23 M) have large enough masses (MBH ∼ 20 M) to form through direct collapse, and therefore these BHs receive no natal kicks. Binaries (potential BH-NS progenitors) with such stars always survive BH formation. On the other hand, for the delayed supernova engine, BHs begin forming just above MZAMS ∼ 15 M. In the mass range MZAMS ∼ 15 − 23 M, these BHs form with low mass (MBH ≲ 10 M) and they tend to receive high natal kicks that can easily disrupt binary progenitors of BH-NS mergers. This leads to the BH-NS merger rate difference between models that employ rapid and delayed supernova engines. There are also some differences between the models for higher initial star masses in terms of BH masses and natal kicks, but they are not that important due to the steep IMF adopted in our calculations. The emergence of a peak in BH mass in the initial-final mass relation associated with the transition between NS and BH formation (MZAMS ∼ 20 M) for the rapid supernova engine, and the continuous gradual BH mass increase with the initial star mass (at least for MZAMS ≲ 35 M) for the delayed supernova engine are the major reasons behind the rate difference. These particular features of BH mass dependence on initial star mass were explained in the context of underlying hydrodynamical simulations of supernova explosions that were used to create both supernova engine models (Fryer et al. 2012; Belczynski et al. 2012b).

thumbnail Fig. 12.

Top: initial (ZAMS) star mass vs. final (Mremnant) mass relation for single stars with a metallicity of Z = 0.001 for two supernova engines: rapid employed (e.g., in model M230) and delayed employed (e.g., in model M280). The blue strip at the bottom marks the NS mass range adopted in our models (MNS = 1 − 2.5 M). Bottom: fallback decreased natal kicks received by remnants (NS or BH) formed out of single stars (Z = 0.001) for a given initial (ZAMS) mass for the two supernova engines. The shape of the adopted initial mass function (IMF) is also shown.

In previous sections we presented all of our results in terms of fractions of the BH-NS merger population. The rates listed in Table 10 can be used to convert these fractions to merger rate densities. In some cases, a high merger rate can compensate for a low fraction of BH-NS mergers with desired properties, or viceversa. For example, we estimated that the highest fraction of LIGO/Virgo detected BH-NS mergers will be associated with kilonovae for models M233 and M383, while other models provide ∼3 − 16 times smaller fractions (see Table 9). However, models M233 and M383 have much smaller (by ≲10) BH-NS merger rate densities than models M230 and M280.

The merger rate density of BH-NS systems was estimated from the O1 and O2 LIGO/Virgo observations to be < 610 Gpc−3 yr−1 (Abbott et al. 2019c). Using one BH-NS merger candidate from the O3 LIGO/Virgo observations (GW190814), Belczynski et al. (2020) estimated the rate density of BH-NS mergers to be in the range of 1.6 − 60 Gpc−3 yr−1. However, it was noted by the LIGO/Virgo Collaboration that this may be a BH-BH merger (Abbott et al. 2020). More recently, based on two BH-NS merger candidates (GW200105 and GW200115), the LIGO/Virgo Collaboration estimated the BH-NS merger rate to be in the range of 7.4 − 320 Gpc−3 yr−1 (90% credible level; The LIGO Scientific Collaboration 2021).

Our calculated BH-NS merger rate densities among the tested models are: 0.2 − 6.7 Gpc−3 yr−1 (z = 0), 0.4 − 13.2 Gpc−3 yr−1 (z = 0.5), and 0.7 − 24.9 Gpc−3 yr−1 (z = 1). These predicted rates are somewhat small compared to the most recent LIGO/Virgo estimates, especially since LIGO/Virgo can currently detect BH-NS mergers only to a redshift of z ∼ 0.1. There are a number of factors that have a significant impact on merger rates and that we have not tested here. The star formation rate and associated cosmic metallicity evolution have nontrivial impacts on derived merger rates (Chruslinska et al. 2019; Chruslinska & Nelemans 2019; Tang et al. 2020; Boco et al. 2021; Santoliquido et al. 2021). New developments to understand the physics of the formation of merging double compact objects may also lead to significant changes in merger rates. For example, in this study we limited our modeling to the formation of BH-NS mergers through the CE isolated binary evolution channel. However, a non-CE channel of isolated binary evolution may produce BH-NS mergers as well, altering the merger rates and BH-NS properties (van den Heuvel et al. 2017; Stevenson et al. 2019; Neijssel et al. 2019; Shao & Li 2021; Olejak et al. 2021). We plan to dedicate a separate study to assessing the impact of these processes on BH-NS merger populations.

8. Conclusions

We performed a suite of binary evolution calculations to model a population of BH-NS mergers. Such mergers are gravitational wave sources and can potentially also produce EM emission (e.g., kilonova, short GRB). They are also interesting in the context of existence (or inexistence) of the FMG, the dearth of compact objects (whether they are NSs or BHs) in the mass range ∼2 − 5 M. At present, it seems that EM observations are beginning to place some compact objects in this mass range (Cromartie et al. 2020; Thompson et al. 2019), and the same is found for GW LIGO/Virgo detections (Abbott et al. 2020). Our study may be summarized as follows:

  1. We predict that only a very small fraction of BH-NS systems will be accompanied by mass ejection during the merger process. The fraction of such systems within the overall BH-NS merger population depends sensitively (0.6 − 9%; see Table 6) on uncertain input physics. For our best guesses of input physics (model M280), we predict that ≲1% of BH-NS mergers are accompanied by mass ejection and by kilonovae bright enough to be detectable with a large 8-m class telescope. This result was obtained for a specific MPA1 EoS (see Fig. 1) and for rather optimistic assumptions on kilonova detectability (see Sect. 2.4). Since the mass ejection (and EM detectability) depends sensitively on BH spin, the future detections (or lack thereof) of BH-NS mergers in EM may distinguish among several existing models of angular momentum transport in massive stars that sets the BH natal spin (Bavera et al. 2020; Belczynski et al. 2020).

  2. If more compact objects within the FMG are found, this would place limits such as those assumed in the “delayed” explosion scenario from Fryer et al. (2012), with weaker explosions with material accretion (e.g., fallback) after the launch of the explosion. A qualitatively similar result is also obtained by Zevin et al. (2020). In our scenarios with a delayed supernova engine, we find that a significant fraction of BH-NS mergers may host a FMG compact object: ≳30% (models M280, M283, M383, and M483; see Table 3).

  3. Our models indicate that extreme BH-NS mass ratio systems may form in classical isolated binary evolution. For example, it is possible for systems to form with a mass ratio (NS mass to BH mass) as low as q = 0.02 with ∼1 M NS and ∼50 M BH. However, such systems are only a very small part of the overall BH-NS merger population, and typical BH-NS mergers are expected to form with mass ratios of q = 0.1 − 0.2.

Our current binary evolution models with BH-NS mergers (and also with BH-BH and NS-NS) mergers are publicly available for further exploration online4.


3

By reverse cumulative distribution we mean the fraction of systems with a higher value than the current argument

Acknowledgments

We would like to thank the anonymous referee for their very useful and detailed comments and suggestions that helped to improve our manuscript. We would like to thank Ben Farr, Martyna Chruslinska, Coen Neijssel, Zoheyr Doctor and all the participants of the “Astrophysics of LIGO/Virgo sources in O3 era” conference who suggested research improvements. R.O.S. would like to thank Erika Holmbeck for her helpful discussions. K.B. and P.D. acknowledge support from the Polish National Science Center (NCN) grant Maestro (2018/30/A/ST9/00050). T.B. was supported by the FNP grant Team/2016-3/19 and NCN grant UMO-2017/26/M/ST9/00978. R.O.S. acknowledges support from NSF AST-1909534.

References

  1. Abbott, B. P., Abbott, R., Abbott, T., et al. 2017, ApJ, 850, L39 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abbott, B., Abbott, R., Abbott, T., et al. 2018a, Phys. Rev. Lett., 121, 161101 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abbott, B. P., Abbott, R., Abbott, T., et al. 2018b, Liv. Rev. Rel., 21, 3 [NASA ADS] [CrossRef] [Google Scholar]
  4. Abbott, B., Abbott, R., Abbott, T., et al. 2019a, Phys. Rev. X, 9, 031040 [NASA ADS] [Google Scholar]
  5. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2019b, ApJ, 882, L24 [Google Scholar]
  6. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2019c, Phys. Rev. X, 9, 031040 [Google Scholar]
  7. Abbott, R., Abbott, T. D., Abraham, S., et al. 2020, ApJ, 896, L44 [Google Scholar]
  8. Abbott, R., Abbott, T. D., & Abraham, S. 2021, ApJ, 915, L5 [NASA ADS] [CrossRef] [Google Scholar]
  9. Ade, P. A., Aghanim, N., Arnaud, M., et al. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Akmal, A., Pandharipande, V., & Ravenhall, D. 1998, Phys. Rev. C, 58, 1804 [Google Scholar]
  11. Alvarez-Castillo, D., Ayriyan, A., Benic, S., et al. 2016, Eur. Phys. J. A, 52, 69 [NASA ADS] [CrossRef] [Google Scholar]
  12. Bailyn, C. D., Jain, R. K., Coppi, P., & Orosz, J. A. 1998, ApJ, 499, 367 [NASA ADS] [CrossRef] [Google Scholar]
  13. Barbieri, C., Salafia, O., Colpi, M., Ghirlanda, G., & Perego, A. 2021, A&A, 654, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Bavera, S. S., Fragos, T., Qin, Y., et al. 2020, A&A, 635, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Belczynski, K., Kalogera, V., & Bulik, T. 2002, ApJ, 572, 407 [NASA ADS] [CrossRef] [Google Scholar]
  16. Belczynski, K., Taam, R. E., Kalogera, V., Rasio, F. A., & Bulik, T. 2007, ApJ, 662, 504 [Google Scholar]
  17. Belczynski, K., Kalogera, V., Rasio, F. A., et al. 2008a, ApJS, 174, 223 [Google Scholar]
  18. Belczynski, K., Taam, R. E., Rantsiou, E., & Van Der Sluys, M. 2008b, ApJ, 682, 474 [NASA ADS] [CrossRef] [Google Scholar]
  19. Belczynski, K., Bulik, T., Fryer, C. L., et al. 2010, ApJ, 714, 1217 [NASA ADS] [CrossRef] [Google Scholar]
  20. Belczynski, K., Bulik, T., & Bailyn, C. 2011, ApJ, 742, L2 [NASA ADS] [CrossRef] [Google Scholar]
  21. Belczynski, K., Bulik, T., & Fryer, C. L. 2012a, ArXiv e-prints [arXiv:1208.2422] [Google Scholar]
  22. Belczynski, K., Wiktorowicz, G., Fryer, C. L., Holz, D. E., & Kalogera, V. 2012b, ApJ, 757, 91 [NASA ADS] [CrossRef] [Google Scholar]
  23. Belczynski, K., Bulik, T., Mandel, I., et al. 2013, ApJ, 764, 96 [NASA ADS] [CrossRef] [Google Scholar]
  24. Belczynski, K., Holz, D. E., Bulik, T., & O’Shaughnessy, R. 2016a, Nature, 534, 512 [NASA ADS] [CrossRef] [Google Scholar]
  25. Belczynski, K., Repetto, S., Holz, D. E., et al. 2016b, ApJ, 819, 108 [NASA ADS] [CrossRef] [Google Scholar]
  26. Belczynski, K., Klencki, J., Fields, C. E., et al. 2020, A&A, 636, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Beldycki, B., & Belczynski, K. 2016, Acta Astron., 66, 347 [NASA ADS] [Google Scholar]
  28. Bhattacharya, M., Kumar, P., & Smoot, G. 2019, MNRAS, 486, 5289 [NASA ADS] [CrossRef] [Google Scholar]
  29. Boco, L., Lapi, A., Chruslinska, M., et al. 2021, ApJ, 907, 110 [NASA ADS] [CrossRef] [Google Scholar]
  30. Chruslinska, M., & Nelemans, G. 2019, MNRAS, 1974 [Google Scholar]
  31. Chruslinska, M., Nelemans, G., & Belczynski, K. 2019, MNRAS, 482, 5012 [Google Scholar]
  32. Côté, B., Belczynski, K., Fryer, C. L., et al. 2017, ApJ, 836, 230 [CrossRef] [Google Scholar]
  33. Cromartie, H., Fonseca, E., Ransom, S. M., et al. 2020, Nat. Astron., 4, 72 [NASA ADS] [CrossRef] [Google Scholar]
  34. de Mink, S. E., & Belczynski, K. 2015, ApJ, 814, 58 [Google Scholar]
  35. Dominik, M., Belczynski, K., Fryer, C., et al. 2012, ApJ, 759, 52 [Google Scholar]
  36. Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [Google Scholar]
  37. Farr, W. M., Sravan, N., Cantrell, A., et al. 2011, ApJ, 741, 103 [NASA ADS] [CrossRef] [Google Scholar]
  38. Fontes, C., Fryer, C., Hungerford, A., Wollaeger, R., & Korobkin, O. 2020, MNRAS, 493, 4143 [CrossRef] [Google Scholar]
  39. Foucart, F., Hinderer, T., & Nissanke, S. 2018, Phys. Rev. D, 98, 081501 [NASA ADS] [CrossRef] [Google Scholar]
  40. Fragione, G., & Loeb, A. 2019, MNRAS, 486, 4443 [NASA ADS] [CrossRef] [Google Scholar]
  41. Fragione, G., Loeb, A., & Rasio, F. A. 2020, ApJ, 895, L15 [NASA ADS] [CrossRef] [Google Scholar]
  42. Fryer, C. L., & Kalogera, V. 2001, ApJ, 554, 548 [Google Scholar]
  43. Fryer, C. L., Belczynski, K., Wiktorowicz, G., et al. 2012, ApJ, 749, 91 [Google Scholar]
  44. Fuller, J., & Ma, L. 2019, ApJ, 881, L1 [Google Scholar]
  45. Fuller, J., Piro, A. L., & Jermyn, A. S. 2019, MNRAS, 485, 3661 [NASA ADS] [Google Scholar]
  46. Gao, H., Ai, S.-K., Cao, Z.-J., et al. 2020, Front. Phys., 15, 1 [Google Scholar]
  47. Gomez, S., Hosseinzadeh, G., Cowperthwaite, P. S., et al. 2019, ApJ, 884, L55 [NASA ADS] [CrossRef] [Google Scholar]
  48. Greiner, J., Burgess, J. M., Savchenko, V., & Yu, H. F. 2016, ApJ, 827, L38 [NASA ADS] [CrossRef] [Google Scholar]
  49. Gupta, A., Gerosa, D., Arun, K., Berti, E., & Sathyaprakash, B. 2020, Phys. Rev. D, 101, 103036 [NASA ADS] [CrossRef] [Google Scholar]
  50. Hinderer, T., Nissanke, S., Foucart, F., et al. 2019, Phys. Rev. D, 100, 063021 [NASA ADS] [CrossRef] [Google Scholar]
  51. Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [Google Scholar]
  52. Hobbs, G., Lorimer, D. R., Lyne, A. G., & Kramer, M. 2005, MNRAS, 360, 974 [Google Scholar]
  53. Hotokezaka, K., & Piran, T. 2017, ApJ, 842, 111 [Google Scholar]
  54. Ivanova, N., Justham, S., Chen, X., et al. 2013, A&ARv, 21, 59 [Google Scholar]
  55. Kawaguchi, K., Kyutoku, K., Shibata, M., & Tanaka, M. 2016, ApJ, 825, 52 [NASA ADS] [CrossRef] [Google Scholar]
  56. Klencki, J., Moe, M., Gladysz, W., et al. 2018, A&A, 619, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Korobkin, O., Wollaeger, R., Fryer, C., et al. 2021, ApJ, 910, 116 [NASA ADS] [CrossRef] [Google Scholar]
  58. Kreidberg, L., Bailyn, C. D., Farr, W. M., & Kalogera, V. 2012, ApJ, 757, 36 [NASA ADS] [CrossRef] [Google Scholar]
  59. Kroupa, P. 2001, MNRAS, 322, 231 [NASA ADS] [CrossRef] [Google Scholar]
  60. Kyutoku, K., Ioka, K., Okawa, H., Shibata, M., & Taniguchi, K. 2015, Phys. Rev. D, 92, 044028 [NASA ADS] [CrossRef] [Google Scholar]
  61. Lee, W. H., & Kluzniak, W. 1995, Acta Astron., 45, 705 [NASA ADS] [Google Scholar]
  62. Li, L.-X., & Paczyński, B. 1998, ApJ, 507, L59 [NASA ADS] [CrossRef] [Google Scholar]
  63. Liu, K., Eatough, R. P., Wex, N., & Kramer, M. 2014, MNRAS, 445, 3115 [CrossRef] [Google Scholar]
  64. MacLeod, M., Antoni, A., Murguia-Berthier, A., Macias, P., & Ramirez-Ruiz, E. 2017, ApJ, 838, 56 [Google Scholar]
  65. Madau, P., & Fragos, T. 2017, ApJ, 840, 39 [Google Scholar]
  66. Metzger, B. D. 2017, Liv. Rev. Rel., 20, 3 [NASA ADS] [CrossRef] [Google Scholar]
  67. Miller, J. M., Ryan, B. R., Dolence, J. C., et al. 2019a, Phys. Rev. D, 100, 023008 [NASA ADS] [CrossRef] [Google Scholar]
  68. Miller, M. C., Lamb, F. K., Dittmann, A. J., et al. 2019b, ApJ, 887, L24 [Google Scholar]
  69. Miyazaki, S., Komiyama, Y., Sekiguchi, M., et al. 2002, PASJ, 54, 833 [Google Scholar]
  70. Mondal, S., Belczyński, K., Wiktorowicz, G., Lasota, J.-P., & King, A. R. 2020, MNRAS, 491, 2747 [Google Scholar]
  71. Neijssel, C. J., Vigna-Gómez, A., Stevenson, S., et al. 2019, MNRAS, 490, 3740 [Google Scholar]
  72. Nissanke, S., Kasliwal, M., & Georgieva, A. 2013, ApJ, 767, 124 [Google Scholar]
  73. Olejak, A., Belczynski, K., & Ivanova, N. 2021, A&A, 651, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Özel, F., Psaltis, D., Narayan, R., & McClintock, J. E. 2010, ApJ, 725, 1918 [CrossRef] [Google Scholar]
  75. Paczyński, B. 1991, Am. Inst. Phys. Conf. Ser., 265, 144 [Google Scholar]
  76. Pavlovskii, K., Ivanova, N., Belczynski, K., & Van, K. X. 2017, MNRAS, 465, 2092 [Google Scholar]
  77. Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Podsiadlowski, P., Joss, P., & Hsu, J. 1992, ApJ, 391, 246 [NASA ADS] [CrossRef] [Google Scholar]
  79. Postnov, K., Kuranov, A., & Simkin, I. 2019, Astron. Lett., 45, 728 [NASA ADS] [CrossRef] [Google Scholar]
  80. Qin, Y., Fragos, T., Meynet, G., et al. 2018, A&A, 616, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Rosswog, S. 2005, ApJ, 634, 1202 [NASA ADS] [CrossRef] [Google Scholar]
  82. Ruffert, M., Janka, H. T., Takahashi, K., & Schaefer, G. 1997, A&A, 319, 122 [Google Scholar]
  83. Sana, H., De Mink, S., de Koter, A., et al. 2012, Science, 337, 444 [NASA ADS] [CrossRef] [Google Scholar]
  84. Santoliquido, F., Mapelli, M., Giacobbo, N., Bouffanais, Y., & Artale, M. C. 2021, MNRAS, 502, 4877 [NASA ADS] [CrossRef] [Google Scholar]
  85. Sedda, M. A. 2020, Commun. Phys., 3, 1 [NASA ADS] [CrossRef] [Google Scholar]
  86. Shao, Y., & Li, X.-D. 2021, ApJ, 920, 81 [NASA ADS] [CrossRef] [Google Scholar]
  87. Spruit, H. C. 1999, A&A, 349, 189 [NASA ADS] [Google Scholar]
  88. Stevenson, S., Sampson, M., Powell, J., et al. 2019, ApJ, 882, 121 [Google Scholar]
  89. Tang, P. N., Eldridge, J. J., Stanway, E. R., & Bray, J. C. 2020, MNRAS, 493, L6 [Google Scholar]
  90. The LIGO Scientific Collaboration, The Virgo Collaboration,& The KAGRA Scientific Collaboration, 2021. ArXiv e-prints [arXiv:2111.03634] [Google Scholar]
  91. Thompson, T. A., Kochanek, C. S., Stanek, K. Z., et al. 2019, Science, 366, 637 [NASA ADS] [CrossRef] [Google Scholar]
  92. van den Heuvel, E. P. J., Portegies Zwart, S. F., & de Mink, S. E. 2017, MNRAS, 471, 4256 [Google Scholar]
  93. Vink, J. S., de Koter, A., & Lamers, H. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Webbink, R. 1984, ApJ, 277, 355 [NASA ADS] [CrossRef] [Google Scholar]
  95. Wiringa, R. B., Fiks, V., & Fabrocini, A. 1988, Phys. Rev. C, 38, 1010 [Google Scholar]
  96. Wollaeger, R. T., Korobkin, O., Fontes, C. J., et al. 2018, MNRAS, 478, 3298 [NASA ADS] [CrossRef] [Google Scholar]
  97. Wollaeger, R. T., Fryer, C. L., Fontes, C. J., et al. 2019, ApJ, 880, 22 [NASA ADS] [CrossRef] [Google Scholar]
  98. Zappa, F., Bernuzzi, S., Pannarale, F., Mapelli, M., & Giacobbo, N. 2019, Phys. Rev. Lett., 123, 041102 [NASA ADS] [CrossRef] [Google Scholar]
  99. Zevin, M., Spera, M., Berry, C. P. L., & Kalogera, V. 2020, ApJ, 899, L1 [Google Scholar]
  100. Zhang, N.-B., & Li, B.-A. 2019, ApJ, 879, 99 [NASA ADS] [CrossRef] [Google Scholar]
  101. Zhu, Y., Wollaeger, R. T., Vassh, N., et al. 2018, ApJ, 863, L23 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Merger rates

In Table A.1 we list the intrinsic merger rate density of BH-NS systems with compact objects in the FMG for all models and for the local Universe (z ∼ 0). Data in this table are the combination of Table 3 and the first column of Table 10. In other words, it is assumed that fractions of BH-NS mergers with various combinations of compact objects in the FMG (Table 3) are constant in the redshift range z = 0 − 1.

Table A.1.

Intrinsic merger rate density [Gpc−3 yr−1] of BH-NS systems in which both (FMGboth) components, only a BH (FMGBH) component, only an NS (FMGNS) component, or no (FMGnone) components are within the FMG.

In Table A.2 we list the merger rate density of LIGO/Virgo detectable BH-NS systems with extreme and small mass ratios for all the models. Data in this Table are the combination of Table 5 and the first column of Table 10 (local Universe z ∼ 0).

Table A.2.

Merger rate density [Gpc−3 yr−1] of LIGO/Virgo detectable (S/R > 8) BH-NS systems with extreme and small mass ratios.

In Table A.3 we list the merger rate density of LIGO/Virgo detectable BH-NS systems with detectable kilonovae for all the models. Data in this Table are the combination of Table 9 and the first column of Table 10 (local Universe z ∼ 0).

Table A.3.

Merger rate density [Gpc−3 yr−1] of BH-NS systems with detectable kilonova among LIGO/Virgo detectable (S/R > 8) events for our assumed exposure times.

All Tables

Table 1.

LIGO/Virgo BH-NS merger candidates.

Table 2.

Calculated models, Supernova engine: R/D–rapid/delayed.

Table 3.

Intrinsic fraction of (z < 1) BH-NS mergers η1 in which both (FMGboth) components, only a BH (FMGBH) component, only an NS (FMGNS) component, or no (FMGnone) components are within the FMG.

Table 4.

Fraction of LIGO/Virgo (S/R > 8) BH-NS systems η3 where both (FMGboth) components, only a BH (FMGBH) component, only an NS (FMGNS) component, or no (FMGnone) components are within the FMG.

Table 5.

Fraction of LIGO/Virgo detectable (S/R > 8) BH-NS mergers with extremely small ratios and extremely small mass ratios.

Table 6.

Intrinsic fraction of BH-NS mergers (z < 1) η2 with any mass ejecta (Mej > 0.001 M) and with significant mass ejecta (Mej > 0.01 M).

Table 7.

Intrinsic BH-NS merger population with any mass ejecta (Mej > 0.001 M, z < 1).

Table 8.

LIGO/Virgo detectable BH-NS mergers (S/R > 8) with any mass ejecta Mej > 0.001 M.

Table 9.

Fraction of BH-NS mergers with detectable kilonova among LIGO/Virgo detectable (S/R > 8) events for assumed exposure times.

Table 10.

Intrinsic BH-NS merger rate density Rd (Gpc−3 yr−1).

Table A.1.

Intrinsic merger rate density [Gpc−3 yr−1] of BH-NS systems in which both (FMGboth) components, only a BH (FMGBH) component, only an NS (FMGNS) component, or no (FMGnone) components are within the FMG.

Table A.2.

Merger rate density [Gpc−3 yr−1] of LIGO/Virgo detectable (S/R > 8) BH-NS systems with extreme and small mass ratios.

Table A.3.

Merger rate density [Gpc−3 yr−1] of BH-NS systems with detectable kilonova among LIGO/Virgo detectable (S/R > 8) events for our assumed exposure times.

All Figures

thumbnail Fig. 1.

Comparison of the NS mass-radius relation of four different equations of state (EoS s) with the LIGO/Virgo constraints from NS-NS merger GW170817 (90% credible limits shown by green and blue boxes: Abbott et al. 2018a). We employ the MPA1 EoS, which reaches ∼2.5 M (our adopted maximum NS mass) and agrees with the LIGO/Virgo constraints. Other relations are WFF1 (Wiringa et al. 1988), APR4 (Akmal et al. 1998), and DD2 (Alvarez-Castillo et al. 2016).

In the text
thumbnail Fig. 2.

Mass distribution of compact objects in BH-NS mergers in the low-redshift Universe (z < 1) for models that employ a rapid supernova engine (M230, M233). There is a very clear mass gap between NSs and BHs. It should be noted that our data become noisy for masses above ≳20 M. However, apparent decline in merger rate density for these higher masses is real and expected in our models.

In the text
thumbnail Fig. 3.

Mass distribution of compact objects in BH-NS mergers in the low-redshift Universe (z < 1) for models that employ a delayed supernova engine (M280, M283). There is no mass gap between NSs and BHs. It should be noted that our data become noisy for masses above ≳20 M, although apparent decline in merger rate density for these higher masses is real and expected in our models.

In the text
thumbnail Fig. 4.

Mass ratio distribution for BH-NS mergers in the low-redshift Universe (z < 1). For each BH-NS system, we calculated the mass ratio q = MNS/MBH ∈ [0, 1]. Values indicate merger rate densities estimated from the systems within each q bin, that is to say dRate dt−1dV−1dq−1Δq, with Δq = 0.01. For transparency purposes, we use simple merger rate density as the vertical axis label. Most mergers have rather small mass ratios q ∼ 0.1 − 0.2, but there are also extreme mass ratio systems such as 56 M BH + 1.3 M NS (q = 0.023). For models with the delayed supernova engine (M280 and M283), there is a secondary peak for high mass ratios q > 0.4, which consists mostly of double compact objects with both components in the FMG (see Sect. 4).

In the text
thumbnail Fig. 5.

Cumulative distribution of LIGO/Virgo detectable (S/R > 8) BH-NS mergers in terms of mass ratio (see also Table 5).

In the text
thumbnail Fig. 6.

Evolutionary route leading to the formation of the 38.9 M + 1.3 M (q = 0.033) BH-NS merger for metallicity Z = 0.0001 (model M280). The evolution from ZAMS to the formation of the BH-NS system takes ∼30.4 Myr, and then it takes ∼2.29 Gyr for the system to in-spiral (angular momentum loss due to the emission of gravitational waves), leading to the final merger of the two compact objects. (CHeB: core helium burning star, He star: naked helium star, HG: Hertzsprung gap star). No tidal spin-up is expected due to large orbital separations when binary components are naked helium stars: a ∼ 25 R (Porb = 2.1d; star A) and a ∼ 200 R (Porb = 60d; star B; see Sect. 2.1).

In the text
thumbnail Fig. 7.

Evolution of two massive stars that leads to the formation of the moderate mass ratio (q = 0.12) BH-NS (14.7 + 1.8 M) merger (model M230). The top panel shows the mass evolution of the primary and secondary stars, while the bottom panel shows the changes in the semimajor axis and eccentricity. Evolution from ZAMS to the BH-NS system formation takes 10.5 Myr, and then it takes another 2.7 Gyr for the two compact objects to merge and produce gravitational wave radiation that is potentially detectable by LIGO/Virgo. No mass ejection (no kilonova or short GRB) is expected in this case.

In the text
thumbnail Fig. 8.

Evolution of binary system that forms a 3.5 + 2.1  M BH-NS merger (mass ratio of q = 0.6). We note that both compact objects are within the FMG and we expect that this merger will not lead to any mass ejection due to a high NS mass.

In the text
thumbnail Fig. 9.

Intrinsic distribution of ejecta mass in BH-NS mergers with z < 1. The structures at Mej ≳ 0.03 M are shot noise.

In the text
thumbnail Fig. 10.

Intrinsic reverse cumulative distribution of ejecta mass (> Mej) in BH-NS mergers with z < 1.

In the text
thumbnail Fig. 11.

Kilonova brightness distribution for LIGO/Virgo detectable (S/R > 8) BH-NS mergers.

In the text
thumbnail Fig. 12.

Top: initial (ZAMS) star mass vs. final (Mremnant) mass relation for single stars with a metallicity of Z = 0.001 for two supernova engines: rapid employed (e.g., in model M230) and delayed employed (e.g., in model M280). The blue strip at the bottom marks the NS mass range adopted in our models (MNS = 1 − 2.5 M). Bottom: fallback decreased natal kicks received by remnants (NS or BH) formed out of single stars (Z = 0.001) for a given initial (ZAMS) mass for the two supernova engines. The shape of the adopted initial mass function (IMF) is also shown.

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.