Top quark electroweak couplings at future lepton colliders

We perform a comparative study of the reach of future \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-$$\end{document}e+e- collider options for the scale of non-resonant new physics effects in the top quark sector, phrased in the language of higher-dimensional operators. Our focus is on the electroweak top quark pair production process \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^- \rightarrow Z^*/\gamma \rightarrow t\bar{t} $$\end{document}e+e-→Z∗/γ→tt¯, and we study benchmark scenarios at the ILC and CLIC. We find that both are able to constrain mass scales up to the few \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {TeV}$$\end{document}TeV range in the most sensitive cases, improving by orders of magnitude on the forecast capabilities of the LHC. We discuss the role played by observables such as forward–backward asymmetries, and making use of different beam polarisation settings, and highlight the possibility of lifting a degeneracy in the allowed parameter space by combining top observables with precision Z-pole measurements from LEP1.


Introduction
Studies of the couplings of the top quark form one of the cornerstones of the analysis program pursued at the LHC, where top quark related processes, most notably tt production, are ubiquitous. Whilst improving the precision of Standard Model predictions is important in its own right, a key motivation for the detailed study of the top quark is to understand its potential role in the physics behind electroweak symmetry breaking, which has yet to be elucidated. As the only fermion with a Yukawa coupling y t 1, the top is a pertinent place to look for the signatures of new physics at the TeV scale. Motivated scenarios range from compositeness theories (e.g. [1][2][3]) all the way to supersymmetry (see [4] for a review).
If the new physics that couples to the top quark is heavy and/or weakly coupled, an appropriate formulation of the effects of new physics is to view the SM as an effective thea e-mail: christoph.englert@glasgow.ac.uk b e-mail: m.russell.2@research.gla.ac.uk ory, where the non-standard couplings are parameterised by operators of dimension D > 4. The leading effects for collider observables typically enter at D = 6, L eff = L SM + i c (6) i O (6) i 2 + · · · (I.1) It is also typically justified to truncate the EFT expansion 1 at this order. Global fits of top quark couplings to data from ATLAS and CMS are already well established (see e.g. Refs. [6][7][8][9][10][11][12] for a representative sample as well as [13] for a discussion of electroweak precision effects), and a similar effort is well under way in Higgs physics, where the associated limits on Wilson coefficients are typically statistics limited. Still, despite the impressive statistical sample of top quark data that enters these fits, the subsequent direct bounds on the operator coefficients, and, by extension, the scale of new physics that would generate those operators, are rather weak [7]. There are few top quark measurements at the LHC that can be considered "precision observables" (helicity fractions in top decays are an exception [14]). By dimensional analysis, the strength of the interference of these operators with energy √ s typically scales as s/ 2 , and the areas of phase space that are most sensitive are plagued by correlated experimental and theoretical systematics. Moreover, the associated weak limits translate into values of that are probed by the high energy bins of the measurement, bringing into question the validity of the truncated EFT description [5] and care needs to be taken when combining measurements of different exclusive energy ranges of a binned distribution [15]. Inclusive cross-sections, being typically dominated by the threshold region √ s ∼ (2)m t , are under more theoretical control, but bring far less sensitivity.
Lepton colliders are not vulnerable to either of these problems. Firstly, there is excellent control over the hard scale of the interaction √ s, so one can always ensure that the limits on the D = 6 operators are consistent with a wellbehaved EFT expansion. 2 Secondly, the theoretical uncertainties from Standard Model calculations are much smaller: there are no PDFs, and the current state of the art precision for tt production is N 3 LO QCD at fixed-order [16], and NNLO+NNLL including threshold resummation, which bring SM scale uncertainty variation bands to the percent level [17].
The physics case for a e + e − collider is by now well established. The principal motivation is to perform a detailed precision study of the couplings of the Higgs boson in the much cleaner environment that a lepton collider affords, which will bring Higgs coupling measurements to an accuracy that will not be challenged by the LHC, even after it collects 3 ab −1 of data [18,19].
The electroweak couplings of the top quark are also clearly within the remit of such a collider. Currently, the only handle on top quark electroweak couplings from the LHC is through the associated production pp → tt V where V ∈ {Z , W, γ }. Whilst measurements of these processes are now approaching the 5σ level, the pull that they have on a global fit is small [7]. Measurements of electroweak single top production bring stronger bounds, but they are sensitive to a smaller subset of operators.
At a lepton collider, on the other hand, the process e + e − → Z * /γ → tt is extremely sensitive to top electroweak couplings. While the overall rate is more modest than at the LHC due to the parametric α EW /α s and s-channel suppression, the process is essentially background free, and it would constitute the first true precision probe of the electroweak sector of the top quark, and open up a new avenue for top quark couplings, complementary to the well-studied top QCD interactions.
Several studies of the prospects for improvement of top measurements at future colliders have already been undertaken (see for example Refs. [20][21][22][23][24][25]), in particular for the proposed International Linear Collider (ILC), but none have explicitly quantified the gain in the constraints on the top electroweak sector of the SMEFT, nor provided a comparative study of different collider options. The purpose of this paper is to provide such a study.
In Sect. 2 we discuss the e + e − → Z * /γ → tt process in the Standard Model EFT, and the D = 6 operators that generate interference with the Standard Model. To motivate the study, in Sect. 4 we discuss the present status of constraints on top electroweak couplings, and discuss projections for the lifetime of the LHC. The rest of the paper is devoted to lepton collider projections. In Sect. 5 we discuss ILC constraints, based on the H-20 running scenario. In Sect. 6, we focus on prospects from CLIC, comparing the potential of the two future collider projections. In Sect. 7 we analyse the complementarity between the bounds derived from tt measurements at a future collider and Z -pole measurements, before summarising in Sect. 8.

Top electroweak couplings
In the Standard Model, the electroweak tt Z coupling is given by the vector-axial-vector coupling To capture effects beyond the SM in this Lagrangian there are two approaches: one can write down anomalous couplings for the tt Z vertex, such that L ttZ receives a term where q = p t − pt . While this has the advantage of elucidating the various spin structures that can impact the tt Z vertex, it has the drawback that it does not allow for a simple power counting of which anomalous couplings would have the strongest effect. For example, the coefficient C 2 A is zero in the Standard Model, so that any corrections to it come solely from new physics contributions, which should be smaller than couplings that have SM interference.
To augment this description, one can instead supplement Eq. (II.1) with higher-dimensional operators. At leading order in the SMEFT, the list of operators that generate modifications to the tt Z vertex is, expressed in the basis and notation of Ref. [26]: The dictionary between the D = 6 operators of Eq. (II.4) and the anomalous couplings of Eq. (II.3) is where the superscript 33 denotes that we are considering the third generation only in the fermion bilinears of Eq. (II.4). Since c (3) ϕq and c (1) ϕq only appear with an overall opposite sign, we can only constrain the operator O (3) ϕq − O (1) ϕq ≡ O ϕq from tt Z couplings. We will discuss a method for bounding the two operators independently later in the paper.
C 2 A is generated by a CP-odd combination of operators, therefore it does not interfere with SM amplitudes and so its effects are expected to be smaller. Since in this study we are more interested in the absolute mass scales of these operators, we set all Wilson coefficients to be real, 3 however, we note that CP-sensitive observables such as angular distributions can also distinguish the CP character of the Wilson coefficients. We also assume that the new physics solely impacts the tt Z vertex, so we do not consider operators which modify the Zeē vertex, nor four-fermion operators which can contribute to the qq → tt or e + e − → tt processes, though we note that including these would in general weaken our bounds (see e.g. Refs. [7,11] for constraints on the former and Ref. [13] for the latter).

Analysis setup
We implemented the operators of Eq. (II.4) in a Feyn-Rules [28] model file, which was exported via UFO [29] to MadGraph [30], in order to generate parton-level events. For the analytic results discussed later, the model file was exported to FeynArts/FormCalc [31,32], which we used to calculate the full dependence of each observable used on the operators of Eq. (II.4). To generate lepton collider constraints, we minimise the χ 2 between the full 5D analytic expression and SM pseudodata using iMinuit, a Python/C++ implementation of the MINUIT [33] algorithm. For hadron collider constraints, we approximate the full 5D dependence of the cross-section on the operators using Professor; an interpolation-based method described in Ref. [34]. We crosscheck our analytic results against Monte Carlo scans using a UFO implementation of the full basis of Ref. [26] and find good agreement.

Hadron collider constraints
To appreciate the impact of the operators of Eq. (II.4), in Fig. 1 we plot the ratio of the full tt Z cross-section with each operator switched on individually, to the NLO SM prediction, taken from Ref. [35]. For ease of interpretation, we split up the cross-section into the contribution from the interference term and the quadratic term. We see firstly that the operators O uW and O u B have the strongest impact on the total crosssection, but this comes purely from the squared term (this was also noted in Ref. [8]). The remaining operators have a milder effect on the cross-section, but their interference term dominates. We also see that the operators O The LHC bounds on the coefficients of these operators from 8 TeV tt Z production cross-sections have been presented in Ref. [7]. The current constraints are weak. Since then, ATLAS and CMS have presented measurements using 13 TeV collision data, with measured values 0.9 ± 0.3 pb [36] and 0.7 ± 0.21 pb [37], respectively. The constraints on the operators using these two measurements are shown in Fig. 2, where the coefficients are normalised to the 'bar' notationc i = c i v 2 / 2 , and the operators are switched on individually.  We see that the current constraints are still quite weak, mainly due to the large (∼30%) experimental uncertainties. These measurements are currently statistics dominated, so it is instructive to ask what the expected improvement is over the lifetime of the LHC. Using a constant systematic uncertainty of 10% based on the current estimate, we also plot in Fig. 2 the constraints using 300 fb −1 and 3 ab −1 of SM pseudodata. We see that there will be an improvement by factors of 1.5 to 2 by the end of Run III, but after this the measurement is saturated by systematics.
To highlight the benefits of improving the theory description in tandem, we also show in Fig. 2 the projected constraints if theory uncertainties are improved to 1% from the current O(10%) precision, which does not seem unreasonable over the timescales we are considering. We see again that there will be no subsequent improvement after 300 fb −1 unless experimental systematics are reduced.
Finally, it should be noted that as more data becomes available, it may be possible to measure tt Z cross-sections differentially in final state quantities. Since cuts on the final state phase space can enhance sensitivity to the region where naïve power counting says D = 6 operators become more important, differential distributions could substantially improve the fit prospects, as has already been demonstrated for tt production. [7,38].
To illustrate this, in Fig. 3 we plot the distributions for the Z boson transverse momentum and top pair transverse momentum, both for the SM only case and with each operator switched on to a value ofc i 0.3; approximately the maximum allowed by current constraints in Fig. 2. We see that extra enhancement in the tail is visible for the field strength tensor operators O uW and O u B , due to the extra momentum dependence in the numerator from the field strength tensor. For the ϕ-type operators, since the interference is solely proportional to ϕ † ϕ → v 2 / 2 , there is no extra enhancement at high p T .
We do not estimate the improvement of the fit by taking these distributions into account, since this would require proper estimates of experimental systematics and tracking the nontrivial correlations between the kinematic quantities in the massive 3-body final state. Here, we merely comment that it may be an avenue worth pursuing as more data becomes available.

ILC constraints
Going beyond the LHC, currently, the most mature proposal is for a linear e + e − collider with a centre of mass energy ranging from 250 GeV to up to 1 TeV. There are several scenarios for integrated luminosity and CM energy combinations. The most-studied is the so-called H-20 option, which involves running at 500 GeV for 500 fb −1 of data, followed by 200 fb −1 of data at the tt threshold to perform detailed measurements of the top quark mass, and 300 fb −1 of data at √ s = 250 GeV to maximise the machine's Higgs potential with high precision. After a luminosity upgrade, a further 3.5 ab −1 is gathered at √ s = 500 GeV, followed by another √ s = 250 GeV run at 1.5 ab −1 . Since we are most interested in the ILC mass reach for new physics, in this study we focus on the 500 GeV ILC running. 4 An important parameter for lepton colliders is the energy spread of electron and positron beams (see e.g. [39]). In order to estimate the effect on our results, we use the results of [39] to calculate the expected change in the cross-section by including the effects of initial state radiation, beam spread and beamstrahlung. We find that, for the typical beam profile, the associated uncertainty is not a limiting factor and we neglect these effects in the following.

The tt total cross-section
Top pair production has a more modest rate here than at a hadron collider. The state-of-the-art Standard Model calculations for (unpolarised) e + e − → tt production are at N 3 LO QCD in the narrow-width approximation [16,17], and at NLO QCD including off-shell effects (which have been demonstrated to be important [40]) and at NLO EW [41] (with partial NNLO results in Ref. [42]) and predict a cross- section σ 0.57 pb. The conventional scale variation gives a QCD uncertainty at the per-mille level. While this rate is more than a factor of a thousand smaller than at the 13 TeV LHC, the process is essentially background free. Thus, after even 500 fb −1 of data the statistical uncertainty will be approximately 0.2%, and so completely subdominant to the systematics. We can thus repeat the exercise of extracting the bounds on the coefficients of the operators of Eq. (II.4) using SM pseudodata. As a guide for the expected numerical constraints, we also plot the ratio of the total cross-section in the presence of the operators to the SM prediction, this time using the total (unpolarised) cross-section at the 500 GeV ILC. This is shown on the left of Fig. 4.
We see again that the operators O uW and O u B are the strongest, however, unlike the case of tt Z production the interference term dominates at small c i / 2 . The result of this is that there is a cancellation between the interference and quadratic terms at approximately c i / 2 −3 TeV −2 , leading to a SM-like cross-section and a second, degenerate minimum in the χ 2 . The constraints obtained from a one-ata-time fit of these operators to the SM pseudodata is shown in the red bars on the right of Fig. 6.
The operators O uW and O u B are very tightly constrained, due to their much stronger impact on the cross-section stem- GeV, the operator colour coding is the same as Fig. 1. We also show a 5% uncertainty band around the SM prediction, to estimate the expected constraints ming from the extra momentum dependence flowing through the vertex. The ϕ-type operators are more weakly constrained, but on the whole the constraints are typically 10-100 times stronger than for the LHC tt Z production projections in Sect. 4, which is unsurprising giving the difference in precision.
Individual constraints are less useful in practice, however. Firstly, in a plausible UV scenario that would generate these operators one would typically expect more than one to be generated at once, so that one-at-a-time constraints cannot be straightforwardly linked to a specific 'top-down' model. Secondly, there can in general be cancellations between different operators for a given observable that can yield spurious local minima and disrupt the fit. This would not be visible in the individual constraints, and so would obscure degeneracies in the operator set that could be broken by considering different observables. Therefore, we also consider constraints where we marginalise over the remaining three coefficients in the fit. These are shown in the blue bars on the right of Fig. 6. 5 We see that, with the exception of O uW and O u B , marginalising over the full operator set wipes out the constraints. This is because even for large values of coefficients, the pull that a particular operator has on the cross-section can easily be cancelled by another operator. We can conclude that, despite the impressive precision that can be achieved in extracting the cross-section, it has limited use in constraining new physics in a simultaneous global fit of several operators. It is worthwhile to make use of other measurements.

Polarised beams
One of the principal strengths of lepton colliders is that the polarisation of the incoming beams can be finely controlled, so that the relative contributions between different subprocesses to a given final state can be tuned. Moreover, because the dependence of top observables on the operators of Eq. (II.4) depends strongly on the initial state polarisation, varying the settings increases the number of independent measurements that can be used to place bounds in a global fit.
To emphasise this point, we study the forward-backward asymmetry, defined as where θ t is the polar angle between the top quark and the incoming electron, for three incoming beam polarisation settings: unpolarised beams, denoted (A t F B ) U ; a fully lefthanded initial polarised electron beam and fully right-handed polarised positron beam, denoted (A t F B ) L ; and vice versa, denoted (A t F B ) R . The SM predictions for these settings at tree 40, 0.37, 0.47}, which agree well with the full NNLO QCD estimates [43,44]. The dependence of these asymmetries on the operators of Eq. (II.4) is shown in Fig. 5.
We see that the dependence on the operators distinctively depends on the initial state polarisations. the departure from the SM prediction is now much stronger for the ϕ-type operators than the total cross-section, which should lead to a sizeable improvement in the final constraints.
To generate these constraints, we consider a global fit of the four operators to six observables: In extracting the constraints, we consider the more realistic ILC polarisation capabilities P e − = ± 0.8, P e + = ∓ 0.3, noting that the cross-section for arbitrary e + e − polarisations is related to the fully polarised one by [45,46] σ P e − P e + = 1 4 where σ RL is the cross-section for fully right-handed polarised electrons and fully left-handed polarised positrons and σ RL is vice versa (the σ RR and σ LL components vanish for p-wave annihilation into spin-1 bosons). Performing a χ 2 fit of the full analytic expression for each observable, using SM pseudodata with 1% experimental error bars (based on studies in Refs. [22,47]) and SM theory uncertainties of 1% (based on the calculations of Refs. [16,17,41,43,44]) the individual and marginalised constraints on these operators are shown in Fig. 6. At the level of individual operators, the constraints are not improved drastically by adding in asymmetry information. For the global fit, however, the constraints lead to much stronger bounds than for fitting to cross-sections [although the marginalisation typically weakens the overall constraints by a factor O(10)].
We see that the constraints are again much stronger for the field strength operators O uW and O u B , where the constraints are at the |c i | 10 −4 level for the individual constraints and |c i | 10 −2 for the marginalised case, corresponding to a mass reach of 10 TeV and 2.16 TeV, respectively, assuming c i 1. The weakest constraints are on the operators O (3) ϕq (O (1) ϕq ), which translate into bounds on of roughly 700 GeV.
While it is encouraging that the bounds are consistent with an EFT formulation, in the sense that √ s, the ILC mass reach for the scale of new physics that would generate these operators is still low. We note, however, that these bounds are on the conservative side, since other observables such as oblique parameters and LEP asymmetries contribute complementary information that will in general tighten them. To keep this fit self-contained, we postpone this discussion until Sect. 7.

CLIC constraints
The Compact Linear Collider (CLIC) project [48,49], with its larger maximum centre-of-mass energy √ s = 3 TeV, will be in a stronger position to discover the effects of some higherdimensional operators, whose effects naïvely scale with the CM energy as s/ 2 . There are two main running scenarios, but both envisage total integrated luminosities of 500 fb −1 at √ s = 500 GeV, 1.5 ab −1 at 1.4 or 1.5 TeV, and 2 ab −1 at 3 (a) (b) Fig. 7 95% confidence ranges for the operators we consider here, from CLIC running at √ s = 3 TeV, assuming 1% theoretical and experimental uncertainties, by fitting to cross-sections, asymmetries, and the combi-nation, with each operator considered individually (a) or in a 5D fit (b). To display both on the same plot, we scale the individual constraints up by a factor of 10, so that the bottom axis is actuallyc i /10 TeV. Again, we focus on the highest energy setting √ s = 3 TeV, to maximise discovery potential for non-resonant new physics through D = 6 operators.
Moving further away from the tt threshold, the total e + e − → Z * /γ → tt rate is smaller than at the ILC; at √ s = 3 TeV it is around 20 fb, which means for the total forecast integrated luminosity at this energy there will be a statistical uncertainty of 0.5%. A total experimental uncertainty of 1% may therefore be too optimistic an estimate once systematics are fully itemised. Nonetheless, for ease of comparison with the ILC figures, we take this as a baseline, and the corresponding constraints, using the same observables and beam settings, are shown in Fig. 7.
We see that, for the individual fit, CLIC constraints are of the same order of magnitude as the ILC ones. 6 Although the direct sensitivity to the operators is enhanced, we see that as we move away from the tt threshold, the interference effect of the ϕ-type operators is much smaller. This is not the case for the operators O uW and O u B , whose contributions stem mainly from the (D = 6) 2 term, as seen on the right of Fig. 4, which receives no suppression. Their individual constraints are close to the ILC values, indicating that energy scale is not the dominant factor driving these limits, but rather the theory and experimental uncertainties which saturate the sensitivity, which we do not vary.
For the more general marginalised fit, we see again that combining cross-section and asymmetry measurements will break blind directions in the fit, leading to much more powerful overall constraints. Unlike for the case of the ILC, however, care must be taken in interpreting these limits in terms of the mass scale of a particular UV model. The marginalised constraint |c ϕu | 0.05, for example, corresponds to a mass scale / √ c 1.1 TeV, which is less than the energy scale probed in the interaction, so that the constraint can only be linked to a particular model if it is very weakly coupled: g * 1.
7 Beyond e + e − → tt: precision observables Obviously the direct constraints that we have focussed on in this work do not exist in a vacuum and the interplay of direct and indirect sensitivity plays an important part in ultimately obtaining the best constraints for a given model (see [51,52]).
To put the expected direct constraints detailed above into perspective we analyse the impact of the considered operators on LEP precision observables. Note, that these Z resonance observables are sensitive to a plethora of other new interactions and a direct comparison is not immediately straightforward [51]. Nonetheless, there is significant discriminative power that is worthwhile pointing out, which we will discuss in the following.

Oblique corrections
The S, T, U parameters [54,55] (see also [56]) are standard observables that capture oblique deviations in the SM electroweak gauge sector from the SM paradigm [57,58] through The operators considered in this work modify these at the oneloop level through diagrams of the type shown in Fig. 8. Throughout we perform our calculation in dimensional regularisation.
The definitions of S, T, U , see [54,55,58], are such that in the SM all divergencies cancel when replacing the renormalised polarisation functions by their bare counterparts. The modifications of Fig. 8, however, induce additional divergencies due to the dimension 6 parts and the introduction of twopoint function counterterms is essential to obtain a UV-finite result; see also [52]. This leads to a regularisation scale μ R dependence of the dimension 6 amplitude parts after renormalisation. It is this part which we focus on as we choose the SM with a 125 GeV Higgs as reference point [53].
As can be seen from Fig. 9, oblique corrections constrain the range of the considered parameter (setting all other contributions to zero). Given that we deal with an effective theory, there is a dependence on the UV cut-off μ R and constraints can be sharpened by requiring that cut-off to lie well above the TeV scale. As such, the oblique corrections are explicitly sensitive to concrete UV completions, and care needs to be taken not to over-emphasise their impact.

Non-oblique corrections
A well-measured quantity at LEP is the Zbb vertex, which enters the prediction of the bottom forward-backward asymmetry A bb F B ; see e.g. [59]. Similar to the operators in Eq. (II.4), in the generic dimension six approach we can expect similar operators for the down-sector of the third fermion family. These will modify the interactions along the same lines as we focussed on above for the top sector. However, due to the different isospin properties, the bottom forward backward asymmetry is now sensitive to the sum c (3) ϕq + c (1) ϕq . This leads ϕq + c (1) ϕq . The exclusion contour is taken from DELPHI collaboration's Ref. [59] for the most constraining measurement at √ s = 91.26 GeV Fig. 11 Allowed 95% confidence regions for the Wilson coefficients c ϕu andc (1) ϕu obtained from combining the information from ILC tt asymmetries and cross-sections (dark blue) and LEP1 bb measurements (lighter blue) to a complementary constraint by the LEP forward backward asymmetry compared to the direct measurements in tt, Fig. 10.
Moreover, the constraints on c (3) ϕq +c (1) ϕq from A bb F B can be combined with the constraints on c (3) ϕq −c (1) ϕq to extract independent bounds on c (3) ϕq and c (1) ϕq . This is shown in Fig. 11. Care should be taken when interpreting these constraints individually, however. We are considering marginalised bounds for the ILC constraint but only one operator combination for the LEP bound. In general, other operators that we do not consider here will impact the Zbb vertex at tree level and in general weaken the bound. This serves as a useful visualisation, however, of the complementarity between past and future colliders in constraining these operators.

Summary and conclusions
Given the unsatisfactory precision of current probes of top quark electroweak couplings from hadron collider measurements, they must be a key priority in the physics agenda of any future linear e + e − collider. By parameterising non-standard top couplings through D = 6 operators, we have analysed the potential for the ILC and CLIC to improve the current precision of the top electroweak sector. Unsurprisingly, if experimental precision would match current estimates, and theory uncertainties can be brought to the same level, the current constraints can be drastically improved by both colliders, with associated bounds on the scale of new physics typically in the 1 TeV to few TeV range, depending on the assumed coupling structures of the underlying model. Using asymmetry measurements as well as cross-sections will be crucial to this endeavour, as will collecting large datasets with several incoming beam polarisations.
We have found that, unlike for the Higgs sector, the large increase in centre-of-mass energy at CLIC does not necessarily offer a competitive advantage over the ILC for bounding new top interactions by the operators we consider, and bounds on the operators that we consider are typically stronger at the latter, though in simultaneous 4D fits the difference is not striking. For some of the operators we consider, the bounds derived from CLIC fits correspond to mass scales smaller than the CM energy that we consider, which can call into question the validity of the EFT description, unless the CLIC sensitivity can exceed the expectations we forecast here.
By combining Z -pole measurements from LEP1 with tt measurements (and future improved electroweak precision measurements), one can in principle break degeneracies in the operator set and disentangle individual operators that could previously be only bounded in combinations. We showed this for the LEP forward-backward asymmetry, this could be improved by fitting other precision elec-troweak observables too. Care must be taken in interpreting the associated constraints, however, as both sets of measurements will in general talk to other operators for which there is no complementarity, and a more systematic approach taking into account EFT loop corrections would have be undertaken before these numerical bounds can be taken at face value.