Consider Einstein-Hilbert action for gravitational field,
S=∫d4x√−g{12κR−1κΛ+LM} where R=Rμνgμν is the curvature scalar, Rμν=Rλμλν is the Racci tensor, κ=8πG/c2.
Einstein-Hilbert action is the most simple Lorentzian invariant action that encodes space-time curvature and gives correct gravity (Newtonian gravity) in weak field limit. A straightforward generalization of Einstein-Hilbert action is the f(R) gravity: R→f(R)=a0+a1R+a2R2+⋯. From the aesthetic point of view, one may argue the higher order terms shouldn't be there. But there terms may arise for a good reason. We know that at high energy ∼TPl, general relativity (GR) will be superseded by quantum gravity. This in fact may even happen before Planck scale. General relativity is merely a low energy effective theory. Effective theories are not always neat. In fact, it is well known that quantized general relativity is non-renormalizable. Therefore, quantum fluctuation would bring infinite terms into the effective Lagrangian.
A natural question is, if these term has any observable effect. Let's restrict ourselves to classical theory only. The question is, how well has general relativity been tested [1]? f(R) gravity and many other competing theories of gravity shares the same foundations with general relativity, except for the amount of gravity produced by the same energy. The relevant test includes deflection of light ray and time delay (Fig. 1), Mercury perihelion shift and spin precession (Table 4), change of Newton's constant (Table 5), mass of graviton (or propagation distance of gravity) etc.
Note that most of these tests are conducted with in solar system. They cannot exclude gravity theory with only large scale effects. It has been known [2-3] R2 and/or RμνRμν term can produce a massive scalar field in addition to graviton. This field couples with matter through Yukawa potential. Its mass can be chosen such than too heavy to propagate a observable distance, but too light ≪MPl to modified gravity. In other words, it is a massive, very weakly interaction particle. Naturally, it is a dark matter candidate, since it only manifests itself in large scales.
Now let's see how it happens. Consider Rn term: Sn=1κ∫d4x√−gRn Now let's take functional derivative with respect to gμν. δ√−gRnδgμν=√−gRn−1(nδRδgμν+R√−gδ√−gδgμν) Now, according to Jacobi's formular, δg=ggμνδgνμ. So 1√−gδ√−gδgμν=−12gμν. δRδgμν=Rμν+gαβ(δΓλαβ;λδgμν−δΓλαλ;βδgμν) Note that δRλαλβ=δΓλαβ;λ−δΓλαλ;β. δΓ's are in fact tensors: δΓλαβ=12gλρ{δgρα;β+δgρβ;α−δgαβ;ρ}.
Hence δRδgμν=Rμν+δgαβ;α;βδgμν−δgα;βα;βδgμν
Therefore, the action δSn=1κ∫d4x√−gδgμνRn−1{nRμν−12gμνR+n−1R(R;μ;ν−gμνR;α;α)+(n−1)(n−2)R2(R;μR;ν−gμνR;αR;α)} yields the equation of motion: ∑nanRn−1{nRμν−12gμνR+n−1R(R;μ;ν−gμνR;α;α)+(n−1)(n−2)R2(R;μR;ν−gμνR;αR;α)}=κ2Tμν
For n=1,a1=12 it reduces to GR. One of the most distinguish feature of this formulation is for n=1 (GR), two extra terms vanish. If we keep terms up to n=2, the equation of motion simplifies as:(1+2a2R)Gμν+a2gμνR2+2a2(R;μ;ν−gμνR;α;α)=κTμν
Now, let's identify the particle contents. Define hμν=gμν−ημν, where ημν is the Minkowski metric. In weak field limit, we can simply take hμν→δgμν and gμν→ημν. But graviton is usually defined as the inverse trace ˉhμν=hμν−12ημνh, where h=hμμ. Note that, Gμν=−12◻ˉhμν and R=ˉhμν,μν−12◻h. The lowest order terms are: S=1κ∫d4x{−14ˉhμν◻ˉhμν+18h(2a2◻2+◻)h−a22ˉhμν∂μ∂ν◻h}
As we can see, there is an additional field besides the conventional graviton field, with mass m2h=−12a2. This field interacts with matter via Yukawa coupling: √−g≃1+δ√−g=1+12h. Therefore, it propagates e−r/mh. Taking R≪m2h≪MPl, R2 gravity produces a massive weakly interacting particle.
[1]: Clifford M. Will, The Confrontation between General Relativity and Experiment , Living Rev. Relativity, 9, (2006), 3, http://www.livingreviews.org/lrr-2006-3
[2]: Jose A. R. Cembranos, Dark Matter from R2 Gravity, Phys. Rev. Lett. 102, 141301 (2009)
[3]: S. Stelle, Gen. Relativ. Gravit. 9, 353 (1978)
.
Einstein-Hilbert action is the most simple Lorentzian invariant action that encodes space-time curvature and gives correct gravity (Newtonian gravity) in weak field limit. A straightforward generalization of Einstein-Hilbert action is the f(R) gravity: R→f(R)=a0+a1R+a2R2+⋯. From the aesthetic point of view, one may argue the higher order terms shouldn't be there. But there terms may arise for a good reason. We know that at high energy ∼TPl, general relativity (GR) will be superseded by quantum gravity. This in fact may even happen before Planck scale. General relativity is merely a low energy effective theory. Effective theories are not always neat. In fact, it is well known that quantized general relativity is non-renormalizable. Therefore, quantum fluctuation would bring infinite terms into the effective Lagrangian.
A natural question is, if these term has any observable effect. Let's restrict ourselves to classical theory only. The question is, how well has general relativity been tested [1]? f(R) gravity and many other competing theories of gravity shares the same foundations with general relativity, except for the amount of gravity produced by the same energy. The relevant test includes deflection of light ray and time delay (Fig. 1), Mercury perihelion shift and spin precession (Table 4), change of Newton's constant (Table 5), mass of graviton (or propagation distance of gravity) etc.
![]() |
Fig. 1: test of deflection of light, δθ=1+γ24GMd and test of time delay experiment δt=2(1+γ)GMln((r⊕+x⊕⋅n)(re−xe⋅n)d2). General relativity predicts γ=1. |

![]() |
Table 5: test of change of Newton constant. |
Now let's see how it happens. Consider Rn term: Sn=1κ∫d4x√−gRn Now let's take functional derivative with respect to gμν. δ√−gRnδgμν=√−gRn−1(nδRδgμν+R√−gδ√−gδgμν) Now, according to Jacobi's formular, δg=ggμνδgνμ. So 1√−gδ√−gδgμν=−12gμν. δRδgμν=Rμν+gαβ(δΓλαβ;λδgμν−δΓλαλ;βδgμν) Note that δRλαλβ=δΓλαβ;λ−δΓλαλ;β. δΓ's are in fact tensors: δΓλαβ=12gλρ{δgρα;β+δgρβ;α−δgαβ;ρ}.
Hence δRδgμν=Rμν+δgαβ;α;βδgμν−δgα;βα;βδgμν
Therefore, the action δSn=1κ∫d4x√−gδgμνRn−1{nRμν−12gμνR+n−1R(R;μ;ν−gμνR;α;α)+(n−1)(n−2)R2(R;μR;ν−gμνR;αR;α)} yields the equation of motion: ∑nanRn−1{nRμν−12gμνR+n−1R(R;μ;ν−gμνR;α;α)+(n−1)(n−2)R2(R;μR;ν−gμνR;αR;α)}=κ2Tμν
For n=1,a1=12 it reduces to GR. One of the most distinguish feature of this formulation is for n=1 (GR), two extra terms vanish. If we keep terms up to n=2, the equation of motion simplifies as:(1+2a2R)Gμν+a2gμνR2+2a2(R;μ;ν−gμνR;α;α)=κTμν
Now, let's identify the particle contents. Define hμν=gμν−ημν, where ημν is the Minkowski metric. In weak field limit, we can simply take hμν→δgμν and gμν→ημν. But graviton is usually defined as the inverse trace ˉhμν=hμν−12ημνh, where h=hμμ. Note that, Gμν=−12◻ˉhμν and R=ˉhμν,μν−12◻h. The lowest order terms are: S=1κ∫d4x{−14ˉhμν◻ˉhμν+18h(2a2◻2+◻)h−a22ˉhμν∂μ∂ν◻h}
As we can see, there is an additional field besides the conventional graviton field, with mass m2h=−12a2. This field interacts with matter via Yukawa coupling: √−g≃1+δ√−g=1+12h. Therefore, it propagates e−r/mh. Taking R≪m2h≪MPl, R2 gravity produces a massive weakly interacting particle.
[1]: Clifford M. Will, The Confrontation between General Relativity and Experiment , Living Rev. Relativity, 9, (2006), 3, http://www.livingreviews.org/lrr-2006-3
[2]: Jose A. R. Cembranos, Dark Matter from R2 Gravity, Phys. Rev. Lett. 102, 141301 (2009)
[3]: S. Stelle, Gen. Relativ. Gravit. 9, 353 (1978)
.