{'fraction_non_alphanumeric': 0.03934680087280918, 'fraction_numerical': 0.013415921728725277, 'mean_word_length': 4.314281439365602, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Webs and Posets

The non-Abelian exponentiation theorem has recently been generalised to correlators of multiple Wilson line operators. The perturbative expansions of these correlators exponentiate in terms of sets of diagrams called webs, which together give rise to colour factors corresponding to connected graphs. The colour and kinematic degrees of freedom of individual diagrams in a web are entangled by mixing matrices of purely combinatorial origin. In this paper we relate the combinatorial study of these matrices to properties of partially ordered sets (posets), and hence obtain explicit solutions for certain families of web-mixing matrix, at arbitrary order in perturbation theory. We also provide a general expression for the rank of a general class of mixing matrices, which governs the number of independent colour factors arising from such webs. Finally, we use the poset language to examine a previously conjectured sum rule for the columns of web-mixing matrices which governs the cancellation of the leading subdivergences between diagrams in the web. Our results, when combined with parallel developments in the evaluation of kinematic integrals, offer new insights into the all-order structure of infrared singularities in non-Abelian gauge theories.


Introduction
Wilson lines continue to generate significant amounts of interest, due to their role in a variety of phenomenological and theoretical applications . For example, they govern the structure of infrared singularities in scattering amplitudes, which lead to large kinematic logarithms which have to be summed up to all orders in perturbation theory [15][16][17]. They also govern scattering in the high energy limit in both gauge theories and gravity [55][56][57][58][59][60][61][62][63][64][65][66][67][68][69][70][71]. For these and other applications, it is important to classify the behaviour of correlators of products of Wilson line operators. Whilst much is known about the case of two Wilson lines meeting at a cusp [72][73][74], the general case of more than two Wilson lines has been studied only recently. The anomalous dimension has been calculated up to two-loop order for both massless [75,76] and massive [77][78][79] particles, where in the former case it was found to have the same colour structure as the one-loop result. This was later explained by constraints on the structure of infrared singularities in massless scattering amplitudes which follow from factorisation and rescalaing symmetry, leading to the formulation of the dipole formula [80][81][82]), a minimal all-order ansatz for the anomalous dimension. Possible corrections to this form have been investigated further in refs. [65,66,[83][84][85][86]. Further progress in this area requires explicit calculations of Wilson line correlators at higher loop orders. In practice this means calculating the exponents of correlators of Wilson lines, from which the relevant anomalous dimensions can be extracted. This itself necessitates the development of efficient computational techniques.
In a recent series of papers [87][88][89][90] (see also [91]), a diagrammatic method has been developed, which allows us to calculate the exponents of Wilson line correlators directly. This provides an efficient framework in which to address higher loop contributions to the multileg soft anomalous dimension, and is a generalisation of the well-known diagrammatic exponentiation of correlators of two Wilson lines (or a Wilson loop) in terms of webs [72][73][74]. As is explained in detail in ref. [87], however, this generalisation is highly non-trivial. Whereas webs in two-parton scattering are single irreducible diagrams, which are each individually free of subdivergences associated with renormalisation of the cusp, in multi-parton scattering webs are closed sets of diagrams related by permutations of gluon emissions on the Wilson lines (we review this material in more detail in section 2). The general form of a given correlator of Wilson line operators Φ i is then given by where the contribution of each web is Here the sums are over all diagrams D, D in the web, whose kinematic and colour parts are given by F(D) and C(D ) respectively. The matrix R DD consists of rational numbers of combinatorial origin, and acts to entangle the colour and kinematic degrees of freedom. Each web has its own web-mixing matrix, and these are known to have special properties. Chief among these is idempotence (proven in ref. [88]), which implies that web mixing matrices are projection operators, with eigenvalues in the set {0, 1}. Each unit eigenvalue is associated with a combination of kinematic factors which survives in the exponent, and which has an associated exponentiated colour factor (ECF), consisting of a superposition of colour factors of individual diagrams in the web. It has recently been shown that all such ECFs correspond to colour factors of diagrams in which all gluons are connected (once the Wilson lines have been removed), thus generalising the so-called non-Abelian exponentiation theorem from the two-line to the multiline case [90].
Other properties of web-mixing matrices are known, such as the fact that their rows sum to zero [87,88], and that their columns appear to sum to zero after contraction with appropriate weight factors [89]. These features are related to two underlying principles: (a) colour factors in the exponent must be fully connected (as described above); (b) the corresponding combinations of kinematic factors are only allowed to contain specific higher-order divergences, which can be cancelled by highly constrained subtraction terms involving subdiagrams of the given web (see [89]). It is therefore clear that a deeper understanding of the combinatorics of web-mixing matrices can provide a valulable insight into the structure of infrared singularities in non-Abelian gauge theories, to all orders in perturbation theory. This is only part of the story, however. A full understanding of higher order singularity structures requires calculation of the necessary kinematic integrals. To this end, new computational techniques are becoming available, which offer the possibility of calculating certain families of kinematic integral to any order in perturbation theory [92,93]. Clearly, analysis of the structure of web mixing matrices and that of kinematic integrals should proceed in parallel.
The aim of this paper is to investigate the combinatorial structure of web mixing matrices in more detail. We will present a useful classification of webs based on the number of external lines they connect before focusing on the special case of webs which connect L + 1 lines with L single gluon exchanges. We will present a simple formula for the rank of any web-mixing matrix in this class and, thus, for the number of independent connected colour factors such a web contributes to. We will then examine the results of ref. [94], a recent study which considered web-mixing matrices from a pure mathematical point of view, and which related their combinatorics to that of order-preserving maps on partially ordered sets (posets), for which a large mathematical literature exists. Using the poset language, we will be able to completely solve for the web-mixing matrices for two special families of webs, at arbitrary order in perturbation theory. The calculation of the corresponding kinematic integrals will be considered elsewhere [92,93]. Finally, we will interpret a recently conjectured column sum rule for web-mixing matrices [89] in the poset language, and look at this in detail for the two particular web families whose web-mixing matrices we are able to exactly solve for.
Although much about web-mixing matrices and their associated kinematic combinations remains unexplored (particularly for webs which contain three and four-gluon vertices away from the Wilson lines), their classification through combinatorial methods offers new insights into the structure of non-Abelian gauge theory. The web-mixing matrix solutions we present here correspond to allorder properties in the exponents of scattering amplitudes, thus correspond to very deep structures in perturbation theory. There is a clear motivation for pushing such investigations to broader web families, in tandem with the calculation of the relevant kinematic integrals.
The structure of the paper is as follows. In section 2 we review the definition of web-mixing matrices, in a self-contained fashion. In particular, we give a combinatorial definition of web mixing matrix elements (first presented in ref. [88]) that will be used frequently throughout the rest of the paper. In section 3 we consider how to systematically classify webs, and provide an expression for the rank of a particular class of webs, as mentioned above. In section 4, we introduce terms and properties relating to posets, which will not necessarily be familiar to a particle physics audience. In section 5 we consider two particular classes of webs, and use the poset language to obtain explicit solutions for their associated web-mixing matrices, before verifying that they encode fully connected colour factors. In section 6 we discuss the interpretation of the weighted column sum rule of ref. [88] in terms of linear extensions, and verify that this holds for the two web families considered in section 5. In section 7 we discuss our results and conclude.

Webs and web-mixing matrices
The aim of this section is to briefly review previous results on webs and their mixing matrices, which will be useful for the rest of the paper. Note that different conventions exist in the literature for what is meant by a multiparton web [87,91] (see in particular appendix A of ref. [88] for a discussion of this point). Throughout this paper, we adopt the terminology of refs. [87][88][89][90].
A multiparton web is a set of diagrams contributing to the exponent of a vacuum expectation value of Wilson line operators meeting at a point (alternatively, to the soft gluon part of a QCD scattering amplitude). An example is shown in figure 1, which illustrates a web connecting three Wilson lines (out of a total of four) at two-loop order. This consists of a pair of diagrams, which are related by permutations of gluon attachments on the external lines. Furthermore, the set is closed under such permutations: permuting the gluons in diagram (a) leads to diagram (b), and vice versa. In general, a web is any such set of diagrams closed under permutations of gluon attachments. Note that the two-loop web of figure 1 consists of distinct gluon exchanges. However, there could also be three and four-gluon vertices off the eikonal lines. This does not affect our definition of a web, which only concerns permutations of gluon attachments associated with the hard external (Wilson) lines. Many examples of webs can be found in refs. [87][88][89][90], and we will see more throughout the present paper. Note that any diagram that consists of a single connected piece (when the Wilson lines are removed) is a web by itself; we will not be interested in such webs in this paper as they do not present any combinatorial problem.
We may label a given web using a notation (n 1 , n 2 . . . n L ), where n i ≥ 0 is the number of gluon attachments on parton line i, and it is assumed that there are L lines in total. Note that this notation does not uniquely specify a given web. In particular, it does not specify possible attachments of gluons off the eikonal lines, via three and four-gluon vertices. This will nevertheless be a convenient notation in what follows, since our interest here is in multiple gluon exchange diagrams with no three or four gluon vertices. Given a web W , we may consider all diagrams D ∈ W . Each diagram D has a kinematic part F(D) and a colour factor C(D). Then the contribution of W to the exponent of the Wilson line correlator is given by eq. (2), where R DD is a web-mixing matrix. The existence of web-mixing matrices was first derived using statistical physics methods (the "replica trick"), which also provides an algorithm for how to calculate them [87]. One may also provide a closed-form combinatorial formula for the element R DD . To this end, one needs the notion of a partition of a web diagram D ∈ W . An m-partition of diagram D is a colouring of D with a number of distinct colours m ≤ n c , where n c is the number of connected subdiagrams 3 in D. This is illustrated in figure 2, which shows an example web diagram at three-loop order together with its partitions (note that the diagram itself is an m-partition with m = 1). Although thinking of partitions as colourings gives an easy way to represent them, we actually use another notation to represent partitions: an m-partition of a diagram D will be represented as a set with m elements, each of which consists of the usual diagram notation, but only including gluons of a particular colour. To illustrate this, consider the 3-partition on the far right of figure 2. We may represent this using the alternative notation of figure 3. Here, each separate colour in the 3-partition gets its own diagram, and there are thus 3 elements in the set. Furthermore, we can then dispense with the colours altogether, as these are no longer needed to separate out the elements of the partition.
The reason for using this alternative notation for partitions is that the diagrams are subject to a multiplication rule, which we have illustrated in figure 4. The rule for combining diagrams is as follows: to form the product, one takes the gluon part of each diagram, and draws them successively emanating outwards from the center of the product diagram. In figure 4(a), for example, the gluon joining external lines 1 and 2 (counting clockwise from top-left) appears closer to the center than the gluon joining external lines 2 and 3 in the product. Another example is shown in figure 4(b).
Here the first diagram on the left-hand side itself has two gluons. In the product diagram, this entire structure is reproduced, followed by the gluon which joins lines 1 and 3. One may also think of this rule as operating separately on the external lines: on each external line, one must order the gluon attachments outwards from the center of the diagram, according to the order in which they appear in the product of individual diagrams on the left-hand side. The reason for considering this multiplication rule is that it describes how the colour factor of a composite diagram is formed from the colour factors of its constitutents. That is, denoting the product of diagrams D 1 and D 2 via D 1 • D 2 , the colour factor of the product is given by Note that this multiplication rule for diagrams is non-commutative, corresponding to the fact that colour matrices associated with gluon emissions are themselves non-commuting. For example, the alternative ordering of the diagrams in figure 4(a) gives the product shown in figure 4(c), which is indeed not the same. Rather, it is a different diagram belonging to the same web.
Having defined the above multiplication rule, we can apply it within m-partitions of a given diagram, by multiplying the elements together. A given m-partition can be used to make m! diagrams (where all elements of the partition are used in each product), according to the number of permutations of the elements. Note that some of the diagrams one forms may be the same. Consider, for example, the 3-partition from figure 3. This can be used to make 6 diagrams, which must belong to the web shown in figure 5. Labelling the elements of figure 3 by [1,2,3], the diagrams that result from each permutation are shown in figure 6. Comparing this with figure 5, we see that diagrams these values is 2+2+1+1=3! as required. One may carry out a similar exercise for any m-partition of any diagram, and this leads us to define the following quantity, first introduced in ref. [88]: The "overlap function of D with respect to P D ", denoted by D |P D , is the number of ways that diagram D is made upon multiplying together all permutations of the elements of the partition P D of diagram D. Here D and D are diagrams in the same closed set (related by gluon permutations).
In the example of figure 6, one finds d|P = a|P = 2 and b|P = c|P = 1, where P is the 3-partition of diagram (a), shown in figure 3. Note that, in general, one has for any fixed partition P D , where n(P D ) is the number of elements in this partition. That is, every permutation of the elements of P D makes some diagram in the set to which D and D belong. Thus, summing over the overlap functions, by their definition, leads to the total number of permutations of the elements, which is n(P D )!.
Given the notion of an overlap function, one may show that the web mixing-matrix element R DD 7 is given by [88] There is a sum over all partitions P D of a given diagram D on the right-hand side, such that the left-hand side still depends on D (and indeed on D which appears in the overlap function in each term on the right-hand side). As already discussed in the introduction, this formula is purely combinatorial in nature, and forms a basis for the formal study of web-mixing matrices [94].
Some general properties of web-mixing matrices have already been established. Firstly, all webmixing matrices are idempotent: This means that web-mixing matrices are projection operators, whose eigenvalues are 0 and 1 (with associated multiplicities). This fact has been proven using statistical physics methods [88], but a combinatorial proof remains elusive. The physical interpretation of the idempotence property was examined in ref. [87] and then in detail in refs. [89,90]. Unit eigenvalues of a web-mixing matrix are associated with certain combinations of kinematic factors, each of which has an accompanying colour factor which corresponds to that of a fully connected gluon diagram [90]. Note in particular that the idempotence property implies that the rank of any web-mixing matrix is the same as its trace. The physical meaning of the rank is that it specifies the number of independent connected colour factors to which a given web contributes.
A second known property is that the row sums of web-mixing matrices are always zero: This was proven using combinatorial methods in ref. [88], and this property expresses the fact that the fully symmetric colour factor that can be obtained from a given web does not contribute to the exponent of the amplitude. Rather, it is generated in the amplitude itself by the exponentiation of lower-order contributions. Also in ref. [88], the sum rule was generalised for planar and non-planar subsets of diagrams in a given web. This provides an interesting link between the combinatorial properties of web-mixing matrices, and the 1/N c expansion, where N c is the number of colours.
Thirdly, all web-mixing matrices are conjectured to obey the weighted column sum rule: where s(D) is the number of ways of independently shrinking connected subdiagrams to the origin, and such that s(D) = 0 for diagrams that are not maximally reducible. For example, in the diagrams of figure 5(b) and (c), there is only one way in each graph of shrinking all three gluons independently to the origin, whereas in (a) and (d), there are two ways, depending on whether one shrinks the upper or lower gluon first. Thus, in this case the vector of multiplicity factors is given by (2,1,1,2). More examples of the s(D) factors were presented in ref. [88], and a justification of the column conjecture given for graphs which do not contain three or four gluon vertices. The physics of this result concerns the renormalisation of the hard interaction vertex at which the Wilson lines meet, and the physical justification of eq. (8) relied upon the fact that individual webs renormalise independently, with no mixing between webs. A combinatorial proof of eq. (8) would provide a basis for generalising the column sum rule, analogous to the sub-row sum rule introduced in ref. [88]. This would provide a link between subleading kinematic divergences, and subleading terms in the 1/N c expansion, as hinted at above. We will return later in the paper to the combinatorial problem of determining the vector s(D) for certain webs, and also see how the weighted column sum rule can be obtained combinatorially in particular cases.
In this section, we have reviewed a number of ideas relating to web mixing matrices, which pave the way for further combinatorial study of their general properties. In the following section, we begin such an investigation, by examining the rank of web-mixing matrices in particular cases.

The rank of web-mixing matrices
Our aim in this section is provide a simple formula for the rank of a given class of web-mixing matrices, and in doing so to provide a systematic classification of webs that is convenient for further study. As explained in the previous section, the rank of a web-mixing matrix has a simple physical meaning: it dictates how many independent connected colour factors a given web contributes to in the exponent of a Wilson line correlator. The general programme for calculating this exponent at a given order in perturbation theory is as follows [90]: 1. Pick a basis of connected colour factors, sufficient to span the space of all possible colour factors. This basis is not unique e.g. an alternative choice of colour factors can be generated by applying Jacobi identities. Reference [90] introduced an effective vertex formalism for calculating web graphs, which gave a convenient prescription for choosing a basis of colour factors.
2. Once a basis has been chosen, one must consider all webs at the given perturbative order. For each web, one must ascertain which connected colour factors it contributes to, and work out what the corresponding combinations of kinematic factors are that accompany each colour factor. This can be done either using web-mixing matrices, whose eigenvectors of unit eigenvalue are in one-to-one correspondence with connected colour factors, or by using the effective vertex formalism discussed in ref. [90].
3. The total contribution to a given colour factor can now be calculated, by carrying out the kinematic integrals found in the previous step. This will be a gauge invariant quantity, consisting of contributions from all webs which contribute to the given colour factor. We call such a quantity a gauge-invariant web, and its calculation may benefit from the use of cleverly designed gauges such as that presented in ref. [24].
A current priority is to calculate the relevant gauge-invariant web combinations at three-loop order, which will extend the state of the art of present-day resummation techniques by a further logarithmic order, leading to significant improvements in the precision of collider physics predictions.
To this end, all required combinations of kinematic factors have been presented, for a given colour factor basis, in ref. [90]. Calculation of the kinematic integrals is also underway [92,93], and involves the use of techniques that potentially generalise to higher orders in perturbation theory. This leads to the exciting possibility that one could calculate kinematic factors for given web families for any number of gluons, which strongly motivates seeking an all-order elucidation of web properties. One of the simpler -and most useful -properties to consider is the rank of web-mixing matrices. A general understanding of which connected colour factors come from which webs is a prerequisite for constructing gauge-invariant webs at a given perturbative order, and the rank provides a first step in this regard.
In investigating the above questions, it is convenient to classify webs according to how many Wilson lines they connect, as was already done in the three-loop analysis of ref. [90]. In general, webs may connect m lines at n-loop order, where 1 ≤ m ≤ n + 1. The case of m = 1 consists of pure selfenergy like contributions, which vanish for massless external lines. The case of m = 2 consists of diagrams which also enter the calculation of the traditional (two-loop) cusp anomalous dimension. The convenience of this classification is that webs can only mix (under gauge transformations) with other webs that connect the same number of lines. Thus, an independent basis of colour factors for m-line webs must be worked out for each value of m separately.
Consider first diagrams which connect n + 1 lines at n-loop order. The possible connected colour factors are shown for the case n = 3 in figure 7, where we use a notation in which the Wilson lines are separated from each other for ease of viewing each gluon diagram (in principle, these Wilson lines all originate from the same hard interaction point). In the connected colour factor, each Wilson line emits a single gluon only, as this is the only way in which the maximal number of lines can be connected. There are three such colour factors, only two of which are independent due to their being related by a Jacobi identity. The notation of figure 7 makes clear that one can view each connected colour factor as a 4-point scattering process, where the external gluons are attached to a Wilson line. Furthermore, each scattering process is at tree-level, which is a direct consequence of the fact that n + 1 lines are connected at n-loop order: one can only form loops by removing a gluon from an external line, such that at most n lines are connected. The number of independent colour factors at n-loop order then amounts to counting the number of independent colour factors in an (n + 1)-point tree-level scattering amplitude. This is known to be (n − 1)! (see e.g. ref. [95]).
For webs which link m < (n + 1) lines at n-loop order, the relevant connected colour factor graphs may have tree or loop-level topologies. That is, starting from the case of n + 1 lines above, every time one disconnects a given external line from the diagram, one may either connect to another line more than once, or create a loop. For example, we show the case of three lines at three-loop order in figure 8. Four possible connected colour factors are shown, three of which are tree-level graphs, and one of which is loop-level (and indeed can be seen to be proportional to the two-loop Y graph) 4 . For webs connecting progressively fewer lines at a given order, we can make the following general observation: Webs which span 1 ≤ m ≤ (n + 1) Wilson lines at n-loop order contribute colour factors corresponding to connected gluon graphs, whose topologies contain up to n + 1 − m loops.
For a given value of m, the question arises of how many independent connected colour factors receive contributions from a given web (namely, the rank of the web-mixing matrix). It is possible to provide a simple formula for the case m = n + 1, namely those webs which connect n + 1 lines with n single gluon emissons, such that there are no three or four-gluon vertices off the Wilson lines (one such example is the (1,2,2,1) web of figure 5): For webs in which n + 1 external lines are connected by n single gluon exchanges, the rank of the web-mixing matrix is given by where n i is the number of gluon emissions on line i.
To show this result, one may use the effective vertex formalism for webs introduced in ref. [90], and which was used there to prove the fully connected colour factor property. As explained in detail in that paper, products of Wilson line operators may be rewritten to involve ordinary (rather than path-ordered) exponentials, thus providing a natural setting in which to use replica arguments (see refs. [54,87]) to prove the connected colour factor property. This introduces an infinite tower of effective vertices, where each "k-vertex" describes the emission of k ≥ 1 gluons from a single Wilson line. Each such vertex consists of a superposition of (k − 1)! independent fully connected colour factors, with associated kinematic coefficients. In this language, diagrams which contribute to the exponent of a Wilson line correlator are produced by linking k-vertices on each Wilson line such that the overall diagram has a single connected piece.
As we have remarked above, the class of webs we are considering gives rise to connected colour factors in which there is at most a single gluon emission from each of the n + 1 external lines. This in turn means that there can be only one effective vertex on each line. As an example, consider the case of a (3,1,1,1) web. In figure 9(a), we show one possible diagram involving a 3-vertex on line 1. This results in a fully connected configuration of gluons. Alternatively, one may consider a 2-vertex and a 1-vertex on line 1 -which gives a disconnected gluon diagram, as shown in figure 9(b). Being disconnected, the latter does not contribute to the exponent.
Having established that for webs connecting n + 1 lines with n gluons each connected colour factor arises from diagrams in which the i th Wilson line contains a single n i -vertex, the total number of independent connected colour factors arising from each web is simply the product of the number of independent connected colour factors in each n i -vertex, which (from above) is (n i − 1)!. This gives the result of eq. (9), as desired. We will see special cases of this result later in the paper, and also connect it to the results of ref. [94].
When fewer than n + 1 lines are linked by gluon exchanges, it is no longer true that connected colour factor diagrams can only be formed by having single effective vertices on each line (plenty of examples can be found in ref. [90]). Nevertheless, it would be interesting to try to extend the above result, albeit after first solving for the total number of independent connected colour factors one expects for m-line webs with m < (n + 1). This appears to be a difficult combinatorial problem in general. Note that ref. [94] derived a result for the (2, 2, . . . , 2) web, namely that the rank of the mixing matrix for n gluon exchanges is given by (n + 1).
In this section, we have provided a simple formula for the rank of the mixing matrices for a special class of webs, namely those consisting of multiple gluon exchanges linking n + 1 Wilson lines at n-loop order. In fact, for certain of these webs one may go much further, and explicitly solve for the entire web-mixing matrix. To this end, and motivated by the pure mathematical study of ref. [94], it is useful to relate web diagrams to partially ordered sets (posets), whose properties we review in the following section.

Webs and posets
In section 2, we have reviewed what is known about web mixing matrices, with a particular emphasis on how they can be obtained combinatorially via overlap functions. In this section, we introduce mathematical ideas regarding partially ordered sets (posets). These were used in ref. [94] to study web mixing matrices, and will be reviewed here as a precursor for presenting the results of that paper. Furthermore, in what follows we will be able to use these poset ideas to obtain an explicit solution for the web-mixing matrix for (1, 2, 2, . . . , 2, 1) webs, and investigate the weighted column sum rule in this context. Given that posets are not necessarily well-known to a physics audience, we will summarise relevant facts about them here in order to make our presentation self-contained. A pedagogical introduction may be found in e.g. ref. [96].
A poset is a set of objects X endowed with an ordering operation ≤ possessing the following properties: This ordering operation applies to pairs of elements in the set X, but such that not all pairs are necessarily ordered with respect to each other (hence the "partially" in "partially ordered set").
A familiar example of a poset is the set of integers, where ≤ has its usual meaning of "less than or equal to". In this case, however, the set is fully ordered, with each element being ordered with respect to every other element.
The structure of a given poset can be represented on a Hasse diagram. In such diagrams, each element of X is represented by a symbol, and a line drawn between distinct elements x and y only if x ≤ y and {z ∈ X : x ≤ z ≤ y} = {x, y}. Then x is drawn below y to indicate the ordering. An example is shown in figure 10, for the set of elements {a, b, c}, such that a ≤ b, a ≤ c but no ordering is implied between elements b and c.
Returning to the web diagrams introduced in section 2, every web diagram D can be identified with a poset. The elements of the poset are the various irreducible subdiagrams that occur in the diagram D. The ordering operation ≤ is defined as follows: given two irreducible subdiagrams x and y, x ≤ y if x lies closer to the origin than y, in the sense that one cannot pull y into the origin without also pulling in x. As an example, consider the diagram shown in figure 11(a). In this diagram, the two gluons labelled A form an irreducible subdiagram: the fact that they are crossed means that one cannot shrink either gluon to the origin independently of the other one. The single gluon exchanges labelled B, C and D are also irreducible subdiagrams. Then the structure of the diagram is encoded by the information that B must be shrunk before A, C or D, and C must be shrunk before D. This corresponds to the Hasse diagram shown in figure 11(b).
Each diagram D ∈ W , for some web W , corresponds to a different Hasse diagram. One may then consider the set of Hasse diagrams instead of the original Feynman diagrams. Note that each Hasse For webs in which all irreducible subdiagrams are connected, there indeed exists a well-defined correspondence between colourings of the web diagrams D ∈ W , and colourings of the vertices of the corresponding Hasse diagrams. In this paper, we will be concerned with special cases of this, in which: (a) web diagrams consist exclusively of irreducible subdiagrams containing multiple gluon exchanges only; (b) any two Wilson lines have at most one gluon exchange between them. For the special cases we consider, the Hasse diagrams have a simple form, allowing a full solution of the relevant web-mixing matrices.
As well as the above notions relating to posets, we will also need the concept of a linear extension in what follows. A linear extension of a given poset P = (X, ≤) is defined as a permutation of the elements of X, such that all pairwise ordering relations remain satisfied. A simple example is provided by the poset whose Hasse diagram is shown in figure 10. This is specified by the set this ordering, shown in figure 12, and thus the number of linear extensions for this poset is equal to 2. In general the classification of linear extensions is complicated -it is at least known that the problem of counting the linear extensions of a poset is P -complete [97]. Given a poset P , we denote the set of linear extensions of P by L(P ). It can then be proven quite generally that the diagonal elements of any web-mixing matrix are given by ref. [94] Here π is a given linear extension of the poset corresponding to diagram D, and p the number of members in the set (i.e. the number of vertices in the corresponding Hasse diagram). As stated above, π corresponds to a permutation of the elements of the set, and one may label this by a permutation of [12 . . . p]. For example, we may label the two Hasse diagrams in figure 12 ( Note that the process of sequentially shrinking subdiagrams to the origin, discussed above from the combinatorial perspective, directly relates to the ultraviolet subdivergences present in these diagrams. This was observed already in refs. [87,91] and analysed in detail in ref. [89]. We will return to this point in section 6.
Having introduced various concepts relating to posets, and their relationship to web diagrams, we next consider the application of these ideas to specific web families. This is the subject of the following section.

General solutions for web-mixing matrices
In the previous section, we presented the notion of a partially ordered set, and explained how this concept relates to web diagrams. Methods from the theory of posets were used in ref. [94] to prove general properties of web-mixing matrices for two special cases, namely (1, 1, . . . , 1, n) and (1, 2, 2, . . . , 2, 1) webs, where we have used the notation introduced in section 2. An expression for the rank of mixing matrices for (2, 2, . . . 2) webs was also given. Our aim here is not to reproduce the proofs contained in ref. [94], to which we refer the reader for all relevant details. Rather, our goal in the present study is to use the poset ideas from the previous section to go further than ref. [94] in solving for the web-mixing matrices of interest explicitly, to all orders in perturbation theory.

5.1
The (1, 2, 2, . . . , 2, 1) web In this section we consider webs in which at least one, and at most two attachments occur on any given leg. One may choose to label the external lines such that the pattern of exchanges is labelled by (1, 2, . . . , 2, 1). This notation does not fully specify the family of webs of interest: we will further assume that all consecutive lines are connected by gluon exchanges. An example is shown in figure 13(a), for the case of seven gluons connecting eight Wilson lines, and where all external lines are shown emanating from the same hard interaction vertex. It is also convenient to use an alternative representation of these diagrams, in which all external lines are shown parallel to each other. That is, one may redraw the diagram of figure 13(a) as shown in figure 13(b). Shrinking gluons to the origin in the former amounts to sinking gluons in the latter, and the hard interaction can be thought of as a horizontal line at the bottom of the diagram.
The (1, 2, 2, . . . 1) web family is of physical interest because it connects n + 1 external lines with n gluons. This is the maximum possible number of external lines which can be entangled at a given loop order, and corresponds to the class of webs whose rank was discussed in section 3. A given number of gluon exchanges n specifies a unique (1, 2, 2, . . . 1) web (other similar webs at a given loop order can be obtained by permuting the external lines). The number of diagrams in this web is equal to the dimension of the permutation group which acts on the gluon emission vertices, namely We collect explicit examples from n = 2 to n = 4 in appendix A, together with the relevant mixing matrices.  Figure 14: (a) The (1,2,2,1) web redrawn as in figure 13(b); (b) corresponding Hasse diagrams.
As stated above, this family has been considered in ref. [94], where some general properties of the web-mixing matrices for all n have been proven. Of relevance for the current study is that, for these webs, one has Tr R (1,2,...,2,1) = 1.
Given the rank is equal to the trace for an idempotent matrix, this provides a special case of eq. (9), which dictates that the mixing matrix of any (1, 2, 2, . . . , 2, 1) web has unit rank. Thus, there is only one linearly independent row or column. Here we use this property, together with the structure of the Hasse diagrams corresponding to each web diagram, to write down a solution for the web-mixing matrix for any n, extending the results of ref. [94].
The first thing to do is to relate the diagrams in a (1, 2, 2, , . . . , 2, 1) web to Hasse diagrams. We can do this by drawing the diagrams in the form shown in figure 13(b), and then associating a vertex • with each gluon exchange. This is possible because all irreducible subdiagrams correspond to a single gluon exchange. An example is shown in figure 14 for the (1,2,2,1) web. Note that all the Hasse diagrams correspond to a single chain, which may or may not have kinks in. This is due to the special nature of the (1, 2, . . . , 2, 1) web, namely that each gluon exchange has only two neighbouring gluons (on adjacent external lines), and no other gluons exchanged between its own pair of external lines. This allows us to immediately write down the Hasse diagrams for any number of gluons n: one simply draws all possible kink assignments in a chain with n − 1 links. As a further example, the case of n = 4 is shown in figure 15.
The special form of these Hasse diagrams provides a basis for further study of the web-mixing matrix, and also the shrinking factors s(D) which enter the weighted column sum rule. In order to investigate the former, we need to be able to evaluate the overlap functions introduced here in section 2. That is, the number of ways of forming a diagram in a (1, 2, . . . , 2, 1) web from a given partition of some other diagram in the same web. To this end, one may first see what partitions look like in terms of Hasse diagrams. As discussed in section 2, a partition corresponds to a colouring of the connected subdiagrams in a given web diagram. In the present case, this amounts to a colouring of gluons, namely of the vertices of the Hasse diagram. A given m-partition separates the Hasse  diagram into m groups of vertices. Furthermore, the multiplication rule for web diagrams easily translates into a multiplication rule for the elements of a partitioned Hasse diagram. Examples with n = 4 are shown in figure 16, which shows two different 2-partitions of the same web diagram, and the diagrams obtained by taking products of the partition elements. One sees that the effect of the colouring is to break links between consecutive vertices if these have different colours; otherwise links are left intact. In each product, the grey element of the partition can be placed above or below the black element, due to the fact that the product of the related Feynman diagrams corresponds to placing subgraphs further or nearer to the hard interaction. It is straightforward to generalise this example to higher numbers of gluons (n), and number of colours (m).
We can now use these ideas to solve explicitly for a particular column of the web-mixing matrix. We start by considering the Hasse diagram consisting of a chain with no kinks, such that each vertex is above the vertex to its left. For the (1,2,2,1) case, this is the diagram in the upper-right of figure 14(b); for n = 4, this is the first diagram in figure 15. It is clear that such a diagram exists for any value of n, and we denote this diagram generically by ↑. It is then possible to find exactly the overlap functions ↑ |P D , namely the number of ways of making diagram ↑ from an arbitrary partition P D of diagram D in the same web. One may then, in turn, solve for the values R D↑ , namely the entire column of the web matrix corresponding to the diagram ↑.  figure 15, and the diagrams formed by taking products of the partition elements.
The first step is to note that if ↑ can be made from the partition P D , there is only one possible way to do this. This is because the various vertices of the Hasse diagram must be multiplied together in the unique sequence which generates the ascending pattern of ↑. Thus, in general the overlap function ↑ |P D can only be one or zero. Furthermore, it is clear that ↑ cannot be made out of partitions in which vertices of the same colour are separated by one or more vertices of a different colour. An example of this is shown in the upper diagram of figure 16, which has two grey vertices separated by a black vertex. In multiplying the partition elements together, the two grey vertices must occur at the same level, and either above or below the black vertices. This rules out the possibility of obtaining the diagram ↑, and this argument clearly generalises to higher numbers of gluons and colours.
It follows that ↑ |P D can only be non-zero for partitions in which groups of consecutive vertices are painted the same colour. Furthermore, in order to make ↑, none of these groups can contain a descending link. An example is shown in figure 17, in which we consider 2-partitions of the upper-rightmost diagram in figure 15. In the upper partition, a descending link is preserved by the colouring, which must necessarily survive in the diagrams formed by taking products of the partition elements. In the lower partition, however, the descending link is broken by choosing two different colours on either side of the link. The partition elements then consist solely of groups of ascending vertices, and it is then obvious that one can form ↑ from such a partition, as shown explicitly in figure 17.
Let us now consider a diagram D whose Hasse diagram has r descending links, and an m-partition of this diagram (equivalently, a colouring of D with m colours). In order to make ↑, the partition must be such that all descending links in the Hasse diagram are removed. This requires at least r + 1 colours, such that we must have m > r for ↑ |P D to be non-zero. After removing the descending links, any assignment of the remaining colours amongst the ascending links will give a partition which can reproduce ↑. Each of these assignments corresponds to a separate partition P D , and the number of such partitions amounts to the number of ways of choosing m − 1 − r links (the number of remaining colours) out of n − 1 − r (the number of ascending links, which is the total number of links (n − 1) minus r descending ones). Thus, one has P D ,n(P D )=m where Θ(n) is the Heaviside function, such that Θ(0) = 1. Substituting this result into eq. (5), one finds where we have replaced the sum over individual partitions with a sum over the number of elements in a partition m ≡ n(P D ). One may now set p = n − 1 − r (the number of ascending links in the Hasse diagram for D) to get where we have set s = n − m and used the symmetry property of binomial coefficients We may now use the rather quirky identity (see e.g. eq. (5) of ref. [98]) to get This is an important formula: it tells us the value of any element of the column of the web-mixing matrix corresponding to diagram ↑, for any number of gluons! Furthermore, the only important quantity is seen to be the number of ascending links p in the Hasse diagram of diagram D. Regarding the column as a whole, we can write a simple expression for the multiplicity n p of the element given by eq. (17). This will simply be the number of Hasse diagrams with p ascending links. Given that there are (n − 1) links in each Hasse diagram in total, this corresponds to the number of ways of choosing p links from (n − 1), namely A consistency check of this result is that the sum over all multiplicities is

20
where we have used a well-known summation formula for binomial coefficients (see e.g. ref. [96]). The right-hand side is the dimension of the web-mixing matrix for the (1, 2, . . . , 2, 1) web with n gluon exchanges, as it should be.
We can now go further this. The fact that all columns are proportional to the column R D↑ (as follows from eq. (11)) means that we may write the elements of the mixing matrix as where a E are coefficients to be determined. A full solution for the web-mixing matrix now corresponds to finding the latter. To this end, we may use the result of ref. [90] -namely, that left-eigenvectors of the web-mixing matrix are in one-to-one correspondence with connected colour factors. The (1, 2, 2, . . . , 2, 1) web is special in that there is only one such connected colour factor [94], and it is straightforward to see how this arises diagramatically: the only way to form a fully connected colour factor is to form a commutator of the colour generators on each line which has two gluon emissions 5 .
One may show that the left-eigenvector c D with unit eigenvalue of the web-mixing matrix is given by: where p D is the number of ascending links in the Hasse diagram corresponding to D. To see this, one may use eqs. (17), (20) and (21) to find We may convert the sum over diagrams D into a sum over the number of ascending links p in the corresponding Hasse diagrams, to get In the second step, we have used the result for the multiplicity factor n p from eq. (18). Equation (23) tells us that c D of eq. (21) is the left eigenvector we are looking for, provided Another way to describe this left eigenvector is as follows: The left eigenvector for the web-mixing matrix for a (1, 2, . . . , 2, 1) web is proportional to c E = (−1) t(E) , where t(E) is the number of permutations of gluon emissions that takes diagram ↑ to diagram E. This is consistent with the fact that the single connected colour factor for this web has a commutator of gluon generators on each external line that has two gluon emissions: interchanging any pair of gluon emissions leads to a minus sign. An example is the (1,2,2,1) case of figure 5: where Not only has the above analysis determined the left eigenvector of the web-mixing matrix with unit eigenvalue -it has also fixed the coefficients a E in eq. (20). We have thus explicitly obtained a general solution for the web-mixing matrix: where, as above, p D is the number of ascending links in the Hasse diagram of diagram D. We may check that this solution satisfies a number of constraints. Firstly, there is that of idempotence of the web-mixing matrix. The squared web-mixing matrix is given by and we may simplify the right-hand side by observing Inserting this on the right-hand side of eq. (28) yields as required. The zero sum row rule is also satisfied, given that as follows from eq. (21). The reader may easily verify that the general solution reproduces the mixing matrices for the cases n = 2, n = 3 and n = 4, which are collected in appendix A. For higher values of n, we may also compare eq. (27) with the general combinatorial formula for diagonal elements of web-mixing matrices, presented in eq. (10). A precise comparison at all orders requires a more involved understanding of linear extensions of the kinked chain than has been presented here (although we return to this point in section 6). Nevertheless, one may note that eq. (10) also involves the inverse binomial coefficients.
One may also check that our solution satisfies the weighted column sum rule of eq. (8). However, this is not so straightforward given that one must first elucidate the structure of the s(D) values for each web. We postpone a discussion of this point to section 6.
As will be detailed elsewhere [92], it is possible to develop systematic methods for carrying out the integration of the kinematic factors for a general number of gluon exchanges. Combined with the above results, this would provide a complete solution of this web family to all orders in perturbation theory, and strongly motivates the investigation of more general webs. (1, 1, 1, . . . , 1, n) web

The
In the previous section, we have discussed the (1, 2, 2, . . . , 2, 1) web at length, and given a full solution for the web-mixing matrix for an arbitrary number of gluon exchanges. In this section,  1,1,1,3) web.
we present a solution for a different web family, namely that of (1, 1, 1, . . . , n). This is in fact a shorter exercise than the previous case, as most of the necessary results have already been derived in ref. [94].
The Hasse diagrams for the (1, 1, . . . , 1, n) case consist of n! vertical chains, each corresponding to a permutation of the gluon emissions on the line which has n gluon emissions. An example is shown, for the n = 3 case, in figure 18. Consequently, one can label each diagram by a permutation π of [12 . . . n]. A general element of the web-mixing matrix can then be denoted R DD ≡ R πσ , where π and σ are the permutations corresponding to diagrams D and D respectively. We then have the following result 6 : The web-mixing matrix element R πσ is given by where r πσ is the number of ascents in the permutation that takes σ to π.
The definition of the latter quantity, and how to find the web-mixing matrix element, is best clarified by a simple example. To this end, we consider the web of figure 21, whose web-mixing matrix is shown in eq. (50). The ordering of diagrams shown in the figure can be labelled by the permutations 123, 132, 312, 321, 213 and 231, where these denote the ordering of the gluons (moving outwards from the hard interaction) on the Wilson line which has three gluon emissions. Consider now the (3,2) element of the matrix, corresponding to R (312),(132) . The permutation that takes 312 to 132 has 1 ascent (in other words, one consecutive pair of increasing numbers; to see this, note that it is the same permutation as that which takes 123 to 213), and thus from eq. (37) one expects as indeed is observed. As a second example, consider the (2,6) element of the matrix, corresponding to R (132,231) . The permutation which takes 132 to 231 would take 123 to 321, and thus has 0 ascents. Thus one expects as observed.
Some comments are in order. First, the web-mixing matrices for this family are symmetric i.e. R πσ = R σπ . This follows from the definition of r πσ , and the fact that the number of ascents in the permutation σ → π is the same as the number of descents in the inverse permutation π → σ. That is, Swapping r πσ → r σπ thus leaves the right-hand side of eq. (37) invariant, via the reflection property of binomal coefficients. Secondly, there is a remarkable similarity between the expression of eq. (37), and the formula for the possible values in the (1, 2, 2, . . . , 1) case, eq. (27). This is because the Hasse diagrams in the present case are governed by permutations of [12. . . . n], where these permutations can be classified by the number of ascents and descents, and thus mapped to kinked chain diagrams. The multiplicities of the various values in this case, however, are different, as we discuss below.
Equation (37) has a number of corollaries. Firstly, the permutation that takes a diagram to itself has n − 1 ascents (i.e. it takes [12 . . . n] to itself). Thus, the diagonal elements are given by This can clearly be observed in the examples in appendix A. It is also possible to derive this result directly from eq. (10), which gives a general formula for the diagonal elements of any web-mixing matrix (in our case, the number of vertices in the Hasse diagram is p = n). The poset for each web diagram has a fully ordered Hasse diagram (e.g. figure 18), thus has only one linear extension. Furthermore, the number of descents in the permutation which takes a diagram to itself is zero, and thus eq. (39) must correspond to eq. (10) with des(π) = 0. This is indeed the case.
Knowledge of the diagonal elements allows us to obtain a simple solution for the trace: Tr R (1,1,...,1,n) = (n − 1)!, which follows from the fact that each diagonal element is 1/n, and there are n! of them. Again the trace corresponds to the rank, and thus one derives a special case of the general result of eq. (9).
Having solved for the possible elements of the web-mixing matrix, it is also convenient to note their multiplicities in each row and column. These are given by the following result [94]: Identifying r ≡ r πσ , the values in a given row or column of the web-mixing matrix are given by eq. (37) with all possible values of r. The multiplicity of each value is given by the Eulerian number n r .
Given that the Eulerian numbers may not be generally familiar to the reader, we briefly introduce their properties in appendix B. As clarification of their role as multiplicities in web-mixing matrices, consider the n = 4 case, whose web-mixing matrix is shown in eq. (51). In any given row or column, the possible values are given by eq. (37) with r = 0, 1, 2, 3, which gives -1/4, 1/12, -1/12 and 1/4 respectively. Furthermore, these should have multiplicities 1, 11, 11 and 1, as given by the relevant Eulerian numbers from table 2. This is indeed observed.
Note that the sum of multiplicities in any given row or column is where we have used a known identity of Eulerian numbers (in words: summing over r the number of permutations with r ascents must give the total number of permutations, which is n!). The right-hand side of eq. (41) is the dimension of the web-mixing matrix, as it should be. Note that the sum over elements in any given row is where we have used eq. (52). This is an explicit realisation of the zero sum row rule of eq. (7), which has been proven to be true fully generally in ref. [88]. A similar sum also ensures that for the (1, 1, . . . , 1, n) web family, the weighted column sum rule is satisfied. We return to this point in section 6.
In this section, we have presented solutions for two specific families of web-mixing matrices, for any number of gluon exchanges, extending the results of ref. [94]. Next, we consider the weighted column sum rule of ref. [89] in the poset language.

The weighted column sum rule
In this section, we return to a specific combinatorial property of web matrices, namely the weighted column sum rule of eq. (8). This was first conjectured in ref. [89], and is crucially related to the renormalisation of the hard interaction vertex at which the external eikonal lines meet. Based on this relation, it was shown in ref. [89] that for webs which consist solely of multiple gluon exchanges, the sum rule must hold. However, the sum rule appears to be more general than this, and was indeed satisfied in all examples considered in ref. [89], including those webs which contain multiple gluon vertices off the eikonal lines. This perhaps is not surprising, given the results of ref. [90], which show that webs containing multiple gluon exchanges are intimately linked with webs which have multiple gluon vertices, by virtue of sharing the same connected colour factors. They are also linked by gauge transformations.
Given the above remarks, it is useful to study the weighted column sum rule from a purely combinatorial point of view. Whilst a full combinatorial proof for arbitrary webs remains elusive, we will relate some key concepts which arise from the poset language of section 4, and also examine the sum rule for the particular cases of (1, 2, . . . , 2, 1) and (1, 1, . . . , n) webs.

s(D) and linear extensions
Our first step is to tighten up the combinatorial definition of the number of ways of shrinking irreducible subdiagrams to the origin in a given web diagram. This may be stated as follows, for 26 all classes of web: Let P be the poset corresponding to a given web diagram D. Then the number of ways of shrinking all irreducible subdiagrams to the origin is the number of linear extensions of P .
This follows from the definition of a linear extension which, as discussed in section 4, corresponds to a permutation of the elements of a poset which preserves the ordering requirements. The Hasse diagram for a web diagram D containing m irreducible subdiagrams contains m vertices, one for each subdiagram. Each possible shrinking sequence must preserve the ordering of these vertices, and corresponds to a permutation of the vertices of the Hasse diagram. Thus, each shrinking sequence can be uniquely identified with a linear extension of the poset of D, such that the total number of shrinking sequences is equal to the number of linear extensions. As an example, consider the diagram of figure 11(a), whose Hasse diagram is shown in figure 11 This is easily seen from the fact that the Hasse diagram for any (1, 1, . . . , 1, n) web diagram is fully ordered (e.g. figure 18), and thus all gluon exchanges can only be shrunk in a single sequence, beginning with the gluon closest to the hard interaction, and working outwards.
where we have used eq. (52). Note that this is analagous to the proof of the zero sum row property for these webs, as performed in eq. (42). This is due to the fact that the same values occur with the same multiplicities in any given row or column of the web-mixing matrix, and also crucially depends on the fact that s(D) is the same for all diagrams.
The fact that the row and column rules appear to be related for this web family (a consequence in this special case of symmetry of the web-mixing matrix) is interesting by itself. The row rule is known to be related to colour properties of web diagrams, including the 1/N c expansion [88].
On the other hand, the column rule is related to kinematic information, namely the structure of singularities [89]. A link between the row and column rules is then a form of colour-kinematic duality, which may be related to other such dualities in the literature. It would be interesting to investigate the possible relationship between the row and column rules in more detail and for more general webs. Another possible question is that of whether there is a sub-column rule analagous the sub-row rule introduced in ref. [89] via an expansion about the planar limit.
Here we have seen explicitly that the weighted column sum rule is satisfied, for this particular family, at the level of individual webs. Furthermore, this follows purely from combinatorial considerations. This is consistent with the fact that webs renormalise individually, with no mixing between webs. Next, we consider the (1, 2, . . . , 2, 1) case.
6.3 Column sum rule for (1, 2, . . . , 2, 1) webs This web family makes an interesting counterpoint to that of the previous section. It is nontrivial in that the s(D) values for each web are different for different diagrams D. This in turn implies that the weighted column sum rule is different to the zero sum row rule.
An explicit solution for a single column of the web-mixing matrix for the (1, 2, . . . , 2, 1) family (corresponding to diagram ↑) has been given in eq. (17), where each value occurs with a multiplicity given by the binomial coefficient of eq. (18). All other columns are proportional to this column as a consequence of eq. (11) (which together with idempotence implies that all web-mixing matrices for this family have unit rank). Hence, if the column sum rule is satisfied for the column R E↑ , it is satisfied for all columns.
As stated in section 6.1, s(D) for a maximally reducible web diagram is the number of linear extensions of the poset corresponding to diagram D. In the present case, this corresponds to linear extensions of posets whose Hasse diagrams correspond to kinked chains (see e.g. figure 15). We collect the corresponding s(D) values up to nine loop order, together with their multiplicities, in appendix C. As we describe in the appendix, a closed form solution for these values is not known. However, one may check the weighted column sum rule in particular cases, for the column R D↑ . One must first identify which s(D) values correspond to which diagrams. The value of R D↑ for each diagram is then given by eq. (17). As an example, the s(D) values for the n = 3 case and corresponding to the diagrams in figure 14 are 2, 1, 1 and 2. The corresponding elements of R D↑ are -1/6, +1/3, +1/3 and -1/6. One then finds D s(D)R D↑ = 0 as required.
As a highly non-trivial example, one may consider also the case n = 5. The Hasse diagrams corresponding to the 16 web diagrams are shown in figure 19, and the corresponding s(D) factors and R D↑ values in table 1. Note that we have multiplied the web-mixing matrices by 5! = 120 in order to make the numbers easier to read. Contracting the s(D) values with the web-mixing  = 0, as expected if the weighted column sum rule is true. One sees from this particular example that the nature of the cancellations involved is highly non-trivial, involving an intimate interplay between the numbers of eq. (17), and the linear extensions of the kinked chain.
A dramatic simplification occurs for even values of n. Consider a diagram D and its conjugate diagram D * whose Hasse diagram is obtained from that of D by reflection about the horizontal. The shrinking factors for these two diagrams are related by which follows from the fact that the shrinking sequences for diagram D * can be obtained by reversing those for D. The Hasse diagram of D * will have n − 1 − p ascending links if that of D has p ascending links, and thus from eq. (21) one may write Hence, for even n values, conjugate diagrams (which have the same s(D) values) have column entries in the web-mixing matrix which are equal in magnitude but opposite in sign. The weighted column sum rule is then automatically satisfied, and thus proven for the even (1, 2, 2, . . . , 2, 1) webs.
In this section, we have expressed the factors s(D) occuring in the weighted column sum rule of eq. (8) in poset language, namely as the number of linear extensions of the poset corresponding to web diagram D, where this is maximally reducible. In the (1, 1, . . . , 1, n) case, we see that the sum rule can indeed be proven combinatorially. The (1, 2, 2, . . . , 2, 1) case is more intricate, and involves combinatorial quantities (linear extensions of kinked chain posets) about which little is currently known. However, we have proven that the sum rule is satisfied for individual webs, for an even number of gluon exchanges. Whilst a formal proof of the weighted column sum rule (also including other web families) remains elusive, we hope that the ideas of this section will be useful in the further study of this property.

Conclusion
The study of infrared singularities in non-Abelian gauge theories continues to be of interest, not least as this is intimately linked to the improvement of collider physics predictions through resummation of large logarithms. In this paper we have continued the development of a diagrammatic ("web") approach to soft gluon exponentiation or, equivalently, correlators of multiple Wilson lines. The web approach to multiparton scattering identifies closed sets of diagrams, whose kinematic and colour parts mix according to web-mixing matrices. This systematic organisation greatly simplifies the calculation of multiloop contributions, including the structure of subdivergences. A full understanding of the exponent necessitates the development of methods for carrying out the necessary kinematic integrals -and these are underway [92,93] -but also the classification of the mixing matrices themselves. Here we have focused on the latter. Their pure combinatorial origin, already identified in ref. [88], has allowed us to relate the study of web-mixing matrices to that of partially ordered sets, a connection which was already explored in a pure mathematical context in ref. [94].
For a special class of webs, namely those consisting of n + 1 Wilson lines linked by n gluon exchanges, we have presented a simple formula for the rank of the mixing matrix. This tells us how many independent colour factors receive contributions from a given web, and it is an interesting combinatorial problem in its own right to extend this to more general webs.
Using the poset language of ref. [94], we have gone much further than studying the rank in being able to provide explicit solutions for two particular web families, namely the (1, 2, 2, . . . , 2, 1) and (1, 1, . . . , 1, n) webs. It is interesting to note that these are extremal cases of the special class mentioned above -they contribute 1 and (n − 1)! independent colour factors respectively, where the total number of colour factors at a given order is known to be (n − 1)!. The formulae for the mixing matrices involve intriguing combinations of inverse binomial coefficients, matching the similar behaviour observed for diagonal elements of more general webs in ref. [94]. Combined with general methods for carrying out kinematic integrals, details of which will be presented elsewhere [92,93], these results are potentially very powerful, allowing unprecedented insights into the all-order structure of the exponents of amplitudes. This provides strong motivation for generalising the results presented here, particularly to webs involving three and four-gluon vertices off the Wilson lines. However, it is also worth stressing that there may be gauges [24] in which the effect of some of these graphs can be absorbed in webs consisting of multiple gluon exchanges only, making the web-mixing matrices for the latter more physically relevant.
The poset language also allows some understanding of the weighted column sum rule for webmixing matrices presented in ref. [89], which implements the cancellation of subdivergences. While a general combinatorial proof of this property remains elusive, the poset language seems to provide the right tools for such a proof. A deeper understanding of the columns of web-mixing matrices may provide insights into a perceived duality between subleading colour corrections, and subleading kinematic divergences, given the zero sub-row sum rule observed in ref. [88].
Further use of the poset langauge presented here and in ref. [94] (and possibly other combinatorial methods) can be used to generalise the results to other web families. Our ultimate aim is to probe the all-order structure of infrared singularities and the associated physical consequences, by elucidating where possible the structure of web-mixing matrices on the one hand, and kinematic integrals on the other. Progress on both of these fronts is ongoing.

A Examples of web-mixing matrices
In this appendix, we collect explicit examples of web-mixing matrices for the (1, 2, 2, . . . 2, 1) and (1, 1, . . . , 1, n) family, where n is the number of gluon exchanges in each case.

2, 1) webs
The web diagrams corresponding to the n = 2 and n = 3 cases are shown in figures 1 and 5 respectively. The corresponding web mixing matrices are: R (1,2,1) = 1 2 and R (1,2,2,1) = 1 6 The (1,2,2,2,1) web is shown in figure 20, where we label each diagram according to a notation first introduced in ref. [87], and in which gluon exchanges associated with the same connected subdiagram are labelled by integers i. The m-tuple [n 1 , n 2 . . . , n m ] associated with each external line then labels the gluon emissions on that line, working inwards towards the hard interaction vertex. In this case the mixing matrix is Figure 21: Web whose mixing matrix is given by eq. (50).
A.2 (1, 1, 1, . . . 1, n) webs The n = 2 case can be obtained by relabelling external lines in the (1, 2, 1) web shown in figure 1, whose mixing matrix is given in eq. (47). The n = 3 case is shown in figure 21, and the mixing matrix is with diagrams ordered as in the figure.
For the n = 4 case (1, 1, 1, 1, 4), we do not reproduce the 24 web diagrams here. We may label each diagram concisely by a permutation of (1234), where this denotes the ordering of the gluons on the parton line with 4 attachments (away from the hard interaction), such that gluon i joins this line with external line i. One then has the mixing matrix where the ordering of the diagrams is given by

B Eulerian numbers
In this appendix, we briefly summarise the properties of Eulerian numbers, which occur in the solution of the web-mixing matrix for (1, 1, . . . , 1, n) webs.
The Eulerian number n r counts the number of permutations of n objects with r ascents (consecutive pairs of rising numbers), where 0 ≤ r ≤ n − 1. Results for the first few values of n and r are shown in table 2.
A number of identities regarding Eulerian numbers are known. For the present paper, the following