We present the first complete next-to-next-to-leading order (NNLO) QCD predictions for differential distributions in the top-quark pair production process at the LHC. Our results are derived from a fully differential partonic Monte Carlo calculation with stable top quarks which involves no approximations beyond the fixed-order truncation of the perturbation series. The NNLO corrections improve the agreement between existing LHC measurements [V. Khachatryan et al. (CMS Collaboration), Eur. Phys. J. C 75, 542 (2015)] and standard model predictions for the top-quark transverse momentum distribution, thus helping alleviate one long-standing discrepancy. The shape of the top-quark pair invariant mass distribution turns out to be stable with respect to radiative corrections beyond NLO which increases the value of this observable as a place to search for physics beyond the standard model. The results presented here provide essential input for parton distribution function fits, implementation of higher-order effects in Monte Carlo generators, as well as top-quark mass and strong coupling determination.

Four years ago, one of us introduced a novel subtraction scheme [1] for the evaluation of double-real radiation contributions to cross sections at next-to-next-to-leading order (NNLO) in QCD. This approach, named SecToR Improved Phase sPacE for Real radiation (STRIPPER), has already found several nontrivial applications. In particular, it has allowed for the determination of NNLO corrections to hadronic top-quark pair production, fully differential top-quark decays, inclusive semileptonic charmless b-quark decays, associated Higgs boson and jet production in gluon fusion, muon decay spin asymmetry, and tchannel single-top production. Common to these calculations was the use of conventional dimensional regularization (CDR). In this publication, we present a complete formulation of the subtraction scheme for arbitrary processes with any number of colored partons in the final state, and up to two partons in the initial state. Furthermore, we modify the integrated subtraction terms of the double-real radiation to enable the introduction of the 't Hooft-Veltman version of dimensional regularization (HV), in which resolved states are four-dimensional. We demonstrate the correctness of our approach on the example of top-quark pair production in the gluon fusion channel.

In this work we present for the first time predictions for top-quark pair differential distributions at the LHC at NNLO QCD accuracy and including EW corrections. For the latter we include not only contributions of O(α 2 s α), but also those of order O(α s α 2 ) and O(α 3 ). Besides providing phenomenological predictions for all main differential distributions with stable top quarks, we also study the following issues. 1) The effect of the photon PDF on top-pair spectra: we find it to be strongly dependent on the PDF set used -especially for the top p T distribution. 2) The difference between the additive and multiplicative approaches for combining QCD and EW corrections: with our scale choice, we find relatively small differences between the central predictions, but reduced scale dependence within the multiplicative approach.3) The potential effect from the radiation of heavy bosons on inclusive top-pair spectra: we find it to be, typically, negligible.

We calculate all major differential distributions with stable top-quarks at the LHC. The calculation covers the multi-TeV range that will be explored during LHC Run II and beyond. Our results are in the form of high-quality binned distributions. We offer predictions based on three different parton distribution function (pdf) sets. In the near future we will make our results available also in the more flexible fastNLO format that allows fast re-computation with any other pdf set. In order to be able to extend our calculation into the multi-TeV range we have had to derive a set of dynamic scales. Such scales are selected based on the principle of fastest perturbative convergence applied to the differential and inclusive cross-section. Many observations from our study are likely to be applicable and useful to other precision processes at the LHC. With scale uncertainty now under good control, pdfs arise as the leading source of uncertainty for TeV top production. Based on our findings, true precision in the boosted regime will likely only be possible after new and improved pdf sets appear. We expect that LHC top-quark data will play an important role in this process.

We demonstrate that a purposefully normalised NNLO m tt differential spectrum can have very small theoretical uncertainty and, in particular, a small sensitivity to the top quark mass. Such observable can thus be a very effective bump-hunting tool for resonances decaying to tt events during LHC Run II and beyond. To illustrate how the approach works, we concentrate on one specific example of current interest, namely, the possible 750 GeV di-gamma excess resonance Φ. Considering only theoretical uncertainties, we demonstrate that it is possible to distinguish pp → Φ → tt signals studied in the recent literature [Hespel, Maltoni and Vryonidou, arXiv:1606.04149] from the pure SM background with very high significance. Alternatively, in case of non-observation, a strong upper limit on the decay rate Φ → tt can be placed.