Two-loop diagrams in non-relativistic QCD with elliptics
B. A. Kniehl, A. V. Kotikov, A. I Onishchenko, O. L. Veretin
TL;DR
This work analyzes two-loop NRQCD-relevant diagrams with elliptic structure arising from two masses $m$ and $M$ under threshold kinematics. It develops a master-parameter approach that reduces two-loop diagrams to integrals over one-loop building blocks with an $s$-dependent effective mass, and provides three distinct representations: a mass-ratio expansion in $x= m^2/M^2$, integral representations, and generalized hypergeometric function expressions, valid to ${\rm O}(\varepsilon)$ in $d=4-2\varepsilon$. In the equal-mass limit, the results are recast in terms of elliptic sums and elliptic integrals, revealing the elliptic nature of these diagrams beyond polylogarithms. The findings deliver practical analytic tools for NRQCD matching, parapositronium decays, and near-threshold $t\bar t$ production, and establish a foundation for further exploration of elliptic polylogarithms and related special functions in multi-loop computations.
Abstract
In this paper we consider two-loop two-, three- and four-point diagrams with elliptic structure in the case of two different masses $m$ and $M$. The latter diagrams generally arise within NRQCD matching procedures and are relevant for parapositronium decay and top pair production at threshold. We present the obtained results in several different representations: series solution with binomial coefficients, integral representation and representation in terms of generalized hypergeometric functions. The results are valid up to $\mathcal{O}(\varepsilon)$ terms in $d=4-2\varepsilon$ space-time dimensions. In the limit of equal masses $m=M$ the obtained results are written in terms of elliptic constants with explicit series representation.
