Explaining the Muon g-2 Anomaly in Deflected AMSB for NMSSM

We propose to embed the General NMSSM (Next-to-Minimal Supersymmetric Standard Model) into the deflected AMSB (Anomaly Mediated Supersymmetry Breaking) mechanism with Yukawa/gauge deflection contributions. After integrating out the heavy messenger fields, the analytical expressions of the relevant soft SUSY breaking spectrum for General NMSSM at the messenger scale can be calculated. We find that successful EWSB (Electroweak Symmetry Breaking) and realistic low energy NMSSM spectrum can be obtained in some parameter regions. In addition, we find that the muon $g-2$ anomaly and electron $g-2$ anomaly (for positive central value electron $g-2$ experimental data) can be jointly explained to $1\sigma$ and $2\sigma$ range, respectively. The $Z_3$ invariant NMSSM, which corresponds to $\xi_F=0$ in our case, can also jointly explain the muon and electron anomaly to $1\sigma$ and $2\sigma$ range, respectively.


Introduction
The Standard Model (SM) of particle physics is very successful in explaining the vast experimental measurements up to the electroweak (EW) scale, including the discovery of Higgs boson by the Large Hadron Collider (LHC) [1,2]. However, it still has many theoretical and aesthetic problems, for example, the quadratic divergence of the fundamental Higgs scalar mass, the dark matter (DM) puzzle, and the origin of baryon asymmetry in the universe. In addition, it has been known for a long time that the theoretical prediction of the muon anomalous magnetic moment a µ ≡ (g − 2) µ /2 for SM has subtle deviations from the experimental values. In fact, combining the recent reported E989 muon g − 2 measurement with the previous BNL result [3,4], the updated world average experimental value of a µ is given by [5] a FNAL+BNL µ = (11659206.2 ± 4.1) × 10 −10 , (1.1) which has a 4.2σ deviation from the SM prediction [6] ∆a FNAL+BNL µ = (25.1 ± 5.9) × 10 −10 .
(1. 2) In addition to the reported muon g − 2 anomaly, the experimental data on electron g − 2 also reported some deviations from the SM predictions. From the measurement of the fine structure constant α em (Cs) by the Berkeley experiment using 133 Cs atoms [7], the experimental value on electron g − 2 by [8] has a 2.4σ deviation from the SM prediction [9] ∆a Exp with the center value taking the negative sign. Such a result is not consistent with the most accurate 2020 measurement using 87 Rb atoms [10], which reported a 1. with the center value taking the positive sign. Although the deviation of electron g−2 is still controversial, its possible theoretical implications should not be overlooked. The previous problems and anomalies strongly indicate that SM should not be the whole story and it only acts as the low energy effective theory of some new physics beyond the SM. Various new physics models had been proposed to deal with the problems that bother the SM. Among them, low-energy supersymmetry (SUSY) is the most attractive one, which can accommodate almost all the solutions of such problems together in a single framework. In particular, the discovered 125 GeV Higgs scalar lies miraculously in the small '115 − 135' GeV window predicted by the low energy SUSY, which is a strong hint of weak scale SUSY. So, if low energy SUSY is indeed the new physics beyond the SM, it should account for the new muon g − 2 anomaly. SUSY explanations of the muon g − 2 anomaly can be seen in the literatures . Furthermore, there have already been several discussions offering combined explanations of the experimental results for electron and muon g − 2 anomaly in the SUSY framework [39][40][41][42][43][44][45][46][47]. On the other hand, it is rather non-trivial for the low-energy SUSY breaking spectrum to be consistent with recent experimental bounds, for example, the LHC exclusion bounds and DM null search results. As the low-energy soft SUSY breaking parameters are fully determined by the predictive SUSY breaking mechanism, their intricate structures can be a consequence of the mediation mechanism of SUSY breaking in the UV theory, for example, the well-motivated anomaly mediated SUSY breaking (AMSB) [48,49] mechanism.
Minimal AMSB, which is determined solely by the F-term VEV of the compensator superfield F φ (its value is approximately equal to the gravitino mass m 3/2 ), is insensitive to its UV theory [50] and predicts a flavor conservation soft SUSY breaking spectrum. Unfortunately, negative slepton squared masses are predicted and the minimal scenario must be extended. Although there are many possible ways to tackle such tachyonic slepton masses problems, the most elegant solution from an aesthetic point of view is the deflected AMSB (dAMSB) [51][52][53][54][55] scenario, which adopts non-trivially an additional messenger sector to deflect the AMSB trajectory and push the negative slepton squared masses to positive values by additional gauge mediation contributions.
Next-to-Minimal Supersymmetric Standard Model (NMSSM) [56,57] can elegantly solve the µ problem that bothers the Minimal Supersymmetric Standard Model (MSSM) with an additional singlet sector. In addition, with additional tree-level contributions or through doublet-singlet mixing, NMSSM can accommodate easily the discovered 125 GeV Higgs boson mass. However, soft SUSY breaking parameters of low-energy NMSSM from a typical SUSY breaking mechanism, such as gauge-mediated SUSY breaking (GMSB) [58,59], are 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 [60]. Such difficulties always persist in ordinary AMSB-type scenarios. We find that the phenomenologically interesting NMSSM spectrum with successful EWSB can be successfully generated by combining both AMSB and GMSB-type contributions. In addition, the discrepancy between the theoretical predictions for the muon (and electron) anomalous magnetic momentum and the experiments can be explained in such model. This paper is organized as follows. In Section 2, we propose our model and discuss the general expression for soft SUSY parameters. The soft SUSY parameters for General NMSSM are given in our scenario. The relevant numerical results are studied in Section 3. Section 4 contains our conclusions.

Soft SUSY Breaking Parameters of NMSSM from Deflected AMSB
We know that NMSSM is well motivated theoretically to solve the µ − Bµ problem. By imposing the discrete Z 3 symmetry, a bare µ term is forbidden in the NMSSM. An effective µ ef f parameter can be generated by after the S field acquires the VEV s , which lies of order the electroweak scale and breaks the Z 3 symmetry. The Z 3 invariant superpotential couplings are given by [56,57] The soft SUSY breaking parameters are given as The breaking of Z 3 discrete symmetry in the scale-invariant NMSSM is bothered with the domain-wall problem, creating unacceptably large anisotropies of the CMB and spoiling successful BBN predictions. In addition, discrete global symmetry cannot be exact at the Planck scale and may be violated by Planck scale suppressed gravitational interactions. Although one can find complicated solutions to such difficulties, it is interesting to go beyond the Z 3 invariant scheme and adopt the general NMSSM (GNMSSM) to evade the previous problems. In fact, effective GNMSSM can also be generated from scale invariant NMSSM if one introduces, in the Jordan frame, a χH u H d term in the frame function that couples to the curvature.
GNMSSM is the most general single gauge singlet chiral superfield S extension of MSSM, including the most general renormalizable couplings in the superpotential and the corresponding soft SUSY breaking terms in L sof t . On the other hand, as noted previously, GNMSSM does not adopt the Z 3 symmetry, which can contain the Z 3 breaking terms [56,57] So, the general superpotential of GNMSSM contains with the corresponding new terms of soft SUSY breaking parameter Such low-energy soft SUSY breaking parameters are determined by the corresponding SUSY breaking mechanism in the UV-completion theory [58,[61][62][63]. However, it is rather non-trivial to generate the phenomenological desirable low-energy soft SUSY breaking spectrum from UV theory. In fact, null search results of sparticles with 139 f b −1 of data at the (13 TeV) LHC by the ATLAS and CMS collaborations [64] suggest that the low-energy SUSY spectrum should have an intricate pattern. For example, the first two generation squarks need to be heavy to avoid the stringent constraints from LHC. To explain the muon g − 2 anomaly, light electroweakinoes and sleptons are always preferred. On the other hand, the discovered 125 GeV Higgs boson by both the ATLAS and CMS collaborations of LHC may indicate stop masses of order 5∼10 TeV or TeV scale stops with large trilinear coupling A t . As light stops can be preferred to keep electroweak naturalness, large trilinear coupling A t is, therefore, needed to accommodate the 125 GeV Higgs. We know that minimal GMSB predicts vanishing A t at the messenger scale and introducing messenger-matter interactions will sometimes be bothered with additional flavor constraints. So, we propose to generate the phenomenological desirable low-energy NMSSM soft spectrum in the predictive AMSB-type framework with additional Yukawa mediation contributions.
As the slepton sector of NMSSM is the same as that of MSSM, the prediction of the NMSSM soft spectrum from AMSB is still bothered with the tachyonic slepton problem. The most elegant solution to such a problem is to add an additional messenger sector, which can deflect the AMSB trajectory properly so that the slepton squared masses can be pushed to positive values. The minimal version of deflected AMSB in NMSSM, which is a straightforward extension of ordinary minimal deflected AMSB for MSSM, adopts no additional interactions other than that between the spurion field and the messengers. However, it is still problematic in realizing successful EWSB and triggering a non-vanishing s , as the necessary condition A 2 κ 9m 2 S [56,57] can not be satisfied easily. In addition, our numerical results indicate that the corresponding SUSY contributions to muon g − 2 anomaly is very small for the few survived parameter points. So, additional coupling terms in the superpotential involving the singlet S and messengers can be introduced to bring new Yukawa deflection contributions to the soft SUSY breaking spectrum and alter those parameters relevant for EWSB.
To evade unwanted mixing between the singlet chiral field S and the spurion field X, we adopt the following form of superpotential at the messenger scale with the coupling form first proposed in [65] for GMSB embedding of NMSSM. The soft SUSY breaking parameters in deflected AMSB can be given by mixed anomaly mediation and gauge mediation contributions. When the renormalization group equation (RGE) evolves down below the messenger thresholds, the anomaly mediation contributions will receive threshold corrections after integrating out the messenger fields. The wavefunction renormalization approach [66,67] can also be used in deflected AMSB [68,69] to keep track of supersymmetry-breaking effects with a manifestly supersymmetric formalism and obtain the soft SUSY breaking spectrum. In deflected AMSB, the messenger threshold M mess and the RGE scale µ in the wavefunction superfield and gauge kinetic superfield can be replaced by the following combinations involving the spurious chiral fields X which, after substituting the F-term VEV of X and φ, can obtain the soft SUSY breaking parameters. Below the messenger scale, the messenger superfields will be integrated out and will deflected the trajectory from ordinary AMSB trajectory. The superpotential for pseudo-moduli superfield W (X) can be fairly generic and leads to a deflection parameter of either sign given by After integrating out the heavy messengers, we can obtain the soft SUSY breaking spectrum at the messenger scale. The soft gaugino masses at the messenger scale can be given by The trilinear soft terms can be determined by the wavefunction renormalization factors In our convention, the anomalous dimension are expressed in the holomorphic basis [70] G (2.14) and ∆G ≡ G + −G − the discontinuity of anomalous dimension across the messenger threshold. Here, 'G + (G − )' denote the values of anomalous dimension above (below) the messenger threshold, respectively. The soft scalar masses are given by (2.16) From the previous analytic expressions, we can calculate the soft SUSY breaking parameters for MSSM at the messenger scale. The gaugino masses are calculated to be with the corresponding beta function and the changes of beta function for the gauge couplings 2N, 2N, 2N ), (2.18) with N representing the family of messengers in (2.8).
The trilinear soft terms are calculated to be 20) and the discontinuity of the Yukawa beta functions across the messenger threshold The sfermion masses can be calculated to be For the third generation, the sfermion masses are given by within which we also include the top, bottom and tau Yukawa contributions. Such contributions can not be neglected at large values of tan β.
Expression of the soft m 2 S receives both anomaly mediation, Yukawa mediation, and interference contributions with the pure anomaly mediation part being 26) and the Yukawa deflection part (including the interference terms) being 3 Joint Explanation of Muon and Electron g − 2 Anomaly We would like to give a joint explanation of muon and electron g − 2 anomalies in the NMSSM framework from deflected AMSB, taking into account all other low-energy experimental data/exclusion bounds. The SUSY contributions to muon g − 2 are dominated by the chargino-sneutrino and the neutralino-smuon loop, and the corresponding Feynman diagrams are shown in Figure 1 [72]. At the leading order of tan β and m W /m SU SY , they are evaluated as [73] ∆a µ (μ L ,μ R ,B) = α Y 4π Here, m µ is the muon mass, m SU SY the SUSY breaking masses and µ the Higgsino mass, respectively. The loop functions are defined as which are monochromatically increasing for x > 0, y > 0 and satisfy 0 ≤ f C,N (x, y) ≤ 1. They satisfy f C (1, 1) = 1/2 and f N (1, 1) = 1/6 in the limit of degenerate masses. The SUSY contributions to the muon g − 2 will be enhanced for small soft SUSY breaking masses and large value of tan β. Similar expressions are evident for electron g − 2 after properly replacing the couplings and mass parameters for muon by those for electrons. The inclusion of the singlino component in NMSSM will not give sizable contributions to ∆a µ because of the suppressed coupling of singlino to the MSSM sector. However, the lightest neutral CP-odd Higgs scalar could give non-negligible contributions to a µ if it is quite light [74]. The positive two-loop contribution is numerically more important for a light CP-odd Higgs at approximately 3 GeV and the sum of both one loop and two loop contributions is maximal around m a 1 ∼6 GeV. It is non-trivial to jointly explain both muon and electron g − 2 anomalies in a single framework. The joint explanation of muon and electron g − 2 anomalies with a negative center value for ∆a e needs either large non-universal trilinear A-terms [47,75] or flavor violating off-diagonal elements in the slepton mass matrices [40,47]. Without explicit flavor mixings, the two anomalies can also be explained by arranging the bino-slepton and chargino-sneutrino contributions differently between the electron and muon sectors, requiring heavy left-hand smuon [39] or light selectrons, wino, and heavy higgsino [42]. Without large flavor violation in the lepton sector, we anticipate that the new physics contributions to the leptonic g − 2 will, in general, scale with the corresponding lepton square masses. Although such scaling solutions cannot explain the electron g − 2 data in (1.3) with negative central value to 2σ range, it, however, can be consistent with the most accurate data in (1.4) with a positive central value. In our model, universal soft SUSY breaking parameters are predicted for the muon and the electron sector at the messenger scale. Two loop RGE running will only split slightly the low energy spectrum for the two sectors at the SUSY scale. Therefore, we anticipate that the scaling solution will approximately hold in our case.
We use NMSSMTools 5.6.2 [61,76] to scan the whole parameter space to find the desired parameter regions that can account for the muon and electron g −2 anomaly. In our numerical calculations, the choices N = 1 and λ X;a = λ X , λ S;a = λ S are adopted for simply. In addition, although there are three new free parameters in the GNMSSM superpotential in comparison to Z 3 invariant NMSSM, we adopt the most predictive choice with only vanishing ξ F and other parametersμ = µ = 0. Non-vanishing ξ F will lead to additional tadpole terms for S and alter the values of s , h u , and h d for the EWSB minima.
The ranges of the input parameters at the messenger scale M mess are chosen to satisfy with F φ M mess . Below the messenger scale, the heavy messenger fields are integrated out and the corresponding RGE trajectory reduces to ordinary GNMSSM type, which can be seen to be deflected from the ordinary AMSB trajectory [48,51] and can, possibly, push the tachyonic slepton masses to positive values. So, the renormalization group equations of GNMSSM [61] are used to evolve the soft SUSY parameters from the messenger scale to the EW scale.
We also impose the following constraints from low energy experimental data other than that already encoded in the NMSSMTools package (i) The CP-even component S 2 in the Goldstone-'eaten' combination of H u and H d doublets corresponds to the SM Higgs. Such an S 2 dominated CP-even scalar should lie in the combined mass range for the Higgs boson, 122GeV < M h < 128GeV [1,2]. Note that the uncertainty is 3 GeV instead of default 2 GeV because large λ may induce additional O(1) GeV correction to M h at the two-loop level [77].
(ii) Direct search bounds for low mass and high mass resonances at LEP, Tevatron, and LHC by using the package HiggsBounds-5.5.0 [78].
A sufficient condition to ensure vacuum stability at the EW scale is the requirement that EW vacuum is the global minimum (true vacuum) of the scalar potential. If the EW vacuum is a local minimum (false vacuum), the relevant parameter regions can still be allowed if the false EW vacuum is meta-stable with a lifetime longer than the age of the universe.
(vi) The relic density of the dark matter should satisfy the Planck result Ω DM h 2 = 0.1199±0.0027 [90] in combination with the WMAP data [91] (with a 10% theoretical uncertainty).
We have the following discussions on our numerical results.
• Although it is fairly non-trivial to realize successful EWSB in NMSSM from predictive UV-completion models, for example in ordinary GMSB, numerical scan indicates that some parameter points can still survive the EWSB conditions in our case. In fact, additional couplings in the superpotential involving the singlet S and messengers can change the AMSB predictions of m 2 S and A λ , A κ so as that the necessary condition A 2 κ 9m 2 S for s = 0 and other EWSB conditions can be satisfied. The values of tan β at the EW scale can be obtained iteratively after we minimize the scalar potential to obtain s . The allowed values of λ, κ and the corresponding µ ef f are shown in the left panel of Figure 2. We can see that the allowed values of λ and κ are always not large. The dependence of ξ F versus the low scale tan β are also shown in the right panel of Figure 2. An interesting observation is that successful EWSB can still be allowed with ξ F = 0, which is just the Z 3 -invariant NMSSM case. Survived points that can satisfy the constraints (i-vi). The allowed ranges of λ versus µ ef f , κ are shown in the left panel while the values of ξ F versus ∆a µ , tan β are shown in the right panel.
• From our numerical results, we can see in the right panel of Figure 3 that the muon g−2 anomaly can be explained to 1σ range. As noted previously, small flavor violation in the lepton sector will predict that ∆a e and ∆a µ satisfy the scaling relation which predicts the same sign of ∆a e as that of ∆a µ . An explanation of the muon g − 2 anomaly can also lead to the explanation of the electron g − 2 anomaly in 2σ range for positive central value electron g − 2 experimental data in (1.3). Figure 3 shows a scatter plot of ∆a e and ∆a µ with the corresponding SM-like Higgs masses in different colors, as the SM-like Higgs mass always exclude a large portion of otherwise allowed parameter regions. From (3.11), due to their dependences on the square of the corresponding lepton masses, ∆a e can be seen to be of order 10 −14 when ∆a µ ∼ O(10 −9 ). However, ∆a e ∼ O(10 −14 ) leads to an apparent horizontal line, when the plot for ∆a e versus ∆a µ is shown in Figure 3. The NMSSM specific contributions to ∆a µ,e are dominantly given by the Barr-Zee type two-loop contributions involving the lightest CP-odd scalar a 1 . However, our numerical results indicate that the a 1 relevant NMSSM specific contributions to ∆a µ,e are always small and subdominant.
The plot of the SUSY contributions to muon anomalous magnetic momentum ∆a µ versus the ξ F parameter are shown in the right panel of Figure 2. It can also be seen from the panel that the Z 3 -invariant NMSSM case, which corresponds to ξ F = 0, can also explain the muon g − 2 anomaly to 1σ range (and the electron g − 2 anomaly to 2σ range by scaling relations).
• The left panel of Figure 4 shows the plot of the SUSY contributions to muon anomalous magnetic momentum ∆a µ versus the gluino mass Mg. In AMSB-type scenarios, the F φ parameter determines the mass scales of all the soft SUSY breaking parameters. The larger the value of F φ , the heavier the sfermion and the gaugino masses. We know that light sleptons and electroweakinoes with masses below 0.5∼1 TeV are preferred to explain the muon/electron g − 2 anomaly via chargino-sneutrino and the neutralino-smuon loops. So, the SUSY explanations of muon/electron g − 2 anomalies prefer smaller F φ , consequently imposing an upper bounds on the low-energy sparticle masses. Our numerical results show that the gluino masses are bounded to lie 3.5 TeV ≤ Mg ≤ 6.0 TeV if the muon g − 2 anomaly is explained upon 3σ range. Gluino masses upon 3.5 TeV can possibly be discovered in the future 100 TeV FCC-hh collider.
In mSUGRA type models with universal gaugino masses at the GUT scale or GMSB type models, the gaugino ratios at the EW scale are always given by  g−2 anomaly rather hard. In our case, the gaugino mass ratios change approximately to M 1 :M 2 :M 3 ≈ (6.6 − 2d):2 (1 − 2d):6 (−3 − 2d) at the EW scale. Therefore, with a proper range of deflection parameter d, the gluino mass can be heavy without contradicting the requirements of light electroweakinoes by the explanation of the muon g − 2 anomaly. We should note that a positive deflection parameter d is always favored to solve the tachyonic slepton problem for few messenger species in deflected AMSB. To tune the slepton squared masses to small positive values, the range of d are constrained to lie in a small range. In fact, our numerical results indicate that the deflection parameters, which parameterizes the relative size between the anomaly mediation contributions and the gauge/Yukawa mediation contributions, are constrained to lie 0.55 < d < 0.9 (see the right panel of Figure 4), allowing the gluino to be heavier than 4 TeV for O(100) GeV wino.
• It can be seen from the previous figures that the observed SM-like 125 GeV Higgs can be accommodated easily in our model. Additional tree-level contributions to SM-like Higgs mass from NMSSM in general allow much lighter stop masses in comparison to MSSM. In addition, the trilinear coupling A t are always predicted to be large in deflected AMSB-type models, which are welcome to give sizeable contributions to the SM-like Higgs mass. Light stops and large A t can also improve the naturalness measurements of the theory. On the other hand, the positive value of A t tends to decrease to zero and further to large negative values when it RGE evolves down from high input scale to EW scale [92]. So, the values of A t at the EW scale may not be large for a mildly large messenger scale M mess , making the A t contribution to the SM-like Higgs mass not important for some range of M mess . Therefore, the stop masses are always not light because the allowed values of λ are small, leading to small tree-level contributions to the SM-like Higgs mass. It can be seen from Figure 4 that the Higgs mass can be as high as 123.7 GeV (124.5 GeV) and the messenger scale M mess is constrained to be larger than 3 × 10 11 GeV ( 1.2 × 10 10 GeV) if the muon g − 2 anomaly is explained upon 2σ (3σ) level, respectively.
As a comparison, the Higgs mass is upper bounded to be 118 GeV (120 GeV) when the muon g − 2 anomaly is explained at 2σ (3σ) level in the CMSSM/mSUGRA, because light sleptons also indicate light stops (with an universal m 0 input at GUT scale), leading to small loop contributions to Higgs masses. So, our deflected AMSB realization of NMSSM is much better in solving the muon g − 2 anomaly than that of minimal gravity mediation realization of MSSM.
• Our numerical results indicate that the lightest neutralino DM is always wino-like, which can annihilate very efficiently and lead to the under abundance of DM unless the DM particle mass is heavier than 3 TeV. The NMSSM-specific singlino component is negligibly small, which, therefore, will not play an important role in DM annihilation processes. Our numerical results indicate that the DM particle is constrained to be lighter than 500 GeV. Therefore, additional DM components, such as the axino, are always needed to provide enough cosmic DM. We also check (see the figures in Figure 5) that the Spin-Independent (SI) and Spin-Dependent (SD) DM direct detection constraints, for example, the LUX [93], XENON1T [94,95], and PandaX-4T [96,97], can be satisfied for a large portion of survived points.

Conclusions
Realistic deflected AMSB-type models are always predictive and can easily lead to light slepton masses when the tachyonic slepton problems that bother them are solved. In addition, the predicted gaugino mass ratios at the EW scale in deflected AMSB are given by M 1 :M 2 :M 3 ≈ (6.6 − 2d):2 (1 − 2d):6 (−3 − 2d), which are different to the ratios M 1 :M 2 :M 3 ≈ 1:2:6 that appeared in mSUGRA/CMSSM and GMSB. Therefore, much lighter electroweakinoes can still be consistent with the LHC 2.2 TeV gluino lower mass bound. Therefore, the soft SUSY breaking spectrum from deflected AMSB type scenarios with light sleptons and light electroweakinoes are always favored to solve the muon and electron g − 2 anomalies. As the AMSB type scenarios always predict insufficient winolike DM, additional DM species, such as the NMSSM-specific singlino component, can be welcome to provide additional contributions to the DM relic density. Therefore, the embedding of NMSSM framework into realistic AMSB-type UV theory is fairly interesting and well-motivated.
In this paper, we propose to provide a joint explanation of electron and muon g − 2 anomalies in UV SUSY models in the framework of the anomaly mediation of SUSY breaking. We embed the General NMSSM into the deflected AMSB mechanism with Yukawa/gauge deflection contributions and obtain the relevant soft SUSY breaking spectrum for General NMSSM. After integrating out the heavy messenger fields, the analytical expressions of the relevant soft SUSY breaking spectrum for General NMSSM at the messenger scale can be calculated, which can be RGE evolved to EW scale with GNMSSM RGE equations. Our numerical scan indicates that successful EWSB and realistic lowenergy NMSSM spectrum can be obtained in some parameter regions. Furthermore, we find that, adopting the positively central value electron g − 2 experimental data, it is possible to jointly explain the muon g − 2 anomaly and the electron g − 2 anomaly within a range of 1σ and 2σ, respectively. The Z 3 invariant NMSSM, which corresponds to ξ F = 0 in our case, can also jointly explain the muon and electron anomaly to 1σ and 2σ range, respectively.