Solving The Muon g-2 Anomaly With Natural NMSSM From Generalized Deflected AMSB

We propose to realize (natural) NMSSM spectrum from deflected AMSB with new messenger-matter interactions. With additional messenger-matter interactions involving ${\bf 10}\oplus{\bf \overline{10}}$ representation messengers, the muon g-2 anomaly can be solved at $2\sigma$ (or $3\sigma$) level with the corresponding gluino mass range $2.8~{\rm TeV}<m_{\tilde{g}}<5.4~ {\rm TeV}$ (or $2.6 ~{\rm TeV}<m_{\tilde{g}}<7.3~ {\rm TeV}$). Besides, our scenario is fairly natural within which the involved fine tuning can be as low as 47. So, in the framework of AMSB-type scenarios, NMSSM can be advantageous in explaining the muon g-2 anomaly in compare with MSSM.


Introduction
Low energy supersymmetric (SUSY) theories, which was regarded for a long time as one of the most appealing extensions of the Standard Model(SM), are well motivated to understand the large hierarchy between the electroweak(EW) scale and the Planck scale. The observed 125 GeV scalar discovered by both the ATLAS [1] and CMS collaborations [2] of LHC, although be compatible with a SM Higgs boson within experimental and theoretical uncertainties, can also be interpreted as the CP-even Higgs boson of the minimal supersymmetric standard model (MSSM). Besides, there are many other hints of low energy SUSY: the observed Higgs mass falls miraculously within the narrow 115 − 135 GeV ′ window ′ required by the MSSM; the mass of top quark mass also lies within the range required by SUSY to radiatively break the EW gauge symmetry; gauge coupling unification, which is not exact in SM, can be successfully realized in the framework of low energy SUSY; the lightest SUSY particle (LSP), which is stable under R-parity, can serve as a natural dark matter candidate with the right order of DM relic density.
It is well known that, to achieve the 125 GeV Higgs in the MSSM, we need either top squarks of order 10 TeV with small mixing or TeV-scale top squarks with large mixing, reintroducing some amount of fine-tuning in the MSSM. This difficulty can be ameliorated in the framework of the Next-to-Minimal Supersymmetric Standard Model(NMSSM). The NMSSM, which is the simplest gauge singlet extension of the MSSM [3], has various preferable features and had drawn a lot of attention in recent years. It can elegantly solve the µ-problem in the MSSM by generating an effective µ-term once the singlet scalar acquires a vacuum expectation value (VEV). Furthermore, due to possible new tree level contributions to the Higgs mass, the NMSSM can easily accommodate the 125 GeV Higgs boson, ameliorating the fine-tuning required.
However, the lack of significant hints in EW precision tests and the null observation of conclusive signals for superparticles at the LHC seem to make the possible discovery of weak scale SUSY doubtable. Recent analysis by LHC have set rather strong constraints on the sparticle masses within the context of various simplified SUSY models: the gluino mg > 1.5 ∼ 1.9 TeV [4] and the lighter top squark mt 1 > 0.85 TeV [5] with even stronger limits on the first generation squarks. Given the rapid accumulation of luminosity by LHC, a 2 TeV gluino is expected to be probed in some SUSY model. Besides, in order to account for the anomaly of muon anomalous magnetic moment in Brookheaven experiments, the weak scale SUSY spectrum should display an intricate structure. More than 3σ discrepancy between the SM prediction [6] and the E821 measurement on [7] g µ − 2 requires the existence of relatively light sleptons and EW gaugino [8]. So supersymmetry breaking mechanism, which determines the whole low energy SUSY spectrum, is thus very crucial to understand such non-generic spectrum.
There are many interesting ways to mediate the SUSY breaking effects in the hidden sector to the visible MSSM sector. The anomaly mediation [9] SUSY breaking mechanisms (AMSB), which is determined by one parameter F φ ≃ m 3/2 , is insensitive to the UV theory [10] and predicts a flavor conservation soft SUSY breaking spectrum. Unfortunately, the minimal AMSB predicts negative slepton masses and must be extended. The most el-egant solution to tachyonic slepton is the deflected AMSB [11] scenario in which additional messenger sectors are introduced to deflect the AMSB trajectory and lead to positive slepton mass with additional gauge mediation contribution. On the other hand, N ≥ 4 species are needed to give positive slepton masses and 125 GeV Higgs with naturally small negative deflection parameters, possibly leading to strong gauge couplings below GUT scale (or Landau pole below Planck scale). So less messenger species are preferable in deflect AMSB. We propose a way to reduce the number of messenger species by introducing general messengermatter interactions in deflected AMSB [12] with small deflection parameter. The low energy NMSSM from typical SUSY breaking mechanism, such as GMSB, is always bothered by the requirement to achieve successful EWSB with suppressed trilinear couplings A κ , A λ and m 2 S , rendering the model building non-trivial [13]. Such difficulty can be ameliorated in AMSB-type scenarios with enhanced trilinear couplings. We find that the phenomenological interesting NMSSM spectrum can be successfully generated with generalized deflected AMSB mechanism involving messenger-matter interactions.
This paper is organized as follows. In Sec 2, we propose our model and discuss the general expressions for soft SUSY parameters from deflected AMSB with general messengermatter interactions. Relevant numerical results are studied in Sec 3. Sec 4 contains our conclusions.

The NMSSM From Deflected AMSB With Messenger-Matter Interactions
It is known that standard GMSB or AMSB can not give a viable NMSSM spectrum. Lacking gauge interactions for S, the A λ , A κ -terms are typically very suppressed in GMSB. In AMSB, large A λ , A κ needs large λ and κ so as to induce large positive m 2 S , suppressing the singlet VEV [14]. By introducing typical form of the messenger sector with messengermatter interactions, negative m 2 S and large A λ , A κ can be generated in our scenario. The framework of the NMSSM includes an extra gauge singlet superfield S in addition to the MSSM sector. In the Z 3 -invariant NMSSM, we have a scale invariant superpotential without the linear, bilinear terms for S and the explicit µ-term. Superpotential of the Z 3 -invariant NMSSM is given by The soft SUSY breaking parameters are given as We propose to generate the soft SUSY breaking parameters of the NMSSM in the framework of deflected AMSB mechanism. The form of the superpotential at the GUT scale (upon the messenger threshold) includes the messenger sectors and superpotential W (X) for pseudo-moduli ′ X ′ with Here we introduce ′ N ′ family of vector-like messengers in the 10 ⊕ 10 representations of SU (5) GUT group Vector-like 5 ⊕5 messengers can also be possibly introduced. We fit the standard model matter contents in the5 a and 10 a representation of SU (5). In the superpotential, the messenger-matter interactions forẼ c L,a sleptons are introduced to give new contributions to the soft masses for such sleptons. Note that it is also possible to introduce the interactions forL L;a (with interactions for5 messengers) or both. Possible massive messenger fields can be included in the superpotential without affecting the AMSB trajectory upon the threshold determined by the pseudo-moduli. It is well known that possible domain wall problems will arise in the Z 3 -symmetric NMSSM. We assume that the Z 3 symmetry is broken by some higher dimensional operators.
Deflection from AMSB trajectory can be used to solve the tachyonic slepton problem. The deflection parameter, which characterize the deviation from the ordinary AMSB trajectory, will depend on the form of the pseudo-moduli superpotential W (X). Such superpotential can be fairly generic and lead to a deflection parameter ′ d ′ of either sign given by In many circumstances, positive deflection parameter [15] could be welcome since less messenger species are needed. Positive deflection parameter in AMSB can be realized by a carefully chosen superpotential [15] or from strong dynamics [16]. After Renormalization Group Equation(RGE) running to the messenger scale set by the pseudo-moduli VEV, the messenger sector reduces For simply, we adopt the simplest case with one family of 10 ⊕ 10 messengers in our scenario. It is known that one additional vector-like 10 representation messengers will not spoil perturbative gauge coupling unification. We also adopt the case with ′ a = 1, 2, 3 ′ which means that the messenger-matter interactions exist for all the three generations.
The messenger-matter interactions involving the gauge singlet S could change its AMSB trajectory and at the same time avoid possible mixing between S and X. The soft SUSY breaking parameters in deflected AMSB can be given by a typical mixed anomaly and gauge mediation contribution. After RGE running down below the messenger threshold, the anomaly part will receive threshold corrections after integrating out the messenger fields. We use the wavefunction renormalization approach [17,18] and replace the messenger threshold M mess and the RGE scale µ by spurious chiral fields X with in the wavefunction superfield and gauge kinetic superfield to obtain the soft SUSY parameters. We reformulate some of the previous results below [12,19,20]. The soft gaugino mass is given at the messenger scale by The trilinear soft terms will be determined by the superpotential after replacing canonical normalized superfields and given by In our convention, the anomalous dimension are expressed in the holomorphic basis [21] (2.12) with ∆G ≡ G + − G − the discontinuity across the messenger threshold. Here G + (G − ) denote respectively the value upon (below) the messenger threshold. The soft scalar masses are given by Details of the expression involving the derivative of ln X can be found in our previous works [12,19,20]. We can see that there are several origins of the contributions in addition to the AMSB contributions: • The gauge-anomaly interference contributions: • Pure gauge mediation contributions With additional interactions for the gauge singlet superfield S, the soft SUSY breaking scalar masses m 2 S and trilinear coupling A κ , A λ will receive new contributions so as to ameliorate the apparent conflict between the requirements of negative m 2 S and large A κ , A λ in ordinary AMSB.
From the previous analytic results for the general deflected AMSB scenario, the soft SUSY breaking parameters in the NMSSM at the messenger scale after integrating out the messenger fields can be given explicitly. The relevant tedious expressions can be found in the appendix.

Numerical Results
EWSB condition is one of the most important constraints for the top-down realization of the NMSSM spectrum. From the minimization conditions of the NMSSM Higgs potential, one of the EWSB condition in the NMSSM can be written as

It is clear from such EWSB condition that m 2
Hu , µ ef f should be light and of order M 2 Z /2 to avoid large fine tuning. On the other hand, the most important corrections Σ u u to m 2 Hu come from the one loop diagrams involving the top squarks and top quark. So, in order to guarantee small loop corrections to m 2 Hu , relatively light stops are needed to accommodate EW naturalness. In the MSSM, light stops should excess 1 TeV (or even heavier without stop maximal mixing) to accommodate the observed 125 GeV Higgs. In the NMSSM, large loop corrections involving stops are not necessary [22] for sizable λ 0.5 and low tan β. The main constraints on light stops come from the LHC direct search constraints.
The messenger-matter interactions in our scenario are chosen to be universal for simply In principle, their precise values should be given by RGE running from the GUT scale to the messenger scale. However, the total messenger contents upon the messenger scale can be model dependent because of the ′ decoupling theorem ′ [10] in AMSB which state that simple messenger threshold (by pure mass term) will not deflect the AMSB trajectory at leading order. By assigning different origin for messenger thresholds (determined by moduli VEV or pure mass term), effects of certain messengers on low energy AMSB spectrum can decouple although they can still contribute to gauge coupling(Yukawa) unification. Because of the freedom to add messengers (with pure mass thresholds) upon messenger scale, we neglect the RGE effects for messenger-matter interactions in our numerical study for simply. The free parameters in our scenario are: We need to check if successful EWSB condition is indeed fulfilled. In fact, the soft SUSY mass m 2 Hu , m 2 H d , m 2 S can be reformulated into µ, tan β, M 2 Z by the minimum condition of the scalar potential. Usually, M A can be used to replace A κ by We should note that κ is a free parameter while tan β is not. This choice is different to ordinary numerical setting in the NMSSM in which tan β is free while κ is a derived quantity [23]. Our choice can be convenient for those predictable NMSSM models from topdown approach. A guess of tan β is made to obtain the relevant Yukawa y t , y b couplings at the EW scale. After RGE evolving up to the messenger scale as the theory inputs, the whole soft SUSY breaking parameters at the messenger scale can be obtained. Low energy tan β can be obtained iteratively from such spectrum with minimization condition for the Higgs potential. The purpose of deflection in AMSB is to solve the notorious tachyonic slepton problem. So non-tachyonic slepton should be obtained at the SUSY scale. Positive slepton masses can be realized by introducing either sign of deflection parameter d with proper mattermessenger interactions.
We use the package NMSSMTools 5.01 [24] to scan the whole parameter space. We will interest in relatively large values of λ in order to increase the tree-level mass of the 125 GeV CP-even Higgs boson. The parameters are chosen to satisfy: Besides, we require that λ 2 + κ 2 0.7 to satisfy the perturbative bounds. In our scan, we impose the following constraints: • The lower bounds from current LHC constraints on SUSY particles [4,5]: -Light stop mass: mt 1 0.85 TeV.
• The lower bounds for 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.
• Flavor constraints from the rare decays of B mesons. We adopt the recent experimental results [25]: • The CP-even component S 2 in the Goldstone-′ eaten ′ combination of H u and H d doublets corresponds to the SM Higgs. The S 2 dominated CP-even scalar should lie in the combined mass range for the Higgs boson: 122GeV < m h 0 < 127GeV from ATLAS and CMS data. Note that the uncertainty is 3 GeV instead of default 2 GeV because large λ may induce additional O(1) GeV correction to m h 0 at two-loop level [26], which is not included in the NMSSMTools.
• The relic density of DM should satisfy the Planck data Ω DM = 0.1199 ± 0.0027 [27] in combination with the WMAP data [28](with a 10% theoretical uncertainty). In our scenario, only the upper bound for DM relic density is used.
• The g µ −2 discrepancy should be solved in our scenario. The E821 experimental result on the muon anomalous magnetic moment at the Brookhaven AGS [29] is given as which is larger than the SM prediction The deviation is about 3σ ∆a µ (expt − SM) = (255 ± 80) × 10 −11 . (3.8) Recent discussions on g µ − 2 anomaly can be seen in [30].
We have the following discussions on our numerical results: • Successful EWSB condition in the NMSSM always set stringent constraints on the constrained low energy inputs, especially when such low energy inputs are determined by certain UV-completed theory. In our scenario, the (deflected) AMSB-type mechanism with few free parameters determines the whole SUSY spectrum at the messenger scale. So successful EWSB condition will rule out much parameter space. Random scan in our scenario indicates that many points can still survive the EWSB conditions. The allowed range for the characteristic NMSSM parameters κ and λ can be seen in fig.1. We can see that our numerical results give the lower bound for κ and λ with κ low ≈ 0.43 and λ low ≈ 0.35 in addition to the upper bounds κ, λ ∼ 0.62. • The left plot of fig.2 shows the SUSY contributions to g µ − 2 anomalous magnetic momentum ∆a µ versus the gluino mass mg. It is obvious that the g µ − 2 anomaly can be successfully solved at both the 2σ(green points) and 3σ (grey points) level. The MSSM is known to be able to yield sizable contributions to the ∆a µ which come dominantly from the chargino-sneutrino and the neutralino-smuon loops. The g µ − 2 anomaly, of order  The inclusion of singlino in the NMSSM will not give sizable contributions to ∆a µ because of the suppressed coupling of singlino to the MSSM sector. Two loop contributions involving the Higgs is negligible in SM. In fact, the SM one-loop Higgs/muon diagram is, about four orders of magnitude below the sensitivity of the g µ − 2 experiment. However, the new Higgs bosons could have an important impact on a µ if the lightest neutral CP-odd Higgs scalar can be quite light [31]. As noted in that paper, the positive two-loop contribution is numerically more important for a light CP-odd Higgs heavier than 3 GeV and the sum of both one loop and two loop contributions is maximal around m a 1 ∼ 6 GeV. We show the CP-even and CP-odd Higgs masses in fig.3 for the allowed parameter space by ∆a µ . From the left plot of fig.3, we can see that the lightest CP-odd Higgs a 1 , which always lies above 62 GeV, is not light enough to give sizable contributions to ∆a µ . The main contributions are similar to that in the MSSM. We should also note that the 125 GeV SM-like Higgs always corresponds to the lightest CP-even Higgs boson for our survived points. The possibility that the 125 GeV Higgs being the second lightest CP-even scalar is not given in our scenario. The LHC discovery prospects for the scalar sector of the NMSSM can be seen in [32][33][34][35][36].
• In most of the allowed parameter spaces, the lightest CP-odd scalar a 1 and the second lightest CP-even scalar h 2 are singlet dominated. The corresponding singlet components of the a 1 and h 2 are displayed in the lower panels of fig.3. If kinetically allowed, the h 2 scalar can decay primarily into a 1 scalar pairs and vector boson pairs.
The cross section of h 2 production from dominant gluon fusion is small in compare with the SM-like Higgs because of its dominant singlet component. In case the pri- mary h 2 → a 1 a 1 decay channel is open, the W + W − , ZZ final states from gluon fusion produced h 2 scalar decay are typically suppressed by order 10 −5 in compare with the SM-like Higgs of the same mass, making it difficult to give a signal of statistical significance at the LHC. On the other hand, in a small region of parameter space within which the h 2 → a 1 a 1 channel is not open for light h 2 and heavy a 1 , the h 2 can decay primarily into W + W − pairs with llνν(l = e, µ), lνjj final states. Because of the reduced production rate and large background from pp → W + W − and tt → bbW + W − , they could hardly be discovered by recent LHC experiments. The possibility to detect the singlet-dominated CP-even scalar in diphoton channel had been discussed in [32]. However, in our scenario, the branch ratio into diphoton is very small in compare with that for a SM Higgs boson of the same mass, making the diphoton final states not very promising to search for at √ s = 13TeV LHC.
From our numerical results, the CP-odd scalar a 1 should be heavier than 62 GeV (half of the SM-like Higgs mass). If this scalar is lighter than half of the SM-like Higgs boson, the SM-like Higgs boson could decay into a 1 pairs. Because the width of the SM-like Higgs boson is quite narrow, such an exotic decay may have a sizable branching ratio and in turn suppress greatly the visible signals of the SM-like Higgs boson at the LHC. So the LHC measurements, in particular those that are put on the Higgs exotic decay modes, will give stringent constraints on such light CP-odd scalar.
The a 1 scalar can be produced from heavy scalar cascade decay or direct production in association with a bottom-quark pair [37]. Each a 1 from h 2 → a 1 a 1 will decay primarily to bb and τ + τ − , yielding final states that will typically have large backgrounds at the LHC. Requiring one a 1 to decay to µ + µ − pairs can allow one to use the dimuon invariant mass as a discrimination between the signal and the background. However, the σ × Br cross section is rather small. Possible 4γ signal from both the diphoton decay modes of a 1 will not be significant enough due to the absence of a possible enhancement of the diphoton branching fraction. Two benchmark points are shown in Table 1 with the decay modes for light scalars and the h 2 production cross section from gluon fusion. • It is also clear from fig.2 that the gluino should lie between 2.6 TeV to 7.3 TeV to solve the observed g µ −2 anomaly. The upper bound on gluino masses can be understood as follows: The scale of the low energy soft SUSY parameters from anomaly mediation is determined by the value of F φ . The lightness ofμ,ν µ ,B,W to account for the g µ − 2 anomaly sets an upper bound on F φ , which subsequently sets the mass scale of the gluino.
The observed SM-like Higgs can also be explained in our scenario as anticipated (see the right plot of fig.2) because additional contributions from the NMSSM can increase the predicted SM-like Higgs mass. Besides, in order to solve the g µ − 2 anomaly at 2σ(3σ) level, the Higgs mass can be as high as 125.6GeV (126.7 GeV). As a comparison, the Higgs mass is upper bounded by 118 GeV(120 GeV) if the g µ − 2 anomaly is solved at 2σ(3σ) level in the CMSSM. So, our scenario is much better in solving the g µ − 2 anomaly and satisfying the Higgs mass measurement.
In the framework of (extended) anomaly mediation, the NMSSM seems to be advan-tageous in compare with the MSSM in which the gluino masses are bounded to below 2.0 TeV(2.5 TeV) if the g µ − 2 anomaly is solved at 2σ (3σ) level with complete GUT multiplet messengers [20]. In fact, the latest LHC search results had already set a lower bound 1.5 ∼ 1.9 TeV on gluino mass. In most of the allowed parameter space, the MSSM from (deflected) AMSB can not solve the g µ − 2 anomaly at 2σ level unless incomplete GUT multiplet messengers are introduced. Although the two scenarios in the MSSM and the NMSSM introduce different types of deflections in their messenger sectors, an important constraint on the parameters comes from the 125 GeV Higgs mass.
In the MSSM, the observed 125GeV Higgs mass will set stringent constraints on the soft SUSY breaking parameters, especially on the trilinear coupling A t and stop masses mt 1 , mt 2 . So the parameter F φ ≃ m 3/2 , which determines the whole soft SUSY spectrum, is bounded to lie in a small region (in combing with other constraints). Thus, the gaugino masses, which are fully determined by F φ and d, are constrained to lie below 2.5 TeV. On the other hand, the observed 125 GeV Higgs can easily be accommodated in the NMSSM and will not set very stringent constraints on the parameter F φ . So the gaugino masses, which are predicted by the value of F φ , can be relatively heavy in the NMSSM.
The deflection with non-trivial messenger-matter interactions involvingμ c L is crucial in obtaining sizable contributions to ∆a µ . We can see from fig.4 that 1.4 < λ 0 < 2 with positive ′ d ′ in the range 0.45 < d < 0.85 is required for all allowed points in our scenario. Positive deflection in this range can cure the tachyonic slepton problem of ordinary AMSB and give small positive slepton masses to account for sizable ∆a µ . Without such new messenger-matter interactions in AMSB type scenario, the two problems can not be solved simultaneously. More 10 ⊕ 10 messenger species will change the allowed range of d. Similar conclusion holds for scenarios with messengermatter interactions involving the doublet (ν µ;L ,μ L ).
• As noted earlier, relatively light stops mt 1 ,t 2 1.5 TeV are needed by naturalness with EW fine tuning ∆ EW 10 in the (radiative) natural MSSM [39]. Heavy stops are not necessarily required to interpret the 125 GeV Higgs in the NMSSM because of the additional tree-level contributions. So it seems natural to achieve (electroweak) naturalness in the NMSSM instead of the MSSM. Although new LHC search results had already imposed stringent constraints on the mass of top squarks to excess 1 TeV, EW naturalness with O (TeV) stops can still be possible with large |A t |. Detailed studies on the radiative natural NMSSM will appear soon in our subsequent studies [40].
The fine-tuning with respect to certain input parameter a is defined as while the total fine-tuning is defined to be ∆ = max a (∆ a ) with {a} the set of param- eters defined at the input scale.
The fine-tuning of our scenario within the allowed parameter space by ∆a µ , the LHC results and the EW precision measurement is shown in fig.5. In the allowed region, the fine tuning ∆ 360. The fine tuning can be as low as 47 in certain regions which can be seen to generate a natural NMSSM spectrum. All the trilinear couplings A t are negative with 1.6TeV < mt 1 < 5.8TeV. As expected, the low fine-tuning region corresponds to lightt 1 . Besides, for a given A t , lightert 1 always corresponds to less fine tuning. Similar conclusion in the MSSM can be seen in [39].
• The lightest neutralino can act as an excellent DM candidate in the NMSSM. Our numerical scan (see fig.6) indicates that such neutralino DM mass should lie in the range 100GeV < m LSP < 300GeV to account for the g µ − 2 anomaly. The lightest neutralino in the NMSSM should be the mixing of the bino(B), wino(W ), higgsino(H 1,2 ) and singlino (S) 10) with N 1i the matrix elements to diagonalize the neutralino mass matrix. In our scenario, we find that the wino is the dominant component of LSP (see left-plot of Figure 5. Survived points with their lightest stop mass mt 1 versus the trilinear couplings for top Yukawa A t and also their corresponding fine tuning involved. Note that the trilinear coupling A t of all points are negative. fig.6). This result can be partially understood from the gaugino ratios at the EW Figure 6. The left plot shows the matrix elements to diagonalize the neutralino mass matrix N 1i versus the lightest neutralino mass(mχ1 0 ) while the right plot shows the effective µ parameter for higgsino mass µ ef f ≡ λ s versus the lightest neutralino mass(mχ1 0 ). All survived points satisfy the LHC and EW precision measurement constraints as well as the required value of ∆a µ .
scale. The gaugino mass ratios at the weak scale are given approximately by (3.11) Knowing the range of the deflection parameter 0.45 < d < 0.85, the lightest gaugino can be identified from the above formula. The plot for effective µ ≡ λ s versus m LSP can be seen in the right plot of fig.6. We can see that the scale of effective µ is not too large in compare with mχ1 0 . So the higgsino components in neutralino LSP can be sizable. The NMSSM-specific singlino component is negligibly small which will not play an important role in DM annihilation processes. It is well known that wino DM will always lead to under abundance of DM relic abundance unless mW ∼ 3 TeV. This result holds in our scenario (see fig.7). The wino-like DM can not provide full relic abundance with its mass of order O(100) GeV. Other DM components, for example, axion or axino, are necessary to provide enough cosmic DM. We also check that the spin-independent DM direct detection constraints, for example, the LUX2016 [41] and the PandaX [42], are satisfied for all survived points.

Conclusion
We propose to realize the (natural) NMSSM spectrum from deflected AMSB with new messenger-matter interactions. With additional messenger-matter interactions involving 10 ⊕ 10 representation messengers, the muon g-2 anomaly can be solved at 2σ (or 3σ) level with the corresponding gluino mass range 2.8 TeV < mg < 5.4 TeV (or 2.6 TeV < mg < 7.3 TeV). Besides, our scenario is fairly natural within which the involved fine tuning can be as low as 47. So, in the framework of AMSB-type scenarios, the NMSSM can be advantageous to explain the muon g-2 anomaly in compare with the MSSM.
Appendix: Soft SUSY Parameter For the NMSSM From the previous analytic expressions, we can obtain the soft SUSY breaking parameters for the NMSSM at the messenger scale. The gaugino masses are given as The trilinear soft terms are given by , Expressions for scalars are rather complicate, they can be parameterized as the sum of each contributions • Pure deflected anomaly mediation contribution δ d without new yuakwa couplings Here N = 3 in our scenario.