The second LIGO-Virgo catalog of gravitational-wave (GW) transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes that impact the observational properties of mergers. Given our set of formation models, we find significant differences between the branching fractions of the underlying and detectable populations, and that the diversity of detections suggests that multiple formation channels are at play. A mixture of channels is strongly preferred over any single channel dominating the detected population: an individual channel does not contribute to more than 70% of the observational sample of binary black holes. We calculate the preference between the natal spin assumptions and common-envelope efficiencies in our models, favoring natal spins of isolated black holes of 0.1, and marginally preferring common-envelope efficiencies of 2.0 while strongly disfavoring highly inefficient common envelopes. We show that it is essential to consider multiple channels when interpreting GW catalogs, as inference on branching fractions and physical prescriptions becomes biased when contributing formation scenarios are not considered or incorrect physical prescriptions are assumed. Although our quantitative results can be affected by uncertain assumptions in model predictions, our methodology is capable of including models with updated theoretical considerations and additional formation channels.

Context. After years of scientific progress, the origin of stellar binary black holes is still a great mystery. Several formation channels for merging black holes have been proposed in the literature. As more merger detections are expected with future gravitational-wave observations, population synthesis studies can help to distinguish between them. Aims. We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of observables such as black-hole spins, masses and cosmological redshifts of mergers. Methods. To achieve this aim, we used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. Results. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the following three main gravitational-wave observables: the effective inspiral spin parameter χeff, the chirp mass Mchirp and the cosmological redshift of merger zmerger. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect approximately 80% of events with χeff <  0.1, while the remaining 20% of events with χeff ≥ 0.1 are split into ∼10% with Mchirp <  15 M⊙ and ∼10% with Mchirp ≥ 15 M⊙. Moreover, we find that Mchirp and χeff distributions are very weakly dependent on the detector sensitivity. Conclusions. The favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase.

We study the impact of mass-transfer physics on the observable properties of binary black hole populations that formed through isolated binary evolution. We used the POSYDON framework to combine detailed MESA binary simulations with the COSMIC population synthesis tool to obtain an accurate estimate of merging binary black hole observables with a specific focus on the spins of the black holes. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of the effective inspiral spin parameter χeff, chirp mass Mchirp, and binary mass ratio q. We find that low common envelope efficiency translates to tighter orbits following the common envelope and therefore more tidally spun up second-born black holes. However, these systems have short merger timescales and are only marginally detectable by current gravitational-wave detectors as they form and merge at high redshifts (z ∼ 2), outside current detector horizons. Assuming Eddington-limited accretion efficiency and that the first-born black hole is formed with a negligible spin, we find that all non-zero χeff systems in the detectable population can come only from the common envelope channel as the stable mass-transfer channel cannot shrink the orbits enough for efficient tidal spin-up to take place. We find that the local rate density (z ≃ 0.01) for the common envelope channel is in the range of ∼17–113 Gpc−3 yr−1, considering a range of αCE ∈ [0.2, 5.0], while for the stable mass transfer channel the rate density is ∼25 Gpc−3 yr−1. The latter drops by two orders of magnitude if the mass accretion onto the black hole is not Eddington limited because conservative mass transfer does not shrink the orbit as efficiently as non-conservative mass transfer does. Finally, using GWTC-2 events, we constrained the lower bound of branching fraction from other formation channels in the detected population to be ∼0.2. Assuming all remaining events to be formed through either stable mass transfer or common envelope channels, we find moderate to strong evidence in favour of models with inefficient common envelopes.

AM CVn binaries consist of a WD accreting from a hydrogen-deficient star (or WD) companion (Warner, 1995;Solheim, 2010). In their formation history (Fig. 1.6 and Section 1.3.1.1), AM CVns form after at least one CE phase of their progenitor system. The current RLO is initiated, due to orbital damping caused by GW radiation, at orbital periods of typically 5−20 min (depending on the nature and the temperature of the companion star), and the mass-transfer rate is determined

Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed subpopulation of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ≈10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ≈85% and ≈60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction fB ∼ 0.05 our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift z ∈ [0, 2.5], a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances.