Binary Black Hole Mergers Research Articles

Constraining extended cosmologies with GW×LSS cross-correlations

The rapid development of gravitational wave astronomy provides the unique opportunity of exploring the dynamics of the Universe using clustering properties of coalescing binary black hole mergers. Gravitational wave data, along with information coming from future galaxy surveys, have the potential of shedding light about many open questions in Cosmology, including those regarding the nature of dark matter and dark energy. In this work we explore which combination of gravitational wave and galaxy survey datasets are able to provide the best constraints both on modified gravity theories and on the nature of the very same binary black hole events. In particular, by using the public Boltzmann code Multi_CLASS, we compare cosmological constraints on popular ΛCDM extensions coming from gravitational waves alone and in conjunction with either deep and localized or wide and shallow galaxy surveys. We show that constraints on extensions of General Relativity will be at the same level of existing limits from gravitational waves alone or one order of magnitude better when galaxy surveys are included. Furthermore, cross-correlating both kind of galaxy survey with gravitational waves datasets will allow to confidently rule in or out primordial black holes as dark matter candidate in the majority of the allowed parameter space.

Supernova calibration by gravitational waves

Hubble tension is one of the most important problems in cosmology. Although the local measurements on the Hubble constant with Type Ia supernovae (SNe Ia) are independent of cosmological models, they suffer the problem of zero-point calibration of the luminosity distance. The observations of gravitational waves (GWs) with space-based GW detectors can measure the luminosity distance of the GW source with high precision. By assuming that massive binary black hole mergers and SNe Ia occur in the same host galaxy, we study the possibility of re-calibrating the luminosity distances of SNe Ia by GWs. Then we use low-redshift re-calibrated SNe Ia to determine the local Hubble constant. We find that we need at least 7 SNe Ia with their luminosity distances re-calibrated by GWs to reach a 2% precision of the local Hubble constant. The value of the local Hubble constant is free from the problems of zero-point calibration and model dependence, so the result can shed light on the Hubble tension.

Random Convolutional Kernels for Space-Detector Based Gravitational Wave Signals

Neural network models have entered the realm of gravitational wave detection, proving their effectiveness in identifying synthetic gravitational waves. However, these models rely on learned parameters, which necessitates time-consuming computations and expensive hardware resources. To address this challenge, we propose a gravitational wave detection model tailored specifically for binary black hole mergers, inspired by the Random Convolutional Kernel Transform (ROCKET) family of models. We conduct a rigorous analysis by factoring in realistic signal-to-noise ratios in our datasets, demonstrating that conventional techniques lose predictive accuracy when applied to ground-based detector signals. In contrast, for space-based detectors with high signal-to-noise ratios, our method not only detects signals effectively but also enhances inference speed due to its streamlined complexity—a notable achievement. Compared to previous gravitational wave models, we observe a significant acceleration in training time while maintaining acceptable performance metrics for ground-based detector signals and achieving equal or even superior metrics for space-based detector signals. Our experiments on synthetic data yield impressive results, with the model achieving an AUC score of 96.1% and a perfect recall rate of 100% on a dataset with a 1:3 class imbalance for ground-based detectors. For high signal-to-noise ratio signals, we achieve flawless precision and recall of 100% without losing precision on datasets with low-class ratios. Additionally, our approach reduces inference time by a factor of 1.88.

The most luminous AGN do not produce the majority of the detected stellar-mass black hole binary mergers in the local Universe

ABSTRACT Despite the increasing number of gravitational wave (GW) detections, the astrophysical origin of binary black hole (BBH) mergers remains elusive. A promising formation channel for BBHs is inside accretion discs around supermassive black holes, that power active galactic nuclei (AGN). In this paper, we test for the first time the spatial correlation between observed GW events and AGN. To this end, we assemble all sky catalogues with 1,412 (242) AGN with a bolometric luminosity greater than 1045.5erg s−1 ($10^{46}\, {\rm erg\, s}^{-1}$) with spectroscopic redshift of z ≤ 0.3 from the Milliquas catalogue, version 7.7b. These AGN are cross-matched with localization volumes of BBH mergers observed in the same redshift range by the LIGO and Virgo interferometers during their first three observing runs. We find that the fraction of the detected mergers originated in AGN brighter than $10^{45.5}\, {\rm erg\, s}^{-1}$ ($10^{46}\, {\rm erg\, s}^{-1}$) cannot be higher than 0.49 (0.17) at a 95 per cent credibility level. Our upper limits imply a limited BBH merger production efficiency of the brightest AGN, while most or all GW events may still come from lower luminosity ones. Alternatively, the AGN formation path for merging stellar-mass BBHs may be actually overall subdominant in the local Universe. To our knowledge, ours are the first observational constraints on the fractional contribution of the AGN channel to the observed BBH mergers.

Disappearing thermal X-ray emission as a tell-tale signature of merging massive black hole binaries

ABSTRACT The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localization, multiwavelength, time-dependent EM models are needed to identify the right host galaxy. We studied merging MBHBs embedded in a circumbinary disc (CBD) using high-resolution two-dimensional simulations, with a Γ-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the CBD. We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multiband EM signature: we find that the MBHB produces strong thermal X-ray emission until 1–2 d prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days, regardless of any post-merger effects such as GW recoil or mass-loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring.

Trajectories of Bright Stars and Shadows around Supermassive Black Holes as Tests of Gravity Theories

General relativity (GR), created more than a century ago, has been checked in various experimental and observational tests. At an early stage of its development, GR predictions were tested in problems where the gravitational field is weak and relativistic corrections can be considered as small perturbations of the Newtonian theory of gravity. However, in recent years due to the progress of new technologies it turned out to be possible to verify the predictions of GR in the limit of a strong gravitational field, as it was done to verify predictions about the profile of the X-ray line of iron Kalpha , estimates of the gravitational wave signal during the mergers of binary black holes and/or neutron stars and during the reconstruction of the shadows of black holes in Sgr A* and M87*. Groups of astronomers using the Keck and VLT (GRAVITY) telescopes confirmed the GR predictions for the redshift of the spectral lines of the S2 star near the passage of its pericenter (these predictions were done in the first post-Newtonian approximation). It is expected that in the near future, observations of bright stars using large telescopes VLT (GRAVITY), Keck, E-ELT and TMT will allow us to verify the predictions of GR in the strong gravitational field of supermassive black holes. Observations of bright stars in the vicinity of the Galactic Center and reconstructions of the shadows of black holes allow not only to verify the predictions of the GR, but also to obtain restrictions on alternative theories of gravity.

Evolutionary Origins of Binary Neutron Star Mergers: Effects of Common Envelope Efficiency and Metallicity

The formation histories of compact binary mergers, especially stellar-mass binary black hole mergers, have recently come under increased scrutiny and revision. We revisit the question of the dominant formation channel and efficiency of forming binary neutron star (BNS) mergers. We use the stellar and binary evolution code MESA and implement a detailed method for common envelope and mass transfer. We perform simulations for donor masses between 7 M ⊙ and 20 M ⊙ with a neutron star (NS) companion of 1.4 M ⊙ and 2.0 M ⊙ at two metallicities, using varying common envelope efficiencies and two different prescriptions to determine if the donor undergoes core collapse or electron capture, given their helium and carbon–oxygen cores. In contrast to the case of binary black hole mergers, for an NS companion of 1.4 M ⊙, all BNS mergers are formed following a common envelope phase. For an NS mass of 2.0 M ⊙, we identify a small subset of mergers following only stable mass transfer if the NS receives a natal kick sampled from a Maxwellian distribution with velocity dispersion σ = 265 km s−1. Regardless of the supernova prescription, we find more BNS mergers at subsolar metallicity compared to solar.

Astrophysical uncertainties in the gravitational-wave background from stellar-mass compact binary mergers

ABSTRACT We investigate the stochastic gravitational-wave background (SGWB) produced by merging binary black holes (BBHs) and binary neutron stars (BNSs) in the frequency ranges of Laser Interferometer Gravitational-Wave Observatory (LIGO)/Virgo/Kagra and Laser Interferometer Space Antenna (LISA). We develop three analytical models, which are calibrated to the measured local merger rates, and complement them with three population synthesis models based on the cosmic code. We discuss the uncertainties, focusing on the impact of the BBH mass distribution, the effect of the metallicity of the progenitor stars, and the time delay distribution between star formation and compact binary merger. We also explore the effect of uncertainties in binary stellar evolution on the background. For BBHs, our analytical models predict ΩGW in the range [4 × 10−10 to 1 × 10−9] (25 Hz) and [1 × 10−12 to 4 × 10−12] (3 mHz), and between [2 × 10−10 to 2 × 10−9] (25 Hz) and [7 × 10−13 to 7 × 10−12] (3 mHz) for our population synthesis models. This background is unlikely to be detected during the LIGO/Virgo/Kagra O4 run, but could be detectable with LISA. We predict about 10 BBH and no BNS mergers that could be individually detectable by LISA for a period of observation of 4 yr. Our study provides new insights into the population of compact binaries and the main sources of uncertainty in the astrophysical SGWB.

The Division between Weak and Strong Explosions from Failed Supernovae

Some massive stars likely fail to produce core-collapse supernovae, but these failed supernovae (FSNe) can generate an electromagnetic outburst prior to the disappearance of the star, as the mass lost to neutrinos during the stellar core collapse results in the formation and breakout of a second shock. We show that, when the mass lost to neutrinos is sufficiently small, there are two self-similar solutions that describe the propagation of a weak shock into a hydrodynamically expanding envelope that simultaneously yield accretion onto the black hole. The larger Mach number solution is unstable and yields the minimum Mach number that a shock must have to strengthen into the energy-conserving regime. Above a critical mass loss, there are no weak-shock solutions, implying that there are only strong explosions if the neutrino mass loss is above a critical value, and this value is a few percent of the mass of the star (and is physically achievable) for typical parameters. Our results imply that the fate of the explosion from an FSN—weak with little to no mass ejection or strong with the expulsion of the majority of the envelope—is a sensitive function of the stellar properties and the neutrino mass loss. We also show that there is a second type of self-similar solution for the shock that results in the settling of the gas near the compact object, which may be applicable to nonterminal stellar eruptions and the response of a gaseous disk to gravitational-wave induced mass loss from a binary black hole merger.

Things That Might Go Bump in the Night: Assessing Structure in the Binary Black Hole Mass Spectrum

Several features in the mass spectrum of merging binary black holes (BBHs) have been identified using data from the Third Gravitational Wave Transient Catalog (GWTC-3). These features are of particular interest as they may encode the uncertain mechanism of BBH formation. We assess if the features are statistically significant or the result of Poisson noise due to the finite number of observed events. We simulate catalogs of BBHs whose underlying distribution does not have the features of interest, apply the analysis previously performed on GWTC-3, and determine how often such features are spuriously found. We find that one of the features found in GWTC-3, the peak at ∼35 M ☉, cannot be explained by Poisson noise alone: peaks as significant occur in 1.7% of catalogs generated from a featureless population. This peak is therefore likely to be of astrophysical origin. The data is suggestive of an additional significant peak at ∼10 M ☉, though the exact location of this feature is not resolvable with current observations. Additional structure beyond a power law, such as the purported dip at ∼14 M ☉, can be explained by Poisson noise. We also provide a publicly available package, GWMockCat, that creates simulated catalogs of BBH events with correlated measurement uncertainty and selection effects according to user-specified underlying distributions and detector sensitivities.

Pulsational pair-instability supernovae in gravitational-wave and electromagnetic transients

ABSTRACT Current observations of binary black hole (BBH) merger events show support for a feature in the primary BH-mass distribution at $\sim \, 35 \ \mathrm{M}_{\odot }$, previously interpreted as a signature of pulsational pair-instability supernovae (PPISNe). Such supernovae are expected to map a wide range of pre-supernova carbon–oxygen (CO) core masses to a narrow range of BH masses, producing a peak in the BH mass distribution. However, recent numerical simulations place the mass location of this peak above $50 \ \mathrm{M}_{\odot }$. Motivated by uncertainties in the progenitor’s evolution and explosion mechanism, we explore how modifying the distribution of BH masses resulting from PPISN affects the populations of gravitational-wave (GW) and electromagnetic (EM) transients. To this end, we simulate populations of isolated BBH systems and combine them with cosmic star formation rates. Our results are the first cosmological BBH-merger predictions made using the binary_c rapid population synthesis framework. We find that our fiducial model does not match the observed GW peak. We can only explain the $35 \ \mathrm{M}_{\odot }$ peak with PPISNe by shifting the expected CO core-mass range for PPISN downwards by $\sim {}15 \ \mathrm{M}_{\odot }$. Apart from being in tension with state-of-the art stellar models, we also find that this is likely in tension with the observed rate of hydrogen-less super-luminous supernovae. Conversely, shifting the mass range upward, based on recent stellar models, leads to a predicted third peak in the BH mass function at $\sim {}64 \ \mathrm{M}_{\odot }$. Thus we conclude that the $\sim {}35 \ \mathrm{M}_{\odot }$ feature is unlikely to be related to PPISN.

GRMHD study of accreting massive black hole binaries in astrophysical environment: A review

We present an overview of recent numerical advances in the theoretical characterization of massive binary black hole (MBBH) mergers in astrophysical environments. These systems are among the loudest sources of gravitational waves (GWs) in the universe and particularly promising candidates for multimessenger astronomy. Coincident detection of GWs and electromagnetic (EM) signals from merging MBBHs is at the frontier of contemporary astrophysics. One major challenge in observational efforts searching for these systems is the scarcity of strong predictions for EM signals arising before, during, and after merger. Therefore, a great effort in theoretical work to-date has been to characterize EM counterparts emerging from MBBHs concurrently to the GW signal, aiming to determine distinctive observational features that will guide and assist EM observations. To produce sharp EM predictions of MBBH mergers it is key to model the binary inspiral down to coalescence in a full general relativistic fashion by solving Einstein’s field equations coupled with the magnetohydrodynamics equations that govern the evolution of the accreting plasma in strong-gravity. We review the general relativistic numerical investigations that have explored the astrophysical manifestations of MBBH mergers in different environments and focused on predicting potentially observable smoking-gun EM signatures that accompany the gravitational signal.

Precessing binary black holes as better dark sirens

ABSTRACT Gravitational waves (GWs) from binary black hole mergers provide unique opportunities for cosmological inference such as standard sirens. However, the accurate determination of the luminosity distance of the event is limited by the correlation between the distance and the angle between the binary’s orbital angular momentum and the observer’s line of sight. In the letter, we investigate the effect of precession on the distance estimation of binary black hole events for the third-generation (3G) GW detectors. We find that the precession can enhance the precision of distance inference by one order of magnitude compared to the scenario where precession is absent. The constraint on the host galaxies can be improved due to the improved distance measurement, therefore the Hubble constant can be measured with higher precision and accuracy. These findings underscore the noteworthy impact of precession on the precision of distance estimation for 3G ground-based GW detectors, which can serve as highly accurate probes of the Universe.

Revisiting the properties of GW190814 and its formation history

ABSTRACT GW190814 was reported during LIGO’s and Virgo’s third observing run with the most asymmetric component masses (an ∼23 M⊙ black hole and an ∼2.6 M⊙ compact object). Under the assumption that this event is a binary black hole (BBH) merger formed through the isolated binary evolution channel, we reanalyse the publicly released data of GW190814 with the modified astrophysical priors on the effective spin χeff, and further explore its formation history using detailed binary modelling. We show that GW190814 is likely to have been formed through the classical common envelope channel. Our findings show that the properties inferred using the modified astrophysical priors are consistent with those inferred by the uniform priors. With the newly inferred properties of GW190814, we perform detailed binary evolution of the immediate progenitor of the BBH (namely a close binary system composed of a BH and a helium star) in a large parameter space, taking into account mass-loss, internal differential rotation, supernova kicks, and tidal interactions between the helium star and the BH companion. Our findings show that GW190814-like events could be formed in limited initial conditions just after the common envelope phase: an ∼23 M⊙ BH and a helium star of MZamsHe ∼ 8.5 M⊙ at solar metallicity (∼ 7.5 M⊙ at 10 per cent solar metallicity) with an initial orbital period at around 1.0 d. Additionally, the inferred low spin of the secondary indicates that the required metallicity for reproducing GW190814-like events should not be too low (e.g. Z ≳ 0.1 Z⊙).

On the Likely Dynamical Origin of GW191109 and Binary Black Hole Mergers with Negative Effective Spin

With the growing number of binary black hole (BBH) mergers detected by LIGO/Virgo/KAGRA, several systems have become difficult to explain via isolated binary evolution, having components in the pair-instability mass gap, high orbital eccentricities, and/or spin–orbit misalignment. Here we focus on GW191109_010717, a BBH merger with component masses of and M ⊙ and an effective spin of , which could imply a spin–orbit misalignment of more than π/2 rad for at least one of its components. Besides its component masses being in the pair-instability mass gap, we show that isolated binary evolution is unlikely to reproduce the proposed spin–orbit misalignment of GW191109 with high confidence. On the other hand, we demonstrate that BBHs dynamically assembled in dense star clusters would naturally reproduce the spin–orbit misalignment and masses of GW191109 and the rates of GW191109-like events if at least one of the components were to be a second-generation BH. Finally, we generalize our results to all events with a measured negative effective spin, arguing that GW200225 also has a likely dynamical origin.

Massive binary black holes from Population II and III stars

ABSTRACT Population III stars, born from the primordial gas in the Universe, lose a negligible fraction of their mass via stellar winds and possibly follow a top-heavy mass function. Hence, they have often been regarded as the ideal progenitors of massive black holes (BHs), even above the pair instability mass gap. Here, we evolve a large set of Population III binary stars (metallicity Z = 10−11) with our population-synthesis code sevn, and compare them with Population II binary stars (Z = 10−4). In our models, the lower edge of the pair-instability mass gap corresponds to a BH mass of ≈86 (≈91) M⊙ for single Population III (II) stars. Overall, we find only mild differences between the properties of binary BHs (BBHs) born from Population III and II stars, especially if we adopt the same initial mass function and initial orbital properties. Most BBH mergers born from Population III and II stars have primary BH mass below the pair-instability gap, and the maximum secondary BH mass is <50 M⊙. Only up to ≈3.3 per cent (≈0.09 per cent) BBH mergers from Population III (II) progenitors have primary mass above the gap. Unlike metal-rich binary stars, the main formation channel of BBH mergers from Population III and II stars involves only stable mass transfer episodes in our fiducial model.

The role of stellar expansion on the formation of gravitational wave sources

ABSTRACT Massive stars are the progenitors of black holes and neutron stars, the mergers of which can be detected with gravitational waves (GW). The expansion of massive stars is one of the key factors affecting their evolution in close binary systems, but it remains subject to large uncertainties in stellar astrophysics. For population studies and predictions of GW sources, the stellar expansion is often simulated with the analytic formulae from Hurley et al. (2000). These formulae need to be extrapolated and are often considered outdated. In this work, we present five different prescriptions developed from 1D stellar models to constrain the maximum expansion of massive stars. We adopt these prescriptions to investigate how stellar expansion affects mass transfer interactions and in turn the formation of GW sources. We show that limiting radial expansion with updated 1D stellar models, when compared to the use of Hurley et al. (2000) radial expansion formulae, does not significantly affect GW source properties (rates and masses). This is because most mass transfer events leading to GW sources are initialized in our models before the donor star reaches its maximum expansion. The only significant difference was found for the mass distribution of massive binary black hole mergers (Mtot > 50 M⊙) formed from stars that may evolve beyond the Humphreys–Davidson limit, whose radial expansion is the most uncertain. We conclude that understanding the expansion of massive stars and the origin of the Humphrey–Davidson limit is a key factor for the study of GW sources.

Understanding Binary Systems—a Comparison between COSMIC and MESA

We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, as they are promising candidates for the progenitor systems behind energetic, longer lasting, radio bright GRB jets. These systems also may be key contributors to the rates of binary black hole mergers throughout our Universe.

Symmetry Breaking in Merging Binary Black Holes from Young Massive Clusters and Isolated Binaries

Properties of the to-date observed binary black hole (BBH) merger events suggest a preference toward spin–orbit aligned mergers. Naturally, this has caused widespread interest and speculations regarding implications on various merger formation channels. Here we show that (i) not only the BBH merger population from isolated binaries but also (ii) BBH population formed in young massive clusters (YMCs) would possess an asymmetry in favor of aligned mergers, in the distribution of the events’ effective spin parameter (χ eff). In our analysis, we utilize BBH merger outcomes from state-of-the-art N-body evolutionary models of YMCs and isolated binary population synthesis. We incorporate, for the first time in such an analysis, misalignments due to both natal kicks and dynamical encounters. The YMC χ eff distribution has a mean (an antialigned merger fraction) of 〈χ eff〉 ≤ 0.04 (f X − ≈ 40%), which is smaller (larger) than but consistent with the observed asymmetry of 〈χ eff〉 ≈ 0.06 (f X − ≈ 28%) as obtained from the population analysis by the LIGO–Virgo–KAGRA collaboration. In contrast, isolated binaries alone tend to produce a much stronger asymmetry; for the tested physical models, 〈χ eff〉 ≈ 0.25 and f X − ≲ 7%. Although the YMC χ eff distribution is more similar to the observed counterpart, none of the channels correctly reproduce the observed distribution. Our results suggest that further extensive model explorations for both isolated binary and dynamical channels as well as better observational constraints are necessary to understand the physics of “the symmetry breaking” of the BBH merger population.

Event Rate of Strongly Lensed Gravitational Waves of Stellar Binary Black Hole Mergers Produced by Dynamical Interactions

Gravitational waves emitted from stellar binary black hole (sBBH) mergers can be gravitationally lensed by intervening galaxies and detected by future ground-based detectors. A great amount of effort has been put into the estimation of the detection rate of lensed sBBH originating from the evolution of massive binary stars (EMBS channel). However, sBBHs produced by the dynamical interaction in dense clusters (dynamical channel) may also be dominant in our universe and their intrinsic distribution of physical properties can be significantly different from those produced by massive stars, especially mass and redshift distribution. In this paper, we investigate the event rate of lensed sBBHs produced via the dynamical channel by Monte Carlo simulations and the number is yr−1 for the Einstein telescope and yr−1 for Cosmic Explorer, of which the median is about ∼2 times the rate of sBBHs originating from the EMBS channel (calibrated by the local merger rate density estimated for the dynamical and the EMBS channel, i.e., and , respectively). Therefore, one may constrain the fraction of both the EMBS and dynamical channels through the comparison of the predicted and observed number of lensed sBBH events statistically.