Road map through the desert: unification with vector-like fermions
RRoad map through the desert: unificationwith vector-like fermions
Kamila Kowalska and Dinesh Kumar
National Centre for Nuclear ResearchPasteura 7, 02-093 Warsaw, Poland
Abstract
In light of null results from New Physics searches at the LHC, we look at unifica-tion of the gauge couplings as a model-building principle. As a first step, we considerextensions of the Standard Model with vector-like fermions. We present a comprehen-sive list of spectra that feature fermions in two distinct SU (3) C × SU (2) L × U (1) Y representations, in which precise gauge coupling unification is achieved. We derive up-per and lower limits on vector-like masses from proton decay measurements, runningof the strong gauge coupling, heavy stable charged particle searches, and electroweakprecision tests. We demonstrate that due to a particular hierarchy among the massparameters required by the unification condition, complementarity of various experi-mental strategies allows us to probe many of the successful scenarios up to at least10 TeV. ∗ [email protected] † [email protected] a r X i v : . [ h e p - ph ] D ec Introduction
Unification of three fundamental forces of the Standard Model (SM) into a single gaugeinteraction has been an enticing idea since the mid 1970s [1, 2, 3, 4, 5]. It emerged as a naturalcontinuation of intellectual efforts that, in merging apparently unrelated phenomena, soughtthe key to a deeper understanding of nature, first by combining electricity and magnetism intoa unified description, later leading to the establishment of the electroweak theory. Althoughthe concept of unification as an underlying organizing principle stems to some extent froma sense of aesthetics, it finds a more robust justification in the fact that the renormalizedgauge couplings of the SM, while evaluated at higher and higher energies, seem to convergetowards a common value. This behavior might be understood as a manifestation of a new,unified description of fundamental interactions known as Grand Unified Theory (GUT).Precise gauge coupling unification, however, is not really achieved in the SM as discrep-ancies among the GUT-scale values of the SM couplings reach several percent. To make itwork, the particle spectrum needs to be extended in order to modify the renormalizationgroup (RG) running of the couplings below the GUT scale. Supersymmetry (SUSY) hasthe advantage of leading to gauge unification in a quite natural way, yet no experimentalevidence of the low-scale SUSY has been found so far. While this fact does not undercut itcompletely as a theoretical framework, it is timely to ask to what extent unification of thegauge couplings is a unique property among various extensions of the SM. In other words,how many different beyond-the-SM (BSM) scenarios can be found whose particle spectrumdiffer quantitatively from the one of SUSY, and still allow for precise unification.In addressing this question we would like to remain as generic as possible. On theother hand, a truly comprehensive study of all imaginable SM extensions would be a highlychallenging (if not impossible) task. For that reason our approach will be incremental: we aregoing to begin with a relatively simple BSM setup, which will then be gradually extended toencompass more complex structures. It is in this spirit that we regard the issue of unificationas a long term research project, a road map that would guide the model building throughthe desert between the electroweak (EW) and GUT scales.We begin with defining the common framework for any unification analysis that we aregoing to undertake. The most important requirement is that the SM gauge symmetry persistsup to the unification scale. It means, we will not consider Pati-Salam [2] or trinification [11,12] type of GUTs as they do not require simultaneous unification of all three SM gaugecouplings. For theoretical consistency, we also demand perturbativity of the renormalizedmodel parameters up to the GUT energies.Now we are going to make several additional assumptions, which on the one hand willsubstantially simplify the analysis, on the other will restrict the types of BSM scenariosthat will be considered. Therefore, such assumptions may be dropped in the future stud-ies. The additional requirements we impose are the following: (i) any extension of the SMmust be anomaly free; (ii) scenarios with low scale unification, M GUT ∼ < GeV, are notallowed (i.e. general dimension-6 operators leading to proton decay are not forbidden by In principle, the GUT-scale values of the couplings can also be modified by high-scale threshold correc-tions [6, 7, 8, 9]. These corrections, however, are strongly model dependent and for a certain range of theGUT-particle masses they become negligibly small [10]. Therefore, we neglect the effects of GUT thresholdcorrections throughout this study. SU (3) C × SU (2) L × U (1) Y charge assignements for BSM fermions and scalars. Relativelyrecently the first attempt has been made in Ref. [28] to systematically study all possible VLextensions of the SM, in which BSM matter multiplets form incomplete representations of SU (5). Scenarios with two distinct representations (and no more than six VL pairs in eachof them) were considered, while independent VL masses were limited to 5 TeV. In the present work, we build on the findings of Ref. [28] and extend their analysis inseveral different directions. First of all, we boost the allowed mass range of VL fermionsup to 10 TeV. While it may seem far beyond the reach of modern colliders, we will showthat due to a particular hierarchy among the VL spectra allowing for unification, as wellas to complementarity of various experimental strategies, one is actually able to deriveexclusion lower bounds even on multi-TeV masses. Secondly, we do not a priori limit themaximum number of VL pairs in each representation. It turns out that, when this extracondition is discarded, novel solutions with respect to Ref. [28] can be found. Finally, wethoroughly discuss a variety of experimental methods that allow one to test the successfulunification scenarios. We derive upper and lower limits on VL masses from proton decaymeasurements, running of the strong gauge coupling, heavy stable charged particle searches,and EW precision tests. We demonstrate that by combining independent experimentalresults we manage in many cases to probe (and to exclude) essentially the whole parameterspace of a given model.The paper is organized as follows. In Sec. 2 we define the fundamental building blocks ofour BSM scenarios in terms of the transformation properties of VL fermions under the SMgauge symmetry group. Sec. 3 presents the main results of the study: a comprehensive listof representations that allow for unification of the SM gauge couplings. In Sec. 4 we discussin detail experimental bounds that constrain the parameter space of the successful scenarios.We present our conclusion in Sec. 5. Technical details of the analysis are collected in twoappendices.
We begin our discussion with constructing a set of generic extensions of the SM that satisfythe requirements defined in the introduction, i.e. no extra gauge symmetry is imposed, Results for 3- and 4-representation scenarios and the VL mass fixed at 1 TeV were shown in Ref. [28] aswell. However, as we will demonstrate in the present study, mass hierarchy among various VL representationsis one of the main factors of the successful gauge coupling unification. For this reason, the fixed-mass analysescan not be considered comprehensive. R ∞
24 4 4
28 8 2 2 Table 1: Maximal number of VL fremions withthe mass of 10 TeV, which allows for perturba-tive gauge couplings below 10 GeV. Y F = 0 isassumed. R R
13 56 23 Table 2: Maximal value of the hypercharge,which allows for perturbative gauge couplings be-low 10 GeV ( N F = 2 and the VL mass is set at10 TeV). the only new particles in the spectrum are VL fermions, and the renormalized couplingsremain perturbative at the unification scale. We additionally assume for the purpose of thisstudy that any Yukawa interactions generated by the BSM sector and allowed by the gaugesymmetry are negligible. We thus introduce N F i copies of new fermionic fields, which transform under the SU (3) C × SU (2) L × U (1) Y gauge group as VL multiplets( R F i , R F i , Y F i ) ⊕ ( ¯ R F i , ¯ R F i , − Y F i ) . (1)Note that we count separately over both components of the pair, so N F i can only assumeeven values (with an exception of fermions that transform in an adjoint representation of anon-abelian gauge symmetry group). The index i runs over the number of distinct represen-tations. At this stage both i and N F i are unconstrained.Upper bounds on the dimension of possible VL representations and on the number offermions that transform accordingly are provided by perturbativity condition. Let us firstconsider the SM extended by one representation of VL fermions and assume Y F = 0. Forvarious combinations of R and R and increasing number of VL copies, we run the SMgauge couplings from the low-energy scale, which we identify with the top quark mass M t and at which the couplings assume the following values [37] g ( M t ) = 1 . , g ( M t ) = 0 . , g Y ( M t ) = 0 . , (2)up to 10 GeV. The VL mass is fixed at 10 TeV as this is the largest allowed value of thisparameter considered in the present study. RG equations (RGEs) for the gauge couplings ina general quantum field theory are well known [38] and we summarize their explicit two-loopform in Appendix A.In Table 1 we show the maximal allowed number of VL fermions, N F max , for which thegauge couplings remain perturbative ( g i (cid:46) π ) up to 10 GeV. One can see that color octetsand electroweak quadruplets are the highest representations possible, and that a total numberof 11 different combinations of SU (3) C and SU (2) L charges is allowed. Non-zero hypercharge The impact of non-gauge interactions on unification of the gauge couplings will be discussed elsewhere. Note that we consider this particular value as a rough estimate of the limit imposed on the value of theGUT scale by proton decay measurement. The actual experimental bounds will be discussed in Sec. 4 N F max . Thus, the requirement of perturbativity up to around theunification scale reduces the possible number of VL fermions that transform under SU (3) C and SU (2) L . In some cases, however, N F max exceeds the maximum number of 12 VL copiesadpoted in Ref. [28]. We will demonstrate in Sec. 3 that several novel solutions with respectto those presented in Ref. [28] can be found if that somewhat arbitrary assumption is relaxed.So far our discussion was quite generic as the conclusions regarding the properties of theallowed VL representations resulted merely from the requirement of perturbativity, indepen-dently on what happens at the unification scale and how the expected GUT symmetry isrealized. It is, however, not the case for the hypercharge. This particular quantum num-ber is much more difficult to deal with in a general manner, as in principle it can assumecontinuous values. Additionally, hypercharge normalization is not unique as it depends on aparticular embedding of the SM into a GUT gauge group [39], and different normalizationsmay lead to different predictions regarding the gauge coupling unification. For these reasonswe have to depart at this point from an entirely model-independent approach.We assume from now on that at the unification scale the SU (5) symmetry is restoredand VL fermions are embedded into multiplets of SU (5) just like it is the case for the SMfields. This seems to be the most natural choice since SU (5) not only can play the role ofa self-contained unified gauge symmetry [1], but also shows up in breaking chains of largerGUT groups. In Appendix B the decomposition of the irreducible SU (5) representations ofincreasing dimensions into irreducible representations of the SM gauge group is summarized.It is enough to consider representations up to dimension 75, since the larger ones decomposeeither to representations that have already appeared in Eq. 30, or to representations whosedimensions exceed the limits presented in Table 1. Additionally, in Table 2 we provideinformation about the maximal value of hypercharge for a single pair of VL fermions ( N F =2), for which the gauge couplings remain perturbative up to 10 GeV. When combined,perturbativity bounds in Tables 1 and 2 eliminate some of representations listed in Eq. 30.Eventually, we are left with a set of 24 distinct non-singlet SU (3) C × SU (2) L × U (1) Y representations:color singlets : ( , , , ( , , − , (cid:0) , , (cid:1) , (cid:0) , , − (cid:1) , ( , , , ( , , , (3) (cid:0) , , (cid:1) , (cid:0) , , − (cid:1) , color triplets : (cid:0) , , − (cid:1) , (cid:0) ¯3 , , − (cid:1) , (cid:0) ¯3 , , (cid:1) , (cid:0) ¯3 , , − (cid:1) , (cid:0) , , (cid:1) , (cid:0) ¯3 , , (cid:1) , (cid:0) ¯3 , , − (cid:1) , (cid:0) , , − (cid:1) , (cid:0) ¯3 , , − (cid:1) , color sextets : (cid:0) ¯6 , , − (cid:1) , (cid:0) , , − (cid:1) , (cid:0) ¯6 , , (cid:1) , (cid:0) , , (cid:1) , color octets : ( , , , ( , , , (cid:0) , , (cid:1) . These are the fundamental building blocks of the VL unification scenarios we are going toanalyze in the next section.
We are now in the position to perform a comprehensive analysis of the SM extensions withVL fermions that could potentially lead to precise unification of three SM gauge couplingsat the energies in the range [10 − ] GeV. With this goal in mind, we would like5o proceed in a systematic way, gradually increasing the complexitity of the constructedmodels. The simplest scenarios could be engineered by adding to the SM one of the VLrepresentations listed in Eq. 3. It is, however, a known fact [28] that precise unification is notpossible within such a framework. The next possibility is then to consider two different VLrepresentations with an arbitrary number of fermion copies within each of them. Theoreticaland phenomenological properties of such models are the main objective of the present study.Adding three (or more) independent VL representations makes our task more and morechallenging. Note that 276 distinctive combinations of the VL representations listed in Eq. 3need to be considered in the two-representation case. This figure increases to 2024 whenthree, and to 10626 when four different representations are considered. Each combinationrequires to scan over the numbers of VL fermions in each representation (see Table 1), as wellas on their masses, which we always assume to be uncorrelated. This means that numericalcomplexity of the problem grows exponentially with every independent representation added.For this reason we focus in this study on the simplest case, leaving more complicated SMextensions for future work.The numerical procedure employed in our analysis is the following. We use the 2-loopSM RGEs from M t up to the scale M , at which the lightest of VL fermions show up inthe spectrum. We assume for simplicity that all N F copies of the same representation havea common mass. Above M we switch to the 2-loop RGEs for a generic BSM scenario,Eq. 15-17. At the scale M the effects of heavier VL fermions need to be taken into account.Finally, we define a unification scale, M GUT , as the scale at which all three gauge couplingsacquire a common value, g GUT , g GUT = g ( M GUT ) = g ( M GUT ) = g ( M GUT ) . (4)We require that the unified coupling is perturbative, i.e. g GUT ≤ π .Precision of the gauge coupling unification can be quantified by a set of three mismatchparameters, (cid:15) , , . For each combinations ( i, j ), where i, j = 1 , ,
3, we determine a twocoupling unification scale from the condition g i ( M ij GUT ) = g j ( M ij GUT ) = g ij . We then define adeviation of the third coupling from g ij as (cid:15) k = g k ( M ij GUT ) − g ij g ij (5)and determine the true unification scale by requiring (cid:15) GUT = min( (cid:15) , (cid:15) , (cid:15) ). In the SM (cid:15) SMGUT = 7 . (cid:15) MSSMGUT = 1 .
1% whenall sparticle masses set at 1 TeV. Therefore, we define the precise gauge unification (PGU)by a condition (cid:15) < N , N , and their masses, M , M , and for each point in the 4-dimensionalparameter space we determine (cid:15) GUT and M GUT . The parameters M and M are variedbetween 0.25 TeV and 10 TeV. The main results of the analyses are summarized in Table 3.We found 13 different scenarios that allow for the PGU at the scale 10 − GeV. Nine ofthem have been previously identified in [28], while the scenarios F2, F3, F6 and F13 present6 cenario R F R F N N VL mass/GUT scaleF1 (cid:0) , , (cid:1) (cid:0) , , (cid:1)
12 2 Fig. 1(a)F2 (cid:0) , , (cid:1) (cid:0) , , (cid:1)
20 4 Fig. 1(b)F3 (cid:0) , , (cid:1) (cid:0) , , (cid:1)
22 4 Fig. 1(c)F4 (cid:0) , , (cid:1) ( , ,
0) 8 1 Fig. 1(d)F5 (cid:0) , , (cid:1) ( , ,
0) 12 2 Fig. 1(e)F6 (cid:0) , , (cid:1) ( , ,
0) 14 2 Fig. 1(f)F7 ( , , (cid:0) , , − (cid:1) , , (cid:0) , , − (cid:1) , , (cid:0) , , − (cid:1) (cid:0) , , (cid:1) (cid:0) , , − (cid:1) (cid:0) , , − (cid:1) (cid:0) , , (cid:1) (cid:0) , , (cid:1) (cid:0) , , (cid:1) (cid:0) , , (cid:1) (cid:0) , , (cid:1) (cid:15) GUT ≤ − GeV. The VL masses vary between 0.25 TeV and10 TeV. In columns 2 and 3 transformation properties of both representations with respect to the SM gaugesymmetry group are given. It is understood that, if applicable, R F also encompasses its own complexconjugation. Columns 4 and 5 display number of VL fermions in each representation. In the last column wedirect the reader to a corresponding figure illustrating the allowed VL mass ranges and the iso-contours ofthe unification scale. completely novel solutions, characterized by either more than 12 copies of VL fermions inone of the representations, or VL masses larger than 5 TeV. Note that some of the successfulcombinations of representations allow for various choices of fermion numbers. This is, forexample, the case for the scenarios F1, F2 and F3, in which VL fields transforming as (cid:0) , , (cid:1) and (cid:0) , ¯2 , − (cid:1) of SU (3) C × SU (2) L × U (1) Y can show up in 12, 20 and 22 copies, while thosetransforming as (cid:0) , , (cid:1) and (cid:0) ¯6 , , − (cid:1) in 2 and 4 copies. Thus, only 7 combinations of R and R are really unique.Among 24 various representations listed in Eq. 3, only 9 can contribute to the successfulunification. These are adjoint representations of both SU (3) C and SU (2) L gauge groups,( , ,
0) and ( , , (cid:0) , , (cid:1) and (cid:0) , , (cid:1) ; fundamental representations of the SM right-handed up and down quarks, (cid:0) , , − (cid:1) , (cid:0) , , − (cid:1) ; and three exotic representations that arenot realized by the ordinary matter. Incidentally, the resemblance of the quantum numberscharacterizing VL fermions that allow for the PGU to the ones of the SM particles canhave important phenomenological consequences once the Yukawa-driven mixing with theSM fermions is allowed.It turns out that whether the gauge coupling unification is possible in a given modelhinges strongly on hierarchy among the VL fermion masses. To illustrate this dependence,in Fig. 1 and Fig. 2 we present a distribution of the mismatch parameter (cid:15) GUT as a functionof M and M for all 13 scenarios summarized in Table 3. In red the region of the PGUis indicated, which will be of main interest for our further phenomenological analysis. Asan additional information, we show in different shades of yellow departure from the precise7 a) F1 (b) F2 (c) F3(d) F4 (e) F5 (f) F6(g) F7 (h) F8 (i) F9Figure 1: Distribution of the mismatch parameter (cid:15) GUT as a function of M and M for the scenarios F1 - F9of Table 3. Red area corresponds to the PGU ( (cid:15) GUT ≤ (cid:15) GUT . Isocontours of theunification scale (in GeV) are indicated as dashed black curves. a) F10 (b) F11(c) F12 (d) F13Figure 2: Distribution of the mismatch parameter (cid:15) GUT as a function of M and M for the scenarios F10- F13 of Table 3. Red area corresponds to the PGU ( (cid:15) GUT ≤ (cid:15) GUT . Isocontoursof the unification scale (in GeV) are indicated as dashed black curves. unification condition as quantified by the increasing values of (cid:15)
GUT . The width of the colorbands can give one some idea on how easy the unification is, or, in other words, what is therequired degree of fine-tuning among the mass parameters.The successful PGU scenarios can be divided in three distinctive categories, dependingon the required mass hierarchy among the VL fermions. We will be referring to them laterusing the following labels:H0 : M ∼ M , scenarios F1 , F3 , F4 , (6)H1 : M (cid:29) M , scenarios F6 , F8 , F9 , F11 , F12 , H2 : M (cid:28) M , scenarios F2 , F5 , F7 , F10 , F13 . In Sec. 4 we will demonstrate that the mass hierarchy characterizing a given scenario iscrucial for the way the scenario can be tested experimentally.Another quantity that significantly differentiates among the VL combinations listed in9 ow Medium HighH0
F1, F4 F3F6, F11 H1 F9 F12 F8F2, F5 H2 F7, F10 F13
Table 4: Properties of the PGU scenarios in terms of the VL mass hierarchy and the unification scale andtheir susceptibility to various experimental search strategies. Light blue indicates the models that are testedby the proton decay measurements. Those highlighted in light green can be subject to color searches: R -hadrons and running of the strong gauge coupling. Light red corresponds to the scenarios tested throughEW interaction in lepton-like HSCP searches and EW precision tests. Table 3 is the unification scale. In Fig. 1 and Fig. 2 its isocontours are indicated as dashedblack curves. Depending on the order of magnitude of M GUT , we divide our scenarios inthree categories:Low : M GUT (cid:39) GeV , scenarios F1 , F4 , F9 , (7)Medium : M GUT ∼ GeV , scenarios F2 , F5 , F6 , F7 , F10 , F11 , F12 , High : M GUT ∼ GeV , scenarios F3 , F8 , F13 . We summarize the characteristics of the PGU scenarios in terms of the VL mass hierarchyand the unification scale in Table 4. The color code refers to experimental techniques thatcan be employed in order to test the available parameter space of a model. We will discussthem in details in Sec. 4.Before closing this section, we would like to comment on the fate of the identified PGUscenarios when the VL masses are pushed to energies much higher than 10 TeV. To analyzethis issue, we repeated the numerical procedure of Sec. 3 extending the scanning ranges of M and M up to 10 TeV. We found that all the models listed in Table 3, except F1 andF4, remain valid at higher energies, as could be anticipated from the shape of the red areas inFig. 1 and Fig. 2. Additionally, several new combinations of two VL representations becomeavailable.The high energy behavior of the PGU scenarios may seem somehow discouraging, asin principle one may put the VL fermions well above the reach of any existing colliderexperiment and still achieve the gauge coupling unification. There is, however, one importantremark to be made. The main factor that decides whether a given scenario is accepted or not,is the value of the unification scale, which we require to stay in the range 10 − GeV.When the mass of VL fermions increases, the unification scale decreases, as confirmed bythe shape of M GUT isocontours in Fig. 1 and Fig. 2. As a consequence, experimental boundsfrom the proton decay measurements may at some point come into play. In the next sectionwe will show that this is indeed the case and that the allowed parameter space of the PGUscenarios is limited from above. 10
Experimental tests of the PGU scenarios
In the previous section we identified all possible combinations of two VL fermion represen-tations that allow for precise unification of the three SM gauge couplings. In the followingwe will focus on phenomenological properties of the PGU scenarios and discuss in detailsvarious experimental ways of testing the available parameter space. We remind the readerthat we assume negligible Yukawa couplings among the BSM sector and the SM, thereforeour VL fermions can only be produced via gauge interactions. Experimental signatures ofthe models with a large number of VL fermions have been discussed by one of us in Ref. [40]and we follow closely its approach. We will demonstrate the complementarity among thebounds provided by various experimental searches, resulting from the fact that each of themaim at constraining particular sets of color and electroweak quantum numbers. In combina-tion with the specific mass hierarchies required by the PGU (see Table 4), it will allow us toderive strong lower bounds on the VL masses.
We assume that at the unification scale the symmetry group SU (3) C × SU (2) L × U (1) Y is embedded into a larger GUT group. Since in the unified framework the SM quarks andleptons belong to the same GUT multiplets, interactions are generated, mediated by heavygauge bosons, that violate both the baryon and lepton number conservation. Proton decay isthen a generic prediction of such scenarios. In non-SUSY models the dominant contributionto the proton decay width comes from dimension-6 gauge operators of the common structure QQQL . The exact form of these countributions highly depends on the realization of the GUTsymmetry. However, a rough estimation of the proton lifetime can be made [41] τ p = (cid:18) πg (cid:19) (cid:18) M GUT
GeV (cid:19) × . × − years , (8)as a simple function of the unification scale M GUT and the value of the unified gauge coupling g GUT .Proton decay has been experimentally searched for since the early 1990s by Super-Kamiokande (SK) underground water Cherenkov detector. The strongest lower bound onthe proton lifetime is set by the decay channel p → e + π and reads τ p > . × years [42].Hyper-Kamiokande (HK), a next generation machine, will be able to extend the limit by atleast one order of magnitude, up to ∼ × years [43].The present (solid blue line) and projected (dashed blue line) limits from the protondecay as a function of the VL masses M and M are shown in Figs. 3 and 4. The shadedblue area above the lines is disfavored. Scenarios F1 and F4 are already entirely excluded bythe SK measurements. Scenario F9, on the other hand, is going to be entirely tested by HK.Scenarios that fall within the reach of the current proton decay experiments are also markedin Table 4 in light blue. As expected, all of them belong to the category “low unificationscale”.Proton decay provides a unique experimental way of testing the PGU scenarios char-acterized by the BSM sector at the energy scales far above the reach of any present-daycollider experiment. In fact, for all but one models from Table 3 it provides upper bounds on11 roton decay Running g R -hadrons HSCP EWPO SummaryModel M max1 M max2 M min1 M min2 M min1 M min2 M min1 M min2 M min1 M min2 plotF1 Excluded - 0.7 - 1.8 0.8 - 1.7 - Fig. 3(a)F2 25 180 - 1.1 - 1.8 0.8 (6.0) 2.0 - Fig. 3(b)F3 350 200 - 1.1 (2.2) 1.8 0.8 - 2.2 - Fig. 3(c)F4 Excluded - 0.4 (1.2) 2.0 0.8 - 1.2 - Fig. 3(d)F5 10 50 - 0.8 - 2.0 0.8 (3.0) 1.5 - Fig. 3(e)F6 500 50 (9.0) 0.8 ( >
10) 2.0 0.8 - 1.7 - Fig. 3(f)F7 20 100 - 0.5 - 1.7 1.1 - 1.2 - Fig. 3(g)F8 2 × × (3.0) 0.8 (6.0) 1.7 1.1 - 1.2 - Fig. 3(h)F9 Excluded HK (4.5) 0.7 ( >
10) 1.8 1.1 - 1.5 - Fig. 3(i)F10 250 1000 - 1.1 - 1.8 1.2 (3.0) 2.0 - Fig. 4(a)F11 600 200 0.2 0.2 (5.0) 1.8 - - 0.5 1.0 Fig. 4(b)F12 6 ×
400 0.2 0.6 (4.0) 1.8 - - 0.7 1.5 Fig. 4(c)F13 - 2 × >
10) - - 0.7 1.7 Fig. 4(d)Table 5: Exclusion bounds on the VL masses M and M provided by different experiments (all in TeV).In column 2 we indicate the models that are exluded by the measurement of the proton lifetime by Super-Kamiokande [42], as well as the upper bounds provided by the projected Hyper-Kamiokande measure-ment [43]. In column 3 limits from the running strong coupling constant measurement by CMS [44] areshown. The numbers in parentheses indicate indirect limits whose derivation is described in the text. Incolumns 4 and 5 bounds from the ATLAS 13 TeV HSCP searches [45] are presented, for colored and non-colored particles, respectively. 100 TeV projections for the EWP tests [46] are shown in column 6. the allowed VL masses. We report them in the second column of Table 5 for the projectedreach from HK (the corresponding current bounds from SK are approximately one orderof magnitude weaker). One can see that for several scenarios that belong to the mediumGUT-scale category the upper bounds on VL masses are of the order of “only” several-tensTeV. This feature opens up an exiting possibility of entirely probing those scenarios in the(however distant) future. The RG running of the strong gauge coupling constant has been tested experimentally upto the energies of around 1 . √ s = 8 TeV with an integrated luminosityof 19 . / fb by the CMS Collaboration [44]. The value of the running coupling is extractedfrom the data as a function of the energy scale at which it is evaluated. The measurementis consistent with the predictions of the SM and as such poses a constraint on the minimalmass of any exotic colored particle.In the third column of Table 5 we summarize the lower bounds on the VL masses M and M in each PGU scenario. The same limits are also depicted in Figs. 3 and 4 as dark greensolid lines. Obviously, only the representations that transform non-trivially under SU (3) C can be directly constrained by the data.There is, however, an interesting observation to be made. In the scenarios with the mass12ierarchy H1, characterized by the color VL fermions much lighter than the non-colored ones,the running strong coupling allows one to indirectly put very strong lower bounds on massesof the fermions that are SU (3) C singlets and would be otherwise not affected by the CMSmeasurement. We indicate them in Table 5 as numbers in parentheses. As an example, let usconsider scenario F6, whose parameter space is subject to various experimental constraintspresented in Fig. 3(f). The direct lower bound on M from the running of g reads in thiscase 0 . M (cid:29) M . As a result, it is almost entirely probed by CMSand an indirect bound on the mass M can be derived, which reads in this case M ∼ > In the absence of Yukawa interactions with the SM quarks and leptons, the VL fermionsare stable and can be experimentally looked for at colliders through heavy stable chargedparticle (HSCP) searches. Dedicated analyses performed both by the ATLAS and CMScollaborations utilize observables related to the ionization energy loss ( dE/dx ) and time offlight (ToF), which allow to distinguish massive and non-relativistic HSCPs from the lightSM particles traveling with velocities close to the speed of light.Two categories of signals are usually considered, depending on the type of charges carriedby HSCP: one that consists of particles interacting strongly, and another in which HSCPsare lepton-like color singlets. The two would differ both by the production mechanism atthe LHC and by the size of the production cross section. In the following we will discussthem separately. Let us first assume that a heavy stable particle can interact strongly. If the lifetime of sucha colored HSCP is longer than typical hadronization time scale, it can form colorless QCDbound states with the SM quarks and gluons, the so-called R -hadrons.The most recent ToF and dE/dx based analyses have been performed by ATLAS using adata sample corresponding to 36 fb − of proton-proton collisions at √ s = 13 TeV [45], andby CMS using 2.5 fb − of data at the same energy [47]. Since in both cases no significantdeviations from the expected SM background have been observed, a model-independent 95%confidence level (C.L.) upper bound on the R -hadron production cross section can be derived. The presence of a stable charged particle at cosmological scales may be problematic from the point ofview of dark matter properties. A way out is to introduce Yukawa interaction with the SM, small enoughnot to affect the RG running but large enough to allow the charged particle to decay. Note, however, thatin the case of representations ( , , , , ) and ( , , ) it is not possible to construct a decay operatorwith the SM matter. One would need to introduce, for example, additional scalars charged under SU (3) C .This provides a motivation to extend the current analyses in the future to include the scalar fields. a) Scenario F1 (b) Scenario F2 (c) Scenario F3(d) Scenario F4 (e) Scenario F5 (f) Scenario F6(g) Scenario F7 (h) Scenario F8 (i) Scenario F9Figure 3: Summary of experimental bounds on VL masses M and M for scenarios F1-F9. In gray the PGUregion is indicated. The area below and left to the solid green line is excluded by the measurement of therunning strong coupling constant by CMS [44]. The limits from the 13 TeV ATLAS R -hadrons search [45]are indicated as a dashed green line. The corresponding lepton-like HSCP search excludes the area left to thedashed red line. 100 TeV projections for the EWP tests [46] are depicted as red dotted lines. Blue shadedregion marks the exlusion by the proton decay measurement at Super-Kamiokande [42]. A projected reachof Hyper-Kamiokande [43] is shown as a blue dashed line. a) Scenario F10 (b) Scenario F11(c) Scenario F12 (d) Scenario F13Figure 4: The same as Fig. 3 but for scenarios F10-F13. Such a result can then be translated into a lower bound on the BSM fermion mass withinan arbitrary framework. An example usually considered by the collaborations is the gluino,the SUSY partner of the gluon and a benchmark for a BSM fermion with the SM charges( , , SU (3) C quantum numbers. We calculated the pp → ¯ QQ cross section at the leading order (LO) using MadGraph5 aMC@NLO , and then rescaled it withthe k -factor of 2 [48] to account for higher-order QCD corrections and to reproduce thecross section quoted in [45] for gluino pair-production. We then compared the result withthe observed exclusion limit on the gluino derived by ATLAS. The corresponding exclusionbounds applied to parameters M and M are indicated in Figs. 3 and 4 as dashed greenlines. We also summarize them in the fourth column of Table 5. The cross section changes by over three orders of magnitude over the VL mass range considered in [45].The resulting exclusion bound is, therefore, very mildly sensitive to higher-order order corrections to thecross section. R -hadron searches probe the parameter space in the same directionas the measurement of the running strong coupling constant from Sec. 4.2. Therefore, themass parameter of the colored VL representations is constrained. As before, indirect boundson the non-colored representations can be derived in the case of the type H1 mass hierarchy.The effect is particularly visible for scenarios F6 and F9, which turn out to be excluded bythe R -hadron searches up to at least 10 TeV. In several other scenarios M receives a stronglower bound as well, which we indicate in Table 5 as a number in parentheses.Note that in the presence of non-zero Yukawa interactions the VL colored fermions maydecay before a bound state is formed. In such a case, the limits from the R -hadron searcheswill no longer apply. For that reason we indicate them in Figs. 3 and 4 with dashed lines,as contrasted with the running strong coupling constraints that are model-independent oncethe SU (3) C charges are fixed. If a charged HSCP does not interact hadronically, it will be produced through Drell-Yan(DY) processes and will predominantly lose energy via ionization inside the detector. Inthe analysis [45] with 36 fb − of data ATLAS interpreted the model-independent results ina benchmark model that assummes DY production of charginos. The corresponding lowerlimits on the HSCP mass reads 1090 GeV. In the analogous study by CMS [47] based on2.5 fb − dataset, bounds on the mass of a generic lepton-like fermions with a unit electriccharge were derived at 550 GeV.To set lower bounds on the masses of VL fermions that are SU (3) C singlets, we usedthe chargino-dedicated search by ATLAS. The LO production cross sections were calculatedwith MadGraph5 aMC@NLO , but no k -factor was added. The corresponding exclusion boundsin the ( M , M ) plane are depicted in Figs. 3 and 4 as dashed red lines. We also summarizethem in the fifth column of Table 5.The limits from the lepton-like HSCP searches allow to probe the scenarios with themass hierarchy H2, in which the non-colored VL fermions are lighter than the colored ones.Enhanced susceptibility to electroweak searches is marked in Table 4 in light red. In general,the limits are significantly weaker than the corresponding color-based bounds due to thelower production cross section. On the other hand, indirect lower bounds on the mass M can be derived in scenarios F2, F5, F10, which turns out to be much stronger than the directbounds from the R -hadron searches or the running strong coupling measurement. A complementary way to study properties of the VL fermions is to look at the processesbelow the M VL mass threshold. Such an approach can result particularly important if VLfermions are too heavy to be directly produced in the colliders, or not long-lived enough fordedicated HSCP searches to be effective. In this regard, high-energy measurements of DYprocesses at the LHC offer a promising way to indirectly look for VL fermions by testingdepartures from the SM predictions in electroweak precision (EWP) observables [46].In the VL extensions of the SM considered in this paper, the BSM contributions canmanifest themselves in two oblique parameters [49, 50] that are sensitive to the presence of16 cenario Current status Experimental testF1 Excluded proton decayF2 M > . M > . M > . M > . R -hadronsF4 Excluded proton decayF5 M > . M > . R -hadronsF7 M > . M > . R -hadrons, HSCPF8 M > . M > . R -hadronsF9 To be tested by HK proton decayF10 M > . M > . M > . M > . R -hadronsF12 M > . M > . R -hadronsF13 Excluded up to 10 TeV R -hadronsTable 6: Summary of current experimental status of the successful PGU scenarios. states charged under the EW gauge symmetry, W and Y . The experimental bounds on W and Y are derived from the measuerements of charged and neutral currents DY at hadroncolliders. The VL fermion contributions to the parameters W and Y are directly related tothe corresponding beta functions and given by [51] W, Y = g , π m W M × ∆ B BSM2 , . (9)Here, ∆ B BSM2 , denote the pure BSM contributions to the one-loop coefficients B , , Eq. 19and Eq. 20, respectively.The most up-to-date EWP experimental limits have been presented in [46], includingdata from LEP [52] and LHC 8 TeV measurements by ATLAS [53] and CMS [54]. Wechecked that they do not provide any bounds on the parameter space of the PGU scenariosunder study. However, since the effects of W and Y on DY processes grow with energy,the present experimental bounds can be significantly improved at the future colliders, byroughy two orders of magnitude at the projected 100 TeV machine [46]. The correspondingprojections with 3 ab − are depicted in Figs. 3 and 4 as dotted red lines. They are alsosummarized in the sixth column of Table 5.As expected from the size of the corresponding gauge couplings and group-theoreticalfactors, the constraints on W are stronger than those on Y . Therefore, the projected EWPbounds are particularly powerful for VL representations with the non-trivial SU (2) L charges.As a consequence, in most cases it is the mass parameter M that can be directly constrainedby the EWP tests. The only exceptions are scenarios F11-F13, in which both VL representa-tions can be constrained. Note also that in those three cases the projected EWP bounds canactually be competitive with the present day measurement of the running strong gauge cou-pling constant. Finally, it is worth to stress that similarly to what we observed in Sec. 4.3.2for the lepton-like HSCP searches, in the PGU scenarios with the mass hierarchy H2 indirectlower bounds on the mass of the SU (2) L singlet representations M can be obtained.17o summarize the findings of this section, we collect in Table 6 information about thecurrent experimental status of the successful PGU scenarios. In this regards, we can dividethem into three distinct categories. The first one encompasses scenarios F1 and F4, whichare already excluded by the proton decay measurements, and scenario F9, which will beentirely tested by Hyper-Kamiokande. The second category corresponds to those scenarios(F6 and F13) that are excluded up to at least 10 TeV, but which become allowed once higherVL masses are considered. The remaining eight scenarios feature the parameter space thatstill evades experimental bounds for VL masses in the multi- TeV regime. It should be noted,however, that some of them (F2, F5, F8 and F11) could be in the future and with more dataentirely tested within the considered mass range by the HSCP searches, while two others(F7 and F12) can be tested for the most part. Scenarios F3, and F10, on the other hand,will remain more challenging to explore. In light of null results from New Physics searches at the LHC, we look at unification ofthe gauge couplings as a model-building principle and classify possible SM extensions thatfeature this property.As a first step, we considered in this study extensions of the SM with two distinctrepresentations of VL fermions. We analyzed all their possible combinations with the numberof fermions in each representation limited only by perturbativity of the gauge couplings atthe unification scale. We found 13 different combinations of two representations that allowfor precise gauge unification at energies higher than 10 GeV, and for VL masses in therange 0 . −
10 TeV.Interdependence between types of spectra required by the unification condition and theirsusceptibility to experimental tests is the main characteristics of successful PGU scenarios.We showed that the effectiveness of a given search in probing the allowed parameter spaceof a model is directly related to its two features: mass hierarchy among VL fermions andthe value of the unification scale. Scenarios in which the colored fermions are much lighterthan the non-colored ones may be almost entirely tested by the measurement of the runningstrong coupling and by the LHC R -hadron searches. And vice versa , if non-colored fermionsare much lighter, HSCP searches and EW precision tests become very effective. On theother hand, scenarios in which both VL masses are of the same order remain beyond thereach of present-days colliders. In this case, however, null outcome from the proton decayexperiments allows to exclude those models that feature the low unification scale.The results presented in this study clearly highlight the importance of combining differentexperimental strategies in order to derive the most robust constraints on the PGU parameterspace. In this regard, proton decay measurements play a particular role, as they offer theonly mean of probing the VL spectra above the multi- TeV regime. There is also a greatpotential in the direct HSCP searches at the LHC. We hope that our results will proveuseful for experimental collaborations in choosing benchmark BSM scenarios for their futureanalyses.The current study can be extended in different directions. First of all, one may considermore complex (and more realistic) BSM scenarios, featuring for example more than two VL18epresentations or extra scalars. Secondly, the effects of non-gauge interactions (Yukawa andscalar types) should be taken into account, as they are bound to affect the phenomenologyof PGU scenarios. After all, the desert seems like an interesting place to explore. ACKNOWLEDGMENTS
We would like to thank Enrico Sessolo for his comments on the manuscript. The use of theCIS computer cluster at the National Centre for Nuclear Research in Warsaw is gratefullyacknowledged. This work is supported by the National Science Centre (Poland) under theresearch Grant No. 2017/26/E/ST2/00470.
A Group invariants and beta functions
General two-loop beta functions for a system of gauge couplings g i of a direct-product sym-metry group G i × · · · read [38] β i = dgd ln µ = g i (4 π ) (cid:20) − C ( G i ) + 23 S ( R F i ) + 13 S ( R Si ) (cid:21) (10)+ g i (4 π ) (cid:20) − C ( G i ) + (cid:16) C ( R F i ) + 103 C ( G i ) (cid:17) S ( R F i ) + (cid:16) C ( R Si ) + 23 C ( G i ) (cid:17) S ( R Si )+ k (cid:88) j =1 g j (cid:16) C ( R F j ) S ( R F i ) + 4 C ( R Sj ) S ( R Si ) (cid:17)(cid:35) , where G i , R F i and R Si denote contributions from gauge bosons, Weyl fermions, and complexscalars respectively. C ( R ) is a quadratic Casimir invariant, S ( R ) a Dynkin index of arepresentation R , and the sum is meant in both S ( R F i ) and S ( R Si ) over all fermion andscalar representations transforming nontrivially under G i .The quadratic Casimir operator for the representation R of a symmetry group G is definedas C ( R ) δ ij = ( t A t A ) ij = d (cid:88) A =1 t A t A , (11)where t A are the generators of G in the representation R . The Dynkin index of a represen-tation R is instead given by S ( R ) δ A B = Tr (cid:8) t A t B (cid:9) . (12)The two are related through the dimensions of the representation R , d ( R ), and of the adjoint, d (Adj), S ( R ) d (Adj) = C ( R ) d ( R ) . (13)It is convenient to parameterize the quadratic Casimir operator, the Dynkin index, andthe dimension of the representation through the weights ( p, q ) for irreducible SU (3) repre-sentations R , and, similarly, through the highest weight (cid:96) for SU (2) representations R ,19 ( R ) = ( p + 1)( q + 1)( p + q + 2) ,C ( R ) = p + q + ( p + q + pq ) , with p, q = 0 , · · · ,d ( R ) = 2 (cid:96) + 1 ,C ( R ) = (cid:96) ( (cid:96) + 1) , with (cid:96) = 0 , , · · · . (14)The two-loop beta functions for the SM augmented with N F fermions in the representa-tion ( R F , R F , Y F ) are straightforwardly derived from10, and read β = g (4 π ) B + g (4 π ) (cid:0) C g + C g + C g (cid:1) , (15) β = g (4 π ) B + g (4 π ) (cid:0) C g + C g + C g (cid:1) , (16) β = g (4 π ) B + g (4 π ) (cid:0) C g + C g + C g (cid:1) , (17)with the one-loop coefficients determined as B = − N F S ( R F ) d ( R F ) , (18) B = −
196 + 23 N F S ( R F ) d ( R F ) , (19) B = 4110 + 25 N F d ( R F ) d ( R F ) Y F . (20)The two-loop coefficients are given by C = −
26 + N F S ( R F ) d ( R F ) (cid:0) C ( R F ) + 10 (cid:1) , (21) C = 92 + 2 N F S ( R F ) C ( R F ) d ( R F ) , (22) C = 1110 + 65 N F S ( R F ) d ( R F ) Y F , (23) C = 12 + 2 N F S ( R F ) C ( R F ) d ( R F ) , (24) C = 356 + N F S ( R F ) d ( R F ) (cid:0) C ( R F ) + 203 (cid:1) , (25) C = 910 + 65 N F S ( R F ) d ( R F ) Y F , (26) C = 444 + 65 N F C ( R F ) d ( R F ) d ( R F ) Y F , (27) C = 2710 + 65 N F C ( R F ) d ( R F ) d ( R F ) Y F , (28)20 = 19950 + 1825 N F d ( R F ) d ( R F ) Y F . (29) B Decomposition of the irreducible SU (5) representa-tions In this Appendix we collected the branching rules for the embedding SU (5) ⊃ SU (3) × SU (2) × U (1) [55], = (cid:0) , , (cid:1) ⊕ (cid:0) , , − (cid:1) , (30) = ( , , ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ (cid:0) , , (cid:1) , = ( , , ⊕ (cid:0) , , (cid:1) ⊕ (cid:0) , , − (cid:1) , = ( , , ⊕ ( , , ⊕ ( , , ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) ¯3 , , (cid:1) , = (cid:0) , , − (cid:1) ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ (cid:0) ¯6 , , (cid:1) ⊕ ( ¯10 , , , = (cid:0) , , − (cid:1) ⊕ (cid:0) , , (cid:1) ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ ( , , ⊕ (cid:0) ¯6 , , (cid:1) , = (cid:0) , , (cid:1) ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) ¯3 , , (cid:1) ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ (cid:0) ¯6 , , − (cid:1) ⊕ (cid:0) , , (cid:1) , = ( , , − ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ (cid:0) ¯6 , , − (cid:1) ⊕ (cid:0) , , (cid:1) ⊕ (cid:0) , , (cid:1) , = (cid:0) , , (cid:1) ⊕ (cid:0) , , (cid:1) ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) ¯3 , , (cid:1) ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) , , (cid:1) ⊕ (cid:0) , , − (cid:1) , (cid:48) = ( , , − ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ (cid:0) ¯6 , , − (cid:1) ⊕ (cid:0) , , (cid:1) ⊕ (cid:16) (cid:48) , , (cid:17) , (cid:48) = ( , , ⊕ (cid:0) , , (cid:1) ⊕ (cid:0) , , − (cid:1) ⊕ (cid:0) ¯3 , , − (cid:1) ⊕ (cid:0) ¯3 , , (cid:1) ⊕ (cid:0) ¯6 , , − (cid:1) ⊕ (cid:0) , , (cid:1) ⊕ ( , , ⊕ ( , , . References [1] H. Georgi and S. L. Glashow, “Unity of All Elementary Particle Forces,”
Phys. Rev.Lett. (1974) 438–441.[2] J. C. Pati and A. Salam, “Lepton Number as the Fourth Color,” Phys. Rev.
D10 (1974)275–289. [Erratum: Phys. Rev.D11,703(1975)].[3] R. N. Mohapatra and J. C. Pati, “Left-Right Gauge Symmetry and an IsoconjugateModel of CP Violation,”
Phys. Rev.
D11 (1975) 566–571.[4] H. Fritzsch and P. Minkowski, “Unified Interactions of Leptons and Hadrons,”
AnnalsPhys. (1975) 193–266.[5] H. Georgi, “The State of the ArtGauge Theories,” AIP Conf. Proc. (1975) 575–582.[6] G. P. Cook, K. T. Mahanthappa, and M. A. Sher, “Effects of Heavy Colored HiggsScalars on SU(5) Predictions,” Phys. Lett. (1980) 398–400.217] V. V. Dixit and M. Sher, “The Futility of High Precision SO(10) Calculations,”
Phys.Rev.
D40 (1989) 3765.[8] P. Langacker and N. Polonsky, “Uncertainties in coupling constant unification,”
Phys.Rev.
D47 (1993) 4028–4045, arXiv:hep-ph/9210235 [hep-ph] .[9] K. Hagiwara and Y. Yamada, “Grand unification threshold effects in supersymmetricSU(5) models,”
Phys. Rev. Lett. (1993) 709–712.[10] J. Schwichtenberg, “Gauge Coupling Unification without Supersymmetry,” Eur. Phys.J.
C79 no. 4, (2019) 351, arXiv:1808.10329 [hep-ph] .[11] Y. Achiman and B. Stech,
New Phenomena in Lepton-Hadron Physics . edited by D.E.C.Fries and J. Wess, Plenum, New York, 1979, p. 303.[12] A. de R´ujula, H. Georgi, and S. L. Glashow,
Fifth Workshop on Grand Unification .edited by K. Kang, H. Fried, and P. Frampton, World Scientific, Singapore, 1984, p. 88.[13] P. Fileviez Perez and S. Ohmer, “Low Scale Unification of Gauge Interactions,”
Phys.Rev.
D90 no. 3, (2014) 037701, arXiv:1405.1199 [hep-ph] .[14] T. G. Rizzo, “Desert guts and new light degrees of freedom,”
Phys. Rev.
D45 (1992)3903–3905.[15] B. Zhang and H.-Q. Zheng, “Top quark, heavy fermions and the composite Higgs boson,”
Commun. Theor. Phys. (2001) 162–166, arXiv:hep-ph/0003065 [hep-ph] .[16] D. Choudhury, T. M. P. Tait, and C. E. M. Wagner, “Beautiful mirrors and precisionelectroweak data,” Phys. Rev.
D65 (2002) 053002, arXiv:hep-ph/0109097 [hep-ph] .[17] L.-F. Li and F. Wu, “Coupling constant unification in extensions of standard model,”
Int. J. Mod. Phys.
A19 (2004) 3217–3224, arXiv:hep-ph/0304238 [hep-ph] .[18] D. E. Morrissey and C. E. M. Wagner, “Beautiful mirrors, unification of couplingsand collider phenomenology,”
Phys. Rev.
D69 (2004) 053001, arXiv:hep-ph/0308001[hep-ph] .[19] G. F. Giudice and A. Romanino, “Split supersymmetry,”
Nucl. Phys.
B699 (2004)65–89, arXiv:hep-ph/0406088 [hep-ph] . [Erratum: Nucl. Phys.B706,487(2005)].[20] I. Dorsner and P. Fileviez Perez, “Unification without supersymmetry: Neu-trino mass, proton decay and light leptoquarks,”
Nucl. Phys.
B723 (2005) 53–76, arXiv:hep-ph/0504276 [hep-ph] .[21] D. Emmanuel-Costa and R. Gonzalez Felipe, “Minimal string-scale unification of gaugecouplings,”
Phys. Lett.
B623 (2005) 111–118, arXiv:hep-ph/0505257 [hep-ph] .[22] R. Shrock, “Variants of the Standard Model with Electroweak-Singlet Quarks,”
Phys.Rev.
D78 (2008) 076009, arXiv:0809.0087 [hep-ph] .2223] I. Gogoladze, B. He, and Q. Shafi, “New Fermions at the LHC and Mass of the HiggsBoson,”
Phys. Lett.
B690 (2010) 495–500, arXiv:1004.4217 [hep-ph] .[24] R. Dermisek, “Insensitive Unification of Gauge Couplings,”
Phys. Lett.
B713 (2012)469–472, arXiv:1204.6533 [hep-ph] .[25] R. Dermisek, “Unification of gauge couplings in the standard model with extra vectorlikefamilies,”
Phys. Rev.
D87 no. 5, (2013) 055008, arXiv:1212.3035 [hep-ph] .[26] I. Dorsner, S. Fajfer, and I. Mustac, “Light vector-like fermions in a minimal SU(5)setup,”
Phys. Rev.
D89 no. 11, (2014) 115004, arXiv:1401.6870 [hep-ph] .[27] M.-L. Xiao and J.-H. Yu, “Stabilizing electroweak vacuum in a vectorlike fermionmodel,”
Phys. Rev.
D90 no. 1, (2014) 014007, arXiv:1404.0681 [hep-ph] . [Adden-dum: Phys. Rev.D90,no.1,019901(2014)].[28] B. Bhattacherjee, P. Byakti, A. Kushwaha, and S. K. Vempati, “Unification with Vector-like fermions and signals at LHC,”
JHEP (2018) 090, arXiv:1702.06417 [hep-ph] .[29] N. Arkani-Hamed and S. Dimopoulos, “Supersymmetric unification without low en-ergy supersymmetry and signatures for fine-tuning at the LHC,” JHEP (2005) 073, arXiv:hep-th/0405159 [hep-th] .[30] V. Barger, J. Jiang, P. Langacker, and T. Li, “String scale gauge coupling unificationwith vector-like exotics and non-canonical U(1)(Y) normalization,” Int. J. Mod. Phys.
A22 (2007) 6203–6218, arXiv:hep-ph/0612206 [hep-ph] .[31] V. Barger, N. G. Deshpande, J. Jiang, P. Langacker, and T. Li, “Implications of Canon-ical Gauge Coupling Unification in High-Scale Supersymmetry Breaking,”
Nucl. Phys.
B793 (2008) 307–325, arXiv:hep-ph/0701136 [hep-ph] .[32] L. Calibbi, L. Ferretti, A. Romanino, and R. Ziegler, “Gauge coupling unification, theGUT scale, and magic fields,”
Phys. Lett.
B672 (2009) 152–157, arXiv:0812.0342[hep-ph] .[33] I. Donkin and A. Hebecker, “Precision Gauge Unification from Extra Yukawa Cou-plings,”
JHEP (2010) 044, arXiv:1007.3990 [hep-ph] .[34] C. Liu and Z.-h. Zhao, “ θ and the Higgs mass from high scale supersymmetry,” Com-mun. Theor. Phys. (2013) 467–471, arXiv:1205.3849 [hep-ph] .[35] S. Zheng, “Effective Higgs Theories in Supersymmetric Grand Unification,” Eur. Phys.J.
C77 no. 9, (2017) 588, arXiv:1706.01071 [hep-ph] .[36] S. Zheng, “Minimal Vectorlike Model in Supersymmetric Unification,” arXiv:1904.10145 [hep-ph] .[37] D. Buttazzo, G. Degrassi, P. P. Giardino, G. F. Giudice, F. Sala, A. Salvio, andA. Strumia, “Investigating the near-criticality of the Higgs boson,”
JHEP (2013)089, arXiv:1307.3536 [hep-ph] . 2338] M. E. Machacek and M. T. Vaughn, “Two Loop Renormalization Group Equations in aGeneral Quantum Field Theory. 1. Wave Function Renormalization,” Nucl. Phys.
B222 (1983) 83–103.[39] A. Perez-Lorenzana and W. A. Ponce, “GUTs and string GUTs,”
Europhys. Lett. (2000) 296–301, arXiv:hep-ph/9911540 [hep-ph] .[40] A. D. Bond, G. Hiller, K. Kowalska, and D. F. Litim, “Directions for model buildingfrom asymptotic safety,” JHEP (2017) 004, arXiv:1702.01727 [hep-ph] .[41] P. Nath and P. Fileviez Perez, “Proton stability in grand unified theories, in strings andin branes,” Phys. Rept. (2007) 191–317, arXiv:hep-ph/0601023 [hep-ph] .[42]
Super-Kamiokande
Collaboration, K. Abe et al. , “Search for proton decay via p → e + π and p → µ + π in 0.31 megatonyears exposure of the Super-Kamiokandewater Cherenkov detector,” Phys. Rev.
D95 no. 1, (2017) 012004, arXiv:1610.03597[hep-ex] .[43]
Hyper-Kamiokande
Collaboration, K. Abe et al. , “Hyper-Kamiokande Design Re-port,” arXiv:1805.04163 [physics.ins-det] .[44]
CMS
Collaboration, V. Khachatryan et al. , “Measurement and QCD analysis of double-differential inclusive jet cross sections in pp collisions at √ s = 8 TeV and cross sectionratios to 2.76 and 7 TeV,” JHEP (2017) 156, arXiv:1609.05331 [hep-ex] .[45] ATLAS
Collaboration, M. Aaboud et al. , “Search for heavy charged long-lived particlesin the ATLAS detector in 36.1 fb − of proton-proton collision data at √ s = 13 TeV,” Phys. Rev.
D99 no. 9, (2019) 092007, arXiv:1902.01636 [hep-ex] .[46] M. Farina, G. Panico, D. Pappadopulo, J. T. Ruderman, R. Torre, and A. Wulzer,“Energy helps accuracy: electroweak precision tests at hadron colliders,”
Phys. Lett.
B772 (2017) 210–215, arXiv:1609.08157 [hep-ph] .[47]
CMS
Collaboration, V. Khachatryan et al. , “Search for long-lived charged particlesin proton-proton collisions at √ s = 13 TeV,” Phys. Rev.
D94 no. 11, (2016) 112004, arXiv:1609.08382 [hep-ex] .[48]
Particle Data Group
Collaboration, M. Tanabashi et al. , “Review of ParticlePhysics,”
Phys. Rev.
D98 no. 3, (2018) 030001.[49] R. Barbieri, A. Pomarol, R. Rattazzi, and A. Strumia, “Electroweak symmetry breakingafter LEP-1 and LEP-2,”
Nucl. Phys.
B703 (2004) 127–146, arXiv:hep-ph/0405040[hep-ph] .[50] G. Cacciapaglia, C. Csaki, G. Marandella, and A. Strumia, “The Minimal Set of Elec-troweak Precision Parameters,”
Phys. Rev.
D74 (2006) 033011, arXiv:hep-ph/0604111[hep-ph] . 2451] D. S. M. Alves, J. Galloway, J. T. Ruderman, and J. R. Walsh, “Running Elec-troweak Couplings as a Probe of New Physics,”
JHEP (2015) 007, arXiv:1410.6810[hep-ph] .[52] A. Falkowski and K. Mimouni, “Model independent constraints on four-lepton opera-tors,” JHEP (2016) 086, arXiv:1511.07434 [hep-ph] .[53] ATLAS
Collaboration, G. Aad et al. , “Measurement of the double-differential high-mass Drell-Yan cross section in pp collisions at √ s = 8 TeV with the ATLAS detector,” JHEP (2016) 009, arXiv:1606.01736 [hep-ex] .[54] CMS
Collaboration, V. Khachatryan et al. , “Measurements of differential and double-differential Drell-Yan cross sections in proton-proton collisions at 8 TeV,”
Eur. Phys.J.
C75 no. 4, (2015) 147, arXiv:1412.1115 [hep-ex] .[55] R. Slansky, “Group theory for unified model building,”
Physics Reports no. 1, (1981)1 – 128.