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.     

Steady shocks around black holes produced by sub-Keplerian flows with negative energy

We discuss a special case of formation of axisymmetric shocks in the accretion flow of ideal gas on to a Schwarzschild black hole: when the total energy of the flow is negative. The result of our analysis enlarges the parameter space for which these steady shocks are exhibited in the accretion of gas rotating around relativistic stellar objects. Since Keplerian discs have negative total energy, we guess that, in this energy range, the production of the shock phenomenon might be easier than in the case of positive energy. So our outcome reinforces the view that sub-Keplerian flows of matter may significantly affect the physics of the high energy radiation emission from black hole candidates. We give a simple procedure to obtain analytically the position of the shocks. The comparison of the analytical results with the data of one-dimensional (1D) and two-dimensional (2D) axisymmetric numerical simulations confirms that the shocks form and are stable.     

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.     

Shock oscillation model for quasi-periodic oscillations in stellar mass and supermassive black holes

We numerically examine centrifugally supported shock waves in 2D rotating accretion flows around a stellar mass  (10 M_⊙)  and a supermassive  (10⁶ M_⊙)  black holes over a wide range of input accretion rates of     . The resultant 2D shocks are unstable with time and the luminosities show quasi-periodic oscillations (QPOs) with modulations of a factor of 2–3 and with periods of a tenth of a second to several hours, depending on the black hole masses. The shock oscillation model may explain the intermediate frequency QPOs with 1–10 Hz observed in the stellar mass black hole candidates and also suggest the existence of QPOs with the period of hours in active galactic nuclei. When the accretion rate     is low, the luminosity increases in proportion to the accretion rate. However, when     greatly exceeds the Eddington critical rate     , the luminosity is insensitive to the accretion rate and is kept constantly around  ∼3 L _E  . On the other hand, the mass-outflow rate     increases in proportion to     and it amounts to about a few per cent of the input mass-flow rate.     

The runaway instability of thick discs around black holes I. The constant angular momentum case

We present results from a numerical study of the runaway instability of thick discs around black holes. This instability is an important issue for most models of cosmic gamma-ray bursts, where the central engine responsible for the initial energy release is such a system consisting of a thick disc surrounding a black hole. We have carried out a comprehensive number of time-dependent simulations aimed at exploring the appearance of the instability. Our study has been performed using a fully relativistic hydrodynamics code. The general relativistic hydrodynamic equations are formulated as a hyperbolic flux-conservative system and solved using a suitable Godunov-type scheme. We build a series of constant angular momentum discs around a Schwarzschild black hole. Furthermore, the self-gravity of the disc is neglected and the evolution of the central black hole is assumed to be that of a sequence of exact Schwarzschild black holes of varying mass. The black hole mass increase is thus determined by the mass accretion rate across the event horizon. In agreement with previous studies based on stationary models, we find that by allowing the mass of the black hole to grow the disc becomes unstable. Our hydrodynamical simulations show that for all disc-to-hole mass ratios considered (between 1 and 0.05), the runaway instability appears very fast on a dynamical time-scale of a few orbital periods, typically a few 10 ms and never exceeding 1 s for our particular choice of the mass of the black hole (2.5 M) and a large range of mass fluxes  ( m 10^-3 M s^-1)  . The implications of our results in the context of gamma-ray bursts are briefly discussed.     

The evolution of misaligned accretion discs and spinning black holes

In this paper, we consider the process of alignment of a spinning black hole and a surrounding misaligned accretion disc. We use a simplified set of equations, that describe the evolution of the system in the case where the propagation of warping disturbances in the accretion disc occurs diffusively, a situation likely to be common in the thin discs in active galactic nuclei (AGN). We also allow the direction of the hole spin to move under the action of the disc torques. In such a way, the evolution of the hole–disc system is computed self-consistently. We consider a number of different situations and we explore the relevant parameter range, by varying the location of the warp radius R _w and the propagation speed of the warp. We find that the dissipation associated with the twisting of the disc results in a large increase in the accretion rate through the disc, so that AGN accreting from a misaligned disc are likely to be significantly more luminous than those accreting from a flat disc. We compute explicitly the time-scales for the warping of the disc and for the alignment process and compare our results with earlier estimates based on simplified steady-state solutions. We also confirm earlier predictions that, under appropriate circumstances, accretion can proceed in a counter-aligned fashion, so that the accreted material will spin-down the hole, rather than spinning it up. Our results have implication in a number of different observational features of AGN such as the orientation and shape of jets, the shape of X-ray iron lines and the possibility of obscuration and absorption of X-ray by the outer disc as well as the general issue of the spin history of black holes during their growth.     

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 structure of black hole magnetospheres I. Schwarzschild black holes

We introduce a multipolar scheme for describing the structure of stationary, axisymmetric, force-free black hole magnetospheres in the '3+1' formalism. We focus here on Schwarzschild spacetime, giving a complete classification of the separable solutions of the stream equation. We show a transparent term-by-term analogy of our solutions with the familiar multipoles of flat-space electrodynamics. We discuss electrodynamic processes around disc-fed black holes in which our solutions find natural applications: (i) 'interior' solutions in studies of the BlandfordZnajek process of extracting the rotational energy of holes, and of the formation of relativistic jets in active galactic nuclei and 'microquasars'; (ii) 'exterior' solutions in studies of accretion disc dynamos, disc-driven winds and jets. On the strength of existing numerical studies, we argue that the poloidal field structures found here are also expected to hold with good accuracy for rotating black holes, except for the cases of the maximum possible rotation rates. We show that the closed-loop exterior solutions found here are not in contradiction with the MacdonaldThorne theorem, as these solutions, which diverge logarithmically on the horizon of the hole , only apply to those regions that exclude .