Soft bootstrap and supersymmetry

The soft bootstrap is an on-shell method to constrain the landscape of effective field theories (EFTs) of massless particles via the consistency of the low-energy S-matrix. Given assumptions on the on-shell data (particle spectra, linear symmetries, and low-energy theorems), the soft bootstrap is an efficient algorithm for determining the possible consistency of an EFT with those properties. The implementation of the soft bootstrap uses the recently discovered method of soft subtracted recursion. We derive a precise criterion for the validity of these recursion relations and show that they fail exactly when the assumed symmetries can be trivially realized by independent operators in the effective action. We use this to show that the possible pure (real and complex) scalar, fermion, and vector exceptional EFTs are highly constrained. Next, we prove how the soft behavior of states in a supermultiplet must be related and illustrate the results in extended supergravity. We demonstrate the power of the soft bootstrap in two applications. First, for the N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 1 and N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 ℂℙ1 nonlinear sigma models, we show that on-shell constructibility establishes the emergence of accidental IR symmetries. This includes a new on-shell perspective on the interplay between N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 supersymmetry, low-energy theorems, and electromagnetic duality. We also show that N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 supersymmetry requires 3-point interactions with the photon that make the soft behavior of the scalar O(1) instead of vanishing, despite the underlying symmetric coset. Second, we study Galileon theories, including aspects of supersymmetrization, the possibility of a vector-scalar Galileon EFT, and the existence of higher-derivative corrections preserving the enhanced special Galileon symmetry. The latter is addressed both by soft bootstrap and by application of double-copy/KLT relations applied to higher-derivative corrections of chiral perturbation theory.


Introduction
Effective field theories (EFTs) encode the low-energy dynamics of the light degrees of freedom in a physical system. The general principle of EFTs is to include all possible local interaction terms permissible by symmetries up to a certain order in the derivative expansion. Irrelevant operators are suppressed by powers of the UV cutoff and have dimensionless Wilson coefficients that parameterize the (possibly unknown) UV physics. Of particular interest, both for formal and phenomenological applications, are the EFTs describing the low-energy interactions of Goldstone modes of spontaneously broken symmetries. Traditionally, such effective actions are constructed explicitly from the underlying symmetry breaking pattern using the method of nonlinear realization [1][2][3]. However, constructing effective actions one by one is not an efficient approach to the problem of classifying such models and studying the properties of the associated scattering amplitudes. Similar to gauge and gravity theories, the Lagrangian description of EFTs has an enormous redundancy in the form of nonlinear field redefinitions which are completely invisible in the S-matrix [4,5]. The modern on-shell approach completely avoids both the redundant description and the associated process of calculating observables from explicitly given Lagrangians. Instead one uses the required physical and mathematical properties of the on-shell scattering amplitudes to constrain the underlying models and directly calculate the physical scattering amplitudes.
The effective actions for Goldstone modes typically have the unusual property that while there may be an infinite number of gauge invariant local operators at a fixed order in the derivative expansion, the associated infinite set of Wilson coefficients is determined in terms of a finite number of independent parameters. How can this be understood in purely on-shell terms? The traditional explanation is that the spontaneously broken symmetries are nonlinearly realized on the fundamental fields and therefore mix operators in the effective action of different valence. From a more physical perspective, the spontaneously broken symmetries manifest themselves on the physical observables via low-energy or soft theorems. The non-independence of the Wilson coefficients is required to produce a cancellation between Feynman diagrams that ensures the low-energy theorem to hold. This

JHEP01(2019)195
is a redundant statement: while the number of independent parameters required to specify the effective action at a given order is reparametrization invariant, the actual Wilson coefficients are not. As we will see, from a purely on-shell perspective the collapse from an infinite number of free parameters to a finite number is a symptom of the underlying recursive constructiblility of the S-matrix, which itself can be understood as a consequence of the low-energy theorems.
It is instructive to consider an explicit example that illustrates these ideas. Consider a flat 3-brane in 5d Minkowski space. There is a Goldstone mode φ associated with the spontaneous breaking of translational symmetry in the direction transverse to the brane, and it is well-known that the leading low-energy dynamics is governed by the Dirac-Born-Infeld (DBI) action. In static gauge, it takes the form where Λ 4 is the brane tension. The action trivially has a constant shift symmetry φ → φ+c which implies that the DBI amplitudes have vanishing single-soft limits. In particular, when one of its momentum lines is taken soft, the Feynman vertex it sits on goes to zero as O( ). There are no cubic interactions, so propagators remain finite. Hence, every tree-level Feynman diagram goes to zero as O( ). What may be surprising is that a cancellation occurs between Feynman diagrams such that the soft behavior of any tree-level DBI n-point amplitude is enhanced to O( 2 ). For example for the 6-point amplitude, the O( )-contributions of the pole diagrams cancel against those of the 6-point contact term, leaving an overall O( 2 ) soft behavior: The cancellation of the O( )-contributions requires the coefficients of the 4-and 6-particle interactions (∂φ) 4 and (∂φ) 6 to be uniquely related. Interestingly we can invert the logic of this argument. Begin with the most general effective action constructed from the operators present in the DBI action, but now with a priori independent Wilson coefficients c i , schematically Imposing that the amplitudes of this model satsify O( 2 ) low-energy theorems generates an infinite set of relations among the c i . Up to non-physical ambiguities related to field redefinitions, the unique solution to these constraints is the DBI action. In that sense, DBI is the unique leading-order 4d real single-scalar theory with O( 2 ) low-energy theorems [6].

JHEP01(2019)195
The cancellation of the O( )-terms in the DBI amplitudes is a manifestation of a less obvious symmetry of the action. The broken Lorentz transformations transverse to the brane induce an enhanced shift symmetry on the brane action of the form φ → φ+c µ x µ +. . ., where the "+ . . ." stand for field-dependent terms. A theory with interaction terms built from scalar fields with at least two derivatives on every field would trivially have the enhanced shift symmetry that leads to the O( 2 ) soft behavior, but this is not the case for DBI. Therefore DBI is in a class of EFTs that have been described in previous work as exceptional [6]. This example illustrates the Lagrangian-based description of what is meant by an exceptional EFT: a local field theory of massless particles with shift symmetries that lead to an enhanced soft behavior of the scattering amplitudes beyond what is obvious from simple counting of derivatives on the fields. 1 The on-shell significance of the exceptional EFTs was first described in [7,8]. It was shown, for the case of scalar effective field theories, that the class of exceptional EFTs as defined above coincides precisely with the class of EFTs for which there exists a valid method of on-shell recursion. On-shell recursion for scattering amplitudes in the form of BCFW [9,10] or those based on various types of multi-line shifts [11][12][13][14] have been around for several years now, but they are often not valid in EFTs. Technically, this is because higher-derivative interactions tend to give 'bad' large-z behavior of the amplitudes under the complex momentum shifts and as a result there are non-factorizable contributions from a pole at z = ∞. A more physical reason is that in order for a recursive approach to have a chance, it has to be given information about how higher-point terms are possibly connected to the lower-point interactions. Standard recursion relations basically only 'know' gaugeinvariance, so in the DBI example they have no opportunity to know about any relation between the couplings of (∂φ) 4 and (∂φ) 6 . So, naturally, a recursive approach to calculate amplitudes in exceptional EFTs needs to know about the low-energy theorems, since -as illustrated for DBI -this is what ties the higher-point interactions to the lower-point ones. This is the additional input used in the early work of [15][16][17] in which vanishing soft limit arguments were used to construct the amplitudes of pion scattering without an explicit action. The idea was more recently formalized in the form of the soft subtracted recursion relations presented in [7]; they provide a tool to calculate the leading (and possibly nextto-leading) order contribution to the S-matrix of an exceptional EFT. 2 The existence of valid recursion relations gives us our sought-after on-shell characterization of the relation among the Wilson coefficients of Goldstone EFTs. The infinite set of a priori independent local operators at leading order in the derivative expansion determine the leading-order part of the S-matrix. For a generic EFT, the presence of independent operators of valence n corresponds to the appearance of independent coefficients on contact 1 This definition is a little imprecise. In standard usage, an EFT is defined by some physical data including the spectrum of particles and associated symmetries and corresponds to an effective action with operators at all orders in the derivative expansion. The defining property of an exceptional EFT however is typically only valid at leading or next-to-leading order. The equivalent on-shell statement is that the scattering amplitudes of the EFT are only recursively constructible at the same order in the expansion. 2 A recent paper [18] derives Ward identities from the enhanced shift symmetry in expectional EFTs and discusses a different type of recursion relations.

JHEP01(2019)195
contributions for amplitudes with n external particles. If the scattering amplitudes are recursively constructible at a given order, then no such independent coefficients can appear since the entire amplitude must be determined by factorization into amplitudes with fewer external particles. Furthermore, the recursion must take as its input a finite set of seed amplitudes that depend on only a finite number of parameters. Beyond being an efficient method for calculating explicit scattering amplitudes in known models, the subtracted recursion relations can be implemented as a numerical algorithm to explore and classify the landscape of possible EFTs. We term this program the soft bootstrap due to the structural similarity of the method with the conformal bootstrap [19,20]. The method is described in detail in section 3.5, here we give a simplified description. We consider EFTs as defined by a set of on-shell soft data: a spectrum of massless states, linearly realized symmetries and low-energy theorems. We use general ansätze for scattering amplitudes of low valence and low mass dimension, consistent with the assumed spectrum and linear symmetries, as input for subtracted recursion. If the ansätze satisfy a certain criterion guaranteeing the validity of the subtracted recursion relations and if the assumed soft data corresponds to a valid EFT, then the output of the recursion should correspond to a physical scattering amplitude. Here valid EFT means the existence of the assumed EFT as a local, unitary, Poincaré invariant quantum field theory.
For tree-level scattering amplitudes this includes the requirement that the only singularities of the amplitude correspond to factorization on a momentum channel. Conversely if no such valid EFT exists, or equivalently if the assumed soft data is inconsistent, then the output of the recursion generically will not correspond to a physical scattering amplitude and this may be detected through the presence of non-physical or spurious singularities. In practice, the ansätze are parametrized by a finite number of coefficients, and the removal of spurious singularities often places constraints on these coefficients.
The soft bootstrap program was initiated in [8], where it was used to explore the landscape of real scalar EFTs with vanishing low-energy theorems. The results are reviewed and extended in section 4. This paper should be understood as a continuation and generalization of this program, incorporating richer soft data including spinning particles and linearly realized supersymmetry. In section 1.1 we provide a brief overview of exceptional EFTs studied in this paper before summarizing our main results in section 1.2 that also provides an outline of the paper.

Overview of EFTs
In this paper, we extend the application of the soft bootstrap from real scalars to any massless helicity-h particle and we derive a precise criterion for the validity of the soft subtracted recursion relations. By the new validity criterion, the on-shell characterization of an exceptional EFT will precisely be that its amplitudes are constructible using soft recursion.
Our work requires a precise definition of the degree of softness of the amplitude. This is given in section 3.1. For now, let us simply introduce the soft weight σ as where S (0) n = 0. Table 1 summarizes the soft weights for various known cases of spontaneous symmetry breaking. The earlier example of DBI corresponds to the case of spontaneously broken higher-dimensional Poincaré symmetry; only the breaking of the translational symmetry actually gives rise to a Goldstone mode [21] and it will have σ = 2.
Here follows a brief overview of exceptional EFTs that appear in this paper. We include the connection between their soft behavior and Lagrangian shift symmetries: • DBI can be extended to a complex scalar Dirac-Born-Infeld theory and coupled supersymmetrically to a fermion sector described by the Akulov-Volkov action of Goldstinos from spontaneous breaking of supersymmetry. In extended supersymmetric DBI, the vector sector is Born-Infeld (BI) theory. The soft weights are σ Z = 2 for the complex scalars Z of DBI, σ ψ = 1 for the fermions of Akulov-Volkov, and σ γ = 0 for the BI photon. The soft behaviors can be associated with shift symmetries Z → Z + c + v µ x µ and ψ → ψ + ξ, where ξ is a constant Grassmann-number. 3 N = 1 supersymmetric Born-Infeld couples the BI vector to the Goldstino of Akulov-Volkov.
• Nonlinear sigma models (NLSM) describe the Goldstone modes of sponteneously broken internal symmetries and have scalars with constant shift symmetries that give σ = 1 soft weights in the low-energy theorems. A common example of an NLSM is chiral perturbation theory in which the scalars live in a coset space U(N ) × U(N )/U(N ).
The complex scalar CP 1 NLSM can be supersymmetrized with a fermion sector that is Nambu-Jona-Lasinio (NJL) model. The complex scalars have shift symmetry Z → Z + c and σ Z = 1 while the fermions have no shift symmetry and σ ψ = 0. We

JHEP01(2019)195
study both the N = 1 and 2 supersymmetric CP 1 NLSM. 4 • A NLSM can have a non-trivial subleading operator that respects the shift symmetry and hence also the low-energy theorems with σ = 1. This operator is known as the Wess-Zumino-Witten (WZW) term and has a leading 5-point interaction.
• Galileon scalar EFTs arise in various contexts and have the extended shift symmetry φ → φ + c + v µ x µ that gives low-energy theorems with σ = 2. As such they can be thought of as subleading operators of the DBI action, and are called DBI-Galileons. They can also be decoupled from DBI (at the cost of having no UV completion).
In 4d there are two independent Galileon operators: the quartic and quintic Galileon.
(By a field redefinition, the cubic Galileon is not independent from the quartic and quintic.) When decoupled from DBI, the quartic Galileon has an even further enhanced shift symmetry φ → φ + c + v µ x µ + s µν x µ x ν that gives low-energy theorems with soft weight σ = 3 and is then called the Special Galileon [8,22].
• The quartic Galileon has a complex scalar version with σ Z = 2 (but it cannot have σ Z = 3). It has an N = 1 supersymmetrization [23,24] in which the fermion sector trivially realizes a constant shift symmetry that gives σ ψ = 1.
• There is evidence [24] that the quintic Galileon may have an N = 1 supersymmetrization. This involves a complex scalar whose real part is a Galileon with σ = 2 and imaginary part is an R-axion with σ = 1.
We now summarize the main results obtained in this paper.

Outline of results
• In section 2 a brief review is given of the Wilsonian effective action. The notion of the reduced dimension of an operator is defined and the relevance to power-counting in the derivative expansion is explained.
• In section 3 we present a review and elaboration on the method of soft subtracted recursion. The asymptotic (large-z) behavior of a scattering amplitude under the momentum deformation is determined using a novel method exploiting the properties of tree amplitudes of massless particles under complex scale transformations. This result is then used to formulate a precise constructibility criterion (3.10) for the applicability of the method. The failure of an EFT (at some order in the derivative expansion) to satisfy the criterion is shown to be equivalent to the existence of independent local operators which are "trivially" invariant under an extended shift symmetry. The systematics of the soft bootstrap algorithm for constraining EFTs is described.
• In section 4 several numerical applications of the soft bootstrap are presented. The landscape of constructible EFTs with simple spectra consisting of a single massless JHEP01(2019)195 complex scalar, Weyl fermion, or vector boson is exhaustively explored. In particular, our analysis shows that there can be no vector Goldstone bosons with vanishing soft theorems. A similar result follows from an algebraic analysis that appeared around the same time as this paper [25].
• In section 5 we describe the interplay between soft behavior and supersymmetry. From the supersymmetry Ward identities we show that the soft weights of the states in an N = 1 multiplet can differ by at most one. Implications for superconformal symmetry breaking and constraints on low-energy theorems in extended supergravity are presented as examples.
• In section 6, we apply recursion to construct the scattering amplitudes of the N = 1, 2 CP 1 nonlinear sigma models at leading (two-derivative) order. For the N = 1 case, it is shown that recursive constructibility together with the conservation of U(1) charges by the seed amplitudes implies that (at two-derivative order) all tree amplitudes of this model conserve an additional accidental U(1) charge. For the N = 2 model, recursive constructibility is non-trivial due to the presence of 3-point interactions and non-vanishing scalar soft limits, but can be achieved using the supersymmetry Ward identities (see appendix D). Using this, we show that all tree amplitudes satisfy the Ward identities of SU(2) R and conserve an additional U(1) R under which the vector bosons are charged. (A detailed inductive proof of the SU(2) R Ward identities is given in appendix C.) The connection between the existence of such chiral charges for vector bosons and known results about special Kähler geometry are described, in particular we highlight the emergence of electric-magnetic duality. Finally, an explicit form of the singular low-energy theorem for the vector bosons of the N = 2 model is presented.
• Section 7 contains brief comments on supersymmetrizations of DBI and Born-Infeld.
• In section 8 various applications of the soft bootstrap algorithm to Galileon-like models are presented. Previous results on the N = 1 supersymmetrization of the quartic-and quintic-Galileon are elaborated upon, in particular the various possible soft weight assignments to the states in the multiplet are described in detail.
The existence of an extension of the special Galileon with non-trivial couplings to a massless vector is considered and evidence is given in favor of the existence of such a model. The soft bootstrap algorithm is applied to the problem of classifying higherderivative corrections to the special Galileon effective action that preserve the lowenergy theorem via the associated on-shell matrix elements. Compatible amplitudes are classified up to couplings of dimension −12 for quartic interactions and −17 for quintic interactions. These results are compared with the output of the double-copy in the form of the field theory KLT relations as applied to chiral perturbation theory. These two constructions are found to agree for quartic interactions but not for quintic.
• In appendix B many explicit forms of calculated amplitudes for various models considered in this paper are presented.

JHEP01(2019)195
2 Structure of the effective action The low-energy dynamics of a physical system can be described by a Wilsonian effective action containing a set of local quantum fields for each of the on-shell asymptotic states with all possible local interactions allowed by the assumed symmetries: (2.1) Here S 0 denotes the free theory, i.e. the kinetic terms, Λ is a characteristic scale of the problem, and c O are dimensionless constants. The sum is over all local Lorentz invariant operators O(x) of the schematic form where A, . . . , D are integer exponents. In this paper we focus on EFTs in which the operators O are manifestly gauge invariant. 5 We assign the following quantities to a local operator

JHEP01(2019)195
to determine when two operators of this form produce tree-level diagrams with couplings of the same mass dimension. Morally ρ is the same as the reduced dimension∆ [O]. The latter is the natural generalization of ρ to operators containing particles of all spins.
The quantity∆ is useful for clarifying the notion of what it means for an interaction to be leading order in an EFT with only irrelevant interactions. In the deep IR, the relative size of the dimensionless Wilson coefficients in the effective action is unimportant since lower dimension operators will always dominate over higher dimension operators. It is therefore only necessary to isolate the contributions that are leading in a power series expansion of the amplitudes in the inverse UV cutoff scale Λ −1 . The dominant interactions in the deep IR are generated by operators that minimize this quantity. As an illustrative example, consider an effective action for scalars with interaction terms of the form The reduced dimensions∆ are 2 and 5/3 for the quartic and quintic interactions respectively. The quintic interaction should therefore dominate over the quartic in the deep IR.
To see this explicitly we have to compare amplitudes with the same number of external states, so we compare the contributions from tree-level Feynman diagrams to the 8-point amplitude: This confirms that the diagrams arising from the quintic interaction dominate the 8-point amplitude.
It is useful to introduce the notion of fundamental interactions (or fundamental operators) in an EFT. These are the lowest dimension operator(s) whose on-shell matrix elements can be recursed to define all matrix elements of the theory at leading order in the low-energy expansion.
Consider the DBI action. The leading interaction comes from an operator of the form 1 Λ 4 ∂ 4 φ 4 and as discussed in the introduction, with the associated 4-point amplitude as input, all other n-point amplitudes in DBI can be constructed with soft subtracted recursion relations. If the action had contained an interaction term of the form c 5 Λ 5 ∂ 5 φ 4 , then 1 Λ 4 ∂ 4 φ 4 would not be sufficient to determine dominating contributions at n-point order, i.e. both interactions would need to be considered fundamental for soft recursion.
The operators immediately subleading to DBI in the brane-effective action are encoded in the DBI-Galileon. In 4d, there are two such independent couplings, 6 namely for a quartic interaction of the schematic form b 4 Λ 6 ∂ 6 φ 4 and a quintic interaction of the form b 5 Λ 9 ∂ 8 φ 5 ; these both have∆ = 3 whereas DBI has∆ = 2. Thus the DBI-Galileon has a total of three fundamental operators: the 4-point DBI interaction and the 4-and 5-point Galileon interactions.

Holomorphic soft limits and low-energy theorems
We rely on the 4d spinor helicity formalism (for reviews, see [28][29][30][31]) in which a massless onshell momentum is written p = −|p [p|. This presents an ambiguity in how to take the soft limit (1.
so the choice of soft limit is simply reflected in a helicity-dependent overall scaling factor. We choose to minimize the power of in the soft limit by letting the choice depend on the sign of the helicity of the particle: specfically, we take p soft → p soft = − |s [s| holomorphically for any state with non-negative helicity: 7 For a negative-helicity particle, we use the anti-holomorphic prescription |s] → |s]. For scalars, it makes no difference which choice is made. We characterize the soft behavior of amplitudes of massless particles in terms of a holomorphic soft weight σ (or, for brevity, just soft weight). It is defined in terms of the holomorphic soft limit (3.2) as where S (0) n = 0. This way of taking the soft limit is closely correlated with the shifts introduced for the soft subtracted recursion relations in the following.

Review of soft subtracted recursion relations
We consider complex momentum deformations of the form The label i = 1, 2, . . . , n runs over the n massless particles in the scattering amplitude. The shifted momentap i are on-shell by virtue of p 2 i = 0 and satisfy momentum conservation 7 Taking the soft limit as simply as in (3.2) is not compatible with overall momentum conservation.
To stay on the algebraic locus of momentum conservation in momentum space, we take the limit with appropriate shifts in a subset of the n − 1 other momentum variables. The precise prescription can be found in equation (6) of [32]. The details will not affect the main line of the discussion in this paper, but we note that all calculations are done manifestly on-shell, including the soft limits.

JHEP01(2019)195
when the shift coefficients a i satisfy the condition in (3.4). (We discuss the solutions to this condition in section 3.5.) When evaluated on the shifted momentap i , an n-point amplitude becomes a function of z and we write it asÂ n (z). The subtracted recursion relations for an n-point tree-level amplitude A n are derived from the Cauchy integral dz zÂ where the contour surrounds all the poles at finite z and the function F is defined as The vanishing of the integral in (3.5) requires absence of a simple pole at z = ∞. We derive a sufficient criterion for this behavior in section 3.3. The shift (3.4) is implemented on the spinor helicity variables according to the sign of the helicity h i of particle i as (3.7) The limit z → 1/a i is then precisely the soft limitp i → 0 of the ith particle in the deformed amplitude. Hence, if the amplitude satisfies low-energy theorems of the form (3.3) with weights σ i for each particle i, the integral (3.5) will not pick up any non-zero residues from poles arising from the function F when it is chosen as in (3.6). Therefore the only simple poles in (3.5) arise from z = 0 and factorization channels in the deformed tree amplitude. They occur where internal momenta go on-shell,P 2 I = 0. The residue theorem then states that the residue at z = 0 equals minus the sum of all such residues, and factorization on these poles gives . (3.8) The sums are over all factorization channels I, the two solutions z ± I toP 2 I = 0, and all possible particle types |ψ (I) that can be exchanged in channel I. These recursion relations are called soft subtracted recursion relations. When F = 1, the recursion is called unsubtracted.
The expression for the solutions z ± I to the quadratic equationP 2 I = 0 involves square roots, but those must cancel since the tree amplitude is a rational function of the kinematic variables. On channels where the amplitude factorizes into two local lower-point amplitudes (meaning that they have no poles), the cancellations of the square roots can be made manifest. This is done by a second application of Cauchy's theorem, which for each channel I converts the sum of residues at z = z ± I to the sum of the residues at z = 0 and z = 1/a i for all i. Details are provided in appendix A, here we simply state the result: if A (3.9) Note that this form of the recursion relation is typically only valid at low points since it requires that the amplitude factorizes into a form where all subamplitudes are local. The recursion relation in the form (3.9) is manifestly rational in the kinematic variables, and we will be using (3.9) for the applications in this paper. Note that only the first term in (3.9) has poles. Therefore the sum of the 1/a i residues over all channels must be a local polynomial in the momenta.

Validity criterion
The purpose of including F (z) in (3.5) is to improve the large-z behavior of the integrand so that one can avoid a pole at z = ∞. This is necessary in EFTs, where the large-z behavior of the amplitude typically does not allow for unsubtracted recursion relations with F (z) = 1 to be valid without a boundary term from z = ∞. A sufficient condition for absence of a simple pole at infinity is that the deformed amplitude vanishes as z → ∞. Below we show that for a theory with a single fundamental interaction (see section 2) of valence v and coupling of mass-dimension [g v ] the criterion for validity of the subtracted recursion relations is Here s i is the spin (not helicity) of particle i and σ i is its soft behavior (3.3). Alternatively, one can write the constructibility criterion in terms of the reduced dimension∆, introduced in (2.3), as The criterion generalizes to theories with more than one fundamental coupling by replacing n−2 v−2 [g v ] in (3.10) by the sum over all couplings contributing to the diagrammatic expansion of the amplitude in question; the precise criterion is given in (3.19).
Proof of the criterion (3.10). To avoid a pole at infinity in the Cauchy integral (3.5), it is sufficient to requireÂ n (z)/F (z) → 0 as z → ∞. To start with, we determine the large-z behavior of the deformed amplitudeÂ n (z).
Generically, in a theory of massless particles with couplings g k , a tree-level amplitude takes the form The mass dimension is also fixed by simple dimensional analysis to be since an n-point scattering amplitude in 4d has to have mass-dimension 4 − n.
It is useful to consider a modified scale transformation defined as The effect of this scaling can be obtained from the uniform scaling (3.13) via a little group transformation (3.1) on all momenta with t = λ 1/2 . Therefore under (3.15), M j scales as where s i is the spin (not helicity) of particle i. For the case of a theory with a single fundamental interaction of valence v with coupling g v , the number of couplings appearing in an n-point amplitude is n−2 v−2 , and therefore we have under the modified scale transformation (3.15). Under the momentum shift (3.7), the deformed tree amplitudeÂ n (z) can be written where the subscripts ± refer to the sign of the helicity of each particle. In the last line we used the behavior (3.16) under the modified scaling (3.15). At large z, the amplitude in the last line of (3.17) is the original unshifted amplitude evaluated at a momentum configuration with q i = −a i p i . These momenta are all on-shell and satisfy, via (3.4), momentum conservation. The only way the tree amplitude could have a singularity at this momentum configuration would be if an internal line went onshell. This can always be avoided for generic momenta. 8 Thus we conclude from (3.17) that for large z, the deformed amplitude behaves aŝ where D is given in (3.16). The inequality allows for the possibility that A n could have a zero at q i = −a i p i .

JHEP01(2019)195
Our mission was to find a criterion forÂ n (z)/F (z) → 0 as z → ∞. By the definition (3.6), we have F (z) → z i σ i for large z. From our analysis of the large-z behavior of A n (z), we can therefore conclude that, at worst,Â n (z)/F (z) → z D− i σ i . The sufficient criterion for absence of a pole at infinity, and hence for validity of the subtracted recursion relation, is then D − i σ i < 0. This is precisely the condition (3.10). This concludes the proof.
It is straightforward to generalize the constructibility criterion to EFTs with more than one fundamental interaction, Recall that in effective field theories, the couplings have negative mass-dimension. This means that the constructibility criterion tends to be dominated by the fundamental interactions associated with operators of the highest mass-dimension that can contribute to the n-point amplitude.
Example 1. Let us once again return to the example of DBI. The action has a fundamental quartic vertex g 4 (∂φ) 4 with a coupling of mass-dimension [g 4 ] = −4. The constructibility criterion (3.10) for the n-scalar amplitude is n(1−σ S ) < 0, where σ S is the soft behavior of the scalar φ. Since σ S = 2 in DBI, all DBI tree amplitudes are constructible via the subtracted soft recursion relations, as claimed in the introduction. The failure of the constructibility criterion for σ S = 1 is simply the statement that an EFT whose interactions are built from powers of (∂φ) 2 trivially has a constant shift symmetry and hence σ S = 1, so there are no constraints from shift symmetry on the coefficients of (∂φ) 2k in terms of that of (∂φ) 4 and then one has no chance of recursing A 4 to get all-point amplitudes.

Non-constructibility = triviality
We have derived a constructibility criterion, but what does it mean? The answer is quite simple: if an n-point amplitude can be constructed recursively from lower-point on-shell amplitudes, there cannot exist a local gauge-invariant n-field operator that contributes to the amplitude without modifying its soft behavior. We define a trivial operator to be one with at least 4 fields whose matrix elements manifestly have a given soft weight σ. Let us now assess what it takes to make an operator of scalar, fermion, and vector fields trivial.
Scalars. Operators with at least m derivatives on each scalar field will trivially have single-soft scalar limits with σ S = m.
Fermions. We have chosen the soft limit (3.2) according to the helicity such that the fermion wavefunctions do not generate any soft factors of . Thus a trivial soft behavior must come from derivatives on each fermion field in the Lagrangian. We conclude that the trivial soft behavior σ F = smallest number of derivatives on each fermion field.
Photons. Gauge invariance tells us that we should construct the interaction terms using the field strength F µν . 9 When associated with an external photon, the Feynman rule for F µν gives p µ ν − p ν µ . Naively, it may seem to be linear in the soft momentum, but under the holomorphic soft shift (3.7) it is actually O( 0 ). Recall that in spinor helicity formalism, a positive helicity vector polarization takes the form µ where q is a reference spinor. Hence, for a positive helicity photon we have This is explicitly independent of the reference spinor q because F µν is gauge invariant. For a positive helicity particle, we take the soft limit holomorphically as |p → |p (while |p] → |p]), so we explicitly see that We conclude that an operator with photons has trivial soft behavior that is determined by the smallest number of derivatives on each field strength F µν .
In an EFT where photon interactions are built only from the field strengths, the matrix elements are O(1) when a photon is taken soft. This, for example, is exactly the case for Born-Infeld theory in which the photons have σ = 0.
Constructibility. Suppose we study an n-particle amplitude with n s scalars, n f fermions, and n γ photons in an EFT whose fundamental v-particle interactions all have couplings of the same mass-dimension [g v ]. The criterion (3.10) for constructibility via subtracted soft recursion relations can be written as is the number of vertices needed at n-point.
Non-constructibility = triviality. Let us assess if there can be a local contact term for an n-particle amplitude with n s scalars, n f fermions, and n γ photons and soft behaviors σ s , σ f , and σ γ , respectively. As discussed above, a contact term that has such trivial soft behavior takes the form Or covariant derivatives Dµ = ∂µ + igAµ. In this paper, we focus on scalars and fermions that do not transform under any gauge-U(1), therefore photons must couple via Fµν .

JHEP01(2019)195
(for brevity we have not distinguished ψ andψ). In 4d, the mass-dimension of the coupling g n is easily computed as Using n = n s + n f + n γ , we can rewrite this as Compare this with (3.21); we note that the constructibility criterion is simply that or maybe more intuitively, that g n has more negative mass-dimension than n v g vvertices. So, when constructibility holds, the n-particle amplitude constructed from the n v v-valent vertices cannot be influenced by a contact term that trivially has the soft behavior: such a contact term would be too high order in the EFT due to all the derivatives needed to trivialize the soft behavior. That of course makes sense; were there such an independent local contact term, it could be added to the result of recursion with any coefficient without changing any of the properties of the amplitude. Hence recursion cannot possibly work in that case. (This is analogous to the example in [28,29] for constructibility in scalar-QED via BCFW; the difference here is that the subtracted soft recursion relations "know" about the soft behavior in addition to gauge-invariance.) The argument is easily extended to the case where the theory has fundamental vertices of different valences and mass-dimensions. We conclude that the constructibility criterion (3.10) is equivalent to the non-existence of local n-particle operators with couplings of the same mass-dimension and trivial soft behavior: Non-constructibility = Triviality.

Implementation of the subtracted recursion relations
Here we present details relevant for the practical implementation of the soft subtracted recursion relations.
Solving the shift constraints. Conservation of the momentum for the shifted momentâ p i (3.4) requires the shift variables a i to satisfy In 4d, the l.h.s. can be viewed as a 4 × n matrix p µ i of rank 4 (if n ≥ 5) multiplying a n-component vector a i . Hence the valid choices of parameters a i form a vector space given by the kernel of the matrix p µ i . For n ≥ 5 any subset of four momenta are generically linearly independent, so the p µ i -matrix has full rank. By the rank-nullity theorem, the dimension of the kernel is therefore n − 4. However, there is always a trivial solution which consists of all a i 's equal, hence non-trivial solutions to (3.25) exist only when n ≥ 6.
Practically, the linear system of equations is solved by dotting in p j , i.e. we have i s ji a i = 0 for j = 1, 2, . . . , n . (3.26) The symmetric n × n-matrix with entries s ji has rank 4, so the linear system (3.26) can be solved for say a 1 , a 2 , a 3 , and a 4 in terms of the n − 4 other a i 's.

JHEP01(2019)195
Soft bootstrap. Subtracted recursion relations can be used to calculate tree amplitudes in EFTs of Goldstone modes in theories we already know well, such as DBI, Akulov-Volkov etc. However, the soft subtracted recursion relations can also be used as a tool to classify and assess the existence of exceptional EFTs with a given spectrum of massless particles and low-energy theorems with given weights σ.
The approach to the classification of special EFTs is as follows: (1) Model input: the spectrum of massless particles and the coupling dimensions of the fundamental interactions in the model.
(2) Symmetry assumptions: the n-particle amplitudes have soft behavior with weight σ i for the ith particle.
If the constructibility criterion (3.10) is not satisfied, the assumptions (1) and (2) are trivially satisfied and we cannot constrain the couplings in the EFTs; it is not exceptional. If the constructibility criterion (3.10) is satisfied for input (1) and (2), one can use the soft subtracted recursion relations to test whether a theory can exist with the above assumptions. One proceeds as follows.
The fundamental vertices give rise to local amplitudes which must be polynomials 10 in the spinor helicity brackets, and it is simple to construct the most general such ansatz for the local input amplitudes. One can further restrict this ansatz by imposing on it the soft behaviors associated with the assumed symmetries. The result of recursing this input from the fundamental vertices is supposed to be a physical amplitude and therefore it must necessarily be independent of the n − 4 parameters a i that are unfixed by (3.25). If that is not the case for any ansatz of the fundamental input amplitudes (vertices), we learn that there cannot exist a theory with the properties (1) and (2) above. On the other hand, an a iindependent result is evidence (but not proof) of the existence of such a theory. It may well be that a i -independence requires some of the free parameters in the input amplitudes to be fixed in certain ways and this can teach us important lessons about the underlying theory. The test of a i -independence can be done efficiently numerically, and this way one can scan through theory-space to test which symmetries are compatible with a given model input.
Additionally, one can impose further constraints from unbroken global symmetries, for example, one can restrict the input from the fundamental amplitudes by imposing the supersymmetry Ward identities. We shall see examples of this in later sections.
4d and 3d consistency checks. There is a subtlety that must be addressed for n = 6. In that case, the solution space is 2-dimensional, but one solution is the trivial one with all a i equal. Furthermore, one can rescale all a i . This means that if the recursed result for the amplitude depends on the a i only through ratios of the form

JHEP01(2019)195
it will appear to be a i -independent numerically, but the result will nonetheless have spurious poles. To detect this problem numerically, we dimensionally reduce the recursed result to 3d. 11 Then the space of solutions to (3.25) is (n − 3)-dimensional, so there are non-trivial solutions and a numerical 3d test will reveal dependence on ratios such as (3.27) for n = 6. We refer to the consistency checks of a i -independence as 4d and 3d consistency checks, respectively, or simply as n-point tests when applied to construction of n-point amplitudes.
In this paper, we use 6-, 7-and 8-point tests. In section 4, we present an overview of the resulting space of exceptional pure real and complex scalar, fermion, and vector EFTs.
Special requirements for non-trivial 5-point interactions. Consider 5-particle interactions which are non-trivial with respect to a given soft behavior. This could for example be the Wess-Zumino-Witten (WZW) term, which with 4 derivatives on 5 scalars has a non-trivial σ = 1 soft behavior. Or the 5-point Galileon, which with 8 derivatives on 5 scalars has a non-trivial σ = 2. Constructibility tells us that one must be able to calculate such 5-point amplitudes from soft recursion relations via factorization, i.e. (3.28) However, there are no 3-point amplitudes available that could possibly make this work. The reason is that the only 3-scalar interaction with a non-zero on-shell amplitude is φ 3 , which gives rise to amplitudes with σ = −1 [32]. So we appear to have a contradiction: the constructibility criterion tells us that these 5-particle amplitudes are recursively constructible, but it is obviously impossible to construct them from lower-point input.
What goes wrong is that at 5-points, there are no non-trivial choices of the a i parameters that give valid recursion relations in 4d. So we have to go to 3d kinematics to resolve this issue. The above contradiction persists in 3d, so the only resolution is that these non-trivial constructible 5-point amplitudes must vanish in 3d kinematics.
Indeed they do: for WZW term and the quintic Galileon, the 5-point matrix elements are The Levi-Civita contraction makes it manifest that these amplitudes vanish in 3d. We conclude that any non-trivial (in the sense of soft behavior) 5-particle interaction must vanish in 3d. Thus, it is no coincidence that the WZW and quintic Galileon 5-point amplitudes are proportional to Levi-Civita contractions.

Soft bootstrap
We now turn to examples of how the soft recursion relations can be used to examine the existence of exceptional EFTs. The landscape of real scalar theories was previously studied in [6][7][8]14]. We outline it briefly below for completeness, but otherwise focus on

JHEP01(2019)195
new results, in particular for complex scalars, fermions, and vectors. This section considers only theories with one kind of massless particle. One can of course also couple scalars, fermions, and vectors in EFTs, and this is discussed in sections 6, 7, and 8.

Pure scalar EFTs
Consider an EFT with a single real scalar field φ. There can only be non-vanishing 3-point amplitudes in φ 3 -theory and this gives amplitudes with soft weight σ = −1. Focusing on EFTs with soft weights σ ≥ 0, the lowest-point amplitude is 4-point.
The on-shell factorization diagrams that contribute in the recursion relations (3.9) for A 6 (1 φ 2 φ 3 φ 4 φ 5 φ 6 φ ) are composed of a product of two 4-point amplitudes, for example the 123-channel diagram is 12 One sums over the 10 independent permutations corresponding to the 10 distinct factorization channels. 13 For complex scalars, we assume that the input 4-point amplitudes are of the form A 4 (1 Z 2Z 3 Z 4Z); 14 one can also consider more general input but it would not be compatible with supersymmetry, so in the present paper we do not discuss such options. At 6-point, there is only one type of amplitude that can arise from such 4-point input via recursion, and that is To get the full amplitude, one must sum over all factorization channels: In the following we consider real and complex scalar theories with 4-and 5-point fundamental vertices. 12 The momenta in the hatted amplitudes are shifted; for simplicity, we do not write the hats on the momentum variables explicitly. Note that in particular P φ should really be understood asP φ withP 2 φ = 0. 13 We do not consider color-ordering in this section. With color-ordering, one only includes the factorization diagrams from cyclic permutations of the external lines.
14 There is no color-ordering implied in any of the amplitudes here. We simply alternate Z andZ states as odd/even numbered momentum lines. In later sections, other helicity states are grouped similarly, in particular for supersymmetric cases, states that belong to the positive helicity sector sit on odd-numbered lines and negative helicity sector states on even-numbered lines. This is convenient for the practical implementation but should not be misunderstood as an indication of color-ordering.

Fundamental 4-point Interactions
Consider a theory of a single real scalar with fundamental 4-point interactions. We parameterize A ansatz 4 as the most general polynomial in the Mandelstam variables s, t, u (with s + t + u = 0) and full Bose symmetry. We subject the recursed result for A 6 to the test of a i -independence, as described in section 3.5. The result is In the table, we list the coupling dimension [g] of the fundamental quartic couplings along with the most general ansatz for the corresponding 4-point amplitude. The dash, −, indicates that the constructibility criterion (3.10) fails; this means "triviality" in the sense described in section 3.4). "F" indicates that the soft recursion fails to give an a i -independent result, and hence no such theory can exist with the given assumptions. When a case passes the 6-point test, we are able to uniquely identify which theory it is. In the above table, the non-trivial theories that pass the 6-point test are: φ 4 -theory, DBI, and the quartic Galileon. The latter automatically has σ = 3 (which is called the Special Galileon) and passes 6-point test for both σ = 2 and σ = 3. The analysis for complex scalars proceeds similarly and the results are The non-trivial theories are |Z| 4 -theory, the CP 1 NLSM (which is studied in further detail in section 6), and the complex scalar versions of DBI and the quartic Galileon. Note that there does not exist a complex scalar version of the Special Galileon with σ = 3. The results for the 6-point amplitudes of each of the theories with σ > 0 can be found in appendix B. (2) Bose symmetry, and (3) assumed soft behavior σ. In many cases, these constraints on the 5-point input amplitudes are sufficient to rule out such theories (assuming no other interactions) without even applying soft recursion. As discussed at the end of section 3.5, non-trivial 5-point amplitudes must vanish in 3d kinematics, so they are naturally written using the Levi-Civita tensor, as in the two cases of WZW and the quintic Galileon (3.29).
• 1 complex scalar. We assume input amplitudes of the form A 5 (1 Z 2Z3 Z 4Z5 Z ). Two cases pass the 8-point test.
The complex-scalar version of the quintic Galileon with [g 5 ] = −9 and σ Z = 2. The 5-point amplitude is same as for the real-scalar quintic Galileon. The fact that it passes the 8-point test is somewhat trivial: because of the two explicit factors of momentum for 4 out of 5 particles, the residues at 1/a i vanish identically for each factorization channel. The same is true for the real Galileon, so the 8-point test is not really effective as an indicator of whether such a theory may exist.
Suppose the putative complex-scalar quintic Galileon is coupled to the complex scalar DBI. Then we can conduct a 7-point test based on factorization into a quantic Galileon and a quartic DBI subamplitude. The test of a i -independence requires the coupling constant g 5 to vanish. This means that the DBI-Galileon with a complex scalar cannot have a 5-point interaction.
At [g 5 ] = −9, there is a 6-parameter family of 5-point amplitudes with σ Z = 1. The EFT with such amplitudes is generally non-constructible. However, a 1-parameter sub-family is compatible with the constraints of supersymmetry. As discussed in [24] and further in section 8.1 this may be a candidate for a supersymmetric quintic Galileon with a limited sector of constructible amplitudes.

Pure fermion EFTs
Let us now consider EFTs with only fermions and fundamental interactions of the form ∂ 2m ψ 2ψ2 . This is not the only choice, but it is the option compatible with supersymmetry. Moreover, we have found that couplings of "helicity violating" 4-point interactions in the fermion sector must vanish by the 6-point test in all pure-fermion cases we tested. The calculations proceed much the same way as for scalars, except that one must be more careful with signs when inserting fermionic states on the internal line. The diagrams needed for

JHEP01(2019)195
the recursive calculation of the 6-fermion amplitude are just like those in the scalar case (4.1), but now the permutations have to be taken with a sign: The input 4-point amplitudes are fixed by little group scaling to be 24 [13] times a Mandelstam polynomial of degree m − 1 that must be symmetric under s ↔ u to ensure Fermi antisymmetry for identical fermions. The most general input amplitudes for low values of m are summarized in the table below that also shows the result of the recursive 6-point test: We comment briefly on these results: • The NJL model has the fundamental 4-fermion interactionψ 2 ψ 2 and the result of recursing it to 6-point is given in appendix B.1. The relevance of this model will for our purposes be as part of the supersymmetrization of the NLSM (see section 6).
• Akulov-Volkov theory of Goldstinos is the only non-trivial EFT with coupling of mass-dimension −4. The Goldstinos in this theory have low-energy theorems with σ = 1. The 6-fermion amplitude is given in (B.14) in appendix B.2.
• There are no constructible purely fermionic EFTs with fundamental quartic coupling [g 4 ] = −6. Nonetheless, as was shown in [24], the quartic Galileon has a supersymmetrization with a 4-fermion fundamental interaction, however, the fermion has σ = 1, so the all-fermion amplitudes in that theory are not constructible by soft recursion: one needs additional input from supersymmetry. We refer the reader to [24] and present some further details in section 8.1.
• For [g] = −8 and σ = 2, the 6-point numerical test is passed in 4d kinematics without constraints on g and g ; that is because the recursed result depends only on ratios (3.27). When the 3d consistency check is employed, we learn that we must set g = 0 to ensure a i -independence. (This is not a strong test since the particular form of the interaction, stu, ensures that all 1/a i -poles cancel in each factorization individual diagram.) Hence, the theory that passes the 6-point test with σ = 2 has A 4 (1 ψ , 2ψ, 3 ψ , 4ψ) = g 24 [13]stu. The subtracted recursion relations fail at n > 6, which means that at 8-point and higher, this model is not uniquely determined by its symmetries. The Lagrangian construction of this theory has been studied as a fermionic generalization of the scalar Galileon [34].

Pure vector EFTs
Pure abelian vector EFTs consist of interaction terms built from F µν -contractions, possibly dressed with extra derivatives. In 4d, the Cayley-Hamilton relations imply that theories built from just field strengths F µν can be constructed from two types of index-contractions, namely (see for example [35]) If one assumes parity, the Lagrangian can only contain even powers of g. One can then write an ansatz for the Lagrangian as As established in section 3.4, a model with photon interactions built of F µν -contractions only have soft behavior σ = 0. The simplest 4-photon interactions may naively look like the vector equivalent of the constructible φ 4 scalar EFT. However, that is not the case. For the scalar, the 6-particle operator 1 Λ 2 φ 6 is subleading to the pole contributions with two φ 4 -vertices. However, for photons the pole terms with two 1 Λ 4 F 4 -vertices are exactly the same order as 1 Λ 8 F 6 . Therefore amplitudes in a theory with F n interactions and σ = 0 are non-constructible, in other words it is trivial to have σ = 0 for any choice of coefficients b i . One may ask if it is possible to choose the parameters b i in (4.9) such that the amplitudes have enhanced soft behavior σ > 0. The 6-point soft recursive test shows that this is impossible, i.e. no models exist with Lagrangians of the form (4.9) and σ > 0.
Nonetheless, the class of theories with pure F n -interactions do include one particularly interesting case, namely Born-Infeld (BI) theory. The BI Lagrangian can be written in 4d as (4.10) Upon expansion, the Lagrangian will take the form (4.9) with some particular coefficients b i . As noted, those particular coefficients do not change the single-soft behavior of amplitudes, the BI photon also has σ = 0. Nonetheless, BI theory does have the distinguishing feature of being the vector part of a supersymmetric EFT. In particular, N = 1 supersymmetric Born-Infeld theory couples the BI vector to a Goldstino mode whose self-interactions are described by the Akulov-Volkov action. One can also view Born-Infeld as the vector part of the N = 2 or N = 4 supersymmetrization of DBI. It was argued recently [35] that supersymmetry ensures BI amplitudes to vanish in certain multi-soft limits. Based on that, the BI amplitudes can be calculated unambiguously using on-shell techniques [35]. Alternatively, one can show that the N = 1 supersymmetry Ward identities uniquely fix the BI amplitudes in terms of amplitudes with Goldstinos; we discuss this briefly in section 7 and in further detail in the context of partial breaking of supersymmetry in a forthcoming paper. Next, one can consider EFTs in which the field strengths are dressed with derivatives, for example

JHEP01(2019)195
Theories with fundamental 4-point interactions are non-constructible for σ = 0 and fail the soft recursion a i -independence 6-point test for σ > 0. One implication of this is that there can be no vector Goldstone bosons with vanishing low-energy theorems. This conclusion was also reached in [25], but from a very different algebraically-based analysis. A second implication is that the pure vector sector of an N ≥ 2 Galileon model is non-constructible with the basic soft recursion, and other properties (such as supersymmetry) have to be specified in order to determine those amplitudes recursively.
There are other interesting vector EFTs: we study in detail the N = 2 supersymmetric NLSM in section 6. Furthermore, massive gravity [36][37][38] motivates the existence of a vector-scalar theory coupling Galileons to a vector field; we explore this in section 8.2.

Soft limits and supersymmetry
For models with unbroken supersymmetry, the on-shell amplitudes satisfy a set of linear relations known as the supersymmetry Ward identities [39,40]. (For recent reviews and results, see [28,29,41].) In this section, we use N = 1 supersymmetry to derive general consequences for the soft behavior for massless particles in the same supermultiplet. It is not assumed that these particles are Goldstone or quasi-Goldstone modes; the results apply to all N = 1 supermultiplets of massless particles. The consequences for extended supersymmetry are directly inferred from the N = 1 constraints.

N = 1 supersymmetry Ward identities
We consider N = 1 chiral and vector supermultiplets. We use the following shorthand for the action of the supercharges on individual particles with momentum label i: for chiral multiplets where Z is a complex scalar and ψ is a Weyl fermion. The superscripts ± refer to the helicity of the particle. Q † raises helicity by 1/2 while Q lowers it by 1/2. The prefactor is what goes outside the amplitude when the supercharge acts on it, e.g.
Due to the Grassmann nature of the supercharges, there is a minus sign for each fermion that the supercharge has to move past to get to the ith state.

JHEP01(2019)195
Similarly for a vector multiplet: where ψ is a Weyl fermion and γ is a vector boson. In this notation, the supersymmetry Ward identities are equivalent to the statement that the following action of the supercharges annihilates the amplitude [28,29,41] where L i is equal to the number of fermions to the left of Q ( †) · i and the factors P i = 0 or 1 correspond to the additional minus signs associated with the spinor prefactors as described in tables (5.1) and (5.3). Note that the action of the supercharges always changes the number of fermions by ±1, but that amplitudes are non-vanishing only if the number of fermions is even. So to get an interesting relation among amplitudes on the right-hand-side, the amplitude on the left-hand-side must vanish identically.

Soft limits and supermultiplets
We consider the chiral multiplet and vector multiplet separately and then extend the results to enhanced supersymmetry. Chiral multiplet. Define the soft factors S (i) n as the momentum dependent coefficients in the holomorphic soft expansion taken here for simplicity on the first particle The soft weights are σ Z and σ ψ for the scalar and fermion, respectively. To see how supersymmetry forces relations among the soft weights and soft factors we use (5.4) to write A n 1 + ψ , . . . , Q · i, . . . , n , where the arbitrary X-spinor cannot be proportional to |1 or |1].

JHEP01(2019)195
Taking the holomorphic soft expansion on the right-hand-side of these expressions, in the second line only, an extra power of appears in the denominator and we find The leading power of on the right-hand-side must match the leading power on the left. It is possible that cancellations among the terms on the right-hand-side may effectively increase the leading power but never decrease it. This then gives the following inequalities for which there are only two solutions These two options have different consequences for the soft factors. For σ Z = σ ψ +1, we have while for σ φ = σ ψ , we have In addition there will be an infinite number of similar relations which come from matching higher powers in .
Vector multiplet. We define the soft factors as The analysis of the supersymmetry Ward identities proceeds similarly to that of the chiral multiplet and results in only two options for the soft weights:

JHEP01(2019)195
The consequences for the soft factors are for σ ψ = σ γ + 1 13) and for σ γ = σ ψ (5.14) Note that we have made no assumptions about the sign of σ, so the relations derived here are totally general. Also, the supersymmetry Ward identities hold at all orders in perturbation theory, so the relations among the soft behaviors remain true at loop-level.
Extended supersymmetry. Relations between the soft weights of particles in the same massless supermultiplets in extended supersymmetry follow directly from the N = 1 results above, since the supersymmetry Ward identities take the same form for each pair of (s, s + 1 2 )-multiplets. In particular, the soft weights of the boson (σ B ) and fermion (σ F ) in a (s, s + 1 2 )-multiplet are related as These relations will be useful in later applications in this paper. For now, we make a small aside and demonstrate the application of (5.15) to the case of spontaneously broken superconformal symmetry and for unbroken extended supergravity.

Application to superconformal symmetry breaking
The breaking of conformal symmetry gives rise to a single Goldstone mode [21], often called the dilaton. It has been established in the literature [42][43][44] that this dilaton obeys lowenergy theorems with σ = 0. In a superconformal theory, breaking of conformal invariance must be accompanied by breaking of the superconformal symmetries. This follows from the algebra: Assuming Q-supersymmetry to be unbroken, the dilaton will be joined by a Goldstone mode from the broken R-symmetry to form a complex scalar Z with σ Z = 0. 15 It follows from our general analysis that the fermionic partner of Z will have σ = 0 or σ = −1. For

JHEP01(2019)195
the latter, Yukawa-interactions are necessary [32] and supersymmetry then requires cubic scalar interactions Z|Z| 2 + h.c. which would imply σ = −1 for the dilaton. Since σ Z = 0, σ = −1 is not possible for the dilaton and we conclude that the Goldstino mode associated with the breaking of the superconformal fermionic symmetries generated by S and S † must have low-energy theorems with soft weight σ = 0.
An example is N = 4 SYM on the Coulomb branch with the simplest breaking pattern. 16 The R-symmetry is broken from SO(6) to SO(5) and the five broken generators give rise to five Goldstone modes which join the dilaton of the conformal breaking to be the 6 real scalars of an N = 4 massless multiplet. The supermultiplet also contains the 4 Goldstinos associated with the four broken superconformal generators. The supermultiplet is capped off by a U(1) vector whose soft weight, by the above analysis, must be either σ = 0 or −1. The states that are charged under this U(1) are the massive W -multiplets and in their presence, one can have σ = −1, otherwise σ = 0 for the vector.

Application to supergravity
It is well-known that gravitons have a universal soft behavior [47]: when the soft limit (3.2) is applied to a single graviton, the amplitude diverges as 1/ 3 , i.e. the soft weight is σ 2 = −3.
(In this section, we use a subscript on the soft weight to indicate the spin of the particle.) Applying (5.15) shows that the gravitino can have σ 3/2 = −2 or −3. However, unitarity and locality constraints show [32] that amplitudes cannot be more singular than 1/ 2 for a single soft gravitino, so it must be that σ 3/2 = −2. This must be true in any supergravity theory.
Consider now a graviphoton in N ≥ 2 supergravity. Its supersymmetry Ward identities with the gravitino imply σ 1 = −2 or σ 1 = −1. The σ 1 = −2 behavior requires the graviphoton, and by supersymmetry also the gravitino, to interact with a pair of electrically charged particles via a dimensionless coupling; however, for the gravitino such a coupling is inconsistent with unitarity and locality [32]. So there is only one option, namely σ 1 = −1.
In pure N ≥ 3 supergravity, we also have spin-1 2 fermions in the graviton supermultiplet. By (5.15) and the previous results, they can have either σ 1/2 = −1 or 0. The analysis in [32] shows that σ 1/2 = −1 requires a dimensionless coupling of the spin-1 2 particle with two other particles, for example via a Yukawa coupling. Since there are no dimensionless couplings in pure supergravity, it follows from [32] that the amplitude has to be O( 0 ) or softer. This leaves only one option, namely that σ 1/2 = 0 in pure supergravity.
In pure N ≥ 4 supergravity, the scalars in the supermultiplet can have σ 0 = 0 or σ 0 = 1. If we focus on the MHV sector, the supersymmetry Ward identities give It is known that the scalar cosets of N ≥ 4 pure supergravity theories in 4d are symmetric, and therefore lead to σ 0 = 1 vanishing low-energy theorems. But at the level of the on-shell amplitudes, this conclusion does not follow from the supersymmetry Ward identities alone: as we have seen, they give σ 0 = 1 or σ 0 = 0. That analysis has to remain true at all loop-orders. In N = 4 supergravity, for example, the anomaly of the U(1) Rsymmetry can be expected to affect the soft behavior at some order. Our arguments show that it cannot happen in the MHV sector, but does not rule it out beyond MHV; this is what the σ 0 = 0 accounts for. Furthermore, one can add higher-derivative operators to the supergravity action such that supersymmetry is preserved but the low-energy theorems are not. Indeed, string theory does this in the α -expansion by adding to the N = 8 tree-level action a supersymmetrizable operator α 3 e −6φ R 4 . This operator does not affect the soft behavior of MHV amplitudes, but it is known that it does result in non-vanishing single soft scalar limits for 6-particle NMHV amplitudes at order α 3 [48,49].
The results for N = 8 supersymmetry are summarized in table 2.

MHV classification and examples of supersymmetry Ward identities
For later convenience, we state here the explicit form of the supersymmetry Ward identities (5.4) for a few particularly useful cases. We focus on the chiral multiplet, but similar results apply to the vector multiplet. First we make the simple observation that amplitudes with all Z's or only oneZ and rest Z's vanish:

JHEP01(2019)195
This follows from the supersymmetry Ward identities such as Dotting in [2| gives (5.17). Similarly A n (1 + ψ 2 − ψ 3 Z 4 Z . . . n Z = 0 and so on. In the context of gluon scattering, the equivalent statements are that amplitudes with helicity structure + + + . . . + or − + + . . . + vanish. These helicity configurations are often called "helicity violating". The simplest non-vanishing amplitudes are often denoted MHV (Maximally Helicity Violating) in the context of gluon scattering and we adapt the same nomenclature here. MHV amplitudes obey the simplest supersymmetry Ward identities in that they are just linear proportionality relations. For example, it follows from Similarly, one finds that the MHV amplitude with four fermions is proportional to the one with two fermions. To summarize, MHV amplitudes satisfy

(5.20)
The second-simplest class of supersymmetric Ward identities relate amplitudes in the NMHV class. In this paper, the 6-particle amplitudes play a central role, so we write down the 6-point NMHV supersymmetry Ward identities explicitly: We now turn to applications of these results.

JHEP01(2019)195 6 Supersymmetric non-linear sigma model
Perhaps the simplest and most familiar class of models that exhibit both linearly realized supersymmetry and interesting low-energy theorems are the supersymmetric non-linear sigma models. Of particular interest are the coset sigma models for which the target manifold is a homogeneous space G/H. At lowest order, the coset sigma model captures the universal low-energy behavior of the scalar Goldstone modes of a spontaneous symmetry breaking pattern G → H, where G and H are the isometry and isotropy groups of the target manifold respectively. If the target manifold is additionally a symmetric space and there are no 3-point interactions, then the off-shell Ward-Takahashi identities for the spontaneously broken currents imply σ = 1 vanishing low-energy theorems for the Goldstone scalars. An interesting recent perspective on coset sigma models can be found in [50]. At leading order it is fairly straightforward to calculate the on-shell scattering amplitudes for such a model from the (two-derivative) non-linear sigma model effective action. Using the methods of on-shell recursion, the use of an effective action is unnecessary. Instead, we may assume low-energy theorems and on-shell Ward identities of the isotropy group H as the on-shell data that defines the model. Using the procedure of the soft bootstrap described in section 3.5, we may apply subtracted recursion to construct the contributions to the S-matrix at leading order.
A particularly simple and well-studied example of such a construction has previously been given for the U(N )×U(N ) U(N ) coset sigma model [7,26]. There are several nice features of this model which make it an appealing toy-model to study on-shell. As will be discussed in section 8.4, at leading order (∆ = 1 or equivalently two-derivative) the isotropy U(N ) symmetry allows for the construction of flavor -ordered partial amplitudes with only (n−3)! independent amplitudes for the scattering of n Goldstone scalars.
The situation is somewhat less straightforward for models describing the low-energy dynamics of the Goldstone modes of internal symmetry breaking with some amount of linearly realized supersymmetry. 17 There are several interesting consequences of this combination of symmetries. The states must form mass degenerate multiplets of the supersymmetry algebra, which in this case means that the Goldstone scalars must always transform together with additional massless spinning states. As discussed in section 5.2, the lowenergy theorems of each of the particles in these Goldstone multiplets are not independent.
It is well-known in the literature of supersymmetric field theories that to construct a supersymmetric action, the massless scalar modes must parametrize a target space manifold with Kähler structure for N = 1 supersymmetry [51]. For N = 2 supersymmetry the target space manifold must have the structure as a non-Kähler submanifold and includes additional directions in field space or equivalently includes additional massless quasi-Goldstone scalars [53]. In general there is no unique way to extend the symmetry breaking coset to a Kähler manifold, because in any given example the spectrum of quasi-Goldstone modes depends on the details of the UV physics. Correspondingly, the quasi-Goldstone scalars do not satisfy the kind of universal low-energy theorems necessary for us to construct the scattering amplitudes recursively.
Instead, in this section we will study the interplay of low-energy theorems and supersymmetry by considering the simplest symmetric coset that is both Kähler and special-Kähler and therefore should admit both an N = 1 and N = 2 supersymmetrization. Our assumption here is that the target manifold is the coset manifold and therefore the massless spectrum should contain only two real scalar degrees of freedom, both Goldstone modes. They form a single complex scalar field Z, Z which carries a conserved charge associated with the isotropy U(1). These properties uniquely determine the Goldstone multiplets as an N = 1 chiral and N = 2 vector multiplet respectively. The main results of this section are (1) the demonstration that both the N = 1 and N = 2 CP 1 non-linear sigma models are constructible on-shell using recursion without the need to explicitly construct an effective action. And (2) this construction gives a new onshell perspective on the relationship between the linearly realized target space isotropies of M V and electric-magnetic duality transformations of the associated vector bosons.
• Scattering amplitudes satisfy isotropy U(1) Ward identities under which Z,Z are charged.
Using the approach of the soft bootstrap, we begin by constructing the most general onshell amplitudes at lowest valence that are consistent with the above data and minimizẽ ∆. There are no possible 3-point amplitudes consistent with the assumptions and so we must begin at 4-point. A |Z| 4 interaction, corresponding to∆ = 0, is consistent with U(1)

JHEP01(2019)195
conservation but violates the assumed low-energy theorem. The next-to-lowest reduced dimension interactions correspond to∆ = 1 and have a unique 4-point amplitude consistent with the assumptions Note that at 4-point, the conservation of the U(1)-charge for the complex scalar is automatically enforced as a consequence of the supersymmetry Ward identitites. We will see that this implies the conservation of the U(1) charge for amplitudes with arbitrary number of external particles corresponding to∆ = 1. Note that this is not automatic for higher order (∆ > 1) corrections and must be imposed as a separate constraint. Using (5.20) the remaining 4-point amplitudes are completely determined by supersymmetry; it is convenient to summarize the component amplitudes in a single superamplitude [54] Here we have introduced two chiral superfields Φ + and Φ − that contain the positive and negative helicity fields of the N = 1 chiral multiplet as η is the Grassmann coordinate of N = 1 on-shell superspace and η i denotes the η-coordinate of the i th superfield. We can obtain all the component amplitudes by projecting out components of the superfield. For example, the all-fermion amplitude can be derived as follows It is useful to note that the expression (6.4) is manifestly local. It follows that all component amplitudes are free of factorization singularities, indicating the absence of 3-point interactions in this theory. Note also that the pure fermion sector is exactly the NJL model detected by the soft bootstrap in section 4.2. Next, we use these 4-point amplitudes to recursively construct n-point amplitudes. Following the discussion in section 5, we note that the soft weight of the fermion must be either σ ψ = 0 or σ ψ = 1. Making the conservative choice σ ψ = 0, we evaluate the constructibility criterion on the above on-shell data, where n f is the number of external fermion states of the n-point amplitude and n s = n−n f is the number of external scalar states. For n > 4, this condition is satisfied for all npoint amplitudes. We find that recursively constructing the 6-point amplitudes yields an a i -independent expression. All the 6-point amplitudes can be found in appendix B.1.

JHEP01(2019)195
If however we make the stronger assumption σ ψ = 1, the recursively constructed 6point amplitude is a i -dependent and therefore fails the consistency checks. As a result we conclude that the true soft weight of the fermion of our theory is σ ψ = 0 and this is sufficient to construct the S-matrix at leading order from the 4-point seed amplitudes (6.4).
The recursive constructibility of the S-matrix has non-trivial consequences for the possible conserved additive quantum numbers. In a recursive model the only non-zero amplitudes are those which can be constructed by gluing together lower-point on-shell amplitudes where the states X,X on either side of the factorization channel I have CP conjugate quantum numbers. As discussed further in appendix C, if an additive quantum number is conserved by all seed amplitudes then it must be conserved by all recursively constructible amplitudes.
For example, in the present context the seed amplitudes conserve two independent U(1) charges: We know to expect the existence of an isotropy U(1) under which the scalars are charged, but from our on-shell construction it is unclear whether this should be U(1) A or a combination of U(1) A and U(1) B . We have presented the charges as two independent R-symmetries but more correctly we should consider them as a single global U(1) and a U(1) R . The presence of a second conserved quantum number is not part of the definition of the CP 1 non-linear sigma model but is instead an emergent or accidental symmetry at lowest order in the EFT. In general one would expect U(1) A × U(1) B to be explicitly broken to the isotropy U(1) by higher dimension operators.
• Scattering amplitudes satisfy isotropy U(1) Ward identities under which Z,Z are charged.
Note that, importantly, we do not impose the soft weight of the scalars σ Z = σZ = 1. As we will explain further below, no model with the above properties and vanishing scalar soft limits exists.
To proceed, interactions with reduced dimension∆ = 0 (such as Yukawa interactions) are incompatible with N = 2 supersymmetry for a single vector multiplet. Thus, the minimal value is∆ = 1; that is of course also the value for the N = 1 model. It is curious to note that N = 2 supersymmetry is sufficient to uniquely construct the S-matrix at this order in∆. As we show in the following, without assuming vanishing scalar soft limits, the restriction of the external states to a single chiral multiplet (Z,Z, ψ 1+ , ψ − 1 ) reproduces the N = 1 CP 1 sigma model.
As in the previous section, for∆ = 1 the 4-point scalar amplitude takes the form (6.3). All 4-point component amplitudes are uniquely fixed by the 4-scalar amplitudes by the N = 2 supersymmetry Ward identities and they can be encoded compactly into superamplitudes using two chiral superfields [54] Here η 1 and η 2 are the Grassmann coordinates of N = 2 on-shell superspace. The R-indices on ψ a are raised and lowered using ab , so ψ − 2 = 21 ψ 1− = ψ 1− and ψ − 1 = 12 ψ 2− = −ψ 2− . In terms of the superfields, the 4-point superamplitude can be expressed as ij η ia η ja . (6.9) We use η ia to denote the a th Grassmann coordinate of the i th external superfield. In contrast to (6.4), the superamplitude (6.9) generates component amplitudes that are not local due to the factorization singularity at P 2 13 → 0. For example, consider the following component amplitude [14] 24 [24] . (6.10) Locality and unitarity imply that this 4-point amplitude must factorize into 3-point amplitudes on the singularity at P 2 13 → 0. Denoting the helicity of the exchanged particle h,

JHEP01(2019)195
the amplitude factorizes as The contribution to the residue on the singularity takes the form with the 3-point amplitudes completely determined by Poincaré invariance and little group scaling. Comparing with the explicit form of the residue calculated from (6.10) we find that h = 1/2 and g 1 g 2 = −1. The exchanged particle of helicity h = 1/2 can be either ψ 1+ or ψ 2+ . The locality of the ) tells us that they do not factorize on the (P 13 ) 2 → 0 pole. We conclude that We carry out a similar exercise with A 4 (1 + γ 2 − γ 3 + γ 4 − γ ) for a particle of helicity h in the P 2 13 → 0 factorization channel. Comparing with the 4-point amplitude (6.9) fixes h = 0. This could correspond to either Z orZ exchange. The absence of a P 2 14 → 0 pole in 14) where g 3 g 4 = 1. Demanding that all non-local 4-point amplitudes factorize correctly fixes −g 1 = g 2 = g 3 = g 4 = −1. The 3-point superamplitudes are ij η ia η ja , (6.15) where 2 a=1 f a is defined as f 1 f 2 .

JHEP01(2019)195
It is interesting to observe that even though the N = 0, 1 and 2 CP 1 NLSM have the pure scalar 4-point amplitude in common, in the latter case the extended supersymmetry together with locality require the presence of 3-point amplitudes.
We are now in a position to address the constructibility of general n-point amplitudes. Since we are not assuming vanishing soft limits as part of our on-shell data, we are not able to make use of subtracted recursion. This is only problematic for a subset of the amplitudes in this model, at least at leading order. The unsubtracted constructibility criterion for this model reads 4 < n f + 2n v , (6.16) where n f and n v are the number of fermions and vector bosons respectively. It turns out that the amplitudes that do not satisfy this criterion can be determined from the N = 2 supersymmetry Ward identities in terms of those that do; explicit formulae are given in appendix D. Remarkably, without making any strong assumptions about the structure of low-energy theorems for the scalars, which usually characterize the sigma model coset structure, the N = 2 supersymmetry is sufficient at leading order to both construct the entire S-matrix and reproduce the amplitudes of the N = 1 and N = 0 models as special cases. This same statement can be made in the perhaps more familiar language of local field theory. At this order in the EFT expansion, the S-matrix elements should be calculable from some effective action, the bosonic sector of which should be described by a twoderivative Lagrangian of the general form where P (|Z| 2 ) and Q(|Z| 2 ) are some functions analytic around Z ∼ 0. Insisting that the Smatrix elements satisfy the on-shell N = 2 supersymmetry Ward identities is equivalent to requiring the existence of off-shell N = 2 supersymmetry transformations under which the effective action is invariant. The on-shell uniqueness result is equivalent to the statement that the off-shell N = 2 supersymmetry uniquely (up to field redefinitions) determines the form of the two-derivative effective action. In particular, the function P (|Z| 2 ) is uniquely determined to be corresponding to the Fubini-Study metric on CP 1 . Since the entire S-matrix is determined, we can explicitly demonstrate how the presence of the vector bosons modifies the structure of the low-energy theorems from the naive vanishing soft limits suggested by the coset structure. Consider the following relation among 5-point amplitudes given by the N = 2 supersymmetry Ward identities The amplitude on the right-hand-side satisfies (6.16) and therefore is constructible using unsubtracted recursion. This gives the non-constructible amplitude on the left-handside as

JHEP01(2019)195
The soft limits on particles 1, 2, 3 and 4 vanish, as expected. The soft limit on particle 5, however, is O(1), contrary to the expected soft behavior for a Goldstone mode of a symmetric coset. Explicitly It is interesting that the coupling to the photons, required by N = 2 supersymmetry, results in non-vanishing soft scalar limits for a theory with a symmetric coset. In principle, this amplitude could have had a contact contribution of the form ∝ [12] 2 , but our calculation shows that such a term would be incompatible with N = 2 supersymmetry. The maximal R-symmetry group that this model can realize is U(2) R = U(1) R × SU(2) R . We will now verify that the SU(2) R symmetry Ward identities hold for the seed amplitudes, the U(1) R we will address separately. To do this we choose a basis for the generators of SU(2) R . The scalars and vectors both transform as SU(2) singlets. The positive helicity fermion species ψ 1,2+ will transform in the fundamental representation under The negative helicity fermions transform in the anti-fundamental withT i = −T † i . This tells us that the T 0 -Ward identity is satisfied as long as the fermion species appear in pairs of (a) different helicity, same species or (b) same helicity, different species. This is true of all the non-zero amplitudes in this model. The action of T + and T − are We find that all 3-point and 4-point amplitudes in this model satisfy the SU(2) R Ward identities, for example [24] (s + t + u) = 0 . (6.24) As discussed above, we conclude that at leading order the SU(2) R Ward identities are satisfied by all amplitudes in the N = 2 model. Following the same approach as described for the N = 1 model, conservation laws satisfied by the seed amplitudes imply that the same quantities are conserved by all leadingorder amplitudes if they are recursively constructible (see appendix C). This result extends to non-Abelian symmetries, which in the on-shell language correspond to Ward identities for non-diagonal generators; this is shown for SU(2) in appendix C. The amplitudes that

JHEP01(2019)195
are not constructible using recursion are fixed by supersymmetry in terms of those that are. Therefore, they will also respect the conservation laws and non-Abelian symmetries of the seed amplitudes.
This model also conserves a separate U(1) R charge. We know to expect the conservation of the charge associated with the U(1) isotropy group. In the N = 1 case we found that the scattering amplitudes conserve an R-charge U(1) A assigned only to the complex scalar but it was consistent with the existence of U(1) B that the isotropy U(1) might also assign a charge to the fermion or even to assign equal charges in the form of a global symmetry. In the present context we also have two independent U(1) symmetries. The first is the U(1) ⊂ SU(2) R which assigns opposite charges to the fermions ψ 1+ and ψ 2+ . The second assigns charges to each of the states which, up to overall normalization can be deduced from the 3-and 4-point seed amplitudes and are summarized in the following table: These are the only linear symmetries compatible with the seed amplitudes. The isotropy U(1) must therefore be identified with some linear combination of U(1) R and U(1) ⊂ SU(2) R . This is perhaps surprising, it tells us that the massless vector boson must also be charged under the isotropy U(1). Just as for the fermions, the vector charges are chiral meaning that the positive and negative helicity states have opposite charges. Such charges for vectors are associated with electric-magnetic duality symmetries.
Such an extra U(1) R symmetry is possible because the maximal outer-automorphism group of the N = 2 supersymmetry algebra is U(2) R . The assignment of the associated charges is, up to normalization, fixed by the charge of the highest helicity state in the multiplet. It is interesting to observe that in the present context, knowledge of the non-vanishing 4-point amplitudes is insufficient to determine the U(1) R charge assignments. It is only from considering the 3-point amplitudes that we find the assignment of a non-zero chiral charge for the vector bosons unavoidable. Consider for example the amplitudes (6.14). Since the scalar is required to be charged under the isotropy U(1), which in this case must be the U(1) R since there are no other symmetries under which the scalar is charged, we see that the vector must also be charged and satisfy 2q[γ + ] = −q[Z]. The existence of fundamental 3-point interactions in this model was deduced by demanding that the singularities of the 4-point amplitudes be identified with physical factorization channels. From an on-shell point of view, it is therefore an unavoidable consequence of locality, unitarity

JHEP01(2019)195
and supersymmetry that the M V isotropy group of an N = 2 non-linear sigma model acts on the vector bosons as an electric-magnetic duality transformation.
The necessary existence of the fundamental 3-point amplitudes (6.13) and (6.14) has a further interesting consequence for the low-energy behavior of the vector boson. In [32] it was shown that singular low-energy theorems arise from the presence of certain 3-point amplitudes. In the notation used in [32] the 3-point amplitudes (6.13) and (6.14) are classified as a = 1 in the soft limit of a positive helicity vector boson. Therefore a vector boson present in amplitudes which contain at least one of the following other particles: Z, ψ a+ or γ + has soft weight σ γ = −1. Using the general formalism developed in [32], we can write down the low-energy theorem of the vector bosons in this subclass of amplitudes [sk] sk A n (1, 2, . . . , F + · k, . . . , n) + O 0 .
(6.25) Here we are using a notation similar to [55] with the introduction of an operator F + which acts on the one-particle states as and annihilates the states of the negative helicity multiplet. A similar operator F − can be defined for the soft limit of a negative helicity vector. Using equation (5.13) in conjunction with the soft behavior (6.25) of the (n + 1)-point amplitude results in the following identity for the residual n-point amplitudes where here P i = 0 or 1 corresponds to the additional signs associated with the prefactors of both the supersymmetry Ward identities and the operator F + given in table (6.26). Note that the action of Q 1 and F + commute on all physical states, so there is no ambiguity when i = j in the sums. Moreover, rearranging the order of the sums, it becomes clear that for each fixed j, the sum over i expresses a supersymmetry Ward identity for the npoint amplitudes. As such, the identity (6.27) does not impose further constraints beyond supersymmetry.

Super Dirac-Born-Infeld and super Born-Infeld
In the soft bootstrap analysis of section 4, we encountered three theories with a fundamental quartic interaction whose couplings are of mass-dimension −4: DBI, Akulov-Volkov, and Born-Infeld. These EFTs can all be related by supersymmetry. We will discuss them in further detail in future work, so for now we simply note the following:

JHEP01(2019)195
• The N = 1 supersymmetric Dirac-Born-Infeld model has as its pure scalar sector the complex scalar DBI theory with σ Z = 2 and as its pure fermion sector Akulov-Volkov theory with σ ψ = 1. All amplitudes are constructible with soft subtracted recursion. We present the expressions for the 4-and 6-point amplitudes in appendix B.2.
• The N = 1 supersymmetric Born-Infeld model combines Akulov-Volkov theory with Born-Infeld theory with σ γ = 0. All amplitudes are constructible with the soft subtracted recursion relations of section 3, except the pure vector ones, but they are uniquely fixed by the supersymmetry Ward identities. The 4-and 6-point amplitudes are given in appendix B.3.
• Extended supersymmetry binds BI, Akulov-Volkov, and DBI into one supersymmetric exceptional EFT. In particular, N = 4 supersymmetry binds together DBI, Akulov-Volkov, and Born-Infeld theory and will be discussed further in forthcoming work. The N = 4 super-DBI amplitudes can be constructed using the CHY approach [56].

Galileons
Galileons are scalar effective field theories that arise in a multitude of contexts and as a result can be defined in different ways. In 4d, Galileons are 1. Higher-derivative scalar field theories with second-order equations of motion and absence of Ostrogradski ghosts. These theories have three free parameters: the cubic, quartic and quintic interaction coupling constants. A field redefinition removes the cubic interaction in favor of a linear combination of the quartic and quintic. The scattering amplitudes are of course invariant under the field redefinition, so for the purpose of studying perturbative scattering amplitudes, we consider only the quartic and quintic Galileons.
2. The non-linear realization of the algebra Gal(4, 1) which is anİnönü-Wigner contraction of the ISO(4, 1) symmetry algebra [57]. Truncated to leading order in the reduced dimension∆, this gives an effective field theory of a real massless scalar φ with σ = 2 vanishing soft limits and coupling dimensions [g 4 ] = −6 and [g 5 ] = −9 for the quartic and quintic interactions respectively.
3. Subleading contributions to the low-energy effective action on a 3-brane embedded in a 5d Minkowski space. The leading contribution to this EFT is the DBI action and including the Galileon terms, the model is often called the DBI-Galileon. In the limit of infinite brane tension, the Galileons decouple from DBI. The non-Z 2symmetric cubic and quintic interactions arise from considering the effective action on an end-of-the-world brane.
4. Scalar effective field theories that arise from the massless decoupling limit of Fierz-Pauli-type massive gravity [36,37] and from the decoupling limit of Proca theories.

JHEP01(2019)195
It is not obvious if these definitions are equivalent. The equivalence between Definitions 2 and 3 is straightforward since ISO(4, 1) is the Poincaré symmetry of the 5d embedding space. In the brane picture of Definition 3, the DBI-Galileon scalar is a Goldstone boson that arises from the spontaneous breaking of translational symmetry transverse to the brane, with the contraction of the 5d Poincaré algebra equivalent to the non-relativistic limit of the fluctuations of the brane into the extra dimension [58].
In an approach based on scattering amplitudes, it is natural to use the second definition of Galileon theories, based on their soft weight σ = 2 and fundamental coupling dimension. This is what we do in the following, however, we do comment on the connections to the other definitions. In section 8.1, we briefly review our recent results about the supersymmetrization of (DBI-)Galileon theories in 4d and cover some details that were left out in [24]. Motivated by Definition 4, we investigate the possibility of a scalar-vector Galileon theory in section 8.2. In sections 8.3 and 8.4, we focus our attention on the Special Galileon. In section 8.3 we address the question of subleading operators respecting the enhanced σ = 3 soft behavior. In section 8.4, we approach the same question from a double-copy construction.

Galileons and supersymmetry
This section reviews and expands on the results of [24] for N = 1 supersymmetrization of Galileon models. Two approaches to forming a complex scalar Z = φ + iχ are considered: (a) Both φ and χ are Galileons so that the complex scalar Z has soft weight σ Z = 2, or (b) φ is a Galileon but χ only has constant shift symmetry; then σ φ = 2 and σ χ = 1, and hence σ Z = 1. A natural interpretation of χ is as an R-axion.
Both options were considered in [24].
Option (a): σ Z = 2. Consider first the quartic Galileon. As discussed in section 5.5, to be compatible with supersymmetry, the 4-point complex scalar amplitudes must have two Z's and twoZ's; such an amplitude is in the MHV class. It is also clear from the table of "soft bootstrap" results in (4.4) that there is a unique complex scalar quartic Galileon theory 18 with σ Z = 2 based on the 4-point interaction with A 4 (1 Z 2Z 3 Z 4Z) = g 4 stu. The other 4-point amplitudes in a supersymmetric theory are fixed by A 4 (1 Z 2Z 3 Z 4Z) using the supersymmetry Ward identity (5.20). By (5.8), the soft behavior of the fermion must be either σ ψ = 1 or 2. The all-fermion amplitudes are constructible when σ ψ = 2, and our soft bootstrap results for fermion theories (4.7) show that no such theory exists. Therefore, the fermions in a supersymmetric Galileon theory with σ Z = 2 must have σ ψ = 1.
In a supersymmetric quartic Galileon theory with σ Z = 2 and σ ψ = 1, the constructibility criterion (3.10) for n-point amplitudes with n s scalars and n f fermions is n f < 4. Thus at 6-point, we can only use soft subtracted recursion to compute the amplitudes with at JHEP01(2019)195 most two fermions. However, as discussed in [24], two of the six supersymmetry Ward identities (5.21)-(5.23) uniquely determine the 4-and 6-fermion amplitudes. The remaining four identities in (5.21)-(5.23) are used as consistency checks. The expressions for the 6point amplitudes of the supersymmetric quartic Galileon can be found in appendix B.4. We have checked that the recursively constructed 4-and 6-point amplitudes match those that we calculate from the Lagrangian superspace construction of the quartic Galileon in [23].
The supersymmetry Ward identities at 8-point and higher do not uniquely determine the non-constructible amplitudes of the supersymmetric quartic Galileon. We therefore suspect that the quartic Galileon fails to be unique at 8-point and higher [24].
The quintic Galileon does not admit a supersymmetrization with σ Z = 2 for the complex scalar. As discussed at the end of section 4.1.2, there are no obvious obstructions from the soft-recursion tests to a complex scalar decoupled quintic Galileon with However, it is not compatible with the 5-point supersymmetry Ward identities. It follows that the cubic Galileon also cannot be supersymmetrized with σ Z = 2.
Option (b): σ Z = 1. Consider a quartic complex scalar theory where the real part of the complex scalar Z is the Galileon φ and the imaginary part is an R-axion χ. The constructibility criterion with σ φ = 2 and σ χ = σ ψ = 1 is 2n χ + n f < 4, so there are only two mixed amplitudes to check; they do not restrict the 2-parameter family of input amplitudes [24]. We have checked that the constructible 6-point amplitudes are compatible with DBI. For a quintic Galileon with σ Z = 1, we found [24] a unique solution to the supersymmetry Ward identities The amplitudes A 5 (1Z 2 Z 3Z 4 Z 5Z), A 5 (1Z 2 Z 3Z 4 ψ 5ψ), and A 5 (1Z 2 ψ 3ψ 4 ψ 5ψ) follow from conjugation of the above. 19 It is interesting that the fermions in these 5-point amplitudes automatically have σ ψ = 1.
To test consistency of a supersymmetric quintic Galileon with σ φ = 2, σ χ = 1, and σ ψ = 1, we consider the 7-point and 8-point amplitudes in the decoupled Galileon theory. In both cases, the constructibility criterion is 2n χ +n f < 4. The (few) non-trivial constructible amplitudes pass the soft subtraction recursive tests of a i -independence. We have also tested compatibility with the supersymmetric DBI interactions: at 7-point the constructibility criterion is 2n χ + n f < 8 and again the constructible 7-point amplitudes pass the test. This indicates that there may indeed be a supersymmetric brane-theory with both quartic

JHEP01(2019)195
and quintic terms subleading to DBI. The scalar φ is the Goldstone mode of the broken transverse translational symmetry whereas the scalar χ is an R-axion. The fermion ψ is a genuine Goldstino of partial broken supersymmetry. We discuss such scenarios further in forthcoming work.

Vector-scalar special Galileon
It is known that scalar Galileon theories arise in certain limits of massive gravity [36,37] (for a review, see [38]). An on-shell massive graviton in 4d has 5 polarization states and the decoupling limit gives one real massless scalar (the Galileon) and a massless photon in addition to the massless graviton. So we expect there to be an EFT of a real Galileon scalar coupled to vector. 20 The vector couples quadratically to the scalar and was consistently truncated off in [37]. Some subsequent studies have discussed the photon-scalar coupling of Galileons, see for example [59]. Here, we use soft recursion to give some definitive results about the possible scattering amplitudes in such a theory.
If the scalar has σ φ = 2, only the scalar amplitudes are constructible, and we are not able to say anything about the vector sector and its couplings to the scalar. If however the couplings are tuned in such a way that the cubic and quintic Galileon interactions are set to zero then in the scalar sector the soft weight of the scalar is enhanced to σ φ = 3, the special Galileon scenario. At present it is unknown whether this enhancement of symmetry can be understood in some natural way from the decoupling limit of some model of massive gravity. Moreover, it is not a priori clear if the σ φ = 3 enhancement can survive coupling to other particles.
We use the power of the soft bootstrap to construct the most general amplitudes consistent with the special Galileon low-energy theorem. We use the 6-point test to exclude EFTs with a special Galileon coupled non-trivially to a photon with σ γ > 0. For the model with σ φ = 3 and σ γ = 0, we find that the soft recursion 6-point test reduces the most general 6 real-parameter ansatz for the scalar and scalar-vector interactions to a 3 real-parameter family: The couplings of the pure vector sector are unconstrained; the most general ansatz is The most interesting feature of the above result is the relation between the coefficients of the amplitudes A 4 (1 φ 2 φ 3 φ 4 φ ) and A 4 (1 − γ 2 φ 3 φ 4 + γ ). The former is the familiar quartic 20 The decoupling of these interactions from the graviton is not clear [38].

JHEP01(2019)195
Galileon, while the latter would arise from an operator of the form where F ± are as defined in and below (3.20) The relation between the couplings strongly indicates the existence of a non-linear symmetry which mixes the scalar and vector modes. Describing the action of this symmetry and its consequences is left for future work.

Higher derivative corrections to the special Galileon
The real quartic Galileon has low-energy theorems with σ = 3 soft weight. Being agnostic about the origin of the special Galileon, from an EFT perspective, one should write a Lagrangian with all possible operators that respect the symmetries of the theory in a derivative expansion. The authors of [60] found that among a specific subclass of Lagrangian operators, namely those with the schematic form ∂ 4 φ 4 , ∂ 6 φ 4 and ∂ 8 φ 5 , the special Galileon is the unique choice that can give enhanced soft limits with σ = 3 soft weight. In this section, we investigate much more exhaustively the possible higher-derivative quartic and quintic operators compatible with σ = 3 soft behavior. This is done using soft-subtracted recursion relations to calculate the 6-and 7-point scattering amplitudes of the model. Taking into account Bose symmetry, the most general ansatz one can write down for the 4-point matrix element of local operators is .
The leading term with coupling c 0 /Λ 6 is the usual quartic Galileon. The terms suppressed by higher powers of the UV cutoff Λ encode all possible higher-derivative quartic operators of the scalar field up to order Λ −14 .
We apply the 6-point test with σ = 3 and find that consistency requires c 1 = c 3 = 0 in the ansatz (8.7). The 4-point amplitude then becomes

JHEP01(2019)195
From this, we understand that there cannot exist an 8-derivative Lagrangian operator that preserves the special Galileon symmetry. Additionally, at 6-, 10-and 12-derivative order there exist unique quartic operators compatible with σ = 3. In section 8.4, we show explicitly that the result (8.8) can also be obtained from an application of the BCJ double-copy. Next we examine the possible existence of quintic operators compatible with σ = 3. We combine input from the quartic Galileon with the most general possible ansatz for the 5-point matrix elements and use the 7-point test to assess compatibility with σ = 3. The soft subtracted recursion relations at 7 points are valid if Since the 7-point coupling constant is the product of a quartic (with mass dimension −6 or lower) and a quintic coupling, the latter must then satisfy With Bose symmetry and the requirement that the ansatz for the 5-point amplitude must have soft weight σ = 3, we are left with In the above, (1234) = µνρσ p µ 1 p ν 2 p ρ 3 p σ 4 , the sum i<j means 4 i=1 5 j=i+1 , while the sum P is over all permutations of {1, 2, 3, 4, 5}, (−1) |P | is the signature of the permutation and P i is its ith element. There are no contributions to the amplitude that have less than 14 derivatives. The 1/Λ 14 -term satisfies the constructibility criterion and vanishes in 3d kinematics, in agreement with the discussion of section 3.5. Two of the 1/Λ 17 -terms also vanish in 3d kinematics, but this was not a priori expected since they are too high order to satisfy constructibility.
The 7-point test implies no constraints on the coefficients d 1 , d 2 , d 3 and d 4 . This is evidence in favor of the existence of four 5-point operators that preserve the special Galileon symmetry. Next, in section 8.4, we investigate whether this result can be obtained from a double-copy prescription, similar to the 4-point case.

Comparison with the field theory KLT relations
The significance of the special Galileon extends well beyond the contraction limit of the 3brane effective field theory and the decoupling limit of massive gravity. The enhancement of the soft behavior to σ = 3 (which degenerates to σ = 2 when the DBI interactions are re-introduced) or correspondingly the extension of the non-linearly realized symmetry algebra suggests that this model has a fundamental significance of its own that is at present JHEP01(2019)195 only partially understood. Perhaps one of the deepest and least understood aspects of the special Galileon is its role in the (field theory) KLT algebra as the product of two copies of the U(N )×U(N ) non-linear sigma model. For N = 2, 3 this coset sigma model has been intensively studied as a phenomenological model of the lightest mesons under the name Chiral Perturbation Theory (χPT). Henceforth we will use this name to avoid confusion with the CP 1 non-linear sigma model discussed in section 6.
The double-copy relation between χPT and the special Galileon was first understood in the CHY auxilliary world-sheet formalism [61]. Specifically, it was shown in the CHY formalism that the leading order contribution to scattering in the special Galileon model can be obtained from the KLT product where α, β index the (n − 3)! independent color(flavor)-orderings. 21  where P(2, 3, 4) denotes the sum of all permutations of legs 2, 3 and 4. For the formulae (8.12) and (8.13) to even be well-defined, the color-ordered amplitudes on the right-hand-side must satisfy a number of non-trivial relations to reduce the number of independent partial amplitudes to (n − 3)! for the scattering of n particles. The existence of a color-ordered representation is itself non-trivial and not guaranteed to be satisfied in all models with color structure [62]. In all known cases where the double-copy relations (8.12) give a sensible, physical output, the reduction to a reduced basis of size (n − 3)! is accomplished by two sets of identities among the partial amplitudes, namely the Kleiss-Kuijf and fundamental Bern-Carrasco-Johansson relations. That these identites obtain for amplitudes calculated in the leading two-derivative action of χPT was first established in [63] using semi-on-shell recursion techniques developed in [64].
Our goal in this section is to connect two (possibly discrepant) definitions of the special Galileon model: 1. The special Galileon is the most general effective field theory of a real massless scalar with σ = 3 vanishing soft limits.
2. The special Galileon is the double-copy of two copies of χPT. 21 We use square brackets for the arguments of a color-ordered amplitude.

JHEP01(2019)195
What we have described above is the known fact that these definitions agree at the lowest non-trivial order. In the previous section we used soft subtracted recursion to construct the most general 4-and 5-point amplitudes consistent with the first definition up to order Λ −12 and Λ −17 respectively. To determine if these results agree with the second definition we must first construct the most general 4-and 5-point amplitudes in χPT compatible with the requirements of the double-copy. Here we are following the approach of [62] and making the most conservative possible assumptions. Specifically we assume that both the explicit form of the double-copy (8.13)  The superscript k counts both the mass dimension of the function and the number of derivatives in the underlying effective operator. In this language, the double-copy-compatibility conditions take the form C: We make a general parametrization of the polynomial functions as      A few comments about this result. As expected, the leading 2-derivative contribution is compatible with the conditions (8.18). Surprisingly, there are no compatible contributions from 4-derivative operators, but there are unique contributions at 6-and 8-derivative order. Moreover, the structure of the result here agrees with the 4-point amplitude of Abelian Z-theory [65]. The Z-theory model is a top-down construction which gives open string scattering amplitudes as the field theory double-copy of Yang-Mills and a higher-derivative extension of χPT. The Z-amplitudes are by construction guaranteed to satisfy the doublecopy-compatibility conditions but with Wilson coefficients g i having precise values calculated from the known string amplitudes. The method of this section can be understood as the bottom-up converse of the Z-theory construction, and at 4-point we find agreement. To summarize, we have shown that up to 8-derivative order there is a 3-parameter family of operators that generate 4-point matrix elements compatible with the conditions required for the double-copy to be well-defined. We could continue this to higher order, but our ability to compare with the methods of section 8.3 are bounded above at this order by the constructibility criterion.
To construct the associated amplitudes in the special Galileon model (according to the second definition described above) we use the first relation in (8.13). The result is (8.21) in precise agreement with the special Galileon amplitude (8.8).
As an additional check to the results obtained above, we calculate the 6-point amplitudes of both χPT and the special Galileon. Up to order O(Λ −6 ) the χPT amplitude can be calculated using soft subtracted recursion with (8.20) as input. Note that only three factorization channels contribute to this calculation because the rest do not preserve color ordering. The resulting amplitude, We are now in position to calculate the 6-point special Galileon amplitude with two different methods. We can either use the 6-point KLT relation in (8.13) or use soft subtracted recursion with (8.21) as input. The results of these calculations match perfectly up to order O(Λ −18 ), which is the furthest the recursive calculation can go.

JHEP01(2019)195
Shifting our focus to 5-point amplitudes, we find that it is not possible to reproduce (8.11) as a double-copy of two (identical or non-identical) color-ordered scalar amplitudes, despite the perfect agreement at 4-and 6-points. Starting from a general ansatz for the scalar color-ordered amplitude, we find that the leading contribution that satisfies all C, KK and BCJ constraints is O(Λ −15 ) corresponding to a valence 5 scalar-field operator with 14 derivatives. The existence of such an operator at all is interesting since there are apparently no odd point amplitudes in Z-theory [65]! At this order we find that the kinematic structure of Z-theory does not coincide with the most general possible double-copycompatible higher-derivative extension of χPT. Or perhaps said differently, just like string theory fixes the Wilson coefficients in the 4-point result (8.20) to take particular (non-zero) values, it appears to fix the Wilson coefficients of the odd-point amplitudes to be zero.
When we use the second relation of (8.13) with this result, we obtain a 5-point scalar amplitude of order O(Λ −33 ), which is significantly subleading to the amplitude (8.11) we calculated in the previous section for the special Galileon.

Outlook
There are several interesting questions that remain unanswered in this work. In section 4 we applied the soft bootstrap to classes of models with simple spectra consisting of a single particle of a particular spin. Furthermore, we gave a limited examination of classes of models with linearly realized supersymmetry with spectra consisting of a single multiplet. There is a potentially vast landscape of constructible models with more complicated spectra and possible futher interesting linearly realized symmetries.
We have already seen examples of this; in section 6 further symmetry (in this case electromagnetic duality symmetry) emerges as an unavoidable consequence of the combination of low-energy theorems and linear N = 2 supersymmetry. Similarly we should expect the soft bootstrap to reveal models with complicated non-linear symmetries. In section 8.2 we have given evidence in favor of the existence of such a symmetry underlying a vector-scalar extension of the special Galileon.
Our results also suggest two additional applications for the soft bootstrap. The first is to the classification of higher-derivative operators. The method applied in sections 8.3 and 8.4 to the special Galileon and χPT is generalizable to a large class of EFTs with manifest advantages over traditional methods. The second is a useful cross-check on results concerning exceptional EFTs obtained via the double copy. In section 8.4 we found the puzzling result that there exist valence 5 operators invariant under the special Galileon symmetry which apparently cannot be constructed as the double copy of subleading χPT operators.
It would be reasonable to expect further, similarly rich and unexpected, phenomena to be present throughout the landscape of constructible EFTs. Physics, UC Santa Barbara, for hospitality during the 'Scattering Amplitudes and Beyond' program which was supported under Grant No. NSF PHY17-48958 to the KITP. This work was supported in part by the US Department of Energy under Grant No. DE-SC0007859. CRTJ was supported by a Leinweber Graduate Fellowship and MH by a Rackham Predoctoral Fellowship from the University of Michigan.
A Derivation of (3.9) In this appendix, we derive the manifestly local form (3.9) of the subtracted recursion relations. For a given factorization channel, consider from the recursion relations (3.8) the expression where the contour surrounds only the two poles z ± I . The second equality is non-trivial and deserves clarification. In the second expression, the subamplitudesÂ Let us assess the large-z behavior of the integrand in (A.1). The L and R subamplitudes have couplings g L and g R such that g L g R = g n , with g n the coupling of A n . Their massdimensions are related as [g L ] + [g R ] = [g n ]. Hence, using n L + n R = n + 2 and (3.16), we find that the numerator behaves at large z aŝ where s P denotes the spin of the particle exchanged on the internal line and D is the large z behavior of the A n which we know satisfies D − n i=1 σ i < 0, by the assumption that the amplitude A n is recursively constructible by the criterion (3.10). We therefore conclude that the integrand in (A.1) behaves as z D−1− n i=1 σ i −2s P , i.e. it goes to zero as 1/z 2 or faster. Hence, there is no simple pole at z → ∞.
If we deform the contour, we get the sum over all poles z = z ± I inÂ where F (z) = n i=1 (1 − a i z) σ i . This form of the recursion relation is manifestly rational in the momenta.
Note that only the z = 0 residues give pole terms in A n . Therefore the sum of the 1/a i residues over all channels must be a local polynomial in the momenta. For example, it is valid for the reconstruction of the 6-point scalar amplitude of NLSM, but not for the reconstruction of an 8-point amplitude.

B Explicit expressions for amplitudes
In this appendix, we present expressions for the 4-and 6-point amplitudes of the theories discussed in the main text. The 6-point amplitudes were reconstructed with the 4-point ones as input, by means of the subtracted recursion relations and the supersymmetry Ward identities also discussed in the main text.

B.1 Supersymmetric CP 1 NLSM
Below, we list the amplitudes for the CP 1 N = 1 supersymmetric NLSM. This model is discussed in section 6 as an illustration of our methods.
The 4-point amplitudes are: They serve as the input for computing the 6-point amplitudes recursively:

B.2 Supersymmetric Dirac-Born-Infeld theory
The amplitudes of N = 1 supersymmetric Dirac-Born-Infeld theory are all recursively constructible. The 4-point amplitudes are  In this subsection, we list the amplitudes of Born-Infeld theory. This theory is the leading order contribution to the effective field theory of a Goldstone N = 1 vector multiplet. The 4-point amplitudes are In this case, only A 6 (1 + γ 2 − γ 3 + γ 4 − γ 5 + ψ 6 − ψ ) has local terms.

B.4 Supersymmetric quartic Galileon theory
Below, we list the amplitudes of an N = 1 supersymmetric quartic Galileon. This model was discussed in detail in [24] and reviewed in section 8. The 4-point amplitudes are

JHEP01(2019)195
Now consider Ward identities generated by elements of a semi-simple Lie algebra. In the root space decomposition of the algebra, we can choose a triplet of generators: raising operators T + , lowering operators T − , and "diagonal" T 0 generators, for each positive root that satisfy the algebra In order for representations of this algebra to be physical, CPT must be an algebra automorphism. The CPT charge conjugation generator C must also flip the sign of the additive T 0 -charge. So we determine the action of C to be C · T 0 · X = −T 0 · C · X = −T 0 ·X , where X is a physical state and we have defined the conjugate stateX to be the charge conjugate of X, i.e.X = C · X.
If the S-matix is recursively constructible (at some order in the derivative expansion) then each n-point amplitude is given as a sum over factorization singularities with residues given in terms of a product of amplitudes with fewer external states where I labels all possible factorization channels and X the exchanged internal states. Since T 0 is diagonal, the Ward identity generated by T 0 works just like in the Abelian case -charges can be assigned to the physical states and recursion preserves this charge in any n-point amplitude. More complicated are the non-diagonal generators T ± . For simplicity, we present the argument explicitly for SU(2) R Ward identities as they apply to the N = 2 NLSM described in section 6.2. For SU(2) R , the action of T + on the fermion helicity states is given in (6.23). The scalar and vectors are singlets under SU(2) R . The statement of the SU(2) R Ward identity is that T + ·A n (1, . . . , n) = 0. The inductive assumption is that this holds true for the lower-point amplitudes in the recursive expression for A n (1, . . . , n). We already know from section 6.2 that SU(2) R is a symmetry of the 3and 4-point amplitudes, so that provides the basis of induction.
The action of T + on the recursive expression for an n-point amplitude is T + · A n (1, . . . , n) ≡ where P i = 0 or 1 corresponds to the additional signs in the prefactors for the action of T + as given in table (6.23). We now prove that this expression vanishes channel by channel. Without loss of generality, we will show that the contribution from the (1 . . . k) ± channel vanishes independently, where + means the contribution from the z ± residue. The argument follows for all other factorization channels by replacing (1 . . . k) ± by I ± . For the (1 . . . k)-channel, the relevant part of (C.4) that we want to show vanishes is and similarly forÂ R . Using this relation and splitting the sum over particles X allows us to rewrite (C.6) as − X (−1) P X Â L (1, . . . , k, T + · X)Â R (X, k + 1, . . . , n) − X (−1) PX Â L (1, . . . , k, X )Â R (T + ·X , k + 1, . . . , n) .
(C. 10) where Q X refers to the prefactors for the action of T − as given in table (6.23). This vanishes when T + · T − ·X =X and P T − ·X + QX = 0 for any state X such that T + · X = 0. For SU(2) R , we can check explicitly that these conditions are satisfied. The only states