Binary Black Hole Mergers Research Articles

Searching for Gravitational-wave Counterparts Using the Transiting Exoplanet Survey Satellite

In 2017, the LIGO and Virgo gravitational-wave (GW) detectors, in conjunction with electromagnetic (EM) astronomers, observed the first GW multimessenger astrophysical event, the binary neutron star (BNS) merger GW170817. This marked the beginning of a new era in multimessenger astrophysics. To discover further GW multimessenger events, we explore the synergies between the Transiting Exoplanet Survey Satellite (TESS) and GW observations triggered by the LIGO–Virgo–KAGRA Collaboration (LVK) detector network. TESS's extremely wide field of view (∼2300 deg2) means that it could overlap with large swaths of GW localizations, which often span hundreds of square degrees or more. In this work, we use a recently developed transient detection pipeline to search TESS data collected during the LVK’s third observing run, O3, for any EM counterparts. We find no obvious counterparts brighter than about 17th magnitude in the TESS bandpass. Additionally, we present end-to-end simulations of BNS mergers, including their detection in GWs and simulations of light curves, to identify TESS's kilonova discovery potential for the LVK's next observing run (O4). In the most optimistic case, TESS will observe up to one GW-found BNS merger counterpart per year. However, TESS may also find up to five kilonovae that did not trigger the LVK network, emphasizing that EM-triggered GW searches may play a key role in future kilonova detections. We also discuss how TESS can help place limits on EM emission from binary black hole mergers and rapidly exclude large sky areas for poorly localized GW events.

Neural Importance Sampling for Rapid and Reliable Gravitational-Wave Inference.

We combine amortized neural posterior estimation with importance sampling for fast and accurate gravitational-wave inference. We first generate a rapid proposal for the Bayesian posterior using neural networks, and then attach importance weights based on the underlying likelihood and prior. This provides (1)a corrected posterior free from network inaccuracies, (2)a performance diagnostic (the sample efficiency) for assessing the proposal and identifying failure cases, and (3)an unbiased estimate of the Bayesian evidence. By establishing this independent verification and correction mechanism we address some of the most frequent criticisms against deep learning for scientific inference. We carry out a large study analyzing 42 binary black hole mergers observed by LIGO and Virgo with the SEOBNRv4PHM and IMRPhenomXPHM waveform models. This shows a median sample efficiency of ≈10% (2 orders of magnitude better than standard samplers) as well as a tenfold reduction in the statistical uncertainty in the log evidence. Given these advantages, we expect a significant impact on gravitational-wave inference, and for this approach to serve as a paradigm for harnessing deep learning methods in scientific applications.

Constraining the number of spacetime dimensions from GWTC-3 binary black hole mergers

In modified gravity models that allow for additional noncompact spacetime dimensions, energy from gravitational waves can leak into these extra spacetime dimensions, leading to a reduction in the amplitude of the observed gravitational waves, and thus are a source of potential systematics in the inferred luminosity distances to gravitational wave sources. Since binary black hole (BBH) mergers are standard sirens, we use the pair-instability supernova mass gap and its predicted features to determine a mass scale in order to break the mass-redshift degeneracy and thus infer the redshift of the source. We simultaneously fit the BBH population and the extra spacetime dimensions parameters from gravitational leakage models using BBH observations from the recently released GWTC-3 catalog. We set constraints on the number of spacetime dimensions and find that $D=3.9{3}_{\ensuremath{-}0.05}^{+0.08}$ at 68% C.L. for models that are independent of a screening scale, finding that the GWTC-3 constraint is as competitive as that set from GW170817 and its electromagnetic counterpart. For models where gravity leaks below a certain screening scale ${R}_{c}$, we find $D=4.4{9}_{\ensuremath{-}0.87}^{+1.63}$ and ${\mathrm{log}}_{10}\text{ }{R}_{c}/\mathrm{Mpc}=4.7{5}_{\ensuremath{-}0.99}^{+0.86}$ with a transition steepness ${\mathrm{log}}_{10}\text{ }n=0.8{5}_{\ensuremath{-}0.86}^{+0.76}$ for the leakage, which for the first time are constrained jointly with the BBH population at cosmological distances. These constraints are consistent with General Relativity (GR), where gravitational waves propagate in $D=3+1$ spacetime dimensions. Using the BBH population to probe modifications to standard cosmological models provides an independent test of GR that does not rely on any electromagnetic information but purely on gravitational wave observations.

Model exploration in gravitational-wave astronomy with the maximum population likelihood

Hierarchical Bayesian inference is an essential tool for studying the population properties of compact binaries with gravitational waves. The basic premise is to infer the unknown prior distribution of binary black hole and/or neutron star parameters such component masses, spin vectors, and redshift. These distributions shed light on the fate of massive stars, how and where binaries are assembled, and the evolution of the Universe over cosmic time. Hierarchical analyses model the binary black hole population using a prior distribution conditioned on hyperparameters, which are inferred from the data. However, a misspecified model can lead to faulty astrophysical inferences. In this paper we answer the question: given some data, which prior distribution---from the set of all possible prior distributions---produces the largest possible population likelihood? This distribution (which is not a true prior) is $\overline{)\ensuremath{\pi}}$ (pronounced ``pi stroke''), and the associated maximum population likelihood is $\overline{)\mathcal{L}}$ (pronounced ``L stroke''). The structure of $\overline{)\ensuremath{\pi}}$ is a linear superposition of delta functions, a result which follows from Carath\'eodory's theorem. We show how $\overline{)\ensuremath{\pi}}$ and $\overline{)\mathcal{L}}$ can be used for model exploration/criticism. We apply this $\overline{)\mathcal{L}}$ formalism to study the population of binary black hole mergers observed in the LIGO-Virgo-KAGRA Collaboration's third gravitational-wave transient catalog. Based on our results, we discuss possible improvements for gravitational-wave population models.

Measuring dark matter spikes around primordial black holes with Einstein Telescope and Cosmic Explorer

Future ground-based gravitational wave observatories will be ideal probes of the environments surrounding black holes with masses $1\ensuremath{-}10{M}_{\ensuremath{\bigodot}}$. Binary black hole mergers with mass ratios of order $q={m}_{2}/{m}_{1}\ensuremath{\lesssim}{10}^{\ensuremath{-}3}$ can remain in the frequency band of such detectors for months or years, enabling precision searches for modifications of their gravitational waveforms with respect to vacuum inspirals. As a concrete example of an environmental effect, we consider here a population of binary primordial black holes which are expected to be embedded in dense cold dark matter spikes. We provide a viable formation scenario for these systems compatible with all observational constraints and predict upper and lower limits on the merger rates of small-mass-ratio pairs. Given a detected signal of one such system by either Einstein Telescope or Cosmic Explorer, we show that the properties of the binary and of the dark matter spike can be measured to excellent precision with one week's worth of data, if the effect of the dark matter spike on the waveform is taken into account. However, we show that there is a risk of biased parameter inference or missing the events entirely if the effect of the predicted dark matter overdensity around these objects is not properly accounted for.

Study of the intermediate mass ratio black hole binary merger up to 1000:1 with numerical relativity

We explicitly demonstrate that current numerical relativity techniques are able to accurately evolve black hole binaries with mass ratios of the order of 1000:1. This proof of principle is relevant for future third generation gravitational wave detectors and space mission LISA, as by purely numerical methods we would be able to accurately compute gravitational waves from the last stages of black hole mergers, as directly predicted by general relativity. We perform a sequence of simulations in the intermediate to small mass ratio regime, , with the small hole starting from rest at a proper distance . We compare these headon full numerical evolutions with the corresponding semianalytic point particle perturbative results finding an impressive agreement for the total gravitational radiated energy and linear momentum as well as for the waveform spectra. We display numerical convergence of the results and identify the minimal numerical resolutions required to accurately solve for these very low amplitude gravitational waves. This work represents a first step towards the considerable challenge of applying numerical-relativity waveforms to interpreting gravitational-wave observations by LISA and next-generation ground-based gravitational-wave detectors.

Precision cosmology with primordial GW backgrounds in presence of astrophysical foregrounds

The era of Gravitational-Wave (GW) astronomy will grant the detection of the astrophysical GW background from unresolved mergers of binary black holes, and the prospect of probing the presence of primordial GW backgrounds. In particular, the low-frequency tail of the GW spectrum for causally-generated primordial signals (like a phase transition) offers an excellent opportunity to measure unambiguously cosmological parameters as the equation of state of the universe, or free-streaming particles at epochs well before recombination. We discuss whether this programme is jeopardised by the uncertainties on the astrophysical GW foregrounds that coexist with a primordial background. We detail the motivated assumptions under which the astrophysical foregrounds can be assumed to be known in shape, and only uncertain in their normalisation. In this case, the sensitivity to a primordial signal can be computed by a simple and numerically agile procedure, where the optimal filter function subtracts the components of the astrophysical foreground that are close in spectral shape to the signal. We show that the degradation of the sensitivity to the signal in presence of astrophysical foregrounds is limited to a factor of a few, and only around the frequencies where the signal is closer to the foregrounds. Our results highlight the importance of modelling the contributions of eccentric or intermediate-mass black hole binaries to the GW background, to consolidate the prospects to perform precision cosmology with primordial GW backgrounds.

Search for neutrino counterparts to the gravitational wave sources from LIGO/Virgo O3 run with the ANTARES detector

Since 2015 the LIGO and Virgo interferometers have detected gravitational waves from almost one hundred coalescences of compact objects (black holes and neutron stars). This article presents the results of a search performed with data from the ANTARES telescope to identify neutrino counterparts to the gravitational wave sources detected during the third LIGO/Virgo observing run and reported in the catalogues GWTC-2, GWTC-2.1, and GWTC-3. This search is sensitive to all-sky neutrinos of all flavours and of energies > 100 GeV, thanks to the inclusion of both track-like events (mainly induced by νμ charged-current interactions) and shower-like events (induced by other interaction types). Neutrinos are selected if they are detected within ± 500 s from the GW merger and with a reconstructed direction compatible with its sky localisation. No significant excess is found for any of the 80 analysed GW events, and upper limits on the neutrino emission are derived. Using the information from the GW catalogues and assuming isotropic emission, upper limits on the total energy E tot,ν emitted as neutrinos of all flavours and on the ratio fν = E tot,ν /E GW between neutrino and GW emissions are also computed. Finally, a stacked analysis of all the 72 binary black hole mergers (respectively the 7 neutron star-black hole merger candidates) has been performed to constrain the typical neutrino emission within this population, leading to the limits: E tot,ν < 4.0 × 1053 erg and fν < 0.15 (respectively, E tot,ν < 3.2 × 10^53 erg and fν < 0.88) for E -2 spectrum and isotropic emission. Other assumptions including softer spectra and non-isotropic scenarios have also been tested.

4-OGC: Catalog of Gravitational Waves from Compact Binary Mergers

We present the fourth Open Gravitational-wave Catalog (4-OGC) of binary neutron star (BNS), binary black hole (BBH), and neutron star–black hole (NSBH) mergers. The catalog includes observations from 2015 to 2020 covering the first through third observing runs (O1, O2, O3a, and O3b) of Advanced LIGO and Advanced Virgo. The updated catalog includes seven BBH mergers that were not previously reported with high significance during O3b for a total of 94 observations: 90 BBHs, 2 NSBHs, and 2 BNSs. The most confident new detection, GW200318_191337, has component masses and its redshift of (90% credible interval) may make it the most distant merger so far. We estimate the merger rate of BBH sources, assuming a power-law mass distribution containing an additive Gaussian peak, to be Gpc−3 yr−1 at a redshift of z = 0 (0.2). For BNS and NSBH sources, we estimate a merger rate of Gpc−3 yr−1 and Gpc−3 yr−1, respectively, assuming the known sources are representative of the total population. We provide reference parameter estimates for each of these sources using an up-to-date model accounting for instrumental calibration uncertainty. The corresponding data release also includes our full set of subthreshold candidates.

Coalescing black hole binaries from globular clusters: mass distributions and comparison to gravitational wave data from GWTC-3

ABSTRACTWe use our cluster population model, cBHBd, to explore the mass distribution of merging black hole binaries formed dynamically in globular clusters. We include in our models the effect of mass growth through hierarchical mergers and compare the resulting distributions to those inferred from the third gravitational wave transient catalogue. We find that none of our models can reproduce the peak at m1 ≃ 10 M⊙ in the primary black hole mass distribution that is inferred from the data. This disfavours a scenario where most of the sources are formed in globular clusters. On the other hand, a globular cluster origin can account for the inferred secondary peak at m1 ≃ 35 M⊙, which requires that the most massive clusters form with half-mass densities $\rho _{\rm h,0}\gtrsim 10^4~{\rm M}_{\odot} \, {\rm pc}^{-3}$. Finally, we find that the lack of a high-mass cut-off in the inferred mass distribution can be explained by the repopulation of an initial mass gap through hierarchical mergers. Matching the inferred merger rate above ≃50 M⊙ requires both initial cluster densities $\rho _{\rm h,0}\gtrsim 10^4~{\rm M}_{\odot} \, {\rm pc}^{-3}$, and that black holes form with nearly zero spin. A hierarchical merger scenario makes specific predictions for the appearance and position of multiple peaks in the black hole mass distribution, which can be tested against future data.

Population of Merging Compact Binaries Inferred Using Gravitational Waves through GWTC-3

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc−3 yr−1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc−3 yr−1, assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc−3 yr−1 at a fiducial redshift (z=0.2). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to (1+z)κ with κ=2.9−1.8+1.7 for z≲1. Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2−0.2+0.1 to 2.0−0.3+0.3M⊙. We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3−0.5+0.3 and 27.9−1.8+1.9M⊙. While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60M⊙, which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χi≈0.25. While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidence of antialigned spins among the binary population. We observe an increase in spin magnitude for systems with more unequal-mass ratio. We also observe evidence of misalignment of spins relative to the orbital angular momentum.22 MoreReceived 4 February 2022Revised 28 October 2022Accepted 19 December 2022DOI:https://doi.org/10.1103/PhysRevX.13.011048Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasGravitational wave sourcesGravitational wavesTechniquesGravitational wave detectionGravitation, Cosmology & Astrophysics

Outlook for detecting the gravitational-wave displacement and spin memory effects with current and future gravitational-wave detectors

Gravitational-wave memory effects arise from nonoscillatory components of gravitational-wave signals, and they are predictions of general relativity in the nonlinear regime that have close connections to the asymptotic properties of isolated gravitating systems. There are many types of memory effects that have been studied in the literature. In this paper we focus on the ``displacement'' and ``spin'' memories, which are expected to be the largest of these effects from sources such as the binary black hole mergers which have already been detected by LIGO and Virgo. The displacement memory is a change in the relative separation of two initially comoving observers due to a burst of gravitational waves, whereas the spin memory is a portion of the change in relative separation of observers with initial relative velocity. As both of these effects are small, LIGO, Virgo, and KAGRA can only detect memory effects from individual events that are much louder (and thus rarer) than those that have been detected so far. By combining data from multiple events, however, these effects could be detected in a population of binary mergers. In this paper, we present new forecasts for how long current and future detectors will need to operate in order to measure these effects from populations of binary black hole systems that are consistent with the populations inferred from the detections from LIGO and Virgo's first three observing runs. We find that a second-generation detector network of LIGO, Virgo, and KAGRA operating at the O4 (``design'') sensitivity for 1.5 years and then operating at the O5 (``plus'') sensitivity for an additional 1.5 years can detect the displacement memory. For Cosmic Explorer, we find that displacement memory could be detected for individual loud events, and that the spin memory could be detected in a population after 5 years of observation time.

Searching for clues of past binary supermassive black hole mergers in nuclear star clusters

ABSTRACTGalaxy mergers are common processes in the Universe. As a large fraction of galaxies hosts at their centres a central supermassive black hole (SMBH), mergers can lead to the formation of a supermassive black hole binary (SMBHB). The formation of such a binary is more efficient when the SMBHs are embedded in a nuclear star cluster (NSC). NSCs are dense and massive stellar clusters present in the majority of the observed galaxies. Their central densities can reach up to $10^7\, {\rm M_{\odot }}\,{\rm pc^{-3}}$ and their masses can be as large as a few $10^7\, {\rm M_{\odot }}$. The direct detection of an SMBHB is observationally challenging. In this work, we illustrate how the large-scale structural and dynamical properties of an NSC can help to identify nucleated galaxies that recently went through a merger that possibly led to the formation of a central SMBHB. Our models show that the merger can imprint signatures on the shape, density profile, rotation, and velocity structure of the NSC. The strength of the signatures depends on the mass ratio between the SMBHs and on the orbital initial conditions of the merger. In addition, the number of hypervelocity stars produced in the mergers is linked to the SMBHB properties. The merger can also contribute to the formation of the nuclear stellar disc of the galaxy.

Numerical relativity higher order gravitational waveforms of eccentric, spinning, nonprecessing binary black hole mergers

We use the open source, community-driven, numerical relativity software, the Einstein Toolkit to study the physics of eccentric, spinning, nonprecessing binary black hole mergers with mass-ratios $q=\{2, 4, 6\}$, individual dimensionless spin parameters $\chi_{1z}=\pm0.6$, $\chi_{2z}=\pm0.3$, that include higher order gravitational wave modes $\ell\leq4$, except for memory modes. Assuming stellar mass binary black hole mergers that may be detectable by the advanced LIGO detectors, we find that including modes up to $\ell=4$ increases the signal-to-noise of compact binaries between $3.5\%$ to $35\%$, compared to signals that only include the $\ell=|m|=2$ mode. We use two waveform models, TEOBResumS and SEOBNRE, which incorporate spin and eccentricity corrections in the waveform dynamics, to quantify the orbital eccentricity of our numerical relativity catalog in a gauge-invariant manner through fitting factor calculations. Our findings indicate that the inclusion of higher order wave modes has a measurable effect in the recovery of moderately and highly eccentric black hole mergers, and thus it is essential to develop waveform models and signal processing tools that accurately describe the physics of these astrophysical sources.

Reconstruction of power spectrum of primordial curvature perturbations on small scales from primordial black hole binaries scenario of LIGO/VIRGO detection

As a candidate bound for the Binary Black Hole (BBH) merger events detected by LIGO/Virgo, Primordial Black Holes (PBHs) provide a useful tool to investigate the primordial curvature perturbations on small scales. Using the GWTC-1 to GWTC-3 catalogs, under the scenario that PBHs originate from large primordial curvature perturbations on small scales during inflationary epoch, we for the first time reconstruct the power spectrum of primordial curvature perturbations on small scales. It is found that the value of the amplitude of the primordial power spectrum is enhanced to $\mathcal{O}(10^{-2})$ on scales $\mathcal{O}(1)$ pc. This may imply the validity of PBH as a possible BBH merger candidate.

Binary black hole mergers from Population III stars: uncertainties from star formation and binary star properties

Binary black hole mergers from Population III stars: uncertainties from star formation and binary star properties

Binary black hole mergers from Population III stars: uncertainties from star formation and binary star properties

Binary black hole mergers from Population III stars: uncertainties from star formation and binary star properties

Inferring eccentricity evolution from observations of coalescing binary black holes

The origin and formation of stellar-mass binary black holes remains an open question that can be addressed by precise measurements of the binary and orbital parameters from their gravitational-wave signal. Such binaries are expected to circularize due to the emission of gravitational waves as they approach merger. However, depending on their formation channel, some binaries could have a non-negligible eccentricity when entering the frequency band of current gravitational-wave detectors. In order to measure eccentricity in an observed gravitational-wave signal, accurate waveform models that describe binaries in eccentric orbits are necessary. In this work we demonstrate the efficacy of the improved TEOBResumS waveform model for eccentric coalescing binaries with aligned spins. We first validate the model against mock signals of aligned-spin binary black hole mergers and quantify the impact of eccentricity on the estimation of other intrinsic binary parameters. We then perform a fully Bayesian reanalysis of GW150914 with the eccentric waveform model. We find (i) that the model is reliable for aligned-spin binary black holes and (ii) that GW150914 is consistent with a non-eccentric merger although we cannot rule out small values of initial eccentricity at a reference frequency of $20$ Hz. Finally, we present a systematic method to measure the eccentricity and its evolution directly from the gravitational-wave posterior samples. Such an estimator is useful when comparing results from different analyses as the definition of eccentricity may differ between models. Our scheme can be applied even in the case of small eccentricities and can be adopted straightforwardly in post-processing to allow for direct comparison between models.

The appearance of a merging binary black hole very close to a spinning supermassive black hole

ABSTRACTThe mass and distance of a binary black hole (BBH) are fundamental parameters to measure in gravitational wave (GW) astronomy. It is well-known that the measurement is affected by cosmological redshift, and recent works also showed that Doppler and gravitational redshifts could further affect the result if the BBH coalesces close to a supermassive black hole (SMBH). Here we consider the additional lensing effect induced by the nearby SMBH on the measurement. We compute the null geodesics originating within 10 gravitational radii of a Kerr SMBH to determine the redshift and magnification of the GWs emitted by the BBH. We find a positive correlation between redshift and demagnification, which results in a positive correlation between the mass and distance of the BBH in the detector frame. More importantly, we find a higher probability for the signal to appear redshifted and demagnified to a distant observer, rather than blueshifted and magnified. Based on these results, we show that a binary at a cosmological redshift of zcos = (10−2–10−1) and composed of BHs of $(10\!-\!20)\, \mathrm{M}_\odot$ could masquerade as a BBH at a redshift of zcos ∼ 1 and containing BHs as large as $(44\!-\!110)\, \mathrm{M}_\odot$. In the case of extreme demagnification, we also find that the same BBH could appear to be at zcos &amp;gt; 10 and contain subsolar-mass BHs. Such an effect, if not accounted for, could bias our understanding of the origin of the BHs detected via GWs.

On the gravitational lensing interpretation of three gravitational wave detections in the mass gap by LIGO and Virgo

ABSTRACT We search for gravitational wave (GW) events from LIGO-Virgo’s third run that may have been affected by gravitational lensing. Gravitational lensing delays the arrival of GWs, and alters their amplitude – thus biasing the inferred progenitor masses. This would provide a physically well-understood interpretation of GW detections in the ‘mass gap’ between neutron stars and black holes, as gravitationally lensed binary neutron star (BNS) mergers. We selected three GW detections in LIGO-Virgo’s third run for which the probability of at least one of the constituent compact objects being in the mass gap was reported as high with low latency – i.e. candidate lensed BNS mergers. Our observations of powerful strong lensing clusters located adjacent to the peak of their sky localization error maps reached a sensitivity $\rm AB\simeq 25.5$ in the z′ band with the GMOS instruments on the Gemini telescopes, and detected no candidate lensed optical counterparts. We combine recent kilonova light-curve models with recent predictions of the lensed BNS population and the properties of the objects that we followed up to show that realistic optical counterparts were detectable in our observations. Further detailed analysis of two of the candidates suggests that they are a plausible pair of images of the same low-mass binary black hole merger, lensed by a local galaxy or small group of galaxies. This further underlines that access to accurate mass information with low latency would improve the efficiency of candidate lensed BNS selection.