Pith. Sign in

REVIEW

On the evaluation of a certain class of Feynman diagrams in x-space: Sunrise-type topologies at any loop order

High Energy Physics - Phenomenology · Submitted 28 Jun 2005 · 196 references

Pith's verdict: A solid, useful methods review of the authors' own x-space program, with a handful of new results attached; the only real soft spot is the unchecked five-loop finite part.

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 · 10 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0506286 (version 1)
Pith IDEW67NLX5, Pith's permanent record of this paper and its review
CitationA.A. Pivovarov, J.G. K\"orner, S. Groote. "On the evaluation of a certain class of Feynman diagrams in x-space: Sunrise-type topologies at any loop order." arXiv:hep-ph/0506286 (2005). https://pith.science/paper/EW67NLX5
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

Sunrise Feynman diagrams at any loop order reduce to a single one-dimensional integral when written in configuration space, with divergent parts analytic and finite parts numeric.

The reading

The paper argues that a whole family of multi-loop Feynman diagrams — the sunrise (or sunset/banana) topologies, where $n+1$ propagators connect two external points — becomes tractable at every loop order if one refuses to work in momentum space. In configuration space the diagram is just a product of propagators with no integration at all, and the Fourier transform that would normally require $n$ entangled loop integrals collapses into one radial integral against a Bessel-function kernel. On that representation the authors build a complete practical toolkit: the divergent parts come out analytically for arbitrary masses, arbitrary space-time dimension and any number of loops, while the finite parts require only a single numerical quadrature. They also show that the spectral density (equivalently, the multi-particle phase space) follows from an inverse Bessel transform along a contour, and they work out threshold, large-momentum, heavy-mass and odd-dimensional limits. The claim matters because it turns a class of diagrams that grows combinatorially harder with each loop into a routine one-dimensional computation, and because it gives a uniform handle on phase space for multi-particle final states.

What carries the argument

The load-bearing object is the Bessel-transform representation of the Fourier-transformed diagram, $2π^{λ+1}∫_0^∞ (px/2)^{-λ}J_λ(px)∏_i D(x,m_i)x^{2λ+1}dx$, together with its inversion: the spectral density is $ρ(s)=(2π)^{λ}/(i s^{λ/2})∫_{c-i∞}^{c+i∞}I_λ(x√s)Π(x)x^{λ+1}dx$, a contour integral that stays finite because the contour bypasses small $x$. Renormalization is handled by adding local counterterms $C_i$ times $δ(x)$ and its derivatives, with the coefficients extracted by expanding the propagators at small $x$ under a Gaussian damping factor and integrating term by term with Gamma functions. For near-threshold behaviour the authors split the kernel into Hankel functions $H^{±}$: the $H^{+}$ piece maps onto a vacuum bubble with one extra line and is analytic at threshold, while the $H^{-}$ piece carries the cut and, after double asymptotic expansion, reduces to $∑_j A_j Γ(1-a-j+2ε)Δ^{a+j-1-2ε}$ with $Δ=M-p$.

Extended reading notes

Core claim

The central claim is that for sunrise-type topologies the configuration-space representation makes the loop count almost irrelevant. Writing each propagator as $D(x,m)=(mx)^{λ}K_{λ}(mx)/((2π)^{λ+1}x^{2λ})$ with $D=2λ+2$, the $n$-loop diagram is the plain product $Π(x)=∏_{i}D(x,m_i)$, and its Fourier transform becomes the single integral $2π^{λ+1}∫_0^∞ (px/2)^{-λ}J_λ(px)Π(x)x^{2λ+1}dx$. From this one line the authors derive: analytic pole parts for any mass configuration and any loop order (via small-$x$ expansion with a damping factor $e^{-μ^2x^2}$ and Euler Gamma integrals); finite parts as one numerical integration after momentum subtraction; a contour-integral formula for the spectral density from the inverse $K$-transform; closed-form results in odd dimensions where the propagators degenerate to exponentials; and a reduction scheme for high propagator powers that funnels everything into two master integrals. Worked checks run from the two-loop genuine sunrise through three-, four- and five-loop vacuum bubbles, reproducing known momentum-space results where those exist and producing new ones where they do not.

Load-bearing premise

The finite-part numbers are trustworthy only if the small-$x$ subtraction is accurate enough that the analytically integrated and numerically integrated pieces cancel cleanly; the authors themselves note in Sec. 3.8 that for the five-loop bubble the cancellation is severe and their subtraction choice is non-optimal, and no uncertainty estimate is attached to the quoted digits.

Editorial extensions

If this is right

  • Any sunrise-type diagram, at any loop order and with any mass assignment, can be renormalized with counterterms obtained analytically and a finite part obtained from one numerical integral, so the cost of adding loops is a constant factor rather than a combinatorial explosion.
  • The spectral density of an $(n+1)$-line sunrise diagram is exactly the $(n+1)$-particle phase space, so multi-body decay phase space volumes become computable by the same one-dimensional contour integral, including in non-integer dimensions.
  • In odd space-time dimensions the propagators become exponentials and the whole diagram evaluates in closed form; for $D=3$ the spectral density is $ρ(s)=(\sqrt{s}-M)^{n-1}\theta(s-M^2)/(2(4π)^n(n-1)!\sqrt{s})$.
  • Near threshold the expansion coefficients are generated mechanically, reproducing the momentum-space region-separation coefficients for the equal-mass three-line case and extending them to four lines and beyond.
  • When one mass is much lighter than the others the plain threshold series fails at $E\sim 2m_0$; keeping that one propagator exact and expanding only the rest extends validity up to $E\sim M$ in terms of hypergeometric functions.

Reading between the lines

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

  • The same representation should let one assemble the spectral densities of nested topologies by convolution, since the phase-space convolution $Φ_n=\int Φ_k Φ_p Φ_1$ stays inside the class of functions the authors already control.
  • Because the finite parts are one-dimensional quadratures over smooth integrands, the method is a natural cross-check engine for automated momentum-space reducers: any disagreement between the two pipelines localizes a bug in the reduction tables rather than in the physics.
  • The stated cancellation problem between the analytic and numerical pieces suggests the practical bottleneck is not the loop order but the choice of subtraction point; a subtraction matched to the dominant small-$x$ asymptotics, rather than a generic Taylor subtraction at $p=0$, would likely recover many more digits at fixed working precision.
  • The odd-dimensional closed forms give rigorous bracketing for even-dimensional results, since $K_{3/2}>K_1>K_{1/2}$ on the positive axis; this could be turned into certified error bars for four-dimensional numerics.
Original abstract

We review recently developed new powerful techniques to compute a class of Feynman diagrams at any loop order, known as sunrise-type diagrams. These sunrise-type topologies have many important applications in many different fields of physics and we believe it to be timely to discuss their evaluation from a unified point of view. The method is based on the analysis of the diagrams directly in configuration space which, in the case of the sunrise-type diagrams and diagrams related to them, leads to enormous simplifications as compared to the traditional evaluation of loops in momentum space. We present explicit formulae for their analytical evaluation for arbitrary mass configurations and arbitrary dimensions at any loop order. We discuss several limiting cases of their kinematical regimes which are e.g. relevant for applications in HQET and NRQCD.

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, 4,291 input tokens, $0.0002 reported.

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

The abstract, sentence by sentence

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

We review recently developed new powerful techniques to compute a class of Feynman diagrams at any loop order, known as sunrise-type diagrams. These sunrise-type topologies have many important applications in many different fields of physics and we believe it to be timely to discuss their evaluation from a unified point of view. The method is based on the analysis of the diagrams directly in configuration space which, in the case of the sunrise-type diagrams and diagrams related to them, leads to enormous simplifications as compared to the traditional evaluation of loops in momentum space. We present explicit formulae for their analytical evaluation for arbitrary mass configurations and arbitrary dimensions at any loop order. We discuss several limiting cases of their kinematical regimes which are e.g. relevant for applications in HQET and NRQCD.

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 review recently developed new powerful techniques to compute a class of Feynman diagrams at any loop order, known as sunrise-type diagrams. contribution 0.68 0.20 0.14 asserted 0.97
s1 These sunrise-type topologies have many important applications in many different fields of physics and we believe it to be timely to discuss their evaluation from a unified point of view. background 0.96 0.11 0.07 asserted 0.84
s2 The method is based on the analysis of the diagrams directly in configuration space which, in the case of the sunrise-type diagrams and diagrams related to them, leads to enormous simplifications as compared to the traditional evaluation of loops in momentum space. method 0.99 0.44 0.25 asserted 1.00
s3 We present explicit formulae for their analytical evaluation for arbitrary mass configurations and arbitrary dimensions at any loop order. result 0.67 0.93 0.60 asserted 1.00
s4 We discuss several limiting cases of their kinematical regimes which are e.g. relevant for applications in HQET and NRQCD. contribution 0.85 0.80 0.06 asserted 0.79

What kind of paper the classifier reports this is: review 0.70 · evidence from proof or derivation 0.97.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
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. Genuine sunrise (a) and general topology of the class of s
Figure 2. The spectral density for the four-line sunrise-type diag
Figure 3. Three-loop massive sunrise-type bubble. A dot on a line cor
Figure 4. Comparison of functions fn(z) for different values of n, plotted on a logarithmic scale Using these relations, we obtain K1−ε(z) = − X∞ n=0 ε 2n (2n)! d dz + ε z ! f2n(z). (65) The family of functions fk(z) is rather close to the original set of Bessel functions Kν(z) (cf
Figure 5. Integrands for numerical integration in case of the four
Figure 6. A five loop vacuum bubble loop case since there exist no results on the five-loop bubble in the literature. The integral representation of the five-loop bubble with equal masses m is given by Π˜ 5(p 2 = 0) = Z D(x, m) 6 d Dx. (82) Evaluating numerically the analytical part one obtains N5Π˜ ana 5 (0) = m8  3ε −5 + 16.5ε −4 + 51.95833ε −3 + 125.6715ε −2 + 259.9876ε −1 −1360392.5934 − 16888723.177ε − 111392297.46ε 2 −518606741.1ε 3 + O(ε 4 )  (83) while the numerical integration of the nonsingular part gives N5Π˜ num 5 (0) = m8  1360739.9485 + 16886269.683ε
Figure 7. Example of a multi-body decay of a quark in the Standard Mo
Figure 8. Leading order contribution to the two-point correlator in
Figure 9. Representation of Π˜ +(p) as a standard vacuum bubble with an additional line. The cross on the line denotes an arbitrary number of derivatives of that line. The two parts Π˜ ±(p) of the polarization function Π( ˜ p) have a completely different behaviour near threshold which allows one to analyze them independently. This observation makes the subsequent analysis straightforward. We first consider the contribution of the Π˜ +(p) part which reduces to a regular sunrise-type diagram. Indeed, using the relation given by Eq. (A14), Kλ(z) = πi 2 e iλπ/2H + λ (iz) (96) between Bessel…
Figure 10. Various results for the spectral density for
Figure 11. Various solutions for the one-loop spectral density for
Figure 12. A non-factorizable diagram at NLO for B0 − B¯0 mixing m1 q1 m2 q2 mn-1 qn-1 mn qn
Figure 13. Sunrise-type diagram with additional external momenta
Figure 14. The diagrams for BN (0, 0, 2, 2, 2, 1), BN (0, 0, 2, 2, 2, 2), and BN (0, 0, 2, 3, 3, 4). D(x, m) we have to use D (1)(x, m) = Z d Dp (2π) D e ipµx µ (p 2 + m2 ) 2 = (x/m) 1−λ 2(2π) λ+1 Kλ−1(mx) = (x/m) ε 2(2π) 2−ε K−ε(mx). (185) One obtains (ξ = mx) BN (0, 0, 2, 2, 2, 2) = 2 1−2ε (Γ(1 + ε))3Γ(2 − ε) Z ∞ 0 (K−ε(ξ))4 ξ 3+2ε dξ. (186) We now can use the general formula [135] " ∂Kν (z) ∂ν # ν=±n = ± 1 2 n! nX−1 k=0  z 2 k−n Kk(z) k!(n − k) , n ∈ {0, 1, . . .} (187) to expand the Bessel function in a series with respect to its index. In case of
Figure 15. spectacle topology diagram in two different representat
Figure 16. Irreducible loop replacing one of the propagators

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

References
How later papers used it

49 later papers cite it

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

Pith read 8 of the 61 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 1 of 2 papers with this label
How the literature uses it

48 arXiv papers cite this in 60 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 · 24 sentences without a confident label

Uses its method or tool 2 papers
Compares with it 1 paper
Cites it as background 32 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.