Binary Black Hole Coalescence Research Articles

Perturbations of spinning black holes in dynamical Chern-Simons gravity: Slow rotation equations

The detection of gravitational waves resulting from the coalescence of binary black holes by the LIGO-Virgo-Kagra Collaboration has inaugurated a new era in gravitational physics. These gravitational waves provide a unique opportunity to test Einstein’s general relativity and its modifications in the regime of extreme gravity. A significant aspect of such tests involves the study of the ringdown phase of gravitational waves from binary black hole coalescence, which can be decomposed into a superposition of various quasinormal modes. In general relativity, the spectra of quasinormal modes depend on the mass, spin, and charge of the final black hole, but they can also be influenced by additional properties of the black hole spacetime, as well as corrections to the general theory of relativity. In this work, we focus on a specific modified theory known as dynamical Chern-Simons gravity. We employ the modified Teukolsky formalism developed in a previous study and lay down the foundations to investigate perturbations of slowly rotating black holes admitted by the theory. Specifically, we derive the master equations for the Ψ0 and Ψ4 Weyl scalar perturbations that characterize the radiative part of gravitational perturbations, as well as the master equation for the scalar field perturbations. We employ metric reconstruction techniques to obtain explicit expressions for all relevant quantities. Finally, by leveraging the properties of spin-weighted spheroidal harmonics to eliminate the angular dependence from the evolution equations, we derive two, radial, second-order, ordinary differential equations for Ψ0 and Ψ4, respectively. These two equations are coupled to another radial, second-order, ordinary differential equation for the scalar field perturbations. This work is the first attempt to derive a master equation for black holes in dynamical Chern-Simons gravity using curvature perturbations. The master equations we obtain can then be numerically integrated to obtain the quasinormal mode spectrum of slowly rotating black holes in this theory, making progress in the study of ringdown in dynamical Chern-Simons gravity. Published by the American Physical Society 2024

The Evolution of Accreting Binaries: From Brown Dwarfs to Supermassive Black Holes

Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios (q ≡ M 2/M 1) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios H/r ∈ {0.1, 0.05, 0.03} and vary both the magnitude and spatial dependence of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with 0.01 ≤ q ≤ 0.04 outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outward. However, in a number of cases, the disks accreting onto binaries with mass ratios ∼0.07 fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the ∼10% level even for mass ratios of 0.01. We demonstrate that the spatial dependence of the viscosity (e.g., α vs. constant ν) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at q ∼ 0.01 and sustain deep, asymmetric cavities.

Unveiling the Merger Structure of Black Hole Binaries in Generic Planar Orbits.

The precise modeling of binary black hole coalescences in generic planar orbits is a crucial step to disentangle dynamical and isolated binary formation channels through gravitational-wave observations. The merger regime of such coalescences exhibits a significantly higher complexity compared to the quasicircular case, and cannot be readily described through standard parametrizations in terms of eccentricity and anomaly. In the spirit of the effective one body formalism, we build on the study of the test-mass limit, and introduce a new modeling strategy to describe the general-relativistic dynamics of two-body systems in generic orbits. This is achieved through gauge-invariant combinations of the binary energy and angular momentum, such as a dynamical "impact parameter" at merger. These variables reveal simple "quasi-universal" structures of the pivotal merger parameters, allowing us to build an accurate analytical representation of generic (bounded and dynamically bounded) orbital configurations. We demonstrate the validity of these analytical relations using 311 numerical simulations of bounded noncircular binaries with progenitors from the RIT and SXS catalogs, together with a custom dataset of dynamical captures generated using the Einstein Toolkit, and test-mass data in bound orbits. Our modeling strategy lays the foundations of accurate and complete waveform models for systems in arbitrary orbits, bolstering observational explorations of dynamical formation scenarios and the discovery of new classes of gravitational wave sources.

Non-stationary astrophysical stochastic gravitational-wave background: a new probe to the high-redshift population of binary black holes

ABSTRACT The astrophysical stochastic gravitational-wave background (SGWB) originates from the mergers of compact binary objects that are otherwise undetected as individual events, along with other sources such as supernovae, magnetars, etc. The individual gravitational-wave (GW) signal is time-varying over a time-scale that depends on the chirp mass of the coalescing binaries. Another time-scale that plays a role is the time-scale at which the sources repeat, which depends on the merger rate. The combined effect of these two leads to a breakdown of the time translation symmetry of the observed SGWB and a correlation between different frequency modes in the signal covariance matrix of the SGWB. Using an ensemble of SGWB due to binary black hole coalescence, calculated using simulations of different black hole mass distributions and merger rates, we show how the structure of the signal covariance matrix varies. This structure in the signal covariance matrix brings additional information about the sources on top of the power spectrum. We show that there is a significant improvement in the figure of merit by using this additional information in comparison to only power spectrum estimation for the LIGO–Virgo–KAGRA (LVK) network of detectors with the design sensitivity noise with 2 yr of observation. The inclusion of the off-diagonal correlation in the covariance of the SGWB in the data analysis pipelines will be beneficial in the quest for the SGWB signal in LVK frequency bands as well as in lower frequencies and in getting an insight into its origin.

Universality in binary black hole dynamics: An integrability conjecture

The waveform of a binary black hole coalescence appears to be both simple and universal. In this essay we argue that the dynamics should admit a separation into “fast and slow” degrees of freedom, such that the latter are described by an integrable system of equations, accounting for the simplicity and universality of the waveform. Given that Painlevé transcendents are a smoking gun of integrable structures, we propose the Painlevé-II transcendent as the key structural element threading a hierarchy of asymptotic models aiming at capturing different (effective) layers in the dynamics. Ward’s conjecture relating integrable and (anti-)self-dual solutions can provide the avenue to encode background binary black hole data in (nonlocal) twistor structures.

Importance nested sampling with normalising flows

Abstract We present an improved version of the nested sampling algorithm nessai in which the core algorithm is modified to use importance weights. In the modified algorithm, samples are drawn from a mixture of normalising flows and the requirement for samples to be independently and identically distributed (i.i.d.) according to the prior is relaxed. Furthermore, it allows for samples to be added in any order, independently of a likelihood constraint, and for the evidence to be updated with batches of samples. We call the modified algorithm i-nessai. We first validate i-nessai using analytic likelihoods with known Bayesian evidences and show that the evidence estimates are unbiased in up to 32 dimensions. We compare i-nessai to standard nessai for the analytic likelihoods and the Rosenbrock likelihood, the results show that i-nessai is consistent with nessai whilst producing more precise evidence estimates. We then test i-nessai on 64 simulated gravitational-wave signals from binary black hole coalescence and show that it produces unbiased estimates of the parameters. We compare our results to those obtained using standard nessai and dynesty and find that i-nessai requires 2.68 and 13.3 times fewer likelihood evaluations to converge, respectively. We also test i-nessai of an 80 s simulated binary neutron star signal using a reduced-order-quadrature basis and find that, on average, it converges in 24 min, whilst only requiring 1.01 × 10 6 likelihood evaluations compared to 1.42 × 10 6 for nessai and 4.30 × 10 7 for dynesty. These results demonstrate that i-nessai is consistent with nessai and dynesty whilst also being more efficient.

Inferring the astrophysical population of gravitational wave sources in the presence of noise transients

ABSTRACT The global network of interferometric gravitational wave (GW) observatories (LIGO, Virgo, KAGRA) has detected and characterized nearly 100 mergers of binary compact objects. However, many more real GWs are lurking sub-threshold, which need to be sifted from terrestrial-origin noise triggers (known as glitches). Because glitches are not due to astrophysical phenomena, inference on the glitch under the assumption it has an astrophysical source (e.g. binary black hole coalescence) results in source parameters that are inconsistent with what is known about the astrophysical population. In this work, we show how one can extract unbiased population constraints from a catalogue of both real GW events and glitch contaminants by performing Bayesian inference on their source populations simultaneously. In this paper, we assume glitches come from a specific class with a well-characterized effective population (blip glitches). We also calculate posteriors on the probability of each event in the catalogue belonging to the astrophysical or glitch class, and obtain posteriors on the number of astrophysical events in the catalogue, finding it to be consistent with the actual number of events included.

LISA science performance in observations of short-lived signals from massive black hole binary coalescences

LISA science performance in observations of short-lived signals from massive black hole binary coalescences

Measuring the ringdown scalar polarization of gravitational waves in Einstein-scalar-Gauss-Bonnet gravity

We model the scalar waves produced during the ringdown stage of binary black hole coalescence in Einstein scalar Gauss-Bonnet (EsGB) gravity, using numerical relativity simulations of the theory in the decoupling limit. Through a conformal coupling of the scalar field to the metric in the matter-field action, we show that the gravitational waves in this theory can have a scalar polarization. We model the scalar quasi-normal modes of the ringdown signal in EsGB gravity, and quantify the extent to which current and future gravitational wave detectors could observe the spectrum of scalar radiation emitted during the ringdown phase of binary black hole coalescence. We find that within the limits of the theory's coupling parameters set by current theoretical and observational constraints, the scalar ringdown signal from black hole remnants in the $10^1 - 10^3 \, M_{\odot}$ mass range is expected to be well below the detectability threshold with the current network of gravitational-wave detectors (LIGO-Virgo-KAGRA), but is potentially measurable with next-generation detectors such as the Einstein Telescope.

Analyses of overlapping gravitational wave signals using hierarchical subtraction and joint parameter estimation

ABSTRACT In the coming years, third-generation detectors such as Einstein Telescope and Cosmic Explorer will enter the network of ground-based gravitational-wave detectors. Their current design predicts a significantly improved sensitivity band with a lower minimum frequency than existing detectors. This, combined with the increased arm length, leads to two major effects: the detection of more signals and the detection of longer signals. Both will result in a large number of overlapping signals. It has been shown that such overlapping signals can lead to biases in the recovered parameters, which would adversely affect the science extracted from the observed binary merger signals. In this work, we analyse overlapping binary black hole coalescences with two methods to analyse multisignal observations: hierarchical subtraction and joint parameter estimation. We find that these methods enable a reliable parameter extraction in most cases and that joint parameter estimation is usually more precise but comes with higher computational costs.

Rapid search for massive black hole binary coalescences using deep learning

The coalescences of massive black hole binaries are one of the main targets of space-based gravitational wave observatories. Such gravitational wave sources are expected to be accompanied by electromagnetic emissions. Low latency detection of the massive black hole mergers provides a start point for a global-fit analysis to explore the large parameter space of signals simultaneously being present in the data but at great computational cost. To alleviate this issue, we present a deep learning method for rapidly searching for signals of massive black hole binaries in gravitational wave data. Our model is capable of processing a year of data, simulated from the LISA data challenge, in only several seconds, while identifying all coalescences of massive black hole binaries with no false alarms. We further demonstrate that the model shows robust resistance to a wide range of generalization cases, including various waveform families and updated instrumental configurations. This method offers an effective approach that combines advances in artificial intelligence to open a new pathway for space-based gravitational wave observations.

Follow-up Survey for the Binary Black Hole Merger GW200224_222234 Using Subaru/HSC and GTC/OSIRIS

The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (also known as S200224ca) and classified as a binary-black hole coalescence, on 2020 February 24. Given its relatively small localization skymap (71 deg2 for a 90% credible region; revised to 50 deg2 in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the r2 and z bands. Observations were conducted on 2020 February 25 and 28 and March 23, with the first epoch beginning 12.3 hr after the GW detection. The survey covered the highest-probability sky area of 56.6 deg2, corresponding to a 91% probability. This was the first deep follow-up (m r ≳ 24, m z ≳ 23) for a binary-black hole merger covering >90% of the localization. By performing image subtraction and candidate screening including light-curve fitting with transient templates and examples, we found 22 off-nucleus transients that were not ruled out as the counterparts of GW200224_222234 with our Subaru/HSC data alone. We also performed GTC/OSIRIS spectroscopy of the probable host galaxies for five candidates; two are likely to be located within the 3D skymap, whereas the others are not. In conclusion, 19 transients remain as possible optical counterparts of GW200224_222234; but we could not identify a unique promising counterpart. If there are no counterparts in the remaining candidates, the upper limits of the optical luminosity are erg s−1 and erg s−1 in the r2 and z bands, respectively, at ∼12 hr after GW detection. We also discuss improvements in the strategies of optical follow-ups for future GW events.

Nonlinear Effects in Black Hole Ringdown.

We report evidence for nonlinear modes in the ringdown stage of the gravitational waveform produced by the merger of two comparable-mass black holes. We consider both the coalescence of black hole binaries in quasicircular orbits and high-energy, head-on black hole collisions. The presence of nonlinear modes in the numerical simulations confirms that general-relativistic nonlinearities are important and must be considered in gravitational-wave data analysis.

Generation and propagation of nonlinear quasinormal modes of a Schwarzschild black hole

In the analysis of a binary black hole coalescence, it is necessary to include gravitational self-interactions in order to describe the transition of the gravitational wave signal from the merger to the ringdown stage. In this paper we study the phenomenology of the generation and propagation of nonlinearities in the ringdown of a Schwarzschild black hole, using second-order perturbation theory. Following earlier work, we show that the Green's function and its causal structure determines how both first-order and second-order perturbations are generated, and hence highlight that both of these solutions share some physical properties. In particular, we discuss the sense in which both linear and quadratic quasi-normal modes (QNMs) are generated in the vicinity of the peak of the gravitational potential barrier (loosely referred to as the light ring). Among the second-order perturbations, there are solutions with linear QNM frequencies (whose amplitudes are thus renormalized from their linear values), as well as quadratic QNM frequencies with a distinct spectrum. Moreover, we show using a WKB analysis that, in the eikonal limit, waves generated inside the light ring propagate towards the black hole horizon, and only waves generated outside propagate towards an asymptotic observer. These results might be relevant for recent discussions on the validity of perturbation theory close to the merger. Finally, we argue that even if nonlinearities are small, quadratic QNMs may be detectable and would likely be useful for improving ringdown models of higher angular harmonics and future tests of gravity.

Upgraded waveform model of eccentric binary black hole based on effective-one-body-numerical-relativity for spin-aligned binary black holes

Effective one body numerical relativity waveform models for spin-aligned binary black holes (SEOBNR) are based on the effective one body theoretical framework and numerical relativity simulation results. SEOBNR models have evolved through version 1 to version 4. We recently extended SEOBNRv1 model to SEOBNRE (Effective One Body Numerical Relativity waveform models for Spin-aligned binary black holes along Eccentric orbit) model which is also valid for spin-aligned binary black hole coalescence along eccentric orbit. In this paper, we update our previous SEOBNRE model to make it consistent to SEOBNRv4 which is the most widely used SEOBNR waveform model. This upgraded SEOBNRE model improves accuracy compared to previous SEOBNRE model, especially for highly spinning black holes. For spin-aligned binary black holes with mass ratio [Formula: see text], dimensionless spin [Formula: see text] and orbital eccentricity [Formula: see text] at reference frequency [Formula: see text] ([Formula: see text] is the total mass of the binary black hole, [Formula: see text]Hz), the upgraded SEOBNRE model can always fit numerical relativity waveform better than 98.2%. For most cases, the fitting factor can even be better than 99%.

Merging black holes in dwarf galaxies: calculating binary black hole coalescence time-scales from simulations forLISAdetection

ABSTRACTSupermassive black holes (SMBHs) merging in dwarf galaxies will be detectable by the Laser Interferometer Space Antenna (LISA) in the mid-2030s. Previous cosmological hydrodynamic simulations have shown the prediction of massive BHs merging in dwarf galaxies, but these simulations are limited by their resolution and cannot follow BH pairs all the way to coalescence. We calculate the delay time between BH pairing and merger based on the properties of the BHs and their host galaxies, and use these properties to calculate gravitational wave strains for eleven different binary BHs that merge inside dwarf galaxies from eight cosmological simulations. This delay time calculation accounts for dynamical friction due to gas and stars, loss-cone scattering, and hardening of the binary due to gravitational radiation. Out of the eleven BH mergers in the simulations, five BH pairs will merge within 0.8–8 Gyr of forming a close pair and could be observed by LISA, and the remaining six are unresolved due to resolution limitations of the simulation. As all five of the resolved close pairs merge within a Hubble time, we make the broad estimate that close SMBH pairs in dwarf galaxies will merge and be detectable by LISA, but this estimate depends on either the presence of gas during orbital decay or a solution to the dynamical buoyancy problem in cored potentials.

Optimizing the placement of numerical relativity simulations using a mismatch predicting neural network

Gravitational wave observations from merging compact objects are becoming commonplace, and as detectors improve and gravitational wave sources become more varied, it is increasingly important to have dense and expansive template banks of predicted gravitational waveforms. Since numerical relativity is the only way to fully solve the nonlinear merger regime of general relativity for comparably massed systems, numerical relativity simulations are critical for gravitational wave detection and analysis. These simulations are computationally expensive, with each simulation placing one point within the high dimensional parameter space of binary black hole coalescences. This makes it important to have a method of placing new simulations in ways that use our computational resources optimally while ensuring sufficient coverage of the parameter space. Accomplishing this requires predicting the impact of a new set of parameters before performing the simulation. To this effect, this paper introduces a neural network to predict the mismatch between the gravitational waves of two binary systems. Using this network, we then show how we can propose new numerical relativity simulations that will provide the most benefit. We also use the network to identify gaps in existing public catalogs and identify degeneracies in the binary black hole parameter space.

Novel Ringdown Amplitude-Phase Consistency Test.

The ringdown signal emitted during a binary black hole coalescence can be modeled as a linear superposition of the characteristic damped modes of the remnant black hole that get excited during the merger phase. While checking the consistency of the measured frequencies and damping times against theKerr BH spectrum predicted by general relativity (GR) is a cornerstone of strong-field tests of gravity, the consistency of measured excitation amplitudes and phases have been largely left unexplored. For a nonprecessing, quasicircular binary black hole merger, we find that GR predicts a narrow region in the space of mode amplitude ratio and phase difference, independently of the spin of the binary components. Using this unexpected result, we develop a new null test of strong-field gravity which demands that the measured amplitudes and phases of different ringdown modes should lie within this narrow region predicted by GR. We call this the amplitude-phase consistency test and introduce a procedure for performing it using information from the ringdown signal. Lastly, we apply this test to the GW190521 event, using the multimodal ringdown parameters inferred by Capano etal. [arXiv:2105.05238]. While ringdown measurements errors for this event are large, we show that GW190521 is consistent with the amplitude-phase consistency test. Our test is particularly well suited for accommodating multiple loud ringdown detections as those expected in the near future, and can be used complementarily to standard black-hole spectroscopy as a proxy for modified gravity, compact objects other than black holes, binary precession and eccentricity.

Eccentricity or spin precession? Distinguishing subdominant effects in gravitational-wave data

ABSTRACT Eccentricity and spin precession are key observables in gravitational-wave astronomy, encoding precious information about the astrophysical formation of compact binaries together with fine details of the relativistic two-body problem. However, the two effects can mimic each other in the emitted signals, raising issues around their distinguishability. Since inferring the existence of both eccentricity and spin precession simultaneously is – at present – not possible, current state-of-the-art analyses assume that either one of the effects may be present in the data. In such a setup, what are the conditions required for a confident identification of either effect? We present simulated parameter inference studies in realistic LIGO/Virgo noise, studying events consistent with either spin precessing or eccentric binary black hole coalescences and recovering under the assumption that either of the two effects may be at play. We quantify how the distinguishability of eccentricity and spin precession increases with the number of visible orbital cycles, confirming that the signal must be sufficiently long for the two effects to be separable. The threshold depends on the injected source, with inclination, eccentricity, and effective spin playing crucial roles. In particular, for injections similar to GW190521, we find that it is impossible to confidently distinguish eccentricity from spin precession.

Pseudospectrum and binary black hole merger transients

The merger phase of binary black hole coalescences is a transient between an initial oscillating regime (inspiral) and a late exponentially damped phase (ringdown). In spite of the non-linear character of Einstein equations, the merger dynamics presents a surprisingly simple behaviour consistent with effective linearity. On the other hand, energy loss through the event horizon and by scattering to infinity renders the system non-conservative. Hence, the infinitesimal generator of the (effective) linear dynamics is a non-selfadjoint operator. Qualitative features of transients in linear dynamics driven by non-selfadjoint (in general, non-normal) operators are captured by the pseudospectrum of the time generator. We propose the pseudospectrum as a unifying framework to thread together the phases of binary black hole coalescences, from the inspiral-merger transition up to the late quasinormal mode ringdown.