Gluon scattering on self-dual radiative gauge fields

We present all-multiplicity formulae, derived from first principles in the MHV sector and motivated by twistor string theory for general helicities, for the tree-level S-matrix of gluon scattering on self-dual radiative backgrounds. These backgrounds are chiral, asymptotically flat gauge fields characterised by their free radiative data, and their underlying integrability is captured by twistor theory. Tree-level gluon scattering scattering amplitudes are expressed as integrals over the moduli space of holomorphic maps from the Riemann sphere to twistor space, with the degree of the map related to the helicity configuration of the external gluons. In the MHV sector, our formula is derived from the Yang-Mills action; for general helicities the formulae are obtained using a background-coupled twistor string theory and pass several consistency tests. Unlike amplitudes on a trivial vacuum, there are residual integrals due to the functional freedom in the self-dual background, but for scattering of momentum eigenstates we are able to do many of these explicitly and even more is possible in the special case of plane wave backgrounds. In general, the number of these integrals is always less than expected from standard perturbation theory, but matches the number associated with space-time MHV rules in a self-dual background field, which we develop for self-dual plane waves.


Introduction
There are many physical scenarios which are well approximated by perturbative quantum field theory (QFT) in strong background fields: fixed, non-trivial solutions to the classical equations of motion which are treated non-perturbatively [1][2][3][4]. In the context of gauge theories, strong Maxwell or Yang-Mills background fields play an important role in the study of laser physics (cf., [5][6][7]), the colour-glass condensate effective theory of heavy ion collisions (cf., [8,9]) and rapidity evolution of Wilson lines modeling high-energy hadron collisions (cf., [10][11][12]). However, the computation of physical observables such as scattering amplitudes in these scenarios is notoriously difficult due to the complicated nature of perturbative Feynman rules defined by the background field Lagrangian. Consequently, there are no tree-level calculations beyond four external states for QFTs in even the simplest strong gauge field backgrounds 1 .
By contrast, recent decades have seen the development of myriad all-multiplicity formulae for the scattering amplitudes of gauge theory in trivial backgrounds at tree-level and beyond. Crucially, these all-multiplicity results are not derived from standard spacetime perturbation theory. This raises the question: is it possible to derive all-multiplicity results in strong field Yang-Mills theory by exploiting alternatives to background perturbation theory on space-time?
In four space-time dimensions, one powerful example of such an alternative is twistor theory; an algebro-geometric formalism that (very roughly) encodes physical fields in geometric data on a complex projective variety called twistor space. Twistor theory enables gluon scattering to be formulated as a perturbative expansion around the integrable selfdual sector of gauge theory. This can be operationalized in terms of a 'twistor string' governing holomorphic maps from the Riemann sphere to twistor space, with the degree of the map corresponding to the helicity configuration of the external gluons in a scattering process (i.e., how 'far away' from self-duality the configuration is) [19,20]. As such, twistor theory underpins one of the most well-known (and earliest) examples of an all-multiplicity formula for scattering amplitudes: the Roiban-Spradlin-Volovich-Witten (RSVW) formula for the full tree-level S-matrix of Yang-Mills theory (in a trivial background) [21].
A priori, it might seem unlikely that a mathematical formalism like twistor theory could say anything about strong field gauge theory, a topic usually studied with the bruteforce methods of background-dressed Feynman rules on space-time. However, the underlying chirality of twistor theory means that it can be used to essentially trivialize the self-dual sector of Yang-Mills theory, which means that it should be applicable to scattering on self-dual backgrounds in much the same way as a trivial background.
In this paper, we present formulae for all tree-level gluon scattering amplitudes in any self-dual (SD) radiative Yang-Mills background. These backgrounds are chiral gauge fields which are characterized by their free radiative data and are asymptotically flat (or 'almost everywhere' asymptotically flat); as such, they admit a well-defined S-matrix. An important special case of such SD radiative backgrounds are the self-dual plane waves, which can be viewed as a coherent superposition of positive helicity gluons valued in a Cartan subalgebra of the gauge group.
Despite having unconstrained functional freedom, these SD radiative backgroundsand gluon perturbations propagating on them -admit an elegant description in terms of twistor theory. For the maximal helicity violating (MHV) sector (i.e., with two negative helicity external gluons and the remainder positive helicity), this allows us to derive -from first principles -momentum space expressions for the corresponding tree amplitudes at arbitrary multiplicity; these include as a special case the formula of [22] for SD plane wave backgrounds. We also provide formulae for all-multiplicity tree-level gluon scattering in any helicity configuration on a SD radiative background which are motivated by twistor string theory. While the formulae for these N k MHV amplitudes are conjectural for k ≥ 2, they are easily seen to pass several fundamental consistency tests.
Although the non-trivial background obviously introduces new features, our formulae share a surprising amount of structure with their well-known antecedents on a trivial background, and manifest remarkable simplicity in comparison to the expectations of space-time perturbation theory. For instance, in the special case where the SD radiative background is Cartan-valued, the MHV colour-ordered partial amplitude is given by: r s 4 1 2 2 3 · · · n − 1 n n 1 d 4 x exp n i=1 (i k i · x + e i g(x, κ i )) , (1.1) where the external gluons are parametrized by their on-shell asymptotic 4-momenta k αα i = κ α iκα i and charges e i with respect to the background, and gluons r, s are negative helicity with all others positive helicity. In this expression, the ratio of angle brackets (constructed from the un-dotted momentum spinors, i j = κ α i κ j α ) is familiar from the Parke-Taylor formula [23] for MHV gluon scattering in a trivial background. But 4-dimensional momentum conservation is absent (as Poincaré invariance is broken by the background), replaced by a single space-time integral against an exponential factor, where the Cartan-valued function g(x, κ) is determined by the SD background.
For a N d−1 MHV amplitude on a generic SD radiative background, our central conjecture is the formula (presented with N = 4 supersymmetry for simplicity) where the integral is over the moduli of a degree d holomorphic map from the Riemann sphere to twistor space, the a i are twistor wavefunctions for the gluon multiplets, and H i are certain holomorphic frames which encode the information about the background. In a trivial background, H i = id, and this expression immediately reduces to the RSVW formula. In general, there are 4d residual (bosonic) integrals in this formula, but when the external gluons are represented by momentum eigenstates, 2(d + 1) integrations can be done explicitly against delta functions (see (4.15)), yielding a formula supported on 2(n − d − 1) of the refined/polarized scattering equations known from ambitwistor strings in trivial backgrounds [24]. In the special case of a self-dual plane wave background, the number of residual integrals can be reduced to 2d−1 (see (4.29)) due to enhanced symmetry. Such residual integrals are a generic feature of amplitudes in strong fields, as the background generically breaks Poincaré invariance and introduces functional degrees of freedom. However, even in the most general case, the 4d residual moduli integrals of our formula are less than the number of integrations expected from space-time perturbation theory in every helicity configuration for n > 3. This number of integrations is precisely what would be generated for a space-time version of the MHV formalism [25] in a nontrivial background, and in section 5 we give the full off-shell propagator in MHV axial gauge in a self-dual plane wave background. Combined with our MHV amplitudes, this yields an MHV formalism on such backgrounds.
The paper is organized as follows. Section 2 begins with an overview of SD radiative gauge fields and their twistor theory. Section 3 derives the MHV gluon amplitude on a SD radiative background from a first-principles using a space-time generating functional. The resulting formula is shown to pass basic consistency tests arising from the weak-field background limit. In section 4, we present a formula for the complete gluon tree-level S-matrix on a SD radiative gauge field background, motivated by twistor string theory, which also passes these basic consistency tests.
Throughout, the special cases of Cartan-valued backgrounds and self-dual plane waves are used as illustrative and explicit examples. Appendix A includes Feynman diagram computations in space-time perturbation theory which match our formulae at 3-and 4points, while appendix B evaluates the formula for the tree-level S-matrix in the special case of an impulsive self-dual plane wave background, where interesting simplifications emerge. Lastly, appendix C develops MHV rules on SD plane wave backgrounds which provide an explanation for the number of residual integrals in our formulae.
2 Twistor theory of self-dual radiative gauge fields A (non-linear) gauge field in Minkowski space is defined to be purely radiative if it is sourcefree, asymptotically flat at both past and future conformal null infinity I − ∪ I + = I and completely determined by the free characteristic data at either I − or I + . In four dimensions, this radiative data is encoded in two functions on I valued in the adjoint representation of the gauge group. The natural observable of Yang-Mills theory in such a configuration is the S-matrix, defined perturbatively on the radiative background by the scattering map of data from I − to that at I + . The reader uninterested in precise definitions of radiative self-dual gauge fields can skip section 2.1, proceeding to section 2.2 which developes the twistor theory of self-dual radiative backgrounds.
Complexifying the data and background field, one of the free radiative functions on I can be set to zero to obtain a self-dual (SD) solution of the Yang-Mills equations. Such solutions can only be real in split or Euclidean signature, but we will not impose such reality conditions; complex circularly polarized solutions have meaning in the quantum theory. In four dimensions, the SD sector is famously integrable: SD solutions can be constructed from complex analytic data with twistor theory [26]. For SD radiative gauge fields, there is a twistor construction which manifests their integrability in terms of the characteristic data at I [27]. In this section, we review this twistorial description of SD radiative gauge fields and demonstrate how gluon perturbations on these backgrounds are encoded in their twistor theory. Throughout, we highlight the concrete example of self-dual plane waves (SDPWs); while not strictly asymptotically flat, as the fields do not fall off fast enough in the null symmetry direction, they are sufficiently 'almost asymptotically flat' for scattering to be well-defined [28,29].
In four-dimensions, a 2-form F decomposes into self-dual (SD) and anti-self-dual (ASD) parts as: where * is the Hodge star of the four-dimensional Minkowski metric. In the 2-spinor formalism, the decomposition of a field strength into its SD and ASD parts is given by with the SD field strength encoded byFαβ =F (αβ) and ASD field strength encoded by F αβ = F (αβ) . Thus, the self-duality condition F − = 0 is F αβ = 0 in the 2-spinor formalism.
The conformal null infinity of Minkowski space, I , is obtained in the usual fashion by conformal compactification. To preserve Lorentz invariance, we work with retarded coordinates (u, r, λ α ,λα), where λ α = (λ 0 , λ 1 ) are holomorphic homogeneous coordinates on P 1 , subject to the equivalence relation (u, r, λ α ,λα) ∼ (|b| 2 u, |b| −2 r, b λ α ,bλα) , ∀b ∈ C * . (2.5) The standard non-homogeneous retarded coordinates (u, r, ζ,ζ) are recovered by choosing a constant future-pointing time-like vector t αα = diag (1/ √ 2, 1/ √ 2) obeying t 2 = 1, and definingλ α := √ 2 t ααλα . Rescaling so that λλ = √ 2 gives the original (u, r) and With such choices, x αα = u t αα + r λ αλα , so that λ αλα is the null vector aligned along the outward going light cones from the origin. Conformal compactification is achieved through the inversion R = r −1 , with future null infinity I + brought to the finite hypersurface R = 0 by a conformal rescaling given by (2.6) The natural conformal equivalence class of degenerate metrics on I + is then written in homogeneous coordinates as with (u, λ α ,λα) the projective Bondi coordinates on I + ∼ = R × S 2 , (λ,λ) being homogeneous coordinates on S 2 and u the coordinate on R [31][32][33]. This projective formulation [27,[34][35][36][37] gives rise to the line bundles O(p, q) → I + whose sections are represented by homogeneous functions f (u, λ,λ) obeying f (|b| 2 u, bλ,bλ) = b pbq f (u, λ,λ) for any b ∈ C * . The metric (2.7) takes values in O(2, 2), so to fix a conformal scale -for instance, to give the round metric on S 2 -one must divide by λλ 2 . More generally, functions with non-trivial spin and conformal weight (s, w) on I + in the standard (non-projective) description correspond to sections of O(p, q) with s = p−q 2 and w = p+q 2 . Translations by a αα on space-time act at I + via u → u + a αα λ αλα . Thus I + is the total space of an affine bundle O R (1, 1) → P 1 , whose sections are real-valued functions f (λ,λ) which obey f (bλ,bλ) = |b| 2 f (λ,λ). In curved space-time, the choice of zero section is acted on more generally by supertranslations u → u+f where f is any section of O R (1, 1); thus, there is no preferred choice of zero-section for the u-coordinate.
Points x ∈ M are realized at I + by their light cone cuts: the intersection of the light cone from x with I + . To describe this cut, parametrize the light cone of x by x αα + r λ αλα with λλ = √ 2 and identify it with u t αα + r ′ λ αλα for some u, r ′ . Contracting with λ αλα gives u = x αα λ αλα , (2.8) as the light cone cut of x at I + . Although our focus here and elsewhere will be I + , identical considerations apply to I − . Indeed, the projective description of I − has fibre coordinate on the total space of O R (1, 1) → P 1 given by the advanced time coordinate v = t + r rather than u.
Radiative gauge fields: A linear, negative helicity spin-1 field is described by its linearised ASD field strength φ αβ (x), which obeys the zero-rest-mass equation ∂ αα φ αβ = 0. This field 'peels' as r → ∞ according to [33,38] where λ α is aligned and covariantly constant along the light rays of constant (u, ζ,ζ) (i.e., with tangent λ αλα ). The object φ 2 = φ 2 (u, λ,λ) is the radiation data, defining a section of O(−3, −1) on I in the homogeneous formalism. Assuming suitable globality, it is straightforward to see that linear spin-1 fields are determined by their radiation data at I by means of the Kirchoff-d'Adhémar integral formulae adapted to I in the unphysical (i.e., conformally rescaled) space-time. This is an integral over the light-cone cut [30,38,39] It is easy to see by differentiating under the integral sign that this solves the zero-restmass equation, and corresponds to the given radiation data when the field is defined and differentiable over I + including the vertex i + (i.e., future time-like infinity) [30]. The positive helicity spin-1 fieldφαβ is given by complex conjugation. These assumptions do not apply in the presence of sources (e.g., when the field has Coulombic parts), and these are therefore not encoded in φ 2 . This motivates our definition of radiative gauge fields: Definition 2.1 A radiative gauge field in four-dimensional Minkowski space-time is one that is completely characterized by its radiation data, so that φ 2 is the free characteristic data for the field.
This definition also applies to non-linear gauge fields, where the gauge potential is a more natural starting point than the radiation field itself.
Self-dual radiative gauge fields: Let A be a gauge field in four-dimensional Minkowski space-time, viewed as a 1-form valued in the Lie algebra g of a compact gauge group. In retarded Bondi coordinates, an asymptotically flat gauge field in temporal gauge A u = 0 can be restricted to I + and pulled back to projective coordinates in the conformally compactified space-time as [40][41][42][43]: Thus, A 0 and A 0 have spin weights −1 and +1, respectively, and both have conformal weight −1.
For the gauge field (2.10), the 2-form field strength F at I + has leading (in the conformally rescaled space-time) SD and ASD parts (cf., [41,44,45]) so A 0 ,Ā 0 respectively determine the leading asymptotic ASD and SD field strengths, sometimes known as the broadcasting function of the gauge field 2 , to be with values in O(−3, −1)⊗ g. The broadcasting function, and thus A 0 , actually determines the free radiative data for the fully non-linear gauge field [40,47,48] by (2.9) with suitable regularity assumptions. However, we also wish to include plane waves in the class of radiative fields, as they are essentially pure radiation. Thus, we allow for A 0 to have isolated singular points on the S 2 at I + (but generally extended in u); such gauge fields will be 'almost everywhere' asymptotically flat. For Lorentz-real gauge fields, there are no purely self-dual (SD) or chiral radiation fields since A 0 andĀ 0 are complex conjugates. However, complex gauge fields have independent data A 0 andÃ 0 , so it is possible to define a purely SD radiative gauge field: Definition 2.2 A self-dual radiative gauge field is a complex radiative gauge field that is asymptotically flat almost everywhere and determined by its radiative dataÃ 0 , with A 0 = 0. Example: Self-dual gauge theory plane waves An important example of a SD radiative gauge field -to which we repeatedly return -is that of a self-dual plane wave (SDPW). This is a highly-symmetric example, being invariant under translations orthogonal to a preferred null direction n αα = ι αια , and is furthermore valued in a Cartan subalgebra h ⊂ g [29,49]. In lightfront coordinates the metric on M becomes (2.14) With n · x = x − the SDPW gauge field may be written where f (x − ) is an arbitrary function of the lightfront coordinate x − , taking values in the Cartan subalgebra h. Choosing a spin frame so that o αια dx αα = dz and n αα = ι αια (i.e., The SDPW gauge field (2.15) is clearly invariant under translations of x + , z andz but it also has two other (less obvious) null rotation symmetries that together form a five-dimensional Heisenberg algebra whose centre is n (cf., [29,50]). Wheṅ f and hence the field strength is compactly supported in x − , the SDPW is an example of a sandwich plane wave, which admits a well-defined S-matrix [28,29]; we assume this sandwich property from now on. These solutions do not satisfy standard asymptotically flat fall-off conditions as r → ∞. On one hand, the field is compactly supported along every light ray tangent to λ αλα when λ ι = 0. On the other hand, it is constant and has no fall-off along light rays with λ α = ι α . We can nevertheless make sense of the radiation data as a distribution supported at λ ι = 0. In the projective coordinates for I + the free radiative data is given bỹ where the Cartan-valued function 3 F is the first antiderivative of f (Ḟ = f ), and the holomorphic delta function (technically, a (0, 1)-distribution) is defined as for any complex variable z. Thus,Ã 0 has a double delta function singularity at the generator of I + where λ ι = 0 and so is only almost everywhere asymptotically flat. The SDPW broadcasting function is therefore 4 These identifications can be verified by means of (2.9).

Twistor theory of SD radiative gauge fields
The twistor space PT of complexified Minkowski space M ∼ = C 4 is an open subset of P 3 . Using homogeneous coordinates Z A = (µα, λ α ), define 5 Points x ∈ M correspond to holomorphic, linearly embedded Riemann spheres X ∼ = P 1 in PT [52] given by the incidence relation For such finite x αα , the twistor variables λ α are identified with the homogeneous coordinate λ α on each twistor 'line' X ∼ = P 1 . It is easy to see that two twistor lines X, Y ⊂ PT intersect if and only if the corresponding points x, y ∈ M are null separated. Conversely, holding Z A constant, (2.20) corresponds to a totally null ASD two-plane in complexfied Minkowski space M. Selecting a real slice of M with Euclidean R 4 or split R 2,2 signature gives rise to a complex conjugation on PT that is quaternionic for Euclidean signature [53,54], or that fixes the RP 3 on which all components of Z A are real for split signature. Our focus in this paper will be on Lorentz signature R 1,3 where the reality structure induces a pseudo-Hermitian structure of signature (2, 2) on PT given by Points of real Minkowski space correspond to lines X that lie in PN = {Z ∈ PT|Z ·Z = 0}, which conversely is the space of real light rays in R 1,3 .
Since each light ray in R 1,3 makes it out to I + , there is a natural projection given by where the reality of u follows from Z ·Z = 0. It follows from the formula in (2.22) that p naturally extends to p : where I + C is a partial complexification of I + obtained by taking u ∈ C but retaining λα =λα inside the full complexification. Thus p is holomorphic in µα but of course not in λ α . The map p projects X given by (2.20) to the light cone cut of x in I + (2.8). 4 That the curvature componentφ2 at I + should be identified with the gauge dependent f is perhaps a puzzle. This is resolved by the fact that the radiative assumption gives the preferred choice for f that arises by integrating the gauge invariantḟ from i + so that i + is regular. 5 We follow the conventions of [51].
Characteristic data on I + C can be pulled back using p to PT (or PN if it has no analytic extension in u). On PT the Kirchoff-d'Adhémar integral formula (2.9) can be re-interpreted [55] as a gauge fixed version of the Penrose transform [56,57]. The Penrose transform is the correspondence In the Dolbeault representation of the cohomology group, an element f ∈ H 1 (PT, O(−4)) is a∂-closed (0, 1)-form on twistor space modulo exact forms of homogeneity weight −4. The correspondence operates via the integral formula where f | X denotes restriction to the twistor line X ∼ = P 1 via the incidence relations (2.20). The Kirchoff-d'Adhémar integral formula (2.9) follows as a special gauge fixing within the Dolbeault equivalence class via the identification which is easily seen to be homogeneous of weight −4 in Z A and∂-closed. For positive helicity, the Penrose transform reads By setting f =Ã 0 Dλ, this can be identified with the conjugated version of (2.9).
Twistor description of SD radiative gauge fields: The positive helicity Penrose transform has nonlinear extensions for SD gauge fields [58], with wide application to the construction of monopoles and instantons, and indeed to the general theory of integrable systems (cf., [26,59,60]). For the SD Yang-Mills equations, the relevant construction is: There is a one-to-one correspondence between: • SD gauge fields on M, and • holomorphic vector bundles E → PT with E| X topologically trivial for every twistor line X ∼ = P 1 corresponding to x ∈ M.
Further restrictions on E → PT give details of the gauge group on space-time and reality structure; for example, for gauge group SU(N ) in split signature, E has rank N , is equipped with a positive definite real form on the real slice of PT, and det E is trivial [61]. The bundle E → PT is often referred to as the Ward bundle.
Theorem 1 can also be understood in terms of radiative data at I + via a Dolbeault presentation. This expresses the complex structure on the Ward bundle E → PT in terms of a∂-operatorD taking sections of E to (0, 1)-forms that satisfiesD 2 = 0. SD Yang-Mills is conformally invariant and the bundle and connection naturally extend to I + as a bundle with the connection given in (2.10). The Ward bundle for a SD radiative gauge field then arises simply by pulling back this bundle with∂-operator arising from (2.10) to PT [27,41,45] with the (0, 1)-form a valued in End E. ThatD 2 = 0 follows trivially, as the only antiholomorphic dependence of a is inλ, and Dλ ∧ Dλ = 0.
Proof of Theorem 1: In the non-abelian case, we can work on I + following [27,41,45]. Upon restriction to any line X ⊂ PT, it follows from the assumptions 6 in Theorem 1 that E is holomorphically trivial, which is the statement that we can find a frame H(x, λ,λ) : This relation is sometimes referred to as the 'Sparling equation,' and by Liouville's theorem any solution will be unique up to H → Hg(x) for some matrix function g independent of λ. Given such an H, it follows that the integrability of which implies for F αβ the ASD field strength of the gauge connection on M. Thus, integrability of the Lax pair is equivalent to the SD Yang-Mills equations F αβ = 0 for the connection D αα . ✷ In general, it is hard to solve for H explicitly. However, whenÃ 0 takes values in a Cartan subalgebra h ⊂ g, theD in (2.27) defines a cohomology class a ∈ H 0,1 ∂ (PT, O ⊗ h), and the SD radiative gauge field on M can be recovered using (2.26). In this case, a| X ∈ H 0,1 ∂ (P 1 , O ⊗ h); but this cohomology group vanishes, so for some Cartan-valued g(x, λ) of homogeneity zero in λ. This can be obtained explicitly from a Green's function on the sphere as where the choice of the spinor o α amounts to an arbitrary gauge choice. With this we see thatDH| X = 0 is solved by setting where exp(· · · ) denotes the exponential map of the gauge group. Now (2.28) gives whereφ 2 = ∂ uÃ 0 is the SD broadcasting function. From this one obtains the integral formula for A itself: expressed equivalently in terms of the complex structure on the Ward bundle or the radiative data.
Background coupled gluons: A linearised gluon field a αα (x) propagating in a background gauge field A αα has a linearised field strength f ab = D [a a b] which can be decomposed into its SD and ASD parts:fαβ and f αβ , respectively. However, on a non-trivial background there is no longer a natural global decomposition a b = a + b + a − b of the potential into SD and ASD parts, as an infinitesimal gauge transformation a a = D a γ leads to f αβ = [F αβ , γ] where F is the curvature of the background field, so such a decomposition cannot be gauge invariant. When the background gauge field is SD, it is gauge invariant to characterise a positive helicity gluon by f αβ = 0. Although it is not possible to write down a gauge invariant potential for a negative helicity gluon, it can nevertheless be characterised asymmetrically by a curvature 2-form f αβ satisfying D αα f αβ = 0.
For a positive helicity gluon (with suitable analyticity) on a SD radiative gauge field background, the Penrose transform represents the gluon in terms of cohomology in twistor space [56,57,59]: where H 0,1 D denotes the Dolbeault cohomology with respect toD. To see how this isomorphism works in the direction from twistor space to space-time, let a ∈ H 0,1 D (PT, O⊗End E).
Upon restriction to the twistor line X, E| X is holomorphically trivial with trivialization H, so a| X takes values in a trivial bundle on X. In this holomorphic frame a| X becomes , and the triviality of this cohomology group implies that the Green's function used in (2.32) can be used to find some j(x, λ) such that for j(x, λ) valued in the adjoint and of homogeneity zero in λ.
Acting on both sides of (2.37) with λ α D αα and using (2.29) gives∂| X (λ α D αα j) = 0. By Liouville's theorem, this implies the existence of an adjoint-valued gluon field a αα (x) such that Contracting both sides of (2.38) with λ β Dα β (and recalling that the background is SD and j is essentially a perturbation to H) shows that this is a positive helicity gluon (i.e., f αβ = 0), as desired.
For a negative helicity gluon in a SD radiative background gauge field, the Penrose transform is This integral produces a well-defined adjoint-valued object on M as H −1 b| X H takes values in H 0,1 (P 1 , O(−4) ⊗ g). Using (2.29), it follows that as required.
Momentum eigenstates on SD backgrounds: For the computation of scattering amplitudes, it is useful to have cohomological representatives corresponding to momentum eigenstates. That is, what are the twistor representatives corresponding to gluons which have null momentum k αα = κ ακα at I − C (i.e., before passing through the non-trivial SD radiative background)? Remarkably, the representation (2.27) for the partial connection D ensures that standard twistor representatives for momentum eigenstates (in a trivial background) provide the answer.
For a positive helicity gluon, consider the representative where T a is a generator of the adjoint representation of the gauge group and the holomorphic delta functionδ has support only where both components of its argument vanish, i.e., where λ α = sκ α . The second equality in (2.42) follows on the support of the delta function only whenκ =κ, but then serves to give a(Z) in terms of its distributional characteristic data for a momentum eigenstate. It is straightforward to show that a isD-closed as a consequence of the fact that a is proportional to Dλ and only depends on anti-holomorphic variables throughλ.
To obtain the corresponding gluon wavefunction, one follows (2.37) -(2.38); on a general non-abelian SD radiative background, it is difficult to solve for H explicitly, so let us make the simplifying assumption that the background is Cartan-valued. In order to keep track of the colour structure, assume that T a is an eigenstate of the Cartan h so that the colour structure of the representative is not deformed when conjugating by the holomorphic frame. Thus T a defines a 'charge' e i with respect to the background, defined as the root e = (e i ) 1≤i≤rank(g) , with e i defined by Upon performing the scale integral in (2.42) and applying (2.33), this means that where ξ α is an arbitrary constant spinor (on which the expression does not depend). From (2.43) and (2.37), it follows that is easily seen to solve the Hamilton-Jacobi equation To determine the positive helicity gluon field itself, it helps to rewrite the relation (2.34) as for gα(x, λ) two Cartan valued functions which are homogeneous of degree −1 in λ. Using the Schouten identity, these can be explicitly found from (2.32) and (2.35) to be where the spinorKα (x) :=κα + e gα(x, κ) , (2.49) defines the on-shell null momentum of the gluon field K αα (x) := ∂ αα φ − e A αα = κ αKα (x) at any point; self-duality of the background means that only the dotted momentum spinor is 'dressed' by the background.
For a negative helicity gluon with momentum κ ακα at I − C , take the representative , and feeding this into the integral formula (2.40) gives As expected by the self-duality of the background, this differs from the linearised field strength of a negative helicity gluon in a trivial background only by the appearance of φ(x) in the exponential.
Example: SDPWs For the special case of the self-dual plane wave (SDPW) given by (2.15), the partial connection on the associated Ward bundle E → PT is given by Plugging this into (2.35), we immediately recover the SDPW gauge potential A αα (x) = −f (x − ) o αια , and from (2.46), upon using the Schouten identity or equivalently (2.47).
In lightfront coordinates (2.13), the components of an on-shell 4-momentum are and combined with (2.53) this yields the associated solution to the Hamilton-Jacobi equation (2.45): As expected, this is the chiral truncation of the solution to the Hamilton-Jacobi equation in a general plane wave background [7,29,65]. The dressed on-shell momentum at any lightfront time also matches the chiral projection of a dressed gluon momentum in a generic plane wave background to the SD sector [66,67].
Making use of the null symmetry of the SDPW background, the lightfront gauge condition n · a = ι αια a αα = 0 can be imposed on both positive and negative helicity gluon fields. In this positive helicity case of (2.48), this simply fixes ξ α = ι α , while in the negative helicity case it allows explicit determination of the gluon wavefunction: This transverse-lightfront gauge is the only gauge for which the negative helicity gluon has no SD linearised field strength.

MHV scattering on gauge theory backgrounds
We now turn to tree-level MHV scattering of gluons in a SD radiative gauge field background. The space-time generating functional for these amplitudes is lifted to twistor space and perturbatively expanded to yield an all-multiplicity expression for MHV scattering which displays much of the simplicity of its flat background cousin, the Parke-Taylor formula [23].

The MHV generating functional
The chiral formulation of gauge theory in (complexified) Minkowski space M is given by the Chalmers-Siegel action [68] S where F αβ is the ASD part of the field strength of the gauge field A, B αβ is an auxiliary ASD 2-form valued in the adjoint of the gauge group, and g is the dimensionless coupling constant. The field equations are so B αβ can be integrated out, resulting in a Lagrangian which differs from the standard Yang-Mills Lagrangian by a multiple of the topological term tr(F ∧F ). Thus, the Chalmers-Siegel action is perturbatively equivalent to Yang-Mills theory in Minkowski space, and will produce the same gluon amplitudes. Furthermore, it is apparent from (3.2) that in this formulation, g = 0 corresponds to the SD sector (i.e., F αβ = 0). Thus, small g ≪ 1 perturbation theory corresponds to a perturbative expansion around the integrable SD sector, with the auxiliary field B αβ encoding ASD fluctuations on a SD background. This formalism is particularly useful for giving a geometric perspective on MHV scattering amplitudes. The MHV helicity configuration (on any perturbative background) involves two negative helicity gluons and an arbitrary number of positive helicity gluons. When the perturbative background is trivial or a SD radiative gauge field, this configuration can be realized as a two-point function of the negative helicity gluons on a SD background [69]. Expanding this composite background into a non-linear SD radiative gauge field and a superposition of positive helicity external gluons then recovers the desired MHV amplitude.
The chiral Chalmers-Siegel formulation is well-adapted to this setup. Let R denote the infinite-dimensional phase space of fields (A, B) ∈ Ω 1 (g) ⊕ Ω 2 − (g) which solve the field equations (3.2), and P the Yang-Mills phase space obtained as the quotient of R by gauge transformations. Let (A, 0) ∈ R be a SD solution that can be perturbed to define a superposition of a SD radiative gauge field and an arbitrary number of positive helicity gluons on that SD background. The boundary term of (3.1) defines a symplectic form on P: where C ⊂ M is a three-dimensional Cauchy surface and δ is the exterior derivative on R. It is standard that ω is independent of the choice of C and descends to become nondegenerate on P after quotienting by gauge transformations.
where β γδ = β (γδ) are the components of the linearised ASD 2-form: β = β γδ dx γα ∧ dx δα . The linearised equations induce a splitting of U by a short exact sequence where U + is the space of linearised SD solutions modulo gauge: As mentioned in §2.2, the SD field strength of the background means that the space of linearised ASD solutions must be defined asymmetrically, so U − is given by fluctuations in the ASD curvature, but for which there is no canonical SD-free potential. We can, however, define what it means to be purely ASD at a hypersurface using the symplectic form as follows. For It is clear that for A 1 , A 2 ∈ U + we have ω(A 2 , A 1 ) = 0 so that U + is a Lagrangian subspace of U . This also gives duality pairing between (α 1 , β 1 = 0) ∈ U + and β 2 ∈ U − by This is not hypersurface independent when β 1 = 0. Nevertheless, it can be used to define a splitting of (3.5), U = U + ⊕ U − , at any given C following [69] by: It is easily checked that at I + So assuming falloff as u → ∞, the definition for A 1 to be ASD at I + is equivalent to settingĀ 0 1 = 0, as expected. Thus, one can prepare a positive frequency ASD field at I − and measure how much SD part it has acquired from the SD background by the time it gets to I + . This can be measured by pairing with a negative frequency ASD field A 2 at I + : this pairing will be the generating functional for the MHV amplitude. So let A 1 be ASD at I − and A 2 be ASD at I + ; evaluating the pairing at I + and using definition 3.1 gives the generating function for the MHV amplitude as upon using Stokes' theorem for ∂M = I + − I − , the fact that A 1 is ASD at I − and the linearised field equations (3.4). Note that in the final line, we have normalised β 1,2 → g −2 β 1,2 for later convenience.
This ω(A 2 , A 1 ) is the generating functional for MHV amplitudes, which we wish to expand perturbatively in terms of positive helicity gluons on a fixed SD radiative gauge field background. The right hand side depends implicitly on the background A as the β 1,2 satisfy background-coupled field equations. The MHV amplitude will be obtained by iteratively perturbing this dependence on A. This would be challenging on space-time, but we will see that on lifting (3.11) to twistor space, where integrability of the SD sector is manifest, it becomes straightforward, leading to a simple extension of the standard MHV amplitude.
Lift to twistor space: The lift of (3.11) to twistor space is easily achieved using the tools in Section 2.2. Observing that the generating functional is equivalent to and partial connection associated to the SD background, respectively. Consequently, the generating functional is lifted to twistor space using (2.40): where H is the global holomorphic frame for E| X . The integral is over two copies of the twistor line X 1,2 ∼ = P 1 ⊂ PT, followed by an integration over the moduli space of lines (i.e., M itself). The notation H 1 stands for H(x, λ 1 ), b 1 stands for b| X 1 so that b 1 = b 1 (x, λ 1 ), and so forth. 8

Perturbative expansion and MHV amplitude
We now take the SD background in the generating functional (3.13) to be constructed from a fixed, classical SD radiative background (treated non-perturbatively) and a finite number of positive helicity gluons propagating on this background. On twistor space, this gives thē ∂-operator on E → PT as∂ + a + a :=D + a , (3.15) where a is defined by the characteristic data of the SD radiative background via (2.27), and a is simply a (0, 1)-form that will later be expanded to finite order as a = i a i with each a i ∈ H 0,1 D (PT, O ⊗ End E) encoding SD perturbations/positive helicity gluons on the background.
In order to expand the generating functional (3.13) in a around the backgroundD, note that dependence on a arises only through the holomorphic frame H for the∂-operator The kernel of (D + a) −1 acting on sections of O(−1) ⊗ E| X is given by 9 The generating functional (3.13) also follows from a Lagrangian perspective by lifting the fully off-shell Chalmers-Siegel action (3.1) to twistor space [70] as: and D 3 Z stands for the canonical section of Ω 3,0 (PT, O (4)). This twistor action is equivalent to that on space-time; solutions to its equations of motion modulo gauge are in one-toone correspondence with solutions to (3.2), and the actions agree when evaluated on these corresponding extrema [70]. The two-point amplitude of two negative helicity gluon fields evaluated on a SD background is given by the bi-linear portion of the action evaluated on the (on-shell) external states. When b1, b2 ∈ H 0,1 D (PT, O(−4) ⊗ End E) are the two ASD external states withD on E describing the background, the bi-linear part of the action is then the generating functional (3.13). Thus the MHV amplitude follows from background perturbation theory in twistor space itself (cf., [71][72][73]). 9 This is the 'holomorphic Wilson line' forD + a along the complex line X of [63,64,74].
For H satisfyingD| X H = 0, write We can now expand U X using the geometric series about a = 0 of (D + a) −1 , which gives with the identification λ 0 ≡ λ.
In terms of U X , the MHV generating functional is: Here, an overall numerical factor has been dropped, the label n + 1 ≡ 1, and the external wavefunctions have been relabeled after performing the expansion so that the negative helicity states are in an arbitrary position in the colour trace.
To obtain the momentum space expression for the MHV amplitude, we insert momentum eigenstate representatives (2.42) and (2.50) into (3.21). As usual, it is possible to be more explicit when the SD radiative background is valued in a Cartan subalgebra; in this case, let each external gluon come with colour vector T i in some root eigenspace of the Cartan so that [a, T i ] = e i a for some root valued charges e i . The P 1 integrals in (3.21) are performed against the holomorphic delta functions in the momentum eigenstate representatives, setting λ i to the corresponding κ i momentum spinor. Further, it follows that H −1 i a i H i = e e i g(x,κ i ) a i , where g(x, κ i ) is the Cartan-valued function defined by the background (2.31) -(2.33).
Thus, the resulting formula for the colour-ordered MHV amplitude in a Cartan-valued SD radiative background is: where charge conservation must hold in order for the total trace to be non-zero, as the trace is invariant under the action of the Cartan subalgebra. As expected, the SD background does not affect the ratio of ASD (angle bracket) momentum spinor contractions familiar from the Parke-Taylor formula on a flat background. However, 4-momentum conservation is replaced by an integral over Minkowski space which cannot be performed explicitly, since the full functional freedom of the background enters through the function g appearing in the exponential. The replacement of momentum conservation by such residual space-time integrals is a familiar feature of background perturbation theory; a generic strong background breaks Poincaré invariance. However, the expectation from standard background perturbation theory is that an n-point tree-level gluon amplitude should have n − 2 residual space-time integrals in a strong background. It is therefore remarkable that (3.22) contains only a single space-time integral regardless of the number of positive helicity gluons in the MHV sector; in this framework, the amplitude behaves like a single contact term.
While surprising from the perspective of the standard Yang-Mills Lagrangian, this is not unexpected from the perspective of the MHV formalism, whereby the Feynman rules of gauge theory are entirely composed of MHV vertices and scalar propagators linking positive and negative helicity legs [25]. The MHV formalism can be obtained from a particular gauge-fixing of the twistor action (3.14) [75][76][77], which can be implemented non-perturbatively. The MHV generating functional is equal to the non-local term in this twistor action, which contains only a single space-time integral. Another, less general, perspective is provided by a non-local field redefinition which can be used on space-time to recast the perturbative Yang-Mills Lagrangian in a form where MHV interactions are captured by a single vertex at arbitrary multiplicity [78,79], although this is only possible for a subset of SD radiative backgrounds which have commuting null and translational symmetries (such as SDPWs).
Evaluation on SDPWs: Substituting (2.53) into (3.22) enables further simplification of the MHV amplitude on a SDPW background. In this case, the integrand depends only on the lightfront coordinate x − . This allows us to perform the integrals over three of the coordinates -x + , z,z -to produce a 3-momentum conserving delta function, denoted by δ 3 +,⊥ to signify that the z,z coordinates are transverse to the wave. The remaining lightfront integral is left intact, in keeping with the genericity of the wave-profile f (x − ). The final answer can be written compactly as where F n (x − ) is the Volkov exponent, a common building block of amplitudes in strong electromagnetic fields and laser physics (cf., [7,65]). It is defined by choosing any set of n − 1 distinct external gluons, say {i 1 , . . . , i n−1 } ⊂ {1, . . . , n}, and setting for K ia the background-dressed null momenta (2.57). Then F n is given by with the equivalence of these expressions following from the support of 3-momentum conservation and the fact that K 2 i = 0 for each external gluon. This is precisely the MHV amplitude for gluons on a SDPW gauge field background presented in [22], derived here as a special case of the MHV amplitude on a SD radiative background.

Consistency checks
It is straightforward to show that the MHV amplitude (3.22) passes basic consistency checks. In appendix A, the formulae are explicitly matched against Feynman diagram computations in background perturbation theory at 3-and 4-points. This becomes unwieldy for higher points n ≥ 5 (indeed, even the n = 4 calculation is fairly complicated). Here, we show that the formula has the correct weak background limit.
For the flat background, (3.22) reduces to the standard Parke-Taylor formula for MHV gluon scattering [23]: in a trivial background g(x, κ i ) becomes independent of κ i so the second term in the exponential of (3.22) vanishes due to charge conservation. The integral over M can then be performed against the remaining exponential and gives momentum conservation, which is the Parke-Taylor formula. In the weak field perturbative limit of (3.22), the background is treated to first order and we should recover the Parke-Taylor formula with an additional positive helicity gluon corresponding to the linearised background. Perturbatively, a general SD radiative background is generated by SDPWs, so the perturbative limit corresponds to expanding (3.24) to first order in f (x − ). This enters only through the Volkov exponent (3.26) to give (3.28) To evaluate this perturbative limit, the background is represented by a Fourier mode for q αα := ω ι αια the on-shell momentum of the background, which can now be viewed as a single positive helicity gluon. Note that if the perturbative background is a pure photon, then e j = 0 for all j = 1, . . . , n and (3.28) vanishes, as required by the photon-decoupling identity [80].
Recalling the definition of the colour-charges, it follows that 30) for t i the generator of the Cartan associated with the background and T a j the generator associated to the external gluon. This means that the sum over j = 1, . . . , n appearing in (3.28) is equivalent to summing over the ways of inserting t i into the colour-ordering of the original external gluons: where we have performed the ds integral and the sum over j is understood to be cyclic (i.e., j = n + 1 = 1). The relative minus sign between the terms in parentheses on the first line arises from the commutator in (3.30). The dx − integral in (3.31) can now be performed straightforwardly to give κ iκi r s 4 1 2 2 3 · · · n − 1 n n 1 n j=1 j j + 1 j ι ι j + 1 , (3.32) after using the Schouten identity and dropping an irrelevant overall factor of −2πi ω −1 . The j th term in this sum is which is the Parke-Taylor formula for (n+1)-point MHV scattering in a trivial background, where the additional positive helicity gluon has momentum q αα and is inserted between gluons j and j + 1 in the colour ordering. This is precisely the expected perturbative limit, with the remaining terms in (3.32) giving all other locations of this additional gluon in the colour ordering.

N k MHV amplitudes
We now consider tree-level amplitudes for gluons and gravitons on self-dual backgrounds in a generic N k MHV helicity configuration. In a trivial background, twistor-string theory [19,20,81,82] leads to remarkable formulae for N k MHV amplitudes in gauge theory and gravity as integrals over the moduli of holomorphic maps from the Riemann sphere to twistor space; the degree d of the map is related to the helicity configuration by d = k + 1 [21].
As we have seen, the twistor correspondence works just as well with a SD radiative gauge field, so the same principle can be applied to amplitudes on such backgrounds. This again leads to a well-defined formula that is much simpler than might be expected from space-time perturbation theory. There are 4(k + 1) twistor moduli integrals at N k MHV. Using space-time Feynman diagrams the number of integrals is given by the space-time dimension multiplied by the number of vertices: 4(n − 2). Thus, the twistor moduli space counting is as one might expect for a space-time analogue of an MHV vertex expansion in which there are k + 1 MHV vertices, a significant reduction for general amplitudes when n > 3. Restricting to a SDPW, there are significant further simplifications arising from the additional symmetries.
In this section, we work with N = 4 super-Yang-Mills (SYM). The inclusion of supersymmetry is just a tool to simplify the resulting formulae; since the expressions are at tree-level, individual particle content is recovered by extracting components from an expansion in the supermomenta. Unlike the MHV amplitudes, these N k MHV formulae are conjectural: we do not yet have a systematic way to derive them from first principles. However, we show that they pass several non-trivial consistency checks, such as having the correct flat background and perturbative limits.

N = 4 super-Yang-Mills & the RSVW formula: Supersymmetric twistor space PT
is an open subset of P 3|4 with homogeneous coordinates Z I = (µα, λ α , χ a ); this is the target space of twistor string theory for N = 4 super-Yang-Mills (SYM). On a trivial background, the tree-level S-matrix of N = 4 SYM is computed by worldsheet correlators in twistor string theory governing holomorphic maps from the Riemann sphere P 1 to PT [19,20,81,82]. The resulting Roiban-Spradlin-Volovich-Witten (RSVW) formula [21] for the N = 4 SYM tree-level S-matrix is obtained as an integral over the moduli space of degree d holomorphic maps into twistor space for amplitudes in the N d−1 MHV helicity sector.
Such a degree d map can be represented as where σ a = (σ 0 , σ 1 ) are homogeneous coordinates on the P 1 worldsheet, the parameters U I a(d) are the map moduli and a(d) is a multi-index denoting a totally symmetric object with d indices. There are 4(d + 1) bosonic and 4(d + 1) fermionic degrees of freedom in the moduli. The RSVW formula for the colour-ordered tree-level amplitudes at NMHV degree d − 1 is then [19,21] where (i j) = ǫ ab σ i b σ j a is the SL(2, C)-invariant inner product on the homogeneous coordinates of P 1 and Dσ = (σ dσ). Division by the (infinite) volume of GL(2, C) is understood in the Faddeev-Popov sense, acting on all integrated parameters, and accounts for the C * × SL(2, C) degeneracy in the description of the holomorphic map (4.1). The a i ∈ H 0,1 (PT, O) ⊗ g are twistor wavefunctions for the external gluon multiplets, usually taken as momentum eigenstates (cf., [21,74,83]). The veracity of this formula is established by showing that it has the correct factorization properties [84][85][86][87]. An on-shell N = 4 SYM multiplet is characterised in a momentum eigenstate representation by an on-shell supermomentum (κ ακα , κ α η a ), with corresponding twistor representatives (cf., [74]): With such momentum eigenstates inserted in (4.2), the moduli of the µα and χ a components of the map (4.1) can be integrated out to obtain delta functions [21]: so isolating the appropriate coefficient is simply an algebraic exercise. For example, in the case of the MHV amplitude (d = 1), the required coefficient is Amplitudes on a SD radiative background: The formula (4.2) naturally extends to any SD background, described on PT by a holomorphic vector bundle E → PT with partial connectionD =∂ + a obeyingD 2 = 0. On restriction to a rational curve Z(σ), this bundle will (generically) be holomorphically trivial, with trivialization H(U, σ) defined byD| Z(σ) H = 0. Equipped with this holomorphic trivialization, factors of (i i + 1) −1 appearing in (4.2) due to Wick contractions in the twistor string worldsheet current algebra are replaced by H −1 i H i+1 (i i + 1) −1 , where H i := H(U r , σ i ). This leads to a formula for the colour-ordered N d−1 MHV tree amplitudes on the background: with the external gluon wavefunctions a i ∈ H 0,1 D (PT, O ⊗ End E). For a general, nonabelian background, it is difficult to determine the holomorphic frame H(U, σ) explicitly, but for a SD radiative background valued in a Cartan h ⊂ g, H can be written as an integral formula.
For such a SD radiative background, there is a higher-degree analogue of the splitting (2.31): on restriction to the curve Z(σ) there exists some g(U, σ) such that a| Z(σ) =∂| Z(σ) g, with for ξ ∈ P 1 an arbitrarily chosen point needed to define the inverse of the∂-operator acting on sections of O over P 1 . With this, it follows that the holomorphic frame is defined using the exponential map of g: For an external gluon with colour T a i and charge e i with respect to the background, it thus follows that With this, the formula for the colour-ordered amplitude on a Cartan-valued SD radiative background is Remarkably, all background dependence is captured by a simple exponential factor. Unlike the RSVW formula (4.4), it is impossible to perform any of the moduli integrals in (4.10) on a generic background, even with momentum eigenstates for the external gluon multiplets. This is because g(U, σ) introduces generic dependence on the moduli for a general SD radiative background. Nevertheless, the formula (4.10) has the remarkable property of posessing fewer integrals than expected from space-time perturbation theory. On a generic SD radiative background, one expects 4(n − 2) (bosonic) integrals for a npoint tree-level gluon amplitude. Evaluated on momentum eigenstates (4.3) -now valued in End E -our formula has 4d residual integrals that are not fixed by delta functions; but d ≤ n − 3 for all n > 3 and thus the background-dressed RSVW formula has fewer residual integrals than space-time perturbation theory in every helicity sector.
If one extracts only external gluons from the N = 4 multiplets in (4.10), it is possible to reduce the number of moduli integrations through a judicious reparametrization of the map moduli. Recall that for a N k MHV amplitude, k = d + 1 of the external gluons are negative helicity; without loss of generality, let gluons r = 1, . . . , d + 1 be negative helicity, with the remainder positive helicity. The parametrization of Z I (σ) in terms of the moduli U I a(d) can now be traded for where the (d + 1) points in twistor space {Z I r } parametrize the new moduli. Note that in this parametrization, the map has the property that Z I (σ r ) = Z I r . Substituting this parametrization into (4.10), the fermionic components of the Z I r appear linearly in the exponentials of the negative helicity gluons, and can be integrated out directly to leave only bosonic integrals. This gives the N d−1 MHV formula as Dσ j (j j + 1) , (4.13) with the product corresponding to whatever colour-ordering is desired. Now insert momentum eigenstates with the form of the negative helicity states following from the Z I (σ r ) = Z I r property of the parametrization (4.11).
The moduli integrations in d 2 λ r can now be performed against the delta functions in the negative helicity wavefunctions to give: Thus, the g i appearing in the exponents are Cartan-valued functions of {κ r , µ r , u r , σ r } as well as σ i . In this formula, we obtain 2(n − d − 1) of the refined/polarized four-dimensional scattering equations 10 [24,88], written in this parametrization as The remaining 2(d + 1) refined scattering equations would have arisen from the d 2 µ r integrals, but these cannot be performed analytically due to the non-trivial background. In (4.15) there are 2(n + d + 1) integrations against 2(n − d − 1); taking into account the GL(2, C) quotient this leaves 4d free integrations. Of course, this matches the number of residual integrals in the general formula (4.7), but this parametrization enables us to explicitly remove 2|4(d + 1) of the integrals against delta functions. A further 2(n − d − 1) integrals can be explicitly removed against the remaining delta functions in (4.15) by solving (4.17) for σ i in terms of the {u r , σ r , κ r , κ i }, but the resulting formulae are not so enlightening. Note that the same computations can be performed starting from (4.7) without assuming that a lies in a Cartan subalgebra, as momentum eigenstates are given by the expression (4.3) on any SD radiative background.
Evaluation on SDPWs: The formula (4.10) can be further simplified when the background is restricted to be a SDPW. In this case, recall that the background partial connection pulled back to the curve Z(σ) is (4.18) From this, we see that the only contributions to (4.8) arise when ι λ(σ) = 0, so it makes sense to parametrize the moduli of the map λ α (σ) in a way that makes this manifest. Now, ι λ(σ) is a homogeneous function of degree d on P 1 ; in an affine patch σ α = (1, z) this means that ι λ(z) is a degree d polynomial in z. Let {b 1 , . . . , b d } denote the d roots of this polynomial, and b 0 encode the overall scale so that: These {b 0 , b 1 , . . . , b d } provide an alternative parametrization of the moduli of ι λ(σ) ; written in homogeneous coordinates this parametrization is simply with each b r ∈ P 1 constituting only a single degree of freedom in the moduli measure. This parametrization also induces a natural choice of basis for H 0 (P 1 , O(d)): Finally, the µα(σ) components of (4.1) can be decomposed as with {m a(d) , m a(d) } the 2d + 2 moduli of the map. Feeding these parametrizations of the map into (4.8), one obtains . (4.25) The d + 1 moduli m a(d) of [ι µ(σ)] can be conveniently reparametrized by x r ν r s r (σ) . This parametrization makes it clear that g depends on only the {x r } moduli of the µα map: with dependence on the moduli of the λ α map implicit on the left-hand side. 11 Inserting (4.28) into the general formula and evaluating on momentum eigenstates, the moduli integrals inm a(d) and y as well as all of the fermionic moduli for χ a can be performed to yield delta functions. The result is: where J(b, ν) is the Jacobian resulting from the various reparametrizations performed on the map moduli: for |b 0 b 1 · · · b d | the Vandermonde determinant, and is the generalized Volkov exponential associated with degree-d maps into twistor space. Compared to (4.10) for a general SD radiative background, there are fewer residual moduli integrals in a SDPW. It is easy to see that the collection of delta functions in (4.29) imply momentum conservation in the x + , z andz-directions, as expected for scattering amplitudes in a plane wave background. Taking into account this 3-momentum conservation and the GL(2, C) quotient, there are 2d − 1 residual integrals in A n,d that are not saturated by delta functions. This counting matches the number of integrations expected for a N d−1 MHV amplitude computed with MHV diagrams in a SDPW background: one integral for each MHV vertex and one integral for each scalar propagator.
While (4.10) -and consequently (4.29) -are conjectural, it is easy to see that these formulae pass several non-trivial tests of their validity. Indeed, background gauge invariance and the correct trivial background and perturbative limits follow immediately by straightforward generalisations of the arguments of section 3.3. Furthermore, when d = 1 the formulae are equal to the MHV expressions (3.22) and (3.24). In the latter case, this follows after using the GL(2, C) quotient to fix b 0 = 1, b α 1 = ι α , ν 0 = 0 and ν 1 = −1, and setting ξ α = o α in the d = 1 version of (4.29).
Finally, we re-emphasize that N = 4 supersymmetry is playing a trivial role here: pure gluon amplitudes are extracted from these formulae algebraically, in exactly the same way as in a flat background, by expanding the fermionic delta functions as described in (4.5).

Further discussion
In this paper, we proved the gauge theory formulae on a strong self-dual plane wave background announced in [22] from first principles and extended it to an arbitrary self-dual radiative background. In addition, we gave conjectural formulae for the full tree-level gluon S-matrix on any self-dual radiative background; these require many fewer integrations than their corresponding counterparts computed via standard space-time Feynman diagrams in background perturbation theory. An explanation for why our formulae contain so many fewer integrations than what is naïvely expected is provided by a background-dressed version of the MHV formalism. In a trivial background, the MHV formalism is an alternative set of Feynman rules for perturbative Yang-Mills theory in which the only vertices are MHV vertices, where legs of opposite helicity are connected by a scalar propagator that includes an off-shell prescription for the un-dotted momentum spinors [25,89].
Applying background perturbation theory to the twistor action (3.14) of Yang-Mills theory defines the MHV vertices for such a background-dressed MHV formalism, and fixing an appropriate axial gauge (to eliminate the cubic MHV vertex) makes these the only vertices in the theory. In the special case of a SDPW background, this can be done explicitly and a formula for the propagator on twistor space can be found (see appendix C): where φ(x) is given by (2.56), having taken k off-shell, and the various delta functions are defined byδ This is structurally equivalent to the trivial background propagator, with the axial gauge defined by the spinorια associated with the SDPW [75,76], but with the extra phase factor e i(φ(x 1 )−φ(x 2 )) . It would be interesting to develop the MHV formalism on any SD radiative background, and to compare the formulae arising from such a 'disconnected' prescription to the 'connected' formulae presented here; Green's functions on more general self-dual backgrounds were already obtained explicitly in [90], albeit not in MHV gauge. If an exact matching were possible, then this would be a route towards proving our formulae; however, even on a trivial background the connected and disconnected prescriptions have only been partially linked (cf., [91]) -even through each is known to be correct. More generally, the usual unitarity-based arguments which prove the original RSVW formula cannot be directly applied to scattering amplitudes on a SD radiative background, as the functional freedom of the background means that even tree-level amplitudes are not rational functions of the kinematic data. Of course, unitarity and its consequences must still leave their imprint on the structure of amplitudes in a background, but it is not immediately clear how these can be operationalized in a way to constrain or prove whether a formula for a given scattering amplitude is correct. Nevertheless, it was recently shown that the analytic structure of 4-point amplitudes in strong field QED can be constrained by gauge invariance [92]; extending these findings to non-abelian gauge theory could provide a powerful new tool to check our formulae for N k≥1 MHV amplitudes on a SD radiative background.
One restriction in our work was to radiative backgrounds. This excludes charged background fields or instantons. The main reason for this was technical in the sense that such solutions have nontrivial topology. Nevertheless, the underlying twistor theory is well-developed for such solutions so that it wouldn't be too hard to extend our methods here. The application to amplitudes on instanton backgrounds was already addressed in [62]. Charged backgrounds might provide a starting point for applications to gravitational waves.
An obvious goal for future research is to find all-multiplicity expressions for the treelevel S-matrix of Yang-Mills theory on any radiative background, rather than the chiral cases considered here. An obvious route towards this is ambitwistor string theory [93], which is a non-chiral formalism that produces compact, all-multiplicity formulae for gauge theory scattering at tree-level and beyond on a trivial background [94]. Ambitwistor string theory can be coupled to a fixed background gauge field and is anomaly free (at genus zero) precisely when the background obeys the Yang-Mills equations [95]; furthermore, the fixed vertex operators for gluon perturbations on any gauge field background are known [96], and have been used to explicitly calculate 3-point amplitudes on a (real, non-chiral) plane wave in any number of space-time dimensions [97]. To extend these computations beyond 3-points an understanding of integrated vertex operators (or equivalently, the scattering equations) in a background gauge field is required.
Nevertheless, we observe that there are potential phenomenological applications of our formulae -even though the backgrounds are chiral. In particular, back-reaction effects from probes on strong gauge fields (cf., [98,99]) can be studied by considering scattering on a complexified background gauge field [100]. It would be very interesting to explore if our results can be of use in providing high-multiplicity data in the study of back-reaction in strong field QED or QCD.
Finally, in [22] we also presented all-multiplicity formulae for MHV scattering of gravitons in a strong self-dual gravitational plane wave. The proof and generalizations of those formulae will appear in a separate upcoming publication; the theory is quite distinct from the gauge theory case studied here.

Acknowledgments
TA is supported by a Royal Society University Research Fellowship. AS is supported by a Mathematical Institute Studentship, Oxford. LJM is partially supported by STFC grant ST/T000864/1.

A 3-and 4-particle checks on SD plane waves
In this appendix, we compute the 3-and 4-particle gluon and graviton MHV amplitudes in SDPWs directly from Feynman diagrams in background perturbation theory. This serves as a check on the MHV formulae derived in sections 3, which are shown to match the Feynman diagram computations for n = 3, 4.
The Feynman rules for perturbative Yang-Mills theory on a general gauge theory plane wave background (and in generic space-time dimension) can be found in [18,29], and are easily specialized to the case of a SDPW background. Vertices are simply read off from the background field Lagrangian, while the gluon propagator G ab ab (x, y) obeys where D a = ∂ a − i A a is the background covariant derivative, A a is the SDPW gauge field, F ab is its field strength, and e is the charge of the propagated gluon with respect to the Cartan-valued background. The Green's function solving this equation is given by where N = −πi/(2π) 4 is a convenient normalisation factor. The phase factors φ(x) are defined by (2.56), now understood to be evaluated at the off-shell momentum k. The tensor structure D ab (x − , y − ) is The tree-level 3-point gluon amplitude is read off from the cubic interaction in the background field Lagrangian [29]: where E i are the background-dressed polarizations of the external gluons, K i are the dressed momenta and F 3 (x − ) is the 3-point Volkov exponent of (3.26). It is straightforward to see that A 3 vanishes when all three external gluons have the same helicity; for the MHV configuration we take gluons 1 and 2 to be negative helicity and gluon 3 to be positive helicity. Making a judicious choice for the reference spinor ξ α in the positive helicity polarization (2.48) (i.e., a judicious gauge choice) substantially simplifies the calculation; we use ξ α = κ 2 α so that In this case, E 1 · E 2 = E 2 · E 3 = 0, and the integrand in A 3 reduces to Now, in general, 3-momentum conservation in the x + , z,z-directions can be written in the equivalent forms using charge conservation n i=1 e i = 0. For n = 3 this implies that 1 3 [ι 3] = − 1 2 [ι 2], which in turn simplifies the 3-point MHV amplitude to which matches the formula (3.24) at n = 3.
For completeness, we note that the 3-point MHV is also easily computed from the background field Lagrangian, and is again simplified on the support of 3-momentum conservation: is not the parity conjugate of A MHV 3 , as the former contains background dressed spinor contractions. This is a result of the chirality of the SDPW background itself.
The tree-level 4-gluon amplitude on a SDPW background recieves contributions from four Feynman diagrams [18]: where the first three terms represent the s, t and u channels, and the last term is the contact diagram coming from the quartic interaction. The simplest of these is the contact contribution: (A.10) The exchange diagrams contain a propagator insertion; for example, the s-channel is given by where the propagator tensor structure D s ab is defined with respect to exchanged momentum k 1 + k 2 and charge e 1 + e 2 . The integral is over two lightfront coordinates x − and y − , with the integration measure There are similar formulae for the t-and u-channels [18].
A priori, these expressions look a long way off our MHV formula (3.24), which contains only a single lightfront integral at arbitrary multiplicity. However, intricate cancellations in the integrands enable the application of some non-trivial integration-by-parts (IBP) identities (cf., [18]) which reduce all space-time expressions to a single lightfront integral. Let particles 1, 2 be negative helicity and 3, 4 be positive helicity, with dressed polarizations given by This gauge choice eliminates the contact contribution to the 4-point amplitude completely. The s-channel diagram evaluates to (A.14) To reduce this to a single lightfront integral, consider the following manipulation: Now, let G ab (x, y) be the propagator for an adjoint-valued scalar field of charge e with respect to the background: This identity can be used to replace the δ 4 (x − y) in (A.15), following which we integrate the charged Laplacian D 2 by parts. The crux of the resulting IBP identity is the relation is simplified to follows by exchanging particle labels 3 and 4. Note that it is highly non-trivial to obtain integrands to which these IBP identities can be applied; this is certainly not guaranteed by any naïve property of the Feynman diagrammatics. The colour-ordered 4-gluon amplitude is now obtained by using the Jacobi identity

B Impulsive SDPWs
It is illustrative to consider our formula for the tree-level S-matrix in the limit where the SD radiative background gauge field becomes an impulsive SDPW. This means that the background field strength has support only along a single lightfront hypersurface, representing an instantaneous burst of chiral colour radiation. This setting enables the x r moduli integrals of (4.29) to be performed explicitly. An impulsive SDPW located on the lightfront x − = 0 is described by a gauge field A = −α Θ(x − ) dz, whereα is a constant vector in a Cartan subalgebra of the gauge group and Θ(x − ) is the Heavyside step function. It is straightforward to show that the function F(x r ) appearing in g(x, σ) through (4.28) is given by F(x r ) =α x r Θ(x r ) for each r = 1, . . . , d. Each dx r integral in the formula (4.29) for A n,d can now be performed by splitting +∞ −∞ dx r (· · · ) = 0 −∞ dx r (· · · ) + +∞ 0 dx r (· · · ) , where an appropriate iε-prescription is assumed to ensure convergence of the integrals. These integrals can now be performed straightforwardly, with the result: where the impulsive factor I n,d is defined by: In this expression, the role of the small parameter ε is to regulate the flat background limit of the amplitudes. Observe that in the weak field limit, whereα → 0, the impulsive factor I n,d obeys

C Twistor action and MHV formalism
The MHV formalism can be derived from the twistor action (3.14) description of perturbative Yang-Mills theory. For our purposes, we need a version of this action coupled to a SD background gauge field on space-time, which is given by shifting a → a + a, whereD =∂ + a gives an integrable background complex structure on the Ward bundle E → PT, and H solves (3.16). The MHV vertices can be computed by expanding the second term of (C.1) as a series in a around the background a. This calculation was already done from a generating functional viewpoint to arrive at (3.21). Thus, to furnish MHV rules, we are left to derive the MHV propagator. In this appendix, we illustrate its derivation for a SDPW background (cf., [62,77] for the derivation in a trivial background).
The following analysis is best carried out in Euclidean signature on space-time, R 4 , so that the twistor space is identified with its projective spinor bundle PT ≃ R 4 × P 1 . In coordinates (x, λ) on this space, one can re-express the complex structure of twistor space as [54]∂ =ē 0∂ 0 +ēα∂α , where, introducing quaternionic conjugationλ α = (λ 1 , −λ 0 ), one defines a standard basis of (0, 1)-forms,ē 0 = Dλ as well as a basis of (0, 1)-vector fields, The Ward bundle also gets identified with the pullback of the Yang-Mills bundle by the natural projection PT → R 4 . The gauge connection A a on R 4 lifts to a partial connection on E → PT, 12 a = −i λ α A ααēα .

(C.5)
This defines our background coupled complex structureD =ē 0∂ 0 +ēα λ α D αα , where D a = ∂ a − i A a is the space-time covariant derivative.
A similar solution can be found for b in momentum space, having accounted for homogeneities. With (C.11), (C.12) it is straightforward to construct the MHV propagator as a Green's function forD, giving (5.1).