QCD in the Color-Flow Representation

For many practical purposes, it is convenient to formulate unbroken non-abelian gauge theories like QCD in a color-flow basis. We present a new derivation of SU(N) interactions in the color-flow basis by extending the gauge group to U(N)xU(1)' in such a way that the two U(1) factors cancel each other. We use the quantum action principles to show the equivalence to the usual basis to all orders in perturbation theory. We extend the known Feynman rules to exotic color representations (e.g. sextets) and interactions (e.g. $\epsilon_{ijk}$). We discuss practical applications as they occur in automatic computation programs.


Introduction
The analysis of particle physics experiments at colliders depends on reliable theoretical predictions for cross sections of scattering processes. In the LHC era, hard processes with more external partons than ever have become accessible, and their analysis is essential for unveiling the physics at the Terascale. This results in a two-sided challenge: on one hand, complex Standard Model processes must be computed with unprecedented precision and on the other hand, the parameter space for many models for physics beyond the Standard Model must be scanned with sufficient accuracy. This situation has prompted the development of computer programs that can automatically compute differential cross sections and sample them efficiently on phase space starting from the specification of a Lagrangian or, equivalently, a set of Feynman rules.
In QCD, quarks and antiquarks come in three colors and gluons in eight. Thus, there are many amplitudes with different color quantum numbers that must be computed and summed in quadrature, eventually. In Feynman diagram based calculations, the contribution of each diagram factorizes into the dependence on color and the dependence on flavor, polarization and momenta. This allows to compute the color factors once and for all. However, for processes with many external particles, the number of Feynman diagrams grows factorially, and more efficient, e. g. recursive, algorithms must be used that can take advantage of the cancellations among diagrams. Unfortunately, the color dependence does not factorize for the whole scattering amplitude and new expansions, like color ordered amplitudes must be employed to allow a separate computation of color factors.
It is therefore worthwhile to investigate efficient representations of the color dependence that work for complete amplitudes, i. e. sums of Feynman diagrams. It turns out that expressing everything in reducible tensor products of the fundamental representation and its conjugate instead of higher irreducible representations is beneficial both for computation and for interfacing to parton showers, fragmentation and hadronization. If there are no interaction vertices with exotic color structures, the resulting amplitudes can always be decomposed into weighted sums of products of Kronecker-δs, so called color flows 1 .
The adjoint representation of U(N) is isomorphic to the tensor product N ×N of the fundamental representation and its conjugate. Thus, the color coordinate a of a U(N)-gauge boson 1 If there are couplings like a totally antisymmetric vertex 3 i,j,k=1 ǫ ijk φ i φ ′ j φ ′′ k in SU (3), the corresponding invariant tensors can also appear in the final result.
in the adjoint representation can be represented equivalently by a pair of (i, j) using the decomposition A j i = A a [T a ] j i . This equivalence is employed in 't Hooft's double line notation [1], in which each U(N)-gauge boson line is replaced by a pair of matter and anti-matter lines for the purpose of computing color factors: It was soon realized that the diagrammatic rules could be extended to accomodate the tracelessness of the generators of SU(N) [2]. Indeed, an efficient diagrammatical algorithm for computing arbitrary color factors in any gauge group was constructed [3] using a generalization of the double line notation to higher representations 2 .
However, the color-flow representation is more than a useful calculational device: Feynman diagrams in the double line notation bear a striking resemblance to dual diagrams [6,7] and indeed the N → ∞-limit is related to a theory of open strings [1]. It is therefore no accident that the successful semiclassical models for the fragmentation of quarks and gluons that are based on string dynamics [8] require a specification of the parton cross sections in the color-flow basis [9]. The earliest examples of Monte-Carlo event generators that combined the partonic cross sections for the hard interactions with QCD-inspired fragmentation models [10] provided a template for the interfacing of hard interactions with fragmentation. It became the standard HEPEVT common block [11], which includes color flow information. Originally developed for interfacing Monte-Carlo event generators for hard scattering matrix elements with libraries or external programs for parton showers, fragmentation and hadronization at LEP1, it was widely adopted and has grown into the Les Houches Accords (LHA) [13].
Therefore the color-flow basis is particularly suited for the computation of partonic scattering amplitudes that will be interfaced with the other components of a Monte-Carlo event generator. Indeed, most generators for hard multi-parton scattering can provide the corresponding information. One can, however, go one step further and not only represent the final result in a color-flow basis, but use Feynman rules in a color-flow basis already for the computation. A lucent derivation of these Feynman rules was given in [14]. Since the derivation in [14] is based on the QCD Lagragian, it is not necessary to use Feynman rules at all, but one can implement modern recursive algorithms [15,16,17] directly in the color-flow basis [18,19,20].
However, the derivation [14] of the color-flow representation is incomplete in two directions: firstly, there is no consideration of interactions with more exotic color structures, in particular beyond QCD with fermionic matter in the fundamental representation, and secondly the discussion is deliberately confined to tree level amplitudes. These two limitations are in fact related, and overcoming them not only of theoretical interest: the most important light Higgs production channel at LHC involves the dimension-5 operator H Tr (F µν F µν ), which arises from a loop and corresponds to an octet-octet-singlet coupling that cannot be described straightforwardly in the framework provided by [14]. Indeed, our generalization of the color-flow representation was prompted by the implementation of the effective Higgs-gluon-gluon vertex in WHIZARD's [18] optimized matrix element generator O'Mega [17]. While it was simple to "fudge" the Feynman rules for a single insertion, it turned out that these prescriptions would give incorrect results for multiple insertions.
This paper is organized as follows: we start with a motivation and description of our implementation of color-flow QCD in section 2 and discuss quantum field-theoretical aspects in more detail in section 3. In section 4, we show that our description is valid to all orders in perturbation theory. In the following two sections we discuss some applications in detail: tree level amplitudes in section 5 and the already mentioned effective interactions from loops in section 6. Before concluding, we describe the implementation of interactions with exotic color structures in section 7. For reference, we repeat the color-flow Feynman rules for QCD with fermionic matter in the fundamental representation in appendix A.

The Model
In this section, we construct a variant of QCD which contains two extra gauge degrees of freedom. The model is set up not for its physics content, but for its usefulness in practical perturbative calculations. In fact, as we will show in later sections, in physical quantities the extra degrees of freedom cancel, so all predictions are identical to ordinary QCD.

QCD Preliminaries
We consider QCD with a single fermionic matter species in the N-dimensional fundamental representation of SU(N). While we are really interested in the N = 3 case, it is sometimes useful to explicitly keep the dependence on N, as we will do through most of this paper.
Furthermore, we will not assume that the matter representation is vector-like as in QCD, although we do follow the QCD notation with Dirac fermions in the Lagrangian.
Let us first recall basic facts of perturbative QCD. The perturbation series is derived from a Lagrangian which splits into two parts, a gauge-invariant part L inv and a gauge-fixing part L gf . The gauge-invariant part is given by 3 The gauge-fixing part depends on the chosen gauge-fixing procedure. In a manifestly covariant formulation with linear gauge-fixing, it takes the form Since the matter fields are in the fundamental (defining) N-dimensional representation of SU(N), we choose to represent the Lie-algebra valued fields G, A, B, c,c by traceless N × N matrices. 4 The field-strength tensorĜ can be expressed in terms of the gauge potentialÂ aŝ Furthermore, the gauge-fixing term involves a Nakanishi-Laudrup (NL) auxiliary fieldB, a Fadeev-Popov ghost fieldĉ, and an antighost fieldĉ.
We may assume a manifestly gauge-invariant renormalization procedure such as MS, so that the Lagrangian retains the form (2,3) in each order of the perturbative expansion. 5 The fields in the Lagrangian are understood to be renormalized, order by order, such that the model contains only two (renormalized) real parameters, the gauge-coupling g and the gauge-fixing parameter ξ. The latter drops out of physical quantities.
The gauge group SU(N) is a subgroup of the general linear group GL(N), therefore the fields (e. g.,Â), in the N × N matrix representation, obey algebraic constraints: The constraints are automatically satisfied if we introduce the usual basis for the Lie algebra representation, T a (a = 1, . . . N 2 − 1), and writeÂ = a A a T a . The basis elements obey the hermiticity and trace conditions:

The Color-Flow Representation
In the following, we will eliminate the basis T a , and instead work with the individual matrix elements of the gauge fields in the chosen Lie algebra representation [5]. This approach is known as a color-flow representation. For concreteness, we consider the gauge potentialÂ. Without algebraic constraints, a N × N matrix field contains N 2 complex degrees of freedom.
We have to implement the algebraic contraints (5). The hermiticity condition iŝ In effect, each unordered index combination (ij) carries only one, instead of two, complex degree of freedom. We can understand this as if each ordered index combination (i, j) carried Similarly, matter propagators are represented by a single directed color line. External states are represented by a terminating line (matter) or two opposing terminating lines (gauge). As long as we consider only the interactions following from the Lagrangian as quoted above, the rule of matrix multiplication trivially ensures that color is conserved at each vertex, i. e., lines end only at external states.
While this graphical approach is straightforward, we also have to enforce the trace condition from (5), explicitly which reduces the number of degrees of freedom by one. The implementation of this constraint amounts to additional modified color-flow patterns in Feynman graphs [14]. While the construction of these patterns is easy to understand, they are, by definition, a non-local modification of the naive color-flow expansion. This is a complication in automatic, in particular non-diagrammatic, algorithms, which we want to avoid.
Instead, we follow a field-theoretical approach and set up a modified QCD theory. This theory has the algebraic constraints inherent in its field content, such that they are automatically satisfied by the generated color-flow amplitudes. In essence, it generates (i) terms that complete the A field to a U(N) matrix so that it has no trace condition and naive color-flow Feynman rules apply, and separately (ii) terms that subtract the spurious contributions without complicating the color-flow rules. To make the diagrammatic expansion unambiguous, the two fields have to be formally independent of each other.

Singlet and Phantom Gluons
We denote the field which complements the gluon field matrix as the singlet gluon A 0 . This field is associated with a "zeroth" generator T 0 which is not traceless but satisfies and let a run from 0 to N 2 − 1. Explicitly, (We choose upper indices for the fundamental and lower indices for the anti-fundamental color representation.) The matrix field resulting from combiningÂ with A 0 has N 2 independent components, the color-flow gluons A i j . Analogously, for the NL and ghost fields, we replace the traceless fields by unconstrained fields B, c,c and singlet fields B S , c S ,c S , respectively. The field strength G is modified accordingly.
For the subtraction terms, we introduce a phantom gluon fieldÃ. The phantom gluon is an independent U(1) gauge boson with the wrong sign in the propagator. Like the other gluons, it couples to the matter fields via the QCD coupling g. It does not couple to gluon fields. We also introduce the corresponding NL fieldB. (We may also introduce phantom ghost fields, but these can be omitted with linear gauge fixing, since the U(1) gauge group is abelian.)

Color-Flow QCD
We can now set up the Lagrangian for the color-flow version of QCD, that includes the singlet gluon (as a component of the gluon matrix A i j ) and the phantom gluonÃ, together with their associated field strength, NL and ghost fields. where and the decoupling ghosts for the phantom gluon have been omitted.
A priori, this theory is different from QCD. We will denote it as "color-flow QCD" and we will show its equivalence to QCD in section 3. It contains N 2 instead of N 2 − 1 independent degrees of freedom in the A field, plus the additional U(1) ′ gauge fieldÃ with wrong-sign propagator. Thus, it can be regarded as a U(N) × U(1) ′ or, equivalently, SU(N) × U(1) × U(1) ′ gauge theory. In a path-integral formulation, there are independent integrations over all N 2 + 1 gauge components. The singlet field A 0 is present, but hidden in the gauge-field matrix.
In the Lagrangian (12), we have encoded the algebraic properties of the color-flow decomposition as a (perturbative) quantum field theory on its own. Therefore, we can expect that amplitudes calculated from this field theory manifestly exhibit the advantages of this approach.
In calculations, we therefore retain the decomposition of the matrix-valued fields in terms of their matrix elements, A i j . As explained above, these non-hermitean vector bosons are represented by double lines, one for the color and one for the anticolor flow, flowing in opposite directions. The vertex color factors are Kronecker deltas, which diagrammatically correspond to color lines continuing through the vertex. The phantom gluon, on the other hand, is analogous to a photon. In particular, it carries no color flow. Color lines always terminate at external states. Hence, for any tree diagram, the color factor is simply N p , where p is the number of distinct color lines in the diagram (after squaring). Since this is determined by the colors of the external states, the whole color algebra becomes trivial. The price for this is a proliferation of extra Feynman diagrams (which have identical kinematics, however, and thus do not blow up the calculation.) We list the Feynman rules that follow from (12) in Appendix A.

Example
Let us look at the exchange of a gluon between fermion lines. In ordinary QCD, this is represented by the diagram In color-flow QCD, we get two diagrams, one for gluon and one for phantom exchange Usually, this color-flow decomposition is obtained by applying the relation (14) backwards. The Feynman rules derived from the color-flow Lagrangian generate it directly.
Also at the squared-amplitude level, the color-flow representation leads to a reorganization of the color-flow calculation, which considerably simplifies the algorithm. There is a further benefit: common parton-shower and hadronization algorithms expect the color connections to be expressed in a color-flow basis. We obtain both the exact amplitude (at the given order) and, by omitting the interference terms and all terms containing external phantoms, the projection on the possible color-flow patterns in a single step.
In summary, we work with a trivial version of color algebra since the complete color flow information is represented on a graphical basis. Subtractions in the color factors are accounted for by extra diagrams which involve the phantom gluon. No approximation is involved.

Field-theoretical Considerations
In the following, we demonstrate that color-flow QCD as described above is equivalent to the original SU(N) theory, i. e., A 0 andÃ cancel each other in physical quantities. This equivalence may appear trivial, at least in tree-level amplitudes. However, it bears some subtleties, in particular if it is to be extended to all orders. To promote the approach from a technical device to a consistent field theory, we have to prove that the equivalence can be maintained, by suitable renormalization conditions, to all orders in perturbation theory.
This situation is similar to the covariant quantization of QED, where the scalar and longitudinal photon components, which are both unphysical, cancel in physical quantities. In a nonabelian theory, this requires the additional introduction of ghost fields [21,22]. In the Fock space, the unphysical fields lead to indefinite-norm states, but it is possible to build a Hilbert space out of equivalence classes of physical states in which the S matrix is unitary. As we will show below, the same argument applies in the present case: the phantom gluon generates negative-norm states, but it is always produced coherently with the singlet gluon (which is now hidden in the N × N gluon matrix), and therefore a unitary S matrix is well defined.

Algebraic Identities at Lowest Order
A consistent formulation of a quantum field theory is based on the symmetries that the effective action, including the integrated Lagrangian as its tree approximation, has to satisfy order by order in a perturbative expansion. These are representable as algebraic identities [22].
In the present subsection, we look only at the lowest order, so we temporarily ignore the complications that necessitate the replacement of gauge symmetry by BRST symmetry, ghost fields, etc. Given the fields A i j ,Ã, B i j ,B,ψ, ψ, we observe that the lowest-order action Γ 0 = L from (12), with the ghost term omitted, satisfies the Ward identity of U(N), where the gauge transformations of the fields are This Ward identity decomposes into separate identities for the SU(N) and U(1) factors. We are interested only in the U(1) part. We can project this out by evaluating the Ward identity (15) for ω = ω 0 1, This is accompanied by the gauge-fixing condition for the singlet Furthermore we define, for any matrix-valued field M ∈ {A, G, . . .}, the linear combinations and the dual functional operators Using (19), we can rewrite terms appearing in the Lagrangian as Computing the derivative of the lowest-order action with respect to A andÃ and combining the results, we obtain the phantom equation It expresses the cancellation between singlet and phantom terms. Since A + is an abelian gauge field, applying δ/δA ν + to the ghost term (13) yields zero, so the phantom equation holds for the complete lowest-order Lagrangian in the BRST formalism.

Fock-Space Cancellation
To clarify the physical implications of the phantom equation (22), let us use (21) to express the singlet and phantom fields ( 1 N Tr A,Ã) by (A + , A − ). The remainder of the matrix A is traceless, this is the original SU(N) gluon with all of its interactions. We must show that the extra fields A + and A − do not contribute to observable quantities.
We eliminate the B fields via their equations of motion. The terms in the resulting Lagrangian that depend on A + or A − take the form The propagator interchanges A + and A − , but only A − couples to matter. Couplings to ordinary gluons are absent in (23), because G µν As a result, whenever a singlet/phantom gluon is created by a matter current, it cannot be annihilated, and vice versa. In short, A + and A − do not introduce any observable interactions.
In fact, we may freely add any term to the effective action that depends on A − at least linearly, but has no dependence on A + . Any such interaction is unobservable. This pattern will continue to hold to all orders, unless a loop diagram induces an interaction of some current with A + . However, the phantom equation (22) prohibits this: it tells that the interactions of A + are completely accounted for by the leading-order Lagrangian, up to a renormalization of the coupling constant. Therefore, if we can maintain the phantom equation on the effective action to all orders, observables computed from the color-flow QCD theory are identical to those from ordinary QCD.
In actual applications, we do not express the Lagrangian in terms of A + and A − . Instead, we work with the U(N) field A and the U(1) ′ phantomÃ, as we did in the previous sections.
Hence, the Fock-space cancellation is not entirely trivial since it involves cancellations among graphs that belong to different gauge groups, at least superficially. We have to enforce the phantom equation explicitly in order to guarantee that no A + interaction arises in the effective action.

Symmetries
If the phantom equation (22) is combined with the Ward identity (17) and the gauge-fixing condition (18), we obtain an analogous Ward identity and gauge-fixing condition for the phantom gluon. The complete symmetry, at lowest order, is U(N) × U(1) ′ , with the additional condition that all gauge couplings are equal (to the SU(N) gauge coupling).
In the previous subsection, we introduced linear combinations auf the extra gauge fields, A + and A − . Correspondingly, we may identify gauge groups U(1) + and U(1) − which are orthogonal combinations of the original extra U(1) and U(1) ′ gauge groups. These gauge invariances are described by the Ward identities which follow from the Ward identity (17) and the phantom equation (22).
The corresponding gauge-fixing conditions are Note that the matter fields are neutral under U(1) + . This is in line with the observation that However, by the argument in the previous subsection, we see that a U(1) − anomaly will not contribute to observables, so it does not invalidate unitarity of the S matrix projected onto physical states.

Renormalization
In this section we outline a proof that the color-flow version of QCD can be extended consistently to all orders in perturbation theory. We follow the algebraic renormalization procedure, as explained in detail in Ref. [22].

Conditions Imposed on the Effective Action
with the linear gauge-fixing equation These identities involve the gauge, gauge-fixing, and ghost fields A i j , B i j , c i j ,c i j , expressed in the color-flow basis, the matter fields ψ,ψ, and the sources for the BRST variations, ρ A , ρ ψ ,ρ ψ , ρ c . Furthermore, the effective action depends on the phantom fieldÃ and its associated NL fieldB. The dependence on those fields is accounted for by the phantom equation (22), which we state here in the form and by the gauge-fixing condition forÃ, The algebraic conditions are supplemented by renormalization conditions which fix the field normalizations and a minimal set of interaction parameters. In the present model, these are just the gauge couplings.
Alternatively, in a minimal-subtraction scheme such as MS, we read off renormalization conditions after regularizing the loop integrals and subtracting the divergencies. This can also be done here, but with a caveat: we may use any renormalization scheme for the SU(N) coupling and for the ordinary SU(N) fields, but in order to keep the simplicity of the approach, we impose the condition that the singlet U (1)

Lowest-Order Solution
The next step in the renormalization procedure is to identify the lowest-order effective action.
We can verify that the Lagrangian (12) is the unique solution of the constraints of the preceding section (with correct gauge group and representation) which is an integrated local polynomial in the fields of dimension equal to four.
The BRST invariance of the theory, and of the lowest-order approximation in particular, is encoded in the ST identity (27) and leads to lowest-order BRST transformations of the form The U(1) part decouples. The corresponding ghost fields Tr c and Trc are free and can be dropped, such that this part of the ST identity can be replaced by the U(1) Ward identity (17).

Inductive Renormalization
For an inductive proof of renormalizability, we can now assume that we have found a Lagrangian which generates a renormalized effective action at order n that satisfies all conditions. The The only non-standard part is the phantom equation (29). We do not know whether the right-hand side is zero at order n + 1. The quantum action principles [23,24,25,26] tell us, however, that the equation must take the form where the effective action Γ (n) is evaluated from the order-n Lagrangian, including terms in the result up to order n + 1, i. e., diagrams with n + 1 loops. The operator X (n+1) ν on the right-hand side is a local polynomial in fields and derivatives of dimension three with the same quantum numbers as the left-hand side.
Taking the divergence of (32), we obtain This is precisely the Ward identity for the U(1) + symmetry (25) with an anomaly on the right-hand side. We have to verify that the right-hand side is in fact zero.
It is well known that the only possible obstruction to such a U(1) Ward identity is the Adler-Bell-Jackiw (ABJ) anomaly [12], i. e., an operator corresponding to a triangle diagram with three external gauge bosons. In a vector-like theory such as QCD, this vanishes. Nevertheless, we may consider a chiral SU(N) theory where this is not necessarily the case.
With the hypothesis that the phantom equation (29) is valid for the renormalized fields and interactions at order n, we can apply the change of basis (19) to the order-n renormalized fields and rewrite this equation as in complete analogy with the tree-level equation (22). Taking the second derivative with respect to matter fields of this relation, we immediately see that A (n) + does not interact at all. In particular, the interaction (ψ / A + ψ) (n) vanishes. Hence, there is no triangle diagram involving A + that can contribute to the ABJ anomaly. Integrating this zero, we conclude that the right-hand side of (33) vanishes at order n + 1. The U(1) + symmetry is preserved. 7 Since the right-hand side of (33) vanishes, we have ∂ ν X (n+1) ν = 0. A divergenceless vector must be the derivative of an antisymmetric tensor, thus Y (n+1) is a SU(N)-invariant tensor of dimension two, a local polynomial of the fields. The only possibility is with constants a (n+1) and b (n+1) . If we now modify the Lagrangian by the finite counterterms and re-evaluate Eq. (33), this term is cancelled. In summary, by adding local counterterms to the Lagrangian, we are able to satisfy all required conditions to order n + 1. In particular, we establish the vanishing of all nontrivial interactions of A + , i.e., the exact cancellation of the singlet and phantom interactions.
We note that the cancellation does not apply to the orthogonal combination A − , so in the presence of chiral matter, an anomaly operator involving A − is allowed. However, as argued In practice, adding those counterterms is necessary and natural. In the usual MS renormalization scheme, the singlet gluon propagator will receive a renormalization different from the octet gluon propagator, because only the octet has self-couplings. In order to keep the colorflow scheme simple, we require finite counterterms that restore the equality of the propagator residues, and thus of the SU(N) and U(1) couplings, at the next order. Then, we also have to renormalize the phantom propagator, i.e., the U(1) ′ coupling, to the very same value. This renormalization coincides with (37). We can do this freely since the values of the U(1) and U(1) ′ couplings do not enter any observable quantity.

Applications: Tree-Level Amplitudes
The color-flow approach is particularly useful for the automatic calculation of tree-level (squared) amplitudes. In this section, we choose a few simple examples that show how this works in practice. Instead of starting with ordinary QCD diagrams, an implementation may construct the amplitude diagrammatically from scratch, treating U(N) gluons and quarks of definite color and phantoms as distinct, ordinary particle species. The cross section is computed by squaring diagrams, including all interference terms. We recall that U(N) gluons do not interfere with phantom gluons.

Algorithm
Each squared or interference diagram has a color weight W which is simply a signed integer power of N, given by In automatic calculations, the color-flow approach has the advantage that the combinatorics of constructing Feynman graphs is already implemented, so the proliferation of diagrams does not raise a bookkeeping problem. (The algorithm should avoid to compute identical kinematics twice, however.) Counting distinct colors is rather simple and can be done, for tree graphs, by looking at the external state. In fact, common conventions [13] require an event-generating program to classify the external state in terms of color connections. On the other hand, computing color factors algebraically in a T a -f abc basis requires some additional infrastructure, and the transformation to the color-connection basis has to be done explicitly.
Algorithms that do not expand an amplitude in diagrams but compute off-shell wave function recursively [15,16,17] benefit even more from the color-flow approach. Assigning color factors to diagrams is not applicable there, so a straightforward implementation of QCD color algebra would require keeping a color degree of freedom in each off-shell wave function while constructing the amplitude, and applying a color-matrix multiplication at each vertex. By contrast, in the color-flow approach each color line can be understood as labeling an independent particle species, hence the algorithm need not know about color at all. It just has to distinguish particle species and their respective Feynman rules.
In the following examples, we verify the color-flow result against the equivalent color-algebra result, where we evaluate the appropriate trace of T a and f abc matrices directly. We are not interested in the kinematical part of the calculation that involves propagators, Lorentz factors, and integrations over momenta, so we omit them. We also leave out the QCD coupling which is attached to each vertex, and the color-averaging factor for the initial state. We just quote color weights.

Quark-Antiquark Scattering
We are looking at the process in pure QCD, so at tree level there is a single diagram, gluon exchange.
To eliminate all open color indices, we square the amplitude and sum over colors. In the diagrams, squaring connects the final state of the amplitude and its complex conjugate; this is indicated by a vertical dashed line. Similarly, the initial states of the amplitude and its complex conjugate, the open lines at the left and right margins are also understood to be pairwise connected.
⇐⇒ + + + For this squared amplitude, standard color algebra yields as its color weight.
Color-flow QCD gives four diagrams for the squared amplitude, which differ only in the color factor. We recall that there is a factor N for each closed color line, and a factor −1/N for each phantom propagator. Including 1/ √ 2 for each vertex, we obtain We observe that there is some redundancy in this case which could be eliminated before computing the result. In any case, we can read off the color factor directly without using color algebra.

Four-Jet Production in e + e −
Let us now consider where we can ignore the colorless initial state when drawing color-flow diagrams. For the amplitude, we have three QCD diagram structures which decompose into ten color-flow diagrams: We now look at the squared color-flow diagrams. Symmetries between diagrams simplify the calculation. Summing over all diagrams with common kinematics, we recover the QCD color weights as linear combinations of powers of N. We show three terms, the remaining ones are equivalent regarding their color flow: In the color-algebra formalism, we compute the color weights which agree with the color-flow results, as required. We are particularly interested in loop amplitudes that can be inserted as effective vertices in tree-level diagrams. We absorb the color sum over internal closed color lines in the corresponding vertex factor. As a result, we again can deduce the remaining overall color factor of any amplitude merely looking at the external lines.
An important example is the coupling of Higgs and electroweak bosons to gluons, which occurs first at one-loop level. As in the previous section, we are interested only in color flow, so we ignore all kinematical and coupling factors.
Here, there are no interferences. To get the correct factors, we only have to remember that each phantom gluon comes with a factor −1/N, while the singlet gluon carries color, which yields a factor N when summed over. As in the tree-level case, we just have to count color lines crossing the cut to obtain the color weights of the squared diagrams.

gg → HH
Again, the correct factor is recovered, albeit many cancellations are involved.
At this point, we digress somewhat and discuss a possible way to remove this redundancy before the sum is computed. Depending on the method by which the squared amplitude is computed, this may be useful for improving efficiency.
We would like to apply the direct cancellation between a phantom and a singlet line, where it appears obvious. If the amplitude is completely expanded in terms of diagrams, explicit singlets appear always attached to one of the effective vertices that involve a closed color loop, including the loop pertaining to the effective non-QCD interaction that we are considering here. In any case, if we apply this argument to the gg → HH amplitude before squaring, we obtain the much simpler result When we compute the square, we have to take into account that singlets, being contained in the U(N) gluon matrix, interfere with U(N) gluons. This interference acts as a projection operator which turns a U(N) gluon into a singlet, on the other side of the cut: Graphically, we can again apply the cancellation procedure and obtain the final diagrammatic There is no redundancy left. Note that, effectively, the cancellation procedure has switched the sign of the phantom-loop graph, so we would have obtained this result if we just had included this graph with switched sign, but no graphs that explicitly involve singlets. This observation might be generalized and incorporated into the algorithm.

H → ggg
The effective one-loop vertex H → ggg has two color structures, f abc and d abc , which originate from the difference and sum of the possible loop orientations, respectively. In the color-flow basis, the Feynman rules are The square of the f term is straightforward: The interference vanishes as expected. The square of the d term yields Eliminating the external singlet lines as before would allow us to combine the last two terms, i. e., which is slightly simpler. under discussion provide extra colored fields and interactions that violate this rule. Automatic computation programs should be able to deal with these extensions, which is possible in the color-flow approach. We consider the two simplest cases, epsilon-tensor couplings and colorsextet fields.

Epsilon Tensor Couplings
The SU(N) algebra has a second invariant besides the Kronecker delta δ i j , the epsilon tensor ǫ i 1 ...i N with N upper (or lower) indices. Standard QCD does not contain any interaction proportional to the epsilon tensor -such interactions are not compatible with the Standard-Model spectrum and symmetries up to operator dimension four. However, they may be introduced together with new exotic particles, e. g., as remnants of GUT spectra. For instance, an isosinglet scalar color-triplet D with down-quark charge could become a diquark-type particle via a Yukawa-type coupling ǫ ijk D i d j u k .
In a tree-level amplitude with one epsilon-tensor vertex, there are N open color lines that end (or begin) within the diagram without a matching anticolor. The coupling induces an ordering of those lines; exchanging any two of them switches the overall sign (for discussion of the large N limit of such terms cf. [27] The latter color structure occurs, e. g., when considering single diquark production and decay at a hadron collider.
If we want to cast these observations into an algorithm, we have to deal with several issues.
We have to keep track of the ordering of color lines, because this determines the overall sign.
More importantly, in the context of algorithms that do not expand everything in terms of Feynman diagrams, combining interactions as shown above is a nonlocal procedure that cannot easily be incorporated.
As a possibility for circumventing these difficulties, we propose to treat epsilon-tensor cou- By matching the topology with one of the situations shown above, we obtain the correct color factor.

Color-Sextet Particles
Color-sextet particles or, more generally, particles in the symmetric bi-fundamental representation of SU(N), are not present in the Standard Model. However, they might exist at mass scales accessible at the LHC. A possible source are extended Higgs (super-)multiplets in GUTs.
For definiteness, we consider a color-sextet scalar particle σ. An explicit representation is a complex symmetric matrix in color space with six independent entries σ 1,2,3 andσ 1,2,3 , with components σ ij , where σ ij = σ ji . We also need the antiparticle σ ij . This particle couples to gluons and gluon pairs. It may also couple to colorless states (e. g., Higgs), and linearly to quark pairs. The Lagrangian is The covariant derivative, in color-flow QCD, is given by We derive the Feynman rules. The propagator has two components which symmetrize the color flow. The same symmetrization is needed for external sextet particles: Gluons interact symmetrically with both color lines: pµ + qµ pµ + qµ −2(pµ + qµ) Couplings to colorless particles are simply iλ (78) and the coupling to quark pairs has the Feynman rules We list some color weights for simple processes, without dividing by in-state multiplicity:

Conclusions
The expansion of QCD amplitudes in a color-flow basis has been known as a useful device in various contexts of perturbative and non-perturbative calculations. In this paper, we have demonstrated that it can also be understood as a field theory of its own, a priori different from, but equivalent to, standard QCD. We have shown that this field theory is well defined, renormalizable, unitary, and yields observable predictions identical to ordinary QCD.
For practical applications, color-flow QCD has advantages in the context of automatic calculations, where it systematically generates color-connected amplitudes that can be matched to parton-shower and hadronization algorithms, avoiding some extra effort that pertains to color algebra and change of bases. In particular, it is useful for algorithms which do not work with a diagrammatic expansion of the amplitude.
We have implemented color-flow QCD in the O'Mega matrix element generator, which is the tool for tree-level amplitude generation contained in the WHIZARD event generator package [18] together with its parton shower generator [28]. The color-flow approach also provides a convenient way of incorporating exotic color interactions, such as epsilon-tensor couplings and color-sextet particles. If desired, the methods presented in the present paper can readily be extended to other exotics, e. g., color-decuplet fields.
For higher-order calculations, color-flow QCD can be taken at face value, provided the renormalization procedure is properly implemented (cf. Sec. 4), which is essentially trivial.
Concerning mixed QCD-electroweak processes at loop level, the examples in Sec. 6 show that there are minor technical complications, which however do not invalidate the procedure. A comprehensive treatment of NLO calculations would also require a color-flow exposition of (dipole) subtraction and parton splitting kernels. This is not discussed in the present paper, but straightforward and part of a different publication [29].
Finally, we remark that the methods developed in the present paper could be generalized and applied to other gauge groups, including exceptional groups, using the completeness relations from [3]. Propagators: