Related papers: Why and how to use a differential equation method …
This article displays a proof of concept of the mixed analytical/numerical method, presented in previous publications, to compute two-loop functions with up to five massive propagators in a scalar theory having three- and four-leg vertices…
The differentiable programming paradigm is a cornerstone of modern scientific computing. It refers to numerical methods for computing the gradient of a numerical model's output. Many scientific models are based on differential equations,…
The two-loop box contributions to massive Bhabha scattering may be reduced to two-loop box master integrals (MIs) with five, six, and seven internal lines, plus vertices and self energies. The self-energy and vertex MIs may be solved…
We present the analytic calculation of two-loop master integrals that are relevant for $tW$ production at hadron colliders. We focus on the integral families with only one massive propagator. After choosing a canonical basis, the…
We describe methods for evaluating one-loop integrals in $4-2\e$ dimensions. We give a recursion relation that expresses the scalar $n$-point integral as a cyclicly symmetric combination of $(n-1)$-point integrals. The computation of such…
Sunset integrals are among the simplest of two-loop integrals that appear in perturbative quantum field theories and possess up to four distinct mass scales. By means of integration by parts identities, they can be written in terms of four…
In this talk, we review recent developments towards the calculation of multi-loop scattering amplitudes. In particular, we discuss how the colour-kinematics duality can provide new integral relations at one-loop level via the Loop-Tree…
We present a method for solving a class of initial valued, coupled, non-linear differential equations with `moving singularities' subject to some subsidiary conditions. We show that this type of singularities can be adequately treated by…
We present the analytic calculation of the Master Integrals necessary to compute the planar massive QCD corrections to Di-photon (and Di-jet) production at hadron colliders. The masters are evaluated by means of the differential equations…
In this study, a recursive solution technique in conjunction with generalized integrating factors is presented and applied to address first and second order linear differential equations. This approach demonstrates practical utility in…
The Monte Carlo method is a thriving and mathematically beautiful numerical technique used extensively, nowadays, to deal with many demanding problems in diverse fields. Here, we present an iterative Monte Carlo algorithm to work out very…
We discuss the extension of the maximal-unitarity method to two loops, focusing on the example of the planar double box. Maximal cuts are reinterpreted as contour integrals, with the choice of contour fixed by the requirement that integrals…
We present analytical results for one-loop five-point master integrals with up to three off-shell legs. The method of canonical differential equations along with the Simplified Differential Equations approach is employed. All necessary…
We introduce a method for calculating individual elements of matrix functions. Our technique makes use of a novel series expansion for the action of matrix functions on basis vectors that is memory efficient even for very large matrices. We…
The $\varepsilon$-form of a system of differential equations for Feynman integrals has led to tremendeous progress in our abilities to compute Feynman integrals, as long as they fall into the class of multiple polylogarithms. It is…
A method for the evaluation of the epsilon expansion of multi-loop massless Feynman integrals is introduced. This method is based on the Gegenbauer polynomial technique and the expansion of the Gamma function in terms of harmonic sums.…
We describe the analytic calculation of the master integrals required to compute the two-mass three-loop corrections to the $\rho$ parameter. In particular, we present the calculation of the master integrals for which the corresponding…
A general procedure for the calculation of a class of two-loop Feynman diagrams is described. These are two-point functions containing three massive propagators, raised to integer powers, in the denominator, and arbitrary polynomials of the…
Loop diagram calculations typically rely on reduction to a finite set of master integrals in $4 - 2\epsilon$ dimensions. It has been shown that for any problem, the masters can be chosen so that their coefficients are finite as $\epsilon…
Estimating the trace of the inverse of a large matrix is an important problem in lattice quantum chromodynamics. A multilevel Monte Carlo method is proposed for this problem that uses different degree polynomials for the levels. The…