Reconstruction from scalar-tensor theory and the inhomogeneous equation of state in f(T) Gravity

General relativity (GR) characterizes gravity as a geometric properly exhibited as curvature on spacetime. Teleprallelism describes gravity through torsional properties, and can reproduce GR at the level of equations. Similar to f(R) gravity, on taking a generalization, f(T) gravity can produce various modifications its gravitational mechanism. The resulting field equations are inherently distinct to f(R) gravity in that they are second order. In the present work, f(T) gravity is examined in the cosmological context with a number of solutions reconstructed by means of an auxiliary scalar field. To do this, various forms of the Hubble parameter are considered with an f(T) lagrangian emerging for each instance. In addition, the inhomogeneous equation of state (EoS) is investigated with a particular Hubble parameter model used to show how this can be used to reconstruct the f(T) lagrangian. Observationally, both the auxiliary scalar field or exotic terms in the FRW field equations give the same results, meaning that the variation in the Hubble parameter may be interpreted as the need to reformulate gravity in some way as is done in f(T) gravity.


I. INTRODUCTION
Modified gravity is one of the two direct approaches for reproducing the late time acceleration observed in the Universe [1][2][3]. Additionally, there are other consistency problems that must eventually be tackled in the general relativity (GR) approach to gravity [4,5]. The question then becomes what reformulation of gravity should be adopted, or whether we should take an extension of GR as our starting position. This is indeed a popular approach to the problem where many works consider an extension to the GR Einstein-Hilbert action in which the predominant attempt can be represented by the f (R) gravity formulations [6,7].
One proposition that has gained interest in recent years is that of teleparallel gravity [8][9][10][11]. The concept was initiated by Einstein himself shortly after his introduction of GR [12]. At first this was simply an equivalent reformulation of GR, called the teleparallel equivalent of general relativity (TEGR), which replaces the Einstein-Hilbert lagrangian with an analogue teleparallel torsion quantity, T (known as the torsion scalar).
As with the Einstein-Hilbert action and the f (R) approach, the simplest generalization is to * jackson.said@um.edu.mt an arbitrary function of the torsion scalar, f (T ) [8]. Similarly, the TEGR formulation is equivalent to GR at the level of equations however its generalization is distinctly different from f (R) in a number of important ways. The first being that the resulting field equations continue to be second order while those of f (R) are fourth order. This has led to interesting insights both cosmologically and otherwise [8,9,11,13,14].
Any alternative to GR should answer some of the open questions in GR. From the cosmological perspective, it would be interesting to put together a consistent theory that predicts the expansion behavior of both the early and latetime universe. In fact this is one of the primary motivations behind the ΛCDM modification to Einstein's GR formulation. While the initiation of the need for modifications to Einstein's GR comes from the observation in 1998 of the accelerating expansion of the universe [1,2], there is now overwhelming observational evidence for the need to modify the Einstein formulation [15].
The goal of the present work is to produce known and new lagrangian models within the f (T ) gravity context using several Hubble parameter models [16], using an auxiliary scalar field as a conduit to perform this reconstruction. This scalar field, essentially, takes on the role of dark energy, and would be the mediator through which these cosmological effects would take over.
Unless stated otherwise, geometric units are used where G = 1 = c. Also Latin indices are used to refer to local inertial coordinates while Greek ones are used to refer to global coordinates. The paper is divided as follows, in section II, f (T ) gravity is introduced with some focus on the distinction between f (R) and f (T ) gravity, and cosmological effects. In section III, the reconstruction work is presented along with the particular lagrangians that reproduce the Hubble parameters being considered. The inhomogeneous EoS is then derived in section IV where a particular example is given along with the general approach to determining f (T ) lagrangians. Finally the results are summarized and discussed in section V.

II. f (T ) GRAVITY AND COSMOLOGY
GR and its variants are largely based on the metric tensor, g µν , in terms of fundamental dynamical variables. This tensor acts as a potential quantity while the curvature is in turn represented through the Levi-Civita connection (which is torsion-free), Γ λ µν . In teleparallelism this connection is replaced by the Weitzenböck connection,Γ λ µν . This new connection is curvature-free and is based on two fundamental dynamical variables, namely the tetrads (or vierbein) and the spin connection. The tetrads, e a µ , are four orthonormal vectors that relate inertial and global frames in that they build the metric up from the Minkowski metric by means of an application of this transformation. In real terms, they represent the physical observer and can be related to the metric by means of where η ab = diag(1, −1, −1, −1). The tetrads obey the following inverse relations On the other hand, the spin connection, ω b aµ , is not a tensor and is dependent on the particular system under consideration, i.e. it accounts for the coordinate system in a way retains the covariance of f (T ) gravity [9]. Indeed in the literature there is a division between tetrads. There are those that produce a vanishing spin connection, ω b aµ = 0, called pure tetrads while in the impure tetrad frames the spin connection takes on some nonzero values.
In this work we take the flat FRW metric where x i represent Cartesian coordinates. The natural choice of tetrad for this metric is e a µ = diag(1, a(t), a(t), a(t)).
Since this tetrad has vanishing spin connection components [9], no contributions needs to be included in what follows. In turn, the Weitzenböck connection takes on the form,Γ λ µν = e λ a ∂ µ e a ν [10]. This naturally leads to the torsion tensor The difference between the Weitzenböck and the Levi-Civita connections is characterized by the contortion tensor Lastly the superpotential tensor is introduced which is produced purely for mathematical convenience [8]. These tensors can be combined to form the torsion scalar, T = T a µν S µν a , which is the lagrangian for TEGR.
It is at this point that the distinction between curvature and this torsion-based theory of gravity can be made clearer. In GR we adopt the Ricci scalar, R, as the gravitational lagrangian [17] whereas in the TEGR case the torsion scalar is taken. The difference between the two quantities obviously lies in a boundary term since they produce the same theory at the level of equations [8]. The distinction can be quantified through [18] where B = 2 e ∂ µ eT λ λµ = 2∇ µ T λ λµ is the boundary term. Therefore, taking the lagrangian −T + B will exactly reproduce the Ricci scalar. However, this similarity is lost once the theory is generalized to an arbitrary function thereof, that is, unless the function is f (R) = f (−T + B). In all other cases, f (R) and f (T ) will not coincide in any meaningful way. Moreover, the ensuing field equations are unique in that, out of the three possible quantities involved, namely R, T , and B, f (T ) is the only lagrangian that produces second order field equations [18,19]. Now, using the tetrad in Eq.(4), the torsion scalar turns out to be [9] Generalizing the action to an arbitrary function of the torsion scalar, f (T ), results in where κ = 4πG and e = det e a µ . On taking a variation with respect to the tetrad the following field equations emerge [8] , f T and f T T denote the first and second derivatives of f (T ) with respect to T , and L m is the matter lagrangian.
The resulting Friedmann equations turn out to be [9] and where dots denote derivatives with respect to cosmic time.
With these governing equations in hand, an auxiliary scalar field can now be introduced to reconstruct the f (T ) lagrangian.

III. RECONSTRUCTION IN f (T ) GRAVITY USING SCALAR-TENSOR THEORY
In general it is difficult to determine specific cosmologies directly from the modified Friedmann equations in Eqs. (13,14) due to the general nature of these relations. For this reason a technique first developed for f (R) gravity in Ref. [20] is investigated in the f (T ) setting. Following Refs.( [20][21][22]) an auxiliary scalar field, φ, without a kinetic term is explored. This approach has shown to be very effective in determining solutions to otherwise insoluble equations. In this context consider the action where P (φ) and Q(φ) are arbitrary functions. Taking a variation with respect to the tetrads, e a µ , gives The benefit of taking this scalar field form is that a second variation can be taken, one with respect to the scalar field itself. This results in the simple relation where primes denote derivatives with respect to the scalar field, φ. This second relation can in turn be used to determine a relation between the scalar field and the torsion scalar, φ = φ(T ). Replacing the result in the action in Eq. (15) gives a way to resolve the lagrangian in the original action in Eq.(10) The scalar field governing relations in Eq.(16) can be taken for any setting, in particular we take them for the metric in Eq.(3) which results which are the Friedmann equations for the scalar field.
The scalar field can easily be taken to be the coordinate time, φ = t, and setting the EoS to be p m = ω m ρ m . The stress-energy tensor has not been altered in any significant way so the same conservation equation naturally follows, namelẏ ρ m + 3H(1 + ω m )ρ m = 0. Solving this relation naturally leads to the cosmic time dependence of the energy density The solution for ρ m is important because with this in hand, the Friedmann equations for the scalar field can be used to determine the scalar field action components as a function of the Hubble parameter. However, in this case the Hubble parameter will be a function of the scalar field, H = g(φ), which is where the relationship between the torsion scalar and the scalar field comes in. Using Eqs. (16) and taking the coupling parameters as κ = 1 = ρ m,0 for convenience, the scalar field functions can be determined. Since the theory is second order, the scalar field turns out to be soluble in terms of the Hubble parameter. After solving the coupled differential equations where k is some constant, and For any given scalar field functional dependence, P = P (φ), Eq.(21) returns a Lambert-W function so this is not considered. Instead various Hubble parameter profiles are taken in order to investigate the reconstruction advantages of this approach. In the following these Hubble parameter model are introduced along with the associated motivation, the Lagrangian functions are then worked out.

Model 1
First a two term Hubble parameter which takes different forms at early and late times is considered. Take [16,20] where t s is an arbitrary characteristic time that represents the time at which the big rip would occur. At early times this Hubble parameter leads to a period of accelerated expansion similar to inflation which is followed directly by a decelerated period of expansion. As t → t s , the expansion of the universe becomes super-accelerated leading to the big rip. This is one straightforward way to reconcile the early-and late-times behavior of the universe.
For early times, t << t s , it straightforwardly follows that Similar to the concept of inflation, for a time t close to zero, the acceleration parameter is positive,ä a > 0. The model then naturally enters a period of deceleration [16] while the condition t << t s is retained. The scalar field functions can then be determined which will turnout to fix the f (T ) lagrangian. Firstly, by Eq.(21) it is found that and by Eq. (22) it follows that By using Eq.(9) the scalar field can be expressed in terms of the torsion scalar, φ 2 = −12H 1 /T . In turn, this leads to a functional form for the original f = f (T ) lagrangian through the scalar field in the lagrangian of Eq.(15), giving This represents the active part of the lagrangian that would be needed to account for the early time segment of the universe in f (T ) gravity. Now for very late times, approaching the big rip, i.e t → t s resulting in Hubble and acceleration parameters As the big rip time is approached the universe enters a phase of very rapid expansion again [16]. As shown above, the scalar field functions can be determined through Eq.(21) and Eq.(22) and analogously Inverting the torsion scalar in Eq.(9) results in a reconstructed function for the torsion lagrangian In this case the additional component of the lagrangian is again exponential however in this case it is positive since it plays a role at late times. The intermediary phase of the universe would still be accounted for by the TEGR component. The end result is a lagrangian for f (T ) gravity given by where the a i and A i are positive constants. The resulting theory would satisfy the cosmological requirements. However more work would have to be done to understand the behavior of the theory on smaller scales such as for galaxies and the solar system. One way to get a more constrained lagrangian from this approximate form would be to compare supernova data and reduce the arbitrariness of the free parameters. Model 2 In model 1 the very late-time rapid acceleration takes on a phantom nature. Another approach is to consider [16] where H 0 and H 1 > 0 are constants, and n is a positive integer. The n > 1 region is considered since it gives three well defined regions corresponding to an early acceleration phase (interpreted as inflation), then deceleration, and finally late-time acceleration. The corresponding acceleration parameter turns out to bë For the model under consideration we do not need to take limits to determine the f (T ) lagrangian, so the scalar functions can straightfor-wardly be determined giving and Inverting the torsion scalar relation in Eq.(9) results in a scalar field The arbitrary function f (T ) can then be determined in general As in the first model the TEGR term is still included to govern the intermediate stage of evolution of the universe, resulting in the gen- . The square root term actually emerges naturally in other models [11,13,14] however when it appears on its own it does not contribute to the eventual Friedmann equations whereas in this case it is implicit in another function. In the current context it is interesting for it to turn out to form part of the lagrangian solution. To a much lesser extent this is not dissimilar to the Gauss-Bonnet generalizations [23].

Model 3
Lastly, the following model is considered [16] where H i , H l , c and α are positive constants. In model 2 additional scalar fields have to be introduced to suppress issues with inflation whereas in the current context a clear link with observation is attainable. For the early and late universe the hubble parameter tends to the constants H i and H l respectively. In this way, H i would drive inflation and H l would take on the small cosmological constant for late times.
The corresponding acceleration parameter turns out to bë a a = Putting this into the scalar field functions in the lagrangian in Eq. (21) and Eq.(22) results in and Q(φ) = 10k As in the previous case, the relationship between the scalar field and the torsion scalar follows straightforwardly Through the mechanics of the scalar field, φ, Eq.(18) is utilized to produce the teleparallel lagrangian Incorporating TEGR again gives a relatively simple lagrangian that contains different lagrangian elements for the different epochs of the cosmological history. In this instance this gives the model f ( which has two effective Hubble parameters for early and late times. This last model can work as a bridge between the first two model since it contains both terms in √ −T and also the exponential term in T . Altogether a wide variety of potential lagrangian terms emerge from considering this collection of Hubble parameter ansatz.

IV. THE INHOMOGENEOUS EQUATION OF STATE
In the previous section different Hubble functions are used to reconstruct the gravitational lagrangian using a scale factor approach as the vehicle for reconstruction. In this section the additional terms that f (T ) offers are interpreted as an extra dark fluid. The resulting EoS of the dark fluid depends on the Hubble parameter and its derivatives, in an inhomogeneous way.
Taking the modified Friedmann equations from Eqs. (13,14) 3H where no other matter contributions are admitted. This would represent a cosmology dominated by modified gravity terms. On comparison with the standard Friedmann equations from GR, the energy and pressure densities contributions can be easily identified as where a perfect fluid context is being assumed. That is, the modified terms can be interpreted as forming part of a cosmic perfect fluid that takes on the role of dark energy. Following this line of thought a dark fluid EoS naturally emerges as which can be represented as where the second term is dependent only on the Hubble parameter (κ = 1). This is the inhomogeneous EoS for the dark fluid which has been investigated in a number of works [20,[24][25][26][27]. Following Refs. [20,28], this equation can be generalized to where g H,Ḣ,Ḧ, .
fT . On combining the Friedmann equations and using the inhomogeneous EoS gives the following differential equation in f (T ) where κ is taken as unity for convenience in what follows. This equation can be solved only once a Hubble parameter is assigned since it is essentially a linear first order differential equation in f (T ).
As an example the following Hubble parameter is considered [28] where H 0 and H 1 are positive constants. The Hubble model gives an early decelerating period followed up an accelerating period, which are studied in in turn as limiting cases. Another reason why this is done is because the relationship between the torsion scalar and cosmic time is non-invertible in the general case. This is needed to eliminate any appearances of cosmic time as t = t(T ) in the eventual f (T ) lagrangian. For the initial decelerating period, H(t) ∼ H 1 /t, giving a time dependence and solving Eq.(50) gives where λ is an integration constant. It is not uncommon for f (T ) gravity to produce lagrangians that cannot be put in closed form [8,29,30]. Thus by using Eq.(49), the inhomogenous EoS equation that turns out to reproduce this kind of behavior is given by On the other hand, for late-times the Hubble parameter takes the form H(t) ∼ H 0 t, and combining with Eq.(9) results in the time dependence (55) Again solving the differential equation Eq.(50) turns out to give a lagrangian whereλ is an integration constant. This general type of lagrangian has shown some promise in cosmological settings [8,29]. Taking derivatives gives the general inhomogeneous EoS equation as This indeed reproduces the desired late time behavior.

V. DISCUSSION AND CONCLUSION
In this work a number of different solutions of f (T ) theory have emerged from investigating the cosmological scenario within the context of both early and late time epochs. This was achieved through the vehicle of using an auxiliary scalar similar to the approach [20]. This method works very well for generalized theories because it gives an second way of solving the Friedmann equations.
Secondly, the inhomogeneous EoS was also derived using a particular Hubble parameter, with both early and late time reconstructions. Together these two approaches offer two methods of using an auxiliary scalar field to reconstruct the f (T ) lagrangian. Moreover, both the auxiliary scalar field and exotic modified gravity terms interpreted as a perfect fluid EoS result in the same state parameter behavior.
f (T ) theories may provide a suitable description of the expansion history of the universe without the need to include exotic components to the stress-energy part. The models would need to be constrained at lower scales, such as in the galactic and solar system regime. Following this the models should next be compared with cosmological data such as supernovae luminosity distance and the CMB profile, similar to the initiative in f (R) theory, however this is out of the scope of this paper.
Beyond the models considered for the first method another two were attempted, namely g(φ) = φ α and g(φ) = h 2 0 1 t 2 0 −φ 2 + 1 t 2 1 +φ 2 however in either case the problem becomes intractable very fast. Other approaches may include perturbative analyses.
For the models that do emerge, firstly they are predominantly exponential functions but some are not. In a follow-up study, it would be interesting to understand better the transition between the different epochs in terms of the lagrangian terms presented here. This would naturally involve an analysis of the stability of each model.