We present the version 2.0 of the program package GoSam for the automated calculation of one-loop amplitudes. GoSam is devised to compute one-loop QCD and/or electroweak corrections to multi-particle processes within and beyond the Standard Model. The new code contains improvements in the generation and in the reduction of the amplitudes, performs better in computing time and numerical accuracy, and has an extended range of applicability. The extended version of the "Binoth-Les-Houches-Accord" interface to Monte Carlo programs is also implemented. We give a detailed description of installation and usage of the code, and illustrate the new features in dedicated examples.

We show that the evaluation of scattering amplitudes can be formulated as a problem of multivariate polynomial division, with the components of the integration-momenta as indeterminates. We present a recurrence relation which, independently of the number of loops, leads to the multi-particle pole decomposition of the integrands of the scattering amplitudes. The recursive algorithm is based on the Weak Nullstellensatz Theorem and on the division modulo the Gröbner basis associated to all possible multi-particle cuts. We apply it to dimensionally regulated one-loop amplitudes, recovering the well-known integrand-decomposition formula. Finally, we focus on the maximum-cut, defined as a system of on-shell conditions constraining the components of all the integration-momenta. By means of the Finiteness Theorem and of the Shape Lemma, we prove that the residue at the maximum-cut is parametrized by a number of coefficients equal to the number of solutions of the cut itself.

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.

We present a semi-analytic method for the integrand reduction of one-loop amplitudes, based on the systematic application of the Laurent expansions to the integranddecomposition. In the asymptotic limit, the coefficients of the master integrals are the solutions of a diagonal system of equations, properly corrected by counterterms whose parametric form is known a priori. The Laurent expansion of the integrand is implemented through polynomial division. The extension of the integrand-reduction to the case of numerators with rank larger than the number of propagators is discussed as well.

We report on the calculation of the cross section for Higgs boson production in association with three jets via gluon fusion, at next-to-leading-order (NLO) accuracy in QCD, in the infinite top-mass approximation. After including the complete NLO QCD corrections, we observe a strong reduction in the scale dependence of the result, and an increased steepness in the transverse momentum distributions of both the Higgs and the leading jets. The results are obtained with the combined use of GoSam, Sherpa, and the MadDipole/MadEvent framework.

Elaborating on the four-dimensional helicity scheme, we propose a pure four-dimensional formulation (FDF) of the d-dimensional regularization of one-loop scattering amplitudes. In our formulation particles propagating inside the loop are represented by massive internal states regulating the divergences. The latter obey Feynman rules containing multiplicative selection rules which automatically account for the effects of the extra-dimensional regulating terms of the amplitude. We present explicit representations of the polarization and helicity states of the four-dimensional particles propagating in the loop. They allow for a complete, four-dimensional, unitarity-based construction of d-dimensional amplitudes. Generalized unitarity within the FDF does not require any higher-dimensional extension of the Clifford and the spinor algebra. Finally we show how the FDF allows for the recursive construction of d-dimensional one-loop integrands, generalizing the fourdimensional open-loop approach.

We present the integrand reduction via multivariate polynomial division as a natural technique to encode the unitarity conditions of Feynman amplitudes. We derive a recursive formula for the integrand reduction, valid for arbitrary dimensionally regulated loop integrals with any number of loops and external legs, which can be used to obtain the decomposition of any integrand analytically with a finite number of algebraic operations. The general results are illustrated by applications to two-loop Feynman diagrams in QED and QCD, showing that the proposed reduction algorithm can also be seamlessly applied to integrands with denominators appearing with arbitrary powers. arXiv:1307.5832v2 [hep-ph]

We present the calculation of the NLO QCD corrections to the associated production of a Higgs boson and two jets, in the infinite top-mass limit. We discuss the technical details of the computation and we show the numerical impact of the radiative corrections on several observables at the LHC. The results are obtained by using a fully automated framework for fixed order NLO QCD calculations based on the interplay of the packages GoSam and Sherpa. The evaluation of the virtual corrections constitutes an application of the d-dimensional integrand-level reduction to theories with higher dimensional operators. We also present first results for the one-loop matrix elements of the partonic processes with a quark-pair in the final state, which enter the hadronic production of a Higgs boson together with three jets in the infinite top-mass approximation.