A Flexible Gaussian Process Reconstruction Method and the Mass Function of the Coalescing Binary Black Hole Systems
Yin-Jie 银杰 Li 李, Yuan-Zhu 远瞩 Wang 王, Ming-Zhe 明哲 Han 韩, Shao-Peng 少鹏 Tang 唐, Qiang 强 Yuan 袁, Yi-Zhong 一中 Fan 范, Da-Ming 大明 Wei 韦
Abstract
Abstract We develop a new method based on Gaussian processes 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 located at 20–30 M ⊙ 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 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi mathvariant="italic"></mml:mi> <mml:mo><</mml:mo> <mml:mn>20</mml:mn> <mml:mspace width="0.25em"/> <mml:msub> <mml:mrow> <mml:mi>M</mml:mi> </mml:mrow> <mml:mrow> <mml:mo>⊙</mml:mo> </mml:mrow> </mml:msub> </mml:math> may be favored by the data. Assuming a mass-independent mass ratio distribution of p ( q ) ∝ q 1.4 , we further obtain a distribution of primary mass and find that there is a feature located in the range of (30, 40) M ⊙ , which can be related to Broken Power Law and Power Law + Peak distributions described in Collaboration et al. Besides, the merger rate of BBHs is estimated to be <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi mathvariant="italic"></mml:mi> <mml:mo>=</mml:mo> <mml:msubsup> <mml:mrow> <mml:mn>26.29</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>8.96</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>14.21</mml:mn> </mml:mrow> </mml:msubsup> <mml:mspace width="0.25em"/> <mml:msup> <mml:mrow> <mml:mi>Gpc</mml:mi> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>3</mml:mn> </mml:mrow> </mml:msup> <mml:mspace width="0.25em"/> <mml:msup> <mml:mrow> <mml:mi>yr</mml:mi> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> </mml:math> , supposing that there is no redshift evolution.