Hints of spin-magnitude correlations and a rapidly spinning subpopulation of binary black holes

The complex astrophysical processes leading to the formation of binary black holes and their eventual merger are imprinted on the spins of the individual black holes. We revisit the astrophysical distribution of those spins based on gravitational waves from the third gravitational wave transient catalog GWTC-3, (Abbott et al. 2023a), looking for structure in the two-dimensional space defined by the dimensionless spin magnitudes of the heavier ($\chi_1$) and lighter ($\chi_2$) component black holes. We find support for two distinct subpopulations with greater than $95\%$ credibility. The dominant population is made up of black holes with small spins, preferring $\chi_1 \approx 0.2$ for the primary and $\chi_2 \approx 0$ for the secondary; we report signs of an anticorrelation between $\chi_1$ and $\chi_2$, as well as as evidence against a subpopulation of binaries in which both components are nonspinning. The subdominant population consists of systems in which both black holes have relatively high spins and contains $20^{+18}_{-18}\%$ of the binaries. The binaries that are most likely to belong in this subpopulation are massive and slightly more likely to have spin-orientations aligned with the orbital angular momentum--potentially consistent with isolated binary formation channels capable of producing large spins, like chemically homogeneous evolution. This hint of a rapidly spinning subpopulation hinges on GW190517, a binary with large and well-measured spins. Our results, which are enabled by novel hierarchical inference methods, represent a first step towards more descriptive population models for black hole spins, and will be strengthened or refuted by the large number of gravitational wave detections expected in the next several years.