Long-lived colored scalars at the LHC

We study the collider signatures of a long-lived massive colored scalar transforming trivially under the weak interaction and decaying within the inner sections of a detector such as ATLAS or CMS. In our study, we assume that the colored scalar couples at tree-level to a top quark and a stable fermion, possibly arising from a dark sector or from supersymmetric extensions of the Standard Model. After implementing the latest experimental searches for long-lived colored scalars, we observe a region of parameter space consistent with a colored electroweak-singlet scalar with mass between ∼200\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sim }200$$\end{document}–350 GeV and a lifetime between 0.1–1 mm/c\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {mm}/c$$\end{document} together, with a nearly degenerate dark fermion that may be probed at the s=13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=13$$\end{document} TeV LHC. We show that a search strategy using a combination of cuts on missing transverse energy and impact parameters can exclude regions of parameter space not accessed by prompt searches. We show that a region of parameter space within our simplified model may naturally arise from the light-stop window regime of supersymmetric extensions of the Standard Model, where a light mostly right-handed stop has a mass slightly larger than the lightest neutralino and decays through a four-body process.


Motivation
Searching for long-lived massive particles is an active program at the LHC. The ATLAS and CMS collaborations have carried out an extensive list of analysis targeting different topologies for the decay products of long-lived particles within the inner detector as well as particles with large enough lifetimes to escape the detector. In particular, both collaborations have carried out an analysis sensitive to very long-lived charged particles using time-of-flight (TOF) and energy loss, dE/dx information [1,2]. This class of searches is highly model independent since a main requirement for a e-mail: apuente@physics.carleton.ca b e-mail: szynkman@fisica.unlp.edu.ar an event to pass the experimental cuts is the appearance of a muon-like object in the muon spectrometer matched to a track in the inner detector. These searches rule out very longlived gluinos, stops, sbottoms, and staus with masses below ∼1250, 900, 850, 500 GeV, respectively. Searches that target signatures in the inner detector have also been dealt with at the LHC. In particular, both CMS and ATLAS have conducted an analysis looking for disappearing tracks [3,4]. This search is particularly interesting since it targets supersymmetric (SUSY) models where a chargino is near mass-degenerate with the lightest neutralino, rendering the chargino with a long lifetime. Within this scenario, the chargino decays to a neutralino and a charged pion with very low p T , low enough for its track to pass reconstruction algorithms. The limits on the chargino mass are on masses below ∼270 GeV and proper lifetimes of ∼0.2 ns. Other searches look for interesting topologies that arise within SUSY scenarios with Rparity violation. In [5,6], the CMS and ATLAS collaborations search for long-lived neutral particles that decay to quarkantiquark pairs and a muon. These long-lived exotic particles are pair produced via the decay of a non-SM Higgs boson or the decay of a scalar quark to the lightest neutralino and a light quark. In the latter case the neutralino is long-lived and decays to a muon and a quark-antiquark pair through a small R-parity-violating coupling. This search is sensitive to long-lived particles with a proper decay length between 2-40 cm and production cross sections below 0.5−3 fb. A similar search is carried out by CMS [7], and targets scalar quark decays to a long-lived neutralino that subsequently decays to a same-flavor lepton pair and a neutrino with proper decay lengths in the range 0.01-100 mm. The analysis is able to set an upper limit on the production cross section in the range 0.2-5 fb for scalar quarks masses above 350 GeV. In addition, the CMS collaboration has a very rigid search for displaced supersymmetry with a dilepton final state [8], in this case one muon and one electron with opposite charges. The search is sensitive to the decay of a long-lived scalar top quark to a lepton and a b-jet in models of R-parity violation. The search has a highest mass exclusion of scalar top quark with mass below 790 GeV and a proper lifetime of 2 cm/c. The ATLAS collaboration has performed a very comprehensive study of displaced vertices, targeting long-lived particles decaying into two leptons or five or more charged particles [9]. Limits are set on models of R-parity violation, split SUSY, and gauge mediated SUSY breaking. In addition, searches have been conducted to analyze scenarios where a long-lived particle is stopped in the calorimeters and decays out of time [10,11]. This class of analyses are sensitive to colored particles that hadronize and interact with matter, and strongly depend on the hadronization models implemented as well as the modeling of the hadron's interactions with matter. Displaced vertices are well motivated scenarios of models beyond the Standard Model. They appear within the context of SUSY breaking [12][13][14], models of weak-scale R-parity violation [15][16][17], Hidden valleys [18][19][20][21], baryogengesis [22,23], dark QCD [24] and late decays of righthanded neutrinos within SUSY extensions of the Standard Model [25]. In addition, there is an increasing activity in the recasting front, to better constrain a wide range of these models using the existing LHC displaced vertex analyses [23,24,26].
The analyses described above, which target long-lived massive particles in the inner detector, do not implement hard cuts on the missing transverse energy, / E T , to identify the signal and discriminate the background. This is due to the fact that either, one has a very compressed spectrum where there is very little missing energy carried away by a stable particle, or the decay products of the long-lived massive particles can all be reconstructed within the tracker. In addition, they are mostly driven by very specific models. In this work we would like to address whether a long-lived charged massive particle can leave a signature that can be reconstructed within the inner detector in association with a large amount of missing transverse energy, / E T . Furthermore, we wish to address this question in a model independent manner that can then be used as a template to compare experimental analyses to different theoretical frameworks. Due to the large number of experimental signatures and classes of long-lived massive particles, we focus only in the case where a colored electroweaksinglet scalar decays to a right-handed top quark and a dark fermion that is stable throughout the detector length. The dark fermion is a true source of missing energy, and depending on the parameters, it can carry away a large amount of energy. In our analysis we focus on pair production of the colored electroweak-singlet scalar and the semileptonic decay mode of the top quark. Our experimental signature would consist of a lepton pair with opposite charges, jets, and missing energy; and in this way, our experimental signature is similar to the analysis on displaced SUSY by the CMS collaboration [8], but we allow for the possibility of a large amount of / E T .
Light colored scalar particles with the quantum numbers of an up-type quark are well motivated scenarios. This is the case within SUSY extensions of the Standard Model such as the Minimal Supersymmetric Standard Model (MSSM). Within this theoretical framework scalar top quarks are used to address the electroweak hierarchy problem of the SM. In addition, a light scalar top quark can be used to explain the observed baryon asymmetry, in particular, if its mass lies below that of the top quark. Recently, light stops have been studied within the framework of various SUSY breaking mechanisms and flavor structures [27,28], and various analyses have been dedicated to the feasibility of electroweak baryogenesis (EWB) [29][30][31][32]. In addition, various phenomenological studies have looked at the possible collider signals of a light stop at the LHC. The work by the authors in [33] addresses all possible stop decay modes with special emphasis on the four-body and flavor-violating channels, while the work by the authors in [34] addresses primarily the four-body decay mode. Recently the authors in [35] have proposed a monojet-like search optimized for the four-body decay mode of the stop with an additional bjet requirement. These studies significantly complement the search efforts by the CMS [36,37] and ATLAS [38,39] collaborations to address the light-stop window regime with the full data set at center of mass energies of 8 TeV. In our study we aim to further complement the analysis of the light-stop window, that is, with slightly displaced vertices in the range 0.1-10 mm.
The summary of our study is as follows: In Sect. 2 we introduce the model and in Sect. 3 we study experimental constraints from direct searches at colliders, both prompt and long-lived. In Sect. 4 we present a search strategy that can be used to set limits on our framework at the LHC with 13 TeV center of mass energies. In Sect. 5 we discuss a model that may lead to our experimental signature and in Sect. 6 we provide concluding remarks.

Toy model
In this section we introduce a toy model aimed at parametrizing the decay of a long-lived colored electroweak-singlet scalar, φ, to a top quark and a new hidden Majorana fermion, X . The Majorana fermion communicates with the SM via the colored scalar. This new Majorana fermion is stable at scales larger than those of a modern detector and will be a source of true missing energy. We can write the coupling of the hidden sector to the SM top quark via the following Lagrangian: Models featuring a Majorana fermion, singlet under the SM gauge group and a non-trivially charged scalar have been studied in [40] in the context of minimal decaying dark matter. In particular, they have explored the possibility that the stability of the dark matter candidate, X , is not guaranteed. The study covers a wide range of parameter space in accordance with a freeze-in mechanism [41], arising from the decays of the colored scalars in thermal equilibrium, and a superWIMP mechanism [42], from the decay of the colored scalar after it has undergone freeze out, to account for the observed dark matter relic abundance. In our study we limit our analysis to a stable Majorana fermion, stable throughout the length of a modern particle collider detector such as ATLAS or CMS and look further into its identity in Sect. 5. In the limit where the mass of the colored electroweaksinglet scalar, m φ , is larger than m t + m X , the leading decay mode is given by whereas for masses below m t + m X , the three-body decay mode through an off-shell top and the four-body mode through and off shell W open up. An analytic expression for the three-body decay mode in the context of the MSSM where a stop decays to a W gauge boson, a b-quark and a neutralino has been studied in [43]. The four-body decay width can be calculated using the approximate formula introduced in [34] and its validity is for small M = m φ − m X . In particular, for masses m φ − m X < m W the four-body decay mode branching fractions are given by Furthermore, in this class of theories, a flavor-violating decay mode may exist and compete with the four-body decay channel. SUSY theories are one classic example where the colored electroweak-singlet scalar identified with a scalar top quark decays into a charm quark and the lightest neutralino. This coupling can be parametrized by a mixing angle, θ t,c , which depends on the flavor structure of the model. However, the authors in [33,34] show that this region is very complex as the four-body and flavor-violating channels dominate in different regions but become comparable for a mixing parameter of order 10 −5 or when the stop is mostly right-handed. The study in [33] carries out an in-depth analysis of the sensitivity that the LHC has to the different decay modes of a light stop within the minimal SUSY Standard Model and within the context of a promptly decaying stop. The authors in [34] instead, focus primarily on a prompt mostly right-handed stop decaying to leptons through a four-body decay channel.
In both cases it appears that a light stop is only viable if its mass lies above 200 GeV. Within our framework, two colored electroweak-singlets can be pair produced, each decaying to a b-quark, a lepton and missing energy. For small couplings, λ, and mass difference M = m φ −m X , φ may be long-lived and decay within the inner detector, particularly for proper decay lengths in the range 0.1-100 mm. The leptons can be reconstructed and the long lifetime of φ may be inferred by the displacement of the leptons from the primary vertex through their impact parameters. This is due to the large correlation that exists between the unstable particle's lifetime and the impact parameters of the decay products. The analysis by the CMS collaboration [8] is especially sensitive to such a signature with the exception that the analysis is performed without the use of hard cuts on the hadronic activity or missing transverse energy. In the following we study the sensitivity that the LHC has to our simplified model, in particular, to scenarios with large / E T .

Prompt searches at LHC 8 TeV
Collider signatures from long-lived particles are relatively free from SM backgrounds. However, conventional searches are limited to impact parameters below O(1) mm. This setup is mainly motivated to probe decays from heavy flavor quarks and τ -leptons. In addition, most searches incorporate track quality cuts with requirements on the impact parameters of charged tracks [44,45] and use impact parameter-base tagging [46][47][48][49] for heavy flavor quarks and τ 's. The latter poses a problem for long-lived particles decaying to b-jets, since this scenario will lead to displaced vertices from B decays with impact parameters above ∼1 mm; thus leading to a large suppression on the b-tagging efficiency. In order to study the properties of jets as a function of the average lifetime of a decaying colored scalar we use Delphes [50] for jet reconstruction and tracking information. We simulate the detector response for events for a slightly compressed scenario of the model discussed in the previous section using m φ = 150 GeV and m X = 20 GeV. In Fig. 1 we show in green the efficiency that all jets in the events contain at least one well reconstructed charged track as a function of the particle's lifetime. We use the ATLAS guidelines for tracks [51] that place standard quality cuts on the transverse impact parameter, |d 0 | < 2 mm, and the difference between the longitudinal impact parameter and the position of the primary vertex, |z 0 − z v | < 10 mm. We observe a decrease in the efficiency for decay lengths above 1 mm and an increase above 1 m. Above 10 mm, a jet will still contain a number of tracks but many will have impact parameters not consistent with the quality criteria. At and above 1 m, the long-lived particle decays within the calorimeter. Most tracks found inside this jet will mostly come from activity in the inner detector; these jets will have a higher probability of satisfying the quality criteria. The b-tagging efficiency is depicted by the blue line within Fig. 1 while the two b-tagging efficiency is depicted by the red line. In Delphes, the default b-tagging efficiency is of 40 % for b-jets with p T > 20 GeV and |η| < 2.5. We are using a compressed spectrum, where the b-tagging efficiency for b's originating from a prompt decay is of 20 %. In addition, Delphes does not implement any quality criteria on the impact parameters on the tracks that arise from the decay of a B-hadron. We apply the quality criteria in [52,53] and observe that for b's from prompt decays the efficiency is on the order of 10 % with a slight decrease on the b-tagging efficiency above 1 mm. The loss in efficiency is even greater for the two b-tagging efficiency. Therefore, to probe long-lived particle decays, especially within the inner detector, analyses must incorporate dedicated triggers and beyond the accustomed particle reconstruction algorithms. The latter should have the potential to pin down the origin of many particles (tracks) to one single displaced vertex [5,6,9].
Our simplified model has new colored degrees of freedom that can be pair produced. The colored scalars decay to a dark fermion that is a source of missing transverse energy and to a top quark, and depending on the spectrum and coupling strength, two-, three-, and four-body decay modes may be open. The latter is important in the region where the scalar and dark fermion have a small mass gap and the decay leads to soft jets or leptons that may escape detection since they may not be fully reconstructed. Therefore, most of the SUSY searches with zero leptons will have a potential reach to long-lived colored scalars that decay with lifetimes above 0.1 mm/c. For this class of particles the decay products, jets, will have impact parameters within the thresholds used in prompt SUSY searches. However, a significant fraction of leptons will not be identified, since tracks with impact parameters above ∼1 mm are not considered in conventional collider searches. This can be seen qualitatively in Fig. 2, where we show the fraction of events as a function of the electron's impact parameter for decay lengths in the range of 0.1-100 mm.
In order to analyze the reach of prompt searches at the 8 TeV LHC to our simplified scenario, we use the analyses validated for phenomenological use by the CheckMATE collaboration [54]. To achieve the largest sensitivity, we analyze particle lifetimes in the range 0.1-10 mm/c using the following collider searches by the ATLAS and CMS collaborations: In Fig. 3 we show the region of parameter space in the m φ -m X plane excluded with the ATLAS and CMS analyses described above for a colored scalar with a 0.1 mm decay length. In particular, we observe that the CMS analysis using the α T variable [60] yields the largest exclusion in the nondegenerate region while the monojet [39] analysis constrains the degenerate window, where most of the signal is due to a highp T jet from initial state radiation. The monojet con- Excluded regions in the m φ -m X plane with the prompt searches included in this sections at 8 TeV center of mass energies. We assume the exclusion applies to lifetimes between 0.1 and 10 mm/c since no additional suppression on the signal acceptance is expected from the displaced jets in the event straint is very insensitive to the particle's lifetime. However, the region excluded by searches incorporating large number of jets can be suppressed as the particle's lifetime increases.
Recently, the ATLAS collaboration has re-cast a number of prompt searches to constrain a long-lived gluino [10]. In particular, these searches require large number of jets and no leptons. For example, the ATLAS search for squarks and gluinos [61] requires a leading jet with p T above 130 GeV in all of its signal regions, and a second jet with p T > 60 GeV. Many ATLAS analyses reject events containing a jet with p T > 100 GeV and |η| < 2 if the charged p T fraction satisfies, f ch < 0.02 or both f ch < 0.05 and an electromagnetic fraction, f em , above 0.9 to reject cosmics or detector malfunctions [62,63]. In addition, jets fail reconstruction if they appear to originate from additional collisions, that is, jets with p T > 50 GeV and |η| < 2.5 are required to have at least one well reconstructed track as described above. In Fig. 4 we show the acceptance rate obtained after demanding between one and two jets with p T > 100 GeV and satisfying the above quality criteria. The acceptance rate is flat for lifetimes below ∼100 mm/c. If the highest efficiency (for small lifetimes) corresponds to the efficiency of current prompt searches, then prompt searches will lose sensitivity for lifetimes above 100 mm/c. In particular, for decay lengths between 100 mm and 10 m, sensitivity of prompt searches would be suppressed by roughly 10-60 %, respectively. This region of parameter space is then very difficult to probe with prompt searches as well as displaced vertex searches since in the latter, a full reconstruction of the secondary vertex is required, which can only happen at displacements from the interaction point below 600 mm. Above ∼1-10 m, the longlived particle decays within the calorimeter and most tracks found within the jet will most likely come from activity in the inner detector. This is depicted by the flat efficiency found above 1-10 m. We conclude that in our particular simplified scenario, prompt searches with lepton vetoes are incredibly effective for lifetimes below 10 mm/c. In light of this result, which falls short of an in-depth analysis of a collider detector, we encourage the collaborations to continue their studies regarding the applicability of prompt searches to lifetimes above 0.1 mm/c in order to better pinpoint the region where more specialized triggers must be used; to both address interesting physics and the issue of pile-up at a stronger hadron collider.
In the following section we analyze a long-lived search by the CMS collaboration that can further probe the parameter space of this model, in particular, regions where displaced leptons from the decays of colored scalar decays appear with impact parameters between 0.01 and 2 cm.

Long-lived searches at LHC 8 TeV
In this section we focus on an analysis by the CMS collaboration that searches for long-lived charged massive particles through their decays into leptons with 19.5 fb −1 of data at 8 TeV center of mass energies [8]. In particular, this search is sensitive to pair produced long-lived stops that decay to a bquark and a lepton through an R-parity-violating vertex [17]. The search is designed to be as model independent as possible and focuses exclusively on a displaced isolated dilepton signature. In particular, the final state consists of exactly one electron and one muon, oppositely charged, and isolated within a cone R = 0.3, 04, respectively. The leptons are required to have a p T > 25 GeV and lie within the detector acceptance of |η| < 2.5. The two leptons are required to be separated in the η-φ plane by more than 0.5. Events are rejected if a jet with p T > 10 lies within R < 0.5 of either selected lepton.
The main background leading to a displaced lepton in association with jets is from Z production with decay into a τ + τ − final state. QCD multijet and tt production are also dominant sources of background since a displaced lepton may appear from the misidentification of a jet and from the decay of a heavy flavor jet. B and D meson decays lead to displaced leptons with a mean displacement of ∼500 µm while leptons from τ decays have a mean lifetime of ∼87 µm. Other sources of background come from single top production, diboson, and W + jets production. In order to validate our analysis, we simulate the Z j (Z → τ + τ − ) and the tt backgrounds, and use the data driven QCD multijet background results quoted by the CMS collaboration. A Monte Carlo approach for the QCD component will require a very large sample to properly handle the multijet background. We use MadGraph 5 [64] in order to simulate the backgrounds and the pair production of a colored scalar with model files generated with FeynRules [65]. The parton showering and hadronization are carried out by Pythia [66], as well as the decays of the colored scalar which are implemented by a decay table that takes into account the two-, three-, and fourbody decay branching fractions. In our analysis we use the Delphes [50] fast detector simulator to reconstruct jets with the anti-k T algorithm using a parameter size of R = 0.5, and to isolate and reconstruct tracks. The lepton identification is performed by matching a track to a truth-level lepton. Furthermore, in the simulation of the backgrounds we further match the leptons to the decay of a τ or a heavy flavor jet. The electron reconstruction efficiencies are implemented using data from the CMS collaboration on electron performance [67] and we implement a muon reconstruction efficiency of 95 %. The isolation, I l , criterion is calculated by summing the p T of all tracks within a cone of size R = 0.3 and 0.4 around the electron and muon, respectively. We require that Events passing the preselection criteria described above are further categorized by their transverse impact parameter. The transverse impact parameter is calculated using a track's coordinate of closest approach to the primary vertex in the transverse plane, x d , y d , and the track's transverse momentum: In order to exclude contamination from promptly decaying particles, both leptons are required to have an impact parameter, |d 0 | above 100 µm. This requirement eliminates almost completely subdominant sources of backgrounds that yield an electron or a muon from the decay of a W , Z , and a top quark. Lastly, the analysis is constrained within a region where |d 0 | < 2. cm. In order to properly validate our analysis and to accurately constrain the parameter space of the model introduced in this work we analyze the dependence of the reconstruction and isolation efficiencies as a function of the impact parameter in a pure sample of displaced leptons arising from the decay of a colored scalar with mass m φ = 150 GeV. The results are shown in Fig. 5. The blue dashed line corresponds to a long-lived decaying particle with a mean lifetime of τ = 0.0001 mm, while the red and green dashed lines correspond to mean lifetimes of 0.1 and 100 mm/c, respectively. The plot shows a clear relationship between a particle's lifetime and the impact parameter of the decay product and it is consistent with what we observe in Fig. 2. We observe in Fig. 5 that even though leptons from decaying particles with lifetimes of 100 mm/c are likely to appear with very small impact parameters, the efficiency for proper isolation and reconstruction is smaller than leptons arising from decays of shorter lived particles at these small values of the impact parameters.
We simulate the signal pertaining to the model described in the previous section for three colored electroweak-singlet scalar decay lengths: 0.1, 1, and 10 mm in the range of masses m φ < 350 GeV and m X < 300 GeV as well as the Z j (Z → τ + τ − ) and tt backgrounds applying the preliminary cuts described above but only demanding that there is either one In order to constrain the parameter space of our simplified model we implement the CMS analysis which is per- formed in three signal regions parametrized by the transverse impact parameter: The most exclusive region requires that both the electron and the muon have an impact parameter above 0.1 cm. This region is labeled S R 3 . The intermediate region, S R 2 is populated with events not in S R 3 but with both leptons above 0.05 cm. S R 1 is populated with events failing the S R 2 criteria but with both leptons above 0.02 cm. We carry out a channel by channel exclusion using a 95 % CL number of excluded signal events. In Table 1 we show the number of observed and expected events calculated in the CMS analysis as well as the number of events to make an exclusion at 95 % confidence level. We fully recast the Z → τ τ background in the window of impact parameters 0.01 < |d 0 | < 2 cm after implementing all preliminary cuts and obtain 82 events vs. the 98 obtained by CMS. However, to make the exclusion we make use of the experimental results due to our inability to fully recast the tt and QCD backgrounds. In addition, we simulate the RPV decay mode Table 1 Number of expected and observed events given by the CMS analysis on displaced SUSY [8]. The third row corresponds to the number of excluded signal events at 95 % CL s . In addition, we have re-cast the model analyzed in the CMS search and show the number of events within each signal region for the three benchmark points in our simulation and in parentheses those of CMS t → b l analyzed by the CMS collaboration and show that our results are consistent with the experiment. In Fig. 8a (top panel) we show our exclusion region in the m φ -m X plane using the CMS results given in Table 1 for a colored scalar with a proper lifetime of 0.1 mm/c, while Fig. 8b, c correspond to lifetimes of 1 and 10 mm/c, respectively. We observe stronger exclusions in regions where the SM background is low; that is, impact parameters between 0.1-2 cm and decay lengths above 1 mm. In particular, for lifetimes above 10 mm, the CMS search for displaced leptons will reach beyond the prompt SUSY exclusions. Run II at the LHC with center of mass energies of 13 TeV will extend the reach in areas that require a larger cross section or luminosity. Regions where a larger amount of missing transverse energy is present will require an analysis that incorporates signal regions using a hard cut on / E T to reduce the SM background. This is particularly true for colored scalar masses between ∼170-250 GeV and a lifetime below 1 mm/c as well as dark fermion masses above ∼45-180 GeV.
The authors in [26] have carried out an in-depth study to recast existing long-lived searches by the CMS and ATLAS collaborations for various general scenarios of SUSY with hadronic decays. In particular, they have explored the region with light scalar top quarks, motivated by naturalness, decaying to a top quark and a gravitino in the context of gauge mediation of SUSY breaking. In this scenario, a light stop decays to a gravitino a W gauge boson and a b-quark through an off-shell top. They find that for a nearly massless gravitino, the displaced dijet analysis by CMS [5] and the displaced muon plus tracks analysis by ATLAS [6] can exclude colored scalars with masses between 100 and 200 GeV with lifetimes all the way down to 0.1 mm. However, for larger dark fermion masses, the gravitino in their context, the fourbody decay mode of the colored scalar mode dominates leading to softer jets and a lower value of H T . In addition, the analysis in [26] places a very conservative bound on the lifetime of the colored scalar of 1 mm when the decay is prompt, but as it was shown in the previous section, a lifetime below 1 mm is not completely ruled out, in particular if the colored scalar has a mass above ∼200 GeV and the dark fermion lies above ∼80 GeV.
In the following, we analyze the LHC reach to our simplified model using √ s = 13 TeV. We propose a search strategy more sensitive to colored scalars with lifetimes below 1 mm/c.

Long-lived colored scalars at LHC 13 TeV
In this section we discuss a very simple search strategy that can be used to probe long-lived colored scalars with a lifetime below 1 mm/c at the LHC. As we have shown in the previous section, for lifetimes above 1 mm/c, the CMS search [8] rules out most of the parameter space up to colored scalar masses of 350 GeV. This is due to the very little SM background present there. Even though a lifetime below ∼1. mm/c is usually considered prompt, prompt searches have not yet ruled out the entire parameter space and long-lived searches are still weak due to the large SM background. We thus outline how a search for leptons and missing energy can be used to complement the exclusion carried out by prompt searches if we further apply a cut on the lepton's impact parameter. In addition, the strategy is simple as it can be used to completely eliminate the tt and QCD multijet backgrounds. We do this through a combination of both an upper bound on the sum of the jet p T in the event and a lower bound on the missing energy. The latter is used to suppress the QCD multijet background, as most of the MET arises from the misreconstruction of jets. In addition, we only trigger on a displaced muon arising form the decay of a long-lived particle within the inner detector. Other sources of backgrounds, such as diboson and W +jets are easily suppressed after requiring a transverse impact parameter for the muon above 0.01 cm.
We use MadGraph 5 [64] in order to simulate the backgrounds and the pair production of a colored scalar with model files generated with FeynRules [65]. Furthermore, as in the analysis described in Sect. 3.2, we only simulate the Z j (Z → τ + τ − ) and tt backgrounds. The parton showering and hadronization are carried out by Pythia [66], as well as the decays of the colored scalar which are implemented by a decay table that takes into account the two-, three-, and four-body decay branching fractions. In our analysis we use the Delphes [50] fast detector simulator to reconstruct jets with p T > 10 GeV using the anti-k T algorithm with parameter R = 0.5, as well as to isolate and reconstruct tracks. The muon identification is performed by matching a track to a truth-level muon with p T > 25 GeV and |η| < 2.5. We implement a muon reconstruction efficiency of 95 % and introduce an isolation criterion, I l , by summing the p T of all tracks within a cone of size R = 0.4 around the muon and requiring that Lastly, we reject events with H T > 100 GeV in order to suppress the QCD and tt backgrounds [68], but advise that one may define different signal regions in H T to enhance the  Table 2 Signal and background cross sections after a series of cuts for a colored scalar with mass, m φ = 200 GeV and a dark fermion of mass, m X = 150 GeV. The average lifetime for the scalar is of 0.1 mm/c. In the last column we show the significance for luminosities of 10, 30, and 100 fb −1 with an overall upper bound on the sum of the jet p T given by reach, especially at collider energies above 8 TeV. The analysis is restricted to a muon's transverse impact parameters, defined in Eq. (5), between 0.01-0.1 cm to further reduce backgrounds where the lepton originates from the decay of a heavy flavor quark. In Fig. 9a, b we show the / E T and |d 0 | μ distributions for the Z j (Z → τ + τ − ) and tt backgrounds together with the signal from a colored scalar with a proper lifetime of 0.1 mm/c. In the figure, the black solid line corresponds to m φ = 180 GeV and m X = 121 GeV. The dashed back line and the solid yellow line correspond to values of m φ (m X ) = 180 (141), 180 (161) GeV, respectively.
In Table 2 we show the LHC reach at 13 TeV for a colored scalar with mass m φ = 200 GeV and a lifetime of 0.1 mm/c together with a dark fermion with mass, m X = 150 GeV. To probe the sensitivity to this simplified model we define the significance of a signal using the following statistical estimator: where N s and N b denote the number of signal background events. Most of the variation depends mostly on the amount of missing momentum and the ability to b-tag the event (to further suppress the Z → τ τ in exchange of a significant drop in the b-tagging efficiency). From Table 2 we observe that in the degenerate window; the most important cut is on the number of b-jets; since this region of parameter space is already characterized by a large amount of / E T . In the non-degenerate window, a stronger reach can be achieved by keeping events with lower / E T . Therefore; we would like to point out that b-tagging is only relevant in regions where small MET is required to increase the number of signal events and to significantly suppress the background. Requiring two b-jets will only destroy any possibility to extract the signal. In Fig. 10 we show the entire light colored scalar region in the m φ -m X plane that the LHC Run II may be able to probe with 10, 30, and 300 fb −1 , respectively.

Realizations
The framework studied in this work fits naturally into the light stop window scenario studied in [34]. In that work, the MSSM is consistent with collider data for stop masses above 200 GeV, naturally emerges from renormalization group (RG) evolution, predicts the correct dark matter relic abundance and it is consistent with flavor constraints. This scenario can be accommodated through a light mostly righthanded stop with mass between 200-400 GeV (which corresponds to φ in this work), a heavy mostly left-handed stop, a gluino with mass below 1.5 TeV and a light neutralino (here, X ) with nearly degenerate with the light stop.
Given an LSP neutralino with a mass very close to the stop mass, m ≡ mt 1 − m X ≈ 20-50 GeV, besides the two-body flavor-violating decayt 1 → c X, the four-body decays are also relevant.
where θ tc is the stop-scharm mixing angle, which in [34] is estimated to be of order O (10 −5 ). We see that these decays yield decay lengths between t ≈ 0.01-0.1 mm. However,  Fig. 11a, we show the shaded region consistent with decay lengths in the range 0.001-1 mm in the absence of a flavor-violating decay mode. Larger lifetimes will require an additional suppression of the stop four-body decay mode. This suppression may exist in extensions of the MSSM where the lightest neutralino is mostly a SM gauge singlet. In Fig. 11b, the shaded regions denote decay lengths in the range 0.1-100 mm after suppressing the stop-neutralino coupling by an O(0.01) coupling. This region of parameter space is then a driving force motivating an analysis for displaced vertices such as the one discussed in Sect. 4 to further probe this light-stop window of the MSSM for dark fermion masses above ∼50 GeV.

Conclusions
We have studied a simplified scenario where a colored electroweak-singlet scalar couples to a dark sector through the right-handed top quark. Within our framework the colored scalar can decay after traveling through the inner detector leading to displaced leptons and missing energy. The amount of missing energy strongly depends on the spectrum of the model, but it can be sufficiently large to escape current bounds on displaced vertices that mostly target the parameter space of R-parity-violating scenarios of the MSSM. After implementing a recent search by the CMS collaboration we observe that a color scalar with lifetime in the range ∼0.1-1 mm/c and masses between ∼200-350 GeV may have escaped detection. We have proposed a search strategy at the 13 TeV LHC that can probe a large enough amount of the parameter space with as little as 10 fb −1 by implementing a cut on the lepton's impact parameter and requiring a large amount of missing transverse energy. However, a statistical significant probe will require between 100-300 fb −1 .
In addition, we observe that our simplified model naturally fits into light-stop window scenarios of supersymmetric extensions of the SM, where we observe that flavor-violating mixing angles below 10 −5 can lead to purely displaced vertices between 0.1-10 mm for colored scalar masses between 100-250 GeV. This is incredibly exciting since it opens the possibility that a slightly long-lived right-handed stop could be spotted at the LHC.