Distinguishing f(R) theories from general relativity by gravitational lensing effect

The post-Newtonian formulation of a general class of f(R) theories is set up in a third-order approximation. It turns out that the information of a specific form of f(R) gravity is encoded in the Yukawa potential, which is contained in the perturbative expansion of the metric components. Although the Yukawa potential is canceled in the second-order expression of the effective refraction index of light, detailed analysis shows that the difference of the lensing effect between the f(R) gravity and general relativity does appear at the third order when f′′(0)/f′(0)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{f''(0)/f'(0)}$$\end{document} is larger than the distance d0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_0$$\end{document} to the gravitational source. However, the difference between these two kinds of theories will disappear in the axially symmetric spacetime region. Therefore only in very rare case the f(R) theories are distinguishable from general relativity by gravitational lensing effect in a third-order post-Newtonian approximation.


Introduction
Recently modified gravity theories have received increasingly attention in issues related to "dark energy" [1][2][3], "dark matter" [4][5][6][7], as well as non-trivial tests on gravity beyond general relativity (GR) [8]. Historically, Einstein's GR is the simplest relativistic theory of gravity with correct Newtonian limit. To pursue new physics, Weyl and Eddington even began to consider modifying GR just after it was established [9,10]. From the viewpoints of perturbutive quantum gravity, GR is non-renormalizable [11][12][13][14], while higherorder gravity theories might alleviate the problem. From the phenomenological viewpoints, there are many ways to modify GR, and some empirical approaches seem to have promising prospect, such as Dvali-Gabadadze-Porrati gravity [15], tensor-vector-scalar theory [16] and Einstein-Aether a e-mail: liu.hongguang@cpt.univ-mrs.fr b e-mail: mayg@bnu.edu.cn theory [17]. Among such extended theories, particular attention has been devoted to the so-called f (R)-gravity. This kind of theories is based on a generalization of the Einstein Hilbert Lagrangian to nonlinear functions f (R) of the Ricci scalar [18]. f (R)-gravity covers a lot of characteristics of higher-order gravity and is convenient to be operated. Hence, f (R) theories provide an ideal tool to study the possible extension of GR. f (R) theories of gravity can also be nonperturbatively quantized by loop quantum gravity approach [19,20].
To confront f (R)-gravity with observations in Solar System, one can get constraints on the theories from different measurements, such as the EötWash experiment [21], the geodesic precession of gyroscopes measured by Gravity Probe B [22] and the precession of the binary pulsars PSR J0737-3039 [23]. At cosmological scales one would expect to employ f (R) theories to account for the problems of "dark energy" [24][25][26] and "dark matter" [27][28][29] needed in GR. If f (R) gravity could account for dark matter, besides matching the rotation curves of galaxy clusters, it should also match the measurements on gravitational lensing effect [30]. However, it is shown in [31] that, in a second-order post-Newtonian approximation, a rather general class of f (R) theories is indistinguishable from GR in gravitational lensing effect. Nevertheless, we will show in this paper that a class of f (R) theories is indeed distinguishable from GR in gravitational lensing effect in a third-order post-Newtonian approximation. However, the possibility to account for the dark matter problem with f (R) theory in lensing effect is highly suppressed due to this tiny third-order difference.
This paper is organized as follows. In Sect. 2, we briefly review the field equations of metric f (R)-gravity. In Sect. 3 the post-Newtonian approximation of a class of f (R) theories is formulated to the desired order. In Sect. 4 we introduce the gravitational lensing effect in metric theories of gravity and show how the f (R) gravity can be distinguishable from GR at third-order post-Newtonian approximation. The dif-ference of the lensing reflection indices is discussed in an example. Finally, conclusions and remarks are given in Sect. 5. Throughout the paper, the metric tensor g μν takes the signature (−, +, +, +).

Field equations of f (R) theory
In metric f (R) theories of gravity, the action of gravity coupled to matter fields is given by where g is the determinant of the metric tensor g μν , χ = 8π G/c 4 with G and c being the Newtonian gravitational constant and the vacuum speed of light, respectively, R = g μν R μν is the Ricci scalar, f (R) is a nonlinear function and S M is the standard matter action. The variation of action (2.1) with respect to the metric g μν yields the Euler-Lagrange equations where ∇ μ is the covariant derivative for g μν , g := ∇ μ ∇ μ , and T μν = (−2c/ √ −g)(δS M /δg μν ) is the energy momentum tensor of matter. Taking the trace of Eq. (2.2) we can get where T is the trace of T μν . Using Eq. (2.3), we can rewrite Eq. (2.2) as

Post-Newtonian expansion
The matter constituents in the universe are usually well approximated by a perfect fluid with mass density ρ and pressure p [32]. Hence we assume that the Newtonian potential U of the mass distribution, the typical velocities v and the pressure of the fluid obey such approximation respectively.
In the post-Newtonian approximation, we can further expand the dynamical variables in the field equations perturbatively in powers of 1/c, since we have the following order relation [33][34][35]: where is the ratio of the energy density to the rest-mass density.
We consider the case that the gravitational field is weak and assume that in absence of a gravitational field the background space-time is flat [30]. We also assume that f (0) = 0, which neglects the contribution of a possible cosmological constant and excludes some form of f (R) theories, e.g, f (R) = 1/R. Note that actually the contribution of a possible cosmological constant can be equivalently substituted by the corresponding contribution of an energy-momentum tensor. Moreover, the f (R) form which is unable to get weak field solution is useless here. In a weak field regime the metric tensor can be expanded about the Minkowski metric η μν in its Lorentzian coordinate system as where |h μν | 1. Up to third order the components of the metric tensor can be written as [31,33,36]: where the left upper index (n) means the order O(n). Using Eq. (3.3) we can get the components of the Ricci tensor: Assuming f (R) to be analytic at R = 0, to the second order the Ricci scalar and thus f (R) and f (R) read, respectively, To the leading order the components of the energy-momentum tensor of matter fields read 3) yields at the second order (2) R = −χ (2) T 00 / f (0), which is consistent with the equation of GR at the same order. Thus in this approach, GR is nothing else but the first term of the Taylor expansion of a more general f (R) theory. As one can see from the action (2.1), we ask f (R) to carry the same dimension as R's. Thus both f (R) and the term f (0) (2) R in Eq. (3.5) are dimensionless. Since the term f (0) (2) R is required to be of order O(2), the expansion of Eq. (3.5) will break down if (2) To derive neat equations in the post-Newtonian approximation, we impose the gauge conditions [31,36] With the gauge conditions, we get from Eqs. (2.3) and (2.4) For the sake of physics and simplicity, we consider the case of f (0) > 0 and f (0) > 0 and define Note that in this case the constant 1/ f (0) can be absorbed into the gravitational constant G if necessary. Then from Eq. (3.9) we can get [34,36] ( with the Yukawa potential (3.14) Note that the information of a specific form of f (R) gravity is encoded in the parameter α in the potential V , and we only consider the solution with α > 0. It should be noticed that, for the other solution with α < 0, the potential V would tend to be divergent at infinity.
It is easy to show by using Eq. (3.9) that Eq. (3.10) can be written as Thus we get Using Eq. (3.13), the solution of Eq. (3.16) can be given by where the Newtonian potential U reads It is obvious that the Newtonian potential U remains unchanged for different forms of f (R) gravity. Similarly, the solution of Eq. (3.11) reads From Eqs. (3.10) and (3.11), we have Using Eqs. (3.17) and (3.19), it is easy to see that Thus from Eq. (3.13) we have Then Eq. (3.17) can be written as where the potential ψ is defined as such that ∇ 2 ψ = −2U . Hence Eq. (3.12) can be written as The solution of Eq. (3.25) reads (3.27) So up to third-order post-Newtonian approximation, the final form of the metric components reads (3.28) In contrast, the metric components to the same order approximation in GR reads [33] g 00 = −1 + 2 c 2 U, (3.29) Hence the difference between the f (R) gravity and GR comes from the Yukawa-like potential V and Z ,0i . In the limit f (0) → 0, we get α → ∞ and V → 0. Then the solution (3.28) of f (R) gravity goes back to (3.29) of GR. On the other hand, it is straightforward to see that, in the limit f (0) → ∞, we have α → 0 and hence get the most obvious departure of f (R) gravity from GR.

Gravitational lensing
A gravitational lens refers to a distribution of plates (such as a cluster of galaxies) between a distant source (a background galaxy) and an observer that is capable of bending the light from the source, as it travels towards the observer. The lensing effect can magnify and distort the image of the background source [37]. According to Fermat's principle, the world line of a light should extremize its arrival time T with respect to an observer under the variation of γ . In metric theories of gravity, this principle implies that the world line of a light coincides with a null geodesic in the spacetime.
In the Lorentzian coordinate system of the flat background spacetime, let dl 2 = δ i j dx i dx j be the spatial Euclidean line element. Up to a constant, the travel time of light on a null geodesic γ is given by where we defined the effective refraction index of light as Then Eq. (4.1) takes the form similar to that of the propagation of a light through a medium in Newtonian space and time.

Second-order expansion
At the second-order post-Newtonian approximation, the only nonzero perturbative metric components are h 00 and h i j in Eqs. (3.17) and (3.19). For a null geodesic, we have [31,38,39] Hence we can get where h is defined such that Eq. (3.19) can be written as (2) h i j = (2) hδ i j . Using Eq. (4.2), we can obtain the effective index of refraction as which at the second order reads which illustrates that at second order, the effective refraction index n is only determined by the Newtonian potential U . As shown in Eqs. (3.28) and (3.29), the difference between f (R) gravity and GR comes from the potential V rather than U . Hence one cannot distinguish f (R) theories from GR at second-order approximation by the gravitational lensing effect [31].

Third-order expansion
We now consider the third-order post-Newtonian approximation which is needed in dealing with light rays in spacetime [33]. At the third-order expansion, the line element of the metric can be written as For a null geodesic, by using Eq. (4.2) we have Thus we obtain where In the third-order approximation, Eq. (4.11) can be expressed as F( (2) h 00 , (2) h, (3) where (4) E( (2) h 00 , (2) h) represents the expansion terms at fourth order. It is obvious that the − sign in front of the function F in Eq. (4.10) should be neglected, since otherwise the refraction index n would become negative. Hence Eq.
(4.10) becomes It is easy to see that Eq. (4.13) can be solved: where u i = dx i /dt is the components of the coordinate speed of light, and n 2 := 1 + 2U (x, t)/c 2 is the refraction index at second order. Therefore, in a third-order post-Newtonian approximation, the effective refraction index of light is obviously dependent on the third-order metric components h 0i . From Eqs. (3.27), (3.28) and (3.29) one can see that, in contrast to the case of GR, in f (R) gravity h 0i is effected also by the Yukawa potential V . Hence, f (R) theories are in principle distinguishable from GR by gravitational lensing effect in a third-order post-Newtonian approximation.

Differences: an example
Although the difference of lensing effect between f (R) gravity and GR is encoded in the third-order terms, it is still unclear whether the difference can actually be detected at this order and in which case the departure become most obvious.
To answer these questions, we first recall from Eq. (4.14) that the difference at third-order effect is contained in the difference of the metric components h 0 j between f (R) gravity and GR, which reads A straightforward calculation leads to are two monotone decreasing functions. Hence the lensing refraction indices of f (R) gravity and GR will have the biggest departure in the limit of α → 0, which reads One may notice that this is nothing else but the potential ψ ,0i appearing in the third-order post-Newtonian approximation of GR. Thus, in the case of the highest departure, the difference of the lensing refraction indices is at the same third order of GR. By noticing that the functions g and h satisfy the relation where d = |x − x |, Eq. (4.16) can be further simplified. In terms of cylindrical coordinates {r, θ, z}, Eq. (4.16) can be written as where we used the identity d θ |x − x | = d θ r 2 + r 2 − 2rr cos(θ ) + (z − z ) 2 = rr sin(θ )/|x − x | and θ = 0.
In an axially symmetric spacetime, it is reasonable to consider the case that the velocities v of the gravitational sources are all tangent to the r -θ plane. Then from Eq. (4.20) one gets Therefore, in an axially symmetric spacetime region, which coincides in most cases with those of galaxies and compact objects, one cannot distinguish f (R) theories from GR by the lensing correction to the third-order term. This result suggests that there are few opportunities to distinguish f (R) theory from GR even in a third-order post-Newtonian approximation. However, the difference n will not vanish in a nonaxially symmetrical spacetime region, thus one could detect the difference in principle. Since the potential appearing in Eq. (4.20) without function g are third-order post-Newtonian terms, the order of n is determined by the order of the functions g(αd) (or g(αd) and h(αd), with the same order as shown in Fig. 1. It is shown that, for α ≤ 1 d , n will be around 10 −1 times the 3rd GR terms, and thus keeps the same third order. However, for α ≥ 10 d , n will be less than 10 −2 times the third-order terms, and hence is indistinguishable with the fourth-order term. This estimation would approach the exact result for the spacetime region far away from the matter center. Then the functions g and h could be approximated by the values g(αd 0 ) and h(αd 0 ), where d 0 is the distance of the position to the matter center. Thus one could write the first-order approximation of (4.16) as Fig. 1 The evaluation of the functions g(αd) and h(αd) respect to αd in which the whole integration is approximated by the integration of the region around the center where most part of the matter locate. Therefore, in highly non-axially symmetric spacetime region, it is possible to distinguish the lensing of GR from those of the f (R) theories satisfying α ≤ 1 d 0 . It also requires that the measurement can approach the 10 −1 precision of the third-order effect.

Concluding remarks
In this paper, the post-Newtonian approximation of a general class of f (R) theories is formulated up to third order. In the third-order expansion, the metric components contain not only the Newtonian potential U but also the Yukawa potential V together with the third order potentials. Note that f (R) theories can be transformed into generalized Brans-Dicke theories by suitable conformal transformations. Since the post-Newtonian formulation of Brans-Dicke gravity has been well studied [33], one can check the consistency of the post-Newtonian formulations between the two kinds of theories. It turns out that, in the limit of α → 0, our result (3.28) of f (R) gravity coincides with the result of Brans-Dicke gravity given in [33]. The proof will be presented in Appendix A.
In our post-Newtonian formulation, the information of a specific form of f (R) theories is contained in the Yukawa potential. While the Yukawa potential does not show in the second-order expression of the effective refraction index n of light, it does appear in the third-order expression of n. Therefore in principle we could distinguish f (R) gravity and GR. Moreover, detailed analysis shows that a series of f (R) forms, more specific, whose parameter 1/α ∼ f (0)/ f (0) is larger than the distance to the massive center, are distinguishable from GR by the gravitational lensing effect in a third-order post-Newtonian approximation. It should be noted that the conclusion that f (R) theories can lead to the gravitational lensing effect different from that of GR can also be obtained by the approach of Minkowski functionals [40]. Moreover, it should be pointed out that the third-order perturbations are also distinguishing f (R) theories from GR in view of the Birkhoff theorem [41]. However, it is shown in this paper that, in the axially symmetrical spacetime region, the gap term between these two kinds of theories vanishes and hence they are indistinguishable at third order.
One of the motivations for developing modified gravity theories is to account for the observed mass profiles in galaxies as well as clusters of galaxies without the inclusion of dark matter. The existence of dark matter in GR is confirmed by the observational data not only from the dynamical analysis, such as rotation curves in spiral galaxies [42] and velocity dispersions in early-type systems [43,44], but also from gravitational lensing observations [45,46]. Observations indicate that we need to take into account almost the same large amount of dark matter to explain the gravitational lensing effect as that for the dynamical data like the velocity dispersion or the temperature profile of the X-ray emitting intracluster medium [47,48] in galaxy clusters or spiral galaxies. Up to now, certain f (R) theories are tested by the dynamical data in galaxy clusters and spiral galaxies [27][28][29][49][50][51]. However, concerning the gravitational lensing observations, our results here show a disfavor of the attempts in this direction. For any f (R) form which could be weakly expanded, the lensing effect correction due to the f (R) will be at most of third order, which is at most 10 −2 times of the leading order, i.e., the second-order post-Newtonian effect. Moreover, the fact that in axially symmetrical spacetime region there is no difference in lensing effect between these two kinds of theories strongly indicates that most of the lensing observations will not show the difference even at third order. Thus it is impossible to explain the lensing observations in the pure f (R) theories that we are considering without any dark matter involved. Actually there is already some evidence implying that f (R) theories without dark matter behave badly for galaxy clusters [52].
It is still possible to determine the parameter α 2 := f (0)/ (3 f (0)) though the precise observational results in non-axially symmetric system. Thus in near future, precise observations of lensing effect would be useful to distinguish certain f (R) theories from GR. It should be remarked that our result is only valid for the f (R) forms which could be weakly expanded. It is interesting to further study whether the dark matter content can be replaced by other unexpandable f (R) theories or other kinds of modified gravity.
ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .