A purely numerical method, Direct ComputationMethod is applied to evaluate Feynman integrals. This method is based on the combination of an efficient numerical integration and an efficient extrapolation. In addition, high-precision arithmetic and parallelization technique can be used in this method if required. We present the recent progress in development of this method and show results such as one-loop 5-point and two-loop 3-point integrals.
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.
A new approach to compute Feynman Integrals is presented. It relies on an integral representation of a given Feynman Integral in terms of simpler ones. Using this approach, we present, for the first time, results for a certain family of non-planar five-point two-loop Master Integrals with one external off-shell particle, relevant for instance for $H+2$ jets production at the LHC, in both Euclidean and physical kinematical regions.
We present numerical results which are needed to evaluate all non-trivial master integrals for four-loop massless propagators, confirming the recent analytic results of[1]and evaluating an extra order in $ep$ expansion for each master integral.
We present pySecDec, a new version of the program SecDec, which performs the factorisation of dimensionally regulated poles in parametric integrals, and the subsequent numerical evaluation of the finite coefficients. The algebraic part of the program is now written in the form of python modules, which allow a very flexible usage. The optimization of the C++ code, generated using FORM, is improved, leading to a faster numerical convergence. The new version also creates a library of the integrand functions, such that it can be linked to user-specific codes for the evaluation of matrix elements in a way similar to analytic integral libraries.
A method for reducing Feynman integrals, depending on several kinematic variables and masses, to a combination of integrals with fewer variables is proposed. The method is based on iterative application of functional equations proposed by the author. The reduction of the one-loop scalar triangle and box integrals with massless internal propagators to simpler integrals is described in detail. The triangle integral depending on three variables is represented as a sum over three integrals depending on two variables. By solving the dimensional recurrence relations for these integrals, an analytic expression in terms of the $_2F_1$ Gauss hypergeometric function and the logarithmic function was derived. By using the functional equations, the one-loop box integral with massless internal propagators, which depends on six kinematic variables, was expressed as a sum of 12 terms. These terms are proportional to the same integral depending only on three variables different for each term. For this integral with three variables, an analytic result in terms of the $F_1$ Appell and $_2F_1$ Gauss hypergeometric functions was derived by solving the recurrence relation with respect to the spacetime dimension $d$. The reduction equations for the box integral with some kinematic variables equal to zero are considered.