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.