We use hexagonalization to compute four-point correlation functions of long BPS operators with special R-charge polarizations. We perform the computation at weak coupling and show that at any loop order our correlators can be expressed in terms of single value polylogarithms with uniform and maximal transcendentality. As a check of our computation we extract nine-loop OPE data and compare it against sum rules of (squared) structures constants of non-protected exchanged operators described by hundreds of Bethe solutions.

We present the full form of a four-point correlation function of large BPS operators in planar N = 4 Super Yang-Mills to any loop order. We do this by following a bootstrap philosophy based on three simple axioms pertaining to (i) the space of functions arising at each loop order, (ii) the behaviour in the OPE in a double-trace dominated channel and (iii) the behaviour under a double null limit. We discuss how these bootstrap axioms are in turn strongly motivated by empirical observations up to nine loops unveiled through integrability methods in our previous work [9] on this simplest correlation function.

We initiate the study of four-point functions of large BPS operators at any value of the coupling. We do it by casting it as a sum over exchange of superconformal primaries and computing the structure constants using integrability. Along the way, we incorporate the nested Bethe ansatz structure to the hexagon formalism for the three-point functions and obtain a compact formula for the asymptotic structure constant of a non-BPS operator in a higher rank sector.