The planar two-loop scalar Feynman integrals contributing to the massive NNLO QCD corrections for W-boson pair production via quark-antiquark annihilation can be classified into three family branches, each of which is reduced to a distinct set of master integrals (MIs), totaling 27, 45 and 15, respectively. These MIs are analytically calculated using the method of differential equations, with solutions expanded as Taylor series in the dimensional regulator ϵ. For the first two family branches, the differential systems can be successfully transformed into canonical form by adopting appropriate bases of MIs. This enables the MIs of these family branches to be expressed either as Goncharov polylogarithms (GPLs) or as one-fold integrals over GPLs, up to O(ϵ4). In contrast, the differential system for the third family branch can only be cast into a form linear in ϵ due to the presence of elliptic integrals. The solution to this linear-form differential system is expressed in an iterated form owing to the strictly lower-triangular structure of the coefficient matrices at ϵ = 0. Our analytic expressions for these MIs are verified with high accuracy against the numerical results from the AMFlow package.
Read full abstract