No Arabic abstract
We develop a new method based on Gaussian process to reconstruct the mass distribution of binary black holes (BBHs). Instead of prespecifying the formalisms of mass distribution, we introduce a more flexible and nonparametric model with which the distribution can be mainly determined by the observed data. We first test our method with simulated data, and find that it can well recover the injected distribution. Then we apply this method to analyze the data of BBHs observations from LIGO-Virgo Gravitational-Wave Transient Catalog 2. By reconstructing the chirp mass distribution, we find that there is a peak or a platform locating at $20-30,M_{odot}$ rather than a single-power-law-like decrease from low mass to high mass. Moreover, one or two peaks in the chirp mass range of $mathcal{M}<20,M_{odot}$ may be favored by the data. Assuming a mass-independent mass ratio distribution of $p(q)propto q^{1.4}$, we further obtain a distribution of primary mass, and find that there is a feature locating in the range of $(30, 40),M_{odot}$, which can be related to textsc{Broken Power Law} and textsc{Power Law + Peak} distributions described in The LIGO Scientific Collaboration et al. (2020). Besides, the merger rate of BBHs is estimated to $mathcal{R} = 26.29^{+14.21}_{-8.96}~{rm Gpc^{-3}~yr^{-1}}$ supposing there is no redshift evolution.
With the black hole mass function (BHMF; assuming an exponential cutoff at a mass of $sim 40,M_odot$) of coalescing binary black hole systems constructed with the events detected in the O1 run of the advanced LIGO/Virgo network, Liang et al.(2017) predicted that the birth of the lightest intermediate mass black holes (LIMBHs; with a final mass of $gtrsim 100,M_odot$) is very likely to be caught by the advanced LIGO/Virgo detectors in their O3 run. The O1 and O2 observation run data, however, strongly favor a cutoff of the BHMF much sharper than the exponential one. In this work we show that a power-law function followed by a sudden drop at $sim 40,M_odot$ by a factor of $sim $a few tens and then a new power-law component extending to $geq 100M_odot$ are consistent with the O1 and O2 observation run data. With this new BHMF, quite a few LIMBH events can be detected in the O3 observation run of advanced LIGO/Virgo. The first LIMBH born in GW190521, an event detected in the early stage of the O3 run of advanced LIGO/Virgo network, provides additional motivation for our hypothesis.
We analyze the LIGO/Virgo GWTC-2 catalog to study the primary mass distribution of the merging black holes. We perform hierarchical Bayesian analysis, and examine whether the mass distribution has a sharp cutoff for primary black hole masses below $65 M_odot$, as predicted in pulsational pair instability supernova model. We construct two empirical mass functions. One is a piece-wise function with two power-law segments jointed by a sudden drop. The other consists of a main truncated power-law component, a Gaussian component, and a third very massive component. Both models can reasonably fit the data and a sharp drop of the mass distribution is found at $sim 50M_odot$, suggesting that the majority of the observed black holes can be explained by the stellar evolution scenarios in which the pulsational pair-instability process takes place. On the other hand, the very massive sub-population, which accounts for at most several percents of the total, may be formed through hierarchical mergers or other processes.
Many proposed scenarios for black hole (BH) mergers involve a tertiary companion that induces von Zeipel-Lidov-Kozai (ZLK) eccentricity cycles in the inner binary. An attractive feature of such mechanisms is the enhanced merger probability when the octupole-order effects, also known as the eccentric Kozai mechanism, are important. This can be the case when the tertiary is of comparable mass to the binary components. Since the octupole strength [$propto (1-q)/(1+q)$] increases with decreasing binary mass ratio $q$, such ZLK-induced mergers favor binaries with smaller mass ratios. We use a combination of numerical and analytical approaches to fully characterize the octupole-enhanced binary BH mergers and provide analytical criteria for efficiently calculating the strength of this enhancement. We show that for hierarchical triples with semi-major axis ratio $a/a_{rm out}gtrsim 0.01$-$0.02$, the binary merger fraction can increase by a large factor (up to $sim 20$) as $q$ decreases from unity to $0.2$. The resulting mass ratio distribution for merging binary BHs produced in this scenario is in tension with the observed distribution obtained by the LIGO/VIRGO collaboration, although significant uncertainties remain about the initial distribution of binary BH masses and mass ratios.
On June 8, 2017 at 02:01:16.49 UTC, a gravitational-wave signal from the merger of two stellar-mass black holes was observed by the two Advanced LIGO detectors with a network signal-to-noise ratio of 13. This system is the lightest black hole binary so far observed, with component masses $12^{+7}_{-2},M_odot$ and $7^{+2}_{-2},M_odot$ (90% credible intervals). These lie in the range of measured black hole masses in low-mass X-ray binaries, thus allowing us to compare black holes detected through gravitational waves with electromagnetic observations. The sources luminosity distance is $340^{+140}_{-140}$ Mpc, corresponding to redshift $0.07^{+0.03}_{-0.03}$. We verify that the signal waveform is consistent with the predictions of general relativity.
We study the impact of mass-transfer physics on the observable properties of binary black hole populations formed through isolated binary evolution. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of $chi_{eff}$, $M_{chirp}$ and $q$. We find that low common envelope efficiency translates to tighter orbits post common envelope and therefore more tidally spun up second-born black holes. However, these systems have short merger timescales and are only marginally detectable by current gravitational-waves detectors as they form and merge at high redshifts ($zsim 2$), outside current detector horizons. Assuming Eddington-limited accretion efficiency and that the first-born black hole is formed with a negligible spin, we find that all non-zero $chi_{eff}$ systems in the detectable population can come only from the common envelope channel as the stable mass-transfer channel cannot shrink the orbits enough for efficient tidal spin-up to take place. We find the local rate density ($zsimeq 0.01$) for the common envelope channel is in the range $sim 17-113~Gpc^{-3}yr^{-1}$ considering a range of $alpha_{CE} in [0.2,5.0]$ while for the stable mass transfer channel the rate density is $sim 25~Gpc^{-3}yr^{-1}$. The latter drops by two orders of magnitude if the mass accretion onto the black hole is not Eddington limited because conservative mass transfer does not shrink the orbit as efficiently as non-conservative mass transfer does. Finally, using GWTC-2 events, we constrain the lower bound of branching fraction from other formation channels in the detected population to be $sim 0.2$. Assuming all remaining events to be formed through either stable mass transfer or common envelope channels, we find moderate to strong evidence in favour of models with inefficient common envelopes.