NMSSM from alternative deflection in generalized deflected anomaly mediated SUSY breaking

We propose a new approach to generate messenger–matter interactions in deflected anomaly mediated SUSY breaking mechanism from typical holomorphic messenger–matter mixing terms in the Kahler potential. This approach is a unique feature of AMSB and has no analog in GMSB-type scenarios. New coupling strengths from the scaling of the (already known) Yukawa couplings always appear in this approach. With messenger–matter interactions in deflected AMSB, we can generate a realistic soft SUSY breaking spectrum for next-to-minimal supersymmetric standard model (NMSSM). Successful electroweak symmetry breaking conditions, which is not easy to satisfy in NMSSM for ordinary AMSB-type scenario, can be satisfied in a large portion of parameter space in our scenarios. We study the relevant phenomenology for scenarios with (Bino-like) neutralino and axino LSP, respectively. In the case of axino LSP, the SUSY contributions to Δaμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta a_\mu $$\end{document} can possibly account for the muon g-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g-2$$\end{document} discrepancy. The corresponding gluino masses, which are found to below 2.2 TeV, could be tested soon at LHC.


Introduction
Weak scale supersymmetry (SUSY) is one of the most interesting candidates for new physics beyond the standard model a e-mail: feiwang@zzu.edu.cn (SM). Its low energy phenomenology is determined mainly by the relevant soft SUSY breaking parameters which are required to preserve flavor and CP with a good accuracy. Such soft SUSY breaking parameters can be predicted by the mediation mechanism of SUSY breaking. So it is crucial to understand various well-motivated mediation mechanisms, for example, the gauge mediated SUSY breaking (GMSB) [1][2][3][4][5][6][7], the anomaly mediated SUSY breaking (AMSB) [8,9] mechanisms.
Minimal AMSB, which is determined solely by the parameter F φ m 3/2 , is insensitive to the UV theory [10] and predicts a flavor conservation soft SUSY breaking spectrum. Unfortunately, negative slepton squared masses will appear and the minimal scenario must be extended. Although there are many possible ways to tackle such tachyonic slepton problem, the most elegant solution from aesthetical point of view is the deflected AMSB [11][12][13][14] (dAMSB) scenario. In deflected AMSB, additional messenger sectors are introduced to deflect the AMSB trajectory and additional gauge mediation contributions can possibly push the negative slepton squared masses to positive values [15,16]. On the other hand, N ≥ 4 species are always needed to give positive slepton squared masses with naturally negative deflection parameters, possibly leading to strong gauge couplings below GUT scale or Landau pole below Planck scale. Besides, electroweak naturalness and the discovered 125 GeV Higgs boson by both the ATLAS [17] and CMS collaborations [18] of LHC, may indicate respectively the lightness of stop and large trilinear coupling A t ; null search results of sparticles with 36 f b −1 of data at the (13 TeV) LHC by the ATLAS and CMS collaborations [19][20][21][22] suggest that the low energy SUSY spectrum should display an intricate pattern, for example, the first two generation squarks need to be heavy to avoid the stringent constraints from LHC. However, spectrum of such a type will in general not be predicted in ordinary (d)AMSB scenarios.
We had proposed to introduce general messenger-matter interactions in dAMSB to solve the previous problems [23] which can be advantageous in various aspects. The previous formalism in [23,24], within which the messenger-matter interactions are introduced in the superpotential, has an analog in GMSB [25][26][27][28][29][30][31]. In this paper, we propose an alternative approach to including messenger-matter interactions in (d)AMSB from typical holomorphic terms in the Kahler potential. Such an approach, with the messenger scales set by the compensator vacuum expectation value (VEV), can only be possible in AMSB type scenarios.
Next-to-minimal supersymmetric standard model (NMSSM) [32,33] is a singlet extension of MSSM that has various advantages. The SUSY preserving μ parameter, which (however) needs to lie near the soft SUSY breaking scale to trigger EWSB, could be generated with the correct scale once the singlet scalar acquires a VEV. Besides, with additional tree-level contributions or through doubletsinglet mixing, NMSSM can easily accommodate the discovered 125 GeV Higgs boson without much fine-tuning. 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 [34,35]. Such difficulty always persists in ordinary (d)AMSB-type scenarios. We find that phenomenological interesting NMSSM spectrum can be successfully generated from non-trivial holomorphic messenger-matter mixing terms in the Kahler potential. Besides, the discrepancy between the theoretical predictions for the muon anomalous magnetic momentum and the experiments, can possibly be solved in our scenario. This paper is organized as follows. In Sect. 2, we propose our new scenario and discuss the general methods to obtain the soft SUSY parameters. The soft SUSY parameters for NMSSM are also given. The relevant numerical results are discussed in Sect. 3. Section 4 contains our conclusions.

Messenger-matter interactions from Kahler potential
There are two possible ways to deflect the AMSB trajectory with the presence of messengers, by pseudo-moduli field or holomorphic terms (for messengers) in the Kahler potential. In our previous work [23], messenger-matter interactions are introduced in the superpotential involving the pseudomoduli field. We find that messenger-matter interactions can also be consistently generated by holomorphic messengermatter mixing terms in the Kahler potential. In order to show the most general features of this approach, we discuss the relevant soft parameters in the framework of NMSSM.
As mentioned previously, the low energy soft SUSY breaking spectrum of NMSSM obtained from typical SUSY breaking mechanism, such as GMSB and ordinary extended (d)AMSB, is always bothered by the requirement to achieve successful electroweak symmetry breaking (EWSB), rendering the model building non-trivial [34,35]. In fact, EWSB within NMSSM in general requires a large VEV for the singlet. This prefers a negative m 2 S and/or large A-terms for the singlet superpotential interactions, A λ and A κ [36]. AMSB always predicts a positive value for m 2 S and A λ , A κ can be sizable only if λ, κ are also large which would induce a larger m 2 S , suppressing the singlet VEV. Such difficulty can possibly be ameliorated in AMSB-type scenarios with enhanced trilinear couplings. We find that phenomenological interesting NMSSM spectrum can be successfully generated by introducing holomorphic messenger-matter mixing terms in the Kahler potential.
We propose to introduce new holomorphic term in the Kahler potential in addition to ordinary canonical Kahler kinetic terms and Z 3 invariant NMSSM superpotential Here φ is the conformal compensator field and S 1 , S 2 , T are gauge singlet superfields;P m,a and Q m,a are the standard model matter superfields in the5 and 10 representations of SU (5) with 'a = 1, 2, 3' the family index. Additional vectorlike messengers in theP, P 5 ⊕ 5) andQ, Q (10 ⊕ 10) representations of SU (5) are introduced to solve the tachyonic slepton problems and at the same time deflect the ordinary AMSB trajectory. In this paper, we assume c P , c Q , c P,Q m,a are real and c P,Q m,a = 0 only for a = 3 for simply. Such mixing between the third generation fermions and additional vectorlike fermions always appear in new physics models, such as top (top-bottom ) seesaw model or extra dimension models. We should note that c P,Q m,a = 0 for the first two generations a = 1, 2 are also possible and such possibility will be discussed subsequently. Besides, the choice of superpotential W M is not unique. For example, we can adopt couplings between S and messengers used in [37] or [38] (with double messenger species).
With only c P,Q m,3 = 0, the holomorphic terms in the Kahler potential reduce to after rescaling each superfield with the compensator field φ, namely → φ . With the F-term VEVs of compensator φ = 1 + F φ θ 2 , we have the potential for the singlets where the first term are obtained by picking out the F φ terms.
We thus arrive at the mass matrix for the scalar components It can be seen that the mass matrix has vanishing determinant.
We can define c S ≡ c 2 S,1 + c 2 S,2 and redefine the fields The mixing angle can be given as (2.6) A minus sign for the angle is kept for future convenience. The zero eigenvalue of the scalar mass matrix corresponds to the combinationS 1 . The fermionic component ofS 1 , which is orthogonal toS 0 , can also be seen to be massless from the Kahler potential. The non-vanishing mass eigenstates for scalar matrix are given by (2.7) Such expressions are analog to that of the GMSB with T,S 0 the messenger-like fields. After integrated out the messengers, we can obtain their contributions to the low energy soft SUSY breaking spectrum. S 1 can be written as the combination which, after substituting into the superpotential, can lead to couplings between the massless fieldsS 1 and heavy massive messenger-type fieldsS 0 . Similar to the gauge singlet case, we can define the massive combinationsK 1 , K 1 for P, Q within Eq. (2.2) and their massless orthogonal combinationsK 2 , K 2 as 3 ) can be written as the combination of massive stateK 1 (K 1 ) and massless stateK 2 (K 2 ) P m,3 = sin ψ 1K1 + cos ψ 1K2 , The SUSY breaking effects from compensator F-term VEVs can be taken into account by introducing a spurion superfields R W = d 2 θ c P RP P + c Q RQ Q + · · · , (2.12) with the spurion VEV as which gives the deflection parameter (2.14) The spurion messenger-matter interactions will affect the AMSB RGE trajectory after integrating out the heavy modes.
The superpotential involving the matter and singlet superfields at the F φ scale are given as The terms containing the 10,5 representations of SU (5) reduce to the sum of their components below the GUT scale and should be understood as the abbreviation of this sum at low energy a=1,2,3 The same holds for terms containing K a ,K a . Besides, terms involving the triplet components of H,H are integrated out by assuming proper doublet-triplet splitting mechanism to generate heavy triplet Higgs masses.
After integrated out all the heavy modes including T,S 0 ;K 1 , K 1 ; P,Q etc, the low energy theory will reduce to NMSSM. The effects of integrating out the messengers can be taken into account by using Giudice-Rattazi's wavefunction renormalization approach [39]. The messenger threshold M R can be further promoted to the other chiral spurion field X with M R = √ X † X . The superfieldS 1 will act as the singlet S appearing in ordinary NMSSM superpotential and K 2 , K 2 as the third generation superfields. Note that there is a scaling of various couplings appeared in NMSSM λ =λ sin θ, κ =κ sin 3 θ, y b,τ = y D,L 33 cos ψ 1 cos ψ 2 , In the subsequent studies, we must ensure that correct scaled (or unscaled) couplings are used in the expressions. New flavor dependent interactions involving the messengers from y D,L i jP m,i Q m, jH5 and y U i j Q m,i Q m, j H 5 are not dangerous because these new flavor dependent interactions are aligned with the MSSM Yukawa coupling. Diagonalizing the MSSM Yukawa couplings will simultaneously diagonalize these additional Yukawa couplings.
We should briefly discuss the most general case with c P m;1,2 = 0. Similar to Eq. (2.9), we can define the mixing matrix within which the three orthogonal combinationsK a(a = 1, 2, 3) are determined up to an arbitrary rotation under SO(3) transformation O P . The scalar components of the three combinationsK a(a = 1, 2, 3) correspond to the massless eigenvalues of the 5 × 5 sfermion mass matrix for (P * ,P,P m 1 ,P m 2 ,P m 3 ) that can be identified to be the three generation squark/sleptons in MSSM. Similar conclusions hold for Q m a . Besides, the relations between y t , y b , y τ and y U 33 , y D 33 , y L 33 will be non-trivial and depends on various parameters in the mixing matrix After fixing the rotation matrix O P and O Q , we can obtain precisely the relations between SM Yukawa couplings and the couplings in the superpotential. For new flavor dependent interactions involving one messengers, such as (y U ) 0bK Q;0KU c L ;b H u , the coupling can be seen to satisfy   [40,41] type mechanism with additional horizontal flavor symmetry.
After integrating out the messenger fields, the wavefunction will depend on the messenger threshold set by the spurion superfield R. The soft gaugino masses are given at the messenger scale by Because of the non-renormalization of the superpotential, the trilinear soft terms will be determined by the wavefunction normalization as The anomalous dimension are expressed in the holomorphic basis [42] as (2.25) with G ≡ G + −G − the discontinuity across the messenger threshold. Here G + (G − ) denotes the value above (below) the messenger threshold, respectively. The soft SUSY breaking scalar masses are given by The dAMSB soft scalar masses can be divided into several parts, namely the gauge-anomaly interference part, the pure gauge mediation part as well as the ordinary anomaly mediation part.

Numerical results
From the previous general formulas for soft SUSY breaking parameters, we can obtain the analytical expressions in our scenario at the scale F φ after integrating out the messenger fields. Some of the lengthy expressions are given explicitly in the appendix. We can see that new contributions will be given to the trilinear couplings A κ .A λ , m 2 S which could be helpful to trigger EWSB in NMSSM.
We use NMSSMTools5.1.2 [43,44] to scan the whole parameter space. The free parameters of our inputs are chosen to satisfy with the range of the mixing angle −π/2 < θ, ψ 1 , ψ 2 < π/2 and we require that λ 2 + κ 2 0.7 to satisfy the perturbative bounds. The soft SUSY mass Z by the minimization conditions of the scalar potential. Usually, the parameter M A can be used to replace A κ by We note that κ is a free parameter while the tan β ≡ v u /v d is not. This choice is different to ordinary numerical setting from top-down in which tan β is free while κ is a derived quantity [45]. Such choice is convenient for those predictable NMSSM models with a UV completion. Firstly, we need to guess a value for tan β to obtain the relevant Yukawa y t , y b couplings at the electroweak (EW) scale. After renormalization group equation (RGE) evolution of g i , y t , y b couplings from EW scale 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 the minimization condition of the Higgs potential [46]. The purpose of deflection in AMSB is to solve the notorious tachyonic slepton problem. So non-tachyonic slepton should be obtained at the EW scale.
In our scan, we also impose the following collider constraints  A possible hint of new physics beyond the SM is the muon g − 2 anomaly. The E821 experimental result of the muon anomalous magnetic moment at the Brookhaven AGS [55] is given by which is larger than the SM prediction [56,57] a SM μ = 116591834(49) × 10 −11 . (3.5) We adopt the conservative estimation 4.7 × 10 −10 a μ 52.7 × 10 −10 in our scenario. If possible, we would like our theory to explain such muon g-2 anomaly.
To address the strong CP problem in SUSY, the axino, which is the fermionic superpartner of axion, can be predicted. If the existence of axino could be confirmed by experiments, the two theoretical hypotheses, which are designed to solve two respective hierarchy problems: the strong CP problem by a very light axion and the gauge hierarchy problem by SUSY, could be validated.
The SUSY version of the Kim-Shifman-Vainshtein-Zakharov (KSVZ) [58,59] axion model, which introduces a SM singlet and a pair of extra vector-like quarks that carry U (1) P Q charges while keeps the SM fermions and Higgs fields neutral under U (1) P Q symmetry, could be introduced in NMSSM with additional singlets and messengers. The decoupling theorem [60] of anomaly mediation, which states that pure mass threshold will not deflect the AMSB trajectory, guarantee that such new messengers, with their masses determined by the singlet VEVs between 10 10 GeV S 10 12 GeV (to be compatible with the experimental bound on f a [61]) will not change the original AMSB predictions upon F φ .
The axino mass, which is strongly model dependent, can be much smaller or much larger than the SUSY scale. Some models predict mã of the order of gravitino mass m 3/2 [61] while the spontaneously broken global SUSY predict it to lie of order mã ∼ O(M 2 SU SY / f a ) [62]. So one usually treats axino mass as a free parameter and assumes the axino interactions to be given by the U (1) P Q symmetry. If axino is heavy, the lightest neutralino will be DM candidate; if the axino is the LSP and acts as DM particle, colored or charged NLSP can decay into axino and non-thermally generate the observed axino DM relic density.
Depending on the nature of the LSP, we have the following discussions from our numerical results • Scenario I-The lightest neutralino χ 0 1 as the LSP: Successful EWSB conditions impose stringent constrains on the input parameters of NMSSM, especially when such low energy inputs are determined by a UV-completed theory. Random scan in this case indicates that still some points can survive the EWSB conditions by leading to a iteratively stable value of tan β. The allowed range for characteristic NMSSM parameters κ and λ can be seen in the upper left panel of Fig. 1. Besides, the discovered 125 GeV Higgs can be successfully interpreted in our scenario. The mass scales of the soft spectrum are uniquely determined by the gravitino mass F φ in AMSB-type scenarios. The plot of F φ versus Higgs boson mass is shown in the upper right panel of  Fig. 1. In the NMSSM, the SM-like Higgs can be pushed to 125 GeV by additional tree-level contributions for relatively large λ and small tan β. If the lightest CP-even scalar is mostly singlet-like, mixing between it and the SM-like Higgs can lead to an increase of the SM-like Higgs mass. Although the possibility of SM-like Higgs being the second lightest CP-even scalar is attractive, we find that almost all the survived points predict the lightest CP-even scalar as the 125 GeV Higgs. Numerical results indicate that the muon g − 2 anomaly cannot be explained in this scenario. The SUSY contributions to a μ are shown in different color in the upper right panel of Fig. 1. This result can be understood because the required a μ can be achieved only if the relevant sparticles (μ,ν μ ,B,W ,H ) are lighter than 600 ∼ 700 GeV for tan β ∼ 10 in MSSM [63][64][65][66][67][68]. The inclusion of singlino in NMSSM will not give sizable contributions to a μ because the coupling of singlino to MSSM sector is suppressed. Although the two loop contributions involving the Higgs is negligible in SM, the new higgs bosons in NMSSM could have an important impact on a μ if the lightest neutral CP-odd Higgs scalar is very light [69]. As noted there, 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. In our scenario, the lightest CP-odd Higgs a 1 is not light enough to give sizable contributions to a μ . So the main contributions are similar to that in the MSSM and are not large enough to account for the discrepancy. In the case of χ 0 1 as the LSP, the dark matter (DM) relic density is required to satisfy the bounds set by the Planck data [70] in combination with the 5σ WMAP data It is well known that the Bino LSP will lead to overabundance of DM unless co-annihilation with sfermions or s-channel resonance increase the annihilation rates. In our scenario, the Bino-like DM always co-annihilates with the lightest stau or sbottom and leads to the right DM relic density. The mass differences betweenτ 1 (b 1 ) NLSP andχ 0 1 LSP are typically of order O(GeV ). Interactions between Bino dark matter and the nucleons are primarily mediated by t-channel scalar Higgses or by s-channel squarks. The Higgs exchange diagrams dominate the Spin-Independent (SI) cross section of χ 0 1 on nucleon for sufficiently heavy squarks. The Spin-Dependent (SD) cross section, which is dominated by the Z-boson exchange diagram as long as the squarks are sufficiently heavy, can provide a good probe of the gaugino and Higgsino parameters. Large SD cross sections can appear only if χ 0 1 DM contains a large Higgsino component. Although the SD direct detection cross section is generally larger than the SI cross section, it is much more difficult to probe in the experiment, as the SD cross section does not scale directly with the mass of the nuclei. In fact, current bounds on the neutron SD cross section are less stringent by a factor of 10 6 . The SI and SD direct detection results are shown in the lower right panel of Fig. 1. In our scenario, the neutralino DM is almost pure Bino with very tiny Higgsino and Wino contents. Such (almost pure) Bino-like DM can survive the direct detection constraints from LUX [72] and PANDAX [73] because of the suppressed Higgsino component.
• Scenario II-The axino as the LSP: As mentioned before, the axino can be the LSP and act as the DM candidate. So the previously forbidden charged or colored LSP region can be revived. In some of the allowed parameter space, numerical scan indicate that lightest ordinary supersymmetric particle (LOSP) will be the lightest sbottom, stau, tau-sneutrino etc. Even in this scenario, successful EWSB conditions can impose stringent constrains on the input parameters of NMSSM, especially for characteristic NMSSM parameters κ and λ. Their allowed range can be seen in the upper left panel of Fig. 2. The plot of F φ versus Higgs boson mass is shown in the upper right panel of Fig. 2. We can see that the 125 GeV Higgs mass can also be easily accommodated in this scenario. The muon g − 2 discrepancy, however, prefer this scenario. The most stringent constraint, as noted in our previous papers [74], is the gluino mass bound. We can see that the SUSY contributions to a μ can reach 14×10 −10 in the case of slepton (τ 1 ,ẽ R ) LOSP with mg ≈ 1.8 TeV.
The mass bound of gluino can be understood as follows: as the whole low energy spectrum is determined by the value of F φ , the mass scales ofμ etc determine the upper bound of F φ , which, on the other hand, sets a bound on gluino mass. We find that the gluino mass is bounded to less than 2.2 TeV if we adopt the lower limit of the required a μ . Such light gluino will soon be tested by the future LHC searching results. The explanation of the muon g−2 discrepancy prefers the effective μ parameter to be light, which could reduce the fine tuning (FT) involved in our theory. Besides, the stop mass can also be predicted to lie just upon the experimental bounds from recent LHC and be fairly natural [75]. We show the Barberi-Giudice (BG) FTs [76] for the survived points in the figures. We can see that the FT involved in this scenario can be as low as 50. DM particle in this scenario is the axino from LOSP decay (neglecting the effects of saxion condensate [77]). As the non-thermally produced relic density of axino can be related to that of NLSP via the relic density of DM will not impose stringent constraints on our parameter space. We show the relic density of NLSP in Fig. 2. The thermal production of axino, which may give important contributions (or even dominate over the non-thermally production) to the observed relic density and be model dependent, will impose strong bound on the reheating temperature after inflation by requiring such production should not overclose the universe. 1 Relevant discussions on the reheating temperature and cosmological consequences of axino DM can be found in [79,80]. The extremely weak interaction strength of axino makes its detection in DM direct detection experiments, as well as at collider experiments, rather hopeless. The hint of the axino DM may show up from the properties of the LOSP. In the case of axino LSP, the LOSP typically decays with a lifetime of less than one second and practically be stable inside the collider detector. The (electrically) charged or colored particle as a LOSP would appear as a stable particle inside the detector. The injection of high-energetic hadronic and electromagnetic particles, produced from late decays of an NLSP to axino (with lifetime less than one second), will not affect the abundance of light elements produced during Big Bang Nucleosynthesis (BBN).

Conclusions
We propose a new approach to generate messenger-matter interactions in deflected anomaly mediated SUSY breaking mechanism from typical holomorphic messenger-matter mixing terms in the Kahler potential. This approach is a unique feature of AMSB and has no analog in GMSBtype scenarios. New coupling strengths obtained from the scaling of the (already known) Yukawa couplings, such as y t tan ψ 2 , y t tan 2 ψ 2 , y b tan ψ 2 tan ψ 1 , κ cot θ, κ cot 2 θ , always appear in this approach, which is a salient feature of this scenario. With messenger-matter interactions in deflected AMSB, we can generate a realistic soft SUSY breaking spectrum for next-to-minimal supersymmetric standard model (NMSSM). Successful electroweak symmetry breaking conditions, which is not easy to satisfy in NMSSM for ordinary AMSB-type scenario, can be satisfied in a large portion of parameter space in our scenarios. We study the relevant phenomenology for scenarios with (Bino-like) neu-tralino and axino LSP, respectively. In the case of Bino-like LSP, most of parameter space can survive the recent SI and SD DM direct detection experiment; in the case of axino LSP, the SUSY contributions to a μ can possibly account for the muon g-2 discrepancy. The corresponding gluino masses (for those points that can explain the muon g-2 anomaly), which are found to below 2.2 TeV, could be tested soon at LHC.

Appendix: The soft SUSY breaking parameters
We list the relevant soft SUSY breaking parameters obtained from the analytical formulas (2.22) to (2.26). For gaugino masses, we have The trilinear couplings are given as with the deflection parameter d = −2.