Abstract

Our universe is found to expand sometimes with a rapidly accelerated speed, sometimes being retarded by gravity and somewhere with a steady accelerating nature. American astronomer Henrietta Swan Leavitt, in 1912, established a relation between the regular luminosity period of a Cepheid star with its absolute magnitude. Using this as a tool, distances of nearby galaxies hosting Cepheid variables were measured. One contemporary scientist of her time, Vesto Slipher noticed redshifts in far away galaxies.

Meanwhile, two publications from Russia(1922) and Belgium(1927) came in front (by Alexander Friedmann and Georges Lemaître respectively) which have shown different upcoming possibilities for the cosmos, especially a tendency to grow a nonstatic future which was contrary to the idea of Einstein’s theory.

In 1929, redshifts of several spiral galaxies were measured by Milton Humason. Humason, with his collaborator Edwin Hubble, soon found the distances of Cepheid variables located in these galaxies. Standard luminosity of Cepheid variables allows to measure both the distance and redshift of host galaxies which reside in “nearby” region. A linear fit to such a redshift vs. distance data was plotted and the corresponding analysis became popular as Hubble’s law which supported the expansion of cosmos. However, expectations that the speed of the universe will slow down with time were formed.

Another kind of standard candles, type 1a supernovae residing in different farthest galaxies were observed by two teams, captained by astronomers Adam Riess, Saul Perlmutter and Brian Schmidt around 1998 which predicted the late time cosmic acceleration [1, 2, 3, 4, 5]. Justifications towards the question why our universe is expanding faster over time were tried to build via two major pathways. Modification of gravity and inclusion of exotic stress energy.

The first mentioned way led to different modified gravities. General Relativity (GR) was able to perfectly support solar system gravity precision tests , e.g., gravitational redshift, gravitational lensing of distant stars’ light, Mercury’s perihelion precision, Shapiro’s time delay effect etc. But in the outer world, some cataclysmic events like black hole-black hole merger [6] have pointed that a quantum correction of GR to justify such strong gravitational effects is required.

A \( d\) -dimensional manifold \( \cal{M}\) invariant under diffeomorphism, equipped with a metric \( g_{\mu \nu}\) , constructs the base of GR as a field theory. The action of this case looks like

Here \( \tilde{R}\) is scalar curvature of the concerned manifold, \( \tilde{\Lambda}\) is the cosmological constant, \( g=det(g_{\mu \nu})\) and \( \cal{L_M}\) denotes a generic matter Lagrangian. Generalised Newton’s constant \( G_N^d\) has the energy dimension \( [G_N^{(d)}]=2-d\) which reduces to \( [G_N^{(4)}]=-2\) in \( 4D\) .

Variating equation (1) with respect to the metric \( g_{\mu \nu}\) , we obtain the Einstein field equation (\( c=1\) )

here the energy momentum tensor \( T_{\mu\nu} = -\frac{2}{\sqrt{|g|}} \frac{\delta    }{\delta g^{\mu \nu}} \left(\sqrt{|g|} \cal{L_M}\right)\) . Modification in the Lagrangian led us to different modified gravities.

In the next section (i.e., in \( Section 2\) ), we will briefly state about a particular type of black hole solution in Horava Lifshitz(HL) gravity, popularly known as Kehagias Sfetsos solution. In \( Section 3\) , some special kind of dark energy models will be discussed. In \( Section 4\) dark energy accretion model will be built for KS black hole. In the last section, we will discuss the results in brief.

The concept of Horava-Lifshitz \( (HL)\) gravity incorporates the manifold \( \cal{M}\) with a foliation structure, which lets the metric to be written in the form [7]

where \( N\) , \( N^i\) and \( \gamma_{ij}\) are respectively the lapse, shift and spatial leaf metric.

The said foliation preserves diffeomorphism \( t\rightarrow t'(t)\) , \( x^i \rightarrow x^{i' t}(t,x)\) . For this transformation, \( N\) , \( N^i\) and \( \gamma_{i j}\) change as

Extrinsic curvature \( K_{i j}\) covariantly transforms under these changes as

With this, the modified action for \( HL\) gravity turns to be

where \( \sqrt{|g|}=N\sqrt{\gamma}\) , \( D=d-1\) is the dimension of spatial leafs, \( \lambda\) is a dimensionless coupling and \( \cal{V}\) contains only spatial derivatives. \( K_{ij}\) is the only covariant tensor which involves time derivatives. This makes the action second order in time and let to avoid the presence of Orstrogradsky ghosts for all observers. If \( {\cal V} = R - 2\Lambda\) and \( \lambda=1\) , the action noted reproduces the \( D+1\) decomposition of the Einstein-Hilbert action with the symmetry group enhanced to full diffeomorphisms.

For a detailed balance conditioned \( HL\) gravity, the action takes the explicit form

where \( w\) , \( \kappa\) and \( \lambda\) symbolize coupling constants which are dimensionless. \( \Lambda\) is the positive cosmological constant related to the Infrared limit and posessing mass dimension \( 2\) . An antisymmetric tensor \( \epsilon^{lmt}\) and a parameter \( \mu\) with mass dimension \( 1\) are introduced. The extrinsic curvature \( K_{lm}\) is now written with respect to the metric tensor

Cotton tensor \( C^{lm}= \epsilon^{ltq} \nabla_t (R^m_q - \frac{1}{4} \delta^m_q).\) Using the projectibility condition of \( HL\) gravity [8] and choosing the Friedmann-Robertson-Walker metric as \( N=1\) , \( g_{lm}=a^2(t) \gamma_{lm}\) , \( N^l=0\) and

where the curvature parameter \( K\) can opt the values \( -1,  0\) \( \& \) \( 1\) when it specifies open, flat or closed universe respectively. \( a(t)\) is the scale factor.

Choosing dark matter to take the form of a perfect fluid, the related Friedmann equation transforms into

Here \( H=\frac{\dot{a}}{a}\) is the Hubble parameter and \( \rho_{tot}\) is usual energy density \( =\rho_m + \rho_{\phi}\) . A Schwarzschild like solution was obtained in HL gravity by Kehagias and Sfetsos [9]. Considering \( \lambda =1\) the fluid equation turns

when a flat space is chosen, i.e., the curvature parameter \( k=0\) , higher order derivative terms will give no contribution in action. However, for a nonzero \( k\) significant changes are expected with small value of the scale factor \( a(t)\) .

The KS solution is given by [10]

where \( A(r)\) is given as

where \( w = \frac{16 \mu^2}{{\cal{\kappa}}^2}.\)

In section \( Section 3\) , we will now recapitulate regarding some dark energy candidates for which equation of states are functions of redshift.

We will shift our eye towards the modification of stress energy side of field equation. This is commonly known as the involvement of dark energy (DE). Several candidates are proposed. ‘Quintessence’ is the first kind of candidate where equation of state looks like [11]

where \( \phi\) and \( V\) are respectively related scalar field and potential.

‘Phantom’ energy is another kind of model [12] which generates quick acceleration leading to future Cosmic singularity named “Big Rip”.

Combining quintessence and phantom, another hypothetical scenario regarding quintom DE model is proposed. For low redshift, EoS looks like,

which, for high redshift, turns

Free parameters of this model are determined as \( w_0 =1.06 \pm 0.14\) and \( w_1=0.36 \pm 0.62\) [13]. This model can not account for why cosmic single fluid model or a scalar field model fails to realise available Quintom model in FLRW universe. To justify, either a non-minimal coupling is to be introduced or extra degrees of freedom will be required [14].

Next generation DE models are proposed as

Amount of acceleration will be high if \( \left. \left(\frac{d w}{dz} \right) \right|_{z=0}\) is negative. Assuming a flat universe and applying data from SDSS, dimensionless density for dark matter, \( \Omega_{DM} = 0.3 \pm 0.1\) , EoS is found to be \( w(z) \sim -1\) . Linear Redshift Parameterization (LRP) [15] has the equation of state

In the article [16], Chevallier-Polarski-Linder (CPL) Redshift parameterization model is proposed which is able to reproduce the cosmic picture for \( z \leq 2\) . At \( z=0\) , \( \Lambda\) CDM is almost perfectly mimicked. A special care has been taken where Cosmic Microwave Background (CMB) peak is formed. EoS is given as

Jassal-Bagala-Padmanabhan (JBP) Redshift Parameterization [17] proposed a model with the EoS

This model is shown to restrict huge variation of energy density at low redshift when SNeIa observations with the constraints from WMAP observation are applied.

In the article [18], Barboza-Alcaniz (BA) redshift parameterization EoS is proposed as,

while constraining free parameters, we find for \( 68.3\%\) confidence, \( w_0= -1.14_{-0.24}^{+0.31}\) , \( w_1 = 0.84_{-1.59}^{0.65}\) and \( \Omega_m^0=0.27 \pm 0.003\) . Cosmos depicted in this model is found to transit from deceleration to acceleration at \( z_t \simeq 0.58\) .

Barboza-Alcaniz-Zhu-Silva (BAZS) [19] have considered another parameterization

i.e.,

We have constrained three free parameters, we have found three kinds of confidence contour graphs and the best fit values of \( w_0^{BAZS}\) , \( w_{\beta}\) and \( \beta\) .

We will solve different fluid’s densities using the cosmological equation of continuity

Energy density of all the said DE models can be determined. In the next section, we will study the accretion of different DE candidates upon Kehagias Sfetsos black holes.

Assuming the fluid to follow perfect fluid nature and recalling the energy momentum tensor of DE

\( \rho\) and \( p\) respectively are the energy density and the pressure of DE. For the present work, we will choose \( \rho=\rho_{\phi}^{i}\) , \( p=p^i_{\phi}\) , where ‘\( i\) ’ will tell exactly which DE model is chosen. \( u^\mu\equiv \frac{dx^\mu}{d\tau}\) is the four speed for the \( \mu^{th}\) coordinate, \( x^\mu\) . \( \tau\) is proper time. 4-velocity follows the norm property, \( u_\mu u^\nu=1\) . We also suppose that the black hole’s spherical symmetry is not influenced by the accreted DE.

Without any loss of generality, we assume \( u^1=u<0\) (as this is radially infalling) and \( u^2=u^3=0\) due to symmetry and velocity components are derived as

To construct the continuity equation, we consider the components of \( T^{\nu}_{\mu}\) as

Along with the 4-current, conservation of mass flux \( \frac{\partial J^\mu}{\partial x^\mu} = 0\) can be written as

and opting the \( 0^{th}\) component of \( T^{\nu}_{\mu ; \nu} = 0\) , we obtain

where \( \lambda_1\) is another arbitrary constant of integration.

Next we will consider the energy momentum conservation \( u^\mu T^\nu_{\mu ; \nu} =0\) , i.e., the projection of energy momentum conservation law on one fluid’s 4-velocity and we obtain [20]

\( \lambda_2 > 0\) symbolizes another integrating constant related the energy flux [20, 21, 22, 23]. \( \rho\) and \( \rho_{\infty}\) represent the values of the density at finite distance and at infinity respectively.

We choose a dimensionless radial distance parameter \( x=\frac{r}{M}\) , \( M\) is the mass of central gravitating object. Hence equation (23) reduces to

Now we choose a formal arbitrary function \( n\) such that

Here the medium is chosen to possess only individual conserved particles [22] and the corresponding number density is taken as \( n\) .

Integrating equation (24) we obtain

where \( n_\infty\) is called the concentration of DE [22].

From equation (22), we obtain

This converts equation (27) into

To determine \( \lambda_2\) , which has the dimension of energy flux, we evaluate different fluid parameters at the critical point. Let us consider a parameter having the same dimension as the velocity given by

We can conclude that \( V^2\) is equivalent to the effective sound speed squared of the concerned accreting medium. We symbolize and define it as

Now differentiating equations (24) and (28) and using (26), we obtain

Clearly the radial inward speed gradient \( \frac{du}{dx}\) is now expressed as a numerator to denominator ratio. It is followed, in the interval \( 0\leq u \leq \) speed of light, the denominator turns zero at some critical distance \( x=x_c\) .

for a physical flow, the numerator should vanish at the same critical point \( x_c\) as well. Hence at \( x=x_c\) , we obtain two relations

and

where \( '\equiv \frac{d}{dx}\) denotes the differentiations with respect to \( x\) .

Now, using equations (32) and (33), from equation (27), we derive,

Wherever the critical point \( x_c\) takes place, we can calculate the density \( \rho_c = \rho(x_c)\) there. Again if a particular value of \( \rho_c\) is given, we can obtain corresponding \( x_c\) and \( u_c\) .

Hence it will be very straight forward to calculate

So, for different redshift parametrizations, we can obtain the constants \( \lambda_2.\)

The rate of change of mass of the black hole is calculated as

Different references [10, 24, 25, 26, 27] suggests that the rate of change of mass can be calculated at every pair of \( \left(\rho,  p(\rho)\right)\) , which does not satisfy the dominant energy condition, whenever we apply the equation of state \( p=w\rho\) .

Hence

Sign of \( \dot{M}\) is solely determined by that of \( (\rho+p)\) . For quintessence \( -1 < \frac{p}{\rho} < -\frac{1}{3}\) and hence \( \dot{M}>0\) . But once phantom barrier \( (\frac{p}{\rho}= -1)\) is crossed, \( p+\rho< 0\) and hence \( \dot{M}<0\) .

here \( M_0 (M(z=0))\) is the present day mass of the black hole. Present time energy density \( \rho(z=0)=\rho_0\) is constituted of two parts : present time matter density \( \rho_m\) and same for dark energy, i.e., \( \rho_{DE_0}\) .

Hence we obtain the equation for 4-D Kehagias-Sfetsos(KS) black hole in Horava-Lifshitz (HL) gravity as

We consider two main contributories to the density and pressure, viz., DE density \( (\rho_{DE})\) and matter density \( (\rho_m)\) .

Hence,

and

Matter density is solved from (18) as [28]

where \( w_m = \frac{p_m}{\rho_m}\) is the EoS for matter contribution and \( \rho_{m_0} (=\rho_m(z=0))\) is the present time value of matter density. In terms of dimensionless density parameter, \( \Omega_{m_0} = \frac{\rho_{m_0}}{3 H_0^2}\) , \( \rho_m\) can be written as

With these details, we are going to study the accretion properties of four different DE models in the next four subsections.

4.1 Accretion of Dark Energy Following Chevallier-Polarski-Linder Parametrization onto Kehagias-Sfetsos Black Hole

Solving the equation of continuity (18) for the CPL EoS (14), we obtain the energy density of CPL type DE as

\[ \begin{equation} \rho^{CPL}_{\phi} (z) = \rho^{CPL}_{\phi 0}~ exp\bigg\{-3w_1^{CPL}~\frac{z}{1+z}\bigg\}(1+z)^{3\left(1+w^{CPL}_0+w^{CPL}_1\right)}~~~~~~~~, \end{equation} \]

(42)

where \( \rho_{\phi 0}^{CPL}\) is the DE energy density at \( z=0\) . Hence the total energy density of universe filled up by CPL type of fluid is evolved as

\[ \begin{equation} \rho^{CPL}_{tot}(z) = 3H_0^2 \left[\Omega_m (1+z)^3 + \Omega^{CPL}_{\phi 0}~ exp\bigg\{-3w_1^{CPL}~\frac{z}{1+z}\bigg\}(1+z)^{3\left(1+w^{CPL}_0+w^{CPL}_1\right)}\right]~~~~. \end{equation} \]

(43)

So, utilizing the equation (43) in the equation (39), where \( \Omega_{\phi 0}^{CPL}= \frac{\rho_{\phi 0}^{CPL}}{3H_0^2}\) is the dimensionless fractional density for CPL, we get for CPL parametrization,

\[ \begin{equation} M^{CPL}=\frac{M_0^{CPL}}{1+\frac{4 \pi \lambda_2^{CPL} M_0}{3} \int _{\rho_{\phi_0}^{CPL}}^{\rho_{\phi}^{CPL}}\frac{d\left(3H_0^2 \left[\Omega_m (1+z)^3 + \Omega^{CPL}_{\phi 0}~ exp\bigg\{-3w_1^{CPL}~\frac{z}{1+z}\bigg\}(1+z)^{3(1+w^{CPL}_0+w^{CPL}_1)}\right]\right)}{\sqrt{\frac{\kappa^2}{12}\left[\left(3H_0^2 \left[\Omega_m (1+z)^3 + \Omega^{CPL}_{\phi 0}~ exp\bigg\{-3w_1^{CPL}~\frac{z}{1+z}\bigg\}(1+z)^{3(1+w^{CPL}_0+w^{CPL}_1)}\right]\right)-6k\mu^4 (1+z)^2-\frac{3k\kappa^2 \mu^2 (1+z)^4}{16}\right]}}}~~~~. \end{equation} \]

(44)

The constant \( \lambda_2^{CPL}\) is calculated using the equation (14) to derive the expression

\[ \frac{d\rho'}{\rho'+p(\rho')}=\frac{d\rho'}{\left\{1+w_0+w_1-w_1a(\rho')\right\}\rho'}~~~~. \]

We solve the scale factor \( a(\rho')\) numerically as the function of density over a large range. For a Bondi type accretion, we have followed that \( x_c\) is likely to occur somewhere inbetween 4.5 to 10 Schwarzschild radius. [29], [30] Performing numerical integration from a large distant point to such a limited critical value, we have evaluated \( \lambda_2^{CPL}\) . We have plotted growing mass to zero redshift mass ratio in figure 1a to 1c for different sets of best fit values of the free parameters.

In figure 1a, where the DE free parameters are chosen as \( w_0=-1\) and \( w_1=0.81\) [31], matter equation of state is varied from stiff fluid (\( \frac{p}{\rho}=1\) ) to quintessence barrier (\( \frac{p}{\rho}=-\frac{1}{3}\) ). As \( z\) decreases, mass is followed to increase first slowly and then rapidly. Mass at a particular redshift is higher if \( w_m\) is low. For a better insight, we have taken logarithm. Shifting to figure 1b and 1c, we observe the same pattern is kept. But near the \( z=0\) point, this ratio is higher as the values of \( w_1\) is reduced.

4.2 Accretion onto Kehagias-Sfetsos Black hole with Dark Energy Following Jassal-Bagala-Padmanabhan (JBP) Redshift Parameterization

For this candidate, the fractional density contribution turns (using equation (15) in the equation (18)).

\[ \begin{equation} \rho^{JBP}_{\phi}(z) = \rho^{JBP}_{\phi 0}~ exp\bigg\{\frac{3w^{JBP}_1}{2(1+z)}\bigg\}(1+z)^{3(1+w^{JBP}_0)}~~~~~~~~, \end{equation} \]

(45)

where \( \rho^{JBP}_{\phi 0} = \rho^{JBP}_{\phi}\) (at \( z=0\) ), comes as the constant of integration.

\[ \begin{equation} \rho^{JBP}_{tot}(z) = 3H_0^2 \left[\Omega_{m0} (1+z)^{3(1+w_m)} + \Omega^{JBP}_{\phi 0} exp\bigg\{ \frac{3w^{JBP}_1}{2(1+z)}\bigg\}(1+z)^{3(1+w^{JBP}_0)} \right]~~~~~~~, \end{equation} \]

(46)

in terms of dimensionless density parameter \( \Omega^{JBP}_{\phi 0}=\frac{\rho^{JBP}_{\phi 0}}{3 H_0^2}\) . So, puttting equation (46) in the equation (39), we will get for JBP parameterization,

\[ \begin{equation} M^{JBP}=\frac{M_0^{JBP}}{1+\frac{4 \pi \lambda_2^{JBP} M_0}{3} \int _{\rho_{\phi_0}^{JBP}}^{\rho_{\phi}^{JBP}}\frac{d\left(3H_0^2 \left[\Omega_{m0} (1+z)^3 + \Omega^{JBP}_{\phi 0} exp\bigg\{ \frac{3w^{JBP}_1}{2(1+z)}\bigg\}(1+z)^{3(1+w^{JBP}_0)} \right] \right)}{\sqrt{\frac{\kappa^2}{12}\left[\left(3H_0^2 \left[\Omega_{m0} (1+z)^3 + \Omega^{JBP}_{\phi 0} exp\bigg\{ \frac{3w^{JBP}_1}{2(1+z)}\bigg\}(1+z)^{3(1+w^{JBP}_0)} \right] \right)-6k\mu^4 (1+z)^2-\frac{3k\kappa^2 \mu^2 (1+z)^4}{16}\right]}}}~, \end{equation} \]

(47)

\( M^{JBP}\) is the BH mass in JBP surroundings at an arbitrary redshift \( z\) (\( M^{JBP}\) ) and \( M_0 (M(z=0))^{JBP}\) is the present mass of the black hole.

We can calculate the constant \( \lambda_2^{JBP}\) using

\[ \frac{d\rho'}{\rho'+p(\rho')}=\frac{d\rho'}{\left\{1+w_0+w_1a(\rho')\left(1-a(\rho')\right)\right\}\rho'}~~~~. \]

Though the inverse function \( a(\rho)\) can not be explicitly determined. We have calculated \( \lambda^{JBP}_2\) by the help of numerical integration. \( log_{10}\left(\frac{M^{JBP}}{M_0^{JBP}}\right)\) is plotted in figure 2a-c.

For JBP, near \( z=0\) , accretion rates are higher than the CPL case. Here if \( w_1\) is less, mass at the same redshift is less. Hence accretion is powerful if \( w_1\) is high. As CPL, in this case also, when matter part is following lesser value of equation of state accretion is powerful. Even we observe the accretion for \( w_m=-\frac{1}{3}\) and \( w_m=1\) differ by \( \sim 10^2\) of order around \( z\sim 3\) .

4.3 Accretion upon Kehagias-Sfetsos Black hole as Accreting Fluid Barboza-Alcaniz (BA) Redshift Parameterization

Incorporating the equation of state (16) in the cosmological equation of continuity, the fractional distribution of BA parameterization model towards the density turns

\[ \begin{equation} \rho^{BA}_{\phi}(z) = \rho^{BA}_{\phi 0} \bigg\{\frac{1+z^2}{(1+z)^2}\bigg\}^{\frac{3w^{BA}_1}{2}} \times (1+z)^{3(w^{BA}_0 + w^{BA}_1)}~~~~~~~~. \end{equation} \]

(48)

where \( \rho_{\phi 0}^{BA}\) is obtained at \( z=0\) . Hence the density of the concerned model turns \( \Omega^{BA}_{\phi 0}=\frac{\rho^{BA}_{\phi 0}}{3 H_0^2}\) .

\[ \begin{equation} \rho^{BA}_{tot}(z) = 3H_0^2 \left[\Omega_{m0} (1+z)^{3(1+w_m)} + \Omega^{BA}_{\phi 0} \bigg\{\frac{1+z^2}{(1+z)^2}\bigg\}^{\frac{3w^{BA}_1}{2}} \times (1+z)^{3(w^{BA}_0 + w^{BA}_1)}\right]~~~~~~~~, \end{equation} \]

(49)

So, incorporating equation (49) in the equation (39), we will get for BA parameterization,

\[ \begin{equation} M^{BA}=\frac{M_0^{BA}}{1+\frac{4 \pi \lambda^{BA}_1 M_0}{3} \int _{\rho_{\phi_0}^{BA}}^{\rho_{\phi}^{BA}}\frac{d\left(3H_0^2 \left[\Omega_{m0} (1+z)^3 + \Omega^{BA}_{\phi 0} \bigg\{\frac{1+z^2}{(1+z)^2}\bigg\}^{\frac{3w^{BA}_1}{2}} \times (1+z)^{3(w^{BA}_0 + w^{BA}_1)} \right] \right)}{\sqrt{\frac{\kappa^2}{12}\left[\left(3H_0^2 \left[\Omega_{m0} (1+z)^3 + \Omega^{BA}_{\phi 0} \bigg\{\frac{1+z^2}{(1+z)^2}\bigg\}^{\frac{3w^{BA}_1}{2}} \times (1+z)^{3(w^{BA}_0 + w^{BA}_1)} \right] \right)-6k\mu^4 (1+z)^2-\frac{3k\kappa^2 \mu^2 (1+z)^4}{16}\right]}}}~~, \end{equation} \]

(50)

\( M^{BA}\) is the black hole mass embedded in BA at an arbitrary redshift (\( M^{BA}\) ) and \( M_0^{BA} (M^{BA}(z=0))\) is the present mass of the black hole in BA.

\( \lambda_1^{BA}\) is determined by numerical integration and another numerical integration from a large distance to \( x_c\) helped us to plot \( log_{10}\left(\frac{M}{M_0}\right)\) in fig 3a-c. BAZS accretion curves are found different than CPL and JBP cases. Clearly, two phases of accretion are found. Upto a transition value of redshift, \( z_t\) , mass increases slowly. As we enter in the area \( z<z_t\) , a rapid increase is followed. This rapidness is again larger once the matter contribution possesses matter EoS \( w_m<0\) .

4.4 Accretion of Kehagias-Sfetsos Black hole with Dark Energy Following Barboza-Alcaniz-Zhu-Silva Redshift Parameterization (BAZS) Parameterization

Solving the equation of continuity (18) for the BAZS EoS (17), we find the energy density of BAZS type DE as,

\[ \begin{equation} \rho^{BAZS}_{\phi}(z) = \rho^{BAZS}_{\phi 0}(1+z)^{\frac{3\{w_{\beta}+\beta(1+w^{BAZS}_0)\}}{\beta}}exp\left[\frac{3w_{\beta}}{\beta^2}\{(1+z)^{-\beta} -1\}\right]~~~~, \end{equation} \]

(51)

where \( \rho_{\phi 0}^{BAZS}\) is the DE density at present epoch. Hence the total energy density of universe filled up by BAZS type candidate is calculated as,

\[ \begin{equation} \rho^{BAZS}_{tot}(z) = 3H_0^2 \left[\Omega_m (1+z)^3 + \Omega^{BAZS}_{\phi 0}(1+z)^{\frac{3\{w_{\beta}+\beta(1+w^{BAZS}_0)\}}{\beta}}exp\left[\frac{3w_{\beta}}{\beta^2}\{(1+z)^{-\beta} -1\}\right]\right]~~~~. \end{equation} \]

(52)

So, putting the equation (52) in the equation (39) we get for BAZS parameterization,

\[ M^{BAZS}=~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ \]

\[ \begin{equation} \frac{M_0^{BAZS}}{1+\frac{4 \pi \lambda^{BAZS}_1 M_0}{3} \int _{\rho_{\phi_0}^{BAZS}}^{\rho_{\phi}^{BAZS}}\frac{d\left(3H_0^2 \left[\Omega_m (1+z)^3 + \Omega^{BAZS}_{\phi 0}(1+z)^{\frac{3\{w_{\beta}+\beta(1+w^{BAZS}_0)\}}{\beta}}exp\left[\frac{3w_{\beta}}{\beta^2}\{(1+z)^{-\beta} -1\}\right]\right] \right)}{\sqrt{ \frac{\kappa^2}{12}\left[\left(3H_0^2 \left[\Omega_m (1+z)^3 + \Omega^{BAZS}_{\phi 0}(1+z)^{\frac{3\{w_{\beta}+\beta(1+w^{BAZS}_0)\}}{\beta}}exp\left[\frac{3w_{\beta}}{\beta^2}\{(1+z)^{-\beta} -1\}\right]\right] \right)-6k\mu^4 (1+z)^2-\frac{3k\kappa^2 \mu^2 (1+z)^4}{16} \right]}}}~~~, \end{equation} \]

(53)

\( M^{BAZS}\) is the BH mass in BAZS at an arbitrary redshift (\( M^{BAZS}\) ) and \( M_0 (M(z=0))^{BAZS}\) is the present mass of the black hole in BAZS.

We plot \( log_{10}\left(\frac{M}{M_0}\right)\) for BAZS accretion as the previous cases. But the plots 4a-c are much different than previous cases. If \( w_{\beta}\) is positive with a positive \( \beta\) , mass plots for different \( w_m\) are followed not to differ much. They stayed only \( 10^{-6}\) to \( 10^{-6.5}\) orders of \( M_0\) at \( z\sim 3\) . Once negative \( w_{\beta}\) is chosen with small negative or positive \( \beta\) , huge differences for \( w_m\) are noticed. For \( w_0=-0.54\) , \( w_{\beta}=-3.58\) , \( \beta=3.5\) and \( w_m=-\frac{1}{3}\) , we follow a point of inflection to form around \( z\sim 0.5\) . This says \( w_m=-\frac{1}{3}\) is an exotic curve which led the mass of the accreting black hole to reach a peak but failed to accumulate more mass at that stage. As time rolled, DE grew more dominating and we again follow hikes in \( M\) as compared to \( M_0\) .

In figure 5a and 5b, we extrapolate the mass growth for CPL and BA parameterization universe. For this, we observe the mass of the central gravitating object to fall in future. Around \( z\sim -1\) , mass fall abruptly towards almost an order \( \sim 10^{-25}\) of its present mass. Universe filled up with BA type fluid also forces a black hole to loose its mass and around \( z\sim -1\) , \( 10^{-8}\) th order of present time black hole mass is left. These two results support the fate of black hole around a future singularity named “Big Rip” as mentioned by different articles [21].

While modification of gravity is almost a fashion by changing the Lagrangian of Einstein Hilbert action, Petr Horava proposed a quantum gravity where the ‘time’ is treated to be more fundamental quantum quantity than the general relativity. If a detailed balanced Horava gravity is chosen with projectibility condition, a black hole solution is given in this gravity theory by Kehagias and Sfetsos. We have studied a Michel-like accretion around such a black hole. We have considered the universe to be filled of dark energy which homogeneously fill the cosmos all over and exerts negative pressure to generate the late time cosmic acceleration as speculated by different type Ia supernovae observations since 1998. Among several kinds of dark energy candidates, redshift parametrization models to support the \( \Lambda CDM\) model at zero redshift and mimic other observations at different redshift models. In this article, four models - Chevallier-Polarski-Linder, Jassal-Bagala-Padmanabhan, Barboza-Alcaniz and Barboza-Alcaniz-Zhu-Silva are picked for accretion studies.

When a fluid is accreted towards a black hole, transonic accretion may get formed if the fluid speed equals the sonic speed for that concerned accreting fluid at some finite radial distance. Dark energy candidate, being re-pulsars, nullifies attractive power of black hole and hence critical points are shifted to nearer points of the gravitating core. Now different candidate possesses different kind of repulsion. That too is different for different sets of free parameters. Hence not only for separate dark energy candidate but also for a particular candidate, every individual set of free parameters generates a new model and hence a distinct critical distance. This gives us different values of an important constant \( \lambda_2\) which determines the growth of mass of the central gravitating object. Hence the mass accretion rate equation for each separate dark energy model is established and solved numerically to find the growing mass. To understand the situation deeply, we plot logarithm of the ratio of mass to present time mass for these models.

First, in Chevallier-Polarski-Linder universe we follow a Horava black hole to achieve \( \sim 10^{-6.5}\) to \( 10^{-5}\) th part of its present day mass at around \( z\sim 3\) . Mass growth has two distinct rates, first slowly and then faster. For Jassal-Bagala-Padmanabhan filled universe, the mass growth is slower than Chevallier-Polarski-Linder. At present time, in both Chevallier-Polarski-Linder and Jassal-Bagala-Padmanabhan cases mass growth does not depend much on the dark matter’s equation of state. In past (\( z\sim 3\) ) \( w_m=1\) and \( w_m=-\frac{1}{3}\) cases differ large. Among them difference for Chevallier-Polarski-Linder is even larger than Jassal-Bagala-Padmanabhan. Black holes in Barboza-Alcaniz type dark energy filled universe is followed to grow steeply in the neighbourhood of \( z\sim 0\) . For Barboza-Alcaniz-Zhu-Silva with negative \( w_0\) and \( w_{\beta}\) with positive \( \beta=3.5\) , for a redshift range (in recent past) mass growth stops if \( w_m=-\frac{1}{3}\) .

Next question arose in our mind regarding the fate of such black holes in future. We follow, if our universe is supposed to be filled up by dark fluids of type Chevallier-Polarski-Linder or Barboza-Alcaniz, mass of a black hole will reduce in future due to accretion of surrounding dark fluids. At \( z\sim -1\) mass is speculated to fall around \( -25^{th}\) order of the present day mass of it. This indicates that in future, black holes sitting in Horava Lifshitz Gravity and a dark energy background evaporates to almost an infinitesimal mass especially for Chevallier-Polarski-Linder type of dark energy model. Barboza-Alcaniz model supports evaporation but not upto minute remnant. For Chevallier-Polarski-Linder case, before we may reach the “Big Rip”, black holes are almost evaporated out such that the conflict regarding the fate of compact objects at future singularity even will not arise.

Acknowledgment: RB and FR thank IUCAA, Pune, India for providing Visiting Associateship. FR is also thankful to SERB, DST, DST, Govt. of India for financial support.

Data Availability Statement : No data was generated or analyzed in this study.

Conflict of Interest : There are no conflicts of interest.

Funding Statement : There is no funding to report for this article.