We present the first systematic study of strong binary-single and binary-binary black hole (BH) interactions with the inclusion of general relativity. By including general relativistic effects in the ...equations of motion during strong encounters, the dissipation of orbital energy from the emission of gravitational waves (GWs) can lead to captures and subsequent inspirals with appreciable eccentricities when entering the sensitive frequency ranges of the LIGO and Virgo GW detectors. It has been shown that binary-single interactions significantly contribute to the rate of eccentric mergers, but no studies have looked exclusively into the contribution from binary-binary interactions. To this end, we perform binary-binary and binary-single scattering experiments with general relativistic dynamics up through the 2.5 post-Newtonian order included, both in a controlled setting to gauge the importance of non-dissipative post-Newtonian terms and derive scaling relations for the cross section of GW captures, as well as experiments tuned to the strong interactions from state-of-the art globular cluster (GC) models to assess the relative importance of the binary-binary channel in facilitating GW captures and the resultant eccentricity distributions of inspiral from channel. Although binary-binary interactions are 10-100 times less frequent in GCs than binary-single interactions, their longer lifetime and more complex dynamics leads to a higher probability for GW captures to occur during the encounter. We find that binary-binary interactions contribute 25%-45% of the eccentric mergers that occur during strong BH encounters in GCs, regardless of the properties of the cluster environment. The inclusion of higher multiplicity encounters in dense star clusters therefore have major implications on the predicted rates of highly eccentric binaries potentially detectable by the LIGO/Virgo network. Because gravitational waveforms of eccentric inspirals are distinct from those generated by merging binaries that have circularized, measurements of eccentricity in such systems would highly constrain their formation scenario.
On 2019 August 14, the LIGO and Virgo detectors observed GW190814, a gravitational-wave signal originating from the merger of a black hole (BH) with a compact object. GW190814's compact-binary source ...is atypical both in its highly asymmetric masses and in its lower-mass component lying between the heaviest known neutron star (NS) and lightest known BH in a compact-object binary. If formed through isolated binary evolution, the mass of the secondary is indicative of its mass at birth. We examine the formation of such systems through isolated binary evolution across a suite of assumptions encapsulating many physical uncertainties in massive-star binary evolution. We update how mass loss is implemented for the neutronization process during the collapse of the proto-compact object to eliminate artificial gaps in the mass spectrum at the transition between NSs and BHs. We find it challenging for population modeling to match the empirical rate of GW190814-like systems while simultaneously being consistent with the rates of other compact binary populations inferred from gravitational-wave observations. Nonetheless, the formation of GW190814-like systems at any measurable rate requires a supernova engine model that acts on longer timescales such that the proto-compact object can undergo substantial accretion immediately prior to explosion, hinting that if GW190814 is the result of massive-star binary evolution, the mass gap between NSs and BHs may be narrower or nonexistent.
Abstract
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 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.
The formation of merging binary black holes can occur through multiple astrophysical channels such as, e.g., isolated binary evolution and dynamical formation or, alternatively, have a primordial ...origin. Increasingly large gravitational-wave catalogs of binary black-hole mergers have allowed for the first model selection studies between different theoretical predictions to constrain some of their model uncertainties and branching ratios. In this work, we show how one could add an additional and independent constraint to model selection by using the stochastic gravitational-wave background. In contrast to model selection analyses that have discriminating power only up to the gravitational-wave detector horizons (currently at redshifts
z
≲ 1 for LIGO–Virgo), the stochastic gravitational-wave background accounts for the redshift integration of all gravitational-wave signals in the Universe. As a working example, we consider the branching ratio results from a model selection study that includes potential contribution from astrophysical and primordial channels. We renormalize the relative contribution of each channel to the detected event rate to compute the total stochastic gravitational-wave background energy density. The predicted amplitude lies below the current observational upper limits of GWTC-3 by LIGO–Virgo, indicating that the results of the model selection analysis are not ruled out by current background limits. Furthermore, given the set of population models and inferred branching ratios, we find that, even though the predicted background will not be detectable by current generation gravitational-wave detectors, it will be accessible by third-generation detectors such as the
Einstein
Telescope and space-based detectors such as LISA.
The formation and evolution of binary stars are critical components of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling or merging compact ...binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite, COSMIC, which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signals observable by the Laser Interferometer Space Antenna.
Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic ...population of dynamically formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1 – 2 Gpc−3 yr−1 to the binary merger rate in the local universe, increasing to ≳ 10 Gpc−3 yr−1 at z ∼ 3 . Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture.
Abstract
One of the goals of gravitational-wave astrophysics is to infer the number and properties of the formation channels of binary black holes (BBHs); to do so, one must be able to connect ...various models with the data. We explore benefits and potential issues with analyses using models informed by population synthesis. We consider five possible formation channels of BBHs, as in Zevin et al. (2021b). First, we confirm with the GWTC-3 catalog what Zevin et al. (2021b) found in the GWTC-2 catalog, i.e., that the data are not consistent with the totality of observed BBHs forming in any single channel. Next, using simulated detections, we show that the uncertainties in the estimation of the branching ratios can shrink by up to a factor of ∼1.7 as the catalog size increases from 50 to 250, within the expected number of BBH detections in LIGO–Virgo–KAGRA's fourth observing run. Finally, we show that this type of analysis is prone to significant biases. By simulating universes where all sources originate from a single channel, we show that the influence of the Bayesian prior can make it challenging to conclude that one channel produces all signals. Furthermore, by simulating universes where all five channels contribute but only a subset of channels are used in the analysis, we show that biases in the branching ratios can be as large as ∼50% with 250 detections. This suggests that caution should be used when interpreting the results of analyses based on strongly modeled astrophysical subpopulations.
Abstract
We study the population properties of merging binary black holes in the second LIGO–Virgo Gravitational-Wave Transient Catalog assuming they were all formed dynamically in gravitationally ...bound clusters. Using a phenomenological population model, we infer the mass and spin distribution of first-generation black holes, while self-consistently accounting for hierarchical mergers. Considering a range of cluster masses, we see compelling evidence for hierarchical mergers in clusters with escape velocities ≳100 km s
−1
. For our most probable cluster mass, we find that the catalog contains at least one second-generation merger with 99% credibility. We find that the hierarchical model is preferred over an alternative model with no hierarchical mergers (Bayes factor
>
1400
) and that GW190521 is favored to contain two second-generation black holes with odds
>
700
, and GW190519, GW190602, GW190620, and GW190706 are mixed-generation binaries with
>
10
. However, our results depend strongly on the cluster escape velocity, with more modest evidence for hierarchical mergers when the escape velocity is ≲100 km s
−1
. Assuming that all binary black holes are formed dynamically in globular clusters with escape velocities on the order of tens of km s
−1
, GW190519 and GW190521 are favored to include a second-generation black hole with odds
>
1
. In this case, we find that 99% of black holes from the inferred total population have masses that are less than 49
M
⊙
, and that this constraint is robust to our choice of prior on the maximum black hole mass.
GW190412 is the first observation of a black hole binary with definitively unequal masses. GW190412's mass asymmetry, along with the measured positive effective inspiral spin, allowed for inference ...of a component black hole spin: the primary black hole in the system was found to have a dimensionless spin magnitude between and (90% credible range). We investigate how the choice of priors for the spin magnitudes and tilts of the component black holes affect the robustness of parameter estimates for GW190412, and report Bayes factors across a suite of prior assumptions. Depending on the waveform family used to describe the signal, we find either marginal to moderate ( :1- :1) or strong ( :1) support for the primary black hole being spinning compared to cases where only the secondary is allowed to have spin. We show how these choices influence parameter estimates, and find the asymmetric masses and positive effective inspiral spin of GW190412 to be qualitatively, but not quantitatively, robust to prior assumptions. Our results highlight the importance of both considering astrophysically motivated or population-based priors in interpreting observations and considering their relative support from the data.
Abstract Gravitational-wave observations provide the unique opportunity of studying black hole formation channels and histories—but only if we can identify their origin. One such formation mechanism ...is the dynamical synthesis of black hole binaries in dense stellar systems. Given the expected isotropic distribution of component spins of binary black holes in gas-free dynamical environments, the presence of antialigned or in-plane spins with respect to the orbital angular momentum is considered a tell-tale sign of a merger’s dynamical origin. Even in the scenario where birth spins of black holes are low, hierarchical mergers attain large component spins due to the orbital angular momentum of the prior merger. However, measuring such spin configurations is difficult. Here, we quantify the efficacy of the spin parameters encoding aligned-spin ( χ eff ) and in-plane spin ( χ p ) at classifying such hierarchical systems. Using Monte Carlo cluster simulations to generate a realistic distribution of hierarchical merger parameters from globular clusters, we can infer mergers’ χ eff and χ p . The cluster populations are simulated using Advanced LIGO-Virgo sensitivity during the detector network’s third observing period and projections for design sensitivity. Using a “likelihood-ratio”-based statistic, we find that ∼2% of the recovered population by the current gravitational-wave detector network has a statistically significant χ p measurement, whereas no χ eff measurement was capable of confidently determining a system to be antialigned with the orbital angular momentum at current detector sensitivities. These results indicate that measuring spin-precession through χ p is a more detectable signature of hierarchical mergers and dynamical formation than antialigned spins.