Total: 1
We present a basis of special functions and an efficient numerical implementation for the planar two-loop master integrals that enter the leading-colour amplitudes for $pp\to t\bar{t}H$ with a closed light-quark loop. We provide analytic results for the special functions at transcendental weight one and two in terms of logarithms and dilogarithms, while higher weights are evaluated as one-fold integrals. We focus on the nested roots that arise, studying their analytic structure on the physical region and constructing numerically stable representations of letters and $\mathrm{d}\log$ forms that involve them. We find that a number of $\mathrm{d}\log$ forms have an antiderivative that cannot be represented by a single logarithm over the physical region and tackle them numerically with adaptive Chebyshev approximation. We implement the functions in \texttt{C++} as an extension of the \texttt{PentagonFunctions++} library. We study the performance of our implementation over physical phase-space points and find that the evaluation is accurate and fast enough for phenomenological application. This provides an important building block for NNLO studies of $t\bar{t}H$ production at the LHC.