Cosmological Analysis of Reconstructed $\mathcal{F}(T,T_\mathcal{G})$ Models

In this paper, we analyze cosmological consequences of the reconstructed generalized ghost pilgrim dark energy $\mathcal{F}(T,T_\mathcal{G})$ models in terms of redshift parameter $z$. For this purpose, we consider power-law scale factor, scale factor for two unified phases and intermediate scale factor. We discuss graphical behavior of the reconstructed models and examine their stability analysis. Also, we explore the behavior of equation of state as well as deceleration parameters and $\omega_{\Lambda}-\omega_{\Lambda}^{'}$ as well as $r-s$ planes. It is found that all models are stable for pilgrim dark energy parameter $2$. The equation of state parameter satisfies the necessary condition for pilgrim dark energy phenomenon for all scale factors. All other cosmological parameters show great consistency with the current behavior of the universe.


Introduction
The accelerated cosmic expansion phenomenon is undoubtedly the biggest achievement of the twentieth century. The source behind this expansion is said to be a repulsive force called dark energy (DE) having large negative pressure. This energy is evenly scattered in the universe but we do not know much about its nature as well as composition. It is suggested by WMAP experiment that the universe has budget as 73% DE, 23% dark matter and 4% baryonic matter. The cosmic expansion goes through different stages of dark matter and DE normally characterized by the equation of state (EoS) parameter. These ranges include ω < −1 for phantom, ω = −1 for vacuum (cosmological constant (Λ)) and −1 < ω < 1 3 for quintessence DE dominated eras. The matter dominated eras corresponding to ω = 0, 1 and 1 3 represent cold dust, stiff and radiation dominated eras, respectively.
The cosmological constant is the best ingredient to discuss the DE mystery but it has issues like coincidence and fine-tuning. This motivates researchers to find some alternatives to describe the DE nature. The most appealing proposals in this scenario are either to modify matter or gravitational side of the Einstein-Hilbert action. The matter modification provides different DE models like Chaplygin gas, phantom, quintessence, k-essence, holographic and pilgrim dark energy (PDE) etc [1]- [5]. On the other hand, the gravitational modification leads to modified theories. Among these theories, there is a modification based on torsional formation of general relativity dubbed as teleparallel theory. In this theory, the basic entity is torsion instead of curvature.
The generalization of teleparallel theory is known as f(T) theory in which torsion scalar T is replaced by an arbitrary function f(T) in the action. Recently, an extension of this theory is proposed by introducing teleparallel equivalent Gauss-Bonnet term T G known as F (T, T G ) theory. The reason behind this generalization is to develop an action that includes higher torsion correction terms. Kofinas and Saridakis [6] presented this distinct theory by evaluating a torsion equivalent of Gauss-Bonnet term without using curvature formalism. They analyzed several observable like EoS, DE density and matter density parameters by assuming two specific F (T, T G ) models [7] and found this theory as explaining the cosmic evolution.
Kofinas et al. [8] discussed dynamical analysis of spatially flat FRW metric by considering a particular F (T, T G ) model and concluded that the universe can exhibit various DE dominated solutions such as cosmological constant, quintessence or phantom like solutions that depend upon the values of corresponding model parameters. Waheed and Zubair [9] investigated energy bounds with perfect fluid using Hubble, deceleration, jerk and snap parameters. Zubair and Jawad [10] explored laws of thermodynamics at apparent horizon of FRW metric. Jawad [11] studied energy conditions for FRW universe analytically.
Various DE models have been developed in the context of quantum gravity as well as general relativity. One of the them is the Veneziano ghost DE model defined as ρ T = µH, where µ is a constant [12]. This model is interesting as it does not involve any new degree of freedom or new parameter. The Veneziano ghost energy density has the form H + O(H 2 ) that provides enough amount of vacuum energy to explore the expansion phenomenon. However, ghost DE model involves only the term H in its energy density. Therefore, Cai et al. [13] added the term H 2 in the ghost DE model as ρ T = µH + νH 2 , where ν is another constant, known as generalized ghost DE model. Fernandez [14] discussed ghost DE models along with scalar field whereas Malekjani [15] established different f(R) models considering ghost as well as generalized ghost DE models.
Wei [16] presented another DE model dubbed as pilgrim DE (PDE) motivated by the fact that phantom like DE has enough strength to prevent black hole formation rather than other types of DE. The PDE also encouraged this fact due to its same repulsive nature. The generalized ghost DE density is further established by involving PDE parameter as  [20] investigated ghost DE model in F (T, T G ) gravity and examined its cosmological consequences through the reconstructed model. Sharif and Nazir [21] reconstructed GGPDE F (T, T G ) models and discussed their corresponding EoS parameters versus PDE parameter.
In this paper, we study cosmological behavior of the reconstructed models [21] versus redshift parameter z and discuss their stability through squared speed of sound parameter. We investigate these reconstructed models through EoS parameter, deceleration parameter, ω Λ − ω ′ Λ analysis and r − s plane. The paper is arranged as follows. Next section provides basic introduction of F (T, T G ) gravity. In section 3, we briefly describe the well-known scale factors. Section 4 analyzes the evolution trajectories via cosmological parameters. In the last section, we summarize the results.

F(T, T G ) Gravity
In this section, we provide a concise review of F (T, T G ) gravity in the background of FRW geometry. The tetrad field e A (x α ) has a fundamental role in f (T ) as well as F (T, T G ) gravity. Trivial tetrad is the simplest one expressed as e A = ∂ α δ α A and e B = ∂ α δ α B , where δ α a is the Kronecker delta. These are not commonly used because they provide zero torsion. The non-trivial tetrad have different behavior, so they are more supportive in describing teleparallel theory. These tetrad can be represented as The metric tensor can also be expressed in the product of tetrad fields as The Weitzenböck connection ω A B (x α ) that describes parallel transportation, has the following form The structure coefficients C C AB are defined as . Similarly, we can express the torsion as well as curvature tensors as The contorsion tensor is defined by Finally, the torsion scalars T and T G take the form The action for F (T, T G ) gravity is proposed by Kofinas and Saridakis [6] where L m is the matter Lagrangian and κ 2 = 1. The teleparallel equivalent to general relativity is obtained by substituting F (T, T G ) = −T . We can also have Gauss-Bonnet theory when F (T, T G ) = αT G − T , where α represents Gauss-Bonnet coupling. The F (T, T G ) field equations can be obtained by varying the action as where and Now we discuss cosmological significance of F (T, T G ) theory by considering flat FRW universe model as where a(t) is the scale factor. There exist infinite possible tetrad fields for each metric, thus we choose a common tetrad field for FRW metric as where dual is defined as The corresponding torsion scalars are Here, H =˙a a defines Hubble parameter and dot indicates time derivative. Substituting the above values in Eq.(2), we obtain We can rewrite the above equations in usual form as where the energy density and pressure for DE sector are The energy conservation equations in terms of dark matter and DE arė

Cosmic Scale Factors
Here, we briefly describe some scale factors through which we explore the cosmological behavior of our reconstructed models [21].

• Power-Law Scale Factor
This scale factor is defined as [22] a(t) = a 0 t n , where n > 0, a 0 > 0. This form of scale factor provides a great consistency for flat FRW metric with the supernova data. For n > 1, it gives an accelerating universe. Using this scale factor, we obtain the corresponding values as and Eq.(1) becomes

• Scale Factor for Unified Phases
The following scale factor unifies matter as well as DE dominated phases. The Hubble parameter takes the form as [22]- [24] which leads to the following form of scale factor as When t is very small, we obtain H(t) ∼ H 2 t which exhibits the presence of perfect fluid with ω Λ = 2 3 H −1 2 − 1. Moreover, when t is very large, H → H 1 yielding constant Hubble parameter which leads to de Sitter universe. This type of Hubble parameter yields a transition from matter to DE dominated phases. The corresponding values of torsion scalars are Using Eqs. (1) and (10), it follows that The scale factor and the corresponding Hubble parameter are defined as [25] a(t) = exp(bt m ), where 0 < m < 1 and b is an arbitrary constant. This scale factor is much useful in cosmological analysis as it has great consistency with astrophysical observations. Both torsion scalars for this scale factor take the form The corresponding energy density of GGPDE is

Cosmological Analysis Via Well-Known Scale Factors
In this section, we explore the behavior of the reconstructed models and investigate their stability. We discuss EoS as well as deceleration parameters and ω Λ − ω ′ Λ as well as r − s planes by using the above three scale factors. We consider the reconstructed F (T, T G ) models [21] obtained by equating the corresponding energy densities of GGPDE model and F (T, T G ) gravity. For this purpose, we equate Eqs.(1) and (7), i.e., ρ Λ = ρ T as We can determine solution of the above equation only for a particular choice of the scale factor. Thus we consider all the above three scale factors and analyze their behavior.
• Now, we investigate stability of the reconstructed GGPDEF (T, T G ) model through the squared speed of sound v 2 s defined as The positive value (v 2 s > 0) indicates stability of the model whereas its negative value (v 2 s < 0) corresponds to instability of the model. Using Eqs. (7), (8) and (13)    • The evolutionary behavior of EoS parameter for the reconstructed F (T, T G ) model is analyzed by evaluating ω Λ through Eqs. (7) and (8) as follows Substituting Eq.(13) in the above expression, we obtain the EoS parameters for u = 2 (Figure 3 left) and u = −2 (right) in terms of z. We consider same values of the corresponding constants as taken earlier. We investigate the evolution of EoS parameter in the interval −0.9 ≤ z ≤ 2 for n 1 = 3.2, n 2 = 4 and n 3 = 5. Figure 3 (left plot) shows that ω Λ starts from phantom region, cuts the phantom divide line and at the end, it becomes zero. This means that EoS parameter shows quintom behavior for u = 2. Similarly, the right plot represents that ω Λ starts from dust like matter era, passes via quintessence as well as vacuum DE eras and finally enters in the phantom era. Hence, ω Λ behaves like quintom for u = −2. In both cases, the reconstructed F (T, T G ) model satisfies PDE phenomenon.
• The deceleration parameter q is described as Its positive value indicates decelerating behavior, q = 0 expresses constant expansion and negative value corresponds to accelerating universe. Substituting the value of ω Λ in the above equation, we obtain deceleration parameter (Figure 4). The left plot for u = 2 indicates that q attains negative values in the range −0.9 ≤ z < 0.15, hence represents accelerating universe in this interval. At z = 0.15, it becomes zero showing constant behavior and for z > 0.15, it leads to decelerating universe. In the right plot of Figure 4 (u = −2), the deceleration parameter exhibits negative values in the interval −0.9 ≤ z ≤ 1.5 showing accelerating behavior.
• The plane ω Λ − ω Λ ′ is developed for examining different DE models. Initially, Caldwell and Linder [28] used this method to study the behavior of quintessence DE model. They suggested that the covered area in the phase plane corresponds to two regions, thawing region (ω Λ < 0, ω Λ ′ > 0) as well as freezing region (ω Λ < 0, ω Λ ′ < 0). Here, we discuss ω Λ − ω Λ ′ plane corresponding to u = 2, −2 for three different values of n. Figure 5 (u = 2) shows that ω Λ − ω ′ Λ plane represents thawing regions for all three values of n. Similarly, all the curves exhibit the same behavior for u = −2 in the interval −0.9 ≤ z ≤ 1 as shown in Figure 6. Hence, ω Λ − ω ′ Λ plane shows consistency with the current behavior of the universe for n = 3.2, 4 and 5.
• Many DE models have been suggested to understand the phenomenon of DE that ultimately explain the current behavior of the universe. Some of them provide same values of the Hubble and deceleration parameters. Thus it is necessary to determine which one gives better information about acceleration of the expanding universe. For this purpose, Sahni et al. [29] presented two dimensionless parameters named as statefinder parameters and are defined as .
We can also write the parameter r in terms of q as r = 2q 2 +q−q ′ . These parameters describe the well-known cosmic regions such as (r, s) = (1, 1) indicating CDM (cold dark matter) limit and (r, s) = (1, 0) showing ΛCDM limit. The region s > 0, r < 1 describes phantom and quintessence eras while s < 0, r > 1 represents Chaplygin gas model. Here, we establish r − s planes corresponding to our reconstructed GG-PDEF (T, T G ) models for the same three values of n. We assume the same values of corresponding parameters as in the previous section in the range −0.5 ≤ z ≤ 5. Figures 7 and 8 show that both reconstructed models for u = 2 and −2 provide the regions of quintessence and phantom DE eras as s > 0 and r < 1 for all three values of n.

Scale Factor for the Unified Phases
Now, we investigate the behavior ofF (T, T G ) models through the scale factor for unified phases. For this purpose, we substitute Eq.(10) in (12) with a = a 0 (1 + z) −1 to get a differential equation. We obtain the reconstructed  GGPDEF (T, T G ) model in terms of z by solving this differential equation numerically. For both u = 2 as well as −2, we consider H 1 = 0.9 and three different values of H 2 = 2.7, 2.75, 2.8. Figure 9 shows that both reconstructed models represent increasing behavior as the value of z increases.
• Figure 10 represents the behavior of squared speed of sound parameter for u = 2 and −2. Both plots show that positive values of v 2 s for z > 1.02 confirm the stability ofF(T, T G ) models.
• Figure 11 indicates that EoS parameter starts from matter dominated era initially (ω Λ = 0) for u = 2 as well as −2. At z = 0, ω Λ crosses the phantom divide line for u = 2 except when H 2 = 2.7. For u = −2, the EoS parameter remains in the matter dominated era for all values of H 2 and represents phantom dominated era for z > 2.2. However, in both cases, the EoS parameter shows consistency with the current expanding behavior of the cosmos as the value of z increases.
• Figure 12 shows deceleration parameter in terms of redshift parameter z. This is zero in the interval −0.9 ≤ z ≤ 1.1 representing the constant cosmic expansion in both cases u = 2 as well as −2. However, for z > 1.1, negative values of q give rise to accelerating universe.  • The behavior of ω Λ −ω ′ Λ is shown in Figure 13. The left plot represents freezing region for H 2 = 2.7 and 2.75 while it corresponds to thawing region for H 2 = 2.8. In the right plot, ω Λ −ω ′ Λ expresses thawing region for H 2 = 2.75 and 2.8.
• The behavior of statefinder parameters is shown in Figure 14. In the left plot (u = 2), we notice that s < 0, r > 1 showing Chaplygin gas model while the right plot (u = −2) indicates s > 0, r < 1 implying phantom and quintessence eras of the universe.

Intermediate Scale Factor
For this scale factor, we obtain reconstructed GGPDEF (T, T G ) models numerically by using Eq.(11) in (12) as shown in Figure 15. We investigate the behavior of our models by assuming three different values of m as m 1 = 0.33,      • Figure 16 confirms the stability of our corresponding model in the range −0.9 ≤ z ≤ 1.15 for u = 2, while for u = −2, the model shows instability.
• Figure 17 indicates that for both values of the PDE parameter, the EoS parameter exhibits transition from phantom towards matter dominated era by crossing the phantom divide line for all three values of m.
• Figure 18 implies that both plots attain negative values in the range −0.9 ≤ z ≤ 1.8 which describe the accelerating behavior of the universe.
• Figure 19 shows the plots for intermediate ω Λ − ω ′ Λ plane that correspond to thawing regions.
• In Figure 20 for u = 2 (left), we attain the point (r, s) = (1, 0) which shows ΛCDM limit. Also, the statefinder parameters represent phantom and quintessence regions for 0 ≤ z ≤ 1. On the other hand, these parameters indicate Chaplygin gas model for 1 ≤ z ≤ 2.5. Similarly, the right plot for u = −2 shows Chaplygin gas model.

Concluding Remarks
In this paper, we have investigated the behavior of reconstruction models [21] in F (T, T G ) gravity along with GGPDE model as well as three scale factors.
For this purpose, we have considered two values of PDE parameter, i.e., u = 2 and −2. We have explored the role of some cosmological parameters versus redshift parameter z in this scenario. We have observed that our models show stability for all the scale factors except power-law and intermediate for u = −2. The equation of state parameter in both cases (u = 2, −2) represents quintom behavior for all the scale factors. It is found that reconstructed GGPDEF (T, T G ) models fulfill the condition of PDE phenomenon. The plot of deceleration parameter versus z exhibits accelerated expansion of the universe.
We have observed that ω Λ − ω ′ Λ plane displays thawing region for powerlaw as well as intermediate scale factor. The scale factor for two unified phases provides freezing region for H 2 = 2.8. The r − s plane shows phantom and quintessence regions for power-law model (u = 2) as well as for the scale factor of unified phases (u = −2). For intermediate case, we have achieved ΛCDM limit. We have found that both the planes, i.e., ω Λ − ω Λ ′ as well as r − s are consistent with the current cosmic behavior.
Sharif and Jawad [17] analyzed GGPDE model by investigating its cosmological consequences in general relativity. Our results are consistent with them for non-interacting case with µ = 1.55 and ν = 1.91. Jawad and Rani [18] investigated the reconstructed models, their stability and EoS parameter in the modified Horava-Lifshitz f (R) gravity. Our results are in great agreement with their work. Sharif and Nazir [19] discussed the same cosmological parameters for GGPDE model in f(T) gravity. We have also compared our results with [19] and found that the EoS parameter and cosmological planes for u = 2 represent consistency with the same values of model parameters. We have noticed that EoS parameter is also consistent with the observational data [30] given as ω Λ = −1.13 +0.24