Pith. sign in

REVIEW 4 major objections 4 minor 24 references

Angular phase-space integrals with four denominators through Mellin--Barnes

T0 review · 4 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Four-denominator angular phase-space integrals are claimed to admit closed analytic forms in Goncharov polylogarithms through the finite order in the dimensional regulator, with multi-mass cases reduced to the single-mass case by partial fr

desk verdict Massless four-denominator angular integrals in GPL form look solid, but the single-massive result is under-verified: Table 2 contradicts the 'perfect agreement' claim, so the paper needs a fix before I'd trust it. read the letter →

arxiv 2508.15952 v1 pith:WSSPL2JI submitted 2025-08-21 hep-ph hep-thmath-phmath.MP

Taushif Ahmed , Syed Mehedi Hasan , Andreas Rapakoulias This is my paper
classification hep-phhep-thmath-phmath.MP
keywords angularphase-spaceintegralsMellin–BarnesGoncharovpolylogarithmsdimensionalregularizationpartialfractiondecompositionhigher-orderQCDcorrectionsreal-emissioniterated
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

In higher-order QCD predictions, real-emission phase-space integrals factorize into a process-dependent radial part and a universal angular part. This paper targets the four-denominator angular integral, the next case beyond two- and three-denominator results, and claims that its massless and single-massive versions can be evaluated through the finite term in the dimensional regulator as explicit Goncharov polylogarithms (GPLs) — iterated logarithmic integrals that are analytically complete and numerically fast. Double-, triple-, and quartic-massive cases are then obtained by a partial-fraction decomposition that rewrites each multi-mass angular integral as a linear combination of single-mass ones. Getting this object into GPL form matters because the angular factor can then be combined with the radial factor to produce fully analytic phase-space integrals for high-order scattering processes. The paper also highlights the solution of six- and seven-fold Mellin–Barnes integrals as a first in the literature.

What carries the argument

The Mellin–Barnes representation (2.1) turns the angular integral into a multi-fold contour integral over gamma-function products. The load-bearing steps are: expansion around epsilon=0 with analytic continuation; conversion of balanced MB integrals — those with equal gamma-function exponent sums on both sides of each contour — into real integrals via Euler beta functions and delta-function identities; and an iterative GPL manipulation that uses the integral definition of GPLs to move integration variables to the rightmost argument, regularizing spurious singularities along the way. For the multi-mass cases, a two-point splitting parameter lambda_ij supplies the partial-fraction decompositio

What would settle it

Compute the single-massive angular integral at Set 1 of Table 2 (v11=0.11, v12=0.45, v13=0.40, v14=0.35, v23=0.30, v24=0.35, v34=0.40) to high precision by an independent method; if the finite term settles at -16.5783 rather than -16.4429, the GPL expression has an error at O(epsilon^0).

Watch

Extended reading notes

Core claim

The central claim is that the massless four-denominator angular integral I^(0)_{1,1,1,1} and the single-massive analogue I^(1)_{1,1,1,1} can be expressed up to O(epsilon^0) in terms of Goncharov polylogarithms, with highest transcendental weight two and a single 1/epsilon pole from a collinear singularity. The paper obtains these expressions by converting six- and seven-fold Mellin–Barnes integrals into real integrals, then rearranging the resulting iterated integrals into GPL form. It further derives explicit partial-fraction identities that express double-, triple-, and quartic-massive four-denominator angular integrals as combinations of single-massive ones. The paper presents this as the

Load-bearing premise

The load-bearing premise is that the two in-house GPL-reduction algorithms—described only schematically in Section 4—correctly handle non-linear GPL weights; Table 2's Set-1 finite-term mismatch (-16.4429 vs -16.5783) is direct evidence that premise is not yet fully verified.

Editorial extensions

If this is right

  • Massless and single-massive four-denominator angular integrals become available as closed GPL expressions, reducing evaluation at a phase-space point from about half an hour by direct MB numerics to about a second.
  • All double-, triple-, and quartic-massive configurations are linear combinations of single-mass results, so only one genuinely new integral needs to be solved for the whole four-denominator family.
  • The GPL form can be combined with the radial component of a phase-space integral, opening the way to fully analytic cross-section predictions rather than one-dimensional numerical integration.
  • The same algorithmic MB-to-GPL conversion extends beyond angular integrals, making high-fold Mellin–Barnes integrals amenable to analytic evaluation in other settings.
  • Recursion relations reduce integrals with higher denominator powers to a finite master set, an angular analogue of integration-by-parts identities.

Reading between the lines

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

  • If the algorithm scales as claimed, five-denominator angular integrals are the immediate next target; the likely obstruction is the same non-linear GPL-weight reduction that the paper only sketches.
  • The Table 2 Set-1 mismatch suggests a residual error in one of the fifteen real integrals or in the GPL reduction, so the single-mass finite term should be treated as provisional until independently confirmed.
  • The final function space being only weight-2 GPLs, with square roots entering only through the letters, hints at a kinematic structure worth mapping: it may predict which higher-denominator cases stay within the GPL algebra.
  • Releasing the reduction code would turn this from a stated result into a reproducible method that others can apply to their own Mellin–Barnes integrals.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The manuscript presents a Mellin--Barnes (MB) treatment of the four-denominator angular phase-space integral (1.1), expanding it in the dimensional regulator and reducing the resulting six-fold (massless) and seven-fold (single-massive) MB integrals to expressions in terms of Goncharov polylogarithms up to O(epsilon^0). For the double-, triple-, and quartic-massive configurations, the paper does not perform the MB reductions directly but instead extends a partial-fraction decomposition (Eqs. (6.9), (6.11), (6.12)) that relates these integrals to the single-massive result. It reports pole terms matching refs. [10,17] and claims the first GPL representation of such high-fold MB integrals.

Significance. If the results are correct, this is a useful and non-trivial contribution: explicit GPL expressions for the massless and single-massive four-denominator integrals, a new algorithmic reduction of high-fold MB integrals, and a partial-fraction route to multi-mass configurations that avoids solving even more complicated MB integrals. The massless finite part passes the numerical checks in Table 1, and the pole terms agree with known results. However, the main new quantity, the single-massive finite part, is only checked in three phase-space points, and one of those three checks fails at a level far above the displayed precision. Since the single-massive expression is supplied only as an ancillary file and the reduction algorithm of Section 4 is described only schematically, the central claim is not yet independently verifiable. The paper is therefore promising but currently falls short of full support for its main result.

major comments (4)
  1. [Section 6.1, Table 2, Eq. (6.1), ancillary file <one_mass_ep0.m>] The text states that the finite term of the single-massive integral has been verified against direct MB evaluation and that the agreement is 'perfect'. However, Table 2, Set 1 shows -16.4429 for the authors' expression versus -16.5783 for direct MB evaluation, a difference of 0.1354. This is many orders of magnitude larger than the six-digit precision reported in the table. Because this discrepancy concerns the sole direct numerical check of the new single-massive GPL result, it is load-bearing. Please correct either the expression/file or the direct evaluation, and show the numerical error/uncertainty of both results.
  2. [Section 4 and Section 6.1] The final GPL expressions are not displayed in the paper; they are only given in Mathematica-readable ancillary files. The reduction from real integrals to GPLs is described schematically ('in-house algorithms', 'we iterate this process', 'carefully regulate spurious singularities'), and no code is provided. Thus a reader cannot verify the central step of the calculation. Please provide the analytic expressions, the code, or a complete proof-of-reduction for at least one nontrivial case, including how non-linear rational weights are handled and how spurious singularities are regulated.
  3. [Section 6.2, Eqs. (6.9)-(6.12) and variable definitions after Eq. (6.12)] The partial-fraction decomposition formulas are central to the multi-mass results but they contain apparent typos, e.g. the expression for v57 reads 'lambda-bar12 lambda-bar23 v13 + lambda-bar12 lambda23 v13 + ...', i.e. the same v13 term appears twice, and v35 is written as 'lambda-bar12 v13 + lambda12^2 v33', which is dimensionally/kinematically suspicious. The paper also refers to 'explicit forms defined in (7.1)' although Eq. (7.1) is the recursion relation, not the list of new scalar products. These formulas must be corrected and verified, since the double-, triple-, and quartic-massive results inherit their correctness from these identities.
  4. [Section 6.2, Tables 1-2] No numerical check is presented for the partial-fraction identities (6.9), (6.11), or (6.12). The claims for double-, triple-, and quartic-massive cases currently rest entirely on the algebra of the decomposition and the single-massive result. At least one phase-space point for each of these configurations should be compared with direct MB evaluation, or with an independent method, before the multi-mass extension can be considered established.
minor comments (4)
  1. [Section 5, Eq. (5.4)] In the definition of I8, the factor 'v2 2 v12' is typeset ambiguously; it appears to have a missing subscript. Please clarify.
  2. [Section 6.1, after Eq. (6.6)] The text refers to 'square roots inside X^(1)_{1,2}', while Eq. (6.7) defines the quantity Z^(1)_{1,2}. Please make the notation consistent.
  3. [Section 6.2] The auxiliary vectors p5, p6, ... are introduced without an explicit statement that they are massless and without a consistency check of their Gram matrix. Because the decomposition identities are algebraic, please state the defining properties of these vectors explicitly.
  4. [General] The paper claims agreement with ref. [17] but does not show the comparison. A brief table or a statement of the numerical values would increase confidence, especially given the discrepancy in Table 2.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the GPL results are derived from an exact MB representation and numerically cross-checked; self-citations are methodological only.

full rationale

No circular dependency was found. The MB representation (Eq. 2.1) is an exact rewriting of the angular integral; the subsequent expansion, analytic continuation and residue extraction are standard operations, and the final GPL expressions for the massless and single-massive cases are checked against direct numerical evaluation of the MB integrals (Tables 1 and 2) and against published pole formulas (refs. [10,17]). The multi-mass results are obtained by explicit partial-fraction identities (Eqs. 6.9, 6.11, 6.12) that reduce to the single-massive case algebraically; no fitted parameter is renamed as a prediction. The self-citations to refs. [3,4] are for the general MB methodology and prior three-propagator application, not for the four-propagator GPL result itself, so they are not load-bearing. The numerical disagreement in Table 2 for the single-massive finite term (Set 1: -16.4429 vs -16.5783) is a validation discrepancy that the paper glosses over as 'perfect agreement', but it is a correctness risk rather than a circularity, since the analytic result is not constructed from the MB comparison values.

Assumptions & free parameters 0 free parameters · 4 assumptions · 1 invented entities

No free parameters are fitted. The central results are produced under standard dimensional-regularization and Mellin-Barnes assumptions, plus two ad hoc additions: the private GPL reduction algorithm and the auxiliary-vector construction for multi-mass decompositions. These are the least externally supported ingredients.

assumptions (4)
  • standard math The Mellin-Barnes representation (2.1) correctly represents the angular integral for Re(j_k)>0 and by analytic continuation for the cases considered.
    Used in Sections 2, 5, and 6; standard technique in the field, with v_kl != 0 assumed (massless diagonal elements omit their MB variable).
  • domain assumption After analytic continuation and expansion in epsilon, the gamma functions have positive real parts so balanced MB integrals can be rearranged into Euler beta functions and the order of integration exchanged (Sections 3, Eqs. (3.3)-(3.5)).
    Crucial for converting six- and seven-fold MB integrals to real integrals; not proven for all contours, but backed by numerical checks.
  • ad hoc to paper The new auxiliary momenta p5, p6, ... introduced by partial fraction decomposition exist as massless vectors with the scalar products defined in Section 6.2.
    Eqs. (6.9)-(6.12) and the list after (6.12) define scalar products like v15, v25, etc., but the paper does not demonstrate that a consistent set of massless vectors realizes them. Load-bearing for the multi-mass claims.
  • ad hoc to paper The GPL reduction algorithms in Section 4 correctly handle non-linear rational weights and spurious singularities.
    The final massless and single-massive expressions depend on these private algorithms; no code or complete derivation is included.
invented entities (1)
  • Auxiliary massless vectors p5...p10 in partial fraction decomposition
    purpose: Express double-, triple-, and quartic-massive four-denominator integrals as combinations of single-massive integrals.
    Mathematical bookkeeping constructed from existing scalar products; no physical interpretation or independent falsifiable handle.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Angular phase-space integrals with four denominators through Mellin--Barnes." pith.science (2026). https://pith.science/paper/WSSPL2JI

@misc{pith2026250815952,
  author       = {Pith},
  title        = {Pith review of: Angular phase-space integrals with four denominators through Mellin--Barnes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WSSPL2JI}},
  note         = {Machine review of arXiv:2508.15952}
}
read the original abstract

We compute four-denominator angular phase-space integrals using the Mellin--Barnes (MB) technique in dimensional regularisation. Independent of the scattering process, an angular integral can be categorised based on the nature of the momenta appearing in the denominators. We address all scenarios involving fully massless and massive momenta. We present a partial fraction decomposition that relates angular integrals with multiple massive momenta to those with a single massive momentum. By solving six- and seven-fold MB integrals, we express the final results up to the finite order in the dimensional regulator in terms of Goncharov polylogarithms.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

24 extracted references · 8 canonical work pages

  1. [12]

    Haug, V.A

    J. Haug, V.A. Smirnov and F. Wunder, On multi-propagator angular integrals , 2508.00693

    work page arXiv

  2. [17]

    Salvatori, The Tropical Geometry of Subtraction Schemes , 2406.14606

    G. Salvatori, The Tropical Geometry of Subtraction Schemes , 2406.14606

    arXiv

  3. [1]

    Somogyi, Angular integrals in d dimensions , J

    G. Somogyi, Angular integrals in d dimensions , J. Math. Phys. 52 (2011) 083501 [1101.3557]

    arXiv 2011

  4. [2]

    Anastasiou, C

    C. Anastasiou, C. Duhr, F. Dulat and B. Mistlberger, Soft triple-real radiation for Higgs production at N3LO , JHEP 07 (2013) 003 [ 1302.4379]

    arXiv 2013

  5. [3]

    Ahmed, S.M

    T. Ahmed, S.M. Hasan and A. Rapakoulias, Phase-space integrals through Mellin-Barnes representation, 2410.18886

    arXiv

  6. [4]

    Ahmed, S

    T. Ahmed, S. Goyal, S.M. Hasan, R.N. Lee, S.-O. Moch, V. Pathak et al., NNLO phase-space integrals for semi-inclusive deep-inelastic scattering , Phys. Rev. D 112 (2025) 014020 [2412.16509]

    arXiv 2025

  7. [5]

    Schellekens, PERTURBATIVE QCD AND LEPTON PAIR PRODUCTION (6, 1981)

    A.N. Schellekens, PERTURBATIVE QCD AND LEPTON PAIR PRODUCTION (6, 1981)

    work page 1981

  8. [6]

    van Neerven, Dimensional Regularization of Mass and Infrared Singularities in Two Loop On-shell Vertex Functions , Nucl.Phys

    W. van Neerven, Dimensional Regularization of Mass and Infrared Singularities in Two Loop On-shell Vertex Functions , Nucl.Phys. B268 (1986) 453

    work page 1986

Show all 24 references
  1. [7]

    Beenakker, H

    W. Beenakker, H. Kuijf, W.L. van Neerven and J. Smith, QCD Corrections to Heavy Quark Production in p anti-p Collisions , Phys. Rev. D 40 (1989) 54

    1989

  2. [8]

    Lyubovitskij, F

    V.E. Lyubovitskij, F. Wunder and A.S. Zhevlakov, New ideas for handling of loop and angular integrals in D-dimensions in QCD , JHEP 06 (2021) 066 [ 2102.08943]

    2021 arXiv

  3. [9]

    Wunder, Asymptotic behavior of angular integrals in the massless limit , Phys

    F. Wunder, Asymptotic behavior of angular integrals in the massless limit , Phys. Rev. D 109 (2024) 076022 [ 2403.09773]

    2024 arXiv

  4. [10]

    Smirnov and F

    V.A. Smirnov and F. Wunder, Expansion by regions meets angular integrals , JHEP 08 (2024) 138 [ 2405.13120]

    2024 arXiv

  5. [11]

    Haug and F

    J. Haug and F. Wunder, Angular integrals with three denominators via IBP, mass reduction, dimensional shift, and differential equations , JHEP 03 (2025) 141 [ 2410.18177]

    2025

  6. [13]

    Anderle, D

    D. Anderle, D. de Florian and Y. Rotstein Habarnau, Towards semi-inclusive deep inelastic scattering at next-to-next-to-leading order , Phys. Rev. D 95 (2017) 034027 [ 1612.01293]

    2017 arXiv

  7. [14]

    Wang, J.O

    B. Wang, J.O. Gonzalez-Hernandez, T.C. Rogers and N. Sato, Large Transverse Momentum in Semi-Inclusive Deeply Inelastic Scattering Beyond Lowest Order , Phys. Rev. D 99 (2019) 094029 [1903.01529]

    2019 arXiv

  8. [15]

    Vollinga and S

    J. Vollinga and S. Weinzierl, Numerical evaluation of multiple polylogarithms , Comput.Phys.Commun. 167 (2005) 177 [ hep-ph/0410259]. – 17 –

    2005 arXiv

  9. [16]

    Naterop, A

    L. Naterop, A. Signer and Y. Ulrich, handyG —Rapid numerical evaluation of generalised polylogarithms in Fortran, Comput. Phys. Commun. 253 (2020) 107165 [ 1909.01656]

    2020 arXiv

  10. [18]

    Smirnov, Analytical result for dimensionally regularized massless on shell double box , Phys

    V.A. Smirnov, Analytical result for dimensionally regularized massless on shell double box , Phys. Lett. B 460 (1999) 397 [ hep-ph/9905323]

    1999 arXiv

  11. [19]

    Tausk, Nonplanar massless two loop Feynman diagrams with four on-shell legs , Phys.Lett

    J. Tausk, Nonplanar massless two loop Feynman diagrams with four on-shell legs , Phys.Lett. B469 (1999) 225 [ hep-ph/9909506]

    1999 arXiv

  12. [20]

    Czakon, Automatized analytic continuation of Mellin-Barnes integrals , Comput

    M. Czakon, Automatized analytic continuation of Mellin-Barnes integrals , Comput. Phys. Commun. 175 (2006) 559 [ hep-ph/0511200]

    2006 arXiv

  13. [21]

    Duhr and F

    C. Duhr and F. Dulat, PolyLogTools — polylogs for the masses , JHEP 08 (2019) 135 [1904.07279]

    2019 arXiv

  14. [22]

    Duhr, Mathematical aspects of scattering amplitudes , 1411.7538

    C. Duhr, Mathematical aspects of scattering amplitudes , 1411.7538

    arXiv

  15. [23]

    Chetyrkin and F.V

    K.G. Chetyrkin and F.V. Tkachov, Integration by parts: The algorithm to calculate β-functions in 4 loops , Nucl. Phys. B 192 (1981) 159

    1981

  16. [24]

    Laporta, High precision calculation of multiloop Feynman integrals by difference equations , Int

    S. Laporta, High precision calculation of multiloop Feynman integrals by difference equations , Int. J. Mod. Phys. A 15 (2000) 5087 [ hep-ph/0102033]. – 18 –

    2000 arXiv

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.