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.     

New self-similar solutions of polytropic gas dynamics

We explore semicomplete self-similar solutions for the polytropic gas dynamics involving self-gravity under spherical symmetry, examine behaviours of the sonic critical curve and present new asymptotic collapse solutions that describe 'quasi-static' asymptotic behaviours at small radii and large times. These new 'quasi-static' solutions with divergent mass density approaching the core can have self-similar oscillations. Earlier known solutions are summarized. Various semicomplete self-similar solutions involving such novel asymptotic solutions are constructed, either with or without a shock. In contexts of stellar core collapse and supernova explosion, a hydrodynamic model of a rebound shock initiated around the stellar degenerate core of a massive progenitor star is presented. With this dynamic model framework, we attempt to relate progenitor stars and the corresponding remnant compact stars: neutron stars, black holes and white dwarfs.     

Virial shocks in galactic haloes?

We investigate the conditions for the existence of an expanding virial shock in the gas falling within a spherical dark matter halo. The shock relies on pressure support by the shock-heated gas behind it. When the radiative cooling is efficient compared with the infall rate, the post-shock gas becomes unstable; it collapses inwards and cannot support the shock. We find for a monatomic gas that the shock is stable when the post-shock pressure and density obey     . When expressed in terms of the pre-shock gas properties at radius r it reads as  ρ r Λ( T )/ u ³ < 0.0126  , where ρ is the gas density, u is the infall velocity and Λ( T ) is the cooling function, with the post-shock temperature   T ∝ u ²  . This result is confirmed by hydrodynamical simulations, using an accurate spheri-symmetric Lagrangian code. When the stability analysis is applied in cosmology, we find that a virial shock does not develop in most haloes that form before   z ∼ 2  , and it never forms in haloes less massive than a few  10¹¹ M_⊙  . In such haloes, the infalling gas is not heated to the virial temperature until it hits the disc, thus avoiding the cooling-dominated quasi-static contraction phase. The direct collapse of the cold gas into the disc should have non-trivial effects on the star formation rate and on outflows. The soft X-ray produced by the shock-heated gas in the disc is expected to ionize the dense disc environment, and the subsequent recombination would result in a high flux of Lα emission. This may explain both the puzzling low flux of soft X-ray background and the Lα emitters observed at high redshift.     

On the angular momentum transfer on to compact stars in binary systems

Results of three-dimensional numerical simulations of the gas transfer in close binary systems show that, in addition to the formation of a tidally induced spiral shock wave, it is also possible for accretion streams to be produced, having low specific angular momentum in a region close to the accreting star. These streams are mainly placed above the orbital disc but are also unevenly present in the equatorial plane. The relevance of such flows is related to formation of hot coronae or bulges in regions very close to the accretor centre. The eventual formation of such bulges and shock-heated flows is interesting in the context of advection-dominated solutions and for the explanation of spectral properties of the black hole candidates in binary systems.     

The post-explosion shock propagation in the central region of our Galaxy

Immediate consequences of nuclear explosions on the structure and physical state of a galactic disk are considered in this paper. Explosions in the nucleus of a Galaxy generate strong shock waves which, when propagating onward heat and condensing the gas, form thin dense ring-like gaseous features behind it. Such rings and dense gaseous complexes have been observed in the central region of the Galaxy. These features have been treated here as the remnants of galactic shocks generated by nuclear explosions. We have estimated the time elapsed since the corresponding explosion, the energy released by explosion and the initial temperature and the velocity of the shock wave thus generated. The cooling of the gas heated by strong shocks has also been considered. The time taken by shock-heated gas to cool to its original temperature has been estimated to be of the order of 10⁵ to 10⁶ yr, according to the initial shock temperature which is about 9×10⁶ K or 6.4×10⁷ K. The rate of emission of energy and the total amount of energy dissipated away in the form of radiation in the cooling process, have been calculated for different values of initial shocktemperatures and also for different field intensities. The high-energy radiation emitted in the cooling process is suggested here as a source for the heating of dust grains, which ultimately are radiated in the infrared spectrum. Thus, a part of the infrared radiation, as measured by many authors, in the central region of the Galaxy, may originate ultimately from the cooling of the shock-heated gas there.     

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.     

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.     

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 effect of a central supermassive black hole on gas fuelling

When a supermassive black hole exists in the centre of a galaxy, an additional inner Lindblad resonance (ILR) exists inside the usual ILRs. We study gas dynamics in a weakly barred potential with a central supermassive black hole by using 2D numerical simulations, and we investigate the effect of the additional ILR on the fuelling of gas into nuclear starburst regions or active galactic nuclei (AGNs). Our numerical results show that strong trailing spiral shocks are formed at the resonance region, and that the gas in the shock region is rapidly fuelled into a central region and makes a nuclear gas ring. As a result, a large amount of gas is concentrated in the nuclear region beyond the ILR in a dynamical time-scale.     

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.     

Identification of a blue object with a newly detected X-ray source

An approximate analytical solution for self-similar flow behind a spherically-symmetric magnetogasdynamic strong shock wave is investigated using the method of Laumbach and Probstein (1969). The total energy of flow is taken to be dependent on the shock radius obeying a power law. The shock is propagating into a perfect gas at rest with non-uniform density and magnetic field.     

Self-similar propagation of axisymmetric radiative magnetogasdynamic shocks with increasing energy,II

In this paper self-similar solutions have been investigated for the propagation of axisymmetric radiative gasdynamic shocks caused by an explosion into an inhomogeneous ideal gas permeated by a current free azimuthal magnetic field. The effects of radiation flux and magnetic field together have been seen in the region of interest on the other flow variables. The total energy of the flow between the inner expanding surface and the shock is taken to be dependent on shock radius obeying a power law. The radiative pressure and energy have been neglected.     

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.     

General polytropic magnetofluid under self-gravity: Voids and shocks

We study the self-similar magnetohydrodynamics (MHD) of a quasi-spherical expanding void (viz. cavity or bubble) surrounding the centre of a self-gravitating gas sphere with a general polytropic equation of state. We show various analytic asymptotic solutions near the void boundary in different parameter regimes and obtain the corresponding void solutions by extensive numerical explorations. We find novel void solutions of zero density on the void boundary. These new void solutions exist only in a general polytropic gas and feature shell-type density profiles. These void solutions, if not encountering the magnetosonic critical curve (MCC), generally approach the asymptotic expansion solution far from the central void with a velocity proportional to radial distance. We identify and examine free-expansion solutions, Einstein–de Sitter expansion solutions, and thermal-expansion solutions in three different parameter regimes. Under certain conditions, void solutions may cross the MCC either smoothly or by MHD shocks, and then merge into asymptotic solutions with finite velocity and density far from the centre. Our general polytropic MHD void solutions provide physical insight for void evolution, and may have astrophysical applications such as massive star collapses and explosions, shell-type supernova remnants and hot bubbles in the interstellar and intergalactic media, and planetary nebulae.     

Approximate solutions for isothermal flows behind strong spherical shocks with variable energy

Self-similar unsteady flows with zero temperature gradient behind strong spherical shocks propagating in non-uniform perfect gas at rest are investigated. The total energy of the flow is assumed to be varying with the shock radius obeying a power law. Approximate solutions in a closed analytical form are obtained using the integral method. Also these solutions are shown to be useful to calculate easily and quickly the shock temperature, X-ray surface brightness and luminosity which are important in astrophysical problems. It is found that these approximate solutions are in close agreement with numerical solutions.     

The influence of the Mach number on the stability of radiative shocks

We study the stability properties of hydrodynamic shocks with finite Mach numbers. The linear analysis supplements previous analyses which took the strong shock limit. We derive the linearized equations for a general specific heat ratio as well as temperature and density power-law cooling functions, corresponding to a range of conditions relevant to interstellar atomic and molecular cooling processes. Boundary conditions corresponding to a return to the upstream temperature  ( R = 1)  and to a cold wall  ( R = 0)  are investigated. We find that for Mach number   M > 5  , the strong shock overstability limits are not significantly modified. For   M < 3  , however, shocks are considerably more stable for most cases. In general, as the shock weakens, the critical values of the temperature power-law index (below which shocks are overstable) are reduced for the overtones more than for the fundamental, which signifies a change in basic behaviour. In the   R = 0  scenario, however, we find that the overstability regime and growth rate of the fundamental mode are increased when cooling is under local thermodynamic equilibrium. We provide a possible explanation for the results in terms of a stabilizing influence provided downstream but a destabilizing effect associated with the shock front. We conclude that the regime of overstability for interstellar atomic shocks is well represented by the strong shock limit unless the upstream gas is hot. Although molecular shocks can be overstable to overtones, the magnetic field provides a significant stabilizing influence.     

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.     

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.     

Self-similar problem with increasing energy in radiative magneto-gas dynamics

In the present paper self-similar solutions have been investigated for the propagation of piston driven, radiative gas-dynamic shocks into an inhomogeneous ideal gas permeated by a current free azimuthal magnetic field for spherical symmetry. The effects of radiation flux and magnetic field together have been seen in the region of interest on the other flow variables. The total energy of the flow between the piston and the shock is taken to be dependent on the shock radius obeying a power law. The radiative pressure and energy have been neglected. This problem is more general than the others done so far. The word piston implies some means to drive plasma radially onwards.     

Dynamic evolution of a quasi-spherical general polytropic magnetofluid with self-gravity

In various astrophysical contexts, we analyze self-similar behaviours of magnetohydrodynamic (MHD) evolution of a quasi-spherical polytropic magnetized gas under self-gravity with the specific entropy conserved along streamlines. In particular, this MHD model analysis frees the scaling parameter n in the conventional polytropic self-similar transformation from the constraint of n+γ=2 with γ being the polytropic index and therefore substantially generalizes earlier analysis results on polytropic gas dynamics that has a constant specific entropy everywhere in space at all time. On the basis of the self-similar nonlinear MHD ordinary differential equations, we examine behaviours of the magnetosonic critical curves, the MHD shock conditions, and various asymptotic solutions. We then construct global semi-complete self-similar MHD solutions using a combination of analytical and numerical means and indicate plausible astrophysical applications of these magnetized flow solutions with or without MHD shocks.