Numerical scattering amplitudes with pySecDec

Abstract

We present a major update of the program pySecDec, a toolbox for the evaluation of dimensionally regulated parameter integrals. The new version enables the evaluation of multi-loop integrals as well as amplitudes in a highly distributed and flexible way, optionally on GPUs. The program has been optimised and runs up to an order of magnitude faster than the previous release. A new integration procedure that utilises construction-free median Quasi-Monte Carlo rules is implemented. The median lattice rules can outperform our previous component-by-component rules by a factor of 5 and remove the limitation on the maximum number of sampling points. The expansion by regions procedures have been extended to support Feynman integrals with numerators, and functions for automatically determining when and how analytic regulators should be introduced are now available. The new features and performance are illustrated with several examples. Program summaryProgram Title: pySecDecCPC Library link to program files:https://doi.org/10.17632/dnrkf5jxzh.4Developer's repository link:https://github.com/gudrunhe/secdec, https://secdec.readthedocs.io (Online documentation)Licensing provisions: GNU Public License v3Programming language: Python, Form, C++, CudaExternal routines/libraries:GSL [1], NumPy [2], SymPy [3], Nauty [4], Cuba [5], Form [6], GiNaC and CLN [7], Normaliz [8], GMP [9].Journal reference of previous version: Comput. Phys. Commun. 273 (2022) 108267 [1].Does the new version supersede the previous version?: Yes.Nature of problem: Scattering amplitudes at higher orders in perturbation theory are typically represented as a linear combination of coefficients — containing the kinematic invariants and the space-time dimension — multiplied with loop integrals which contain singularities and whose analytic representation might be unknown.Solution method: Extraction of singularities in the dimensional regularization parameter as well as in analytic regulators for potential spurious singularities is done using sector decomposition. The combined evaluation of the integrals with their coefficients is performed in an efficient way.Restrictions: Depending on the complexity of the problem, limited by memory and CPU/GPU time.