Dirac gauginos and the 125 GeV Higgs

We investigate the mass, production and branching ratios of a 125 GeV Higgs in models with Dirac gaugino masses. We give a discussion of naturalness, and describe how deviations from the Standard Model in the key Higgs search channels can be simply obtained. We then perform parameter scans using a SARAH package upgrade, which produces SPheno code that calculates all relevant quantities, including electroweak precision and flavour constraint data, to a level of accuracy previously impossible for this class of models. We study three different variations on the minimal Dirac gaugino extension of the (N)MSSM.


JHEP06(2013)073 1 Introduction
The LHC experiments have claimed the discovery of a new particle near 125 GeV [1,2]. Its production and decays make it a good candidate for the Higgs boson. However, more experimental data is needed for a precise determination of its quantum numbers and interactions to know if the Higgs sector is the Standard Model one, or an extended version, for instance as required in supersymmetric models. A clear indication of non-minimality would be a significant excess or deficit in at least one decay channel, but the Higgs mass itself can also be regarded as such, since it is somewhat high for the most natural version of the MSSM. This is one motivation for this work, in which we study the main properties of such a Higgs boson in models with Dirac gaugino masses . These have numerous virtues compared to their Majorana counterparts, not least that they allow for increased naturalness [6,23,34,[37][38][39] which is particularly important in the light of recent LHC SUSY searches [40][41][42]; but also that the direct production of gluinos is suppressed, loosening the bounds from direct LHC searches [29,31]; they allow relaxation of flavour constraints such as due to Br(B → sγ) [43]; they preserve R-symmetry (allowing for simpler supersymmetry-breaking models) and can be motivated from higher-dimensional theories as being derived from an N = 2 supersymmetry in the gauge sector at some scale.
To give gauginos Dirac masses, we must add new chiral superfields in the adjoint representation of each gauge group: a singlet S = S + √ 2θχ S + · · · , an SU(2) triplet T = O +· · · . Of these, the triplet and singlet may have new renormalisable couplings with the Higgs which allow more possibilities to obtain the desired mass range than in the MSSM. More precisely, the lightest Higgs mass in these models is determined by opposite competing effects. On one hand, the presence of new couplings in the extended scalar sector leads to an enhancement of this mass by allowing new contributions to the quartic Higgs coupling. On the other hand, the supersoft operators that include the Dirac mass induce new D-term couplings, which increase Higgs mixing and thus tend to reduce the lightest Higgs mass. However, this is only potentially problematic if the triplet scalar soft mass is small; but we shall demonstrate in section 4 -along with a general discussion of naturalness in Dirac gaugino models -that it may be naturally large enough to avoid this problem. Higgs mixing involving the singlet induced by a Dirac Bino mass, however, then presents an intriguing opportunity: it allows the decays of the Higgs to b quarks and τ leptons to be suppressed while preserving the rate to W s and Zs, and an enhancement of the diphoton rate. We give an explanation of this and a discussion of the Higgs production and decay rates in section 5.
Previous attempts to quantitatively study the Higgs sectors of Dirac gaugino models have been hampered by the lack of numerical tools to do so; until now only one-loop effective potential calculations of the Higgs mass have been possible. This is in contrast to the MSSM, where the leading corrections are known to three-loop order. However, with an upgrade of the SARAH package described in section 3, it is now possible to study models with arbitrary gauge groups and matter content: it can generate SPheno code which calculates all one-loop pole masses and tadpoles, which allows a much more accurate determination of the Higgs mass. We implement a minimal Dirac-gaugino extension of the (N)MSSM

JHEP06(2013)073
which lends itself to four particularly interesting sub-classes of models, which we review in section 2. One such class of models is the "MSSM in disguise," where all the new scalars are too heavy to observe or mix substantially with the Higgs; this is a good toy scenario to use to test of the code, and we do just that.
The SPheno code produced by SARAH can also calculate the branching ratios and production cross-sections of the Higgs, as well as electroweak precision observables such as ∆ρ and flavour constraints such as Br(B → sγ). We take full advantage of the latter in investigating the "MSSM without µ term" [7] in section 6, comparing the results for ∆ρ with approximations given in [7]. Unfortunately those constraints in addition to those on chargino masses and the Higgs mass yield that model problematic. However, we propose a new model, which we call "dynamical µ models," in which the singlet obtains a substantial expectation value -we show in model scans in section 7 that this not only alleviates all of the problems of the MSSM without µ term, but also allows for Higgs mixing and thus interesting deviations in the Higgs production rates and branching ratios that may better fit the current data than the standard model.

One model into four
We shall consider the framework of e.g. [18,23] of adding minimal fields to the MSSM to allow Dirac gaugino masses. As described above, this requires three adjoint chiral superfields S, T, O and we include the supersoft operators However, adding Dirac gaugino masses to the MSSM can introduce many extra parameters: not just the Dirac masses themselves, but also new superpotential couplings and soft terms. By making certain assumptions about the new parameters we can then arrive at different limits of the model with different phenomenology; we shall consider four such limits.
The most general renormalisable superpotential that we can write down is The usual scalar soft terms are Similarly there are the A-terms for W 2 . One limit of the above model would be to allow all parameters, including Majorana gaugino masses, to be significant and non-vanishing. Such models may have virtues due to the extra Higgs couplings and extra gaugino states (the charginos could contribute, for example, to enhancing the Higgs-to-diphoton decay channel) but we shall leave the exploration of this to future work. Instead, we shall explore models where the gaugino masses are entirely Dirac, taking as motivation the possibility of preserving R-symmetry in some sector of the theory (but not in all: it must ultimately be broken in any case). Of particular interest to us are: • MSSM in disguise: here we shall allow a µ-term, and assume that the only source of R-symmetry violation arises in the supersymmetry-breaking sector, but permit only a B µ term. This assumption will be preserved by renormalisation group running and so it can be justified by high-energy boundary conditions; non-zero A-terms would lead to Majorana gaugino masses. Since it is the "MSSM in disguise" we shall take the scalar singlet and triplet to be very massive (several TeV). We perform some scans over models of this type in section 3.
• MSSM without µ term ( / µSSM): this is the scenario of [7], similar to the MSSM in disguise but taking µ = 0. Here we shall insist that the singlet vev is small, so the chargino mass must be generated by the coupling λ T . We investigate this scenario in section 6.
• Dynamical µ models: in this scenario, we again take µ = 0 but allow a substantial non-zero expectation value for the singlet. In this way, the vev and the coupling λ S lead to an effective µ-term. Models of this type are very natural and interesting from the point of view of Higgs mixing: we perform scans over these models in section 7.
• Dynamical µ and Bµ models: this is the scenario of [23] where we allow a non-zero κ, breaking R-symmetry in the visible sector, but allowing µ and B µ to be generated via a non-zero singlet vev -so we can set all R-symmetry-breaking parameters in the supersymmetry breaking sector to zero. It is somewhat similar to the NMSSM, but the Dirac masses lead to some interesting differences. Models of this type can be very natural, but we leave scans of their parameter space to future work. We review the tree-level properties of the generic case in appendix A; see also [23]. However, common to all of the above scenarios is the assumption that Dirac gaugino masses JHEP06(2013)073 dominate over Majorana ones. For this to be natural, A-terms must also be small, so in our searches we shall set (unless otherwise stated) Some notation. We now introduce some notation relevant for the following: we redefine the singlet and triplet scalars in terms of real components S ≡ 1 The expectation values v S , v T are associated with new non-trivial potential minimisation conditions, which we give in equation (A.7). The scalars s R , T R mix with the Higgs fields, with a 4 × 4 mass matrix given in equation (A.12); the mixing will be important in section 7.

SUSY search constraints
This paper is substantially concerned with the Higgs data from the LHC (and the Tevatron), which is summarised in appendix C, along with the relevant electroweak precision and rare decay constraints. However, the LHC (along still with LEP) has also put several limits on supersymmetric particles, which we must also take into account in model searches. Since these will be relevant for all of the following sections, we shall summarise here the relevant constraints on our models due to recent searches at the LHC. Motivated by naturalness, our models typically have several neutralinos and charginos in the hundred to a few hundred GeV mass range; these then influence the bounds on the squarks and Higgs.

Squarks and gluinos
In our models the gluino is a Dirac particle, and this will therefore affect the bounds on its mass -and also those of the squarks. One scenario would be to take a very heavy gluino, where the bounds on the squarks soften [29,31]; however, in the opposite limit of a lighter gluino with heavy squarks, the bounds on the gluino mass actually strengthen. In the Majorana case, the bounds on first two generation squarks can be as high as 1.5 TeV [44] and the gluino about 1.24 TeV [45]. In the absence of an updated study for the Dirac case, we shall conservatively take the bound on first two generation squarks to be the same as in the Majorana case, and the gluino to be heavier than 1.5 TeV. However, it should be stressed that the gluino mass, coming from a supersoft operator, does not affect the other masses to the same extent as in the Majorana case.

Stops
Stops are now excluded up to 660 GeV for a massless neutralino [46], and indeed this bound remains roughly constant for a lightest neutralino mass up to about 200 GeV. Hence in our models, when we take the stops to be light we shall set them at 700 GeV.

Charginos and heavy neutralinos
The bounds on charginos and neutralinos heavier than the lightest one (the LSP) depend on whether the decays proceed via sleptons (rather than electroweak gauge bosons) and the mass of the neutralino [47]. If the sleptons are light so that the decays proceed via them, then the bounds are very strong, excluding charginos and heavy neutralinos up to JHEP06(2013)073 600 GeV for LSP masses up to about 300 GeV (after which the bounds dramatically fall away). However, if the sleptons are much heavier than the electroweak scale (and indeed the charginos/heavy neutralinos) then the bounds are much softer -up to 315 GeV if the LSP has mass 110 GeV. If the LSP is heavier than this, the bounds essentially vanishand we are merely left with the LEP bound of 105 GeV for the chargino mass. Hence we shall simply take this to imply that the sleptons should be relatively heavy compared to the charginos and the LSP should be above about 110 GeV.

Implementation in SARAH and SPheno
For our numerical studies we have extended the Mathematica package SARAH [48-51] to support Dirac gauginos. 1 This support includes an automatic derivation of all specific D-and fermion mass terms in the context of Dirac gauginos as well as the full two-loop RGEs using the results of ref. [39]. Both aspects are also covered in the presence of several Abelian gauge groups and kinetic mixing using the approach of ref. [52].
Afterwards we implemented our model and used the possibility of SARAH to create source code for SPheno [53,54]. For this implementation we used the most general setup allowing all possible superpotential and soft-breaking terms of eqs. (2.2)-(2.5). Also Majorana as well as Dirac mass terms could be present at the same time. The different limits of the model were obtained by a convenient choice of the input parameters in our scans. As mentioned above, it is also possible to create a SPheno version with SARAH which includes the two-loop RGEs. In this case it is possible to study also constrained GUT models. However, we leave this for further analysis and treated nearly all parameters at the SUSY scale as independent. The only relation between parameters stems from the necessity to be at the minimum of the vacuum: since there are four non-trivial vacuum minimisation conditions (given in equations (A.5), (A.6) and (A.7)) we must use these to eliminate four parameters. From a top-down perspective, it would be preferable to specify the soft masses m S , m T as was done in [23], and derive from that v S , v T . However, the equations are nonlinear in these, and so it is preferable to instead take v S , v T , tan β, µ as input parameters in the code, and treat m 2 S , m 2 T , m 2 Hu , m 2 H d as output parameters (calculated including oneloop tadpoles). As further inputs for our studies we the soft-breaking terms as well as the superpotential couplings at the SUSY scale: the standard model gauge and Yukawa couplings are calculated at M SU SY using two-loop standard model RGEs from M Z [55].
The obtained code including these boundary conditions was compiled with SPheno 3.1.2 and provides a fully functional spectrum calculator: the entire mass spectrum is calculated at one-loop using the full dependence on the momentum of the external particle. These calculations include the entire CP and flavor structure and, of course, also the mixing between MSSM fields and the new states. For a detailed discussion of the calculation of the loop corrected mass spectrum using SARAH and SPheno for extensions of the MSSM we refer the interested reader to ref. [56]. We note that, in the absence of the theoretical JHEP06(2013)073 calculation, the code cannot include the two-loop corrections to the Higgs mass involving (Dirac) gluino masses that can be important in the MSSM and NMSSM. In those cases, the Higgs mass is usually increased by 2 to 4 GeV. Throughout we shall be conservative and allow a variation of ±4 GeV for the mass of the Higgs in the scans, but we expect that typically the shift will be positive.
In addition, the produced SPheno version provides routines to calculate the decay widths and branching ratios for all SM superpartners and Higgs fields. In general, these are pure tree-level calculation. However, in the case of the Higgs particles the loop induced decays into two photons and two gluons are calculated. Here, not only the known SM or MSSM contributions are included but SPheno takes automatically the loops arising from all charged particles in the model into account. This calculation is comparable to the analytical results given in section 5, but also the NLO corrections due to quarks and squarks are added as given in ref. [57]. Similarly, for the Higgs decays into quarks the dominant QCD corrections due to the gluon are added [58]. Finally, the SPheno version for Dirac gauginos also calculates the observables b → sγ and ∆ρ with the same precision as done in the MSSM by SPheno 3.1.1 including all contributions due to the new states present in the considered model. The details of these calculations can be found in ref. [54] and references therein. With version 3.3.0 the SPheno output of SARAH provides also of a full one-loop calculation of B 0 s,d → l + l − [59]. Of course, even if the used for this study mostly the SPheno output of SARAH, models with Dirac gauginos are supported also in all other interfaces of SARAH. Hence, the user can use SARAH to obtain model files for CalcHep/CompHep [60,61], WHIZARD/OMEGA [62,63], FeynArts/FormCalc [64,65] or in the UFO format [66] which is supported by MadGraph 5 [67]. The model files for CalcHep can also be used together with MicrOmegas for dark matter studies [68].

Comparison with effective potential
It is possible to calculate an approximation for the Higgs mass via the effective potential technique. In the decoupling limit of large Bµ, this yields where the first line is the tree-level mass, the second line is the usual contribution from stops with mt 1,2 the physical masses and A t set to zero; the λ i are the coefficients of the terms in the most general CP-conserving gauge-invariant potential up to quartic order (e.g. λ 2 is the coefficient of the |H u | 4 operator) about zero vev [69]. This is a good approximation when the Higgs vev is small compared to the energy scales being integrated out; in this Approx m h GeV Figure 1. Comparison of effective potential method (labelled "approx m h ") with a parameter scan using SPheno code produced by SARAH. Parameters are given in the text; a searched was performed using the SPheno code, restricting to m h = 125 ± 4 GeV, and from the same input parameters for each model the one-loop Higgs mass is shown (note that this has a larger range). Only models with relatively low couplings λ S , λ T (roughly λ 2 case, we shall consider integrating out the adjoint scalars, which is appropriate for the MSSM in disguise. We give the full expression for the coefficients λ i in appendix B, but an interesting limit is to consider B S = B T = 0 and neglect v, v S and the Dirac masses m Di ; then the contribution from the singlet and triplet scalars is [23] where in the last line we show the limit that the scalar masses become equal. By performing a scan over parameters we can compare this approximation with the more accurate results from SPheno using the code produced by SARAH; we show the results with the choice of tan β = 50, m D2 = 600 GeV, first two generation sfermion mass squareds of 4 × 10 7 (GeV) 2 , third generation sfermion mass squareds 4 × 10 6 (GeV) 2  varied between 10 6 and 10 10 (GeV) 2 . As can be seen from the figure, there is good agreement between the two, although of course the approximate effective potential calculation exhibits a wider variation of masses; that there is no apparent correlation is unsurprising, essentially reflecting the error margin in the effective potential method.

JHEP06(2013)073 4 Dirac gauginos and natural SUSY
One of the increasingly attractive features of Dirac gaugino models is that the supersoft operators allow for increased naturalness [6,23,34,[37][38][39][40][41]; they do not appear in the RGEs for the soft masses and so they only affect the stop mass via a UV-finite correction, allowing for heavy gluinos with light stops. On the other hand, the Higgs potential is corrected at tree-level by two competing effects: one is the enhancement of the Higgs mass (at low tan β) by the new couplings λ S , λ T , clearly evident in equation (3.1); the other is a reduction in the effective D-term Higgs quartic coupling due to the Dirac mass terms; if we ignore the superpotential couplings and integrate out the heavy scalars then [6,23] This latter effect manifests itself as mixing terms in the Higgs mass matrix (A.12) ∆ hs , ∆ ht ; provided that the other soft masses are large enough the suppression can be avoided. We may then ask how large these soft masses can be while preserving naturalness, to which we consider the radiative corrections to m 2 H u,d : where Λ is the UV cutoff of the theory. Using ∆ ≡ The latter value is particularly important, because we also have the constraint that the triplet scalar expectation value must be small in order to avoid a large ρ parameter. We have (see appendix) which is satisfied for any value of m D2 for m T 1.4 TeV; but interestingly is also satisfied for any value of m T if m D2 1.4 TeV. An interesting choice for the parameters λ S , λ T would be for them to take their N = 2 supersymmetric values at some scale, 2, although corrections due to the running would break this exact relation. However, even if the N = 2 scale is intermediate (such as 10 12 GeV) assuming a desert we would still find λ S 0.2, λ T 0.4 at the low energy scale. In this case, in order to preserve naturalness we would require In order to satisfy the tree-level ∆ρ constraint, we would either need to work at small c 2β , or take m D2 1.4 TeV and ensure that the ensuing Higgs mixing allows a large enough Higgs mass. We leave exploring this interesting possibility to future work: in this paper we shall take large values of m T and small values of λ T to keep the loop-level corrections to ∆ρ small. One final naturalness-related issue in these models is that the requirement of small A-terms typically diminishes stop mixing; in the case of large tan β the mixing is almost eliminated. Hence the contribution from the stops to ∆ρ is [70][71][72]: which is similar to the experimental value. In "natural SUSY" MSSM models, to survive the current direct search constraints the stops must be about 700 GeV (or be near the top mass), and so in the absence of large A-terms (as we are assuming) the stops by themselves will not be able to lift the Higgs mass to 125 GeV, but have a large impact on the electroweak precision corrections. On the other hand, in the context of λSUSY [72] the model can remain natural for stops as heavy as 1.5 TeV [72,73] because the relative contribution of stops to the Higgs mass compared to the tree-level effect is small. Clearly λ S in our models is the same as the λ in λSUSY. Hence the simplest natural scenario is to take small tan β, small λ T , and λ S ∼ 1 with m T ∼ few TeV and m S TeV. The Dirac gaugino masses m D3 , m D2 can be naturally large, but the Dirac Bino mass will be bounded above by the amount of mixing that it induces between the singlet and lightest Higgs (thus reducing the Higgs mass) to be a few hundred GeV. As a result of this, an amusing feature is that natural Dirac gauginos will lead to a Majorana neutralino, since there will be non-negligible mixing between the Bino and the neutral Higgsinos.

Higgs production and branching ratios
It is now important to consider the production cross-sections and branching ratios of the Higgs. In our Dirac gaugino models there is a singlet scalar which may mix with the lightest Higgs state, so we shall consider the branching ratios into each channel taking into account the mixing, providing some approximate expressions and then comparing them to the output of the SPheno code created by SARAH.
We shall use the standard definitions where A ii is the amplitude.

JHEP06(2013)073
To take Higgs mixing into account, consider the rotation of the states h, H, s R via a matrix S so that  We shall assume throughout that the lightest Higgs field h 1 has mass 125 GeV, i.e. there is no additional lighter singlet. In this notation, we can then calculate how the production and decay channels are modified. At 125 GeV the production cross-sections (as listed on the CERN yellow pages [74]) are Of the initial eigenstates, only h couples at tree level to the vector bosons. The coupling to gluons and in the ttH process is proportional to the top quark coupling squared, so (noting that S ij are the elements of the above matrix S, so that S 11 denotes the amount of the original state h in the physical lightest Higgs, and S 21 is the amount of H) Hence we can write We can use the same approach for the decay branching ratios: Since the photon couples at one loop to the Higgs, and the singlet couples to charged fields, the expression for the coupling to the photon will be more complicated. At m h = 125 GeV, the standard model Higgs branching ratio (as listed on the CERN yellow pages [74]) is dominated by So we can write To approximate the diphoton channel, we require the expression where the functions A s where s is the spin of the field in the loop are standard and given in [57]; τ X ≡ 4m 2 X /m 2 h 1 . For fermions and scalars they are well approximated by the large mass limits of 4/3 and 1/3 respectively, but for the W boson and top quark the values are A 1 (τ W ) −8.32, (4/3)A 1/2 (τ t ) 1.84. If we consider that the singlet couples to some charged Dirac fermion we can write for the couplings g h 1 ii We can then consider enhancements through various fields remaining light, taking the mixing into account.

Charginos
The chargino mass matrix is expanded by new charged states from the triplet: in the basis

JHEP06(2013)073
Since the loop function A 1/2 varies very little between the lower bound on chargino masses from LEP (105 GeV, or 92 GeV with caviats) and infinite mass, it is very well approximated by 4/3 over all cases of interest. Then we can well approximate In the limit M 2 = M T = λ ST = 0, this simplifies to The scenario of light charginos is particularly appropriate for the MSSM in disguise; since the adjoint scalars are all massive the main phenomenological difference with the MSSM will be the charginos and neutralinos. Hence this represents one useful limit of this formula, where there is negligible mixing between the Higgs states (i.e. S 11 1, S i1 0 ∀ i = 1). Then it is clear that for an appreciably large coupling λ T and large c 2β the γγ channel can be enhanced without affecting the other channels. We performed at scan at λ T = 1, tan β = 50 varying µ, and plot the contours of the Higgs to diphoton branching ratio in figure 2, showing also the effect on ∆ρ. This analysis is similar in spirit to models adding a triplet to the MSSM [75], except instead of including a Majorana Wino and supersymmetric triplet mass we have included a Dirac Wino mass.
An alternative application of formula (5.13) is to allow appreciable Higgs mixing but take the large m D2 limit, leaving light Higgsinos. We then find (5.14) Hence the diphoton production rate can be enhanced for suitably large vλ S /µ and S 31 . This is applicable for the NMSSM too, and is particularly interesting because by varying the Higgs mixing terms we can simultaneously enhance µ γγ , decrease the µ bb and µ ττ while maintaining µ W W and µ ZZ roughly unchanged if so desired. This is similar to singlet extensions of the MSSM without Dirac gauginos [76] and can be easily understood from equations (5.5) and (5.7): by having a small positive admixture of H we enhance the coupling of the Higgs to tops, and hence to gluons; by reducing S 11 we reduce the coupling to bottoms and taus, which reduces the total width of the Higgs -both of these can compensate for the reduction in coupling to W s. In figure 3 we show how µ bb and µ τ τ are affected by the mixing. The above effect of the mixing and charginos on µ γγ is then illustrated in figures 4 and 5.
JHEP06(2013)073  Figure 2. Results of a scan over a subspace of the "MSSM in disguise" with λ T = 0.7, tan β = 50, varying m D2 and µ in order to vary the lightest chargino mass which is restricted to be Dirac Winolike (above the contour the lightest chargino is Higgsino-like). The stop and sbottom masses were taken to be equal and varied in order to keep the Higgs mass at 125±4 GeV (within reasonable accuracy of the results of the SARAH-produced SPheno code, which cannot include two-loop effects due to the Dirac gluino). Other non-zero soft parameters were: m D1 = 200 GeV, m D3 = 1000 GeV, λ S = 0.1 (the gluino mass is technically just too low for current limits, but we have verified that raising the gluino mass above current bounds does not noticeably affect the result); the slepton and first two generations of squark masses squared were 4×10 7 (GeV) 2 ; the singlet and triplet scalar masses were approximately 5000 GeV; Bµ was varied over [300, 5 × 10 5 ] (GeV) 2 . Left: contours of µ γγ (GeV). Right: µ γγ vs ∆ρ; orange points have the lightest chargino mass greater than 105 GeV, blue ones smaller. The value 0.0008, usually taken as a (2σ) upper bound, is shown as the red dashed line.

Charged Higgs
The charged Higgs fields can also contribute to the diphoton rate in these models. The full charged Higgs matrix involves not just the H ± of the MSSM, but also the charged triplet scalars; however, from the ρ-parameter constraint we know that the triplet scalars must be heavy (see section 4), so we shall neglect their contribution. Hence we can approximate A large contribution to the diphoton rate from charged Higgs loops can then arise when mixing between the lightest Higgs and the singlet is substantial. Note that light charged Higgs fields in the limit of very small tan β also often demand light stops and charginos to cancel large contributions to b → sγ.  3, and the red contour is a rough 95% confidence-level preferred region according to the data given in appendix C (excluding electroweak precision data, which is dependent on the spectra of any other light particles).

Stops and staus
Of the squarks in the theory, the stops and staus -having the strongest couplings to the light Higgs -have been studied in the MSSM as candidates to modify the Higgs to dipho-JHEP06(2013)073 ton rate [77,78]. In models with Dirac gauginos, the new D-term contributions to the potential modify the squark masses; for the stop and stau their mass matrices are given in appendix equations (A. 19) and (A.20). Neglecting the normal D-term components, v T and the A-terms gives approximately It is straightforward to derive expressions for the the couplings of the stops and staus to the Higgs eigenstates, but since the full expressions are lengthy we give here simplified formulae neglecting subleading terms proportional to M Z and setting the A-terms and v T to zero: Figure 6. Here mt i , mτ i are the masses of the stop and stau eigenstates respectively. The Dirac mass enters into the above by shifting the mases (e.g. m 2

MSSM without µ-term
In the presence of Dirac gaugino masses it is possible to remove the Higgsinos' mass term from the superpotential of the MSSM [7]. This is another way of curing the intrinsic µ-problem of the MSSM. Furthermore, an approximate U(1) R symmetry naturally guarantees that tan β is large, explaining the top/bottom quark mass hierarchy. In contrast to its appealing theoretical aspact, the / µSSM is under substantial pressure from experimental data. First, LEP put a lower limit on the mass of the lightest chargino of 94 GeV [79]. The chargino mass eq. (5.10) reads in this limit It is known that it is possible to fulfill this bound by a careful choice of λ T and m 2D : a large value of λ T as well as m 2D around 107 GeV is needed to maximize the mass of the lightest chargino, see figure 6. The highest mass which can be reached at tree-level is about 110 GeV. This is also not improved at the one-loop level because the loop corrections due to the heavy triplets even tend to decrease the mass. Hence, the highest mass we could find in our scans calculating the full one-loop spectrum was 103 GeV. This is already under heavy pressure using the latest results for electroweakino searches at LHC and LEP, see section 2.1.3. Hence, one would need special mass configurations like compressed spectra in order to soften the mass limit and somehow keep this model alive. However, a more severe problem is that the large values of λ T lead often to a huge contribution to ∆ρ. The reason is that this coupling breaks the custodial SU(2) L symmetry present in the SM and MSSM: especially, the chargino and neutralino-loops contribute differently to the W and Z self-energies. The large impact was already pointed out in ref. [7] using an approximation for the resulting contributions to the T -parameter. We repeat this analysis with the full numerical evaluation of ∆ρ. We find a strong correlation in this model between the mass of the lightest chargino and the smallest possible value of ∆ρ as depicted in figure 7.

JHEP06(2013)073
Even if the full calculations leads to somewhat smaller values of ∆ρ than the approximate one for λ T > 1, all points which fulfill the limit of mχ+ 1 > 94 GeV suffer from a large ∆ρ of at least 0.003. Many points are even above 0.01. Note that ∆ρ is very quickly increasing with the chargino mass. Therefore, demanding ∆ρ < 0.0008 would rule out all points with chargino masses above 20 GeV. Furthermore, it can also be seen that in general the points with a Higgs mass between 122 and 128 GeV lead independently from the chargino mass to ∆ρ > 0.001. The reason is that a sizable contribution of λ T to the tree-level Higgs mass is needed. This could, of course, be circumvented to some extent by allowing for even larger values of T top . However, this is in contradiction to the approximate U(1) R -symmetry which suppresses in general the trilinear terms. Therefore, if we restrict ourself to moderate values of the squared squark mass parameters and trilinear soft-term, this model is always in conflict with ∆ρ: even if it would be somehow possible to find kinematical configurations to significantly reduce the LEP limits on the chargino mass, the now existing bound on the Higgs mass still predicts too large values of ∆ρ.

Dynamical µ models
Although the "MSSM without µ term" may be severely challenged, by allowing a substantial expectation value for the singlet the various problems can be cured. The Higgs potential will always lead to a non-zero value for v S as can be seen from the minimisation condition (A.16), but in order for this to be significant we can allow a non-zero negative JHEP06(2013)073 value for B S and/or a non-zero tadpole term t S . Both of these are generically present and do not break R-symmetry. In this scenario, as in the MSSM without µ term, we take the only significant source of R-symmetry breaking to be a Bµ term.
This scenario is particularly interesting from the perspective of Higgs mixing, since the singlet adjoint scalar will typically be light -we see from the minimisation conditions that the parameterm 2 S in the tree-level mass-squared matrix (A.12) is and, expecting from naturalness and RGE running [39] There may thus be substantial mixing between the singlet and original "h" eigenstate; the size of Bµ term then controls the amount of "H" in the lightest Higgs mass eigenstate.
Moreover, the singlet couples to the gauginos via the coupling λ S , which, if m D2 is not small, will be predominantly Higgsino-like. This then offers the possibility of realising the scenario considered in section 5.1. We have therefore conducted a scan over a portion of the parameter space of these models, concentrating on models with a small component of mixing between h and H but substantial S 11 and S 31 components, using the SPheno code produced by SARAH. tan β was taken to be 1.5 and λ S was varied from a negative initial value in order to fix the Higgs mass at 125 ± 4 GeV -recall that this is a rather conservative error range. The other parameters varied were m D1 ∈ [−800, 800] GeV, v S ∈ [130, 430] GeV, Bµ ∈ [312, 90312] (GeV) 2 , while the non-zero fixed soft parameters were λ T = 0.021, B S = −5 × 10 5 (GeV) 2 , t S = −1.5 × 10 7 GeV 3 , m 2 T = 2.5 × 10 7 (GeV) 2 , m 2 O = 9 × 10 6 (GeV) 2 , m D2 = 600 GeV, m D3 = 1500 GeV. The slepton and first two generations of squark masses squared were 4 × 10 7 (GeV) 2 while the third generation squark masses squared were 1.5 × 10 6 (GeV) 2 . Figure 8 shows the µ W W and µ τ τ values versus µ γγ (µ ZZ and µ bb being, as per the approximate formulae, almost identical to µ W W and µ τ τ respectively) while also giving the values of ∆ρ and BR(b → sγ). As can be seen from the plots, there are many experimentally viable model points. To further elucidate the comparison with our predicted scenario, a plot of the mixing parameters S 11 and S 21 with µ γγ revealed by the colour of the points is shown in figure 9.
It is worthwhile to pick out one example point; this has λ S = −1.5, Bµ = 7.57 × Models with light stops. In the previous scan, we held the third generation masses to be heavy to be above search bounds and to diminish their contribution to ∆ρ whilst still remaining natural. However, it is also straightforward to find models of the above class that have light stops, which would be natural even for the MSSM but also interesting for

Conclusions
Dirac gaugino models are gaining increased interest as non-minimal supersymmetric standard models with enhanced naturalness compared to the MSSM and an enhanced Higgs mass that can also relax bounds on direct superpartner searches. With the latest update to the SARAH package, it is now possible to study such models quantitatively using modern numerical tools, and this work is a first step in exploring phenomenologically the lowenergy parameter space. We have discussed the properties of three different Dirac gaugino scenarios that are subclasses of the minimal Dirac gaugino extension of the (N)MSSM: the "MSSM in disguise," the "MSSM without µ term" of [7] and a new scenario involving a dynamical µ term. While the first of these is phenomenologically very similar to the MSSM with higher-dimensional operators, we found that the second is unfortunately severely challenged by the current data. The third scenario, on the other hand, can be particularly natural and also has many characteristics appropriate to allow Higgs mixing and thus modifications of the Higgs production and decay rates; in particular, it is possible for example to enhance the diphoton signal, suppress the bottom and tau signals, while leaving the Z and W channels roughly the same as the Standard Model. We have performed a first examination of its parameter space but clearly it would be interesting to examine it further, particularly as new Higgs data becomes available.
There are now many interesting directions for future work. One will be to compare specific models directly with collider data, particularly in the context of models with light JHEP06(2013)073 stops. In addition, it would be interesting to see how embedding the models we have discussed in particular high-energy completions affects the discussion of naturalness. Furthermore, constraints due to dark matter (assuming a thermal history of the universe or otherwise) can now also be applied. On the technical side, to further refine the precision of the Higgs mass, the leading two-loop corrections involving the Dirac gluinos should now be calculated. This work is therefore one step on the increasingly attractive path of bringing the phenomenology of Dirac gauginos closer to the level of understanding of the (N)MSSM.

Acknowledgments
We thank Jong Soo Kim, Nicolas Bernal and Werner Porod for fruitful discussions. MDG was supported by ERC advanced grant 226371. He would like to additionally thank Veronica Sanz and Michael Trott for helpful discussions, and Chris Wymant for pointing out some typos. KB is supported in part by the European contract "UNILHC" PITN-GA-2009-237920.

A Tree-level parameters of the model
In this section we summarise the tree-level parameters of the model; see also [23]. Here we add the tadpole term for the singlet and expressions for the stop and stau mass matrices including the new Dirac gaugino D-term corrections.

A.1 Higgs potential
It will be useful to introduce the following effective mass parameters:

A.1.1 Equations of motion for the CP-even neutral fields
The scalar potential for the CP-even neutral fields is given by: where the effective masses for the real parts of the S and T fields read: There is no restriction on the sign of the different mass parameters m 2 S and B S at this stage.

JHEP06(2013)073
The imaginary parts of the fields have been dropped as their vevs are vanishing due to the assumed CP conservation [15]. The coefficients of the corresponding quadratic terms: do not, in contrast to the CP-even partners, receive contributions from D-terms proportional to the Dirac masses.
As is customarily done for the (N)MSSM, the minimization of the scalar potential allows here also to expressμ andBµ as a function of the other parameters: and The new equations are We can use these to solve for the masses in terms of the vevs. However, since the vev of T contributes to the W boson mass, the electroweak precision data give important bounds on the parameters of the model. For instance, using ∆ρ = (4.2 ± 2.7) × 10 −4 [80][81][82][83][84], we require: which is satisfied for v T 4 GeV. For large triplet masses we have

A.1.2 Masses of the CP even neutral scalars
Introducting the notatioñ

JHEP06(2013)073
the mass matrix for the CP even scalars in the basis {h, H, S R , T 0 R } is:  where we have defined: which vanishes when λ S and λ T take their N = 2 values [10]. We denote non-diagonal elements describing the mixing of S R and T 0 R states with the light Higgs h: stand for the corresponding mixing with heavier Higgs, H.
Then using the minimisation conditions we can write Clearly ∆ hs M 2 Z is necessary to prevent a see-saw reduction in the lightest Higgs mass. We havem

A.2 Squark masses
The squark masses are modified by the Dirac mass terms via the D-term contribution. We give here the expressions for the mass matrices for stops, sbottoms and staus: The entries of the sbottom mass matrix read The entries of the stau mass matrix read

B One-loop effective potential
The standard general form of the Higgs potential up to quartic order is In Dirac gaugino models where we integrate out the heavy singlet and triplet scalars, we find λ 3 and λ 4 have a tree level contribution, so we write λ 3 = 2λ 2 T + λ 3 , λ 4 = λ 2 S − λ 2 T + λ 4 where λ 3 , λ 4 are the loop corrections to the potential. Then we find at one loop λ 4 = λ 5 = λ 6 = λ 7 = 0, with the remaining contributions from the singlet and triplet scalars given by  1.14 ± 0.43 µ τ τ 1.1 ± 0.41) 0.8 ± 0.7 1.02 ± 0.35 Table 2. Table of latest production cross-section times branching ratios over standard model values for CMS [85], ATLAS [86] (both using nearly 20 fb −1 ) and the combined Tevatron search [88].

JHEP06(2013)073
To generate the crude 95% confidence level limits in the text, we simply take the observables from table 2, combine with the above values for ρ and Br(B → X s γ), compute the value for χ 2 , and compare to the 95% limit from the χ 2 distribution with the corresponding number of variables. Since we are doing a crude analysis, where there are duplicate values for the same observable we simply take the weighted mean and add the errors in inverse quadrature (i.e. 1 σ 2 = 1 , with one exception: since this paper is partly concerned with methods of obtaining excesses from the Higgs diphoton rate, we exclude that data from the fit, since the combined CMS and ATLAS value now barely shows an excess. Since the two values differ substantially compared to their errors, we take the opinion that it is not reasonable to combine them and that the CMS value may in principle be revised again. Open Access. This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.