Going beyond Top EFT
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 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 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 () much smaller than their mass (), and typically their performance quickly stops at 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 ( and 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, . 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 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,
where . 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), :
and the coupling between right-handed tops and right-handed quarks,
In Ref. Ellis:2020unq , we showed that the operator was mostly constrained by Run 1 and Run 2 Higgs observables plus the TeVatron and LHC datasets. Moreover, among the four-fermion operators, only the (= , or ) operators were constrained, predominantly by the top data and, to a lesser extent, by 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) |
---|---|---|
0.36 | ||
-0.4 | 5. | |
-0.45 | 4.0 | |
-1.0 | -0.42 |
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 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 (), singlet under , and a singlet fermion (), which is a Dark Matter candidate. Under the imposed symmetry the BSM fields are odd and the SM are even, so the renormalizable BSM lagrangian becomes:
(1) |
with , so the DM candidate is stable. The only viable decay channel for is , where the top is off-shell if . Note that the interactions between the Dark Matter candidate and the SM are fully controlled by the 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 coupling is fixed by the LSP composition and it is of the order of the EW couplings (). In the scenario discussed here we assume to be a free parameter, which can be as large as allowed by perturbativity, . 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: and as long as the value of is properly chosen. In particular, large coupling values, e.g. , 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 () 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:
(2) | |||||
where represent any light quark flavor, represents the 3rd generation quark doublet and is the (on-shell) top mass.
The connection with the SMEFT operators in the Warsaw basis described in Sec. 2 and the two operators is as follows
where is the top Yukawa. Unlike the general SMEFT framework, the and coefficients are correlated and determined by the underlying UV parameters. These coefficients were computed using Matchete Fuentes-Martin:2022jrf and are given by:
(3) | |||||
(4) | |||||
(5) |
where is the top mass and . For example, taking the limit (), one finds that this model produces a particular pattern in the SMEFT parameter space:
(6) | ||||
(7) |
We point out that for the results discussed later we only consider BSM contributions up to order and the operator will be ignored. Furthermore, the 4-top operator 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 scaling, the bounds on our benchmark would be weaker than those from final state, due to the kinematic suppression due to a top radiating three tops from and the inherent loop suppression. Note, also, that a large four-top operator is constrained by perturbativity and, in particular, by the absence of bound states. As it was first discussed in Ref. Bardeen:1992mhn , bound states would be formed when four-fermion interactions like grow above some critical value, estimated to be , where represents the scale at which confinement occurs. Comparing with the matching in Eq. 6 and approximating , non perturbativity would require . Therefore, for the results presented below, we impose .
The size of the and coefficients and their ratio is illustrated in Fig. 1. As we can see, for BSM masses around 1 TeV and , the coefficients are TeV-2. Furthermore, we see that is typically times larger than and is always negative. As discussed in Sec. 2, usual SMEFT analysis constrain these coefficients to 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 , only the sub-TeV region of parameter space can be tested.
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:
(8) |
where the and 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 and . 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.
While contains quite a large number of terms and an involved tensor structure, the expression for the form factor can be written in a compact form using the triangular Passarino-Veltmann loop functions Passarino:1978jh ; tHooft:1978jhc :
(9) |
where , are the top and anti-top momenta and . All the loop integrals are functions of and 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 () and neglecting the top mass () they simplify considerably, resulting in:
(10) |
where:
(11) | ||||
In Fig. 3 we show the form factor dependence on coming from the factor defined in Eq.( 11). As we can see, peaks around or , displaying a very broad resonant behavior. We also see that for the form factor is dominated by its imaginary part, resulting in a negative interference term, as discussed below.
3.2 Differential distributions in
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 and the top transverse momentum . 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 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: .
The distributions were computed at leading order in and (), which corresponds to the Born plus the interference terms:
(12) |
In the EFT approach this is equivalent to keeping only the terms.222We have verified that the contribution from the term () is always subdominant for perturbative values of the BSM coupling, . 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 () is only affected by , while the gluon process () depends on both form factors, it is interesting to investigate the individual contributions from each process. For instance, in the EFT regime, we have , resulting in a negative interference contribution from .
In Fig. 4 we show the (normalized) distributions for two sets of BSM masses and . 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 GeV. This is expected, since the EFT approximation is clearly not valid for . 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 . The behavior of the 1-loop distribution is also negative for most values of the invariant mass, except for , where the distribution resembles a broad resonance and its contribution becomes positive. The contribution from the quark initiated process (), 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 (relative to the SM), while the 1-loop distribution presents a very broad enhancement around TeV. In addition, for large invariant mass values ( TeV), the 1-loop 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, presents a broad bump behaviour near and a dominance of the imaginary part for large values. As a result, the full 1-loop distribution displays an excess for the invariant mass bins close to 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 . 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 , as expected. The broad resonant behavior of the 1-loop distributions is once again present, but now it only starts to appear at 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.
In addition to the top pair invariant mass, the top 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 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 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 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 . It is also important to point out that while the EFT approximation is valid for , it is not so easy to identify a similar condition for the distribution. This can be seen in Fig.5 (right), where the EFT approximation fails even at values much smaller than .
Boosted Top Phase-Space | |
---|---|
Jet Cuts | |
GeV | |
Fat jet Cuts | |
GeV | |
120 GeV GeV | |
contains one | |
Lepton Cuts | |
GeV | |
GeV | |
Cut | GeV |
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 production, and ii) indirect searches, i.e. measurements of 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: or . For instance, for sufficiently large masses we expect the EFT results to be valid. Also, for very small , the loop contributions to top pair production are suppressed and direct searches will be more sensitive.
4.1 Direct Searches
If is sufficiently light ( TeV), it will be copiously produced at the LHC, since it is colored. The leading order diagrams for production are shown in Figure 6. The signatures generated by production and decay strongly depend on the mass difference . For sufficiently large , the signature corresponds to on-shell tops plus missing energy (), while in the compressed scenario () we have plus -jets and additional (soft) leptons and light jets. The compressed scenario tends to be more challenging, resulting in weaker constraints on . We also point out that the signal does not dependent on the BSM coupling , except for the scalar width. Although for sufficiently small widths the scalar can become long-lived, in the following we assume large enough so 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 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 can become relevant. Therefore, in addition to the stop searches, we have considered the CMS jets plus 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 production cross-section at leading order (LO), a constant k-factor was used to approximate the NNLO+NNLL result Beenakker:2016gmf .
In Figure 7 we show the 95% C.L. excluded region in the vs plane. For large mass differences (), 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 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 GeV are excluded, while for large the exclusion goes up to 1.3 TeV.
4.2 Indirect searches in
As discussed in Sec. 4.1, the constraints from direct searches for depend strongly on the scalar-DM mass difference and exclude masses up to TeV. These constraints, however, do not depend on the BSM coupling (), except for , which could render the colored scalar long-lived. On the other hand, if , BSM loop contributions to 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 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:
(13) |
where () is the background (signal) prediction computed at LO using MadGraph5_aMC@NLO. The reweighting factors obtained through this procedure are typically for the CMS invariant mass bins and for the ATLAS transverse momentum bins. With the above expressions and the covariance matrices () provided by the experimental collaborations, we can then make use of the measured distributions () to compute limits on the BSM coupling . Following the same procedure used by the experimental collaborations we define a function as:
(14) |
where , so the 95% C.L. limit on corresponds to .
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
We first consider the CMS measurement CMS-TOP-20-001 of the differential cross-sections at TeV using the full Run 2 luminosity, 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 GeV and GeV and correspond to a compressed scenario currently beyond the reach of direct searches. The BSM coupling is considerably large, , but still within the perturbative regime. As we can see, the 1-loop distribution significantly deviates from data for 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 at 95% C.L. using the 1-loop calculation, while the EFT approximation results in . Note that the error band shown in the bottom subplot of Fig. 8 corresponds only to . The full correlations, however, are essential for computing the limits and provide stronger constraints than assuming uncorrelated bin uncertainties.
ATLAS
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 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 , 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 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 and "light" BSM masses, GeV, GeV, while the second one shows the same distributions, but for heavier masses ( TeV and TeV) and a larger BSM coupling, . 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 , 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 coupling. In particular, for the signal shown in the left plot of Fig. 9, we obtain at 95% C.L. using the 1-loop calculation, while 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.
Note that, when compared to the invariant mass distributions from Fig. 8, the measurement seems to be more sensitive to the BSM signal than the invariant mass distribution, since a larger excess is seen in the intermediate 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, , and for large BSM couplings, . In order to compare the constraints from direct and indirect searches, we scan over the BSM masses and compute the limits on obtained from the CMS measurement and the ATLAS distribution. In Figure 10 we show the region excluded by direct searches for and the exclusion curves from indirect searches for and . The regions to the left of the blue (red) curves are excluded at 95% C.L. by the CMS invariant mass (ATLAS ) measurement. We also show the corresponding curves obtained assuming the EFT approximation. As expected from the results in the previous Section, the measurement is more sensitive to the BSM signal and excludes masses up to GeV in the highly compressed region, if we take . 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, , the exclusion curves obtained using the invariant mass distribution fall inside the direct search excluded region and are not shown. Once we consider , the measurement becomes competitive with direct searches in the compressed region and exclude masses up to GeV. But the limits obtained from the distribution are still stronger, excluding up to 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 , 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 would likely be excluded by DM Direct Detection searches Garny:2018icg . In this case, the 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.
All the results presented in Fig.10 rely on 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 TeV-2, translating to 400 GeV for , while the corresponding constraint we have obtained using the EFT approximation is 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 and four-fermion 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 production at the LHC at leading order in the BSM coupling . 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 . The counter-terms diagrams needed for regularizing the divergent diagrams are not shown.
Appendix B Form Factors and EFT matching
In order to verify the expression for the form factor defined in Eq.(9) and the EFT coefficients defined in Eqs.(4) and (4), we compute the amplitude for the quark initiated process, , 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, , we should obtain the same result for both approaches, thus validating our implementation. A similar exercise could in principle be done for the process. However, due to the large number of diagrams (see Figure 11) and the complexity of the form factor, the validation of has only been done numerically.
We start with the 1-loop calculation using form factors. Note that for the process only the top-top-gluon form factor () is relevant and the top momenta appearing in are on-shell, . In this case the form factor defined in Eq.(9) can be simplified, since and similarly . Using these results we obtain:
(15) |
where are the triangular scalar loop functions and 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 . Similar expressions can be found for the compressed region and have been used in our results for scenarios with .:
(16) |
where and . Both the divergent () and the terms in cancel with the corresponding terms from , which is the only divergent loop integral in Eq.(15).
Using the results above, the amplitude for the first diagram in Figure 11 becomes:
(17) |
where is the quark (top) spinor and is the anti-quark (anti-top) spinor. The color indices for the Gell-Mann matrices are contracted with the spinors and are not explicitly shown for simplicity.
The amplitude can be greatly simplified using the on-shell relations:
(18) |
and . Applying these relations to Eq.(17) we obtain:
(19) |
Finally, using momentum conservation, , where are the initial state momenta and the on-shell relations for the massless incoming quarks we have: . Hence:
(20) |
where we have used .
The above result corresponds to the full 1-loop calculation for the 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:
(21) |
where 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:
(22) |
In order to verify the above relations we must expand the loop functions and the counter-terms to leading order in , which gives:
(24) |
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 and ATLAS measurements following the procedure outlined in Sec.4.2, but now assuming . The results are shown in Figure 12, where we display the exclusion curves at 95% C.L. assuming and using the 1-loop calculation. As we can see the exclusion on is reduced by 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.
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 = 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 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 differential cross sections in proton-proton collisions at = 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 jet production cross-sections in collisions at TeV 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 = 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 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 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 = 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 = 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 , 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 plus missing transverse momentum final state at 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 production cross sections in the full kinematic range using lepton+jets events from proton-proton collisions at = 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].