The NUHM2 after LHC Run 1

We make a frequentist analysis of the parameter space of the NUHM2, in which the soft supersymmetry (SUSY)-breaking contributions to the masses of the two Higgs multiplets, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m^2_{H_{u,d}}$$\end{document}mHu,d2, vary independently from the universal soft SUSY-breaking contributions \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m^2_0$$\end{document}m02 to the masses of squarks and sleptons. Our analysis uses the MultiNest sampling algorithm with over \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4 \times 10^8$$\end{document}4×108 points to sample the NUHM2 parameter space. It includes the ATLAS and CMS Higgs mass measurements as well as the ATLAS search for supersymmetric jets + \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${/\!\!E}_T$$\end{document}/ET signals using the full LHC Run 1 data, the measurements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{BR}(B_s \rightarrow \mu ^+\mu ^-)$$\end{document}BR(Bs→μ+μ-) by LHCb and CMS together with other B-physics observables, electroweak precision observables and the XENON100 and LUX searches for spin-independent dark-matter scattering. We find that the preferred regions of the NUHM2 parameter space have negative SUSY-breaking scalar masses squared at the GUT scale for squarks and sleptons, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_0^2 < 0$$\end{document}m02<0, as well as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m^2_{H_u} < m^2_{H_d} < 0$$\end{document}mHu2<mHd2<0. The tension present in the CMSSM and NUHM1 between the supersymmetric interpretation of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(g-2)_\mu $$\end{document}(g-2)μ and the absence to date of SUSY at the LHC is not significantly alleviated in the NUHM2. We find that the minimum \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2 = 32.5$$\end{document}χ2=32.5 with 21 degrees of freedom (dof) in the NUHM2, to be compared with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2/\mathrm{dof} = 35.0/23$$\end{document}χ2/dof=35.0/23 in the CMSSM, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2/\mathrm{dof} = 32.7/22$$\end{document}χ2/dof=32.7/22 in the NUHM1. We find that the one-dimensional likelihood functions for sparticle masses and other observables are similar to those found previously in the CMSSM and NUHM1.


Introduction
Supersymmetric (SUSY) models are among the bestmotivated extensions of the Standard Model (SM) that might be discovered at the Large Hadron Collider (LHC). They stabilise the electroweak hierarchy [1,2] and facilitate grand unification [3][4][5][6][7], and the lightest supersymmetric particle (LSP) provides a natural candidate for the cosmological dark matter [8,9]. However, the absence of a signal in direct searches for SUSY particles in Run 1 of the LHC [10,11] sets strong constraints on supersymmetric models, as do the measurement of the mass and properties of the Higgs boson [12,13] and precision measurements of rare decays such as B s → μ + μ − [14][15][16].
Gaining a fully accurate picture of the effects of these constraints requires that they be combined in global statistical fits within specific supersymmetric models. Particularly well-motivated and simplified versions of the minimal supersymmetric Standard Model (MSSM) [17,18] are derived from grand unified theory (GUT) model-building considerations. There have been a number of analyses  of the constraints imposed by LHC Run 1 data on the parameter spaces of such models, particularly the constrained MSSM (CMSSM) [43][44][45][46][47][48][49][50][51][52][53][54][55][56][57][58][59][60], whose parameters are the soft supersymmetry (SUSY)-breaking masses m 0 , m 1/2 and A 0 that are universal at the GUT scale, and tan β, the ratio of the two vacuum expectation values of the two Higgs doublets. There have also been some studies of the LHC constraints on the NUHM1 [61][62][63][64], in which the soft SUSY-breaking contri-butions to the masses of the electroweak Higgs multiplets, m 2 H d ,H u , are equal but non-universal. However, these models have become very constrained by the recent data. The anomalous magnetic moment of the muon (g − 2) μ [65][66][67][68][69][70][71][72][73][74] is a particular source of tension, as has been reinforced by the recent convergence in the Standard Model (SM) calculations of (g − 2) μ based on τ decays and different sets of e + e − annihilation data [75,76]. As is well known, the ∼3.5 σ discrepancy between the observed value and SM prediction can be reduced by SUSY contributions due to relatively light electroweakly interacting superpartners. In the simple GUT-based models mentioned above, direct searches and the Higgs mass force the coloured superpartners to be so heavy that, due to the universality of the soft SUSY-breaking parameters m 0 and m 1/2 at the GUT scale that leads also to relatively heavy electroweak superpartners, these models cannot remove the (g −2) μ anomaly [77]. Also for this reason, the LHC searches for leptons and electroweak inos do not impact significantly the parameter spaces of these GUT-based models.
A related extension of these models which a priori might be able to alleviate this tension is the NUHM2 [78,79], in which m 2 H d = m 2 H u = m 2 0 in general, 1 but the soft SUSYbreaking parameters m 0 , m 1/2 and A 0 are still universal at the GUT scale. An equivalent formulation of the NUHM2 is to treat the pseudoscalar mass M A and supersymmetric Higgs mass term μ as free parameters, which could lead to interesting phenomenology associated with light higgsinos and/or a light pseudoscalar Higgs. Moreover, new terms in the renormalisation group equations (RGEs) associated with the scalar-mass non-universality in the NUHM2 may lead to lighter left-handed sleptons, offering further avenues for ameliorating the tension with (g − 2) μ (see Sect. 3.1 for details).
We find that the NUHM2, despite its freedom in the choices of M A and μ, is unable to alleviate significantly the tension between the absence to date of SUSY at the LHC and the supersymmetric interpretation of (g − 2) μ that had been found previously in the CMSSM and NUHM1. We find that the minimum χ 2 /dof = 32.5/21 in the NUHM2, to be compared with χ 2 /dof = 35.0/23 in the CMSSM and χ 2 /dof = 32.7/22 in the NUHM1. A novel feature of the best NUHM2 fit is that the preferred regions of the NUHM2 parameter space have negative SUSY-breaking scalar masses squared for squarks and sleptons, m 2 0 < 0, as well as It is quite possible that SUSY-breaking scalar masses are negative at the GUT scale and yet, when run down to the weak scale, there are no tachyonic scalars in the theory. This extension to the CMSSM parameter space in the context of a gravitino LSP was considered in [85,86], in mirage-mediated models in [87], and in recent post-Higgs gauge mediation constructions [88]. These models as well as NUHM2 models with m 2 H u < 0 and m 2 H d < 0 are potentially problematic due to the presence of charge-and colour-breaking minima, particularly along F-and D-flat directions [89,90]. However, so long as the standard electroweak vacuum is long-lived, the relevance of other vacua becomes a cosmological question related to our position in field space after inflation. For a discussion of cosmological issues associated with such tachyonic soft SUSY-breaking mass parameters; see [91].
As an output of our analysis, we compare the onedimensional likelihood functions for sparticle masses and other observables in the NUHM2 with those found previously in the CMSSM and NUHM1. The 95 % CL lower limits on the gluino, squark, stop and stau masses are not very different in the NUHM2 from those found previously in the CMSSM and NUHM1. However, the distinction found in those models between low-and high-mass regions of their respective parameter spaces is largely lost in the NUHM2 because of its greater flexibility in satisfying the dark-matter constraint. In addition to sparticle masses, we also present NUHM2 predictions for BR(B s → μ + μ − ) and the spinindependent dark-matter scattering cross section, σ SI p .

Analysis procedure
We follow closely the procedure described in [19]. Our treatment of the non-LHC constraints is identical with the treatment in [19], and we treat the top quark mass and the strong coupling as nuisance parameters with Gaussian priors: m t = 173.2 ± 0.9 GeV and α s (M Z ) = 0.1185 ± 0.0006. We again use the MultiNest algorithm to sample the NUHM2 parameter space, just as we did previously for the CMSSM and NUHM1 models. As mentioned in the Introduction, we use a NUHM2 sample comprising ∼4 × 10 8 points, with the aim of sampling adequately features of the six-dimensional NUHM2 parameter space , tan β}, ensuring in particular that all high-likelihood regions are identified and well characterised. We sample the ranges −1333 GeV < m 0 < 4000 GeV, 0 < m 1/2 < 4000 GeV, −5 × 10 7 GeV 2 < m 2 H u , m 2 H d < 5 × 10 7 GeV 2 , −8000 GeV < A 0 < 8000 GeV and 2 < tan β < 68. (Here and subsequently, negative values of m 0 should be understood as m 0 ≡ Sign(m 2 0 ) |m 2 0 | < 0, and we use analogous definitions for negative values of m H u and m H d .) The parameter ranges are scanned by dividing the range of m 0 into 4 segments, and the ranges of m 1/2 , m H u and m H d into 3 segments each, yielding a total of 108 boxes. Their boundaries are smeared using a Gaussian function so as to sample the NUHM2 parameter space smoothly, which also provides some information beyond the nominal sampling range, as we discuss later in the case of m H u and m H d .
We merge this dedicated sample of the NUHM2 parameter space with the samples of the CMSSM and NUHM1 parameter spaces used in [19]. The latter are subspaces of the full NUHM2 parameter space, and the CMSSM and NUHM1 points provide supplementary sampling of the likelihood function of the NUHM2.
We construct a global likelihood function that receives contributions from the usual electroweak precision observables, as well as B-decay measurements such as BR(b → sγ ), BR(B u → τ ν τ ) and BR(B s → μ + μ − ). Bounds on their experimental values as well as those on the cosmological dark matter density, the cross section for spin-independent dark-matter scattering from the LUX experiment and the LHC search for supersymmetric signals are given in [77], with updates detailed in [92]. The observables we use, as well as the values and errors we assume, are given in Table 1, with references to their sources.
into FeynHiggs [97][98][99][100][101] yield an upward shift of M h for scalar top masses in the (multi-)TeV range and reduce the theoretical uncertainty in the Higgs mass calculation [139], which is nevertheless significantly larger than the variations in the best-fit Higgs mass since its discovery and the differences between the values reported by ATLAS and CMS. We incorporate here the public results of the search for jets + / E T events without leptons using the full ATLAS Run 1 data set of ∼20/fb at 8 TeV [10], which has greater sensitivity to the relevant parts of the NUHM2 parameter space than searches including leptons and/or b quarks. 5 Experimental searches for jets + / E T events are typically analysed within the framework of the CMSSM for some fixed A 0 and tan β. The applicability of these analyses to other A 0 and tan β values, as well as to constraining the NUHM1,2, requires some study and justification. One issue is that, for any specific set of values of m 0 , m 1/2 , A 0 and tan β, the sensitivities of ATLAS and CMS to jets + / E T events might depend on the degree of non-universality in the NUHM1,2. A second issue is that the range of m 0 in the NUHM2 that is consistent with theχ 0 1 LSP requirement depends on the degrees of non-universality. Specifically, this requirement is compatible with m 2 0 < 0 in the NUHM2, a possibility that is absent for the CMSSM, but can occur in the NUHM1 for m 1/2 2000 GeV when m 2 H d = m 2 H u < 0 and dominates over m 2 0 in the renormalisation-group evolution. In the NUHM2 it is even easier to obtain m 2 0 < 0 and remain compatible with a neutralino LSP, because a combination of soft supersymmetry-breaking parameters known as S (defined below) may be non-zero.
Since the ATLAS experiment quotes limits only for the CMSSM with m 2 0 > 0, we rely on a previous dedicated study of the jets + / E T search at 7 TeV [92], made using the Delphes [140] generic simulation package with a 'card' to emulate the performance of the ATLAS detector, that showed that the LHC results could be extrapolated to m 2 0 < 0. As shown in Figs. 2, 3 and 4 of [92], his study confirmed that / E T constraints in the (m 0 , m 1/2 ) plane of the CMSSM are relatively insensitive to tan β and A 0 , as stated in [141], and that the / E T constraints are also quite insensitive to the degrees of non-universality in the NUHM1,2, with any variations in the sensitivity being smaller than the uncertainties in our simulation. Specifically, it was found that the 95 % CL bounds in the (m 0 , m 1/2 ) plane of the CMSSM were approximately independent of A 0 and tan β, as also stated by CMS [142]; the same was true for m 2 The same is expected to be true for the 8-TeV ATLAS jets + / E T search 5 Searches with leptons and/or b quarks are important at large m 0 , but they extend only to values of m 1/2 that are disfavoured already by other constraints.
[10] used here, which uses a similar event selection to the ATLAS 7-TeV data studied in [92]. Finally, we also incorporate here the most recent constraints on A/H production from ATLAS and CMS [124,125], using the same approach as in [19].

Scalar-mass parameters and renormalisation
Before discussing our results for the NUHM2, we briefly review another important difference between this model and its more constrained relatives. When is non-zero. In both the CMSSM and NUHM1, S = 0 and is a fixed point of the RGEs at the one-loop level and remains zero at any scale [144,145]. However, in the NUHM2, with m 2 H u = m 2 H d , S = 0 at the GUT scale, as seen in (1), which can cause the low-energy spectrum to differ significantly from that in the CMSSM or NUHM1. For example, consider the renormalisation-group equation for the τ R mass: When S < 0, the evolution of m 2 τ R receives a positive contribution as it runs down from the GUT scale to the electroweak scale. As a result, ensuring a neutralino LSP becomes a generic possibility even when m 2 0 < 0. 6 Furthermore, the masses of left-handed sleptons may run to lighter values than their right-handed counterparts, allowing for new coannihilation channels to regulate the neutralino relic density [78,79], or larger contributions to (g − 2) μ .

Model parameter planes
The (m 0 , m 1/2 ) plane: We first present results for the (m 0 , m 1/2 ) plane shown in Fig. 1. We denote the best-fit point by a filled green star and the χ 2 = 2.30 and 5.99 contours by solid red and blue contours, respectively. These would correspond to 68 and 95 % CL contours if the errors were Gaussian. In the upper 6 In the NUHM1, the flexibility to allow m 2 H d = m 2 Hu < 0 with a different value from m 2 0 can also affect the running to ensure a neutralino LSP when m 2 0 < 0, but only in a restricted region when m 1/2 is large enough.  Fig. 1 we also show the best-fit points in the NUHM1 and CMSSM (shaded and open green stars), and the 68 and 95 % CL contours in these models (dashed and dotted red and blue contours, respectively). It is apparent from Fig. 1 that the extents of these contours that the χ 2 function for the NUHM2 is quite shallow, and we emphasise that the best-fit point and other details of the χ 2 function should not be over-interpreted.
We see that the 68 % CL NUHM2 region in the upper left panel of Fig. 1 extends in a lobe down to m 1/2 ∼ 300 to 2000 GeV for −500 GeV m 0 2000 GeV, whereas m 0 is relatively unrestricted for m 1/2 2500 GeV. At the 95 % CL we find m 1/2 500 GeV for m 0 0. The best-fit point in the NUHM2 has m 0 ∼ −500 GeV and m 1/2 ∼ 1800 GeV. The LHC / E T search with the most impact on the parameter space is that with jets and zero leptons, which constrains the NUHM2 parameter space most when m 0 1500 GeV. As already mentioned, we have verified previously [92] that this constraint is approximately independent of the other NUHM2 parameters in the (m 0 , m 1/2 ) region of interest. Searches for events with b-jets and/or leptons have greater sensitivity when m 0 1500 GeV, but are important only outside the 95 % CL contour, at lower m 1/2 , so we have not studied in detail their sensitivity to the model parameters.
In the case of the NUHM1, the range of m 0 where low values of m 1/2 2000 GeV are allowed at the 68 % CL (within the dashed red contour in Fig. 1) is much smaller, being limited to 200 GeV m 0 1000 GeV. The case of the CMSSM is much more restrictive, with only a small part of the 68 % CL region (within the dotted red contour in Fig. 1) with 300 GeV m 0 1300 GeV appearing when m 1/2 1800 GeV. Moreover, this case has a largest allowed value of m 1/2 ∼ 2500 GeV at the 95 % CL, whereas we observe no upper bound on m 1/2 in either the NUHM1 or the NUHM2.

The dark-matter constraint:
The dark-matter density constraint is less restrictive in the NUHM2 than in the NUHM1 and, particularly, the CMSSM. In the regions of interest, the dark-matter density is generally brought down into the range allowed by cosmology through enhancement of (co-)annihilation processes due to particular properties of the spectrum. In the other panels of Fig. 1 we use different colours of shading to visualise the impacts of these processes, by displaying areas of the 95 % CL regions in the (m 0 , m 1/2 ) planes where the following conditions are satisfied: χ ± 1 coannihilation (green): each of which is surrounded by a dotted contour. Regions where more than one of these conditions are satisfied are shaded in darker colours. We have also explored the focuspoint [146][147][148] criterion |μ/mχ0 1 − 1| < 0.3, and found that it is not relevant in the displayed portions of the (m 0 , m 1/2 ) planes. We note that the criteria above are approximate, being intended only to serve as guides to the different regions in the (m 0 , m 1/2 ) planes.
We see in the upper right panel of Fig. 1 that the low-mass region of the CMSSM is in the stau coannihilation region [149][150][151][152][153][154][155][156][157] (pink shading) and its high-m 0 region (blue shading) is in the funnel region where the LSPs annihilate rapidly through the s-channel heavy Higgs resonances A/H [43][44][45][46][47]. The best-fit CMSSM point now lies in the stau coannihilation region: the difference from the low-mass best-fit point found in [19] is due to using the updated ATLAS jets + / E T constraint [10]. The current CMSSM best-fit point is very similar to the previous local best fit in the low-mass region. We also see for 1000 GeV m 0 2000 GeV and m 1/2 2000 GeV (shaded purple) a CMSSM region where both the stau-coannihilation and funnel criteria are satisfied.
In the NUHM1, as seen in the lower left panel of Fig. 1 it is possible to satisfy the χ h 2 constraint for larger values of m 1/2 than are possible in the CMSSM, thanks to the extra degree of freedom associated with the soft SUSYbreaking contribution to the Higgs masses. In the low-mass NUHM1 region, the relic density is again determined by stau coannihilation (pink shading), whereas at large m 0 and m 1/2 2500 GeV the rapid annihilation via the A/H fun-nel (blue shading) is important. The NUHM1 best-fit point is in a high-mass region where χ h 2 is determined by coannihilations of nearly degenerateχ 0 1 ,χ ± 1 andχ 0 2 [157-161] (green shading), since μ m 1/2 and the LSP is nearly a pure higgsino. There is also a region where both stau andχ ± 1 coannihilations are important (dark green shading).
In the case of the NUHM2, all four of the mechanisms (3) come into play, as we see in the lower right panel of Fig. 1. As in the cases of the CMSSM and NUHM1, there are regions where stau coannihilation (pink), rapid annihilation via A/H bosons (blue) andχ ± 1 coannihilations (green) are important, as well as a region where both stau andχ ± 1 coannihilations are important (dark green). We also see two small bands with (m 0 , m 1/2 ) ∼ (2000, 600) GeV where stop coannihilation [162][163][164][165][166][167] is important.
Our best-fit point for the NUHM2 has m 2 0 < 0 in the pink region where the relic density is fixed by stau coannihilation. As can be seen in Fig. 2, the LSP and the lighter stau are indeed very nearly degenerate at this point, with the other sleptons only slightly heavier but the other sparticles significantly more massive. Also, M A 2mχ0 1 , so there is no significant enhancement of LSP annihilations via direct-channel resonances. We emphasise, however, that the NUHM2 spectrum is poorly determined, and that this and other processes play important roles in other parts of the NUHM2 parameter space. We find M h = 124.8 GeV at the best-fit point. For comparison, the lower panels of Fig. 2 display the best-fit spectra in the CMSSM (left) and the NUHM1 (right). In the case of the CMSSM, the best-fit point is also in the stau coannihilation region, whereas the best NUHM1 fit is in a region where both stau andχ ± 1 coannihilations are important.
Other parameter planes: Figure 3 displays the (m 0 , tan β) plane (left) and the (tan β, m 1/2 ) plane (right) in the NUHM2, CMSSM and NUHM1. In both panels, we see that a large range 5 tan β 60 is allowed at the 68 % CL (solid red contour). 7 The range of tan β within the 68 % CL region is restricted to values 40 for the lower-mass lobe in Fig. 1 where m 0 1000 GeV and m 1/2 2500 GeV. Once again, we see that the additional freedom of being able effectively to choose μ and M A independently allows solutions with the correct relic density over a wider range of the parameters m 0 , m 1/2 and tan β. The region of the (m 0 , tan β) plane with |m 0 | 1000 GeV is generally in the stau coannihilation region, whereas in the region at larger m 0 and tan β 40 χ ± 1 coannihilation is important. The prominent horizontal lobe in the left-hand plot at tan β ∼ 50 is associated with the A-funnel region.  . We see again that the best-fit point has m 0 < 0, and that both m 2 H u,d < 0 are favoured, with a preference for m 2 H u < m 2 H d . 8 The reason for this preference can be understood from (1). 8 However, the exact locations of the CL contours near the best-fit point in the right panel of Fig. 4 Figure 6 displays the (m 0 , A 0 ) plane (left) and the (A 0 , m 1/2 ) plane (right) for the NUHM2. The fit does not exhibit any overall preference for a sign of A 0 . However, we see that negative values of A 0 are generally preferred when m 0 and m 1/2 are large, whereas the low-mass lobe in Fig. 1 is generally associated with positive values of A 0 . 9 This tendency is driven by the value of M h measured at the LHC. Figure 7 displays the (M A , tan β) plane in the NUHM2 (solid lines), CMSSM (dashed lines) and NUHM1 (dotted lines). In the NUHM2 we see a 95 % CL lower limit on M A that increases from ∼200 GeV when tan β ∼ 5 to 1000 GeV when tan β ∼ 50, which is essentially determined by the H/A → τ τ constraint [124,125], with cut-outs due to the χ 2 penalties as different mechanisms for satisfying the χ h 2 constraint come into play or become ineffective. The bestfit value of M A ∼ 2500 GeV, but the global χ 2 function is very flat, and this model parameter is not well determined, and could be as low as 500 GeV at the 68 % CL. We find a 95 % CL lower limit tan β 4, which is quite insensitive to the value of M A . We find thatχ ± 1 coannihilation is generally important for M A 2000 GeV, whereas stau coannihilation is important for M A 2000 GeV. The A/H funnel becomes important for M A ∼ 2000 GeV, and also for tan β 50. Table 2 summarises our results for the our best-fit points in a global fit to the NUHM2, compared with fits in the NUHM1 and the CMSSM using the same post-LHC Run 1 data set. We see that the total χ 2 in the best NUHM2 fit is lowered by only χ 2 = 0.2 from the best NUHM1 fit, so the extra parameter in the NUHM2 does not provide a significant advantage. According to the F-test, there is a 77 % chance that the data are represented better by the NUHM1 than by the CMSSM, whereas there is only a 28 % chance that the NUHM2 is an improvement on the NUHM1, and a 78 % chance that the 9 We recall that we use the same convention for the sign of A 0 as in [19,92], which is opposite to the convention used in, e.g., SoftSUSY.  NUHM2 represents the data better than the CMSSM. None of these can be regarded as significant.

Summary of NUHM2 global fit
We note that the NUHM2 best-fit value of m 0 is small and negative, and that it is accompanied by values of m 2 H u and m 2 H d that are also negative and larger in magnitude. 10 We have checked the vacuum stability of the best-fit point using the Vevacious code [168], finding that it is metastable. The best-fit value of m 1/2 in the NUHM2 lies significantly beyond the direct lower limit from sparticle searches at the LHC. We also find that a positive value of A 0 is preferred, in contrast to the NUHM1 and the CMSSM which have much larger values of m 0 and m 1/2 at their best fit points. That said, we repeat that the likelihood functions are extremely shallow, and the 68 % ranges very large, so the best fit point should not be over-interpreted. 10 The best-fit value of m 2 H d lies slightly outside the nominal range of our scan, but within a region that is sampled with ∼50 % of the density within the range. We have verified using MINUIT that although the global χ 2 function is quite flat, it has a well-defined minimum at this point.

Predictions for physical observables
We now turn to the predictions for physical observables that emerge from our frequentist analysis of the NUHM2 parameter space, and compare them with corresponding predictions from our previous analyses of the CMSSM and NUHM1 parameter spaces [19]. Since the CMSSM is a subset of the NUHM1, which is itself a subset of the NUHM2, However, this is not immediately visible in the plots below, in which we plot the difference χ 2 from the minimum value of χ 2 in that model shown in the Table. In general, after falling from high values at low masses, the χ 2 are generally flat at high masses. However, there are some features associated with, for example, transitions between different mechanisms for bringing the relic density into the allowed range, which we comment on in the discussion below.

Sparticle masses
In the left panel of Fig. 8 we display the χ 2 function in the NUHM2 (solid line) as a function of mg. We see that mg 1.5 TeV is preferred at the 95 % CL, 11 as was the case in the CMSSM and NUHM1, and that the χ 2 function is quite flat for mg 2.5 TeV. The lower limit on mg is mainly due to the ATLAS jets + / E T constraint, counteracted to some extent by (g − 2) μ : the LHC M h measurement plays no role. The best-fit point has mg ∼ 3670 GeV as seen also in Fig. 2. At low masses, the χ 2 function is similar to that for the CMSSM (dotted line), and also to the NUHM1(dashed line) when mg 2 TeV. Above this mass, the difference between the χ 2 functions for the NUHM2 and the NUHM1 is largest for 3 TeV mg 5 TeV, where the extra freedom Table 2 The best-fit points found in global fits in the CMSSM, the NUHM1 and the NUHM2, using the same experimental constraints (and their theoretical interpretations): the difference in the CMSSM best-fit from that found in [19] is due to using the updated ATLAS jets + / E T constraint [10]. We note that the overall likelihood functions in all the models are quite flat, so that the precise locations of the bestfit points are not very significant, and for this reason we do not quote uncertainties. The p-values quoted would have the interpretations of probabilities if the likelihood functions given by the χ 2 statistics were Gaussian Model  The right panel of Fig. 8 displays the χ 2 likelihood as a function of mq R , defined here to be the average of the spartners of the right-handed components of the four lightest quarks. We see that mq 1.5 TeV at the 95 % CL in the NUHM2, driven essentially by the ATLAS jets + / E T constraint, with a best-fit value mq R ∼ 3080 GeV as seen also in Fig. 2, and that the χ 2 function in this model is very similar to those in the NUHM1 and CMSSM for mq R 2 TeV. However, the χ 2 functions in these models differ quite significantly for 2 TeV mq R 4.5 TeV, reflecting the fact visible in Fig. 1 that the separation between the low-and high-mass regions becomes less pronounced as the Higgs mass universality is progressively relaxed. This can be traced back to the broader range of options for bringing the cold dark-matter density into the range preferred by cosmology.
In the left panel of Fig. 9 we display the χ 2 likelihood as a function of mt 1 . In this case the lower-mass limit is not driven by the ATLAS jets + / E T search. On the other hand, the χ 2 functions for these models are quite different at both larger and smaller mt 1 : lower masses are not so strongly disfavoured in the NUHM2, and the features found in the CMSSM at mt 1 ∼ 1 TeV and ∈(2, 3) TeV are not found in the NUHM2, whose χ 2 function falls almost monotonically as mt 1 increases. This reflects again the fact that the low-and high-mass regions are less distinct in the NUHM2, whereas in the CMSSM the stau coannihilation region is quite separated from the H/A funnel region at high masses. There are also some stop coannihilation points at low mt 1 . The bestfit point has mt 1 ∼ 3420 GeV as seen also in Fig. 2. The right panel of Fig. 9 displays the χ 2 functions in the NUHM2, NUHM1 and CMSSM as functions of mτ 1 . At low mass, we see that the χ 2 functions are almost identical in the three models, giving a lower bound mτ 1 300 GeV at the 95 % CL, driven by the ATLAS jets + / E T search. At intermediate masses, the χ 2 functions in the NUHM1 and NUHM2 are reduced by the operation of extra dark matter density reduction mechanisms, which are operative in the NUHM2 also at higher masses, but not in the NUHM1 The structures seen in the χ 2 functions for the NUHM1 (dashed line) and CMSSM (dotted line) are absent for the NUHM2, whose χ 2 function (solid line) has a shallow minimum at mτ 1 ∼ 780 GeV.
The left panel of Fig. 10 displays the dependences of the χ 2 functions in the NUHM2, NUHM1 and CMSSM on M A . We see that the χ 2 function for the NUHM2 is quite flat above ∼500 GeV, following a steep rise at lower masses and a 95 % CL lower limit M A 200 GeV, which is largely due to the H/A → τ τ constraint [124,125] as mentioned previously. The best-fit point has M A ∼ 2470 GeV as seen also in Fig. 2. The right panel of Fig. 10 displays the corresponding χ 2 function for μ. Like M A , this extra degree of freedom in the NUHM2 is poorly constrained by current data. are quite similar at low masses, being largely driven by the ATLAS jets + / E T constraint, and we find that mχ0 1 250 GeV at the 95 % CL. The χ 2 function in the NUHM2 (solid line) then has a shallow minimum for mχ0 1 ∈ (600, 1000) GeV, with a best-fit value ∼760 GeV. As already mentioned, the NUHM2 best-fit point is in the stau coannihilation region, with mτ 1 − mχ0 1 ∼ 18 GeV and the other sleptons slightly heavier, as also seen in Fig. 2. In the case of mχ± 1 , the NUHM2 χ 2 function has a 95 % CL lower bound 500 GeV and a shallow minimum for mχ± 1 ∈ (1000, 1500) GeV and a best-fit value ∼1430 GeV as also seen in Fig. 2. The extra degree of freedom in the NUHM2 compared to the NUHM1 does not relax significantly the lower bounds on theχ ± 1 andχ 0 1 masses. The left panel of Fig. 12 displays the χ 2 functions for R μμ (defined here as B R(B s → μ + μ − )/B R(B s → μ + μ − ) SM ) in the NUHM2, NUHM1 and CMSSM. We see that they are almost identical, and that all three models allow no scope for R μμ to fall significantly below the SM value within the 95 % confidence level range. For R μμ above the Standard Model value, the χ 2 functions all rise in the same way as the contribution from the experimental constraint on R μμ (red line), implying that the other constraints do not impose significant constraints on R μμ above the Standard Model value. The fact that the CMSSM appears to have slightly larger freedom for R μμ is related to the fact the total χ 2 is larger than in the other models. Shifting the CMSSM curve in the right panel of Fig. 12 to account for that difference, the CMSSM region would be fully contained in the NUHM1,2 regions, as expected because of the stronger restrictions in the CMSSM.

The anomalous magnetic moment of the muon
The right panel of Fig. 12 displays the χ 2 functions for the difference from the SM: in the NUHM2, NUHM1 and CMSSM, as blue solid, dashed and dotted lines, respectively. Also shown, as a solid red line, is the (g − 2) μ contribution to the χ 2 function. As is well known, the other constraints, principally those from the LHC, do not allow a large SUSY contribution to (g − 2) μ within the NUHM1 (dashed line) or the CMSSM (dotted line). We find that in the NUHM2 the most important role is played by the LHC M h measurement. As we also see in the right panel of Fig. 12, there is significantly more flexibility in the NUHM2 contribution to (g − 2) μ (solid line). However, even in this case the model is unable to reduce the discrepancy between the theoretical prediction and the central experimental value much below the χ 2 ∼ 9 level. We find the (g − 2) μ contributions to the global χ 2 to be 9.2 (10.5) (8.8) in the CMSSM (NUHM1) (NUHM2). A reduction of the minimum value of the global χ 2 function w.r.t. the SM [19] is found at the level of χ 2 ∼ 4.0, with a best-fit value of , σ SI p ) plane in the NUHM2, with results in the CMSSM and NUHM1 shown for comparison. The star and contours have the same significations as in Fig. 1 values of σ SI p ∼ 4 × 10 −49 cm 2 , within the range where the atmospheric neutrino background dominates, are slightly favoured relative to the range σ SI p ∼ 10 −45 cm 2 . In all the three models, the steep rise in the χ 2 function at low values of σ SI p is due to the contribution from Higgs exchange via the small Higgsino component in theχ 0 1 .

Summary and conclusions
In this paper we have presented the results of a frequentist global fit of the NUHM2 model. Previous analyses of the CMSSM and NUHM1 models [19] have shown those models to be very constrained by available data. One might have wondered whether the extra degrees of freedom in the Higgs sector in the NUHM2 scenario would alleviate this tension, but we found that this was not the case. Our fit employed ∼4 × 10 8 points in the NUHM2 parameter space, and we paid particular attention to the part of the NUHM2 parameter space where m 2 0 < 0. Applying the ATLAS constraints on jets + / E T to the NUHM1,2 (and especially to m 2 0 < 0) required an extrapolation from the published results, which we previously validated for 7 TeV limits using an implementation of the Delphes collider detector simulation code set to emulate the ATLAS detector.
The minimum value of χ 2 /dof was 32.5/21, to be compared with the values χ 2 /dof ∼ 32.7/22 and 35.0/23 found in our previous analyses of the NUHM1 and CMSSM, respectively. We found that ranges of m 2 H u < m 2 H d < m 2 0 < 0 are favoured. We find similar tension between (g − 2) μ and the LHC Higgs and jets + / E T constraints in the NUHM2 as in the NUHM1 and CMSSM. The best-fit values of mg and mq R in the NUHM2 are ∼3 TeV, with χ 2 functions that are quite flat for masses 2 TeV. The freedom effectively to vary μ and M A in the NUHM2 does not suffice to provide a better fit to (g − 2) μ and suggests that if this anomaly persists then some non-universality among the SUSY-breaking scalar masses may be required.
On the one hand, it is encouraging that the results of this NUHM2 analysis are relatively similar to those found previously for the NUHM1 and the CMSSM, suggesting that the type of frequentist analysis presented here is robust with respect to simple expansions of the CMSSM parameter space. On the other hand, this analysis suggests that it would be interesting to study models in which the GUT universality assumptions are further relaxed, with a corresponding increase in the number of parameters. Such models may offer the prospect of a significant reduction in χ 2 if they can relax the tension between (g − 2) μ and the LHC constraints. Similarly, models based on a phenomenological definition of lowenergy soft supersymmetry-breaking parameters, variants of the pMSSM [170][171][172][173][174][175][176][177], may also ameliorate the tension. This may offer another path of extension beyond the well-studied CMSSM, NUHM1 and NUHM2 scenarios.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3 / License Version CC BY 4.0.