Quadratic Gravity: Action and Implications

• Quadratic gravity is a framework that adds R², Ricci, and Riemann squared terms to the conventional Einstein–Hilbert action, thereby capturing higher-order UV corrections.
• It produces fourth-order field equations that introduce new propagating degrees of freedom, including a massive spin-2 ghost and a massive scalar mode.
• The theory’s improved renormalizability motivates investigations into nonlocal modifications and alternative formulations to reconcile quantum effects with gravity.

Quadratic Action of Gravity

A quadratic action of gravity is a generally covariant field theory in which the Lagrangian includes all possible terms that are quadratic in the curvature tensors—specifically, in the scalar curvature $R$, Ricci tensor $R_{\mu\nu}$, and Riemann tensor $R_{\mu\nu\lambda\sigma}$—in addition to the conventional Einstein–Hilbert term. Such theories are motivated by ultraviolet (UV) renormalizability, effective field theory corrections to General Relativity (GR), and attempts to achieve a deeper ultraviolet completion of quantum gravity. The quadratic action modifies both the classical field equations and the quantum properties of gravity, introducing characteristic new propagating degrees of freedom, higher-derivative dynamics, and intricate issues regarding unitarity and stability.

1. Structure of the Quadratic Gravitational Action

The most general parity-even, diffeomorphism-invariant Lagrangian in four dimensions involving curvatures up to mass dimension four is

$$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$

where $\kappa^2 = 32\pi G$, with $G$ Newton's constant, and $\alpha, \beta, \gamma$ are dimensionless couplings. The cosmological constant $\Lambda$ is often set to zero for theoretical analysis. Utilizing the Gauss–Bonnet theorem, which in $d=4$ makes the combination

$$E = R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} - 4 R_{\mu\nu}R^{\mu\nu} + R^2$$

a total derivative, one typically reduces the three quadratic invariants to two independent ones, often choosing $R_{\mu\nu}$0 and the Weyl tensor squared $R_{\mu\nu}$1 as the basis. The action can then be recast as

$R_{\mu\nu}$2

with $R_{\mu\nu}$3 real, dimensionless parameters controlling the strengths of the quadratic terms (Donoghue et al., 2021, Salvio, 2018, Alvarez-Gaume et al., 2015, Silveravalle, 2022).

2. Classical Field Equations and Their Properties

Varying the action with respect to $R_{\mu\nu}$4 produces field equations of up to fourth order in derivatives: \begin{align*} \frac{1}{\kappa^2}(R_{\mu\nu} - \tfrac12 g_{\mu\nu} R) &+ 2\alpha \left[ R R_{\mu\nu} - \tfrac12 g_{\mu\nu} R^2 + g_{\mu\nu}\Box R - \nabla_\mu \nabla_\nu R \right] \ &+ \beta \left[ \Box R_{\mu\nu} + \tfrac12 g_{\mu\nu} \Box R - \nabla_\mu \nabla_\nu R - 2 R_{\mu\rho\nu\sigma}R^{\rho\sigma} + \tfrac12 g_{\mu\nu} R_{\rho\sigma} R^{\rho\sigma} \right] \ &+ \gamma \left[ -\tfrac12 g_{\mu\nu} R_{\alpha\beta\rho\sigma}R^{\alpha\beta\rho\sigma} + 2 R_{\mu\alpha\beta\gamma}R_\nu{}^{\alpha\beta\gamma} - 4\Box R_{\mu\nu} + 2g_{\mu\nu} \Box R + 4 \nabla_\mu \nabla_\nu R \right] \ &= 0 \end{align*} Given the presence of up to four derivatives of the metric, the theory propagates more degrees of freedom than GR and exhibits richer solution spaces—including, for example, black holes with non-Schwarzschild exteriors, wormholes, and solutions with different horizon structures or Yukawa-type potentials (Silveravalle, 2022).

When written with nonlocal or analytic functions of the d'Alembertian acting on curvature, as in

$R_{\mu\nu}$5

the variation rules require careful treatment of operators acting on fields and their variations, with higher-derivative structures handled by repeated integration by parts and the use of algebraic identities (Biswas et al., 2013).

3. Spectrum and Propagators: New Degrees of Freedom and Ghosts

The linearized spectrum around flat spacetime reveals the physical content. Perturbing $R_{\mu\nu}$6, the kinetic terms can be diagonalized using spin projection operators.

The physical spectrum consists of:

• A massless spin-2 graviton (as in GR)
• A massive spin-2 "ghost" excitation with mass $R_{\mu\nu}$7 and negative residue (Ostrogradsky–Lee–Wick ghost)
• A massive spin-0 scalar mode with mass $R_{\mu\nu}$8 if $R_{\mu\nu}$9 (healthy if sign chosen appropriately)

The propagator in momentum space reveals the characteristic $R_{\mu\nu\lambda\sigma}$0 fall-off at high energies: $R_{\mu\nu\lambda\sigma}$1 where $R_{\mu\nu\lambda\sigma}$2 and $R_{\mu\nu\lambda\sigma}$3 are the Barnes–Rivers spin projector operators (Donoghue et al., 2021, Salvio, 2018, Kubo et al., 5 Feb 2025).

For ghost-free nonlocal theories, entire functions $R_{\mu\nu\lambda\sigma}$4 in the propagator denominator can suppress new poles, avoiding physical ghosts if $R_{\mu\nu\lambda\sigma}$5 has no zeros in the complex plane (example: $R_{\mu\nu\lambda\sigma}$6) (Biswas et al., 2013).

A unique feature arises in the pure $R_{\mu\nu\lambda\sigma}$7 theory: on flat backgrounds, only the scalar propagates, while the graviton mode can be absent; but on de Sitter/anti–de Sitter backgrounds or with an added Einstein term, the massless spin-2 returns (Alvarez-Gaume et al., 2015).

4. Renormalization, UV Behavior, and Asymptotics

Quadratic gravity is perturbatively renormalizable: all divergences up to dimension-4 operators can be absorbed into the couplings of $R_{\mu\nu\lambda\sigma}$8, $R_{\mu\nu\lambda\sigma}$9, the cosmological constant, and the Newton constant. The superficial degree of divergence is reduced due to the $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$0 propagator scaling, greatly improving the UV behavior over GR (Donoghue et al., 2021, Salvio, 2018).

One-loop $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$1-functions for the dimensionless couplings, in the pure gravity and generic matter sector, show: $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$2 where $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$3 accounts for the number and type of matter fields. The theory is thus asymptotically free in the $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$4 coupling, while the $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$5 coupling typically grows towards the UV.

If all matter couplings flow to a fixed point, the theory can approach conformal gravity at infinite energy, evading the Landau pole (“agravity” scenario) (Salvio, 2018).

Nonlocal, ghost-free variants can achieve UV softening and even asymptotic freedom while preserving unitarity, at the cost of introducing nonlocal operators (Biswas et al., 2013).

5. Unitarity, the Ghost Problem, and Physical Interpretations

The ghost-like massive spin-2 excitation (Ostrogradsky mode) is a central challenge. At the perturbative level:

• The ghost has a negative-norm or negative-energy contribution in the Hilbert space (signaled by the wrong sign of the kinetic pole).
• Stability is spoiled, as the Hamiltonian is unbounded below.

Proposed remedies include:

• Lee–Wick prescription: Interpreting the ghost as an unstable resonance, choosing analytic continuations (prescriptions in the propagator) such that only unitary asymptotic states contribute, enforcing unitarity for physical scattering amplitudes.
• Indefinite-metric or Dirac–Pauli quantization: Quantizing the ghost sector with modified inner product (possibly including complexified phase-space variables), restoring probabilistic interpretation in the sense that all transition probabilities are nonnegative and sum to one (Salvio, 2024).
• In nonlocal quadratic models, choosing form factors so that no new poles appear in the propagator—thus avoiding physical ghosts (Biswas et al., 2013). However, this shifts the issue to the nature of nonlocal modification and possible UV acausality.

Open problems include nonperturbative construction of the physical Hilbert space, identification of a true ghost-free sector, and mathematical control of quantum probabilities.

6. Physical Consequences and Applications

Table: New Physical Features from Quadratic Action Terms

Property	Description	References
Renormalizability	Power counting UV behavior improved; all divergences absorbed	(Donoghue et al., 2021, Salvio, 2018)
New modes	Propagation of a massive spin-2 ghost and a massive scalar	(Silveravalle, 2022, Alvarez-Gaume et al., 2015)
Modified Newtonian potential	Yukawa corrections, modified black hole solutions	(Silveravalle, 2022)
Black holes/wormholes/naked singularities	Exotic spherically symmetric solutions beyond Schwarzschild	(Silveravalle, 2022)
Inflationary cosmology	Starobinsky-like inflation from $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$6 term	(Kubo et al., 5 Feb 2025, Alvarez-Gaume et al., 2015)
UV completeness (potential)	Well-defined up to infinite energy for special RG flows	(Salvio, 2018)
Quantum path integrals	Formulation of background-independent/Euclidean measures	(Salvio, 2024, Belokurov et al., 2024)
Unification scenarios	Emerges naturally in string compactifications or Weyl–DBI frameworks	(Ghilencea, 14 Aug 2025, Alvarez-Gaume et al., 2015)

• Astrophysical implications: The quadratic terms lead to modifications in the post-Newtonian dynamics of binary inspirals, impacting gravitational waveforms. Current experiments constrain the couplings of quadratic terms (Kim et al., 2019).
• Cosmology: The $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$7 term supports inflationary dynamics (Starobinsky inflation), with predictions for primordial gravitational waves and consistency relations in the inflationary tensor spectrum (Kubo et al., 5 Feb 2025).
• Static objects: The Schwarzschild metric is no longer unique; black holes with additional hair, wormhole solutions, or naked singularities can be realized, controlled by the extra Yukawa charges associated with the scalar and ghost masses (Silveravalle, 2022).

7. Alternative Formulations and Extensions

• First order (Palatini, metric-affine) formalism: Quantizing both the metric and independent connection fields leads to an inequivalent renormalizable theory. In this approach, all propagators fall as $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$8, and no explicit spin-2 ghost arises, at least at the linearized level (Alvarez et al., 2017, Alvarez et al., 2017). The theory exhibits UV conformal invariance and allows for the Planck scale to be dynamically generated by spontaneous symmetry breaking via matter couplings.
• Nonlocal/infinite-derivative gravity: Nonlocal analytic form factors $$S = \int d^4x\,\sqrt{-g} \left[ \frac{1}{2\kappa^2} R + \alpha\,R^2 + \beta\,R_{\mu\nu}R^{\mu\nu} + \gamma\,R_{\mu\nu\lambda\sigma}R^{\mu\nu\lambda\sigma} + \Lambda \right]$$9 can regulate the spectrum, eliminating ghosts and rendering the theory both UV-soft and unitary for suitable analytic choices (Biswas et al., 2013).
• Weyl-Dirac-Born-Infeld action: In unified gauge-theoretic frameworks, a Born–Infeld-type action in Weyl geometry automatically yields quadratic gravity at leading order in a dimensionless small parameter, with gravity self-regularizing UV divergences and providing a natural matching to the Einstein–Hilbert term via dynamical symmetry breaking (Ghilencea, 14 Aug 2025).
• Torsionful extensions: Incorporating quadratic invariants in the torsion sector introduces further scalar and axial-vector modes, modifying Friedmann dynamics, effectively providing alternative "inflationary" or acceleration branches without explicit scalar fields (Iosifidis et al., 2021).

---

Quadratic gravity, in its various forms—local, nonlocal, metric-affine, or embedded in larger symmetry frameworks—remains a theoretically compelling, technically nontrivial, and actively investigated framework for both classical and quantum gravitational physics. The delicate balance between renormalizability, unitarity, novel phenomenology, and compatibility with existing observations continues to drive research at the intersection of gravity, field theory, and cosmology (Donoghue et al., 2021, Salvio, 2018, Biswas et al., 2013, Kubo et al., 5 Feb 2025).