Observational Constraints on $f(T)$ gravity from varying fundamental constants

We use observations related to the variation of fundamental constants, in order to impose constraints on the viable and most used $f(T)$ gravity models. In particular, for the fine-structure constant we use direct measurements obtained by different spectrographic methods, while for the effective Newton's constant we use a model-dependent reconstruction, using direct observational Hubble parameter data, in order to investigate its temporal evolution. We consider two $f(T)$ models and we quantify their deviation from $\Lambda$CDM cosmology through a sole parameter. Our analysis reveals that this parameter can be slightly different from its $\Lambda$CDM value, however the best-fit value is very close to the $\Lambda$CDM one. Hence, $f(T)$ gravity is consistent with observations, nevertheless, as every modified gravity, it may exhibit only small deviations from $\Lambda$CDM cosmology, a feature that must be taken into account in any $f(T)$ model-building.

On the other hand, in some modified cosmological scenarios one can obtain a variation of the fundamental constants, such as the fine structure and the Newton's constants. Such a possibility has been investigated in the literature since Dirac [23] and Milne and Jordan [24] times. Later on, Brans and Dicke proposed the time variation of the Newton's constant, driven by a dynamical scalar field coupled to curvature [25], while Gamow triggered subsequent speculations on the possible variation of the fine structure constant [26]. Similarly, in recent modified gravities, which involve extra degrees of freedom comparing to general relativity, one may obtain such a variation of the fundamental constants [27,28]. However, since experiments and observations give strict bounds on these variations [29][30][31][32], one can use them in order to constrain the theories at hand.
In the present work we are interested in investigating the constraints on f (T ) gravity by observations related to the variation of fundamental constants. In particular, since f (T ) gravity predicts a variation of the finestructure and Newton's constants, we will use the recent observational bounds of these variations in order to constrain the f (T ) forms as well as the range of the involved parameters. The plan of the work is the following: In Section II we give a brief review of f (T ) gravity and cosmology. In Section III we investigate the constraints on specific f (T ) gravity models arising from the observational bounds of the fine-structure constant variation, while in Section IV we study the corresponding constraints that arise from the observational bounds of the Newton's constant variation. Finally, in Section V we summarize our results.

II. f (T ) GRAVITY AND COSMOLOGY
In this section we provide a short review of f (T ) gravity and cosmology. We use the tetrad fields e µ A , which form an orthonormal base at each point of the tangent space of the underlying manifold (M, g µν ), where g µν = η AB e A µ e B ν is the metric tensor defined on this manifold (we use Greek indices for the coordinate space and Latin indices for the tangent one). Furthermore, instead of the torsionless Levi-Civita connection which is used in the Einstein-Hilbert action, we use the curvatureless [10]. Hence, the gravitational field in such a formalism is described by the following torsion tensor: Subsequently, the Lagrangian of the Teleparallel equivalent of general relativity, namely the torsion scalar T , is constructed by contractions of the torsion tensor as [10] T One may consider generalized theories in which the Lagrangian T is extended to an arbitrary function f (T ), similarly to the f (R) extension of curvature-based gravity. In particular, such gravitational action will read as where e = det(e A µ ) = √ −g, and G N is the Newton's constant. Additionally, along the gravitational action (3) we consider the matter sector, and hence the total action writes as where L m (e A µ , Ψ M ) is the total matter Lagrangian including the electromagnetic field. Finally, variation in terms of the tetrad fields give rise to the field equations as where f T = ∂f /∂T , f T T = ∂ 2 f /∂T 2 , and with T (m) ρ ν the total matter energy-momentum tensor. In the above equation we have inserted for convenience the "superpotential" tensor S µν Applying f (T ) gravity in a cosmological framework we consider a spatially flat FLRW universe with line element ds 2 = −dt 2 + a 2 (t)[dr 2 + r 2 dθ 2 + sin 2 θ dφ 2 ], which arises from the diagonal tetrad e A µ = diag(1, a(t), a(t), a(t)), with a(t) the scale factor. In this case, the field equations (5) become where ρ and p are respectively the total matter energy density and pressure, and where ρ T , p T are the effective dark-energy energy density and the pressure of gravitational origin, given by In the above expressions we have used that which arises straightforwardly from (2) in the FLRW universe. Finally, from equations (8), (9) we can define the effective dark-energy equation of state (EoS) as where as usual we use the redshift z = a0 a − 1, as the independent variable, and for simplicity we set a 0 = 1. Clearly, w has a dynamical nature.
In the following we focus on two well-studied, viable f (T ) models, which correspond to a small deviation from ΛCDM cosmology, and which according to [20][21][22] are the ones that fit the observational data very efficiently.
• The first scenario is the power-law model (hereafter where θ, b are the two free model parameters, out of which only one is independent. Inserting this f (T ) form into the first Friedmann equation (6) at present time, i.e. at redshift z = 0, one may derive that where Ω m0 = 8πGρm0 is the corresponding density parameter at present. Hence, and using additionally that ρ m = ρ m0 (1 + z) 3 , equation (6) for this model can be written as Lastly, we mention that the above model for b = 0 reduces to ΛCDM cosmology, while for b = 1/2 it gives rise to the Dvali-Gabadadze-Porrati (DGP) model [33].
• The second scenario is the square-root-exponential (hereafter f 2 CDM) of [12], with in which β and p the two free model parameters out of which only one is independent. Inserting this f (T ) form into (6) at present time, one obtains that Finally, the first Friedmann equation (6) for this model can be written as Lastly, note that this model reduces to ΛCDM cosmology for p → +∞. Hence, in the following, for this model it will be convenient to set b ≡ 1/p, and hence ΛCDM cosmology is obtained for b → 0 + .

III. OBSERVATIONAL CONSTRAINTS FROM FINE-STRUCTURE CONSTANT VARIATION
In this section we will use observational data of the variation of the fine-structure constant α, in order to constrain f (T ) gravity. Let us first quantify the α-variation in the framework of f (T ) cosmology. In general, in a given theory the fine-structure constant is obtained using the coefficient of the electromagnetic Lagrangian. In the case of modified gravities, this coefficient generally depends on the new degrees of freedom of the theory [28,34]. Even if one starts from the Jordan-frame formulation of a theory, with an uncoupled electromagnetic Lagrangian, and although the electromagnetic Lagrangian is conformally invariant, and it is not affected by conformal transformations between the Jordan and Einstein frames, thus it will acquire a dependence on the extra degree(s) of freedom due to quantum effects [35]. In particular, if φ is the extra degree of freedom that arises from the conformal transformationg µν = Ω 2 g µν from the Jordan to the Einstein frame, then quantum effects such as the presence of heavy fermions (note that this does not necessarily require new physics till the Planck scale) will induce a coupling of φ to photons, namely [35] where F µν is the electromagnetic tensor, g bare the bare coupling constant, and with M pl = 1/(8πG N ) the Planck mass and β γ = O(1) a constant (we have assumed that β γ φ M pl ). Hence, the scalar coupling to the electromagnetic field will imply a dependence of the fine structure constant of the form [34,35] 1 where the subscripts denote the Einstein and Jordan frames respectively, or equivalently The above factor is in general time-(i.e redshift-) dependent. Therefore, it proves convenient to normalize it in order to have ∆α = 0 at present (z = 0), which in Although the above procedure is straightforward in cases where a conformal transformation from the Jordan to the Einstein frame exists, it becomes more complicated for theories where such a transformation is not known.
In case of f (T ) gravity, it is well known that a conformal transformation does not exist in general, since transforming the metric asg µν = Ω 2 g µν , with Ω 2 = f T a smooth non-vanishing function of spacetime coordinates, one obtains Einstein gravity plus a scalar field Lagrangian, plus the transformed matter Lagrangian, plus a non vanishing term 2Ω −6∂µ Ω 2T ρ ρµ [36]. This additional term forbids the complete transformation to the Einstein frame, and hence in every application one has indeed to perform calculations in the more complicated Jordan one.
In order to avoid performing calculations in the Jordan frame we will make the reasonable assumption that f (T ) = T +const.+corrections, which has been shown to be the case according to observations [19][20][21][22], and holds for the two forms considered in this work, namely (12) and (15), too. Hence, Ω 2 = 1 + corrections, and theñ ∂ µ Ω 2 is negligible, which implies that the above extra term can be neglected. In the end of our investigation, we will verify the validity of the above assumption. Thus, we can indeed obtain an approximate transformation to the Einstein frame, and in particular the introduced degree of freedom reads as φ = − √ 3/f T [36]. Hence, inserting this into (19) we acquire and thus inserting into relation (22), we can easily extract the variation of the fine-structure constant as where f T 0 = f T (z = 0). Lastly, since β γ = O(1), the above relation becomes Hence, for a general f (T ), the ratio ∆α/α indeed depends on z, through the f T (z) function (we remind that according to (10), T (z) = −6H 2 (z)), while in the case of standard ΛCDM cosmology, where f (T ) = T + Λ, ∆α/α becomes zero.
z ∆α/α(ppm) Ref.  In the following, we confront relation (25) with observations of the fine-structure constant variation, in order to impose constraints on f (T ) gravity (it proves that the neglected term between (24) and (25) imposes an error of the order of 10 −9 and hence our approximation is justified). We use direct measurements of the finestructure constant that are obtained by different spectrographic methods, summarized in Table I. Additionally, along with these data sets, and in order to diminish the degeneracy between the free parameter of the models, we use 580 Supernovae data (SNIa) from Union 2.1 compilation [43], as well as data from BAO observations, adopting the three measurements of A(z) obtained in [44], and using the covariance among these data given in [45].
In the following two subsections, we analyze two viable models, namely f 1 CDM of (12) and f 2 CDM of (15), separately.
For the power-law f 1 CDM model of (12), we easily acquire where we have used also (10). Inserting (26) into (25) we can derive the evolution of ∆α/α as (27) where the ratio H 2 (z)/H 2 0 is given by (14). We mention that while analyzing the model for the data set of ∆α/α of Table I, we have marginalized over Ω m0 , and thus the statistical information focuses only on the parameter b. For the fittings ∆α/α + SNIa and ∆α/α + SNIa + BAO, we have considered Ω m0 as a free parameter, and we have found that Ω m0 = 0.23±0.13 (for ∆α/α + SNIa) and Ω m0 = 0.293 ± 0.023 (for ∆α/α + SNIa + BAO) at 1σ confidence level. Finally, in Fig. 1 we  (12), obtained from the joint analysis ∆α/α + SNIa + BAO. The cross marks the best-fit value.
present the 68.27% and 95.45% confidence regions in the plane Ω m0 − b, considering the observational data ∆α/α + SNIa + BAO. Note that these results are in qualitative agreement with those of different observational fittings [20,21], and show that ΛCDM cosmology (which is obtained for b = 0) is inside the obtained region. In fact, one may notice from Table II that the reduced χ 2 for ∆α/α + SNIa and ∆α/α + SNIa + BAO data are very close to 1, while for single data from ∆α/α its value slightly exceeds 1 although not significantly. Additionally, in order to examine the late-time asymptotic behavior of the scenario at hand, in Fig. 2 we depict the evolution of the equation-of-state parameter given in (11), applying a reconstruction at 1σ confidence level via error propagation using the joint analysis ∆α/α + SNIa + BAO. As we can see, w at late times acquires values very close to '−1', as expected. For a more detailed  investigation of the late-time asymptotics up to the far future one must apply the method of dynamical system analysis as it was done in [46], where it was thoroughly shown that the universe will end in a de Sitter phase. The evolution of the equation-of-state parameter given in (11), for the f1CDM power-law model of (12), applying a reconstruction at 1σ confidence level via error propagation using the joint analysis ∆α/α + SNIa + BAO.
In summary, it is clear that f 1 CDM model, under all the above three different combinations of statistical data sets, remains close to ΛCDM cosmology as expected. Lastly, note that this is a self-consistent verification for the validity of our assumption that f (T ) = T + const. + corrections, which allowed us to work in the Einstein frame.
For the square-root-exponential f 2 CDM model of (15), we easily obtain Inserting (28) into (25) we can derive the evolution of ∆α/α as where the ratio H 2 (z)/H 2 0 is given by (17). We mention that while analyzing the model for the data set of ∆α/α of Table I, we have marginalized over Ω m0 , and thus the statistical information focuses only on the parameter b. For the fittings ∆α/α + SN Ia and ∆α/α+SN Ia+BAO we have taken Ω m as a free parameter, and we note that Ω m0 = 0.277 ± 0.019 (for ∆α/α + SNIa) and Ω m0 = 0.283 ± 0.016 (for ∆α/α + SNIa + BAO) at 1σ confidence level. Finally, in Fig. 3 we  (15), obtained from the joint analysis ∆α/α + SNIa + BAO. The cross marks the best-fit value.
present the 68.27% and 95.45% confidence regions in the plane Ω m0 − b, considering the observational data ∆α/α + SNIa + BAO (we have taken b 0.001 in order to avoid divergences in the function H(z) at high redshifts). Note that these results are in qualitative agreement with those of different observational fittings [20,21], and show that ΛCDM cosmology (which is obtained for b → 0 + ) is inside the obtained region. Furthermore, and similarly to the f 1 CDM model, from Table III we deduce that although the data from ∆α/α alone show a slightly deviating nature (reduced χ 2 = 1.1) from ΛCDM scenario, but for ∆α/α + SNIa and ∆α/α + SNIa + BAO data it is implied that the model is very close to ΛCDM cosmology. This is also a self-consistent verification for the validity of our assumption that f (T ) = T + const. + corrections, which allowed us to work in the Einstein frame.  Lastly, in order to examine the late-time asymptotic behavior of f 2 CDM model, in Fig. 4 we depict the evolution of the equation-of-state parameter given in (11), applying a reconstruction at 1σ confidence level via error propagation using the joint analysis ∆α/α + SNIa + BAO, where one can see that w at late times acquires values very close to '−1', as expected. Similarly to the previous model, for a more detailed investigation of the late-time asymptotics one must apply the method of dynamical system analysis [46], where it can thoroughly be shown that the universe will end in a de Sitter phase. The evolution of the equation-of-state parameter given in (11), for the f2CDM square-root-exponential model of (15), applying a reconstruction at 1σ confidence level via error propagation using the joint analysis ∆α/α + SNIa + BAO.

IV. OBSERVATIONAL CONSTRAINTS FROM THE NEWTON'S CONSTANT VARIATION
In this section we will directly use the observational constraints imposed on f (T ) models in order to examine the variation of the gravitational constant G N . Let us first quantify the G N -variation in the framework of f (T ) cosmology. As is well known, a varying effective gravitational constant is one of the common features in many modified gravity theories [1]. In case of f (T ) gravity, the effective Newton's constant G ef f can be straightforwardly extracted as [21,47] Hence, for the f 1 CDM power-law model of (12), and using (26), we obtain where the ratio H 2 (z)/H 2 0 is given by (14) (clearly, for b = 0 we have that G ef f (z) = G N = const.). Similarly, for the f 2 CDM square-root-exponential model of (15), and using (28), we acquire where the ratio H 2 (z)/H 2 0 is given by (17) Let us now use the above expressions for G ef f (z) and confront them with the observational bounds of the Newton's constant variation. We use observational Hubble parameter data in order to investigate the temporal evolution of the function G ef f (z), since such a compilation is usually used to constrain cosmological parameters, due to the fact that it is obtained from model-independent direct observations. We adopt 37 observational Hubble parameter data in the redshift range 0 < z ≤ 2.36, compiled in [48], out of which 27 data points are deduced from the differential age method, whereas 10 correspond to measures obtained from the radial baryonic acoustic oscillation method.
We apply the following methodology: Firstly, we estimate the error in the measurements associated with the function G ef f /G N , for both models of (31) and (32), via the standard method of error propagation theory, namely and we fix the free parameters of the two models within the values obtained in the joint analysis of [22] and of Section III of the current work. Then, the measurements of G ef f /G N are calculated directly for each redshift defined in the adopted compilation.
In the left graph of Fig. 5 we depict the 1σ confidencelevel estimation of the function G ef f (z)/G N from 37 Hubble data points, in the case of the f 1 CDM powerlaw model of (12). Additionally, in the right graph of Fig. 5  from the observational Hubble parameter data. When we perform the analysis within the known range of the parameter b for this model (from [22] as well as from Fig. 1 above), we find that G ef f /G N ≈ 1. Nevertheless, a minor deviation is observed for the fixed value of b = 0.01 (see the left graph of Fig. 5). For in-stance, note that G ef f (z = 0.07)/G N = 0.992 ± 0.004 and G ef f (z = 2.36)/G N = 0.99923 ± 0.00005, for the first and the last data points of the redshift interval [0, 2.36], respectively.
In Fig. 6 we present the corresponding graphs for the f 2 CDM square-root-exponential model of (15). When we perform the analysis within the known range of the parameter b ≡ 1/p for this model (from [22] as well as from Fig. 3 above), we find that G ef f /G N ≈ 1, similarly to the case of f 1 CDM model.
In summary, from the analysis of this section, we verify the results of the previous section, namely that the parameter b that quantifies the deviation of both f 1 CDM and f 2 CDM models from ΛCDM cosmology, is very close to zero. These results are in qualitative agreement with previous observational constraints on f (T ) gravity, according to which only small deviations are allowed, with ΛCDM paradigm being inside the allowed region [19][20][21][22].

V. CONCLUSIONS
In the present work we have used observations related to the variation of fundamental constants, in order to impose constraints on the viable and most used f (T ) gravity models. In particular, since f (T ) gravity predicts a variation of the fine-structure constant, we used the recent observational bounds of this variation, from direct measurements obtained by different spectrographic methods, along with standard probes such as Supernovae type Ia and baryonic acoustic oscillations, in order to constrain the involved model parameters of two viable and wellused f (T ) models.
For both the f 1 CDM power-law model, as well as for the f 2 CDM square-root-exponential model, we found that the parameter that quantifies the deviation from ΛCDM cosmology can be slightly different than its ΛCDM value, nevertheless the best-fit value is very close to the ΛCDM one. Additionally, since f (T ) gravity predicts a varying effective gravitational constant, we quantified its temporal evolution with the use of the previously constrained model parameters. For both f 1 CDM and f 2 CDM models, we found that the deviation from ΛCDM cosmology, is very close to zero.
These results are in qualitative agreement with previous observational constraints on f (T ) gravity [20][21][22], however they have been obtained through completely independent analysis. In summary, f (T ) gravity is consistent with observations, and thus it can serve as a candidate for modified gravity, although, as every modified gravity, it may have only small deviation from ΛCDM cosmology, a feature that must be taken into account in any f (T ) model-building.