Compact Binary Mergers Research Articles

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.

Compact Binary Merger Rate with Modified Gravity in Dark Matter Spikes

In this study, we investigate the impact of modified gravity (MG) on the merger rate of compact binaries within dark matter spikes surrounding supermassive black holes (SMBHs). Specifically, we calculate the binary merger rates involving primordial black holes (PBHs) and/or neutron stars (NSs) in Hu–Sawicki f(R) gravity and the normal branch of Dvali–Gabadadze–Porrati (nDGP) gravity, with three SMBH mass functions: Benson, Vika, and Shankar. The results show consistently higher merger rates predicted for PBH–PBH and PBH–NS binaries in these gravity models compared to general relativity (GR), in particular at lower SMBH masses and for steeper dark matter spike density profiles. The predicted merger rates are compared to the LIGO–Virgo–KAGRA observations in constraining the parameters of the theory. In particular, we find steeper dark matter spike density profiles in the MG scenarios compared to GR. When compared to current observational constraints on PBH abundance, the mass ranges allowed by Hu–Sawicki f(R) models are found to be wider than those allowed by nDGP models, for given merger rates. The results are highly dependent on the choice of SMBH mass function, with the Vika and Shankar mass functions predicting lower abundances. The considerable sensitivity of the results to the assumed gravity scenario and SMBH mass function demonstrates the necessity of incorporating the corresponding theoretical uncertainties when making relatively robust predictions on compact binary merger rates and, as a result, on PBH properties.

Multispectral Sirens: Gravitational-wave Cosmology with (Multi-) Subpopulations of Binary Black Holes

The cosmic expansion rate can be directly measured with gravitational-wave (GW) data of the compact binary mergers by jointly constraining the mass function of the population and the cosmological model via the so-called spectral sirens. Such a method relies on the features in the mass functions, which may originate from some individual subpopulations and hence become blurred/indistinct due to the superposition of different subpopulations. In this work, we propose a novel approach to constrain the cosmic expansion rate with subpopulations of GW events, named multispectral sirens. The advantage of the multispectral sirens compared to the traditional spectral sirens is demonstrated by the simulation with the mock data. The application of this approach to the GWTC-3 data yields H0=73.3−25.6+29.9Mpc−1kms−1 (median and symmetric 68.3% credible interval), which is about 19% tighter than the result inferred with the traditional spectral sirens utilizing a powerlaw+peak mass function. The incorporation of the bright standard siren GW170817 with a uniform prior in [10, 200] (log-uniform prior in [20,140]) Mpc−1 km s−1 gives H0=71.1−7.5+15.0(70.3−7.1+12.9)Mpc−1kms−1 (68.3% confidence level), corresponding to an improvement of ∼26% (23%) with respect to the measurement from sole GW170817.

Correction to: Rapid source classification and distance estimation for compact binary mergers with PyCBC live

Correction to: Rapid source classification and distance estimation for compact binary mergers with PyCBC live

Compact object populations over cosmic time II. Compact object merger rates and masses over redshift from varying initial conditions

ABSTRACT We perform a first study of the impact of varying two components of the initial conditions in binary population synthesis of compact binary mergers – the initial mass function, which is made metallicity- and star formation rate-dependent, and the orbital parameter (orbital period, mass ratio, and eccentricity) distributions, which are assumed to be correlated – within a larger grid of initial condition models also including alternatives for the primary mass-dependent binary fraction and the metallicity-specific cosmic star formation history. We generate the initial populations with the sampling code bossa and evolve them with the rapid population synthesis code compas. We find strong suggestions that the main role of initial conditions models is to set the relative weights of key features defined by the evolution models. In the two models we compare, black hole–black hole (BHBH) mergers are the most strongly affected, which we connect to a shift from the common envelope to the stable Roche lobe overflow formation channels with decreasing redshift. We also characterize variations in the black hole–neutron star (BHNS) and neutron star–neutron star (NSNS) final parameter distributions. We obtain the merger rate evolution for BHBH, BHNS, and NSNS mergers up to $z=10$, and find a variation by a factor of $\sim 50\textnormal {--}60$ in the local BHBH and BHNS merger rates, suggesting a more important contribution from initial conditions than previously thought, and calling for a complete exploration of the initial conditions model permutations.

Unveiling the Progenitors of a Population of Likely Peculiar Gamma-Ray Bursts

Traditionally, gamma-ray bursts (GRBs) are classified as long and short GRBs, with T 90 = 2 s being the threshold duration. Generally, long-duration GRBs (LGRBs; T 90 > 2 s) are associated with the collapse of massive stars, and short-duration (SGRBs; T 90 < 2 s) are associated with compact binary mergers involving at least one neutron star. However, the existence of a population of so-called “peculiar GRBs”—i.e., LGRBs originating from mergers or long Type I GRBs, and SGRBs originating from collapsars or short Type II GRBs—has challenged the traditional paradigm of GRB classification. Finding more peculiar GRBs may help to give more insight into this issue. In this work, we analyze the properties of machine-learning-identified long Type I GRB and short Type II GRB candidates, long GRBs-I and short GRBs-II (the so-called “peculiar GRBs”). We find that long GRBs-I almost always exhibit properties similar to Type I GRBs, which suggests that mergers may indeed produce GRBs with T 90 > 2 s. Furthermore, according to the probability given by the redshift distribution, short GRBs-II almost exhibit properties similar to Type II GRBs. This suggests that the populations of short Type II GRBs are not scarce and that they are hidden in a large number of samples without redshifts, which is unfavorable for the interpretation that the jet progression leads to a missed main emission.

Diversity in Fermi/GBM Gamma-Ray Bursts: New Insights from Machine Learning

Classification of gamma-ray bursts (GRBs) has been a long-standing puzzle in high-energy astrophysics. Recent observations challenge the traditional short versus long viewpoint, where long GRBs are thought to originate from the collapse of massive stars and short GRBs from compact binary mergers. Machine learning (ML) algorithms have been instrumental in addressing this problem, revealing five distinct GRB groups within the Swift Burst Alert Telescope (BAT) light-curve data, two of which are associated with kilonovae (KNe). In this work, we extend our analysis to the Fermi Gamma-ray Burst Monitor catalog and identify five clusters using unsupervised ML techniques, consistent with the Swift/BAT results. These five clusters are well separated in the fluence-duration plane, hinting at a potential link between fluence, duration, and complexities (or structures) in the light curves of GRBs. Further, we confirm two distinct classes of KN-associated GRBs. The presence of GRB 170817A in one of the two KN-associated clusters lends evidence to the hypothesis that this class of GRBs could potentially be produced by binary neutron star mergers. The second KN-associated GRB cluster could potentially originate from neutron star–black hole mergers. Future multimessenger observations of compact binaries in gravitational waves and electromagnetic waves can be paramount in understanding these clusters better.

Inferring Binary Properties from Gravitational-Wave Signals

This review provides a conceptual and technical survey of methods for parameter estimation of gravitational-wave signals in ground-based interferometers such as Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo. We introduce the framework of Bayesian inference and provide an overview of models for the generation and detection of gravitational waves from compact binary mergers, focusing on the essential features that are observable in the signals. Within the traditional likelihood-based paradigm, we describe various approaches for enhancing the efficiency and robustness of parameter inference. This includes techniques for accelerating likelihood evaluations, such as heterodyne/relative binning, reduced-order quadrature, multibanding, and interpolation. We also cover methods to simplify the analysis to improve convergence, via reparameterization, importance sampling, and marginalization. We end with a discussion of recent developments in the application of likelihood-free (simulation-based) inference methods to gravitational-wave data analysis.

Using equation of state constraints to classify low-mass compact binary mergers

Using equation of state constraints to classify low-mass compact binary mergers

New results on the gamma-ray burst variability--luminosity relation

At the dawn of the gamma--ray burst (GRB) afterglow era, a Cepheid-like correlation was discovered between the time variability $V$ and the isotropic-equivalent peak luminosity $L_ iso $ of the prompt emission of about a dozen long GRBs with measured redshift available at that time. Soon afterwards, the correlation was confirmed in a sample of about 30 GRBs, even though it was affected by significant scatter. Unlike the minimum variability timescale (MVT), $V$ measures the relative power of short-to-intermediate timescales. We aim to test the correlation using about 200 long GRBs with spectroscopically measured redshift, detected by Swift Fermi and Konus/ WIND for which both observables can be accurately estimated. The variability for all selected GRBs was calculated according to the original definition using the 64 ms background-subtracted light curves of Swift /BAT ( Fermi /GBM) in the 15--150 (8--900) keV energy passband. Peak luminosities were either taken from the literature or derived from modelling broad-band spectra acquired with either Konus/ WIND or Fermi /GBM. The statistical significance of the correlation has weakened to $ 2$&lt;!PCT!&gt;, mostly due to the appearance of a number of smooth and luminous GRBs that are characterised by a relatively small $V$. At odds with most long GRBs, three out of four long-duration merger candidates have high $V$ and low iso The luminosity is more tightly connected with shortest timescales measured by MVT than the short to intermediate timescales measured by $V$. We discuss the implications for internal dissipation models and the role of the $e^ pm $ photosphere. We identified a few smooth GRBs with a single broad pulse and low $V$ that might have an external shock origin, in contrast with most GRBs. The combination of high variability $(V 0.1)$, low luminosity iso $ erg s$^ $, and short MVT ($ 0.1$ s) could be a good indicator for a compact binary merger origin.

Gauss-Bonnet Cosmology: large-temperature behaviour and bounds from Gravitational Waves

We provide a transparent discussion of the high temperature asymptotic behaviour of Cosmology in a dilaton-Einstein-Gauss-Bonnet (dEGB) scenario of modified gravity with vanishing scalar potential. In particular, we show that it has a clear interpretation in terms of only three attractors (stable critical points) of a set of autonomous differential equations: w = -1/3, w = 1 and 1 < w < 7/3, where w ≡ p/ρ is the equation of state, defined as the ratio of the total pressure and the total energy density. All the possible different high-temperature evolution histories of the model are exhausted by only eight paths in the flow of the set of the autonomous differential equations. Our discussion clearly explains why five out of them are characterized by a swift transition of the system toward the attractor, while the remaining three show a more convoluted evolution, where the system follows a meta-stable equation of state at intermediate temperatures before eventually jumping to the real attractor at higher temperatures. Compared to standard Cosmology, the regions of the dEGB parameter space with w = -1/3 show a strong enhancement of the expected Gravitational Wave stochastic background produced by the primordial plasma of relativistic particles of the Standard Model. This is due to the very peculiar fact that dEGB allows to have an epoch when the energy density ρ rad of the relativistic plasma dominates the energy of the Universe while at the same time the rate of dilution with T of the total energy density is slower than what usually expected during radiation dominance. This allows to use the bound from Big Bang Nucleosynthesis (BBN) to put in dEGB a constraint T RH ≲ (108 – 109) GeV on the reheating temperature of the Universe T RH. Such BBN bound is complementary to late-time constraints from compact binary mergers.

The Pristine Inner Galaxy Survey (PIGS)

The most metal-poor stars provide valuable insights into the early chemical enrichment history of a system, carrying the chemical imprints of the first generations of supernovae. The most metal-poor region of the Sagittarius dwarf galaxy remains inadequately observed and characterised. To date, only ∼4 stars with [Fe/H] &lt; −2.0 have been chemically analysed with high-resolution spectroscopy. In this study, we present the most extensive chemical abundance analysis of 12 low-metallicity stars with metallicities down to [Fe/H] = −3.26 and located in the main body of Sagittarius. These targets, selected from the Pristine Inner Galaxy Survey, were observed using the MIKE high-resolution spectrograph at the Magellan-Clay telescope, which allowed us to measure up to 17 chemical species. The chemical composition of these stars reflects the imprint of a variety of type II supernovae (SNe II). A combination of low- to intermediate-mass high-energy SNe and hypernovae (∼10 − 70 M⊙) is required to account for the abundance patterns of the lighter elements up to the Fe-peak. The trend of the heavy elements suggests the involvement of compact binary merger events and fast-rotating (up to ∼300 km s−1) intermediate-mass to massive metal-poor stars (∼25 − 120 M⊙) that are the sources of rapid and slow processes, respectively. Additionally, asymptotic giant branch stars contribute to a wide dispersion of [Ba/Mg] and [Ba/Eu]. The absence of an α−knee in our data indicates that type Ia supernovae did not contribute in the very metal-poor region ([Fe/H] ≤ −2.0). However, they might have started to pollute the interstellar medium at [Fe/H] &gt; −2.0, given the relatively low [Co/Fe] in this metallicity region.

Impact of Higher Harmonics of Gravitational Radiation on the Population Inference of Binary Black Holes

Templates modeling just the dominant mode of gravitational radiation are generally sufficient for the unbiased parameter inference of near-equal-mass compact binary mergers. However, neglecting the subdominant modes can bias the inference if the binary is significantly asymmetric, very massive, or has misaligned spins. In this work, we explore if neglecting these subdominant modes in the parameter estimation of nonspinning binary black hole mergers can bias the inference of their population-level properties such as mass and merger redshift distributions. Assuming the design sensitivity of the advanced LIGO-Virgo detector network, we find that neglecting subdominant modes will not cause a significant bias in the population inference, although including them will provide more precise estimates. This is primarily because asymmetric binaries are expected to be rarer in our detected sample, due to their intrinsic rareness and the observational selection effects. The increased precision in the measurement of the maximum black hole mass can help in better constraining the upper mass gap in the mass spectrum.

GRB 211211A: The Case for an Engine-powered over r-process-powered Blue Kilonova

The recent gamma-ray burst (GRB) GRB 211211A provides the earliest (∼5 hr) data of a kilonova (KN) event, displaying bright (∼1042 erg s−1) and blue early emission. Previously, this KN was explained using simplistic multicomponent fitting methods. Here, in order to understand the physical origin of the KN emission in GRB 211211A, we employ an analytic multizone model for r-process-powered KNe. We find that r-process-powered KN models alone cannot explain the fast temporal evolution and the spectral energy distribution (SED) of the observed emission. Specifically, (i) r-process models require high ejecta mass to match early luminosity, which overpredicts late-time emission, while (ii) red KN models that reproduce late emission underpredict early luminosity. We propose an alternative scenario involving early contributions from the GRB central engine via a late low-power jet, consistent with plateau emission in short GRBs and GeV emission detected by Fermi-LAT at ∼104 s after GRB 211211A. Such late central engine activity, with an energy budget of ∼a few percent of that of the prompt jet, combined with a single red KN ejecta component, can naturally explain the light curve and SED of the observed emission, with the late-jet–ejecta interaction reproducing the early blue emission and r-process heating reproducing the late red emission. This supports claims that late low-power engine activity after prompt emission may be common. We encourage early follow-up observations of future nearby GRBs and compact binary merger events to reveal more about the central engine of GRBs and r-process events.

Observing black hole mergers beyond the pair-instability mass gap with next-generation gravitational wave detectors

Stellar evolution predicts the existence of a mass gap for black hole remnants produced by pair-instability supernova dynamics, whose lower and upper edges are very uncertain. We study the possibility of constraining the location of the upper end of the pair-instability mass gap, which is believed to appear around mmin∼130M⊙, using gravitational wave observations of compact binary mergers with next-generation ground-based detectors. While high metallicity may not allow for the formation of first-generation black holes on the “far side” beyond the gap, metal-poor environments containing population III stars could lead to such heavy black hole mergers. We show that, even in the presence of contamination from other merger channels, next-generation detectors will measure the location of the upper end of the mass gap with a relative precision close to Δmmin/mmin≃4%(Ndet/100)−1/2 at 90% CL, where Ndet is the number of detected mergers with both members beyond the gap. These future observations could reduce current uncertainties in nuclear and astrophysical processes controlling the location of the gap. Published by the American Physical Society 2024

Convolutional neural networks for signal detection in real LIGO data

Searching the data of gravitational-wave detectors for signals from compact binary mergers is a computationally demanding task. Recently, machine-learning algorithms have been proposed to address current and future challenges. However, the results of these publications often differ greatly due to differing choices in the evaluation procedure. The Machine Learning Gravitational-Wave Search Challenge was organized to resolve these issues and produce a unified framework for machine-learning search evaluation. Six teams submitted contributions, four of which are based on machine-learning methods, and two are state-of-the-art production analyses. This paper describes the submission from the team TPI FSU Jena and its updated variant. We also apply our algorithm to real O3b data and recover the relevant events of the GWTC-3 catalog. Published by the American Physical Society 2024

Unresolved stochastic background from compact binary mergers detectable by next-generation ground-based gravitational-wave observatories

Unresolved stochastic background from compact binary mergers detectable by next-generation ground-based gravitational-wave observatories

Alternative statistical interpretation for the apparent plateaus in the duration distributions of gamma-ray bursts

Context. The existence of a plateau in the short-duration tail of the observed distribution of cosmological long-soft gamma ray bursts (LGRBs) is posited to be the first direct evidence of collapsars. A similar plateau in the short-duration tail of the observed duration distribution of short-hard gamma ray bursts (SGRBs) has been suggested as evidence of compact binary mergers. Aims. We present an equally plausible alternative interpretation for this evidence, based on a purely statistical approach. Methods. Specifically, we show that the observed plateau in the short-duration tail of the duration distribution of LGRBs can naturally occur in the statistical distributions of strictly positive physical quantities, exacerbated by the effects of mixing with the duration distribution of SGRBs, observational selection effects, and data aggregation (e.g., binning) methodologies. The observed plateau in the short-duration tail of the observed distributions of SGRBs may likewise result from a combination of sample incompleteness and inhomogeneous binning of data. We further confirm the impact of these factors on the observation of a plateau in the duration distributions of GRBs through extensive numerical Monte Carlo simulations. Results. This analysis corroborates and strengthens a purely statistical and sample-incompleteness interpretation of the observed plateau in the duration distribution of LGRBs and SGRBs, without invoking the physics of collapsars or jet-propagation through the stellar envelope.

Can we identify primordial black holes? Tidal tests for subsolar-mass gravitational-wave observations

The detection of a subsolar object in a compact binary merger is regarded as one of the smoking gun signatures of a population of primordial black holes (PBHs). We critically assess whether these systems could be distinguished from stellar binaries, for example composed of white dwarfs or neutron stars, which could also populate the subsolar mass range. At variance with PBHs, the gravitational-wave signal from stellar binaries is affected by tidal effects, which dramatically grow for moderately compact stars as those expected in the subsolar range. We forecast the capability of constraining tidal effects of putative subsolar neutron star binaries with current and future LIGO-Virgo-KAGRA (LVK) sensitivities as well as next-generation experiments. We show that, should LVK O4 run observe subsolar neutron-star mergers, it could measure the (large) tidal effects with high significance. In particular, for subsolar neutron-star binaries, O4 and O5 projected sensitivities would allow measuring the effect of tidal disruption on the waveform in a large portion of the parameter space, also constraining the tidal deformability at O(10%) level, thus excluding a primordial origin of the binary. Vice versa, for subsolar PBH binaries, model-agnostic upper bounds on the tidal deformability can rule out neutron stars or more exotic competitors. Assuming events similar to the subthreshold candidate SSM200308 reported in LVK O3b data are PBH binaries, O4 projected sensitivity would allow ruling out the presence of neutron-star tidal effects at ≈3σ CL, thus strengthening the PBH hypothesis. Future experiments would lead to even stronger (&gt;5σ) conclusions on potential discoveries of this kind. Published by the American Physical Society 2024

First Constraints on Compact Binary Environments from LIGO-Virgo Data.

The LIGO-Virgo analyses of signals from compact binary mergers observed so far have assumed isolated binary systems in a vacuum, neglecting the potential presence of astrophysical environments. We present here the first investigation of environmental effects on each of the events of GWTC-1 and two low-mass events from GWTC-2. We find no evidence for the presence of environmental effects. Most of the events decisively exclude the scenario of dynamical fragmentation of massive stars as their formation channel. GW170817 results in the most stringent upper bound on the medium density (≲21 g/cm^{3}). We find that environmental effects can substantially bias the recovered parameters in the vacuum model, even when these effects are not detectable. We forecast that the Einstein Telescope and B-DECIGO will be able to probe the environmental effects of accretion disks and superradiant boson clouds on compact binaries.