We report on the progress in constructing contracted one-loop tensors. Analytic results for rank R=4 tensors, cross-checked numerically, are presented for the first time.
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 dimen
sion $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$.
In this paper, we describe a numerical approach to evaluate Feynman loop integrals. In this approach the key technique is a combination of a numerical integration method and a numerical extrapolation method. Since the computation is carried out in a
fully numerical way, our approach is applicable to one-, two- and multi-loop diagrams. Without any analytic treatment it can compute diagrams with not only real masses but also complex masses for the internal particles. As concrete examples we present numerical results of a scalar one-loop box integral with complex masses and two-loop planar and non-planar box integrals with masses. We discuss the quality of our numerical computation by comparisons with other methods and also propose a self consistency check.
The long-standing problem of representing the general massive one-loop Feynman integral as a meromorphic function of the space-time dimension $d$ has been solved for the basis of scalar one- to four-point functions with indices one. In 2003 the solut
ion of difference equations in the space-time dimension allowed to determine the necessary classes of special functions: self-energies need ordinary logarithms and Gauss hypergeometric functions $_2F_1$, vertices need additionally Kamp{e} de F{e}riet-Appell functions $F_1$, and box integrals also Lauricella-Saran functions $F_S$. In this study, alternative recursive Mellin-Barnes representations are used for the representation of $n$-point functions in terms of $(n-1)$-point functions. The approach enabled the first derivation of explicit solutions for the Feynman integrals at arbitrary kinematics. In this article, we scetch our new representations for the general massive vertex and box Feynman integrals and derive a numerical approach for the necessary Appell functions $F_1$ and Saran functions $F_S$ at arbitrary kinematical arguments.
We introduce an algebro-geometrically motived integration-by-parts (IBP) reduction method for multi-loop and multi-scale Feynman integrals, using a framework for massively parallel computations in computer algebra. This framework combines the compute
r algebra system Singular with the workflow management system GPI-Space, which is being developed at the Fraunhofer Institute for Industrial Mathematics (ITWM). In our approach, the IBP relations are first trimmed by modern algebraic geometry tools and then solved by sparse linear algebra and our new interpolation methods. These steps are efficiently automatized and automatically parallelized by modeling the algorithm in GPI-Space using the language of Petri-nets. We demonstrate the potential of our method at the nontrivial example of reducing two-loop five-point nonplanar double-pentagon integrals. We also use GPI-Space to convert the basis of IBP reductions, and discuss the possible simplification of IBP coefficients in a uniformly transcendental basis.
We discuss briefly the first numerical implementation of the Loop-Tree Duality (LTD) method. We apply the LTD method in order to calculate ultraviolet and infrared finite multi-leg one-loop Feynman integrals. We attack scalar and tensor integrals wit
h up to six legs (hexagons). The LTD method shows an excellent performance independently of the number of external legs.