Two-loop $gg\to AA$ amplitudes near threshold and the double parton scattering singularity

Abstract

The two-loop amplitudes for pseudoscalar pair production in gluon fusion, $gg\to AA$, with each pseudoscalar emitted from its own effective gluon vertex, are known in closed form, but large cancellations near threshold obscure their analytic structure. We expand the IR-finite amplitudes in the pseudoscalar velocity ${\beta }_A$ through ${\beta }_A^6$ at fixed scattering angle. All inverse powers of ${\beta }_A$ and the odd powers through ${\beta }_A^5$ cancel, and at leading power the angular dependence enters only through $\ln (2p_T/m_A)$. We identify this logarithm with the double parton scattering (DPS) singularity: with the overall constant fixed at threshold, the factorisation formula of Gaunt and Stirling reproduces the velocity dependence of its one-loop coefficient in the forward limit to $3\times 10^{-15}$, and a one-loop $gg\to HH$ analysis confirms the corresponding CP-even weight. At two loops the equal-helicity amplitude grows as ${\ln }^3{p}_T$, with the leading logarithmic ratio of planar $\mathcal{N}=4$ super Yang-Mills theory, while the opposite helicity amplitude, which has no one-loop DPS logarithm, grows as ${\ln }^2{p}_T$. In the finite virtual correction, relative to the squared tree amplitude, the highest logarithm at order ${\alpha }_s^2$ is $-\frac{7}{64}{\alpha }_s^2{N}_c^2{\ln }^2(2p_T/m_A)$. For ${\beta }_A\le 0.2$ the series reproduce the unexpanded amplitudes with relative errors of at most $6.77\times 10^{-5}$ at the tested points; they provide this two-loop contribution to the hard function for soft-gluon resummation of the pair invariant-mass distribution near $m_{AA}=2m_A$ and are released with code for their numerical evaluation.

arXiv comment39 pages, 6 figures, 8 tables; the threshold series, code for their numerical evaluation, construction code and validation records are provided as ancillary files