We present Fuchsia -an implementation of the Lee algorithm, which for a given system of ordinary differential equations with rational coefficients

No abstract

No abstract

Among the unitarity cuts of 4-loop massless propagators two kinds are currently fully known: the 2-particle cuts with 3 loops corresponding to form-factors, and the 5-particle phase-space integrals. In this paper we calculate master integrals for the remaining ones: 3-particle cuts with 2 loops, and 4-particle cuts with 1 loop. The 4-particle cuts are calculated by solving dimensional recurrence relations. The 3-particle cuts are integrated directly using 1→3 amplitudes with 2 loops, which we also re-derive here up to transcendentality weight 7. The results are verified both numerically, and by showing consistency with previously known integrals using Cutkosky rules. We provide the analytic results in the space-time dimension 4 − 2ε as series in ε with coefficients being multiple zeta values up to weight 12. In the ancillary files we also provide dimensional recurrence matrices and SummerTime files suitable for numerical evaluation of the series in arbitrary dimensions with any precision.

We report on the recent progress in reducing differential equations for Feynman master integrals to canonical form with the help of a method proposed by Roman Lee. For the first time, we present Fuchsia -our open-source implementation of the Lee algorithm written in Python using mathematical routines of a free computer algebra system SageMath. We demonstrate Fuchsia by reducing differential equations for NLO contributions to splitting functions in QCD, which contain both loops and legs integrals.