“…Furthermore, we implement our expressions in a publicly available C++ library that allows one to evaluate these contributions efficiently and in a stable manner, making our results readily available for a future full calculation of the corresponding NNLO QCD corrections. 1 Our work builds on a number of recent technical advances that have paved the way to the calculation of two-loop scattering amplitudes for 2 → 3 processes with phenomenological applications [32][33][34][35][36][37]. We employ the known pure two-loop master integrals [21,26,38,39], which can be written in terms of a special set of multi-valued transcendental functions -the so-called pentagon functions [40,41].…”