Evidence for a Correlation between Binary Black Hole Mass Ratio and Black Hole Spins

The astrophysical origins of the binary black hole systems seen with gravitational waves are still not well understood. However, features in the distribution of black hole masses, spins, redshifts, and eccentricities provide clues into how these systems form. Much has been learned by investigating these distributions one parameter at a time. However, we can extract additional information by studying the covariance between pairs of parameters. Previous work has shown preliminary support for an anticorrelation between mass ratio q ≡ m ₂/m ₁ and effective inspiral spin χ _eff in the binary black hole population. In this study, we test for the existence of this anticorrelation using updated data from the third gravitational-wave transient catalog and improve our copula-based framework to employ a more robust model for black hole spins. We find evidence for an anticorrelation in (q, χ _eff) with 99.7% credibility. This may imply high common-envelope efficiencies, stages of super-Eddington accretion, or a tendency for binary black hole systems to undergo mass-ratio reversal during isolated evolution. Covariance in (q, χ _eff) may also be used to investigate the physics of tidal spinup as well as the properties of binary black hole-forming active galactic nuclei.