Seeking Spinning Subpopulations of Black Hole Binaries via Iterative Density Estimation

Attempts to understand the formation of binary black hole (BBH) systems detected via gravitational wave (GW) emission are affected by many unknowns and uncertainties, from both the observational and theoretical (astrophysical modelling) sides. Binary component spins have been proposed as a means to investigate formation channels, however obtaining clear inferences is challenging, given the apparently low magnitude of almost all merging BH spins and their high measurement uncertainties. Even for the effective aligned spin $\chi_{\mathrm{eff}}$ which is more precisely measured than component spins, specific model assumptions have been required to identify any clear trends. Here, we reconstruct the joint component mass and $\chi_{\mathrm{eff}}$ distribution of BBH mergers with minimal assumptions using the GWTC-3 catalog, using an iterative kernel density estimation (KDE)-based method. We reproduce some features seen in previous analyses, for instance a small but preferentially positive $\chi_{\mathrm{eff}}$ for low-mass mergers; we also identify a possible subpopulation of higher-spin BBH with $|\chi_{\mathrm{eff}}|$ up to $\sim\! 0.75$ for primary masses $m_1 \gtrsim 40\,M_\odot$, in addition to the bulk of the distribution with $|\chi_{\mathrm{eff}}| \lesssim 0.2$. This finding is consistent with previous studies indicating a broader spin distribution at high mass, suggesting a distinct origin for the high-spin systems. We also identify a new potential trend of low-mass BBHs: the \emph{derivative} of $\chi_{\mathrm{eff}}$ with respect to $m_1$ ($m_2$) is positive (negative) over the $10$--$15\,M_\odot$ range. This apparent structure may be related to} a previously reported anticorrelation between mass ratio and $\chi_{\mathrm{eff}}$.