# Analysis off(R) Theory Corresponding to NADE and NHDE

Authors: M. Sharif; M. Zubair
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

We develop the connection off(R) theory with new agegraphic and holographic dark energy models. The function f(R) is reconstructed regarding the f(R) theory as an effective description for these dark energy models. We show the future evolution of f and conclude that these functions represent distinct pictures of cosmological eras. The cosmological parameters such as equation of state parameter, deceleration parameter, statefinder diagnostic, and w−w′ analysis are investigated which assure the evolutionary paradigm of f.


---

## Body

## 1. Introduction

Expanding paradigm of the universe has been affirmed by the contemporary observational data [1–5]. The prime source behind this dramatic change in the evolution of the universe is said to be dark energy (DE). Dark energy is a strange type of gravitationally repulsive energy component, spreading over 72% of the contents in the universe. The nature of DE is still a question mark, and various representations have been proposed in general theory of relativity to understand it. Holographic dark energy (HDE) appeared as one of the most eminent candidates to address the issue of cosmic acceleration. The density of the HDE has been proposed by incorporating the mathematical form of the holographic principle as [6, 7]
(1)ρϑ=3c2Mp2L-2,
where c is a constant, Mp-2=8πG is the reduced Planck mass, and L is the infrared (IR) cutoff. Though Hubble horizon H-1 is the natural possibility for L, it does not imply the cosmic acceleration [8]. Li [7] suggested that the future event horizon is the most appropriate choice for IR cutoff, which seems to be consistent with recent measurements.The modification of the IR cutoff in HDE has been reported in different scenarios, such as introducing new time scale and consideringL as a function of the Ricci scalar in both original and generalized forms. Wei and Cai [9] suggested a new model of agegraphic DE by introducing conformal time as the time scale for the FRW universe and is known as new agegraphic DE (NADE). Wu et al. [10] discussed the evolution of the new agegraphic quintessence DE models in ω-ω′ phase plane both with and without interaction. The NADE model has been formulated in the context of alternative theories, such as Brans-Dicke theory [11] and Hor˘ava-Lifshitz gravity [12]. Granda and Oliveros [13, 14] proposed a new IR cutoff for HDE in terms of H and H˙ and discussed the correspondence of new HDE (NHDE) with models of scalar fields. This work has been extended for interacting case in nonflat universe [15].The modification of the Einstein-Hilbert action is another promising approach to explain the fact of cosmic acceleration. In this regard, there are various theories of gravity, such asf(R) [16–19], f(R,T), where T is the trace of the energy-momentum tensor [20–26], Gauss-Bonnet gravity [27, 28], and so forth. The action of f(R) theory with matter Lagrangian ℒM is defined as
(2)ℐ=∫‍dx4-g[Mp22f(R)+ℒ(M)].
In the literature [29–34], people have discussed the cosmological reconstruction of f(R) theory according to the class of HDE models. Capozziello et al. [29] developed an effective numerical scheme for reconstructing f(R) from Hubble parameter of a given DE model and applied this scheme to the quintessence DE model and Chaplygin gas.Following [29], Wu and Zhu [30] reconstructed f(R) according to HDE and explored the future evolution for different values of the parameter c. Feng [31] analyzed the effect of parameter α on reconstructed f(R) corresponding to Ricci DE. The explicit functions of f(R) in FRW universe can also be obtained from the reconstruction procedure according to the given DE model. Setare [32, 33] obtained f(R) functions corresponding to HDE and NADE by assuming an ansatz for the scale factor. In [34], reconstruction has been executed for both ordinary and entropy corrected models of holographic and NADE. In a recent work [25], we have reconstructed f(R,T) models according to holographic and NADE and found that the said models can represent the quintessence/phantom regimes of the universe.Here, we regard the NADE and NHDE as promising models and apply the numerical scheme for reconstructingf(R) without introducing any additional DE factor. The future evolution of f(R) is presented for different values of the essential parameters. We assure the evolution of f(R) by analyzing the corresponding behavior of cosmographic parameters in particular DE models. The paper has the following format: in Section 2, we reconstruct the f(R) theory according to NADE and discuss the future evolution. Section 3 provides the evolution of f(R) corresponding to NHDE. In Section 4, we summarize our findings.


## 2. Reconstruction from NADE

We consider the NADE density of the form [9]
(3)ρϑ=3n2Mp2η2,
where the factor 3n2 is inserted to parameterize some uncertainties, namely, the specific forms of cosmic quantum fields, the role of spacetime curvature, and so forth. And η is the conformal time in FRW background
(4)η=∫‍dta(t)=∫‍daHa2.For the flat FRW geometry comprising of matter component and NADE, the Friedmann equation is given by(5)3H2=ρM+ρϑ,
where ρM=ρM0e-3x from the energy conservation equation of matter. By defining the fractional matter and DE densities Ωϑ=ρϑ/ρcri,ΩM=ρM/ρcri with ρcri=3Mp2H2, the Hubble parameter H(x) is obtained as
(6)H(x)=H0(ΩM0e-3x1-Ωϑ)1/2.Differentiatingρϑ with respect to x and making use of DE conservation equation, the equation of state (EoS) parameter in NADE is obtained as
(7)ωϑ=-1+2Ωϑe-x3n.Using (3) and (4) with the relation Ωϑ=n2/H2η2, we obtain
(8)Ωϑ′=Ωϑ(1-Ωϑ)(3-2nΩϑe-x),
where prime denotes derivative with respect to x=lna. The initial condition on Ωϑ can be set from (5) as
(9)Ωϑ0+ΩM0=1.One can determineΩϑ using (8) and (9), and hence, the evolution of the universe in NADE can be executed.The field equations off(R) theory can be found by varying action (2) with respect to the metric
(10)Gαβ=Mp-2Tαβ(M)+Tαβ(curv),
where
(11)Tαβ(curv)=1fR[12gαβ(f-RfR)+fR;μν(gαμgβν-gαβgμν)]
originates from the curvature contribution to the effective energy-momentum tensor; the subscript R denotes derivative with respect to the scalar curvature and Tμν(M)=T^μν(M)/f(R). T^μν(M) is the standard matter energy-momentum tensor. For the flat FRW geometry, the respective field equations together with the conservation equation are given by
(12)3Mp2H2=ρT,-Mp2(2H˙+3H2)=pT,ρ˙T+3H(ρT+ωT)=0,
where ρT=ρM+ρcurv and pT=pM+pcurv. In this discussion, we consider the pressureless matter without any curvature-matter interaction. Equation (12) can be combined to single equation
(13)H˙=-12f[3H02ΩM0e-3x+(R¨-HR˙)fR+R˙2fRR].
Employing the relation d/dt=Hd/dx, we replace t by x, and hence, (13) can be translated into 3rd-order differential equation of f(x) as
(14)ℬ3(x)d3fdx3+ℬ2(x)d2fdx2+ℬ1(x)dfdx=-3H02ΩM0e-3x,
where ℬi are functions of H(x) and its derivatives given by (A.1).We aim to solve this equation to obtainf[R(x)] using the Hubble parameter H(x). For this purpose, we set the boundary conditions of the form [29]
(15)(dfdx)x=0=(dRdx)x=0,(d2fdx2)x=0=(d2Rdx2)x=0,f(x=0)=f(R0)=6H02(1-ΩM0)+R0.If the functionH(x) is known, then the coefficients ℬi and hence the function f(R) can be evaluated corresponding to the given DE model. In case of NADE, we do not have explicit form of H(x) whereas H(x) and its derivatives can be represented in terms of Ωϑ(x). Therefore, after lengthy calculations, the coefficients ℬi are interpreted in the form of Ωϑ(x) and Ωϑ′(x). We solve numerically the system of (8) and (14) together with the conditions (9) and (15).In [35], the constraint on parameter n is developed from the cosmological data for a flat universe consisting of DE and matter component, and the best fit value is found to be n=2.76-0.109+0.111. For the nonflat universe, Zhang et al. [36] found that the most appropriate measure of n from the WMAP 7-year observations is n=2.673-0.077-0.151-0.222+0.053+0.127+0.199. In this study, we set n=2.3,2.8,3.3,3.8 and ΩM0=0.27. For this choice of parameters, the function f(R) is plotted against R as shown in Figure 1. It is clear that functions appear distinct if |R| (or z=e-x-1) is large, while these functions seem to coincide for small |R|. The evolution of f shown in Figure 1 is quite similar to that for HDE [30]. We also plot these functions on lf-lR plane as shown in Figure 2, where lf=ln(-f) and lR=ln(-R). Our results are consistent with that in [29], and the parameter n appears as fundamental element in identifying the nature of NADE.Figure 1
Reconstructedf(R) for NADE with 0⩽z⩽10.Figure 2
Reconstructedf(R) for NADE in lf-lR plane with 0⩽z⩽10.To explore the effect ofn further, let us see the future evolution of |R|. Figure 3 shows that future variation of |R| is alike for different values of n and almost favors the cosmological constant. In Figures 4 and 5, we present the future evolution of f for n=2.3,2.8,3.3,3.8. For n=2.8, the curve seems to be similar to that for c=1 in HDE, but in the late stage of the universe, there would be a sudden change. The evolution of f is similar for other values of n which is shown in Figure 5. These plots indicate distinct characteristics for the NADE.Figure 3
Future evolution ofR in NADE with -1≲z⩽1.Figure 4
Future evolution off(R) for n=2.8 with -1≲z⩽2.Figure 5
Future evolution off(R) for different values of n with -1≲z⩽2.We extend our discussion and explore the evolution of NADE for the cosmographic parameters such as EoS parameter, deceleration parameter, statefinder diagnostic, andωϑ-ωϑ′ analysis. It is clear from (7) that NADE does not permit the crossing of phantom divide line ωϑ=-1, also if z→-1 and Ωϑ→1, then ωϑ→-1 in future evolution of NADE. Figure 6(a) shows that for our choice of parameter n, the EoS parameter of NADE favors the quintessence era, and in the late time, it mimics the cosmological constant regime. We also show the evolution of ωϑ′ in ωϑ-ωϑ′ plane for different values of n in NADE. Figure 6(b) depicts that the ωϑ-ωϑ′ plane represents the ΛCDM model (ωϑ=-1,ωϑ′=0) when z→-1 (or x→∞). The present values of ωϑ and ωϑ′ are denoted by dots on each curve. For n=2.3,2.8,3.3,3.8, the present values of (ωϑ,ωϑ′) are given by (-0.752,-0.218),(-0.796,-0.175),(-0.827,-0.147), and (-0.850,-0.126), respectively.Evolution trajectories in NADE for (a)ωϑ versus z and (b) ωϑ and ωϑ′ in ωϑ-ωϑ′ plane with -1≲z⩽2.
(a)
(b)Sahni et al. [37, 38] defined the statefinder diagnostic parameters {r,s} of the form
(16)r=a⃛aH3,s=(r-1)3(q-1/2).Introducing the EoS parameter and dimensionless density of DE, (16) is transformed as
(17)r=1-32Ωϑ[ωϑ′-3ωϑ(1+ωϑ)],s=-13ωϑ[ωϑ′-3ωϑ(1+ωϑ)].The deceleration parameterq in terms of ωϑ and Ωϑ is given by
(18)q=12(1+3ωϑΩϑ).The variation of deceleration parameterq with z for the NADE without interaction is shown in Figure 7(a). The transition of the universe from decelerating epoch to the accelerated era can be seen, and it will end up with q=-1 representing the de Sitter model. The sign flip of q depends on the selection of n; the era of cosmic acceleration starts earlier for small values of n as compared to larger values.Evolution trajectories in NADE for (a)q versus z and (b) the statefinder diagnostic in s-r plane with -1≲z⩽2.
(a)
(b)The plots of statefinder parameters in thes-r plane for n=2.3,2.8,3.3,3.8 are shown in Figure 7(b). The dots in the diagram correspond to present day values of statefinder parameters (r0,s0), which are denoted as (0.150,0.627)(red), (0.129,0.660)(black), (0.113,0.691)(gray), and (0.100,0.720) (blue). The evolution trajectories of the statefinder diagnostic are represented for the future evolution, and these will end up to star symbol {r=1,s=0}, the ΛCDM model. We also plot the statefinder diagnostic in q-r plane for our selection of parameter n together with the flat ΛCDM model. It can be seen from Figure 8 that evolution trajectories for NADE in q-r plane commence from the fix point (q=0.5,r=1), which represents the standard cold dark matter regime (SCDM). These curves end at (q=-1,r=1), the de Sitter model in future evolution of the universe. The past and future eras of the universe and present day values of (q,r) are represented by stars and dots, respectively.Figure 8
The statefinder plot for NADE inq-r plane. The green line represents the ΛCDM with -1≲z⩽2.


## 3. Reconstruction from NHDE

The energy density of HDE with Granda-Oliveros cutoff is given by [13, 14]
(19)ρϑ=3Mp2(μH2+υH˙),
where μ and υ are positive constants. Using this value of ρϑ, (5) can be written in the form
(20)E2(x)=22(1-μ)+3υΩM0e-3x+g0e2x(1-μ)/υ,
where E(x)=H(x)/H0, H0 is the present day value of Hubble parameter and g0 is the constant of integration which can be obtained using the condition E(x=0)=1 as
(21)g0=1-22(1-μ)+3υΩM0.Following [13, 14], the NHDE density is expressed as
(22)ρϑ=3H02[2μ-3υ2(1-μ)+3υΩM0e-3x+g0e2x(1-μ)/υ].The pressure of NHDE can be obtained using this value ofρϑ in the conservation equation of DE
(23)pϑ=-ρϑ-13dρϑdx=-3H02(2(1-μ)+3υ3υ)g0e2x(1-μ)/υ.Manipulating (22) and (23), the EoS parameter of NHDE turns out to be
(24)ωϑ=-[2(1-μ)+3υ]2g0e2x(1-μ)/υ3υ(2μ-3υ)ΩM0e-3x+[2(1-μ)+3υ]g0e2x(1-μ)/υ.Proceeding in a similar fashion as in the case of NADE, we obtain(25)𝒟3(x)d3fdx3+𝒟2(x)d2fdx2+𝒟1(x)dfdx=-ΩM0e-3x,
where 𝒟i are functions of E(x) and its derivatives, see the appendix. For the HDE with Granda-Oliveros cutoff, the expression for H(x) is directly useable in numerical computations. If we substitute (20) and constraint (21) in differential equation (25), then the resulting equation can be solved numerically under the boundary conditions (15).In [13, 14], the best fit values of parameters μ and υ are suggested as μ≈0.93 and υ≈0.5 to keep NHDE consistent with the theory of big-bang nucleosynthesis. Wang and Xu [39] developed the best fit values of parameters (μ,υ) in both flat and nonflat NHDE models from the current observational data. They found the best fit parameters for the flat model as μ=0.8502-0.0875-0.1064+0.0984+0.1299 and υ=0.4817-0.0773-0.0955+0.0842+0.1176. In this study, we select the parameters (μ=0.85,υ=0.48), (μ=0.93,υ=0.56), and (μ=1,υ=0.63). The plot of the function f(R) versus R for NHDE is shown in Figure 9(a). It shows that reconstructed f(R) for NHDE is the same for different values of parameters (μ,υ). We also plot these results on lf-lR plane as shown in Figure 9(b). To be more definite about the behavior of f, we draw plot on lf-x plane presented in Figure 10. This shows that different values of parameters do not affect the shapes of curves unless x>-0.3, and aftermath these curves would depict different picture.Reconstructedf(R) in (a) f-R plane and (b) lf-lR plane for NHDE with 0⩽z⩽10.
(a)
(b)Figure 10
Reconstructedf(R) for NHDE in lf-x plane with 0⩽z⩽10, where x=ln(1+z)-1.In order to explore the distinctive effect of parameters(μ,υ), we get insight of future evolution. First, we investigate the future evolution of R versus red shift which is shown in Figure 11. For (μ,υ)<(1,0.63), |R| would take infinitely large values in future which indicate the phantom era with ωϑ<-1 dominating over the matter part, leading to the big rip singularity. When (μ,υ)=(1,0.63), there is a slight variation in |R| and diagram assures the DE model with ωϑ=-1, the cosmological constant. The difference in the selected parameters (μ,υ) can be seen more effectively in the future evolution of f reconstructed according to the NHDE as shown in Figures 12, 13, and 14. For μ=0.85 and υ=0.48 in Figure 12, the curve shows that initially |R| decreases before reaching the present epoch (x=0), which changes its direction and it would increase leading to the phantom DE. In this scenario, |R| keeps growing, whereas f initially decreases and then attains positive value approaching to +∞. In fact, in phantom DE models, the point of reversion is a common characteristic because the DE components succeed in their competition with matter contents of the universe. For μ=0.93 and υ=0.56, we have almost identical picture as in Figure 12, but here the growing rate in f is comparatively large. For μ=1 and υ=0.63, we have a linear dependence of f on R leading to constant which is in accordance to the de Sitter model, where f(R)=R+Λ.Figure 11
Future evolution ofR in NHDE with -1⩽z⩽1.Figure 12
Future evolution off(R) for (μ,υ)=(0.85,0.48) with -1≲z⩽2. The dot denotes the present day and is the point of reversion for DE dominated model.Figure 13
Future evolution off(R) for (μ,υ)=(0.93,0.56) with -1≲z⩽2. The dot denotes the present day and is the point of reversion for DE dominated model.Figure 14
Future evolution off(R) for (μ,υ)=(1,0.63) with -1≲z⩽2. The dot denotes the present day and is the point of reversion for DE dominated model.Now, we discuss the evolution of the NHDE for the selected parameters(μ,υ) and interpret the behavior of EoS parameter, deceleration parameter, and statefinder diagnostic. The plot of EoS parameter for future evolution in NHDE is shown in Figure 15. It shows that the NHDE represents the de Sitter phase of the universe for (μ,υ)=(1,0.63). For (μ,υ)<(1,0.63), the EoS parameter intersects the phantom divide line (ωϑ=-1) and behaves as quintom model of DE [40]. In this perspective, ωϑ ends up with phantom era which may lead to cosmic doomsday when all the astronomical objects will be ripped apart. It is evident that domain of ωϑ in NHDE is consistent with the observational data of WMAP5 which establishes range of -1.11<ωϑ<-0.86 [41].Figure 15
Evolution trajectories of EoS parameter in NHDE for different values of parameters(μ,υ) with -1≲z⩽2. The dot denotes the present day, and star represents the ΛCDM model.The evolution ofq is represented in Figure 16 which confirms the behavior of ωϑ. The curve for (μ,υ)=(1,0.63) assures the ΛCDM model with q=-1. The transition from deceleration to accelerated epoch can be seen from this plot, and the values of redshift at the transition point are consistent with the observational results [42, 43]. The plot of statefinder diagnostic in NHDE for different values of parameters (μ,υ) is shown in Figure 17. The dot represents the fix point r=1,s=0 (i.e., the de Sitter phase), and all the curves pass through this point.Figure 16
Evolution trajectories ofq in NHDE for different values of parameters (μ,υ) with -1≲z⩽2. The dot denotes the present day, and star represents the ΛCDM model.Figure 17
Evolution trajectories of statefinder diagnostic in NHDE for different values of parameters(μ,υ) with -1≲z⩽2.


## 4. Conclusions

Thef(R) theory stands as one of the prosperous contexts to describe the cosmic evolution and the present day observational consequences. This theory appears to be a potential candidate in explaining the late time accelerated expansion. A profound model which can explain the cosmic evolution in a definite way is still under consideration. The cosmological reconstruction of f(R) gravity has been explored in [29–34], and the issue of which approach should be used is still alive. In [32–34], f(R) function corresponding to a class of HDE models has been constructed by assuming some ansatz for the scale factor in FRW background. A more effective scheme to reconstruct f(R) theory from the given evolution history H(z) is developed by Capozziello et al. [29]. In this scheme, the significant thing is that we can develop the correspondence of f(R) theory to the given DE model by using the expression of respective H(z). As a result, one can find the f(R) theory which explains the same dynamics (i.e., the cosmic evolution) as predicted by the given DE model. Now, it is of interest to consider the modified HDE models and address the intrinsic degeneracy among the f(R) theory and DE models. We find that predictions of both candidates (f(R) theory and DE models) reconcile as they represent distinct features of the same picture.In this work, we have reconstructed the functionf(R) according to NADE and NHDE in flat FRW geometry. The numerical reconstruction scheme is applied to obtain the evolution trajectories of f in different scenarios. In this reconstruction procedure, the Hubble parameter plays a significant role as we transform all quantities in terms of H and H˙. We summarize our results as follows.(i)
For NADE,H is given in terms of Ωϑ, so we solve the system of evolution equations for both Ωϑ and f. The results are shown in Figures 1 and 2 which are consistent with the constructed functions in the literature [29–31]. In comparison with HDE, the future variation of |R| and f shows identical behavior for different values of n. We can say that the behavior of f suggests the de Sitter phase in late time evolution of the universe. The cosmological parameters have been explored in NADE for n=2.3,2.8,3.3,3.8 to make sure the evolution of f. Figures 5–8 evidently show that NADE favors the quintessence regime, and in future evolution, it may end up with the de Sitter phase. Thus, our results for reconstructed f(R) are consistent with the independent evolution of NADE. We would like to emphasize that we have taken significantly different values of n, but all of these contribute similar results.(ii)
In case of  HDE with Granda-Oliveros cutoff, the Hubble parameter in the formE(x)=H(x)/H0 is directly used in numerical calculations. We have shown the function f in f-R and lf-lR planes in Figure 9. These curves seem to be identical for different values of parameters (μ,υ), and slight difference is found for lf-x plane which is shown in Figure 10. Further, we probe the future evolution of |R| and obtain distinct variations accordingly as (μ,υ)⩽(1,0.63). The future evolution of f in Figures 12–14 evidently shows the role of parameters (μ,υ). These plots represent distinct features of f which have been later confirmed by the evolution trajectories of ωϑ and q in Figures 15 and 16. For (μ,υ)=(0.85,0.48) and (0.93,0.56), it can be seen that f depicts the phantom DE era, and in such case ωϑ<-1 and q<-1. For (μ,υ)=(1,0.63), we have f representing the de Sitter phase with ωϑ=-1 and q=-1. Thus, our results for the function f corresponding to NHDE coincide with that of cosmographic parameters.It is to be noted that NADE and NHDE models are developed in the context of general relativity rather than any modified theory such asf(R) gravity. We have reconstructed f(R) by considering the curvature part as an effective description of these DE models. We also emphasize that this work is more comprehensive when comparing with previous ones as it involves the analysis of cosmological parameters to ensure the evolution of reconstructed function f(R).



---

*Source: 790967-2013-05-19.xml*