We extend the auxiliary-mass-flow (AMF) method originally developed for Feynman loop integration to calculate integrals involving also phase-space integration. Flow of the auxiliary mass from the boundary ($infty$) to the physical point ($0^+$) is obtained by numerically solving differential equations with respective to the auxiliary mass. For problems with two or more kinematical invariants, the AMF method can be combined with traditional differential equation method by providing systematical boundary conditions and highly nontrivial self-consistent check. The method is described in detail with a pedagogical example of $e^+e^-rightarrow gamma^* rightarrow tbar{t}+X$ at NNLO. We show that the AMF method can systematically and efficiently calculate integrals to high precision.
One approach to the calculation of cross sections for infrared-safe observables in high energy collisions at next-to-leading order is to perform all of the integrations, including the virtual loop integration, by Monte Carlo numerical integration. In a previous paper, two of us have shown how one can perform such a virtual loop integration numerically after first introducing a Feynman parameter representation. In this paper, we perform the integration directly, without introducing Feynman parameters, after suitably deforming the integration contour. Our example is the N-photon scattering amplitude with a massless electron loop. We report results for N = 6 and N = 8.
We show that direct Feynman-parametric loop integration is possible for a large class of planar multi-loop integrals. Much of this follows from the existence of manifestly dual-conformal Feynman-parametric representations of planar loop integrals, and the fact that many of the algebraic roots associated with (e.g. Landau) leading singularities are automatically rationalized in momentum-twistor space---facilitating direct integration via partial fractioning. We describe how momentum twistors may be chosen non-redundantly to parameterize particular integrals, and how strategic choices of coordinates can be used to expose kinematic limits of interest. We illustrate the power of these ideas with many concrete cases studied through four loops and involving as many as eight particles. Detailed examples are included as ancillary files to this works submission to the arXiv.
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 computer 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.
In this manuscript, which is to appear in the proceedings of the conference MathemAmplitude 2019 in Padova, Italy, we provide an overview of the module intersection method for the the integration-by-parts (IBP) reduction of multi-loop Feynman integrals. The module intersection method, based on computational algebraic geometry, is a highly efficient way of getting IBP relations without double propagator or with a bound on the highest propagator degree. In this manner, trimmed IBP systems which are much shorter than the traditional ones can be obtained. We apply the modern, Petri net based, workflow management system GPI-Space in combination with the computer algebra system Singular to solve the trimmed IBP system via interpolation and efficient parallelization. We show, in particular, how to use the new plugin feature of GPI-Space to manage a global state of the computation and to efficiently handle mutable data. Moreover, a Mathematica interface to generate IBPs with restricted propagator degree, which is based on module intersection, is presented in this review.
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.