Massive Galileon Positivity Bounds

The EFT coefficients in any gapped, scalar, Lorentz invariant field theory must satisfy positivity requirements if there is to exist a local, analytic Wilsonian UV completion. We apply these bounds to the tree level scattering amplitudes for a massive Galileon. The addition of a mass term, which does not spoil the non-renormalization theorem of the Galileon and preserves the Galileon symmetry at loop level, is necessary to satisfy the lowest order positivity bound. We further show that a careful choice of successively higher derivative corrections are necessary to satisfy the higher order positivity bounds. There is then no obstruction to a local UV completion from considerations of tree level 2-to-2 scattering alone. To demonstrate this we give an explicit example of such a UV completion.


Introduction
Low energy effective field theories of scalar fields are part and parcel of cosmological model building. They are a near essential ingredient in inflationary theories, and form the basis of most theories of, or alternatives to, dark energy. In many proposed models, the scalar is an assumed low energy field in an otherwise unknown high energy (UV) completion. In the absence of explicit UV guidance, effective field theories can be constructed according to the standard principle that every operator consistent with the underlying symmetries and field content is included in the Lagrangian. The form of the scalar low-energy effective field theory (LEEFT) is then significantly controlled by the assumed symmetry, be it exact or approximate. A special class of such LEEFT are the Galileon models [1] where the assumed global symmetry for the scalar field π is the spacetime dependent transformation π → π + c + v µ x µ . Theories of this type were discovered in the context of massive theories of gravity, originally in the Dvali-Gabadadze-Porrati model [2], where π describes the degree of freedom associated with the helicity zero mode of the massive graviton [3].
In order for a gapped (i.e. massive) scalar theory to admit a standard Wilsonian UV completion, the 2-to-2 scattering amplitude must satisfy a number of so-called positivity bounds [4,5]. These are derived based on the cherished assumptions that the scattering amplitude is Lorentz invariant, unitary, polynomially bounded in momenta, crossing symmetric and is analytic in the complex energy plane modulo certain poles and branch cuts. Of these conditions, the latter two are tied to locality and causality. Polynomial boundedness, the statement that the scattering amplitudes do not grow faster than a given polynomial (or slightly more generally a linear exponential) of complex momenta is necessary so that Fourier transforms are well defined to ensure that the amplitudes can be given meaning in real space. This is tied to locality of real space correlation functions. Analyticity is motivated by causality, and in the special cases where it can be derived rigorously, analyticity follows from the properties of the real space retarded Green's functions which are used to determine the S-matrix amplitude. Although no rigorous proof of full analyticity of the S-matrix has ever been given, it is straightforward to show that to any order in perturbation theory the scattering amplitude remains analytic, and it is generally argued that the singularities (position of poles and branch cuts) on the physical sheet in the full S-matrix are the same as those seen in perturbation theory [6].
The lowest order forward limit positivity bounds were previously used to argue that the massless Galileon [7] had no standard UV completion, because the coefficient of s 2 (s being the square of the center of mass energy) in the (pole subtracted) scattering amplitude which must necessarily be positive definite, was found to be zero [4]. This unusual behaviour is a consequence of the special soft scattering properties of Galileons, which are in turn tied to the Galileon symmetry π → π + c + v µ x µ [8,9]. This argument, however, relies on a subtle procedure which introduces a mass as an IR regulator and sends m → 0 at the end of the calculation.
The standard positivity bounds are only well-defined in the presence of a mass gap. The reason is two-fold: on the one hand a massless theory can violate the Froissart bound [10,11] which affects the number of subtractions necessary. If for example three subtractions were needed, due to the cross section growing faster than σ ∼ s (something which is technically possible for a massless theory) then it is impossible to place a bound on the sign of the coefficient of the s 2 term. Secondly, the mass gap is necessary to have an analytic region for the scattering amplitude which connects the upper and lower half complex s planes. In the massless limit, the amplitude at t = 0 is ill-defined due to a pole and branch point 1 .
Fortunately these problems are easily dealt with by considering a massive Galileon theory. In its original realization, the Galileon arose as the helicity-zero mode of a massive spin-2 resonance, and thus was neither strictly massless nor massive. It is only in a specific decoupling limit that it becomes massless. On the other hand, a massive Galileon sector is more naturally embedded in interacting spin-2 theories such as Ghost-free massive gravity [12][13][14], (see [15,16] for a recent review of massive gravity)-the scalar mode in this massive gravity theory is massive away from the decoupling limit and corresponds to the helicity-0 component of the massive graviton. Despite naïvely breaking the Galileon symmetry, a mass term leaves intact the Galileon non-renormalization theorem [17]. Furthermore, no additional operators violating the Galileon symmetry are generated at any order in loops by the presence of a mass term. It is therefore very natural to consider a massive Galileon theory, and we will show here that such a theory can indeed satisfy the positive bounds required for an analytic Lorentz invariant UV completion, providing its higher order derivative terms have suitably chosen coefficients.
Our results are consistent with the recent findings that forward limit positivity bounds are satisfied for various IR extensions of the Galileon (massive gravity [18], and pseudo-linear massive gravity / Proca theory [19]). While one might conclude from this that the non-scalar modes in Ghost-free massive gravity play an important role in satisfying the positivity requirements, based on the massive Galileon result we see that these requirements are already satisfied for the massive scalar sector alone. In fact, as we are able to apply our positivity bounds away from the forward limit [5], our requirements on the massive Galileon are stronger than equivalent bounds in the literature.
Starting with the hypothetical existence of a local, analytic Wilsonian UV completion to a Lorentz invariant massive scalar LEEFT, we have derived previously a series of positivity bounds that the (pole subtracted) 2 → 2 scattering amplitude and its derivatives with respect to the Mandelstam variables should satisfy, see Ref. [5]. These represent an infinite number of requirements that place constraints on all the coefficients of the scattering amplitude when Taylor expanded in terms of the Mandelstam variables. This Taylor expansion is always well defined given the analyticity of the pole subtracted amplitude in the Mandelstam triangle 0 ≤ s, t, u ≤ 4m 2 . These bounds are valid at and away from the forward scattering limit and as we shall see later, the more stringent bounds are often away from the forward limit. These bounds apply to any scalar EFT with a mass gap, which includes the massive Galileon on flat spacetime. The infinite number of bounds derived in [5] places constraints on the EFT coefficients, not only the leading Galileon operators, but also on all higher derivative operators that enter the massive Galileon LEEFT.
We begin by quickly reviewing the infinite number of scalar positivity bounds in section 2 before focusing on the massive Galileon Wilsonian action in section 3 which reviews a key similarity between the massive and massless Galileon, namely that they share a nonrenormalization theorem which guarantees that only Galileon invariant local terms are generated in the 1PI effective action. We then apply the positivity bounds to the leading Galileon operators in section 4 and infer a bound on the coefficients that relates the cubic and the quartic Galileon operators. We then turn to the leading higher order derivative operators in section 5, and derive their respective bounds. We prove that no local, analytic and Lorentz invariant Wilsonian UV completion could ever lead to the lowest order Galileon operators without also involving some higher derivative operators in the LEEFT, besides those being generically generated by quantum corrections,. To illustrate our results we provide in section 6 a simple yet explicit example of a Lorentz invariant and renormalizable UV theory that leads to a special massive Galileon LEEFT in the sense that the LEEFT preserves the Galileon symmetry (up to the mass term) and enjoys the Galileon non-renormalization theorem. We conclude in section 7.
We also supplement our results with four appendices: We generalize our results beyond the Galileon LEEFT to include the leading bounds for any massive and Lorentz invariant scalar LEEFT (without any assumption of symmetry) in Appendix A. In appendix B we prove that the positivity bounds are strictly positive definite (and can never be positive semidefinite in an interacting theory). In appendix C we provide an alternative (but equivalent) formulation for the coefficients and functions that enter the positivity bounds. In appendix D we provide the explicit diagonalization used in the UV example provided in section 6.

Review of the Scalar Positivity Bounds
In [5], using the unitary and analytic properties of the scalar scattering amplitude, we have derived an infinite number of positivity bounds for the derivatives of pole subtracted amplitudes. In the following, we will review the main ingredients of the proof and results. The following discussion is not specific to Galileons, but we shall apply it to the massive Galileon LEEFT in sections 4, 5 and 6 and to a general massive scalar field LEEFT in Appendix A.

Pole Subtracted Dispersion Relation
In what follows we shall be interested in the 2 → 2 scattering amplitude A(s, t), for a single scalar species of mass m which can be expanded into partial waves as where s, t, u are the usual Mandelstam variables and θ is the scattering angle in the center of mass frame, cos θ = 1 + 2t s−4m 2 . Making use of the partial wave unitarity bound Im a ≥ |a | 2 > 0 and the properties of the Legendre polynomial P , one can infer that in the physical region s > 4m 2 ∂ n ∂t n ImA(s, t = 0) ≥ 0, for n ≥ 0 .
In appendix B we prove that taken together with the assumption of analyticity this is in fact a strict positivity ∂ n ∂t n ImA(s, t = 0) > 0, for n ≥ 0 . By some weak assumptions on the analyticity on the Mandelstam complex plane, it has been shown [20,21] that A(s, t) is analytic in the twice cut s plane for fixed t and in the disk |t| < 4m 2 for fixed s (excluding the obvious poles of s and t) 2 . This leads to ∂ n ∂t n ImA(s, t) > 0, for n ≥ 0 and for all 0 ≤ t < 4m 2 , s ≥ 4m 2 .

(2.4)
On the other hand, making use of the Froissart-Martin bound [22][23][24], one can arrive at the twice subtracted dispersion relation where a(t) is some unknown function and for scalar field theories λ is independent of s and t.
In what follows we shall be working in terms of the variable v rather than the center of mass energy s, and denote by B(v, t) the pole subtracted dispersion relation (2.8) The derivatives of B(v, t) are designated by B (N,M ) , 9) and are consequently being evaluated at s = 2m 2 − t/2. Provided 0 ≤ t < 4m 2 , then 0 ≤ s ≤ 2m 2 , which lies in the Mandelstam triangle 0 ≤ s, t, u < 4m 2 in which the pole subtracted amplitudes are known to be analytic, and hence all the derivatives are well defined.

Positivity Bounds
Combining the previous expressions together with the bound (2.4) that follows from unitarity and analyticity, we have shown in [5] that the following quantity must necessarily be positive definite if the theory is to have an analytic and Lorentz invariant UV completion, where Y (2N,M ) (t) is defined by the following recurrence relation where M 2 = Min(µ+t/2−2m 2 ) = 2m 2 +t/2, and the coefficients c r and β k defined recursively by Alternative but fully equivalent expressions for these coefficients are also given in appendix C.
In what follows we shall see how to apply those bounds to the massive Galileon LEEFT.

Tree versus Loop Bounds
The bounds Y (2N,M ) (t) > 0 derived in [5] are true for the full all-loop scattering amplitude. We can however, apply them directly to the tree level LEEFT. If we compute a scattering amplitude to tree level in the low energy effective theory, then there will be no imaginary parts in the region 4m 2 ≤ µ < Λ 2 th where Λ th is the threshold to produce new heavy states, i.e. the mass of the lightest state that lies outside of the low energy effective field theory. This allows us to apply the bounds Y In the application of the tree level bounds, we must however be careful in how to interpret the bounds on the higher derivative terms. To see the problem let us first consider the forward scattering limit. If in a given theory the tree level scattering amplitude takes the form while it is clear that the forward scattering limit bounds impose c 1 > 0, we cannot further declare c 2 > 0 without further specifying how we separate trees and loops (the renormalization prescription), since on computing a one-loop diagram, we will obtain renormalization prescription sensitive local terms that will contribute at the same order 1/Λ 8 . The problem arises if there is a single scale in the problem, i.e. Λ, at which the theory is strongly coupled, then it is no longer possible to separate the tree and loop contributions that arise at the same order in a power expansion in s/Λ 2 since the loop expansion itself breaks down at s ∼ Λ 2 .
Weak coupling: As noted in [4], this problem is resolved if it is assumed that the UV theory is weakly coupled, in which case there exists an additional small parameter g for which the tree scattering amplitude takes the form then the one-loop contribution will be of order g 2 and so if g 1 we can safely put the bound on all higher derivatives of A tree . So if we assume a weakly coupled UV completion, we can consistently impose Y (2N,M ) tree (t, Λ th ) > 0 for all N and M . Interestingly, the explicit example we give for a UV completion of a massive Galileon in section 6 falls into this category.
Massive Galileon: In the case of the massive Galileon however we can do better. The tree amplitude for a massive Galileon takes the form 16) and the loop contributions 3 begin at Given the assumption m Λ, the loop corrections to the tree level coefficients computed up to and including 1/Λ 10 corrections are negligible relative to the existing tree level contribution. This holds even though the coefficient of s 2 is already m 2 /Λ 2 suppressed. As a result we may take seriously all of our tree level bounds Y (2N,M ) tree (t, Λ th ) > 0 that will be computed in sections 4 and 5 applied up to and including order 1/Λ 10 , i.e. the finite number of them that receive only contributions from the tree level scattering amplitude expanded to this order. This already allows us to put non-trivial bounds on the higher derivative coefficients without needing to assume that the UV completion is weakly coupled, and is precisely what we shall do in what follows.
Including Loops: We can also go beyond this as follows. Given the same assumption of the hierarchy m Λ, then to a given order in 1/Λ 2 there are only a finite number of loops that contribute significantly. Specifically, if we compute the scattering amplitude to order 1/Λ 2K , then we need only compute loops to order N loop = Floor[K/3]. We can then impose those loop level bounds Y (2N,M ) (t) > 0 (with M 2 = Min(µ + t/2 − 2m 2 ) = 2m 2 + t/2), that include only contributions from the scattering amplitude coefficients up to order 1/Λ 2K . These can be strengthened by using the knowledge of the light loops to that order, to compute their contribution to B (N,M ) (t) in the region in which they can be computed perturbatively. This is achieved as follows: assuming that perturbation theory can be trusted up to a scale Λ m where 1, then we may define (see [5] for the origin of this combination) Following the arguments of [5] we may then show that Y It is understood that if the amplitude is computed to order 1/Λ 2K , then only those bounds Y (2N,M ) Λ (t) > 0 that include only contributions from derivatives that arise up to this order are taken seriously.
The previous arguments ensure that for any LEEFT, we may apply the positivity bounds to any desired order in the EFT expansion, provided at least we include the light loops to the desired order. The strongest form of the bound will then be obtained by subtracting off the known contribution from the light loops to define Y (2N,M ) Λ (t). In the case of an assumed weakly coupled UV completion, which our explicit example in section 6 falls into, we can impose all orders in N, M of tree level bounds. In the special case of a massive Galileon, it is sufficient to work at tree level up to and including order 1/Λ 10 which is what we do in the following.

Massive Galileon Effective Field Theory
In what follows we start by considering the Lagrangian for a massive extension to the Galileon [7] which, in four flat spacetime dimensions, takes the form where we have used the notation Π µν = ∂ µ ∂ ν π, square brackets represent the trace of a tensor, and antisymmetrization is defined without 1/n!, e.g.
It is convenient to use the standard canonical normalization for the scalar field g 2 = 1/2. The dimensionless coefficient g 3 (or g 4 ) could be absorbed in the definition of Λ but we keep it separate for later convenience. L Galileon [π] contains the distinguished Galileon operators whose equations of motion are second order in derivatives.
Since the bounds derived in section 2 only deal with the 2 → 2 scattering amplitude, and since we will mainly focus on tree-level for π, we do not need consider higher than quartic interactions.

Wilsonian Action
When viewed as an effective field theory, the Galileon must be supplemented by an infinite number of higher derivative operators which also respect the Galileon symmetry. As a result, the full Wilsonian action for this massive Galileon LEEFT S W [π] is then, where L h.d. is a scalar function constructed from all the possible tensor combinations of two or more derivatives acting on the field. The precise form of L h.d. depends strongly on the renormalization scheme employed to compute loops since these operators are expected to receive order unity renormalizations 4 .
Stated differently, the massive Galileon may be defined as any effective field theory for which the action transforms under the transformation π → π + c + v µ x µ as The massive Galileon has several important properties that put it on the same footing as its massless counterpart 5 : 1. Quantum corrections preserve the Galileon symmetry (provided that the Galileon couples to all other fields through Galileon invariant interactions), 4 We may for instance define SW in the Euclidean as the effective action in which trees and loops of the heavy fields for all momenta, and low energy field π for momenta above the scale Λ th , are integrated out, so that the cutoff for the remaining light loop integrals is Λ th . However, a more practical definition is to define SW as the Wilsonian action for which the remaining π loops can be computed in dimensional regularization.
There is no loss in generality in this approach since the terms discarded in dimensional regularization are precisely the local operators already included in L h.d. 5 For a recent discussion on non-renormalization theorems of this type see [25].
2. Quantum corrections do not renormalize the coefficients of the leading Galileon operators in L mGal , i.e. they neither renormalize the coefficients g n , nor the mass scale m.
To demonstrate this, it is simplest to consider the expression for the one-particle irreducible effective action Γ[π], (see Ref. [17]), The 1PI effective action can be split into its classical and quantum parts, and performing a similar split in the integration measure of the path integral π = π + √ χ we then have As is well known the path integral may be evaluated to determine Γ q [π] as an expansion in powers of . The key observation is that although the addition of a mass term to S W breaks the Galilean symmetry of S W , it nevertheless leaves invariant δ n δπ n S W [π] for n ≥ 2. This is manifest since and for all higher functional derivatives the mass term does not enter, e.g.
We may thus conclude that F [π] and hence Γ q [π] is manifestly Galileon invariant This argument is easily generalized to include the Galileon interacting with other fields, provided that the interactions to the other fields are themselves manifestly invariant under the Galileon symmetry.
The second part of the non-renormalization theorem states that the Galileon operators S Galileon and the mass term are not renormalized [1,3]. This follows simply from the observation that δ 2 δπ(x)δπ(y) S W [π] and all higher derivatives depend on π only in the combination of functions of Π µν and its derivatives. For instance, for the leading massive Galileon operators 12) are explicitly only dependent on π through powers of the invariant combination Π µν . Similarly since S h.d. itself only depends on functions of Π µν and its derivatives then the same clearly holds for all functional derivatives with respect to π. This implies that the local counter terms generated in Γ q [π] will contain only local functions of Π µν and its derivatives. However, S mGal [π] contains a dependence on π through fewer than two derivatives. Thus no local counterterm generated in Γ q [π] can renormalize S mGal [π], although they will renormalize S h.d.
[π] which is precisely why the latter terms are included. Once again, the addition of a mass term does not affect this property since the mass only arises as a constant, π independent term in Mass Parameter and Pole of the Propagator: A clear product of these results is that 'mass' of the Galileon defined as the non-derivative coefficient of π 2 term in Γ[π], i.e.
, will also receive no quantum corrections. However the physical mass, defined as the pole of the propagator will receive finite quantum corrections from higher derivative terms. These occur because loops can generate terms in Γ q [π] of the Galileon invariant form π 1+n π/Λ 2n (with n ≥ 1) which will shift the physical pole m phys away from m, A small Galileon mass is technically natural as long as m Λ, which is a manifestation of the t'Hooft naturalness argument.
In practice however, when computing the scattering matrix, it is always possible to ignore higher derivative contributions in S W to the quadratic action the reason being that all such terms can be removed by a field redefinition of the form that has the virtue of preserving the Galileon symmetry. The field redefinition will have the affect of adjusting the coefficients in the Galileon invariant interactions and we can take the point of view that this has already been done from the outset.

Constraints on the Massive Galileon
Now let us apply the positivity bound to the tree level 2 → 2 scattering amplitude of the massive Galileon. Excluding the higher derivative operators, the 2 → 2 scattering amplitude for the massive Galileon (3.2) in the centre of mass frame is given by with The pole subtracted amplitude is then B(s, t) = a 00 + a 10 x + a 01 y, Significantly, for order unity g 3 and g 4 , we see that a 10 ∼ m 2 Λ 2 a 01 . This is related to the fact that the massless Galileon has enhanced soft behaviour due to the Galileon symmetry [8,9], and so in the limit m → 0 we need a 10 to vanish.
Note that these leading Galileon interactions (3.2) have given a scattering amplitude truncated at sixth order in energy, and so from the results of [4,5] we have two independent positivity bounds, Y (2,0) : a 10 + a 01t > 0 (4.6) where as mentioned earlier, Λ th is the threshold scale at which new physics must necessarily enter to restore analyticity and unitarity (since we are dealing with the tree-level amplitude). For the massive Galileon LEEFT to make sense at all, Λ 2 th should lie above the scale 4m 2 and ideally around or above the scale Λ 2 .
We can distinguish between 3 different scenarios, depending on the ratio g 4 /g 2 3 , 1. If g 4 /g 2 3 ≤ 3/4, then both bounds (4.6, 4.7) are satisfied for any value of 0 ≤ t < 4m 2 , and without any restriction neither on Λ th , nor on the mass (this implies that the Galileon mass can be taken to be arbitrarily small, without violation of these bounds).
2. For 3/4 < g 4 /g 2 3 ≤ 7/8, analyticity imposes the following upper bound on Λ th , Interestingly the strongest form of this bound arises not in the forward scattering limit t → 0, but rather in the opposite limit t → 4m 2 and the above upper bound corresponds to t = 4m 2 . This illustrates the power of extending the constraints beyond the forward scattering limit.
Even though the ratio g 4 /g 2 3 can in theory be larger than 3/4, in practise, since the LEEFT only makes sense if Λ 2 th m 2 , the ratio g 4 /g 2 3 can never get much larger than 3/4. The allowed region of parameter space is shown in figure 1. We emphasise that there is no condition imposed on g 5 , which does not contribute to the 2 → 2 tree-level amplitude.

Strong Coupling Scale
It is interesting to relate this to the scale for which perturbative unitarity breaks down for 2 → 2. This occurs when any of the partial waves violate the optical theorem, |a (s)| < 1. In this case, the largest multipole moment is at = 0, 32π s s − 4m 2 a 0 (s) = 1 −1 d cos θ P 0 (cos θ)A(s, t) (4.9) focusing again on tree-level contributions since our aim here is to compare with the bounds obtained previously from tree-level 2 → 2 scattering. Then generically the strong coupling scale implied by this process is unless we artificially tune (g 4 − 3g 2 3 /4) to be small (i.e. of order m 2 /Λ 2 or smaller). This combination of parameters has a natural explanation. In the absence of a mass term, we can use the Galileon duality transformation [26,27] to remove the cubic interaction and put it in the quartic interaction (and quintic one, which is irrelevant to this discussion). The combination g 4 = g 4 − 3g 2 3 /4 is precisely the new coefficient of the quartic Galileon operator Figure 1. Summary of constraints imposed on coefficients of cubic (g 3 ) and quartic (g 4 ) Galileon operators to respect known bounds imposed by (1). the existence of a hypothetical local, analytic UV completion, (2). a cutoff which is above the Galileon mass, and (3). the existence of a static and spherically symmetric Vainshtein mechanism. The boundary between no analytic Wilsonian UV completion and a potential UV completion with an unacceptably low cutoff is at g 4 = 7/8g 2 3 . For g 4 < 3/4g 2 3 there is to date no known obstruction for the potential existence of an analytic UV completion.
after this transformation has been performed. Thus tuning g 4 − 3g 2 3 /4 = g 4 to be small is in effect artificially making the strong coupling scale (related to tree-level 2 → 2 scattering) large by switching off interactions. Given this it is more natural to define Λ (which has so far remained a free parameter) as the strong coupling scale in the m → 0 limit, which at the level of the 2 → 2 scattering amplitude amounts to setting |g 4 − 3g 2 3 /4| = 1. With this convention we see that the bound on Λ th in the region 4/3 < g 4 /g 2 3 ≤ 7/8 is and with the usual reasonable assumption that g 3 is of order unity we would find at best Λ th ∼ m, which renders the LEEFT inconsistent. This is an example of where, as pointed out in [5], it is possible that analyticity acts as a stronger constraint on the cutoff of an effective field theory than perturbative unitarity alone.

Analyticity vs Vainshtein Mechanism
At this stage it is interesting to compare how these bounds overlap with the requirement that the Galileon exhibits a Vainshtein mechanism [7]. Although not central to their use as effective field theories, for the majority of phenomenological applications it is necessary that Galileons are in the Vainshtein screened region which suppresses their contributions to fifth forces, evading otherwise strong solar system constraints on gravity.
For instance, for a spherically symmetric source, such as a star, the Vainshtein mechanism requires that we can find a real regular solution to the static spherically symmetric Galileon equations. When working with the leading Galileon operators L Galilon , since the quintic Galileon vanishes in 3d, i.e. vanishes in any static configuration, then the quintic Galileon (or the coefficient g 5 ) does not enter this argument. From [7], see eqns.(57-60), we see that the Vainshtein only works in that static and spherically symmetric configuration if (4.14) In the special case where g 3 = 0, the positivity bound Eqn. (4.8) requires that which then excludes any possibility that a low energy effective field theory could have a unitary and analytic UV completion, with an active phenomenological Vainshtein mechanism. More generally the combined requirements imply

Positivity bounds on Higher Derivative Terms
Having shown that the massive Galileon is consistent with the leading order positivity constraints, we may look to higher order bounds Y (2N,M ) > 0. Based on the tree level computation done so far involving the operators in L mGal , we have that B (2N,M ) = 0 for all 2N +M ≥ 4, which automatically violates the higher order positivity bounds. Whilst loops from the light fields will generate a non-zero contribution to these B (2N,M ) , the tree level bounds are putting constraints on the unknown heavy physics which UV completes the Galileon. This heavy physics which has been integrated out gives rise to higher derivative operators in the EFT, previously denoted as L h.d. in (3.3), and will contribute to the tree level amplitudes and to B (2N,M ) . Consequently analyticity and unitarity of the unknown UV physics will impose constraints on the coefficients of these higher derivative operators.
For instance, to account for the leading order contributions to the 2 → 2 scattering amplitudes from this heavy physics (i.e., O(1/Λ 8 )), we must include the following cubic and quartic higher derivative interactions Each of these terms will give rise to contributions to the 2 → 2 scattering that scale as 1/Λ 8 . This is clear for the quartic interactions, and for the cubic it arises for diagrams for which one vertex is a 1/Λ 3 interaction (from the cubic Galileon) and the second is a 1/Λ 5 interaction. Hence on dimensional grounds alone, these tree level interactions will give a contribution to scattering amplitude of the form s 4 /Λ 8 which will show up as a finite contribution to the higher order Y ( 3) These are bounded by, These are the only independent bounds at this order. The higher order bounds with 2N +M > 4 cannot be computed without a knowledge of the amplitude beyond O(1/Λ 8 ). The first 2 bounds are the same as the bounds of Eq (4.6) and (4.7), but now include small corrections of O(1/Λ 8 ) from higher order derivative terms. The bound of Y (4,0) yields This bound is of course easily satisfied as we have (superficially) 4 new parameters that enter unsuppressed. In practice not all of these parameters are independent because of the ability to do field redefinitions, however the combination d 1 2 + d 4 + g 3 4 (3c 1 + 2c 2 ) is automatically invariant under field redefinitions. Crucially it is not possible to set d 1 = d 4 = c 1 = c 2 = 0.
Thus the existence of a local UV completion requires that the LEEFT has non-zero higher derivative operators. From an EFT point of view this is not too surprising since these operators will inevitably be generated from loops of the heavy fields. The new input is that already at tree level it is necessary to include these operators, i.e. it is not possible to tune the theory so that all these higher derivative terms vanish at some scale. They necessarily arise from integrating out the heavy fields that UV complete the theory. Once again, if light loops are computed in dimensional regularization, then they will make only (m/Λ) 4 suppressed contributions to the coefficients at this order.
This trend will continue if we look at higher order contributions to the amplitude, as more indices will come in and thus more possible operators at each order. For example, up to O(1/Λ 10 ) contribution, we need to include operators that schematically are of the form At each new order we will obtain new bounds, but the increase in the number of new coefficients will adequately compensate this. As we have already discussed in section 2.3, once we reach the order 1/Λ 12 then in the absence of a weak coupling parameter the one-loop of the light field contribute at the same order. The higher loops remain suppressed as long as m Λ. It is then necessary to either apply the exact version of the bounds, or follow the method discussed in section 2.3 and impose the bounds Y

Manifestly Galileon Invariant Formulation
In this section, we consider a simple UV completion of a massive Galileon, obtained via the introduction of a single heavy field H of mass M H . After integrating out the heavy field, we obtain manifestly Galileon invariant interactions for the light Galileon field π. Remembering that a massive Galileon can be defined as a theory for which under the Galileon transformation π → π + v µ x µ the Lagrangian transforms as δL = −m 2 πv µ x µ up to total derivatives, then it is straightforward to see that the following renormalizable theory respects this symmetry: Here we require |α| < 1 to avoid a ghost instability and M H m to set an appropriate EFT hierarchy that allows us to integrate out the heavy field. The Wilsonian effective action for the massive Galileon is defined via the path integral e iS W [π] = DH e iS UV [π,H] , (6.2) and will take the form of an expansion in loops of the heavy field where n counts the number of heavy loops (i.e. of loops of the heavy field H).
Explicitly integrating out the heavy field H to determine S (0) W corresponds to solving the classical equation of motion for H to give H tree , and then substituting back in the Lagrangian. This leads to As mentioned above, the higher derivative quadratic terms can be removed by a field redefinition at the price of redefining the coefficients of the interactions. To this order however the resulting interactions are relatively uninteresting since for example the operator ( π) 4 , although Galileon invariant, can be field redefined into m 8 π 4 + . . . , and so at tree level in the heavy fields there are no truly higher derivative interactions.
The situation is different if we include loops from the heavy fields. For instance, at one-loop the action picks up a contribution S (1) Expanding this we will, for example, obtain terms of the form and since H tree depends on π, H tree ∼ −α π/M 2 H + . . . then this corresponds to interactions Once again, these interactions are manifestly Galileon invariant, as required by the nonrenormalization theorem, and can be field redefined into H π 2 n π 2 + . . . ,  Figure 2.ππ →ππ diagrams up to 1-loop in the heavy field. We emphasize that from the LEEFT picture these are all tree-level diagrams of the light field.
corresponding to genuinely non-trivial higher derivative interactions. These show up in the scattering amplitude as s dependent contributions which we will see explicitly in the exact form of the scattering amplitude given below.

Diagonalized Formulation
In practice, to calculate the scattering amplitude in the UV theory, it is easier to work with an action in which both the kinetic term and mass terms are diagonalized: The explicit form of the diagonalization transformations is given in appendix D. The salient point is that for In this form the Galileon symmetry is realized in the sensẽ π →π + v µ x µ , (6.11) H →H − βv µ x µ , (6.12) Even though the heavy field shifts in this representation, since the shift is linear, the nonrenormalization theorem remains unaffected.
We are interested in calculating the 2 → 2 scattering amplitude between the light fields ππ →ππ. According to Eq (D.11), the π andπ fields are slightly different, but as a result of the equivalence theorem and the LSZ formalism the S-matrix for the two sets of asymptotic states evaluated on-shell are the same. Up to one loop in the heavy field, we only have the diagrams given in Fig. 2 for the 2 → 2 scattering ofπ. The amplitude is given by where computed using dimensional regularization (λ → µ 4−dλ ) in the MS subtraction scheme.
The closed form of the integral A X depends on the value of X. To make use of our positivity bounds, we can focus on the range 0 ≤ X < 4m 2 4M 2 , within which we have classical field configurations need to satisfy ∂ M H . Nevertheless, from this argument the weakly coupled UV completion is not in conflict with the possibility of a Vainshtein mechanism, provided the bounds derived in section 4.3 are satisfied and the gradients are under control.

Discussion
Since its first appearance within the context of (soft) massive gravity theories, the consistency of Galileon LEEFTs has remained a matter of much debate over the past decade [4]. While a mass would technically break the Galileon symmetry it does so in a way that preserves the non-renormalization theorem and all the essential features of the Galileon. Moreover in most of its known realizations, i.e. within the context of massive gravity theories, the Galileon appears as the helicity-0 mode of the graviton in a particular decoupling limit and is massive away from that decoupling limit. It is therefore natural to include the mass as part of the Galileon LEEFT.
In parallel, assuming the existence of any local, Lorentz-invariant and analytic Wilsonian UV completion imposes an infinite number of positivity bounds on the 2 → 2 scattering amplitude and its derivatives of any scalar LEEFT with a mass gap, [4,5] and are hence directly applicable to the massive Galileon LEEFT. Using all the tree-level positivity bounds known so far (both those previously derived in the literature as well as the new ones very recently derived in [5]), we have shown the existence of an entire region of parameter space which shows no obstruction (at tree-level) to the potential existence of a standard Wilsonian UV completion. A direct consequence of the positivity bounds derived in [5] is that higher derivative operators (that also respect the Galileon symmetry) are necessarily required to be present if the LEEFT is to have a standard Wilsonian UV completion.
Further requiring the existence of an active Vainshtein mechanism for static and spherically symmetric configuration does reduce this region but not entirely. However we emphasize that the analysis performed in this work has nothing to say about the validity of the Vainshtein regime where the field is strongly coupled (see [29] for a discussion on this point).
The absence of direct obstructions to the potential existence of a standard Wilsonian UV completion from 2 → 2 tree-level considerations alone, are by no means to be taken as an indication that such a UV completion will definitively exist. In the context of the leading Galileon operators (that do not get renormalized), finding such a UV completion would certainly be a success in itself, which is of course well-beyond the scope of this work. However for a particular massive Galileon LEEFT we were able to provide an explicit example of UV completion, which illustrates the fact that a Wilsonian UV completion is indeed possible and even explicitly constructible in some of these massive Galileon LEEFT.
In the context of the specific UV complete example we have found, we can manifestly see that in the massless limit, the Galileon LEEFT either becomes a trivial free theory or its cutoff vanishes. This realization is fully consistent with the results found in [4] arguing for the absence of standard Wilsonian UV completion for the massless Galileon. For that case, alternatives to the usual Wilsonian picture seem to remain as the only possibility [30][31][32][33][34][35][36].
This gives rise to the pole subtracted amplitude B(s, t) = a 00 + a 10 x + a 01 y, Note that a 10 is no longer m 2 /Λ 2 suppressed, which means that for this case the t dependence of the positivity bounds can be viewed as a small effect, as for 0 ≤ t < 4m 2 this dependence is suppressed by m 2 /Λ 2 . Truncating the amplitude to O(Λ −6 ), the Y (2,0) (t = 0) and Y (2,1) (t = 0) bounds are respectively 4c 44 − c 32 g 3 > 0 (A.10) Including a nonzero c 44 or c 32 makes it much easier to satisfy the positivity bounds for a wide range of (g 3 , g 4 ). This is not surprising, the source of tension between Galileon theories and positivity has always been that Galileon symmetry seemed to prevent operators from contributing to B (2,0) , and so discarding the Galileon symmetry naturally eases this tension.
Although these bounds are far fewer in number than the coefficients in the effective Lagrangian, this is a reflection of the fact that many of these operators are redundant to this order since they may be removed by further field redefinition. For example, where we will work to order φ 4 only. Performing this field redefinition on (A.1) we find that, again up to total derivatives and φ operates which can be removed by a further redefinition. This transformation preserves S matrix elements, and indeed we find that the amplitude (A.3) is invariant. A special case of transformations of this form is the Galileon duality [26? ] which has the additional property that it forms a group, two duality transformations taken consecutively are equivalent to a single duality transformation. The ability to perform field redefinitions represents a degeneracy in our EFT parameters, seemingly different (c nm , g n ) are describing identical theories. This degeneracy can be removed by fixing a choice of the four coefficients d nm in (A.12). For example, one can use d 20 and d 22 to set c 32 = c 44 = 0, and then d 30 or d 32 to set c 42 = 0. This leaves us with the massive Galileon (3.2), plus a π 3 and π 4 vertex, which do not contribute to the large s behaviour of the amplitude, and hence do not affect the unitarity of a UV completion. Indeed, with this choice of d mn , we find that the positivity bounds reduce to those found in the main text (4.8).
Therefore for any scalar field theory on flat space, regardless of whether or not it has Galileon symmetry, the leading order positivity bounds on the four point function can be written as the bounds (4.8) after an appropriate field redefinition.
B Proof that ∂ n t ImA(s, 0) > 0 in the physical region.
In this appendix we prove that while the optical theorem implies a semi-definite bound on the imaginary part of the coefficients to the partial wave expansion of the 2 → 2 scattering amplitude, we necessarily have a definitive positive bound on the imaginary part of the amplitude and its derivatives, rather than a semi-definite bound.

C Equivalent Expressions for the Positivity Bound
In [5], we have proven that the combinations Y ( These coefficients can also be expressed in terms of the Euler numbers E 2k and Bernoulli numbers B k , Alternatively, we may also notice that these coefficients are simply the coefficients of the Taylor expansion of sech and tan: sech(x/2) = This allows the recursive definition of Y (2N,M ) to be solved solely in terms ofB(v, t), As a function, the function D can be defined as follows, with F (x) = [tanh(x/2) + tan(x/2)]/(2x), but as an operator,D 2N is to be understood as the Taylor series expansion in t of the previous function,