We present a systematic algorithm for the perturbative computation of soft functions that are defined in terms of two light-like Wilson lines. Our method is based on a universal parametrisation of the phase-space integrals, which we use to isolate the singularities in Laplace space. The observable-dependent integrations can then be performed numerically, and they are implemented in the new, publicly available package SoftSERVE that we use to derive all of our numerical results. Our algorithm applies to both SCET-1 and SCET-2 soft functions, and in the current version it can be used to compute two out of three NNLO colour structures associated with the so-called correlated-emission contribution. We confirm existing two-loop results for about a dozen e + e − and hadron-collider soft functions, and we obtain new predictions for the C-parameter as well as thrust-axis and broadening-axis angularities.

We present compact integral representations for the calculation of two-loop anomalous dimensions for a generic class of soft functions that are defined in terms of two light-like Wilson lines. Our results are relevant for the resummation of Sudakov logarithms for e + e − event-shape variables and inclusive hadron-collider observables at next-to-next-toleading logarithmic accuracy within Soft-Collinear Effective Theory (SCET). Our formalism applies to both SCET-1 and SCET-2 soft functions and we clarify the relation between the respective soft anomalous dimension and the collinear anomaly exponent. We confirm existing two-loop results for about a dozen dijet soft functions and obtain new predictions for the angularity event shape and the soft-drop jet-grooming algorithm.

We extend our algorithm for automating the calculation of two-loop dijet soft functions to observables that do not obey the non-Abelian exponentiation theorem, i.e. to those that require an independent calculation of the uncorrelated-emission contribution. As the singularity structure of uncorrelated double emissions differs substantially from the one for correlated emissions, we introduce a novel phase-space parametrisation that isolates the corresponding divergences. The resulting integrals are implemented in SoftSERVE 1.0, which we release alongside of this work, and which we supplement by a regulator that is consistent with the rapidity renormalisation group framework. Using our automated setup, we confirm existing results for various jet-veto observables and provide a novel prediction for the soft-drop jet-grooming algorithm.