Current population models of binary black hole distributions are difficult to interpret because standard population inferences hinge on modeling choices, which can mask or mimic real structure. The maximum population likelihood `` formalism'' provides a means to investigate and interpret features in the distribution of binary black holes using only data -- without specifying a population model. It tells us if features inferred from current population models are truly present in the data or if they arise from model misspecification. It also provides guidance for developing new models by highlighting previously unnoticed features. In this study, we utilize the formalism to examine the binary black hole population in the LIGO--Virgo--KAGRA (LVK) fourth Gravitational-Wave Transient Catalog (GWTC-4). Our analysis supports the existence of a gap around 45\, M_ in the secondary black hole mass distribution and identifies a widening in the distribution of the effective inspiral spin parameter χₑff near this mass as recently reported by Tong et al. (2025). Similar to earlier studies, we find support for an anti-correlation between χₑff and mass ratio. However, we argue that this may be a spurious correlation arising from misspecification of the joint distribution of black hole masses. Furthermore, we identify support for dimensionless black hole spin magnitudes at approximately χ 0. 2 and χ0. 7. The data support the existence of a correlation between the spin magnitudes χ₁ and χ₂, though subsequent study is required to determine if this feature is statistically significant. The accompanying data release includes samples, which can be used to compare theoretical predictions to LVK data and to assess assumptions in parameterised models.
Guttman et al. (Thu,) studied this question.