Unification with Vector-like fermions and signals at LHC
Biplob Bhattacherjee, Pritibhajan Byakti, Ashwani Kushwaha, Sudhir K, Vempati
TL;DR
This paper explores minimal Standard Model extensions with vector-like fermions that achieve gauge coupling unification, include neutrino mass mechanisms or dark matter candidates, and predict new particles detectable at the LHC.
Contribution
It introduces new models with vector-like fermions for gauge unification, some of which incorporate neutrino masses or dark matter, and discusses their LHC signatures.
Findings
Models achieve gauge coupling unification with vector-like fermions.
Some models include Type III seesaw for neutrino masses.
LHC searches already probe these particles around 1 TeV.
Abstract
We look for minimal extensions of Standard Model with vector like fermions leading to precision unification of gauge couplings. Constraints from proton decay, Higgs stability and perturbativity are considered. The simplest models contain several copies of vector fermions in two different (incomplete) representations. Some of these models encompass Type III seesaw mechanism for neutrino masses whereas some others have a dark matter candidate. In all the models, at least one of the candidates has non-trivial representation under . In the limit of vanishing Yukawa couplings, new QCD bound states are formed, which can be probed at LHC. The present limits based on results from 13 TeV already probe these particles for masses around a TeV. Similar models can be constructed with three or four vector representations, examples of which are presented.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
Figure 11
Figure 12
Figure 13
Figure 14
Figure 15
Figure 16
Figure 17
Figure 18
Figure 19
Figure 20
Figure 21
Figure 22
Figure 23
Figure 24
Figure 25
Figure 26| Parameter | Value | Description | |
|---|---|---|---|
| 80.384 0.015 GeV | Pole mass of W boson Agashe:2014kda | ||
| 91.1876 0.0021 GeV | Pole mass of Z boson Agashe:2014kda | ||
| 125.09 0.21 0.11 GeV | Pole mass of Higgs boson Aad:2015zhl | ||
| 173.34 0.76 0.3 GeV | Pole mass of top quark ATLAS:2014wva | ||
| 0.1184 0.0007 | gauge SU(3)c coupling Bethke:2012zza |
| (6,1,1/3) | 1 | ||||
|---|---|---|---|---|---|
| (6,1,2/3) | 1 | ||||
| (8,1,0) | 1 | ||||
| (3,1,1/3) | 1 | ||||
| (3,1,2/3) | 1 |
| Model1 | Rep2 1(6,1,1/3) | ||
|---|---|---|---|
| Model2 | Rep2 2(8,1,0) | ||
| Model3 | Rep2 4(3,1,1/3) | ||
| Model4 | Rep1 2(3,1,2/3) | ||
| Rep2 2(3,2,1/6) | |||
| Model5 | Rep2 1(6,1,2/3) | ||
| Model6 | Rep2 2(6,1,2/3) | ||
| Model7 | Rep1 1(3,1,1/3) | ||
| Rep2 1(3,2,1/6) | |||
| Model8 | Rep2 1(8,1,0) | ||
| Model9 | Rep2 6(3,1,1/3) |
| S.No. | SM Rep | Source | Dynkin Indices | S.No. | SM Rep | Source | Dynkin Indices |
|---|---|---|---|---|---|---|---|
| 1 | 10 | 21 | , | ||||
| 2 | 22 | 45, 70 | |||||
| 3 | 23 | , | |||||
| 4 | 24 | 70 | |||||
| 5 | 25 | ||||||
| 6 | 5, 45, 70 | 26 | |||||
| 7 | 40 | 27 | 15 | ||||
| 8 | 24 | 28 | 50 | ||||
| 9 | 15 | 29 | 35, 40 | ||||
| 10 | 70 | 30 | 75 | ||||
| 11 | 35 | 31 | 70 | ||||
| 12 | 32 | ||||||
| 13 | 33 | 24 | |||||
| 14 | 34 | 40 | |||||
| 15 | 5, 45, 50, 70 | 35 | 45, 50, 70 | ||||
| 16 | 10, 40 | 36 | 75 | ||||
| 17 | 45 | 37 | 35 | ||||
| 18 | 75 | 38 | |||||
| 19 | 10, 15, 40 | 39 | 70 | ||||
| 20 | 24, 75 | 40 |
| Vector-Like Fermion | Couples to |
|---|---|
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
aainstitutetext: Centre for High Energy Physics, Indian Institute of Science, Bangalore 560012bbinstitutetext: Department of Physics, Pandit Deendayal Upadhyaya Adarsha Mahavidyalaya (PDUAM) Eraligool, Karimganj, 788723, Assam, India
Unification with Vector-like fermions and signals at LHC
Biplob Bhattacherjee b
Pritibhajan Byakti a
Ashwani Kushwaha a
Sudhir K Vempati
Abstract
We look for minimal extensions of Standard Model with vector like fermions leading to precision unification of gauge couplings. Constraints from proton decay, Higgs stability and perturbativity are considered. The simplest models contain several copies of vector fermions in two different (incomplete) representations. Some of these models encompass Type III seesaw mechanism for neutrino masses whereas some others have a dark matter candidate. In all the models, at least one of the candidates has non-trivial representation under . In the limit of vanishing Yukawa couplings, new QCD bound states are formed, which can be probed at LHC. The present limits based on results from 13 TeV already probe these particles for masses around a TeV. Similar models can be constructed with three or four vector representations, examples of which are presented.
1 Introduction
For the past few decades, the path to Beyond Standard Model (BSM) physics has been dictated mostly by solutions to the hierarchy problem Csaki:2016kln . However, with no experimental evidence to support this endeavor, from either LEP, Tevatron or the LHC so far, one might wish to explore alternate paths which do not contain a solution to the hierarchy problem. Furthermore, there could solutions to the hierarchy problem which do not introduce any new particles all the way up to GUT scales. The relaxion idea and its variant for example, propose a cosmological solution to the hierarchy problem without introducing any new physics at the weak scale Graham:2015cka ; Espinosa:2015eda ; Choi:2015fiu ; Kaplan:2015fuy .
One of the guiding principles for these alternate paths is the unification of gauge coupling constants. Popular models like split supersymmetry ArkaniHamed:2004yi ; Giudice:2004tc ; ArkaniHamed:2004fb have been proposed which have part of the MSSM particle spectrum at the weak scale and rest (scalar spectrum) at an intermediate scale. The current limits on the stable, long lived R-hadrons which are a prediction of these models are about 1.5- 1.61 TeV Aaboud:2016uth ; Khachatryan:2016sfv . However, this framework depends crucially on the underlying MSSM framework. Generalization without supersymmetry are important to explore.
With this view point, we revisit extensions of the Standard Model with vector-like fermions which lead to precision gauge coupling unification (for earlier works in this direction, please see Liu:2012qua ; Hartmann:2014fya ; Rizzo:1991tc ; Choudhury:2001hs ; Morrissey:2003sc ; Barger:2004sf ; EmmanuelCosta:2005nh ; Barger:2006fm ; Barger:2007qb ; Calibbi:2009cp ; Donkin:2010ta ; Hall:2009nd ; Dermisek:2012ke ; Li:2003zh ; Shrock:2008sb ; Dorsner:2005fq ; Chkareuli:1994ng ; Gogoladze:2010in ; FileviezPerez:2007nh ; Dermisek:2012as ; Dorsner:2014wva ; Xiao:2014kba ). There are several virtues of these models:
(i) They have minimal constraints from electroweak precision parameters, especially from S and T parameters delAguila:2008pw ; Aguilar-Saavedra:2013qpa ; Ishiwata:2015cga ; delAguila:2000rc ; Bizot:2015zaa ; Lavoura:1992np ; Ellis:2014dza , as long as the mixing between vector-like fermions and SM fermions is small.
(ii) They do not lead to any anomalies as they are vector in nature.
(iii) They can be tested directly at the collider experiments like LHC. The kind of signals depend whether on the amount of mixing they have with the Standard Model fields.
(iv) If they have mixing with SM quarks, it is possible that they can be probed indirectly in flavour physics.
To our knowledge, there has not been a recent survey of models containing vector fermions leading to gauge coupling unification. An earlier analysis was done in Ref. Rizzo:1991tc with the available LEP data at that time. We have updated where those models stand in the Appendix, F. In addition to improvements in the gauge coupling measurements and theoretical threshold calculations which are now available at NNLO, an important role is played by the experimental discovery and the (almost) precise determination of the Higgs mass. It has been shown that the Higgs potential becomes unstable from scales close to GeV Degrassi:2012ry , depending on the exact values of the top mass and . Thus a Grand Unified Theory should not only lead to gauge coupling unification but also keep Higgs potential stable all the way up to the GUT scale.
In the models presented here as we will see the Higgs potential naturally remains stable all the way up to the GUT scale. In the view that the primary existence of these vector particles is unification of gauge couplings, we dub them “unificons”. However, as we will see later, these models do not restrict themselves only to unification. In some models, we find solutions with a provision for Type III seesaw mechanism for neutrino masses, and in some others there is a WIMP (Weakly Interacting Massive Particle) dark matter candidate. Thus “unificon” models can indeed have wide phenomenological reach solving other problems in Standard Model like neutrino masses and dark matter.
As a search for all possible models with extra vector-like fermions would be a herculean task, we resort to minimality. We assume unification of gauge couplings á la SU(5). Additional vector-like particles appear as incomplete representations of SU(5). We have looked at all possible incomplete decompositions emanating from SU(5) representations up to dimension 75. The number of copies in each representation is taken to be which is an integer between 1 and 6. The mass range of these additional vector-like fermions is chosen to be TeV, where is a number taken to be approximately between to .
There are no solutions with successful gauge coupling unification as long as the vector-like fermions come in one single representation. This holds true even if increase the number of copies all the way to six, the maximum we have allowed per representation111We have considered the Yukawa couplings of the extra vector-like fermions and the mixing with the SM fermions to be negligible. This can be arranged by imposing discrete symmetry.. The minimal set of successful models with two different representations each with varied number of copies is listed in Table 2. All these models satisfy constraints from proton decay and the stability of the Higgs potential. Both representations come in several copies. Some solution allows for degeneracy between the fermions of the different representations, where as in some cases require non-degeneracy of the fermions in representation 1 and fermions in representation 2.
Interestingly all models have at least one representation with non trivial colour quantum numbers which makes them attractive from LHC point of view. In the limit of negligible Yukawa couplings, these colour states in SU(3) representations of the type 3, 6, 8 form bound states and are produced at LHC. The present limits on these bound states from 13 TeV run of LHC are already touching the 1 TeV mark, depending on the decay mode and the final states. We provide in detail limits on the relevant SU(3) representation bound states.
We also looked for solutions with three and four different representations. Unlike the two representation case, we considered degenerate spectrum for all the vector-like fermions in these two cases. Several solutions are found which are listed in Appendix B and Appendix C. The rest of the paper is organised as follows: In the next Section 2 we recap the essential RG required for gauge coupling unification and stability of the Higgs potential. In Section 3 we present the results for two fermion different representation case. In Section 4 we present the properties of each successful model. In Section 5, we discuss the bound state formalism of the colour vector-like fermions and limit from LHC. We close with a conclusion and outlook. In Appendix D we have tabulated all forty representations of coming from representations upto dimension 75 Slansky:1981yr with their Dynkin index. In Appendix E constraints on mixing between SM fermions with vector-like quark is summarized. Appendix G summaries the two-loop RG equation of Standard Model.
2 Recap of essential RG
2.1 One loop gauge unification
It is well known that gauge couplings do not unify precisely in the Standard Model. If one insists on unification of the guage couplings at the GUT scale, the required is 0.204 (for one loop beta functions) instead of the current experimental value of = 0.23129 0.00005 1674-1137-40-10-100001 . As argued in the introduction, in the present work, we look for additional vector-like matter fermions, close to the weak scale, which can compensate the deviation and lead to successful gauge coupling unification. At the 1-loop level, the beta functions for the three gauge couplings are given as
[TABLE]
where is runs over all the three gauge groups. The functions have the general form:
[TABLE]
Here is quadratic Casimir and is Dynkin index of representation R. and represents vector, Weyl fermion and complex scalar field respectively. For U(1) group and are
[TABLE]
For SU(N) group is defined as follows:
[TABLE]
The following are the list of Dynkin Indices for lower dimensional Representations:
[TABLE]
More complete list on quadratic casimirs can be found in Aranda:2009wh . Within SM, the beta functions take the value
[TABLE]
In the presence of a vector-like fermion at the scale greater than weak scale, given the gauge coupling unification at , the (Eq. (1)) take the form:
[TABLE]
where and capture effect of addition of vector-like fermions at the scale . The parameter is an useful measure of unification of gauge couplings. It is defined as
[TABLE]
Where the second line can be derived from Eq. (6) assuming unification at . The parameters are defined as . In the absence of new vector-like particles, is independent of the running scale . In their presence however, there is a dependence but it is typically mild. For the case where the new particles are close to weak scale TeV, and when , the log factor, is close to one. In this case, the expression for unified theories is given by
[TABLE]
Note that the Eq. (7) can purely be determined from experiments at . It’s value is given by
[TABLE]
In the SM, if we insist on unified gauge couplings at , at the weak scale, takes the value 0.5 clearly in conflict with experiments. In MSSM, turns out to be . Of course, these arguments are valid only at one loop. There is deviation in Eq. (9) when higher loops are considered. In our analysis, most of the successful models have a of to . The above discussion can be easily generalised for more than one Vector field at scales . It has the following general form at the 1-loop level.
[TABLE]
where we assumed the hierarchy of the scales as etc.
2.2 Two loop RG evolution of gauge couplings
To improve the precision in unification of gauge couplings, we consider two loop beta functions. At the two loop level, the beta functions involve Yukawa couplings which makes them model dependent. Vector-like fermions which typically have “bare” mass terms in the Lagrangian, can also mix with the Standard Model fermions through Yukawa interactions if allowed by the gauge symmetry. However, this mixing is subject to strong phenomenological constraints delAguila:2008pw ; Aguilar-Saavedra:2013qpa ; Ishiwata:2015cga ; delAguila:2000rc ; Bizot:2015zaa . A detailed Discussion on the mixing constraints can be found in Appendix E.
In the present analysis, we restrict ourselves to models with minimal or zero vector-like fermion and SM mixing through the Higgs mechanism. With this assumption, we can safely neglect the Yukawa contribution from the new sector to the gauge coupling unification. The RG equations at the two loop level are given by Jones:1981we ; Arason:1991ic ; Luo:2002ey :
[TABLE]
where the first term in the right hand side is due to one-loop which was discussed in the previous subsection. The second term is purely from gauge interactions whereas the third terms involves the Yukawa terms where the suffixes mean the up-type,down-type and lepton-type couplings. The expression for the coefficients appearing in the second term of the above equation are as follows Jones:1981we ; Luo:2002ti :
[TABLE]
where means dimension of the representation and other factors C(R) and T(R) are already defined in Eq. (2).
For the Standard Model, the values of are as follows:
[TABLE]
In the third term of Eq. (12), we have the coefficients and for the standard model particles it has the following form:
[TABLE]
As we are considering the Yukawa couplings between the vector-like fermions with Higgs boson to be negligible222This can be organised by imposing discrete symmetries distinguishing SM partners from vector-like fermions, the contribution of vector-like particles to coefficient can be taken as zero. On the other hand , where is used to indicate contribution from additional vector-like fermions. We’ll give explicit values of for each of the viable models in Section 4.
Two-loop RG running for the Yukawa couplings is given as
[TABLE]
The SM RG for these Yukawa couplings are shown in Appendix G. Here we will address the effect of new fermion fields in RG of Yukawa couplings Luo:2002ti ; Machacek:1984zw ; Machacek:1983fi ; Machacek:1983tz . The one loop beta functions of these couplings are not be affected by new matter(fermion) fields because we considered the Yukawa couplings between the vector-like fermions with Higgs boson to be negligible. Two loop beta functions get contributions from the diagrams shown in Fig. (1), which results in the following terms:
[TABLE]
2.3 Evolution of Higgs Self coupling
The modification of the gauge beta functions in the presence of additional vector-like particles can have implications on the evolution of the Higgs self coupling. At the outset, one might consider that since there are no new large Yukawa couplings333Firstly by assumption, As we will see in the next section, this is automatic in most models as Yukawa couplings with new vector-like fermions are not gauge invariant., the evolution of the Higgs self coupling might be in the safe region. While this is true, the evolution of the SM Yukawa couplings is itself modified in these models as seen in the previous sub-section. It is thus worthwhile to check explicitly the stability of Higgs self coupling along with gauge coupling unification.
To check the Higgs stability we follow EliasMiro:2011aa ; Degrassi:2012ry ; Buttazzo:2013uya who have checked for the stability using three loop beta functions and NNLO matching conditions. We use the beta function of the at the two loop and put a condition that is always positive at all scales of evolution. Two-loop RG running for the Higgs quartic coupling are shown below
[TABLE]
where beta functions for SM Higgs quartic couplings are defined in Appendix G. The effect of new fermion fields in RG of Higgs quartic couplings are:
[TABLE]
To solve the RG equations we need boundary values of the coupling constants and masses at the top mass () scale. The quantities of interest are Higgs quartic coupling (), Yukawa couplings and gauge coupling, which can be calculated in terms of physical observables W-boson mass (), Z-boson mass (), Higgs mass () and at the two loop level. The input parameters are calculated in the -scheme. More detailed can be found in Buttazzo:2013uya . For the RG running we use the central value of Top mass. The input values of SM parameters and couplings are listed in Table 1.
Values of the relevant couplings at scale are as follows:
[TABLE]
where all the parameters with mass dimension has written in GeV. Central values of the above couplings are calculated upto NNLO (Degrassi:2012ry for ) order for all of them except the for which we considered NNNLO Chetyrkin:1999ys ; Chetyrkin:1999qi ; Melnikov:2000qh . The value of , is extracted from the global fit of Ref. Bethke:2012zza in the effective SM with 5 flavours. Including RG running from to at 4 loops in QCD and at 2 loops in the electroweak gauge interactions, and 3 loop QCD matching at to the full SM with 6 flavours, the strong gauge couplig is calculated. The contribution of the bottom and tau Yukawa couplings, are computed from the b-quark mass ( GeV) and Tau mass ( GeV) Xing:2007fb .
Threshold Corrections at GUT Scale
One of the main concerns which remains now is the possible effect of threshold corrections at the GUT scale, which can be quite significant. These corrections are highly model dependent. In some GUT models, with no extra matter at the weak scale (other than the Standard Model particle content), it is possible to achieve gauge coupling unification through large threshold corrections at the GUT scale Lavoura:1993su . While such extreme situations are no longer valid due to the constraint on the stability of the Standard Model Higgs potential, it is still possible that GUT scale threshold corrections could play an important role. To study the impact of threshold corrections on gauge coupling unification, we define the following parameters: and . Note that coincides with when all , at the scale . In the presence of threshold corrections, one could allow for deviations in in terms of at the GUT scale444Another model independent parameterisation for the threshold corrections was presented in Ellis:2015jwa . Defining , we see that is as large as in the Standard Model. In our survey of models below, we have allowed for variations in up to . A more conservative set of models is tabulated in Appendix A which have of 3.
Proton Decay
Models studied in this work can lead to proton decay mediated by the gauge bosons at GUT scale. The lifetime of proton decay is extremely sensitive to the heavy gauge bosons (). For these models, using the simple decay width formulae, we estimate the life time of the proton, where the current experimental value is of order years Takhistov:2016eqm .
3 Gauge coupling unification with vector-like fermions
As mentioned in Introduction, in our search for successful models with gauge coupling unification, we focus on vector-like matter in incomplete representations of SU(5). We have considered (incomplete) representations Slansky:1981yr up-to dimension 75, which contains a 15 of SU(3) of QCD as the largest component. The full list of incomplete representations is presented in Appendix D. As can be seen from the Table 8, there are 40 representations which we have considered. Note that representations 4, 5 in Table 8 do not come as incomplete representations of SU(5) instead they are singlet representations of SU(5). Our search strategy is start with copies of representation , with all the copies degenerate in mass, and look for unification of the gauge couplings. The maximum number of copies is taken to be 10. The number of representation types considered simultaneously is restricted up to four. An important constraint comes from proton decay, which restricts the scale of unification to lie above (at least) GeV. As mentioned above, in addition to unification, we also consider that the Higgs potential should be stable all the way up to the GUT scale. In the computations, we have also varied the input parameters to lie within their two sigma regions. The masses of new vector-like are assumed to lie between 250 GeV - 5 TeV.
For we searched for the mass of the vector-like fermion from 250 GeV - 5TeV, considering number of vector-like fermions . These masses for copies have been considered degenerate for simplicity and no successful model was observed. The simplest solutions we found contain at least two different representation content each with a different number of copies. We call these solutions “minimal unificon models”. These are listed in Table 2. We now explain the notation used in the Table. The two representations considered are called Rep1 and Rep2. The representation is described as , where introduced earlier is the number of copies of the representation, is the representation of the field under the gauge group G of the SM.
Furthermore, in the above, we mentioned only one part of the representation instead of the complete vector multiplet for brevity. For example, actually means . Colored representations like may mean two possibilities: (a) and (b) . On the other hand, the real representations like and are not short-hand notations. In the second last column, the entries are written in units of GeV. Thus except the first model, all the models have unification scale larger than GeV. All models appeared as the solution of one loop RG equation. In the third and fifth columns, we show the mass range of the vector-like fields. One can see that if we increase the mass of one representation, the mass of the other field also increases (as shown in Fig. (2)(b), Fig. (3)(b), Fig. (4)(b), Fig. (5)(b), Fig. (6)(b), Fig. (7)(b), Fig. (8)(b), Fig. (9)(b) and Fig. (10)(b)).
Solutions with three types of representations are also possible. These are listed in Table 6 of Appendix B. Here we made a restricted choice that all the representations and their copies are degenerate in mass of about 1 TeV. As can been seen from the Table, the minimum number of extra vector-like fermions required is seven over the three representations, where as the maximum number is eighteen. All of them have unification scale less than GeV, which puts them at risk with Proton decay. The life time of the proton in these models is of order years which in contrast with the experimental value years Takhistov:2016eqm . The maximum number of representations we have chosen simultaneously is four. Searching for models with different masses for each copy and each representation is computationally very intensive. Thus, we have considered all the four representations and their copies to be degenerate in mass at 1 TeV. The list of successful models is given in Table 7 of Appendix C. The minimum number of vector-like particles required over all representations is five and the maximum is twenty. As with the three representation case, we find that the Unification scale is smaller than GeV with the exception of one model (Model No 17 of Table 7). As before from the arguments of Proton decay, these models can have potentially small proton life times in conflict with experiment. We do not address this issue here.
4 Minimal unificon models
In this section we concentrate on the Minimal vector-like fermion Unification Models. The list of such of models is given in Table 2. Several interesting features are evident from the Table 2.
(a) Except for the first and eigth model, all the models have unification scale above and thus are safe with proton decay. (b) The minimalist model is model 7, with only two vector-like fermions one with a mass range of 0.250-5 TeV and another within a mass range of 250-5000 GeV. This model might have constraints from direct searches of vector-like quarks at LHC and elsewhere if there is significant mixing with SM particles. In its absence, as we assumed here, the bound will be different. We will discuss it in the next section. (c) The maximum number of vector-like fermions needed is nine in Model 9.
We now discuss each of these models in detail.
4.1 Model 1
In this model555We have cross-checked our Two Loop RG equation of this model with the publicly code SARAH Staub:2008uz for consistency., we have six copies of , which we called Rep1, with mass range between 250 GeV to 5000 GeV and one copy of , called Rep2, with mass range from 250 GeV to 5000 GeV. Rep1 field is lepton doublet like field and thus it can interact with right handed electron and the Higgs field through Yukawa interactions. This field mainly decays to gauge bosons like boson and . For the sake of simplicity of the two loop gauge coupling RG running, we impose appropriate symmetries to these fermion doublets. This symmetry cut-down all the Yukawa terms involving these fields at the renormalisable level and only gauge couplings are allowed. Lightest neutral component of these fermions can be a dark matter candidate. This type of dark matter is called inert fermion doublet dark matter Arina:2011cu ; Arina:2012aj . Rep2 is more exotic and at the renormalisation level, it can interact with the gauge bosons only. It cannot decay to any standard model particles. Thus they form bound states. Phenomenology of this is studied in detail in the next Section 5.
For most points in this model vector-like fermions in Rep1 can be degenerate with vector-like fermions in Rep2 () as shown in Fig. (2)(b). However, there could be points in which either of the and are possible. The change in the beta functions in these three possibilities are as follows:
(a)
(I)
[TABLE]
[TABLE]
[TABLE]
(b)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
(c)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
A sample unification point is shown in Fig. (2)(a), six copies of lepton like vector fermions with degenerate mass of 1210 GeV and one copy of Rep2 with a mass of 1260 GeV is considered. The figure shows unification clearly. The running of and are also shown. The panel Fig. (2)(b) has the mass distribution in Rep1-Rep2 mass plane. The model clearly prefers degeneracy of Rep1 and Rep2 for successful unification.
4.2 Model 2
We got six copies of in mass range between 250 GeV to 2000 GeV and two copies of with mass range from 500 GeV to 5 TeV. Similar to the previous model, Rep1 field is lepton like field and thus all the comments are applicable here. Rep2 is gluino like and at the renormalisation level, it can interact with the gluons only and does not have any decay chain. Possibility of any higher dimension decaying operators and its collider phenomenology are studied in the next Section 5.
In the model, is always less than . The change in the beta functions in the two thresholds are as follows:
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
A sample unification point is shown in Fig. (3)(a), six copies of lepton like vector fermions with degenerate mass of 620 GeV and two copy of Rep2 with a mass of 4310 GeV is considered. The figure shows unification of gauge couplings as well as running of and . Mass distribution in Rep1-Rep2 mass plane is shown in Fig. (3)(b).
4.3 Model 3
In this model, we got two copies of and four copies of . The mass ranges of Rep1 and Rep2 are (250 GeV,5 TeV) and (250 GeV,5 TeV) respectively. Rep1 can be a viable candidate of type III Foot:1988aq ; Ma:1998dn seesaw model with fermion mass of M. The neutrino masses are generically given by a factor , where is the vacuum expectation value of the Higgs field. For large M (of the order of GeV), small neutrino masses are generated even for Yukawa couplings of . On the other hand, either smaller Yukawa couplings (which would not effect the RG running) or extended seesaw mechanisms, such as those of the inverse seesaw models delAguila:2008hw , are required to obtain small neutrino masses while keeping close to a few hundreds of GeV. However, we can also impose appropriate symmetries. This symmetry removes all the Yukawa terms involving these fields at the renormalisable level and only gauge couplings are allowed666Seesaw requires Yukawa couplings, our model does not have a seesaw mechanism for neutrino masses.. Neutral component of these fermions is a viable dark matter candidate. This type of dark matter are referred as wino like dark matter and have been discussed in Cirelli:2005uq ; Bhattacherjee:2014dya ; Cirelli:2014dsa ; Moroi:2013sla .
Rep2 has same representation like the down quark. This colour vector-like fermion can form a bound state and annihilate to diphoton, dijet etc. event, which we studied in Section 5.
For most points in this model vector-like fermions in Rep1 can be degenerate with vector-like fermions in Rep2 () as shown in Fig. (2)(b). However, there could be points in which either of the , and are possible. The change in the beta functions in these possibilities are as follows:
(a)
(I)
[TABLE]
[TABLE]
[TABLE]
(b)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
(c)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
In Fig. (4)(a), a sample of gauge couplings, and running is shown with two copies of weak-isospin triplet vector-like fermions with degenerate mass of 800 GeV and four copies of bottom like vector quark with a mass of 3030 GeV. Fig. (4)(b) shows the mass distribution in Rep1-Rep2 mass plane.
4.4 Model 4
This model is interesting as representations of the vector-like matter are like up quarks (Rep1) and left handed quark (Rep2). They appear in two copies for each and their mass ranges are (500 GeV,5 TeV) and (250 GeV,4.5 TeV) respectively. These vector-like quark can be probed at LHC as a bound state, which is studied in Section 5.
In the model, is greater than . However, some parameters of vector like fermion mass have and . The change in the beta functions in these three possibilities are as follows:
(a)
(I)
[TABLE]
[TABLE]
[TABLE]
(b)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
(c)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
The running of gauge couplings, and are shown in Fig. (5)(a), considering two copies of top like vector fermions with degenerate mass of 3175 GeV and two copies of left handed vector-like quark with a mass of 730 GeV. Fig. (5)(b) shows the mass distribution in Rep1-Rep2 mass plane.
4.5 Model 5
This Model consist of 3 copies of vector-like fermion , which is triplet under SU(2) representation (Rep1) and one copy of vector-like fermion which is sextet under SU(3) representation (Rep2). The mass range of Rep1 and Rep2 are (1.8 TeV to 5TeV) and (250 GeV to 950 GeV) respectively. The possible scenarios of Rep1 has been discussed in Model 3 and Rep2 has been mentioned in Model 1 with hypercharge 2/3. In the model, is greater than .
The change in the beta functions in the two thresholds are as follows:
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
A sample unification point is shown in Fig. (6)(a), three copies of weak isospin triplet vector-like fermions with degenerate mass of 4.16 TeV and one copy of color sextet vector-like fermion with a mass of 280 GeV is considered. The Fig. (6)(a) shows unification clearly. Fig. (6)(b) shows the mass distribution in Rep1-Rep2 mass plane.
4.6 Model 6
This model consist of one copy of Rep1= and two copies of Rep2=. The mass range for Rep1 and Rep2 are (250 GeV to 2 TeV) and (1 TeV to 5 TeV) respectively. The Rep1 is fourplet under SU(2) representation and has been studied under minimal dark matter in Ref. Cirelli:2005uq . To our Knowlegde this is a first time it appeared in the unification of gauge coupling. Rep2 is exotic sextet under SU(3), which we discussed in Model 5. In the model, is greater than . The change in the beta functions in the two thresholds are as follows:
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
Gauge coupling unification and running of and are also shown in Fig. (7)(a), with one copy of weak isospin fourplet vector-like fermions with degenerate mass of 1.51 TeV and two copies of color sextet vector-like fermion with a mass of 4.81 TeV. The Fig. (7)(b) has the mass distribution in Rep1-Rep2 mass plane.
4.7 Model 7
This model consist of one copy of Rep1= and two copies of Rep2=. The mass range for Rep1 and Rep2 are (250 GeV to 5 TeV) and (250 GeV to 5 TeV) respectively. Representation one has been discussed in Model 4 with Rep1 having hypercharge 2/3. The difference can been studied with their bound state decay to diphoton channel, as shown in Section 5. In this model, there could be points in which either of the , and are possible. The change in the beta functions in the three conditions are as follows:
(a)
(I)
[TABLE]
[TABLE]
[TABLE]
(b)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
(c)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
A sample unification point is shown in Fig. (8)(a), one copy of bottom like vector fermions with degenerate mass of 4.65 TeV and one copy of left handed quark like vector fermion with a mass of 309 GeV is considered. The figure shows unification clearly. The running of and are also shown. The panel Fig. (8)(b) has the mass distibution in Rep1-Rep2 mass plane.
4.8 Model 8
This model consist of four copies of Rep1= and one copy of Rep2=. The mass range for Rep1 and Rep2 are (300 GeV to 5 TeV) and (300 GeV to 5 TeV) respectively. This representation has been discussed in Model 2 with different number of particles for each represenatation. The difference can been studied with their bound state decay to diphoton channel and dijet, as shown in Section 5.
For most points in this model vector-like fermions in Rep1 can be degenerate with vector-like fermions in Rep2 () as shown in Fig. (9)(b). However, there could be points in which either of the and are possible.
The change in the beta functions in the two thresholds are as follows:
(a)
(I)
[TABLE]
[TABLE]
[TABLE]
(b)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
(c)
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
A sample unification point is shown in Fig. (9)(a), four copies of lepton like vector fermions with degenerate mass of 1.86 TeV and one copy of gluion like vector fermion with a mass of 1.38 TeV is considered. The figure shows unification clearly. The running of and are also shown. The panel Fig. (9)(b) has the mass distibution in Rep1-Rep2 mass plane.
4.9 Model 9
This model consist of three copies of Rep1= and six copies of Rep2=. The mass range for Rep1 and Rep2 are (1.1 TeV to 5 TeV) and (250 GeV to 1.8 TeV) respectively. This representation has been discussed in Model 3 with different number of particle for each represenatation. The difference can been studied with their bound state decay to diphoton channel and dijet, as shown in Section 5. In the model, is greater than . The change in the beta functions in the two thresholds are as follows:
(I)
[TABLE]
[TABLE]
[TABLE]
(II)
[TABLE]
[TABLE]
[TABLE]
A sample unification point is shown in Fig. (10)(a), three copies of weak-isospin triplet vector-like fermions with degenerate mass of 4.6 TeV and six copies of bottom like vector fermion with a mass of 1.6 GeV is considered. The figure shows unification clearly. The running of and are also shown. The panel Fig. (10)(b) has the mass distibution in Rep1-Rep2 mass plane.
5 Collider Signature of Minimal vector-like fermion models
The models listed in Table 2 have several exotic states lying close to electroweak scale, which can be probed at LHC. Models have exotic lepton like states (uncoloured) mostly in doublet, triplet and fourplet representation of SU(2). These states are produced at the LHC through Drell-Yan process and typically have cross-section of the order 10 fb Franceschini:2008pz (roughly slepton production or exotic lepton production). These particles decay through Yukawa interaction to lighter SM leptons. In the limit of vanishing Yukawa couplings, these particles can manifest as missing energy and disappearing charge track at LHC and limits from monojets and disappearing tracks could apply to our model. The LHC at 14 TeV with integrated luminosity 3000 fb*-1* is only sensitive to mass of order 400 GeV Low:2014cba . In the following we will concentrate on the strongly interacting exotic sector; which appears in all the successful models.
5.1 Decay Operators
The models tabulated in the above has exotic fields and some of these fields don’t have renormalisation level decay operators. These fields are (i) , (ii) and (iii) . Now question is whether we can have higher dimensional operators or not. Note that if there exists any higher dimensional operator then there must be some new fields which got integrated out in some higher scales. Now this scale has to be high (close to the GUT scale) as otherwise unification will be disturbed. These higher dimensional operators are suppressed as
[TABLE]
where is the dimension of the operator . Six-dimensional operators are suppressed by square of the GUT scale and thus life-time of the particle is expected to be High ( years). Thus we are focusing only on the five dimensional operators. Any five dimensional operator for decay of such particle must have the forms:
[TABLE]
where in the place of the Higgs and SM fermions fields one can use their conjugate fields. Thus colour charge of the exotic field has to be neutralized by SM fermion to form a five dimensional operator involving the Higgs. In the SM, there is no such field and hence possibility (1) is not possible. For the second case, colour representation of the exotic times that of the SM fermion field must transform as any one of 1, 8, 10 and 27 dimensional representation. However we don’t have SM field with above representation hence, this second possibility is also ruled out. These exotic fields can form a bound state and in the next subsection we’ll discuss this in details.
5.2 Formalism for Bound state
In this section we investigate the possibility of producing bound states of the colour vector-like fermions. The idea of bound state has been studied, in understanding bottom and charm quark through their bound states. For the formation of bound state, we assume the new vector-like fermion () is long lived so that it has time to form a bound state prior to decaying. This condition is easily satisfied in our case, as the Yukawa coupling between the new vector-like fermions and SM particle is assumed to be negligible. The bound state formalism has been studied in Kats:2012ym ; Kats:2009bv , where they focus on pair-produced colour particle Beyond the Standard Model by the observation of diphoton, dijet etc. resonances arising from QCD bound state.
We assume that the only interaction that contribute to the production of bound state is the Standard Model SU(3) colour gauge interaction. We estimate the annihilation rates and parton-level cross-section at leading order, along with NLO MSTW parton distribution functions Martin:2009iq , to compute the LHC signals for = 8 TeV, 13 TeV and 14 TeV evaluated at scale . The production cross section of colour singlet spin zero bound state from constituent vector-like fermion with colour representation 3, 6, 8 are shown in Fig. (11) and Fig. (15). As pointed out in Ref. Younkin:2009zn , NLO corrections to cross-section can increase the diphoton resonance arising from stoponium by 25%. Therefore, large uncertainties are expected in our result of factor of two or so. This still can allow us to constraints minimal vector-like fermion model.
Further uncertainty in our results arises because of limits extracted from ATLAS and CMS result, which is obtained for a fixed spin and production channel. Signal shape have some dependence on the acceptance, intrinsic width and whether a jet is due to parton-level gluon or quark, this adds to some uncertainties.
A pair of near threshold can form a QCD bound state, which we defined as . If the decay width of is smaller than its respective binding energy, it can be observed as a resonance which annihilates to SM particles. For particles () of mass , the Bohr radius of relevant bound state is much smaller than the QCD scale and the velocity of its constituents is non relativistic, we can estimate bound state as modified hydrogenic approximation. For a particle in the colour representation , the potential between and depends on the colour representation of the pair through the casimirs of and as
[TABLE]
where is defined as the running coupling at the scale of the average distance between the two particle in the corresponding hydrogenic state, which is order of the Bohr radius (for which we used Ref. Pineda:1997hz ). The binding energy of the wave functions at the origin for the ground state are given by
[TABLE]
The quantum number of determines the production as well as the decay modes of bound state particle . The cross-section for the bound state to be produced by initial-state partons and is given as
[TABLE]
where is the mass of the bound state, is the velocity of or in CM frame. The production cross-section of any narrow resonance of mass and spin from parton and , and the decay rate of bound state to and , are related by
[TABLE]
where denotes the colour representation of particle .
In the next subsection we will strict ourself to study the colour singlet and spin zero (J=0) bound state system. Assuming the production cross-section of is dominated by gluon fusion. The gluon fusion partonic production cross-section of bound state is given by
[TABLE]
Depending on the quantum number of , bound state can decay to diphoton, dijet, , and channels. The production of preceding pair events produced in proton-proton collisions in LHC can be predicted as .
Here we will identify the channels in which the bound state resonance would be most easily measurable and compute the corresponding cross-section as a function of the mass, colour representation and charge of the constituent particles. The promising final states that we analyzed are diphoton and dijet channels. In the case of SU(2) multiplet the large mass splitting is constrained by Electroweak precision test, which modifies the oblique parameter T and S Lavoura:1992np , hence we have analysed our results in degenerate mass scenario.
5.3 Signals
5.3.1 , , , channel
Any spin half particle can be produced in pairs (in gg collisions ) in an S-wave colour singlet bound state, which can decay as typically narrow , , resonance. The decay width of the , , signal due to spin bound state is given as PhysRevD.35.3366
[TABLE]
[TABLE]
[TABLE]
where , is the third component of the weak isospin for the left and right handed state of the fermion, is the charge of particle, and .
Model No. 4 of minimal fermion model contains constituent of vector-like fermion (3,2,1/6) with SU(2) doublet. This can also decay to channel, which is comparable to channel. The decay width for is given as PhysRevD.35.3366 ,
[TABLE]
where , .
The branching fraction of the isoweak singlet fermions which satisfied the gauge coupling unification and vacuum stability are tabulated in Table 3.
Model No. 4 with vector-like fermion constituent (3,2,1/6), can decay to or , , and channels. With charge Q=-1/3 the branching fraction at mass TeV is , , , and respectively and for Q = 2/3 is , , , respectively. We observed that in a large isoweak SU(2) represenatation the total decay width can be larger than its width into gg.
Both ATLAS and CMS have performed a search of resonant production of photon pairs for scalar particle (J=0). ATLAS ATLAS-CONF-2016-059 analysis is based on data corresponding to an integrated luminosity of 15.4 fb*-1* at =13 recorded in 2015 and 2016. CMS Khachatryan:2016yec data sample correspond to luminosity 12.9 fb*-1* at =13 in 2016, combined statistically with the previous data of 2012 and 2015 at =8 and =13 respectively, with luminosity of 19.7 and 3.3 fb*-1*.
5.3.2 Dijet channel
S-wave bound state with spin can be produced via and annihilating mostly to . For j=1/2 there is also a comparable contribution from S-wave colour octet bound states produced via and annihilating to , which we will not discuss here.
The decay width of signal due to spin colour singlet bound state is,
[TABLE]
Search for narrow resonances decaying to dijet final states in proton-proton collision has been performed by the ATLAS and CMS collaborations using the LHC run data at 8 TeV as well as 13 TeV. CMS Khachatryan:2016ecr study has been performed with integrated luminosity 18.8 fb*-1* at = 8 TeV using a novel technique called data scouting. ATLAS Aad:2014aqa has studied with = 8 TeV using full integrated luminosity of 20.3 fb*-1* masses upto 4.5 TeV.
In run-II, ATLAS Aaboud:2017yvp with centre-of-mass energy = 13 has studied the dijet search using the data collected in 2015 and 2016 with luminosity 3.5 fb*-1* and 33.5 fb*-1* respectively and CMS CMS-PAS-EXO-16-056 has presented a data with luminosity 36 fb*-1* considering masses above 600 GeV.
5.4 Limits on Signals from CMS and ATLAS
In next section we examine the constraints on masses of bound state from dijet and diphoton bounds considering one copy of constituent vector-like fermions. We have used the recent limits of ATLAS and CMS for diphoton resonance at centre of energy =13 TeV from 2015 and as well as 2016 data. Dijet bounds has been considered for centre of energy =8 and 13 TeV from both ATLAS and CMS.
As we have n number of copies of vector-like fermions described in the in Section 4 for two fermions representation, we will give the exclusion limits of vector-like fermion particle occurring in different models with n number of copies in the Table 4.
5.4.1 Dijet Bounds
In Fig. (12)(a)(b) we present the as a function of the mass of the resonance considering one copy of constituent vector-like fermions. The black line is the upper limit on this cross-section from ATLAS Aad:2014aqa 8 TeV and blue line is from CMS Khachatryan:2016ecr 8 TeV data in Fig. (12)(a). Fig. (12)(b) shows the dijet limits from ATLAS(black) ATLAS-CONF-2016-069 13 TeV and CMS(blue) Sirunyan:2016iap 13 TeV data. We can clearly say that the dijet limits are not strong enough to rule any of the models, if they have only one copy of constituent fermions.
In the Fig. (15), we have plotted (black solid line) the projected limit for 14 TeV LHC at 3000 fb*-1* for the dijet cross Section Yu:2013wta . Assuming model, 14 TeV limits on mass of and coupling between gauge field with quark has been calculated in Ref. Yu:2013wta . Using this limit, we have calculated 14 TeV projected limit on dijet cross-section. We have found that mass of vector-like fermion with colour representation six can be excluded up to 800-900 GeV at the HL-LHC.
5.4.2 Diphoton Bounds
The diphoton channel has played a very important role in discovering the Higgs Boson. It can be a very important channel to look at BSM physics. We present the production of diphoton channel as a function of the resonance mass considering one copy of constituent vector-like fermions in Fig. (13). Black line is the upper limit on this cross-section from ATLAS ATLAS-CONF-2016-059 13 TeV and blue line is from CMS Khachatryan:2016yec 13 TeV data. It can be observed that the upper limits on cross-section can give stringent bound on the masses of vector-like fermions ().
There has been searches in , and resonances from these bound states. ATLAS Aad:2015ipg has performed a combination of individual searches in all-leptonic, and all hadronic final states to search for heavy bosons decaying to and with integrated luminosity of 20.3 fb*-1* at 8 TeV. The sensitivity is weaker than channel for ATLAS Aad:2015mna at 8 TeV by around 1000. Both CMS CMS-PAS-EXO-16-025 and ATLAS Aad:2014fha have performed a resonance decaying to at centre-of-mass energy of 8TeV at integrated luminosity 20.3 and 19.7 fb*-1* respectively. Where sensitivity is weaker than diphoton channel is weaker by order 10.
CMS Khachatryan:2016odk has performed a searches in resonance in leptonic channel final decay state at centre of mass energies of 8 and 13 TeV. The bounds are weaker than diphoton bounds by factor of 200. ATLAS ATLAS-CONF-2016-082 has searched for heavy resonance decaying to and pair decaying to leptonic and hadronic channels at a centre of mass energy 13 TeV with total integrated luminosity 13.2 fb*-1*. The sensitivity is still weaker by factor 1000 with respect to diphoton channel.
6 Summary and Outlook
Unification of gauge couplings is one of the most important signatures of a successful Grand Unified Theory beyond the electroweak scale. We look for models with extra vector-like fermions at the weak scale which can lead to successful unification of gauge couplings. With two representation, we find a class of nine models leading to successful unification of gauge couplings. An interesting aspect of these is that all of them contain coloured vector-like fermion in the spectrum. The coloured set of the vector-like fermions can be probed at LHC by looking for bound states formed by them and their probable decays. We have already listed the present bounds from LHC for each successful model. The future runs of LHC are sensitive to further mass ranges of these particles. Finally, it would be interesting to look for complete GUT models with this particle spectrum.
Acknowledgments
The authors thank Tao Han, Rohini Godbole and Xerxes Tata for collaboration in the initial stages of the project. We thank Xerxes Tata for comments on the manuscript and a reference. SKV thanks visits to University of Pittsburgh and University of Hawaii where this project was started. SKV acknowledges support from IUSTFF Grant:JC-Physics Beyond Standard Model/23-2010 during the visit PITT PACC, Department of physics and Astronomy, University of Pittsburgh, USA and Department of Physics and Astronomy, University of Hawaii, Honolulu, USA. The work of B.B. is supported by the Department of Science and Technology, Government of India, under the Grant Agreement number IFA13-PH-75 (INSPIRE Faculty Award). P.B. thanks Palash B. Pal for giving a clue to calculate 6 dimensional representations of generators.
Appendix A Two Representation Case
Here we enlist the models which satisfy gauge coupling unification and positivity of higgs potenatial for Two fermion representation model, with .
Appendix B Three Representation Case
Here we enlist the models which satisfy gauge coupling unification and positivity of higgs potenatial for three fermion representation model. Unlike Two Representation case, we made a restricted choice that all the representations and their copies are degenerate in mass of about 1 TeV, with up to ten copies in each representation. All of the models have unification scale less than GeV, which does not satisfy with Proton decay constraint. The models are listed below in Table 6
Appendix C Four Representation Case
Here we enlist the models which satisfy gauge coupling unification and stable higgs potenatial upto grand unified scale for four fermion representation model. Here also we restricted representations and their copies are degenerate in mass of about 1 TeV. We have allowed for up to ten copies in each model. Except one model, all of the models have unification scale less than GeV, which does not satisfy with Proton decay constraint. The models are listed below in Table 7
Appendix D Representations and Dynkin indices
We considered all the representations coming from representations upto dimension 75. In Table 8, we listed those forty representations Slansky:1981yr with their contribution to beta function (i.e. Dynkin index) considering them as scalar fields. One can straight-forwardly derive corresponding Dynkin indices if the fileds are vector-like fermion just by multiplying the tabulated value with 2 if the representation is real and by multiplying with 4 if the considered representation is complex.
Appendix E Mixing between SM particle with vector-like fermion
In this section we will assume that the new vector-like fermions interact with the SM fermions via Yukawa interactions. SM contains , and , and Higgs doublet, . It can be easily be understood that, among the vector-like fermions considered in this work, new vector-like fermions coupling to the SM ones with renormalisable couplings can only appear in top and bottom partner gauge-covariant multiplets, and in lepton and neutrino partner with definite quantum numbers, which has been studied in Cacciapaglia:2015ixa ; delAguila:2008pw ; Bizot:2015zaa ; delAguila:2000rc ; Cacciapaglia:2010vn ; Gopalakrishna:2011ef ; Aguilar-Saavedra:2013qpa ; Fujikawa:1994we and some of them tabulated in Table 9.
Here we will briefly overview the leading order constraints coming from EW precision tests, direct searches at colliders and Higgs physics. It is reasonable to assume that, only third family of SM fermions have sizable contribution from new vector-like fermions.
E.1 Vector like quarks
Due to mixing of the SM top and bottom quark with vector-like fermions partners, the resulting physical up and down type quark mass eigenstates and may contain non-zero and components, leading to a deviation in their couplings to Z and W bosons. In this case, the relation between weak and mass eigenstates for up quark can be parameterized as two matrices ,
[TABLE]
Similar unitary matrices can be written for down sector. The mixing angles in the left and right sectors are not independent, but have a relation (see also Dawson:2012di ; Fajfer:2013wca ; Atre:2011ae )
[TABLE]
where and are the mass of SM fermion and vector-like fermion respectively.
This mixing gives new contributions to the oblique parameters S and T Peskin:1990zt , which is precisely measured at LEP and SLC. The contributions to S, T in models with T, B singlets and (T B) doublets are studied in Lavoura:1992np ; Aguilar-Saavedra:2013qpa ; Bizot:2015zaa ; Angelescu:2015kga , which would give a constraints in mixing parameters between SM and their vector-like fermions partners. For singlet B quark, the constraints from is strong, which gives upper bound on mixing to be 0.04. For singlet T quark upper bound of is 0.15 to 0.10 for mass range 600 GeV to 2 TeV respectively, from S and T parameter. For (T B) doublet, the constraints from EW precision gives upper bound on to be 0.06 and, between 0.13 to 0.09 for mass range 600 GeV to 2 TeV respectively, considering the splitting between and of 2 GeV.
Direct Searches
A full model of vector-like Quark decaying to SM particles and search strategies to discover at LHC has been studied in Ref. AguilarSaavedra:2009es ; Contino:2008hi ; Aguilar-Saavedra:2013qpa and Ref. within. The singlet Quark decays as,
[TABLE]
The singlet quark decays are
[TABLE]
doublet assuming that they couple to the third generation, are the same as for singlets,
[TABLE]
We would summaries the mass constraints coming from direct searches of VLQ at the LHC.
For Integrated luminosity of 19.5 fb*-1* at TeV CMS Chatrchyan:2013uxa experiment at the Large Hadron Collider searched for the quark decaying into three different final states, bW, tZ, and tH. The search is carried out using events with at least one isolated lepton. The lower limits are set on the T quark mass at 95 confidence level between 687 and 782 GeV for all possible values of the branching fractions into the three different final states assuming strong production.
A search in CMS Khachatryan:2015oba is performed in five exclusive channels: a single- lepton channel, a multilepton channel, two all-hadronic channels optimized either for the bW or the tH decay, and one channel in which the Higgs boson decays into two photons. A statistical combination of these results is performed and lower limits on the T quark mass are set. Depending on the branch- ing fractions, lower mass limits between 720 and 920 GeV at 95 confidence level are found. A search similar to Top like vector quark, heavy B quark vec- tor couplings to W, Z, and H bosons, is carried out by CMS experiment Khachatryan:2015gza . The B quark is assumed to be pair produced and to decay in one of three ways: to tW, bZ, or bH. The search is carried out in final states with one, two, and more than two charged leptons, as well as in fully hadronic final states.Each of the channels in the exclusive final-state topologies is designed to be sensitive to specific combinations of the B quark-antiquark pair decays. A statistical combination of these results gives lower limits on the B quark mass between 740 GeV and 900 GeV with 95 confidence level, depending on the values of the branching fractions of the B quark to tW, bZ, and bH.
ATLAS has also searched for exotic quark, heavy X quark with decaying to tW gives a lower bound of mass 840 GeV Aad:2015mba with 95 C.L. . Quark Y with chagre decaying to Wb gives lower bound of mass 770 GeV Aad:2015kqa with 95 C.L.. The experimental searches assume pair production via strong interactions and prominent decays in the indicated channels.
E.2 Vector like leptons
In this section we discuss new colourless fermions. Weak iso-triplet with zero hyper-charge vector-like fermion can couple to left handed SM fermions and higgs as:
[TABLE]
where the matrix notation of is as follows
[TABLE]
The contribution of to the EW precision parameters is vanishingly small delAguila:2008pw , since the mixing angle are suppressed by and the loop induced mass splitting between the = GeV Ibe:2012sx . In the limit we can realize it as a type III seesaw model Foot:1988aq with neutrino mass .
In the limit , this can be realized as a wino like dark matter Cirelli:2005uq .
SM fermions can also couple to four different possible vector-like leptons, a weak singlet E, a weak doublet L or , a weak triplet . The effect of these vector-like leptons on modification on the Higgs decays, anomalous magnetic moment to the muon and lepton flavour violation decays are studied in Refs. Falkowski:2013jya ; Altmannshofer:2013zba ; Kumar:2015tna ; Chen:2016lsr ; Fujikawa:1994we ; Ishiwata:2015cga .
Direct search
The limits on strongly depend on the SM generation that couples to the heavy leptons. The limits on doublet L, couplings only to the third generation is 270 GeV and coupling with and gives bound of 450 GeV, Ref. Falkowski:2013jya , while the LEP limit remains more constraining in the case of the singlet E, 100 GeV. For the exotic doublet with a doubly-charged component, Ref. Altmannshofer:2013zba reports 320 GeV.
Appendix F Earlier Scan of models By Tom Rizzo
In this section we update the work done in Ref. Rizzo:1991tc . They studied the grand unified theories in context of additional degree of freedom at electroweak scale. S (F) indicates that the quantum numbers following it refer to a complex scalar (vector like fermion) representation. () is the number of fields of type A (B) in the scenario.
[TABLE]
Appendix G Two loop Beta Function
For Standard Model, in Yukawa sector the beta function are Jones:1981we ; Arason:1991ic ; Luo:2002ey
[TABLE]
where one loop contribution are given as
[TABLE]
with
[TABLE]
the two-loop contribution are given as
[TABLE]
[TABLE]
[TABLE]
with
[TABLE]
and
[TABLE]
In Higgs sector we present functions for the quartic coupling:
[TABLE]
where the one loop contribution is given as,
[TABLE]
with
[TABLE]
and the two loop contribution is given as:
[TABLE]
where is the number of generation of fermions in SM.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) C. Csáki and P. Tanedo, Beyond the Standard Model , in Proceedings, 2013 European School of High-Energy Physics (ESHEP 2013): Paradfurdo, Hungary, June 5-18, 2013 , pp. 169–268, 2015. 1602.04228 . DOI . · doi ↗
- 2(2) P. W. Graham, D. E. Kaplan and S. Rajendran, Cosmological Relaxation of the Electroweak Scale , Phys. Rev. Lett. 115 (2015) 221801 , [ 1504.07551 ]. · doi ↗
- 3(3) J. R. Espinosa, C. Grojean, G. Panico, A. Pomarol, O. Pujolàs and G. Servant, Cosmological Higgs-Axion Interplay for a Naturally Small Electroweak Scale , Phys. Rev. Lett. 115 (2015) 251803 , [ 1506.09217 ]. · doi ↗
- 4(4) K. Choi and S. H. Im, Realizing the relaxion from multiple axions and its UV completion with high scale supersymmetry , JHEP 01 (2016) 149 , [ 1511.00132 ]. · doi ↗
- 5(5) D. E. Kaplan and R. Rattazzi, Large field excursions and approximate discrete symmetries from a clockwork axion , Phys. Rev. D 93 (2016) 085007 , [ 1511.01827 ]. · doi ↗
- 6(6) N. Arkani-Hamed, S. Dimopoulos, G. F. Giudice and A. Romanino, Aspects of split supersymmetry , Nucl. Phys. B 709 (2005) 3–46 , [ hep-ph/0409232 ]. · doi ↗
- 7(7) G. F. Giudice and A. Romanino, Split supersymmetry , Nucl. Phys. B 699 (2004) 65–89 , [ hep-ph/0406088 ]. · doi ↗
- 8(8) N. Arkani-Hamed and S. Dimopoulos, Supersymmetric unification without low energy supersymmetry and signatures for fine-tuning at the LHC , JHEP 06 (2005) 073 , [ hep-th/0405159 ]. · doi ↗