
discussed in [70, 71].
This work will investigate the projection of the t¯
tsystem’s transverse momentum with respect to a reference
unit vector ~τ on the azimuthal plane, more explicitly, qτ≡ |~qτk| ≡ |~qT·~τ|. In contrast to the traditional
transverse momentum spectrum dσt¯
t/dqT,qT=|~qT|, where both components of ~qTare measured and thereby
constrained, the present observable dσt¯
t/dqτonly concerns the projected piece ~qτk, leaving the perpendicular
part ~qτ⊥unresolved and, hence, it should be integrated out. As will be demonstrated in this paper, the act
of integrating out the perpendicular components will introduce new and distinguishing features to the qτ
spectrum, particularly in regards to the treatment on the azimuthal asymmetric contributions [66,67,72,73].
Induced by the soft and collinear radiation, the fixed-order calculation of the qτdistribution exhibits substan-
tial higher-order corrections in the small qτregion. This, thus, necessitates a resummation of the dominant
contributions in this regime to all orders to stabilise the perturbative predictions. In order to accomplish
this target, one of the prerequisite conditions is to determine the dynamic regions driving the asymptotic
behaviour in the limit qτ→0. For the classic transverse momentum resummation, this analysis was first
presented for Drell-Yan production in [74] by means of inspecting the power laws of a generic configura-
tion on the pinch singularity surface [75–77]. It was proven that the leading singular behaviour in the
small qTdomain was well captured by the beam-collinear, soft, and hard regions. However, this conclusion
cannot be straightforwardly applied onto the qτresummation in top-pair production, as the deep recoil
configuration |~qτ⊥| ∼ Mt¯
tqτ, which stems from the integral over the perpendicular component, was kine-
matically excluded in [74]. Therefore, for delivering an honest and self-consistent study on dσt¯
t/dqτ, this
work will reappraise the scalings of the relevant configurations, comprising both the qT-like configuration
Mt¯
t |~qτ⊥| ∼ qτand the asymmetric one |~qτ⊥| ∼ Mt¯
tqτ.
To this end, we will exploit the strategy of expansion by regions [78–81] to motivate the momentum modes
governing the low qτregime, which will cover the beam-collinear, soft, central-jet, and hard regions. Then,
the soft-collinear effective theory (SCET) [82–91] and the heavy quark effective theory (HQET) [92–95]
are used to embody those dynamic modes, thereby calculating the effective amplitudes and the respective
differential cross sections. After carrying out a multipole expansion, the results constructed by those dynamic
regions all reflect the unambiguous scaling behaviors, which can be determined from the power prescriptions
of the relevant effective fields. From the outcome, we point out that the leading asymptotic behavior of
dσt¯
t/dqτis still resultant of the symmetric configuration Mt¯
t |~qτ⊥| ∼ qτ, which is in practice dictated by
the beam-collinear, soft, and hard momenta, akin to the case of dσt¯
t/dqT, whereas the contributions from
the |~qτ⊥| ∼ Mt¯
tqτpattern are suppressed by at least one power of λτ≡qτ/Mt¯
t.
Upon the identification of the leading regions, we make use of the decoupling properties of the soft modes [58,
82] to derive the factorisation formula for dσt¯
t/dqτ. Owing to the integration over ~qτ⊥, the impact space
integrals herein are all reduced from 2D to 1D. Thus, the azimuthal asymmetric contributions, which in
principle contain divergent terms in the asymptotic regime in the general dσt¯
t/d~qTcross section after com-
pleting the inverse Fourier transformations, do not contribute any divergences to the qτspectrum. This
is the second ~qT-based observable free of asymmetric singularities in addition to the azimuthally averaged
spectra dσt¯
t/dqT[64, 65, 68].
To implement the resummation, we employ the renormalisation group equations (RGE) and the rapidity
renormalisation group equations (RaGE) to evolve the intrinsic scales in the respective ingredients and in turn
accomplish the logarithmic exponentiations [96–99]. Alternative approaches can also be found in [100–109].
For assessing the resummation accuracy, we take the logarithmic counting rule LM∼α−1
s∼λ−1
Lthroughout
and organise the perturbative corrections to the relevant sectors in line with the following prescription,
dσt¯
t
dqτ∼σBorn
t¯
texp LMf0(αsLM)
| {z }
(LL)
+f1(αsLM)
| {z }
(NLL,NLL0)
+αsf2(αsLM)
| {z }
(N2LL,N2LL0)
+α2
sf3(αsLM)
| {z }
(N3LL,N3LL0)
+. . .
×1(LL,NLL) + αs(NLL0,N2LL) + α2
s(N2LL0,N3LL) + α3
s(N3LL0,N4LL) + . . . .
(1.1)
Therein, the desired precisions of the anomalous dimensions are specified between the square brackets in the
exponent as for a given logarithmic accuracy, while the according requirements on the fixed-order elements
are presented within the curly brackets. In this work, we will evaluate and compare the resummed qτspectra
on the next-to-leading-logarithmic (NLL), N2LL, and approximate N2LL0(aN2LL0) levels.
The paper is structured as follows. In Sec. 2, we will utilise the strategy of expansion by dynamic re-
gions [78–81] and effective field theories, i.e. SCETII [89–91] and HQET [92–95], to derive the factorisation
3