LPT Orsay 1163
L2C 11
LUPM 11002
Naturalness and Fine Tuning in the NMSSM: Implications of Early LHC Results
Ulrich Ellwanger, Grégory EspitalierNoël and Cyril Hugonie
LPT, UMR 8627, CNRS, Université de Paris–Sud, 91405 Orsay, France
L2C, UMR 5221, CNRS, Université de Montpellier II, 34095 Montpellier, France
LUPM, UMR 5299, CNRS, Université de Montpellier II, 34095 Montpellier, France
We study the fine tuning in the parameter space of the semiconstrained NMSSM, where most soft Susy breaking parameters are universal at the GUT scale. We discuss the dependence of the fine tuning on the soft Susy breaking parameters and , and on the Higgs masses in NMSSM specific scenarios involving large singletdoublet Higgs mixing or dominant HiggstoHiggs decays. Whereas these latter scenarios allow a priori for considerably less fine tuning than the constrained MSSM, the early LHC results rule out a large part of the parameter space of the semiconstrained NMSSM corresponding to low values of the fine tuning.
1 Introduction
The first motivation for supersymmetric extensions of the Standard Model (SM) stems from the solution of the naturalness or fine tuning problem in the Higgs sector of the SM [1, 2, 3, 5, 4] (besides the unification of the gauge couplings and the possibility to explain dark matter): In the SM with an ultraviolet cutoff much larger than the electroweak scale , the bare Higgs mass squared must satisfy roughly . Hence must be of the order , but must be finetuned relative to with a precision of the order . This fine tuning is enormous for of the order of a GUT scale, but only a tuning of is considered as natural. Within supersymmetric (Susy) extensions of the SM with Susy breaking terms of the order , the necessary tuning between the parameters is of the order of , and hence independent from an ultraviolet cutoff .
First results of searches for Susy by the ATLAS and CMS collaborations at the LHC, based on fb of data taken at 7 TeV centerofmass energy, imply lower bounds on Susy breaking gluino and up/down squark masses in the 1 TeV range [6, 7] (for the latest publications, see the ATLAS and CMS notes on the web pages [8, 9]). These bounds reduce the phenomenologically viable range of parameters in Susy extensions of the SM. However, an obviously interesting question is the impact of these negative results on the necessary tuning between the parameters in the remaining parameter region.
It is well known that the nonobservation of a Higgs boson at LEP [10] implies already a “little fine tuning problem” in the Minimal Supersymmetric Standard Model (MSSM) (see e.g. [11, 12, 13, 14]), where the field content in the Higgs sector is as small as possible, but large radiative corrections are required in order to lift the mass of the lightest neutral CPeven Higgs boson above the lower LEP bound. Large radiative corrections require relatively large Susy breaking top squark masses compared to the electroweak scale. Via loop diagrams, large top squark masses lead to relatively large soft Susy breaking Higgs mass terms, which require some tuning among the parameters of the MSSM such that the Higgs vacuum expectation values (vevs) are of , well below the scale of the Higgs mass terms. The required tuning among the parameters of the MSSM is typically of the order of a few %.
The “little fine tuning problem” of the MSSM, originating from LEP constraints, is alleviated in certain regions of the parameter space of the NexttoMinimal Supersymmetric Standard Model (NMSSM). The NMSSM is the simplest Susy extension of the SM with a scale invariant superpotential, i.e. where the only dimensionful parameters are the soft Susy breaking terms. No supersymmetric Higgs mass term is required as in the MSSM, since it is generated dynamically by the vacuum expectation value of a gauge singlet superfield (see [15, 16] for recent reviews). Together with the neutral components of the two SU(2) doublet Higgs fields and of the MSSM, one finds three neutral CPeven and two CPodd Higgs states in this model.
The additional coupling of to and can lead to a larger mass of the SMlike neutral Higgs boson , and to mixings of the physical CPeven Higgs bosons in terms of the weak eigenstates , and , implying reduced couplings of the physical eigenstates to the boson. Both phenomena make it easier to satisfy the LEP bounds [17], and allow to alleviate the little finetuning problem [18]. (See [19] for an evaluation of the upper bound on of about 140 GeV if is required to remain perturbative below the GUT scale; if is allowed to be larger, can be even larger [20].)
Moreover, can decay preferably into lighter NMSSMspecific singletlike Higgs bosons. In this case LEP bounds on are lower, and again the required fine tuning can be considerably smaller than in the MSSM [21, 22, 23, 18, 24]. Consequently it becomes important to study the impact of the recent bounds from the LHC on the fine tuning within the NMSSM, which is the purpose of this paper.
The most frequently used quantitative measure for fine tuning is the maximum of the logarithmic derivative of with respect to all fundamental parameters (if the fundamental Lagrangian is given at the GUT scale) [25, 26, 27, 28, 11, 29, 30, 31, 32, 33]:
(1) 
(See [34, 35] for alternatives; sometimes instead of is used in the argument of the logarithm, which leads to an obvious additional factor of 2. Subsequently we prefer to study linear relations between all masses and couplings.) depends on the point in parameter space and is, roughly speaking, inversely proportional to the required fine tuning (as discussed above) between the parameters . Hence, for a given point in parameter space, should be as small as possible, preferably of . Preferred regions in the parameter space spanned by are those where is minimal (denoted by ). In practice, the value of depends on the choice of independent fundamental parameters defining the model, and on the implementation of phenomenological constraints as the dark matter relic density, the anomalous magnetic moment of the muon etc..
Including WMAP constraints on the dark matter relic density (but leaving aside the top Yukawa coupling in the list of ), has been studied recently within the constrained MSSM (cMSSM, with universal soft Susy breaking terms at the GUT scale) in [30, 31, 32, 33]. First investigations of the impact of the early LHC results on the Susy parameter space in the cMSSM have been performed in [36, 37, 38, 39], in the cMSSM with nonuniversal sfermion masses in [40], and within the general MSSM in [41, 38, 39].
Compared to alternative procedures as likelihood scans and/or Bayesian techniques (see [42, 43] for studies within constrained versions of the NMSSM, and [44] for a recent discussion) a disadvantage of is that it does not allow to marginalise (i.e. to integrate over) parts of the parameters which, in turn, would allow to determine “most likely” values for given quantities as masses of specific particles, given all present experimental constraints. (Clearly, such predictions for “most likely” masses seem to be of limited use; for instance, they would have failed miserably if applied to the SMlike Higgs mass in the preLEP era.)
One can also leave aside the issue of quantitative fine tuning, content oneself with the fact that the fine tuning in Susy is always much smaller than in the SM with a large cutoff , and determine “most likely” values for parameters exclusively from best fits to data from electroweak precision experiments. The impact of recent LHC bounds on such best fits has been studied recently in [45, 46, 48, 47, 49, 43].
On the other hand, the constraints from recent or future LHC results on the quantitative fine tuning measure can contribute to the discussion on the impact of LHC results on Supersymmetry in general. Hence we will compare these constraints within the parameter space of the semiconstrained sNMSSM (where singletspecific soft terms are allowed to be nonuniversal) to the cMSSM, obtained from the sNMSSM in the limit . Therefore we study the dependence of on the universal soft Susy breaking parameters and (gaugino and scalar masses, respectively) and on the gluino and up/down squark masses. In particular we investigate the relevance of NMSSM specific scenarios in the Higgs sector for fine tuning, as scenarios with large singlet/doublet mixing and scenarios with dominant decays [21, 22, 23, 24, 16, 43].
In the next Section we define the model, the procedure for the determination of , and discuss some specific properties of in the NMSSM. Our results and their discussion are given in Section 3, and we conclude with a summary in Section 4.
2 Fine Tuning in the NMSSM and the MSSM
The NMSSM differs from the MSSM by the presence of the gauge singlet superfield . The Higgs mass term in the superpotential of the MSSM is replaced by the coupling of and to and a selfcoupling , hence the superpotential is scale invariant in this simplest invariant version of the NMSSM:
(2) 
where we have confined ourselves to the Yukawa couplings of and to the quarks and leptons , , and of the third generation; a sum over the three generations is implicitly assumed. (In (2), for the first and the last time, the fields denote superfields.) Once assumes a vev , the first term in generates an effective term
(3) 
The soft Susy–breaking terms consist of mass terms for the gaugino, Higgs and sfermion fields
(4)  
as well as trilinear interactions between the sfermion and the Higgs fields, including the singlet field
(5)  
(Again, an effective MSSMlike parameter is generated.)
All parameters in the above Lagrangian depend on the energy scale via the corresponding renormalization group (RG) equations, which account for the dominant radiative corrections involving large logarithms. In the constrained NMSSM, one imposes unification of the soft Susy–breaking gaugino masses, sfermion and Higgs masses as well as trilinear couplings at the grand unification (GUT) scale :
(6) 
Since the singlet superfield could play a special role (its couplings to a hidden sector, responsible for supersymmetry breaking, could differ from the MSSMlike fields), we will allow for nonuniversal singletspecific soft terms and at the grand unification scale. This is the socalled semiconstrained NMSSM, denoted by sNMSSM subsequently. Including the top quark Yukawa coupling due to its influence on the RG equations, the Lagrangian of the sNMSSM depends on eight parameters at the GUT scale:
(7) 
The calculation of defined in (1) proceeds in two steps: First, we compute the variations
(8) 
with respect to the parameters at the Susy scale (the Susy scale is defined to be of the order of the soft Susy breaking terms). Subsequently these variations are contracted with the Jacobian
(9) 
which takes care of the renormalization group running of the parameters between the Susy and the GUT scales. Then we obtain
(10) 
The parameters at the Susy scale are those appearing in the (effective) Higgs potential, whose minimization determines the vevs of , of and of . These vevs determine, in turn, the quantities
(11) 
Including the dominant top quark/squark induced radiative corrections, the three minimisation equations are given by [16]
(12) 
where denotes an average value of the top squark masses. (It is not necessary to be more precise here, in contrast to the radiative corrections to the physical Higgs masses.) It is straightforward to express the vevs , and in terms of , and with the help of these equations.
Hence the relevant parameters at the Susy scale are given by (leaving aside the electroweak gauge couplings and , as well as inside the logarithm)
(13) 
In order to compute the variations (see (8)) with respect to these parameters, we use
(14) 
for . Since all partial derivatives of the equations can be computed explicitely, the three equations (14) can be solved for (and, separately, for and ) as function of all , which allows to determine the variations in (8).
At this stage it is useful to recall the origin of the “little fine tuning problem” in the MSSM. Neglecting the radiative corrections, the minimisation equations (12) of the Higgs potential imply, with in the MSSM,
(15) 
In the absence of fine tuning, all terms on the right hand side of (15) should be of comparable magnitude, and no large cancellations should occur; hence both and should not be much larger than . However, from the RG equations one typically obtains , which is often required to be much larger (in absolute value) than : At least within the MSSM, the SMlike Higgs scalar mass increases proportionally to due to top/stop induced radiative corrections. Then, large values for are unavoidable in order to satisfy the LEP bound. Albeit large stop masses are consistent with the nonobservation of top squarks, they would generate an uncomfortably large value for which has to be cancelled by in (15).
For large one finds for , following (8) with or ,
(16) 
Accordingly large values for (leading, generally, to large values for ) reflect well the necessary fine tuning if and hence are large.
In the NMSSM is replaced by . For large , the above reasoning remains essentially unchanged: For (valid in most of the parameter space), in (12) gives
(17) 
Replacing and (17) for in (15), one finds again from (8)
(18) 
in the NMSSM. Hence, quite obviously, large values for are unnatural as well. However, due to the NMSSM specific contributions to the Higgs masses and mixings or NMSSM specific Higgs decays, LEP bounds on the Higgs sector can be satisfied for smaller top/stop induced radiative corrections, hence for smaller values of , allowing for smaller values for and .
It remains to express the variations of the parameters at the Susy scale in terms of variations of the parameters at the GUT scale, i.e. to compute the Jacobian in (9), via the integration of the RG equations for the parameters. In the cMSSM (with boundary conditions at the GUT scale as in (6)) one can always write
(19) 
and
(20) 
where the coefficients and depend on the gauge and Yukawa couplings.
In the typical case where all in (19) satisfy , one can verify that all variations , with , are less than 1. At first sight, due to for these parameters , this seems to reduce the necessary fine tuning in the MSSM. However, due to , remains always large. Moreover, should generally be included in the list of parameters [50, 27], and the corresponding variation can be large. This holds particularly in the socalled focus point region of the MSSM where , and in (19) for specific values of (but a large derivative of with respect to ).
In the sNMSSM, additional terms and appear on the right hand side of (19), which have little impact in practice. Instead of (20), the RG equations for , and will now play some role since, replacing (17) for in , depends on these parameters (apart from a dependence on ). In fact one can verify that, in the MSSM limit of the NMSSM where , , satisfies the same RG equation as . All in all we cannot expect dramatic effects on the fine tuning from the somewhat different running of the parameters between the Susy and the GUT scale in the NMSSM.
In practice our computation of the different variations in the space of parameters (7), in order to find its maximum as function of (see (1)) at a specific point in the parameter space, is performed as follows: For each such point, the code NMSPEC [51] inside NMSSMTOOLs [52, 53] is used in order to compute the Higgs and sparticle spectrum including radiative corrections as described in these references. Constraints from LEP, Bphysics and the anomalous magnetic moment of the muon are taken care of according to the latest updates given on the web site http://www.th.upsud.fr/NMHDECAY/nmssmtools.html. (No constraints on the dark matter relic density are imposed; however, in most cases these could be satisfied by giving up the bino mass unification at the GUT scale, i.e. chosing an appropriate mass for the lightest neutralinolike Susy particle (LSP) without impact on the results relevant here.)
For phenomenologically acceptable points, is computed from the three minimization equations in (12), following the procedure described above. The Jacobian , i.e. the variations of the parameters at the Susy scale in terms of variations of the parameters at the GUT scale, is computed numerically from the two loop RG equations. This allows to obtain the necessary quantities , whose maximum with respect to all parameters (7) defines .
3 Results for the cMSSM and the sNMSSM
To start with, we apply our procedure to the cMSSM, which allows for comparisons with the sNMSSM and the available literature. The relevant parameters at the GUT scale are
(21) 
As usual, and (and hence and ) are determined by and , but they still contribute to the definition of . Next we scan over a grid of values of and . For each set of these values, we scan over and . Keeping only sets of parameters consistent with all phenomenological constraints, we look for values of and which minimize as defined in (1) for fixed and . The resulting minimal values of can be represented in the plane , or in the plane (where refers to squarks of the first generation) in Figs. 1.
In order to guide the eye, we have indicated lower bounds from ATLAS and CMS notes on analyses of jets and missing , based on an integrated luminosity fb: In the plane , lower bounds from ATLAS [6] (interpreted within the cMSSM with , ) are shown as a black line, and lower bounds from CMS [7] are shown as a red line. In the plane , lower bounds from ATLAS [6] are shown as as a black line. (The latter bounds from ATLAS are obtained in a simplified model where squarks decay only into quarks + a neutralino with a branching ratio of 100 %.)
In the white regions in Figs. 1, phenomenological constraints cannot be satisfied for any values of and : Either a stau would be the LSP (left hand side of the left panel), or a charged slepton, chargino, neutralino or a CPeven Higgs boson is excluded by LEP2/Tevatron (lower part of the left panel), or squarks are excluded by CDF/D0 (left hand side of the right panel) or are theoretically unaccessible (right hand side of the right panel).
Note that, for GeV and GeV in the left panel, decreases hardly with decreasing Susy breaking parameters and ; the minimal value of in the preLHC allowed region is about . In this region of the parameter space we observe the ”little fine tuning problem” of the MSSM due to the LEP bound on the SMlike Higgs mass, which hardly depends on and . The impact of the LEP bound becomes clear once we minimize the fine tuning for fixed lightest Higgs mass (without imposing LEP constraints on the Higgs sector) as a function not only of and , but also of and . The result is shown in Fig. 2.
We see the strong increase of with for GeV. Accordingly the LEP constraint GeV implies in agreement with Figs. 1 (for low values of and ). In fact, is always just above 114 GeV in the entire planes in Fig. 1, once and are chosen such that is minimized for fixed and , but LEP constraints are applied.
Hence, for low values of and (in the region GeV and GeV), is determined by the LEP constraints on (implying the little fine tuning problem of the MSSM) and not much affected by the lower bounds on gluino and squark masses from early LHC searches: In the region above the ATLAS/CMS bounds, the minimal value of increases only to .
We also see in Fig. 2 that does not decrease systematically with , once lower bounds on sparticle masses from LEP and the Tevatron are imposed: In order to minimize , preferred values of would be in the range GeV which is excluded in the MSSM, but not in the NMSSM (see below).
For larger values of or , the origin of the required fine tuning is different: Here it is simply the fact that the weak scale (determined essentially by ) is small compared to the Susy breaking scale, which requires some tuning between the parameters. Since at the Susy scale is closely related to the squark masses, increases rapidly with (for TeV) as it is visible in the right panel in Figs. 1.
The fine tuning in the cMSSM has recently been analysed in [30, 31, 32, 33]. The procedure and precision in these papers is similar to ours, except that constraints on the dark matter relic density are applied in [30, 31, 32, 33], but contributions to from the top Yukawa coupling are left aside. From Fig. 7d in [31] we find, once LEP constraints are applied, a minimal value of for not too large values of . Given that in [31] is twice as large as our defined in (1), this coincides well with the left panel of Figs. 1 for moderate values of . However, for GeV, decreases to in [31], whereas increases with in the left panel of Figs. 1. In fact, for larger values of , is dominated by contributions from whose absence in [30, 31, 32, 33] explains the different results for in this region.
Next we turn to the sNMSSM. In various regions of the parameter space of the sNMSSM, unconventional properties of the Higgs sector allow to alleviate the LEP bounds, lowering the minimal possible values of . We found it interesting to study these lower bounds on separately for different scenarios in the NMSSM Higgs sector (see also [18]), since these will have very different implications for future Higgs searches at the LHC. Hence we distinguish subsequently the following two scenarios:
(1) The lightest CP even Higgs boson has a large singlet component (H/S mixing). This implies a reduced coupling to the boson, and allows for masses well below 114 GeV [10].
(2) A CPeven Higgs boson decays dominantly into a pair of lighter CPodd bosons, , allowing again for masses well below 114 GeV [10]. (We omit the index of in the following.)
The search for the minimal fine tuning in each of these scenarios is performed similar to the procedure in the cMSSM: Again we scan over a grid of values of and . Now, for each set of these values, we scan over , , and ( and are determined by and , but included in the definition of ) using a Monte Carlo Markov Chain (MCMC) technique. Keeping only sets of parameters consistent with all phenomenological constraints, we look for values of , , and which minimize as defined in (1) for fixed and , allowing us to represent the resulting minimal values of in the plane , or in the plane . (In order to distinguish the scenarios above, we require essentially for scenario (1), but for scenario (2).)
In the scenario (1) (H/S mixing), the corresponding results for are shown in Figs. 3. Now the constraint on the left hand side in the left panel from the absence of a stau LSP has disappeared, since a singlinolike neutralino can be the LSP. From here onwards, the bounds from ATLAS and CMS are indicative only, since the signals for supersymmetry in the NMSSM can be different notably in the case of a singlinolike LSP [54].
Compared to the fine tuning in the cMSSM in Figs. 1, we see that can be considerably smaller for not too large values of the Susy breaking parameters and [18]: A large singlet component of (in the range ) allows for lighter masses compatible with LEP constraints, which reduces the required fine tuning, see Fig. 4 below. The mass of the second mostly SMlike CPeven Higgs boson is always just above 114 GeV. This is now easier to satisfy than in the MSSM, since the doublet/singlet mixing shifts the mass of the mostly doubletlike Higgs boson upwards. The values of the NMSSMspecific coupling do not have to be large to this end; its value is always .
Respecting just the preLHC phenomenological constraints, the fine tuning measure can be as small as 14 for low values of and in this scenario, but this is precisely the region which is constrained by the first unsuccessful searches for Susy at the LHC. (As stated above, these constraints depend on the decay properties of the u/dsquarks and gluinos. Additional bino neutralino decay processes can reduce signature in all sparticle decay cascades. Applying nevertheless the bounds for the cMSSM scenarios studied by the ATLAS and CMS collaborations to the sNMSSM, we find that the smallest admissible value of becomes for GeV, GeV, similar to the smallest admissible value of in the cMSSM.)
For larger values of and , the fine tuning within this scenario (1) of the sNMSSM becomes similar to the one within the cMSSM: As explained above, the origin of the fine tuning is now the smallness of the weak scale with respect to and not the LEP constraints on the Higgs mass; hence the possibility to alleviate the LEP constraints within the sNMSSM is less relevant.
In Fig. 4 we show as function of the mass of the dominantly singletlike Higgs state, minimizing as a function of , , , taking into account LEP constraints from Higgs and sparticle searches. (The irregular structures originate from the LEP constraints on , which lead to irregular upper bounds on the coupling of to the boson as function of .) Again, does not decrease systematically with decreasing given the phenomenological constraints on sparticle masses, but the behaviour is somewhat different from the dependence of on the SMlike Higgs mass in the MSSM. Here we find that is minimal for below GeV which coincides with the Higgs mass range where a 2.3 excess in is observed at LEP 2 [10]. (This excess could be explained here since, although is dominantly a gauge singlet, it still has a nonvanishing – but reduced – coupling to the boson, and will decay dominantly into .)
Next we turn to the sNMSSM scenario (2) involving light pseudoscalars , into which the SMlike CPeven Higgs boson can decay. Again the required fine tuning is reduced [21, 22], since LEP constraints allow for smaller masses than 114 GeV. These constraints depend on :
a) For GeV, will dominantly decay into a pair of quarks. This signature has been studied by the OPAL and DELPHI groups at LEP [55, 56] implying GeV if has SMlike couplings to the boson [10]. Still, this lower bound on allows for lower values of .
b) For GeV, decays dominantly into a pair of leptons. The signature has recently been reanalysed by the ALEPH group [57] for GeV, implying again lower limits on . However, for in the range GeV, can and will mix strongly with the CPodd states [58] which implies a considerable reduction of the [59]. Hence bounds on hardly constrain in this case; now we find that the dominant constraints result from the lower limits on depending on remaining subdominant [10].
Our results for the minimal fine tuning in the sNMSSM scenario (2) are shown in Figs. 5 in the same planes as before. Similar to the sNMSSM scenario (1), respecting just the preLHC constraints, the fine tuning measure can be as small as 9 for low values of and . Applying naively the bounds for the cMSSM scenarios studied by the ATLAS and CMS collaborations to the sNMSSM, we find that the smallest admissible value of becomes for GeV, GeV, hence below the smalles admissible value in the cMSSM.
In the region of low , the is larger than 80 %. Again, the values of the NMSSMspecific coupling do not have to be large in order to favour decays once these are kinematically allowed; the value of varies between 0.02 and 0.16 in Figs. 5.
In Fig. 6 we show the minimal value of as function of the mass of the now SMlike state , imposing LEP constraints. Again the preferred value of is not always as small as possible; now the minimal fine tuning is obtained for in the range GeV, and increases again strongly for larger values of .
It is also interesting to study the dependence of the fine tuning on . At first sight, depends hardly on within the sNMSSM, where differs from at the GUT scale: A variation of from to GeV (still allowing for decays) corresponds to a variation of from GeV to GeV at the GUT scale, which has practically no effect on the fine tuning measure . However, different values for correspond to different lower LEP bounds on ; lower bounds , in turn, affect as in Fig. 6. The resulting minimal values of as function of are shown in Fig. 7.
The structures in Fig. 7 can be explained as follows: For GeV, LEP bounds on (here with , and ) apply, see Figs. 3 in [10]. The corresponding lower bounds on are somewhat weaker for GeV than for GeV, leading to a lower fine tuning for GeV where GeV. (For GeV, the lower bound on has a dip implying a dip in . For GeV, increases since GeV is required for kinematic reasons; here this scenario is obviously not preferred.) For GeV, strong lower bounds on from ALEPH [57] on the apply at first sight. However, due to a reduced branching ratio for form mixing [59], the window GeV is hardly affected by these constraints, allowing for to assume values for which is minimal according to Fig. 6. (For GeV, parameters have to be tuned in order to satisfy constraints from .) Hence we find two distinct regions for where decays allow for a considerable reduction of the fine tuning.
4 Conclusions
We have studied the amount of fine tuning in the parameter space of the semiconstrained NMSSM, and compared it to the cMSSM. Representing the minimal fine tuning in the plane allowed us to study the impact of early LHC results. First we verified quantitatively, to what extend the NMSSMspecific scenarios in the Higgs sector allow to alleviate the LEP constraints with respect to the MSSM. We found indeed that a considerable reduction of the fine tuning is possible in scenarios where the lightest CPeven Higgs state is dominantly singlet like, or in scenarios where decays are possible, notably for GeV and .
If one applies naively the early LHC constraints on the cMSSM in the plane to the sNMSSM, the impact on the minimal fine tuning is stronger in the sNMSSM since these constraints affect in particular the region of smaller Susy breaking terms where the fine tuning in the sNMSSM can be relatively small. Still, the necessary fine tuning in the sNMSSM in scenarios where decays are possible is smaller than in the cMSSM, but only if the Susy breaking terms are not too large (hence if sparticles are not too heavy). Otherwise, if the Susy breaking terms are larger, the fine tuning does not originate from LEP constraints on the Higgs sector, but from the smallness of the weak scale with respect to the Susy scale. This hardly depends on details of the Higgs sector of the Susy model, as long as the fundamental parameters are defined at the GUT scale and influence each other through the renormalization group running between the GUT and the weak scale. Hence the same question should be reanalysed in the NMSSM with gauge mediated supersymmetry breaking [60], where the Susy breaking parameters can originate from much lower scales, and where the Susy breaking Higgs mass terms can differ considerably from squark mass terms.
Acknowledgements
U. E. acknowledges support from the French ANR LFVCPVLHC. C. H. acknowledges support from the French ANR TAPDMS (ANR09JCJC0146).
References
 E. Witten, Nucl. Phys. B 188 (1981) 513.
 S. Dimopoulos and H. Georgi, Nucl. Phys. B 193 (1981) 150.
 E. Witten, Phys. Lett. B 105 (1981) 267.
 N. Sakai, Z. Phys. C 11 (1981) 153.
 R. K. Kaul and P. Majumdar, Nucl. Phys. B 199 (1982) 36.
 ATLAS Collaboration, ATLCOMPHYS2011981, presentation at the 2011 Europhysics Conference On High Energy Physics, Grenoble, France, 21  27 July 2011
 CMS Collaboration, “Search for supersymmetry in allhadronic events with ”, CMSPASSUS11003, presentation at the 2011 Europhysics Conference On High Energy Physics, Grenoble, France, 21  27 July 2011
 https://twiki.cern.ch/twiki/bin/view/AtlasPublic/SupersymmetryPublicResults
 https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsSUS
 S. Schael et al. [ALEPH and DELPHI and L3 and OPAL Collaborations and LEP Working Group for Higgs Boson Searches], Eur. Phys. J. C 47 (2006) 547 [arXiv:hepex/0602042].
 P. H. Chankowski, J. R. Ellis and S. Pokorski, Phys. Lett. B 423 (1998) 327 [arXiv:hepph/9712234].
 R. Barbieri and A. Strumia, Phys. Lett. B 433 (1998) 63 [arXiv:hepph/9801353].
 G. L. Kane and S. F. King, Phys. Lett. B 451 (1999) 113 [arXiv:hepph/9810374].
 L. Giusti, A. Romanino and A. Strumia, Nucl. Phys. B 550 (1999) 3 [arXiv:hepph/9811386].
 M. Maniatis, Int. J. Mod. Phys. A 25 (2010) 3505 [arXiv:0906.0777 [hepph]].

U. Ellwanger, C. Hugonie and A. M. Teixeira,
Phys. Rept. 496 (2010) 1
[arXiv:0910.1785 [hepph]].  M. BasteroGil, C. Hugonie, S. F. King, D. P. Roy and S. Vempati, Phys. Lett. B 489 (2000) 359 [arXiv:hepph/0006198].
 R. Dermisek and J. F. Gunion, Phys. Rev. D 77 (2008) 015013 [arXiv:0709.2269 [hepph]].
 U. Ellwanger and C. Hugonie, Mod. Phys. Lett. A 22 (2007) 1581 [arXiv:hepph/0612133].
 R. Barbieri, L. J. Hall, Y. Nomura and V. S. Rychkov, Phys. Rev. D 75 (2007) 035007 [arXiv:hepph/0607332].
 R. Dermisek and J. F. Gunion, Phys. Rev. Lett. 95 (2005) 041801 [arXiv:hepph/0502105].
 R. Dermisek and J. F. Gunion, Phys. Rev. D 73 (2006) 111701 [arXiv:hepph/0510322].
 R. Dermisek and J. F. Gunion, Phys. Rev. D 76 (2007) 095006 [arXiv:0705.4387 [hepph]].
 R. Dermisek, Mod. Phys. Lett. A 24 (2009) 1631 [arXiv:0907.0297 [hepph]].
 J. R. Ellis, K. Enqvist, D. V. Nanopoulos and F. Zwirner, Mod. Phys. Lett. A 1 (1986) 57.
 R. Barbieri and G. F. Giudice, Nucl. Phys. B 306 (1988) 63.
 B. de Carlos and J. A. Casas, Phys. Lett. B 309 (1993) 320 [arXiv:hepph/9303291].
 P. Ciafaloni and A. Strumia, Nucl. Phys. B 494 (1997) 41 [arXiv:hepph/9611204].
 P. H. Chankowski, J. R. Ellis, M. Olechowski and S. Pokorski, Nucl. Phys. B 544 (1999) 39 [arXiv:hepph/9808275].
 S. Cassel, D. M. Ghilencea and G. G. Ross, Phys. Lett. B 687 (2010) 214 [arXiv:0911.1134 [hepph]].
 S. Cassel, D. M. Ghilencea and G. G. Ross, Nucl. Phys. B 835 (2010) 110 [arXiv:1001.3884 [hepph]].
 S. Cassel, D. M. Ghilencea, S. Kraml, A. Lessa and G. G. Ross, JHEP 1105 (2011) 120 [arXiv:1101.4664].
 S. Cassel and D. M. Ghilencea, “A review of naturalness and dark matter prediction for the Higgs mass in MSSM and beyond,” arXiv:1103.4793 [hepph].
 G. W. Anderson and D. J. Castano, Phys. Lett. B 347 (1995) 300 [arXiv:hepph/9409419].
 K. L. Chan, U. Chattopadhyay and P. Nath, Phys. Rev. D 58 (1998) 096004 [arXiv:hepph/9710473].
 A. Strumia, JHEP 1104 (2011) 073 [arXiv:1101.2195].
 S. Akula, N. Chen, D. Feldman, M. Liu, Z. Liu, P. Nath and G. Peim, Phys. Lett. B 699 (2011) 377 [arXiv:1103.1197].
 J. A. Conley, J. S. Gainer, J. L. Hewett, M. P. Le and T. G. Rizzo, “Supersymmetry Without Prejudice at the 7 TeV LHC,” arXiv:1103.1697 [hepph].
 M. Farina, M. Kadastik, D. Pappadopulo, J. Pata, M. Raidal and A. Strumia, “Implications of XENON100 results for Dark Matter models and for the LHC,” arXiv:1104.3572 [hepph].
 K. Sakurai and K. Takayama, “Constraint from recent ATLAS results to nonuniversal sfermion mass models and naturalness,” arXiv:1106.3794 [hepph].
 S. Scopel, S. Choi, N. Fornengo and A. Bottino, Phys. Rev. D 83 (2011) 095016 [arXiv:1102.4033].
 D. E. LopezFogliani, L. Roszkowski, R. R. de Austri and T. A. Varley, Phys. Rev. D 80 (2009) 095013 [arXiv:0906.4911].
 J. F. Gunion, D. E. LopezFogliani, L. Roszkowski, R. R. de Austri and T. A. Varley, “NexttoMinimal Supersymmetric Model Higgs Scenarios for Partially Universal GUT Scale Boundary Conditions,” arXiv:1105.1195.
 F. Feroz, K. Cranmer, M. Hobson, R. Ruiz de Austri and R. Trotta, JHEP 1106 (2011) 042 [arXiv:1101.3296].
 B. C. Allanach, Phys. Rev. D 83 (2011) 095019 [arXiv:1102.3149].
 P. Bechtle et al., “What if the LHC does not find supersymmetry in the sqrt(s)=7 TeV run?,” arXiv:1102.4693 [hepph].
 B. C. Allanach, T. J. Khoo, C. G. Lester and S. L. Williams, JHEP 1106 (2011) 035 [arXiv:1103.0969].
 O. Buchmueller et al., Eur. Phys. J. C 71 (2011) 1634 [arXiv:1102.4585].
 P. Bechtle et al., “Present and possible future implications for mSUGRA of the nondiscovery of SUSY at the LHC,” arXiv:1105.5398 [hepph].
 G. G. Ross and R. G. Roberts, Nucl. Phys. B 377 (1992) 571.
 U. Ellwanger and C. Hugonie, Comput. Phys. Commun. 177 (2007) 399 [arXiv:hepph/0612134].
 U. Ellwanger, J. F. Gunion and C. Hugonie, JHEP 0502 (2005) 066 [arXiv:hepph/0406215].
 U. Ellwanger and C. Hugonie, Comput. Phys. Commun. 175 (2006) 290 [arXiv:hepph/0508022].
 U. Ellwanger, A. Florent and D. Zerwas, JHEP 1101 (2011) 103 [arXiv:1011.0931 [hepph]].
 G. Abbiendi et al. [OPAL Collaboration], Eur. Phys. J. C 37 (2004) 49 [arXiv:hepex/0406057].
 J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 38 (2004) 1 [arXiv:hepex/0410017].
 S. Schael et al. [ALEPH Collaboration], JHEP 1005 (2010) 049 [arXiv:1003.0705].
 F. Domingo, U. Ellwanger and M. A. SanchisLozano, Phys. Rev. Lett. 103 (2009) 111802 [arXiv:0907.0348 [hepph]].
 F. Domingo and U. Ellwanger, JHEP 1106 (2011) 067 [arXiv:1105.1722].
 U. Ellwanger, C. C. JeanLouis and A. M. Teixeira, JHEP 0805 (2008) 044 [arXiv:0803.2962].