---
title: Higher-Order Curvature-Scalar Gravity
url: https://www.emergentmind.com/topics/higher-order-curvature-scalar-gravity
type: topic
---

# Higher-Order Curvature-Scalar Gravity

In the literature surveyed here, higher-order curvature-scalar gravity functions as an umbrella description for gravitational theories in which the Einstein-Hilbert action is extended by nonlinear curvature invariants and, in many formulations, by explicit scalar fields or scalar degrees of freedom generated through duality transformations. Typical ingredients include \(R^2\), \(R^3\), \(R\square R\), \(RR_{\mu\nu}R^{\mu\nu}\), Gauss-Bonnet couplings, and general \(\mathcal{L}(\text{Riemann})\) contractions; their applications range from inflation and nonsingular bounces to black holes, holography, braneworld localization, and discrete quantum-gravity constructions [1610.08519] [2510.17965] [2311.12104].

## 1. Scope and representative structures

A broad formulation is provided by general \(\mathcal{L}(\text{Riemann})\) theories in arbitrary dimension,
\[
S=\int_{\mathcal{M}} d^D x \sqrt{|g|}\,\left[\mathcal{L}(R_{\mu\nu\rho\sigma},g^{\alpha\beta})+L_{\rm matter}\right],
\]
whose field equations can be written in terms of
\[
\mathcal{E}_{\mu\nu}\equiv P_{\mu}^{\ \sigma\rho\lambda} R_{\nu\sigma\rho\lambda} - \frac{1}{2}g_{\mu\nu}\mathcal{L} - 2\nabla^{\alpha}\nabla^{\beta} P_{\mu\alpha\beta\nu} = \frac{1}{2} T_{\mu\nu},
\]
with \(P^{\mu\nu\rho\sigma}\equiv \partial \mathcal{L}/\partial R_{\mu\nu\rho\sigma}\). Bueno, Cano, Min, and Visser showed that such theories can be classified on maximally symmetric backgrounds by four constants \((a,b,c,e)\), yielding a spectrum that may contain a massless graviton, a ghost-like massive spin-2 mode, and a scalar mode. The corresponding masses are
\[
\kappa_{\rm eff} = \frac{1}{4e-8\Lambda(D-3)a}, \qquad
m_s^2 = \frac{e(D-2) - 4\Lambda [a+b D(D-1) + c(D-1)]}{2a + Dc + 4b(D-1)},
\]
\[
m_g^2 = \frac{-e + 2\Lambda(D-3)a}{2a+c},
\]
so the absence of the massive spin-2 ghost requires \(2a+c=0\), while removal of the scalar requires \(2a + Dc + 4b(D-1)=0\) [1610.08519].

Within this broad class, several concrete higher-order curvature-scalar frameworks recur across the literature.

| Framework | Representative action | Scalar-sector role |
|---|---|---|
| General \(\mathcal{L}(\text{Riemann})\) gravity | \(\int \sqrt{|g|}\,\mathcal{L}(R_{\mu\nu\rho\sigma},g^{\alpha\beta})\) | May propagate a scalar mode depending on \((a,b,c,e)\) |
| \(\mathcal{F}(R,T)\) bounce model | \(\mathcal{F}(R,T)=R+\lambda_1 R^2+2\lambda_2 T\) | Effective matter coupling through \(T\) [2206.12423] |
| \(f(R)\)-dGB theory | \(\frac{1}{16\pi}\!\int\! \sqrt{-g}\{f(R)-(\nabla_\mu\phi)(\nabla^\mu\phi)+2\eta(\phi)\mathcal{G}\}\) | Two nonminimally coupled scalar fields after scalarization [2510.17965] |
| HOCG black-hole framework | \(\int d^4x \sqrt{-g}\, f(R,G,\phi)+S_{\rm matter}\) | Scalar coupled to higher-curvature invariants [2509.11985] |

The scalar sector is not always inserted by hand. In some theories it is explicit, as in \(f(R,\mathcal G,\phi)\) models or Einstein-dilaton-Gauss-Bonnet-type constructions; in others it emerges through conformal or Legendre transformations of higher-curvature terms. This suggests that “curvature-scalar” is best understood as a structural relation rather than a single Lagrangian ansatz.

## 2. Scalarization, dual frames, and formulation dependence

For \(f(R)\) theories, the scalar-tensor correspondence is mathematically explicit. In the regular-spacetime analysis of \(f(R)\) gravity, a conformal transformation
\[
\mathbf{g}_{ab}=\Omega^2 g_{ab}, \qquad \Omega^{3}=f'(R),
\]
together with
\[
\kappa_5 \phi = \frac{2}{\sqrt{3}}\ln f'(R),
\]
maps the Jordan-frame action to Einstein gravity plus a scalar field with potential
\[
V(\phi)=\frac{f'(R)R-f(R)}{2\kappa_5^2 [f'(R)]^{5/3}}.
\]
The same work shows that the equivalence extends from the action to the field equations, and even to their projection on a lower-dimensional hypersurface, provided the conformal transformation is regular [1604.05301].

A five-dimensional braneworld realization of this correspondence uses \(F(R)=R+\alpha R^2\). After introducing an auxiliary field and performing the conformal transformation \(G_{MN}\to \tilde G_{MN}=e^\sigma G_{MN}\), the theory becomes Einstein gravity plus a scalar field \(\Phi\). In that setting, fermion localization properties are identical in the \(F(R)\) and scalar-tensor descriptions: increasing \(\alpha\) localizes massless chiral modes toward the TeV brane, while massive KK fermions localize toward the Planck brane, and the radion-fermion coupling for massive KK modes is strongly suppressed on the visible brane [1707.06532].

Formulation dependence is especially sharp in Palatini and metric-affine theories. In the Palatini treatment of Higgs inflation with higher-order gravity, the Jordan-frame action
\[
S = \int d^4 x \sqrt{-g} \left\{ \frac12 \left[ \left(1+\xi h^2 \right) R + \beta R^2 + \alpha R_{(\mu\nu)} R^{(\mu\nu)} \right] - \frac12 g^{\mu\nu} \partial_\mu h \partial_\nu h - V(h) \right\}
\]
is restricted to terms built from the symmetric part of the Ricci tensor. In this sector there are no new gravitational degrees of freedom; the scalar perturbation spectrum is unchanged, while the tensor perturbation spectrum is suppressed by the higher-order curvature couplings. The Einstein-frame potential becomes
\[
\tilde{U}(\phi)=\frac{U(\chi(\phi))}{1 + (2\alpha + 8\beta) U(\chi(\phi))},
\]
and the tensor-to-scalar ratio is
\[
r = \frac{16\epsilon_U}{1 + (2\alpha + 8\beta) U}
\]
[2106.09438].

That simplification is not generic. In metric-affine Ricci-based gravity, the claim that higher-order curvature theories are generically ghost-free is explicitly rejected. When projective symmetry is broken and the action depends on the full Ricci tensor rather than only \(R_{(\mu\nu)}\), the non-projectively invariant sector propagates ghost-like degrees of freedom. The pathology can be avoided either by imposing projective symmetry,
\[
\Gamma^\lambda_{\mu\nu}\rightarrow \Gamma^\lambda_{\mu\nu} + \delta^\lambda_\mu \xi_\nu,
\]
or by enforcing a torsion-free constraint \(T^\lambda_{[\mu\nu]}=0\) [1901.08988]. A common misconception is therefore that “Palatini” or “metric-affine” automatically removes higher-derivative pathologies; the cited results show that this is true only in specially constrained sectors.

## 3. Cosmology: second-order FLRW sectors, inflation, and bounces

A major line of work isolates higher-curvature theories whose FLRW equations remain second order. “Cosmological Gravities” are defined as higher-curvature theories whose equations of motion for FLRW configurations are of second order in time derivatives for the scale factor, just as in Einstein gravity. An all-orders family in \(D\geq 3\) is
\[
\mathcal{L}_{\text{cosmo}} = -2\Lambda + R + \sum_{n=2}^\infty \mu_n \ell^{2n-2} \mathcal{C}^{(n)},
\]
with a unique nontrivial equivalence class at each curvature order. In these theories, the linearized equations for scalar cosmological perturbations contain no more than two time derivatives, and explicit four-dimensional examples were given up to fifth order in the curvature [2311.12104].

An earlier FLRW construction based on Lovelock tensors reached a related conclusion from a different direction. By constructing higher-order terms linear in \(\ddot a\), the resulting Lagrangians can be reduced by integration by parts to dependence only on \(a\) and \(\dot a\), so no scalar degrees of freedom appear in the cosmological model. The resulting dynamics takes the form of \(f(H^2)\) cosmology, with
\[
f(H^2)=3p,
\]
and inflationary behavior of the scale factor can arise without introducing ad hoc inflaton fields [1202.2428].

Inflationary model building in higher-order curvature-scalar gravity has proceeded along several distinct routes. One route extends Starobinsky inflation by including all scalar-curvature corrections up to second order:
\[
S = \frac{M_{\mathrm{Pl}}^2}{2} \int d^4x \sqrt{-g}\left( R + \frac{1}{2\kappa_0} R^2 + \frac{\alpha_0}{3\kappa_0^2} R^3 - \frac{\beta_0}{2\kappa_0^2} R \Box R \right).
\]
In the Einstein frame this becomes a two-field model. Phase-space analysis exhibits a slow-roll attractor and a stable spiral at the origin for \(0\leq \beta_0\leq 3/4\), scalar and isocurvature perturbations can be separated into \(Q_\sigma\) and \(Q_s\), and observational agreement with Planck 2018, BICEP3/Keck, and BAO is obtained for approximately \(\alpha_0 \lesssim 10^{-4}\) and \(\beta_0\lesssim 10^{-2}\) [2207.02103].

A dimension-six extension focused on
\[
R+R^2+RR_{\mu\nu}R^{\mu\nu}
\]
derives a four-dimensional autonomous system in the Jordan frame and linearizes in the cubic parameter \(\gamma_0\). The slow-roll potential becomes
\[
V(\varphi) \approx \frac{1}{6} \left( 1 - e^{-\varphi} - \frac{13}{144}\gamma_0 e^{\varphi} \right)^2,
\]
with inflationary observables
\[
n_s \approx 1 - \frac{2}{N} \left( 1 + \frac{52}{243}\gamma_0 N^2 \right), \qquad
r \approx \frac{12}{N^2 \left( 1 - \frac{52}{243}\gamma_0 N^2 \right)}.
\]
Compatibility with Planck, BICEP/Keck, and BAO requires \(|\gamma_0|\lesssim 10^{-3}\), while negative \(\gamma_0\) restores compatibility with recent ACT, Planck, and DESI results [2509.09167].

A more radical possibility is inflation in a higher-curvature dominated regime. In Geometric Inflation models,
\[
3F(H)=\frac{1}{M_{\rm Pl}^2}\rho,\qquad -H'\frac{dF(H)}{dH}=\frac{1}{M_{\rm Pl}^2}(\rho+P),
\]
with
\[
F(H)=H^2+\Lambda_*^{-2}\sum_{n=3}^{\infty}(-1)^n c_n (\Lambda_* H)^{2n}.
\]
If \(F(H)\) has a GR regime, a transition or flat regime, and a steep inflating regime, then a minimally coupled scalar with quadratic potential can drive more than 100 e-folds with \(\phi,\Delta\phi\ll M_{\rm Pl}\) in explicit examples [2007.07651].

Bounce cosmology provides the complementary use of higher-curvature terms. In \(\mathcal{F}(R,T)\) gravity with
\[
\mathcal{F}(R,T)=R+\lambda_1 R^2 + 2\lambda_2 T,
\]
and a power-law Hubble parameter
\[
H(t)=\alpha+\beta t + \gamma t^n,
\]
the conditions \(H(t_b)=0\), \(\dot H(t_b)>0\), and \(a(t_b)>0\) can be satisfied, the equation-of-state parameter crosses the phantom divide \(\omega=-1\), and the bounce avoids the cosmological initial singularity [2206.12423]. A related analysis in quadratic curvature gravity minimally coupled to a scalar field distinguishes two energy-momentum tensors. If only the scalar field is counted, NEC, WEC, and DEC remain satisfied while SEC is violated during the bounce. If an effective tensor is constructed to absorb higher-curvature corrections, all four energy conditions are violated near the bounce, emphasizing that the non-Einsteinian behavior resides in the gravitational sector rather than in exotic matter [2603.23564].

## 4. Spherically symmetric solutions and black-hole phenomenology

Higher-order curvature-scalar gravity also generates a large class of static and spherically symmetric solutions. A six-derivative scalar-curvature theory with
\[
S = \frac{M_P^2}{2}\int d^4x\,\sqrt{-g}\left[ R + \frac{1}{2\kappa_0}R^2 + \frac{\alpha_0}{3\kappa_0^2}R^3 - \frac{\beta_0}{2\kappa_0^2}R\square R \right] + \int d^4x\,\sqrt{-g}\mathcal{L}_m
\]
yields weak-field potentials of Yukawa type,
\[
\phi = -\frac{GM}{r}\left(1 + a_+ e^{-r/l_+} + a_- e^{-r/l_-} \right),
\]
with
\[
l_{\pm} = \sqrt{\frac{3}{2\kappa_0}\left( 1 \pm \sqrt{1 - \frac{4\beta_0}{3}} \right)}.
\]
Stability restricts the parameters to
\[
\kappa_0>0,\qquad 0\leq \beta_0 < \frac34,
\]
so only Yukawa-type corrections are physically allowed. Non-Schwarzschild black holes might exist mathematically, but for reasonable values of the parameters their horizon radii are extremely small, \(r_h\lesssim 10^{-31}\,\text{m}\), making macroscopic deviations from Schwarzschild unfeasible [2004.04878].

The strong-field situation is more restrictive in Einsteinian cubic gravity. For the action
\[
\mathcal{S} = \frac{M_\mathrm{pl}^2}{2}\int d^4x\sqrt{-g}\left(R+\alpha \mathcal P\right),
\]
odd-parity perturbations around static and spherically symmetric black holes reveal three propagating degrees of freedom rather than the single odd-parity mode of GR. At least one is always a ghost, and one mode has
\[
c_{\Omega 3}^2 = -\frac12,
\]
so static and spherically symmetric vacuum black holes supported by unsuppressed cubic curvature terms are excluded by ghost and Laplacian instabilities [2305.07217].

A distinct HOCG black-hole solution employs a scalar field coupled to higher-order curvature terms. In one formulation,
\[
S = \int d^D x \sqrt{-g} \left\{ \frac{1}{2\kappa^2} f(R) - \frac{1}{2}\partial_\mu \phi_1 \partial^\mu \phi_1 + V(\phi_1) + H(\phi_1) G \right\},
\]
and the static line element is built from
\[
S(r) = 1 - \frac{c}{r} + \frac{\digamma}{r^2}, \qquad
S_1(r) = 1 - \frac{c}{r} + \frac{c \digamma^{3/2}}{r^4}.
\]
This solution supports up to two horizons \(r_\pm = \frac{c \pm \sqrt{c^2-4\digamma}}{2}\), admits Wald and Barrow entropies, exhibits Davies-type phase transitions where the heat capacity diverges, and has positive quasi-local energy and Gibbs free energy in the thermodynamic analysis. Geodesic deviation yields an explicit stability condition outside the event horizon [2509.02643].

The phenomenology of a related static metric,
\[
A(r,\xi)=1-\frac{2M}{r}+\frac{\xi}{r^2},\qquad
B(r,\xi)=1-\frac{2M}{r}+\frac{2M\xi^{3/2}}{r^4},
\]
has been analyzed in detail. Increasing \(\xi\) decreases the event horizon radius and increases the Cauchy horizon radius; the photon sphere is shifted to
\[
r_{\text{photon}}=\frac12\left(3M+\sqrt{9M^2-8\xi}\right),
\]
the shadow radius decreases, and quasinormal modes for scalar, vector, tensor, and spinorial perturbations become less damped. Weak-field and strong-deflection lensing are both modified, and comparison with EHT and Solar System data constrains \(\xi\): for M87* one finds \(\xi/M^2\lesssim 0.091\), while the tightest Solar System bound comes from light deflection, \(|\xi|\lesssim 10^{13}\,\mathrm{m}^2\) [2509.11985].

Combining \(f(R)\) gravity with Einstein-dilaton-Gauss-Bonnet interactions produces a bi-scalar theory with second-order field equations,
\[
S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ f(R) - (\nabla_\mu \phi)(\nabla^\mu \phi) + 2\eta(\phi) \mathcal{G} \right\}.
\]
In that framework black holes are modified by \(f(R)\) terms, retain the EdGB features of a minimum mass and multiple branches, and display a nontrivial suppression of Ricci-scalar divergence in the interior. Even so, the singularity structure and elliptic regions inside the horizon remain similar to those of pure EdGB gravity, so adding individual higher-order terms does not resolve the underlying ill-posedness [2510.17965].

## 5. Entropy, holography, and matter localization

Black-hole thermodynamics in higher-curvature gravity departs sharply from the pure Wald picture once dynamical processes are considered. For arbitrary higher-curvature gravity with scalar and electromagnetic fields, the entropy satisfying the linearized second law is
\[
S = \int_\mathcal{S} \sqrt{y}\,\Big[ p_W - 2 P^{ij} C_{ij} + N(\nabla_b \phi) \Big],
\]
where \(p_W\) is the Wald entropy density, \( -2P^{ij}C_{ij}\) is the higher-curvature correction already known from previous work, and \(N(\nabla_b\phi)\) is a new correction sourced by the minimal coupling interaction between gravity and the scalar field. The electromagnetic field does not contribute to this entropy correction at linear order. A standard misconception is therefore that only non-minimal curvature couplings modify the entropy relevant for the second law; the cited result shows that minimally coupled scalars can also do so [2211.06882].

In holography, three-dimensional new massive gravity,
\[
S = \frac{1}{16\pi G} \int d^3 x \sqrt{-g}\left[ \sigma R + \frac{2}{\ell^2} + \frac{1}{m^2}\left(R_{\mu\nu}R^{\mu\nu} - \frac{3}{8}R^2\right) \right],
\]
provides an explicit AdS/BCFT example in which higher-curvature terms can be absorbed into effective couplings. The effective Newton constant is
\[
G_{\mathrm{eff}} = \frac{G}{\sigma + \frac{1}{2m^2 L^2}},
\]
the central charge becomes
\[
c = \frac{3L}{2G}\left(\sigma+\frac{1}{2m^2L^2}\right),
\]
and the boundary entropy and entanglement entropy retain the Einstein-gravity form after the substitution \(G\to G_{\mathrm{eff}}\) or, equivalently, \(c\to c_{\rm NMG}\). The holographic \(g\)-theorem continues to hold under the null energy condition [1201.1988].

Matter localization in warped extra dimensions shows a different manifestation of higher-order curvature-scalar effects. In a five-dimensional \(F(R)=R+\alpha R^2\) braneworld, the bulk geometry becomes \(y\)-dependent rather than constant-curvature, massless left- and right-chiral fermion zero modes are driven toward the TeV brane as \(\alpha\) increases, and massive KK modes move toward the Planck brane. The effective radion-fermion coupling for zero modes increases with \(\alpha\), but the ratio \(\lambda_{L,R}/m_{\rm rad}\) decreases, suppressing collider visibility; the same localization pattern persists in the scalar-tensor dual description [1707.06532].

## 6. Structural issues, discrete extensions, and unresolved questions

A recurring theme is that higher-order curvature-scalar gravity is highly non-generic: consistency depends on carefully chosen combinations, symmetry restrictions, or background truncations. In general \(\mathcal{L}(\text{Riemann})\) theories, the massive spin-2 sector is ghost-like unless the coefficients satisfy the Einstein-like constraints, and only specially tuned combinations such as Einsteinian cubic gravity, its quartic extensions, \(f(\text{Lovelock})\), or the new cubic term \(\mathcal Y\) remove the ghost-like spin-2 mode at linear order [1610.08519]. On the cosmological side, second-order FLRW dynamics can be engineered systematically, but this property does not automatically extend to all perturbative or strong-field sectors [2311.12104].

Discrete quantum-gravity work provides an independent route to higher-order invariants. In causal set theory, the Benincasa-Dowker operator
\[
B(\phi(x)) = \frac{1}{\ell^2}\left(\alpha_d \phi(x) + \beta_d \sum_{i=1}^{n_d} C_i^{(d)} \sum_{y\in L_i(x)} \phi(y)\right)
\]
has the continuum limit
\[
\lim_{\ell\to 0} B(\phi(x)) = \left(\Box - \frac12 R(x)\right)\phi(x).
\]
Iterating the operator yields
\[
\lim_{\ell\to 0} B^2(4)=R(x)^2-2\Box R(x),
\]
and, more generally,
\[
\langle B^n(\phi(x))\rangle \to \left(\Box-\frac12 R(x)\right)^n \phi(x)+\mathcal O(\ell^2),
\]
so higher-order Ricci-scalar invariants emerge directly from the discrete causal structure [2301.13525].

The main controversies are therefore not about whether higher-curvature terms can be written down, but about when they define physically viable theories. Metric-affine theories generically remain pathological without projective symmetry or additional constraints [1901.08988]. Unsuppressed cubic corrections can exclude black-hole solutions through ghost and Laplacian instabilities [2305.07217]. In higher-curvature-dominated inflation, perturbation theory and the calculation of cosmological observables remain technically difficult and are explicitly identified as open questions [2007.07651]. In combined \(f(R)\)-Gauss-Bonnet theories, second-order field equations are retained, but singularity resolution and hyperbolicity problems are not achieved by adding finitely many individual higher-order terms [2510.17965].

Taken together, these results indicate that higher-order curvature-scalar gravity is not a single theory but a research program organized around a tension: higher-curvature operators and scalar sectors greatly enlarge the space of cosmological and black-hole solutions, yet only specially structured models avoid ghost modes, gradient instabilities, entropy-law violations, or strong-field ill-posedness. This suggests that the enduring problems in the subject are not the construction of new invariants, but the identification of nonperturbative sectors in which the enlarged dynamics remains mathematically controlled and observationally viable.

Source: https://www.emergentmind.com/topics/higher-order-curvature-scalar-gravity