Pith. Sign in

REVIEW

Two-Loop Quark and Gluon Form Factors in Dimensional Regularisation

High Energy Physics - Phenomenology · Submitted 5 Jul 2005 · 22 references

Pith's verdict: The exact crossed-triangle integral A6 is the real deliverable; the ε² gluon form factor is new but only numerically checked, and the check is reported in one acknowledgement sentence.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 1 major objection · 8 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0507061 (version 3)
Pith IDZCPYU36F, Pith's permanent record of this paper and its review
CitationD. Maitre, T. Gehrmann, T. Huber. "Two-Loop Quark and Gluon Form Factors in Dimensional Regularisation." arXiv:hep-ph/0507061 (2005). https://pith.science/paper/ZCPYU36F
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

The two-loop quark and gluon form factors are obtained in closed form to all orders in the dimensional-regularisation parameter, after the last missing two-loop master integral is expressed through generalised hypergeometric functions of单位.

The reading

In dimensional regularisation the dimension of space-time is continued away from four, $d=4-2epsilon$, and every loop integral becomes a Laurent series in $epsilon$. Physical predictions need only the first few terms of that series, but higher-order terms are exactly what infrared factorisation formulas consume: the pole operators at three loops act on the $epsilon^2$ pieces of the two-loop form factors. The obstacle was that one of the three two-loop vertex master integrals, the crossed triangle $A_6$, had no exact expression — only a series known up to the finite term. The authors give $A_6$ in closed form as $Gamma$-functions plus generalised hypergeometric functions of unit argument, expand it with a purpose-built package, and then write out the one-loop form factors to $epsilon^4$ and the two-loop quark and gluon form factors to $epsilon^2$. If the expansions are right, the two-loop form factors are no longer limited by integral evaluation at all: any order in $epsilon$ can be produced on demand.

What carries the argument

The load-bearing object is the exact expression for $A_6$ in Eq. (5): a sum of four terms built from $Gamma$-functions of $epsilon$-dependent arguments times ${}_3F_2(1,-4epsilon,-2epsilon;1-3epsilon,1-2epsilon;1)$, ${}_3F_2(1,1,1+2epsilon;2+epsilon,2+2epsilon;1)$ and ${}_4F_3(1,1-epsilon,-4epsilon,-2epsilon;1-3epsilon,1-2epsilon,1-2epsilon;1)$. Generalised hypergeometric functions of unit argument have parameters that here carry the regulator $epsilon$, so expanding in $epsilon$ means expanding in the parameters of the function, not just in its argument; that is what the HypExp package does, and it is the mechanism that converts a closed-form integral into a Laurent series to any desired order. The second piece of machinery is the reduction of the full two-loop form factors onto the master integrals $A_{2,LO}$, $A_3$, $A_4$ and $A_6$ using integration-by-parts and Lorentz-invariance identities, which is carried out exactly in $d$ so that no accuracy is lost before the final expansion.

Extended reading notes

Core claim

The paper's central claim is that the two-loop crossed-triangle integral $A_6 = int frac{d^dk}{(2pi)^d}frac{d^dl}{(2pi)^d} frac{1}{k^2 l^2 (k-p_1-p_2)^2 (k-l)^2 (k-l-p_2)^2 (l-p_1)^2}$ admits an exact expression in $d$ dimensions as a combination of $Gamma$-functions and a ${}_3F_2$ and a ${}_4F_3$ of unit argument (Eq. 5). With this, every master integral entering the massless quark form factor $gamma^* to qbar q$ and the gluon form factor $H to gg$ is known to all orders in $epsilon$, so the form factors themselves are known to all orders once the reduction coefficients are inserted. The authors give the eighth-order Laurent expansion of $A_6$ (Eq. 6), the renormalised one-loop form factors through $epsilon^4$ (Eqs. 15 and 16), and the renormalised two-loop form factors through $epsilon^2$ (Eqs. 17 and 18), the latter being the order required for the ultraviolet renormalisation and infrared factorisation of the corresponding three-loop quantities. A note added at the end states that expanding the unrenormalised two-loop quark form factor to $epsilon^2$ reproduces Eq. (B.1) of the concurrent work by Moch, Vermaseren and Vogt.

Load-bearing premise

Everything rests on the correctness of the HypExp-generated Laurent expansions of the ${}_3F_2$ and ${}_4F_3$ functions in Eq. (5), and the only validation stated for those expansions is an acknowledgement of independent numerical checks by sector decomposition, with no numbers, error estimates or comparison tables shown.

Editorial extensions

If this is right

  • The $epsilon^2$ coefficients of the two-loop quark and gluon form factors supply exactly the input that the three-loop infrared singularity operators project onto, so the previously unknown hard radiation constants in the three-loop infrared factorisation formula for these form factors become fixable.
  • Because the two-loop vertex master integrals reappear as subtopologies in three-loop reductions whenever a loop is pinched off, the $epsilon^2$ terms given here can be inserted directly into three-loop form-factor calculations.
  • Any other integral that reduces to $A_6$ can now be expanded to arbitrary order in $epsilon$ rather than only to the finite term, since the closed form plus HypExp removes the order-by-order bottleneck.
  • The paper argues the same technique should apply to three- and four-loop two-point functions and to multi-particle phase-space integrals, which are cuts of such two-point functions, for fully inclusive observables.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The methodological bet worth watching is that exact-in-$d$ master integrals plus a parameter-expansion package beat order-by-order numerical sector decomposition once high $epsilon$ orders are needed; if that holds generally, the bottleneck in multi-loop work shifts from evaluating integrals to reducing diagrams onto them.
  • The quark half of the result is cross-checked against an independent published expansion, but the gluon half is not: the paper offers no analogous external comparison for Eq. (18), so the gluon $epsilon^2$ coefficients are the more exposed part of the claim.
  • The appearance of the multiple zeta value $zeta_{5,3}$ already at order $epsilon^4$ in $A_6$ hints that pushing to higher orders will require a working calculus of multiple zeta values, which is a real constraint on any 'expand to any order' promise.
  • Since the form factors are now known as exact functions of $d$, one could ask whether the $epsilon to 0$ limit of the product of two-loop form factors and the infrared operators exponentiates in a way that is visible before expansion, which would turn the order-by-order matching into a structural statement rather than a coefficient hunt.
Original abstract

We compute the two-loop corrections to the massless quark form factor $\gamma^* \to q\bar q$ and gluon form factor $H\to gg$ to all orders in the dimensional regularisation parameter $\epsilon=(4-d)/2$. The two-loop contributions to the form factors are reduced to linear combinations of master integrals, which are computed in a closed form, expressed as $\Gamma$-functions and generalised hypergeometric functions of unit argument. Using the newly developed HypExp-package, these can be expanded to any desired order, yielding Laurent expansions in $\epsilon$. We provide expansions of the form factors to order $\epsilon^2$, as required for ultraviolet renormalisation and infrared factorisation of the three-loop form factors.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 3,723 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.17 · human participants 0.01 · controlled comparison 0.01 · randomized assignment 0.01 · states quantitative result 0.36 · compares with prior results 0.05 · mentions code or data availability 0.23.

The abstract, sentence by sentence

arXiv listing abstract v3. Results and the main claim are underlined; .

We compute the two-loop corrections to the massless quark form factor $\gamma^* \to q\bar q$ and gluon form factor $H\to gg$ to all orders in the dimensional regularisation parameter $\epsilon=(4-d)/2$. The two-loop contributions to the form factors are reduced to linear combinations of master integrals, which are computed in a closed form, expressed as $\Gamma$-functions and generalised hypergeometric functions of unit argument. Using the newly developed HypExp-package, these can be expanded to any desired order, yielding Laurent expansions in $\epsilon$. We provide expansions of the form factors to order $\epsilon^2$, as required for ultraviolet renormalisation and infrared factorisation of the three-loop form factors.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We compute the two-loop corrections to the massless quark form factor $\gamma^* \to q\bar q$ and gluon form factor $H\to gg$ to all orders in the dimensional regularisation parameter $\epsilon=(4-d)/2$. contribution 1.00 0.96 0.41 asserted 1.00
s1 The two-loop contributions to the form factors are reduced to linear combinations of master integrals, which are computed in a closed form, expressed as $\Gamma$-functions and generalised hypergeometric functions of unit argument. method 0.76 0.94 0.22 asserted 1.00
s2 Using the newly developed HypExp-package, these can be expanded to any desired order, yielding Laurent expansions in $\epsilon$. method 0.99 0.73 0.07 asserted 0.96
s3 We provide expansions of the form factors to order $\epsilon^2$, as required for ultraviolet renormalisation and infrared factorisation of the three-loop form factors. result 0.73 0.95 0.34 asserted 1.00

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 0.88.

What was examined, and what was not

sources read
  • arXiv listing abstract v3
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References
How later papers used it

94 later papers cite it

Evidence: substantial. 86 of the citing papers share no author with it, and 32 of those carry a confident label from Pith's citation classifier.

Pith read 11 of the 121 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its data or values 1 of 22 papers with this label
How the literature uses it

93 arXiv papers cite this in 120 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 72 sentences without a confident label

Uses its method or tool 6 papers
  • The back-to-back asymptotics of the EEC can be described using Soft and Collinear Effective Theory (SCET) [ref] via the following factorization theorem [ref] [equation] In EEC_fact_thm_q, $J_0$ is the Bessel function arising from the Fourier transform due to the azymuthal symmetry of the EEC measurement, $H_qq̅$ is the quark color singlet SCET hard function, which is related to the IR finite part of the quark form factors [ref] and can be extracted up to 4 loops from the recent result of Lee:2022nhh, and $_q$ is the quark EEC jet function which is known up to N$^3$LO [ref] .

    The Four-Loop Rapidity Anomalous Dimension and Event Shapes to Fourth Logarithmic Order · EEC in the back-to-back limit · 0.87

  • Amongst the resulting terms are integrals in $(D+2n)$ dimensions (where $n ∈ℤ$) and with raised powers of propagators. * Higher-dimensional scalar integrals are shifted back to $D$ dimensions using dimensional recurrence relations [ref] . * Using the Mathematica package LiteRed [ref] , the resulting $D$-dimensional scalar integrals are reduced to a basis of masters using integration-by-parts relations (IBPs) [ref] . * Known expressions for the four-point two-loop master integrals [ref] are inserted.

    Infrared and transcendental structure of two-loop supersymmetric QCD amplitudes · Integration & transcendental weight · 1.00

  • In this way the three-loop vertex integrals are reduced to two-, and three-loop vacuum integrals which are treated with MATAD [ref] , and to massless one- and two-loop vertex contributions which are reduced to master integrals (see, e.g., Ref. [ref] ) with the help of FIRE [ref] and an independent program based of the Laprota-method.

    Three-loop vertex corrections · Virtual NNLO corrections to Higgs production in gluon fusion · 1.00

Uses its data or values 21 papers
Compares with it 4 papers
Cites it as background 12 papers
  • The three-loop form factor master integrals were first computed numerically to $𝒪(ϵ^2)$ in reference [ref] using dimensional recurrence relations [ref] but, with the help of the celebrated PSLQ algorithm [ref] , the authors were able to recover the analytical solutions.[At the time, it was conjectured that zeta and multiple zeta values of, at most, weight eight would appear in the higher-order results required to expand the three-loop form factors through to $𝒪(ϵ^2)$.] The explicit higher-order results for the three-loop masters are, regrettably, scattered over three different articles [ref] .

    On the Computation of Form Factors in Massless QCD with Finite Master Integrals · Notation and Conventions · 0.81

  • The two- and three-loop QCD corrections to $F_q$ and $F_g$ are available from Refs. [ref] , At four loops, only partial results have been obtained so far.

    Quark and gluon form factors in four-loop QCD · Introduction · 0.84

  • Within QCD such form factors are known up to three loops for some years now [ref] , where a vector current with no association to the gauge group has been used for the fermions.

    Infrared Anomalous Dimensions at Three-Loop in the SM from Conserved Currents · Introduction · 0.89

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.