Invariant-mass distribution of top-quark pairs and top-quark mass determination

In this study, we investigate the invariant-mass distribution of top-quark pairs near the 2mt threshold, which strongly influences the determination of the top-quark mass mt. Higher-order non-relativistic corrections lead to large contributions, which are not included in the state-of-the-art theoretical predictions. A factorization formula is derived to resum such corrections to all orders in the strong-coupling, and necessary ingredients are calculated to perform the resummation at next-to-leading power. We combine the resummation with fixed-order results and present phenomenologically relevant numerical results. The resummation effect significantly increases the differential cross-section in the threshold region and makes the theoretical prediction more compatible with experimental data. We estimate that using our prediction in the determination of mt will lead to a value closer to the direct measurement result.

Top-quark pair production is among the most important scattering processes at the Large Hadron Collider (LHC).The experimental measurements of its differential cross-sections have achieved remarkably high resolutions [1][2][3][4] in the Run 2 of the LHC with a center-of-mass energy .In the meantime, the theoretical modeling of the kinematic distributions in quantum chromodynamics (QCD) was significantly improved with the emergence of the next-to-nextto-leading order (NNLO) predictions [5][6][7] supplemented with the resummation of large logarithms to the next-tonext-to-leading logarithmic accuracy (NNLL') [8][9][10].The complete next-to-leading order (NLO) corrections incorporating electroweak (EW) effects [11,12] have also been combined with the higher order QCD corrections, resulting in the state-of-the-art standard model (SM) pre-dictions for this process [13,14].Despite these developments, small discrepancies between the high precision theoretical and experimental results are persistent.A notable discrepancy is the invariant-mass distribution near the threshold, measured by the CMS collaboration at the 13 TeV LHC [2,3], where is the top-quark mass, and the invariant mass is defined as .In Fig. 1, we depict the averaged distribution in the [300, 380] GeV range.The CMS result [2] is shown as the green band, which reflects the combined statistical and systematic uncertainty.The central values of various theoretical predictions (NNLO from [6], NNLO+EW from [13], and NNLO+NNLL' from [9]) are shown for comparison.There is a clear gap between the experimental and theoretical results.The theoretical predictions in this region are strongly dependent on .Consequently, this discrepancy has a strong impact on the extraction of from kinematic distributions [15], which favors a value of that is significantly lower than the current world average.Therefore, it is of interest to investigate whether the top quark is indeed lighter than expected, or whether there are missing contributions in this region that are not incorporated in the most up-to-date theoretical predictions.
In this study, we investigate a class of higher order QCD corrections in the threshold region .These originate from exchanges of Coulomb-type gluons as well as soft gluons between the top and anti-top quarks.The leading contributions take the form , where is the velocity of the top and anti-top quarks in the rest frame.The region corresponds to , where the top and anti-top quarks are slowly moving with respect to each other.In the small region, the contributions (as well as those of ) are enhanced and lead to large corrections beyond the NNLO(+NNLL') results.These corrections are divergent as , signaling a breakdown of the perturbative expansion in this limit.At NNLO (order ), this does not pose a severe problem, as the integration over is still convergent, which yields a finite prediction for the averaged differential cross-section, as shown in Fig. 1.However, at order and beyond, even the integration across becomes divergent, and one has to perform an all-order resummation to restore the validity of theoretical predictions.Furthermore, at fixed orders in perturbation theory, due to the phase-space constraint.This is not physical, because bound-state effects can make the mass of the system smaller than , where Q is a heavy quark.For charm quarks and bottom quarks, such bound-state effects are genuinely non-perturbative.For top-quark pairs, the bound-state effects can be (to a good approximation) studied in perturbation the- t ory, as long as we integrate in a large enough range around the threshold .
In the following part of this study, we present a factorization formula for the partonic differential cross-sections that is valid in the region , which schematically takes the form .The hard function H describes exchanges and emissions of hard gluons with typical momenta of , while the potential function J resums exchanges of Coulomb and soft gluons between the top and anti-top quarks.Notably, we do not take the limit here, where is the partonic center-of-mass energy.Consequently, extra hard radiations are allowed and incorporated in our framework into the hard function H. Therefore, our factorization formula is intrinsically different from those in [16][17][18], where the approximation was employed.Our formalism is similar to that of [19,20]; however, it takes into account the kinematic dependence of the renormalization and factorization scales, and differs in the treatment of subleading corrections in .Subsequently, we present numeric results that are relevant to LHC phenomenology based on our factorization formula.We show that the resummation effects indeed enhance the distribution significantly in the threshold region, and make the theoretical prediction more compatible with the experimental data.

Factorization near threshold t t
We begin with the invariant mass distribution written as a convolution of the partonic luminosity functions and partonic differential cross-sections where denote partons within the colliding protons; , , and is the factorization scale.The parton luminosity functions are defined by where is the parton distribution function (PDF) of the parton i in the proton p. Notably, in writing down Eq. ( 1), we have taken into account that the factorization scale may depend on kinematic variables (which we collectively denote as ) other than and .This is necessary, as we combine our result with the NNLO result of [6], where the scales are correlated with the variable , where and are the transverse momenta of the top and anti-top quarks, respectively.
We are concerned with QCD corrections to the partonic differential cross-sections.At the leading order  (LO) in , only the and gg channels yield non-vanishing contributions where is the renormalization scale, and is the scattering angle of the top quark in the rest frame.The coefficients , with labelling the color configuration of the system, are given by β 3) is proportional to due to the two-body phase space.This indicates that the LO differential cross-section approaches zero at the threshold .This is no longer true at higher orders in , essentially because of the terms mentioned in the Introduction.More precisely, the NLO differential cross-section approaches a positive constant in the limit, at which the NNLO and higher order ones diverge.This behavior makes the perturbative expansion badly convergent in the threshold region.In the following, we analyze the properties of the differential cross sections and resum the large and corrections to all orders of .This is the only way to restore the predictive power of perturbative QCD near threshold.
In the threshold limit or , there are large hierarchies among the energy scales , , and .Using the method of regions, we identify the following momentum modes in the rest frame [21] hard : O(M tt ) The top and anti-top quarks move very slowly and are non-relativistic objects, whose residue momenta correspond to the potential mode.They can interact with each other through soft, ultrasoft, and potential gluons, which are described by the effective Lagrangian of potential non-relativistic QCD (pNRQCD) [22,23].The hard mode describes fluctuations with typical momenta of , resulting in Wilson coefficients of effective operators in pNRQCD.We note that there are no collinear modes in our setup.This is related to the fact that we are not taking the limit , such that the extra emissions are not constrained.Consequently, we do not need to employ the soft-collinear effective theory, as opposed to [17,18].Using the method of effective field theory (EFT), we derive a factorization formula for the partonic cross-sections in the threshold region.The power expansion of pNRQCD is according to the counting .Up to the next-to-leading power (NLP), the factorization formula reads where are hard functions describing hard gluon exchanges and emissions, which depend on the transverse momentum and the rapidity Y of the pair; (E) with are potential functions describing exchanges of potential, soft and ultrasoft gluons between the top and anti-top quarks.At NLP, ultrasoft gluon exchanges among initial-state partons and the pair do not contribute.Consequently, there is no soft function in our factorization formula.This contribution could be present at higher powers in and , which is an interesting subject of study for the future.The prefactors in Eq. ( 6) are chosen such that the LO cross-sections in the EFT reproduce the exact ones in Eq. (3).These prefactors take into account subleading corrections in at LO and have significant impacts at higher orders in .This, along with the explicit kinematic dependence of the hard function, makes our factorization formula substantially different from that in [19,20].
The hard functions are calculated order-by-order in in terms of differential cross-sections in the limit .For resummation at NLP, we must find the hard functions up to the NLO.The LO hard functions are given by 1 − z Q T We analytically calculated the NLO corrections to the hard functions, which were hereto unknown in the literature.They receive contributions from both virtual exchanges and real emissions of gluons.The virtual and real contributions are separately infrared divergent.The divergences cancel when combined with the PDF counterterms.After the cancellation, we obtain finite NLO hard functions involving singular distributions of , , and Y.The results are more complicated than those in Chinese Physics C Vol. 44, No. 9 (2020) 091001 091001-3 the literature [19,20], as we must maintain the dependencies on kinematic variables in order to use the -based renormalization and factorization scales.The potential function is related to the imaginary part of the pN-RQCD Green function of the pair at the origin.It receives contributions from potential, soft, and ultrasoft modes.The explicit expressions for the Green function up to NLP can be found in, e.g., [18,24].The finite width of the top quark is considered by the replacement of in the potential function.This treatment is valid at NLP; however, additional considerations are required at higher powers (we count ) [25][26][27][28][29].

M tt
Combining the hard and potential functions, we can now produce NLP resummed predictions for the distribution near threshold.These can be further combined with fixed-order ones via the matching formula where we used "nLO" and "nnLO" to denote fixed-order expansions of the NLP resummed result ( 6) up to the second and third order in .Their differences with respect to the full NLO and NNLO results are higherpower terms in , which we wish to consider through the matching procedure.

Numeric results
We present the numerical results based on our factorization formula (6) and the matching formula (8).Throughout our calculation, we take , , and use the NNPDF3.1 NNLO PDFs with .Because the small-resummation is valid only at low , we restrict ourselves to GeV in this work, which is the region most sensitive to .Following [9,13], we choose the default values of the renormalization scale and the factorization scale to be , the corresponding results are regarded as central values.We then vary the scales simultaneously up and down by a factor of two to estimate the remaining theoretical uncertainties.With such a procedure, the error bands are not necessarily symmetric around the central values.First, in Fig. 2 we show the exact NLO correction computed using [30], in comparison with its approximations in various kinematic limits.Such a comparison is important to assess the validity of the approximations within the phase-space region of interest.Fig. 2 has a few significant implications, which we like to elaborate on.1) The result corresponds to the second term in the fixed-order expansion of our NLP formula (6).It provides an excellent description of the full NLO correction in the range [340,380] GeV.This demonstrates that our small-approximation is indeed valid and supports the resummation performed in this work.2) The result corresponds to the "soft " or "threshold " limit considered in e.g., [8,31].It should be stressed that the concept of "soft" or "threshold" has a completely different meaning than those used in this work.We observe that, while not perfect, the limit provides a reasonable approximation to the full result in the low region.It is known that this limit works better in the high region [8,31].3) Finally, the double limit and corresponds to the "soft " limit considered in [16,20].This approximation does not capture the dominant contribution at NLO.This is essentially the reason why we do not consider such a "soft " resummation in this work.
We now turn to compare the NLP resummed result Eq. ( 6) and its fixed-order expansion.As in Eq. ( 8), we label the expansion up to the i-th term ( ) as n i LO.Here, the lower-case 'n' implies that these are approximations to the full N i LO results in the limit.The n 0 LO result equals exactly the full LO result because of the prefactors in Eq. (6).We show such a comparison in Fig. 3, from which we can draw several important conclusions.1) The fixed-order expansion converges rather quickly when is significantly larger than .However, when approaches the threshold, the behavior becomes out-of-control.In particular, the n 3 LO result tends to , while the n 4 LO tends to in the limit ( ).
2) The singularity at is regularized by the resummation effects, and we obtain a finite prediction near with the NLP resummation formula.In the resummed result, the region is allowed due to bound-state effects.The shape of the NLP curve for depends crucially on the top quark width.However, we have verified that the integrated cross-section in the [300,380] GeV range is insensitive to the value of .3) As a final implication of Fig. 3, we note that the NLP curve almost overlaps with the nLO and n 2 LO curves for .This indicates that in this region, the matched (N)NLO+NLP results of Eq. ( 8) are governed by the fixed-order (N)NLO calculations.Effectively, this shows that we do not apply the small-resummation to regions where the break down of the EFT description is a possibility.Indeed, we have checked that the dominant beyond-NNLO correction comes from the region , where and pNRQCD is perfectly applicable.Finally, we match our resummed calculation to the NLO and NNLO results according to Eq. ( 8).The NLO results are computed using and the NNLO results are obtained from [6,7,32,33].The NLO+NLP and NNLO+NLP results are shown in Fig. 4 and Table 1 along with the NLO and NNLO results, compared against the CMS measurement.We find that the resummation effects enhance the NNLO differential cross-section by about 9%, and make the theoretical prediction more compatible with experimental data.As discussed in the Introduction, this enhancement will have a significant influence on the extraction of the top-quark mass from kinematic distributions.Although we cannot repeat the analysis of [15], we can estimate the impact by studying the -dependence of the NLO (used in the fit of [15]) and NNLO+NLP predictions for the [300, 380] GeV range.For example, the central value of the NNLO+NLP differential cross-section for is given by 1.434 pb/GeV.To achieve the same differential cross-section using the NLO calculation, one must lower the value of to approximately 171 GeV.This roughly explains the outcome of [15] and shows that using our prediction in the fit will lead to a shift of about 1.5 GeV resulting in a value closer to the world average.We studied the invariant-mass distribution near the threshold.In this region, a small gap exists between the most up-to-date theoretical predictions and the experimental measurements by the CMS collaboration at the LHC.This leads to a discrepancy between the values of determined from direct measurements and from fitting kinematic distributions in production.We show that higher-order non-relativistic effects lead to large corrections to the differential cross-section in the threshold region, which were not included in the state-of-the-art theoretical predictions.We derive a factorization formula to resum such corrections to all orders in strong coupling, and calculate necessary ingredients to perform the resummation at next-to-leading power.We combine the resummation with NLO and NNLO results, and present numeric results relevant for LHC phenomenology.We find that the resummation effect increases the differential crosssection in the range by about 9%.This makes the theoretical prediction more compatible with experimental data, and leads to a shift of about 1.5 GeV on the extracted top-quark mass from kinematic distributions, resulting in a value that is more consistent with the world average.Our results can be easily combined with the soft gluon resummation of [8][9][10]31] and with electroweak corrections of [11][12][13][14].This allows a precision prediction for the invariant-mass spectrum across the entire phase space.Our framework can be extended to study the t t + jets production process, which is also used to extract the top quark mass [34,35].

L. L. Yang would like to thank A. Mitov for useful discussions.
tt

Fig. 1 .
Fig. 1. (color online) Averaged invariant mass distribution in [300, 380] GeV range.The CMS result [2] is shown as the green band, while the central-values of various theoretical predictions are shown in comparison.

Table 1 .
Differential cross-sections (in units of pb/GeV) entering