We report on the first calculation of next-to-next-to-leading order (NNLO) QCD corrections to the inclusive production of ZZ pairs at hadron colliders. Numerical results are presented for pp collisions with centre-of-mass energy ($\sqrt{s}$) ranging from 7 to 14 TeV. The NNLO corrections increase the NLO result by an amount varying from $11\%$ to $17\%$ as $\sqrt{s}$ goes from 7 to 14 TeV. The loop-induced gluon fusion contribution provides about $60\%$ of the total NNLO effect. When going from NLO to NNLO the scale uncertainties do not decrease and remain at the $\pm 3\%$ level.Comment: Reference added, version published on Physics Letters

Charged gauge boson pair production at the Large Hadron Collider allows detailed probes of the fundamental structure of electroweak interactions. We present precise theoretical predictions for on-shell W+ W- production that include, for the first time, QCD effects up to next to next to leading order in perturbation theory. As compared to next to leading order, the inclusive W+ W- cross section is enhanced by 9% at 7 TeV and 12% at 14 TeV. The residual perturbative uncertainty is at the 3% level. The severe contamination of the W+ W- cross section due to top-quark resonances is discussed in detail. Comparing different definitions of top-free W+ W- production in the four and five flavor number schemes, we demonstrate that top-quark resonances can be separated from the inclusive W+ W- cross section without a significant loss of theoretical precision.

We introduce a class of iterated integrals, defined through a set of linearly independent integration kernels on elliptic curves. As a direct generalisation of multiple polylogarithms, we construct our set of integration kernels ensuring that they have at most simple poles, implying that the iterated integrals have at most logarithmic singularities. We study the properties of our iterated integrals and their relationship to the multiple elliptic polylogarithms from the mathematics literature. On the one hand, we find that our iterated integrals span essentially the same space of functions as the multiple elliptic polylogarithms. On the other, our formulation allows for a more direct use to solve a large variety of problems in high-energy physics. We demonstrate the use of our functions in the evaluation of the Laurent expansion of some hypergeometric functions for values of the indices close to half integers.

Abstract:We compute the full set of two-loop Feynman integrals appearing in massless two-loop four-point functions with two off-shell legs with the same invariant mass. These integrals allow to determine the two-loop corrections to the amplitudes for vector boson pair production at hadron colliders, qq → V V , and thus to compute this process to next-tonext-to-leading order accuracy in QCD. The master integrals are derived using the method of differential equations, employing a canonical basis for the integrals. We obtain analytical results for all integrals, expressed in terms of multiple polylogarithms. We optimize our results for numerical evaluation by employing functions which are real valued for physical scattering kinematics and allow for an immediate power series expansion.

Abstract:We elaborate on the method of differential equations for evaluating Feynman integrals. We focus on systems of equations for master integrals having a linear dependence on the dimensional parameter. For these systems we identify the criteria to bring them in a canonical form, recently identified by Henn, where the dependence of the dimensional parameter is disentangled from the kinematics. The determination of the transformation and the computation of the solution are obtained by using Magnus and Dyson series expansion. We apply the method to planar and non-planar two-loop QED vertex diagrams for massive fermions, and to non-planar two-loop integrals contributing to 2 → 2 scattering of massless particles. The extension to systems which are polynomial in the dimensional parameter is discussed as well.