QCD amplitudes with 2 initial spacelike legs via generalised BCFW recursion

We complete the generalisation of the BCFW recursion relation to the off-shell case, allowing for the computation of tree level scattering amplitudes for full High Energy Factorisation (HEF), i.e. with both incoming partons having a non-vanishing transverse momentum. We provide explicit results for color-ordered amplitudes with two off-shell legs in massless QCD up to 4 point, continuing the program begun in two previous papers. For the 4-fermion amplitudes, which are not BCFW-recursible, we perform a diagrammatic computation, so as to offer a complete set of expressions. We explicitly show and discuss some plots of the squared 2 → 2 matrix elements as functions of the differences in rapidity and azimuthal angle of the final state particles.


JHEP02(2017)009
This paper is organized as follows. In section 2 we formally prove that it is possible to apply the BCFW recursion to tree level QCD scattering amplitudes with at least one gluon with up to two off-shell legs, using straightforward diagrammatic arguments. In section 3 we compute all the 3-point amplitudes with two off-shell legs, which are the fundamental starting point of the recursion and can be obtained via BCFW itself starting from the known on-shell results, as well as diagrammatically. Section 4 presents the explicit derivations of the 4 point amplitudes for 2 → 2 scattering in High Energy Factorization. Their BCFW calculation is worked out. For the cases with 4 fermions, when BCFW does not work, a diagrammatic evaluation is performed, so as to explicitly compute them all. Section 5 collects the all-leg results for MHV amplitudes. Finally, section 6 shows an explicit analytical study of the squared matrix elements and how these compare to the Hybrid Factorization case for the 4-gluon scattering. Some more technical details and cross checks are presented in the appendices.

The BCFW relation for amplitudes with two off-shell legs
We always consider scattering amplitudes with all particles outgoing or incoming and refer the reader to [34,35] for the details about our conventions as well as an introduction to the BCFW recursion. General informations about the notations are contained in appendix C.

Recursibility with two off-shell legs
In order to apply the BCFW recursion, one has to make the scattering amplitude A a function of the auxiliary complex variable A(z) in such a way that lim z→∞ A(z) = 0 . (2.1) We refer to [12,13] for the original proofs in the on-shell case respectively for only gluons and with fermions as well. In [35] a thorough discussion of allowed shifts for amplitudes with only one off-shell leg can be found. In this section we show that amplitudes with two off-shell legs and any number of fermion pairs are always BCFW-recursible, as long as there is at least one gluon among the external particles. By this we mean that they can always be expressed in terms of lower point amplitudes and of amplitudes of the same rank but with one less off-shell leg, according to the classification of residues provided in [34,35]. An argument leading to the corresponding conclusions for on-shell amplitudes has been around for some time [37]. It employs the background field method, previously used in [36] to give a physical motivation of the existence of the BCFW recursion relation. Its conclusion is that BCFW is viable for gauge theory and gravity amplitudes coupled to fermions and/or scalars as long as at least 1 gluon or graviton is an external leg [37]. Here it will be sufficient for us to argue in terms of individual diagrams.
In the gluon case, one can always shift the two off-shell legs and thus we get a factor 1/z 2 from the external lines, as in the case which is diagrammatically examined in [12], when the shifted gluons have helicities (h i , h j ) = (−, +). In fact, as off-shell gluons do not have corresponding polarisation vectors, we deliberately include their propagators in the amplitude in order to ensure proper behaviour for z → ∞. This implies the existence of additional poles, dubbed 'C' and 'D' terms, whose residues are just the amplitude with the leg exhibiting the pole being properly taken on-shell, as detailed in [34,35]. Another kind of poles featured only by off-shell amplitudes are due to vanishing denominators in eikonal propagators.
From the original argument in [12] for all-gluon amplitudes, we know that in individual diagrams there can be at most one more three-gluon vertex than there are propagators. This implies that, as long as we have a factor 1/z 2 from external legs, all diagrams behave at least as well as ∼ 1/z asymptotically. Also, from the observations above we understand that, as long as there is one gluon among the external particles, we can get a factor 1/z, whether it is on-shell or off-shell. The question is whether one can explicitly prove that this is enough in all cases. Fermions interact with the gauge field via a 3-point momentum-independent vertex and their propagators approach a constant for z → ∞. The key observation is that n-point amplitudes with k fermion pairs and m gluons (n = 2 k + m) can be thought of as n-point gluon amplitudes made by diagrams obtained by turning k couples of external gluons into fermions and changing accordingly vertices and propagators along the diagram, as depicted in the example in figure 1.
Clearly, we can do this switch only for the diagrams for which it is possible to make the fermion number flow along the diagram without meeting a 4-point vertex, but this simply means that the topologies of diagrams for n-gluon scattering for which this is not possible do not contribute to the amplitudes describing scattering of n − 2k gluons and 2k fermions. For instance, in the case of 4-point amplitudes, we can obtain the color ordered amplitudes for ggψψ from those for 4-gluon scattering with this prescription and, clearly, the diagram made by the 4-gluon vertex will not generate any contribution. This is illustrated in figure 2, JHEP02(2017)009 Figure 2. How to get the planar contributions to the scattering of 2-gluon and 2 fermions from the 2 → 2 gluon scattering.
where the three planar diagrams contributing to the 4-gluon amplitude are turned into the two diagrams representing a scattering process involving two fermions and two gluons. Now, for every possible Feynman diagram, the numbers of fermions and gluons, N f /g , are related to the number of propagators and vertices by the relations where the meaning of subscripts is obvious. Let us suppose that an n-gluon amplitude is turned into an amplitude with n − 2k gluons and k fermion pairs switching vertices and propagators according to P g → P f and V 3g → V f along the flow of the fermion number. Then we can express n in two ways for diagrams of the two amplitudes, where primes of course refer to the amplitude with fermions. As, obviously, V 4g = V 4g and where we have used (2.2) and ∆ is the difference between non-primed and primed quantities. The last equation means that the difference between gluon vertices and gluon propagators is bigger for the all-gluon amplitude than for the amplitude with 2k fermions. Consequently, the asymptotic behaviour of internal lines of the latter as a function of z will always be better, i.e. at least as good as O(z 0 ) As we can always shift one gluon leg which goes as 1/z for z → ∞, the desired asymptotic condition is fulfilled.

JHEP02(2017)009
We conclude by pointing out that eikonal propagators behave as gluon propagators for z → ∞, as the direction p does not shift for off-shell particles, whereas eikonal vertices do not depend on the momenta of the eikonal particles (see the Feynman rules in [35]), so that our arguments hold for diagrams containing such terms too.
Finally, in the event that one has to shift one gluon and one fermion line, it is always possible to do that without getting a z factor from the fermion spinor. For instance, it is known in the on-shell case that it is not allowed to shift one fermion and one gluon with the same helicity sign. This happens because an amplitude with, for instance, one fermion pair is given by If one shifts the fermion q − and one gluon with negative helicity g − , then the shift vector guaranteeing a 1/z factor from the gluon is e µ = g|γ µ |q]/2, for which the large z behaviours of the legs cancel out. On the other hand, g|γ µ |q]/2 is manifestly a good shift vector. A completely similar argument holds for A(q − , q + , g 1 , . . . , g n ) = [q| . . . |q . Then, when working with amplitudes with on-shell fermions, one can always choose the shift vector in such a way that the external fermion spinor does not shift. In [35] we called this the original LW prescription because it was first discussed in [13]. Finally, if the fermion is off-shell, then the spinor is not shifted anyway, so that our argument is complete. This concludes our proof.

Residues
The classification of the possible residues was already presented in [34,35]. The general formula which generalizes the BCFW recursion relation to the off-shell case is where the index s refers to the particle species, which can be gluon (g) or fermion (f ), and h refers to the helicity. We recall the nature of these contributions. The A s p,h residues are due to the poles which appear also in the original BCFW recursion for on-shell amplitudes, whereas the B s i terms are due to the poles appearing in the propagator of auxiliary eikonal quarks [25,26]. If the i-th particle is on-shell or off-shell but shifted, these terms are not present. The case when one off-shell leg is right next to a shifted leg leads to the appearance of B residues, which imply 2-point functions defined only for shifted momenta. Their proper treatment boils down to the the same final prescription, no matter whether the off-shell particle is a gluon or a fermion. In the case it is a gluon, this problem was already addressed in [34]. Here we collect the results in order to have them all available in the same place. Notice that the k-lines can represent an off-shell gluon as well as an off-shell fermion and the associated eikonal propagators. We explicitly distinguish the cases in which the i-th particle in the 2-point function is on-shell, as well as the shift vector; instead, the double j-line can be either on-shell or off-shell, but the resulting different shifts The bottom line is that, in these cases, in the (n−2)-point amplitude in the B-residue there is one more shifted leg, according to the prescriptions reported above. As the situation with the 2-point amplitudes is rather specific, in order to indicate them during the derivations reported in the following we use the symbol A 2 (K * ,ĵ). Finally, C s and D s denote the same kind of residue, which appears respectively when the shifted i-th or j-th particle is off-shell. They are due to the vanishing shifted momentum in the propagator of the off-shell particle: k 2 i (z) = 0 or k 2 j (z) = 0. These contributions are given by an amplitude with equal number of legs evaluated with the i-th or j-th leg onshell but with momentum spinors shifted according tables 1 and 2 times proper prefactors. For gluons the prefactors are shown in the figure, for fermions the only difference is x → √ x.
Let us recall that it was explicitly proven in [35] that, if the shift vector is chosen in such a way as to match the original LW prescription in the on-shell limit, then the C f and D f terms vanish identically. This was enough there, where only one particle was off-shell; now in principle both shifted particles (gluon and fermion) can be off-shell, so the statement is better recast in a slightly different way: the C f and D f residues are zero if the fermion spinor in the shift vector is the same as the one used in the diagrammatic construction of the amplitude. With our conventions (see also eq. (2.6)), the correspondence is This statement can be checked to be equivalent to the previous one through the explicit analysis presented in [35]. As we will see first in the section dedicated the to 3-point amplitudes, also if both fermions are off-shell it is necessary to completely distinguish the two spinor choices. This statement will also be explicitly illustrated with examples. This was not necessary with only one off-shell particle, because in that case the spinor of the off-shell fermion is completely determined by the spinor used for the on-shell one, whose choice is in turn dictated by (2.8).
Before presenting explicit results, one might ask whether the BCFW algorithm is necessarily the most efficient one when it comes to calculate multi-leg tree level amplitudes. In particular, it is known from [38,39] that, at least for all-gluon amplitudes, it numerically outperforms the Berends-Giele off-shell recursion [6] only for multiplicities n < 9. Besides, JHEP02(2017)009 (h g , hq) = ( * ∨ −, −) × Table 2. Shifted quantities needed for the evaluation of C f and D f residues in the case in which the antiquark is shifted. A completely similar table holds for the case in which q is shifted (q → q).
BCFW can be made more efficient by the use of hard-wired lower-point amplitudes, so that the recursion terminates with these instead of the 3-point amplitudes. There is no special reason to expect that fermion pairs should change much. As far as the sole computation of the amplitudes is concerned, such a limiting multiplicity would be relevant for processes with 7 particles final state or, maybe more realistically for the time being, for NNLO corrections to 5-jet production, which is the present benchmark multiplicity for NLO calculations achieved by the BLACK-HAT and NJET collaborations [40,41], if one assumes that the hurdles to 2-loop computations with final-state multiplicities higher than 2 is to be overcome any time soon. Therefore, we can safely state that BCFW bears a twofold advantage: it employs, at every step, even in this modified version, only gauge-invariant quantities and is numerically competitive with other methods.

JHEP02(2017)009 3 The 3-point amplitudes
All the 3-point amplitudes with 1 off-shell leg have been provided in the first appendix of [35]. There we showed their explicit diagrammatic construction, but also stressed the complete independence of the BCFW approach from any Feynman-diagram calculation; indeed, we also showed how off-shell 3-point amplitudes can be inferred via the BCFW procedure from their on-shell counterparts, which are in turn determined solely by symmetry requirements. The same holds for this case, with the difference that the different kinematics implies that we must take two residues into account, one for each off-shell leg. This means that each of the amplitudes below is the sum of two terms, obtained by going on-shell with respect to both shifted off-shell legs in turn. We do not provide any diagrammatic construction here but we just list the complete set of results.

Two off-shell gluons
The 3 gluon amplitudes with two off-shell legs were given for the first time in [34] and are [13] . (3.1)

Off-shell gluon and off-shell fermion
Concerning amplitudes with one off-shell gluon and one off-shell fermion, simple calculations give the results , .
We explicitly show how to derive the first of (3.2) through BCFW. Let e µ = 1|γ µ |q]/2 be the shift vector. Asq cannot have but a negative helicity, such a shift vector does not respect the LW prescription. Despite the fact that for an off-shell fermion it is still viable, the residue associated with the fermion external momentum going on-shell is not vanishing anymore.
From tables 1 and 2 and taking prefactors properly into account, we see that

JHEP02(2017)009
where the right side follows from the results derived in [35]. Replacing the hatted spinors according to the results for the fermion D-pole in the same tables, we finally get the first line of (3.2).

Two off-shell fermions
This is the most complicated case and it deserves further clarifications than the previous two. First we report the results and then we comment on the meaning of the notation. , .
The superscripts accompanying the off-shell fermions must be explained. A fermion pair in an amplitude with all outgoing or all incoming particles must feature opposite helicities and thus there are only two possibilities, namely where the spinors assignment is uniquely determined by the fermion helicity. If only one of the fermions in the pair is off-shell, then its helicity is uniquely fixed by its partner's, which in turn determines the corresponding spinors, as in (3.5). But if both fermions in the pair are off-shell, then one has to consider both options, The meaning of the superscripts should now be clear: for each off-shell pair and for every possible configuration of the other particles, there are two different off-shell amplitudes. These amplitudes are distinguished by the two possible different assignments of the fermion helicities.
In appendix D we derive the first amplitude in two ways, first by explicit diagrammatic computation and then via BCFW recursion, in order to fully illustrate the point for the interested reader.

The 4-point amplitudes
In this section we report the new computational results: we present the analytical expressions for the four-point amplitudes with two off-shell particles. In the cases with at least one external gluon, we employed the BCFW recursion and cross-checked with the numerical results obtained with the help of AVHLIB [42], which constitutes by itself a very non trivial cross-check. All expressions for the amplitudes, except the four-fermion JHEP02(2017)009 cases, have been implemented for numerical evaluation in amp4hef. 1 Finally, all the fully squared matrix elements, averaged (summed) over initial (final) state colors and spins, were successfully cross-checked with the results reported in [19].

Amplitudes with 4 gluons and no fermions
These amplitudes were essentially all given already in [34]. They are The first four amplitudes in the first two lines are Maximally Helicity Violating (MHV) and are given by the general formulas (5.1).
Concerning the other four amplitudes, the following two are calculated in [34]: . (4.1) From cyclic invariance it is apparent that A(1 * , 2 + , 3 * , 4 − ) can be obtained from the first amplitude in (4.1) just by the relabelings (1 ↔ 3, 2 ↔ 4). In order to get A(1 * , 2 − , 3 + , 4 * ) from the second, instead, it is necessary to assume that the color-order reversed relation holds, before performing the relabelings (1 ↔ 4, 2 ↔ 3). We checked explicitly that this is actually the case for this amplitude and in appendix B we recall the short and beautiful proof of (4.2) presented in [43] and show that it carries over directly to our amplitudes with off-shell legs as well.

2 off-shell gluons
There are two such amplitudes, namely A(2 * ,q ± , q ∓ , 1 * ). So we compute the first one explicitly in appendix D, the other one being analogous. Here we just report the final JHEP02(2017)009 results, .

1 off-shell gluon and 1 off-shell fermion
This is the most complicated case, since these amplitude also depend on the position of the off-shell gluon relative to the fermion pair, a situation already met in [35].
There are 8 independent such amplitudes, the other 8 being parity conjugated and thus given simply by their adjoints. The independent ones can conveniently be grouped into MHV and NMHV amplitudes, i.e.
The general formulas for MHV amplitudes were already given in section 5.
As for the NMHV ones, we have to compute 4 of them. We explicitly show the calculation of the first one in appendix D, providing only the results for the remaining three and for all the adjoint amplitudes.

2 off-shell fermions
The final set of amplitudes we can evaluate with BCFW recursion are the ones with two on-shell gluons and one off-shell fermion pair. There are 8 of them, 4 of the kind A(2 ± ,q * , q * , 1 ± ) and 4 of the kind A(2 ± ,q * , q * , 1 ∓ ), because the fermion pair, as detailed in section 3, can come in two ways, depending on the way the spinors are assigned. The formulas for the MHV amplitudes A(2 ± ,q * , q * , 1 ± ) were given in section 5. The 4 non MHV amplitudes require some work, which is explicitly done for two of them which are non parity-conjugated in the appendix D.
The final explicit results are  Figure 3. The three contributions to A(ū * + , u * − ,d + , d − ) and A(ū * + , u * − ,d − , d + ). In order to compute this amplitudes, we can no longer apply the BCFW recursion, as the condition of having at least one external gauge boson is not fulfilled. So, as we want to provide a full set of 4-point amplitudes for high-energy factorization, we have to compute them via Feynman diagrams expansion. The case with two quark pairs of different flavours is simpler and we will start with it. Let us use the letters 'u' and 'd' for different flavours. Modulo parity conjugation, the independent amplitudes are We start with the case with an off-shell quark-anti-quark pair of the same flavour, which is diagrammatically represented in figure 3. Using the Feynman rules, one can calculate quite easily that Next is the configuration with two off-shell antiquarks, which is made by the four contributions depicted in figure 3. Applying the aforementioned Feynman rules to the four contributions shown in figure 4, we arrive at (4.9)

Amplitudes with 2 equal-flavour fermion pairs
We come to the last two sets of amplitudes, which we need to describe the processes q * q * → qq andq * q * →q q.

JHEP02(2017)009
Concerning the first process, we have the four independent contributions , , As we are working in an all-incoming setup, the letter u denotes the off-shell quarks and u stands for the on-shell antiquarks, which become outgoing quarks when switching their momenta signs. The particles with the label 1 together make the t-channel.
For the second process instead, the four independent contributions are , . (4.11) Here the labels 1 and 2 stand respectively for the off-shell and on-shell quark-antiquark pair, whereas the t-channel is made by the momenta of the two (anti-) quarks.

JHEP02(2017)009
The interested reader can reproduce these results using the Feynman rules and applying them to Feynman diagrams which are similar to the ones shown in the previous section. Specifically, the amplitudes A(u * + 1 ,ū − 1 ,ū − 2 , u * + 2 ) A(u * + 1 ,ū − 1 ,ū + 2 , u * − 2 ) are schannel processes and are thus represented by diagrams similar to the set in figure 3.

The four amplitudes
are u-channel scatterings. Each of these later six processes is made of four contributions analogous to those depicted in figure 4, with the only different that the gluon exchange must happen in the corresponding t or u channel instead than in the s one.

The general formulas for MHV amplitudes
In this section we derive formulas for amplitudes which are Maximally Helicity Violating (MHV). First of all we recall the MHV formulas derived in [34] for the all-gluon case, [21] .
Then we set up to find out the general structure of the MHV amplitudes in the case we have 1 off-shell gluon and 1 off-shell fermion and two off-shell fermions.

Squared 2 → 2 matrix elements
In this section we sketch the kinematical properties which are specific to matrix elements in HEF and may be relevant for further phenomenological investigations. We also square our matrix elements by averaging (summing) over the initial (final) state colours and spins, explaining how to express them in terms of the rapidity difference ∆y and the azimuthal angle difference ∆φ between the final state partons. Only for this section, we adopt the (physical) convention that the two off-shell particles are incoming and the other two are outgoing. This means that momentum conservation reads In doing this, we generalise the discussion of the collinear case presented in [44] Working in light-cone variables, the q ± momenta for a general 4-vector q are defined as q ± = q 0 ± q 1 , (6.2) and x 1 is the beam direction. The rapidity variable is defined as

JHEP02(2017)009
which naturally allows to express the particle energy and longitudinal momentum of massless particles as The incoming momenta k a/b and the outgoing momenta p a/b are thus found to be respectively p a = p a⊥ e ya , p a⊥ e −ya ; p a⊥ = p a⊥ e ya , e −ya ; cos φ a , sin φ a Here, differently from the collinear case, k a/b⊥ = 0, which implies imbalance for the final state transverse momenta, p a⊥ = − p b⊥ ; k a/b⊥ ≡ | k a/b⊥ |, p a/b⊥ ≡ | p a/b⊥ |, S is the proton center of mass energy, y a/b are the rapidities of the final states particles and the azimuthal angle difference is defined as the angle between their directions ∆φ = arccos (p a⊥ ·p b⊥ ), whereas θ a is the azimuthal angle for the transverse momentum of the first incoming offshell particle. By projecting momentum conservation (6.2) onto the light-cone components, one easily gets the further constraints By using (6.6), we finally end up with the parameterizations k a = p a⊥ e ya + p b⊥ e ya+∆y , 0; k a⊥ cos θ a , k a⊥ sin θ a , p a = p a⊥ e ya , e −ya ; cos φ a , sin φ a , p b = p b⊥ e ya+∆y , e −ya−∆y ; cos(φ a + ∆φ), sin(φ a + ∆φ) , after introducing the rapidity difference ∆y. We are thus left with 8 independent parameters on which the squared matrix elements depend: k a⊥ , θ a , p a⊥ , p b⊥ , y a , φ a , ∆y, ∆φ. At the end of this paper, we present some plots of squared matrix elements as functions of (∆y, ∆Φ) for the following given values of the first 6 variables, k a⊥ = p a⊥ = p b⊥ = 20 GeV, θ a = 0.2, y a = −2, φ a = 2.9. We choose to show the processes g * g * → gg, g * g * → qq and q q → qq plotted in figures 7-8. In figures 5-6, instead, we have inserted the plots of the squared amplitudes for 2-gluon scattering with two and one on-shell legs respectively. As momentum conservation in the collinear case enforces φ = π, the second variable is chosen to be the final state transverse momentum of each jet, so that the plot of an actual process is to be thought as a one dimensional section perpendicular to the k T axis. It is nevertheless useful for illustrative purposes.  Finally, in figures 9 and 10 we show explicitly the ratio of the squared matrix elements for g * g * → gg and g * g → gg for the aforementioned as well as one more set of input values, which differ from the previous ones only for the value of φ a = 0.5; this is done to stress the non trivial difference between the two configurations and the sensitivity of the difference between the two cases to initial conditions. In all the plots the strong coupling constant is set α s = 0.2, so as to facilitate the comparison with the plots in [45]. The partial amplitudes making up the matrix element for g * g → gg, which are all MHV, can be found in [34,35].

JHEP02(2017)009
For the sake of clarity, we will limit our discussion to the only-gluon case. What is striking in figures 6 and 7 with respect to figure 5 is the possibility of having the configuration of final state jets away from back-to-back configuration. These configurations are precisely what High Energy Factorization is supposed to produce, thanks to the transverse JHEP02(2017)009  momentum imbalance due to the initial state partons non vanishing k T . In the collinear approach, such a dispersion for 2 → 2 LO scattering process is generated at the cross section level by initial state parton shower effects. The spike in the center, which diverges for vanishing ∆φ and ∆y, corresponds to the physical situation in which the two jets originating from the final state partons are unresolved and become one single jet. As for the comparison between the hybrid factorization case (only one off-shell leg) and the full High Energy Factorization one (two off-shell legs), we limit ourselves to showing three plots -figures 9-11 -featuring the ratios of only-gluon squared matrix elements on a coinciding subset of independent variables (i.e. the final state transverse momenta), while keeping the other variables (θ 1 , φ 1 , y 1 and p T 1 ) fixed. As we scan three possible differences of final state transverse momenta, in ascending order, we clearly observe that JHEP02(2017)009   the rapidity gradient in the ratio plots increases. A more complete scan of the various possible kinematic configurations would require too much space for a paper like the present one, which is concerned mainly with theoretical aspects of the problem, but remains an interesting problem to be addressed separately.

Summary and perspectives
We have studied the problem of calculating amplitudes with 2 off-shell legs for High Energy Factorization via BCFW recursion. We found a full generalisation of the recursibility criteria devised in [34,35], proving independently that, exactly as in the fully on-shell case, amplitudes with at least one gluon leg can be recursed, whereas, when there are only external fermions, one still has to rely on the ordinary Feynman diagram expansion. Other types of generalisations of BCFW are underway in the on-shell case [46][47][48][49], which amount JHEP02(2017)009 to devise ways to pin down the boundary contributions which show up when the asymptotic condition on the complexified amplitude is not satisfied. The off-shell implementation of such methods is, however, well beyond the scope of the present work and we will not discuss it further. For the first time, a complete set of exact analytical expressions of color ordered 4point amplitudes with 2 off-shell legs for a Yang-Mills theory with massless fermions has been provided, together with MHV formulas for any number of legs. Also a discussion of their properties and differences w.r.t. the collinear case is proposed. All expressions for the amplitudes, except the cases with two quark-antiquark pairs, have been implemented in amp4hef [42] for numerical evaluation.
Concerning the potential future developments, it would of course be good to push the program to higher multiplicities, especially now that the theoretical framework to perform such calculations has been completely defined. Considering the present applicability limits mentioned at the end of section 2.2, it would be nice to arrive at expressions for 6-or 7-point amplitudes. Nevertheless, performing such calculations by hand for every case can be quite demanding, so the development of a dedicates software package is envisaged to this goal.

A The 3-point amplitudes with less off-shell legs
In order to make our paper self-consistent, here we list the 3-point amplitudes used in the construction of the 4-point amplitudes presented in this paper. They are either on-shell or have only 1 off-shell leg.

JHEP02(2017)009
where the inductive hypothesis is used between the first and second line and which completes the proof. For gluon amplitudes with up to 2 off-shell legs, one can always choose to shift the off-shell gluons in order to avoid B-poles, so that the only difference w.r.t. the on-shell case is the presence of additional n-point amplitudes with one less off-shell leg, because of the C and D residues; the latter amplitudes can be BCFW recursed themselves until all off-shell legs disappear and the color-order reversed relation holds because of (B.1). For the A-terms, instead, the procedure in (B.1) works fine, thus completing the sought generalization.

C.1 Momenta
We repeat here the parameterisation of the momenta used in [34,35], which reads k µ 1 + k µ 2 + · · · + k µ n = 0 momentum conservation (C.1) where for every k µ i there is a corresponding orthogonal, light-like direction p µ i . With the help of an auxiliary light-like four-vector q µ , the momentum k µ can be decomposed in terms of its light-like direction p µ , satisfying p · k = 0, and a transversal part, following If the momentum k is on-shell, then it coincides with its associated direction, p = k. The coefficients κ and κ * do not depend on the auxiliary momentum q µ (see appendix C.3) and We focus on the fundamental color-ordered or dual amplitudes, where the gauge-group factors have been stripped off and all particles are massless. By construction, these amplitudes contain only planar Feynman graphs and are constructed with color-stripped Feynman rules.
The polarisation vectors for gluons can be expressed as where q is the auxiliary light-like vector and g stands for the gluon momentum. In the expressions of the amplitudes, gluons will be denoted by the number of the corresponding JHEP02(2017)009 particle, whereas we will always explicitly distinguish quarks and antiquarks with q andq respectively. Finally, the polarization vectors for the auxiliary photons coming in pairs with the off-shell quarks are where f denotes the fermion momentum spinor and q is again the auxiliary vector. When shifting momenta with an auxiliary complex variables, as required by the BCFW procedure, we pick up two particles, say i and j, and we require their light-like directions to be the reference vector of each other, so that their momenta -with a generally non vanishing transverse component -become (C.9) In our notation |p a ](|p a ) ≡ |a](|a ), where g a is any gluon, whereas for quarks |p q ](|p q ) ≡ |q](|q ) and |pq](|pq ) ≡ |q](|q ). Just as in [34], we choose the shift vector to be which satisfies p i · e = p j · e = e · e = 0 . (C.11) Shifted quantities are denoted by a hat symbol. The most general case shifted momenta are thusK Momentum conservation and the eikonal conditions p i ·K i = 0 and p j ·K j = 0 are manifestly preserved. So, for an off-shell momentum, only one of the scalar coefficients (κ or κ * ) in the transverse momentum shifts, whereas the directions do not shift. We report the shifts induced by i|γ µ |j]/2, i on-shell: |î] = |i] + z|j] j off-shell:κ * i = κ * j + z ij j on-shell: |ĵ = |j − z|i (C. 13) JHEP02(2017)009

C.2 Spinors
The left-and right-handed spinors used in this paper can be defined in the following way: The "dual" spinors are defined as The definition of the "dual" spinors does not involve complex conjugation and all spinors are well defined for complex momenta.
Defining the Pauli vector as our conventions for the Dirac matrices are so that γ 5 is given by The dyadic product satisfies the identity The following spinor identities hold and are extensively used in the paper, with the convention that p µ , q µ and r µ are always light like, whereas k µ does not have to be light-like  Figure 12. Feynman diagrams contributing to A(1 + ,q * + , q * − ).
The last two relations trivially reduce to 0 = 0 for odd n, and the one before reduces to 0 = 0 for even n.

C.3 Schouten identity
A very insightful way to express the Schouten identity is to write it as a completeness relation. For any p µ , r µ with p 2 = r 2 = 0 and p·r = 0 we have where we used the fact that p|K /|p] = 2p·k = 0. The same can be shown for κ * .

D Some explicit calculations
This appendix is devoted to preset some explicit BCFW derivations which were omitted in the text for the sake of shortness.
The first calculation we present is the evaluation of the first of the three-point amplitudes in (3.4). We do this in two ways, diagrammatically (see figure 12) and via BCFW, in order to provide an explicit cross-check. There are three Feynman diagrams contributing, shown in figure 12, where all momenta are taken as incoming. In the first one the two off-shell fermions are coming from the QED vertex between the auxiliary fermion and photon, in the remaining two auxiliary particles pair are bent open in turn and auxiliary fermions appear in internal lines as eikonal propagators, denoted with a dashed line according to the Feynman rules [35]. Notice that, as pointed out in [26], auxiliary photons interact only with the corresponding auxiliary fermion, not with any other charged fermion.  Figure 13. The only BCFW residue for A(1 + ,q * + , q * − ) with shift vector e µ = q|γ µ |1]/2.
The implicit expression for our amplitude that one gets from the Feynman diagrams is thus (D.1) where the helicity of the auxiliary photons is enforced by their position, because (see eq. (C.8)) the other one would give just zero. We choose the fermion momenta p q and pq to be auxiliary vectors for each other's auxiliary photon. One degree of freedom still at our disposal in the computation is the choice of the auxiliary vector for the polarization of the gluon. From eq. (C.7) one can check that choosing p q annihilates the first and third contribution, so that we are left with after exploiting the definitions of κ q and κ * q and momentum conservation. Now we compute the same amplitude by BCFW recursion using q|γ µ |1]/2 as shift vector, which respects the LW prescription in the on-shell limit, so that we only have a B-pole, due to the eikonal propagator of the auxiliary eikonal fermion used forq going onshell, as depicted in figure 13. The following relations are derived by requiring the location of the pole to be pq ·K q = 0, In particular, they implyˆ µ 1 + = µ 1 + , giving On the left of the unshifted eikonal propagator is the 2-point amplitude of the gluon and theq fermion current, whereas on its right is the product of the two auxiliary photon vertices with the q propagator ( evaluated for shifted momentum) in the middle. One can visualise these "unconventional" 2-point amplitudes thinking about them as the parts of the second Fenman diagram in figure 13 on the left and right side of the propagator and evaluated for shifted momenta.