Probing Non-holomorphic MSSM via precision constraints, dark matter and LHC data

In this analysis we explore the phenomenological constraints of models with non-holomorphic soft SUSY breaking terms in a beyond the MSSM scenario having identical particle content. The model referred as NHSSM shows various promising features like the possibility of a strong reduction in electroweak fine-tuning even for a scenario of a heavy higgsino type of LSP, a fact that is unavailable in pMSSM models. The other important aspect is satisfying the muon $g-2$ data even for a small $\tan\beta$ via a small value of coupling $A_\mu'$ associated with the tri-linear non-holomorphic soft term. Thus, a large SUSY contribution to muon $g-2$ is possible even for a significantly large smuon mass $m_{\tilde {\mu_1}}$. The Higgs mass radiative corrections are contributed by both the holomorphic and non-holomorphic trilinear soft parameters $A_t$ and $A_t'$, thus diluting the requirement to have a larger $A_t$ to satisfy the Higgs mass data. The model also provides with valid parameter space satisfying the constraint of $B \rightarrow X_s +\gamma$ for large values of $\tan\beta$, a scenario unfavourable in pMSSM.


Introduction
The discovery of the Higgs Boson at the ATLAS [1] and the CMS [2] experiments of the Large Hadron Collider (LHC) marks the completion of particle searches within the realm of the Standard Model (SM) [3]. The SM is quite successful in explaining electroweak and strong interactions and the associated Higgs mechanism is found to be a viable method for generating masses for fermions and electroweak gauge bosons. Despite its success in explaining most of the observed experimental results, there are many theoretical issues and experimental facts that cannot be addressed while staying within the SM. The gauge hierarchy problem, baryogenesis, the fact that neutrinos have masses, the absence of a dark matter candidate, are a few of the important issues that motivate us to explore Beyond the SM (BSM) scenarios. Models involving Supersymmetry (SUSY) such as the Minimal Supersymmetric Standard Model (MSSM) [4][5][6][7] are prominent candidates for BSM physics.
However, the fact remains that even after the first few years of running of the Large Hadron Collider (LHC), SUSY is yet to be found. This has obviously put serious constraints on various models of low energy SUSY. In the post-Higgs discovery years, the lighter Higgs boson of MSSM to have a mass of m h ∼ 125 GeV [8] translates into large radiative corrections [9]. This demands a heavier top squark sector. A large fine-tuning is to be accepted. Furthermore, LHC has pushed up the lower limits of masses of the first two generations of squarks as well as gluino beyond a TeV. At the same time, SUSY models are increasingly being constrained via B-physics related measurements at LHCb. On the other hand, regarding dark matter (DM) [10], the measurements from WMAP/PLANCK [11,12] for the DM relic density or LUX [13] experiment for DM direct detection have put significant limits.
We further emphasize that, the data from the Brookhaven experiment for the anomalous magnetic moment of muon or (g − 2) µ points out a significant deviation (3.2σ) from its SM based evaluation while we note that various uncertainties of the SM contributions to (g − 2) µ are being reduced over the last few years. This leads to stringent constraints on the scalar and gaugino sectors of SUSY models. A combined requirement for satisfying the relic density range from WMAP or PLANCK experiments apart from satisfying the LHC derived sparticle mass [14] bounds particularly creates tension so as to have a reasonably large a SUSY µ . Here, a SUSY µ refers to the SUSY contribution to the theoretical evaluation of the muon anomaly a µ ≡ 1 2 (g − 2) µ which is supposed to be equal to the difference between the experimental value and the SM evaluation of the observable. The Higgsinos and Wino are typically required to be heavy for becoming candidates of dark matter. On the other hand, satisfying (g − 2) µ with dominant contributions from the chargino-sneutrino loops demands non-decoupling higgsinos or wino along with a light sneutrinoν µ . We emphasize that in models with gaugino mass universality like minimal supergravity (mSUGRA) [15]/constrained MSSM (CMSSM) [5], the chargino-sneutrino loops dominate in a SUSY µ [16,17]. Thus in the present scenario of sparticle mass limits a SUSY µ is not large enough to explain the observed deviation. In the parameter space of MSSM that is consistent with the Higgs mass data, dark matter relic density, collider limits for scalar and gaugino (electroweakino) masses, it turns out that the primary contribution to a SUSY µ comes from the loop diagrams containing neutralinos and smuons, particularly from the L-R mixing terms that scale with M 1 µ tan β [18]. In order to have a larger a SUSY µ so as to account for the deviation the above quantity needs to be large or in other words this restricts the smuon mass mμ 1 to become large. A comprehensive analysis in a model that identifies the valid region of parameter space satisfying limits from flavor physics such as that from Br(B → X s + γ) and Br(B s → µ + µ − ) (which we would collectively refer as B-physics constraints), dark matter constraints, while also having a moderate degree of fine-tuning and most importantly that would easily accommodate the (g − 2) µ limits even for a small tan β all at one go, is undoubtedly important. In this analysis, keeping ourselves contained within the MSSM particle setup we would like to explore whether a consideration of non-holomorphic (NH) soft SUSY breaking terms may be able to reduce the stringency arising out of the (g − 2) µ constraint in particular apart from satisfying all the above mentioned phenomenological requirements.
Away from MSSM, particularly in models with singlet scalars NH soft breaking terms potentially fall in the class of terms that may cause hard SUSY breaking [19][20][21][22][23][24]. Considering, for example, a hidden sector SUSY breaking scenario like supergravity, one may generically consider a spontaneous SUSY breaking due by the vacuum expectation value of an auxiliary field F belonging to a chiral superfield X. This causes appearance of soft terms in the Lagrangian that are associated with the coupling of X with another chiral superfield Φ or a gauge field strength superfield W a α . In a supergravity framework where the mass scale M is large (typically the Planck mass) one obtains the following [20].
Here, the parameter µ is introduced for the soft term so as to follow closely with the usual MSSM notation [5]. Considering the vacuum expectation value X = θθ < F > and denoting < F > simply by F one has the usual soft terms of MSSM namely the gaugino mass term, the cubic and the analytic scalar squared mass terms and a non-analytic scalar mass term coming out of the above D-term contribution as given below.
We note that F/M should refer to a weak scale mass which we consider here as the W-boson mass M W . Apart from the above contributions, there can be other D-term contributions that are classified as "maybe soft" [20] such as the ones arising out of 1 M . These terms in a broader sense may cause quadratic divergence, thus become hard SUSY breaking terms, in scenarios where the visible sector contains a singlet superfield. Nevertheless, the terms are generally highly suppressed in a supergravity type of scenario. Indeed this is why such NH contributions to the Lagrangian are traditionally ignored while discussing models with high scale based SUSY breaking. We must, however, note that, as pointed out in [25], in absence of any gauge singlet field in the visible sector such suppression may not be possible if the supersymmetry breaking effect is communicated to the visible sector at a lower energy. For example, this may happen in scenarios with gauge-mediated supersymmetry breaking [25] 1 . Thus having no gauge singlet superfields MSSM may as well include such NH terms that are soft SUSY breaking in nature. The terms can hardly be ignored in the most general sense. This was discussed or at least pointed out in several works [20-22, 24, 27-31].
Apart from the references shown above related to the issue of absence of quadratic divergence in MSSM in presence of NH soft terms and the possible origin of the terms in relation to a hidden sector SUSY breaking model, we will now briefly refer to a few specific works related to phenomenology.
A general analysis with NH soft terms that incorporated renormalization group evolutions with or without using R-parity violation was presented in Ref. [29]. This was followed by a study [30] with NH SUSY breaking terms in an essentially Constrained Minimal Supersymmetric Model (CMSSM) 1 For scenarios with F -term SUSY breaking leading to 1 M 3 type of suppression in nonstandard supersymmetry breaking terms including Dirac gaugino mass terms see ref. [26]. setup with all input parameters including the NH ones being given at the gauge coupling unification scale while exploring relevant phenomenological constraints like Br(B → X s + γ). Analyses that particularly focussed on the Higgs sector and constraints like Br(B → X s + γ) are Refs. [31,32].
However, in contrast to Ref. [31] that used input parameters given entirely at the unification scale, Ref. [32] analyzed in a mixed set-up where the NH parameters were given at the electroweak scale in an otherwise CMSSM type of setup. Similar mixed input parameters were used in Ref. [28] that discussed reparametrization invariance in special circumstances of the choice of parameters 2 , spectra as well as the effect on fine-tuning. Ref. [33] may be seen for its emphasis on the Higgs sector while considering CP-violating phases. Nonstandard SUSY breaking was also used in phenomenological studies with R-parity violating NH soft SUSY breaking terms in MSSM framework in Refs. [29,34,35].
Unlike all the previous analyses where universal models of CMSSM type were considered with NH parameters being given either at the unification scale or at the electroweak scale, our work on the Non-Holomorphic Supersymmetric Standard Model (NHSSM) will entirely use electroweak scale input parameters similar to what is considered in phenomenological MSSM (pMSSM) model [36] in relation to MSSM. In this phenomenological NHSSM (pNHSSM) framework, we will explore the extent NH parameters influence on satisfying the (g − 2) µ constraint apart from the effect on electroweak fine-tuning and m h via the associated radiative corrections. The Higgs mass limit is achieved for smaller values of |A t |, the trilinear coupling parameter corresponding to the top-quark in comparison with what is required for MSSM. Furthermore, as we will see such NH parameters do not affect the fine-tuning measure since the Higgs scalar potential would not have any dependence on such parameters 3 .
We will additionally focus on the low energy processes like Br(B → X s +γ) and Br(B s → µ + µ − ) that also receive a significant amount of contributions from the NH terms. Particularly these constraints may indeed be quite severe for large values of tan β in MSSM. However, the NH terms are able to alter the above branching ratios and a large region of parameter space that would be excluded in MSSM is restored. We will further see that there are valid regions of pNHSSM parameter space that is consistent with correct relic abundance for dark matter (DM) and these may be probed for the Direct and Indirect Detection limits of DM.
Our paper is organized as follows. In Sec.2 we discuss NHSSM and particularly explain the impact of NH terms on the i) Higgs and other scalar sectors, ii) charginos and neutralinos, electroweak fine-tuning and iii) phenomenological aspects related to the constraints coming from dark matter, muon g − 2 as well as Br(B → X s + γ) and Br(B s → µ + µ − ). We present the results of our analysis in Sec.3. Finally, we conclude in Sec.4.

Non-Holomorphic Supersymmetric Standard Model
We remind that MSSM is considered to have only holomorphic soft SUSY breaking terms. The trilinear soft terms, in particular, are given by as follows [5] We have only shown here the dominant terms involving the third generations of fermions. It was shown that in the absence of any Standard Model gauge singlet it is possible to extend the SUSY breaking soft sector by including NH soft SUSY breaking terms, without aggravating any quadratic divergence [19,22]. Thus the NH soft terms of the NHSSM in general that include trilinear coupling terms as well as a coupling term involving Higgsinos are given by [31,32], We will now discuss the effect of involving NH terms on several sectors of MSSM sparticle spectra, particularly in the Higgs sector, the squark and slepton sectors, the electroweakinos (charginos and neutralinos) apart from its effect on a fine-tuning measure. We will also include low energy and cosmological observables like dark matter relic density while also taking into account the LHC bounds for sparticles and the Higgs mass data [8].

Influence of non-holomorphic terms on the scalar sector
The NH trilinear coupling parameters may cause a significant amount of change in the masses of squarks and sleptons. For example, the mass matrix for the up type of scalar quark, in general, is given by [31,32].
Similar matrices for a slepton or a down-type of squark may be written as given below.
Clearly, µ of MSSM that contributes to L-R mixing of squark is replaced by µ + A u in NHSSM.
The contributions of the NH terms will thus be more effective for i) low tan β in the case of up type of squarks and ii) large tan β in case of down type of squarks or sleptons.
An effect on top-squark sector is transmitted to the Higgs mass radiative corrections. We remind that in the MSSM framework the discovery of the Higgs boson with a mass of 125.09 ± 0.24GeV [8] is translated into a large radiative corrections to the mass of the lighter neutral CP-even Higgs boson h. 4 . The above requirement pushes up the masses of the top-squarks in MSSM which in turn indicates the need for a large value of |A t | so as to have a larger Left-Right mixing.
Considering the top-stop loops which constitute the most contributing terms, the above radiative corrections in NHSSM read as follows.
Herem t refers to the running top-quark mass that includes corrections from the electroweak, QCD and SUSY QCD effects. The maximal mixing scenario refers to X t = √ 6M S [5,6] where M S = √ mt 1 mt 2 . We note that unlike MSSM, with a suitable combination of the signs of A t and A t , it is possible to limit |A t | so as to satisfy the Higgs mass constraint.
We shall now discuss the effect of considering NHSSM on the Electroweakino (Chargino-Neutralino) sector and fine tuning.

Electroweakinos in NHSSM
The NH parameter µ modifies both the chargino and neutralino mass matrices. Essentially the charginos and neutralinos now have higgsino components corresponding to a higgsino mass of |µ − µ |. Thus in NHSSM the neutralino mass matrix reads [31,32]: Similarly for the chargino matrix we have [31,32]: We note that the LEP bound on the lighter chargino mass will essentially apply to |µ−µ | instead of µ. However, as we will see below the fine-tuning measure is still dependent on µ rather than µ and this provides with a quite unique signature of a possibility of having low fine-tuning irrespective of the nature of dark matter considered in NHSSM. For example, one can have a bino-like LSP even for very low µ.

Electroweak fine tuning in pNHSSM
The nonholomorphic trilinear parameters are associated with charged or colored scalars whereas the parameter µ is associated with fermions (higgsinos). Thus with no influence of the nonholomorphic 4 Precision measurements indicate that the discovered Higgs Boson is consistent with being SM-like [8]. With the above in perspective we consider the decoupling regime of MSSM Higgs characterized by M 2 soft breaking terms on the neutral scalar potential, the latter is same as that of MSSM as given below.
where µ is the bilinear Higgs mixing parameter of the superpotential, m Hu , m H d are scalar mass parameters and b is the Higgs mixing parameter within the SUSY breaking soft sector. After minimization with respect to the vacuum expectation values (vevs) of the neutral Higgs scalars, one finds, the well-known relations for electroweak symmetry breaking [5,6] and, (2.10) Clearly, at the tree-level one requires a fine cancellation between the two terms in the right-hand side of Eq.2.9 coming from SUSY breaking parameters m H d , m Hu and supersymmetry preserving parameter µ for obtaining the left-hand side namely m 2 Z /2 where m Z refers to the measured value of 91.2 GeV. The degree of cancellation broadly indicates a measure of the fine-tuning at the electroweak scale. We consider here a general definition of electroweak fine tuning (EWFT) that uses log derivatives [46][47][48], namely where p i ≡ {µ 2 , b, m Hu , m H d } are the parameters that determine the tree-level value of Z boson mass. The total EWFT at low scale is given by, It turns out that with large tan β the most important terms are ∆(µ) 4µ 2 Z tan 2 β . The above expressions show that even for a moderately large tan β, a small value of EWFT demands a lower value of µ. However, very small values of µ are excluded in MSSM due to LEP bound of lighter chargino mass and we will see that this has an important significance in relation to the fine tuning in NHSSM in Sec.3. We however point out that for small tan β and very small µ (typically much smaller than the above chargino mass limit) situation may arise where ∆(m Hu ) and ∆(m H d ) become larger than ∆(µ) [49]. As a result even for negligible values of µ one may obtain finite EWFT as we will see in Sec.3.3.
In our discussion on EWFT we must, however, remember that, we should include the principal corrections due to one-loop radiative effects due to top-stop loops and this leads to [49], Albeit this depends on the choice of the cut-off scale Λ. The requirement of a large A t in the post Higgs@ 125 GeV scenario increases δm 2 Hu which specially shows the need to include the one-loop corrections shown above. This obviously enhances the EWFT, although we will not include this effect in our EWFT measure. We must, however, point out that minimizing the Higgs potential near a scale where the logarithmic term vanish in Eq.2.13 may reduce the requirement of inclusion of the one-loop radiative corrections significantly as it was discussed in the context of hyperbolic branch and electroweak fine-tuning Ref. [50].
As noted before, a small value of µ is consistent with smaller EWFT. However, the lighter chargino mass bound from LEP limits µ hence to the EWFT measure not to become too small. In NHSSM the higgsino content of electroweakinos of Eqs.2.6 and 2.7 depends on |µ − µ |, thus EWFT may become small irrespective of the mass of electroweakinos.

Low Energy Constraints viz. (B
It is known that rare B-decays within the SM like (B → X s + γ), (B s → µ + µ − ) that are helicity suppressed may have large contributions from the radiative corrections due to superpartners in the loops. The SM contributions to (B → X s + γ) almost saturate the experimental data. Thus any BSM correction should be low enough to accommodate the difference of the SM and the experimental results, which is albeit small. SUSY parameter space is thus strongly constrained via cancellation of dominantly contributing diagrams, while the contributions individually may be large. One must also remember that there are next to leading order (NLO) contributions that can also be quite significant, specially for large values of tan β. In the SM the dominant radiative corrections come from t − W loops. In MSSM the significantly contributing diagrams involve t − H ± andt −χ ± loops. The contributions from the former loops share the same sign with the t − W loop contributions of the SM. In NHSSM, the soft terms from the NH trilinear coupling A t and the bilinear higgsino coupling µ have significant influences on the SUSY diagrams associated with the above flavor related processes. In contrast to MSSM, loops involving left-right mixing of top-squarks in NHSSM are associated with the factor A t − (µ + A t ) cot β [31,32], whereas for the contributions involving the higgsino loops, µ is replaced by the difference (µ − µ ) [31,32].
Regarding the constraint from Br(B s → µ + µ − ) that is typically stringent for large tan β and small pseudoscalar mass in MSSM models, it turns out that the available parameter space that survives after imposing the Br(B → X s +γ) constraint, is not affected much when one imposes the constraint from Br(B s → µ + µ − ).

Results
The focus of our analysis in this section would be the important features of NHSSM in relation with MSSM, particularly the effects of considering NHSSM on i) Higgs boson mass, ii) flavor violating processes like Br(B → X s + γ) and Br(B s → µ + µ − ) iii) SUSY contributions to muon g − 2, iv) dark matter relic density and v) electroweak fine-tuning. Regarding the last two points we will particularly demonstrate the fact that NHSSM allows to have a higgsino-like DM with a sufficiently low electroweak fine-tuning, a feature almost impossible to get in MSSM. We will additionally show a few benchmark points consistent with the various above constraints.

Impact of non-holomorphic soft parameters on m h
The effects of nonholomorphic parameters particularly A t on radiative corrections to the CP-even lighter Higgs boson mass m h as enumerated in Eq.2.5 is particularly prominent for smaller tan β.
Keeping this in mind, we choose tan β = 10 and show the extent of variation of m h due to varying A t and A t . This is displayed in Fig.1(a). We assume a 3 GeV window in m h leading to the following range [51].    Fig.1(a) shows the variation of m h against A t for the scanning ranges of Eq.3.2 for tan β = 10.
The magenta and cyan colored regions correspond to NHSSM and MSSM respectively. Fig.1(b) is same as In order to probe NHSSM signatures on the Higgs boson mass we plot both the NHSSM and MSSM specific parameter points in Fig.1(a). The magenta colored points in Fig.1(a) correspond to the NHSSM scenario where variation due to relevant holomorphic and non-holomorphic parameters are as referred in Eq.3.2. We isolate the MSSM specific parameter points in the cyan colored region by a choice of A t = µ = 0. Clearly, focusing on the non-maximal region of m h and a given value of A t we note that the lighter Higgs boson mass may have a 2−3 GeV amount of enhancement/decrease, a signature of NHSSM. Additionally, compared to MSSM, NHSSM is able to provide with correct ranges of m h for a significantly lower value of |A t |. We must, however, note that for smaller |A t | less than a TeV or so, the contribution from A t to the radiative corrections to the Higgs mass is hardly large enough so as to satisfy the lower limit of Eq.3.1. Needless to mention a choice of a heavier third generation of squark would easily enhance m h close to its upper bound of Eq.3.1. Now, we would like to focus on tan β = 40. Since the contribution of A t is suppressed by tan β (Eq.2.3), we expect only a marginal impact on m h . This is evident from the appearance of only a small spread around the cyan region of Fig.1(b). Although we have considered A b to be vanishing, it may be noted that a significantly larger value of A b may lead to non-negligible contribution towards m h via the effect of sbottom loops. This effect is indeed enhanced by tan β and depends 7 . Similar contribution arises also from the stau loops in the presence of large A τ in the large tan β regime. However for our analysis of Fig.1 with µ < 1 TeV, we hardly expect any significant contribution to m h from the sbottom loops since the prefactor of tan β in the off-diagonal sbottom mass matrices is not too large with the given range of parameter regions considered in this analysis. This is also true for stau loops. Apart from the effect of NH trilinear parameters particularly A t on m h via top-squarks in the loops, we must remember that the other NH parameter µ may play an important role via the chargino loop contributions to the Higgs boson mass [57,58]. The latter contributions that are intrinsically negative are essentially independent of tan β [57,58]. As a result, in the region near A t = 0 in Fig.1(b) where A t is not able to influence on X t because of suppression via tan β, we find a spread of the magenta points toward the smaller direction of m h . That the effect does not depend on tan β is manifested in the similar region of Fig.1(a). Thus, larger values of µ − µ may cause a decrease in m h for both values of tan β as used in the figure. Although enhancement of m h due to NH parameters is not very significant in the large tan β limit, the impact of the above parameters on particularly the low energy phenomenological constraints like Br(B → X s + γ) is extremely important as we will see in Sec.3.2.

Effects of non-holomorphic parameters on SUSY contributions to Br(B → X s + γ)
and Br(B s → µ + µ − ) In this subsection, we would like to discuss the results of including the constraints of Br(B → X s +γ) and Br(B s → µ + µ − ) on the NHSSM parameter space. The experimental limits on Br(B → X s + γ) 7 See [56] and references therein along with references for the NHSSM and the MSSM.
The region with large A t along with µA t < 0 (µ is scanned over positive value as in Eq.3.2) is discarded via the lower bound of Eq 3.3 whereas the region with µA t > 0 is disallowed via the upper bound of Eq 3.3 9 . Thus, a large |A t | regions become unavailable in MSSM which in turn causes m h to go below the lower limit of Eq.3.1. Certainly, m h can be increased via increasing the third generation of scalar mass that would enhance the Higgs mass radiative corrections.
In contrast to MSSM, we find that the magenta region corresponding to NHSSM includes parameter points that satisfy the Higgs mass bounds in addition to the B-physics constraints.
Thus large values of A t (with preference to negative region) correspond to valid parameter zones simply because of the fact that the role played by A t is effectively replaced by A t − A t cot β in NHSSM (see Sec.2.4). Thus a scan over A t even for a large value of tan β in NHSSM is able to accommodate appreciably large values of |A t | consistent with the Higgs mass as well as B-physics constraints.

Electroweak fine-tuning and higgsino dark matter
Typically a higgsino dominated dark matter with a mass around a few hundred GeV produces extremely large annihilation cross section. Apart from the LSP pair annihilation there is a substantial amount of χ 0 1 −χ 1 ± coannihilation. A larger higgsino content in a primarily bino dominated χ 0 1 such as what one obtains in the focus point [62,63]/hyperbolic branch scenario [50] may produce the right relic abundance satisfying the experimental constraint of DM relic density. However, this is highly constrained by the direct detection of DM experiments like LUX [13]. In MSSM, a highly higgsino dominated LSP satisfies the DM relic density limits as given by Eq.3.5 from PLANCK [12] data for an LSP mass of ∼ 1 TeV [64,65].
Certainly, with a 1 TeV higgsino mass the electroweak fine-tuning estimate ∆ T otal of Eq.2.12 becomes generally large 10 . One, on the other hand, is able to isolate significantly the EWFT measure from LSP mass in NHSSM simply because of the fact that the Higgs potential (Eq.2.8) does not depend on the nonholomorphic Higgsino parameter µ , whereas the higgsino content of the LSP is determined via the difference of µ and µ .
We now scan the NHSSM parameter space focusing on a higgsino dominated LSP. We select a fixed value of 3 TeV for all the following masses namely M 1 ,M 2 and m A along with all the squark and slepton mass parameters in a scenario of varying µ − µ as mentioned below. 10 See however Ref. [66] where the authors considered specific GUT scenarios with non-universal gaugino masses.
We note that the LEP limit on lighter chargino mass is translated in NHSSM as |µ − µ | > ∼ 100 GeV and we probe an LSP mass zone up to 1.5 TeV so that χ 0 1 remains sufficiently higgsino dominated in its composition whereas we do not include the s-channel A-boson annihilation region. Furthermore, in order to have the Higgs mass in the correct range, we vary the trilinear parameters A t and A t as given below 11 .
The parameter ranges of Eq.3.6 and 3.7 apply to the discussion of the present subsection only. In Fig.3 we plot ∆ T otal (Eq.2.12) vs m χ 0 1 for MSSM and NHSSM cases drawn in blue and brown colors respectively. The MSSM part of the analysis corresponds to vanishing µ and A t and appears as a thin blue line in the middle. The relic density limits are satisfied in the vertical strip shown in magenta or green corresponding to NHSSM and MSSM cases. The lighter chargino mass limit disallows the LSP mass to go below 100 GeV (white region in the left). Fig.3(a) shows that ∆ T otal for NHSSM can either be larger or smaller than the MSSM specific values. The larger ∆ T otal region, of course, occurs when |µ| is larger while the lower region corresponds to smaller |µ|. The LSP mass is essentially same as |µ − µ |. Fig.3(b) shows a similar result for tan β = 40. The lowest value of ∆ T otal satisfying the DM relic density constraint of Eq.3.5 for MSSM is about 500, whereas the same for NHSSM for tan β = 10 is about 50 rather than being vanishingly small 12 . On the other hand, ∆ T otal for tan β = 40 can indeed approach zero for vanishingly small µ. The above difference of small µ behavior of ∆ T otal is indeed consistent with the discussion of Sec. 2.3. We have used micrOMEGAs [67] for relic density computation.
We note that a higgsino type of LSP generally satisfies the LUX data [65] and it may be probed via XENON1T [68].  12 We essentially agree with the analysis of Ref. [37].  Thus NHSSM can easily accommodate the stringent muon g−2 constraint even with a small amount of NH trilinear coupling A µ by allowing larger smuon masses. Apart from the above in relation to the effect of the combined constraints from B-physics and Muon g − 2, one finds that NHSSM can accommodate the large tan β regimes (that naturally increases a SUSY µ ) easily in comparison with the MSSM scenario. Table 1 compares MSSM and NHSSM spectra for two benchmark points in detail. Both the points for NHSSM satisfy all the relevant constraints, whereas the corresponding MSSM points do not necessarily satisfy the same.  The Higgs mass constraint of Eq.3.1 is however not imposed. The blue, green and brown regions correspond to satisfying the muon g − 2 constraint within 1σ, 2σ and 3σ limits respectively.
Finally, we will comment on the possible effect of considering a negative sign of µ. In regard to Sec.3.1 the part involving the radiative corrections to Higgs boson mass is controlled by X t = A t − (µ + A t ) cot β. Hence for large tan β the contribution coming from the part involving µ and A t is suppressed. On the other hand, for smaller values of tan β (= 10) the above contribution is relatively larger. However, since A t as well as A t both are scanned with a larger range (−3 TeV

Conclusion
In MSSM the superpotential is a holomorphic function of superfields and one considers soft SUSY breaking terms that are also holomorphic function of fields. However a SUSY theory devoid of an SM gauge singlet allows non-holomorphic soft SUSY breaking terms in the Lagrangian, and this has been used in various beyond the MSSM scenarios analyzing neutrino physics, leptogenesis, CP violation etc. In this analysis, we focus on the relevant phenomenological constraints while considering non-holomorphic soft SUSY breaking terms in a beyond the MSSM scenario with identical particle content as that of MSSM. Our work on Non-Holomorphic Supersymmetric Standard Model (NHSSM) uses electroweak scale input parameters similar to what is considered in the pMSSM model. This is unlike the previous analyses where only the non-holomorphic parameters were given at the electroweak scale while other soft parameters belonged to the grand unification scale or all the input parameters were given at the aforesaid scale.
We particularly analyze NHSSM specific effects on the Higgs mass radiative corrections, electroweak fine-tuning, electroweakino spectra, the constraint due to Br(B → X s +γ) and Br(B s → µ + µ − ) and the novel signature of NHSSM that enhances a SUSY µ so that it can easily accommodate the Muon g−2 limits even for larger smuon masses or small tan β. In the context of radiative corrections to the Higgs boson mass it is seen that for tan β = 10 and a given value of A t , the variation of NHSSM parameters, particularly due to A t , the trilinear nonholomorphic parameter may easily cause a change in m h by 2-3 TeV in either direction positive or negative. In other words, NHSSM does not necessarily require large values of |A t | in order to produce the right amount of radiative corrections to Higgs mass, since both A t and A t contribute toward the corrections. The radiative contributions to m h due to A t is suppressed by tan β leading to a quite small effect for tan β = 40. However, for large tan β NHSSM makes parameter space with large |A t | to become valid via its effects on Br(B → X s + γ). We note that with squark and gluino masses assuming just above the LHC limits and tan β = 40 the constraints from Br(B → X s + γ) in MSSM eliminates the large |A t | zones altogether so that m h goes below 122.1 GeV, the lower limit of Higgs mass considered in the analysis.
The non-holomorphic terms in NHSSM through their contributions toward Br(B → X s + γ) allows a significant part of the large |A t | region to become valid.
It is known that a higgsino type of LSP in MSSM that is supposed to satisfy the PLANCK data on DM relic density has a mass of around 1 TeV. In pMSSM this obviously increases the electroweak fine-tuning due to the sufficiently large value of µ. In contrast, NHSSM is able to produce a drastic reduction of the electroweak fine-tuning measure even for such a large mass of higgsino.
The dependence of electroweak fine-tuning on µ rather than µ , the bilinear Higgs nonholomorphic parameter whereas the fact that electroweakino masses are related to the difference of µ and µ indeed isolates the two sectors 13 . The electroweak fine-tuning can either decrease or increase depending on the relative contributions of µ and µ to the difference µ − µ . 13 As mentioned before, we agree with the result of the recent analysis of ref. [37] in this regard.
Regarding the (g − 2) µ constraint, NHSSM is able to significantly enhance a SUSY µ even for a small tan β via a small value of the associated trilinear coupling parameter A µ . This is true even for a significantly large smuon mass mμ 1 . This is indeed a very novel feature of NHSSM. Just a small amount A µ like 50 GeV even for a small tan β (= 10) may significantly alter the MSSM predictions on lighter smuon mass that would satisfy the (g − 2) µ constraint at 1σ level . This has the potential to cause a significant change in predictions involving SUSY models in general.

Acknowledgment
A.D. would like to thank the Council of Scientific and Industrial Research, Government of India and Indian Association for the Cultivation of Science for support and computational facilities. UC would like to thankfully acknowledge a correspondence with Florian Staub.