Cubic Metric-Affine Gravity (MAG)

• Cubic Metric-Affine Gravity (MAG) is a four-dimensional extension of GR where both the metric and affine connection are dynamic, incorporating torsion and nonmetricity.
• The inclusion of cubic invariants enriches the theory, allowing ghost-free propagation and stable vector, axial, and massive tensor modes through rigorous algebraic constraints.
• Novel solutions, such as modified black holes and gravitational waves with non-Einsteinian polarizations, highlight the potential observational signatures of cubic MAG.

Cubic Metric-Affine Gravity (MAG) is a four-dimensional gauge-theoretic extension of General Relativity wherein both the metric $g_{\mu\nu}$ and independent affine connection $\tilde\Gamma^\lambda{}_{\mu\nu}$ are dynamical, with their associated non-Riemannian features—torsion $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$ and nonmetricity $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$—entering the action through quadratic and, crucially, cubic invariants constructed from the curvature, torsion, and nonmetricity tensors. The inclusion of cubic terms results in a classical theory with a significantly richer geometric and dynamical structure, including new propagating degrees of freedom and exact solutions absent in the purely quadratic case, and most notably provides a consistent ghost-free framework for the propagation of vector, axial, and massive tensor modes that evade numerous no-go theorems for higher-spin fields (Bahamonde et al., 2024, Bahamonde et al., 5 Nov 2025).

1. Action Principle and Structural Features

The general cubic Metric-Affine Gravity action is given by

$$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$

where $-R$ is the usual Einstein–Hilbert term, and the Lagrangian densities $\mathcal{L}_{\rm Quad}$, $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$, and $\mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)}$ systematically incorporate all parity-even, metric-compatible quadratic and cubic invariants in the curvature, torsion, and nonmetricity.

Explicitly, the structure up to cubic order is

$\tilde\Gamma^\lambda{}_{\mu\nu}$0

$\tilde\Gamma^\lambda{}_{\mu\nu}$1

where $\tilde\Gamma^\lambda{}_{\mu\nu}$2, $\tilde\Gamma^\lambda{}_{\mu\nu}$3, $\tilde\Gamma^\lambda{}_{\mu\nu}$4, and $\tilde\Gamma^\lambda{}_{\mu\nu}$5 are free coupling constants detailed in the referenced appendices. These terms generate, respectively, second- and third-order invariants in the field strengths, with the cubic sector collecting all algebraically inequivalent invariants mixing curvature with up to two explicit torsion or nonmetricity tensors.

A distinct feature of cubic MAG is the large parameter space: the cubic sector alone can be parametrized by up to 209 couplings, subject to algebraic conditions required for physical consistency.

2. Field Equations and Irreducible Decomposition

Variation of the action with respect to the independent fields yields three sets of field equations:

• For the metric $\tilde\Gamma^\lambda{}_{\mu\nu}$6 (the “metric equation”):

$\tilde\Gamma^\lambda{}_{\mu\nu}$7

• For the torsion $\tilde\Gamma^\lambda{}_{\mu\nu}$8:

$\tilde\Gamma^\lambda{}_{\mu\nu}$9

• For the nonmetricity $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$0:

$$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$1

Key analytical steps for solution include separation of the connection into its Levi-Civita and distortion components, grouping terms according to variational derivatives, and exploiting the irreducible decomposition of torsion and nonmetricity. The latter can be summarized as:

• Torsion: vector $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$2, axial $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$3, tensor $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$4
• Nonmetricity: Weyl $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$5, vector $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$6, axio-tensor $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$7, trace-free $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$8

The explicit constitutive tensors and their index contractions, as functions of the quadratic and cubic couplings, govern the propagation and mixing of these modes.

3. Stability of Vector and Axial Modes

Quadratic MAG is well known for propagating unhealthy vector and axial degrees of freedom—manifesting as either Ostrogradsky ghosts or Hamiltonians unbounded from below—due to mixing terms such as $$T^\lambda{}_{\mu\nu} = 2\tilde\Gamma^\lambda{}_{[\mu\nu]}$$9 or $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$0, and cross-derivative couplings between inequivalent field irreps. The cubic invariants, when suitably tuned, can cancel these problematic interactions.

Precisely, the presence of 209 independent cubic couplings $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$1 allows imposition of 39 algebraic stability constraints, yielding a 194-dimensional “safe” parameter hypersurface. On this subspace, the kinetic matrix for the set $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$2 is block diagonal: $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$3 where the $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$4 are explicit linear combinations of the original couplings.

Ghost-free propagation and the absence of tachyons across all backgrounds require: $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$5

$$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$6

These algebraic conditions guarantee stability under linear perturbations for all irreducible vector and axial modes.

4. Black Hole and Wave Solutions in Cubic MAG

Cubic MAG admits novel exact solutions with no analog in either GR or quadratic MAG. Two notable classes are:

Reissner–Nordström-like Black Holes

The static, spherically symmetric ansatz allows non-vanishing irreducible torsion and nonmetricity components alongside the metric: $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$7

$$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$8

Upon solving the 123 algebraic stability constraints in addition to the field equations, the general solution features four independent charges:

• Mass $$Q_{\lambda\mu\nu} = \tilde\nabla_\lambda g_{\mu\nu}$$9
• Spin $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$0
• Dilation $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$1
• Shear $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$2

The metric function is modified vs. GR,

$$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$3

All vector/axial modes remain massless, while irreducible tensors $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$4 acquire universal, $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$5-independent mass terms proportional to a parameter $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$6.

Gravitational Wave (pp-wave) Solutions

New families of exact gravitational wave solutions exist in cubic MAG, exhibiting dynamical torsion and nonmetricity content (Bahamonde et al., 5 Nov 2025). The general pp-wave Ansatz is

$$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$7

with all torsion and nonmetricity components functions of $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$8 and Lie-dragged by the null vector $$S = \frac{1}{16\pi} \int d^4x \sqrt{-g} \left\{ -R + \mathcal{L}_{\rm Quad} + \mathcal{L}_{\rm curv\!-\!tors}^{(3)} + \mathcal{L}_{\rm curv\!-\!nonm}^{(3)} + \mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)} \right\}$$9.

Depending on which post-Riemannian sectors are present:

• In the Riemann–Cartan branch ($-R$0), a torsion function $-R$1 contributes an explicit term to $-R$2,
• In the Weyl–Cartan branch ($-R$3 but $-R$4), both $-R$5 and a harmonic scalar $-R$6 can be present,
• In the general metric-affine case, a traceless nonmetricity vector $-R$7 with at most linear $-R$8 dependence, alongside $-R$9 and $\mathcal{L}_{\rm Quad}$0, affects the metric function. The solution reads: $\mathcal{L}_{\rm Quad}$1 where the coefficients $\mathcal{L}_{\rm Quad}$2 depend on the fundamental couplings.

5. Polarization Structure and Physical Implications

The algebraic classification of pp-wave solutions in MAG, by Type N conditions on the field strengths, parallels the Petrov–Segre scheme in GR but must incorporate irreducible post-Riemannian components. The physical content is manifest in the geodesic deviation equation for test particles, where, in addition to the standard tensor polarizations ($\mathcal{L}_{\rm Quad}$3) of helicity $\mathcal{L}_{\rm Quad}$4, cubic MAG generically produces a breathing scalar mode ($\mathcal{L}_{\rm Quad}$5): $\mathcal{L}_{\rm Quad}$6 Whenever post-Riemannian amplitudes $\mathcal{L}_{\rm Quad}$7, the breathing mode is nonvanishing; thus, cubic MAG predicts a genuine helicity-0 polarization in gravitational radiation, a phenomenological signature in principle observable by gravitational wave interferometers sensitive to non-Einsteinian modes (e.g., LIGO/Virgo bounds on $\mathcal{L}_{\rm Quad}$8, $\mathcal{L}_{\rm Quad}$9). Present limits only exclude dominant scalar components, but advanced detectors may probe the subdominant regime ($\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$0–$\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$1) relevant for cubic-MAG waves.

Stability of these solutions is safeguarded by the same algebraic ghost-free constraints originating from the action-level analysis.

6. Massive Tensors, No-Go Theorems, and Theoretical Significance

Cubic MAG circumvents long-standing obstacles for consistently interacting higher-spin fields that afflict quadratic theories. In quadratic MAG, irreducible massless tensor fields associated with torsion and nonmetricity generically lead to strong-coupling, causality, and positivity pathologies, violating constraints such as the Weinberg–Witten theorem and the Velo–Zwanziger condition.

In the cubic framework, the same tensor modes ($\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$2, $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$3, $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$4) acquire genuine Proca-type mass terms proportional to the parameter $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$5:

• They propagate with healthy, bounded-from-below Hamiltonians,
• Avoid causality (Velo–Zwanziger) issues,
• Evade restrictions forbidding massless higher-spin conserved currents.

Consequently, in stable cubic MAG models, all propagating post-Riemannian degrees of freedom are either

• spin-0 (scalars: $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$6, $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$7, $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$8),
• spin-1 (vectors: $\mathcal{L}_{\rm curv\!-\!tors}^{(3)}$9, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$0, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$1, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$2),
• or massive spin-2–like fields ($\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$3, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$4, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$5);

no massless higher-spin degrees of freedom remain, sidestepping the classic no-go theorems and providing a self-consistent, stable low-energy classical theory for metric-affine gravity (Bahamonde et al., 2024, Bahamonde et al., 5 Nov 2025).

---

Sector	Degrees of Freedom	Status in Stable Cubic MAG
Vector/axial (Torsion, Nonmetricity)	$\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$6, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$7, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$8, $\mathcal{L}_{\rm curv\!-\!nonm}^{(3)}$9	Ghost-free, massless, block-diagonal kinetic
Irreducible tensors	$\mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)}$0, $\mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)}$1, $\mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)}$2	Massive, positive-definite Hamiltonians
Metric (spin-2, $\mathcal{L}_{\rm curv\!-\!tors\!-\!nonm}^{(3)}$3)	Graviton	Einsteinian in quadratic limit

A plausible implication is that observational signatures from non-Einsteinian gravitational waves or black holes carrying post-Riemannian fluxes would serve as a probe of the extended dynamical content of cubic metric-affine gravity.