Pith. Sign in

REVIEW

An algebraic/numerical formalism for one-loop multi-leg amplitudes

High Energy Physics - Phenomenology · Submitted 27 Apr 2005 · 88 references

Pith's verdict: A careful 2005 methods paper whose algebraic core holds; the caveats are stated by the authors rather than hidden, and the stress-test's target is real but smaller than presented.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 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 →

· 3 major objections · 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/0504267 (version 3)
Pith IDJPQPGAIJ, Pith's permanent record of this paper and its review
CitationC.Schubert, E.Pilon, G.Heinrich, J.Ph.Guillet, T.Binoth. "An algebraic/numerical formalism for one-loop multi-leg amplitudes." arXiv:hep-ph/0504267 (2005). https://pith.science/paper/JPQPGAIJ
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 paper sets out a reduction scheme that writes any one-loop Feynman diagram, with any number of external legs and with massive or massless internal lines, as a fixed short list of triangle and box integrals, without the spurious Gram-mim

The reading

The paper is trying to establish a bookkeeping scheme for one-loop Feynman diagrams that works for an arbitrary number $N$ of external legs, for massive as well as massless particles, and that never produces the spurious inverse Gram determinants which make conventional momentum-space reduction numerically unstable near degenerate kinematics. The authors argue three things at once: that a subtraction identity reorganises each tensor integral into a finite remainder plus a lower-rank, lower-point integral, so that all infrared poles end up inside three-point functions and can be handled analytically; that for $N$ at least 5 no integral in more than $n+4$ dimensions is ever needed; and that by choosing Feynman-parameter integrals with non-trivial numerators as the endpoints, the form factors can be written with no inverse Gram determinant anywhere. If this holds, a next-to-leading-order calculation for a multi-jet or multi-particle process becomes a mechanical lookup of form factors plus a numerical evaluation of a small basis, which is precisely the bottleneck that keeps multi-leg NLO cross sections from being automated. The paper pairs the algebra with a numerical fallback, a multi-dimensional contour deformation in Feynman parameter space, so that the same basis can be evaluated even where the analytic expressions degenerate.

What carries the argument

The central object is the shift-invariant difference vector $(r_i - r_j)$, which replaces the loop momentum in the numerator and makes the whole reduction invariant under $k$ to $k + r_0$; iterating the reduction therefore never requires redefining the loop momentum or proliferating terms. Around it sit the kinematic matrix $S$ and its Moore-Penrose pseudo-inverse, used so that the defining equation still has a solution when $S$ is singular (that is, for $N$ at least 7 or at exceptional kinematics), the reduction coefficients $b_i$ and $C^{j}_{a}$, and the special five-point identity $T^{}_{ab} = 2 V^{}_a V^{}_b B$ that converts the would-be $1/B$ singularities into a factor of order $O(1-n/4)$. The basis the reduction lands on is deliberately chosen: $n$- and $(n+2)$-dimensional three-point functions with up to three Feynman parameters in the numerator, and $(n+2)$- and $(n+4)$-dimensional four-point functions with up to three and one parameters respectively.

Extended reading notes

Core claim

The load-bearing result is the subtraction identity that defines the reduction coefficients $C^{j}_{a}$ through $S_{ij} C^{j}_{a} = (r_i - r_j)^{}_{a}$, where $S_{ij} = (r_i-r_j)^2 - m_i^2 - m_j^2$ is the kinematic matrix. Solving it splits a rank-$r$ $N$-point tensor integral into an infrared-finite piece plus a sum of rank-$(r-1)$, $(N-1)$-point integrals. For $N$ at least 6 the two auxiliary quantities $V^{}_b$ and $T^{}_{ab}$ both vanish in four dimensions, so only the pinched terms survive and the reduction of any $N$-point integral beyond five legs is trivial. For $N=5$ the identity $T^{}_{ab} = 2 V^{}_a V^{}_b B$ (with $B$ the sum of reduction coefficients) makes the coefficient of the higher-dimensional pentagon of order $O(1-n/4)$, and appendix C proves by induction that the accompanying $1/B$ terms cancel analytically, leaving form factors free of both higher-dimensional pentagons and inverse Gram determinants. The paper states this explicitly as a representation not previously given in the literature, and supplies the full form factor lists for $N=1$ through $N=5$.

Load-bearing premise

The claim that the numerical contour-deformation evaluation stays reliable at degenerate kinematics is demonstrated on a single two-to-two scattering example with one interpolation parameter, and is extrapolated from there to forward and backward scattering and to other mass configurations.

Editorial extensions

If this is right

  • Any one-loop $N$-point amplitude can be expressed as a linear combination of a fixed basis of two-point functions, $n$- and $(n+2)$-dimensional three-point functions, and $(n+2)$- and $(n+4)$-dimensional four-point functions with Feynman parameters in the numerator.
  • Infrared poles appear only inside three-point functions, so the divergent part of an amplitude can be extracted analytically and separated from the finite remainder before any numerical integration is attempted.
  • Tensor integrals with six or more external legs reduce iteratively to five-point integrals plus reduction coefficients, so no form factors beyond the five-point case need ever be derived.
  • A practical evaluator can be built as a hybrid: the fast analytic representation wherever the Gram determinant is comfortably non-zero, and the contour-deformation numerics only in the thin shell of phase space where kinematics degenerate.

Reading between the lines

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

  • The guarantee that higher-dimensional pentagons drop out is stated for rank up to five; the paper itself notes that rank-six pentagons are ultraviolet divergent, so the order-$\epsilon$ term can no longer be discarded and an extra rational piece appears in eq. (C.105). Processes with effective Higgs-gluon couplings or higher-rank numerators therefore sit outside the clean version of the claim.
  • Avoiding inverse Gram determinants removes only the spurious denominators. The paper concedes that $\delta V \cdot K \cdot \delta V$ and $\det S$ can still vanish at genuine soft, collinear, or threshold configurations, so the stability gain is confined to the class of artefacts rather than to physical singularities.
  • The decision rule the paper proposes, a threshold $B\Lambda^2 \ge 1$ switching between analytic and numerical evaluation, is effectively a design specification for a software package: the expensive numerical path is invoked only in a small fraction of phase space, which is what makes the cost acceptable.
  • The helicity-decomposition relations in appendix E, where the six-point reduction coefficients collapse to spinor-trace ratios, suggest that a dedicated massless six-leg implementation could sidestep the general machinery entirely for that case.
Original abstract

We present a formalism for the calculation of multi-particle one-loop amplitudes, valid for an arbitrary number N of external legs, and for massive as well as massless particles. A new method for the tensor reduction is suggested which naturally isolates infrared divergences by construction. We prove that for N>4, higher dimensional integrals can be avoided. We derive many useful relations which allow for algebraic simplifications of one-loop amplitudes. We introduce a form factor representation of tensor integrals which contains no inverse Gram determinants by choosing a convenient set of basis integrals. For the evaluation of these basis integrals we propose two methods: An evaluation based on the analytical representation, which is fast and accurate away from exceptional kinematical configurations, and a robust numerical one, based on multi-dimensional contour deformation. The formalism can be implemented straightforwardly into a computer program to calculate next-to-leading order corrections to multi-particle processes in a largely automated way.

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, 5,599 input tokens, $0.0002 reported.

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

The abstract, sentence by sentence

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

We present a formalism for the calculation of multi-particle one-loop amplitudes, valid for an arbitrary number N of external legs, and for massive as well as massless particles. A new method for the tensor reduction is suggested which naturally isolates infrared divergences by construction. We prove that for N>4, higher dimensional integrals can be avoided. We derive many useful relations which allow for algebraic simplifications of one-loop amplitudes. We introduce a form factor representation of tensor integrals which contains no inverse Gram determinants by choosing a convenient set of basis integrals. For the evaluation of these basis integrals we propose two methods: An evaluation based on the analytical representation, which is fast and accurate away from exceptional kinematical configurations, and a robust numerical one, based on multi-dimensional contour deformation. The formalism can be implemented straightforwardly into a computer program to calculate next-to-leading order corrections to multi-particle processes in a largely automated way.

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 present a formalism for the calculation of multi-particle one-loop amplitudes, valid for an arbitrary number N of external legs, and for massive as well as massless particles. contribution 1.00 0.94 0.28 asserted 1.00
s1 A new method for the tensor reduction is suggested which naturally isolates infrared divergences by construction. contribution 0.66 0.90 0.13 hedged possibility 0.86
s2 We prove that for N>4, higher dimensional integrals can be avoided. result 0.91 0.94 0.17 asserted 1.00
s3 We derive many useful relations which allow for algebraic simplifications of one-loop amplitudes. contribution 0.93 0.96 0.11 asserted 1.00
s4 We introduce a form factor representation of tensor integrals which contains no inverse Gram determinants by choosing a convenient set of basis integrals. contribution 0.59 0.95 0.12 asserted 1.00
s5 For the evaluation of these basis integrals we propose two methods: An evaluation based on the analytical representation, which is fast and accurate away from exceptional kinematical configurations, and a robust numerical one, based on multi-dimensional contour deformation. method 0.94 0.88 0.09 asserted 0.81
s6 The formalism can be implemented straightforwardly into a computer program to calculate next-to-leading order corrections to multi-particle processes in a largely automated way. contribution 0.45 0.49 0.23 asserted 0.82

What kind of paper the classifier reports this is: method or tool 1.00 · evidence from proof or derivation 1.00.

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

Figures
Figure 1. General N-point one-loop graph with momentum and propagator labelling. pi are defined as incoming. Momentum conservation implies X N i=1 pi = 0 . (3) For future reference we label each propagator qi by the number i as shown in
Figure 2. Graphical representation of pinch integrals. Each topolo
Figure 3. Real and imaginary parts of the basis integrals
Figure 4. Same as Fig. 7.2 but for the basis integrals
Figure 5. Reduction of N-point tensor integrals to basis integrals. The indicated equations can be implemented directly into an algebraic computer program. I n N (1|j1|j1, j2|j1, j2, j3) denotes the integral I n N with zero, one, two or three Feynman parameters in the numerator. 36
Figure 6. Evaluation of the basis integrals: the triangle case,
Figure 7. Evaluation of the basis integrals: the box case,

Figures reproduced from arXiv:hep-ph/0504267 by the authors.

References
How later papers used it

56 later papers cite it

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

Pith read 20 of the 95 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 method or tool 6 of 18 papers with this label · 2 by the same authors

By the same authors

Compares with it 1 of 2 papers with this label
How the literature uses it

52 arXiv papers cite this in 91 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 · 51 sentences without a confident label

Builds on it 1 paper
Uses its method or tool 16 papers
Compares with it 2 papers
  • We have compared the results with an independent implementation using FeynArts and FormCalc [ref] to generate and simplify the diagrams, where the tensor integrals are algebraically reduced to scalar integrals using the procedure described in [ref] .

    NLO Cross Sections for the LHC using GOLEM: Status and Prospects · NLO Results for $qq̅→bb̅bb̅$ CITEMARK22Z · 0.89

  • Finally, we mention that, always with the help of HELAC-1loop, we successfully compared our predictions for six-quark 1-loop amplitudes (with three different flavours) with the results produced by the GOLEM group [ref] .

    Feynman Rules for the Rational Part of the QCD 1-loop amplitudes · Results · 1.00

Cites it as support 1 paper
  • Although higher-dimensional pentagons ($I_5^n+2$, $I_5^n+4$, …) formally appear in the reduction, too, it has been shown in [ref] that the coefficients of these integrals are always of order $𝒪(ε)$ and therefore drop out in phenomenological calculations.

    Automated Evaluation of One-Loop Six-Point Processes for the LHC · Basis Integrals · 0.98

Cites it as background 12 papers

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.