Accretion onto Some Well-Known Regular Black Holes
In this work, we discuss the accretion onto static spherical symmetric regular black holes for specific choices of equation of state parameter. The underlying regular black holes are charged regular black hole using Fermi-Dirac Distribution, logistic distribution, nonlinear electrodynamics, respectively and Kehagias-Sftesos asymptotically flat regular black hole. We obtain the critical radius, critical speed and squared sound speed during the accretion process near the regular black holes. We also study the behavior of radial velocity, energy density and rate of change of mass for each regular black holes.
At present, type 1a supernova , Cosmic microwave background (CMB) radiation  and large scale structure [3, 4] have shown that our universe is currently in accelerating expansion period. Dark energy is responsible for this acceleration and it has strange property that violates the null energy condition (NEC) and weak energy condition (WEC) [5, 6] and produces strong repulsive gravitational effects. Recent observations suggests that approximately 74% of our universe is occupied by dark energy and the rest 22% and 4% is of dark matter and ordinary matter respectively. Nowadays dark energy is the most challenging problem in astrophysics. Many theories have been proposed to handle this important problem in last two decade. Dark energy is modeled using the relationship between energy density and pressure by perfect fluid with equation of state (EoS) . The candidates of dark energy are phantom like fluid , quintessence , cosmological constant . Other models are also proposed for explanation of dark energy like k-essence, DBI-essence, Hessence, dilation, tachyon, chaplygin gas etc -.
On the other hand, existence of essential singularities (which leads to various black holes (BHs)) is one of the major problems in general relativity (GR) and it seems to be a common property in most of the solutions of Einstein field equations. To avoid these singularities, regular BHs (RBHs) have been developed. These BHs are solutions of Einstein Equation with no essential singularity hence their metric is regular everywhere. Strong energy condition (SEC) is violated by these RBHs somewhere in space time [17, 18] while some of these satisfy the WEC. However, it is necessary for those RBHs to satisfy WEC having de sitter center. The study of RBHs solutions is very important to understand the gravitational collapse. Since Penrose cosmic censorship conjecture claims that is singularities predicted by GR [19, 20] occur, they must be explained by event horizons. Bardeen  has done a poineer work in this way by presenting the first RBH known as ”Bardeen Black Hole”, which satisfy the WEC.
The discussion about the properties of BHs have led many interesting phenomenon. Accretion onto the BHs is one of them. When a massive condensed object (e.g. black holes, neutron stars, stars etc.) try to capture a particle of the fluid from its surroundings, then the mass of condensed object has been effected. This process is known as accretion of fluid by condensed object. Due to accretion the planets and star form inhomogeneous regions of dust and gas. Supermassive BHs exist at center of giant galaxies which suggests that they could have formed through accretion process. It is not necessary that the mass of BH increases due to accretion process, sometimes in falling matter is thrown away like cosmic rays . First of all, the problem of accretion on compact object was investigated by Bondi using Newtonian theory of gravity . After that many researchers such as Michel , Babichev et al. [25, 26], Jamil  and Debnath  have discussed the accretion on Schwarzschild BH under different aspects. Kim et al  and Madrid et al.  studied accretion of dark energy on static BH and Kerr-Newman BH. Sharif and Abbas  discussed the accretion on stringy charged BHs due to phantom energy.
Recently, framework of accretion on general static spherical symmetric BHs has been presented by Bahamonde and Jamil . We have extended this general formalism for some RBHs. We analyze the effect of mass of RBH by choosing different values of EoS parameter. This paper is established as follows: In section 2, we derive general formalism for spherically static accretion process. In section 3, we discuss some RBHs and for each case, we explain the critical radius, critical points, speed of sound, radial velocities profile, energy density and rate of change of BH mass. In the end, we conclude our results.
2 General Formalism For Accretion
The generalized static spherical symmetry is characterized by the following line element
where , and are the functions of only. The energy-momentum tensor is considered as perfect fluid which is isotropic and inhomogeneous and defined as follows
where p is pressure, is energy density and is the four velocity which is given by
where is the proper time. We have and both equal to zero due to spherical symmetry restrictions. Here pressure, energy density and four velocity components are only the functions of . The normalization condition of four velocity must satisfy , we get
where , can be negative or positive due to square root which represents the backward or forward in time conditions. However, is required for accretion process otherwise for any outward flows . Both inward and outward flows are very important in astrophysics. One can assume that the fluid is dark energy or any kind of dark matter. For spherical symmetric BH, the proper dark energy model could be obtained by generalizing Michel’s theory. In dark energy accretion, Babichev et al.  have introduced the above generalization on Schwarzschild black hole. Similarly, some authors [22, 31] have extended this procedure for generalized static spherically symmetric BH. In these works, equation of continuity plays an important role which turns out to be
where is the constant of integration. Using , we obtain continuity (or relativistic energy flux) equation
Furthermore, assuming a certain EoS in this case. After some calculations, the above equation becomes
here prime represents the derivative with respect to . By integrating the last equation, we obtain
where is another constant depends upon and . Moreover, the equation of mass flux yields
In addition, we have introduced the variable
If the bracketed terms in Eq.(12) vanishes, we obtain the critical point (where speed of sound equal to speed of flow) which is located at . Hence at critical point, we get
and Eq.(12) turns out to be
Also, is the critical speed of flow evaluated at critical value . We can decoupled the above two equations and obtain
The speed of sound is evaluated at as follows
Obviously, and can never be negative and hence
Moreover, the rate of change of BH mass can be defined as follows 
Here dot is derivative with respect to time. We can observe that the mass of BH will increase for the fluid and hence the accretion occurs outside the BH. Otherwise, for like fluid, the mass of BH will decrease. The mass of BH cannot remain fixed because it will decrease in hawking radiation while it will increase in accretion. If we consider the time dependence of BH mass, then first assume that it will not change the geometry and symmetry of space time. Hence the space time metric remain static spherical symmetric .
3 Spherical Symmetric Metrics with Charged RBHs
In this section, we discuss the spherically symmetric metrics with charged RBHs in which . For this assumption, Equation (16) give
Although, our focus on charged RBHs metrics with event horizons, the present analysis is forbidden for horizon space time. In many cases, we concern on critical values (critical radius), critical velocities, speed of sound in fluid, behavior of energy density of fluid, radial velocity and rate of change of mass of accreting objects. So horizon does not involve anywhere .
3.1 Charged RBH Using Fermi-Dirac Distribution
The said RBH solution has the following metric functions 
where the Fermi-Dirac distribution function is
By replacing , we can obtain the distribution function as
with normalization factor is . Also the distribution function satisfies
where . Hence the metric functions turn out to be
If we set and , we obtain
In both equations, the difference of factor must be noted .
It is possible to integrate the conversation laws and obtain analytical expressions of the physical parameters. For simplicity, we will study the barotropic case where the fluid has an equation . Using (5) and (11), we obtain
The velocity profile for different values of is shown in Figure 1. Here refer to stiff, dust and cosmological constant respectively and and refer to quintessence and phantom energy. It can be seen that for the radial velocity of the fluid is negative and it is positive for . If the flow is outward then is not allowed and vice versa. In the case of the fluid is at rest at . Figure 2 represents the behavior of energy density of fluids in the surrounding area of BH. Obviously the WEC and DEC satisfied by dust, stiff and quintessence fluids. When phantom fluid () moves towards BH then energy density decreases and reverse will happen for dust, stiff and quintessence fluids (). Asymptotically at infinity for while it approaches to maximum at and near the BH.
Figure 3 represents the change in BH mass for different values of . The mass of the BH will increase near it and at for respectively. On the other hand, mass of BH decrease near it and at for . Hence the mass of BH increases due to accretion of quintessence, dust and stiff matter while it decreases due to accretion phantom like fluids.
The critical values, critical velocities and speed of sound are obtained for different values of the EoS parameter in Table 1. Critical radius is shifting to left when increases. Thus, the infalling the fluid acquires supersonic speeds closer to BH. Same critical radius is obtained for and with same critical velocities but opposite direction. We get negative speed of sound at and positive speed of sound for remaining critical radius. Also, the speed of sound increases near the RBH. For this metric, we find that
Also, the condition (18) yields
3.2 Charged RBH Using Logistic Distribution
The Logistic distribution function is 
in which replace then we obtain the distribution function
with normalization factor is . Also the distribution function satisfies
where . The horizons can be obtained for where . The metric function can be written as
If we set then we obtain the Schwarzschild BH and if we set we get
It is noteworthy that this metric function corresponds to an Ay´on-Beato and Garcýa BH .
The velocity profile for different values of is shown in Figure 4. It can be observed that for the radial velocity of the fluid is negative and it is positive for . If the flow is inward then is not allowed and vice versa. In the case of the fluid is at rest at . Figure 5, represents the behavior of energy density of fluids in the surrounding area of BH. Obviously the WEC and DEC are satisfied by dust, stiff and quintessence fluids. When phantom like fluid () moves towards BH then energy density decreases and reverse will happen for dust, stiff and quintessence fluids ().
The of RBH for distinct EoS parameters is obtained by using (19)
Figure 6 represents the change in BH mass against x. It is evident that the mass of BH increases due to quintessence, dust and stiff fluids and it decreases due to phantom fluids.
The critical radius, critical velocities and speed of sound are obtained for different values of EoS parameter in Table 2. Critical radius is shifting to right when increases. Thus the infalling the fluid acquires supersonic speeds closer to BH. For phantom like fluid, quintessence, dust and stiff matter the critical radius and critical velocities are explained in the above table. Same critical radius is obtained for and with same critical velocities but differ in sign. We obtained negative speed of sound at and positive speed of sound at . Near the BH speed of sound will increase. For this metric we find that
Also, the condition (18) yields
3.3 Charged RBH from Nonlinear Electrodynamics
Using the line element
Here the function
and its associated electric field source is
where q and M represent electric charge and mass respectively . The solution elaborate RBH and its global structure is like R-N BH. The asymptotic behavior of the solution is
So the metric function
The radial velocity and energy density for this metric are given by
The absolute value of the velocity profile for different values of is shown in Figure 7. It can be observed that for the radial velocity of the fluid is negative and it is positive for . If the flow is inward then is not allowed and vice versa. In the case of the fluid is at rest at . Figure 8 represents the energy density of fluids in the region of BH. It is apparent that the WEC and DEC is satisfied by phantom fluids. When phantom fluids moves towards the BH the energy density increases on the other hand it decreases for dust and stiff matter.
The rate of change of mass is given by
The rate of change of in the BH mass against x is plotted in Figure 9. Due to accretion of dust and stiff matter the mass of the BH will increase for small values of x and vice versa for phantom fluids. It is also noted that the maximum rate of RBH mass increases due to followed by .
The critical values, critical velocities and speed of sound are obtained for different values of EoS parameter in Table 3. Critical radius is shifting to right when increases. Speed of sound is negative at and near the BH the speed of sound will increase. For this RBH we find that
Also, the condition (18) yields
3.4 Kehagias - Sftesos asymptotically flat BH
KS studied the following BH metric
In the frame work of Horava theory, where m is the mass, b is the positive constant related to coupling constant of theory. The metric asymptotically behaves the usual Schwarzschild BH 
for . The KS metric have two horizons at
The radial velocity and energy density are given by
The radial velocity for different values of is shown in Figure 10. The radial velocity is negative for phantom like fluid and positive for quintessence, dust and stiff matter. The evolution of energy density of fluids in the surrounding area of RBH is plotted in Figure 11. The energy density for phantom fluids is negative while the energy density for stiff, dust and quintessence fluids is positive.
For this RBH, rate of change of mass becomes
Figure 11 represents the rate of change in RBH mass against x. We see that the RBH mass will increase for and it will decrease for .
The critical values, critical velocities and speed of sound for different values of is presented in the Table 4. For quintessence matter, we obtain very large critical radius. Similarly as before, we obtain the same critical radius for dust and phantom like fluids and same critical velocities but differ in sign. If we increase the EoS parameter then the critical radius is shifted near RBH. It is evident that the critical velocity is negative for phantom like fluid and positive for quintessence, dust and stiff matter. The speed of sound is negative at and positive for remaining critical radius. For this metric, we find that
The condition (18) becomes
4 Concluding Remarks
In this work, we have investigated the accretion onto various RBHs (such as RBH using Fermi-Dirac Distribution, RBH using logistic distribution, RBH using nonlinear electrodynamics and Kehagias-Sftesos asymptotically flat RBH) which are asymptotically leads to Schwarzschild and Reissner-Nordstrom BHs (most of them satisfy the WEC). We have followed the procedure of Bahamonde and Jamil  and obtained the critical points, critical velocities and the behavior of sound speed for chosen RBHs. Moreover, we have analyzed the behavior of radial velocity, energy density and rate of change of mass for RBHs for various EoS parameters. For calculating these quantities, we have assumed the barotropic EoS and found the relationship between the conservation law and barotropic EoS. We have found that the radial velocity () of the fluid is positive for stiff, dust and quintessence matter and it is negative for phantom-like fluids. If the flow is inward then is not allowed and is not allowed for outward flow. Also, we have obtained that the energy density remains positive for quintessence, dust and stiff matter while becomes negative for phantom-like fluid near RBHs.
In addition, the rate of mass of BH is dynamical quantity, so the analysis of the nature its mass in the presence of various dark energy models may become very interesting in the present scenario. Also, the sensitivity (increasing or decreasing) of BHs mass depend upon the nature of fluids which accretes onto it. Therefore, we have considered the various possibilities of accreting fluids such as dust and stiff matter, quintessence and phantom. We have found that the rate of change of mass of all RBHs increases for dust and stiff matter, quintessence-like fluid since these fluids do not have enough repulsive force. However, the mass of all RBHs decreases in the presence of phantom-like fluid (and the corresponding energy density and radial velocity becomes negative) because it has strong negative pressure. This result shows the consistency with several works [22, 31, 35, 36, 37, 38, 39]. Also, this result favor the phenomenon that universe undergoes the big rip singularity, where all the gravitationally bounded objects dispersed due to phantom dark energy.
Although, we have assumed the static fluid which may be extended for non-static fluid without assuming any EoS and can be obtained more interesting results. This is left for future considerations.
-  Perlmutter, et al., S.: Supernova Cosmology Project Collaboration. Astrophys. J. 517, 565 (1999).
-  Spergel, et al., D.N.: WMAP Collaboration. Astrophys. J. Suppl. 170, 377 (2007).
-  Eisenstein, et al., D.J.: SDSS Collaboration. Astrophys. J. 633, 560 (2005).
-  Riess, et al., A.G.: Supernova Search Team Collaboration. Astron. J. 116, 1009 (1998).
-  Johri, V.B.: Phys. Rev. D 70, 041303 (2004).
-  Lobo, F.S.N.: Phys. Rev. D 71, 084011 (2005).
-  Nojiri, S. and Odintsov, S.: Phys. Rep. 505, 59144 (2011).
-  Armendariz-Picon, C., Mukhanov, V.F. and Steinhardt, P.J.: Phys. Rev. Lett. 85, 4438 (2000).
-  Gasperini, et al., M.: Phys. Rev. D 65, 023508 (2002).
-  Gumjudpai, B. and Ward, J.: Phys. Rev. D 80, 023528 (2009).
-  Martin, J. and Yamaguchi, M.: Phys. Rev. D 77, 123508 (2008).
-  Wei, H., Cai, R.G. and Zeng, D.F.: Class. Quant. Grav. 22, 3189 (2005).
-  Sen, A.: JHEP 0207, 065 (2002).
-  Caldwell, R.R.: Phys. Lett. B 545, 23 (2002).
-  Kamenshchik, A.Y., Moschella, U. and Pasquier, V.: Phys. Lett. B 511, 265 (2001).
-  Copeland, E.J., Sami, M. and Tsujikawa,S.: Int. J. Mod. Phys. D 15, 1753 (2006).
-  Elizalde, E. and Hildebrandt, S. R.: Phys. Rev. D 65, 124024 (2002).
-  Zaslavskii, O. B.: Phys. Lett. B 688, 278 (2010).
-  Hawking, S.W. and Ellis, G.F.: The Large Scale Structure of SpaceTime (Cambridge Univ. Press 1973).
-  Senovilla, J.M.M.: Gen. Relat. Gravt. 30, 701 (1998).
-  Bardeen, J.: presented at GR5, Tiflis, U.S.S.R., and published in the conference proceedings in the U.S.S.R. (1968).
-  Bahamonde, S. and Jamil, M.: Eur. Phys. J. C 75, 508 (2015).
-  Bondi, H.: Mon. Not. R. Astron. Soc. 112, 195 (1952).
-  Michel, F.C.: Astrophys. Space Sci. 15, 153 (1972).
-  Babichev, et al., E.: Phys. Rev. Lett. 93, 021102 (2004).
-  Babichev, E., Dokuchaev, V. and Eroshenko, Y.: J. Exp. Theor. Phys. 100, 528538 (2005).
-  Jamil, M.: Eur. Phys. J. C 62, 609 (2009).
-  Sharif, M. and Abbas, G.: Chin. Phys. Lett. 29, 010401 (2012).
-  Kim, S.W. and Kang, Y.: Int. J. Mod. Phys. Conf. Ser. 12, 320 (2012).
-  Jimenez, Madrid, J.A. and Gonzalez-Diaz, P.F.: Grav. Cosmol. 14, 213 (2008).
-  Debnath, U.: Eur. Phys. J. C 75, 129 (2015).
-  Leonardo, Balart and Elias, C. Vagenas: Phys. Rev. D 90, 124045 (2014).
-  Sharif, M. and Jawad, A.: Mod. Phys. Lett. A 25, 3241-3250(2010).
-  Culetu, H.: Astrophys. Space Sci. 2, 360(2015).
-  Debnath, U.: Eur. Phys. J. C 75, 449 (2015).
-  Wei, H.: Class. Quantum Grav. 29(2012)175008.
-  Lobo, F.S.N.: Phys. Rev. D 71(2005)124022; Lobo, F.S.N.: Phys. Rev. D 71(2005)084011; Sushkov, S.: Phys. Rev. D 71(2005)043520.
-  Babichev, E., Dokuchaev, V. and Eroshenko, Y.: Phys. Rev. Lett. 93(2004)021102.
-  Sharif, M. and Abbas, G.: Chin. Phys. Lett. 28(2011)090402; Martin-Moruno, P.: Phys. Lett. B 659(2008)40; Jamil, M., Rashid, M.A. and Qadir, A.: Eur. Phys. J. C 58(2008)325; Babichev, E. et al.: Phys. Rev. D 78(2008)104027; Jamil, M.: Eur. Phys. J. C 62(2009)325; Jamil, M. and Qadir, A.: Gen. Rel. Grav. 43(2011)1069; Bhadra, J. and Debnath, U.: Eur. Phys. J. C 72(2012)1912.