# Long-Lived Heavy Quarks: A Review

Authors: Mathieu Buchkremer; Alexander Schmidt
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

We review the theoretical and experimental situation for long-lived heavy quarks, or bound states thereof, arising in simple extensions of the Standard Model. If these particles propagate large distances before their decay, they give rise to specific signatures requiring dedicated analysis methods. In particular, vector-like quarks with negligible couplings to the three known families could have eluded the past experimental searches. While most analyses assume prompt decays at the production vertex, novel heavy quarks might lead to signatures involving displaced vertices, new hadronic bound states, or decays happening outside of the detector acceptance. We perform reinterpretations of existing searches for short- and long-lived particles, and give suggestions on how to extend their
reach to long-lived heavy quarks.


---

## Body

## 1. Introduction

Over the past decades, we have discovered that nature consists of a given number of elementary particles, deeply connected with the known fundamental forces driving our universe. The recent observation of a new particle resembling the long-sought Higgs boson at the Large Hadron Collider (LHC) now provides us with strong evidence for the validity of the Standard Model (SM) [1, 2]. Yet, while it is generally acknowledged that the latter comprises three generations of chiral quarks and leptons, various fundamental problems do not find their answers within this framework. For instance, long-standing issues such as the origin of the fermion mass hierarchy or the nature of the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix hint at the possible need for new physics. While the SM is expected to lose its predictive power as the experimental searches now reach the TeV scale, models suggesting new heavy fermions might offer solutions to these problems in the near future.Although a fourth sequential family of quarks is now disfavoured by the recent results of the searches for the Higgs boson [3], other models predicting new heavy particles beyond the top quark are still consistent with the current experimental measurements. In particular, the possibility for vector-like quarks, that is, quarks having their left- and right-handed components transforming identically under the electroweak gauge group, is a common feature in many scenarios going beyond the SM, for example, extra dimensional models, Little Higgs models, grand unified theories, and so forth [4]. Nonchiral quarks also have the peculiarity to decouple in the heavy mass limit, leading to SM-like signals. Should the mixings of these with the light SM fermions be suppressed, the Higgs production rates would not be easily distinguishable from the SM expectations [5]. These new quarks could be sufficiently long-lived to behave like effectively stable particles, evading the current searches as they propagate over sizeable distances.From the collider searches, it is clear that if new states with masses less than a hundred GeV had existed, they would have been observed. On the other hand, the experimental reach for detecting novel heavy stable particles above the TeV scale depends on how readily identifiable such states are at the LHC. While particles with nanoseconds lifetimes could fly away from the primary interaction point before they decay to ordinary particles and lead to displaced vertices, they can also hadronise, allowing for a possibly rich spectrum of new exotic and heavy bound states. Such stable massive particles are anticipated in many new physics models, either due to the presence of a new conserved quantum number (e.g.,R parity in supersymmetric models) or because the decays are suppressed by kinematics or small couplings (see [6] and references therein for an exhaustive review). Although the majority of the past studies focused on stable states in supersymmetry, the possibility for new stable quarks received few attention at the LHC. While the CMS and ATLAS experiments are now setting strong limits on long-lived gluino, stop and stau pair production (see, for instance, [7, 8]), scenarios involving quasi-stable heavy quarks have been barely investigated at the LHC.In the following, we review the phenomenology of new, long-lived heavy quarks. Assuming a general parameterisation, Section2 first examines their production and decay modes at the LHC, considering both the chiral and vector-like scenarios. The possible signatures for displaced vertices and stable massive hadrons are then covered within the context of negligible mixings with the SM fermions. The experimental aspects of the collider searches for new long-lived quarks are presented in Section 3. We review the situation for prompt decay searches, displaced vertices, and Heavy Stable Charged Particles (hereafter, HSCPs) and highlight some of their limitations. Reinterpretations of the existing searches are then presented for short- and long-lived particles. Considering the possibility to improve the sensitivity of the existing analyses to long-lived quarks beyond the TeV scale, some alternative search topologies are finally described. Our conclusions are given in Section 4.


## 2. Long-Lived Quarks Signatures

### 2.1. Production

As we review the possibility for long-lived quarks in general, two different scenarios will be distinguished in this work. New chiral fermions, for which the left- and right-handed chiralities have different charges under the electroweak gauge group, gain their masses from the electroweak symmetry breaking mechanism. If heavy, a sequential fourth generation quark doublet(t′,b′)L can therefore induce large corrections to loop observables. Such a scenario is now severely constrained from the electroweak precision data and the Higgs search results. Given that chiral fermions couple to the Higgs boson with a strength proportional to their Yukawa couplings, they do not decouple from its production, as the corresponding rate should increase by a factor of about 9 due to the additional fermion loops occurring in gluon fusion. The unobserved enhancement thus now strongly disfavours a fourth generation Standard Model [3].Nonchiral quarks, on the other hand, still provide a viable extension to the SM and certainly require a dedicated discussion. Considering the introduction of new vector-like heavy partners mixing with the SM quarks, additional parameters are allowed fromSU(2)L×U(1)Y invariant Yukawa interactions and Dirac mass terms in the SM Lagrangian. The light quark couplings to the Higgs and electroweak bosons are consequently modified, along with the introduction of new couplings between the heavy and the SM quarks.Still, it is important to emphasise that such new states cannot have arbitrary quantum numbers, since they can only mix with the SM quarks through a limited number of gauge-invariant couplings. Classifying them intoSU(2)L multiplets, their Yukawa terms only allow for seven distinct possibilities, that is, the two singlets, the three doublets, and the two triplets displayed in Table 1.Table 1
Vector-like multiplets allowed to mix with the SM quarks through Yukawa couplings. The electric charge is the sum of the third component of isospinT3 and of the hypercharge Y. Adapted from [9].
Q
q
T
2
/
3
B
-
1
/
3
(
X
5
/
3
T
2
/
3
)
(
T
2
/
3
B
-
1
/
3
)
(
B
-
1
/
3
Y
-
4
/
3
)
(
X
5
/
3
T
2
/
3
B
-
1
/
3
)
(
T
2
/
3
B
-
1
/
3
Y
-
4
/
3
)
T
3
0
0
1
/
2
1
/
2
1
/
2
1
1
Y
2
/
3
-
1
/
3
7
/
6
1
/
6
-
5
/
6
2
/
3
-
1
/
3The production of such new heavy quarks, either chiral or vector-like, is usually assumed to proceed dominantly at the LHC through gluon fusion,gg→QQ¯. Nevertheless, depending on the model at hand, electroweak single production can also provide an alternative mechanism as it is not affected by the large phase-space suppression from which pair production suffers. In particular, new heavy quarks can be produced singly in flavour-changing processes via the electroweak interaction through qiqj(-)→V*→qkQ, where V=W,Z. A comparison between the dominant production modes can be found in [10, 11] for fourth generation and in [12] for vector-like quarks. Assuming order unity couplings, benchmark cross-sections have been evaluated for various mass scenarios and reproduced in Table 2.Table 2
Electroweak single production cross-sections (in pb) as a function of the massmQ of the heavy quark, assuming order unity couplings as computed in [12]. In contrast, strong pair production pp→QQ- yields 3.61 (0.63) fb for mQ=900 and 1800 GeV/c2 at s=7 and 14 TeV.
Channel
s
=
7 TeV
s
=
14 TeV
m
Q
=
900 GeV/c2
m
Q
=
1800 GeV/c2
p
p
→
W
*
→
j
X
1.4
0.36
p
p
→
W
*
→
j
X
-
0.037
0.0092
p
p
→
W
*
→
j
T
0.61
0.16
p
p
→
W
*
→
j
T
-
0.052
0.013
p
p
→
Z
*
→
j
T
0.43
0.11
p
p
→
Z
*
→
j
T
-
0.025
0.0064
p
p
→
W
*
→
j
B
0.69
0.18
p
p
→
W
*
→
j
B
-
0.089
0.022
p
p
→
Z
*
→
j
B
0.18
0.047
p
p
→
Z
*
→
j
B
-
0.034
0.0088
p
p
→
W
*
→
j
Y
0.29
0.074
p
p
→
W
*
→
j
Y
-
0.12
0.031As described in [12], the leading contributions to vector-like T and B single production arise from ud→W*→jT, du→W*→jB, qd→Z*→jB and qu→Z*→jT, where q denotes a valence quark parton and j a generic light quark jet. Unlike QCD pair production, the above processes, however, scale with the Q-q quark coupling squared and can be strongly suppressed if the associated mixings are negligibly small. On the other hand, Figures  2–4 in [12] indicate that σ(pp→QQ¯) falls off faster than the single production cross-sections as soon as mQ reaches the TeV scale, while the electroweak production channels displayed in Table 2 are enhanced by a factor O(mQ2/mV2), originating from the longitudinal polarisation of the gauge boson V. Furthermore, the relative rates of the jQ and jQ¯ channels strongly depend on the valence quark density in the initial state. Given the difference in the PDFs of valence and sea quarks in the initial states, new heavy quarks are produced singly at a much higher rate than antiquarks for increasing mQ. While the u quark contributes more than the d quark, the neutral current contributions are weaker than the charged current ones. Finally, we underline the absence of the tree-level processes pp→Z*→jQ for vector-like quarks involving exotic charges, given that they only interact with other states via tree-level charged currents [12]. Although most of the above statements are model dependent, we conclude that the large running energy of the LHC makes single electroweak production a promising discovery channel when considering new heavy quark searches with mQ≳1 TeV/c2.


### 2.2. Decay

In the search of new fermions, the associated decay modes require careful attention. The exact partial width for a new sequential heavy quarkQ decaying on-shell to a light quark q through a charged current can be written as
(1)Γ(Q⟶qW)=GFmQ38π2|κQq|2f2(mqmQ,mWmQ),(2)f2(α,β)=[(1-α2)2+β2(1+α2)-2β4]×[1-(α+β)2][1-(α-β)2],
where κQq denotes the generic Q-q quark coupling, equal to the CKM mixing matrix element VQq when considering a new sequential family of quarks. Interestingly, the width (1) holds for both chiral and vector-like singlet quarks, given that all heavy-to-light charged current quark decays occur to be pure V-A processes with identical rates [13]. If mQ≫mq, the above partial width can be shown to reach the asymptotic form
(3)Γ(Q⟶qW)≃170.5MeV×|κQq|2×mQ3mW3.
Considering a new heavy quark Q decaying exclusively via the above 2-body decay mode, we evaluate in Table 3 the corresponding decay lengths, with cτ≃2×10-10 GeV/Γ(GeV) μm and |κQq|=1. With increasing mass, the probability for a new heavy quark Q to decay weakly thus becomes more and more important as its lifetime decreases. However, small couplings to the lighter SM quarks strongly affect this statement as the corresponding widths are suppressed by |κQq|2. As we will see in Section 2.4, this could lead to a very different phenomenology as the quarks Q, if very heavy, would form bound states and possibly decay hadronically.Table 3
Partial widths and decay lengths for a new heavy quarkQ decaying to light quarks via on-shell charged currents Q→qW, with mq=0 and 173.2 GeV/c2 and |κQq|=1, as computed from (1).
m
Q (GeV/c2)
Γ (GeV)
c
τ
×
1
0
12 (μm)
300
8.73 (2.39)
22.91 (83.85)
400
20.90 (11.14)
9.57 (17.96)
500
40.94 (27.86)
4.89 (7.18)
600
70.81 (54.52)
2.82 (3.67)
700
112.50 (93.06)
1.78 (2.15)
800
167.97 (147.43)
1.19 (1.38)
900
239.18 (213.57)
83.62 (9.36) ×10-1
1000
328.12 (299.46)
6.10 (6.68) ×10-1
1100
436.74 (405.05)
4.58 (4.94) ×10-1
1200
567.02 (532.31)
3.53 (3.76) ×10-1
1300
720.93 (683.21)
2.77 (2.92) ×10-1
1400
900.43 (859.72)
2.22 (2.33) ×10-1
1500
1107.5 (1063.79)
1.81 (1.88) ×10-1
1600
1344.1 (1297.4)
1.49 (1.54) ×10-1
1700
1612.2 (1562.5)
1.24 (1.28) ×10-1
1800
1913.8 (1861.1)
1.05 (1.07) ×10-1Several conclusions can already be drawn at this stage. As far as a heavy fourth quark family(t′,b′)L is concerned, we highlight that the ratios of the charged current decay rates to the lighter families only depend on the off-diagonal mixing elements. Assuming that the extended 4×4 CKM matrix is unitary, the CKM mixing elements |VQq| imply nonunity branching ratios for fourth generation quarks in general.The possibility for tree-level Flavour Changing Neutral Currents leads to similar conclusions for models involving new vector-like quarks interacting with the SM families. Depending on their Yukawa couplings, they can undergo tree-level transitions to lighter quarks andW, Z, and Higgs bosons. Considering vector-like quarks Q=T,B with mQ≫mq, the neutral currents Q→qZ and Q→qH can occur with a rate of the same magnitude as given by (3), if they are allowed to mix with the light SM quarks q. If no assumption is made on the hierarchy of the Q-q coupling strengths, no exclusive decay mode to the first, second, or third family should be preferred.For nonchiral quarks decaying toWq, Zq, and Hq, long lifetimes can be expected in the |κQq|≃0 limit. Additionally, it is interesting to notice that Γ(Q→qZ)/Γ(Q→qH)≃1 holds in the large mass limit mQ≫mq for most scenarios, while the charged-to-neutral current ratios Γ(Q→qW)/Γ(Q→q(Z,H)) approach either 2 or 0 depending on the vector-like representation [14].If more than one new heavy quark is considered, and if they are allowed to decay into each other, the partial width (3) becomes inaccurate due to the possible competition with the heavy-to-heavy transitions Q1→Q2V(*), where V=W,Z,H, depending on the model at hand. Indeed, if the decay rates of the daughter particles are significant with respect to the mass of the parent, the decay of an unstable particle is allowed through threshold effects, even if kinematically forbidden. For both chiral and vector-like quarks, real and virtual V emission near threshold proceed with the rate
(4)Γ(Q1→Q2V(*))=GF2mQ15192π3|κQ1Q2|2f3(mQ12mV2,mQ22mV2,ΓV2mV2),f3(α,β,γ)=2∫0(1-β)2dx‍×[(1-β)2+x(1+β)-2x2]λ(1,x,β)[(1-αx)2+γ2],
where λ(a,b,c)=a2+b2+c2-2ab-2bc-2ac [4]. While the Q→qV rates essentially depend on the Q-q quark couplings, the heavy quark mass splittings |mQ1-mQ2| drive the strength of the heavy-to-heavy transitions [11, 15]. If the quarks Q1 and Q2 belong to the same weak multiplet with a large mass splitting, the Q1→Q2V decay could be very rapid and dominate the other decay modes if there is no |κQ1Q2| suppression. On the other hand, the lightest partner can only decay into an SM family quark and is thus plausibly long-lived for small mixings. Should this be the case, we will see in Sections 2.4 and 2.5 that the latter could hadronize, leading to a different decay phenomenology. This occurs specifically for SU(2)L vector-like singlets, which free quark decays can only proceed via nonzero mixing into the lighter generations.As far as the doublet and triplet representations are concerned, new vector-like partners should be very close in mass if one assumes that isospin conservation is respected. Although it may be broken by higher order corrections, such a degeneracy has been shown in [16] to require a mass difference between 70 and 110MeV for nonchiral doublets, allowing for nanoseconds lifetimes. If their couplings to the lighter quarks are negligible, the partial widths for X5/3, T2/3, B-1/3, and Y-4/3 quarks are thus plausibly small.Summarising these results, the range of the observable couplings to which long-lived particle searches are sensitive at the LHC can be evaluated. From the approximation (3), the partial decay length for a vector-like quark reads
(5)λ≃2cm×(10-8|κQq|)2(500GeV/c2mQ)3.
Assuming mQ≲1 TeV/c2 and an experimental resolution of at least 25 μm, the heavy quark decays are prompt only if their associated couplings are above the 10-7 level. If |κQq|≲10-7, the corresponding decay lengths could be observed over a few tens of centimeters. In the limit of such small mixings with the SM fermions, it is thus natural to motivate searches for displaced vertices.


### 2.3. Classification of Signatures

We now detail the possible situations giving rise to the observable signatures of long-lived quarks at the LHC.The classification given in Table4 summarises the long-lived quark scenarios corresponding to the small mixing scenarios treated in this work. We emphasise that if mixing with all SM quark families is permitted, direct searches should aim at being sensitive to all light quarks q in the observed final states. Given that their interactions are allowed through arbitrary Yukawa couplings, the searches for long-lived vector-like quarks should be as inclusive as possible in order to cover the full spectrum of possibilities at the LHC. Such a program is important to set new constraints on the Q→qV(*) decay modes, independently of any assumptions on the Q-q mixing sector. The relative importance of the Q→(u,c,d,s)V channels is most relevant compared to the Q→(t,b)V charged currents, as no exclusive decay mode should be preferred for V=W,Z,H.Table 4
Possible decay signatures for new long-lived quarks (mQ=1 TeV/c2), as discussed in the text.
(i)
(ii)
(iii)
Short-lived
Intermediate
Long-lived
|
κ
Q
q
|
≳

1
0
-
7
≲

1
0
-
7
≳

1
0
-
9
≲

1
0
-
9
Γ (GeV)
≳

1
0
-
12
≲

1
0
-
12
≳

1
0
-
16
≲

1
0
-
16
τ (s)
≲

1
0
-
13
≳

1
0
-
13
≲

1
0
-
9
≳

1
0
-
9As illustrated in Figure1, three cases for new long-lived quarks can be distinguished depending on their mixing parameters with the SM fermions. (i)
Firstly, if all couplings with the known SM fermions are larger than the10-7 level, all heavy-to-light decays Q→qV are prompt, over distances smaller than 25 μm. This defines the short-lived scenario (i), with no experimental difference compared to the direct searches.(ii)
The second scenario (ii) arises for intermediate decay lengths ranging between a few microns and centimetric distances. If there is at least one light quarkq such that 10-7≲|κQq|≲10-9, displaced events could be observed, yet with a partial width suppressed by |κQq|2. Consequently, if all the couplings but one are below the 10-7 level, a single exclusive channel would lead to a prompt decay signature corresponding to a short-lived heavy quark. A possible exception, though, is the case for new heavy multiplets with sizeable mass splittings, as allowed in extra generation models and possible extensions [11, 17–19]. If mQ1≳mQ2+mV, the heaviest quark Q1 can be short-lived and decay semiweakly to Q2V(*), while the lightest partner is likely stable if all its decay modes suffer severe suppression. If mQ1≃mQ2, on the other hand, all heavy-to-heavy transitions are suppressed and both quarks could be long-lived on detector scales.(iii)
For particles with decay lengths larger than the detector dimensions lies the long-lived region (iii), the details of which will be discussed in Section2.4. If all heavy quark couplings with the SM fermions are below the 10-9 level, the stable case becomes a relevant scenario, possibly in conjunction with displaced events.Figure 1
Lifetimeτ in seconds (γβ≃1) for a sequential heavy quark Q decaying to the lighter families by emitting a real W boson, for mQ=300 GeV/c2 (in blue) and mQ=1800 GeV/c2 (in red). The three regions of interest refer to the conventions given in Table 4.We emphasise that the scenarios for which prompt decays and displaced vertices could be observed only call for one of the allowed decay modes to fulfill the associated conditions. The stable scenario, however, requires all decay lengths and partial widths to lie in the ranges given in Table4. As we will see in Section 2.4, if all Q-q quark couplings to the SM families are smaller than 10-2, such new heavy fermions could hadronise. As a result, the annihilation decays and hadronic transitions between the formed bound states would dominate, while the single quark decay events might be suppressed.


### 2.4. Signatures for Heavy Stable Particles: Heavy Quarkonia

Previously, we have discussed how displaced events could occur in the single quark decay transitions in the case of smallQ-q couplings. In this section, we describe the possibility for new heavy quarks to bind into hadrons if their lifetime is large enough, providing us with an interesting alternative for new signal searches at the LHC.While new coloured particles with masses heavier than the top quark are usually assumed to decay as free particles, their possible formation into baryons and mesons is an important issue to consider in the case of small couplings with the SM fermions.Should the partial widthsQ→q(W,Z,H) be suppressed, the quark Q might be stable and hadronise. Assuming that the binding force in a fourth generation bound state is of Coulombic nature, the seminal condition
(6)mQ<125GeV(100GeV)×|VQq|-2/3
was derived in [20] for QQ¯ (Qq¯) formation. According to (6), new heavy quarks form hadronic states if their mixing with known quarks is sufficiently small. Mixings roughly smaller than 0.2 are typically required for 300 GeV/c2 new quarks to form bound states, which remains perfectly consistent with the current bounds from the electroweak precision observables. As shown on Figure 2 for a new up-like fourth generation quark t′, the electroweak precision fits typically restrict the mixing matrix element |Vt′b|≃|Vtb′| to be smaller than 10-1 for mt′>700 GeV/c2. If lighter, the t′ total width can be smaller than 10MeV [11], thus allowing new chiral quarks to form bound states with nanosecond lifetimes.Figure 2
Reproduction of theχ2 fit of the oblique parameters S and T as given from Figure 2 in [11], considering a chiral fourth generation quark t′. The upper limits at 68% (darker blue) and 95% (lighter blue) CL on the t′ mass are depicted as a function of the CKM mixing |Vt′b|, freely varied over [300,800] GeV/c2 and [0,0.3], respectively. The solid (dashed) curve denotes the upper bound mt′(|Vt′b|) under which t′t¯′ and t′q¯ (t′¯q) bound states can form, as defined by relation (6).As far as nonchiral fermions are concerned, the electroweak precision constraints provide upper bounds on the mixing parameters, but are not as restrictive given that their effects decouple in the limit of large masses [21]. Although vector-like quarks are usually considered to couple dominantly to the heaviest third quark family (see [12, 22] for some exceptions), they are allowed, in principle, to mix with all up- (down-) type quarks [14].In any case,|κQq|≲10-2 is used as a rule of thumb in most of the allowed vector-like representations [14]. Furthermore, corrections to the CKM matrix elements are required to be small for SU(2)L singlets [23], doublets, and triplets [24]. In this framework, the hadronic production of new bound states involving sequential or nonchiral quarks is a conceivable possibility at the LHC. Assuming that (6) is satisfied due to small Q-q couplings, their lifetime can be longer than the orbital period of the corresponding QQ¯ (Qq¯) bound state, allowing to build up new quarkonium resonances, open-flavour mesons, and baryons [20].If their formation is governed by the same strong interactions that are responsible for the existence of the ordinary hadronic spectroscopy, the associated mass spectrum can be determined from the known properties of low-energy hadron physics [25]. Considering a simplified model of QCD-like hadrons, new heavy quarks would form QQ¯ bound states (hereafter defined as ηQ) from induced, strongly attractive potentials. Conjointly, quarkonium formation can arise from Higgs boson exchange, as the corrections from Yukawa-type forces cannot be neglected for large quark masses [26, 27]. In [28], Hung detailed the numerical resolution of a general nonrelativistic Higgs exchange potential, with interesting consequences on the associated scalar spectrum. More recently, Enkhbat et al. gave a preliminary discussion of the fourth generation Yukawa bound states phenomenology at the LHC, considering an ultraheavy degenerate new quark family [29]. Although binding energies of Yukawa origin certainly provide a more consistent framework for a dedicated analysis, we restrict the present discussion to the context of Coulomb-like potential models and assume that the short-distance behaviour of the quark-antiquark potential dominates (we refer the reader to [30, 31] for previous and more exhaustive phenomenological reviews). In this context, the effects of the quarkonium wave function must be taken into account as we consider the corresponding production and decay modes. In perturbative QCD, the dominant contribution to ηQ production cross-section in pp collisions reads [30]
(7)σ(pp⟶gg⟶ηQ)=π2τ8mη3Γ(ηQ⟶gg)[τdLdτ]gg,(8)Γ(ηQ⟶gg)=8αs2(mη2)3mη2|RS(0)|2,
where mη=2mQ, and [τdL/dτ]gg denotes the gluon luminosity with τ=mη2/s. The full NLO partonic cross-sections for the gg, qg, and qq¯ initiated reactions are provided in [32] and have been updated in [33]. Interestingly, they are all proportional to the square of the QQ¯ radial wave function RS at the origin, which might drive the heavy quarkonium production cross-section σ(gg→ηQ) to be substantially smaller than σ(gg→QQ¯). Indeed, the quarkonium production rate typically amounts to a few percent of the pair production cross-section [32–34].Depending on whether the decay rate of the constituent particles,ΓQ, is larger or smaller than the bound state annihilation rate, two different cases that can be considered as the bound states can either undergo single quark decays or annihilate hadronically. As detailed in [35], the strength of the annihilation decay signal is enhanced when the intrinsic width of the heavy quark Q is decreased. Should its rate be of the same order of magnitude as the binding energy of the ηQ state, the quarkonium can be broad and display little evidence of any resonance behaviour over the continuum. On the other hand, if its 2-body decays are suppressed due to kinematical suppression or small couplings as discussed in Section 2.2, Q can only decay through off-shell intermediate states. In such a case, the signature would correspond to an annihilation signal, namely, a hadronic final state markedly distinguishable from pair production and decay.We now give a brief and qualitative description of some of the expectedηQ hadronic signatures at the LHC, considering either a chiral or vector-like heavy quark Q. While similar results hold for all S- and P-wave quarkonium states, we limit our analysis to the case of a stable, neutral JPC=0-+ pseudoscalar state ηQ. In particular, an1S0 quarkonium state formed from down-type fourth generation b′ quarks is known to provide a possible candidate if the condition (6) is satisfied [33]. The latter is mainly produced through the gluon fusion process (8) with a production cross-section two orders of magnitude larger than the JPC=1-- vector state [33]. If b′ is the lightest fourth generation quark, a ηb′ bound state can decay either through q-b′ family mixing with q=u,c,t, or to boson pairs. Whether the single quark decays would compete with the ηb′ hadronic modes depends on the heavy quark masses and mixings. If mb′>1 TeV/c2, a fourth-generation down-type quark dominantly decays into fermion-antifermion pairs if |Vqb′|≳10-2. For mb′<1 TeV/c2, the ηb′ total width lies below 1 GeV/c2, and the quarkonium state decays proceed dominantly via the annihilation diagrams depicted in Figure 3 [33]. Decays to gauge boson pairs, fermion-antifermion pairs, gauge-Higgs, and Higgs boson pairs can occur through γ,Z, and H exchange in the s-channel or proceed via quark mixing in the t- and u-channels, if allowed. Decays to Higgs boson pairs (ηQ→HH) are forbidden by CP conservation in the case of JPC=0-+ pseudoscalars. While the processes involving charged currents are suppressed in the case of small quark couplings, the ηb′→WW mode is allowed at loop-level via the exchange of the SU(2)L partner t′, even if |Vtb′|≃0. Decays to WW pairs are not allowed for ηT,B bound states involving singlet vector-like quarks.First-order diagrams for theηQ→p1p2 decays, where p1 and p2 can be a gauge boson pair, a fermion-antifermion pair, a gauge-Higgs boson pair, or a pair of Higgs bosons. The strong decay mode ηQ→gg is also allowed but not shown above. The (a), (b), and (c) diagrams denote decays via γ, Z, or H exchange in the s-channel, while (d) and (e) correspond to the p1p2 decays involving quark exchange in the t- and u-channels. Figure from [30].
(a)
(b)
(c)
(d)
(e)An exhaustive study of the production and subsequent decays of theS- and P-wave quarkonium states of a heavy chiral quark is available in [30, 36], where all the allowed production mechanisms and decay patterns have been thoroughly investigated. As an illustration, we compare in Figure 4 the branching ratios of the allowed decay modes for quarkonium states formed by a chiral fourth generation quark b′ and a vector-like isosinglet B quark [36]. The quarkonium annihilation decay rates depend markedly on the mass scale.Branching ratios as given in [36] for ηQ quarkonium, considering a chiral b′ quark (a) and a down-type singlet vector-like B quark (b), with mH=125 GeV/c2. |Vtb′|≃1 is assumed in the fourth generation case.
(a)
(b)In the chiral case, the main decay modes areZH, gg, WW, ZZ, Zγ, and γγ, respectively. The fermionic decays ηb′→qq¯ are also allowed if |Vqb′|≠0. Substantial decay rates to qq¯ fermion pairs are allowed in case of large Q-q couplings. For pseudoscalar ηb′ masses lighter than 450 GeV/c2, the dominant annihilation mode proceeds via the strong interaction, that is, ηb′→gg. For larger masses, the width into two gluons is decreasing, whereas the annihilation mode into ZH pairs takes over. This follows from the fact that the decay rates into Higgs and longitudinal gauge bosons are enhanced by the large Yukawa coupling of the heavy b′ quark, so that the ZH branching ratio becomes sizeable for increasing masses. The latter channel then leads to a salient signature for heavy fourth generation bound states.While fourth generation pseudoscalar quarkonia decays allow for striking signals,TT¯ and BB¯ bound states of vector-like quark singlets can only decay to gg, γγ, γZ, and ZZ pairs, as the ηQ→ZH mode is absent [36]. This follows from the fact that the ZH mode proceeds through the axial part of the neutral current coupling, which in turn is proportional to the third component of the weak isospin. Only γγ, Zγ, and ZZ would then signal the visible hadronic modes with branching ratios of about 10-4. TT¯ and BB¯ bound states of SU(2)L vector-like singlet quarks are thus hardly observable at the LHC as the gluon mode dominates. On the other hand, T and B vector-like doublet (triplet) partners with nonvanishing weak isospin can have large branching ratios to ZH via γ, Z, or H exchange in the s-channel. The ηX and ηY quarkonia provide a noticeable exception, given that the exotic quarks X5/3 and Y-4/3 do not couple to neutral currents at tree level. Since they only mix with the other states via charged currents, the annihilation modes including γ, Z, and H exchange are thus forbidden for ηX and ηY, which then mainly decay to gg and WW bosons. Although it is suppressed below the percent level, the distinctive diphoton decay mode also provides a possible golden channel for discovery. If new stable, coloured particles are produced at the LHC, the resulting final states could indeed allow for resonant signatures including photons and leptons, from which their quantum numbers and masses can be identified independently of their decay modes. As the decay rate for such bound states might be small, model-independent searches for new resonances could be competitive with the direct searches [35].


### 2.5. Signatures for Heavy Stable Particles: Open-Flavour Mesons

WhileQQ¯ quarkonium resonances might be challenging to observe, the possibility for “open-flavour” hadrons Qq¯ (Q¯q) and Qqq (Q¯q¯q¯), with q being a light SM quark, also covers a broad spectrum of new heavy mesons and baryons to search for. Should they be produced at the LHC, such states can form by capturing a light quark partner from the vacuum. As they pass through the detector, they can transform into various slowly moving heavy states. Assuming that they hadronise with u, d, and s quarks after being produced, Table 5 lists the corresponding allowed stable mesons and baryons.Table 5
Q
q
- (Q-q) mesons and Qqq (Q-q-q-) (anti-)baryons involving Q=X5/3, T2/3, B-1/3, and Y-4/3 vector-like quarks. The states in bold font are hadrons of which yields are expected to be substantial at the LHC, as predicted in [25] for penetration depths between 0 and 3 meters. For simplicity, only the neutral and positively charged hadrons are displayed.
Charge
Mesons
Baryons
Q
=
0
T
u
-
,
T
-
u
,
B
d
-
,
B
-
d
,
B
s
-
,
B
-
s
T
d
d
,
T
d
s
,
T
s
s
,
B
u
d
,
B
-
u
-
d
-
,
B
u
s, B-u-s-,Yuu,Y-u-u-
Q
=
1
X
u
-
,
T
d
-
,
T
s
-
,
B
-
u
,
Y
-
d
,
Y
-
s
X
d
d
,
X
d
s
,
X
s
s
,
T
u
d
,
T
u
s
,
B
u
u, B-d-d-,B-d-s-,B-s-s-,Y-u-d-,Y-u-s-
Q
=
2
X
d
-
,
X
s
-
,
Y
-
u
X
u
d
,
X
u
s
,
T
u
u
,
Y
-
d
-
d
-
,
Y
-
d
-
s
-
,
Y
-
s
-
s
-
Q
=
3
—
X
u
uThe properties of new hypothetical fourth generation quarkonium and open-flavour mesons have been examined previously in [37] considering large-N expansion techniques. In [38], their masses and decay constants have been estimated from the experimental measurements of the ordinary bb¯ (cc¯) hadron masses, using the QCD sum rules. Interestingly, both approaches indicate that the corresponding spectra share numerous features with most of the scenarios predicting stable bound states of heavy quarks beyond the Standard Model, depending on their masses, charges, and angular momenta.Interestingly, the interactions of such open-flavour mesons with the material occur to be similar to those ofR-hadrons, that is, stable hadrons composed of a supersymmetric particle and at least one SM quark (see [39, 40] and references therein for related works). The 1/2 difference in spin with respect to stop hadrons has little impact on the search strategies, which depend mostly on the interactions between the lighter quarks and the matter of the detector. This is in substance the description given by the spectator model [41, 42], for which the allowed transitions of heavy and light quarks within a given bound state are known to be flavour and spin independent [43, 44]. In this context, the heavy quark Q can be considered as a source of kinetic energy of which sole function is to give mass and momentum to the underlying hadron. Such new heavy bound states then behave as rather passive objects, consisting of a noninteracting heavy component, accompanied by lighter constituents (u and d quarks, mostly). While the partons are being scattered within the detector, their cross-sections vary with their inverse mass in perturbative QCD [6, 25]. A heavy Qq¯ (Qq¯) bound state is then expected to suffer small energy losses when interacting with the material, given that only its light quark half is responsible for the hadronic interactions with the detector.We list in Table6 the various allowed 2→2 and 2→3 transitions involving meson-to-meson and meson-to-baryon conversions. In general, these interactions lead to simultaneous pion emission as the initial mesons convert into slowly moving baryons.Table 6
Meson-to-Meson and Meson-to-Baryon processes forQu- and Qd- mesonic states, where Q=X5/3,T2/3,B-1/3,Y-4/3 and q=u,d. Similar transitions can be obtained for the charge conjugated processes involving Q-q. Assuming that the cross-sections for all the transitions are of the same order of magnitude, Meson-to-Baryon processes to Qud states are more likely than those involving Quu and Qdd in the final state.
(a)
Meson-to-Meson2→2 processes
Charge exchange
Q
u
-
+
p
↔
Q
d
-
+
n
Elastic scattering
Q
q
-
+
p
/
n
→
Q
q
-
+
p
/
n
(b)
Meson-to-Baryon2→2 processes
Baryon exchange
Q
u
-
+
n
→
Q
u
d
+
π
-
Q
u
-
+
n
→
Q
d
d
+
π
0
Q
u
-
+
p
→
Q
u
u
+
π
-
Q
u
-
+
p
→
Q
u
d
+
π
0
Q
d
-
+
n
→
Q
u
d
+
π
0
Q
d
-
+
n
→
Q
d
d
+
π
+
Q
d
-
+
p
→
Q
u
u
+
π
0
Q
d
-
+
p
→
Q
u
d
+
π
+
(c)
Meson-to-Meson2→3 processes
Inelastic scattering
Q
u
-
+
n
→
Q
d
-
+
n
+
π
-
Q
u
-
+
n
→
Q
u
-
+
n
+
π
0
Q
u
-
+
n
→
Q
u
-
+
p
+
π
-
Q
u
-
+
p
→
Q
d
-
+
p
+
π
-
Q
u
-
+
p
→
Q
u
-
+
p
+
π
0
Q
u
-
+
p
→
Q
u
-
+
n
+
π
+
Q
u
-
+
p
→
Q
d
-
+
n
+
π
0
Q
d
-
+
n
→
Q
u
-
+
n
+
π
+
Q
d
-
+
n
→
Q
d
-
+
n
+
π
0
Q
d
-
+
n
→
Q
u
-
+
p
+
π
0
Q
d
-
+
n
→
Q
d
-
+
p
+
π
-
Q
d
-
+
p
→
Q
u
-
+
p
+
π
+
Q
d
-
+
p
→
Q
d
-
+
p
+
π
0
Q
d
-
+
p
→
Q
d
-
+
n
+
π
+
(d)
Meson-to-Baryon2→3 processes
Baryon exchange
Q
u
-
+
n
→
Q
u
d
+
π
-
+
π
0
Q
u
-
+
n
→
Q
d
d
+
π
-
+
π
+
Q
u
-
+
n
→
Q
d
d
+
π
0
+
π
0
Q
u
-
+
p
→
Q
u
u
+
π
-
+
π
0
Q
u
-
+
p
→
Q
u
d
+
π
-
+
π
+
Q
u
-
+
p
→
Q
u
d
+
π
0
+
π
0
Q
u
-
+
p
→
Q
d
d
+
π
+
+
π
0
Q
d
-
+
n
→
Q
u
d
+
π
-
+
π
+
Q
d
-
+
n
→
Q
u
d
+
π
0
+
π
0
Q
d
-
+
n
→
Q
d
d
+
π
0
+
π
+
Q
d
-
+
p
→
Q
u
u
+
π
-
+
π
+
Q
d
-
+
p
→
Q
u
u
+
π
0
+
π
0
Q
d
-
+
p
→
Q
u
d
+
π
0
+
π
+While they interact with the nuclear matter, most of the new states suffer multiple scattering and convert into baryons, allowing forQuu, Qud, and Qdd states, as new heavy mesons are kinematically favoured to increase their baryon number by emitting one or more pions [25]. Due to the lack of these in the material, the reverse reaction is known to be less favoured.Interestingly, new mesons which convert at the beginning of the scattering chain could generate tracks with “dashed-lines,” signalling possible baryonic or electric charge exchange. Indeed,Qq¯ and Q¯q bound states traversing a medium composed of light quarks likely flip their electric charge, frequently interchanging their parton constituents with those of the material nuclei. Eventually, they could leave the detector as heavy stable neutral particles and hence not be observable in muon detectors. The modeling of the nuclear interactions of new heavy hadrons traveling through matter, a survey of which is given in [25, 45], actually favours scenarios with significant charge suppression if heavy open-flavour mesons have a sizeable probability to transform into neutral particles while traversing the detector. This poses a serious challenge for the experimental searches, since the most conservative scenarios assume a complete charge suppression where 100% of the produced hadrons become neutral before reaching the muon detectors. Taking a specific example, we consider pair-produced T2/3 quarks hadronising into the heavy states Tq¯ and T¯q immediately after production. For convenience, we assume that a majority (≳90%) of the new heavy quarks initially hadronise into mesons, while a smaller amount (≲10%) form baryons [6]. As they interact farther in the detector, most of the Tq¯ mesons eventually exit the detector as Tu¯ and Td¯ mesonic states, or as Tud baryons. Baryon-to-meson conversions are allowed as well, as these states can also annihilate back into Tu¯ and Td¯ mesons, yet with a smaller rate. T¯q mesons of vector-like antiquarks, on the other hand, unlikely give rise to antibaryons, but can still flip their charge through the exchange of u and d quarks with the material. Possibly large flavour fractions for T¯u and T¯d can thus be expected in the detector, with a negligible amount of T¯q¯q¯ states. If such antibaryons were to form, they would quickly annihilate by baryon-to-meson interactions, producing pions [46].After travelling through the detectors a few nuclear lengths away from the production vertex, the novel hadrons transform in roughly 100% to the positively charged baryonTud, whereas one half of the T¯-states thus transform into T¯u mesons, and the other half in T¯d. The corresponding hadronic flavour decompositions, evaluated as a function of the penetration depth in the detector, can be read from Figure 5 of [25] for R-hadrons involving stable scalar top and antitop quarks. We notice that similar results are obtained when replacing stop t~ squarks by stable quarks T of equal mass. The Tud, Tu¯ (T¯u), and Td¯ (T¯d), states are then expected to retain the largest flavour composition fractions when considering penetration depths larger than 3 meters and are thus the most likely observable states, which we emphasised in bold font in Table 5.Except forYuu and Y¯u¯u¯, it is interesting to notice that none of the states listed in Table 5 is neutral if involving vector-like quarks with exotic charges. The corresponding bound states might give rise to possibly large fractions of slowly moving charged particles Xud and Yud with Q=+2e and -e, respectively, accompanied by a small amount of Xu¯, Xd¯, Yu¯, and Yd¯ mesons (and charge conjugates). These particles are all electrically charged as they go across the detector, emitting charged and neutral pions in the process together with losing small amounts of energy. Given the small expected interaction cross-section from the Meson-to-Baryon transitions, the processes shown in Table 6 leave small energy deposits in the calorimeters, of the order of O(1) GeV [6]. On the other hand, such stable hadrons would lead to observable tracks due to the ionisation energy losses, allowing for signatures similar to slow-moving muons with high transverse momentum. Searches for stable charged particles, if adapted to such a case, thus provide a promising strategy to rule out the possibility for novel exotic quarks with large lifetimes. Given these very specific signals, the aforementioned signatures can certainly be discriminated at the LHC.


## 2.1. Production

As we review the possibility for long-lived quarks in general, two different scenarios will be distinguished in this work. New chiral fermions, for which the left- and right-handed chiralities have different charges under the electroweak gauge group, gain their masses from the electroweak symmetry breaking mechanism. If heavy, a sequential fourth generation quark doublet(t′,b′)L can therefore induce large corrections to loop observables. Such a scenario is now severely constrained from the electroweak precision data and the Higgs search results. Given that chiral fermions couple to the Higgs boson with a strength proportional to their Yukawa couplings, they do not decouple from its production, as the corresponding rate should increase by a factor of about 9 due to the additional fermion loops occurring in gluon fusion. The unobserved enhancement thus now strongly disfavours a fourth generation Standard Model [3].Nonchiral quarks, on the other hand, still provide a viable extension to the SM and certainly require a dedicated discussion. Considering the introduction of new vector-like heavy partners mixing with the SM quarks, additional parameters are allowed fromSU(2)L×U(1)Y invariant Yukawa interactions and Dirac mass terms in the SM Lagrangian. The light quark couplings to the Higgs and electroweak bosons are consequently modified, along with the introduction of new couplings between the heavy and the SM quarks.Still, it is important to emphasise that such new states cannot have arbitrary quantum numbers, since they can only mix with the SM quarks through a limited number of gauge-invariant couplings. Classifying them intoSU(2)L multiplets, their Yukawa terms only allow for seven distinct possibilities, that is, the two singlets, the three doublets, and the two triplets displayed in Table 1.Table 1
Vector-like multiplets allowed to mix with the SM quarks through Yukawa couplings. The electric charge is the sum of the third component of isospinT3 and of the hypercharge Y. Adapted from [9].
Q
q
T
2
/
3
B
-
1
/
3
(
X
5
/
3
T
2
/
3
)
(
T
2
/
3
B
-
1
/
3
)
(
B
-
1
/
3
Y
-
4
/
3
)
(
X
5
/
3
T
2
/
3
B
-
1
/
3
)
(
T
2
/
3
B
-
1
/
3
Y
-
4
/
3
)
T
3
0
0
1
/
2
1
/
2
1
/
2
1
1
Y
2
/
3
-
1
/
3
7
/
6
1
/
6
-
5
/
6
2
/
3
-
1
/
3The production of such new heavy quarks, either chiral or vector-like, is usually assumed to proceed dominantly at the LHC through gluon fusion,gg→QQ¯. Nevertheless, depending on the model at hand, electroweak single production can also provide an alternative mechanism as it is not affected by the large phase-space suppression from which pair production suffers. In particular, new heavy quarks can be produced singly in flavour-changing processes via the electroweak interaction through qiqj(-)→V*→qkQ, where V=W,Z. A comparison between the dominant production modes can be found in [10, 11] for fourth generation and in [12] for vector-like quarks. Assuming order unity couplings, benchmark cross-sections have been evaluated for various mass scenarios and reproduced in Table 2.Table 2
Electroweak single production cross-sections (in pb) as a function of the massmQ of the heavy quark, assuming order unity couplings as computed in [12]. In contrast, strong pair production pp→QQ- yields 3.61 (0.63) fb for mQ=900 and 1800 GeV/c2 at s=7 and 14 TeV.
Channel
s
=
7 TeV
s
=
14 TeV
m
Q
=
900 GeV/c2
m
Q
=
1800 GeV/c2
p
p
→
W
*
→
j
X
1.4
0.36
p
p
→
W
*
→
j
X
-
0.037
0.0092
p
p
→
W
*
→
j
T
0.61
0.16
p
p
→
W
*
→
j
T
-
0.052
0.013
p
p
→
Z
*
→
j
T
0.43
0.11
p
p
→
Z
*
→
j
T
-
0.025
0.0064
p
p
→
W
*
→
j
B
0.69
0.18
p
p
→
W
*
→
j
B
-
0.089
0.022
p
p
→
Z
*
→
j
B
0.18
0.047
p
p
→
Z
*
→
j
B
-
0.034
0.0088
p
p
→
W
*
→
j
Y
0.29
0.074
p
p
→
W
*
→
j
Y
-
0.12
0.031As described in [12], the leading contributions to vector-like T and B single production arise from ud→W*→jT, du→W*→jB, qd→Z*→jB and qu→Z*→jT, where q denotes a valence quark parton and j a generic light quark jet. Unlike QCD pair production, the above processes, however, scale with the Q-q quark coupling squared and can be strongly suppressed if the associated mixings are negligibly small. On the other hand, Figures  2–4 in [12] indicate that σ(pp→QQ¯) falls off faster than the single production cross-sections as soon as mQ reaches the TeV scale, while the electroweak production channels displayed in Table 2 are enhanced by a factor O(mQ2/mV2), originating from the longitudinal polarisation of the gauge boson V. Furthermore, the relative rates of the jQ and jQ¯ channels strongly depend on the valence quark density in the initial state. Given the difference in the PDFs of valence and sea quarks in the initial states, new heavy quarks are produced singly at a much higher rate than antiquarks for increasing mQ. While the u quark contributes more than the d quark, the neutral current contributions are weaker than the charged current ones. Finally, we underline the absence of the tree-level processes pp→Z*→jQ for vector-like quarks involving exotic charges, given that they only interact with other states via tree-level charged currents [12]. Although most of the above statements are model dependent, we conclude that the large running energy of the LHC makes single electroweak production a promising discovery channel when considering new heavy quark searches with mQ≳1 TeV/c2.


## 2.2. Decay

In the search of new fermions, the associated decay modes require careful attention. The exact partial width for a new sequential heavy quarkQ decaying on-shell to a light quark q through a charged current can be written as
(1)Γ(Q⟶qW)=GFmQ38π2|κQq|2f2(mqmQ,mWmQ),(2)f2(α,β)=[(1-α2)2+β2(1+α2)-2β4]×[1-(α+β)2][1-(α-β)2],
where κQq denotes the generic Q-q quark coupling, equal to the CKM mixing matrix element VQq when considering a new sequential family of quarks. Interestingly, the width (1) holds for both chiral and vector-like singlet quarks, given that all heavy-to-light charged current quark decays occur to be pure V-A processes with identical rates [13]. If mQ≫mq, the above partial width can be shown to reach the asymptotic form
(3)Γ(Q⟶qW)≃170.5MeV×|κQq|2×mQ3mW3.
Considering a new heavy quark Q decaying exclusively via the above 2-body decay mode, we evaluate in Table 3 the corresponding decay lengths, with cτ≃2×10-10 GeV/Γ(GeV) μm and |κQq|=1. With increasing mass, the probability for a new heavy quark Q to decay weakly thus becomes more and more important as its lifetime decreases. However, small couplings to the lighter SM quarks strongly affect this statement as the corresponding widths are suppressed by |κQq|2. As we will see in Section 2.4, this could lead to a very different phenomenology as the quarks Q, if very heavy, would form bound states and possibly decay hadronically.Table 3
Partial widths and decay lengths for a new heavy quarkQ decaying to light quarks via on-shell charged currents Q→qW, with mq=0 and 173.2 GeV/c2 and |κQq|=1, as computed from (1).
m
Q (GeV/c2)
Γ (GeV)
c
τ
×
1
0
12 (μm)
300
8.73 (2.39)
22.91 (83.85)
400
20.90 (11.14)
9.57 (17.96)
500
40.94 (27.86)
4.89 (7.18)
600
70.81 (54.52)
2.82 (3.67)
700
112.50 (93.06)
1.78 (2.15)
800
167.97 (147.43)
1.19 (1.38)
900
239.18 (213.57)
83.62 (9.36) ×10-1
1000
328.12 (299.46)
6.10 (6.68) ×10-1
1100
436.74 (405.05)
4.58 (4.94) ×10-1
1200
567.02 (532.31)
3.53 (3.76) ×10-1
1300
720.93 (683.21)
2.77 (2.92) ×10-1
1400
900.43 (859.72)
2.22 (2.33) ×10-1
1500
1107.5 (1063.79)
1.81 (1.88) ×10-1
1600
1344.1 (1297.4)
1.49 (1.54) ×10-1
1700
1612.2 (1562.5)
1.24 (1.28) ×10-1
1800
1913.8 (1861.1)
1.05 (1.07) ×10-1Several conclusions can already be drawn at this stage. As far as a heavy fourth quark family(t′,b′)L is concerned, we highlight that the ratios of the charged current decay rates to the lighter families only depend on the off-diagonal mixing elements. Assuming that the extended 4×4 CKM matrix is unitary, the CKM mixing elements |VQq| imply nonunity branching ratios for fourth generation quarks in general.The possibility for tree-level Flavour Changing Neutral Currents leads to similar conclusions for models involving new vector-like quarks interacting with the SM families. Depending on their Yukawa couplings, they can undergo tree-level transitions to lighter quarks andW, Z, and Higgs bosons. Considering vector-like quarks Q=T,B with mQ≫mq, the neutral currents Q→qZ and Q→qH can occur with a rate of the same magnitude as given by (3), if they are allowed to mix with the light SM quarks q. If no assumption is made on the hierarchy of the Q-q coupling strengths, no exclusive decay mode to the first, second, or third family should be preferred.For nonchiral quarks decaying toWq, Zq, and Hq, long lifetimes can be expected in the |κQq|≃0 limit. Additionally, it is interesting to notice that Γ(Q→qZ)/Γ(Q→qH)≃1 holds in the large mass limit mQ≫mq for most scenarios, while the charged-to-neutral current ratios Γ(Q→qW)/Γ(Q→q(Z,H)) approach either 2 or 0 depending on the vector-like representation [14].If more than one new heavy quark is considered, and if they are allowed to decay into each other, the partial width (3) becomes inaccurate due to the possible competition with the heavy-to-heavy transitions Q1→Q2V(*), where V=W,Z,H, depending on the model at hand. Indeed, if the decay rates of the daughter particles are significant with respect to the mass of the parent, the decay of an unstable particle is allowed through threshold effects, even if kinematically forbidden. For both chiral and vector-like quarks, real and virtual V emission near threshold proceed with the rate
(4)Γ(Q1→Q2V(*))=GF2mQ15192π3|κQ1Q2|2f3(mQ12mV2,mQ22mV2,ΓV2mV2),f3(α,β,γ)=2∫0(1-β)2dx‍×[(1-β)2+x(1+β)-2x2]λ(1,x,β)[(1-αx)2+γ2],
where λ(a,b,c)=a2+b2+c2-2ab-2bc-2ac [4]. While the Q→qV rates essentially depend on the Q-q quark couplings, the heavy quark mass splittings |mQ1-mQ2| drive the strength of the heavy-to-heavy transitions [11, 15]. If the quarks Q1 and Q2 belong to the same weak multiplet with a large mass splitting, the Q1→Q2V decay could be very rapid and dominate the other decay modes if there is no |κQ1Q2| suppression. On the other hand, the lightest partner can only decay into an SM family quark and is thus plausibly long-lived for small mixings. Should this be the case, we will see in Sections 2.4 and 2.5 that the latter could hadronize, leading to a different decay phenomenology. This occurs specifically for SU(2)L vector-like singlets, which free quark decays can only proceed via nonzero mixing into the lighter generations.As far as the doublet and triplet representations are concerned, new vector-like partners should be very close in mass if one assumes that isospin conservation is respected. Although it may be broken by higher order corrections, such a degeneracy has been shown in [16] to require a mass difference between 70 and 110MeV for nonchiral doublets, allowing for nanoseconds lifetimes. If their couplings to the lighter quarks are negligible, the partial widths for X5/3, T2/3, B-1/3, and Y-4/3 quarks are thus plausibly small.Summarising these results, the range of the observable couplings to which long-lived particle searches are sensitive at the LHC can be evaluated. From the approximation (3), the partial decay length for a vector-like quark reads
(5)λ≃2cm×(10-8|κQq|)2(500GeV/c2mQ)3.
Assuming mQ≲1 TeV/c2 and an experimental resolution of at least 25 μm, the heavy quark decays are prompt only if their associated couplings are above the 10-7 level. If |κQq|≲10-7, the corresponding decay lengths could be observed over a few tens of centimeters. In the limit of such small mixings with the SM fermions, it is thus natural to motivate searches for displaced vertices.


## 2.3. Classification of Signatures

We now detail the possible situations giving rise to the observable signatures of long-lived quarks at the LHC.The classification given in Table4 summarises the long-lived quark scenarios corresponding to the small mixing scenarios treated in this work. We emphasise that if mixing with all SM quark families is permitted, direct searches should aim at being sensitive to all light quarks q in the observed final states. Given that their interactions are allowed through arbitrary Yukawa couplings, the searches for long-lived vector-like quarks should be as inclusive as possible in order to cover the full spectrum of possibilities at the LHC. Such a program is important to set new constraints on the Q→qV(*) decay modes, independently of any assumptions on the Q-q mixing sector. The relative importance of the Q→(u,c,d,s)V channels is most relevant compared to the Q→(t,b)V charged currents, as no exclusive decay mode should be preferred for V=W,Z,H.Table 4
Possible decay signatures for new long-lived quarks (mQ=1 TeV/c2), as discussed in the text.
(i)
(ii)
(iii)
Short-lived
Intermediate
Long-lived
|
κ
Q
q
|
≳

1
0
-
7
≲

1
0
-
7
≳

1
0
-
9
≲

1
0
-
9
Γ (GeV)
≳

1
0
-
12
≲

1
0
-
12
≳

1
0
-
16
≲

1
0
-
16
τ (s)
≲

1
0
-
13
≳

1
0
-
13
≲

1
0
-
9
≳

1
0
-
9As illustrated in Figure1, three cases for new long-lived quarks can be distinguished depending on their mixing parameters with the SM fermions. (i)
Firstly, if all couplings with the known SM fermions are larger than the10-7 level, all heavy-to-light decays Q→qV are prompt, over distances smaller than 25 μm. This defines the short-lived scenario (i), with no experimental difference compared to the direct searches.(ii)
The second scenario (ii) arises for intermediate decay lengths ranging between a few microns and centimetric distances. If there is at least one light quarkq such that 10-7≲|κQq|≲10-9, displaced events could be observed, yet with a partial width suppressed by |κQq|2. Consequently, if all the couplings but one are below the 10-7 level, a single exclusive channel would lead to a prompt decay signature corresponding to a short-lived heavy quark. A possible exception, though, is the case for new heavy multiplets with sizeable mass splittings, as allowed in extra generation models and possible extensions [11, 17–19]. If mQ1≳mQ2+mV, the heaviest quark Q1 can be short-lived and decay semiweakly to Q2V(*), while the lightest partner is likely stable if all its decay modes suffer severe suppression. If mQ1≃mQ2, on the other hand, all heavy-to-heavy transitions are suppressed and both quarks could be long-lived on detector scales.(iii)
For particles with decay lengths larger than the detector dimensions lies the long-lived region (iii), the details of which will be discussed in Section2.4. If all heavy quark couplings with the SM fermions are below the 10-9 level, the stable case becomes a relevant scenario, possibly in conjunction with displaced events.Figure 1
Lifetimeτ in seconds (γβ≃1) for a sequential heavy quark Q decaying to the lighter families by emitting a real W boson, for mQ=300 GeV/c2 (in blue) and mQ=1800 GeV/c2 (in red). The three regions of interest refer to the conventions given in Table 4.We emphasise that the scenarios for which prompt decays and displaced vertices could be observed only call for one of the allowed decay modes to fulfill the associated conditions. The stable scenario, however, requires all decay lengths and partial widths to lie in the ranges given in Table4. As we will see in Section 2.4, if all Q-q quark couplings to the SM families are smaller than 10-2, such new heavy fermions could hadronise. As a result, the annihilation decays and hadronic transitions between the formed bound states would dominate, while the single quark decay events might be suppressed.


## 2.4. Signatures for Heavy Stable Particles: Heavy Quarkonia

Previously, we have discussed how displaced events could occur in the single quark decay transitions in the case of smallQ-q couplings. In this section, we describe the possibility for new heavy quarks to bind into hadrons if their lifetime is large enough, providing us with an interesting alternative for new signal searches at the LHC.While new coloured particles with masses heavier than the top quark are usually assumed to decay as free particles, their possible formation into baryons and mesons is an important issue to consider in the case of small couplings with the SM fermions.Should the partial widthsQ→q(W,Z,H) be suppressed, the quark Q might be stable and hadronise. Assuming that the binding force in a fourth generation bound state is of Coulombic nature, the seminal condition
(6)mQ<125GeV(100GeV)×|VQq|-2/3
was derived in [20] for QQ¯ (Qq¯) formation. According to (6), new heavy quarks form hadronic states if their mixing with known quarks is sufficiently small. Mixings roughly smaller than 0.2 are typically required for 300 GeV/c2 new quarks to form bound states, which remains perfectly consistent with the current bounds from the electroweak precision observables. As shown on Figure 2 for a new up-like fourth generation quark t′, the electroweak precision fits typically restrict the mixing matrix element |Vt′b|≃|Vtb′| to be smaller than 10-1 for mt′>700 GeV/c2. If lighter, the t′ total width can be smaller than 10MeV [11], thus allowing new chiral quarks to form bound states with nanosecond lifetimes.Figure 2
Reproduction of theχ2 fit of the oblique parameters S and T as given from Figure 2 in [11], considering a chiral fourth generation quark t′. The upper limits at 68% (darker blue) and 95% (lighter blue) CL on the t′ mass are depicted as a function of the CKM mixing |Vt′b|, freely varied over [300,800] GeV/c2 and [0,0.3], respectively. The solid (dashed) curve denotes the upper bound mt′(|Vt′b|) under which t′t¯′ and t′q¯ (t′¯q) bound states can form, as defined by relation (6).As far as nonchiral fermions are concerned, the electroweak precision constraints provide upper bounds on the mixing parameters, but are not as restrictive given that their effects decouple in the limit of large masses [21]. Although vector-like quarks are usually considered to couple dominantly to the heaviest third quark family (see [12, 22] for some exceptions), they are allowed, in principle, to mix with all up- (down-) type quarks [14].In any case,|κQq|≲10-2 is used as a rule of thumb in most of the allowed vector-like representations [14]. Furthermore, corrections to the CKM matrix elements are required to be small for SU(2)L singlets [23], doublets, and triplets [24]. In this framework, the hadronic production of new bound states involving sequential or nonchiral quarks is a conceivable possibility at the LHC. Assuming that (6) is satisfied due to small Q-q couplings, their lifetime can be longer than the orbital period of the corresponding QQ¯ (Qq¯) bound state, allowing to build up new quarkonium resonances, open-flavour mesons, and baryons [20].If their formation is governed by the same strong interactions that are responsible for the existence of the ordinary hadronic spectroscopy, the associated mass spectrum can be determined from the known properties of low-energy hadron physics [25]. Considering a simplified model of QCD-like hadrons, new heavy quarks would form QQ¯ bound states (hereafter defined as ηQ) from induced, strongly attractive potentials. Conjointly, quarkonium formation can arise from Higgs boson exchange, as the corrections from Yukawa-type forces cannot be neglected for large quark masses [26, 27]. In [28], Hung detailed the numerical resolution of a general nonrelativistic Higgs exchange potential, with interesting consequences on the associated scalar spectrum. More recently, Enkhbat et al. gave a preliminary discussion of the fourth generation Yukawa bound states phenomenology at the LHC, considering an ultraheavy degenerate new quark family [29]. Although binding energies of Yukawa origin certainly provide a more consistent framework for a dedicated analysis, we restrict the present discussion to the context of Coulomb-like potential models and assume that the short-distance behaviour of the quark-antiquark potential dominates (we refer the reader to [30, 31] for previous and more exhaustive phenomenological reviews). In this context, the effects of the quarkonium wave function must be taken into account as we consider the corresponding production and decay modes. In perturbative QCD, the dominant contribution to ηQ production cross-section in pp collisions reads [30]
(7)σ(pp⟶gg⟶ηQ)=π2τ8mη3Γ(ηQ⟶gg)[τdLdτ]gg,(8)Γ(ηQ⟶gg)=8αs2(mη2)3mη2|RS(0)|2,
where mη=2mQ, and [τdL/dτ]gg denotes the gluon luminosity with τ=mη2/s. The full NLO partonic cross-sections for the gg, qg, and qq¯ initiated reactions are provided in [32] and have been updated in [33]. Interestingly, they are all proportional to the square of the QQ¯ radial wave function RS at the origin, which might drive the heavy quarkonium production cross-section σ(gg→ηQ) to be substantially smaller than σ(gg→QQ¯). Indeed, the quarkonium production rate typically amounts to a few percent of the pair production cross-section [32–34].Depending on whether the decay rate of the constituent particles,ΓQ, is larger or smaller than the bound state annihilation rate, two different cases that can be considered as the bound states can either undergo single quark decays or annihilate hadronically. As detailed in [35], the strength of the annihilation decay signal is enhanced when the intrinsic width of the heavy quark Q is decreased. Should its rate be of the same order of magnitude as the binding energy of the ηQ state, the quarkonium can be broad and display little evidence of any resonance behaviour over the continuum. On the other hand, if its 2-body decays are suppressed due to kinematical suppression or small couplings as discussed in Section 2.2, Q can only decay through off-shell intermediate states. In such a case, the signature would correspond to an annihilation signal, namely, a hadronic final state markedly distinguishable from pair production and decay.We now give a brief and qualitative description of some of the expectedηQ hadronic signatures at the LHC, considering either a chiral or vector-like heavy quark Q. While similar results hold for all S- and P-wave quarkonium states, we limit our analysis to the case of a stable, neutral JPC=0-+ pseudoscalar state ηQ. In particular, an1S0 quarkonium state formed from down-type fourth generation b′ quarks is known to provide a possible candidate if the condition (6) is satisfied [33]. The latter is mainly produced through the gluon fusion process (8) with a production cross-section two orders of magnitude larger than the JPC=1-- vector state [33]. If b′ is the lightest fourth generation quark, a ηb′ bound state can decay either through q-b′ family mixing with q=u,c,t, or to boson pairs. Whether the single quark decays would compete with the ηb′ hadronic modes depends on the heavy quark masses and mixings. If mb′>1 TeV/c2, a fourth-generation down-type quark dominantly decays into fermion-antifermion pairs if |Vqb′|≳10-2. For mb′<1 TeV/c2, the ηb′ total width lies below 1 GeV/c2, and the quarkonium state decays proceed dominantly via the annihilation diagrams depicted in Figure 3 [33]. Decays to gauge boson pairs, fermion-antifermion pairs, gauge-Higgs, and Higgs boson pairs can occur through γ,Z, and H exchange in the s-channel or proceed via quark mixing in the t- and u-channels, if allowed. Decays to Higgs boson pairs (ηQ→HH) are forbidden by CP conservation in the case of JPC=0-+ pseudoscalars. While the processes involving charged currents are suppressed in the case of small quark couplings, the ηb′→WW mode is allowed at loop-level via the exchange of the SU(2)L partner t′, even if |Vtb′|≃0. Decays to WW pairs are not allowed for ηT,B bound states involving singlet vector-like quarks.First-order diagrams for theηQ→p1p2 decays, where p1 and p2 can be a gauge boson pair, a fermion-antifermion pair, a gauge-Higgs boson pair, or a pair of Higgs bosons. The strong decay mode ηQ→gg is also allowed but not shown above. The (a), (b), and (c) diagrams denote decays via γ, Z, or H exchange in the s-channel, while (d) and (e) correspond to the p1p2 decays involving quark exchange in the t- and u-channels. Figure from [30].
(a)
(b)
(c)
(d)
(e)An exhaustive study of the production and subsequent decays of theS- and P-wave quarkonium states of a heavy chiral quark is available in [30, 36], where all the allowed production mechanisms and decay patterns have been thoroughly investigated. As an illustration, we compare in Figure 4 the branching ratios of the allowed decay modes for quarkonium states formed by a chiral fourth generation quark b′ and a vector-like isosinglet B quark [36]. The quarkonium annihilation decay rates depend markedly on the mass scale.Branching ratios as given in [36] for ηQ quarkonium, considering a chiral b′ quark (a) and a down-type singlet vector-like B quark (b), with mH=125 GeV/c2. |Vtb′|≃1 is assumed in the fourth generation case.
(a)
(b)In the chiral case, the main decay modes areZH, gg, WW, ZZ, Zγ, and γγ, respectively. The fermionic decays ηb′→qq¯ are also allowed if |Vqb′|≠0. Substantial decay rates to qq¯ fermion pairs are allowed in case of large Q-q couplings. For pseudoscalar ηb′ masses lighter than 450 GeV/c2, the dominant annihilation mode proceeds via the strong interaction, that is, ηb′→gg. For larger masses, the width into two gluons is decreasing, whereas the annihilation mode into ZH pairs takes over. This follows from the fact that the decay rates into Higgs and longitudinal gauge bosons are enhanced by the large Yukawa coupling of the heavy b′ quark, so that the ZH branching ratio becomes sizeable for increasing masses. The latter channel then leads to a salient signature for heavy fourth generation bound states.While fourth generation pseudoscalar quarkonia decays allow for striking signals,TT¯ and BB¯ bound states of vector-like quark singlets can only decay to gg, γγ, γZ, and ZZ pairs, as the ηQ→ZH mode is absent [36]. This follows from the fact that the ZH mode proceeds through the axial part of the neutral current coupling, which in turn is proportional to the third component of the weak isospin. Only γγ, Zγ, and ZZ would then signal the visible hadronic modes with branching ratios of about 10-4. TT¯ and BB¯ bound states of SU(2)L vector-like singlet quarks are thus hardly observable at the LHC as the gluon mode dominates. On the other hand, T and B vector-like doublet (triplet) partners with nonvanishing weak isospin can have large branching ratios to ZH via γ, Z, or H exchange in the s-channel. The ηX and ηY quarkonia provide a noticeable exception, given that the exotic quarks X5/3 and Y-4/3 do not couple to neutral currents at tree level. Since they only mix with the other states via charged currents, the annihilation modes including γ, Z, and H exchange are thus forbidden for ηX and ηY, which then mainly decay to gg and WW bosons. Although it is suppressed below the percent level, the distinctive diphoton decay mode also provides a possible golden channel for discovery. If new stable, coloured particles are produced at the LHC, the resulting final states could indeed allow for resonant signatures including photons and leptons, from which their quantum numbers and masses can be identified independently of their decay modes. As the decay rate for such bound states might be small, model-independent searches for new resonances could be competitive with the direct searches [35].


## 2.5. Signatures for Heavy Stable Particles: Open-Flavour Mesons

WhileQQ¯ quarkonium resonances might be challenging to observe, the possibility for “open-flavour” hadrons Qq¯ (Q¯q) and Qqq (Q¯q¯q¯), with q being a light SM quark, also covers a broad spectrum of new heavy mesons and baryons to search for. Should they be produced at the LHC, such states can form by capturing a light quark partner from the vacuum. As they pass through the detector, they can transform into various slowly moving heavy states. Assuming that they hadronise with u, d, and s quarks after being produced, Table 5 lists the corresponding allowed stable mesons and baryons.Table 5
Q
q
- (Q-q) mesons and Qqq (Q-q-q-) (anti-)baryons involving Q=X5/3, T2/3, B-1/3, and Y-4/3 vector-like quarks. The states in bold font are hadrons of which yields are expected to be substantial at the LHC, as predicted in [25] for penetration depths between 0 and 3 meters. For simplicity, only the neutral and positively charged hadrons are displayed.
Charge
Mesons
Baryons
Q
=
0
T
u
-
,
T
-
u
,
B
d
-
,
B
-
d
,
B
s
-
,
B
-
s
T
d
d
,
T
d
s
,
T
s
s
,
B
u
d
,
B
-
u
-
d
-
,
B
u
s, B-u-s-,Yuu,Y-u-u-
Q
=
1
X
u
-
,
T
d
-
,
T
s
-
,
B
-
u
,
Y
-
d
,
Y
-
s
X
d
d
,
X
d
s
,
X
s
s
,
T
u
d
,
T
u
s
,
B
u
u, B-d-d-,B-d-s-,B-s-s-,Y-u-d-,Y-u-s-
Q
=
2
X
d
-
,
X
s
-
,
Y
-
u
X
u
d
,
X
u
s
,
T
u
u
,
Y
-
d
-
d
-
,
Y
-
d
-
s
-
,
Y
-
s
-
s
-
Q
=
3
—
X
u
uThe properties of new hypothetical fourth generation quarkonium and open-flavour mesons have been examined previously in [37] considering large-N expansion techniques. In [38], their masses and decay constants have been estimated from the experimental measurements of the ordinary bb¯ (cc¯) hadron masses, using the QCD sum rules. Interestingly, both approaches indicate that the corresponding spectra share numerous features with most of the scenarios predicting stable bound states of heavy quarks beyond the Standard Model, depending on their masses, charges, and angular momenta.Interestingly, the interactions of such open-flavour mesons with the material occur to be similar to those ofR-hadrons, that is, stable hadrons composed of a supersymmetric particle and at least one SM quark (see [39, 40] and references therein for related works). The 1/2 difference in spin with respect to stop hadrons has little impact on the search strategies, which depend mostly on the interactions between the lighter quarks and the matter of the detector. This is in substance the description given by the spectator model [41, 42], for which the allowed transitions of heavy and light quarks within a given bound state are known to be flavour and spin independent [43, 44]. In this context, the heavy quark Q can be considered as a source of kinetic energy of which sole function is to give mass and momentum to the underlying hadron. Such new heavy bound states then behave as rather passive objects, consisting of a noninteracting heavy component, accompanied by lighter constituents (u and d quarks, mostly). While the partons are being scattered within the detector, their cross-sections vary with their inverse mass in perturbative QCD [6, 25]. A heavy Qq¯ (Qq¯) bound state is then expected to suffer small energy losses when interacting with the material, given that only its light quark half is responsible for the hadronic interactions with the detector.We list in Table6 the various allowed 2→2 and 2→3 transitions involving meson-to-meson and meson-to-baryon conversions. In general, these interactions lead to simultaneous pion emission as the initial mesons convert into slowly moving baryons.Table 6
Meson-to-Meson and Meson-to-Baryon processes forQu- and Qd- mesonic states, where Q=X5/3,T2/3,B-1/3,Y-4/3 and q=u,d. Similar transitions can be obtained for the charge conjugated processes involving Q-q. Assuming that the cross-sections for all the transitions are of the same order of magnitude, Meson-to-Baryon processes to Qud states are more likely than those involving Quu and Qdd in the final state.
(a)
Meson-to-Meson2→2 processes
Charge exchange
Q
u
-
+
p
↔
Q
d
-
+
n
Elastic scattering
Q
q
-
+
p
/
n
→
Q
q
-
+
p
/
n
(b)
Meson-to-Baryon2→2 processes
Baryon exchange
Q
u
-
+
n
→
Q
u
d
+
π
-
Q
u
-
+
n
→
Q
d
d
+
π
0
Q
u
-
+
p
→
Q
u
u
+
π
-
Q
u
-
+
p
→
Q
u
d
+
π
0
Q
d
-
+
n
→
Q
u
d
+
π
0
Q
d
-
+
n
→
Q
d
d
+
π
+
Q
d
-
+
p
→
Q
u
u
+
π
0
Q
d
-
+
p
→
Q
u
d
+
π
+
(c)
Meson-to-Meson2→3 processes
Inelastic scattering
Q
u
-
+
n
→
Q
d
-
+
n
+
π
-
Q
u
-
+
n
→
Q
u
-
+
n
+
π
0
Q
u
-
+
n
→
Q
u
-
+
p
+
π
-
Q
u
-
+
p
→
Q
d
-
+
p
+
π
-
Q
u
-
+
p
→
Q
u
-
+
p
+
π
0
Q
u
-
+
p
→
Q
u
-
+
n
+
π
+
Q
u
-
+
p
→
Q
d
-
+
n
+
π
0
Q
d
-
+
n
→
Q
u
-
+
n
+
π
+
Q
d
-
+
n
→
Q
d
-
+
n
+
π
0
Q
d
-
+
n
→
Q
u
-
+
p
+
π
0
Q
d
-
+
n
→
Q
d
-
+
p
+
π
-
Q
d
-
+
p
→
Q
u
-
+
p
+
π
+
Q
d
-
+
p
→
Q
d
-
+
p
+
π
0
Q
d
-
+
p
→
Q
d
-
+
n
+
π
+
(d)
Meson-to-Baryon2→3 processes
Baryon exchange
Q
u
-
+
n
→
Q
u
d
+
π
-
+
π
0
Q
u
-
+
n
→
Q
d
d
+
π
-
+
π
+
Q
u
-
+
n
→
Q
d
d
+
π
0
+
π
0
Q
u
-
+
p
→
Q
u
u
+
π
-
+
π
0
Q
u
-
+
p
→
Q
u
d
+
π
-
+
π
+
Q
u
-
+
p
→
Q
u
d
+
π
0
+
π
0
Q
u
-
+
p
→
Q
d
d
+
π
+
+
π
0
Q
d
-
+
n
→
Q
u
d
+
π
-
+
π
+
Q
d
-
+
n
→
Q
u
d
+
π
0
+
π
0
Q
d
-
+
n
→
Q
d
d
+
π
0
+
π
+
Q
d
-
+
p
→
Q
u
u
+
π
-
+
π
+
Q
d
-
+
p
→
Q
u
u
+
π
0
+
π
0
Q
d
-
+
p
→
Q
u
d
+
π
0
+
π
+While they interact with the nuclear matter, most of the new states suffer multiple scattering and convert into baryons, allowing forQuu, Qud, and Qdd states, as new heavy mesons are kinematically favoured to increase their baryon number by emitting one or more pions [25]. Due to the lack of these in the material, the reverse reaction is known to be less favoured.Interestingly, new mesons which convert at the beginning of the scattering chain could generate tracks with “dashed-lines,” signalling possible baryonic or electric charge exchange. Indeed,Qq¯ and Q¯q bound states traversing a medium composed of light quarks likely flip their electric charge, frequently interchanging their parton constituents with those of the material nuclei. Eventually, they could leave the detector as heavy stable neutral particles and hence not be observable in muon detectors. The modeling of the nuclear interactions of new heavy hadrons traveling through matter, a survey of which is given in [25, 45], actually favours scenarios with significant charge suppression if heavy open-flavour mesons have a sizeable probability to transform into neutral particles while traversing the detector. This poses a serious challenge for the experimental searches, since the most conservative scenarios assume a complete charge suppression where 100% of the produced hadrons become neutral before reaching the muon detectors. Taking a specific example, we consider pair-produced T2/3 quarks hadronising into the heavy states Tq¯ and T¯q immediately after production. For convenience, we assume that a majority (≳90%) of the new heavy quarks initially hadronise into mesons, while a smaller amount (≲10%) form baryons [6]. As they interact farther in the detector, most of the Tq¯ mesons eventually exit the detector as Tu¯ and Td¯ mesonic states, or as Tud baryons. Baryon-to-meson conversions are allowed as well, as these states can also annihilate back into Tu¯ and Td¯ mesons, yet with a smaller rate. T¯q mesons of vector-like antiquarks, on the other hand, unlikely give rise to antibaryons, but can still flip their charge through the exchange of u and d quarks with the material. Possibly large flavour fractions for T¯u and T¯d can thus be expected in the detector, with a negligible amount of T¯q¯q¯ states. If such antibaryons were to form, they would quickly annihilate by baryon-to-meson interactions, producing pions [46].After travelling through the detectors a few nuclear lengths away from the production vertex, the novel hadrons transform in roughly 100% to the positively charged baryonTud, whereas one half of the T¯-states thus transform into T¯u mesons, and the other half in T¯d. The corresponding hadronic flavour decompositions, evaluated as a function of the penetration depth in the detector, can be read from Figure 5 of [25] for R-hadrons involving stable scalar top and antitop quarks. We notice that similar results are obtained when replacing stop t~ squarks by stable quarks T of equal mass. The Tud, Tu¯ (T¯u), and Td¯ (T¯d), states are then expected to retain the largest flavour composition fractions when considering penetration depths larger than 3 meters and are thus the most likely observable states, which we emphasised in bold font in Table 5.Except forYuu and Y¯u¯u¯, it is interesting to notice that none of the states listed in Table 5 is neutral if involving vector-like quarks with exotic charges. The corresponding bound states might give rise to possibly large fractions of slowly moving charged particles Xud and Yud with Q=+2e and -e, respectively, accompanied by a small amount of Xu¯, Xd¯, Yu¯, and Yd¯ mesons (and charge conjugates). These particles are all electrically charged as they go across the detector, emitting charged and neutral pions in the process together with losing small amounts of energy. Given the small expected interaction cross-section from the Meson-to-Baryon transitions, the processes shown in Table 6 leave small energy deposits in the calorimeters, of the order of O(1) GeV [6]. On the other hand, such stable hadrons would lead to observable tracks due to the ionisation energy losses, allowing for signatures similar to slow-moving muons with high transverse momentum. Searches for stable charged particles, if adapted to such a case, thus provide a promising strategy to rule out the possibility for novel exotic quarks with large lifetimes. Given these very specific signals, the aforementioned signatures can certainly be discriminated at the LHC.


## 3. Experimental Aspects

In this section, we review the past and current experimental searches for new long-lived quarks. Limitations of these searches are then discussed along with possible extensions to enhance sensitivity to heavy quarks. Reinterpretations of direct searches for heavy quarks as well as general searches for long-lived particles are given in the context of long-lived heavy quarks.


### 3.1. Previous Searches at Tevatron

The searches undertaken at the Tevatron resulted in already stringent mass bounds on long-lived heavy quarks, yet not without assumptions. Looking for long-lived parents of theZ boson in displaced vertices from pp¯ collisions at s=1.8 TeV, CDF set limits on the cross-section of a fourth generation charge -1/3e quark as a function of its lifetime with an integrated luminosity of 90 pb-1 [47]. Isolated electron-positron pairs originating from Z decays were searched for in the exclusive channel b′→bZ with mb′>mZ. Finding no evidence for new long-lived particles, CDF excluded mb′<148 GeV/c2 for cτ= 1 cm at 95% CL. This limit drops to mZ+mb≃96 GeV/c2 if cτ> 22 cm or cτ< 0.009 cm. Previously, the D0 collaboration already ruled out the range mZ/2<mb′<mZ+mb from b′→bγ searches in [48] for all proper lifetimes, while a b′ quark with a mass lower than mZ/2 was previously dismissed by the LEP direct searches [49]. Within a 90 pb-1 data sample of pp¯ collisions recorded during 1994-1995, CDF performed a search for low velocity massive charged stable particles leaving large amounts of energy in the calorimeters [50].Assuming a muonlike penetration and searching for an anomalously high ionisation energy loss signature, the data was found to agree with background expectations, and upper limits of the order of 1 pb were derived on the production cross-section. Sensitive to long-lived fourth generation quarks scenarios, the lower boundsmb′>190 GeV/c2 and mt′>220 GeV/c2 have been obtained for q=-1/3e and q=2/3e stable quarks, respectively, with no observed excess over background. More recently, D0 studied a 1.1 fb-1 data sample and looked for Z→e-e+ decays assumed to follow from a long-lived b′ quark parent with BR(b′→bZ)=100% [51]. Assuming that the electromagnetic showers were originating from the same vertex, the analysis found no hint away from the pp¯ interaction point. D0 excluded mb′<190 GeV/c2 at the 95% confidence level for decay lengths between 3.2 mm and 7 m. With the same integrated luminosity, CDF excluded mb′<268 GeV/c2 at 95% CL, considering a long-lived b′ quark decaying exclusively into a Z boson and a b jet [52]. However, it has been emphasised in [53] that the assumption BR(b′→bZ)=100% is inaccurate for mb′>255 GeV/c2, given that the b′ decay mode b′→tW should take over if mb′>mt+mW. Furthermore, the processes b′→(u,c)W were hinted to proceed with nonnegligible rates, leading to an even lower branching ratio. Indeed, even if the b′→tW decay dominates, the aforementioned mass limits depend sensitively on the CKM mixing elements between the fourth and the first three generations. If nonunity couplings between the fourth and the lighter quarks are allowed, the CDF lower bound on mb′ can be significantly affected.Additionally, it is known that the above CDF limits do not apply for long-lived heavy quarks decaying between roughly 1 cm and 3 m within the detector [53]. Should their couplings to SM quarks lie in the corresponding range 4.5×10-9–7.8×10-8, there exist no bounds on the t′ and b′ masses in this uncovered region. For decay lengths larger than the detector dimensions, the lower bounds drop to mt′>220 GeV/c2 and mb′>190 GeV/c2, as obtained from the stable quark searches [50, 51].


### 3.2. LHC

#### 3.2.1. Limitations of Direct Searches

While dedicated results for long-lived bound states of vector-like heavy quarks are still missing at the LHC, various searches for heavy quarks of zero lifetime have already been performed by the ATLAS and CMS experiments. In this section, we discuss how the results of searches for prompt production could be reinterpreted for lifetimes larger than10-10 s.The main aspects of reinterpreting these results are the branching ratios to the investigated final states, which may be different in alternative models, and of course the lifetime of the heavy particle.Considering here a specific example, the current best limit onb′ production, published by the CMS collaboration, excludes production cross-sections of σ>0.1 pb at 95% CL with an integrated luminosity of 4.9 fb-1 at s=7 TeV [54], assuming a branching ratio BR(b′→tW)=100%. As we have seen in Section 2.2, the assumption of exclusive branching ratios is not generally valid for new heavy quarks. Nevertheless, the fraction to which a hypothetical signature occurs can be used to recalculate the cross-section limit in the most naïve approximation
(9)σtrue=σexcluded×BRassumedBRtrue,
where σexcluded is the excluded cross-section under the assumption of BRassumed and BRtrue is the branching ratio in the given model.While this statement is a trivial estimation, the reinterpretation for longer lifetimes requires more thought. In particular, the implicit assumption of heavy quarks decaying promptly at the production vertex is valid only for very short lifetimes. If these particles form bound states and propagate certain distances before their decay, they can potentially escape the direct searches. In such a scenario, the published limits could be considerably weaker, depending on the particle’s lifetimes, simply because they would fly too far to be detected at the primary vertex.For intermediate lifetimes with displaced decay vertices within the detector volume (cf. the region (ii) discussed in Section2.3), a recalculation of the published limits can be attempted. Assuming an exponential decay function, a fraction of the decays will always happen in the vicinity of the beam line, so that the prompt searches will pick them up. Based on this fraction, one can recalculate the limits as a function of the heavy quark lifetime. However, in order to do this correctly, one requires the exact selection efficiency as a function of the displacement from the beam line. Unfortunately, such information has not been made available by the experiments. Still, it can be estimated from an educated guess under specific assumptions. For instance, the analysis in [54] applies b-jet tagging, which requires high-quality tracks to originate within the inner detector volume. The CMS innermost part, the pixel detector consists of three barrel layers with the innermost layer at a distance of 4.4cm from the beam line. This represents the technical limitation to this analysis, so that we can assume that the selection efficiency vanishes at a transverse displacement of around 4cm.Even in analyses without application ofb-jet tagging, stringent quality cuts are usually required to select reconstructed detector objects originating from the primary vertex. One of the motivations for these requirements is the rejection of pileup. Jets from displaced decays which do not point to the production vertex are mostly removed by the pileup cleaning procedures. We therefore make the assumption that events are kept in the direct searches if the decay happens at less than 4cm from the primary interaction vertex, and lost otherwise. The fraction of the lost events can be determined from the distance that these particles travel before their decay. This distance d depends on the mass, the momentum, and the proper lifetime τ. It is defined in the laboratory frame as d=γ·β·c·τ. The fraction is obtained by convoluting d with an exponential decay function and the momentum spectrum which determines γ and β.We take the momentum spectrum of pair-produced heavy quarks as predicted by the MadGraph [55] Monte Carlo generator. The assumed center of mass energy is s=8 TeV. The momentum is needed for a precise calculation of the decay distance on an event by event basis because the velocity may be significantly smaller than c. The pT distribution for various quark masses is shown in Figure 5. The distribution of the velocity in units of c is shown in Figure 6.Figure 5
Transverse momentum distribution of pair-produced heavy quarks as predicted by MadGraph.Figure 6
Velocity distribution of pair-produced heavy quarks in units ofβ=v/c as predicted by MadGraph.Our simulation assumes that the production mechanism and kinematics are independent of the lifetime of the particle. Possible limitations of that assumption are discussed in Section3.2.4. Calculating the distribution of decay vertices in the detector frame, we obtain the fraction of decays outside of the geometrical acceptance of 4cm as a function of the proper lifetime τ. Table 7 summarises our results for four benchmark points in the parameter space. The results are displayed as a function of the proper lifetime in Figure 7.Table 7
Fraction of rejected decays of heavy quarks outside of the geometrical acceptance of4 cm around the primary interaction vertex for three different lifetimes and four different heavy quark masses.
Lifetime
1
0
-
10 s
1
0
-
9 s
1
0
-
8 s
m
=
300 GeV/c2
32%
86.5%
98.5%
m
=
500 GeV/c2
26%
84.1%
98.3%
m
=
700 GeV/c2
21%
82.3%
97.9%
m
=
1000 GeV/c2
16%
79.4%
97.5%Figure 7
Fraction of rejected decays of heavy quarks outside of the geometrical acceptance of4cm around the primary interaction vertex as a function of the proper lifetime τ.We conclude that direct searches are only valid for lifetimes which are considerably shorter than about10-10 s; otherwise the particles would propagate too far and be rejected by the selection criteria. We also stress that this result is relatively independent of the particle masses within our approximations.The published limits on the heavy quark production cross-sections can be simply reinterpreted in the region10-8 s <τ<10-10 s by
(10)σtrue=σexcluded×(1-f),
where f is the fraction of lost events given in Table 7 and Figure 7. The recalculation of the mass limits can then be done by comparing the recalculated cross-sections with predicted cross-sections as a function of the mass.


#### 3.2.2. Displaced Vertices

The immediate question that arises from the previous results is whether lifetimes larger than10-10 s are covered by dedicated searches for long-lived particles, either through displaced vertices or signatures of stable particles propagating through the full detector.Displaced topologies have been searched for at LHC and published, for instance, in [56, 57]. These searches are indirectly applicable to heavy quarks and will be discussed in this section.The CMS analysis [56] is looking for massive long-lived spinless neutral χ bosons, produced in decays of Higgs bosons (H→χχ). The χ bosons are decaying to dileptons χ→l+l-. This analysis is feasible thanks to the CMS track reconstruction algorithm which is able to identify very displaced tracks not originating from the primary interaction vertex. By fitting tracks of two oppositely charged leptons to a common vertex, transverse displacements up to 50cm from the beam line can be reconstructed. This analysis is very sensitive, and it is able to put limits on production cross-sections of the order of σ×BR<10-3 pb at 95% CL, depending on the assumed masses of the χ bosons and their mother particles. The decays of heavy quarks such as Q→tW, for instance, can lead to displaced dilepton vertices as well, because both W bosons can decay leptonically.In this analysis, the momentum of the vertex is required to be parallel to the vector pointing from the primary vertex to theχ boson decay vertex in the transverse plane. This “collinearity” cut is 0.2 (0.8) radians for muon (electron) final states. A direct interpretation of the published limits in context of long-lived heavy quarks is therefore not straightforward. Due to the presence of the neutrinos and of the b-quark jet in Q→tW decays, the momentum of the dilepton vertex will not be parallel to the vector pointing from the primary to the secondary vertex. However, we can attempt to estimate the efficiency of this collinearity cut based on simulated events. We use again our MadGraph simulation from Section 3.2.1 and calculate the angle between the momentum of a heavy quark Q and the dimuon momentum in leptonic Q→tW decays. It is found that the distribution of this collinearity angle is fairly independent of the mass of the heavy quark. The collinearity cut has an efficiency of 68% for electrons and 25% for muons in leptonic Q→tW decays which represents a moderate decrease. A reinterpretation of the published results, taking these efficiencies into account, will be considered in the following.A central component that would be necessary to derive trustworthy limits on heavy quark production from displaced vertex searches is the Monte Carlo simulation of the signal process. The interactions of these heavy particles with the detector material and the efficiency of their reconstruction and selection need to be estimated in order to facilitate the calculation of limits. Difficulties may arise from missing implementations of certain exotic models in the simulation software. Existing simulations of heavy stable particles such as those ofR-hadrons may provide good approximations though [25].We see in [56] that different assumptions for the masses mH and mχ lead to different exclusion limits due to the kinematical properties of the final state. A direct matching of the mass of a heavy quark to mH and mχ is not straightforward. The assumption that the results from [56] are applicable for heavy quarks can therefore only be understood as an approximation. In the following we assume that the limits for higher mH and mχ values are more applicable to heavy quarks, so we only consider those mass points with mH≥400 GeV/c2 and mχ≥50 GeV/c2. Fortunately, the resulting exclusion limits vary only in a small window between 5·10-4 pb and 2·10-3 pb in the dimuon channel for 1cm <cτ<10cm when scanning the various mass points for mH and mχ. Following the strategy of a worst case scenario we use the worst limit of 2·10-3 pb which weakens to 8·10-3 pb taking the effect of the collinearity cut into account. Now the tW final state has only a 1% branching ratio into dimuons which results in an excluded cross-section of σ·BR(Q→tW)>0.8 pb for 1cm <cτ<10cm. In the context of QCD pair production of heavy quarks, the value of 0.8 pb corresponds to a mass of mQ=435 GeV/c2 at s=7 TeV [58]. This can therefore be assumed to be the exclusion limit under the assumption of a branching fraction BR(Q→tW)=100%. Comparing this to the searches at CDF and D0 quoted in Section 3.1, we see a clear improvement of the limits, at least in the range 1cm <cτ<10cm. While the reinterpretation is difficult due to the discussed reasons, much stronger bounds would certainly be possible by including the heavy quark signatures into the displaced vertex analysis.A displaced vertex analysis has also been performed by the ATLAS collaboration [59], following a slightly different strategy. Displaced vertices are reconstructed in an inclusive way, using all displaced tracks in the event as potential seeds. In contrast to CMS, where only opposite-sign dimuon vertices have been used, the ATLAS approach requires at least five tracks at the vertex and an invariant mass of the vertex of at least 10 GeV/c2. To ensure a good fit quality, only vertices within the fiducial barrel pixel detector volume are considered, which means that transverse displacements up to 18cm are considered. In addition, one muon with pT>50 GeV/c is required to be associated with the signal vertex. This analysis has been carried out in the context of an R-parity violating supersymmetric scenario, deriving limits on σ·BR2, where BR2 is the square of the branching ratio for produced squark pairs to decay via long-lived neutralinos to muons and quarks. The requirement that the triggering muon candidate is associated with the displaced vertex ensures that the selection efficiency for each neutralino is independent of the rest of the event. This facilitates a reinterpretation for scenarios with different numbers of long-lived neutralinos in the event.Also in this analysis the limits depend on masses and kinematics of the hypothetical long-lived particles and their mother particles. In the most pessimistic scenario the excluded cross-sections are approximatelyσ·BR2>10-2 pb for lifetimes of 3 mm <cτ<102 mm using an integrated luminosity of 4.4 fb-1 at s=7 TeV. For even larger lifetimes of 102 mm <cτ<103 mm the limits are better than σ·BR2>10-1 pb.The interpretation of these results for decays such as theQ→tW or even Q→bW is possible because these signatures would certainly give rise to decay vertices with very high invariant mass and high track multiplicity. The tW final states have at least one muon in 19% of the cases which means that the excluded values are σ·BR(Q→tW)>0.05 pb for 3 mm <cτ<102 mm. This is about an order of magnitude better than our reinterpretation of the CMS result. The value of 0.05 pb corresponds to mQ=650 GeV/c2 [60].


#### 3.2.3. Heavy Stable Charged Particles

For very long-lived scenarios (cf. the region (iii) discussed in Section2.3), the new heavy states do not decay, but they travel through the full detector. In the following we review recent results by CMS [61] along with an assessment of their relevance for HSCPs.The main identification criteria for HSCPs are high momentum, high ionization, and long time of flight (TOF). The energy lossdE/dx along the track is calculated from the charge deposits in the silicon tracker, while the TOF is obtained from the arrival time in the muon system. These quantities are uncorrelated for SM particles. This noncorrelation is used to estimate the background yields in signal-depleted regions.Models of charge suppression due to interaction of HSCPs with the detector material are considered as well. If charged particles become neutral while propagating through the detector they do not reach the muon chambers. The TOF criterion cannot be used in this case. Therefore, the results have also been estimated without the TOF requirement at the cost of weaker exclusion limits.The obtained limits for a given mass can be very different depending on the assumed particle type. This is due to the predicted kinematics of the particle and its expected detector signal. For instance, the limit on pair production of a scalar top quark with a mass between300 and 800 GeV/c2 is σ< 3 fb for an integrated luminosity of 5 fb-1 at s=7 TeV.Limits for long-lived heavy quarks have not been considered explicitly, unfortunately. However, as an approximation, we make the assumption that a long-lived bound state of a stop would behave similar or equal to a bound state of a heavy up-type quark [25]. This concerns both the interaction with the detector material and the production mechanism which is assumed to proceed via the strong interaction. It is then possible to use the stop exclusion limits directly for stable heavy quarks. Limitations in that assumption are discussed in Section 3.2.4.The predicted production cross-sections of heavy quarks are quite large for low masses (Section2.1). Even with a large systematic safety margin, these cross-sections can be considered to be excluded by the HSCP searches. The heavy quark pair production cross-section for mQ=800 GeV/c2 (the highest mass value considered in [61]), computed at NLO, is σ=9.7 fb at s=7 TeV [58], which can therefore be assumed to be excluded for sufficiently long lifetimes. Heavy quark masses of the order of mQ≈1 TeV/c2 are not yet excluded though, neither in direct searches for prompt decays nor in long-lived searches.To make a quantitative statement about the excluded lifetimes by the HSCP searches, we repeat our simulation from Section3.2.1. To be detected by these analyses, the particles have to be produced at the primary vertex and they have to traverse the full tracking devices, and the muon chambers for the combined tracker+TOF analysis (or the tracker for the tracker-only analysis). The CMS tracker has a length of about 6m and a diameter of 2.8m [62]. We can calculate the fraction of decays inside the tracker volume and assume that these decays will not be reconstructed in the HSCP analyses because of missing hits in the outermost tracker layers.Table8 and Figure 8 summarise our estimations for a few benchmark points and as a function of the proper lifetime. The transition between the two extrema of being fully efficient and losing all decays happens within two orders of magnitudes of the lifetime between 10-9 s and 10-7 s. From these results we conclude that heavy quark masses reaching mQ=800 GeV/c2 can be excluded at 95% CL for lifetimes longer than 10-7 s. As the selection efficiency of the HSCP searches is about 50% for a lifetime of τ=10-8 s, the published limit by the HSCP searches weakens by a factor of two, which is still sufficient to exclude the 800 GeV/c2. The efficiency drops quickly for even shorter lifetimes.Table 8
Fraction of rejected decays of heavy quarks within the geometrical acceptance of the CMS tracking detector. Decays within the tracker are rejected because the particle's trajectory does not reach the outermost tracker layers.
Lifetime
2
·
1
0
-
9 s
1
0
-
8 s
1
0
-
7 s
m
=
300 GeV/c2
90%
47%
7%
m
=
500 GeV/c2
94%
53%
8%
m
=
700 GeV/c2
96%
58%
10%
m
=
1000 GeV/c2
98%
64%
11%Figure 8
Fraction of rejected decays of heavy quarks within the volume of the CMS tracking detector as a function of the proper lifetimeτ.To conclude this section, we would like to mention one additional very interesting experimental strategy. As discussed in detail in [63], heavy metastable particles may be stopped in the detectors and give out-of-time decays. This kind of signature has been investigated by the CMS collaboration [7] using gaps of no collisions between the proton bunches. Particles of extremely long lifetimes (between 10 μs and 1000 s) may produce jets which can be triggered by the calorimeters during periods of no proton collisions. Using a dataset in which CMS recorded an integrated luminosity of 4 fb-1 and a search interval corresponding to 246 hours of trigger live time, 12 events have been observed with a background prediction of 8.6±2.4 events. This result has been interpreted in the context of long-lived gluinos (R-hadrons) and scalar top quarks. The best limit is obtained for lifetimes larger than 10-6 s where the upper limit on the production of stop quarks has been found to be σ(pp→t~t~)·BR(t~→tχ~0)>0.7 pb. When interpreting this for pair production of heavy quarks, the cross-section of 0.7 pb corresponds to mQ=445 GeV/c2 [58]. These results are clearly less stringent than those of the stable HSCP searches discussed above; however, they are able to probe signatures which can potentially escape the HSCP searches in case the particle propagates with an extremely low velocity or does not reach the outermost tracker layers before it stops.


#### 3.2.4. Discussion of the Results

In this section we summarise our findings from Sections3.2.1 to 3.2.3, followed by a discussion of possible shortcomings and alternative search strategies.The following list gives an overview of the different types of searches for heavy or long-lived particles along with the results of our reinterpretations.(i)
Prompt decays of heavy quarks: by far, the largest number of published searches for heavy quarks assumes their prompt production and decay at the proton-proton interaction vertex. We estimated that this assumption is valid for lifetimes up toτ<10-10 s (cτ<3cm). The selection efficiency decreases for longer lifetimes and drops to 10%–20% for τ>10-9 s (cτ>30cm).(ii)
Displaced vertices: our reinterpretation of the ATLAS results [59] for tW final states rules out cross-sections of σ·BR(Q→tW)>0.05 pb for lifetimes of 10-11 s <τ<3·10-10 s (3 mm <cτ<10cm), corresponding to an excluded heavy quark mass of mQ<650 GeV/c2.(iii)
Heavy Stable Charged Particles: these searches require the heavy state to propagate through the full tracking detectors to be detected. Our reinterpretation of the CMS search [61], assuming similar behaviour as for scalar long-lived top quarks, excludes cross-sections of σ>3 fb for very long lifetimes τ>10-7 s (cτ>30m). The mass limit mQ<800 GeV/c2 is still valid for lifetimes τ>10-8 s (cτ>3m).This is obviously still a patchwork of different sensitivities for different lifetimes and different models and methods. We also see that a small geometrical window of inefficiency remains. This is the region beyond displacements of more than about 50 cm but within the volume of the tracking devices. In this region, both the searches for displaced vertices and the searches for HSCP do not reconstruct any events. The exact size of this geometrical region depends on the experiment. We estimate the fraction of decays of heavy quarks within this geometrical region based on our simulation from the previous sections. The result is shown in Figure9. We see that in the worst case scenario only 50% of the decays are lost for a lifetime of τ=3·10-9 s. This represents a moderate decrease of the selection efficiency for a narrow window, so that we can expect a full coverage of the lifetime spectrum in the future.Figure 9
Fraction of decays of heavy quarks between the acceptance of the displaced vertex analysis and HSCP analyses.Yet, some of the assumptions made in this section may have certain shortcomings. First of all, we considered a general framework in which the new heavy quarks are pair-produced, with kinematics independent of the lifetime of the particle. Our motivation was that the current searches for long-livedR-hadrons formed from stable stops generally assume that strong pair production dominates. In supersymmetric models, however, the situation changes drastically when the squark masses become closer to the TeV scale. As pointed out in [64], the squark pair production from light quarks can dominate the gluon fusion mechanism for large masses due to the decreasing probability for the gluons to carry large fractions of the proton momentum. For scenarios where either the squark mass or the mass of the gluino propagator reaches 1 TeV, Higgs and squark weak production starts to dominate. The production channels then become more and more model dependent for increasing masses, arguing for dedicated MC simulations above such mass values.Secondly, we assumed that the heavy quark pair production cross-sections were of the same magnitude as the rate of particles hadronising into bound states, while we stressed in Section2.4 that this is not the case in general. In particular, the ηQ quarkonium production rate is expected to be smaller than the QQ¯ pair production cross-section, due to the possibly large suppression arising from the bound state wave function, which also leads to a different decay phenomenology. Precise calculations would be needed to refine our analysis and the above discussion might need to be adjusted in case of strongly bound states.As a possible suggestion to extend the reach of the above analyses, another search topology may be the case in which a heavy quark decays within the tracker volume. The signature to be searched for may then be one or more displaced jets, which trigger the event, and originate from the heavy quark decay. In addition, a short track with a few hits and largedE/dx pointing to the decay vertex could be required. This signature may have relatively low backgrounds as the jets should have a large transverse momentum. The tracks within these jets are not required to arise from the primary interaction vertex which makes track reconstruction difficult. To circumvent this, one could perform a regional search for short tracklets pointing towards the energy deposit in the calorimeter.


## 3.1. Previous Searches at Tevatron

The searches undertaken at the Tevatron resulted in already stringent mass bounds on long-lived heavy quarks, yet not without assumptions. Looking for long-lived parents of theZ boson in displaced vertices from pp¯ collisions at s=1.8 TeV, CDF set limits on the cross-section of a fourth generation charge -1/3e quark as a function of its lifetime with an integrated luminosity of 90 pb-1 [47]. Isolated electron-positron pairs originating from Z decays were searched for in the exclusive channel b′→bZ with mb′>mZ. Finding no evidence for new long-lived particles, CDF excluded mb′<148 GeV/c2 for cτ= 1 cm at 95% CL. This limit drops to mZ+mb≃96 GeV/c2 if cτ> 22 cm or cτ< 0.009 cm. Previously, the D0 collaboration already ruled out the range mZ/2<mb′<mZ+mb from b′→bγ searches in [48] for all proper lifetimes, while a b′ quark with a mass lower than mZ/2 was previously dismissed by the LEP direct searches [49]. Within a 90 pb-1 data sample of pp¯ collisions recorded during 1994-1995, CDF performed a search for low velocity massive charged stable particles leaving large amounts of energy in the calorimeters [50].Assuming a muonlike penetration and searching for an anomalously high ionisation energy loss signature, the data was found to agree with background expectations, and upper limits of the order of 1 pb were derived on the production cross-section. Sensitive to long-lived fourth generation quarks scenarios, the lower boundsmb′>190 GeV/c2 and mt′>220 GeV/c2 have been obtained for q=-1/3e and q=2/3e stable quarks, respectively, with no observed excess over background. More recently, D0 studied a 1.1 fb-1 data sample and looked for Z→e-e+ decays assumed to follow from a long-lived b′ quark parent with BR(b′→bZ)=100% [51]. Assuming that the electromagnetic showers were originating from the same vertex, the analysis found no hint away from the pp¯ interaction point. D0 excluded mb′<190 GeV/c2 at the 95% confidence level for decay lengths between 3.2 mm and 7 m. With the same integrated luminosity, CDF excluded mb′<268 GeV/c2 at 95% CL, considering a long-lived b′ quark decaying exclusively into a Z boson and a b jet [52]. However, it has been emphasised in [53] that the assumption BR(b′→bZ)=100% is inaccurate for mb′>255 GeV/c2, given that the b′ decay mode b′→tW should take over if mb′>mt+mW. Furthermore, the processes b′→(u,c)W were hinted to proceed with nonnegligible rates, leading to an even lower branching ratio. Indeed, even if the b′→tW decay dominates, the aforementioned mass limits depend sensitively on the CKM mixing elements between the fourth and the first three generations. If nonunity couplings between the fourth and the lighter quarks are allowed, the CDF lower bound on mb′ can be significantly affected.Additionally, it is known that the above CDF limits do not apply for long-lived heavy quarks decaying between roughly 1 cm and 3 m within the detector [53]. Should their couplings to SM quarks lie in the corresponding range 4.5×10-9–7.8×10-8, there exist no bounds on the t′ and b′ masses in this uncovered region. For decay lengths larger than the detector dimensions, the lower bounds drop to mt′>220 GeV/c2 and mb′>190 GeV/c2, as obtained from the stable quark searches [50, 51].


## 3.2. LHC

### 3.2.1. Limitations of Direct Searches

While dedicated results for long-lived bound states of vector-like heavy quarks are still missing at the LHC, various searches for heavy quarks of zero lifetime have already been performed by the ATLAS and CMS experiments. In this section, we discuss how the results of searches for prompt production could be reinterpreted for lifetimes larger than10-10 s.The main aspects of reinterpreting these results are the branching ratios to the investigated final states, which may be different in alternative models, and of course the lifetime of the heavy particle.Considering here a specific example, the current best limit onb′ production, published by the CMS collaboration, excludes production cross-sections of σ>0.1 pb at 95% CL with an integrated luminosity of 4.9 fb-1 at s=7 TeV [54], assuming a branching ratio BR(b′→tW)=100%. As we have seen in Section 2.2, the assumption of exclusive branching ratios is not generally valid for new heavy quarks. Nevertheless, the fraction to which a hypothetical signature occurs can be used to recalculate the cross-section limit in the most naïve approximation
(9)σtrue=σexcluded×BRassumedBRtrue,
where σexcluded is the excluded cross-section under the assumption of BRassumed and BRtrue is the branching ratio in the given model.While this statement is a trivial estimation, the reinterpretation for longer lifetimes requires more thought. In particular, the implicit assumption of heavy quarks decaying promptly at the production vertex is valid only for very short lifetimes. If these particles form bound states and propagate certain distances before their decay, they can potentially escape the direct searches. In such a scenario, the published limits could be considerably weaker, depending on the particle’s lifetimes, simply because they would fly too far to be detected at the primary vertex.For intermediate lifetimes with displaced decay vertices within the detector volume (cf. the region (ii) discussed in Section2.3), a recalculation of the published limits can be attempted. Assuming an exponential decay function, a fraction of the decays will always happen in the vicinity of the beam line, so that the prompt searches will pick them up. Based on this fraction, one can recalculate the limits as a function of the heavy quark lifetime. However, in order to do this correctly, one requires the exact selection efficiency as a function of the displacement from the beam line. Unfortunately, such information has not been made available by the experiments. Still, it can be estimated from an educated guess under specific assumptions. For instance, the analysis in [54] applies b-jet tagging, which requires high-quality tracks to originate within the inner detector volume. The CMS innermost part, the pixel detector consists of three barrel layers with the innermost layer at a distance of 4.4cm from the beam line. This represents the technical limitation to this analysis, so that we can assume that the selection efficiency vanishes at a transverse displacement of around 4cm.Even in analyses without application ofb-jet tagging, stringent quality cuts are usually required to select reconstructed detector objects originating from the primary vertex. One of the motivations for these requirements is the rejection of pileup. Jets from displaced decays which do not point to the production vertex are mostly removed by the pileup cleaning procedures. We therefore make the assumption that events are kept in the direct searches if the decay happens at less than 4cm from the primary interaction vertex, and lost otherwise. The fraction of the lost events can be determined from the distance that these particles travel before their decay. This distance d depends on the mass, the momentum, and the proper lifetime τ. It is defined in the laboratory frame as d=γ·β·c·τ. The fraction is obtained by convoluting d with an exponential decay function and the momentum spectrum which determines γ and β.We take the momentum spectrum of pair-produced heavy quarks as predicted by the MadGraph [55] Monte Carlo generator. The assumed center of mass energy is s=8 TeV. The momentum is needed for a precise calculation of the decay distance on an event by event basis because the velocity may be significantly smaller than c. The pT distribution for various quark masses is shown in Figure 5. The distribution of the velocity in units of c is shown in Figure 6.Figure 5
Transverse momentum distribution of pair-produced heavy quarks as predicted by MadGraph.Figure 6
Velocity distribution of pair-produced heavy quarks in units ofβ=v/c as predicted by MadGraph.Our simulation assumes that the production mechanism and kinematics are independent of the lifetime of the particle. Possible limitations of that assumption are discussed in Section3.2.4. Calculating the distribution of decay vertices in the detector frame, we obtain the fraction of decays outside of the geometrical acceptance of 4cm as a function of the proper lifetime τ. Table 7 summarises our results for four benchmark points in the parameter space. The results are displayed as a function of the proper lifetime in Figure 7.Table 7
Fraction of rejected decays of heavy quarks outside of the geometrical acceptance of4 cm around the primary interaction vertex for three different lifetimes and four different heavy quark masses.
Lifetime
1
0
-
10 s
1
0
-
9 s
1
0
-
8 s
m
=
300 GeV/c2
32%
86.5%
98.5%
m
=
500 GeV/c2
26%
84.1%
98.3%
m
=
700 GeV/c2
21%
82.3%
97.9%
m
=
1000 GeV/c2
16%
79.4%
97.5%Figure 7
Fraction of rejected decays of heavy quarks outside of the geometrical acceptance of4cm around the primary interaction vertex as a function of the proper lifetime τ.We conclude that direct searches are only valid for lifetimes which are considerably shorter than about10-10 s; otherwise the particles would propagate too far and be rejected by the selection criteria. We also stress that this result is relatively independent of the particle masses within our approximations.The published limits on the heavy quark production cross-sections can be simply reinterpreted in the region10-8 s <τ<10-10 s by
(10)σtrue=σexcluded×(1-f),
where f is the fraction of lost events given in Table 7 and Figure 7. The recalculation of the mass limits can then be done by comparing the recalculated cross-sections with predicted cross-sections as a function of the mass.


### 3.2.2. Displaced Vertices

The immediate question that arises from the previous results is whether lifetimes larger than10-10 s are covered by dedicated searches for long-lived particles, either through displaced vertices or signatures of stable particles propagating through the full detector.Displaced topologies have been searched for at LHC and published, for instance, in [56, 57]. These searches are indirectly applicable to heavy quarks and will be discussed in this section.The CMS analysis [56] is looking for massive long-lived spinless neutral χ bosons, produced in decays of Higgs bosons (H→χχ). The χ bosons are decaying to dileptons χ→l+l-. This analysis is feasible thanks to the CMS track reconstruction algorithm which is able to identify very displaced tracks not originating from the primary interaction vertex. By fitting tracks of two oppositely charged leptons to a common vertex, transverse displacements up to 50cm from the beam line can be reconstructed. This analysis is very sensitive, and it is able to put limits on production cross-sections of the order of σ×BR<10-3 pb at 95% CL, depending on the assumed masses of the χ bosons and their mother particles. The decays of heavy quarks such as Q→tW, for instance, can lead to displaced dilepton vertices as well, because both W bosons can decay leptonically.In this analysis, the momentum of the vertex is required to be parallel to the vector pointing from the primary vertex to theχ boson decay vertex in the transverse plane. This “collinearity” cut is 0.2 (0.8) radians for muon (electron) final states. A direct interpretation of the published limits in context of long-lived heavy quarks is therefore not straightforward. Due to the presence of the neutrinos and of the b-quark jet in Q→tW decays, the momentum of the dilepton vertex will not be parallel to the vector pointing from the primary to the secondary vertex. However, we can attempt to estimate the efficiency of this collinearity cut based on simulated events. We use again our MadGraph simulation from Section 3.2.1 and calculate the angle between the momentum of a heavy quark Q and the dimuon momentum in leptonic Q→tW decays. It is found that the distribution of this collinearity angle is fairly independent of the mass of the heavy quark. The collinearity cut has an efficiency of 68% for electrons and 25% for muons in leptonic Q→tW decays which represents a moderate decrease. A reinterpretation of the published results, taking these efficiencies into account, will be considered in the following.A central component that would be necessary to derive trustworthy limits on heavy quark production from displaced vertex searches is the Monte Carlo simulation of the signal process. The interactions of these heavy particles with the detector material and the efficiency of their reconstruction and selection need to be estimated in order to facilitate the calculation of limits. Difficulties may arise from missing implementations of certain exotic models in the simulation software. Existing simulations of heavy stable particles such as those ofR-hadrons may provide good approximations though [25].We see in [56] that different assumptions for the masses mH and mχ lead to different exclusion limits due to the kinematical properties of the final state. A direct matching of the mass of a heavy quark to mH and mχ is not straightforward. The assumption that the results from [56] are applicable for heavy quarks can therefore only be understood as an approximation. In the following we assume that the limits for higher mH and mχ values are more applicable to heavy quarks, so we only consider those mass points with mH≥400 GeV/c2 and mχ≥50 GeV/c2. Fortunately, the resulting exclusion limits vary only in a small window between 5·10-4 pb and 2·10-3 pb in the dimuon channel for 1cm <cτ<10cm when scanning the various mass points for mH and mχ. Following the strategy of a worst case scenario we use the worst limit of 2·10-3 pb which weakens to 8·10-3 pb taking the effect of the collinearity cut into account. Now the tW final state has only a 1% branching ratio into dimuons which results in an excluded cross-section of σ·BR(Q→tW)>0.8 pb for 1cm <cτ<10cm. In the context of QCD pair production of heavy quarks, the value of 0.8 pb corresponds to a mass of mQ=435 GeV/c2 at s=7 TeV [58]. This can therefore be assumed to be the exclusion limit under the assumption of a branching fraction BR(Q→tW)=100%. Comparing this to the searches at CDF and D0 quoted in Section 3.1, we see a clear improvement of the limits, at least in the range 1cm <cτ<10cm. While the reinterpretation is difficult due to the discussed reasons, much stronger bounds would certainly be possible by including the heavy quark signatures into the displaced vertex analysis.A displaced vertex analysis has also been performed by the ATLAS collaboration [59], following a slightly different strategy. Displaced vertices are reconstructed in an inclusive way, using all displaced tracks in the event as potential seeds. In contrast to CMS, where only opposite-sign dimuon vertices have been used, the ATLAS approach requires at least five tracks at the vertex and an invariant mass of the vertex of at least 10 GeV/c2. To ensure a good fit quality, only vertices within the fiducial barrel pixel detector volume are considered, which means that transverse displacements up to 18cm are considered. In addition, one muon with pT>50 GeV/c is required to be associated with the signal vertex. This analysis has been carried out in the context of an R-parity violating supersymmetric scenario, deriving limits on σ·BR2, where BR2 is the square of the branching ratio for produced squark pairs to decay via long-lived neutralinos to muons and quarks. The requirement that the triggering muon candidate is associated with the displaced vertex ensures that the selection efficiency for each neutralino is independent of the rest of the event. This facilitates a reinterpretation for scenarios with different numbers of long-lived neutralinos in the event.Also in this analysis the limits depend on masses and kinematics of the hypothetical long-lived particles and their mother particles. In the most pessimistic scenario the excluded cross-sections are approximatelyσ·BR2>10-2 pb for lifetimes of 3 mm <cτ<102 mm using an integrated luminosity of 4.4 fb-1 at s=7 TeV. For even larger lifetimes of 102 mm <cτ<103 mm the limits are better than σ·BR2>10-1 pb.The interpretation of these results for decays such as theQ→tW or even Q→bW is possible because these signatures would certainly give rise to decay vertices with very high invariant mass and high track multiplicity. The tW final states have at least one muon in 19% of the cases which means that the excluded values are σ·BR(Q→tW)>0.05 pb for 3 mm <cτ<102 mm. This is about an order of magnitude better than our reinterpretation of the CMS result. The value of 0.05 pb corresponds to mQ=650 GeV/c2 [60].


### 3.2.3. Heavy Stable Charged Particles

For very long-lived scenarios (cf. the region (iii) discussed in Section2.3), the new heavy states do not decay, but they travel through the full detector. In the following we review recent results by CMS [61] along with an assessment of their relevance for HSCPs.The main identification criteria for HSCPs are high momentum, high ionization, and long time of flight (TOF). The energy lossdE/dx along the track is calculated from the charge deposits in the silicon tracker, while the TOF is obtained from the arrival time in the muon system. These quantities are uncorrelated for SM particles. This noncorrelation is used to estimate the background yields in signal-depleted regions.Models of charge suppression due to interaction of HSCPs with the detector material are considered as well. If charged particles become neutral while propagating through the detector they do not reach the muon chambers. The TOF criterion cannot be used in this case. Therefore, the results have also been estimated without the TOF requirement at the cost of weaker exclusion limits.The obtained limits for a given mass can be very different depending on the assumed particle type. This is due to the predicted kinematics of the particle and its expected detector signal. For instance, the limit on pair production of a scalar top quark with a mass between300 and 800 GeV/c2 is σ< 3 fb for an integrated luminosity of 5 fb-1 at s=7 TeV.Limits for long-lived heavy quarks have not been considered explicitly, unfortunately. However, as an approximation, we make the assumption that a long-lived bound state of a stop would behave similar or equal to a bound state of a heavy up-type quark [25]. This concerns both the interaction with the detector material and the production mechanism which is assumed to proceed via the strong interaction. It is then possible to use the stop exclusion limits directly for stable heavy quarks. Limitations in that assumption are discussed in Section 3.2.4.The predicted production cross-sections of heavy quarks are quite large for low masses (Section2.1). Even with a large systematic safety margin, these cross-sections can be considered to be excluded by the HSCP searches. The heavy quark pair production cross-section for mQ=800 GeV/c2 (the highest mass value considered in [61]), computed at NLO, is σ=9.7 fb at s=7 TeV [58], which can therefore be assumed to be excluded for sufficiently long lifetimes. Heavy quark masses of the order of mQ≈1 TeV/c2 are not yet excluded though, neither in direct searches for prompt decays nor in long-lived searches.To make a quantitative statement about the excluded lifetimes by the HSCP searches, we repeat our simulation from Section3.2.1. To be detected by these analyses, the particles have to be produced at the primary vertex and they have to traverse the full tracking devices, and the muon chambers for the combined tracker+TOF analysis (or the tracker for the tracker-only analysis). The CMS tracker has a length of about 6m and a diameter of 2.8m [62]. We can calculate the fraction of decays inside the tracker volume and assume that these decays will not be reconstructed in the HSCP analyses because of missing hits in the outermost tracker layers.Table8 and Figure 8 summarise our estimations for a few benchmark points and as a function of the proper lifetime. The transition between the two extrema of being fully efficient and losing all decays happens within two orders of magnitudes of the lifetime between 10-9 s and 10-7 s. From these results we conclude that heavy quark masses reaching mQ=800 GeV/c2 can be excluded at 95% CL for lifetimes longer than 10-7 s. As the selection efficiency of the HSCP searches is about 50% for a lifetime of τ=10-8 s, the published limit by the HSCP searches weakens by a factor of two, which is still sufficient to exclude the 800 GeV/c2. The efficiency drops quickly for even shorter lifetimes.Table 8
Fraction of rejected decays of heavy quarks within the geometrical acceptance of the CMS tracking detector. Decays within the tracker are rejected because the particle's trajectory does not reach the outermost tracker layers.
Lifetime
2
·
1
0
-
9 s
1
0
-
8 s
1
0
-
7 s
m
=
300 GeV/c2
90%
47%
7%
m
=
500 GeV/c2
94%
53%
8%
m
=
700 GeV/c2
96%
58%
10%
m
=
1000 GeV/c2
98%
64%
11%Figure 8
Fraction of rejected decays of heavy quarks within the volume of the CMS tracking detector as a function of the proper lifetimeτ.To conclude this section, we would like to mention one additional very interesting experimental strategy. As discussed in detail in [63], heavy metastable particles may be stopped in the detectors and give out-of-time decays. This kind of signature has been investigated by the CMS collaboration [7] using gaps of no collisions between the proton bunches. Particles of extremely long lifetimes (between 10 μs and 1000 s) may produce jets which can be triggered by the calorimeters during periods of no proton collisions. Using a dataset in which CMS recorded an integrated luminosity of 4 fb-1 and a search interval corresponding to 246 hours of trigger live time, 12 events have been observed with a background prediction of 8.6±2.4 events. This result has been interpreted in the context of long-lived gluinos (R-hadrons) and scalar top quarks. The best limit is obtained for lifetimes larger than 10-6 s where the upper limit on the production of stop quarks has been found to be σ(pp→t~t~)·BR(t~→tχ~0)>0.7 pb. When interpreting this for pair production of heavy quarks, the cross-section of 0.7 pb corresponds to mQ=445 GeV/c2 [58]. These results are clearly less stringent than those of the stable HSCP searches discussed above; however, they are able to probe signatures which can potentially escape the HSCP searches in case the particle propagates with an extremely low velocity or does not reach the outermost tracker layers before it stops.


### 3.2.4. Discussion of the Results

In this section we summarise our findings from Sections3.2.1 to 3.2.3, followed by a discussion of possible shortcomings and alternative search strategies.The following list gives an overview of the different types of searches for heavy or long-lived particles along with the results of our reinterpretations.(i)
Prompt decays of heavy quarks: by far, the largest number of published searches for heavy quarks assumes their prompt production and decay at the proton-proton interaction vertex. We estimated that this assumption is valid for lifetimes up toτ<10-10 s (cτ<3cm). The selection efficiency decreases for longer lifetimes and drops to 10%–20% for τ>10-9 s (cτ>30cm).(ii)
Displaced vertices: our reinterpretation of the ATLAS results [59] for tW final states rules out cross-sections of σ·BR(Q→tW)>0.05 pb for lifetimes of 10-11 s <τ<3·10-10 s (3 mm <cτ<10cm), corresponding to an excluded heavy quark mass of mQ<650 GeV/c2.(iii)
Heavy Stable Charged Particles: these searches require the heavy state to propagate through the full tracking detectors to be detected. Our reinterpretation of the CMS search [61], assuming similar behaviour as for scalar long-lived top quarks, excludes cross-sections of σ>3 fb for very long lifetimes τ>10-7 s (cτ>30m). The mass limit mQ<800 GeV/c2 is still valid for lifetimes τ>10-8 s (cτ>3m).This is obviously still a patchwork of different sensitivities for different lifetimes and different models and methods. We also see that a small geometrical window of inefficiency remains. This is the region beyond displacements of more than about 50 cm but within the volume of the tracking devices. In this region, both the searches for displaced vertices and the searches for HSCP do not reconstruct any events. The exact size of this geometrical region depends on the experiment. We estimate the fraction of decays of heavy quarks within this geometrical region based on our simulation from the previous sections. The result is shown in Figure9. We see that in the worst case scenario only 50% of the decays are lost for a lifetime of τ=3·10-9 s. This represents a moderate decrease of the selection efficiency for a narrow window, so that we can expect a full coverage of the lifetime spectrum in the future.Figure 9
Fraction of decays of heavy quarks between the acceptance of the displaced vertex analysis and HSCP analyses.Yet, some of the assumptions made in this section may have certain shortcomings. First of all, we considered a general framework in which the new heavy quarks are pair-produced, with kinematics independent of the lifetime of the particle. Our motivation was that the current searches for long-livedR-hadrons formed from stable stops generally assume that strong pair production dominates. In supersymmetric models, however, the situation changes drastically when the squark masses become closer to the TeV scale. As pointed out in [64], the squark pair production from light quarks can dominate the gluon fusion mechanism for large masses due to the decreasing probability for the gluons to carry large fractions of the proton momentum. For scenarios where either the squark mass or the mass of the gluino propagator reaches 1 TeV, Higgs and squark weak production starts to dominate. The production channels then become more and more model dependent for increasing masses, arguing for dedicated MC simulations above such mass values.Secondly, we assumed that the heavy quark pair production cross-sections were of the same magnitude as the rate of particles hadronising into bound states, while we stressed in Section2.4 that this is not the case in general. In particular, the ηQ quarkonium production rate is expected to be smaller than the QQ¯ pair production cross-section, due to the possibly large suppression arising from the bound state wave function, which also leads to a different decay phenomenology. Precise calculations would be needed to refine our analysis and the above discussion might need to be adjusted in case of strongly bound states.As a possible suggestion to extend the reach of the above analyses, another search topology may be the case in which a heavy quark decays within the tracker volume. The signature to be searched for may then be one or more displaced jets, which trigger the event, and originate from the heavy quark decay. In addition, a short track with a few hits and largedE/dx pointing to the decay vertex could be required. This signature may have relatively low backgrounds as the jets should have a large transverse momentum. The tracks within these jets are not required to arise from the primary interaction vertex which makes track reconstruction difficult. To circumvent this, one could perform a regional search for short tracklets pointing towards the energy deposit in the calorimeter.


## 3.2.1. Limitations of Direct Searches

While dedicated results for long-lived bound states of vector-like heavy quarks are still missing at the LHC, various searches for heavy quarks of zero lifetime have already been performed by the ATLAS and CMS experiments. In this section, we discuss how the results of searches for prompt production could be reinterpreted for lifetimes larger than10-10 s.The main aspects of reinterpreting these results are the branching ratios to the investigated final states, which may be different in alternative models, and of course the lifetime of the heavy particle.Considering here a specific example, the current best limit onb′ production, published by the CMS collaboration, excludes production cross-sections of σ>0.1 pb at 95% CL with an integrated luminosity of 4.9 fb-1 at s=7 TeV [54], assuming a branching ratio BR(b′→tW)=100%. As we have seen in Section 2.2, the assumption of exclusive branching ratios is not generally valid for new heavy quarks. Nevertheless, the fraction to which a hypothetical signature occurs can be used to recalculate the cross-section limit in the most naïve approximation
(9)σtrue=σexcluded×BRassumedBRtrue,
where σexcluded is the excluded cross-section under the assumption of BRassumed and BRtrue is the branching ratio in the given model.While this statement is a trivial estimation, the reinterpretation for longer lifetimes requires more thought. In particular, the implicit assumption of heavy quarks decaying promptly at the production vertex is valid only for very short lifetimes. If these particles form bound states and propagate certain distances before their decay, they can potentially escape the direct searches. In such a scenario, the published limits could be considerably weaker, depending on the particle’s lifetimes, simply because they would fly too far to be detected at the primary vertex.For intermediate lifetimes with displaced decay vertices within the detector volume (cf. the region (ii) discussed in Section2.3), a recalculation of the published limits can be attempted. Assuming an exponential decay function, a fraction of the decays will always happen in the vicinity of the beam line, so that the prompt searches will pick them up. Based on this fraction, one can recalculate the limits as a function of the heavy quark lifetime. However, in order to do this correctly, one requires the exact selection efficiency as a function of the displacement from the beam line. Unfortunately, such information has not been made available by the experiments. Still, it can be estimated from an educated guess under specific assumptions. For instance, the analysis in [54] applies b-jet tagging, which requires high-quality tracks to originate within the inner detector volume. The CMS innermost part, the pixel detector consists of three barrel layers with the innermost layer at a distance of 4.4cm from the beam line. This represents the technical limitation to this analysis, so that we can assume that the selection efficiency vanishes at a transverse displacement of around 4cm.Even in analyses without application ofb-jet tagging, stringent quality cuts are usually required to select reconstructed detector objects originating from the primary vertex. One of the motivations for these requirements is the rejection of pileup. Jets from displaced decays which do not point to the production vertex are mostly removed by the pileup cleaning procedures. We therefore make the assumption that events are kept in the direct searches if the decay happens at less than 4cm from the primary interaction vertex, and lost otherwise. The fraction of the lost events can be determined from the distance that these particles travel before their decay. This distance d depends on the mass, the momentum, and the proper lifetime τ. It is defined in the laboratory frame as d=γ·β·c·τ. The fraction is obtained by convoluting d with an exponential decay function and the momentum spectrum which determines γ and β.We take the momentum spectrum of pair-produced heavy quarks as predicted by the MadGraph [55] Monte Carlo generator. The assumed center of mass energy is s=8 TeV. The momentum is needed for a precise calculation of the decay distance on an event by event basis because the velocity may be significantly smaller than c. The pT distribution for various quark masses is shown in Figure 5. The distribution of the velocity in units of c is shown in Figure 6.Figure 5
Transverse momentum distribution of pair-produced heavy quarks as predicted by MadGraph.Figure 6
Velocity distribution of pair-produced heavy quarks in units ofβ=v/c as predicted by MadGraph.Our simulation assumes that the production mechanism and kinematics are independent of the lifetime of the particle. Possible limitations of that assumption are discussed in Section3.2.4. Calculating the distribution of decay vertices in the detector frame, we obtain the fraction of decays outside of the geometrical acceptance of 4cm as a function of the proper lifetime τ. Table 7 summarises our results for four benchmark points in the parameter space. The results are displayed as a function of the proper lifetime in Figure 7.Table 7
Fraction of rejected decays of heavy quarks outside of the geometrical acceptance of4 cm around the primary interaction vertex for three different lifetimes and four different heavy quark masses.
Lifetime
1
0
-
10 s
1
0
-
9 s
1
0
-
8 s
m
=
300 GeV/c2
32%
86.5%
98.5%
m
=
500 GeV/c2
26%
84.1%
98.3%
m
=
700 GeV/c2
21%
82.3%
97.9%
m
=
1000 GeV/c2
16%
79.4%
97.5%Figure 7
Fraction of rejected decays of heavy quarks outside of the geometrical acceptance of4cm around the primary interaction vertex as a function of the proper lifetime τ.We conclude that direct searches are only valid for lifetimes which are considerably shorter than about10-10 s; otherwise the particles would propagate too far and be rejected by the selection criteria. We also stress that this result is relatively independent of the particle masses within our approximations.The published limits on the heavy quark production cross-sections can be simply reinterpreted in the region10-8 s <τ<10-10 s by
(10)σtrue=σexcluded×(1-f),
where f is the fraction of lost events given in Table 7 and Figure 7. The recalculation of the mass limits can then be done by comparing the recalculated cross-sections with predicted cross-sections as a function of the mass.


## 3.2.2. Displaced Vertices

The immediate question that arises from the previous results is whether lifetimes larger than10-10 s are covered by dedicated searches for long-lived particles, either through displaced vertices or signatures of stable particles propagating through the full detector.Displaced topologies have been searched for at LHC and published, for instance, in [56, 57]. These searches are indirectly applicable to heavy quarks and will be discussed in this section.The CMS analysis [56] is looking for massive long-lived spinless neutral χ bosons, produced in decays of Higgs bosons (H→χχ). The χ bosons are decaying to dileptons χ→l+l-. This analysis is feasible thanks to the CMS track reconstruction algorithm which is able to identify very displaced tracks not originating from the primary interaction vertex. By fitting tracks of two oppositely charged leptons to a common vertex, transverse displacements up to 50cm from the beam line can be reconstructed. This analysis is very sensitive, and it is able to put limits on production cross-sections of the order of σ×BR<10-3 pb at 95% CL, depending on the assumed masses of the χ bosons and their mother particles. The decays of heavy quarks such as Q→tW, for instance, can lead to displaced dilepton vertices as well, because both W bosons can decay leptonically.In this analysis, the momentum of the vertex is required to be parallel to the vector pointing from the primary vertex to theχ boson decay vertex in the transverse plane. This “collinearity” cut is 0.2 (0.8) radians for muon (electron) final states. A direct interpretation of the published limits in context of long-lived heavy quarks is therefore not straightforward. Due to the presence of the neutrinos and of the b-quark jet in Q→tW decays, the momentum of the dilepton vertex will not be parallel to the vector pointing from the primary to the secondary vertex. However, we can attempt to estimate the efficiency of this collinearity cut based on simulated events. We use again our MadGraph simulation from Section 3.2.1 and calculate the angle between the momentum of a heavy quark Q and the dimuon momentum in leptonic Q→tW decays. It is found that the distribution of this collinearity angle is fairly independent of the mass of the heavy quark. The collinearity cut has an efficiency of 68% for electrons and 25% for muons in leptonic Q→tW decays which represents a moderate decrease. A reinterpretation of the published results, taking these efficiencies into account, will be considered in the following.A central component that would be necessary to derive trustworthy limits on heavy quark production from displaced vertex searches is the Monte Carlo simulation of the signal process. The interactions of these heavy particles with the detector material and the efficiency of their reconstruction and selection need to be estimated in order to facilitate the calculation of limits. Difficulties may arise from missing implementations of certain exotic models in the simulation software. Existing simulations of heavy stable particles such as those ofR-hadrons may provide good approximations though [25].We see in [56] that different assumptions for the masses mH and mχ lead to different exclusion limits due to the kinematical properties of the final state. A direct matching of the mass of a heavy quark to mH and mχ is not straightforward. The assumption that the results from [56] are applicable for heavy quarks can therefore only be understood as an approximation. In the following we assume that the limits for higher mH and mχ values are more applicable to heavy quarks, so we only consider those mass points with mH≥400 GeV/c2 and mχ≥50 GeV/c2. Fortunately, the resulting exclusion limits vary only in a small window between 5·10-4 pb and 2·10-3 pb in the dimuon channel for 1cm <cτ<10cm when scanning the various mass points for mH and mχ. Following the strategy of a worst case scenario we use the worst limit of 2·10-3 pb which weakens to 8·10-3 pb taking the effect of the collinearity cut into account. Now the tW final state has only a 1% branching ratio into dimuons which results in an excluded cross-section of σ·BR(Q→tW)>0.8 pb for 1cm <cτ<10cm. In the context of QCD pair production of heavy quarks, the value of 0.8 pb corresponds to a mass of mQ=435 GeV/c2 at s=7 TeV [58]. This can therefore be assumed to be the exclusion limit under the assumption of a branching fraction BR(Q→tW)=100%. Comparing this to the searches at CDF and D0 quoted in Section 3.1, we see a clear improvement of the limits, at least in the range 1cm <cτ<10cm. While the reinterpretation is difficult due to the discussed reasons, much stronger bounds would certainly be possible by including the heavy quark signatures into the displaced vertex analysis.A displaced vertex analysis has also been performed by the ATLAS collaboration [59], following a slightly different strategy. Displaced vertices are reconstructed in an inclusive way, using all displaced tracks in the event as potential seeds. In contrast to CMS, where only opposite-sign dimuon vertices have been used, the ATLAS approach requires at least five tracks at the vertex and an invariant mass of the vertex of at least 10 GeV/c2. To ensure a good fit quality, only vertices within the fiducial barrel pixel detector volume are considered, which means that transverse displacements up to 18cm are considered. In addition, one muon with pT>50 GeV/c is required to be associated with the signal vertex. This analysis has been carried out in the context of an R-parity violating supersymmetric scenario, deriving limits on σ·BR2, where BR2 is the square of the branching ratio for produced squark pairs to decay via long-lived neutralinos to muons and quarks. The requirement that the triggering muon candidate is associated with the displaced vertex ensures that the selection efficiency for each neutralino is independent of the rest of the event. This facilitates a reinterpretation for scenarios with different numbers of long-lived neutralinos in the event.Also in this analysis the limits depend on masses and kinematics of the hypothetical long-lived particles and their mother particles. In the most pessimistic scenario the excluded cross-sections are approximatelyσ·BR2>10-2 pb for lifetimes of 3 mm <cτ<102 mm using an integrated luminosity of 4.4 fb-1 at s=7 TeV. For even larger lifetimes of 102 mm <cτ<103 mm the limits are better than σ·BR2>10-1 pb.The interpretation of these results for decays such as theQ→tW or even Q→bW is possible because these signatures would certainly give rise to decay vertices with very high invariant mass and high track multiplicity. The tW final states have at least one muon in 19% of the cases which means that the excluded values are σ·BR(Q→tW)>0.05 pb for 3 mm <cτ<102 mm. This is about an order of magnitude better than our reinterpretation of the CMS result. The value of 0.05 pb corresponds to mQ=650 GeV/c2 [60].


## 3.2.3. Heavy Stable Charged Particles

For very long-lived scenarios (cf. the region (iii) discussed in Section2.3), the new heavy states do not decay, but they travel through the full detector. In the following we review recent results by CMS [61] along with an assessment of their relevance for HSCPs.The main identification criteria for HSCPs are high momentum, high ionization, and long time of flight (TOF). The energy lossdE/dx along the track is calculated from the charge deposits in the silicon tracker, while the TOF is obtained from the arrival time in the muon system. These quantities are uncorrelated for SM particles. This noncorrelation is used to estimate the background yields in signal-depleted regions.Models of charge suppression due to interaction of HSCPs with the detector material are considered as well. If charged particles become neutral while propagating through the detector they do not reach the muon chambers. The TOF criterion cannot be used in this case. Therefore, the results have also been estimated without the TOF requirement at the cost of weaker exclusion limits.The obtained limits for a given mass can be very different depending on the assumed particle type. This is due to the predicted kinematics of the particle and its expected detector signal. For instance, the limit on pair production of a scalar top quark with a mass between300 and 800 GeV/c2 is σ< 3 fb for an integrated luminosity of 5 fb-1 at s=7 TeV.Limits for long-lived heavy quarks have not been considered explicitly, unfortunately. However, as an approximation, we make the assumption that a long-lived bound state of a stop would behave similar or equal to a bound state of a heavy up-type quark [25]. This concerns both the interaction with the detector material and the production mechanism which is assumed to proceed via the strong interaction. It is then possible to use the stop exclusion limits directly for stable heavy quarks. Limitations in that assumption are discussed in Section 3.2.4.The predicted production cross-sections of heavy quarks are quite large for low masses (Section2.1). Even with a large systematic safety margin, these cross-sections can be considered to be excluded by the HSCP searches. The heavy quark pair production cross-section for mQ=800 GeV/c2 (the highest mass value considered in [61]), computed at NLO, is σ=9.7 fb at s=7 TeV [58], which can therefore be assumed to be excluded for sufficiently long lifetimes. Heavy quark masses of the order of mQ≈1 TeV/c2 are not yet excluded though, neither in direct searches for prompt decays nor in long-lived searches.To make a quantitative statement about the excluded lifetimes by the HSCP searches, we repeat our simulation from Section3.2.1. To be detected by these analyses, the particles have to be produced at the primary vertex and they have to traverse the full tracking devices, and the muon chambers for the combined tracker+TOF analysis (or the tracker for the tracker-only analysis). The CMS tracker has a length of about 6m and a diameter of 2.8m [62]. We can calculate the fraction of decays inside the tracker volume and assume that these decays will not be reconstructed in the HSCP analyses because of missing hits in the outermost tracker layers.Table8 and Figure 8 summarise our estimations for a few benchmark points and as a function of the proper lifetime. The transition between the two extrema of being fully efficient and losing all decays happens within two orders of magnitudes of the lifetime between 10-9 s and 10-7 s. From these results we conclude that heavy quark masses reaching mQ=800 GeV/c2 can be excluded at 95% CL for lifetimes longer than 10-7 s. As the selection efficiency of the HSCP searches is about 50% for a lifetime of τ=10-8 s, the published limit by the HSCP searches weakens by a factor of two, which is still sufficient to exclude the 800 GeV/c2. The efficiency drops quickly for even shorter lifetimes.Table 8
Fraction of rejected decays of heavy quarks within the geometrical acceptance of the CMS tracking detector. Decays within the tracker are rejected because the particle's trajectory does not reach the outermost tracker layers.
Lifetime
2
·
1
0
-
9 s
1
0
-
8 s
1
0
-
7 s
m
=
300 GeV/c2
90%
47%
7%
m
=
500 GeV/c2
94%
53%
8%
m
=
700 GeV/c2
96%
58%
10%
m
=
1000 GeV/c2
98%
64%
11%Figure 8
Fraction of rejected decays of heavy quarks within the volume of the CMS tracking detector as a function of the proper lifetimeτ.To conclude this section, we would like to mention one additional very interesting experimental strategy. As discussed in detail in [63], heavy metastable particles may be stopped in the detectors and give out-of-time decays. This kind of signature has been investigated by the CMS collaboration [7] using gaps of no collisions between the proton bunches. Particles of extremely long lifetimes (between 10 μs and 1000 s) may produce jets which can be triggered by the calorimeters during periods of no proton collisions. Using a dataset in which CMS recorded an integrated luminosity of 4 fb-1 and a search interval corresponding to 246 hours of trigger live time, 12 events have been observed with a background prediction of 8.6±2.4 events. This result has been interpreted in the context of long-lived gluinos (R-hadrons) and scalar top quarks. The best limit is obtained for lifetimes larger than 10-6 s where the upper limit on the production of stop quarks has been found to be σ(pp→t~t~)·BR(t~→tχ~0)>0.7 pb. When interpreting this for pair production of heavy quarks, the cross-section of 0.7 pb corresponds to mQ=445 GeV/c2 [58]. These results are clearly less stringent than those of the stable HSCP searches discussed above; however, they are able to probe signatures which can potentially escape the HSCP searches in case the particle propagates with an extremely low velocity or does not reach the outermost tracker layers before it stops.


## 3.2.4. Discussion of the Results

In this section we summarise our findings from Sections3.2.1 to 3.2.3, followed by a discussion of possible shortcomings and alternative search strategies.The following list gives an overview of the different types of searches for heavy or long-lived particles along with the results of our reinterpretations.(i)
Prompt decays of heavy quarks: by far, the largest number of published searches for heavy quarks assumes their prompt production and decay at the proton-proton interaction vertex. We estimated that this assumption is valid for lifetimes up toτ<10-10 s (cτ<3cm). The selection efficiency decreases for longer lifetimes and drops to 10%–20% for τ>10-9 s (cτ>30cm).(ii)
Displaced vertices: our reinterpretation of the ATLAS results [59] for tW final states rules out cross-sections of σ·BR(Q→tW)>0.05 pb for lifetimes of 10-11 s <τ<3·10-10 s (3 mm <cτ<10cm), corresponding to an excluded heavy quark mass of mQ<650 GeV/c2.(iii)
Heavy Stable Charged Particles: these searches require the heavy state to propagate through the full tracking detectors to be detected. Our reinterpretation of the CMS search [61], assuming similar behaviour as for scalar long-lived top quarks, excludes cross-sections of σ>3 fb for very long lifetimes τ>10-7 s (cτ>30m). The mass limit mQ<800 GeV/c2 is still valid for lifetimes τ>10-8 s (cτ>3m).This is obviously still a patchwork of different sensitivities for different lifetimes and different models and methods. We also see that a small geometrical window of inefficiency remains. This is the region beyond displacements of more than about 50 cm but within the volume of the tracking devices. In this region, both the searches for displaced vertices and the searches for HSCP do not reconstruct any events. The exact size of this geometrical region depends on the experiment. We estimate the fraction of decays of heavy quarks within this geometrical region based on our simulation from the previous sections. The result is shown in Figure9. We see that in the worst case scenario only 50% of the decays are lost for a lifetime of τ=3·10-9 s. This represents a moderate decrease of the selection efficiency for a narrow window, so that we can expect a full coverage of the lifetime spectrum in the future.Figure 9
Fraction of decays of heavy quarks between the acceptance of the displaced vertex analysis and HSCP analyses.Yet, some of the assumptions made in this section may have certain shortcomings. First of all, we considered a general framework in which the new heavy quarks are pair-produced, with kinematics independent of the lifetime of the particle. Our motivation was that the current searches for long-livedR-hadrons formed from stable stops generally assume that strong pair production dominates. In supersymmetric models, however, the situation changes drastically when the squark masses become closer to the TeV scale. As pointed out in [64], the squark pair production from light quarks can dominate the gluon fusion mechanism for large masses due to the decreasing probability for the gluons to carry large fractions of the proton momentum. For scenarios where either the squark mass or the mass of the gluino propagator reaches 1 TeV, Higgs and squark weak production starts to dominate. The production channels then become more and more model dependent for increasing masses, arguing for dedicated MC simulations above such mass values.Secondly, we assumed that the heavy quark pair production cross-sections were of the same magnitude as the rate of particles hadronising into bound states, while we stressed in Section2.4 that this is not the case in general. In particular, the ηQ quarkonium production rate is expected to be smaller than the QQ¯ pair production cross-section, due to the possibly large suppression arising from the bound state wave function, which also leads to a different decay phenomenology. Precise calculations would be needed to refine our analysis and the above discussion might need to be adjusted in case of strongly bound states.As a possible suggestion to extend the reach of the above analyses, another search topology may be the case in which a heavy quark decays within the tracker volume. The signature to be searched for may then be one or more displaced jets, which trigger the event, and originate from the heavy quark decay. In addition, a short track with a few hits and largedE/dx pointing to the decay vertex could be required. This signature may have relatively low backgrounds as the jets should have a large transverse momentum. The tracks within these jets are not required to arise from the primary interaction vertex which makes track reconstruction difficult. To circumvent this, one could perform a regional search for short tracklets pointing towards the energy deposit in the calorimeter.


## 4. Summary and Conclusions

The current searches for heavy quarks at the LHC mostly ignore the option of long lifetimes. Although a new chiral family of fermions is now strongly disfavoured from the recent Higgs search results, such a scenario remains of prime importance when considering vector-like quarks beyond the Standard Model. As nonchiral fermions decouple in the limit of vanishing mixing with the lighter generations, they can be long-lived and form bound states. Dedicated searches for novel heavyQQ¯ quarkonia, Qqq baryons, and Qq¯ (Q¯q) open-flavour mesons might provide promising strategies for future investigations.As we have detailed, our reinterpretation of the HSCP searches indicates that quark masses of800 GeV/c2 can be excluded at 95% CL for lifetimes longer than 10-8 s. This result, however, assumes that the interactions of the heavy quarks with the detector material are similar to that of stop particles. As discussed in this paper, there are potential restrictions to such an assumption, and the forthcoming searches should be interpreted in a model-specific context. Additionally, it is usually considered that the bound states production cross-sections are comparable to that of pair-produced heavy quarks. This is not necessarily the case for the reasons we presented.Finally, searches for prompt decays are still valid for long lifetimes to a certain extent, but the assumptions about branching ratios and the potential loss of efficiency need to be considered. We showed that these analyses start to become insensitive for lifetimes longer than10-10 s, corresponding to Q-q quark couplings below the 10-9 level. For lifetimes beyond 10-10 s, searches for displaced vertices in the tracker volume have been shown to play an important role.



---

*Source: 690254-2013-04-04.xml*