---
title: 'GLPV Theories: Beyond Horndeski Gravity'
url: https://www.emergentmind.com/topics/gleyzes-langlois-piazza-vernizzi-glpv
type: topic
---

# GLPV Theories: Beyond Horndeski Gravity

Gleyzes–Langlois–Piazza–Vernizzi (GLPV) theories are a class of scalar–tensor modified-gravity models that extend Horndeski theory beyond the requirement of manifestly second-order field equations while still propagating only three physical degrees of freedom, namely two tensor polarizations and one scalar mode [1408.0670]. In unitary gauge they are naturally formulated in ADM variables, where they retain a constrained, degenerate structure that removes the would-be Ostrogradsky mode. In the later taxonomy of scalar–tensor theories, GLPV occupies a distinguished beyond-Horndeski corner of the broader DHOST framework, and it has become a central reference point for questions about Hamiltonian consistency, disformal frame transformations, screening, compact-object solutions, and cosmological observables [1907.02072].

## 1. Emergence from Horndeski and effective field theory

The immediate conceptual background of GLPV is the effective field theory of a single extra scalar degree of freedom in cosmology. In the EFT formulation of dark energy and modified gravity, the starting point is the most general unitary-gauge Lagrangian built from the lapse and the extrinsic and intrinsic curvature of constant-time hypersurfaces. Focusing on linear perturbations, one can identify seven operators that lead to equations of motion containing at most two derivatives, while the homogeneous background is controlled by only three time-dependent coefficients [1304.4840]. Within that framework, Horndeski theory occupies a six-operator subset satisfying a specific relation between the coefficients conventionally denoted \(m_4^2\) and \(\tilde m_4^2\), whereas allowing \(m_4^2\neq \tilde m_4^2\) produces a broader class of single-field theories with second-order linear dynamics. That broader class is the EFT precursor of what is now called GLPV [1304.4840].

This EFT viewpoint is significant because it reframes the criterion for consistency. Horndeski achieves ghost freedom by arranging the covariant equations to be manifestly second order. GLPV instead allows higher-order equations at the Lagrangian level but preserves the physical spectrum through degeneracy and constraints. This distinction became foundational for later DHOST developments, where the absence of extra modes is characterized by degeneracy of the kinetic structure rather than by manifest second-order equations alone [1907.02072].

In the cosmological EFT language, GLPV is also associated with the appearance of the beyond-Horndeski parameter \(\alpha_{\rm H}\). Horndeski corresponds to \(\alpha_{\rm H}=0\), whereas GLPV allows \(\alpha_{\rm H}\neq 0\), together with related time-dependent EFT coefficients such as \(\alpha_{\rm K}\), \(\alpha_{\rm B}\), \(\alpha_{\rm M}\), and \(\alpha_{\rm T}\) [2004.07021]. This parameterization is especially useful for linear cosmological perturbations and for connecting GLPV models to large-scale-structure observables.

## 2. ADM formulation and relation to Horndeski

GLPV theories are most compactly defined in unitary gauge, assuming the scalar field has a timelike gradient and choosing
\[
\phi=\phi(t).
\]
With the ADM decomposition
\[
ds^2=-N^2dt^2+h_{ij}(dx^i+N^i dt)(dx^j+N^j dt),
\]
the GLPV action can be written as
\[
S=\int dt\,d^3x\,\sqrt{h}\sum_{n=2}^{5}L_n,
\]
with
\[
L_2=A_2(t,N),\qquad L_3=A_3(t,N)K,
\]
\[
L_4=A_4(t,N)K_2+B_4(t,N)R,
\]
\[
L_5=A_5(t,N)K_3+B_5(t,N)K_{ij}G^{ij},
\]
where
\[
K_2=K^2-K_{ij}K^{ij},\qquad
K_3=K^3-3K K_{ij}K^{ij}+2K^i{}_j K^j{}_k K^k{}_i
\]
and \(R\), \(G_{ij}\) are the intrinsic three-curvature scalar and Einstein tensor of \(h_{ij}\) [1408.0670]. The functions \(A_2,A_3,A_4,A_5,B_4,B_5\) depend on \(t\) and \(N\), or equivalently on \((\phi,X)\) in covariant language.

Horndeski is recovered by imposing specific relations among these functions. In the quartic sector,
\[
A_4=2X B_{4,X}-B_4,
\]
and in the quintic sector,
\[
A_5=-\frac13 X B_{5,X}.
\]
GLPV is defined by retaining the same ADM operator structure while dropping these relations [1412.6210]. In covariant notation, the quartic beyond-Horndeski deviation is often encoded through
\[
- F_4 X^2 \equiv B_4 + A_4 - 2X B_{4,X},
\]
so that Horndeski corresponds to \(F_4=0\), while GLPV allows \(F_4\neq 0\) [1508.06364].

A standard dimensionless measure of this departure is
\[
\alpha_{\rm H}\equiv \frac{2X B_{4,X}-B_4}{A_4}-1=\frac{F_4 X^2}{A_4}.
\]
Pure Horndeski has \(\alpha_{\rm H}=0\). In GLPV, \(\alpha_{\rm H}\neq 0\) measures the deformation of the scalar–curvature sector away from the Horndeski structure, including the mismatch between the \(A_4\) and \(B_4\) pieces that would combine into the Einstein–Hilbert form in general relativity [1508.06364].

This ADM representation is not merely notational. It makes explicit which combinations of lapse, extrinsic curvature, and intrinsic curvature are retained, which constraints survive, and how GLPV generalizes Horndeski while staying within the same basic unitary-gauge architecture.

## 3. Hamiltonian structure and degree-of-freedom counting

The defining nontrivial claim about GLPV is that it propagates only three physical degrees of freedom even though its equations are generically higher order. A full Hamiltonian proof was given for a simple but nontrivial GLPV subclass with \(A_5=0\) by Lin, Mukohyama, Namba, and Saitou [1408.0670].

In that subclass, the canonical momenta conjugate to \(h_{ij}\) are linear in the extrinsic curvature,
\[
\pi^{ij}=\sqrt{h}\,\big[A_3 h^{ij}+2A_4(h^{ij}K-K^{ij})+B_5 G^{ij}\big],
\]
which permits inversion for \(K_{ij}\) provided \(A_4\neq 0\) [1408.0670]. Because the action contains no \(\dot N\) or \(\dot N^i\), the primary constraints are
\[
\pi_N\approx 0,\qquad \pi_i\approx 0.
\]
Preservation in time generates secondary constraints
\[
\mathcal{H}_i\approx 0,\qquad \mathcal{C}\approx 0,
\]
where \(\mathcal{C}=-\partial_N\mathcal{H}\) plays the role of a lapse-related Hamiltonian constraint [1408.0670].

A key point of the analysis is that the naive momentum constraint must be corrected. The generator of spatial diffeomorphisms is not \(\mathcal{H}_i\) alone but
\[
\mathcal{H}^{\text{(tot)}}_i=\mathcal{H}_i+\pi_N\,\partial_i N.
\]
With this correction, \(\pi_i\) and \(\mathcal{H}^{\text{(tot)}}_i\) form the first-class set associated with spatial diffeomorphisms, while \((\pi_N,\mathcal{C})\) form a second-class pair, provided
\[
\frac{\partial^2\mathcal{H}}{\partial N^2}\neq 0.
\]
The lapse is therefore fixed by constraints rather than generating an extra gauge symmetry [1408.0670].

The degree-of-freedom count is then straightforward. The phase space contains 20 variables:
\[
(N,N^i,h_{ij};\pi_N,\pi_i,\pi^{ij}),
\]
namely 10 configuration variables and 10 conjugate momenta. After gauge fixing the three spatial diffeomorphisms, one has 14 second-class constraints in total, leaving a 6-dimensional physical phase space:
\[
20-14=6.
\]
This corresponds to
\[
\frac{6}{2}=3
\]
configuration-space degrees of freedom, interpreted as two tensor modes and one scalar [1408.0670].

This result is central for the subject. It establishes that “beyond Horndeski” does not mean “beyond three degrees of freedom.” In GLPV, higher derivatives in the equations of motion do not automatically imply an extra Ostrogradsky ghost; the relevant question is whether the Hamiltonian system is degenerate and constrained.

## 4. Disformal transformations, Einstein frame, and DHOST embedding

Disformal transformations play a structural role in GLPV:
\[
g_{\mu\nu}\to \hat g_{\mu\nu}
=\Omega^2(\phi)\,g_{\mu\nu}+\Gamma(\phi,X)\,\partial_\mu\phi\,\partial_\nu\phi.
\]
In unitary gauge, this transformation preserves the GLPV ADM form: a GLPV action is mapped into another GLPV action with transformed coefficients \(\hat A_i,\hat B_i\) [1506.08561]. This closure under disformal maps sharply distinguishes GLPV from a merely ad hoc extension of Horndeski.

For cosmological perturbations on flat FLRW backgrounds, both the curvature perturbation and the tensor perturbation are invariant under the disformal transformation in unitary gauge:
\[
\hat\zeta=\zeta,\qquad \hat\gamma_{ij}=\gamma_{ij}.
\]
The quadratic scalar and tensor actions retain the same form, with transformed coefficients such as
\[
\hat q_s=\frac{q_s}{\Omega^3},\qquad \hat c_s^2=\Omega^2 c_s^2,\qquad
\hat q_t=\frac{q_t}{\Omega^3},\qquad \hat c_t^2=\Omega^2 c_t^2,
\]
and the inflationary scalar and tensor power spectra are invariant up to next-to-leading order in slow-roll [1412.6210].

A particularly useful frame is the Einstein frame defined by choosing \(\Omega\) and \(\Gamma\) so that the tensor quadratic action takes the same form as in general relativity. In that frame,
\[
\hat q_t=\frac{M_{\rm pl}^2}{8\hat{\bar N}},\qquad
\hat c_t^2=\hat{\bar N}^2,
\]
and the effective gravitational potential seen by matter perturbations separates into a GR-like contribution plus terms induced by the disformal coupling [1506.08561]. When the transformed action belongs to a Horndeski subclass, one has
\[
\hat\Psi+\hat\Phi=-\hat\alpha_{\rm H}\,\frac{\hat{\delta N}}{\hat{\bar N}},
\]
so \(\hat\alpha_{\rm H}=0\) implies
\[
\hat\Phi=-\hat\Psi,
\]
i.e. no anisotropic stress in the Einstein frame [1506.08561].

In the broader scalar–tensor classification, GLPV is naturally embedded in DHOST. In the notation of shift-symmetric quadratic-plus-cubic DHOST theories,
\[
F_0=G_2,\quad F_1=G_3,\quad F_2=G_4,\quad F_3=G_5,
\]
with the GLPV identifications
\[
A_1=-A_2=2G_{4X}+Xf_4,\qquad A_3=-A_4=2f_4,\qquad A_5=0,
\]
\[
3B_1=-B_2=\frac32 B_3=G_{5X}+3Xf_5,
\]
\[
-2B_4=B_5=2B_6=-B_7=6f_5,\qquad B_8=B_9=B_{10}=0
\]
[1907.02072]. In this language, pure quartic GLPV belongs to quadratic DHOST, while pure quintic GLPV belongs to cubic DHOST. The black-hole analysis of shift-symmetric DHOST theories indicates that fully general quartic-plus-quintic GLPV with both \(f_4\neq 0\) and \(f_5\neq 0\) cannot simultaneously be degenerate and admit the exact constant-\(X\) Schwarzschild or Schwarzschild–(A)dS solutions considered there; combined with gravitational-wave speed constraints, this pushes viable GLPV models toward purely quartic, effectively quadratic-DHOST sectors [1907.02072].

## 5. Compact objects, singularities, and black-hole solutions

A major nonperturbative issue for GLPV is the behavior of static, spherically symmetric solutions. For a metric
\[
ds^2=-e^{2\Psi(r)}dt^2+e^{2\Phi(r)}dr^2+r^2(d\theta^2+\sin^2\theta\,d\varphi^2)
\]
and a regular central field profile with \(\phi'(0)=0\), one finds that if
\[
\alpha_{\rm H}(r\to 0)\to \alpha_{{\rm H},c}\neq 0,
\]
then the Ricci scalar behaves as
\[
R\simeq -\frac{2\alpha_{\rm H}}{r^2}
\]
near the origin. The resulting geometry is a conical, or solid-angle-deficit, singularity [1508.06364]. The same conclusion holds both in vacuum Schwarzschild-like solutions and inside compact bodies under the simplifying assumption \(\phi'(r)=0\): a nonzero constant central limit of \(\alpha_{\rm H}\) produces a curvature divergence that cannot be removed by tuning integration constants [1508.06364].

This pathology is not generic to every GLPV model. The singularity disappears if
\[
\alpha_{\rm H}\to 0\qquad \text{as}\qquad r\to 0.
\]
A sufficient model-building strategy is to require that the constant-in-\(X\) piece of \(\alpha_{\rm H}\) vanish identically, which in the notation
\[
A_4=-\frac{M^2}{2}F_1(\phi)+f_1(X),\qquad
B_4=\frac{M^2}{2}F_2(\phi)+f_2(X)
\]
amounts to setting
\[
F_1(\phi)=F_2(\phi)
\]
and ensuring that \(f_i(X)\) have no constant term [1508.06364]. Then \(\alpha_{\rm H}\) is proportional to positive powers of \(X\), and since \(X\propto \phi'^2\to 0\) at the center, the singularity is avoided.

An explicit regular example is an extension of the covariant Galileon with a dilatonic coupling,
\[
F(\phi)=e^{-2q\phi/M_{\rm pl}},\qquad
f_1(X)=a_4 X^n,\qquad f_2(X)=b_4 X^n,
\]
with \(F_1=F_2=F\). In that case \(\alpha_{\rm H}\to 0\) as \(r\to 0\), the Ricci scalar remains finite, and a Vainshtein mechanism operates so that general-relativistic behavior is recovered inside a Vainshtein radius [1508.06364].

The central singularity issue is also illuminating by contrast. Beyond-generalized Proca theories reduce to shift-symmetric GLPV in the scalar limit \(A^\mu\to \nabla^\mu\chi\), yet in the full vector theory the temporal vector component regularizes the center and removes the solid-angle-deficit singularity that would appear in the scalar GLPV limit. Near \(r=0\), the vector theory enforces regular expansions
\[
\Psi(r)=\Psi_2 r^2+O(r^4),\qquad
\Phi(r)=\Phi_0+\Phi_2 r^2+O(r^4),
\]
and curvature invariants remain finite [1608.08390]. This comparison suggests that the conical singularity is tied specifically to the scalar-only realization of the beyond-Horndeski quartic structure rather than to the formal presence of an \(f_4\)-type interaction alone.

Black-hole solutions form a separate line of inquiry. Shift-symmetric GLPV theories admit exact stealth Schwarzschild and Schwarzschild–(A)dS solutions with a linearly time-dependent scalar and constant kinetic term in those subclasses satisfying specific algebraic conditions on the coupling functions. However, the DHOST analysis indicates that quintic and full quartic-plus-quintic GLPV sectors are strongly restricted by degeneracy and gravitational-wave constraints, leaving pure quartic GLPV as the most plausible surviving corner for such exact GR-metric black holes [1907.02072].

## 6. Cosmological perturbations, non-Gaussianity, and observational probes

At linear order, GLPV cosmology differs from Horndeski not only through background functions but through the coupling between tensor and scalar sectors. In a simple GLPV dark-energy model, a slight deviation of the tensor propagation speed squared \(c_{\rm t}^2\) from 1 can induce a large modification of the scalar propagation speed squared \(c_{\rm s}^2\) whenever the scalar kinetic energy density is much smaller than the matter density, which is the generic early-time situation in dark-energy models [1503.06539]. A scaling solution with approximately constant \(\rho_X/\rho_m\) avoids this amplification. In the same setup, if the oscillating mode of scalar perturbations is initially suppressed, the anisotropy parameter
\[
\eta=-\frac{\Phi}{\Psi}
\]
can deviate significantly from 1; with generic initial conditions, the deviation of \(c_{\rm t}^2\) from 1 instead produces large oscillations of \(\Psi\) with a frequency related to \(c_{\rm s}^2\) [1503.06539]. These behaviors feed directly into CMB and weak-lensing observables.

In large-scale structure, GLPV effects are often encoded phenomenologically in the second-order density kernel
\[
F_2(\mathbf{k}_1,\mathbf{k}_2,z)
=\kappa_s(z)\,\alpha_s(\mathbf{k}_1,\mathbf{k}_2)
-\frac{2}{7}\lambda_s(z)\,\gamma(\mathbf{k}_1,\mathbf{k}_2),
\]
with
\[
\alpha_s=1+\frac12\frac{\mathbf{k}_1\cdot \mathbf{k}_2}{k_1^2 k_2^2}(k_1^2+k_2^2),
\qquad
\gamma=1-\frac{(\mathbf{k}_1\cdot\mathbf{k}_2)^2}{k_1^2 k_2^2}.
\]
In the parameterization used for weak-lensing higher-order spectra, GR has \(\kappa_s=\lambda_s=1\), Horndeski has \(\kappa_s=1\) with \(\lambda_s\neq 1\), and GLPV allows both \(\kappa_s\neq 1\) and \(\lambda_s\neq 1\) [2004.07021]. This changes the projected cumulant correlators and therefore the skew-spectrum, kurt-spectra, and related higher-order weak-lensing observables. The same modified kernel enters the integrated bispectrum of convergence maps, making squeezed weak-lensing observables a direct probe of GLPV departures from GR [1902.04877].

During inflation, GLPV fits naturally into the unified EFT of inflation. The quadratic action retains the universal form
\[
S_2=\int d^4x\,a^3 Q\left[\dot\zeta^2-\frac{c_s^2}{a^2}(\partial\zeta)^2\right],
\]
with \(Q>0\) and \(c_s^2>0\) required for stability, while the cubic action can be written in a basis of operators controlled by time-dependent EFT coefficients [1804.07741]. A striking result is that GLPV does not generate new scalar bispectrum shapes beyond those already present in Horndeski; rather, it changes the time dependence and amplitudes of the same operator basis. The squeezed-limit consistency relation remains valid even when slow-roll is transiently violated, provided the single-clock assumptions of the EFT remain in force [1804.07741].

An analogous statement holds for primordial tensor non-Gaussianity. In Gao’s unifying ADM framework, the GLPV generalization does not introduce any new tensor cubic operators relative to Horndeski; it only loosens the relations among the coefficients of the same two basic tensor interactions. Four genuinely new tensor cubic operators appear only in theories beyond GLPV [1512.01380]. When the modification of the tensor dispersion relation is small, there is only a single cubic term generating squeezed tensor non-Gaussianity, and it is precisely the term already present in general relativity; the other cubic interactions peak at equilateral configurations [1512.01380].

A more recent development concerns scalar-induced gravitational waves in GLPV theories that are disformally disconnected from Horndeski. In that case, the tensor action contains a new scalar–scalar–tensor interaction controlled by the combination
\[
A_5 C_4^{\rm H}-A_4 C_5^{\rm H},
\]
which vanishes in Horndeski and in GLPV sectors reachable from Horndeski by allowed disformal transformations. This interaction introduces third derivatives in the source for induced tensor modes and enhances the gravitational-wave signal. For a scale-invariant primordial scalar spectrum, the induced gravitational-wave spectral density scales as
\[
\Omega_{\rm GW}(f)\propto f^5,
\]
a much steeper rise than in GR or Horndeski-related theories [2509.05027]. This suggests that early-universe GLPV phases disformally disconnected from Horndeski may be distinguishable through a genuinely new gravitational-wave signature.

Taken together, these results define the modern status of GLPV. The class is broader than Horndeski but narrower than completely general higher-derivative scalar–tensor theories; it is best understood as a constrained, degenerate extension whose consistency is established by Hamiltonian structure, whose relation to Horndeski is organized by disformal transformations, and whose phenomenology is sharply shaped by compact-object regularity, screening, gravitational-wave constraints, and higher-order cosmological observables [1408.0670].

Source: https://www.emergentmind.com/topics/gleyzes-langlois-piazza-vernizzi-glpv