Solar system tests for realistic f(T) models with non-minimal torsion–matter coupling

In the previous paper, we have constructed two f(T) models with non-minimal torsion–matter coupling extension, which are successful in describing the evolution history of the Universe including the radiation-dominated era, the matter-dominated era, and the present accelerating expansion. Meantime, the significant advantage of these models is that they could avoid the cosmological constant problem of Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM. However, the non-minimal coupling between matter and torsion will affect the tests of the Solar system. In this paper, we study the effects of the Solar system in these models, including the gravitation redshift, geodetic effect and perihelion precession. We find that Model I can pass all three of the Solar system tests. For Model II, the parameter is constrained by the uncertainties of the planets’ estimated perihelion precessions.


Introduction
General relativity (GR), whose centenary was recently celebrated [1,2], uses the metric tensor as the fundamental dynamical variable, and chooses the torsionless Levi-Civitá connection to describe the gravitation field. In the teleparallel equivalent of general relativity (TEGR) [3][4][5], the curvatureless Weizenböck connection is chosen instead, and the tetrad field plays the role of fundamental variable. TEGR provides the gauge structure of gravitation: the gauge theory for the group of translation on the tangent bundle of the spacetime. In this sense, teleparalellism may be considered as an approach to including gravity into the unification of gauge theories and maybe, further, to quantize the gravity.
Under the framework of TEGR, different modifications of gravity theory have been proposed for better understanding of the currently accelerating Universe [6,7]. One possible scheme is f (T ) theory, in which the torsion scalar T in the a e-mail: zhaixh@shnu.edu.cn b e-mail: kychz@shnu.edu.cn gravitation Lagrangian is replaced by an arbitrary function f (T ) [8][9][10][11][12][13][14]. Compared to f (R) theory, the field equations are second order instead of fourth one, which is an important advantage of f (T ) theory. As a further extension, and an analogy to that of f (R) theories [15][16][17][18][19][20][21][22][23][24][25][26][27][28][29], non-minimal torsion-matter coupling f (T ) gravity has been proposed and studied [30]. In the previous research [31], using the observation data of type Ia supernovae (SNeIa), cosmic microwave background (CMB), and baryon acoustic oscillations (BAO), we established two concrete f (T ) models with non-minimal torsion-matter coupling extension, and we found that they are successful in describing the observation of the Universe and its large-scale structure and evolution. The joint fitting led to m0 = 0.255 ± 0.010, b0 h 2 = 0.0221 ± 0.0003 and H 0 = 68.54 ± 1.27 for model I and m0 = 0.306 ± 0.010, b0 h 2 = 0.0225±0.0003 and H 0 = 60.97±0. 44 for model II at 1σ confidence level, where H 0 is the Hubble parameter at present, h = H 0 /(100km/s/Mpc), and m0 and b0 are the density parameters of dark matter and baryon matter. In a recent work, we also studied the structure formation of these two models [32].
On the other hand, any modified gravity theory must confront the Solar system tests which have been passed in GR [33][34][35][36][37]. In fact, Solar system effects in some modified gravities including minimal coupling f (T ) models have been considered [38][39][40][41][42][43][44]. In Ref. [42], the authors have investigated the perihelion precession, light bending, Shapiro time delay and gravitation redshift in minimal coupling f (T ) gravity with f (T ) = T + αT n , and constrained the parameter α for n = 2, 3. However, non-minimal geometry-matter coupling theories predict that the movements of massive particles are no longer geodesic, which will have influences on some of the Solar system effects. Similar cases have been studied in some non-minimal coupling f (R) theories (see e.g. [16,45,46]).
In this paper, we consider the Solar system tests including the gravitation redshift, geodetic effect, and perihelion precession in the two models we have constructed in Ref. [31]. The gravitation redshift is one representative effect that does not involve massive particle. Since the two models we are considering assume that the coupling between radiation and torsion is still minimal for simplicity, this kind of effects is practically the same as the minimal coupling theory. The geodetic effect and perihelion precession, however, involve the movement of massive objects, and thus should be affected by the non-minimal coupling between matter and torsion. Judging from this, it is necessary to combine the cosmological restriction with Solar system effects for the parameters in two models. We have gotten the cosmological restriction using the observation data of SNeIa, CMB and BAO as follows [31]: A = 0.188 ± 0, 048, B = 0.510 ± 0.060 (for Model I); A = 0.633 ± 0.012 and B is a free parameter (for Model II). It is our task in this paper to check whether or not these values of parameters will pass the Solar system tests.
The paper is organized as follows: In Sect. 2 we briefly review the f (T ) theory with power law torsion-matter coupling extension and our two cosmologically fitted models. In Sect. 3, we calculate the weak field limit of the spherically symmetric background metric. In Sects. 4, 5, and 6 we consider the gravitation redshift, geodetic effect and perihelion precession in our two models, respectively. We conclude our research in Sect. 7.

Brief review of the theory
On the spacetime differentiable manifold M, one can generally find dual pairs of linearly independent frame fields {e a , e a }, such that e a (e b ) = δ a b . One particular pair would be the gradients {∂ a } of the tangent space coordinates {x a } and their covectors {dx a }. One such frame field, also known as the tetrad or vierbein field, forms a base for the vectors on the tangent space T p M at each point p ∈ M. On the overlap area of different patches of M, the members of one base can be written in terms of members of the other, e.g. e a = e μ a ∂ μ or e a = e a μ dx μ .
The spacetime metric g is related to the tangent space metric η by or conversely And hence the determinant TEGR uses these vierbein fields as dynamical variables to define the Weitzenböck connection [47] λ νμ ≡ e λ a ∂ μ e a ν = −e a ν ∂ μ e λ a .
The torsion and contorsion tensors are then given by and Utilizing the tensor one can define the torsion scalar TEGR [3][4][5] This means a c = e a μ e μ c is one of the entries of a Lorentz transformation. That is, different vierbein fields satisfying the Lorentz transformation relation will give the same spacetime metric. Though this does not make much difference in TEGR, it does have some significant effects on f (T ) theories, one of which is that the f (T ) theories are no longer Lorentz invariant theories [48]. The fact that different vierbein fields (related by a Lorentz group) lead to the same metric calls for an appropriate choice of vierbein for f (T ) theories and their extensions, especially when spherical coordinates is considered [49].
One of the extensions of f (T ) theories is the non-minimal torsion-matter coupling f (T ) gravity, proposed by Harko et al. [30]. The action of the non-minimal torsion-matter coupling f (T ) gravity can be written as where L m is the Lagrangian density of matter. With the action principle applied on Eq. (14) with respect to the vierbein, the equation of motion is then given by where

and
T α β is the energy-momentum tensor of matter given by And it takes the usual form for a perfect fluid, where p and ε are the pressure and energy density of the matter, respectively, and u μ is the 4-velocity.
The covariant derivative (related to the Levi-Civitá connection) of Eq. (15) gives This suggests that the energy-momentum tensor is no longer conservative. However, contracting Eq. (18) with u μ , we have If we take the matter Lagrangian density to be L m = ε [30,[50][51][52][53], then as discussed in Ref. [31], in cosmological cases, i.e. with the Friedmann-Lemaître-Robertson-Walker (FLRW) metric g μν = diag(1, −a, −a, −a), Eq. (19) will lead back to the usual conservation law of matteṙ In Ref. [31], we have constructed two concrete models with We have found that these two models are successful in describing the observation of the Universe and its large-scale structure and evolution. In order to see whether these two models can pass the Solar system tests, we will consider the gravitation redshift, geodetic effect, and perihelion precession of these two models in the following sections.

Spherically symmetric weak field solution
Firstly, we calculate the background metric. Since we are considering the spherically symmetric vacuum solution, Eq. (15) can be rewritten as with Therefore, we are actually dealing with the minimal coupling f (T ) theory since L m vanishes at outside of the central mass. As mentioned before, one should be careful when it comes to choosing vierbein fields in spherical coordinates for f (T ) theories. The general spherically symmetric metric can be written as It is found [49] that the vierbein field is a viable choice of vierbein for f (T ) theories in that it does not impose extra constraints on the form of f (T ) and it preserves the Birkhoff theorem, that is, The torsion scalar is then given by where the prime indicates a derivative with respect to r . Thus, Eq. (22) can reduce to In the weak field limit ξ, ζ 1, these equations can be solved up to linear approximation as for Model I, and for Model II, where the integral constant is chosen to be the central spherically symmetric mass M as in the Schwarzschild solution. These results are consistent with the general weak field solution for minimal coupling power law f (T ) theory given in Refs. [41,54].
With the (vacuum) background vierbein field solved, we are ready to calculate the Solar system effects.

Gravitation redshift
The gravitation redshift for a radiation signal is given by where ν, ν 0 are the frequencies of the signal measured at the position of r, r 0 , respectively. Since this effect does not involve any massive test particle, the background solutions (28) and (29) can be used directly. Thus, up to linear approximation, the gravitation redshift is for Model I and for Model II. The common terms in each case are the prediction of GR, and the last terms are the modifications from the models. The Hydrogen-Maser in Gravity Probe A redshift experiment [55] launched a spacecraft nearly vertically upward to 10,000 km to measure this effect. Using the cosmologically best-fit values (21), we have the gravitation redshift modifications δ( ν ν 0 ) I ≤ 9.045 × 10 −78 for Model I and δ( ν ν 0 ) II = 2.331B × 10 −20 for Model II. The redshift experiment reaches a δ rs = 10 −14 accuracy [55]. It is obvious that the modifications are within the error for both models if B < 10 6 for Model II.

Geodetic effect
Consider a gyroscope described by a spin 4-vector s μ orbiting around the central mass. Its 4-velocity is given by u μ = (u t , 0, 0, u φ ) where u φ is the orbiting angular velocity, and u μ u μ = 1. s μ is spacelike such that u μ s μ = 0. It is found that as long as the direct effects of tidal forces can be neglected, s μ is Fermi-Walker transported along u μ [37,56], that is, Contracting Eq. (18) with the induced metric h μλ = g μλ − u μ u ν , we have Combining Eqs. (33) and (34), and suppressing the hydrodynamical term ∇ ν p, which is not related to f 2 , we have where μ νλ is the Levi-Civitá connection and we have taken L m = ε [30,[50][51][52]. The r, θ, φ components of Eq. (35) are The first and third equations vanish because u r , u θ = 0. Since s μ u μ = 0, let θ = π/2, we can put into Eq. (36) and then we have The general solution of Eq. (38) can be written as where the angular frequency is given by On the other hand, from Eq. (34), we have Using u μ u μ = 1, we have . (42) Thus, the precessing rate of the gyroscope ω = u φ − . For the two models being considered, up to linear approximation, it is , The first term 3M/2r is the prediction of GR, the rest terms are the modifications from the models. Thus the relative modifications are , The Gravity Probe B mission [57] tested this effect 1 with four gyroscopes at the typical altitude of 642 km. Analysis of the data resulted in a relative error of geodetic precessing rate less than 0.4%. With the cosmologically best-fit value 1 It measured also the gravitomagnetic spin precession induced in GR by the Earth's rotation; see, e.g., [58]. (21) and the mass of Earth (G M ⊕ 22487.9 eV −1 , with the gravitational constant G recovered), the relative modifications of our models are | ω I /ω GR | ≤ 3.005 × 10 −31 and ω II /ω GR = 2.152B × 10 −11 . It is obvious that the modifications from both models are within the experiment error if B < 10 8 for Model II.

Perihelion precession
For a massive test particle orbiting around the central mass, the equation of motion can be obtained from Eq. (41): Assume that the orbit lies in the θ = π/2 plane and dθ = 0. The t and φ components of the equation of motion are which can be integrated: where k, h are integral constants. Now from the line element we have Up to linear approximation, Eq. (49) can be expressed as for Model I and for Model II. The first line of each model represents the GR case. Eq. (51) can be compared to the classical equation of motion where V, E are the potential and total energy, and L is the angular momentum, and m 0 is the mass of the test particle. In our cases, both L and E are constants. Thus in the first lines of Eqs. (50) and (51), we can identify V Newt. = − M r as the Newtonian potential, E = k 2 −1 2 as the conserved energy per unit mass, l as the conserved angular momentum per unit mass, and V GR = − l 2 M r 3 as the potential from GR. The rest lines of Eqs. (50) and (51) can be seen as additional potentials from the modifications of the models. Thus, the equation of motion can be re-expressed as for both models and for Model I and for Model II. V 2 and V 4,II are constants, which do not contribute to the precession. All the other potentials are power law functions of r . Perihelion precession for such kind of potentials were calculated in, e.g., [59][60][61]. For a potential V (r ) in the form V (r ) = γ r λ , it is found that the extra perihelion precession per orbital period is where e ecc. and r a are the eccentricity and the semi major axis of the orbit, respectively, and I (e ecc. , λ) = 2π 0 cos θ dθ (1 + e ecc. cos θ ) λ+1 . For a small eccentricity e ecc. , I (e ecc. , λ) −e ecc. π(1 + λ).
Hence, the prediction of the extra perihelion precession from GR is which is the well-known result. The extra modifications of perihelion precession are for Model I and for Model II. The extra perihelion precessions of the planets (from Mercury to Saturn) in the Solar system estimated from the recent ephemerides INPOP10a [62] and EPM2011 [63] (see also [36]) are listed in Table 1.
While it is not very conclusive whether or not some of these planets have anomalous perihelion precessions [36], the preliminary bounds on the two models should be obtained by comparing Eqs. (60) and (61) with the uncertainties of the perihelion precessions in Table 1. Note that Eqs. (60) and (61) give the theoretical predictions of extra changes of the orbit after one revolution, while Table 1 gives the extra precession rates of the planets estimated from the observation. Divided by the orbital periods of the planets, the theoretical predictions from Eqs. (60) and (61) can be compared with observation. Using the semi major axes r a , eccentricities e ecc. and the orbital periods T 0 of the planets' orbits, as well as the Solar mass (G M 7.483 × 10 9 eV −1 , with the gravitational constant recovered), we can calculate the extra precession rates of Models I and II with the cosmologically best-fit values of A, B (21), which are listed in Table 2.
One can see that with the best-fit parameters, the extra perihelion precessions from Model I are less than the uncertainties in Table 1 by about 10 orders of magnitude. Thus, if the unproved anomalies are overlooked it can be concluded that Model I with the best-fit parameters can pass the test of perihelion precession. In the same spirit, if we require that for Model II the extra precessions from the model |δφ II /T 0 | listed in Table 2 are less than 10 −4 mas/year, we can obtain the upper limit of the parameter B < 7.41 × 10 −7 .

Conclusion and discussion
In this paper, we have studied the gravitation redshift, geodetic effect and perihelion precession in the two concrete models of f (T ) theory with non-minimal torsion-matter coupling extension, which have been previously constructed and investigated for cosmology. It is found that Model I can pass all three Solar system tests. For Model II, if it is required that the theoretical extra perihelion precessions are indistinguishable under the current observation accuracy and less than 10 −4 mas/year, the free parameter B should be <7.41×10 −7 , otherwise the theoretical values will surpass the bounds of the uncertainties of the planets' perihelion precessions.
Neither the gravitation redshift nor Shapiro time delay involves massive particle. Since the coupling between radiation and torsion is still minimal in Models I and II, these effects are practically the same as the minimal coupling theory [38,39,42]. Furthermore, the test of Shapiro time delay has been omitted in this paper because its constraint is looser than that of the gravitation redshift for the parameters A and B in the Models I and II.
In a nutshell, we have shown that Models I and II successfully describe not only the observation of the Universe and its large-scale structure and evolution, but also the Solar system effects of gravitation. Therefore, these models are realistic ones. However, we find no acceptable H 0 and m0 when f 1 (T ) and f 2 (T ) are both proportional to T 2 . Further research is needed as to whether or not there is a better f (T ) model with non-minimal torsion-matter coupling extension. For instance, one can consider the case of f 1 (T ) ∝ T m and f 2 (T ) ∝ T n and take optimum of (m, n) using observation at the large and small scales. We will consider this question in the future work.