Likelihood analysis of supersymmetric SU(5) GUTs

We perform a likelihood analysis of the constraints from accelerator experiments and astrophysical observations on supersymmetric (SUSY) models with SU(5) boundary conditions on soft SUSY-breaking parameters at the GUT scale. The parameter space of the models studied has seven parameters: a universal gaugino mass \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{1/2}$$\end{document}m1/2, distinct masses for the scalar partners of matter fermions in five- and ten-dimensional representations of SU(5), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_5$$\end{document}m5 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{10}$$\end{document}m10, and for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {5}$$\end{document}5 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{\mathbf{5}}}$$\end{document}5¯ Higgs representations \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{H_u}$$\end{document}mHu and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{H_d}$$\end{document}mHd, a universal trilinear soft SUSY-breaking parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_0$$\end{document}A0, and the ratio of Higgs vevs \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tan \beta $$\end{document}tanβ. In addition to previous constraints from direct sparticle searches, low-energy and flavour observables, we incorporate constraints based on preliminary results from 13 TeV LHC searches for jets + events and long-lived particles, as well as the latest PandaX-II and LUX searches for direct Dark Matter detection. In addition to previously identified mechanisms for bringing the supersymmetric relic density into the range allowed by cosmology, we identify a novel \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\tilde{u}_R}/{\tilde{c}_R} - \tilde{\chi }^{0}_{1}$$\end{document}u~R/c~R-χ~10 coannihilation mechanism that appears in the supersymmetric SU(5) GUT model and discuss the role of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\tilde{\nu }}_\tau }$$\end{document}ν~τ coannihilation. We find complementarity between the prospects for direct Dark Matter detection and SUSY searches at the LHC.

It is interesting to explore also models that are less simplified than the CMSSM, but not as agnostic as the pMSSM, in that they incorporate a limited number of simplifying assumptions. GUTs motivate the assumption that the gaugino masses are universal, and constraints on flavour-changing neutral interactions suggest that the soft SUSY-breaking masses for scalars with identical quantum numbers are also universal. However, there is no compelling phenomenological reason why the soft SUSY-breaking masses for scalars with different quantum numbers should be universal.
Specific GUTs may also provide some guidance in this respect. For example, in an SO (10) GUT the scalar masses of all particles in a given generation belonging to a single 16 representation of SO(10) would be universal, as would those for the 5 and5 SU(5) Higgs representations that belong to a single 10 of SO (10) and break electroweak symmetry, as in the NUHM1. In contrast, the SU(5) framework is less restrictive, allowing different masses for scalars in5 and 10 representations [72], and also for the 5 and5 Higgs representations. Thus it is a 1-parameter extension of the NUHM2. In this paper we explore the theoretical, phenomenological, experimental and cosmological constraints on this SU(5)based SUSY GUT model. This relaxation of universality is relevant for the evaluation of several different constraints from both the LHC and elsewhere. For example, the most powerful LHC constraints on the CMSSM, NUHM1 and NUHM2 are those from the classic / E T searches [1][2][3][6][7][8]. These constrain principally the right-handed squarks, whose decays are dominated by theq R → qχ 0 1 channel that maximizes the / E T signature. On the other hand, the decay chains of left-handed squarks are more complicated, typically involving theχ ± 1 , resulting in a dilution of the / E T signature and more importance for final states including leptons. In a SUSY SU(5) GUT, the lefthanded squarks and the right-handed up-type squarks appear in 10 representations whereas the right-handed down-type squarks appear in5 representations, with independent soft SUSY-breaking masses. Hence the impacts of the LHC / E T and other constraints need to be re-evaluated.
The possible difference between the soft SUSY-breaking contributions to the masses of the squarks appearing in a 10 of SU(5), i.e., up-type squarks and left-handed downtype squarks, and those appearing in a5 of SU(5), i.e., right-handed down-type squarks, offers a new avenue for compressing the stop spectrum. Also, as we shall see, with m 5 = m 10 there is the possibility that mũ R ,c R are much smaller than the other squark masses, leading to another type of compressed spectrum. 1 In principle, the constraints from flavour observables may also act differently when m 5 = m 10 . For example, the soft SUSY-breaking masses of the left-and right-handed charge +2/3 quarks are independent, and flavour observables such as BR(b → sγ ) and BR(B s → μ + μ − ) depend on both of them, in general.
Another experimental constraint whose interpretation may be affected by the non-universality of scalar masses is (g − 2) μ . A priori, a SUSY explanation of the discrepancy between the Standard Model (SM) prediction and the experimental measurement of (g − 2) μ requires relatively light smuons, either right-and/or left-handed, which are in 10 and 5 representations, respectively. It is interesting to investigate to what extent the tension between a SUSY interpretation of (g − 2) μ and the LHC constraints on squarks that is present in more constrained SUSY models could be alleviated by the extra degree of freedom afforded by the5-10 disconnect in SU (5).
Finally, we recall that in large parts of the regions of the CMSSM, NUHM1 and NUHM2 parameter spaces favoured at the 68% CL the relicχ 0 1 density is brought into the range allowed by Planck [74] and other data via coannihilation with the stau and other sleptons [75][76][77][78][79][80][81][82]. In an SU(5) GUT, the left-and right-handed sleptons are in different representations,5 and 10, respectively. Hence they have different masses, in general, providing more flexibility in the realization of coannihilation. Specifically, as mentioned above, the freedom to have m 5 = m 10 allows the possibility that the right-handed up-and charm-flavour squarks,ũ R andc R , are much lighter than the other squarks, opening up the novel possibility ofũ R /c R −χ 0 1 coannihilation, as we discuss below. Our analysis of the available experimental constraints largely follows those in our previous studies of other variants of the MSSM [10][11][12][13][14][15][16][17], the main new feature being that we incorporate the constraints based on the preliminary results from LHC searches for jets + / E T events with ∼13/fb of data at 13 TeV [9]. For this purpose, we recast available results for simplified models with the mass hierarchies mg mq and vice versa. We also include the preliminary constraints from LHC searches in 13-TeV data for the heavy MSSM Higgs bosons and long-lived charged particles, and incorporate in combination the recent PandaX [83] and LUX [84] data.
The SUSY SU(5) GUT model we study is set up in Sect. 2, and our implementations of constraints and analysis procedure are summarized in Sect. 3. Section 4 describes how we characterize different Dark Matter (DM) mechanisms, including the novelũ R /c R −χ 0 1 coannihilation mechanism, ν τ coannihilation and a hybrid possibility. Section 5 contains our results in several model parameter planes, and Sect. 6 describes various one-dimensional likelihood functions including those for several sparticle masses, (g − 2) μ and various other observables. Higgs boson branching ratios (BRs) are presented in Sect. 7, followed by a comparison of the SU(5) with the NUHM2 results in Sect. 8. The possibility of a long-livedτ 1 is discussed in Sect. 9, and the prospects for direct DM detection are discussed in Sect. 10. Finally, Sect. 11 presents a summary and some conclusions.

Supersymmetric SU(5) GUT model
We assume a universal, SU(5)-invariant gaugino mass parameter m 1/2 , which is input at the GUT scale, as are the other SUSY-breaking parameters listed below.
We assume the conventional multiplet assignments of matter fields in the minimal superymmetric GUT: where the subscript i = 1, 2, 3 is a generation index. The only relevant Yukawa couplings are those of the third generation, particularly that of the t quark (and possibly the b quark and the τ lepton) that may play an important role in generating electroweak symmetry breaking. In our discussion of flavour constraints, we assume the MFV scenario in which generation mixing is described by the Cabibbo-Kobayashi-Maskawa (CKM) model. This is motivated by phenomenological constraints on low-energy flavour-changing neutral interactions, as is our assumption that the soft SUSYbreaking scalar masses for the different 10 i and5 i representations are universal in generation space, and are denoted by m 10 and m 5 , respectively. In contrast to the CMSSM, NUHM1 and NUHM2, we allow m 5 = m 10 . We assume a universal soft trilinear SUSY-breaking parameter A 0 . We assume the existence of two Higgs doublets H u and H d in 5 and5 representations that break electroweak symmetry and give masses to the charge +2/3 and charge −1/3 and −1 matter fields, respectively. It is well known that this assumption gives a (reasonably) successful relation between the masses of the b quark and the τ lepton [85][86][87], but not for the lighter charge −1/3 quarks and charged leptons. We assume that whatever physics resolves this issue is irrelevant for our analysis, as would be the case, for instance, if corrections to the naive SU(5) mass relations were generated by higher-dimensional superpotential terms [88]. In the absence of any phenomenological constraints, we allow the soft SUSY-breaking contributions to the H u and H d masses, m H u and m H d , to be different from each other, as in the NUHM2, as well as from m 5 and m 10 . As in the CMSSM, NUHM1 and NUHM2, we allow the ratio of Higgs vacuum expectation values, tan β, to be a free parameter.
In addition to these electroweak Higgs representations, we require one or more Higgs representations to break the SU(5) GUT symmetry. The minimal possibility is a single 24 representation , but we do not commit ourselves to this minimal scenario. It is well known that this scenario has problems with rapid proton decay 2 and GUT threshold effects on gauge coupling unification. We assume that these issues are resolved by the appearance of additional fields at or around the GUT scale that are otherwise irrelevant for TeV-scale phenomenology. The effective low-energy Higgsino mixing coupling μ is a combination of an input bilinear H u H d coupling and possible trilinear and higher-order couplings to GUT-scale Higgs multiplets such as H u H d . We assume that these combine to yield μ = O(1) TeV and positive, without entering into the possibility of some dynamical mechanism, and commenting below only briefly on the case μ < 0.

Implementations of constraints and analysis procedure
Our treatments in this paper of many of the relevant constraints follow very closely the implementations in our previous analyses of other supersymmetric models [10][11][12][13][14][15][16].
For the convenience of the reader, we summarise the constraints in Table 1. In the following subsections we review our implementations, highlighting new constraints and instances where we implement constraints differently from our previous work.

Electroweak and flavour constraints
We treat as Gaussian constraints all electroweak precision observables, all B-physics and K -physics observables except for BR(B s,d → μ + μ − ). The χ 2 contribution from BR(B s,d → μ + μ − ), combined here in the quantity R μμ [13], is calculated using a combination of the CMS [89] and LHCb [90] results described in [91] with the more recent result from ATLAS [92]. We extract the corresponding χ 2 contribution in Table 1 by applying to the two-dimensional likelihood provided by the combination of these experiments the minimal flavour violation (MFV) assumption that applies in the SU(5) model. We calculate the elements of the CKM matrix using only experimental observables that are not included in our set of flavour constraints. We have updated our implementations of all the flavour constraints, and now use the current world average value of m t [93]. These and all other constraints whose implementations have been changed are indicated by arrows and boldface in Table 1.

Higgs constraints
We use the combination of ATLAS and CMS measurements of the mass of the Higgs boson: M h = 125.09 ± The χ 2 contributions of 85 Higgs search channels from the LHC and the Tevatron are evaluated using HiggsSignals, see [136,137], where a complete list of references can be found. The χ 2 contributions from the limits from searches for the heavy neutral MSSM Higgs bosons in the H/A → τ + τ − channels are evaluated using the code HiggsBounds [133,[138][139][140], which incorporates the results of CMS searches [132,133] with ∼25 fb −1 of 8 TeV data. The contributions from the two possible production modes, gg → H/A and bb → H/A, are combined in a consistent manner, depending on the MSSM parameters. The results from HiggsBounds have been compared with the published CMS analysis, and are in very good agreement [133]. The corresponding χ 2 contribution is labelled "2D likelihood" in Table 1. For the corresponding constraint with 13 fb −1 of 13 TeV data, we implement an approximate treatment of the χ 2 contribution using the preliminary result of ATLAS [134], as we describe in more detail below. Limits from other Higgs boson searches are not relevant for the investigation in this paper and are therefore not included.

LHC / E T constraints at 13 TeV
ATLAS and CMS have recently announced preliminary results from / E T searches with ∼13/fb of data at 13 TeV, using simplified models for gluino and squark pair production [5,9]. These searches assume mg mq and mq mg, respectively, and 100% BRs for the decaysg → ffχ 0 1 ( f = q, b, t) andq → qχ 0 1 , respectively, which maximize the possible corresponding / E T signatures. Neither of these assumptions is valid in the SUSY SU(5) GUT model: as we will see in more detail later, the mg and mq masses are quite similar in much of the favoured region of parameter space, 4 and in general other decay modes dilute the / E T signature, although larger-multiplicity final states may compensate through an increase in transverse energy H T [141]. These other decay modes populate other search channels including leptons, which we do not consider in this paper as they were of limited importance in our previous analyses of the CMSSM, NUHM1 and NUHM2, having impact only for relatively large squark masses and small m 1/2 . Figure 1 displays the ratios of thegg cross section (left panel) and theqq +qq cross section (right panel) that we find in ranges of mq and mg that are representative of those favoured in our analysis before implementing the LHC 13-TeV / E T constraint, relative to the cross sections found in the simplified models with mg mq and mq mg, respectively. We have used NLL-fast-3.1 [142] to obtain the cross section at NLO + NLL level. In both plots a large area at higher squark masses is visible, as well as a thin strip at ∼500 GeV. The latter corresponds to lighterũ R and c R discussed below. We see that thegg cross section (left panel) is generally smaller than in the corresponding simplified model by a factor >2 due to the destructive interference between the s-channel gluon exchange diagram and the t-channel squark exchange diagram in qq →gg, thus weakening the LHC constraints as discussed below. On the other hand, theqq +qq cross section (right panel) is generally a factor 10 larger than in the simplified model, except in theũ R /c R −χ 0 1 coannihilation strip at small mũ R , mc R , mχ0 1 ∼ 500 GeV and m 1/2 ∼ 2500 GeV, to which we return later. The enhancement of the squark cross section is due to the fact that in the squark-neutralino simplified model there is no production mode with total baryon number B = 2/3, qq →qq, because gluinos are assumed to be absent. On the other hand, in our model mg ∼ min(mq ), and qq →qq (with t-channelg exchange) becomes the dominant squark production mode in the large mq region, due to the valence quark-parton dominance in the proton in the large x regime. Figure 2 displays the CMS 95% confidence limits in the (mg, mχ0 1 ) plane from a hadronic jets plus / E T search [9] within a simplified model assuming that the decay modẽ g → qqχ 0 1 occurs with 100% BR (solid black lines). These limits are compared with the best-fit points (green stars) and the regions in the fits that are preferred at χ 2 = 2.30 and χ 2 = 5.99 (red and blue contours, respectively). Here and in the following analogous parameter planes, we use the χ 2 = 2.30 and χ 2 = 5.99 contours as proxies for the boundaries of the 68 and 95% CL regions in the fit.
In addition, within the 95% CL region in Fig. 2 we have indicated the dominant (>50%)g decays found in our analysis. We note that many model points do not have any decay mode with BR >50% within the 95% CL region and that, for those that do, the dominant decays are two-bodyg →qq modes that were not considered in [9]. Because of this and the fact that thegg cross section is always smaller than in the gluino simplified model by a factor > 2 (see the left panel of Fig. 1), the LHC 13-TeV / E T constraint from the gluino simplified model has only negligible impact. Our LHC 13-TeV / E T constraint on the gluino mass actually comes indirectly from the squark mass constraint estimated using the squark simplified model discussed below, since the squark and gluino masses are related via renormalization group evolution in the SU(5) model. The left panel in Fig. 2 was obtained before implementing the LHC 13-TeV / E T 95% confidence limit on gluino and squark pair production, while in the right panel this constraint is included. We note that the simplified-model exclusion in this analysis extended to mg 1650 GeV, below the gluino mass at the pre-LHC 13 TeV best-fit point, and barely reaching the 68% CL contour (solid red line). Figure 3 contains an analogous set of planes for CMS / E T searches for squarks, where the CMS limit assuming a simplified model with heavy gluino and 100% BRs for q → qχ 0 1 is displayed (black lines): the solid lines assume that all the squarks of the first two generations are degenerate, the dashed lines assume two degenerate squarks, and the dotted lines assume just one squark. The planes in the ) obtained without (with) the CMS 13-TeV constraint. The dominant (>50%)g decays into first-and secondgeneration quarks and squarksq L ,R and third-generation quarks and squarkst/b 1,2 found in the SUSY SU(5) model are colour-coded as indicated upper panels display mχ0 1 and the masses of the first-and second-generation right-handed up-type squarks (here commonly denotedũ R ), while the planes in the lower panels are for the down-type squarks (here commonly denotedd R ). The main decay modes of theũ R (upper) and thed R (lower) are indicated over much of the preferred parameter space, and we note that the dominant (>50%) decay modes of both right-handed up-and down-type squarks are indeed into the corresponding quark flavour +χ 0 1 for nearly the whole 68% CL regions, as assumed in the squark simplifiedmodel search. This is, however, not the case for the left-handed up-and down-type squarks (not shown), whose dominant decays are intoχ ± 1 and electroweak doublet partner quark flavours. Furthermore, within the displayed 95% CL regions there are also large areas where decays into gluinos, not considered in the simplified model, are dominant.
Because theq R → qχ 0 1 decays are important, and also because theqq +qq cross section in our sample is much larger than that found at large mq forqq in the simplified model with mq mg, as seen in the right panel of Fig. 1, we have implemented a recast of this search in our global analy- , assuming a simplified model with heavy gluinos and 100% BR forq → qχ 0 1 : the solid lines assume that all the squarks of the first two generations are degenerate, the dashed lines assume two degenerate squarks, and the dotted lines assume just one squark. All panels show the best-fit point (green star), 68 and 95% CL contours (red and blue lines, respectively) for mχ0 1 and the masses of the first-and second-generation right-handed up-type squarksũ R (upper panels) and the down-type squarksd R (lower panels). In both cases, the left panels were obtained without the CMS 13-TeV constraint, and the right panels include it. The dominant (>50%)q decays found in the SUSY SU(5) model are colour-coded as indicated sis, 5 and the comparison between the left panels (without this contribution) and the right panels (with this contribution) in Fig. 3 shows the importance of this constraint.
Our implementation of the LHC 13-TeV / E T constraint is based on [9]. In this analysis, the CMS Collaboration provides a map of the 95% CL cross-section upper limit as a function of mq and mχ0 1 assuming pp →qq and 100% BR 5 Theũ R /c R −χ 0 1 coannihilation strip visible in the upper panels of Fig. 3 at mũ R mχ0 1 ∼ 500 GeV is the subject of a later dedicated discussion.
forq → qχ 0 1 . This is indeed the dominant production and decay mode in most parts of the 68% CL regions of the considered model, as can be seen in Figs. 1 and 3. For each point we compare our calculation of (σq¯q +σqq ) BR 2 q→qχ 0 1 with the CMS 95% CL upper limit on the cross section: σ UL (mq ,χ 0 1 ).
We model the χ 2 penalty as so that the CMS 95% CL upper limit corresponds to χ 2 (/ E T ) = 5.99 and χ 2 scales as the square of the number of signal events, N sig , which gives the right scaling. We have checked that our implementation (2) reproduces the ±1 σ band in the two-dimensional exclusion limit provided by CMS [9], with a discrepancy that is much smaller than the width of the ±1 σ band.
The aforementioned CMS analysis [9] also looks at three simplified gluino models assuming 100% BR forg → ffχ 0 1 with f = q, b, t, respectively, and provides corresponding cross-section upper limit maps as a function of mg and mχ0 1 .
We implement these constraints by defining χ 2 by analogy with Eq. (2).
We also consider the pp →qg process, treating it as follows. This process is only relevant when mq ∼ mg. In this regime, if mq > mg (mg > mq ),q (g) tends to decay intog (q), radiating soft jets. If these soft jets are ignored, we are left with thegg (qq) system. In this approximation, the impact of pp →qg can therefore be estimated by adding an extra contribution σqg B Rq →qg (σqg B Rg →qq ) to σgg (σqq + σq¯q ). In general, SUSY searches are designed to look for high p T objects, and one loses a small amount of sensitivity by ignoring soft jets. We therefore believe that our implementation of the pp →qg process is conservative.
Finally, we estimate the total χ 2 penalty from the 6

Constraints on long-lived charged particles
We also include in our analysis LHC constraints from searches for heavy long-lived charged particles (HLCP) that are, in general, relevant to coannihilation regions where the mass difference between the lightest SUSY particle (LSP) and the next-to-lightest SUSY particle (NLSP) may be small and the NLSP may therefore be long-lived. As we discuss below, important roles are played in our analysis byτ 1 , χ ± 1 andũ R /c R coannihilation, but only in theτ 1 case is the NLSP -LSP mass difference small enough to offer the possibility of a long-lived charged particle. We implement in our global analysis the preliminary CMS 13-TeV result [130] using tracking and time-of-flight measurements, based 6 One could be concerned that summing up the χ 2 contributions from different simplified-model limits would overestimate the exclusion limit, since these signal regions are not necessarily independent. This would be indeed the case if the same event sample were confronted with multiple overlapping signal regions. In our case, however, the signal sample is divided into statistically independent sub-samples, corresponding to the simplified-model topologiesg → bbχ 0 1 , → ttχ 0 1 , etc., and these sub-samples are confronted with the corresponding simplified-model limits only once. In such a case the χ 2 (/ E T ) estimate (2) provides a conservative limit when there is no significant excess in the data. on the recipe and the efficiency map as a function of the pseudo-rapidity and velocity of the HLCP given in [131]. We use Pythia 8 [143] and Atom [144][145][146][147][148] to generate and analyse the events, and assume that the efficiencies for detecting slow-movingτ 1 s are similar at 8 and 13 TeV. 7 The efficiency contains a lifetime-dependent factor ∝ exp(−dm/ pτ ), where d is a distance d 10 m that depends on the pseudorapidity, and m, p and τ are the mass, momentum and lifetime of the long-lived particle. This factor drops rapidly for particles with lifetimes 10 ps, corresponding to mτ 1 − mχ0 1 1.6 GeV.

Constraints on heavy neutral Higgs bosons from Run II
Concerning the production of heavy neutral Higgs bosons, in addition to the 8 TeV constraints on H/A → τ + τ − provided by HiggsBounds, we also take into account the preliminary exclusion limits obtained by ATLAS from searches for generic spin-0 bosons φ in the τ τ final state with an integrated luminosity of 13.3 fb −1 at 13 TeV that were presented at the ICHEP 2016 conference and described in [134] (see also the CMS results in [135]). Upper bounds on σ × BR(φ → τ τ ) are reported for each M φ separately for the gluon fusion production channel and for production in association with a bb pair assuming there is no contamination between the modes, assuming a single resonance. We compute the cross sections and the BRs in the MSSM using FeynHiggs, adding the contributions for φ = H and φ = A, using the average of the two masses, which are degenerate within the experimental resolution. This result is compared with the upper limit from the corresponding channel neglecting contamination. This approach leads to a conservative limit since we underestimate the signal yield in each channel by neglecting the contamination (the events from the other production mode). As in Eq. (2), the χ 2 penalties are modelled as where . Finally we take the stronger χ 2 rather than combining them, in order to be on the conservative side 8 :

Other constraints
The most important other constraint update is that on spinindependent DM scattering. We incorporate in our global fit the recent result published by the PandaX-II experiment [83], which we combine with the new result from the LUX Collaboration [84], as discussed in more detail in Sect. 8.

Sampling procedure
As discussed in the previous Section, the SUSY SU(5) GUT model we study has seven parameters: m 1/2 , m 5 , m 10 , m H u , m H d , A 0 and tan β. The ranges of these parameters that we scan in our analysis are listed in Table 2. The quoted negative values actually correspond to negative values of m 2 5 , m 2 10 , m 2 H u and m 2 H d : for convenience, we use the notation sign(m 2 ) × |m 2 | → m. The negative values of m 5 and m 10 that are included in the scans may be compatible with early-Universe cosmology [155], and yield acceptable tachyon-free spectra. In the portions of the scans with negative values of m H u and m H d , although the effect of the top quark Yukawa coupling in the renormalization group equations is important, it may not be the mechanism responsible for generating electroweak symmetry breaking, since m H u and m H d are negative already at the input scale. We sample this parameter space using MultiNest v2.18 [156][157][158], dividing the seven-dimensional parameter space into 108 boxes, as also described in Table 2. This has two advantages: it enables us to run MasterCode on many nodes in parallel, and it enables us to probe more efficiently for local features in the likelihood function. For each box, we choose a prior such that 80% of the sample has a flat distribution within the nominal range, while 20% of the sample is in normally distributed tails outside the box. Our resultant total sample overlaps smoothly between boxes, avoiding any spurious features at the box boundaries. The total number of points in our sample is ∼125 × 10 6 , of which ∼8 × 10 6 have χ 2 < 10.

Dark matter mechanisms
The relic density of the LSP, assumed here to be the lightest neutralino,χ 0 1 , which is stable in supersymmetric SU(5) because of R-parity, may be brought into the narrow range allowed by the Planck satellite and other measurements [74] via a combination of different mechanisms. It was emphasized previously [16] in studies of the CMSSM, NUHM1 and NUHM2 that simple annihilations of pairs of LSPs into conventional particles would not have been sufficient to bring the relicχ 0 1 density down into the Planck range for values of mχ0 1 compatible with the LHC search limits and other constraints on these models. Instead, there has to be some extra mechanism for suppressing the LSP density. Examples include enhanced, rapid annihilation through direct-channel resonances such as Z , h, H/A. Another possibility is coannihilation with some other, almost-degenerate sparticle species [18,[159][160][161][162]: candidates for the coannihilating species identified in previous studies include thẽ τ 1 ,μ,ẽ,ν,t 1 andχ ± 1 . We introduced in [16] measures on the sparticle mass parameters that quantify the mass degeneracies relevant to the above-mentioned coannihilation and rapid annihilation processes, of which the following are relevant to our analysis of the SUSY SU(5) GUT model 9 : We also indicate above the colour codes used in subsequent figures to identify regions where each of these degeneracy conditions applies. We have verified in a previous study [16] that CMSSM, NUHM1 and NUHM2 points that satisfy the DM density constraint fulfill one or more of the massdegeneracy conditions, and that they identify correctly the mechanisms that yield the largest fractions of final states, which are usually 50% [14,169].
In much of the region satisfying theτ 1 degeneracy criterion above, theν τ has a similar mass, and can contribute to coannihilation [71]. We highlight the parts of the sample where sneutrino coannihilation is important by introducing a shading for regions where theν τ is the next-to-lightest sparticle (NLSP), and it obeys the degeneracy conditioñ We discuss later the importance of this supplementary DM mechanism.
As we discuss in this paper, a novel possibility in the SU(5) SUSY GUT is coannihilation with right-handed uptype squarks,ũ R andc R , which may be much lighter than the other squarks in this model, as a consequence of the freedom to have m 5 = m 10 . We quantify the relevant mass degeneracy criterion bỹ As we shall see in the subsequent figures, this novel degeneracy condition can play an important role when m 5 m 10 . The existence of this new coannihilation region was verified using SSARD [129], an independent code for calculating the supersymmetric spectrum and relic density.
We also distinguish in this analysis 'hybrid' regions where theτ 1 coannihilation and H/A funnel mechanisms may be relevant simultaneously: also with the indicated colour code.

Parameter planes
We display in Fig. 4 features of the global χ 2 function for the SUSY SU(5) GUT model in the (m 5 , m 1/2 ) plane (left panel) and the (m 10 , m 1/2 ) plane (right panel), profiled over the other model parameters. 10 Here and in subsequent parameter planes, the best-fit point is shown as a green star, the 68% CL regions are surrounded by red contours, and the 95% CL regions are surrounded by blue contours (as mentioned above, we use the χ 2 = 2.30 and χ 2 = 5.99 contours as proxies for the boundaries of the 68 and 95% CL regions in the fit). The regions inside the 95% CL contours are shaded Table 3 Parameters of the best-fit point in the SUSY SU(5) GUT model, with mass parameters given in GeV units. The numbers in parentheses in the bottom row are for a fit that does not include the LHC 13-GeV constraints and the recent PandaX-II and LUX constraints on DM scattering. Note that we use the same convention for the sign of A 0 as in [10][11][12][13][14][15][16], which is opposite to the convention used in, e.g., SoftSUSY, and that we use the notation sign(m 2 ) × |m 2 | → m for m 5 (11) according to the dominant DM mechanisms discussed in the previous section, see the criteria (4, 6, 7). In the (relatively limited) unshaded regions there is no single dominant DM mechanism.
As we see in Fig. 4, the best-fit point is at relatively small values of m 5 , m 10 and m 1/2 , close to the lower limit on m 1/2 , whereas the 68% CL region extends to much larger values of m 5 , m 10 and m 1/2 . The values of the model parameters at the best-fit point are listed in Table 3. 11 The upper row of numbers are the results from the current fit including the latest LHC 13-TeV and PandaX-II/LUX constraints, and the numbers in parentheses in the bottom row were obtained using instead the previous LHC 8-TeV and XENON100 constraints, but the same implementations of the other constraints. The most significant effect of the new LHC data has been to increase the best-fit value of m 1/2 by ∼160 GeV: the changes in the other fit parameters are not significant, in view of the uncertainties. As we discuss in more detail later, the favoured fit regions are driven by the (g − 2) μ constraint towards the boundary of the region excluded by the / E T constraint. Away from this boundary, the global χ 2 function is quite flat.
The best-fit point and much of the 68% CL region lie within the pink shaded region whereτ 1 −χ 0 1 coannihilation is the dominant DM mechanism. At larger values of m 5 and m 10 we encounter a blue shaded region where rapid annihilation via direct-channel H/A poles is dominant. We also see darker shaded hybrid regions whereτ 1 and H/A annihilation are important simultaneously. At larger values of m 1/2 3000 GeV, in the green shaded regions, the dominant DM mechanism isχ ± 1 −χ 0 1 coannihilation. There is also a band in the (m 10 , m 1/2 ) plane with m 10 1500 GeV and m 1/2 ∼ 1000 GeV, allowed at the 95% CL, whereν NLSP τ coannihilation is important.
We also note the appearance within the 95% CL region at m 1/2 ∼ 1000 GeV, and m 10 ∼ −1000 GeV of the novelũ R /c R −χ 0 1 coannihilation region (shaded yellow). 11 The SLHA files for the best-fit point and other supplementary material can be found in [17].
To understand the origin of this novelty, consider the oneloop renormalization-group equations for the states in the 10 representations of SU(5), namely (q L , u c L , e c L ) i , above the highest MSSM particle mass (all masses are understood to be scalar fermion masses, and we suppress subscripts L ): where t ≡ ln(Q/Q 0 ) with Q the renormalization scale and Q 0 some reference scale, and where the trace in S sums over the generations. Theũ R /c R − χ 0 1 coannihilation mechanism becomes important in a region of the SUSY SU(5) GUT parameter space where m 2 5 is very large and positive (∼27 TeV 2 ), m 2 10 is small and negative (∼−1.4 TeV 2 ), m 2 H u is very large and negative (∼−23 TeV 2 ), and m 2 H d is very large and positive (∼50 TeV 2 ). In this region, therefore, X t is very large and negative (∼ − 35 TeV 2 ), X b and X τ are suppressed because of small Yukawa couplings (tan β is not large in this region), and S is also very large and negative (∼ − 73 TeV 2 ), since m 2 H u − m 2 H d is large and negative and Tr(m 2 q −m 2 l −2m 2 u c +m 2 d c +m 2 e c ) vanishes at the GUT scale. Inspection shows that the X t terms in (8) and (9) drive the stop and sbottom masses upwards, and the S terms in (8) and (10) drive the left-handed squark and right-handed slectron masses upwards. On the other hand, the S term in (9) drives the right-handed squark masses downwards. Since there are no counteracting X terms for theũ R andc R , these have lower masses than the other sfermions, opening the way to aũ R /c R −χ 0 1 coannihilation region. 12 As discussed in more detail later, we used the Atom [144][145][146][147][148] simulation code for a dedicated verification that points in this region escape all the relevant LHC constraints. These points avoid exclusion by the LHC constraints through a combination of a strong mass degeneracy, mũ R /c R − mχ0 1 50 GeV, leading to strong suppression of the standard / E T signature, and the reduction of the production rate compared to the simplified model that assumes mass degeneracy of all eight light flavour squarks (see Fig. 1). These effects are clearly visible in Fig. 18 of [4]. Figure 5 displays the corresponding information in the (m 5 , m 10 ) plane of the SUSY SU(5) GUT model. As already reported in Table 3, here we see directly that the best-fit point has very small (and slightly negative) m 5 , and that m 10 is somewhat larger, exploiting the possibility that m 5 = m 10 , which is offered in this model. We also see again that the 68% CL region extends to values of m 5 and m 10 beyond thẽ τ 1 coannihilation region. We also note that in most of the rest of this planeχ ± 1 −χ 0 1 coannihilation is dominant, with only scattered regions where rapid H/A annihilation is important, even in combination withτ 1 coannihilation.
Projections of our results in the (tan β, m 1/2 ), (tan β, m 5 ) and (tan β, m 10 ) planes are shown in Fig. 6. We see that values of tan β 4 are allowed at the 95% CL, that the range tan β ∈ (8, 57) is favoured at the 68% CL, and that there is no phenomenological upper limit on tan β at the 95% CL. 13 The best-fit point has tan β = 13, as also reported in Table 3.
The pinkτ 1 −χ 0 1 coannihilation region is very prominent in the (tan β, m 1/2 ) projection shown in the upper panel of Fig. 6, as is the blue rapid H/A annihilation region and the purpleτ 1 −χ 0 1 coannihilation + H/A funnel hybrid region 13 The RGE evolution of the Yukawa couplings blows up for tan β 60.
at large tan β and m 1/2 . While the H/A funnel appears in the CMSSM only when tan β 45 for μ > 0 [18][19][20][21][22], in the SU(5) SUSY GUT model, it is found at significantly lower tan β, due to the separation of m H u and m H d from m 5 and m 10 , effectively making m A (and μ) free parameters as in the NUHM2. There is also a region in the (tan β, m 1/2 ) plane with tan β 10 and m 1/2 ∼ 1000 GeV whereν NLSP  [116]. Moreover, the Dark Matter mechanisms do not exhibit any preference for values of M h above or below the nominal central value. Thus, there is no apparent tension between this LHC measurement and the other constraints on the SUSY SU(5) GUT model, with the notable exception of (g − 2) μ .
As is well known, the calculation of M h in the MSSM is particularly sensitive to the value of the trilinear soft SUSYbreaking parameter A 0 as well as the stop squark masses. The latter depend in the SUSY SU(5) GUT model on m 10 and m 1/2 , but are insensitive to m 5 .
The (m H u , m H d ) plane is shown in Fig. 8. We see that the best-fit point lies in the quadrant where both m H u and m H d are negative, and that the 68% CL region extends also to the quadrant where m H d is negative and m H u is positive, as does theτ 1 −χ 0 1 coannihilation region. On the other hand, theχ ± 1 −χ 0 1 coannihilation region lies in the upper quadrants where m H d > 0. There is also an intermediate region, characterized by the H/A funnel mechanism and its hybridization withτ 1 coannihilation, part of which is also allowed at the 68% CL. There is also a region with The left panel in Fig. 9 displays the (M A , tan β) plane in the supersymmetric SU(5) GUT model. We see that M A 800(1000) GeV at the χ 2 = 5.99 (2.30) level, corresponding to the 95 (68) % CL, which is largely due to the interplay of the indirect constraints on (M A , tan β) such as M h (see also [135]) as well as the direct constraints from the LHC heavy MSSM Higgs searches. Even for large tan β, where these constraints impose the strongest lower  Fig. 4 limit on M A , it is much weaker than our global limit, which is M A 2800(>4000) GeV at the 95 (68) % CL. 14 We observe the same behaviour in the right panel of Fig. 9, where the one-dimensional likelihood profile for M A is shown. Indeed, the lightest pseudoscalar mass allowed at χ 2 = 4 is ∼920 GeV. The best-fit point in the global fit has (M A , tan β) (1600 GeV, 13): this is considerably beyond the present and projected LHC reach, though poorly determined. 14 We note, however, that this constraint might weaken with a larger parameter sample.

One-dimensional likelihood functions
We now discuss the one-dimensional χ 2 functions for various observable quantities.  of the SU(5) sample with m 5 /m 10 ∈ [0.9, 1.1], which we discuss in more detail later. 15 The current SU(5) fit exhibits minima of χ 2 at masses 2.5 TeV: mg 2600 GeV, common squark mass mq 2200 GeV, mũ R , md R , mt 1 2200 GeV and mτ 1 540 GeV, followed by a rise at higher mass towards a plateau with χ 2 2. The minimum is relatively sharp for mg, mq and mτ 1 , whereas it is broader for mt 1 . The exact values are listed in Table 4 and depicted in Fig. 11. In this figure we also indicate decay branching ratios (BRs) exceeding 20% by dashed lines, which are thicker for more important BRs. Figure 12 displays the one-dimensional 68 and 95% CL ranges 15 The χ 2 functions for the NUHM2 subsample are calculated relative to its minimum χ 2 , which is ∼0.4 higher than the minimum χ 2 for the full SU(5) sample.
for the Higgs and sparticle masses in the supersymmetric SU(5) model as darker and lighter coloured bands, with the best-fit values shown as blue lines.
Concerning future e + e − colliders, one can see that the best-fit masses of the lightest neutralino and stau are ∼ 500 GeV, and some other 68% CL ranges go down to 500 GeV, offering the possibility of pair production at a collider with √ s ∼ 1 TeV, as envisaged for the final stage of the ILC [172,173]. Going to higher centre-of-mass energies, e.g., √ s 3 TeV [173,174] as anticipated for CLIC, significant fractions of the 68% CL ranges of electroweak sparticle masses can be covered.
As already noted, a novel feature of the SUSY SU(5) GUT model with (m 5 = m 10 ) is that theũ R andc R may be much lighter than the other squarks. This leads to the possibility of aũ R /c R −χ 0 1 coannihilation strip where mũ R and mc R ∼ 500 GeV, which is visible as a second local minimum of χ 2 with χ 2 > 4 in the centre right panel of Fig. 10.
We have checked specifically whether this strip is allowed by the available LHC constraints. To this end, we verified using the Atom simulation code that points along this strip are consistent with the published constraints from the LHC 8-TeV data. We have also checked that this strip is consistent with the preliminary simplified-model search forqq +qq at 13 TeV reported by CMS. The left panel of Fig. 13 displays as a solid/dashed blue line the one-dimensional χ 2 function  Fig. 4 for mũ R − mχ0 1 including/omitting the 13-TeV data (the corresponding lines for mc R − mχ0 1 are very similar), and the right panel of Fig. 13 shows the region of the (mũ R , mχ0 1 ) plane where χ 2 < 5.99, i.e., allowed at the 95% CL. We find that σ (qq +qq) < 0.1 pb in this region, whereas the cross-section upper limit as given in [9] is 1 pb. We conclude that this simplified-model search does not affect the likelihood in thisũ R /c R −χ 0 1 coannihilation strip region. However, it will be explored further by future LHC data with increased luminosity.
Finally, we comment on the impact of the constraints from mono-jet searches [175][176][177]. These searches are designed to be sensitive to the highly compressed mass region by limiting the multiplicity of the high p T jets. In theũ R /c R −χ 0 1 coannihilation region, the mass difference is mildly compressed, mũ R /c R − mχ0 1 ∼ 40 GeV, and the jets fromũ R /c R decays are still resolvable from the background. Such extra jets will spoil the characteristic of the mono-jet event and reduce the efficiency. The degradation of the sensitivity for the mildly compressed region is clearly seen for example in Fig. 5 of [175]. For this reason, the mono-jet searches lose sensitivity to theũ R /c R −χ 0 1 coannihilation region, compared to the jets + / E T analysis [9], and we do not consider them in this paper.
Another novel feature of the SUSY SU(5) GUT model is visible in Table 4 and Fig. 11. Having m 5 = m 10 allows the possibility of strong mixing between theτ R in the 10 representation and theτ L in the5 representation. For example, at the best-fit point theτ 1 is an almost equal mixture ofτ L and τ R : This large mixing explains the level repulsion m 200 GeV between theτ 1 andτ 2 seen in Table 4, which is much larger than the splitting m 50 GeV between the almost unmixedẽ 1 ∼ẽ R andẽ 2 ∼ẽ L , which is also seen in Table 4.
We show in Fig. 14 the contribution to the global χ 2 function of (g − 2) μ (in teal), as a function of m 5 (left panel), m 10 (middle panel) and m 1/2 (right panel). In each case, there is a well-defined minimum that is lower than the plateau at large mass values by χ 2 2. In contrast, the contribu-tions to the global χ 2 function of the other observables are relatively featureless over large ranges of m 5 , m 10 and m 1/2 , with the exception of the contribution from the LHC 13-TeV data (mainly due to the / E T constraint), which rises sharply at low m 1/2 , as shown in the stacked red histogram in the right panel of Fig. 14. Because we profile over the other parameters, this does not have much impact on the dependence of χ 2 on m 5 and m 10 , as seen in the left and middle panels. The well-defined minima seen in the (g −2) μ contributions in the left and middle panels of Fig. 14 occur at quite small values of m 5 and m 10 , reflecting the fact that (g − 2) μ is sensitive to the soft symmetry-breaking contributions to the masses of bothμ L andμ R . These are m 5 and m 10 , respectively, so maximizing the SUSY contribution to (g − 2) μ and thereby minimizing the (g −2) μ contribution to χ 2 prefers small values of both m 5 and m 10 . Similarly, the SUSY contribution to (g − 2) μ is suppressed for large gaugino masses, explaining the aversion to large m 1/2 seen in the right panel of Fig. 14.
The principal contributions to the global χ 2 function at the best-fit point for the SUSY SU(5) GUT model are given in Table 5, and the corresponding pulls at the best-fit point are displayed graphically in Fig. 15. Apart from (g −2) μ , the other contributions deserving of comment include the following. The large contribution from HiggsSignals reflects the large number of channels considered, and has negligible variation for most of the points in our sample. We note that A FB (b) makes a contribution that is not much smaller than that of (g − 2) μ at the best-fit point, and that A e LR and σ 0 had also make relatively large contributions to the global χ 2 function. These observables reflect the residual tensions in the electroweak precision observables at the Z peak, which are present in the SM and the SUSY SU(5) GUT model is unable to mitigate.
In order to compare the quality of the SU(5) fit to the results of previous MasterCode analyses of competing models [15], we follow the prescription used there of subtracting from the total χ 2 given in Table 5 and Fig. 15, namely 100.34, the χ 2 contributions originating from Fig. 12 The one-dimensional 68 and 95% CL ranges of masses we obtain for the current fit in the supersymmetric SU(5) model, shown in dark and light orange, respectively. The best-fit point is represented by blue lines   Table 5 The principal χ 2 contributions of observables at the best-fit point in the SUSY SU(5) GUT model, together with the total χ 2 function 3.40 0.35 0.78 6.79 0.82 0.08  HiggsSignals [136,137], which dominate the global χ 2 function and would bias the analysis. Fig. 15 lists 36 separate contributions to the total χ 2 function. The first three (m t , M Z , and α (5) had (M Z )) are treated as nuisance parameters and the two LHC MET constraints at 8 and 13 TeV are applied as a single constraint. Omitting the HiggsSignals constraints in our determination of the number of degrees of freedom leaves 30 constraints, with seven parameters for the SU(5) model and hence 23 degrees of freedom.The χ 2 contributions from the relevant constraints sum to 32.39, corresponding to a χ 2 probability of 9%. This can be compared with the χ 2 probability values of 11, 12, 11 and 31% found in [15] for the CMSSM, NUHM1, NUHM2 and pMSSM10, respectively, using LHC Run 1 constraints. However, as in [15], we stress that these χ 2 probabilities are only approximate since, for example, they neglect correlations between the observables. A more complete treatment using toys, as done in the last reference of [37][38][39][40][41][42][43][44][45], is beyond the scope of this work.
There are a couple of important corollaries to this observation, one concerning mt 1 . It is sensitive to A 0 as well as the soft SUSY-breaking contributions to thet L andt R mass parameters (which are both given by m 10 in the SUSY SU(5) GUT model). Since A 0 is relatively poorly determined, the χ 2 minimum for mt 1 is relatively shallow, as seen in the lower left panel of Fig. 10.
The second observation concerns the sign of μ. All our analysis has been for μ > 0, which is the sign capable of mitigating the discrepancy between the experimental value of (g − 2) μ and the SM prediction. For μ < 0, the largemass plateau would have a similar height as in Fig. 14, but the χ 2 function would rise monotonically at low values of m 5 , m 10 and m 1/2 , instead of featuring a dip. Thus, the μ < 0 possibility would be disfavoured by χ 2 2, and the global minimum would lie at large masses and be ill defined.
The χ 2 distributions for some more observables are shown in Fig. 16, We see that the minima for mχ0 1 (upper left panel) and mχ± 1 (upper right panel) are quite well defined, mirroring the structure in the χ 2 function for mτ 1 shown in the lower right panel of Fig. 10. The preference for a (very) smallτ 1 − χ 0 1 mass difference is seen in the lower left panel of Fig. 16, and reflects the fact, commented on in connection with many previous figures, that the best-fit point and much of the 68% CL region lies in theτ 1 −χ 0 1 coannihilation region. On the other hand, a small mt 1 −mχ0 1 mass difference is disfavoured, as seen in the lower right plot of Fig. 16, reflecting the fact that stop coannihiliation does not play a significant role.
Theχ ± 1 −χ 0 1 coannihilation region is prominent in the previous figures, and also contains parameter sets that are preferred at the 68% CL. Hence a smallχ ± 1 −χ 0 1 mass difference is also allowed at the χ 2 1 level, as seen in the left panel of Fig. 17, although the best-fit point has mχ± 1 − mχ0 1 ∼ 470 GeV. However, values of theχ ± 1 lifetime that are allowed at the 95% CL are all too short to provide a long-lived particle signal, as seen in the right panel of Fig. 17. 16 We now discuss the one-dimensional likelihood functions for electroweak precision observables and observables in the flavour sector shown in Fig. 18, together with the current experimental measurements and their uncertainties shown as dotted grey lines. The upper left panel displays (g − 2) μ , and we see that the global minimum occurs for (g−2) μ 0.4× 10 −9 , with χ 2 −2 compared to the case (g − 2) μ = 0. We see again that the SUSY SU(5) GUT model is able to mitigate slightly the discrepancy between the SM and the measurement of (g − 2) μ , although it does not provide a substantial improvement over the SM prediction.
As for M h , as shown in the upper right panel of Fig. 18 the χ 2 function is minimized close to the nominal experimental value, and is quite symmetric, showing no indication of any tension in the SUSY SU(5) GUT model fit. Likewise, the best-fit value of M W (lower left panel of Fig. 18) is highly compatible with the experimental measurement, and that for BR(B s,d → μ + μ − ) (lower right panel) is very close to the SM prediction, and hence also compatible with the experimental measurement. We note that, whereas values of BR(B s,d → μ + μ − ) that are slightly larger than the SM value are possible, smaller values are strongly disfavoured in the SUSY SU(5) GUT model.

Higgs branching ratios
We present in Fig. 19 the one-dimensional likelihood functions for the ratios of supersymmetric SU(5) and SM predictions for the BRs of h → γ γ (left panel), h → Z Z * (middle panel) 17 and h → gg decays (right panel). We see that in   each case the preferred region in the fit corresponds to a prediction in the SU(5) model that deviates from the SM case by at most a few %, whereas the present experimental uncertainties in the different coupling modifiers squared (employing some theory assumptions) are typically O(30)% [179], and a precision of O(5−10%) (with the same theory assumptions) can be reached by the end of the LHC programme. On the other hand, future e + e − colliders such as the ILC, CLIC or FCC-ee anticipate a precision at the percent level for couplings to fermions and at the permille level for couplings to massive gauge bosons [173,180]. This offers the possiblity that deviations from the SM in the SUSY SU(5) GUT model can be measured in the future.

Comparison with previous results
In previous papers we have studied the CMSSM, NUHM1 and NUHM2 using the LHC 8-TeV results and earlier DM scattering constraints. None of these models are directly comparable to the supersymmetric SU(5) model studied here, which has four different soft SUSY-beaking scalar mass parameters, m 5 , m 10 , m H u and m H d . The most similar is the NUHM2, which has the three parameters m 0 = m 5 = m 10 , m H u and m H d . Here we compare the supersymmetric SU(5) results found in this paper using LHC 13-TeV data with 'fake' NUHM2 results obtained by selecting a subset of this SU(5) sample with m 5 /m 10 ∈ [0.9, 1.1] (which were also displayed as grey lines in Fig. 10) and with previous NUHM2 results [14]. Figure 20 compares the one-dimensional χ 2 likelihood functions for mg (upper left), mq R (upper right), mt 1 (lower left) and mτ 1 (lower right) found in the SU(5) model including LHC 13-TeV constraints (solid blue lines) with the restricted fake NUHM2 sample (solid grey lines) and, for comparison, results from our previous NUHM2 analysis that used only the LHC 7-and 8-TeV constraints (dashed grey lines) [14].
We see here and in Fig. 10 that the restricted 'fake' NUHM2 sample exhibits, in general, best-fit masses that are similar to those found in the full SU(5) sample. The most noticeable differences are that lower masses are disfavoured in the restricted sample relative those in the full SU(5) model, indicating that the latter has some limited ability to relax the NUHM2 lower bounds on sparticle masses, e.g., at the 95% CL. The previous NUHM2 analysis [14] also yielded similar best-fit masses but, as could be expected, gave 95% CL lower limits on sparticle masses that were further relaxed. Similar features can also be observed in Figs. 16, 17, 18 and 19, where we have also included the 'fake' NUHM2 subsample.
Restricting further our SU(5) to mimic the NUHM1, let alone the CMSSM, is not useful because of the increased sampling uncertainties in such restricted samples. However, we showed in [14] that our NUHM2 LHC 7-and 8-TeV results for the exhibited sparticle masses were broadly similar to those for the NUHM1 and the CMSSM [13], and we expect the impacts of the LHC 13-TeV data on these models to be comparable to that in the NUHM2.
Finally, we ask whether or not there is a significant improvement in the SU(5) fit compared to that in the NUHM2 subsample, thanks to the additional parameter (m 5 and m 10 replacing m 0 ). The NUHM2 subsample has a total χ 2 = 100.8, which is reduced to 32.8 when we remove the contributions from HiggsSignals, as discussed earlier. It should be noted that the NUHM2 subsample is statistically significantly smaller than that of the SU(5) sample. The quoted NUHM2 χ 2 represents only an upper bound on the χ 2 of the best-fit point that would be found in a more complete sample of the NUHM2. Since the NUHM2 model has one less parameter than the SU(5) model, it has 24 degrees of freedom, and its χ 2 probability is 11%. According to the Wilks test [181], the probability that the data are represented better by the SU(5) model than by the NUHM2 subsample is 50%, while the F test [182] yields a 40% probability. Therefore we conclude that there is no evidence that the extra parameter of SU(5) provides a significant improvement.

The possibility of a long-livedτ 1
The possibility of a very smallτ 1 −χ 0 1 mass difference opens up the possibility that theτ 1 might have a long lifetime, as discussed in the contexts of the CMSSM, NUHM1 and NUHM2 in [16]. This would occur if mτ 1 − mχ0 1 < m τ . As seen in the lower left panel of Fig. 16, the best-fit point has a mass difference ∼20 GeV, outside this range, but mτ 1 − mχ0 1 < m τ is allowed with χ 2 ∼ 1. In Fig. 21 we analyze the lifetime of theτ 1 . We see in the upper left panel of Fig. 21 that there is essentially no χ 2 penalty for 10 −9 s ττ 1 10 −2 s, with lifetimes ∼10 −10 s and 10 3 s allowed with χ 2 1. Distinguishing a separated-vertex signature at the LHC would be challenging for smaller values of ττ 1 , and there would be significant disruption of the successful conventional Big Bang nucleosynthesis calculations for ττ 1 10 3 s [183][184][185][186][187][188][189][190][191][192].
The upper right plot of Fig. 21 compares theτ 1 lifetime with its mass. The plane is characterized by a strip with 800 GeV τ 1 1200 GeV allowed at the 68% CL, while the 95% CL region is significantly wider, ranging from mτ 1 ∼ 500 GeV to mτ 1 ∼ 2000 GeV.
The lower panels of Fig. 21 display the regions of the (m 5 , m 1/2 ) (left) and (m 10 , m 1/2 ) (right) planes in the SUSY SU(5) GUT model where the lowest-χ 2 points have 10 −10 s < ττ 1 < 10 3 s. The colour-coding indicates the lifetimes of these points, as indicated in the legends. The contours for χ 2 < 2.30(5.99) relative to the best-fit point in our sample are shown as solid red and blue lines, respectively. One can see that larger lifetimes occur all over the displayed parameter space, with a slight preference for larger m 5 or m 10 values. Lower panels the (m 5 , m 1/2 ) and (m 10 , m 1/2 ) planes, coloured according to the values of ττ 1 . The 68 and 95% CL contours in these three planes are coloured red and blue, respectively 10 Direct dark matter detection As already mentioned, the PandaX-II experiment [83] has recently published results from its first 98.7 days of data, which currently provide the most stringent upper limits on the spin-independent DM scattering cross section on protons, σ SI p . In parallel, the LUX Collaboration [84] has presented preliminary constraints on σ SI p from 332 days of data. We have combined these two constraints on σ SI p into a single experimental likelihood function, which we have then convoluted with an estimate of the theoretical uncertainty in the calculation of σ SI p , as described in [16], to constrain the SUSY SU(5) GUT parameter space. This constraint has been used in obtaining the global fit whose results we have presented in the previous Sections. Here we discuss the future prospects for direct DM detection in light of our global fit. , σ SI p ) plane seen in Fig. 22. We see that there are regions favoured at the 68% CL that lie relatively close to this boundary, whereas the main 68% CL region and the best-fit point have smaller values of σ SI p . We also note that the H/A funnel andχ ± 1 −χ 0 1 DM mechanisms favour values of σ SI p that are relatively close to the PandaX-II/LUX boundary, whereas theτ 1 −χ 0 1 mechanism and its hybridization with the H/A funnel favour smaller values of σ SI p . The upcoming XENON1T [193] experiment will be able to probe the wholẽ χ ± 1 coannihilation region and a substantial part of the H/A funnel region.
We also display in Fig. 22 the projected 95% exclusion sensitivity of the future LUX-Zeplin (LZ) and XENONnT experiments (solid purple and dashed blue lines respec- , σ SI p ) plane in the SUSY SU(5) GUT model. The solid green line is the 95% CL upper limit from the XENON100 experiment, and the dashed black solid line is the new 95% CL upper limit from the LUX experiment. The solid black line shows the 95% CL exclusion contour for our combination of the PandaX-II and LUX experiments, the solid purple line shows the projected 95% exclusion sensitivity of the LUX-Zeplin (LZ) experiment, the solid and dashed blue lines show the projected 95% sensitivities of the XENON1T and XENONnT experiments, respectively, and the dashed orange line shows the astrophysical neutrino 'floor', below which astrophysical neutrino backgrounds dominate (yellow region). The other line colours and shadings within the 68 and 95% CL regions are the same as in Fig. 4 tively) [193,194], and the astrophysical neutrino 'floor' (dashed orange line) [195,196], below which astrophysical neutrino backgrounds dominate (yellow region). We see that much of theτ 1 −χ 0 1 coannihilation region and the region of its hybridization with the H/A funnel lie below the projected sensitivities of the LZ and XENONnT experiments, and substantial portions of them also lie below the neutrino 'floor'. On the bright side, however, we recall that theτ 1 −χ 0 1 region, in particular, lies at relatively small values of m 5 , m 10 and m 1/2 , offering greater prospects for detection at the LHC than, e.g., theχ ± 1 −χ 0 1 region, so there is complementarity in the prospects of the LHC and direct DM experiments for probing the SUSY SU(5) GUT model, as was noted previously for other SUSY models [16].

Summary and conclusions
We have explored in this paper the experimental, phenomenological, astrophysical and cosmological constraints on the minimal SUSY SU(5) GUT model. In this scenario the GUT-scale universal soft SUSY-breaking scalar mass m 0 is replaced by independent masses for the 10 and5 sfermions. This flexibility introduces some features that are novel compared to the GUT-universal CMSSM, NUHM1 and NUHM2.
In general we observe that many best-fit values of the coloured particles are within the reach of the HL-LHC, but that the preferred regions clearly extend beyond the reach of the final stage of the LHC. On the other hand, the bestfit masses of some electroweakly interacting particles are ∼500 GeV, offering the possibility of pair production at a collider with √ s ∼ 1 TeV, as envisaged for the final stage of the ILC. Going to higher centre-of-mass energies, √ s 3 TeV as anticipated for CLIC, significant fractions of the 68% CL ranges of electroweak sparticle masses can be covered.
One novelty is the appearance of aũ R /c R −χ 0 1 coannihilation region that appears where m 2 5 is large and positive, m 2 10 is small and negative, and m 2 H u and m 2 H d are large and negative. On the other hand, we find thatt 1 −χ 0 1 coannihilation is not important in the SUSY SU(5) GUT model, nor are the focus-point region and rapidχ 0 1χ 0 1 annihilation via direct-channel h and Z poles. We have checked that theũ R /c R −χ 0 1 coannihilation region is not yet excluded by searches for / E T events at the LHC, because the production rate is reduced compared to the case where all eight squarks are mass degenerate and the smallũ R /c R −χ 0 1 mass difference suppresses this signature. However, this region may be accessible with future LHC runs.
We have also highlighted the possibility that aν τ NLSP might have an important coannihilation role. Another novelty is the composition of theτ 1 NLSP in a significant region of the model parameter space. In the GUT-universal CMSSM, NUHM1 and NUHM2 models, the universality of m 0 and the greater renormalization for SU(2) doublets impose a substantial mass difference between theτ 2 and theτ 1 , with the latter being predominantly aτ R . However, in the SUSY SU(5) GUT model with m 5 = m 10 , theτ R andτ L may have similar masses, and the off-diagonal entries in theτ mass matrix may cause large mixing and repulsion between theτ 1 andτ 2 masses.
On the other hand, one experimental signature that is shared by the SUSY SU(5) GUT model and GUTuniversal models is the possible appearance of a long-lived (metastable)τ 1 . This is a feature of a significant fraction (but not all) of theτ 1 −χ 0 1 coannihilation region. The prospects for direct DM detection are mixed: they are relatively good in theχ ± 1 −χ 0 1 coannihilation region, but less promising in the rapid H/A annihilation and hybrid regions, though potentially detectable in the planned LUX-Zeplin experiment. On the other hand, theτ 1 −χ 0 1 coannihilation region probably lies beyond the reach of this experiment, as does part of the hybrid region. Indeed, portions of these regions lie below the neutrino 'floor'. On the other hand, sub-stantial parts of these regions are accessible to LHC searches for long-lived particles and / E T .