Model dependence of transonic properties of accretion flows around black holes

We analytically study how the behaviour of accretion flows changes when the flow model is varied. We study the transonic properties of the conical flow, a flow of constant height and a flow in vertical equilibrium, and show that all these models are basically identical, provided that the polytropic constant is suitably changed from one model to another. We show that this behaviour is extendible even when standing shocks are produced in the flow. The parameter space where shocks are produced remains roughly identical in all these models when the same transformation among the polytropic indices is used. We present applications of these findings.     

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.     

Mini-discs around spinning black holes

Accretion on to black holes in wind-fed binaries and in collapsars forms small rotating discs with peculiar properties. Such 'mini-discs' accrete on the free-fall time without the help of viscosity and nevertheless can have a high radiative efficiency. The inviscid mini-disc model was previously constructed for a non-rotating black hole. We extend the model to the case of a spinning black hole, calculate the structure and radiative efficiency of the disc and find their dependence on the black hole spin. If the angular momenta of the disc and the black hole are anti-aligned, a hydrodynamic analogue of Penrose process takes place.     

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.     

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.     

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.     

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.     

Global solutions of adiabatic accretion flows with isothermal shocks in Kerr black hole geometry

This paper presents global solutions of adiabatic accretion flows with isothermal shocks in Kerr black hole geometry. It is known that in the previously studied cases, where the flow including the shock is either entirely adiabatic or entirely isothermal, there can be no more than one stable shock solution, and the solution can only be of α –x type. However, the solution topology in the present case shows remarkable new characteristics: for the same flow parameters there can be two stable shock solutions satisfying physical boundary conditions, and the solution can be of three types, namely α– x , x –α and α–α type. In addition, shocks in the present case occur for a parameter region different from that for Rankine–Hugoniot shocks. These results greatly increase the possibilities of shock formation in astrophysical flows. It is also significant that the effects of frame-dragging of a rapid Kerr black hole on the shock formation are discovered. Finally, a brief comparison is made between shocked inviscid flows and two types of shock-free viscous flows, namely those of Shakura & Sunyaev and Narayan & Yi, and some comments are made about the fact that numerous authors who have studied transonic global solutions of accretion flows have found no shocks.     

Warp and eccentricity propagation in discs around black holes

We consider the inward propagation of warping and eccentric disturbances in discs around black holes under a wide variety of conditions. In our calculations, we use secular theories of warped and eccentric discs and assume the deformations to be stationary and propagating in a disc model similar to regions (a) and (b) of Shakura & Sunyaev discs. We find that the propagation of deformations to the innermost regions of the disc is facilitated for low viscous damping and high accretion rate. We relate our results to the possible excitation of trapped inertial modes, and to the observations of high-frequency quasi-periodic oscillations (QPOs) in black hole systems in the very high spectral state.     

Super-Eddington accretion discs around Kerr black holes

We calculate the structure of the accretion disc around a rapidly rotating black hole with a super-Eddington accretion rate. The luminosity and height of the disc are reduced by the advection effect. In the case of large viscosity parameter, α>0.03, the accretion flow deviates strongly from thermodynamic equilibrium and overheats in the central region. With increasing accretion rate, the flow temperature steeply increases, reaches maximum, and then falls off. The maximum is achieved in the advection-dominated regime of accretion. The maximum temperature in the disc around a massive black hole of M =10⁸ M⊙ with α=0.3 is of order 3×10⁸ K. The discs with large accretion rates can emit X-rays in quasars as well as in galactic black hole candidates.     

A massive thick disc around a massive black hole and its runaway instability

We investigate the runaway instability of configurations consisting of a massive dense but non-self-gravitating thick disc/torus surrounding a massive black hole (MBH). We limit our model parameters to values that result in a self-consistent thick disc around an MBH. We identify, analytically, the index of the angular momentum distribution that will form a thick disc as an initial equilibrium state, and obtain the mass ratio of the disc to the central black hole for which the disc is dominated by the radiation pressure. In our theoretical framework we find that a self-consistent thick disc with constant angular momentum leads to a runaway instability on a dynamical time-scale. However, even a slight increase of the specific angular momentum outwards has a strong stabilizing effect on the accretion process. Finally, we discuss our results and present possible applications to high-energy emission.     

Studies of accretion flows around rotating black holes – II. Standing shocks in the pseudo-Kerr geometry

Standing, propagating or oscillating shock waves are common in accretion and winds around compact objects. We study the topology of all possible solutions using the pseudo-Kerr geometry. We present the parameter space spanned by the specific energy and angular momentum and compare it with that obtained from the full general relativity to show that the potential can work satisfactorily in fluid dynamics also, provided the polytropic index is suitably modified. We then divide the parameter space depending on the nature of the solution topology. We specifically study the nature of the standing Rankine–Hugoniot shocks. We also show that as the Kerr parameter is increased, the shock location generally moves closer to the black hole. In future, these solutions can be used as guidelines to test numerical simulations around compact objects.     

Smoothed particle hydrodynamic simulations of viscous accretion discs around black holes

Viscous Keplerian discs become sub-Keplerian close to a black hole since they pass through sonic points before entering into it. We study the time evolution of polytropic viscous accretion discs (both in one- and two-dimensional flows) using smoothed particle hydrodynamics. We discover that for a large region of the parameter space spanned by energy, angular momentum and polytropic index, when the flow viscosity parameter is less than a critical value, standing shock waves are formed. If the viscosity is very high then the shock wave disappears. In the intermediate viscosity, the disc oscillates very significantly in the viscous time-scale. Our simulations indicate that these centrifugally supported high density regions close to a black hole play an active role in the flow dynamics, and consequently, the radiation dynamics.     

Emission-line profiles from self-gravitating toroids around black holes

We have computed line profiles from self-gravitating toroids around black holes. The specific angular momentum of the toroids is assumed to be constant in space. The images of the toroids show peculiar features in the far sides of the black holes. Concerning the line profiles, the red wing extends to the very low frequency region because the location of the inner edge is rather near the event horizon of the black hole and consequently the velocity of the inner edges of the toroids can be greater than that of the Kepler discs.     

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.