Leptoquark mechanism of neutrino masses within the grand unification framework

We demonstrate the viability of the one-loop neutrino mass mechanism within the framework of grand unification when the loop particles comprise scalar leptoquarks (LQs) and quarks of the matching electric charge. This mechanism can be implemented in both supersymmetric and non-supersymmetric models and requires the presence of at least one LQ pair. The appropriate pairs for the neutrino mass generation via the up-type and down-type quark loops are S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_3$$\end{document}–R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_2$$\end{document} and S1,3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_{1,\,3}$$\end{document}–R~2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{R}_2$$\end{document}, respectively. We consider two distinct regimes for the LQ masses in our analysis. The first regime calls for very heavy LQs in the loop. It can be naturally realized with the S1,3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_{1,\,3}$$\end{document}–R~2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{R}_2$$\end{document} scenarios when the LQ masses are roughly between 1012\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{12}$$\end{document} and 5×1013\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$5 \times 10^{13}$$\end{document} GeV. These lower and upper bounds originate from experimental limits on partial proton decay lifetimes and perturbativity constraints, respectively. Second regime corresponds to the collider accessible LQs in the neutrino mass loop. That option is viable for the S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_3$$\end{document}–R~2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{R}_2$$\end{document} scenario in the models of unification that we discuss. If one furthermore assumes the presence of the type II see-saw mechanism there is an additional contribution from the S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_3$$\end{document}–R2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_2$$\end{document} scenario that needs to be taken into account beside the type II see-saw contribution itself. We provide a complete list of renormalizable operators that yield necessary mixing of all aforementioned LQ pairs using the language of SU(5). We furthermore discuss several possible embeddings of this mechanism in SU(5) and SO(10) gauge groups.


Introduction
Leptoquarks (LQs) are colored states that couple quarks to leptons. They can thus yield novel physical processes such as proton decay or help explain experimentally observed phenomena that cannot be successfully addressed within the Standard Model (SM) of elementary particle physics. For example, neutrino masses of Majorana nature can be generated through the one-loop level processes if one introduces at least two particular scalar LQ multiplets [1,2] to the SM particle content. It is our intention to investigate the viability of this particular mechanism within a context of grand unification. This is where the LQs first emerged after all [3][4][5][6]. For exhaustive lists of references on the LQ phenomenology one can consult reviews on the subject [7][8][9][10] or turn to the numerous studies of specific aspects of the LQ related physics [11][12][13][14][15][16][17]. The one-loop contributions towards neutrino masses that we study have been considered extensively in the literature [1,2,[18][19][20][21][22]. Our intention, in contrast to the existing studies, is to analyse possibilities to have a more fundamental origin of this mechanism and to provide several realistic examples.
The idea to have radiatively induced neutrino masses in the grand unified theory framework has been around for a very long time [23]. There are several explicit implementations of this approach that one can find in the literature within both the SU (5) [24][25][26][27][28][29][30] and the SO (10) [23,31,32] contexts. What sets our study apart from the existing work is that we exclusively use scalar LQs to generate neutrino masses at the one-loop level.
To start, we present an overview of the most salient features of the LQ neutrino mass mechanism. Only then do we proceed to discuss two distinct implementations of this approach to address the issue of neutrino mass within the grand unification framework. We list the transformation Table 1 Transformation properties of scalar LQs under the SM gauge group. The list of the most relevant SU (5) (SO (10)) representations that accommodate them is presented in the third (fourth) column. We assume the standard embedding of U (1) charges in SO (10) LQ (SU (3), SU (2), U (1)) SU (5) SO (10)  properties of scalar LQs under the SM gauge group in Table 1.
We adopt symbolic notation to represent LQ multiplets [14]. We also denote a given representation with the associated dimensionality whenever possible. To single out a particular electric charge eigenstate from a given LQ multiplet we use superscripts [10]. For example, S 3 comprises three electric charge eigenstates that we label S . This fixes the hypercharge normalization we use throughout the manuscript.
The mechanism we want to study, in its minimal form, requires the presence of one scalar multiplet that transforms asR 2 and another one that has the transformation properties of either S 1 or S 3 in addition to the SM particle content. The following two features are crucial if one is to generate neutrino mass(es) at the one-loop level. Firstly,R whereỹ RL 2 , y L L 1 , y L L 3 , and y D are 3 × 3 matrices in flavor space. 1 H (≡ (1, 2, 1/2)) is the Higgs boson of the SM, τ k , k = 1, 2, 3, are Pauli matrices, and a, b, c = 1, 2 are the SU (2) group space indices. The couplings ofR Secondly,R 2 can mix with either S 1 or S 3 through the Higgs boson. In fact, the LQ pairs S 1 -R should mix in order for the mechanism to work. In the latter case the states S −2/3 3 andR 2/3 * 2 also mix. The relevant parts of the scalar interactions are 1 The chiralities of the quark-lepton pair that the LQ couples to are denoted with the superscript labels ofỹ RL 2 , y L L 1 , and y L L 3 .
where λ 1 and λ 3 are dimensionful parameters that we take to be real for simplicity. We denote the squared-masses of the two physical LQs of the 1/3 electric charge with m 2 LQ 1 and m 2 LQ 2 regardless of whether these states originate from the S 1 -R −1/3 * 2 or S 1/3 combination. The angle that diagonalizes the 2 × 2 squared-mass matrix m 2 1 (m 2 3 ) for the S 1 -R ) pair is labeled θ 1 (θ 3 ). The squared-mass matrices m 2 1 and m 2 3 take the form where H represents a vacuum expectation value (VEV) of electrically neutral component of the SM Higgs field. Here, The mechanism is very economical since the same scalar field H , upon the electroweak symmetry breaking, provides masses for the SM charged fermions and introduces a mixing term for the LQs. The particles that propagate in the loop that generates neutrino Majorana mass(es) are the down-type quarks and scalar LQs of the matching electric charge. The associated one-loop Feynman diagrams are presented in the upper panel of Fig. 1.
The effective neutrino mass matrix in the basis of the physical down-type quarks and LQs reads [18] Fig. 1 The one-loop neutrino mass diagrams for the S 1, 3 -R 2 (S 3 -R 2 ) scenario in the upper (lower) panel. See text for full details (m N ) αβ = 3 sin 2θ 1, 3 32π 2 where α, β, δ = 1, 2, 3 are flavor indices, 33 ) are the down-type quark masses, and Before we proceed we have one specific comment with regard to the previous discussion. It concerns a possibility that the fermions that propagate in the neutrino mass loop are the up-type quarks instead of the down-type quarks. This seems to be a viable possibility if one starts with the R 2 -S 3 combination. The most essential Yukawa interactions for this scenario are where y RL 2 and y U are 3 × 3 matrices in flavor space. The couplings of R  but only if all three multiplets, i.e., R 2 ,R 2 , and S 3 , are present in the set-up [18]. Third option is to have one additional scalar S(≡ (1, 3, 1)) that acquires a VEV. The tree-level mixing of R 2/3 2 with S −2/3 * 3 is then possible and the off-diagonal entry of the relevant 2 × 2 squared-mass matrix is proportional to a product of the VEVs of neutral fields in S and H . The scalar interactions that are needed to implement the second and third option are where λ 2 is a dimensionful parameter, whereas κ 1 and κ 2 are both dimensionless parameters. One can trivially adapt Eqs. (4), (5), and (6) to the up-type quark scenario in order to find the associated neutrino mass matrix m N . Let us denote with θ 2 the mixing angle between where S represents the VEV of electrically neutral component of S and κ = κ 1 + κ 2 . All one needs to do is to first evaluate θ 2 by replacing λ 1, 3 H with 2κ H S in Eq. (4) and then substitute parameters θ 1, 3 ,ỹ RL 2 , and y L L 1, 3 with θ 2 , y RL 2 , and − √ 2(V * CKM y L L 3 ), respectively, in Eqs. (5) and (6). The down-type quark masses in Eq. (6) also need to be replaced with the masses of the up-type quarks, i.e., (m 1 , m 2 , pair is shown in the lower panel of Fig. 1. We will make further comments on this potentially important contribution towards neutrino masses later on. Our aim is to implement the one-loop neutrino mass mechanism in the framework of grand unification. We accordingly investigate viability of two distinct regimes in Sect. 2 using mainly the language of SU (5) gauge group. First regime corresponds to a scenario where the LQs behind the neutrino mass generation reside at a very high energy scale. This possibility is discussed in Sect. 2.1. Second regime corresponds to a scenario where the neutrino masses are generated with the Large Hadron Collider (LHC) accessible scalar LQs. We demonstrate viability of that scenario in Sect. 2.2. The summary of our findings is presented in Sect. 3.

Grand unification vs. one-loop neutrino mass
Let us proceed with a realistic implementation of the oneloop neutrino mass mechanism with scalar LQs in the grand unification framework. We primarily use the language of the SU (5) gauge group in what follows. The SM fermions reside in 10 α and 5 α of SU (5), where α(= 1, 2, 3) is a flavor index [6]. The exact decompositions of 5 α (10 α ) under the SM reads The relevant embeddings of scalar LQs in the SU (5) representations are presented in Table 1. We clearly need to have either one 10-or one 15-dimensional scalar representation in order to introduce oneR 2 multiplet in any SU (5) model. Relevant contraction that yields y RL 2d R ν LR −1/3 2 term whenR 2 is part of 10-dimensional (15dimensional) representation is y αβ 5 α 5 β 10 (y αβ 5 α 5 β 15). We identify (ỹ RL 2 ) αβ to be −y αβ / √ 2, where y αβ are elements of an antisymmetric (symmetric) complex matrix in the case whenR 2 originates from 10-dimensional (15-dimensional) representation.
The mass mechanism that we discuss can also be implemented in the SO(10) framework. See Table 1 for the standard embedding of scalar LQs in the SO (10) representations. In particular, ifR 2 originates from 120-dimensional (126-dimensional) representation of SO(10) the relevant couplings to the SM fermions will be antisymmetric (symmetric) in flavor space. These properties thus closely mirror the SU (5) flavor structure of theR 2 couplings. The associated SO (10) operators are y αβ 16 α 16 β 120 and y αβ 16 α 16 β 126, where we assume that one 16-dimensional SO (10) representation comprises one generation of the SM fermions and one right-chiral neutrino.
The origin of the term y L L 3d C L ν L S 1/3 3 is unique in SU (5) as can be seen from Table 1. Namely, S 3 resides in a 45dimensional representation and the relevant contraction that generates aforementioned couplings is y 45 αβ 10 α 5 β 45. One can thus identify y L L 3 with y 45 / √ 2, where y 45 is related to the masses of the charged fermions and down-type quarks as we show in the next paragraph. The situation with R 2 seems more involved since R 2 can be found in either 45- where all the VEVs are taken to be real. The VEV normalization yields v 2 5 /2 + 12v 2 45 = v 2 , where v(= 246 GeV) is the electroweak VEV [34]. The SU (5) symmetry thus dictates that y 45 45 ). The term y L L 1d C L ν L S 1 originates from y 5 αβ 10 α 5 β 5 and y 45 αβ 10 α 5 β 45 for S 1 ∈ 5 and S 1 ∈ 45, respectively. In the former (latter) case one can identify y L L 1 with −y 5 / √ 2 (y 45 /2). Finally, one needs to provide the mixing term for at least one of the relevant LQ pairs in order to complete the neutrino mass loop. There are two very different regimes for the scalar LQ masses that we can envisage with this in mind. First option is that the LQs behind the neutrino mass generation reside at a very high energy scale. This could provide compliance of the set-up with the experimental bounds on proton decay. The main issue with this regime could be associated with the size of the relevant lepton-quark-LQ couplings. Namely, these couplings might need to be unrealistically large in order to (re)produce neutrino mass scales that are compatible with experimental observations. It turns out that this is not the case and we accordingly demonstrate in Sect. 2.1 why and how this particular scenario can be realized within the grand unification frameworks.
The second option is that the scalar LQs are very light. That scenario is especially appealing since the LHC accessible LQs could also affect flavor physics observables. The main difficulty with this particular set-up is to explain observed levels of matter stability. 2 Namely, S 1 and S 3 can both have "diquark" couplings that, in combination with the lepton-quark-LQ couplings that are needed to generate neutrino masses, destabilize protons and bound neutrons. 3 To avoid conflict with stringent limits on proton lifetime one would need to either forbid or substantially suppress these "diquark" operators. This might be very difficult from the model building point of view since unification of matter multiplets dictates common origin of both types of couplings. One would also need to prevent mixing between these LQs and any other LQ in the theory that has "diquark" couplings to ensure stability of matter. This might also represent a challenge since one needs to mix specific LQ multiplets in order to generate neutrino masses in the first place. We show that both of these issues can be successfully addressed for the S 3 -R 2 and S 3 -R 2 scenarios in Sect. 2.2. The S 1 -R 2 option, on the other hand, is problematic due to difficulty with suppression of the S 1 "diquark" couplings in the simplest of models and we opt not to discuss it in the light LQ regime.

Heavy leptoquark regime
Let us turn our attention to a scenario where the LQs are heavy. We assume in what follows that all the LQ masses need to be at or exceed 10 12 GeV to ensure proton stability. This is a very conservative estimate since it is certainly above a lower bound that can be extracted from the latest data on proton stability within the SU (5) framework [36]. We show that the one-loop neutrino masses can be realized in this part of phenomenologically available parameter space if the fermions in the neutrino mass loop are exclusively the down-type quarks.   16 GeV. The necessary mixing between S 1 (∈ 5) andR 2 (∈ 15) can be generated through contractions of the form 5 i 5 j 15 i j and 5 i 5 j 15 jk 24 i k . These, again, yield an angle θ 1 that is comparable in strength to our estimate for θ 3 . We can furthermore safely assume that the m b (≈ 1 GeV) contribution dominates the sum in Eq. (5). Putting all this together implies that where we suppress flavor indices and assume that the mass splitting between LQs is not substantial, i.e., we take that ln(m 2 LQ 2 /m 2 LQ 1 ) ≈ 1. The approximation of Eq. (13) shows that the entries in the product (ỹ RL 2 y L L 1, 3 ) do not have to be very large to correctly describe the neutrino mass scale. For example, in the non-degenerate normal hierarchy case of the neutrino masses the largest entry on the left side of Eq. (13) needs to be at the level of 5 × 10 −2 eV which would imply that (ỹ RL 2 y L L 1, 3 ) ∼ 5 × 10 −3 . 4 The back-of-the-envelope estimate we present clearly demonstrates the viability of this option. Note that there is an upper bound on the heavier of the two LQs in this set-up if one demands perturbativity of the Yukawa coupling entries inỹ RL 2 and y L L 1, 3 matrices. We find it to be roughly at 5 × 10 13 GeV. This implies that the two LQs must reside in relatively narrow mass window from 10 12 to 5 × 10 13 GeV in order to accommodate all the relevant constraints. One can furthermore infer that ln(m 2 LQ 2 /m 2 LQ 1 ) < 5, which is in agreement with our initial assumption.
This particular possibility to generate neutrino masses, in our view, has been overlooked in the literature on grand unification. For example, there are two non-supersymmetric models that already have all the necessary ingredients to incorporate this particular scenario. The first model [30] introduces one 10-dimensional scalar representation on top of 5, 24, and 45 in order to generate neutrino masses through the Zee mechanism [24]. The second model [38] resorts to one 15dimensional scalar representation in addition to 5, 24, and 45 in order to generate neutrino masses through the type II see-saw mechanism [39,40]. Again, both of these models can accommodate the one-loop mechanism we discuss.
The heavy LQ regime is also tailor-made for the SO(10) framework. This could especially be beneficial in the scenarios that fail to accommodate neutrino masses in satisfactory manner. Clearly, it is sufficient to have either a 120-or a 126dimensional representation to introduce LQs that transform as S 1 , S 3 , andR 2 . This means that the relevant LQ couplings to the SM matter are always in place if one assumes standard embedding of the SM fermions in SO (10). The only remaining element, i.e., the LQ mixing, depends on the exact scalar sector of the SO(10) theory. We opt to show only one example due to existence of several distinct ways one can realistically break SO (10) For example, if we introduce one 210-dimensional representation to break SO (10) there is an operator of the form 210 10 126 that exists regardless of whether the theory is supersymmetric or not that yields a mixing between S 1 (∈ 10) andR 2 (∈ 126), Here, 10 and 126 are scalar representation that generate masses of the SM charged fermions.

Light leptoquark regime
To demonstrate that the collider accessible LQ scenario is a viable option to generate neutrino masses one needs to address the issue of the LQ mixing. Namely, if the genuine LQ states mix with the states that have "diquark" couplings it is hard to imagine that matter stability holds at the experimentally observed levels. We focus exclusively on a scenario whenR 2 originates from 15-dimensional representation. The analysis for the 10-dimensional representation case is completely analogous as we show in Appendix A. The SU (5) scenario under consideration comprises the following scalar representations: 5, 15, 24, and 45. We note that R 2 ,R 2 , and S 3 do not have "diquark" couplings [41] at renormalizable level if the charged fermion mass relations are given with Eqs. There are ten non-trivial operators that mix the LQ states of the same electric charge if the only VEVs present are the ones proportional to v 24 , v 45 , and v 5 . Nine (four) of these contractions affect the 1/3 (2/3) electric charge states. There are no contractions that mix LQs of the 4/3 electric charge through these VEVs. The complete list of relevant SU (5) contractions is relegated to Appendix A. It turns out that one can write a 4 × 4 squared-mass matrix for the 1/3 electric charge LQs in a block diagonal form where the relevant two blocks are of dimension 2 × 2 each. The basis for this matrix is , where we explicitly denote the origin of LQ multiplets that transform as S 1 under the SM gauge group. The mixing term between S 1/3 * 3 andR −1/3 2 we referred to previously as λ 3 H is proportional to a product of v 24 with v 5 . Since the LQs of the 4/3 electric charge do not mix the associated 2 × 2 squared-mass matrix has only diagonal entries. These findings guarantee the matter stability even in the presence of the mixing that is needed to generate neutrino masses. Components of S 3 andR 2 can thus be very light and the resulting neutrino mass matrix is correctly described through the expression of Eq. (5) due to a block diagonal form of the relevant LQ squared-mass matrix. We briefly postpone the discussion of the mixing between the LQ states with electric charge of 2/3, since these originate from R 2 ,R 2 , and S 3 multiplets that have no "diquark" couplings in this set-up and consequently do not directly affect matter stability.
Let us summarize the main features of the light LQ set-up.R 2 (S 3 ) originates from 15 (45) of SU (5). Again, R 2 could instead originate from 10-dimensional representation. The SU (5) symmetry is broken by the VEV of 24 down to SU (3) × SU (2) × U (1). The Higgs field VEVs that complete the electroweak symmetry breaking reside in both 5 and 45. The light LQ states are components of R 2 and S 3 and they help generate neutrino masses. Three out of six LQs of the model -S 1 (45), S 1 (5), andS 1mediate proton decay and need to be heavy. R 2 can in principle be of an arbitrary mass. Finally, the mass matrix for the up-type quarks is symmetric in accordance with Eq. (12) which has implications for the gauge-mediated proton decay [42]. We plan to pursue the phenomenology of this set-up in future work. In this respect, the state S 3 with mass close to the LHC reach has been proven to play a beneficial role in addressing hints of lepton flavor universality violation in b → s and b → c ν processes [43,44].
We have, in our analysis, neglected possible VEVs of electrically neutral fields in 15-and 24-dimensional representations. The former (latter) field resides in the (1, 3, 1) ((1, 3, 0)) component of 15 (24). We normalize these additional VEVs of 15 ≡ 15 i j and 24 ≡ 24 i j to be 15 55 = v 15 and 24 4 4 = − 24 5 5 = v S , respectively. The presence of these VEVs introduces seven additional SU (5) operators that one needs to include in the analysis of the LQ mixing. We list these operators in Appendix A.
The one-loop mechanism we discuss is not the only possible contribution towards neutrino masses in the light LQ regime. Note that the VEV of the 15-dimensional representation can generate neutrino mass(es) of Majorana nature through the type II see-saw mechanism [39,40]. 5 More importantly, the up-type quarks can also contribute towards neutrino mass generation since the scalars R 2/3 2 , R 2/3 2 , and S −2/3 * 3 mix with or without the VEV of the 15dimensional representation [18]. In the latter case we find that the up-type quark contribution is completely negligible. In the former case the mixing angle θ 2 between R 2 ) that the relevant off-diagonal entries 12, 13, and 23 for the symmetric squared-mass matrix of the 2/3 electric charge LQs are proportional to v 15 v 5 , v 24 v 5 , and v 45 v 5 , respectively. This can increase the maximum allowed value of v 15 but only by a factor of 10. The leading neutrino mass contributions due to propagation of the up-type quarks and the down-type quarks are thus proportional to O(10 −3 )m t and O(1)m b , respectively, and can be comparable in strength in some parts of the available parameter space. A self-consistent study of the neutrino mass(es) should take into account all these contributions ifR 2 originates from 15-dimensional representation and the VEV proportional to v 15 is turned on. IfR 2 originates from the 10-dimensional representation the only relevant contribution in this regime is due to the down-type quark loop.

Conclusions
The one-loop neutrino mass mechanism with scalar LQs in the loop can be embedded within the framework of grand unification, regardless of whether the scenario is supersymmetric or not. There exist two distinct regimes for the LQ masses.
One option is to have heavy LQs in the loops that generate neutrino masses. This option can be naturally realized with the S 1, 3 -R 2 combinations of LQs. The type II see-saw mechanism contribution could also be present and important in some parts of the accessible parameter space. The nice feature of the heavy LQ limit is that the masses of the LQs in the loop can only be between 10 12 and 5 × 10 13 GeV in order to simultaneously avoid experimental limits on partial proton decay lifetimes and still satisfy perturbativity constraints on the lepton-quark-LQ couplings.
The other option is to have collider accessible LQs in the loop. That particular limit can be realized via the loops that contain the down-type quarks and scalars of the matching electric charge that are the mixture of S 3 andR 2 multiplets. The S 1 -R 2 combination is not a viable option in this limit due to existence of "diquark" couplings of S 1 in the minimal set-up. If the theory also contains an SU (2) triplet scalar (1, 3, 1) that gets the VEV one needs to take into account two additional neutrino mass contributions. One is the type 6 Note that v 15 is bounded from above due to the existing electroweak precision measurements of the so-called ρ parameter [47]. This bound can be avoided if one judiciously adjusts v 15 and v S to be approximately equal [38]. 7 For the latest direct bounds on LQ masses from the LHC data see, for example, Refs. [48,49]. II see-saw contribution and the other one is the one-loop contribution due to propagation of the up-type quarks and the scalar states of the same electric charge that originate from the mixture of S 3 and R 2 multiplets. These three mechanisms can coexist and be of equal importance in some parts of the available parameter space.
We discuss possible origins of scalar LQs that are needed to complete the neutrino mass generating loops using the language of SU (5). We also provide a list of all SU (5) contractions that generate the LQ mixing terms. We furthermore argue that all the necessary ingredients to implement the oneloop neutrino mass mechanism are present in any SO (10) theory with the standard embedding of the matter fields that generates charged fermion masses through renormalizable contractions.
To obtain a scenario comprising 5-, 10-, 24-, and 45dimensional scalar representations one should replace 15dimensional representation with 10-dimensional one wherever possible. Note that some of the contractions that one obtains with the simple substitution yield zero due to the antisymmetric nature of 10 i j (= −10 ji ) in the SU (