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.     

Gravitational waves from resolvable massive black hole binary systems and observations with Pulsar Timing Arrays

Massive black holes are key components of the assembly and evolution of cosmic structures, and a number of surveys are currently on going or planned to probe the demographics of these objects and to gain insight into the relevant physical processes. Pulsar Timing Arrays (PTAs) currently provide the only means to observe gravitational radiation from massive black hole binary systems with masses  ≳10⁷ M_⊙  . The whole cosmic population produces a stochastic background that could be detectable with upcoming PTAs. Sources sufficiently close and/or massive generate gravitational radiation that significantly exceeds the level of the background and could be individually resolved. We consider a wide range of massive black hole binary assembly scenarios, investigate the distribution of the main physical parameters of the sources, such as masses and redshift, and explore the consequences for PTAs observations. Depending on the specific massive black hole population model, we estimate that on average at least one resolvable source produces timing residuals in the range  ∼5–50 ns  . PTAs, and in particular the future Square Kilometre Array, can plausibly detect these unique systems, although the events are likely to be rare. These observations would naturally complement on the high-mass end of the massive black hole distribution function future surveys carried out by the Laser Interferometer Space Antenna .     

Gravitational wave bursts from collisions of primordial black holes in clusters

The rate of gravitational wave bursts from the mergers of massive primordial black holes in clusters is calculated. Such clusters of black holes can be formed through phase transitions in the early Universe. The central black holes in clusters can serve as the embryos of supermassive black holes in galactic nuclei. The expected burst detection rate by the LISA gravitational wave detector is estimated.     

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.     

Non-thermal transient sources from rotating black holes

Rotating black holes can power the most extreme non-thermal transient sources. They have a long-duration viscous time-scale of spin-down, and produce non-thermal emissions along their spin-axis, powered by a relativistic capillary effect. We report on the discovery of exponential decay in Burst and Triensient Source Experiment (BATSE) light curves of long gamma-ray bursts (GRBs) by matched filtering, consistent with a viscous time-scale, and identify ultra-high energy cosmic rays (UHECRs) about the Greisen–Zatsepin–Kuzmin (GZK) threshold with linear acceleration of ion contaminants along the black hole spin-axis, consistent with black hole masses and lifetimes of Fanaroff–Riley type II (FR II) active galactic nuclei (AGN). We explain the absence of UHECRs from BL Lac objects due to UHECR emissions preferably at appreciable angles away from the black hole spin-axis. Black hole spin may be the key to unification of GRBs and their host environments, and to AGN and their host galaxies. Our model points to long-duration bursts in radio from long GRBs without supernovae and gravitational waves from all long GRBs.     

Gravitational wave emission from a companion black hole in the presence of an accretion disc around a super-massive Kerr black hole

Gravitational wave signal characteristics from a binary black hole system in which the companion moves through the accretion disc of the primary are studied. We chose the primary to be a super-massive  ( M = 10⁸ M_⊙)  Kerr black hole and the companion to be a massive black hole  ( M = 10⁵ M_⊙)  to clearly demonstrate the effects. We show that the drag exerted on the companion by the disc is sufficient to reduce the coalescence time of the binary. The drag is primarily due to the fact that the accretion disc on a black hole deviates from a Keplerian disc and becomes sub-Keplerian due to inner boundary condition on the black hole horizon. We consider two types of accretion rates on to the companion. The companion is deeply immersed inside the disc and it can accrete at the Bondi rate which depends on the instantaneous density of the disc. However, an accretion disc can also form around the smaller black hole and it can accrete at its Eddington rate. Thus, this case is also studied and the results are compared. We find that the effect of the disc will be significant in reducing the coalescence time and one needs to incorporate this while interpreting gravitational wave signals emitted from such a binary system.     

A simple model for lensing by black holes in galactic nuclei

The lensing properties of the Plummer model with a central point mass and external shear are derived, including the image multiplicities, critical curves and caustics. This provides a simple model for a flattened galaxy with a central supermassive black hole. For the Plummer model with black hole, the maximum number of images is four, provided the black hole mass is less than an upper bound which is calculated analytically. This introduces a method to constrain black hole masses by counting images, thus applicable at cosmological distance. With shear, the maximum number of images is six and we illustrate the occurrence of an astroid caustic and two metamorphoses.     

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.     

Studies of accretion flows around rotating black holes – III. Shock oscillations and an estimation of the spin parameter from QPO frequencies

In the present communication of our series of papers dealing with the accretion flows in the pseudo-Kerr geometry, we discuss the effects of viscosity on the accretion flow around a rotating black hole. We find the solution topologies and give special attention to the solutions containing shocks. We draw the parameter space where standing shocks are possible and where the shocks could be oscillating and could produce quasi-periodic oscillations (QPOs) of X-rays observed from black hole candidates. In this model, the extreme locations of the shocks give the upper limits of the QPO frequencies  (ν_QPO)  which could be observed. We show that both the viscosity of the flow and the spin of the black hole a increase the QPO frequency while, as expected, the black hole mass reduces the QPO frequencies. Our major conclusion is that the highest observed frequency gives a strict lower limit of the spin. For instance, a black hole exhibiting  ν_QPO∼ 400  and  700 Hz  must have the spin parameters of   a > 0.25  and  >0.75  , respectively, provided viscosity of the flow is small. We discuss the implications of our results in the light of observations of QPOs from black hole candidates.     

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.     

Spectral variability in transonic discs around black holes

Transonic discs with accretion rates relevant to intrinsically bright Galactic X-ray sources ( L ≈10³⁸–10³⁹ erg s⁻¹) exhibit a time-dependent cyclic behaviour due to the onset of a thermal instability driven by radiation pressure. In this paper we calculate radiation spectra emitted from thermally unstable discs to provide detailed theoretical predictions for observationally relevant quantities. The emergent spectrum has been obtained by solving self-consistently the vertical structure and radiative transfer in the disc atmosphere. We focus on four particular stages of the disc evolution, the maximal evacuation stage and three intermediate stages during the replenishment phase. The disc is found to undergo rather dramatic spectral changes during the evolution, emitting mainly in the 1–10 keV band during outburst and in the 0.1–1 keV band off-outburst. Local spectra, although different in shape from a blackbody at the disc effective temperature, may be characterized in terms of a hardening factor f . We have found that f is more or less constant, both in radius and in time, with a typical value ∼ 1.65.     

The runaway instability of self-gravitating tori with non-constant specific angular momentum around black holes

We discuss the runaway instability of axisymmetric tori with non-constant specific angular momentum around black holes, taking into account self-gravity of the tori. The distribution of specific angular momentum of the tori is assumed to be a positive power law with respect to the distance from the rotational axis. By employing the pseudo-Newtonian potential for the gravity of the spherical black hole, we have found that self-gravity of the tori causes a runaway instability if the amount of the mass which is transferred from the torus to the black hole exceeds a critical value, i.e. 3 per cent of the mass of the torus. This has been shown by two different approaches: (1) by using equilibrium models and (2) by dynamical simulations. In particular, dynamical simulations using an SPH code have been carried out for both self-gravitating and non-self-gravitating tori. For non-self-gravitating models, all tori are runaway stable. Therefore we come to the conclusion that self-gravity of the tori has a stronger destabilizing effect than the stabilizing effect of the positive power-law distribution of the angular momentum.     

Behaviour of dissipative accretion flows around black holes

We investigate the behaviour of dissipative accreting matter close to a black hole, as this provides important observational features of galactic and extragalactic black hole candidates. We find a complete set of global solutions in the presence of viscosity and synchrotron cooling. We show that advective accretion flow can have a standing shock wave and the dynamics of the shock is controlled by the dissipation parameters (both viscosity and cooling). We study the effective region of the parameter space for standing as well as oscillating shock. We find that the shock front always moves towards the black hole as the dissipation parameters are increased. However, viscosity and cooling have opposite effects in deciding the solution topologies. We obtain two critical cooling parameters that separate the nature of the accretion solution.     

On an excitation mechanism for trapped inertial waves in discs around black holes

According to one model, high-frequency quasi-periodic oscillations (QPOs) can be identified with inertial waves, trapped in the inner regions of accretion discs around black holes due to relativistic effects. In order to be detected, their amplitudes need to reach large enough values via some excitation mechanism. We work out in detail a non-linear coupling mechanism suggested by Kato, in which a global warping or eccentricity of the disc has a fundamental role. These large-scale deformations combine with trapped modes to generate 'intermediate' waves of negative energy that are damped as they approach either their corotation resonance or the inner edge of the disc, resulting in amplification of the trapped waves. We determine the growth rates of the inertial modes, as well as their dependence on the spin of the black hole and the properties of the disc. Our results indicate that this coupling mechanism can provide an efficient excitation of trapped inertial waves, provided the global deformations reach the inner part of the disc with non-negligible amplitude.     

Gravitational instability of discs with dissipative coronae around supermassive black holes

We study the dynamical structure of a self-gravitating disc with coronae around a supermassive black hole. Assuming that the magnetorotational instability responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has a major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscosity coefficient, the mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic Centre, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic Centre, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. The existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.