Black hole spin dependence of general relativistic multi-transonic accretion close to the horizon

We introduce a novel formalism to investigate the role of the spin angular momentum of astrophysical black holes in influencing the behavior of low angular momentum general relativistic accretion. We propose a metric independent analysis of axisymmetric general relativistic flow, and consequently formulate the space and time dependent equations describing the general relativistic hydrodynamic accretion flow in the Kerr metric. The associated stationary critical solutions for such flow equations are provided and the stability of the stationary transonic configuration is examined using an elegant linear perturbation technique. We examine the properties of infalling material for both prograde and retrograde accretion as a function of the Kerr parameter at extremely close proximity to the event horizon. Our formalism can be used to identify a new spectral signature of black hole spin, and has the potential of performing the black hole shadow imaging corresponding to the low angular momentum accretion flow.     

Standing shocks in magnetized dissipative accretion flow around black holes

We explore the global structure of the accretion flow around a Schwarzschild black hole where the accretion disc is threaded by toroidal magnetic fields. The accretion flow is optically thin and advection dominated. The synchrotron radiation is considered to be the active cooling mechanism in the flow. With this, we obtain the global transonic accretion solutions and show that centrifugal barrier in the rotating magnetized accretion flow causes a discontinuous transition of the flow variables in the form of shock waves. The shock properties and the dynamics of the post-shock corona are affected by the flow parameters such as viscosity, cooling rate and strength of the magnetic fields. The shock properties are investigated against these flow parameters. We further show that for a given set of boundary parameters at the outer edge of the disc, accretion flow around a black hole admits shock when the flow parameters are tuned for a considerable range.     

Model dependence of the multi-transonic behaviour,stability properties and the corresponding acoustic geometry for accretion onto rotating black holes

Stationary, multi-transonic, integral solutions of hydrodynamic axisymmetric accretion onto a rotating black hole have been compared for different geometrical configurations of the associated accretion disc structures described using the polytropic as well as the isothermal equations of state. Such analysis is performed for accretion under the influence of generalised post Newtonian pseudo Kerr black hole potential. The variations of the stationary shock characteristics with black hole spin have been studied in details for all the disc models and are compared for the flow characterised by the two aforementioned equations of state. Using a novel linear perturbation technique it has been demonstrated that the aforementioned stationary solutions are stable, at least upto an astrophysically relevant time scale. It has been demonstrated that the emergence of the horizon related gravity like phenomena (the analogue gravity effects) is a natural consequence of such stability analysis, and the corresponding acoustic geometry embedded within the transonic accretion can be constructed for the propagation of the linear acoustic perturbation of the mass accretion rate. The analytical expression for the associated sonic surface gravity κ has been obtained self consistently. The variations of κ with the black hole spin parameter for all different geometric configurations of matter and for various thermodynamic equations of state have been demonstrated.     

Application of a relativistic accretion disc model to X-ray spectra of LMC X-1 and GRO J1655-40

We present a general relativistic accretion disc model and its application to the soft-state X-ray spectra of black hole binaries. The model assumes a flat, optically thick disc around a rotating Kerr black hole. The disc locally radiates away the dissipated energy as a blackbody. Special and general relativistic effects influencing photons emitted by the disc are taken into account. The emerging spectrum, as seen by a distant observer, is parametrized by the black hole mass and spin, the accretion rate, the disc inclination angle and the inner disc radius.
We fit the ASCA soft-state X-ray spectra of LMC X-1 and GRO J1655-40 by this model. We find that, having additional limits on the black hole mass and inclination angle from optical/UV observations, we can constrain the black hole spin from X-ray data. In LMC X-1 the constraint is weak, and we can only rule out the maximally rotating black hole. In GRO J1655-40 we can limit the spin much better, and we find 0.68 a 0.88 . Accretion discs in both sources are radiation-pressure dominated. We do not find Compton reflection features in the spectra of any of these objects.     

Spectral signatures of dissipative standing shocks and mass outflow in presence of Comptonization around a black hole

Accretion flows having positive specific energy are known to produce outflows and winds which escape to a large distance. According to Two Component Advective Flow (TCAF) model, centrifugal pressure dominated region of the flow just outside the black hole horizon, with or without shocks, acts as the base of this outflow. Electrons from this region are depleted due to the wind and consequently, energy transfer rate due to inverse Comptonization of low energy photons are affected. Specifically, it becomes easier to cool this region and emerging spectrum is softened. Our main goal is to show spectral softening due to mass outflow in presence of Compton cooling. To achieve this, we modify Rankine-Hugoniot relationships at the shock front when post-shock region suffers mass loss due to winds and energy loss due to inverse Comptonization. We solve two-temperature equations governing an accretion flow around a black hole which include Coulomb exchange between protons and electrons and other major radiative processes such as bremsstrahlung and thermal Comptonization. We then compute emitted spectrum from this post-shock flow. We also show how location of standing shock which forms outer boundary of centrifugal barrier changes with cooling. With an increase in disc accretion rate \((\dot{m}_{d})\) , cooling is enhanced and we find that the shock moves in towards the black hole. With cooling, thermal pressure is reduced, and as a result, outflow rate is decreased. We thus directly correlate outflow rate with spectral state of the disc.     

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.     

Slim accretion discs: a model for ADAF–SLE transitions

We numerically construct slim, global, vertically integrated models of optically thin, transonic accretion discs around black holes, assuming a regularity condition at the sonic radius and boundary conditions at the outer radius of the disc and near the black hole. In agreement with several previous studies, we find two branches of shock-free solutions, in which the cooling is dominated either by advection or by local radiation. We also confirm that the part of the accretion flow where advection dominates is in some circumstances limited in size: it does not extend beyond a certain outer limiting radius. New results found in our paper concern the location of the limiting radius and the properties of the flow near to it. In particular, we find that beyond the limiting radius the advective-dominated solutions match on to Shapiro, Lightman &38; Eardley (SLE) discs through a smooth transition region. Therefore, the full global solutions are shock-free and unlimited in size. There is no need to postulate an extra physical effect (e.g. evaporation) for triggering the ADAF–SLE transition. It occurs as a result of standard accretion processes described by the classic slim disc equations.     

Constraints on black hole spins with a general relativistic accretion disk corona model

The peaks in the spectra of the accretion disks surrounding massive black holes in quasars are in the far-UV or soft X-ray band, which are usually not observed. However, in the disk corona model, soft photons from the disk are Comptonized to high energy in the hot corona, and the hard X-ray spectra(luminosity and spectral shape) contain information on the incident spectra from the disk. The values of black hole spin parameter a*are inferred from the spectral fitting, which are spread over a large range, ~-0.94 to 0.998. We find that the inclination angles and mass accretion rates are well determined by the spectral fitting, but the results are sensitive to the accuracy of black hole mass estimates. No tight constraints on the black hole spins are achieved, if the uncertainties in black hole mass measurements are a factor of four,which are typical for the single-epoch reverberation mapping method. Recently, the accuracy of black hole mass measurement has been significantly improved to 0.2- 0.4 dex with the velocity resolved reverberation mapping method. The black hole spin can be well constrained if the mass measurement accuracy is50%. In the accretion disk corona scenario, a fraction of power dissipated in the disk is transported into the corona, and therefore the accretion disk is thinner than a bare disk for the same mass accretion rate,because the radiation pressure in the disk is reduced. We find that the thin disk approximation, H/R0.1,is still valid if 0.3 m 0.5, provided half of the dissipated power is radiated in the corona above the disk.     

Models for jet power in elliptical galaxies: a case for rapidly spinning black holes

The power of jets from black holes is expected to depend on both the spin of the black hole and the structure of the accretion disc in the region of the last stable orbit. We investigate these dependencies using two different physical models for the jet power: the classical Blandford–Znajek (BZ) model and a hybrid model developed by Meier. In the BZ case, the jets are powered by magnetic fields directly threading the spinning black hole while in the hybrid model, the jet energy is extracted from both the accretion disc as well as the black hole via magnetic fields anchored to the accretion flow inside and outside the hole's ergosphere. The hybrid model takes advantage of the strengths of both the Blandford–Payne and BZ mechanisms, while avoiding the more controversial features of the latter. We develop these models more fully to account for general relativistic effects and to focus on advection-dominated accretion flows (ADAFs) for which the jet power is expected to be a significant fraction of the accreted rest mass energy.
We apply the models to elliptical galaxies, in order to see if these models can explain the observed correlation between the Bondi accretion rates and the total jet powers. For typical values of the disc viscosity parameter  α∼ 0.04 –0.3  and mass accretion rates consistent with ADAF model expectations, we find that the observed correlation requires   j ≳ 0.9  ; that is, it implies that the black holes are rapidly spinning. Our results suggest that the central black holes in the cores of clusters of galaxies must be rapidly rotating in order to drive jets powerful enough to heat the intracluster medium and quench cooling flows.     

Study of accretion disk dynamics in presence of cooling

We examine the behaviour of accretion flow around a rotating black hole in presence of cooling. We obtain global flow solutions for various accretion parameters that govern the accreting flow. We show that standing isothermal shock wave may develop in such an advective accretion flow in presence of cooling. This shocked solution has observational consequences as it successfully provides the possible explanations of energy spectra as well as generation of outflows/jets of various galactic and extra-galactic black hole candidates. We study the properties of isothermal shock wave and find that it strongly depends on the cooling efficiency. We identify the region in the parameter space spanned by the specific energy and specific angular momentum of the flow for standing isothermal shock as a function of cooling efficiencies and find that parameter space gradually shrinks with the increase of cooling rates. Our results imply that accretion flow ceases to contain isothermal shocks when cooling is beyond its critical value.     

Transfer of energy and angular momentum in the magnetic coupling between a rotating black hole and the surrounding accretion disc

The transfer of energy and angular momentum in the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc is discussed based on a mapping relation derived by considering the conservation of magnetic flux with two basic assumptions: (i) the magnetic field on the horizon is constant, (ii) the magnetic field on the disc surface varies as a power law with the radial coordinate of the disc. The following results are obtained: (i) the transfer direction of energy and angular momentum between the BH and the disc depends on the position of a co-rotation radius relative to the MC region on the disc, which is eventually determined by the BH spin; (ii) the evolution characteristics of a rotating BH in the MC process without disc accretion are depicted in a parameter space, and a series of values of the BH spin are given to indicate the evolution characteristics; (iii) the efficiency of converting accreted mass into radiation energy of a BH–disc system is discussed by considering the coexistence of disc accretion and the MC process; (iv) the MC effects on disc radiation and the emissivity index are discussed and it is concluded that they are consistent with the recent XMM–Newton observation of the nearby bright Seyfert 1 galaxy MCG–6-30-15 with reference to a variety of parameters of the BH–disc system.     

Modelling the behaviour of accretion flows in X-ray binaries

We review how the recent increase in X-ray and radio data from black hole and neutron star binaries can be merged together with theoretical advances to give a coherent picture of the physics of the accretion flow in strong gravity. Both long term X-ray light curves, X-ray spectra, the rapid X-ray variability and the radio jet behaviour are consistent with a model where a standard outer accretion disc is truncated at low luminosities, being replaced by a hot, inner flow which also acts as the launching site of the jet. Decreasing the disc truncation radius leads to softer spectra, as well as higher frequencies (including quasi periodic oscillations, QPOs) in the power spectra, and a faster jet. The collapse of the hot flow when the disc reaches the last stable orbit triggers the dramatic decrease in radio flux, as well as giving a qualitative (and often quantitative) explanation for the major hard–soft spectral transition seen in black holes. The neutron stars are also consistent with the same models, but with an additional component due to their surface, giving implicit evidence for the event horizon in black holes. We review claims of observational data which conflict with this picture, but show that these can also be consistent with the truncated disc model. We also review suggested alternative models for the accretion flow which do not involve a truncated disc. The most successful of these converge on a similar geometry, where there is a transition at some radius larger than the last stable orbit between a standard disc and an inner, jet dominated region, with the X-ray source associated with a mildly relativistic outflow, beamed away from the disc. However, the observed uniformity of properties between black holes at different inclinations suggests that even weak beaming of the X-ray emission may be constrained by the data. After collapse of the hot inner flow, the spectrum in black hole systems can be dominated by the disc emission. Its behaviour is consistent with the existence of a last stable orbit, and such data can be used to estimate the black hole spin. By contrast, these systems can also show very different spectra at these high luminosities, in which the disc spectrum (and probably structure) is strongly distorted by Comptonization. The structure of the accretion flow becomes increasingly uncertain as the luminosity approaches (and exceeds) the Eddington luminosity, though there is growing evidence that winds may play an important role. We stress that these high Eddington fraction flows are key to understanding many disparate and currently very active fields such as ULX, Narrow Line Seyfert 1’s, and the growth of the first black holes in the Early Universe.     

Hysteresis effects and diagnostics of the shock formation in low angular momentum axisymmetric accretion in the Kerr metric

The secular evolution of the purely general relativistic low angular momentum accretion flow around a spinning black hole is shown to exhibit hysteresis effects. This confirms that a stationary shock is an integral part of such an accretion disc in the Kerr metric. The equations describing the space gradient of the dynamical flow velocity of the accreting matter have been shown to be equivalent to a first order autonomous dynamical systems. Fixed point analysis ensures that such flow must be multi-transonic for certain astrophysically relevant initial boundary conditions. Contrary to the existing consensus in the literature, the critical points and the sonic points are proved not to be isomorphic in general, they can form in a completely different length scales. Physically acceptable global transonic solutions must produce odd number of critical points. Homoclinic orbits for the flow possessing multiple critical points select the critical point with the higher entropy accretion rate, confirming that the entropy accretion rate is the degeneracy removing agent in the system. However, heteroclinic orbits are also observed for some special situation, where both the saddle type critical points of the flow configuration possesses identical entropy accretion rate. Topologies with heteroclinic orbits are thus the only allowed non-removable degenerate solutions for accretion flow with multiple critical points, and are shown to be structurally unstable. Depending on suitable initial boundary conditions, a homoclinic trajectory can be combined with a standard non-homoclinic orbit through an energy preserving Rankine-Hugoniot type of stationary shock, and multi-critical accretion flow then becomes truly multi-transonic. An effective Lyapunov index has been proposed to analytically confirm why certain class of transonic flow cannot accommodate shock solutions even if it produces multiple critical points.     

Is the Bardeen–Petterson effect responsible for the warping and precession in NGC 4258?

Strong evidence for the presence of a warped Keplerian accretion disc in NGC 4258 (M 106) has been inferred from the kinematics of water masers detected at subparsec scales. Assuming a power-law accretion disc and using constraints on the disc parameters derived from observational data, we have analysed the relativistic Bardeen–Petterson effect driven by a Kerr black hole as the potential physical mechanism responsible for the disc warping. We found that the Bardeen–Petterson radius is comparable to or smaller than the inner radius of the maser disc (independent of the allowed value for the black hole spin parameter). Numerical simulations for a wide range of physical conditions have shown that the evolution of a misaligned disc due to the Bardeen–Petterson torques usually produces an inner flat disc and a warped transition region with a smooth gradient in the tilt and twist angles. Since this structure is similar to that seen in NGC 4258, we propose that the Bardeen–Petterson effect may be responsible for the disc warping in this galaxy. We estimated the time-scale necessary for the disc inside of the Bardeen–Petterson radius to align with the black hole's equator, as a function of the black hole spin. Our results show that the Bardeen–Petterson effect can align the disc within a few billion years in the case of NGC 4258. Finally, we show that if the observed curvature of the outer anomalous arms in the galactic disc of NGC 4258 is associated with the precession of its radio jet/counterjet, then the Bardeen–Petterson effect can provide the required precession period.     

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.     

Dissipative accretion flows around a rotating black hole

We study the dynamical structure of a cooling dominated rotating accretion flow around a spinning black hole. We show that non-linear phenomena such as shock waves can be studied in terms of only three flow parameters, namely the specific energy     , the specific angular momentum (λ) and the accretion rate     of the flow. We present all possible accretion solutions. We find that a significant region of the parameter space in the     plane allows global accretion shock solutions. The effective area of the parameter space for which the Rankine–Hugoniot shocks are possible is maximum when the flow is dissipation-free. It decreases with the increase of cooling effects and finally disappears when the cooling is high enough. We show that shock forms further away when the black hole is rotating compared to the solution around a Schwarzschild black hole with identical flow parameters at a large distance. However, in a normalized sense, the flow parameters for which the shocks form around the rotating black holes are produced shocks closer to the black hole. The location of the shock is also dictated by the cooling efficiency in that higher the accretion rate     , the closer is the shock location. We believe that some of the high-frequency quasi-periodic oscillations may be due to the flows with higher accretion rate around the rotating black holes.     

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.     

A model of magnetically induced disc–corona for black hole binaries

We propose a model of magnetic connection (MC) of a black hole with its surrounding accretion disc based on large-scale magnetic field. The MC gives rise to transport of energy and angular momentum between the black hole and the disc, and the closed field lines pipe the hot matter evaporated from the disc, and shape it in the corona above the disc to form a magnetically induced disc–corona system, in which the corona has the same configuration as the large-scale magnetic field. We numerically solve the dynamic equations in the context of the Kerr metric, in which the large-scale magnetic field is determined by dynamo process and equipartition between magnetic pressure and gas pressure. Thus we can obtain a global solution rather than assuming the distribution of large-scale magnetic field beforehand. The main MC effects lie in three aspects. (1) The rotational energy of a fast-spinning black hole can be extracted, enhancing the dissipation in the accretion disc, (2) the closed field lines provide a natural channel for corona matter escaping from disc and finally falling into black hole and (3) the scope of the corona can be bounded by the conservation of magnetic flux. We simulate the high-energy spectra of this system by using Monte Carlo method, and find that the relative hardness of the spectra decreases as accretion rate or black hole spin a _* increases. We fit the typical X-ray spectra of three black hole binaries  (GRO J1655−40, XTE 1118+480 and GX 339−4)  in the low/hard or very high state.