Parameterized post Newtonian approximation in a teleparallel model of dark energy with a boundary term

We study the parameterized post-Newtonian approximation in teleparallel model of gravity with a scalar field. The scalar field is non-minimally coupled to the scalar torsion as well as to the boundary term introduced in [1]. We show that, in contrast to the case where the scalar field is only coupled to the scalar torsion, the presence of the new coupling affects the parameterized post-Newtonian parameters. These parameters for different situations are obtained and discussed.


Introduction
In teleparallel model of gravity, instead of the torsionless Levi-Civita connections, curvatureless Weitzenböck connections are used [2]. Teleparallel equivalent of general relativity was first introduced in [3] as an attempt for unification of electromagnetism and gravity. This theory is considered as an alternative theory of usual general relativity and has been recently employed to study the late time acceleration of the Universe [4]. This can be accomplished by considering modified f (T ) models [5], where T is the torsion scalar, or by introducing exotic field such as quintessence. Assuming a nonminimal coupling between the scalar field and the torsion opens new windows in studying the cosmological evolution [6], and can be viewed as a promising scenario for late time acceleration and super-acceleration [7].
A non minimally coupled scalar field, like scalar-tensor model, may alter Newtonian potential. So it is necessary to check if the model can pass local gravitational tests such as solar system observations. This can be done in the context of the parameterized post-Newtonian formalism [10]. In [8], and [9] it was shown that when the scalar field is only coupled to the scalar torsion, there is no deviation from general relativity in the parameterized post-Newtonian (PPN) parameters and the theory is consistent with gravitational tests and solar system observations.
Recently a new coupling between the scalar field and a boundary term B, corresponding to torsion divergence B ∝ ∇ µ T µ , was introduced in [1], where the cosmological consequences of such a coupling for some simple power law scalar field potential, and the stability of the model were discussed. There was found that that the system evolves to an attractor solution, corresponding to late time acceleration, without any fine tuning of the parameters. In this framework, the phantom divide line crossing is also possible. Thermodynamics aspects of this model were studied in [11]. This model includes two important subclasses, i.e. non minimally coupled quintessence to the Ricci scalar and non minimally coupled quintessence to the scalar torsion. Another important feature of this model is its ability to describe the present cosmic acceleration in the framework of Z 2 symmetry breaking by alleviating the coincidence problem [12].
In this paper, we aim to investigate whether this new boundary coupling may affect the Newtonian potential and PPN parameters: γ(r) and β(r).
The scheme of the paper is as follows: In the second section we introduce the model and obtain the equations of motion. In the third section, we obtain the weak field expansion of the equations in the PPN formalism and obtain and discuss their solutions for spherically symmetric metric. We show that the PPN parameters may show deviation from general relativity. We consider different special cases and derive explicit solutions for the PPN parameters in terms of the model parameters and confront them with observational data.

The model and the field equations
In our study we use vierbeins e a = e a µ ∂ µ whose duals, e a µ , are defined through e a µ e a ν = δ ν µ . The metric tensor is given by g µν = η ab e a µ e b ν , η = diag(−1, 1, 1, 1). e = det(e a µ ) = det √ −g.
where k 2 = 8πG N , and G N is Newtonian gravitational constant. The torsion scalar is defined by and the boundary term is [13] where T µ = T λ λ µ . The Weitzenböck torsion, and connection are given by and respectively. S ρ µν is defined according to Note that R = −T + B, where R is the Ricci scalar curvature. Hence for χ = −ǫ the model reduces to a quintessence model coupled non-minimally to the scalar curvature, while for χ = 0, we recover the quintessence model coupled non-minimally to the scalar torsion.
By variation of the action (1) with respect to the vierbeins we obtain τ ν µ is the energy momentum tensor of matter. The trace of (7), multiplied by −δ ν µ /2, is By combining (8) and (7) we obtain Note that the trace of the energy momentum tensor is τ = g µν τ µν . In the same way, variation of the action with respect the scalar field gives Eqs. (9) and (10) are the main equations that we will work with in the following parts.

post-Newtonian formalism
To investigate the post-Newtonian approximation [10] of the model, the perturbation is specified by the velocity of the source matter | v| such that e.g. O(n) ∼ | v| n . The matter source is assumed to be a perfect fluid obeying the post-Newtonian hydrodynamics: where ρ is energy density, p is the pressure and Π is the specific internal energy. u µ is the four-vector velocity of the fluid. The velocity of the source matter is v i = u i u 0 . The orders of smallness of energy momentum tensor ingredients are [10] ρ where U is the Newtonian gravitational potential. The components of the energy momentum tensor are given by We expand the metric around the Minkowski flat background as [9,8] Note h (1) µν = 0 [10]. Accordingly, the vierbeins may be expanded as [9] Note B (1) a µ = 0. In our analysis we need non zero components of the metric up to order 4, i.e : h 00 . We also use the notation B µν = η µσ B σ ν and δ a σ B a ν = B σ ν . By comparing (14) and (15) we derive (like [9], and [8], B ij is assumed to be diagonal) We introduce two functions A, and γ (which is one of the PPN parameters) through [8] B The scalar field is expanded as where and φ 0 is a constant cosmological background. φ 0 is of order O(0) and may evolve in times of order of the Hubble time, so in solar system tests we assume that it is static. The time derivative, ∂ 0 = ∂ ∂t , of the other fields are weighted with order O(1) [10].
The potential around the background is Defining After these preliminaries, let us solve the eqs. (9) and (10) order by order in the PPN formalism. At zeroth order (9) and (10) imply 0-0 component of (9) gives which at order 2 reduces to: resulting in where the potential is given by To obtain (25), we have used . (27) By taking the trace of i-j component of (9), at order 2, we obtain: which reduces to At the second order perturbation, the boundary term B, defined in (3), is derived as Hence from (10) the equation of motion of the scalar field becomes Eqs. (25), (29), and (31) are our three main equations for determining A, γ, and ψ (2) . Using these three equations, for a given U , A is derived as and ψ (2) is obtained as γ is determined by the equation where Γ := γ−1 γ+1 . (34) is a nonhomogeneous screened Poisson equation whose solution is where Equation (32) allows us to take where G is defined through So one can define an effective G ef f. through (7) at the third order gives which by using reduces to To simplify computations one may employ gauge condition which determines B (3) 0 j in terms of second order parameters. This gauge is compatible with eqs. (25) and (29). Using one can find that (23) at the fourth order gives Also, the scalar field equation at the fourth order is By using and (2) .
(45), and (49), are our main results in the fourth order. These equations together with (42) and (43) in the third order, and (32), (33), (34), in the second order must be solved to give the post-Newtonian parameters.
To solve these complicated equations, we consider solutions specified by U = U (r) which results in The gauge (43) implies ∂ j B (3) 0 j = 0. Therefore (42) reduces to For v i = 0, (51) gives B (3) 0 i = 0 (by assumption that perturbation terms vanish at large distance). In this situation eqs. (45) and (49) and respectively. To obtain post-Newtonian parameters we must obtain A, ψ (2) , and γ(r). By inserting them in (52) and (53), we obtain solutions for B (4) 0 0 . To do so we consider a spherically symmetric metric with a point source.

φ 0 = 0
For χ = 0, we may have also a situation with no deviation in PPN parameters from general relativity, this occurs for φ 0 = 0. For example for potentials and (22)) leads to φ 0 = 0 which by using (59-61) results in γ = 1, G = 1, and β = 1. Therefore in this case too, there is no deviation from general relativity for the PPN parameters.
In this experiment the gravitational interaction, in terms of Astronomical Unit takes place at r ≃ 7.44 × 10 −3 AU [16] . The parameter β is determined by lunar laser ranging experiments via the Nordtvedt effect [17]. This test indicates the bound [15] at a gravitational interaction distance r = 1AU [16]. (79) and (80) restrict the parameters of our model. For V (φ) = 0 , (79)and (68) give In the limiting cases (71) and (72) we find and respectively.
For λr ≫ 1, we have and (80) restricts our parameters as

conclusion
The teleparallel model of gravity with a quintessence (nonminimally) coupled to the torsion and also to a boundary term (proportional to the torsion divergence) was considered (see (1)). Although the model shows some interesting aspects in cosmology and in describing the late time acceleration of the Universe, but it must also pass local gravitational and solar system tests. So we studied the parameterized post-Newtonian (PPN) approximation of the model. We obtained the equations of motion (see the second section), and solve them order by order to obtain PPN parameters (see the third section). Explicit expression for PPN parameters in spherically symmetric metric were obtained and different possible situations were discussed. Our results show that PPN parameters, except for some special cases, i.e. in the absence of boundary term and also zero scalar field background, differ from general relativity. So we conclude that, coupling of the scalar field to the boundary term generally makes the model deviate from general relativity in the PPN limit.
Since T and B are not invariant under local Lorentz transformations, the teleparallel model with the boundary term is not invariant under Lorentz transformations unless one takes χ = −ǫ. Despite this, in spacetimes with spherical symmetry like Schwarzschild spacetime and so on, it is possible to choose good or preferred tetrads to solve this issue [18]. In scalar-tetrad theories of gravity the preferred tetrads cannot be detected via measuring the metric components [19]. Similarly, in our model, PPN parameters in the standard post-Newtonian formalism do not identify the effect of preferred tetrads. To include these effects one must generalize the post-Newtonian approach, as was pointed out in [8].