Phase-space integrals through Mellin-Barnes representation

Abstract: We compute angular phase-space integrals with three and four denominators analytically, working within dimensional regularisation via the Mellin-Barnes (MB) representation. The approach converts multifold MB integrals into real parametric integrals and expresses all results in terms of Goncharov polylogarithms (GPLs). For three denominators, all-massless results are obtained to $\mathcal{O}(ε^2)$ and the single-massive case to $\mathcal{O}(ε)$; for four denominators, both the massless and single-massive cases are solved to $\mathcal{O}(ε^0)$. Integrals with multiple massive momenta follow from a partial fraction decomposition reducing them to the single-massive case. Recursion relations relating integrals with higher denominator powers to master integrals are derived. These are essential ingredients to solving full phase-space integrals.

• The paper presents a systematic Mellin-Barnes method that translates multi-fold integrals into explicit Goncharov polylogarithm representations for phase-space integrals in QCD.
• It details a robust ε-expansion and analytic continuation pipeline that resolves singularities and reduces complex integrals to real parametric forms.
• The method achieves significant computational speedup and scalability, enabling high-precision analytic evaluations of NNLO real-emission contributions.

Mellin-Barnes Representation of Multi-Denominator Phase-Space Integrals

Introduction and Context

The analytic computation of real-emission phase-space (PS) integrals is a critical yet comparatively underdeveloped aspect of perturbative QCD, especially at higher orders where complex multi-loop virtual corrections and real-emission contributions must be treated simultaneously. This paper systematically develops a Mellin-Barnes (MB) based framework for evaluating angular PS integrals with three and four denominators in the context of dimensional regularization, focusing on the full infrared singularity structure and finite terms that enter physical cross sections.

The angular component of PS integrals, when properly isolated, carries universal collinear/infrared pole structure and is fully expressible in terms of the reference momenta. For small numbers of propagator denominators ($n=1,2$), closed analytic forms in terms of hypergeometric or Lauricella functions are standard; however, for $n\geq 3$ the all-orders $\epsilon$ expansion is structurally complex and was previously unsolved in generality. This work leverages the MB representation to systematically translate these multi-fold contour integrals into real parametric integrals, producing results in terms of Goncharov polylogarithms (GPLs), which are amenable to further analytic manipulation.

Mellin-Barnes Pipeline and Reduction Strategy

The authors construct MB representations for angular integrals, explicitly encoding kinematic dependencies in products of $\Gamma$-functions over integration variables associated with the scalar products of external momenta. The expansion and evaluation pipeline comprises:

1. Analytic Continuation: Systematic tracking and resolution of contour-crossing poles as $\epsilon\to 0$ ensures a well-defined Laurent expansion. Residue extraction and iterative pole handling reduce MB integrals to lower-fold cases where necessary.
2. $\epsilon$-Expansion: Following analytic continuation, the integrand is safely expanded in $\epsilon$. Each term in the Laurent series corresponds to simplified MB integrals.
3. MB to Real Parametric Integrals: Balancing $\Gamma$-function products enables reduction to Euler beta functions, leading to real integrals via delta-function constraints.
4. Expression in GPLs: The resulting integrands are integrated iteratively, yielding GPLs with alphabets determined by the kinematic structure; rationalization techniques manage square roots when necessary.

This pipeline is fully algorithmic and scalable to large numbers of denominators, with computational complexity as the primary bottleneck.

Results: Three-Denominator Integrals

All-Massless Case

For $n=3$ with massless external momenta, the MB approach yields explicit GPL representations up to $\mathcal{O}(\epsilon^2)$, with three scalar invariants and a threefold MB integral. The collinear singularity structure is confirmed: a single $n\geq 3$0 pole emerges, and no double pole is present, consistent with expectations from phase-space geometry. The absence of square roots in the integration variables strongly facilitates iterated integration, and the expressions obtained are directly suitable for convolution with process-dependent radial components.

Single-Massive Case

The addition of a massive external momentum increases the MB integral to fourfold. Eleven distinct MB integrals arise at $n\geq 3$1, ten of which are standard, while the eleventh involves a quadratic radicand requiring rationalization. The resulting expressions are fully converted to GPLs, with a systematic algorithm for extension to additional massive legs.

Multiple Massive Legs and Recursion Relations

Partial fraction decomposition, exploiting linearity in the coordinate representation, enables systematic reduction of multi-massive integrals to sums over single-massive configurations. Recursion relations generated by kinematic differentiation provide full reduction to a master integral basis, analogous to IBP identities in multi-loop Feynman integrals. This framework greatly simplifies the computation and reuse of integrated results in higher-order calculations.

Combination with Radial Part and Singularities

The convolution with the radial component is enabled by the closed GPL form. For massless integrals, soft singularities appear in individual $n\geq 3$2 coefficients, but these are shown to factorize and exponentiate, allowing an explicit subtraction scheme that renders the angular component safely expandable order by order.

Results: Four-Denominator Integrals

All-Massless and Single-Massive Cases

The MB approach generalizes to $n\geq 3$3, yielding six- and seven-fold integrals for massless and single-massive configurations, respectively. The $n\geq 3$4 pole in the massless case is manifestly $n\geq 3$5 symmetric, determined entirely by collinear singularity, while the finite part is expressed in GPLs of weight 2 with an alphabet including 17 letters (rational and square-root valued). These letters encode complicated branch point structures analogous to Gram determinants in multi-loop computations.

In the single-massive case, symmetry reduction is observed, and the expansion produces a distinct set of 11 GPL letters. The computational efficiency improvement is stark: numerical evaluation of the GPL expressions is $n\geq 3$61800 times faster than direct MB numerical integration.

Extension to Multi-Massive and Higher Denominator Cases

Extensions to two, three, and four massive external momenta follow by iterated PF decomposition, with explicit expressions provided. Recursion relations for shifting propagator powers are generated analogously to the $n\geq 3$7 case, allowing for efficient reduction.

Algorithmic Scalability

The MB/GPL method is fully algorithmic with no conceptual limit on the number of denominators, only practical computational challenges as the MB fold increases.

Implications and Outlook

From a practical standpoint, the explicit GPL representations for angular PS integrals facilitate analytic and high-precision numerical evaluation of real-emission contributions in higher-order QCD calculations. This is particularly notable for processes such as semi-inclusive deep-inelastic scattering at NNLO and beyond, where fully analytic phase-space integration remains a significant challenge.

The theoretical significance is twofold: (1) the MB method provides a systematic reduction to real integrals, and (2) GPLs afford closure under iterated integration, a decisive benefit over representations involving Clausen functions. This enables analytic convolution with the radial component, ensuring that final cross sections may, in many cases, be obtained in explicit analytic form.

For future developments, the methodology is directly scalable to $n\geq 3$8 denominators, relevant to multi-particle final states in high-multiplicity perturbative corrections and multi-loop phenomenology. In cases where function spaces exceed GPLs (e.g., elliptic or modular forms), the MB pipeline still outputs results as iterated integrals, retaining analytic accessibility for convolutions. The computational infrastructure developed herein paves the way for precision QCD at electron-ion collider and other high-energy facilities.

Conclusion

This paper establishes a robust, analytic, and algorithmic Mellin-Barnes-based framework for evaluating angular phase-space integrals with three and four denominators in dimensional regularization, expressing all results as Laurent expansions in $n\geq 3$9 with complete GPL representations. The recursive structure, PF decomposition, and singularity management strategies developed herein constitute essential tools for the analytic computation of real-emission phase-space integrals in multi-loop QCD. The closure under iterated integration and significant computational speedup demonstrate both theoretical and practical utility, facilitating advances in precision perturbative phenomenology and providing a foundation for further generalizations in advanced scattering calculations (2604.01505).