Next Article in Journal
Optimization of Energetic Train Cooperation
Next Article in Special Issue
A Confront between Amati and Combo Correlations at Intermediate and Early Redshifts
Previous Article in Journal
Path Mapping Approach for Network Function Virtualization Resource Allocation with Network Function Decomposition Support
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamical Properties of Dark Energy Models in Fractal Universe

1
CAMS, Business School, Faculty of Management Studies, University of Central Punjab, Lahore 54000, Pakistan
2
Department of Mathematics, Govt. College University, Faisalabad 38000, Pakistan
3
Department of Mathematics, COMSATS University Islamabad, Lahore-Campus, Lahore 54000, Pakistan
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(9), 1174; https://doi.org/10.3390/sym11091174
Submission received: 30 July 2019 / Revised: 8 September 2019 / Accepted: 9 September 2019 / Published: 16 September 2019
(This article belongs to the Special Issue Beyond the Standard Cosmological Model in the Multi-messenger Era)

Abstract

:
In this paper, we consider the flat FRW spacetime filled with interacting dark energy and dark matter in fractal universe. We work with the three models of dark energy named as Tsallis, Renyi and Sharma–Mittal. We investigate different cosmological implications such as equation of state parameter, squared speed of sound, deceleration parameter, statefinder parameters, ω e f f ω e f f (where prime indicates the derivative with respect to ln a , and a is cosmic scale factor) plane and Om diagnostic. We explore these parameters graphically to study the evolving universe. We compare the consistency of dark energy models with the accelerating universe observational data. All three models are stable in fractal universe and support accelerated expansion of the universe.

1. Introduction

The presence of networks of voids and structures in distribution of galaxies in space is one of the remarkable findings in the last two decades. The earlier pervasive examination of galaxy distribution with focus only on angular coordinates was unable to locate such high scale structures. Such findings refuse to accommodate the standard cosmological theories, in which homogeneity at large scale is of immense importance [1]. The existence of such a large scale at which the matter distribution approaches an average density is a matter of intense debate [2,3,4]. At small scales, it is well established that the distribution of galaxies is fractal, although its deviation from homogeneity is still debatable. Some statistical methods have been used to analyze the three dimensional data with some disagreements [3,4]. The fractal scenario is still resisted due to its incompatibility with the standard theories and isotropy of microwave back ground radiation [5,6]. The fractal distribution for matter was suggested in the framework of an expanding universe, the Friedmann cosmology [7].
A fractal is an identical and fundamental altering distribution of points at all scales, excluding its gravitational dynamics in framework of FRW solutions to general relativity [8,9]. The problem arises because of incompatibility of a fractal with the Cosmological Principle, which requires matter distribution to be isotropic and homogeneous [5]. In the case of non-analytical structure, such as fractal, all points are statistically equivalent, satisfying the conditional cosmological principle [8,9,10]. Systems whose effective dimensionality changes with scale may have fractal behavior. To understand the fractal effects, the present work is an effort to analyze cosmic acceleration under influence of fractal effects.
Many pieces of observational evidence such as different luminosity type Ia supernovae (SNIa) [11], Large Scale Structure [12] and Cosmic Microwave Background (CMB) radiation [13] strongly support the accelerated expansion of our universe. The candidates responsible for this expanding scenario provoked a new type of matter violating strong energy condition, i.e., ρ + p < 0 . Such exotic matters with negative energy density are named as Dark Energy (DE) [14]. One amongst the proposed candidates of DE appealing researchers the most is Holographic Dark energy (HDE). The idea behind this model is the Holographic principle, which emphasizes that for a physical system the number of degrees of freedom should be scaled by its bounding area instead of its volume [15,16] constrained by infrared cutoff [17]. By using this idea, Li [18] recommended the following constraint on its energy density: ρ d 3 C 2 M p 2 / L 2 , where the equality sign holds only for saturated holographic bound. Here, C 2 is the dimensionless constant, L denotes the IR-cutoff radius and M p 2 = ( 8 π G ) 1 is the reduced Plank mass. Li used three options for IR cutoffs, namely Hubble horizon, the future event horizon and the particle horizon. According to Li, only the future event horizon has the capability of accommodating the sufficient acceleration of universe. Fischler and Susskind proposed HDE models based on the cosmological state of Holographic principal [19]. The derivation of holographic energy density ρ d = 3 C 2 M p 2 / L 2 is based on the entropy–area relationship S ~A ~ L 2 of black holes. Here, A = 4 π L 2 is area of horizon. A more general representation of HDE model with NO cutoff was presented Nojiri and Odinstov, who studied its cosmological and thermodynamical aspects [20,21].
Sheykhi [22] developed HDE model by using Hubble horizon as IR cutoff and believed that this model is capable of illustrating present state of universe by considering interaction of DE and Dark matter (DM). Hu and Ling examined the interacting HDE model and cosmological parameters through observational constraints [23]. In the low red-shift region, the HDE model supports the present observations. For a specific choice of interacting term, they attempted to reduce the cosmic coincidence problem by using different options for time rate of change of ratio of DM to HDE densities. By using a number of recent observational schemes, Ma et al. [24] observed big rip singularity for different parameters while exploring interacting and non-interacting HDE with DM. They found that HDE models are more compatible with observations as compared to Λ CDM Model.
The DE and horizon entropy has an apparent influence on each other in thermodynamics. The exotic nature of space time and the Bekenstein entropy being non-extensive entropy measure, the generalized entropies, i.e., Tsallis and Renyi entropies, have been used to examine cosmological and gravitational phenomenon [25,26,27,28]. Amongst several generalized entropy formalisms to study cosmological phenomena, Tsallis and Renyi entropies develop more precise pictures of the universe. Another generalized entropy called the Sharma–Mittal measure [25] has been used to propose a new model of HDE [29]. In this model, the role of IR-cutoff is played by the Hubble horizon. It is stable and compatible with universe expansion history.
In the present work, we study the cosmological implications of HDE in fractal universe. We work in a flat FRW universe filled with interacting DM and DE. We use three different HDE models, namely Tsallis Holographic dark energy (THDE), Renyi Holographic dark energy (RHDE) and Sharma–Mittal Holographic dark energy (SMHDE), for this purpose. These models of HDE can also be generated from extended generalized form of EoS as presented by [21] and also presented in [30]. For analysis, we take help of EoS parameter, square speed of sound and deceleration parameter. Further, we make use of ω e f f ω e f f , statefinders { r , s } and Om diagnostics. In Section 2, the field equations under selected scenario and expressions for all the used HDE models are given. Section 3 contains discussion regarding deceleration parameter for the three HDE models in fractal universe. Similarly, Section 4, Section 5, Section 6, Section 7 and Section 8 contain the views about EoS parameter, square speed of sound, statefinder parameter, Om diagnostic and ω e f f ω e f f , respectively. The last section comprises conclusion after detailed discussions of the previous sections.

2. Field Equations with Solutions

In a Fractal universe, time and space coordinates scale identically, i.e., [ x μ ] = 1 , μ = 0 , 1 , , D 1 . The non-trivial measure (as in Lebesgue–Stieltjes integrals): d D x d ρ ( x ) , [ ρ ] = D α D replaces the standard measure in action. Here, D is the topological dimension of the embedding space-time. During the evolution of the system, for a given time, a preserved fraction of states is represented roughly by the parameter α > 0 . Such structures are obtainable by varying effective dimensionality of universe at different scales. This attribute can be acquired by introducing fractal action [31,32]. The general Boral probability measure ρ is assumed as measure on Fractal set. Therefore, the metric-spacetime M furnished with measure ρ is represented by ( M , ρ ) in D dimensions. In this case, ρ is continuous with d ρ = d D x ν ( x ) , some multidimensional Lebesgue measure, and ν is the weight function (fractal function). The total action of Einstein gravity in fractal spacetime can be given as
S = S g + S m .
The gravitational part of action is
S g = M p 2 d D x ν ( x ) g ( R 2 Λ ω μ ν μ ν ) ,
and
S m = d D x ν ( x ) g L m ,
is action of matter part minimally coupled with gravity [31,32,33]. Here, M p = ( 8 π G ) 1 2 is the reduced Planck mass and ω represents fractal parameter, R is the Ricci scalar.
In homogeneous, isotropic FLRW spacetime model, the Friedmann equations in a fractal universe can be obtained by taking variation of action with respect to metric tensor as given below
H 2 + k a 2 + H ν ˙ ν ω 6 ν ˙ 2 = 1 3 ( ρ d + ρ m ) + Λ 3 .
Here, ρ d and ρ m are energy densities of DE and DM, respectively, H = a ˙ a is Hubble parameter. k is the curvature constant with k = 1 , 0 , 1 for open, flat and closed universe, respectively. We take time-like fractal profile ν = a γ with assumptions Λ = 0 , k = 0 and γ = constant. Now, the Equation (2) takes the following form
H 2 1 γ ω 6 ν 2 γ 2 = 1 3 ( ρ d + ρ m ) .
The energy conservation equation in Fractal universe is given as
ρ ˙ + ( ρ + p ) ( 3 H + ν ˙ ν ) = 0 .
After substituting ν = a γ , Equation (4) can be written as
ρ ˙ + H ( ρ + p ) ( 3 γ ) = 0 .
We introduce the phenomenological term Q to represent interaction between dark matter and dark energy model. Here, the decay rate should be proportional to the present value of Hubble parameter H as determined by Supernova and CMB data. The interaction term Q describes the energy flow between the two fluids. Observations from [34,35] give clear picture of interaction between DE and DM. The supernova data together with CMB and large scale structure [36] describe that such an interaction could be inferred from the expansion history of the universe. The relevant studies on the interaction between DE and DM have been carried out in [37,38,39]. The transition from matter domination to DE domination is given by an appropriate energy exchange rate. An interaction term having decay rate proportional to present value of Hubble parameter is assumed to be a good fit to the expansion history of the universe. The evidence for such interacting DE has also been provided by Baryon Oscillation Spectroscopic Survey (BOSS) [40]. Some authors [41] have also discussed various DE models in different frameworks and found interesting results. Amongst the different forms of Q available in the literature, we choose the following form of Q as in [42]
Q = 3 ξ H ρ m .
Here, ξ (may be positive or negative) is the interaction parameter. Now the conservation equation takes the following form
ρ d ˙ + H ( ρ d + p d ) ( 3 γ ) = Q , ρ m ˙ + H ( ρ m ) ( 3 γ ) = Q ,
with p m = 0 being pressureless fluid. From Equation (7), we can easily find
ρ m = ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) .
Next, we discuss the energy densities for HDE models under consideration.

2.1. THDE

The Tsallis definition of entropy plays a vital role in the study of gravitational and cosmological systems. The Tsallis entropy content of system is the power law function of system area confirmed by the quantum gravity considerations. The horizon entropy of black hole can be modified as shown by Tsallis and Cirto [43].
S δ = γ A δ .
Here, δ is the non additivity parameter [43] and γ is unknown constant and A = 4 π L 2 representing the area of horizon. The Bekenstein entropy is recoverable at the limit δ = 1 and γ = 1 4 G with = k b = c = 1 . At this limit, the system is describable by ordinary distribution of probability. By using Equation (9) and relation L 3 Λ 3 ( S δ ) 3 4 described by Cohen et al. [17], we get Λ 4 γ ( 4 π ) δ L 2 δ 4 . This gives the following form to energy density of Tsallis holographic dark energy (THDE) [44].
ρ d = B H 4 2 δ .
By using Equations (3) and (7), we can calculate H ˙ and p d as
H ˙ = ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ ,
p d = 3 ξ ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) B H 4 2 δ + B ( 4 2 δ ) H 2 2 δ γ 3 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ .

2.2. RHDE

Renyi (S) and Tsallis ( S T ) are two well-known one-parameter generalized entropies defined as
S = 1 δ l n Σ i = 1 W P i 1 δ ,
S T = 1 δ Σ i = 1 W ( P i 1 δ P i ) .
Here, δ = 1 Q ; combining the above one-parametric entropy measures with each other, we can get the following relation as in [45,46,47,48]
S = 1 δ l n ( 1 + δ S T ) .
Using the above equation and assumption ρ d d V T d S , Renyi holographic dark energy (RHDE) can be defined as
ρ d = 3 C 2 H 2 8 π ( 1 + δ π H 2 ) .
The following expression for H ˙ and p d can be obtained by using Equations (3) and (7):
H ˙ = ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ,
p d = 1 H ( 3 γ ) ( 3 ξ H ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + 3 C 2 H 3 ( 3 γ ) 8 π ( 1 + δ π H 2 ) + 2 F H 3 ( H 2 + 2 δ π ) ( H 2 + δ π ) 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ) .

2.3. SMHDE

Another generalized entropy measure introduced by Sharma and Mittal defined as [49,50]
S S M = 1 1 r Σ i = 1 W P i 1 δ 1 r δ 1
Using Equations (13)–(15), one can easily write
S S M = 1 R 1 + δ S T R δ 1 ,
with R = 1 r . As Bekenstein–Hawking entropy is supposed to be an appropriate candidate for Tsallis entropy, we can replace S T in Equation (20) by S B = A 4 with A being the horizon area and get
S S M = 1 R 1 + δ A 4 R δ 1 ,
Using Equation (21) and Λ 4 S L 4 as described above, we approach SMHDE density given below
ρ d = 3 C 2 H 4 8 π R 1 + π δ H 2 R δ 1 ,
with C 2 being free parameter. With the help of Equations (3), (7) and (22), we derive the following expressions for H ˙ and p d
H ˙ = ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ,
p d = 1 H ( 3 γ ) ( 3 ξ H ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) 3 C 2 H 5 ( 3 γ ) 8 π R ( 1 + δ π H 2 ) β 1 3 C 2 8 π R ( 4 H 3 1 + δ π H 2 β 1 ) 2 δ π β H 1 + δ π H 2 β ) ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) ( C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 ) 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 .

3. Deceleration Parameter

To determine the cosmological acceleration of any cosmological model, the deceleration parameter plays a vital role. A negative value of deceleration parameter q < 0 exhibit cosmic acceleration of model under discussion, while positive q > 0 shows a decelerating universe. The deceleration parameter can be calculated as
q = 1 H ˙ H 2 .

3.1. THDE

In the case of THDE, by using the expression of H ˙ from Equation (11) in Equation (25), we get the following expression for deceleration parameter q
q = 1 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H 3 ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 5 2 δ .
In Figure 1, we plot the variation of deceleration parameter with respect to redshift parameter with in the range 0.9 z 0.1 . Here, we take ω = 0.9 , ξ = 2 , B = 1 , ρ 0 = 1 , γ = 2 , α = 1 and δ = 2 . From the plot, we see that the deceleration parameter remains negative before z = 0.15 , exhibiting the cosmic acceleration. After z = 0.15 , the curve representing q becomes positive. For the redshift parameter 0.1 z 0.15 , q represents decelerated phase and then for 0.15 z 0.9 , q exhibits accelerated expansion. From the curve, it is evident that the universe is evolving from early decelerated phase towards late accelerated phase.

3.2. RHDE

The following expression of deceleration parameter for RHDE can be obtained by substituting value of H ˙ from Equation (17) into Equation (25)
q = 1 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ν 2 6 H 3 ( 1 γ ω 6 ν 2 γ 2 ) 2 F H 5 H 2 + δ π 2 H 2 + 2 δ π .
Figure 2 contains the plot of the expression for deceleration parameter versus redshift parameter 0.9 z 0.1 for RHDE. Here, we use same values of constants as discussed above. We observe that the deceleration parameter q remains negative before z = 0.14 and becomes positive after z = 0.14 . It represents accelerated universe before about z = 0.14 and decelerated phase of universe after z = 0.14 . This means our universe is going through an accelerated phase of expansion.

3.3. SMHDE

To get deceleration parameter expression in SMHDE, we replace H ˙ in Equation (25) by Equation (23) and achieve the following
q = 1 H 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 .
The plot of the expression for deceleration parameter versus redshift parameter 0.9 z 0.1 for SMHDE is given in Figure 3. Under the same above-mentioned values for constants along with C = 1 and β = 800 , the behavior of deceleration parameter is identical with those in THDE and RHDE. It shows decelerated phase of the universe in the early epoch and eventually enters the era of accelerated expansion of the universe.

4. Eos Parameter

EoS parameter is mainly used for categorizing DE models, i.e., for non-relativistic matter we have ω e f f < 0 . Likewise, for quintessence, we must have 1 < ω e f f < 1 3 . The cosmological model is represented by ω e f f = 1 and the phantom barrier holds for ω e f f < 1 . Similarly, ω e f f > 1 corresponds to quintom behavior. To achieve the EoS parameter, we use the following equation
ω e f f = p e f f ρ e f f .

4.1. THDE

By substituting values from Equations (10) and (12) into Equation (29), we get the following expression for square speed of sound.
ω e f f = ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + B H 4 2 δ 1 ( 3 ξ ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) γ 3 B H 4 2 δ + ( B ( 4 2 δ ) H 2 2 δ γ 3 ) ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ ) .
Figure 4 gives the plot of EoS parameter ω e f f for THDE for redshift parameter 0.9 z 0.2 with all the constants bearing same values as mentioned above. Following the classification described previously, we interpret the plot as below. The EoS parameter attains the value 1 from z = 0.9 to z = 0.7 , showing that THDE is behaving as cosmological constant. Afterwards, it behaves as quintessence as it takes values within the range 1 < ω e f f < 0.33 when redshift parameter is in range 0.7 < z < 0.27 . For z > 0.27 , THDE shows quintom behavior.

4.2. RHDE

To get an expression for EoS parameter in RHDE, we substitute values into Equation (29) from Equations (16) and (18)
ω e f f = 1 H ( 3 γ ) ( 3 ξ H ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + H ( 3 γ ) F H 4 ( H 2 + δ π ) 1 + 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H ω ν 2 γ 2 H 3 2 H 1 γ ω 6 ν 2 γ 2 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 1 ) ( F H 2 ( H 2 + δ π ) 1 + ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) ) 1 .
The plot of EoS parameter ω e f f for RHDE against redshift parameter 0.9 z 0.3 is given by Figure 5. Here, all the constants have same values as mentioned above. Here, RHDE behaves as cosmological constant as the EoS parameter attains the value 1 from z = 0.9 to z = 0.8 . Then, for redshift parameter z > 0.8 , the curve lies in the range 1 < ω e f f < 0.33 , showing the quintessence behavior.

4.3. SMHDE

To obtain the expression for EoS parameter in the case of SMHDE as given below we substitute values from Equations (22) and (24) into Equation (29)
ω e f f = 1 H ( 3 γ ) ( 3 ξ H ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) 3 C 2 H 5 ( 3 γ ) 8 π R ( 1 + δ π H 2 ) β 1 3 C 2 8 π R ( 4 H 3 1 + δ π H 2 β 1 ) 2 δ π β H 1 + δ π H 2 β ) ( ρ 0 3 ( 3 ξ + γ 3 ) H ( 1 + z ) ( 3 ξ + γ 3 ) ω ν 2 γ 3 H 3 3 ) ( 2 H ( 1 γ ω ν 2 γ 2 6 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 3 C 2 H 4 8 π R ( 1 + δ π H 2 ) β 1 + ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) 1 .
The graph for EoS parameter in case of SMHDE is given in Figure 6. Here, we use same previously described vales of constant for redshift parameter 0.35 z 0.26 . For redshift parameter 0.35 z 0.29 , SMHDE acts as quintessence and after z = 0.29 it shows quintom-like behavior.

5. Square Speed Of Sound

For stability analysis of HDE models, we make use of the square speed of sound C s 2 . We calculate it as follows
C s 2 = p e f f ρ e f f .
Here, p e f f = p d + p m and ρ e f f = ρ d + ρ m . For a HDE model to be stable, we must have C s 2 > 0 .

5.1. THDE

The expression for square speed of sound in THDE is as given below
C s 2 = ( 3 ξ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) γ 3 H + B ( 4 2 δ ) H 2 2 δ 3 ( γ 3 ) ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) 5 H 2 ω γ 3 ( 1 + z ) 2 γ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ + ρ 0 ( 3 ξ + γ 3 ) 2 ( 1 + z ) ( 3 ξ + γ 3 ) H 2 3 + 2 ω γ 4 H 4 ( 1 + z ) 2 γ 3 2 1 γ ω 6 ( 1 + z ) 2 γ γ 2 B 3 ( 4 2 δ ) ( 3 2 δ ) H 2 2 δ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ 2 ) 2 H 1 γ ω 6 ( 1 + z ) 2 γ γ 2 B 3 ( 4 2 δ ) ( 3 2 δ ) H 2 2 δ 1 + ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ 2 B ( 4 2 δ ) ( 2 2 δ ) γ 3 H 1 2 δ B ( 4 2 δ ) H 3 2 δ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ ) ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H + B ( 4 2 δ ) H 3 2 δ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ ) 1 .
In Figure 7, we plot the expression for C s 2 versus redshift parameter 0.7 z 0 with same values of constants. The square speed of sound remains positive. This exhibits the stability of interacting THDE model in fractal universe.

5.2. RHDE

For RHDE, the square speed of sound has the following expression
C s 2 = 3 ξ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H γ 3 2 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) + 1 γ 3 ( H 1 2 H ( 1 γ ω 6 ν 2 γ 2 ) 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 1 2 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ( ρ 0 3 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) 2 ω 3 H 4 γ 4 ν 2 + 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 ( F H 2 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 4 3 F H 4 ( H 2 + δ π ) 3 ( H 2 + 2 δ π ) + 2 3 F H 4 ( H 2 + δ π ) 2 ( 1 γ ω 6 ν 2 γ 2 ) ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 3 H 2 ω γ 3 ν 2 H 2 ω γ 3 ν 2 ρ 0 3 ( 3 ξ + γ 3 ) 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) + 2 F ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ) 2 ( 3 H 2 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 4 H 4 ( H 2 + δ π ) 3 ( H 2 + 2 δ π ) + 2 H 4 ( H 2 + δ π ) 2 ) ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 F H ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 ) .
The square speed of sound for RHDE is given in Figure 8 against the redshift parameter. Here, we take ξ = 2 , α = 1 and the other constants are same as described above. It is evident from plot that square speed of sound remains positive for 0.46 z 0.454 . Hence, the RHDE model in fractal universe is stable.

5.3. SMHDE

The square speed of sound has the following expression in the case of SMHDE:
C s 2 = 1 3 γ ( 3 ξ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H 1 3 C 2 8 π R ( 12 H 2 1 + δ π H 2 β 1 10 δ π β 1 + δ π H 2 β 1 + 4 δ 2 π 2 β ( β 1 ) 1 + δ π H 2 β 1 H 2 ) ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( ( 1 γ ω 6 ν 2 γ 2 ) 2 H C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) 3 + 3 C 2 8 π R 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ( 2 H ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ) 1 ( 2 3 ω γ 4 ν 2 H 4 + 1 3 ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) ( ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) + ρ 0 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) 2 ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ( 12 H 2 1 + δ π H 2 β 1 10 δ π β 1 + δ π H 2 β 1 + 4 δ 2 π 2 β ( β 1 ) 1 + δ π H 2 β 1 H 2 ) ( ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R ( 4 H 3 ( 1 + δ π H 2 ) β 1 ) 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) 2 5 3 ω γ 3 ν 2 H 2 ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) + ( H 2 3 + γ ) 3 C 2 8 π R ( 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) ( 3 C 2 8 π R 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R ( 4 H 3 ( 1 + δ π H 2 ) β 1 ) 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) 1 .
Figure 9 contains the plot of the square speed of sound for SMHDE, with all values of constants as mentioned above and ξ as in the previous plot with β = 800 . Here, the trajectory shows the stability of SMHDE model from z = 0.8 to z = 0.5 .

6. Statefinder Parameters

To obtain a diagnostic tool involving higher derivative of scale factor, Sahni et al. [51] introduced the state finder pair { r , s } . It is a geometrical diagnostic being depending on scale factor. Mainly, it is used to recognize DE models. Distinct evolutionary trajectories are shown in { r , s } plane by different DE models. Trajectories in s-r plane relevant to different cosmological models show qualitatively different behaviors. By using statefinder diagrams, we can explicitly diagnose the spatial curvature contribution in different DE models. This statefinder parameter { r , s } in terms of Hubble parameter can be defined as
r = 1 + 3 H ˙ H 2 + H ¨ H 3
and
s = 3 H H ˙ + H ¨ 3 H ( 2 H ˙ + 3 H 2 ) .

6.1. THDE

By using Equation (11) in Equations (37) and (38), we get the following expressions for r and s:
r = 1 + 3 H 2 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ + 1 H 3 ( ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) 5 H 2 ω γ 3 ( 1 + z ) 2 γ 3 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ + ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H 2 + 2 ω γ 4 3 ( 1 + z ) 2 γ ( 2 ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 B 3 ( 4 2 δ ) ( 3 2 δ ) H 2 2 δ ) ) ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ 2 ) 2 1 γ ω 6 ( 1 + z ) 2 γ γ 2 B 3 ( 4 2 δ ) ( 3 2 δ ) H 2 2 δ 1 ,
s = 1 3 ( 3 H 2 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ + 1 H 3 ( ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) 5 H 2 ω γ 3 ( 1 + z ) 2 γ 3 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ + ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H 2 + 2 ω γ 4 3 ( 1 + z ) 2 γ ( 2 ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 B 3 ( 4 2 δ ) ( 3 2 δ ) H 2 2 δ ) ) ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ 2 ) 2 1 γ ω 6 ( 1 + z ) 2 γ γ 2 B 3 ( 4 2 δ ) ( 3 2 δ ) H 2 2 δ 1 ) ( 1 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ 6 H 3 ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 5 2 δ ) 1 .
Figure 10 shows variation of statefinder parameter s against r in the case of THDE. It is evident from the plot that we have r = 1 for s = 0 . This corresponds to the standard Λ CDM cosmological model for accelerating universe.

6.2. RHDE

In the case of RHDE, Equations (37) and (38) can be modified using Equation (17) to get the following expressions:
r = 1 + 3 H 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) + H 3 2 H ( 1 γ ω 6 ν 2 γ 2 ) 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 1 ( ρ 0 3 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) 2 ω 3 H 4 γ 4 ν 2 + 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 ( F H 2 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 4 3 F H 4 ( H 2 + δ π ) 3 ( H 2 + 2 δ π ) + 2 3 F H 4 ( H 2 + δ π ) 2 ( 1 γ ω 6 ν 2 γ 2 ) ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 3 H 2 ω γ 3 ν 2 H 2 ω γ 3 ν 2 ρ 0 3 ( 3 ξ + γ 3 ) 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) ,
s = 1 3 ( 3 H 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) + H 3 2 H ( 1 γ ω 6 ν 2 γ 2 ) 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 1 ( ρ 0 3 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) 2 ω 3 H 4 γ 4 ν 2 + 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 ( F H 2 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 4 3 F H 4 ( H 2 + δ π ) 3 ( H 2 + 2 δ π ) + 2 3 F H 4 ( H 2 + δ π ) 2 ( 1 γ ω 6 ν 2 γ 2 ) ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 3 H 2 ω γ 3 ν 2 H 2 ω γ 3 ν 2 ρ 0 3 ( 3 ξ + γ 3 ) 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) ) 3 2 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ν 2 6 H 3 ( 1 γ ω 6 ν 2 γ 2 ) 2 F H 5 H 2 + δ π 2 H 2 + 2 δ π 1 .
Figure 11 gives plot of statefinder parameter s against r in the case of RHDE. It is clear from the plot that we have r = 1 for s = 0 . This corresponds to the standard Λ CDM cosmological model for accelerating universe.

6.3. SMHDE

For SMHDE, we customize Equations (37) and (38) using Equation (23) to attain the following expressions:
r = 1 + 3 H 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 + H 3 ( 2 H ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ) 1 ( 2 3 ω γ 4 ν 2 H 4 + 1 3 ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 + ρ 0 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) 2 ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R ( 12 H 2 1 + δ π H 2 β 1 10 δ π β 1 + δ π H 2 β 1 + 4 δ 2 π 2 β ( β 1 ) 1 + δ π H 2 β 1 H 2 ) ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) 2 5 3 ω γ 3 ν 2 H 2 ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) ,
s = 1 3 ( 3 H 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 + H 3 ( 2 H ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ) 1 ( 2 3 ω γ 4 ν 2 H 4 + 1 3 ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 β 1 ) 1 + ρ 0 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) 2 ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R ( 12 H 2 1 + δ π H 2 β 1 10 δ π β 1 + δ π H 2 β 1 + 4 δ 2 π 2 β ( β 1 ) 1 + δ π H 2 β 1 H 2 ) ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) 2 5 3 ω γ 3 ν 2 H 2 ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) ) ( H 2 ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) ( C 2 8 π R ( 4 H 3 ( 1 + δ π H 2 ) β 1 ) ) 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 3 2 ) 1 .
The variation of statefinder parameter s against r in the case of SMHDE is given in Figure 12. This corresponds to the standard Λ CDM cosmological model for accelerating universe as it is clear from the plot that we have r = 1 for s = 0 .

7. Om Parameter

Regarding observations, Om diagnostic appeared to be more simple as compared to { r , s } statefinder. One among its salient features is its capability of discriminating dynamical DE models from Λ -CDM without any information about matter density. For Λ -CDM, quintessence and phantom models, Om has zero negative and positive curvatures, respectively.
Figure 13 gives the plot of O m -diagnostic for THDE versus redshift parameter z. Within the range 0.9 z 0.1 , the trajectory shows positive curvature. This leads to phantom behavior of THDE model in the fractal universe. The plot of O m -diagnostic for RHDE against redshift parameter z is given in Figure 14. For range 0.9 z 0.1 , the plot represents positive curvature leading to phantom behavior of the RHDE model in the fractal universe. Following the same pattern, SMHDE also crosses the phantom line. It is evident in Figure 15, which shows positive curvature of trajectory.

8. ω e f f - ω e f f Plane Analysis

The dynamical property of DE models can be analyzed by using ω e f f - ω e f f plane proposed by Linder and Caldwell [52]. Here, ω e f f = p e f f ρ e f f is the EoS parameter and ω e f f is the evolution parameter for EoS parameter. To get this evolution parameter, we differentiate EoS parameter ω e f f with respect to l n a . The ω e f f - ω e f f plane has been divided into two parts, namely “thawing” and “freezing”. The conditions for thawing and freezing parts are ω e f f < 0 , ω e f f > 0 and ω e f f < 0 , ω e f f < 0 , respectively.

8.1. THDE

By differentiating ω e f f for THDE with respect to l n a , we get the following expression for ω e f f in the case of THDE:
ω e f f = H 1 ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + B H 4 2 δ 2 ( ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + B H 4 2 δ ( 3 ξ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H γ 3 B ( 4 2 δ ) H 3 2 δ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ + B γ 3 ( 4 2 δ ) ( H 2 2 δ ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) 3 5 H 2 ω γ 3 ν 2 3 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ + ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H 2 3 + ( B ( 4 2 δ ) ( 3 2 δ ) H 2 2 δ 3 2 ( 1 γ ω 6 ν 2 γ 2 ) ) ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ 2 ) 2 H 1 γ ω 6 ν 2 γ 2 B 3 ( 4 2 δ ) H 3 2 δ 1 + ( 2 2 δ ) H 1 2 δ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ 2 ) ) ( 3 ξ ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) γ 3 B H 4 2 δ + B ( 4 2 δ ) H 2 2 δ γ 3 ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ ) ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H + B ( 4 2 δ ) H 3 2 δ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ( 1 + z ) 2 γ ( γ 2 ) 6 H ( 1 γ ω 6 ( 1 + z ) 2 γ γ 2 ) B ( 4 2 δ ) H 3 2 δ ) .
In the case of THDE, the plot for ω e f f - ω e f f is given by Figure 16. We use the same values of constant as described above. It is quite clear from the plot that the Eos parameter ω e f f < 0 and the evolutionary parameter ω e f f < 0 as well. This shows THDE represents the freezing region.

8.2. RHDE

The expression for ω e f f in the case of RHDE takes the following form:
ω e f f = H 1 ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + F H 2 ( H 2 + δ π ) 1 2 ( ( ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + F H 2 ( H 2 + δ π ) 1 ) ( 3 ξ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H ( γ 3 ) + H 1 ( γ 3 ) 2 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ) 1 ( ρ 0 3 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) 2 ω 3 H 4 γ 4 ν 2 + 2 ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 ( F H 2 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 4 3 F H 4 ( H 2 + δ π ) 3 ( H 2 + 2 δ π ) + 2 3 F H 4 ( H 2 + δ π ) 2 ( 1 γ ω 6 ν 2 γ 2 ) ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 3 H 2 ω γ 3 ν 2 H 2 ω γ 3 ν 2 ρ 0 3 ( 3 ξ + γ 3 ) 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) + 2 F ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ) 2 ( 3 H 2 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 4 H 4 ( H 2 + δ π ) 3 ( H 2 + 2 δ π ) + 2 H 4 ( H 2 + δ π ) 2 ) 2 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 2 1 H 2 ( γ 3 ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 4 π H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) + 1 ) + 1 H ( γ 3 ) ( 3 ξ H ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) + H ( 3 γ ) F H 4 ( H 2 + δ π ) 1 + 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ν 2 ) 2 H ( 1 γ ω 6 ν 2 γ 2 ) 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) 1 ) ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω γ 3 ν 2 ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) 2 3 F H 3 ( H 2 + δ π ) 2 ( H 2 + 2 δ π ) ) 1 ) .
The plot for ω e f f - ω e f f in the case of RHDE is given by Figure 17. Clearly, it is representing freezing region, being Eos parameter ω e f f < 0 and the evolutionary parameter ω e f f < 0 .

8.3. SMHDE

For SMHDE, we get ω e f f as given below:
ω e f f = H 1 3 C 2 H 4 8 π R ( 1 + δ π H 2 ) β 1 + ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) 2 ( 3 C 2 H 4 8 π R ( 1 + δ π H 2 ) β 1 + ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) ( 1 3 γ ( 3 ξ ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H 1 3 C 2 8 π R ( 12 H 2 1 + δ π H 2 β 1 10 δ π β 1 + δ π H 2 β 1 + 4 δ 2 π 2 β ( β 1 ) 1 + δ π H 2 β 1 H 2 ) ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( ( 1 γ ω 6 ν 2 γ 2 ) 2 H C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) 3 + 3 C 2 8 π R 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ( 2 H ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ) 1 ( 2 3 ω γ 4 ν 2 H 4 + 1 3 ( ρ 0 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) ( ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) + ρ 0 ( 3 ξ + γ 3 ) 2 H 2 ( 1 + z ) ( 3 ξ + γ 3 ) ) 2 ( 1 γ ω 6 γ 2 ν 2 ) C 2 8 π R ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ( 12 H 2 1 + δ π H 2 β 1 10 δ π β 1 + δ π H 2 β 1 + 4 δ 2 π 2 β ( β 1 ) 1 + δ π H 2 β 1 H 2 ) ( ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R ( 4 H 3 ( 1 + δ π H 2 ) β 1 ) 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) 2 5 3 ω γ 3 ν 2 H 2 ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) + ( H 2 3 + γ ) 3 C 2 8 π R ( 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) ) 1 H ( 3 γ ) ( 3 ξ H ρ 0 ( 1 + z ) ( 3 ξ + γ 3 ) 3 C 2 H 5 ( 3 γ ) 8 π R ( 1 + δ π H 2 ) β 1 3 C 2 8 π R ( 4 H 3 1 + δ π H 2 β 1 ) 2 δ π β H 1 + δ π H 2 β ) ( ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) ( C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 ) 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 3 C 2 8 π R ( 4 H 3 1 + δ π H 2 β 1 2 δ β H 1 + δ π H 2 β 1 ) ρ 0 3 ( 3 ξ + γ 3 ) ( 1 + z ) ( 3 ξ + γ 3 ) H H 3 ( ω ν 2 γ 3 3 ) ( 2 H ( 1 γ ω 6 ν 2 γ 2 ) C 2 8 π R 4 H 3 ( 1 + δ π H 2 ) β 1 2 δ π β H ( 1 + δ π H 2 ) β 1 ) 1 ) .
In the case of SMHDE, the plot for ω e f f - ω e f f is given by Figure 18. Clearly, it is representing thawing region, being the EoS parameter ω e f f < 0 and the evolutionary parameter ω e f f > 0 .

9. Conclusions

In this paper, we have considered the flat FRW universe in fractal spacetime. We have investigated some cosmological implications such as EoS parameter, deceleration parameter, Om parameter, State-finder parameter and ω e f f - ω e f f plane. For this purpose we have taken three DE models, namely Tsallis DE, Renyi DE and Sharma mittal DE. We have selected an appropriate interacting situation for HDE and DM to avoid coincidence. We have derived the expressions for all above-mentioned cosmological parameters in terms of red shift parameter z. We have fixed the constants as ξ = 2 , ρ = 1 , δ = 1 , B = 1 , γ = 2 and α = 1 generally. For some cases of RHDE and SMHDE, we have taken ξ = 2 and β = 800 . The deceleration parameter q describes the evolution of the universe. The three cosmological models, i.e., THDE, RHDE and SMHDE, of the evolving universe transit from early decelerating phase ( q > 0 ) to current accelerating phase ( q < 0 ). On the basis of the EoS parameter ω e f f , we can say THDE exhibit quintom-like nature as it shows transition of the universe from Cosmological constant towards quintessence. RHDE shows conversion of the universe from Cosmological constant towards quintessence. SMHDE shows quintom-like behavior. For the stability analysis of THDE, RHDE and SMHDE, we investigate the graphical behavior of square speed of sound. All three interacting HDE models under consideration in fractal universe are stable. From ω e f f ω e f f plane, we can say that THDE and RHDE represent the freezing region and reinforce the accelerated expansion of the universe, whereas the third model, SMHDE, represents the thawing region. In the same manner, we have { r , s } = { 1 , 0 } representing standard Λ C D M cosmological model of the accelerating universe for THDE, RHDE and SMHDE. The plots for Om diagnostic of THDE, RHDE and SMHDE are almost identical. The plot shows positive curvature leading to phantom behavior of selected HDE models. It is noted that the constraints on effective EoS are also compatible with recent results [53].
Finally, it is observed that EoS parameter of all DE models lies within the following constraints at the present epoch (Table 1).

Author Contributions

M.U.S. has done analysis as well as proof reading of the manuscript, A.I. has completed the manuscript with calculations and figures while A.J. has contributed in the improvement of the manuscript and prepared the final draft.

Funding

This research received no external funding.

Acknowledgments

Abdul Jawad is thankful to HEC, Islamabad, Pakistan for its financial support under the grant No: 5412/Federal/NRPU/R&D/HEC/2016.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peebles, P.J.E. Principles of Physical Cosmology; Princeton University Press: Princeton, NJ, USA, 1993. [Google Scholar]
  2. Davis, M. Is the Universe Homogeneous on Large Scales? In Critical Dialogues in Cosmology; Turok, N., Ed.; World Scientific Publishing Company: Singapore, 1997; p. 13. [Google Scholar]
  3. Wu, K.K.S.; Lahav, O.; Rees, M. The Large-Scale Smoothness of the Universe. Nature 1999, 397, 225. [Google Scholar] [CrossRef]
  4. Joyce, M.; Montuori, M.; Sylos, L.F.; Pietronero, L. Comment on the paper “the eso slice project galaxy redshift survey: v. evidence for a d=3 sample dimensionality”. Astron. Astrophys. 1999, 344, 387. [Google Scholar]
  5. Coles, P. Cosmology: An unprincipled Universe? Nature 1998, 391, 120. [Google Scholar] [CrossRef]
  6. Smoot, G.; Bennett, C.L.; Kogut, A.; Wright, E.L.; Aymon, J.; Boggess, N.W.; Cheng, E.S.; de Amici, G.; Gulkis, S.; Hauser, M.G.; et al. Structure in the COBE differential microwave radiometer first-year maps. Astrophys. J. 1992, 396, L1. [Google Scholar] [CrossRef]
  7. Joyce, M.; Anderson, P.W.; Montuori, M.; Sylos Labini, F.; Pietronero, L. Fractal Cosmology in an Open Universe. Europhys. Lett. 2000, 49, 416. [Google Scholar] [CrossRef]
  8. Mandelbrot, B. The Fractal Geometry of Nature. Am. J. Phys. 1983, 51, 286. [Google Scholar] [CrossRef]
  9. Mandelbrot, B.B. Fractality, Lacunarity, and the Near-Isotropic Distribution of Galaxies. In Current Topics in Astrofundamental Physics: Primordial Cosmology; Sanchez, N., Zichichi, A., Eds.; NATO ASI Series (Series C: Mathematical and Physical Sciences); Springer: Dordrecht, The Netherlands, 1998; Volume 511, p. 583. [Google Scholar]
  10. Coleman, P.H.; Pietronero, L. The fractal structure of the universe. Phys. Rep. 1992, 213, 311. [Google Scholar] [CrossRef]
  11. Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.; Goobar, A.; Groom, D.E.; et al. Measurements of Ω and Λ from 42 high-redshift supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  12. Abazajian, K.; Adelman-McCarthy, J.K.; Agüeros, M.A.; Allam, S.S.; Anderson, K.S.; Anderson, S.F.; Annis, J.; Bahcall, N.A.; Baldry, I.K.; Bastian, S.; et al. The second data release of the sloan digital sky survey. Astron. J. 2004, 128, 502. [Google Scholar] [CrossRef]
  13. Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.; Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. Astron. J. 1998, 116, 1009. [Google Scholar] [CrossRef]
  14. Bamba, K.; Capozziello, S.; Nojiri, S.; Odintsov, S.D. Dark energy cosmology: the equivalent description via different theoretical models and cosmography tests. Astro. Phys. Space Sci. 2012, 342, 155–228. [Google Scholar] [CrossRef] [Green Version]
  15. Hooft, G. Dimensional Reduction in Quantum Gravity. arXiv 2009, arXiv:gr-qc/9310026v2. [Google Scholar]
  16. Susskind, L. The World as a Hologram. J. Math. Phys. 1995, 36, 6377. [Google Scholar] [CrossRef]
  17. Cohen, A.; Kaplan, D.; Nelson, A. Effective Field Theory, Black Holes, and the Cosmological Constant. Phys. Rev. Lett. 1999, 82, 4971. [Google Scholar] [CrossRef]
  18. Li, M. A Model of Holographic Dark Energy. Phys. Lett. B 2004, 603, 1. [Google Scholar] [CrossRef]
  19. Fischler, W.; Susskind, L. Holography and Cosmology. arXiv 1998, arXiv:hep-th/9806039. [Google Scholar]
  20. Nojiri, S.; Odintsov, S.D. Unifying phantom inflation with late-time acceleration: Scalar phantom non-phantom transition model and generalized holographic dark energy. Gen. Relativ. Gravit. 2006, 38, 1285–1304. [Google Scholar] [CrossRef]
  21. Nojiri, S.; Odintsov, S.D. Inhomogeneous equation of state of the universe: Phantom era, future singularity, and crossing the phantom barrier. Phys. Rev. D 2005, 72, 023003. [Google Scholar] [CrossRef] [Green Version]
  22. Sheykhi, A. Holographic scalar field models of dark energy. Phys. Rev. D 2011, 84, 107302. [Google Scholar] [CrossRef]
  23. Hu, B.; Ling, Y. Interacting dark energy, holographic principle, and coincidence problem. Phys. Rev. D 2006, 73, 123510. [Google Scholar] [CrossRef] [Green Version]
  24. Ma, Y.Z.; Gong, Y.; Chen, X. Features of holographic dark energy under combined cosmological constraints. Eur. Phys. J. C 2009, 60, 303. [Google Scholar] [CrossRef]
  25. Masi, M. A step beyond Tsallis and Renyi entropies. Phys. Lett. A 2005, 338, 217. [Google Scholar] [CrossRef]
  26. Tsallis, C. The Nonadditive Entropy Sq and Its Applications in Physics and Elsewhere: Some Remarks. Entropy 2011, 13, 1765. [Google Scholar] [CrossRef]
  27. Renyi, A. Probability Theory; North-Holland: Amsterdam, The Netherlands, 1970. [Google Scholar]
  28. Tsallis, C. Possible generalization of Boltzmann-Gibbs statistics. J. Stat. Phys. 1988, 52, 479. [Google Scholar] [CrossRef]
  29. Jahromi, A.S.; Moosavi, S.A.; Moradpour, H.; Graça, J.M.; Lobo, I.P.; Salako, I.G.; Jawad, A. Generalized entropy formalism and a new holographic dark energy model. J. Phys. Lett. B. 2018, 2, 52. [Google Scholar]
  30. Nojiri, S.; Odintsov, S.D.; Saridakis, E.N. Modified cosmology from extended entropy with varying exponent. Eur. Phys. J. C 2019, 79, 242. [Google Scholar] [CrossRef]
  31. Calcagni, G. Quantum field theory, gravity and cosmology in a fractal universe. JHEP 2010, 1003, 120. [Google Scholar] [CrossRef]
  32. Calcagni, G. Fractal universe and quantum gravity. Phys. Rev. Lett. 2010, 104, 251301. [Google Scholar] [CrossRef]
  33. Calcagni, G. Multi-scale gravity and cosmology. JCAP 2013, 12, 41. [Google Scholar] [CrossRef]
  34. Salvatelli, V.; Said, N.; Bruni, M.; Melchiorri, A.; Wands, D. Indications of a Late-Time Interaction in the Dark Sector. Phys. Rev. Lett. 2014, 113, 181301. [Google Scholar] [CrossRef] [Green Version]
  35. Yang, T.; Guo, Z.K. Reconstructing the interaction between dark energy and dark matter using Gaussian processes. Phys. Rev. D 2015, 91, 123533. [Google Scholar] [CrossRef] [Green Version]
  36. Feng, C.; Wang, B.; Su, R.-K. Testing the viability of the interacting holographic dark energy model by using combined observational constraints. JCAP 2007, 9, 005. [Google Scholar] [CrossRef]
  37. Zimdahl, W.; Pavon, D.; Chimento, L.P. Interacting Quintessence. Phys. Lett. B 2001, 521, 133. [Google Scholar] [CrossRef]
  38. Wang, B.; Zang, J.; Lin, C.-Y.; Abdalla, E.; Micheletti, S. Interacting dark energy and dark matter: Observational constraints from cosmological parameters. Nucl. Phys. B 2007, 778, 69–84. [Google Scholar] [CrossRef]
  39. Amendola, L.; Tsujikawa, S.; Sami, M. Phantom damping of matter perturbations. Phys. Lett. B 2006, 632, 155–158. [Google Scholar] [CrossRef] [Green Version]
  40. Ferreira, E.G.M.; Quintin, J.; Costa, A.A.; Abdalla, E.; Wang, B. Evidence for interacting dark energy from BOSS. arXiv 2017, arXiv:1412.2777. [Google Scholar] [CrossRef]
  41. Gupta, G.; Saridakis, E.N.; Sen, A.A. Nonminimal quintessence and phantom with nearly flat potentials. Phys. Rev. D 2009, 79, 123013. [Google Scholar] [CrossRef]
  42. Santos, L. Constraining interacting dark energy with CMB and BAO future surveys. Phys. Rev. D 2017, 96, 103529. [Google Scholar] [CrossRef] [Green Version]
  43. Tsallis, C.; Cirto, L.J.L. Black hole thermodynamical entropy. Eur. Phys. J. C 2013, 73, 2487. [Google Scholar] [CrossRef]
  44. Tavayef, M.; Sheykhi, A.; Bamba, K.; Moradpour, H. Tsallis Holographic Dark Energy. arXiv 2018, arXiv:1804.02983v1. [Google Scholar] [CrossRef]
  45. Majhi, A. Non-extensive Statistical Mechanics and Black Hole Entropy From Quantum Geometry. Phys. Lett. B 2017, 775, 32–36. [Google Scholar] [CrossRef]
  46. Komatsu, N. Cosmological model from the holographic equipartition law with a modified Rényi entropy. Eur. Phys. J. C 2017, 77, 229. [Google Scholar] [CrossRef]
  47. Moradpour, H.; Bonilla, A.; Abreu, E.M.C.; Neto, J.A. Accelerated cosmos in a nonextensive setup. Phys. Rev. D 2017, 96, 123504. [Google Scholar] [CrossRef] [Green Version]
  48. Czinner, V.G.; Iguchi, H. Black hole entropy and the zeroth law of thermodynamics. Int. J. Mod. Phys. D 2015, 24, 1542015. [Google Scholar] [CrossRef]
  49. Sharma, B.D.; Mitta, P. New nonadditive measures of entropy for discrete probability distributions. J. Math. Sci. 1975, 10, 28–40. [Google Scholar]
  50. Sharma, B.D.; Mitta, P. New non-additive measures of relative information. J. Cmbin Inform. Syst. Sci. 1977, 2, 122–132. [Google Scholar]
  51. Sahni, V. Statefinder—A new geometrical diagnostic of dark energy. JETP Lett. 2003, 77, 201–206. [Google Scholar] [CrossRef]
  52. Caldwell, R.R.; Linder, E.V. The Limits of Quintessence. Phys. Rev. Lett. 2005, 95, 141301. [Google Scholar] [CrossRef]
  53. Aviles, A.; Gruber, C.; Luongo, O.; Quevedo, H. Cosmography and constraints on the equation of state of the Universe in various parametrizations. Phys. Rev. D 2012, 86, 123516. [Google Scholar] [CrossRef]
  54. Ade, P.A.R. Planck 2015 results. XIII. Cosmological parameters. A&A 2016, 594, A13. [Google Scholar]
  55. Hinshaw, G.; Larson, D.; Komatsu, E.; Spergel, D.N.; Bennett, C.L.; Dunkley, J.; Nolta, M.R.; Halpern, M.; Hill, R.S.; Odegard, N.; et al. Nine-year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results. Astrophys. J. Suppl. 2013, 208, 19. [Google Scholar] [CrossRef]
Figure 1. Deceleration parameter.
Figure 1. Deceleration parameter.
Symmetry 11 01174 g001
Figure 2. Deceleration parameter.
Figure 2. Deceleration parameter.
Symmetry 11 01174 g002
Figure 3. Deceleration parameter.
Figure 3. Deceleration parameter.
Symmetry 11 01174 g003
Figure 4. Eos parameter.
Figure 4. Eos parameter.
Symmetry 11 01174 g004
Figure 5. EoS parameter.
Figure 5. EoS parameter.
Symmetry 11 01174 g005
Figure 6. EoS parameter.
Figure 6. EoS parameter.
Symmetry 11 01174 g006
Figure 7. Square speed of sound.
Figure 7. Square speed of sound.
Symmetry 11 01174 g007
Figure 8. Square speed of sound.
Figure 8. Square speed of sound.
Symmetry 11 01174 g008
Figure 9. Square speed of sound.
Figure 9. Square speed of sound.
Symmetry 11 01174 g009
Figure 10. Statefinder parameter.
Figure 10. Statefinder parameter.
Symmetry 11 01174 g010
Figure 11. Statefinder parameter.
Figure 11. Statefinder parameter.
Symmetry 11 01174 g011
Figure 12. Statefinder parameter.
Figure 12. Statefinder parameter.
Symmetry 11 01174 g012
Figure 13. Om diagnostic parameter.
Figure 13. Om diagnostic parameter.
Symmetry 11 01174 g013
Figure 14. Om diagnostic parameter.
Figure 14. Om diagnostic parameter.
Symmetry 11 01174 g014
Figure 15. Om diagnostic parameter.
Figure 15. Om diagnostic parameter.
Symmetry 11 01174 g015
Figure 16. Plot of EoS parameter ω e f f = X and its evolution parameter ω e f f = G .
Figure 16. Plot of EoS parameter ω e f f = X and its evolution parameter ω e f f = G .
Symmetry 11 01174 g016
Figure 17. Plot of EoS parameter ω e f f = X and its evolution parameter ω e f f = G .
Figure 17. Plot of EoS parameter ω e f f = X and its evolution parameter ω e f f = G .
Symmetry 11 01174 g017
Figure 18. Plot of EoS parameter ω e f f = X and its evolution parameter ω e f f = G .
Figure 18. Plot of EoS parameter ω e f f = X and its evolution parameter ω e f f = G .
Symmetry 11 01174 g018
Table 1. Summary of the observational data on ω D E .
Table 1. Summary of the observational data on ω D E .
ω e f f Observational SchemesRefrences
1.1313 0.25 + 0.24 Planck+WP+BAO[54]
1.09 ± 0.17 Planck+WP+Union 2.1[54]
1.13 0.14 + 0.13 Planck+WP+SNLS[54]
1.24 0.19 + 0.18 WMAP+eCMB+BAO+ H 0 +SNe Ia[54]
1.073 0.089 + 0.090 WMAP+eCMB+BAO+ H 0 [55]
1.084 ± 0.063 WMAP+eCMB+BAO+ H 0 +SNe[55]

Share and Cite

MDPI and ACS Style

Shahzad, M.U.; Iqbal, A.; Jawad, A. Dynamical Properties of Dark Energy Models in Fractal Universe. Symmetry 2019, 11, 1174. https://doi.org/10.3390/sym11091174

AMA Style

Shahzad MU, Iqbal A, Jawad A. Dynamical Properties of Dark Energy Models in Fractal Universe. Symmetry. 2019; 11(9):1174. https://doi.org/10.3390/sym11091174

Chicago/Turabian Style

Shahzad, Muhammad Umair, Ayesha Iqbal, and Abdul Jawad. 2019. "Dynamical Properties of Dark Energy Models in Fractal Universe" Symmetry 11, no. 9: 1174. https://doi.org/10.3390/sym11091174

APA Style

Shahzad, M. U., Iqbal, A., & Jawad, A. (2019). Dynamical Properties of Dark Energy Models in Fractal Universe. Symmetry, 11(9), 1174. https://doi.org/10.3390/sym11091174

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop