Higgs alignment from extended supersymmetry

We consider the effective type-II two-Higgs doublet model originating from Dirac gaugino models with extended supersymmetry in the gauge sector, which is automatically aligned in the simplest realisations. We show that raising the scale at which the extended supersymmetry is manifest and including quantum corrections actually improves the alignment. Using an effective field theory approach including new threshold corrections and two-loop RGEs, plus two-loop corrections to the Higgs mass in the low-energy theory, we study the implications from the Higgs mass and other experimental constraints on the scale of superpartners. We contrast the results of the minimal Dirac gaugino model, where alignment is automatic, with the hMSSM and the MRSSM, where it is not, also providing an hMSSM-inspired analysis for the new models.


Introduction
In the absence of signals of strongly-coupled particles at the LHC, it has become important to study the possibility of new particles that couple to standard model (SM) states only via couplings of electroweak strength. The bounds on such particles are still relatively weak but with much luminosity to arrive there is still a substantial parameter space to explore, and such theories perhaps represent now the best chance for discoveries. Among such theories, one that has received significant and now increasing attention is the Two Higgs Doublet Model (THDM); see e.g. [1][2][3][4] and references therein. It is important to ask the question: "does the Higgs sector just consist of one doublet?" because the answer will give profound information about nature. If there are indeed additional fundamental scalars that mix with the Higgs boson, then this dramatically worsens the Hierarchy problem and a e-mail: kbenakli@lpthe.jussieu.fr b e-mail: goodsell@lpthe.jussieu.fr c e-mail: swilliamson@lpthe.jussieu.fr would necessitate a rethinking of our ideas of naturalness. On the other hand, such sectors naturally appear in the context of supersymmetry (SUSY) and it is conceivable that a second Higgs doublet could be the harbinger of a full SUSY theory.
However, the measurements of the Higgs boson's couplings already place significant constraints on the amount of mixing that it can suffer. It is for this reason that there has been much interest in the idea of alignment in the Higgs sector, i.e. that the mass eigenstates align with the vacuum expectation value, because in this case the couplings would be exactly SM-like.
To quantify this, consider two Higgs doublets 1 , 2 which mix, and then rotate their neutral components as follows: where we shall throughout use the notation c β ≡ cos β, s β ≡ sin β, t β ≡ tan β.
In this basis, we can write the mass matrix as where the quantities Z 1 , Z 5 , Z 6 are functions of the quartic couplings and mixing angles only; we shall give explicit expressions for this relationship later, in equation (2.20). Clearly the mass eigenstates are onlyh,H if and this is the condition for alignment, because the fields align with the electroweak vacuum expectation value. On the other hand, if Z 6 = 0, we must make a further rotation which is conventionally parameterised by an angle α as where now h, H are the two mass eigenstates. We shall assume throughout that h is the lightest eigenstate. In terms of the masses of the physical bosons m h,H this gives . (1.4) In both the type-I and type-II THDM, there is a Higgs eigenstate that couples to the up-type quarks, and we define this eigenstate to be 2 . This means that the ratio of the h coupling to all up-type quarks compared to the SM Higgs' value is while the ratio of the coupling to vector bosons to the SM value is also determined entirely by the mixing (neglecting loop effects from the rest of the extended Higgs sector): However, there is a combined ATLAS + CMS bound [5] on the ratio of these: (1. 6) This is enough to constrain where the latter bound comes from the value t β = 7.3, and the bound is much more stringent for large or small t β . For m H somewhat above m h this is a rather weak constraint, only becoming relevant when the two states approach degeneracy. However, in the type-II THDM, there is another constraint from the ratio of the ratio of the neutral Higgs coupling to all down-type quarks compared to its SM value (1.8) and, since from the previous constraint we know that the denominator is nearly equal to one, we have −0.04 t β t β−α 0.2 (1.9) which in turn implies t β−α t β and so s β−α c β−α This leads to a sensible constraint; for example, for m H = 600 GeV and t β = 5 it leads to Z 6 0.2. So we see that either we should take the mass m H to be large, in which case we have decoupling, or we keep it light in order to possibly detect it at the LHC, in which case we need alignment without decoupling (see e.g. [6]). However, as we have seen this is non-trivial; as the LHC measurements become more precise, the constraints will tighten further, and it is in this spirit that it is important to consider models where the alignment is natural rather than ad hoc.
The problem for the different types of THDM is that alignment without decoupling is not generic when we choose the masses -or equivalently quartic couplings -from the bottom up. Hence it is logical to derive the couplings of the THDM from some higher-energy theory and look for cases where alignment arises naturally. For example, [7][8][9] proposed models which lead to a natural alignment condition, based on additional bosonic symmetries. Here, on the other hand, we shall show how alignment arises automatically in a class of supersymmetric models, in contrast to the MSSM or NMSSM [10], with the additional benefits of (greatly) increasing the naturalness of the model and being able to predict the scale of new superpartners. Moreover, we shall show that quantum corrections actually improve the alignment!
The class of models that we shall consider have a gauge sector which is enhanced to N = 2 supersymmetry at a (potentially high) scale M N =2 . This fits into the framework of Dirac gaugino models, which have been well-studied in, for example . In particular, the idea of N = 2 supersymmetry in the gauge sector only and the consequences for the Higgs sector were first explored in [16] and recently studied in [72,73]. In general, though, this was either taken to be at the same scale as the other superpartners [73], or only a rough estimate of the main contribution of the chiral sector was included [16], while we shall show that increasing M N =2 improves alignment and increases naturalness! In Sect. 2 we will describe our theory and how it leads to natural alignment at tree level. In Sect. 3 we will outline the effect of radiative corrections. In Sect. 4 we perform a precision study of the model using an EFT approach to obtain the parameters at low energies, give predictions for the scale of new physics from the value of the Higgs mass, and explore the consequences for alignment. In Sect. 5 we consider all of the relevant constraints on the model space, including the latest LHC search for decays to τ pairs, b → sγ searches and electroweak precision constraints, and show how this affects our model. In the appendix we give all of the one-loop threshold corrections for our model at the scale of supersymmetry. Finally, in Sect. 6 we briefly consider the case of the MRSSM. To endow gauginos with a Dirac mass, at a minimum we need to add chiral fermions in the adjoint representation of each gauge group, which means adding adjoint chiral superfields: a singlet S, an SU (2) triplet T, and an SU (3) octet O. If we add just these fields, then we have the simplest Diracgaugino extension of the MSSM whose Higgs sector has been well studied [21,26,32,44,69]. However, we can then choose the superpotential according to the symmetries that we want to preserve. One motivation for the adjoint fields is as the additional degrees of freedom from an N = 2 supersymmetric gauge multiplet, and then the H u , H d fields become an N = 2 hypermultiplet; in this work we shall assume that N = 2 supersymmetry in the gauge/Higgs sector only is valid above some scale M N =2 . In this case, we can immediately write down the superpotential which contains the only interactions compatible with N = 2 SUSY and includes a central role for the R-symmetry. Indeed, under the R-symmetry of the N = 1 theory the adjoint scalars must have zero charge, and this prevents couplings of the form S 2 , S 3 etc which would otherwise be permitted by the gauge symmetry. The condition of N = 2 supersymmetry imposes (where g Y , g 2 are the hypercharge and SU (2) gauge couplings) at the scale M N =2 , which we shall in general take to be greater than the N = 1 SUSY scale. We must also add supersymmetry-breaking terms, and these do not necessarily need to respect the same symmetries as supersymmetric terms. The most general choice that we can make for the Higgs and adjoint scalar sector for the standard soft terms is where W iα are the supersymmetric gauge field strengths.
Since we are interested in Dirac gaugino masses and their attractive theoretical and phenomenological properties, we should expect that the terms that violate R-symmetry should be small: this includes the Majorana gaugino masses; A S , A T ; but also B μ . However, we require that the Rsymmetry is broken at some scale, since we believe that global symmetries cannot be exact; but also, in this model, the Higgs must carry R-charge and so the absence of an Raxion requires it. Indeed, the R-axion is essentially the Higgs pseudoscalar, whose mass is controlled by the B μ term. We therefore, as in earlier works, take B μ to have a small but non-zero value. We can also take motivation from models of gauge mediation of supersymmetry [20,68], where the trilinears are all small, and we shall mostly neglect them in the following (although they do not significantly affect the analysis).
On the other hand, in gauge-mediated models the adjoint scalars are typically the heaviest states. Taking large m S , m T , m O then motivates integrating them out of the light spectrum. Interestingly, since B μ should remain small due to the approximate R-symmetry, if we were to tune the Higgs masses such that only one remains light, then we would have very large tan β, and would have trouble obtaining the correct Yukawa couplings for the down-type quarks and leptons. This implies that a second Higgs should be taken to be somewhat light, and motivates studying the two-Higgs doublet limit of the model.
Finally, we note that this model does not have gaugecoupling unification. If we wish to naturally restore gauge coupling unification, we can add additional vector-like lepton fields, as was done in [55,64]. Since they are vector-like, we could also allow them to be hypermultiplets of the N = 2 at M N =2 , but their inclusion will little change the discussion in this paper so for sake of generality we shall neglect them.

The MRSSM
Another very popular realisation of Dirac gaugino models is the MRSSM [17,54,65,66,70]. In this model, we preserve an exact continuous R-symmetry by including some R-Higgs doublet superfields which couple to the Higgs bosons but do not obtain an expectation value, allowing the Higgs doublets H u , H d to have zero R-charge. The Higgs superpotential becomes 1 where the difference in sign is explained by the different charges of the hypermultiplets. 2 R-symmetry then limits the possible soft-supersymmetry breaking terms to consist of only the supersoft operator, squark/slepton masses and The terms on the last line are usually neglected, but there is no symmetry that forbids them (even if we expect them to be small in e.g. gauge mediation models).

Two-Higgs doublet model limit
The Higgs sectors of the models in the previous subsection have been comprehensively studied. However, here we wish to map them onto the two Higgs doublet model once To map our supersymmetric model onto this, we choose to make the identification from which we can write down The parameters λ i were given at tree-level and with some loop corrections in [21,44] in the limit of neglecting μ and m DY , m D2 . However, when we integrate out the adjoint scalars and retain these terms, there are corrections due to the presence of trilinear couplings; setting the parameters A S , A T to zero, we find for the minimal model: (2.11) Here we have defined In fact, the terms suppressed by m S R , m T P all have the effect of suppressing the Higgs quartic coupling: in the limit of large Dirac gaugino masses so that we can neglect m 2 S , B S , m 2 T , B T we find This simply corresponds to the well-known fact (see e.g. [14]) that the adjoint scalars eliminate the D-term potential of the Higgs, because they couple via the D-term. Writing φ i for (anti)fundamental scalars and for adjoint scalars, we have where T a are the generators of the gauge group with coupling g, and we see that the above will always be zero when we integrate out . For the MRSSM, for simplicity again neglecting A S , A T -for completeness we give the full corrections in Appendix D.4 -we find In this case, the supersoft limit is even worse, because in that limit all of the λ i vanish. However, even with the addi-tions of λ S and λ T in the minimal model, the potential is not stable in this limit -for example if H d or H u are set to zero the quartic terms vanish -and so we would require loop corrections to prevent runaway vacua. An investigation of whether this is even viable is beyond the scope of this paper: instead, since we do not want to substantially reduce the Higgs quartic coupling at low scales we shall consider instead that |m DY | m S , |m D2 | m T . As is also well known (see e.g. [32,44]) and we shall later discuss, this limit is also imposed on us by electroweak precision tests. In this limit we have instead at tree-level Hence for the rest of the paper we shall consider our lowenergy theory to be a type-II two Higgs doublet model with an additional (Dirac) bino and wino (the gluino must remain heavy due to LHC constraints -currently of the order of 2 TeV). We shall fix the boundary conditions at high energies and find some interesting conclusions.

Tree-level alignment
In [16,73] the Higgs sector of Dirac gaugino models was investigated in the limit that the couplings λ S , λ T took their N = 2 supersymmetric values at the low energy scale. However, they also pointed out that alignment in the Higgs sector would be broken by quantum corrections to the (2, 2) element of the Higgs mass matrix. In this section we shall consider just the potential at tree-level, and in Sect. 4 consider loop corrections, contrasting our results with theirs.
To begin with, the mass-matrices for the CP-even neutral scalars in the two-Higgs doublet model can be parametrised in the alignment basis where is (see e.g. [1][2][3]) where, using λ 345 ≡ λ 3 + λ 4 + λ 5 we have The parameter m A is the pseudoscalar mass, given by while the charged Higgs mass is The neutral Higgs masses are For our minimal model we have and (2.27) The Higgs mass matrix is diagonalised to find the physical Higgs masses and the mixing angle α. From the identification of the 2HDM parameters in (2.10) we obtain (2.28) The condition for alignment is the diagonalisation of M 2 i.e. Z 6 → 0. From Eq. (2.26) we see this amounts at tree-level to having In other words, when the couplings respect their N = 2 values, the Higgs doublets are automatically aligned! From Eqs. (2.28, 2.29) we find that in this alignment limit, c β−α → 0 and s β−α → 1, therefore the heavy CP-even neutral scalar doest not take part in electroweak symmetry breaking while h is a Standard Model Higgs-like boson. The tree-level masses of the two neutral CP-even Higgs bosons are while the charged Higgs boson mass is given by correcting the expression given in [16,73]. Hence, at treelevel, the model exhibits alignment for any value of tan β and the tree-level Higgs mass is independent of tan β (which was already noted in [16,73]). On the other hand, for the MRSSM there is no automatic alignment, because the Higgs sector at tree-level closely resembles that of the MSSM once the adjoint scalars and R-Higgs fields are decoupled; this can be seen just by putting λ S = λ T = 0 in the above equations. In the following we shall therefore mostly focus on the minimal Dirac gaugino model (with some further comments about the MRSSM).

Radiative corrections to alignment
As mentioned above, the perfect alignment obtained at treelevel is not preserved when the radiative corrections to the scalar effective potential are taken into account. In addition to the corrections already present in the MSSM, there are two new sources for this misalignment. The first is due to the appearance of chiral fields, quarks and leptons, at a scale M N =2 . This scale can be identified with the fundamental scale of the theory, or an intermediate scale where a partial breaking N = 2 → N = 1 is achieved (while an explicit realisation of this partial supersymmetry breaking remains unknown for a chiral theory, there is not a no-go theorem showing it to be impossible). The second large contribution comes from the mass splitting between fermonic and bosonic components of all of the superfields, i.e. coming from the N = 2 → N = 0 (or N = 1 → N = 0) breaking. We will discuss them here in turn.

Misalignment from
When we run our couplings from the N = 2 scale M N =2 to the scale of the N = 1 supersymmetric superparticles (which we shall call M SUSY ) there will be a splitting induced of λ S and λ T relative to the N = 2 SUSY relations. This in turn will lead to misalignment at M SUSY via a non-zero Z 6 : To obtain an estimate of the magnitude of this splitting, we can integrate over the difference in the beta functions for λ S and λ T to leading order: These equations are only useful for small M N =2 /M SUSY , because for large ratios the top Yukawa coupling can change by a factor of two or more, but it gives an indication of the amount of misalignment: even for log M N =2 M SUSY ∼ O(10) we find This is a small deviation from alignment indeed, and very encouraging. We shall investigate this quantitatively in Sect. 4, and will find that due to the diminishing Yukawa couplings at high energies, 3 the actual splitting is smaller than this naive estimate. As an aside, a similar conclusion is reached if we extend our Dirac Gaugino theory by including additional fields to the Minimal Dirac Gaugino Supersymmetric standard Model (MDGSSM) to restore gauge coupling unification.

Misalignment from
More significantly, there is the potential misalignment induced from the threshold corrections at M SUSY and then the running between M SUSY and the scale of the THDM; let us take the matching scale (as commonly done) to be the electroweak vev v. These can both be approximated at one loop by corrections to the δλ i . In the approximation that the singlet and triplet scalars are degenerate with mass M , and the stop squarks are degenerate with mass mt and neglecting the splitting between the couplings λ S,T and their N = 2 values we find, matching at a scale μ: using y t , y b , y τ to denote the top, bottom and τ Yukawa couplings. We give full (updated) expressions in the limit m DY , m D2 m S , m T in Appendix 1. We then find the remarkable result that the singlet/triplet scalar contributions to Z 6 exactly cancel out! We then find that the dominant contribution to Z 6 is that coming from the stops: where m t is the top quark mass. Although the magnitude of this is the same as the loop contribution to Z 6 in the MSSM, the misalignment thus induced is much smaller, because (a) there is no tree-level contribution, and (b) it is also proportional to the stop correction to the Higgs mass, which is smaller than in the MSSM due to the tree-level boost to the Higgs mass. To investigate the misalignment in this model further, however, we shall in the next section perform a precision study using numerical tools, where we shall use the logic of the hMSSM [74]/h2MSSM [73] to show that the misalignment in the model is even smaller than the above naive estimate.

Precision study
To precisely study the quantum corrections to alignment in our minimal model, we implemented the low-energy model consisting of the THDM supplemented by a Dirac bino and a Dirac Wino into the package SARAH. We describe the couplings of the model in detail in Appendix 1. We then modified the code to implement the boundary at a supersymmetry scale M SUSY and use the two-loop supersymmetric RGEs for the minimal Dirac gaugino extension of the MSSM [46,75] as generated by SARAH [76][77][78]. While this theory does not fit into a GUT and has no gauge coupling unification, we implemented an N = 2 supersymmetry scale where the couplings λ S , λ T take their N = 2 supersymmetric values. By running all the way from a low-scale Q (which we take to be the scale of the Dirac gauginos and Heavy Higgses, but could equally be m top ) up to M SUSY and then M N =2 and back down, iterating until the results converge, we were able to find consistent values of the parameters. At the scale Q, the threshold corrections are those that are included in SARAH by default: • One-loop matching of Yukawa couplings to the Standard Model values, plus two-loop strong corrections to the top Yukawa. • One-loop gauge threshold corrections.
We employed the two-loop RGEs for this model up to M SUSY , and then at the scale M SUSY , we implemented the following thresholds: • Tree-level correction to the λ i from Dirac gaugino masses given in (2.11), even if we are otherwise neglecting the Dirac gaugino masses. We take M SUSY to be a common mass of left-and righthanded stops, and assume that other MSSM particles have masses at this scale; we allow the singlet at triplet scalars to be heavier at a scale M . We eliminate all R-symmetryviolating terms (such as squark trilinear couplings) and assume that This means that we neglect squark mixing, which greatly simplifies the thresholds. The thresholds for supersymmetric particles that we include are then nearly complete in this limit: the gauge and Yukawa threshold corrections vanish for the MSSM couplings, and we only neglect the corrections to the gauge/Yukawas induced by the adjoint scalars -since their effect is in general very small; we leave the calculation and implementation of these for future work. However, we do include their contribution to the Higgs quartic couplings. Furthermore, we know that in the limit of zero squark mixing the two-loop corrections to the Higgs quartic couplings are also small or even vanishing [85], and so we are justified in neglecting them.
To perform a more general scan over the parameter space including trilinear scalar couplings, general masses and allowing μ, m DY , m D2 to be of the order of M SUSY we would need to compute the additional threshold corrections. While we expect that the effect of μ, m DY , m D2 on our results will be very small, it would nonetheless be interesting to compute these in the future.
We performed scans over the values of tan β and varied M SUSY to obtain a light Higgs mass of 125.15 GeV. For the other values we take by imposing As we shall later see, these are compatible with all current experimental constraints. Note, on the other hand, that we shall not discuss collider limits on the electroweakinos because the effect of changing their masses is tiny.
In the scans we see little deviation between m tree A and the mass of the heavy/charged Higgses because the mixing is small; indeed the results are not especially sensitive to m tree A as a result.

Running from the N = 2 scale
At the scale Q = 400 GeV, we find g Y = 0.37, g 2 = 0.64 ± 0.01. These are barely different at the SUSY scale and vary little with M N =2 , but we do find some dependence of the ratios √ 2λ S /g Y , √ 2λ T /g 2 on this scale, which we give in Fig. 1. The values in the plot were taken with a common supersymmetric scale of M SUSY = 3 TeV and have essentially no dependence on m A .
An alternative way of visualising this information is in the quantity Z 1 evaluated at the SUSY scale. Since our model is always very near alignment, this gives the "tree-level" Higgs mass and so in Fig. 2 . We see that for M N =2 = M SUSY this is always essentially M Z , while as we increase M N =2 we obtain a further enhancement to the Higgs mass at small tan β 1.5.
If we were to include no further corrections, then the value of Z 6 at M SUSY would be given by 10 10 GeV and 10 16 GeV 10 10 GeV and 10 16 GeV, which corresponds to the "tree-level" value of the Higgs mass before we take running from M SUSY (or equivalently the SUSY corrections at M Z ) into account (we take v = 246 GeV in the figure). We see that increasing M N =2 increases the Higgs mass, particularly for small tan β > 1. 5 for the value of Z 6 at the SUSY scale almost exactly correspond to the above equation, which we plot in Fig. 3. The differences (particularly the tiny difference from zero for the N = 2 scale equal to M SUSY ) come from the tree-level and loop-level shifts.

Running below M SUSY
Once we include the two-loop running below M SUSY , the picture changes substantially. This is dominated by the effects of the stops via their absence from the RGEs; we plot the results of Z 6 for the same scan as in Fig. 3 at the scale of our low-energy theory in Fig. 4.
Interestingly, the results can be understood by following the reasoning of the hMSSM [74]/h2MSSM [73] treatment. In that framework, the quantum corrections to the Higgs mass are assumed to be dominated by the (2, 2) component -and further that we can neglect the contributions to the other components compared to the tree-level ones. We shall first review what happens in the hMSSM and then apply the analysis to our case.

(Lack of) alignment in the hMSSM
In the hMSSM [74], we have λ 2 = Now let us suppose that we tune the values to obtain alignment. We then have which leads to i.e. it is impossible to achieve alignment without decoupling or going to the large tan β limit with these approximations. If we do not neglect the other contributions to Z 6 , in the case of exact alignment we then have Since we expect λ 1 , λ 345 λ 2 , and for t β > 1 we have s 2 β > c 2 β , s 2 β > |c 2β |, this is still impossible to satisfy. 4 However, we will find that for our scenario things are somewhat better.

Alignment in the Dirac-gaugino model
Using the expressions (2.16) for the quartic couplings, we can rewrite This leads to The corrections δλ i can be interpreted as either coming from running the couplings between the scale M SUSY and Q, or alternatively from integrating out the supersymmetric particles at the scale Q. In the latter case we can obtain an estimate of their values from the expressions (3.5) and see that they are typically suppressed relative to /v 2 by a numerical factor and also the ratio of the electroweak gauge coupling to the strong gauge coupling or top Yukawa, and we find that we can therefore continue with the hMSSM approximation and neglect them. However, the effect from the running of λ S , λ T is non-negligible: eliminating in exchange for the Higgs mass and defininĝ (4.10) 4 In the full MSSM the radiative corrections to λ 345 and λ 7 -which we are neglecting -can be large enough to allow alignment for tan β 10 [6,86]. However, this also requires significant stop mixing, which we do not have in our model.
we have We shall later give the expressions for eliminating λ 2 and calculating Z 6 in any THDM with general λ i , i = 1 . . . 4 in Eqs. (6.3) and (6.4).
A comparison of the above formula with the curves in Fig. 4 shows that this gives a reasonable fit. In the case of M N =2 = M SUSY the expression is particularly simple, but in the other cases we need to take account of the varation of √ 2λ S (M SUSY ), √ 2λ T (M SUSY ) with tan β that can be seen in Fig. 1.
The main conclusion that can be drawn from the above formula is that the misalignment coming from the squark corrections required to enhance the Higgs mass can be compensated by the effect of running λ S , λ T . Indeed, we see from Fig. 4 that for M N =2 = 10 16 GeV, Z 6 is essentially vanishing for tan β 3. From the curves in the figure, we see that increasing the N = 2 scale causes a partial or total cancellation of the misalignment contributions, meaning that the Higgs boson is accidentally very Standard-Model-like, independent of the mass of the heavy Higgs! This is the main result of the paper.

Higgs mass bounds on the SUSY scale
Finally we consider the effect of the loop corrections in the low-energy theory on the Higgs mass (i.e. those coming from the Higgs sector itself, the top and the electroweakinos). In Fig. 5 we show the tree-level and one-loop values for the Higgs mass as we vary tan β (with M SUSY fixed to ensure m h = 125.15 GeV at two loops). We find a significant upward shift of about 7 GeV at one-loop, and then a downward shift of about 1 or 2 GeV from one to two loops. Note that we can interpret the "tree-level" Higgs mass as the loop-level Higgs mass in the full Dirac gaugino model including the effects of the stops and gluinos (which in the EFT formalism appear via the RGEs, rather than fixed-order diagrams).
In Fig. 6 we show the final curve of tan β against M SUSY , for different values of the N = 2 scale between M SUSY and 10 16 GeV.
The plot shows that there is a minimum for M SUSY around tan β 2 or 3, particularly for larger values of M N =2 , which  can be understood in terms of the splitting of λ T from its N = 2 value and the consequent boost to the Higgs mass, which can be clearly seen in Fig. 2.
The results in Fig. 6 contrast starkly with the MSSM case matched onto the 2HDM as shown in e.g. [87]: due to the enhancement to the Higgs mass from the new couplings already seen in Fig. 2 we have a much lower SUSY scale. On the other hand, there are significant differences from the values quoted in [73] which are most closely related to the case M N =2 = M SUSY ; here of course we have light electroweakinos, although the largest difference is the significantly more accurate EFT calculation employed here.

Experimental constraints
Since our model realises excellent alignment, the light Higgs couplings are very nearly Standard-Model-like across the whole parameter space, and so there is no significant constraint from those -this is in contrast to e.g. the hMSSM scenario, where for low tan β the Higgs couplings provided until recently the most important lower bound on the Heavy Higgs mass. However, there are still significant constraints on the parameter space coming from electroweak precision tests, flavour and direct searches, as we detail below.

Electroweak precision corrections
There are two contributions to the electroweak precision parameters: those coming from the high-energy theory, and those coming from the low-energy theory. In the high-energy theory there will be contributions at tree-level from the triplet scalars: they should obtain a vacuum expectation value, and in our EFT this manifests itself as generating effective operators.
In the limit of zero CP violation, and neglecting the terms A S , A T we can write the effective operator arising from integrating out the triplet as quite simply where we understand summation on the index a and When we give a vacuum expectation value to the Higgs, this translates into the constraint from the expectation value of the triplet: For μ = 500 GeV and an approximately N = 2 value for λ T , with small tan β insisting that this contribution does not exceed the experimental bound by 3σ gives m T P > 1500 GeV (5.5) while simply saturating without exceeding the central best-fit value would limit instead m T P > 2 TeV. On the other hand, we also have a contribution from the electroweakinos at loop level, which increases as the Dirac mass/μ-term become smaller. Hence they cannot be arbitrarily light. In Fig. 7 we plot the value of ρ calculated in the low-energy theory for the scan values (4.1) and find that they are below the experimental limit across the whole parameter space.

Bounds on tan β and m A
The most stringent constraints on the parameter space of our model come from the searches for pp → H/A → τ τ at the LHC; and the decay B → sγ determined in [89], which bounds the charged Higgs mass to be heavier than 580 GeV independent of the value of tan β (which in turn bounds the mass of the pseudoscalar Higgs to be above around 568 GeV).
The bounds from run 1 of the LHC were rather mild on the hMSSM: they restricted tan β < 8 for low m A (see e.g. [86,91]). In [73] it was claimed that in the h2MSSM these bounds would apply unaltered; while it is true that the couplings to the pseudoscalar are the same in the h2MSSM and hMSSM, the "heavy" Higgs does have altered couplings at small m A and tan β -since it is more aligned. Since the production of the Heavy Higgs is dominated at small tan β by gluon fusion, and at large tan β by the bbH process, then we would expect some differences at small tan β. However, recently, ATLAS produced a much enhanced bound [90] on gluon fusion and bbH production and then decay to τ pairs; they also interpreted this in terms of the hMSSM. To compare to our model we computed Higgs production using SusHi [92][93][94][95][96][97][98] and rescaled the production cross-sections according to b-quark and gluon couplings computed in our SARAH/SPheno code, then multiplied by the tau decay branching fraction, and combined the bound assuming that the signals from H/A production overlap for small mass differences. We show the result in Fig. 8, where we also show the bound from [90] on the hMSSM. We find almost Fig. 7 ρ calculated at one-loop in the low-energy theory, for different values of M N =2 given in the legend. We see that the magnitude is roughly equal to the experimental error, and we are always well within 3σ of the experimental central value (which is anyway above the Standard Model value by 1.6σ )

Fig. 8 Bounds from
pp → H/A → τ + τ − (blue region) and B → sγ (red region, m A 568 GeV) interpreted in the m A / tan β plane for the hMSSM (taken from [90]) and our model no difference, except that the bound on our model is very slightly weaker once decays to the electroweakinos are permitted. However, the branching ratio to electroweakinos in that region is never significant enough to reduce the τ decay fraction.

Alignment in the MRSSM
For completeness we now discuss the case of the MRSSM in the same limit as for the DG-MSSM. Since the tree-level THDM parameters are the same as those of the MSSM in the limit of large adjoint scalar and R-Higgs masses, 5 there is no contribution to Z 6 from the running of the parameters 5 So therefore Z 1 , Z 5 , Z 6 are the same as in the MSSM case, i.e. Eqs. (2.24)-(2.26) with λ S = λ T = 0. λ S u,d , λ T u,d . We can first write the neutral Higgs mass matrix as If we consider the loop corrections due to λ S u,d , λ T u,d to be small, then the analysis of alignment is identical to the MSSM case, and we can apply the hMSSM logic. However, if we instead take them to be non-negligible -such as in [53,54,65,66] -then the contributions to λ 2 no longer dominate, and the hMSSM reasoning may no longer apply. On the other hand, the largest contribution from the other particles will still be to λ 2 , and so we can assume that To eliminate , we eliminate λ 2 in terms of the Higgs mass, which for general λ i , i = 1 . . . 4 (and λ 5 = λ 6 = λ 7 = 0): this can then be substituted into the expression for Z 6 : We give the loop corrections to the λ i from the adjoint scalars in Appendix C.3, but in the simplified case of m T + = m T − = m S R = m S I = M and g Y = g 2 = 0 we have, for matching at a scale μ: If we then take (as in [54,65,66]) λ S u = −λ S d ≡ λ, λ T u = λ T d ≡ , and allow an additional contribution /v 2 to λ 2 from the stops, then we have We see that when the couplings λ, are large enough, the alignment will always be improved compared to the MSSM, because the enhancement to Z 1 is always greater than that to Z 6 . We note three cases of particular interest: 1. If we increase the contributions from the adjoint scalars to the point that we can neglect those from the stops, then we see that for small tan β we will easily have alignment (in contrast to the MSSM case).
2. Alternatively, we could enhance the contributions from λ rather than , since the former coupling does not contribute to Z 6 . 3. On the other hand, if we take the N = 2 supersymmetric limit we find, using the expressions in Appendix C.3 (and no longer neglecting the gauge couplings): giving δλ 345 = δλ 1 ≡ δλ, (6.9) so there is no shift to Z 6 from the adjoint scalars, but we do have a shift to Z 1 , i.e If the mass of the adjoint scalars is comparable to the mass of the stops, then this will however never be significant.
On the other hand, if we take the adjoint scalars to be very heavy, then this indicates that we can have improved alignment relative to the MSSM. To quantify this, we can use our above expression for Z 6 (6.4): which is the numerator for the MSSM case. In the case that m 2 A m 2 h (which corresponds to our case of interest -even though we would like m A small enough to not entirely be in the decoupling limit), we therefore find and so the deviation of Z 6 from the MSSM value due to the adjoint scalars should be less than 4%. On the other hand, as we shall see below, they can still have a significant effect on the SUSY scale.
Therefore, from the analysis above, in all three cases of interest, the alignment will never be as good as for our minimal Dirac gaugino model, because of the tree-level contribution to misalignment: we shall illustrate this for the N = 2 case in the next subsection.

Numerical analysis of an N = 2 MRSSM
To compare with our previous analysis of the DG-MSSM, here we present a simplified numerical analysis for an N = 2 MRSSM, as defined in point 3 above and Eq. (6.7). From our estimations above, the alignment should only differ from the MSSM when relatively extreme values are taken for the adjoint scalar masses, and so to perform a precise analysis we would need to have a tower of effective field theories and the appropriate threshold corrections. Instead we decided to neglect all loop-level threshold corrections other than those from the adjoint scalars (although we use 2-loop RGEs throughout) and performed a simple analysis where the low energy model was approximated by the Standard Model and type-II two-Higgs doublet model. In this way we should obtain an idea of how the adjoint scalar masses cause the SUSY scale and alignment to vary from the predictions of the MSSM.

Procedure
Two-loop Standard Model matching values were implemented at m t for the standard model gauge, Yukawa, and Higgs quartic couplings from [99] and a two-loop Standard Model running was performed up to an intermediate scale Q = 600 GeV, where the λ i (Q) couplings were given approximate values to be determined through future iterations between the scales Q and M N =2 . The two-Higgs doublet model 2-loop running was implemented up to the supersymmetry breaking scale defining the leading squark masses, M SU SY , where guesses were made for the inputs of the parameters λ S,T u,d . The MRSSM was then run to 2-loops to some high scale M N =2 where the N = 2 boundary conditions (6.7) were implemented. All two-loop beta functions were generated in SARAH, and the value of m tree A = 600 GeV was taken as in the minimal Dirac gaugino case. In this simplified model, as the electroweakinos are not taken to be light, the intermediate scale Q is taken to match the choice of heavy Higgs mass. Indeed, with these choices we should understand the Dirac gaugino masses m DY , m D2 , m D3 and the higgsino mass to be at M SUSY , and also the masses of the R-Higgs fields R u,d should be at that scale, because we do not implement any threshold corrections from those fields (leaving these to future work).
On the run down, λ i (M SU SY ) were matched to the 1-loop threshold corrections coming from the heavy S, T scalars as given in Appendix C.3, taking the adjoint scalars to be degenerate with mass M . This process was iterated, re-matching the gauge and Yukawa couplings onto their 2-loop Standard Model running values at the scale Q, while the λ i and λ S,T u,d couplings were matched to the outputs from the previous running until their values converged. Finally, the λ i parameters were mapped back onto the Higgs quartic coupling using λ(Q) = Z 1 (Q) and the Standard Model couplings were run back down to m t . The correct Higgs mass was selected from the criterion λ(m t ) = 0.252±0.002, corresponding to a pole Higgs mass of m h = 125 ± 0.5 GeV.
This process was executed for scans over the values tan β ∈ [2,20] Here the modulus of the ratio is plotted, since the λ S d ratio is negative to respect the N = 2 supersymmetry relations. As expected, the model is closest to the alignment limit when the N = 2 scale is closer to the M SU SY scale. It can be seen that the Higgs mass is boosted to a greater extent by the down-type couplings than the up-type, where the ratio √ 2λ T d /g has the largest effect, especially for higher values of N = 2 scale.    Figure 12 shows little deviation in the results for Z 6 (Q), regardless of M and M N =2 . Indeed, as anticipated above, the results are almost indistinguishable from the MSSM case, since the adjoint scalars in the MRSSM never give a large boost to the quartic couplings even for the extreme cases we  have taken. Exceptionally, the couplings in the case of very heavy scalars and very high M N =2 are considerably enhanced and deviate from the N = 2 relations, making the alignment in this case just marginally worse. While the adjoint scalars give only a very small boost to the Higgs mass, on the other hand it is enough to cause noticeable effects in the predicted M SU SY scale, shown in Fig. 13, because of the logarithmic nature of the contributions from other SUSY states. Figure 13 shows the values of M SU SY against tan β over the parameter scan producing a Higgs mass corresponding to m h = 125 ± 0.5 GeV: this margin is reflected in the enclosed transparent area. 6 For tan β < 4, M SU SY is required to be, at the very least, 20 TeV for the highest values of M N =2 and very heavy M , and is closer to ∼ 100 TeV for lower values. M SU SY stabilises around tan β = 10 for all values of M and M N =2 , where at this point M SU SY can be as low as several hundred GeV for M N =2 = 10 16 and very heavy scalars. In this final extreme case (which is of course excluded experimentally, but we give as an indication of the possible effects) the logarithms being resummed in the RGEs become smaller, it is possible that any neglected threshold effects could make a significant difference and the results become unreliable, but we leave this additional analysis to future work.

Conclusions
We have considered the the consequences for the simplest realisation of Dirac gaugino models when we impose N = 2 supersymmetric boundary conditions for the Higgs/gauge sector at some energy scale. We found that the model naturally realises alignment in the Higgs mass matrix, and that surprisingly this is preserved even by quantum corrections. Even more interestingly, the departure from N = 2 relations due to running of the couplings actually leads to both an enhanced Higgs mass (and thus lower SUSY scale/more natural model) and also improved alignment when we take the effects of the squarks into account.
We have provided the most accurate calculation to date for the SUSY scale for a Dirac gaugino model by employing the effective field theory approach, with one-loop boundary conditions at the high scale and two loops at the THDM scale. This leads to the prediction that the scale of coloured superpartners should be above 3 TeV (when we allow a very high scale for the breaking of the approximate N = 2 SUSY) but across most of the parameter space it is below 10 TeV. While this is not encouraging for the detection of stops/gluinos at the LHC, this is well within the reach of a future 100 TeV collider. On the other hand, the LHC or a future e + e − collider should be able to explore the electroweak sector of the model, including the Higgs sector and the electroweakinos (if they are light).
There are many possible avenues for future work: improving the accuracy of the matching at M SUSY (as noted recently, matching at two-loop order is often necessary for accuracy of the loop expansion to include all non-logarithmic corrections [100], although in this class of models as we have discussed all of the missing corrections are believed to be small) and including the effects of the electroweakinos in the matching at one loop, so that we can consider the model with m DY ∼ m D2 ∼ M SUSY ; also with the full set of thresholds we could perform an estimate of the error in the calculation (which, again, should already be small -see e.g. the estimates for the MSSM case in [101]); or including the effects of possible R-symmetry violating terms. On the other hand, it would also be interesting to more fully explore the consequences for different Dirac gaugino models, such as the MRSSM, where we have only performed a preliminary analysis.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

A THDM with light electroweakinos
The limit that we are interested in has the electroweakinos much lighter than the singlet and triplet scalars; in order to avoid washing out the tree-level Higgs quartic coupling and generating a large contribution to ρ they should be light. At energies below the stop/sbottom masses, then, we have an effective theory of the two-Higgs doublet model augmented by light electroweakinos. This looks a little like Split supersymmetry or the scenario of [102] (which considered a split scenario with both Higgs doublets light), except that our electroweakinos have Dirac masses and our gluino is heavy; here are therefore new Yukawa couplings between the Higgs doublets i , the left and right binoB i and winoW a i for i = 1, 2, and the higgsinosh u,d : This gives neutral and charged fermion mass matrices where the bases are Note that M B , M W are 2 × 2 matrices. At the SUSY scale, we match the above to the corresponding couplings in the Dirac gaugino theory: We choose to make the definition meaning H u · H d ↔ − † 1 2 , which leads to the identifications These are, however, given in terms of the DR parameters: making the conversion to MS we find

B Threshold corrections
In this section we give the one-loop threshold corrections to the couplings in our theory. Throughout we use the definitions where μ is the renormalisation scale at which the quantities are evaluated.

B.1 Conversion from MS to DR
The conversion of the gauge couplings from the MS to DR renormalisation scheme is given by For the Yukawa couplings, we retain only the strong gauge coupling dependence: For the Higgs quartic couplings, we define and then If we express the quartic couplings in terms of the MS gauge couplings at tree level, then we have a further shift from the shift to g 2 of + κ 6 g 4 2 for λ 1,2 and − κ 6 g 4 2 for λ 3 .

B.2.1 Matching at the SUSY scale
In the limit that we take in the body of the paper, all of the threshold corrections coming from squarks vanish at the matching scale. However, to extend the results of [87,103] to our model, we have computed the corrections coming from stops, sbottoms and staus to the quartic couplings allowing non-zero squark trilinears and μ. They are given by where the δ (1) th λ i contributions are those from bubble, triangle and box diagrams and are unchanged from the MSSM case given in [87,103], while the δ (1) λ i are the wavefunction corrections that are modified for our model:

B.2.2 Matching at a general scale
If the squarks are not degenerate or we integrate them out at a scale other than a common SUSY scale, then in our limit we have If we consider all squarks to be at a common SUSY scale M S , then these simplify to (B.9)

B.3 Contributions from the S, T scalars
Here we present the contributions to the quartic couplings coming from the adjoint scalars S, T , in the limit m DY , m D2 m S , m T , B S , B T and assuming no CPviolation. The scalars have masses The loop corrections to the quartic couplings are: These results update those previously given in the literature by including the electroweak contributions.

C One-loop RGEs
For our numerical study we use two-loop RGEs throughout, as generated by SARAH. They are too long to put into print; however, for illustration we provide here the one-loop expressions for the low-energy theory of the THDM with electroweakinos, after making the simplification that: • Only third generation Yukawa couplings are included.
• No CP-violation, hence all couplings real.
• Once we respect the matching conditions (A.5-A.8), the beta functions for the couplings that are zero at the supersmmetry scale are zero along the flow, and hence we set the couplings g 11 1d , g 21 1u , g 22 1d , g 12 1u , g 11 2d , g 21 2u , g 22 2d , g 12 2u to zero.

D MRSSM corrections
Here we collect the tree-level and leading one-loop threshold corrections to the THDM paramters in the MRSSM.

D.4 Tree-level
The tree-level λ i are given by (D.1) The shifts from integrating out the adjoint scalars give The one-loop corrections from the adjoint scalars in the limit that we can neglect the Dirac gaugino masses are given by: