Constraining Quirky Tracks with Conventional Searches
Marco Farina, Matthew Low

TL;DR
This paper reinterprets existing collider searches and proposes a new method to constrain quirks, hypothetical particles with unique confinement properties, using data from the CMS experiment, achieving TeV-scale mass limits.
Contribution
It introduces a novel reinterpretation of collider limits and proposes a new quirk search strategy in magnetic-field-less CMS data, extending constraints on quirks.
Findings
Mass constraints in the TeV range depending on quantum numbers.
Reinterpretation of monojet and stable charged particle searches.
Proposal of a new quirk search method in CMS data.
Abstract
Quirks are particles that are both charged under the standard model and under a new confining group. The quirk setup assumes there are no light flavors of the new confining group so that while the theory is in a confining phase, the distance between quirk-antiquirk pairs can be macroscopic. In this work, we reinterpret existing collider limits, those from monojet and heavy stable charged particle searches, as limits on quirks. Additionally, we propose a new search in the magnetic-field-less CMS data for quirks and estimate the sensitivity. We focus on the region where the confinement scale is roughly between 1 eV and 100 eV and find mass constraints in the TeV-range, depending on the quirk's quantum numbers.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3| cut | HSCP | monojet | 0T |
|---|---|---|---|
| GeV | GeV | GeV | |
| m |
| spin | charge | ||||
|---|---|---|---|---|---|
| ( 1 eV) | ( 100 eV) | ( eV) | |||
| fermion | 2 | 1.6 TeV | 1.0 TeV | 500 GeV | |
| scalar | 2 | 1.3 TeV | 700 GeV | 350 GeV | |
| fermion | 2 | 650 GeV | 150 GeV | ||
| scalar | 2 | 350 GeV | 60 GeV | ||
| fermion | 5 | 1.8 TeV | 1.1 TeV | 600 GeV | |
| scalar | 5 | 1.4 TeV | 850 GeV | 450 GeV | |
| fermion | 5 | 800 GeV | 200 GeV | 30 GeV | |
| scalar | 5 | 450 GeV | 80 GeV |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Constraining Quirky Tracks with Conventional Searches
Marco Farina
New High Energy Theory Center, Department of Physics, Rutgers University,
136 Frelinghuisen Road, Piscataway, NJ 08854, USA
Matthew Low
School of Natural Sciences, Institute for Advanced Study,
Einstein Drive, Princeton, NJ 08540, USA
Abstract
Quirks are particles that are both charged under the standard model and under a new confining group. The quirk setup assumes there are no light flavors of the new confining group so that while the theory is in a confining phase, the distance between quirk-antiquirk pairs can be macroscopic. In this work, we reinterpret existing collider limits, those from monojet and heavy stable charged particle searches, as limits on quirks. Additionally, we propose a new search in the magnetic-field-less CMS data for quirks and estimate the sensitivity. We focus on the region where the confinement scale is roughly between 1 eV and 100 eV and find mass constraints in the TeV-range, depending on the quirk’s quantum numbers.
Introduction — The Large Hadron Collider (LHC) has now been running for several years and continues to be our most direct probe of electroweak-scale physics. The primary directions of phenomenological studies have been naturalness-driven and signature-driven. Along the signature-driven direction, only relatively small developments have been made in the study of unusual particle tracks. Track reconstruction at colliders relies on the simple assumption that all particles follow simple helical trajectories characteristic of the motion of charged particles in a magnetic field. There are new physics scenarios, however, that transcend that assumption and give rise to much stranger types of tracks at the LHC. Some examples of these track signatures include tracks that abruptly change direction (kinked tracks), tracks that begin partway through the detector (appearing tracks), tracks with anomalous deposits of energy, and tracks with unusual curvature (see Burdman et al. (2008); Harnik and Wizansky (2009); Asai et al. (2011); Harnik et al. (2011); Meade et al. (2012); Fok and Kribs (2011); Jung and Lee (2015); Knapen et al. (2016) for past and related theory studies). The latter case is typical of quirks, which will be the focus of this work.
Quirks are particles that are both charged under the standard model (SM) and under a new confining group Okun (1980a, b); Kang and Luty (2009). The quirk setup assumes there are no light flavors of the new confining group so that while the theory is in a confining phase, the distance between quirk-antiquirk pairs can be macroscopic. This leads to an interesting array of collider signatures based on the length, , of the flux tube, or string, between the quirks.
It is when becomes comparable to the length scales relevant for detectors that quirk tracks exhibit unusual curvature. Due to the challenges in identifying such tracks, there have been very few dedicated searches performed for quirks. A search from DZero sets the only bound on quirks which is when (where the individual quirks are not resolved) Abazov et al. (2010).
In this work, we will show for the first time that strong bounds can be set on quirks, at collider-relevant scales, using entirely standard LHC searches with no modifications to tracking algorithms. These searches are sensitive for macroscopic string lengths. In addition to reinterpreting existing searches, we propose a new search that can be performed in the magnetic-field-less “0T” data of CMS (still using standard tracking algorithms). In the 0T data all known charged particles are expected to leave straight tracks making this dataset a nearly background-free sample for certain types of tracks with anomalous curvature. While we propose a specific search for quirks, we are optimistic that the use of such a dataset can be extended to other scenarios beyond quirks.
The remainder of the paper is organized as follows: first, we briefly describe quirk models and their collider phenomenology. Then, we reinterpret monojet searches and heavy stable charged particle (HSCP) searches in the quirk parameter space, leveraging the fact that hadron colliders automatically scan a range of largely due to the velocity distribution of particles produced. Next, we suggest a novel use of the 0T data from CMS to constrain quirks. Finally, we conclude with projected results and comments on the remaining quirk parameter space.
Quirks (at the LHC) — We now introduce the minimal ingredients for a quirk model. To the SM gauge group we add a new “infracolor” gauge group that is assumed to be asymptotically-free with a confinement scale and to the SM particle content we add a new species, , with mass and infracolor representation size . The particle is called a quirk when it is much heavier than the confinement scale (). Since is assumed to be the lightest infracolored particle, there are no particles lighter than that can form “hadrons” and the only hadronic states are glueballs with masses a little above . When quirks are pair produced, for instance at the LHC, there are no light hadrons to break the infracolor flux tube between the quirks. This flux tube, or string, connecting the pair can be macroscopic in size and its tension results in an attractive force between the two particles.
If the quirks are colored they hadronize, via QCD, forming color-neutral states with SM quarks or gluons, similar to the well-studied case of -hadrons in supersymmetry (we will adopt the name -hadron for the color-neutral quirk-parton state). The two -hadrons are still connected by the infracolor flux tube and for all practical purposes the effect of QCD hadronization on the quirk dynamics can be neglected as . On the other hand, the electric charge of the -hadron, does affect its detection; we will return to this point in the discussion of our results. Clearly, tracking methods cannot be used for neutral -hadrons.
We can now study the trajectory of quirks at the LHC. The equations of motion are given by the Nambu-Goto action with point masses on the ends of the string in an electromagnetic background Kang and Luty (2009); Luscher and Weisz (2002). The equations of motion for a single quirk are
[TABLE]
Above is a unit vector that points towards the other quirk and is used to define , , and . The quirk charge is and the magnetic field is . For both the HSCP and monojet searches we use the CMS magnetic field of while for the 0T search we use . The tension is given by which is proportional to . There have been estimates that in QCD Deur et al. (2015), but we take for simplicity (as the difference is simply a rescaling of the parameter space).
In the absence of external forces and when the quirks are back-to-back the maximum distance between them can be calculated to be
[TABLE]
where is the boost factor. While the true string length, , can be different, we use as a simple approximation. The factor follows from dimensional analysis and the factor plays a relevant role in collider searches. Numerically, one has
[TABLE]
From Eq. (3) one can map different types of searches to the appropriate range of . For the quirks only separate a microscopic distance, comparable to the typical tracking resolution () so that the quirk pair is observed as a single highly-ionizing straight track Abazov et al. (2010). For one finds that is macroscopic and leads, in general, to oddly curved tracks.111Note that the range roughly corresponds to the length range . The distance is still relevant for the LHC because the sagitta of an LHC track is roughly where is the radius of the tracker. Taking the sagitta to be comparable to the tracker resolution , one finds that there is sensitivity up to . Finally, for the effective string length is megascopic and does not play a role in collider searches, leaving the quirks to appear as stable charged particles.
By oddly curved, we simply mean that the quirk tracks deviate from the standard helix that curves in the -plane. Quirk tracks can exhibit a wide variety of strange behaviors, for example, curving outside the -plane, reversing direction, or passing through the same detector layer more than once. Crucially, however, they can also closely approximate a standard helix trajectory (at least inside the detector volume). The factor in Eq. (3) allows for a range of values for a given .
Reinterpreting Existing Searches — As mentioned above, for a wide range of there is a non-zero probability that a quirk track would be reconstructed at the LHC. When this happens the quirk will appear simply as a heavy stable (or long-lived) charged particle. In a collider, such particles are found by looking for tracks with large deposits of energy and/or a long time of flight (as compared to muons). When both tracks fail to be reconstructed monojet searches will have sensitivity, provided that the quirks have been produced in association with a sufficiently energetic jet (through either initial or final state radiation). Monojet searches look for large missing transverse energy that results from a jet recoiling against undetected particles.
The probability to reconstruct a track is characterized by the track efficiency and is shown in Fig. 1 for HSCP searches (red) and for monojet searches (blue). The track efficiency is computed by applying a series of track selection cuts (as used by CMS in their HSCP analysis CMS (2016a)) listed in Table 1 and described below.
First, the quirk propagation is found by solving Eq. (1). We use the straight string approximation throughout. Each time a quirk passes through a tracker layer it registers as a hit with an efficiency that we assume to be 100%. To account for the fact that in practice the hit efficiency is less than 100% we increase the requirement to from the typical . The requirement specifies the minimum number of layers a track must hit.
We model the tracker geometry following the specifications of the CMS tracker,222We use the CMS tracker to maximize the accuracy of our results for the HSCP and 0T searches. which consists of a combination of barrel layers and endcap layers that cover the range where a particle will pass through layers depending on its trajectory Chatrchyan et al. (2008). Details of our tracker model are given in the appendix.
The tracks are then fit to a helix, which corresponds to the trajectory of a charged particle in the longitudinal magnetic field of the detector. A helix is given by
[TABLE]
where is the parameter along the curve, is the radius, is the initial azimuthal direction, and is the dip angle. For a completely general helix there are 3 additional parameters specifying the initial position of the helix, but we set this to the origin for simplicity. The depends on the charge of the particle.
Next, we perform a fit to Eq. (4) and assume that quirks with would be reconstructed as tracks CMS (2016a). We use a spatial resolution of for each hit Chatrchyan et al. (2008). The cut is applied to the measured of the track (computed from the best fit value) rather than the true of the particle.
The first search that we reinterpret is the HSCP search. The most sensitive HSCP searches have been performed at 13 TeV by CMS with 12.9 fb*-1* CMS (2016a) and by ATLAS with 3.2 fb*-1* Aaboud et al. (2016a) and are presented as mass limits on stable particles of different charges and SU(3) representations. We follow the event selection of the former analysis which primarily consists of a cut on the track and a cut of . We only consider the sample that would be selected by the muon trigger which adds the additional requirement that in order that the trigger be efficient. We generate quirk pair events using Madgraph 5 v2.3.3 Alwall et al. (2014). The tracking efficiency is shown in Fig. 1 (red). We define the track efficiency as the number of tracks passing all track requirements divided by the number of tracks passing the cut. As expected, the efficiency drops as the string length becomes comparable to the detector scale, while it approaches 100% when the string tension is small compared to the Lorentz force. The asymptotic value at low is determined by primarily by the cut.
The second search we reinterpret is the monojet search. From CMS the strongest search is at 13 TeV and uses 12.9 fb*-1* CMS (2016b) and from ATLAS it is at 13 TeV using 3.2 fb*-1* Aaboud et al. (2016b). While the CMS search has better reach, we use the ATLAS limits because they are presented as limits on compressed stop squarks which is kinematically more similar to quirk events than the setups in the CMS search. We also show limits scaling the ATLAS result up to 12.9 fb -1. We generate quirk pair events along with a radiated jet of and follow the monojet event selection in Aaboud et al. (2016b). The track selection in this case still identifies when a good track would be selected, however, opposite to the HSCP case, this means the event would be rejected. For this reason we plot the quantity in Fig. 1 (blue). Note that even for small the efficiency does not go to zero because some tracks still fail the selection criteria.
Using both HSCP and monojet searches provides a very complementary approach so that while our simplified tracker parametrization could differ from a full simulation it still captures the crucial features. In particular, the complementarity ensures that our results reliably cover the full range of we study. When one efficiency degrades, the other reaches its maximum (and at the maximum the tracking reconstruction details are less important). Another important feature captured in Fig. 1 is the weak dependence on the quirk mass (since is approximately linear in mass). In a similar manner, the quirk’s electric charge has a minor impact on the efficiencies. Indeed, while we use a fermionic color triplet quirk with quantum numbers under the SM gauge group as a case study, we will provide limits for a few other cases. The results are shown in Fig. 3 which will be further explained in the following section.
Using the 0T Data — In addition to the reinterpreted searches, we propose an entirely new search whose sensitivity is maximal in the challenging region near . This search makes use of the 0.6 fb*-1* of 13 TeV data with a 0T magnetic field Khachatryan et al. (2016). Without a magnetic field, all SM particles, both neutral and charged, travel in straight lines. Quirks, on the other hand, still curve due to the string tension. This means one can simply count the number of curved tracks in the 0T data and accordingly set a limit or make an observation of quirks. Operationally, this would entail running the tracking algorithm on the 0T data while pretending there is a magnetic field and should not require any modifications to the tracking algorithm itself.
The efficiency for identifying a track in the 0T data is shown in Fig. 2. For the 0T search we use the monojet trigger CMS (2016c) which uses an analysis-level cut of for the leading jet and where excludes muons (and would also exclude the quirk tracks). While in principle one could use the muon trigger as in done in the HSCP searches, in practice the monojet trigger is more effective. The reason is that without sizable initial radiation for the quirk system to recoil against, the quirks will be almost back-to-back in the transverse plane. This means there is no curvature in the -plane so that the quirk trajectories cannot be reconstructed as non-straight helices. Therefore, at least some initial state radiation is required for a non-zero efficiency.
We generate events with a single jet with and apply the track cuts in Table 1. These closely follow the selection from the HSCP search with the exception that we add a requirement that the fitted radius must be . We estimate the value from the sagitta of a track
[TABLE]
where is the chord length, corresponding to the radius of the tracker. We take Chatrchyan et al. (2008) and set the sagitta to the single hit resolution of Chatrchyan et al. (2008) and require a single hit fluctuation to find . Note that estimate assumes 3 hits. Since a straight track faking a quirk requires hits, the fake rate is much lower than indicated by the requirement. The cut is responsible for the decrease in track efficiency at .
We assume that the fake rate is sufficiently low and that multiple scattering effects faking a curved track are sufficiently rare to treat the analysis as zero-background. In principle, if a track with non-zero curvature is discovered the event could be inspected and checked for the presence of a second curved track, providing a smoking gun of the signal. A limit is projected corresponding to observing events Feldman and Cousins (1998).
Discussion of Results — The results are shown, for a fermion with , in Fig. 3. The shaded red region shows the limits from HSCP searches which drive the limits for . The shaded green region shows the limits from the ATLAS monojet search that used 3.2 fb*-1* and the unshaded green limit scales up the limit to a dataset of 12.9 fb*-1* (the amount used in the CMS monojet search). Our projection for the 0T search is given by the shaded blue region which uses 0.6 fb*-1*. The unshaded blue line shows a hypothetical dataset of 20 fb*-1* with no magnetic field and is the minimum amount of data required to probe parameter space beyond what is covered by HSCP and monojet searches. The HSCP and 0T bounds are cut off at because our statistics there are insufficient for a reliable estimate.
Regarding QCD hadronization, the HSCP searches at the LHC use two different models CMS (2016a); Aaboud et al. (2016a). The first model Kraan (2004); Mackeprang and Rizzi (2007) assumes that the heavy hadrons can be charged or neutral when exiting the calorimeter while the second model Mackeprang and Milstead (2010) assumes the all heavy hadrons are neutral when exiting the calorimeter. The 0T search only uses information from the tracker and therefore does not depend on this assumption. We take the fraction of -hadrons that are charged to be from Pythia 8 Sjöstrand et al. (2015).
For the 0T and HSCP searches we allow for 1 or 2 identified tracks while for the monojet search we require 0 identified tracks. The overall efficiency includes an acceptance factor ( in the relevant region) that was not used in Figs. 1 and 2. We find that HSCP can constrain higher masses than monojet searches or the 0T search. This is not surprising as both the monojet and 0T searches require an additional radiated jet and have a larger background or smaller dataset.
In Table 2 we report the limits for various other quantum numbers using , , and as benchmark points.
The gray lines in Fig. 3 show contours of constant . The used to compute is the mean of the velocity distribution at 13 TeV. The contours give an idea of the length scales where each search is most effective and conversely, show where in parameter space other types of quirk searches would lie.
For , smaller than shown in Fig. 3, the 0T limit quickly goes away but the HSCP limit stays constant at 1.6 TeV because simply gets larger. On the other end, for , the monojet limit stays constant until where the quirk system will appear as a single straight track. Here, HSCP searches might have some sensitivity again. For large enough eventually quirk-antiquirk annihilation can become prompt and searches for various resonances, like or can be relevant Kang and Luty (2009).
Before concluding we note that when a quirk pair crosses, soft radiation could be emitted as hadrons or glueballs, leading to a loss of energy and modification of the quirk trajectories, as well as extra activity in the detector. In most of the parameter space shown, however, we study so that the quirks typically do not cross each other.
Finally, one can ask how useful would modifications or extensions to tracking techniques be? If it is possible to trigger on such anomalous tracks, one expects to set bounds competitive with HSCP searches and extend limits above . For quirks, however, one limitation at the LHC is that the string length should be larger than the average distance traveled by the quirk between two layers of the tracker, which is . If not, then the sampling of the track would lose information on the peculiar “periodic” structure of the trajectory. It is therefore in the region where substantial improvements are possible. This corresponds to an order of magnitude in and we believe that a more detailed study is warranted, not only at ATLAS and CMS, but also at more specialized detector such as MoEDAL Acharya et al. (2014).
Outlook — In this paper, we demonstrated that while quirk dynamics can result in very exotic tracks, they can also result in very standard looking tracks allowing for standard searches to constrain a substantial region of parameter space. In particular, reinterpreting HSCP and monojet searches allows one to set limits in the regions and respectively, which corresponds to effective string lengths of and .
For colored quirks the limits range from 1.0 TeV to 1.6 TeV with . The limits get correspondingly higher as is increased. Limits on uncolored quirks were found to range from 150 GeV to 650 GeV.
We then proposed a novel use of the 0T data from CMS which involved looking for curved tracks in the dataset. Given the small size of the dataset, 0.6 fb*-1*, we predict that no curved tracks (at least due to quirks) should be observed in the data as HSCP searches already rule out the parameter space. Amusingly, the 0T search could outdo the current HSCP limits if it had only of data.
We chose a few sample quantum numbers, in Table 2, but it would be interesting to see limits for a larger variety of quantum numbers. On the experimental side, it would interesting to see if dedicated quirk searches can be done and would they compare to the monojet and HSCP searches.
Finally, we argued that in the region there is an opportunity to go well beyond the sensitivity of current searches by developing novel tracking techniques. Such techniques could then be applied to other physics cases, for instance, kinked or appearing tracks. Given the simplicity of our 0T analysis, we hope that this work can serve as motivation for moving towards more involved tracking modifications in order to fully exploit the LHC’s potential for unusual tracks.
Acknowledgements.
The authors would like to thank Raffaele Tito D’Agnolo, Markus Luty, Duccio Pappadopulo, Joshua Ruderman, Andreas Salzburg, and Kris Sigurdson for useful discussions and Yuri Gershtein, Philip Harris, Simon Knapen, Scott Thomas, and Nhan Tran for helpful discussions and reading the manuscript. M.F. is supported in part by the DOE Grant DE-SC0010008. M.L. is supported by a Frank and Peggy Taplin membership at the Institute for Advanced Study. This work was partly completed at KITP, which is supported in part by the National Science Foundation under Grant No. NSF PHY11-25915 and at the Aspen Center for Physics, which is supported by National Science Foundation grant PHY-1066293.
Appendix A: Tracker Model Details
We model the tracker following the specifications in Chatrchyan et al. (2008). The tracker is comprised of a cylindrical barrel that surrounds the beam pipe and two disks (end caps) on each side of the barrel. We consider each layer to be of negligible thickness. In cylindrical coordinates we specify barrel layers by their position and extent and disk layers by their position and extent.
The tracker consists of the pixel detector:
- •
Pixel barrel (3 layers)
- •
Pixel end cap (2 layers)
and the strip tracker:
- •
Tracker inner barrel (4 layers)
- •
Tracker outer barrel (6 layers)
- •
Tracker inner disk (3 layers)
- •
Tracker end caps (9 layers)
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Burdman et al. (2008) G. Burdman, Z. Chacko, H.-S. Goh, R. Harnik, and C. A. Krenke, Phys. Rev. D 78 , 075028 (2008) , ar Xiv:0805.4667 [hep-ph] . · doi ↗
- 2Harnik and Wizansky (2009) R. Harnik and T. Wizansky, Phys. Rev. D 80 , 075015 (2009) , ar Xiv:0810.3948 [hep-ph] . · doi ↗
- 3Asai et al. (2011) S. Asai, Y. Azuma, M. Endo, K. Hamaguchi, and S. Iwamoto, JHEP 12 , 041 (2011) , ar Xiv:1103.1881 [hep-ph] . · doi ↗
- 4Harnik et al. (2011) R. Harnik, G. D. Kribs, and A. Martin, Phys. Rev. D 84 , 035029 (2011) , ar Xiv:1106.2569 [hep-ph] . · doi ↗
- 5Meade et al. (2012) P. Meade, M. Papucci, and T. Volansky, Phys. Rev. Lett. 109 , 031801 (2012) , ar Xiv:1103.3016 [hep-ph] . · doi ↗
- 6Fok and Kribs (2011) R. Fok and G. D. Kribs, Phys. Rev. D 84 , 035001 (2011) , ar Xiv:1106.3101 [hep-ph] . · doi ↗
- 7Jung and Lee (2015) S. Jung and H.-S. Lee, (2015), ar Xiv:1503.00414 [hep-ph] .
- 8Knapen et al. (2016) S. Knapen, S. Pagan Griso, M. Papucci, and D. J. Robinson, (2016), ar Xiv:1612.00850 [hep-ph] .
