Revisiting (cid:12)ne-tuning in the MSSM

: We evaluate the amount of (cid:12)ne-tuning in constrained versions of the minimal supersymmetric standard model (MSSM), with di(cid:11)erent boundary conditions at the GUT scale. Speci(cid:12)cally we study the fully constrained version as well as the cases of non-universal Higgs and gaugino masses. We allow for the presence of additional nonholomorphic soft-terms which we show further relax the (cid:12)ne-tuning. Of particular impor-tance is the possibility of a Higgsino mass term and we discuss possible origins for such a term in UV complete models. We point out that loop corrections typically lead to a reduction in the (cid:12)ne-tuning by a factor of about two compared to the estimate at tree-level, which has been overlooked in many recent works. Taking these loop corrections into account, we discuss the impact of current limits from SUSY searches and dark matter on the (cid:12)ne-tuning. Contrary to common lore, we (cid:12)nd that the MSSM (cid:12)ne-tuning can be as small as 10 while remaining consistent with all experimental constraints. If, in addition, the dark matter abundance is fully explained by the neutralino LSP, the (cid:12)ne-tuning can still be as low as (cid:24) 20 in the presence of additional non-holomorphic soft-terms. We also discuss future prospects of these models and (cid:12)nd that the MSSM will remain natural even in the case of a non-discovery in the foreseeable future.


Introduction
The discovery of the Higgs with a mass of about 125 GeV [1,2] is, to date, the biggest success of the large hadron collider (LHC). In contrast, there has not been any evidence for new physics. This puts very strong constraints on the masses of new coloured particles as predicted, for instance, by supersymmetry (SUSY); working with very simplified assumptions, it is possible to exclude gluinos and first/second generation squarks nearly up to 2 TeV [3][4][5][6]. These experimental results raise the question of how natural are the simplest versions of supersymmetry such as the minimal supersymmetric standard model (MSSM). It has become common lore that a strong tension between the LHC results and 'natural SUSY' is now the case; for example in the Constrained MSSM it was found, because of the rather heavy Higgs mass, that one needs a fine-tuning (FT) of 300 or more to accommodate correct electroweak symmetry breaking (EWSB) [7]. This has led to increasing interest in JHEP03(2017)021 less constrained versions of the MSSM [8] and also in other SUSY models such as those with singlet extensions which can increase the Higgs mass. It was found that these models can still be considered as mildly fine-tuned even in fully constrained scenarios, see e.g. refs. [9][10][11] and references therein.
However, even in the MSSM, the situation is not as clear as often claimed; the main source of FT is usually the µ-term in the superpotential of the MSSM which gives mass to the Higgsinos and simultaneously contributes to the Higgs masses. As a result natural SUSY is usually associated with light Higgsinos, not heavier than a few hundred GeV [12][13][14][15][16][17][18][19][20][21]. It has been recently pointed out that this is not strictly correct because of the possibility of an additional source of Higgsino mass [22]; if the Higgsinos gain mass via a soft-breaking term µ , which does not affect the FT very much, it is possible to get 1 TeV Higgsinos without a large FT penalty. We show that such a large mass can arise through large radiative corrections involving Higgs portal couplings to SUSY breaking in a hidden sector. It is also possible in sequestered SUSY breaking that radiative effects generate a natural cancellation between the µ term and the soft Higgs mass terms that leaves only an Higgsino mass [23]. Moreover, if one relaxes the requirement of the unification of the gaugino mass terms at the GUT scale, regions of parameter space exist with only mild fine-tuning [24,25]. This is a natural possibility in many GUT models where the gaugino mass terms get split by the breaking of the enhanced gauge group and also in string theory models (see [26] and [25] and references therein).
The aim of this work is to perform a study of the fine tuning in constrained versions of the MSSM extended by non-holomorphic Higgsino mass and other soft-terms and to discuss the impact of gaugino mass hierarchies at the GUT scale. In the course of this work we found that the loop-corrections to the FT can be very important and typically reduce the FT prediction by a factor of 2. 1 This disagrees with the treatment of loop corrections in ref. [14] and with a series of recent papers that follow this treatment. We discuss the origin of this disagreement and give the full prescription of how to deal with loops when calculating FT. To take account of these corrections we have added a section revising the results for the MSSM. This paper is organised as follows: in section 2.1 we discuss the loop corrected FT calculation and compute the improvement compared to the previous calculations. In section 3 we give a short introduction to the MSSM extended by non-holomorphic soft-terms and we analyse possible origins of a large Higgsino mass term, both of holomorphic and non-holomorphic origin. In section 4 we present the results of exhaustive numerical scans for the FT in the MSSM, constrained by the latest experimental results. We conclude in section 5.

JHEP03(2017)021
2 Revisiting the fine-tuning calculation in the MSSM including loop corrections 2.1 Fine-tuning measure at tree-level Let us start our discussion of the MSSM fine tuning at tree-level, in order to gain some intuition and to connect to most results in the literature. After electroweak symmetry breaking (EWSB) the Higgs potential in the MSSM is given by from which one derives the minimum conditions (tadpole equations) All SUSY parameters are understood as running DR quantities at the SUSY scale. As usual, we define tan β = vu v d and note that B µ is proportional to the mass M A of the pseudo-scalar Higgs boson. In the decoupling limit (M A → ∞) and for tan β → ∞ one finds the simple relation which is often presented in the form This makes the origin of the (little) hierarchy problem within the MSSM apparent: the r.h.s. contains terms which are naturally O(M SUSY ), the SUSY breaking scale. Thus, in order to obtain the measured value of M Z there must be a cancellation between these terms which demands a certain level of tuning. There are different measures to quantify the amount of fine-tuning ∆ FT . A widely used one is the sensitivity measure proposed in refs. [28,29] ∆ ≡ max Abs ∆ p , Here, p are the independent parameters of the model, and the quantity ∆ −1 gives a measure of the accuracy to which independent parameters must be tuned to get the correct electroweak breaking scale. 2 Applying this measure to eq. (2.6), one finds The fine tuning measure can be related to a factor in a likelihood fit to the data. Interpreting this in a probabilistic sense fine tuning of O(50) adds about 0.5 to the value of χ 2 /degree of freedom in a SUSY fit to data [30].

JHEP03(2017)021
Using p 2 = {µ 2 , m H 2 u } the very naive estimate for the fine-tuning is found to be Thus, a small FT needs moderately small |µ| and |m 2 hu | at the low scale. In the following, we do not discuss the FT as function of the parameters at the SUSY scale, but consider a high-scale model in which the fundamental parameters are defined at the scale of grand unification (GUT scale). The reason is that in UV complete models it may be incorrect to treat all parameters as independent at the SUSY scale because correlations among parameters are usually present. These correlations can significantly affect the calculated fine-tuning. For instance if SUSY breaking leads to degenerate soft scalar masses there is a cancellation between the tree level and radiative contributions to the Higgs mass that leads to a reduction of the sensitivity of the Higgs mass to the initial scalar masses, the so-called 'focus point' [31][32][33][34]. As a result the fine-tuning measure is significantly reduced. It has been pointed out in refs. [25,26,[35][36][37][38] that one finds also 'gaugino focus point' scenarios if a hierarchy among the gaugino masses is assumed at the GUT scale; for example such hierarchies can result from the breaking of a larger gauge group or from an underlying string theory [25,26].
Before we turn to a detailed discussion of the FT in the MSSM assuming such a GUT model let us evaluate the effect of loop corrections on the inferred FT.

Fine tuning including loop corrections
The discussion of the impact of radiative corrections to the FT has a long history, see for instance the early works refs. [39,40], where loop corrections were found to be potentially important. More recent papers however, such as ref. [43], claim that the impact of loop corrections on the inferred FT are negligible, in particular for the FT with respect to µ. To resolve this discrepancy consider the method proposed by ref. [14] to incorporate the loop corrections by re-writing the loop corrected tadpoles as with Σ uu = ∂∆V ∂v 2 and where all loop effects are absorbed into Σ uu . This method was afterwards applied in many other works, see for instance refs. [18,19,21,. Crucially, when calculating the fine tuning, these papers all treat Σ uu as independent of v and thus do not find a correlation between the tree and loop level terms when calculating ∂ ln v 2 ∂ ln µ . As a result the FT is found to be insensitive to loop corrections. However in general Σ uu does depend on the electroweak VEV v (see below for an explicit example), implying immediately that this treatment is incorrect.
To estimate the effect of loop corrections on the FT we start with the loop corrected tadpole equation, working in the decoupling limit for simplicity (our numerical results in later sections don't rely on this simplification)

JHEP03(2017)021
where Σ u = ∂∆V ∂v . In order to take account of the correlations we use the general parametri- Using eqs. (2.11) and (2.12) one may express the loop-level fine tuning in the form which reproduces the tree level result, eq. (2.9) in the limit (Σ i = 0). The Σ i depend only very weakly on µ and for simplicity we have assumed ∂Σ i ∂µ = 0 as in ref. [43]. Thus, the dominant change in the FT measure does not come from the variation in Σ, but is due to the overall size of Σ 2,3 .

Example: stop corrections without mixing
Let us discuss the simple example of stop corrections without mixing and with degenerate stop masses to illustrate our points. The loop corrections to the scalar potential can be written as t is a universal stop soft-breaking mass, Y t is the DR top Yukawa coupling and Q is the renormalisation scale. For this, one finds Clearly Σ uu does depend on the electroweak vev v. Using our parametrisation from above, the coefficients in eq. (2.12) can be expressed as For this simple example, the change in the fine-tuning from tree-to loop-level can be approximated as follows In the following we assume that the Σi are independent of v. This is not strictly speaking correct because of the appearance of v in the logs in the loop corrections. This effect is neglected here in the analytical discussion but taken into account in our numerical checks.

JHEP03(2017)021
where we set Q 2 = m 2 t and used mt m t . In figure 1 we plot r FT as function of the SUSY scale M SUSY showing that one can expect that the FT improves by about a factor of two if the loop corrections are correctly included. We also show the FT with respect to µ as a function of µ as well as the simple analytical estimate for the FT. Given that the reduction in FT is about a factor of 2 (see also [41,42] for a similar result), we propose to use when estimating the FT with respect to µ.

The new soft-terms
In our discussion of the estimate of the FT we have seen that, despite the reduction in fine tuning due to loop effects, small FT prefers a moderately small µ-term. Therefore the Higgsino mass mH , which is mainly given by µ in the MSSM, is strongly limited by the requirement of small FT. However, this constraint can be softened if there are other sources of Higgsino mass. This leads us to consider the possibility of additional soft SUSY breaking terms including a non-holomorphic contribution to the Higgsino mass, In the context of the MSSM these terms are still 'soft' in the sense that they do not lead to quadratic divergences at radiative order [83], such quadratic divergences only appearing in models with additional gauge singlet fields. Usually the non-holomorphic Higgsino

JHEP03(2017)021
mass term is simply not written down, as it can be reabsorbed into the other soft-terms and superpotential parameters. As shown in [22] however, this term can be crucial when considering FT and should therefore not be absorbed. 4

Impact and origin of large µ
Of particular relevance to the Higgsino mass bound is the possibility of the Higgsino mass µ in eq. (3.1). While the soft-terms in eq. (3.1) keep the minimum conditions in eqs. (2.3)-(2.4) unchanged, the µ term shift the Higgsino mass to It has been pointed out in ref. [22] that a large µ can have very important consequences on the FT. This is in particular the case when the constraints from the dark matter relic are included. The correct abundance of Higgsino dark matter is naturally obtained for Higgsino masses of around 1 TeV. However, this would immediately correspond to a FT of at least 120 as can be seen from eq. (2.20). On the other hand, if the Higgsino mass stems to a large extent from µ , 1 TeV Higgsinos could still be consistent with a FT well below 50 [22]. We are going to extend the analyses of ref. [22] in section 4 by performing exhaustive parameter scans for a variety of different GUT boundary conditions. First, however, we discuss schemes in which a large µ term can arise in a UV complete model.
As discussed by Martin [91], in superspace co-ordinates the µ term has the form In the MSSM the Higgsino gets mass via its coupling to the Majorana masses of the gauginos. For gaugino masses of O(M SUSY ), much greater than the EW breaking scale this mass is of O(v 2 /M SUSY ) and so is much smaller than the SUSY breaking scale. However this suggests that large Higgsino masses may be generated by coupling to heavy gauginos in the SUSY breaking sector. As a concrete example consider a U(1) hidden sector with gaugino mass generated by the term 1 M d 2 θXW α W α when F X = 0. In addition to the MSSM Higgs fields H d,u we add Higgs portal mediator fields Φ 3,4 , Φ, Φ and a further field Y that also acquires an F-term. In table 1 we exhibit the transformation properties of the fields under the Standard Model gauge group SU(2) L × U(1) Y and the additional hidden sector U(1) gauge group. With the charge assignments shown the superpotential has the form where Y obtains a scalar VEV Y = M as well as an F-term VEV F Y . With these couplings the graph of figure 2 generates an Higgsino mass mh uhd ∝ where the JHEP03(2017)021  Figure 2. Graph generating Higgsino mass via gaugino mediation. messenger masses in the denominator depend on the details of the model. One may readily check that there is no equivalent one-loop scalar mass term for the MSSM Higgs fields. 5 5 Another way to obtain an Higgsino mass without the SUSY related Higgs scalar mass term has been pointed out by Perez, Roy and Schmaltz [23]. It results from renormalisation group (RG) scaling from approximate strong dynamics in a hidden sector. The authors show that, starting with a supersymmetric µ term for the Higgs and an independent soft Higgs scalar mass, the RG running drives the Higgs scale soft mass squared to −µ 2 , naturally cancelling the supersymmetric contribution to the Higgs mass squared term where the equality of coefficients is a natural result of the RG running. However, although this gives a Higgsino mass without a Higgs scalar mass, this does not come from the operator in eq. (3.3). As a result, even in models with Standard Model singlet superfields such as the NMSSM, there are no quadratic divergences and so the mechanism can be applied to these models too. An analysis of this case is beyond the scope of this paper. A Higgsino mass without its scalar counterpart can also be generated via the Scherk-Schwarz mechanism (see [86]). More recently constructions of natural SUSY models relying on the Scherk-Schwarz mechanism have received some interest [87][88][89]. In this case the mass is comparable to the scale of compactification and the new states associated with the extra dimension spoil the simple gauge coupling unification considered in this paper (see however [90] for an attempt to achieve gauge coupling unification also in these kinds of set-ups).

Numerical setup
For our analysis we have used the Mathematica package SARAH [92][93][94][95][96][97] in order to generate a spectrum generator based on SPheno [98,99]. This provides a full-fledged spectrum generator including the RGEs for the standard SUSY terms up to two loops, based on refs. [100], while the non-holomorphic soft-terms are fully included at one-loop level based on the results of refs. [101,102]. The SUSY mass spectrum is calculated at full one-loop, while for the Higgs states the important two-loop corrections are covered in the effective potential approach [103,104]: this includes also the one-and two-loop corrections to the Higgs mass via the non-holomorphic soft-terms. All important flavour constraints such as b → sγ, ∆M Bq , B → ll are checked via the FlavorKit functionality [105]. The routine to calculate the fine-tuning via SARAH/SPheno has been improved during this work to include also the loop corrections as discussed in section 2.2. 6 The numerical procedure works now as follows: for each independent parameter (e.g. m 0 , m 1/2 , A 0 , µ, B µ ) the following steps are performed for a given parameter point This iteration also ensures that the VEV-dependence of all logs in the loop function is fully included. The overall fine-tuning is taken to be the maximal value of the fine-tuning calculated in this way, ∆ = max ∆ p

Boundary conditions
In the following we numerically study the amount of FT in constrained versions of the MSSM with and without non-holomorphic Higgsino masses. We consider both the fully constrained versions where both the gaugino masses and all scalar masses unify at the GUT scale as well as the case of non-universal gaugino and Higgs masses. In all cases we assume universal squark and slepton masses The changes are public with version 4.9.2 of SARAH. and we parametrise the trilinear soft-terms as

JHEP03(2017)021
As usual we take the trilinear A terms to be proportional to the corresponding superpotential couplings and for simplicity we assume the same relation for the A 0 term. Our results are rather insensitive to this choice. The case of non-universal gaugino masses we parametrise via coefficients a and b at the GUT scale: In total we study 6 different boundary conditions, which are defined in table 2. All parameters discussed above as well as the value of the Higgsino mass, µ , are defined at the GUT scale. In addition, we fix the values of µ and B µ by requiring correct EW symmetry breaking (EWSB) and use tan β as well as the phase of µ as input. The fine-tuning is calculated with respect to (m 0 , m 1/2 , A 0 , A 0 , µ, µ , B µ ) and for the case of non-universal Higgs masses also with respect to (m 2 hu , m 2 h d ) We do not include any FT penalty for a and b assuming that they are fixed in the underlying theory. 7

Constraints
A dedicated re-casting of the collider limits from SUSY searches is beyond the scope of this work. Therefore we require a simplistic cut on the gluino mass depending on the neutralino mass based on ref. [3,5]. Of course, in realistic SUSY models, the limits on the different masses can be significantly weaker compared to the results of simplified models [106,107]. In addition we require all charged electroweak states such as charginos to have masses larger than 100 GeV. In addition, we require that the Higgs mass should be in the range 122 GeV < m h < 128 GeV. Since the pseudoscalar and heavy Higgs are usually much heavier for all points which survive the cuts on the SUSY masses, the Higgs coupling rates are very SM-like and no additional constraints come from their measurements. 7 As discussed in e.g. [25] there are a number of high-scale scenarios which predict fixed ratios among the gaugino masses, in which case an infinitesimal variation of a and b is not meaningful.

JHEP03(2017)021
In a second step we assume that the lightest neutralino makes up (part of) the dark matter and we assume standard thermal freeze out as the main production mechanism. To have a consistent scenario we require that either the neutralino dark matter does not over-close the universe, i.e. Ωh 2 < 0.13 (weak requirement), or that it explains the entire DM abundance in the universe, 8 0.11 < Ωh 2 < 0.13 (strong requirement). We also require that the neutralino passes the latest direct detection limits from LUX [109]. To evaluate the LUX constraint we rescale the bound by Ωh 2 in the case of an under-abundance.
When discussing the dark matter properties of a given scenario, one should not forget that there is (at least) a gravitino in the particle spectrum beyond the usual SM superpartners, with a mass depending on the underlying supersymmetry breaking mechanism. If the gravitino is the LSP, it is a good dark matter candidate and the lightest neutralino becomes the NLSP and will decay into the gravitino. In this case the bounds from the relic abundance don't apply. If the gravitino is heavier, it can be an additional source of dark matter production via late decays, depending on the cosmological evolution of the Universe. With this in mind it is obvious that the possible cuts due to the dark matter properties are on a different footing than the limits from collider searches and hence we show the results of these cuts individually below.

The fully constrained MSSM
Let us start our discussion of the fine-tuning in the fully constrained MSSM with (CNHSSM) and without (CMSSM) the new soft terms. Although FT in the CMSSM has been widely discussed, we redo the analysis here to take account of the loop corrections discussed above that have often been missed. Specifically we perform a random scan in the following parameter space: 0 < m 0 < 3000 GeV 100 < m 1/2 < 3000 GeV 10 < tan β < 50 −3000 GeV < A 0 < 3000 GeV sign µ ± 1 For the CNHSSM we in addition scan over The main results of these scans are shown in figures 3 and 4. In figure 3 we show the fine tuning against the SM-like Higgs mass and against the lightest neutralino mass for the CMSSM and the CNHSSM. Here and in the following we distinguish four cases: JHEP03(2017)021 . Light gray points are before any cuts. Medium gray points are after bounds from collider searches have been taken into account. The light coloured points in addition pass the weak dark matter cuts while the dark coloured points include the lower bound on Ωh 2 and hence require the correct relic abundance. It can be seen that the lightest neutralino in the CMSSM is always a bino, while in the CNHSSM also a Higgsino LSP is possible (see text for colour coding). Bottom: FT vs. the mass of the lightest neutralino.
• Points with a well defined spectrum without any further cuts (light gray) • Points passing the cuts on the Higgs mass and LHC limits (medium gray) • Points passing the weak dark matter cuts (medium colour) • Points passing the strong dark matter cuts (dark colour) The coloured points contain additional information about the composition of the lightest neutralino, with red for a bino, green for a Wino and blue for a Higgsino. Mixtures between these different states are shown via RGB colour coding.
As expected from the discussion in section 2.1, compared to a tree-level analysis, the FT becomes smaller by about a factor of 2 if the loops are included properly. If only the collider bounds are applied, the smallest FT we find within the CMSSM is about 130 JHEP03(2017)021 (cf. table 3 for a summary of our findings for all considered cases) and quickly rises with the Higgs and neutralino masses. For a unified gaugino mass m 1/2 at the high scale it is well known that at the low scale the ratios of gaugino mass parameters roughly scale as M 3 : M 2 : M 1 ∼ 6 : 2 : 1. Hence the bino is always the lightest gaugino for a universal m 1/2 . In addition Higgsinos turn out to be typically significantly heavier, such that the lightest neutralino is an almost pure bino, which can also be inferred from the colour coding in figure 3. Due to its rather small couplings a bino LSP has a small annihilation cross section and is therefore typically overproduced implying that only a small number of points pass all the cuts. This is because one needs to sit in fine-tuned co-annihilation or resonance regions, which is a potential tuning beyond the one in the electroweak sector.
For the case of the CNHSSM, the picture after collider constraints have been taken into account is rather similar as far as the FT is concerned, with the smallest FT of about 110. The large difference however concerns the composition of the lightest neutralino: due to the extra contribution to the Higgsino mass light Higgsinos are well possible, and it is much easier to obtain the correct relic abundance, as can be seen from figure 3. However, only JHEP03(2017)021 the usual µ-term enters the electroweak symmetry breaking conditions and hence there is no difference in this respect compared to the CMSSM: the usual µ-term is still large and while the extra term µ can make the Higgsinos much lighter, the FT is still sizeable.
Future prospects. In the coming years the LHC will continue to collect large amounts of data and simultaneously dark matter direct detection experiments will improve the bounds on the dark matter nucleon scattering cross section by about 2 orders of magnitude. Given these expected improvements it is interesting to ask how naturalness will be constrained in a given scenario in the foreseeable future. The answer can be inferred from figure 4. In the top panel we show the FT in the m χ -mg plane to estimate the sensitivity of future collider searches. In the bottom panel we show the FT in the m χ -σ p plane, where σ p is the DM nucleon scattering cross section. Here we rescale the DM nucleon scattering cross section for the case of an LSP underabundance in order to be able to sensibly compare to direct detection bounds, σ p → σ p · Ω th DM /Ω obs DM . For both cases we show the current bounds in the given plane, with the other collider bounds taken into account in both cases and in addition the dark matter cuts applied to the bottom panels. In the bottom panels we also show the expected sensitivity of the planned XENONnT or LZ [110] experiments. We see that for the CMSSM the valid points in the m χ -mg plane are on a line with m χ ∼ 1/6mg as expected for a bino LSP due to the gaugino mass relations. For the CNHSSM the LSP can be lighter than 1/6mg if there is a sizeable Higgsino component, but of course not heavier (a heavier Higgsino would no longer be the LSP). In both cases the FT significantly increases with increasing gluino mass, with a minimal FT for mg ≥ 3 TeV (≥ 5 TeV) of about 230 (700) in both cases (cf . table 3). Constraining naturalness even more may be a more precise (theoretical) determination of the Higgs mass, given the steep slope in FT. Direct detection experiments are typically not sensitive with respect to the FT within the CMSSM, due to the small couplings of the bino LSP, while for the CNHSSM the remaining parts of parameter space with FT about 100 will be tested by future direct detection experiments.

Non-universal Higgs masses
Non-universal Higgs masses are often discussed as a first step beyond the CMSSM. One of the motivations is that in many GUT models different masses for the Higgs fields are possible, as they usually originate from a different GUT multiplet compared to the matter fields. We therefore briefly discuss the results for the MSSM with non-universal Higgs mass terms with and without the new soft terms (denoted NHSSM-NUHM and MSSM-NUHM respectively). In order to minimise bias in our results we solve the tadpole equations not only with respect to µ and Bµ as before but also with respect to the soft Higgs masses m 2 h d and m 2 hu and overlay the different scans. Similarly as before our results are summarised in figures 5 and 6. It turns out that after collider cuts have been taken into account the resulting FT is even slightly worse than in the fully constrained versions. This is likely due to the fact that the approximate focus point behaviour of a universal m 0 is no longer present. On the other hand the LSP is typically a mixture of bino and Higgsino and a number of points have the correct relic abundance. Regarding future prospects this JHEP03(2017)021 scenario is again rather similar to the previous case, with FT increasing significantly with gluino mass.

Non-universal gaugino masses
Let us finally come to the case of non-universal gaugino masses, i.e. we relax the assumption that all gaugino masses unify at the GUT scale and allow for the following ratios in the gaugino masses at the high scale: Non-universal gaugino masses have attracted quite some attention due to the existence of a gaugino focus point, allowing for a significant reduction in FT. This is illustrated in figure 7 where we show the FT as function of the gaugino mass ratios a, b. Note that these ratios can still be consistent with an underlying GUT structure, see e.g. [26] for more details.
Similarly as before the main results are summarised in figures 8 and 10. We see that for both cases, the MSSM-NUGM as well as the NHSSM-NUGM, a large reduction in FT JHEP03(2017)021  is possible. Also the increase of FT with m h is significantly milder than in the previous cases. After the collider constraints are taken into account the FT in both cases can be as small as 10. From figure 8 also the flexibility within the gaugino sector is evident, allowing for bino (red), Higgsino (blue) and Wino (green) LSPs. Higgsinos as well as Winos annihilate rather efficiently in the early Universe, leading to an underabundance in the DM relic density if their masses are below a TeV. Therefore requiring the weak dark JHEP03(2017)021 matter cut with no lower bound on the relic abundance, the FT stays as small as 10 in both cases (requiring however an alternative DM source such as axions or axinos or a different cosmological history). When applying the strong DM cut to explain the dark matter via the neutralino LSP, the Higgsino (Wino) mass gets pushed to about 1 TeV (2 TeV) as can be seen in figure 8. 9 In the MSSM this requires a µ term of about 1 TeV, resulting in significant tuning (of about 120). In the NHSSM large Higgsino masses can be obtained due to the extra µ contribution without a large FT penalty, allowing for FT of about 20 with the correct relic abundance. For clarity we show black lines in figure 8 to illustrate this. The full black line corresponds to the expectation ∆ tree = 2 from the MSSM where the Higgsino mass is up to small corrections given by the µ-term. The dotted black line is a rough estimate for the loop corrections using ∆ loop = 1 2 ∆ tree . One can see that this approximation works well in the usual MSSM while in the presence of an extra contribution to the Higgsino mass the FT can be much smaller, since the Higgsino mass and µ term JHEP03(2017)021 are now independent. To illustrate this further, we show in figure 9 the minimal FT in the (µ, µ )-plane after the collider cuts (left) and after the additional requirement of the correct relic abundance (right). It can be seen that before the lower bound on the relic abundance is imposed there is no preference for µ = 0, which is however strongly preferred once this bound is implemented.

Future prospects
Let us finally discuss how future collider and direct detection experiments test the natural range of parameters for the case of non-universal gaugino masses. In figure 10 we again show the FT in the in the m χ -mg plane to estimate the sensitivity of future collider searches (top panel) as well as the FT in the m χ -σ p plane (bottom panel). Interestingly, once we impose the mass of the SM-like Higgs boson to be in the measured range, the smallest fine tuning is found for gluino masses above the current bounds and the region with small FT extends to very large gluino masses. In fact for gluino masses as heavy as 5 TeV the fine tuning can still be between 20 and 30, as can also be inferred from table 3. This observation applies both to the case with and without the extra non-holomorphic contribution to the Higgsino mass, although the case with extra terms is slightly preferred. What is different however is the dependence on the LSP mass, as expected: without the extra contribution there is a large FT penalty when going to large LSP masses, while with this contribution the increase in fine tuning is very mild. Turning to the dark matter properties of these models, it can be seen that the region with small FT extends down to very small direct detection cross section. Indeed the cross sections are often below the expected sensitivity of LZ and XENONnT and can be as small as the neutrino background so that it seems impossible to probe the entire natural parameter space of the model.
To summarise we present in table 3 the minimum FT ∆ min for all considered cases after applying the various cuts discussed in the text, both for the current bounds and estimating the sensitivity of future experiments.

Summary and conclusions
In this work we have reassessed the FT in the MSSM for a number of different GUT scale SUSY breaking boundary conditions. Specifically we considered the fully constrained MSSM with and without an additional non-holomorphic Higgsino mass term. We further relaxed the condition of universality of Higgs and gaugino masses at the GUT scale, see table 2 for all considered boundary conditions. We find that due to the proper inclusion of loop corrections the FT is typically smaller by about a factor of two compared to a tree-level analysis of the sensitivity measure for the well-known case of the CMSSM. The inclusion of the non-holomorphic Higgsino mass does not appreciably change the FT for this case, but allows for more flexibility in the neutralino sector, making it possible to satisfy the dark matter relic abundance constraint much more easily. Non-universal Higgs masses don't improve on the FT -in fact the FT is typically even worse due to the loss of the scalar focus point. The situation is very different for non-universal gaugino masses, which allow for a significant reduction in FT, assuming that the ratio of gaugino masses is fixed within the underlying theory. In this case we find points with FT as small as 10 which pass all the collider and dark matter limits, if one allows for a dark matter underabundance. If we require the correct dark matter density, which typically requires a Higgsino LSP with mass of about 1 TeV, the FT increases significantly, becoming larger than 100 due to the large µ term when no additional contribution to the Higgsino mass is taken into account. With JHEP03 (2017) Table 3. The best fine-tuning ∆ min found after applying the different cuts in the MSSM for all considered boundary conditions. The first three columns correspond to the current limits with cuts as described in the text. The last four columns are sensitivity estimates for a lower bound on the gluino mass mg of 3 and 5 TeV respectively. The columns denoted by DD in addition take into account a bound on the LSP nucleon cross section as expected for the future direct detection experiments XENONnT and LZ.
such an additional contribution however the FT can stay as small as ∼ 20 even for a 1 TeV Higgsino.
We also evaluated the prospects of testing the naturalness of these different models within the next 10-20 years, taking into account more stringent cuts on the gluino mass as well as future direct detection experiments. We find that while the fully universal case as well as the case with non-universal Higgs masses will be significantly more constrained, the case of non-universal gauginos allows for very large gluino masses without a large penalty on the FT measure. In fact when taking into account the known SM-like Higgs mass, the region with smallest FT is for gluino masses above current bounds. Even taking into account future direct detection experiments such as XENONnT and LZ, these low fine JHEP03(2017)021 tuned regions survive, offering an interesting way to evade all bounds while still maintaining one of the original motivations of SUSY, i.e. fully solving the hierarchy problem of the SM.