Untangling the merger history of massive black holes with LISA

Binary black hole coalescences emit gravitational waves that will be measurable by the space-based detector LISA to large redshifts. This suggests that LISA may be able to observe black holes grow and evolve as the Universe evolves, mapping the distribution of black hole masses as a function of redshift. An immediate difficulty with this idea is that LISA measures certain redshifted combinations of masses with good accuracy: if a system has some mass parameter m , then LISA measures  (1+ z ) m   . This mass–redshift degeneracy makes it difficult to follow the mass evolution. In many cases, LISA will also measure the luminosity distance D of a coalescence accurately. Since cosmological parameters (particularly the mean density, the cosmological constant and the Hubble constant) are now known with moderate precision, we can obtain z from D and break the degeneracy. This makes it possible to untangle the mass and redshift and to study the mass and merger history of black holes. Mapping the black hole mass distribution could open a window on to an early epoch of structure formation.     

Laser Interferometer Space Antenna double black holes: dynamics in gaseous nuclear discs

We study the inspiral of double black holes, with masses in the Laser Interferometer Space Antenna ( LISA ) window of detectability, orbiting inside a massive circumnuclear, rotationally supported gaseous disc. Using high-resolution smoothed particle hydrodynamics simulations, we follow the black hole dynamics in the early phase when gas-dynamical friction acts on the black holes individually, and continue our simulation until they form a close binary. We find that in the early sinking the black holes lose memory of their initial orbital eccentricity if they corotate with the gaseous disc. As a consequence, the massive black holes bind forming a binary with a low eccentricity, consistent with zero within our numerical resolution limit. The cause of circularization resides in the rotation present in the gaseous background where dynamical friction operates. Circularization may hinder gravitational waves from taking over and leading the binary to coalescence. In the case of counter-rotating orbits, the initial eccentricity (if present) does not decrease, and the black holes may bind forming an eccentric binary. When dynamical friction has subsided, for equal mass black holes and regardless their initial eccentricity, angular momentum loss, driven by the gravitational torque exerted on the binary by surrounding gas, is nevertheless observable down to the smallest scale probed (≃1 pc). In the case of unequal masses, dynamical friction remains efficient down to our resolution limit, and there is no sign of formation of any ellipsoidal gas distribution that may further harden the binary. During inspiral, gravitational capture of gas by the black holes occurs mainly along circular orbits; eccentric orbits imply high relative velocities and weak gravitational focusing. Thus, the active galactic nucleus activity may be excited during the black hole pairing process and double active nuclei may form when circularization is completed, on distance scales of tens of parsecs.     

The evolution of massive black hole seeds

We investigate the evolution of high-redshift seed black hole masses at late times and their observational signatures. The massive black hole seeds studied here form at extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating dark matter haloes with seeds formed in this way, we follow the mass assembly of these black holes to the present time using a Monte Carlo merger tree. Using this machinery, we predict the black hole mass function at high redshifts and at the present time, the integrated mass density of black holes and the luminosity function of accreting black holes as a function of redshift. These predictions are made for a set of three seed models with varying black hole formation efficiency. Given the accuracy of present observational constraints, all three models can be adequately fitted. Discrimination between the models appears predominantly at the low-mass end of the present-day black hole mass function which is not observationally well constrained. However, all our models predict that low surface brightness, bulgeless galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts has a direct influence on the black hole occupation fraction in galaxies at   z = 0  . This effect is more pronounced for low-mass galaxies. This is the key discriminant between the models studied here and the Population III remnant seed model. We find that there exist a population of low-mass galaxies that do not host nuclear black holes. Our prediction of the shape of the M _BH–σ relation at the low-mass end is in agreement with the recent observational determination from the census of low-mass galaxies in the Virgo cluster.     

The imprint of massive black hole formation models on the LISA data stream

The formation, merging and accretion history of massive black holes (MBHs) along the hierarchical build-up of cosmic structures leaves a unique imprint on the background of gravitational waves (GWs) at mHz frequencies. We study here, by means of dedicated simulations of black hole build-up, the possibility of constraining different models of black hole cosmic evolution using future GW space-borne missions, such as LISA . We consider two main scenarios for black hole formation, namely, one where seeds are light (  ≃10² M_⊙  , remnant of Population III stars) and one where seeds are heavy (  ≳10⁴ M_⊙  , direct collapse). In all the models we have investigated, MBH binary coalescences do not produce a stochastic GW background, but rather, a set of individual resolved events. Detection of several hundreds merging events in a 3-yr LISA mission will be the sign of a heavy seed scenario with efficient formation of black hole seeds in a large fraction of high-redshift haloes. At the other extreme, a low event rate, about a few tens in 3 yr, is peculiar of scenarios where either the seeds are light, and many coalescences do not fall into the LISA band, or seeds are massive, but rare. In this case a decisive diagnostic is provided by the shape of the mass distribution of detected events. Light binaries  ( m < 10⁴ M_⊙)  are predicted in a fairly large number in Population III remnant models, but are totally absent in direct collapse models. Finally, a further, helpful diagnostic of black hole formation models lies in the distribution of the mass ratios in binary coalescences. While heavy seed models predict that most of the detected events involve equal-mass binaries, in the case of light seeds, mass ratios are equally distributed in the range 0.1–1.     

Implications of primordial black holes on the first stars and the origin of the super-massive black holes

If the cosmological dark matter has a component made of small primordial black holes (BHs), they may have a significant impact on the physics of the first stars and on the subsequent formation of massive BHs. Primordial BHs would be adiabatically contracted into these stars and then would sink to the stellar centre by dynamical friction, creating a larger BH which may quickly swallow the whole star. If these primordial BHs are heavier than  ∼10²² g  , the first stars would likely live only for a very short time and would not contribute much to the reionization of the Universe. They would instead become  10–10³ M_⊙  BHs which (depending on subsequent accretion) could serve as seeds for the super-massive BHs seen at high redshifts as well as those inside galaxies today.     

The gravitational wave background from star–massive black hole fly-bys

Stars on eccentric orbits around a massive black hole (MBH) emit bursts of gravitational waves (GWs) at periapse. Such events may be directly resolvable in the Galactic Centre. However, if the star does not spiral in, the emitted GWs are not resolvable for extragalactic MBHs, but constitute a source of background noise. We estimate the power spectrum of this extreme mass ratio burst background (EMBB) and compare it to the anticipated instrumental noise of the Laser Interferometer Space Antenna (LISA). To this end, we model the regions close to an MBH, accounting for mass segregation, and for processes that limit the presence of stars close to the MBH, such as GW inspiral and hydrodynamical collisions between stars. We find that the EMBB is dominated by GW bursts from stellar mass black holes, and the magnitude of the noise spectrum  ( fS _GW)^1/2  is at least a factor of ∼10 smaller than the instrumental noise. As an additional result of our analysis, we show that LISA is unlikely to detect relativistic bursts in the Galactic Centre.     

The black hole mass – spheroid luminosity relation

The differing   M _bh– L   relations presented in McLure & Dunlop, Marconi & Hunt and Erwin et al. have been investigated. A number of issues have been identified and addressed in each of these studies, including but not limited to the removal of a dependency on the Hubble constant, a correction for dust attenuation in the bulges of disc galaxies, the identification of lenticular galaxies previously treated as elliptical galaxies and the application of the same ( Y ∣ X ) regression analysis. These adjustments result in relations which now predict similar black hole masses. The optimal K -band relation is  log( M _bh/M_⊙) =−0.37(±0.04)( M _K + 24) + 8.29(±0.08)  , with a total (not intrinsic) scatter in log M _bh equal to 0.33 dex. This level of scatter is similar to the value of 0.34 dex from the     relation of Tremaine et al. and compares favourably with the value of 0.31 dex from the   M _bh– n   relation of Graham & Driver. Using different photometric data, consistent relations in the B and R band are also provided, although we do note that the small  ( N = 13)  R -band sample used by Erwin et al. is found here to have a slope of −0.30 ± 0.06. Performing a symmetrical regression on the larger K -band sample gives a slope of ∼−0.40, implying M _bh∝ L ^1.00. Implications for galaxy–black hole co-evolution, in terms of dry mergers, are briefly discussed, as are the predictions for intermediate mass black holes. Finally, as noted by others, a potential bias in the galaxy sample used to define the   M _bh– L   relations is shown and a corrective formula provided.     

Dual black holes in merger remnants – I. Linking accretion to dynamics

We study the orbital evolution and accretion history of massive black hole (MBH) pairs in rotationally supported circumnuclear discs up to the point where MBHs form binary systems. Our simulations have high resolution in mass and space which, for the first time, makes it feasible to follow the orbital decay of a MBH either counter- or corotating with respect to the circumnuclear disc. We show that a moving MBH on an initially counter-rotating orbit experiences an 'orbital angular momentum flip' due to the gas-dynamical friction, i.e. it starts to corotate with the disc before a MBH binary forms. We stress that this effect can only be captured in very high resolution simulations. Given the extremely large number of gas particles used, the dynamical range is sufficiently large to resolve the Bondi–Hoyle–Lyttleton radii of individual MBHs. As a consequence, we are able to link the accretion processes to the orbital evolution of the MBH pairs. We predict that the accretion rate is significantly suppressed and extremely variable when the MBH is moving on a retrograde orbit. It is only after the orbital angular momentum flip has taken place that the secondary rapidly 'lights up' at which point both MBHs can accrete near the Eddington rate for a few Myr. The separation of the double nucleus is expected to be around ≲10 pc at this stage. We show that the accretion rate can be highly variable also when the MBH is corotating with the disc (albeit to a lesser extent) provided that its orbit is eccentric. Our results have significant consequences for the expected number of observable double active galactic nuclei at separations of ≲100 pc.     

The tidal disruption rate in dense galactic cusps containing a supermassive binary black hole

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes, the gravitational potential is dominated by the more massive primary black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power-law cusp. We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted by the larger black hole. This effect is analogous to the so-called Kozai effect well known in celestial mechanics. We develop an analytical theory for the secular evolution of the stellar orbits and calculate the rate of tidal disruption. We compare our analytical theory with a simple numerical model and find very good agreement.
Our results show that for a primary black hole mass of  ∼10⁶–10⁷ M_⊙  , the black hole mass-ratio   q > 10⁻²  , cusp size ∼1 pc, the tidal disruption rate can be as large as  ∼10⁻²–1 M_⊙ yr⁻¹  . This is at least 10²–10⁴ times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical-friction time-scale, and it is rather short: ∼10⁵ yr. The dependence of the tidal disruption rate on the mass ratio, and on the size of the cusp, is also discussed.     

Gravitational waves from scattering of stellar-mass black holes in galactic nuclei

Stellar-mass black holes (BHs) are expected to segregate and form a steep density cusp around supermassive black holes (SMBHs) in galactic nuclei. We follow the evolution of a multimass system of BHs and stars by numerically integrating the Fokker–Planck energy diffusion equations for a variety of BH mass distributions. We find that the BHs 'self-segregate', and that the rarest, most massive BHs dominate the scattering rate closest to the SMBH  (≲10⁻¹ pc)  . BH–BH binaries form out of gravitational wave emission during BH encounters. We find that the expected rate of BH coalescence events detectable by Advanced LIGO is  ∼1–10² yr⁻¹  , depending on the initial mass function of stars in galactic nuclei and the mass of the most massive BHs. We find that the actual merger rate is likely ∼10 times larger than this due to the intrinsic scatter of stellar densities in many different galaxies. The BH binaries that form this way in galactic nuclei have significant eccentricities as they enter the LIGO band (90 per cent with   e > 0.9  ), and are therefore distinguishable from other binaries, which circularize before becoming detectable. We also show that eccentric mergers can be detected to larger distances and greater BH masses than circular mergers, up to  ∼700 M_⊙  . Future ground-based gravitational wave observatories will be able to constrain both the mass function of BHs and stars in galactic nuclei.     

The black hole mass–galaxy age relation

We present evidence that there is a significant correlation between the fraction of the mass of a galaxy that lies in its central black hole and the age of the galactic stellar population. Since the absorption-line indices that are used to estimate the age are luminosity-weighted, they essentially measure the time since the last significant episode of star formation in the galaxy. The existence of this correlation is consistent with several theories of galaxy formation, including the currently favoured hierarchical picture of galaxy evolution, which predicts just such a relation between the black hole mass and the time since the last burst of merger-induced star formation. It is not consistent with models in which the massive black hole is primordial, and hence uncoupled from the stellar properties of the galaxy.     

Star clusters around recoiled black holes in the Milky Way halo

Gravitational wave emission by coalescing black holes (BHs) kicks the remnant BH with a typical velocity of hundreds of  km s⁻¹  . This velocity is sufficiently large to remove the remnant BH from a low-mass galaxy but is below the escape velocity from the Milky Way (MW) galaxy. If central BHs were common in the galactic building blocks that merged to make the MW, then numerous BHs that were kicked out of low-mass galaxies should be freely floating in the MW halo today. We use a large statistical sample of possible merger tree histories for the MW to estimate the expected number of recoiled BH remnants present in the MW halo today. We find that hundreds of BHs should remain bound to the MW halo after leaving their parent low-mass galaxies. Each BH carries a compact cluster of old stars that populated the core of its original host galaxy. Using the time-dependent Fokker–Planck equation, we find that the present-day clusters are  ≲1 pc  in size, and their central bright regions should be unresolved in most existing sky surveys. These compact systems are distinguishable from globular clusters by their internal (Keplerian) velocity dispersion greater than 100 km s⁻¹ and their high mass-to-light ratio owing to the central BH. An observational discovery of this relic population of star clusters in the MW halo would constrain the formation history of the MW and the dynamics of BH mergers in the early Universe. A similar population should exist around other galaxies and may potentially be detectable in M31 and M33.