# GUT Physics in the Era of the LHC

^{1}TRIUMF Theory Group, Vancouver, BC, Canada^{2}Department of Physics, University of Oslo, Oslo, Norway^{3}ARC Centre of Excellence for Particle Physics at the Tera-Scale, School of Physics and Astronomy, Monash University, Melbourne, VIC, Australia^{4}Department of Physics and Astronomy, University College London, London, United Kingdom^{5}Max-Planck-Institut für Kernphysik, Heidelberg, Germany^{6}Department of Physics, University of Ljubljana, Ljubljana, Slovenia^{7}Jožef Stefan Institute, Ljubljana, Slovenia

Grand Unified Theories (GUTs) are one of the most interesting high-energy completions of the Standard Model, because they provide a rich, powerful and elegant group-theoretical framework able to resolve a variety of problems remaining in our current understanding of particle physics. They usually act as motivators for many low energy BSM theories, such as left-right symmetric or supersymmetric models, and they serve to fill the gap between the experimentally reachable low energies and the physics in the ultraviolet. In recent years, however, they have fallen slightly from the spotlight, in favor of “simplified” models with more specific phenomenological predictions. The aim of this review is to summarize the state of the art on GUTs and argue for their importance in modern physics. Recent advances in experiments permit to test the predictions of GUTs at different energy scales. First, as GUTs can play a role in the inflationary dynamics of the early Universe, their imprints could be found in the CMB observations by the Planck satellite. Remarkably enough, GUTs could manifest themselves also in terrestrial tests; several planned experiments aim to probe the proton stability and to establish order of magnitude higher bounds on its lifetime. Moreover, the predictions of specific GUT models could be tested even at the LHC thanks to its high energy reach, via searches for exotic states or additional contributions to flavor anomalies.

## 1. Introduction

The Standard Model (SM) [1–3] of particle physics is an incredible successful theory of subatomic physics. It describes the electroweak and strong interactions of fundamental particles with surprising accuracy up to the energy scales of modern day experiments. Further supported by the discovery of the Higgs boson [4, 5], it stands as one of the best evidences that symmetries and the mechanism of spontaneous symmetry breaking play a critical role on the Universe at the smallest scales [6–8].

In spite of its success at explaining with astonishing precision most of the observed phenomena, the SM cannot be the ultimate theory of particle physics. The Higgs quartic coupling in the SM becomes negative at scales ≳10^{10} GeV, rendering the vacuum state of the theory unstable at high energies [9]. New physics must then surface below or around that scale to stabilize the vacuum. Furthermore, there is a continuously increasing amount of observations that are in tension with the predictions of the SM. From the discovery of neutrino oscillations [10, 11] to the recent anomalies in the flavor sector [12–14], these phenomena cannot be explained with the SM alone and contributions from beyond the SM (BSM) physics may be required to accommodate them.

Grand Unified Theories (GUTs) [15–19] are well motivated extensions of the SM that can address several of its outstanding issues. As the SM does for electromagnetism and weak nuclear decays, GUTs exploit the power of symmetries and group theory to unify the electroweak and strong interactions into a single force. As can be noticed in Figure 1, the flow of the SM gauge couplings already hints at a possible unification at a high scale, thereby providing further motivation for GUTs as appealing BSM models.

The predicted unification of forces provides an explanation for the *ad hoc* nature of the SM charge assignments and their accidental anomaly cancellation [20, 21]. Through the introduction of new fields and symmetries, GUTs can resolve many of the issues of the SM: they can provide an explanation for the lightness of neutrino masses, as well as introducing additional contributions that can accommodate some of the observed flavor anomalies. In addition, GUTs can also live alongside other BSM models, such as Supersymmetry (SUSY), with both theories complementing each other and on the whole becoming a better candidate for a successful BSM theory [22].

Naturally living at high energies, it is expected that GUTs have some observable consequences for the cosmological evolution of the Universe, as they can play a role during the inflationary epoch and their phase transitions may be the source for matter-antimatter asymmetry or gravitational waves [23]. With or without SUSY, GUTs also make predictions that can be tested at low energy experiments such as particle colliders, which can search for new exotic states or deviations on flavor observables. Some of its high energy repercussions can also be probed with precision experiments, with nucleon decay measurements at the forefront. In short, Grand Unified Theories have observable effects in many fronts and can be probed by current and upcoming experiments in the near future.

Therefore, throughout this review we will focus on the description of GUT models and their observable consequences. We will introduce the basic concepts and summarize some of the modern research in GUTs. Out of all the possible observable probes of GUTs, we will focus on a subset of them. In the cosmological front we will outline the advances on inflationary GUTs, detail their observable gravitational wave signatures, from cosmic strings and phase transitions, and their relation with mechanisms for baryo and leptogenesis. The low energy front will cover collider searches for supersymmetry, leptoquarks and exotic states. Lastly, a number of precision tests of unification will be detailed, including nucleon decay, flavor observables and neutrinoless double beta decay.

As ultraviolet (UV) completions of the Standard Model that live at very high energies, GUTs are often closely related to theories of gravity, such as string theory. In fact, many unified theories arise as four-dimensional compactifications in some realization of superstring theory [24, 25]. However, throughout this review we will not concern ourselves with these string theory realizations. For a review on embedding GUTs in the heterotic string and outcomes of string compactification for unified theories see Raby [22].

## 2. Basics on GUT Models

Grand Unified Theories are extensions of the SM with larger symmetry groups. Strictly speaking, GUTs require that the unified group be a simple group, e.g., *SU*(5), *SO*(10), or *E*_{6}. However, here we use the term GUT more loosely, referring to any extension of the SM symmetries including product groups such as *SU*(5) × *U*(1) and *SU*(4) × *SU*(2) × *SU*(2). Candidate groups for a realistic GUT model must satisfy two conditions: they must contain the SM group as a subgroup and they must have complex representations that reproduce the chiral structure of the SM. Although groups with pseudoreal representations have been studied as candidates for unified theories, *E*_{7} [26], we will not consider them here.

### 2.1. *SU*(5)

The first appearance of a GUT in the literature dates back to 1974 when H. Georgi and S. Glashow proposed the unification of the SM gauge group into a simple group, *SU*(5) [15]. In their proposal all the left-handed fermions of a single generation fell into two representations of the group, $\overline{5}$ and **10**, in the following way

and the gauge and Higgs sector of the theory were embedded into the **24** and **5** representations, respectively. In addition to the SM Higgs boson present in the representation **5**, often a scalar $\overline{5}$ representation is also present, which contains a second *SU*(2)_{L} doublet, necessary for UV completions of two-Higgs doublet models (2HDM) [27].

The Georgi-Glashow (GG) model was the first attempt of a fully-unified model for particle physics, and it provided a neat explanation for the hypercharge quantization in the SM. The traceless hypercharge generator can be written as [21]

which when acting upon the representations of *SU*(5) results in the specific hypercharge assignments of the SM fields, i.e., for α = 1, *Y*(*Q*) = 1/6, *Y*(*L*) = −1/2, *Y*(*u*^{c}) = −2/3, *Y*(*d*^{c}) = 1/3 and *Y*(*e*^{c}) = 1. In unified theories one often uses the “GUT normalization” of the hypercharge, which corresponds simply to a rescaling of the charges and gauge couplings as ${Y}_{GUT}=\sqrt{3/5}\text{}Y$ and ${g}_{1}=\sqrt{5/3}{\text{}g}^{\prime}$ [28].

Spontaneous symmetry breaking of *SU*(5) happens when a scalar field in a non-trivial representation of the group acquires a vacuum expectation value (vev). The minimal representation of *SU*(5) that can achieve this goal while keeping the SM phase unbroken is the **24** [15, 29]. This vev provides a mass to the off-diagonal *SU*(5) gauge bosons while the SM gauge bosons remain massless.

By virtue of the unification into a single gauge group, the GG model requires strict unification of the SM gauge couplings, which is hinted at but not really achieved in the SM [28, 29], as can be seen in Figure 1, as well as that of Yukawa couplings for each of the two representations, a difficult task in its minimal version [30, 31].

The minimal *SU*(5) version suffers from other afflictions beyond the failed gauge and Yukawa unification mentioned above. One of these is the introduction of an artificial hierarchy, known as doublet-triplet splitting [32, 33], in the components to the Higgs representation **5**. The colored components must be quite heavy to avoid rapid proton decay whereas the uncolored components must be relatively light, for they correspond to the SM Higgs doublet, now know to have a mass of *m*_{h} = 125.18 GeV [34]. Solutions to this problem in several *SU*(5) models have been proposed, such as the “missing partner mechanism” [33, 35] or the “double missing partner mechanism” [36, 37].

Another case where the minimal *SU*(5) falls short is the lack of a mechanism for the generation of neutrino masses. Extended scalar sectors can generate neutrino masses in the type-I [38] and type-III [39] seesaw mechanisms^{1}, but the resulting theories are often non-renormalisable. Renormalisable *SU*(5) models can also be constructed where the neutrino masses are generated via a mixture of type-I and type-III seesaw [40] or the Zee mechanism [41, 42].

Worst of all, however, is the fact that the vanilla *SU*(5) model predicts rapid proton decay. The lifetime of the proton can be naively estimated as [30]

with *m*_{p} the mass of the proton and *M*_{X} the mass of the mediator field at the scale of unification. The apparent unification of gauge couplings happens at an energy scale μ ~ 10^{15} GeV, which gives a half life for the proton of the order of 10^{31} years, far below the experimental bound from the Super-Kamiokande experiment of 1.6 × 10^{34} years [43]. Particular choices of the Higgs sector of the *SU*(5) model, however, avoid this issue, rendering non-minimal *SU*(5) models viable candidates [44–46]. Furthermore, *SU*(5) models with vector-like fermions can be consistent with current limits and even predict an upper bound on the lifetime of proton decay [47].

One of the fundamental issues with GUT models, which remains as a concern today, is the gauge hierarchy problem. The large hierarchy between the mass scale of unification and the electroweak scale poses a problem since it causes large loop corrections to the Higgs mass [48]. Supersymmetry (SUSY) was proposed as a solution to this issue [49] and even acted as a motivation for unified theories since some of its minimal realizations, such as the MSSM, predicted the unification of gauge couplings, as can be seen in Figure 2.

Supersymmetric GUTs are in fact rather popular and have in some cases been proven to be more successful at representing nature than regular GUTs [50]. In SUSY *SU*(5) theories the scale of unification is typically larger than in non-supersymmetric models^{2}, around μ ~ 2 × 10^{16} GeV as can be seen in Figure 2. This has two advantageous consequences for these models: the larger mass scale for the gauge mediators imposes a further suppression on nucleon decay processes, consistent with experimental measurements [50]; and pushes the unification scale beyond the scale of inflation, which helps to dilute the magnetic monopoles naturally present in the theory [52]. Another issue in vanilla *SU*(5) models that can be addressed in its supersymmetric version is the doublet-triplet splitting, where the Higgs doublets are made light via cancellations of the superpotential parameters [32, 33, 53, 54].

### 2.2. Flipped *SU*(5)

An alternative solution to the issues of *SU*(5) models, without supersymmetry, is what is now known as *flipped* *SU*(5) [55, 56]. The flipped version differs from regular *SU*(5) in its gauge group, extended to *SU*(5) × *U*(1), and the manner in which the SM particle fields are embedded into representations of the group. In contrast to Equation (1), the matter representations in the flipped *SU*(5) model are

where ν^{c} labels the right-handed neutrino field, whose presence provides a mechanism for neutrino mass generation, which was absent in vanilla *SU*(5).

With these different embeddings of the SM fields, the hypercharge operator is no longer contained in *SU*(5), as in Equation (2), but it is rather a combination of diagonal generators of both *SU*(5) and *U*(1). With standard normalization the hypercharge operator can be written as [57]

where *T*_{24} is a diagonal generator of *SU*(5) and *X* the *U*(1) charge.

Due to the extended gauge sector and modified unified conditions of flipped *SU*(5), proton decay does not become an issue [57]. In addition, in flipped *SU*(5) magnetic monopoles cannot be created since the supergroup containing the charge operator is not a simple group [52, 57].

As was the case with regular *SU*(5) models, flipped *SU*(5) can be extended with the help of supersymmetry. The combination of solutions to the issues of the *SU*(5) model that both SUSY and flipped *SU*(5) offer makes SUSY flipped *SU*(5) one of the most popular GUTs in the literature [58, 59], in spite of not realizing a full unification of gauge couplings.

Flipped *SU*(5) models are also well motivated from their UV completions, since they can easily be derived naturally from weakly-coupled string theory. As opposed to vanilla *SU*(5), which undergoes symmetry breaking via a **24**-dimensional representation, the flipped *SU*(5) model does not require large dimensional representations, as it breaks via a **10**_{1}, and it is therefore easier to obtain from a manifold compactification of string theory [60, 61].

### 2.3. Pati-Salam and the Left-Right Symmetry

Around the same time that the *SU*(5) model was proposed, R. Pati and A. Salam suggested another extension of the SM [16]. They proposed a rearrangement of the SM fields into a different group configuration, *SU*(4)_{c} × *SU*(2)_{L} × *SU*(2)_{R}. Though not really a fully unified theory, it provided a partial unification of leptons and quarks into a single color group, *SU*(4)_{c}, while at the same time introducing another copy of *SU*(2) for the right-handed sector of the theory. This model automatically contains a right-handed neutrino, embedded into a *SU*(2)_{R} doublet with the right-handed charged lepton. Thus, the SM fields fall into two representations of the group in the following way

One of the major successes of the Pati-Salam (PS) model was being the first appearance of a left-right symmetric model, with a right-handed sector *SU*(2)_{R} and a heavy right-handed gauge boson *W*_{R} [62]. It was also the original proposal for the idea of quark-lepton complementarity. As an amalgamation of the two ideas, the PS group maximally contains the left-right symmetry group, *SU*(3)_{c} × *SU*(2)_{L} × *SU*(2)_{R} × *U*(1)_{B−L} [17, 63], as well as the quark-lepton unified group, *SU*(4)_{c} × *SU*(2)_{L} × *U*(1)_{R} [64, 65].

Left-right symmetric models, à* la* Pati-Salam or of the type *SU*(3)_{c} × *SU*(2)_{L} × *SU*(2)_{R} × *U*(1)_{B−L}, are fairly popular because they naturally include a right-handed neutrino and can generate light neutrino masses via some type of seesaw mechanism [66, 67]. Similar to PS, left-right symmetric (LR) models are not fully unified theories, yet they can be an intermediate step on the breaking chain of a PS model [68] or some other unified theory [69].

Symmetry breaking in the PS model can happen through a number of different paths, depending on the vev of the scalar fields in the theory. The most compelling paths preserve either the LR symmetry, with the LR group *SU*(3)_{c} × *SU*(2)_{L} × *SU*(2)_{R} × *U*(1)_{B−L} as an intermediate step, or quark-lepton unification, with *SU*(4)_{c} × *SU*(2)_{L} × *U*(1)_{R} as a subgroup [70]. Further symmetry breaking from the LR symmetry model happens when either a pair of *SU*(2) doublets (one left-handed and one right-handed) [62, 71], or a pair of *SU*(2) triplets (left and right-handed) [68] acquire a vev. In both PS and LR theories, the hypercharge operator is written as a linear combination of the diagonal *SU*(2)_{R} generator and the *B*−*L* charge [*U*(1)_{B−L} generator embedded in *SU*(4)_{c} in PS] as

As opposed to the case of *SU*(5) the proton is often stable in PS and LR models. This occurs because the gauge sector of the theory preserves *B* and *L* number independently and the only scalar fields that can mediate the transition are in antisymmetric representations, rarely seen in PS or LR models [68].

The addition of supersymmetry to PS and left-right symmetric models [72] is not as straightforward as with other GUT models. The simplest scenario with both SUSY and LR symmetry was shown to fail to achieve spontaneous symmetry breaking [73]. In order to circumvent this issue one must either add extra fields, such as a parity-odd singlet^{3} [74] or an extra Higgs field [75], or alternatively supplement the Lagrangian with non-renormalizable operators [76]. One of the main advantages of SUSY LR models, and the reason why so much effort is put on resolving the SSB issue, is that they naturally preserve *R*-parity, since *B* − *L* is a gauge symmetry of the theory, which forbids the dangerous baryon and lepton number violating operators that appear in the MSSM, thereby making the lightest SUSY particle stable [77].

### 2.4. *SO*(10)

Although the GG and PS models seem quite distinct in their approach to unification, they have a common ancestor. Both *SU*(5) × *U*(1) and *SU*(4) × *SU*(2) × *SU*(2) are maximal subgroups of another Lie group of larger dimension, *SO*(10). This was first realized by Fritzsch and Minkowski [18], and independently by Georgi [19], who proposed a model of unification with several intermediate steps. *SO*(10) has since been the most popular choice as a unification group, since it provides a vast display of options for field configurations and symmetry breaking patterns.

One of the many key features of *SO*(10) models is that they fully unify a generation of SM fermions into a single representation of the group. Thus, the 16 Weyl fermions, including right-handed neutrinos, can be embedded into the fundamental **16** representation of *SO*(10) as

Due to the transformation properties of the *SO*(10) group, the spinor representation **16** is a complex representation, thereby satisfying chirality as observed in the Standard Model. Additionally *SO*(10) is a “safe algebra” [78], it does not suffer from anomalies as, for example, the *SU*(5) case above, which makes model building in *SO*(10) easier for it does not rely on some specific field configurations to cancel the gauge anomalies [79].

Despite the large amount of *SO*(10) models in the literature, a common feature is that the gauge fields are embedded in the adjoint representation of the group, that is **45**,

which contains the SM gauge bosons as well as off-diagonal components which, as happened in *SU*(5), can mediate quark-lepton transitions, known as leptoquarks. The Yukawa sector in *SO*(10) models is often also quite recurrent because, at the renormalizable level, it can only be of the form [80]

where Y is the matrix of Yukawa couplings, *C*_{L} and *C*_{10} the charge conjugation matrices in the Poincaré and *SO*(10) groups, Γ_{i} the generators of *SO*(10) in the spinor representation and Φ^{i}, Φ^{ijk}, and Φ^{ijklm} are scalar fields in the **10, 120**, and $\overline{126}$ representations, respectively. The SM Higgs field is, therefore, some linear combination of these fields and hence the SM fermion masses directly follow from the Yukawa matrix **Y** and the vacuum expectation values of the scalar fields. The particular choice of the scalar sector is typically guided by the principle of minimality. While the minimal regular (non-SUSY) *SO*(10) model with SSB driven by the **45** and **126** Higgs representations has been revived and still represents a phenomenologically viable scenario [81–85], this is not the case of the minimal SUSY *SO*(10) model [86, 87].

Symmetry breaking in *SO*(10) models can occur through one of many different paths. Since both *SU*(5) × *U*(1) and *SU*(4) × *SU*(2) × *SU*(2) are subgroups of *SO*(10), they can be an intermediate step on the symmetry breaking path toward the Standard Model, as can be any of their respective subgroup [70, 88, 89]. Alternatively *SO*(10) can be broken directly to the SM group, without intermediate steps [90]. All the possible breaking paths from *SO*(10) can be seen in Figure 3. The particular symmetry breaking scenario that is realized in a *SO*(10) model depends exclusively on the scalar sector of the theory and the configuration of the vacuum, and it constitutes one of the major differences among *SO*(10) models in the literature.

**Figure 3**. Patterns of symmetry breaking from *SO*(10) to the SM group [21].

Regular *SO*(10) models may suffer from some of the same issues as regular *SU*(5), namely rapid proton decay can occur with a low unification scale. The main solution to this problem, as it was with *SU*(5), is the addition of supersymmetry. SUSY *SO*(10) [90, 91] theories are rather popular and given the large number of degrees of freedom they possess, such as symmetry breaking pattern, field content, etc., they can easily avoid many of the traditional issues in unified theories.

Alike to the *SU*(5) model, it is possible to construct alternative embeddings of the SM fermions into representations of the group via the addition of an Abelian group. In the *flipped* *SO*(10) × *U*(1) model [92] the SM fermion content is not fully embedded into the **16** representation of the group, but rather into the direct product **16**_{1}⊕**10**_{−2}⊕**1**_{4}. This model loses its unified nature in favor of more degrees of freedom for the Yukawa and symmetry breaking sectors of the theory, which are no longer constrained by the statements above [80].

### 2.5. *E*_{6}

The GUT models described so far have unification groups that are part of the infinite series *SU*(*n*) or *SO*(2*n*), which means that for each successful model with particular *n* there is an infinite number of alternatives with order larger than *n*. For instance, the *SO*(18) group has been studied as a candidate for gauge and family unification [93]. The exceptional algebras, however, are unique so they are more aesthetically appealing candidates as theories of unification [20, 94]. Among all exceptional algebras, only *E*_{6} is large enough to contain the SM as a subgroup and admits complex representations.

The fermionic content in the *E*_{6} theory is embedded in the fundamental **27**-dimensional representation of the group, which contains the SM fermions plus exotic fields. The particular allocations of SM fermions into the **27** representation depends on the subgroup of *E*_{6} that breaks into after SSB, which is typically either the trinification group, *SU*(3)_{c} × *SU*(3)_{L} × *SU*(3)_{R} [95–97] or *SO*(10) × *U*(1) [98, 99]. The decomposition of the fundamental **27** into these subgroups is

The field content in *E*_{6} models is quite vast. There are 78 gauge bosons, of which only 45 survive at lower scales if *SO*(10) is the breaking path, or even fewer in the case of *SU*(3) × *SU*(3) × *SU*(3), just 24. The minimal scalar content needs at least a scalar field in the combination $\overline{27}+\overline{351}+\overline{35{1}^{\prime}}$, which contains the SM Higgs, and a scalar responsible for SSB of *E*_{6}, which is dependent upon the pattern of symmetry breaking, e.g., a **78** for breaking to *SO*(10) × *U*(1).

One of the main motivations for *E*_{6} as a unification group comes from superstring theory, where it was shown to emerge as a four-dimensional compactification of the heterotic *E*_{8} × *E*_{8} superstring theory [24, 25]. In fact, the presence of compactified extra dimensions near the scale of *E*_{6} breaking can trigger symmetry breaking of the *E*_{6} group via the Hosotani mechanism [100] straight into *SU*(3) × *SU*(2) × *U*(1) × *U*(1) × *U*(1).

Most of the research in *E*_{6} theories has been typically considered only within the scope of supersymmetry, which ties in with their motivation as low-energy limits of superstring theory where spacetime supersymmetry emerges naturally after compactification. Being a simple Lie group, *E*_{6} benefits from the prediction of gauge coupling unification in supersymmetry, as did *SU*(5) and *SO*(10), which strengthens the motivation. In addition to the rich phenomenology of the MSSM, the *E*_{6} model adds quite a few phenomenological predictions on its own, from exotic fermion states to new heavy gauge bosons [101, 102].

## 3. Selected Topics in GUTs

Model building in unified theories involves more than the selection of the group theory properties as introduced in section 2. There are a few outstanding issues that need to be addressed in order to construct a realistic model. Gauge coupling unification is typically one of the most pressing issues, which can often be resolved by intermediate steps in the breaking chain or by the addition of supersymmetry to the theory. In this section we describe the interplay between SUSY and GUTs, as well as other topics such as the generation of light neutrino masses and the unification of the Yukawa couplings.

### 3.1. Supersymmetry and GUTs

Supersymmetry is a very appealing theory on its own right. It is one of the most aesthetically pleasing extensions of the Standard Model and it has an extremely rich phenomenology that can be readily tested at colliders and other experiments. SUSY GUTs [22] are a conglomerate of the numerous advantages of unified theories and the predicting power of supersymmetry. One of the most attractive features of SUSY is that it can stabilize the electroweak scale against quantum corrections, the so called *hierarchy problem* [103, 104] and provides a mechanism for dynamic electroweak symmetry breaking [105, 106].

In addition, if *R*-parity is conserved [107] the lightest supersymmetric particle (LSP) is stable. Therefore, SUSY automatically predicts the existence of a Dark Matter candidate and can easily produce scenarios with the correct relic abundance [108, 109].

As previously mentioned, one of the major motivations for SUSY GUTs is that the minimal MSSM model predicts gauge coupling unification at some high scale ~ 10^{16} [110]. As was seen in Figure 2, just taking the one-loop RGE flow of the gauge couplings, the unification at the GUT scale is fairly successful. The RGEs for the gauge couplings at one-loop have an analytic solution of the form [111]

where *i* = 1, 2, 3 labels the coupling of the *U*(1), *SU*(2), and *SU*(3) subgroups of the SM gauge group, and *b*_{i} are parameters that depend on the field content. For the MSSM these are *b*_{i} = (33/5, 1, −3). With a degenerate sparticle spectrum at 100 GeV, these one-loop RGEs unify at ${M}_{GUT}~2.5\times 1{0}^{16}$ GeV with α_{GUT} ~ 0.0388.

This picture, however, relies on a light and almost degenerate supersymmetric spectrum. For heavier or split spectra the situation changes drastically, often spoiling gauge unification altogether. A unification measure can be defined to assess how the unification of gauge couplings changes with the masses of the supersymmetric particles as

where μ_{ij} is the energy scale at which ${\alpha}_{i}^{-1}$ and ${\alpha}_{j}^{-1}$ unify. Figure 4 shows how the unification measure varies with respect to the SUSY scale for an MSSM model with degenerate SUSY masses (blue). One can distinctly see that for larger sparticle masses, the unification of gauge couplings significantly worsens, from a 70% unification for *M*_{SUSY} ~ 100 GeV to <30% at *M*_{SUSY} ~ 1 TeV. Consequently, in addition to solving the little hierarchy problem without too much fine tuning, a light sparticle spectrum is clearly preferred to achieve gauge coupling unification.

**Figure 4**. Evolution of the unification measure with the SUSY scale for the MSSM with almost degenerate masses (blue) and scenario with lighter electroweakinos (green), calculated at one-loop. Shaded regions include threshold corrections over the solid lines.

In spite of the appeal of SUSY GUTs, the combined effort of several collider experiments has not found any clear evidences of SUSY particles^{4}. Hence, minimal and light versions of the MSSM are in tension with experimental evidence and that makes achieving gauge coupling unification much harder. This tension relaxes slightly once the mass degeneracy condition is forgone. If the sparticle masses vary considerably across the spectrum, it is possible to evade experimental bounds for those masses more strongly constrained (e.g., gluinos, squarks, etc.) while keeping part of the spectrum light. Mass splittings across the supersymmetric spectrum can be taken into account by the contribution of threshold corrections at the SUSY scale, which are of the type [112, 113]

with ϕ the scalar fields in the MSSM (sfermions), *m*_{ϕ} their masses and *S*(ϕ) their Dynkin indices; and ψ the fermions in the MSSM (gauginos and Higgsinos), *m*_{ψ} and *S*(ψ) their masses and Dynkin indices. The shaded blue region in Figure 4 corresponds to MSSM models with slightly non-degenerate masses. Although these models exhibit the same trend as before, decreasing the unification measure as *M*_{SUSY} increases, some of these achieve a better gauge coupling unification than the degenerate case, with up to 60% unification for *M*_{SUSY} ~ 1 TeV.

Supersymmetric models with large splittings in the particle spectrum can modify this picture significantly. A special case, known as Split Supersymmetry [114, 115], has all the sfermions decoupled at the GUT scale and only gauginos and Higgsinos remain light, protected by chiral symmetry. This model is very well motivated within the context of unified theory, because the decoupled fields form full multiplets of *SU*(5), so the unification of the gauge couplings is not affected [114]. Hence, the behavior of the unification measure for these models is identical to the semi-degenerate MSSM case from above (blue line and shaded region in Figure 4), but has the advantage of allowing a lighter spectrum since some of the strongest experimental constraints are on the squark masses, which are decoupled from the spectrum. These Split-SUSY models, however, predict the existence of a light gluino, which is unfortunately strongly constrained by experiments. Alternative versions of this model with light electroweakinos (~ 100 GeV) and heavy gluinos (~ 5 TeV) have been studied [116], but these fail to provide successful gauge coupling unification for a light spectrum, slowly improving at larger scales, as can be noticed from the green solid line and shaded region (threshold corrections) in Figure 4. Split-SUSY and the light electroweakino model are just two extreme cases, the former requiring very light spectrum for successful unification and the latter a heavy spectrum. A number of models can be constructed with different spectra that have intermediate predictions for gauge coupling unification. In fact, with a precise analyses of threshold corrections, a number of supersymmetric models with large mass splittings have been shown to achieve exact unification, with a relatively light spectrum [117].

### 3.2. Neutrino Masses

The observed oscillations of neutrino flavors [10, 118] require any successful extension of the SM to incorporate non-zero masses for at least two neutrino species. Effectively, these masses are generated by the 5-dimensional Weinberg operator

where *c*_{W} is the Wilson coefficient, Λ denotes the operator's cut-off scale and *L* and *H* are the lepton and Higgs doublets, respectively. A typical UV-completion of this effective operator is some kind of seesaw mechanism [66, 119], which allows to satisfy elegantly the requirement of tiny neutrino mass size. Generally, these neutrino mass schemes assume a presence of new, heavy degrees of freedom, which are ideally motivated by other BSM physics. As has been shown, a number of different seesaw set-ups can be very naturally incorporated within the GUT framework. Provided that only a single type of new particle is added to the SM field content, there are three basic tree-level seesaw types [120].

#### 3.2.1. Seesaw Type I

In the original and simplest seesaw mechanism of type I the right-handed neutrino singlets must be added to the model [66, 121–123]. As the current experimental data require only two neutrinos to be massive, the minimal scenario must include two right-handed neutrino states. This extension then allows to write both Dirac and Majorana neutrino mass terms

where *y*^{ν} is the matrix of neutrino Yukawa couplings and *M*_{M} denotes the Majorana mass matrix. Hence, taking ${m}_{D}={y}^{\nu}v$, with *v* being the SM Higgs vev, the neutrino mass matrix can be written in the usual form

The block-diagonalisation of this matrix leads to the light mass of the oscillating neutrinos

as the Majorana mass parameter can be chosen to be arbitrarily large. Considering the neutrino Yukawa couplings of order one and the Majorana mass around 10^{14} GeV, the desired neutrino mass sizes of order *m*_{ν} ≈ 0.1 eV are obtained. The type-I seesaw mechanism can be implemented in the GUT framework. Particularly, it arises very naturally in *SO*(10) GUT, where the right-handed neutrino singlet can be accommodated together with all the other fermions in a single **16**_{F} spinor representation.

#### 3.2.2. Seesaw Type II

The second possibility to construct a seesaw mechanism is to assume a heavy scalar *SU*(2)_{L}-triplet Δ_{L},

which allows to write the following Lagrangian terms

Diagonalisation of the type-II seesaw mass matrix [67, 124, 125] then generates neutrino mass

and for *M*_{Δ}≫*v* the required suppression is obtained.

Also this seesaw can be responsible for neutrino mass generation in GUTs. For instance, in *SO*(10) unification the left-handed scalar triplet is contained by the **126** Higgs field, which is usually considered to be present in the scalar sector of the theory. It has been shown that type-II seesaw can be the dominant neutrino mass scheme within both SUSY [126] and non-SUSY [127] *SO*(10) GUTs.

#### 3.2.3. Seesaw Type III

The third option for a UV-completion of the Weinberg operator in Equation (14) is to introduce new fermionic *SU*(2)_{L} triplets ${\text{T}}_{F}^{I}$ [128]^{5}. Their interaction with the SM content is analogous to the type I seesaw, namely,

The neutrino mass matrix for type III seesaw then reads

and for ${M}^{{T}_{F}}\gg {y}^{{T}_{F}}v$ the smallness of neutrino masses is ensured.

The incorporation of the type-III seesaw mechanism into GUTs has been studied in literature [39, 120, 129]. When implemented within *SU*(5) models, type-III seesaw comes automatically in hand with the type-I seesaw, as both fields responsible for these mechanisms share the same adjoint representation **24**_{F}.

#### 3.2.4. Inverse Seesaw

At low energies the light neutrino masses can be generated at tree level via the so called inverse seesaw mechanism. This string theory motivated [25] scheme can be constructed when a non-minimal lepton content of a given theory is assumed. Namely, extra singlet leptons must be added to the model, which is in general allowed for any gauge theory [119]. The minimalistic extension of the SM particle content leading to inverse seesaw requires a pair of left-handed two-component lepton singlets *N*^{c} and *S* [98]. Taking three generations of these new singlet fields, one can write the 9 × 9 mass matrix of the neutral leptons in the basis $\left\{{\nu}_{L}^{\ell},{{N}^{I}}^{c},{S}^{A}\right\}$ (with *A* = *a, b, c*) as

where M and μ are the mass matrices corresponding to the *SU*(2)_{L} singlets, while *m*_{D} is the Dirac neutrino mass matrix as usual. As predicted by some string models, the Majorana mass entries corresponding to ν_{L} and *N* are zero. Thus, the only Majorana mass parameter is the matrix μ, which corresponds to the extra singlet *S*. This entry is then responsible for lepton number violation. If μ is set to be zero, the *B* − *L* symmetry is restored, the matrix M^{IS} degenerates and the three oscillating neutrinos become massless.

On the other hand, for non-vanishing μ such that μ≪*m*_{D}≪*M* the resulting mass matrix of the light neutrino eigenstates reads

The main difference from the standard seesaw scenarios is that in the present case neutrinos become light for μ → 0, not for large values of Majorana mass parameter. This is also the reason why one talks about “inverse” seesaw. As vanishing μ enhances the symmetry of the theory, the assumption of its small value can be considered to be natural [130, 131].

#### 3.2.5. Linear Seesaw

A particularly interesting realization of the inverse seesaw mechanism can be constructed within the *SO*(10) GUTs framework with broken D-parity [132]. The so called linear seesaw mechanism consists in extending the minimal fermionic content of the *SO*(10) model, contained by three copies of the **16**_{F} representation, by three gauge singlets *S*^{A}. The original version of this scheme was designed within the supersymmetric *SO*(10) framework; however, it can be constructed also in non-supersymmetric scenarios. The mass matrix for the neutral fermions in the basis $\left\{{\nu}_{L}^{\ell},{{N}^{I}}^{c},{S}^{A}\right\}$ has the following form

Here, m_{D} denotes the Dirac neutrino mass, M is the heaviest Dirac neutrino mass term mixing *N*-*S* and m_{L} stands for the small term mixing ν-*S*, which breaks the (*B* − *L*) symmetry. The light neutrino masses are then given by the expression

which depends linearly on m_{D} (and therefore also on corresponding Yukawa couplings). In the present scenario it is the large parameter *M* given by the unification scale what ensures the smallness of neutrino masses. Hence, the lightness of neutrinos is independent of the (*B* − *L*) symmetry breaking scale, which can consequently lie at low, experimentally accessible energies without spoiling the desired size of neutrino masses or the unification.

#### 3.2.6. Other Neutrino Mass Models

Despite the success of seesaw mechanisms, one can think of a number of alternative neutrino mass generation schemes. From the phenomenological point of view, these can be even more interesting, as they often predict (unlike the three usual seesaws) a low-energy origin of neutrino masses. The light neutrino masses are obtained using a small lepton-number-violating parameter (similarly as in the inverse seesaw), or they can be suppressed by loops and small Yukawa couplings. While the former option can be realized e.g., within supersymmetric models with R-parity breaking [133], the latter possibility refers to the models of neutrino mass generation via calculable radiative corrections (i.e., the Zee mechanism) [41, 134]. A two-loop mechanism generating neutrino masses within a minimal *SO*(10) GUT was identified by Witten [135] and the same scheme can be constructed also in the flipped *SU*(5) context [136–138].

### 3.3. Yukawa Coupling Unification and Fermion Masses

In fully unified theories, such as *SU*(5) and *SO*(10), the gauge couplings must unify at some high energy scale. This is typically achieved automatically in SUSY GUTs, as mentioned above, due to the RGE flow of the MSSM gauge couplings (c.f. Figure 2), but it can also be achieved through the addition of new scalar representations [139] or with a multi-step symmetry breaking pattern [91].

Along gauge coupling unification, *SU*(5) and *SO*(10) models also require the unification of the Yukawa couplings. The largest hierarchy on the fermion masses happens in the third generation where *m*_{t}/*m*_{b} ~ 40 and *m*_{b}/*m*_{τ} ~ 2.3, hence Yukawa unification in GUTs is always studied within the third generation only. In *SU*(5) the charged leptons live in the same representation as the down-type quarks, $\overline{5}$
, and as such it is expected that at the GUT scale *y*_{b} = *y*_{τ}, whereas in *SO*(10) all SM fermions are embedded into the same 16-dimensional representation, so the unification condition becomes *y*_{t} = *y*_{b} = *y*_{τ}.

Although a natural prediction of GUTs, Yukawa unification is not easily achieved in the vanilla *SU*(5) and *SO*(10) models [29–31]. As can be seen in the left-side plot of Figure 5, the Yukawa couplings in the SM are far from unification. In spite of this, a few successful attempts to solve the unification of *y*_{b} and *y*_{τ} in *SU*(5) inspired models have been performed, either by including large scalar *SU*(5) representations to the field content [31, 139, 140] or by adding Planck scale suppressed interaction of the Higgs field to the SM fermions [141].

**Figure 5**. One loop renormalisation group flow of the SM **(left)** and MSSM **(right)** Yukawa couplings, with *m*_{0} = 2 TeV, *m*_{1/2} = 3 TeV, *A*_{0} = 0 and tan β = 40 (solid), tan β = 30 (dashed) and tan β = 15 (dotted).

In SUSY GUTs, however, Yukawa unification can often be achieved in some regions of the full supersymmetric parameter space. As can be seen in the right-side picture in Figure 5, the Yukawa couplings in the MSSM tend to run toward convergence at high scales, and they can be seen to almost unify for large values of tan β [142–145]. This occurs because the third generation fermion masses depend on tan β in the following way [146]

which can realize the hierarchy *m*_{t}≫*m*_{b}, *m*_{τ} even in *SO*(10) or *E*_{6} models where one expects *y*_{t} ~ *y*_{b} ~ *y*_{τ}. These solutions with large tan β can spoil radiative EW symmetry breaking in unified models, since $B\mu \simeq \frac{{M}_{A}^{2}}{tan\beta}\simeq 0$ implies that ${m}_{{H}_{d}}^{2}-{m}_{{H}_{u}}^{2}>{m}_{Z}^{2}$ [147], contrary to the unified picture where ${m}_{{H}_{u}}^{2}={m}_{{H}_{d}}^{2}$. This issue can often be resolved either by *ad hoc* splittings of the Higgs soft masses at the GUT scale, or by considering the effect of *D*-terms in the boundary conditions at the GUT scale [148], which naturally imposes a splitting of ${m}_{{H}_{d}}^{2}-{m}_{{H}_{u}}^{2}=4{m}_{D}^{2}$.

In addition to satisfying *m*_{t} ≫ *m*_{b}, *m*_{τ} for unified Yukawa couplings, one can lift the hierarchy between *m*_{b} and *m*_{τ} with the inclusion of radiative corrections on the *b* mass. At one loop the *b* quark couples to *H*_{u} via a gluino or Higgsino loop [146], as can be seen in Figure 6, which adds a correction to *m*_{b} of the type [149]

Though similar corrections appear for *m*_{t} and *m*_{τ}, they are negligible compared to δ*m*_{b}. The correction on *m*_{t} is not proportional to tan β, which is required to be large to satisfy *t*−*b*−τ unification. Further, δ*m*_{τ} does not have a gluino loop and the Higgsino contribution is inversely proportional to ${m}_{\stackrel{~}{{\nu}_{t}}}$ which is typically much larger than ${m}_{\stackrel{~}{t}}$, and therefore the contribution is small. These radiative corrections on *m*_{b} are proportional to tan β and therefore can be significant, up to 50% [148], which can spoil the hierarchy *m*_{t} ≫ *m*_{b}. There are regions of the SUSY parameter space, however, where it is possible to reduce δ*m*_{b} while keeping tan β large [148–150], thereby successfully predicting *t*−*b*−τ unification, even factoring in LHC searches [151].

In SUSY *SU*(5) models the more straightforward boundary condition *y*_{b} = *y*_{τ} is imposed. It was found that, in addition to the large tan β scenarios from above, *b*−τ unification can also be achieved in a region of parameter space with low tan β [152, 153]. However, a sufficiently low tan β might struggle to lift sufficiently the tree level Higgs mass to the observed value, and hence there remains some tension between unified *b*−τ models of low tan β and the observed Higgs mass [154].

A number of other mechanisms have been proposed to satisfy the Yukawa unification conditions. Intermediate breaking steps, such as the Pati-Salam group, can modify the Yukawas RGEs in a favorable manner achieving quasi-unification [155, 156]. Non-canonical seesaw mechanisms in neutrino models require *b*−τ unification to match the observed neutrino mixings [157]. Or the inclusion of certain higher dimensional operators can successfully yield Yukawa unification [158].

Beyond the unification of the Yukawa couplings for each of the families of SM fermions, the mass hierarchies among the different families remains an open question. Although GUTs by themselves do not make predictions on the nature of this hierarchy, they often include a fair amount of parameters and mixing matrices that are unconstrained and can fit the fermion masses. Additionally GUTs are often extended with family symmetries, continuous or discrete, which can, with a smaller set of parameters, accurately predict the fermion mass hierarchies, as well as their mixings and CP phases encoded in the CKM and UPMNS matrices [159–163]. We will not discuss family symmetries any further since they fall beyond the scope of this work.

## 4. Modern Day GUTs

Since their first appearance in the late 70s, a large number of GUT models have been proposed. These vary according to the symmetry group employed, the symmetry breaking mechanism and field content among others. Some of them were driven by the experimental results of the time and other by new theoretical insights. In this section we attempt to outline a small, non-exhaustive, subset of GUT models, aiming to explore those with strong phenomenological consequences and some that have been in the spotlight in recent years. We thus focus on left-right symmetric models, SUSY *SO*(10), trinification models and E_{6}SSM.

### 4.1. Left-Right Symmetric Models

One of the minimal extensions of the SM is the earlier mentioned left-right symmetric model [16, 17, 62, 63, 67]. Despite not being real GUTs, LR models can very conveniently play the role of an intermediate symmetry restored between the electroweak scale and the GUT scale [69, 164, 165]. The LR framework has attracted a lot of attention particularly in connection with the LHC [166–178], as it typically predicts new physics at energies that can be probed by the collider searches.

The fermionic particle content of LR models is given by a straightforward LR symmetric extension of the SM content, i.e., the right-handed doublets are introduced

As a result, right-handed neutrinos are naturally included making the left-handed neutrinos acquire mass in the LR models, which is a highly desirable feature of a BSM model. The presence of the right-handed neutrino partners is also essential for cancellation of the *B* − *L* gauge anomaly.

The Higgs sector of LR symmetric theories can vary. The minimal scenarios mostly include a scalar bi-doublet^{6}

containing the SM Higgs, which subsequently gives masses to quarks and leptons. The corresponding vev reads

where $v\equiv \sqrt{{v}_{1}^{2}+{v}_{2}^{2}}$ and it mixes the left-handed and right-handed gauge bosons as described below.

Besides the bi-doublet, typically a pair of scalar triplets

or doublets

must be added to the Higgs sector in order to break the LR gauge group to the SM. In fact, the right-handed scalar is enough to do so, but inclusion of the left-handed triplet (or doublet) preserves the LR symmetry (so called “manifest LR symmetry”), i.e., the *SU*(2)_{L} and *SU*(2)_{R} gauge couplings are equal: *g*_{L} = *g*_{R}.

If no additional fermions besides the SM fermionic content are considered, at least two bi-doublets must be present in the scalar sector to account for the correct SM flavor physics [164]. In a model with a single bi-doublet the Yukawa Lagrangian implies that the up-quark mass matrix is proportional to the down-quark mass matrix (independently of the vev structure); thus, the CKM matrix becomes trivial *V*_{CKM} = 1.

Consequently, the LR symmetry breaking takes place in two steps. At first, the neutral component of right-handed scalar triplet (or doublet) gets the vev *v*_{R} and breaks the LR gauge group to the SM gauge group. Subsequently, the bi-doublet acquiring its vev breaks the SM gauge group to *SU*(3)_{C} ⊗ *U*(1)_{Q}. Based on the observations it can be assumed that *v*_{R}≫*v*_{Φ1}, *v*_{Φ2}.

Depending on the scalar content of a particular LR model, different ways of light neutrino mass generation can be employed. Having right-handed neutrino singlets means that type-I seesaw is always the option. In general, the neutrino mass matrix can take the form

where m_{D} denotes the Dirac mass matrix, while M_{M, L} and M_{M, R} are the Majorana mass matrices corresponding to the left-handed and right-handed neutrinos, respectively.

The Yukawa couplings in LR models include the scalar bi-doublet,

where $\stackrel{~}{\Phi}={\sigma}^{2}{\Phi}^{*}{\sigma}^{2}$. Then the Dirac neutrino mass matrix and the mass matrix of charged leptons are in the broken phase given by

In case that the right-handed scalar triplet Δ_{R} is responsible for the LR symmetry breaking, one can write also the Yukawa couplings for the right-handed lepton doublet in the form

where Δ_{R} = **Δ**_{R}·**τ**. After the triplet acquires its vev

the LR symmetry is broken and the right-handed neutrino receives Majorana mass ${\text{M}}_{M,R}={\text{y}}^{{\Delta}_{R}}{v}_{R}\gg v$, which allows for type-I seesaw mechanism.

When the Higgs sector contains also the left-handed scalar triplet Δ_{L} with vev

it generates the left-handed Majorana mass matrix ${\text{M}}_{M,L}={\text{y}}^{{\Delta}_{L}}{v}_{L}$ switching on type-II seesaw mechanism.

In principle, the type-I and type-II seesaws can be combined giving the “full” seesaw matrix (35). The resulting light neutrino mass matrix reads

Specifically, if *v*_{Φ2} = 0 is assumed for simplicity, then the formula (42) can be rewritten in terms of the parameters of the LR models as

Hence, if the hierarchy *v*_{R}≫*v*_{Φ1}≫*v*_{L} is satisfied, the neutrino masses become small.

In models with the LR symmetry breaking driven by the right-handed doublet χ_{R} instead of the triplet Δ_{R} the light neutrino masses can be obtained employing the inverse [98, 179–181] and/or linear [175, 182, 183] seesaw mechanisms, provided that a singlet fermion {**1, 1, 1**, 0} is added to the model particle content. Alternatively, it is also possible to construct the type-III seesaw mechanism, if a left-handed or right-handed fermionic triplet is present within the LR model [184, 185]. Lastly, neutrino mass generation in LR models via the Zee mechanism can be achieved with the addition of a charged scalar boson {**1, 1, 1**, 2} [186].

### 4.2. SUSY *SO*(10) Models

Supersymmetric *SO*(10) models are rather appealing GUTs, for they combine together the advantages of SUSY, Pati-Salam and *SU*(5) models. As briefly outlined in section 2.4, *SO*(10) models unify all fermions of a generation in the SM into a single representation, of dimension **16**. This decomposes into the maximal subgroups as

As mentioned before, at the renormalizable level the Yukawa sector of *SO*(10) includes the Higgs representations **10, 120**, and $\overline{126}$, which are promoted to superfields in SUSY *SO*(10). Hence, the superpotential of the Yukawa sector is

where Y_{i} are matrices of Yukawa couplings in family space. One of the most remarkable features of SUSY *SO*(10) is that, starting from a Yukawa unified scenario, Y_{10} and Y_{126} are sufficient to reproduce the full mass spectrum of SM fermions, along with the measured values of mixings in the quark and neutrino sector [187–197]. The mass matrices of SM fermions M_{i} can be written as [86]

where M_{D}, M_{L}, and M_{R} are the Dirac and Majorana masses in types I and II seesaw (c.f. section 3.2), and the *v*'s are the various vacuum expectation values of **10**, $\overline{126}$ and the left and right-handed *SU*(2) triplets.

The minimal SUSY *SO*(10) model therefore contains the Higgs superfields **10** and $\overline{126}$, responsible for fermion masses, and a pair of representations **126** and **210** which trigger the symmetry breaking of *SO*(10) [198–200]. Although quite appealing due to its minimal set of model parameters, this model does not achieve the right level of gauge coupling unification and suffers from rapid proton decay [87, 201].

Many solutions have been implemented to resurrect minimal *SO*(10) models. The spectrum of soft masses in the low energy MSSM strongly affects the outcome of gauge coupling unification, as was seen in section 3.1, hence modifications on the hierarchy of soft masses can help toward solving the issues with SUSY *SO*(10) models [202–206]. Additionally, extended scalar sectors, either containing a **120** [207, 208] or a **54** [209] representation, can increase the unification scale through strong threshold effects, thereby alleviating the constraint of nucleon decay. Recently it has been shown that a combination of new Higgs representations with a modified spectrum of soft masses can accommodate gauge coupling unification and nucleon decay constraints, while still being able to predict a suitable spectrum of fermion masses [210].

As with many GUT models, SUSY *SO*(10) makes predictions that can be tested in a number of different fronts. Collider searches at the LHC [211, 212] as well as dark matter searches [213] can discover the predicted light, TeV-scale, states. Precision tests such as nucleon decays [214, 215], lepton flavor violation [216] and flavor observables [217] can probe the validity of the models at high scales. For more details on probing SUSY *SO*(10) and GUTs in general see sections 5–7.

### 4.3. Trinification

As a maximal subgroup of *E*_{6}, the trinification gauge group *SU*(3)_{c} × *SU*(3)_{L} × *SU*(3)_{R} is an alternative approach to SUSY *SO*(10) on the road to *E*_{6} unification. The matter content of trinification models per generation typically looks like [218]

where *h*_{ij} are the components of two Higgs doublets, ϕ a SM singlet field and *D* and *D*^{c} color-triplets. An additional *Z*_{3} symmetry is often considered to make the gauge couplings unify at the GUT scale, *g*_{c} = *g*_{L} = *g*_{R}.

Given the presence of exotic fields in the matter multiplets, trinification models struggle to trigger spontaneous symmetry breaking without making the matter content impossibly heavy. Additional Higgs multiplets [219–222] and/or non-renormalizable operators [223–226] are usually introduced to alleviate this issue. Unfortunately these models tend to produce tension with current limits on proton decay [227] and collider searches [228, 229].

The fundamental challenge behind the issues of trinification is the complicated vacuum structures and the large number of parameters they depend on. Solutions to this problem involve the addition of family symmetries which reduce the number of parameters and facilitate the study of the symmetry breaking vacuum [230], further simplified by embedding the theory into larger dimensional groups such as *E*_{8} [231, 232].

### 4.4. E_{6}SSM

The Exceptional Supersymmetric Standard Model (ESSM or E_{6}SSM) [233, 234] is an extension of the MSSM motivated as a low energy effective theory from a *E*_{6} unified GUT model at high scales. At low scales it has the gauge group *SU*(3)_{c} × *SU*(2)_{L} × *U*(1)_{Y} × *U*(1)_{N}, where the additional *U*(1)_{N} factor is leftover from the symmetry breaking of *E*_{6}. All the superfields in the theory are contained in three copies of the **27** representation of *E*_{6}, which decompose under the *SU*(5) × *U*(1)_{N} subgroup as [235]

where ${10}_{1}^{i}$ and ${\overline{5}}_{2}^{i}$ are the matter multiplets for all three generations, with the standard embeddings of matter fields in *SU*(5) (c.f. section 2.1). The superfields ${\overline{5}}_{-3}^{i}$ and ${5}_{-2}^{i}$ contain the two Higgs doublets of the MSSM, *H*_{u} and *H*_{d}, plus two copies of pairs of exotic doublets, ${H}_{u}^{1,2}$ and ${H}_{d}^{1,2}$ and three copies of exotic triplets *D*_{i} and ${\overline{D}}_{i}$. Lastly, the singlets ${1}_{5}^{i}$ and ${1}_{0}^{i}$ correspond to exotic singlet fields *S*_{i}, responsible for *U*(1)_{N} breaking at low scales, and right-handed neutrino fields *N*_{i}, respectively.

Anomaly cancellation of the *U*(1)_{N} factor in the E_{6}SSM model is guaranteed so long as the only decoupled state is the singlet neutrino field. *N*_{i} can be as large as necessary to provide light neutrino masses through type I seesaw mechanism and generate the baryon assymmetry of the Universe via leptogenesis [236, 237]. The remaning fields of the **27**^{i} multiplets charged under *U*(1)_{N} remain at energies below the breaking of *U*(1)_{N} and hence anomalies cancel. Light colored states have dangerous consequences, however, for they can mediate baryon and lepton number violating interactions leading to rapid proton decay. In order to avoid that, the original E_{6}SSM model postulates the existence of an approximate *Z*_{2} symmetry that forbids those interactions. An exact *Z*_{2} symmetry can also be considered [238, 239], but in such a case additional exotic states must be introduced to ensure that the exotic quarks are not stable.

Gauge coupling unification in the E_{6}SSM model requires the addition of incomplete multiplets of *E*_{6} at low scales, since full multiplets do not modify the unification properties of the RGE flow. A pair of fields *H* and $\overline{H}$ in conjugate representations are added, to ensure no anomalies are reintroduced. Alternatively, a Pati-Salam intermediate step has been postulated that achieves gauge coupling unification without the need of additional superfields. This “minimal” E_{6}SSM model, however, predicts unification at the Planck scale so quantum gravity corrections may play a role and affect the outcome of unification [240].

The E_{6}SSM has a rather rich phenomenology since most of the predicted exotic states live at low energies. The constrained E_{6}SSM (cE_{6}SSM) is a version of the full E_{6}SSM that exploits the properties of unification of *E*_{6} and assumes universal scalar and gaugino soft masses at the GUT scale [241, 242]. Predictions of the cE_{6}SSM include contributions to the Higgs mass and rare decays [243, 244] together with light exotic states, such as the ${Z}_{N}^{\prime}$ associated with the *U*(1)_{N} broken symmetry, and the color triplet fermions *D* and $\overline{D}$, all of which can be probed at the LHC [245]. Lastly, as in the MSSM the lightest supersymmetric particle is stable, so it is a valid candidate for dark matter. In contrast with the regular neutralino dark matter in the MSSM, the dark matter candidate in the E_{6}SSM is a mixture of binos, winos and higgsinos, as well as the inert singlinos and higgsinos in ${H}_{u,d}^{i}$ and *S*_{i} [246–250].

## 5. Cosmology and the Early Universe

### 5.1. Inflation and GUTs

Cosmic inflation plays an important role in theories of Grand Unification, as it is needed to dilute relics such as magnetic monopoles, which are produced ubiquitously in GUT models^{7}. The requirement to dilute these relics therefore determines the scale of inflation in specific models [251]. Moreover, since generic inflation models are associated with a scale ${\Lambda}_{\text{inflation}}~1{0}^{16}$ GeV, it becomes attractive to associate the inflaton with a GUT scalar.

To agree with observations, inflation models need to predict a large number of observable e-folds ($N={\int}_{{t}_{0}}^{{t}_{e}}Hdt\gtrsim 60$), as well as small spectrum density of fluctuations δ*ρ*/*ρ* ~ 10^{−5}. For an effectively single field model, this can be illustrated by the tension between the Lyth bound (a measure of the field excursion necessary to solve the problems inflation was invented to solve) given in terms of the number of e-foldings *N* [252],

and the amplitude of the Cosmic Microwave Background (CMB) anisotropies, which implies [253, 254],

It is clear, then, that inflation requires a very flat scalar potential. Hence, it is attractive to consider inflation models in which the potential is dynamically generated [255–258].

A successful example of such a model was realized as early as 1983 [259]. This paper considered a potential of a Coleman-Weinberg form [255],

Such a potential can only be made compatible with CMB constraints if *A* is very small (presently, *A* ≲ 10^{−14} [260]). Therefore, the potential in Equation (51) cannot be due to loops of *SU*(5) gauge bosons. Instead, [259] considered a model in which the inflaton transformed as a singlet of *SU*(5), couples weakly to the adjoint and fundamental Higgs fields, and therefore obtains a vacuum expectation value when *SU*(5) breaks into the SM. The original CW-GUT inflation model [259] predicts primordial gravitational waves with tensor-to-scalar index 0.02 ≤ *r* ≤ 0.1 [261]. Although this is currently not in tension with the CMB-constraints [253, 254], several modifications have been proposed which predict smaller *r* [260, 262–264].

An alternative class of GUT inflation models are based on no-scale supergravity. It was realized in 2013 [265] that particular realizations of no-scale supergravity (SUGRA) models of inflation can be equivalent to the Starobinsky model [266], in which inflation is realized from a non-minimal Einstein-Hilbert action $S=\frac{1}{2}\int dx\sqrt{-g}(R+{R}^{2}/6{M}^{2})$. The correspondence can be seen by a conformal transformation, such that the model is equivalent to canonical gravity plus a scalar field [265, 267]. The scalar potential then becomes

Starobinsky-like models are attractive candidates for inflation models, as they make viable predictions for inflationary observables without the need to introduce a large set of finely tuned parameters. Examples of no-scale SUGRA inflation models include sneutrino inflation, which can be consistently implemented in flipped *SU*(5) SUSY GUTs [268–272]. Other no-scale GUT inflation models identify the inflaton with the Higgs boson, and circumvent the problems associated with conventional Higgs inflation [273, 274].

A phenomenological approach was taken by Hertzberg and Wilczek [275]. Here it was assumed that inflation is driven by the vacuum energy associated with unification. It was shown that several examples of large-field (Δϕ ~ *M*_{p}) models of inflation give predictions consistent with the CMB-constraints [253, 254].

GUT-inflation has also been studied in combination with other cosmological scenarios, most importantly with non-thermal leptogenesis and the seesaw mechanism for neutrino masses [276–278]^{8}. Models of sneutrino inflation are well suited for this purpose [268–272].

### 5.2. Cosmological Constraints on Cosmic Strings

Cosmic strings are generic cosmological predictions of many GUT theories [279–281]. Field theories with broken gauge symmetries may have a vacuum state that is not unique, such that different points in physical space may have distinct (but degenerate) vacuum configurations. By continuity of the field, the interpolating field values must be taken on in between these points, which gives rise to an energetic object called a topological defect, or (in the one-dimensional case) a string.

The simplest description of cosmic strings after their formation approximates the fundamental *Nambu-Goto* strings. Nambu-Goto strings are characterized by the dimensionless string tension *Gμ*, where μ is the mass per unit length and *G* is Newton's constant. Strings produced at the GUT scale typically have a mass per unit length of the order of μ ~ 10^{21}kgm^{−1} and a thickness of 10^{−24}m, such that the tension is *Gμ* ~ 10^{−6} [282]. For comparison, strings formed around the electroweak scale are expected to have much smaller tensions, *Gμ* ~ 10^{−34}. As the expansion of the Universe stretches strings, while the string tension stays constant and in the absence of a decay mechanism, ρ_{strings} would grow with the scale factor and eventually dominate the energy density of the Universe. Cosmic strings can decay into gravitational or scalar radiation, however. In the presence of such a decay channel an attractor scaling solution is reached, in which the strings maintain a constant fraction of the energy budget.

Cosmic strings could be detected through gravitational lensing and anisotropies in the CMB [283, 284]. Cosmic strings imprint on the CMB as line-like discontinuities, caused by a boost of photons toward the observer as a string moves across the line of sight [285, 286]. Planck data constrains the Nambu-Goto string tension *Gμ* < 1.8 × 10^{−7} [287], the non-local string tension *Gμ* < 10.6 × 10^{−7} [287] and the Abelian Higgs string model *Gμ* < 2.0 × 10^{−7} [288].

If the strings decay gravitationally, such radiation can be observed in Gravitational Wave (GW) experiments [289–293]. Strong gravitational radiation bursts may be produced by cusps [294–298]: the LIGO/VIRGO collaboration reported an experimental upper limit of *Gμ* < 10^{−8} in some regions of the parameter space, in which gravitational backreaction determines the size of the loops [299]. Pulsar Timing Arrays (PTAs) potentially give more stringent bounds, as they can already probe the stochastic GW background; depending on the model, $G\mu <{O}(1{0}^{-12}-1{0}^{-11})$ [300–302]. However, the relative importance of the gravitational decay channel has been the source of some disagreement in the literature. Simulations of Nambu-Goto strings [303–308] and full field-theoretic simulations of the Abelian Higgs model [309–311] differ in the typical scale on which the strings form loops. Loops of the order of the string width *r*_{s} can radiate heavy particles (as the natural mass of coupled particles is $m~{r}_{s}^{-1}$); loops of the typical inter-string spacing ξ are expected to decay gravitationally [312]. Recent field-theory simulations [311] suggest a mechanism to transport energy from large to small scales, which is not currently understood. Numerical results [313–315] also indicate that the simplest analytical models are due an update, when non-Abelian gauge groups are considered. Furthermore, different groups simulating Nambu-Goto strings differ in the distribution of the loop size. Simulations in which the gravitational radiation back-reacts on the string have smoother features, which hinders the formation of smaller loops [303–305]. In these simulations, the PTA constraints on the stochastic background and LIGO/VIRGO constraints on burst are stronger [301].

The shape of the fractional energy density Ω_{GW} ≡ *f*/*ρ*_{c}(*dρ*_{GW}/*df*) power spectrum from cosmic strings is expected to be nearly scale-invariant, with an amplitude defined by the characteristic string tension *Gμ*, such that if it has a large enough amplitude, the signal would be seen in frequency windows of different experiments [316]. This distinguishes the power spectrum from other sources. In particular, an observation of GW at pulsar timing arrays, if coming from cosmic strings, will draw attention to interferometer searches for this source.

### 5.3. Gravitational Waves From Phase Transitions

Grand unification models can accommodate a rich scalar sector, which can result in a complicated cosmological history involving exotic phase transitions. Some GUT inspired possibilities are: a color breaking phase transition where color is broken and restored when leptoquarks acquire a vacuum expectation value in an intermediate transition, which can catalyse baryogenesis [317, 318], *B*−*L* and *L* violating phase transitions [319–322], hidden sector phase transitions [323–326], and a Pati-Salam transition [327]. If any such phase transition occurs through bubble nucleation, an observable relic gravitational wave spectrum can be seen today, for a review see Mazumdar and White [328], Caprini and Figueroa [329], and Weir [330]. Furthermore, GUT models often require the existence of extra singlets. For example, the E_{6}SSM model studied in Athron et al. [241] had 3 generations of singlets and such singlets can, in principle, catalyse the electroweak phase transition to be strongly first order as well [331–333].

The gravitational wave spectrum generated from a cosmic phase transition has three contributions: a contribution from the collision of scalar shells, and potentially long-lasting contributions from sound waves and turbulence in the plasma. The total spectrum can thus be written as,

Although much uncertainty remains about the precise form of these spectra, all three are controlled by four thermal parameters, which can be computed by first principles [328, 334]: the latent heat released during the transition (conventionally normalized by the radiation energy density), denoted α, the nucleation rate (conventionally normalized to the Hubble parameter) *β*/*H*, the temperature at which the transition occurs *T*_{*} and the velocity of the bubble wall *v*_{w}. The collision term is expected to be sub-dominant for transitions associated with a broken gauge group [335]. The sound wave contribution to the power spectrum is [334] and [336]

where ${\overline{U}}_{f}^{2}~(3/4){\kappa}_{f}{\alpha}_{T}$ is the rms fluid velocity, Γ ~ 4/3 is the adiabatic index, κ_{f} is the efficiency of converting the latent heat into gravitational waves and *g*_{*} the number of relativistic degrees of freedom. The frequency dependence is captured by the spectral state

The other notable, albeit sub-dominant, contribution is the contribution from magneto-hydrodynamic turbulence in the plasma. The power spectrum from this contribution is given by [337],

where ϵ is the fraction of the energy that contributes to turbulence, typically taken to be in the range (0.05, 0.1) [334]. In this case the spectral form is a function of two time scales,

where *h*^{*} is the Hubble rate at the transition temperature^{9}.

For a single scalar field transition, without a tree-level barrier between the true and the false vacuum, β/*H* tends to be *O*(10^{3}) or greater [232]. The transition temperature is the same order of magnitude as the mass of the scalar. Therefore, transitions with scalar masses *O*(10^{5}) GeV can be probed by ground-based interferometers such as the Einstein Telescope [338], Kagra [339] and cosmic explorer [340], whereas space-based LISA will probe transitions at the electroweak scale [334]. The former can be more directly related to studies of GUTs-we show benchmark examples for a Pati-Salam phase transition are shown in Figure 7. The visibility of the spectrum tends to grow with the ratio *v*/*m*, the gauge coupling constant *g*, the rank of the (sub) group being broken and the number of other particles acquiring a mass during the transition [323]. Furthermore, it was found in Croon et al. [323] that some non-trivial model discrimination is possible if one observes a primordial power spectra due to the increase in visibility as well as moderate correlations between thermal parameters, shown in Figure 8 for *SU*(*N*)/*SU*(*N* − 1) cosets.

**Figure 7**. Plane of thermal parameters with contours of constant sound wave (in black) and turbulence (in blue) peak amplitude. The thicker line shows the peak sensitivity of the Einstein telescope [338]. Points denote thermal parameters produced by a Pati-Salam phase transition with a Pati-Salam scale of ${M}_{\text{PS}}=1{0}^{5}$ GeV. The peak frequency for the sound wave spectrum is indicated by the color scaling. Figure taken from Croon et al. [327].

**Figure 8**. Thermal parameters from a renormalizable potential for a scalar field in the fundamental representation of *SU*(*N*_{G}) which gets broken to *SU*(*N*_{G} − 1) **(left). Right** panel is the same aside from the inclusion of a non-renormalizable sextet term in the potential and the gauge coupling constant is fixed to unity. In the above *N*_{G} denotes the order of the group, *N*_{f} is the number of fermions in the fundamental representation that are strongly coupled to the scalar field, in direct analogy with the SM save that the Yukawa couplings are set to unity. Note that ξ is the ratio of latent heat to radiation energy density more commonly denoted α. Contours of constant peak amplitude for the sound wave spectrum are shown with the darker line corresponding to LISA sensitivity range for a power law spectrum that has been integrated over frequency [341]. Note in the above *v*_{w} = 0.5 for the left plot and *v*_{w} = 1 for the right plot using the efficiency terms in Espinosa et al. [342]. Figure taken from Croon et al. [323].

If multiple scalar fields are involved in a transition the barrier between the true and false vacuum can persist at zero temperature due to triscalar or non renormalizable operators [333, 343]. In such a case significantly more supercooling is possible and the transition temperature is no longer confined to be the same order of magnitude as the scalar mass. This implies that β/*H* can be quite small and the latent heat can be large, increasing the visibility of the gravitational wave and reducing the peak frequency. A caveat to this is that recent work found that phase transitions that involve a large amount of supercooling may fail to complete due to the onset of inflation [344]. Regardless, the thermal parameter space in the case of multifield phase transitions is broader, which minimizes model discrimination somewhat, though not completely [323].

### 5.4. Baryo-/Leptogenesis

The existence of a baryon asymmetry in the Universe (BAU) is one of the central problems of modern cosmology [345, 346]. At the same time, the concordance between different measurements of the primordial baryon asymmetry is a triumph of modern cosmology with BBN and CMB limits giving [347, 348]

respectively. Any explanation for the baryon asymmetry must satisfy the three Sakharov conditions [349]

• Baryon number *B* violation

• C and CP violation

• A departure from thermal equilibrium.

Early attempts at generating the BAU focused on *B* violating decays (for a review see [350]). Such decays typically violate *B*+*L* while conserving *B* − *L* (for an exception see [351]). For example, *SU*(5) GUTs are invariant under changes to a global phase conjugate to *B* − *L* number, whereas *SO*(10) has a local *U*(1)_{B − L} symmetry. However, any primordial *B*+*L* asymmetry is washed out by *B*+*L* violating electroweak sphalerons. Therefore, only a primordial *B* − *L* asymmetry will be preserved unless sphalerons are quenched.

Leptogenesis allows for a *B* − *L* violating operator, $m{\overline{\nu}}_{R}^{c}{\nu}_{R}$, that is also responsible for a light neutrino mass via type-I seesaw mechanism (see section 3.2). CP violating decays of such sterile neutrinos ensure a net *B* − *L* asymmetry which electroweak sphalerons convert to a baryon asymmetry. Electroweak baryogenesis by contrast uses the sphalerons themselves to generate a net *B*+*L* asymmetry which cannot be washed out before the sphalerons are quenched [345, 346]. More specifically, if the electroweak phase transition is strongly first order, bubbles of electroweak broken phase populate a medium of symmetric phase with sphalerons quenched only inside such bubbles. Particles can have CP violating interactions with the bubble wall which biases the sphalerons to produce a net *B*+*L* asymmetry. Some of this asymmetry is swept up in the expanding bubble wall where it is preserved.

GUTs are only relevant to electroweak baryogenesis if the GUT model motivates some light BSM states. Recent work on electroweak baryogenesis in the *E*_{6}-SSM utilized three generations of singlet superfields to motivate a CPV source involving singlino-Higgsino interactions with the bubble wall [331].

A feature of leptogenesis during GUTs is the possibility of new CP violating decay channels due to the presence of leptoquarks [352, 353]. This allows a lower minimum mass for the lightest sterile neutrino than in the minimal scenario [236]^{10}. Much of the recent focus on baryogenesis within GUTs involves leptogenesis with some intriguing concordance achieved in the case of *SU*(5) [355] and *SO*(10) GUTs [271, 278, 356–360]. A generic feature of *SO*(10) GUTs is normal ordering of neutrino masses and a negative Dirac phase [361], both of which are favored by current observational limits [360]. Furthermore, many GUTs, including *SO*(10), predict a Dirac neutrino mass matrix that is not too different from the up quark mass matrix and *SO*(10) leptogenesis also achieves agreement in the currently observed atmospheric mixing angle in the first octant. Realistic models with two right-handed neutrinos can emerge in Δ(27) × *SO*(10) models and *A*4 × *SU*(5) supersymmetric models [362, 363]. The third right-handed neutrino can either decouple because it is very heavy or because its Yukawa coupling is very small. The latter case implies a stable particle that can play the role of dark matter [364].

## 6. Direct Collider Searches for GUTs

### 6.1. Searches for Supersymmetry

As we have seen in section 3.1 supersymmetry plays a rather important role on many unified theories and motivates the unification of gauge couplings at large scales. However, both its solution to the hierarchy problem and gauge coupling unification often rely on a light sparticle spectrum, around or below the TeV scale. Thus, searches for supersymmetric particles has been part of the research programme in collider physics for the last few decades, from searches at LEP and the Tevatron, to the recent results of the LHC, and it is still part of the proposed physics programme for future colliders, e.g., CLIC, ILC, or VLHC [365–368].

In R-parity conserving SUSY the lightest supersymmetric particle (LSP) is stable. This has strong consequences for SUSY searches, for the LSP will escape the collider in the form of missing transverse energy (MET)^{11} [369]. In addition, R-parity requires that sparticles are pair-produced in colliders, hence the different searches for supersymmetry are classified according to the particle that is produced in pairs. The production cross sections of the different sparticle species are rather different and often determine the exclusion or detection power of a particular channel. For instance, the strongest exclusion limits at the LHC across the sparticle spectrum are on first and second generation squarks and gluinos which, as can be seen in Figure 9, have the largest production cross sections [370, 371].

**Figure 9**. Pair production cross sections as a function of the mass of sparticle $\stackrel{~}{\varphi}$ at center of mass energy $\sqrt{s}=13$ TeV.

Squarks and gluinos are produced in pairs at the LHC in the combinations $\stackrel{~}{g}\stackrel{~}{g}$, $\stackrel{~}{q}\stackrel{~}{q}$ and $\stackrel{~}{q}\stackrel{~}{g}$ and their main decay channels are $\stackrel{~}{q}\to q{\stackrel{~}{\chi}}_{1}^{0}$ and $\stackrel{~}{g}\to q\overline{q}{\stackrel{~}{\chi}}_{1}^{0}$, with the neutralino LSP escaping the collider. Hence the typical signature for these processes has multiple jets and large missing energy. The decay topologies for these signatures are depicted in Figure 10. Other decay modes for squark and gluinos involve the production of charginos or heavier neutralinos, $\stackrel{~}{q}\to q{\stackrel{~}{\chi}}_{2}^{0}$, $\stackrel{~}{q}\to {q}^{\prime}{\stackrel{~}{\chi}}_{1}^{\pm}$, $\stackrel{~}{g}\to q\overline{q}{\stackrel{~}{\chi}}_{2}^{0}$ and $\stackrel{~}{g}\to q\overline{q}{\stackrel{~}{\chi}}_{1}^{\pm}$, which then decay to *W* and *Z* bosons and ${\stackrel{~}{\chi}}_{1}^{0}$. The final state signatures depend on the decay modes of the gauge bosons, and can have (0–4) leptons, jets and MET. ATLAS and CMS have reported results from the last run of the LHC at 13 TeV and 36 fb^{−1} for searches with jets and MET final states [372–374], with one lepton, jets and MET [375–377], same and opposite-sign dilepton pairs, jets and MET [378, 379], two or three leptons, jets and MET [380, 381], 3rd generation squarks (with and without Higgs reconstruction) and MET [382–384] and hadronic τ decays, jets and MET [385], among others. These searches set a lower limit for a range of simplified models on the mass of the gluino of ${m}_{\stackrel{~}{g}}\gtrsim 2.1$ TeV and the mass of the first and second generation squarks of ${m}_{\stackrel{~}{q}}\gtrsim 1.5$ TeV.

**Figure 10**. Most simple topologies for the production of squarks and gluinos at the LHC with decays to jets + MET.

The next strongest production cross section is that of stop and sbottom pairs. The main decay channel for stops and sbottoms is $\stackrel{~}{t}\to t{\stackrel{~}{\chi}}_{1}^{0}$ and $\stackrel{~}{b}\to b{\stackrel{~}{\chi}}_{1}^{0}$, respectively. This topology is similar to the decay of first and second generation squarks, with the added complexity that neither *t* or *b* produce a clean jet, but rather have many decay channels that can result in numerous jets, leptons and, of course, MET. Secondary decay channels for $\stackrel{~}{t}$ and $\stackrel{~}{b}$ involve decays to a chargino, $\stackrel{~}{t}\to b{\stackrel{~}{\chi}}_{1}^{\pm}$ and $\stackrel{~}{b}\to t{\stackrel{~}{\chi}}_{1}^{\pm}$, with subsequent decays involving *W* bosons, or decays into a heavy neutralino, $\stackrel{~}{t}\to t{\stackrel{~}{\chi}}_{2}^{0}$ and $\stackrel{~}{b}\to b{\stackrel{~}{\chi}}_{2}^{0}$, which in turn decays into a *Z* or a Higgs boson and ${\stackrel{~}{\chi}}_{1}^{0}$. The latest searches of the LHC experiments for pair-produced stops and sbottoms target final states with jets and MET [374, 386, 387], b-jets and MET [388, 389], one lepton, jets and MET [390, 391], two and three leptons, jets and MET [379–381, 392, 393] and final states with a *h* or a *Z* boson and MET [394], among others. These searches exclude masses of stops and sbottoms up to ${m}_{\stackrel{~}{t}}~1$ TeV and ${m}_{\stackrel{~}{b}}~900$ GeV for some simplified models.

In the cases where the colored sector of a supersymmetric model has large masses, the direct production of chargino, neutralino and slepton pairs dominate. A pair of directly produced sleptons decay typically like $\stackrel{~}{l}\to l{\stackrel{~}{\chi}}_{1}^{0}$. Neutralinos and charginos are produced in pairs in a number of different combinations, the most commonly studied of which are ${\stackrel{~}{\chi}}_{0}^{2}{\stackrel{~}{\chi}}_{1}^{\pm}$ and ${\stackrel{~}{\chi}}_{1}^{\pm}{\stackrel{~}{\chi}}_{1}^{\pm}$. The decays of heavy neutralinos and charginos produce *W*, *Z* or Higgs bosons and the lightest neutralino. Further decay of *W* and *Z* sets the final states targeted by ATLAS and CMS searches, such as the final state with two leptons and MET [395–397], many leptons and MET [397–400], leptons, jets and MET [393], taus and MET [401, 402], and b-jets plus MET [403], among others. Due to their low production cross sections, the exclusion limits on slepton masses from direct production are quite weak and they only reach up to around ${m}_{\stackrel{~}{l}}~500$ GeV. Stronger limits on slepton masses can be inferred from neutralino/chargino production with sleptons in the cascade, reaching up to ~850 GeV. The limits on electroweakinos (neutralinos and charginos) are very sensitive to the parameter choices for the simplified model analyses performed by the experiments, hence the exclusion limits on ${\stackrel{~}{\chi}}_{2}^{0}$ and ${\stackrel{~}{\chi}}_{1}^{\pm}$ vary from search to search and from signal region to signal region, roughly in the range ${m}_{{\stackrel{~}{\chi}}_{0}^{2}},{m}_{{\stackrel{~}{\chi}}_{\pm}^{1}}\in (500\mathrm{\text{GeV}},1.1\mathrm{\text{TeV}})$. Similarly the limits on the lightest neutralino varies in the range ${m}_{{\stackrel{~}{\chi}}_{1}^{0}}\in (200,700)$ GeV. In addition, there is a hint of an excess in one of the two and three lepton final state analyses by the ATLAS collaboration in the low ${m}_{{\stackrel{~}{\chi}}_{1}^{0}}$ region, with a reported significance of 2 and 3σ in the 2*l* and 3*l* channels, respectively [399].

In addition to sparticle searches, SUSY can also be probed through searches for heavy and charged Higgs bosons. The MSSM predicts the existence of two CP-even scalars, *h* and *H*, one CP-odd pseudoscalar, *A*, and a charged scalar *H*^{±}. The lightest CP-even scalar, *h*, is said to be “SM-like,” as its mass and couplings are aligned with the Higgs boson discovered at the LHC [4, 5], the so called *alignment limit*. Neutral heavy Higgses can be produced at the LHC in the same manner as the SM Higgs, that is by gluon fusion, vector boson fusion (VBF) and associated production, with a *t* and/or *b* quark. Thus, the same mechanisms that lead to the discovery of the SM Higgs are employed to set exclusion limits on heavier neutral (pseudo)scalars, including signatures where *H* is produced in resonance and decays into two light Higgs bosons *pp* → *H* → *hh*. The final states targeted by these exotic Higgs searches consist of 2–4 leptons, jets and MET from on- and off-shell *W* and *Z* bosons [404, 405], two leptons and MET [406, 407], final states with four b-jets [408, 409] or two b-jets and *WW* [410], γγ [411, 412], ττ [413, 414], μμ [415], or $t\overline{t}$ [416] decays, ditau final states [417, 418] and diphotons, with and without associated *W* bosons [419–421]. Charged Higgs bosons, *H*^{±} can typically be produced with associated resonant and non-resonant top-quark production. Their main decay channels are *H*^{±} → *W*^{±}*Z* [422, 423], *H*^{±} → *t*(*c*)*b* [416, 424] and *H*^{±} → τ^{±}ν [416, 425]. Since no excess has been found for either heavy *H* or *H*^{±}, the experiments set upper limits that strongly depend on the production cross section and, in turn, on tan β. For *H* the limits range from *m*_{H} < 400 GeV for tan β = 2 and production cross sections limits of 0.1 pb for larger masses. For *H*^{±} with ${m}_{{H}^{\pm}}<160$ GeV excluded for all values of tan β and ${m}_{{H}^{\pm}}<1.1$ TeV excluded for tan β = 60.

Many unified theories automatically preserve *R*-parity, such as left-right symmetric or Pati-Salam models, as well as supergroups of them, *SO*(10) or *E*_{6}. This is because they contain a gauged *U*(1)_{B−L} subgroup which effectively makes the LSP stable [76]. Other models, such as *SU*(5), may have *R*-parity violating (RPV) interactions, though in general they will be suppresed since they can lead to rapid proton decay. Since the LSP is no longer stable, collider signatures of RPV typically contain multiple leptons [400], multiple jets [426–429] or both [430, 431] in the final state. These searches often impose strong upper limits on sparticle masses that range from 150 GeV to a few TeV for $\stackrel{~}{t}$, depending on the channel, and from 1 to 2 TeV, for $\stackrel{~}{g}$.

If the LSP is metastable or the lightest chargino and neutralino are almost degenerate, they can live long enough to leave a displaced vertex or a disappearing track on the detector. Detailed searches have been performed by ATLAS [432–436] and CMS [437–439] to search for these long-lived particles, and they have reached exclusion limits comparable to those of the detailed searches above.

Most of the searches described above assumed a neutralino LSP, which is typically the case in gravity mediated SUSY breaking. In gauge mediated SUSY breaking (GSMB) and general gauge mediation (GGM) the LSP is actually a nearly massless gravitino. In these cases new decay channels are open with photons [440–443], Z's [400, 403], Higgses [400, 403, 444], and τs [385, 445] in the final state (see Figure 11).

Although the LHC results are the most recent and, for the most part, they supersede the results of previous collider experiments, such as those at the Tevatron, some experimental limits from LEP still remain relevant today. In particular for models with a significant production of neutralino/chargino or slepton pairs, the limits from ALEPH [446, 447], OPAL [448, 449], L3 [450, 451], and DELPHI [452] on sleptons and electroweakino masses are rather pertinent, as they are largely model independent.

Many of the SUSY searches above are performed using simplified models, e.g., ATLAS jets and MET search assumes a 50–50 split between the decay modes of gluinos [372]. Therefore, the mass and cross section limits obtained are weakened in more complicated models. In order to assess the relevance of many of these exclusion limits on several popular SUSY models, a full global fit of the parameter space of the model is required. Several of these fits have been performed for the CMSSM, NUMH1 and NUMH2 [453–459], phenomenological MSSM models [460–462], SUSY GUT models [463] and electroweak-sector MSSM models [116], by the Zfitter [464]^{12}, SuperBayes [465, 466]^{13}, Fittino [467]^{14}, MasterCode [468]^{15}, and GAMBIT [469–476]^{16} collaborations.

### 6.2. Collider Searches for Leptoquarks

Leptoquarks (LQs) are associated either with the vector (spin 1) particles that correspond to the gauge bosons of the unified gauge group or they can be scalars (spin 0) and belong to a Higgs sector of a unified theory. Vector LQ mass is typically of the order of the unification scale and can only be accessible directly at colliders if the unification scale is low enough (e.g., Pati-Salam models). Scalar representations can also contain light fields, most notably the SM Higgs, but come at the cost of severe fine tuning, as discussed in the section 2 on the example of doublet-triplet splitting problem in 5-dim. representation of *SU*(5). In non-supersymmetric unified models the presence of light colored scalars tends to aid unification (see e.g., [477]). Another important difference is that the scalar LQ interactions can be analyzed without specifying the concrete GUT completion in the ultraviolet. On the other hand, vector (gauge) LQs are sensitive to the mass generation mechanism that is specified in the ultraviolet. Therefore, effective vector LQ models are not renormalizable [478]. Furthermore, even the couplings of vector LQs to the SM gauge sector are not completely fixed by the gauge quantum numbers [479, 480].

Altogether there are six scalar and six vector leptoquarks, listed in Table 1, that couple to the SM matter at the renormalizable level [480–482]. The fermionic number *F* ≡ 3*B*+*L* of leptoquarks that do not couple to diquarks (ϕ*qq*) and are potentially *B* and *L* conserving, must be *F* = 0, whereas LQs with |*F*| = 2 in general destabilize the proton via *B*−*L* conserving decays.

The most important phenomenological characteristic of light LQs (of mass of the order few TeV) is their color triplet nature allowing them to be produced in pairs via strong interactions in a largely model independent manner. In this section we will focus on the on-shell production of LQs at *pp* colliders, since the current mass constraints are dominated by LHC searches. For specific signatures of LQ production in colliders with alternative initial states (see [480]). Pair production of leptoquarks is model independent for the $gg\to \text{LQ}\overline{\text{LQ}}$ partonic process, while the parton level process $q\overline{q}\to \text{LQ}\overline{\text{LQ}}$ is affected also by the *t*-channel lepton exchange diagram (bottom right diagram in Figure 12) that introduces some model dependence when the leptoquark flavor couplings are non-negligible. The partonic cross sections for pair production at leading order are [483–488]:

where ŝ is the partonic center-of-mass energy squared, α_{3} is the strong coupling constant, and $\beta =\sqrt{1-4{m}_{\varphi}^{2}/\u015d}$. A model independent study of weak doublet scalar LQs at the LHC and the interplay with low energy flavor processes was performed in [489], where it was also shown that searches for single LQ production could be more sensitive in the regime of large Yukawas and/or LQ masses. An analysis of pair and single production, along with the corresponding UFO model files LQ_NLO, both for scalar and a vector leptoquark has been presented in Doršner and Greljo [490].

**Figure 12**. Representative diagrams for leptoquark pair production at *pp* colliders. Dots denote the LQ-*q*-ℓ coupling.

On the other hand, single leptoquark production at *pp* colliders is always model dependent (Figure 13). Single leptoquark searches are more effective at larger LQ masses [489].

**Figure 13**. Representative diagrams for single leptoquark production at *pp* colliders. Dots denote the LQ-*q*-ℓ coupling.

On the decay side of the process, experimental searches for pair and single LQ production are targeting a resonance in the *jℓ* channel. The decay width of a scalar leptoquark into a lepton-quark final state is given by [490] and [491]:

Current bounds from dedicated leptoquark pair production have been commonly extracted in the framework that assumed LQ coupling only to a single generation of SM fermions, whereas realistic LQ scenarios could posses richer flavor structure [480, 492]. The experimental upper bounds are given for the product of cross section and the LQ branching fraction probability β^{2}, where β is the probability for LQs to decay to a final state with charged leptons. There have been numerous analyses performed at the LHC for leptoquarks being either of 1st [493], 2nd [494], or 3rd [495–499] generation. More recent studies, motivated by the observed lepton universality violation in *B*-meson decays, allow also for cross-generational couplings, e.g., [500]. Finally, also single LQ production channels are being studied [501]. To conclude, the current lower bounds on leptoquark masses from direct searches at the LHC range from several 100 GeV to above the TeV scale, where the exact bound depends on the size of the flavor couplings.

### 6.3. Other Exotic Searches

Many unified theories include heavy sterile neutrinos that contribute to the mass of the light neutrinos via type-I seesaw mechanism (see section 3.2). These are often associated with the symmetry breaking of a left-right sector of the theory, and thus they are expected to be heavier than the EW scale. Direct searches for heavy neutrinos at colliders often target a decay channel where final state has two same-sign leptons, via *s* or *t*-channel production of a gauge boson that can be left or right handed [502]. Figure 14 shows the Feynman diagram for the golden channel for heavy neutrino searches, *pp* → *W* → *Nl* → *Wll* → *lljj*. The ATLAS and CMS experiments at the LHC have performed searches for heavy neutrinos in LR models with masses *M*_{N} ≈ (20, 1, 600) GeV and have imposed strong limits on the couplings between active and sterile neutrinos [503–505].

Other searches for sterile neutrinos are performed in beam dump experiments [506], where the neutrinos are produced in semileptonic decays of mesons, with masses below 2 GeV [507–510]. For intermediate masses below the *Z* resonance, the strongest limits come from the decay *Z* bosons by the LEP experiments DELPHI and L3 [511, 512].

In addition to singlet fermions and colored leptoquarks, the LHC looks for heavy colorless vector bosons as part of their exotic searches programme. Charged *W*′ and neutral *Z*′ vector bosons are predicted in a number of GUT frameworks and they can often live at low scales, which positions them within the reach of colliders. Clear examples of this are the left-right symmetric models described in section 4.1, that predict light *W*_{R} and *Z*_{R} bosons, or the light *Z*_{N} appearing in E_{6}SSM models.

These states are produced at *pp* colliders through Drell-Yan processes *pp* → *W*′/*Z*′ and subsequently decay into leptons or jets. One of the most targeted processes for *W*′ involve the decay into heavy neutrinos, as in Figure 14, with two same or opposite sign leptons (depending on the Majorana or Dirac nature of the heavy neutrinos) and jets [503, 504]. These searches often use a simplified model where *g*_{R} = *g*_{L} and ${M}_{N}={M}_{W}^{\prime}/2$ resulting in high exclusion limits with ${M}_{{W}^{\prime}}\gtrsim 4.5$ TeV, but it has been shown that these limits weaken somewhat in more general models [173–175]. CMS also reported a search for *W*′ where the vector boson decays to $\tau {\overline{\nu}}_{\tau}$, the τ decaying hadronically [513], with slightly weaker limits.

Narrow resonance searches for *Z*′ have been performed by ATLAS and CMS, targeting final states with two opposite-sign leptons. These searches have yielded model-dependent exclusion limits on ${M}_{{Z}^{\prime}}$. For *E*_{6}-inspired *Z*′, the limits vary around ${M}_{{Z}^{\prime}}\gtrsim (3,3.5)$ TeV, whereas for LR models they are moderately stronger ${M}_{{Z}^{\prime}}\gtrsim 4$ TeV [514, 515].

Lastly, GUTs predict a plethora of different scalar states that can be observed at the LHC if they are light enough, e.g., Δ_{L, R} in LR symmetric models. Searches for neutral and singly charged scalar bosons are identical to the searches for supersymmetric Higgs bosons in section 6.1, so we will not repeat them here. Doubly-charged scalars, such as the ${\delta}_{L}^{\pm \pm}\in {\Delta}_{L}$ in LR models, have been studied by ATLAS and CMS in multilepton final states [516, 517], diboson final states [518] and in long-lived particle studies [519] with model-dependent limits below 1 TeV.

## 7. Precision Tests of Unification

### 7.1. Proton Decay

Unified theories may contain gauge or scalar bosons that mediate transitions between leptons and quarks. These transitions violate baryon *B* and lepton *L* number separately and hence can cause the rapid decay of nucleons [520–529].

In the language of Effective Field Theory, nucleon decay transitions are generated by higher dimensional operators suppressed by the mass of the heavy mediator. The most relevant contribution to nucleon decay comes from dimension-6 operators of the form *qqql*, mediated either by a gauge or scalar boson. In SUSY GUTs, however, dimension 4 and 5 operators can appear, involving *R*-parity violating interactions and mixing among sfermions, respectively. Assuming the conservation of *R* parity and minimal flavor violation (MFV) in the sfermion sector, however, dimension 4 and 5 contributions can be made negligible [530–532]. Therefore, dimension 6 operators dominate the contributions to nucleon decay. These are, in general, model dependent, so calculating the decay width requires knowledge of the flavor structure at the GUT scale and varies among the different decay channels [533, 534]. However, with some simplifying assumptions one can approximate the decay width of the proton as [535]

There are several decay channels for the proton and neutron, each with a different experimental bound. The most stringent of them, known as the *gold channel* for proton decay is *p* → *e*^{+}π^{0}, whose parton-level diagrams can be seen in Figure 15, and with a lower limit on the half-life, set by the Super-Kamiokande, of τ > 1.6 × 10^{34} years [43]. Other processes with slighly lower bounds are τ(*p* → μ^{+}π^{0}) > 7.7 × 10^{33} years [43], τ(*p* → ν*K*^{+}) > 5.9 × 10^{33} years [536] and τ(*nn* → *e*^{±}μ^{∓}) > 4.4 × 10^{33} years [537].

**Figure 15**. Parton level *s* and *t*-channel diagrams for the proton decay channel p → e^{+}π^{0} for a gauge or scalar boson mediator *X*.

The next generation of experiments for nucleon decay has already been proposed. Hyper Kamiokande will take the place of Super-K and has a projected sensitivity in the golden channel *p* → *e*^{+}π^{0} of ~10^{35} years [538]. The Deep Underground Neutrino Experiment (DUNE) [539] expects to improve the limit on *p* → ν*K*^{+} to ~ 3 × 10^{34} years. These increased limits will probe unified teories at the highest scales and could be the smoking gun for them. In case of a positive signal from either of these experiments, precision calculations of proton decay processes with controlled uncertainties will become invaluable [84, 540].

### 7.2. Flavor Phenomenology of Light Leptoquarks

Several flavor couplings of scalar leptoquarks to leptons and quarks, which are in general connected to the GUT contractions of scalar and fermionic representations, allow their virtual effects to be tested in low-energy flavor observables. Such are the decays of hadrons or leptons, precision observables at LEP, and static properties of particles such as dipole moments. On the high-*p*_{T} front, the LHC is also becoming competitive as a flavor probe for virtual effects of particles that are too massive to be produced on-shell.

The correlations between lepton-quark-LQ couplings are determined at low scales by the weak isospin and hypercharge. As an example, consider the weak doublet leptoquark *R*_{2}(3, 2, 7/6) (see Table 1), which can couple to two-types of lepton-quark bilinears:

Here *i, j* are the flavor indices, primed fields are written in the interaction basis, unprimed fields are in the mass basis, except for the neutrinos which are aligned with charged leptons. There are three important features in the above Lagrangian. First, since *R*_{2} is a weak doublet it must couple to left-handed quark doublets, which implies that CKM matrix *V* relates the couplings of up-type and down-type quarks. Second, *R*_{2} couples to both chiralities of quarks and leptons, which leads to scalar and/or tensor effective interactions and could lead to enhanced effects in meson mixing amplitudes, dipole moments, and radiative decays [480]. Third, as *F* = 0 for *R*_{2} we cannot construct interaction term with diquark coupling, implying that proton cannot decay via Δ(*B*−*L*) = 0 process. Generalizing to other LQ states, weak triplets only talk to the left-handed fermions (2⊗2), leading to strictly chiral interaction, whereas singlet LQs can talk to 2⊗2 and 1⊗1 fermion bilinears.

Among the flavor constraints, leptoquarks naturally (at tree-level) contribute to semileptonic effective operators at low scales, therefore the most relevant observables are (semi-)leptonic decays of mesons, baryons, or τ leptons. The most notable charged-current and flavor changing neutral current (FCNC) constraints, and the general framework to address them in leptoquark models, have been spelled out in Leurer [541], Davidson et al. [482], and Doršner et al. [480]. The most constraining are the FCNC observables, where the tree-level LQ contribution can easily stand out of the SM signal, which is 1-loop suppressed in the case of quark FCNC and absent in the case of lepton FCNC. Effective dimension-6 interactions for 4 lepton or 4 quark operators, which drive the ℓ → ℓ′ℓ′ℓ″ (see section 7.4) and meson mixing processes, occur at one-loop [480]. Therefore, meson mixing is in general not among the strictest constraints on LQs (for *B*_{s} mixing see e.g., [542]).

### 7.3. Lepton Flavor Universality

Lepton flavor universality (LFU) ratios, defined as ratios between rates for processes that differ only in lepton flavor, are very well suited to test the validity of the SM. The main advantage is that in the Standard Model LFU is respected by all gauge interactions, the only breaking comes from mass splitting among leptons, which leads to efficient cancellation of hadronic and parametric uncertainties in LFU ratios. Recently, two LFU ratios in *B*-meson decays have been observed

where *l* = *e*, μ. Several experiments found that the ratios ${R}_{{D}^{(*)}}$ are larger than ${R}_{{D}^{(*)}}^{\text{SM}}$. The measurements of *R*_{D} [13, 543, 544] differ by ~ 2 σ with respect to the SM prediction [545] and by ~ 3 σ in the case of ${R}_{{D}^{*}}$ [546–548]. Combined significance reaches 4σ deviation from the SM [549]. The LHCb experiment has also measured ${R}_{{K}^{(*)}}$ LFU ratios, related to the neutral-current process *b* → *sll*, and found them to be lower than expected in the SM. While *R*_{K} was measured in a single kinematical region, *q*^{2} ∈ [1.1, 6]GeV^{2} [12], ${R}_{{K}^{*}}$ was measured also in the ultra-low region *q*^{2} ∈ [0.045, 1.1]GeV^{2} [550]. Each of the ${R}_{{K}^{(*)}}$ measurements is ~ 2.5σ below the SM prediction level [551, 552], and furthermore, there are discrepancies in *b*→*sℓℓ* driven decays that are coherent with the deviation in ${R}_{K}^{(*)}$ if there is ~ 20% reduction in the vector Wilson coefficient *C*_{9} [553, 554].

Light leptoquarks are prime candidates to explain one or both of those puzzles. For the ${R}_{{D}^{(*)}}$ the effective Lagrangian contains four relevant operators:

Model independently it has been shown that ${R}_{{D}^{(*)}}$ can be explained either by rescaling the SM semileptonic operator (*g*_{VL}), by turning on *g*_{T}, or by particular combinations of scalar and tensor operators, *g*_{SL} = ±4*g*_{T} operators, that arise in presence of a non-chiral LQ [555–558]. In order to address ${R}_{{K}^{(*)}}$ one has to modify the vector Wilson coefficient *C*_{9} whereas the axial Wilson coefficient *C*_{10} may also be present in the effective Lagrangian:

Also purely left-handed scenarios with *C*_{9} = −*C*_{10} ≈ −0.6, which are characteristic of LQ weak-singlet or triplet exchange, are in good agreement with ${R}_{{K}^{(*)}}$ and the global fit of *b* → *sℓℓ*. Such left-handed leptoquark solutions have been put forward: triplet scalar *S*_{3}, singlet vector *U*_{1}, triplet vector *U*_{3} [557, 558]. For loop-level explanation of ${R}_{{K}^{(*)}}$ one can also invoke singlet *S*_{1} [559] or doublet *R*_{2} [560], but at the price of large couplings. There are several proposals with scalar leptoquarks that address ${R}_{{K}^{(*)}}$ and/or ${R}_{{D}^{(*)}}$ [559, 561–566], some in the context of unified theories such as *SU*(5) [567–569], left-right symmetry [570], Pati-Salam [571, 572], *SO*(10) [573] and others [574]. Recently it was realized that a singlet vector leptoquark ${U}_{1}^{\mu}(3,1,2/3)$ generates left-handed interactions and partially resolves both LFU puzzles, in many UV frameworks [575–577], including three-flavor extensions of the Pati-Salam model [578–584].

Finally, moderately large leptoquark couplings dictated by the above LFU anomalies can be also observed in processes with virtual LQ exchange, typically in the *t*-channel, resulting in a final state with at least one charged lepton. Inspired by the abovementioned LFU anomalies, processes with final state leptons have been studied, which can probe LQ scenarios for LFU violation observed in *B* meson decays [585–588]. In this case, LQ cannot be produced on-shell and the sensitivity does not deteriorate abruptly with rising LQ mass. Instead there is a smooth transition to the effective theory picture, where heavy LQ is integrated out. Thus, among the LQ induced processes the *t*-channel has the best mass reach for LQs and it is thus complementary to pair and single production [489, 490]. Another recent set of observables at the LHC, targeting the LQ scenarios that are well suited to explain the LFU anomalies ${R}_{{D}^{(*)}}$, are the searches with a single τ lepton in the final state [589]^{17}. Third generation leptoquarks could also be probed in $t\overline{t}$ final states [591].

### 7.4. Lepton Flavor Violation and Dipole Moments

In the SM with massive neutrinos, lepton flavor violation (LFV) can occur via the mixing in the neutrino sector. It is, however, heavily suppressed due to the GIM mechanism [592], as the rate depends on the neutrino masses resulting in an unobservable prediction of order 10^{−55} [593]. Extensions of the SM modify this prediction by introducing additional sources of lepton flavor violation [594–596]. New physics can then be probed by testing the deviations of certain lepton flavor violating processes with respect to the experimental limits.

Charged lepton flavor violating processes are typically of three types: ${l}_{\alpha}^{-}\to {l}_{\beta}^{-}\gamma $, ${l}_{\alpha}^{-}\to {l}_{\beta}^{-}{l}_{\gamma}^{-}{l}_{\delta}^{+}$ (with α ≠ β) and μ −*e* conversion in nuclei [597]. One-loop contributions to the first two processes can occur through a dipole and box diagrams as depicted in Figure 16, with a scalar or vector mediator(s) *X*(*X*′) and a SM or exotic fermion(s) *f*(*f*′) running in the loop. Contributions to μ − *e* conversion follows from the penguin (center) and box (right) diagrams with *l*_{α} = μ, *l*_{β} = *e* and *l*_{γ, δ} = *q*.

**Figure 16**. Diagrams contributing to LFV processes, ${l}_{\alpha}^{-}\to {l}_{\beta}^{-}\gamma $ (left) and ${l}_{\alpha}^{-}\to {l}_{\beta}^{-}{l}_{\gamma}^{-}{l}_{\delta}^{+}$ (center and right) with *X* and *X*′ scalar or vector mediators and *f* and *f*′ fermions.

Unified theories often contain a number of exotic states capable of fulfilling the role of *X* and *f* in Figure 16, violating lepton flavor either by interactions between the leptons and mediators or by mixing in the leptonic sector. The latter case is realized in GUT models with heavy neutrinos (c.f. left-right models in section 4.1), where the mixing between active and sterile neutrinos enhances the LFV contribution [598–600]. The contribution to the branching ratios of the most constraining LFV processes, μ → *eγ*, μ → *eee* and μ −*e* conversion, in these models, with heavy neutrinos of mass *M*_{NI}, active-sterile mixing Θ_{αI}, a right-handed gauge boson *W*_{R} and left and right-handed scalar triplets δ_{L} and δ_{R}, can be written as [601], [602], and [174]

In supersymmetric GUTs there are many possible sources of lepton flavor violation, parametrised by the mixing in the slepton sector of the MSSM, which has contributions to LFV processes of the type [603]

Off-diagonal entries in the slepton mass matrices can be the result of non-minimal flavor violating interactions or non-canonical Yukawa textures at the GUT scale, where the soft masses are supposed to unify [604, 605]. In addition, slepton mixing can be induced in minimal flavor violating (MFV) SUSY via seesaw mechanisms [606–612] or, for Yukawa-unified theories (see section 3.3), it can depend on the CKM matrix at the GUT scale [596]. LFV contributions can also arise in SUSY models where *R*-parity is violated, explicitly or spontaneously, with interaction terms of the type ${l}_{i}{l}_{j}{\stackrel{~}{\nu}}_{k}$ that induce tree-level contributions to *l* → *lll* decays and μ → *e* conversion, as well as new dipole contributions to *l* →* lγ* [613–617].

The anomalous electric, *d*_{i}, and magnetic, *a*_{i}, dipole moments of quarks and leptons follow from processes identical to the diagram on the left in Figure 16, where *l*_{α} and *l*_{β} have the same flavor. Hence contributions from heavy states running in the loops can have a strong effect that can be tested experimentally. As with LFV, the SM contribution to electric dipole moments (EDMs) is tiny, as it is proportional to the CP-violation phase in the CKM matrix [618, 619]. EDMs have not been observed experimentally, so deviations from the SM prediction due to CP-violation in other sectors is strongly constrained [620]. Other sources of CP violation can appear in neutrino mixing [619], phases in fermion-sfermion couplings [621] or extended Higgs sectors [622]. Anomalous magnetic moments (AMM), on the other hand, have been measured with extreme accuracy. In fact, the precision of both the experimental measurement and theoretical prediction for *a*_{μ} has shown a discrepancy of more than 3 standard deviations [34]. New physics contributions have been shown to resolve that tension, particularly in the context of supersymmetry [623].

In the presence of light leptoquarks anomalous dipole moments of leptons or quarks are one-loop processes [480]. A special feature of non-chiral leptoquarks, such as *R*_{2} with couplings (61), is that both *l*_{ij} and *r*_{ij} are non-zero in the interaction Lagrangian ${\overline{q}}^{i}\left[{l}_{ij}{P}_{R}+{r}_{ij}{P}_{L}\right]{\ell}^{j}\varphi $ which then leads to the anomalous moment of the muon:

where *m*_{ϕ} and *Q*_{ϕ} are the charge and mass of the leptoquark and *q* is the flavor of the quark in the loop (see the leftmost diagram in Figure 16). Shown is the leading order contribution in *m*_{q}. The first term increases *a*_{μ} only when *Q*_{ϕ} > 2/3 and it is present for all scalar LQ states that couple to a muon. The second term is relevant for non-chiral LQs and it is chirally enhanced by *m*_{q}/*m*_{ℓ}, possibly leading to large effects with moderate couplings to *b* or *t* quark. Furthermore, the sign of the non-chiral term is adjustable. On the other hand, the same mechanism also enhances dipole LFV transitions, e.g., μ → *eγ*, τ → μγ [560]. Non-chiral LQs may also generate quark or lepton electric dipole moments at 1-loop [624, 625].

Whichever the mechanism, it is clear that LFV and anomalous dipole moments are predicted by many GUT models, with varying strengths for different processes. The full list of processes and their current experimental upper bounds and measurements can be seen in Table 2, where the experiments that have studied the processes are detailed. Furthermore, new experiments are being developed at this moment that will attempt to improve the limits on processes like μ → *eee* (Mu3e [642]) and μ − *e* conversion (COMET [643], Mu2e [644]), with projected limits up to four orders of magnitude stronger than previous studies. Additionally, a new measurement of *a*_{μ} has been performed by the Muon g-2 experiment and it is expected to be released soon [645], which may confirm the deviation observed before, and thus further motivate the need of new physics.

**Table 2**. Upper bounds at 90% C.L on LFV processes and EDMs, and measurements of AMMs, along with the experiments that provided them.

### 7.5. Neutrinoless Double Beta Decay

This rare nuclear process corresponding to a simultaneous conversion of two neutrons to two protons and two electrons within the nucleus [646] is of great interest for particle physics, as it clearly does not conserve lepton number, and thus violates the corresponding accidental Abelian global symmetry of the SM. Consequently, a strong experimental effort is being made to observe this unique process. Unfortunately, its observation will be very difficult, as 0ν*ββ* decay is expected to be extremely rare. There is a number of experiments, some in operation, other being constructed or planned, attempting to measure the decay. An overview of the major collaborations is shown in Table 3. The current experimental lower limits on its half-life are around 10^{26} years [647, 648] and the future searches should reach sensitivities by one or two orders of magnitude higher.

As can be shown, the existence of 0ν*ββ* decay implies the Majorana nature of neutrinos (and vice versa) [649–651] and as such it represents one of the best probes of this BSM hypothesis. Therefore, in GUT models allowing for Majorana neutrino mass generation 0ν*ββ* decay can be in principle always triggered. This, however, does not say anything about the underlying mechanism and the resulting 0ν*ββ* decay rate. The prominent *standard (mass) mechanism* of 0ν*ββ* decay assumes a light Majorana neutrino exchange between the two beta-decaying neutrons. In the SM with light massive neutrinos this process can be depicted as shown in Figure 17, left. Besides the standard scenario a number of *non-standard mechanisms* triggering 0ν*ββ* decay can be constructed. The effective treatment of these exotic mechanisms can be conveniently employed (see e.g., [652–656]).

**Figure 17**. The standard mass mechanism of 0ν*ββ* decay **(left)** and a non-standard contribution that can be constructed in the left-right symmetric models using vector currents of opposite chiralities **(right)**.

As for the UV-complete 0ν*ββ* decay mechanisms, a variety of interesting ones can be constructed within GUTs. For instance, in the left-right symmetric models (where, of course, the standard light neutrino exchange is available) one can think of several exotic mechanisms involving exchange of heavy neutrino as well as light and heavy *W* vector bosons [657]. In the simplest exotic case the light neutrino exchange is substituted by a heavy right-handed neutrino exchange, which means that the involved vector currents and emitted electrons must be also right-handed. Due to the large mass of the propagating neutrino, the interaction can be considered to be contact and we refer to this contribution as to short-range mechanism. Since the right-handed currents are present in left-right symmetric models, it is also possible to draw 0ν*ββ* decay mechanisms, in which the neutrino exchange does not violate chirality. This means that the contribution is not proportional to the neutrino mass and the two outgoing electrons are of opposite chiralities. A possible mechanism of this type is depicted in Figure 17, right. As apparent, the diagram involves one right-handed and one left-handed vector current and since the light neutrino propagator is present, one refers to this contribution as to a long-range 0ν*ββ* decay mechanism.

Leptoquarks, particles appearing prominently in GUTs, can also trigger non-standard 0ν*ββ* decay contributions. It has been described that this is the case, when different leptoquark multiplets mix via a possible leptoquark-Higgs coupling violating lepton number [658, 659]. Diagrams of this type of contributions to 0ν*ββ* decay are shown in Figure 18. The specific helicity structure of the effective four-fermion interaction leads to the fact that this contribution can dominate over the standard mass mechanism. The current lower limits on 0ν*ββ* decay half-life then allow to derive the bounds on corresponding leptoquark parameters.

Neutrinoless double beta decay can be triggered also in supersymmetric theories aspiring for grand unification. In the simplest case, if the MSSM with broken *R* parity is considered, 0ν*ββ* decay diagrams involving *R*-parity-violating couplings and supersymmetric mediators can be drawn [660–664]. An example of a supersymmetric 0ν*ββ* decay mechanism is depicted in Figure 19. Again, from non-observation of 0ν*ββ* decay it is possible to derive limits on the unknown model parameters.

## 8. Outlook and Future Prospects for GUTs

Among the vast landscape of theories beyond the Standard Model, Grand Unified Theories stand out as appealing candidates. As we have seen, GUTs are a collection of ideas from group theory, supersymmetry, neutrino physics, flavor physics and more, which positions them as some of the most complete and attractive theories in the literature. Indeed they are among the few BSM theories capable of simultaneously affecting the highest energy scales, influencing the cosmology of the early Universe, and the low energies, within reach of colliders and terrestrial experiments.

Throughout this review article we have provided a rough sketch of the status of GUTs and some of the associated research in the field. We have described the basic principles behind them and their first appearance in the world of particle physics in the late 70s. A lot of effort was spent in the years after and many of the greatest models were designed at that time. Research in GUTs has continued since, focusing either on specific topics within and alongside them, or on particular models that compiled a few advances in the field. We have thus summarized a non-exhaustive selection of topics and models that are at the forefront of research in unified theories, aiming to provide an overview of the current state of the art.

We are fortunate enough to live in a time where experimental searches are abundant and they cover a rather vast range of fronts. The most cutting-edge technologies have been and are being developed to push the boundaries of our current understanding of particle physics and cosmology. Grand unified theories are and will be put under the microscope by many of these experimental advances, which will confirm, constrain or outright exclude some of the existing models.

The recent observation of gravitational wave signatures opens a new window into the history of the Universe, where events and phenomena that ocurred in the early Universe can be observed with gravitational wave detectors. Cosmic phase transitions associated with patterns of symmetry breaking in unified theories are such events, as they can be the source of stochastic gravitational waves that can be observed today. Transition temperatures above the EW scale, typically associated with the breaking of some intermediate step in a GUT model, can be studied by future gravitational wave experiments such as a LISA, the Einstein telescope, Kagra, the Cosmic Explorer, BBO and DEIGO. Also in the cosmological frontier, GUTs can have a serious impact on the inflationary epoch of the Universe, testable in measurements of the CMB, and can contribute to the baryon asymmetry of the Universe, via baryo and leptogenesis.

At the time of writing we have reached the end of the second run of the LHC, with an outstanding recorded integrated luminosity of about 150 fb^{−1}. Analyses of the accumulated data, however, are still under way and they will probably spill well into the start run 3 in 2021. Many of the analyses already published have strong consequences for the predictions of unified theories, as are direct searches for supersymmetry, leptoquarks or other exotics at ATLAS and CMS. Upcoming results from ongoing and future analyses of the results from the LHC experiments may strengthen the bounds on light states as predicted by SUSY GUTs and other models, or they might show hints of the existence of new particles, whose relevance for GUTs would need to be determined. Upgraded versions of the LHC (HL-LHC, VLHC or FCC) or other future colliders (ILC, CLIC) will certainly boost this programme with increased accuracy and higher energies, which will further probe the low-hanging states predicted by GUTs.

Where colliders search for the low scale predictions of GUTs, precision experiments can explore the intermediate and high scales associated with unification. Nucleon decay limits are often among the strongest probes of fully unified theories and future experiments such as Hyper-Kamiokande and DUNE may set even stronger exclusion limits or perhaps measure signs of proton decay, which would be a smoking gun for GUTs. Furthermore, GUTs can provide contributions to a number of flavor and precision observables, such as LFU, LFV, EDM, AMM, or 0ν*ββ*, some of which are in tension with the SM. Confirmation of these flavor anomalies with more collected data by LHCb and other experiments would be undeniable evidence of the need for new physics models and GUTs are very well suited for that purpose.

To conclude, Grand Unified Theories are still at the vanguard of research in BSM models. They can explain many of the issues of the SM and can accommodate the recent results from the cosmological, precision and collider frontiers with relative ease. Contrary to “simplified” models, GUTs are complete theories that can simultaneously make a large number of testable predictions on the different fronts. Fortunately, these predictions can be explored by upcoming analyses and future experiments, which can set strong exclusion limits in a subset of GUT models. On a more optimistic note, any observation in, for instance, SUSY searches at colliders, gravitational waves signatures or proton decay will stack the odds in favor of some GUT models and will significantly shape the future of the research in particle physics.

## Author Contributions

TG has lead the development of this review. All other authors have contributed equally.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

The authors would like to thank P. Athron, F. Deppisch, A. Kvellestad, and Y. Zhang for helpful discussions. TG was partly funded by the Research Council of Norway under FRIPRO project number 230546/F20 and partly supported by the ARC Centre of Excellence for Particle Physics at the Tera-scale, grant CE110001004. NK was supported by the Slovenian Research Agency under the research core funding grant P1-0035 and No. J1-8137. TRIUMF receives federal funding via a contribution agreement with the National Research Council of Canada and the Natural Science and Engineering Research Council of Canada.

## Footnotes

1. ^See section 3.2 for details on neutrino mass generation through the seesaw mechanism.

2. ^A detailed description of unification in SUSY *SU*(5) can be found in Dorsner et al. [51].

3. ^Although successful in achieving spontaneous symmetry breaking (SSB) in this SUSY LR model, the resulting vacuum state does not preserve the electromagnetic charge.

4. ^See section 6.1 for more details on searches for Supersymmetry.

5. ^Similarly to the right-handed neutrino singlets, only two triplets are necessary, although three (one per flavor) are considered here.

6. ^Here, the representations are labeled the usual way in the order {*SU*(3)_{C}, *SU*(2)_{L}, *SU*(2)_{R}, *U*(1)_{B−L}}.

7. ^However, there are exceptions, most notably the flipped *SU*(5) SUSY GUT theories.

8. ^We expand more on the topic of the baryon assymmetry in subsection 5.4.

9. ^Note that the existence of two time scales in the spectral form means that the peak amplitude for the turbulence contribution cannot be found simply be setting the frequency to either *h*_{*} or *f*_{turb} in Equation 55.

10. ^This limit of course is for the non resonant regime. In case of resonant leptogenesis the masses of the sterile neutrinos can be very low [354].

11. ^In cases where a charged next-to-lightest SUSY particle (NSLP) is stable at detector timescales, no clear MET signal is produced, since the NLSP will decay to the LSP outside the detector.

13. ^www.ft.uam.es/personal/rruiz/superbayes

17. ^See also [590] for constraints on couplings to light quarks.

## References

1. Glashow S. Partial symmetries of weak interactions. *Nucl.Phys.* (1961) **22**:579–88. doi: 10.1016/0029-5582(61)90469-2

4. Aad G, Abajyan T, Abbott B, Abdallah J, Khalek SA, Abdelalim AA, et al. Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC. *Phys Lett.* (2012) **B716**:1–29.

5. Chatrchyan S, Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Aguilo E, et al. Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC. *Phys Lett.* (2012) **B716**:30–61. doi: 10.1016/j.physletb.2012.08.021

7. Englert F, Brout R. Broken symmetry and the mass of gauge vector mesons. *Phys Rev Lett.* (1964) **13**:321–3.

8. Guralnik GS, Hagen CR, Kibble TWB. Global conservation laws and massless particles. *Phys Rev Lett.* (1964) **13**:585–7.

9. Degrassi G, Di Vita S, Elias-Miro J, Espinosa JR, Giudice GF, Isidori G, et al. Higgs mass and vacuum stability in the Standard Model at NNLO. *J High Energy Phys.* (2012) **08**:098.

10. Fukuda Y, Hayakawa T, Ichihara E, Inoue K, Ishihara K, Ishino H, et al. Evidence for oscillation of atmospheric neutrinos. *Phys Rev Lett.* (1998) **81**:1562–7.

11. Ahmad QR, Allen RC, Andersen TC, Anglin JD, Buhler G, Barton JC, et al. Measurement of the rate of ${\nu}_{e}+d\to p+p+{e}^{-}$ interactions produced by ^{8}*B* solar neutrinos at the Sudbury Neutrino Observatory. *Phys Rev Lett.* (2001) **87**:071301. doi: 10.1103/PhysRevLett.87.071301

12. Aaij R, Adeva B, Adinolfi M, Affolder A, Ajaltouni Z, Akar S, et al. Test of lepton universality using *B*^{+}→*K*^{+}ℓ^{+}ℓ^{−} decays. *Phys Rev Lett.* (2014) **113**:151601. doi: 10.1103/PhysRevLett.113.151601

13. Aaij R, Adeva B, Adinolfi M, Affolder A, Ajaltouni Z, Akar S, et al. Measurement of the ratio of branching fractions ${B}({\overline{B}}^{0}\to {D}^{*+}{\tau}^{-}{\overline{\nu}}_{\tau})/{B}({\overline{B}}^{0}\to {D}^{*+}{\mu}^{-}{\overline{\nu}}_{\mu})$. *Phys Rev Lett.* (2015) **115**:111803. doi: 10.1103/PhysRevLett.115.111803

14. Aaij R, Adeva B, Adinolfi M, Ajaltouni Z, Akar S, Albrecht J, et al. Measurement of the ${B}_{s}^{0}\to {\mu}^{+}{\mu}^{-}$ branching fraction and effective lifetime and search for *B*^{0} → μ^{+}μ^{−} decays. *Phys Rev Lett.* (2017) **118**:191801. doi: 10.1103/PhysRevLett.118.191801

17. Mohapatra RN, Pati JC. Left-right gauge symmetry and an isoconjugate model of CP violation. *Phys Rev.* (1975) **D11**:566–71.

18. Fritzsch H, Minkowski P. Unified interaction of leptons and hadrons. *Ann Phys.* (1975) **93**:193–266.

21. Gonzalo Velasco T. *Model Building and Phenomenology in Grand Unified Theories*. University Coll. London (2015). Available online at: http://discovery.ucl.ac.uk/id/eprint/1471160

22. Raby S. Supersymmetric grand unified theories. *Lect Notes Phys.* (2017) **939**:1–308. doi: 10.1007/978-3-319-55255-2

23. Nath P. High energy physics and cosmology at the unification frontier: opportunities and challenges in the coming years. *Int J Mod Phys.* (2018) **A33**:1830017. doi: 10.1142/S0217751X1830017X

24. Candelas P, Horowitz GT, Strominger A, Witten E. Vacuum configurations for superstrings. *Nucl Phys.* (1985) **B258**:46–74.

25. Witten E. Symmetry breaking patterns in superstring models. *Nucl Phys.* (1985) **B258**:75. doi: 10.1016/0550-3213(85)90603-0

27. Branco GC, Ferreira PM, Lavoura L, Rebelo MN, Sher M, Silva JP. Theory and phenomenology of two-Higgs-doublet models. *Phys Rept.* (2012) **516**:1–102.

28. Georgi H, Quinn HR, Weinberg S. Hierarchy of interactions in unified gauge theories. *Phys Rev Lett.* (1974) **33**:451–4.

29. Buras AJ, Ellis JR, Gaillard MK, Nanopoulos DV. Aspects of the grand unification of strong, weak and electromagnetic interactions. *Nucl Phys.* (1978) **B135**:66–92.

30. Chanowitz MS, Ellis JR, Gaillard MK. The price of natural flavor conservation in neutral weak interactions. *Nucl Phys.* (1977) **B128**:506–36.

31. Georgi H, Jarlskog C. A new lepton-quark mass relation in a unified theory. *Phys Lett.* (1979) **86B**:297–300.

33. Masiero A, Nanopoulos DV, Tamvakis K, Yanagida T. Naturally massless higgs doublets in supersymmetric SU(5). *Phys Lett.* (1982) **115B**:380–4.

34. Tanabashi M, Hagiwara K, Hikasa K, Nakamura K, Sumino Y, Takahashi F, et al. Review of particle physics. *Phys Rev.* (2018) **D98**:030001. doi: 10.1088/1674-1137/40/10/100001

35. Grinstein B. A supersymmetric SU(5) gauge theory with no gauge hierarchy problem. *Nucl Phys.* (1982) **B206**:387.

36. Hisano J, Moroi T, Tobe K, Yanagida T. Suppression of proton decay in the missing partner model for supersymmetric SU(5) GUT. *Phys Lett.* (1995) **B342**:138–44.

37. Antusch S, de Medeiros Varzielas I, Maurer V, Sluka C, Spinrath M. Towards predictive flavour models in SUSY SU(5) GUTs with doublet-triplet splitting. *J High Energy Phys.* (2014) **09**:141. doi: 10.1007/JHEP09(2014)141

38. Dorsner I, Fileviez Perez P, Gonzalez Felipe R. Phenomenological and cosmological aspects of a minimal GUT scenario. *Nucl Phys.* (2006) **B747**:312–27. doi: 10.1016/j.nuclphysb.2006.05.006

39. Bajc B, Senjanovic G. Seesaw at LHC. *J High Energy Phys*. (2007) **08**:014. doi: 10.1088/1126-6708/2007/08/014

40. Fileviez Perez P. Renormalizable adjoint SU(5). *Phys Lett.* (2007) **B654**:189–93. doi: 10.1016/j.physletb.2007.07.075

41. Zee A. A theory of lepton number violation, neutrino majorana mass, and oscillation. *Phys Lett.* (1980) **93B**:389.

42. Fileviez Perez P, Murgui C. Renormalizable SU(5) unification. *Phys Rev. D*. (2016) **94**:075014. doi: 10.1103/PhysRevD.94.075014

43. Abe K, Haga Y, Hayato Y, M Ikeda KI, Kameda J, et al. Search for proton decay via *p*→*e*^{+}π^{0} and *p* → μ^{+}π^{0} in 0.31 megaton years exposure of the super-kamiokande water cherenkov detector. *Phys Rev D.* (2017) **95**:012004. doi: 10.1103/PhysRevD.95.012004

44. Dorsner I, Fileviez Perez P. Unification without supersymmetry: neutrino mass, proton decay and light leptoquarks. *Nucl Phys.* (2005) **B723**:53–76. doi: 10.1016/j.nuclphysb.2005.06.016

45. Dorsner I, Fileviez Perez P. Unification versus proton decay in SU(5). *Phys Lett.* (2006) **B642**:248–52. doi: 10.1016/j.nuclphysb.2005.06.016

46. Fornal B, Grinstein B. Grand unified theory with a stable proton. (2018) *Int J Mod Phys.* **A33**:1844013. doi: 10.1142/S0217751X1844013X

47. Fileviez Pérez P, Gross A, Murgui C. Seesaw scale, unification, and proton decay. *Phys Rev D.* (2018) **98**:035032. doi: 10.1103/PhysRevD.98.035032

48. Weinberg S. Implications of dynamical symmetry breaking. *Phys Rev.* (1976) **D13**:974–96. doi: 10.1103/PhysRevD.13.974

50. Nanopoulos DV, Tamvakis K. SUSY GUTS: 4 - GUTS: 3. *Phys Lett.* (1982) **113B**:151–8. doi: 10.1016/0370-2693(82)90413-0

51. Dorsner I, Fileviez Perez P, Rodrigo G. On unification and nucleon decay in supersymmetric grand unified theories based on SU(5). *Phys Lett.* (2007) **B649**:197–205. doi: 10.1016/j.physletb.2007.03.062

52. 't Hooft G. Magnetic monopoles in unified gauge theories. *Nucl Phys.* (1974) **B79**:276–84. doi: 10.1016/0550-3213(74)90486-6

53. Witten E. Mass hierarchies in supersymmetric theories. *Phys Lett.* (1981) **105B**:267. doi: 10.1016/0370-2693(81)90885-6

54. Dimopoulos S, Georgi H. Solution of the Gauge Hierarchy Problem. *Phys Lett.* (1982) **117B**:287–90. doi: 10.1016/0370-2693(82)90720-1

55. De Rujula A, Georgi H, Glashow SL. FLAVOR GONIOMETRY BY PROTON DECAY. *Phys Rev Lett.* (1980) **45**:413. doi: 10.1103/PhysRevLett.45.413

56. Barr SM. A new symmetry breaking pattern for SO(10) and proton decay. *Phys Lett.* (1982) **B112**:219. doi: 10.1016/0370-2693(82)90966-2

57. Derendinger JP, Kim JE, Nanopoulos DV. Anti-SU(5). *Phys Lett.* (1984) **B139**:170. doi: 10.1016/0370-2693(84)91238-3

58. Antoniadis I, Ellis JR, Hagelin JS, Nanopoulos DV. Supersymmetric flipped SU(5) revitalized. *Phys Lett.* (1987) **B194**:231. doi: 10.1016/0370-2693(87)90533-8

59. Ellis JR, Hagelin JS, Kelley S, Nanopoulos DV. Aspects of the flipped unification of strong, weak and electromagnetic interactions. *Nucl Phys.* (1988) **B311**:1. doi: 10.1016/0550-3213(88)90141-1

60. Campbell BA, Ellis JR, Hagelin JS, Nanopoulos DV, Ticciati R. Flipped SU(5) from manifold compactification of the ten-dimensional heterotic string. *Phys Lett.* (1987) **B198**:200–4. doi: 10.1016/0370-2693(87)91496-1

61. Antoniadis I, Ellis JR, Hagelin JS, Nanopoulos DV. GUT model building with fermionic four-dimensional strings. *Phys Lett.* (1988) **B205**:459–65. doi: 10.1016/0370-2693(88)90978-1

62. Senjanovic G, Mohapatra RN. Exact left-right symmetry and spontaneous violation of parity. *Phys Rev*. (1975) **D12**:1502. doi: 10.1103/PhysRevD.12.1502

63. Mohapatra RN, Pati JC. A natural left-right symmetry. *Phys Rev*. (1975) **D11**:2558. doi: 10.1103/PhysRevD.11.2558

64. Smirnov AD. The minimal quark-lepton symmetry model and the limit on Z-prime mass. *Phys Lett*. (1995) **B346**:297–302. doi: 10.1016/0370-2693(95)00015-D

65. Fileviez Perez P, Wise MB. Low scale quark-lepton unification. *Phys Rev D*. (2013) **88**:057703. doi: 10.1103/PhysRevD.88.057703

66. Mohapatra RN, Senjanovic G. Neutrino mass and spontaneous parity violation. *Phys Rev Lett*. (1980) **44**:912. doi: 10.1103/PhysRevLett.44.912

67. Mohapatra RN, Senjanovic G. Neutrino masses and mixings in gauge models with spontaneous parity violation. *Phys Rev*. (1981) **D23**:165. doi: 10.1103/PhysRevD.23.165

68. Mohapatra RN, Marshak RE. Local B-L symmetry of electroweak interactions, majorana neutrinos and neutron oscillations. *Phys Rev Lett*. (1980) **44**:1316–9. doi: 10.1103/PhysRevLett.44.1316

69. Lindner M, Weiser M. Gauge coupling unification in left-right symmetric models. *Phys Lett*. (1996) **B383**:405–14. doi: 10.1016/0370-2693(96)00775-7

70. Aulakh CS, Bajc B, Melfo A, Rasin A, Senjanovic G. SO(10) theory of R-parity and neutrino mass. *Nucl Phys*. (2001) **B597**:89–109. doi: 10.1016/S0550-3213(00)00721-5

71. Senjanovic G. Spontaneous breakdown of parity in a class of gauge theories. *Nucl Phys*. (1979) **B153**:334–64. doi: 10.1016/0550-3213(79)90604-7

72. Antoniadis I, Leontaris GK. A SUPERSYMMETRIC SU(4) x O(4) MODEL. *Phys Lett*. (1989) **B216**:333–5. doi: 10.1016/0370-2693(89)91125-8

73. Kuchimanchi R, Mohapatra RN. No parity violation without R-parity violation. *Phys Rev*. (1993) **D48**:4352–60. doi: 10.1103/PhysRevD.48.4352

74. Cvetic M. Effects of the parity-odd singlet on the *N* = 1 supergravity theory within the left-right symmetric model. *Phys Lett*. (1985) **164B**:55–60. doi: 10.1016/0370-2693(85)90030-9

75. Aulakh CS, Benakli K, Senjanovic G. Reconciling supersymmetry and left-right symmetry. *Phys Rev Lett*. (1997) **79**:2188–91. doi: 10.1103/PhysRevLett.79.2188

76. Aulakh CS, Melfo A, Senjanovic G. Minimal supersymmetric left-right model. *Phys Rev*. (1998) **D57**:4174–8. doi: 10.1103/PhysRevD.57.4174

77. Aulakh CS, Melfo A, Rasin A, Senjanovic G. Supersymmetry and large scale left-right symmetry. *Phys Rev*. (1998) **D58**:115007. doi: 10.1103/PhysRevD.58.115007

78. Georgi H, Glashow SL. Gauge theories without anomalies. *Phys Rev.* (1972) **D6**:429. doi: 10.1103/PhysRevD.6.429

79. Adler SL. Axial vector vertex in spinor electrodynamics. *Phys Rev*. (1969) **177**:2426–38. doi: 10.1103/PhysRev.177.2426

80. Mohapatra RN, Sakita B. SO(2n) grand unification in an SU(N) basis. *Phys Rev*. (1980) **D21**:1062. doi: 10.1103/PhysRevD.21.1062

81. Bertolini S, Di Luzio L, Malinsky M. The quantum vacuum of the minimal SO(10) GUT. *J Phys Conf Ser*. (2010) **259**:012098. doi: 10.1088/1742-6596/259/1/012098

82. Bertolini S, Di Luzio L, Malinsky M. Seesaw scale in the minimal renormalizable SO(10) grand unification. *Phys Rev*. (2012) **D85**:095014. doi: 10.1103/PhysRevD.85.095014

83. Bertolini S, Di Luzio L, Malinsky M. Light color octet scalars in the minimal SO(10) grand unification. *Phys Rev*. (2013) **D87**:085020. doi: 10.1103/PhysRevD.87.085020

84. Kolesova H, Malinsky M. Proton lifetime in the minimal SO(10) GUT and its implications for the LHC. *Phys Rev*. (2014) **D90**:115001. doi: 10.1103/PhysRevD.90.115001

85. Graf L, Malinsky M, Mede T, Susic V. One-loop pseudo-Goldstone masses in the minimal *SO*(10) Higgs model. *Phys Rev D*. (2017) **95**:075007. doi: 10.1103/PhysRevD.95.075007

86. Bertolini S, Schwetz T, Malinsky M. Fermion masses and mixings in SO(10) models and the neutrino challenge to SUSY GUTs. *Phys Rev*. (2006) **D73**:115012. doi: 10.1103/PhysRevD.73.115012

87. Aulakh CS, Garg SK. MSGUT : From bloom to doom. *Nucl Phys*. (2006) **B757**:47–78. doi: 10.1016/j.nuclphysb.2006.07.030

88. Aulakh CS, Girdhar A. SO(10) a la Pati-Salam. *Int J Mod Phys*. (2005) **A20**:865–94. doi: 10.1142/S0217751X0502001X

89. Anastaze G, Derendinger JP, Buccella F. INTERMEDIATE SYMMETRIES IN THE SO(10) MODEL WITH (16+16) + 45 HIGGSES. *Z Phys*. (1983) **C20**:269–73. doi: 10.1007/BF01574861

90. Clark TE, Kuo TK, Nakagawa N. A SO(10) SUPERSYMMETRIC GRAND UNIFIED THEORY. *Phys Lett*. (1982) **115B**:26–8. doi: 10.1016/0370-2693(82)90507-X

91. Aulakh CS, Mohapatra RN. Implications of supersymmetric SO(10) grand unification. *Phys Rev*. (1983) **D28**:217. doi: 10.1103/PhysRevD.28.217

93. Reig M, Valle JWF, Vaquera-Araujo CA, Wilczek F. A Model of Comprehensive Unification. *Phys Lett*. (2017) **B774**:667–70. doi: 10.1016/j.physletb.2017.10.038

94. Langacker P, Wang J. U(1)-prime symmetry breaking in supersymmetric E(6) models. *Phys Rev*. (1998) **D58**:115010. doi: 10.1103/PhysRevD.58.115010

95. Gursey F, Ramond P, Sikivie P. A universal gauge theory model based on E6. *Phys Lett*. (1976) **60B**:177–80. doi: 10.1016/0370-2693(76)90417-2

96. Achiman Y, Stech B. Quark lepton symmetry and mass scales in an E6 unified gauge model. *Phys Lett*. (1978) **77B**:389–93. doi: 10.1016/0370-2693(78)90584-1

97. Shafi Q. E(6) as a unifying gauge symmetry. *Phys Lett*. (1978) **79B**:301–3. doi: 10.1016/0370-2693(78)90248-4

98. Mohapatra RN, Valle JWF. Neutrino mass and baryon number nonconservation in superstring models. *Phys Rev*. (1986) **D34**:1642. doi: 10.1103/PhysRevD.34.1642

99. Buccella F, Miele G. SO(10) from supersymmetric *E*(6). *Phys Lett*. (1987) **B189**:115–7. doi: 10.1016/0370-2693(87)91280-9

100. Hosotani Y. Dynamical mass generation by compact extra dimensions. *Phys Lett*. (1983) **126B**:309–13. doi: 10.1016/0370-2693(83)90170-3

101. Hewett JL, Rizzo TG. Low-energy phenomenology of superstring inspired E(6) models. *Phys Rept*. (1989) **183**:193. doi: 10.1016/0370-1573(89)90071-9

102. London D, Rosner JL. Extra gauge bosons in E(6). *Phys Rev*. (1986) **D34**:1530. doi: 10.1103/PhysRevD.34.1530

103. Dimopoulos S, Raby S, Wilczek F. Supersymmetry and the scale of unification. *Phys Rev*. (1981) **D24**:1681–3. doi: 10.1103/PhysRevD.24.1681

104. Witten E. Dynamical breaking of supersymmetry. *Nucl Phys*. (1981) **B188**:513. doi: 10.1016/0550-3213(81)90006-7

105. Ibanez LE, Ross GG. SU(2)-L x U(1) symmetry breaking as a radiative effect of supersymmetry breaking in guts. *Phys Lett*. (1982) **110B**:215–20. doi: 10.1016/0370-2693(82)91239-4

106. Ellis JR, Nanopoulos DV, Tamvakis K. Grand unification in simple supergravity. *Phys Lett*. (1983) **121B**:123–9. doi: 10.1016/0370-2693(83)90900-0

107. Farrar GR, Fayet P. Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry. *Phys Lett*. (1978) **76B**:575–9. doi: 10.1016/0370-2693(78)90858-4

108. Ellis JR, Hagelin JS, Nanopoulos DV, Olive KA, Srednicki M. Supersymmetric relics from the big bang. *Nucl Phys*. (1984) **B238**:453–76. doi: 10.1016/0550-3213(84)90461-9

109. Jungman G, Kamionkowski M, Griest K. Supersymmetric dark matter. *Phys Rept*. (1996) **267**:195–373. doi: 10.1016/0370-1573(95)00058-5

110. Ellis JR, Kelley S, Nanopoulos DV. Probing the desert using gauge coupling unification. *Phys Lett*. (1991) **B260**:131–7. doi: 10.1016/0370-2693(91)90980-5

111. Martin SP, Vaughn MT. Two loop renormalization group equations for soft supersymmetry breaking couplings. *Phys Rev*. (1994) **D50**:2282. doi: 10.1103/PhysRevD.50.2282

112. Hall LJ. Grand unification of effective gauge theories. *Nucl Phys*. (1981) **B178**:75–124. doi: 10.1016/0550-3213(81)90498-3

113. Weinberg S. Effective gauge theories. *Phys Lett*. (1980) **91B**:51–5. doi: 10.1016/0370-2693(80)90660-7

114. Giudice GF, Romanino A. Split supersymmetry. *Nucl Phys*. (2004) **B699**:65–89. doi: 10.1016/j.nuclphysb.2004.08.001

115. Arkani-Hamed N, Dimopoulos S. Supersymmetric unification without low energy supersymmetry and signatures for fine-tuning at the LHC. *J High Energy Phys*. (2005) **06**:073. doi: 10.1088/1126-6708/2005/06/073

116. Athron P, Balázs C, Buckley A, Cornell JM, Danninger M, Farmer B, et al. Combined collider constraints on neutralinos and charginos. (2018). doi: 10.1140/epjc/s10052-019-6837-x

117. Ellis SAR, Wells JD. Visualizing gauge unification with high-scale thresholds. *Phys Rev*. (2015) **D91**:075016. doi: 10.1103/PhysRevD.91.075016

118. Pontecorvo B. Neutrino experiments and the problem of conservation of leptonic charge. *Sov Phys JETP*. (1968) **26**:984–8. doi: 10.1016/B978-0-12-395657-6.50020-3

119. Schechter J, Valle JWF. Neutrino masses in SU(2) x U(1) theories. *Phys Rev*. (1980) **D22**:2227. doi: 10.1103/PhysRevD.22.2227

120. Ma E. Pathways to naturally small neutrino masses. *Phys Rev Lett*. (1998) **81**:1171–4. doi: 10.1103/PhysRevLett.81.1171

121. Minkowski P. μ → *eγ* at a rate of one out of 10^{9} muon decays? *Phys Lett*. (1977) **67B**:421–8. doi: 10.1016/0370-2693(77)90435-X

123. Gell-Mann M, Ramond P, Slansky R. Complex spinors and unified theories. *Conf Proc*. (1979) **C790927**:315–21.

124. Magg M, Wetterich C. Neutrino mass problem and gauge hierarchy. *Phys Lett*. (1980) **94B**:61–4. doi: 10.1016/0370-2693(80)90825-4

125. Lazarides G, Shafi Q, Wetterich C. Proton lifetime and fermion masses in an SO(10) model. *Nucl Phys*. (1981) **B181**:287–300. doi: 10.1016/0550-3213(81)90354-0

126. Goh HS, Mohapatra RN, Nasri S. SO(10) symmetry breaking and type II seesaw. *Phys Rev*. (2004) **D70**:075022. doi: 10.1103/PhysRevD.70.075022

127. Mohapatra RN, Parida MK. Type II seesaw dominance in non-supersymmetric and split susy SO(10) and proton life time. *Phys Rev*. (2011) **D84**:095021. doi: 10.1103/PhysRevD.84.095021

128. Foot R, Lew H, He XG, Joshi GC. Seesaw neutrino masses induced by a triplet of leptons. *Z Phys*. (1989) **C44**:441. doi: 10.1007/BF01415558

129. Fileviez Perez P. Supersymmetric adjoint SU(5). *Phys Rev*. (2007) **D76**:071701. doi: 10.1103/PhysRevD.76.071701

130. Gonzalez-Garcia MC, Valle JWF. Fast decaying neutrinos and observable flavor violation in a new class of majoron models. *Phys Lett*. (1989) **B216**:360–6. doi: 10.1016/0370-2693(89)91131-3

131. 't Hooft G. Naturalness, chiral symmetry, and spontaneous chiral symmetry breaking. *NATO Sci Ser B*. (1980) **59**:135–57. doi: 10.1007/978-1-4684-7571-5_9

132. Malinsky M, Romao JC, Valle JWF. Novel supersymmetric SO(10) seesaw mechanism. *Phys Rev Lett*. (2005) **95**:161801. doi: 10.1103/PhysRevLett.95.161801

133. Hirsch M, Valle JWF. Supersymmetric origin of neutrino mass. New *J Phys*. (2004) **6**:76. doi: 10.1088/1367-2630/6/1/076

134. Babu KS. Model of ‘Calculable’ majorana neutrino masses. *Phys Lett*. (1988) **B203**:132–6. doi: 10.1016/0370-2693(88)91584-5

135. Witten E. Neutrino masses in the minimal O(10) theory. *Phys Lett*. (1980) **91B**:81–4. doi: 10.1016/0370-2693(80)90666-8

136. Leontaris GK, Vergados JD. The Seesaw mechanism in string models. *Phys Lett*. (1991) **B258**:111–7. doi: 10.1016/0370-2693(91)91217-J

137. Arbeláez Rodríguez C, Kolešová H, Malinský M. Witten's mechanism in the flipped SU(5) unification. *Phys Rev D*. (2014) **89**:055003. doi: 10.1103/PhysRevD.89.055003

138. Harries D, Malinský M, Zdráhal M. Witten's loop in the minimal flipped *SU*(5) unification revisited. *Phys Rev D*. (2018) **98**:095015. doi: 10.1103/PhysRevD.98.095015

139. Giveon A, Hall LJ, Sarid U. SU(5) unification revisited. *Phys Lett*. (1991) **B271**:138–44. doi: 10.1016/0370-2693(91)91289-8

140. Antusch S, King SF, Spinrath M. GUT predictions for quark-lepton Yukawa coupling ratios with messenger masses from non-singlets. *Phys Rev D*. (2014) **89**:055027. doi: 10.1103/PhysRevD.89.055027

141. Ellis JR, Gaillard MK. Fermion masses and higgs representations in SU(5). *Phys Lett*. (1979) **88B**:315–9. doi: 10.1016/0370-2693(79)90476-3

142. Kelley S, Lopez JL, Nanopoulos DV. Yukawa unification. *Phys Lett*. (1992) **B274**:387–92. doi: 10.1016/0370-2693(92)92003-Y

143. Ananthanarayan B, Lazarides G, Shafi Q. Top mass prediction from supersymmetric guts. *Phys Rev*. (1991) **D44**:1613–5. doi: 10.1103/PhysRevD.44.1613

144. Antusch S, Spinrath M. Quark and lepton masses at the GUT scale including SUSY threshold corrections. *Phys Rev*. (2008) **D78**:075020. doi: 10.1103/PhysRevD.78.075020

145. Antusch S, Spinrath M. New GUT predictions for quark and lepton mass ratios confronted with phenomenology. *Phys Rev*. (2009) **D79**:095004. doi: 10.1103/PhysRevD.79.095004

146. Babu KS, Kolda CF. Signatures of supersymmetry and Yukawa unification in Higgs decays. *Phys Lett*. (1999) **B451**:77–85. doi: 10.1016/S0370-2693(99)00204-X

147. Olechowski M, Pokorski S. Electroweak symmetry breaking with nonuniversal scalar soft terms and large tan beta solutions. *Phys Lett*. (1995) **B344**:201–10. doi: 10.1016/0370-2693(94)01571-S

148. Blazek T, Dermisek R, Raby S. Yukawa unification in SO(10). *Phys Rev*. (2002) **D65**:115004. doi: 10.1103/PhysRevD.65.115004

149. Hall LJ, Rattazzi R, Sarid U. The Top quark mass in supersymmetric SO(10) unification. *Phys Rev*. (1994) **D50**:7048–65. doi: 10.1103/PhysRevD.50.7048

150. King SF, Oliveira M. Yukawa unification as a window into the soft supersymmetry breaking Lagrangian. *Phys Rev*. (2001) **D63**:015010. doi: 10.1103/PhysRevD.63.015010

151. Baer H, Raza S, Shafi Q. A heavier gluino from t-b-t Yukawa-unified SUSY. *Phys Lett*. (2012) **B712**:250–4. doi: 10.1016/j.physletb.2012.04.075

152. Langacker P, Polonsky N. The Bottom mass prediction in supersymmetric grand unification: Uncertainties and constraints. *Phys Rev*. (1994) **D49**:1454–67. doi: 10.1103/PhysRevD.49.1454

153. Langacker P, Polonsky N. Implications of Yukawa unification for the Higgs sector in supersymmetric grand unified models. *Phys Rev*. (1994) **D50**:2199–217. doi: 10.1103/PhysRevD.50.2199

154. Baer H, Gogoladze I, Mustafayev A, Raza S, Shafi Q. Sparticle mass spectra from SU(5) SUSY GUT models with *b*−τ Yukawa coupling unification. *J High Energy Phys*. (2012) **03**:047. doi: 10.1007/JHEP03(2012)047

155. Gomez ME, Lazarides G, Pallis C. Yukawa quasi-unification. *Nucl Phys*. (2002) **B638**:165–85. doi: 10.1016/S0550-3213(02)00483-2

156. Gogoladze I, Khalid R, Raza S, Shafi Q. *t*−*b*−τ Yukawa unification for μ < 0 with a sub-TeV sparticle spectrum. *J High Energy Phys*. (2010) **12**:055. doi: 10.1007/JHEP12(2010)055

157. Bajc B, Senjanovic G, Vissani F. b - tau unification and large atmospheric mixing: a Case for noncanonical seesaw. *Phys Rev Lett*. (2003) **90**:051802. doi: 10.1103/PhysRevLett.90.051802

158. Anderson G, Raby S, Dimopoulos S, Hall LJ, Starkman GD. A Systematic SO(10) operator analysis for fermion masses. *Phys Rev*. (1994) **D49**:3660–90. doi: 10.1103/PhysRevD.49.3660

159. Altarelli G, Feruglio F. Discrete flavor symmetries and models of neutrino mixing. *Rev Mod Phys*. (2010) **82**:2701–29. doi: 10.1103/RevModPhys.82.2701

160. Ishimori H, Kobayashi T, Ohki H, Shimizu Y, Okada H, Tanimoto M. Non-abelian discrete symmetries in particle physics. *Prog Theor Phys Suppl*. (2010) **183**:1–163. doi: 10.1143/PTPS.183.1

161. Grimus W, Ludl PO. Finite flavour groups of fermions. *J Phys*. (2012) **A45**:233001. doi: 10.1088/1751-8113/45/23/233001

162. King SF, Luhn C, Stuart AJ. A grand Delta(96) x SU(5) flavour model. *Nucl Phys*. (2013) **B867**:203–35. doi: 10.1016/j.nuclphysb.2012.09.021

163. de Anda FJ, King SF. *SU*(3) × *SO*(10) in 6d. *J High Energy Phys*. (2018) **10**:128. doi: 10.1007/JHEP10(2018)128

164. Arbeláez C, Hirsch M, Malinský M, Romão JC. LHC-scale left-right symmetry and unification. *Phys Rev*. (2014) **D89**:035002. doi: 10.1103/PhysRevD.89.035002

165. Deppisch FF, Gonzalo TE, Graf L. Surveying the SO(10) model landscape: the left-right symmetric case. *Phys Rev*. (2017) **D96**:055003. doi: 10.1103/PhysRevD.96.055003

166. Senjanovic G. Seesaw at LHC through left-right symmetry. *Int J Mod Phys*. (2011) **A26**:1469–91. doi: 10.1142/S0217751X1105302X

167. Nemevsek M, Nesti F, Senjanovic G, Zhang Y. First limits on left-right symmetry scale from LHC data. *Phys Rev*. (2011) **D83**:115014. doi: 10.1103/PhysRevD.83.115014

168. Han T, Lewis I, Ruiz R, Si Zg. Lepton number violation and *W*′ chiral couplings at the LHC. *Phys Rev*. (2013) **D87**:035011. doi: 10.1103/PhysRevD.87.039906

169. Das SP, Deppisch FF, Kittel O, Valle JWF. Heavy neutrinos and lepton flavour violation in left-right symmetric models at the LHC. *Phys Rev*. (2012) **D86**:055006. doi: 10.1103/PhysRevD.86.055006

170. Bambhaniya G, Chakrabortty J, Gluza J, Kordiaczynska M, Szafron R. Left-right symmetry and the charged higgs bosons at the LHC. *J High Energy Phys*. (2014) **05**:033. doi: 10.1007/JHEP05(2014)033

171. Chen CY, Dev PSB, Mohapatra RN. Probing heavy-light neutrino mixing in left-right seesaw models at the LHC. *Phys Rev D*. (2013) **88**:033014. doi: 10.1103/PhysRevD.88.033014

172. Lee CH, Bhupal Dev PS, Mohapatra RN. Natural TeV-scale left-right seesaw mechanism for neutrinos and experimental tests. *Phys Rev D*. (2013) **88**:093010. doi: 10.1103/PhysRevD.88.093010

173. Deppisch FF, Gonzalo TE, Patra S, Sahu N, Sarkar U. Signal of right-handed charged gauge bosons at the LHC? *Phys Rev D*. (2014) **90**:053014. doi: 10.1103/PhysRevD.90.053014

174. Deppisch FF, Gonzalo TE, Patra S, Sahu N, Sarkar U. Double beta decay, lepton flavor violation, and collider signatures of left-right symmetric models with spontaneous *D*-parity breaking. *Phys Rev D*. (2015) **91**:015018. doi: 10.1103/PhysRevD.91.015018

175. Deppisch FF, Graf L, Kulkarni S, Patra S, Rodejohann W, Sahu N, et al. Reconciling the 2 TeV excesses at the LHC in a linear seesaw left-right model. *Phys Rev D*. (2016) **93**:013011. doi: 10.1103/PhysRevD.93.013011

176. Bhupal Dev PS, Mohapatra RN. Unified explanation of the *eejj*, diboson and dijet resonances at the LHC. *Phys Rev Lett*. (2015) **115**:181803. doi: 10.1103/PhysRevLett.115.181803

177. Mitra M, Ruiz R, Scott DJ, Spannowsky M. Neutrino jets from high-mass *W*_{R} gauge bosons in TeV-Scale left-right symmetric models. *Phys Rev D*. (2016) **94**:095016. doi: 10.1103/PhysRevD.94.095016

178. Helo JC, Li H, Neill NA, Ramsey-Musolf M, Vasquez JC. Probing neutrino Dirac Mass in Left-Right Symmetric Models at the LHC and Next Generation Colliders. *Phys Rev D.* (2018) **99**:055042. doi: 10.1103/PhysRevD.99.055042

179. Mohapatra RN. Mechanism for understanding small neutrino mass in superstring theories. *Phys Rev Lett*. (1986) **56**:561–3. doi: 10.1103/PhysRevLett.56.561

180. Dev PSB, Mohapatra RN. TeV scale inverse seesaw in SO(10) and leptonic non-unitarity effects. *Phys Rev*. (2010) **D81**:013001. doi: 10.1103/PhysRevD.81.013001

181. Brdar V, Smirnov AY. Low scale left-right symmetry and naturally small neutrino mass. *J High Energy Phys.* (2018) **1902**:045. doi: 10.1007/JHEP02(2019)045

182. Akhmedov EK, Lindner M, Schnapka E, Valle JWF. Left-right symmetry breaking in NJL approach. *Phys Lett*. (1996) **B368**:270–80. doi: 10.1016/0370-2693(95)01504-3

183. Akhmedov EK, Lindner M, Schnapka E, Valle JWF. Dynamical left-right symmetry breaking. *Phys Rev*. (1996) **D53**:2752–80. doi: 10.1103/PhysRevD.53.2752

184. Fileviez Perez P. Type III Seesaw and Left-Right Symmetry. *J High Energy Phys*. (2009) **03**:142. doi: 10.1088/1126-6708/2009/03/142

185. Duerr M, Fileviez Perez P, Lindner M. Left-right symmetric theory with light sterile neutrinos. *Phys Rev D*. (2013) **88**:051701. doi: 10.1103/PhysRevD.88.051701

186. Fileviez Perez P, Murgui C, Ohmer S. Simple left-right theory: lepton number violation at the LHC. *Phys Rev D*. (2016) **94**:051701. doi: 10.1103/PhysRevD.94.051701

187. Babu KS, Mohapatra RN. Predictive neutrino spectrum in minimal SO(10) grand unification. *Phys Rev Lett*. (1993) **70**:2845–8. doi: 10.1103/PhysRevLett.70.2845

188. Albright CH, Barr SM. Construction of a minimal Higgs SO(10) SUSY GUT model. *Phys Rev*. (2000) **D62**:093008. doi: 10.1103/PhysRevD.62.093008

189. Fukuyama T, Okada N. Neutrino oscillation data versus minimal supersymmetric SO(10) model. *J High Energy Phys*. (2002) **11**:011. doi: 10.1088/1126-6708/2002/11/011

190. Goh HS, Mohapatra RN, Ng SP. Minimal SUSY SO(10) model and predictions for neutrino mixings and leptonic CP violation. *Phys Rev*. (2003) **D68**:115008. doi: 10.1103/PhysRevD.68.115008

191. Goh HS, Mohapatra RN, Ng SP. Minimal SUSY SO(10), b tau unification and large neutrino mixings. *Phys Lett*. (2003) **B570**:215–21. doi: 10.1016/j.physletb.2003.08.011

192. Bertolini S, Frigerio M, Malinsky M. Fermion masses in SUSY SO(10) with type II seesaw: a Non-minimal predictive scenario. *Phys Rev*. (2004) **D70**:095002. doi: 10.1103/PhysRevD.70.095002

193. Matsuda K, Koide Y, Fukuyama T. Can the SO(10) model with two Higgs doublets reproduce the observed fermion masses? *Phys Rev*. (2001) **D64**:053015. doi: 10.1103/PhysRevD.64.053015

194. Matsuda K, Koide Y, Fukuyama T, Nishiura H. How far can the SO(10) two Higgs model describe the observed neutrino masses and mixings? *Phys Rev*. (2002) **D65**:033008. doi: 10.1103/PhysRevD.65.033008

195. Babu KS, Macesanu C. Neutrino masses and mixings in a minimal SO(10) model. *Phys Rev*. (2005) **D72**:115003. doi: 10.1103/PhysRevD.72.115003

196. Dueck A, Rodejohann W. Fits to SO(10) Grand Unified Models. *J High Energy Phys*. (2013) **09**:024. doi: 10.1007/JHEP09(2013)024

197. Deppisch T, Schacht S, Spinrath M. Confronting SUSY SO(10) with updated lattice and neutrino data. *J High Energy Phys*. (2019) **01**:005. doi: 10.1007/JHEP01(2019)005

198. Aulakh CS, Bajc B, Melfo A, Senjanovic G, Vissani F. The Minimal supersymmetric grand unified theory. *Phys Lett*. (2004) **B588**:196–202. doi: 10.1016/j.physletb.2004.03.031

199. Aulakh CS, Girdhar A. SO(10) MSGUT: Spectra, couplings and threshold effects. *Nucl Phys*. (2005) **B711**:275–313. doi: 10.1016/j.nuclphysb.2005.01.008

200. Bajc B, Melfo A, Senjanovic G, Vissani F. The Minimal supersymmetric grand unified theory. 1. Symmetry breaking and the particle spectrum. *Phys Rev*. (2004) **D70**:035007. doi: 10.1103/PhysRevD.70.035007

201. Bajc B, Melfo A, Senjanovic G, Vissani F. Fermion mass relations in a supersymmetric SO(10) theory. *Phys Lett*. (2006) **B634**:272–7. doi: 10.1063/1.2149692

202. Bajc B, Dorsner I, Nemevsek M. Minimal SO(10) splits supersymmetry. *J High Energy Phys*. (2008) **11**:007. doi: 10.1088/1126-6708/2008/11/007

203. De Romeri V, Hirsch M, Malinsky M. Soft masses in SUSY SO(10) GUTs with low intermediate scales. *Phys Rev*. (2011) **D84**:053012. doi: 10.1103/PhysRevD.84.053012

204. Deppisch FF, Desai N, Gonzalo TE. Compressed and Split Spectra in Minimal SUSY SO(10). *Frontin Phys*. (2014) **2**:27. doi: 10.3389/fphy.2014.00027

205. Poh Z, Raby S. Yukawa Unification in an SO(10) SUSY GUT: SUSY on the Edge. *Phys Rev D*. (2015) **92**:015017. doi: 10.1103/PhysRevD.92.015017

206. Ellis SAR, Gherghetta T, Kaneta K, Olive KA. New Weak-Scale Physics from SO(10) with High-Scale Supersymmetry. *Phys Rev D*. (2018) **98**:055009. doi: 10.1103/PhysRevD.98.055009

207. Dutta B, Mimura Y, Mohapatra RN. Suppressing proton decay in the minimal SO(10) model. *Phys Rev Lett*. (2005) **94**:091804. doi: 10.1103/PhysRevLett.94.091804

208. Aulakh CS, Garg SK. The new minimal supersymmetric GUT: spectra, RG analysis and fermion fits. *Nucl Phys*. (2012) **B857**:101–42. doi: 10.1016/j.nuclphysb.2011.12.003

209. Fukuyama T, Ilakovac A, Kikuchi T, Meljanac S, Okada N. SO(10) group theory for the unified model building. *J Math Phys*. (2005) **46**:033505. doi: 10.1063/1.1847709

210. Babu KS, Bajc B, Saad S. Resurrecting minimal yukawa sector of SUSY SO(10). *J High Energy Phys*. (2018) **10**:135. doi: 10.1007/JHEP10(2018)135

211. Bhupal Dev PS, Mohapatra RN. Electroweak symmetry breaking and proton decay in SO(10) SUSY-GUT with TeV W(R). *Phys Rev D*. (2010) **82**:035014. doi: 10.1103/PhysRevD.82.035014

212. Anandakrishnan A, Bryant BC, Raby S, Wingerter A. LHC phenomenology of SO(10) models with Yukawa Unification. *Phys Rev D*. (2013) **88**:075002. doi: 10.1103/PhysRevD.88.075002

213. Baer H, Kraml S, Sekmen S, Summy H. Dark matter allowed scenarios for Yukawa-unified SO(10) SUSY GUTs. *J High Energy Phys*. (2008) **03**:056. doi: 10.1088/1126-6708/2008/03/056

214. Goh HS, Mohapatra RN, Nasri S, Ng SP. Proton decay in a minimal SUSY SO(10) model for neutrino mixings. *Phys Lett*. (2004) **B587**:105–16. doi: 10.1016/j.physletb.2004.02.063

215. Mohapatra RN, Severson M. Leptonic *CP* violation and proton decay in SUSY SO(10). *J High Energy Phys*. (2018) **09**:119. doi: 10.1007/JHEP09(2018)119

216. Masiero A, Vempati SK, Vives O. Seesaw and lepton flavor violation in SUSY SO(10). *Nucl Phys*. (2003) **B649**:189–204. doi: 10.1016/S0550-3213(02)01031-3

217. Dermisek R, Raby S, Roszkowski L, Ruiz de Austri R. Dark matter and B(s) —> mu+ mu- with minimal SO(10) soft SUSY breaking II. *J High Energy Phys*. (2005) **09**:029. doi: 10.1088/1126-6708/2005/09/029

218. Babu KS, He XG, Pakvasa S. Neutrino masses and proton decay modes in SU(3) X SU(3) X SU(3) trinification. *Phys Rev*. (1986) **D33**:763. doi: 10.1103/PhysRevD.33.763

219. Sayre J, Wiesenfeldt S, Willenbrock S. Minimal trinification. *Phys Rev*. (2006) **D73**:035013. doi: 10.1103/PhysRevD.73.035013

220. Stech B. The mass of the Higgs boson in the trinification subgroup of E6. *Phys Rev*. (2012) **D86**:055003. doi: 10.1103/PhysRevD.86.055003

221. Willenbrock S. Triplicated trinification. *Phys Lett*. (2003) **B561**:130–4. doi: 10.1016/S0370-2693(03)00419-2

223. Cauet C, Pas H, Wiesenfeldt S. Trinification, the hierarchy problem and inverse seesaw neutrino masses. *Phys Rev*. (2011) **D83**:093008. doi: 10.1103/PhysRevD.83.093008

224. Dvali GR, Shafi Q. Gauge hierarchy in SU(3)(C) x SU(3)(L) x SU(3)-R and low-energy implications. *Phys Lett*. (1994) **B326**:258–63. doi: 10.1016/0370-2693(94)91319-6

226. Nath P, Arnowitt RL. Symmetry breaking in three generation calabi-yau manifolds. *Phys Rev*. (1989) **D39**:2006.

227. Maekawa N, Shafi Q. Supersymmetric SU^{(3)}3 unification with anomalous U(1)(a) gauge symmetry. *Prog Theor Phys*. (2003) **109**:279–93. doi: 10.1143/PTP.109.279

228. Hetzel J, Stech B. Low-energy phenomenology of trinification: an effective left-right-symmetric model. *Phys Rev*. (2015) **D91**:055026. doi: 10.1103/PhysRevD.91.055026

229. Stech B. Trinification Phenomenology and the structure of Higgs Bosons. *J High Energy Phys*. (2014) **08**:139. doi: 10.1007/JHEP08(2014)139

230. Camargo-Molina JE, Morais AP, Pasechnik R, Wessén J. On a radiative origin of the Standard Model from Trinification. *J High Energy Phys*. (2016) **09**:129. doi: 10.1007/JHEP09(2016)129

231. Camargo-Molina JE, Morais AP, Ordell A, Pasechnik R, Sampaio MOP, Wessén J. Reviving trinification models through an E6 -extended supersymmetric GUT. *Phys Rev*. (2017) **D95**:075031. doi: 10.1103/PhysRevD.95.075031

232. Camargo-Molina JE, Morais AP, Ordell A, Pasechnik R, Wessén J. Scale hierarchies, symmetry breaking and particle spectra in SU(3)-family extended SUSY trinification. *Phys Rev*. (2019) **D99**:035041. doi: 10.1103/PhysRevD.99.035041

233. King SF, Moretti S, Nevzorov R. Theory and phenomenology of an exceptional supersymmetric standard model. *Phys Rev*. (2006) **D73**:035009. doi: 10.1103/PhysRevD.73.035009

234. King SF, Moretti S, Nevzorov R. Exceptional supersymmetric standard model. *Phys Lett*. (2006) **B634**:278–84. doi: 10.1016/j.physletb.2005.12.070

235. Athron P, Hall JP, Howl R, King SF, Miller DJ, Moretti S, et al. Aspects of the exceptional supersymmetric standard model. *Nucl Phys Proc Suppl*. (2010) **200-202**:120–9. doi: 10.1016/j.nuclphysbps.2010.02.074

236. King SF, Luo R, Miller DJ, Nevzorov R. Leptogenesis in the exceptional supersymmetric standard model: flavour dependent lepton asymmetries. *J High Energy Phys*. (2008) **12**:042. doi: 10.1088/1126-6708/2008/12/042

237. Nevzorov R. Leptogenesis as an origin of hot dark matter and baryon asymmetry in the *E*_{6} inspired SUSY models. *Phys Lett*. (2018) **B779**:223–9. doi: 10.1016/j.physletb.2018.02.020

238. Nevzorov R. *E*_{6} inspired supersymmetric models with exact custodial symmetry. *Phys Rev*. (2013) **D87**:015029. doi: 10.1103/PhysRevD.87.015029

239. Nevzorov R. Quasifixed point scenarios and the Higgs mass in the E6 inspired supersymmetric models. *Phys Rev*. (2014) **D89**:055010. doi: 10.1103/PhysRevD.89.055010

240. Howl R, King SF. Planck scale unification in a supersymmetric standard model. *Phys Lett*. (2007) **B652**:331–7. doi: 10.1016/j.physletb.2007.07.035

241. Athron P, King SF, Miller DJ, Moretti S, Nevzorov R. The constrained exceptional supersymmetric standard model. *Phys Rev*. (2009) **D80**:035009. doi: 10.1103/PhysRevD.80.035009

242. Athron P, King SF, Miller DJ, Moretti S, Nevzorov R. Predictions of the constrained exceptional supersymmetric standard model. *Phys Lett*. (2009) **B681**:448–56. doi: 10.1016/j.physletb.2009.10.051

243. Athron P, King SF, Miller DJ, Moretti S, Nevzorov R. Constrained Exceptional Supersymmetric Standard Model with a Higgs Near 125 GeV. *Phys Rev*. (2012) **D86**:095003. doi: 10.1103/PhysRevD.86.095003

244. Athron P, Mühlleitner M, Nevzorov R, Williams AG. Non-Standard Higgs Decays in U(1) Extensions of the MSSM. *J High Energy Phys*. (2015) **01**:153. doi: 10.1007/JHEP01(2015)153

245. Athron P, King SF, Miller DJ, Moretti S, Nevzorov R. LHC signatures of the constrained exceptional supersymmetric standard model. *Phys Rev*. (2011) **D84**:055006. doi: 10.1103/PhysRevD.84.055006

246. Hall JP, King SF. Neutralino dark matter with inert higgsinos and singlinos. *J High Energy Phys*. (2009) **08**:088. doi: 10.1088/1126-6708/2009/08/088

247. Hall JP, King SF. Bino dark matter and big bang nucleosynthesis in the constrained *E*_{6}SSM with massless inert singlinos. *J High Energy Phys*. (2011) **06**:006. doi: 10.1007/JHEP06(2011)006

248. Hall JP, King SF, Nevzorov R, Pakvasa S, Sher M, Nevzorov R, et al. Novel Higgs Decays and Dark Matter in the E(6)SSM. *Phys Rev*. (2011) **D83**:075013. doi: 10.1103/PhysRevD.83.075013

249. Athron P, Harries D, Nevzorov R, Williams AG. *E*_{6} Inspired SUSY benchmarks, dark matter relic density and a 125 GeV Higgs. *Phys Lett*. (2016) **B760**:19–25. doi: 10.1016/j.physletb.2016.06.040

250. Athron P, Harries D, Nevzorov R, Williams AG. Dark matter in a constrained E_{6} inspired SUSY model. *J High Energy Phys*. (2016) **12**:128. doi: 10.1007/JHEP12(2016)128

251. Rocher J, Sakellariadou M. Supersymmetric grand unified theories and cosmology. *J Cosmol Astropart Phys*. (2005) **0503**:004. doi: 10.1088/1475-7516/2005/03/004

252. Lyth DH. What would we learn by detecting a gravitational wave signal in the cosmic microwave background anisotropy? *Phys Rev Lett*. (1997) **78**:1861–3. doi: 10.1103/PhysRevLett.78.1861

253. Ade PAR, Aghanim N, Arnaud M, Arroja F, Ashdown M, Aumont J, et al. Planck 2015 results. XX. Constraints on inflation. *Astron Astrophys*. (2016) **594**:A20. doi: 10.1051/0004-6361/201525898

254. Akrami Y, Arroja F, Ashdown M, Aumont J, Baccigalupi C, Ballardini M, et al. Planck 2018 results. X. Constraints on inflation. *arXiv:1807.06211* (2018).

255. Coleman SR, Weinberg EJ. Radiative corrections as the origin of spontaneous symmetry breaking. *Phys Rev*. (1973) **D7**:1888–910. doi: 10.1103/PhysRevD.7.1888

256. Croon D, Sanz V. Saving natural inflation. *J Cosmol Astropart Phys*. (2015) **1502**:008. doi: 10.1088/1475-7516/2015/02/008

257. Croon D, Sanz V, Setford J. Goldstone inflation. *J High Energy Phys*. (2015) **10**:020. doi: 10.1007/JHEP10(2015)020

258. Croon D, Sanz V, Tarrant ERM. Reheating with a composite Higgs boson. *Phys Rev*. (2016) **D94**:045010. doi: 10.1103/PhysRevD.94.045010

259. Shafi Q, Vilenkin A. Inflation with SU(5). *Phys Rev Lett*. (1984) **52**:691–4. doi: 10.1103/PhysRevLett.52.691

260. Barenboim G, Chun EJ, Lee HM. Coleman-Weinberg Inflation in light of Planck. *Phys Lett*. (2014) **B730**:81–88. doi: 10.1016/j.physletb.2014.01.039

261. Okada N, Senoguz VN, Shafi Q. The observational status of simple inflationary models: an update. Turk *J Phys*. (2016) **40**:150–62. doi: 10.3906/fiz-1505-7

262. Cerioni A, Finelli F, Tronconi A, Venturi G. Inflation and reheating in induced gravity. *Phys Lett*. (2009) **B681**:383–6. doi: 10.1016/j.physletb.2009.10.066

263. Panotopoulos G. Nonminimal GUT inflation after Planck results. *Phys Rev*. (2014) **D89**:047301. doi: 10.1103/PhysRevD.89.047301

264. Karam A, Pappas T, Tamvakis K. Nonminimal coleman–weinberg inflation with an *R*^{2} term. *J Cosmol Astropar Phys.* (2018) **2019**:006. doi: 10.1088/1475-7516/2019/02/006

265. Ellis J, Nanopoulos DV, Olive KA. No-Scale supergravity realization of the starobinsky model of inflation. *Phys Rev Lett*. (2013) **111**:111301. doi: 10.1103/PhysRevLett.111.129902

266. Starobinsky AA. A new type of isotropic cosmological models without singularity. *Phys Lett*. (1980) **B91**:99–102. doi: 10.1016/0370-2693(80)90670-X

267. Whitt B. Fourth order gravity as general relativity plus matter. *Phys Lett*. (1984) **145B**:176–8. doi: 10.1016/0370-2693(84)90332-0

268. Croon D, Ellis J, Mavromatos NE. Wess-zumino inflation in light of planck. *Phys Lett*. (2013) **B724**:165–9. doi: 10.1016/j.physletb.2013.06.016

269. Ellis J, Nanopoulos DV, Olive KA. A no-scale supergravity framework for sub-Planckian physics. *Phys Rev*. (2014) **D89**:043502. doi: 10.1103/PhysRevD.89.043502

270. Ellis J, Gonzalo TE, Harz J, Huang WC. Flipped GUT Inflation. *J Cosmol Astropart Phys*. (2015) **1503**:039. doi: 10.1088/1475-7516/2015/03/039

271. Ellis J, Garcia MAG, Nagata N, Nanopoulos DV, Olive KA. Starobinsky-like inflation and neutrino masses in a No-Scale SO(10) model. *J Cosmol Astropart Phys*. (2016) **1611**:018. doi: 10.1088/1475-7516/2016/11/018

272. Gonzalo TE, Heurtier L, Moursy A. Sneutrino driven GUT inflation in supergravity. *J High Energy Phys*. (2017) **06**:109. doi: 10.1007/JHEP06(2017)109

273. Ellis J, He HJ, Xianyu ZZ. New higgs inflation in a No-Scale supersymmetric SU(5) GUT. *Phys Rev*. (2015) **D91**:021302. doi: 10.1103/PhysRevD.91.021302

274. Ellis J, He HJ, Xianyu ZZ. Higgs inflation, reheating and gravitino production in No-Scale supersymmetric GUTs. *J Cosmol Astropart Phys*. (2016) **1608**:068. doi: 10.1088/1475-7516/2016/08/068

275. Hertzberg MP, Wilczek F. Inflation driven by unification energy. *Phys Rev*. (2017) **D95**:063516. doi: 10.1103/PhysRevD.95.063516

276. Rehman MU, Shafi Q, Wickman JR. GUT inflation and proton decay after WMAP5. *Phys Rev*. (2008) **D78**:123516. doi: 10.1103/PhysRevD.78.123516

277. Boucenna SM, Morisi S, Shafi Q, Valle JWF. Inflation and majoron dark matter in the seesaw mechanism. *Phys Rev*. (2014) **D90**:055023. doi: 10.1103/PhysRevD.90.055023

278. Sravan Kumar K, Vargas Moniz P. Conformal GUT inflation, proton lifetime and non-thermal leptogenesis. *arXiv:1806.09032* (2018).

279. Kibble TWB. Topology of cosmic domains and strings. *J Phys*. (1976) **A9**:1387–98. doi: 10.1088/0305-4470/9/8/029

280. Kibble TWB. Some implications of a cosmological phase transition. *Phys Rept*. (1980) **67**:183. doi: 10.1016/0370-1573(80)90091-5

281. Jeannerot R, Rocher J, Sakellariadou M. How generic is cosmic string formation in SUSY GUTs. *Phys Rev*. (2003) **D68**:103514. doi: 10.1103/PhysRevD.68.103514

282. Achucarro A, Martins CJAP. Cosmic strings. In: Meyers, RA editor. *Encyclopedia of Complexity and Systems Science.* New York, NY: Springer (2009). doi: 10.1007/978-0-387-30440-3_107

283. Vilenkin A. Cosmological density fluctuations produced by vacuum strings. *Phys Rev Lett*. (1981) **46**:1169–72. doi: 10.1103/PhysRevLett.46.1169

284. Brandenberger RH, Turok N. Fluctuations from cosmic strings and the microwave background. *Phys Rev*. (1986) **D33**:2182. doi: 10.1103/PhysRevD.33.2182

285. Gott JR III. Gravitational lensing effects of vacuum strings: exact solutions. *Astrophys J*. (1985) **288**:422–7. doi: 10.1086/162808

286. Kaiser N, Stebbins A. Microwave anisotropy due to cosmic strings. *Nature*. (1984) **310**:391–3. doi: 10.1038/310391a0

287. Ade PAR, Aghanim N, Arnaud M, Ashdown M, Aumont J, Baccigalupi C, et al. Planck 2015 results. XIII. Cosmological parameters. *Astron Astrophys*. (2016) **594**:A13. doi: 10.1051/0004-6361/201525830

288. Lizarraga J, Urrestilla J, Daverio D, Hindmarsh M, Kunz M. New CMB constraints for Abelian Higgs cosmic strings. *J Cosmol Astropart Phys*. (2016) **1610**:042. doi: 10.1088/1475-7516/2016/10/042

289. Vilenkin A. Gravitational radiation from cosmic strings. *Phys Lett*. (1981) **107B**:47–50. doi: 10.1016/0370-2693(81)91144-8

290. Vachaspati T, Vilenkin A. Gravitational radiation from cosmic strings. *Phys Rev*. (1985) **D31**:3052. doi: 10.1103/PhysRevD.31.3052

291. Caldwell RR, Allen B. Cosmological constraints on cosmic string gravitational radiation. *Phys Rev*. (1992) **D45**:3447–68. doi: 10.1103/PhysRevD.45.3447

292. Siemens X, Mandic V, Creighton J. Gravitational wave stochastic background from cosmic (super)strings. *Phys Rev Lett*. (2007) **98**:111101. doi: 10.1103/PhysRevLett.98.111101

293. Sanidas SA, Battye RA, Stappers BW. Projected constraints on the cosmic (super)string tension with future gravitational wave detection experiments. *Astrophys J*. (2013) **764**:108. doi: 10.1088/0004-637X/764/1/108

294. Damour T, Vilenkin A. Gravitational wave bursts from cosmic strings. *Phys Rev Lett*. (2000) **85**:3761–4. doi: 10.1103/PhysRevLett.85.3761

295. Damour T, Vilenkin A. Gravitational wave bursts from cusps and kinks on cosmic strings. *Phys Rev*. (2001) **D64**:064008. doi: 10.1103/PhysRevD.64.064008

296. Damour T, Vilenkin A. Gravitational radiation from cosmic (super)strings: bursts, stochastic background, and observational windows. *Phys Rev*. (2005) **D71**:063510. doi: 10.1103/PhysRevD.71.063510

297. Siemens X, Creighton J, Maor I, Ray Majumder S, Cannon K, Read J. Gravitational wave bursts from cosmic (super)strings: quantitative analysis and constraints. *Phys Rev*. (2006) **D73**:105001. doi: 10.1103/PhysRevD.73.105001

298. Olum KD, Blanco-Pillado JJ. Field theory simulation of Abelian Higgs cosmic string cusps. *Phys Rev*. (1999) **D60**:023503. doi: 10.1103/PhysRevD.60.023503

299. Aasi J, Abadie J, Abbott BP, Abbott R, Abbott T, Abernathy MR, et al. Constraints on cosmic strings from the LIGO-Virgo gravitational-wave detectors. *Phys Rev Lett*. (2014) **112**:131101. doi: 10.1103/PhysRevLett.112.131101

300. Sanidas SA, Battye RA, Stappers BW. Constraints on cosmic string tension imposed by the limit on the stochastic gravitational wave background from the European Pulsar Timing Array. *Phys Rev*. (2012) **D85**:122003. doi: 10.1103/PhysRevD.85.122003

301. Abbott BP, Abbott R, Abbott TD, Acernese F, Ackley K, Adams C, et al. Constraints on cosmic strings using data from the first Advanced LIGO observing run. *Phys Rev*. (2018) **D97**:102002. doi: 10.1103/PhysRevD.97.102002

302. Blanco-Pillado JJ, Olum KD, Siemens X. New limits on cosmic strings from gravitational wave observation. *Phys Lett*. (2018) **B778**:392–6. doi: 10.1016/j.physletb.2018.01.050

303. Lorenz L, Ringeval C, Sakellariadou M. Cosmic string loop distribution on all length scales and at any redshift. *J Cosmol Astropart Phys*. (2010) **1010**:003. doi: 10.1088/1475-7516/2010/10/003

304. Ringeval C, Sakellariadou M, Bouchet F. Cosmological evolution of cosmic string loops. *J Cosmol Astropart Phys*. (2007) **0702**:023. doi: 10.1088/1475-7516/2007/02/023

305. Ringeval C, Suyama T. Stochastic gravitational waves from cosmic string loops in scaling. *J Cosmol Astropart Phys*. (2017) **1712**:027. doi: 10.1088/1475-7516/2017/12/027

306. Blanco-Pillado JJ, Olum KD, Shlaer B. The number of cosmic string loops. *Phys Rev*. (2014) **D89**:023512. doi: 10.1103/PhysRevD.89.023512

307. Blanco-Pillado JJ, Olum KD. Stochastic gravitational wave background from smoothed cosmic string loops. *Phys Rev*. (2017) **D96**:104046. doi: 10.1103/PhysRevD.96.104046

308. Martins CJAP, Shellard EPS. Fractal properties and small-scale structure of cosmic string networks. *Phys Rev*. (2006) **D73**:043515. doi: 10.1103/PhysRevD.73.043515

309. Hindmarsh M, Stuckey S, Bevis N. Abelian higgs cosmic strings: small scale structure and loops. *Phys Rev*. (2009) **D79**:123504. doi: 10.1103/PhysRevD.79.123504

310. Hindmarsh M, Rummukainen K, Tenkanen TVI, Weir DJ. Improving cosmic string network simulations. *Phys Rev*. (2014) **D90**:043539. doi: 10.1103/PhysRevD.90.043539

311. Hindmarsh M, Lizarraga J, Urrestilla J, Daverio D, Kunz M. Scaling from gauge and scalar radiation in Abelian Higgs string networks. *Phys Rev*. (2017) **D96**:023525. doi: 10.1103/PhysRevD.96.023525

312. Brandenberger RH. On the decay of cosmic string loops. *Nucl Phys*. (1987) **B293**:812–28. doi: 10.1016/0550-3213(87)90092-7

313. Hindmarsh M, Rummukainen K, Weir DJ. Numerical simulations of necklaces in SU(2) gauge-Higgs field theory. *Phys Rev*. (2017) **D95**:063520. doi: 10.1103/PhysRevD.95.063520

314. Hindmarsh M, Rummukainen K, Weir DJ. New solutions for non-Abelian cosmic strings. *Phys Rev Lett*. (2016) **117**:251601. doi: 10.1103/PhysRevLett.117.251601

315. Lopez-Eiguren A, Urrestilla J, Achúcarro A, Avgoustidis A, Martins CJAP. Evolution of Semilocal String Networks: II. Velocity estimators. *Phys Rev*. (2017) **D96**:023526. doi: 10.1103/PhysRevD.96.023526

316. Figueroa DG, Hindmarsh M, Urrestilla J. Exact scale-invariant background of gravitational waves from cosmic defects. *Phys Rev Lett*. (2013) **110**:101302. doi: 10.1103/PhysRevLett.110.101302

317. Patel HH, Ramsey-Musolf MJ, Wise MB. Color Breaking in the early universe. *Phys Rev*. (2013) **D88**:015003. doi: 10.1103/PhysRevD.88.015003

318. Ramsey-Musolf MJ, Winslow P, White G. Color breaking baryogenesis. *Phys Rev*. (2018) **D97**:123509. doi: 10.1103/PhysRevD.97.123509

319. Long AJ, Tesi A, Wang LT. Baryogenesis at a lepton-number-breaking phase transition. *J High Energy Phys*. (2017) **10**:095. doi: 10.1007/JHEP10(2017)095

320. Schmitz K. *The B-L Phase Transition: Implications for Cosmology and Neutrinos*. Hamburg U. (2012). Available online at: http://www-library.desy.de/cgi-bin/showprep.pl?thesis12-039

321. Buchmuller W, Domcke V, Schmitz K. Spontaneous B-L breaking as the origin of the hot early universe. *Nucl Phys*. (2012) **B862**:587–632. doi: 10.1016/j.nuclphysb.2012.05.001

322. Buchmüller W, Domcke V, Kamada K, Schmitz K. The gravitational wave spectrum from cosmological *B*−*L* breaking. *J Cosmol Astropart Phys*. (2013) **1310**:003. doi: 10.1088/1475-7516/2013/10/003

323. Croon D, Sanz V, White G. Model discrimination in gravitational wave spectra from dark phase transitions. *J High Energy Phys*. (2018) **08**:203. doi: 10.1007/JHEP08(2018)203

324. Schwaller P. Gravitational waves from a dark phase transition. *Phys Rev Lett*. (2015) **115**:181101. doi: 10.1103/PhysRevLett.115.181101

325. Breitbach M, Kopp J, Madge E, Opferkuch T, Schwaller P. Dark, cold, and noisy: constraining secluded hidden sectors with gravitational waves. *arXiv:1811.11175* (2018).

326. Baldes I, Garcia-Cely C. Strong gravitational radiation from a simple dark matter model. *arXiv:1809.01198* (2018).

327. Croon D, Gonzalo TE, White G. Gravitational waves from a pati-salam phase transition. *J High Energy Phys*. (2019) **02**:083. doi: 10.1007/JHEP02(2019)083

328. Mazumdar A, White G. Cosmic phase transitions: their applications and experimental signatures. *J High Energy Phys.* (2018) doi: 10.1088/1361-6633/ab1f55. [Epub ahead of print].

329. Caprini C, Figueroa DG. Cosmological backgrounds of gravitational waves. *Class Quant Grav*. (2018) **35**:163001. doi: 10.1088/1361-6382/aac608

330. Weir DJ. Gravitational waves from a first order electroweak phase transition: a brief review. *Phil Trans Roy Soc Lond*. (2018) **A376**:20170126. doi: 10.1098/rsta.2017.0126

331. Chao W. Electroweak baryogenesis in the exceptional supersymmetric standard model. *J Cosmol Astropart Phys*. (2015) **1508**:055. doi: 10.1088/1475-7516/2015/08/055

332. Chiang CW, Ramsey-Musolf MJ, Senaha E. Standard model with a complex scalar singlet: cosmological implications and theoretical considerations. *Phys Rev*. (2018) **D97**:015005. doi: 10.1103/PhysRevD.97.015005

333. Profumo S, Ramsey-Musolf MJ, Wainwright CL, Winslow P. Singlet-catalyzed electroweak phase transitions and precision Higgs boson studies. *Phys Rev*. (2015) **D91**:035018. doi: 10.1103/PhysRevD.91.035018

334. Caprini C, Hindmarsh M, Huber S, Konstandin T, Kozaczuk J, Nardini G, et al. Science with the space-based interferometer eLISA. II: gravitational waves from cosmological phase transitions. *J Cosmol Astropart Phys*. (2016) **1604**:001. doi: 10.1088/1475-7516/2016/04/001

335. Bodeker D, Moore GD. Electroweak bubble wall speed limit. *J Cosmol Astropart Phys*. (2017) **1705**:025. doi: 10.1088/1475-7516/2017/05/025

336. Hindmarsh M, Huber SJ, Rummukainen K, Weir DJ. Shape of the acoustic gravitational wave power spectrum from a first order phase transition. *Phys Rev*. (2017) **D96**:103520. doi: 10.1103/PhysRevD.96.103520

337. Caprini C, Durrer R, Servant G. The stochastic gravitational wave background from turbulence and magnetic fields generated by a first-order phase transition. *J Cosmol Astropart Phys*. (2009) **0912**:024. doi: 10.1088/1475-7516/2009/12/024

338. Punturo M, Abernathy M, Acernese F, Allen B, Andersson N, Arun K, et al. The einstein telescope: a third-generation gravitational wave observatory. *Class Quant Grav*. (2010) **27**:194002. doi: 10.1088/0264-9381/27/19/194002

339. Akutsu T, Ando M, Arai K, Arai Y, Araki S, Araya A, et al. KAGRA: 2.5 Generation interferometric gravitational wave detector. *Nat Astron.* (2018) **3**:35–40. doi: 10.1038/s41550-018-0658-y

340. Abbott BP, Abbott R, Abbott TD, Abernathy MR, Ackley K, Adams C, et al. Exploring the Sensitivity of Next Generation Gravitational Wave Detectors. *Class Quant Grav*. (2017) **34**:044001. doi: 10.1088/1361-6382/aa51f4

341. Thrane E, Romano JD. Sensitivity curves for searches for gravitational-wave backgrounds. *Phys Rev*. (2013) **D88**:124032. doi: 10.1103/PhysRevD.88.124032

342. Espinosa JR, Konstandin T, No JM, Servant G. Energy budget of cosmological first-order phase transitions. *J Cosmol Astropart Phys*. (2010) **1006**:028. doi: 10.1088/1475-7516/2010/06/028

343. Grojean C, Servant G, Wells JD. First-order electroweak phase transition in the standard model with a low cutoff. *Phys Rev*. (2005) **D71**:036001. doi: 10.1103/PhysRevD.71.036001

344. Ellis J, Lewicki M, No JM. On the maximal strength of a first-order electroweak phase transition and its gravitational wave signal. *J Cosmol Astropart Phys*. (2018) **04**:003. doi: 10.1088/1475-7516/2019/04/003

345. White GA. *A Pedagogical Introduction to Electroweak Baryogenesis*. IOP Concise Physics. Morgan & Claypool (2016). doi: 10.1088/978-1-6817-4457-5ch1

346. Morrissey DE, Ramsey-Musolf MJ. Electroweak baryogenesis. New *J Phys*. (2012) **14**:125003. doi: 10.1088/1367-2630/14/12/125003

347. Aghanim N, Akrami Y, Ashdown M, Aumont J, Baccigalupi C, Ballardini M, et al. Planck 2018 results. VI. Cosmological parameters. (2018) .

348. Riemer-Sørensen S, Jenssen ES. Nucleosynthesis predictions and high-precision deuterium measurements. *Universe*. (2017) **3**:44. doi: 10.3390/universe3020044

349. Sakharov AD. Violation of CP Invariance, C asymmetry, and baryon asymmetry of the universe. *Pisma Zh Eksp Teor Fiz*. (1967) **5**:32–5.

350. Riotto A. Theories of baryogenesis. In: *Proceedings, Summer School in High-energy physics and cosmology: Trieste, Italy, June 29-July 17, 1998* Trieste (1998). p. 326–436.

351. Gu PH, Sarkar U. SO(10) GUT baryogenesis. *Phys Lett*. (2008) **B663**:80–2. doi: 10.1016/j.physletb.2008.04.001

352. Gehrlein J, Petcov ST, Spinrath M, Zhang X. Leptogenesis in an SU(5) × A_{5} golden ratio flavour model. *Nucl Phys*. (2015) **B896**:311–29. doi: 10.1016/j.nuclphysb.2015.08.019

353. Falcone D, Tramontano F. Leptogenesis with SU(5) inspired mass matrices. *Phys Lett*. (2001) **B506**:1–6. doi: 10.1016/S0370-2693(01)00403-8

354. Pilaftsis A, Underwood TEJ. Resonant leptogenesis. *Nucl Phys*. (2004) **B692**:303–45. doi: 10.1016/j.nuclphysb.2004.05.029

355. Blanchet S, Fileviez Perez P. Baryogenesis via leptogenesis in adjoint SU(5). *J Cosmol Astropart Phys*. (2008) **0808**:037. doi: 10.1088/1475-7516/2008/08/037

356. Chianese M, Di Bari P. Strong thermal *SO*(10)-inspired leptogenesis in the light of recent results from long-baseline neutrino experiments. *J High Energy Phys*. (2018) **05**:073. doi: 10.1007/JHEP05(2018)073

357. Di Bari P, Re Fiorentin M. A full analytic solution of *SO*(10)-inspired leptogenesis. *J High Energy Phys*. (2017) **10**:029. doi: 10.1007/JHEP10(2017)029

358. Di Bari P, Re Fiorentin M. Supersymmetric *SO*(10)-inspired leptogenesis and a new *N*_{2}-dominated scenario. *J Cosmol Astropart Phys*. (2016) **1603**:039. doi: 10.1088/1475-7516/2016/03/039

359. Di Bari P, Marzola L, Re Fiorentin M. Decrypting *SO*(10)-inspired leptogenesis. *Nucl Phys*. (2015) **B893**:122–57. doi: 10.1016/j.nuclphysb.2015.02.005

360. Di Bari P, Marzola L. SO(10)-inspired solution to the problem of the initial conditions in leptogenesis. *Nucl Phys*. (2013) **B877**:719–51. doi: 10.1016/j.nuclphysb.2013.10.027

361. Esteban I, Gonzalez-Garcia MC, Hernandez-Cabezudo A, Maltoni M, Schwetz T. Global analysis of three-flavour neutrino oscillations: synergies and tensions in the determination of θ_{2}3, δ_{C}*P*, and the mass ordering. *J High Energy Phys*. (2018) **01**:106. doi: 10.1007/JHEP01(2019)106

362. Björkeroth F, de Anda FJ, de Medeiros Varzielas I, King SF. Leptogenesis in a Δ(27) × *SO*(10) SUSY GUT. *J High Energy Phys*. (2017) **01**:077. doi: 10.1007/JHEP01(2017)077

363. Björkeroth F, de Anda FJ, de Medeiros Varzielas I, King SF. Towards a complete A_{4}× SU(5) SUSY GUT. *J High Energy Phys*. (2015) **06**:141. doi: 10.1007/JHEP06(2015)141

364. Anisimov A, Di Bari P. Cold dark matter from heavy right-handed neutrino mixing. *Phys Rev*. (2009) **D80**:073017. doi: 10.1103/PhysRevD.80.073017

365. De Roeck A, Ellis J, Grojean C, Heinemeyer S, Jakobs K, Weiglein G, et al. From the LHC to Future Colliders. *Eur Phys J*. (2010) **C66**:525–83. doi: 10.1140/epjc/s10052-010-1244-3

366. Aarons G, Djouadi A, Lykken J, Moenig K, Okada Y, Oreglia M, et al. *International Linear Collider*. Reference Design Report. Vol. 2. physics at the ILC. arXiv:0709.1893 (2007).

367. Baer H, Barklow T, Fujii K, Gao Y, Hoang A, Kanemura S, et al. *The International Linear Collider*. Technical Design Report. Vol. 2. Physics. (2013). doi: 10.2172/1095028

368. Linssen L, Miyamoto A, Stanitzki M, Weerts H. Physics and Detectors at CLIC: CLIC Conceptual Design Report. *arXiv:1202.5940* (2012).

369. Haber HE, Kane GL. The search for supersymmetry: probing physics beyond the standard model. *Phys Rept*. (1985) **117**:75–263. doi: 10.1016/0370-1573(85)90051-1

370. Borschensky C, Krämer M, Kulesza A, Mangano M, Padhi S, Plehn T, et al. Squark and gluino production cross sections in pp collisions at $\sqrt{s}$ = 13, 14, 33 and 100 TeV. *Eur Phys J*. (2014) **C74**:3174. doi: 10.1140/epjc/s10052-014-3174-y

371. Beenakker W, Borschensky C, Krämer M, Kulesza A, Laenen E. NNLL-fast: predictions for coloured supersymmetric particle production at the LHC with threshold and Coulomb resummation. *J High Energy Phys*. (2016) **12**:133. doi: 10.1007/JHEP12(2016)133

372. Aaboud M, Abbott B, Abdinov O, Abeloos B, Abidi SH, AbouZeid O, et al. Search for squarks and gluinos in final states with jets and missing transverse momentum using 36 fb^{−1} of $\sqrt{s}=13$ TeV pp collision data with the ATLAS detector. *Phys Rev*. (2018) **D97**:112001. doi: 10.1103/PhysRevD.97.112001

373. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for natural and split supersymmetry in proton-proton collisions at $\sqrt{s}=13$ TeV in final states with jets and missing transverse momentum. *J High Energy Phys*. (2018) **05**:025. doi: 10.1007/JHEP05(2018)025

374. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for supersymmetry in multijet events with missing transverse momentum in proton-proton collisions at 13 TeV. *Phys Rev*. (2017) **D96**:032003. doi: 10.1103/PhysRevD.96.032003

375. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for squarks and gluinos in events with an isolated lepton, jets, and missing transverse momentum at $\sqrt{s}=13$ TeV with the ATLAS detector. *Phys Rev*. (2017) **D96**:112010. doi: 10.1103/PhysRevD.96.112010

376. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for supersymmetry in events with one lepton and multiple jets exploiting the angular correlation between the lepton and the missing transverse momentum in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *Phys Lett*. (2018) **B780**:384–409. doi: 10.1016/j.physletb.2018.03.028

377. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for Supersymmetry in *pp* Collisions at $\sqrt{s}=13\text{TeV}$ in the Single-Lepton Final State Using the Sum of Masses of Large-Radius Jets. *Phys Rev Lett*. (2017) **119**:151802. doi: 10.1103/PhysRevLett.119.151802

378. Aaboud M, Georges Aad (Marseille C, Abdinov O, Abeloos B, Abidi SH, et al. Search for new phenomena using the invariant mass distribution of same-flavour opposite-sign dilepton pairs in events with missing transverse momentum in $\sqrt{s}=13$ TeV pp collisions with the ATLAS detector. *Eur Phys J*. (2018) **C78**:625. doi: 10.1140/epjc/s10052-018-6081-9

379. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for top squarks and dark matter particles in opposite-charge dilepton final states at $\sqrt{s}=$ 13 TeV. *Phys Rev*. (2018) **D97**:032009. doi: 10.1103/PhysRevD.97.032009

380. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for supersymmetry in final states with two same-sign or three leptons and jets using 36 fb^{−1} of $\sqrt{s}=13$ TeV *pp* collision data with the ATLAS detector. *J High Energy Phys*. (2017) **09**:084. doi: 10.1007/JHEP09(2017)084

381. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for supersymmetry in events with at least three electrons or muons, jets, and missing transverse momentum in proton-proton collisions at $\sqrt{s}=13$ TeV. *J High Energy Phys*. (2018) **02**:067. doi: 10.1007/jhep02(2018)067

382. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for supersymmetry in final states with missing transverse momentum and multiple *b*-jets in proton-proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *J High Energy Phys*. (2018) **06**:107. doi: 10.1007/JHEP06(2018)107

383. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for physics beyond the standard model in events with high-momentum higgs bosons and missing transverse momentum in proton-proton collisions at 13 TeV. *Phys Rev Lett*. (2018) **120**:241801. doi: 10.1103/PhysRevLett.120.241801

384. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for supersymmetry in proton-proton collisions at 13 TeV using identified top quarks. *Phys Rev*. (2018) **D97**:012007. doi: 10.1103/PhysRevD.97.012007

385. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for squarks and gluinos in final states with hadronically decaying τ-leptons, jets, and missing transverse momentum using *pp* collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector. *Phys Rev*. (2018) **99**:012009. doi: 10.1103/PhysRevD.99.012009

386. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for a scalar partner of the top quark in the jets plus missing transverse momentum final state at $\sqrt{s}$=13 TeV with the ATLAS detector. *J High Energy Phys*. (2017) **12**:085. doi: 10.1007/JHEP12(2017)085

387. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for top squark pair production in pp collisions at $\sqrt{s}=13$ TeV using single lepton events. *J High Energy Phys*. (2017) **10**:019. doi: 10.1007/JHEP10(2017)019

388. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for supersymmetry in events with *b*-tagged jets and missing transverse momentum in *pp* collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *J High Energy Phys*. (2017) **11**:195. doi: 10.1007/JHEP11(2017)195

389. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for the pair production of third-generation squarks with two-body decays to a bottom or charm quark and a neutralino in proton-proton collisions at $\sqrt{s}$ = 13 TeV. *Phys Lett*. (2018) **B778**:263–91. doi: 10.1016/j.physletb.2018.01.012

390. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for top-squark pair production in final states with one lepton, jets, and missing transverse momentum using 36 fb^{−1} of $\sqrt{s}=13$ TeV pp collision data with the ATLAS detector. *J High Energy Phys*. (2018) **06**:108. doi: 10.1007/JHEP09(2018)050

391. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for top squarks decaying via four-body or chargino-mediated modes in single-lepton final states in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *J High Energy Phys*. (2018) **09**:065. doi: 10.1007/JHEP09(2018)065

392. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for direct top squark pair production in final states with two leptons in $\sqrt{s}=13$ TeV *pp* collisions with the ATLAS detector. *Eur Phys J*. (2017) **C77**:898.

393. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Searches for pair production of charginos and top squarks in final states with two oppositely charged leptons in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *J High Energy Phys*. (2018) **11**:079. doi: 10.1007/JHEP11(2018)079

394. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for direct top squark pair production in events with a Higgs or *Z* boson, and missing transverse momentum in $\sqrt{s}=13$ TeV *pp* collisions with the ATLAS detector. *J High Energy Phys*. (2017) **08**:006. doi: 10.1007/JHEP11(2017)195

395. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for electroweak production of supersymmetric states in scenarios with compressed mass spectra at $\sqrt{s}=13$ TeV with the ATLAS detector. *Phys Rev*. (2018) **D97**:052010. doi: 10.1103/PhysRevD.97.052010

396. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for supersymmetric partners of electrons and muons in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *Phys Lett*. (2018) **B790**:140–166. doi: 10.1016/j.physletb.2019.01.005

397. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Combined search for electroweak production of charginos and neutralinos in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *J High Energy Phys*. (2018) **03**:160. doi: 10.1007/JHEP03(2018)160

398. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for electroweak production of supersymmetric particles in final states with two or three leptons at $\sqrt{s}=13$TeV with the ATLAS detector. (2018) **C78**:995. doi: 10.1140/epjc/s10052-018-6423-7

399. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for chargino-neutralino production using recursive jigsaw reconstruction in final states with two or three charged leptons in proton-proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *Phys Rev*. (2018) **D98**:092012. doi: 10.1103/PhysRevD.98.092012

400. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for supersymmetry in events with four or more leptons in $\sqrt{s}=13$ TeV *pp* collisions with ATLAS. *Phys Rev*. (2018) **D98**:032009. doi: 10.1103/PhysRevD.98.032009

401. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for the direct production of charginos and neutralinos in final states with tau leptons in $\sqrt{s}=$ 13 TeV *pp* collisions with the ATLAS detector. *Eur Phys J*. (2018) **C78**:154. doi: 10.1140/epjc/s10052-018-5583-9

402. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for supersymmetry in events with a τ lepton pair and missing transverse momentum in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *J High Energy Phys*. (2018) **11**:151. doi: 10.1007/JHEP11(2018)151

403. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for pair production of higgsinos in final states with at least three *b*-tagged jets in $\sqrt{s}=13$ TeV *pp* collisions using the ATLAS detector. *Phys Rev*. (2018) **D98**:092002. doi: 10.1140/epjc/s10052-018-6457-x

404. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for Higgs boson pair production in the *WW*^{(*)}*WW*^{(*)} decay channel using ATLAS data recorded at $\sqrt{s}=13$ TeV. *arXiv:1811.11028* (2018).

405. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for heavy ZZ resonances in the ℓ^{+}ℓ^{−}ℓ^{+}ℓ^{−} and ${\ell}^{+}{\ell}^{-}\nu \overline{\nu}$ final states using proton-proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *Eur Phys J*. (2018) **C78**:293. doi: 10.1140/epjc/s10052-018-5686-3

406. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for heavy resonances decaying into *WW* in the *eνμν* final state in *pp* collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *Eur Phys J*. (2018) **C78**:24. doi: 10.1140/epjc/s10052-017-5491-4

407. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. A search for pair production of new light bosons decaying into muons in proton-proton collisions at 13 TeV. *arXiv:1812.00380* (2018).

408. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for pair production of Higgs bosons in the $b\overline{b}b\overline{b}$ final state using proton-proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *J High Energy Phys*. (2018) **01**:030. doi: 10.1007/JHEP01(2019)030

409. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for beyond the standard model Higgs bosons decaying into a $\text{b}\overline{\text{b}}$ pair in pp collisions at $\sqrt{s}=$ 13 TeV. *J High Energy Phys*. (2018) **08**:113. doi: 10.1007/JHEP08(2018)113

410. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for Higgs boson pair production in the $b\overline{b}W{W}^{*}$ decay mode at $\sqrt{s}=13$ TeV with the ATLAS detector. *J High Energy Phys.* (2018) **04**:092 . doi: 10.1007/JHEP04(2019)092

411. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for Higgs boson pair production in the $\gamma \gamma b\overline{b}$ final state with 13 TeV *pp* collision data collected by the ATLAS experiment. *J High Energy Phys*. (2018) **11**:040. doi: 10.1007/JHEP11(2018)040

412. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for a standard model-like Higgs boson in the mass range between 70 and 110 GeV in the diphoton final state in proton-proton collisions at $\sqrt{s}=$ 8 and 13 TeV. *Phys Lett*. (2018) **B793**:320–47. doi: 10.1016/j.physletb.2019.03.064

413. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for Resonant and Nonresonant Higgs Boson Pair Production in the $b\overline{b}{\tau}^{+}{\tau}^{-}$ Decay Channel in *pp* Collisions at $\sqrt{s}=13$TeV with the ATLAS Detector. *Phys Rev Lett*. (2018) **121**:191801. doi: 10.1103/PhysRevLett.121.191801

414. Khachatryan V, Sirunyan AM, Tumasyan A, Wolfgang Adam EA, Bergauer T, et al. Search for light bosons in decays of the 125 GeV Higgs boson in proton-proton collisions at $\sqrt{s}=8$ TeV. *J High Energy Phys*. (2017) **10**:076. doi: 10.1007/JHEP10(2017)076

415. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for a light pseudoscalar Higgs boson produced in association with bottom quarks in pp collisions at $\sqrt{s}=8$ TeV. *J High Energy Phys*. (2017) **11**:010. doi: 10.1007/JHEP11(2017)010

416. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for charged Higgs bosons decaying into top and bottom quarks at $\sqrt{s}$ = 13 TeV with the ATLAS detector. *J High Energy Phys*. (2018) **11**:085. doi: 10.1007/JHEP11(2018)085

417. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for additional heavy neutral Higgs and gauge bosons in the ditau final state produced in 36 fb^{−1} of pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *J High Energy Phys*. (2018) **01**:055. doi: 10.1007/JHEP11(2018)051

418. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for additional neutral MSSM Higgs bosons in the τ*τ* final state in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *J High Energy Phys*. (2018) **09**:007. doi: 10.1007/JHEP09(2018)007

419. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for new phenomena in high-mass diphoton final states using 37 fb^{−1} of proton–proton collisions collected at $\sqrt{s}=13$ TeV with the ATLAS detector. *Phys Lett*. (2017) **B775**:105–25. doi: 10.1016/j.physletb.2017.10.039

420. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for Higgs boson pair production in the γ*γWW*^{*} channel using *pp* collision data recorded at $\sqrt{s}=13$ TeV with the ATLAS detector. *Eur Phys J*. (2018) **C78**:1007. doi: 10.1103/PhysRevD.98.092002

421. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. A search for pairs of highly collimated photon-jets in *pp* collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector. *Phys Rev*. (2018) **D99**:012008. doi: 10.1103/PhysRevD.99.012008

422. Aad G, Abbott B, Abdallah J, Abdinov O, Aben R, Abolins M, et al. Search for a charged higgs boson produced in the vector-boson fusion mode with decay *H*^{±}→*W*^{±}*Z* Using *pp* collisions at $\sqrt{s}=8$ TeV with the ATLAS experiment. *Phys Rev Lett*. (2015) **114**:231801. doi: 10.1103/PhysRevLett.114.231801

423. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for Charged Higgs Bosons Produced via Vector Boson Fusion and Decaying into a Pair of *W* and *Z* Bosons Using *pp* Collisions at $\sqrt{s}=13\text{TeV}$. *Phys Rev Lett*. (2017) **119**:141802. doi: 10.1103/PhysRevLett.119.141802

424. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for a charged Higgs boson decaying to charm and bottom quarks in proton-proton collisions at $\sqrt{s}=8$ TeV. *J High Energy Phys*. (2018) **11**:115. doi: 10.1007/JHEP11(2018)115

425. Aaboud M, Aad G, Abbott B, Abdallah J, Abdinov O, Abeloos B, et al. Search for charged Higgs bosons produced in association with a top quark and decaying via *H*^{±} → τ*ν* using *pp* collision data recorded at $\sqrt{s}=13$ TeV by the ATLAS detector. *Phys Lett*. (2016) **B759**:555–74. doi: 10.1016/j.physletb.2016.06.017

426. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for R-parity-violating supersymmetric particles in multi-jet final states produced in *p*-*p* collisions at $\sqrt{s}=13$ TeV using the ATLAS detector at the LHC. *Phys Lett*. (2018) **B785**:136–58. doi: 10.1016/j.physletb.2018.08.021

427. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for new phenomena in a lepton plus high jet multiplicity final state with the ATLAS experiment using $\sqrt{s}=13$ TeV proton-proton collision data. *J High Energy Phys*. (2017) **09**:088. doi: 10.1007/JHEP09(2017)088

428. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. A search for pair-produced resonances in four-jet final states at $\sqrt{s}=$ 13 TeV with the ATLAS detector. *Eur Phys J*. (2018) **C78**:250. doi: 10.1140/epjc/s10052-018-5693-4

429. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for *R*-parity violating supersymmetry in pp collisions at $\sqrt{s}=$ 13 TeV using b jets in a final state with a single lepton, many jets, and high sum of large-radius jet masses. *Phys Lett*. (2018) **B783**:114–39. doi: 10.1016/j.physletb.2018.06.028

430. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for B-L R -parity-violating top squarks in $\sqrt{s}$ =13TeV pp collisions with the ATLAS experiment. *Phys Rev*. (2018) **D97**:032003. doi: 10.1103/PhysRevD.97.032003

431. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for resonant production of second-generation sleptons with same-sign dimuon events in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *Eur Phys J*. (2018) **C79**:305. doi: 10.1140/epjc/s10052-019-6800-x

432. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for long-lived particles produced in *pp* collisions at $\sqrt{s}=13$ TeV that decay into displaced hadronic jets in the ATLAS muon spectrometer. *Phys Rev*. (2018) **D99**:052005. doi: 10.1103/PhysRevD.99.052005

433. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for heavy charged long-lived particles in proton-proton collisions at $\sqrt{s}=13$ TeV using an ionisation measurement with the ATLAS detector. *Phys Lett*. (2019) **B788**:96–116. doi: 10.1016/j.physletb.2018.10.055

434. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for long-lived, massive particles in events with displaced vertices and missing transverse momentum in $\sqrt{s}$ = 13 TeV *pp* collisions with the ATLAS detector. *Phys Rev*. (2018) **D97**:052012. doi: 10.1103/PhysRevD.97.052012

435. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for long-lived charginos based on a disappearing-track signature in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *J High Energy Phys*. (2018) **06**:022. doi: 10.1007/JHEP06(2018)022

436. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for long-lived particles in final states with displaced dimuon vertices in *pp* collisions at $\sqrt{s}=$ 13 TeV with the ATLAS detector. *Phys Rev D.* (2018) **99**:012001. doi: 10.1103/PhysRevD.99.012001

437. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for new long-lived particles at $\sqrt{s}=$ 13 TeV. *Phys Lett*. (2018) **B780**:432–54. doi: 10.1016/j.physletb.2018.03.019

438. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for long-lived particles decaying into displaced jets in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *Phys Rev*. (2018) **D99**:032011. doi: 10.1103/PhysRevD.99.032011

439. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for long-lived particles with displaced vertices in multijet events in proton-proton collisions at $\sqrt{s}=$13 TeV. *Phys Rev*. (2018) **D98**:092011. doi: 10.1103/PhysRevD.98.092011

440. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for gauge-mediated supersymmetry in events with at least one photon and missing transverse momentum in pp collisions at $\sqrt{s}=$ 13 TeV. *Phys Lett*. (2018) **B780**:118–43. doi: 10.1016/j.physletb.2018.02.045

441. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for supersymmetry in events with at least one photon, missing transverse momentum, and large transverse event activity in proton-proton collisions at $\sqrt{s}=13$ TeV. *J High Energy Phys*. (2017) **12**:142. doi: 10.1007/JHEP12(2017)142

442. Aaboud M, Aad G, Abbott B, Abdallah J, Abdinov O, Abeloos B, et al. Search for supersymmetry in a final state containing two photons and missing transverse momentum in $\sqrt{s}$ = 13 TeV *pp* collisions at the LHC using the ATLAS detector. *Eur Phys J*. (2016) **C76**:517. doi: 10.1140/epjc/s10052-016-4344-x

443. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for photonic signatures of gauge-mediated supersymmetry in 13 TeV *pp* collisions with the ATLAS detector. *Phys Rev*. (2018) **D97**:092006. doi: 10.1103/PhysRevD.97.092006

444. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for Higgsino pair production in *pp* collisions at $\sqrt{s}$ = 13TeV in final states with large missing transverse momentum and two Higgs bosons decaying via $H\to b\overline{b}$. *Phys Rev*. (2018) **D97**:032007. doi: 10.1103/PhysRevD.97.032007

445. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for top squarks decaying to tau sleptons in *pp* collisions at $\sqrt{s}=13$ TeV with the ATLAS detector. *Phys Rev*. (2018) **D98**:032008. doi: 10.1103/PhysRevD.98.032008

446. Heister A, Schael S, Barate R, Bruneliere R, Bonis ID, Decamp D, et al. Search for charginos nearly mass degenerate with the lightest neutralino in e+ e- collisions at center-of-mass energies up to 209-GeV. *Phys Lett*. (2002) **B533**:223–36. doi: 10.1016/S0370-2693(02)01584-8

447. Heister A, Schael S, Barate R, Bruneliere R, Bonis ID, Decamp D, et al. Absolute lower limits on the masses of selectrons and sneutrinos in the MSSM. *Phys Lett*. (2002) **B544**:73–88. doi: 10.1016/S0370-2693(02)02471-1

448. Abbiendi G, Ainsley C, Akesson PF, Alexander G, Allison J, Amaral P, et al. Search for chargino and neutralino production at s^{(1/2)} = 192-GeV to 209 GeV at LEP. *Eur Phys J*. (2004) **C35**:1–20. doi: 10.1140/epjc/s2004-01758-8

449. Abbiendi G, Ainsley C, Akesson PF, Alexander G, Allison J, Amaral P, et al. Search for anomalous production of dilepton events with missing transverse momentum in e+ e- collisions at s^{(1/2)} = 183-Gev to 209-GeV. *Eur Phys J*. (2004) **C32**:453–73. doi: 10.1140/epjc/s2003-01466-y

450. Acciarri M, Achard P, Adriani O, Aguilar-Benitez M, Alcaraz J, Alemanni G, et al. Search for charginos and neutralinos in *e*^{+}*e*^{−} collisions at $\sqrt{S}$ = 189-GeV. *Phys Lett*. (2000) **B472**:420–33. doi: 10.1016/S0370-2693(99)01388-X

451. Achard P, Adriani O, Aguilar-Benitez M, Alcaraz J, Alemanni G, Allaby J, et al. Search for scalar leptons and scalar quarks at LEP. *Phys Lett*. (2004) **B580**:37–49. doi: 10.1016/j.physletb.2003.10.010

452. Abdallah J, Abreu P, Adam W, Adzic P, Albrecht T, Alderweireld T, et al. Searches for supersymmetric particles in e+ e- collisions up to 208-GeV and interpretation of the results within the MSSM. *Eur Phys J*. (2003) **C31**:421–79. doi: 10.1140/epjc/s2003-01355-5

453. Roszkowski L, Ruiz de Austri R, Trotta R, Tsai YLS, Varley TA. Global fits of the Non-Universal Higgs Model. *Phys Rev*. (2011) **D83**:015014. doi: 10.1103/PhysRevD.83.015014

454. Bertone G, Cerdeno DG, Fornasa M, Ruiz de Austri R, Strege C, Trotta R. Global fits of the cMSSM including the first LHC and XENON100 data. *J Cosmol Astropart Phys*. (2012) **1201**:015. doi: 10.1088/1475-7516/2012/01/015

455. Bechtle P, Bringmann T, Desch K, Dreiner H, Hamer M, Hensel C, et al. Constrained Supersymmetry after two years of LHC data: a global view with Fittino. *J High Energy Phys*. (2012) **06**:098. doi: 10.1007/JHEP06(2012)098

456. Buchmueller O, Cavanaugh R, Roeck AD, Dolan MJ, Ellis JR, Flacher H, et al. The CMSSM and NUHM1 after LHC Run 1. *Eur Phys J*. (2014) **C74**:2922. doi: 10.1140/epjc/s10052-014-2922-3

457. Buchmueller O, Cavanaugh R, Citron M, Roeck AD, Dolan MJ, Ellis JR, et al. The NUHM2 after LHC Run 1. *Eur Phys J*. (2014) **C74**:3212. doi: 10.1140/epjc/s10052-014-3212-9

458. Han C, Hikasa Ki, Wu L, Yang JM, Zhang Y. Status of CMSSM in light of current LHC Run-2 and LUX data. *Phys Lett*. (2017) **B769**:470–76. doi: 10.1016/j.physletb.2017.04.026

459. Athron P, Balázs C, Bringmann T, Buckley A, Chrza̧szcz M, Conrad J, et al. Global fits of GUT-scale SUSY models with GAMBIT. *Eur Phys J*. (2017) **C77**:824. doi: 10.1140/epjc/s10052-017-5167-0

460. Bagnaschi E, Sakurai K, Borsato M, Buchmueller O, Citron M, Costa JC, et al. Likelihood Analysis of the pMSSM11 in Light of LHC 13-TeV Data. *Eur Phys J*. (2018) **C78**:256. doi: 10.1140/epjc/s10052-018-5697-0

461. Athron P, Balázs C, Bringmann T, Buckley A, Chrza̧szcz M, Conrad J, et al. A global fit of the MSSM with GAMBIT. *Eur Phys J*. (2017) **C77**:879. doi: 10.1140/epjc/s10052-017-5196-8

462. Bertone G, Calore F, Caron S, Ruiz R, Kim JS, Trotta R, et al. Global analysis of the pMSSM in light of the Fermi GeV excess: prospects for the LHC Run-II and astroparticle experiments. *J Cosmol Astropart Phys*. (2016) **1604**:037. doi: 10.1088/1475-7516/2016/04/037

463. Bagnaschi E, Costa JC, Sakurai K, Borsato M, Buchmueller O, Cavanaugh R, et al. Likelihood Analysis of Supersymmetric SU(5) GUTs. *Eur Phys J*. (2017) **C77**:104. doi: 10.1140/epjc/s10052-017-4639-6

464. Arbuzov AB, Awramik M, Czakon M, Freitas A, Grunewald MW, Monig K, et al. ZFITTER: A Semi-analytical program for fermion pair production in e+ e- annihilation, from version 6.21 to version 6.42. *Comput Phys Commun*. (2006) **174**:728–58. doi: 10.1016/j.cpc.2005.12.009

465. Ruiz de Austri R, Trotta R, Roszkowski L. A markov chain monte carlo analysis of the CMSSM. *J High Energy Phys*. (2006) **05**:002. doi: 10.1088/1126-6708/2006/05/002

466. Trotta R, Feroz F, Hobson MP, Roszkowski L, Ruiz de Austri R. The Impact of priors and observables on parameter inferences in the Constrained MSSM. *J High Energy Phys*. (2008) **12**:024. doi: 10.1088/1126-6708/2008/12/024

467. Bechtle P, Desch K, Wienemann P. Fittino, a program for determining MSSM parameters from collider observables using an iterative method. *Comput Phys Commun*. (2006) **174**:47–70. doi: 10.1016/j.cpc.2005.09.002

468. Buchmueller O, Cavanaugh R, De Roeck A, Ellis JR, Flacher H, Heinemeyer S, et al. Likelihood Functions for Supersymmetric Observables in Frequentist Analyses of the CMSSM and NUHM1. *Eur Phys J*. (2009) **C64**:391–415. doi: 10.1140/epjc/s10052-009-1159-z

469. Athron P, Balazs C, Bringmann T, Buckley A, Chrza̧szcz M, et al. GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool. *Eur Phys J*. (2017) **C77**:784. doi: 10.1140/epjc/s10052-017-5321-8

470. Athron P, Balázs C, Dal LA, Edsjö J, Farmer B, Gonzalo TE, et al. SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables. *Eur Phys J*. (2018) **C78**:22. doi: 10.1140/epjc/s10052-017-5390-8

471. Bringmann T, Conrad J, Cornell JM, Dal LA, Edsjö J, Farmer B, et al. DarkBit: a GAMBIT module for computing dark matter observables and likelihoods. *Eur Phys J*. (2017) **C77**:831. doi: 10.1140/epjc/s10052-017-5155-4

472. Martinez GD, McKay J, Farmer B, Scott P, Roebber E, Putze A, et al. Comparison of statistical sampling methods with ScannerBit, the GAMBIT scanning module. *Eur Phys J*. (2017) **C77**:761. doi: 10.1140/epjc/s10052-017-5274-y

473. Balázs C, Buckley A, Dal LA, Farmer B, Jackson P, Krislock A, et al. ColliderBit: a GAMBIT module for the calculation of high-energy collider observables and likelihoods. *Eur Phys J*. (2017) **C77**:795. doi: 10.1140/epjc/s10052-017-5285-8

474. Bernlochner FU, Chrzaszcz M, Dal LA, Farmer B, Jackson P, Kvellestad A, et al. FlavBit: A GAMBIT module for computing flavour observables and likelihoods. *Eur Phys J*. (2017) **C77**:786. doi: 10.1140/epjc/s10052-017-5157-2

475. Athron P, Balázs C, Bringmann T, Buckley A, Chrza̧szcz M, Conrad J, et al. Status of the scalar singlet dark matter model. *Eur Phys J*. (2017) **C77**:568. doi: 10.1140/epjc/s10052-017-5113-1

476. Athron P, Balázs C, Beniwal A, Bloor S, Camargo-Molina JE, Cornell JM, et al. Global analyses of Higgs portal singlet dark matter models using GAMBIT. *Eur Phys J*. (2019) **C79**:38. doi: 10.1140/epjc/s10052-018-6513-6

477. Murayama H, Yanagida T. A viable SU(5) GUT with light leptoquark bosons. Mod *Phys Lett*. (1992) **A7**:147–52. doi: 10.1142/S0217732392000070

478. Biggio C, Bordone M, Di Luzio L, Ridolfi G. Massive vectors and loop observables: the *g*−2 case. *J High Energy Phys*. (2016) **10**:002. doi: 10.1007/JHEP10(2016)002

479. Gabrielli E, Marzola L, Raidal M, Veermäe H. Dark matter and spin-1 milli-charged particles. *J High Energy Phys*. (2015) **08**:150. doi: 10.1007/JHEP08(2015)150

480. Doršner I, Fajfer S, Greljo A, Kamenik JF, Košnik N. Physics of leptoquarks in precision experiments and at particle colliders. *Phys Rept*. (2016) **641**:1–68. doi: 10.1016/j.physrep.2016.06.001

481. Buchmuller W, Ruckl R, Wyler D. Leptoquarks in Lepton-Quark Collisions. *Phys Lett*. (1987) **B191**:442–8. doi: 10.1016/0370-2693(87)90637-X

482. Davidson S, Bailey DC, Campbell BA. Model independent constraints on leptoquarks from rare processes. *Z Phys*. (1994) **C61**:613–44. doi: 10.1007/BF01552629

483. Grifols JA, Mendez A. Pair Production of Colored Hypermesons at Very High-energy. *Phys Rev*. (1982) **D26**:324. doi: 10.1103/PhysRevD.26.324

484. Antoniadis I, Baulieu L, Delduc F. Supersymmetric enhancement factor for the one jet cross-section in $p\overline{p}$ Reactions. *Z Phys*. (1984) **C23**:119. doi: 10.1007/BF01557587

485. Eichten E, Hinchliffe I, Lane KD, Quigg C. Super Collider Physics. *Rev Mod Phys*. (1984) **56**:579–707. doi: 10.1103/RevModPhys.56.579

486. Altarelli G, Ruckl R. Searching for a Scalar Top at the CERN $p\overline{p}$ Collider. *Phys Lett*. (1984) **144B**:126–32. doi: 10.1016/0370-2693(84)90190-4

487. Dawson S, Eichten E, Quigg C. Search for supersymmetric particles in hadron-hadron collisions. *Phys Rev*. (1985) **D31**:1581. doi: 10.1103/PhysRevD.31.1581

488. Blumlein J, Boos E, Kryukov A. Leptoquark pair production in hadronic interactions. *Z Phys*. (1997) **C76**:137–53. doi: 10.1007/s002880050538

489. Dorsner I, Fajfer S, Greljo A. Cornering Scalar Leptoquarks at LHC. *J High Energy Phys*. (2014) **10**:154. doi: 10.1007/JHEP10(2014)154

490. Doršner I, Greljo A. Leptoquark toolbox for precision collider studies. *J High Energy Phys*. (2018) **05**:126. doi: 10.1007/JHEP05(2018)126

491. Plehn T, Spiesberger H, Spira M, Zerwas PM. Formation and decay of scalar leptoquarks / squarks in e p collisions. *Z Phys*. (1997) **C74**:611–4. doi: 10.1007/s002880050426

492. Diaz B, Schmaltz M, Zhong YM. The leptoquark Hunter's guide: Pair production. *J High Energy Phys*. (2017) **10**:097. doi: 10.1007/JHEP10(2017)097

493. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for pair production of first-generation scalar leptoquarks at $\sqrt{s}=$ 13 TeV. *Phys Rev*. (2018) **D99**:052002. doi: 10.1103/PhysRevD.99.052002

494. Aaboud M, Aad G, Abbott B, Abdallah J, Abdinov O, Abeloos B, et al. Search for scalar leptoquarks in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS experiment. New *J Phys*. (2016) **18**:093016. doi: 10.1088/1367-2630/18/9/093016

495. Aad G, Abbott B, Abdallah J, Abdinov O, Aben R, Abolins M, et al. Searches for scalar leptoquarks in pp collisions at $\sqrt{s}$ = 8 TeV with the ATLAS detector. *Eur Phys J*. (2016) **C76**:5. doi: 10.1140/epjc/s10052-015-3823-9

496. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Asilar E, Bergauer T, et al. Search for heavy neutrinos or third-generation leptoquarks in final states with two hadronically decaying τ leptons and two jets in proton-proton collisions at $\sqrt{s}=13$ TeV. *J High Energy Phys*. (2017) **03**:077. doi: 10.1007/JHEP03(2017)077

497. Sirunyan AM, Tumasyan A, Adam W, Asilar E, Bergauer T, Brandstetter J, et al. Search for third-generation scalar leptoquarks and heavy right-handed neutrinos in final states with two tau leptons and two jets in proton-proton collisions at $\sqrt{s}=13$ TeV. *J High Energy Phys*. (2017) **07**:121. doi: 10.1007/JHEP07(2017)121

498. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Constraints on models of scalar and vector leptoquarks decaying to a quark and a neutrino at $\sqrt{s}=$ 13 TeV. *Phys Rev*. (2018) **D98**:032005. doi: 10.1103/PhysRevD.98.032005

499. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for third-generation scalar leptoquarks decaying to a top quark and a τ lepton at $\sqrt{s}=$ 13 TeV. *Eur Phys J*. (2018) **C78**:707. doi: 10.1140/epjc/s10052-018-6143-z

500. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for leptoquarks coupled to third-generation quarks in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *Phys Rev Lett*. (2018) **121**:241802. doi: 10.1103/PhysRevLett.121.241802

501. Khachatryan V, Sirunyan AM, Tumasyan A, Adam W, Asilar E, Bergauer T, et al. Search for single production of scalar leptoquarks in proton-proton collisions at $\sqrt{s}=8$ *TeV. Phys Rev*. (2016) **D93**:032005. doi: 10.1103/PhysRevD.93.032005

502. Deppisch FF, Bhupal Dev PS, Pilaftsis A. Neutrinos and Collider Physics. New *J Phys*. (2015) **17**:075019. doi: 10.1088/1367-2630/17/7/075019

503. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for heavy Majorana or Dirac neutrinos and right-handed *W* gauge bosons in final states with two charged leptons and two jets at $\sqrt{s}$ = 13 TeV with the ATLAS detector. *J High Energy Phys*. (2018) **01**:016. doi: 10.1007/JHEP01(2019)016

504. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for a heavy right-handed W boson and a heavy neutrino in events with two same-flavor leptons and two jets at $\sqrt{s}=$ 13 TeV. *J High Energy Phys*. (2018) **05**:148. doi: 10.1007/JHEP05(2018)148

505. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for heavy neutral leptons in events with three charged leptons in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *Phys Rev Lett*. (2018) **120**:221801. doi: 10.1103/PhysRevLett.120.221801

506. Drewes M, Garbrecht B. Combining experimental and cosmological constraints on heavy neutrinos. *Nucl Phys*. (2017) **B921**:250–315. doi: 10.1016/j.nuclphysb.2017.05.001

507. Bergsma F, Dorenbosch J, Allaby JV, Amaldi U, Barbiellini G, Berger C, et al. A Search for Decays of Heavy Neutrinos in the Mass Range 0.5-GeV to 2.8-GeV. *Phys Lett*. (1986) **166B**:473–8. doi: 10.1016/0370-2693(86)91601-1

508. Bernardi G, Carugno G, Chauveau J, Dicarlo F, Dris M, Dumarchez J, et al. FURTHER LIMITS ON HEAVY NEUTRINO COUPLINGS. *Phys Lett*. (1988) **B203**:332–34. doi: 10.1016/0370-2693(88)90563-1

509. Vilain P, Wilquet G, Petrak S, Beyer R, Flegel W, Grote H, et al. Search for heavy isosinglet neutrinos. *Phys Lett*. (1995) **B343**:453–58. doi: 10.1016/0370-2693(94)01422-9

510. Vaitaitis A, Drucker RB, Formaggio J, Adams T, Alton A, Avvakumov S, et al. Search for neutral heavy leptons in a high-energy neutrino beam. *Phys Rev Lett*. (1999) **83**:4943–6. doi: 10.1103/PhysRevLett.83.4943

511. Adriani O, Aguilar-Benitez M, Ahlen S, Alcaraz J, Aloisio A, Alverson G, et al. Search for isosinglet neutral heavy leptons in Z0 decays. *Phys Lett*. (1992) **B295**:371–82.

512. Abreu P, Adam W, Adye T, Azhinenko I, Alekseev GD, Alemany R, et al. Search for neutral heavy leptons produced in Z decays. *Z Phys*. (1997) **C74**:57–71. doi: 10.1007/s002880050370

513. Sirunyan AM, Tumasyan A, Adam W, Ambrogi F, Asilar E, Bergauer T, et al. Search for a W' boson decaying to a τ lepton and a neutrino in proton-proton collisions at $\sqrt{s}=$ 13 TeV. *Phys Lett*. (2018) **B792**:107. doi: 10.1016/j.physletb.2019.01.069

514. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for new high-mass phenomena in the dilepton final state using 36 fb^{−1} of proton-proton collision data at $\sqrt{s}=13$ TeV with the ATLAS detector. *J High Energy Phys*. (2017) **10**:182. doi: 10.1007/JHEP10(2017)182

515. Collaboration C. Search for a high-mass resonance decaying into a dilepton final state in 13 fb^{−1} of pp collisions at $\sqrt{s}=13\text{TeV}$. (2016). Available online at: http://cds.cern.ch/record/2205764

516. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abidi SH, et al. Search for doubly charged Higgs boson production in multi-lepton final states with the ATLAS detector using proton-proton collisions at $\sqrt{s}=13\text{TeV}$. *Eur Phys J*. (2018) **C78**:199. doi: 10.1140/EPJC/S10052-018-5661-Z

517. Collaboration C. A search for doubly-charged Higgs boson production in three and four lepton final states at $\sqrt{s}=13\text{TeV}$. (2017).

518. Aaboud M, Aad G, Abbott B, Abdinov O, Abeloos B, Abhayasinghe DK, et al. Search for doubly charged scalar bosons decaying into same-sign *W* boson pairs with the ATLAS detector. Submitted to: *Eur Phys J*. (2018) **79**:58. doi: 10.1140/epjc/s10052-018-6500-y

519. Aaboud M, Aad G, Abbott B, Abbott DC, Abdinov O, Abeloos B, et al. Search for heavy long-lived multi-charged particles in proton-proton collisions at $\sqrt{s}$ = 13 TeV using the ATLAS detector. (2018) *Phys Rev D*. **99**:052003. doi: 10.1103/PhysRevD.99.052003

520. Abbott LF, Wise MB. The effective hamiltonian for nucleon decay. *Phys Rev*. (1980) **D22**:2208. doi: 10.1103/PhysRevD.22.2208

521. Marciano WJ, Senjanovic G. Predictions of supersymmetric grand unified theories. *Phys Rev*. (1982) **D25**:3092. doi: 10.1103/PhysRevD.25.3092

522. Dimopoulos S, Raby S, Wilczek F. Proton decay in supersymmetric models. *Phys Lett*. (1982) **112B**:133. doi: 10.1016/0370-2693(82)90313-6

523. Nath P, Chamseddine AH, Arnowitt RL. Nucleon Decay in Supergravity Unified Theories. *Phys Rev*. (1985) **D32**:2348–58. doi: 10.1103/PhysRevD.32.2348

524. Arnowitt RL, Chamseddine AH, Nath P. NUCLEON DECAY BRANCHING RATIOS IN SUPERGRAVITY SU(5) GUTS. *Phys Lett*. (1985) **156B**:215–9. doi: 10.1016/0370-2693(85)91512-6

525. Hisano J, Murayama H, Yanagida T. Nucleon decay in the minimal supersymmetric SU(5) grand unification. *Nucl Phys*. (1993) **B402**:46–84. doi: 10.1016/0550-3213(93)90636-4

526. Lucas V, Raby S. Nucleon decay in a realistic SO(10) SUSY GUT. *Phys Rev*. (1997) **D55**:6986–7009. doi: 10.1103/PhysRevD.55.6986

527. Goto T, Nihei T. Effect of RRRR dimension five operator on the proton decay in the minimal SU(5) SUGRA GUT model. *Phys Rev*. (1999) **D59**:115009. doi: 10.1103/PhysRevD.59.115009

528. Murayama H, Pierce A. Not even decoupling can save minimal supersymmetric SU(5). *Phys Rev*. (2002) **D65**:055009. doi: 10.1103/PhysRevD.65.055009

529. Dermisek R, Mafi A, Raby S. SUSY GUTs under siege: proton decay. *Phys Rev*. (2001) **D63**:035001. doi: 10.1103/PhysRevD.63.035001

530. Bajc B, Fileviez Perez P, Senjanovic G. Proton decay in minimal supersymmetric SU(5). *Phys Rev*. (2002) **D66**:075005. doi: 10.1103/PhysRevD.66.075005

531. Bajc B, Fileviez Perez P, Senjanovic G. Minimal supersymmetric SU(5) theory and proton decay: where do we stand? In: *Beyond the Desert: Accelerator, Non-accelerator and Space Approaches Into the Next Millennium. Proceedings, 3rd International Conference on particle physics beyond the standard model, Oulu, Finland, June 2-7, 2002*), Oulu (2002). p. 131–9.

532. Emmanuel-Costa D, Wiesenfeldt S. Proton decay in a consistent supersymmetric SU(5) GUT model. *Nucl Phys*. (2003) **B661**:62–82. doi: 10.1016/S0550-3213(03)00301-8

533. Dorsner I, Fileviez Perez P. How long could we live? *Phys Lett*. (2005) **B625**:88–95. doi: 10.1016/j.physletb.2005.08.039

534. Fileviez Perez P. Fermion mixings versus d = 6 proton decay. *Phys Lett*. (2004) **B595**:476–83. doi: 10.1016/j.physletb.2004.06.061

535. Nath P, Fileviez Perez P. Proton stability in grand unified theories, in strings and in branes. *Phys Rept*. (2007) **441**:191–317. doi: 10.1016/j.physrep.2007.02.010

536. Abe K, Hayato Y, Iyogi K, Kameda J, Miura M, Moriyama S, et al. Search for proton decay via *p*→*vvK*^{+} using 260 kiloton year data of Super-Kamiokande. *Phys Rev*. (2014) **D90**:072005. doi: 10.1103/PhysRevD.90.072005

537. Sussman S, Abe K, Bronner C, Hayato Y, Ikeda M, Iyogi K, et al. Dinucleon and nucleon decay to two-body final states with no hadrons in super-kamiokande. *arXiv:1811.12430* (2018).

538. Abe K, Abe K, Aihara H, Aimi A, Akutsu R, Andreopoulos C, et al. Hyper-kamiokande design report. *arXiv:1805.04163* (2018).

539. Acciarri R, Acero MA, Adamowski M, Adams C, Adamson P, Adhikari S, et al. Long-baseline neutrino facility (LBNF) and deep underground neutrino experiment (DUNE). *arXiv:1512.06148* (2015).

540. Kolešová H, Malinský M. Flavour structure of GUTs and uncertainties in proton lifetime estimates. (2016) *Phys Rev*. **D99**:035005. doi: 10.1103/PhysRevD.99.035005

541. Leurer M. A Comprehensive study of leptoquark bounds. *Phys Rev*. (1994) **D49**:333–42. doi: 10.1103/PhysRevD.49.333

542. Di Luzio L, Kirk M, Lenz A. Updated *B*_{s}-mixing constraints on new physics models for *b*→*sℓ*^{+}ℓ^{−} anomalies. *Phys Rev*. (2018) **D97**:095035. doi: 10.1103/PhysRevD.97.095035

543. Lees JP, Poireau V, Tisserand V, Wu SL, Grauges E, Palano A, et al. Measurement of an excess of $\overline{B}\to {D}^{(*)}{\tau}^{-}{\overline{\nu}}_{\tau}$ decays and implications for charged higgs bosons. *Phys Rev*. (2013) **D88**:072012. doi: 10.1103/PhysRevD.88.072012

544. Hirose S, Iijima T, Adachi I, Adamczyk K, Aihara H, Said SA, et al. Measurement of the τ lepton polarization and *R*(*D*^{*}) in the decay $\overline{B}\to {D}^{*}{\tau}^{-}{\overline{\nu}}_{\tau}$. *Phys Rev Lett*. (2017) **118**:211801. doi: 10.1103/PhysRevLett.118.211801

545. Bailey JA, Bazavov A, Bernard C, Bouchard CM, DeTar C, Du D, et al. B → Dℓv form factors at nonzero recoil and |V_{cb}| from 2+1-flavor lattice QCD. *Phys Rev*. (2015) **D92**:034506. doi: 10.1103/PhysRevD.92.034506

546. Fajfer S, Kamenik JF, Nisandzic I. On the $B\to {D}^{*}\tau {\overline{\nu}}_{\tau}$ Sensitivity to New Physics. *Phys Rev*. (2012) **D85**:094025. doi: 10.1103/PhysRevD.85.094025

547. Bernlochner FU, Ligeti Z, Papucci M, Robinson DJ. Combined analysis of semileptonic *B* decays to *D* and *D*^{*}: *R*(*D*^{(*)}), |*V*_{cb}|, and new physics. *Phys Rev*. (2017) **D95**:115008. doi: 10.1103/PhysRevD.95.115008

548. Bigi D, Gambino P, Schacht S. A fresh look at the determination of |*V*_{cb}| from *B→D ^{*}ℓν. Phys Lett*. (2017)

**B769**:441–5. doi: 10.1016/j.physletb.2017.04.022

549. Amhis Y, Banerjee S, Ben-Haim E, Bernlochner F, Bozek A, Bozzi C, et al. Averages of *b*-hadron, *c*-hadron, and τ-lepton properties as of summer 2016. *Eur Phys J.* (2016) **C77**:895. doi: 10.1140/epjc/s10052-017-5058-4

550. Aaij R, Adeva B, Adinolfi M, Ajaltouni Z, Akar S, Albrecht J, et al. Test of lepton universality with *B*^{0}→*K*^{*0}ℓ^{+}ℓ^{−} decays. *J High Energy Phys*. (2017) **08**:055. doi: 10.1007/JHEP08(2017)055

551. Hiller G, Kruger F. More model-independent analysis of *b*→*s* processes. *Phys Rev*. (2004) **D69**:074020. doi: 10.1103/PhysRevD.69.074020

552. Bordone M, Isidori G, Pattori A. On the Standard Model predictions for *R*_{K} and ${R}_{{K}^{*}}$. *Eur Phys J*. (2016) **C76**:440. doi: 10.1140/epjc/s10052-016-4274-7

553. D'Amico G, Nardecchia M, Panci P, Sannino F, Strumia A, Torre R, et al. Flavour anomalies after the ${R}_{{K}^{*}}$ measurement. *J High Energy Phys*. (2017) **09**:010. doi: 10.1007/JHEP09(2017)010

554. Capdevila B, Crivellin A, Descotes-Genon S, Matias J, Virto J. Patterns of New Physics in *b*→*sℓ*^{+}ℓ^{−} transitions in the light of recent data. *J High Energy Phys*. (2018) **01**:093. doi: 10.1007/JHEP01(2018)093

555. Bhattacharya B, Datta A, London D, Shivashankara S. Simultaneous Explanation of the *R*_{K} and *R*(*D*^{(*)}) Puzzles. *Phys Lett*. 2015) **B742**:370–4. doi: 10.1016/j.physletb.2015.02.011

556. Freytsis M, Ligeti Z, Ruderman JT. Flavor models for $\overline{B}\to {D}^{(*)}\tau \overline{\nu}$. *Phys Rev*. (2015) **D92**:054018. doi: 10.1103/PhysRevD.92.054018

557. Buttazzo D, Greljo A, Isidori G, Marzocca D. B-physics anomalies: a guide to combined explanations. *J High Energy Phys*. (2017) **11**:044. doi: 10.1007/JHEP11(2017)044

558. Angelescu A, Becirevic D, Faroughy DA, Sumensari O. Closing the window on single leptoquark solutions to the *B*-physics anomalies. *J High Energy Phys*. (2018) **10**:183. doi: 10.1007/JHEP10(2018)183

559. Bauer M, Neubert M. Minimal Leptoquark Explanation for the ${\text{R}}_{{D}^{(*)}}$ , R_{K} , and (*g*−2)_{g} Anomalies. *Phys Rev Lett*. (2016) **116**:141802. doi: 10.1103/PhysRevLett.116.141802

560. Becirevic D, Sumensari O. A leptoquark model to accommodate ${R}_{K}^{\text{exp}}<{R}_{K}^{\text{SM}}$ and ${R}_{{K}^{*}}^{\text{exp}}<{R}_{{K}^{*}}^{\text{SM}}$. *J High Energy Phys*. (2017) **08**:104. doi: 10.1007/JHEP08(2017)104

561. Becirevic D, Fajfer S, Košnik N, Sumensari O. Leptoquark model to explain the *B*-physics anomalies, *R*_{K} and *R*_{D}. *Phys Rev*. (2016) **D94**:115021. doi: 10.1103/PhysRevD.94.115021

562. Becirevic D, Košnik N, Sumensari O, Zukanovich Funchal R. Palatable Leptoquark Scenarios for Lepton Flavor Violation in Exclusive *b*→*sℓ*_{1}ℓ_{2} modes. *J High Energy Phys*. (2016) **11**:035. doi: 10.1007/JHEP11(2016)035

563. Chen CH, Nomura T, Okada H. Excesses of muon *g*−2, ${R}_{{D}^{(*)}}$, and *R*_{K} in a leptoquark model. *Phys Lett*. (2017) **B774**:456–64. doi: 10.1016/j.physletb.2017.10.005

564. de Medeiros Varzielas I, King SF. ${R}_{{K}^{(*)}}$ with leptoquarks and the origin of Yukawa couplings. *J High Energy Phys*. (2018) **11**:100. doi: 10.1007/JHEP11(2018)100