Berends-Giele currents in Bern-Carrasco-Johansson gauge for F3- and F4-deformed Yang-Mills amplitudes

We construct new representations of tree-level amplitudes in D-dimensional gauge theories with deformations via higher-mass-dimension operators α′F3 and α′2F4. Based on Berends-Giele recursions, the tensor structure of these amplitudes is compactly organized via off-shell currents. On the one hand, we present manifestly cyclic representations, where the complexity of the currents is systematically reduced. On the other hand, the duality between color and kinematics due to Bern, Carrasco and Johansson is manifested by means of non-linear gauge transformations of the currents. We exploit the resulting notion of Bern-Carrasco-Johansson gauge to provide explicit and manifestly local double-copy representations for gravitational amplitudes involving α′R2 and α′2R3 operators.


JHEP02(2019)078
Contents Recent investigations of scattering amplitudes in gauge theories and gravity revealed a wealth of mathematical structures and surprising connections between different theories. For gravitational theories in D spacetime dimensions, traditional methods for tree amplitudes and loop integrands naively give rise to an exasperating proliferation of terms. Still, the final answers for these quantities across various loop-and leg orders take a strikingly simple form: the dependence on the spin-two polarizations can often be reduced to squares of suitably chosen gauge-theory quantities.
The study of double-copy structures in perturbative gravity originates from string theory, where Kawai Lewellen and Tye (KLT) identified universal relations between openand closed-string tree-level amplitudes [1]. The KLT relations have been later on reformulated in a field-theory framework by Bern, Carrasco and Johansson (BCJ) [2][3][4] such as to flexibly address multiloop integrands. In this way, numerous long-standing questions on the ultraviolet properties of supergravity theories have been resolved [5][6][7][8][9][10], bypassing the spurious explosion of terms in intermediate steps.
This double-copy approach to gravitational amplitudes takes a particularly elegant form once a hidden symmetry of gauge-theory amplitudes is manifested -the duality between color and kinematics due to BCJ [2]. At tree level, the BCJ duality in gauge theories has not only been explained and manifested in string theories [11][12][13][14][15][16] but also extends to various constituents of string-theory amplitudes [17][18][19][20][21][22]. In particular, the following terms in the gauge-field effective action of the open bosonic string 1 in D spacetime dimensions preserve the BCJ duality to the order of α 2 [18], where F µν and α denote the non-abelian field strength and the inverse string tension, respectively. In presence of the effective action (1.1), KLT formulae and BCJ doublecopy representations known from Einstein gravity extend 2 to gravitational tree amplitudes 3 from α R 2 + α 2 R 3 operators [18] involving higher powers in the Riemann curvature R. The schematic notation R 2 and R 3 for operators in the gravitational effective action is understood to comprise additional couplings of a B-field and a dilaton ϕ (such as e −2ϕ R 2 ) known from the low-energy regime of the closed bosonic string [25]. The interplay of higher-mass-dimension operators D 2m F n and D 2m R n in string theories with the BCJ duality and double copy is well understood from the worldsheet de-JHEP02(2019)078 of the NLSM 5 [46] has been recently obtained from the analogous dimensional reduction of α F 3 in a companion paper [47]. In view of these connections, we hope that the notion of BCJ gauge inspires a reformulation of the (YM+F 3 +F 4 )-Lagrangian (1.1) wheresimilar to [44] -the D-dimensional kinematic algebra is manifest. 6 Another source of motivation for this work stems from the renewed interest in the gravitational α R 2 + α 2 R 3 interactions in D = 4 dimensions. While R 3 is well-known to be the first (non-evanescent) two-loop counterterm for pure gravity [52,53], the evanescent one-loop counterterm R 2 was recently found to contaminate dimensional regularization at two loops [54,55]. Moreover, evanescent matrix elements of R 2 are closely related to certain anomalous amplitudes of N = 4 supergravity [56] through the double copy [57]. Finally, when viewed as ambiguities in defining quantum theories, matrix elements of higher dimensional operators can be crucial to restore symmetries when using a non-ideal regulator for loop amplitudes [58]. We hope that our D-dimensional double-copy representations for tree-level amplitudes of (α R 2 + α 2 R 3 )-deformed gravity shed further light into these looplevel topics: either by unitarity or by using the BCJ-gauge currents as building blocks for loop amplitudes that universally represent tree-level subdiagrams. 7

Outline
This work is organized as follows: in section 2, we review the basics of Berends-Giele recursions, the BCJ duality as well as the double copy and establish the associated elements of notation. Section 3 is dedicated to amplitudes of (YM+F 3 +F 4 ) in different types of Berends-Giele representations including a systematic reduction of the rank of the currents. In section 4, an explicit off-shell realization of the BCJ duality is obtained from the Berends-Giele setup. Finally, section 5 relates this realization of the BCJ duality to non-linear gauge freedom and combines the off-shell ingredients from the previous section to manifestly local amplitude representations of (YM+F 3 +F 4 ) and gravity with α R 2 + α 2 R 3 operators. A derivation of the BCJ relations to the order of α 2 from purely kinematic arguments is given in section 5.2.

Review and notation
In this section, we set up notation and review the key ideas and applications of Berends-Giele recursions for tree-level amplitudes in YM theory, in particular • the resummation of Berends-Giele currents to obtain perturbiner solutions to the non-linear field equations • manifestly cyclic Berends-Giele representations of YM amplitudes involving currents of smaller rank than naively expected.

JHEP02(2019)078
We will also review the BCJ duality and the double copy from a perspective which later on facilitates the implementation of these features in tree amplitudes and Berends-Giele currents of (YM+F 3 +F 4 ) as well as gravity with α R 2 + α 2 R 3 operators.

Berends-Giele recursions
An efficient approach to determine the tensor structure of D-dimensional tree amplitudes in pure YM theory has been introduced by Berends and Giele in 1987 [29]. The key idea of the reference is to recursively combine all color-ordered Feynman diagrams involving multiple external on-shell legs and a single off-shell leg. This recursion is implemented via currents J µ 12...p that depend on the polarization vectors e µ i and lightlike momenta k µ i of the external particles i = 1, 2, . . . , p subject to the following on-shell constraints e i · k i = k i · k i = 0 ∀ i = 1, 2, . . . . Currents of arbitrary multiplicity can be efficiently computed from the Berends-Giele recursion [29] The external states have been grouped into multiparticle labels or words P = 12 . . . p. We will represent multiparticle labels by capital letters P, Q, X, Y, . . . and denote their length, i.e. the number of labels in P = 12 . . . p, by |P | = p. The summation over XY = P on the right-hand side of (2.2) instructs to deconcatenate P into non-empty words X = 12 . . . j and Y = j+1 . . . p with j = 1, 2, . . . , p−1 and therefore generates |P |−1 terms. 8 Similarly, XY Z = P encodes 1 2 (|P |−1)(|P |−2) deconcatenations into non-empty words X = 12 . . . j, Y = j+1 . . . l and Z = l+1 . . . p with 1 ≤ j < l ≤ p−1.
Moreover, the right-hand side of (2.2) involves multiparticle momenta k P through Mandelstam invariants or inverse propagators s P k µ P =12...p = k µ 1 + k µ 2 + . . . + k µ p , s P = 1 2 k 2 P . (2.5) Finally, the brackets in (2.3) and (2.4) capture the cubic and quartic Feynman vertices of pure YM theory in Lorenz gauge. As depicted in figure 1, the role of the deconcatenations XY = P and XY Z = P in (2.2) is to connect lower-rank currents J µ X , J ν Y and J λ Z via Feynman vertices in all possible ways that preserve the color order of the on-shell legs in the word P = 12 . . Accordingly, color-ordered on-shell amplitudes at n = p+1 points are recovered by taking the off-shell leg in the rank-p current J µ P on shell: this on-shell limit is implemented by contraction with the polarization vector J µ n = e µ n of the last leg and removing the propagator s −1 12...p in the p-particle channel of J µ P which would diverge by n-particle momentum conservation k 2 12...p → (−k n ) 2 = 0 [29], 9 A YM (1, 2, . . . , n−1, n) = s 12...n−1 J µ 12...n−1 J µ n . (2.6) For instance, the rank-two current due to (2.2) with X = 1 and Y = 2 yields the following representation of the three-point amplitude where cyclicity may be manifested via e 3 ·k 2 = −e 3 ·k 1 by means of on-shell constraints and momentum conservation. Note that Berends-Giele formulae similar to (2.6) have been given

Perturbiners as generating series of Berends-Giele currents
The Berends-Giele construction of the previous section can be related to solutions of the non-linear field equations: generating series of Berends-Giele currents turn out to solve the equations of motion from the action S YM of pure YM theory We use the following conventions in deriving the Lie-algebra valued gluon field A µ and its non-linear field strength F µν from a connection ∇ µ , The relation of tree-level amplitudes with solutions of the field equations via generating series goes back to the "perturbiner" formalism [68][69][70][71][72]. In these references, generating series of MHV amplitudes are derived from self-dual YM theory, see [73] for supersymmetric extensions. The connection between perturbiner solutions and the dimension-agnostic Berends-Giele currents of [29] was established in [27,28] and will now be reviewed. 10 The shuffle product P ¡Q of words P = p1p2 . . . p |P | and Q = q1q2 . . . q |Q| is recursively defined by P ¡∅ = ∅¡P = P , P ¡Q = p1(p2 . . . p |P | ¡Q) + q1(q2 . . . q |Q| ¡P).

JHEP02(2019)078
Lorenz gauge ∂ µ A µ = 0 simplifies the equations of motion [∇ µ , F λµ ] = 0 to the wave equation with the notation = ∂ µ ∂ µ for the d'Alembertian, One can derive formal solutions to (2.13) by means of the perturbiner ansatz (2.14) The summation variables i, j, l, . . . = 1, 2, 3, . . . refer to external-particle labels in an unbounded range, and we have introduced a compact notation P =∅ for sums over nonempty words P = 12 . . . p in passing to the second line. The dependence on the spacetime coordinates x µ enters through plane waves 11 e k P ·x , see (2.5) for the multiparticle momenta k P . The color degrees of freedom in (2.14) are represented through matrix products of the Liealgebra generators t a i whose adjoint indices a 1 , a 2 , . . . are associated with an unspecified gauge group. Upon insertion into the second line of (2.13), the perturbiner ansatz (2.14) can be verified to solve the non-linear field equations [∇ µ , F λµ ] = 0 if its coefficients J µ P obey the Berends-Giele recursion (2.2). Hence, generating series of Berends-Giele currents are formal solutions to the field equations. 12 By the shuffle symmetry (2.8) of the currents J µ P , the matrix products t a i t a j of the Lie-algebra generators on the right-hand side of (2.14) conspire to nested commutators, and the perturbiner solution is guaranteed to be Liealgebra valued [74].
As a convenient reorganization of the Berends-Giele recursion (2.2), one can write the field equations as in the first line of (2.13) and insert a separate perturbiner expansion for the non-linear field strength, The conventional form of plane waves e ik·x with an imaginary unit in the exponent can be recovered by redefining the momenta in this work as k → ik. The equations in the main text follow the conventions where external momenta are purely imaginary in order to keep factors of i from proliferating. 12 Strictly speaking, contributions with several factors of t a j referring to the same external leg j need to be manually suppressed by adding nilpotent symbols to the perturbiner ansatz [68]. For ease of notation, we do not include these symbols into the equations in the main text, and all terms with repeated appearance of a given external leg are understood to be suppressed.  ..p , only cubicvertex diagrams have to be considered in their recursive construction from lower-rank currents. In the depicted example at rank p = 3 with an additional off-shell leg. . . , only two cubic diagrams of s-channel and t-channel type contribute to the four-point amplitude obtained from s 123 J µ 123 J µ 4 .
The trilinear term {J X , J Y , J Z } in (2.4) which represents the quartic vertex of the YM Lagrangian has been absorbed into the non-linear part of the field-strength current B µν P in (2.15). The leftover deconcatenations P = XY in (2.16) can be interpreted as describing cubic diagrams, see figure 1. Let us illustrate this statement with the four-point amplitude s 123 J µ 123 J µ 4 derived from a rank-three current via (2.6): the two deconcatenations (X, Y ) = (12,3) and (1,23) in the recursion (2.16) for J µ 123 can be viewed as the two cubic diagrams in figure 2 where appropriate contributions from the quartic vertex (2.4) are automatically included.
Note that the Lorenz-gauge condition and the field equations imply the relations including transversality of the gluon polarizations for single-particle labels P = i. Moreover, the non-linear gauge symmetry of the action (2.11) under δ Ω A µ = ∂ µ Ω − [A µ , Ω] and δ Ω F µν = −[F µν , Ω] acts on the currents via The scalar currents Ω P are defined by the perturbiner expansion Ω(x) = P =∅ Ω P t P e k P ·x of the gauge scalar in δ Ω . We will later on spell out a choice of gauge-scalar currents Ω P which manifests the BCJ duality at the level of Berends-Giele currents. Another specific choice of Ω P → Ω lin P allows to track the effect of linearized gauge transformations e µ i → k µ i on the i th leg of the Berends-Giele currents in (2.18): one can line up the replacement e µ i → k µ i with a set of gauge transformations that preserves Lorenz gauge. The condition δ Ω lin (∂ µ A µ ) = ∂ µ (δ Ω lin A µ ) = 0 then translates into the recursion [27] which needs to be supplemented with the initial conditions Ω lin j → δ i,j if the linearized gauge transformations e µ i → k µ i only applies to the i th leg. Precursors of the formula (2.19) for linearized gauge transformations of Berends-Giele currents can be found in [65]. Figure 3. Diagrammatic interpretation of the building block M X,Y,Z in (2.20) with multiparticle labels X = x 1 x 2 . . . x p , Y = y 1 y 2 . . . y q and Z = z 1 z 2 . . . z r .

Manifestly cyclic reformulation
Given that the Berends-Giele formula (2.6) for color-ordered amplitudes A YM (1, 2, . . . , n) singles out the last leg n which is excluded from the current J µ 12...n−1 , cyclic invariance in the external legs is obscured. We shall now review a reorganization of the Berends-Giele currents for YM tree amplitudes such that the n th leg enters on completely symmetric footing. Moreover, the subsequent rewritings reduce n-point amplitudes to shorter Berends-Giele currents of rank ≤ n 2 instead of the rank-(n−2) currents in the recursion (2.2) for J µ 12...n−1 . The backbone of the manifestly cyclic Berends-Giele formulae is the building block [28] composed of three currents with multiparticle labels X, Y, Z each of which represents treelevel subdiagrams. The resulting diagrammatic interpretation of M X,Y,Z is depicted in figure 3, and the definition (2.20) along with B µν X = −B νµ X implies permutation antisymmetry M X,Y,Z = −M Y,X,Z and M X,Y,Z = M Y,Z,X expected from the cubic vertex in the figure.
Using k P · J P = 0 and k X + k Y + k Z = 0, it was shown in [28] that the n-point amplitude (2.6) can be rewritten as which will be referred to as "integration by parts" 13  .
Note in particular that the rank of the currents in the manifestly cyclic n-point amplitudes (2.24) is bounded by 14 n 2 rather than n−2 as expected from the recursions (2.2) or (2.16) for J µ 12...n−1 . In section 3, similar expressions with manifest cyclicity and Berends-Giele currents of maximum rank n 2 will be given for the deformed (YM+F 3 +F 4 ) theory.

BCJ duality
The organization of the Berends-Giele recursion (2.16) in terms of cubic-vertex diagrams as exemplified in figure 2 resonates with the BCJ duality between color and kinematics [2]: according to the BCJ duality, scattering amplitudes in non-abelian gauge theories can be represented in a manner such that color degrees of freedom can be freely interchanged with the kinematic variables. While "color" refers to contractions of structure constants f a i a i a k , polarizations and momenta are referred to as "kinematics", and the notion of "freely interchanging" will be shortly made precise. The three-index structure of the contracted structure constants can be visualized via cubic-vertex diagrams with a factor of f a i a i a k for each vertex and contractions of the adjoint indices along the internal edges. Similarly, the kinematic dependence on e µ i , k µ i should also be organized in terms of cubic diagrams to manifest the BCJ duality.
The non-linear Y of the field-strength current B µν P in (2.15) absorbs the contributions from the quartic vertex Tr[A µ , A ν ][A µ , A ν ] in the YM action (2.11). This can be seen from that fact that the non-linear terms have fewer propagators than the rest of (2.16). Hence, the use of field-strength currents amounts to inserting 1 = k 2 P k 2 P such that a quartic vertex is "pulled apart" into two cubic vertices connected by the "fake" 13 This terminology goes back to the fact that the building block (2.20) and the amplitude representation (2.21) descend from ten-dimensional SYM [28,38]: in the setup of these references, (2.22) is a consequence of BRST integration by parts in pure-spinor superspace [39].
14 Earlier examples of such economic and manifestly cyclic Berends-Giele representations have been investigated in [76], but the construction in the reference requires a mixture of quadratic, cubic and quartic combinations of Berends-Giele currents instead of a single building block (2.20). propagator k 2 P . The choice of the channel P in 1 = has to be compatible with the color dressing f abe f ecd of the quartic vertex, where ambiguities arise from the Jacobi relations In figure 4, this situation is visualized in a four-point tree-level context, but there is no limitation to cubic-diagram parametrizations of n-point tree amplitudes as well as multiloop integrands [3,4]. Although the BCJ duality conjecturally applies to loop integrands [3,4], we shall focus on its well-established tree-level incarnation.
Of course, contributions from the higher-order vertices of (Tr F 3 )-and (Tr F 4 )-type can also be cast into a cubic-graph form by repeated insertions of 1 = . For the action (1.1) of (YM+F 3 +F 4 ), the color structure of the F 3 and F 4 operators also boils down to contracted structure constants [18], and the ambiguities due to Jacobi identities (2.25) arise in this situation as well. In the subsequent review of the BCJ duality, the color-dressed tree-level amplitudes Tr(t a 1 t a ρ(2) t a ρ(3) . . . t a ρ(n) )A(1, ρ(2), ρ(3), . . . , ρ(n)) (2.26) may refer to pure YM (A → A YM ), to its (α F 3 + α 2 F 4 )-deformation (A → A YM+F 3 +F 4 ) or to any other generalization that obeys the BCJ duality. Once the kinematic dependence of (2.26) is absorbed into cubic diagrams I, J, K, . . . , one can choose a parametrization [2] M n = where Γ n denotes the set of cubic tree-level graphs with n external legs. The color factors C I represent the contracted structure constants that arise from the traces in (2.26). The kinematic numerators N I are combinations of e µ i and k µ i that can be assembled from the Berends-Giele currents of the theory. Finally, the propagators s −1 e comprise Mandelstam variables (2.5) for the multiparticle momenta in the internal edges e of the graph I.
The parametrization (2.27) is said to manifest the BCJ duality if all the symmetries of the color factors C I carry over to the kinematic numerators N I . More specifically [2]:  Figure 5. Triplets of cubic graphs I, J, K whose color factors C · and kinematic factors N · are both related by a Jacobi identity if the duality between color and kinematics is manifest. The dotted lines at the corners represent arbitrary tree-level subdiagrams and are understood to be the same for all of the three cubic graphs. duality-satisfying representation (2.27), the kinematic numerators exhibit the same antisymmetry properties under flips: (2.28) • For each triplet of graphs I, J, K where the Jacobi identities (2.25) lead to the vanishing of triplets C I + C J + C K , the BCJ duality requires the corresponding triplet of kinematic numerators to vanish as well As visualized in figure 5, such triplets of cubic graphs only differ by a single propagator.
In later sections, we will construct local representatives of the kinematic numerators N I in (2.27) of (YM+F 3 +F 4 ) which do not exhibit any poles in s P and obey the BCJ duality up to and including the order of α 2 . By the Jacobi identities (2.25) of the color factors, the numerators are still far from unique after imposing locality, and generic choices at n ≥ 5 points will fail to obey some of the kinematic Jacobi relations (2.29). Hence, finding a manifestly color-kinematics dual parametrization (2.27) requires some systematics in addressing quartic and higher-order vertices via 1 = k 2 P k 2 P . The additional requirement of locality is particularly restrictive, and we will see that suitable gauge transformations (2.18) of the Berends-Giele currents in (YM+F 3 +F 4 ) give rise to local solutions, generalizing the construction in ten-dimensional SYM [28].
including the order of α 2 [18] by isolating suitable terms in the monodromy relations of open-string tree-level amplitudes [11,12]. For a variety of four-dimensional helicity configurations, kinematic numerators of YM + F 3 subject to Jacobi relations (2.29) can be found in [18]. We will derive generalizations to helicity-agnostic expressions in D dimensions and include the α 2 order of (YM+F 3 +F 4 ). Kinematic antisymmetry relations (2.28) and Jacobi identities (2.29) leave (n−2)! independent instances of N I . A basis of kinematic numerators under these relations can be assembled from the "half-ladder" diagrams depicted in figure 6 which are characterized by a fixed choice of endpoints 1 and n as well as permutations ρ ∈ S n−2 of the remaining legs 2, 3, . . . , n−1. We will denote the basis numerators of the half-ladder diagrams in figure 6 by N 1|ρ(2,3,...,n−1)|n and refer to them as "master numerators".

Double copy
The BCJ duality allows to convert cubic-graph parametrizations (2.27) of gauge-theory amplitudes into gravitational ones: once the gauge-theory numerators N I satisfy the same symmetry properties as the color factors C I (i.e. flip antisymmetry (2.28) and kinematic Jacobi identities (2.29)), then the double-copy formula In case of (YM+F 3 +F 4 )-numerators, the gravitational amplitudes descend from a deformation of the Einstein-Hilbert action by higher-curvature operators of α R 2 + α 2 R 3 [18] as seen in the low-energy effective action of the closed bosonic string [25], see section 5.4 for details. The tilde along with the second copyÑ I of the gauge-theory numerator N I indicates that the i th external gravitational state may arise from the tensor product of different polarization vectors e µ i andẽ µ i .

Perturbiners and Berends-Giele representations for F and F 4
In this section, we apply the Berends-Giele methods of sections 2.1 to 2.3 to the deformed (YM+F 3 +F 4 ) theory known from the low-energy regime of open bosonic strings. The tree-level amplitudes following from the action

JHEP02(2019)078
reproduce the leading orders α 0 , α 1 in the low-energy expansion of bosonic-string amplitudes [80], and a well-defined sector of the α 2 order: the effective action of both bosonic and supersymmetric open strings comprises an operator α 2 ζ 2 Tr F 4 which can be cleanly distinguished from (3.1) by its transcendental prefactor ζ 2 = π 2 6 [81 -83]. Up to and including the order of α 2 , the amplitudes computed from (3.1) obey BCJ relations (2.30) while the α 2 ζ 2 Tr F 4 -operator excluded from (3.1) is incompatible with the BCJ duality [18]. As explained in the reference, these BCJ relations to the order of α 2 rely on the interplay between single-insertions of the α 2 Tr F 4 -operator in (3.1) and double-insertions of α Tr F 3 .
More generally, the accompanying multiple zeta values are instrumental to identify the D 2m F n -operators in string effective actions that admit color-kinematics dual representations. For instance, the entire single-trace gauge sector of the heterotic string obeys BCJ relations [20]. The subsector of open-bosonic-string amplitudes compatible with the BCJ duality was identified in [21], and the amplitude contributions without any zeta-value coefficient were derived from a field-theory Lagrangian [22].
The subsequent Berends-Giele recursions for the amplitudes of (3.1) follow a two-fold purpose: on the one hand, they will be used to generate economic and manifestly cyclic amplitude representations along the lines of section 2.3. On the other hand, they set the stage for • an off-shell realization of the BCJ duality in section 4 • a kinematic proof of the BCJ relations in section 5.2 • a construction of manifestly local gauge-theory numerators subject to kinematic Jacobi relations in section 5.3.
All of these results hold to the order of α 2 and are based on a non-linear gauge transformation of the generating series of Berends-Giele currents similar to (2.18).

Berends-Giele recursions for F 3 and F 4
Our Berends-Giele approach to (YM+F 3 +F 4 ) follows the lines of section 2.2 to derive recursions for the currents from the non-linear equations of motion. The field variation of the action (3.1) is given by 15 and augments (2.11) by α -corrections. In Lorenz gauge ∂ µ A µ = 0, setting (3.2) to zero amounts to a wave equation analogous to (2.13), . The extensions of this lemma to Fµν -dependent quantities X follows straightforwardly from the Leibniz rule.

JHEP02(2019)078
Since we will only be interested in the amplitude contributions up to the order of α 2 , we can simplify (3.3) by dropping terms of order α 3 and higher. At the first order in α , this allows to replace This form of the field equations gives rise to an efficient Berends-Giele recursion: we will study formal solutions of (3.4) modulo α 3 that descend from a perturbiner ansatz where Lorenz gauge and the definition F µν = −[∇ µ , ∇ ν ] of the field strength imply In comparison to the perturbiners (2.14) and (2.15) of undeformed YM theory, the Berends-Giele currents have been renamed as J µ P → A µ P and B µν P → F µν P in order to distinguish these α -dependent quantities from the YM currents in (2.16) and (2.15), In the same way as the Berends-Giele recursion of undeformed YM theory benefits from field-strength currents B µν P , the perturbiner solutions to (3.4) are conveniently expressed in terms of the additional auxiliary currents By their definition in (3.8), the auxiliary currents are determined by A µ P and F µν P , in the same way as the F µν P in (3.6) boil down to the elementary currents A µ P . All the above currents can be shown to obey shuffle symmetry by repeating the arguments for the currents J µ P and F µν P of pure YM theory [27,65]. With the above definitions, the Berends-Giele recursion induced by the field equation (3.4) takes the simple form

JHEP02(2019)078
In analogy to (2.6), the leading orders α ≤2 of the tree amplitudes resulting from the action (3.1) are then given by For instance, the rank-two current due to (3.11) with X = 1 and Y = 2 and the resulting three-point amplitude read The last term of the three-point function illustrates a fundamental difference between the tensor structure of YM amplitudes and their F 3 corrections: contractions of the type (k·e) n do not occur in n-point amplitudes of YM [84][85][86]. Hence, the expressions for A YM+F 3 +F 4 in this work do not belong to the class of D 2m F n amplitudes that are accessible from the open superstring through a combination of color-ordered SYM trees [17]. By the shuffle symmetry (3.10) of the currents, the amplitude representation (3.12) can be used to demonstrate A YM+F 3 +F 4 (1, 2, . . . , n) to also obey the Kleiss-Kuijf relations (2.9). The α 1 -order of our results up to and including n = 6 points has been checked to match the D-dimensional CHY formulae of 16 [26]. The D = 4 helicity components of the CHY expressions 17 in turn have been verified to agree with the results of [88][89][90].
We emphasize that the form of the recursion in (3.11) only involves deconcatenations P = XY into two words X, Y rather than three-word expressions with P = XY Z as seen in (2.2). Hence, the amplitudes (3.12) naturally arise in a cubic-graph parametrization (2.27) as visualized in figure 7. The cubic-graph organization extends to the order of α 2 , although each term of the Tr F 4 vertex in (3.1) involves at least four powers of the A µ field. Still, the quartic-vertex origin of the terms G νµ|ρσ X F ν|ρσ Y in (3.11) is visible through the absence of single-particle currents G µν|λρ i = 0 since there is no deconcatenation XY = i in (3.9). Like this, the part G νµ|ρσ X F ν|ρσ Y of the recursion (3.11) can only contribute at minimum length |X| + |Y | = 3, i.e. to amplitudes (3.12) at multiplicity n ≥ 4.
Note that the equation of motion (3.2) translates into the following expression for the tensor divergence of F µν P :

Manifestly cyclic Berends-Giele representations
This section is dedicated to a reformulation of the Berends-Giele formula (3.12) for A YM+F 3 +F 4 such as to manifest cyclicity and to reduce the maximum rank of the Berends-Giele constituents (3.11) on the right-hand side. This amounts to identifying a deformation 16 We are grateful to Song He and Yong Zhang for providing the analytic expressions. 17 See [87] for a systematic study of the reduction of CHY formulae to D = 4 dimensions. ..p of rank p combine the diagrams and propagators expected in a color-ordered (p+1)-point tree amplitude of (YM+F 3 +F 4 ) with an off-shell leg. . . . Vertices marked with a white and black dot represent the first and second order in α on the right hand side of (3.11), i.e. the cubic-graph parametrization of α F 3 and

JHEP02(2019)078
In analogy with (2.21), the amplitude formula (3.12) can be rewritten as To first order in α , the equivalence with (3.12) is proven to all multiplicities in appendix A.1. At second order in α 2 , (3.16) has been checked analytically to multiplicity n = 6 and numerically up to and including n = 8. In the same way as we are only interested in the orders α ≤2 of the amplitudes (3.16), we will consistently drop terms at the orders α ≥3 in later equations of this work and skip the disclaimer O(α 3 ) for ease of notation. Amplitude representations with manifest cyclicity and lower-rank Berends-Giele currents can be obtained by an integration-by-parts property that takes the same form as (2.22),

JHEP02(2019)078
The rational prefactors of non-prime multiplicities n = 4, 6, 8, 9 avoid overcounting of cubic diagrams when combinations of currents are invariant under less than n cyclic shifts i → i+1. Integration by parts (3.17) can be used to bypass such prefactors in expressions like also see [28,37,38] for the antecedents of these representations in ten-dimensional SYM.
Similarly, the all-multiplicity series of cyclic representations can be imported from ten-dimensional SYM [17].

Gauge algebra of F 3 + F 4 building blocks
The action of non-linear gauge transformations , Ω] is not altered by the higher-mass-dimension operators in the action (3.1). Hence, given perturbiner components Ω P for the gauge scalars Ω, the α -deformed currents of the previous subsection follow the transformations of the YM currents (2.18), One can therefore verify non-linear gauge invariance of the amplitude formula (3.12) by repeating the arguments of the undeformed gauge theory: among the three terms in the gauge variation the first one vanishes by the Lorenz-gauge condition k 12...n−1 · A 12...n−1 = 0 and the second one due to transversality k n · A n = 0 (using momentum conservation k µ 12...n−1 = −k µ n in both cases). The currents in the second line of (3.23) have multiplicity |X|, |Y | ≤ n−2 and are therefore regular as s 12...n−1 → 0, so multiplication with s 12...n−1 causes this term to vanish as well. This rests on the reasonable assumption that the gauge scalars Ω P descend from a perturbiner and can only have poles in s Q for subsets Q ⊆ P .
By the Leibniz property of δ Ω , momentum conservation and the expression (3.14) for contractions of the form k µ P F µν P , one can infer the non-linear gauge transformation of the JHEP02(2019)078 building blocks (3.15). The result is most conveniently expressed in terms of a scalar quantity Ω X,Y,Z,W = Ω [X,Y,Z,W ] which is totally antisymmetric in four multiparticle labels X, Y, Z, W , While the α → 0 limit of (3.24) descends from BRST variations of superspace building blocks in ten-dimensional SYM theory [28], a proof to the first order in α is given in appendix A.3. At the order of α 2 , we have tested (3.24) to the order of |X| + |Y | + |Z| = 6. Based on the gauge algebra (3.24) and permutation antisymmetry Ω X,Y,Z,W = Ω [X,Y,Z,W ] , one can check the non-linear gauge invariance of the manifestly cyclic amplitude representations in (3.18), (3.19) and (3.21). In particular, this will be exploited in later sections to evaluate the M X,Y,Z in BCJ gauge which is tailored to manifest the BCJ duality via local numerators.

Kinematic Jacobi identities in off-shell diagrams
The purpose of this section is to manifest the BCJ duality between color and kinematics in off-shell diagrams of (YM+F 3 +F 4 ). We will construct local solutions to the kinematic Jacobi relations (2.29) in the subdiagram with an off-shell leg drawn in figure 9. This amounts to assigning kinematic numerators to the cubic-vertex diagram in the figure which share the symmetries of the associated color factors The adjoint indices a 1 , a 2 , . . . , a p refer to p on-shell legs, and an off-shell leg is associated with a free adjoint index b carried by the rightmost factor in (4.1). When identifying the dotted off-shell line in figure 9 with an external on-shell leg, we recover the half-ladder diagrams of figure 6 that define the master numerators at n = p+1 points. Accordingly, permutations of figure 9 in 2, 3, . . . , p will be associated with the master numerators in an off-shell setup: cubic diagrams which are not of half-ladder form or do not have leg 1 and the off-shell leg. . . at their endpoints can be reached from (p−1)! permutations of figure 9 through a sequence of Jacobi identities.
In specific examples of (4.1) at rank p ≤ 5, antisymmetry f a i a j a k = f [a i a j a k ] and Jacobi identities (2.25) imply the so-called Lie symmetries for the color factors, The ellipsis in the subscript of each term indicates that lower-rank symmetries in the first labels extend to higher rank. For instance,  Figure 9. This section is dedicated to constructing local and Jacobi-satisfying kinematic representatives for the depicted cubic diagram of (YM+F 3 +F 4 ). The notation a µ 12...p , f µν 12...p , f We will now describe the construction of local kinematic factors for (YM+F 3 +F 4 ) with the same Lie symmetries of (4.2) which imply kinematic Jacobi relations. The recursive procedure presented in this section closely follows the steps of [15], where local superspace building blocks with Lie symmetries have been constructed for ten-dimensional SYM.

Local multiparticle polarizations up to rank three
As we already saw for the Berends-Giele currents of the previous section, each cubic vertex of (YM+F 3 +F 4 ) may introduce powers of α 0 , α 1 or α 2 into the kinematic factors. We no longer distinguish these contributions from the individual vertices (as done by the white and black circles in figure 7 and 8) and collectively refer to all contributions at orders α ≤2 through the off-shell diagram in figure 9. We will start from the numerators in the Berends-Giele recursion (3.11) to construct solutions to the kinematic Jacobi identitiesi.e. realizations of the Lie symmetries in (4.2) -up to the order of α 2 .
Kinematic representatives for the diagram in figure 9 with Lie symmetries will be referred to as multiparticle polarizations and denoted by lowercase parental letters a µ 12...p , f µν 12...p , f µ|νλ 12...p and g µν|λρ 12...p . This notation will help to distinguish the local multiparticle polarizations from the Berends-Giele currents A µ P , F µν P , F µ|νλ P and G µν|λρ P with kinematic poles. In the same way as all the four species of Berends-Giele currents enter the cyclic building block M X,Y,Z in (3.15), we will later on see that the analogous four species of multiparticle polarizations can be combined to Jacobi-satisfying kinematic numerators N I in the sense of section 2.4.
At rank one, the local multiparticle polarizations are defined to match their Berends-Giele counterparts which include the transverse polarization vectors e µ i and do not exhibit any kinematic poles, The simplest multiparticle polarization a µ 12 at rank two is defined by isolating the numerator in the Berends-Giele current (3.13), where the absence of contributions at order α 2 is plausible by the valence of the Feynman vertices from α 2 F 4 . The alternative presentation of (4.4) as a µ 12 = s 12 A µ 12 generalizes to the following two-particle polarizations at higher mass dimension, The local multiparticle polarizations are still proportional to their Berends-Giele counterparts since the latter only describe a single cubic diagram, see figure 10. By the shuffle symmetry 123 . Indeed, antisymmetry a µ 123 = − a µ 213 in the first two indices is inherited from the property (4.6) of the rank-two input. However, the first non-trivial kinematic Jacobi identity for the triplet of cubic diagrams in figure 11 requires a µ 123 + a µ 231 + a µ 312 to vanish, which is not the case. Still, the obstruction takes a special form, where one can factor out the overall momentum k µ 123 and isolate a scalar quantity h 123 that captures the deviation from the Lie symmetries a µ 123 + a µ 231 + a µ 312 = 3k µ 123 h 123 . Amusingly, the explicit form of can be reproduced from the cyclic building block of (3.15). Already the left-hand side of (4. As we will see, the appearance of the overall momentum k µ 123 in the correction (4.10) to a µ 123 is essential to absorb the analogous improvements of Berends-Giele currents into a non-linear gauge transformation (3.22).
Given a multiparticle polarization a µ 12...p at rank p, the construction of its analogues f µν 12...p , f µ|νλ 12...p and g µν|λρ 12...p at higher mass dimension involves contact terms ∼ s ij that preserve the Lie symmetries. For instance, the rank-three generalizations of (4.5) are easily checked to reproduce the symmetries (4.11) of a µ 123 . These contact terms are the local equivalents of the deconcatenation terms in the Berends-Giele currents F µν P , F µ|νλ P and G µν|λρ P in (3.6) and (3.9), see section 5.1 for more details and [15] for superspace analogues.

Local multiparticle polarizations at rank four and five
The numerators of the Berends-Giele recursion (3.11) serve as a starting point to construct higher-rank multiparticle polarizations for the diagrams in figure 9 that satisfy the Lie symmetries (4.2) of the dual color factors. The higher-rank systematics of the redefinition in (4.10) is most conveniently illustrated via examples at p = 4, 5.
Given the multiparticle polarizations at rank three in (4.10) and (4.12), their Lie symmetries imply that the rank-four object (4.14) The pattern of subtractions in (4.14) has been inferred by mimicking BRST transformations in ten-dimensional pure-spinor superspace [15], and it should be possible to give a similar motivation from a study of linearized gauge variations. The key benefit of the redefinition in (4.14) is that the deviation from the rank-four Lie symmetry now takes a factorized form At higher rank, analogous redefinitions in two steps a µ 123...p → a µ 123...p → a µ 123...p will be sufficient to attain Lie symmetries, i.e. there are no additional intermediate steps at p > 4. For instance, the Lie-symmetry satisfying multiparticle polarizations (4.16) and (4.18) at rank four can be used to recursively construct a rank-five quantity subject to the symmetries (4.17) in its first four labels. The rank-five Lie symmetry can be enforced by first performing subtractions analogous to (4.14) and [15],

Local multiparticle polarizations at higher rank
The recursive construction of multiparticle polarizations will now be summarized in terms of all-rank formulae that closely follow their superspace antecedents [15] but incorporate α -corrections. The Lie symmetries of a µ 12...q , f µν 12...q , f of the summation prescription of the form a 1 a 2 . . . a k = X¡Y in (4.27): for a given kparticle label a 1 a 2 . . . a k , the sum runs over all the 2 k pairs of words X and Y whose shuffle product contains a 1 a 2 . . . a k , for instance 20 (X, Y ) = (∅, ∅) in case of k = 0 as well as (X, Y ) = (a 1 , ∅) and (X, Y ) = (∅, a 1 ) in case of k = 1. The outcome of (4.27) still obeys the Lie symmetries in the first p−1 labels and is claimed to factorize k µ 12...p when probing the rank-p Lie symmetry: : p = 2n even .
see the explanation above for the summation prescription j+1, j+2 . . . p = X¡Y . The pattern of contact terms on the right-hand side preserves the Lie symmetries in all the p labels and will be connected with Berends-Giele currents in section 5.1. In the next section, these four Jacobi-satisfying kinematic representatives (4.29) and (4.31) of the offshell diagram in figure 9 will be combined to on-shell numerators of (YM+F 3 +F 4 ).

Berends-Giele currents in BCJ gauge
In this section, we relate the Jacobi-satisfying numerators for cubic off-shell diagrams as constructed in the previous section to gauge-transformed Berends-Giele currents. The idea is to compare the Lorenz-gauge currents A µ P , F µν P , F on the right-hand side of (5.1). The rank-two currents are degenerate with A µ, BCJ 12 = A µ 12 , while the redefinitions of numerators a µ 12...p at ranks p ≥ 3 by h 12...p introduce differences between A µ P from the recursion (3.11) and the alternative currents in (5.1). At rank three, there are two cubic diagrams contributing to A µ, BCJ 123 after dropping the distinction between order-α 0 , α 1 , α 2 vertices in figure 7, and the five cubic diagrams at rank four are depicted in figure 12. The numerator for the last cubic diagram of A µ, BCJ 1234 in the figure with propagators (s 12 s 34 s 1234 ) −1 is defined to be a µ 1234 − a µ 1243 by its relation to half-ladder numerators via Jacobi identities. 21 More generally, each cubic diagram contributing to A µ, BCJ 12...p can be derived from the half-ladder topology via kinematic Jacobi relations. The half-ladder numerators at rank 21 The Lie symmetry a µ 1234 − a µ 1243 = −(a µ 3412 − a µ 3421 ) ensures that this numerator changes sign when trading 12 ↔ 34 by a flip of the central cubic vertex and therefore obeys (2.28).

JHEP02(2019)078
with gauge parameters Ω P → H P . These transformations are checked to apply to all currents up to and including rank five, see appendix B.3 for the explicit form of H 12345 , and the existence of suitable H 12...p is conjectural at higher rank p ≥ 6. As the punchline of (5.7), the local Jacobi-satisfying numerators for off-shell diagrams are related to Lorenz-gauge currents through a non-linear gauge transformation generated by H P t P e k P ·x . (5.8) In the next sections, the local multiparticle polarizations will be used to manifest the BCJ duality between color and kinematics in tree-level amplitudes of (YM+F 3 +F 4 ). Hence, the transformed currents A µ BCJ P , . . . , G µν|λρ, BCJ P related by (5.2) are said to be in BCJ gauge [27,28]. Note that the first perturbiner solutions to the field equations of tendimensional SYM were actually constructed in BCJ gauge [91].

Kinematic derivation of the BCJ relations
As a first application of BCJ-gauge currents, we derive the BCJ relations (2.30) of A YM+F 3 +F 4 amplitudes by inverting their correspondence (5.2) with multiparticle polarizations. The same sequence of arguments has been applied to derive BCJ relations for tree amplitudes of ten-dimensional SYM from superspace currents in BCJ gauge [28], and we adapt the reasoning of the reference to the bosonic amplitudes up to the α 2 -order.
At rank p, inversion of (5.2) relates multiparticle polarizations to BCJ-gauge currents via where S(σ|ρ) 1 denotes the KLT matrix defined in (2.39). The Lie symmetries of the numerators of the BCJ-gauge currents ensure that the Mandelstam invariants from S(σ|ρ) 1 cancel all of their kinematic poles on the right-hand sides of (5.9). However, when repeating these matrix multiplications with Lorenz-gauge currents A µ P , some of the kinematic poles in (m ≥ 3)-particle channels persist 22 Hence, it is a peculiarity of BCJ-gauge currents that local objects are obtained from matrix multiplication with S(σ|ρ) 1 . Similarly, the pole s −1 12...p in the p-particle channel drops out from the following rank-p combinations of BCJ-gauge currents, This can be used to derive BCJ relations among color-ordered (YM+F 3 +F 4 ) amplitudes.
We exploit that the amplitude formula (3.12) is invariant under non-linear gauge transformations -see section 3.3 -and can therefore be written in terms of BCJ-gauge currents, This calls for an all-multiplicity formula for analogous combinations with regular behaviour as s 12...p → 0: the right-hand sides of (5.10) can be generated through the S-map [15,28] A µ, BCJ

JHEP02(2019)078
involving words P = p 1 p 2 . . . p |P | and Q = q 1 q 2 . . . q |Q| . BCJ gauge of the currents implies that the S-map defined in (5.15) removes the pole in s P Q [15] and therefore paves the way for the following form of the BCJ relations [28]  which is equivalent to (2.30) by the KK relations A YM+F 3 +F 4 ((X¡Y ), n) = 0 ∀ X, Y = ∅.

Local Jacobi-satisfying numerators
In this section, we will exploit the multiparticle polarizations of (YM+F 3 +F 4 ) to construct local and Jacobi-satisfying cubic-diagram numerators. The most direct approach is to expand the BCJ-gauge current in the amplitude representation (5.12) via (5.2), . From the remaining propagators in Φ(τ |ρ) 1 , the expressions a µ 1ρ(23...n−1) e µ n will be shown below to take the role of master numerators of the half-ladder diagrams depicted in figure 6. The (n−2)! KKindependent permutations of A YM+F 3 +F 4 in (5.18) incorporate each cubic diagram at least once and therefore define all of the numerators.
(5.19) Consistency with (5.18) can be conveniently checked by expressing Φ(τ |ρ) 1 with τ, ρ ∈ S n−2 as a putative (n+1)-point doubly-partial amplitude −m(1, τ, n, n+1|1, ρ, n+1, n) via (2.38) and (2.42). By its Berends-Giele representation (2.40) [63], the latter can be written as Figure 13. Diagrammatic interpretation of the local building block N X,Y,Z with multiparticle labels X = x 1 x 2 . . . x p , Y = y 1 y 2 . . . y q and Z = z 1 z 2 . . . z r referring to three off-shell half-ladder diagrams that are connected by the central vertex. Note that we no longer distinguish between the orders of α carried by the individual cubic vertices and therefore suppress the white and black dots of figure 8.
Since φ P |Q vanishes unless P is a permutation of Q, the only contribution arises from the deconcatenations with A = n and Y = n leading to Hence, (5.19) at σ = (1, τ, n) reduces to (5.18). For other choices of σ in turn, validity of (5.19) follows from the KK relations of both sides. Hence, by the discussion around (2.36), the cubic-diagram numerators of (5.18) are composed from the masters a µ 1ρ(23...n−1) e µ n as dictated by Jacobi identities. Note that the cubic-diagram numerators in (5.18) and (5.19) are not crossing symmetric, i.e. their functional form in terms of polarizations and momenta depends on the position of the singled-out legs 1 and n in the diagram.
In the same way as the manifestly cyclic representations of section 3.2 assemble npoint amplitudes from Berends-Giele currents of maximum rank n 2 , we will next spell out alternative numerators in terms of lower-rank multiparticle polarizations. In analogy to the cyclic building block M X,Y,Z in (3.15) composed of Lorenz-gauge currents, we define the local combination to describe the cubic diagram in figure 13 (see figure 8 for the analogous diagrammatic interpretation of M X,Y,Z ). There is an ambiguity in relating cubic diagrams to the combinations N X,Y,Z in (5. all of N 123,4,5 , N 12,3,45 and N 1,2,543 describe the same cubic diagram. A valid (n−2)!-set of N X,Y,Z to serve as the master numerators for half-ladder diagrams is given by As a defining property of these master numerators N ... , the central vertex of figure 13 is always chosen to be adjacent to leg n which therefore enters in a single-particle slot. As depicted in figure 14, the numerators in (5.23) describe half-ladder diagrams with endpoints 1 and n−1, where the location of leg n decides about the partition into the three subdiagrams associated with the slots of N X,Y,Z . The remaining labels a 1 , a 2 , . . . , a p , b 1 , b 2 , . . . , b q are a permutation of 2, 3, . . . , n−2 with p+q = n−3. Together with the n−2 different choices of p = 0, 1, . . . , n−3, this exhausts the total of (n−2)! permutations of the larger set 2, 3, . . . , n−2, n. The collection of N 1a 1 a 2 ...ap|n|bq...b 2 b 1 n−1 in (5.23) and figure 14 can be used as an alternative to the master numerators a µ 1ρ(23...n−1) e µ n in (5.19). As a practical advantage of the N ... , their constituents in (5.22) only require multiparticle polarizations of maximal rank n−2 instead of the rank-(n−1) quantities a µ 1ρ(23...n−1) . As demonstrated in appendix C, 23 they yield Jacobi-satisfying amplitude representations of the form (2.36), ρ(2, 3, . . . , j), n, ρ(j+1, . . . , n−2), n−1) . (5.24) Note that the α → 0 order of (5.24) follows from the field-theory limit of the pure-spinor superstring based on the amplitude representations of [14,63] and the superspace gauge described in [27].
Moreover, one can further reduce the maximum rank of the multiparticle polarizations by a generalization of the integration-by-parts relation (3.17). The latter still holds when

(5.26)
The right-hand sides of (5.25) and (5.26) at n ≥ 5 can be assembled from multiparticle polarizations of maximum rank n−3, and the spurious poles in s n−1,n cancel after combining all the terms. The same strategy applies to permutations of N 12...n−2|n|n−1 and N 1|n|23...(n−2)(n−1) in 2, 3, . . . , n−2. Like this, the n-point amplitude representation (5.24) with manifest BCJ duality and local numerators is completely determined by multiparticle polarizations of rank n−3. For instance, the explicit construction of multiparticle polarizations up to rank five in section 4 is sufficient to pinpoint all the eight-point numerators in (5.24). Note that the special footing of legs 1, n−1 and n in (5.24) breaks the crossing symmetry even more heavily than the numerators in (5.19). Still, one can restore crossing symmetry by averaging over all choices of singling out legs i, j, k ∈ {1, 2, . . . , n} instead of 1, n−1 and n.

Relation to string-theory and gravity amplitudes
A major motivation for the construction of (YM+F 3 +F 4 ) numerators with manifest locality stems from their connection with gravitational quantities through the double copy. Following the lines of [18], the double copy of A YM+F 3 +F 4 to amplitudes from highercurvature operators α R 2 +α 2 R 3 can be extracted from the string-theory KLT relations [1] (also see [23,24]

): the leading α -orders of the open-bosonic-string amplitudes
comprise our results for A YM+F 3 +F 4 (and the aforementioned contribution from the supersymmetrizable F 4 -operator which is incompatible with the BCJ duality and can be distinguished by its coefficient ζ 2 [18]). By the interplay with the trigonometric factors in JHEP02(2019)078 the KLT formula, both copies of A super−F 4 (σ) drop out from relevant orders of the closed bosonic string [21], (5.28) where S(ρ|τ ) 1 is the field-theory KLT matrix defined in (2.39) and the permutations ρ, τ act on 2, 3, . . . , n−2. Hence, to the orders considered, the right-hand side of (5.28) describes amplitudes from the low-energy effective action of the closed bosonic string [25] S closed where ϕ denotes the dilaton and H = dB is the field strength of the B-field. In spite of the dilaton admixtures via e −2ϕ , e −4ϕ , the operators along with the first and second order of α are collectively referred to as R 2 and R 3 . While the R 2 operator can be reconciled with up to sixteen supercharges [93], the R 3 operator is not supersymmetrizable [94].
Given the multitude of propagators in KLT formulae of the form (5.28), the locality properties of gravity amplitudes are more transparent in representations involving Jacobisatisfying numerators as in (2.34). For instance, our master numerators for (YM+F 3 +F 4 ) in (5.19) and (5.24) admit a realization of the double-copy structure via where each term has the propagator structure of cubic diagrams. The subscript GR + R 2 + R 3 is just a schematic shorthand for the amplitudes generated by the action (5.29) to the orders of α 2 . As emphasized in [18], the α 2 -order of (5.30) receives contribution from both single-insertions of R 3 operators and double-insertions of R 2 operators. In D = 4 spacetime dimensions, the R 2 contribution to (5.29) is the topological Gauss-Bonnet term. The components at the first order in α of (5.30) with graviton helicities are therefore guaranteed to vanish. Still, the double insertions of R 2 contribute to the α 2 -order of graviton components in four dimensions since the prefactor of e −2ϕ in (5.29) allows for dilaton exchange [18].
On the right-hand side of (5.28) or (5.30), the α 2 order receives both symmetric and asymmetric contributions: terms of the form A YM+F 3 +F 4 α Ã YM+F 3 +F 4 α where both gauge-theory halves contribute a factor of α have been carefully analyzed in D = 4 helicity components [18]. Our results on the α 2 -order of A YM+F 3 +F 4 and its master numerators

JHEP02(2019)078
additionally deliver the contributions to D-dimensional amplitudes of GR + R 2 + R 3 , where both powers of α stem from the same gauge-theory factor. These contributions involving A YM+F 3 +F 4 α 2 explain 25 the departure of M GR+R 2 +R 3 α 2 from the double copy of the first α -order A YM+F 3 +F 4 α which has already been observed in certain D = 4 helicity components [18]. Hence, there is no need to consider higher orders from the α -expansion of the sin( πα 2 k i ·k j ) terms in the string-theory KLT relations as speculated in the reference. At the order of α , one may extract new representations 26 for supersymmetrized matrix elements of R 2 from (5.30) by trading A YM+F 3 +F 4 for color-ordered amplitudes of ten-dimensional SYM and their dimensional reductions. These supersymmetrizations play a key role in recent studies of divergences and duality anomalies of N = 4 supergravity [57,58].

Conclusions and outlook
In this work, we have studied various representations for tree-level amplitudes of Ddimensional gauge theories with α F 3 + α 2 F 4 deformations. Our main results concern the emergence of the BCJ duality between color and kinematics up to the order of α 2 which is driven by the interplay of the F 3 and F 4 -operators [18] and independent on the choice of gauge group. While the BCJ duality has originally been explained by the realization of the α F 3 + α 2 F 4 operators from the open bosonic string, our work takes a different approach by identifying the seeds of the duality in the Berends-Giele currents of (YM+F 3 +F 4 ).
We study the Berends-Giele currents of (YM+F 3 +F 4 ) in the perturbiner formalism [68][69][70][71][72], where non-linear gauge transformations can be mapped to reparametrizations of the scattering amplitudes [27]. First of all, these currents are used to construct manifestly cyclic amplitude representations, where the rank of the contributing currents is systematically reduced. Moreover, we pinpoint a specific non-linear gauge transformation up to the order of five on-shell legs in a single current which rearranges the naive Feynman-diagram output of the action such as to manifest the BCJ duality. Like this, we derive the BCJ relations among color-ordered amplitudes to the order of α 2 from purely kinematic considerations. Furthermore, two kinds of explicit cubic-diagram parametrizations are given for (YM+F 3 +F 4 )-amplitudes where the manifestly local numerators obey kinematic Jacobi relations.
Our construction is inspired by superspace kinematic factors of ten-dimensional SYM [27,28] whose properties were inferred from the conformal-field-theory description of the pure-spinor superstring [14,15,37]. We identify extensions of these superspaceinspired structures to higher orders in α and to operators F 3 , F 4 that do not admit any supersymmetrization. It would be interesting to find a conformal-field-theory derivation of the local multiparticle polarizations that drive our BCJ-duality-satisfying amplitude rep- 25 We are grateful to Johannes Brödel for helpful discussions on this point and checking a representative four-dimensional helicity example. 26 Note that alternative representations with 8 supercharges on both chiral halves can be extracted from the low-energy limit of one-loop string amplitudes in K3 orbifolds [95,96].

JHEP02(2019)078
resentations. One possible starting point is to combine the worldsheet description of the bosonic string with the off-shell techniques of [16]. Alternatively, it might be helpful to identify a vertex-operator origin of the CHY formulae for F 3 amplitudes [26] 27 along with generalizations to higher orders in α . A complementary approach to the α F 3 + α 2 F 4 operators of the bosonic string is suggested by the recent double-copy description of bosonic-string amplitudes [22]: after peeling off the worldsheet integrals that are common with superstring amplitudes, an allorder family of α -corrections of the bosonic string can be traced back to a massive gauge theory dubbed (DF ) 2 + YM. The latter has been constructed in [97] by imposing the BCJ duality on a collection of dimension-six interactions between gauge bosons and massive scalars, and it should reproduce the (YM+F 3 +F 4 )-amplitudes in this work upon lowenergy expansion. It would be interesting to study our results from the (DF ) 2 + YMperspective and to generalize them to arbitrary orders in α by integrating out its massive modes.
Convenient and Jacobi-satisfying representations of tree-level subdiagrams are helpful for loop integrands of string-and field-theory amplitudes, see e.g. [59][60][61][62]. Our results might guide the organization of tensor structures of loop amplitudes in bosonic and heterotic string theories. This in turn could give input on loop integrands of half-maximal supergravity and their interplay with evanescent matrix elements and anomalies [54,55,57,58].

Acknowledgments
We are grateful to the organizers of the PSI Winter School 2018 where this work was initiated. We would like to thank Freddy

A Properties of the cyclic building blocks
In this appendix, we prove several properties of the building block M A,B,C defined in (3.15) to first order in α . These proofs are based on transversality k P · A P = 0 and the truncation 27 Note that the CHY half-integrands Pn in [26] with a puncture zj ∈ C on the Riemann sphere for each external state j = 1, 2, . . . , n and zi,j = zi − zj may be reproduced from the first order in α of Pn(zj, kj, ej) = a µ 12...n−1 e µ n z1,2z2,3 . . . zn−1,nzn,1 (α ) 1 + perm(2, 3, . . . , n−1) .

JHEP02(2019)078
of (3.14) to the first order in α , Moreover, we will use the generalization of (A.1) to higher mass dimension, which follows from the definition [∇ µ , ∇ λ ] = −F µλ and a corollary of the equations of motion, By virtue of (A.1), (A.2) and (A.4), one can rewrite any contraction of F µν P and F µ|νλ P with the corresponding momentum k P in terms of deconcatenations. We will always work to first order in α , but we will split the proofs into different orders for the convenience of the reader.

A.1 Appearance in the amplitudes
The first property we wish to prove is where, as a special case of (3.15), ) . (A.7) We first focus on the zeroth order of α . Notice how the last two terms on the first line of (A.7), when contracted with A µ n , exactly match the first and third terms on the first line JHEP02(2019)078 of (A.6). To see that the remaining terms are equal to each other, we notice that where we have used momentum conservation k µ n = −k µ X − k µ Y and transversality. This matches with the missing terms at the α 0 order of (A.7), upon contraction with A µ n . We now show that the same matching occurs between the terms of order α in both expressions. Equating them leads to When using the Jacobi identity (A.3), two of the terms on the left-hand side combine to cancel exactly the two on the right-hand side: Now we have to make sure that the remaining terms C, D, E, F and G on the left-hand side of (A.9) cancel each other. The first two can be rewritten as XY =12...n−1

JHEP02(2019)078
where the second line follows from (A.2). In passing to the third line, we have again used momentum conservation and exploited antisymmetry F νλ n = F [νλ] n . Hence, the statement (A.5) is proved to first order in α .

A.2 Integration by parts
The second property of M A,B,C we want to prove is the integration-by-parts identity which translates into the following claim at the zeroth order in α , (A.14) We can rewrite the second term ∼ A X A Y in the first line using the antisymmetry of F µν B and the definition (3.6) of F µν A : The analogous sum in the second line of (A.14) has a similar term related by A ↔ B, (A.17) Inserting this into (A.14) and ignoring the O(α )-term in the last line, we conclude that the property (A.13) is indeed satisfied to zeroth order in α . We now show that the property is still valid at the first order in α . In doing that, we have to combine the last term of (A.17) with the O(α )-terms in the M A,B,C of (A.13).

JHEP02(2019)078
Hence, the leftover task is to prove that Let us start by using the Jacobi identity (A.3) to rewrite the following terms: The right-hand sides cancel the first two terms inside the sums in (A.18) over A = XY and B = XY . Then, using the recursive definition of F µ|νλ P in (3.9), we can write and using the same manipulations The terms of the form k µ · F νλ · F ν|µλ · cancel between (A.22), (A.23) and (A.24) by momentum conservation k A +k B = 0. Finally, the contributions of the form A=XY F νλ B F µλ X F µν Y cancel between (A.23), (A.24) and the leftover terms of (A.18). This concludes our proof of the integration-by-parts identity (A.13) to the order of α .

A.3 Gauge algebra
Finally, we want to see how a non-linear gauge transformation (3.22) acts on M X,Y,Z . To zeroth order in α , we get (A. 25) Let us look at the terms which are not inside a sum. We can start by grouping all the ones with the same Ω · coefficient, and use momentum conservation k µ We then rewrite these four terms via

JHEP02(2019)078
and the same identities with (Y ↔ Z). In the first line, we rewrote k µ Y F µν Y using the zeroth order in α of (A.1). The last two lines of (A.27) are based on the antisymmetry of F µν Y and the definition (3.6) of F µν Z . Notice that F µν Y F µν Z cancels in the antisymmetrization w.r.t. Y ↔ Z in (A. 26). Hence, all the leftover terms in (A.25) involve a sum over deconcatenations, either X=AB or one of (X ↔ Y, Z). We collect all the expressions from the cyclic permutations in (A.25) where the sum is X=AB (A. 28) It turns out that the coefficient of each of the Ω's inside the above sum can be identified as some M P,Q,R with various combinations of the three words: The object inside the sum over X = AB is totally antisymmetric in A, B, Y, Z and can be identified as Ω A,B,C,D as defined in (3.24). Hence, the zeroth order of the gauge transformation (A.25) can be written as We now want to extend the proof of (A.30) to the first order in α . First of all, terms of O(α ) have been neglected when inserting (A.1) into the first term of (A.27). Therefore, we carry forward the following terms in δ Ω M X,Y,Z , where we have spelt out all the terms of the same form X=AB as in (A.28). This needs to be combined with the gauge variation of the O(α ) terms in the definition (3.15) of M X,Y,Z : F νλ Y ± perm(X, Y, Z) . (A.32) The only term in (A.32) which is not yet in the form of a deconcatenation sum will now be rewritten via momentum conservation k µ X + k µ Y + k µ Z = 0: (A.33) The first term calls for the relation (A.2), such that all the terms in the quantity G have been expressed via deconcatenation sums: F νλ Z ± perm(X, Y, Z) . δ Ω M X,Y,Z | α 1 = 1 2 In this appendix, we spell out two representations of the local rank-five scalar h 12345 that arises in the redefinition (4.21) towards the multiparticle polarization a µ 12345 . The scalar h 12345 can be expressed in terms of the local building blocks N X,Y,Z defined in (5.22) which are composed from multiparticle polarizations at rank ≤ 3,
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.