Higher-order Sudakov resummation in coupled gauge theories

We consider the higher-order resummation of Sudakov double logarithms in the presence of multiple coupled gauge in. The associated evolution equations depend on the coupled β functions of two (or more) coupling constants α_a and α_b, as well as anomalous dimensions that have joint perturbative series in α_a and α_b. We discuss possible strategies for solving the system of evolution equations that arises. As an example, we obtain the complete three-loop (NNLL) QCD⊗QED Sudakov evolution factor. Our results also readily apply to the joint higher-order resummation of electroweak and QCD Sudakov logarithms.As part of our analysis we also revisit the case of a single gauge interaction (pure QCD), and study the numerical differences and reliability of various methods for evaluating the Sudakov evolution factor at higher orders. We find that the approximations involved in deriving commonly used analytic expressions for the evolution kernel can induce noticeable numerical differences of several percent or more at low scales, exceeding the perturbative precision at N³LL and in some cases even NNLL. Therefore, one should be cautious when using approximate analytic evolution kernels for high-precision analyses.