David Pirtskhalava
email: david.pirtskhalava@sns.it

Luca Santoni
email: luca.santoni@sns.it

Enrico Trincherini
email: enrico.trincherini@sns.it

Filippo Vernizzi
email: filippo.vernizzi@cea.fr

Weakly Broken Galileon Symmetry

Introduction and summary

Over the years, the effective field theory (EFT) approach has proven itself invaluable for organizing physics at different length scales. With the sole knowledge of the relevant physical degrees of freedom and symmetries -exact or approximate -that govern their dynamics, EFTs can be used to infer model-independent physical predictions. The reason for such universality lies in a fundamental property of the great majority of physical systems, according to which the details of physics at short distances do not leave a qualitative impact on their large-distance characteristics, enabling one to analyze theories with various UV structure in a single framework.

In recent years, the EFT approach has offered new insights in various cosmological contexts, including -but not limited to -inflation [START_REF] Creminelli | Starting the Universe: Stable Violation of the Null Energy Condition and Non-standard Cosmologies[END_REF][START_REF] Cheung | The Effective Field Theory of Inflation[END_REF], dark energy [START_REF] Creminelli | The Effective Theory of Quintessence: the w¡-1 Side Unveiled[END_REF][START_REF] Gubitosi | The Effective Field Theory of Dark Energy[END_REF][START_REF] Bloomfield | Dark energy or modified gravity? An effective field theory approach[END_REF] and the largescale structure [START_REF] Carrasco | The Effective Field Theory of Cosmological Large Scale Structures[END_REF]. For example, the effective theory of inflation is based on the observation that the dynamics of the most general theory of 'single-clock' inflation can be universally captured by an EFT nonlinearly realizing time diffeomorphisms (diffs)1 t → t+ξ 0 (t, x) , with spatial diffs x i → x i + ξ i (t, x) realized linearly. The energy scales of interest are those around the inflationary Hubble rate, E ∼ H, i.e. the frequency at which all correlation functions are measured in the CMB. The spectrum of perturbations consists of the two polarizations of the graviton plus the Goldstone boson of time translation symmetry breaking. That the latter mode has to be present in the spectrum is a direct consequence of the symmetry breaking pattern and has little to do with the exact UV details of the microscopic theory of inflation. This Goldstone boson is what we usually refer to as the adiabatic mode. The CMB indicates that precisely this mode is predominantly responsible for the generation of large-scale structures in the universe. Indeed, any UV theory of inflation that does not lead to extra degrees of freedom around the Hubble energies is equivalent to single-clock inflation and is thus captured by the EFT of Refs. [START_REF] Creminelli | Starting the Universe: Stable Violation of the Null Energy Condition and Non-standard Cosmologies[END_REF][START_REF] Cheung | The Effective Field Theory of Inflation[END_REF].

Given that symmetries define the effective theories, it is of interest to explore the possible symmetries of systems consisting of one or more scalar fields coupled to gravitya typical setup for cosmological model-building. For example, in the context of inflation, effective theories with (approximate) shift symmetry play a key role. In this paper, we study yet another possible -and, as we argue below, necessarily approximate -symmetry of cosmological scalar fields: the invariance under internal galileon transformations φ → φ + b µ x µ .

(1.1)

Theories invariant under (1.1) have appeared in various contexts before. To start with, Eq. (1.1) is a symmetry (up to a total derivative) of the simplest possible quantum field theory: that of a free scalar field. The simplest interacting generalization, nontrivially invariant under (1.1), i.e. containing less than two derivatives per field in the Lagrangian, has appeared in Ref. [START_REF] Luty | Strong interactions and stability in the DGP model[END_REF] in the context of the DGP model [START_REF] Dvali | 4-D gravity on a brane in 5-D Minkowski space[END_REF]. The most general scalar theory with the same property -the galileon -has been proposed in [START_REF] Nicolis | The Galileon as a local modification of gravity[END_REF] and subsequently found [START_REF] De Rham | Generalization of the Fierz-Pauli Action[END_REF] to describe the scalar polarization of the ghost-free dRGT massive graviton [START_REF] De Rham | Resummation of Massive Gravity[END_REF].

Galileon theories are interesting in many ways. Their defining property is that, despite being higher-derivative theories, they have no more than second-order equations of motion, thus describing a single propagating scalar mode. Moreover, as a direct consequence of invariance under (1.1), the coefficients of the leading galileon interactions are not renormalized, at least perturbatively [START_REF] Luty | Strong interactions and stability in the DGP model[END_REF], making any tuning of them stable under quantum corrections.

Despite these interesting properties, the invariance under galileon transformations cannot be exact in nature because the couplings of the scalar to gravity, at least, break it explicitly. We will therefore be interested in characterizing the theories that preserve as much as possible the attractive quantum non-renormalization properties of the exactly invariant case. This naturally leads to the notion of weakly broken invariance under (1.1), which we will define in a precise way in what follows. The simplest sufficiently non-trivial theory with weakly broken galileon invariance (and exact shift symmetry) is of the following form

L = - 1 2 (∂φ) 2 + 1 Λ 3 3 (∂φ) 2 φ + 1 Λ 4 2 (∂φ) 4 . (1.2)
While the first two operators in (1.2) are exactly invariant under (1.1) (up to the boundary terms), the quartic operator is a small breaking as far as Λ 2 Λ 3 . In general we would expect other symmetry breaking operators of the form (∂φ) 2n to be generated by quantum corrections at the scale Λ 2 . However, in this case, a stronger result, the remnant of the non-renormalization properties of the invariant action, holds: all the symmetry breaking operators are generated at a scale that is parametrically higher than Λ 2 . This means in particular that the operator (∂φ) 4 gets only small corrections through loop effects. In the presence of gravity, as we will show in the next section, the weak breaking gives rise to even more non-trivial consequences.

Importantly, the theory (1.2) admits a large generalization that retains both its quantum properties and the second-order equations of motion. Perhaps not surprisingly, we will find that the proper generalization fits into the class of the most general Lorentz-invariant theories with second-order field equations, known as Horndeski theories [START_REF] Horndeski | Second-order scalar-tensor field equations in a four-dimensional space[END_REF]. Since, however, a generic Horndeski theory does not have much to do with the invariance under (1.1), we prefer to refer to the subclass that we study here as 'theories with weakly broken galileon (WBG) invariance'. The purpose of this work is to introduce these theories and to set the stage for their detailed phenomenological studies.

The study of galileon-invariant theories in the inflationary context have been initiated in Ref. [START_REF] Burrage | Galileon inflation[END_REF]. The model, referred to as 'galileon inflation', is based on the particular curvedspace extension of the theory that keeps the property of second-order equations of motion intact -the so-called 'covariant galileon' of Ref. [START_REF] Deffayet | Covariant Galileon[END_REF]. It has been pointed out in [START_REF] Burrage | Galileon inflation[END_REF], that such theories enjoy more freedom compared to ghost/DBI inflation-like models as far as phenomenology is concerned. In particular, they can lead to the possibility of lifting the 'large speed of sound/small non-gaussianities' correspondence, characteristic of the latter models. Our results extend the findings of Ref. [START_REF] Burrage | Galileon inflation[END_REF] in several directions. While it has been assumed in [START_REF] Burrage | Galileon inflation[END_REF] that the covariant galileon is the unique model capable of extending the phenomenological virtues of shift-symmetric theories in a radiatively stable way, we find that there is in fact a wider class of models that can achieve this. Unlike the covariant galileon, the theories we will be interested in do not generically reduce to the standard galileon once gravity is turned off. Nevertheless, neither non-renormalization nor the second-order equations of motion of the galileon need to be given up, leading to the possibility of strong -and quantum-mechanically robust -phenomenological differences from slow-roll inflation. This generalizes the well-studied case of DBI models in an interesting way.

The paper is organized as follows. We start in Secs. 2 and 3 by showing that a (shiftsymmetric) subset of Horndeski theories can be derived solely based on the requirement of WBG invariance at the quantum level, a concept that we define along the way. In Sec. 4, we study general cosmological implications of approximate symmetry under galileon transformations and apply our findings to inflationary, as well as late-time cosmologies in Secs. 5 and 6, respectively. Finally, in Sec. 7 we conclude. Various technical calculations that would overwhelm the main body of the text are collected in the two appendices.

Galileon invariance (with and without gravity)

In this section we study the fate of galileon symmetry upon inclusion of gravity, starting out by reviewing flat-space theories exactly invariant under (1.1). If internal galileon invariance is realized in nature, however, it has to be approximate: the coupling to gravity unavoidably breaks this symmetry and, at the same time, defines the sense in which the breaking can be considered weak. Indeed, as we will see below, out of the infinite number of possible ways one can couple the galileon to gravity, there is a unique set of non-minimal couplings that can be qualified as being 'more invariant' under galileon transformations than all the rest. The resulting theory is the covariant galileon of Ref. [START_REF] Deffayet | Covariant Galileon[END_REF]. On the other hand, since the galileon symmetry has to be weakly broken by the couplings to gravity, one should in principle include in the effective theory symmetry-breaking operators of e.g. the form (∂φ) 2n . The latter will be generated by quantum loops with Wilson coefficients bounded from above, so as to be consistent with the approximate symmetries at hand. As we will see, requiring invariance under (1.1) to be only weakly broken will naturally lead to a particular sub-class of Horndeski theories, which is however much more general than just the covariant galileon.

Flat space galileons

Consider a trivial, free theory of a scalar φ. As emphasized above, in addition to more familiar symmetries (such as the ones under constant shifts or conformal transformations), this theory possesses an extra invariance under internal galileon transformations (1.1). The latter leaves the action invariant only up to a boundary term. As shown in [START_REF] Nicolis | The Galileon as a local modification of gravity[END_REF], apart from the free theory -and even a more trivial tadpole term, there are exactly three more interaction terms in four spacetime dimensions that share invariance under internal galileon transformations in a non-trivial way. The corresponding theory can be written as

L = (∂φ) 2 + 5 I=3 c a Λ 3(I-2) 3 L i , (2.1)
with the three interaction terms given by

L 3 = (∂φ) 2 [Φ] , (2.2)
L 4 = (∂φ) 2 [Φ] 2 -[Φ 2 ] , (2.3)
L 5 = (∂φ) 2 [Φ] 3 -3[Φ][Φ 2 ] + 2[Φ 3 ] , (2.4)
where we denote

[Φ] ≡ φ, [Φ 2 ] ≡ ∂ µ ∂ ν φ∂ ν ∂ µ φ, etc.
In addition to being invariant under (1.1), galileon theories share another special property: the associated scalar equations of motion are second order, both in time and in space, despite the higher-derivative interactions in Eqs. (2.2)- (2.4). This guarantees that there are no extra propagating degrees of freedom hidden in φ. Moreover, the particular structure of the Lagrangian (2.1) results in a nonrenormalization theorem that allows to radiatively generate only terms trivially invariant under (1.1), i.e. terms with at least two derivatives acting on φ. In particular, the coefficients c I of galileon interactions are not renormalized by quantum loops in the presence of exact galileon invariance [START_REF] Luty | Strong interactions and stability in the DGP model[END_REF]. In theories of modified gravity with matter, galileons couple to the rest of the degrees of freedom via symmetry-breaking Planck-suppressed operators2 , usually making the renormalization of c I parametrically weak.

Coupling to gravity

It is generally impossible to couple the galileon to gravity3 , while keeping invariance under galileon transformations, or their curved-space generalization, intact. Thus, galileon symmetry is expected to be broken even in the purely scalar sector -i.e. if one sets the metric perturbation to zero 'by hand' -due to loop-generated operators not invariant under (1.1) (these operators of course have to be suppressed by at least one power of the Planck mass).

∂φ ∂ m φ, m ≥ 2 ∂ n h
Figure 1: For every given vertex, a solid straight line corresponds to a possible external leg resulting in a scalar with at most one derivative acting on it. A dashed line corresponds to an external scalar with at least two derivatives, while a wiggly line denotes an external graviton.

For the sake of concreteness, let us concentrate on the operators of the form (∂φ) 2n . After all, it is the absence, or smallness, of these operators that makes the galileon special compared to a generic shift-symmetric theory 4 . These operators can be generated from loops with the appropriate number of insertions of symmetry-breaking vertices, that include at least one graviton line. Interpreting Λ 3 -the smallest scale by which interactions are suppressed in the effective theory5 -as the genuine cutoff, any loop-generated operator can be written as

(∂φ) 2n Λ 4(n-1) k,n , Λ k,n ≡ M k Pl Λ 4n-k-4 3 1 4n-4 , (2.5)
where k denotes some positive integer and we have omitted factors of 16π 2 for simplicity. For fixed k and n sufficiently large, the scale suppressing a given single-derivative operator, Λ k,n , can in principle be arbitrarily close to Λ 3 . If this were true, the resulting theory would by no means be considered as a theory with WBG invariance, since symmetry-breaking operators would be order-one in the units of the cutoff.

We wish to show here that: i) this in fact does not happen for the galileon, even if minimally coupled to gravity; ii) for a very particular, non-minimal coupling, the scale suppressing the symmetry-breaking operators can be made parametrically higher than the analogous scale characterizing all the other ways of introducing gravity into the system. For the latter theories, there is a well-defined separation between the scale suppressing the invariant galileon interactions and the quantum-mechanically generated single-derivative operators: while the former are suppressed by powers of Λ 3 , the scale suppressing the latter is parametrically higher, asymptoting to

Λ 2 ≡ (M Pl Λ 3 3 ) 1/4 (2.6)
for large n in (2.5). Perhaps not surprisingly, we will find that the special type of couplings to gravity required to suppress symmetry-breaking operators is of the Horndeski classthe unique curved-space extension of the (generalized) galileons that leads to second-order equations of motion both for the scalar and the metric.

Let us now prove all of the above-described statements, starting with a few definitions. Consider an arbitrary vertex of the form (∂φ) k (∂ m φ) n ∂ l h p , with m ≥ 2. In terms of Feynman diagrams, we will adopt the convention according to which a solid straight line corresponds to an external leg coming out of this vertex, which results in a scalar with one derivative acting on it. A dashed line corresponds to an external scalar with at least two derivatives, while a wiggly line denotes an exernal graviton, see Fig. 1. For example, a k = 2, n = 1, p = 1 vertex would correspond to the first diagram on the upper line of Fig. 2. In certain cases, the number of solid lines can be less than the number of factors of ∂φ in a Lagrangian interaction term. For example, the special structure of the pure galileon interactions makes them equivalent to vertices with only dashed lines, even though on average these terms contain less than two derivatives per field. This is because a galileon vertex can only lead to external states with at least two spacetime derivatives -the fact that lies at the heart of the non-renormalization theorem associated with these theories [START_REF] Luty | Strong interactions and stability in the DGP model[END_REF]. We now wish to show that there exist more redundant vertices of this kind in the suitable non-minimal extension of the theory to curved space.

Let us have a look at all possible vertices with a single graviton line in the minimally coupled galileon theory, obtained from (2.1) by replacing all derivatives with covariant ones (∂ µ → ∇ µ ). A straightforward inspection yields that there are six of such vertices, shown in Fig. 2. Vertices with three solid lines can in principle arise from picking up a factor of ∂h from covariant derivatives6 , i.e. ∇ 2 φ ∼ ∂ 2 φ + ∂h∂φ, where we denote by ∇ 2 any product of two covariant derivatives applied on φ. This means that, in the minimally coupled galileon theory, the smallest scale by which the operators of the form (2.5) with n = 3a, k = 2a are suppressed, is7 ∼ (M Pl Λ 5

3 ) 1/6 . The latter suppression scale is still too small to be consistent with the generic definition of WBG invariance, which we will introduce and discuss extensively in what follows.

We will now show that it is possible to modify the theory by non-minimal couplings to gravity that result in the elimination of vertices with three solid lines, leaving just a factor of (∂φ) 2 per graviton line. Insertion of these vertices into a generic loop diagram can only lead to symmetry-breaking operators of the form (2.5) suppressed by at least the scale Λ 2 = (M Pl Λ 3

3 ) 1/4 , parametrically higher than (M Pl Λ 5 3 ) 1/6 . Therefore, there is a welldefined sense in which the resulting theory is 'more invariant' under the internal galileon transformations than a generic curved-space extension of the galileon. This defines what we mean by 'theories with WBG invariance' throughout the present work and we will sharpen this definition in the next section, where non-trivial vacua of the resultant theories are considered.

For the cubic covariant galileon, it is easy to show that the first vertex from the lower line of Fig. 2 is absent: the only way to generate it is by picking up a metric perturbation from the covariant laplacian. This yields a term of the form (∂φ) 3 ∂h. Putting the derivative Figure 2: Single-graviton vertices, which can contribute two or three external scalars with one derivative ((∂φ) 2 or (∂φ) 3 ) to a 1PI vertex. The vertices from the second line exactly cancel for the unique (non-minimal) curved-space extension of the galileon that retains second-order equations of motion.

from h on the rest of the fields via partial integration however makes it evident that the corresponding vertex can only have two solid lines, but not three. The case of quartic and quintic galileons is more non-trivial, but straightforward; we show in Appendix A that vertices with three solid lines and one graviton, as well as five solid lines and two gravitons can be removed by suitably adding non-minimal couplings to gravity. The resultant theory is the covariant galileon of Ref. [START_REF] Deffayet | Covariant Galileon[END_REF] -characterized, as a bonus, by second-order equations of motion both for the scalar and the metric.

In summary, the 'most symmetric' generalization of the galileon coupled to gravity consists of the following operators

1 Λ 3 3 L 3 → √ -g 1 Λ 3 3 L min 3 , (2.7)
1 Λ 6 3 L 4 → √ -g 1 Λ 6 3 (∂φ) 4 R -4 L min 4 , (2.8)
1 Λ 9 3 L 5 → √ -g 1 Λ 9 3 (∂φ) 4 G µν ∇ µ ∇ ν φ + 2 3 L min 5 , (2.9)
where L min I denote the galileons (2.2)-(2.4) minimally coupled to gravity. The structure of the full effective theory is such that every pair of external scalars with no more than one derivative on each unavoidably comes with a suppression of at least one power of the Planck scale. With this in mind, quantum-mechanically generated operators of the form (∂φ) 2n can be estimated simply on dimensional grounds to be at most of the following magnitude

(∂φ) 2n M n Pl Λ 3n-4 3 .
(2.10)

The scales (M n Pl Λ 3n-4

3

) 1/(4n-4) suppressing such operators approach Λ 2 only for asymptotically large n, otherwise being parametrically larger.

In the following sections, we will argue that the statement of approximate galileon invariance can be viewed as a statement about non-trivial classical vacua, generically present in the shift-invariant theories at hand. If these vacua are to be insensitive to UV physics, the symmetry-breaking operators in the effective theory can be at most of order (2.10) in magnitude. In such cases we will say that the galileon invariance is broken only weakly by couplings to gravity.

Theories with WBG invariance

Perhaps, the most important lesson that one can draw from the discussion of the previous section is that the covariant galileon is in fact not the most general theory enjoying the abovedescribed quantum properties. Consider the effective theory with the leading operators given by the following subclass of the Horndeski terms

L WBG 2 = Λ 4 2 G 2 (X) , (3.1)
L WBG 3 = Λ 4 2 Λ 3 3 G 3 (X)[Φ] , (3.2)
L WBG 4 = Λ 8 2 Λ 6 3 G 4 (X)R + 2 Λ 4 2 Λ 6 3 G 4X (X) [Φ] 2 -[Φ 2 ] , (3.3)
L WBG 5 = Λ 8 2 Λ 9 3 G 5 (X)G µν Φ µν - Λ 4 2 3Λ 9 3 G 5X (X) [Φ] 3 -3[Φ][Φ 2 ] + 2[Φ 3 ] , (3.4)
where we now extend the notation employed in Eqs. (2.2)-(2.4) by replacing the partial derivatives with covariant ones, i.e.

[Φ] ≡ g µν ∇ µ ∇ ν φ, [Φ 2 ] ≡ ∇ µ ∇ ν φ∇ ν ∇ µ φ, etc. Moreover, G I are arbitrary dimensionless functions of the dimensionless variable X ≡ - 1 Λ 4 2 g µν ∂ µ φ∂ ν φ , (3.5)
and the subscript 'X' means differentiation with respect to this variable 8 . We will then assume that the functions G I are defined through the Taylor expansion

G I = c (0) I -c (1)
I X + c (2)
I X 2 + • • • = c (0) I + c (1) I (∂φ) 2 Λ 4 2 + c (2) I (∂φ) 4 Λ 8 2 + . . . , (3.6)
8

The particular form of the interactions (3.1)-(3.4) (the relative coefficients between two operators in L WBG 4 or L WBG 5 , for example) may appear tuned, and this is sometimes presented as an unfortunate feature in the literature. We stress that this 'tuning', motivated by restoring unitarity in the theory, is in fact natural and stable under quantum corrections. This is very similar to why we work with a gauge-invariant kinetic term -Tr F µν F µν in theories with massive spin-1 particles, despite of there being no gauge invariance in the presence of the mass term.

where c

(n) I are dimensionless, order-one coefficients. Note that by setting all the coefficients c

(n) I to zero except c (1) 2 , c (1) 3 , c (2)
4 and c

(2) 5 we recover the covariant galileon of Eqs. (2.7)-(2.9). As in the example of the introductory section, the theory at hand is characterized by two scales 9 : Λ 2 and Λ 3 Λ 2 . From the EFT standpoint that we are adopting here, all Wilson coefficients g (n) I are measured in the units of the cutoff Λ 3 . The theory (3.1)-(3.4) is the one for which a well-defined class of (canonically normalized) operators -those that respect the galileon symmetry φ → φ + b µ x µ -have order-one coefficients,

g (n) I ∼ O(1) , (3.7)
while all the others are parametrically suppressed with respect to the coefficients c

(n) I by positive powers of the ratio Λ 3 /M Pl . As a consequence, one can define a scaling limit

M Pl → ∞, Λ 2 → ∞, Λ 3 = finite , (3.8)
in which the system recovers the exact internal galileon symmetry.

A comment about the terms c

(1)

4 and c

(1)

5 is in order here 10 . These terms are absent in the covariant galileon, which starts with G 4 , G 5 ∼ (∂φ) 4 . Thus, one can wonder whether they modify our conclusions regarding quantum corrections. By expanding L 4 and L 5 at linear order in metric perturbations, it is straightforward to verify that c 5 do lead to new operators at the scale Λ 3 . However, these operators are just a scalar-tensor generalization of the galileon familiar from the decoupling limit of dRGT massive gravity [START_REF] De Rham | Generalization of the Fierz-Pauli Action[END_REF],

L WBG 4 ∼ h µν ε µαρλ ε νβσλ ∂ α ∂ β φ∂ ρ ∂ σ φ + . . . , L WBG 5 ∼ h µν ε µαρλ ε νβσδ ∂ α ∂ β φ∂ ρ ∂ σ φ∂ λ ∂ δ φ + . . . . (3.9)
The latter interactions obey the same non-renormalization theorem as the galileon [START_REF] De Rham | Nonrenormalization and naturalness in a class of scalar-tensor theories[END_REF], forbidding asymptotic states with less than two derivatives acting on them. In particular, this means that the counting of M -1 Pl factors in loop diagrams of the previous section goes through unaltered: each factor of (∂φ) 2 arising from a generic 1PI vertex still comes with a suppression of at least a factor of M -1

Pl , and our analysis of quantum corrections performed for the covariant galileon also applies for the generalized theory (3.1)-(3.4). 11 Therefore, for simplicity we will disregard these terms and assume that G 4 and G 5 start at least quadratic in X. Indeed, our discussion of the previous chapter guarantees that, once tuned to zero, G 4X and G 5X do not receive any appreciable quantum corrections as a result of the approximate galileon symmetry. 9 The Planck scale arises here from L WBG 5 give inconsequential total derivatives and can be disregarded. 11 That the terms (3.9) can be rewritten as a certain limit of non-minimally coupled scalar-tensor theory has been noticed in [START_REF] Chkareuli | Vainshtein Mechanism In Λ 3 -Theories[END_REF] and the cosmology of that theory has been studied in [START_REF] De Rham | Cosmology of the Galileon from Massive Gravity[END_REF].

Notice that the theories that we propose do not reduce to the galileon for general functions G I , even if one switches off gravity. Rather, they describe a certain generalization thereof, which still has no more than second-order equations of motion. The symmetry (1.1) is thus broken even when one sets h µν = 0 (note that this is different from taking the decoupling limit (3.8)). Nevertheless, the corresponding breaking of invariance under (1.1) is in a well-defined sense 'not stronger' than the one already present once gravity is turned on. In particular, it follows from the very construction of the operators in (2.7)-(2.9) that each pair of factors of ∂φ still comes with a suppression of at least a factor of M Pl in the full effective theory. Thus, the Lagrangians (3.1)-(3.4) define the most general scalar-tensor theory consistent with our definition of WBG invariance.

The peculiar quantum properties characteristic of the theories at hand put loop corrections under complete control for a broad class of physical backgrounds. This can be established by noting that loop-generated operators always have at least one extra factor of M -1

Pl , compared to their 'tree-level' counterparts. Indeed, working in the units in which Λ 3 = 1, each factor of Λ 4 2 becomes equivalent to M Pl , and the non-minimal terms in (3.1)-(3.4) can be schematically written as

M Pl R(∂φ) 2 1 + (∂φ) 2 M Pl + . . . (∇ 2 φ) n , (3.10)
while the similar loop-generated terms can never have the M Pl -enhancement. Likewise, the symmetry-breaking scalar operators generated quantum-mechanically are bound to have the following schematic form

(∂φ) 2n M n Pl ∇ 2 φ m , (3.11)
and the analogous terms in the original action are all enhanced by at least a factor of M Pl compared to this expression. Restoring the factors of Λ 3 , we conclude that if the field expectation values are such that

|X| ∼ < 1, |Y | ∼ < 1 , (3.12)
where X is defined in Eq. (3.5) and

Y ≡ ∇ 2 φ Λ 3 3 , (3.13)
the loop corrections are negligible: they can never compete with the operators from the classical action (3.1)-(3.4) and all predictions based on the latter can be trusted in the full quantum theory. Moreover, not only are the magnitudes of the various operators technically natural, but any possible tuning of these can only be spoiled by corrections of higher order in Λ 3 /M Pl . This generalizes the galileon non-renormalization theorem in the presence of gravity.

On the other hand, if the conditions (3.12) are not satisfied, the background becomes too strong to be trusted: the magnitude of loop corrections becomes pumped up beyond that of the leading operators and galileon invariance alone is useless in organizing the low-energy expansion. One could only trust the classical theory beyond that point if there is an extra structure, leading to a resummation of the series (3.6). An example would be, e.g., the DBI theories [START_REF] Silverstein | Scalar speed limits and cosmology: Acceleration from D-cceleration[END_REF][START_REF] Alishahiha | DBI in the sky[END_REF] (as well as generalizations thereof [START_REF] Hinterbichler | DBI Realizations of the Pseudo-Conformal Universe and Galilean Genesis Scenarios[END_REF][START_REF] Hinterbichler | Dirac-Born-Infeld Genesis: An Improved Violation of the Null Energy Condition[END_REF]), where this happens due to non-linearly realized higher-dimensional spacetime symmetries. However, here we want to remain as general as possible and do not assume any extra symmetry beyond the (weakly broken) galileon one. In the latter case, the requirement of quantum stability expressed by Eq. (3.12) places a strong upper bound on how large the predictions for various physical quantities can be.

De Sitter vacua

Perhaps, the best way to illustrate the physical consequences of the discussion of the previous section is to resort to a concrete example. To this end, we wish to consider the following theory

S = d 4 x √ -g 1 2 M 2 Pl R - 1 2 (∂φ) 2 -V (φ) + 5 I=2 L WBG I + . . . , (4.1)
where L WBG I are the operators defined in Eqs. (3.1)-(3.4) of the previous chapter. We have extracted and canonically normalized the scalar and graviton kinetic terms, so that G 2 , G 4 and G 5 are assumed to start at least quadratic in X (see the discussion below Eq. (3.9)), while G 3 can have a linear piece. In general, one can allow for a small potential,

V (φ) = -λ 3 φ + 1 2 m 2 φ φ 2 + . . . , (4.2)
where the parameters λ, m φ , etc. can be (naturally) arbitrarily small, since they are the only ones that break the scalar shift symmetry (which otherwise is exact even on curved space).

In the next section we will discuss the implications of a weakly broken galileon symmetry for inflation. We will thus be interested in background solutions given by an approximate de Sitter space with the Hubble parameter H. This section is devoted to a preliminary study of such solutions.

For the theory (4.1), the scalar equation of motion on a flat FRW background reads

1 a 3 d dt 2a 3 φF (X, Z) = - dV dφ , (4.3)
where the function F is given by the following expression12

F (X, Z) = 1 2 + G 2X -3ZG 3X + 6Z 2 G 4X X + 2G 4XX + Z 3 3 G 5X X + 2G 5XX , (4.4)
with the variable Z defined so as to roughly coincide in the limit φ H φ with the background value of Y defined in Eq. (3.13),

Z ≡ H φ Λ 3 3 . (4.5)
Moreover, the two Friedmann equations can be written in the following way

3M 2 Pl H 2 = V + Λ 4 2 X 1 2 - G 2 X + 2G 2X -6ZG 3X -6Z 2 G 4 X 2 -4 G 4X X -4G 4XX + 2Z 3 5 G 5X X + 2G 5XX , (4.6)
M 2 Pl Ḣ = - Λ 4 2 XF + M Pl φ(XG 3X -4ZG 4X -8ZXG 4XX -3Z 2 G 5X -2Z 2 XG 5XX ) 1 + 2G 4 -4XG 4X -2ZXG 5X . (4.7)
By choosing appropriate combinations of the coefficients c

(n)
I and in the absence of the potential, one can check that there exist exact linear solutions φ 0 ∝ t to (4.3) with F (X 0 , Z 0 ) = 0, sourcing an exact de Sitter space. Upon turning on a small enough potential for φ, these backgrounds can slightly deviate from de Sitter, Ḣ H 2 . Moreover, if both X 0 and Z 0 are of order one on these solutions, the background curvature is of the order H 2 ∼ Λ 4 2 /M 2 Pl and all terms involving φ in (4.1) contribute equally to the energy density.

More generally, once shift-symmetry is broken, the parameters X and Z measure the departure of quasi de Sitter solutions from slow-roll inflation. In particular, for X 1 and Z

1 one reproduces the standard potential-dominated slow-roll regime. In the opposite limit and for small enough V (φ), the acceleration can be supported by the kinetic part of the action. At the same time, we have argued below Eq. (3.12) that the same two parameters control the magnitude of quantum corrections. Indeed, in the X ∼ < 1 and Z ∼ < 1 regime loop corrections are fully under control, which can be seen by estimating typical loop-generated terms, e.g. of the form

(∇ 2 φ) n Λ 3n-4 3 ∼ Z n Λ 4 3 , (∂φ) 2n M n Pl Λ 3n-4 3 ∼ X n Λ 4 3 . (4.8)
For X and Z less than or of order one, these are parametrically suppressed with respect to the background energy density contributed by the classical Lagrangian, which, as follows from Eq. (4.6), generically satisfies 3M 2 Pl H 2 ∼ > XΛ 4 2 . Saturating either one, or both of the bounds in Eq. (3.12) leads to a non-linear regime where various next-to-leading order terms in the derivative expansion become important, while quantum corrections are still under control. This case will be of our prime interest in the rest of the paper. In particular, in the next section we will explore the properties of the corresponding backgrounds in the context of inflationary physics.

Inflation

The effective field theory for perturbations

The consequences of approximate galileon symmetry for the physics of the early universe can be conveniently studied in the effective field theory of inflation (EFTI) framework [START_REF] Creminelli | Starting the Universe: Stable Violation of the Null Energy Condition and Non-standard Cosmologies[END_REF][START_REF] Cheung | The Effective Field Theory of Inflation[END_REF]. Operationally, the EFT of interest is conveniently constructed in the unitary gauge, i.e. by choosing the special coordinate system where the time slicing coincides with the uniform inflaton hypersurfaces. The spacetime metric can be decomposed into the ADM variables,

ds 2 = -N 2 dt 2 + γ ij N i dt + dx i N j dt + dx j , (5.1)
where N is the laspe and N i the shift, while γ ij is the 3d metric.

The 'breaking' of time diffs allows to include objects that are not invariant under the full 4d diffs, but preserve the reparametrizations of the 3d coordinates on equal-time hypersurfaces. Thus, the possible building blocks of the EFT of interest are homogeneous functions of time, intrinsic and extrinsic curvatures of equal-time hypersurfaces, respectively

R ij and K ij , and temporal components of tensors with upper indices such as N = 1/ -g 00 [START_REF] Creminelli | Starting the Universe: Stable Violation of the Null Energy Condition and Non-standard Cosmologies[END_REF][START_REF] Cheung | The Effective Field Theory of Inflation[END_REF]. To construct the action up to quadratic order in metric fluctuations, we perturb these quantities around a flat FRW background, ds 2 = -dt 2 + a 2 (t)d x 2 , by writing δN ≡ N -1,

δK ij ≡ K ij -Hγ ij (K ≡ K i i ), while (3)
R ij is already at least linear in perturbations. In terms of these quantities, the action reads [START_REF] Creminelli | Starting the Universe: Stable Violation of the Null Energy Condition and Non-standard Cosmologies[END_REF][START_REF] Cheung | The Effective Field Theory of Inflation[END_REF][START_REF] Gleyzes | Essential Building Blocks of Dark Energy[END_REF][START_REF] Piazza | Effective Field Theory of Cosmological Perturbations[END_REF]]

S = d 4 x √ γ N M 2 Pl 2 f (t) (3) R + K ij K ij -K 2 -2 ḟ (t) K N + c(t) N 2 -Λ(t) + M 4 (t) 2 δN 2 -M 3 (t)δKδN - M 2 (t) 2 δK 2 -δK ij δK ij + m2 (t) 2 (3) R δN - M 2 (t) 2 δK 2 + δK ij δK ij + m 1 (t) (3) RδK + . . . . (5.2)
The ellipses stand for Lagrangian terms constructed of space or time derivatives of N , (3) R ij and K ij . Moreover, we have omitted the operator (3) R ij δK ij , since at the quadratic order it can be rewritten as a combination of other operators in (5.2), see [START_REF] Gleyzes | Essential Building Blocks of Dark Energy[END_REF].

As far as inflation is concerned, one can always remove the time-dependent function f by means of a conformal transformation of the metric [START_REF] Cheung | The Effective Field Theory of Inflation[END_REF]. Thus, in this section we will assume f = 1, in which case the time-dependent coefficients c and V in the first line -those that lead to tadpoles -are uniquely determined in terms of the background equations of motion,

c = -M 2 Pl Ḣ , Λ = M 2 Pl (3H 2 + Ḣ) . (5.
3)

The rest of the action contains operators whose coefficients are in principle unconstrained. Depending on which subset of these coefficients is nonzero, the above effective theory describes the linear perturbations in different (and all possible) models of single-field inflation.

Given that a generic model of inflation can be defined via the set of EFT coefficients M 4 (t), M 3 (t), M 2 (t), . . . , it is important to understand what are, if any, the possible hierarchies among these. Indeed, the basic properties of the inflationary background and perturbations such as stability, power spectra, non-gaussianities, etc., crucially depend on the relative magnitudes of these terms. In the case that the problem involves widely separated scales (such as, for instance, the Hubble scale and the UV cutoff of the theory), it is not a priori clear what the quantum-mechanically stable values of these Wilson coefficients are. In particular, treating them as arbitrary parameters can easily lead to results that require tuning, once quantum corrections are taken into account.

The easiest way to arrive at various hierarchies in quantum field theory is using symmetries -exact or approximate. For example, chiral symmetry, although inexact, is the central reason behind the technical naturalness of fermion masses in the standard model of particle physics. In the same spirit, one can use various symmetries in the effective theory (5.2) to constrain the EFT coefficients in different models of inflation. One obvious possible symmetry are constant time translations, t → t + c, which would constrain all coefficients (including the Hubble rate) in the EFT to be time-independent, enforcing the background solution to be an exact de Sitter space. Another symmetry one can impose are arbitrary space-independent reparametrizations of time, t → f (t); see, e.g., [START_REF] Blas | Models of non-relativistic quantum gravity: The Good, the bad and the healthy[END_REF]. In the exact symmetry limit, this would require M, M and m to vanish, leading to interesting theories in one corner of the inflationary theory space [START_REF] Creminelli | Khronon inflation[END_REF]. On the other hand, even if time reparametrizations are not an exact symmetry, their weak breaking (appropriately defined in what follows) would guarantee that any small values of these coefficients are technically natural.

Building on the insights of the previous sections, we wish to investigate the consequences of yet another possible approximate invariance of the action (5.2) -that induced by internal galileon transformations realized on the foliation scalar as in (1.1). If this scalar acquires a linear profile φ ∝ t driving a de Sitter background (similarly to what we discussed in the previous section), the galileon symmetry (1.1) will have the following manifestation in the unitary gauge

t → t + bµ x µ , (5.4)
where bµ and c are a constant four-vector and a scalar, respectively. For more complicated φprofiles, the associated unitary gauge invariance will have a more complicated form. However, in practice it always makes life easier to restore the foliation scalar in the effective action (5.2) in order to explore its approximate symmetries. We will take this route in what follows.

'Ordinary' effective theories

Before turning to the case of models with galileon symmetry, it is instructive to recall what happens in more standard scalar effective theories. There, one starts with an action organized in derivative expansion containing a single scale Λ -the EFT cutoff,

S = S EH + d 4 x √ -g Λ 4 c G 2 (X) -V (φ) + ( φ) 2 Λ 2 c + . . . , (5.5)
where, as before, G 2 is a dimensionless function of the variable X = -(∂φ) 2 /Λ 4 c , parametrizing the leading derivative effects.

As discussed in the previous section, for a flat enough potential the scalar background profile can be approximated by a linear function of time φ 0 = c(t)Λ 2 c t, where c(t) is a slowly varying function. This is true both for slow-roll inflation, where d ln φ/dN ∼ d ln Ḣ/dN is of the order of the slow-roll parameters, and for models such as the ghost condensate, where the field 'velocity' is explicitly constant, i.e. c = const. at the leading order [START_REF] Arkani-Hamed | Ghost inflation[END_REF]. For definiteness, let us concentrate on the latter class of models and assume that c ∼ 1 solves G 2X (c 2 ) = 0. On such solutions, the Hubble rate can be estimated as M 2 Pl H 2 ∼ Λ 4 c and is completely determined by the leading G 2 (X) piece, higher-derivative operators being unimportant for the backgrounds at hand. For example, the quadratic operator shown in Eq. (5.5) can be estimated as ( φ 0 ) 2 /Λ 2 c ∼ Λ 4 c (H/M Pl ) and thus it is negligible for H M Pl , as required for a consistent classical description.

Restricting to the unitary gauge, δφ(t, x) = 0, one straightforwardly finds that all EFT coefficients are determined by the cutoff of the theory

M 4 ∼ Λ 4 c , M 3 ∼ Λ 3 c , M 2 ∼ Λ 2 c , M 2 ∼ Λ 2 c , etc. (5.6)
In this case, the dynamics of small perturbations is fully dominated by the only zeroderivative quadratic operator in the effective theory -δN 2 -and most of the phenomenology is thus determined by the single coefficient M 4 . This leads to interesting characteristic features, such as the possibility of small speed of sound c 2 s of scalar perturbations (if

|M 4 | |M 2
Pl Ḣ|) and, associated to it, large non-gaussianities f NL ∼ 1/c 2 s [START_REF] Chen | Observational signatures and non-Gaussianities of general single field inflation[END_REF][START_REF] Cheung | The Effective Field Theory of Inflation[END_REF]. The contribution of higher-order terms amounts to only slightly correcting the leading results. For example, one can show that the correction to the speed of sound from the operator M 2 (which arises entirely due to mixing with gravity) is of order δc 2 s ∼ M 2 /M 2 Pl ∼ H/M Pl and can be safely ignored. Consequently, one can consistently consider the perturbations of the inflaton field as weakly coupled over a sufficiently broad range of distances encompassing the Hubble scale and straightforwardly apply the derivative expansion. This would correspond to the case of k -inflation [START_REF] Armendariz-Picon | k -inflation[END_REF], ghost inflation [START_REF] Arkani-Hamed | Ghost inflation[END_REF] and other related models.

L WBG I M 4 Λ -4 2 M 3 HΛ -4 2 M 2 H 2 Λ -4 2 m2 H 2 Λ -4 2 I = 2 4X 2 G 2XX × × × I = 3 -6XZ (G 3X + 2XG 3XX ) -2XZG 3X × × I = 4 24XZ 2 (3G 4XX + 2XG 4XXX ) 8XZ 2 G 4X X + 2G 4XX -4Z 2 G 4X M 2 H 2 Λ -4 2 I = 5 2XZ 3 3 G 5X X + 12G 5XX + 4XG 5XXX 2XZ 3 3 G 5X X + 2G 5XX -2Z 3 G 5X M 2 H 2 Λ -4 2
Table 1: Contribution from the Lagrangian terms L WBG I to various unitary-gauge operators, defined in Eq. (5.2). The relation M 2 = m2 holds for all of these. The subscript X means differentiation with respect to X and all derivatives are evaluated on the background. We have assumed that all background quantities obey 'slow-roll', i.e. d/dt H, so that terms involving derivatives of X and Z have been neglected. Moreover, L WBG 4 generically contributes a weakly time-dependent piece to the effective Planck mass through the function f (t) in Eq. (5.2). According to our assumption about the form of G 4 discussed in Sec. 3, this is at least of order X 2 for small X, so that its contribution can be neglected. An analogous piece coming from L WBG 5 is suppressed by a further factor of the slow-roll parameter.

Inflation with WBG symmetry

Going back to the theories characterized by WBG invariance (4.1), one can straightforwardly find the values of the various coefficients in Eq. (5.2), by means of reformulating the discussion of the previous section into the effective theory language [START_REF] Gleyzes | Essential Building Blocks of Dark Energy[END_REF]. The contribution from each Lagrangian L WBG I to different EFTI operators is reported in Tab. 1.

Interestingly, as a generic property of Horndeski theories, only one combination of the operators δK 2 , δK ij δK ij and (3) RδN , given by the following expression

-δK 2 + δK ij δK ij + (3) R δN , (5.7)
appears in the EFT action [START_REF] Gleyzes | A unifying description of dark energy[END_REF]. Thus, the parameters M 2 and m2 are related to each other as M 2 = m2 , while all the others are zero, M 2 = m 1 = . . . = 0. Moreover, it follows from the discussion of Secs. 2 and 3 that this 'tuning' is perfectly stable under loop corrections, as far as X ∼ < 1 and Z ∼ < 1. This implies, for instance, that M 2 M 2 can hold in the full quantum theory due to the (weakly broken) galileon invariance. Moreover, the unique combination (5.7) is in fact a redundant operator and can be removed by a redefinition of the metric [START_REF] Gleyzes | A unifying description of dark energy[END_REF][START_REF] Gleyzes | Effective Theory of Interacting Dark Energy[END_REF] 13 . As a consequence, the single operator δN δK -associated with the cubic galileon -is responsible for all the differences with respect to the more familiar DBI/ghost inflation models, as far as stability and power spectra are concerned.

With this in mind, we derive the quadratic action for the comoving curvature perturbation ζ (defined in the standard way, γ ij = a 2 e 2ζ δ ij ) for the general case of nonzero coefficients M 4 and M 3 in Appendix B. The full result (away from any decoupling/short-distance limit) reads

S ζ = d 4 x a 3 A(t) ζ2 -B(t) ( ∇ζ) 2 a 2 , (5.8)
where

A = M 2 Pl (-4M 4 Pl Ḣ + 3 M 6 + 2M 2 Pl M 4 ) ( M 3 -2M 2 Pl H) 2 , ( 5.9)
B = M 2 Pl (-4M 4 Pl Ḣ + 2M 2 Pl H M 3 -M 6 + 2M 2 Pl ∂ t M 3 ) ( M 3 -2M 2 Pl H) 2 , ( 5.10)
and the speed of sound for short-wavelength perturbations is simply c 2 s = B/A. Depending on the precise values of the parameters X and Z, there are two phenomenologically distinct regimes in a theory of inflation with WBG invariance, which we review in what follows.

The kinetically driven phase

The strongly coupled backgrounds -i.e. the ones for which all terms (3.1)-(3.4) are of the same order -correspond to the function F (X 0 , Z 0 ) vanishing at the leading order, as already noted above. For X and Z strictly constant, one has an exact de Sitter space, which can be made quasi-de Sitter by turning on a small potential for φ. One can see from Eq. (4.4) that for F to vanish requires that at least Z 0 be order-one. Although less evident, it is easy to show that also X 0 has to be order-one, if the potential is to provide a sub-leading contribution to the energy density. Indeed, it follows from Eq. (4.6) that the energy density scales as M 2 Pl H 2 ∼ Λ 4 2 X in this case. Together with the following equality,

Z = √ X M Pl H Λ 2 2 , (5.11)
which simply follows from the very definitions of X and Z, this implies that Z ∼ X. One thus concludes that whenever the theory approaches strong coupling on a (quasi) de Sitter background, Z ∼ 1, all functions G I should be generically resummed. We stress again that this only concerns the case when V (φ) can be neglected at the leading order -a condition we will give up below.

For large backgrounds that nearly saturate the bound (3.12), one consequence of the approximate symmetry under (1.1) is the following radiatively stable hierarchy

M 2 Pl H 2 ∼ M 4 ∼ M 3 H ∼ M 2 H 2 ∼ m2 H 2 m 1 H 3 , M 2 H 2 , . . . , (5.12)
as can be seen from Tab. Pl respectively. For the latter values of these operators, their effects become order-one important for physical quantities such as, for instance, the power spectrum or the speed of sound of the scalar perturbations -see the exact expressions (5.9) and (5.10). In this case, derivative interactions either dominate or contribute significantly to the background curvature, unlike potentially-driven models such as for instance slow-roll or DBI inflation. Remarkably enough, for a sufficiently large hierarchy between Λ 2 and Λ 3 -or, equivalently, between M Pl and Λ 3 , the strong coupling of our backgrounds does not necessarily imply the breakdown of the classical description, as follows from the estimates of Eq. (4.8). For example, in terms of the EFT language, the operator ( φ) 2 /Λ 2 3 leads to the δK

2 term of order δ M 2 ∼ Λ 4 2 /Λ 2 3 ∼ M 2 Pl (H/M Pl ) 2/3
, which is suppressed by a tiny factor compared to the leading contribution 14 .

There is another important characteristic to the 'large' backgrounds X ∼ 1, Z ∼ 1 -or, equivalently, to large values of the EFT coefficients given in Eq. (5.12). Assuming M 2 Pl Ḣ ∼ < M 4 , as required by being close to de Sitter space, there is no short-distance limit in which the metric can be decoupled from the inflaton perturbations. Therefore, mixing with gravity is order-one important at all scales. To see this, let us restore the Goldstone boson π, nonlinearly realizing the spontaneously broken time diffs:

g 00 → g 00 = -N -2 = -1 -2 π -π2 + (∂π) 2 a 2 , δK → δK - ∂ 2 π a 2 + O(π 2 ) .
The relevant part of the π action then reads

-M 2 Pl Ḣ 1 N 2 + 1 2 M 4 (δN ) 2 -M 3 δN δK → π2 c + M 3 M Pl f 2 π δN c ∂ 2 π c a 2 + . . . , (5.13)
where we have defined the 'decay constant' f 2 π ≡ (M 4 /2 -M 2 Pl Ḣ) 1/2 and canonically normalized the Goldstone and the lapse variable as π = π c /f 2 π and δN = δN c /M Pl , respectively. One can see from the last term in this expression, that the mixing with gravity is indeed important for the Goldstone dynamics at all distance scales if the EFT coefficients saturate the strong coupling bound (i.e. if M 3 ∼ M Pl M 2 ). This means that, in this case, one has to perform the full analysis including dynamical gravity in order to extract reliable predictions from the theory. That gravity is order-one important for such backgrounds can also be seen from the exact expressions for the kinetic coefficients of the curvature perturbation (5.9) and (5.10). The calculation of the cubic action/non-gaussianities in this regime is much more involved due to the non-decoupling of gravitational effects and will be reported on elsewhere.

In principle, Z may end up being somewhat smaller than one due to an accident/fine tuning (as discussed above, in our case any tuning of coefficients is respected by quantum corrections in the theories at hand). In this case, one can then define the decoupling limit governed solely by the dynamics of the Goldstone mode or, equivalently, by the comoving curvature perturbation ζ = -Hπ. In this case the analysis beyond the quadratic action becomes significantly simpler. The non-gaussianity in the decoupling-limit has been studied in the context of galileon inflation in Ref. [START_REF] Burrage | Galileon inflation[END_REF] and we will not reproduce the results here. Unlike DBI/ghost inflation, this limit contains cubic operators with four derivatives, such as for instance ∂ 2 π(∂π) 2 . However, all of these operators are redundant and can be reduced on-shell to those with three derivatives [START_REF] Creminelli | Galilean symmetry in the effective theory of inflation: new shapes of non-Gaussianity[END_REF], i.e. π(∂π) 2 and π3 . This means that, at least in the decoupling limit, galileon inflation does not lead to new shapes of non-gaussianity, distinct from those that already arise in DBI models. The genuine novelty is that, unlike DBI theories, more than one operator becomes nonlinear in galileon inflation, leading to a richer parameter space. As a result, the amplitude of equilateral non-gaussianity f equil NL is not necessarily linked to the speed of sound of scalar perturbations and it can be larger than the standard DBI result, f equil NL ∼ 1/c 2 s , see [START_REF] Burrage | Galileon inflation[END_REF] for a more detailed discussion.

Slow-roll inflation with WBG symmetry

A different and particularly interesting case corresponds to a potentially driven universe with a subleading contribution to the energy density due to the derivative part of the action,

M 2 Pl H 2 ∼ V Λ 4 2 X .
(5.14)

In the inflationary context this is just the slow-roll limit.

The slow-roll parameter can be estimated from the Friedmann equations (4.6) and (4.7) to be of order

≡ - Ḣ H 2 ∼ Λ 4 2 X M 2 Pl H 2 , (5.15)
where we have made use of the fact that F is an order-one function and have assumed that F and φ vary slowly with time, Ḟ HF, φ H φ. Using the identity (5.11), the last expression for immediately implies that Z has to be parametrically larger than X, Z ∼ X √ .

(5.16)

Let us consider now a potentially-driven regime, characterized by Z ∼ 1. According to Eq. (5.16), X is parametrically suppressed, being of order √ , and it follows from (5.15), that

M 2 Pl Ḣ = M 2 Pl H 2 ∼ Λ 4 2 X .
(5.17)

The latter scale, denoted by f π in Eq. (5.13), defines the scalar two-point function in ordinary slow-roll models, setting the normalization for the Goldstone mode of broken time diffs,

L = -M 2 Pl Ḣ(∂ µ π) 2 + . . . , (5.18)
where ellipses denote all sorts of additional terms, completely unimportant in slow-roll inflation, including those describing the mixing of π with gravity.

Let us now have a look at the values of the additional operators present in the case of interest and listed in Tab. 1. A crucial observation is that, when Z is of order one, the contribution to M 4 and M 3 from all L WBG

I with I = 2 scales as M 4 ∼ Λ 4 2 X ∼ M 2 Pl Ḣ and M 3 ∼ M 2
Pl Ḣ/H respectively. The M 4 operator that arises from L WBG 2 is down by an additional factor of X ∼ √ compared to this, so that, as long as perturbations are concerned, higher-derivative terms are more important than the next-to-leading order terms in the derivative expansion of the form (∂φ) 2n . Such a scaling of extra EFT operators is quite remarkable, since it implies that the perturbation Lagrangian (5.18) is modified at order one compared to the usual slow-roll case, while the background evolution is still fully governed by the potential 15 .

This can also be seen from the exact expressions for the quadratic scalar perturbation Lagrangian. Taking into account the magnitudes of the two EFTI operators at hand, let us denote

M 4 = -αM 2 Pl Ḣ = αM 2 Pl H 2 , M 3 = βM 2 Pl Ḣ H = -βM 2 Pl H , (5.19)
where α and β are constants of order one. Then, from Eqs. (5.9) and (5.10) one finds

A = 1 + α 2 M 2 Pl + O( 2 ) , B = 1 - β 2 M 2 Pl + O( 2 ) , c 2 s = 2 -β 2 + α . (5.20)
This explicitly shows that, by suitably choosing the order-one coefficients α and β, which themselves are functions of the free parameters of the theory, the speed of sound of the scalar perturbations end up anywhere in the allowed region 0 ≤ c 2 s ≤ 1. Moreover, tuning β to be close to 2 yields a regime with highly subluminal scalar perturbations, potentially implying strong coupling and large 'accidental' non-gaussianities [34], all in a technically natural way (recall that any tuning in the theories at hand is protected by WBG invariance).

Last but not least, we note that in the potentially driven regime, one can generically define the decoupling limit -unlike the kinetically dominated case considered above. This can be seen by noting that the term describing the mixing of the (canonically normalized) Goldstone π c with gravity is parametrically suppressed on backgrounds close to de Sitter. Indeed, using the expressions for the magnitudes of EFT coefficients in Eq. (5.19), one can readily evaluate the coefficient of the mixing term in Eq. (5.13) to be of order

M 3 M Pl f 2 π ∼ √ . (5.21)
During inflation, this is a parametrically small number meaning that the Goldstone mode can be studied independently from the gravitational degrees of freedom, at frequencies of order of the Hubble scale. The same conclusion can be made regarding all other EFT operators present in the theory.

One can also see from Table 1 that, while generalizing galileon inflation, the most general inflationary theory defined by (4.1) possesses a much broader parameter space. While three parameters in Tab. 1 characterize galileon inflation (G 3X , G 4XX , G 5XX ), a generic theory with WBG invariance has six more parameters 16 at the level of the two-point scalar function (G 4X , G 5X , G 2XX , G 3XX , G 4XXX , G 5XXX ). Higher-point functions will introduce even more freedom. Moreover, arbitrary tuning of the parameters is stable under radiative corrections and taking advantage of this fact can lead to interesting theories in various corners of the EFT of inflation theory space. The complete analysis of possible physical predictions of the theory (4.1) is beyond the scope of the present discussion. 6 Late-time universe

Dark energy

The theories characterized by WBG symmetry can also be applied in the context of the late-time cosmic acceleration. For simplicity, we will neglect the scalar potential V (φ) in this section, contrarily to the case discussed in the previous section. Moreover, we reabsorb the canonical scalar kinetic term in the definition of the function G 2 , so that the action of interest reads as follows

S = d 4 x √ -g 1 2 M 2 Pl R + 5 I=2 L WBG I . (6.1)
The homogeneous equations on a flat FRW background obtained by varying this action are similar to those, discussed in Sec. 5. However, in this section we choose to rewrite them in a slightly different form and notation.

In terms of the two variables X and Z defined in Eqs. (3.5) and (4.5), the homogeneous scalar evolution equation reads [START_REF] Kobayashi | G-inflation: Inflation driven by the Galileon field[END_REF]

d dt a 3 φF DE (X, Z) = 0 , (6.2)
where F DE (X, Y ) is defined as

F DE (X, Z) ≡ G 2X -3ZG 3X + 6Z 2 G 4X X + 2G 4XX + Z 3 3 G 5X X + 2G 5XX . (6.3)
We assume that all matter couples minimally to the metric g µν , so that the Friedmann equations can be written as

3M 2 * H 2 = ρ m + ρ DE , (6.4)
-2M 2 * Ḣ = ρ m + p m + ρ DE + p DE , (6.5)
where ρ m and p m respectively denote the energy density and pressure of matter. The effective Planck mass squared, M 2 * , is generically time-dependent due to the non-minimal couplings of the scalar in L WBG 4 and L WBG 5 , and is explicitly given by the following expression

M 2 * ≡ M 2 Pl (1 + 2G 4 -4XG 4X -2ZXG 5X ) . (6.6)
In terms of the EFT parameters of Eq. (5.2),

M 2 * = M 2 Pl f + M 2 .
Finally, the unperturbed effective energy density and pressure read

ρ DE = Λ 4 2 X - G 2 X + 2G 2X -6ZG 3X +12Z 2 G 4X X + 2G 4XX + 4Z 3 G 5X X + G 5XX , (6.7)
p DE + ρ DE = 2Λ 4 2 XF DE + 2M Pl φX × G 3X -4Z G 4X X + 2G 4XX -Z 2 3 G 5X X + 2G 5XX . (6.8)
For Horndeski theories, the phenomenological deviations from the ΛCDM model at the level of linear perturbations can be conveniently parametrized in terms of four timedependent parameters [START_REF] Gleyzes | Essential Building Blocks of Dark Energy[END_REF]. It is convenient to use the dimensionless functions introduced in [START_REF] Bellini | Maximal freedom at minimum cost: linear large-scale structure in general modifications of gravity[END_REF], which can be expressed in terms of the EFT coefficients of Eq. (5.2) as [START_REF] Gleyzes | A unifying description of dark energy[END_REF]

α K = 2c + M 4 M 2 * H 2 , α B = M 2 Pl ḟ -M 3 2M 2 * H , α M = M 2 Pl ḟ + Ṁ 2 M 2 * H , α T = - M 2 M 2 * , (6.9)
where the function c is fixed by the background Friedmann equations derived from action

(5.2), 2c = p DE + ρ DE + M 2 Pl (H ḟ -f ) . (6.10)
For the theory of interest (6.1), their explicit expressions are reported in App. C. The contributions from the Lagrangians L WBG I to each of these parameters is given in Tab. 2.

Let us first briefly comment on the physical meaning of the α-parameters and their observational consequences (see, e.g. [START_REF] Bellini | Maximal freedom at minimum cost: linear large-scale structure in general modifications of gravity[END_REF][START_REF] Gleyzes | A unifying description of dark energy[END_REF] for a more detailed discussion). • The first of these functions, α K , parametrizes the kinetic energy of scalar perturbations, induced by the four Lagrangian terms L WBG

L WBG I M 2 * H 2 Λ -4 2 α K M 2 * H 2 Λ -4 2 α B M 2 * H 2 Λ -4 2 α M M 2 * H 2 Λ -4 2 α T I = 2 2X(G 2X + 2XG 2XX ) × × × I = 3 -12XZ(G 3X + XG 3XX ) XZG 3X × × I = 4 12XZ 2 G 4X X + 8G 4XX + 4XG 4XXX -4XZ 2 G 4X X + 2G 4XX - φ Λ 3 3 α B 4Z 2 G 4X I = 5 4XZ 3 3 G 5X X + 7G 5XX + 2XG 5XXX -XZ 3 3 G 5X X + 2G 5XX 2 ḢZ 3 H 2 G 5X -2 φ Λ 3 3 α B 2Z 2 Z - φ Λ 3 3 G 5X

I

. In terms of the EFT for perturbations (5.2), the contributions to α K arise from the operators ∝ δN 2 , described by the coefficients c and M 4 . This parameter is enough to describe linear perturbations in the minimally coupled quintessence and k-essence models. In the minimal case, dark energy fluctuations behave as those of a perfect fluid with the speed of propagation

c 2 s = 3(1 + w DE )Ω DE /α K ,
where Ω DE is the energy density of φ in the units of the critical one. A well-studied example corresponds to the limit α K 1, which leads to zero sound speed [START_REF] Creminelli | The Effective Theory of Quintessence: the w¡-1 Side Unveiled[END_REF][START_REF] Creminelli | Spherical collapse in quintessence models with zero speed of sound[END_REF].

• The second function, α B , parametrizes the mixing between metric and scalar field fluctuations [START_REF] Deffayet | Imperfect Dark Energy from Kinetic Gravity Braiding[END_REF]. It can be induced by the Lagrangians L WBG I with 3 ≤ I ≤ 5, or by the operator containing δN δK in the EFT for perturbations, which describes kinetic mixing with gravity. This operator induces a typical scale (approximately given by k B /a 3H Ω m /2 for α M = α T = 0) [START_REF] Bellini | Maximal freedom at minimum cost: linear large-scale structure in general modifications of gravity[END_REF], below which dark energy can cluster with energy density fluctuations δρ DE 2δρ m α 2 B / ((α 2 K + 6α 2 B )c 2 s ). Above this scale, dark energy fluctuations simply behave as those of a perfect fluid. . It parametrizes the non-minimal coupling in scalar-tensor theories such as Brans-Dicke [START_REF] Brans | Mach's principle and a relativistic theory of gravitation[END_REF] and, as such, induces a slip between the gravitational potentials. In terms of the EFT operators, α M can be generated by time variation of the parameters f and M 2 . In particular, in f (R) theories, α M = -α B .

• Finally, the last function, α T , parametrizes the deviation of the speed of sound of tensor perturbations from unity in the presence of the Lagrangian L WBG 4 and L WBG 5 [START_REF] Kobayashi | G-inflation: Inflation driven by the Galileon field[END_REF], and also induces a slip between the gravitational potentials. It is generated by the EFT operator, multiplied by the coefficient M 2 .

A particularly interesting case corresponds to a constant X, where the scalar field grows linearly with time, φ 0 ∝ t. In this case, the functions G I , and all their derivatives are time-independent when evaluated on the background, while Z changes proportionally to the Hubble rate. One can see from the form of the scalar equation, Eq. (6.2), that the φ 0 ∝ t profile can only be a solution if the following relations are satisfied

G 2X = 0 , G 3X = 0 , G 4X + 2XG 4XX = 0 , 3G 5X + 2XG 5XX = 0 . (6.11)
It then follows from Eqs. (6.7) and (6.8) that the effective energy density on such backgrounds is given by

ρ DE = -Λ 4 2 (G 2 + 2Z 3 G 5X
) and the equation of state is that of the cosmological constant, with w DE ≡ p DE /ρ DE = -1. Moreover, while the dimensionless parameters α K , α M and α T in Tab. 2 are generally nonzero and time dependent, Eq. (6.11) implies that α B = 0. In this case, in order to avoid ghost instabilities one must require that α K ≥ 0 [START_REF] Gleyzes | Essential Building Blocks of Dark Energy[END_REF][START_REF] Bellini | Maximal freedom at minimum cost: linear large-scale structure in general modifications of gravity[END_REF]. Note that since Z ∝ H, ρ DE is time dependent, even though w DE = -1. This is because the effective Planck mass varies with time and the dark energy does not follow the usual conservation equation. Demanding that the energy density remains finite at early times requires it to be constant, i.e. G 5X = 0 and thus that also α M = 0, see Tab. 2. In this case, at the background level this model exaclty behaves as a cosmological constant. At the level of linear perturbations, the dimensionless parameters α K and α T do not vanish and contribute to the sound speed of scalar and tensor fluctuations, respectively as c 2 s = -2α T /α K and c 2 T = 1 + α T (see [START_REF] Bellini | Maximal freedom at minimum cost: linear large-scale structure in general modifications of gravity[END_REF][START_REF] Gleyzes | A unifying description of dark energy[END_REF]). One can avoid gradient instabilities by requiring the positivity of both these speed of fluctuations squared, which implies -1 ≤ α T ≤ 0. This is the simplest application of WBG symmetry to the late time-acceleration. The assumption that the background profile of the field is linear with time, φ 0 ∝ t, and that the background expansion history is the same as in ΛCDM considerably restricts the values that the parameters α a can consistently take. These solutions are not possible for the covariant galileon [START_REF] Deffayet | Covariant Galileon[END_REF]. In this case one usually assumes a tracker solution with Z = const. [START_REF] De Felice | Cosmology of a covariant Galileon field[END_REF]; this imposes a particular expansion history, which has been shown to lead to observations that are disfavoured with respect to ΛCDM (see [START_REF] Barreira | The observational status of Galileon gravity after Planck[END_REF] for a recent analysis and e.g. [START_REF] Appleby | The Paths of Gravity in Galileon Cosmology[END_REF][START_REF] Neveu | Experimental constraints on the uncoupled Galileon model from SNLS3 data and other cosmological probes[END_REF] for analysis that do not assume the tracker solution). More sophisticated examples can be constructed using the Lagrangians (3.1)- (3.4), by allowing the background solution for φ to be different from a profile linear in time. In this case, from Tab. 2 one generically expects that the parameters α a may assume any value smaller than unity, α a ∼ X ∼ < 1.

Static sources and the Vainshtein mechanism

Any theory with gravitationally coupled scalars should be able to pass the solar system constraints in order to be able to qualify as an acceptable alternative to general relativity. In the context of theories with WBG invariance, one can imagine coupling φ to the trace of the matter stress tensor in the following way

L mat ∼ 1 M Pl φT m . (6.12)
It is then necessary to make sure that the exchange of φ by realistic sources does not modify the Newtonian potential by any detectable amount. There are a number of ways to suppress the contribution from φ to the gravitational potential of astrophysical objects (stars, planets, etc.) in modified gravity. In the case of the theories with approximate galileon invariance we are interested in, what guarantees that φ is screened beyond the observable values is the Vainshtein mechanism [START_REF] Vainshtein | To the problem of nonvanishing gravitation mass[END_REF][START_REF] Deffayet | Nonperturbative continuity in graviton mass versus perturbative discontinuity[END_REF]. Indeed, in the decoupling limit (3.8), the theory reduces to the Galileon, which is known to naturally incorporate the screening below the Vainshtein redius [START_REF] Nicolis | The Galileon as a local modification of gravity[END_REF]. For a souce of mass M source , this radius reads

r V = M source M Pl Λ 3 3 1/3
. (6.13)

The decoupling limit, on the other hand, does capture all relevant astrophysical scales and all results obtained in that limit can be fully trusted. To see this directly, one can evaluate the quantity X on a pure Galileon solution of Ref. [START_REF] Nicolis | The Galileon as a local modification of gravity[END_REF] and check that it is extremely small everywhere in space 17 , so that all the results of Ref. [START_REF] Nicolis | The Galileon as a local modification of gravity[END_REF] apply to our case without any modification.

Conclusions and outlook

Galileon-symmetric scalar theories have played an important role in cosmological modelbuilding in the context of both the early and the late-time universe. The exact symmetry is only possible in the absence of gravity and it leads to remarkable properties -both at the classical and quantum levels. However, upon inclusion of gravity galileon invariance necessarily has to be broken, raising the question regarding the fate of its attractive flatspace properties.

In this paper, we have introduced the notion of weakly broken galileon invariance, which allows to naturally generalize the galileon to curved space (see Sec. 2). Requiring that the quantum non-renormalization properties of the flat-space theory is preserved to the maximal possible extent also in the presence of gravity has led us, quite uniquely, to a particular subclass of the ghost-free Horndeski theories, presented in Sec. 3. The resulting effective theory is characterized by two scales -or, more appropriately, by parametrically small coefficients of symmetry-breaking operators in the units of the effective field theory cutoff -allowing to retain the quantum non-renormalization properties of the galileon for a broad range of physical backgrounds.

As we have shown, the latter properties guarantee that de Sitter solutions that these theories generically possess are insensitive to loop corrections. This gives a lot of freedom to implement these solutions for building technically natural models of both the inflationary, as well as the late-time universe. One of the purposes of this work has been to introduce these models and set the stage for their detailed phenomenological studies.

In Sec. 5, we have studied two relevant regimes in the context of inflation. The first one describes genuinely strongly coupled physics, characterized by X = φ2 /Λ 4 2 ∼ 1 and Z = H φ/Λ 3 3 ∼ 1 in the notation of Sec. 3. It generalizes galileon inflation of Ref. [START_REF] Burrage | Galileon inflation[END_REF] and saturates the bound on the 'strength' of the backgrounds for which the classical analysis can still be trusted. The dynamics of perturbations on these backgrounds is completely dominated by higher-derivative operators. Yet, the mechanism is fully predictive, quantum loops giving subleading corrections to the classical results. In a generic theory with weakly broken galileon invariance, there are six additional parameters (two of which have been set to zero in this paper) compared to the particular case of galileon inflation at the level of the scalar two-point function. More freedom arises when higher n-point functions are concerned.

The other, particularly interesting, regime corresponds to partially strongly coupled backgrounds, with Z ∼ 1 and X ∼ √ , resembling slow-roll inflation. In this scenario, inflation happens mostly due to the potential; however, weakly broken galileon invariance makes it possible for the higher-derivative operators to become relevant for the dynamics of scalar perturbations, modifying their action at order one. In particular, this allows for a significantly subluminal speed of sound for the curvature perturbation ζ as a result of a possible parameter tuning -which is, importantly, respected by loop corrections to the leading order in M -1 Pl . The small speed of sound generically translates into strong coupling of scalar perturbations, leading to large non-gaussianities. Such 'accidentally' moderate/strong non-gaussian perturbations can never arise in ordinary slow roll models. Thus, this represents an interesting possibility that the theories with weakly broken galileon invariance can offer. These matters will be studied in detail in an upcoming work [34].

The theories introduced in Sec. 3 can also be used to drive the current accelerated expansion of the universe, as discussed in Sec. 6. The situation is richer than the covariant galileon case studied in the literature. In particular, while one can recover a background expansion history closed to the one currently observed, at the level of linear perturbations these theories can lead to modifications of gravity of different types. Expressed in terms of the dimensionless effective field theory parameters α a defined for Horndeski theories, one generically finds that these cannot be too large, α a ∼ < 1, leaving the possibility of modifications which are compatible with current constraints but possibly testable by future observations. We leave a more thorough study of concrete scenarios for the future.

As a final remark, we note that we have not carried out a complete analysis of quantum corrections in the theories with weakly broken galileon invariance; we have rather made the most conservative estimates of their magnitudes, so that the statements made above should remain true for the most general completion of these theories above the scale Λ 3 . Whether or not our findings regarding the structure of the low-energy effective theory (4.1) imply anything useful about its UV completion, is still to be seen. We postpone the detailed study of this question for future work.

implying that the two contributions to the possible vertex with one graviton and three solid lines exactly cancel for our resulting generalized theory.

The way we have chosen to couple the galileon to gravity has been dictated by our desire to raise the scale suppressing loop-generated symmetry-breaking operators. However, there is something more to it: the resulting non-minimal theory can be recognized as that of Horndeski type, leading to second-order equations both for the metric and the scalar. In the retrospect, this is not that surprising: the 'three solid line' vertices, if present, are the only ones that would introduce higher-order equations for the metric in the theory of minimally coupled galileon, as can be easily seen from e.g. Eq. (A.2). Our non-minimal term exactly cancels this contribution. We thus observe something similar to what happens in massive gravity [START_REF] De Rham | Generalization of the Fierz-Pauli Action[END_REF]: there, eliminating terms that lead to higher-order equations (or, equivalently, to the Boulware-Deser ghost) automatically raises the quantum cutoff of the theory; in our case, the cutoff (Λ 3 ) remains unaltered, but a similar procedure raises the scale of breaking for the galileon symmetry.

The same result can be straightforwardly extended to the case of the quintic galileon. A brute-force expansion of the minimally-coupled term, using (A.1), yields

2 3 L min 5 ∼ (∂φ) 2 ∂ ρ φ ∂ ρ h ( φ) 2 -(∂ µ ∂ ν φ) 2 -2∂ ρ h µν (∂ µ ∂ ν φ φ -∂ µ ∂ α φ∂ ν ∂ α φ) -2∂ µ h ρµ ( φ) 2 -(∂ α ∂ β φ) 2 + 4∂ µ h ρν (∂ µ ∂ ν φ φ -∂ µ ∂ α φ∂ ν ∂ α φ) . (A.4)
It is impossible to manipulate the last expression by partial integration or otherwise, so as to be left with less than three factors of ∂φ. This would lead to the last vertex with three solid lines in Fig. 2, and therefore higher-order equations for the metric. One can however show that this term can be eliminated by adding a non-minimal piece with an exact coefficient corresponding to the Horndeski theory. To this end, we note that to the linear order in the metric perturbation, the following relation holds

√ -g(∂φ) 4 G µν ∇ µ ∇ ν φ = - 1 2 (∂φ) 4 µαρλ νβσλ ∂ α ∂ β h ρσ ∂ µ ∂ ν φ + O(h 2 ) , (A.5)
where denotes the totally antisymmetric symbol with 1234 = 1. Expanding the antisymmetric product and integrating by parts, we obtain

(∂φ) 4 G µν ∇ µ ∇ ν φ ∼ (∂φ) 2 ∂ ρ φ ∂ ρ h (∂ µ ∂ ν φ) 2 -( φ) 2 + 2∂ ρ h µν (∂ µ ∂ ν φ φ -∂ µ ∂ α φ∂ ν ∂ α φ) .
This exactly cancels the first line of the minimal term's contribution (A.4). Partially integrating the second line of (A.4) on the other hand, gives

-2(∂φ) 2 ∂ ρ φ∂ ν h ρµ ∂ µ ∂ α φ∂ ν ∂ α φ -∂ µ ∂ ν φ φ + ∂ µ φ∂ ν φ -∂ µ ∂ ν ∂ α φ∂ α φ . (A.6)
One can now check, that whatever stands to the right of (∂φ) 2 ∂ ρ φ in the last expression, is a total derivative (∂ ν of a local operator). This means that one final partial integration can get rid of one solid line in the corresponding vertex, reducing the number of solid lines to two. This proves our statement regarding elimination of vertices with one graviton and more than two solid lines.

One last statement we wish to prove regards the absence of vertices with two gravitons and more than four solid lines. The only possible obstruction comes from the non-minimal coupling obtained above, (∂φ) 4 G µν ∇ µ ∇ ν φ, in which we take the Einstein tensor linear in the metric perturbation, and pick up another factor of h from the covariant derivative 19 . One can show however, that the corresponding term can always be put into the form in which there are no more than four scalars with a single derivative acting on them. Indeed, using the expression for G µν in terms of the antisymmetric symbol, we have

(∂φ) 4 G µν ∇ µ ∇ ν φ ∼ 1 4 (∂φ) 4 µαρλ νβσλ ∂ α ∂ β h ρσ (2∂ µ h γν -∂ γ h µν )∂ γ φ . (A.7)
This form makes it straightforward to convince oneself, that the factors involving two gravitons in the above expression collect into either a total derivative, or a total derivative up to corrections involving at most four factors of ∂φ:

µαρλ νβσλ ∂ α ∂ β h ρσ ∂ µ h γν = ∂ µ ( µαρλ νβσλ ∂ α ∂ β h ρσ h γν ) (∂φ) 4 ∂ γ φ µαρλ νβσλ ∂ α ∂ β h ρσ ∂ γ h µν ∼ 1 2 (∂φ) 4 ∂ γ φ ∂ γ ( µαρλ νβσλ ∂ α ∂ β h ρσ h µν ) . (A.8)
At this point, all that remains is to do a single partial integration to complete the proof of the statement, made in the beginning of the present section.

B Spectrum of (near-) de Sitter space Exact de Sitter space saturates the Null Energy Condition (NEC), which seems to be one of the most robust ways of covariantly expressing the requirement of energy positivity in general relativity 20 . This seems to be supported by the fact, that at least in the simplest effective theories -those that obey the standard derivative expansion with non-degenerate leading-order effects -NEC-violation inevitably implies either ghost or gradient instability at short distances [START_REF] Dubovsky | Null energy condition and superluminal propagation[END_REF]. Moreover, the exact dS limit of the corresponding backgrounds can only be stable over a finite time, necessarily suffering from Jeans-like instabilities [START_REF] Creminelli | Starting the Universe: Stable Violation of the Null Energy Condition and Non-standard Cosmologies[END_REF]. Furthermore, there is additional, purely quantum evidence against the existence of exact de Sitter spacetime [START_REF] Dvali | Quantum Exclusion of Positive Cosmological Constant?[END_REF]. In view of these facts, it is important to understand the most general conditions for stability of dS as well as its small deformations.

In this appendix, we study this question for the effective theory governed by WBG invariance. Since these theories lead to second-order equations of motion, the dispersion relation for the scalar perturbations will have a ω 2 ∝ k 2 form. Moreover, as noted above, the operators M 2 and m2 are redundant, and can be removed by a disformal transformation on the tensor mode; one can therefore concentrate on the only remaining quadratic operators: δN 2 and δN δK in the effective theory (5.2). It is convenient to fix the unbroken spatial diffs so as to put the 3D metric in the following form [START_REF] Maldacena | Non-Gaussian features of primordial fluctuations in single field inflationary models[END_REF]

g ij = a 2 (t)[(1 + 2ζ)δ ij + γ ij ], ∂ i γ ij = γ ii = 0 , (B.1)
where ζ and γ capture the physical scalar and tensor perturbations. Modes of different helicity do not mix on the homogeneous and isotropic backgrounds considered here, so that we can discard the helicity-1 part of the shift altogether, setting N i ≡ δ ij ∂ j β. In the effective theory for perturbations (5.2), both the lapse and the shift can be integrated out from their respective equations of motion21 :

M 2 Pl (3) R + K 2 -K ij K ij + 2 N 2 Ḣ -2(3H 2 + Ḣ) + 2M 4 δN -2 M 3 (δK -3HδN ) = 0 , (B.2) ∇ i M 2 Pl K i j -δ i j K -δ i j M 3 δN -M 2 δK i j -M 2 δKδ i j = 0 , (B.3)
where the lapse perturbation is denoted by δN ≡ N -1 and all EFT coefficients except for M 4 and M 3 are set to zero for the reasons outlined above. Solving these equations for δN and ψ to linear order and substituting the solutions back into (5.2), one finds the following quadratic action for the curvature perturbation The speed of sound of short-wavelength scalar perturbations is given by c 2 s = B/A. The last two expressions agree with those of Ref. [START_REF] Creminelli | Galilean Genesis: An Alternative to inflation[END_REF] up to a change of operator basis.

S ζ = d 4 x

4 in

 4 the cosmological constant, which for simplicity we tune to zero, c normalization of the scalar and graviton kinetic terms respectively, while c

•

  The function α M is defined as α M ≡ d ln M 2 * / ln a, and is thus related to the time variation of the effective Planck mass M 2 * induced by the Lagrangians L WBG 4 and L WBG 5

  1. Here, ellipses denote all other coefficients of operators of various canonical dimension, properly normalized with the help of the appropriate number of powers of the Hubble scale. The strongly-coupled case at hand corresponds to the largest possible values of the coefficients M 4 , M 3 and { M 2 , m2 }, i.e. of order M 2 Pl H 2 , M 2 Pl H and M 2

Table 2 :

 2 Contribution from the Lagrangian terms L WBG I to the various effective α-parameters defined in Eq. (6.9), denoting the deviations from the ΛCDM model in linear perturbation theory. Their complete expressions are reported in App. C, see Eqs. (C.1)-(C.4).

In what follows, we will at times abuse the nomenclature by referring to this redundancy as 'time translations'.

In massive gravity, where φ describes the helicity-0 polarization of the graviton, this simply follows from the equivalence principle.

Unless gravity is massive, see[START_REF] Gabadadze | A Covariant Master Theory for Novel Galilean Invariant Models and Massive Gravity[END_REF].

Moreover, the size of these operators is an avatar for technical naturalness of non-trivial de Sitter vacua of the classical theory, as we will show in the next section.

We always assume in this paper that Λ 3 M Pl .

Moving the derivative on h to the rest of the fields in the vertex generically does not reduce the number of solid lines for the quartic and quintic galileons, since this extra derivative can act on the factor ∂ 2 φ, representing the dashed line.

This can be seen by inserting enough number of vertices with one graviton and three solid lines into a generic 1PI loop diagram.

Note that both G 4X /X and G 5X /X in this equation are finite even in the X → 0 limit due to our assumption that G 4 and G 5 start at least quadratic in X, see the discussion below Eq. (3.9).

This redefinition corresponds to transforming into the frame where the tensor modes propagate at the unit speed of sound[START_REF] Creminelli | Resilience of the standard predictions for primordial tensor modes[END_REF].

Note that quantum corrections generically introduce higher time derivatives (e.g. δ Ṅ 2 ) at the same scales as the higher space derivatives. These lead to loss of unitarity associated with ghost excitations with masses above the EFT cutoff (cured by whatever completes the theory in the UV). The scales associated with the higher space-and time-derivative operators can only be disentangled when the speed of sound is small, c s 1, in which case the latter scale becomes enhanced by powers of c -1 s .

This is analogous to what occurs in DBI inflation[START_REF] Silverstein | Scalar speed limits and cosmology: Acceleration from D-cceleration[END_REF], even though there are crucial differences. Here we have in mind a scenario closer to slow-roll inflation than DBI models, in that e.g. the Hubble friction is important for the inflationary dynamics of φ and sustaining inflation requires a flat enough potential -or, in other words, that the usual relation ∼ (M Pl V /V ) 2 holds. This is not necessarily true in DBI inflation.

Assuming G 4X and G 5X are zero, which we did in this paper.

In the case of the cubic Galileon, X becomes of order one only at the Schwarzschild radius of the source r S , while, in the most general case, it is a constant within the Vainshtein radius, of order of the tiny ratio (r S H 0 ) 2/3 , where H 0 is the current Hubble rate.

Although parametrically higher than Λ 3 , this symmetry-breaking scale is still low enough to disrupt the de Sitter spectrum (5.12).

A vertex with more than four solid lines coming from the minimally coupled quintic galileon would involve at least three gravitons.

This is unlike e.g. the Strong Energy Condition (SEC), which is known to be violated both in the early and in the late universe.

For the purposes of deriving the quadratic action for ζ only the linearized version of these equations matters. This is because the contribution of higher-order terms in δN and β to the free action for the curvature perturbation vanishes on-shell.

Here the parameter α B is defined as -1/2 the one defined in Ref.[START_REF] Bellini | Maximal freedom at minimum cost: linear large-scale structure in general modifications of gravity[END_REF].

Acknowledgements

It is a pleasure to thank Paolo Creminelli and Greg Gabadadze for valuable discussions. The work of D.P. is supported in part by MIUR-FIRB grant RBFR12H1MW and by funds provided by Scuola Normale Superiore through the program "Progetti di Ricerca per Giovani Ricercatori". The work of E.T. is supported in part by MIUR-FIRB grant RBFR12H1MW. F.V. acknowledges kind hospitality by Scuola Normale Superiore, where part of this project was completed.

A Quantum corrections

In this appendix, we extend the discussion of Sec. 2 to show that the vertices with three solid lines and one graviton, as well as five solid lines and two gravitons (see Fig. 2) can be removed in a suitable curved-space extension of the galileon terms.

We have considered the case of the cubic galileon in the main text. One has to work a little more to understand the case of quartic minimally coupled galileon. The contribution to the second vertex from the lower line of Fig. 2 can be extracted by picking up a factor of ∂h from the covariant derivative acting on one of the scalars. The relevant term reads,

Here (and in the rest of the present section), by '∼' we mean 'equals up to a total derivative and up to terms with fewer factors of ∂φ', and we do not distinguish between upper and lower indices for simplicity (everything is contracted with the flat-space metric). Using (A.1), we find after a little bit of algebra

where L min 4 denotes the Lagrangian term, obtained by minimally covariantizing (2.3). Inserting the latter vertex into a generic 1PI loop diagram would induce galileon symmetrybreaking operators at the scale M 1/6 Pl Λ 5/6 3 , parametrically exceeding Λ 2 -as discussed above 18 . Fortunately, it turns out possible to raise it by adding non-minimal couplings, capable of eliminating the vertices with three solid lines coming from the minimally coupled theory. For the quartic galileon, the right coupling is √ -g(∂φ) 4 R . Indeed, expanding the Ricci scalar to the linear order in h and integrating by parts, we obtain √ -g(∂φ

C Effective parameters for dark energy

Deviations from ΛCDM for WBG invariant theories can be parametrized in terms of the following time-dependent functions [START_REF] Bellini | Maximal freedom at minimum cost: linear large-scale structure in general modifications of gravity[END_REF] 22