Coalescence Of Supermassive Black Holes Research Articles

From 100 MHz to 10 GHz: Unveiling the spectral evolution of the X-shaped radio galaxy in Abell 3670

Context. X-shaped radio galaxies (XRGs) are characterised by two pairs of misaligned lobes: active lobes hosting radio jets and the wings. None of the formation mechanisms proposed thus far are able to exhaustively reproduce the diverse features observed among XRGs. Emerging evidence has proposed the existence of sub-populations of XRGs forming via different processes. Aims. The brightest cluster galaxy (BCG) in Abell 3670 (A3670) is a dumbbell system hosting the XRG MRC 2011-298. The morphological and spectral properties of this interesting XRG were first characterised based on Karl G. Jansky Very Large Array (JVLA) data at 1–10 GHz. In the present work, we follow up on MRC 2011-298 with the upgraded Giant Metrewave Radio Telescope (uGMRT) at 120–800 MHz to further constrain its properties and origin. Methods. We carried out a detailed spectral analysis sampling different spatial scales. Integrated radio spectra, spectral index maps, radio colour-colour diagrams, and radiative age maps of both the active lobes and prominent wings were employed to test the origin of the source. Results. We confirm a progressive spectral steepening from the lobes to the wings. The maximum radiative age of the source is ~80 Myr, with the wings being older than the lobes by ≳30 Myr in their outermost regions. Conclusions. The observed properties are in line with an abrupt reorientation of the jets by ~90 deg from the direction of the wings to their present position. This formation mechanism is further supported by the comparison with numerical simulations in the literature, which additionally highlight the role of hydrodynamic processes in the evolution of large wings such as those of MRC 2011-298. It is plausible that the coalescence of supermassive black holes could have triggered the spin-flip of the jets. Moreover, we show that the S-shape of the radio jets is likely driven by precession with a period of P ~ 10 Myr.

RABBITS – II. The impact of AGN feedback on coalescing supermassive black holes in disc and elliptical galaxy mergers

ABSTRACT In this study of the ‘Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations’ (RABBITS) series, we investigate the orbital evolution of supermassive black holes (SMBHs) during galaxy mergers. We simulate both disc and elliptical galaxy mergers using the ketju code, which can simultaneously follow galaxy (hydro-)dynamics and small-scale SMBH dynamics with post-Newtonian corrections. With our SMBH binary subgrid model, we show how active galactic nuclei (AGNs) feedback affects galaxy properties and SMBH coalescence. We find that simulations without AGN feedback exhibit excessive star formation, resulting in merger remnants that deviate from observed properties. Kinetic AGN feedback proves more effective than thermal AGN feedback in expelling gas from the centre and quenching star formation. The different central galaxy properties, which are a result of distinct AGN feedback models, lead to varying rates of SMBH orbital decay. In the dynamical friction phase, galaxies with higher star formation and higher SMBH masses possess denser centres, become more resistant to tidal stripping, experience greater dynamical friction, and consequently form SMBH binaries earlier. As AGN feedback reduces gas densities in the centres, dynamical friction by stars dominates over gas. In the SMBH hardening phase, compared to elliptical mergers, disc mergers exhibit higher central densities of newly formed stars, resulting in accelerated SMBH hardening and shorter merger time-scales (i.e. $\lesssim 500$ Myr versus $\gtrsim 1$ Gyr). Our findings highlight the importance of AGN feedback and its numerical implementation in understanding the SMBH coalescing process, a key focus for low-frequency gravitational wave observatories.

Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers

ABSTRACTWe introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the ketju code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi–Hoyle–Lyttleton accretion into the binary phase and incorporates preferential mass accretion on to the secondary SMBH, which is motivated by results from small-scale hydrodynamical circumbinary disc simulations. We perform idealized gas-rich disc galaxy merger simulations using pure thermal or pure kinetic active galactic nuclei (AGNs) feedback. Our binary accretion model provides more physically motivated SMBH mass ratios, which are one of the key parameters for computing gravitational wave (GW) induced recoil velocities. The merger time-scales of our simulated SMBH binaries are in the range tmerge ∼ 10–400 Myr. Prograde in-plane equal-mass galaxy mergers lead to the shortest merger time-scales, as they experience the strongest starbursts, with the ensuing high stellar density resulting in a rapid SMBH coalescence. Compared to the thermal AGN feedback, the kinetic AGN feedback predicts longer merger time-scales and results in more core-like stellar profiles, as it is more effective in removing gas from the galaxy centre and quenching star formation. This suggests that the AGN feedback implementation plays a critical role in modelling SMBH coalescences. Our model will be useful for improving the modelling of SMBH mergers in gas-rich galaxies, the prime targets for the upcoming LISA GW observatory.

The role of mergers and gas accretion in black hole growth and galaxy evolution

We use a semi-analytic galaxy formation model to study the co-evolution of supermassive black holes (SMBHs) with their host galaxies. Although the coalescence of SMBHs is not important, the quasar-mode accretion induced by mergers plays a dominant role in the growth of SMBHs. Mergers play a more important role in the growth of SMBH host galaxies than in the SMBH growth. It is the combined contribution from quasar mode accretion and mergers to the SMBH growth and the combined contribution from starburst and mergers to their host galaxy growth that determine the observed scaling relation between the SMBH masses and their host galaxy masses. We also find that mergers are more important in the growth of SMBH host galaxies compared to normal galaxies which share the same stellar mass range as the SMBH host galaxies.

Evolution of gas disc–embedded intermediate mass ratio inspirals in theLISAband

ABSTRACTAmong the potential milliHz gravitational wave (GW) sources for the upcoming space-based interferometer LISA are extreme- or intermediate-mass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range 105M⊙ ≲ M ≲ 107M⊙ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGNs), where torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform. Here we use a suite of 2D hydrodynamical simulations with the moving-mesh code disco to present a systematic study of disc torques. We measure torques on an inspiralling BH and compute the corresponding waveform deviations as a function of the binary mass ratio q ≡ M2/M1, the disc viscosity (α), and gas temperature (or equivalently Mach number; $\mathcal {M}$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends non-trivially on the combination of these parameters. The gas imprint is detectable by LISA for binaries embedded in AGN discs with surface densities above $\Sigma _0\ge 10^{4-6} \rm \, g cm^{-2}$, depending on q, α, and $\mathcal {M}$. Deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GW-driven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of disc–satellite interaction.

The effect of cooling on the accretion of circumprimary discs in merging supermassive black hole binaries

ABSTRACT At the final stages of a supermassive black hole coalescence, the emission of gravitational waves will efficiently remove energy, and angular momentum from the binary orbit, allowing the separation between the compact objects to shrink. In the scenario where a circumprimary disc is present, a squeezing phase will develop, in which the tidal interaction between the disc and the secondary black hole could push the gas inwards, enhancing the accretion rate on to the primary and producing what is known as an electromagnetic precursor. In this context, using 3D hydrodynamic simulations, we study how an adiabatic circumprimary accretion disc responds to the varying gravitational potential as the secondary falls on to the more massive object. We included a cooling prescription controlled by the parameter β = Ωtcool, which will determine how strong the final accretion rate is: a hotter disc is thicker, and the tidal interaction is suppressed for the gas outside the binary plane. Our main results are that for scenarios where the gas cannot cool fast enough (β ≥ 30), the disc becomes thick and renders the system invisible, while for β ≤ 10 the strong cooling blocks any leakage on to the secondary’s orbit, allowing an enhancement in the accretion rate of two orders of magnitude stronger than the average through the rest of the simulation.

Accelerated orbital decay of supermassive black hole binaries in merging nuclear star clusters

ABSTRACT The coalescence of supermassive black holes (SMBHs) should generate the strongest sources of gravitational waves (GWs) in the Universe. However, the dynamics of their coalescence is the subject of much debate. In this study, we use a suite of N-body simulations to follow the merger of two nuclear star clusters (NSCs), each hosting an SMBH in their centre. We find that the presence of distinct star clusters around each SMBH has important consequences for the dynamical evolution of the SMBH binary: (i) The separation between the SMBHs decreases by a few orders of magnitude in the first few Myrs by the combined effects of dynamical friction and a drag force caused by tidally stripped stars. In fact, this is a significant speedup for equal mass ratio binaries, and becomes extreme for unequal mass ratios, e.g. 1:10 or 1:100, which traditional dynamical friction alone would not permit to bind. (ii) The subsequent binary hardening is driven by the gravitational slingshots between the SMBH binary and stars, and also depends on the mass ratio between the SMBHs. Thus, with this additional drag force, we find that all SMBHs in our suite coalesce within a Hubble time. Given that about 50 per cent of Milky Way-sized galaxies host NSCs, our results are encouraging for upcoming GW observations with the Laser Interferometer Space Antenna – LISA – which will detect SMBH coalescence in the 104–107 M⊙ mass range.

Supermassive black holes coalescence mediated by massive perturbers: implications for gravitational waves emission and nuclear cluster formation

A large fraction of galactic nuclei is expected to host supermassive black hole binaries (BHB), likely formed during the early phase of galaxies assembly and merging. In this paper, we use a large set of state-of-art numerical models to investigate the interplay between a BHB and a massive star cluster (GCs) driven toward the galactic centre by dynamical friction. Varying the BHB mass and mass ratio and the GC orbit, we show that the reciprocal feedback exerted between GCs and the BHB shapes their global properties. We show that, at GC-to-BHB mass ratios above 0.1, the GC affects notably the BHB orbital evolution, possibly boosting its coalescence. This effect is maximized if the GC moves on a retrograde orbit, and for a non-equal mass BHB. We show that the GC debris dispersed around the galactic nucleus can lead to the formation of a nuclear cluster, depending on the BHB tidal field, and that the distribution of compact remnants resulting from the GC disruption can carry information about the BHB orbital properties. We find that red giant stars delivered by the spiralling GC can be disrupted at a rate of $\simeq (0.7-7)\times 10^{-7}$ yr$^{-1}$ for BHB masses $\sim 10^7{\rm M}_\odot$, while tens to hundreds of stars can be possibly observed in the galactic halo as high-velocity stars, with velocities up to $\sim 2000$ km s$^{-1}$, depending on the BHB orbital properties.

Post-Newtonian Dynamical Modeling of Supermassive Black Holes in Galactic-scale Simulations

Abstract We present KETJU, a new extension of the widely used smoothed particle hydrodynamics simulation code GADGET-3. The key feature of the code is the inclusion of algorithmically regularized regions around every supermassive black hole (SMBH). This allows for simultaneously following global galactic-scale dynamical and astrophysical processes, while solving the dynamics of SMBHs, SMBH binaries, and surrounding stellar systems at subparsec scales. The KETJU code includes post-Newtonian terms in the equations of motions of the SMBHs, which enables a new SMBH merger criterion based on the gravitational wave coalescence timescale, pushing the merger separation of SMBHs down to ∼0.005 pc. We test the performance of our code by comparison to NBODY7 and rVINE. We set up dynamically stable multicomponent merger progenitor galaxies to study the SMBH binary evolution during galaxy mergers. In our simulation sample the SMBH binaries do not suffer from the final-parsec problem, which we attribute to the nonspherical shape of the merger remnants. For bulge-only models, the hardening rate decreases with increasing resolution, whereas for models that in addition include massive dark matter halos, the SMBH binary hardening rate becomes practically independent of the mass resolution of the stellar bulge. The SMBHs coalesce on average 200 Myr after the formation of the SMBH binary. However, small differences in the initial SMBH binary eccentricities can result in large differences in the SMBH coalescence times. Finally, we discuss the future prospects of KETJU, which allows for a straightforward inclusion of gas physics in the simulations.

X-Shaped Radio Galaxies and the Nanohertz Gravitational Wave Background

AbstractCoalescence of supermassive black holes (SMBHs) in galaxy mergers is potentially the dominant contributor to the low frequency gravitational wave background (GWB). It was proposed by Merritt & Ekers that X-shaped radio galaxies are signposts of such coalescences and that their abundance might be used to predict the magnitude of the GWB. Cheung identified a sample of 100 candidate X-shaped radio galaxies using the NRAO FIRST survey; these are small-axial-ratio extended radio sources with off-axis emission. In Roberts et al. we made radio images of 52 of these sources with resolution of about 1 arcsecond using archival Very Large Array data. Fifty-one of the 52 were observed at 1.4 GHz, seven were observed at 1.4 and 5 GHz, and one was observed only at 5 GHz. Our higher resolution VLA images along with FIRST survey images of the sources in the sample reveal that extended extragalactic radio sources with small axial ratios are largely (60%) cases of double radio sources with twin lobes that have off-axis extensions, usually with inversion-symmetric structure. The available radio images indicate that at most 20% of sources might be genuine X-shaped radio sources that could have formed by a restarting of beams in a new direction following an interruption and axis flip. The remaining 20% are in neither of these categories.These images indicate that at most a small fraction of the candidates might be genuine X-shaped radio sources that were formed by a restarting of beams in a new direction following a major merger, or by spin drift caused by BH-BH interaction. This suggests that fewer than 1.3% of extended radio sources appear to be candidates for genuine axis reorientations (“spin flips”), or 2.2% if possible “axis drift” sources are included, much smaller than the 7% suggested by Leahy & Parma. Thus, the associated GWB may be substantially smaller than previous estimates. These results can be used to normalize detailed calculations of the SMBH coalescence rate and the GWB.

ULTRAMASSIVE BLACK HOLE COALESCENCE

Although supermassive black holes (SMBHs) correlate well with their host galaxies, there is an emerging view that outliers exist. Henize 2-10, NGC 4889, and NGC1277 are examples of SMBHs at least an order of magnitude more massive than their host galaxy suggests. The dynamical effects of such ultramassive central black holes is unclear. Here, we perform direct N-body simulations of mergers of galactic nuclei where one black hole is ultramassive to study the evolution of the remnant and the black hole dynamics in this extreme regime. We find that the merger remnant is axisymmetric near the center, while near the large SMBH influence radius, the galaxy is triaxial. The SMBH separation shrinks rapidly due to dynamical friction, and quickly forms a binary black hole; if we scale our model to the most massive estimate for the NGC1277 black hole, for example, the timescale for the SMBH separation to shrink from nearly a kiloparsec to less than a parsec is roughly 10 Myr. By the time the SMBHs form a hard binary, gravitational wave emission dominates, and the black holes coalesce in a mere few Myr. Curiously, these extremely massive binaries appear to nearly bypass the 3-body scattering evolutionary phase. Our study suggests that in this extreme case, SMBH coalescence is governed by dynamical friction followed nearly directly by gravitational wave emission, resulting in an rapid and efficient SMBH coalescence timescale. We discuss the implications for gravitational wave event rates and hypervelocity star production.

RECOILING SUPERMASSIVE BLACK HOLES: A SEARCH IN THE NEARBY UNIVERSE

The coalescence of a binary black hole can be accompanied by a large gravitational recoil due to anisotropic emission of gravitational waves. A recoiling supermassive black hole (SBH) can subsequently undergo long-lived oscillations in the potential well of its host galaxy, suggesting that offset SBHs may be common in the cores of massive ellipticals. We have analyzed HST archival images of 14 nearby core ellipticals, finding evidence for small ($\lesssim 10$ pc) displacements between the AGN (locating the SBH) and the center of the galaxy (the mean photocenter) in 10 of them. Excluding objects that may be affected by large-scale isophotal asymmetries, we consider six galaxies to have detected displacements, including M87, where a displacement was previously reported by Batcheldor et al. 2010. In individual objects, these displacements can be attributed to residual gravitational recoil oscillations following a major or minor merger within the last few Gyr. For plausible merger rates, however, there is a high probability of larger displacements than those observed, if SBH coalescence took place in these galaxies. Remarkably, the AGN-photocenter displacements are approximately aligned with the radio source axis in four of the six galaxies with displacements, including three of the four having relatively powerful kpc-scale jets. This suggests intrinsic asymmetries in radio jet power as a possible displacement mechanism, although approximate alignments are also expected for gravitational recoil. Orbital motion in SBH binaries and interactions with massive perturbers can produce the observed displacement amplitudes but do not offer a ready explanation for the alignments.

Prompt tidal disruption of stars as an electromagnetic signature of supermassive black hole coalescence

A precise electromagnetic measurement of the sky coordinates and redshift of a coalescing black hole binary holds the key for using its gravitational wave (GW) signal to constrain cosmological parameters and to test general relativity. Here we show that the merger of ~10^{6-7}M_sun black holes is generically followed over a period of years by multiple electromagnetic flares from tidally disrupted stars. The sudden recoil imparted to the merged black hole by GW emission promptly fills its loss cone and results in a tidal disruption rate of stars as high as ~0.1 per year. The prompt disruption of a star within a single galaxy over a short period provides a unique electromagnetic flag of a recent black hole coalescence event, and sequential disruptions could be used on their own to calibrate the expected rate of GW sources for pulsar timing arrays or the proposed Laser Interferometer Space Antenna (LISA).

DETECTION, LOCALIZATION, AND CHARACTERIZATION OF GRAVITATIONAL WAVE BURSTS IN A PULSAR TIMING ARRAY

Efforts to detect gravitational waves by timing an array of pulsars have focused traditionally on stationary gravitational waves: e.g., stochastic or periodic signals. Gravitational wave bursts --- signals whose duration is much shorter than the observation period --- will also arise in the pulsar timing array waveband. Sources that give rise to detectable bursts include the formation or coalescence of supermassive black holes (SMBHs), the periapsis passage of compact objects in highly elliptic or unbound orbits about a SMBH, or cusps on cosmic strings. Here we describe how pulsar timing array data may be analyzed to detect and characterize these bursts. Our analysis addresses, in a mutually consistent manner, a hierarchy of three questions: \emph{i}) What are the odds that a dataset includes the signal from a gravitational wave burst? \emph{ii}) Assuming the presence of a burst, what is the direction to its source? and \emph{iii}) Assuming the burst propagation direction, what is the burst waveform's time dependence in each of its polarization states? Applying our analysis to synthetic data sets we find that we can \emph{detect} gravitational waves even when the radiation is too weak to either localize the source of infer the waveform, and \emph{detect} and \emph{localize} sources even when the radiation amplitude is too weak to permit the waveform to be determined. While the context of our discussion is gravitational wave detection via pulsar timing arrays, the analysis itself is directly applicable to gravitational wave detection using either ground or space-based detector data.

Brightening of an Accretion Disk due to Viscous Dissipation of Gravitational Waves during the Coalescence of Supermassive Black Holes

Mergers of supermassive black hole binaries release peak power of up to approximately 10(57) erg s(-1) in gravitational waves (GWs). As the GWs propagate through ambient gas, they induce shear and a small fraction of their power is dissipated through viscosity. The dissipated heat appears as electromagnetic (EM) radiation, providing a prompt EM counterpart to the GW signal. For thin accretion disks, the GW heating rate exceeds the accretion power at distances farther than approximately 10(3) Schwarzschild radii, independently of the accretion rate and viscosity coefficient.

The evolution of black hole mass and spin in active galactic nuclei

Observations show that the central black hole in galaxies has a mass M of only ∼10−3 of the stellar bulge mass. Thus, whatever process grows the black hole also promotes star formation with far higher efficiency. We interpret this in terms of the generic tendency of active galactic nucleus (AGN) accretion discs to become self-gravitating outside some small radius Rsg∼ 0.01–0.1 pc from the black hole. We argue that mergers consist of sequences of such episodes, each limited by self-gravity to a mass ΔMepisode∼ 10−3M, with angular momentum characteristic of the small part of the accretion flow which formed it. In this picture, a major merger with ΔMmerger∼M gives rise to a long series of low-mass accretion disc episodes, all chaotically oriented with respect to one another. Thus, the angular momentum vector oscillates randomly during the accretion process, on mass-scales ∼ 103 times smaller than the total mass accreted in a major merger event. We show that for essentially all AGN parameters, the disc produced by any accretion episode of this type has lower angular momentum than the hole, allowing stable co- and counter-alignment of the discs through the Lense–Thirring effect. A sequence of randomly oriented accretion episodes as envisaged above then produces accretion discs stably co- or counter-aligned with the black hole spin with almost equal frequency. Accretion from these discs very rapidly adjusts the hole's spin parameter to average values (the precise range depending slightly on the disc vertical viscosity coefficient α2) from any initial conditions, but with significant fluctuations (Δa∼±0.2) about these. We conclude that (i) supermassive black holes (SMBH) should on average spin moderately, with the mean value decreasing slowly as the mass increases; (ii) SMBH coalescences leave little long-term effect on ⁠; (iii) SMBH coalescence products in general have modest recoil velocities, so that there is little likelihood of their being ejected from the host galaxy; (iv) black holes can grow even from stellar masses to ∼5 × 109M⊙ at high redshift z∼ 6; and (v) jets produced in successive accretion episodes can have similar directions, but after several episodes the jet direction deviates significantly. Rare examples of massive holes with larger spin parameters could result from prograde coalescences with SMBHs of similar mass, and are most likely to be found in giant ellipticals. We compare these results with observation.

Publisher’s Note: Higher signal harmonics, LISA’s angular resolution, and dark energy [Phys. Rev. D76, 104016 (2007)

It is generally believed that the angular resolution of the Laser Interferometer Space Antenna (LISA) for binary supermassive black holes (SMBH) will not be good enough to identify the host galaxy or galaxy cluster. This conclusion, based on using only the dominant harmonic of the binary SMBH signal, changes substantially when higher signal harmonics are included in assessing the parameter estimation problem. We show that in a subset of the source parameter space the angular resolution increases by more than a factor of 10, thereby making it possible for LISA to identify the host galaxy/galaxy cluster. Thus, LISA's observation of certain binary SMBH coalescence events could constrain the dark energy equation of state to within a few percent, comparable to the level expected from other dark energy missions.

Higher signal harmonics, LISA’s angular resolution, and dark energy

It is generally believed that the angular resolution of the Laser Interferometer Space Antenna (LISA) for binary supermassive black holes (SMBH) will not be good enough to identify the host galaxy or galaxy cluster. This conclusion, based on using only the dominant harmonic of the binary SMBH signal, changes substantially when higher signal harmonics are included in assessing the parameter estimation problem. We show that in a subset of the source parameter space the angular resolution increases by more than a factor of 10, thereby making it possible for LISA to identify the host galaxy/galaxy cluster. Thus, LISA’s observation of certain binary SMBH coalescence events could constrain the dark energy equation of state to within a few percent, comparable to the level expected from other dark energy missions.

Finding the Electromagnetic Counterparts of Cosmological Standard Sirens

The gravitational waves (GWs) emitted during the coalescence of supermassive black holes (SMBHs) in the mass range ~104-107 M?/(1 + z) will be detectable out to high redshifts with the future Laser Interferometer Space Antenna (LISA). The distance and direction to these standard sirens can be inferred directly from the GW signal, with a precision that depends on the masses, spins, and geometry of the merging system. In a given cosmology, the LISA-measured luminosity distance translates into a redshift shell. We calculate the size and shape of the corresponding three-dimensional error volume in which an electromagnetic counterpart to a LISA event could be found, taking into account errors in the background cosmology (as expected by the time LISA flies), weak gravitational lensing (de)magnification due to inhomogeneities along the line of sight, and potential source-peculiar velocities. Weak-lensing errors largely exceed other sources of uncertainties (by a factor of ~7 for typical sources at z = 1). Under the plausible assumption that SMBH-SMBH mergers are accompanied by gas accretion leading to Eddington-limited quasar activity, we then compute the number of quasars that would be found in a typical three-dimensional LISA error volume, as a function of BH mass and event redshift. Low redshifts offer the best opportunities to identify quasar counterparts to cosmological standard sirens. For mergers of ~4 ? (105-107) M? SMBHs, the LISA error volume will typically contain a single near-Eddington quasar at z ~ 1. If SMBHs are spinning rapidly, the error volume is smaller and may contain a unique quasar out to redshift z ~ 3. This will allow a straightforward test of the hypothesis that GW events are accompanied by bright quasar activity and, if the hypothesis proves correct, will guarantee the identification of a unique quasar counterpart to a LISA event, with a B-band luminosity of LB ~ (1010-1011) L?. Robust counterpart identifications would allow unprecedented tests of the physics of SMBH accretion, such as precision measurements of the Eddington ratio. They would clarify the role of gas as a catalyst in SMBH coalescences and would also offer an alternative method to constrain cosmological parameters.

Realistic event rates for detection of supermassive black hole coalescence by LISA

The gravitational waves generated during supermassive black hole (SMBH) coalescence are prime candidates for detection by the satellite LISA. We use the extended Press—Schechter formalism combined with empirically motivated estimates for the SMBH—dark matter halo mass relation and SMBH occupation fraction to estimate the maximum coalescence rate for major SMBH mergers. Assuming efficient binary coalescence, and guided by the lowest nuclear black hole mass inferred in local galactic bulges and nearby low-luminosity active galactic nuclei (≈105 M⊙), we predict approximately 15 detections per year at a signal-to-noise ratio greater than 5, in each of the inspiral and ringdown phases. Rare coalescences between SMBHs having masses in excess of 107 M⊙ will be more readily detected via gravitational waves from the ringdown phase.