Warped Massive Gravity: Theory & Applications

• Warped massive gravity is a framework that unifies massive graviton models with warped extra dimensions across six-, five-, and three-dimensional settings.
• In six dimensions, warping generates an infinite tower of Kaluza–Klein modes that yield Yukawa corrections to the four-dimensional Newtonian potential, exemplifying computable model dynamics.
• Lower-dimensional realizations enhance our understanding of holography and asymptotic symmetries, providing explicit insights into black hole physics and cosmological behavior.

Warped massive gravity denotes a family of constructions in which a massive graviton sector is combined with a warped geometry. In the literature represented here, the term covers at least three technically distinct settings: a warped six-dimensional world with an extra two-sphere, where four-dimensional gravity receives a tower of Kaluza–Klein Yukawa corrections (Kokado et al., 2018); a four-dimensional dRGT theory embedded on a brane in a warped five-dimensional AdS bulk, where the bulk induces an additional kinetic term for the longitudinal mode and suppresses both the vDVZ discontinuity and low-scale strong coupling (Gabadadze et al., 2019); and three-dimensional massive-gravity theories whose distinguished solutions are warped AdS$_3$ or warped BTZ geometries with $SL(2,\mathbb R)\times U(1)$ structure and warped-CFT asymptotics (0906.1243, 0902.4634). A recent cosmological formulation uses a five-dimensional ghost-free massive graviton with a brane-localized four-dimensional massive gravity potential and studies FLRW evolution on the brane (Garcia-Saenz et al., 11 Sep 2025).

1. Core geometries and uses of the term

A first higher-dimensional realization is a warped $6$D geometry with line element

$$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$

so that a four-dimensional Minkowski sector is warped by a function of the internal angular coordinate while the extra space is a round $S^2$ of fixed radius $a$. In that model the positive-energy condition fixes

$\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$

with constraints $4a\le 1$ and $\Lambda a^2<-2$ (Kokado et al., 2018).

A second usage places nonlinear four-dimensional massive gravity on a brane inside a warped AdS$_5$ bulk with

$SL(2,\mathbb R)\times U(1)$0

and splits the action as $SL(2,\mathbb R)\times U(1)$1. In this construction the bulk curvature is not merely background structure: it generates a kinetic term for the four-dimensional helicity-0 mode and changes the decoupling-limit dynamics of the massive graviton (Gabadadze et al., 2019).

A third usage, dominant in the three-dimensional literature, concerns massive gravity theories admitting warped AdS$SL(2,\mathbb R)\times U(1)$2 backgrounds and warped black holes. In those works the relevant symmetry is reduced from $SL(2,\mathbb R)\times U(1)$3 to $SL(2,\mathbb R)\times U(1)$4, and the asymptotic symmetry algebra is typically a Virasoro algebra semidirectly extended by a $SL(2,\mathbb R)\times U(1)$5 current algebra rather than the Brown–Henneaux $SL(2,\mathbb R)\times U(1)$6 structure (0906.1243, Yekta, 2016). This suggests that “warped massive gravity” is best treated as a class of related frameworks rather than a single universal Lagrangian.

2. Six-dimensional warped worlds and the KK correction to Newtonian gravity

In the six-dimensional model, the bulk stress tensor is taken in the general diagonal form

$SL(2,\mathbb R)\times U(1)$7

with

$SL(2,\mathbb R)\times U(1)$8

This ansatz does not commit to a particular microscopic bulk matter sector, and conservation follows from the Bianchi identity. With the warp factor above, explicit expressions are obtained for $SL(2,\mathbb R)\times U(1)$9, such as

$6$0

The Einstein equation is written as

$6$1

or equivalently

$6$2

where the factor $6$3 is specific to $6$4D (Kokado et al., 2018).

After linearization and gauge fixing, the graviton fluctuation is separated as

$6$5

with periodic $6$6-dependence $6$7. For small warp parameter $6$8, the $6$9-equation reduces to a spheroidal differential equation,

$$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$0

with $$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$1. The KK decomposition then yields four-dimensional modes satisfying

$$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$2

and, for small warp parameter,

$$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$3

The normalization condition is

$$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$4

These facts make the KK spectrum explicit rather than purely formal (Kokado et al., 2018).

The four-dimensional potential between masses $$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$5 and $$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$6 is

$$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$7

with $$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$8 and $$ds^2 = g_{\mu\nu}(x,\theta)\,dx^\mu dx^\nu + a^2\left(d\theta^2+\sin^2\theta\, d\phi^2\right), \qquad g_{\mu\nu}=\phi^2(\theta)\,\eta_{\mu\nu},$$9. The potential is therefore a Newtonian term plus an infinite tower of Yukawa corrections. At short distance, the large-$S^2$0 tail satisfies $S^2$1 as $S^2$2, so the correction behaves as $S^2$3 near the origin (Kokado et al., 2018). In this sense, the model is a computable warped realization of massive-gravity corrections to four-dimensional gravity.

3. Warped AdS$S^2$4 embeddings, longitudinal dynamics, and amplitude structure

The five-dimensional AdS embedding starts from the flat-space problem of ordinary four-dimensional massive gravity, where the helicity-0 mode $S^2$5 obtains its dynamics only through mixing with the helicity-2 field. In the schematic decoupling-limit brane action,

$S^2$6

diagonalization produces a matter coupling

$S^2$7

and the flat-space theory becomes strongly coupled at

$S^2$8

That low scale underlies the usual strong-coupling problem, while the non-decoupling of the scalar polarization in the limit $S^2$9 yields the vDVZ discontinuity (Gabadadze et al., 2019).

The warped solution is to couple the brane theory to a weakly coupled AdS$a$0 bulk. In the bulk, the longitudinal Stückelberg vector has action

$a$1

so curvature generates a mass for the Goldstone vector. In the $a$2D decoupling limit, the bulk longitudinal scalar $a$3 obeys

$a$4

which induces on the brane a nonlocal kinetic term for $a$5,

$a$6

At low energy this becomes approximately local,

$a$7

The scalar is then canonically normalized as

$a$8

and its matter coupling is suppressed by $a$9. The theory can remain weakly coupled up to energies of order $\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$0, and the linear theory smoothly approaches GR, so the vDVZ discontinuity disappears (Gabadadze et al., 2019).

The same construction has a holographic interpretation in which a four-dimensional nonlinear massive gravity mixes with a non-conserved symmetric tensor of a cutoff CFT, and that mixing generates the large kinetic term for the longitudinal mode. In this dual picture, the CFT continuum supplies the missing dynamics of $\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$1, weakens its self-interactions, and weakens its coupling to the matter stress tensor (Gabadadze et al., 2019).

An amplitude-based realization of warped KK massive gravity appears in compactified five-dimensional warped gauge and gravity theories on an RS1 background,

$\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$2

After KK reduction, the five-dimensional graviton yields a four-dimensional tensor $\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$3, vector $\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$4, and scalar $\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$5. In $\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$6 gauge, the gravitational equivalence theorem connects helicity-$\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$7 KK gravitons to scalar Goldstones and helicity-$\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$8 KK gravitons to vector Goldstones at high energy, and the paper further constructs an extended double-copy for massive KK amplitudes: exact for full three-point amplitudes at tree level and valid at leading order in the high-energy expansion for four-point amplitudes (Hang et al., 2024). This provides an on-shell formulation of warped massive gravity in which the KK graviton sector is controlled by the same equivalence-theorem logic that organizes spontaneously broken gauge theories.

4. Cosmology in warped massive gravity

In the cosmological formulation, warped massive gravity is defined by a five-dimensional ghost-free massive graviton in the bulk together with a four-dimensional dRGT potential localized on the brane. The action contains a $\phi(\theta)=\epsilon\, e^{a\sin^2\theta},$9D Einstein–Hilbert term with $4a\le 1$0, a $4a\le 1$1D dRGT potential $4a\le 1$2, a $4a\le 1$3D induced Einstein term, and a $4a\le 1$4D dRGT potential $4a\le 1$5. The fiducial geometry is AdS in flat slicing,

$4a\le 1$6

and the construction is explicitly motivated by the same idea that the warped AdS structure upgrades the effective strong-coupling scale of the four-dimensional massive-gravity sector (Garcia-Saenz et al., 11 Sep 2025).

Because flat and closed FLRW geometries are generally inconsistent with a flat fiducial metric in massive gravity, the analysis uses the open FLRW branch. The bulk ansatz is

$4a\le 1$7

with induced brane metric

$4a\le 1$8

The Stückelberg equation gives

$4a\le 1$9

and nontrivial cosmology requires

$\Lambda a^2<-2$0

A useful consequence is

$\Lambda a^2<-2$1

so the four-dimensional dRGT contribution acts like an effective vacuum-energy term on the brane (Garcia-Saenz et al., 11 Sep 2025).

The cosmological system simplifies into two distinguished classes. In the Neumann model, one imposes

$\Lambda a^2<-2$2

The brane dynamics is then governed by a modified Raychaudhuri equation containing the term

$\Lambda a^2<-2$3

Numerically, this model exhibits Big Bang type solutions, Big Crunch type solutions, and a broad class with a non-singular bounce. Near the bounce, if $\Lambda a^2<-2$4 is small,

$\Lambda a^2<-2$5

so the bounce is smooth and does not require exotic matter. The same model can also admit a Big Rip behavior in some parameter regions (Garcia-Saenz et al., 11 Sep 2025).

In the special model, one instead tunes the massive-gravity coefficients so that

$\Lambda a^2<-2$6

This tuning implies $\Lambda a^2<-2$7, removes the problematic $\Lambda a^2<-2$8 term, and permits a first integral of the Raychaudhuri equation. For $\Lambda a^2<-2$9, the resulting Friedmann equation has a DGP-like two-branch structure. At late times it behaves like GR plus a cosmological constant, although the early-time evolution can differ substantially; special parameter choices can even give unconventional late-time scalings such as $_5$0 (Garcia-Saenz et al., 11 Sep 2025).

The observational analysis in that work specializes to $_5$1 and uses Planck 2018 CMB likelihoods together with PantheonPlusSH0ES supernovae. The special model is reported to be generally compatible with current data at a level comparable to $_5$2CDM, and for priors allowing negative $_5$3 it can reduce or even remove the Hubble tension, with the strongest case giving consistency of the CMB-inferred and SH0ES/PantheonPlus-inferred $_5$4 values at about $_5$5. The authors also stress several limitations: only background cosmology is considered, perturbative stability is not analyzed, the induced fiducial metric is assumed flat, the simplest boundary condition $_5$6 is imposed, $_5$7 symmetry across the brane is assumed, and the observational study covers only a restricted subset of the special-model parameter space (Garcia-Saenz et al., 11 Sep 2025).

5. Warped AdS$_5$8 sectors in three-dimensional massive gravity

Three-dimensional massive gravity provides a separate and highly developed notion of warped massive gravity. In cosmological new massive gravity, the action

$_5$9

admits both BTZ black holes and warped AdS$SL(2,\mathbb R)\times U(1)$00 black holes. Stationary rotationally symmetric solutions are constructed by reducing the metric to a mechanical system for a Minkowski three-vector $SL(2,\mathbb R)\times U(1)$01 and using the quadratic ansatz $SL(2,\mathbb R)\times U(1)$02. The warped branch is characterized by a parameter $SL(2,\mathbb R)\times U(1)$03 fixed by the couplings, and the causally regular black-hole range is

$SL(2,\mathbb R)\times U(1)$04

At the critical point $SL(2,\mathbb R)\times U(1)$05, the regular solution is instead $SL(2,\mathbb R)\times U(1)$06 with vanishing entropy, mass, and angular momentum (0902.4634).

The same geometric sector appears in BHT massive gravity and in generalized minimal massive gravity. In GMMG, the spacelike stretched warped AdS$SL(2,\mathbb R)\times U(1)$07 black hole is written in ADM form,

$SL(2,\mathbb R)\times U(1)$08

with $SL(2,\mathbb R)\times U(1)$09, $SL(2,\mathbb R)\times U(1)$10, and $SL(2,\mathbb R)\times U(1)$11 fixed by the warp parameter $SL(2,\mathbb R)\times U(1)$12 and the horizon radii $SL(2,\mathbb R)\times U(1)$13. The geometry has isometry $SL(2,\mathbb R)\times U(1)$14, and the auxiliary one-forms are chosen as

$SL(2,\mathbb R)\times U(1)$15

so that the field equations reduce to algebraic conditions on the constants $SL(2,\mathbb R)\times U(1)$16 and the couplings $SL(2,\mathbb R)\times U(1)$17 (Setare et al., 2017).

The warped AdS$SL(2,\mathbb R)\times U(1)$18 phase space is larger than the stationary black-hole family. In ghost-free three-dimensional massive gravity including TMG, NMG, and MMG contributions,

$SL(2,\mathbb R)\times U(1)$19

there exist exact Kerr–Schild deformations of warped AdS$SL(2,\mathbb R)\times U(1)$20 black holes describing evanescent gravitons. Their profile satisfies

$SL(2,\mathbb R)\times U(1)$21

and yet the conserved gravitational energy is independent of the wave profile and equals that of the undeformed black hole. The higher-curvature terms modify the allowed decay exponent $SL(2,\mathbb R)\times U(1)$22, so these configurations are not locally equivalent to empty warped AdS$SL(2,\mathbb R)\times U(1)$23 once the full theory is taken into account (Giribet et al., 2015). A plausible implication is that warped massive gravity in three dimensions should be viewed as a genuine dynamical sector, not merely as a family of quotients of homogeneous backgrounds.

6. Boundary conditions, asymptotic charges, and warped holography

The foundational asymptotic analysis in topologically massive gravity introduces boundary conditions adapted to spacelike and timelike warped AdS$SL(2,\mathbb R)\times U(1)$24, rather than Brown–Henneaux asymptotics. For the spacelike case, the asymptotic symmetry generators $SL(2,\mathbb R)\times U(1)$25 and $SL(2,\mathbb R)\times U(1)$26 obey

$SL(2,\mathbb R)\times U(1)$27

and the corresponding charge algebra contains a Virasoro algebra together with a current algebra,

$SL(2,\mathbb R)\times U(1)$28

with

$SL(2,\mathbb R)\times U(1)$29

For the timelike case, the sign of the current-algebra central term flips. The charges are finite, integrable, and conserved, but the energy contribution of current descendants is sector-dependent: in the spacelike case they lower the energy, whereas in the timelike case they raise it (0906.1243).

In BHT massive gravity, the canonical Hamiltonian analysis in first-order form leads to an improved differentiable generator $SL(2,\mathbb R)\times U(1)$30, with surface term

$SL(2,\mathbb R)\times U(1)$31

The asymptotic symmetry algebra is a semidirect sum of a Virasoro algebra and a $SL(2,\mathbb R)\times U(1)$32 Kac–Moody algebra,

$SL(2,\mathbb R)\times U(1)$33

with

$SL(2,\mathbb R)\times U(1)$34

A Sugawara construction yields two commuting Virasoro algebras, and the resulting charges satisfy the first law, the Smarr relation, and the Cardy entropy formula for warped black holes (Yekta, 2016).

GMMG exhibits the same basic holographic pattern. The asymptotic conserved charges realize a Virasoro $SL(2,\mathbb R)\times U(1)$35 algebra with central extensions $SL(2,\mathbb R)\times U(1)$36 and $SL(2,\mathbb R)\times U(1)$37, a Sugawara-like construction reconstructs two commuting Virasoro generators, and both the ordinary Cardy formula and the warped Cardy formula reproduce the gravitational entropy exactly. On that basis the dual theory is argued to be a warped CFT rather than a conventional two-dimensional CFT with full left-right symmetry (Setare et al., 2017).

A more recent NMG analysis adapts “quadratic ensemble” boundary conditions to warped AdS$SL(2,\mathbb R)\times U(1)$38. In Fefferman–Graham gauge it imposes

$SL(2,\mathbb R)\times U(1)$39

and finds an asymptotic solution space with two constants $SL(2,\mathbb R)\times U(1)$40, $SL(2,\mathbb R)\times U(1)$41 and two chiral functions $SL(2,\mathbb R)\times U(1)$42, $SL(2,\mathbb R)\times U(1)$43. The asymptotic algebra is $SL(2,\mathbb R)\times U(1)$44, or after central extension $SL(2,\mathbb R)\times U(1)$45. The surface charges are finite but not integrable on the full phase space; integrability is restored by restricting to

$SL(2,\mathbb R)\times U(1)$46

Within that subsector, warped BTZ black holes are embedded naturally, and the warped Cardy formula reproduces the thermodynamic entropy (Sajadi et al., 8 Oct 2025).

The holographic interpretation is reinforced by hidden conformal symmetry. For the self-dual warped AdS$SL(2,\mathbb R)\times U(1)$47 black hole in TMG, the massive scalar wave equation becomes hypergeometric after the small angular momentum condition

$SL(2,\mathbb R)\times U(1)$48

and can be rewritten as the eigenvalue equation of an $SL(2,\mathbb R)\times U(1)$49 Casimir,

$SL(2,\mathbb R)\times U(1)$50

Because the geometry is obtained by a periodic identification in the angular direction, one hidden conformal sector is broken to $SL(2,\mathbb R)\times U(1)$51 while the other remains unbroken. The Cardy formula with

$SL(2,\mathbb R)\times U(1)$52

reproduces the Bekenstein–Hawking entropy, and the scalar absorption cross section and quasinormal frequencies take the expected CFT form (Li et al., 2010).

7. Perturbations, stability, and thermodynamic-chaotic probes

Exact perturbative probes provide some of the sharpest tests of warped massive gravity. For the self-dual warped AdS$SL(2,\mathbb R)\times U(1)$53 black hole in TMG, scalar, vector, and spinor perturbations all reduce to hypergeometric equations after separation of variables. Imposing ingoing behavior at the horizon and vanishing Dirichlet boundary conditions at infinity yields quantized frequencies of the universal form

$SL(2,\mathbb R)\times U(1)$54

with appropriate conformal weights for each spin. The resulting quasinormal spectrum matches the pole structure of retarded Green’s functions in the dual chiral CFT and thus supplies a quantitative test of warped AdS/CFT (Li et al., 2010).

For spacelike stretched warped AdS$SL(2,\mathbb R)\times U(1)$55 black holes in TMG, the massive scalar equation

$SL(2,\mathbb R)\times U(1)$56

is also exactly solvable in terms of hypergeometric functions. The effective potential approaches a constant $SL(2,\mathbb R)\times U(1)$57 at infinity rather than the asymptotically flat value $SL(2,\mathbb R)\times U(1)$58, and superradiance occurs when the usual frequency condition is satisfied. Nevertheless, the quasinormal and bound-state spectra show that all physically admissible modes satisfy

$SL(2,\mathbb R)\times U(1)$59

Even after placing a stationary mirror outside the horizon and imposing Dirichlet boundary conditions, no superradiant instability appears because the geometry does not provide the trapping well needed for repeated amplification (Ferreira, 2013). This result is frequently contrasted with Kerr, whose rotating superradiant modes can become unstable in the presence of confinement.

Warped massive gravity has also been analyzed through chaos and thermodynamic phase structure. In warped BTZ/WCFT holography, the appropriate diagnosis uses effective sector-dependent temperatures rather than a single global temperature. For rotating BTZ,

$SL(2,\mathbb R)\times U(1)$60

and the same logic extends to warped BTZ black holes and WCFTs, where the Virasoro and $SL(2,\mathbb R)\times U(1)$61 sectors feel different effective temperatures. Within that framework, the apparent violation of the chaos bound is resolved, the grand canonical ensemble is identified as the physical ensemble for warped AdS$SL(2,\mathbb R)\times U(1)$62/warped BTZ thermodynamics, and the phase structure depends sensitively on whether the bulk action is TMG or NMG. The same work also proposes that boundary modular scrambling modes are related to bulk curvature invariants and syzygies (Ghodrati, 2022).

Taken together, these perturbative and thermodynamic analyses show that warped massive gravity is not only a background classification problem. It is also a framework in which KK spectra, longitudinal-mode dynamics, asymptotic symmetry algebras, exact quasinormal modes, evanescent excitations, and even chaos diagnostics can be computed explicitly. Across its higher-dimensional and three-dimensional realizations, the recurring structural theme is that warping reshapes the infrared behavior of massive gravity without removing the central role of extra polarizations, boundary conditions, and mode mixing.