Skip to content

f(R) Gravity and Its Scalar Representation

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

An f(R)f(R) action changes Einstein gravity by replacing a term linear in the Ricci scalar with a nonlinear function. In the metric variational problem this introduces one scalar degree of freedom: the trace of the metric equation is a differential equation for fR≡df/dRf_R\equiv df/dR. An auxiliary-field transformation makes that scalar explicit, but only on a branch where the relevant Legendre map is regular and invertible. Varying an independent connection instead gives a different theory whose scalar is algebraically, rather than dynamically, determined by matter.

Required background. Gravitational Actions and Boundary Terms supplies the variational boundary conditions, while Curvature supplies the Ricci identities and curvature contractions used below.

Helpful background. Einstein Equations supplies the trace and Bianchi checks; Formulations of General Relativity clarifies metric and independent-connection variations; Quantum Gravity provides the neighboring EFT interpretation without turning this classical equivalence into a microscopic claim.

Use the site signature (−,+,+,+)(-,+,+,+) and curvature convention, and let matter couple minimally to the metric. A useful Dirichlet problem fixes both the induced metric and fRf_R at a non-null boundary, with the scalar-weighted Gibbons–Hawking–York completion. Fixing the metric alone is generally insufficient for this higher-derivative action. With those data, or with compactly supported bulk variations, take

S[g,ψ]=116πG∫Md4x−g f(R)+Sm[g,ψ]+S∂M.S[g,\psi] =\frac{1}{16\pi G}\int_M d^4x\sqrt{-g}\,f(R) +S_m[g,\psi]+S_{\partial M}.

The two identities needed in the bulk variation are

δ−g=−12−g gabδgab,δR=Rabδgab+gab□δgab−∇a∇bδgab.\delta\sqrt{-g}=-\frac12\sqrt{-g}\,g_{ab}\delta g^{ab}, \qquad \delta R=R_{ab}\delta g^{ab} +g_{ab}\Box\delta g^{ab}-\nabla_a\nabla_b\delta g^{ab}.

Integrating the last two terms by parts gives

fRRab−12fgab+(gab□−∇a∇b)fR=8πGTab.f_RR_{ab}-\frac12fg_{ab} +\left(g_{ab}\Box-\nabla_a\nabla_b\right)f_R =8\pi G T_{ab}.

This derivation and its boundary assumptions are developed in Sotiriou and Faraoni 2010, § II.A, equations (5)–(12). Taking the trace exposes the extra equation,

3□fR+fRR−2f=8πGT.3\Box f_R+f_RR-2f=8\pi G T.

For f(R)=R−2Λf(R)=R-2\Lambda, fR=1f_R=1 and the derivative term vanishes, so both equations reduce to their Einstein counterparts. For a genuinely nonlinear ff, the trace evolves the scalar quantity fRf_R. Covariant divergence gives zero by the Bianchi identity and the commutator acting on fRf_R; this is an independent conservation check, not another equation of motion.

Introduce an auxiliary field χ\chi before solving any field equation:

S[g,χ,ψ]=116πG∫d4x−g[f(χ)+f′(χ)(R−χ)]+Sm+S∂M.S[g,\chi,\psi] =\frac{1}{16\pi G}\int d^4x\sqrt{-g} \left[f(\chi)+f'(\chi)(R-\chi)\right]+S_m+S_{\partial M}.

Variation with respect to χ\chi gives

f′′(χ)(R−χ)=0.f''(\chi)(R-\chi)=0.

On a branch with f′′(χ)≠0f''(\chi)\ne0, this fixes χ=R\chi=R. If in addition Φ=f′(χ)\Phi=f'(\chi) is one-to-one on the branch, define the Legendre potential

V(Φ)=χ(Φ)Φ−f ⁣(χ(Φ)).V(\Phi)=\chi(\Phi)\Phi-f\!\left(\chi(\Phi)\right).

The action becomes

S[g,Φ,ψ]=116πG∫d4x−g[ΦR−V(Φ)]+Sm+S∂M,S[g,\Phi,\psi] =\frac{1}{16\pi G}\int d^4x\sqrt{-g} \left[\Phi R-V(\Phi)\right]+S_m+S_{\partial M},

which is the Jordan-frame Brans–Dicke action with ω=0\omega=0. Here the action contains −V-V; the potential called VV on the general scalar–tensor action page appears as −2V-2V and equals half the Legendre potential used here. Its equations are

ΦGab=8πGTab+∇a∇bΦ−gab□Φ−12Vgab,R=V,Φ.\begin{aligned} \Phi G_{ab} &=8\pi G T_{ab}+\nabla_a\nabla_b\Phi-g_{ab}\Box\Phi -\frac12Vg_{ab},\\ R&=V_{,\Phi}. \end{aligned}

Eliminating RR with the metric trace yields

3□Φ+2V−ΦV,Φ=8πGT,3\Box\Phi+2V-\Phi V_{,\Phi}=8\pi G T,

the same differential content as the f(R)f(R) trace. The complete equivalence and the role of f′′≠0f''\ne0 are stated in Sotiriou and Faraoni 2010, § III.A, equations (50)–(63).

An Einstein-frame kinetic term can be obtained locally with gabE=Φgabg^E_{ab}=\Phi g_{ab} and a logarithmic scalar coordinate, but this requires Φ>0\Phi\gt0. Matter then couples to Φ−1gabE\Phi^{-1}g^E_{ab}, so the map is a field redefinition only when that transformed matter action, its clocks, and its sources are retained.

Metric and Palatini theories are not the same

Section titled “Metric and Palatini theories are not the same”

In the Palatini problem, gabg_{ab} and a torsion-free connection Γ^abc\hat\Gamma^a{}_{bc} are varied independently:

S[g,Γ^,ψ]=116πG∫d4x−g f(R)+Sm[g,ψ],S[g,\hat\Gamma,\psi] =\frac{1}{16\pi G}\int d^4x\sqrt{-g}\,f(\mathcal R) +S_m[g,\psi],

where R=gabRab(Γ^)\mathcal R=g^{ab}\mathcal R_{ab}(\hat\Gamma) and the displayed matter action does not depend on Γ^\hat\Gamma. The equations are

fRR(ab)−12fgab=8πGTab,∇^c ⁣(−g fRgab)=0.f_{\mathcal R}\mathcal R_{(ab)}-\frac12fg_{ab}=8\pi G T_{ab}, \qquad \hat\nabla_c\!\left(\sqrt{-g}\,f_{\mathcal R}g^{ab}\right)=0.

The second equation makes Γ^\hat\Gamma the Levi-Civita connection of hab=fRgabh_{ab}=f_{\mathcal R}g_{ab} on a nonsingular branch. Crucially, the trace is algebraic,

fRR−2f=8πGT.f_{\mathcal R}\mathcal R-2f=8\pi G T.

The scalar representation therefore has Brans–Dicke parameter ω=−3/2\omega=-3/2 and no independent scalar wave equation. The derivation and the metric-affine qualification appear in Sotiriou and Faraoni 2010, §§ II.B–II.C and III.B–III.C, equations (14)–(22) and (64)–(70). If matter couples to the independent connection, even this Palatini representation changes; it must not be inferred from the metric theory.

Let f(R)=R+aR2f(R)=R+aR^2. Then

Φ=1+2aχ,V(Φ)=(Φ−1)24a.\Phi=1+2a\chi, \qquad V(\Phi)=\frac{(\Phi-1)^2}{4a}.

In metric f(R)f(R) gravity, linearizing the trace about Minkowski space with T=0T=0 gives

(□−m02) δR=0,m02=16a.(\Box-m_0^2)\,\delta R=0, \qquad m_0^2=\frac{1}{6a}.

Thus a>0a\gt0 gives a non-tachyonic scalar in this background, and a→0a\rightarrow0 sends its mass above any fixed low-energy scale. This is a controlled Einstein limit only for fields and frequencies that remain below that growing mass and on the branch Φ>0\Phi\gt0. At f′′=0f''=0, at a turning point of f′f', or where fRf_R changes sign, the scalar coordinate or Einstein metric is singular; the equivalence cannot be continued by declaration.

This chapter analyzes classical actions, modes, and EFT validity. Empirical likelihoods and numerical bounds remain separate in Tests of Gravity, while claims about a microscopic completion remain in Quantum Gravity.

Call every f(R) formulation one theory. Metric, Palatini, and metric-affine variations have different equations and matter couplings. The choice of independent fields is part of the definition.

Use the scalar map through a singular point. The auxiliary action is globally equivalent only where f′′≠0f''\ne0 and the selected branch of f′f' is invertible; an Einstein-frame map further requires fR>0f_R\gt0.

For f(R)=R+aR2f(R)=R+aR^2, perform the auxiliary-field transformation and derive the linearized scalar mass. Then repeat the trace calculation in the Palatini theory and explain why no Klein–Gordon equation results.

Solution

The auxiliary equation gives χ=R\chi=R for a≠0a\ne0, while Φ=1+2aχ\Phi=1+2a\chi and V=(Φ−1)2/(4a)V=(\Phi-1)^2/(4a). The metric trace is 6a□R−R=06a\Box R-R=0 in vacuum, hence (□−1/(6a))R=0(\Box-1/(6a))R=0. In the Palatini theory the trace instead reads (1+2aR)R−2(R+aR2)=−R=8πGT(1+2a\mathcal R)\mathcal R-2(\mathcal R+a\mathcal R^2) =-\mathcal R=8\pi GT. It fixes R=−8πGT\mathcal R=-8\pi GT algebraically, so there is no independently propagating scalar.

  • Sotiriou, Thomas P., and Valerio Faraoni. 2010. “f(R)f(R) Theories of Gravity.” Reviews of Modern Physics 82: 451–497. DOI.