Binary black holes in which both spins are aligned with the binary's orbital angular momentum do not precess. However, the up-down configuration, in which the spin of the heavier (lighter) black hole is aligned (anti-aligned) with the orbital angular momentum, is unstable to spin precession at small orbital separations. We first cast the spin precession problem in terms of a simple harmonic oscillator and provide a cleaner derivation of the instability onset. Surprisingly, we find that following the instability, up-down binaries do not disperse in the available parameter space but evolve toward precise endpoints. We then present an analytic scheme to locate these final configurations and confirm them with numerical integrations. Namely, unstable up-down binaries approach mergers with the two spins coaligned with each other and equally misaligned with the orbital angular momentum. Merging up-down binaries relevant to LIGO/Virgo and LISA may be detected in these endpoint configurations if the instability onset occurs prior to the sensitivity threshold of the detector. As a by-product, we obtain new generic results on binary black hole spin-orbit resonances at 2nd~post-Newtonian order. We finally apply these findings to a simple astrophysical population of binary black holes where a formation mechanism aligns the spins without preference for co- or counteralignment, as might be the case for stellar-mass black holes embedded in the accretion disk of a supermassive black hole.

Multiple approaches are required to study the evolution of black-hole binaries. While the post-Newtonian approximation is sufficient to describe the early inspiral (even from infinitely large orbital separation), only numerical relativity can capture the full complexity of the dynamics near merger. We combine multi-timescale post-Newtonian integrations with numerical-relativity surrogate models, thus mapping the entire history of the binary from its asymptotic configuration at past-time infinity to the post-merger remnant. This approach naturally allows us to assess the impact of the precessional and orbital phase on the properties - mass, spin, and kick - of the merger remnant. These phases introduce a fundamental uncertainty when connecting the two extrema of the binary evolution.

Two of the dominant channels to produce the black-hole binary mergers observed by LIGO and Virgo are believed to be the isolated evolution of stellar binaries in the field and dynamical formation in star clusters. Their relative efficiency can be characterized by a "mixing fraction." Pair instabilities prevent stellar collapse from generating black holes more massive than about $45 M_\odot$. This "mass gap" only applies to the field formation scenario, and it can be filled by repeated mergers in clusters. A similar reasoning applies to the binary's effective spin. If black holes are born slowly rotating, the high-spin portion of the parameter space (the "spin gap") can only be populated by black hole binaries that were assembled dynamically. Using a semianalytical cluster model, we show that future gravitational-wave events in either the mass gap, the spin gap, or both can be leveraged to infer the mixing fraction between the field and cluster formation channels.