Physics-based phenomenological modeling of binary black hole hierarchical formation 1: Synthetic universes from globular cluster simulations for GWTC

Turning simulations into synthetic universes

Gravitational-wave detectors are building a census of merging black holes. The challenge is to work out what that census says about where the binaries formed. The paper arXiv:2609.05996v1 introduces a way to make that comparison more physical: instead of fitting an abstract distribution alone, it builds synthetic universes from simulations of globular clusters and compares them with the GWTC-5.0 catalog.

The framework is designed to let the underlying physical parameters vary continuously. Cluster mass, metallicity, formation redshift, compactness, and the spins black holes are born with can all change within the population model. That makes it possible to ask a sharper question: which combinations of physical conditions produce a catalog resembling the one we observe?

What happens inside a dense cluster?

In a globular cluster, black holes can meet, merge, and leave behind a heavier remnant that merges again. These hierarchical mergers provide a natural route to the high-mass part of the catalog. In the model, they also produce a wider spread of effective spins, and the abstract reports that this behavior is reproduced when black holes are born with zero natal spin.

The cluster population is compared alongside a phenomenological field, or isolated-binary, channel. This second population supplies lower-mass systems with preferentially aligned spins. Including both channels changes the inferred cluster-compactness picture: ordinary dense globular-cluster birth radii can fit the data without requiring environments as extreme as nuclear clusters.

The model also includes an intermediate-mass, isotropic component to address residual tension between 15 and 30 solar masses. The authors interpret it as field remnants reprocessed in clusters, while stressing that this is a hypothesis—not a detection—and that current numerical support does not allow a reliable evidence comparison.

Predictions that grow more useful with every event

After physical normalization, the model gives a high-mass globular-cluster fraction of about 0.39%. It reports local rates of roughly 9.1 Gpc^-3 yr^-1 for the cluster channel and 16.0 Gpc^-3 yr^-1 for the field channel. It also predicts linked features: mass-spectrum breaks near 35 and 70 solar masses, a mass-ratio feature around q ≈ 0.5 from first-plus-second-generation pairings, and a symmetric effective-spin distribution that broadens sharply above 45 solar masses.

These correlations are the point of the exercise. As the gravitational-wave catalog grows, observations can test not just one preferred mass or spin distribution, but detailed physical models whose parameters move continuously through plausible formation scenarios. That turns a growing list of black-hole mergers into a progressively more discriminating probe of how stellar environments build them.