Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
aainstitutetext: Centro de Ciências Naturais e Humanas, Universidade Federal do ABC, Santo André, 09210-580 SP, Brazilbbinstitutetext: Instituto de Física Corpuscular (IFIC), Universidad de Valencia-CSIC, E-46980 Valencia, Spain

Going beyond Top EFT

André Lessa b    and Verónica Sanz andre.lessa@ufabc.edu.br veronica.sanz@uv.es
Abstract

We present a new way to interpret Top Standard Model measurements going beyond the SMEFT framework. Instead of the usual paradigm in Top EFT, where the main effects come from tails in momenta distributions, we propose an interpretation in terms of new physics which only shows up at loop-level. The effects of these new states, which can be lighter than required within the SMEFT, appear as distinctive structures at high momenta, but may be suppressed at the tails of distributions. As an illustration of this phenomena, we present the explicit case of a UV model with a 𝒵2subscript𝒵2\mathcal{Z}_{2}caligraphic_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT symmetry, including a Dark Matter candidate and a top-partner. This simple UV model reproduces the main features of this class of signatures, particularly a momentum-dependent form factor with more structure than the SMEFT. As the new states can be lighter than in SMEFT, we explore the interplay between the reinterpretation of direct searches for colored states and Dark Matter, and Top measurements, made by ATLAS and CMS in the differential tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG final state. We also compare our method with what one would expect using the SMEFT reinterpretation, finding that using the full loop information provides a better discriminating power.

1 Introduction

At the LHC, searches for new phenomena in the Standard Model Effective Field Theory (SMEFT) framework are now commonplace, see e.g. Refs. CMS:2023ixc ; CMS:2022hjj ; CMS:2022uox ; CMS:2020gsy ; CMS:2021cxr ; ATLAS:2022xyx ; ATLAS:2021jgw ; ATLAS:2023jma for recent experimental results. They provide a way to re-interpret Standard Model (SM) measurements which exploits their full kinematic range and can guide combinations of different channels. The most striking signatures of SMEFT show up at tails in energy-momentum distributions Ellis:2014dva , extreme kinematic regions where the SM contribution is scarce and the new phenomena more visible. A similar story can be told for light axion-like particles (ALPs), whose derivative couplings also induce prominent effects in the tails of distributions from SM measurements Gavela:2019cmq ; Folgado:2020utn .

On the other hand, direct searches for new physics are based on on-shell production of the new states, which can then decay, interact with the detector or escape detection. Direct searches can be based on signatures with very low SM background and/or searches for excesses in specific channels and phase space regions. For instance, in the case of bump-hunt searches, one would scan for deviations from a smooth SM background in a range of resonance masses. Nowadays, these searches are sensitive to very high masses (well above TeV) for traditional channels (e.g., dijet or dilepton) and the LHC experimental collaborations are continuously broadening the coverage for possible final states. Despite their impressive sensitivity, resonance searches have an intrinsic limitation: they make sense for narrow states, with widths (ΓΓ\Gammaroman_Γ) much smaller than their mass (m𝑚mitalic_m), and typically their performance quickly stops at Γ/mgreater-than-or-equivalent-toΓ𝑚absent\Gamma/m\gtrsimroman_Γ / italic_m ≳ 0.3, see e.g. Refs. CMS:2021ctt ; CMS:2019gwf for recent experimental analyses with variable widths 111A theoretical proposal to broaden the scope of these searches has been presented in Ref. Chivukula:2017lyk .. In addition, due to trigger requirements, soft final states which can appear in compressed Beyond the Standard Model (BSM) scenarios can also be very challenging for direct searches, which must rely on initial or final state radiation for triggering.

Somewhere in between the SMEFT and on-shell paradigms lies the proposal of this paper, namely the exploration of scenarios that are just beyond the reach of direct searches, but are not correctly described by the Effective Field Theory (EFT) limit. As we will show, these scenarios can be probed by SM measurements, but their signal can be very distinct from what would be expected from SMEFT. In particular, if the new states contribute to SM measurements through loop diagrams, their effect in the differential distributions would be localized in a kinematic region, resembling a very broad bump. In the limit that the states running in the loop are very heavy, this localization would shift towards high invariant masses, reaching the SMEFT limit.

As discussed in Ref. Cepedello:2022pyx , scenarios which contribute at loop-level at leading order are a good testing bed for the interplay between direct searches and indirect probes, as typically the new states can be lighter than in scenarios with tree-level contributions. And among the set of loop-induced UV models, those with a Dark Matter (DM) candidate are particularly interesting Cepedello:2023yao . Moreover, we will focus on scenarios with a special relation with the top sector, which will allow us to draw a comparison with the current efforts on the Top EFT searches and provide an alternative to those. A similar approach was considered in Ref. Drozd:2015kva , but within the context of Higgs couplings.

The paper is structured as follows. In Sec. 2 we review the main concepts in the Top sector of the SMEFT. In the section 3, we present a minimal loop-induced scenario with a DM particle and a heavy top-partner, discuss the analytical behaviour of their contribution to top observables, and explore the connection to the Top EFT. In Sec. 4 the limits from direct searches are reviewed, as well as the limits from precise SM top observables (mtt¯subscript𝑚𝑡¯𝑡m_{t\bar{t}}italic_m start_POSTSUBSCRIPT italic_t over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT and pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT distributions). Those direct and indirect probes are placed together and compared with the Top EFT limit in Sec. 4.3. Finally, in Sec. 5, we conclude. Auxiliary information concerning the loop calculation, matching to the EFT regime and the limit setting are given in Appendices A, B and C.

The datasets used to obtain all the results presented here as well as additional details are available in the Zenodo zenodoRepo and GitHub repositories.

2 Top EFT

One of the goals of this paper is to show a different way to interpret SM measurements, namely to search for new states which can be relatively light when compared to the LHC energy scale and only contribute to higher dimensional operators at the loop level. To illustrate this point, we will show results from a UV extension of the SM with a singlet fermionic Dark Matter candidate and a scalar colored state, a partner of the right-handed top, tRsubscript𝑡𝑅t_{R}italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT. For energy scales sufficiently smaller than the BSM masses, this scenario and its phenomenology can be matched to a reduced set of Top EFT operators.

The top EFT has been described in many works, e.g. Buckley:2015lku ; Etesami:2017ufk ; Maltoni:2019aot ; Brivio:2019ius ; Miralles:2021dyw ; Kassabov:2023hbm , and it is a subset of the SMEFT Lagrangian. In this paper we will use the Warsaw Grzadkowski:2010es convention to classify the independent operators. Note that if we allowed for the most general flavour structure, we would find that there are 2499 different types of operators which contribute at dimension-six Alonso:2013hga , but this number is drastically reduced once we assume some type of flavour structure in the UV completions, an assumption well motivated by the obstinate absence of anomalies in flavour observables. In particular, when focusing on the top physics, it is common to consider a flavour SU(3)5𝑆𝑈superscript35SU(3)^{5}italic_S italic_U ( 3 ) start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT symmetry, leading to a top-specific scenario. The details of this scenario, including a classification of the relevant operators and their limits from a global fit can be found in Ref. Ellis:2020unq . As we will show below, the particular UV model considered here induces a new interaction of two tops with gluons,

𝒪tG=(Q¯σμνTAtR)φ~GμνA,subscript𝒪𝑡𝐺¯𝑄superscript𝜎𝜇𝜈superscript𝑇𝐴subscript𝑡𝑅~𝜑superscriptsubscript𝐺𝜇𝜈𝐴\mathcal{O}_{\scriptscriptstyle tG}=(\bar{Q}\sigma^{\mu\nu}T^{A}t_{R})% \widetilde{\varphi}\,G_{\mu\nu}^{A},caligraphic_O start_POSTSUBSCRIPT italic_t italic_G end_POSTSUBSCRIPT = ( over¯ start_ARG italic_Q end_ARG italic_σ start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) over~ start_ARG italic_φ end_ARG italic_G start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT ,

where σμν=i2[γμ,γν]superscript𝜎𝜇𝜈𝑖2superscript𝛾𝜇superscript𝛾𝜈\sigma^{\mu\nu}=\frac{i}{2}\left[\gamma^{\mu},\gamma^{\nu}\right]italic_σ start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT = divide start_ARG italic_i end_ARG start_ARG 2 end_ARG [ italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT , italic_γ start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT ]. This gluon-top coupling modifier will be accompanied by a set of four-fermion interactions between the right-handed tops and the quark doublets (of any flavor), qLsubscript𝑞𝐿q_{L}italic_q start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT:

𝒪tq(8)=(q¯LγμTAqL)(t¯RγμTAtR),superscriptsubscript𝒪𝑡𝑞8subscript¯𝑞𝐿subscript𝛾𝜇superscript𝑇𝐴subscript𝑞𝐿subscript¯𝑡𝑅superscript𝛾𝜇superscript𝑇𝐴subscript𝑡𝑅\mathcal{O}_{\scriptscriptstyle tq}^{\scriptscriptstyle(8)}=(\bar{q}_{L}\gamma% _{\mu}T^{A}q_{L})(\bar{t}_{R}\gamma^{\mu}T^{A}t_{R}),caligraphic_O start_POSTSUBSCRIPT italic_t italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT = ( over¯ start_ARG italic_q end_ARG start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_q start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) ,

and the coupling between right-handed tops and right-handed quarks,

𝒪tt(8)=(t¯RγμTAtR)(t¯RγμTAtR),𝒪t(u/d)(8)=(t¯RγμTAtR)((u¯R/d¯R)γμTA(uR/dR)).formulae-sequencesuperscriptsubscript𝒪𝑡𝑡8subscript¯𝑡𝑅subscript𝛾𝜇superscript𝑇𝐴subscript𝑡𝑅subscript¯𝑡𝑅superscript𝛾𝜇superscript𝑇𝐴subscript𝑡𝑅superscriptsubscript𝒪𝑡𝑢𝑑8subscript¯𝑡𝑅subscript𝛾𝜇superscript𝑇𝐴subscript𝑡𝑅subscript¯𝑢𝑅subscript¯𝑑𝑅superscript𝛾𝜇superscript𝑇𝐴subscript𝑢𝑅subscript𝑑𝑅\mathcal{O}_{\scriptscriptstyle tt}^{\scriptscriptstyle(8)}=(\bar{t}_{R}\gamma% _{\mu}T^{A}t_{R})(\bar{t}_{R}\gamma^{\mu}T^{A}t_{R}),\,\mathcal{O}_{% \scriptscriptstyle t(u/d)}^{\scriptscriptstyle(8)}=(\bar{t}_{R}\gamma_{\mu}T^{% A}t_{R})((\bar{u}_{R}/\bar{d}_{R})\gamma^{\mu}T^{A}(u_{R}/d_{R})).caligraphic_O start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT = ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) , caligraphic_O start_POSTSUBSCRIPT italic_t ( italic_u / italic_d ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT = ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) ( ( over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT / over¯ start_ARG italic_d end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT ( italic_u start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT / italic_d start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) ) .

In Ref. Ellis:2020unq , we showed that the operator 𝒪tGsubscript𝒪𝑡𝐺\mathcal{O}_{\scriptscriptstyle tG}caligraphic_O start_POSTSUBSCRIPT italic_t italic_G end_POSTSUBSCRIPT was mostly constrained by Run 1 and Run 2 Higgs observables plus the TeVatron and LHC tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG datasets. Moreover, among the four-fermion operators, only the 𝒪tX(8)superscriptsubscript𝒪𝑡𝑋8\mathcal{O}_{\scriptscriptstyle tX}^{\scriptscriptstyle(8)}caligraphic_O start_POSTSUBSCRIPT italic_t italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT (X𝑋Xitalic_X= q𝑞qitalic_q, u𝑢uitalic_u or d𝑑ditalic_d) operators were constrained, predominantly by the top data and, to a lesser extent, by tt¯V𝑡¯𝑡𝑉t\bar{t}Vitalic_t over¯ start_ARG italic_t end_ARG italic_V measurements. The current limits on these four operators are shown in Table 1, and one can see that the difference between the individual and marginalised limits is quite dramatic, particularly for the four-fermion operators. The reason is that marginalised limits correspond to a global fit to many operators, beyond these shown here, which contribute to the same set of observables. Hence, the inclusion of more operators tend to weaken the limits for each operator.

Operator Individual fit (TeV-2) Marginalised fit (TeV-2)
𝒪tGsubscript𝒪𝑡𝐺\mathcal{O}_{\scriptscriptstyle tG}caligraphic_O start_POSTSUBSCRIPT italic_t italic_G end_POSTSUBSCRIPT 0.010.1+0.086subscriptsuperscript0.010.0860.1-0.01^{+0.086}_{-0.1}- 0.01 start_POSTSUPERSCRIPT + 0.086 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.1 end_POSTSUBSCRIPT 0.360.6+0.12subscriptsuperscriptabsent0.120.6{}^{+0.12}_{-0.6}start_FLOATSUPERSCRIPT + 0.12 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 0.6 end_POSTSUBSCRIPT
𝒪tq(8)superscriptsubscript𝒪𝑡𝑞8\mathcal{O}_{\scriptscriptstyle tq}^{\scriptscriptstyle(8)}caligraphic_O start_POSTSUBSCRIPT italic_t italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT -0.40.85+0.06subscriptsuperscriptabsent0.060.85{}^{+0.06}_{-0.85}start_FLOATSUPERSCRIPT + 0.06 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 0.85 end_POSTSUBSCRIPT 5.13+2.2subscriptsuperscriptabsent2.213{}^{+2.2}_{-13}start_FLOATSUPERSCRIPT + 2.2 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 13 end_POSTSUBSCRIPT
𝒪tu(8)superscriptsubscript𝒪𝑡𝑢8\mathcal{O}_{\scriptscriptstyle tu}^{\scriptscriptstyle(8)}caligraphic_O start_POSTSUBSCRIPT italic_t italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT -0.451.1+0.23subscriptsuperscriptabsent0.231.1{}^{+0.23}_{-1.1}start_FLOATSUPERSCRIPT + 0.23 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 1.1 end_POSTSUBSCRIPT 4.0 11+19subscriptsuperscriptabsent1911{}^{+19}_{-11}start_FLOATSUPERSCRIPT + 19 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 11 end_POSTSUBSCRIPT
𝒪td(8)superscriptsubscript𝒪𝑡𝑑8\mathcal{O}_{\scriptscriptstyle td}^{\scriptscriptstyle(8)}caligraphic_O start_POSTSUBSCRIPT italic_t italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT -1.02.5+0.38subscriptsuperscriptabsent0.382.5{}^{+0.38}_{-2.5}start_FLOATSUPERSCRIPT + 0.38 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 2.5 end_POSTSUBSCRIPT -0.4212+11subscriptsuperscriptabsent1112{}^{+11}_{-12}start_FLOATSUPERSCRIPT + 11 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 12 end_POSTSUBSCRIPT
Table 1: Limits on the relevant subset of Top operators at 95% C.L., from Ref. Ellis:2020unq .

3 Beyond Top EFT

The Top EFT operators discussed in the previous section are useful for describing new physics effects on energy scales well below the BSM masses. However, once we consider the TeV energies probed by the LHC and BSM particles with masses around 1 TeV, the validity of the EFT regime is not guaranteed. In this case we need to go beyond the Top EFT and consider the UV extension of the SM. Motivated by Dark Matter, we consider a BSM scenario with a 𝒵2subscript𝒵2\mathcal{Z}_{2}caligraphic_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT parity, which ensures the stability of the DM candidate. An important consequence of this assumption is that it forbids linear couplings of the new states to two SM particles and the SMEFT operators are only induced at the one-loop level Cepedello:2022pyx . In addition we will build in this model a special connection to the top sector, a possibility that has been partly explored in the context of DM relic abundance and collider phenomenology, see Refs. Haisch:2015ioa ; Plehn:2017bys ; Delgado:2016umt ; Garny:2018icg .

3.1 An explicit example: a UV extension with Dark Matter and a top partner

In order to incorporate the main features described above and be minimal, we consider the simple case of a scalar top partner (φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT), singlet under SU(2)L𝑆𝑈subscript2𝐿SU(2)_{L}italic_S italic_U ( 2 ) start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT, and a singlet fermion (χ𝜒\chiitalic_χ), which is a Dark Matter candidate. Under the imposed 𝒵2subscript𝒵2\mathcal{Z}_{2}caligraphic_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT symmetry the BSM fields are odd and the SM are even, so the renormalizable BSM lagrangian becomes:

BSM=χ¯(i∂̸12mχ)χ+|DμφT|2mT2|φT|2(yDMφTχ¯tR+h.c.)\mathcal{L}_{BSM}=\bar{\chi}\left(i\not{\partial}-\frac{1}{2}m_{\chi}\right)% \chi+|D_{\mu}\varphi_{T}|^{2}-m_{T}^{2}|\varphi_{T}|^{2}-\left(y_{\mathrm{DM}}% \varphi_{T}^{\dagger}\bar{\chi}t_{R}+h.c.\right)caligraphic_L start_POSTSUBSCRIPT italic_B italic_S italic_M end_POSTSUBSCRIPT = over¯ start_ARG italic_χ end_ARG ( italic_i ∂̸ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT ) italic_χ + | italic_D start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_y start_POSTSUBSCRIPT roman_DM end_POSTSUBSCRIPT italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over¯ start_ARG italic_χ end_ARG italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT + italic_h . italic_c . ) (1)

with mT>mχsubscript𝑚𝑇subscript𝑚𝜒m_{T}>m_{\chi}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT > italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT, so the DM candidate is stable. The only viable decay channel for φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT is φTχ+tsubscript𝜑𝑇𝜒𝑡\varphi_{T}\to\chi+titalic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT → italic_χ + italic_t, where the top is off-shell if ΔM=mTmχ<mtΔ𝑀subscript𝑚𝑇subscript𝑚𝜒subscript𝑚𝑡\Delta M=m_{T}-m_{\chi}<m_{t}roman_Δ italic_M = italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT - italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT < italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT. Note that the interactions between the Dark Matter candidate and the SM are fully controlled by the yDMsubscript𝑦DMy_{\mathrm{DM}}italic_y start_POSTSUBSCRIPT roman_DM end_POSTSUBSCRIPT coupling. Although this scenario can be phenomenologically similar to the minimal supersymmetric standard model (MSSM) with a Bino LSP and a right stop, we point out that in the supersymmetric case the yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT coupling is fixed by the LSP composition and it is of the order of the EW couplings (yDM0.11similar-tosubscript𝑦𝐷𝑀0.11y_{DM}\sim 0.1-1italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ∼ 0.1 - 1). In the scenario discussed here we assume yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT to be a free parameter, which can be as large as allowed by perturbativity, yDM4πless-than-or-similar-tosubscript𝑦𝐷𝑀4𝜋y_{DM}\lesssim 4\piitalic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≲ 4 italic_π. The Dark Matter implications of this scenario were studied in Ref. Garny:2018icg , where it has been shown that the correct Dark Matter relic density can be achieved for a wide range of mass values: 10 GeVmχ<50 TeVless-than-or-similar-to10 GeVsubscript𝑚𝜒50 TeV10\mbox{ GeV}\lesssim m_{\chi}<50\mbox{ TeV}10 GeV ≲ italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT < 50 TeV and ΔM500 GeVless-than-or-similar-toΔ𝑀500 GeV\Delta M\lesssim 500\mbox{ GeV}roman_Δ italic_M ≲ 500 GeV as long as the value of yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT is properly chosen. In particular, large coupling values, e.g. yDM3greater-than-or-equivalent-tosubscript𝑦𝐷𝑀3y_{DM}\gtrsim 3italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≳ 3, are needed in some regions of parameter space. In this work we do not impose any Dark Matter constraints, since these can be modified by the presence of additional (heavy) BSM states and/or a non-standard cosmological evolution.

The above model can lead to several implications at the LHC and low energy observables. In this work we are mostly interested in the complementarity between direct searches for the top scalar and constraints from top pair production observables. A full study of the direct and indirect constraints on the BSM model is left for a future work.

3.1.1 EFT Limit

In the heavy mass limit (mT,mχmt,sformulae-sequencemuch-greater-thansubscript𝑚𝑇subscript𝑚𝜒subscript𝑚𝑡𝑠m_{T},m_{\chi}\gg m_{t},\sqrt{s}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT ≫ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , square-root start_ARG italic_s end_ARG) the BSM contributions for the model defined in Sec. 3.1 can be described by an effective field theory, where the colored scalar and dark fermion have been integrated out. In this case we have the following dimension-six effective Lagrangian:

EFT=mtCgGμνA(t¯TAσμνt)subscript𝐸𝐹𝑇subscript𝑚𝑡subscript𝐶𝑔superscriptsubscript𝐺𝜇𝜈𝐴¯𝑡superscript𝑇𝐴superscript𝜎𝜇𝜈𝑡\displaystyle\mathcal{L}_{EFT}=m_{t}C_{g}\;G_{\mu\nu}^{A}\left(\bar{t}T^{A}% \sigma^{\mu\nu}t\right)caligraphic_L start_POSTSUBSCRIPT italic_E italic_F italic_T end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_G start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT ( over¯ start_ARG italic_t end_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_σ start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT italic_t ) +\displaystyle++ Cq(t¯RTAγμtR)(Q¯LTAγμQL+u¯RTAγμuR+d¯RTAγμdR)subscript𝐶𝑞subscript¯𝑡𝑅superscript𝑇𝐴superscript𝛾𝜇subscript𝑡𝑅subscript¯𝑄𝐿superscript𝑇𝐴superscript𝛾𝜇subscript𝑄𝐿subscript¯𝑢𝑅superscript𝑇𝐴superscript𝛾𝜇subscript𝑢𝑅subscript¯𝑑𝑅superscript𝑇𝐴superscript𝛾𝜇subscript𝑑𝑅\displaystyle C_{q}\left(\bar{t}_{R}T^{A}\gamma^{\mu}t_{R}\right)\left(\bar{Q}% _{L}T^{A}\gamma^{\mu}Q_{L}+\bar{u}_{R}T^{A}\gamma^{\mu}u_{R}+\bar{d}_{R}T^{A}% \gamma^{\mu}d_{R}\right)italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) ( over¯ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT + over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT + over¯ start_ARG italic_d end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_d start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) (2)
+\displaystyle++ Cq(t¯RTAγμtR)(Q¯3,LTAγμQ3,L)subscript𝐶𝑞subscript¯𝑡𝑅superscript𝑇𝐴superscript𝛾𝜇subscript𝑡𝑅subscript¯𝑄3𝐿superscript𝑇𝐴superscript𝛾𝜇subscript𝑄3𝐿\displaystyle C_{q}\left(\bar{t}_{R}T^{A}\gamma^{\mu}t_{R}\right)\left(\bar{Q}% _{3,L}T^{A}\gamma^{\mu}Q_{3,L}\right)italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) ( over¯ start_ARG italic_Q end_ARG start_POSTSUBSCRIPT 3 , italic_L end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 3 , italic_L end_POSTSUBSCRIPT )
+\displaystyle++ CtR(t¯RTAγμtR)(t¯RTAγμtR)subscript𝐶𝑡𝑅subscript¯𝑡𝑅superscript𝑇𝐴superscript𝛾𝜇subscript𝑡𝑅subscript¯𝑡𝑅superscript𝑇𝐴superscript𝛾𝜇subscript𝑡𝑅\displaystyle C_{tR}\left(\bar{t}_{R}T^{A}\gamma^{\mu}t_{R}\right)\left(\bar{t% }_{R}T^{A}\gamma^{\mu}t_{R}\right)italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) ( over¯ start_ARG italic_t end_ARG start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT )

where u,d,Q𝑢𝑑𝑄u,d,Qitalic_u , italic_d , italic_Q represent any light quark flavor, Q3,Lsubscript𝑄3𝐿Q_{3,L}italic_Q start_POSTSUBSCRIPT 3 , italic_L end_POSTSUBSCRIPT represents the 3rd generation quark doublet and mtsubscript𝑚𝑡m_{t}italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is the (on-shell) top mass.

The connection with the SMEFT operators in the Warsaw basis described in Sec. 2 and the two operators Cg,qsubscript𝐶𝑔𝑞C_{g,q}italic_C start_POSTSUBSCRIPT italic_g , italic_q end_POSTSUBSCRIPT is as follows

Cgsubscript𝐶𝑔\displaystyle C_{g}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT =\displaystyle== yt1CtGΛ2,superscriptsubscript𝑦𝑡1subscript𝐶𝑡𝐺superscriptΛ2\displaystyle y_{t}^{-1}\frac{C_{tG}}{\Lambda^{2}},italic_y start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT divide start_ARG italic_C start_POSTSUBSCRIPT italic_t italic_G end_POSTSUBSCRIPT end_ARG start_ARG roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,
Cqsubscript𝐶𝑞\displaystyle C_{q}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT =\displaystyle== Ctq(8)Λ2=Ct(u/d)(8)Λ2superscriptsubscript𝐶𝑡𝑞8superscriptΛ2superscriptsubscript𝐶𝑡𝑢𝑑8superscriptΛ2\displaystyle\frac{C_{tq}^{(8)}}{\Lambda^{2}}=\frac{C_{t(u/d)}^{(8)}}{\Lambda^% {2}}divide start_ARG italic_C start_POSTSUBSCRIPT italic_t italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT end_ARG start_ARG roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = divide start_ARG italic_C start_POSTSUBSCRIPT italic_t ( italic_u / italic_d ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT end_ARG start_ARG roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG
CtRsubscript𝐶𝑡𝑅\displaystyle C_{tR}italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT =\displaystyle== Ctt(8)Λ2,superscriptsubscript𝐶𝑡𝑡8superscriptΛ2\displaystyle\frac{C_{tt}^{(8)}}{\Lambda^{2}},divide start_ARG italic_C start_POSTSUBSCRIPT italic_t italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT end_ARG start_ARG roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,

where ytsubscript𝑦𝑡y_{t}italic_y start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is the top Yukawa. Unlike the general SMEFT framework, the Cgsubscript𝐶𝑔C_{g}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT and Cqsubscript𝐶𝑞C_{q}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT coefficients are correlated and determined by the underlying UV parameters. These coefficients were computed using Matchete Fuentes-Martin:2022jrf and are given by:

Cgsubscript𝐶𝑔\displaystyle C_{g}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT =\displaystyle== gsyDM2384π21mT21(1x)4[16x+3x2+2x36x2log(x)]subscript𝑔𝑠superscriptsubscript𝑦𝐷𝑀2384superscript𝜋21superscriptsubscript𝑚𝑇21superscript1𝑥4delimited-[]16𝑥3superscript𝑥22superscript𝑥36superscript𝑥2𝑥\displaystyle-\frac{g_{s}y_{DM}^{2}}{384\pi^{2}}\frac{1}{m_{T}^{2}}\frac{1}{% \left(1-x\right)^{4}}\left[1-6x+3x^{2}+2x^{3}-6x^{2}\log(x)\right]- divide start_ARG italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 384 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG ( 1 - italic_x ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG [ 1 - 6 italic_x + 3 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 6 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_log ( italic_x ) ] (3)
Cqsubscript𝐶𝑞\displaystyle C_{q}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT =\displaystyle== gs2yDM2576π21mT21(1x)4[29x+18x211x3+6x3log(x)]superscriptsubscript𝑔𝑠2superscriptsubscript𝑦𝐷𝑀2576superscript𝜋21superscriptsubscript𝑚𝑇21superscript1𝑥4delimited-[]29𝑥18superscript𝑥211superscript𝑥36superscript𝑥3𝑥\displaystyle\frac{g_{s}^{2}y_{DM}^{2}}{576\pi^{2}}\frac{1}{m_{T}^{2}}\frac{1}% {\left(1-x\right)^{4}}\left[2-9x+18x^{2}-11x^{3}+6x^{3}\log(x)\right]divide start_ARG italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 576 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG ( 1 - italic_x ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG [ 2 - 9 italic_x + 18 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 11 italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 6 italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_log ( italic_x ) ] (4)
CtRsubscript𝐶𝑡𝑅\displaystyle C_{tR}italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT =\displaystyle== yDM4128π21mT21(1x)3[1x2+2xlog(x)]superscriptsubscript𝑦𝐷𝑀4128superscript𝜋21superscriptsubscript𝑚𝑇21superscript1𝑥3delimited-[]1superscript𝑥22𝑥𝑥\displaystyle-\frac{y_{DM}^{4}}{128\pi^{2}}\frac{1}{m_{T}^{2}}\frac{1}{\left(1% -x\right)^{3}}\left[1-x^{2}+2x\log(x)\right]- divide start_ARG italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 128 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG ( 1 - italic_x ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG [ 1 - italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_x roman_log ( italic_x ) ] (5)

where mtsubscript𝑚𝑡m_{t}italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is the top mass and xmχ2/mT2𝑥superscriptsubscript𝑚𝜒2superscriptsubscript𝑚𝑇2x\equiv m_{\chi}^{2}/m_{T}^{2}italic_x ≡ italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. For example, taking the limit x1𝑥1x\to 1italic_x → 1 (mTmχsimilar-to-or-equalssubscript𝑚𝑇subscript𝑚𝜒m_{T}\simeq m_{\chi}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≃ italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT), one finds that this model produces a particular pattern in the SMEFT parameter space:

Cgsimilar-to-or-equalssubscript𝐶𝑔absent\displaystyle C_{g}\simeqitalic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ≃ 12gsyDM2384π21mT2 , Cq32gs2yDM2576π21mT2 , CtR13yDM4128π21mT2similar-to-or-equals12subscript𝑔𝑠superscriptsubscript𝑦𝐷𝑀2384superscript𝜋21superscriptsubscript𝑚𝑇2 , subscript𝐶𝑞32superscriptsubscript𝑔𝑠2superscriptsubscript𝑦𝐷𝑀2576superscript𝜋21superscriptsubscript𝑚𝑇2 , subscript𝐶𝑡𝑅similar-to-or-equals13superscriptsubscript𝑦𝐷𝑀4128superscript𝜋21superscriptsubscript𝑚𝑇2\displaystyle-\frac{1}{2}\frac{g_{s}y_{DM}^{2}}{384\pi^{2}}\frac{1}{m_{T}^{2}}% \mbox{ , }C_{q}\simeq\frac{3}{2}\frac{g_{s}^{2}y_{DM}^{2}}{576\pi^{2}}\frac{1}% {m_{T}^{2}}\mbox{ , }C_{tR}\simeq-\frac{1}{3}\frac{y_{DM}^{4}}{128\pi^{2}}% \frac{1}{m_{T}^{2}}- divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 384 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ≃ divide start_ARG 3 end_ARG start_ARG 2 end_ARG divide start_ARG italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 576 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT ≃ - divide start_ARG 1 end_ARG start_ARG 3 end_ARG divide start_ARG italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 128 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (6)
Cq=absentsubscript𝐶𝑞absent\displaystyle\Rightarrow C_{q}=⇒ italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = 2gsCg2subscript𝑔𝑠subscript𝐶𝑔\displaystyle-2\,g_{s}\,C_{g}- 2 italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT (7)

We point out that for the results discussed later we only consider BSM contributions up to order yDM2superscriptsubscript𝑦𝐷𝑀2y_{DM}^{2}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and the CtRsubscript𝐶𝑡𝑅C_{tR}italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT operator will be ignored. Furthermore, the 4-top operator CtRsubscript𝐶𝑡𝑅C_{tR}italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT contribution to the top pair production at the LHC is negligible, as it corresponds to a two-loop contribution, and its inclusion would not change the results presented in this work. Nonetheless, this four-top operator could be searched for in the four-top final state, which has recently been observed at the LHC CMS:2023ftu ; ATLAS:2023ajo . Despite the strong yDM4superscriptsubscript𝑦𝐷𝑀4y_{DM}^{4}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT scaling, the bounds on our benchmark would be weaker than those from tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG final state, due to the kinematic suppression due to a top radiating three tops from CtRsubscript𝐶𝑡𝑅C_{tR}italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT and the inherent loop suppression. Note, also, that a large four-top operator is constrained by perturbativity and, in particular, by the absence of tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG bound states. As it was first discussed in Ref. Bardeen:1992mhn , bound states would be formed when four-fermion interactions like CtRsubscript𝐶𝑡𝑅C_{tR}italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT grow above some critical value, estimated to be CtR8π23Λ2similar-to-or-equalssubscript𝐶𝑡𝑅8superscript𝜋23superscriptΛ2C_{tR}\simeq\frac{8\pi^{2}}{3\Lambda^{2}}italic_C start_POSTSUBSCRIPT italic_t italic_R end_POSTSUBSCRIPT ≃ divide start_ARG 8 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG, where ΛΛ\Lambdaroman_Λ represents the scale at which confinement occurs. Comparing with the matching in Eq. 6 and approximating ΛmTsimilar-to-or-equalsΛsubscript𝑚𝑇\Lambda\simeq m_{T}roman_Λ ≃ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, non perturbativity would require yDM18greater-than-or-equivalent-tosubscript𝑦𝐷𝑀18y_{DM}\gtrsim 18italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≳ 18. Therefore, for the results presented below, we impose yDM10subscript𝑦𝐷𝑀10y_{DM}\leq 10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≤ 10.

The size of the Cqsubscript𝐶𝑞C_{q}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT and Cgsubscript𝐶𝑔C_{g}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT coefficients and their ratio is illustrated in Fig. 1. As we can see, for BSM masses around 1 TeV and yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5, the coefficients are Cg,q102similar-tosubscript𝐶𝑔𝑞superscript102C_{g,q}\sim 10^{-2}italic_C start_POSTSUBSCRIPT italic_g , italic_q end_POSTSUBSCRIPT ∼ 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT TeV-2. Furthermore, we see that Cqsubscript𝐶𝑞C_{q}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT is typically 22.5similar-toabsent22.5\sim 2-2.5∼ 2 - 2.5 times larger than |Cg|subscript𝐶𝑔|C_{g}|| italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT | and Cgsubscript𝐶𝑔C_{g}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT is always negative. As discussed in Sec. 2, usual SMEFT analysis constrain these coefficients to 0.110.110.1-10.1 - 1 TeV-2. Although these constraints can not be directly applied to our scenario, since both coefficients are present and correlated, one would still expect that, for yDM5less-than-or-similar-tosubscript𝑦𝐷𝑀5y_{DM}\lesssim 5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≲ 5, only the sub-TeV region of parameter space can be tested.

Refer to caption
Refer to caption
Refer to caption
Figure 1: Values for the EFT coefficients Cgsubscript𝐶𝑔C_{g}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT and Cqsubscript𝐶𝑞C_{q}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT computed according to Eqs.( 3) and ( 4). The BSM and strong couplings were taken as yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5 and αs=0.13subscript𝛼𝑠0.13\alpha_{s}=0.13italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.13. The shaded region corresponds to mT<mχsubscript𝑚𝑇subscript𝑚𝜒m_{T}<m_{\chi}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT < italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT, which corresponds to a stable colored scalar. The bottom plot shows the values taken by the ratio of both coefficients.

3.1.2 1-Loop Form Factors

The EFT approach discussed in Sec. 3.1.1 is only valid for energies well below the BSM masses. As we will show in Sec. 4, at the LHC it is possible to probe distributions at energies up to a few TeV. Therefore the EFT validity is not guaranteed when using such measurements to look for new physics. In this case we need to compute the full loop contributions to the relevant observables, which are valid at any scale. For sufficiently high values of BSM masses, the loop contributions should reproduce the EFT results.

In order to compute the 1-loop contributions to top pair production distributions, we have computed form factors for the effective top-top-gluon and top-top-gluon-gluon couplings induced by the loop diagrams shown in Fig. 2. The form factors can be then written as effective, momentum dependent couplings:

FF=π2gsyDM2Gμt¯[μ(pt,pt¯)]t+π2gs2yDM2GμGνt¯[μν(pg,pt,pt¯)]tsubscript𝐹𝐹superscript𝜋2subscript𝑔𝑠superscriptsubscript𝑦𝐷𝑀2subscript𝐺𝜇¯𝑡delimited-[]superscript𝜇subscript𝑝𝑡subscript𝑝¯𝑡𝑡superscript𝜋2superscriptsubscript𝑔𝑠2superscriptsubscript𝑦𝐷𝑀2subscript𝐺𝜇subscript𝐺𝜈¯𝑡delimited-[]superscript𝜇𝜈subscript𝑝𝑔subscript𝑝𝑡subscript𝑝¯𝑡𝑡\mathcal{L}_{FF}=\pi^{2}g_{s}y_{DM}^{2}G_{\mu}\bar{t}\left[\mathcal{F}^{\mu}% \left(p_{t},p_{\bar{t}}\right)\right]t+\pi^{2}g_{s}^{2}y_{DM}^{2}G_{\mu}G_{\nu% }\bar{t}\left[\mathcal{F}^{\mu\nu}\left(p_{g},p_{t},p_{\bar{t}}\right)\right]tcaligraphic_L start_POSTSUBSCRIPT italic_F italic_F end_POSTSUBSCRIPT = italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_G start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG [ caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ) ] italic_t + italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_G start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_G start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG [ caligraphic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT ( italic_p start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ) ] italic_t (8)

where the μsuperscript𝜇\mathcal{F}^{\mu}caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT and μνsuperscript𝜇𝜈\mathcal{F}^{\mu\nu}caligraphic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT form factors contain the full momenta dependence as well as the Dirac and color structures, which we suppress for simplicity. In order to determine these functions, all the diagrams shown in Fig. 2 were computed using FeynArts Hahn:2000kx and FeynCalc Shtabovenko:2020gxv ; Shtabovenko:2016sxi ; Mertig:1990an and the results were then used to extract μsuperscript𝜇\mathcal{F}^{\mu}caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT and μνsuperscript𝜇𝜈\mathcal{F}^{\mu\nu}caligraphic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT. The form factors also include the counter-terms required for renormalizing the top self-energy and the top-top-gluon vertex, which were computed using NLOCT Degrande:2014vpa under the on-shell renormalization scheme.

Refer to caption
Figure 2: Loop diagrams used to compute the form factors defined in Eq.( 8).

While μνsuperscript𝜇𝜈\mathcal{F}^{\mu\nu}caligraphic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT contains quite a large number of terms and an involved tensor structure, the expression for the gtt¯𝑔𝑡¯𝑡g-t-\bar{t}italic_g - italic_t - over¯ start_ARG italic_t end_ARG form factor can be written in a compact form using the triangular Passarino-Veltmann loop functions Passarino:1978jh ; tHooft:1978jhc :

μ=iTAsuperscript𝜇𝑖superscript𝑇𝐴\displaystyle\mathcal{F}^{\mu}=iT^{A}caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = italic_i italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT {1[p1μ(C1+2C11)p2μ(C1+2C12)]\displaystyle\biggl{\{}\not{p}_{1}\left[p_{1}^{\mu}\left(C_{1}+2C_{11}\right)-% p_{2}^{\mu}\left(C_{1}+2C_{12}\right)\right]{ italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ) - italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) ]
+2[p2μ(C2+2C22)p1μ(C2+2C12)]subscriptitalic-p̸2delimited-[]superscriptsubscript𝑝2𝜇subscript𝐶22subscript𝐶22superscriptsubscript𝑝1𝜇subscript𝐶22subscript𝐶12\displaystyle+\not{p}_{2}\left[p_{2}^{\mu}\left(C_{2}+2C_{22}\right)-p_{1}^{% \mu}\left(C_{2}+2C_{12}\right)\right]+ italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT [ italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 22 end_POSTSUBSCRIPT ) - italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) ]
+2γμ(C00+δR)}PR+2iTAγμδLPL\displaystyle+2\gamma^{\mu}\left(C_{00}+\delta_{R}\right)\vphantom{\left[\left% (p_{1}\right)\right]}\biggr{\}}P_{R}+2iT^{A}\gamma^{\mu}\delta_{L}P_{L}+ 2 italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) } italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT + 2 italic_i italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT (9)

where PR,L=12(1±γ5)subscript𝑃𝑅𝐿12plus-or-minus1superscript𝛾5P_{R,L}=\frac{1}{2}\left(1\pm\gamma^{5}\right)italic_P start_POSTSUBSCRIPT italic_R , italic_L end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( 1 ± italic_γ start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT ), p1,p2subscript𝑝1subscript𝑝2p_{1},p_{2}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are the top and anti-top momenta and s=(p1+p2)2𝑠superscriptsubscript𝑝1subscript𝑝22s=(p_{1}+p_{2})^{2}italic_s = ( italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. All the loop integrals Ci,Cijsubscript𝐶𝑖subscript𝐶𝑖𝑗C_{i},C_{ij}italic_C start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_C start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT are functions of (p12,s,p22)superscriptsubscript𝑝12𝑠superscriptsubscript𝑝22(p_{1}^{2},s,p_{2}^{2})( italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_s , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) and δL,Rsubscript𝛿𝐿𝑅\delta_{L,R}italic_δ start_POSTSUBSCRIPT italic_L , italic_R end_POSTSUBSCRIPT are the counter-terms obtained using the on-shell renormalization scheme (see Appendix B for more details). Although the analytical expressions for the loop integrals are quite involved in the general case, in the mass degenerate limit (mχmTsimilar-to-or-equalssubscript𝑚𝜒subscript𝑚𝑇m_{\chi}\simeq m_{T}italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT ≃ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT) and neglecting the top mass (mtmT,mχ,smuch-less-thansubscript𝑚𝑡subscript𝑚𝑇subscript𝑚𝜒𝑠m_{t}\ll m_{T},m_{\chi},\sqrt{s}italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≪ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT , square-root start_ARG italic_s end_ARG) they simplify considerably, resulting in:

μ=superscript𝜇absent\displaystyle\mathcal{F}^{\mu}=caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = iTA32π4[2(2p1μ+1p2μ)sγμ]PR(s,mT2) (mχmTmt0),𝑖superscript𝑇𝐴32superscript𝜋4delimited-[]2subscriptitalic-p̸2superscriptsubscript𝑝1𝜇subscriptitalic-p̸1superscriptsubscript𝑝2𝜇𝑠superscript𝛾𝜇subscript𝑃𝑅𝑠superscriptsubscript𝑚𝑇2 (mχmTmt0)\displaystyle i\frac{T^{A}}{32\pi^{4}}\left[2\left(\not{p}_{2}p_{1}^{\mu}+\not% {p}_{1}p_{2}^{\mu}\right)-s\gamma^{\mu}\right]P_{R}\mathcal{F}(s,m_{T}^{2})\;% \mbox{ ($m_{\chi}\to m_{T}$, $m_{t}\to 0$)},italic_i divide start_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT end_ARG start_ARG 32 italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG [ 2 ( italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT + italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ) - italic_s italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ] italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT caligraphic_F ( italic_s , italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ( italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT → italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT → 0 ) , (10)

where:

(s,mT2)𝑠superscriptsubscript𝑚𝑇2absent\displaystyle\mathcal{F}(s,m_{T}^{2})\equivcaligraphic_F ( italic_s , italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ≡ 1s2[mT2L2(s,mT2)+s(s4mT2)L(s,mT2)+3s] and1superscript𝑠2delimited-[]superscriptsubscript𝑚𝑇2superscript𝐿2𝑠superscriptsubscript𝑚𝑇2𝑠𝑠4superscriptsubscript𝑚𝑇2𝐿𝑠superscriptsubscript𝑚𝑇23𝑠 and\displaystyle\frac{1}{s^{2}}\left[m_{T}^{2}L^{2}(s,m_{T}^{2})+\sqrt{s\left(s-4% m_{T}^{2}\right)}L(s,m_{T}^{2})+3s\right]\mbox{ and }divide start_ARG 1 end_ARG start_ARG italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_s , italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + square-root start_ARG italic_s ( italic_s - 4 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG italic_L ( italic_s , italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) + 3 italic_s ] and (11)
L(s,mT2)𝐿𝑠superscriptsubscript𝑚𝑇2absent\displaystyle L(s,m_{T}^{2})\equivitalic_L ( italic_s , italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ≡ log[(s(s4mT2)+2mT2s)/(2mT2)]𝑠𝑠4superscriptsubscript𝑚𝑇22superscriptsubscript𝑚𝑇2𝑠2superscriptsubscript𝑚𝑇2\displaystyle\log\left[\left(\sqrt{s\left(s-4m_{T}^{2}\right)}+2m_{T}^{2}-s% \right)/(2m_{T}^{2})\right]roman_log [ ( square-root start_ARG italic_s ( italic_s - 4 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG + 2 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_s ) / ( 2 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ]
Refer to caption
Figure 3: Dependence of form factor on the top pair invariant mass s=m(tt¯)2𝑠𝑚superscript𝑡¯𝑡2s=m(t\bar{t})^{2}italic_s = italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for mχ=mTsubscript𝑚𝜒subscript𝑚𝑇m_{\chi}=m_{T}italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT and mt0subscript𝑚𝑡0m_{t}\to 0italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT → 0. The curves correspond to the real, imaginary and absolute values of the factor \mathcal{F}caligraphic_F defined in Eq.(11).

In Fig. 3 we show the form factor dependence on s𝑠sitalic_s coming from the (s,mT2)𝑠superscriptsubscript𝑚𝑇2\mathcal{F}(s,m_{T}^{2})caligraphic_F ( italic_s , italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) factor defined in Eq.( 11). As we can see, \mathcal{F}caligraphic_F peaks around s5mT2similar-to-or-equals𝑠5superscriptsubscript𝑚𝑇2s\simeq 5m_{T}^{2}italic_s ≃ 5 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT or m(tt¯)5mTsimilar-to-or-equals𝑚𝑡¯𝑡5subscript𝑚𝑇m(t\bar{t})\simeq\sqrt{5}m_{T}italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ≃ square-root start_ARG 5 end_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, displaying a very broad resonant behavior. We also see that for smT2much-greater-than𝑠superscriptsubscript𝑚𝑇2s\gg m_{T}^{2}italic_s ≫ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT the form factor is dominated by its imaginary part, resulting in a negative interference term, as discussed below.

3.2 Differential distributions in tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG

Since the BSM model discussed here mostly couples to the top quark, it can impact top pair production and its measured distributions. In order to illustrate these effects, we consider the top pair invariant mass m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) and the top transverse momentum pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ). We have implemented the lagrangians defined in Eqs.(2) and (8) in the UFO Degrande:2011ua format, which allows us to generate events using either the full 1-loop calculation or the EFT approximation. For all the results discussed below we have produced 150k MC events for pptt¯𝑝𝑝𝑡¯𝑡pp\to t\bar{t}italic_p italic_p → italic_t over¯ start_ARG italic_t end_ARG at parton level using MadGraph5_aMC@NLO Alwall:2014hca ; Frederix:2018nkq with an invariant mass bias, so we can appropriately describe the high energy tail of the distributions. The top quarks were then decayed using MadSpin MadSpin . We have used the PDF set NNPDF23_nlo_as_0119 and the factorization and normalization scales were set to the top transverse mass: μF=μR=mt2+pT2subscript𝜇𝐹subscript𝜇𝑅superscriptsubscript𝑚𝑡2superscriptsubscript𝑝𝑇2\mu_{F}=\mu_{R}=\sqrt{m_{t}^{2}+p_{T}^{2}}italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = square-root start_ARG italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG.

The distributions were computed at leading order in αssubscript𝛼𝑠\alpha_{s}italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT (𝒪(αs,αsyDM2)𝒪subscript𝛼𝑠subscript𝛼𝑠superscriptsubscript𝑦𝐷𝑀2\mathcal{O}(\alpha_{s},\alpha_{s}y_{DM}^{2})caligraphic_O ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )), which corresponds to the Born plus the interference terms:

||2=|SM|2+2Re(SMBSM)superscript2superscriptsubscriptSM22ResubscriptsuperscriptSMsubscriptBSM|\mathcal{M}|^{2}=|\mathcal{M}_{\rm SM}|^{2}+2\mathrm{Re}\left({\mathcal{M}^{*% }_{\rm SM}\mathcal{M}_{\rm BSM}}\right)| caligraphic_M | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = | caligraphic_M start_POSTSUBSCRIPT roman_SM end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 roman_R roman_e ( caligraphic_M start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_SM end_POSTSUBSCRIPT caligraphic_M start_POSTSUBSCRIPT roman_BSM end_POSTSUBSCRIPT ) (12)

In the EFT approach this is equivalent to keeping only the 1/Λ21superscriptΛ21/\Lambda^{2}1 / roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT terms.222We have verified that the contribution from the 1/Λ41superscriptΛ41/\Lambda^{4}1 / roman_Λ start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT term (|BSM|2superscriptsubscriptBSM2|\mathcal{M}_{\rm BSM}|^{2}| caligraphic_M start_POSTSUBSCRIPT roman_BSM end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT) is always subdominant for perturbative values of the BSM coupling, yDM4πless-than-or-similar-tosubscript𝑦𝐷𝑀4𝜋y_{DM}\lesssim 4\piitalic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≲ 4 italic_π. We also point out that the interference term can be negative or positive depending on the behavior of the form factors at distinct energy scales. Since the quark initiated process (qq¯tt¯𝑞¯𝑞𝑡¯𝑡q\bar{q}\to t\bar{t}italic_q over¯ start_ARG italic_q end_ARG → italic_t over¯ start_ARG italic_t end_ARG) is only affected by μsuperscript𝜇\mathcal{F}^{\mu}caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT, while the gluon process (ggtt¯𝑔𝑔𝑡¯𝑡gg\to t\bar{t}italic_g italic_g → italic_t over¯ start_ARG italic_t end_ARG) depends on both form factors, it is interesting to investigate the individual contributions from each process. For instance, in the EFT regime, we have Cg<0subscript𝐶𝑔0C_{g}<0italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT < 0, resulting in a negative interference contribution from ggtt¯𝑔𝑔𝑡¯𝑡gg\to t\bar{t}italic_g italic_g → italic_t over¯ start_ARG italic_t end_ARG.

Refer to caption
Refer to caption
Figure 4: Distributions for the tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG invariant mass. The SM (Born) contribution is shown by the filled histogram, while the solid histograms show the BSM contributions (interference term) using the full 1-loop calculation and the EFT approximation. The left plot shows the distributions for low BSM masses (mT=500subscript𝑚𝑇500m_{T}=500italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 500 GeV and mχ=400subscript𝑚𝜒400m_{\chi}=400italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 400 GeV) and the right plot for higher masses (mT=1000subscript𝑚𝑇1000m_{T}=1000italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 1000 GeV and mχ=900subscript𝑚𝜒900m_{\chi}=900italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 900 GeV). The bottom subplots show the ratio of the expected number of events for the BSM contributions to the SM one.

In Fig. 4 we show the (normalized) m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) distributions for two sets of BSM masses and yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT. The filled histogram shows the SM LO distribution, while the BSM interference term from Eq.( 12) is shown by the solid histograms. The dark red and green histograms correspond to gluon and quark initiated processes computed using the full 1-loop form factors, while the light red and green ones show the results using the EFT approximation. Since the interference can be negative, the upper plots show the absolute value of the distributions, while the lower subplots show the ratio of each BSM contribution to the SM result.

For the light mass case (left plot) we see that the EFT and 1-loop curves start to differ around m(tt¯)600similar-to𝑚𝑡¯𝑡600m(t\bar{t})\sim 600italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ∼ 600 GeV. This is expected, since the EFT approximation is clearly not valid for s=m(tt¯)mT,mχformulae-sequence𝑠𝑚𝑡¯𝑡greater-than-or-equivalent-tosubscript𝑚𝑇subscript𝑚𝜒\sqrt{s}=m(t\bar{t})\gtrsim m_{T},m_{\chi}square-root start_ARG italic_s end_ARG = italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ≳ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT. First we point out that within the EFT approximation the gluon initiated process always follows very closely the SM distribution, thus simply resulting in a rescaling of the total cross-section. For the scenario investigated here, this contribution is always negative, since Cg<0subscript𝐶𝑔0C_{g}<0italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT < 0. The behavior of the 1-loop distribution is also negative for most values of the invariant mass, except for m(tt¯)2mTsimilar-to𝑚𝑡¯𝑡2subscript𝑚𝑇m(t\bar{t})\sim 2m_{T}italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ∼ 2 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, where the distribution resembles a broad resonance and its contribution becomes positive. The contribution from the quark initiated process (qqtt¯𝑞𝑞𝑡¯𝑡qq\to t\bar{t}italic_q italic_q → italic_t over¯ start_ARG italic_t end_ARG), on the other hand, is always harder than the SM one for both the EFT and 1-loop distributions. Within the EFT approximation, however, the BSM contribution always increases with m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) (relative to the SM), while the 1-loop distribution presents a very broad enhancement around m(tt¯)=2mT=1𝑚𝑡¯𝑡2subscript𝑚𝑇1m(t\bar{t})=2m_{T}=1italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) = 2 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 1 TeV. In addition, for large invariant mass values (m(tt¯)2.5greater-than-or-equivalent-to𝑚𝑡¯𝑡2.5m(t\bar{t})\gtrsim 2.5italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ≳ 2.5 TeV), the 1-loop qqtt¯𝑞𝑞𝑡¯𝑡qq\to t\bar{t}italic_q italic_q → italic_t over¯ start_ARG italic_t end_ARG contribution becomes negative. These features can be traced back to the discussion on the form factor in the previous section. As seen in Fig. 3, \mathcal{F}caligraphic_F presents a broad bump behaviour near 5mT5subscript𝑚𝑇\sqrt{5}\,m_{T}square-root start_ARG 5 end_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT and a dominance of the imaginary part for large s𝑠sitalic_s values. As a result, the full 1-loop distribution displays an excess for the invariant mass bins close to 2mT2subscript𝑚𝑇2m_{T}2 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT and the "intermediate" bins would be the most sensitive to BSM contributions. This behavior would not be expected if we (wrongly) assumed the EFT approximation to hold, since its distributions tend to always grow with m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ). We also see that the EFT approximation considerably overestimates the signal at the tail of the distribution.

Once we consider higher BSM masses, as shown in the right plot of Fig. 4, the EFT and 1-loop distributions agree fairly well up to smT,mχsimilar-to𝑠subscript𝑚𝑇subscript𝑚𝜒\sqrt{s}\sim m_{T},m_{\chi}square-root start_ARG italic_s end_ARG ∼ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT, as expected. The broad resonant behavior of the 1-loop distributions is once again present, but now it only starts to appear at m(tt¯)5mT2.2similar-to-or-equals𝑚𝑡¯𝑡5subscript𝑚𝑇similar-to-or-equals2.2m(t\bar{t})\simeq\sqrt{5}m_{T}\simeq 2.2italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ≃ square-root start_ARG 5 end_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≃ 2.2 TeV. In this example the higher bins would be the most sensitive to the BSM contributions and the constraints are stronger than the ones expected from the EFT approximation, since the 1-loop distribution is clearly larger than the EFT one at the tail of the distribution. In addition, the gluon and quark initiated processes in the 1-loop calculation are both positive at the tail, while the gluon curve is always negative if we assume the EFT approximation, thus reducing the total BSM EFT signal.

Refer to caption
Refer to caption
Figure 5: Same as in Figure 4 for the top transverse momentum. The boosted top phase space was select according to the cuts defined in Table. 2

In addition to the top pair invariant mass, the top pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT can also be used to constrain new physics contributions ATLAS-TOPQ-2019-23 . Since below we will consider the ATLAS measurement from Ref. ATLAS-TOPQ-2019-23 and it includes fiducial phase-space cuts, in Figure 5 we show the pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT distributions after applying the ATLAS selection. All cuts applied at the particle level are listed in Table 2 and target events where one top decays hadronically and the other leptonically. The distributions shown in Fig.5 correspond to the transverse momentum of the hadronically decaying top. Although the differences between the 1-loop calculation and the EFT approximation are not so dramatic as in the m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) case, we also notice relevant differences between the two methods. In particular, the gluon initiated contributions from the 1-loop results are positive for a wide range of pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT values, while it is always negative within the EFT approximation. We also see that total 1-loop result is larger than the EFT approximation up to pT3mTsimilar-tosubscript𝑝𝑇3subscript𝑚𝑇p_{T}\sim 3m_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ∼ 3 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT. It is also important to point out that while the EFT approximation is valid for s^=m(tt¯)mT,mχformulae-sequence^𝑠𝑚𝑡¯𝑡much-less-thansubscript𝑚𝑇subscript𝑚𝜒\sqrt{\hat{s}}=m(t\bar{t})\ll m_{T},m_{\chi}square-root start_ARG over^ start_ARG italic_s end_ARG end_ARG = italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ≪ italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT, it is not so easy to identify a similar condition for the pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT distribution. This can be seen in Fig.5 (right), where the EFT approximation fails even at pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT values much smaller than mT,mχsubscript𝑚𝑇subscript𝑚𝜒m_{T},m_{\chi}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT.

Boosted Top Phase-Space
Jet Cuts n(j)>0𝑛𝑗0n(j)>0italic_n ( italic_j ) > 0
pT>36subscript𝑝𝑇36p_{T}>36italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT > 36 GeV
|η|<2.5𝜂2.5|\eta|<2.5| italic_η | < 2.5
Fat jet Cuts n(j)>0𝑛𝑗0n(j)>0italic_n ( italic_j ) > 0
pT>355subscript𝑝𝑇355p_{T}>355italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT > 355 GeV
|η|<2.0𝜂2.0|\eta|<2.0| italic_η | < 2.0
120 GeV <m<220absent𝑚220<m<220< italic_m < 220 GeV
contains one b𝑏bitalic_b
Lepton Cuts n(l)=1𝑛𝑙1n(l)=1italic_n ( italic_l ) = 1
pT>27subscript𝑝𝑇27p_{T}>27italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT > 27 GeV
|η|<2.5𝜂2.5|\eta|<2.5| italic_η | < 2.5
ΔR(l,b)<2.0Δ𝑅𝑙𝑏2.0\Delta R(l,b)<2.0roman_Δ italic_R ( italic_l , italic_b ) < 2.0
m(l,b)<180𝑚𝑙𝑏180m(l,b)<180italic_m ( italic_l , italic_b ) < 180 GeV
ΔR(l,j)>0.4Δ𝑅𝑙𝑗0.4\Delta R(l,j)>0.4roman_Δ italic_R ( italic_l , italic_j ) > 0.4
ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT Cut ETmiss>20superscriptsubscript𝐸𝑇miss20E_{T}^{\rm miss}>20italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT > 20 GeV
Table 2: Fiducial phase-space cuts applied at parton level to reproduce the boosted top phase-space considered by ATLAS in Ref. ATLAS-TOPQ-2019-23 .

4 LHC Constraints

As shown in Sec. 3.2, the BSM model considered here can have an impact on the differential top distributions. However, for masses smaller than a few TeV, the BSM states can be produced on-shell at the LHC. Therefore this scenario can potentially be constrained by: i) direct searches, i.e. searches for on-shell φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT production, and ii) indirect searches, i.e. measurements of pptt¯𝑝𝑝𝑡¯𝑡pp\to t\bar{t}italic_p italic_p → italic_t over¯ start_ARG italic_t end_ARG distributions. In this Section we mostly aim to address the following questions:

  • Can indirect searches be complementary to direct searches?

  • What is the impact of (wrongly) assuming the EFT approximation when constraining the model?

Clearly the above answers depend on the model parameters: (mT,mχ,yDM)subscript𝑚𝑇subscript𝑚𝜒subscript𝑦𝐷𝑀(m_{T},m_{\chi},y_{DM})( italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT , italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ) or (mT,ΔM=mTmχ,yDM)formulae-sequencesubscript𝑚𝑇Δ𝑀subscript𝑚𝑇subscript𝑚𝜒subscript𝑦𝐷𝑀(m_{T},\Delta M=m_{T}-m_{\chi},y_{DM})( italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , roman_Δ italic_M = italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT - italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT , italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ). For instance, for sufficiently large masses we expect the EFT results to be valid. Also, for very small yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT, the loop contributions to top pair production are suppressed and direct searches will be more sensitive.

Refer to caption
Figure 6: Leading order diagrams for the on-shell production of the colored scalar.

4.1 Direct Searches

If φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT is sufficiently light (mTless-than-or-similar-tosubscript𝑚𝑇absentm_{T}\lesssimitalic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≲ TeV), it will be copiously produced at the LHC, since it is colored. The leading order diagrams for φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT production are shown in Figure  6. The signatures generated by φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT production and decay strongly depend on the mass difference ΔM=mTmχΔ𝑀subscript𝑚𝑇subscript𝑚𝜒\Delta M=m_{T}-m_{\chi}roman_Δ italic_M = italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT - italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT. For sufficiently large ΔMΔ𝑀\Delta Mroman_Δ italic_M, the signature corresponds to on-shell tops plus missing energy (ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT), while in the compressed scenario (ΔMmtmuch-less-thanΔ𝑀subscript𝑚𝑡\Delta M\ll m_{t}roman_Δ italic_M ≪ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT) we have ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT plus b𝑏bitalic_b-jets and additional (soft) leptons and light jets. The compressed scenario tends to be more challenging, resulting in weaker constraints on mTsubscript𝑚𝑇m_{T}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT. We also point out that the signal does not dependent on the BSM coupling yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT, except for the scalar width. Although for sufficiently small widths yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT the scalar can become long-lived, in the following we assume yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT large enough so φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT always have prompt decays.

Since the LHC signatures are the same employed on stop searches, we make use of SModelS Kraml:2013mwa ; Alguero:2020grj ; Alguero:2021dig ; MahdiAltakach:2023bdn to reinterpret the ATLAS and CMS constraints on stop-neutralino simplified models and identify the most relevant analyses. The compressed scenario is particularly challenging and for very small ΔMΔ𝑀\Delta Mroman_Δ italic_M the decay products can be very soft and missed by most event selection criteria. In this region of parameter space searches for Dark Matter production, which target initial state radiation (ISR) jets plus ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT can become relevant. Therefore, in addition to the stop searches, we have considered the CMS jets plus ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT search CMS-EXO-20-004 which targets scenarios with hard jets coming from ISR. We have recast this analysis and used MadGraph5_aMC@NLO, Pythia 8.306 Pythia8 and Delphes Delphes to reinterpret the CMS results for the BSM scenario from Sec. 3.1. Although we have computed the φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT production cross-section at leading order (LO), a constant k-factor k=1.5𝑘1.5k=1.5italic_k = 1.5 was used to approximate the NNLO+NNLL result Beenakker:2016gmf .

Refer to caption
Figure 7: 95 % C.L. exclusions in the mTsubscript𝑚𝑇m_{T}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT versus ΔM=mTmχΔ𝑀subscript𝑚𝑇subscript𝑚𝜒\Delta M=m_{T}-m_{\chi}roman_Δ italic_M = italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT - italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT plane from direct searches for on/off-shell tops and missing energy. The red curve corresponds to the CMS search CMS-EXO-20-004 for ISR jets plus ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT (EXO-20-004), the green curve shows the excluded region from the CMS search CMS-SUS-20-002 for 0,1010,10 , 1 and 2 leptons plus ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT (CMS-SUS-20-002), while the orange and blue curves show the excluded regions from the ATLAS (ATLAS-SUSY-2018-12) and CMS (CMS-SUS-19-006) searches for hadronic tops plus ETmisssuperscriptsubscript𝐸𝑇missE_{T}^{\rm miss}italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_miss end_POSTSUPERSCRIPT from Refs. ATLAS-SUSY-2018-12 and  CMS-SUS-19-006 , respectively.

In Figure 7 we show the 95% C.L. excluded region in the ΔMΔ𝑀\Delta Mroman_Δ italic_M vs mTsubscript𝑚𝑇m_{T}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT plane. For large mass differences (ΔM>mWΔ𝑀subscript𝑚𝑊\Delta M>m_{W}roman_Δ italic_M > italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT), the most relevant analyses are the CMS combined stop search CMS-SUS-20-002 and the CMS search for jets and missing energy CMS-SUS-19-006 . These searches lose sensitivity to scenarios with small ΔMΔ𝑀\Delta Mroman_Δ italic_M and limits are not provided in this case, leading to the sharp cut-off seen on the green and orange curves in Fig. 7. The ATLAS search for hadronic tops plus missing energy ATLAS-SUSY-2018-12 , however, also targets the compressed scenario and is the most sensitive search in this region of parameter space, excluding mass differences down to 10-20 GeV. For even smaller mass differences, the decay products are very soft and the CMS search for ISR jets becomes relevant, as shown by the red curve in Fig. 7. Overall we see that in the highly compressed scenario scalar masses up to mT500similar-to-or-equalssubscript𝑚𝑇500m_{T}\simeq 500italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≃ 500 GeV are excluded, while for large ΔMΔ𝑀\Delta Mroman_Δ italic_M the exclusion goes up to 1.3 TeV.

4.2 Indirect searches in tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG

As discussed in Sec. 4.1, the constraints from direct searches for φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT depend strongly on the scalar-DM mass difference and exclude masses up to mT1.3similar-to-or-equalssubscript𝑚𝑇1.3m_{T}\simeq 1.3italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≃ 1.3 TeV. These constraints, however, do not depend on the BSM coupling (yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT), except for yDM1much-less-thansubscript𝑦𝐷𝑀1y_{DM}\ll 1italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≪ 1, which could render the colored scalar long-lived. On the other hand, if yDM1greater-than-or-equivalent-tosubscript𝑦𝐷𝑀1y_{DM}\gtrsim 1italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≳ 1, BSM loop contributions to pptt¯𝑝𝑝𝑡¯𝑡pp\to t\bar{t}italic_p italic_p → italic_t over¯ start_ARG italic_t end_ARG are enhanced and can become sizeable. The top distributions have been measured at high accuracy both by CMS and ATLAS and found to be in good agreement with the SM predictions. Here we follow closely the approach developed in Ref. Esser:2023fdo , which considered the top pair invariant mass measured by CMS CMS-TOP-20-001 and the top transverse momentum measured by ATLAS ATLAS-TOPQ-2019-23 to constrain the couplings of axion-like particles to the top quark.

It is well known that the tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG distributions can be significantly modified by NLO and NNLO QCD corrections Catani:2019hip ; Catani:2019iny . Therefore, for the results below, we use the corresponding SM predictions at NNLO quoted by ATLAS or CMS. However, it is beyond the scope of this work to compute the BSM signal to this level of accuracy. Nonetheless, we approximate the impact of higher order corrections on the BSM contribution (interference with the SM) using a bin-dependent reweighting factor:

ki=NSMi(NNLO)NSMi(LO)NBSMi(NNLO)kiNBSMi(LO),subscript𝑘𝑖superscriptsubscript𝑁SM𝑖NNLOsuperscriptsubscript𝑁SM𝑖LOsuperscriptsubscript𝑁BSM𝑖NNLOsimilar-to-or-equalssubscript𝑘𝑖superscriptsubscript𝑁BSM𝑖LOk_{i}=\frac{N_{\rm SM}^{i}({\rm NNLO})}{N_{\rm SM}^{i}({\rm LO})}\Rightarrow N% _{\rm BSM}^{i}({\rm NNLO})\simeq k_{i}N_{\rm BSM}^{i}({\rm LO}),italic_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = divide start_ARG italic_N start_POSTSUBSCRIPT roman_SM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_NNLO ) end_ARG start_ARG italic_N start_POSTSUBSCRIPT roman_SM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_LO ) end_ARG ⇒ italic_N start_POSTSUBSCRIPT roman_BSM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_NNLO ) ≃ italic_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT roman_BSM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_LO ) , (13)

where NSMi(LO)superscriptsubscript𝑁SM𝑖LON_{\rm SM}^{i}({\rm LO})italic_N start_POSTSUBSCRIPT roman_SM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_LO ) (NBSMi(LO)superscriptsubscript𝑁BSM𝑖LON_{\rm BSM}^{i}({\rm LO})italic_N start_POSTSUBSCRIPT roman_BSM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_LO )) is the background (signal) prediction computed at LO using MadGraph5_aMC@NLO. The reweighting factors obtained through this procedure are typically ki1.31.6similar-to-or-equalssubscript𝑘𝑖1.31.6k_{i}\simeq 1.3-1.6italic_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≃ 1.3 - 1.6 for the CMS invariant mass bins and ki0.61.8similar-to-or-equalssubscript𝑘𝑖0.61.8k_{i}\simeq 0.6-1.8italic_k start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≃ 0.6 - 1.8 for the ATLAS transverse momentum bins. With the above expressions and the covariance matrices (Cijsubscript𝐶𝑖𝑗C_{ij}italic_C start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT) provided by the experimental collaborations, we can then make use of the measured distributions (NObsisuperscriptsubscript𝑁Obs𝑖N_{\rm Obs}^{i}italic_N start_POSTSUBSCRIPT roman_Obs end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT) to compute limits on the BSM coupling yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT. Following the same procedure used by the experimental collaborations we define a χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT function as:

χ2(yDM)=i,j=binsΔiCij1Δj,superscript𝜒2subscript𝑦𝐷𝑀subscript𝑖𝑗binssubscriptΔ𝑖subscriptsuperscript𝐶1𝑖𝑗subscriptΔ𝑗\chi^{2}(y_{DM})=\sum_{i,j={\rm bins}}\Delta_{i}C^{-1}_{ij}\Delta_{j},italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_i , italic_j = roman_bins end_POSTSUBSCRIPT roman_Δ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_C start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT roman_Δ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , (14)

where Δi=[NObsiNSMi(NNLO)yDM2NBSMi(NNLO)]subscriptΔ𝑖delimited-[]superscriptsubscript𝑁Obs𝑖superscriptsubscript𝑁SM𝑖NNLOsuperscriptsubscript𝑦𝐷𝑀2superscriptsubscript𝑁BSM𝑖NNLO\Delta_{i}=\left[N_{\rm Obs}^{i}-N_{\rm SM}^{i}({\rm NNLO})-y_{DM}^{2}N_{\rm BSM% }^{i}({\rm NNLO})\right]roman_Δ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = [ italic_N start_POSTSUBSCRIPT roman_Obs end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT - italic_N start_POSTSUBSCRIPT roman_SM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_NNLO ) - italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_N start_POSTSUBSCRIPT roman_BSM end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT ( roman_NNLO ) ], so the 95% C.L. limit on yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT corresponds to Δχ2=3.84Δsuperscript𝜒23.84\Delta\chi^{2}=3.84roman_Δ italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 3.84.

Below we present results for both the EFT approach discussed in Sec. 3.1.1 and the full 1-loop form factors from Sec. 3.1.2. Although the former should only be valid at high masses, it is interesting to compare both approaches and quantify how the EFT constraints deviate from the full 1-loop calculation.

CMS m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG )

We first consider the CMS measurement CMS-TOP-20-001 of the differential m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) cross-sections at s=13𝑠13\sqrt{s}=13square-root start_ARG italic_s end_ARG = 13 TeV using the full Run 2 luminosity, =137137\mathcal{L}=137caligraphic_L = 137 fb-1. The measurement includes the full kinematic range and extends up to invariant masses of 3.5 TeV. CMS has unfolded the measured distributions and provided measurements at the parton level, which can be used to constrain BSM contributions. In addition, the covariance matrix of the measurements has been provided, which allows us to include systematic and statistical uncertainty correlations. For the SM predictions we have considered the NNLO prediction quoted by CMS computed using MATRIX Grazzini:2017mhc , but since the covariance matrix for the predictions was not given, we have not included it when computing the limits on the BSM signal.

In Fig. 8 we show the measured invariant mass distribution, the SM prediction and the SM plus BSM prediction using the 1-loop form factors or the EFT approximation. The bottom subplot shows the ratio of predictions to data. The BSM masses are mT=700subscript𝑚𝑇700m_{T}=700italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 700 GeV and mχ=690subscript𝑚𝜒690m_{\chi}=690italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 690 GeV and correspond to a compressed scenario currently beyond the reach of direct searches. The BSM coupling is considerably large, yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10, but still within the perturbative regime. As we can see, the 1-loop distribution significantly deviates from data for m(tt¯)2mT=1.4greater-than-or-equivalent-to𝑚𝑡¯𝑡2subscript𝑚𝑇1.4m(t\bar{t})\gtrsim 2m_{T}=1.4italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) ≳ 2 italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 1.4 TeV, as expected from the behavior discussed in Sec. 3.2. We also note that for the intermediate bins the EFT contribution underestimates the signal, while for the last two bins it is close to the 1-loop calculation. However, since the uncertainty in the last bins is quite high, the BSM signal is mostly constrained by the intermediate bins. As a result, the EFT approximation significantly underestimates the constraints. In particular, for the point shown in Fig. 8, we obtain yDM<7.8subscript𝑦𝐷𝑀7.8y_{DM}<7.8italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT < 7.8 at 95% C.L. using the 1-loop calculation, while the EFT approximation results in yDM<10.7subscript𝑦𝐷𝑀10.7y_{DM}<10.7italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT < 10.7. Note that the error band shown in the bottom subplot of Fig. 8 corresponds only to Ciisubscript𝐶𝑖𝑖\sqrt{C_{ii}}square-root start_ARG italic_C start_POSTSUBSCRIPT italic_i italic_i end_POSTSUBSCRIPT end_ARG. The full correlations, however, are essential for computing the limits and provide stronger constraints than assuming uncorrelated bin uncertainties.

Refer to caption
Figure 8: Invariant mass distribution, m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ), for mT=700subscript𝑚𝑇700m_{T}=700italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 700 GeV, mχ=690subscript𝑚𝜒690m_{\chi}=690italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 690 GeV and yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10. The data points show the unfolded distribution measured by CMS CMS-TOP-20-001 , while the filled histogram shows the SM prediction at NNLO CMS-TOP-20-001 . The solid histogram shows the SM plus BSM distribution computed using the 1-loop form factors, while the dashed histogram shows the same distribution within the EFT approximation. The bottom subplot shows the ratio of the measured distribution to the SM, SM plus 1-loop and SM plus EFT distributions. The band shows the uncertainties ignoring correlations, i.e. Ciisubscript𝐶𝑖𝑖\sqrt{C_{ii}}square-root start_ARG italic_C start_POSTSUBSCRIPT italic_i italic_i end_POSTSUBSCRIPT end_ARG.

ATLAS pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t )

As discussed in Sec. 3.2, the top transverse momentum can also be impacted by BSM contributions. In order to constrain the signal, we consider the ATLAS measurement ATLAS-TOPQ-2019-23 of the top pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT for semi-leptonic decaying tops. Unlike the CMS measurement discussed in the previous Section, ATLAS considers the fiducial phase-space for boosted tops, which can be approximated by the cuts listed in Table 2. The measured distribution is then unfolded to the parton level. For the SM prediction we consider the values quoted by ATLAS obtained using MadGraph5_aMC@NLO and Pythia 8 after NNLO reweighting ATLAS-TOPQ-2019-23 . The covariance matrix is also provided for the unfolded measurement, but it only includes the statistical uncertainties. However, the total systematical uncertainties for each bin is also given and we include them as a diagonal contribution to Cijsubscript𝐶𝑖𝑗C_{ij}italic_C start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT, which means we ignore correlations of systematical uncertainties333Although this is clearly an approximation, there is not enough information publicly available to properly include the correlations of the systematical uncertainties..

The measured pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT distribution is shown in Fig. 9 along with the SM prediction (filled histogram) and the total distribution (SM plus BSM) computed using the 1-loop form factors and the EFT approximation. The first plot shows the distribution for yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5 and "light" BSM masses, mT=600subscript𝑚𝑇600m_{T}=600italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 600 GeV, mχ=590subscript𝑚𝜒590m_{\chi}=590italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 590 GeV, while the second one shows the same distributions, but for heavier masses (mT=1subscript𝑚𝑇1m_{T}=1italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 1 TeV and mχ=0.9subscript𝑚𝜒0.9m_{\chi}=0.9italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 0.9 TeV) and a larger BSM coupling, yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10. For the lighter BSM masses we see that the 1-loop distribution deviates more strongly from data in the intermediate bins, as expected from the behavior seen in Fig. 5. We also see that the EFT underestimates the signal in the intermediate bins, while overestimates it for the highest bin. Due to the large statistical uncertainties at large pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ), the ATLAS measurement is mostly sensitive to deviations in the low to intermediate bins. As a result, the EFT approximation results in weaker constraints to the yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT coupling. In particular, for the signal shown in the left plot of Fig. 9, we obtain yDM<4.7subscript𝑦𝐷𝑀4.7y_{DM}<4.7italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT < 4.7 at 95% C.L. using the 1-loop calculation, while yDM<7.4subscript𝑦𝐷𝑀7.4y_{DM}<7.4italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT < 7.4 if we assume the EFT approximation. A similar behavior is also seen at larger masses, as shown by the right plot in Fig. 9. In this case, however, the EFT distribution is smaller than the 1-loop one for all bins and once again underestimates the sensitivity to new physics.

Refer to caption
Refer to caption
Figure 9: Transverse momentum distribution, pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ), for mT=600subscript𝑚𝑇600m_{T}=600italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 600 GeV, mχ=590subscript𝑚𝜒590m_{\chi}=590italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 590 GeV, yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5 (left) and mT=1subscript𝑚𝑇1m_{T}=1italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT = 1 TeV, mχ=0.9subscript𝑚𝜒0.9m_{\chi}=0.9italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT = 0.9 TeV, yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10 (right). The data points show the unfolded distribution measured by ATLAS ATLAS-TOPQ-2019-23 , while the filled histogram shows the SM prediction at NNLO from Ref. ATLAS-TOPQ-2019-23 . The solid histogram shows the SM plus BSM distribution computed using the 1-loop form factors, while the dashed histogram shows the same distribution within the EFT approximation. The bottom subplot shows the ratio of the measured distribution to the SM, SM plus 1-loop and SM plus EFT distributions. The band shows the uncertainties ignoring correlations, i.e. Ciisubscript𝐶𝑖𝑖\sqrt{C_{ii}}square-root start_ARG italic_C start_POSTSUBSCRIPT italic_i italic_i end_POSTSUBSCRIPT end_ARG.

Note that, when compared to the invariant mass distributions from Fig. 8, the pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ) measurement seems to be more sensitive to the BSM signal than the invariant mass distribution, since a larger excess is seen in the intermediate pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT bins. In addition, both measured distributions display under-fluctuations in several bins with respect to the SM prediction. As a result, the constraints on the BSM signal are stronger than expected (see Appendix C for more details).

4.3 Results

The discussion in the previous Sections showed that top measurements can be sensitive to new physics and complementary to direct searches, specially in the compressed region, ΔMmtmuch-less-thanΔ𝑀subscript𝑚𝑡\Delta M\ll m_{t}roman_Δ italic_M ≪ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, and for large BSM couplings, yDM5greater-than-or-equivalent-tosubscript𝑦𝐷𝑀5y_{DM}\gtrsim 5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≳ 5. In order to compare the constraints from direct and indirect searches, we scan over the BSM masses and compute the limits on yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT obtained from the CMS m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) measurement and the ATLAS pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ) distribution. In Figure 10 we show the region excluded by direct searches for φTsubscript𝜑𝑇\varphi_{T}italic_φ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT and the exclusion curves from indirect searches for yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5 and yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10. The regions to the left of the blue (red) curves are excluded at 95% C.L. by the CMS invariant mass (ATLAS pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT) measurement. We also show the corresponding curves obtained assuming the EFT approximation. As expected from the results in the previous Section, the pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ) measurement is more sensitive to the BSM signal and excludes masses up to mT630similar-to-or-equalssubscript𝑚𝑇630m_{T}\simeq 630italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≃ 630 GeV in the highly compressed region, if we take yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5. This exclusion goes slightly beyond the masses probed by direct searches in the compressed scenario and illustrate the complementarity between the two types of searches. For the same value of the BSM coupling, yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5, the exclusion curves obtained using the invariant mass distribution fall inside the direct search excluded region and are not shown. Once we consider yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10, the m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) measurement becomes competitive with direct searches in the compressed region and exclude masses up to mT800similar-to-or-equalssubscript𝑚𝑇800m_{T}\simeq 800italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≃ 800 GeV. But the limits obtained from the pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ) distribution are still stronger, excluding up to mT1.1similar-to-or-equalssubscript𝑚𝑇1.1m_{T}\simeq 1.1italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≃ 1.1 TeV and are complementary to direct searches even beyond the compressed region.

In Figure 10 we also display the exclusion curves obtained assuming the EFT approximation. As discussed in Sec. 3.2, the EFT regime is not valid for the range of BSM masses considered here and the energies probed by the top measurements. Indeed we see that the EFT calculation considerably underestimates the excluded regions. In particular, for yDM=5subscript𝑦𝐷𝑀5y_{DM}=5italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 5, the region excluded by top measurements falls completely inside the region already excluded by direct searches if we assume the EFT distributions.

Note that if we assume this minimal model fully explains the Dark Matter relic abundance, values of yDM3greater-than-or-equivalent-tosubscript𝑦𝐷𝑀3y_{DM}\gtrsim 3italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT ≳ 3 would likely be excluded by DM Direct Detection searches Garny:2018icg . In this case, the tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG measurements discussed here are not yet competitive to other searches. Nonetheless, once more LHC data is collected, the constraints obtained using the full 1-loop calculations could become relevant. The same is not true for the EFT analysis, since its constraints are too conservative. Therefore considering the 1-loop results is indeed essential for properly assessing the impact of top measurements to Dark Matter models.

Refer to caption
Figure 10: Comparison of the exclusion curves from direct stop searches and from top measurements. The solid lines show the exclusion obtained using the full 1-loop calculation, while the dashed lines correspond to the EFT approximation for distinct values of the BSM coupling yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT. The limits obtained using the CMS invariant mass measurement from Ref. CMS-TOP-20-001 are shown by the blue curves, while the ones obtained from the ATLAS pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT measurement ATLAS-TOPQ-2019-23 are shown in red.

All the results presented in Fig.10 rely on tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG distributions and do not exploit the full range of LHC and low energy data. On the other hand, global fits of SMEFT operators include a large number of observables. Although the validity of SMEFT is limited to the high mass region, it is still interesting to compare the results from the global fit in Ref. Ellis:2020unq and shown in Table 1 with the limits from Fig.10. First we point out that the strongest limit from the global fit for a single operator gives |CtG|/Λ2=|Cg|<0.1subscript𝐶𝑡𝐺superscriptΛ2subscript𝐶𝑔0.1|C_{tG}|/\Lambda^{2}=|C_{g}|<0.1| italic_C start_POSTSUBSCRIPT italic_t italic_G end_POSTSUBSCRIPT | / roman_Λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = | italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT | < 0.1 TeV-2, translating to mTgreater-than-or-equivalent-tosubscript𝑚𝑇absentm_{T}\gtrsimitalic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≳ 400 GeV for yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10, while the corresponding constraint we have obtained using the EFT approximation is mTgreater-than-or-equivalent-tosubscript𝑚𝑇absentm_{T}\gtrsimitalic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ≳ 800 GeV. Note, though, that by considering a single Wilson coefficient in the SMEFT analysis we are not making use of the correlations between the gluon-top 𝒪tGsubscript𝒪𝑡𝐺\mathcal{O}_{\scriptscriptstyle tG}caligraphic_O start_POSTSUBSCRIPT italic_t italic_G end_POSTSUBSCRIPT and four-fermion 𝒪tX(8)superscriptsubscript𝒪𝑡𝑋8\mathcal{O}_{\scriptscriptstyle tX}^{\scriptscriptstyle(8)}caligraphic_O start_POSTSUBSCRIPT italic_t italic_X end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 8 ) end_POSTSUPERSCRIPT operators, which explains the weaker bound. Second, we have shown that the EFT result tends to underestimate the constraints. Therefore these two factors considerably enhances the SM measurements sensitivity to new physics when compared to a global SMEFT analysis.

5 Conclusions

In this paper we have presented a new way to re-interpret SM top measurements which can access light new physics scales not suitable for the SMEFT framework, providing a path to go beyond the Top EFT approach. In particular, we propose to consider new physics scenarios which produce loop-induced signatures in SM final states, motivated by the existence of Dark Matter. In order to properly assess the sensitivity of top measurements to the BSM signal, we have computed the leading one-loop BSM contributions to top pair production through the use of form-factor effective couplings, which fully capture the BSM one-loop effects and can differ significantly from the behavior of SMEFT operators. We have found that, while the use of Top EFT operators produce an excess in the tail of distributions, the behavior of the loop calculation can be similar to a very broad bump. Therefore, in order to fully capture the effects of new physics one cannot simply re-interpret current SMEFT analyses, but should extend them making use of different kinematic regions other than distribution tails. Our work shows that this re-interpretation is possible and would make use of the same differential measurements considered by the experimental collaborations, just focusing in a different region of phase space. This procedure could be generalized to other SMEFT sectors, like diboson and Higgs observables.

In addition, using the information of the full form-factor, SM measurements can be more sensitive to new physics than expected by the usual SMEFT analysis. We have shown that, for sufficiently large BSM couplings, this enhancement in sensitivity renders the SM measurements complementary to direct searches, allowing to extend the excluded region of the BSM parameter space.

We point out that the results presented here are mostly intended to illustrate the potential of considering the one-loop form factors when modeling new physics contributions to SM measurements as opposite to the standard EFT approach. A precise determination of the constraints from indirect searches would require the inclusion of other top measurements, the calculation of QCD NLO corrections to the signal and more detailed information about the SM uncertainties and their correlations, which we leave for future work. Finally, the broad resonant behavior of the signal and the negative interference at high energy bins could also be better exploited to constrain new physics. In particular, ratios of intermediate and high invariant mass bins could in principle reduce the systematical uncertainties and enhance the signal sensitivity.

Acknowledgements.
We would like to thank Maeve Madigan for her help with the reinterpretation of the ATLAS and CMS top analyses and Johnathon Gargalionis for his help with the EFT matching procedure. A.L. is supported by FAPESP grants no. 2018/25225-9 and 2021/01089-1. The research of VS is supported by the Generalitat Valenciana PROMETEO/2021/083 and the Ministerio de Ciencia e Innovacion PID2020-113644GB-I00.

Appendix A Feynman Diagrams

In this Section we list all the relevant Feynman diagrams required for computing the BSM contributions to tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG production at the LHC at leading order in the BSM coupling yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT. The diagrams are shown in Figure 11 and are divided into quark initiated (first row) and gluon initiated (second, third and fourth rows) processes. In the left column we show the diagrams used for the 1-loop calculation, while the right column shows the equivalent ones used within the EFT approximation. The diagrams which correspond to a top/anti-top permutation of other diagrams are not explicitly shown, but are indicated by (t¯t)¯𝑡𝑡\left(\bar{t}\leftrightarrow t\right)( over¯ start_ARG italic_t end_ARG ↔ italic_t ). The counter-terms diagrams needed for regularizing the divergent diagrams are not shown.

Refer to caption
Figure 11: Loop diagrams and their EFT equivalents contributing to t¯t¯𝑡𝑡\bar{t}tover¯ start_ARG italic_t end_ARG italic_t production at leading order in the BSM coupling yDMsubscript𝑦𝐷𝑀y_{DM}italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT.

Appendix B Form Factors and EFT matching

In order to verify the expression for the form factor μsuperscript𝜇\mathcal{F}^{\mu}caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT defined in Eq.(9) and the EFT coefficients Cg,Cqsubscript𝐶𝑔subscript𝐶𝑞C_{g},C_{q}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT defined in Eqs.(4) and (4), we compute the amplitude for the quark initiated process, qq¯tt¯𝑞¯𝑞𝑡¯𝑡q\bar{q}\to t\bar{t}italic_q over¯ start_ARG italic_q end_ARG → italic_t over¯ start_ARG italic_t end_ARG, using both the 1-loop form factors and the EFT operators. The corresponding diagrams are shown by the first row in Figure 11. In the EFT limit, mT,mχmt,sformulae-sequencemuch-greater-thansubscript𝑚𝑇subscript𝑚𝜒subscript𝑚𝑡𝑠m_{T},m_{\chi}\gg m_{t},\sqrt{s}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT ≫ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , square-root start_ARG italic_s end_ARG, we should obtain the same result for both approaches, thus validating our implementation. A similar exercise could in principle be done for the ggtt¯𝑔𝑔𝑡¯𝑡gg\to t\bar{t}italic_g italic_g → italic_t over¯ start_ARG italic_t end_ARG process. However, due to the large number of diagrams (see Figure 11) and the complexity of the form factor, the validation of μνsuperscript𝜇𝜈\mathcal{F}^{\mu\nu}caligraphic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT has only been done numerically.

We start with the 1-loop calculation using form factors. Note that for the qq¯tt¯𝑞¯𝑞𝑡¯𝑡q\bar{q}\to t\bar{t}italic_q over¯ start_ARG italic_q end_ARG → italic_t over¯ start_ARG italic_t end_ARG process only the top-top-gluon form factor (μsuperscript𝜇\mathcal{F}^{\mu}caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT) is relevant and the top momenta appearing in μsuperscript𝜇\mathcal{F}^{\mu}caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT are on-shell, p12=p22=mt2superscriptsubscript𝑝12superscriptsubscript𝑝22superscriptsubscript𝑚𝑡2p_{1}^{2}=p_{2}^{2}=m_{t}^{2}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. In this case the form factor defined in Eq.(9) can be simplified, since C1=C1(p12,s,p22)=C1(mt2,s,mt2)=C2subscript𝐶1subscript𝐶1superscriptsubscript𝑝12𝑠superscriptsubscript𝑝22subscript𝐶1superscriptsubscript𝑚𝑡2𝑠superscriptsubscript𝑚𝑡2subscript𝐶2C_{1}=C_{1}(p_{1}^{2},s,p_{2}^{2})=C_{1}(m_{t}^{2},s,m_{t}^{2})=C_{2}italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_s , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_s , italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and similarly C22=C11subscript𝐶22subscript𝐶11C_{22}=C_{11}italic_C start_POSTSUBSCRIPT 22 end_POSTSUBSCRIPT = italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT. Using these results we obtain:

μ=superscript𝜇absent\displaystyle\mathcal{F}^{\mu}=caligraphic_F start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT = iTAPR[(C1+2C11)(1p1μ+2p2μ)(C1+2C12)(1p2μ+2p1μ)]𝑖superscript𝑇𝐴subscript𝑃𝑅delimited-[]subscript𝐶12subscript𝐶11subscriptitalic-p̸1superscriptsubscript𝑝1𝜇subscriptitalic-p̸2superscriptsubscript𝑝2𝜇subscript𝐶12subscript𝐶12subscriptitalic-p̸1superscriptsubscript𝑝2𝜇subscriptitalic-p̸2superscriptsubscript𝑝1𝜇\displaystyle iT^{A}P_{R}\left[\left(C_{1}+2C_{11}\right)\left(\not{p}_{1}p_{1% }^{\mu}+\not{p}_{2}p_{2}^{\mu}\right)-\left(C_{1}+2C_{12}\right)\left(\not{p}_% {1}p_{2}^{\mu}+\not{p}_{2}p_{1}^{\mu}\right)\right]italic_i italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT [ ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ) ( italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT + italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ) - ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) ( italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT + italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ) ]
+2iTAPRγμ(C00+δR)+2iTAPLγμδL2𝑖superscript𝑇𝐴subscript𝑃𝑅superscript𝛾𝜇subscript𝐶00subscript𝛿𝑅2𝑖superscript𝑇𝐴subscript𝑃𝐿superscript𝛾𝜇subscript𝛿𝐿\displaystyle+2iT^{A}P_{R}\gamma^{\mu}\left(C_{00}+\delta_{R}\right)+2iT^{A}P_% {L}\gamma^{\mu}\delta_{L}+ 2 italic_i italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ( italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) + 2 italic_i italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT (15)

where Ci,Cijsubscript𝐶𝑖subscript𝐶𝑖𝑗C_{i},C_{ij}italic_C start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_C start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT are the triangular scalar loop functions and δL,Rsubscript𝛿𝐿𝑅\delta_{L,R}italic_δ start_POSTSUBSCRIPT italic_L , italic_R end_POSTSUBSCRIPT are the counter-terms. Using NLOCT Degrande:2014vpa , dimensional regularization and the on-shell renormalization conditions we obtain444The expressions in Eq.(16) for the counter-term are valid for mT>mχ+mtsubscript𝑚𝑇subscript𝑚𝜒subscript𝑚𝑡m_{T}>m_{\chi}+m_{t}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT > italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT + italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT. Similar expressions can be found for the compressed region and have been used in our results for scenarios with mT<mχ+mtsubscript𝑚𝑇subscript𝑚𝜒subscript𝑚𝑡m_{T}<m_{\chi}+m_{t}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT < italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT + italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT.:

δL=subscript𝛿𝐿absent\displaystyle\delta_{L}=italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT = 1128π4xt2λ{2λ[xt(xt+2x+log(x)2)(x1)2log(x)]\displaystyle\frac{1}{128\pi^{4}x_{t}^{2}\sqrt{\lambda}}\biggl{\{}2\sqrt{% \lambda}\left[x_{t}\left(x_{t}+2x+\log(x)-2\right)-(x-1)^{2}\log(x)\right]divide start_ARG 1 end_ARG start_ARG 128 italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT square-root start_ARG italic_λ end_ARG end_ARG { 2 square-root start_ARG italic_λ end_ARG [ italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + 2 italic_x + roman_log ( italic_x ) - 2 ) - ( italic_x - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_log ( italic_x ) ]
4[xt(xt+x2+x2)(x1)3]log(λxt+x+12x)}\displaystyle\qquad-4\left[x_{t}\left(x_{t}+x^{2}+x-2\right)-(x-1)^{3}\right]% \log\left(\frac{\sqrt{\lambda}-x_{t}+x+1}{2\sqrt{x}}\right)\biggr{\}}- 4 [ italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_x - 2 ) - ( italic_x - 1 ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ] roman_log ( divide start_ARG square-root start_ARG italic_λ end_ARG - italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + italic_x + 1 end_ARG start_ARG 2 square-root start_ARG italic_x end_ARG end_ARG ) }
δR=subscript𝛿𝑅absent\displaystyle\delta_{R}=italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = (x1xt)128π4xt2λ{2[(xt1)2+(x2)x]log(λxt+x+12x)\displaystyle\frac{\left(x-1-x_{t}\right)}{128\pi^{4}x_{t}^{2}\sqrt{\lambda}}% \biggl{\{}2\left[\left(x_{t}-1\right)^{2}+(x-2)x\right]\log\left(\frac{\sqrt{% \lambda}-x_{t}+x+1}{2\sqrt{x}}\right)divide start_ARG ( italic_x - 1 - italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) end_ARG start_ARG 128 italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT square-root start_ARG italic_λ end_ARG end_ARG { 2 [ ( italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_x - 2 ) italic_x ] roman_log ( divide start_ARG square-root start_ARG italic_λ end_ARG - italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + italic_x + 1 end_ARG start_ARG 2 square-root start_ARG italic_x end_ARG end_ARG )
+λ[2xt(xt+x1)log(x)]}164π4[1ϵ+log(μ2mT2)]\displaystyle\qquad+\sqrt{\lambda}\left[2x_{t}-\left(x_{t}+x-1\right)\log(x)% \right]\biggr{\}}-\frac{1}{64\pi^{4}}\left[\frac{1}{\epsilon}+\log\left(\frac{% \mu^{2}}{m_{T}^{2}}\right)\right]+ square-root start_ARG italic_λ end_ARG [ 2 italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - ( italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + italic_x - 1 ) roman_log ( italic_x ) ] } - divide start_ARG 1 end_ARG start_ARG 64 italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG italic_ϵ end_ARG + roman_log ( divide start_ARG italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ] (16)

where x=mχ2/mT2,xt=mt2/mT2formulae-sequence𝑥superscriptsubscript𝑚𝜒2superscriptsubscript𝑚𝑇2subscript𝑥𝑡superscriptsubscript𝑚𝑡2superscriptsubscript𝑚𝑇2x=m_{\chi}^{2}/m_{T}^{2},x_{t}=m_{t}^{2}/m_{T}^{2}italic_x = italic_m start_POSTSUBSCRIPT italic_χ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and λ=1+x2+xt22xxt2xt2x𝜆1superscript𝑥2superscriptsubscript𝑥𝑡22𝑥subscript𝑥𝑡2subscript𝑥𝑡2𝑥\lambda=1+x^{2}+x_{t}^{2}-2xx_{t}-2x_{t}-2xitalic_λ = 1 + italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_x italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - 2 italic_x start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - 2 italic_x. Both the divergent (1/ϵ1italic-ϵ1/\epsilon1 / italic_ϵ) and the log(μ2/mT2)superscript𝜇2superscriptsubscript𝑚𝑇2\log(\mu^{2}/m_{T}^{2})roman_log ( italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) terms in δRsubscript𝛿𝑅\delta_{R}italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT cancel with the corresponding terms from C00subscript𝐶00C_{00}italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT, which is the only divergent loop integral in Eq.(15).

Using the results above, the amplitude for the first diagram in Figure 11 becomes:

loop=subscriptloopabsent\displaystyle\mathcal{M}_{\rm loop}=caligraphic_M start_POSTSUBSCRIPT roman_loop end_POSTSUBSCRIPT = iπ2yDM2gs21sTAv¯q¯γμuqu¯t(iμ)vt¯𝑖superscript𝜋2superscriptsubscript𝑦𝐷𝑀2superscriptsubscript𝑔𝑠21𝑠superscript𝑇𝐴subscript¯𝑣¯𝑞superscript𝛾𝜇subscript𝑢𝑞subscript¯𝑢𝑡𝑖subscript𝜇subscript𝑣¯𝑡\displaystyle-i\pi^{2}y_{DM}^{2}g_{s}^{2}\frac{1}{s}T^{A}\bar{v}_{\bar{q}}% \gamma^{\mu}u_{q}\bar{u}_{t}\left(-i\mathcal{F}_{\mu}\right)v_{\bar{t}}- italic_i italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_s end_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( - italic_i caligraphic_F start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ) italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT
=\displaystyle== iπ2yDM2gs21sTATA{2v¯q¯γμuq[(C00+δR)u¯tγμPRvt¯+δLu¯tγμPLvt¯]\displaystyle-i\pi^{2}y_{DM}^{2}g_{s}^{2}\frac{1}{s}T^{A}\,T^{A}\biggl{\{}2% \bar{v}_{\bar{q}}\gamma^{\mu}u_{q}\left[\left(C_{00}+\delta_{R}\right)\bar{u}_% {t}\gamma_{\mu}P_{R}v_{\bar{t}}+\delta_{L}\bar{u}_{t}\gamma_{\mu}P_{L}v_{\bar{% t}}\right]- italic_i italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_s end_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT { 2 over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ ( italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ]
+v¯q¯1uq[(C1+2C11)u¯t1PRvt¯(C1+2C12)u¯t2PRvt¯]subscript¯𝑣¯𝑞subscriptitalic-p̸1subscript𝑢𝑞delimited-[]subscript𝐶12subscript𝐶11subscript¯𝑢𝑡subscriptitalic-p̸1subscript𝑃𝑅subscript𝑣¯𝑡subscript𝐶12subscript𝐶12subscript¯𝑢𝑡subscriptitalic-p̸2subscript𝑃𝑅subscript𝑣¯𝑡\displaystyle\quad+\bar{v}_{\bar{q}}\not{p}_{1}u_{q}\left[\left(C_{1}+2C_{11}% \right)\bar{u}_{t}\not{p}_{1}P_{R}v_{\bar{t}}-\left(C_{1}+2C_{12}\right)\bar{u% }_{t}\not{p}_{2}P_{R}v_{\bar{t}}\right]+ over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT - ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ]
+v¯q¯2uq[(C1+2C11)u¯t2PRvt¯(C1+2C12)u¯t1PRvt¯]}\displaystyle\quad+\bar{v}_{\bar{q}}\not{p}_{2}u_{q}\left[\left(C_{1}+2C_{11}% \right)\bar{u}_{t}\not{p}_{2}P_{R}v_{\bar{t}}-\left(C_{1}+2C_{12}\right)\bar{u% }_{t}\not{p}_{1}P_{R}v_{\bar{t}}\right]\biggr{\}}+ over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT - ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ] } (17)

where uq(t)subscript𝑢𝑞𝑡u_{q(t)}italic_u start_POSTSUBSCRIPT italic_q ( italic_t ) end_POSTSUBSCRIPT is the quark (top) spinor and vq¯(t¯)subscript𝑣¯𝑞¯𝑡v_{\bar{q}(\bar{t})}italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG ( over¯ start_ARG italic_t end_ARG ) end_POSTSUBSCRIPT is the anti-quark (anti-top) spinor. The color indices for the Gell-Mann matrices TAsuperscript𝑇𝐴T^{A}italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT are contracted with the spinors and are not explicitly shown for simplicity.

The amplitude can be greatly simplified using the on-shell relations:

up=mtup,italic-p̸subscript𝑢𝑝subscript𝑚𝑡subscript𝑢𝑝\displaystyle\not{p}u_{p}=m_{t}u_{p},italic_p̸ italic_u start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , u¯p=mtu¯psubscript¯𝑢𝑝italic-p̸subscript𝑚𝑡subscript¯𝑢𝑝\displaystyle\quad\bar{u}_{p}\not{p}=m_{t}\bar{u}_{p}over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_p̸ = italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT
vp=mtvp,italic-p̸subscript𝑣𝑝subscript𝑚𝑡subscript𝑣𝑝\displaystyle\not{p}v_{p}=-m_{t}v_{p},italic_p̸ italic_v start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = - italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , v¯p=mtv¯psubscript¯𝑣𝑝italic-p̸subscript𝑚𝑡subscript¯𝑣𝑝\displaystyle\quad\bar{v}_{p}\not{p}=-m_{t}\bar{v}_{p}over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_p̸ = - italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT (18)

and PR=PLitalic-p̸subscript𝑃𝑅subscript𝑃𝐿italic-p̸\not{p}P_{R}=P_{L}\not{p}italic_p̸ italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = italic_P start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_p̸. Applying these relations to Eq.(17) we obtain:

loop=subscriptloopabsent\displaystyle\mathcal{M}_{\rm loop}=caligraphic_M start_POSTSUBSCRIPT roman_loop end_POSTSUBSCRIPT = iπ2yDM2gs21sTATA{2v¯q¯γμuq[(C00+δR)u¯tγμPRvt¯+δLu¯tγμPLvt¯]\displaystyle-i\pi^{2}y_{DM}^{2}g_{s}^{2}\frac{1}{s}T^{A}\,T^{A}\biggl{\{}2% \bar{v}_{\bar{q}}\gamma^{\mu}u_{q}\left[\left(C_{00}+\delta_{R}\right)\bar{u}_% {t}\gamma_{\mu}P_{R}v_{\bar{t}}+\delta_{L}\bar{u}_{t}\gamma_{\mu}P_{L}v_{\bar{% t}}\right]- italic_i italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_s end_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT { 2 over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ ( italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ]
+mt(C11C12)v¯q¯(1+2)uqu¯t(PRPL)vt¯subscript𝑚𝑡subscript𝐶11subscript𝐶12subscript¯𝑣¯𝑞subscriptitalic-p̸1subscriptitalic-p̸2subscript𝑢𝑞subscript¯𝑢𝑡subscript𝑃𝑅subscript𝑃𝐿subscript𝑣¯𝑡\displaystyle\qquad+m_{t}\left(C_{11}-C_{12}\right)\bar{v}_{\bar{q}}\left(\not% {p}_{1}+\not{p}_{2}\right)u_{q}\bar{u}_{t}\left(P_{R}-P_{L}\right)v_{\bar{t}}+ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT - italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT - italic_P start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT
+mt(C1+C11+C12)v¯q¯(12)uqu¯tvt¯}\displaystyle\qquad+m_{t}\left(C_{1}+C_{11}+C_{12}\right)\bar{v}_{\bar{q}}% \left(\not{p}_{1}-\not{p}_{2}\right)u_{q}\bar{u}_{t}v_{\bar{t}}\biggr{\}}+ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT } (19)

Finally, using momentum conservation, 1+2=1+2subscriptitalic-p̸1subscriptitalic-p̸2subscriptitalic-k̸1subscriptitalic-k̸2\not{p}_{1}+\not{p}_{2}=\not{k}_{1}+\not{k}_{2}italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_k̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_k̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, where k1,2subscript𝑘12k_{1,2}italic_k start_POSTSUBSCRIPT 1 , 2 end_POSTSUBSCRIPT are the initial state momenta and the on-shell relations for the massless incoming quarks we have: v¯q¯(1+2)uq=0subscript¯𝑣¯𝑞subscriptitalic-p̸1subscriptitalic-p̸2subscript𝑢𝑞0\bar{v}_{\bar{q}}\left(\not{p}_{1}+\not{p}_{2}\right)u_{q}=0over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = 0. Hence:

loop=subscriptloopabsent\displaystyle\mathcal{M}_{\rm loop}=caligraphic_M start_POSTSUBSCRIPT roman_loop end_POSTSUBSCRIPT = iπ2yDM2gs21sTATA{\displaystyle-i\pi^{2}y_{DM}^{2}g_{s}^{2}\frac{1}{s}T^{A}\,T^{A}\biggl{\{}- italic_i italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_s end_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT {
v¯q¯γμuq[2(C00+δRδL)u¯tγμPRvt¯+2δLu¯tγμvt¯]subscript¯𝑣¯𝑞superscript𝛾𝜇subscript𝑢𝑞delimited-[]2subscript𝐶00subscript𝛿𝑅subscript𝛿𝐿subscript¯𝑢𝑡subscript𝛾𝜇subscript𝑃𝑅subscript𝑣¯𝑡2subscript𝛿𝐿subscript¯𝑢𝑡subscript𝛾𝜇subscript𝑣¯𝑡\displaystyle\bar{v}_{\bar{q}}\gamma^{\mu}u_{q}\left[2\left(C_{00}+\delta_{R}-% \delta_{L}\right)\bar{u}_{t}\gamma_{\mu}P_{R}v_{\bar{t}}+2\delta_{L}\bar{u}_{t% }\gamma_{\mu}v_{\bar{t}}\right]over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ 2 ( italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT - italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT + 2 italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ]
+mt(C1+C11+C12)v¯q¯(12)uqu¯tvt¯}\displaystyle+m_{t}\left(C_{1}+C_{11}+C_{12}\right)\bar{v}_{\bar{q}}\left(\not% {p}_{1}-\not{p}_{2}\right)u_{q}\bar{u}_{t}v_{\bar{t}}\biggr{\}}+ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT } (20)

where we have used PL=1PRsubscript𝑃𝐿1subscript𝑃𝑅P_{L}=1-P_{R}italic_P start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT = 1 - italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT.

The above result corresponds to the full 1-loop calculation for the qq¯tt¯𝑞¯𝑞𝑡¯𝑡q\bar{q}\to t\bar{t}italic_q over¯ start_ARG italic_q end_ARG → italic_t over¯ start_ARG italic_t end_ARG process. In order to compare it to the EFT approximation, we compute the same process, but now using the EFT lagrangian defined in Eq.(2) and the diagrams shown in the first row (right column) of Figure 11. The EFT amplitude in this case is simply:

EFT=subscriptEFTabsent\displaystyle\mathcal{M}_{\rm EFT}=caligraphic_M start_POSTSUBSCRIPT roman_EFT end_POSTSUBSCRIPT = 2igsmtsTATACg[2mt(v¯q¯γμuq)(u¯tγμvt¯)+v¯q¯(21)uq(u¯tvt¯)]2𝑖subscript𝑔𝑠subscript𝑚𝑡𝑠superscript𝑇𝐴superscript𝑇𝐴subscript𝐶𝑔delimited-[]2subscript𝑚𝑡subscript¯𝑣¯𝑞superscript𝛾𝜇subscript𝑢𝑞subscript¯𝑢𝑡subscript𝛾𝜇subscript𝑣¯𝑡subscript¯𝑣¯𝑞subscriptitalic-p̸2subscriptitalic-p̸1subscript𝑢𝑞subscript¯𝑢𝑡subscript𝑣¯𝑡\displaystyle-2ig_{s}\frac{m_{t}}{s}T^{A}\,T^{A}C_{g}\left[2m_{t}\left(\bar{v}% _{\bar{q}}\gamma^{\mu}u_{q}\right)\left(\bar{u}_{t}\gamma_{\mu}v_{\bar{t}}% \right)+\bar{v}_{\bar{q}}\left(\not{p}_{2}-\not{p}_{1}\right)u_{q}\left(\bar{u% }_{t}v_{\bar{t}}\right)\right]- 2 italic_i italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT divide start_ARG italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT [ 2 italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ) ( over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ) + over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ) ]
iTATACq(v¯q¯γμuq)(u¯tγμPRvt¯)𝑖superscript𝑇𝐴superscript𝑇𝐴subscript𝐶𝑞subscript¯𝑣¯𝑞superscript𝛾𝜇subscript𝑢𝑞subscript¯𝑢𝑡subscript𝛾𝜇subscript𝑃𝑅subscript𝑣¯𝑡\displaystyle-iT^{A}\,T^{A}C_{q}\left(\bar{v}_{\bar{q}}\gamma^{\mu}u_{q}\right% )\left(\bar{u}_{t}\gamma_{\mu}P_{R}v_{\bar{t}}\right)- italic_i italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ) ( over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ) (21)

where Cq,Cgsubscript𝐶𝑞subscript𝐶𝑔C_{q},C_{g}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT , italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT are the EFT coefficients defined in Eqs.(4) and (4). Before we can compare the EFT result with the 1-loop calculation it is useful to write the above amplitude as:

EFT=subscriptEFTabsent\displaystyle\mathcal{M}_{\rm EFT}=caligraphic_M start_POSTSUBSCRIPT roman_EFT end_POSTSUBSCRIPT = iπ2yDM2gs21sTATA{\displaystyle-i\pi^{2}y_{DM}^{2}g_{s}^{2}\frac{1}{s}T^{A}\,T^{A}\biggl{\{}% \vphantom{\left[\left(\not{p}_{1}\right)\right]}- italic_i italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_s end_ARG italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT {
v¯q¯γμuq[2(sCq2π2gs2yDM2)u¯tγμPRvt¯+2(2mt2Cgπ2gsyDM2)u¯tγμvt¯]subscript¯𝑣¯𝑞superscript𝛾𝜇subscript𝑢𝑞delimited-[]2𝑠subscript𝐶𝑞2superscript𝜋2superscriptsubscript𝑔𝑠2superscriptsubscript𝑦𝐷𝑀2subscript¯𝑢𝑡subscript𝛾𝜇subscript𝑃𝑅subscript𝑣¯𝑡22superscriptsubscript𝑚𝑡2subscript𝐶𝑔superscript𝜋2subscript𝑔𝑠superscriptsubscript𝑦𝐷𝑀2subscript¯𝑢𝑡subscript𝛾𝜇subscript𝑣¯𝑡\displaystyle\bar{v}_{\bar{q}}\gamma^{\mu}u_{q}\left[2\left(\frac{sC_{q}}{2\pi% ^{2}g_{s}^{2}y_{DM}^{2}}\right)\bar{u}_{t}\gamma_{\mu}P_{R}v_{\bar{t}}+2\left(% \frac{2m_{t}^{2}C_{g}}{\pi^{2}g_{s}y_{DM}^{2}}\right)\bar{u}_{t}\gamma_{\mu}v_% {\bar{t}}\right]over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT italic_γ start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ 2 ( divide start_ARG italic_s italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT + 2 ( divide start_ARG 2 italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT ]
+mt(2Cgπ2gsyDM2)v¯q¯(12)uqu¯tvt¯}\displaystyle+m_{t}\left(\frac{-2C_{g}}{\pi^{2}g_{s}y_{DM}^{2}}\right)\bar{v}_% {\bar{q}}\left(\not{p}_{1}-\not{p}_{2}\right)u_{q}\bar{u}_{t}v_{\bar{t}}\biggr% {\}}+ italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( divide start_ARG - 2 italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) over¯ start_ARG italic_v end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_p̸ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_p̸ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_u start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT over¯ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT over¯ start_ARG italic_t end_ARG end_POSTSUBSCRIPT } (22)

Comparing Eqs.(20) and (22) we see that, in the EFT limit, we must have:

C1+C11+C12=subscript𝐶1subscript𝐶11subscript𝐶12absent\displaystyle C_{1}+C_{11}+C_{12}=italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT = (2Cgπ2yDM2gs)2subscript𝐶𝑔superscript𝜋2superscriptsubscript𝑦𝐷𝑀2subscript𝑔𝑠\displaystyle\left(\frac{-2C_{g}}{\pi^{2}y_{DM}^{2}g_{s}}\right)( divide start_ARG - 2 italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG )
C00+δRδL=subscript𝐶00subscript𝛿𝑅subscript𝛿𝐿absent\displaystyle C_{00}+\delta_{R}-\delta_{L}=italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT - italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT = (sCq2π2gs2yDM2)𝑠subscript𝐶𝑞2superscript𝜋2superscriptsubscript𝑔𝑠2superscriptsubscript𝑦𝐷𝑀2\displaystyle\left(\frac{sC_{q}}{2\pi^{2}g_{s}^{2}y_{DM}^{2}}\right)( divide start_ARG italic_s italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG )
δL=subscript𝛿𝐿absent\displaystyle\delta_{L}=italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT = (2mt2Cgπ2yDM2gs)2superscriptsubscript𝑚𝑡2subscript𝐶𝑔superscript𝜋2superscriptsubscript𝑦𝐷𝑀2subscript𝑔𝑠\displaystyle\left(\frac{2m_{t}^{2}C_{g}}{\pi^{2}y_{DM}^{2}g_{s}}\right)( divide start_ARG 2 italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG ) (23)

In order to verify the above relations we must expand the loop functions Ci,Cijsubscript𝐶𝑖subscript𝐶𝑖𝑗C_{i},C_{ij}italic_C start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_C start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT and the counter-terms δL,Rsubscript𝛿𝐿𝑅\delta_{L,R}italic_δ start_POSTSUBSCRIPT italic_L , italic_R end_POSTSUBSCRIPT to leading order in s,mt𝑠subscript𝑚𝑡s,m_{t}italic_s , italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, which gives:

C1+C11+C12=1192π41mT21(x1)4[16x+3x2+2x36x2log(x)]+𝒪(1mT4)subscript𝐶1subscript𝐶11subscript𝐶121192superscript𝜋41superscriptsubscript𝑚𝑇21superscript𝑥14delimited-[]16𝑥3superscript𝑥22superscript𝑥36superscript𝑥2𝑥𝒪1superscriptsubscript𝑚𝑇4\displaystyle C_{1}+C_{11}+C_{12}=\frac{1}{192\pi^{4}}\frac{1}{m_{T}^{2}}\frac% {1}{(x-1)^{4}}\left[1-6x+3x^{2}+2x^{3}-6x^{2}\log(x)\right]+\mathcal{O}(\frac{% 1}{m_{T}^{4}})italic_C start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 11 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 192 italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG ( italic_x - 1 ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG [ 1 - 6 italic_x + 3 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 6 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_log ( italic_x ) ] + caligraphic_O ( divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG )
C00+δRδL=11152π4smT21(x1)4[29x+18x211x3+6x3log(x)]+𝒪(1mT4)subscript𝐶00subscript𝛿𝑅subscript𝛿𝐿11152superscript𝜋4𝑠superscriptsubscript𝑚𝑇21superscript𝑥14delimited-[]29𝑥18superscript𝑥211superscript𝑥36superscript𝑥3𝑥𝒪1superscriptsubscript𝑚𝑇4\displaystyle C_{00}+\delta_{R}-\delta_{L}=\frac{1}{1152\pi^{4}}\frac{s}{m_{T}% ^{2}}\frac{1}{(x-1)^{4}}\left[2-9x+18x^{2}-11x^{3}+6x^{3}\log(x)\right]+% \mathcal{O}(\frac{1}{m_{T}^{4}})italic_C start_POSTSUBSCRIPT 00 end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT - italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 1152 italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_s end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG ( italic_x - 1 ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG [ 2 - 9 italic_x + 18 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 11 italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + 6 italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT roman_log ( italic_x ) ] + caligraphic_O ( divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG )
δL=1192π4mt2mT21(x1)4[16x+3x2+2x36x2log(x)]+𝒪(1mT4)subscript𝛿𝐿1192superscript𝜋4superscriptsubscript𝑚𝑡2superscriptsubscript𝑚𝑇21superscript𝑥14delimited-[]16𝑥3superscript𝑥22superscript𝑥36superscript𝑥2𝑥𝒪1superscriptsubscript𝑚𝑇4\displaystyle\delta_{L}=\frac{-1}{192\pi^{4}}\frac{m_{t}^{2}}{m_{T}^{2}}\frac{% 1}{(x-1)^{4}}\left[1-6x+3x^{2}+2x^{3}-6x^{2}\log(x)\right]+\mathcal{O}(\frac{1% }{m_{T}^{4}})italic_δ start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT = divide start_ARG - 1 end_ARG start_ARG 192 italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_m start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 1 end_ARG start_ARG ( italic_x - 1 ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG [ 1 - 6 italic_x + 3 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_x start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 6 italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_log ( italic_x ) ] + caligraphic_O ( divide start_ARG 1 end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) (24)

Finally, comparing the above results with the expressions for Cqsubscript𝐶𝑞C_{q}italic_C start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT and Cgsubscript𝐶𝑔C_{g}italic_C start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT given in Eqs.(3) and (4), we see that the matching relations in Eq.(23) are indeed satisfied, as expected.

Appendix C Indirect Searches - Expected Limits

All the indirect searches constraints presented in Sec.4.3 make use of the unfolded top measurements discussed in Sec.4.2. As shown in Figs.8 and 9, both the invariant mass and transverse momentum distributions display a few bins where the SM prediction is above the measured values. Although these differences are within two standard deviationsCMS-TOP-20-001 , they result in stronger limits than expected. In order to illustrate and quantify this difference we compute the expected limits on the signal using the CMS m(tt¯)𝑚𝑡¯𝑡m(t\bar{t})italic_m ( italic_t over¯ start_ARG italic_t end_ARG ) and ATLAS pT(t)subscript𝑝𝑇𝑡p_{T}(t)italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_t ) measurements following the procedure outlined in Sec.4.2, but now assuming NObsi=NSMi(NNLO)subscriptsuperscript𝑁𝑖Obssubscriptsuperscript𝑁𝑖SMNNLON^{i}_{\rm Obs}=N^{i}_{\rm SM}(\rm NNLO)italic_N start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Obs end_POSTSUBSCRIPT = italic_N start_POSTSUPERSCRIPT italic_i end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_SM end_POSTSUBSCRIPT ( roman_NNLO ). The results are shown in Figure 12, where we display the exclusion curves at 95% C.L. assuming yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10 and using the 1-loop calculation. As we can see the exclusion on mTsubscript𝑚𝑇m_{T}italic_m start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT is reduced by 200300similar-toabsent200300\sim 200-300∼ 200 - 300 GeV when we compare the observed and expected exclusions. A better treatment of the SM predictions and their uncertainties will likely bring the observed results within better agreement with the expected exclusion. Nonetheless, our overall conclusions about the potential complementarity between direct and indirect searches and the lack of validity of the EFT approximation still hold.

Refer to caption
Figure 12: Comparison of the exclusion curves from direct stop searches and from top measurements. The solid (dotted) curves show the observed (expected) exclusions. All the limits were computed using the 1-loop form factors with yDM=10subscript𝑦𝐷𝑀10y_{DM}=10italic_y start_POSTSUBSCRIPT italic_D italic_M end_POSTSUBSCRIPT = 10. The limits obtained using the CMS invariant mass measurement from Ref.CMS-TOP-20-001 are shown by the blue curves, while the ones obtained from the ATLAS pTsubscript𝑝𝑇p_{T}italic_p start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT measurementATLAS-TOPQ-2019-23 are shown in red.

References

  • (1) CMS collaboration, Search for new physics in top quark production with additional leptons in the context of effective field theory using 138fb-1 of proton-proton collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 13 TeV, Tech. Rep. CMS-PAS-TOP-22-006, CERN, Geneva (2023).
  • (2) CMS collaboration, Search for new physics using effective field theory in 13 TeV pp collision events that contain a top quark pair and a boosted Z or Higgs boson, Phys. Rev. D 108 (2023) 032008 [2208.12837].
  • (3) CMS collaboration, Constraints on anomalous Higgs boson couplings to vector bosons and fermions from the production of Higgs bosons using the τ𝜏\tauitalic_ττ𝜏\tauitalic_τ final state, Phys. Rev. D 108 (2023) 032013 [2205.05120].
  • (4) CMS collaboration, Combined Higgs boson production and decay measurements with up to 137 fb-1 of proton-proton collision data at sqrts = 13 TeV, Tech. Rep. CMS-PAS-HIG-19-005, CERN, Geneva (2020).
  • (5) CMS collaboration, Measurement of Wγ±superscript𝛾plus-or-minus{}^{\pm}\gammastart_FLOATSUPERSCRIPT ± end_FLOATSUPERSCRIPT italic_γ differential cross sections in proton-proton collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 13 TeV and effective field theory constraints, Phys. Rev. D 105 (2022) 052003 [2111.13948].
  • (6) ATLAS collaboration, Combined effective field theory interpretation of Higgs boson and weak boson production and decay with ATLAS data and electroweak precision observables, Tech. Rep. ATL-PHYS-PUB-2022-037, CERN, Geneva (2022).
  • (7) ATLAS collaboration, Measurements of W+W+1limit-fromsuperscript𝑊superscript𝑊1W^{+}W^{-}+\geq 1~{}italic_W start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_W start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT + ≥ 1jet production cross-sections in pp𝑝𝑝ppitalic_p italic_p collisions at s=13𝑠13\sqrt{s}=13~{}square-root start_ARG italic_s end_ARG = 13TeV with the ATLAS detector, JHEP 06 (2021) 003 [2103.10319].
  • (8) ATLAS collaboration, Top EFT summary plots September 2023, Tech. Rep. ATL-PHYS-PUB-2023-027, CERN, Geneva (2023).
  • (9) J. Ellis, V. Sanz and T. You, Complete Higgs Sector Constraints on Dimension-6 Operators, JHEP 07 (2014) 036 [1404.3667].
  • (10) M.B. Gavela, J.M. No, V. Sanz and J.F. de Trocóniz, Nonresonant Searches for Axionlike Particles at the LHC, Phys. Rev. Lett. 124 (2020) 051802 [1905.12953].
  • (11) M.G. Folgado and V. Sanz, On the Interpretation of Nonresonant Phenomena at Colliders, Adv. High Energy Phys. 2021 (2021) 2573471 [2005.06492].
  • (12) CMS collaboration, Search for resonant and nonresonant new phenomena in high-mass dilepton final states at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 13 TeV, JHEP 07 (2021) 208 [2103.02708].
  • (13) CMS collaboration, Search for high mass dijet resonances with a new background prediction method in proton-proton collisions at s=𝑠absent\sqrt{s}=square-root start_ARG italic_s end_ARG = 13 TeV, JHEP 05 (2020) 033 [1911.03947].
  • (14) R.S. Chivukula, P. Ittisamai, K. Mohan and E.H. Simmons, Broadening the Reach of Simplified Limits on Resonances at the LHC, Phys. Rev. D 96 (2017) 055043 [1707.01080].
  • (15) R. Cepedello, F. Esser, M. Hirsch and V. Sanz, Mapping the SMEFT to discoverable models, JHEP 09 (2022) 229 [2207.13714].
  • (16) R. Cepedello, F. Esser, M. Hirsch and V. Sanz, SMEFT goes dark: Dark Matter models for four-fermion operators, JHEP 09 (2023) 081 [2302.03485].
  • (17) A. Drozd, J. Ellis, J. Quevillon and T. You, Comparing EFT and Exact One-Loop Analyses of Non-Degenerate Stops, JHEP 06 (2015) 028 [1504.02409].
  • (18) A. Lessa and V. Sanz, Beyond Top EFT - Code and Data, 2023. 10.5281/zenodo.10277973.
  • (19) A. Buckley, C. Englert, J. Ferrando, D.J. Miller, L. Moore, M. Russell et al., Constraining top quark effective theory in the LHC Run II era, JHEP 04 (2016) 015 [1512.03360].
  • (20) S.M. Etesami, S. Khatibi and M. Mohammadi Najafabadi, Study of top quark dipole interactions in tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG production associated with two heavy gauge bosons at the LHC, Phys. Rev. D 97 (2018) 075023 [1712.07184].
  • (21) F. Maltoni, L. Mantani and K. Mimasu, Top-quark electroweak interactions at high energy, JHEP 10 (2019) 004 [1904.05637].
  • (22) I. Brivio, S. Bruggisser, F. Maltoni, R. Moutafis, T. Plehn, E. Vryonidou et al., O new physics, where art thou? A global search in the top sector, JHEP 02 (2020) 131 [1910.03606].
  • (23) V. Miralles, M.M. López, M.M. Llácer, A. Peñuelas, M. Perelló and M. Vos, The top quark electro-weak couplings after LHC Run 2, 2107.13917.
  • (24) Z. Kassabov, M. Madigan, L. Mantani, J. Moore, M. Morales Alvarado, J. Rojo et al., The top quark legacy of the LHC Run II for PDF and SMEFT analyses, JHEP 05 (2023) 205 [2303.06159].
  • (25) B. Grzadkowski, M. Iskrzynski, M. Misiak and J. Rosiek, Dimension-Six Terms in the Standard Model Lagrangian, JHEP 10 (2010) 085 [1008.4884].
  • (26) R. Alonso, E.E. Jenkins, A.V. Manohar and M. Trott, Renormalization Group Evolution of the Standard Model Dimension Six Operators III: Gauge Coupling Dependence and Phenomenology, JHEP 04 (2014) 159 [1312.2014].
  • (27) J. Ellis, M. Madigan, K. Mimasu, V. Sanz and T. You, Top, Higgs, Diboson and Electroweak Fit to the Standard Model Effective Field Theory, JHEP 04 (2021) 279 [2012.02779].
  • (28) U. Haisch and E. Re, Simplified dark matter top-quark interactions at the LHC, JHEP 06 (2015) 078 [1503.00691].
  • (29) T. Plehn, J. Thompson and S. Westhoff, Dark Matter from Electroweak Single Top Production, Phys. Rev. D 98 (2018) 015012 [1712.08065].
  • (30) A. Delgado, A. Martin and N. Raj, Forbidden Dark Matter at the Weak Scale via the Top Portal, Phys. Rev. D 95 (2017) 035002 [1608.05345].
  • (31) M. Garny, J. Heisig, M. Hufnagel and B. Lülf, Top-philic dark matter within and beyond the WIMP paradigm, Phys. Rev. D 97 (2018) 075002 [1802.00814].
  • (32) J. Fuentes-Martín, M. König, J. Pagès, A.E. Thomsen and F. Wilsch, A proof of concept for matchete: an automated tool for matching effective theories, Eur. Phys. J. C 83 (2023) 662 [2212.04510].
  • (33) CMS collaboration, Observation of four top quark production in proton-proton collisions at s=13TeV, Phys. Lett. B 847 (2023) 138290 [2305.13439].
  • (34) ATLAS collaboration, Observation of four-top-quark production in the multilepton final state with the ATLAS detector, Eur. Phys. J. C 83 (2023) 496 [2303.15061].
  • (35) W.A. Bardeen and C.T. Hill, Top quark condensates and the symmetry breaking of the electroweak interactions, Adv. Ser. Direct. High Energy Phys. 10 (1992) 649.
  • (36) T. Hahn, Generating Feynman diagrams and amplitudes with FeynArts 3, Comput. Phys. Commun. 140 (2001) 418 [hep-ph/0012260].
  • (37) V. Shtabovenko, R. Mertig and F. Orellana, FeynCalc 9.3: New features and improvements, Comput. Phys. Commun. 256 (2020) 107478 [2001.04407].
  • (38) V. Shtabovenko, R. Mertig and F. Orellana, New Developments in FeynCalc 9.0, Comput. Phys. Commun. 207 (2016) 432 [1601.01167].
  • (39) R. Mertig, M. Bohm and A. Denner, FEYN CALC: Computer algebraic calculation of Feynman amplitudes, Comput. Phys. Commun. 64 (1991) 345.
  • (40) C. Degrande, Automatic evaluation of UV and R2 terms for beyond the Standard Model Lagrangians: a proof-of-principle, Comput. Phys. Commun. 197 (2015) 239 [1406.3030].
  • (41) G. Passarino and M.J.G. Veltman, One Loop Corrections for e+ e- Annihilation Into mu+ mu- in the Weinberg Model, Nucl. Phys. B 160 (1979) 151.
  • (42) G. ’t Hooft and M.J.G. Veltman, Scalar One Loop Integrals, Nucl. Phys. B 153 (1979) 365.
  • (43) C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer and T. Reiter, UFO - The Universal FeynRules Output, Comput. Phys. Commun. 183 (2012) 1201 [1108.2040].
  • (44) J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079 [1405.0301].
  • (45) R. Frederix, S. Frixione, V. Hirschi, D. Pagani, H.S. Shao and M. Zaro, The automation of next-to-leading order electroweak calculations, JHEP 07 (2018) 185 [1804.10017].
  • (46) P. Artoisenet, R. Frederix, O. Mattelaer and R. Rietkerk, Automatic spin-entangled decays of heavy resonances in Monte Carlo simulations, JHEP 03 (2013) 015 [1212.3460].
  • (47) ATLAS collaboration, Measurements of differential cross-sections in top-quark pair events with a high transverse momentum top quark and limits on beyond the Standard Model contributions to top-quark pair production with the ATLAS detector at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 13 TeV, JHEP 06 (2022) 063 [2202.12134].
  • (48) S. Kraml, S. Kulkarni, U. Laa, A. Lessa, W. Magerl, D. Proschofsky-Spindler et al., SModelS: a tool for interpreting simplified-model results from the LHC and its application to supersymmetry, Eur. Phys. J. C 74 (2014) 2868 [1312.4175].
  • (49) G. Alguero, S. Kraml and W. Waltenberger, A SModelS interface for pyhf likelihoods, Comput. Phys. Commun. 264 (2021) 107909 [2009.01809].
  • (50) G. Alguero, J. Heisig, C.K. Khosa, S. Kraml, S. Kulkarni, A. Lessa et al., Constraining new physics with SModelS version 2, JHEP 08 (2022) 068 [2112.00769].
  • (51) M. Mahdi Altakach, S. Kraml, A. Lessa, S. Narasimha, T. Pascal and W. Waltenberger, SModelS v2.3: Enabling global likelihood analyses, SciPost Phys. 15 (2023) 185 [2306.17676].
  • (52) CMS collaboration, Search for new particles in events with energetic jets and large missing transverse momentum in proton-proton collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 13 TeV, JHEP 11 (2021) 153 [2107.13021].
  • (53) C. Bierlich et al., A comprehensive guide to the physics and usage of PYTHIA 8.3, 2203.11601.
  • (54) DELPHES 3 collaboration, DELPHES 3, A modular framework for fast simulation of a generic collider experiment, JHEP 02 (2014) 057 [1307.6346].
  • (55) W. Beenakker, C. Borschensky, R. Heger, M. Krämer, A. Kulesza and E. Laenen, NNLL resummation for stop pair-production at the LHC, JHEP 05 (2016) 153 [1601.02954].
  • (56) CMS collaboration, Combined searches for the production of supersymmetric top quark partners in proton–proton collisions at s=13TeV𝑠13TeV\sqrt{s}=13\,\text{Te}\text{V}square-root start_ARG italic_s end_ARG = 13 italic_Te italic_V, Eur. Phys. J. C 81 (2021) 970 [2107.10892].
  • (57) ATLAS collaboration, Search for a scalar partner of the top quark in the all-hadronic tt¯𝑡¯𝑡t{\bar{t}}italic_t over¯ start_ARG italic_t end_ARG plus missing transverse momentum final state at s=13𝑠13\sqrt{s}=13square-root start_ARG italic_s end_ARG = 13 TeV with the ATLAS detector, Eur. Phys. J. C 80 (2020) 737 [2004.14060].
  • (58) CMS collaboration, Search for supersymmetry in proton-proton collisions at 13 TeV in final states with jets and missing transverse momentum, JHEP 10 (2019) 244 [1908.04722].
  • (59) F. Esser, M. Madigan, V. Sanz and M. Ubiali, On the coupling of axion-like particles to the top quark, JHEP 09 (2023) 063 [2303.17634].
  • (60) CMS collaboration, Measurement of differential tt¯𝑡¯𝑡t\bar{t}italic_t over¯ start_ARG italic_t end_ARG production cross sections in the full kinematic range using lepton+jets events from proton-proton collisions at s𝑠\sqrt{s}square-root start_ARG italic_s end_ARG = 13  TeV, Phys. Rev. D 104 (2021) 092013 [2108.02803].
  • (61) S. Catani, S. Devoto, M. Grazzini, S. Kallweit and J. Mazzitelli, Top-quark pair production at the LHC: Fully differential QCD predictions at NNLO, JHEP 07 (2019) 100 [1906.06535].
  • (62) S. Catani, S. Devoto, M. Grazzini, S. Kallweit, J. Mazzitelli and H. Sargsyan, Top-quark pair hadroproduction at next-to-next-to-leading order in QCD, Phys. Rev. D 99 (2019) 051501 [1901.04005].
  • (63) M. Grazzini, S. Kallweit and M. Wiesemann, Fully differential NNLO computations with MATRIX, Eur. Phys. J. C 78 (2018) 537 [1711.06631].