We present the new version of OpenLoops, an automated generator of tree and one-loop scattering amplitudes based on the open-loop recursion. One main novelty of OpenLoops 2 is the extension of the original algorithm from NLO QCD to the full Standard Model, including electroweak (EW) corrections from gauge, Higgs and Yukawa interactions. In this context, among several new features, we discuss the systematic bookkeeping of QCD-EW interferences, a flexible implementation of the complex-mass scheme for processes with on-shell and off-shell unstable particles, a special treatment of on-shell and off-shell external photons, and efficient scale variations. The other main novelty is the implementation of the recently proposed on-the-fly reduction algorithm, which supersedes the usage of external reduction libraries for the calculation of tree-loop interferences. This new algorithm is equipped with an automated system that avoids Gram-determinant instabilities through analytic methods in combination with a new hybrid-precision approach based on a highly targeted usage of quadruple precision with minimal CPU overhead. The resulting significant speed and stability improvements are especially relevant for challenging NLO multi-leg calculations and for NNLO applications. only little user intervention. Moreover, OpenLoops is used as a building block of Matrix [50] for the calculation of NNLO QCD observables. In this context, the automation of EW corrections in OpenLoops 2 opens the door to ubiquitous NLO QCD+NLO EW simulations in Sherpa [51, 52] and NNLO QCD+NLO EW calculations in Matrix [53].The OpenLoops 2 code is publicly available on the Hepforge webpage https://openloops.hepforge.org and via the Git repository https://gitlab.com/openloops/OpenLoops. It consists of a processindependent base code and a process library that covers several hundred partonic processes, including essentially all relevant processes at the LHC. The desired processes can be easily accessed through an automated download mechanism. The set of available processes is continuously extended, and possible missing processes can be promptly generated by the authors upon request.The paper is organised as follows. Section 2 presents the structure of the original open-loop recursion and the new on-the-fly reduction algorithm. Numerical instabilities and the new hybrid-precision system are discussed in detail. Section 3 deals with general aspects of NLO calculations and their automation in OpenLoops. This includes the bookkeeping of towers of terms of variable order α p s α q , the treatment of input parameters, optimal couplings for external photons, the renormalisation of the SM at O(α s ) and O(α), the on-shell and complex-mass schemes, and the I-operator. Section 4 provides instructions on how to use the program, starting from installation and process selection, and including the various interfaces for the calculation of matrix elements, colour/spin correlators, and tree amplitudes in colour space. Technical benchmarks concerning the speed and numerical sta...

We analytically compute the dominant contributions to the β-functions for the top-Yukawa coupling, the strong coupling and the Higgs self-coupling as well as the anomalous dimensions of the scalar, gluon and quark fields in the unbroken phase of the Standard Model at threeloop level. These are mainly the QCD and top-Yukawa corrections. The contributions from the Higgs self-interaction which are negligible for the running of the top-Yukawa and the strong coupling but important for the running of the Higgs self-coupling are also evaluated.Recent exciting evidence from several SM-like Higgs search channels at both the CERN Large Hadron Collider and the Fermilab Tevatron [6-8] point to the possibility of a SM Higgs boson with a mass in the vicinity of 125 GeV which is in truly remarkable agreement with the aforementioned prediction. 1 This calls for more precise calculations, in particular of β-functions, in the SM. In the present paper we are particularly interested in the evolution of the Higgs self-coupling as well as the top-Yukawa coupling in the SM.1 Note that the boundary condition λ(M P lanck ) = 0, leading to the prediction of the Higgs mass close the the experimental evidence, has been also discussed recently in [9].

Building on the open-loop algorithm we introduce a new method for the automated construction of oneloop amplitudes and their reduction to scalar integrals. The key idea is that the factorisation of one-loop integrands in a product of loop segments makes it possible to perform various operations on-the-fly while constructing the integrand. Reducing the integrand on-the-fly, after each segment multiplication, the construction of loop diagrams and their reduction are unified in a single numerical recursion. In this way we entirely avoid objects with high tensor rank, thereby reducing the complexity of the calculations in a drastic way. Thanks to the on-the-fly approach, which is applied also to helicity summation and for the merging of different diagrams, the speed of the original open-loop algorithm can be further augmented in a very significant way. Moreover, addressing spurious singularities of the employed reduction identities by means of simple expansions in rank-two Gram determinants, we achieve a remarkably high level of numerical stability. These features of the new algorithm, which will be made publicly available in a forthcoming release of the OpenLoops program, are particularly attractive for NLO multi-leg and NNLO real-virtual calculations.

The advent of efficient numerical algorithms for the construction of one-loop amplitudes has played a crucial role in the automation of NLO calculations, and the development of similar algorithms at two loops is a natural strategy for NNLO automation. Within a numerical framework the numerator of loop integrals is usually constructed in four dimensions, and the missing rational terms, which arise from the interplay of the (D − 4)-dimensional parts of the loop numerator with 1/(D − 4) poles in D dimensions, are reconstructed separately. At one loop, such rational terms arise only from UV divergences and can be restored through process-independent local counterterms. In this paper we investigate the behaviour of rational terms of UV origin at two loops. The main result is a general formula that combines the subtraction of UV poles with the reconstruction of the associated rational parts at two loops. This formula has the same structure as the R-operation, and all poles and rational parts are described through a finite set of processindependent local counterterms. We also present a general formula for the calculation of all relevant two-loop rational counterterms in any renormalisable theory based on onescale tadpole integrals. As a first application, we derive the full set of two-loop rational counterterms for QED in the R ξ-gauge.

We analytically calculate higher order corrections to coefficient functions of the operator product expansion (OPE) for the Euclidean correlator of two energy-momentum tensors in massless QCD. These are the three-loop contribution to the coefficient C_0 in front of the unity operator O_0=1 and the one and two-loop contributions to the coefficient C_1 in front of the gluon "condensate" operator O_1=-1/4 G^{\mu \nu} G_{\mu \nu}. For the correlator of two operators O_1 we present the coefficient C_1 at two-loop level.Comment: v2: comments and appendix added, results are now available in computer readable form; v3: JHEP version, extended discussion of the method of projector