KCL-PH-TH/2015-42, LCTS/2015-29, CERN-PH-TH/2015-229,

UMN–TH–3502/15, FTPI–MINN–15/41, IPMU15–0174, CETUP2015-021

Beyond the CMSSM without an Accelerator:

Proton Decay and Direct Dark Matter Detection

John Ellis, Jason L. Evans, Feng Luo, Natsumi Nagata,

Keith A. Olive and Pearl Sandick

Theoretical Physics and Cosmology Group, Department of Physics,

King’s College London, Strand, London WC2R 2LS, UK

TH Division, Physics Department, CERN, CH-1211 Geneva 23, Switzerland

William I. Fine Theoretical Physics Institute, School of Physics and Astronomy,

University of Minnesota, Minneapolis, MN 55455, USA

Kavli IPMU (WPI), UTIAS, University of Tokyo, Kashiwa, Chiba 277-8583, Japan

Department of Physics and Astronomy, University of Utah, Salt Lake City, Utah 84112, USA

Abstract

We consider two potential non-accelerator signatures of generalizations of the well-studied constrained minimal supersymmetric standard model (CMSSM). In one generalization, the universality constraints on soft supersymmetry-breaking parameters are applied at some input scale below the grand unification (GUT) scale , a scenario referred to as ‘sub-GUT’. The other generalization we consider is to retain GUT-scale universality for the squark and slepton masses, but to relax universality for the soft supersymmetry-breaking contributions to the masses of the Higgs doublets. As with other CMSSM-like models, the measured Higgs mass requires supersymmetric particle masses near or beyond the TeV scale. Because of these rather heavy sparticle masses, the embedding of these CMSSM-like models in a minimal SU(5) model of grand unification can yield a proton lifetime consistent with current experimental limits, and may be accessible in existing and future proton decay experiments. Another possible signature of these CMSSM-like models is direct detection of supersymmetric dark matter. The direct dark matter scattering rate is typically below the reach of the LUX-ZEPLIN (LZ) experiment if is close to , but may lie within its reach if GeV. Likewise, generalizing the CMSSM to allow non-universal supersymmetry-breaking contributions to the Higgs offers extensive possibilities for models within reach of the LZ experiment that have long proton lifetimes.

September 2015

## 1 Introduction

Supersymmetry remains a favored extension of the Standard Model (SM), despite its non-appearance during Run 1 of the LHC[1, 2]. Indeed, the discovery of a 125-GeV Higgs boson at the LHC [3] has supplemented the traditional arguments for supersymmetry, including the naturalness of the electroweak scale, the unification of the fundamental interactions and the existence of a cold dark matter candidate (if -parity is conserved). The minimal supersymmetric extension of the SM (MSSM) predicts the existence of a Higgs boson with mass GeV, and is a prime example of new physics capable of stabilizing the electroweak vacuum for GeV [4]. Moreover, global fits in the framework of simple supersymmetric models suggest that the couplings of the lightest supersymmetric Higgs boson should be very similar to those of the Higgs boson in the SM, as is indicated by the ATLAS and CMS experiments [5, 6]. When the supersymmetric particle masses are large, which is the case we consider, the Higgs couplings become even more like the SM couplings.

If these arguments are valid, the following questions must be answered: which supersymmetric model is found in Nature, and how may it be tested? To begin to answer these questions, we focus here on the MSSM, and more specifically on constrained versions in which the soft supersymmetry-breaking scalar masses , gaugino masses and trilinear terms are assumed to have universal values at some high input mass scale . Typically, is chosen to be at the grand-unified-theory (GUT) scale, a scenario called the constrained MSSM (CMSSM) [7, 8, 9, 10, 11, 12], in which the ratio of Higgs vacuum expectation values, , is a free parameter.

In order to find models with less-constrained dark matter scenarios and simultaneously a sufficiently long lifetime for the proton, we focus here on two one-parameter extensions of the CMSSM: ‘sub-GUT’ models [13] in which is free, and the NUHM1 [14, 15], in which the two Higgs soft masses are equal at the input scale, , but are allowed to differ from . We will also discuss ‘sub-GUT’ models obtained from minimal supergravity (mSUGRA), which are more constrained than the CMSSM, since the gravitino mass and the trilinear and bilinear soft supersymmetry-breaking terms are related: [16, 17]. Since mSUGRA models have one fewer parameter, is no longer free. Sub-GUT mSUGRA models have the same number of free parameters as in the CMSSM, but viable models can readily be found.

Although the standard CMSSM with is still viable, there remain only restricted regions of the parameter space of the CMSSM (and, a fortiori, of the more restrictive mSUGRA model) in which a successful prediction for can be reconciled with the measured cold dark matter density [11, 12, 18, 19, 20, 21, 22, 23, 24, 25, 26]. The parameter spaces of these models become more restricted when they are embedded in an SU(5) GUT, because they tend to have a proton lifetime which is shorter than the current experimental limits [27, 28, 29], even if the supersymmetric sparticle masses are rather heavy.

These problems can be avoided, however, if is identified with
some scale lower than the typical GUT scale. An effective scale
of supersymmetry breaking significantly below the GUT scale,
, is not without theoretical motivation. For example,
mirage unification models [30] and other scenarios such
as [31] give exactly such boundary conditions for the soft supersymmetry-breaking
parameters. Phenomenologically, these sub-GUT models with
have been shown to have an appropriate cold dark matter density
in a considerably larger parameter space [21]. As
could be expected, a significant part of this parameter space contains points
that are compatible with the LHC measurement of and other
phenomenological constraints, such as the non-detection of
supersymmetric particles at the LHC [1, 2] ^{1}^{1}1Note
that we do not impose any constraint from the anomalous magnetic moment
of the muon, [32].. This reduced tension in
sub-GUT models is due to the reduced running of the soft
masses, which leads to a sparticle spectrum that is, in
general, more compressed than models with . Moreover, this compression of the spectrum leads to more avenues for
coannihilation [33], which is effective in reducing the relic neutralino density into the range allowed by
cosmology.

Similarly, if the Higgs soft masses are allowed to differ from the soft masses of the matter scalars, as in the NUHM, there are more viable options for dark matter. In both the CMSSM and in the NUHM, the Higgs mixing mass, , and the pseudoscalar mass, , are determined through the minimization of the Higgs potential. However, either and/or can be traded for the Higgs soft mass, which can be calculated using the minimization of the Higgs potential.

Here we examine two potential non-accelerator observables in the contexts of these less-constrained models: the proton lifetime and the elastic scattering cross section for the direct detection of dark matter.

In [34], the proton lifetime was computed by renormalization group (RG) running the gauge couplings up to the GUT scale, defined to be where the two electroweak couplings are equal. The imperfection in the unification of the electroweak couplings with the strong coupling was then used to determine the size of the color-triplet Higgs threshold, which then determined the color-triplet Higgs mass [35, 36, 37]. Using this procedure with sub-TeV stops and Higgsinos and decoupled first- and second-generation sfermions, it was shown that the lifetime of the proton was shorter than the experimental constraints. Since the experimental constraints are now stronger, this problem has become even worse.

This problem can be avoided in many ways. One particularly simple way is to include an additional pair of and Higgs boson supersymmetric multiplets that do not couple to any of the SM fields. Below the SU(5) breaking scale, the colored and flavored Higgs mass become free parameters. If the portion of the Higgs supersymmetric multiplet that has SU(3) charges is lighter than the portion with SU(2) charges, the thresholds in the couplings will be different from those of minimal SU(5), and the colored Higgs masses can be made sufficiently heavy that the proton decay constraints can be met [34]. Other possibilities for alleviating this problem include forbidding the dimension-five operator leading to proton decay using extra dimensions [38], more complicated Higgs sectors [39], flipped SU(5) [40], or a Peccei-Quinn symmetry [41]. The problem is also alleviated in models with scalar masses that are TeV [42, 43, 44], as in pure gravity mediation [45].

The LHC constraints on sfermion masses and the observed Higgs mass of GeV motivate us to consider once again the decoupling limit as an explanation for the long lifetime of the proton. The decoupling limit was unsuccessful in [34] due to the assumption of a light third generation. However, if the third generation is also decoupled, the proton lifetime is extended. Since a heavier third generation is favoured by the GeV mass of the Higgs, we find this to be a reasonable scenario for suppressing proton decay. However, the real challenge in this scenario is to find regions of parameter space that combine a viable dark matter candidate with an acceptably long proton lifetime.

We will find that the minimal supersymmetric grand unified theory based on SU(5) [46] with a CMSSM spectrum does have a very small region of parameter space that combines a GeV Higgs, a sufficiently long proton lifetime, and a viable dark matter candidate. This may occur either in the focus-point region [47] or in a region where the dark matter density is obtained by stop coannihilation with the bino [48, 24]. However, as we show below, sub-GUT and NUHM1 models are less restricted by proton decay constraints. The proton lifetime is longer in sub-GUT models, in general, because the stop masses are larger due to reduced RG running. Since the lifetime of the proton scales as a power of the stop mass, this enhances the proton lifetime. Moreover, these sub-GUT models have an acceptable dark matter density in regions where the bino and the lighter stau coannihilate [49]. In the NUHM1, TeV Higgsinos are possible for any value of and and, if the Higgsino mass is (TeV), the Higgsino can be a thermal relic dark matter candidate [50]. If and are large, the proton lifetime is greatly enhanced because of the large stop mass. In all cases, compatibility with minimal SU(5) requires relatively low values of .

We also examine whether such models are compatible with present experimental constraints on the direct detection of dark matter through spin-independent elastic scattering, as provided, e.g., by the LUX experiment [51], and whether they can be probed by the next generation of such experiments, e.g., XENON1T [52] and LUX-ZEPLIN (LZ) [53].

For the purpose of our study, we use FeynHiggs [54] to calculate as a function of the model input parameters. Since one expects an uncertainty GeV in this calculation, we assume that any model yielding a prediction GeV may be acceptable. Even with this theoretical uncertainty, we find that the measurement generally gives a stronger constraint on the model parameters than do the direct LHC searches for supersymmetric particles published so far. As we also discuss, another important constraint is provided by the experimental search for decay [55], particularly at large . Since we do not impose any constraint, and the Higgs and other LHC constraints exclude small values of , the impact of the constraint [56] is reduced. We use SSARD [57] to calculate the particle spectrum, proton lifetimes, and elastic scattering cross sections.

The layout of this paper is as follows. In Section 2 we summarize the features of the CMSSM and mSUGRA models that are relevant for our analysis, and introduce their extensions to sub-GUT and NUHM1 models. In Section 3, we discuss the basics of our calculations of the proton life-time and elastic scattering cross sections in CMSSM-like models. Section 4 displays our results. Finally, Section 5 summarizes our conclusions.

## 2 The Models

### 2.1 mSUGRA

Minimal supergravity (mSUGRA) models have a quadratic Kähler potential for the chiral superfields, and the effective scalar potential is [58, 59, 16]

(1) |

where is the superpotential for the matter superfields ^{2}^{2}2We use the same notation
for the chiral superfields and their spin-zero components., which takes the following form in the MSSM:

(2) |

We denote by and (, , and ) the left- (right-)handed matter superfields, the Yukawa couplings are denoted by the , and are the pair of MSSM Higgs doublets with superpotential mixing coefficient . The in Eq. (1) are the trilinear terms in the superpotential, is the bilinear part corresponding to the term, and is the gravitino mass. In mSUGRA one finds scalar mass universality with , and there is a relation between the tri- and bilinear supersymmetry-breaking terms:

(3) |

These conditions apply at an input renormalization scale, , which may or may not be identified with the grand unification scale . If the gauge kinetic function is minimal, there is also gaugino mass universality, with a common mass that we assume to apply at the same input scale .

The two electroweak vacuum conditions are

(4) |

and

(5) |

where the soft supersymmetry-breaking Higgs masses denoted by are here evaluated at the weak scale, and and are loop corrections [60]. An mSUGRA model has just three continuous parameters: , and . The conditions (4, 5) can be used to determine as well as the magnitude of , but the sign of is undetermined. We consider in this paper both signs of in selected cases.

### 2.2 The CMSSM

The CMSSM is effectively a one-parameter generalization of mSUGRA, in which the relation (3) between and is dropped, which allows to be taken as an extra free parameter. In addition, in general, which is possible in SUGRA models only if the supergravity Kähler potential has non-minimal kinetic terms. Thus the CMSSM is specified by four parameters, , , , , the sign of . Here we assume that is sufficiently large to be irrelevant.

### 2.3 Sub-GUT versions of the CMSSM and mSUGRA

Generalizations of both mSUGRA and the CMSSM are possible if the input scale for universality of the supersymmetry-breaking terms differs from . We concentrate here on ‘sub-GUT’ models with [13]: a value of above the GUT scale [61, 62] would introduce many more GUT parameters, requiring a separate in-depth study. Sub-GUT versions of mSUGRA have four parameters: , and , whereas sub-GUT versions of the CMSSM have as an extra parameter (assuming again that is irrelevantly large). We found in [21] that sub-GUT mSUGRA models are phenomenologically viable in a relatively restricted range of that straddles the Polonyi value [63, 16].

### 2.4 The NUHM1

Another one- or two-parameter generalization of the CMSSM is the NUHM [64, 9], in which the values of the soft supersymmetry-breaking contributions to the Higgs masses and at the input scale are allowed to differ from the universal scalar mass . In the NUHM1 considered here, it is assumed that at the input scale [14, 15]. One may choose either or (through its relation to ) as a free parameter, and use the minimization conditions (4 and 5) to solve for . The examples shown here treat as a free parameter, since this displays more readily the interesting results in this scenario.

The two-parameter extension known as the NUHM2 [65, 15] drops the requirement that at the input scale. In this case, both and are allowed to be free input parameters. Alternatively, one can choose both and at the weak scale as free input parameters and use the minimization conditions (4 and 5) to solve for and . We do not study the NUHM2 in this paper, as most of the interesting aspects of the NUHM2 are contained in NUHM1 scans.

## 3 Calculations

### 3.1 Proton Decay Lifetimes

In this subsection we describe how we calculate proton decay rates in the minimal supersymmetric SU(5) GUT model [46]: for further discussions and detailed formulae, see [42, 44, 45, 66]. This model is the simplest supersymmetric extension of the original Georgi-Glashow model [67], in which the MSSM matter superfields are embedded into a representation of SU(5) for each generation, and the MSSM Higgs superfields and are incorporated in a pair of and 5 superfields, respectively. The SU(3) components of the and 5 are called the color-triplet Higgs fields. The dominant contribution to proton decay in this model is given by the exchange of these color-triplet Higgs fields [68], which induce dimension-five baryon-number violating operators, whereas the exchanges of SU(5) gauge bosons yield dimension-six operators. In this case, the dominant proton decay channel is the mode, and we focus on the partial decay rate for this channel in the following.

We obtain the low-energy effective Lagrangian below the GUT scale by first integrating out the color-triplet Higgs fields. In superfield notation, the effective Lagrangian is given by

(6) |

with the effective operators and defined by

(7) |

where are SU(3) color indices, are generation indices, and is the totally antisymmetric tensor. The Wilson coefficients of the above operators, and at the GUT scale are evaluated from the tree-level color-triplet Higgs exchange diagrams, with the results given in Appendix A. As shown in Eq. (21) of the appendix, these Wilson coefficients include up-type quark Yukawa couplings and down-type quark/lepton Yukawa couplings , which should be unified at the GUT scale in the minimal SU(5) GUT. Note, however, that although the minimal SU(5) GUT relation between the bottom and masses [69] is approximately consistent with the experimental values, this is not the case for the strange and masses, nor for the down and masses. The GUT Higgs couplings must therefore be more complicated than in the minimal SU(5) GUT, e.g., with additional higher-dimensional Higgs representations [70] and/or contributions to the fermion masses from higher-dimensional superpotential terms [71]. In practice, the ambiguity in choosing whether the down-type quark or lepton Yukawa couplings sets the scale of proton decay results in about a factor of 20 uncertainty in the proton decay calculation [45], which represents our ignorance of GUT-scale physics in the Yukawa sector. In the following calculation, we take down-type quark Yukawa couplings as the GUT-scale Yukawa couplings, which leads to longer proton decay lifetimes and thus gives rather conservative bounds on the model parameter space.

In addition, the GUT-scale Yukawa couplings introduce two extra phase factors [72], which give rise to additional uncertainty [66]. It turns out, however, that the effects of these unknown phases are actually negligible, as shown in Appendix A. Thus, we neglect these effects in our analysis.

After integrating out the color triplet Higgs boson, the GUT-scale Wilson coefficients are then evolved down to
the sfermion mass scale according to the renormalization group equations
(RGEs), which are also presented in
Appendix A. At the sfermion mass threshold,
sfermions in the external lines of the dimension-five effective
operators are integrated out via the one-loop diagrams illustrated in
Fig. 1, yielding
dimension-six four-fermion operators [73, 74]. In the absence of flavor violation in
the sfermion sector ^{3}^{3}3The effects on proton decay of possible flavor violation in the
sfermion sector are discussed in [44]., only the operators, and with , , and
give sizable contributions to proton decay; the contributions
of the other operators are suppressed by small Yukawa couplings and/or the
off-diagonal CKM matrix elements. The one-loop diagrams in
Fig. 1 yield the following effective Lagrangian below the
sfermion mass scale:

(8) |

with the operators composed of SM fermion fields,

(9) |

corresponding, respectively, to the operators and in [74], where are SU(2) indices, is the totally antisymmetric tensor, and , , and . The coefficients of the operators in (8) are also given in Appendix A. Note that, since a chirality flip in the internal wino/Higgsino propagator is required in the processes shown in Fig. 1, the operator coefficients contain factors of if , where is the wino or Higgsino mass and is the mass of a sfermion running in the loop. As a result, if the magnitude of the Higgsino mixing term is much smaller than that of the wino mass , the wino contribution dominates the Higgsino contribution, and vice versa if .

The coefficients in (8) are then run down to the electroweak scale using the one-loop RGEs given in [75]. At the electroweak scale, we transform to an operator basis in the low-energy theory, and the operator coefficients are evolved to the hadron scale, GeV, using the two-loop RGEs obtained in [76]. Finally, using the hadron matrix elements of the operators at GeV, we obtain the partial decay width of the channel. These matrix elements are evaluated using the QCD lattice simulation performed in [77]. This procedure, as well as the relevant formulae, is also summarized in Appendix A.

As can be seen from (21), the proton decay rate depends on the mass of the color-triplet Higgs field . Thus, to evaluate the proton lifetime, we need to determine the size of . To that end, we use the method discussed in [35, 36, 37]. In this method, the GUT-scale threshold corrections to the gauge coupling constants are used to estimate the masses of the GUT particles. The GUT-scale matching conditions for the gauge coupling constants at one-loop level in the scheme [78] in the minimal SU(5) GUT are given as follows [79, 80, 81], assuming no additional GUT-scale physics:

(10) |

where , , , and are the gauge coupling constants of U(1), SU(2), SU(3), and SU(5), respectively, with related to the hypercharge gauge coupling through , and and are the masses of the heavy gauge bosons and the adjoint Higgs fields, respectively. Note that these conditions do not include scale-independent terms since we use the scheme for the renormalization. These equations then yield

(11) |

and the upper relation allows one to evaluate from the coupling constants of the SM gauge interactions at the GUT scale determined using the RGEs [35, 36, 37].

Before concluding this subsection, we discuss the qualitative dependence of the proton decay lifetime on the MSSM parameters. As already mentioned above, the loop functions for the diagrams in Fig. 1 give rise to a factor of . Therefore, the proton lifetime becomes longer if the sfermion masses are taken to be larger. In addition, as can be seen from (21), the decay amplitude contains both the up- and down-type Yukawa couplings, which leads to a factor of . Moreover, the Higgsino-exchange contribution also has an extra factor of . As a result, the proton decay rate is strongly enhanced for moderate/large values of . For these reasons, large sfermion masses and small are favorable for evading the proton decay constraints.

### 3.2 Elastic Scattering Cross Sections

Next, we review the calculation of the neutralino-nucleus elastic scattering cross sections that we use in the following analysis. There are two types of interactions that induce dark matter-nuclei scattering: spin-independent (SI) and spin-dependent (SD). Since there is no interference between these two interactions, we can evaluate the SD and SI scattering cross sections separately.

We first consider SI scattering. The SI elastic scattering cross section of the neutralino LSP with a nucleus is expressed in terms of the SI neutralino-nucleon effective coupling () as follows:

(12) |

where and are the masses of the neutralino LSP and the target nucleus, respectively, and () denotes the atomic (mass) number of the target nucleus.

The SI neutralino-nucleon scattering matrix elements are induced by the exchange of squarks and neutral Higgs bosons. To evaluate the effective coupling , we first obtain the neutralino-quark/gluon effective operators by integrating out the squarks and Higgs bosons. Then, using the nucleon matrix elements of these effective operators, we can calculate the effective coupling . For more details, see [82, 83, 84, 85, 86, 87]. As a result, is expressed in terms of the coefficients of the neutralino-quark effective scalar interactions, (where and denote the neutralino LSP and quarks, respectively), as [84, 85]

(13) |

Here the are the quark masses,
is the nucleon mass, the are the nucleon matrix elements of the
light-quark mass operators, and denotes the gluon contribution to the nucleon mass.
We extract the values of from the pion-nucleon
-term MeV and MeV
[88] ^{4}^{4}4This choice corresponds to , which is larger than the values obtained in
lattice QCD simulations [89]. We note that the
uncertainties in these quantities significantly affect the resultant
scattering cross sections [85].. Analytic
expressions for the are presented in
[84, 85]. The second term of the right-hand side in
Eq. (13) represents the (long-distance) contribution of heavy
quarks to the neutralino-gluon interactions, which can be related to the
quark couplings via the triangle diagrams associated with
the trace anomaly of the energy-momentum tensor
[90]. In Eq. (13), we neglect the effects of the
twist-2 operators [82] as well as the short-distance
contribution of quarks to the gluon operators [86, 87], since they are numerically small when squarks are
rather heavy, which is the case we discuss below.

In the models we study, the dominant contribution to is given by the exchange of neutral Higgs bosons, since the squarks tend to be heavy. Moreover, in a wide range of parameter space, the Higgs sector is close to the decoupling limit, and the LSP is a bino-Higgsino mixed state with , , . In this case, the expression for is approximated by

(14) |

for the bino LSP case, while

(15) |

for the Higgsino LSP case. Here, is the weak mixing angle, () is the mass of the lighter (heavier) neutral Higgs boson, and and for up- and down-type quarks, respectively. As can be seen from these expressions, the neutralino-nucleon scattering cross sections decrease when the difference between and gets large. In addition, it is found that the SI effective coupling depends on the sign of and, in particular, when is negative the coupling can be significantly suppressed due to cancellations (this feature is sometimes called the “blind spot” [84, 91, 92]).

We next discuss SD scattering, for which the neutralino-nucleus scattering cross section is given by

(16) |

where is the Fermi constant, is the total spin of the target nucleus, and

(17) |

with () being the expectation value of the total spin of protons (neutrons) in the target nucleus. The SD neutralino-nucleon effective coupling, , is expressed as

(18) |

where we use given in [85], and denotes the SD neutralino-quark couplings, which are induced by the exchange of -boson and squarks. The analytic formula for is again given in [84, 85]. As in the SI case, the squark contribution is suppressed compared with the boson contribution in the parameter region we are interested in. Furthermore, when , , , is approximated by

(19) |

for a bino-like LSP, whilst for a Higgsino-like LSP we have

(20) |

where denotes the third component of the SU(2) generators.

As seen above, the neutralino-nucleus scattering cross sections are suppressed when gauginos/Higgsinos are heavy. In such cases, electroweak-loop contributions may dominate the tree-level Higgs and boson contributions [86, 93]. It turns out, however, that in the case of a bino-Higgsino LSP, the electroweak loop contributions are quite small [94], and thus we neglect them in our calculation.

Because of the coherent nature of the SI neutralino-nucleus scattering as shown in Eq. (12), the current and future direct detection experiments are much more sensitive to the SI scattering compared to SD scattering. For this reason, we mainly discuss SI scattering in the following.

## 4 Results

### 4.1 Cmssm

In view of the discussion in Section 3.1, in our study of the proton lifetime we focus on relatively small values of , and have chosen in Fig. 2. For larger values of , the proton lifetime becomes smaller than the current experimental bound, and minimal SU(5) is not viable.

We show in Fig. 2 four examples of planes in the CMSSM with . In the left panels we choose , whereas in the right panels we choose . We take in the upper panels and in the lower panels. Higgs mass contours are shown as red dot-dashed curves labelled by in GeV in 1 GeV intervals starting at 122 GeV. We recall that to calculate we use FeynHiggs [54], which carries a roughly 1.5 GeV uncertainty. In the left panels, the light mauve shaded region in the upper part of the figure is excluded because there are no solutions to the Higgs minimization conditions: along this boundary . Electroweak symmetry breaking (EWSB) fails here because the Higgs soft masses at the GUT scale are large and the RG running to the weak scale does not suppress the Higgs soft masses sufficiently for EWSB to occur. Because large gluino masses can assist electroweak symmetry breaking effectively at two loops, the value of that is allowed increases for increasing . Just below the region where EWSB fails, there is a dark blue shaded region where the relic density falls within the range determined by CMB experiments [95]. Since the relic density of dark matter is now determined quite accurately (), for the purpose of visibility we have shown the strip for which the relic density lies in the range . This strip is in the focus-point region [23, 47] where the Higgsinos are much lighter than the stops. The correct dark matter density is realized either by coannihilation of the Higgsino with the bino when is smaller, or by Higgsino annihilations when the Higgsino mass is of order a TeV for larger values of . The TeV-scale Higgsino dark matter region continues well beyond the bounds of the figure. We note also that the brown shaded regions at the bottoms of the panels are excluded because there the LSP is the lighter charged stau lepton. The planes also feature stau-coannihilation strips close to the boundary of the brown shaded region. They extend to about TeV, but are very difficult to see on the scale of this plot, even with our enhancement of the relic density range. We note that for this value of there are no relevant constraints from rare B decays.

Contours of the proton lifetime using down-type Yukawa couplings (see the discussion given in Sec. 3.1) are shown as solid black curves that are labelled in units of yrs. Thus the limit of yrs would exclude everything below the curve labelled 0.05. For the nominal value of GeV, neglecting the theoretical uncertainties in the calculation of , we see that in the upper left plane of Fig. 2 the Higgs contour intersects the focus-point region where yrs. Much of the focus-point strip in this figure may be probed by future proton decay experiments. Changing the sign of has almost no effect on the proton lifetime, as seen in the lower left panel of Fig. 2, but the calculated Higgs mass is smaller by GeV, which is less than the uncertainty in the FeynHiggs calculation of .

In the right panels of Fig. 2, since increasing drives a larger splitting between the two stops, there are excluded regions shaded in brown in the upper halves of the panels, where the lighter stop becomes the LSP. Close to this boundary (but again difficult to see) there is a narrow blue strip where the dark matter density is brought into the allowed range by coannihilation with the lighter stop. The relatively large value of leads to large Higgs masses in most of the plane, but the Higgs mass is somewhat smaller along this strip for . We recall that the Higgs mass is sensitive to the off-diagonal element in the stop mass matrix, which is proportional to and peaks when is roughly 2.5 times the geometric mean of the two stop masses. For positive , is relatively large and we are past the peak where is maximized. In contrast, for , is smaller (as there is some cancellation between the two terms) and we are closer to the peak and is larger. This effect is pronounced along the upper left edge because the stop is much lighter in this region. For , the GeV contour, which is consistent with the experimental value when uncertainties in the theoretical calculation are considered, intersects the strip at TeV where yrs. Note that we terminate the larger contours where the calculated value becomes unreliable: near the endpoints of these curves, the uncertainty in the FeynHiggs calculation of is GeV. For , the GeV contour intersects the stop coannihilation strip when TeV and the proton lifetime is significantly smaller ( in these units). For slightly lower than the value shown in these panels of Fig. 2, large uncertainties in from FeynHiggs appear when yrs. When the stop coannihilation strip is no longer present. On the other hand, when the central value of the Higgs mass along the stop strip drops below 122 GeV when , which is unacceptably small.

We show in Fig. 3 the spin-independent cross section, , as a function of the neutralino mass for the two upper panels in Fig. 2 with . The points in each panel represent results of a scan of the parameter space. In the upper panels, darker points fall within 3 of the dark matter relic density that fits best the Planck data. Lighter points have smaller relic densities and should not be excluded. However, whenever the relic density is below the central value determined by Planck, we scale the cross section downwards by the ratio of the calculated density to the Planck density. From these panels, we find that the cases give relatively large SI scattering cross sections, while those for are significantly suppressed. In the case of , the values of and are close to each other, and thus the LSP is a well-mixed bino-Higgsino state. This leads to a large SI scattering cross section, as can be seen from Eqs. (14) and (15). The set of darkly shaded points with good relic density are found mostly at GeV due to the fact that these points are mainly Higgsino LSPs. As the bino mass is increased, the scattering cross section decreases. However, the points sampled here all have GeV which produces the lower boundary of the points displayed. Because of the constraints coming from the Higgs mass, the scattering cross sections for all dark matter candidates in these models are accessible at LZ. Along the focus point, the LSP mass varies downward as is decreased and the cross section is maximal at around 3 pb. Due to the small uncertainty in the Planck relic density, we find very few darkly shaded points in this region. Points with smaller cross section are found between the focus point strip and the no-EWSB boundary where the relic density is below the Planck density. On the other hand, for the cases, the LSP is almost pure bino and squarks are quite heavy; for these reasons, we obtain very small SI cross sections. The solid curve in Fig. 3 corresponds to the current LUX limit [51] and thus some models (including the focus point models) are excluded by this limit when . The thin black dashed curve corresponds to the projected LZ sensitivity [53, 96] and almost all of the points sampled when are therefore testable. The thick orange dashed line corresponds to the irreducible neutrino background [97, 96]. All of the points sampled when fall below the neutrino background and probing them would require a directional recoil detector [98].

In the lower panels of Fig. 3, we see the same points, now colored to show the Higgs mass ranges. The darkest points have the calculated Higgs mass in the range 124–126 GeV, medium shaded points have in the ranges 123–124 GeV or 126–127 GeV, and the lightest points have in the ranges 122–123 GeV or 127–128 GeV. All of these are compatible with the experimental measurement, within twice the FeynHiggs uncertainty.

In Fig. 4, we show the points corresponding to the lower panels of Fig. 2 with . There is relatively little change in the scattering cross sections for . When , the cross sections are in general somewhat lower due to the cancellation mentioned in Sec. 3.2. For the points have moved to lower , but remain for the most part unobservable. We again see that the Higgs mass constraint puts the bulk of the points within reach of LZ.

Finally, in Fig. 5 we show the spin-dependent cross sections, , for the upper panels in Fig. 2 when . The points and shadings are identical to those in the previous two figures. Here the thick black solid curve is the upper limit from PICO [99] and the thin curves are obtained from IceCube [100] limits based on annihilations into pairs (solid) or pairs (dashed). For the focus-point models, annihilations proceed primarily into electroweak gauge bosons, or final states with some non-negligible contributions from , for which the may be applicable. Models with lie just below the current bounds again because of the highly-mixed nature of the LSP, whilst the models with predict cross section far below these bounds.

### 4.2 mSUGRA

In [21], mSUGRA models were considered with (the Polonyi [63] value) and for comparison. The computed value of is generally . As a result, proton lifetime limits for these cases would be short, in violation of the experimental bounds unless a non-minimal version of SU(5) is adopted.

In the case of the Polonyi model, viable regions of the parameter space contain a gravitino LSP, which would give negligible signals in direct detection experiments. However, for the larger value there are some regions of parameter space with a bino LSP with the relic density held in check by stau coannihilations, but for GeV and there are strong constraints from B-physics observables in that case. This model is highly constrained and we do not discuss it further here, though we return later to mSUGRA models with universality imposed below the GUT scale.

### 4.3 subGUT

In the left panels of Fig. 6, we show examples of planes with , , GeV and (upper panel) and (lower panel). In both planes, one finds three distinct brown shaded regions where the LSP is no longer neutral and/or uncolored. At the left, at low and TeV the lighter stop is the LSP, at low for all the lighter stau is the LSP, and along a diagonal strip that rises from the stau LSP region the lighter chargino is the LSP. When is approximately equal to the bino mass, , two of the neutralino mass eigenstates are strongly mixed bino-Higgsino states. In this case, 1-loop corrections to these masses can differ significantly from the correction to the second Higgsino (which is nearly identical to the correction to the lighter chargino) and cause the chargino to become the LSP. To the left of this region, the LSP is the bino and the relic density gets too large. To the right of this region, the LSP is the Higgsino and, because and are large here, the relic density again gets too large.

There are also three distinct regions in these panels of Fig. 6 where the relic density is consistent with the Planck constraint. Somewhat offset from the stop LSP region, we see a curved band which is produced by stop coannihilation. This region is much broader than in the typical CMSSM case, due to the increased degeneracy of the SUSY particles. Then, just above the stau LSP region we see the familiar stau coannihilation strip. At this value of it extends to far greater values of than it would in the CMSSM with : this is generally possible for sufficiently small , and is again due to the degeneracy of the SUSY particles. We note that, to the right of the chargino LSP region, stau coannihilation occurs between the stau and a Higgsino LSP instead of the more usual bino. Finally, between the stau and stop strips, we see a narrow funnel region where . For , the funnel extends to TeV, whilst for , it extends past the end of the plot to TeV. Looking now at the red dashed contours of , we see that in the stop coannihilation strip and the funnel region the Higgs mass is somewhat too small: GeV for , whereas the uncertainty from FeynHiggs is GeV. This is to be expected, since the stop masses are light in this region and so the corrections to the Higgs mass are small. On the other hand, the Higgs mass exceeds 124 GeV for , and hence much of the stau coannihilation region is acceptable. For (lower left panel of Fig. 6), the values of are somewhat higher, and parts of the stop coannihilation strip may be acceptable.

As in the previous subsections, the solid black lines are contours of the proton decay lifetime. For (where the Higgs mass is acceptable), the lifetime exceeds the experimental bound of yrs. As discussed earlier, in order to obtain a sufficiently long lifetime, we are forced to relatively small values of . For larger than the value shown here, the lifetime along the stau strip drops below the experimental bound, as seen in the right panels of Fig. 6 where is chosen, and one would need to abandon minimal SU(5), as the proton lifetime is less than yrs over much of the plane.

Qualitatively, we see similar features for the LSP and relic density in the right panels. Smaller () is possible, but one needs to go to higher values of to ensure a sufficiently heavy Higgs boson. For smaller , the extent of the stau strip is reduced, making it difficult to obtain both a heavy enough Higgs and a long proton lifetime. This reduction in the stau coannihilation strip is due to a reduction in the Higgsino masses as is reduced. The stau coannihilation band in this figure is actually assisted by several other supersymmetric particles with masses similar to the stau, including the charged and neutral Higgsinos. The Higgsino masses are roughly set by ( is set by EWSB conditions which is related to the stop mass and hence to ). Since the bino and stau masses continue to grow as is increased, eventually the Higgsino becomes the LSP with no potential coannihilation partners. At this point, the stau coannihilation band disappears and the Higgsino becomes the LSP. Lowering reduces the value of for which the Higgsino becomes the LSP, and so reduces the size of the stau coannihilation strip. However, if we rely on non-minimal SU(5) to lengthen the lifetime of the proton, we can go to larger values of , as seen in the right panels of Fig. 6. This allows one to obtain simultaneously a large enough Higgs mass and a small enough relic density.

Results for the SI cross section for the subGUT CMSSM cases displayed in Fig. 6 are shown in Figs. 7 (for ) and 8 (for ). We see in the upper panels of Fig. 7 that the SI cross sections for models with a relic density compatible with the Planck range (darker blue points) are generally below the current upper limit but within reach of the LZ experiment. The dark shaded points originate from what appears to be the stau coannihilation strip. As noted above, obtaining the correct relic density at such large LSP masses requires additional coannihilation and mass degeneracies. In this case, the Higgsinos also have masses comparable to the bino mass. Along this horizontal strip of points in the left panel of Fig. 7, the mass difference () is relatively constant and hence from Eq. (14), when is small, we obtain a cross section which is relatively constant as well. The same is true for most of the under-dense (paler blue) points for (upper right panel), but under-dense points for may have SI cross sections below the LZ sensitivity though above the neutrino background level. The lower panels of Fig. 7 show that many of the points within reach of the LZ experiment have values of close to the experimental value. Fig. 8 (for ) exhibits somewhat lower values of the SI cross section, in general. Note that the cluster of points with GeV correspond to bino LSPs in the stop coannihilation region. These points do not appear when as that region has GeV and hence not included in our scan for elastic cross sections. Consequently, most models with a relic density compatible with the Planck range are beyond the LZ sensitivity, and some of the points are below the neutrino background level. The same is true a fortiori for the points with under-dense relic neutralinos. Finally we note in passing, that the 3 nearly horizontal points below the neutrino background (for ) originate from the funnel region (there are similar points when but more difficult to discern in the figure).

Turning now to subGUT mSUGRA models, we consider only the Polonyi model, i.e., ,
with .
The left panel of Fig. 9 is for GeV.
The mauve shaded regions in the upper left
and lower right parts of the plane are where electroweak symmetry breaking fails
( in the upper left and a diverging Yukawa coupling due to an excessive value for in the lower right), and in the
central brown shaded region the stau is the LSP. Above this region, various processes contribute to bringing
the relic density into the Planck range. Over much of this plane, the LSP is mostly Higgsino which is nearly degenerate with the next lightest superparticle (NLSP) which is a chargino in this case as well as the 2nd Higgsino.
In the blue shaded area above the stau LSP region,
in addition to neutralino coannihilations, stau coannihilation also enhances the cross section in this strip.
In the wide blue shaded region above the stau strip (recall that we are here showing regions where the relic
density lies between .06 and 0.2) the Higgs funnel (lower part of this strip and conventional focus-point region (upper part) have merged.
Below the stau LSP region, the gravitino becomes the LSP and can be dark matter. The values of
are quite large for all points in the left panel of Fig. 9. For this reason, the proton lifetime in
minimal SU(5) models is much too short, and some non-minimal model must be considered ^{5}^{5}5We note also
that the constraint excludes a region at small and (shaded green), and that the
constraints is also relevant in much of the allowed region of the plane..
The right side of Fig. 9 is for GeV, and shares the qualitative features
of the electroweak symmetry breaking and stau LSP constraints. The values of are somewhat
smaller than in the GeV case, but still much too large to obtain a sufficiently long proton lifetime in minimal SU(5) models.
The dark matter constraint is satisfied in a focus-point strip close to the electroweak symmetry breaking boundary,
which has now demerged from the funnel and stau coannihilation strip. The rapid-annihilation funnel is now
clearly visible as a separate well-defined region.
In the stau strip and in the funnel, the LSP is once again a bino, though the masses of the Higgsinos are
not much larger. There is also
a gravitino dark matter region below the stau LSP region.

Fig. 10 displays results for the SI cross section in these subGUT mSUGRA models. We see in the upper panels that the SI cross section is generally between the current LUX upper limit and the prospective LZ sensitivity, though some models (particularly for GeV) have cross sections above the LUX limit and a few under-dense models have SI cross sections below the LZ sensitivity. For GeV, we clearly see the pile of points with the Planck relic density at GeV corresponding to a Higgsino LSP near the broad intersection of the stau strip and focus point swath. The dark blue points in this figure continue to higher Higgsino masses along the stau coannihilation strip. Very low mass points (all lightly shaded) correspond to regions in Fig. 9a that are to the left of the blue shaded region. In the white region to the left, the relic density is small, and in the white region to the right (between the stau strip and funnel) the relic density is too high. For GeV, we see two very distinct regions in Fig. 10. The region with lower masses ( GeV and cross section between and pb) originate from the focus point region. The remainder of the points come from either the funnel or the stau strip and can be more easily distinguished by the lower panels showing the Higgs mass ranges. The relative paucity of dark shaded blue points stems from the fact that the true Planck strips are quite thin in this case. Note also that there is no pile up of points at 1100 GeV as the LSP is most bino rather than Higgsino at the higher value of . We see in the lower panels of this figure that models with GeV lie in the intersection region for GeV, and as noted above the dark brown shaded points for GeV at low masses come from the focus point where as we now see that the middle group around GeV originate in the funnel, and the group at larger masses lie in the stau strip. All of the dark shaded points lie within reach of the LZ experiment (though some GeV models are excluded already by the LUX upper limit).

We conclude this subsection by showing results for the SD cross sections in these subGUT mSUGRA models in Fig. 11. We see in the upper panels that the SD cross sections are generally smaller than the PICO bound [99], and also below the IceCube upper limits [100] for both and final states (which are likely to be more similar to the model final states). There is a handful of GeV models (most of them under-dense) whose predictions lie close to the IceCube limit, but most model predictions are significantly below it. We see in the lower right panel of Fig. 11 that many of the models close to the IceCube limit have FeynHiggs values close to the experimental value. Comparing this figure with Fig. 10, it seems that there are better prospects for discovering SI scattering in these subGUT mSUGRA scenarios. However, we recall that these models yield proton lifetimes that are too short in minimal SU(5), pointing to the need for some non-minimal model.

### 4.4 Nuhm1

As mentioned earlier, in the NUHM1 one has the freedom to treat as a free parameter, and in the following we show planes for some representative choices of , and . In fact, as long as is small enough to obtain an acceptable proton lifetime, the qualitative behaviour of the parameter space is relatively insensitive to , though there is some dependence of the Higgs mass contours on , as could be expected.

We show in the upper left panel of Fig. 12 the plane for , and GeV, which exhibits a small stau LSP region at low and . Since is fixed, the composition of the LSP changes as is increased. At small the LSP is mainly bino and the relic density is too high. As is increased, the Higgino component increases and the relic density passes through the Planck range across a relatively narrow, near-vertical transition strip. (Note that, in all four panels of this figure, the blue region corresponds to just the 3 Planck range rather than the extended range used in previous figures.) At larger the LSP is a Higgsino with a mass of about 1050 GeV which is slightly low for a Higgsino LSP and, as a result, the relic density is somewhat too small when TeV. In this panel, we see that we obtain an acceptable Higgs mass ( GeV) when TeV. The proton lifetime is sufficiently large ( yrs) for this value of .

In the upper right panel of Fig. 12, we have increased slightly to 1050 GeV. The most striking feature is that the dark matter region fills the right part of the plane: indeed, it extends infinitely far to the right towards large gaugino masses. In this case, when the gaugino mass is large, the LSP is a nearly pure Higgsino, as is the NLSP. This near-degeneracy facilitates coannihilation that brings the relic density within the acceptable range, with being determined predominantly by [50]. The Higgsino mass in this case is very close to 1100 GeV, which remains constant at large . Thus there is a very large (infinite) area where the relic density matches the Planck result. At low , the relic density is too large and drops monotonically as the gaugino mass is increased and asymptotes to the Planck density at very large . For TeV, when , the Higgs mass contours are nearly vertical and the value of was chosen to maximize the area with good relic density and Higgs masses. The area between TeV and 9 TeV has between 124 and 126 GeV, and increasing (decreasing) by 0.5 would raise (lower) by roughly 1 GeV. Much of this region has yrs: requiring yrs implies either TeV for small or TeV for TeV. As is increased past 1050 GeV, the left edge of the blue shaded region moves quickly to the right and the relic density would be too large over much of the plane. The relic density would now asymptote to a value in excess of the Planck density. The Higgs mass is independent of for small but the Higgs mass contours bend to the left as is increased, so that the Higgs mass becomes large at larger .

In the lower panels of Fig. 12, we have taken . In the left panel, and the transition strip from bino to Higgsino dark matter is much narrower and occurs at much lower GeV. Had we chosen as in the previous plots, the Higgs mass would be far too small. This can be compensated in this panel by choosing larger , and we have chosen in this panel. As in the CMSSM, there is now a shaded region where the LSP is a stop in the upper left of the panel. There is a barely visible stop coannihilation strip that runs close to the stop LSP boundary, from the transition strip down to smaller and . There is also a narrow stau coannihilation strip running on top of the stau LSP region at low . In the right panel, we have again taken and increased to 1000 GeV. The relic density region resembles that in the upper left panel of the same figure, though the Higgs masses are now notably larger. The transition strip is now centered on GeV, which is compatible within the experimental measurement within the theoretical uncertainties. We note also that the proton lifetime is far smaller in the lower panels due to the larger value of . Indeed, in the lower left panel is always below 0.01 yrs.

The elastic scattering cross sections for the four panels of Fig. 12 are shown in Figs. 13 and 14. The left panels of Fig. 13 correspond to the NUHM1 model with , , and GeV. Viable points (with the correct relic density or less) have gaugino masses of around 3 TeV (for the correct relic density) or greater (less than the Planck density). In either case, the LSP mass is just over 1 TeV, which explains why all the points line up vertically. Most of the points (though not all) lie below the current LUX limit and all of them lie above the LZ projected reach. Note that, in principle, this vertical strip could extend further down, into the neutrino background, if we continued to sample points at higher . Our sampling of points includes only points with between 122 and 128 GeV. Concerning the right panels with GeV, we recall that much of the plane contains a Higgsino LSP with the desired relic density. In that region, the mass of the LSP is always very close to 1100 GeV and that fact is readily seen in the right panels of Fig. 13, where all the points stack vertically at GeV. All of these points lie below the current LUX bound, but most of them are within the projected reach of LZ. As in the previous example, the vertical strip of points could go lower if we sampled to larger where GeV. As seen in the lower right panel of this figure, the points with between 124 and 126 GeV are all accessible to LZ.

In Fig. 14, we show the spin-independent cross sections for . In the left panels with GeV, we again see a pile-up at a fixed LSP mass, GeV, corresponding in this case to the transition strip and the region to its right. Since the strip is so narrow in this case, there are very few dark-shaded points, and these have cross sections that exceed the LUX bound. There are also a few points at lower that originate near the stop or stau coannihilation strips. All of the points shown lie within the LZ projected reach. One also sees in the lower left panel that many of the points have a Higgs mass in the 124–126 GeV range. Indeed, for an underdense Higgsino-like LSP with GeV (to the right of the transition strip in Fig. 12), there are also many points with GeV, though these points are eclipsed in Fig. 14 by those with more favorable . In the right panels of Fig. 14, for GeV, all points in the transition strip and to the right of the strip have a narrow range of LSP masses fixed by the value of . These points lie just below the current LUX bound. Points that are most compatible with Planck dark matter abundance and have GeV have the largest SI cross sections, within an order of magnitude of the current LUX limit. All points considered here, with 122 GeV GeV, should be accessible with LZ.

## 5 Discussion

Large parts of the CMSSM parameter space are excluded by the absence (so far) of proton decay, if the CMSSM is embedded within the minimal SU(5) GUT. There are regions of parameter space with and several TeV that are still allowed, however. Typically, these models predict a spin-independent dark matter-nucleon scattering cross section that falls below the current LUX upper limit but could be accessible to the planned LZ experiment. In fact, it is the constraints coming from the Higgs boson mass which exclude the bulk of the model points which are beyond the reach of LZ. The prospects for direct detection of spin-independent dark matter scattering are reduced for and for , and the cross sections for spin-dependent dark matter scattering are generally substantially below the current bounds from PICO and IceCube.

In view of the limited perspectives within the CMSSM, we have explored in this paper the prospects for probing other MSSM scenarios via proton decay and dark matter detection. In one class of scenarios, called subGUT models, universality of the soft supersymmetry-breaking masses is retained, but is imposed at some scale . Within the subGUT CMSSM there are more possibilities for bringing the supersymmetric relic density within the range allowed by Planck and other experiments even if and are each several TeV, thanks in particular to the more compressed spectrum and consequently the greater possibilities for coannihilation processes that bring the dark matter density down into the allowed range. However, small values of are preferred again, as in the GUT-scale CMSSM, unless one adopts a non-minimal GUT structure. For GeV and or 10, we find spin-independent dark matter-nucleon scattering cross sections that are well within the range allowed by LUX, and the spin-independent cross section may fall below the neutrino background level, particularly for .

In mSUGRA models, the possibilities are very limited if , but open up for . On the other hand, in mSUGRA models is no longer a free parameter, and the electroweak vacuum conditions typically require large values that give severe problems within the minimal SU(5) GUT framework. That said, spin-independent dark matter scattering may again lie within reach of the LZ experiment, whereas spin-dependent scattering cross sections generally lie below the PICO and IceCube upper limits.

In the NUHM1 one may regard the Higgs mixing parameter as an extra free parameter compared to the CMSSM. This freedom opens up new possibilities for models that respect the dark matter, Higgs mass, and proton decay constraints. In particular, since varying varies the Higgsino component of the LSP, there is the possibility of a ‘well-tempered’ transition region as well as the more familiar stau and stop coannihilation possibilities for bringing the relic neutralino density into (or below) the Planck range. Moreover, for some values as seen in the upper right panel of Fig. 12, in particular, the relic density may lie within the Planck range up to indefinitely high values of and . In this case, the proton lifetime may certainly be long enough to survive the present experimental lower limit whereas, as seen in Fig. 13, the spin-independent dark matter scattering cross section is likely to be within reach of the planned LZ experiment.

In conclusion, the examples studied in this paper show that there are certainly interesting possibilities for probing supersymmetric models beyond the CMSSM via searches for proton decay and direct dark matter scattering.

## Acknowledgements

The work of J.E. was supported in part by the London Centre for Terauniverse Studies (LCTS), using funding from the European Research Council via the Advanced Investigator Grant 267352 and from the UK STFC via the research grant ST/J002798/1. The work of F.L. was also supported by the European Research Council Advanced Investigator Grant 267352. The work of J.L.E., N.N. and K.A.O. was supported in part by DOE grant DE-SC0011842 at the University of Minnesota. The work of N.N. was also supported by Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists. The work of P.S. was supported in part by NSF Grant No. PHY-1417367. P.S. would also like to thank CETUP* (Center for Theoretical Underground Physics and Related Areas) for its hospitality and partial support during the 2015 Summer Program.

## Appendix

## Appendix A Wilson Coefficients and RGEs for Proton Decay

In this Appendix, we summarize the matching conditions and RGEs used in the proton decay calculation discussed in Sec. 3.1. The Wilson coefficients and in Eq. (6) are given at the GUT scale by

(21) |

where and are up- and down-type Yukawa couplings at the GUT scale, respectively, is the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and denote the extra phases appearing in the GUT Yukawa couplings. They are taken such that they satisfy , and thus there are two independent degrees of freedom [72]. For the defi