Einsteinian Cubic Gravity Overview

• Einsteinian Cubic Gravity is a higher-derivative extension of general relativity featuring a unique cubic curvature invariant that propagates only the massless graviton.
• ECG yields tractable second-order equations for static, rotating, and charged black holes, leading to distinctive thermodynamic properties and stability features.
• Observable signatures in gravitational lensing, black hole shadows, and gravitational waves position ECG as a promising candidate for testing modifications to GR in strong-field regimes.

Einsteinian Cubic Gravity (ECG) is a higher-derivative extension of general relativity characterized by the addition of a parity-even, dimension-independent cubic curvature invariant to the gravitational action. The unique feature of ECG is that, when linearized about maximally symmetric backgrounds, its spectrum contains only the massless, transverse graviton, with all potential ghost or scalar modes decoupling. This property yields second-order equations for static, spherically symmetric solutions, ensuring tractable and physically meaningful modifications to black hole, cosmological, and compact object structure. ECG occupies a privileged position among cubic extensions of Einstein gravity, being neither trivial nor topological in four dimensions and unifying desirable features—nontrivial propagation, universal couplings, and absence of linearized ghosts—previously known only separately in Lovelock or topological gravities.

1. Action, Structure, and Linearized Spectrum

The ECG action in four dimensions is defined as

$$S = \frac{1}{16\pi G} \int d^4x\,\sqrt{-g}\bigl[ R - 2\Lambda + \lambda\,\mathcal{P} \bigr],$$

where $\lambda$ is a dimensionful (mass$^{-4}$) or dimensionless coupling depending on conventions, and

$$\mathcal{P} = 12\,R_{a}{}^{c}{}_{b}{}^{d}\, R_{c}{}^{e}{}_{d}{}^{f}\, R_{e}{}^{a}{}_{f}{}^{b} + 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}$$

is the unique, dimension-independent parity-even cubic curvature invariant selected to ensure that only the massless graviton propagates in the linearized theory (Bueno et al., 2016, Bueno et al., 2023).

Linearization about maximally symmetric backgrounds shows the absence of scalar and massive spin-2 ghosts, a property that distinguishes ECG from generic cubic gravity theories. The relative coefficients in $\mathcal{P}$ are chosen to be universal in dimension and recover the Gauss–Bonnet combination at quadratic order. Unlike Lovelock gravities, which are topological or trivial in four dimensions, ECG remains nontrivial and leads to genuine modifications of black hole and cosmological solutions (Bueno et al., 2016).

2. Black Hole Solutions: Static, Rotating, and Charged Cases

Static Spherically Symmetric Black Holes

ECG admits single-function generalizations of Schwarzschild and RN(A)dS metrics, the metric taking the form

$$ds^2 = -f(r)\,dt^2 + \frac{dr^2}{f(r)} + r^2\,d\Omega_2^2,$$

where $f(r)$ satisfies a non-linear ODE derived from varying the action. For vanishing cosmological constant,

$$-(f-1)\,r - \lambda\left[ \frac{1}{3}f^{\prime3} + \frac{1}{r}f^{\prime2} - \frac{2}{r^2}f(f-1)\,f' - \frac{1}{r}f\,f''\bigl(r f' - 2(f-1)\bigr)\right] = 2M$$

(Bueno et al., 2016, Hennigar et al., 2018, Li et al., 2024). The black hole's thermodynamics are given exactly in terms of horizon data (radius $r_h$, temperature $T$, and entropy $S$), with the first law $dM = T\,dS$ verified independently of the cubic coupling (Bueno et al., 2016).

Small ECG black holes display positive specific heat below a critical mass,

$M_{\rm crit} = \frac{16}{27}G^{-1/2}\lambda^{1/4},$

and are thus thermodynamically stable, in contrast to Schwarzschild (where all solutions have negative specific heat) (Bueno et al., 2016, Hennigar et al., 2018).

Rotating Black Holes

Slowly rotating black holes in ECG admit a metric of the form

$$ds^2 = -f(r)\,dt^2 + \frac{dr^2}{f(r)} + 2a\,r^2 p(r)\sin^2\theta\,dt\,d\phi + r^2(d\theta^2 + \sin^2\theta\,d\phi^2),$$

where $f(r)$ and $p(r)$ obey two coupled second-order ODEs, tractable numerically for arbitrary coupling parameter (Adair et al., 2020). ECG corrections modify the horizon radius, angular velocity, ISCO, photon sphere, and shadow in calculable fashion, with all quantities acquiring corrections $\sim\lambda/M^4$ at leading order (Adair et al., 2020, Burger et al., 2019).

Charged Black Holes and Nonuniqueness

Electrically and magnetically charged solutions exist, with the master equation

$$2GM-\frac{GQ^2}{r} =-(f-1)\,r - G^2\lambda\left[ 4f'^3+12\frac{f'^2}{r} -24f(f-1)\frac{f'}{r^2} -12\,f\,f''\bigl(f'-\tfrac{2(f-1)}{r}\bigr) \right].$$

Above a threshold coupling ($\lambda > M^4/768 G^2$), two regular, asymptotically flat black holes can share the same mass and charge—an explicit, pathology-free violation of the black hole uniqueness theorem in four dimensions (2002.04071). All charged ECG black holes lack a Cauchy horizon, evading mass-inflation instabilities and trivially satisfying strong cosmic censorship.

3. Stability, Spectrum, and EFT Constraints

Linearized perturbations (odd-parity Regge–Wheeler modes) show that, aside from the massless spin-2 graviton, ECG generates two additional branches with high-mass (or tachyonic) character: $$\omega^2 = c_r^2 k_r^2 + c_\theta^2 f k_\theta^2 + f\mathcal{M}^2,$$ with masses $\mathcal{M}^2 \sim r^3/(\lambda r_h)\gg 1/\ell^2$ for horizon radius $r_h \gg \ell$ (with EFT cutoff set by the scale $\ell$). These extra modes always lie far above the EFT window and do not contribute to low-energy dynamics (Bueno et al., 2023).

Within the EFT regime ($r_h \gg \ell,\, k \ll 1/\ell$), only the ordinary graviton propagates, and all metric corrections are under perturbative control. Any "pathology" observed (superluminality, instability) by probing $k\sim 1/\ell$ is an artifact of exceeding the EFT domain (Bueno et al., 2023).

4. Cosmology, Braneworlds, and Holography

Cosmological ECG (CECG) (Quiros et al., 2020) modifies the Friedmann–Raychaudhuri system: $$3H^2(1+16\beta H^4) = \rho_m + \Lambda, \qquad 2\dot H (1+48\beta H^4) = -(\rho_m + p_m).$$ The global phase-space admits three attractors: the "inflationary Big Bang" (source), standard matter-dominated (saddle), and a de Sitter manifold (late-time attractor with $\Lambda>0$). The vacuum spectrum remains the massless graviton; ghost and scalar modes decouple.

In five-dimensional warped braneworld scenarios (Lessa et al., 2023), ECG corrections are incorporated via the unique cubic term $P$, yielding second-order equations for the warp factor and ensuring localization and stability of the four-dimensional gravity on the brane, provided the cubic coupling stays below a certain bound.

In holography, ECG provides a well-defined AdS$_4$/CFT$_3$ toy model, where the black-hole phase diagram, entanglement entropy, and transport coefficients (notably, the shear viscosity to entropy ratio) can be computed exactly as non-analytic functions of the ECG coupling. No violation of the KSS bound occurs in the physical coupling domain (Bueno et al., 2018).

5. Observational Signatures and Tests

Geodesics, Shadows, and Lensing

ECG modifies marginally bound and innermost stable circular orbits (MBO/ISCO), photon spheres, and shadow radii, inducing observable corrections $\sim\lambda/M^4$ (Hennigar et al., 2018, Li et al., 2024, Adair et al., 2020). Periodic orbits, periastron precession, zoom-whirl structures, and gravitational waveforms show quantitative shifts as a function of the ECG coupling; best-fit to S2-star data around Sgr A* yields current limits $\lambda \lesssim 5\times10^{29} M_\odot^4$ (Li et al., 2024).

Gravitational lensing, including time delays and image positions, exhibits miliarcsecond-level differences from GR predictions for black holes at maximal allowed $\lambda$, potentially observable with current VLBI and interferometric techniques (Poshteh et al., 2018).

Regular Black Holes and Wormholes

Perturbative ECG corrections preserve regular, magnetically charged black holes sourced by nonlinear electrodynamics, with de Sitter-like cores, remnant horizon radius, and modified heat capacity phase structure (Lessa et al., 2023). Traversable wormhole solutions exist with analytic shape functions threaded by ordinary (weak-energy-condition-respecting) matter, as the cubic curvature stress makes up the required exotic energy near the throat (Mehdizadeh et al., 2019).

Accretion, Blandford–Znajek Process, and Spin Tests

ECG modifies the Blandford–Znajek jet power at second relative order in black hole spin,

$$P_{\rm BZ} = \Phi^2 \left[ \frac{1}{6\pi} (M\Omega_H) + \frac{\alpha + \beta\zeta}{6\pi} (M\Omega_H)^3 + \cdots \right],$$

where the cubic coupling appears in the coefficient $\beta$ at quartic order in spin. Rapidly spinning holes (e.g., M87) can thus break the GR-ECG degeneracy observationally (Peng et al., 2023).

Static, horizonless ECG solutions admit unique, static stable circular orbits, yielding non-monotonic ZAMO velocities (the Aschenbach effect) and Doppler-invariant accretion rings—signatures distinguishable from GR (Zhao et al., 26 Jan 2026).

6. Weyl-Gauged ECG and Further Extensions

Weyl–gauged ECG (WECG) (Dengiz, 2024) is constructed in Weyl geometry, with a local compensating scalar and Weyl vector, and is free from any dimensionful couplings. The bare ECG emerges as the low-energy limit after spontaneous or radiative breaking of Weyl symmetry. The model admits AdS or flat vacua (but not dS), with the Planck scale set by the scalar expectation value. The possibility of embedding the cubic invariants in Weyl geometry and the role of nonperturbative effects (instantons, tachyon condensation) in vacuum stability remain open research directions.

7. Scattering Amplitudes, Effective Potentials, and Theoretical Bounds

On-shell amplitude techniques efficiently compute the semi-classical two-body potentials in ECG, with new non-dispersive terms at $\sim 1/r^6$ and quantum $\sim 1/r^7$ corrections. The leading singularity method bypasses full loop integrations and Tensors, providing direct access to classical GR and ECG modifications in the scattering angle, impulse, and black hole metric (Emond et al., 2019, Burger et al., 2019).

No current observation constrains the cubic coupling to meaningful theoretical bounds for astrophysical objects; solar-system limits are many orders of magnitude weaker than those accessible in strong-field regimes near compact objects.

---

In summary, Einsteinian Cubic Gravity is the unique, universal, ghost-free cubic extension of GR that admits well-defined static, rotating, and charged black holes, regular wormholes, cosmologies, and braneworlds, all with distinctive corrections calculable in terms of a single cubic curvature coupling. Its observable signatures—thermodynamics, geodesics, lensing, accretion physics, and gravitational waves—remain under active investigation in theory and experiment. The consistent EFT regime, holographic implications, and extensions (such as Weyl–gauged ECG and higher-dimensional formulations) highlight ECG as the leading candidate for nontrivial cubic modifications to four-dimensional gravity (Bueno et al., 2016, Bueno et al., 2023, Bueno et al., 2016, Adair et al., 2020, Bueno et al., 2018, Li et al., 2024, Poshteh et al., 2018, Zhao et al., 26 Jan 2026).