Proceedings of the Fifth Annual LHCP July 27, 2019
Extra Quarks Decaying to Dark Matter
[0.25cm] Beyond the Narrow Width Approximation
Hugo Prager^{†}^{†}†Speaker., Stefano Moretti, Dermot O’Brien, Luca Panizzi
School of Physics and Astronomy, University of Southampton, Southampton SO17 1BJ, UK
[0pt] Particle Physics Department, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, UK
[0pt] Physics Department, CERN, CH1211, Geneva 23, Switzerland
[0pt] Dipartimento di Fisica, Università di Genova and INFN, Sezione di Genova,
via Dodecaneso 33, 16146 Genova, Italy
ABSTRACT
We explore the effects induced by a finite width in processes of pair production of a heavy topquark partner and its subsequent decay into a bosonic Dark Matter (DM) candidate – either scalar or vector – and a SM uptype quark at the Large Hadron Collider (LHC). We discuss the configurations of masses, widths and couplings where this phenomenology can be important in a simple model with just one such objects. Finally, we emphasise the correct definition of signal and background to be adopted as well as stress the importance of new dedicated experimental searches.
PRESENTED AT
The Fifth Annual Conference
on Large Hadron Collider Physics
Shanghai Jiao Tong University, Shanghai, China
May 1520, 2017
1 Introduction
Understanding the nature of potential eXtra Quarks (XQs) and evident DM may possibly mean to study the two sides of the same coin, as these states appear together in a variety of Beyond the SM (BSM) scenarios, e.g., in Universal Extra Dimensions [1, 2]), where indeed the DM candidate (scalar or vector) can be produced via the decay of an XQ. In this kind of scenarios, we get signatures with significant missing transverse energy (), hence similar to the case of Supersymmetry (SUSY) with conserved parity, an aspect which makes it possible to interpret SUSY results in the Narrow Width Approximation (NWA) in terms of limits on (narrow) XQs [3, 4]. The main benefit of this is of course that this procedure, based on the NWA assumption, makes the approach model independent, as the Branching Ratio (BR) of the XQ into DM (and a standard quark) is the only unknown BSM parameter, apart from the XQ mass, since the dominant production channel of XQs at the LHC is QCD induced pair production. However, this approach is rather limiting in terms of the XQ parameters that can be probed, as the XQ can have a large decay width. In fact, here we want to evaluate the effects of a large width of the XQ in the determination of its crosssection and, hence, in the (re)interpretation of bounds from experimental searches. Recall that, in these circumstances, i.e., when the XQ is offshell, in order to preserve gauge invariance, additional diagrams should be included in the definition of the signal, as we shall do below. These will include topologies with one XQ only. As immediate consequence of this, we notice that such an approach becomes necessarily model dependent. It is the purpose of this writeup to discuss the salient features of this phenomenology, borrowing results from [5].
Lagrangian terms
The interaction between a singlet DM (scalar or vector) and the XQ for a coupling with third generation SM quark is described by the following lagrangian terms^{*}^{*}*The interaction terms for a coupling with first and second generation quark are completely analogous.
where and represent XQ singlets and a XQ doublet, while and represent scalar and vector DM. All the new particles are odd under a symmetry which is needed to make the DM stable under the assumption that all the SM states are even under the same symmetry. Notice that depending on the values of the coupling constants the XQ can be either vectorlike or chiral[4].
Our results have been obtained considering the vectorlike component of a doublet, although we checked that singlet and chiral XQs produce analogous results.
Observables and conventions
Two different processes leading to the same final state are considered:

QCD pair production and onshell decay: in the NWA regime the crosssection is given by and this process is usually considered in experimental searches for XQ pair production,

Full signal: all the topologies containing at least one XQ are taken into account, including some which are missing in the NWA, and the crosssection will be labelled by ,
and we study the ratio to provide a quantitative measure of the difference between the two processes.
Analysis tools and experimental searches
We perform our study by scanning over three free parameters: the masses , and the width . We consider scenarios where the DM state has masses = 10 GeV and 500 GeV^{†}^{†}†The case of DM with mass of 1 TeV has also been considered in [5] but it will not be presented here. and where the XQ has mass , with (such that its onshell decay is kinematically allowed) up to = 2500 GeV, which is the maximal value to obtain at least 1 signal event for pair production at LHC@13TeV with 100/fb integrated luminosity. The width ranges from 0% (NWA) to 40% of the mass.
Our numerical results at partonic level are obtained using MadGraph5 [6, 7] and a model we implemented in Feynrules [8]. Events are simulated using the PDF set cteq6l1 [9] and then passed to Pythia 8 [10, 11], which takes care of the hadronisation and parton showering.
To analyse our results and perform a recasting with a set of 13 TeV analyses considering final states compatible with our scenarios, we employ CheckMATE 2 [12], which uses Delphes 3 [13] for the emulation of detector effects. We consider all the ATLAS and CMS analyses at 13 TeV available within the CheckMATE database, the most relevant ones in our case being ATLAS 1604.07773 [14], ATLAS 1605.03814 [15], ATLASCONF2016050 [16].
2 Extra quark interacting with DM and the SM top quark
The possible decay channels for a XQ coupling to third generation SM quarks are , i.e. .
Large width effects at parton level
We show in Fig.2 the ratio for an LHC energy of 13 TeV for two representative cases.
Relative difference between full signal and NWA crosssections. Left: vector DM of mass 10 GeV; right: scalar DM of mass 500 GeV. \parpic
As expected, the offshell contributions are negligible in the NWA and become more and more relevant as increases. The increase is especially fast when the mass approaches the kinematic limit, and this can be explained by a nontrivial combination of factors: as the gap between the masses decreases, a larger width opens a larger phase space for decay. Moreover, in some regions becomes similar to even for large values of the width; nevertheless, the differential distributions of the signal in the NWA and finite width are different and therefore, it is not possible to describe the process through the NWA approximation.
Large width effects at detector level
We show in Fig.2 the results obtained with the searches implemented in CheckMATE for a DM particle of mass 10 GeV and 500 GeV.
CheckMATE results for a DM particle of mass 10 GeV and 500 GeV. We show in black (grey) the exclusion line for the scalar (vector) DM scenario. \parpic[right]
The width dependence of the bounds is very weak in both cases, and the numerical values of the bounds for scalar and vector DM are similar, despite the fact that the signal crosssections are different. The weak dependence of the bounds on the width is determined by a combination of crosssection and selection efficiencies: the cuts on and have a strong effect on the signal and result in smaller efficiencies when increases.
Left: Signal crosssection; right: efficiency of the most relevant signal region. \parpic
We show in Fig.2 the signal crosssection and the efficiency for the most relevant signal region together with the exclusion line for a scalar DM of mass 10 GeV. The plot clearly shows that the increase of crosssection is compensated by the decrease of efficiency. Therefore the NWA hypothesis made for experimental searches does not overestimate the bound, even if the has a large width.
3 Extra quark interacting with DM and the SM up quark
\parpic[right]
We now study the case of XQs coupling a DM candidate and first generation SM quarks. The possible final states are and , i.e. .
In this case new topologies containing collinear divergences (due to the gluon splitting) are present, such as the one showed on the right. These diagrams affect the signal but not the QCD pair production, therefore giving rise to a large increase of the ratio . For this reason we use logarithmic plots in the following.
Large width effects at parton level
We show in Fig.3 the relative difference between full signal and NWA crosssections for representative scenarios.
Relative difference between full signal and NWA crosssections. Left: vector DM of mass 10 GeV; right: scalar DM of mass 500 GeV.
\parpic
The ratio is much larger than for third generation coupling and correctly approaches zero in the NWA. There are no regions where cancellations make similar to for large width, but the crosssection ratio still decreases in regions that are very similar to the ones observed for 3rd generation. This means that even if a cancellation of opposite effects still takes place, the effect of the additional diagrams is dominant.
Large width effects at detector level
In Fig.3 we observe a sizable width dependence of the bounds, which are still similar for scalar and vector DM in the NWA, as expected. Furthermore, it is possible to distinguish scalar and vector DM.
The shape of the bounds is again fetermined by a combination of crosssection and efficiencies contributions. The crosssection effect is dominant in this case: for both scalar and vector DM, the bounds track the different behaviours in the scaling of . In Fig.3 we have shown examples with one of the two DM spins.
The NWA hypothesis made by experimentalist may underestimate the bounds if the width is large, as the bounds may reach values larger than 2 TeV.
4 Exclusion limits in the plane
The scenarios we are considering are described by three parameters: the masses of the and of the DM and the width of the . The 2 exclusion bound identifies a 3D surface in the space defined by such parameters and therefore it is instructive to study the projections of the surface on the plane identified by the masses of and DM for different values of the ratio. This allows us to directly compare bounds on and bosonic DM with analogous results in other models, such as supersymmetry. Exclusion limits for SUSY searches of stops are often presented in the plane, where is the neutralino. Fig.1 shows the bounds in the plane for specific values of : the NWA case, 20% and 40%. Results for a quark coupling to DM and the charm quark are included in this figure.
The assumptions about the couplings of with SM quarks strongly affect the qualitative behaviours of the exclusion limits. In the NWA the bounds exhibit the same shape than the SUSY ones whereas if the has a large width the region where the is almost degenerate with the DM is also excluded, even if the masses are large, above the TeV. The dependence on the width is stronger when the couples with the SM up quark, and in this case it is also possible to distinguish scalar from vector DM in the large width regime.
In all cases the bounds in the NWA are less stringent than those in the large width regime. The qualitative difference in the width dependence of the bounds, depending on the different assumptions about the couplings, is however remarkable. The design of new experimental searches could definitely take into account these effects to try achieve not only the discovery of new signals in channels with but also a characterisation of the discovered signal.
5 Conclusions
In conclusion, results from modeldependent LHC searches for XQs decaying to DM (spin 0 or 1) and a SM quark (which we have taken here to be both of type) do not account for effects induced by either the large XQ width, additional XQ topologies (required by gauge invariance) or both, yet, these can be very large even in a simplified model with only one XQ, as we shown here. Therefore, in the attempt to make use of existing experimental results in such circumstances, one should rescale observed limits to the actual ones upon accounting for such effects (as we have done here) or else attempt deploying new ones adopting different selection strategies which minimise (in the case of exclusion) or indeed exalt (in the case of discovery) such effects. The reward of this approach, albeit being modeldependent, can be tantalising. On the one hand, more sensitivity can be gained to XQs and DM. On the other hand, the spin nature of the latter can be accessed. This is generally more possible in the case of XQ decays to light quarks (first and second generation) than to heavy ones (third generation).
References
 [1] I. Antoniadis, A Possible new dimension at a few TeV, Phys.Lett. B246 (1990) 377–384.
 [2] G. Cacciapaglia, A. Deandrea and J. LlodraPerez, A Dark Matter candidate from Lorentz Invariance in 6D, JHEP 1003 (2010) 083, [0907.4993].
 [3] G. Cacciapaglia, A. Deandrea, J. Ellis, J. Marrouche and L. Panizzi, LHC MissingTransverseEnergy Constraints on Models with Universal Extra Dimensions, Phys. Rev. D87 (2013) 075006, [1302.4750].
 [4] S. Kraml, U. Laa, L. Panizzi and H. Prager, Scalar versus fermionic top partner interpretations of searches at the LHC, JHEP 11 (2016) 107, [1607.02050].
 [5] S. Moretti, D. O’Brien, L. Panizzi and H. Prager, Production of extra quarks decaying to Dark Matter beyond the Narrow Width Approximation at the LHC, 1705.07675.
 [6] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer and T. Stelzer, MadGraph 5 : Going Beyond, JHEP 1106 (2011) 128, [1106.0522].
 [7] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer et al., The automated computation of treelevel and nexttoleading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079, [1405.0301].
 [8] A. Alloul, N. D. Christensen, C. Degrande, C. Duhr and B. Fuks, FeynRules 2.0  A complete toolbox for treelevel phenomenology, Comput. Phys. Commun. 185 (2014) 2250–2300, [1310.1921].
 [9] J. Pumplin, D. Stump, J. Huston, H. Lai, P. M. Nadolsky et al., New generation of parton distributions with uncertainties from global QCD analysis, JHEP 0207 (2002) 012, [hepph/0201195].
 [10] T. Sjostrand, S. Mrenna and P. Z. Skands, A Brief Introduction to PYTHIA 8.1, Comput. Phys. Commun. 178 (2008) 852–867, [0710.3820].
 [11] T. Sjostrand, S. Mrenna and P. Z. Skands, PYTHIA 6.4 Physics and Manual, JHEP 0605 (2006) 026, [hepph/0603175].
 [12] D. Dercks, N. Desai, J. S. Kim, K. Rolbiecki, J. Tattersall and T. Weber, CheckMATE 2: From the model to the limit, 1611.09856.
 [13] DELPHES 3 collaboration, J. de Favereau et al., DELPHES 3, A modular framework for fast simulation of a generic collider experiment, JHEP 1402 (2014) 057, [1307.6346].
 [14] ATLAS collaboration, M. Aaboud et al., Search for new phenomena in final states with an energetic jet and large missing transverse momentum in collisions at ââTeV using the ATLAS detector, Phys. Rev. D94 (2016) 032005, [1604.07773].
 [15] ATLAS collaboration, M. Aaboud et al., Search for squarks and gluinos in final states with jets and missing transverse momentum at 13 TeV with the ATLAS detector, Eur. Phys. J. C76 (2016) 392, [1605.03814].
 [16] ATLAS collaboration, “Search for top squarks in final states with one isolated lepton, jets, and missing transverse momentum in = 13 TeV pp collisions with the ATLAS detector.” http://cds.cern.ch/record/2206132, Aug, 2016.