We briefly sketch the methods for a numerically stable evaluation of tensor one-loop integrals that have been used in the calculation of the complete electroweak one-loop corrections to $PepPemto4 $fermions. In particular, the improvement of the new methods over the conventional Passarino--Veltman reduction is illustrated for some 4-point integrals in the delicate limits of small Gram (and other kinematical) determinants.
We show how to evaluate tensor one-loop integrals in momentum space avoiding the usual plague of Gram determinants. We do this by constructing combinations of $n$- and $(n-1)$-point scalar integrals that are finite in the limit of vanishing Gram determinant. These non-trivial combinations of dilogarithms, logarithms and constants are systematically obtained by either differentiating with respect to the external parameters - essentially yielding scalar integrals with Feynman parameters in the numerator - or by developing the scalar integral in $D=6-2e$ or higher dimensions. As an explicit example, we develop the tensor integrals and associated scalar integral combinations for processes where the internal particles are massless and where up to five (four massless and one massive) external particles are involved. For more general processes, we present the equations needed for deriving the relevant combinations of scalar integrals.
We present the computation of a full set of planar five-point two-loop master integrals with one external mass. These integrals are an important ingredient for two-loop scattering amplitudes for two-jet-associated W-boson production at leading color in QCD. We provide a set of pure integrals together with differential equations in canonical form. We obtain analytic differential equations efficiently from numerical samples over finite fields, fitting an ansatz built from symbol letters. The symbol alphabet itself is constructed from cut differential equations and we find that it can be written in a remarkably compact form. We comment on the analytic properties of the integrals and confirm the extended Steinmann relations, which govern the double discontinuities of Feynman integrals, to all orders in $epsilon$. We solve the differential equations in terms of generalized power series on single-parameter contours in the space of Mandelstam invariants. This form of the solution trivializes the analytic continuation and the integrals can be evaluated in all kinematic regions with arbitrary numerical precision.
We present a new algorithm for the reduction of one-loop emph{tensor} Feynman integrals with $nleq 4$ external legs to emph{scalar} Feynman integrals $I_n^D$ with $n=3,4$ legs in $D$ dimensions, where $D=d+2l$ with integer $l geq 0$ and generic dimension $d=4-2epsilon$, thus avoiding the appearance of inverse Gram determinants $()_4$. As long as $()_4 eq 0$, the integrals $I_{3,4}^D$ with $D>d$ may be further expressed by the usual dimensionally regularized scalar functions $I_{2,3,4}^d$. The integrals $I_{4}^D$ are known at $()_4 equiv 0$, so that we may extend the numerics to small, non-vanishing $()_4$ by applying a dimensional recurrence relation. A numerical example is worked out. Together with a recursive reduction of 6- and 5-point functions, derived earlier, the calculational scheme allows a stabilized reduction of $n$-point functions with $nleq 6$ at arbitrary phase space points. The algorithm is worked out explicitely for tensors of rank $Rleq n$.
We construct a specific formalism for calculating the one-loop virtual corrections for standard model processes with an arbitrary number of external legs. The procedure explicitly separates the infrared and ultraviolet divergences analytically from the finite one-loop contributions, which can then be evaluated numerically using recursion relations. Using the formalism outlined in this paper, we are in position to construct the next-to-leading order corrections to a variety of multi-leg QCD processes such as multi-jet production and vector-boson(s) plus multi-jet production at hadron colliders. The final limiting factor on the number of particles will be the available computer power.
A comprehensive study is performed of general massive, tensor, two-loop Feynman diagrams with two and three external legs. Reduction to generalized scalar functions is discussed. Integral representations, supporting the same class of smoothness algorithms already employed for the numerical evaluation of ordinary scalar functions, are introduced for each family of diagrams.