Hypergeometric Series Representations of Feynman Integrals by GKZ Hypergeometric Systems
René Pascal Klausen
TL;DR
The paper establishes a GKZ-hypergeometric framework for Feynman integrals by encoding each graph via the Lee-Pomeransky polynomial $G$ and its Newton polytope $oldsymbol{ riangle_G}$. For regular unimodular triangulations, generalized Feynman integrals $J_ ext{A}(oldsymbol{ u}, z)$ admit explicit Horn-type series representations as finite sums of GKZ $ ext{Gamma}$-series with boundary-determined coefficients, and these representations extend to the $oldsymbol{ u}$-expansion and various parametric forms. It further shows how non-generic limits and $oldsymbol{ u}$-expansions can be managed through limit procedures, leading to a practical algorithm for computing Feynman amplitudes in many kinematic regions, as demonstrated in the full massive sunset example. The study also discusses the role of unimodularity in triangulations, the treatment of non-unimodular cases via subtriangulations or dilations, and outlines open questions about extending unimodularity and applying GKZ to broader parametric representations, with implications for numerical efficiency and analytic structure of amplitudes.
Abstract
We show that almost all Feynman integrals as well as their coefficients in a Laurent series in dimensional regularization can be written in terms of Horn hypergeometric functions. By applying the results of Gelfand-Kapranov-Zelevinsky (GKZ) we derive a formula for a class of hypergeometric series representations of Feynman integrals, which can be obtained by triangulations of the Newton polytope $Δ_G$ corresponding to the Lee-Pomeransky polynomial $G$. Those series can be of higher dimension, but converge fast for convenient kinematics, which also allows numerical applications. Further, we discuss possible difficulties which can arise in a practical usage of this approach and give strategies to solve them.
