Gravitational-wave Observations Suggest Most Black Hole Mergers Form in Triples

The spin-orbit tilt angles $θ_{1(2)}$ of merging stellar-mass black holes provide key insights into their astrophysical origin. Non-parametric population modelling of The LIGO, Virgo, and KAGRA Collaborations (2025a, arXiv:2508.18083) shows that the spin-orbit tilt distribution of mergers in the latest Gravitational-Wave Transient Catalog 4.0 exhibits a global peak at near-perpendicular directions $\cosθ_{1(2)}\approx0$. Here, we recover this feature using hierarchical Bayesian inference with parametric models that are tailored to enhance the diagnostic power about astrophysical formation channels. We find that the spin distribution of the low-mass bulk of the binary black hole merger population ($m_1\lesssim 44.3^{+8.7}_{-4.6}\,\rm M_\odot$) can be well-modelled by a dominant Gaussian component that peaks at $\cosθ_{1(2)}\approx0$, possibly mixed with a subdominant isotropic component. Models that include a component with spins preferentially aligned with the orbit are disfavoured by current data (with Bayes factors $|Δ\ln\mathcal{B}|\approx1$ to $3$) and constrain its contribution to be likely small ($ξ\sim\mathcal{O}(1)\,\%$), although large contributions cannot yet ruled out with certainty. If these findings are reinforced by more detections, they would challenge any major contribution from the traditional isolated-binary formation scenario yielding closely aligned spins. Instead, the dominant component with near-perpendicular spins matches expectations from the evolution of isolated massive stellar triples in the galactic field, where the Lidov-Kozai effect naturally produces a unique overabundance of mergers with $\cosθ_{1(2)}\approx0$.