An Oscillating Holographic Dark Energy in $f(R)$ Gravity

The recent astronomical studies, including High-Redshift Supernovae, SNe Ia, Galaxy clustering, Cosmic microwave background and Large-scale structure, have confirmed that the universe is expanding at an accelerated rate [1, 2, 3, 4, 5]. To explore the accelerating and expanding nature of the universe, one approach is to modify general theory of relativity (GTR). Several modified theories of GTR such as $f(R)$, $f(T)$, $f(R,T)$, $f(G)$ [6, 7, 8, 9, 10, 11] and so on, have been proposed to rectify the problems in cosmology. The simplest modifications to GTR are the $f(R)$ theories, which introduces the arbitrary function of Ricci scalar in Einstein-Hilbert action. Firstly, Buchdahl [12] have proposed the $f(R)$ theory of gravity to generalize Einstein’s GTR. The modified theories have considered as a very interesting theories to describe the dark energy. The new class of modified $f(R)$ gravity have proposed by Nojiri and Odintsov [13] which unify $R^{m}$ early-time inflation with late-time $\Lambda CDM$ epoch. Wu and Zhu [14] has reconstructed the $f(R)$ theory which effectively describes the HDE and reconstruct the function $f(R)$ with the parameter $c<1,c=1$ and $c>1$ respectively. Many features of $f(R)$ theories and various applications to cosmology and gravity, such as dark energy, inflation, local gravity constraints, cosmological perturbations, and spherically symmetric solutions in weak and strong gravitational backgrounds, have been studied by Felice and Tsujikawa [15]. Katore [16] have investigated Bianchi-II, VIII and IX cosmological models with string in the framework of $f(R)$ theory of gravitation. Shamir [17] has investigated an anisotropic LRS Bianchi-I metric in $f(R)$ theory of gravity. Also, Santhi et al. [18] have investigated LRS Bianchi-I with bulk viscous string cosmological model in $f(R)$ gravity and shows that the realistic energy conditions, $\rho\geq 0$ and $\rho_{p}\geq 0$ are satisfied. Recently, An anisotropic cosmological model with magnetized strange quark matter have investigated by Ozdemir and Aktas [19] in $f(R)$ theory of gravity. Also, Katore and Gore [20] have studied $\Lambda$CDM Bianchi-I cosmological models with HDE in $f(R)$ gravity theory.

Another alternative for solving the dark energy problem is HDE, which is based on the Holographic concept initially presented by G’t Hooft [21] in the context of black hole physics. Without decaying into a black hole, this principle should be limited by an infrared cutoff scale $L$ and ultraviolet cutoff scale $\Lambda$. Fischer and Susskind [22] have proposed the cosmological version of holographic principle and have been discussed some of its consequences. The model of HDE is the idea of the short distance cut-off relates to the infrared cut-off has proposed by Li [23]. Setare and Vagenas [24] investigated the cosmic dynamics of the interacting HDE model. Adhav et al. [25] have investigated the homogeneous and anisotropic Bianchi-I cosmological model with interacting dark matter and HDE. Granda and Rojas [26] have studied an analysis of the autonomous system, corresponding critical points and their stability using HDE in the context of $f(R)$ gravity. Recently, Hatkar et al. [27] have investigated the Bianchi-III cosmological model with viscous HDE in Brans-Dicke theory of gravity. Shaikh and his collaborators [28] have studied HDE in $f(G)$ gravity.

The homogeneous and anisotropic Bianchi-II cosmological model plays a vital role in describing the universe at the early stages of its evolution. The homogeneous but anisotropic LRS Bianchi-II perfect fluid cosmological model with constant deceleration parameter have studied by Singh and Kumar [29]. Total energy of the universe has obtained by Aydogdu [30] based on LRS Bianchi-II models by using the Tolman, Papapetrou, and Weinberg energy-momentum complexes. Amirhashchi [31] has explored a novel anisotropic LRS Bianchi-II cosmological model using stiff fluid with variable $\Lambda$. Tiwari et al. [32] have studied the LRS Bianchi-II cosmological model in the presence of perfect fluid with varying $\Lambda$. Adhav [33] has used exponential and power-law expansions to find solutions to field equations for the LRS Bianchi-II cosmological model with anisotropic dark energy. Reddy and Santhi [34] have investigated the LRS Bianchi-II perfect fluid metric in the $f(R,T)$ theory of gravity. Sarkar [35] has investigated the correspondence between LRS Bianchi-II interacting HDE and tachyon scalar field dark energy model with variable deceleration parameter. Maurya et al. [36] have studied anisotropic LRS Bianchi-II DE models with a new approach of time-dependent deceleration parameter. R. L. Naidu et al. [37] have studied anisotropic Bianchi-V dark energy cosmological model with massive scalar field. Mete et al. [38] have examined the LRS Bianchi-II space time with interacting dark matter and HDE for $n=\frac{1}{2}$and $n=\frac{3}{2}$. Recently, Pradhan et al. [39, 40] have studied the LRS Bianchi-II cosmological models in GTR and the modified gravity presented by Saez–Ballester.

With the observation of above discussion, we have studied LRS Bianchi-II HDE cosmological models in the context of $f(R)$ gravity theory and explored the possibilty of bouncing solutions. In section 2, we derived the $f(R)$ gravity field equations for the LRS Bianchi-II cosmological model in the presence of HDE. The scale factor solutions in respect to the field equations determined in section 3. The section 4, devoted to analysed the physical and geometrical parameters. The energy conditions and its graphical interpretaion are discussed in section 5. The last section is concluded with discussion and conclusion.

The simplest and most popular modification of GTR is the $f(R)$ theory of gravity. For the $f(R)$ gravity, the modified action is given as follows:

where $f(R)$ is an arbitrary function of the Ricci scalar $R$ and $L_{m}$ is the matter Lagrangian. It is worth mentioning that the standard of Einstein-Hilbert action can be recovered when $f(R)=R.$

Now, the $f(R)$ gravity field equations are obtained by varying the action with respect to metric $g_{\mu\nu}$ are as follows:

where $F(R)\equiv\frac{df(R)}{dR}$, $\nabla_{\mu}$ denotes covariant differentiation, $T^{m}{}_{\mu\nu}$ and $T^{\Lambda}{}_{\mu\nu}$ are the energy momentum tensor corresponding to the matter and HDE respectively.

We have considered the LRS Bianchi-II space-time in an orthogonal frame given by [31] as

where $A$ and $B$ are functions of time $t$ only.

The scalar curvature $R$ for the metric (3) is given by

where the subscript ’4’ denotes differentiation with respect to time $t$.

The energy momentum tensors for matter and HDE are defined as

where $\rho_{m},$ $\rho_{\Lambda}$ are the energy densities of matter and HDE respectively and $p_{\Lambda}$ is the pressure corresponding to HDE.

The equation of HDE density is given by

where $\alpha,\beta$ are constants and $H$ is the Hubble parameter. This condition, proposed by Granda and Oliveros [41], which must fulfill the limitations imposed by current observations.

The continuity equations for matter and HDE can be obtained as

The equation of continuity for the matter is

The equation of continuity for the HDE is

The equation of state for HDE is defined as

Using equations (10) and (13) in (14), the EoS parameter becomes

The field equations (6)-(9) are a system of four equations in six unknowns $A$, $B$, $f$, $\rho_{m}$, $\rho_{\Lambda}$ and $p_{\Lambda}$. To find the solutions of the field equations, we consider the relations as:

1. 1.

   We suppose that the Shear scalar $\sigma$ is proportional to Expansion scalar $\theta$, resulting in a relation between scale factors $A$ and $B$ as

   $$B=A^{n},$$

   (16)

   where $n$ is an arbitrary real number.

2. 2.

   Secondly, for the simplicity, we assume the technique of Boutrous [42] and Ram [43] and the adhoc relation as

   $$\frac{A_{44}}{A}+\frac{A_{4}{}^{2}}{A^{2}}+\frac{A_{4}B_{4}}{AB}-\frac{B^{2}}{2A^{4}}=\frac{B_{44}}{B}+\frac{A_{4}B_{4}}{AB}+\frac{B^{2}}{2A^{4}}.$$

   (17)

Using equations (16) and (17), we have obtained as

We denote, $A_{4}=S(A)$ which implies that $A_{44}=SS^{\prime},$ where $S^{\prime}=\frac{dS}{dA}.$ Therefore, after solving and integrating above equation leads to

We assume that $n=2$ to find the deterministic solution in cosmic time. Therefore, above equation reduce to

On integrating above equation, we get

Therefore, using equations (16) and (21), we get

where $c_{1}$ and $c_{2}$ are integrating constants.

It is clear from equations (21) and (22) that, both the metric potentials $A$ and $B$ have oscillating nature and has physical meaning with Piao and Zhang [45]. Hence using equations (21) and (22), the LRS Bianchi-II space time for HDE can be written as

We define the Scale factor $a(t)$ and mean Hubble’s parameter $H$ for the LRS Bianchi-II metric (3) may be defined as

where $V=A^{2}B$ is the spatial volume of the universe.

Let us first define the dynamical scalars, such as the expansion scalar $\theta$, shear scalar $\sigma^{2}$, and mean anisotropic parameter $\Delta_{m}$:

where $H_{i}\,\left(i=1,\,2,\,3\right)$ represents the directional Hubble parameter in the directions of $x,y,z$ respectively. The mean anisotropic parameter of the expansion $\Delta_{m}$, plays an important role in determining whether the model is isotropic or anisotropic. It is the measure of deviation from isotropic expansion; when $\Delta_{m}=0,$ the universe expands isotropically.

Also, the deceleration parameter $q$ which is defined as

The cosmological behavior of the universe is predicted by the sign of $q.$ For an accelerating universe, the deceleration parameter $q$ is negative, while for a decelerating universe, it is positive.

Using euations (21) and (22), the value of scale factor $a(t)$ obtained as

The mean Hubble’s parameter $H$ as well as the deceleration parameter $q$ are obtained as

where $c_{2}$ is integrating constant.

It is observed that, for this model the scale factor vanishes throughout each cycle and Hubble parameter diverges $i.e.$ becomes singular at bouncing point. The bounce occurs when for $t=\left(\frac{n\pi}{2}-\frac{1}{2}\right)$, where $n$ is natural number which correspond to Big Crunch or Big bang singularity. Also, the deceleration parameter $q>0$ for any $t\geq 0.$ Hence, the positive sign of $q$ indicates that investigated model undergoes in the phase of deceleration and its behavior is depicted in figure 1.

Also, the physical quantities Shear scalar $\sigma^{2},$ Expansion scalar $\theta$ and mean anisotropic parameter $\Delta_{m}$ are obtained as

It is clear from the above equations that, the mean anisotropic parameter $\Delta_{m}$ is independent of time $t$, $i.e.$ it remains constant throughout the universe’s evolution. As a result, the current model appears to be anisotropic.

Therefore, from equations (33) and (34), we get

The ratio $\frac{\sigma}{\theta}={\rm constant}$ implies that the model under investigation does not approach isotropy. For the given metric (3), the energy densities for matter and HDE are obtained as

where $D,$ $E,$ $c_{2}$ and $c_{3}$ are integrating constants.

The holographic density $(\rho_{\Lambda})$ has oscillates in the positive quadrant as shown in figure 3. The density $\rho_{\Lambda}\geq 0$ is observed, indicating that the investigated model (23) satisfies the weak energy conditions. Our findings are similar to those of Beesham [46] and Rami [47].

Also, the pressure for HDE is obtained as

where $D$ and $E$ are constants of integration.

From figure 3 and 4, it is observed that, density $(\rho_{\Lambda})$ and pressure $(p_{\Lambda})$ of HDE diverges for $t=\left(\frac{n\pi}{2}-\frac{1}{2}\right)$, where $n$ is natural number $i.e.$ becomes singular at bouncing point. This singuarity occurs at bounce when the scale factor tends to zero and density as well as pressure becomes infinite called Type III singularity. The asymptotic behaviour of type III singilarity has observed by Setare et al. [55] for additive $g(T)$ solution.

From equations (14), (38) and (39), the HDE equation of state (EoS) parameter $\omega_{\Lambda}$ becomes

Figure 5 depicts the evolution of EoS parameter $\omega_{\Lambda}$ for HDE, which shows that it starts decreasing until it becomes negative, then increases until it becomes positive after a finite time and reverses back. As a result, the EoS parameter oscillates during the cycle resembles with [56] and varying between -0.778 and 1.016. $i.e.$ for any $t\geq 0\Rightarrow$ $\omega_{\Lambda}\in[-0.778,\,1.016].$ Durrer and Laukenmann [44] have investigated the oscillating universe as a possible alternative to inflation..

Using equations (4), (21) and (22), then Ricci scalar becomes

Solving equations (7), (8), (16) and (41), we get

where $D$ and $E$ are integrating constants.

Also, the matter density parameter $\Omega_{m}$ and HDE density parameter $\Omega_{\Lambda}$ are defined and obtained as

and

Therefore, using equations (43) and (44), the total energy density parameter $\Omega=\Omega_{m}+\Omega_{\Lambda}$ is obtained as

From figure 6, we have observed that the overall density parameter $\Omega$ oscillating during each cycle. During each cycle, it starts decresing from positive value to become negative and then goes incresing and become infinite. As a results, the values of overall density parameter $\Omega<1$, $\Omega>1$ and $\Omega=1$ represents the open, closed and flat universe respectively.

The energy conditions are a set of linear equations derived from Raychaudhuri’s equation that demonstrate that energy density can never be negative and gravity always attracts. The linear combination of pressure and density cannot be negative according to the energy conditions. Therefore, the energy conditions are expressed as follows:

1. (i)

   NEC $\Leftrightarrow$ $\rho_{\Lambda}+p_{\Lambda}\geq 0$

2. (ii)

   SEC $\Leftrightarrow$ $\rho_{\Lambda}+3p_{\Lambda}\geq 0$

3. (iii)

   DEC $\Leftrightarrow$ $\rho_{\Lambda}>\left|p_{\Lambda}\right|\geq 0$ and

4. (iv)

   WEC $\Leftrightarrow$ $\rho_{\Lambda}\geq 0$.

Using the equations (37) and (38), the expression for NEC and SEC are obtain as

The graphical representations of NEC, SEC, DEC and WEC depicted in figure (7) and (8) shows that, they obeys the state of energy conditions. It’s worth noting that our findings are more relevant than those of Ozdemir and Aktas [19] and Sahoo et al. [56], because our model satisfies all of the energy conditions.

We have investigated the HDE cosmological model based on LRS Bianchi-II space time in the context of $f(R)$ gravity. The LRS Bianchi-II cosmological model in $f(R)$ theory of gravity has the solutions from equation (16). It is observed that, our investigated model admits oscillating nature $i.e.$ it has cyclic nature. This investigated model represents the behaviour of a cyclic universe [49, 51] undergo as a sequence of continuous contraction and expansions phase [50, 52, 53, 54]. The mean anisotropic parameter $\Delta_{m}\neq 0$ which indicates that the investigated model is anisotropic throughout the evolution of the universe. Also, we have observed that the cosmological quantities like Shear scalar $\sigma$, Expansion scalar $\theta$ and Hubble parameter $H$ are oscillating infinitely and tends to infinity for $t=\left(\frac{n\pi}{2}-\frac{1}{2}\right)$ and all are vanishes at $t=\left(\frac{(2n-1)\pi}{4}-\frac{1}{2}\right)$, where $n$ is natural number. Also, the deceleration parameter $q>0$ implies that the decelerated expansion of the universe. It’s important to note that our investigation resembles with [16, 46, 56]. Moreover, all of the energy conditions, including the NEC, SEC, DEC, and WEC, are also satisfied.