---
title: Intermediate-Mass Binary Black Hole Mergers
url: https://www.emergentmind.com/topics/intermediate-mass-binary-black-hole-mergers
type: topic
---

# Intermediate-Mass Binary Black Hole Mergers

Searching arXiv for recent papers on intermediate-mass binary black hole mergers and related channels.
Intermediate-mass binary black hole mergers occupy the regime between the stellar-mass binaries now routinely observed by LIGO–Virgo–KAGRA and the massive black hole binaries targeted by LISA. In current usage, the term covers coalescences in which at least one component lies in the intermediate-mass black hole (IMBH) range \(10^2 M_\odot \lesssim M \lesssim 10^6 M_\odot\), but the observational literature also includes heavy systems in the \(100\!-\!450\,M_\odot\) total-mass interval and binaries that produce an IMBH remnant, such as GW190521 with source-frame component masses \(85^{+21}_{-14}\,M_\odot\) and \(66^{+17}_{-18}\,M_\odot\) and a final remnant of \(142^{+28}_{-16}\,M_\odot\) [2406.01687] [1201.5999] [2010.06161].

## 1. Definition, mass scales, and observational scope

The standard mass scale adopted for IMBHs is \(10^2 M_\odot \lesssim M \lesssim 10^6 M_\odot\), while ground-based interferometers are sensitive mainly to the lower end of that interval, from a few tens up to a few thousand solar masses [2406.01687]. Within this broader category, the literature distinguishes at least three overlapping observational classes: comparable-mass IMBH binaries, IMBH–stellar-mass systems in the intermediate mass-ratio inspiral range \(10^{-4}\lesssim q\lesssim 0.1\), and heavy binaries whose components occupy or approach the pair-instability gap and whose merger remnant is itself an IMBH [2507.01083] [2010.06161].

GW190521 sharpened this taxonomy because its component masses lie well inside the “upper black-hole mass gap,” for which first-generation black holes formed by ordinary single-star evolution at low metallicity are expected to have a maximum mass of order \(45\text{–}60\,M_\odot\), with a quoted lower edge near \(\sim 50^{+4}_{-8}\,M_\odot\) [2010.06161]. This made hierarchical assembly a central interpretive theme for intermediate-mass mergers. A plausible implication is that the empirical category “intermediate-mass binary black hole merger” now includes both binaries containing an IMBH and binaries whose observed masses require repeated merger growth across the pair-instability gap.

Search strategies have reflected this mixed observational scope. Initial LIGO–Virgo targeted non-spinning mergers with total mass \(100\!-\!450\,M_\odot\) and mass ratios from \(1{:}1\) to \(4{:}1\) [1201.5999], while more recent LVK searches have explicitly targeted IMRIs with detector-frame masses \(54\,M_\odot < m_1 < 400\,M_\odot\), \(3\,M_\odot < m_2 < q\,m_1\), and \(1/100<q<1/18\) [2507.01083]. This observational spread is one reason the subject links stellar-dynamical, stellar-evolutionary, and relativistic-binary populations rather than a single sharply bounded mass class.

## 2. Formation channels across dense stellar and gaseous environments

A major channel is hierarchical growth in dense stellar systems. In direct \(N\)-body globular-cluster simulations with \texttt{NBODY6-gpu}, one black hole in a representative \(N_0=100{,}000\) cluster undergoes a chain of seven BH–BH mergers within \(6\) Gyr and reaches \(97.8\,M_\odot\), while the broader suite of 58 simulations supports an order-of-magnitude cluster-level probability \(0.01 \lesssim p \lesssim 0.1\) for IMBH formation through hierarchical mergers over a Hubble time in a median-mass globular cluster if gravitational recoil is neglected [2010.06161]. In this picture, first-generation black holes with \(M\lesssim 50\,M_\odot\) mass-segregate into the core, dynamically assemble binaries, and then build second- and higher-generation remnants through repeated merger cycles.

Nuclear star clusters provide a related but deeper-potential variant. Semi-analytic Monte Carlo calculations for nuclear star clusters find IMBH-binary merger rates in the range \(0.01\!-\!10\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}\), depending on IMBH mass, metallicity, seed mass, and initial spin [2204.03745]. Direct \(N\)-body calculations in nucleated dwarf galaxies go further and show that comparable-mass IMBHs merge very efficiently once brought inside the dense nuclear region: binary coalescence times are typically a few hundred million years, hardening rates span \(\sim 1.8\!-\!299\,{\rm pc^{-1}\,Myr^{-1}}\), and the stellar encounters that drive the binary originate predominantly from nuclear star clusters when the NSC-to-IMBH-binary mass ratio is greater than \(10\) [2109.12129].

Other channels populate the same mass regime with different initial conditions. Collisions and mergers of massive globular clusters in a Milky-Way-like disc can pair pre-existing IMBHs; for a globular-cluster collision rate of \(\sim 1.8\) per Gyr and an IMBH occupation fraction \(f_{\rm IMBH}\sim 10\%-20\%\), the expected number of IMBH–IMBH binaries formed per Milky-Way-like galaxy over a Hubble time is \(2\!-\!4\) [1906.05864]. In very massive Population III binaries, IMBH-involving mergers depend sensitively on the common-envelope parameter \(\alpha\lambda\); the study finds that the Einstein Telescope has the potential to detect \(\sim 10\!-\!1000\) BBHs with IMBHs per year [2211.07496]. In dense young clusters with collisional runaways, a very massive star of mass \(m_{\rm run}=f_{\rm run}M_{\rm CL}\) collapses into an IMBH seed, and the subsequent IMBH–BH mergers are expected to be broadly detectable, with some sources observable in multiple bands [2401.06731]. In AGN disks, IMBBHs with total mass \(M\sim 10^2\!-\!10^4\,M_\odot\) are treated as a distinct environmental population in which gas dynamical effects become directly measurable in the waveform [2503.11721].

## 3. Dynamical assembly, recoil, and distinctive waveform signatures

The recurrent dynamical ingredients are mass segregation, three-body binary formation, exchange interactions, binary hardening, and eventual gravitational-wave inspiral. In the globular-cluster hierarchical-merger picture, binaries form dynamically in the core, harden through binary–single and binary–binary encounters, and merge once GW emission dominates at small separation; the heaviest black hole is almost always retained in a binary through exchange bias, which reinforces hierarchical growth [2010.06161]. This same combination of Heggie’s law, exchange preference for the most massive component, and GW-driven inspiral underlies IMBH growth in several cluster channels.

Gravitational recoil is the principal brake on that hierarchy. In the seven-merger “Snowball” chain, post-processing with isotropic spin orientations shows that the second-generation remnant is ejected in \(\approx 85.2\%\) of realizations after the first merger, only \(\sim 0.719\%\) of chains reach length \(L=2\) across the full sample with spin recoil included, \(\sim 0.0112\%\) reach \(L=3\), \(\sim 3.06\times 10^{-4}\%\) reach \(L=4\), and no \(L\ge 5\) chains are seen [2010.06161]. This establishes recoil and spin orientation as gatekeepers of intermediate-mass growth in globular clusters. A plausible implication is that environments with high escape speeds are not merely favorable but often necessary for long hierarchical chains.

Several channels predict substantial eccentricity. In nucleated dwarf galaxies, IMBH binaries attain high eccentricities \(e\sim 0.4\!-\!0.99\), and residual eccentricity may persist to the LISA regime [2109.12129]. In AGN disks, the distinctive feature is not eccentricity but a measurable environmental phase correction from gas dynamical friction. The waveform phase acquires an additional term
\[
\phi_{\rm DF}(f) = -\rho \,\frac{25\pi (3\eta -1) M_c^2}{739328\, \eta^2} \,\gamma_{\rm DF} \left[\pi M_c (1+z) f\right]^{-16/3},
\]
which behaves as an effective \(-5.5\)PN correction and is therefore strongest at low frequencies [2503.11721]. For \(\rho \gtrsim 10^{-12}\,\mathrm{g/cc}\), LISA can detect this imprint; for \(\rho=10^{-12}\,\mathrm{g/cc}\), the density is measurable with order-unity fractional error, and for \(\rho=10^{-10}\!-\!10^{-8}\,\mathrm{g/cc}\) the constraint improves to a few percent [2503.11721]. This provides a direct route to distinguishing AGN-disk IMBBHs from effectively vacuum channels.

## 4. Searches with initial and current ground-based detectors

The first dedicated LIGO–Virgo search for IMBH mergers used a weakly modeled burst strategy on 2005–2007 data, targeting non-spinning systems with total mass \(100\!-\!450\,M_\odot\) and component mass ratios between \(1{:}1\) and \(4{:}1\) [1201.5999]. No plausible signals were observed. The most stringent 90% confidence rate-density upper limit was obtained in the bin centered on \(88+88\,M_\odot\), where the limit is \(0.13\) per \(\mathrm{Mpc}^3\) per Myr for non-spinning sources [1201.5999]. That search established the burst-based methodology for merger–ringdown-dominated signals but did not yet reach the theoretically favored astrophysical rate range.

Matched-filter searches improved this picture in O3. An optimized PyCBC search tailored to detector-frame total masses \(100\!-\!600\,M_\odot\), stricter signal-noise discriminators, and a \(15\) Hz low-frequency cutoff improves the sensitive volume-time product over previous PyCBC-based searches by a factor of \(1.5\) to \(3\), depending on total binary mass, at an inverse false alarm rate of \(100\) years [2106.00193]. Applied to O3a data, it identifies no new significant IMBH binaries but confirms GW190521 with a false alarm rate of \(1\) in \(727\) years [2106.00193]. In parallel, the hierarchical-merger interpretation of GW190521 remains quantitatively plausible in optimistic recoil conditions: a no-recoil globular-cluster estimate yields a volumetric rate \(R_{\rm Snowball}\approx 0.03\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}\), broadly consistent with the inferred GW190521-like rate \(0.13^{+0.30}_{-0.11}\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}\), whereas isotropic spin recoil suppresses the predicted rate to \(3.5\times10^{-12}\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}\) in the same median-mass cluster model [2010.06161].

The dedicated O3 IMRI search extends the observational program to the highly asymmetric sector \(1/100<q<1/18\), using aligned-spin template banks built from SEOBNRv5HM and including higher modes for the first time [2507.01083]. The inclusion of higher modes increases the sensitivity volume by up to \(\sim 500\%\). No significant candidates with IFAR \(>1\) year are found, and the resulting upper limits on the local IMRI merger rate lie between \(\sim 30\) and \(10^3\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}\), depending on the component masses [2507.01083]. This places present LVK searches in a regime where sensitivity to the IMBH–stellar-mass corner is improving rapidly but still remains above many theoretical rate estimates.

## 5. Next-generation detectors, multiband observations, and parameter estimation

Forecasts for next-generation ground-based networks show that IMBBHs should become precision sources provided low-frequency sensitivity reaches a few hertz. For the CE40–CE20–ET network with \(f_{\rm min}=3\) Hz, binaries with component masses \(m_{1,2}\sim 1000\,M_\odot\) have mass errors \(\lesssim 0.1\%\) at \(z=0.5\) and \(\lesssim 1\%\) at \(z=2\), while the source redshift can be measured with percent-level accuracy or better at \(z\lesssim 2\) [2406.01687]. Lighter binaries with \(m_{1,2}\lesssim 300\,M_\odot\) still admit \(O(10)\%\) redshift accuracy even at \(z=10\), and binaries with \(z\lesssim 0.5\) can be localized within \(1\,\mathrm{deg}^2\) for \(m_{1,2}\lesssim 1000\,M_\odot\), with many comparable-mass systems localized within \(0.1\,\mathrm{deg}^2\) [2406.01687]. The same study emphasizes that low-frequency sensitivity is crucial: raising \(f_{\rm min}\) from \(3\) Hz to \(7\!-\!10\) Hz drastically reduces detection horizons and parameter-estimation quality for \(M\sim 10^3\,M_\odot\) binaries [2406.01687].

Space-based observations add both reach and environmental diagnostics. In AGN disks, a single LISA observation of a source with total mass \(10^3\,M_\odot\), mass ratio \(2\), and luminosity distance \(3\) Gpc is sufficient to associate the merger with an AGN disk of density \(\sim 10^{-12}\,\mathrm{g/cc}\), because the density can be estimated with an error bar \(O(100\%)\); joint LISA+ET analyses improve density constraints by a factor \(\sim 2\!-\!5\), and for a \(1000\,M_\odot\), \(q=2\) system the minimum fractional density error reaches \(\sim 0.6\%\) in the examined mass–density range [2503.11721]. More generally, nuclear-star-cluster models predict that a few merger events per year should be detectable with LISA, DECIGO, ET, and LIGO for IMBHs with masses \(\lesssim 1000\,M_\odot\), and a few tens of merger events per year with DECIGO, ET, and LIGO only [2204.03745]. The collisional-runaway cluster models likewise find that most IMBH–BH merger signals could be detected, with some of them being multi-band sources [2401.06731].

These capabilities imply that future catalogs will not merely enlarge the sample; they will separate channels observationally. Mass-gap occupancy, redshift evolution, sky localization, residual eccentricity in the mHz band, higher-mode content in IMRIs, and direct environmental phasing in AGN disks all become measurable diagnostics once IMBBHs are seen across LISA, ET, CE, DECIGO, and improved ground-based networks [2406.01687] [2503.11721].

## 6. Theoretical uncertainties, host dependence, and population-level implications

The principal uncertainty across channels is retention. In globular clusters, recoil treatment, natal-kick prescriptions, the absence or inclusion of primordial binaries, and the lack of post-Newtonian terms in the direct \(N\)-body equations all affect whether hierarchical chains survive long enough to cross the pair-instability gap [2010.06161]. In Population III binaries, the maximum primary BH mass and the survival of IMBH-involving mergers to low redshift depend strongly on the common-envelope parameter \(\alpha\lambda\); for \(\alpha\lambda=0.01\), all IMBH-involving mergers have merged by \(z\sim 3\), while the Einstein Telescope may detect \(\sim 10\!-\!1000\) such systems per year across the allowed model set [2211.07496]. In AGN-disk systems, omitting environmental effects in waveform models biases the chirp mass, mass ratio, and arrival time, which directly affects multiband analyses [2503.11721].

Host structure also controls whether IMBH binaries merge efficiently at all. Nucleated dwarf galaxies produce very efficient pairing and coalescence on timescales of a few hundred million years [2109.12129], whereas in non-nucleated dwarf galaxies the direct \(N\)-body simulations reported in the lowest-density regime find that none of the IMBHs in the simulation suite merge within a Hubble time [2408.14541]. This suggests that the presence or absence of a dense nuclear star cluster is a first-order discriminator between efficient and stalled IMBH assembly in dwarfs.

Intermediate-mass black holes also feed back on the broader black-hole merger ecosystem. In the MOCCA survey, the number of stellar-mass BBH mergers from a globular cluster follows the expected scaling with initial mass and density unless the present-day IMBH mass is more massive than \(\sim 10^4\,M_\odot\) or exceeds about \(1\) per cent of the cluster’s initial total mass; in such cases, binary-black-hole formation and subsequent merger events are significantly reduced, and the under-population of BBHs can reach a factor of \(2\) depending on the clusters’ initial distributions [2008.10823]. A plausible implication is that a mature IMBH population can suppress one observational channel while enhancing another: fewer ordinary cluster BBH mergers, but more IMBH–BH or IMBH-seeded hierarchical events.

Taken together, the current literature presents intermediate-mass binary black hole mergers as a heterogeneous but physically connected population. Dense stellar systems build them through repeated exchange and hardening, gaseous disks imprint environmental phase shifts, dwarf nuclei can merge them efficiently or not at all depending on nuclear structure, and next-generation detectors will measure masses, redshifts, localizations, and in some cases environmental parameters precisely enough to separate these channels empirically [2010.06161] [2406.01687] [2204.03745].

Source: https://www.emergentmind.com/topics/intermediate-mass-binary-black-hole-mergers