High-precision calculation of multi-loop Feynman integrals by difference equations
read the original abstract
We describe a new method of calculation of generic multi-loop master integrals based on the numerical solution of systems of difference equations in one variable. We show algorithms for the construction of the systems using integration-by-parts identities and methods of solutions by means of expansions in factorial series and Laplace's transformation. We also describe new algorithms for the identification of master integrals and the reduction of generic Feynman integrals to master integrals, and procedures for generating and solving systems of differential equations in masses and momenta for master integrals. We apply our method to the calculation of the master integrals of massive vacuum and self-energy diagrams up to three loops and of massive vertex and box diagrams up to two loops. Implementation in a computer program of our approach is described. Important features of the implementation are: the ability to deal with hundreds of master integrals and the ability to obtain very high precision results expanded at will in the number of dimensions.
This paper has not been read by Pith yet.
Forward citations
Cited by 35 Pith papers
-
Efficient AI-Inspired Reduction of Feynman Integrals via Tube Seeding
Machine learning discovers a tube-seeding strategy for IBP reduction of Feynman integrals that scales linearly with numerator power, demonstrated on rank-20 2-loop 5-point integrals.
-
The four-loop non-singlet splitting functions in QCD
Four-loop non-singlet splitting functions in QCD are computed analytically for the first time, with numerical representations provided.
-
Learning to Unscramble Feynman Loop Integrals with SAILIR
A self-supervised transformer learns to unscramble Feynman integrals for online IBP reduction, delivering bounded memory use on complex two-loop topologies while matching Kira's speed on the hardest cases tested.
-
The photon-energy spectrum in $B\to X_s\gamma$ to N$^3$LO: light-fermion and large-$N_{\rm c}$ corrections
N3LO calculation of the B to Xs gamma photon spectrum including complete light-fermion corrections, two massive fermion loops, and large-Nc terms, with improved results in kinetic and MSR mass schemes.
-
Emergence of Calabi-Yau manifolds in high-precision black hole scattering
At 5PM-1SF order, Calabi-Yau three-fold periods emerge in radiation-reacted observables for classical black hole scattering computed with worldline QFT and advanced IBP/DE methods.
-
Gravitational Compton scattering at the fourth post-Minkowskian order
Derives gravitational Compton amplitude at O(G^4) and N-matrix element for scattering phase shift, verified by agreement with black-hole perturbation theory.
-
Taming Symbolic IBP Reduction with Intermediate Bases
An algorithm reconstructs symbolic IBP reduction coefficients via intermediate bases, demonstrated on massive box-triangle and pentagon-triangle integrals using 3289 and 13013 samplings versus over a million unknowns.
-
Analytic results for heavy-quark contributions to charged-current DIS at NNLO
Analytic NNLO partonic coefficient functions for F2, FL, F3 in charged-current DIS with exact charm mass, expressed via Goncharov polylogarithms and Chen iterated integrals.
-
NNLO QCD predictions for $t\bar t W$ production at hadron colliders
NNLO QCD predictions for ttW production at hadron colliders using direct two-loop amplitude computation in the generalised leading-colour limit.
-
Magic Relations and Critical Varieties of Feynman Integrals
Magic relations in Feynman integral families coincide with higher-dimensional critical varieties, enabling a practical test to detect and handle them.
-
Heavy-Quark Condensate and Vacuum Energy Anomalous Dimension at Five Loops
The heavy-quark condensate is computed at five-loop order in QCD with massive quarks, confirming the five-loop vacuum anomalous dimension.
-
SIRENA -- Sum-Integral REductioN Algorithm
SIRENA automates IBP reduction of sum-integrals in finite-temperature QFT, reproduces known results to 3 loops, supplies new 3-loop fermionic reductions, and derives an analytic factorization formula for arbitrary 2-l...
-
Tensor decomposition of $e^+e^-\to\pi^+\pi^-\gamma$ to higher orders in the dimensional regulator
First beyond-NLO tensor decomposition and higher-order analytic one-loop amplitudes for e+e- to pi+pi-gamma, paired with a fast numerical five-point integral evaluator.
-
The OPE Approach to Renormalization: Operator Mixing
OPE-based recursive renormalization for mixed composite operators gives five-loop anomalous dimensions in phi^4 and two-loop in phi^3 models.
-
Feynman integral reduction by covariant differentiation
Covariant differentiation on the dual vector space spanned by master integrals reduces a large class of Feynman integrals to masters, with connections reusable across mass configurations.
-
Feynman integral reduction with intersection theory made simple
Branch representation reduces the variable count for intersection-theory-based Feynman integral reduction to at most 3L-3 for L-loop integrals regardless of leg number.
-
Learning to Unscramble: Simplifying Symbolic Expressions via Self-Supervised Oracle Trajectories
A permutation-equivariant transformer trained on self-supervised oracle trajectories from scrambled expressions achieves near-perfect simplification rates for dilogarithms and 100% success on 5-point gluon scattering ...
-
Resumming Scattering Amplitudes for Waveforms
A new projector-based formalism determines effective potentials from perturbative amplitudes and resums them to compute non-perturbative gravitational waveforms for generic two-body trajectories.
-
Twisted Feynman Integrals: from generating functions to spin-resummed post-Minkowskian dynamics
Twisted Feynman integrals are introduced with graded Symanzik polynomials, classified as exponential periods, and shown to have geometry not inferable from generalized Baikov leading singularities.
-
Double virtual QCD corrections to $t\bar{t}+$jet production at the LHC
Leading-colour two-loop virtual amplitudes for ttbar+jet are extracted analytically via finite-field evaluations and differential equations, then packaged in a C++ library with new numerical integration techniques.
-
Integral Reduction with Kira 2.0 and Finite Field Methods
Kira 2.0 implements finite-field coefficient reconstruction for IBP reductions and improved user-equation handling, yielding lower memory use and faster performance on state-of-the-art problems.
-
A numerical evaluation of planar two-loop helicity amplitudes for a W-boson plus four partons
First numerical evaluation of planar two-loop helicity amplitudes for W-boson plus four partons using finite-field reduction and sector decomposition on a subset of master integrals.
-
A Runway to Dissipation of Angular Momentum via Worldline Quantum Field Theory
The authors introduce static correlators in worldline QFT to compute angular momentum dissipation in black hole scattering, reproducing the known O(G^3) flux and extending the approach to electromagnetism at O(α^3).
-
Planar master integrals for two-loop NLO electroweak light-fermion contributions to $g g \rightarrow Z H$
Analytic expressions for the planar master integrals in two-loop NLO EW light-fermion contributions to gg → ZH are derived via canonical differential equations and expressed using Goncharov polylogarithms or one-fold ...
-
Four-loop Anomalous Dimensions of Scalar-QED Theory from Operator Product Expansion
Four-loop anomalous dimension of φ^Q in scalar-QED computed via OPE, extending prior three-loop results and validating the method in a gauge theory.
-
New algorithms for Feynman integral reduction and $\varepsilon$-factorised differential equations
A geometric order relation in IBP reduction yields a master-integral basis with Laurent-polynomial differential equations on the maximal cut that are then ε-factorized.
-
CHESS: CHEbyshev pSeudo-Spectral transport for Feynman integral differential equations
CHESS package implements Chebyshev-Lobatto spectral collocation for transporting epsilon-factorized differential equations of Feynman master integrals with benchmarks showing rapid convergence and shorter wall times t...
-
Tracing Transcendentality in Protected Correlators of N=4 SYM
Explicit two-loop computations of protected correlators in N=4 SYM yield a universal one-loop term and a planar extrapolation at arbitrary dimension controlled by stress-tensor multiplet count.
-
Weak-field waveforms for generic relativistic orbits
Outlines a Schwinger-Keldysh path-integral framework that derives worldline equations of motion and computes weak-field gravitational waveforms independently for unspecified relativistic orbits.
-
Two-loop QCD corrections to $ H \rightarrow b + \bar{b} + g $ at higher powers in the dimensional regulator
Two-loop QCD corrections to H → b b-bar g amplitude at higher powers of ε in dimensional regularization.
-
Numerical analytical continuation of multivariate hypergeometric functions
A general numerical framework is described for high-precision evaluation and analytic continuation of multivariate hypergeometric functions via Pfaffian systems and the Frobenius method.
-
Three loop QCD corrections to electroweak radiative parameters
Three-loop QCD corrections to electroweak radiative parameters Δρ, Δr, and Δκ are computed, yielding an updated W boson mass prediction relevant for FCC precision targets.
-
One-loop amplitudes for $t\bar{t}j$ and $t\bar{t}\gamma$ productions at the LHC through $\mathcal{O}(\epsilon^2)$
Analytic expressions for one-loop helicity amplitudes in ttj and ttγ production are derived to O(ε²) as linear combinations of pentagon functions with rational coefficients in momentum-twistor variables, obtained via ...
-
Graphical Functions by Examples
Graphical functions, defined as massless three-point position-space integrals, serve as a powerful tool for evaluating multi-loop Feynman integrals, with extensions to conformal field theory and recent algorithmic com...
-
Les Houches 2023 -- Physics at TeV Colliders: Report on the Standard Model Precision Wishlist
The report reviews progress since 2021 in fixed-order computations for LHC applications and identifies processes requiring missing higher-order corrections to match anticipated experimental precision.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.