Massive Black Hole Mergers Research Articles

Detection and prediction of future massive black hole mergers with machine learning and truncated waveforms

Detection and prediction of future massive black hole mergers with machine learning and truncated waveforms

An effective model for the tidal disruption of satellites undergoing minor mergers with axisymmetric primaries

According to the hierarchical formation paradigm, galaxies form through mergers of smaller entities and massive black holes (MBHs), if present at their centers, migrate to the nucleus of the newly formed galaxy, where they form binary systems. The formation and evolution of MBH binaries, and in particular their coalescence timescale, is highly relevant for current and future facilities aimed at detecting the gravitational wave signal produced by the MBHs close to coalescence. While most of the studies targeting this process are based on hydrodynamic simulations, the high computational cost makes a complete parameter space exploration prohibitive. Semianalytic approaches represent a valid alternative, but they require ad hoc prescriptions for the mass loss of the merging galaxies in minor mergers due to tidal stripping, which is not commonly considered or is at best modelled assuming very idealised geometries. In this work we propose a novel, effective model for the tidal stripping in axisymmetric potentials, to be implemented in semi-analytic models. We validated our semi-analytic approach against N-body simulations considering different galaxy sizes, inclinations, and eccentricities, finding only a moderate dependence on the orbit eccentricity. In particular, we find that, for almost circular orbits, our model mildly overestimates the mass loss, and this is due to the adjustment of the stellar distribution after the mass is removed. Nonetheless, the model exhibits a very good agreement with simulations in all the considered conditions, and thus represents an extremely powerful addition to semi-analytic calculations.

In search of the biggest bangs since the big bang

Many galaxies contain SuperMassive Black Holes (SMBHs), whose formation and history raise many puzzles. Pulsar timing arrays have recently discovered a low-frequency cosmological “hum” of gravitational waves that may be emitted by SMBH binary systems, and the JWST and other telescopes have discovered an unexpectedly large population of high-redshift SMBHs. We argue that these two discoveries may be linked, and that they may enhance the prospects for measuring gravitational waves emitted during the mergers of massive black holes, thereby opening the way towards resolving many puzzles about SMBHs as well as providing new opportunities to probe general relativity.

Forming off-center massive black hole binaries in dwarf galaxies through Jacobi capture

It is well established that black holes reside in the central regions of virtually all types of known galaxies. Recent observational and numerical studies however challenge this picture, suggesting that intermediate-mass black holes in dwarf galaxies may be found on orbits far from the center. In particular, constant-density cores minimize orbital energy losses due to dynamical friction, and allow black holes to settle on stable off-center orbits. Using controlled simulations, we study the dynamics of off-center black holes in dwarf galaxies with such cores. We propose a new scenario to describe off-center mergers of massive black holes, starting with a Jacobi capture. We focus on initially circular and co-planar black hole orbits and explore a large parameter space of black hole masses and orbital parameters. We find that Jacobi captures are a complex and chaotic phenomenon that occurs in about 13% of cases in this simplified setup, and we quantify how the likelihood of capture depends on the simulation parameters. We note that this percentage is likely an upper limit of the general case. Nevertheless, we show that Jacobi captures in cored dwarf galaxies can facilitate the formation of off-center black hole binaries, and that this process is sufficiently common to have a substantial effect on black hole growth. While our setup only allows for temporary captures, we expect dissipative forces from baryons and post-Newtonian corrections to maintain the captures over time and to lead to the formation of stable binary systems. This motivates future studies of the effectiveness of such dissipative forces, within stripped nuclei or globular clusters, in forming stable bound systems.

Connecting low-redshift LISA massive black hole mergers to the nHz stochastic gravitational wave background

Pulsar Timing Array (PTA) experiments worldwide recently reported evidence of a nHz stochastic gravitational wave background (sGWB) compatible with the existence of slowly inspiralling massive black hole (MBH) binaries (MBHBs). The shape of the signal contains valuable information about the evolution of z < 1 MBHs above 108 M⊙, suggesting a faster dynamical evolution of MBHBs towards the gravitational-wave-driven inspiral or a larger MBH growth than usually assumed. In this work, we investigate if the nHz sGWB could also provide constraints on the population of merging lower-mass MBHBs (< 107 M⊙) detectable by LISA. To this end, we use the L-Galaxies semi-analytical model applied to the Millennium suite of simulations. We generate a population of MBHs compatible simultaneously with current electromagnetic and nHz sGWB constraints by including the possibility that, in favourable environments, MBHs can accrete gas beyond the Eddington limit. The predictions of this new model for the sGWB show that the global (integrated up to high-z) LISA detection rate is not significantly affected when compared to a fiducial model whose nHz sGWB signal is ∼2 times smaller. In both cases, the global rate yields ∼12 yr−1 and is dominated by systems of 105 − 6 M⊙. The main differences are limited to low-z (z < 3), high-mass (> 106 M⊙) LISA MBHBs. The model compatible with the latest PTA results predicts up to ∼1.6 times more detections, with a rate of ∼1 yr−1. We find that these LISA MBHB systems have 50% probability of shining with bolometric luminosities > 1043 erg s−1. Hence, in case PTA results are confirmed and given the current MBH modelling, our findings suggest there will be higher chances to perform multimessenger studies with LISA MBHB than previously expected.

GRMHD simulations of accretion flows onto massive binary black hole mergers embedded in a thin slab of gas

GRMHD simulations of accretion flows onto massive binary black hole mergers embedded in a thin slab of gas

A multiresolution method for modelling galaxy and massive black hole mergers

ABSTRACT The coalescence of the most massive black hole (MBH) binaries releases gravitational waves (GWs) within the detectable frequency range of pulsar timing arrays (PTAs; 10−9 to 10−6 Hz). The incoherent superposition of GWs from MBH mergers, the stochastic gravitational wave background (GWB), can provide unique information on MBH parameters and the large-scale structure of the Universe. The recent evidence for a GWB reported by the PTAs opens an exciting new window on to MBHs and their host galaxies. However, the astrophysical interpretation of the GWB requires accurate estimations of MBH merger time-scales for a statistically representative sample of galaxy mergers. This is numerically challenging; a high numerical resolution is required to avoid spurious relaxation and stochastic effects, while a large number of simulations are needed to sample a cosmologically representative volume. Here, we present a new multimass modelling method to increase the central resolution of a galaxy model at a fixed particle number. We follow mergers of galaxies hosting central MBHs with the fast multiple method code griffin at two reference resolutions and with two refinement schemes. We show that both refinement schemes are effective at increasing central resolution, reducing spurious relaxation and stochastic effects. A particle number of N ≥ 106 within a radius of five times the sphere of influence of the MBHs is required to reduce numerical scatter in the binary eccentricity and the coalescence time-scale to <30 per cent, a resolution that can only be reached at present with the mass refinement scheme.

GRMHD simulations of accretion flows onto unequal-mass, precessing massive binary black hole mergers

GRMHD simulations of accretion flows onto unequal-mass, precessing massive binary black hole mergers

Massive black hole binaries as sources of low-frequency gravitational waves and X-shaped radio galaxies

ABSTRACT We present the study of multimessenger signatures of massive black hole (MBH) binaries residing in the centres of galaxy merger remnants. In particular, we first focus on the gravitational wave background (GWB) produced by an ensemble of MBH binary inspirals in the frequency range probed by the Pulsar Timing Array (PTA) experiments. The improved estimates of the characteristic strain were obtained with the inclusion of environmental effects on the MBH binary orbital decay within the galaxy merger remnants, added in post-processing to the semi-analytical model of galaxy formation and evolution SHARK. Secondly, we explore two, intriguing in terms of the MBH binary evolution studies, hypotheses aiming to explain the origins of X-shaped radio galaxies – a peculiar type of objects with double lobe structures, constituting approximately 6–10 per cent of known radio loud galaxies. The two considered scenarios involve either an abrupt change in the jet direction after an MBH merger (a spin-flip) or an unresolved close binary, where each of the two components produces a jet. We find that the estimated GWB amplitude at the reference frequency $f_0=1 \, {\rm yr}^{-1}$ is in the range of $A_{\rm { yr^{-1}}} = 1.20\times 10^{-15}{\!-\!}1.46\times 10^{-15}$, which is 50 per cent lower than the strain of the signal detected by the PTA experiments. We also show that the spin-flip scenario considered in gas-poor mergers reproduces the observed properties of X-shaped radio galaxies well in terms of flip angle, redshift, and luminosity distributions.

Signatures of Massive Black Hole Merger Host Galaxies from Cosmological Simulations. I. Unique Galaxy Morphologies in Imaging

Low-frequency gravitational-wave experiments such as the Laser Interferometer Space Antenna and pulsar timing arrays are expected to detect individual massive black hole (MBH) binaries and mergers. However, secure methods of identifying the exact host galaxy of each MBH merger among the large number of galaxies in the gravitational-wave localization region are currently lacking. We investigate the distinct morphological signatures of MBH merger host galaxies, using the Romulus25 cosmological simulation. We produce mock telescope images of 201 simulated galaxies in Romulus25 hosting recent MBH mergers through stellar population synthesis and dust radiative transfer. Based on comparisons to mass- and redshift-matched control samples, we show that combining multiple morphological statistics via a linear discriminant analysis enables identification of the host galaxies of MBH mergers, with accuracies that increase with chirp mass and mass ratio. For mergers with high chirp masses (≳108.2 M ⊙) and high mass ratios (≳0.5), the accuracy of this approach reaches ≳80%, and does not decline for at least ∼1 Gyr after numerical merger. We argue that these trends arise because the most distinctive morphological characteristics of MBH merger and binary host galaxies are prominent classical bulges, rather than relatively short-lived morphological disturbances from their preceding galaxy mergers. Since these bulges are formed though major mergers of massive galaxies, they lead to (and become permanent signposts for) MBH binaries and mergers that have high chirp masses and mass ratios. Our results suggest that galaxy morphology can aid in identifying the host galaxies of future MBH binaries and mergers.

Adapting the PyCBC pipeline to find and infer the properties of gravitational waves from massive black hole binaries in LISA

The laser interferometer space antenna (LISA), due for launch in the mid 2030s, is expected to observe gravitational waves (GWs) from merging massive black hole binaries (MBHBs). These signals can last from days to months, depending on the masses of the black holes, and are expected to be observed with high signal to noise ratios (SNRs) out to high redshifts. We have adapted the PyCBC software package to enable a template bank search and inference of GWs from MBHBs. The pipeline is tested on the LISA data challenge’s Challenge 2a (‘Sangria’), which contains MBHBs and thousands of galactic binaries (GBs) in simulated instrumental LISA noise. Our search identifies all six MBHB signals with more than 92% of the optimal SNR. The subsequent parameter inference step recovers the masses and spins within their 90% confidence interval. Sky position parameters have eight high likelihood modes which are recovered but often our posteriors favour the incorrect sky mode. We observe that the addition of GBs biases the parameter recovery of masses and spins away from the injected values, reinforcing the need for a global fit pipeline which will simultaneously fit the parameters of the GB signals before estimating the parameters of MBHBs.

Implications of the pulsar timing array detections for massive black hole mergers in the LISA band

Implications of the pulsar timing array detections for massive black hole mergers in the LISA band

Too small to fail: characterizing sub-solar mass black hole mergers with gravitational waves

The detection of a sub-solar mass black hole could yield dramatic new insights into the nature of dark matter and early-Universe physics, as such objects lack a traditional astrophysical formation mechanism. Gravitational waves allow for the direct measurement of compact object masses during binary mergers, and we expect the gravitational-wave signal from a low-mass coalescence to remain within the LIGO frequency band for thousands of seconds. However, it is unclear whether one can confidently measure the properties of a sub-solar mass compact object and distinguish between a sub-solar mass black hole or other exotic objects. To this end, we perform Bayesian parameter estimation on simulated gravitational-wave signals from sub-solar mass black hole mergers to explore the measurability of their source properties. We find that the LIGO/Virgo detectors during the O4 observing run would be able to confidently identify sub-solar component masses at the threshold of detectability; these events would also be well-localized on the sky and may reveal some information on their binary spin geometry. Further, next-generation detectors such as Cosmic Explorer and the Einstein Telescope will allow for precision measurement of the properties of sub-solar mass mergers and tighter constraints on their compact-object nature.

Supernova calibration by gravitational waves

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

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

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

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

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

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

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

Multi-messenger study of merging massive black holes in the OBELISK simulation: Gravitational waves, electromagnetic counterparts, and their link to galaxy and black-hole populations

Massive black-hole (BH) mergers are predicted to be powerful sources of low-frequency gravitational waves (GWs). Coupling the detection of GWs with an electromagnetic (EM) detection can provide key information about merging BHs and their environments as well as cosmology. We study the high-resolution cosmological radiation-hydrodynamics simulation OBELISK, run to redshift z = 3.5, to assess the GW and EM detectability of high-redshift BH mergers, modelling spectral energy distribution and obscuration. For EM detectability, we further consider sub-grid dynamical delays in postprocessing. We find that most of the merger events can be detected by LISA, except for high-mass mergers with very unequal mass ratios. Intrinsic binary parameters are accurately measured, but the sky localisation is poor generally. Only ∼40% of these high-redshift sources have a sky localisation better than 10 deg2. Merging BHs are hard to detect in the restframe UV since they are fainter than the host galaxies, which at high redshift are star-forming. A significant fraction, 15–35%, of BH mergers instead outshine the galaxy in X-rays, and about 5 − 15% are sufficiently bright to be detected with sensitive X-ray instruments. If mergers induce an Eddington-limited brightening, up to 30% of sources can become observable. The transient flux change originating from such a brightening is often large, allowing 4 − 20% of mergers to be detected as EM counterparts. A fraction, 1 − 30%, of mergers are also detectable at radio frequencies. Transients are found to be weaker for radio-observable mergers. Observable merging BHs tend to have higher accretion rates and masses and are overmassive at a fixed galaxy mass with respect to the full population. Most EM-observable mergers can also be GW-detected with LISA, but their sky localisation is generally poorer. This has to be considered when using EM counterparts to obtain information about the properties of merging BHs and their environment.

Black hole mergers as tracers of spinning massive black hole and galaxy populations in the OBELISK simulation

Massive black hole (BH) mergers will be key targets of future gravitational wave and electromagnetic observational facilities. In order to constrain BH evolution with the information extracted from BH mergers, one must take into account the complex relationship between the population of merging BHs and the global BH population. We analysed the high-resolution cosmological radiation-hydrodynamics simulation OBELISK, run to redshift z = 3.5, to study the properties of the merging BH population, and its differences with the underlying global BH population in terms of BH and galaxy properties. In post-processing, we calculated dynamical delays between the merger in the simulation at the resolution limit and the actual coalescence well below the resolution scale. We find that merging BHs are hosted in relatively massive galaxies with stellar mass M* ≳ 109 M⊙. Given that galaxy mass is correlated with other BH and galaxy properties, BH mergers tend to also have a higher total BH mass and higher BH accretion rates than the global population of main BHs. These differences generally disappear if the merger population is compared with a BH population sampled with the same galaxy mass distribution as merger hosts. Galaxy mergers can temporarily boost the BH accretion rate and the host’s star formation rate, which can remain active at the BH merger if sub-resolution delays are not taken into account. When dynamical delays are taken into account, the burst has generally faded by the time the BHs merge. BH spins are followed self-consistently in the simulation under the effect of accretion and BH mergers. We find that merging BHs have higher spins than the global population, but similar or somewhat lower spins compared to a mass-matched sample. For our sample, mergers tend to decrease the spin of the final BH remnant.

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.