Singletop production and rare top interactions
Abstract
The study of the top quark’s properties is an important part of the LHC programme. In earlier work, we have studied the rare decay , using effective operators to capture the effects of physics beyond the Standard Model. However top decay is primarily sensitive to new physics in the subTeV energy regime. If this new physics resides at a higher energy scale, then one needs to turn to singletop production. In this paper, we use the channel and channel singletop production cross sections to constrain the new physics parameter space associated with such contact interactions. We also study the net top polarization as a means to distinguish between contributions from operators involving different fermion chiralities and Lorentz structures.
pacs:
14.65.HaI Introduction
The top quark has long been believed to be colluding with new physics (NP). However, intense scrutiny of the top quark’s properties at the Tevatron and at the LHC has so far not revealed any conclusive departures from the Standard Model (SM). For a few years, the Tevatron experiments reported a large forwardbackward asymmetry in top pairproduction. With the accumulation of more statistics and improved calculation of the SM predictions, however, this anomaly disappeared Aaltonen:2017efp (). Nevertheless, the top quark remains a likely suspect – its mass differs from those of other SM fermions by orders of magnitude, so much so that it threatens to push the electroweak vacuum beyond the edge of stability. Moreover, longstanding and freshly surfaced anomalies in the sector continue to fuel speculations that third generation quarks may be the muchsought window to physics beyond the SM.
In earlier work tbbc_sofar (), we proposed a study of rare decay modes of the top quark. Since all the top quark measurements to date have been made in channels involving the dominant decay modes of the top, they would, naturally, have missed signs of new physics that only manifests in the rare decay modes. We examined the sensitivity of the LHC to the rare decay and found that with 3000 fb of data, the LHC would be able to set statistically significant limits on such decays. However, it is evident that top decay would be most sensitive to new physics effects arising at the energy scale of a few hundred GeV at most. If the new physics contributions originate at higher energy scales, the impact on top decay would be too small to be discernible. In order to probe such interactions further, one must turn to singletop production.
In this paper we examine effects that would arise in singletop production due to anomalous couplings between the top, bottom and charm quarks. Assuming such effects to arise at relatively high energy scales, we parametrize them in terms of various fourFermi operators the same operators that we considered in Ref. tbbc_sofar (). A priori, it may seem that the contribution of such operators to singletop production would be diminished by the parton densities of the heavy quarks in the initial state. While this is true in general, the situation is salvaged somewhat by the fact that the competing SM mode is driven by electroweak interactions and not by strong interactions. A detailed numerical study shows that it is possible to set meaningful limits on the parameters of the interactions using existing LHC data. We also present a futuristic scenario in which very stringent limits may be obtained. This possibility, however, is contingent upon the development of reliable techniques to determine the charge of an outgoing quark on an eventbyevent basis. We further examine the possibility of distinguishing between the contributions of new physics operators with different Lorentz structures and fermion chiralities using the polarization of the top quark.
Single top production has received significant attention as a direct probe of physics beyond the SM stp_old (); stp_new (). Since the present analysis focusses on the effects of fourquark operators, it is worth mentioning that such operators have been studied quite extensively, particularly in the context of flavorchanging neutral currents involving the top quark fourfermion ().
The remainder of this paper is organised as follows. In Section II, we discuss the theoretical generalities related to single top production at the LHC and introduce the effective operators that we use to parametrize new physics contributions, though the detailed analytic expressions are consigned to the Appendix. In Section III, we discuss our numerical analyses and results for channel as well as channel singletop production at the LHC at centerofmass energies of 8 TeV and 13 TeV. We conclude in Section IV.
Ii Single top production
At a hadron collider, the dominant production mode for top quarks is . Single top production is subdominant. Nonetheless, it is important as it provides a cleaner way of measuring the electroweak couplings of the top quark. Within the framework of the SM, singletop production at hadron colliders is classified into 3 production channels as shown in Fig. 1, namely, channel, channel and channel. At the LHC, channel production is the dominant mode, followed by associated production. Cross sections for all three channels have been measured at the LHC during the 7TeV and 8TeV runs ATLAS_t7 (); ATLAS_t8 (); ATLAS_s7 (); ATLAS_s8 (); CMS_t7 (); CMS_t8 (); CMS_s7n8 (). The 13TeV run is ongoing and results are already available in some channels ATLAS_t13 (); CMS_t13 (). These are summarised in Fig. 2. The measured cross sections are byandlarge in agreement with the SM predictions, even though certain channels are plagued by large experimental uncertainties.
ii.1 Effective Lagrangian
In Ref. tbbc_sofar (), the contributions from physics beyond the SM to the rare decay were parameterized in terms of the 6dimensional operators given by ,
where
(3)  
Clearly, these operators can also contribute to singletop production. The channel would remain unaffected by these NP contributions. However, final states identical to those produced in channel and channel processes can occur via the contact interactions listed above.
ii.2 Contribution to singletop production from
The operators listed in Eqs. (II.1), (II.1) and (3) can give rise to three possible amplitudes for single top production : , , . In the SM, the first one is an channel process, whereas the second and third are channel processes. In addition, the three final states get contributions from lightquark initial states in the SM. Some of the key features are as follows:

For singletop production due to such operators, the initial states would necessarily consist of bottom and charm quarks. The low densities of these inside the proton tend to suppress the cross section as compared to the SM production rates. This effect is more pronounced in the channel than in the channel as the SM rate for the latter is CKM supressed unless there is a quark in the initial state (see Fig. 1).

The suppression caused by the low initialstate parton densities is compensated to some extent by the growth in the cross section with increase in the parton c.m. energy (). This, of course, is a generic feature of contact interactions.

When new physics is parametrized in terms of effective operators, the implicit assumption is that the operators arise due to physics at a very high energy scale () that is beyond the direct reach of current experiments. When the operators listed in Eqs. (II.1)(3) are applied to top decay, (TeV) is permissible, given that the energy scale of the interaction is 173 GeV. However, in the context of singletop production in collisions at LHC energies, one must necessarily consider significantly larger values of (10 TeV) or higher. In the parameterization used in this work and in our previous work on this topic tbbc_sofar (), . Consequently, correspond to in the subTeV regime. In principle, it should be possible to rule out such new physics scenarios as they could lead to resonances in the channel production mode, causing a spike in the cross section. However, owing to the relatively large uncertainty in the measured channel cross section, no robust conclusions can be drawn at this stage. In the channel, such contributions would manifest, for example, in the form of a harder transverse momentum distribution. However, no such deviations have been found up to 300 GeV ATLAS_t8 ().
In our previous work related to , the 10 couplings from the NP effective Lagrangian were found to appear together in six characteristic combinations. We denoted these six combinations by , defining
(4) 
In addition to depending on the above six quantities, various observables were also found to depend on the real and imaginary parts of and on the combination .^{1}^{1}1The dependence on comes from SMNP interference terms. Terms proportional to are associated with tripleproduct correlations tbbc_sofar ().
In the case of singletop production, once again, NP contributions to the various cross sections of interest can be expressed in terms of these same combinations of NP parameters. Explicit expressions for the matrix elements squared in the different cases may be found in the Appendix. If we restrict our attention to the case in which we sum over the top quark’s spin, we find that the 20dimensional parameter space spanned by the 10 complexvalued parameters is reduced to a 5dimensional parameter space composed of , , , Re() and Im().^{2}^{2}2As we have noted in past work, the situation is somewhat complicated by the fact that the real and imaginary parts of also appear in the parameter . If one wishes to consider top quark polarization effects, there is an additional dependence on , , , and (please refer to the Appendix for details). In the present work we shall assume that and are both real, thereby reducing the size of the parameter space somewhat.^{3}^{3}3Note that plays an important role in partial rate asymmetries tbbc_sofar (); we neglect such effects here. In our numerical work below, we will mostly consider the case in which the polarization of the top quark is ignored (Secs. III.1, III.2 and III.4). A polarizationdependent asymmetry is considered in Sec. III.3.
From the experimental side, we currently have only one input that we can use to constrain this multidimensional parameter space, namely, the production cross section. Therefore, in order to obtain meaningful interpretations of these constraints, we will impose some further restrictions. Throughout the remainder of this paper, we will consider operators with a single Lorentz structure at a time. We will also make certain assumptions about the relative magnitudes of the couplings associated with operators having the same Lorentz structure but different chirality structure.
Iii Numerical Study
In order to extract limits on the ’s (and consequently on ), we start by implementing alongside the Standard Model in MadGraph5 MG5 () using FeynRules FR (). This puts us in a position to calculate the treelevel cross section for singletop production in collisions, for which we use CTEQ6L parton distributions functions (PDFs) CTEQ (), setting both the renormalization and factorization scales to be = 173 GeV. In order to approximate higherorder QCD corrections, we estimate Kfactors as (approx.)/, where (approx.) is obtained from the references listed in Table 1 and is calculated using MadGraph5 in conjuction with MSTW2008LO PDFs MSTW ().^{4}^{4}4MSTW2008LO PDFs are used only for the purpose of determining the Kfactor, to be consistent with the calculations for (approx.). For all subsequent calculations we use CTEQ6L PDFs. We then compute the treelevel cross sections, including both SM and NP effects, and multiply the results by the corresponding Kfactors to obtain estimates of the QCDcorrected values. In the following, “” denotes the SM cross sections obtained in this manner.
channel; 8 TeV  From Ref. NLO_sch_813 ()  From MadGraph5 using MSTW2008LO pdfs  1.74 

channel; 13 TeV  From Ref. NLO_sch_813 ()  From MadGraph5 using MSTW2008LO pdfs  1.73 
channel; 8 TeV  From Ref. NLO_tch_8 ()  From MadGraph5 using MSTW2008LO pdfs  1.06 
channel; 13 TeV  From Ref. CMS_t13 ()  From MadGraph5 using MSTW2008LO pdfs  1.02 
iii.1 channel singletop production
At the LHC, channel processes yield the dominant contribution to singletop production. These processes consist of and where and . Within the SM, the relative magnitudes of the contributions from the different initial states are governed by the densities of the respective partons inside the proton and the CKM factors appearing in the amplitude. Once is introduced, there are additional contributions to and . Figure 3 shows the cross section. As noted above, the operators in , and are considered separately in our analysis. It can be seen that the tensor operators are constrained most tightly, followed by vector and scalar operators. This is expected, given the structure of the s and the fact that large numerical factors accompany wherever they appear. It is intriguing to note that, except in the case of , couplings of (1) are not excluded by the experimental data at 8 TeV. How do we reconcile this with a) our earlier statement that (1) correspond to in the subTeV range, and b) the fact that no exotic physics has been discovered at the LHC so far? We return to this question below, in Sec.III.4.
Figure 4 shows the analogous analysis for the 13TeV data. As compared to the 8TeV data, the 13TeV data appear to be less constraining. At first glance, one is tempted to attribute this to low statistics given that the 8TeV measurement is based on 19.7 fb of data while the 13TeV result is based on 3.2 fb. However, despite the relatively low statistics at 13TeV, it turns out that the largest component of the uncertainty is due to systematics. If future analyses can reduce the systematic uncertainty, then tighter constraints can be expected. Presently, for a more effective comparison between the sensitivities to NP couplings at 8 TeV and 13 TeV, we construct a 10% band around the SM prediction (see Figs. 3 and 4). This gives us an estimate of the improvement in the limits under the assumption that, at both 8 TeV and 13 TeV, the central value of the measurement coincides with the SM prediction and has a 10% uncertainty.
In Figs. 3 and 4, we have allowed all chiral structures associated with a given Lorentz structure to have the same weight. That is, when vector operators are considered (), we have set , and similarly for scalar and tensor operators. In the following, we relax this condition and consider scenarios where , and are smaller than . As expected, this relaxes the constraint on (see Fig. 5). In a UVcomplete scenario, these operators may not all occur simultaneously and there would exist several possibilities for the relative sizes of the corresponding couplings. We illustrate the effect using one such possibility. The same exercise can be carried out for scalar and tensor operators; the results are depicted in Figs. 6 and 7.
iii.2 channel singletop production
Singletop production in the channel is due to the processes , with and . Of these, the largest contribution comes from , since the rest of the processes are CKM suppressed. However, the introduction of enhances the contribution due to . This is shown in Fig. 8.
Since the cross section is small ( 3.7 pb in the Standard Model), CMS CMS_s7n8 () and ATLAS ATLAS_s8 () report the combined cross section due to and . However, in our calculations, we consider only .^{5}^{5}5The NP operators considered in this work do, in fact, contribute to ; we are, however, restricting our attention to cases in which a single top quark is produced. Hence, in this section, the experimental measurements are not indicated on the plots. Instead, the grey band denotes a band of uncertainty around the central value of the corresponding SM cross section. The size of the band is chosen so as to be commensurate with the total uncertainty (statistical + systematic) reported in the actual measurement CMS_s7n8 (); ATLAS_s8 (). In Fig. 8, for example, the size of the band is 40%. While measurements in the channel are yet to be reported for 13 TeV, this data set is expected to be several times larger than the 8TeV data set. Hence we use two bands, of sizes 20% and 10%, to project the limits on in this case (see Fig. 9).
Comparing Fig. 8 with Fig. 3, it is immediately clear that the channel process, even when measured with a lower accuracy, yields more stringent bounds than the channel process. This is easy to understand. is much smaller for the channel process than for the channel process, so even if the relative uncertainty is larger, the absolute deviation allowed is smaller, which leads to tighter constraints on the couplings. A further improvement can be expected to emerge from the channel measurements at 13 TeV (see Fig. 9). Here again we consider scenarios where one chiral structure is dominant. The impacts on the corresponding limits are shown in Figs. 10, 11 and 12.
iii.3 Top polarization
The polarization of the top quark is an important observable at the LHC and a few measurements of it have already been made toppol (); toppol_single (). The net polarization is usually defined as
(5) 
where and , respectively, denote the number of top quarks with spin aligned along or opposite to a chosen direction. The value of depends on the choice of the reference direction. The usual choices for this reference direction are the axis or the momentum of the top itself. In the latter case, and denote the number of top quarks of different helicities. If the top quark is produced in association with another particle (as in the case of singletop production), the momentum of that particle may also be chosen as the reference direction. The utility of lies in the fact that it is sensitive to the chiral structure of the coupling at the production vertex. Fortuitously, due to its large mass, the top quark decays before hadronizing. This allows the polarization information (which would otherwise be lost during hadronization) to be gleaned from the angular distribution of the top quark’s decay products. Very often, the or coming from the top decay is used for this purpose.
In the preceding sections, we have identified regions of parameter space that are compatible with various singletop production cross section measurements. We now examine the reach of as a means to distinguish between the different types of couplings. In particular, we expect to deviate from its Standard Model value when there is an increase in the fraction of in the ensemble, that is, when there are contributions from or . As can be inferred from the expressions in the Appendix, cross section measurements are sensitive to various combinations of sums of ’s. By way of contrast, the spindependent contributions to the amplitude squared are dependent on differences of ’s, such as , for example. For this reason, polarization measurements can yield additional information regarding the types of NP interactions that contribute to the various processes under consideration in this work. In the following analysis, we choose certain values of that yield a cross section within a certain “allowed” range and then compute according to Eq. (5) above, with and denoting the number of top quarks of either helicity.
Table 2 lists channel singletop polarization asymmetries for various combinations of NP contributions at = 8 TeV. In choosing parameters, we have allowed for an enhancement of 20% in the cross section over the SM prediction, noting that the experimental uncertainty is close to 40%. We consider scenarios for which the dominant contribution comes from one operator,^{6}^{6}6This is also more likely from the point of view of a UVcomplete model. keeping the subdominant couplings at approximately onetenth of the dominant one. We find that the deviation in can be considerable. Note that the values of listed in Table 2 are based on the calculation of treelevel cross sections, as described at the beginning of Sec. III. factors cancel out in the calculation of . Estimation of based on higherorder calculations is beyond the scope of this work. However, the dominant higherorder corrections would arise from QCD effects and, as such, they are not expected to alter significantly.
Dominant Contribution  No Contribution  (pb)  

SM  SM  3.7  0.68  
VP1  ,  4.4  0.70  
VP2  ,  4.5  0.70  
VP3  ,  4.5  0.43  
VP4  ,  4.5  0.43  
SP1  ,  4.5  0.73  
SP2  ,  4.5  0.73  
SP3  ,  4.5  0.40  
SP4  ,  4.5  0.41  
TP1  ,  4.5  0.65  
TP4  ,  4.5  0.48 
We have not attempted to estimate the accuracy with which the can be measured at the LHC.^{7}^{7}7We refrain from making such an attempt, since background rejection depends heavily on the specific algorithm used, which, in turn, often involves boosted decision trees toppol_single () and other sophisticated analysis tools developed and trained by experimentalists to get the best from their respective detectors. We note, however, that a recent analysis from CMS toppol_single () was able to make a measurement with approximately 38% uncertainty. We also point out that the results from Ref. toppol_single () cannot be used to compare directly with our results because their choice of reference axis is different from ours. Nonetheless, it is clear that accuracies 10% or better would be needed in order for to be a useful discriminator.
Finally, in Table 3, we list the polarization asymmetry for channel singletop production at = 13 TeV. This time we allow for 10% enhancement in the cross section over the SM prediction. Once again the largest deviations correspond to contributions from and , although the net size of the deviation is smaller. This is largely an artefact of the restriction placed on the overall increase in cross section; allowing for a larger variation in the cross section would, in general, allow for greater deviation in . It would not, however, guarantee a greater deviation in , since the size of depends on the chiral structure of the dominant operator.
Dominant Contribution  No Contribution  (pb)  

SM  SM  7.1  0.67  
VP1  ,  7.8  0.69  
VP2  ,  7.8  0.69  
VP3  ,  7.8  0.55  
VP4  ,  7.8  0.53  
SP1  ,  7.8  0.71  
SP2  ,  7.8  0.71  
SP3  ,  7.8  0.54  
SP4  ,  7.8  0.52  
TP1  ,  7.9  0.65  
TP4  ,  7.9  0.56 
One can, similarly, obtain the polarization asymmetry for channel singletop production. If we follow our earlier strategy of choosing benchmark points from well within the allowed band, we would be considering scenarios where . In that case, the deviation in is 10% and may not be experimentally discernible. However, as one allows for larger , deviations in also begin to increase. Tables similar to Tables 2 and 3 for the channel are not presented here for the sake of brevity.
iii.4 channel singletop production  a Futuristic Analysis
We now return to the futuristic analysis that we briefly alluded to in the Introduction. As discussed in Sec. III.1, channel singletop production gets contributions from all processes of the type and where and . The new physics operators, however, only affect the subprocesses and . If one of these subprocesses could be isolated and studied separately, then it would be possible to obtain far more stringent constraints on the new physics parameters. The reason for the enormous increase in sensitivity is obvious – in the earlier cases, the NP contribution arising from just two subprocesses had to compete with the SM background arising from a multitude of subprocesses. In this case, it would be competing with background arising from just one subprocess. This also resolves the conundrum that we encountered in Sec.III.1 – if of (1) (or, equivalently 1 TeV) are not ruled out, why hasn’t new physics been discovered at the LHC ? It is because the sensitivity of the LHC in this context is limited by the fact that one or two NP amplitudes have to compete against a large number of SM amplitudes.
From the point of view of isolating channel subprocesses, appears at first glance to be more promising than , since events could be isolated from other channel events by the application of an additional tag. In reality, the situation is somewhat complicated.
In order to identify channel singletop production and distinguish it from channel singletop production, experiments already use an additional tag. The idea here is that if a top quark is produced in an channel process, then about 99% of the time, it is accompanied by a bottom quark owing to the strength of . In contrast, a final state is rare in the channel process.^{8}^{8}8A final state is approximately times more likely to occur than a final state in channel singletop production. Since tags do not distinguish between and , both and final states are identified as coming from the channel process. The contamination in channel measurements arising out of such misidenfications is not significant, since the cross section for the (channel) final state is orders of magnitude smaller than the dominant channel contribution.
However, if it were possible to distinguish between and quarks, then one would be able to isolate the process . To estimate the expected improvement in the limits, one only needs to consider the size of the SM background, which, at 13 TeV, is approximately 135 pb for the usual channel production and about 0.014 pb when is isolated. The actual limits are depicted in Fig. 13, which can be compared with Fig. 4 and with the 13TeV plots in Figs. 57. is depicted as a grey band in each of these cases. In Fig. 13(a) and (b), the effect of the interference between the SM and the NP term is discernible, unlike in the corresponding plots in Sec. III.1.
While isolating seems to be an exciting possibility, techniques for distinguishing between flavored quarks and antiquarks in the final state are still at a nascent stage, although some developments in this direction have been reported bcharge (). If such techniques can be improved upon sufficiently so as to become reliable even when the statistics are relatively low, then the channel will become the primary channel of interest.
Iv Summary and Conclusions
Top quark decays are most sensitive to new physics effects at the energy scale of a few hundred GeV. The same effects, if they arise from a higher energy scale, would be more effectively probed in singletop production. In this work, we have focused on NP effects arising from anomalous couplings between the top, bottom and charm quarks. Since we only consider the top’s interactions with heavy quarks, it might at first appear that progress would be thwarted by the low densities of heavy quarks inside the proton. Nevertheless, our detailed study shows that it would be possible to place meaningful constraints on the new physics parameters.
We have considered channel and channel singletop production cross sections to obtain constraints on contact interactions involving , and quarks. Of these two channels, the stronger constraints arise from the channel. This is due to the fact that the Standard Model background cross section is smaller for the channel. Within a given channel ( or ), the limits are most stringent for tensor operators, followed by vector and scalar operators, respectively. This is essentially due to the additional numerical factors that appear from the Dirac traces for each of these operators.
Apart from the total cross section, we have also considered the relative contributions to the cross section from top quarks of different helicities. The polarization asymmetry , which compares the helicity states of the top quark, can be particularly useful in establishing the presence of operators involving , especially since the corresponding Standard Model charged current coupling involves only .
Finally, we have considered a futuristic analysis. This analysis rests on one crucial assumption, namely that quarks in the final state can be distinguished from antiquarks on an eventbyevent basis. Adopting this assumption, we have obtained limits on the NP fourquark operators that are far more stringent than those obtained from the regular channel and channel analyses. The comparison can be seen in Fig. 14. The improvement is indeed startling and perhaps adequate motivation for the pursuit of experimental techniques that will make it possible.
ACKNOWLEDGEMENTS
The authors wish to thank David Cusick and Kristian Stephens for collaboration at early stages of this work, and David London and Bhubanjyoti Bhattacharya for helpful discussions. This work has been partially supported by CONICET and ANPCyT projects PICT20132266 and PICT20160164 (A.S.). K.K. was supported by the U.S. National Science Foundation under Grant PHY1215785. The work of P.S. was supported by the Dept. of Science and Technology, India under the INSPIRE Faculty Scheme Grant IFA14PH105.
Appendix
In this appendix we write down partonlevel expressions for the matrix element squared for the  and channel processes considered in this paper. Throughout, we use the notation and . Also, the spin fourvector for the top quark may be written as st_refs ()
where and represent the top quark’s threemomentum and energy (in a given reference frame), and where is a unit vector defined in the rest frame of the top quark. Note that , as expected. The symbol denotes that spin and color are averaged for the inital state and summed for the final state, except for the spin of the finalstate top quark. In each case,
Finally, we adopt the shorthand notation , taking .
Appendix A :
Appendix B : DcKs ()
Appendix C :
References
 (1) T. A. Aaltonen et al. [CDF and D0 Collaborations], [arXiv:1709.04894].

(2)
K. Kiers, T. Knighton, D. London, M. Russell, A. Szynkman and K. Webster,
Phys. Rev. D 84, 074018 (2011),
[arXiv:1107.0754];
K. Kiers, P. Saha, A. Szynkman, D. London, S. Judge and J. Melendez, Phys. Rev. D 90, 094015 (2014), [arXiv:1407.1724];
P. Saha, K. Kiers, D. London and A. Szynkman, Phys. Rev. D 90, 094016 (2014), [arXiv:1407.1725];
P. Saha, K. Kiers, B. Bhattacharya, D. London, A. Szynkman and J. Melendez, Phys. Rev. D 93, 054044 (2016), [arXiv:1510.00204]. 
(3)
T. Han, M. Hosch, K. Whisnant, B. L. Young and X. Zhang,
Phys. Rev. D 58, 073008 (1998),
[hepph/9806486];
T. M. P. Tait and C.P. Yuan, Phys. Rev. D 63, 014018 (2000), [hepph/0007298];
J. A. AguilarSaavedra, Nucl. Phys. B 804, 160 (2008), [arXiv:0803.3810]. 
(4)
R. Goldouzian and B. Clerbaux,
[arXiv:1612.01795];
D. Pinna, A. Zucchetta, M. R. Buckley and F. Canelli, Phys. Rev. D 96, 035031 (2017), [arXiv:1701.05195];
X. Y. Jiang, J. Z. Han, G. Yang and C. X. Yu, Commun. Theor. Phys. 67, 415 (2017);
I. Ahmed, [arXiv:1711.08348];
P. Pani and G. Polesello, [arXiv:1712.03874]. 
(5)
P. M. Ferreira, O. Oliveira and R. Santos,
Phys. Rev. D 73, 034011 (2006),
[hepph/0510087];
P. M. Ferreira and R. Santos, Phys. Rev. D 73, 054025 (2006), [hepph/0601078];
P. M. Ferreira and R. Santos, Phys. Rev. D 74, 014006 (2006), [hepph/0604144];
R. A. Coimbra, P. M. Ferreira, R. B. Guedes, O. Oliveira, A. Onofre, R. Santos and M. Won, Phys. Rev. D 79, 014006 (2009), [arXiv:0811.1743];
J. A. AguilarSaavedra, Nucl. Phys. B 843, 638 (2011), Erratum: [Nucl. Phys. B 851, 443 (2011)], [arXiv:1008.3562];
G. Durieux, F. Maltoni and C. Zhang, Phys. Rev. D 91, 074017 (2015), [arXiv:1412.7166];
J. A. AguilarSaavedra, C. Degrande and S. Khatibi, Phys. Lett. B 769, 498 (2017), [arXiv:1701.05900].  (6) G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 90, 112006 (2014), [arXiv:1406.7844].
 (7) M. Aaboud et al. [ATLAS Collaboration], Eur. Phys. J. C 77, 531 (2017), [arXiv:1702.02859].
 (8) The ATLAS Collaboration, ATLASCONF2011118.
 (9) G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 756, 228 (2016), [arXiv:1511.05980]; Phys. Lett. B 740, 118 (2015), [arXiv:1410.0647].
 (10) S. Chatrchyan et al. [CMS Collaboration], JHEP 1212, 035 (2012), [arXiv:1209.4533].
 (11) V. Khachatryan et al. [CMS Collaboration], JHEP 1406, 090 (2014), [arXiv:1403.7366].
 (12) V. Khachatryan et al. [CMS Collaboration], JHEP 1609, 027 (2016), [arXiv:1603.02555].
 (13) M. Aaboud et al. [ATLAS Collaboration], JHEP 1704, 086 (2017), [arXiv:1609.03920].
 (14) A. M. Sirunyan et al. [CMS Collaboration], Phys. Lett. B 772, 752 (2017), [arXiv:1610.00678].

(15)
LHC Top Working Group,
https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCTopWGSummaryPlots  (16) J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer and T. Stelzer, JHEP 1106, 128 (2011), [arXiv:1106.0522], http://madgraph.hep.uiuc.edu.
 (17) A. Alloul, N. D. Christensen, C. Degrande, C. Duhr and B. Fuks, Comput. Phys. Commun. 185, 2250 (2014), [arXiv:1310.1921], http://feynrules.irmp.ucl.ac.be.
 (18) J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky and W. K. Tung, JHEP 0207, 012 (2002), [hepph/0201195], http://hep.pa.msu.edu/cteq/public/cteq6.html.
 (19) A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63, 189 (2009), [arXiv:0901.0002], https://mstwpdf.hepforge.org.
 (20) N. Kidonakis, Phys. Rev. D 81, 054028 (2010), [arXiv:1001.5034].
 (21) N. Kidonakis, arXiv:1205.3453; arXiv:1311.0283.

(22)
S. Chatrchyan et al. [CMS Collaboration],
Phys. Rev. Lett. 112, 182001 (2014),
[arXiv:1311.3924];
V. Khachatryan et al. [CMS Collaboration], Phys. Rev. D 93, 052007 (2016), [arXiv:1601.01107];
G. Aad et al. [ATLAS Collaboration], Phys. Rev. Lett. 111, 232002 (2013), [arXiv:1307.6511].  (23) V. Khachatryan et al. [CMS Collaboration], JHEP 1604, 073 (2016), [arXiv:1511.02138].
 (24) V. M. Abazov et al. [D0 Collaboration], Phys. Rev. Lett. 98, 041801 (2007), [hepex/0608044]; G. Aad et al. [ATLAS Collaboration], JHEP 1311, 031 (2013), [arXiv:1307.4568].

(25)
C. Itzykson and J. B. Zuber, Quantum Field Theory, McGrawHill (1980),
ISBN 0070320713;
K. Kiers, A. Szynkman and D. London, Phys. Rev. D 74, 035004 (2006), [hepph/0605123].  (26) David Cusick and Kristian Stephens, private communication.