Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3

Abstract

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mrow><a:mn>1700</a:mn><a:mtext> </a:mtext><a:mtext> </a:mtext><a:msup><a:mrow><a:mi>Gpc</a:mi></a:mrow><a:mrow><a:mo>−</a:mo><a:mn>3</a:mn></a:mrow></a:msup><a:mtext> </a:mtext><a:msup><a:mrow><a:mi>yr</a:mi></a:mrow><a:mrow><a:mo>−</a:mo><a:mn>1</a:mn></a:mrow></a:msup></a:mrow></a:math> and the neutron star–black hole merger rate to be between 7.8 and <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:mrow><c:mn>140</c:mn><c:mtext> </c:mtext><c:mtext> </c:mtext><c:msup><c:mrow><c:mi>Gpc</c:mi></c:mrow><c:mrow><c:mo>−</c:mo><c:mn>3</c:mn></c:mrow></c:msup><c:mtext> </c:mtext><c:msup><c:mrow><c:mi>yr</c:mi></c:mrow><c:mrow><c:mo>−</c:mo><c:mn>1</c:mn></c:mrow></c:msup></c:mrow></c:math>, assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:mrow><e:mn>44</e:mn><e:mtext> </e:mtext><e:mtext> </e:mtext><e:msup><e:mrow><e:mi>Gpc</e:mi></e:mrow><e:mrow><e:mo>−</e:mo><e:mn>3</e:mn></e:mrow></e:msup><e:mtext> </e:mtext><e:msup><e:mrow><e:mi>yr</e:mi></e:mrow><e:mrow><e:mo>−</e:mo><e:mn>1</e:mn></e:mrow></e:msup></e:mrow></e:math> at a fiducial redshift (<g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"><g:mi>z</g:mi><g:mo>=</g:mo><g:mn>0.2</g:mn></g:math>). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:mo stretchy="false">(</i:mo><i:mn>1</i:mn><i:mo>+</i:mo><i:mi>z</i:mi><i:msup><i:mo stretchy="false">)</i:mo><i:mi>κ</i:mi></i:msup></i:math> with <m:math xmlns:m="http://www.w3.org/1998/Math/MathML" display="inline"><m:mi>κ</m:mi><m:mo>=</m:mo><m:mn>2.</m:mn><m:msubsup><m:mn>9</m:mn><m:mrow><m:mo>−</m:mo><m:mn>1.8</m:mn></m:mrow><m:mrow><m:mo>+</m:mo><m:mn>1.7</m:mn></m:mrow></m:msubsup></m:math> for <o:math xmlns:o="http://www.w3.org/1998/Math/MathML" display="inline"><o:mi>z</o:mi><o:mo>≲</o:mo><o:mn>1</o:mn></o:math>. Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from <q:math xmlns:q="http://www.w3.org/1998/Math/MathML" display="inline"><q:msubsup><q:mn>1.2</q:mn><q:mrow><q:mo>−</q:mo><q:mn>0.2</q:mn></q:mrow><q:mrow><q:mo>+</q:mo><q:mn>0.1</q:mn></q:mrow></q:msubsup></q:math> to <s:math xmlns:s="http://www.w3.org/1998/Math/MathML" display="inline"><s:msubsup><s:mn>2.0</s:mn><s:mrow><s:mo>−</s:mo><s:mn>0.3</s:mn></s:mrow><s:mrow><s:mo>+</s:mo><s:mn>0.3</s:mn></s:mrow></s:msubsup><s:msub><s:mi>M</s:mi><s:mo stretchy="false">⊙</s:mo></s:msub></s:math>. We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of <v:math xmlns:v="http://www.w3.org/1998/Math/MathML" display="inline"><v:msubsup><v:mn>8.3</v:mn><v:mrow><v:mo>−</v:mo><v:mn>0.5</v:mn></v:mrow><v:mrow><v:mo>+</v:mo><v:mn>0.3</v:mn></v:mrow></v:msubsup></v:math> and <x:math xmlns:x="http://www.w3.org/1998/Math/MathML" display="inline"><x:msubsup><x:mn>27.9</x:mn><x:mrow><x:mo>−</x:mo><x:mn>1.8</x:mn></x:mrow><x:mrow><x:mo>+</x:mo><x:mn>1.9</x:mn></x:mrow></x:msubsup><x:msub><x:mi>M</x:mi><x:mo stretchy="false">⊙</x:mo></x:msub></x:math>. While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately <ab:math xmlns:ab="http://www.w3.org/1998/Math/MathML" display="inline"><ab:mn>60</ab:mn><ab:msub><ab:mi>M</ab:mi><ab:mo stretchy="false">⊙</ab:mo></ab:msub></ab:math>, which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below <db:math xmlns:db="http://www.w3.org/1998/Math/MathML" display="inline"><db:msub><db:mi>χ</db:mi><db:mi>i</db:mi></db:msub><db:mo>≈</db:mo><db:mn>0.25</db:mn></db:math>. While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum. Published by the American Physical Society 2023