Cubic Galileon: Theory and Applications

• Cubic Galileon is a scalar-tensor model with derivative self-interactions that preserve Galilean symmetry and yield second-order field equations.
• Its Vainshtein screening mechanism suppresses the scalar-mediated fifth force in high-density regions, aligning modified gravity with solar system tests.
• It impacts cosmic acceleration, offers alternatives to dark matter in galactic rotation curves, and supports exotic compact object solutions.

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 (Bellini et al., 2013). 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 (Chan et al., 2018).

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 (Bellini et al., 2013). 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 (Koehn et al., 2013).

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 (Dinda, 2018). In shift-symmetric black-hole and stellar studies the action is commonly written as

$\pi$0

which isolates the cubic Galileon term through the constant $\pi$1 (Ogawa et al., 2019). A pure cubic sector with no standard kinetic term has also been studied,

$\pi$2

in the context of rotating black holes (Grandclément, 2023).

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 (Bellini et al., 2013). 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 (Babichev et al., 2012). 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 (Aslmarand et al., 2021).

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,

$\pi$3

and the ratio $\pi$4 between fifth force and Newtonian force satisfies the algebraic relation

$\pi$5

with $\pi$6 the Vainshtein radius (Chan et al., 2018). For $\pi$7, the fifth force is unscreened and $\pi$8; for $\pi$9, $(\nabla\phi)^2\Box\phi$0, so the scalar force is strongly suppressed (Chan et al., 2018).

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 $(\nabla\phi)^2\Box\phi$1 induces an effective matter–scalar coupling

$(\nabla\phi)^2\Box\phi$2

so that even when the bare coupling $(\nabla\phi)^2\Box\phi$3 is set to zero, matter is sourced by the scalar through the cosmological time evolution of the background field (Babichev et al., 2012). The local solution then depends crucially on asymptotic boundary conditions: Minkowski and de Sitter asymptotics select different branches of the radial solution (Babichev et al., 2012).

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 (Bergliaffa et al., 2016). 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 (Babichev et al., 2012).

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 $(\nabla\phi)^2\Box\phi$4 and $(\nabla\phi)^2\Box\phi$5 (Bellini et al., 2013). Direct $(\nabla\phi)^2\Box\phi$6 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 $(\nabla\phi)^2\Box\phi$7 level, time variation is allowed only at the few-percent level, and the Galileon energy density can contribute only to about $(\nabla\phi)^2\Box\phi$8 of the acceleration energy density, the other $(\nabla\phi)^2\Box\phi$9 being effectively a cosmological constant term (Bellini et al., 2013). Although the best-fit $\phi\to \phi + c + b_\mu x^\mu$0 is slightly smaller than for $\phi\to \phi + c + b_\mu x^\mu$1CDM, the Bayesian evidence disfavors the cubic Galileon at the “Decisive” level, with odds $\phi\to \phi + c + b_\mu x^\mu$2 against the model (Bellini et al., 2013).

A different cosmological direction emphasizes special integrable sectors. For the action

$\phi\to \phi + c + b_\mu x^\mu$3

the choice

$\phi\to \phi + c + b_\mu x^\mu$4

implies $\phi\to \phi + c + b_\mu x^\mu$5 and yields an extra Noether first integral (Giacomini et al., 2017). 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 (Giacomini et al., 2017).

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 $\phi\to \phi + c + b_\mu x^\mu$6CDM or canonical quintessence are potentially detectable with future surveys (Dinda, 2018). 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 $\phi\to \phi + c + b_\mu x^\mu$7CDM, with the marked density power spectrum differing by more than $\phi\to \phi + c + b_\mu x^\mu$8, fewer massive halos, and a far higher number of low-density cells in the cubic Galileon model than in $\phi\to \phi + c + b_\mu x^\mu$9CDM (Zhang et al., 2020). This suggests that voids and other low-density environments provide unusually strong constraining power (Zhang et al., 2020).

A separate variant considers a disformal coupling to Standard Model matter,

$\mathcal{L}_1,\dots,\mathcal{L}_5$0

while keeping the Galileon energy density subdominant (Lawrence et al., 2020). 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 (Lawrence et al., 2020). 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 (Lawrence et al., 2020).

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

$\mathcal{L}_1,\dots,\mathcal{L}_5$1

with $\mathcal{L}_1,\dots,\mathcal{L}_5$2 determined by the Vainshtein-screened algebraic relation above (Chan et al., 2018). The corresponding circular velocity is

$\mathcal{L}_1,\dots,\mathcal{L}_5$3

so the Newtonian rotation curve is boosted by a radius-dependent factor (Chan et al., 2018).

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 $\mathcal{L}_1,\dots,\mathcal{L}_5$4 and best-fit coupling $\mathcal{L}_1,\dots,\mathcal{L}_5$5, reproducing the flat outer profile without invoking a dark matter halo (Chan et al., 2018). 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 $\mathcal{L}_1,\dots,\mathcal{L}_5$6 and $\mathcal{L}_1,\dots,\mathcal{L}_5$7, and the best-fit couplings lie in the range

$\mathcal{L}_1,\dots,\mathcal{L}_5$8

clustered around $\mathcal{L}_1,\dots,\mathcal{L}_5$9 (Chan et al., 2018).

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 (Chan et al., 2018). The same parameter range is consistent with Solar System screening: for $c_4=c_5=0$0 and $c_4=c_5=0$1, the Solar System Vainshtein radius is estimated as $c_4=c_5=0$2, whereas the Solar System size is only $c_4=c_5=0$3, leading to a maximum fractional correction $c_4=c_5=0$4 to Newtonian gravity (Chan et al., 2018).

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 $c_4=c_5=0$5, whereas the cubic Galileon uses a scalar field with derivative self-interactions and Vainshtein screening, with the scale of the modification controlled by $c_4=c_5=0$6, $c_4=c_5=0$7, and $c_4=c_5=0$8 (Chan et al., 2018).

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,

$c_4=c_5=0$9

and the exterior solution matches the corresponding black-hole solution (Ogawa et al., 2019). 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 (Ogawa et al., 2019). These solutions carry an additional integration constant besides the mass of the star, reflecting scalar hair (Ogawa et al., 2019).

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 (Babichev et al., 2016). 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 $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$0 (Babichev et al., 2016). The parameter $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$1 represents black-hole primary hair, because it labels inequivalent solutions with the same asymptotic de Sitter behavior (Babichev et al., 2016).

The pure cubic sector leads to more exotic objects. Fully consistent rotating black holes have been computed numerically in the theory with action

$L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$2

using spectral methods in a 3+1 formalism (Grandclément, 2023). 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 (Grandclément, 2023). 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 (Grandclément, 2023).

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 (Bergliaffa et al., 2016).

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 $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$3, and preserves $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$4 for tensor modes (Aslmarand et al., 2021).

On the other side, several theoretical constructions expose sharp limitations. The most direct is supersymmetrization. In globally $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$5 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 (Koehn et al., 2013). 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 (Koehn et al., 2013).

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 $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$6, completely causal solutions are excluded; for the shift-symmetric branch $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$7, a vacuum completely causal solution exists with $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$8 (Nascimento et al., 2020). 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 (Nascimento et al., 2020).

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 $L_3=-\frac{1}{2}(\partial\phi)^2\Box\phi,$9CDM (Bellini et al., 2013). As a galactic modified-gravity sector, it can fit selected rotation curves without dark matter and remain screened in the Solar System (Chan et al., 2018). 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 (Ogawa et al., 2019). 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.