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.     

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.     

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.     

Magnetic field dragging and the vertical structure of thin accretion discs

The presence of an imposed vertical magnetic field may drastically influence the structure of thin accretion discs. If the field is sufficiently strong, the rotation law can depart from the Keplerian one. We consider the structure of a disc for a given eddy magnetic diffusivity but neglect details of the energy transport. The magnetic field is assumed to be in balance with the internal energy of the accretion flow. The thickness of the disc as well as the turbulent magnetic Prandtl number and the viscosity, α , are the key parameters of our model. The calculations show that the radial velocity can reach the sound speed for a magnetic disc if the thickness is comparable to that of a non-magnetic one. This leads to a strong amplification of the accretion rate for a given surface density. The inclination angle of the magnetic field lines can exceed the critical value 30° (required to launch cold jets) even for a relatively small magnetic Prandtl number of order unity. The toroidal magnetic fields induced at the disc surface are smaller than predicted in previous studies.     

Critical properties of spherically symmetric accretion in a fractal medium

Spherically symmetric transonic accretion of a fractal medium has been studied in both the stationary and the dynamic regimes. The stationary transonic solution is greatly sensitive to infinitesimal deviations in the outer boundary condition, but the flow becomes transonic and stable when its evolution is followed through time. The evolution towards transonicity is more pronounced for a fractal medium than it is for a continuum, and in the former case the static sonic condition is met on relatively larger length scales. The dynamic approach also shows that there is a remarkable closeness between an equation  of motion for a perturbation in the flow, and the metric of an analogue acoustic black hole. The stationary inflow solutions of a fractal medium are as much stable under the influence of linearized perturbations as they are for the fluid continuum.     

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.     

Critical properties of spherically symmetric black hole accretion in Schwarzschild geometry

The stationary, spherically symmetric, polytropic and inviscid accretion flow in the Schwarzschild metric has been set-up as an autonomous first-order dynamical system, and it has been studied completely analytically. Of the three possible critical points in the flow, the one that is physically realistic behaves like the saddle point of the standard Bondi accretion problem. One of the two remaining critical points exhibits the strange mathematical behaviour of being either a saddle point or a centre-type point, depending on the values of the flow parameters. The third critical point is always unphysical and behaves like a centre-type point. The treatment has been extended to pseudo-Schwarzschild flows for comparison with the general relativistic analysis.     

Perturbation of mass accretion rate,associated acoustic geometry and stability analysis

We investigate the stability of stationary integral solutions of an ideal irrotational fluid in a general static and spherically symmetric background, by studying the profile of the perturbation of the mass accretion rate. We consider low angular momentum axisymmetric accretion flows for three different accretion disk models and consider time dependent and radial linear perturbation of the mass accretion rate. First we show that the propagation of such perturbation can be determined by an effective 2 × 2 matrix, which has qualitatively similar acoustic causal properties as one obtains via the perturbation of the velocity potential. Next, using this matrix we analytically address the stability issues, for both standing and travelling wave configurations generated by the perturbation. Finally, based on this general formalism we briefly discuss the explicit example of the Schwarzschild spacetime and compare our results of stability with the existing literature, which instead address this problem via the perturbation of the velocity potential.     

On some transonic aspects of general relativistic spherical accretion on to Schwarzschild black holes

The equations governing general relativistic, spherically symmetric, hydrodynamic accretion of polytropic fluid on to black holes are solved in the Schwarzschild metric to investigate some of the transonic properties of the flow. Only stationary solutions are discussed. For such accretion, it has been shown that real physical sonic points may form even for flow with   γ <4/3  or   γ >5/3  . The behaviour of some flow variables in the close vicinity of the event horizon is studied as a function of specific energy and the polytropic index of the flow.     

The excitation, propagation and dissipation of waves in accretion discs: the non-linear axisymmetric case

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing, and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic  ( γ =5/3)  . Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance  ℳ_r≳0.01  undergo shocks within a distance of order the resonance radius.     

Sonic points and shocks in isothermal accretion and winds in Kerr geometry

This paper is devoted to the study of sonic points and shocks in stationary, axially symmetric, isothermal flows around a Kerr black hole. We first show the dependence of the location of the sonic point with the flow's angular momentum for different isothermal sound speeds. With our selected shock jump conditions, we then discuss the properties of the shock, including the location and the strength. The ambiguity regarding the shock locations is removed by stability analysis. We also find some differences between the shock in isothermal flows and that in adiabatic flows. Subject headings: accretion, accretion disks-black hole physics-hydrodynamics-relativity-shock waves.     

Vertically self-gravitating ADAFs in the presence of toroidal magnetic field

Force due to the self-gravity of the disc in the vertical direction is considered to study its possible effects on the structure of a magnetized advection-dominated accretion disc. We present steady-sate self similar solutions for the dynamical structure of such a type of the accretion flows. Our solutions imply reduced thickness of the disc because of the self-gravity. It also imply that the thickness of the disc will increase by adding the magnetic field strength.     

Atmospheric acoustics of Titan, Mars, Venus, and Earth

Planetary atmospheres are complex dynamical systems whose structure, composition, and dynamics intimately affect the propagation of sound. Thus, acoustic waves, being coupled directly to the medium, can effectively probe planetary environments. Here we show how the acoustic absorption and speed of sound in the atmospheres of Venus, Mars, Titan, and Earth (as predicted by a recent molecular acoustics model) mirror the different environments. Starting at the surface, where the sound speed ranges from ∼200 m/s for Titan to ∼410 m/s for Venus, the vertical sound speed profiles reveal differences in the atmospheres' thermal layering and composition. The absorption profiles are relatively smooth for Mars, Titan, and Earth while Venus stands out with a noticeable attenuation dip occurring between 40 and 100 km. We also simulate a descent module sampling the sound field produced by a low-frequency “event” near the surface noting the occurrence of acoustic quiet zones.     

The radial–azimuthal instability of a hot two-temperature accretion disc with advection

The radial–azimuthal instability of a hot two-temperature accretion disc with advection is examined in this paper. We find that the inclusion of very little advection has significant effects on two acoustic modes for a geometrically thin, cooling-dominated two-temperature disc, but has no effect on acoustic modes for a geometrically slim, cooling-dominated two-temperature disc. We also find that, when azimuthal perturbations are considered, the stability properties of the disc are different from those in the pure radial perturbation case. An increase of the azimuthal wavenumber will stabilize the acoustic modes but make the viscous and thermal modes more unstable for a geometrically thin, cooling-dominated two-temperature disc. It makes the thermal mode more unstable and the acoustic mode more stable, but only affects the instability of the viscous mode for short-wavelength perturbations for a geometrically slim, cooling-dominated two-temperature disc. For a geometrically slim, advection-dominated two-temperature disc, the increase of the azimuthal perturbation makes the I- and O-modes more stable and the thermal mode more unstable, but has no effect on the viscous mode.     

Magnetic field bending in accretion discs with dynamos

We consider the problem of poloidal magnetic field advection and bending of an initially vertical field owing to radial inflow in thin accretion discs. For a ratio of kinematic viscosity to magnetic diffusivity of order unity, significant bending of an externally applied vertical field cannot occur in a disc with no internal dynamo. However, we show that if poloidal field is generated by a dynamo operating near its critical state, then significant field bending may be possible. Our results are of particular relevance to wind launching from accretion discs.     

Two- and three-dimensional numerical simulations of accretion discs in a close binary system

We perform 2D and 3D numerical simulations of an accretion disc in a close binary system using the simplified flux vector splitting (SFS) finite volume method. In our calculations, the gas is assumed to be ideal with γ =1.01, 1.05, 1.1 and 1.2 . The mass ratio of the mass-losing star to the mass-accreting star is unity. Our results show that spiral shocks are formed on the accretion disc in all cases. In 2D calculations we find that the smaller γ is, the more tightly the spiral winds. We observe this trend in 3D calculations as well in a somewhat weaker sense. Mach numbers in our discs are less than 10. These values are lower than the values in observed accretion discs in close binary systems.
Recently, Steeghs, Harlaftis & Horne found the first convincing evidence for spiral structure in the accretion disc of the eclipsing dwarf nova binary IP Pegasi, using the technique known as Doppler tomography. Although the Mach numbers in present calculations are rather low, we may claim that the spiral structure that we discovered in earlier numerical simulations is now found observationally.     

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.     

Three-Dimensional Simulations of Jets from Keplerian Disks: Stability Issues

Three-dimensional simulations of the time-dependent evolution of non-relativistic outflows from the surface of Keplerian accretion disks are presented. We investigate the outflow that arises from a magnetized accretion disk that is initially in hydrostatic balance with its surrounding cold corona. Our simulations show that jets maintain their long-term stability through a self-limiting process wherein the average Alfvénic Mach number within the jet is maintained to order unity. This is accomplished in at least two ways. First, poloidal magnetic field is concentrated along the central axis of the jet forming a `backbone' in which the Alfvén speed is sufficiently high to reduce the average jet Alfvénic Mach number to unity. Second, the onset of higher order Kelvin-Helmholtz `flute' modes (m ≥ 2) reduce the efficiency with which the jet material is accelerated, and transfer kinetic energy of the out flow into the stretched, poloidal field lines of the distorted jet. This too has the effect of increasing the Alfvén speed and thereby reducing the Alfvénic Mach number. The jet is able to survive the onset of the more destructive m=1 mode in this way.     

Fractal features in accretion discs

Fractal concepts have been introduced in the accretion disc as a new feature. Due to the fractal nature of the flow, its continuity condition undergoes modifications. The conserved stationary fractal flow admits only saddle points and centre-type points in its phase portrait. Completely analytical solutions of the equilibrium point conditions indicate that the fractal properties enable the flow to behave like an effective continuum of lesser density, and facilitate the generation of transonicity. However, strongly fractal flows inhibit multitransonicity from developing. The mass accretion rate exhibits a fractal scaling behaviour, and the entire fractal accretion disc is stable under linearized dynamic perturbations.