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.     

A systematic study of variability in a sample of ultraluminous X-ray sources

We present results from a study of short-term variability in 19 archival observations by XMM–Newton of 16 ultraluminous X-ray sources (ULXs). Eight observations (six sources) showed intrinsic variability with power spectra in the form of either a power-law or broken power-law-like continuum and in some cases quasi-periodic oscillations (QPOs). The remaining observations were used to place upper limits on the strength of possible variability hidden within. Seven observations (seven sources) yielded upper limits comparable to, or higher than, the values measured from those observations with detectable variations. These represented the seven faintest sources, all with   f_x < 3 × 10⁻¹² erg cm⁻² s⁻¹  . In contrast, there are four observations (three sources) that gave upper limits significantly lower than both the values measured from the ULX observations with detectable variations, and the values expected by comparison with luminous Galactic black hole X-ray binaries (BHBs) and active galactic nuclei (AGN) in the observed frequency bandpass (10⁻³–1 Hz). This is the case irrespective of whether one assumes characteristic frequencies appropriate for a stellar mass  (10 M_⊙)  or an intermediate mass  (1000 M_⊙)  black hole, and means that in some ULXs the variability is significantly suppressed compared to bright BHBs and AGN. We discuss ways to account for this unusual suppression in terms of both observational and intrinsic effects and whether these solutions are supported by our results.     

Spectral variability of ultraluminous X-ray sources

We study spectral variability of 11 ultraluminous X-ray sources (ULX) using archived XMM–Newton and Chandra observations. We use three models to describe the observed spectra: a power law, a multicolour disc (MCD) and a combination of these two models. We find that seven ULXs show a correlation between the luminosity L _X and the photon index Γ. Furthermore, four out of these seven ULXs also show spectral pivoting in the observed energy band. We also find that two ULXs show an   L _X–Γ  anticorrelation. The spectra of four ULXs in the sample can be adequately fitted with a MCD model. We compare these sources to known black hole binaries (BHB) and find that they follow similar paths in their luminosity–temperature diagrams. Finally, we show that the 'soft excess' reported for many of these ULXs at ∼0.2 keV seems to roughly follow a trend   L _soft∝ T ^−3.5  when modelled with a power law plus a 'cool' MCD model. This is contrary to the   L ∝ T ⁴  relation that is expected from theory and what is seen for many accreting BHBs. The observed trend could instead arise from disc emission beamed by an outflowing wind around a  ∼10 M_⊙  black hole.     

The changing interstellar medium of massive elliptical galaxies and cosmic evolution of radio galaxies and quasars

The recently discovered apparent dramatic expansion in the effective radii of massive elliptical galaxies from   z ≃ 2  to ≃0.1 has been interpreted in terms of either galaxy mergers or the rapid loss of cold gas due to active galactic nuclei (AGN) feedback. In examining the latter case, we have quantified the extent of the expansion, which is uncertain observationally, in terms of the star formation parameters and time of the expulsion of the cold gas. In either case, the large global decrease in stellar density should translate into a major drop in the interstellar medium density and pressure with cosmic epoch. These cosmological changes are expected to have a major influence on the gas accretion mode, which will shift from 'cold' thin disc accretion at high redshifts towards 'hot' Bondi fed Advection Dominated Accretion Flow (ADAF) accretion at low redshifts. The decline of angular momentum inflow would then lead to a spin down of the black hole, for which we have calculated more precise time-scales; a value of about 0.2 Gyr is typical for a  10⁹ M_⊙  central black hole. These results have implications for the different cosmological evolutionary patterns found for the luminosity functions of powerful and weak radio galaxies.     

Quasi-viscous accretion flow – I. Equilibrium conditions and asymptotic behaviour

In a novel approach to studying viscous accretion flows, viscosity has been introduced as a perturbative effect, involving a first-order correction in the α-viscosity parameter. This method reduces the problem of solving a second-order non-linear differential equation (Navier–Stokes equation) to that of an effective first-order equation. Viscosity breaks down the invariance of the equilibrium conditions for stationary inflow and outflow solutions, and distinguishes accretion from wind. Under a dynamical systems classification, the only feasible critical points of this 'quasi-viscous' flow are saddle points and spirals. On large spatial scales of the disc, where a linearized and radially propagating time-dependent perturbation is known to cause a secular instability, the velocity evolution equation of the quasi-viscous flow has been transformed to bear a formal closeness with Schrödinger's equation with a repulsive potential. Compatible with the transport of angular momentum to the outer regions of the disc, a viscosity-limited length-scale has been defined for the full spatial extent over which the accretion process would be viable.     

Self-similar solution of hot accretion flows with ordered magnetic field and outflow

Observations and numerical magnetohydrodynamic (MHD) simulations indicate the existence of outflows and ordered large-scale magnetic fields in the inner region of hot accretion flows. In this paper, we present the self-similar solutions for advection-dominated accretion flows (ADAFs) with outflows and ordered magnetic fields. Stimulated by numerical simulations, we assume that the magnetic field has a strong toroidal component and a vertical component in addition to a stochastic component. We obtain the self-similar solutions to the equations describing the magnetized ADAFs, taking into account the dynamical effects of the outflow. We compare the results with the canonical ADAFs and find that the dynamical properties of ADAFs such as radial velocity, angular velocity and temperature can be significantly changed in the presence of ordered magnetic fields and outflows. The stronger the magnetic field is, the lower the temperature of the accretion flow will be and the faster the flow rotates. The relevance to observations is briefly discussed.     

On the origin of the 511-keV emission in the Galactic Centre

Diffuse 511-keV line emission, from the annihilation of cold positrons, has been observed in the direction of the Galactic Centre for more than 30 yr. The latest high-resolution maps of this emission produced by the SPI instrument on INTEGRAL suggest at least one component of the emission is spatially coincident with the distribution of ∼70 luminous, low-mass X-ray binaries detected in the soft gamma-ray band. The X-ray band, however, is generally a more sensitive probe of X-ray binary populations. Recent X-ray surveys of the Galactic Centre have discovered a much larger population (>4000) of faint, hard X-ray point sources. We investigate the possibility that the positrons observed in the direction of the Galactic Centre originate in pair-dominated jets generated by this population of fainter accretion-powered X-ray binaries. We also consider briefly whether such sources could account for unexplained diffuse emission associated with the Galactic Centre in the microwave (the Wilkinson Microwave Anisotropy Probe 'haze') and at other wavelengths. Finally, we point out several unresolved problems in associating Galactic Centre 511-keV emission with the brightest X-ray binaries.     

Delayed X-ray emission from fallback in compact-object mergers

When double neutron star or neutron star–black hole binaries merge, the final remnant may comprise a central solar-mass black hole surrounded by a  ∼0.01–0.1 M_⊙  torus. The subsequent evolution of this disc may be responsible for short γ-ray bursts (SGRBs). A comparable amount of mass is ejected into eccentric orbits and will eventually fallback to the merger site after ∼0.01 s. In this paper, we investigate analytically the fate of the fallback matter, which may provide a luminous signal long after the disc is exhausted. We find that matter in the eccentric tail returns at a super-Eddington rate and eventually (≳0.1 s) is unable to cool via neutrino emission and accrete all the way to the black hole. Therefore, contrary to previous claims, our analysis suggests that fallback matter is not an efficient source of late-time accretion power and unlikely to cause the late-flaring activity observed in SGRB afterglows. The fallback matter rather forms a radiation-driven wind or a bound atmosphere. In both the cases, the emitting plasma is very opaque and photons are released with a degraded energy in the X-ray band. We therefore suggest that compact binary mergers could be followed by an 'X-ray renaissance', as late as several days to weeks after the merger. This might be observed by the next generation of X-ray detectors.     

Black hole mergers: can gas discs solve the 'final parsec' problem?

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre-existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework.
The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high-photon-energy variability on the relatively long orbital time-scale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary.
We also find that the disc can drive the binary to merger on a reasonable time-scale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e.   a ₀≲ 0.05  pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is ≲0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem for a significant class of mergers.     

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.