pySecDec: A toolbox for the numerical evaluation of multi-scale integrals

By Sophia Borowka, Gudrun Heinrich, Stephan Jahn, Stephen P. Jones, Matthias Kerner, Johannes Schlenk, Tom J. E. Zirke. Computer Physics Communications, 2017

316 citations recorded by the source corpus.

Abstract We present py SecDec , a new version of the program SecDec , which performs the factorization 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. Program summary Program Title: pySecDec Program Files doi: http://dx.doi.org/10.17632/3y8bbz9c9v.1 Licensing provisions: GNU Public License v3 Programming language: python, FORM, C++ External routines/libraries: catch [1], gsl [2], numpy [3], sympy [4], Nauty [5], Cuba [6], FORM [7], Normaliz [8]. The program can also be used in a mode which does not require Normaliz. Journal reference of previous version: Comput. Phys. Commun. 196 (2015) 470–491. Nature of the problem: Extraction of ultraviolet and infrared singularities from parametric integrals appearing in higher order perturbative calculations in quantum field theory. Numerical integration in the presence of integrable singularities (e.g. kinematic thresholds). Solution method: Algebraic extraction of singularities within dimensional regularization using iterated sector decomposition. This leads to a Laurent series in the dimensional regularization parameter ϵ (and optionally other regulators), where the coefficients are finite integrals over the unit-hypercube. Those integrals are evaluated numerically by Monte Carlo integration. The integrable singularities are handled by choosing a suitable integration contour in the complex plane, in an automated way. The parameter integrals forming the coefficients of the Laurent series in the regulator(s) are provided in the form of libraries which can be linked to the calculation of (multi-) loop amplitudes. Restrictions: Depending on the complexity of the problem, limited by memory and CPU time. References: [1] https://github.com/philsquared/Catch/ . [2] http://www.gnu.org/software/gsl/ . [3] http://www.numpy.org/ . [4] http://www.sympy.org/ . [5] http://pallini.di.uniroma1.it/ . [6] T. Hahn, “CUBA: A Library for multidimensional numerical integration,” Comput. Phys. Commun. 168 (2005) 78 [hep-ph/0404043], http://www.feynarts.de/cuba/ . [7] J. Kuipers, T. Ueda and J. A. M. Vermaseren, “Code Optimization in FORM,” Comput. Phys. Commun. 189 (2015) 1 [arXiv:1310.7007], http://www.nikhef.nl/ form/ . [8] W. Bruns, B. Ichim, B. and T. Romer, C. Soger, “Normaliz. Algorithms for rational cones and affine monoids.” http://www.math.uos.de/normaliz/ .

Subfields: Atomic and Molecular Physics, and Optics, Computational Theory and Mathematics, Information Systems

View Paper: https://www.semanticscholar.org/paper/79cc64c75247dcbf135f8b09ec05b83a12dbf2d5

How researchers see this paper

The researchers whose trusted networks rank it highest, strongest advocate first. Every researcher gets a different order.

In Atomic and Molecular Physics, and Optics

  1. Thomas Gehrmann
  2. Gudrun Heinrich
  3. Stefan Kallweit
  4. G. Luisoni
  5. E. Glover
  6. Aude Gehrmann-De Ridder
  7. Stefano Pozzorini
  8. M. Schönherr

In Computational Theory and Mathematics

  1. Gudrun Heinrich
  2. Stefan Kallweit
  3. Thomas Gehrmann
  4. Stefano Pozzorini
  5. Ansgar Denner
  6. I. Moult
  7. G. Luisoni
  8. A. Larkoski

In Information Systems

  1. Gudrun Heinrich

Citations in the corpus