Binary Black Hole Mergers Research Articles

Development of convective envelopes in massive stars. Implications for gravitational wave sources

The structure of stellar envelopes strongly influences the course and outcome of binary mass transfer, in particular of common-envelope (CE) evolution. Convective envelopes can most easily be ejected during CE events, leading to short-period binaries and, potentially, gravitational-wave (GW) sources. Conversely, radiative envelopes are thought to lead to CE mergers and Thorne-\.Zytkow objects (T\.ZOs) or quasi-stars (QSs). Rapid binary models based on Hurley et al. (2000) often assume that any CE event with a Hertzsprung gap donor results in a CE merger, in tension with the literature. We improve on this assumption with a more self-consistent criterion based on the presence of a convective envelope. Using 1D stellar models ( MESA ), we systematically investigated the development of convective envelopes in massive stars. We provided fitting formulae for rapid binary codes and implemented them into the StarTrack population synthesis code to refine the CE treatment and examined the impact on GW sources, T\.ZOs, and QSs. We show that convective envelopes in massive stars are highly sensitive to the treatment of superadiabacity and the mixing length. Our revised CE model significantly reduces (by a factor sim 20) the predicted merger rate of binary black hole (BH-BH) mergers with total masses between sim 20 and 50 $M_ This leads to a bimodal mass distribution with a strong metallicity dependence. We also predict that the current T\.ZO--QS formation rate in the Galaxy (up to sim $ yr$^ $), combined with their predicted lifetimes, makes their detection unlikely. Our study strongly suggests that the role of CE evolution in the formation of BH-BH mergers has been considerably overestimated for BH-BH mergers with $M_ tot $\,geq \,20 $M_ We highlight that any prediction from the CE channel for massive BH-BH mergers (>50 $M_ heavily hinges on our limited understanding of stellar structure and mass loss close to the Eddington limit.

A physically modelled selection function for compact binary mergers in the LIGO-Virgo O3 run and beyond

Abstract Despite the observation of nearly 100 compact binary coalescence (CBC) events up to the end of the Advanced gravitational-wave (GW) detectors' third observing run (O3), there remain fundamental open questions regarding their astrophysical formation mechanisms and environments. Population analysis should yield insights into these questions, but requires careful control of uncertainties and biases. GW observations have a strong selection bias: this is due first to the dependence of the signal amplitude on the source's (intrinsic and extrinsic) parameters, and second to the complicated nature of detector noise and of current detection methods. In this work, we introduce a new physically-motivated model of the sensitivity of GW searches for CBC events, aimed at enhancing the accuracy and efficiency of population reconstructions. In contrast to current methods which rely on re-weighting simulated signals (injections) via importance sampling, we model the probability of detection of binary black hole (BBH) mergers as a smooth, analytic function of source masses, orbit-aligned spins, and distance, fitted to accurately match injection results. The estimate can thus be used for population models whose signal distribution over parameter space differs significantly from the injection distribution. Our method has already been used in population studies such as reconstructing the BBH merger rate dependence on redshift.

PyMerger: Detecting Binary Black Hole Mergers from the Einstein Telescope Using Deep Learning

We present PyMerger, a Python tool for detecting binary black hole (BBH) mergers from the Einstein Telescope (ET), based on a deep residual neural network (ResNet) model. ResNet was trained on data combined from all three proposed subdetectors of ET (TSDCD) to detect BBH mergers. Five different lower-frequency cutoffs (F low)—5 Hz, 10 Hz, 15 Hz, 20 Hz, and 30 Hz—with the match-filter signal-to-noise ratio (MSNR) ranges 4–5, 5–6, 6–7, 7–8, and >8 were employed in the data simulation. Compared to previous work that utilized data from a single subdetector, the detection accuracy from TSDCD has shown substantial improvements, increasing from 60%, 60.5%, 84.5%, 94.5%, and 98.5% to 78.5%, 84%, 99.5%, 100%, and 100% for sources with MSNRs of 4–5, 5–6, 6–7, 7–8, and >8, respectively. The ResNet model is evaluated on the first ET mock data challenge (ET-MDC1) data set, where the model demonstrates strong performance in detecting BBH mergers, identifying 5566 out of 6578 BBH events, with optimal SNRs starting from 1.2 and a minimum and maximum D L of 0.5 Gpc and 148.95 Gpc, respectively. Despite being trained only on BBH mergers without overlapping sources, the model achieves high BBH detection rates. Notably, even though the model was not trained on binary neutron star (BNS) and black hole-neutron star (BHNS) mergers, it successfully detected 11,477 BNS and 323BHNS mergers in ET-MDC1, with optimal SNRs starting from 0.2 and 1, respectively, indicating its potential for broader applicability.

Detection Rate of Galaxy Cluster-lensed Stellar Binary Black Hole Mergers by the Third-generation Gravitational-wave Detectors

Gravitational waves (GWs) from stellar binary black hole (sBBH) mergers can be strongly gravitational lensed by intervening galaxies/galaxy clusters. Only a few works have investigated cluster-lensed sBBH mergers by adopting oversimplified models, while galaxy-lensed ones have been intensively studied. In this paper, we estimate the detection rate of cluster-lensed sBBH mergers with third-generation GW detectors and its dependence on the lens models. We adopt detailed modeling of galaxy cluster lenses by using the mock clusters in the Synthetic Sky Catalog for Dark Energy Science with LSST (cosmoDC2) and/or approximations of the pseudo-Jaffe profile or an eccentric Navarro–Frenk–White dark matter halo plus a bright central galaxy with singular isothermal sphere profile. Considering that the formation of sBBH mergers is dominated by the evolution of the massive binary star (EMBS) channel, we find that the detection rate of cluster-lensed sBBHs is ∼5–84 yr−1, depending on the adopted lens model and uncertainty in the merger rate density, which is about ∼13−2.0+28 yr−1 if adopting relatively more realistic galaxy clusters with central main and member galaxies in the cosmoDC2 catalog, close to the estimated detection rate of sBBH mergers lensed by galaxies. In addition, we also consider the case in which the production of sBBH mergers is dominated by the dynamical interactions in dense stellar systems. We find that the detection rate of cluster-lensed sBBHs from the dynamical channel is about 1.5 times larger than that from the EMBS channel, and the redshift distribution of the former peaks at a higher redshift (∼3) compared with that from the latter (∼2).

Exploring Field-evolution and Dynamical-capture Coalescing Binary Black Holes in GWTC-3

We investigate formation channels for merging binary black holes (BBHs) in GWTC-3, with a dedicated semiparametric population model. The model first describes or excludes a high-spin (with magnitudes of ∼0.7) and high-mass (ranging in ∼20–80M ⊙) subpopulation, which was identified by previous works and can be interpreted as hierarchical mergers. We find that the rest of BBH population can be categorized into two subpopulations with different mass and mass-ratio distributions, as indicated by a Bayes factor of lnB=1.8 . One subpopulation, characterized by nearly aligned spins and consistent with isolated field formation, likely dominates the 10-solar-mass peak in the primary-mass function. The other subpopulation, with isotropic spins and consistent with the dynamical channels, shows a stronger preference for symmetric pairing, and mainly contributes to the 35-solar-mass peak in the primary-mass function. Note that the Bayes factor is not high enough with the currently available data, so that the case of a single population is still acceptable. Additionally, we compare the mass distributions between merging BBHs and black holes (BHs) in high-mass X-ray binaries (HMXBs). We find that the primary mass of the aligned subpopulation is slightly lighter than those of the HMXB BHs, while the isotropic subpopulation is consistent with the HMXB BHs, if the power-law index of its mass function is shifted by 2, as indicated by the dynamical formation channels. However, the spin magnitudes of both subpopulations are significantly smaller than those of the HMXB BHs.

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

Context. 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. Aims. 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 102–105 M⊙. Methods. 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. Results. Our findings indicate that ET could detect a binary black hole (BBH) merger rate of around 2 × 105 Gpc−3 yr−1 for BH masses ranging from 10 to 200 M⊙, with around 100 Gpc−3 yr−1 of this rate specifically attributed to BHs in the 100–200 M⊙ 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 ≈ 2. Conclusions. 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.

Mass Function of Stellar Black Holes as Revealed by the LIGO–Virgo–KAGRA Observations

Ninety gravitational-wave events have been detected by the LIGO–Virgo–KAGRA network and are released in the Gravitational-Wave Transient Catalog. Among these events, 83 cases are definitely binary black hole mergers, since the masses of all the objects involved significantly exceed the upper limit of neutron stars. The black holes in these merger events naturally form two interesting samples, a premerger sample that includes all the black holes before the mergers and a postmerger sample that consists of the black holes generated during the merging processes. The former represents black holes that once existed in the Universe, while the latter represents newly born black holes. Here we present a statistical analysis of these two samples. The nonparametric τ statistic method is adopted to correct for the observational selection effect. The Lynden-Bell C − method is further applied to derive the mass distribution and density function of black holes. It is found that the mass distribution can be expressed as a broken power-law function. More interestingly, the power-law index in the high-mass region is comparable for the two samples. The number density of black holes is found to depend on redshift as ρ(z) ∝ z −2.06—z −2.12 based on the two samples. The implications of these findings on the origin of black holes are discussed.

Black hole-black hole mergers with and without an electromagnetic counterpart: A model for stable tertiary mass transfer in hierarchical triple systems

Triple stars are prevalent within the population of observed stars. Their evolution compared to binary systems is notably more complex and is influenced by unique dynamical, tidal, and mass transfer processes inherent in higher order multiples. Understanding these phenomena is essential for comprehensive insight into multistar evolution and the formation of energetic transients, including gravitational wave (GW) mergers. Our study aims to probe the evolution of triple star systems when the tertiary component fills its Roche lobe and transfers mass to the inner binary. Specifically, we focus on the impact of tertiary mass transfer on the evolution of the inner orbit and investigate whether it could lead to the formation of GW sources with distinct properties. To achieve this, we developed an analytical model that describes the evolution of the inner and outer orbits of hierarchical triples undergoing stable mass transfer from the tertiary component. We have publicly released this model as a python package on Zenodo. Utilising population synthesis simulations, we investigated triples with a Roche-lobe filling tertiary star and an inner binary black hole (BBH). These systems stem from inner binaries experiencing chemically homogeneous evolution (CHE). Our analysis encompasses two distinct populations with metallicities of $Z=0.005$ and $Z=0.0005$, focusing on primary components in the inner binary with initial masses ranging from 20 to $100\ M odot $ and inner and outer orbital separations of up to $40\ R odot $ and $10^5\ R odot $, respectively, targeting the parameter space where chemically homogeneous evolution is anticipated. Our results indicate that for the systems we studied, the mass transfer phase predominantly leads to orbital shrinkage of the inner binary and evolution towards non-zero eccentricities and is accompanied by an expansion of the outer orbit. In the systems where the inner binary components evolve in a chemically homogeneous manner, 9.5<!PCT!> result in mass transfer from the tertiary onto an inner BBH. Within this subset, we predict a high formation efficiency of GW mergers ranging from 85.1<!PCT!> to 100<!PCT!> at $Z=0.005$ and 100<!PCT!> at $Z=0.0005$ with short delay times, partly attributable to the mass transfer phase. Owing to the rarity of triples with a CHE inner binary in the stellar population, we project local merger rates in the range of 0.69 to 1.74 $ Gpc^ $. Of the prospected BBH mergers that enter the LISA and aLIGO frequency band due to GW emission, a fraction is still accreting gas from the tertiary star. This could produce a strong electromagnetic (EM) counterpart to the GW source and maintain high eccentricities as the system enters the frequency range detectable by GW detectors. The occurrence of EM signals accompanying mergers varies significantly depending on model assumptions, with fractions ranging from less than 0.03<!PCT!> to as high as 46.8<!PCT!> of all mergers if the formation of a circumbinary disc is allowed.

Detecting a gravitational wave background from inflation with null energy condition violation: prospects for Taiji

The null energy condition (NEC) is a fundamental principle in general relativity, and its violation could leave discernible signatures in gravitational waves (GWs). A violation of the NEC during the primordial era would imprint a blue-tilted spectrum on the stochastic gravitational wave background (SGWB) at nanohertz frequencies, potentially accounting for the recently detected signal by pulsar timing arrays. Remarkably, models of NEC violation during inflation also predict a nearly scale-invariant GW spectrum in the millihertz frequency range, which could be detectable by upcoming space-based GW detectors such as Taiji. The observation of this distinctive spectrum would provide compelling evidence for new physics beyond the standard cosmological paradigm. In this study, we explore Taiji’s ability to detect an SGWB arising from NEC violation during inflation, considering various foregrounds and noise sources, including an extragalactic foreground from binary black hole mergers throughout the Universe, a galactic foreground from white dwarf binaries, and the intrinsic noise of the Taiji detector. Employing comprehensive Bayesian parameter estimation techniques to analyze simulated Taiji data, we demonstrate a remarkable improvement in precision of three orders of magnitude compared to the NANOGrav 15-year data set for measuring the tensor power spectrum amplitude, PT,2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$P_{T,2}$$\\end{document}, during the second inflationary stage. This substantial enhancement in measurement capabilities underscores Taiji’s potential as a powerful probe for investigating the NEC violation in the early Universe.

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 ($H_0$) from the standard dark siren method using a sample of five well-covered gravitational wave (GW) alerts reported during the first part of the fourth observing run of the Laser Interferometer Gravitational-Wave Observatory (LIGO), the Virgo and Kamioka Gravitational Wave Detector (KAGRA) collaborations (LVK) and with three updated standard dark sirens from third observation run in combination with 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 $H_0$ posterior for 5 new standard sirens with other 10 previous events (using the most recent available data for the five novel events and updated three previous posteriors from O3), finding $H_0 = 70.4^{+13.6}_{-11.7}~{\rm km~s^{-1}~Mpc^{-1}}$ (68 per cent confidence interval) with the catalogue method only. This result represents an improvement of $\sim 23~{{\ \rm per\ cent}}$ comparing the new 15 dark siren constraints with the previous 10 dark siren constraints and a reduction in uncertainty of $\sim 40~{{\ \rm per\ cent}}$ 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 $H_0$ of $68.0^{+4.4}_{-3.8}~{\rm km~s^{-1}~Mpc^{-1}}$, a $\sim 6~{{\ \rm per\ cent}}$ precision from standard sirens, reducing the previous constraint uncertainty by $\sim 10~{{\ \rm per\ cent}}$.

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