2020ApJ...898..152S

Next generation ground-based gravitational-wave detectors will observe binary black hole (BBH) mergers up to redshift z≳10, probing the evolution of compact binary (CB) mergers across cosmic time. Here, we present a new data-driven model to estimate the cosmic merger rate density (MRD) evolution of CBs, by coupling catalogs of CB mergers with observational constraints on the cosmic star formation rate (SFR) density and on the metallicity evolution of the universe. We adopt catalogs of CB mergers derived from recent N-body and population-synthesis simulations, to describe the MRD of CBs formed in young star clusters (hereafter, dynamical CBs) and in the field (hereafter, isolated CBs). The local MRD of dynamical BBHs is R_BBH=64_–20⁺³⁴ Gpc^–3 yr^–1, consistent with the 90% credible interval from the first and second observing runs (O1 and O2) of the LIGO-Virgo collaboration, and with the local MRD of isolated BBHs ( R_BBH=50_–37⁺⁷¹ Gpc^–3 yr^–1). The local MRD of dynamical and isolated black hole-neutron star binaries is R_BHNS=41_–23⁺³³ and 49_–34⁺⁴⁸ Gpc^–3 yr^–1, respectively. Both values are consistent with the upper limit inferred from O1 and O2. Finally, the local MRD of dynamical binary neutron stars (BNSs, R_BNS=151_–38⁺⁵⁹ Gpc^–3 yr^–1) is a factor of two lower than the local MRD of isolated BNSs ( R_BNS=283_–75⁺⁹⁷Gpc^–3 yr^–1). The MRD for all CB classes grows with redshift, reaching its maximum at z∈[1.5,2.5], and then decreases. This trend springs from the interplay between cosmic SFR, metallicity evolution, and delay time of binary compact objects.