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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:hep-ph/0506286 by the authors.
References
How later papers used it
49 later papers cite it
1 later paper by other authors is listed below: 1 used its method or tool.
Uses its method or tool 1 of 2 papers with this label
-
Used the properties of the Bessel function, citing this paper's appendix, to express a derivative in terms of the Bessel function of the third order.
New Renormalization Group Equations and the Naturalness Problem (2009) · 0.95 · paraphrased: the licence does not allow quoting
How the literature uses it
48 arXiv papers cite this in 60 sentences
Papers citing it · 24 sentences without a confident label
Uses its method or tool 2 papers
-
Generally, the discontinuity of $(t-z)^-a_1$ is [ref] [equation] Note that for $a_1=n-ϵ$, where $n$ is a positive integer, the expansion [ref] [equation] and its differential forms are used to calculate the discontinuities of the above integrals.
Analytic calculation of doubly heavy hadron spectral density in coordinate space · Integral representation of hypergeometric function · 0.81
-
Using the properties of the Bessel function (see, e.g., the Appendix in [ref] ) the derivative can be expressed in terms of the Bessel function of the third order.
New Renormalization Group Equations and the Naturalness Problem · The Case of $ϕ^4$ · 0.95
Compares with it 1 paper
-
One can prove that the results of Eq. [ref] numerically agrees with the Bessel representation [ref] for $q^2 < m_0⋯L^2$ : [equation] where $J_0(x)$ is the first kind Bessel function and $K_0(x)$ is the second kind modified Bessel function.
Discontinuities of banana integrals in dispersion relation representation · Introduction · 0.95
Cites it as background 32 papers
-
These are called banana loops and have been extensively studied in the amplitudes literature [ref] ; in the cosmological context, only a handful of analytic expressions exist for the special case of de Sitter space [ref] .
Bananas are Unparticles: Differential Equations and Cosmological Bootstrap · Introduction · 0.98
-
These integrals have long been known to involve elliptic curves [ref] .
Cuts and Isogenies · Introduction · 0.94
-
For two-dimensional QFT, we need Bessel functions $J_0$ and $Y_0$, as well as modified Bessel functions $I_0$ and $K_0$, to define propagators and compute Feynman integrals [ref] .
Wick rotations, Eichler integrals, and multi-loop Feynman diagrams · Background and motivations · 0.97
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.