Papers
Topics
Authors
Recent
Search
2000 character limit reached

Einsteinian Cubic Gravity Overview

Updated 18 March 2026
  • Einsteinian Cubic Gravity is a higher-curvature gravitational theory that extends general relativity with a unique cubic curvature invariant, ensuring only the massless graviton propagates on maximally symmetric backgrounds.
  • It predicts observable deviations in black hole thermodynamics and holography, with modifications in shadow size, ISCO, and angular momentum for static and slowly rotating black holes.
  • The theory also impacts early-universe dynamics with a non-standard inflationary attractor and raises issues of stability due to potential ghost and Laplacian instabilities in non-perturbative regimes.

Einsteinian cubic gravity (ECG) is a higher-curvature extension of general relativity that incorporates a unique cubic-in-curvature invariant in such a way that the linearized spectrum on maximally symmetric backgrounds contains only the massless, transverse graviton. In contrast to Lovelock or other quasi-topological theories, ECG remains nontrivial and non-topological in four dimensions, and exhibits novel dynamics, black-hole physics, and holographic features. This article surveys the structure, solutions, thermodynamics, holography, stability properties, and cosmological dynamics of ECG, with particular emphasis on four-dimensional static and rotating black holes, their observational signatures, and the crucial subtlety of ghost freedom.

1. Gravitational Action and Field Equations

The ECG action in four spacetime dimensions is given by

S=116πGd4xg[R2Λ+λP],S = \frac{1}{16\pi G}\int d^4x \sqrt{-g} \left[ R - 2\Lambda + \lambda \, \mathcal{P} \right],

where RR is the Ricci scalar, Λ\Lambda the cosmological constant (often set to zero in asymptotically flat contexts), and λ\lambda a coupling constant with dimension (length)2^2 or (mass)2^{-2}. The key cubic curvature invariant is

P=12RabcdRcdefRefab+RabcdRcdefRefab12RabcdRacRbd+8RabRbcRca.\mathcal{P} = 12\,R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} + R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} -12\,R_{abcd}R^{ac}R^{bd} +8\,R_{a}{}^{b}R_{b}{}^{c}R_{c}{}^{a}.

The equations of motion are fourth-order PDEs in the metric but, for spherically symmetric ansätze, reduce to a single ordinary differential equation (ODE) for the metric function f(r)f(r): (f1)rλG2[4f3+12f2r24f(f1)fr212ff(f2(f1)r)]=2GM,-(f-1)r - \lambda G^2\left[ 4f'^3 + 12\frac{f'^2}{r} -24f(f-1)\frac{f'}{r^2} -12f f''\big(f' - \frac{2(f-1)}{r}\big)\right] = 2GM, where ()(\,)' denotes RR0 and RR1 is the ADM mass parameter (Bueno et al., 2016). The uniqueness and physicality of ECG stem from the requirement that, upon linearization on maximally symmetric backgrounds, only the Einstein graviton propagates—a constraint that tightly fixes the relative coefficients in RR2 (Bueno et al., 2016).

2. Static, Spherically Symmetric and Rotating Black Holes

For the static spherically symmetric line element,

RR3

ECG admits one-parameter families of asymptotically flat black-hole solutions. The equation for RR4 reduces to Schwarzschild for RR5. For small RR6,

RR7

Near the horizon RR8, RR9, with surface gravity Λ\Lambda0 and Hawking temperature Λ\Lambda1 (Bueno et al., 2016). The Wald entropy is

Λ\Lambda2

In charged cases with Maxwell coupling, the system admits two coexisting black-hole branches above a critical coupling, violating uniqueness yet retaining pathologies such as the disappearance of inner (Cauchy) horizons found in Reissner–Nordström (2002.04071).

The extension to slowly rotating black holes, Λ\Lambda3, leads to coupled ODEs for Λ\Lambda4 and Λ\Lambda5 solvable numerically. ECG corrections shift the horizon angular velocity, ISCO radius, photon sphere, and shadow size, but the deviations become significant only for compact objects with Λ\Lambda6 (Adair et al., 2020).

3. Thermodynamic and Phase Structure

ECG black holes exhibit thermodynamic properties that deviate from the Schwarzschild paradigm. For asymptotically AdS or flat cases, the Hawking temperature and entropy can be computed analytically as functions of Λ\Lambda7 and Λ\Lambda8; the Abbott–Deser mass matches the integration parameter. The first law, Λ\Lambda9, is satisfied for all values of λ\lambda0 and λ\lambda1 (Bueno et al., 2016, Bakhtiarizadeh, 2021). Notably, ECG solutions feature two branches for λ\lambda2: a small λ\lambda3 branch with positive specific heat (thermodynamically stable) and a large λ\lambda4 branch with negative specific heat (unstable), merging at λ\lambda5, signaling a phase transition (Bueno et al., 2016). In AdS, the equation of state is quadratic in λ\lambda6 and exhibits van der Waals–type criticality when λ\lambda7, with the critical ratio λ\lambda8 matching the universal value from classical thermodynamics (Hennigar et al., 2016).

Higher-dimensional extensions require two independent metric functions and admit “super-entropic” black holes that violate the reverse-isoperimetric inequality (Hennigar et al., 2016). The inclusion of gauge fields yields novel features such as the existence of both small and large black holes with the same λ\lambda9, and the complete absence of inner horizons even in over-extremal regimes (2002.04071).

4. Holographic Aspects and Conformal Field Theory Duals

In the AdS context, four-dimensional ECG provides a tractable holographic model for asymptotically locally AdS spacetimes, yielding explicit non-hairy AdS2^20 black holes. The relationship between bulk and boundary data is as follows (Bueno et al., 2018):

  • The graviton two-point function charge is 2^21, with 2^22 determined by 2^23 (here, 2^24 is the dimensionless ECG coupling).
  • The universal entanglement entropy charge 2^25 controls EE across disk regions.
  • Thermal entropy and Rényi entropies of the dual CFT2^26 can be computed using exact or numerically precise bulk black holes.
  • The ratio 2^27 is a non-analytic function of the ECG coupling and is strictly bounded below by 2^28, with positive-energy black holes forbidding any violation of the KSS bound, regardless of the coupling strength.
  • The generalized Gibbons–Hawking–York boundary term enters the Euclidean on-shell action weighted by 2^29, providing a systematic prescription for holographic renormalization and for extracting dual CFT data (Bueno et al., 2018).

5. Stability and Pathology: Ghosts, Laplacian Instabilities, and the EFT Regime

While ECG is constructed to propagate only the massless graviton at the linearized level on maximally symmetric backgrounds, ghost and Laplacian instabilities can arise for non-perturbative values of the coupling. For spherically symmetric black holes with order-unity coupling, odd-parity metric perturbations around the ECG black hole exhibit three dynamical modes (rather than the single pure-spin-2 degree of freedom of GR). One mode always behaves as a ghost (sign mistuned kinetic term), and another acquires a negative sound speed squared for high-2^{-2}0 angular modes, 2^{-2}1, resulting in Laplacian instability on arbitrarily short timescales (Felice et al., 2023).

These pathologies only disappear when restricting ECG to the strictly perturbative regime (EFT validity), with 2^{-2}2, i.e., the cubic term is always subleading compared to the Ricci scalar. Thus, there are no genuinely stable, static, spherically symmetric black-hole solutions in ECG with unsuppressed higher-order curvature terms (Felice et al., 2023).

6. Observational and Phenomenological Signatures

ECG predicts modest but potentially observable deviations from general relativity in strong-field regions (black holes, compact objects), but is essentially indistinguishable from GR in solar-system or weak-field regimes. Key predictions include:

  • Slight increases in the black-hole shadow angular diameter, with 2^{-2}3 (photon sphere) and 2^{-2}4 both shifted slightly outward at linear order in 2^{-2}5 (Hennigar et al., 2018, Li et al., 2024).
  • Enhanced periastron precession (compared to GR) for slightly positive cubic couplings, providing a direct but currently loose bound from S2 star orbits: 2^{-2}6 (Li et al., 2024).
  • Binary inspiral (“zoom–whirl”) waveforms in ECG show suppressed amplitudes and shifted phasing; future LISA or TianQin observations could constrain 2^{-2}7 at the percent level.
  • Gravitational lensing by ECG black holes results in milliarcsecond-level shifts in image positions for supermassive black holes, which are accessible to VLTI/GRAVITY and future EHT-like telescopes (Poshteh et al., 2018).
  • Shadow measurements at 2^{-2}8 precision would constrain the cubic coupling to 2^{-2}9 (Sánchez et al., 3 Feb 2025).

New classes of horizonless compact objects (“frozen gravitational stars,” FGSs) have been numerically constructed; these are naked singularities cloaked by a critical "frozen" surface at the Schwarzschild radius and are observationally indistinguishable from extremal black holes (Wang, 2024).

7. Cosmological Solutions and the Role of the Cubic Term

In cosmology, the cubic curvature term in ECG (and extended "cosmological Einsteinian cubic gravity," CECG) has a pronounced effect on early-universe dynamics. The phase-space analysis of the cosmological Friedmann equations—modified by the cubic term—reveals a universal “inflationary big-bang” attractor (for high curvature, P=12RabcdRcdefRefab+RabcdRcdefRefab12RabcdRacRbd+8RabRbcRca.\mathcal{P} = 12\,R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} + R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} -12\,R_{abcd}R^{ac}R^{bd} +8\,R_{a}{}^{b}R_{b}{}^{c}R_{c}{}^{a}.0), corresponding to non-standard matter-driven inflation that is the generic past behavior for all positive P=12RabcdRcdefRefab+RabcdRcdefRefab12RabcdRacRbd+8RabRbcRca.\mathcal{P} = 12\,R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} + R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} -12\,R_{abcd}R^{ac}R^{bd} +8\,R_{a}{}^{b}R_{b}{}^{c}R_{c}{}^{a}.1 (Quiros et al., 2020). However, the cubic modification alone cannot yield late-time acceleration without a cosmological constant term; the late-time attractor is standard de Sitter only if P=12RabcdRcdefRefab+RabcdRcdefRefab12RabcdRacRbd+8RabRbcRca.\mathcal{P} = 12\,R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} + R_{ab}{}^{cd}R_{cd}{}^{ef}R_{ef}{}^{ab} -12\,R_{abcd}R^{ac}R^{bd} +8\,R_{a}{}^{b}R_{b}{}^{c}R_{c}{}^{a}.2. The CECG model thus provides a purely geometric mechanism for primordial inflation but requires vacuum energy for late-time accelerated expansion.


References

  1. (Bueno et al., 2016, Hennigar et al., 2016, Bueno et al., 2016) — foundational construction, static black holes, thermodynamics, and uniqueness.
  2. (Felice et al., 2023) — instability and ghost analysis beyond EFT regime.
  3. (2002.04071, Bakhtiarizadeh, 2021, Cano et al., 2019, Adair et al., 2020) — rotating, charged black holes and black strings.
  4. (Bueno et al., 2018) — holographic dictionary and thermodynamics; (Hennigar et al., 2018, Sánchez et al., 3 Feb 2025, Poshteh et al., 2018, Li et al., 2024, Wang, 2024) — observational phenomenology, lensing, strong-field tests, and exotic compact objects.
  5. (Feng et al., 2017, Quiros et al., 2020) — cosmological solutions, critical points, and asymptotic analysis.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Einsteinian Cubic Gravity.