Recognition: unknown
All-loop four-quark Bethe-Salpeter kernel
Pith reviewed 2026-05-08 16:22 UTC · model grok-4.3
The pith
The all-loop bare perturbative four-quark Bethe-Salpeter kernel is calculated analytically in the large-flavor limit of QCD.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We analytically calculate the all-loop bare perturbative part of the four-quark Bethe-Salpeter kernel using modern scattering amplitude methods. We work to subleading order in the large number of quark flavors approximation of massless Quantum Chromodynamics, which simultaneously makes an all-loop calculation feasible, is systematically improvable, and preserves asymptotic freedom. It also allows for avoiding the ambiguity of choosing a truncation scheme in Dyson-Schwinger equations. As a byproduct of our calculation, we also provide the result for the gluon and quark propagators.
What carries the argument
The four-quark Bethe-Salpeter kernel, obtained through Integration-By-Parts reduction of Feynman integrals to a basis of master integrals followed by direct integration to generalized polylogarithms.
If this is right
- The kernel can be used as input for nonperturbative formulations of the Bethe-Salpeter equation.
- It provides a way to study bound states without truncation ambiguities in related Dyson-Schwinger equations.
- The method opens a path to phenomenological applications in particle physics.
- The gluon and quark propagators are now known to all loops in this approximation.
Where Pith is reading between the lines
- This result might allow for more precise predictions of tetraquark masses or decay rates when combined with nonperturbative methods.
- Similar techniques could be applied to compute other multi-particle kernels in QCD.
- Comparison with lattice simulations in the large N_f limit could test the perturbative input.
Load-bearing premise
The large number of quark flavors approximation makes the all-loop calculation possible while remaining systematically improvable and preserving asymptotic freedom.
What would settle it
An independent computation of the kernel at a high loop order, such as four loops, using different methods and checking agreement with the analytic expression in terms of generalized polylogarithms.
read the original abstract
We analytically calculate the all-loop bare perturbative part of the four-quark Bethe-Salpeter kernel using modern scattering amplitude methods. We work to subleading order in the large number of quark flavors approximation of massless Quantum Chromodynamics, which simultaneously makes an all-loop calculation feasible, is systematically improvable, and preserves asymptotic freedom. It also allows for avoiding the ambiguity of choosing a truncation scheme in Dyson-Schwinger equations. We exploit state-of-the-art methods in Integration-By-Parts reduction of Lorentz scalar Feynman integrals into a minimal Master Integral basis, and direct integration into Generalized Polylogarithms. As a byproduct of our calculation, we also provide the result for the gluon and quark propagators. We discuss a path towards nonperturbative formulation and potential future phenomenological applications.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims to analytically compute the all-loop bare perturbative four-quark Bethe-Salpeter kernel in massless QCD at subleading order in the large-N_f expansion. It employs modern scattering-amplitude techniques, specifically Integration-By-Parts reduction of Feynman integrals to a minimal master-integral basis followed by direct integration into generalized polylogarithms. Results for the gluon and quark propagators are provided as a byproduct, and a path toward a nonperturbative formulation is outlined.
Significance. An explicit all-loop result for the kernel in a controlled expansion would constitute a notable technical advance, supplying a systematically improvable perturbative input for Bethe-Salpeter studies that avoids truncation-scheme ambiguities. The byproduct propagators would also be of immediate use. The significance hinges on whether the large-N_f framework truly preserves asymptotic freedom and on the verifiability of the final expression.
major comments (2)
- Abstract and introductory discussion of the approximation: the claim that the subleading large-N_f expansion 'preserves asymptotic freedom' is not self-evident. The one-loop beta-function coefficient receives a leading contribution +(2/3)N_f from quark bubbles; for N_f ≫ 1 this drives b0 < 0 and produces a UV Landau pole. The manuscript must explicitly state whether the pure-glue 11 C_A term is retained at leading order while only the kernel diagrams are expanded, or whether a different organization (e.g., fixed finite N_f with 1/N_f corrections) is employed. Without this clarification the central assertion that the result is phenomenologically relevant remains unsubstantiated.
- Main text (calculation section) and any appendices: the manuscript asserts an analytical all-loop result but supplies neither the explicit expression for the four-quark kernel nor the reduction steps to the master-integral basis. To allow verification of the central claim, the final result (in terms of generalized polylogarithms or master integrals) together with the key IBP relations and integration steps must be presented.
minor comments (1)
- The abstract mentions 'state-of-the-art methods in Integration-By-Parts reduction' but does not cite the specific reduction software or algorithm employed; adding the reference would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have revised the manuscript to improve clarity and verifiability.
read point-by-point responses
-
Referee: Abstract and introductory discussion of the approximation: the claim that the subleading large-N_f expansion 'preserves asymptotic freedom' is not self-evident. The one-loop beta-function coefficient receives a leading contribution +(2/3)N_f from quark bubbles; for N_f ≫ 1 this drives b0 < 0 and produces a UV Landau pole. The manuscript must explicitly state whether the pure-glue 11 C_A term is retained at leading order while only the kernel diagrams are expanded, or whether a different organization (e.g., fixed finite N_f with 1/N_f corrections) is employed. Without this clarification the central assertion that the result is phenomenologically relevant remains unsubstantiated.
Authors: We agree that the original wording was insufficiently precise. Our expansion is organized by retaining the full leading-order pure-glue contribution (11 C_A/3) to the beta function while expanding only the four-quark kernel diagrams to subleading order in 1/N_f. This ensures the ultraviolet behavior remains asymptotically free, with the gluonic sector dominating at leading order. We have revised the abstract and the introductory discussion to state this organization explicitly and to clarify the phenomenological relevance of the controlled approximation. revision: yes
-
Referee: Main text (calculation section) and any appendices: the manuscript asserts an analytical all-loop result but supplies neither the explicit expression for the four-quark kernel nor the reduction steps to the master-integral basis. To allow verification of the central claim, the final result (in terms of generalized polylogarithms or master integrals) together with the key IBP relations and integration steps must be presented.
Authors: We accept that the original manuscript omitted sufficient technical detail for independent verification. The complete all-loop expression for the bare four-quark kernel, written in terms of generalized polylogarithms, has been added to a new appendix. We have also included the minimal master-integral basis, the principal IBP reduction identities employed, and an outline of the direct integration procedure that yields the generalized polylogarithms. These additions should permit full reconstruction and checking of the result. revision: yes
Circularity Check
No significant circularity; derivation relies on standard perturbative methods and IBP techniques.
full rationale
The paper presents an analytic all-loop calculation of the four-quark Bethe-Salpeter kernel at subleading large-N_f order using Integration-By-Parts reduction to master integrals and direct integration to generalized polylogarithms. No load-bearing step reduces by construction to a fitted parameter, self-defined quantity, or unverified self-citation chain. The large-N_f approximation is an explicit input choice justified by feasibility and systematic improvability rather than derived from the kernel itself. Byproduct results for propagators follow from the same perturbative expansion without circular renaming or ansatz smuggling. The derivation chain is self-contained against external benchmarks of amplitude methods.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Standard perturbative expansion and Feynman rules of massless QCD
- domain assumption Large-Nf limit renders all-loop summation feasible while preserving asymptotic freedom
Reference graph
Works this paper leans on
-
[1]
D. J. Gross and F. Wilczek,Ultraviolet Behavior of Nonabelian Gauge Theories,Phys. Rev. Lett.30(1973) 1343–1346
1973
-
[2]
H. D. Politzer,Reliable Perturbative Results for Strong Interactions?,Phys. Rev. Lett.30 (1973) 1346–1349
1973
-
[3]
Gross et al.,50 Years of Quantum Chromodynamics,Eur
F. Gross et al.,50 Years of Quantum Chromodynamics,Eur. Phys. J. C83(2023) 1125, [2212.11107]
-
[4]
E. Shuryak and I. Zahed,The Hadron-Parton Bridge, From the QCD Vacuum to Partons, 2601.15085
-
[5]
Jaffe and E
A. Jaffe and E. Witten,Quantum yang-mills theory, 2000. https://www.claymath.org/wp-content/uploads/2022/06/yangmills.pdf
2000
-
[6]
’t Hooft,A Two-Dimensional Model for Mesons,Nucl
G. ’t Hooft,A Two-Dimensional Model for Mesons,Nucl. Phys. B75(1974) 461–470
1974
-
[7]
’t Hooft,A Planar Diagram Theory for Strong Interactions,Nucl
G. ’t Hooft,A Planar Diagram Theory for Strong Interactions,Nucl. Phys. B72(1974) 461
1974
-
[8]
G. F. Sterman,LEADING LOGARITHMIC BEHAVIOR OF ELECTROMAGNETIC FORM-FACTORS IN BOUND STATE MODELS,Phys. Rev. D9(1974) 3188. – 37 –
1974
-
[9]
G. P. Lepage,TWO-BODY BOUND STATES IN QUANTUM ELECTRODYNAMICS, other thesis, 7, 1978
1978
-
[10]
L. von Smekal, A. Hauck and R. Alkofer,A Solution to Coupled Dyson–Schwinger Equations for Gluons and Ghosts in Landau Gauge,Annals Phys.267(1998) 1–60, [hep-ph/9707327]
-
[11]
P. Maris and C. D. Roberts,Dyson-Schwinger equations: A Tool for hadron physics,Int. J. Mod. Phys. E12(2003) 297–365, [nucl-th/0301049]
- [12]
- [13]
-
[14]
G. van Baalen, D. Kreimer, D. Uminsky and K. Yeats,The QCD beta-function from global solutions to Dyson-Schwinger equations,Annals Phys.325(2010) 300–324, [0906.1754]
- [15]
- [16]
-
[17]
Liu,Generic framework for non-perturbative QCD in light hadrons,2501.07776
W.-Y. Liu,Generic framework for non-perturbative QCD in light hadrons,2501.07776
- [18]
-
[19]
F. V. Tkachov,A theorem on analytical calculability of 4-loop renormalization group functions,Phys. Lett. B100(1981) 65–68
1981
-
[20]
K. G. Chetyrkin and F. V. Tkachov,Integration by parts: The algorithm to calculate β-functions in 4 loops,Nucl. Phys. B192(1981) 159–204
1981
- [21]
- [22]
- [23]
- [24]
-
[25]
E. Panzer,Algorithms for the symbolic integration of hyperlogarithms with applications to Feynman integrals,Comput. Phys. Commun.188(2015) 148–166, [1403.3385]
-
[26]
Heinrich,Collider Physics at the Precision Frontier,Phys
G. Heinrich,Collider Physics at the Precision Frontier,Phys. Rept.922(2021) 1–69, [2009.00516]
-
[27]
A. Huss, J. Huston, S. Jones, M. Pellen and R. R¨ ontsch,Les Houches 2023 – Physics at TeV Colliders: Report on the Standard Model Precision Wishlist,2504.06689
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[28]
Beneke,Pole mass renormalon and its ramifications,Eur
M. Beneke,Pole mass renormalon and its ramifications,Eur. Phys. J. ST230(2021) 2565–2579, [2108.04861]. – 38 –
- [29]
-
[30]
Altarelli,Introduction to renormalons, in5th Hellenic School and Workshops on Elementary Particle Physics, pp
G. Altarelli,Introduction to renormalons, in5th Hellenic School and Workshops on Elementary Particle Physics, pp. 221–236, 1996
1996
-
[31]
M. Beneke,Renormalons,Phys. Rept.317(1999) 1–142, [hep-ph/9807443]
work page Pith review arXiv 1999
- [32]
-
[33]
Scattering amplitudes from dispersive iterations of unitarity,
P. Tourkine and A. Zhiboedov,Scattering amplitudes from dispersive iterations of unitarity, JHEP11(2023) 005, [2303.08839]
- [34]
-
[35]
Radozycki,Summing up the perturbation series in the Schwinger model,Eur
T. Radozycki,Summing up the perturbation series in the Schwinger model,Eur. Phys. J. C 6(1999) 549–553, [hep-th/9712041]
- [36]
-
[37]
Weinzierl,Hopf algebras and Dyson-Schwinger equations,Front
S. Weinzierl,Hopf algebras and Dyson-Schwinger equations,Front. Phys. (Beijing)11(2016) 111206, [1506.09119]
-
[38]
Cavalcanti,Renormalons beyond the Borel plane,Phys
E. Cavalcanti,Renormalons beyond the Borel plane,Phys. Rev. D103(2021) 025019, [2011.11175]
-
[39]
All-loop Heavy-Heavy-Light-Light correlators in N= 4 super Yang-Mills theory,
A. Brown, F. Galvagno and C. Wen,All-loop Heavy-Heavy-Light-Light correlators inN= 4 super Yang-Mills theory,JHEP10(2024) 171, [2407.02250]
-
[40]
Tropicalized quantum field theory and global tropical sampling
M. Borinsky,Tropicalized quantum field theory and global tropical sampling,2508.14263
work page internal anchor Pith review Pith/arXiv arXiv
-
[41]
P.-H. Balduf and E. Panzer,Renormalized tropical field theory,2512.21091
-
[42]
L. Clingerman and M. D. Schwartz,Asymptotic Behavior of Diagram Classes,2512.09042
-
[43]
L. Lippstreu,Analytic Properties of Infrared-Finite Amplitudes in Theories with Long-Range Forces,2505.04702. [45]Particle Data Groupcollaboration, P. A. Zyla et al.,Review of Particle Physics,PTEP 2020(2020) 083C01
-
[44]
D. Binosi and J. Papavassiliou,Gauge-invariant truncation scheme for the Schwinger-Dyson equations of QCD,Phys. Rev. D77(2008) 061702, [0712.2707]
- [45]
-
[46]
Banks,Schwinger Dyson Summation of Perturbation Expansions,2410.19683
T. Banks,Schwinger Dyson Summation of Perturbation Expansions,2410.19683
-
[47]
F. Peng and H. Shu,Truncating Dyson-Schwinger equations based on a Lefschetz thimble decomposition and Borel resummation,Phys. Rev. D111(2025) 065003, [2410.13364]
-
[48]
R. Alkofer and L. von Smekal,The Infrared behavior of QCD Green’s functions: Confinement dynamical symmetry breaking, and hadrons as relativistic bound states,Phys. Rept.353(2001) 281, [hep-ph/0007355]. – 39 –
-
[49]
D. Dorigoni,An Introduction to Resurgence, Trans-Series and Alien Calculus,Annals Phys. 409(2019) 167914, [1411.3585]
- [50]
- [51]
-
[52]
S. Caron-Huot, A. Pokraka and Z. Zahraee,Two-point sum-rules in three-dimensional Yang-Mills theory,JHEP01(2024) 195, [2309.04472]
-
[53]
A. Maiezza and J. C. Vasquez,Gluon mass generation from renormalons and resurgence, Phys. Lett. B853(2024) 138697, [2405.01639]
-
[54]
G. V. Dunne,Introductory Lectures on Resurgence: CERN Summer School 2024, 2511.15528
work page internal anchor Pith review Pith/arXiv arXiv 2024
-
[55]
Catani and L
S. Catani and L. Trentadue,Resummation of the QCD Perturbative Series for Hard Processes,Nucl. Phys. B327(1989) 323–352
1989
-
[56]
The Conformal window in QCD and supersymmetric QCD,
E. Gardi and G. Grunberg,The Conformal window in QCD and supersymmetric QCD, JHEP03(1999) 024, [hep-th/9810192]
-
[57]
M. Mari˜ no,Lectures on non-perturbative effects in largeNgauge theories, matrix models and strings,Fortsch. Phys.62(2014) 455–540, [1206.6272]
-
[58]
Witten,Baryons in the 1/n Expansion,Nucl
E. Witten,Baryons in the 1/n Expansion,Nucl. Phys. B160(1979) 57–115
1979
-
[59]
J. M. Maldacena,The LargeNlimit of superconformal field theories and supergravity,Adv. Theor. Math. Phys.2(1998) 231–252, [hep-th/9711200]
work page internal anchor Pith review arXiv 1998
-
[60]
Anti De Sitter Space And Holography
E. Witten,Anti de Sitter space and holography,Adv. Theor. Math. Phys.2(1998) 253–291, [hep-th/9802150]
work page internal anchor Pith review arXiv 1998
-
[61]
Veneziano,Some Aspects of a Unified Approach to Gauge, Dual and Gribov Theories, Nucl
G. Veneziano,Some Aspects of a Unified Approach to Gauge, Dual and Gribov Theories, Nucl. Phys. B117(1976) 519–545
1976
-
[62]
Nogueira,Automatic Feynman graph generation,J
P. Nogueira,Automatic Feynman graph generation,J. Comput. Phys.105(1993) 279–289
1993
- [63]
- [64]
-
[65]
R. N. Lee,LiteRed 1.4: a powerful tool for reduction of multiloop integrals,J. Phys. Conf. Ser.523(2014) 012059, [1310.1145]
work page Pith review arXiv 2014
-
[66]
T. Huber and D. Maˆ ıtre,HypExp, a Mathematica package for expanding hypergeometric functions around integer-valued parameters,Comput. Phys. Commun.175(2006) 122–144, [hep-ph/0507094]
-
[67]
C. Duhr and F. Dulat,PolyLogTools — polylogs for the masses,JHEP08(2019) 135, [1904.07279]
-
[68]
M. Heller and A. von Manteuffel,MultivariateApart: Generalized partial fractions,Comput. Phys. Commun.271(2022) 108174, [2101.08283]
-
[69]
X. Liu and Y.-Q. Ma,AMFlow: A Mathematica package for Feynman integrals computation via auxiliary mass flow,Comput. Phys. Commun.283(2023) 108565, [2201.11669]. – 40 –
-
[70]
Catani,The Singular behavior of QCD amplitudes at two loop order,Phys
S. Catani,The Singular behavior of QCD amplitudes at two loop order,Phys. Lett. B427 (1998) 161–171, [hep-ph/9802439]
-
[71]
Palanques-Mestre and P
A. Palanques-Mestre and P. Pascual,The 1/N −f Expansion of theγand Beta Functions in QED,Commun. Math. Phys.95(1984) 277
1984
- [72]
-
[73]
O. Antipin, N. A. Dondi, F. Sannino, A. E. Thomsen and Z.-W. Wang,Gauge-Yukawa theories: Beta functions at largeN f,Phys. Rev. D98(2018) 016003, [1803.09770]
-
[74]
Radozycki and I
T. Radozycki and I. Bialynicki-Birula,Finite nonperturbative solutions of Dyson-Schwinger equations in QED in the infrared domain,Phys. Rev. D52(1995) 2439–2445
1995
-
[75]
Williams,Schwinger-Dyson equations in QED and QCD: The Calculation of fermion-antifermion condensates, other thesis, 6, 2007
R. Williams,Schwinger-Dyson equations in QED and QCD: The Calculation of fermion-antifermion condensates, other thesis, 6, 2007
2007
-
[76]
G. Eichmann, H. Sanchis-Alepuz, R. Williams, R. Alkofer and C. S. Fischer,Baryons as relativistic three-quark bound states,Prog. Part. Nucl. Phys.91(2016) 1–100, [1606.09602]
-
[77]
Hoyer,Lectures on Bound states, 5, 2016.1605.01532
P. Hoyer,Lectures on Bound states, 5, 2016.1605.01532
-
[78]
Pagels and S
H. Pagels and S. Stokar,The Pion Decay Constant, Electromagnetic Form-Factor and Quark Electromagnetic Selfenergy in QCD,Phys. Rev. D20(1979) 2947
1979
- [79]
-
[80]
D. Binosi and L. Theussl,JaxoDraw: A Graphical user interface for drawing Feynman diagrams,Comput. Phys. Commun.161(2004) 76–86, [hep-ph/0309015]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.