The Redshift Evolution of the Binary Black Hole Mass Distribution from Dense Star Clusters

Gravitational-wave detectors are unveiling a population of binary black hole (BBH) mergers out to redshifts z approximate to 1, and are starting to constrain how the BBH population evolves with redshift. We present predictions for the redshift evolution of the BBH mass and spin distributions for systems originating from dense star clusters. Utilizing a grid of 144 state-of-the-art dynamical models for globular clusters, we demonstrate that BBH merger rates peak at higher redshifts for larger black hole primary masses M 1. Specifically, for M 1 greater than or similar to 40 M circle dot, the BBH merger rate reaches its peak at redshift z approximate to 2.1, while for M 1 less than or similar to 20 M circle dot, the peak occurs at z approximate to 1.1, assuming that the cluster formation rate peaks at z = 2.2. The average BBH primary mass also increases from similar to 10 M circle dot at z = 0 to similar to 30 M circle dot at z = 10. We show that similar to 20% BBHs contain massive remnants from next-generation mergers, with this fraction increasing (decreasing) for larger (smaller) primary masses. This difference is not large enough to significantly alter the effective spins of the BBH population originating from globular clusters, and we find that their effective spin distribution does not evolve across cosmic time. These findings can be used to distinguish BBHs from dense star clusters by future gravitational-wave observations.