Six-Gluon Amplitudes in Planar ${\cal N}=4$ Super-Yang-Mills Theory at Six and Seven Loops

We compute the six-particle maximally-helicity-violating (MHV) and next-to-MHV (NMHV) amplitudes in planar maximally supersymmetric Yang-Mills theory through seven loops and six loops, respectively, as an application of the extended Steinmann relations and using the cosmic Galois coaction principle. Starting from a minimal space of functions constructed using these principles, we identify the amplitude by matching its symmetries and predicted behavior in various kinematic limits. Through five loops, the MHV and NMHV amplitudes are uniquely determined using only the multi-Regge and leading collinear limits. Beyond five loops, the MHV amplitude requires additional data from the kinematic expansion around the collinear limit, which we obtain from the Pentagon Operator Product Expansion, and in particular from its single-gluon bound state contribution. We study the MHV amplitude in the self-crossing limit, where its singular terms agree with previous predictions. Analyzing and plotting the amplitudes along various kinematical lines, we continue to find remarkable stability between loop orders.

These DCI contributions to the amplitude have especially nice properties. In particular, they have been observed to have uniform transcendental weight equal to twice the loop order. While generic amplitudes in this theory may depend on elliptic polylogarithms [44,47,48] and even more complicated [49,50] functions, maximally-helicity-violating (MHV) and nextto-MHV (NMHV) amplitudes are expected [10] to contain only functions drawn from the class of generalized polylogarithms. This class of functions is well understood, and has a Hopf algebra coaction structure that has been exploited to great effect [17,[51][52][53][54][55][56]. (We will often refer to this coaction loosely as a coproduct, even though technically it is not because the two objects a function is mapped to under it are really of different types.) Maximally iterating this coproduct yields an object called the symbol of a function [17]. In six-and seven-particle kinematics (and for higher multiplicity two-loop MHV amplitudes) the letters entering the symbol have intriguing connections to cluster algebras [22,23,25,28,[57][58][59][60][61].
Using these known symbol letters in six-particle kinematics and considering functions with physical branch cuts, it is possible to write down the space of functions within which the amplitude should reside. Moreover, by requiring a general ansatz of these functions to have the right symmetries and kinematic properties, the amplitude can be identified uniquely within this space. This bootstrap procedure was first employed for the three-loop six-particle MHV amplitude [19,24]. Since then, it has been used to fix six-particle amplitudes through five loops [20,26,27,29,31,33] and seven-particle amplitudes through four loops [30,34,62].
Crucial to this progress has been an increasingly refined understanding of the space of functions needed to represent these amplitudes. In particular, the hexagon and heptagon bootstraps have in recent years been augmented by knowledge of the Steinmann relations [33,34,[63][64][65]. In a companion paper [66], we present the next (and probably final) evolution of this space. We find that the Steinmann relations can be applied at any depth in the symbol [45,66]. These conditions apply at the level of the symbol, but there are also important restrictions on the function space that involve multiple zeta value (MZV) constants, which are invisible at the level of the symbol. Once the amplitudes are properly normalized, they belong to a space of functions that satisfies a cosmic Galois coaction principle [67][68][69][70], which provides important restrictions on the MZVs that can appear. Furthermore, the space of functions is minimal through weight seven: it is impossible to eliminate any more functions, because all of them are needed to capture the iterated derivatives (or coproduct entries) of the six and seven loop amplitudes. We call this minimal space of functions H hex .
In this paper, we employ this space of functions to construct the six-particle amplitude through seven loops in the MHV sector and through six loops in the NMHV sector. After constructing an ansatz with the appropriate symmetries, we require that it behaves appropriately under the action of the dual superconformalQ operator [21]-which constrains the symbol letters that can appear in the final entry-and we require the expected leading-power behavior in the collinear limit. To do this, we construct a basis of functions iteratively in the weight, and work out their limiting behavior. At six and seven loops, these constraints leave fewer than 20 free parameters, out of thousands of initial parameters.
To fix the remaining parameters, we consider additional kinematic limits that supply independent information about the amplitude. We first consider the multi-Regge limit, in which outgoing particles are widely separated in rapidity, and the amplitude factorizes in a Fourier-Mellin transformed space [29,40,[71][72][73]. The behavior in this limit is now understood to all orders via integrability techniques [74]. Constraints from this limit yield a unique answer for the amplitude through five loops and for the six-loop NMHV amplitude. However, at both six and seven loops for the MHV amplitude we observe a novel qualitative feature: the appearance of a single function that vanishes in the leading-power collinear and multi-Regge limits. Hence, the free parameter multiplying this function cannot be determined from these limits alone.
To determine this one remaining coefficient in our ansatz, we use a kinematic expansion around the collinear limit, which can be computed in the framework of the Wilson loop (or Pentagon) Operator Product Expansion (OPE) [75][76][77][78][79][80][81][82][83][84][85][86][87]. Quite interestingly, this last step requires going up to second order in the near-collinear expansion, and specifically examining the OPE contribution of the first gluon bound state.
Having determined the amplitudes, we proceed to study their properties, which reveals another novel feature that first appears at six loops: a previously conjectured cross-loop-order relation between MHV amplitudes and NMHV amplitudes [29,33], which held for several loop orders, no longer seems to hold. As part of our analytic study, we present formulas for the values of our amplitudes at various points in the space of kinematics, in order to exhibit their number-theoretic properties. This study includes the limit of physical 2 → 4 and 3 → 3 kinematics where the corresponding Wilson loop approaches a self-crossing configuration [88][89][90][91], and the singular terms can be resummed to arbitrary loop order. We also plot the perturbative amplitudes numerically along various lines in the Euclidean region. The new analytic features we described do not affect a remarkable numerical consistency between loop orders at generic values of the cross ratios. We interpret this consistency as evidence that the perturbative expansion of the amplitudes in these regions has a finite radius of convergence.
This paper is organized as follows. In section 2, we begin by reviewing our minimal space of functions, and discuss how to normalize the six-particle MHV and NMHV amplitudes in order to fit them into this space. Then, in section 3, we describe how to construct the amplitude through seven loops in the MHV sector and at six loops in the NMHV sector, by applying constraints from symmetries and kinematic limits. Section 4 contains the numbertheoretic and numeric exploration of the amplitudes on a variety of kinematic points and lines. In section 5 we examine the self-crossing limit of the MHV amplitude. In section 6 we conclude and discuss directions for future research.
We provide the following ancillary files along with this paper: SixGluonAmpsAndCops.m, SixGluonHPLLines.m, SelfCross.m, SelfCrossSingular.m, hexMRKL1-7.m, WL0-6.m, W1111L0-6.m, WL7.m, and WLOPEblocks.m. These computer-readable files describe results that are too lengthy to place in the text, including (respectively) a coproduct description of the amplitudes, their values on certain Euclidean lines in terms of harmonic polylogarithms [92], the MHV values in the self-crossing limit, a formula for the singular self-crossing behavior to 20 loops, formulas for the amplitudes in the multi-Regge limit, formulas for the near-collinear limit of the framed Wilson loop through six loops for MHV and NMHV and at seven loops for MHV, and individual contributions to the Wilson loop OPE. The files are hosted at [93].

Superamplitudes, kinematic variables and generalized polylogarithms
Let us begin by briefly recalling certain facts about the general structure of six-particle amplitudes in planar N = 4 SYM theory, including their kinematic dependence and the class of functions that encompasses them. Due to the supersymmetry of this theory, we may combine color-ordered amplitudes with different external particles into corresponding superamplitudes, which are defined in an on-shell superspace [94][95][96][97]. The superfield Φ can be written in terms of Grassmann variables η A that transform in the fundamental representation of the SU (4) R-symmetry, where the gluons of each helicity G ± , gluinos of each helicity Γ A and Γ A , and scalars S AB represent the on-shell particle content of the theory. The superamplitude A n (Φ 1 , Φ 2 , . . . , Φ n ) is typically broken into three factors: the BDS ansatz A BDS n [35], the remainder function R n [5], and a ratio function P n [38], giving (2. 2) The BDS ansatz encodes the infrared-divergent part of the amplitude (in dimensional regularization) as well as its non-DCI part. This leaves the remainder function and ratio function, which are finite and respect dual conformal symmetry. The ratio function carries all dependence on the Grassmann variables and can be expanded in Grassmann degree to isolate contributions with different helicity structure: The kinematic dependence of the remainder and ratio function is most conveniently written in terms of dual variables x i and θ i , which are defined in terms of the external momenta k i and the Grassmann variables η A via Here the additional index i on the Grassmann variables associates them with the i th external particle, and λ,λ are two-component spinors. For more background on these dual (super-)coordinates, see for instance ref. [98]. Focusing now on the case of six particles, there are only three helicity configurations: MHV, NMHV, and MHV. The MHV amplitude is parity conjugate to MHV, and thus the entire amplitude is encoded in the MHV and NMHV components. After normalizing them by the BDS ansatz, these amplitudes only depend on the kinematics through dual-conformal invariant cross ratios. Only three algebraically independent cross ratios can be formed for six particles, and traditionally they have been chosen as are squared differences of dual coordinates. These cross ratios can also be expressed in terms of (planar) two-and three-particle Mandelstam invariants using the translation s i,i+1,...,i+n−1 = (k i + k i+1 + · · · + k i+n−1 ) 2 = x 2 i,i+n . The MHV amplitude corresponds to the leading term in the expansion (2.3), and as such it depends only on the remainder function R 6 (u, v, w) defined in eq. (2.2). This function is expected to be a pure (generalized) polylogarithmic function (to be defined more precisely below) of the cross ratios (2.5) to all loop orders, meaning that the kinematic dependence only appears in polylogarithms and not in any rational prefactors multiplying these functions. The NMHV contribution to the ratio function P NMHV is not pure, but it can be written as a sum of pure polylogarithmic functions multiplied by R-invariants.
To define the R-invariants, we first recall the definition of momentum supertwistors [99,100], The momentum twistors Z a , considered as vectors in CP 3 , can be contracted into SL(4) invariants using the Levi-Civita tensor, where in particular x 2 ij ∝ i − 1ij − 1j , and the additional factors (which would make this proportionality exact) cancel out in the fully dual-conformal invariant ratios in eq. (2.5). Similarly, five-brackets of the Z a are dual superconformal invariants known as R-invariants. They are defined to be (f ) ≡ [abcde] = δ 4 χ a bcde + cyclic abcd bcde cdea deab eabc , (2.8) where we denote the five-bracket of legs {a, b, c, d, e} by the remaining leg f . Using these quantities, we can parametrize P NMHV in terms of a parity-even function V and a parity-odd functionṼ : The remainder function, V , andṼ have perturbative expansions in the large-N coupling g 2 defined by where λ is the usual 't Hooft coupling. The coefficients of the perturbative expansion are linear combinations of generalized polylogarithms. These functions are defined as iterated integrals over logarithmic kernels, commonly denoted by with the recursion starting at G(z) = 1. The number n of nested integrations in the above definition is referred to as the weight of the generalized polylogarithm. At L loops in the perturbative expansion in g 2 , the remainder function, V , andṼ and are all pure functions of uniform transcendental weight 2L. In practice this means that the total differential of these functions, and the functions F we will be interested in later, can be written in the form dF = s∈S hex F s d ln s , (2.12) where S hex is the set of nine hexagon symbol letters that will be introduced below. If F has weight n, then each F s is a pure function of weight n − 1.
Because of the motivic structure of generalized polylogarithms [17,51,54,56,68,101], the total differential (2.12) corresponds to the component of the coaction with weight-one functions in its back entry, The map ∆ •,1 can be applied iteratively to the functions in the left factor of the coproduct, breaking down the weight n polylogarithms entering F into n-fold tensor products of logarithms, objects commonly referred to as symbols [17]. The letters s appearing in the arguments of the logarithms belong to the hexagon symbol alphabet: (2.14) Here we have introduced three parity-odd letters, 15) defined in terms of the quantities The cross ratios (2.5) are parity even.) In this paper we will often use the following equivalent symbol alphabet, which is related to the original alphabet S hex through the relations This alphabet makes the Steinmann relations more transparent by isolating the three independent three-particle Mandelstam invariants s i,i+1,i+2 in different letters, a, b and c. In this respect it resembles more closely the usual choice of alphabet for the seven-particle amplitude [30]. At higher loops, expressing these functions in the G-function notation (2.11) becomes overly cumbersome, and it proves more effective to encode the amplitude in terms of its ∆ •,1 coproduct. This amounts to keeping track of the (iterated) derivatives of the amplitude, as well as the integration constants required to reconstruct the amplitude from these derivatives. This way of treating these functions will be reviewed in more depth in a companion paper [66], and has been discussed elsewhere in the literature (see for example ref. [24]). In the next subsection, we describe how we normalize the amplitudes so that they lie in H hex .

Cosmic normalization
In ref. [66] we describe the construction of a minimal space of (extended) Steinman-satisfying hexagon functions. In addition to the extended Steinmann relations, this space is constructed to obey a coaction principle, which is to say it is invariant under the cosmic Galois group [67][68][69][70]. This condition implies that if one studies the coaction of higher-loop amplitudes, the first entry belongs to a stable space whose dimension at fixed weight saturates as the loop order increases. Empirically, this saturation happens around weight L, namely at around half the weight of the corresponding L-loop amplitude. While some of these restrictions have long been understood at symbol level, at the level of functions they result in further restrictions on the transcendental constants that are allowed to appear as free elements: only ζ 4 , ζ 6 , ζ 8 , etc. are needed.
It is highly nontrivial that the six-point amplitudes lie within this space through seven loops. In fact, they only do so once they are properly normalized. First of all, in order to preserve the Steinmann constraints, we must normalize the amplitude by the BDS-like ansatz [31,33,102] (up to a kinematically constant factor). The BDS-like ansatz differs from the BDS ansatz by a factor of where the cusp anomalous dimension is [103] and is the finite, dual conformally invariant part of the one-loop MHV amplitude. However, beyond two loops this normalization is not sufficient to place the amplitudes into H hex , and needs to be adjusted by a kinematical constant. (The Steinmann relations uniquely fix the normalization, up to this constant.) One way to see the need for a new normalization is to inspect the three-loop values of the MHV amplitude E old and NMHV amplitude E old defined in ref. [33] at the point (u, v, w) = (1, 1, 1): Both amplitudes contain (ζ 3 ) 2 with the same nonzero coefficient. The presence of (ζ 3 ) 2 would violate the coaction principle, because of a specific term in its coaction, namely ∆ 3,3 [(ζ 3 ) 2 ] = 2ζ 3 ⊗ ζ 3 . This term contains a ζ 3 in its first entry, but no ζ 3 is allowed there, because the weight-three functions in H hex all vanish at (1, 1, 1) [66]. Therefore, starting at three loops, we modify the BDS-like ansatz A BDS−like 6 by a function of the coupling ρ(g 2 ), in order to fit the amplitudes into a space that is invariant under the cosmic Galois group. In other words, we define the "cosmically normalized" MHV amplitude E in terms of the full, infrared-divergent amplitude A 6 , the BDS-like ansatz, and ρ, as (2.23) The BDS-like ansatz is [31,102] A BDS−like where C (L) is a constant at each loop order, as is f (L) ( ), (2.25) (We suppress the two additional terms in the expansion of f here for simplicity.) Finally, the quantityM 6 ( ) differs from the full one-loop MHV amplitude by an amount proportional to E (1) . It is given bŷ where γ E is the Euler-Mascheroni constant. It is important thatM 6 ( ) only depends on the two-particle invariants s i,i+1 to O( 0 ), so that factoring it out of the amplitude does not affect the Steinmann relations for three-particle invariants [33].
The relation between E and the remainder function R 6 is The coefficient functions that specify the cosmically normalized NMHV amplitude, which we call E andẼ, are defined similarly; in terms of the scheme-independent coefficient functions V andṼ of the ratio function (2.9) they are defined by To pass from the old normalization for E old , E old ,Ẽ old used in ref. [33] to the one used here, simply divide the old functions by ρ(g 2 ).
We compute ρ by requiring that the amplitudes fit into the minimal space that obeys the coaction principle. For example, at three loops we need to have ρ(g 2 ) = 1+8(ζ 3 ) 2 g 6 +O(g 8 ) in order to cancel the two appearances of (ζ 3 ) 2 in eq. (2.22). This criterion does not necessarily determine ρ uniquely. For instance, ρ could be adjusted by terms involving even zeta values ζ 2L at L loops since these have free parameters associated with them (except for ζ 2 at one loop). We choose to fix this ambiguity by not including any ζ 2L term at L loops. Through seven loops, we find that the minimal solution for ρ does not require any genuine MZVs, nor does it contain any factor of ζ 2 : (2.29) This form for ρ is uniquely fixed, given the following assumptions: 1. ρ does not contain MZVs of depth two or higher, e.g. no ζ 5,3 , 2. the coaction principle is satisfied through weight 14 at the point (1, 1, 1), 3. it is also satisfied at an analytic continuation of this point to 3 → 3 self-crossing kinematics, 4. a subspace of the hexagon functions that we can define to all weights, which saturates the space of MZVs in H hex through weight 10, also does so at weight 11.
The last constraint is only needed to uniquely fix ρ (7) . It imposes one additional constraint on the MZVs appearing at weight 11, which in turn implies one fewer allowed MZV at weight 14. Otherwise we would have an ambiguity in ρ at seven loops, because we do not yet know the NMHV amplitude at this order. Once we have fixed ρ at a given loop order, many constraints ensue as we go to higher loop orders. For example, none of the {6, 1, 1} coproducts of the four-loop amplitudes can contain a (ζ 3 ) 2 when evaluated at (1, 1, 1). It is interesting that the same value of ρ works for both the MHV and NMHV amplitudes. Perhaps this fact indicates that ρ is determining a particular infrared regularization scheme.
In fact, ρ seems to be related to the inverse of the cusp anomalous dimension at low loop orders. 1 More precisely, the quantity ρ − 4 g 2 /Γ cusp contains only even zeta values until five loops: Recall that the g 4 , g 6 and g 8 terms could be removed from this relation by redefining ρ, using the fact that the corresponding zeta values are independent, constant elements of H hex . Given how the cusp anomalous dimension is tied to infrared divergences, eq. (2.30) might provide a clue for how ρ defines a "cosmically-preferred" infrared regularization scheme. Through five loops, the values of the MHV amplitude at (1, 1, 1) are, in the new normalization, For the NMHV amplitude (parity even part) they are We will give the six-and seven-loop values in section 4.1.

Bootstrapping the six-particle amplitude
In this section we describe the constraints that we impose in order to uniquely determine the MHV amplitude through seven loops, and the NMHV amplitude through six loops.

Discrete symmetries,Q supersymmetry, and collinear limit
Our starting point is the space of functions described in ref. [66]: generalized polylogarithms with symbol letters drawn from the alphabet S hex of eq. (2.14), that have branch cuts ending where the cross ratios u, v, and w become zero or infinity, and that satisfy the extended Steinmann relations and obey a cosmic Galois coaction principle. The general linear combination of all such functions with transcendental weight 2L forms an initial ansatz for each of the (cosmically normalized) functions E (L) , E (L) , andẼ (L) . These functions inherit a set of discrete symmetries from the dihedral symmetry of the full superamplitude: E is fully symmetric under all permutations of u, v, w, while We also impose the condition [29] which removes an unphysical degree of freedom in the functionẼ, as any cyclically symmetric piece ofẼ drops out of the full amplitude due to the following linear relation between Rinvariants, At this early stage it is convenient to impose the "final entry condition" that follows from theQ equation [21,104,105]. The dual superconformal generatorQ is a first-order differential operator acting on the n-point L-loop N k MHV amplitude, which relates it to an integral over the (n + 1)-point (L − 1)-loop N k+1 MHV amplitude. By choosing appropriate differentials, the latter "source term" can be made to vanish, leading to a set of homogeneous constraints. In the MHV case, they take the form, These relations imply that the final entries of E can only be drawn from a subset of the normal symbol letters, which we may write in the alphabet (2.17) as The NMHV final entry conditions are slightly more involved [31]. The allowed final symbol entries here depend on which R-invariant the polylogarithm multiplies. The allowed combinations are along with their cyclic images. Note that the appearance of products like d ln(m w y u ) means that these constraints link E andẼ; after imposing these final entry conditions, our ansätze for these functions can no longer be considered independent. We apply these constraints first because they are easy to impose, and because they constrain a large number of coefficients. For example, at six loops, of the 3692 undetermined coefficients present in our initial weight 12 ansatz, only 236 and 102 parameters remain in the MHV and NMHV amplitudes, respectively, after imposing the above symmetries and final entry conditions. We next turn to the collinear limit, where the six-point amplitude reduces to a five-point amplitude times a universal splitting function. Since the six-point BDS ansatz captures the correct collinear limits of the six-point amplitude, both the remainder function and the ratio function must vanish at loop level in this limit. E and E do not vanish, but their behavior can be easily found by taking the collinear limit of eq. (2.27). In the u, v, and w variables, the collinear limit corresponds to v → 0 and u → 1 − w, where we have that Similar constraints apply in the cyclically related collinear limits, but are automatically enforced by the discrete symmetry conditions already imposed. Imposing this behavior in the strict collinear limit (which is to say, at leading power) fixes many of the remaining parameters in our ansatz. As we will see in the next section, it leaves us with fewer than 20 free parameters to constrain.

Multi-Regge kinematics
Next we constrain these remaining parameters, using input from additional kinematic limits. The multi-Regge limit, in which the collision of two highly energetic particles produces four particles that are strongly ordered in rapidity, is a prime source of such input. In terms of Mandelstam invariants, with particles 3 and 6 incoming and particles 1,2,4,5 outgoing, the 2 → 4 multi-Regge limit reads In this section, we will use (u 1 , u 2 , u 3 ) instead of (u, v, w) to avoid confusion with the multi-Regge variable w. In the cross ratios the multi-Regge limit becomes, with the ratios held fixed in terms of the complex-valued multi-Regge variable w. The behavior of the remaining letters, may be inferred from the following momentum twistor parametrization, where the columns correspond to (Z 1 , . . . , Z 6 ), the rows correspond to their components, and the limit amounts to √ u 2 u 3 τ → 0. In the Euclidean region, loop corrections to the BDS-normalized amplitude vanish in the multi-Regge limit [40], due to its conformal equivalence to a soft limit [24,106]. Nontrivial behavior in the limit is obtained by analytically continuing into physical 2 → 4 Minkowski kinematics. For particles 3 and 6 in the initial state, the Mandelstam variables s 12 and s 45 in the numerator of u 1 should be continued into the time-like region, which amounts to continuing u 1 → e −2πi u 1 (3.13) before taking the limit. As we will review shortly, it is possible to obtain direct predictions for the behavior of the amplitude in multi-Regge kinematics when the external particles are gluons. For the bosonic MHV amplitude this is true by definition, whereas for the NMHV superamplitude, we need to specialize to its gluonic components. Without loss of generality we may choose the (χ 4 ) 4 component of the ratio function (that is, P NMHV ) [31], which describes the NMHV helicity configuration 2 3 Taking the limits of the R-invariants that the functions E andẼ multiply, we may express the corresponding BDS-normalized gluon amplitude in multi-Regge kinematics as The first factor in this relation takes into account the change from (cosmic) BDS-like normalization back to BDS normalization. Moreover, all transcendental functions should be understood as having been first analytically continued as (3.13), and then evaluated in the limit (3.9). For instance, in this limit, where we also used the shorthand ww * ≡ |w| 2 etc. The amplitude is invariant with respect to the discrete y 2 ↔ 1/y 3 transformation, known as target-projectile symmetry. The latter is a combination of a dihedral flip (y 2 ↔ y 3 or Z i → Z 3−i ) and a parity (y i → 1/y i or Z i → Z i+3 ) transformation. In the multi-Regge limit it becomes equivalent to the inversion w → 1/w, w * → 1/w * . Thus if our ansatz for E,Ẽ already respects target-projectile symmetry in general kinematics, we need only compute the first line of eq. (3.15), and the second line follows from the first by replacing w → 1/w.
The behavior of the amplitude in multi-Regge kinematics may be studied directly within the BFKL approach [29,40,[71][72][73], yielding dispersion relation-type integrals for the amplitude. These integrals have the form, where and Γ cusp is the cusp anomalous dimension, given in eq. (2.20). The remaining building blocks of the integrand are known as the BFKL eigenvalue ω, the regularized (MHV) impact factor Φ Reg and the NMHV form factorH. 3 They were first obtained to the first few orders in perturbation theory by analyzing the effective particles whose exchange dominates in the limit, also by studying the limits of the amplitudes [26], and more recently to all loops with 3 The NMHV form factor is essentially the ratio between the regularized MHV and NMHV impact factors.
the help of integrability [74]. We refer to the latter reference for their expansion to arbitrary loop order, where in the notation therē Finally, the P on the right-hand side of eq. (3.17) denotes the Cauchy principal value of the integral, which has a pole at n = ν = 0. Equivalently, we take half the value of the corresponding residue, when evaluating the integral using Cauchy's residue theorem.
Given that all building blocks in the integrand (3.17) start at O(1), except for ω(ν, n) which starts at O(g 2 ), it is easy to see that at L loops R (L) 6 and P NMHV can be written in the limit as a polynomial in ln(1 − u 1 ), or ln τ , of degree L − 1 (plus power-suppressed terms). The coefficient of ln L−1 τ is referred to as the leading-log (LL) contribution, and similarly the coefficient of ln L−1−k τ corresponds to the (next-to) k -leading logarithmic (N k LL) contribution. We provide the result of evaluating the Fourier-Mellin integral in eq. (3.17), separated by logarithmic order, in the ancillary file hexMRKL1-7.m.
Aside from the rational prefactors 1/(1 + w * ) and w * /(1 + w * ) in the NMHV case, the coefficients of the large logarithms are single-valued harmonic polylogarithms (SVHPLs) in w, w * [107,108]. From eqs. (3.10) and (3.12) it can be seen that these functions have the symbol alphabet {w, 1+w, w * , 1+w * }, and that the first entry can be only ww * or (1+w)(1+ w * ). SVHPLs have the important property that they can be uniquely reconstructed from their holomorphic part, defined to be their value at w * = 0, also removing any divergent ln w * . This observation significantly simplifies the computation of the dispersion integrals (3.17), since the holomorphic part comes only from the residues satisfying iν = −n/2 [109,110] (see also ref. [106]).
We can match this holomorphic part directly to the holomorphic part of the multi-Regge limit of our ansatz for the amplitude. In order to obtain this limit, we iteratively compute the multi-Regge limit of each function in our basis through weight 11. This can be done using the coproduct entries of these functions, which encode the derivative of each function with respect to w. (The multi-Regge limit of all the functions appearing in these coproduct entries is known from the previous iteration in the computation.) Each function's derivative can thus be matched to a basis of SVHPLs of the appropriate weight. (There are 3259 SVHPLs at weight 11, including MZVs multiplied by lower weight SVHPLs.) We fix the constants of integration by first computing all functions on the 2 → 4 self-crossing line where (u, v, w) → (1 − δ, v, v), as described in section 3 of ref. [91] and section 5 of this paper, and then sending v → 0; this limit intersects the w → −1 limit of the multi-Regge limit. In the coproduct representation of our ansatz, matching to holomorphic data means that we only need to consider the final symbol letters w and 1 + w. At seven loops, matching at the level of weight 11 functions means that we compare to the {11, 1, 1, 1} coproduct component. This leaves three integration constants unfixed. However, the constants can be determined as described above, using the self-crossing limit of just the coproduct entries of the amplitude itself, rather than for all weight 12, 13, and 14 functions in H hex . (An analogous procedure was used to evaluate the strict collinear and near-collinear limits, where the boundary conditions are also fixed using the same self-crossing line, but this time on the Euclidean sheet where it intersects the collinear limit at its soft endpoint.) Since ω(ν, n) is zero at tree level, any term containing a square or higher power of ln τ is determined entirely by lower-loop data [33]. That is, at any given loop order such terms may be predicted using only lower-loop amplitudes and the structure of eq. (3.17), without the need for prior knowledge of the precise form of the integrand building blocks from, say, integrability. If all parameters in our ansatz can be fixed by such terms, then our bootstrap is in some sense "pure", in that it does not rely on external data.
Quite interestingly, while the multi-Regge limits at the ln 2 τ level fixed all remaining parameters through five loops [33], they no longer do so at six or seven loops. The NMHV amplitude can be fixed at six loops using just this data, but in the MHV amplitude a single parameter evades determination in the multi-Regge limit at both six and seven loops, even when the independent predictions of ref. [74] are used for the most subleading large logarithms ln τ and ln 0 τ . We will analyze the function accompanying this final six-loop parameter in more detail in subsection 3.4; however, its appearance is not a problem. As we now describe, another limit can fix this coefficient. This limit can, in principle, supply an infinite amount of boundary data for the amplitude.

Near-collinear operator product expansion
Independent predictions for the behavior of the amplitude in an expansion around the collinear limit may be obtained within the framework of the Wilson loop (or Pentagon) Operator Product Expansion [75][76][77][78][79][80][81][82][83][84][85][86][87]. The limit in question is most conveniently described in terms of variables T = e −τ , S = e σ , and F = e iφ that appear in the following parametrization of the momentum twistors Z i [77], where the limit corresponds to τ → ∞ or T → 0. We may express any conformally invariant cross ratio in terms of T , S, and F with the help of the above parametrization, and the definition of the cross ratio in terms of four-brackets (2.7). For example, the cross ratios u, v, and w evaluate to As its name suggests, the natural object in the Wilson loop OPE approach is not the amplitude per se, but the closely related, appropriately normalized (framed) Wilson loop. Focusing momentarily on the MHV case, the relation between E and the corresponding Wilson loop W is where E (1) is defined in eq. (2.21), and X is given by [77] The normalization factor ρ appears in eq. (3.24) to compensate for the fact the Wilson loop is related to the original BDS-like normalized amplitude, not its cosmically normalized cousin. Notice that unlike E, the function X and thus also W are not dihedrally symmetric. The asymmetry is a consequence of the framing of the hexagonal Wilson loop, which requires dividing it by two pentagonal Wilson loops and multiplying back by a quadrilateral one; the orientation of these auxiliary polygons breaks the dihedral symmetry.
The collinear limit corresponds to T → 0, and the Wilson loop OPE is nothing but an expansion of W in powers T E i for different exponents E i , which govern how fast the limit is approached. Both the exponent and the coefficient of these powers can be computed exactly in the coupling, as a consequence of the physical interpretation of each term as an excitation of an integrable flux tube, formed by the Wilson lines of the quadrilateral in the aforementioned framing of the Wilson loop.
At loop order L in the weak coupling regime, which is the focus of this paper, the OPE framework predicts that where f (L) m,n (T, S) is a polynomial (of degree L − 1 for MHV) in ln T , whose coefficients are sums of rational times transcendental functions of S. In the above, m and n correspond respectively to the total particle number and helicity of the flux tube excitation, both of which are good quantum numbers.
Here, we will mostly consider the excitations with m = |n|, and n = ±1 or ±2, whose contributions to the near-collinear limit have the form, |n|,n (T, S) . (3.27) Given that the field content of the theory consists of scalars, fermions and gluons with helicity 0, ±1/2 and ±1 respectively, these contributions only come from gluonic excitations of the same helicity. It can be shown that gluons of the same helicity may also form bound states [78,111]. The flux tube description of these contributions is where [80] W denotes the contribution of a single-particle excitation which is a bound state of |a| gluons with helicity a and rapidity u, and is a superposition of two such excitations, consisting of gluons with helicities a and b such that ab > 0, and rapidities u and v, respectively. The physical quantities p a and γ a = E a − |a| are the momentum and quantum correction to the energy of the flux tube excitation, while µ a and P a|b are the measure and pentagon transition, respectively. The finite-coupling expressions for all these quantities have been provided in ref. [80], together with a method for their systematic weak-coupling expansion. In order to avoid clutter, here we illustrate only the leading terms for ab > 0: We provide expressions for these quantities to eight loops in the ancillary file WLOPEblocks.m. Because the MHV Wilson loop is a bosonic scalar object, it is invariant under parity, which flips the sign of the helicity. This symmetry implies that as can also be seen explicitly in the leading-order expressions above. From the same expressions, it is evident that W [a] will start at O(g 2 ). At L loops, it is a polynomial in ln T of degree L − 1, whereas W [a,b] starts at O(g 8 ), and has degree L − 4. Remarkably, this implies that the single-gluon state has a distinct signature at order T 2 , since at L loops all ln k T contributions to f 2,±2 with L − 1 ≤ k ≤ L − 3 come purely from it. A general algorithm for the weak-coupling evaluation of W [±1] was developed in ref. [112], and later extended to W [±2] [113]. In both cases explicit expressions were obtained up to six loops. 4 In the latter paper, W [±1,±1] was computed at four loops, by observing that although it is represented as a two-fold integral (3.30), it can always be reduced into a sum of products of one-fold integrals, to which the method of ref. [112] is applicable. Because this property is independent of loop order, it was also used [109] to determine W (5) [±1,±1] ; here we have pushed these three OPE contributions to seven loops. We include them in the files WL0-6.m and WL7.m, through six loops and at seven loops respectively.
Thus we have the complete OPE prediction through seven loops for the T F ±1 and T 2 F ±2 terms in the T → 0 limit of the MHV amplitude E (L) , after using the conversion (3.24). At six loops, we compare this prediction to the limit of our ansatz, with the single parameter that survives the multi-Regge limit. Quite interestingly, at six loops we find that the single such parameter is fixed from the T 2 F 2 ln 4 T contribution, namely from the single-gluon bound state alone. (Recall that at six loops, the two-gluon superposition only contributes to terms with at most two powers of ln T .) We have checked that our answer for the amplitude also agrees with all the remaining T 2 F ±2 terms, including those sensitive to the two-gluon superposition. Our answer for the limit of the amplitude also contains the T 2 F 0 terms, which we also include in the ancillary file WL0-6.m. These terms have not been checked against the OPE, although of course we expect them to agree.
At seven loops, we proceed in a slightly different fashion. The combined constraints, up to and including the multi-Regge limit, again leave one parameter free in the ansatz. Because the near-collinear limit is computationally demanding at seven loops, we instead examine the origin, where (u, v, w) → (0, 0, 0). We observe that through six loops the remainder function has only a very mild singularity, quadratic in logarithms of u, v, w. (See section 4.2.) On the other hand, the function multiplying the one free parameter is far more singular at the origin, behaving like ln 6 u ln 4 v ln 4 w + cyclic. Therefore we fix its coefficient by requiring only quadratic behavior in the remainder function at seven loops. We have verified that this choice for the last parameter matches the OPE predictions in the near-collinear limit for all T 1 F ±1 and T 2 F ±2 terms, i.e. to the same accuracy checked at six loops. 5 We have also investigated the near-collinear limit of the final seven-loop ambiguity, and find that its limit contains terms of order T 2 F ±2 ln 6 T , meaning that its coefficient can also be fixed using just the single-gluon bound state. Again, we extract the T 2 F 0 terms in the Wilson loop OPE from the amplitude in this limit, and we include them in the ancillary file WL7.m.
Moving on to the NMHV case, we will focus on taking the near-collinear limit of the (1111) component of the superamplitude, which may be expressed in terms of the E,Ẽ coefficient 4 In particular, the relations between our notations, and those used in the aforementioned papers, are functions as   6 Note that ref. [77] uses ηi, rather than χi, to refer to dual supercoordinates.
are NMHV form factors, responsible for creating OPE excitations that are charged under the R-symmetry of the theory. They are expressible in terms of so-called Zhukowski variables, for which we have also indicated the choice of branch at weak coupling. In contrast to the MHV case, the presence of these form factors implies that OPE contributions with opposite helicity will no longer be equal, e.g. W . Nevertheless, we may also evaluate them using the existing summation algorithms mentioned earlier in this section, after taking into account two minor differences for the case a > 0. First, the presence of inverse powers of the coupling in eq. (3.39) implies that we have to expand the MHV integrand multiplying the form factor to higher order in g 2 than the loop order we wish to compute. And second, particularly for W (1111) [2] , it proves simpler to compute the integral by first defining a reduced integrand f (u) to be the original integrand divided by (u 2 + a 2 /4), then evaluating the integral of the reduced integrand by residues. From that integral we obtain the result for the original integrand by differentiating, 7 In this manner, we have also determined the T F ±1 and T 2 F ±2 terms in the collinear OPE expansion of the (1111) NMHV super-Wilson loop component to six loops. All terms computed agreed with the near-collinear limit of the NMHV amplitude, which is an independent cross check, since all parameters were fixed by the multi-Regge limit in the NMHV case. We include the near-collinear limit of this component of the super-Wilson-loop in the ancillary file W1111L0-6.m, including T 2 F 0 terms which have not yet been checked against the OPE. Table 1 summarizes the number of parameters left after imposing the various constraints we have discussed, through six loops. The left-hand entry in parentheses is the number for the MHV amplitude, while the right-hand one is NMHV. "MRK" refers to multi-Regge kinematics, with the number of N's indicating how many logarithms below the leading logarithms (LL). There are a couple of question marks in the table at six loops, or weight 12, having to do with the number of parity-odd functions that are not allowed by branch cut and other conditions, even though their symbol is allowed by symbol-level constraints [66]. The issue arises because we do not yet have a complete weight 12 basis at function level, but it is unlikely to affect the final number by more than one. There are also some asterisks related to how the constant ρ is fixed. If ρ were known ahead of time at a given loop order, the number of parameters would be exactly the number shown. But in practice it is not, and so one should add the number of asterisks indicated. For example, there are two extra parameters associated with ρ for the six loop MHV amplitude after the strict collinear limits are imposed, in addition to 6 other parameters that are independent of the value of ρ.   Table 1. Remaining parameters in the ansätze for the (MHV, NMHV) amplitude after each constraint is applied, at each loop order. The superscript " * " (" * n") denotes an additional ambiguity (n ambiguities) which arises only due to lack of knowledge of the cosmic normalization constant ρ at the given stage. The "?" indicates an ambiguity about the number of weight 12 odd functions that are "dropouts"; they are allowed at symbol level but not function level. The seven-loop MHV amplitude was constrained in a somewhat different order. As the parameter counts are not directly comparable it is omitted from the table.
The table clearly shows the difference between six loop MHV and all other cases shown, in that one parameter survives all the way through the multi-Regge limit and the one flux tube (T 1 ) OPE constraints. The same is true at seven loops. We omit the seven-loop MHV numbers because we constrained the result somewhat differently, using symbol-level constraints first and then reapplying the constraints at function level, so the numbers are not directly comparable. However, the number of surviving symbol-level parameters after the strict collinear limit is 17. This number is consistent with the general pattern of the number of parameters rising by roughly a factor of 3 per loop.

The fate of inter-loop relations
Once we arrive at the final expressions for the six-particle MHV and NMHV amplitudes through seven and six loops respectively, following the steps we described in the previous subsections, we move on to examine their properties. Most of the quantitative and qualitative analysis in various points and lines of the space of kinematics will be done in the next section, but let us conclude this section by noting another new feature that first appears at six loops.
Based on empirical observations up to L = 4, a relationship between the L-loop MHV amplitude and the (L − 1)-loop NMHV amplitude was conjectured [29]. This relation was then confirmed to hold also at L = 5 [33]. In our notation, which essentially coincides with that of the latter reference, the relation takes the form, (3.42) where F x,y refers to the corresponding component of the ∆ n−2,1,1 "double" coproduct of the function F , and E = E(u, v, w).
Here, we observe that this relation in fact breaks down at six and seven loops. Very interestingly, this breakdown is closely related to the phenomenon described at the end of section 3.1, namely the appearance of functions in H hex that vanish in the near-collinear limits through O(T ) and in the multi-Regge limits. These functions are precisely the ones multiplying the parameters of our ansatz that remain free after applying the corresponding constraints. We saw that there exists one such function at weight 12 that contributes to the MHV amplitude, which we can call Z(u, v, w) for concreteness. A similar function appears at weight 14, which we callZ(u, v, w).
Both Z andZ are totally symmetric under dihedral S 3 transformations. Moreover, their parity-even {2L − 1, 1} coproducts vanish identically, (3.43) so that they are entirely specified by their y u coproduct. The T derivative of any function satisfying eq. (3.43) can be expressed in terms of the coproducts Z yu , Z yv and Z yw , with coefficients that are O(T 0 ) as T → 0. Because Z is parity even, the Z y i are parity-odd, and so they vanish like T 1 as T → 0 (times powers of ln T ). Thus the T derivative of Z vanishes like T 1 ln k T , and so Z itself must vanish like T 2 ln k T (provided that Z is not a constant). In other words, eq. (3.43) alone is enough to ensure that a parity even function Z is undetectable in the near-collinear limit at the level of one flux tube excitation. (It is not as clear to us why such a function has to vanish in the multi-Regge limit.) Why didn't such functions turn up at lower loop orders? The answer is that within our function space there are no parity-even solutions to eq. (3.43) until weight 12! Returning to the connection between the six loop function Z and the relation (3.42), we find that we can satisfy this relation if we shift E (6) → E (6) + αZ(u, v, w) for some rational number α. In other words, the only effect of the MHV-NMHV relation at six loops is to set the coefficient of Z(u, v, w) in E (6) to an incorrect value. This is nontrivial, since eq. (3.42) could alternatively yield no solution at all. We suspect that its validity through five loops is linked to the non-existence of analogs of Z(u, v, w) at lower weight. We should also remark that at seven loops it is not possible to solve eq. (3.42) solely by shifting E (7) → E (7) + βZ(u, v, w) for any β.
Could there be other relations between MHV and NMHV amplitudes and their double coproducts at one loop higher? By surveying our data up to six loops, we find three such relations: (3.44) and its cyclic images. It remains to be seen whether they continue to hold at higher loops, but paralleling our analysis of the breakdown of eq. (3.42), they may be related to the absence or presence of functions that vanish even faster than Z in the near-collinear limit.

Numerics and number theory in kinematic limits
In the ancillary file SixGluonAmpsAndCops.m accompanying this paper, the six-and sevenloop functions E, E, andẼ are expressed in terms of iterated ∆ •,1 coproduct entries, supplemented by boundary data at the point (1,1,1). This is formally equivalent to providing these functions in the generalized polylogarithm notation (2.11) but is far more compact. While at five (and lower) loops, both types of expressions have been published [33,114], the polylogarithmic expressions become too unwieldy by weight 12. In fact, it already proves prohibitively difficult to use these expressions to generate numerics at five loops. However, these functions simplify drastically on a number of codimension two and three subspaces. This allows us to more easily probe the analytic and numerical properties of the amplitudes. In particular, it was observed by several of the authors [24] that the remainder function appears to behave extremely similarly across loop orders for moderate values of the cross ratios (0 u 1). This was interpreted [26,31] as a signature of the finite radius of convergence of the perturbative series in planar N = 4 SYM, a property that is expected from the absence of instantons and renormalons. It can be explicitly demonstrated for some quantities that are known at finite coupling, such as the cusp anomalous dimension [103]. The small-coupling expansion of the cusp anomalous dimension has a radius of convergence of 1/16. Stated differently, at high orders the ratio between the coefficients of successive orders in the coupling approaches −16. 8 From prior observations at lower loops, we expect that for values of u, v, w of order unity, the ratios of successive loop orders for the functions E,Ẽ, and E should approach a similar radius of convergence, and should do so quite rapidly in the loop expansion. Later in this section, we will confirm that this behavior continues to hold at six and seven loops by plotting the ratios of successive loop orders. First we discuss the functions' behavior at the point (1, 1, 1) and at the origin.

The origin
In the limit where u, v, w all approach zero, all hexagon functions become polynomials in the logarithms ln u, ln v, and ln w. We can determine the (transcendental) coefficients of this polynomial for the functions E, E, andẼ by first considering them at the point u = 1 and v, w → 0. This point is in the collinear limit, so we know the values of these three functions here. To translate this information to the origin, we integrate the amplitude along the (u, 0, 0) line, where v, w → 0 but u is generic. Along this line, all hexagon functions become harmonic polylogarithms in u with coefficients that are polynomials in ln v, ln w. A generic weight-2L hexagon function gives rise to a polynomial of total degree 2L in the logarithms. However, from the leading OPE behavior on the double-scaling surface v → 0 one can show that for the L-loop amplitude, this polynomial should have a degree of at most L in ln v, and similarly for the degree in ln u and in ln w [18,19]. 10 That is indeed the maximum logarithmic behavior that we find for the even and odd contributions to the ratio function, V andṼ . Remarkably, for the remainder function through seven loops we find much less singular behavior, at most quadratic in the total power of the logarithms. Also, linear terms in the logarithms do not appear. Because R 6 (u, v, w) is totally 10 The original analyses focused on amplitudes with a one-loop leading OPE contribution. More generally, it can be shown that if the leading OPE contribution an amplitude receives is at k loops, the highest degree of logarithmic divergence is L − k [116], with k = 0 being of course maximal. symmetric, and there are only two totally symmetric quadratic polynomials, we find that (4.10) as (u, v, w) → (0, 0, 0), where P 1 (u, v, w) = ln 2 u + ln 2 v + ln 2 w + ln u ln v + ln v ln w + ln w ln u, (4.11) P 2 (u, v, w) = ln u ln v + ln v ln w + ln w ln u. 2 ) all have vanishing coefficients. Although we do not have a proof that the remainder function is at most quadratic in ln u, ln v, ln w, there is an expectation in the literature that this is indeed the case [102]. Assuming that it holds, it's worth noting that this constraint can be used to eliminate the MHV ambiguity arising from the function Z(u, v, w), which is not fixed by the multi-Regge limit, nor by the OPE at the level of one flux-tube excitation. (These two constraints are apparently closely related by analyticity.) In the limit (u, v, w) → (0, 0, 0), Z is highly singular, with up to 12 total powers of logarithms (but at most 4 powers in any individual logarithm): (4.31) So in principle we could have used the quadratic logarithmic behavior of R 6 (u, v, w) in the (u, v, w) → (0, 0, 0) limit, right after the multi-Regge constraint, to remove the final parameter multiplying Z. In practice, this is how we removed the parameter multiplying the analogous functionZ at seven loops, using the fact that Z(u, v, w) ∼ ln 4 u ln 4 v ln 4 w(ln 2 u + ln 2 v + ln 2 w) . (4.32) After fixing that last parameter, we then verified the OPE limits.

The line (u, u, 1)
Now let us consider the line on which (u, v, w) = (u, u, 1). Here, the hexagon symbol alphabet S hex can be written in terms of just two letters, {u, 1 − u}. This means that functions on this line can be written in terms of harmonic polylogarithms [92] in u, with weight vector Thus, only E(u, u, 1), E(u, 1, u), and E(u, u, 1) represent independent functions. In the ancillary file SixGluonHPLLines.m, we provide these functions through six loops (seven loops for E(u, u, 1)).
We plot E(u, u, 1), E(u, 1, u), and E(u, u, 1) in figures 1, 2, and 3 respectively, through six loops (seven loops for E(u, u, 1)), as ratios of successive loop orders. For u between 10 −2 and 10 2 , remarkably, the ratios flatten out more and more with each additional loop, and they appear to be steadily approaching the cusp value of −16. Near u = 0.1 for E(u, u, 1), and for u between 0.1 and 1 for E(u, u, 1) there is a dip/spike feature, which is simply because each function crosses zero at a slightly different value of u. For u → 0 and u → ∞, the ratios no longer display the expected radius of convergence, either diverging logarithmically at different rates or, for the examples of E(u, u, 1) and E(u, u, 1) below, approaching constant values that do not have the same ratio of −16 between loop orders.  We can also study the line on which (u, v, w) = (u, 1, 1). As in the previous subsection, hexagon functions can here be expressed in terms of harmonic polylogarithms. The symme- Thus, on this line we can express all functions in terms of E(u, 1, 1), E(1, v, 1),Ẽ(u, 1, 1), and E(u, 1, 1). We provide these functions through six loops (seven loops for E(u, 1, 1)) in the ancillary file SixGluonHPLLines.m.

The line (u, u, u) and strong coupling
Unlike the previous two lines considered, projecting to the line u = v = w does not reduce our functions to harmonic polylogarithms. Instead, representing them would require cyclotomic polylogarithms [117]. In practice, we evaluate our functions numerically on this line using series expansions around u = 0, 1, and ∞ with overlapping ranges of convergence. This line is particularly interesting because the remainder function is known here not only for weak coupling, but for strong coupling as well, due to the work of Alday, Gaiotto, and Maldacena (AGM) [102] who represented it in terms of the area of a minimal surface in AdS 5 . In fact, the area can be evaluated numerically for generic kinematics, but on this line an analytic strong-coupling formula is available: (4.39) The last term in eq. (4.39) ensures that R (∞) 6 (u, v, w) vanishes in the collinear limit. In figure 8, we plot our results for the remainder function on this line, alongside the strongcoupling AGM result. We normalize each result by its value at (1, 1, 1) for ease of comparison. This normalization forces all curves to go through unity at u = 1. Once normalized in this way, the functions are almost indistinguishable for u < 1, while for u > 1 their behavior remains rather similar, diverging from each other for large u, similarly to their behavior on other lines.
In ref. [79] Basso, Sever and Vieira (BSV) pointed out that, because certain scalars become very light at strong coupling, the determinant for fluctuations around the minimalarea surface is parametrically of the same order as the area term; that is, it is also proportional to g ∼ √ λ at strong coupling. The net result of this computation is to add an additional constant −π/72 to the AGM result: BSV argued that the collinear limit does not commute with the strong-coupling limit, as a consequence of the constant offset in eq. (4.40). The line (u, u, u) is far from the collinear limit, but we can still ask, how do large perturbative orders behave, compared to eq. (4.39) or eq. (4.40)? To do this, in fig. 9 we re-plot the perturbative results in fig. 8, normalizing everything by R (∞) 6 (u, u, u)/R (∞) 6 (1, 1, 1) and expanding the scale so we can see differences that were indistinguishable in the previous plot. The dashed line shows the strong-coupling BSV prediction in the numerator, and the same R (∞) 6 (u, u, u)/R (∞) 6 (1, 1, 1) in the denominator. Again all curves must go through unity at u = 1. For 0.4 < u < 0.1, all the perturbative results have the same shape as the strong-coupling result R (∞) 6 (u, u, u) to within a few percent. However, a naive extrapolation to yet higher loop orders suggests that the shape for 0.5 < u < 0.1 is starting to resemble the AGM+BSV prediction more closely than the AGM prediction.

The self-crossing limit
A hexagonal Wilson loop can be deformed until two of its lines, on opposite sides of the hexagon, almost cross each other. This limit was first studied for Wilson loops in planar N = 4 super-Yang-Mills theory in refs. [88][89][90]. In ref. [91], it was pointed out that the limit for 2 → 4 kinematics mimics double parton scattering in hadronic collisions, and this limit and an analogous limit for 3 → 3 kinematics were thoroughly explored. Making use of the anomalous dimension matrix for crossing Wilson lines [118], it was argued that the singular parts of E should not depend on the residual kinematic variable characterizing the self-crossing kinematics. By matching to the multi-Regge limit discussed in section 3.2, it was possible to predict the singular terms to high loop orders. In this section, we will test these predictions and provide results for related but nonsingular terms in the self-crossing limit.
In self-crossing kinematics, the cross ratios (u, v, w) approach (1 − δ, v, v) with δ → 0, after performing an analytic continuation from the Euclidean region. In the 2 → 4 case, the analytic continuation to be performed is u → e −2πi u, as in eq. (3.13); the surviving cross ratio v is restricted to 0 < v < 1, and the small parameter δ is positive. In the 3 → 3 case, v is restricted to either v < 0 or v > 1, and δ is negative, so we will express the results in terms of |δ| = −δ. For v < 0, the analytic continuation is u → e +2πi u, v → e +πi v, w → e +πi w. We can go from v < 0 to v > 1 by analytically continuing ln(1 − v) → ln(v − 1) − iπ. In fact, E 3→3 (1 + |δ|, v, v) is smooth around v = ±∞, and so no analytic continuation is necessary after writing it in the right representation in terms of harmonic polylogarithms.
The amplitude is simpler in 3 → 3 self-crossing kinematics than in 2 → 4 kinematics. The reason is that the hexagonal Wilson loop has an alternating structure, incoming-outgoingincoming-· · · , which makes it quasi-Euclidean. As a result, all the singularities as |δ| → 0 are in the imaginary part of the amplitude. In ref. [91], a formula was provided for a nonsingularlyframed Wilson loop, where L = ln(1 − 1/v) and the function g(ln |δ|, Γ cusp ) (not to be confused with the coupling constant g) was given through 7 loops. The relation between W ns and E was also determined, except that the E in ref. [91] differs from the one here by the constant factor ρ given in eq. (2.29). Correcting for that difference, we have Since ρ is a nonsingular constant, the derivative with respect to ln |δ| passes right through it, and we obtain, This formula controls all of the singular terms in E 3→3 as δ → 0. Note that these terms are totally independent of v. In ref. [91], the function g(ln |δ|, Γ cusp ) was evaluated completely through 7 loops, and partially at 8 and 9 loops, using the connection to the multi-Regge limit. However, motivated by the fact that the self-crossing limit is a virtual Sudakov region, where virtual gluons are confined to narrow jets as |δ| → 0, it is possible to give another representation for the singular terms in E 3→3 in the self-crossing limit: where J 1 is the first Bessel function and with γ E the Euler-Mascheroni constant. We have also introduced two additional anomalous dimensions, Γ virt and Γ 3 , which arise in the large-rapidity limit of the flux-tube spectrum [76,77,103,119]. Define the semi-infinite matrix K [103] and vector κ eff by and M = (1 + K) −1 , Q ij = j(−1) j+1 δ ij . Then the three anomalous dimensions are We have Γ virt = −12ζ 3 g 4 + (80ζ 5 + 16ζ 2 ζ 3 )g 6 − (700ζ 7 + 80ζ 2 ζ 5 + 168ζ 3 ζ 4 )g 8 + . . . , (5.10) . . , (5.11) while the first few terms of Γ cusp are given in eq. (2.20). Given these anomalous dimensions, it is straightforward to evaluate eq. (5.4) perturbatively to very high orders (or indeed, nonperturbatively). After expanding the integrand in g 2 , one encounters integrals of the Bessel function J 1 (2ν) multiplied by powers of ln ν + γ E . The integrals can be performed by Taylor expanding the identity [120] around ξ = 0. In an ancillary file, SelfCrossSingular.m, we evaluate eq. (5.4) through 20 loops. (We have multiplied it back by ρ, because we only know ρ to 7 loops.) In the rest of this section, we provide some information about the terms in the MHV 11 amplitude E that are nonsingular as δ → 0. The limit v → ±∞ is smooth, i.e. there are no ln v singularities for E in this limit. We will give the singular terms, as well as some nonsingular constants, in this limit. The following expressions match eqs. (3.28)-(3.33) of ref. [91] through five loops, after accounting for a factor of two difference in the loop expansion parameter, 11 The NMHV amplitude is trickier to evaluate in the self-crossing limit because some of the R-invariants diverge there, necessitating a higher-order expansion of the transcendental functions in δ, which we won't perform here. and (starting from three loops) from the factor of ρ used to define E here: As the point v = 1 is approached from above, i.e. from the 3 → 3 side, the results are also relatively simple. Because the ln |δ|-dependent terms are identical to those presented above for v = ∞, here we give only the finite terms: v → ∞ limit. There do not seem to be many restrictions on the MZVs that can appear at v = ±∞, or in the (singular) imaginary part as v → 1 + . However, the real part as v → 1 + , which is finite as δ → 0, does have MZV restrictions, and a coaction principle operates at this point as well, even though it is not on the Euclidean sheet. In fact, the dimensionality of the space of allowed MZVs here is exactly the same as for the Euclidean point (1, 1, 1) through at least weight 11, although the actual values are different [66].
We have also evaluated the MHV amplitude as a function of v for the 3 → 3 configurations with v < 0 and v > 1, and the 2 → 4 configuration with 0 < v < 1, through seven loops, in terms of harmonic polylogarithms. Because these expressions are rather lengthy at six and seven loops, where they go beyond ref. [91], we provide them in an ancillary file, SelfCross.m.

Conclusions
In this work, we have bootstrapped the six-and seven-loop six-particle MHV amplitudes, as well as the six-loop NMHV amplitude, in planar N = 4 super-Yang-Mills theory. The space of functions H hex that provides the initial ansatz for this procedure is presented in a companion paper [66]. This space of functions takes advantage of constraints following from cosmic Galois theory and the extended Steinmann relations, and appears to be the minimal space required to express the six-point amplitudes in this theory.
In order to fit the amplitude into this space, we have absorbed several constants into a function of the coupling ρ(g 2 ). This function is universal for MHV and NMHV, and it is natural to wonder whether it has physical meaning. Its number-theoretic content appears to be related to the cusp anomalous dimension, but further work will be necessary to establish whether the relation between these quantities can be made precise. Such work might hint at a physical meaning and answer the question: why is it possible to normalize these amplitudes in such a way that they respect a cosmic Galois coaction principle between loop orders?
Collinear and final-entry constraints both powerfully constrain the ansatz formed out of this space of functions. Together with information from the multi-Regge limit, they almost suffice to determine the amplitude completely. A novel feature we encountered, however, was the appearance of one potential contribution to the MHV amplitude at six loops, and again one at seven loops, that vanished in each of these limits. We successfully detected these contributions by their nonvanishing terms in the near-collinear expansion to next-to-leading order, which is governed by the Pentagon Operator Product Expansion, and in particular its first gluon bound state. It is interesting to note that while the multi-Regge limit can be reached from the OPE by an analytic continuation that also includes this state [74], it is nevertheless weaker. In other words, the analytic continuation causes some information to be lost.
Since the Pentagon OPE offers an infinite amount of boundary information, the presence of functions that cannot be detected by their multi-Regge and strict collinear limits is in principle not a problem at higher loops. In practice, accessing increasingly subleading terms in the near-collinear expansion also becomes increasingly difficult. It corresponds to acting with more and more derivatives on the transcendental functions, which in turn requires tabulating many complicated expressions. An alternative would be to partially resum a relatively simple subset of OPE excitations, that spans a more general limit where the amplitude does not vanish when evaluated there in a strict sense. The double scaling limit [80] is an excellent example. It corresponds to sending one cross ratio to zero while holding the other two cross ratios fixed. The resummation of the infinite set of purely gluonic excitations that dominate this limit was initiated in ref. [109]. Indeed, the distinct signature of OPE excitations that were resummed in the latter paper, if pursued to six and seven loops, would be able to fix the final free parameters there.
Better yet, it would be very exciting if we could harness the power of the OPE in order to access the origin in cross ratio space, where we saw that functions which can only be detected from the OPE have unusually divergent behavior, in contrast to the remainder function, which is only quadratic in the logarithms of the cross ratios through seven loops. While the origin is outside of the radius of convergence of the OPE, it certainly lies within the subspace spanned by the aforementioned double scaling limit. We are optimistic that it can be reached by means of an analytic continuation, as was the case with the multi-Regge limit.
We also studied the self-crossing limit in this paper. The singular terms as this limit is approached were computed to seven loops, and found to match previous predictions [91]. A compact, Sudakov-based formula for the singular terms was presented that can be evaluated to very high loop order. In addition we studied the nonsingular terms, which at one point are governed by the coaction principle in a nontrivial way.
We continue to observe that the amplitudes have very numerically consistent ratios between successive loop orders. These ratios approach the radius of convergence previously observed for the cusp anomalous dimension. It is interesting that the ratios of loop orders are already beginning to converge at this order, which is quite atypical for other examples of quantum field theories with finite radii of convergence.
Finally, we are optimistic that these techniques can be applied to other theories without such a high degree of symmetry.