Binary Black Hole Mergers Research Articles

Constraining the AGN formation channel for detected black hole binary mergers up to z=1.5 with the Quaia catalogue

Constraining the AGN formation channel for detected black hole binary mergers up to z=1.5 with the Quaia catalogue

Analytical study of the merging rate of low-mass intermediate-mass black holes in preparation for the future Einstein Telescope

The detection of gravitational wave (GW) signals by Advanced LIGO, Virgo, and KAGRA interferometers opened a new chapter in our understanding of the formation of compact objects. In particular, the detection of GW190521 is observational confirmation of the existence of intermediate-mass black holes (IMBHs); yet more direct observations are needed to better understand the mechanisms behind their formation. In this study, we explore the potential of the next-generation ground-based detector, the Einstein Telescope (ET), to advance our understanding of astrophysics through the detection of GWs emitted by IMBHs. To achieve this, the ET is designed to have improved sensitivity in the low-frequency range of approximately 2-10 Hz, enabling the detection of GWs originating from binary systems containing IMBHs with masses in the range of approximately 10$^2$-10$^5$ $M_ odot We consider black holes (BHs) in the pair-instability form via the hierarchical merger model in galaxies, and approximate the number of events that could be observed by the ET. Our findings indicate that ET could detect a binary black hole (BBH) merger rate of around $2 yr^ $ for BH masses ranging from 10 to 200 M$ odot $, with around 100 $Gpc^ yr^ $ of this rate specifically attributed to BHs in the 100–200 M$ odot $ mass range, which we classify as low-mass IMBHs in this study. This suggests that ET could detect several dozen events similar to GW190521. The exact locations of these BBH mergers are not specified and we count our BH mergers across the entire universe up to a redshift of $z$ approx 2. Observations made with the ET are expected to significantly enhance our comprehension of galactic BH growth, and the existence and characteristics of low-mass IMBHs.

Eccentric mergers in AGN Discs: Influence of the supermassive black-hole on three-body interactions

Abstract There are indications that stellar-origin black holes (BHs) are efficiently paired up in binary black holes (BBHs) in Active Galactic Nuclei (AGN) disc environments, which can undergo interactions with single BHs in the disc. Such binary-single interactions can potentially lead to an exceptionally high fraction of gravitational-wave mergers with measurable eccentricity in LIGO/Virgo/KAGRA. We here take the next important step in this line of studies by performing post-Newtonian N-body simulations between migrating BBHs and single BHs set in an AGN disc-like configuration, with a consistent inclusion of the central supermassive black hole (SMBH) in the equations of motion. With this setup, we study how the fraction of eccentric mergers varies in terms of the initial size of the BBH semi-major axis relative to the Hill sphere, as well as how it depends on the angle between the BBH and the incoming single BH. We find that the fraction of eccentric mergers is still relatively large, even when the interactions are notably influenced by the gravitational field of the nearby SMBH. However, the fraction as a function of the BBH semi-major axis does not follow a smooth functional shape, but instead shows strongly varying features that originate from the underlying phase-space structure. The phase-space further reveals that many of the eccentric mergers are formed through prompt scatterings. Finally, we present the first analytical solution to how the presence of an SMBH in terms of its Hill sphere affects the probability for forming eccentric BBH mergers through chaotic three-body interactions.

The Binary Black Hole Merger Rate Deviates from the Cosmic Star Formation Rate: A Tug of War between Metallicity and Delay Times

Gravitational-wave detectors are now making it possible to investigate how the merger rate of binary black holes (BBHs) evolves with redshift. In this study, we examine whether the BBH merger rate of isolated binaries deviates from a scaled star formation rate density (SFRD)—a frequently used model in state-of-the-art research. To address this question, we conduct population synthesis simulations using COMPAS with a grid of stellar evolution models, calculate their cosmological merger rates, and compare them to a scaled SFRD. We find that our simulated rates deviate by factors up to 3.5 at z ∼ 0 and 5 at z ∼ 9 due to two main phenomena: (i) the formation efficiency of BBHs is an order of magnitude higher at low metallicities than at solar metallicity, and (ii) BBHs experience a wide range of delays (from a few megayears to many gigayears) between formation and merger. The deviations are similar when comparing to a delayed SFRD, and even larger (up to ∼10×) when comparing to SFRD-based models scaled to the local merger rate. Interestingly, our simulations find that the BBH delay time distribution is redshift dependent, increasing the complexity of the redshift distribution of mergers. We find similar results for simulated merger rates of black hole–neutron stars (BHNSs) and binary neutron stars (BNSs). We conclude that the rate of BBH, BHNS, and BNS mergers from the isolated channel can significantly deviate from a scaled SFRD, and that future measurements of the merger rate will provide insights into the formation pathways of gravitational-wave sources.

Dual Active Galactic Nuclei: Precursors of Binary Supermassive Black Hole Formation and Mergers

The presence of dual active galactic nuclei (AGNs) on scales of a few tens of kiloparsecs can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. This study uses the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs (L bol > 1043 erg s−1) and their neighboring SMBHs (any SMBHs closer than 30 pkpc to an AGN) at z ≤ 2. This is our underlying population. We then applied the luminosity threshold of L bol > 1043 erg s−1 to the neighboring SMBHs thereby identifying dual and multiple AGNs. Our findings indicate an increase in the number of both single and dual AGNs from lower to higher redshifts. We also find that the number of dual AGNs with separations of 0.5–4 kpc is twice the number of duals with separations of 4–30 kpc. All dual AGNs in our sample resulted from major mergers. Compared to single AGNs, duals have a lower black hole-to-halo mass ratio. We found that the properties of dual AGN host halos, including halo mass, stellar mass, star formation rate, and gas mass, are generally consistent with those of single AGN halos, albeit tending toward the higher end of their respective property ranges. Our analysis uncovered a diverse array of evolutionary patterns among dual AGNs, including rapidly evolving systems, slower ones, and instances where SMBH mergers are ineffective.

No Evidence for a Dip in the Binary Black Hole Mass Spectrum

Stellar models indicate that the core compactness of a star, which is a common proxy for its explodability in a supernova, does not increase monotonically with the star’s mass. Rather, the core compactness dips sharply over a range of carbon–oxygen core masses; this range may be somewhat sensitive to the star’s metallicity and evolutionary history. Stars in this compactness dip are expected to experience supernovae leaving behind neutron stars, whereas stars on either side of this range are expected to form black holes. This results in a hypothetical mass range in which black holes should seldom form. Quantitatively, when applied to binary stripped stars, these models predict a dearth of binary black holes with component masses ≈10M ⊙–15M ⊙. The population of gravitational-wave signals indicates potential evidence for a dip in the distribution of chirp masses of merging binary black holes near ≈10M ⊙–12M ⊙. This feature could be linked to the hypothetical component mass gap described above, but this interpretation depends on what assumptions are made of the binaries’ mass ratios. Here, we directly probe the distribution of binary black hole component masses to look for evidence of a gap. We find no evidence for this feature using data from the third gravitational-wave transient catalog. If this gap does exist in nature, we find that it is unlikely to be resolvable by the end of the current (fourth) LIGO–Virgo–KAGRA observing run.

A dark standard siren measurement of the Hubble constant following LIGO/Virgo/KAGRA O4a and previous runs

Abstract We present a new constraint on the Hubble constant (H0) from the standard dark siren method using a sample of 5 well-covered gravitational waves (GW) alerts reported during the first part of the fourth LIGO/Virgo/KAGRA observing run and with 3 updated standard dark sirens from third observation run in combination the previous constraints from the first three runs. Our methodology relies on the galaxy catalogue method alone. We use a Deep Learning method to derive the full probability density estimation of photometric redshifts using the Legacy Survey catalogues. We add the constraints from well localized Binary Black Hole mergers to the sample of standard dark sirens analysed in our previous work. We combine the H0 posterior for 5 new standard sirens with other 10 previous events (using the most recent available data for the 5 novel events and updated 3 previous posteriors from O3), finding $H_0 = 70.4^{+13.6}_{-11.7}~{\rm km~s^{-1}~Mpc^{-1}}$ (68% Confidence interval) with the catalogue method only. This result represents an improvement of $\sim 23~{{\%}}$ comparing the new 15 dark siren constrain with the previous 10 dark siren constraint and a reduction in uncertainty of $\sim 40~{{\%}}$ from the combination of 15 dark and bright sirens compared with the GW170817 bright siren alone. The combination of dark and bright siren GW170817 with recent jet constraints yields H0 of $68.0^{+4.4}_{-3.8}~{\rm km~s^{-1}~Mpc^{-1}}$, a $\sim 6~{{\%}}$ precision from Standard Sirens, reducing the previous constraint uncertainty by $\sim 10~{{\%}}$ .

An upper limit on the spins of merging binary black holes formed through isolated binary evolution

As the sensitivity of ground-based gravitational wave detectors progressively increases, observations of black hole mergers will provide us with the joint distribution of their masses and spins. This will be a critical benchmark to validate different formation scenarios. Merging binary black holes formed through the evolution of isolated binary systems require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by the critical rotation limit at their surface, including its deviation from the Keplerian value owing to the outward force provided by radiation. This sets a restriction on their angular momentum content at core-collapse. We aim to determine if this restriction plays a role in the spins of binary black hole mergers. We used detailed calculations of stripped stars with the MESA code at low metallicities ($Z=Z_ $Z_ and $Z_ to determine the dimensionless spins of black holes produced by critically rotating stellar progenitors. To study how such progenitors can arise, we considered their formation through chemically homogeneous evolution (CHE) in binary stars. We used a semi-analytical model to study the physical processes that determine the final angular momentum of CHE binaries, and compared our results against available population synthesis models that rely on detailed binary evolution calculations. We find that above black hole masses of $ 25M_ the dimensionless spin parameter of critically rotating stripped stars ($a=Jc/(GM^2)$) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by isolated binary evolution. CHE can produce binaries where both black holes hit this limit, producing a pileup at the boundary of the excluded region. High-spin black holes arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary black holes detected in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($ eff <0.5$) that merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high-spin population and the existence of a limit set by critical rotation.

Nonlinear studies of modifications to general relativity: Comparing different approaches

Studying the dynamical, nonlinear regime of modified theories of gravity remains a theoretical challenge that limits our ability to test general relativity. Here we consider two generally applicable, but approximate, methods for treating modifications to full general relativity that have been used to study binary black hole mergers and other phenomena in this regime, and compare solutions obtained by them to those from solving the full equations of motion. The first method evolves corrections to general relativity order by order in a perturbative expansion, while the second method introduces extra dynamical fields in such a way that strong hyperbolicity is recovered. We use shift-symmetric Einstein-scalar-Gauss-Bonnet gravity as a benchmark theory to illustrate the differences between these methods for several spacetimes of physical interest. We study the formation of scalar hair about initially nonspinning black holes, the collision of black holes with scalar charge, and the inspiral and merger of binary black holes. By directly comparing predictions, we assess the extent to which those from the approximate treatments can be meaningfully confronted with gravitational wave observations. We find that the order-by-order approach cannot faithfully track the solutions when the corrections to general relativity are non-negligible. The second approach, however, can provide consistent solutions, provided the timescale over which the dynamical fields are driven to their target values is made short compared to the physical timescales. Published by the American Physical Society 2024

A trifecta of modelling tools: a Bayesian binary black hole model selection combining population synthesis and galaxy formation models

ABSTRACT Gravitational waves (GWs) have revealed surprising properties of binary black hole (BBH) populations, but there is still mystery surrounding how these compact objects evolve. We apply Bayesian inference and an efficient method to calculate the BBH merger rates in the Shark host galaxies, to determine the combination of COMPAS parameters that outputs a population most like the GW sources from the LIGO, Virgo, and KAGRA (LVK) transient catalogue. For our COMPAS models, we calculate the likelihood with and without the dependence on the predicted number of BBH merger events. We find strong correlations between hyper-parameters governing the specific angular momentum (AM) of mass lost during mass transfer, the mass-loss rates of Wolf–Rayet stars via winds and the chemically homogeneous evolution (CHE) formation channel. We conclude that analysing the marginalized and unmarginalized likelihood is a good indicator of whether the population parameters distribution and number of observed events reflect the LVK data. In doing so, we see that the majority of the models preferred in terms of the population-level parameters of the BBHs greatly overpredict the number of events we should have observed to date. Looking at the smaller number of models that perform well with both likelihoods, we find that those with no CHE, AM loss occurring closer to the donor during the first mass-transfer event, and/or higher rates of mass-loss from Wolf–Rayet winds are generally preferred by current data. We find these conclusions to be robust to our choice of selection criteria.

Black hole—neutron star binary mergers: the impact of stellar compactness

Abstract Recent gravitational wave (GW) observations include possible detections of black hole—neutron star binary mergers. As with binary black hole mergers, numerical simulations help characterize the sources. For binary systems with neutron star components, the simulations help to predict the imprint of tidal deformations and disruptions on the GW signals. In a previous study, we investigated how the mass of the black hole has an impact on the disruption of the neutron star and, as a consequence, on the shape of the GWs emitted. We extend these results to study the effects of varying the compactness of the neutron star. We consider neutron star compactness in the 0.113–0.2 range for binaries with mass ratios of 3 and 5. As the compactness and the mass ratio increase, the binary system behaves during the late inspiral and merger more like a black hole binary. For the cases with the least compact neutron star, the GWs emitted, in terms of mismatches, are the most distinguishable from those by a binary black hole. The disruption of the star significantly suppresses the kicks on the final black hole. The disruption also affects, although not dramatically, the spin of the final black hole. Lastly, for neutron stars with low compactness, the quasi-normal ringing of the black hole after the merger does not show a clean quasi-normal ringing because of the late accretion of debris from the neutron star.

Binary black hole mergers from Population III star clusters

Binary black holes (BBHs) that are born from the evolution of Population III (Pop. III) stars are one of the main high-redshift targets for next-generation ground-based gravitational-wave (GW) detectors. Their predicted initial mass function and lack of metals make them the ideal progenitors of black holes above the upper edge of the pair-instability mass gap, that is, with a mass higher than ~134 (241) M⊙ for stars that become (or do not become) chemically homogeneous during their evolution. We investigate the effects of cluster dynamics on the mass function of BBHs that are born from Pop. III stars by considering the main uncertainties on the mass function of Pop. III stars, the orbital properties of the binary systems, the star cluster mass, and the disruption time. In our dynamical models, at least ~5% and up to 100% BBH mergers in Pop. III star clusters have a primary mass m1 above the upper edge of the pair-instability mass gap. In contrast, only ≲3% isolated BBH mergers have a primary mass above the gap, unless their progenitors evolved as chemically homogeneous stars. The lack of systems with a primary and/or secondary mass inside the gap defines a zone of avoidance with sharp boundaries in the plane of the primary mass-mass ratio. Finally, we estimate the merger rate density of BBHs. In the most optimistic case, we find a maximum of ℛ ~ 200 Gpc-3 yr-1 at ɀ ~ 15 for BBHs that formed via dynamical capture. For comparison, the merger rate density of isolated Pop. III BBHs is ℛ ≤ 10 Gpc−3 yr−1 for the same model of Pop. III star formation history.

Standard siren cosmology with gravitational waves from binary black hole mergers in active galactic nuclei

Standard siren cosmology with gravitational waves from binary black hole mergers in active galactic nuclei

Greybody factors imprinted on black hole ringdowns. II. Merging binary black holes

Greybody factors imprinted on black hole ringdowns. II. Merging binary black holes

Navigating Unknowns: Deep Learning Robustness for Gravitational-wave Signal Reconstruction

We present a rapid and reliable deep-learning-based method for gravitational-wave (GW) signal reconstruction from elusive, generic binary black hole mergers in LIGO data. We demonstrate that our model, AWaRe, effectively recovers GWs with parameters it has not encountered during training. This includes features like higher black hole masses, additional harmonics, eccentricity, and varied waveform systematics, which introduce complex modulations in the waveform’s amplitudes and phases. The accurate reconstructions of these unseen signal characteristics demonstrate AWaRe's ability to handle complex features in the waveforms. By directly incorporating waveform reconstruction uncertainty estimation into the AWaRe framework, we show that for real GW events, the uncertainties in AWaRe's reconstructions align closely with those achieved by benchmark algorithms like BayesWave and coherent WaveBurst. The robustness of our model to real GW events and its ability to extrapolate to unseen data open new avenues for investigations in various aspects of GW astrophysics and data analysis, including tests of general relativity and the enhancement of current GW search methodologies.

The effect of stellar rotation on black hole mass and spin

Abstract The gravitational wave signature of a binary black hole (BBH) merger is dependent on its component mass and spin. If such black holes originate from rapidly rotating progenitors, the large angular momentum reserve in the star could drive a collapsar-like supernova explosion, hence substantially impacting these characteristics of the black holes in the binary. To examine the effect of stellar rotation on the resulting black hole mass and spin, we conduct a 1D general relativistic study of the end phase of the stellar collapse. We find that the resulting black hole mass at times differs significantly from the previously assumed values. We quantify the dependence of the black hole spin magnitude on the hydrodynamics of the accretion flow, providing analytical relations for calculating the mass and spin based on the progenitor’s pre-collapse properties. Depending on the nature of the accretion flow, our findings have implications for the black hole upper mass gap resulting from pair-instability supernovae, the maximum mass of a maximally rotating stellar black hole, and the maximum effective spin of a BBH formed in a tidally locked helium star - black hole binary.

Measuring kinematic anisotropies with pulsar timing arrays

Recent pulsar timing array (PTA) collaborations show strong evidence for a stochastic gravitational wave background (SGWB) with the characteristic Hellings-Downs interpulsar correlations. The signal may stem from supermassive black hole binary mergers, or early Universe phenomena. The former is expected to be strongly anisotropic, while primordial backgrounds are likely to be predominantly isotropic with small fluctuations. In the case the observed SGWB is of cosmological origin, our relative motion with respect to the SGWB rest frame is a guaranteed source of anisotropy, leading to O(10−3) energy density fluctuations of the SGWB. For such cosmological SGWB, kinematic anisotropies are likely to be larger than the intrinsic anisotropies, akin to the cosmic microwave background (CMB) dipole anisotropy. We assess the sensitivity of current PTA data to the kinematic dipole anisotropy, and we also forecast at what extent the magnitude and direction of the kinematic dipole can be measured in the future with an SKA-like experiment. We also discuss how the spectral shape of the SGWB and the location of the pulsars to monitor affect the prospects of detecting the kinematic dipole with PTA. In the future, a detection of this anisotropy may even help resolve the discrepancy in the magnitude of the kinematic dipole as measured by CMB and large-scale structure observations. Published by the American Physical Society 2024

Differentiable and hardware-accelerated waveforms for gravitational wave data analysis

We propose the use of automatic differentiation through the programming framework for accelerating a variety of analysis tasks throughout gravitational wave (GW) science. Firstly, we demonstrate that complete waveforms which cover the inspiral, merger, and ringdown of binary black holes (i.e., IMRPhenomD) can be written in and demonstrate that the serial evaluation speed of the waveform (and its derivative) is similar to the implementation in . Moreover, allows for graphics processing unit–accelerated waveform calls which can be over an order of magnitude faster than serial evaluation on a CPU. We then focus on three applications where efficient and differentiable waveforms are essential. Firstly, we demonstrate how gradient descent can be used to optimize the ∼200 coefficients that are used to calibrate the waveform model. In particular, we demonstrate that the typical with numerical relativity waveforms can be improved by more than 50%. Secondly, we show that Fisher forecasting calculations can be sped up by ∼3–5× (on a CPU) with no loss in accuracy. This increased speed makes Fisher forecasting for a population of events substantially simpler. Finally, we show that gradient-based samplers like Hamiltonian Monte Carlo lead to significantly reduced autocorrelation values when compared to traditional Monte Carlo methods. Since differentiable waveforms have substantial advantages for a variety of tasks throughout GW science, we propose that waveform developers use to build new waveforms moving forward. Our waveform code, , can be found on GitHub website and will continue to be updated with new waveforms as they are implemented. Published by the American Physical Society 2024

Close Encounters of Wide Binaries Induced by the Galactic Tide: Implications for Stellar Mergers and Gravitational-wave Sources

A substantial fraction of stars can be found in wide binaries with projected separations between ∼102 and 105 au. In the standard lore of binary physics, these would evolve as effectively single stars that remotely orbit one another on stationary Keplerian ellipses. However, embedded in their Galactic environment, the low binding energy of wide binaries makes them exceptionally prone to perturbations from the gravitational potential of the Milky Way and encounters with passing stars. Employing a fully relativistic N-body integration scheme, we study the impact of these perturbations on the orbital evolution of wide binaries along their trajectory through the Milky Way. Our analysis reveals that the torques exerted by the Galaxy can cause large-amplitude oscillations of the binary eccentricity to 1 − e ≲ 10−8. As a consequence, the wide binary members pass close to each other at periapsis, which, depending on the type of binary, potentially leads to a mass transfer or collision of stars or to an inspiral and subsequent merger of compact remnants due to gravitational-wave radiation. Based on a simulation of 105 wide binaries across the Galactic field, we find that this mechanism could significantly contribute to the rate of stellar collisions and binary black hole mergers as inferred from observations of luminous red novae and gravitational-wave events by LIGO/Virgo/Kagra. We conclude that the dynamics of wide binaries, despite their large mean separation, can give rise to extreme interactions between stars and compact remnants.

Unequal-mass highly spinning binary black hole mergers in the stable mass transfer formation channel

Unequal-mass highly spinning binary black hole mergers in the stable mass transfer formation channel