- Physical Review Journals
All JournalsPhysics Magazine

Physical Review D

Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Multiple Mellin-Barnes integrals with polygamma functions

Sumit Banik1,2,* and Samuel Friot3,4,†

  • *Contact author: sumit.banik@physik.uzh.ch
  • †Contact author: samuel.friot@universite-paris-saclay.fr
PDF

Phys. Rev. D 114, 056004 – Published 2 September, 2026

DOI: https://doi.org/10.1103/yfvl-dfk2

Abstract

Mellin-Barnes (MB) integrals appear in various branches of physics and mathematics and are, in particular, used as a standard tool for evaluating multiloop, multiscale Feynman integrals both analytically and numerically. Recent geometric approaches based on conic hulls and triangulations provide a systematic framework for computing multiple MB integrals in terms of multivariate series. These approaches have so far been limited to MB integrals whose integrands are ratios of products of Euler’s gamma functions only. However, in Feynman integral calculus, MB integrals with polygamma functions naturally arise, for instance, after resolving singularities in the dimensional-regularization parameter ε and expanding the MB integrand in powers of ε, as done by the public codes mb.m and mbresolve.m. In this paper, we extend the conic hull and triangulation methods to the computation of MB integrals having polygamma functions in their integrand. We show that the arguments of polygamma functions can be treated in a similar way to the arguments of gamma functions when applying the conic hull and triangulation techniques to identify poles that would contribute to different series solutions. However, since the singularity structure of the polygamma function is different from that of the gamma function, we propose two different ways to compute MB integrals involving polygamma functions, depending on whether the MB integral has straight or nonstraight contours. We have implemented these algorithms in an updated version of the Mathematica package mbconichulls.wl, and we illustrate their use with a set of examples from Feynman integral calculus.

View figure in article

Physics Subject Headings (PhySH)

Article Text

References (79)

  1. V. A. Smirnov, Analytic Tools for Feynman Integrals (Springer, Berlin, Heidelberg, 2012), vol. 250.
  2. I. Dubovyk, J. Gluza, and G. Somogyi, Mellin-Barnes Integrals: A Primer on Particle Physics Applications, Lect. Notes Phys. (Springer, Cham, 2022), Vol. 1008.
  3. R. J. Sasiela, Electromagnetic Wave Propagation in Turbulence: Evaluation and Application of Mellin Transforms, Springer Series on Wave Phenomena (Springer, Berlin, Heidelberg, 1994), Vol. 18.
  4. S. Friot, On Ruby’s solid angle formula and some of its generalizations, Nucl. Instrum. Methods Phys. Res., Sect. A 773, 150 (2015).
  5. D. O. Oriekhov and V. P. Gusynin, RKKY interaction in a doped pseudospin-1 fermion system at finite temperature, Phys. Rev. B 101, 235162 (2020).
  6. P. Appell and J. Kampé de Fériet, Fonctions hypergéométriques et hypersphériques—Polynômes d’Hermite (Gauthier-Villars et Cie, Paris, 1926).
  7. H. Exton, Multiple Hypergeometric Functions and Applications, Ellis Horwood Series in Mathematics and its Applications (Ellis Horwood Ltd., Chichester, England, 1976).
  8. O. I. Marichev, Handbook of Integral Transforms of Higher Transcendental Functions: Theory and Algorithmic Tables, Ellis Horwood Series in Mathematics and its Applications (Ellis Horwood Ltd., Chichester, England, 1983).
  9. H. Srivastava and P. Karlsson, Multiple Gaussian Hypergeometric Series, Ellis Horwood Series in Mathematics and its Applications (E. Horwood, Chichester, England, 1985).
  10. R. B. Paris and D. Kaminski, Asymptotics and Mellin-Barnes Integrals, Encyclopedia of Mathematics and its Applications (Cambridge University Press, Cambridge, England, 2001), Vol. 85.
  11. I. Gonzalez, I. Kondrashuk, V. H. Moll, and L. M. Recabarren, Mellin–Barnes integrals and the method of brackets, Eur. Phys. J. C 82, 28 (2022).
  12. T.-F. Feng, C.-H. Chang, J.-B. Chen, and H.-B. Zhang, GKZ-hypergeometric systems for Feynman integrals, Nucl. Phys. B953, 114952 (2020).
  13. J.-P. Aguilar, C. Coste, and J. O. Korbel, Series representation of the pricing formula for the european option driven by space-time fractional diffusion, Fract. Calc. Appl. Anal. 21, 981 (2018).
  14. J.-P. Aguilar and J. O. Korbel, Option pricing models driven by the space-time fractional diffusion: Series representation and applications, Fractal Fract. 2, 15 (2018).
  15. J. Gluza, K. Kajda, and T. Riemann, AMBRE: A Mathematica package for the construction of Mellin-Barnes representations for Feynman integrals, Comput. Phys. Commun. 177, 879 (2007).
  16. A. V. Belitsky, A. V. Smirnov, and V. A. Smirnov, MB tools reloaded, Nucl. Phys. B986, 116067 (2023).
  17. V. A. Smirnov, Analytical result for dimensionally regularized massless on-shell double box, Phys. Lett. B 460, 397 (1999).
  18. J. B. Tausk, Non-planar massless two-loop Feynman diagrams with four on-shell legs, Phys. Lett. B 469, 225 (1999).
  19. M. Czakon, Automatized analytic continuation of Mellin-Barnes integrals, Comput. Phys. Commun. 175, 559 (2006).
  20. A. V. Smirnov and V. A. Smirnov, On the resolution of singularities of multiple Mellin-Barnes integrals, Eur. Phys. J. C 62, 445 (2009).
  21. M. Kalmykov, V. Bytev, B. A. Kniehl, S.-O. Moch, B. F. L. Ward, and S. A. Yost, Hypergeometric Functions and Feynman Diagrams, in Antidifferentiation and the Calculation of Feynman Amplitudes (Springer, Cham, 2020).
  22. J. Vollinga and S. Weinzierl, Numerical evaluation of multiple polylogarithms, Comput. Phys. Commun. 167, 177 (2005).
  23. I. Dubovyk, J. Gluza, and T. Riemann, Optimizing the Mellin-Barnes approach to numerical multiloop calculations, Acta Phys. Pol. B 50, 1993 (2019).
  24. M. Y. Kalmykov and B. A. Kniehl, Counting the number of master integrals for sunrise diagrams via the Mellin-Barnes representation, J. High Energy Phys. 07 (2017) 031.
  25. M. Y. Kalmykov and B. A. Kniehl, Mellin-Barnes representations of Feynman diagrams, linear systems of differential equations, and polynomial solutions, Phys. Lett. B 714, 103 (2012).
  26. N. I. Usyukina, On a representation for three point function, Teor. Mat. Fiz. 22, 300 (1975).
  27. E. E. Boos and A. I. Davydychev, A Method of evaluating massive Feynman integrals, Theor. Math. Phys. 89, 1052 (1991).
  28. A. I. Davydychev, Some exact results for N point massive Feynman integrals, J. Math. Phys. (N.Y.) 32, 1052 (1991).
  29. A. I. Davydychev, General results for massive N point Feynman diagrams with different masses, J. Math. Phys. (N.Y.) 33, 358 (1992).
  30. N. I. Usyukina and A. I. Davydychev, An Approach to the evaluation of three and four point ladder diagrams, Phys. Lett. B 298, 363 (1993).
  31. A. I. Davydychev and J. B. Tausk, Two loop selfenergy diagrams with different masses and the momentum expansion, Nucl. Phys. B397, 123 (1993).
  32. C. Anastasiou, C. Duhr, F. Dulat, F. Herzog, and B. Mistlberger, Higgs Boson Gluon-Fusion production in QCD at three loops, Phys. Rev. Lett. 114, 212001 (2015).
  33. C. Anastasiou, C. Duhr, F. Dulat, E. Furlan, T. Gehrmann, F. Herzog, A. Lazopoulos, and B. Mistlberger, High precision determination of the gluon fusion Higgs boson cross-section at the LHC, J. High Energy Phys. 05 (2016) 058.
  34. C. Anastasiou, C. Duhr, F. Dulat, E. Furlan, T. Gehrmann, F. Herzog, and B. Mistlberger, Higgs Boson gluon-fusion production beyond threshold in N3LO QCD, J. High Energy Phys. 03 (2015) 091.
  35. I. Dubovyk, A. Freitas, J. Gluza, T. Riemann, and J. Usovitsch, Complete electroweak two-loop corrections to Z boson production and decay, Phys. Lett. B 783, 86 (2018).
  36. P. Banerjee, T. Engel, N. Schalch, A. Signer, and Y. Ulrich, Bhabha scattering at NNLO with next-to-soft stabilisation, Phys. Lett. B 820, 136547 (2021).
  37. J. Davies, G. Mishima, K. Schönwald, M. Steinhauser, and H. Zhang, Higgs boson contribution to the leading two-loop Yukawa corrections to gg→HH, J. High Energy Phys. 08 (2022) 259.
  38. J. Davies, K. Schönwald, M. Steinhauser, and H. Zhang, Next-to-leading order electroweak corrections to gg→HH and gg→gH in the large-mt limit, J. High Energy Phys. 10 (2023) 033.
  39. H. Zhang, Massive two-loop four-point Feynman integrals at high energies with AsyInt, J. High Energy Phys. 09 (2024) 069.
  40. J. Davies, K. Schönwald, M. Steinhauser, and H. Zhang, Analytic next-to-leading order Yukawa and Higgs boson self-coupling corrections to gg→ HH at high energies, J. High Energy Phys. 04 (2025) 193.
  41. J.-P. Aguilar, D. Greynat, and E. De Rafael, Muon anomaly from lepton vacuum polarization and the Mellin-Barnes representation, Phys. Rev. D 77, 093010 (2008).
  42. D. Greynat and E. de Rafael, Hadronic contributions to the muon anomaly in the constituent chiral quark model, J. High Energy Phys. 07 (2012) 020.
  43. J. Charles, E. de Rafael, and D. Greynat, Mellin-Barnes approach to hadronic vacuum polarization and gμ−2, Phys. Rev. D 97, 076014 (2018).
  44. B. Ananthanarayan, J. Bijnens, S. Friot, and S. Ghosh, Analytic representation of FK/Fπ in two loop chiral perturbation theory, Phys. Rev. D 97, 091502 (2018).
  45. Z. Bern, L. J. Dixon, and V. A. Smirnov, Iteration of planar amplitudes in maximally supersymmetric Yang-Mills theory at three loops and beyond, Phys. Rev. D 72, 085001 (2005).
  46. Z. Bern, M. Czakon, L. J. Dixon, D. A. Kosower, and V. A. Smirnov, The four-loop planar amplitude and cusp anomalous dimension in maximally supersymmetric Yang-Mills theory, Phys. Rev. D 75, 085010 (2007).
  47. G. Somogyi, Angular integrals in d dimensions, J. Math. Phys. (N.Y.) 52, 083501 (2011).
  48. T. Ahmed, S. M. Hasan, and A. Rapakoulias, Phase-space integrals through Mellin-Barnes representation, Phys. Rev. D 112, L051903 (2025).
  49. T. Ahmed, S. M. Hasan, and A. Rapakoulias, Angular phase-space integrals with four denominators through Mellin-Barnes, J. High Energy Phys. 11 (2025) 152.
  50. S. Banik, A. Crivellin, L. Naterop, and P. Stoffer, Two-loop anomalous dimensions for baryon-number-violating operators in SMEFT, J. High Energy Phys. 02 (2026) 017.
  51. S. Weinzierl, Feynman Integrals, A Comprehensive Treatment for Students and Researchers, UNITEXT for Physics (Springer, New York, 2022).
  52. A. V. Kotikov, Differential equations method: New technique for massive Feynman diagrams calculation, Phys. Lett. B 254, 158 (1991).
  53. T. Gehrmann and E. Remiddi, Differential equations for two-loop four-point functions, Nucl. Phys. B580, 485 (2000).
  54. J. M. Henn, Multiloop integrals in dimensional regularization made simple, Phys. Rev. Lett. 110, 251601 (2013).
  55. M. Beneke and V. A. Smirnov, Asymptotic expansion of Feynman integrals near threshold, Nucl. Phys. B522, 321 (1998).
  56. A. Pak and A. Smirnov, Geometric approach to asymptotic expansion of Feynman integrals, Eur. Phys. J. C 71, 1626 (2011).
  57. B. Ananthanarayan, A. Pal, S. Ramanan, and R. Sarkar, Unveiling Regions in multi-scale Feynman integrals using singularities and power geometry, Eur. Phys. J. C 79, 57 (2019).
  58. V. A. Smirnov, Expansion by Regions: An Overview (Springer, Cham, 2021).
  59. J. Usovitsch, I. Dubovyk, and T. Riemann, MBnumerics: Numerical integration of Mellin-Barnes integrals in physical regions, Proc. Sci., LL2018 (2018) 046 [arXiv:1810.04580].
  60. B. Ananthanarayan, S. Banik, S. Friot, and S. Ghosh, Massive one-loop conformal Feynman integrals and quadratic transformations of multiple hypergeometric series, Phys. Rev. D 103, 096008 (2021).
  61. K. Alkalaev and S. Mandrygin, Multipoint conformal integrals in D dimensions. Part I. Bipartite Mellin-Barnes representation and reconstruction, J. High Energy Phys. 09 (2025) 118.
  62. K. Alkalaev and S. Mandrygin, Multipoint conformal integrals in D dimensions. Part II. Polygons and basis functions, J. High Energy Phys. 12 (2025) 088.
  63. F. Loebbert, D. Müller, and H. Münkler, Yangian bootstrap for conformal Feynman integrals, Phys. Rev. D 101, 066006 (2020).
  64. B. Ananthanarayan, S. Banik, S. Friot, and S. Ghosh, Double box and hexagon conformal Feynman integrals, Phys. Rev. D 102, 091901 (2020).
  65. S. Choudhury, Notes On de-Sitter Mellin Barnes amplitudes, arXiv:2512.09175.
  66. A. O. Barvinsky, A. E. Kalugin, and W. Wachowski, Functorial properties of Schwinger-DeWitt expansion and Mellin-Barnes representation, Phys. Rev. D 113, 045005 (2026).
  67. B. Ananthanarayan, S. Banik, S. Friot, and S. Ghosh, Multiple series representations of N-fold Mellin-Barnes integrals, Phys. Rev. Lett. 127, 151601 (2021).
  68. S. Banik and S. Friot, Multiple Mellin-Barnes integrals and triangulations of point configurations, Phys. Rev. D 110, 036002 (2024).
  69. S. Banik, MBConicHulls: A Mathematica package for analytic evaluation of N-fold MB integrals. https://github.com/SumitBanikGit/MBConicHulls/.
  70. C. Schneider, Simplifying Multiple Sums in Difference Fields (Springer, Vienna, 2013).
  71. S. Banik and S. Friot, Multiple Mellin-Barnes integrals with straight contours, Phys. Rev. D 107, 016007 (2023).
  72. S. Friot and D. Greynat, On convergent series representations of Mellin-Barnes integrals, J. Math. Phys. (N.Y.) 53, 023508 (2012).
  73. K. J. Larsen and R. Rietkerk, MultivariateResidues: A Mathematica package for computing multivariate residues, Comput. Phys. Commun. 222, 250 (2018).
  74. C. Anastasiou, C. Duhr, F. Dulat, and B. Mistlberger, Soft triple-real radiation for Higgs production at N3LO, J. High Energy Phys. 07 (2013) 003.
  75. B. Eden, The off-shell one- and two-loop box recovered from intersection theory, Nucl. Phys. B1024, 117336 (2026).
  76. B. Ananthanarayan, S. Banik, V. Bernard, S. Friot, S. Ghosh, and U.-G. Meißner, Sunset integrals with up to three mass scales in chiral perturbation theory: A comparative study of the Mellin-Barnes representation technique, Eur. Phys. J. ST 235, 1833 (2026).
  77. A. I. Davydychev and M. Y. Kalmykov, New results for the epsilon expansion of certain one, two and three loop Feynman diagrams, Nucl. Phys. B605, 266 (2001).
  78. S. Banik and S. Friot, Automated computation of multiple Mellin-Barnes integrals having polygamma functions in their integrand, Acta Phys. Pol. B Proc. Suppl. 19, 2-A8 (2026).
  79. J. Rambau, Topcom: Triangulations of point configurations and oriented matroids, Konrad-Zuse-Zentrum für Informationstechnik Berlin (ZIB) Technical Report No. ZIB-Report 02-17, 2002.

Outline

Information

Phys. Rev. D 114, 056004– Published 2 September, 2026
  • Received 5 January 2026
  • Accepted 22 June 2026
Crossmark - Check for updates button
Creative Commons Logo

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP3.

Published by the American Physical Society

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation