NMSSM with generalized deflected mirage mediation

We propose to generate a realistic soft SUSY breaking spectrum for Next-to-Minimal Supersymmetric Standard Model (NMSSM) with a generalized deflected mirage mediation scenario, in which additional Yukawa and gauge mediation contributions are included to deflect the renormalization group equation trajectory. Based on the Wilsonian effective action obtained by integrating out the messengers, the NMSSM soft SUSY breaking spectrum can be given analytically at the messenger scale. We find that additional contributions to mS2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_S^2$$\end{document} can possibly ameliorate the stringent constraints from the electroweak symmetry breaking and 125 GeV Higgs mass. Constraints from dark matter and fine-tuning are also discussed. The Barbieri–Giudice fine-tuning measure and electroweak fine-tuning measure in our scenario can be as low as O(1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal O}(1)$$\end{document}, which possibly indicates that our scenario is natural.


Introduction
Low energy supersymmetry (SUSY), which can solve elegantly the gauge hierarchy problem of the Standard Model (SM) and provide a viable dark matter (DM) candidate, has been regarded by many physicists as one of the most appealing candidates for the TeV-scale new physics. However, the reported data of the Large Hadron Collider (LHC) agree quite well with the SM predictions and no significant deviations have been observed in the electroweak precision measurements and the flavor physics. Besides, the lack of SUSY signals at the LHC [1][2][3] and the difficulty to accommodate the discovered 125 GeV Higgs [4,5] seem to indicate that the low energy SUSY spectrum should display an intricated structure. As the low energy soft SUSY breaking spectrum is determined by the SUSY breaking mechanism, it is interesting to survey the phenomenology related to the SUSY breaking and mediation mechanism from a top-down approach.
SUSY breaking mechanism from flux compactification of type IIB string theory can lead to interesting consequences. In the generalized Kahler-modulus dominated scenarios, the dilaton and complex structure moduli fields can be stabilized by the background NS and RR 3-form fluxes. Such superheavy modes will be integrated out and eliminated from the low energy effective theory. The remaining Kahler moduli fields can be stabilized by non-perturbative effects, such as the gaugino condensation. To generate SUSY breaking in the observable sector and tune the cosmological constant to a tiny positive value, Kachru et al. propose to add an anti-D3 brane at the tip of the Klebanov-Strassler throat to explicitly break SUSY and lift the AdS universe to obtain a dS one [6]. Consequently, the F-component of the light Kahler moduli field could mediate the SUSY breaking effects to the visible sector and result in a mixed modulus-anomaly mediation SUSY breaking scenario [7,8]. It is interesting to note that the involved modulus mediated SUSY breaking contributions can be comparable to that of the anomaly mediation [9,10]. With certain assumptions on the Yukawa couplings and the modular weights, the change of the SUSY breaking contributions (from the inputs) by the renormalization group equation (RGE) evolution and the anomaly mediation contributions could cancel each other at a 'mirage' scale, leading to an apparent pure modulus mediation scenario at such a mirage scale [11]. Such a mixed modulus-anomaly mediation SUSY breaking mechanism is dubbed as 'mirage mediation'.
Anomaly mediation contribution is a crucial ingredient of such a mixed modulus-anomaly mediation mechanism. It is well known that purely anomaly mediation mechanism is bothered by the tachyonic slepton problem [12][13][14][15][16][17][18]. One of its non-trivial extensions with a messenger sector, namely the deflected anomaly mediation SUSY breaking (AMSB) scenario, can elegantly solve such a tachyonic slepton problem through the deflection of the RGE trajectory [19][20][21][22][23][24][25] by additional gauge [26][27][28][29][30][31][32] or Yukawa mediation contributions. Such a messenger sector can also be present in the mirage mediation scenarios and play an important role in generating a preferable low energy SUSY spectrum. In fact, the 'mirage' unification of gaugino masses at TeV scale can be possible in deflected mirage mediation scenarios, even with the simplest minimal KKLT set-up [33][34][35]. In the mixed modulus-anomaly mediation mechanism, a straightforward extension of μ sector to predict a small Bμ term needs some fine tuning [11]. Such a μ − Bμ problem can be solved naturally in the framework of Next-to-Minimal Supersymmetric Standard Model (NMSSM). The realization of NMSSM in TeV mirage mediation scenarios had already been discussed in the literature [36][37][38][39]. However, it was found that only a small portion of the parameter space can be consistent with the electroweak symmetry breaking (EWSB) conditions and at the same time accommodate the 125 GeV Higgs [36,37]. So it is rather interesting to seek new ways to generate a realistic NMSSM spectrum.
Additional gauge or Yukawa mediation contributions in the mirage mediation [33][34][35] can deflect the RGE trajectory and change the low energy soft SUSY predictions. Analytical expressions for the soft SUSY breaking parameters at the messenger scale are given in the Wilsonian effective action approach [40] by one of the authors. Based on the general discussions in [40], we propose to generate the NMSSM spectrum by a generalized deflected mirage mediation mechanism [41] with additional gauge and Yukawa mediation contributions. We find that the inclusion of the messenger sector and non-trivial interactions can possibly alleviate the stringent constraints from the 125 GeV Higgs and the EWSB conditions.
This paper is organized as follows. We briefly review the mirage mediation mechanism in Sect. 2. We state our model and present the analytical expressions for the soft SUSY parameters in Sect. 3. The numerical results are discussed in Sect. 4. Section 5 contains our conclusions.

Brief review of the mirage mediation mechanism
In Type IIB string theory compactified on a Calabi-Yau orientifold, the presence of background fluxes can fix the dilaton and the complex structure moduli, leaving only the Kahler moduli in the four-dimensional Wilsonian effective supergravity action (defined at the boundary scale ) after integrating out the superheavy complex structure moduli and dilaton. The low energy effective Lagrangian in terms of compensator field and a single Kahler modulus that parameterizes the overall size of the compact space [11] is given as with a holomorphic gauge kinetic term The Kahler potential involves the 'no-scale' kinetic form for the Kahler modulus while the superpotential involves the KKLT setup where the first term is generated from the fluxes and the second term from non-perturbative effects, such as gaugino condensation from the non-abelian gauge sector or D3-brane instantons. W 0 denotes the superpotential terms involving the MSSM (NMSSM) sector as well as a possible messenger sector. The modulus T , which is not fixed by the background flux, can be stabilized by such a KKLT-type superpotential with . Note that in the KKLT setup, the flux-induced SUSY breaking is dynamically canceled by the non-perturbative SUSY breaking that stabilizes the Kahler modulus T , leading to a SUSY-preserving AdS solution. In order to obtain a vacuum with a positive cosmological constant and break SUSY, KKLT proposed to add anD3 brane to provide an uplifting operator given by with a positive constant D ∼ O(m 2 3/2 M 2 pl ). The uplifting operator, which represents the low energy consequence of the sequestered SUSY-breaking brane, is independent of visible matter fields and T (with n P = 0) in the minimal KKLT set-up. Explicit SUSY breaking via anti-D3 branes can be replaced by typical D-term or F-term uplifting mechanisms [42][43][44].
With the uplifting low energy effective potential, we have the leading order F-terms of compensator field F φ and Kahler modulus F T [11] The non-zero F-term VEV of the heavy moduli H are given approximately by F H ∼ m 2 3/2 /m U m 3/2 . Therefore, it gives negligible contributions to SUSY breaking [45].
In the mirage mediation, we have m 3/2 ≈ (4π 2 )M 0 numerically, which means that the loop induced anomaly mediation contributions are comparable to the modulus mediation contributions. The importance of the anomaly mediation contributions relative to the modulus mediation contributions can be parameterized by So α 1 corresponds to the limit of pure modulus mediation, while α 1 corresponds to the pure anomaly mediation. Although the minimal KKLT predicts α ≈ 1, other values of α ∼ O(1) can be obtained by assuming proper uplifting operator and proper forms for the potential of Kahler modulus [11]. So we leave the value α as a free parameter in the following discussions.

NMSSM with deflected mirage mediation
The general form of the low energy effective action Eq. (2.1) can be amended to include new ingredients, for example, the NMSSM sector and a messenger sector. The Kahler potential will include not only the 'no-scale' kinetic form for the Kahler modulus but also additional kinetic terms for messenger fields withP i , P i the messenger superfields and the NMSSM superfield. The Kahler metric for matter fields and messengers as well as holomorhic gauge kinetic functions are assumed to depend non-trivially on the Kahler modulus T The values of n X , n , n P , l i depend on the location of the fields on the D3/D7 branes. Besides, universal l i = 1 are adopted in our scenario so that the gauge fields reside on the D7 brane and the gauge coupling unification can be preserved.
The W 0 term within the superpotential Eq. (2.3) is given as with the messenger sector. The '2m' family of messengers can be fitted into 5,5 representation of SU (5) gauge group The VEV of the pseudo-modulus superfield X , which can be determined by W (X ) and other SUSY breaking effects, can set the messenger thresholds through the X P mPm type couplings. The deflection parameter 'd', which characterizes the relative size of contributions between the gauge (Yukawa) mediation and anomaly mediation, can be readily obtained Similar to the deflected AMSB, a positive deflection parameter in the deflected mirage mediation, which can be generated by a carefully chosen superpotential etc. [21][22][23][24][25], may be preferable because less messenger species are needed so that the problem of strong gauge couplings below the GUT scale can be evaded. The purpose to introduce even number of messenger species is to forbid possible kinetic mixing between X and S, otherwise the tadpole term for S would destabilize the weak scale [46]. The discrete Z 3 breaking by EWSB may generate domain walls in the early universe which may lead to an unacceptably large anisotropy of CMB. To avoid such a problem, the Z 3 symmetry is assumed to be broken by some higher dimensional operators. There are two approaches to obtain the low energy SUSY spectrum in the deflected mirage mediation: • In the first approach, the deflected mirage mediation soft SUSY spectrum is given by the expressions in [11] at the boundary scale. Such a spectrum will receive additional contributions towards its RGE running to low energy scale, especially the threshold corrections related to the appearance of messengers [47]. • In the second approach which we will adopt, the soft SUSY spectrum at low energy scale is derived directly from the (low energy) effective action. The SUGRA description in Eq. (2.1) can be seen as the Wilsonian effective action after integrating out the superheavy complex structure moduli and dilaton field. After the pseudomodulus X acquires a VEV and determines the messenger threshold, the messengers can be integrated out to obtain a low energy effective action below the messenger threshold scale. So we anticipate that the Kahler metric Z and gauge kinetic f i should depend non-trivially on the messenger threshold M 2 mess /φ † φ and M mess /φ, respectively. The resulting soft SUSY spectrum below the messenger threshold can be derived from the wavefunction renormalization approach [48][49][50][51][52][53][54]. The analytical expressions for the soft SUSY breaking parameters in the most general form of deflected mirage mediation scenarios are given in [40] by one of the authors.

The modular weight choices for NMSSM superfields
We need to specify the modular weights n i ≡ 1 − m i for NMSSM superfields before we can carry out numerical analysis. In the NMSSM, successful EWSB and the solution to the μ-problem in general require a large VEV for the singlet S. So it is preferable to introduce a negative m 2 S and/or large trilinear terms A λ , A κ for the singlet superpotential interactions. As a negative m 2 S prefers vanishing modulus contributions, we set the modular weight m S = 0. The choice of modulus weight for H u , H d can be understood from the EWSB conditions in NMSSM. From we can see that m 2 H u should be light to avoid a too large finetuning. On the other hand, m 2 The electroweak (EW) naturalness prefers relatively light stops. In the MSSM, light stops below 1 TeV are disfavored because it is difficult to accommodate the observed 125 GeV Higgs. However, large loop corrections involving stops are not necessarily required to interpret the 125 GeV Higgs in NMSSM. So light stops, which is preferable from the criterion of a low EW fine-tuning, are still allowed in NMSSM. Squarks of the first two generations should be heavy to avoid various SUSY CP and flavor constraints. We note that even pure AMSB contributions can guarantee the heaviness of colored SUSY particles. Besides, it is preferable to ameliorate the g μ − 2 discrepancy in the framework of SUSY with light sleptons and electroweakinoes. With the modular weight l i = 1 for gauge fields, a positive deflection parameter 'd' can possibly guarantee the lightness of the electroweakinos.
The notorious tachyonic slepton problem in AMSB can be solved in our scenario. Positive slepton masses can be obtained by introducing a proper deflection parameter 'd' or by adding extra modulus mediation contributions. So we chose the following m i in our scenario: • Modular weights for sleptons are given by • Modular weights for other matter and messenger fields are given by (3.9) • Double messenger species with m = 1 are adopted in our subsequent numerical study.
Note that the messenger modular weights also play a role and contribute to m 2 S . The modular weights n i = 0 correspond to matter fields on D7 branes while n i = 1 on D3 branes. Modular weights n i = 1/2 corresponds to fields on the intersections of the D3-D7 branes.

Analytical expressions for soft SUSY breaking parameters
Now we use the second Wilsonian approach in our analysis. From the analytical expressions for the generalized mirage mediation [40], the soft SUSY breaking parameters in NMSSM at the messenger scale M mess after integrating out the messenger fields can be given explicitly. The gaugino masses are given as with l i = 1 and Within the expression, the relative size between the anomaly and modulus mediation contribution at the messenger scale is determined by the free parameter We define the modular weight q y i jk as (3.14) The trilinear soft terms at the messenger scale are given by (3.19) with new contributions due to non-vanishing G S

20)
The soft SUSY breaking parameters for the scalars can be parameterized by with each part given as follows • The pure modulus contribution part with b i = b i + b i being the beta function upon the messenger thresholds and Expressions for the first two generations can be obtained by setting y t = y b = 0. • The deflected anomaly mediation part with 42) and N = 2 in our scenario for For the third generationQ L ,3 ,Ũ c L , we need to include the 'y t , y b ' Yukawa contributions The contributions to δ IS can be divided into three parts We have the pure anomaly mediation part Besides, the m 2 S term receives new contributions involving λ P because G S is not continuous across the messenger threshold So the Yukawa mediation contribution is The anomaly-gauge (Yukawa) mixing term is given by (3.56) • The interference terms involving the Kahler modulus 'T ': Note that the expressions for sfermions are hold for the third generation, the first two generation can be obtained by setting y t , y b → 0. Within the expressions, modular weight l i = 1 for gauge couplings are used. The previous contributions are just the anomaly-modulus interference part. Possible modulus-gauge interference part will also appear in our scenario. The anomalous dimensions for all fields except S are continuous across the messenger threshold, so their modulus-gauge interference contributions vanish. As noted previously, the anomalous dimension for S is discontinuous across the messenger threshold, so we have the new T, X interference contributions to m 2 S T X (3.66) Here

Numerical results
After fixing the modular weights, the remaining free parameters in our scenario are with F φ /(16π 2 ) ≈ α M 0 and the simplest choice λ D P = λ T P = λ D X = λ T X = λ 0 in our numerical study. Note that for later convenience, the definition of α is four times smaller than α that appears in Eq. (2.8). The ratio α between F φ /(16π 2 ) and M 0 holds in the messenger scale and in general is different from its value at the GUT scale. We choose a positive α in our numerical study. For a negative α, virtual mirage unification at a super-GUT energy scale will appear.
In order to transform m 2 H u , m 2 H d , m 2 S into μ, tan β, M 2 Z , we use the following approximation to obtain μ and tan β iteratively. We use NMSSMTools5.2.0 [55,56] to scan the whole parameter space. The parameters are chosen to satisfy: 10 15 GeV > M mess > 10 5 GeV, In our scan, we impose the following constraints: 122 GeV < M h < 128 GeV from ATLAS and CMS data [4,5]. Note that the uncertainty is 3 GeV instead of the default 2 GeV because a large λ may induce an additional 1 GeV correction to m h at two-loop level [61], which is not included in the NMSSMTools. (III) The relic density of the neutralino dark matter should satisfy the Planck data DM = 0.1199 ± 0.0027 [62] in combination with the WMAP data [63] (with a 10% theoretical uncertainty). (IV) The electroweak precision observables [64] and the lower bounds on neutralinos and charginos, including the invisible decay bounds for Z -boson. The most stringent constraints of LEP require mχ± > 103.5GeV and the invisible decay width (Z →χ 0χ0 ) < 1.71 MeV, which is consistent with the 2σ precision EW measurement non−SM inv < 2.0 MeV [65]. (V) Flavor constraints [66] from B-meson rare decays The deviation is about 3σ a μ (expt − SM) = (255 ± 80) × 10 −11 . (4.12) We adopt a conservative estimation 4.7 × 10 −10 a μ 52.7 × 10 −10 in our numerical results.
We should note that the numerical results depend crucially on whether the 125 GeV Higgs is the lightest CP-even scalar (Type A) or the second lightest CP-even scalar (Type B). We have the following discussions • The low energy soft SUSY breaking spectrum of NMSSM, determined from a top-down approach by a UV-completed theory, is always bothered by the requirement to achieve successful EWSB. As noted previously, EWSB conditions in NMSSM in general require a large VEV for the singlet and consequently prefer a negative m 2 S and/or large A λ , A κ for the singlet potential. However, ordinary mirage mediation scenarios always predict large positive values for m 2 S and not very large A λ , A κ , suppressing the singlet VEV. In our scenario, because of possible negative contributions to m 2 S from new Yukawa interactions, stringent constraints from successful EWSB can be ameliorated. The observed 125 GeV Higgs mass, whether it is the lightest or the second-lightest CP-even scalar in NMSSM, can also be successfully accommodated in our scenario. We can see from Fig. 1 that many samples of (λ, κ) can survive the EWSB conditions in NMSSM. In contrast to the numerical results of TeV mirage mediation in [36,37] within which the allowed (λ, κ) only take values near (0.7,0.1), some portion of (λ, κ) parameter space can survive all constraints in Type A. We should note that such a difference has the following reasons -The choices of the modular weights in [36,37], for example, the values of m Q 3 L and m t c L etc, are different to ours which are given in Eq. (3.9).
-New ingredients, such as the messenger sector which can cause additional deflection of the RGE trajectory, will introduce new free parameters. After all, the mirage mediation scenario can be seen as a special case of our scenario with gauge couplings (Yukawa couplings) being switched off. -The mirage unification scale is set to lie at TeV scale in [36,37]. It is known that simple mirage unification for soft parameters would in general be spoiled in deflected mirage mediation scenario. However, 'mirage' unification for gaugino masses persists which can be proven in our Wilsonian approach (see Appendix A for details). In our scenario, the mirage unification scale for gaugino masses is not constrained to lie at TeV scale. However, we can see from the mirage unification scales shown by different colors in Fig. 1 that even if such scales are required to lie at TeV scale, vast parameter space, which is larger than the numerical result of Refs. ( [36,37]), can survive the constraints from (I-V).
Note that the vacuum stability bounds are also taken into account in our numerical studies, which impose stringent constraints on scenario in [36,37]. In Type B in which the 125 GeV Higgs is the second lightest CP-even scalar, the allowed (λ, κ) parameter space is also much bigger than that in [36,37]. The Higgs mass can be increased by 8 GeV through the mixing with the singlet component for large tan β and λ 0.04. It can be seen from Fig. 2 that the modulus mediation contribution M 0 is bounded to lie between 1 and 8.5 TeV for Type A. A small M 0 always prefers a low messenger scale M mess . A large M 0 , which controls the whole soft SUSY breaking parameters to be heavy, can easily accommodate the SM-like Higgs mass because of large loop corrections from heavy stops in addition to the treelevel contributions involving λ. The value of M 0 is upper bounded to be less than about 5.5 TeV for Type B, which sets an upper bound for the soft SUSY breaking param- The mirage unification scales for gaugino masses are also shown in different colors. All points satisfy the constraints (I-V) eters, especially for the gluino masses. The gluino mass, which is determined by the scale of M 0 , is bounded to below 16 TeV for Type A and below 8 TeV for Type B, respectively. • The Barbieri-Giudice fine-tuning (BGFT) [70] measure is defined as where 'a i ' stands for the set of parameters defined at the input scale. There are two mass scales for the soft SUSY parameters in our scenario, one is the scale that characterize the modulus contributions M 0 and the other is the scale that characterize the anomaly contribution F φ . The latter one is rewritten into a dimensionless quantity α by Eq. (3.13). We can also calculate the electroweak fine-tuning (EWFT) measure EW of the survived points defined in [71,72] EW ≡ max with the relevant terms (see [71,72]) We can see from Fig. 2 that the BGFT measure in our scenario can range from O(1) to O(1000). In fact, the lowest BGFT can reach O(1) for Type B. Such a low fine-tuning possibly indicates that our scenario is natural. We also compare the BGFT measure of the survived points with their corresponding EWFT measure in the last two panels of Fig. 2. It can be seen that the calculated EWFT measure is positively correlated to the corresponding BGFT measure. In most of allowed parameter space, the EWFT and BGFT measures take values of the same order. Besides, the EWFT measure, which can be thought of as providing a lower bound on the electroweak fine-tuning, is always smaller than that of the BGFT measure [73].
In some of the allowed region, the EWFT measure as well as BGFT measure is rather low (being O(1)). As emphasized in [73], EWFT measure is a necessary, albeit not sufficient, measure of electroweak fine-tuning. Low EW need not necessarily mean the model is not finetuned. Rather, it may indicate the possibility that some model might exist with low fine-tuning which might be hidden by the naive application of BG . It is known that low fine-tuning needs light stops as well as a small effective μ, which are naively determined by the dimensional parameter M 0 that controls the whole soft SUSY spectrum with moderate values of α. The lower the M 0 (consequently the lower gluino mass), the lower value of the BGFT(EWFT) measure.
• A positive deflection parameter 'd' is always favored to solve the tachyonic slepton problem in the deflected AMSB for fewer messenger species. In the deflected mirage mediation scenarios, if the modulus contribution is subdominant, a realistic model still prefers a positive deflection parameter 'd' with less messenger species. As the parameter α ≡ α /4 determines the relative size of the contributions between the anomaly mediation and the modulus mediation, a large value of α, which indicates small modulus mediation contributions, needs a large positive deflection parameter 'd' to avoid tachyonic sleptons. We check that large negative values of 'd' are mostly ruled out by the EWSB condition and tachyonic sfermions. It is obvious in the left panel of Fig. 3 that the deflection parameter is constrained to lie at about 1.8 to tune the tachyonic slepton masses to positive values by additional gauge and Yukawa mediation contributions in the region with a large α. It can be seen in Fig. 4 that the corresponding lightest CP-even Higgs mass should lie at a very narrow band centered at about 122.1 GeV with d ≈ 1.8. Besides, the second-lightest CP-odd scalar a 2 is constrained to lie near 6000 GeV if the lightest CPodd scalar a 1 is lighter than 5000 GeV while the lightest CP-odd scalar is constrained to lie near 6000 GeV if the second-lightest CP-odd scalar is heavier than 6000 GeV.
In our scenario, the quantity '4αd' can approximately measure the relative size of deflection contributions (by gauge or Yukawa mediation) to the modulus mediation contributions. We can see from Fig. 3 that the deflection contributions can be dominant in a large portion of the surviving parameter space. Besides, it can also be seen from Fig. 3 that in the modulus mediation dominated regions, that is small α with d = 0, realistic NMSSM spectrums can be possible. This indicates that pure mirage mediation without deflection, which is a special case of our scenarios, can lead to realistic NMSSM spectrum even though it is stringently constrained by EWSB conditions and 125 GeV Higgs. This conclusion agrees with that of [36,37]. Additional deflection with d = 0 from messenger sector can enlarge the possible choice of α in mirage mediation, rendering the mirage-type scenarios more natural. • From Eq. (3.10) in the appendix, we can see the gaugino ratio at the EW scale where g 1 , g 2 , g 3 take values at the GUT scale. On the other hand, the singlino mass is determined by κ and s , which rescales the effective μ e f f parameter by a factor 2κ/λ. In general, a pure singlino-like LSP tends to have a too large relic density due to a comparatively small annihilation cross section because of its small couplings to SM particles.  (A 1 ) couplings. Co-annihilation with heavier χ 0 2 (for m 10GeV) will also efficiently reduce the singlino relic abundance to a proper DM . Besides, the annihilation channels χ 0 1 χ 0 1 → tt, bb can also be important. Similar DM annihilation channels exist for Type B in which the neutralino DM is mainly singlinolike with non-negligible higgsino components. We know that mixed bino-Higgsino DM is severely constrained by direct detection constraints. For a singlinodominant DM, the exchange of a light H 1 can possibly lead to a large direct detection cross section that will be accessible in the present generation of detectors. It can be seen from the lower panels of Fig. 5 that only a small portion of DM parameter space can survive the spinindependent (SI) DM direct detection constraints from the LUX [74], PANDAX [75] and Xenon1T [76]. In fact, direct detection constrained the DM mass to lie in the range [120, 470] GeV for Type A and [50, 400] GeV for Type B, respectively. • Figure 6 shows the SUSY contributions to the muon g − 2. It is known that the required SUSY contributions to a μ can be achieved only if the relevant sparticles(μ,ν μ ,B,W ,H ) are lighter than 600-700 GeV for tan β ∼ 10 in the MSSM [77]. The inclusion of the singlino in the NMSSM can not give sizable contributions to a μ because of the suppressed couplings of singlino to MSSM sector. Although the two loop contributions involving the Higgs are negligible in the SM, new Higgs bosons in the NMSSM could have an important impact on a μ if the lightest neutral CP-odd Higgs scalar is very light [78]. In fact, a positive two-loop contribution is numerically more important for a light CP-odd Higgs being a bit heavier than 3 GeV and the sum of both one-loop and two-loop contributions is maximal around m a 1 ∼ 6 GeV. In our scenario, the lightest CP-odd Higgs a 1 is bounded to be heavier than 40 GeV and give negligible two-loop contributions to a μ in Type A. The main contribution to a μ is thus similar to that in the MSSM. In Type B, the lightest CP-odd Higgs a 1 can lie near 10 GeV and will give important two-loop contributions to increase a μ to values favored by experiments. • Before we finish our discussions on numerical results, we note that the gauge-modulus interference contribution for soft scalar masses will play an important role in phenomenological studies. In our scenario, such a contribution is non-vanishing only for S which is given in Eq. (3.65). This contribution will possibly change the EWSB condition, the dark matter relic density and other collider predictions. We show several benchmark points which are affected by the contributions of T X P (mS 2 ). Table 1 shows several benchmark points, which should not survive the various constraints if the gauge-modulus interference contribution is absent.

Conclusions
We propose to generate a realistic soft SUSY breaking spectrum for Next-to-Minimal Supersymmetric Standard Model with a generalized deflected mirage mediation scenario, in which additional Yukawa and gauge mediation contributions are included to deflect the RGE trajectory. Based on the Wilsonian effective action obtained by integrating out the messengers, the NMSSM soft SUSY breaking spectrum can be given analytically at the messenger scale. We find that additional contributions to m 2 S can possibly ameliorate the