Binary Black Hole Mergers Research Articles

Binary black hole mergers from merged stars in the Galactic field

The majority of massive stars are found in close binaries which: (i) are prone to merge and (ii) are accompanied by another distant tertiary star (triples). Here, we study the evolution of the stellar post-merger binaries composed of the merger product and the tertiary companion. We find that post-merger binaries originating from compact stellar triples with outer semi-major axes $a_{\rm out,init}\lesssim10^1\,-\,10^2\,\rm AU$ provide a new way to form binary black hole mergers in the galactic field. By means of a population synthesis, we estimate their contribution to the total black hole merger rate to be $\mathcal{R}(z=0)=0.3\,-\,25.2\,\rm Gpc^{-3}\,yr^{-1}$. Merging binary black holes that form from stellar post-merger binaries have exceptionally low mass ratios. We identify a critical mass ratio $q\simeq0.5$ below which they dominate the total black hole merger rate in the field. We show that after including their additional contribution, the mass ratio distribution of binary black hole mergers in the galactic field scenario is in better agreement with that inferred from gravitational wave detections.

Signs of higher multipoles and orbital precession in GW151226

We present a reanalysis of GW151226, the second binary black hole merger discovered by the LIGO--Virgo Collaboration. Previous analysis showed that the best-fit waveform for this event corresponded to the merger of a $\sim 14 \, M_\odot$ black hole with a $\sim 7.5 \, M_\odot$ companion, and the posterior distribution in mass ratio ($q \leq 1$) is rather flat. In this work, we perform parameter estimation using a waveform model that includes the effects of orbital precession and higher-order radiative multipole modes, and we find that the source parameters of GW151226 shift towards the low $q$ and high effective spin ($\chi_{\rm eff}$) region and that $q$ is better measured. The new solution has a log likelihood roughly two points higher than when either higher multipoles or orbital precession is neglected and can alter the astrophysical interpretation of GW151226. Additionally, we find it useful to use a flat-in-$\chi{\rm eff}$ prior, which does not penalize the large $|\chi_{\rm eff}|$ region, in order to uncover the higher likelihood region for GW151226. Our solution has several interesting properties: (a) the secondary black hole mass is close to the upper limit of the hypothesized lower mass gap of astrophysical black hole population; and (b) orbital precession is driven by the primary black hole spin, which has a dimensionless magnitude as large as $\sim 0.85$ and is tilted away from the orbital angular momentum at an angle of $\sim 57^\circ$. Since GW151226 is a relatively weak signal, an unambiguous claim of the detection of these effects in the signal cannot be made.

The Gravitational Wave Universe Toolbox

Context:TheGW-Universe Toolboxis a software package that simulates observations of the gravitational wave (GW) Universe with different types of GW detectors, including Earth-based and space-borne laser interferometers and pulsar timing arrays. It is accessible as a website, and can also be imported and run locally as a Python package.Methods:We employ the method used by theGW-Universe Toolboxto generate a synthetic catalogue of detection of stellar-mass binary black hole (BBH) mergers. As an example of its scientific application, we study how GW observations of BBHs can be used to constrain the merger rate as a function of redshift and masses. We study advanced LIGO (aLIGO) and theEinsteinTelescope (ET) as two representatives of the second and third generation GW observatories, respectively. We also simulate the observations from a detector that is half as sensitive as the ET at its nominal designed sensitivity, which represents an early phase of the ET. We used two methods to obtain the constraints on the source population properties from the catalogues: the first uses a parameteric differential merger rate model and applies a Bayesian inference on the parameters; the other is non-parameteric and uses weighted Kernel density estimators.Results:Our results show the overwhelming advantages of the third generation detector over those of the second generation for the study of BBH population properties, especially at redshifts higher than ∼2, where the merger rate is believed to peak. With the simulated aLIGO catalogue, the parameteric Bayesian method can still give some constraints on the merger rate density and mass function beyond its detecting horizon, while the non-parametric method loses the constraining ability completely there. The difference is due to the extra information placed by assuming a specific parameterisation of the population model in the Bayesian method. In the non-parameteric method, no assumption of the general shape of the merger rate density and mass function are placed, not even the assumption of its smoothness. These two methods represent the two extreme situations of general population reconstruction. We also find that, despite the numbers of detected events of the half ET can easily be compatible with full ET after a longer observation duration, and the catalogue from the full ET can still give much better constraints on the population properties due to its smaller uncertainties on the physical parameters of the GW events.

The QCD phase transition behind a PBH origin of LIGO/Virgo events?

The best-motivated scenario for a sizable primordial black hole (PBH) contribution to the LIGO/Virgo binary black hole mergers invokes the QCD phase transition, which naturally enhances the probability to form PBH with masses of stellar scale. We reconsider the expected mass function associated not only to the QCD phase transition proper, but also the following particle antiparticle annihilation processes, and analyse the constraints on this scenario from a number of observations: the specific pattern in cosmic microwave background (CMB) anisotropies induced by accretion onto PBHs, CMB spectral distortions, gravitational wave searches, and direct counts of supermassive black holes (SMBHs) at high redshift. We find that the scenario is not viable, unless an ad hoc mass evolution for the PBH mass function and a cutoff in power-spectrum very close to the QCD scale are introduced by hand. Despite these negative results, we note that a future detection of coalescing binaries involving sub-solar PBHs has the potential to check the cosmological origin of SMBHs at the e ± annihilation epoch, if indeed the PBH mass function is shaped by the changes to the equation of state driven by the thermal history of the universe.

The Gravitational Wave Universe Toolbox

Context.As the importance of gravitational wave (GW) astrophysics increases rapidly, astronomers interested in GWs who are not experts in this field sometimes need to get a quick idea of what GW sources can be detected by certain detectors, and the accuracy of the measured parameters.Aims.The GW-Toolbox is a set of easy-to-use, flexible tools to simulate observations of the GW universe with different detectors, including ground-based interferometers (advanced LIGO, advanced VIRGO, KAGRA, Einstein Telescope, Cosmic Explorer, and also customised interferometers), space-borne interferometers (LISA and a customised design), and pulsar timing arrays mimicking the current working arrays (EPTA, PPTA, NANOGrav, IPTA) and future ones. We include a broad range of sources, such as mergers of stellar-mass compact objects, namely black holes, neutron stars, and black hole–neutron star binaries, supermassive black hole binary mergers and inspirals, Galactic double white dwarfs in ultra-compact orbit, extreme-mass-ratio inspirals, and stochastic GW backgrounds.Methods.We collected methods to simulate source populations and determine their detectability with various detectors. Our aim is to provide a comprehensive description of the methodology and functionality of the GW-Toolbox.Results.The GW-Toolbox produces results that are consistent with previous findings in the literature, and the tools can be accessed via a website interface or as a Python package. In the future, this package will be upgraded with more functions.

Are Binary Black Hole Mergers and Long Gamma-Ray Bursts Drawn from the Same Black Hole Population?

This paper compares the population of binary black hole (BBH) mergers detected by LIGO/Virgo with selected long gamma-ray burst (GRB) world models convolved with a delay function (LGRBs are used as a tracer of stellar-mass BH formation). The comparison involves the redshift distribution and the fraction of LGRBs required to produce the local rate of BBH mergers. We find that BBH mergers and LGRBs cannot have the same formation history, unless BBH mergers have a long coalescence time of several Gyr. This would imply that BHs born during the peak of long GRB formation at redshift z ≈ 2−3 merge within the horizon of current GW interferometers. We also show that LGRBs are more numerous than BBH mergers, meaning that most of them do not end their lives in BBH mergers. We interpret these results as an indication that BBH mergers and LGRBs constitute two distinct populations of stellar-mass BHs, with LGRBs being more frequent than BBH mergers. We speculate that the descendants of LGRBs may resemble galactic high-mass X-ray binaries more than BBH mergers. Finally, we discuss the possible existence of a subpopulation of fast-spinning LGRB descendants among BBH mergers, showing that this population, if it exists, is expected to become dominant beyond redshift z ≈ 1, leading to a change in the observed properties of BBH mergers.

MOCCA-SURVEY data base II – Properties of intermediate mass black holes escaping from star clusters

ABSTRACT In this work, we investigate properties of intermediate-mass black holes (IMBHs) that escape from star clusters due to dynamical interactions. The studied models were simulated as part of the preliminary second survey carried out using the MOCCA code (MOCCA-SURVEY Database II), which is based on the Monte Carlo N-body method and does not include gravitational wave recoil kick prescriptions of the binary black hole merger product. We have found that IMBHs are more likely to be formed and ejected in models where both initial central density and central escape velocities have high values. Most of our studied objects escape in a binary with another black hole (BH) as their companion and have masses between 100 and $\rm 140 {\rm ~M}_{\odot }$. Escaping IMBHs tend to build-up mass most effectively through repeated mergers in a binary with BHs due to gravitational wave emission. Binaries play a key role in their ejection from the system as they allow these massive objects to gather energy needed for escape. The binaries in which IMBHs escape tend to have very high binding energy at the time of escape and the last interaction is strong but does not involve a massive intruder. These IMBHs gain energy needed to escape the cluster gradually in successive dynamical interactions. We present specific examples of the history of IMBH formation and escape from star cluster models. We also discuss the observational implications of our findings as well as the potential influence of the gravitational wave recoil kicks on the process.

Searching for mass-spin correlations in the population of gravitational-wave events: The GWTC-3 case study

One fundamental goal of the newly born gravitational wave astronomy is discovering the origin of the observed binary black hole mergers. Towards this end, identifying features in the growing wealth of data may help in distinguishing different formation pathways. While large uncertainties still affect the binary formation models, spin-mass relations remain characteristic features of specific classes of channels. By focusing on the effective inspiral spin $\chi_\text{eff}$, the best reconstructed spin-related merger parameter, we show that current GWTC-3 data support the hypothesis that a fraction of events may display mass-spin correlations similar to one expected by dynamical formation channels of either astrophysical or primordial nature. We quantify the Bayesian evidence in favour of those models, which are substantially preferred when compared to the Gaussian phenomenological model adopted to describe the distribution of $\chi_\text{eff}$ in the recent LIGO/Virgo/KAGRA population analyses.

Accelerating binary black hole system

Recently, by the end of their second observing run, advanced LIGO and Virgo have detected some binary black hole mergers. In this paper, we show that the acceleration of black holes around each other causes the emergence of two Rindler space–times. The parameters of the black hole in each region change in opposite order to the parameters of the black hole in the other region. Also, the acceleration has a direct relation to the black hole entropies and their observed redshifts.

Flexible and fast estimation of binary merger population distributions with an adaptive kernel density estimator

The LIGO Scientific, Virgo and KAGRA Collaborations recently released the third gravitational wave transient catalog or GWTC-3, significantly expanding the number of gravitational wave signals. To address the---still uncertain---formation channels of the source compact binaries, their population properties must be characterized. The computational cost of the Bayesian hierarchical methods employed thus far scales with the size of the event catalogs, and such methods have until recently assumed fixed functional forms for the source distribution. Here we propose a fast and flexible method to reconstruct the population of LIGO-Virgo merging black hole (BH) binaries without such assumptions. For sufficiently high event statistics and sufficiently low individual event measurement error (relative to the scale of population features) a kernel density estimator (KDE) reconstruction of the event distribution will be accurate. We improve the accuracy and flexibility of KDE for finite event statistics using an adaptive bandwidth KDE (awKDE). We apply awKDE to publicly released parameter estimates for 44 significant (69) BH binary mergers in GWTC-2 (GWTC-3), in combination with a fast polynomial fit of search sensitivity, to obtain a nonparametric estimate of the mass distribution, and compare to Bayesian hierarchical methods. We also demonstrate a robust peak detection algorithm based on awKDE and use it to calculate the significance of the apparent peak in the BH mass distribution around $35\text{ }\text{ }{\mathrm{M}}_{\ensuremath{\bigodot}}$. We find such a peak is very unlikely to have occurred if the true distribution is a featureless power-law (significance of $3.6\ensuremath{\sigma}$ for confident GWTC-2 BBH events, $3.0\ensuremath{\sigma}$ for confident GWTC-3 BBH events).

LIGO–Virgo correlations between mass ratio and effective inspiral spin: testing the active galactic nuclei channel

ABSTRACT Observations by LIGO–Virgo of binary black hole mergers suggest a possible anticorrelation between black hole mass ratio (q = m2/m1) and the effective inspiral spin parameter χeff, the mass-weighted spin projection on to the binary orbital angular momentum. We show that such an anticorrelation can arise for binary black holes assembled in active galactic nuclei (AGNs) due to spherical and planar symmetry-breaking effects. We describe a phenomenological model in which (1) heavier black holes live in the AGN disc and tend to spin-up into alignment with the disc; (2) lighter black holes with random spin orientations live in the nuclear spheroid; (3) the AGN disc is dense enough to rapidly capture a fraction of the spheroid component, but small in radial extent to limit the number of bulk disc mergers; (4) migration within the disc is non-uniform, likely disrupted by feedback from migrators or disc turbulence; (5) dynamical encounters in the disc are common and preferentially disrupt binaries that are retrograde around their centre of mass, particularly at stalling orbits, or traps. Comparisons of predictions in (q, χeff) parameter space for the different channels may allow us to distinguish their fractional contributions to the observed merger rates.

Black hole binaries in cubic Horndeski theories

We study black hole binary mergers in certain cubic Horndeski theories of gravity, treating them fully nonlinearly. In the regime of validity of effective field theory, the mismatch of the gravitational wave strain between Horndeski and general relativity (coupled to a scalar field) can be as large as 10%--13% in the Advanced LIGO mass range. The initial data and coupling constants are chosen such the theory remains in the weakly coupled regime throughout the evolution. In all cases that we have explored, we observe that the waveform in the Horndeski theory is shifted by an amount much larger than the smallness parameter that controls the initial data. This effect is generic and may be present in other theories of gravity involving higher derivatives.

Partial-envelope stripping and nuclear-timescale mass transfer from evolved supergiants at low metallicity

Stable mass transfer from a massive post-main sequence (post-MS) donor is thought to be a short-lived event of thermal-timescale mass transfer (∼10−3 M⊙ yr−1) which within ≲104 yr strips the donor star of nearly its entire H-rich envelope, producing a hot, compact helium star. This long-standing picture is based on stellar models with rapidly expanding Hertzprung gap (HG) donor stars. Motivated by a finding that in low-metallicity binaries, post-MS mass transfer may instead be initiated by donors already at the core-helium burning (CHeB) stage, we used the MESA stellar-evolution code to compute grids of detailed massive binary models at three metallicities: those of the Sun, the Large Magellanic Cloud (LMC, ZFe; LMC/ZFe; ⊙ ≈ 0.36), and the Small Magellanic Cloud (SMC, ZFe; SMC/ZFe; ⊙ ≈ 0.2). Our grids span a wide range in orbital periods (∼3 to 5000 days) and initial primary masses (10 M⊙ to 36 − 53 M⊙, depending on metallicity). We find that metallicity strongly influences the course and outcome of mass-transfer evolution. We identify two novel types of post-MS mass transfer: (a) mass exchange on the long nuclear timescale (ΔTMT ≳ 105 yr, Ṁ ∼ 10−5 M⊙ yr−1) that continues until the end of the CHeB phase, and (b) rapid mass transfer leading to detached binaries with mass losers that are only partially stripped of their envelopes. At LMC and SMC compositions, the majority of binary models with donor masses ≥17 M⊙ follow one of these two types of evolution. In neither (a) nor (b) does the donor become a fully stripped helium star by the end of CHeB. Boundaries between the different types of post-MS mass transfer evolution are associated with the degree of rapid post-MS expansion of massive stars and, for a given metallicity, are sensitive to the assumptions about internal mixing. At low metallicity, due to partial envelope stripping, we predict fewer hot fully stripped stars formed through binary interactions as well as higher compactness of the presupernova core structures of mass losers. Nuclear-timescale post-MS mass transfer suggests a strong preference for metal-poor host galaxies of ultra-luminous X-ray sources with black-hole (BH) accretors and massive donors, some of which might be the immediate progenitors of binary BH mergers. It also implies a population of interacting binaries with blue and yellow supergiant donors. Partially stripped stars could potentially explain the puzzling nitrogen-enriched slowly rotating (super)giants in the LMC.

No slip gravity in light of LISA standard sirens

ABSTRACT Standard sirens (SS) are the gravitational wave analogue of the astronomical standard candles, and can provide powerful information about the dynamics of the Universe up to very high z values. In this work, we generate three mock SS catalogues based on the merger of massive black hole binaries which are expected to be observed in the LISA operating frequency band. Then, we perform an analysis to test modifications of general relativity (GR) inspired by the no slip gravity framework. We find that in the best scenarios, we can constrain the free parameters which quantify deviations from GR to 21 per cent accuracy, while the Hubble parameter can be simultaneously fit to 6 per cent accuracy. In combination with CMB information, we find a 15 per cent accuracy on the modified gravity free parameters and 0.7 per cent accuracy on the Hubble parameter. The SS events at very large cosmological distances to be observed in LISA band will provide a unique way to test nature of gravity, but in the context of the analysis performed here, it will not be possible to distinguish the no slip gravity from GR.

The Imprint of Superradiance on Hierarchical Black Hole Mergers

Ultralight bosons are a proposed solution to outstanding problems in cosmology and particle physics: they provide a dark-matter candidate while potentially explaining the strong charge-parity problem. If they exist, ultralight bosons can interact with black holes through the superradiant instability. In this work we explore the consequences of this instability on the evolution of hierarchical black holes within dense stellar clusters. By reducing the spin of individual black holes, superradiance reduces the recoil velocity of merging binary black holes, which, in turn, increases the retention fraction of hierarchical merger remnants. We show that the existence of ultralight bosons with mass 2 × 10−14 ≲ μ/eV ≲ 2 × 10−13 would lead to an increased rate of hierarchical black hole mergers in nuclear star clusters. An ultralight boson in this energy range would result in up to ≈60% more present-day nuclear star clusters supporting hierarchical growth. The presence of an ultralight boson can also double the rate of intermediate-mass black hole mergers to ≈0.08 Gpc−3 yr−1 in the local universe. These results imply that a select range of ultralight boson masses can have far-reaching consequences for the population of black holes in dense stellar environments. Future studies into black hole cluster populations and the spin distribution of hierarchically formed black holes will test this scenario.

Measuring the Properties of Active Galactic Nuclei Disks with Gravitational Waves

Active galactic nuclei (AGNs) are promising environments for the assembly of merging binary black hole (BBH) systems. Interest in AGNs as nurseries for merging BBHs is rising, following the detection of gravitational waves from a BBH system from the purported pair-instability mass gap, most notably GW190521. AGNs have also been invoked to explain the formation of the high-mass-ratio system GW190814. We draw on simulations of BBH systems in AGNs to propose a phenomenological model for the distribution of black hole spins of merging binaries in AGN disks. The model incorporates distinct features that make the AGN channel potentially distinguishable from other channels, such as assembly in the field and in globular clusters. The model parameters can be mapped heuristically to the age and density of the AGN disks. We estimate the extent to which different populations of mergers in AGNs can be distinguished. If the majority of merging black holes are assembled in AGNs, future gravitational-wave observations may provide insights into the dynamics of AGN disks.

Detectability of a spatial correlation between stellar mass black hole mergers and active galactic nuclei in the local Universe

ABSTRACT The origin of the binary black hole (BBH) mergers detected through gravitational waves (GWs) by the LIGO-Virgo-KAGRA (LVK) collaboration remains debated. One fundamental reason is our ignorance of their host environment, as the typical size of an event’s localization volume can easily contain thousands of galaxies. A strategy around this is to exploit statistical approaches to assess the spatial correlation between these mergers and astrophysically motivated host galaxy types, such as active galactic nuclei (AGNs). We use a likelihood ratio method to infer the degree of GW–AGN connection out to z = 0.2. We simulate BBH mergers whose components’ masses are sampled from a realistic distribution of the underlying population of black holes (BHs). Localization volumes for these events are calculated assuming two different interferometric network configurations. These correspond to the configuration of the third (O3) and of the upcoming fourth (O4) LVK observing runs. We conclude that the 13 BBH mergers detected during the third observing run at z ≤ 0.2 are not enough to reject with a 3σ significance the hypothesis according to which there is no connection between GW and AGNs more luminous than $\approx 10^{44.3}\rm {erg}\ \rm {s}^{-1}$, that have number density higher than 10−4.75 Mpc−3. However, 13 detections are enough to reject this no-connection hypothesis when rarer categories of AGNs are considered, with bolometric luminosities greater than $\approx 10^{45.5}\rm {erg}\ \rm {s}^{-1}$. We estimate that O4 results will potentially allow us to test fractional contributions to the total BBH merger population from AGNs of any luminosity higher than $80{{\ \rm per\ cent}}$.

Evidence of Large Recoil Velocity from a Black Hole Merger Signal.

The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick," reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to v_{f}∼1542_{-1098}^{+747} km/s or v_{f}≳698 km/s (one-sided limit), at 90% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than 0.48% (7.7%) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors.

Gravitational-wave echoes from numerical-relativity waveforms via spacetime construction near merging compact objects

We propose a new approach toward reconstructing the late-time near-horizon geometry of merging binary black holes, and toward computing gravitational-wave echoes from exotic compact objects. A binary black-hole merger spacetime can be divided by a time-like hypersurface into a Black-Hole Perturbation (BHP) region, in which the space-time geometry can be approximated by homogeneous linear perturbations of the final Kerr black hole, and a nonlinear region. At late times, the boundary between the two regions is an infalling shell. The BHP region contains late-time gravitational-waves emitted toward the future horizon, as well as those emitted toward future null infinity. In this region, by imposing no-ingoing wave conditions at past null infinity, and matching out-going waves at future null infinity with waveforms computed from numerical relativity, we can obtain waves that travel toward the future horizon. In particular, the Newman-Penrose $\psi_0$ associated with the in-going wave on the horizon is related to tidal deformations measured by fiducial observers floating above the horizon. We further determine the boundary of the BHP region on the future horizon by imposing that $\psi_0$ inside the BHP region can be faithfully represented by quasi-normal modes. Using a physically-motivated way to impose boundary conditions near the horizon, and applying the so-called Boltzmann reflectivity, we compute the quasi-normal modes of non-rotating ECOs, as well as gravitational-wave echoes. We also investigate the detectability of these echoes in current and future detectors, and prospects for parameter estimation.

Stable Mass Transfer Can Explain Massive Binary Black Hole Mergers with a High-spin Component

Recent gravitational-wave observations showed that binary black hole (BBH) mergers with massive components are more likely to have high effective spins. In the model of isolated binary evolution, BH spins mainly originate from the angular momenta of the collapsing cores before BH formation. Both observations and theories indicate that BHs tend to possess relatively low spins; the origin of fast-spinning BHs remains a puzzle. We investigate an alternative process that stable Case A mass transfer may significantly increase BH spins during the evolution of massive BH binaries. We present detailed binary evolution calculations and find that this process can explain the observed high spins of some massive BBH mergers under the assumption of mildly super-Eddington accretion.