Momentum-space threshold resummation in $tW$ production at the LHC

TL;DR

This paper performs all-order soft-gluon resummation for $tW$ production at the LHC, significantly improving cross section predictions and aligning well with experimental measurements.

Contribution

It introduces a momentum-space threshold resummation method for $tW$ production, enhancing theoretical accuracy beyond NLO.

Findings

NLO+NNLL results increase cross sections by 12-17%.

Resummation improves agreement with LHC measurements.

Provides predictions for 14 TeV LHC.

Abstract

We calculate the soft-gluon corrections for $tW$ production to all orders. The soft limit is defined in the pair invariant mass or one particle inclusive kinematic schemes. We find that at NLO the contribution of the soft-gluon effect dominates in the total cross section or the differential distributions. After resumming the soft-gluon effect to all orders using the renormalization group equation, we find that the NLO+NNLL results increase the NLO cross sections by $12\%\sim 17\%$ depending on the scheme and the collider energy. Our results are in agreement with the measurements at the 8 and 13 TeV LHC. We also provide predictions for the total cross section at the 14 TeV LHC.

Figures

Tables

Table 1. Table 1: The fixed-order total cross section and the power corrections for t ​ W − 𝑡 superscript 𝑊 tW^{-} production with s = 8 𝑠 8 \sqrt{s}=8 TeV. The power corrections are defined as d ​ σ ​ ( NLO ) − d ​ σ ​ ( NLO ​ leading ) 𝑑 𝜎 NLO 𝑑 𝜎 NLO leading d\sigma({\rm NLO})-d\sigma({\rm NLO\ leading}) . The LO results are different in two schemes due to the different choice of the factorization scale. The scale uncertainties are shown.

[pb]	PIM			1PI
	b-veto	DS	DR	b-veto	DS	DR
LO	${6.96}_{-6\%}^{+5\%}$			${7.21}_{-4\%}^{+5\%}$
NLO $bg$	${12.0}_{+2\%}^{-6\%}$			${11.7}_{+4\%}^{-8\%}$
NLO leading	${11.2}_{+2\%}^{-5\%}$			${11.5}_{+0\%}^{-3\%}$
NLO	${9.31}_{-1\%}^{+0\%}$	${9.92}_{-2\%}^{+2\%}$	${10.0}_{-2\%}^{+2\%}$	${9.32}_{+0\%}^{-1\%}$	${10.0}_{-2\%}^{+1\%}$	${10.2}_{-2\%}^{+2\%}$
power corr.	$-{1.87}_{-0.25}^{+0.5}$	$-{1.26}_{-0.40}^{+0.75}$	$-{1.17}_{-0.43}^{+0.76}$	$-{2.14}_{-0.05}^{+0.30}$	$-{1.46}_{-0.21}^{+0.53}$	$-{1.29}_{-0.27}^{+0.63}$

Table 2. Table 2: Total cross sections for t ​ W − 𝑡 superscript 𝑊 tW^{-} production in PIM and 1PI schemes. The NLO cross sections are calculated using DS scheme.

[pb]	PIM		1PI
$\sqrt{s}$	8 TeV	13 TeV	8 TeV	13 TeV
LO	${7.0}_{-6\%}^{+5\%}$	${22.4}_{-2\%}^{+5\%}$	${7.2}_{-4\%}^{+5\%}$	${22.9}_{-1\%}^{+3\%}$
NLO	${9.92}_{-2\%}^{+2\%}$	${32.8}_{-1\%}^{+1\%}$	${10.0}_{-2\%}^{+1\%}$	${33.0}_{-1\%}^{+1\%}$
aNNLO	${11.6}_{-5\%}^{+4\%}$	${37.1}_{-5\%}^{+5\%}$	${11.2}_{-6\%}^{+6\%}$	${35.9}_{-6\%}^{+7\%}$
NLO+NNLL	${11.4}_{-7\%}^{+7\%}$	${36.7}_{-7\%}^{+7\%}$	${11.7}_{-17\%}^{+12\%}$	${37.3}_{-21\%}^{+16\%}$
aNNLO/NLO	1.16	1.13	1.12	1.09
(NLO+NNLL)/NLO	1.15	1.12	1.17	1.13

Equations

$$p(P_1)+p(P_2)\to t(p_3)+W^{-}(p_4)+X(P_X),$$

$$b(p_1)+g(p_2)\to t(p_3)+W^{-}(p_4)+X(p_X).$$

$$\hat{s}=(p_1+p_2{)}^2,{\hat{t}}_1=(p_1-p_3{)}^2-m_t^2,{\hat{u}}_1=(p_2-p_3{)}^2-m_t^2,$$

$${\hat{t}}_1^W=(p_2-p_4{)}^2-M_W^2,{\hat{u}}_1^W=(p_1-p_4{)}^2-M_W^2.$$

$$s=\hat{s}/x_1/x_2,t_1={\hat{t}}_1/x_1,u_1={\hat{u}}_1/x_2,t_1^W={\hat{t}}_1^W/x_2,u_1^W={\hat{u}}_1^W/x_1,$$

$$\frac{d^2{\sigma }^{\mathrm{PIM}}}{d{M}^{2}d\cos \theta }=$$

$$\displaystyle\times\frac{z^{-\eta}}{(1-z)^{1-2\eta}}\tilde{s}_{\rm PIM}\left(\ln\frac{M^{2}(1-z)^{2}}{z\mu_{s}^{2}}+\partial_{\eta},\mu_{s}\right)\frac{e^{-2\gamma_{\rm E}\eta}}{\Gamma(2\eta)}\Bigg{|}_{\eta=(C_{A}+C_{F})a_{\gamma^{\rm cusp}}(\mu_{s},\mu_{f})},$$

$$U_{\mathrm{PIM}}({\mu }_h,{\mu }_s,{\mu }_f)=$$

$$\displaystyle\hskip 68.28644pt+\,a_{\gamma^{h}}(\mu_{s},\mu_{h})+2\,a_{\gamma^{\phi_{q}}}(\mu_{s},\mu_{f})+2\,a_{\gamma^{\phi_{g}}}(\mu_{s},\mu_{f})\Big{]}.$$

$${\gamma }^h=2\left({\gamma }^Q+{\gamma }^q+{\gamma }^g\right)-\frac{C_A}{2}{\gamma }^{\mathrm{cusp}}\ln \frac{m_t^2(-\hat{s})}{{\hat{u}}_1^2}-\left(C_F-\frac{C_A}{2}\right){\gamma }^{\mathrm{cusp}}\ln \frac{m_t^2(-\hat{s})}{{\hat{t}}_1^2}.$$

$$\frac{d^2{\sigma }^{1PI}}{d{p}_T^{2}dy}=$$

$$\times H({\mu }_h)U_{1PI}({\mu }_h,{\mu }_s,{\mu }_f){\tilde{s}}_{1PI}({\partial }_{\eta },{\mu }_s)\frac{1}{s_4}{\left(\frac{s_4}{M_W{\mu }_s}\right)}^{2\eta }\frac{e^{-2{\gamma }_{\mathrm{E}}\eta }}{\Gamma (2\eta )},$$

$$U_{1PI}({\mu }_h,{\mu }_s,{\mu }_f)$$

$$\displaystyle+2\,a_{\gamma^{\phi_{q}}}(\mu_{s},\mu_{f})+2\,a_{\gamma^{\phi_{g}}}(\mu_{s},\mu_{f})+a_{\gamma^{\rm cusp}}(\mu_{s},\mu_{f})\Bigg{(}C_{A}\ln\frac{M_{W}^{2}\mu_{s}^{2}}{\left(\hat{t}_{1}^{\,W}\right)^{2}}+C_{F}\ln\frac{M_{W}^{2}\mu_{s}^{2}}{\left(\hat{u}_{1}^{W}\right)^{2}}\Bigg{)}\Bigg{]}.$$

$$\frac{d}{d\ln \mu }\left(f_{i/p}\otimes f_{j/p}\otimes H\otimes S_{\mathrm{PIM},1PI}\right)=0$$

$$d\sigma (\mathrm{aNNLO})=d\sigma (\mathrm{NNLO}\mathrm{leading})+d\sigma (\mathrm{NLO})-d\sigma (\mathrm{NLO}\mathrm{leading}),$$

$$d\sigma (\mathrm{NLO}+\mathrm{NNLL})=d\sigma (\mathrm{NNLL})+d\sigma (\mathrm{NLO})-d\sigma (\mathrm{NLO}\mathrm{leading})$$

$$\frac{d}{d\ln \mu }{f}_{i/p}(x,\mu )={\int }_y^1\frac{dz}{z}\left[\frac{2C_i{\gamma }^{\mathrm{cusp}}({\alpha }_s)}{(1-z{)}_{+}}+2{\gamma }^{{\varphi }_i}({\alpha }_s)\delta (1-z)\right]f_{i/p}(x/z,\mu ),$$

$$\frac{d}{d\ln \mu }{f}_{q/p}(x_1(s_4),\mu )$$

$$+\left[2C_F{\gamma }^{\mathrm{cusp}}({\alpha }_s)\ln \frac{s_4}{-{\hat{u}}_1^W}+2{\gamma }^{{\varphi }_q}({\alpha }_s)\right]f_{q/p}(x_1(s_4),\mu ),$$

$$\frac{d}{d\ln \mu }{f}_{g/p}(x_2(s_4),\mu )$$

$$+\left[2C_A{\gamma }^{\mathrm{cusp}}({\alpha }_s)\ln \frac{s_4}{-{\hat{t}}_1^W}+2{\gamma }^{{\varphi }_g}({\alpha }_s)\right]f_{g/p}(x_2(s_4),\mu ).$$

$$\displaystyle\tilde{s}_{\rm NLO}(L,\mu_{s})=\frac{\alpha_{s}}{4\pi}\bigg{[}-C_{A}I_{12}-(2\,C_{F}-C_{A})I_{13}-C_{A}I_{23}+C_{F}I_{33}\bigg{]},$$

$$\displaystyle I_{ij}(L)=-\frac{(4\pi\mu)^{2\epsilon}}{\pi^{2-\epsilon}}v_{i}\cdot v_{j}\int d^{d}k\frac{e^{-ik^{0}x_{0}}}{v_{i}\cdot kv_{j}\cdot k}(2\pi)\delta(k^{2})\theta(k^{0})\Bigg{|}_{L_{0}\to-L},$$

$$I_{12}$$

$$I_{13}$$

$$I_{33}$$

$$I_{23}$$

Full text

††preprint:

CERN-TH-2019-014

LA-UR-19-21475

TUM-HEP-1190/19

1 Introduction

The top quark, the most massive elementary particle discovered so far, is playing an important role in testing the standard model (SM) and searching for new physics beyond the SM. At a hadron collider, such as the large hadron collider (LHC), the top quarks can be produced in pairs via strong interaction or in association with a jet or a $W$ boson via weak interactions. The associated production with a $W$ boson offers a particular window to the weak interactions of the top quark and potentially can lead to a direct measurement of the CKM matrix element $V_{tb}$. Besides, $tW$ production is the second largest single top production channel and thus serves as an essential background in search for new physics. So far, the LHC has accumulated a large number of data, based on which the total and differential cross sections of this channel have been measured directly Aad:2012xca ; Aad:2015eto ; Aaboud:2016lpj ; Aaboud:2017qyi ; Chatrchyan:2012zca ; Chatrchyan:2014tua ; Sirunyan:2018lcp .

On the other hand, the precise theoretical predictions provide a valid framework in which important information can be extracted from the experimental data. Since the leading-order (LO) cross section of the $W$ boson associated production is proportional to the strong coupling $\alpha_{s}$ and $|V_{tb}|^{2}$, it is crucial to understand the value of $\alpha_{s}$ in order to extract $V_{tb}$. The inclusion of higher order QCD corrections can help to estimate the factorization and renormalization scale dependence of the cross sections. The next-to-leading order (NLO) QCD corrections have been calculated in refs. Giele:1995kr ; Zhu:2001hw ; Cao:2008af ; Kant:2014oha , and the results including decays of the top quark and the $W$ boson are also available Campbell:2005bb . The parton shower effects in this channel have been studied in refs. Frixione:2008yi ; Re:2010bp ; Jezo:2016ujg .

When considering the higher order corrections for $tW^{-}$ production, there is one subtlety to deal with. Due to the same $tW^{-}b$ final state in both the real correction to $tW^{-}$ and the $t{\bar{t}}$ production with the decay ${\bar{t}}\to W^{-}b$, one needs to find a way to differentiate the two processes or to define the $tW^{-}$ process properly beyond tree level. There are several proposals on the market Tait:1999cf ; Belyaev:2000me ; White:2009yt ; Demartin:2016axk . Here we point out that so far we discuss the $tW^{-}$ production in the five flavor scheme, i.e., the LO is $gb\to tW^{-}$. It is also possible to work in the four flavour scheme in which the LO is $gg\to tW^{-}\bar{b}$ Frederix:2013gra ; Cascioli:2013wga .

The next-to-next-to-leading order (NNLO) QCD corrections to $tW$ production are not available for now. Though the NNLO $N$-jettiness soft function, one of the ingredients for an NNLO calculation using $N$-jettiness subtraction method, has been computed in refs. Li:2016tvb ; Li:2018tsq , the two-loop virtual correction is still the bottleneck because of its dependence on multiple scales. The soft gluon corrections near the threshold have been calculated up to NNNLO based on next-to-next-to-leading logarithms (NNLL) resummation Kidonakis:2006bu ; Kidonakis:2010ux ; Kidonakis:2016sjf , which are considered as an approximation to the full higher order corrections. The three-loop soft anomalous dimension for $tW$ production was calculated by ref. Kidonakis:2019nqa which can be used to study the full NNNLO threshold effects.

In this paper, we will present the soft gluon corrections to $tW$ production using the soft-collinear effective theory Bauer:2000ew ; Bauer:2000yr ; Bauer:2001ct ; Bauer:2001yt ; Beneke:2002ph (see Becher:2014oda for a review) which separates the hard contributions with the large momentum transfer and the soft gluon corrections characterized as low energy contributions. Two different definitions of the soft limit are investigated. One is measured by the threshold variable $1-z=1-M_{tW}^{2}/\hat{s}\to 0$, while the other is given by $s_{4}=(p_{1}+p_{2}-p_{t})^{2}-M_{W}^{2}\to 0$. In principle, these two definitions encode the same soft gluon physics in the threshold limits and they only differ by power suppressed corrections. The threshold contributions up to NNLO are obtained and the resummation is achieved through solving the RG equations of the hard and soft functions. Our results could be taken as an important theoretical input in future experimental analyses.

The paper is organized as follows. In section 2 we show the basic information about the kinematics in this process and the factorization formula of the cross section in the soft limit. The numerical results and relevant discussions are then presented in section 3. We conclude in section 4. The evolution equation of parton distribution functions (PDFs) in the threshold limits and the analytic result of the soft function are given in the appendices.

2 Factorization and resummation formalism

We consider inclusive stable top quark and $W$ boson associated production at the LHC

[TABLE]

where $X$ denotes all the other possible extra radiations in the final states. In the threshold limit at the leading power we only need to consider the partonic channel

[TABLE]

The corresponding LO Feynman diagrams are shown in figure 1. The partonic kinematic variables are defined to be

[TABLE]

The corresponding variables at hadronic level are

[TABLE]

where $x_{1,2}$ are the Bjorken scaling variables.

In the soft limit or threshold limit, the real emissions are highly constrained, only soft gluons allowed in the final state $X$. This limit would be reached if the invariant mass of the final state $M=\sqrt{(p_{3}+p_{4})^{2}}$ approaches the initial partonic center-of-mass energy $\sqrt{\hat{s}}$. As a result, the variable $1-z\equiv 1-M^{2}/\hat{s}\to 0$ in the threshold limit111In the central-of-mass frame the energy of the soft radiation is $E_{g}\approx M(1-z)/2\sqrt{z}$. , and the perturbative expansion of cross section contains a series of large logarithms $\alpha_{s}^{n}\left[\ln^{2n-i}(1-z)/(1-z)\right]_{+}~{}~{}(i=1,2,...,2n)$, which might spoil the convergence of perturbative series. It is our purpose in this work to study the threshold behavior and resum such large logarithms to all orders. Since the soft limit is characterized by the final-state two particles’ invariant mass, it is called the pair invariant mass (PIM) scheme. Besides, there is another scheme, called one particle inclusive (1PI) scheme, in which the soft limit is defined by partonic level $s_{4}\to 0$ with $s_{4}\equiv\hat{s}+\hat{t}_{1}+\hat{u}_{1}+m_{t}^{2}-M_{W}^{2}=(p_{4}+p_{X})^{2}-M_{W}^{2}$. The two schemes measure the soft limit in different ways and the combination of the studies in two schemes provides more complete information on the structure.

In the rest part of this section we briefly show the factorization formula for $tW$ production, which can be derived in a similar way used in the other processes, such as the single top or top quark pair productions Ahrens:2010zv ; Ahrens:2011mw ; Zhu:2010mr ; Wang:2010ue ; Wang:2012dc . In the PIM scheme the cross section in the threshold limit can be factorized to a product of the hard and the soft function 222This factorization is carried out at the leading power of the threshold variable. At next-to-leading power, the threshold factorization becomes more complicated; see the recent paper Beneke:2018gvs . Idilbi:2005ky ; Becher:2007ty , which describes the physics at two different scales, i.e., the large hard scale and the small soft scale, respectively. They contain no large logarithms at their intrinsic scales ($\mu_{h}$ and $\mu_{s}$ respectively) as expected, since they depend only on a single scale there. The resummation of all the large logarithms caused by soft gluon effects is achieved by evolving the two functions from the intrinsic scales to a common factorization scale using their renormalization group (RG) equations. The RG improved differential cross section can be written as

[TABLE]

where polar angle of top quark $\theta$ is defined in the center-of-mass frame, $\lambda=(1-m_{t}^{2}/\hat{s}-M_{W}^{2}/\hat{s})^{2}-4m_{t}^{2}M_{W}^{2}/\hat{s}^{2}$, $\tau=M^{2}/s$ and the other kinematic variables have already been defined in eq. (2). For convenience we suppress the dependence of hard function $H$, soft function $\tilde{s}_{\rm PIM}$ and evolution factor $U_{\rm PIM}$ on the kinematic variables. $\tilde{s}_{\rm PIM}$ denotes the soft function defined in the Laplace space. The evolution factor embodies RG running from hard and soft scale to factorization scale, which is expressed as

[TABLE]

where the definitions of function $S$ and $a_{\gamma}$ and all relevant anomalous dimensions, e.g., $\gamma^{\phi_{q}}$, can be found in appendix A of ref. Ahrens:2010zv . The hard anomalous dimension specific to the $tW$ process is given by Ferroglia:2009ep ; Ferroglia:2009ii

[TABLE]

In the 1PI scheme, the soft radiations are characterized via the threshold variable $s_{4}\approx 2p_{4}\cdot k$ with $k$ the sum of all the momenta of the soft final state Laenen:1998qw ; Becher:2009th ; Ahrens:2011mw ; Wang:2012dc . The analysis is simplified by going to the rest frame of the inclusive final state $W+X$ where $|\vec{p}_{W}|=\mathcal{O}(s_{4}/M_{W})$. In this frame the energy of the soft radiation is $E_{g}\approx s_{4}/2M_{W}$. In the 1PI scheme, the RG improved cross section is

[TABLE]

where the integration range is defined as $x_{1}^{\min}=(M_{W}^{2}-m_{t}^{2}-u_{1})/(s+t_{1})$, $s_{4}^{\rm max}=m_{t}^{2}-M_{W}^{2}+u_{1}+x_{1}(s+t_{1})$, and momentum fraction $x_{2}$ is defined via $x_{1}$ and $s_{4}$ as $x_{2}=(s_{4}+M_{W}^{2}-m_{t}^{2}-x_{1}t_{1})/(u_{1}+x_{1}s)$. The Mandelstam variables are related to the top quark’s transverse momentum $p_{T}$ and rapidity $y$ via $t_{1}=-\sqrt{s}\,m_{\bot}e^{-y}$, and $u_{1}=-\sqrt{s}\,m_{\bot}e^{y}$, with $m_{\bot}=\sqrt{\smash[b]{p}_{T}^{2}+m_{t}^{2}}$. Here the 1PI evolution factor has the form as

[TABLE]

The difference between the RG factors in eq. (2) and eq. (2) arises from the RG equation of the PDFs and soft function. We present the RG evolution of the PDFs in the PIM and 1PI schemes in appendix A. The two-loop anomalous dimensions that govern the evolution of the hard and soft functions and thus determine the scale dependent part are derived from the general structure of the anomalous dimension Ferroglia:2009ep ; Ferroglia:2009ii . The scale independent part of the hard function has been obtained at NLO using modified MadLoop Hirschi:2011pa which makes uses of Ninja Peraro:2014cba , CutTools Ossola:2007ax and OneLOop vanHameren:2009dr packages. We have computed one-loop soft function analytically, which is shown in appendix B. Combining all the ingredients together, we have checked the RG invariance

[TABLE]

in both of the kinematic schemes.

The NLO and NNLO leading power contributions are obtained by setting the scales $\mu_{h},\mu_{s},\mu_{f}$ in Eqs. (2) and (2) equal. In this way, for PIM scheme we capture all the threshold logarithms $\alpha_{s}\left[\ln^{n}(1-z)/(1-z)\right]_{+}$ with $n=1,0$ at NLO and $\alpha_{s}^{2}\left[\ln^{n}(1-z)/(1-z)\right]_{+}$ with $n=3,2,1,0$ at NNLO, as well as the scale dependent logarithms predicted by eq. (2). The similar procedure can be applied to obtain the threshold enhanced logarithms for 1PI scheme. In the following calculations the approximate NNLO (aNNLO) cross section is defined as

[TABLE]

where the NLO power suppressed terms in $1-z$ or $s_{4}$ have been included to give more precise results. We can also match the resummed prediction to the fixed order result by

[TABLE]

with the NNLL result given by eq. (2) and eq. (2).

3 Numerical results

To perform the numerical calculation, the input parameters are set as $m_{t}=173.3$ GeV, $\Gamma_{t}=1.5$ GeV, $M_{W}=80.419$ GeV, $\alpha=1/132.5$ and the Fermi-constant $G_{F}=1.166390\times 10^{-5}$ GeV*-2*. For the LO and NLO calculations we use the CT14 LO and NLO PDF sets Dulat:2015mca as provided by the LHAPDF library Buckley:2014ana , respectively. The aNNLO and resummed predictions are obtained using CT14 NNLO PDF sets. For fixed-order calculations the renormalization scale is set to be the same as the factorization scale. It is natural to set the default hard scale to be the invariant mass of the top quark and $W$ boson, i.e., $\mu_{h}=M$, where the hard function contains no large logarithms. The soft scale is chosen numerically according to the criterion that the perturbative series of the soft function are well behaved Becher:2007ty . Explicitly, we find that the ratio of the soft and hard scale is $0.3\sim 0.4$ in the PIM scheme and $0.3\sim 0.5$ in the 1PI scheme at the 8 TeV LHC. The default factorization scale has been chosen to be $\mu_{f}=M$ and $\mu_{f}=m_{t}+M_{W}$ in the PIM and the 1PI schemes, respectively. The final scale uncertainties are evaluated by varying these scales by a factor of two independently.

As discussed in the introduction there are several methods to deal with the problem of the interference between the real corrections to $tW$ production and $t\bar{t}$ calculation  Tait:1999cf ; Belyaev:2000me ; Campbell:2005bb ; Frixione:2008yi ; White:2009yt ; Demartin:2016axk at NLO, such as Diagram Removal (DR), Diagram Subtraction (DS) and b-jet transverse momentum veto. The differences between these schemes have been discussed a lot before; see e.g. refs. Frixione:2008yi ; Re:2010bp . Since they are only relevant in the power suppressed channels333The problem of the interference exists only in the $gg\to tW^{-}\bar{b}$ or $q\bar{q}\to tW^{-}\bar{b}$ channel which is at subleading power near the threshold., we will not repeat the discussion about their difference in this paper. We notice that in refs. Kidonakis:2010ux ; Kidonakis:2016sjf where the higher order threshold corrections were studied, the power corrections as well as certain leading power logarithmic independent terms are not taken into account. In the rest of this section we only show the predictions of the process $pp\to tW^{-}$. The total cross section for $tW^{-}$ and $\bar{t}W^{+}$ can be obtained by doubling the results, as demonstrated in ref. Frixione:2008yi . The NLO cross sections in the b-jet veto scheme are evaluated using MCFM Campbell:2005bb with $p_{T}^{b-{\rm jet}}<50$ GeV. The cross sections in the DS and DR schemes are calculated by POWHEG-BOX Alioli:2010xd ; Re:2010bp .

Before presenting the resummed result, we firstly investigate the contribution of the leading power terms. From table 1 we can see that the NLO corrections are sizeable, enhancing the LO result by $29\%\sim 44\%$ depending on the different methods to isolate the $tW$ process. These NLO corrections get contributions from all the $bg$, $gg$ and $qq^{\prime}$ channels, though at LO only $bg$ channel exists. Among them, the $bg$ channel dominates or even surpasses the NLO corrections, as indicated in table 1 too. Moreover, the leading power terms of the $bg$ channel can approximate the total result of the $bg$ channel very well, the difference being only $7\%$ and $2\%$ in the PIM and 1PI schemes, respectively. Since the leading power terms can be obtained from the resummed results as discussed above in eq.(11), they can be calculated up to higher orders in $\alpha_{s}$, namely beyond NLO. These make up a major part of the full NNLO corrections and can be taken as an approximation of the latter. The quality of the approximation could be estimated by looking at the power corrections. The NNLO results are still unavailable, so we study the NLO ones which are shown in table 1 as well. It is ready to see that they are negative and around $-20\%\sim-12\%$ in PIM and $-23\%\sim-13\%$ in 1PI kinematic scheme depending on the methods to deal with the interference problem. The contributions of the higher order (in $\alpha_{s}$) power corrections can be obtained by calculating the full NNLO QCD corrections or by making use of the next-to-leading power factorization and resummation, both of which are difficult at the moment and beyond the scope of this paper.

Although the usual way to evaluate the scale uncertainty is to vary the scales by a factor of two, it is also interesting to investigate the factorization scale dependence in a larger region. From the figure 2, we can see that the ratio of the NLO over LO result is insensitive to the factorization scale, always in the region $(1.37,1.53)$, when it is varied from $M/8$ to $4M$. This means that there is no clear choice of the factorization scale to ensure fastest convergence. Moreover, we find that the $bg$ channel is very sensitive to the factorization scale when it is smaller than $M/2$. In order to avoid such a dependence, we have chosen the default factorization scale at $M$. It can also be seen that the NLO leading power terms dominate the $bg$ channel over a large region.

Then we turn to the differential cross sections. We show the $tW$ invariant mass distributions in the PIM scheme in figure 3 and the top quark $p_{T}$ distributions in the 1PI scheme in figure 4. We show results at both the 8 TeV and 13 TeV LHC. It can be seen that the leading power terms are dominant in all the invariant mass or the top $p_{T}$ regions, as in the case of total cross sections. The NNLO leading terms increase the NLO leading cross section by about $10\%$ in most of the region.

Now we present in table 2 the aNNLO and NLO+NNLL result defined in eq. (11) and eq. (12), respectively. The NLO+NNLL (aNNLO) predictions increase the NLO total cross section by $12\%\sim 17\%$ ($9\%\sim 16\%$) depending on the collider energy and the threshold variable, but with larger scale uncertainties. These large uncertainties are mainly from the variation of the factorization scale $\mu_{f}$. At first sight, this is unexpected since we have checked the scale independence near the threshold analytically in eq. (10). However, this is based on the assumption $x_{1,2}\to 1$ as discussed in appendix A. When the kinematics is far away from the threshold limit, this assumption is not valid. The very small scale uncertainties of the NLO results seem like a coincidence because the NLO contributions from $gg$ and $qq^{\prime}$ channels are negative while the contributions from $bg$ channel are positive. Meanwhile they display an opposite behavior under the scale variation; see table 1. Our resummed result or its expansion in $\alpha_{s}$ improves only the result in $bg$ channel. It would be interesting to investigate whether the scale cancellation among different channels happens at higher orders. From table 2 we also find that the total cross sections in the PIM and 1PI scheme are compatible. And the resummed cross sections in PIM kinematics have smaller scale uncertainties.

Lastly, we compare the theoretical results with the measurements of the total cross section for $tW^{-}$ and $\bar{t}W^{+}$ production at the LHC in figure 5. After considering the large experimental uncertainties, the NLO+NNLL predictions are in good agreement with the data at the 8 TeV and 13 TeV LHC. We also give the predictions at the 14 TeV LHC.

4 Conclusion

We have investigated the soft-gluon resummation for $tW$ production in the framework of soft-collinear effective theory. We considered the two different definitions of the threshold limit, $1-M^{2}/\hat{s}\to 0$ and $s_{4}\to 0$, corresponding to the PIM and 1PI kinematic schemes, respectively. We briefly discussed the factorization and resummation formalism in both kinematic schemes. In addition, we have calculated the hard function and soft function at NLO. Expanding the resummed formula in $\alpha_{s}$ gives the leading power terms of the fixed-order results. We found that the NLO leading power contribution is a good approximation to the $bg$ channel at NLO not only for the total cross sections but also for the differential distributions. After resumming the soft gluon effects to all orders using renormalization group equation, we find that the NLO+NNLL results increase the NLO cross sections by about 15(12)% in PIM and 17(13)% in 1PI scheme at the 8(13) TeV LHC, but with large uncertainties which is mostly generated by varying the factorization scale. We compared with the data at the 8 and 13 TeV LHC and found good agreement within uncertainties. We provide the prediction for the 14 TeV LHC.

In future, we can obtain more precise predictions for the $tW$ process by including higher order hard and soft functions in the resummation formalism or by calculating the full NNLO corrections. The latter may be achieved making use of the $N$-jettiness subtraction method Gao:2012ja ; Gaunt:2015pea ; Boughezal:2015dva . The NNLO beam function Gaunt:2014xga ; Gaunt:2014cfa and $N$-jettiness soft function Li:2018tsq for this process have been computed. The only missing part is the two-loop hard function, which requires a huge amount of work. We defer this study to future work.

Acknowledgements

We would like to thank Shi Ang Li for the contribution in the early stage of this work. C.S.L. was supported by the National Nature Science foundation of China, under Grants No. 11875072. H.T.L. was supported by the Los Alamos National Laboratory LDRD program. J.W. was supported by the BMBF project No. 05H15WOCAA and 05H18WOCA1.

Appendix A RG equation of the PDFs near the threshold

In the threshold limit $x_{1,2}\to 1$, the DGLAP evolution for the PDFs can be written as Korchemsky:1992xv ; Moch:2004pa

[TABLE]

where the quadratic Casimir operator $C_{i}$ for the quark is $C_{q}=C_{F}$, and for the gluon is $C_{g}=C_{A}$. In the threshold limit $s_{4}\to 0$ the evolution equations for PDFs are

[TABLE]

A similar derivation for $t\bar{t}$ and single top production can be found in refs. Ahrens:2011mw ; Wang:2012dc .

Appendix B Soft function

The NLO soft function can be written as

[TABLE]

For convenience we evaluate soft integral $I_{ij}$ in the position space, and then transform them into Laplace space as

[TABLE]

with $L_{0}=\ln\left(-\mu^{2}x_{0}^{2}e^{2\gamma_{\rm E}}/4\right)$, and $v_{i}$ are normalized momenta fulfilling on-shell conditions as $v_{1}^{2}=v_{2}^{2}=0$ and $v_{3}^{2}=1$.

In the PIM kinematics the full set of integrals can be found in section III of ref. Li:2014ula . In the 1PI kinematics the integral $I_{12}$ can be obtained from eq. (21) of ref. Ahrens:2011mw by replacing all the kinematics variables to the ones related to $W$. The integrals $I_{13}$, $I_{23}$ and $I_{33}$ are more complicated, and have been first calculated in this paper. The non-vanishing integrals are collected below

[TABLE]

with $\beta_{t}=\sqrt{1-4m_{t}^{2}\hat{s}/(m_{t}^{2}-M_{W}^{2}+\hat{s})^{2}}$, $\beta_{W}=\sqrt{1-4M_{W}^{2}\hat{s}/(m_{t}^{2}-M_{W}^{2}-\hat{s})^{2}}$, $x_{tW}=\sqrt{x_{t}\ x_{W}}$ and $x_{i}=(1-\beta_{i})/(1+\beta_{i})$. In the limit of $M_{W}\to m_{t}$, the integrals reproduce those for $t\bar{t}$ production in ref. Ahrens:2011mw .