Total: 1
Measurements of the binary black hole spin distribution from the growing catalog of gravitational-wave observations can help elucidate the astrophysical processes shaping the formation and evolution of these systems. Spin-orbit resonances are one process of interest, in which the component spin vectors and the orbital angular momentum align into a common plane and jointly precess about the total angular momentum of the system. These resonances, which occur preferentially in systems formed via isolated binary evolution with strong tidal effects, lead to excesses in the distribution of the azimuthal angle between the projections of the component spin vectors onto the orbital plane at ϕ12=0,±π. Previous analyses have demonstrated that this parameter is particularly difficult to constrain for individual binaries. In this work, we conduct the first hierarchical analysis modeling the population-level distribution of ϕ12 simultaneously with the other mass and spin parameters for simulated binary black hole populations to determine whether spin-orbit resonances can be reliably constrained. While we are unlikely to find definitive evidence for spin-orbit resonances with a population of the size expected by the end of the ongoing LIGO-Virgo-KAGRA fourth observing run, we correctly recover the various ϕ12 distributions we simulate within uncertainties. We find that we can place meaningful constraints on the relative excesses at ϕ12=0,±π, which encodes information about mass transfer in the formation of the binary. We can also distinguish between fully isotropic spin angle distributions and those with features in the spin azimuth and tilt distributions. Thus, we show that population-level measurements of the ϕ12 distribution offer a reliable, novel way to probe binary formation channels, dynamics, and mass transfer with gravitational-wave observations.