All-loop four-quark Bethe-Salpeter kernel

We analytically calculate the all-loop bare perturbative part of the four-quark Bethe-Salpeter kernel using modern scattering amplitude methods. We work to subleading order in the large number of quark flavors approximation of massless Quantum Chromodynamics, which simultaneously makes an all-loop calculation feasible, is systematically improvable, and preserves asymptotic freedom. It also allows for avoiding the ambiguity of choosing a truncation scheme in Dyson-Schwinger equations. We exploit state-of-the-art methods in Integration-By-Parts reduction of Lorentz scalar Feynman integrals into a minimal Master Integral basis, and direct integration into Generalized Polylogarithms. As a byproduct of our calculation, we also provide the result for the gluon and quark propagators. We discuss a path towards nonperturbative formulation and potential future phenomenological applications.