W boson mass in minimal Dirac gaugino scenarios

We investigate the conditions for alignment in Dirac Gaugino models with minimal matter content. This leads to several scenarios, including an aligned Dirac Gaugino NMSSM that allows a light singlet scalar. We then investigate the compatibility of minimal Dirac Gaugino models with an enhanced W boson mass, using a new precise computation of the quantum corrections included in the code SARAH 4.15.0.


Introduction
While constraints on heavy Higgs bosons in supersymmetric models are rather stringent at large tan β, excluding masses above a TeV, at small to moderate tan β direct searches do not place significant limits; only indirect constraints from B → sγ limit a heavy charged Higgs boson to be above 568 GeV, roughly independent of tan β [1].On the other hand, in the Minimal Supersymmetric Standard Model (MSSM), this region is likely excluded for an additional neutral Higgs boson below a few hundred GeV due to modifications to the SM-like Higgs boson couplings.This has led to a lot of interest in extensions of the MSSM (or variants of the Two Higgs Doublet Model) where alignment without decoupling is possible [2][3][4][5][6][7][8][9], that is where the mixing between the SM-like Higgs boson and other scalars is minimised so that it aligns with the expectation values and has SM-like couplings.
These models have many interesting phenomenological properties, and have been extensively studied in the literature.They involve, at a minimum, an extension of the MSSM by three adjoint chiral superfields, one for each gauge group; the fermions from these pair with the gauginos to give them a Dirac mass.This means the presence of new scalar fields, in singlet, triplet and octet representations.
Actually, it has not adequately been investigated to what extent the singlet could be light in such models.One condition for this to be the case is that it should not disturb the couplings of the light Higgs -in other words, we should have some amount of alignments without decoupling.ATLAS and CMS both give constraints on the overall signal strength of the Higgs boson µ to be [71,72]: µ =1.06 ± 0.07 (ATLAS), µ = 1.02While this still allows a moderate amount of mixing, the larger the mixing between the flavour eigenstates, the stronger the direct search bounds on the singlet will be.
Therefore in this work we will consider the conditions for an approximate alignment in which the light Higgs mixes neither with the Heavy Higgs nor with the singlet (the triplet being decoupled).This will lead to a scenario that we refer to as the aligned DGNMSSM.
Recently, the CDF experiment reported a new measurement of the mass of the W boson [73].Compared to the SM prediction [73][74][75], this gives as averages (combined Tevatron+LEP [73,[76][77][78][79][80][81][82][83]): If we take the central value of the top quark mass to be 172.89 GeV then the central SM prediction becomes 80352 MeV [84].These differ by 7 standard deviations, although measurements at the LHC [85,86] also differ from the combination of Tevatron+LEP by 4 standard deviations, so at this stage confirmation is required by other experiments.Nevertheless, a modification to the W boson mass is one of the most generic effects of new light particles coupling to the electroweak sector, so such a hint is tantalising.It has generally been assumed in Dirac Gaugino models that the adjoint scalars should be heavy; indeed, the requirement that the triplet scalar vacuum expectation value (vev) must be very small compared to the Standard Model Higgs one -otherwise it would generate a large ρ parameter -is usually ensured by giving the triplet a heavy mass.Amusingly, following the new measurement, the simplest explanation for the enhanced W boson mass is exactly an expectation value for the neutral component of such a triplet.In this work we shall investigate that possibility in minimal Dirac Gaugino models.
Such a triplet scalar also comes along with electroweak fermions, which can modify the quantum corrections too.Therefore a precise computation is required.While a preliminary such computation was performed for the MRSSM [87] using an update to FlexibleSUSY [88,89], and a related computation was performed for the same model in [90], that model lacks a natural enhancement to the W boson quantum corrections.In this work we introduce a similarly precise computation in the package SARAH-4.15.0 and use it to examine the compatibility of our aligned DGNMSSM, along with four other scenarios -the "MSSM without µ term," the MDGSSM, the aligned MDGSSM and the general DGNMSSM -with the new measurement of the W boson mass, or a naive world average value of M world average W = 80411 ± 15 MeV.This work is organised as follows.In section 2 we summarise the essential details of the class of Dirac Gaugino models, including the vacuum minimisation conditions and mass matrices.The conditions for alignment are reviewed and a comparison is made with the cases of the MSSM and NMSSM.We also introduce the different variants we shall consider: the MSSM without µ-term; the MDGSSM, the DGNMSSM and the aligned DGNMSSM.In Section 3 we will study the predictions for all of these classes of models for the W boson mass, examining in particular the effects of a precise computation of the quantum effects.We present our conclusions in section 4.
2 Dirac Gaugino Models with Automatic Treelevel Alignment

Field content and interactions
We shall consider in this work the extension of the MSSM by a minimal matter content to allow Dirac Gaugino masses, as in [18,35].The additional superfields consist of three chiral multiplets, in adjoint representations of the SM gauge group factors (DG-adjoints): a singlet S, an SU (2) W triplet T a , and an SU (3) C octet O a .If we require gauge-coupling unification, even more states should be added to the model.For instance, for an (SU (3)) 3 Grand Unification, the minimal set of chiral multiplets includes also extra Higgs-like doublets R u,d as well as two pairs of vector-like righthanded electron E 1,2 in (1, 1) 1 and Ẽ 1,2 in (1, 1) −1 .We will not consider these states here.
In order to develop an intuition for the different interactions involved, it is helpful to consider a simple picture where the model descends from a supersymmetric theory in D dimensions.The different states can appear in different sectors: some live in the whole D-dimensional bulk, others are localised on four-dimensional hyper-surfaces (branes) at points of the extra dimensions of coordinates The corresponding Lagrangian can be written as where we have not explicitly written the metric factors.The four-dimensional theory arises after a truncation keeping only the compactification zero modes: boundaries . (2.4) A tree-level alignment in the Higgs sector appears in a class of models where the bulk theory leads to a four-dimensional Lagrangian with interactions governed by an N = 2 extended SUSY.In particular, the SM gauge fields and the DG-adjoint fields arise as N = 2 vector supermultiplets, and the two Higgs chiral superfields H d and H d form an N = 2 hypermultiplet, interacting through the superpotential where T ≡1 2 σ a T a , and the dot product is defined as The N = 2 SUSY has a global SU (2) R R-symmetry that rotates between the generators of the two N = 1 supercharges.The scalar components S, T a of S and T a , respectively, are singlets of SU (2) R .This R-symmetry rotates then between the auxiliary fields F a Σ of the adjoint superfields Σ a ∈ {S, T a } and the auxiliary component D a of the corresponding chiral gauge superfields W a i α for U (1) Y and SU (2) W .This implies that form a triplet of SU (2) R .As a consequence, in order that the interactions (2.5) of S and T a with the two Higgs doublets preserve SU (2) R , the couplings λ S and λ T must be related through , and to the couplings g Y and g 2 of the U (1) Y and SU (2) W gauge groups, respectively.Below the scale where the N = 2 SUSY is broken to N = 1, these relationships are spoiled by a small amount through renormalisation group running, so in numerical evaluations we must treat the couplings λ S and λ T as independent parameters.In addition to the SU (2) R R-symmetry which is broken in N = 1 (chiral) sectors, there is a global U (1) R R-symmetry under which the superspace coordinates θ α carry a −1 charge.The U (1) R charges of the H u and H d superfields are R Hu and R H d , respectively.They are arbitrary but subject to the constraint R Hu +R H d = 2.The DGadjoint superfields S, T a , and O a are R-neutral.Below, we shall classify the different N = 1 interactions following whether they preserve or break the U (1) R symmetry.
The boundary Lagrangian can be split into different contributions: (2.9) Here, we denote by L bulk localised kinetic and interaction terms already present in the bulk theory L 4d bulk but appearing with relative coefficients that violate N = 2 supersymmetry.Such terms can a priori be present because the boundary theory preserves only N = 1 SUSY, thus the coefficients of these terms are less constrained.Here, for simplicity, we assume such terms to vanish at tree level, to be only generated by quantum loops after supersymmetry breaking, and will therefore be accounted for in our analysis, at least in part, through the radiative corrections.
Also in (2.9), we have the usual MSSM Yukawa superpotential W Y ukawa with the couplings responsible for the quark and lepton masses: which arises on the brane where the matter field supermultiplets are localised.
In this work, we consider a typical scale for the soft terms, for example squarks or gaugino masses, to be in the phenomenologically interesting range m sof t ∼ TeV.If we denote by Λ a higher scale, for instance related to supersymmetry breaking messenger mass scale or to the Planck scale, then we can consider the relative strength of the diverse SUSY-breaking terms as an expansion in powers of m sof t Λ .We will assume that SUSY-breaking terms in the gauge sector preserve the U (1) R R-symmetry, giving rise to Dirac gaugino masses, while Majorana masses might be generated only by higherorder interaction terms, therefore suppressed by additional powers of hidden-sector couplings and/or where Λ could be the Planck scale (for gravity-induced effects).The effective superpotential for the Dirac gaugino masses reads where O ≡ 1 2 λ a O a , and W i α are the chiral gauge-strength superfields.Finally, the superpotential W N R contains terms that break explicitly the U (1) R R-symmetry: (2.12) The soft SUSY-breaking Lagrangian ∆L sof t can in turn be split in two parts.The first contains the scalar mass and interaction terms that preserve the U (1) R R-symmetry: The second part of ∆L sof t contains the scalar mass and interaction terms that break the U (1) R R-symmetry: as well as the Majorana mass terms M i (with i = 1, 2, 3) for the gauginos.In general, the mechanisms that break R-symmetry and SUSY could be independent of each other, hence in eqs.(2.13) and (2.14) we refrained from defining the soft SUSY-breaking trilinear couplings as proportional to the corresponding superpotential couplings.In this work we assume that the soft SUSY-breaking Higgs-sfermion-sfermion interactions in the second line of eq.(2.14) are suppressed with respect to the R-conserving sfermion mass terms in eq.(2.13).This can be realised in our D-dimensional picture if the quark and lepton superfields are localised on a brane that differs from the one where the breaking of the R-symmetry takes place.Since our study focuses on the electroweak sector of Dirac Gaugino models, we assume for simplicity that the scalar octet O a is heavy and can be integrated out of the theory.To insulate the singlet sector from threshold corrections involving the heavy octet, we also neglect the singlet-octet interaction term proportional to λ SO in the R-violating part of the superpotential, eq.(2.12), as well as the analogous term proportional to T SO in the R-conserving part of the soft SUSY-breaking Lagrangian, eq.(2.13).Similarly, since they cannot appear in some of our scenarios, for simplicity in the following we shall also neglect λ ST , T ST and the tadpole terms t S , ξ S .

The electroweak scalar sector and alignment
We can now discuss the neutral scalar sector of this class of models.The vacuum expectation values 2 (vevs) of the neutral components of the doublets H d and H u are related by v 2 u + v 2 d = v 2 , where v 246 GeV is the electroweak scale, and we define tan β = v u /v d .The neutral singlet and triplet scalars S and T 0 obtain vevs S = v S / √ 2 and T 0 = v T / √ 2, respectively.These lead to effective µ and B µ parameters: The vevs v S and v T are then determined as a solution for the coupled cubic equations: where are effective mass-squared parameters (at zero expectation value) for the real components S R and T 0 R of the neutral singlet and triplet scalars (the analogous masses for the imaginary components S I and T 0 T −B T , respectively).We know that v T must be small -namely, less than a few GeV -to avoid an overlarge tree-level ∆ρ, so to a good approximation we can set v T = 0 in the vacuum minimisation equation for v S , eq. (2.17), and decouple it from the one for v T , eq. (2.18); this would allow the cubic equation for v S to be solved using standard techniques.However, the current state of technology for the computation of loop corrections assumes that we take expectation values as being valid for the true minimum of the full quantum-corrected potential, so in our numerical studies we must take them as inputs.Especially for v T this can lead to complications; see [91] for a recent discussion of this issue.
To discuss the alignment in the Higgs sector, it is now convenient to introduce the so-called Higgs basis for the two doublets, where we defined for convenience two doublets with positive hypercharge, Φ j u ≡ H j u and Φ j d ≡ − ij H * j d .In the Higgs basis the two doublets can be decomposed as i.e., only the neutral component of Φ 1 has a non-zero vev, and the would-be-Goldstone bosons, G ± and G 0 , all lie in Φ 1 .In general, the neutral CP-even fields h and H mix with the neutral CP-even components of the singlet and the triplet.In the basis {h, H, S R , T 0 R }, the tree-level mass matrix reads where (2.30) Exact alignment in the Higgs sector is obtained when one of the eigenstates of the CP-even mass matrix -in this work, we take it to be the lightest one -is aligned in field space with the direction of the SM Higgs vev, and thus has SM-like couplings to gauge bosons and matter fermions.This is equivalent to requiring that h itself be an eigenstate of M 2 , or in other words that M 2 1j = 0 with j = 2, 3, 4. If we make the reasonable assumption that the triplet is heavy, then this can be relaxed to just j = 2, 3.In addition, we will also refer in this work to cases where the singlet can be light without being potentially ruled out by direct searches.In this case we will require the supplementary condition that M 2 23 = 0.
We start our discussion by focusing on the alignment between the two doublets.The use in eq.(2.23) of the N = 2 condition for the singlet and triplet superpotential couplings, see eq. (2.8), implies ∆ h = 0 and M h = M Z , i.e. alignment is automatically realised at the tree level in this class of models, and the tree-level mass of the SM-like Higgs boson is independent of tan β but well below the value observed at the LHC.It is however well known that, in SUSY models, the radiative corrections to the Higgs mass matrix play a crucial role in lifting the prediction for the mass of the SM-like Higgs boson up to the observed value.Moreover, the radiative corrections to the condition in eq.(2.8) for the superpotential couplings of the adjoint superfields can become relevant if the scale M N =2 where the N = 2 SUSY is broken to N = 1 is much larger than the scale where the Higgs mass matrix is computed.All of these corrections inevitably affect also the condition for alignment in the Higgs sector.As was discussed in ref. [66], in DG models the element that mixes the two doublets in the loop-corrected mass matrix can be recast as ) where M 2  11 contains the dominant one-loop contribution from top and stops.The latter consists in a term enhanced by , where y t is the top Yukawa coupling and M SUSY denotes for simplicity a common soft SUSY-breaking mass parameter for the stops.The second term in eq.(2.31) accounts for the deviation of the superpotential couplings from the SU (2) R condition in eq.(2.8), and the ellipses denote one-loop top/stop contributions that are suppressed by small ratios of parameters, one-loop contributions that involve couplings other than y t , and higher-loop contributions.Close to alignment, the loop-corrected mass-matrix element M 2  11 can be empirically identified with the observed mass of the SM-like Higgs boson, M 2 h ≈ 2 M 2 Z .Therefore, eq. ( 2.31) shows that the radiative corrections included in M 2  11 tend to destroy the tree-level alignment in the Higgs sector of DG models.However, when M N =2 is large the evolution of λ S and λ T down to the scale where the Higgs mass matrix is computed makes the second term in eq.(2.31) negative, and partially compensates for the misalignment induced by the top/stop contributions.
It is instructive to compare the condition for doublet alignment in DG models with the analogous conditions in the MSSM and in the NMSSM.In the case of the MSSM, discussed e.g. in refs.[6,92], one finds where A t ≡ T 33 u /y t is the soft SUSY-breaking Higgs-stop-stop interaction parameter, and again the ellipses denote sub-dominant terms.It appears that, in the MSSM, the alignment condition M 2 12 = 0 can be realised radiatively when a large value of tan β suppresses the first term in eq.(2.32), while the parameters M SUSY , A t and µ combine in such a way that the second term is large and negative.In contrast, in DG models the contributions to M 2  12 analogous to the second term in eq.(2.32) are suppressed by the assumption that A t M SUSY , see the comments after eq.(2.14), thus doublet alignment cannot be realised in this way.We remark however that, even with the N = 2 condition for the superpotential couplings, in DG models M 2  12 is smaller by a factor between 2 and 3 -depending on tan β, which we assume to be greater than 1with respect to the case of the MSSM with small A t .
In the case of the NMSSM, discussed e.g. in refs.[7,9], the mixing between h and H is given by where μ ≡ µ+λ S v S / √ 2. Comparing with the case of the MSSM, eq.(2.32), we see that the condition M 2 12 = 0 can be realised even in the absence of a large contribution from the terms proportional to µ eff A t /M 2 SUSY , as long as the singlet-doublet superpotential coupling takes values in the range λ 2 S ≈ (3 − 4) M 2 Z /v 2 , where the exact numerical coefficient depends on the value of tan β.As first pointed out in ref. [7], this condition singles out the region of the NMSSM parameter space where λ S ≈ 0.7 ± 0.05, a much larger value than would be implied by the SU (2) R condition in DG models.
To summarise, the SU (2) R R-symmetry implies exact alignment at the tree level in the Higgs-doublet sector of the DG models, but the alignment is partially spoiled by the radiative corrections that are necessary to obtain a realistic value for the SMlike mass.Alignment in the MSSM can be realised only through radiative corrections, for large tan β and for specific choices of the parameters in the stop sector.Finally, doublet alignment in the NMSSM can be realised even without the help of radiative corrections for an appropriate choice of λ S , which -differently from the DG case with SU (2) R R-symmetry -is treated as a free parameter.
The second condition for Higgs alignment in DG models is M 2 13 = 0, i.e. vanishing mixing between h and S R .Including the dominant contributions from stop loops, we find: where the tree-level mixing term ∆ hs is given in eq.(2.24), µ eff is given in eq.(2.15), and Q is the renormalisation scale at which the parameters entering ∆ hs are expressed.We assumed again a common soft SUSY-breaking mass term M SUSY for the stops, and we neglected terms suppressed by powers of m 2 t /M 2 SUSY .The various terms that contribute to ∆ hs arise from different sectors of the D-dimensional picture discussed earlier in this section: namely, µ and λ S enter the bulk superpotential in eq.(2.5); m 1D enters the R-conserving boundary superpotential in eq.(2.11); M S and κ enter the R-violating boundary superpotential in eq.(2.12); T S enters the R-violating SUSYbreaking Lagrangian in eq.(2.14).Therefore, even if we assume the N = 2 SUSY relation of eq.(2.8) between λ S and g Y , a vanishing M 2  13 can only result from an accidental cancellation between unrelated terms.We also note that, in contrast to the case of M 2  12 , the radiative correction to M 2 13 is not enhanced by tan β with respect to the tree-level part.Thus, its qualitative impact on our discussion of the alignment conditions is limited, as long as the scale Q is not too far from M SUSY .
The minimum conditions of the scalar potential can be exploited to express the mass parameters for the doublets and the singlet in terms of the other Lagrangian parameters and of the vevs v d , v u and v S .In particular, we obtain a relation between the mass parameter m2 SR for the real component of the singlet, see eq. (2.19), and the matrix element M 2  13 given in eq.(2.34): Eq. (2.35) above shows that the condition of vanishing mixing between h and S R , however realised, carries implications for the mass of the singlet.The diagonal element for the singlet in the scalar mass matrix is where we applied the same approximations as in eq.(2.34) for the one-loop correction in the second line.
We now discuss the simplest case in which the global U (1) R R-symmetry is preserved in the superpotential but broken by soft SUSY-breaking terms, in which case we can set κ to zero.Since this implies a vanishing quartic self-coupling for the singlet, the stability of the scalar potential requires that we also assume T κ = 0, even if the trilinear self-coupling of the singlet resides in the R-conserving part of the soft SUSYbreaking Lagrangian.In this scenario, which we shall refer to as the aligned MDGSSM, eq.(2.35) shows that the alignment condition M 2 13 ≈ 0 requires m2 SR ≈ 0 or v S v.The first of these two options implies that the CP-even mass eigenstate that is mostly singlet is relatively light: setting m2 SR = 0 and κ = T κ = 0 in eq.(2.36), and neglecting the small effect of the one-loop correction, we find that the value λ S ≈ 0.7 favored by the alignment condition for the Higgs doublets in the NMSSM, see eq. (2.34), leads to M 2  33 ≈ (122 GeV) 2 , whereas the value λ S ≈ 0.25 implied in our Dirac-gaugino model by the N = 2 SUSY relation of eq.(2.8) leads to M 2  33 ≈ (44 GeV) 2 .We remark that the mixing between S R and the heavier, non-SM-like scalar H, which is controlled by M 2  23 , is suppressed when 33 , and would in any case lower the mass of the singlet-like eigenstate.
The definition in eq.(2.19) shows that even the vanishing of m2 SR requires a cancellation between terms that arise from different sectors of our D-dimensional construction: m 1D from the R-conserving boundary superpotential in eq.(2.11), M S from the R-violating boundary superpotential in eq.(2.12), m 2 S and B S from the R-conserving soft SUSY-breaking Lagrangian in eq.(2.13).If such cancellation is not realised, the alternative requirement for alignment implied by eq.(2.35) when κ = 0 is that v S v.This is not problematic as long as a suitable higgsino mass is provided by the µ term in the bulk superpotential, see eq. (2.5).
Finally, the third condition for Higgs alignment in DG models is M 2 14 = 0, i.e. vanishing mixing between h and T 0 R .The formulas for the relevant mass-matrix elements and for the minimum condition, including the dominant one-loop corrections from top and stop loops, are similar to eqs.(2.34)-(2.36),with the obvious singlet-to-triplet replacements but without terms analogous to those controlled by κ in the singlet case: (2.39) The discussion of the constraints on the triplet mass induced by the condition of doublet-triplet alignment follows the lines of the discussion of singlet-triplet alignment for κ = 0, with the important difference that the condition v T v must in any case be satisfied to avoid an excessive contribution to ∆ρ.As a consequence, it is not necessary to require m2 T R ≈ 0 to obtain approximate alignment.Nevertheless, we remark that the condition of exact alignment M 2 14 = 0 would imply M 2 44 ≈ M 2 W through eqs.(2.37)-(2.39).

The electroweak fermion sector
We now outline the mass spectrum of the electroweak fermions, which will be relevant for our discussion examining the W boson mass.The neutralino mass matrix, in the basis S, B, T 0 , W 0 , H0 d , H0 u reads: , are given by and do not depend on κ, but we have written for completeness the Majorana gaugino masses M 1,2 for the bino and wino respectively.

Scenarios
We have described the general features of minimal Dirac Gaugino models, and explained how some values of the couplings allow Higgs alignment to automatically occur at tree-level.Below, we shall consider the following different scenarios corresponding to specific choices of the model parameters: • General MDGSSM In the general MDGSSM, the only source of R-symmetry violation comes from a small B µ term, which is a radiatively stable condition (the renormalisation group running will not generate other R-symmetry violating terms from a B µ term).This excludes all of W N R ; the only supersymmetric parameters beyond those of the MSSM retained are λ S , λ T , but those are allowed to have any value.Supersoft masses [13] are allowed, and soft supersymmetry-breaking masses are allowed, but squark/sfermion trilinears and T S , T T are not.On the other hand, trilinears involving the adjoint scalars may be allowed (so T κ , T ST , T SO ) but may be argued to be small in typical models.The couplings λ S , λ T enhance the Higgs mass and W mass at the same time; while large values of λ T were previously deemed problematic for the ρ parameter, they are now a virtue.
• MSSM without µ term (RIP) This model, described in [14], is identical to the general MDGSSM execpt that the µ-term is set to zero (although a small v S generates a tiny effective µ).In the MSSM this would yield massless higgsinos, but here the higgsinos obtain a mass through λ T which causes them to mix with the triplet fermion.Thus the µ problem of the MSSM is solved, at the expense of an upper bound on the chargino masses, which, as we shall see, leads to the model being ruled out.
• Aligned MDGSSM By taking λ S , λ T to their their N = 2 values given in Eq. (2.8) in the general MDGSSM we guarantee alignment with the heavy Higgs at tree level.If we further impose that ∆ hs = 0, as described above, then the singlet also has negligible mixing with the SM-like Higgs.We shall refer to this scenario as the aligned MDGSSM.
• DGNMSSM If we instead allow R-symmetry violation in the superpotential and the associated soft-breaking trilinears -in particular for κ, T κ and T S -we can generate µ and B µ terms through a substantial expectation value for the singlet.The model thus resembles the NMSSM, especially if we set the µ, B µ terms (and M S , M T ) to zero; this was proposed in [18].We shall therefore refer to this scenario as the DGNMSSM.

• Aligned DGNMSSM
We can achieve aligment in the DGNMSSM by setting λ S , λ T to their N = 2 values and taking ∆ hs = 0. We shall refer to this scenario as the aligned DGNMSSM; in this work, when we consider the W boson mass, we shall also enforce ∆ Hs = 0 which guarantees that the singlet couplings to SM fields are small, rendering light singlets safe from collider searches.

W mass in Dirac Gaugino models
Dirac gaugino models offer two methods of explaining an enhancement of the W boson mass with respect to the SM: either quantum corrections or a tree-level expectation value for the triplet scalar.In the case of quantum corrections, there are new contributions to the W mass compared to the MSSM, again coming from interactions related to the adjoint triplet; in particular the coupling λ T .
Recall the definition of the ρ parameter is M 2 W ≡ ρc 2 W M 2 Z ; through the presence of the adjoint triplet, DG models contain a tree-level modification to this relation compared to the SM: One could consider a modification to c W instead of ρ as an explanation of an enhanced W mass, but this is discounted based on electroweak precision tests (see e.g.[93]); on the other hand, a triplet expectation value is one of the most generic and acceptable ways of enhancing the W mass at tree level.Naively we could then take the observed value of M Z and the standard value of c W and infer the value of ∆ρ to obtain a given value of M W .However, in the SM and in any BSM theory it is necessary to take certain electroweak observables as input, and in our setup we will take the conventional choice of the Z mass, G F and α.When we modifiy ρ this gives both a modification to M W and a small modification to sin 2 θ W .At tree-level this is so if we want to obtain M W = 80.424 MeV we need ∆M 2 W = 11 (GeV) 2 .Taking s 2 W = 0.23121 this gives ∆ CDF ρ tree =0.0012 Interpreted as a tree-level expectation value for the triplet, this yields which was previously at the upper bound of what was acceptable.
In the minimal Dirac gaugino model, we have from equation (2.39) For a small v T , if we do not tune ∆ ht ≈ 0, then triplets must be heavy.However, we also need heavy winos to evade collider bounds and this implies large m D2 : the connection between electroweakino masses and the expectation value of the triplet is a novel feature of this class of models.In [64] the conclusion was that for winos above 700 GeV there were essentially no constraints on the higgsinos beyond LEP.If the contribution from m D2 dominates ∆ ht this implies which is a natural scale for supersymmetric scalars.Of course, we can have lighter winos provided that the neutralino is not too light, above around 200 to 300 GeV [64].However, the model also contains ample room for quantum corrections to also enhance the W mass.In the following, we shall investigate this for the different scenarios described in section 2.4.

Numerical setup
In order to accurately compute the quantum corrections to the W mass, we use a new EFT approach, closely related to that of [87], implemented in the spectrum-generatorgenerator SARAH.We use the expression: where (M 2 W ) SM is the full two-loop W mass in the SM, as computed in [75], depending on the pole masses of the Higgs boson, top quark, α s and ∆α (5) had .We use the interpolating function from that paper.When we use the average values of the Higgs mass of 125.09GeV and the top quark mass of 172.89 GeV this function gives us M W = 80.354 GeV, just 2 MeV higher than the current world average for the W boson mass in the SM.The expressions computed in the square brackets are the differences between the high-energy theory (HET) and the SM: ∆α are now the gauge threshold corrections between the HET and the SM for the electromagnetic gauge coupling divided by α (so they do not now depend on ∆α (5) had ).In SARAH, we compute the expression in square brackets at the matching scale, which is the mass of the heavy particles.We are therefore ignoring the running from that scale down to the electroweak scale, which is of controllable size, but will nevertheless be included in a future development of the code.One could argue that we should instead perform the matching at the electroweak scale, but then all of the loop functions will contain large logarithms and there can be larger, spurious, running of the couplings of the high-energy theory which can spoil the results.
In order to be a strict one-loop matching between the HET and the SM, the weak mixing angle in the above must be the M S or DR value at the matching scale Q: To extract this, we match the calculations of the Z-boson mass, α(Q) -and we also compute the decay of the muon at the matching scale.In practice, this means that we extract the couplings in the SM at the top mass scale without including the effects of new physics, then run them up to the matching scale.The threshold corrections to α between the two theories are simple to compute since it is unbroken and yield ∆α; the coupings g Y (Q), g 2 (Q) in the high-energy theory are chosen to solve the equation which is done iteratively by progressively running up and down and updating at each step, along with all the other quantities in the high-energy theory.The value for (c 2 W s 2 W ) SM includes a compensatory term for corrections from MS to DR if needed.In the SARAH model file, the couplings λ S , λ T , T S , T T are defined differently to the above: we have lam Hence in our plots and benchmark points we list the values in terms of −λ S , √ 2λ T , −T S , which makes the correspondence with the numerical codes exact.

MSSM without µ term
The "MSSM without µ term" proposed in [14] was an intriguing solution to the µ problem.It was however challenged by the requirement of having a high enough Higgs mass, chargino mass and not too large ρ; indeed the lack of intersection of points satisfying the latter two was demonstrated in [35].It might therefore be tempting to revisit this model in light of the new data about the W mass.However, we shall demonstrate here that it is conclusively ruled out.
Putting aside the Higgs mass constraint, the see-saw effect on the charginos is a problem.LEP put a lower limit on the mass of the lightest chargino of 94 GeV [94][95][96].In that model the chargino mass is It is known that it is possible to fulfil the LEP bound by a careful choice of λ T and m D2 : a large value of λ T as well as m D2 around 107 GeV is needed to maximize the mass of the lightest chargino.It would then be made of a higgsino-wino mixture with two charginos that are light and one (wino-like) somewhat heavier.Unfortunately, subsequent LHC searches are especially sensitive to winos up to about 800 GeV, see [64].
It might be possible to evade this constraint if the light wino and neutralino are close enough in mass so that decays such as χ± → χ0 + W ± are not possible.Most likely this is difficult or impossible to achieve, but without a detailed investigation we cannot exclude the possibility that some region of parameter space might evade direct LHC searches; we can only apply the LEP constraint as a hard lower bound on the chargino mass.
In [35] it was demonstrated that that the corrections to ∆ρ were correlated with the lightest chargino mass; in order to evade the LEP bound, λ T has to be large and this drives large ∆ρ.Here we can give a striking confirmation of this observation by plotting the W mass against the lightest chargino mass for a sample of O(70000) spectra generated using SARAH.We fix the octet scalar, and all squark and slepton soft masses via m 2 O = m 2 q = m 2 l = 10 TeV 2 , and fix the gluino mass to 3 TeV; this ensures that they are beyond all current and near-future bounds.Then we scan over the ranges: We use a Markov Chain Monte-Carlo (MCMC) algorithm to generate points (this helps to obtain points with non-tachyonic spectra with Higgs mass close to the observed value compared to a random scan).There is not intended to be a genuine statistical interpretation of the distribution of the points, but the overall envelope should show where valid sets of parameters exist.The results are shown in figure 1, where the LEP constraint is shown as a vertical green band.It can be clearly seen that it is not possible to both satisfy the LEP constraint and have an acceptable value for the W boson mass; from the W boson mass alone we would predict a chargino of mass below 65 GeV.

W mass in the MDGSSM
In the minimal Dirac Gaugino extension of the Standard Model (MDGSSM) the most typical scenario is to assume that R-symmetry is broken only via a B µ term, i.e. the Higgs sector is special.It is identical to the previous model except that we allow a µ-term.This, however, makes all the difference: now the higgsino mass is not bounded from below, not requiring a large mixing with the winos; and further the enhancement to the Higgs mass is under control.
We perform a new MCMC scan with parameters allowed to vary within the ranges: We fix the octet scalar, and all squark and slepton soft masses via m 2 O = m 2 q = m 2 l = 10 TeV 2 , and fix the gluino mass to 3 TeV.We choose a likelihood function to be a product of a gaussian in the Higgs mass with mean 125 GeV and standard deviation 3 GeV, a gaussian in the W mass with mean 80.413 GeV and standard deviation 20 MeV, and a sigmoid on the constraints (given as the maximum ratio of predicted crosssection to observed, across all channels) from HiggsBounds5 [97-100], which strongly suppresses the likelihood when the observed cross-section ratio is greater than one, but is otherwise close to unity.This choice of likelihood function is merely a device to select desirable points, and the distribution of the points is not meant to have a statistical interpretation in terms of their Bayesian likelihood; in particular, the theory uncertainty on the Higgs mass does not have a statistical interpretation, and a window of 3 GeV is a conservative estimate of the average error, since we use the latest twoloop corrections in the generalised effective potential and gaugeless limit [60,[101][102][103] with pole-mass matching onto the SM [104] (see [105] for a recent review).Such a conservative window of 3 GeV is employed because only a relatively small proportion of points actually generate a spectrum, and in principle a two-stage procedure prefiltering points along those suggested in [64,106] would probably be more efficientor in addition the ability to invert the vacuum minimisation relations and compute v S , v T instead of treating them as inputs, but this is not yet automatically possible in the code in a way that would correctly incorporate the loop corrections to the Higgs masses [91].
We show plots in figure 2 for M W against v T (left plot), and then for M W against λ T (right plot), to show the points that benefit from large quantum corrections as the means of enhancing the W mass: the tree-level expectation just from modifying v T is shown as a solid red curve on the left plot.We also show as a dashed red curve the value of M W that would be obtained by insisting that the shift in ∆ρ only modifies M W without changing sin θ W (if a different method of matching onto the SM parameters were used, for example).
All points shown satisfy all Higgs bounds; have a charged Higgs heavier than 600 GeV (so are safe from B → sγ constraints [1]); have charginos heavier than the LEP limit and winos heavier than 600 GeV.These baseline selections are marked as blue points in the plots; there are about 36000, of which about 10000 have |v T | < 1 GeV.Points shown in yellow further have charginos heavier than 200 GeV, m D2 > 700 GeV; while those marked in green have charginos heavier than 250 GeV and m D2 > 800 GeV (about 1600 points in our sample survive these cuts).The green points are thus almost certainly guaranteed to be safe from current collider bounds (although they may yet be probed in future).The requirement of heavy charged Higgs scalars sets the MSSM-like neutral and pseudoscalar masses to be heavy, and essentially guarantees the safety of all selected points from constraints on the couplings of the SM-like Higgs.Nevertheless, we also filtered the green points with constraints from HiggsSignals [107,108].The  different categories of points show the expected wider range of enhancements to the W mass as the charginos become lighter.We provide a selection of benchmark points, with the input parameters and crucial data, in table 1.
The first clear observation is that the quantum corrections to the W mass are at least as important as the tree-level contribution from the expectation value v T , and the generic contribution to the W mass is positive with no points below the red curve.The asymmetry of the plot with v T is due to the fact that we only take positive Dirac gaugino masses.It is important to note that in the red curves we take the SM value of the W boson mass to be 80.352 GeV, whereas the fitting function of [75] as employed in SARAH gives a value of 80.354 GeV for the central values; and gives 80.356 for a Higgs mass at the lower bound of our permitted range of 122 GeV.However, it is clear that the quantum corrections in our sample are generally more important than ∆ρ tree .Indeed, with the parameter ranges we have chosen, we are not within the range of masses required for decoupling of the quantum corrections (we have checked that the quantum corrections to M W smoothly drop to near zero as the masses of all particles are raised to about 2 TeV or higher).
The second observation is that there is no clear correlation between the W mass and the parameter λ T .Due to our requirement of a large wino mass, the selected points have light neutralinos/charginos of mixed bino/higgsino type.The mass splitings among the higgsinos -and thus the contribution to the W mass -can be driven large by λ S and λ T .These couplings also enhance the Higgs mass at tree-level, so large values are favoured in the scans because we fixed the stop masses at √ 10 TeV.In this class of model there is therefore no particular preference for one or the other coupling.
Since the selected points generally have a mixed bino/higgsino LSP, they have good dark matter candidates, but it is expected that the relic density should be underdense.A detailed investigation of the dark matter-collider complementarity for scenarios satisfying the latest W mass data along the lines of [64] would be an interesting subject for future work, provided that the latest LHC analyses can be recast.As mentioned above, the possibility of a large λ T (and the presence of the singlet scalar/fermion) distinguishes higgsinos in this scenario from those in the MSSM.It is clear that this class of models provides a very natural explanation for an enhancement to the W mass compared to the Standard Model.

W mass in the aligned MDGSSM
In the aligned MDGSSM (where the only source of R-symmetry breaking is the B µ term, and we take T T = 0), we choose the parameters to induce alignment at treelevel in the MDGSSM (so 2 and m DY = c 2β µ eff ) but taking T S = 0 (which means the mixing of the singlet with the heavy Higgs cannot vanish unless it is heavy).To study this scenario we perform a scan with the same strategy as before except that now, since λ S , λ T are fixed, it is necessary to vary the masses of the stops/sbottoms to allow us to find the observed value of the Higgs mass; there is also therefore a preference for models with larger tan β since the tree-level contributions to the Higgs mass at low tan β are not sufficient.We therefore use a common mass for the third generation squarks M SUSY ; we set m  We mostly find points with small µ/m DY and very little enhancement to the W mass because the electroweakinos tend to be light.We give benchmark points in table 2 and plots in figure 3 which demonstrate the lack of enhancement and scarcity of points (790 survived from a scan for one million).

W mass in the aligned DGNMSSM
We turn now to the case of the aligned DGNMSSM described in section 2.4.We set the couplings λ S , λ T to their N = 2 values and then choose m DY and T S to make ∆ hs and ∆ Hs vanish.This leads to .16)This has an interesting consequence because in this model  Then We perform a scan using the same strategy as the previous sections; we use a common mass for the third generation squarks M SUSY ; we set m 2 Q,33 = m 2 U,33 = m 2 D,33 = M 2 SUSY (the other squarks and sleptons we retain fixed at √ 10 TeV).Then we scan with the parameter ranges, using the same likelihood function as before: In our scans we impose that all Higgs searches are satisfied using HiggsBounds and HiggsSignals.We show the results for the W boson mass in figure 4 where the points have the same colour coding as in the previous sections.It is apparent that in this model it is complicated to enhance the W boson mass.This is because we have only a small quantum effect from λ S , λ T , but also because the lightest neutralinos are typically rather light: since µ eff = g Y v S /2 we need a large v S 500 GeV (or 1 TeV for our more stringent points) to have heavy enough higgsinos.Then we need κ negative and not too small to avoid a too-small pseudoscalar/charged Higgs mass (if we neglect κ then M A is bounded by g Y v S , so M A > 600 GeV requires v S 1700 GeV, at the limit of our search range).So this implies that the singlino is generally heavy compared to m DY :  In figure 4 we show the scan results with the same colour coding as in previous sections.It is clear that in this scenario, models which can explain a large W mass are driven by a larger v T with some modest quantum corrections enhancing the mass by O(10) MeV; there is very little spread due to the lack of variation in λ T .However, it is difficult to find points with large enough v T that satisfy other bounds.
Another feature of the selected points is that in almost all cases v T > 0; considering T T = M T = 0 in equation (2.18) means that in order for m2 T R > 0 (so that, at least, the pseudoscalar triplet should be non-tachyonic, since we take B T = 0) and v T < 0 we would need µ eff > |m D2 c 2β |.But we need large tan β to obtain the correct Higgs mass, and therefore g Y v S /2 m D2 ; for our selected points we require a minimum of m D2 > 600 GeV, and so v S would again be beyond our search range.
Since we are interested here in alignment, we may have a light singlet scalar without falling foul of either light Higgs or heavy Higgs searches.In this limit we have T κ of opposite sign to v S allows the singlet to be made light while making M A arbitrarily heavy.In the scan, we do not impose any likelihood bias to search for points with a light singlet, but we show benchmark points passing all constraints which have light singlet masses in table 3.They show a W mass consistent with the SM prediction.

W mass in the general DGNMSSM
Finally we consider the general DGNMSSM, where we allow the values of λ S , λ T to vary and do not fix the values of m DY or T S to require alignment but scan over them.This means that we will not focus on light singlet (or doublet) scalars.Similar to the aligned case, there is still a see-saw effect on the lightest neutralino mass due to the non-zero singlino mass, which can drive down the quantum corrections to the W boson, but a large |λ T | can compensate for this and also help enhance the SM-like Higgs mass.We perform a scan using the strategy as in sections 3. We give plots in figure 5 with the same colour coding as in the previous sections; the difference in the distribution to the previous examples is rather striking.It is clear that in this scenario a large negative λ T and a positive v T is favoured; this gives a tree-level enhancement to the Higgs mass and a loop-level enhancement to the W -boson mass.
The singlino component will mix less with the higgsinos than in the MDGSSM because of the √ 2κv S singlino mass, and thus the effect of λ S on the W mass is diminished.The asymmetry in the signs of λ T and v T can be explained by the fact that we only take positive Dirac gaugino masses in the scans.

Conclusions
We have shown that an aligned Dirac Gaugino NMSSM is possible and compatible with current collider constraints; it can even lead to relatively light singlet scalars that may be of interest to future searches (although would be rather difficult to find directly  as they are difficult to produce).Such a model favours a W boson mass compatible with or just above the SM prediction.We also showed how two different Dirac Gaugino scenarios can easily be compatible with an enhanced W boson mass, including a precise computation of the quantum corrections for the first time, which are now incorporated automatically in the package SARAH.We also used this computation to add more nails to the coffin of the "MSSM without µ term." We have been conservative in our application of collider constraints and concluded that the MDGSSM models would typically contain underdense dark matter densities.However, it would be interesting to examine the issue of dark matter and collider constraints again in all of these classes of models when all the latest searches for electroweakinos have been recast; we have provided ample benchmark points for this purpose.In the DGNMSSM or its aligned version, if we impose strict R-parity or have a heavy gravitino (by no means entirely obvious assumptions), it may be that we require a Higgs funnel to obtain the correct relic density, which would require a sophisticated search strategy to find allowed parameter ranges, along e.g. the lines of [106].However, it is also likely that a light singlino in the aligned DGNMSSM could fulfil the role of the Higgs funnel.We leave these questions to future work.

Figure 1 :
Figure 1: W boson mass vs lightest chargino mass for points generated in the 'MSSM without µ-term.'The horizontal bands show the SM, Tevatron and world average masses for the W boson mass; the vertical green band shows the LEP constraint on the lightest chargino.

Figure 2 :
Figure 2: Left: W boson mass vs v T in the MDGSSM, with the red curve showing the tree-level prediction.Right: W boson mass vs √ 2λ T in the MDGSSM.Colours of the points are described in the text, with those obeying the strictest cuts shown in green.The colourful horizontal bands show the SM range in light orange; the Tevatron+LEP average in purple, and a conservative world average in green.

Figure 3 :
Figure 3: Points in the aligned MDGSSM.Left: M W against triplet expectation value.Right: M W against singlet-like scalar mass.Colour coding as for previous plots. .20)

Figure 4 :
Figure 4: W boson mass in the aligned DGNMSSM.Left: W boson mass against triplet expectation value.Right: W boson mass against singlet-like Higgs mass.Colour coding of points is described in the text.

Figure 5 :
Figure 5: W boson mass in the general DGNMSSM.Left: W boson mass against triplet expectation value.Right: W boson mass against λ T .Colour coding of points as in previous figures.

Table 1 :
Benchmark points for the MDGSSM.Input parameters are given above the double line, and masses of the most important particles below.

Table 2 :
Benchmark points with a light singlet in the Aligned MDGSSM

Table 3 :
Benchmark points for the "aligned DGNMSSM".Input parameters are given above the double line, and masses of the most important particles below.

Table 4 :
Benchmark points for the general DGNSSM.Input parameters are given above the double line, and masses of the most important particles below.