---
title: 'Cubic Galileon: Theory and Applications'
url: https://www.emergentmind.com/topics/cubic-galileon
type: topic
---

# Cubic Galileon: Theory and Applications

Searching arXiv for recent and foundational papers on the cubic Galileon to ground the article in the literature.
The cubic Galileon is a scalar–tensor sector in which a scalar field, usually denoted \(\phi\) or \(\pi\), supplements the Einstein–Hilbert dynamics through a canonical kinetic term and the cubic derivative self-interaction \((\nabla\phi)^2\Box\phi\). In flat space it is associated with the Galilean shift \(\phi\to \phi + c + b_\mu x^\mu\), and in its covariant realizations it belongs to the Horndeski or Kinetic Gravity Braiding subclass, retaining second-order field equations despite higher-derivative interactions in the Lagrangian [1306.1262]. The cubic Galileon has been studied as a mechanism for Vainshtein screening, late-time acceleration, modified galaxy dynamics, scalar hair around compact objects, and as a building block in more elaborate massive-gravity or disformally coupled constructions [1803.01455].

## 1. Formal structure and variants

In four dimensions the Galileon family consists of five independent operators \(\mathcal{L}_1,\dots,\mathcal{L}_5\); the cubic Galileon is the truncation with \(c_4=c_5=0\), so that the scalar sector keeps the tadpole or linear potential term, the canonical kinetic term, and the cubic interaction [1306.1262]. In the standard flat-space normalization, the basic cubic operator is
\[
L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,
\]
and its special status is that varying it yields equations of motion containing no derivatives higher than second order [1302.0840].

Different subliteratures adopt different but equivalent parameterizations. In cosmological applications one often writes a covariant Horndeski action with \(G_2\) and \(G_3\) functions, so that the cubic sector is realized through a \(G_3(\phi,X)\Box\phi\) interaction [1801.01741]. In shift-symmetric black-hole and stellar studies the action is commonly written as
\[
S=\int d^4x\sqrt{-g}\,\Big[\zeta(R-2\Lambda)-\eta(\partial\phi)^2+\gamma\,\Box\phi\,(\partial\phi)^2+\mathcal{L}_m\Big],
\]
which isolates the cubic Galileon term through the constant \(\gamma\) [1911.01669]. A pure cubic sector with no standard kinetic term has also been studied,
\[
S[g,\phi]=\int \left[R+\gamma(\partial\phi)^2\Box\phi\right]\sqrt{-g}\,d^4x,
\]
in the context of rotating black holes [2308.11245].

A broad distinction runs through the literature. One line of work treats the cubic Galileon as a dark-energy or modified-gravity sector at cosmological scales [1306.1262]. Another studies it as a local or astrophysical fifth-force sector, often with asymptotically flat or asymptotically de Sitter boundary conditions around stars and black holes [1212.1394]. A third embeds the cubic interaction into larger frameworks, such as quasi-dilaton massive gravity, where the scalar sector remains Horndeski-like but is coupled to a dRGT graviton potential [2107.09919].

## 2. Screening, effective coupling, and local dynamics

The defining physical mechanism of the cubic Galileon is Vainshtein screening. In the static, spherically symmetric, non-relativistic limit relevant for galaxies, the scalar obeys a modified Poisson-like equation,
\[
\frac{1}{r^2}\frac{d}{dr}\!\left(r^2\frac{d\phi}{dr}\right)
+\frac{1}{2\beta\Lambda}\frac{1}{r^2}\frac{d}{dr}\!\left[r\left(\frac{d\phi}{dr}\right)^2\right]
=8\pi\beta G\rho,
\]
and the ratio \(x(r)\equiv F_5/F_N\) between fifth force and Newtonian force satisfies the algebraic relation
\[
x+\left(\frac{r_V}{r}\right)^3\frac{x^2}{2\beta^2}=2\beta^2,
\qquad
r_V(r)=\left[\frac{GM(r)}{\Lambda^2}\right]^{1/3},
\]
with \(r_V\) the Vainshtein radius [1803.01455]. For \(r\gg r_V\), the fifth force is unscreened and \(x\approx 2\beta^2\); for \(r\ll r_V\), \(x\propto (r/r_V)^{3/2}\ll 1\), so the scalar force is strongly suppressed [1803.01455].

In covariant, time-dependent, spherically symmetric backgrounds the same physics appears through kinetic braiding between the scalar and the metric. For the cubic covariant Galileon in a cosmological background, a linearly time-dependent scalar \(\varphi(t,r)=\varphi(r)+\dot\varphi_c t+\varphi_0\) induces an effective matter–scalar coupling
\[
\alpha_{\rm eff}=\alpha+\frac{k_3}{M^2}\dot\varphi_\infty,
\]
so that even when the bare coupling \(\alpha\) is set to zero, matter is sourced by the scalar through the cosmological time evolution of the background field [1212.1394]. The local solution then depends crucially on asymptotic boundary conditions: Minkowski and de Sitter asymptotics select different branches of the radial solution [1212.1394].

Perturbations around such backgrounds propagate in an effective metric rather than the spacetime metric itself. In the accretion problem on Schwarzschild, the cubic Galileon background defines a sonic horizon outside which perturbations carry a positive-definite energy and whose flux makes the perturbation energy monotonically decrease; steady-state cubic Galileon accretion is therefore linearly stable [1610.00067]. Inside matter, the same kinetic braiding leads to a friction force proportional to the cosmological time derivative of the scalar, producing efficient damping of scalar perturbations [1212.1394].

## 3. Cosmology, late-time acceleration, and structure formation

As a dark-energy model, the cubic Galileon has been analyzed most extensively in a covariant background with \(c_1,c_2,c_3\neq 0\) and \(c_4=c_5=0\) [1306.1262]. Direct \(H(z)\) measurements constrain this model very strongly. For the cubic Galileon with a linear potential, the best-fitting equation of state is indistinguishable from that of a cosmological constant at the \(0.1\%\) level, time variation is allowed only at the few-percent level, and the Galileon energy density can contribute only to about \(10\%\) of the acceleration energy density, the other \(90\%\) being effectively a cosmological constant term [1306.1262]. Although the best-fit \(\chi^2\) is slightly smaller than for \(\Lambda\)CDM, the Bayesian evidence disfavors the cubic Galileon at the “Decisive” level, with odds \(>100:1\) against the model [1306.1262].

A different cosmological direction emphasizes special integrable sectors. For the action
\[
S=\int d^{4}x\sqrt{-g}\left(
\frac12 R-\frac12\partial^\mu\phi\partial_\mu\phi-V(\phi)-\frac12 g(\phi)\partial^\mu\phi\partial_\mu\phi\Box\phi
\right),
\]
the choice
\[
V(\phi)=V_0 e^{-\lambda\phi},
\qquad
g(\phi)=g_0 e^{\lambda\phi},
\]
implies \(V(\phi)g(\phi)=\alpha=\text{const}\) and yields an extra Noether first integral [1703.05860]. This integrable cubic Galileon admits exact power-law solutions and new de Sitter lines, but many critical points are eliminated by ghost or Laplacian instabilities; viable sectors survive only in restricted parameter regions [1703.05860].

At the level of large-scale structure, the cubic Galileon leaves relatively modest linear signatures but more distinctive nonlinear and environmental ones. Weak-lensing analyses of the minimally coupled cubic Galileon with linear potential find that the model can produce stable late-time acceleration while keeping tensor propagation luminal, and that convergence power spectrum and bispectrum differences from \(\Lambda\)CDM or canonical quintessence are potentially detectable with future surveys [1801.01741]. N-body simulations sharpen this picture: the matter power spectrum, halo mass function, galaxy-galaxy weak lensing signal, marked density power spectrum, and counts-in-cells all differ from \(\Lambda\)CDM, with the marked density power spectrum differing by more than \(10\%\), fewer massive halos, and a far higher number of low-density cells in the cubic Galileon model than in \(\Lambda\)CDM [2004.12659]. This suggests that voids and other low-density environments provide unusually strong constraining power [2004.12659].

A separate variant considers a disformal coupling to Standard Model matter,
\[
\tilde g_{\mu\nu}=g_{\mu\nu}+\frac{1}{M_D^4}\partial_\mu\phi\,\partial_\nu\phi,
\]
while keeping the Galileon energy density subdominant [2007.16052]. In that setting, the cubic Galileon is not used as dark energy; instead it is a relic sector whose background contribution is “screening-like” during matter domination because the quadratic and cubic operators nearly cancel in the energy density [2007.16052]. This suppresses the Integrated Sachs–Wolfe signature that otherwise constrains self-accelerating cubic Galileons, and the induced difference between electromagnetic and gravitational propagation is too small to be constrained by GW170817 or by electromagnetic dispersion data [2007.16052].

## 4. Galactic phenomenology and the missing-mass problem

A distinct application treats the cubic Galileon as an alternative to particle dark matter at galactic scales. In the Cubic Galileon Gravity model of galactic rotation curves, matter feels the Newtonian force plus a fifth force
\[
\mathbf{F}_5=-\beta\nabla\phi,
\qquad
F_{\rm tot}=(1+x)F_N,
\]
with \(x(r)\) determined by the Vainshtein-screened algebraic relation above [1803.01455]. The corresponding circular velocity is
\[
v^2(r)=(1+x(r))\frac{GM(r)}{r},
\]
so the Newtonian rotation curve is boosted by a radius-dependent factor [1803.01455].

Applied to the Milky Way using the baryonic SLFC model of Flynn et al. with the dark-matter halo removed, the model fits the observed rotation curve with \(\Lambda=H_0=67.8\ {\rm km\,s^{-1}\,Mpc^{-1}}\) and best-fit coupling \(\beta\approx 20\), reproducing the flat outer profile without invoking a dark matter halo [1803.01455]. Applied to the SPARC sample, after cuts on Hubble type, distance, and distance uncertainty, 28 galaxies are used. Their baryonic masses are reconstructed from the SPARC disk and bulge velocity fields with stellar mass-to-light ratios \(\Upsilon_{\rm disk}\) and \(\Upsilon_{\rm bulge}\), and the best-fit couplings lie in the range
\[
\beta\approx 12-27,
\]
clustered around \(\beta\sim 20\) [1803.01455].

Within that study, the cubic Galileon reproduces rising inner rotation curves that track the baryonic distribution and approximately flat outer rotation curves, and is therefore presented as an alternative to dark matter at galactic scales [1803.01455]. The same parameter range is consistent with Solar System screening: for \(\Lambda=H_0\) and \(\beta\sim 20\), the Solar System Vainshtein radius is estimated as \(r_V\sim 100\ {\rm pc}\), whereas the Solar System size is only \(\sim 100\ {\rm AU}\approx 10^{-3}\ {\rm pc}\), leading to a maximum fractional correction \(x\sim 10^{-6}\) to Newtonian gravity [1803.01455].

The galactic cubic Galileon is often compared with MOND. Both explain rotation curves with baryons plus modified gravity, but the mechanisms are different: MOND introduces an acceleration scale \(a_0\), whereas the cubic Galileon uses a scalar field with derivative self-interactions and Vainshtein screening, with the scale of the modification controlled by \(r_V(r)\), \(\beta\), and \(\Lambda\) [1803.01455].

## 5. Compact objects, black holes, and relativistic stars

The cubic Galileon has a particularly rich strong-gravity sector. In asymptotically de Sitter space, static spherically symmetric stars can be supported by a scalar with linear time dependence,
\[
\phi(t,r)=q t+\int dr\,\frac{\chi(r)}{h(r)},
\]
and the exterior solution matches the corresponding black-hole solution [1911.01669]. Due to the Vainshtein mechanism, the stellar structure is indistinguishable from that of General Relativity with the same central density as long as the stellar radius is shorter than the Vainshtein radius; outside that radius the scalar field is no longer suppressed [1911.01669]. These solutions carry an additional integration constant besides the mass of the star, reflecting scalar hair [1911.01669].

Earlier work on “Black holes in a cubic Galileon universe” found analytic three-dimensional BTZ-like black holes with a non-trivial scalar field that modifies the effective cosmological constant, and numerical four-dimensional asymptotically flat and asymptotically de Sitter solutions [1605.07438]. In the de Sitter case there are three branches associated with the effective cosmological constant, and for two of them the theory admits families of spherically symmetric black holes parameterized by the scalar velocity \(q\) [1605.07438]. The parameter \(q\) represents black-hole primary hair, because it labels inequivalent solutions with the same asymptotic de Sitter behavior [1605.07438].

The pure cubic sector leads to more exotic objects. Fully consistent rotating black holes have been computed numerically in the theory with action
\[
S[g,\phi]=\int \left[R+\gamma(\partial\phi)^2\Box\phi\right]\sqrt{-g}\,d^4x,
\]
using spectral methods in a 3+1 formalism [2308.11245]. These solutions confirm the vanishing of the ADM and Komar masses while retaining finite angular momentum, and the surface gravity is not constant over the horizon, so the zeroth law of black-hole thermodynamics is violated [2308.11245]. The same study finds negative energy density everywhere, hence violation of the weak and dominant energy conditions, and concludes that these solutions do not seem to be a valid alternative to astrophysical black holes, especially with a vanishing mass [2308.11245].

Black-hole accretion furnishes a complementary strong-field test. For a cubic Galileon accreting steadily onto Schwarzschild, the scalar equation admits two branches matched at a sonic horizon, and linear perturbations propagate in an effective metric whose energy flux through that horizon is negative, implying linear stability of the accreting configuration [1610.00067].

## 6. Theoretical tensions, extensions, and status

The cubic Galileon occupies an unusual position in modified-gravity theory. On one side, it is one of the simplest Horndeski sectors with Vainshtein screening and, in minimally coupled form, luminal tensor propagation. This has motivated extensions such as cubic Galileon massive gravity, in which a quasi-dilaton scalar with cubic Galileon kinetic structure couples to a dRGT massive graviton. In that setting the theory admits self-accelerating solutions, fits the Union2 sample of 557 Type Ia supernovae with \(\chi^2_{\rm min}=543.579\), and preserves \(c_T^2=1\) for tensor modes [2107.09919].

On the other side, several theoretical constructions expose sharp limitations. The most direct is supersymmetrization. In globally \(N=1\) supersymmetry, the cubic Galileon scalar must be embedded into a chiral multiplet, so the real scalar is accompanied by a second real scalar. The unique supersymmetric extension that reduces to the standard cubic Galileon when the second scalar is set to zero produces generically third-order equations of motion and therefore an Ostrogradsky ghost [1302.0840]. The theory can still be interpreted as an effective field theory below the ghost scale, but it loses the status of a fundamental ghost-free Galileon [1302.0840].

Global spacetime structure raises separate issues. In Gödel-type metrics, the covariant cubic Galileon with matter content given by a perfect fluid and an electromagnetic field admits consistent space-time homogeneous solutions only under specific coupling relations. For the branch with \(2c_2^2=c_1c_3\), completely causal solutions are excluded; for the shift-symmetric branch \(c_1=0\), a vacuum completely causal solution exists with \(m^2=4\omega^2\) [2009.13242]. This indicates that the cubic Galileon can either preserve or eliminate closed timelike curves depending on the branch and on how the scalar profile enters the geometry [2009.13242].

The overall status is therefore sharply scale- and context-dependent. As a cosmological dark-energy replacement, the cubic Galileon is strongly squeezed toward behavior nearly indistinguishable from \(\Lambda\)CDM [1306.1262]. As a galactic modified-gravity sector, it can fit selected rotation curves without dark matter and remain screened in the Solar System [1803.01455]. As a nonlinear scalar-tensor theory, it supports screened relativistic stars, scalar-haired black holes, stable accretion flows, and highly non-GR rotating black holes [1911.01669]. This suggests that “cubic Galileon” names not a single phenomenological scenario, but a distinctive derivative self-interaction whose viability depends on which sector of the theory is activated and which asymptotic branch is chosen.

Source: https://www.emergentmind.com/topics/cubic-galileon