Pith. sign in

REVIEW 4 major objections 4 minor 3 cited by

A generating function with two auxiliary vectors reduces every tensor integral of the two-loop sunset to seven master integrals through explicit recurrences.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-04 15:40 UTC pith:GDEJPXLB

load-bearing objection First two-loop generating-function reduction is a real technical step, but the paper's completeness and linear-independence claims are unproven and the recurrences are not independently verified; worth refereeing, not trusting yet. the 4 major comments →

arxiv 2509.18730 v2 pith:GDEJPXLB submitted 2025-09-23 hep-th

Tensor Reduction of Sunset by Generating Function

Xinhe Chen , Bo Feng , Liang Zhang This is my paper
classification hep-th
keywords generating functionsunset diagramtensor reductiontwo-loop Feynman integralsPV reductionsyzygy equationsrecurrence relationsmaster integrals
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper attempts to establish that the generating-function reduction method, previously tested only at one loop, works for a genuine two-loop integral: the sunset diagram. The sunset is the simplest two-loop topology, with three propagators connecting two loop momenta and one external momentum, and its general tensor integrals are a natural first stress test. The authors insert exponentials of two auxiliary vectors into the sunset integrand, so that every tensor integral appears as a coefficient in the expansion of the resulting generating function. They construct five differential equations for these coefficients—three PV-type natural equations and two syzygy equations free of doubled propagators—and solve them by series expansion. The result is a complete set of recurrence relations that reduce any high-rank tensor coefficient to four boundary coefficients and the seven master integrals.

Core claim

The authors' claim is that the complete tensor reduction of the two-loop sunset is encoded in a finite system of differential equations for a generating function, and that this system can be solved by series expansion. Every tensor integral of the sunset is a coefficient α_{abnmk} in the expansion of the generating function in the five scalar combinations built from the auxiliary vectors R1, R2 and the external momentum K. The first three differential equations yield recurrences that lower the indices n and m and reduce any α_{abnmk} to coefficients α_{a'b'000}; the two syzygy-derived equations then provide recurrences that solve α_{ab000} in terms of α_{0b'000} and α_{1b'000}. With four bou

What carries the argument

The central object is the generating function I_gen = ∫ d^D l1 d^D l2 e^{l1·R1+l2·R2}/(D1 D2 D3), where D1, D2, D3 are the three sunset propagators and R1, R2 are auxiliary vectors. Its expansion coefficients in the five independent contractions K·R1, K·R2, K^2 R1^2 − (K·R1)^2, K^2 R2^2 − (K·R2)^2, and K^2 R1·R2 − (K·R1)(K·R2) are exactly the tensor integrals of the sunset. The paper combines three natural PV-type differential equations with two syzygy equations—IBP identities constructed so that doubled propagators cancel—to form a complete system of five differential equations. Expanding these equations in the five scalar variables turns them into algebraic recurrences that lower the indic

Load-bearing premise

The load-bearing premise is that the five differential equations form a complete system and that the five equations used in the final step to solve α_{0b000} and α_{1b000} are linearly independent; if either fails at some spacetime dimension or kinematic point, the recurrences do not close.

What would settle it

Evaluate the 5×5 coefficient matrix of Step Six at, say, D=4 and b=2 or b=3; if its determinant is identically zero over the masses and K^2, the recurrences (3.24)–(3.25) degenerate and the claimed completeness fails. An independent check would generate a rank-4 sunset coefficient from the recurrences and compare it with the coefficient obtained by direct solution of the integration-by-parts identities at a generic mass point.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Any rank-r tensor integral of the sunset diagram can be reduced to the seven master integrals by repeated application of the recurrences, without solving a large linear system for each rank.
  • The generating-function method is shown to pass from one-loop integrals to a nontrivial two-loop topology, making the same differential-equation-plus-syzygy construction a plausible route for other multiloop topologies.
  • The number of boundary conditions needed to solve the system equals the number of top-sector master integrals, giving a new way to determine the number of master integrals from the completeness of the differential system.
  • The reduction coefficients are recovered as explicit functions of spacetime dimension D, the masses, and K^2, uniformly in tensor rank, giving analytic control over the coefficients.
  • Subsector integrals, such as the tadpole and two-propagator integrals, are treated as predetermined inputs and propagate through the recurrences, so improving subsector reduction automatically improves the sunset reduction.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Editor's inference: the same 'natural equations plus syzygy equations' construction should extend to other two-loop topologies, with the number of syzygy equations set by the defect of the PV equations; the main obstruction is finding simple syzygies that avoid doubled propagators in each topology.
  • Editor's inference: denominators such as D+2k−5, (a−1)a(−3+a+D), and b(b−1) in the recurrences indicate special spacetime dimensions where the recurrences need separate treatment; checking these limits could expose resonances requiring a different basis choice or new identities among master integrals.
  • Editor's inference: since the recurrences are rational functions of D and the kinematical variables, they could be combined with modular or finite-field reconstruction to produce analytic reduction coefficients for higher ranks where direct symbolic expansion is slow.
  • Editor's inference: the boundary-counting observation suggests a general consistency criterion—the dimension of the solution space of the completed differential system should match the number of top-sector master integrals—which could serve as a check when applying the method to new diagrams.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The manuscript proposes a generating-function method for the tensor reduction of the two-loop sunset diagram. It constructs a system of five differential equations for the reduction coefficients of the generating function: three 'natural' Passarino–Veltman-type equations (Sec. 2.1) and two syzygy-based IBP equations (Sec. 2.2). After expanding the coefficients in the five Lorentz scalars, the paper derives recurrence relations (Secs. 3.2–3.3) that are claimed to reduce every sunset tensor integral to the seven master integrals of Eq. (2.3). The beta sub-sector coefficients are computed independently from one-loop tadpole generating functions in Appendix A. The central claim is that the resulting recurrences are complete and efficient for arbitrary rank.

Significance. If the recurrences are correct, the paper would be a substantive step toward extending generating-function reduction beyond one loop, with no fitted parameters and an explicit avoidance of doubled-propagator IBP relations. The derivation is self-contained and the beta sub-sectors are obtained from independent tadpole inputs, which is a strength. However, the paper's central claim rests on two algebraic facts that are asserted but not proved: the completeness of the five differential equations and the linear independence of the five equations used in Step Six. A numerical or symbolic cross-check against an established IBP reduction is also absent. The method is promising, but the current verification is not yet commensurate with the strength of the claims.

major comments (4)
  1. [Sec. 3.3, Step Six] The recurrences (3.24) and (3.25) for alpha_{0b000} and alpha_{1b000} are the keystone of the whole reduction, but they follow from the unproved assertion that the five equations built from (3.21) and (3.22) 'can be considered linearly independent.' No determinant is given, and the coefficient matrix is not displayed. If this 5x5 matrix is singular for generic D or for some b (e.g., at the physical D=4 for a particular b), the recurrences do not follow and the claimed completeness fails. Please supply the determinant or an independent proof of nonsingularity, and state the exceptional values explicitly.
  2. [Sec. 2.2, end] The sentence 'Having obtained the complete set of differential equations' asserts completeness of the five-equation system (three PV equations plus two syzygy equations) without proof. This completeness is load-bearing: without it, the subsequent series solution cannot be guaranteed to determine all reduction coefficients. A proof that the five differential operators span the relevant solution space, or a check that the system has the correct number of independent solutions (matching the four initial conditions of Sec. 3.1), is needed.
  3. [Secs. 3.2–3.3, denominators] Several recurrence denominators vanish for values of the indices or dimension that are not excluded. Examples: (3.17)/(3.18) contain D+2k+2n-3 and D+2k+2m-3; (3.20) contains D+2k-5 and 12+D^2-11k+2k^2+D(3k-7); (3.23) contains (a-1)a(-3+a+D); (3.24) contains (b-1)b; (3.25) contains (b-1)b(-6+2b+3D). The paper only notes invalidity for b=0,1 in (3.24)-(3.25), but not for the other factors, some of which vanish at integer D (e.g., -3+a+D at D=2 for a=1). The status of these singular cases must be discussed; otherwise the recurrences cannot be claimed for all D and ranks.
  4. [General verification] The manuscript contains no independent check of the final recurrence system against a standard IBP reduction (e.g., LiteRed, FIRE, or Kira) for even one nontrivial tensor integral. Given the length and complexity of Eqs. (3.16)–(3.25), such a check—for instance, a rank-3 or rank-4 sunset tensor integral at D=4 or D=4-2eps—would substantially raise confidence. The absence of this check, together with the two unproved algebraic assertions above, leaves the central claim under-verified.
minor comments (4)
  1. [Eq. (2.3)] There is a typographical error: after the definition of I6 there is a double comma ', ,' before I7.
  2. [Sec. 3.2, Step Four] The text says four equations are used to solve for alpha_{ab11k}, but the equations themselves are not listed; the reader must reconstruct them from (3.13) and (3.14). A brief display or supplementary material would improve reproducibility.
  3. [App. A, Eq. (A.2)] The double factorial notation (r-1)!! is used without definition; consider adding a short note, especially since r can be odd/even and the expression implicitly vanishes for odd r.
  4. [Throughout] The symbol s0 is introduced only in Eq. (3.10); earlier references to K^2 would be clearer if s0 were defined at first use in Sec. 2.

Circularity Check

0 steps flagged

No circular reduction: recurrences are derived from IBP/syzygy differential equations with independent one-loop inputs; the asserted linear independence in Step Six is an unproved algebraic gap, not a circular step.

full rationale

The derivation chain is not circular. The reduction coefficients alpha are defined by Igen = I·alpha (Eq. 2.4) and solved from five differential equations: three PV equations (2.7)-(2.9) and two syzygy equations (2.17),(2.19), all derived from explicit IBP identities with the generating function. The source terms beta_{ij} are subsector reductions computed in Appendix A from one-loop tadpole generating functions; these are independent, externally checkable inputs, not fits and not the sunset tensor reduction being derived. No parameter is fitted to the target data and no 'prediction' is renamed input: the recurrences of Sec. 3 express higher alpha coefficients in terms of lower ones and terminate at the four initial conditions (3.4)-(3.7). The self-citations (e.g., [38,42,46]) supply standard one-loop results and the generating-function/PV framework, but the sunset-specific syzygy computation and recurrence derivation are carried out in this paper, so the self-citations are not load-bearing in the circular sense. The two structural concerns identified by a skeptical reader are real but are not circularity: (i) at the end of Sec. 2.2 the completeness of the five differential equations is asserted, not proved; (ii) in Sec. 3.3 Step Six the paper states 'these 5 equations can be considered linearly independent' and in Sec. 3.2 Step Four similar 4x4 solvability is used without determinant checks, so (3.24),(3.25) and the claimed completeness are conditional on these algebraic facts. Since neither step assumes the reduction it derives, these are verification/correctness gaps, not self-referential reductions. Score 1 reflects minor non-load-bearing self-citations and the conditional-but-not-circular completeness gap; no circular step is exhibited.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

The paper does not fit any parameters or introduce new physics. It relies on standard IBP identities and the completeness of the master basis. The main additional assumptions are the completeness of the syzygy-augmented differential system and the linear independence of the final five equations, both asserted without proof.

axioms (5)
  • domain assumption Integration-by-parts identities hold in dimensional regularization (Eq. 2.10).
    The entire derivation is based on IBP identities applied to the sunset integral.
  • domain assumption The generating function admits a convergent (formal) power series in the five Lorentz scalars y1, y2, x11, x22, x12 (Eq. 3.1).
    The series expansion and recurrence solution assume this expansion captures all tensor structures.
  • domain assumption The seven integrals in Eq. (2.3) form a complete basis for the sunset topology.
    The reduction to seven master integrals is standard for the general sunset.
  • ad hoc to paper The five differential equations (three PV, two syzygy) form a complete system for the reduction coefficients.
    Stated at the end of Sec. 2.2 without proof; the solving procedure is intended to demonstrate it.
  • ad hoc to paper The five equations selected in Step Six (Sec. 3.3) are linearly independent.
    Asserted: 'these 5 equations can be considered linearly independent'; no determinant or proof is given.

pith-pipeline@v1.3.0-alltime-deepseek · 23657 in / 15897 out tokens · 111393 ms · 2026-08-04T15:40:50.841198+00:00 · methodology

0 comments
read the original abstract

Recently, the generating function has been proposed as an alternative reduction method. This method has been tested at the one-loop level, including the tensor reduction and propagators with higher powers. In this work, we initiate the study of the method for higher loops by focusing on the sunset diagram, which is the simplest nontrivial two-loop integral. By employing PV reduction equations together with syzygy equations, we construct a complete system of differential equations. Through series expansion, we derive a complete set of recurrence relations, which can efficiently reduce any high-rank tensor structure.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Efficient AI-Inspired Reduction of Feynman Integrals via Tube Seeding

    hep-ph 2026-06 unverdicted novelty 8.0

    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.

  2. An Algorithm for the Symbolic Reduction of Multi-loop Feynman Integrals via Generating Functions

    hep-ph 2026-05 unverdicted novelty 7.0

    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...

  3. Twisted Feynman Integrals: from generating functions to spin-resummed post-Minkowskian dynamics

    hep-th 2025-12 unverdicted novelty 7.0

    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.

Reference graph

Works this paper leans on

55 extracted references · 49 linked inside Pith · cited by 3 Pith papers

  1. [1]

    V. A. Smirnov,Analytic tools for Feynman integrals, vol. 250. 2012, 10.1007/978-3-642-34886-0

    doi:10.1007/978-3-642-34886-0 2012

  2. [2]

    Weinzierl,Feynman Integrals

    S. Weinzierl,Feynman Integrals. A Comprehensive Treatment for Students and Researchers, UNITEXT for Physics. Springer, 2022, 10.1007/978-3-030-99558-4, [2201.03593]

    Pith/arXiv arXiv 2022

  3. [3]

    F. V. Tkachov,A theorem on analytical calculability of 4-loop renormalization group functions,Phys. Lett. B100(1981) 65

    1981

  4. [4]

    K. G. Chetyrkin and F. V. Tkachov,Integration by parts: The algorithm to calculate β-functions in 4 loops,Nucl. Phys. B192(1981) 159

    1981

  5. [5]

    Passarino and M

    G. Passarino and M. J. G. Veltman,One Loop Corrections for e+ e- Annihilation Into mu+ mu- in the Weinberg Model,Nucl. Phys. B160(1979) 151. – 18 –

    1979

  6. [6]

    Ossola, C

    G. Ossola, C. G. Papadopoulos and R. Pittau,Reducing full one-loop amplitudes to scalar integrals at the integrand level,Nucl. Phys. B763(2007) 147 [hep-ph/0609007]

    Pith/arXiv arXiv 2007

  7. [7]

    Ossola, C

    G. Ossola, C. G. Papadopoulos and R. Pittau,CutTools: A Program implementing the OPP reduction method to compute one-loop amplitudes,JHEP03(2008) 042 [0711.3596]

    Pith/arXiv arXiv 2008

  8. [8]

    Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower,One loop n point gauge theory amplitudes, unitarity and collinear limits,Nucl. Phys. B425(1994) 217 [hep-ph/9403226]

    Pith/arXiv arXiv 1994

  9. [9]

    Britto, F

    R. Britto, F. Cachazo and B. Feng,Generalized unitarity and one-loop amplitudes in N=4 super-Yang-Mills,Nucl. Phys. B725(2005) 275 [hep-th/0412103]

    Pith/arXiv arXiv 2005

  10. [10]

    Anastasiou, R

    C. Anastasiou, R. Britto, B. Feng, Z. Kunszt and P. Mastrolia,D-dimensional unitarity cut method,Phys. Lett. B645(2007) 213 [hep-ph/0609191]

    Pith/arXiv arXiv 2007

  11. [11]

    Britto and B

    R. Britto and B. Feng,Unitarity cuts with massive propagators and algebraic expressions for coefficients,Phys. Rev. D75(2007) 105006 [hep-ph/0612089]

    Pith/arXiv arXiv 2007

  12. [12]

    Liu, Y.-Q

    X. Liu, Y.-Q. Ma and C.-Y. Wang,A Systematic and Efficient Method to Compute Multi-loop Master Integrals,Phys. Lett. B779(2018) 353 [1711.09572]

    Pith/arXiv arXiv 2018

  13. [13]

    Liu and Y.-Q

    Z.-F. Liu and Y.-Q. Ma,Automatic computation of Feynman integrals containing linear propagators via auxiliary mass flow,Phys. Rev. D105(2022) 074003 [2201.11636]

    Pith/arXiv arXiv 2022

  14. [14]

    Badger, H

    S. Badger, H. Frellesvig and Y. Zhang,Hepta-Cuts of Two-Loop Scattering Amplitudes,JHEP04(2012) 055 [1202.2019]

    Pith/arXiv arXiv 2012

  15. [15]

    Zhang,Integrand-Level Reduction of Loop Amplitudes by Computational Algebraic Geometry Methods,JHEP09(2012) 042 [1205.5707]

    Y. Zhang,Integrand-Level Reduction of Loop Amplitudes by Computational Algebraic Geometry Methods,JHEP09(2012) 042 [1205.5707]

    Pith/arXiv arXiv 2012

  16. [16]

    K. J. Larsen and Y. Zhang,Integration-by-parts reductions from unitarity cuts and algebraic geometry,Phys. Rev. D93(2016) 041701 [1511.01071]

    Pith/arXiv arXiv 2016

  17. [17]

    Laporta,High-precision calculation of multiloop Feynman integrals by difference equations,Int

    S. Laporta,High-precision calculation of multiloop Feynman integrals by difference equations,Int. J. Mod. Phys. A15(2000) 5087 [hep-ph/0102033]

    Pith/arXiv arXiv 2000

  18. [18]

    Anastasiou and A

    C. Anastasiou and A. Lazopoulos,Automatic integral reduction for higher order perturbative calculations,JHEP07(2004) 046 [hep-ph/0404258]

    Pith/arXiv arXiv 2004

  19. [19]

    R. N. Lee,LiteRed 1.4: a powerful tool for reduction of multiloop integrals,J. Phys. Conf. Ser.523(2014) 012059 [1310.1145]

    Pith/arXiv arXiv 2014

  20. [20]

    A. V. Smirnov and F. S. Chukharev,FIRE6: Feynman Integral REduction with modular arithmetic,Comput. Phys. Commun.247(2020) 106877 [1901.07808]

    Pith/arXiv arXiv 2020

  21. [21]

    A. V. Smirnov and M. Zeng,FIRE 6.5: Feynman integral reduction with new simplification library,Comput. Phys. Commun.302(2024) 109261 [2311.02370]

    Pith/arXiv arXiv 2024

  22. [22]

    von Manteuffel and C

    A. von Manteuffel and C. Studerus,Reduze 2 - Distributed Feynman Integral Reduction,1201.4330. – 19 –

    Pith/arXiv arXiv

  23. [23]

    Maierh¨ ofer, J

    P. Maierh¨ ofer, J. Usovitsch and P. Uwer,Kira—A Feynman integral reduction program,Comput. Phys. Commun.230(2018) 99 [1705.05610]

    Pith/arXiv arXiv 2018

  24. [24]

    Klappert, F

    J. Klappert, F. Lange, P. Maierh¨ ofer and J. Usovitsch,Integral reduction with Kira 2.0 and finite field methods,Comput. Phys. Commun.266(2021) 108024 [2008.06494]

    Pith/arXiv arXiv 2021

  25. [25]

    Ruijl, T

    B. Ruijl, T. Ueda and J. A. M. Vermaseren,Forcer, a FORM program for the parametric reduction of four-loop massless propagator diagrams,Comput. Phys. Commun.253(2020) 107198 [1704.06650]

    Pith/arXiv arXiv 2020

  26. [26]

    Peraro,FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs,JHEP07(2019) 031 [1905.08019]

    T. Peraro,FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs,JHEP07(2019) 031 [1905.08019]

    Pith/arXiv arXiv 2019

  27. [27]

    Z. Wu, J. Boehm, R. Ma, H. Xu and Y. Zhang,NeatIBP 1.0, a package generating small-size integration-by-parts relations for Feynman integrals,Comput. Phys. Commun.295(2024) 108999 [2305.08783]

    Pith/arXiv arXiv 2024

  28. [28]

    X. Guan, X. Liu, Y.-Q. Ma and W.-H. Wu,Blade: A package for block-triangular form improved Feynman integrals decomposition,Comput. Phys. Commun.310 (2025) 109538 [2405.14621]

    Pith/arXiv arXiv 2025

  29. [29]

    Chen,Semi-automatic calculations of multi-loop Feynman amplitudes with AmpRed,Comput

    W. Chen,Semi-automatic calculations of multi-loop Feynman amplitudes with AmpRed,Comput. Phys. Commun.312(2025) 109607 [2408.06426]

    Pith/arXiv arXiv 2025

  30. [30]

    von Manteuffel and R

    A. von Manteuffel and R. M. Schabinger,A novel approach to integration by parts reduction,Phys. Lett. B744(2015) 101 [1406.4513]

    Pith/arXiv arXiv 2015

  31. [31]

    Peraro,Scattering amplitudes over finite fields and multivariate functional reconstruction,JHEP12(2016) 030 [1608.01902]

    T. Peraro,Scattering amplitudes over finite fields and multivariate functional reconstruction,JHEP12(2016) 030 [1608.01902]

    Pith/arXiv arXiv 2016

  32. [32]

    Liu and Y.-Q

    X. Liu and Y.-Q. Ma,Determining arbitrary Feynman integrals by vacuum integrals, Phys. Rev. D99(2019) 071501 [1801.10523]

    Pith/arXiv arXiv 2019

  33. [33]

    X. Guan, X. Liu and Y.-Q. Ma,Complete reduction of integrals in two-loop five-light-parton scattering amplitudes,Chin. Phys. C44(2020) 093106 [1912.09294]

    Pith/arXiv arXiv 2020

  34. [34]

    von Hippel and M

    M. von Hippel and M. Wilhelm,Refining Integration-by-Parts Reduction of Feynman Integrals with Machine Learning,JHEP05(2025) 185 [2502.05121]

    Pith/arXiv arXiv 2025

  35. [35]

    Song, T.-Z

    Z.-Y. Song, T.-Z. Yang, Q.-H. Cao, M.-x. Luo and H. X. Zhu,Explainable AI-assisted Optimization for Feynman Integral Reduction,2502.09544

    Pith/arXiv arXiv

  36. [36]

    Ablinger, J

    J. Ablinger, J. Bl¨ umlein, C. Raab, C. Schneider and F. Wißbrock,Calculating Massive 3-loop Graphs for Operator Matrix Elements by the Method of Hyperlogarithms,Nucl. Phys. B885(2014) 409 [1403.1137]

    Pith/arXiv arXiv 2014

  37. [37]

    D. A. Kosower,Direct Solution of Integration-by-Parts Systems,Phys. Rev. D98 (2018) 025008 [1804.00131]

    Pith/arXiv arXiv 2018

  38. [38]

    Feng,Generation function for one-loop tensor reduction,Commun

    B. Feng,Generation function for one-loop tensor reduction,Commun. Theor. Phys. 75(2023) 025203 [2209.09517]. – 20 –

    Pith/arXiv arXiv 2023

  39. [39]

    X. Guan, X. Li and Y.-Q. Ma,Exploring the linear space of Feynman integrals via generating functions,Phys. Rev. D108(2023) 034027 [2306.02927]

    Pith/arXiv arXiv 2023

  40. [40]

    C. Hu, T. Li, J. Shen and Y. Xu,An explicit expression of generating function for one-loop tensor reduction,JHEP02(2024) 158 [2308.13336]

    Pith/arXiv arXiv 2024

  41. [41]

    T. Li, Y. Song and L. Zhang,Solving arbitrary one-loop reduction via generating function,Eur. Phys. J. C85(2025) 142 [2404.04644]

    Pith/arXiv arXiv 2025

  42. [42]

    C. Hu, B. Feng, J. Shen and Y. Zhang,General one-loop generating function by integration-by-part relations,Phys. Rev. D111(2025) 076026 [2403.16040]

    Pith/arXiv arXiv 2025

  43. [43]

    C. Hu, Y. Hu and J. Shen,Direct expression for one-loop tensor reduction with Lorentz indices via a generating function,Phys. Rev. D111(2025) 096005 [2501.11150]

    Pith/arXiv arXiv 2025

  44. [44]

    Hu, W.-D

    C. Hu, W.-D. Li and X. Li,Generating Function of Loop Reduction by Baikov Representation,2504.02573

    arXiv

  45. [45]

    C. Hu, T. Li and X. Li,One-loop Feynman integral reduction by differential operators,Phys. Rev. D104(2021) 116014 [2108.00772]

    Pith/arXiv arXiv 2021

  46. [46]

    Feng and T

    B. Feng and T. Li,PV-reduction of sunset topology with auxiliary vector,Commun. Theor. Phys.74(2022) 095201 [2203.16881]

    Pith/arXiv arXiv 2022

  47. [47]

    B. Feng, T. Li, H. Wang and Y. Zhang,Reduction of general one-loop integrals using auxiliary vector,JHEP05(2022) 065 [2203.14449]

    Pith/arXiv arXiv 2022

  48. [48]

    B. Feng, J. Gong and T. Li,Universal treatment of the reduction for one-loop integrals in a projective space,Phys. Rev. D106(2022) 056025 [2204.03190]

    Pith/arXiv arXiv 2022

  49. [49]

    Li,Nontrivial one-loop recursive reduction relation,JHEP07(2023) 051 [2209.11428]

    T. Li,Nontrivial one-loop recursive reduction relation,JHEP07(2023) 051 [2209.11428]

    Pith/arXiv arXiv 2023

  50. [50]

    Zhang,Tensor loop reduction via the Baikov representation and an auxiliary vector,Phys

    L. Zhang,Tensor loop reduction via the Baikov representation and an auxiliary vector,Phys. Rev. D110(2024) 016013 [2309.00930]

    Pith/arXiv arXiv 2024

  51. [51]

    Gluza, K

    J. Gluza, K. Kajda and D. A. Kosower,Towards a Basis for Planar Two-Loop Integrals,Phys. Rev. D83(2011) 045012 [1009.0472]

    Pith/arXiv arXiv 2011

  52. [52]

    R. M. Schabinger,A New Algorithm For The Generation Of Unitarity-Compatible Integration By Parts Relations,JHEP01(2012) 077 [1111.4220]

    Pith/arXiv arXiv 2012

  53. [53]

    Zhang,Integration-by-parts identities from the viewpoint of differential geometry, in19th Itzykson Meeting on Amplitudes 2014, 8, 2014,1408.4004

    Y. Zhang,Integration-by-parts identities from the viewpoint of differential geometry, in19th Itzykson Meeting on Amplitudes 2014, 8, 2014,1408.4004

    Pith/arXiv arXiv 2014

  54. [54]

    Zhang,Lecture Notes on Multi-loop Integral Reduction and Applied Algebraic Geometry, 12, 2016,1612.02249

    Y. Zhang,Lecture Notes on Multi-loop Integral Reduction and Applied Algebraic Geometry, 12, 2016,1612.02249

    Pith/arXiv arXiv 2016

  55. [55]

    Singular4-4-0 — A computer algebra system for polynomial computations

    W. Decker, G.-M. Greuel, G. Pfister and H. Sch¨ onemann, “Singular4-4-0 — A computer algebra system for polynomial computations.” http://www.singular.uni-kl.de, 2024. – 21 –

    2024