A precise determination of the top-quark pole mass

TL;DR

This paper uses the Principle of Maximum Conformality to accurately determine the top-quark pole mass by comparing theoretical predictions with collider measurements, confirming QCD consistency.

Contribution

It applies PMC scale-setting to precisely extract the top-quark pole mass from collider data, improving the accuracy and consistency of QCD predictions.

Findings

Top-quark pole mass determined as approximately 174 GeV.

Predictions agree with direct measurement averages.

Validates the use of PMC in high-energy QCD analyses.

Abstract

The Principle of Maximum Conformality (PMC) provides a systematic way to eliminate the renormalization scheme and renormalization scale uncertainties for high-energy processes. We have observed that by applying PMC scale-setting, one obtains comprehensive and self-consistent pQCD predictions for the top-quark pair total cross-section and the top-quark pair forward-backward asymmetry in agreement with the measurements at the Tevatron and LHC. As a step forward, in the present paper, we determine the top-quark pole mass via a detailed comparison of the top-quark pair cross-section with the measurements at the Tevatron and LHC. The results for the top-quark pole mass are $m_t=174.6^{+3.1}_{-3.2}$ GeV for the Tevatron with $\sqrt{S}=1.96$ TeV, $m_t=173.7\pm1.5$ GeV and $174.2\pm1.7$ GeV for the LHC with $\sqrt{S} = 7$ TeV and $8$ TeV, respectively. Those predictions agree with the average, $173.34\pm0.76$ GeV, obtained from various collaborations via direct measurements. The consistency of the pQCD predictions using the PMC with all of the collider measurements at different energies provides an important verification of QCD.

Tables

Table 1. Table 1: The NNLO top-quark pair production cross-sections for the Tevatron and LHC (in unit: pb), comparing conventional versus PMC scale settings. Here all production channels have been summed up. Three typical choices for the initial renormalization scales μ r = m t / 2 subscript 𝜇 𝑟 subscript 𝑚 𝑡 2 \mu_{r}=m_{t}/2 , m t subscript 𝑚 𝑡 m_{t} and 2 ​ m t 2 subscript 𝑚 𝑡 2m_{t} have been adopted.

	Conventional			PMC
${\mu }_r$	$m_t/2$	$m_t$	$2m_t$	$m_t/2$	$m_t$	$2m_t$
${\sigma }_{\mathrm{Tevatron}}^{1.96\mathrm{TeV}}$	7.54	7.29	7.01	7.43	7.43	7.43
${\sigma }_{\mathrm{LHC}}^{7TeV}$	172.07	167.67	160.46	174.97	174.98	174.99
${\sigma }_{\mathrm{LHC}}^{8TeV}$	244.87	239.03	228.94	249.16	249.18	249.19
${\sigma }_{\mathrm{LHC}}^{13TeV}$	792.36	777.72	746.92	807.80	807.83	807.86

Table 2. Table 2: The coefficients c 0 , 1 , 2 , 3 subscript 𝑐 0 1 2 3 c_{0,1,2,3} as determined from the PMC predictions for the top-quark pair cross-section by varying the top-quark pole mass from 160 GeV to 190 GeV. The notation [ σ th ​ ( m t ) + Δ ​ σ th + ​ ( m t ) ] delimited-[] subscript 𝜎 th subscript 𝑚 𝑡 Δ subscript superscript 𝜎 th subscript 𝑚 𝑡 [\sigma_{\rm th}(m_{t})+\Delta\sigma^{+}_{\rm th}(m_{t})] indicates that the coefficients are determined using the maximum cross section within its allowable parameter range, and [ σ th ​ ( m t ) − Δ ​ σ th − ​ ( m t ) ] delimited-[] subscript 𝜎 th subscript 𝑚 𝑡 Δ subscript superscript 𝜎 th subscript 𝑚 𝑡 [\sigma_{\rm th}(m_{t})-\Delta\sigma^{-}_{\rm th}(m_{t})] corresponds to the minimum cross section.

		$c_0\left[\mathrm{pb}\right]$	$c_1\left[\mathrm{pb}\right]$	$c_2\left[\mathrm{pb}\right]$	$c_3\left[\mathrm{pb}\right]$
	${\sigma }_{\mathrm{th}}$	$7.6181$	$-0.06140$	$2.1135\times {10}^{-4}$	$-1.9319\times {10}^{-6}$
Tevatron	${\sigma }_{\mathrm{th}}+\Delta {\sigma }_{\mathrm{th}}^{+}$	$7.7580$	$-0.06261$	$2.1711\times {10}^{-4}$	$-1.9923\times {10}^{-6}$
	${\sigma }_{\mathrm{th}}-\Delta {\sigma }_{\mathrm{th}}^{-}$	$7.4796$	$-0.06019$	$2.0572\times {10}^{-4}$	$-1.8750\times {10}^{-6}$
	${\sigma }_{\mathrm{th}}$	$179.2422$	$-1.2311$	$4.7155\times {10}^{-3}$	$-3.3920\times {10}^{-5}$
LHC${|}_{7TeV}$	${\sigma }_{\mathrm{th}}+\Delta {\sigma }_{\mathrm{th}}^{+}$	$182.8195$	$-1.2590$	$4.8479\times {10}^{-3}$	$-3.5338\times {10}^{-5}$
	${\sigma }_{\mathrm{th}}-\Delta {\sigma }_{\mathrm{th}}^{-}$	$175.7093$	$-1.2037$	$4.5866\times {10}^{-3}$	$-3.2489\times {10}^{-5}$
	${\sigma }_{\mathrm{th}}$	$255.0975$	$-1.5718$	$5.2644\times {10}^{-3}$	$-4.2394\times {10}^{-5}$
LHC${|}_{8TeV}$	${\sigma }_{\mathrm{th}}+\Delta {\sigma }_{\mathrm{th}}^{+}$	$260.1779$	$-1.6078$	$5.4191\times {10}^{-3}$	$-4.4240\times {10}^{-5}$
	${\sigma }_{\mathrm{th}}-\Delta {\sigma }_{\mathrm{th}}^{-}$	$250.0801$	$-1.5364$	$5.1128\times {10}^{-3}$	$-4.0618\times {10}^{-5}$
	${\sigma }_{\mathrm{th}}$	$825.5955$	$-3.2873$	$5.1997\times {10}^{-3}$	$-1.0274\times {10}^{-4}$
LHC${|}_{13TeV}$	${\sigma }_{\mathrm{th}}+\Delta {\sigma }_{\mathrm{th}}^{+}$	$841.9260$	$-3.3675$	$5.4202\times {10}^{-3}$	$-1.0730\times {10}^{-4}$
	${\sigma }_{\mathrm{th}}-\Delta {\sigma }_{\mathrm{th}}^{-}$	$809.4638$	$-3.2084$	$4.9863\times {10}^{-3}$	$-9.8738\times {10}^{-5}$

Table 3. Table 3: Top-quark pole mass (in unit GeV) determined by D0 collaboration  Abazov:2009si ; Abazov:2011cq ; Abazov:2011mi ; Abazov:2016ekt , where the theoretical predictions for top-quark pair production is based on conventional (Conv.) scale-setting. Our predictions using PMC scale-setting are presented as a comparison.

	dilepton		lept.+jets	lept.+jets+dilepton
Conv.	${171.5}_{-8.8}^{+9.9}$  Abazov:2009si	$171.6\pm 4.3$  Abazov:2011cq ; Beneke:2012wb	${166.7}_{-4.5}^{+5.2}$  Abazov:2011mi ; Abazov:2011pta	${172.8}_{-3.2}^{+3.4}$  Abazov:2016ekt
PMC	${174.0}_{-9.8}^{+8.5}$	${172.7}_{-4.3}^{+4.1}$	$171.1\pm 4.9$	${174.6}_{-3.2}^{+3.1}$

Table 4. Table 4: Top-quark pole mass (in unit GeV) determined by CMS and ATLAS collaborations at S = 7 𝑆 7 \sqrt{S}=7 TeV  Chatrchyan:2012bra ; Chatrchyan:2013haa ; Aad:2014kva ; Khachatryan:2016mqs , where the theoretical predictions for top-quark pair production is based on conventional (Conv.) scale-setting. Our predictions using PMC scale-setting are presented as a comparison.

	dilepton	dilepton-e$\mu$
Conv.	${177.0}_{-3.3}^{+3.6}$ Chatrchyan:2012bra ; Chatrchyan:2013haa	$171.4\pm 2.6$ Aad:2014kva	${174.1}_{-2.4}^{+2.2}$ Khachatryan:2016mqs
PMC	$177.5\pm 2.4$	$171.8\pm 1.6$	$173.7\pm 1.5$

Table 5. Table 5: Top-quark pole mass (in unit GeV) determined by CMS and ATLAS collaborations at S = 8 𝑆 8 \sqrt{S}=8 TeV  Aad:2014kva ; Khachatryan:2016mqs , where the theoretical predictions for top-quark pair production is based on conventional (Conv.) scale-setting. Our predictions using PMC scale-setting are presented as a comparison.

	dilepton-e$\mu$
Conv.	$174.1\pm 2.6$ Aad:2014kva	${174.6}_{-2.5}^{+2.3}$ Khachatryan:2016mqs
PMC	$174.3\pm 1.7$	$174.2\pm 1.7$

Equations

$${\sigma }_{H_1{H}_2\to t\bar{t}X}$$

$${\mathcal{L}}_{ij}(s,S,{\mu }_f)=\frac{1}{S}\underset{s}{\overset{S}{\int }}\frac{d\hat{s}}{\hat{s}}{f}_{i/H_1}\left(x_1,{\mu }_f\right)f_{j/H_2}\left(x_2,{\mu }_f\right),$$

$${\hat{\sigma }}_{ij}$$

$$f(m_t)$$

$$f_{\mathrm{th}}(\sigma |m_t)=\frac{1}{\sqrt{2\pi }\Delta {\sigma }_{\mathrm{th}}(m_t)}\exp \left[-\frac{{\left(\sigma -{\sigma }_{\mathrm{th}}(m_t)\right)}^2}{2\Delta {\sigma }_{\mathrm{th}}^2(m_t)}\right].$$

$${\sigma }_{\mathrm{th}}(m_t)$$

$$f_{\exp }(\sigma |m_t)=\frac{1}{\sqrt{2\pi }\Delta {\sigma }_{\exp }(m_t)}\exp \left[-\frac{{\left(\sigma -{\sigma }_{\exp }(m_t)\right)}^2}{2\Delta {\sigma }_{\exp }^2(m_t)}\right],$$

$$m_t{|}_{\mathrm{Tevatron},\sqrt{S}=1.96\mathrm{TeV}}$$

$$m_t{|}_{\mathrm{LHC},\sqrt{S}=7\mathrm{TeV}}$$

$$m_t{|}_{\mathrm{LHC},\sqrt{S}=8\mathrm{TeV}}$$

$$m_t^{\overline{\mathrm{MS}}}(m_t){|}_{\mathrm{Tevatron},\sqrt{S}=1.96\mathrm{TeV}}$$

$$m_t^{\overline{\mathrm{MS}}}(m_t){|}_{\mathrm{LHC},\sqrt{S}=7\mathrm{TeV}}$$

$$m_t^{\overline{\mathrm{MS}}}(m_t){|}_{\mathrm{LHC},\sqrt{S}=8\mathrm{TeV}}$$

$$m_t=174.0\pm 1.1\mathrm{GeV}\mathrm{and}{m}_t^{\overline{\mathrm{MS}}}(m_t)=163.4\pm 1.0\mathrm{GeV}.$$

Full text

I Introduction

The top-quark is the heaviest particle of the Standard Model (SM), and its mass is one of the fundamental parameters within the SM. The large top-quark mass implies a strong top-quark Yukawa coupling to the Higgs boson, playing a special role in testing the electroweak symmetry breaking mechanism and for the search of new physics beyond the SM. The top quark decays before hadronization, and one can determine its mass by directly measuring its decay products Dalitz:1992np . Such measurements allow for the direct extraction of the top-quark mass ($m_{t}$), which however, relies heavily on the detailed reconstruction of the kinematics and reconstruction efficiency Chatrchyan:2013haa ; Aaboud:2016000 . In 2014, a combination of measurements of the top-quark mass performed by the CDF and D0 experiments at the Tevatron collider and the ATLAS and CMS experiments at the Large Hadron Collider (LHC) gives ATLAS:2014wva , $m_{t}=173.34\pm 0.76$ GeV. The direct measurements are based on analysis techniques which use top-pair events provided by Monte Carlo (MC) simulation for different assumed values of the top quark mass. Applying those techniques to data yields a mass quantity corresponding to the top quark mass scheme implemented in the MC, thus it is referred as the “MC mass”.

Another important approach for extracting the top-quark mass is done by using detailed comparisons of the pQCD predictions with the corresponding measurements; this method is indirect, but it provides complementary information on the top quark compared to direct measurements. Theoretical arguments suggest that the top-quark MC mass is within $\sim 1$ GeV of its pole mass Buckley:2011ms , and thus its use has a negligible effect on the determination of pole mass Khachatryan:2016mqs ; Fleming:2007qr 111The position of the pole in the quark propagator is defined as its pole mass, and the on-shell quark propagator has no infrared divergences in perturbation theory, it thus provides a perturbative definition of the quark mass Tarrach:1980up ; Kronfeld:1998di .. Thus in our present calculations, we shall only extract the top-quark pole mass and as usual directly take the determined top-quark MC mass by the experimental groups as the value of the top-quark pole mass. Recently, such indirect extractions of $m_{t}$ from the top-quark pair production channels by various experimental collaborations have been performed, giving the pole value, $m_{t}=173.8^{+1.7}_{-1.8}$ GeV from CMS  Khachatryan:2016mqs , $m_{t}=172.8^{+3.4}_{-3.2}$ GeV from D0 Abazov:2016ekt , and $m_{t}=172.9^{+2.5}_{-2.6}$ GeV from ATLAS Aad:2014kva .

A key goal for the indirect determinations is to have a precise theoretical prediction for the top-quark pair production cross-section in order to provide maximal constraints on $m_{t}$. Practically, one can first set an arbitrary initial renormalization scale to do the pQCD calculation, whose value only need to be large enough to ensure the applicability of the pQCD theory. Under conventional scale-setting, the renormalization scale is fixed to its initial value, which is usually chosen as the typical momentum flow of the process or the one to eliminate large logs in the perturbative series. More explicitly, it is conventional to take the renormalization scale in those predictions as the top-quark mass $m_{t}$ to eliminate the large logarithmic terms such as $\ln(\mu_{r}/m_{t})$; one then varies the renormalization scale over an arbitrary range such as $[m_{t}/2,2m_{t}]$ to ascertain the uncertainty. At sufficiently high order, a small renormalization scale-dependent prediction may be achieved for global quantities such as the total cross-section. However, such small renormalization scale dependence of the resulting prediction is due to cancelations among different orders; the renormalization scale uncertainty for each order is still uncertain and could be very large. In fact, when one applies conventional scale-setting, the renormalization scheme- and initial renormalization scale- dependence is introduced at any fixed order.

The Principle of Maximum Conformality (PMC) Wu:2013ei ; Brodsky:2011ta ; Brodsky:2011ig ; Mojaza:2012mf ; Brodsky:2013vpa provides a systematic way to eliminate renormalization scheme-and-scale ambiguities. It generalizes the BLM scale setting procedure Brodsky:1982gc to all orders. As in QED GellMann:1954fq , one shifts the argument of the running coupling at each order in the pQCD series to absorb all occurrences of the $\beta$-function. In addition, a convergent pQCD series without factorial renormalon divergence can be obtained. The PMC predictions are renormalization-scheme independent at each order in $\alpha_{s}$, since all of the scheme-dependent $\{\beta_{i}\}$-terms in the QCD perturbative series have been resummed into the running couplings. The PMC satisfies renormalization group invariance and satisfies all of the self-consistency conditions of the renormalization group Brodsky:2012ms ; Wu:2014iba , and it reduces in the $N_{C}\to 0$ Abelian limit Brodsky:1997jk to the standard Gell-Mann-Low method GellMann:1954fq . A number of PMC applications are summarized in the review Wu:2015rga ; in each case the PMC works successfully and leads to improved agreement with experiment.

By applying PMC scale-setting, we can achieve optimal renormalization scales of the process and thus obtain precise predictions for the top-quark pair production cross-section without conventional renormalization scale uncertainty Brodsky:2012rj ; Brodsky:2012sz ; Brodsky:2012ik ; Wang:2014sua . Because of the uncalculated high-order terms, there is residual renormalization scale dependence for the PMC prediction. However such residual renormalization scale dependence is generally small either due to the perturbative nature of the PMC scales or due to the fast convergence of the conformal pQCD series; e.g. we have found that the residual renormalization for top-pair production is negligibly small at the NNLO level. The PMC predictions for the top-quark pair forward-backward asymmetry are also in agreement with the corresponding CDF and D0 measurements Wang:2015lna , since it correctly assigns different renormalization scales in the one- and two- gluon exchange amplitudes.

In subsequent sections, we will determine the top-quark pole mass from a detailed comparison of the top-quark pair production cross-section predicted by applying the PMC with the measured values obtained by the Tevatron and LHC experiments.

II Top-quark pair production at the hadron colliders and the determination of top-quark pole mass

The hadronic cross-section for the top-quark pair production can be written as the convolution of the factorized partonic cross-section $\hat{\sigma}_{ij}$ with the parton luminosities ${\cal L}_{ij}$:

[TABLE]

where

[TABLE]

$x_{1}={\hat{s}/S}$ and $x_{2}={s/\hat{s}}$. Here $S$ denotes the hadronic center-of-mass energy squared, and $s=x_{1}x_{2}S$ is the subprocess center-of-mass energy squared. The parameter $\mu_{r}$ denotes the (initial) renormalization scale and $\mu_{f}$ denotes the factorization scale. The choice of $\mu_{r}$ is arbitrary, which is only need to be in pQCD region ($\gg\Lambda_{\rm QCD}$) and usually people set its value as the typical momentum flow of the process; and for this process, $\mu_{r}$ is usually chosen as $m_{t}$. The function $f_{i/H_{\alpha}}(x_{\alpha},\mu_{f})$ ($\alpha=1$ or $2$) describes the probability of finding a parton of type $i$ with a light-front momentum fraction between $x_{\alpha}$ and $x_{\alpha}+dx_{\alpha}$ in the proton $H_{\alpha}$.

The partonic subprocess cross-section $\hat{\sigma}_{ij}$ up to NNLO level can be expanded as a power series of $\alpha_{s}$:

[TABLE]

where $\rho=4m_{t}^{2}/s$. In the literature, the perturbative coefficients up to NNLO level have been calculated by various groups, e.g. Refs.Nason:1987xz ; Nason:1989zy ; Beenakker:1988bq ; Beenakker:1990maa ; Czakon:2008ii ; Moch:2008qy ; Beneke:2011mq ; Kidonakis:2010dk ; Baernreuther:2012ws ; Czakon:2012pz ; Czakon:2013goa . The LO, NLO and NNLO coefficients $f^{0}_{ij}$, $f^{1}_{ij}$ and $f^{2}_{ij}$ can be explicitly read from the HATHOR program Aliev:2010zk and the Top++ program Czakon:2011xx , where $(ij)=\{(q\bar{q}),(gg),(gq),(g\bar{q})\}$ stands for the four production channels, respectively. By carefully identifying the $n_{f}$-terms specifically associated with the $\{\beta_{i}\}$-terms in $f^{0}_{ij}$, $f^{1}_{ij}$ and $f^{2}_{ij}$, and by using the degeneracy pattern of the renormalization group equation in a recursive way, one can determine the $\beta$ terms and thus the correct arguments of the strong couplings at each perturbative order. The remaining $n_{f}$ terms arise from quark loop contributions which are ultraviolet finite. A detailed determination of the PMC scales for $\hat{\sigma}_{ij}$ up to NNLO level, including a careful treatment of the separate renormalization scales of the Coulomb-type rescattering corrections appearing in the threshold region, have been presented in Refs.Brodsky:2012rj ; Brodsky:2012sz . We shall not repeat these formulae here; the interested readers may turn to those two references for details.

In doing the numerical analysis, we will first take the top-quark pole mass as $m_{t}=173.3$ GeV toppole and choose the parton distribution functions (PDF) as the CT14 version of the CTEQ collaboration Dulat:2015mca . The NNLO $\alpha_{s}$-running is adopted with its normalization fixed in $\overline{\rm MS}$-scheme using $\alpha_{s}(M_{Z})=0.118$.

The setting of the factorization scale $\mu_{f}$ is a separate, important issue 222We have found that the factorization scale dependence is suppressed after applying the PMC Wang:2014sua ; Wang:2016wgw ; this can be explained by the fact that the pQCD series behaves much better after applying the PMC.; however, a possible determination can be based on the light-front holographic QCD Brodsky:2014yha . It determines a scale $Q_{0}$ at the interface between nonpertubative and perturbative QCD. In the analysis given here, we will take $\mu_{f}=m_{t}$.

We present the NNLO top-quark pair production cross-section at the hadronic colliders Tevatron and LHC for both conventional and PMC scale settings in Table 1, where three typical initial renormalization scales are adopted. The results shown in Table 1 show that if one uses conventional scale-setting, the renormalization scale dependence of the NNLO cross-section is still about $6\%-7\%$ for $\mu_{r}\in[m_{t}/2,2m_{t}]$. If one analyzes the pQCD series in detail, one finds that the dependence of the NNLO cross-section on the guess of the renormalization scale using conventional scale-setting is due to cancelations among different orders, and the renormalization scale dependence of each perturbative term is rather large Wang:2015lna . Thus computing a finite number of additional higher-order terms could soften the renormalization scale dependence for the total cross-section to a certain degree, but it does not eliminate the dependence on the choice of the initial renormalization scale, especially when the detailed dependence on the renormalization scales at each order is also important.

When PMC scale-setting is used, the renormalization scales are fixed by using the renormalization group equation recursively, thus fixing the arguments of the strong couplings at each order. There is residual renormalization scale dependence due to unknown NNNLO and higher-order contributions, for example, we need to known the $\beta$-terms at the NNNLO level to fix the PMC scale of the NNLO-terms. Table 1 shows that the residual renormalization scale dependence of the NNLO total cross-section is negligibly small for $\mu_{r}\in[m_{t}/2,2m_{t}]$, which is less than $0.1\%$ even when taking a quite large initial renormalization scale range $\mu_{r}\in[m_{t}/4,20m_{t}]$. The PMC scales are distinct at different orders, as in QED. Since the PMC scales are determined from perturbative input, any renormalization scale uncertainty of the pQCD series is transferred at finite order to the small uncertainty of the PMC scales.

If setting $\mu_{r}=m_{t}/2$ for conventional scale setting, the pQCD convergence is better than the cases of $\mu_{r}=m_{t}$ and $\mu_{r}=2m_{t}$, whose total cross-section is also close to the PMC prediction. Thus, for conventional scale setting, the best choice of an effective renormalization scale for top-quark pair production is $\mu_{r}\sim m_{t}/2$ other than the conventional suggested $m_{t}$ Wang:2014sua . The choice of $\mu_{r}\sim m_{t}/2$ is also suggested in Ref.Czakon:2016dgf by using the principle of fastest perturbative convergence.

Table 1 shows the PMC predictions for the top-quark pair total cross-section: $\sigma^{1.96\rm TeV}_{\rm Tevatron}=7.43^{+0.14}_{-0.13}~{}\rm pb$ at the Tevatron, $\sigma^{7\rm TeV}_{\rm LHC}=175.0^{+3.5}_{-3.5}~{}\rm pb$, $\sigma^{8\rm TeV}_{\rm LHC}=249.2^{+5.0}_{-4.9}~{}\rm pb$, and $\sigma^{13\rm TeV}_{\rm LHC}=807.8^{+16.0}_{-15.8}~{}\rm pb$ at the LHC for $\sqrt{S}=7$, 8 and 13 TeV, respectively. Those predictions agree with the Tevatron and LHC measurements Aaltonen:2013wca ; Chatrchyan:2013ual ; Aad:2012vip ; Chatrchyan:2013kff ; Aad:2015dya ; Chatrchyan:2012vs ; Aad:2012qf ; Chatrchyan:2016abc ; Aad:2014kva ; Khachatryan:2016mqs ; Khachatryan:2015fwh ; Khachatryan:2014loa ; Aad:2015pga ; Khachatryan:2016kzg ; Khachatryan:2015uqb ; Aaboud:2015AAAA ; Aaboud:2016pbd . A comparison of the PMC prediction for the top-quark pair production cross-section with the LHC measurements is shown in Fig.(1) for $\sqrt{S}=7$ TeV and $8$ TeV. As in Ref.Khachatryan:2016mqs , the theoretical error bands in Fig.(1) is estimated by using the CT14 error PDF sets Dulat:2015mca with the range of $\alpha_{s}(M_{Z})\in[0.117,0.119]$.

It is important to study the ratio of total cross sections $R^{8/7}=(\sigma^{8\rm TeV}_{\rm LHC})/(\sigma^{7\rm TeV}_{\rm LHC})$, since the experimental uncertainties, which are correlated between the two analyses (at $\sqrt{S}$ = 7 or 8 TeV) cancel out, leading to an improved precision in comparison to the individual measurements. The predicted cross-section ratio by the PMC is $R^{8/7}|_{\rm PMC}=1.42\pm 0.04$, which shows excellent agreement with the latest CMS measurement $R^{8/7}|_{\rm CMS}=1.43\pm 0.04\pm 0.07\pm 0.05$ Chatrchyan:2016abc .

As we have shown above, the PMC provides a comprehensive and self-consistent pQCD explanation for the top-quark pair production cross-section as well as the top-quark pair forward-backward asymmetry. The behavior of the top-quark pair production cross section allows a direct determination of the top-quark pole mass by comparing the pQCD prediction with the data.

Following the method of Ref.Aaboud:2016000 , we define a likelihood function

[TABLE]

Here $f_{\rm th}(\sigma|m_{t})$ is the normalized Gaussian distribution, which is defined as

[TABLE]

The top-quark pair production cross-section is a function of the top-quark pole mass $m_{t}$, decreases with increasing $m_{t}$. It can be parameterized as Beneke:2011mq

[TABLE]

where all masses are given in units of GeV. $\Delta\sigma_{\rm th}(m_{t})$ stands for the maximum error of the cross-section for a fixed $m_{t}$; it is estimated by using the CT14 error PDF sets Dulat:2015mca with range of $\alpha_{s}(M_{Z})\in[0.117,0.119]$. The determined coefficients $c_{0,1,2,3}$ are given in Table 2.

In order to determine the precise values for the coefficients $c_{0,1,2,3}$, we have used a wide range of the top-quark pole mass, i.e. $m_{t}\in[{\rm 160\;GeV},{\rm 190\;GeV}]$. We define $\sigma_{\rm th}(m_{t})$ as the cross-section at a fixed $m_{t}$, where all input parameters are taken at their central values, $[\sigma_{\rm th}(m_{t})+\Delta\sigma^{+}_{\rm th}(m_{t})]$ is the maximum cross-section within the allowable parameter range, and $[\sigma_{\rm th}(m_{t})-\Delta\sigma^{-}_{\rm th}(m_{t})]$ is the minimum value. Similarly, $f_{\rm exp}(\sigma|m_{t})$ is the normalized Gaussian distribution

[TABLE]

where $\sigma_{\rm exp}(m_{t})$ is the measured cross-section, and $\Delta\sigma_{\rm exp}(m_{t})$ is the uncertainty for $\sigma_{\rm exp}(m_{t})$.

We present the top-quark pair NNLO production cross-section (5) versus the top-quark pole mass at different hadron-hadron collision energies in Figs.(2, 3). The coefficients $c_{0,1,2,3}$ are determined by the PMC predictions. In these figures, the experimental measurements are presented for comparison, where the thinnest shaded bands are for the PMC predictions and the thickest shaded bands are for the combined experimental results respectively. The agreement of the PMC predictions with the measurements, as shown by Figs.(2, 3), makes it possible to achieve reliable predictions for top-quark pole mass. A precise range of values for the pole mass can thus be achieved in comparison with pQCD predictions based on conventional scale-setting. In the following, we will determine the top-quark pole mass such that the maximum value of the likelihood function (3) is achieved.

The D0 collaboration determined the top-quark pole mass by comparing the theoretical predictions based on conventional scale-setting with the measurements of the top-quark pair production cross-sections at the Tevatron Abazov:2009si ; Abazov:2011cq ; Abazov:2011mi ; Abazov:2016ekt . The results for various production channels are presented in Table 3. As a comparison, we present our predictions using PMC scale-setting in Table 3. For the calculation of the likelihood function (3), we have used the experimental measurements in these references as the input for $f_{\rm exp}(\sigma|m_{t})$.

Table 3 shows that the top-quark pole mass determined from the dilepton channel which are measured at the Tevatron Run I stage possesses the largest uncertainty Abazov:2009si . It will be improved by more precise data for the dilepton and the lepton + jets channels obtained at the Run II stage Abazov:2011cq ; Abazov:2011mi ; Abazov:2016ekt .

We present the likelihood function defined in Eq.(3) at the Tevatron in Fig.(4), where the measured combined inclusive top-quark pair cross-section of Ref.Abazov:2016ekt are adopted as the experimental input. By evaluating the likelihood function, we obtain $m_{t}=174.6^{+3.1}_{-3.2}$ GeV, where the central value is extracted from the maximum of the likelihood function, and the error ranges are obtained from the $68\%$ area around the maximum. As indicated by Figs.(2, 3), due to the elimination of renormalization scale uncertainty. The PMC predictions have less uncertainty compared to the predictions by using conventional scale-setting. Thus the uncertainty of the precision of top-quark pole mass is dominated by the experimental errors. For example, the PMC determination for the pole mass via the combined dilepton and the lepton + jets channels data is about $1.8\%$, which is almost the same as that of the recent determination by the D0 collaboration, $172.8^{+3.4}_{-3.2}$ GeV Abazov:2016ekt whose error is $\sim 1.9\%$.

The CMS and ATLAS collaborations have determined the top-quark pole mass by using measurements of top-quark pair production cross-sections at the LHC Chatrchyan:2012bra ; Chatrchyan:2013haa ; Aad:2014kva ; Khachatryan:2016mqs together with the theoretical predictions derived from conventional scale-setting; the results for various production channels are presented in Tables 4 and 5 for $\sqrt{S}=7$ and 8 TeV, respectively. As a comparison, we also present our predictions using PMC scale-setting in the two Tables. Similarly, for calculating the likelihood function (3), we use the experimental measurements in those references as the input for $f_{\rm exp}(\sigma|m_{t})$.

By using the measured cross section $\sigma_{\rm exp}(m_{t})$ together with its error $\Delta\sigma_{\rm exp}(m_{t})$ from the latest CMS measurement Khachatryan:2016mqs , we present the likelihood functions at the LHC in Fig.(5). Because the experimental uncertainty at the LHC is smaller than that of Tevatron, the determined top-quark pole mass by using the LHC data has better precision in comparison with the analysis using the Tevatron data. By evaluating the likelihood functions, we obtain $m_{t}=173.7\pm 1.5$ GeV for $\sqrt{S}=7$ TeV, and $m_{t}=174.2\pm 1.7$ GeV for $\sqrt{S}=8$ TeV. The precision of the top-quark pole masses determination is improved to be (${\pm 1.5}$) for $\sqrt{S}=7$ TeV and (${\pm 1.7}$) for $\sqrt{S}=8$ TeV.

By evaluating the likelihood function (3) using the corresponding measurements of the latest Tevatron and LHC collaborations, we obtain the following predictions for the top-quark pole mass,

[TABLE]

By using the relation between the pole mass and the $\overline{\rm MS}$ mass up to four-loop level Marquard:2015qpa ; Kataev:2015gvt , we can convert the top-quark pole mass to the $\overline{\rm MS}$ definition. For $\mu_{r}=m_{t}$, we obtain

[TABLE]

The weighted average of those predictions then leads to

[TABLE]

We summarize the top-quark pole masses determined at both the Tevatron and LHC in Fig.(6), where our PMC predictions and previous predictions from other collaborations Abazov:2009ae ; Abazov:2011pta ; Abazov:2016ekt ; Aaboud:2016000 ; Aad:2014kva ; Aad:2015waa ; Chatrchyan:2013haa ; Khachatryan:2016mqs are presented. For reference, the combination of Tevatron and LHC direct measurements of the top-quark mass is presented as a shaded band, giving $m_{t}=173.34\pm 0.76$ GeV ATLAS:2014wva . It shows that our new top-quark pole mass determined by PMC agree with the combination of Tevatron and LHC direct measurements.

III Summary

We have achieved precise predictions for the top-quark pair production cross-section with minimal dependence on the choice of the initial renormalization scale by using PMC. The resulting predictions are in agreement with measurements done by both the Tevatron and the LHC Collaborations. We have given a new determination of the top-quark pole mass by comparing the PMC prediction for the top-quark pair cross-sections with the latest measurements; a detailed comparison of previous determinations given in the literature has also been presented. Our new determination of the top-quark pole masses provide complementary information compared to direct measurements.

The determined top-quark pole masses are cross-checked by other determinations used different techniques. Typically, the mass $m_{t}=175.8^{+2.7}_{-2.4}$ GeV from an electroweak fits Baak:2012kk , the mass $m_{t}=173.5\pm 3.0\pm 0.9$ GeV reconstructed from lepton + $J/\psi$ from b-jet Khachatryan:2016pek , the mass $m_{t}=173.2\pm 1.6$ GeV from dilepton kinematic distributions2017mlepton and the best direct measurement results $m_{t}=172.99\pm 0.85$ GeV from ATLAS Aaboud:2016igd and $m_{t}=172.35\pm 0.51$ GeV from CMS Khachatryan:2015hba . The consistency of the pQCD predictions using the PMC with all of the collider measurements at different energies and different techniques provides an important verification of QCD.

The PMC provides a systematic, rigorous method for eliminating renormalization scheme-and-scale ambiguities at each order in perturbation theory. As we have shown in our previous papers, the PMC is applicable to a wide variety of perturbatively calculable processes. In each case, the ad hoc renormalization scale uncertainty conventionally assigned to the pQCD predictions can be eliminated. The residual renormalization scale dependence due to uncalculated high-order terms are usually small due to a more convergent pQCD series. The PMC, with its solid physical and rigorous theoretical background, thus will greatly improve the precision of tests of the Standard Model.

Acknowledgements: The authors would like to thank Hui-Lan Liu for helpful discussions. This work was supported in part by the Natural Science Foundation of China under Grant No.11547010, No.11625520, No.11705033 and No.11325525; by the Project of Guizhou Provincial Department of Science and Technology under Grant No.2016GZ42963 and the Key Project for Innovation Research Groups of Guizhou Provincial Department of Education under Grant No.KY[2016]028 and No.KY[2017]067; and by the Department of Energy Contract No.DE-AC02-76SF00515. SLAC-PUB-16934.