New algorithms for Feynman integral reduction and varepsilon-factorised differential equations
Pith reviewed 2026-05-17 20:56 UTC · model grok-4.3
The pith
A new geometric order on integration-by-parts identities produces master integrals with Laurent-polynomial differential equations in ε on the maximal cut.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using a new geometric order relation in the integration-by-parts reduction produces a basis of master integrals whose differential equations on the maximal cut are Laurent polynomials in the regularisation parameter ε and compatible with a filtration. This reduction step is performed entirely with rational functions. A subsequent procedure then ε-factorises the resulting differential equations.
What carries the argument
Geometric order relation applied to integration-by-parts identities, which selects a preferred basis making the maximal-cut differential equations Laurent polynomials in ε and compatible with a filtration.
If this is right
- The reduction to master integrals proceeds using only rational arithmetic.
- Maximal-cut differential equations take a simple Laurent polynomial form in ε.
- The selected basis respects a filtration that simplifies subsequent analysis.
- The two-step procedure applies to integral families of varying complexity.
Where Pith is reading between the lines
- The rational-function step could be automated inside existing computer-algebra systems to reduce manual basis choices.
- Extending the ordering criterion beyond the maximal cut might simplify the construction of full differential equations for scattering amplitudes.
- The approach may connect to other geometric selection methods used in algebraic-geometry treatments of Feynman integrals.
Load-bearing premise
The newly introduced geometric order relation on the integration-by-parts identities produces a basis whose maximal-cut differential equations are automatically Laurent polynomials in ε and compatible with the required filtration.
What would settle it
A concrete Feynman integral family for which the geometric ordering fails to produce differential equations on the maximal cut that are Laurent polynomials in ε or that violate the filtration condition.
Figures
read the original abstract
In this paper, we give a detailed account of the algorithm outlined in [1] for Feynman integral reduction and $\varepsilon$-factorised differential equations. The algorithm consists of two steps. In the first step, we use a new geometric order relation in the integration-by-parts reduction to obtain a basis of master integrals, whose differential equations on the maximal cut are of a Laurent polynomial form in the regularisation parameter $\varepsilon$ and compatible with a filtration. This step works entirely with rational functions. In a second step, we provide a method to $\varepsilon$-factorise the aforementioned Laurent differential equations. The second step may introduce algebraic and transcendental functions. We illustrate the versatility of the algorithm by applying it to different examples with a wide range of complexity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a two-step algorithm for Feynman integral reduction and ε-factorised differential equations. The first step introduces a geometric order relation on integration-by-parts identities to select a basis of master integrals such that the maximal-cut differential equations are Laurent polynomials in the regularisation parameter ε and compatible with a filtration; this step operates entirely with rational functions. The second step provides a method to ε-factorise the resulting Laurent equations, which may introduce algebraic or transcendental functions. The algorithm is illustrated on examples spanning a range of complexities.
Significance. If the geometric order relation reliably produces the claimed Laurent-polynomial and filtration-compatible form for arbitrary integral families, the work would offer a constructive, rational-first route to preparing master-integral bases for ε-factorised systems, which is useful for multi-loop phenomenology. The explicit separation of the rational reduction step from the subsequent factorisation step is a clear methodological strength, and the constructive nature of the procedure (starting from standard IBP identities without parameter fitting) is a positive feature.
major comments (2)
- [Abstract and §3] Abstract and §3 (geometric order construction): the claim that the newly defined geometric order on IBP identities automatically yields maximal-cut DEs that are Laurent polynomials in ε and filtration-compatible for arbitrary families is not supported by a general argument; the reduction matrix's preservation of the precise ε-pole structure is shown only on selected examples, leaving open the possibility that additional ε-dependent denominators in some sectors violate the Laurent form.
- [§4] §4 (examples): explicit verification that the output bases satisfy the Laurent-polynomial and filtration properties is not provided for all cited cases, which is required to substantiate that the first step works entirely with rational functions without sector-specific adjustments.
minor comments (2)
- [§2] Clarify the precise definition of the geometric order relation with a short worked example in the main text (rather than deferring all details to an appendix) to improve readability.
- [§5] Add a brief discussion of how the method scales with the number of propagators or loops, even if only qualitatively.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and the constructive comments. The distinction between the rational reduction step and the subsequent factorisation is indeed a key feature of the approach, and we welcome the opportunity to strengthen the presentation of the geometric order construction and the supporting verifications.
read point-by-point responses
-
Referee: [Abstract and §3] Abstract and §3 (geometric order construction): the claim that the newly defined geometric order on IBP identities automatically yields maximal-cut DEs that are Laurent polynomials in ε and filtration-compatible for arbitrary families is not supported by a general argument; the reduction matrix's preservation of the precise ε-pole structure is shown only on selected examples, leaving open the possibility that additional ε-dependent denominators in some sectors violate the Laurent form.
Authors: The geometric order is defined directly on the space of integrals using a filtration that is independent of ε. Because the IBP relations are generated over the field of rational functions in the kinematic variables and ε, and the ordering selects a basis in which no integral is reduced using a relation that would place ε in a denominator, the resulting reduction matrix remains in the Laurent polynomial ring. This property follows from the construction rather than from case-by-case verification. We will revise §3 to spell out this reasoning explicitly, including a short argument that any potential ε-dependent denominator would contradict the maximality of the chosen order. While the manuscript demonstrates the procedure on families of increasing complexity, we do not assert a universal theorem covering every conceivable integral family; the algorithmic guarantee is tied to the rational character of the reduction. revision: partial
-
Referee: [§4] §4 (examples): explicit verification that the output bases satisfy the Laurent-polynomial and filtration properties is not provided for all cited cases, which is required to substantiate that the first step works entirely with rational functions without sector-specific adjustments.
Authors: We agree that explicit checks improve clarity. In the revised version we will add, either in the main text or in an appendix, a concise verification for each example that the differential equations on the maximal cut are indeed Laurent polynomials in ε and respect the filtration induced by the geometric order. These checks will be performed directly on the reduction matrices obtained in the first step, confirming that only rational operations are used. revision: yes
Circularity Check
New geometric order relation in IBP reduction yields desired DE form via constructive rational procedure without self-referential reduction.
full rationale
The paper presents a two-step algorithm: first, a newly introduced geometric order relation on IBP identities produces a master-integral basis whose maximal-cut differential equations are Laurent polynomials in ε and filtration-compatible, all using rational functions; second, ε-factorisation may introduce algebraic/transcendental functions. This is described as a constructive procedure starting from standard IBP identities. No equations or steps reduce the target basis property to a fitted parameter, renamed input, or self-citation chain that is itself unverified. The reference to the outline in [1] is acknowledged but the present manuscript supplies the detailed account and examples; the central claim does not collapse to a prior self-citation by construction. The derivation is therefore self-contained against the stated inputs and does not exhibit the enumerated circularity patterns.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Integration-by-parts identities generate a finite basis of master integrals for any given Feynman integral family.
- ad hoc to paper A geometric order relation on the integrals can be defined such that the resulting differential equations on the maximal cut are Laurent polynomials in ε.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the pole order o ... number of non-zero consecutive residues r ... filtrations W•, F•_geom and F•_comb
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 3 Pith papers
-
The spectrum of Feynman-integral geometries at two loops
Two-loop Feynman integrals involve Riemann spheres, elliptic curves, hyperelliptic curves of genus 2 and 3, K3 surfaces, and a rationalizable Del Pezzo surface of degree 2.
-
An Algorithm for the Symbolic Reduction of Multi-loop Feynman Integrals via Generating Functions
A new generating-function framework turns IBP relations into differential equations in a non-commutative algebra, yielding an iterative algorithm that derives symbolic reduction rules and checks completeness for topol...
-
Integrand Analysis, Leading Singularities and Canonical Bases beyond Polylogarithms
Feynman integrals selected for unit leading singularities in complex geometries satisfy epsilon-factorized differential equations with new transcendental functions corresponding to periods and differential forms in th...
Reference graph
Works this paper leans on
-
[1]
Bree et al., (2025), arXiv:2506.09124
ε -collaboration, I. Bree et al., (2025), arXiv:2506.09124
-
[2]
E. D’Hoker, M. Hidding, and O. Schlotterer, Phys. Rev. Le tt. 133, 021602 (2024), arXiv:2308.05044. 77
-
[3]
R. Marzucca, A. J. McLeod, B. Page, S. Pögel, and S. Weinzi erl, Phys. Rev. D 109, L031901 (2024), arXiv:2307.11497
-
[4]
L. de la Cruz and P . V anhove, Lett. Math. Phys. 114, 89 (2024), arXiv:2401.09908
-
[5]
Schottky–Kronecker forms and hyperelliptic polylogarithms,
K. Baune, J. Broedel, E. Im, A. Lisitsyn, and F. Zerbini, J . Phys. A 57, 445202 (2024), arXiv:2406.10051
-
[6]
Higher-genus Fay-like identities from meromorphic generating functions,
K. Baune, J. Broedel, E. Im, A. Lisitsyn, and Y . Moeckli, ( 2024), arXiv:2409.08208
- [7]
- [8]
- [9]
- [10]
- [11]
- [12]
-
[13]
Relating flat connections and polylogarithms on higher genus Riemann surfaces,
E. D’Hoker, B. Enriquez, O. Schlotterer, and F. Zerbini , (2025), arXiv:2501.07640
- [14]
-
[15]
Feynman integrals, elliptic integrals and two-parameter K3 surfaces
C. Duhr and S. Maggio, (2025), arXiv:2502.15326
-
[16]
Duhr, Modular forms for three-loop banana integrals , 2502.15325
C. Duhr, (2025), arXiv:2502.15325
-
[17]
M. Becchetti, C. Dlapa, and S. Zoia, Phys. Rev. D 112, L031501 (2025), arXiv:2503.03603
-
[18]
Aspects of canonical differential equations for Calabi-Yau geometries and beyond
C. Duhr et al., JHEP 06, 128 (2025), arXiv:2503.20655
- [19]
-
[20]
F. V . Tkachov, Phys. Lett. B 100, 65 (1981)
work page 1981
-
[21]
K. G. Chetyrkin and F. V . Tkachov, Nucl. Phys. B 192, 159 (1981)
work page 1981
-
[22]
High-precision calculation of multi-loop Feynman integrals by difference equations
S. Laporta, Int. J. Mod. Phys. A 15, 5087 (2000), arXiv:hep-ph/0102033
work page internal anchor Pith review Pith/arXiv arXiv 2000
- [23]
-
[24]
A. V . Smirnov, Comput. Phys. Commun. 189, 182 (2015), arXiv:1408.2372
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[25]
Kira - A Feynman Integral Reduction Program
P . Maierhöfer, J. Usovitsch, and P . Uwer, Comput. Phys. Commun. 230, 99 (2018), arXiv:1705.05610
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [26]
- [27]
-
[28]
A. V . Kotikov, Phys. Lett. B 254, 158 (1991)
work page 1991
-
[29]
A. V . Kotikov, Phys. Lett. B 267, 123 (1991), [Erratum: Phys.Lett.B 295, 409–409 (1992)]
work page 1991
-
[30]
Remiddi,Differential equations for Feynman graph amplitudes,Nuovo Cim
E. Remiddi, Nuovo Cim. A 110, 1435 (1997), arXiv:hep-th/9711188
-
[31]
Differential Equations for Two-Loop Four-Point Functions
T. Gehrmann and E. Remiddi, Nucl. Phys. B 580, 485 (2000), arXiv:hep-ph/9912329
work page internal anchor Pith review Pith/arXiv arXiv 2000
-
[32]
J. M. Henn, Phys. Rev. Lett. 110, 251601 (2013), arXiv:1304.1806
work page internal anchor Pith review Pith/arXiv arXiv 2013
- [33]
-
[34]
X. Liu and Y .-Q. Ma, Comput. Phys. Commun. 283, 108565 (2023), arXiv:2201.11669
-
[35]
A Systematic and Efficient Method to Compute Multi-loop Master Integrals
X. Liu, Y .-Q. Ma, and C.-Y . Wang, Phys. Lett. B 779, 353 (2018), arXiv:1711.09572
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[36]
Z.-F. Liu and Y .-Q. Ma, Phys. Rev. Lett. 129, 222001 (2022), arXiv:2201.11637. 78
-
[37]
J. Henn et al., Phys. Rev. Lett. 135, 031601 (2025), arXiv:2501.01847
-
[38]
The $\varepsilon$-form of the differential equations for Feynman integrals in the elliptic case
L. Adams and S. Weinzierl, Phys. Lett. B781, 270 (2018), arXiv:1802.05020
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[39]
I. Hönemann, K. Tempest, and S. Weinzierl, Phys. Rev. D 98, 113008 (2018), arXiv:1811.09308, [Erratum: Phys.Rev.D 110, 059901 (2024)]
- [40]
- [41]
- [42]
- [43]
-
[44]
L. Görges, C. Nega, L. Tancredi, and F. J. Wagner, JHEP 07, 206 (2023), arXiv:2305.14090
- [45]
- [46]
- [47]
- [48]
- [49]
-
[50]
N. Schwanemann and S. Weinzierl, SciPost Phys. 18, 172 (2025), arXiv:2412.07522
- [51]
-
[52]
M. Becchetti, F. Coro, C. Nega, L. Tancredi, and F. J. Wag ner, JHEP 06, 033 (2025), arXiv:2502.00118
- [53]
-
[54]
J. Chen, L. L. Y ang, and Y . Zhang, (2025), arXiv:2503.23 720
work page 2025
- [55]
-
[56]
S. Pögel, X. Wang, and S. Weinzierl, JHEP 04, 117 (2023), arXiv:2212.08908
- [57]
-
[58]
S. Pögel, X. Wang, and S. Weinzierl, Phys. Rev. Lett. 130, 101601 (2023), arXiv:2211.04292
- [59]
- [60]
-
[61]
Calabi-Yau Feynman integrals in gravity:ε-factorized form for apparent singularities
H. Frellesvig, R. Morales, S. Pögel, S. Weinzierl, and M . Wilhelm, JHEP 02, 209 (2025), arXiv:2412.12057
-
[62]
On canonical differential equations for Calabi-Yau multi-scale Feynman integrals
S. Maggio and Y . Sohnle, JHEP 10, 202 (2025), arXiv:2504.17757
- [63]
-
[64]
S. Pögel, T. Teschke, X. Wang, and S. Weinzierl, (2025), arXiv:2507.23594
-
[65]
Moser, Mathematische Zeitschrift 1, 379 (1959)
J. Moser, Mathematische Zeitschrift 1, 379 (1959)
work page 1959
-
[66]
R. N. Lee, JHEP 04, 108 (2015), arXiv:1411.0911
work page internal anchor Pith review Pith/arXiv arXiv 2015
- [67]
-
[68]
epsilon: A tool to find a canonical basis of master integrals
M. Prausa, Comput. Phys. Commun. 219, 361 (2017), arXiv:1701.00725
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[69]
Fuchsia: a tool for reducing differential equations for Feynman master integrals to epsilon form
O. Gituliar and V . Magerya, Comput. Phys. Commun. 219, 329 (2017), arXiv:1701.04269
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[70]
Algorithmic transformation of multi-loop master integrals to a canonical basis with CANONICA
C. Meyer, Comput. Phys. Commun. 222, 295 (2018), arXiv:1705.06252. 79
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[71]
D. Chicherin and V . Sotnikov, JHEP 20, 167 (2020), arXiv:2009.07803
-
[72]
D. Chicherin, V . Sotnikov, and S. Zoia, JHEP 01, 096 (2022), arXiv:2110.10111
-
[73]
M. Hidding, Comput. Phys. Commun. 269, 108125 (2021), arXiv:2006.05510
-
[74]
T. Armadillo, R. Bonciani, S. Devoto, N. Rana, and A. Vic ini, Comput. Phys. Commun. 282, 108545 (2023), arXiv:2205.03345
- [75]
-
[76]
P . Petit Rosàs and W . J. Torres Bobadilla, JHEP 09, 210 (2025), arXiv:2507.12548
-
[77]
Cluster adjacency properties of scattering amplitudes
J. Drummond, J. Foster, and Ö. Gürdo ˘gan, Phys. Rev. Lett. 120, 161601 (2018), arXiv:1710.10953
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[78]
D. Chicherin, J. M. Henn, and G. Papathanasiou, Phys. Re v. Lett. 126, 091603 (2021), arXiv:2012.12285
-
[79]
A. Pokraka, M. Spradlin, A. V olovich, and H.-C. Weng, (2 025), arXiv:2506.11895
-
[80]
A. Kristensson, M. Wilhelm, and C. Zhang, Phys. Rev. Let t. 127, 251603 (2021), arXiv:2106.14902
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.