Threshold Drop in Accretion Density if Dark Energy is Accreting onto a Supermassive Black Hole

TL;DR

This paper investigates how dark energy, modeled as viscous modified Chaplygin gas, affects accretion density profiles around supermassive black holes, revealing a threshold drop in accretion when viscosity and dark energy effects are considered.

Contribution

It introduces a novel analysis of viscous dark energy's impact on accretion densities near black holes using an alpha-disc model and compares results with observational data.

Findings

Wind strength increases with dark energy as accreting agent.

Accretion density exhibits a threshold drop with specific viscosity and dark energy effects.

Results align with existing observational data.

Abstract

Galactic structures are supposed to be formed out of dark matter clustering. Some examples of supermassive black holes in the central regions of high redshift galaxies say that the concerned supermassive black holes have completed their constructions in a time less than it generally should be. To justify such discrepancies, we are forced to model about existences of black hole mimickers and exotic phenomena acting near the supermassive black holes. Motivated by these we study the natures of exotic matters, especially dark energy near the black holes. We choose modified Chaplygin gas as dark energy candidate. Again, the descriptions of gravitational waves or the attenuations of them when they are tunnelling through cosmological distances help us to measure the shear viscosity of the medium through which the waves have been travelled. Delayed decaying models of dark matters also suggest that dark energy and viscosity may come up as a byproduct of such decays or interactions. We consider the viscous nature of the medium, i.e., the dark energy. To do so, we choose an alpha-disc model as proposed by Shakura and Sunyaev. We study the variations of densities through accretion and wind branches for a different amount of viscosity regulated by the Shakura-Sunyaev's alpha parameter, spin parameter and different properties of accreting fluids, viz, the properties of adiabatic fluid and modified Chaplygin gas. We compare these results with each other and some existing density profiles drawn from observational data-based simulations. We follow that our result supports the data observed till date. Specifically, we see the wind to get stronger for dark energy as accreting agent. Besides, we see the accretion to have a threshold drop if the viscosity is chosen along with the repulsive effects of dark energy.

Equations

$${\mathcal{F}}_g\left(\mathcal{X},j\right)=\frac{{\left({\mathcal{X}}^2-2j\sqrt{\mathcal{X}}+j^2\right)}^2}{{\mathcal{X}}^3\{\sqrt{\mathcal{X}}\left(\mathcal{X}-2\right)+j{\}}^2},$$

$$p_{MCG}={\alpha }_{MCG}{\rho }_{MCG}-\frac{{\beta }_{MCG}}{{\rho }_{MCG}^{n_{MCG}}},$$

$$\frac{d}{d\mathcal{X}}\left(\mathcal{X}u\Sigma \right)=0,$$

$$u\frac{du}{d\mathcal{X}}+\frac{1}{\rho }\frac{dp}{d\mathcal{X}}-\frac{{\lambda }^2}{{\mathcal{X}}^3}+{\mathcal{F}}_g\left(\mathcal{X},j\right)=0,$$

$$u\frac{d\lambda }{d\mathcal{X}}=\frac{1}{\mathcal{X}\Sigma }\frac{d}{d\mathcal{X}}\left[{\mathcal{X}}^2{\alpha }_{SS}\left(p+\rho u^2\right)h\left(\mathcal{X}\right)\right],$$

$$h\left(\mathcal{X}\right)=c_s\sqrt{\frac{\mathcal{X}}{{\mathcal{F}}_g}}.$$

$$c_s^2=\frac{\partial p_{MCG}}{\partial {\rho }_{MCG}}={\alpha }_{MCG}+\frac{{\beta }_{MCG}{n}_{MCG}}{{\rho }_{MCG}^{n_{MCG}+1}}.$$

$$2c_s\frac{d{c}_s}{d\mathcal{X}}=-\left(n_{MCG}+1\right)\frac{{\beta }_{MCG}{n}_{MCG}}{{\rho }^{n_{MCG}+2}}\frac{d\rho }{d\mathcal{X}}=-\left(n_{MCG}+1\right)\left(c_s^2-{\alpha }_{MCG}\right)\times \frac{1}{\rho }\frac{d\rho }{d\mathcal{X}}$$

$$\frac{1}{\rho }\frac{d\rho }{d\mathcal{X}}=-\frac{2c_s}{\left(n_{MCG}+1\right)\left(c_s^2-{\alpha }_{MCG}\right)}\frac{d{c}_s}{d\mathcal{X}}.$$

$$\frac{1}{\rho }\frac{d\rho }{d\mathcal{X}}=\frac{1}{\rho }\frac{dp}{d\rho }\frac{d\rho }{d\mathcal{X}}=c_s^2\left(\frac{1}{\rho }\frac{d\rho }{d\mathcal{X}}\right).$$

$$\frac{1}{{\rho }_{MCG}}\frac{d{p}_{MCG}}{d\mathcal{X}}=-\frac{2c_s^3}{\left(n_{MCG}+1\right)\left(c_s^2-{\alpha }_{MCG}\right)}\frac{d{c}_s}{d\mathcal{X}}=-\frac{2c_s}{\left(n_{MCG}+1\right)}\frac{d{c}_s}{d\mathcal{X}}-\frac{\alpha }{\left(n_{MCG}+1\right)}\frac{d}{d\mathcal{X}}\left\{\ln \left(c_s^2-{\alpha }_{MCG}\right)\right\}$$

$$=-\frac{1}{n_{MCG}+1}\frac{d}{d\mathcal{X}}\left(c_s^2\right)-\frac{{\alpha }_{MCG}}{n_{MCG}+1}\frac{d}{d\mathcal{X}}\{ln\left(c_s^2-{\alpha }_{MCG}\right)\}.$$

$$\dot{\mathcal{M}}=\Theta \rho c_s\frac{{\mathcal{X}}^{\frac{3}{2}}}{{\mathcal{F}}_g^{\frac{1}{2}}}u,$$

$$\frac{d{c}_s}{d\mathcal{X}}=\left(\frac{3}{2\mathcal{X}}-\frac{1}{2{\mathcal{F}}_g}\frac{d{\mathcal{F}}_g}{d\mathcal{X}}+\frac{1}{u}\frac{du}{d\mathcal{X}}\right)\left\{\frac{\left(n_{MCG}+1\right)c_s\left(c_s^2-{\alpha }_{MCG}\right)}{\left(1-n_{MCG}\right)c_s^2+{\alpha }_{MCG}\left(n_{MCG}+1\right)}\right\}$$

$$\begin{gathered} \frac{d\lambda }{d\mathcal{X}}=\frac{\mathcal{X}{\alpha }_{SS}}{u}[\frac{1}{2}\left(\frac{5}{\mathcal{X}}-\frac{1}{{\mathcal{F}}_g}\frac{d{\mathcal{F}}_g}{d\mathcal{X}}\right)\left\{\frac{\left(n_{MCG}+1\right){\alpha }_{MCG}-c_s^2}{n_{MCG}}+u^2\right\} \\ +2u\frac{du}{d\mathcal{X}}+\left\{\left(\frac{\left(n_{MCG}+1\right){\alpha }_{MCG}-c_s^2}{n_{MCG}}+u^2\right)\frac{1}{c_s}-\left(c_s^2+u^2\right)\left(\frac{1}{n_{MCG}+1}\frac{2c_s}{c_s^2-{\alpha }_{MCG}}\right)\right\}\frac{d{c}_s}{d\mathcal{X}}]. \end{gathered}$$

$$\frac{du}{d\mathcal{X}}=\frac{\frac{{\lambda }^2}{{\mathcal{X}}^3}-{\mathcal{F}}_g\left(\mathcal{X}\right)+\left(\frac{3}{\mathcal{X}}-\frac{1}{{\mathcal{F}}_g}\frac{d{\mathcal{F}}_g}{d\mathcal{X}}\right)\frac{c_s^4}{\{\left(1-n_{MCG}\right)c_s^2+{\alpha }_{MCG}\left(n_{MCG}+1\right)\}}}{u-\frac{2c_s^4}{u\{\left(1-n_{MCG}\right)c_s^2+{\alpha }_{MCG}\left(n_{MCG}+1\right)\}}}.$$

$$\mathcal{A}{\left(\frac{du}{d\mathcal{X}}\right)}_{\mathcal{X}={\mathcal{X}}_c}^2+\mathcal{B}{\left(\frac{du}{d\mathcal{X}}\right)}_{\mathcal{X}={\mathcal{X}}_c}+\mathcal{C}=0.$$

$$\frac{GMm}{R^2}=F_{grav}=F_{rad}=P_{rad}\kappa m=\frac{L}{c}\frac{1}{4\pi R^2}\kappa m=\frac{L}{c}\frac{1}{4\pi R^2}\frac{{\sigma }_T}{m_p}m$$

$$\dot{\mathcal{M}}=\frac{4\pi GM{m}_p}{ϵc{\sigma }_T}=7.03592\times 10^{-16}\times \left(10^7{M}_{\odot }\right)se{c}^{-1}.$$

$$\rho =\frac{2.285\times 10^{-21}\times \sqrt{\frac{{\mathcal{F}}_g}{{\mathcal{X}}^3}}}{u{c}_s}.$$

$$\frac{\rho \left(\mathcal{X}{r}_g\right)}{{\rho }_{crit}}=\frac{{\delta }_c}{\mathcal{X}{\left(1+\mathcal{X}\right)}^2},$$

$$\rho (\mathcal{X}{r}_g)={\rho }_{s}exp\left[-b_n\left\{(\mathcal{X}{)}^{\frac{1}{n}}-1\right\}\right],$$

Full text

**Keywords : Black Hole Accretion Disc, Density of Accreting Fluid, Dark Energy.

**

PACS: 95.30.Sf,95.36.+x, 95.35.+d, 98.80.Cq, 98.80.-k

1 Introduction

One of the very important predictions of Einstein’s General Relativity (GR hereafter) is about the existence of a particular type of compact objects, the name of which is popularly coined as “Black Holes” (BHs hereafter). BHs are nothing but such a dense object, the escape velocity from the surface of which is more than the speed of light in vacuum, $c$ (for simplest example of such objects known as Schwarzschild BH, calculated by Schwarzschild in 1916, $v_{escape}=\sqrt{\frac{2GM}{R^{2}}}\geq c$, i.e., $R\leq\frac{2GM}{c^{2}}$, $G$ is Newton’s gravitational constant, $M$ is the mass of the central gravitating object, $R$ is the radius of the object). There exist many indirect proofs of the existences of BHs drawn from the peculiar behaviours and high energy physics of the infalling matters flowing towards the steep gravitational well around the concerned BH. With the help of observations of the “Event Horizon Telescope”, we are able to figure out the first-ever image of a BH’s event horizon [1]. To state more correctly, what has actually been observed is the BH’s shadow (a combination of event horizon and gravitational lensing of photons) not the event horizon itself. In any case, allow us to bring to readers’ attentions on the work [2] where the authors have studied the possibility that the object observed might not be a Kerr BH but a more exotic object such as a superspinar, i.e., a BH spinning above the Kerr limit, which however would not necessarily leads to a naked singularity since quantum gravity effects could conspire to “cover” the singularity (in fact it is claimed that such objects arise quite naturally in string theory, and some have even claimed that a detection of superspinning BHs could be a “proof” of string theory).

This particular supermassive BH is located at the central region of the elliptic galaxy M87. The mass of this central BH, span of the accretion disc, the angular speed of the disc and accretion rate onto the BHare found to be roughly equal to $\left(6.5\pm 0.2_{stat}\pm 0.7_{phys}\right)\times 10^{9}M_{\odot},~{}~{}2.5\times 10^{4}AU$, of the order $10^{3}kms^{-1}$ and $0.1M_{\odot}yr^{-1}$ respectively. Besides, recent gravitational wave detections[3] of BH pair mergers can be treated to be proofs of supermassive BHs as well.

Galactic structures through all over the universe possess two common features: the presence of a supermassive BH (SMBH hereafter) at the center and dark matter(DM) embeddings and surroundings all over the galaxy. The causes of these two phenomena, however, are not entirely clear to us. We are able to observe SMBH of mass $\sim 800$ million $M_{\odot}$ even at redshift $z=7.54$ which is to be formed within a span less than a billion years only [4]. To justify such existences, alternative models of SMBHs are proposed. In the background of classical GR, different models for extended objects are considered, eg. gravastars and boson stars. Besides, existences of exotic cases like naked singularity etc are theoretically used to mimic SMBHs. Motivations of these concepts were to consider only the gravitational effects of the BH alternatives and to distinguish between the observational evidences drawn from the BHs and the BH mimickers. Authors of the reference [5] have shown, for Milky way, these two components’ effects can not be distinguished in the inward region of a sphere persuing radius $\sim 100AU$ distance from the center. If considerable amounts of pulsars are present in the vicinity of a SMBH, then only we might be able to distinguish between the BH shadow and the evidences of BH mimickers present near the BH. These researches strongly establish the presence of exotic matter/energy near the core regions of galaxies. The variation of densities for the model considered in the reference [5] is shown in figure (0). The density falls abruptly as we move far from the central region of the galaxy.

“Dark energy” (DE hereafter) is a nomenclature to a hypothetical fluid which permeates all over the cosmos, exerting negative pressure to accelerate the late time expansion of the universe[17]. Interactive DM-DE models also support that if DM clusters are found in galactic centers, DE will grow from delayed decay of the DM. Along with this, the bulk viscosity is also supposed to grow in these regions [6, 7].

Accretion flow of interstellar plasmas surrounding a relativistic star, i.e., a Schwarzschild singularity was first so ever studied by Michel [8]. Before Michel, Bondi had worked on accretion properties of more or less Newtonian type and onto a nonrelativistic celestial object. Bondi did not consider the transonic properties of accretion disc either [9]. Michel has tried to plot the radial inward velocity($v_{x}$) of accretion flow for adiabatic fluids and has derived the energy flux of the flow. He has assumed the total internal energy of the infalling matter to get radiated out just before the infall. In Michel’s work, the role of viscosity was neglected.

Shakura and Sunyaev [10] tried to estimate the quantity of infalling matter and to make simpler modelling of the angular momentum transport which builds up the thin disc-like structure. They have predicted that the efficiency of the magnetic field etc grown in a disc can not exceed the thermal energy of the matter which again can be quantified through the speed of sound ($c_{s}$) in the said accreting fluid. Compared to the scale of radial distance from the center, the turbulence is homogeneous and isotropic. Calculations on molecular viscosity, laminar flow properties and magnetic field-angular momentum transport connection led them to the replacement of tangential stresses by the simplified term: $-{\cal W}_{r\phi}\sim\alpha_{SS}\rho c_{s}^{2}$, where $\alpha_{SS}$ is the well known Shakura Sunyaev parameter and $\rho$ is the density of the accreting fluid. Again, Michel or Shakura-Sunyaev neglected the GR effects on the curvature near BH except for the truncation of the thin disc inside the innermost stable orbit. A complete GR view of accretion, for the first time, was established by Novikov and Thorne [11].

After we have done with the simplifications of the complex natures of accretion disc viscosity, the very next problem to get solved was to simplify the highly nonlinear GR equations pointed out by Novikov and Thorne. To consider a stationary flow ($\frac{\partial}{\partial t}\equiv 0$), i.e., to simplify the analysis we must replace the relativistic attracting force offered by the central BH with a pseudo Newtonian force. We also consider an axisymmetric system with coordinates $r$, $\phi$ and $z$ with $\frac{\partial}{\partial\phi}\equiv 0$. Besides, Shakura and Sunyaev’s work was for thin accretion disc which must not sustain anymore if the accretion rate is enough high. Paczynski and Witta [12] have taken a way out where flat disc was considered, inner edge of which is extended somewhere between marginally bound circular orbit and last stable circular orbit controlled under a pseudo Newtonian potential, $\psi=-\frac{GM}{r-r_{g}}$, $r_{g}=\frac{2GM}{c^{2}}$, $r$ is the radial distance from the center of the gravitating object. Later, this potential which was constructed for Schwarzschild metric was generalized for Kerr metric by B. Mukhopadhyay as [13]:

[TABLE]

where ${\cal X}=\frac{r}{r_{g}}$, $j$ is dimensionless specific angular velocity for rotating BH. As the rotational parameter, $j$ is incorporated in the pseudo-Newtonian potential, we can find different accretion studies [14] which analyzed the effects of the spin parameter on the accretion disc. Viscous accretion study with this potential is constructed by the same group of researchers [15] and analysis of sonic point, density, temperature and abundance of different light elements in the disc and the outflow are studied.

Besides, the scenario of ongoing studies of adiabatic accretion properties around BHs, Babichev E. et al have started the study [16] of DE accretion. There are several DE candidates which have been proposed and worked with to support to late-time cosmic acceleration. Among them, Modified Chaplygin Gas (MCG hereafter) is a particular one. Variations of its characterizing parameters turn MCG to mimic radiation era to phantom epoch [18, 19, 20]. The Equation of State (EoS hereafter) of MCG is given by

[TABLE]

where $\alpha_{MCG}$, $\beta_{MCG}$ and $n_{MCG}$ are MCG parameters, $\rho_{MCG}$ is the energy density of MCG. Like other DE candidates, different properties of accretion of MCG onto BHs have also been studied[21, 22, 23]. In literature, we can find works on Chaplygin gas accretion. Biswas, R. et. al. [21] have speculated that wind becomes stronger for DE accretion and matter is thrown out with a speed equal to that of light at a finite distance from the BH. This implies that the repulsive power of DE makes the accretion phenomena weaker. Rotation of the central engine, in addition, makes this fainting out process more efficient. This work was followed by Biswas, R. [22]. In this article, it was speculated that the nonviscous accretion’s density profile is weaker if MCG accretion is considered. On the other hand, the wind branch’s density is very high near to the BH. Recently, Dutta, S. and Biswas, R. [24] have shown that, if the viscosity is encountered, the wind branch achieves speed equal to that of light at a distance which is lesser than the nonviscous case. Again BH spin enhances this tendency. In this present letter, we will study the variations of density in accretion and wind branches. The variations of density of accreting exotic matter have studied by different authors. Density variations in the accretion of DM halo is a well-opted topic as DM, unlike DE, is preferable to form clusters[25, 26].

In this letter, we will consider viscous accretion of DE on SMBHs and study the density profiles of both accretion and wind flows. We have taken both the non-rotating and rotating BHs as the central engines. MCG is chosen as the DE model. Every result has been compared with corresponding adiabatic cases.

We wish to study the effects of viscosity, the negative pressure of DE and rotation on the accretion, precisely, on the density of the accretion and wind flow around SMBH. Study of such density profiles may reveal the natures of accretion when an exotic matter is intended to accrete in.

The next part of this letter comprises of the mathematical formulation of viscous/nonviscous flow and adiabatic/MCG gas upon rotating/non-rotating BHs. After that, in section 3, we will study the variations of density graphically. In the end, we will briefly discuss the results and conclude.

2 Mathematical Formulation of Viscous Dark Energy Accretion Model

To construct the mathematical model of the problem, we have taken the equation of continuity for cylindrical accretion model as

[TABLE]

where $\Sigma$ is vertically integrated density given by, $\Sigma=I_{c}\rho_{e}h\left({\cal X}\right),$ $I_{c}$= constant (related to equation of state of fluid)= $1$ for simplicity, $\rho_{e}$=density of the accreting fluid at equitorial plane, $h\left({\cal X}\right)$= half thickness of the disc. $u=u_{\cal X}=\frac{v_{\cal X}}{c}$, $v_{\cal X}$ is the radially inward speed of accretion. Next three equations will be the three components of Navier Stokes equation :

(a) Radial momentum balance equation:

[TABLE]

(b) Azimuthal momentum balance equation:

[TABLE]

where $\lambda$ is specific angular momentum.

Assuming the vertical equilibrium from the vertical component of Navier Stokes equation we get the expression for $h({\cal X})$ as

[TABLE]

Last four equations together will complete the present model for three dependent variables $u$, $c_{s}$ and $\lambda$. Differentiating equation (2) with respect to $\rho$ we have,

[TABLE]

On differentiation with ${\cal X}$ and using equation (7) we obtain

[TABLE]

which yields

[TABLE]

Now, we will try to replace the pressure gradient over density term, i.e., the second term of the equation (4), $\frac{1}{\rho}\frac{dp}{d{\cal X}}$, as

[TABLE]

Combining (8) and (9), and rearranging we have

[TABLE]

[TABLE]

Now integrating the equation (3), we get the mass conservation equation

[TABLE]

where $\Theta$ is geometrical constant.

Replacing the value of $\rho$ from equation (7) in (11) and differentiating the whole term we get a differential equation for $c_{s}$ as

[TABLE]

and from the equations (5) and (2) we get the specific angular momentum gradient as

[TABLE]

Thus we get all the required values to replace $\frac{1}{\rho}\frac{dp}{d{\cal X}}$ in the equation (4), hence we obtain the radial velocity gradient as

[TABLE]

It is clear to observe that the denominator turns zero as $u=\frac{\sqrt{2}c_{s}^{2}}{\left\{(1-n_{MCG})c_{s}^{2}+\alpha_{MCG}(1+n_{MCG})\right\}^{\frac{1}{2}}}$ and this value lies in the interval $[0,~{}1]$ for some chosen ${\cal X}={\cal X}_{c}$(as $c_{s}$ lies in $[0,~{}1]$). Vanishing denominator signifies an unphysical flow. But as we consider the flow to be a continuous one, the numerator should simultaneously vanish at ${\cal X}={\cal X}_{c}$.

Applying the L’Hospital’s rule in equation (14) at the critical point (say ${\cal X}={\cal X}_{c}$) we get a quadratic equation of $\frac{du}{d{\cal X}}$ in the form,

[TABLE]

Where

${\cal A}=1+\frac{1}{c_{sc}^{2}}-\frac{4\left(n_{MCG}+1\right)\left(c_{sc}^{2}-\alpha_{MCG}\right)}{{\left\{\left(1-n_{MCG}\right)c_{sc}^{2}+\alpha_{MCG}\left(n_{MCG}+1\right)\right\}}^{2}}$,

${\cal B}=\frac{4\lambda u}{{\cal X}_{c}^{3}}+\frac{2u_{c}c_{sc}\left(1-n_{MCG}\right)}{\left\{\left(1-n_{MCG}\right)c_{sc}^{2}+\alpha_{MCG}\left(n_{MCG}+1\right)\right\}}-\frac{2u_{c}\left(n_{MCG}+1\right)\left(c_{sc}^{2}-\alpha_{MCG}\right)}{{\left\{\left(1-n_{MCG}\right)c_{sc}^{2}+\alpha_{MCG}\left(n_{MCG}+1\right)\right\}}^{2}}\left\{\frac{3}{{\cal X}_{c}}-\frac{1}{{\cal F}_{g}}\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)_{{\cal X}={\cal X}_{c}}\right\}+\\ \left[\left(c_{sc}^{2}+u_{c}^{2}\right)\left(\frac{1}{n_{MCG}+1}\frac{2c_{sc}}{c_{sc}^{2}-\alpha_{MCG}}\right)-\left\{\frac{3}{{\cal X}_{c}}-\frac{1}{{\cal F}_{g}}\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)_{{\cal X}={\cal X}_{c}}\right\}\frac{4u_{c}}{c_{sc}{\left\{\left(1-n_{MCG}\right)c_{sc}^{2}+\alpha_{MCG}\left(n_{MCG}+1\right)\right\}}}-\left(\frac{\left(n_{MCG}+1\right)\alpha_{MCG}-c_{sc}^{2}}{n_{MCG}}+u^{2}\right)\frac{1}{c_{sc}}\right]\\ \times\left\{\frac{u_{c}\left(n_{MCG}+1\right)\left(c_{sc}^{2}-\alpha_{MCG}\right)}{2c_{sc}^{3}}\right\}$

and ${\cal C}={\cal D}+{\cal E}+{\cal F}$.

The values of ${\cal D}$, ${\cal E}$ and ${\cal F}$ are,

${\cal D}=\left[\left(c_{sc}^{2}+u_{c}^{2}\right)\left(\frac{1}{n_{MCG}+1}\frac{2c_{sc}}{c_{sc}^{2}-\alpha_{MCG}}\right)-\left\{\frac{3}{{\cal X}_{c}}-\frac{1}{{\cal F}_{g}}\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)_{{\cal X}={\cal X}_{c}}\right\}\frac{4u_{c}}{c_{sc}{\left\{\left(1-n_{MCG}\right)c_{sc}^{2}+\alpha_{MCG}\left(n_{MCG}+1\right)\right\}}}\right.\\ \left.-\left(\frac{\left(n_{MCG}+1\right)\alpha_{MCG}-c_{sc}^{2}}{n_{MCG}}+u^{2}\right)\frac{1}{c_{sc}}\right]\left\{\frac{3}{2{\cal X}_{c}}-\frac{1}{2{\cal F}_{g}}\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)_{{\cal X}={\cal X}_{c}}\right\}\left\{\frac{\left(n_{MCG}+1\right)c_{sc}\left(c_{sc}^{2}-\alpha_{MCG}\right)}{\left(1-n_{MCG}\right)c_{sc}^{2}+\alpha_{MCG}\left(n_{MCG}+1\right)}\right\}$,

${\cal E}=\left\{\frac{3}{{\cal X}_{c}}-\frac{1}{{\cal F}_{g}}\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)_{{\cal X}={\cal X}_{c}}\right\}\frac{u_{c}c_{sc}\left(1-n_{MCG}\right)}{{\left\{\left(1-n_{MCG}\right)c_{sc}^{2}+\alpha_{MCG}\left(n_{MCG}+1\right)\right\}}}+\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)_{{\cal X}={\cal X}_{c}}+\left\{\frac{1}{{\cal F}_{g}^{2}}\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)^{2}_{{\cal X}={\cal X}_{c}}-\frac{1}{{\cal F}_{g}}\left(\frac{d^{2}{\cal F}_{g}}{d{\cal X}^{2}}\right)_{{\cal X}={\cal X}_{c}}-\frac{3}{{\cal X}_{c}^{2}}\right\}\frac{u_{c}^{2}}{2}$,

${\cal F}=\frac{\lambda\alpha_{ss}}{{\cal X}_{c}^{2}u_{c}}\left\{\frac{1}{{\cal F}_{g}}\left(\frac{d{\cal F}_{g}}{d{\cal X}}\right)_{{\cal X}={\cal X}_{c}}-\frac{5}{{\cal X}_{c}}\right\}\left\{\frac{\left(n_{MCG}+1\right)\alpha_{MCG}-c_{sc}^{2}}{n_{MCG}}+u_{c}^{2}\right\}$.

Here $c_{sc}=c_{s}\left({\cal X}={\cal X}_{c}\right)$. We can obtain the accretion and wind fluid speeds by solving equations (12), (13) and ((14)) starting from the critical point with proper initial conditions. Initial relation and hence the values of fluid speed and sound speed may be obtained by solving $N\left({\cal X}_{c}\right)=D\left({\cal X}_{c}\right)=0$ in the equation (14) due to the requirement of transonic smooth flow. $\lambda({\cal X}_{c})=\lambda_{c}$ has been chosen artificially.

We will try to find out density variations through the accretion and wind profiles. To do that we will use the concept of Eddington mass accretion rate. To keep an accretion process running, the considered central object should possess a luminosity up to a maximum limit beyond which the radiation pressure of that object will overcome gravity so that material outside the object will be swiped away from it, rather than falling inwards. This maximum luminosity is known as Eddington luminosity. To obtain this limit, we will balance the gravitational force with radiation force to obtain

[TABLE]

where $\kappa$ is the opacity, $\sigma_{T}$ is Thompson scattering cross-section and $m_{p}$ is the mass of a proton. So $L_{Edd}=\frac{4\pi GMcm_{p}}{\sigma_{T}}$. If $\epsilon$ fraction of the mass accretion rate is radiated out as energy, i.e., $L=\epsilon\dot{{\cal M}}c^{2}$ and we get Eddington mass accretion rate for a $10^{7}M_{\odot}$ object as

[TABLE]

From the equation (11) we get the density in $gcm^{-3}$ as an equation of dimensionless quantities ${\cal X}$, $c_{s}$ and $u$ as

[TABLE]

3 Variations of Density in Accretion and Wind Branch : Solutions

The DM distribution in galaxies and clusters of galaxies is well described by $\Lambda$CDM numerical n-body simulations. This theory also discussed cosmological evolution and the formation of astrophysical objects. These numerical simulations produce radial density profiles over a few decades in a radius near universal DM halo [27, 28]. Generally, Navarro-Frenk-White (NFW)[29] profile supplies us a good fit to density profiles of equilibrium DM halo [30, 31, 32, 33, 34, 35, 36]. But some discrepancies do occur between simulations and observations. As an example, we can look up to the core problem [37] of N-body simulations which gives cuspy central density profiles [38] which can fit well by NFW profile, while observation provides core like centers in DM dominated systems like dwarf galaxies [39, 40].

Usually, radial density profile is supposed to have a monotonically changing logarithmic slope such as NFW or Einasto profile [54]. There are two different classes which exhibit non-monotonic changes in their density profile slopes [55]. The first class consists of the systems that have density oscillation defined through their distribution function $f(E)$, or differential energy distribution $N(E)$. Examples of this class are DARKexp profiles [49], Polytropic [50, 51], Isothermal sphere(IS) [52], King Profile [53]. On the other hand, the members of the second class are the systems which are the product of dynamical evolution, like observed and simulated galaxies and pure DM halos. NFW profile and Einasto profile are the examples of this class.

Here, in our letter, we only consider NFW profile and Einasto profile as these two profiles give a good fit with the density profile of DM in the center of galaxy or cluster of galaxies. Navarro, Frenk and White used high-resolution simulations to study the formation of CDM halos with mass spanning about $4th$ orders of magnitude, ranging from dwarf galaxy halos to rich galaxy clusters. Their work showed that the equilibrium density profile of CDM halos of all masses can be accurately fitted by the simple formula [29].

[TABLE]

where, $\rho_{crit}=\frac{3H^{2}}{8\pi G}$ is the critical density and $\delta_{c}$ is a characteristic (dimensionless) density. The NFW simulation suggests that the density profile of an isolated equilibrium halo can be specified by giving two parameters, the mass of the halo and the characteristic density of the halo. In some particular hierarchical model, these parameters are related in such a way that the characteristic density is proportional to the mean cosmic density. Thus the characteristic density reflects the density of the universe at the collapse time of the objects, which merge to form the halo core. Now using our model with the equation (17) we get the density $4.3\times 10^{-23}gcm^{-3}$ at a distance of six Schwarzschild radius distance from the center of the central engine. It is to be noted that the order of the density we get $(10^{-23})$ which is higher than the density of the universe (approximately of order $10^{-30}$).

Now we will study the Einasto profile. The formula for this profile is given by

[TABLE]

where $\rho_{s}=\rho_{crit}\times\delta_{c}$, $n$= polytropic index, $b_{n}$ is a parameter depending on polytropic index. Einasto profile implies that the spherically averaged density profile of DM haloes are not universal but depends on the slope of DM power spectrum. The haloes, which are slowly growing through merger and accretion have steeper and more centrally concentrated density profiles with extended outer envelopes. Whereas the faster-growing haloes have shallower inner profiles and steep outer ones. We can obtain the density at $x=6$ with this formula both for adiabatic case and MCG. For adiabatic case, the density is $4.38\times 10^{-21}gcm^{-3}$ and for MCG it is $1.40\times 10^{-21}gcm^{-3}$.

3.1 Graphical Interpretations of the Density Profiles

We will plot density vs radial distance. Fig 1.1.1.a & 1.1.2.a describes the density vs radial distance curves for adiabatic fluid without viscosity with NFW profiles and Einasto profile respectively. Besides, the concerned central engine is taken to be a non-rotating, i.e., a Schwarzschild one. The solid lines in the graphs are the plots of accretion density whereas the dashed curve is for the wind density. It is seen along with the accretion curve density increases as we go nearer to the BH. This is quite obvious due to the pressure felt in the steep gravitational field.

For wind, density at far is lower than that in the accretion curve and increases more rapidly than the accretion density and near the BH it is enough high. Near the BH, mass is piled up creating a throwback or outward shock for the infalling matter creating higher wind. This may raise the wind density near the BH. As we go far, outward throwing force decreases and density also decreases.

Fig 1.1.1.b & 1.1.2.b show the density vs radial distance curve for inviscid MCG accretion for non-rotating BH with NFW profiles and Einasto profile respectively. We see the wind is highly dense when it is near to BH. As it goes far from the BH, density falls very quickly. High density, even if for MCG produces less negative pressure. But as the density falls, negative pressure becomes high and the disc must have fainted.

Now we will plot the density vs radial distance curve for adiabatic flow with viscosity $\alpha_{ss}=10^{-4}$ in fig 1.2.1.a & 1.2.2.a for non-rotating BH and viscosity $\alpha_{ss}=10^{-2}$ in Fig 1.3.1.a & 1.3.2.a with NFW profiles and Einasto profile respectively. We see that the wind branch is less steeper, it has the highest value near the central engine and while we go far from the central engine the density decreasing spontaneously. We can also observe that the accretion branch changes its properties with increasing viscosity. At a certain distance, the density of accretion jumped over suddenly and the accretion disc no longer exists there.

In fig 1.3.1.b & 1.3.2.b we plot the same graph with viscosity $\alpha_{ss}=10^{-4}$ we show that the density of accretion has no notable change.

Fig 1.3.1.b & 1.3.2.b show the density vs radial distance curve for MCG accretion with viscosity $\alpha_{ss}=10^{-2}$ for non-rotating BH with NFW profiles and Einasto profile respectively. We see the wind is highly dense when it is near to BH. As it goes far from the BH, density falls very quickly. High density, even if for MCG produces less negative pressure. But as the density falls, negative pressure becomes high and the disc must have fainted.

In figure $2.1.1.a$ we have plotted density profile for non-viscous adiabatic accretion for a rotating BH with $j=0.5$. It is observed that for adiabatic accretion the wind density very near to BH is $10^{4}$ times lighter than the non-rotating case if the viscosity is not working. Even far from the BH, when we are thousand Schwarzchild radius distance apart from the BH, the wind density for adiabatic fluid is ten times lighter than the non-rotating case. However, the accretion density is more or less identical to the non-rotating case.

In $2.1.1.b$ we have plotted the density for accreting Chaplygin gas onto a rotating BH with $j=0.5$. No viscosity is considered. Contrary to the adiabatic case we see that the wind density is hundred times heavier than the non-rotating case when we are very near to the BH. Rotation makes the disc radius shorter. This points out to the fact that dark energy’s repulsive effect is inhibited by the rotation of BH.

Fig $2.1.2.a$ and $2.1.2.b$ are the curves for Einasto profiles. The basic nature match with $2.1.1.a$ and $2.1.1.b$ respectively.

In fig $2.2.1.a$ to $2.2.2.b$ we have plotted the density profiles for rotating BH ($j=0.5$). Viscosity started to act here ($\alpha_{ss}=10^{-4}$). The order of density in the unit $g/cm^{3}$ is equal to the non-viscous case. But for both the adiabatic and MCG case we see the wind branch ended very near to the BH (than the non-viscous rotating case). In addition, the accreting branch of the adiabatic case is formed to grow high as we go far from the BH (at a particular distance). For MCG if we go far from BH the accretion density increases abruptly at a certain point and then becomes constant. This signifies that the accreting MCG’s density is constant at far. Then as it goes near to the BH the density falls at a particular distance (say ${\cal X}_{fall}$). Then it increases very slowly. So if viscosity and BH rotation both are working as a catalyst and an inhibitor, not only the wind is strengthened, but also accretion is been weakened at a certain distance.

In the figures $2.3.1.a$ to $2.3.2.b$ the viscosity is increased ($\alpha_{ss}=10^{-2}$), keeping the rotation same. For the MCG accretion case we see ${\cal X}_{fall}|_{\alpha_{ss}=10^{-2}}>{\cal X}_{fall}|_{\alpha_{ss}=10^{-4}}$. This means if viscosity is high, the weakening of accretion happens at further distance. This supports the conclusion that viscosity added with MCG decreases the accretion density as well.

4 Brief Discussions and Conclusions

In this letter, our motivation was to study the variations of density of accreting and wind fluids where DE is assumed to accrete onto SMBHs. Primarily, it is observed that for adiabatic flow if we are at a distance $x>x_{sonic}$ then $\rho_{wind}<\rho_{accretion}$. But once we are at the region $x<x_{sonic}$ then $\rho_{wind}>\rho_{accretion}$. Near to the BH, the wind density is very high. But if we consider DE to accrete in, the wind density increases remarkably. But the addition of viscosity gives us two distinct features. Firstly, the accretion disc’s radius is more shortened with the viscous flow. The second important thing is the sudden reduction of accreting density. As we increase the viscosity, we observe that this sudden falling point goes far from the BH. This signifies if DE is treated to be the necessary parameter, more viscous is the flow, more shortened is the length from BH at which the wind velocity becomes equal to that of light. Besides this strong outward wind created by the repulsive nature of DE, we observe that the accretion density is also reduced at a distance. A threshold drop of density in accretion profile can be pointed out clearly. So DE not only increases the power of wind it also reduces the power of accretion if viscosity and rotation are working together with it. To match the results found by us is compared with figure 0. We observe that our results are staying in the range of density predicted by the reference [5].

**Achkowledgement : ** This research is supported by the project grant of Government of West Bengal, Department of Higher Education, Science and Technology and Biotechnology (File no:- $ST/P/S\&T/16G-19/2017$). RB thanks IUCAA for providing Visiting Associateship. RB also thanks Prof. Banibrata Mukhopadhyay, Department of Physics, IISc, Bangalore, India-560012 for a detailed discussion regarding this problem earlier. SD thanks to Department of Higher Education, Government of West Bengal for research fellowship.