NLO QCD corrections to off-shell top–antitop production with semi-leptonic decays at lepton colliders

The study of top-quark properties will be a central aspect of the physics programme of any future lepton collider. In this article, we investigate the production of top-quark pairs in the semi-leptonic decay channel in e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\textrm{e}}^{+}\textrm{e}^{-}$$\end{document} collisions, whose experimental signature is one charged lepton, jets, and missing energy. We present for the first time fiducial cross sections and differential distributions at next-to-leading-order accuracy in QCD for the full off-shell process. We find that the QCD corrections for the considered process are strongly dependent on the beam energies and range from few per cent up to more than 100%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100\%$$\end{document} (near threshold and above 1TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\,\text {TeV} $$\end{document}). We focus, in particular, on two scenarios: one close to threshold (365GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$365\,\text {GeV} $$\end{document}), dominated by top-pair production, and one at the TeV scale (1.5TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.5\,\text {TeV} $$\end{document}), for which irreducible-background contributions become relevant. An assessment of polarised-beam effects is also provided.


Introduction
What large-scale collider experiment will come after the end of the Large Hadron Collider (LHC) is currently an open question.At present, several options are being considered which include lepton colliders such as the International Linear Collider (ILC) [1][2][3], the FCC-ee [4], or the Compact Linear Collider (CLIC) [5].In all cases, the study of top-quark properties will play a central role in the physics programme of those facilities.
The main advantage of lepton colliders over hadron ones is the possibility to tune very precisely the centreof-mass (CM) energy of the experiment.Thus, one can perform a scan of energies covering the threshold region for the production of a pair of top quarks.This provides a very clean access to key properties of the top quark such as its mass and width [6].
On the experimental side, significant prospective work has been done [7][8][9] to estimate the potential gain in performing such measurements and to assess their experimental limitations.On the theory side, great efforts have been put in providing precise predictions using non-relativistic QCD and resummation techniques at threshold [10][11][12][13].Differential predictions including also the transition to the continuum described by fixed-order QCD have been obtained in Ref. [14].
Above threshold, several predictions at fixed order in QCD have been provided for the on-shell production of a top-antitop pair, i.e. e + e − → t t, reaching next-to-next-tonext-to-leading-order (N 3 LO) accuracy for the inclusive cross section [15,16] and next-to-next-to-leading-order (NNLO) accuracy at the differential level [17][18][19][20].For the off-shell top-antitop pair production with leptonic decays, i.e. e + e − → j b j b − + ν¯ ν , which is well-defined both below and above threshold, several next-to-leading-order (NLO) QCD predictions have been provided [21][22][23].Regarding electroweak (EW) corrections, NLO accuracy for the inclusive cross section has been achieved long ago [24][25][26] and later supplemented with O α 2 ISR effects [27].Recently, the QED ISR effects at NLL in collinear factorisation have been matched to NLO EW corrections for on-shell production [28].
It is worth emphasising that for off-shell predictions the fully leptonic final state has been usually considered in the literature, with the exception of some sensitivity and background studies relying on LO off-shell simulations in the lepton-plus-jets channel [29][30][31][32].

arXiv:2302.04188v1 [hep-ph] 8 Feb 2023
In particular, the NLO QCD corrections for the semileptonic final state, i.e. e + e − → j b j b j j µ + ν µ are still unknown. 1The lepton+jets channel has the advantage to possess a larger cross section owing to the larger W-decay branching ratio.It also allows to fully reconstruct the momenta of the top quarks.
In the present work, we fill this gap by computing for the first time NLO QCD corrections for the process e + e − → j b j b j j µ + ν µ .In particular, we discuss phenomenological results in the case where all final-state particles are well separated, which corresponds to a so-called resolved topology as opposed to the case where light jets are allowed to be clustered in a large-radius b jet (boosted topology).We provide cross sections and differential distributions for different CM energies.
A further advancement of this calculation concerns the implementation of the FKS subtraction scheme [34] in the Monte Carlo integration code MoCaNLO.Among others, the present calculation served to validate the implementation of the FKS subtraction terms for processes with only final-state soft and collinear singularities.
This article is organised as follows.In Section 2, the process under investigation is presented (Section 2.1), the input parameters and event selections are listed (Section 2.2), and several remarks are provided regarding our implementation (Section 2.3).Section 3 discusses numerical results for the fiducial cross section and differential distributions.Finally, in Section 4 the main results obtained are summarised.

Calculation details 2.1 Definition of the process
In the present work, we consider the production of a topantitop pair in e + e − collisions in the semi-leptonic decay channel, e + e − → j b j b j j µ + ν µ , at NLO QCD accuracy.All final-state particles (quarks and leptons) are considered massless, and no quark mixing is taken into account (unit CKM matrix).With such a choice, the NLO corrections of order O(α s α 6 ) are genuine QCD corrections to the leading-order (LO) EW [O(α 6 )] cross section, as the EW corrections to the LO interference [O(α s α 5 )] vanish thanks to colour algebra.The real corrections are made of all possible gluon emissions from any of the coloured particles.The virtual corrections consist in the interference of Born amplitudes with one-loop ones, which are obtained by the insertion of a gluon in the tree-level matrix elements.As illustrated in Figure 1 for the leading order, all possible non-resonant and off-shell contributions are accounted for.In the top row, on the left-hand side, the typical production of a pair of top quarks and their semileptonic decay is depicted.The middle diagram shows a Higgs-strahlung type contribution where the Higgs boson decays into a pair of W bosons and the Z boson into a 1 For the LHC, the semi-leptonic final state has been computed by some of us few years ago [33].
bottom-antibottom pair.The diagram on the right-hand side shows a contribution to the same final state that does not involve any resonant top or antitop quark or Higgs boson.In the second row, a tri-boson (left) and a single-top (right) contribution are shown.We do not consider initialstate-radiation (ISR) and beam-strahlung effects of QED type as we restrict ourselves to QCD corrections.

Input parameters and kinematic selections
The computation is carried out in the five-flavour scheme, therefore m b = 0 is assumed throughout.The on-shell weak-boson masses and decay widths are fixed as [35] M OS W = 80.379 GeV, and then converted into the pole values [36].The Higgsboson and top-quark pole masses are chosen as [35] M H = 125 GeV, While the Higgs-boson width is taken from Ref. [37], the numerical value of the top-quark width is obtained by applying relative QCD corrections from Ref. [38] to the LO top-quark width computed following Ref.[39].All unstable particles are treated within the complex-mass scheme [40][41][42][43].
The EW coupling constant α is computed within the G µ scheme [44] with the Fermi constant set to The running of the strong coupling α s is carried out at two loops using the Recola program [45], assuming α s (M Z ) = 0.118.Finally, the renormalisation scale is set to µ R = m t , and the scale uncertainty is obtained by varying µ R by a factor 2 up and down.
In the following, we consider e + e − collisions at several CM energies.In addition to a scan of the integrated cross sections between 300 GeV and 2 TeV, shown in Section 3.1, we focus on two particular CM energies.Specifically, in Section 3.2 we provide differential results for 365 GeV, i.e. the highest collision energy envisioned for the FCC-ee [4].A similar energy is planned for the first operating scenario of CLIC [8,46], targeting the production of t t pairs above threshold.In Section 3.3, we show differential results for 1.5 TeV, i.e. the second operating stage of CLIC [9].
In our setup, the jet-clustering is carried out with the k T algorithm [47] and a resolution radius R = 0.4 [7,8].Note that a generalised version of the k T algorithm has been used in the t t study in the fully-leptonic decay channel [23].The clustering algorithm is applied on partons (quarks and gluons) with a minimum angle of 0.7721 • , which corresponds to a rapidity of 5. We use the following selection cuts, which are inspired by event selections applied in CLIC and FCC-ee studies [7][8][9].In our calculation, the events are required to have: • a minimum missing transverse-momentum p T,miss > 20 GeV, which is defined as the transverse momentum of the neutrino; • a minimum transverse momentum p T > 20 GeV for the antimuon, the light jets, and the two b-tagged jets; • an angular acceptance of 10 • < θ < 170 • for the antimuon, the light jets, and the two b-tagged jets; • a minimum rapidity-azimuthal-angle distance between the antimuon and the jets, ∆R j , ∆R j b > 0.4; • an invariant-mass cut on the system formed by the two hardest visible light jets, the charged lepton, and the neutrino of M jjµ + νµ > 130 GeV.
We require at least two light jets that fulfil all requirements on the transverse momentum, the angular acceptance, and the rapidity-azimuthal-angle distance to leptons (visible jets).Out of these jets, we select the two hardest ones (according to transverse momentum) for the distributions shown below.The condition on the invariant mass M jjµ + νµ selects a kinematic region that excludes the Higgs-boson decay into two leptons and two jets.
The two b jets present in the final state can be associated with the leptonically decaying top quark (j ).This is achieved by finding the maximum of a likelihood function that is the product of two Breit-Wigner distributions (of the top and antitop quark), as done in Ref. [48].The likelihood function mimics the top-and antitop-quark propagators, assuming three-body decays after recombination (both at LO and at NLO QCD), and reads with p abc = p a + p b + p c .The combination of bottom jets {j b,i , j b,j } that maximises L ij defines the two bottom jets originating from the leptonic and hadronic top quarks.Note that in Eq. (5) all possible combinations of light jets and b jets are considered.This includes also light jets with a minimum angle of 0.7721 • that do not fulfil the transverse-momentum, angular, and rapidity-azimuthalangle distance requirements for visible jets.Note further that the neutrino momenta are extracted from Monte Carlo truth, assuming that the hardscattering CM energy is exactly the one of the e + e − collision, i.e. neglecting ISR and beam-strahlung effects.

Implementation and validation
To carry out the present calculation, we have employed the Monte Carlo program MoCaNLO.In the past, MoCaNLO has been successfully used for several topassociated computations at NLO QCD and/or EW accuracy at hadron colliders [33,[48][49][50][51][52][53][54].The present work is the first application of MoCaNLO to a lepton-collider process.The program uses phase-space mappings similar to those of Refs.[40,55,56] and has shown to be particularly efficient for NLO calculations for high-multiplicity processes (up to 2 → 8).The tree-level and one-loop matrix elements are obtained from Recola [45,57] using the integral library Collier [58].For the subtraction of infrared divergences, the original code relies on the Catani-Seymour subtraction scheme [59][60][61].For the present calculation, we have implemented the FKS scheme [34] following closely Refs.[62,63].
To validate our implementation of the FKS scheme, we have compared our results against those obtained with the well-tested Catani-Seymour scheme for several NLO QCD calculations at lepton colliders including di-jet production, di-boson production in the semi-leptonic channel, off-shell top-antitop production in the fully leptonic channel, and the process considered in this work.In all cases, we have found perfect agreement within the Monte Carlo uncertainty at the level of both fiducial cross sections and differential distributions.For the √ s = 365 GeV setup considered in this article, the fiducial cross section obtained at NLO QCD with the two subtraction schemes reads, respectively, showing perfect agreement.Agreement has also been found at the differential level, where the numerical differences between the NLO distributions obtained with the two schemes are well within integration uncertainties bin by bin.We provide in Figure 2 a comparison of the differential results obtained with the two subtraction schemes for two observables, namely an angular one and a transverse-momentum one.The integration errors displayed in the lower inset of the plots are the combined integration errors of both NLO calculations.We observe that the bin-wise agreement is within the integration uncertainties over the whole spectrum.An analogous picture has been found for all other observables that we have computed.
In addition to the comparison against the dipole formalism, our implementation of the FKS subtraction scheme has allowed for further tests: • The cancellation of infrared poles in the n-body contribution to the NLO QCD cross section has been verified by evaluating the virtual contribution V and the integrated FKS counterterm I at different values of the infrared scale µ IR (from 10 −8 GeV to 10 8 GeV) and checking that the sum V + I is independent of the µ IR choice.This has been carried out for a large number of phase-space points finding agreement up to 12 digits.• The cancellation of phase-space singularities between the real matrix element and the FKS subtraction counterterm has been verified by constructing real-phasespace points that approach the soft, collinear, and softcollinear regions by means of a rescaling of the radiation variables.• The FKS-subtraction parameters ξ c and δ [34,62,63] which define the integration boundaries for the soft and collinear regimes, have been varied, confirming that the sum of the subtraction counterterm and its integrated counterpart is independent of them.Selected results regarding this subtraction test are shown in Table 1, where the reader can observe a rather strong impact of the soft parameter ξ c on the size and sign of subtracted real and virtual contributions, compared to a milder effect of the collinear one δ.• The evaluation of FKS sector functions and their sum rules have been checked in the subtracted-real contribution by means of a variation of the exponents a, b that enter the sector functions, which are defined in Eq. (5.11) of Ref. [63].For a fixed choice of the FKS parameters (ξ c = δ = 0.01) we have calculated the subtracted-real contributions for different choices of such exponents, e.g.
finding perfect agreement within integration errors.
The default values used for the results presented here are ξ c = δ = 0.01 and a = b = 1.

Fiducial cross sections
In this section, we report results for the fiducial cross section in the setup defined in Section 2. In Table 3, we provide the results at LO and NLO QCD accuracy for several choices of the CM energy, including regimes below and above the t t threshold.Since e + e − → j b j b j j µ + ν µ is a purely EW process at LO, there is no scale dependence at this order.At NLO, the QCD uncertainty comes from three-point variations of the renormalisation scale, µ R /µ The NLO QCD corrections are strongly dependent on the CM energy as can be seen in Figure 3, where the fiducial cross section is provided as a function of the CM energy in the range from 300 GeV to 2 TeV.The largest QCD corrections are observed slightly below the top-antitop threshold.In this regime, the presence of the Coulomb singularity renders the NLO QCD corrections divergent for on-shell top quarks [21,64], while the inclusion of decay effects makes them finite though still very large, reaching almost 170% of the LO cross section.In the case of the fully leptonic top-quark decays, a similar behaviour is observed at and around threshold [21,23].Above threshold the corrections turn negative.For the semi-leptonic process, the NLO QCD corrections become positive for energies above 1 TeV.For example at 1.5 TeV, the corrections are very large of the order of 60%.This is in contrast with the fully leptonic case where the corrections stay negative at high energies.This difference is due to the specific event selection.At LO, the cross section is suppressed in the semi-leptonic channel by the jet clustering (with R = 0.4) which effectively forbids boosted W bosons decaying into two quarks.At NLO QCD, in the presence of real gluon radiation, this constraint is lifted for sufficiently hard gluons.This part of the phase space therefore opens up and leads to relatively large corrections.This effect is specific to the semi-leptonic final state as in the fully leptonic case e + e j b j b jj + + X @ NLO QCD, unpolarised beams  e + e j b j b jj + + X @ NLO QCD, unpolarised beams  there are no cuts preventing boosted W bosons.We have verified this explanation by running the calculation for different jet radii.Smaller relative QCD corrections are found for a smaller jet-clustering radius (R = 0.1, allowing configurations with more boosted W bosons).For instance, at √ s = 1 TeV we find: Finally, we mention that the off-shell calculation embeds t t contributions as well as irreducible-background contributions that become more and more important with increasing CM energy, as observed in ILC sensitivity studies at 500 GeV [29,31,32].We have checked numerically that background contributions become indeed relevant.
In particular, we have found that single-top topologies [e.g. e + e − → bW − t, bW + t, as shown in Figure 1(e)] are the largest of these contributions.Tri-boson topologies [e + e − → W + W − Z, as shown in Figure 1(d)] also contribute but to a lesser extent.It is therefore interesting to realise that at very high energy the final state under investigation is not only made of top-antitop topologies but also of many others, rendering the reconstruction of top-antitop pairs difficult.

Differential distributions at 365 GeV
In this section, several differential distributions are presented at LO and NLO QCD accuracy.While the upper panels of the plots contain the absolute predictions, the lower ones show the corresponding K factors.In the following, when an observable refers to either the leptonically or the hadronically decaying top quark, their definition follows from the maximisation of the likelihood function in Eq. (5).
In Figure 4, several transverse-momentum distributions are shown.The first two are for the reconstructed bottom quark from the leptonically and hadronically decaying top quarks, respectively.At LO, both distributions are almost identical with a pronounced drop around 95 GeV.This results from a hard cut for production of on-shell top-quark pairs decaying into on-shell W bosons that can be evaluated to For √ s = 365 GeV, this amounts to p T,j b 94.4 GeV.At NLO QCD, the picture changes as the hadronic top quark receives significantly more corrections at high transverse momenta.This is related to additional real-radiation jets that are mis-reconstructed as top-decay jets.It is worth mentioning that the reconstructed bottom-quark distributions are very close to those obtained with Monte Carlo truth.For what concerns the reconstructed top quarks (leptonic and hadronic), both at LO and NLO QCD accuracy the two distributions show very similar qualitative behaviours.On-shell production of top-antitop quarks is restricted to p T,t < ( √ s/2) 1 − 4m 2 t /s ≈ 58 GeV, leading to a sharp drop of the distribution above this value.In the off-shell region, the NLO QCD corrections are somewhat larger for the hadronically decaying top quark.For p T,t 127 GeV the recoiling system cannot contain a resonant W boson anymore explaining the shoulder near 120 GeV.We note that in both cases, around 100 GeV and above, the top-transverse-momentum distributions become numerically unstable.Besides the lower statistics, this is simply due to the fact that at such energies, the process is not exclusively made of top-antitop topologies as explained in the previous section.
In Figure 5, several invariant-mass distributions are displayed.The first two are for the invariant masses of the leptonically [Figure 5(a)] and hadronically [Figure 5(b)] decaying top quarks, respectively.It is interesting to observe that, as for the transverse-momentum distributions, the LO predictions are essentially identical while at NLO they significantly differ.This is due to the fact that the hadronically decaying top quark possesses three partons in the final state as opposed to only one for the leptonically decaying one, leading to more final-state radiations in the hadronic case.As a consequence, more events are moved from the resonance or above to below the resonance owing to final-state radiation that is not reconstructed with the decay products of the top/antitop quark forming hence a large radiative tail (see, for instance, Refs.[33,65]).In the case of the hadronically decaying top quark, the effect is so large that the NLO cross section becomes negative above the resonance.Such a behaviour has already been observed for the same final state at a hadron collider [33] and requires the inclusion of higher-order corrections for a proper description of this observable.
A radiative tail also appears in the invariant-mass distribution of the two hardest light jets [Figure 5(c)], which at LO reconstruct a W boson. Again, the effects are extremely large with K factors reaching ten below the resonance.The distribution in the invariant mass of the system formed by the reconstructed bottom quark from the leptonically decaying top quark and the antimuon [Figure 5(d)] has been found to be very sensitive to the top-quark mass as it possesses an on-shell edge at 65,66].While the relative corrections are flat in the on-shell region, they strongly increase above the on-shell edge.
Finally, Figures 5(e)-5(f) show distributions in the cosine of the production angle of the reconstructed top quarks.Both distributions are relatively similar up to a reflection of the directions.Indeed, given that, as opposed to the LHC, the initial state is asymmetric, the top and antitop quarks have preferred directions while generally ending up in a back-to-back configuration.The NLO QCD corrections are flat and reproduce those of the fiducial cross section.

Differential distributions at 1.5 TeV
In Figure 6 we show the results obtained for a few selected observables in e + e − collisions at 1.5 TeV CM energy.We stress that the technique used in the 365 GeV analysis to reconstruct the top and antitop quarks is not performing well at 1.5 TeV, owing to the presence of sizeable irreducible backgrounds not involving a top-antitop pair in the fiducial volume.Therefore, we do not show topreconstructed observables in this section as their physical interpretation is unclear.The higher CM energy shifts the most populated transverse-momentum and invariant-mass regions to higher values than at 365 GeV and enhances irreducible backgrounds that are suppressed for CM energies around the top-pair threshold.
In Figure 6(a) we show the differential results in the transverse momentum of the second-hardest light jet.At LO and for on-shell W bosons, this observable is characterised by a kinematic cut-off which, assuming small angles between the two jets, is given by p T,j 2 ,max ∼ m jj,max /∆R jj,min ∼ M W /0.4 ∼ 200 GeV [33].The NLO QCD corrections, rather flat in the most populated region (p T,j 2 150 GeV), fill the kinematic regime that is suppressed at LO, driven by real corrections with a third jet from gluon radiation that are tagged as the second-hardest jet.The light-jet kinematics is strongly affected by QCD corrections at NLO, as can also be observed in Figure 6(b) where the invariant mass of the two hardest light jets is considered.The typical Breit-Wigner shape coming from the LO W-boson hadronic decay is distorted by QCD corrections both below and above the W-boson pole mass.The QCD corrections at the peak are negative and qualitatively similar to the ones at 365 GeV in Figure 5(b).The radiative tail below M W is similar to the one observed at 365 GeV, while the enhanced K factor found for M jj 200 GeV originates from hard gluon radiation.
The distribution in the invariant mass of the b-jet pair is considered in Figure 6(c).A clear peak at M Z highlights the contributions with a resonant Z boson decaying into two b jets, produced in association with two W bosons.
The Z-boson peak sits on top of contributions coming both from the t t-production process and other sizeable backgrounds as single-top production.The QCD corrections are large and negative at the Z-boson peak, similarly to what can be observed in Figure 6(b) for the hadronic Wboson decay.Below M Z , the QCD corrections are large and positive, due to the LO suppression and the presence of unclustered gluon radiation.Above M Z the corrections are positive and diminish in size from 100% to a few percent around 1 TeV.For an invariant mass close to the maximal possible value for on-shell production of about 1170 GeV [twice the transverse momentum resulting from Eq. ( 10)] the leading order is suppressed and the K factor increases.
Rather flat relative QCD corrections are found for the distribution in the invariant mass of the bottom-antimuon system, shown in Figure 6(d).As in Figure 5(d), the bottom momentum is obtained from Monte Carlo truth.The edge at m2 t − M 2 W observed in the 365 GeV scenario is present also in the high-energy scenario, but the drop of the LO cross section around this threshold is less severe at TeV-scale energies due to the increased irreduciblebackground contributions which do not embed a resonant top quark.This leads to a QCD K factor that does not increase for M µ + b > 153 GeV.

Polarised-beam effects
In the baseline FCC-ee scenarios, the beams are planned to be unpolarised [4,67], while the CLIC and ILC facilities are envisoned to collide 80%-polarised electrons and possibly 30%-polarised positrons [1,5].It has been claimed [3] that polarised beams at lepton colliders are beneficial to enhance the sensitivity to EW parameters and possible new-physics effects, increase the signal-to-background ratio for several signatures, and keep systematics under control.Assuming a partial polarisation along the beam axis (often dubbed longitudinal polarisation) for both the positron (P e + ) and the electron (P e − ), the cross section for a given process reads, + (1 where σ XY is the cross section for a positron with helicity X and an electron with helicity Y , and L(R) stands for left(right)-handed helicity.Note that in annihilation processes, like the one we consider in this work, the Standard-Model dynamics only allows for a combined angular momentum equal to 1. Therefore, the initial-state leptons cannot carry the same helicity, i.e. σ LL = σ RR = 0.
In Table 2 we show the integrated cross sections for a number of beam-polarisation choices and for both energy scenarios considered in this work.Close to the threshold, the configuration with a right-handed positron and a left-handed electron (RL) gives a fiducial cross section approximately 2.5 times larger than the opposite helicity configuration (LR).The QCD corrections relative to the corresponding LO cross section are almost identical in all pure or mixed helicity configurations, in agreement with the results of Ref. [23] for undecayed top quarks.At high energy, the ratio between the RL and LR cross sections sizeably increases at LO, being ≈ 4 at 1.5 TeV, while the relative NLO QCD corrections are larger for the LR state (+80%) than for the RL one (+50%).This difference between the two pure polarisation states is driven by the real-radiation contributions which open up new helicity configurations in the final state that are suppressed at LO.When irreducible-background effects become relevant, the trivial factorisation of QCD corrections from the initial-state helicity configuration does not hold anymore [23].
The effects of beam polarisation are maximal in differential distributions, mostly for angular observables.In Figures 7-9 we show differential results at fixed initialstate helicities (LR, RL) for both the 365 GeV and the 1.5 TeV energy scenarios.Differential results for partially polarised beams, though not considered here, can be estimated easily performing a bin-by-bin combination of the RL and LR distributions according to Eq. (11).At energies well above the t t threshold, selecting the RL helicity state is expected to enhance the ratio of single-top and non-resonant contributions over the top-antitop ones [30].
In Figure 7 we consider the distribution in the cosine of the polar angle of the second hardest light jet.The LR and RL shapes at 365 GeV are related by an almost perfect mirroring about cos θ j 2 = 0, up to the different overall normalisation.This holds both at LO and at NLO QCD.In the RL state, the antitop quark typically goes forward, as shown in Figure 5(f), and is mostly right handed [68].According to the helicity structure of the tree-level on-shell top-decay amplitude, the W − boson from the decay of a mostly right-handed antitop quark is produced backward with respect to the antitop direction (in the antitop rest frame), therefore giving light jets that are typically produced with cos θ j < 0. In fact, this is the case for the second-hardest jet in Figure 7(a) but also for the hardest jet, though with different distribution shapes.The same reasoning applies to the LR initial state, with a flip of sign in the left-right asymmetry in the decay matrix element, therefore motivating the almost exact mirror symmetry between the LR and RL states.The QCD effects are slightly different for the LR and RL states, but in both cases they are larger where the corresponding LO are more kinematically suppressed.The situation is significantly different at 1.5 TeV.The RL state features asymmetric peaks at the distribution endpoints, the LR one peaks in forward regions, while being suppressed in backward regions.This results from the fact that at high energies the antitop quarks have high energies, while still going mostly in the forward direction for the RL state.As a consequence, their decay products also often end up in forward direction.The QCD corrections mostly fill the region with negative cos θ j 2 for both helicity states.Especially for the LR state, the QCD K factor is between 5 and 6 in the region that is mostly suppressed at LO. Shape differences between the LR and RL helicity states are not observed for all angular observables.In Figure 8 we show distributions in the cosine of the polar angle of the hardest b jet.For this variable, the LR and RL distributions are characterised by very similar LO shapes and QCD K factors at 365 GeV, with the most populated region being the central one.At 1.5 TeV the central region is suppressed, while the forward-backward ones are favoured by both the t t and the single-top contributions. 2 The inclusion of QCD corrections enhances mostly the suppressed central regions.In spite of a much  larger K factor for the LR state, the NLO shapes are quite similar for the two pure helicity states.The rationale is that an approximate factorisation of QCD corrections w.r.t. the beam helicity is fulfilled in the fiducial volume for a beam energy around the top-mass threshold, while it is broken for TeV-scale CM energies mostly due to the LO suppression of some helicity configurations and the opening of new ones thanks to gluon emission.The difference between different helicity selections can be appreciated not only in angular distributions, but also in the invariant mass of the b-jet pair considered in Figure 9.Note that we focus here on masses below 200 GeV for both CM energies, while the unpolarised distribution in Figure 6(c) clearly shows that the relevant and forward directions and thus at high energies also their decay products.
range for the 1.5 TeV scenario is much larger.The contributions involving the Z decay into two b jets give a peak at the Z mass that is sizeable only in the RL shape, while its presence is hardly visible for the LR state, since e + e − → W + W − Z contributes basically only for the former owing to the purely left-handed coupling of the W boson.While this holds for both energies, the tri-boson contribution to the cross section is small compared to the dominant t t-production process at 365 GeV, but its contribution is larger at 1.5 TeV.The difference between the LR and RL states at 1.5 TeV is propagated to the QCD corrections which are sizeable and negative at the Z peak for the RL state, while they are less pronounced for the LR one.On the contrary, at 365 GeV, the QCD K factors are almost identical for the two helicity states, apart from a narrow region around the Z peak and close to the kinematic  Same structure as in Figure 7.
boundary for on-shell top production M j b j b 189 GeV given by Eq. ( 10) and M j b j b < 2p T,j b .
As shown above, the beam-polarisation effects are relevant for some observables but completely irrelevant for others.These effects are in general enhanced for higher CM energies, where the interplay among different subprocesses enhances spin configurations in the final state depending on the polarisation of the incoming beams.In particular, the appearance of hard-gluon radiation at NLO can allow for different spin configurations than at LO.A due comment concerns the flavour of the final-state lepton.Since we have considered e + e − → j b j b j j µ + ν µ , with different lepton flavours in the initial and final state no t-channel or boson-fusion topology appears.If a positron was considered in the final state, contributions from the RR and LL initial-state configurations would be non vanishing, and more resonant structures would be present making the spin structure of the process even more involved.

Conclusion
Top-quark physics will play a central role in any of the potential future lepton colliders.It is therefore very important to provide precise theoretical predictions.So far, only the off-shell production of a top-antitop pair with fully leptonic decays had been computed at NLO QCD.
Nonetheless, the semi-leptonic channel offers many advantages such as the larger cross section and the possibility to fully reconstruct the momenta of the top quarks.In this paper, we have provided the first calculation of the NLO QCD corrections to the full process e + e − → j b j b j j µ + ν µ .To that end, we have implemented the FKS subtraction scheme in the Monte Carlo integrator MoCaNLO.We have successfully validated this implementation against the already existing Catani-Seymour scheme in our Monte Carlo program.
At the level of fiducial cross sections, the QCD corrections strongly depend on the collision energy, ranging from huge positive values at the top-pair threshold to negative values above threshold but lower than 1 TeV and increasingly positive values beyond 1 TeV.This dependence results from the Coulomb singularity in the threshold region and the selection of the decay jets that suppresses specific kinematic regimes at high energies at LO but not at NLO because of hard gluon radiation.
The behaviour of QCD corrections becomes even more striking at differential level, with huge K factors in regions where the LO cross section is suppressed and the real corrections open up new kinematic topologies.This holds for both collision energies we have considered (365 GeV, 1.5 TeV).The radiative corrections become large especially for invariant-mass and transverse-momentum distributions.The most extreme K factors are found at a 1.5 TeV collision energy where the QCD effects are enhanced far from the top-pair threshold.
The structure of LO contributions and QCD radiative corrections can be understood in more detail when selecting pure helicity states for the electron and positron beams.The factorisation of higher-order QCD effects with respect to the spin state of the leptonic initial state is typically confirmed at collision energies around the toppair threshold.It is clearly broken at 1.5 TeV, owing to enhanced irreducible background processes and the suppression of certain helicity configurations at high energy.

Fig. 1 .
Fig.1.Sample tree-level Feynman diagrams for the process e + e − → j b j b j j µ + νµ.
t lep b ) and the hadronically decaying antitop quark (j t had b

Fig. 3 .
Fig. 3. Fiducial cross sections for e + e − → j b j b j j µ + νµ at LO and NLO QCD at various CM energies √ s.Numerical values are shown in the table (left), where the digit in parentheses indicates the Monte Carlo statistical error, while the sub-and superscripts in per cent indicate the renormalisation-scale uncertainties calculated with three-point scale variations.The integrated cross section is shown in the figure (right) at LO (blue) and NLO QCD (red) as a function of the CM energy.The red-shaded band is obtained by means of three-point renormalisation-scale variations.

Fig. 4 .
Fig. 4. Differential distributions for a 365 GeV CM energy in the transverse momentum of the reconstructed bottom quark from the leptonically decaying top quark (top left), the reconstructed bottom quark from the hadronically decaying top quark (top right), the reconstructed leptonically decaying top quark (bottom left), and the reconstructed hadronically decaying top quark (bottom right).The red-shaded band is obtained by means of three-point renormalisation-scale variations and the lower panel displays the K factor.

Fig. 5 .
Fig.5.Differential distributions for a 365 GeV CM energy in the invariant masses of the reconstructed leptonically decaying top quark (a), the reconstructed hadronically decaying top quark (b), the system of the two hardest light jets and the bottom quark from the leptonically decaying top quark with the antimuon (d), as well as in the cosine of the production angle of the reconstructed leptonically decaying top quark (e) and the reconstructed hadronically decaying top quark (f).Note that for Figure5(d), the Monte Carlo truth momenta are used.Same structure as in Figure4.

Fig. 6 .
Fig. 6.Differential distributions for a 1.5 TeV CM energy: transverse momentum of the second hardest light jet (a), invariant masses of the system formed by the two hardest light jets (b), the two b jets (c), and the bottom and the antimuon (d).The bottom momentum used in Figure 6(d) is taken from Monte Carlo truth.Same structure as in Figure 4.

RLRLFig. 7 .
Fig.7.Differential distributions in the cosine of the polar angle of the second hardest light jet for CM energies 365 GeV (left) and 1.5 TeV (right).The top panels show the absolute differential cross sections at LO (blue) and at NLO QCD (red) for the RL (solid) and LR (dashed) helicity combinations of the initial-state leptons.The red-shaded band is obtained by means of three-point renormalisation-scale variations and the lower panel displays the K factor.

RLFig. 8 .
Fig. 8. Differential distributions in the cosine of the polar angle of the hardest b jet for CM energies 365 GeV (left) and 1.5 TeV (right).Same structure as in Figure 7.

RLRLFig. 9 .
Fig. 9. Differential distributions in the invariant mass of the b-jet pair for CM energies 365 GeV (left) and 1.5 TeV (right).Same structure as in Figure 7.

Table 1 .
Differential distributions in the cosine of the angular separation between the hardest light jet and the antimuon (left) and in the transverse momentum of the hardest b jet (right).The NLO QCD distributions obtained with the Catani-Seymour (dubbed CS) and FKS subtraction schemes are plotted in the top panel besides the LO results.In the bottom panel the ratio between the FKS and CS results (solid line) and the corresponding Monte Carlo integration error (shaded band) are shown.. Comparison of subtracted virtual (V subtr

Table 2 .
Fiducial cross sections for e + e − → j b j b j j µ + νµ at LO and NLO QCD for several beam polarisations and for CM energies 365 GeV and 1.5 TeV.The digit in parentheses indicates the Monte Carlo statistical error, while the sub-and super-scripts in per cent indicate the scale variation.