Cosmological Consequences of New Dark Energy Models in Einstein-Aether Gravity

Abstract

In this paper, we reconstruct various solutions for the accelerated universe in the Einstein-Aether theory of gravity. For this purpose, we obtain the effective density and pressure for Einstein-Aether theory. We reconstruct the Einstein-Aether models by comparing its energy density with various newly proposed holographic dark energy models such as Tsallis, Rényi and Sharma-Mittal. For this reconstruction, we use two forms of the scale factor, power-law and exponential forms. The cosmological analysis of the underlying scenario has been done by exploring different cosmological parameters. This includes equation of state parameter, squared speed of sound and evolutionary equation of state parameter via graphical representation. We obtain some favorable results for some values of model parameters

1. Introduction

Nowadays, it is believed that our universe is undergoing an accelerated expansion with the passage of cosmic time. This cosmic expansion has been confirmed through various observational schemes such as supernova type Ia (SNIa) [1,2,3,4] and the cosmic microwave background (CMB) [5,6,7,8,9]. The source behind the expansion of the universe is a mysterious force called dark energy (DE) and its nature is still ambiguous [10,11,12,13]. The current Planck data shows that DE accounts for $68.3\%$ of the total energy contents of the universe. The first candidate for describing the DE phenomenon is the cosmological constant but it has fine tuning and cosmic coincidence problems. Due to this reason, different DE models as well as theories of gravity with modifications have been suggested. The dynamical DE models include a family of Chaplygin gas as well as holographic DE models, scalar field models such as K-essence, phantom, quintessence, ghost, etc. [14,15,16,17,18,19,20,21,22,23,24,25,26].

One of the DE model is the holographic DE (HDE) model which becomes a favorable technique now-a-days to study the DE mystery. This model is established in the framework of holographic principle which corresponds to the area instead of volume for the scaling of the number of degrees of freedom of a system. This model is an interesting effort in exploring the nature of DE in the framework of quantum gravity. In addition, the HDE model gives the relationship between the energy density of quantum fields in vacuum (as the DE candidate) to the cutoffs (infrared and ultraviolet). Cohen et al. [27] provided a very useful result about the expression of the HDE model density which is based on the vacuum energy of the system. The black hole mass should not be overcome by the maximum amount of the vacuum energy. Taking into account the nature of spacetime along with long term gravity, various entropy formalisms have been used to discuss the gravitational and cosmological setups [28,29,30,31,32,33]. Recently, some new HDE models are proposed, like Tsallis HDE (THDE) [31], Rényi HDE model (RHDE) [32] and Sharma-Mittal HDE (SMHDE) [33].

The examples of theories with modification setups include $f\left(R\right),f\left(T\right)$, $f(R,\mathcal{T}),f\left(G\right)$ etc., where R shows the Ricci scalar representing the curvature, T means the torsion scalar, $\mathcal{T}$ is the trace of the energy-momentum tensor and G goes as the invariant of Gauss–Bonnet [34,35,36,37,38,39,40,41,42,43,44,45]. For recent reviews in terms of DE problems including modified gravity theories, see, for instance [46,47,48,49,50,51,52]. The Einstein-Aether theory is one of the modified theories of gravity [53,54] and accelerated expansion phenomenon of the universe has also been investigated in this theory [55]. Meng et al. have also discussed the current cosmic acceleration through DE models in this gravity [56,57]. Recently, Pasqua et al. [58] have made versatile studies on cosmic acceleration through various cosmological models in the presence of HDE models.

In the present work, we will develop the Einstein-Aether gravity models in the presence of modified HDE models and well-known scale factors. For these models of modified gravity, we will extract various cosmological parameters. In the next section, we will give a brief review of the Einstein-Aether theory. In Section 3, we present the basic cosmological parameters as well as well-known scale factors. We will discuss the cosmological parameters for modfied HDE models in Section 4, Section 5 and Section 6. In the last section, we will summarize our results.

2. Einstein-Aether Theory

As our universe is full with many of the natural occurring phenomenons. One of them is transfer of light from one place to another and second is how gravity acts. To explain these kinds of phenomenons, many of the physicists were used the concept of Aether in many of the theories. In modern physics, Aether indicates a physical medium that is spread homogeneously at each point of the universe. Hence, it was considered that it is a medium in space that helps light to travel in a vacuum. According to this concept, a particular static frame reference is provided by Aether and everything has absolute relative velocity in this frame. That is suitable for Newtonian dynamics extremely well. But, when Einstein performed different experiments on optics in his theory of relativity, then Einstein rejected this ambiguity. When CMB was introduced, many of the people took it a modern form of Aether. Gasperini has popularized Einstein-Aether theories [59]. This theory is said to be covariant modification of general relativity in which unit time like vector field(aether) breaks the Lorentz Invariance (LI) to examine the gravitational and cosmological effects of dynamical preferred frame [53]. Following is the action of Einstein-Aether theory [60,61].

$$S=\int d^{4}x\sqrt{-g}\left(\frac{R}{4\pi G}+L_{EA}+L_m\right),$$

(1)

where $L_{EA}$ represents the Lagrangian density for the vector field and $L_m$ indicates Lagrangian density of matter field. Further, $g,R$ and G indicate determinant of the metric tensor $g^{\mu \nu }$, Ricci scalar and gravitational constant respectively. The Lagrangian density for vector field can be written as

$$\begin{array}{ccc}\hfill L_{EA}& =& \frac{M^2}{16\pi G}F\left(K\right)+\frac{1}{16\pi G}\lambda (A^a{A}_a+1),\hfill \end{array}$$

(2)

$$\begin{array}{ccc}\hfill K& =& M^{-2}{K}_{cd}^{ab}{\nabla }_a{A}^c{\nabla }_b{A}^{d}c\hfill \end{array}$$

(3)

$$\begin{array}{ccc}\hfill K_{cd}^{ab}& =& c_1{g}^{ab}{g}_{cd}+c_2{\delta }_c^a{\delta }_d^b+c_3{\delta }_d^a{\delta }_c^b,a,b=0,1,2,3.\hfill \end{array}$$

(4)

where $\lambda$ represents a Lagrangian multiplier, dimensionless constants are denoted by $c_i,M$ referred as coupling constant parameter and $A^a$ is a tensor of rank one, that is a vector. The function $F\left(K\right)$ is any arbitrary function of K. We obtain the Einstein field equations from Equation (1) for the Einstein-Aether theory as follows

$$\begin{array}{ccc}\hfill G_{ab}& =& T_{ab}^{EA}+8\pi G{T}_{ab}^m,\hfill \end{array}$$

(5)

$$\begin{array}{ccc}\hfill {\nabla }_a\left(\frac{dF}{dK}{J}_b^a\right)& =& 2\lambda A_b,\hfill \end{array}$$

(6)

where $J_b^a=-2K_{bc}^{ad}{\nabla }_d{A}^c,T_{ab}^{EA}$ shows energy momentum-tensor for vector field and $T_{ab}^m$ indicates energy-momentum tensor for mater field. These tensors are given as

$$\begin{array}{ccc}\hfill T_{ab}^m& =& (\rho +p)u_{a}ub+p{g}_{ab},\hfill \end{array}$$

(7)

$$\begin{array}{ccc}\hfill T_{ab}^{EA}& =& {\displaystyle \frac{1}{2}}{\nabla }_d\left((J_a^d{A}_b-J_a^d{A}_b-J_{\left(ab\right)}{A}^d){\displaystyle \frac{dF}{dK}}\right)-Y_{\left(ab\right)}{\displaystyle \frac{dF}{dK}}\hfill \\ & +& {\displaystyle \frac{1}{2}}{g}_{ab}{M}^{2}F+\lambda A_a{A}_b,\hfill \end{array}$$

(8)

where p and $\rho$ represent energy density and pressure of the matter respectively. Furthermore, $u_a$ expresses the four-velocity vector of the fluid and given as $u_a=(1,0,0,0)$ and $A_a$ is time-like unitary vector and is defined as $A_a=(1,0,0,0)$. Moreover $Y_{ab}$ is defined as

$$Y_{ab}=c_1\left(\left({\nabla }_d{A}_a\right)\left({\nabla }^d{A}_b\right)-\left({\nabla }_a{A}_d\right)\left({\nabla }_a{A}^d\right)\right),$$

(9)

where indices $\left(ab\right)$ show the symmetry.

The Friedmann equations modified by the Einstein-Aether gravity are given as follows

$$\begin{array}{ccc}\hfill ϵ\left(\frac{F}{2K}-\frac{dF}{dK}\right)H^2+\left(H^2+\frac{k}{a^2}\right)& =& \left(\frac{8\pi G}{3}\right)\rho ,\hfill \end{array}$$

(10)

$$\begin{array}{ccc}\hfill ϵ\frac{d}{dt}\left(H\frac{dF}{dK}\right)+\left(-2\dot{H}+\frac{2k}{a^2}\right)& =& 8\pi G(p+\rho ).\hfill \end{array}$$

(11)

Here K becomes $K=\frac{3ϵH^2}{M^2}$, where $ϵ$ is a constant parameter. The energy density of Einstein-Aether theory is denoted by ${\rho }_{EA}$ and called the effective energy density, while the effective pressure in the Einstein-Aether gravity is given by $p_{EA}$. So, we can rewrite Equations (10) and (11) as

$$\begin{array}{ccc}\hfill \left(H^2+\frac{k}{a^2}\right)& =& \left(\frac{8\pi G}{3}\right)\rho +\frac{1}{3}{\rho }_{EA},\hfill \end{array}$$

(12)

$$\begin{array}{ccc}\hfill \left(-2\dot{H}+\frac{2k}{a^2}\right)& =& 8\pi G(p+\rho )+({\rho }_{EA}+p_{EA}),\hfill \end{array}$$

(13)

where

$$\begin{array}{ccc}\hfill {\rho }_{EA}& =& 3ϵH^2\left(\frac{dF}{dK}-\frac{F}{2K}\right),\hfill \end{array}$$

(14)

$$\begin{array}{ccc}\hfill p_{EA}& =& -3ϵH^2\left(\frac{dF}{dK}-\frac{F}{2K}\right)-ϵ\left(\dot{H}\frac{dF}{dK}+H\frac{d\dot{F}}{dK}\right).\hfill \end{array}$$

(15)

$$\begin{array}{ccc}& =& {\rho }_{EA}-\frac{{\dot{\rho }}_{EA}}{3H}.\hfill \end{array}$$

(16)

The equation of state (EoS) parameter for the Einstein-Aether can be obtained by using Equations (14) and (15), and it is given by

$${\omega }_{EA}=\frac{p_{EA}}{{\rho }_{EA}}=-1-\frac{\dot{H}\frac{dF}{dK}+H\frac{d\dot{F}}{dK}}{3H^2(\frac{dF}{dK}-\frac{F}{2K})}.$$

(17)

3. Cosmological Parameters

To understand the geometry of the universe, the following are some basic cosmological parameters.

3.1. Equation of State Parameter

In order to categorize the different phases of the evolving universe, the EoS parameter is widely used. In particular, the decelerated and accelerated phases contain DE, DM, radiation dominated eras. This parameter is defined in terms of energy density $\rho$ and pressure p as $\omega =\frac{p}{\rho }$.

• In the decelerated phase, the radiation era $0<\omega <\frac{1}{3}$ and cold DM era $\omega =0$ are included.
• The accelerated phase of the universe has following eras: $\omega =-1\Rightarrow$ cosmological constant, $-1<\omega <\frac{-1}{3}\Rightarrow$ quintessence and $\omega <-1\Rightarrow$ phantom era of the universe.

3.2. Squared Speed of Sound

To examine the behavior of DE models, there is another parameter which is known as squared speed of sound. It is denoted by $v_s^2$ and is calculated by the following formula

$$v_s^2=\frac{\dot{p}}{\dot{\rho }}.$$

(18)

The stability of the model can be checked by this relation. If its graph is showing negative values then we may say that model is unstable and in case of non-negative values of the graph, it represents the stable behavior of the model.

3.3. $\omega$-${\omega }^{\prime }$ Plane

There are different DE models which have different properties. To examine their dynamical behavior, we use $\omega$-${\omega }^{\prime }$ plane, where prime denotes the derivative with respect to $lna$ and subscript $\Lambda$ indicates DE scenario. This method was developed by Caldwell and Linder [62] and divides $\omega$-${\omega }^{\prime }$ plane into two parts. One is the freezing part in which evolutionary parameter gives negative behavior for negative EoS parameter, i.e., ${\omega }^{\prime }<0,\omega <0$, while for positive behavior of evolutionary parameter corresponding to negative EoS parameter yields the thawing part (${\omega }^{\prime }>0,\omega <0$) of the evolving universe.

3.4. Scale Factor

The scale factor is the measure of how much the universe has expanded since a given time. It is represented by $a\left(t\right)$. Since the latest cosmic observations have shown that the universe is accelerating so $a\left(t\right)>0$. As Einstein-Aether is one of the modified theory which may produce the accelerated expansion of the universe, by using this theory, we can reconstruct various well-known DE models. In order to do this, we take some modified HDE models such as THDE, RHDE and SMHDE models. Since $F\left(K\right)$ is a function that the Einstein-Aether theory contains, which can be determined by comparing the densities with the above DE models. For this purpose, we use some well-known forms of the scale factor, $a\left(t\right)$. We consider two forms of scale factors $a\left(t\right)$ in terms of power and exponential terms. These are

(i)

Power-law form: $a\left(t\right)=a_0{t}^m,m>0$, where $a_0$ is a constant which indicates the value of scale factor at present-day [63,64]. From this scale factor, we get $H,\dot{H},K$ as follows

$$H=\frac{m}{t},\dot{H}=-\frac{m}{t^2},K=\frac{3ϵm^2}{M^2{t}^2}.$$

(19)

(ii)

Exponential form: $a\left(t\right)=e^{\alpha t^{\theta }}$ where $\alpha$ is a positive constant and $\theta$ lies between 0 and 1. This scale factor gives

$$H=\alpha \theta t^{\theta -1},\dot{H}=\alpha \theta (\theta -1)t^{\theta -2},K=\frac{3ϵ{\alpha }^2{\theta }^2{t}^{2(\theta -1)}}{M^2}.$$

(20)

4. Reconstruction from the Tsallis Holographic Dark Energy Model

The energy density of THDE model is given by [31]

$${\rho }_D=B{L}^{2\delta -4},$$

(21)

where B is an unknown parameter. Taking into account Hubble radius as IR cutoff L, that is $L=\frac{1}{H}$, we have

$${\rho }_D=B{H}^{-2\delta +4}.$$

(22)

In order to construct a DE model in the framework of Einstein-Aether gravity with THDE model, we compare the densities of both models, (i.e., ${\rho }_{EA}={\rho }_D$). This yields

$$\frac{dF}{dK}-\frac{F}{2K}=\frac{B}{3ϵ}{H}^{-2\delta +2},$$

(23)

which results the following form

$$F\left(K\right)=\frac{2B{M}^{-2\delta +2}{K}^{-\delta +2}}{(-2\delta +3){\left(3ϵ\right)}^{-\delta +2}}+C_1\sqrt{K}.$$

(24)

Power-law form of scale factor:

Using the expression of $F\left(K\right)$ along with Equation (19) in (14), we obtain the energy density and pressure as

$$\begin{array}{ccc}\hfill {\rho }_{EA}& =& {\displaystyle \frac{m^2\left(3^{\delta }2B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }-9{ϵ}^2{C}_1\right)}{6K^{\frac{3}{2}}{t}^2ϵ}},\hfill \\ \hfill p_{EA}& =& m\left(4B{M}^{-2\delta }{ϵ}^{\delta }\right(-3^{\delta }{K}^{1-\delta }{M}^2(4-9m-2\delta +6m\delta )-12M^2(-2\hfill \\ & +& \delta )(-1+\delta ){\left({\displaystyle \frac{m^2ϵ}{M^2{t}^2}}\right)}^{1-\delta }){(-3+2\delta )}^{-1}+({\displaystyle \frac{9(-1+6m){ϵ}^2}{K^{\frac{3}{2}}}}+{\displaystyle \frac{3\sqrt{3}}{m^4}}\hfill \\ & \times & M^4{t}^4\sqrt{{\displaystyle \frac{m^2ϵ}{M^2{t}^2}}}\left)C_1\right){\left(36t^2ϵ\right)}^{-1}.\hfill \end{array}$$

(25)

Using these expressions of energy density and pressure, we find the values of some cosmological parameters in the following. The EoS parameter takes the following form

$$\begin{array}{ccc}\hfill {\omega }_{EA}& =& K^{\frac{3}{2}}\left(4B{M}^{-2\delta }{ϵ}^{\delta }\right(-3^{\delta }{K}^{1-\delta }{M}^2(4-9m-2\delta +6m\delta )-12M^2(-2\hfill \\ & +& \delta )(-1+\delta ){\left({\displaystyle \frac{m^2ϵ}{M^2{t}^2}}\right)}^{1-\delta }){(-3+2\delta )}^1+({\displaystyle \frac{9(-1+6m){ϵ}^2}{K^{\frac{3}{2}}}}+{\displaystyle \frac{3\sqrt{3}{M}^4}{m^4}}\hfill \\ & \times & t^4\sqrt{{\displaystyle \frac{m^2ϵ}{M^2{t}^2}}}\left)C_1\right){\left(6m\left(23^{\delta }B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }-9{ϵ}^2{C}_1\right)\right)}^{-1}.\hfill \end{array}$$

(26)

The derivative of EoS parameter with respect to $lna$ is given by

$$\begin{array}{ccc}\hfill {\omega }_{EA}^{\prime }& =& K^{\frac{3}{2}}t({\displaystyle \frac{96B{m}^2{M}^{-2\delta }(1-\delta )(-2+\delta )(-1+\delta ){ϵ}^{1+\delta }{\left(\frac{m^2ϵ}{M^2{t}^2}\right)}^{-\delta }}{t^3(-3+2\delta )}}+(-{\displaystyle \frac{3\sqrt{3}}{m^2}}\hfill \\ & \times & {\displaystyle \frac{M^{2}tϵ}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}}}+{\displaystyle \frac{12\sqrt{3}{M}^4{t}^3\sqrt{\frac{m^2ϵ}{M^2{t}^2}}}{m^4}}\left)C_1\right){\left(6m^2\left(23^{\delta }B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }-9{ϵ}^2{C}_1\right)\right)}^{-1}.\hfill \end{array}$$

(27)

We plot the EoS parameter versus z using the relation $t=\frac{1}{{(1+z)}^{\frac{1}{m}}}$ taking values of constants as $B=5,M=5,\delta =1.8,ϵ=1$ and $C_1=2$. We plot ${\omega }_{EA}$ for three different values of scale factor parameter m as $m=2,3,4$ as shown in Figure 1. All the three trajectories represent the phantom behavior of the universe related to the redshift parameter. In Figure 2, we plot ${\omega }_{EA}^{\prime }-\omega EA$ plane taking same values of the parameters for $-1\le z\le 1$. For $m=2$ and 3, the evolving EoS parameter shows negative behavior with respect to negative EoS parameter which indicates the freezing region of the universe. The trajectory of ${\omega }_{EA}^{\prime }$ for $m=4$ represents the positive behavior for negative EoS parameter and expresses the evolving universe in thawing region of the universe.

Also the squared speed of sound in the underlying scenario becomes

$$\begin{array}{ccc}\hfill v_s^2& =& \left({\left({\displaystyle \frac{m^2ϵ}{M^2{t}^2}}\right)}^{-\delta }\right(8B{m}^4{ϵ}^{\delta }(-12K^{\frac{3}{2}+\delta }{m}^2{(-2+\delta )}^2(-1+\delta )ϵ+3^{\delta }{K}^{\frac{5}{2}}{M}^2\hfill \\ & \times & t^2(4-9m-2\delta +6m\delta ){\left({\displaystyle \frac{m^2ϵ}{M^2{t}^2}}\right)}^{\delta })+3K^{\delta }{M}^{2\delta }{t}^2(-3+2\delta ){\left({\displaystyle \frac{m^2ϵ}{M^2{t}^2}}\right)}^{\delta }\hfill \\ & \times & \left(6(1-6m)m^4{ϵ}^2+\sqrt{3}{K}^{\frac{3}{2}}{M}^4{t}^4\sqrt{{\displaystyle \frac{m^2ϵ}{M^2{t}^2}}}\right)C_1\left)\right)/(12m^5{t}^2(-3+2\delta )\hfill \\ & \times & \left(-23^{\delta }B{K}^{\frac{5}{2}}{M}^2{ϵ}^{\delta }+9K^{\delta }{M}^{2\delta }{ϵ}^2{C}_1\right)).\hfill \end{array}$$

(28)

Figure 3 shows the plot of $v_s^2$ versus z to check the behavior of the Einstein-Aether model for the THDE and power-law scale factor for same values of parameters. The trajectories represent the negative behavior of the model which indicated the instability of the model.

Exponential form of scale factor:

Following the same steps for the exponential form of the scale factor, we get energy density and pressure as

$$\begin{array}{ccc}\hfill {\rho }_{EA}& =& {\displaystyle \frac{t^{-2+2\theta }{\alpha }^2{\theta }^2\left(23^{\delta }B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }-9{ϵ}^2{C}_1\right)}{6K^{\frac{3}{2}}ϵ}},\hfill \\ \hfill p_{EA}& =& {\displaystyle \frac{1}{12ϵ}}{t}^{-2+\theta }\alpha \theta ({\displaystyle \frac{2t^{\theta }\alpha \theta \left(-23^{\delta }B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }+9{ϵ}^2{C}_1\right)}{K^{\frac{3}{2}}}}-{ϵ}^2(-1\hfill \\ & +& \theta \left)\right(8B{M}^{-2\delta }(-2+\delta ){ϵ}^{-2+\delta }(3^{\delta }{K}^{1-\delta }{M}^2-6M^2(-1+\delta )({\displaystyle \frac{t^{-2+2\theta }}{M^2}}\hfill \\ & \times & {\alpha }^2ϵ{\theta }^2{)}^{1-\delta }){\left(3(-3+2\delta )\right)}^{-1}+(-{\displaystyle \frac{3}{K^{\frac{3}{2}}}}+{\displaystyle \frac{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}{{\alpha }^4{ϵ}^2{\theta }^4}}\sqrt{3}{M}^4\hfill \\ & \times & t^{4-4\theta }\left)C_1\right)).\hfill \end{array}$$

(29)

Now, by using the above density and pressure, we obtain the EoS parameter and its derivative for the Einstein-Aether gravity as follows

$$\begin{array}{ccc}\hfill {\omega }_{EA}& =& \left(K^{\frac{3}{2}}{t}^{-\theta }\right(\frac{2t^{\theta }\alpha \theta \left(-23^{\delta }B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }+9{ϵ}^2{C}_1\right)}{K^{\frac{3}{2}}}-{ϵ}^2(-1+\theta )\hfill \\ & \times & \left(8B{M}^{-2\delta }(-2+\delta ){ϵ}^{-2+\delta }\right(3^{\delta }{K}^{1-\delta }{M}^2-6M^2(-1+\delta )(\frac{t^{-2+2\theta }{\alpha }^2}{M^2}\hfill \\ & \times & ϵ{\theta }^2{)}^{1-\delta }){\left(3(-3+2\delta )\right)}^{-1}+\left(-\frac{3}{K^{\frac{3}{2}}}+\frac{\sqrt{3}{M}^4{t}^{4-4\theta }\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}{{\alpha }^4{ϵ}^2{\theta }^4}\right)\hfill \\ & \times & C_1\left)\right))/\left(2\alpha \theta \left(23^{\delta }B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }-9{ϵ}^2{C}_1\right)\right),\hfill \\ \hfill {\omega }_{EA}^{\prime }& =& \frac{1}{6{\alpha }^2\theta \left(23^{\delta }B{K}^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }-9{ϵ}^2{C}_1\right)}{K}^{\frac{3}{2}}{t}^{-2(1+\theta )}(-1+\theta )(8B{K}^{-\delta }\hfill \\ & \times & M^{-2\delta }(-2+\delta ){ϵ}^{\delta }{\left(\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}\right)}^{-\delta }(-6K^{\delta }{t}^{2\theta }{\alpha }^2(-1+\delta )ϵ(2+2\delta \hfill \\ & \times & (-1+\theta )-\theta )\theta +3^{\delta }K{M}^2{t}^2{\left(\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}\right)}^{\delta }){\left(-3+2\delta \right)}^{-1}-3\hfill \\ & \times & t^{2-4\theta }\left(3t^{4\theta }{\alpha }^4{ϵ}^2{\theta }^5+\sqrt{3}{K}^{\frac{3}{2}}{M}^4{t}^4(3-4\theta )\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}\right)C_1(K^{\frac{3}{2}}{\alpha }^4\hfill \\ & \times & {\theta }^5{)}^{-1}).\hfill \end{array}$$

The squared speed of sound for the second form of the scale factor is given by

$$\begin{array}{ccc}\hfill v_s^2& =& K^{\frac{3}{2}}{t}^{-\theta }\left(8B{M}^{-2\delta }{ϵ}^{\delta }\right(-6M^2(-2+\delta )(-1+\delta )(4+2\delta (-1+\theta )-3\theta )\hfill \\ & \times & {\left(\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}\right)}^{1-\delta }-3^{\delta }{K}^{1-\delta }{M}^2((-2+\delta )(-2+\theta )+3t^{\theta }\alpha (-3+2\hfill \\ & \times & {\delta \left)\theta \right))(-3+2\delta )}^{-1}+3t^{-4\theta }(3t^{4\theta }{\alpha }^4{ϵ}^2(-2+\theta ){\theta }^4+36t^{5\theta }{\alpha }^5{ϵ}^2{\theta }^5+\sqrt{3}\hfill \\ & \times & K^{\frac{3}{2}}{M}^4{t}^4\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}(-1+2\theta )\left)C_1{\left(K^{\frac{3}{2}}{\alpha }^4{\theta }^4\right)}^{-1}\right)\left(12\alpha \theta \right(23^{\delta }B\hfill \\ & \times & K^{\frac{5}{2}-\delta }{M}^{2-2\delta }{ϵ}^{\delta }-9{ϵ}^2{C}_1\left)\right).\hfill \end{array}$$

Figure 4 represents the graph of the EoS parameter versus z for the exponential form of the scale factor taking $B=5=M,\delta =1.8,ϵ=1,C_1=-0.5,\theta =0.5$ and scale factor parameter $\alpha =2,3,4$. This parameter represents the phantom behavior of the universe for $\alpha =3$ and after a transition from quintessence to phantom era for $\alpha =2$. For $\alpha =4$, the trajectory of the EoS parameter corresponds to the $\Lambda$-CDM model ${\omega }_{EA}=-1$. In Figure 5, the graph is plotted between ${\omega }_{EA}^{\prime }$ and ${\omega }_{EA}$. The graph represents initially freezing region and then indicates the thawing region of the evolving universe. As we increase the value of $\alpha$, the trajectories indicate the thawing region only. However, the graph of $v_s^2$ versus z as shown in Figure 6 shows the unstable behavior.

5. Reconstruction from Rényi Holographic Dark Energy Model

The energy density of the RHDE model is [33]

$${\rho }_D=\frac{3C^2{L}^{-2}}{8\pi \left(1+\frac{\delta \pi }{H^2}\right)}.$$

(30)

For the Hubble horizon, it takes the form

$${\rho }_D=\frac{3C^2{H}^2}{8\pi \left(1+\frac{\delta \pi }{H^2}\right)}.$$

(31)

Now we compare the Einstein-Aether model energy density with the RHDE model density (i.e., ${\rho }_{EA}={\rho }_D$) in order to get the reconstructed equation,

$$\frac{dF}{dK}-\frac{F}{2K}=\frac{C^2}{8\pi ϵ(1+\frac{\delta \pi }{H^2})}.$$

(32)

The solution of this equation is given by

$$F\left(K\right)=\frac{C^2}{4\pi Mϵ}\left(KM-\sqrt{3K\delta ϵ\pi }arctan\left(\frac{M\sqrt{K}}{\sqrt{3ϵ\pi \delta }}\right)\right)+C_2\sqrt{K}.$$

(33)

Power-law form of the scale factor:

Inserting all the corresponding values into Equations (14) and (15), we get density and pressure of the Einstein-Aether gravity model as follows

$$\begin{array}{ccc}\hfill {\rho }_{EA}& =& \frac{1}{8\pi t^2}3C^2{m}^2\left(1-\frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right).\hfill \\ \hfill p_{EA}& =& \frac{1}{24\pi t^5}{C}^{2}m(9m{t}^3\left(-1+\frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right)-\frac{1}{M^4}ϵ(3M^3\hfill \\ & \times & t^3(M\left(-2+\frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right)+\frac{\mathrm{ArcTan}\left(\frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta ϵ\pi }}\right)}{K}\sqrt{3}\hfill \\ & \times & \sqrt{K\delta ϵ\pi }-\frac{C_2}{\sqrt{K}}){\left(ϵ\right)}^{-1}+(m^2(-3m^{2}t{\delta }^2{ϵ}^2{\pi }^2\left(m^2-t^2\delta \pi \right)-3\hfill \\ & \times & m{\delta }^{3/2}{ϵ}^2{\pi }^{\frac{3}{2}}{\left(m^2+t^2\delta \pi \right)}^2\mathrm{ArcTan}\left(\frac{m}{t\sqrt{\delta }\sqrt{\pi }}\right)+\sqrt{3}Mt\sqrt{\delta ϵ\pi }mϵ\hfill \\ & \times & \sqrt{\frac{\delta \pi }{M^2{t}^2}}{\left(m^2+t^2\delta \pi \right)}^2{C}_2\left)\right)/\left({\left(\frac{m^2ϵ}{M^2{t}^2}\right)}^{3/2}\sqrt{\delta ϵ\pi }\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}\right(m^2\hfill \\ & +& t^2\delta \pi {)}^2\left)\right)).\hfill \end{array}$$

In this case, the EoS parameter takes the form

$$\begin{array}{ccc}\hfill {\omega }_{EA}& =& (9m{t}^3\left(-1+\frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right)-\frac{1}{M^4}ϵ(3M^3{t}^3\left(M\right(-2\hfill \\ & +& \frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)})+\frac{\sqrt{3}\sqrt{K\delta ϵ\pi }\mathrm{ArcTan}\left(\frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta ϵ\pi }}\right)}{K}-\frac{C_2}{\sqrt{K}})\hfill \\ & \times & {\left(ϵ\right)}^{-1}+\left(m^2\right(-3m^{2}t{\delta }^2{ϵ}^2{\pi }^2\left(m^2-t^2\delta \pi \right)-3m{\delta }^{\frac{3}{2}}{ϵ}^2{\pi }^{\frac{3}{2}}(m^2+t^2\hfill \\ & \times & \delta \pi {)}^2\mathrm{ArcTan}\left(\frac{m}{t\sqrt{\delta }\sqrt{\pi }}\right)+\sqrt{3}Mt\sqrt{\delta ϵ\pi }\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}{\left(m^2+t^2\delta \pi \right)}^2\hfill \\ & \times & C_2\left)\right)/\left({\left(\frac{m^2ϵ}{M^2{t}^2}\right)}^{\frac{3}{2}}\sqrt{\delta ϵ\pi }\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}{\left(m^2+t^2\delta \pi \right)}^2\right)\left)\right)/\left(9t^3\right(m\hfill \\ & -& \frac{3m\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\left)\right),\hfill \end{array}$$

and ${\omega }_{EA}^{\prime }$ is given as follows

$$\begin{array}{ccc}\hfill {\omega }_{EA}^{\prime }& =& 3(6-\frac{4\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}\sqrt{\delta ϵ\pi }}+m\left(-9+\frac{27\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right)+\sqrt{\delta ϵ\pi }\hfill \\ & \times & \left(-\frac{9\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}+\frac{4m^2\sqrt{\frac{m^2ϵ}{M^2{t}^2}}\left(m^4+4m^2{t}^2\delta \pi +t^4{\delta }^2{\pi }^2\right)}{\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}{\left(m^2+t^2\delta \pi \right)}^3}\right)\hfill \\ & +& \frac{4t\sqrt{\delta }\sqrt{\pi }\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}\mathrm{ArcTan}\left(\frac{m}{t\sqrt{\delta }\sqrt{\pi }}\right)}{m\sqrt{\frac{m^2ϵ}{M^2{t}^2}}\sqrt{\delta ϵ\pi }}-\frac{3\sqrt{3}\sqrt{K\delta ϵ\pi }\mathrm{ArcTan}\left(\frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta ϵ\pi }}\right)}{KM})\hfill \\ & +& \frac{\left(9-\frac{4\sqrt{3}\sqrt{K}}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}}\right)C_2}{\sqrt{K}M}{\left(9m\left(m-\frac{3m\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right)\right)}^{-1}.\hfill \end{array}$$

The expression for squared speed of sound turns out as

$$\begin{array}{ccc}\hfill v_s^2& =& -3(-4+\frac{\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}\sqrt{\delta ϵ\pi }}+m\left(6-\frac{18\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right)+\sqrt{\delta ϵ\pi }\hfill \\ & \times & \left(\frac{6\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}-\frac{m^2\sqrt{\frac{m^2ϵ}{M^2{t}^2}}\left(m^4+4m^2{t}^2\delta \pi +11t^4{\delta }^2{\pi }^2\right)}{\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}{\left(m^2+t^2\delta \pi \right)}^3}\right)\hfill \\ & -& \frac{t\sqrt{\delta }\sqrt{\pi }\sqrt{\frac{m^2\delta {ϵ}^2\pi }{M^2{t}^2}}\mathrm{ArcTan}\left(\frac{m}{t\sqrt{\delta }\sqrt{\pi }}\right)}{m\sqrt{\frac{m^2ϵ}{M^2{t}^2}}\sqrt{\delta ϵ\pi }}+\frac{2\sqrt{3}\sqrt{K\delta ϵ\pi }\mathrm{ArcTan}\left(\frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta ϵ\pi }}\right)}{KM})\hfill \\ & +& \frac{\left(-6+\frac{\sqrt{3}\sqrt{K}}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}}\right)C_2}{\sqrt{K}M}{\left(18\left(m-\frac{3m\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}\right)\right)}^{-1}.\hfill \end{array}$$

The plot of the EoS parameter is shown in Figure 7 with respect to z. All the trajectories of the EoS parameter represent the quintessence phase of the universe. Figure 8 shows the graph of ${\omega }_{EA}-{\omega }_{EA}^{\prime }$ plane for same range of z. The trajectories of ${\omega }_{EA}^{\prime }$ describe the negative behavior for all ${\omega }_{EA}<0$ give the freezing region of the universe. To check the stability of the underlying model, Figure 9 shows the unstable behavior of the model. However, for $m=2$, we get some stable points for $z<-0.475$.

Exponential form of scale factor:

Taking into account second scale factor Equation (15) along with $F\left(K\right)$, we get the following energy density and pressure

$$\begin{array}{ccc}\hfill {\rho }_{EA}& =& \frac{3C^2{t}^{-2+2\theta }{\alpha }^2{\theta }^2\left(-3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }+\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)\right)}{8\sqrt{K}\pi \left(K{M}^2+3\delta ϵ\pi \right)}\hfill \end{array}$$

(34)

$$\begin{array}{ccc}\hfill p_{EA}& =& {\displaystyle \frac{1}{24\pi }}{C}^2{t}^{-3+\theta }\alpha \theta (9t^{1+\theta }\alpha \theta \left(-1+{\displaystyle \frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}}\right)-\frac{1}{M^2}ϵ(-1\hfill \\ & +& \theta \left)\right(3Mt(M\left(2-{\displaystyle \frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}}\right)-\frac{\mathrm{ArcTan}\left({\displaystyle \frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta ϵ\pi }}}\right)}{K}\hfill \\ & \times & \sqrt{3}\sqrt{K\delta ϵ\pi }+{\displaystyle \frac{C_2}{\sqrt{K}}}){\left(ϵ\right)}^{-1}+(\delta \pi \left(3t^{\theta }\alpha {\delta }^{\frac{3}{2}}{ϵ}^2\theta {\pi }^{\frac{3}{2}}\right(t^{1+\theta }\alpha \sqrt{\delta }\theta \sqrt{\pi }\hfill \\ & \times & \left(t^{2\theta }{\alpha }^2{\theta }^2-t^2\delta \pi \right)+{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^2\mathrm{ArcTan}\left({\displaystyle \frac{t^{-1+\theta }\alpha \theta }{\sqrt{\delta }\sqrt{\pi }}}\right))\hfill \\ & -& \sqrt{3}Mt\sqrt{\delta ϵ\pi }\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^2{C}_2\left)\right)/(\alpha \theta \hfill \\ & \times & \sqrt{{\displaystyle \frac{t^{-2+2\theta }ϵ}{M^2}}}{\left(\delta ϵ\pi \right)}^{\frac{3}{2}}\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^2\left)\right)).\hfill \end{array}$$

(35)

The EoS parameter is obtained from the above energy density and pressure. This parameter with its derivative are given by

$$\begin{array}{ccc}\hfill {\omega }_{EA}& =& \left(\sqrt{K}{t}^{-1-\theta }\left(K{M}^2+3\delta ϵ\pi \right)\right(9t^{1+\theta }\alpha \theta \left(-1+{\displaystyle \frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}}\right)\hfill \\ & -& \frac{1}{M^2}ϵ(-1+\theta )\left(3Mt\right(M\left(2-{\displaystyle \frac{3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)}}\right)-\sqrt{3}\sqrt{K\delta ϵ\pi }\hfill \\ & \times & \frac{\mathrm{ArcTan}\left({\displaystyle \frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta ϵ\pi }}}\right)}{K}+\frac{C_2}{\sqrt{K}}){\left(ϵ\right)}^{-1}+(\delta \pi \left(3t^{\theta }\alpha {\delta }^{\frac{3}{2}}{ϵ}^2\theta {\pi }^{\frac{3}{2}}\right(t^{1+\theta }\alpha \sqrt{\delta }\hfill \\ & \times & \theta \sqrt{\pi }\left(t^{2\theta }{\alpha }^2{\theta }^2-t^2\delta \pi \right)+{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^2\mathrm{ArcTan}\left({\displaystyle \frac{t^{-1+\theta }\alpha \theta }{\sqrt{\delta }\sqrt{\pi }}}\right))\hfill \\ & -& \sqrt{3}Mt\sqrt{\delta ϵ\pi }\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^2{C}_2\left)\right)/(\alpha \theta \sqrt{ϵ}\hfill \\ & \times & \sqrt{{\displaystyle \frac{t^{-2+2\theta }}{M^2}}}{\left(\delta ϵ\pi \right)}^{\frac{3}{2}}\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^2\left)\right)\left)\right)/(9\hfill \\ & \times & \alpha \theta \left(-3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }+\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)\right)).\hfill \end{array}$$

(36)

$$\begin{array}{ccc}\hfill {\omega }_{EA}^{\prime }& =& \left(t^{-5\theta }(-1+\theta )\left(K{M}^2+3\delta ϵ\pi \right)\right(3(-{\displaystyle \frac{3\sqrt{K}M{t}^{3\theta }{\alpha }^3{\theta }^4\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{K{M}^2+3\delta ϵ\pi }}\hfill \\ & +& {\displaystyle \frac{1}{{ϵ}^2{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^3}}K{M}^3{t}^{2+\theta }\alpha \theta \sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}(2t^{6\theta }{\alpha }^6ϵ{\theta }^7\alpha \theta \hfill \\ & \times & \sqrt{{\displaystyle \frac{t^{-2+2\theta }ϵ}{M^2}}}+t^6{\delta }^2{\pi }^2(2\delta ϵ\theta \sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}\pi +(1-2\theta )\sqrt{\delta ϵ\pi }\alpha ϵ\theta \hfill \\ & \times & \sqrt{{\displaystyle \frac{t^{-2+2\theta }\delta \pi }{M^2}}})+2t^{4+2\theta }{\alpha }^2\delta {\theta }^2\pi (3\delta ϵ\theta \sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}\pi -4(-1\hfill \\ & +& \theta )\sqrt{\delta ϵ\pi }\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}})+t^{2+4\theta }{\alpha }^4{\theta }^4(6\delta ϵ\theta \sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}\pi \hfill \\ & +& (-1+2\theta )\sqrt{\delta ϵ\pi }\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}\left)\right)-{\delta }^{3/2}{\pi }^{3/2}(\sqrt{3}{t}^{3\theta }{\alpha }^3{ϵ}^{3/2}{\theta }^4\hfill \\ & \times & \sqrt{K\delta ϵ\pi }\mathrm{ArcTan}\left({\displaystyle \frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta }\sqrt{ϵ}\sqrt{\pi }}}\right)+K{M}^3{t}^3(1-2\theta )\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}\hfill \\ & \times & \sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}\mathrm{ArcTan}\left({\displaystyle \frac{t^{-1+\theta }\alpha \theta }{\sqrt{\delta }\sqrt{\pi }}}\right)\left){\left({\left(\delta ϵ\pi \right)}^{3/2}\right)}^{-1}\right)+\alpha \theta (3\hfill \\ & \times & \sqrt{K}{t}^{3\theta }{\alpha }^2ϵ{\theta }^3-\sqrt{3}K{M}^2{t}^{2+\theta }\sqrt{{\displaystyle \frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}(-1+2\theta )\left)C_2{\left(ϵ\right)}^{-1}\right)\hfill \\ & \times & )/\left(9\sqrt{K}M{\alpha }^5{\theta }^5\left(-3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }+\sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)\right)\right).\hfill \end{array}$$

(37)

The correspond expression for $v_s^2$ is given by

$$\begin{array}{ccc}\hfill v_s^2& =& \left(t^{-\theta }\left(K{M}^2+3\delta ϵ\pi \right)\right(3({\displaystyle \frac{3\sqrt{K}M\left(-2+\theta +6t^{\theta }\alpha \theta \right)\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }}{K{M}^2+3\delta ϵ\pi }}\hfill \\ & +& KM(4-{\displaystyle \frac{t^6{\delta }^3{\pi }^3\sqrt{\frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}\sqrt{\delta ϵ\pi }{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^3}}+t^{6\theta }{\alpha }^6{\theta }^6\sqrt{\delta ϵ\pi }\hfill \\ & \times & {\displaystyle \frac{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{M^2{\left(\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}\right)}^{3/2}{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^3}}+8M^2{t}^6\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}\hfill \\ & \times & {\left(\delta ϵ\pi \right)}^{\frac{3}{2}}{\displaystyle \frac{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{{ϵ}^3{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^3}}+2\theta (-1-3t^{\theta }\alpha -4M^2{t}^6\alpha \theta {\left(\delta ϵ\pi \right)}^{\frac{3}{2}}\hfill \\ & \times & \sqrt{{\displaystyle \frac{t^{-2+2\theta }ϵ}{M^2}}}{\displaystyle \frac{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2\delta {ϵ}^2{\theta }^2\pi }{M^2}}}{{ϵ}^3{\left(t^{2\theta }{\alpha }^2{\theta }^2+t^2\delta \pi \right)}^3}}\left)\right)+K{t}^{-1+\theta }\alpha \theta \sqrt{\delta ϵ\pi }\alpha ϵ\theta \sqrt{{\displaystyle \frac{t^{-2+2\theta }\delta \pi }{M^2}}}\hfill \\ & \times & {\displaystyle \frac{\mathrm{ArcTan}\left(\frac{t^{-1+\theta }\alpha \theta }{\sqrt{\delta }\sqrt{\pi }}\right)}{M\sqrt{\delta }{\left(\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}\right)}^{3/2}\sqrt{\pi }}}+\sqrt{3}(-2+\theta )\sqrt{K\delta ϵ\pi }\mathrm{ArcTan}\left({\displaystyle \frac{\sqrt{K}M}{\sqrt{3}\sqrt{\delta ϵ\pi }}}\right))\hfill \\ & +& \sqrt{K}\left(6-3\theta -{\displaystyle \frac{\sqrt{3}\sqrt{K}}{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}}\right)C_2\left)\right)/\left(18\sqrt{K}M\alpha \theta \right(-3\sqrt{\delta ϵ\pi }\sqrt{K\delta ϵ\pi }\hfill \\ & +& \sqrt{K}\left(K{M}^2+3\delta ϵ\pi \right)\left)\right).\hfill \end{array}$$

(38)

Figure 10 represents the graph of the EoS parameter versus z for the RHDE model taking an exponential form of the scale factor. For $\alpha =2$, initially the trajectory expresses the transition from decelerated phase to accelerated phase and then crosses the phantom divide line and gives the phantom phase of the universe. For higher values of the $\alpha$, that is for $\alpha =3,4$, the trajectories of the EoS parameter represents the quintessence phase. In Figure 11, we plot the graph of evolution parameter of EoS versus EoS parameter which gives the freezing region of the universe. Figure 12 shows the graph of $v_s^2$ for stability analysis of the model. Initially the graph gives the stability and then for decreasing z, the model becomes unstable. As we increase the value of $\alpha$, the trajectories give more stable points.

6. Reconstruction from the Sharma-Mittal Holographic Dark Energy Model

Sharma-Mittal introduced two parametric entropy and defined it as [32]

$$S_{SM}=\frac{1}{1-r}\left({\left({\Sigma }_{i=1}^n{P}_i^{1-\delta }\right)}^{1-r/\delta }-1\right),$$

(39)

where r is a new free parameter. The expression of the SMHDE model for the Hubble horizon is given by

$${\rho }_D=\frac{3ϵH^4}{8\pi R}\left({(1+\frac{\delta \pi }{H^2})}^{\frac{R}{\delta }}-1\right).$$

(40)

By comparing the energy densities of the SMHDE model and Einstein-Aether gravity model, we find

$$\frac{dF}{dK}-\frac{F}{2K}=\frac{K{M}^2}{24ϵ\pi R}\left({\left(1+\frac{3ϵ\pi \delta }{K{M}^2}\right)}^{\frac{R}{\delta }}-1\right),$$

(41)

which leads us to the following solution

$$F\left(K\right)=\frac{K^2{M}^2\left(-1+{}_2{F}_1(-\frac{3}{2},-\frac{R}{\delta },\frac{-1}{2},\frac{-3\pi \delta ϵ}{K{M}^2})\right)}{36\pi Rϵ}+C_3\sqrt{K}.$$

(42)

Power-law form of scale factor:

For this scale factor, we obtain

$$\begin{array}{ccc}\hfill {\rho }_{EA}& =& \frac{K{m}^2{M}^2\left(-1+{\left(1+\frac{3\pi \delta ϵ}{K{M}^2}\right)}^{\frac{R}{\delta }}\right)}{8\pi R{t}^2}.\hfill \end{array}$$

(43)

$$\begin{array}{ccc}\hfill p_{EA}& =& {\displaystyle \frac{1}{72t^4}}m\left({\displaystyle \frac{1}{\pi R\left(m^2+\pi t^2\delta \right)}}\right(3m^2(-8\left(m^2+\pi t^2\delta \right)\hfill \\ & +& 3{\left(1+{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)}^{\frac{R}{\delta }}\left(3m^2+\pi t^2(-2R+3\delta )\right))ϵ+K{M}^2\hfill \\ & \times & t^2\left(m^2+\pi t^2\delta \right)\left(-4+9m-3(-1+3m){\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{\frac{R}{\delta }}\right)\hfill \\ & -& \left(m^2+\pi t^2\delta \right)(3m^2ϵ{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)\hfill \\ & -& K{M}^2{t}^2{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)\left)\right)\hfill \\ & -& {\displaystyle \frac{12t^2ϵ\left(-3+\frac{\sqrt{3}\sqrt{K}}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}}\right)C_3}{\sqrt{K}}}).\hfill \end{array}$$

(44)

The cosmological parameters are given by

$$\begin{array}{ccc}\hfill {\omega }_{EA}& =& {\displaystyle \frac{1}{9Km{M}^2{t}^2\left(-1+{\left(1+\frac{3\pi \delta ϵ}{K{M}^2}\right)}^{\frac{R}{\delta }}\right)}}\pi R\left({\displaystyle \frac{1}{\pi R\left(m^2+\pi t^2\delta \right)}}\right(3m^2\hfill \\ & \times & \left(-8\left(m^2+\pi t^2\delta \right)+3{\left(1+{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)}^{\frac{R}{\delta }}\left(3m^2+\pi t^2(-2R+3\delta )\right)\right)\hfill \\ & \times & ϵ+K{M}^2{t}^2\left(m^2+\pi t^2\delta \right)(-4+9m-3(-1+3m){\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{\frac{R}{\delta }}\hfill \\ & \times & )-\left(m^2+\pi t^2\delta \right)(3m^2ϵ{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)\hfill \\ & -& K{M}^2{t}^2{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)\left)\right)-12t^2ϵ\hfill \\ & \times & {\displaystyle \frac{\left(-3+\frac{\sqrt{3}\sqrt{K}}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}}\right)C_3}{\sqrt{K}}}),\hfill \end{array}$$

(45)

$$\begin{array}{ccc}\hfill {\omega }_{EA}^{\prime }& =& m^4ϵ(-3{\left(1+{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)}^{\frac{R}{\delta }}(5m^4+2m^2\pi t^2(-3R+5\delta )+{\pi }^2{t}^4(4R^2\hfill \\ & -& 10R\delta +5{\delta }^2\left)\right)-{\left(m^2+\pi t^2\delta \right)}^2(-16\hfill \\ & +& {}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)\left)\right)-4\sqrt{3}{M}^2\pi R{t}^4\hfill \\ & \times & {\left(m^2+\pi t^2\delta \right)}^2\sqrt{{\displaystyle \frac{m^2ϵ}{M^2{t}^2}}}{C}_3\left(3K{m}^4{M}^2{t}^2{\left(m^2+\pi t^2\delta \right)}^2\right(-1+(1\hfill \\ & +& {\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}{)}^{\frac{R}{\delta }}\left)\right),\hfill \end{array}$$

(46)

$$\begin{array}{ccc}\hfill v_s^2& =& {\displaystyle \frac{1}{18K^{\frac{3}{2}}{m}^3{M}^2{t}^2{\left(m^2+\pi t^2\delta \right)}^2\left(-1+{\left(1+\frac{3\pi \delta ϵ}{K{M}^2}\right)}^{R/\delta }\right)}}(\sqrt{K}{m}^2\hfill \\ & \times & \left(3m^2\right(-32{\left(m^2+\pi t^2\delta \right)}^2+3{\left(1+{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)}^{R/\delta }(11m^4+2m^2\pi t^2\hfill \\ & \times & (-5R+11\delta )+{\pi }^2{t}^4\left(4R^2-14R\delta +11{\delta }^2\right)\left)\right)ϵ-2K{M}^2{t}^2(m^2\hfill \\ & +& \pi t^2\delta {)}^2\left(4-9m+3(-1+3m){\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{R/\delta }\right)-{\left(m^2+\pi t^2\delta \right)}^2\hfill \\ & \times & (3m^2ϵ{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^2\delta }{m^2}}\right)-2K{M}^2{t}^2\hfill \\ & \times & {}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)\left)\right)+12\pi R{t}^2(m^3\hfill \\ & +& m\pi t^2\delta {)}^2ϵ\left(6-{\displaystyle \frac{\sqrt{3}\sqrt{K}}{\sqrt{\frac{m^2ϵ}{M^2{t}^2}}}}\right)C_3),\hfill \end{array}$$

(47)

We plot the EoS parameter for the SMHDE model with respect to the redshift parameter as shown in Figure 13 for the power-law scale factor. For $m=3$ and 4, the trajectories represent the transition from quintessence to phantom phase while $m=2$ indicates the phantom era throughout for z. The plot of this parameter with its evolution parameter is given in Figure 14, which shows the freezing region of the evolving universe. However, for higher values of m, we may get thawing region $({\omega }_{EA}^{\prime }>0)$. Figure 15 gives the graph of squared speed of sound versus redshift. The trajectory for $m=2$ shows the stability of the model as redshift parameter decreases while other trajectories describe the unstable behavior of the model.

Exponential form of scale factor:

Following the same steps, we obtain the following expressions for energy density, pressure and parameters for exponential scale factor. These are:

$$\begin{array}{ccc}\hfill {\rho }_{EA}& =& \frac{K{M}^2{t}^{-2+2\theta }{\alpha }^2\left(-1+{\left(1+\frac{3\pi \delta ϵ}{K{M}^2}\right)}^{R/\delta }\right){\theta }^2}{8\pi R},\hfill \end{array}$$

(48)

$$\begin{array}{ccc}\hfill p_{EA}& =& {\displaystyle \frac{1}{24}}{t}^{-4+\theta }\alpha \theta (-{\displaystyle \frac{3K{M}^2{t}^{2+\theta }\alpha \left(-1+{\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{R/\delta }\right)\theta }{\pi R}}-ϵ(-1+\theta )\hfill \\ & \times & \left(t^2\right(K^{3/2}{M}^2(-4+3{\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{R/\delta }\hfill \\ & +& {}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right))+36\pi RϵC_3)\hfill \\ & \times & {\left(3\sqrt{K}\pi Rϵ\right)}^{-1}+t^{2\theta }{\alpha }^2{\theta }^2(-8-3{\left(1+{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right)}^{R/\delta }(\pi t^2(-2R\hfill \\ & +& 3\delta )+3t^{2\theta }{\alpha }^2{\theta }^2){(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2)}^{-1}\hfill \\ & +& {}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right){\left(\pi R\right)}^{-1}-4\sqrt{3}ϵ\hfill \\ & \times & {\displaystyle \frac{C_3}{M^2{\left(\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}\right)}^{3/2}}}\left)\right)),\hfill \end{array}$$

(49)

$$\begin{array}{ccc}\hfill {\omega }_{EA}& =& {\displaystyle \frac{1}{3K{M}^2\alpha \left(-1+{\left(1+\frac{3\pi \delta ϵ}{K{M}^2}\right)}^{R/\delta }\right)\theta }}\pi R{t}^{-2-\theta }(-3K{M}^2{t}^{2+\theta }\alpha (-1\hfill \\ & +& {\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{R/\delta })\theta {\left(\pi R\right)}^{-1}-ϵ(-1+\theta )(t^2\left(K^{3/2}{M}^2\right(-4\hfill \\ & +& 3{\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{R/\delta }+{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right))\hfill \\ & +& 36\pi RϵC_3){\left(3\sqrt{K}\pi Rϵ\right)}^{-1}+t^{2\theta }{\alpha }^2{\theta }^2(-8-3{\left(1+{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right)}^{R/\delta }(\pi t^2\hfill \\ & \times & (-2R+3\delta )+3t^{2\theta }{\alpha }^2{\theta }^2){(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2)}^{-1}\hfill \\ & +& {}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right){\left(\pi R\right)}^{-1}-4\sqrt{3}ϵ\hfill \\ & \times & {\displaystyle \frac{C_3}{M^2{\left(\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}\right)}^{3/2}}}\left)\right)),\hfill \end{array}$$

(50)

$$\begin{array}{ccc}\hfill {\omega }_{EA}^{\prime }& =& {\displaystyle \frac{1}{9K{M}^2{\alpha }^2\left(-1+{\left(1+\frac{3\pi \delta ϵ}{K{M}^2}\right)}^{R/\delta }\right){\theta }^2}}\pi R{t}^{-2(1+\theta )}(-1+\theta )\hfill \\ & \times & \left({\displaystyle \frac{1}{\pi R{\left(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2\right)}^2}}\theta \right(K{M}^2{t}^2\left(-4+3{\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{R/\delta }\right)\hfill \\ & \times & {\left(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2\right)}^2+3t^{2\theta }ϵ\theta (8(-2+\theta ){\left(\pi t^2\alpha \delta +t^{2\theta }{\alpha }^3{\theta }^2\right)}^2-3{\alpha }^2\hfill \\ & \times & {\left(1+{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right)}^{R/\delta }(t^{4\theta }{\alpha }^4{\theta }^4(-5+2\theta )+2\pi t^{2+2\theta }{\alpha }^2{\theta }^2(R(3-2\theta )\hfill \\ & +& \delta (-5+2\theta ))+{\pi }^2{t}^4\left(2R\delta (5-4\theta )+4R^2(-1+\theta )+{\delta }^2(-5+2\theta )\right)\left)\right)\hfill \\ & +& {\left(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2\right)}^2(K{M}^2{t}^2{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)\hfill \\ & -& 3t^{2\theta }{\alpha }^2ϵ\theta {(-1+2\theta )}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right)\left)\right)\hfill \\ & +& 12t^2ϵ\left({\displaystyle \frac{3\theta }{\sqrt{K}}}+{\displaystyle \frac{\sqrt{3}(1-2\theta )}{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}}\right)C_3),\hfill \end{array}$$

(51)

$$\begin{array}{ccc}\hfill v_s^2& =& {\displaystyle \frac{1}{18K{M}^2\alpha \left(-1+{\left(1+\frac{3\pi \delta ϵ}{K{M}^2}\right)}^{R/\delta }\right)\theta }}\pi R{t}^{-2-\theta }({\displaystyle \frac{1}{\pi R{\left(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2\right)}^2}}\hfill \\ & \times & \left(K{M}^2{t}^2{\left(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2\right)}^2\right(-8+4\theta +18t^{\theta }\alpha \theta -3{\left(1+{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)}^{R/\delta }\hfill \\ & \times & \left(-2+\theta +6t^{\theta }\alpha \theta \right))+3t^{2\theta }ϵ(8(-4+3\theta ){\left(\pi t^2\alpha \delta \theta +t^{2\theta }{\alpha }^3{\theta }^3\right)}^2-3{\alpha }^2\hfill \\ & \times & {\left(1+{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right)}^{R/\delta }{\theta }^2(t^{4\theta }{\alpha }^4{\theta }^4(-11+8\theta )+2\pi t^{2+2\theta }{\alpha }^2{\theta }^2(R(5-4\theta )\hfill \\ & +& \delta (-11+8\theta ))+{\pi }^2{t}^4\left(2R\delta (7-6\theta )+4R^2(-1+\theta )+{\delta }^2(-11+8\theta )\right)\left)\right)\hfill \\ & -& {\left(\pi t^2\delta +t^{2\theta }{\alpha }^2{\theta }^2\right)}^2(K{M}^2{t}^2(-2+\theta ){}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{3\pi \delta ϵ}{K{M}^2}}\right)\hfill \\ & +& 3t^{2\theta }{\alpha }^2ϵ{\theta }^2{}_2{F}_1\left(-{\displaystyle \frac{3}{2}},-{\displaystyle \frac{R}{\delta }},-{\displaystyle \frac{1}{2}},-{\displaystyle \frac{\pi t^{2-2\theta }\delta }{{\alpha }^2{\theta }^2}}\right)\left)\right)+12t^2ϵC_3\hfill \\ & \times & {\displaystyle \frac{\left(6-3\theta -\frac{\sqrt{3}\sqrt{K}}{\sqrt{\frac{t^{-2+2\theta }{\alpha }^2ϵ{\theta }^2}{M^2}}}\right)}{\sqrt{K}}}).\hfill \end{array}$$

(52)

For the exponential scale factor for the SMHDE model, the plot of the EoS parameter in Figure 16 represents the phantom behavior initially, but converges to the cosmological constant behavior for $\alpha =2,3$, as z decreases. For $\alpha =4$, the EoS parameter gives the phantom behavior. Figure 17 represents the graph of the ${\omega }_{EA}^{\prime }-{\omega }_{EA}$ plane, which shows the positive behavior of ${\omega }_{EA}^{\prime }$ versus negative ${\omega }_{EA}$ expressing thawing region of the universe. The squared speed of sound graph gives unstable behavior of the SMHDE model in the framework of the Einstein-Aether theory of gravity as shown in Figure 18.

7. Summary

In this work, we have discussed Einstein-Aether gravity and utilized its effective density and pressure. We have developed the Einstein-Aether models by using some holographic dark energy models. In the presence of a free function $F\left(K\right)$, we have treated the affective density and pressure as DE. From the modified HDE models such as the THDE, RHDE and SMHDE models, we have formed the unknown function $F\left(K\right)$ for the Einstein-Aether theory by considering the power-law form and exponential forms of scale factor. We have discussed some cosmological parameters, like the EoS parameter with its evolutionary parameter and squared speed of sound to check the stability of the reconstructed models for this theory.

The remaining results have been summarized as follows:

• EoS parameter for power-law scale factor:
  • THDE ⇒ phantom behavior,
  • RHDE ⇒ quintessence phase,
  • SMHDE ⇒ transition from quintessence to phantom phase for $m=3,4$, phantom era for $m=2$.
• EoS parameter for exponential scale factor:
  • THDE ⇒ transition from quintessence to phantom era for $\alpha =2$, phantom behavior for $\alpha =3$, $\Lambda$CDM model for $\alpha =4$,
  • RHDE ⇒ phantom phase for $\alpha =2$, quintessence phase for $\alpha =3,4$
  • SMHDE ⇒ cosmological constant behavior for $\alpha =2,3$, phantom behavior for $m=4$.
• ${\omega }^{\prime }$-$\omega$ plane for power-law scale factor:
  • THDE ⇒ freezing region for $m=2,3$, thawing region for $m=4$,
  • RHDE ⇒ freezing region,
  • SMHDE ⇒ freezing region.
• ${\omega }^{\prime }$-$\omega$ plane for exponential scale factor:
  • THDE ⇒ freezing region to thawing region,
  • RHDE ⇒ freezing region,
  • SMHDE ⇒ thawing region.
• Squared speed of sound for power-law scale factor:
  • THDE ⇒ unstable,
  • RHDE ⇒ unstable,
  • SMHDE ⇒ stable for $m=2$, unstable for $m=3,4$.
• Squared speed of sound for exponential scale factor:
  • THDE ⇒ unstable,
  • RHDE ⇒ stability for higher values and instability for lower values,
  • SMHDE ⇒ unstable.

It is mentioned here that for $m=2$ for the power-law form of the scale factor in the case of the SMHDE model, we obtain a phantom region with stable behavior in the freezing region which leads to the most favorable result within the current cosmic expansion scenario.