Circular orbits and accretion disk around AdS black holes surrounded by dark fluid with Chaplygin-like equation of state

In the present work we study the geodesic motion and accretion process of a test particle near an Anti-de Sitter (ADS) BH surrounded by a dark fluid with a Chaplygin-like equation. Within the defined paradigm, we investigate on the equatorial plane and examine circular geodesics along with their features related to stabilities, radiation energy flux, oscillations and orbits. The general form of the fluid accretion onto the AdS BH through dynamical analysis and mass expansion also has discussed in a depth. Additionally, a few more interesting topics, e.g. the effective potential, angular momentum, specific energy, radiation energy and epicyclic frequencies have also been examined thoroughly. All the attributes are physically acceptable within the observational signatures and ranges.


Introduction
As a theory of gravitation Einstein's General Relativity (GR) provides an in-depth and wide infrastructural facilities to understand all kinds of phenomena involved in the field of astrophysics and cosmology.Remarkably, Schwarzschild [1] first found out the solution to GR in connection to stellar model which can be treated as a toy model devoid of most of the real situations and physical aspects.Nevertheless, this historical solution provided an interesting feature of GR that this theory can conceive the idea of black hole (BH) in the spherically symmetric static metric.Here one may face an awkward situation for the condition either R = 2M or r → 0, where r, R and M , respectively are the radial coordinate, the radius of the stellar system and mass of the star.This odd and unexpected feature in the Schwarzschild [1] is the well-known singularity of a BH.At this juncture it may be a funny information that Einstein did not believe even in the existence of black holes where he refutes it by commenting that "Schwarzschild singularities" do not exist in physical reality [2,3].
However, the historical timeline how the "Schwarzschild singularities" gradually have taken shape over the time are as follows: (i) in 1939, Oppenheimer along with Volkoff worked on the collapsing star and possibility of the formation of a BH [4], (ii) again in the same year Oppenheimer and Snyder proposed a BH scenario in their paper 'On Continued Gravitational Contraction' [5], (iii) in 1958, Finkelstein pointed out the Schwarzschild surface as an event horizon [6], (iv) in 1960, an extension of the Schwarzschild model was proposed by Kruskal [7], (v) in 1963, Kerr obtained an exact solution for a rotating BH [8], (vi) in 1965, Newman [9] investigated on the axisymmetric solution for a BH with rotation and charge, (vii) in 1965, Penrose [10] found generic solutions to Einstein's GR which are involved in singularities, (viii) in 1967-1977, through a series of research papers the no-hair theorem emerged out which attributes that a stationary BH solution can completely be described by the three physical parameters, viz.mass, angular momentum and electric charge [11][12][13][14], (vii) (vii) in the period 1973-1975, a thorough investigations by Bardeen, Carter and Hawking [15] and Bekenstein [16][17][18] formulated the BH thermodynamics, (ix) in 1974, Hawking [19] raised an novel question whether "Black hole explosions?" takes place which eventually gave birth to the so-called Hawking radiation and (x) innumerable theoretical studies on BH by hundreds of scientists, recently studied a few of them are as follows [20][21][22][23][24][25][26][27][28].Most recently, some new black hole solutions have been calculated in the background of different aspect including gravitational decoupling approach and quantum induced [29][30][31][32][33][34][35][36][37][38][39][40][41][42][43] On the other hand, based on the observational evidence, in 1975, a galactic X-ray source named as Cygnus X-1, has been recognized as the first black hole [44,45].However, very recently on 11 February 2016, Abbott et al. [46] announced the first direct detection of gravitational waves which represents the observation of a BH merger.Again on 10 April 2019, the first direct image of a BH and its vicinity was published following observations made by the Event Horizon Telescope (EHT) in 2017 of the supermassive black hole in Messier 87's galactic centre [47,48].Till date there have been collective evidences on the real existence of BH [49].
In connection to major properties and structures, whether it is of primordial or stellar or galactic center-based BH, there are a lot of works available in the literature, viz.(i) the Hawking radiation which has already been mentioned [19] is a unique combination of GR and quantum mechanical effect, (ii) the Lense-Thirring effect which becomes indicative due to the angular momentum related frame dragging by the gravitomagnetic field [50], (iii) an event horizon which acts as an one-way boundary membrane in spacetime and allows matter as well as light-quanta to pass inside [51][52][53], (iv) the innermost stable circular orbit is a special feature of GR for which any infinitesimal inward perturbations to a circular orbit will lead to spiraling into the black hole [54] and (v) an accretion disk which is a structural formation around a BH by diffused material in order to maintain angular momentum [55].The impact of a cosmological constant on the test-particle motion in BH backgrounds was examined by Stuchlík [56].Stuchlík computed the circular orbits and the characteristics of purely radial trajectories around the charged black hole while accounting for the cosmological constant in the same study.In a different study, Schwarzschild-de Sitter and Schwarzschild-anti-de Sitter spacetime properties were examined by Stuchlík and his collaborator Hledík [57].They made use of the embedding diagrams, photon escape cones, and effective potential for the motion of test particles and photons.They demonstrate that in the Schwarzschild-de Sitter spacetimes, coincidence does not hold.In a recent review article, Stuchlík et al., [58] examined the effects of cosmic repulsion and external magnetic fields on jets connected to rotating BHs and accretion disks rotating around them against the background of the cosmological constant.
In the present work related to BH we are interested to focus on the two specific physical properties of it, viz.circular orbits and accretion disks.Regarding the circular orbits, the basic idea is that in GR there does exist an innermost stable circular orbit contrary to the Newtonian gravity where a test particle can orbit steadily in reference of a central massive object.Depending on the situations, i.e. energy of the test particle and spin of the BH, a circular orbit may lead to spiraling effect [54].As far as inward spiraling of a particle is concerned one can have a proximity to an accretion disc which is a resultant of accumulative particles due to tremendous gravitational attraction acting on the gaseous matter.
Therefore, based on the above mentioned literature and physical status of some of the astrophysical aspects, our attempt in the present article is to carry out investigation on some features of BH, especially circular orbits and accretion disks.
To reach at this goal we have presented first an AdS BH solution by considering the spherically symmetric spacetime in Section 2. In the next Section 3 we have provided geodesic motion regime under which circular motions (3.1), radiation energy flux (3.2) and oscillations (3.3) have been discussed.For the accretion disks related discussion we have chosen Section 4 where we have presented dynamic parameters (4.1) and mass expansion (4.2).We provide a brief overview of the circular equatorial geodesics 5 under which we specifically presented epicyclic frequencies (5.1).Finally, Section 6 is devoted to for discussion and concluding remarks.

AdS black hole surrounded by dark fluid with Chaplygin-like equation of state
To investigate the connection between the BH spacetime and the field when the BH is static spherically symmetric and its atmosphere is made up of a field with an explicit Lagrangian, it is practical to solve the gravitational field equations and the equation of motion of the relevant field simultaneously.When we are unable to understand the nature of the matter around the BH solutions due the influence of quintessence dark energy and Chaplygin gas, then it is crucial to think about how to investigate the relationship between matter and spacetime curvature using just the fluid matter's equation of state.For the current analysis, We consider the following form for the metric to describe our static spherically symmetric spacetime.
In order to matter source, we start with the perfect fluid, which is given as where p and ρ mention the isotropic pressure and energy density respectively.As pioneer in the working field, Semiz [59] has investigated the equation of state (EOS) p = ωρ (ω is a constant) for ideal fluid sources (dust, phantom energy, radiation, or dark energy) by using Einstein's equations as a basis.Further, Kiselev [60] calculated a novel BH spacetime by considering the ambient quintessence matter as an anisotropic fluid, which quickly turned into an extraordinarily popular toy model.For the Chaplygin dark fluid (CDF), from a field theoretical perspective, even now its background producing processes remain unclear [61,62], while there are a few plausible options in the framework of string theory [63,64].The CDF is typically represented by adding a self-interacting potential and a scalar field φ with the Lagrangian L φ = − 1 2 ∂ a φ∂ a φ − U (φ) [65,66].In this study, we assume that the CDF is anisotropic and that a covariant form of its stress-energy tensor may be expressed as [67] T where p r and p t are defining the radial and the tangential pressure components, u a defines the four-velocity of fluid and k a represents a unit space-like vector, which is orthogonal to u a , with u a and k a with satisfying u a u a = −1, k a k a = 1 and u a k a = 0. Π ab = g ab + u a u b − k a k b is a tensor projection onto two orthogonal surfaces u a and k a .When one operates within the fluid's comoving frame, they can achieve that u a = (− √ f , 0, 0, 0) and k a = (0, 1/ √ g, 0, 0).
Finally, the stress-energy tensor through Eq. (2.3) can be provided as The anisotropic factor is defined as the difference between radial and tangential pressures p r − p t .The stress-energy tensor simplifies to the conventional isotropic at p r = p t .Even if a cosmological fluid has anisotropy in the gravitational field generated by a BH, its equation of state should appear as p = p(ρ) at cosmological scale, allowing one to constrain the tangential pressure p t by taking an isotropic average over the angles and and some extra conditions ⟨T i j ⟩ = p(ρ)δ i j , one can get the following relation where the relation ) for quintessence matter, from the Eq. ( 2.5), one can use p t = 1 2 (1 + 3ω)ρ, which is compatible with the radial pressure p r = −ρ [60].The non-linear EOS for the CDF in our scenario is p = − B ρ , where B is a positive constant.The tangential pressure produces p t = 1 2 ρ − 3B 2ρ for p r = −ρ.As a result, the CDF's stress-energy tensor can be written as ) At cosmological scale, the anisotropy of the CDF vanishes and the EOS yields p = −B/ρ.Moreover, as an important condition, i.e., T t t = T r r is required, the relationship between the metric components g(r) = f (r) can be accomplished through a suitable time rescaling without sacrificing generality.The Einstein tensor components are calculated as ) (2.9) Now, by using the Eqs.(2.6)-(2.7)and Eqs.(2.8)-(2.9),one can get the following relation We now examine the existence of the cosmological constant.The two above differential equations allow us to calculate the two unknown functions, i.e., f (r) and ρ(r), analytically.Now, one can simply the energy density of CDF by solving the above set of differential equations (2.10)-(2.11)as: where a normalization factor q > 0 denotes the CDF's intensity.Additionally, the conservation law for the stress-energy tensor ∇ b T ab = 0 directly yields Eq. (2.12).It is evident that the CDF energy density can be approximated by r 6 ≪ q 2 /B as demonstrating that the CDF behaves as though its energy density were dependent on r −3 for matter content.Now, by using Eq.(2.12) into Eq.(2.10), we have the final analytical solution for f (r) [68] as where M is the BH mass, however in the current case, the BH mass is regarded as a point mass BH.Given r → ∞, f (r) can be adopted asymptotically The solution of the asymptotic behavior is examined by both the CDF parameter B and cosmological parameter Λ.In this work, we apprehension the AdS BH, so we constrain  To study the horizon structure, one can take f (r h ) = 0, and find the event horizon radius r h , with the mass of the BH can be expressed as The effect of parameters on the lapse f (r) is shown in Fig. (2).Some important insights regarding horizon structure are listed below: • In the upper plot by the left position has three solution curves black,red and green, in which black curve has no horizon for M = 0.45.Red curve has an event horizon for M = 0.53 and green curve has two horizon for M = 0.60 at the specific values of other parameters B = 1, q = 1 and Λ = −1.
• The upper right plot has three solution curves black,red and green, in which black curve has no horizon for B = 0.01.Red curve has an event horizon for B = 0.10 and green curve has two horizon for B = 0.50 at the specific values of other parameters M = 1, q = 1 and Λ = −1.• The bottom plot has three solution curves black,red and green, in which black curve has no horizon for q = 1.0.Red curve has an event horizon for q = 1.50 and green curve has two horizon for q = 2.0 at the specific values of other parameters M = 1, B = 1 and Λ = −1.
We have noted the interesting behavior of AdS BH, all the solution curves inclines towards the event horizon as all the parameters are increased in value of M , B and q.

Geodesic motion regime
In this Section, we compute the overall findings for particle geodesic motion by using the underlying symmetric and static space-time.For the analysis, we make a choice of two Killing vectors, say ε t = ∂ t and ε ϕ = ∂ ϕ under the effect of two constants of motion like E (conserved energy) and other one is L (angular momentum) with the following formulas: with u a = dx a dτ = (u t , u r , u θ , u ϕ ) are the four-velocity vectors for the moving particles, which fulfil the normalization condition, i.e., u a u a = −1, one can get the following relation: Veff max Along the equatorial plane, i.e., (θ = π 2 ), from Eqs. (3.1) and (3.2), we have the following relations From Eq. (3.3), the conserved energy equation for the particles motion can be obtained with an effective potential V ef f as The effective potential analysis of AdS BH plays a crucial role in geodesic motion, which clearly depends on parameter f (r) and the angular momentum L. The local extremum of the effective potential can be used to identify the position of the circular orbits.The following are the salient features of Fig. 3: • In the upper plot by the left position, the black curve has no extremum for L = 0.The first extremum is detected in red curve at V ef f = 0.99 for L = 2.2 √ 3.In this plot V ef f is increased for increasing L. In this plot, the values of other parameters are M = 1, β = 1, q = 1 and Λ = −1.
• In the left right plot, the black curve has an extremum for B = 1.0.The V ef f is increased for decreasing B. In this plot, the values of other parameters are M = 1, L = 5, q = 1 and Λ = −1.
• In the bottom plot, the black curve has an extremum for q = 1.0.The solution curves of V ef f are very small increased for increasing q.In this plot, the values of other parameters are M = 1, β = 1, L = 5 and Λ = −1.
We investigated the position of the orbits, which is represented by the point of extremum around the AdS BH at a distance of r ≈ 2, where the potential increases for increasing the angular momentum L and parameter q decreases for reducing the parameter B.

Circular motions
An equatorial plane can be used to investigate a particles circular motion.In this case, the radial component of circular motion r must be constant, i.e. u r = ur = 0. We get V e f f = E 2 and d dr V e f f = 0 from Eq. (3.4).These relationships allow us to write the particular energy E, specific angular momentum L, angular velocity Ω ϕ , and angular momentum L as: ) ) To achieve the specific energy and angular momentum, we have the following condition r 2f (r) − rf ′ (r) > 0. (3.10) One can acquire the limiting area of a circular orbit by solving the above condition.Additionally, two relations E 2 < 1 and E 2 = 1 are taken for the bound and marginally bound orbits respectively.Marginally bound orbits are easily discovered by using Eq.(3.6).Thus The photon sphere is characterized by this relation.The significant points of the specific energy, which is shown in Fig. 4 are as follows: • The specific energy is decreased gradually by increasing the parameter B. The curves are shifted inwards to the singularity.
• The specific energy is decreased gradually by decreasing the parameter q and vice versa.
The curves are shifted inwards to the singularity until it hits the singularity at r = 1.5 and go up to the maximum radius.
The angular momentum was described in Fig. 5, and its salient features are as follows: • The angular momentum is decreased gradually by increasing the parameter B. The curves are shifted onward to the singularity.However, the black and red curves are quickly fall onto the singularity and decreasing the radius.
• The angular momentum is decreased gradually by increasing the parameter q and vice versa.All the curves quickly fall onto the singularity and decreasing the radius and vice versa.

Radiation energy flux
The most explosive astrophysical events can be caused by radiation that is created when falling elements that are falling from rest at infinity accrete onto the black hole.This process releases gravitational energy from the falling elements.The radiation energy flux over the accretion disk is caused by the radiant energy associated with the specific energy E, angular momentum L and angular velocity Ω ϕ that were examined by Kato et al. [69] and can be provided as where Ṁ is an accretion rate, K is a radiation flux and g is another important parameter, which is defined as On the equatorial plane, one can obtain the following relation where W = (−rf ′ (r) + 2f (r)), X = (−f ′ (r) + rf ′′ (r)), Y = (rf ′ (r) + 2f (r)) and Now, by considering the relation K(r) = σT 4 (r), between the radiation and temperature.Here, σ is Stefan's constant.The detailed scenario of this Stefan's constant was discussed by [70] and further studied by [71].
• The black curve has a maximum radiation energy flux for B = 0.70 and q = 1.0 and it is decreasing for decreasing values of B and increasing q.The radius is decreased as the radiation flux increased by the parameters.
• The same behavior is seen in the picture of temperature and the radiation parameter Z as seen in the radiation energy flux.An additional essential evaluation of the accretion disk is the efficiency of the accreting particle.Consequently, to achieve the highest accreting efficiency, we use the relation η = 1 − E. We described the accreting efficiency in Fig. 8 and has the following key insights: • The efficiency is increased for increasing both the parameters B and q.
• Increasing efficiency increased the radius and the curves are go up to the maximum radius.
• We see the increasing efficiency and vice versa in left plot for taking specific value of q = 1 and variation in B.
• We see the increasing efficiency and vice versa in right plot for taking specific value of B = 1 and variation in q.

Oscillations
In general, three frequencies are needed to discuss the particle motion in an accretion disk.Harmonic motion within the vertical frequency, Ω θ , harmonic motion within the radial frequency, Ω r , and harmonic motion within the orbital frequency Ω ϕ .Here, the radial and vertical motions over the circular equatorial plane can be calculated.In this regard, we examine the following relations 1 2 dr dt ef f respectively.By using the Eq. ( 3.1), one can get radial motion u θ = 0 and vertical motion, u r = 0. Further, using u r = dr dτ = dr st u t and u θ = dθ dτ = dθ st u t , we have where . Then from Eq. (3.17), we obtain and

Accretion Regime
Herein, we shall formulate some general results of accretion.For these results, one can use the energy-momentum tensor of a perfect fluid with the following four-velocity vectors where τ represents the particles geodesic motion proper time.For the analysis of normalization criterion for spherically symmetric steady-state flow is u a u a = −1.We consider the four velocities of the present fluid as: The inward flow u r < 0 generates the negative velocity of the fluid and accretion.Whereas for the outward flow u r > 0, the fluid velocity should be positive.Therefore, for the present accretion analysis, we rely on the conservation law of energy and momentum tensor, which is further defined as: Now, by using the second kind of Christoffel symbol Γ and the covariant derivative ;, with √ −g = r 2 sin 2 θ, one can get the following relation where D 0 is the constant of integration.Now, by taking the befits from the relation of four-velocity and conservation law and through u a T ab ;b = 0, we have By the supposition u a u a = −1 and g ab ;b = 0, the above relation reduces to Taking only non-zero components, we obtain By integration where D 1 is the new constant of integration.Further, the minus sign is taken on the right-hand side that is u r < 0, we get where D 2 is the again constant of integration.The above setup produces the mass flux, which is further defined as: and it is also revised as Therefore, we have the following conservation law for mass as where D 3 is the constant like D 1 and D 2 .

Dynamical parameters
In this Subsection, we shall suppose equation of state p = kρ, where k is the state parameter for the isothermal fluid.In is necessary to mention that in an isothermal fluids, the temperature will be constant throughout the flow.The sound speed parameters should be remained constant for such kind of fluids p ∝ ρ throughout the accretion.Now, from Eqs. (4.9), (4.10) and (4.13), we have In the above equation D 4 represents the integration constant.Now, by taking p = kρ in above equation, one can obtain the general form of radial velocity as Therefore, the final form of radial velocity for current considered AdS BH is given as By considering Eq. (4.13), the general form of the density for the relevant fluid is defined as Now, the energy density for strong and weak fields is provided as We characterized the velocity of fluid in Fig. 9 and has the following key points: • The fluid velocity is increased for increasing B. The bound radius is decreased when the velocity increases.
• The fluid velocity is increased for increasing q.Similarly, the bound radius is decreased when the velocity increases.
We characterized the density behavior in Fig. 10 and has the following key points: • The fluid density is decreased for increasing B, until it goes up near to the singularity r = 2.0 at B = 1.0.
• The fluid density is increased for decreasing q, and all the solution curves quickly fall onto the singularity.We analyzed the mass accretion rate for AdS BH in Fig. 11 and has the following key points: • The maximum accretion rate of AdS BH occurs for (B = 0.7) and (q = 1.0).The radius is decreased gradually as increasing the accretion rate.
• Accretion rate is increased for increasing B and decreasing q at different distances.

Mass expansion
The mass of BH per unit of time can be measured by the mass accretion rate, which is equal to the area times flux at the boundary of BH.Further, it is usually denoted by Ṁ , and it can be determined by the metric parameters as well as the characteristics of the accreting fluid.It is defined as: where In order to investigate the mass expansion, we should supposed that our boundary at infinity, i.e., r = r ∞ .The peculiar infinity where the massive particles fall into the BH is indicated by the radius at infinity.Also, ρ ∞ is ρ energy density at r ∞ and p ∞ means p pressure at r ∞ .Hence, the above relation can be revised as: where On integrating, on can get where M i is the initial mass, M t is the mass of BH with the critical accretion time and is the critical accretion of time evolution.Therefore, the final expression for BH mass accretion rate is expressed as: (4.23)

Circular equatorial geodesics
The explicit form of the effective potential, specific energy, specific angular momentum of a moving particle is important for the physical analysis.From Eqs. (3.5), (3.6) and (3.7), we have and L 2 = 2M +

Epicyclic frequencies
The moving particles in a circular orbit achieve a small oscillations in the direction of vertical and radial frequencies.These oscillations are the special effects of perturbation on the moving particles in a circular orbit.Hence, the required frequencies are given by Br 6 r 2 +Br 3 + q arcsinh q √ Br 3 . (5.5) We characterized the radial and vertical frequencies and their ratios in Fig. 12 with the following key points: • The vertical frequency Ω θ increased by increasing B among very small change in the radius.
• The vertical frequency Ω θ increased by decreasing q and decreasing the radius.
• The radial frequency Ω r increased by decreasing B among very small change in the radius.
• The radial frequency Ω r increased by decreasing q and decreasing the radius.
• The ratio Ω θ Ωr transformed the solution curves from minimum radius to maximum radius.

Conclusion
The geodesic motion and accretion process of a test particle near an AdS BH surrounded by a dark fluid with a Chaplygin-like equation of state have investigated in the current analysis.Under this specified theme, we examined the equatorial plane and circular geodesics along with their stabilities, radiations energy flux, oscillations and orbits.The general form of the fluid accretion onto the AdS BH through dynamical analysis and mass expansion also has discussed.Some other interesting features like the effective potential, angular momentum, specific energy, radiation energy and epicyclic frequencies have been investigated.All the important highlights of the current analysis are listed as follows: • The effect of parameters on the horizon via metric function f (r) has provided in Fig. 2. We have observed that AdS BH has intriguing behavior; with larger values of all the parameters M , B and q, all of the solution curves are moved inward toward an event horizon.
• The effective potential for the considered AdS BH has presented in Fig. 3.As effective potential is very important for the geodesic motion to find the location of the circular orbits by the local extremum.The stability of the orbits is also depended on the lapse function and the angular momentum L. The point of extremum around the AdS BH at a distance of r ≈ 2 represents the orbits locations, which we have investigated.
Here, an increase in potential corresponds to an increase in angular momentum L and parameter q, whereas a reduction in parameter B results in a decrease in potential.
• The effect of involved parameters on specific energy E and angular momentum L have shown in Figs. 4 and 5 respectively.
• The graphical behavior of integrated factor Z from radiation parameter, radiation energy K(r) and temperature T with Stefan's constant of the considered AdS BH has given in Figs. 6 and 7 respectively.Further, accreting efficiency parameter has provided in Fig. 8.It has noted that the process releases gravitational energy in the falling elements.The radiating energy corresponding to the specific energy E, angular momentum L, and angular velocity is what causes the radiation energy flux over the accretion disk.
• The behavior of radial velocity, energy density, and mass accretion rate of the proposed, which has surrounded by dark fluid with Chaplygin-like equation of state have given in Figs. 9, 10 and 11 respectively.The radial velocity has increased for the increasing values of B. Interestingly, the bound radius has decreased when the velocity increased.
The fluid velocity has increased for increasing in q.Similarly, the bound radius has decreased when the velocity increased.It has observed that the maximum accretion rate of AdS BH occurs for (B = 0.7) and (q = 1.0).The radius has decreased gradually as increased the accretion rate.
• The effect of parameters B and q through different values on epicyclic frequencies of the AdS BH versus r has shown graphically in Fig. 12.
In conclusion, the current analysis of accretion disks around AdS black holes surrounded by dark fluid through a Chaplygin-like EOS has offered a significant and insightful observations regarding the behavior of these type of systems.It can be noticed that the stability, structure and energy transport mechanisms of the accretion disks are modified by the emergence of new dynamics and attributes brought about by the existence of a dark fluid with an EOS exhibiting the Chaplygin gas.

Fig. ( 1 )
shows how the metric function f (r) depends on the parameters B and q.

Figure 1 .
Figure 1.An illustration of the BH mass varying with q and B. We set Λ = -1.A physical BH mass needs to match the extremal BH mass at least, for the specific values of B, q, and Λ.

Figure 2 .
Figure 2.An illustration of horizons of the AdS BH (f (r)) versus r for different values of M , B and q.

Figure 3 .
Figure 3.An illustration of V ef f of the AdS BH versus r for different values of B and q.

Figure 4 .Figure 5 .
Figure 4.The illustration of the specific energy E of the AdS BH versus r for different values of B and q.

Figure 6 .
Figure 6.An illustration of Z of the AdS BH versus r for different values of B and q.

Figure 7 .
Figure 7.An illustration of the radiation energy K(r) and the temperature of the AdS BH versus r for different values of B and q.

Figure 8 .
Figure 8.An illustration of the efficiency of the AdS BH versus r for different values of B and q.

Figure 9 .Figure 10 .
Figure 9.An illustration of the radial velocity of the AdS BH versus r for different values of B and q.

Figure 11 .
Figure 11.An illustration of mass accretion rate of the AdS BH versus r for different values of B and q.

Figure 12 .
Figure 12.An illustration of epicyclic frequencies of the AdS BH versus r for different values of B and q.