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.     

The role of flow geometry in influencing the stability criteria for low angular momentum axisymmetric black hole accretion

Using mathematical formalism borrowed from dynamical systems theory, a complete analytical investigation of the critical behaviour of stationary flows in low angular momentum axisymmetric black hole accretion, provides significant insight about the nature of the phase trajectories corresponding to transonic accretion in the steady state, without taking recourse to any explicit numerical method commonly reported in the literature on multi-transonic black hole accretion discs and related astrophysical phenomena. Investigation of an accretion process around a non-rotating black hole, forming different geometrical configurations of the flow structure under the influence of various pseudo-Schwarzschild potentials, reveals that the general profile of the parameter space divisions describing multi-critical accretion, is roughly equivalent for various flow geometries. However, a mere variation of the polytropic index of the flow cannot map a critical solution from one flow geometry to another, since the numerical domain of the parameter space responsible for producing multi-critical accretion does not undergo a continuous transformation in multi-dimensional parameter space. The stationary configuration used to demonstrate the aforementioned findings is shown to be stable under time-dependent linearised perturbations for all kinds of flow geometries, driven by any pseudo-Schwarzschild potential, and using a standard equation of state. Finally, the structure of the acoustic metric corresponding to the propagation of the linear perturbation is discussed for various flow geometries used.     

Variation of the gas and radiation content in the sub-Keplerian accretion disk around black holes and its impact to the solutions

We investigate the variation of the gas and the radiation pressure in accretion disks during the infall of matter to the black hole and its effect to the flow. While the flow far away from the black hole might be non-relativistic, in the vicinity of the black hole it is expected to be relativistic behaving more like radiation. Therefore, the ratio of gas pressure to total pressure (β) and the underlying polytropic index (γ) should not be constant throughout the flow. We obtain that accretion flows exhibit significant variation of β and then γ, which affects solutions described in the standard literature based on constant β. Certain solutions for a particular set of initial parameters with a constant β do not exist when the variation of β is incorporated appropriately. We model the viscous sub-Keplerian accretion disk with a nonzero component of advection and pressure gradient around black holes by preserving the conservations of mass, momentum, energy, supplemented by the evolution of β. By solving the set of five coupled differential equations, we obtain the thermo-hydrodynamical properties of the flow. We show that during infall, β of the flow could vary up to ∼300%, while γ up to ∼20%. This might have a significant impact to the disk solutions in explaining observed data, e.g. super-luminal jets from disks, luminosity, and then extracting fundamental properties from them. Hence any conclusion based on constant γ and β should be taken with caution and corrected.     

Mass and spin co-evolution during the alignment of a black hole in a warped accretion disc

In this paper, we explore the gravitomagnetic interaction of a black hole (BH) with a misaligned accretion disc to study BH spin precession and alignment jointly with BH mass M _BH and spin parameter a evolution, under the assumption that the disc is continually fed, in its outer region, by matter with angular momentum fixed on a given direction     . We develop an iterative scheme based on the adiabatic approximation to study the BH–disc co-evolution: in this approach, the accretion disc transits through a sequence of quasi-steady warped states (Bardeen–Petterson effect) and interacts with the BH until the spin   J _BH  aligns with     . For a BH aligning with a corotating disc, the fractional increase in mass is typically less than a few per cent, while the spin modulus can increase up to a few tens of per cent. The alignment time-scale     is of  ∼10⁵–10⁶ yr  for a maximally rotating BH accreting at the Eddington rate. BH–disc alignment from an initially counter-rotating disc tends to be more efficient compared to the specular corotating case due to the asymmetry seeded in the Kerr metric: counter-rotating matter carries a larger and opposite angular momentum when crossing the innermost stable orbit, so that the spin modulus decreases faster and so the relative inclination angle.