A Combinatoric Shortcut to Evaluate CHY-forms

In \cite{Chen:2016fgi} we proposed a differential operator for the evaluation of the multi-dimensional residues on isolated (zero-dimensional) poles.In this paper we discuss some new insight on evaluating the (generalized) Cachazo-He-Yuan (CHY) forms of the scattering amplitudes using this differential operator. We introduce a tableau representation for the coefficients appearing in the proposed differential operator. Combining the tableaux with the polynomial forms of the scattering equations, the evaluation of the generalized CHY form becomes a simple combinatoric problem. It is thus possible to obtain the coefficients arising in the differential operator in a straightforward way. We present the procedure for a complete solution of the $n$-gon amplitudes at one-loop level in a generalized CHY form. We also apply our method to fully evaluate the one-loop five-point amplitude in the maximally supersymmetric Yang-Mills theory; the final result is identical to the one obtained by Q-Cut.


Introduction
Scattering equations are derived at tree level for the high-energy behavior of string theory in [2,3] and have drawn attention from theoretical physicists in diversified contexts [4,5]. Incorporating scattering equations, Cachazo, He and Yuan [6][7][8] propose a closed formula for arbitrary n-point tree amplitudes in a variety of massless quantum field theories. This form is proven in [9] for Yang-Mills theory in arbitrary dimensions and a polynomial form of the scattering equations is obtained by the same group in [10].
The scattering equations are generalized to loop levels in a number of contexts, such as open string theory, pure spinor formalism of superstring, and ambitwistor string theory [11][12][13]. The CHY expressions are subsequently extended to one and two loops for the bi-adjoint scalar theory, gauge theory and gravity in [14][15][16][17]. In addition to the ambitwitor approach, the one-loop generalized CHY forms are obtained from tree-level ones in one higher spatial dimension in [18,19] for scalar and gauge theories, and higher-loop CHY forms for scalar theory are constructed in [20].
Besides constructing the generalized CHY forms, it is also a challenge to evaluate such multi-dimensional CHY integrals. Indirect methods, such as "building block method" [21][22][23][24] and "Integration rules" [25][26][27][28][29][30][31][32], evaluate the generalized CHY integrals reductively. Direct approaches are also explored. At tree level, the n-point CHY form is evaluated for gauge and gravity theories in special kinematics in [33] and the scattering equations are solved in four dimensions up to six points in [34][35][36]. Elimination theory is exploited in solving the scattering equations up to seven points in [37,38] and a general prescription based on elimination theory is proposed in [39]. A direct evaluation of the CHY form for the MHV tree amplitude is given in [40]. Algebraic geometry based methods-the companion matrix method [41], Bezoutian matrix method [42], and polynomial reduction techniques [43,44]-are employed to evaluate the CHY expressions, without solving the underlying scattering equations.
In our previous work [1] we proposed a conjecture that enables us to compute multidimensional residues on isolated poles (0-dimension) by a differential operator.
Here we briefly recall that conjecture.
Suppose f 1 , f 2 , ..., f k are homogeneous polynomials in complex variables z 1 , z 2 , ... z k of degrees d 1 , d 2 , ..., d k respectively. And we assume that the common zeros of f 1 , f 2 , ..., f k consists of a single isolated point p. Let R(z i ) be a holomorphic function in a neighborhood of p. Then the conjecture states that the residue of R at p can be computed by a differential operator D as follows, where D is of the form D = {r i } a r 1 ,r 2 ,··· ,r k (∂ 1 ) r 1 (∂ 2 ) r 2 · · · (∂ k ) r k . (1.2) Here coefficients a r 1 ,r 2 ,...,r k are z-independent constants and ∂ i = ∂ ∂z i , i = 1, . . . , k. The sum is done over all solutions of the equation k i=1 r i = k h=1 d h − k. Moreover, it is conjectured that D is uniquely determined by two conditions respectively from 1) the local duality theorem [45,46] and 2) the intersection number of the divisors D i = (f i ).
This conjecture is verified numerically to be widely applicable for computing both degenerate and non-degenerate multi-dimensional residues, as long as the poles are isolated. In the application to evaluating the generalized CHY integrals, the linear equations for a r 1 ,··· ,r k arising from the local duality theorem and the intersection number requirement are the structures of interest. In this paper, we study these structures by introducing a tableau representation for each a r 1 ,··· ,r k . Using the tableau representation, we explain how to obtain the differential operator efficiently, especially for the polynomial scattering equations, which makes the evaluation of the generalized CHY forms a simple combinatoric problem. In this way, the evaluation of CHY integrals in their prepared forms (which we define later in this paper) is easily achieved for any number of external lines at one loop.
This leads straightforwardly to a full evaluation of the n-gon one-loop integrand in the generalized CHY form. Furthermore, many a term in the integrands of super-Yang-Mills expressed in the generalized CHY form can also be re-casted into the prepared form, which are then calculated effortlessly. We apply our method on the one-loop five-point amplitude in SYM and find that the long expression in the Pfaffian can be transformed into the prepared form by cross-ratio identities [32]. Upon evaluating the prepared forms, the final result is identical to the ones obtained in [47] by Q-Cuts [48].

Prescription for determining differential operators
In this section we introduce our new method to determine the differential operators that is more efficient than the approach taken in [1]. We begin with a warmup example for the four point one loop SYM integrand. After that we present the notations used in obtaining the coefficients in the differential operators. Finally we use this method to get the complete solutions of the coefficients for a class of CHY integrands.

Toy model: one-loop four-point SYM integrand
The one-loop scattering amplitude for four external gluons in N = 4 SYM has been revisited in various contexts in literature [47,49]. The generalized CHY integral for its loop integrand is given in [15,16]. A detailed analysis of the evaluation of the generalized CHY integral has already been presented in [1]. In this section we reconsider this example as a toy model to illustrate the main ideas of the method proposed in this paper. Details of and more general discussions on the method are to appear in later sections.
Firstly let us recall the settings. The four-point generalized CHY [15,16] form of the integrand reads, where P T 4 is the well-known Parke-Taylor factor that reads .
where S 4 is the cyclic permutation group of four objects. For illustrative purpose, here we only consider the first term in P T 4 which leads to the following integral where we fix the gauge σ 4 = 1 and the polynomials h i 's are and we have used the following notations, As shown in [1], using the global residue theorem, the integral can be rewritten as the following whereh i is obtained from h i with the replacement σ 4 → σ ′ 4 . From now on in this section we will simply write σ 4 for σ ′ 4 for notational convenience. The main idea of the method presented in [1] is to replace the contour integration with a differential operator D, which should be of the third-order in this particular case and can be written as D = 0≤r i ≤3 , r 1 +r 2 +r 3 +r 4 =3 a r 1 ,r 2 ,r 3 ,r 4 ∂ ∂σ 1 where the coefficients a's are to be fixed by the local duality theorem and the intersection number equation. One of the main points of the current paper is to provide an efficient method to determine those coefficients without solving the equations by brute force.
Now we briefly describe how the shortcut method works out in this particular example. There are 20 coefficients a r 1 ,r 2 ,r 3 ,r 4 to be determined. Before turning to the equations from local duality and intersection number, we firstly classify the coefficients as follows. For each coefficient a we define its rank to be R(a r 1 ,r 2 ,r 3 ,r 4 ) = r 1 !r 2 !r 3 !r 4 !. Furthermore, we associate a tableau with each a. All tableaux (up to subscript permutation) for our example are shown in Fig 1. The tableau in Fig 1(a) (including the permutations of the red tiles) denotes a 0,1,1,1 and all those with permuted subscripts. In Fig 1(b), the corresponding variables are a 0,0,1,2 and all those with permuted subsripts. In Fig 1(c), the corresponding variables are a 0,0,0,3 and all the permutations ρ of the index {0, 0, 0, 3}. With the definition of the rank of coefficients and the tableaux representation, solving the equations of the a's can be carried out in a very specific way that is much more efficient than treating it as a general equation solving problem. Recall that the local duality theorem yields, We will see that the equation D (σ i σ j σ 1 ) = 0 simply leads to the vanishing of all a r 1 ,_,_,_ with r 1 ≥ 1. So this means, for a non-zero a, the number of tiles in the first row in the corresponding tableau is zero. On the other hand, the equation Dh 3 = 0 means the vanishing of a's of rank 1. Thus a's corresponding to the tableau in Fig. 1 (a) are all vanishing. The equation D σ i σ ih1 will give us a relation between a's of rank 6 and a's of rank 2, e.g a 0,0,0,3 = −2! 3!ℓ 4 (ℓ 1 a 1,0,0,2 + ℓ 2 a 0,1,0,2 + ℓ 3 a 0,0,1,2 ) .
This equation can actually be read off directly from the tableau in Fig. 1 (c). By permuting the red tiles of a 0,0,0,3 's tableau among the rows, we could get the tableaux corresponding to a 1,0,0,2 , a 0,1,0,2 and a 0,0,1,2 respectively. In terms of the coefficients, this exactly corresponds to (2.7). Namely, a 0,0,0,3 multiplied by its rank and l 4 equals the minus sum of a 1,0,0,2 , a 0,1,0,2 and a 0,0,1,2 multiplied by their ranks and l i , where i indicates the row index of the red tile of the corresponding tableau.
The remaining unknown a's are those corresponding to the tableau Fig. 1 (b) with rank 2. They are all the variables present in the intersection number equation (2.8) Other equations involving variables of rank two are D σ i σ jh1 with i = j and D σ ih2 .
For example D σ 3 σ 4h1 = 0 is 2ℓ 4 a 0,0,1,2 + 2ℓ 3 a 0,0,2,1 = 0, (2.9) where we have used the fact that a 0,1,1,1 and a 1,0,1,1 are zero. This equation is represented by Fig. 2  Similarly this equation can also be read off from the tableau. For that purpose, we paint red the two tiles as shown in Fig. 2(b) and permute the red tiles among the four rows. To get Eq. (2.10), we simply take the sum of the products of the corresponding a's with their ranks and corresponding l i,j 's, where i, j are the row indices of the red tiles of respective tableaux. Accordingly, it is also seen that the first term and the

Preliminary: notations & tableaux
In the previous discussion, the one-loop SYM integrand for four external particles is massaged and turned into a sum of several terms, all of which are of the following form, where σ j (j = 1, · · · , n − 1) are the (n − 1) variables that need to be integrated out (with the gauge choice σ n = 1) and σ ′ n is the auxiliary parameter that comes into play in the homogenization of the scattering equations (for details, see Section 4.2 and 4.3 of [1]). h j (j = 1, · · · , n − 1) denote the (n − 1) polynomial scattering equations (3.6) that are necessary to capture the behavior of an n-point scattering and h ′ n = σ i (i = 1, · · · , n − 1) or h ′ n = σ ′ n . In this form, the divisors are generated either by one of the polynomial scattering equations or by the simplest monomial possible. As reviewed in the introduction, such contour integrals are associated with a differential operator with a certain number of coefficients fixed by the local duality theorem and the intersection number. These constraints are all linear and especially the local duality theorem gives rise a sparse constraint matrix, which in general can be solved by various sparse matrix methods. In the case of an integral in the form of (2.12), thanks to the beautiful mathematical structures of scattering equations, these constraints can indeed be solved analytically. In this section we present the method that reduces considerably the size of the constraint matrix and obtains quickly the analytic solutions of the coefficients. For this reason, we would like to call an expression of form (2.12) prepared. Moreover, as we will observe in more examples, the CHY-like representations for amplitudes/integrands can often be reduced to the prepared form or similar ones with a slightly modified h ′ n , even though the original seem much more complicated. For an expression with a modified h ′ n which still remains a monomial, we expect our algorithm can be generalized.
As shown in the discussion of the four-point one-loop integrand, to demonstrate our method in a more intuitive way, a graphical representation of the coefficients in the operator comes in handy. Let us take a coefficient a 0,1,2,3,4 as an example and address a few details of its corresponding tableau depicted in Fig.3  tableau corresponds to the i-th index of the coefficient while the number of the tiles in the i-th row is equal to the value of the i-th index. The total number of tiles is the order of the corresponding differential operator. As mentioned above, we define the rank of as coefficient a r 1 ,r 2 ,··· ,rn as the following integer number, R(a r 1 ,r 2 ,··· ,rn ) = r 1 ! r 2 ! · · · r n ! , (2.13) and this number is written on top of the corresponding tableau. Sometimes we drop the row labels in the tableau as shown in Figure.3(b), and this tableau represents the class of a's whose indices are related by permutations. For instance, the tableau in Figure.3(b) is corresponding to the coefficients {a ρ(0,0,1,2,3) } ρ∈S 5 where S 5 is the permutation group of five items and ρ(0, 0, 1, 2, 3) means a permuation of the five digits (0, 0, 1, 2, 3).
As shown in Figure.3, we have painted some tiles red in the tableaux. For a given unpainted tableau, we associate with it a class of colored tableaux which are obtained in the following way: in each row only the rightmost tile can be painted red, and we can choose to paint it or not. So altogether we have 2 J colored tableaux for each given unpainted one, where J is the number of nonempty rows in the unpainted tableau. From each colored tableau, we can get several new tableaux (without considering their coloring) by permuting the red tiles among the rows. For example, in Figure.4 are depicted the new tableaux resulted from the tile moving of the first colored tableau. In this example, we can obtain a 0,0,1,2,3 , a 0,0,2,1,3 , a 0,1,1,1,3 and a 1,0,1,1,3 from a 0,0,1,1,4 . So far we have played with the tableaux. Now we demonstrate how the tile moving can acquire actual meanings from the local duality theorem. Recall that the polynomial scattering equation of degree m for n particles takes the following form (throughout this paper we adopt the gauge choice σ n = 1 and all the polynomial scattering equations, if not specified otherwise, are homogenized with an auxiliary parameter σ ′ n ; for notational compactness, we drop the prime and simply denote the auxiliary as σ n from now on), (2.14) The order of the differential operator D associated with the prepared form for n points is Since the differential operator D acting on a monomial in σ i 's does nothing but picking out the exponents, the local duality theorem yields, where r j 's are non-negative integers satisfying the Frobenius equation r 1 + r 2 + · · · r n = M − m. Here we see that the rank of a r 1 ,··· ,rn is merely the factorial product that arises when taking derivatives multiple times. There are C m n coefficients a r 1 ,··· ,rn 's on the righthand side above and their corresponding tableaux all have r j black tiles in the row j and m red tiles that sit in the m different rows {i 1 , i 2 , · · · , i m }. The tableaux corresponding to these a's are exactly those related by permuting the red tiles as mentioned above! In other words, given a tableau that has r j black tiles in the j-th row and the rows {i 1 , · · · , i m } have one red tile each, there is a unique equation following from the local duality theorem that relates this tableau and those obtained by permuting the red tiles. Furthermore, scanning over the monomials of the correct degrees to construct the local duality equations is equivalent to exhausting all the legitimate colorings and writing down the corresponding equations.
As mentioned above, the equations arising from the local duality theorem are not independent from each other and the number of equations is much larger than the number of coefficients. With the correspondence between the colorings and the local duality constraints established, we are ready to present a method that systematically selects just enough independent local duality equations to fix the coefficients completely.

Coefficient generating algorithm
In this section, we spell out our method that analytically solves the local duality and the intersection number constraints associated with a prepared form and display the general solutions of the coefficients in the corresponding differential operator.
Although there are several possible colored tableaux associated with a given uncolored tableau, there is a preferred one that leads (by shuffling the red tiles among the rows) to new tableaux whose ranks are either lower than or equal to that of the original. This particular colored tableau is constructed as follows. For convenience let us assume a r 1 ,r 2 ,··· ,rn is a coefficient with indices such that r i 1 r i 2 · · · r in . In the prefered colored tableau, we paint red the last tile in row i n (the top row). Now for all subsequent rows, if the row i p has one red tile and r ip − r i p−1 1, r i p−1 > 0, we should also paint red the last tile of the row r i p−1 ; otherwise all the tiles in the row r i p−1 remain black and the coloring process stops here.
There are two types of the coefficients a r 1 ,r 2 ,··· ,rn : let (j 1 , j 2 , · · · j n ) be a particular permutation of the indices such that r j 1 r j 2 · · · r jn ; if r j+1 − r j 1 for j = 1, · · · n − 1, a r 1 r 2 ···rn is called elementary; otherwise it is called non-elementary. It is easy to observe that, for an elementary coefficient, the aforementioned preferred coloring relates it with those of lower or equal ranks. A non-elementary coefficient gets turned into only lower-rank ones through the preferred coloring. Moreover, for a particular prepared form for n points, the corresponding tableaux are all of M = (n − 1)(n − 2)/2 tiles in total. Up to row permutation, there is a unique tableau with the highest rank among the elementary ones, that is the tableau associated with a 0,0,1,2,3,...,n−2 whose rank is R e = (n − 2)!(n − 1)! · · · 1!.
The above observations about the prepared form actually provide us a natural solving order in which the local duality constraints can be conveniently solved. In the following discussion, we will show that in the operator associated with a prepared form, any coefficients with ranks lower than R e must vanish while the elementary ones with rank R e can be easily worked out, which leaves us only non-elementary ones with ranks greater than R e . Those higher rank coefficients can however be obtained by an inductive method: suppose we have obtained all the coefficients -both elementary and non-elementary -up to some rank R 0 R e and the next possible rank is R 1 , then a non-elementary coefficient of rank R 1 whose tableau is painted in the preferred way can be directly read off from the relation associated with the colored tableau, i.e. it can be represented in terms of other coefficients of lower ranks only. This procedure continues all the way up to the highest-rank coefficients and thus in this way we could determine all the coefficients analytically one after another.
Trivial coefficients In general, the differential operator D may be very complicated. However, it is easy to see that many coefficients in D are actually zero, i.e trivial, when D is associated with a form that is prepared (see (2.12)).
Firstly, for a prepared form with h n = σ i , the coefficients a r 1 ,··· ,rn with r i 1 must vanish, because the local duality theorem for h n yields, where s j are non-negative solutions to r 1 + · · · + r n = M − 1.
Secondly, the coefficients a r 1 ,··· ,rn must vanish if R(a r 1 ,··· ,rn ) < R e . To see this, let us consider a prepared form with h n = σ i . Since r 1 + · · · + r n = M, the lowest-rank coefficients must be elementary and have r i 1 when n 5 (when n = 4 the lowestrank coefficient has the i-th index being 0 and all the rest 1), which obviously vanish. Now we only need to show that every elementary coefficient of the form a r 1 ,··· ,r i−1 ,0,r i+1 ···rn (where r 1 ! · · · r n ! < R e ) is zero, since the non-elementary ones can be written in terms of coefficients with lower ranks only. Let {i 1 , · · · , i i−1 , i i+1 , · · · , i n } be a permutation of {1, · · · , i − 1, i + 1, · · · , n} such that r i 1 · · · r in and r i j+1 − r i j 1, r i+1 − r i−1 1. Then we must have r in n − 3 in order to guarantee that the rank is lower than R e . 1 On the other hand, we must also have r i 1 = 1 because r 2 + · · · + r n = M = ord(D). 2 Hence the tableau for a r 1 ,··· ,r i−1 ,0,r i+1 ,··· ,rn has only the i-th row empty and all the other rows filled. The aforementioned preferred coloring leads to (n − 1) red tiles in total and the local duality relation represented by this colored tableau involves a r 1 ,··· ,r i 1 ,0,r i+1 ,··· ,rn and other coefficients that all have r i = 1, i.e. a r 1 ,··· ,r i−1 ,0,r i+1 ,··· ,rn is rewritten as a linear combination of several vanishing a's and therefore is zero. By induction, all the non-elementary coefficients with ranks lower than R e are also zero. 1 If r in = n − 2 and all the other r j 's saturate their lower bounds, we have 2 Likewise, if r 1 = 0 and all the other r j 's saturate their upper bounds, we have r 2 + · · · + r n = r i2 + · · · + r in = 0 + 1 + · · · + (n − 3) + (n − 3) < (n − 1)(n − 2) 2 = M .
Therefore we must have r i1 = 1 and r in n − 3.
Non-trivial elementary coefficients Now we turn to the elementary coefficients of rank R e . Recall that the intersection number constraint reads where h n = σ i and the scattering equations h j (j = 1, · · · n − 1) are polynomials in σ i 's of degree j respectively. Using the explicit forms of h's, we can see that only coefficients of ranks no greater than R e appear in (2.17).

(2.23)
Now we consider another two permutations, where the right-hand sides are related by two permutation actions and the corresponding equation is given by the local duality theorem for h 2 which reads, where m j is defined the same way as above and for this case we have m v (ρ 2 ) 1,2,··· ,n−3,n−2 ℓ 1,··· ,i−1,i+1··· ,n .
(2.28) Summary of the coefficient generating method Now let us summarize our method that computes the coefficients a r 1 ,··· ,rn in the differential operator D for a prepared CHY-type expression for n points with h n = σ i .
The relations corresponding to the painted tableaux in Figure.6 are independent from each other, since every relation involves a new coefficient. These relations are obviously complete as well, since they are linear and fix all the coefficients uniquely.

One-loop five-point SYM amplitude
In this and the next sections, we further illustrate the aforementioned method by discussing more examples, the one-loop MHV amplitude for five external particles in N = 4 SYM and the n-gon amplitude for an arbitrary number of particles. At first sight, the one-loop SYM integrand [15] might appear considerably more complicated than the four-point case due to the existence of the Pfaffian. Fortunately for us, supersymmetries, which enter the story in the GSO projection, simplify the integrand/amplitude by a great deal and cast the integrand into a form that is no longer intimidating at all. We first demonstrate this process using the relations among the Jacobi theta functions studied in [50] and the decomposition of the polarization vectors. In the end we arrive at an expression that is a sum of prepared forms and the evaluation of this expression becomes straightforward.

Simplify the one-loop generalized CHY integrand
The generalized CHY form for the one-loop 5-point integrand in SYM is first conjectured as the following by Geyer et al in [15,16] and checked against the corresponding Q-cut expressions up to five points, where P T 5 is the well-known Parke-Taylor factor for five points and its general expression reads P T n = ρ∈Sn 1 σ ρ(1) (σ ρ(1) − σ ρ(2) )(σ ρ(2) − σ ρ(3) ) · · · (σ ρ(n−1) − σ ρ(n) ) , (3.2) in which S n is the cyclic permutation group of n elements. The shorthand notation F denotes the summation of the Pfaffians over the GSO sectors that have even spin structures, namely F = α c α Z α P f (W α ) | τ →∞ where τ is regarded simply as a parameter in the SYM computations, but can be interpreted as the moduli parameter of the corresponding one-loop worldsheet in string theory. 3 Z α denotes the theory-specific partition functions for the GSO sectors and in the case of N = 4 super Yang-Mills reads where η(τ ) is the Dedekind eta function and α labels the GSO sectors with even spin structures, that is, (R+), (NS+), (NS−) respectively. The coefficients c α take care of the GSO projection, and thus take the values {c 2 , c 3 , c 4 } = {+1, −1, +1}. The matrix in the Pfaffian takes the following skew-symmetric form, where the polarization vectors and the momenta of the external gluons are denoted as ǫ's and k's (without loss of generality, we choose the helicities of ǫ 1 and ǫ 2 to be negative and the rest positive throughout this section). S α denotes the fermionic propagator in the RNS formalism of string theory, with α labeling the sector, whose explict form will not be used in this paper. We have also re-parameterized the variables as σ = e 2πi(z− τ 2 ) and the scattering equations in terms of the complex variables z i read θ 1 (z ij |τ ) , with G denoting the bosonic propagator. Now we begin to simplify the Pfaffian (3.4), which by definition is a polynomial in the elements of the matrix (3.4). The terms in this polynomial are all of the following form with coefficients containing the kinetic informations only, is a permutation of (i 1 , · · · , i m ). The diagonal terms C ll of the off-diagonal blocks in (3.4) are the only ones that contribute factors other than S α . The union of all the indices above is simply {1, 2, 3, 4, 5} with each number appearing twice while the sets of indices {i 1 , · · · , i m } and {j 1 , · · · j 5−m } have no overlap. The coefficients c({ǫ}, {k}) and the entries C ll are evidently universal for all GSO sectors and can be pulled out from the the summation over the spin structures; the only factors relevant for the GSO summation are the products of S α 's. Such summations are worked out explicitly at one loop for the N = 4 partition functions in [50], exploiting the properties of the Jacobi theta functions. For m ≤ 3 the weighted sum simply vanishes and thus the surviving terms in the five-point Pfaffian are those with m = 4 or m = 5, which sum up to the following simple results when z i 1 i ′ 1 + · · · z imi ′ m = 0 (this condition is satisfied trivially in our case), Using these relations, the sum of the Pfaffians corresponding to the MHV amplitude over the even spin structures can be directly written down as the following, where ρ{345} and ρ{2345} are the permutation groups S 3 /S 2 4 and S 4 /S 2 × S 2 respectively. This expression does not yet put all the external particles on the same footing. In the N = 4 SYM case, it is possible to make the cyclic symmetry manifest at the level of integrand, by projecting the polarization vectors ǫ µ i onto the external momenta k µ i 's and eliminating the loop momentum ℓ using the scattering equations (3.6). The decomposition of the polarization vector reads, 3 v ρ 4 4 and ǫ ρ 1 ρ 2 ρ 3 ρ 4 denotes the Levi-Civita tensor and v ρ i i denotes the ρ i -th component of the vector v j . Now all the ǫ i · ℓ terms can be rewritten in terms of R i,j 's using the scattering equations, leading to the following form of (3.9) that treats all the external lines democratically, where γ [ij] is defined in [47] as with ρ being a permutation element in S 5 /(S 2 × S 3 ).
Substituting (3.11) into (3.1) and taking the limit τ → 0, the one-loop MHV amplitude for five points reads

Transformations to the prepared CHY integrand
So far we have reduced the Pfaffian in the SYM integrand to a simple expression. Substituting the explicit expressions for the Parke-Taylor factor and applying the linear transformations that turn the original scattering equations to polynomials, we write down the following expressions for the five-point integrand, The factor . The former term above is already of the prepared form (2.12) up to the global residue theorem (which only introduces an overall minus sign to the result), while the latter, although not yet of the form we want, can be easily massaged into one, using the following cross ratio identities derived in [32], where i ∈ [1,5]. The denominator on the right-hand side contains the factors (σ i − σ j ) of non-adjacent pairs (ij) only and these factors are canceled out by the Vandermonde determinant that was introduced by the linear transformations of the scattering equations, leaving the monomial σ i in the denominator, which is just the case of the prepared form up to the global residue theorem.
Collecting all the prepared forms, the integrand (3.14) is rearranged as follows, where and (3.16)

Evaluating the prepared form of the five point integrand
In order to evaluate (3.14), we first construct the differential operators that capture the residues corresponding to the prepared expressions obtained previously. Without loss of generality, we take the following prepared form with h 5 = σ 1 as an example to demonstrate the detailed evaluation (all the other terms can be computed the same way), where H(σ) is a holomorphic function of the σ i 's. The differential operator for this residue is of order 6 and has 210 coefficients a r 1 ,r 2 ,r 3 ,r 4 ,r 5 to be determined.
With all the coefficients fixed for the case of (3.17), it is now straightforward to evaluate the P 1 and B [1,2],1 terms in (3.14), by acting with the differential operator D.
Explicitly, the P 1 term becomes = 3!β 1,2,3,4,5 a 0,0,1,2,3 = β 1,2,3,4,5 2ℓ 5 ℓ 4,5 ℓ 3,4,5 ℓ 2, 3,4,5 . (3.21) The integrand term containing B [1,2],1 is holomorphically equivalent to [12] s 12 ℓ 5 ℓ 4,5 ℓ 3,4,5 (3.22) All other terms are similar. The final results for all the terms are shown in Table 1 and 2, where the correspondence between each term and a particular forward-limit channel in the Q-cut analysis [48] is displayed. Table 1: Final results for the pentagon diagrams. The first line lists the terms while the second line their analytic expressions respectively. The last line shows the forwardlimit channels corresponding to each term, where the cut lines in the q-cut analysis are depicted in Fig. 9. For a given cut line, its momentum is chosen to be ℓ.   Table.1, the first and the second lines show the terms and their respective final expressions, while in the last line their corresponding forward-limit channels are given and the cut lines are depicted in Figure.10. Again, for a given cut line in the q-cut analysis, its momentum is chosen to be ℓ. Figure 10: Box Diagram.
This concludes the computation of the residues at finite poles and at infinity (the computation of the residue at infinity is transformed to evaluating a prepared form with h n = σ n in the same fashion as in [1]).

Analysis for the spurious poles
In this subsection we investigate the spurious poles that arise from the extra solutions to the polynomial scattering equations. Such contribution should vanish in consistent supersymmetric theories and we will show this is indeed true in this particular case. It is easy to observe that the extra solutions all coincide at the position σ 1 = σ 2 = · · · = σ n−1 = 1. Naturally, we prefer to deal with poles at the origin and shift the spurious pole there by the parameter transformation σ i → σ i + 1. Therefore the residue at the spurious pole is associated with the meromorphic form, where the factor H 5 is read off from (3.14) and The product i+1<j (σ i − σ j ) contains all non-adjacent pairs of (σ i − σ j ) with the gauge choice being implicitly σ 5 = 1. The shifted polynomial scattering equations read the following, Note that the leading order terms above do not give rise to a zero-dimensional intersection and we have to apply transformations on these polynomials such that the divisors generated by them coincide at an isolated point. The explicit expressions for such transformations will not be needed and we simply display the new polynomials here,ĥ According to the transformation law, we also need to take care of the determinant of the transformations, which is simply 1 here.
In principle we need to construct a different differential operator following Section 4.4 of [1]; and to this end the hatted polynomials are homogenized in a way such that the higher-order terms are lowered to the same degree as the leading order ones, leading to the order of this operator being 4. Luckily the actual form of the operator is not necessary in this case, since (3.24) is holomorphic at the shifted origin and the lowest term in its numerator is of degree 5. The action of the fourth-order differential operator on (3.24) must vanish when evaluated at the shifted origin.
Since the spurious poles do not contribute, the total residue obtained previously is the final integrand for the five-point SYM amplitude at one loop which simply reads, where the expressions for I P i 5 and I are given in Table 1 and Table 2. Just like its four-point counterpart studied in [1], our result presents a clear one-to-one correspondence with the forward-limit channels in the Q-cut analysis. We expect this nice property continues to hold for higher points.
The method we have illustrated so far generalizes naturally to the one-loop integrands for any number of external particles. For higher-point cases, we expect the integrands can be decomposed into several residues that are associated with the prepared form, or similar expressions with slightly more complicated h n 's.
Similar to the analysis in Section.2, when h n is a monomial, the corresponding local duality theorem leads to the vanishing of some coefficients.
The local duality equations arising from the scattering equations are universal despite the choice of h n , that is to say, the fact that the coefficients in the corresponding differential operator can be related to those with lower or equal ranks only holds for any number of points, and most of the relations among these coefficients remain the same. The construction of the differential operator always boils down to solving for the elementary coefficients. Since the number of the elementary coefficients is quite limited, these coefficients can be solved efficiently.

One-loop n-gon amplitude: A direct evaluation
In this section we discuss the evaluation of the generalized CHY form for one-loop n-gon amplitude. The evaluation can be done directly since all n-gon amplitudes are already of the prepared form. The one-loop n-gon amplitude is conjectured to be [18] where P T 4 is the Parke-Taylor factor. The terms in (4.1), after homogenization, are of form where i||j denotes pairs of neighboring indices in the cyclic order. In terms of the differential operator representation, we have where v (ρ −1 ) means the result of performing ρ −1 permutation on v = {0, 0, 1, 2, 3, · · · , n− 3, n − 2}, and subscripts of v indexes the components of v (ρ −1 ) .These a's have already been obtained in Eq. (2.28). On the other hand, as we will see from later discussion, here the spurious pole does not contribute. Thus the sum over ρ ∈ S n of (4.3) is the final result of the general one-loop n-gon amplitudes.
According to our conjecture, the differential operator D for the residue at the shifted origin is (n − 4)(n − 1) 2 + 2, where n 4. Then it is easy to find that the spurious residue are all vanishing for n-gon with n 4.

Conclusion and outlook
In this paper we use our previously proposed differential operator to compute the residues on the solutions of scattering equations. We find that the coefficients in the differential operator can be determined in a combinatoric way. We present an analytical solution for the differential operator for the prepared forms of the generalized CHY integrands for n-point scattering amplitudes. We then use said differential operator to evaluate the one-loop CHY integrand for five points in N = 4 SYM which is casted into the prepared forms, and the one-loop CHY integrands of the n-gon amplitudes which are naturally of the prepared forms for any number of external lines. In both examples, our final results are identical with the ones obtained through the Q-cut analysis respectively.
Although all the examples studied in this paper can be massaged to the prepared form, more complicated denominators may turn up in the integrands when we investigate higher-point one-loop amplitudes in super Yang-Mills. An immediate followup is to generalize our method for differential operators corresponding to these new ingredients. Once such generalization is worked out, the evaluation of any one-loop integrand in SYM will become straightforward and economical.
At the moment of writing, the construction of generalized CHY forms in super Yang-Mills theory has been promoted to two loops for four external lines. Another natural direction to take is to study the combinatoric structures of the two-loop scattering equations and lift the coefficient-generating method to two loops, which in turn will provide a useful tool in constructing the generalized CHY form for any number of external lines at two-loop level.
Moreover, the method we propose is theory-independent and thus has the potential to evaluate any CHY-type expressions in pure Yang-Mills and Gravity theory. This, in turn, will make the generalized CHY forms an efficient method for computing amplitudes.