Large muon (g − 2) with TeV-scale SUSY masses for tan β → ∞

The muon anomalous magnetic moment aμ is investigated in the MSSM for tan β → ∞. This is an attractive example of radiative muon mass generation with completely different qualitative parameter dependence compared to the MSSM with the usual, finite tan β. The observed, positive difference between the experimental and Standard Model values can only be explained if there are mass splittings, such that bino contributions dominate over wino ones. The two most promising cases are characterized either by large Higgsino mass μ or by large left-handed smuon mass mL. The required mass splittings and the resulting aμSUSY are studied in detail. It is shown that the current discrepancy in aμ can be explained even in cases where all SUSY masses are at the TeV scale. The paper also presents useful analytical formulas, approximations for limiting cases, and benchmark points.

(1.1) It can be explained by a variety of new physics models below or at the TeV scale. Generally, the new physics contributions a NP µ are strongly correlated with the loop contributions to the muon mass, δm NP µ , and they are suppressed by two powers of the typical new physics mass scale M NP [7,22]. This relation can be written as where C NP = O(δm NP µ /m µ ) is given by the model-dependent relative contribution to the muon mass.
Of special interest are, therefore, models of radiative muon mass generation: in these models, δm NP µ amounts to the entire physical muon mass and C NP = O(1). Excluding fine-tuning in the muon mass, these models yield the largest a NP µ , compared to other ones 1 The experimental result has been obtained at BNL [1]; a fourfold improvement in precision is expected from the new experiments [2][3][4][5]. For reviews of the theory prediction see refs. [6,7]; recent theory progress has been achieved on the QED [8][9][10][11], electroweak [12], and hadronic contributions [13][14][15][16][17]; for specific reviews of the hadronic light-by-light contributions and expected future improvements see [18][19][20][21]. The value quoted in the text is based on the hadronic contributions from ref. [13].

JHEP10(2015)026
with the same new physics scale M NP . As the estimate a NP µ = O(m 2 µ /M 2 NP ) holds, the observed deviation can in principle be explained by a new physics scale of the order 2 TeV.
The general idea of radiative muon mass generation can be realized in the renormalizable and calculable framework of the minimal supersymmetric standard model (MSSM). As discussed below in section 2 there are two distinct possibilities, one of which has already been studied in refs. [23][24][25]. The possibility considered in the present paper is the limit with the up-and down-type Higgs vacuum expectation values v u,d . In this limit v d , the muon mass, and all other down-type lepton and quark masses vanish at tree level. The masses arise from finite loop diagrams generating non-holomorphic couplings of down-type fermions to the "wrong" Higgs doublet H u . These loop diagrams are also important for finite tan β, and they have been discussed extensively in the literature, often in the context of B-physics [26][27][28][29][30][31][32][33][34][35], but also in the context of the muon magnetic moment [36]. Because of these loop diagrams, the limit tan β → ∞ exists and is phenomenologically viable [37,38].
We consider a µ in the MSSM for tan β → ∞ for two reasons. On the one hand, this scenario is a calculable realization of the generally interesting idea of radiative muon mass generation, and the results can be indicative of the more general situation. In particular, it is no special case of the model-independent analysis of ref. [39], where simplified models with only two relevant particle masses have been considered. We will see that all cases of interest here involve at least three relevant masses at the TeV scale.
On the other hand, the tan β → ∞ limit opens up an intriguing area of supersymmetry (SUSY) parameter space where a µ behaves qualitatively very differently from the standard MSSM case (with moderate tan β). This standard MSSM case is well known (for reviews see [40][41][42]; recent works are [43][44][45][46][47][48][49][50][51][52][53][54][55][56][57]), and it requires SUSY masses in the few-hundred GeV range in order to explain the deviation (1.1); the LHC experiments, however, start to exclude parts of the relevant parameter space [43,58]. Thus it is well motivated to ask whether SUSY can explain the deviation even if all SUSY masses are much higher, at the TeV scale. The answer can be expected to be provided by SUSY radiative muon mass scenarios such as the one considered here.
Our approach is a low-energy phenomenological one, with the aim to answer whether and how TeV-scale SUSY masses are compatible with the deviation (1.1) in the MSSM limit tan β → ∞, also taking into account other experimental data and the internal consistency of the theory. Nevertheless it is important to briefly take a top-down perspective and review how infinite tan β might result from a more fundamental theory. One appealing possibility is an unbroken continuous R-symmetry, which forbids the Bµ soft breaking term. Such a symmetry has been used in ref. [59] to construct a "one Higgs doublet model", in which v d = 0. R-symmetric models such as this one however contain non-MSSM degrees of freedom, and therefore our MSSM study does not directly apply to them. A mechanism which can naturally generate vanishing/small Bµ and infinite/large tan β in the MSSM is gauge-mediated SUSY breaking. In its pure form [60] it requires Bµ = 0 at the messenger scale, and a non-vanishing Bµ is only generated by renormalization-group running at lower JHEP10(2015)026 scales. The connection between gauge-mediation and the large-tan β MSSM has been studied in detail in ref. [38], however with several simplifying restrictions which allow for very large but not infinite tan β. The promising results explained below provide motivation for further model building to actually realize infinite tan β in a more fundamental theory, but this is beyond the scope of the present paper.
We now give a brief preview of the MSSM limit tan β → ∞ to set the stage for the remainder of the paper. In the simple case that all relevant SUSY masses are equal to M SUSY and tan β is moderate, the SUSY contribution to a µ is approximately given by the one-loop diagrams as For large tan β it has already been shown in ref. [36] that higher-order terms leading in tan β can become important and change the linear dependence. These higher-order terms can be resummed as where, still in the case of equal SUSY masses, the self-energy leads to Hence, in the desired limit of infinite tan β the SUSY contribution becomes demonstrating that the magnitude of the contribution is very large, even for M SUSY in the multi-TeV range. The proportionality to tan β has been replaced by a constant behaviour, and the dependence on sign(µ) has disappeared. However, the sign is predicted to be wrong, so such a scenario with infinite tan β and equal SUSY masses is definitely excluded.
If the SUSY masses are not all equal, the above approximations do not apply, and in the remainder of the paper we will give a full investigation of the five-dimensional parameter space. In particular we will characterize the regions that lead to a positive SUSY contribution to a µ , which agrees with the observed deviation.
In section 2 we define the model, provide analytical results and useful approximation formulas, and briefly review constraints from other sectors. Comprehensive numerical analyses and their physical interpretation for different mass parameter regions are presented in section 3. We also discuss relevant constraints on the parameter space. Further discussion and the conclusions are given in section 4. In addition, we present relevant benchmark points and scan plots showing the possible values of the lightest SUSY mass parameter for which the discrepancy (1.1) can be explained.

Analytical results and approximations
The MSSM tree-level muon mass is given by the product of the muon Yukawa coupling y µ and the down-type Higgs vacuum expectation value v d . Radiative muon mass generation in the MSSM requires the tree-level mass to vanish. There are two generic possibilities: either y µ = 0 or v d = 0. In case y µ = 0, the muon mass can be generated from loop diagrams involving binos and smuons and the non-standard soft supersymmetry breaking term A µ . The value of a µ in this scenario has already been studied in refs. [23,24], and the observed deviation in a µ can indeed be explained for bino and smuon masses around 1 TeV. A similar study has been performed with the holomorphic trilinear coupling in ref. [25].
Here we consider the second possibility v d = 0, or equivalently the limit of eq. (1.3). In the following we explain the relevant formulas governing the muon mass, Yukawa coupling, and magnetic moment in this limit.
The physical muon (pole) mass m µ is given by the on-shell muon self energy 2 The explicit result can be found in the appendix. This equation is used to determine the value of the muon Yukawa coupling y µ . Eq. (2.2) also defines the "reduced" self energy factor ∆ red µ (which satisfies ∆ µ = tan β ∆ red µ with the usual definition for ∆ µ [36] and which agrees with the quantity µ of ref. [38]). The definition highlights the physics behind the muon mass generation: chiral symmetry (under which left-and right-handed muon transform with different phases) is broken by the Yukawa coupling, and the muon mass is generated by loop-induced couplings of the muon to the vacuum expectation value of the "wrong" Higgs doublet v u . Hence the factor y µ v u can be pulled out of all contributions to the muon mass. (Another factor of the Higgsino mass µ could also be pulled out of all contributions, because µ has to appear in all couplings of the muon to v u due to Peccei-Quinn symmetry.) It should be noted that the coupling constants, mixing matrices, and mass eigenvalues entering the self energy partially depend on the muon Yukawa coupling, so ∆ red µ has a small residual dependence on y µ . The combination is thus an important quantity. In the standard case with moderate tan β it is strictly equal to m tree µ tan β and also often identified with m µ tan β, but the latter identification is only possible if either higher-order effects can be ignored or on-shell renormalization is used.
Here we have to distinguish m µ from m tree µ and tan β eff from tan β. 2 All the following formulas are given for the limit tan β → ∞. The formulas valid exactly also for arbitrary (small, large or infinite) tan β can easily be reconstructed by replacing ∆ red

JHEP10(2015)026
The physics behind the generation of a µ is the same as for the muon mass. To highlight this similarity we write introducing another dimensionless "reduced" quantity a red µ . In the standard case, we simply have a SUSY µ ≈ tan β a red µ . Again, a red µ has only a small residual dependence on y µ , and its explicit result can be found in the appendix.
Combining the equations for the muon mass generation and the magnetic moment, we obtain the MSSM prediction in the limit tan β → ∞: We record here further useful relations between the quantities introduced so far: In contrast to eqs. (1.4), (1.6) from the Introduction, the exact result eq. (2.5) depends on five independent SUSY mass parameters: the Higgsino mass µ, the gaugino (bino and wino) masses M 1,2 , and the left-and right-handed smuon soft mass parameters m L,R . The result has interesting symmetry properties due to cancellations between numerator and denominator: on the one hand, a SUSY Now, useful mass-insertion approximations for the MSSM muon self energy and the SUSY contribution to the anomalous magnetic moment of the muon are provided. Their validity for the parameter ranges of interest will be discussed at the end of the subsequent section. To maximize the magnitude of a µ and minimize constraints from CP -violating observables we restrict our considerations to real mass parameters µ, M 1,2 . Taking advantage of the above-mentioned symmetries we choose, without loss of generality, µ and M 1 to be positive. Only the sign of M 2 remains arbitrary. Under the assumption M Z M SUSY , the SUSY contributions are given by the five mass-insertion diagrams in figure 1 [41,42,61]. The self energy factor ∆ red µ can then be decomposed as the sum of where g 1,2 are the U(1) and SU(2) gauge couplings, and the loop function is given by This loop function is dimensionful and scales as 1/M 2 , where M denotes the largest of the three mass arguments. The result for the self energy corresponds to the one from ref. [36] if the limit tan β → ∞ is taken. Likewise, for the case M Z M SUSY , the SUSY contribution to the anomalous magnetic moment of the muon can be approximated by the sum of the expressions corresponding to the mass-insertion diagrams in figure 1, with an external photon coupling to the charged internal line. Via eq. (2.4), this result can be related to the one for finite tan β quoted e.g. in refs. [42,47], where the loop functions have been introduced. The functions F C 2 and F N 2 can be found in the appendix.

JHEP10(2015)026
We now provide seminumerical approximations which allow to directly read off numerical orders of magnitude and signs and which facilitate the phenomenological discussion. The dimensionless loop functionŝ are especially useful if the sign and size of the different contributions shall be compared. The first two of these functions have been introduced in refs. [37] and [43] respectively. In the special case of three equal masses, we getÎ(1, 1) = 1/2,K W (1, 1) = 5/6 and K N (1, 1) = 1/6. In general,Î andK N lie between 0 and 1;K W lies between 0 and 1.155. Furthermore, in the case of equal arguments,Î(x, x) andK N (x, x) increase monotonically, whileK W (x, x) increases up to a value of about 1.155 at x ≈ 11.15 and then decreases towards a limit of 1 for x → ∞. By inserting all numerically known Standard Model quantities we obtain a seminumerical version of eqs. (2.7) and the negative MSSM muon self energy, which should be equal to the muon mass [see eq. (2.2)], can be written as where the order of terms is the same as in eqs. (2.7), and the first term combines the contributions from loops with charged or neutral winos. Similarly, by plugging in numbers into eqs. (2.9), we obtain the seminumerical massinsertion approximation for a SUSY µ in the tan β → ∞ limit as

JHEP10(2015)026
which corresponds to the result from ref. [43] in the limit tan β → ∞ and combines the contributions that are created by a charged or neutral internal wino in the first term. The approximate expressions (2.12) and (2.13) have several noteworthy features.
• They are linear in y µ because the residual dependence of ∆ red µ and a red µ on y µ vanishes for M Z /M SUSY → 0.
• All terms involve the factor µM 1 or µM 2 . This explains the symmetry of a SUSY µ [given by eq. (2.5)] under a sign change either of µ or of M 1 and M 2 .
• Under sign change of M 2 alone, the wino contributions change their signs relative to the bino contributions.
• The signs of the contributions are related: each mass-insertion diagram contributes with equal sign to ∆ red µ and a red µ , except for the chargino diagram, where the signs are opposite.
• Thus if the chargino diagram (a) in figure 1 dominates in both ∆ red µ and a red µ , we get a SUSY µ < 0. In contrast, for neutralino dominance a SUSY µ becomes positive and, therefore, has the correct sign with regard to an explanation of the discrepancy (1.1).
• The contributions (2.7b) and (2.9b) from the diagram with an internal neutral wino can never dominate since they are always smaller than the ones from the diagram with an internal charged wino. In ∆ red µ , the charged and neutral wino contribution add constructively, in a red µ destructively.

Experimental constraints from other sectors
Before we present our numerical analyses, we briefly discuss constraints on parameter space from other observables and sectors and show that the arising constraints can be satisfied without restricting the five parameters µ, M 1,2 , m L,R relevant for a µ . Clearly, an observable strongly related to a µ is the lepton flavour violating decay µ → eγ. The correlation of the two MSSM predictions has been discussed extensively in the literature [62][63][64][65][66]. If the prediction for a SUSY µ is fixed, the SUSY prediction for µ → eγ can be estimated quite well up to the unknown lepton flavour violating parameters. In particular the one-loop amplitudes for the two quantities share the same tan β enhancement; the correlation between the two does not strongly depend on tan β. E.g. figures 13-16 of ref. [65] remain valid also for tan β → ∞. Therefore, like for moderate tan β, the decay µ → eγ is compatible with experimental bounds [67] for sufficiently small lepton flavour violating parameters.
Further constraints from B-physics have been discussed extensively in ref. [38]; taking the limits of those results as tan β → ∞ shows that agreement between theory and experiment can be achieved. The parameters mainly constrained by B-physics are the heavy Higgs boson masses as well as the stop trilinear coupling A t . Since these parameters are not relevant for the discussion of a µ , we give only a brief account of how they are constrained. The new physics contribution to BR(B + → τ + ν) can be suppressed below 2σ with σ being the experimental error, by having the charged Higgs mass M H + around a few

JHEP10(2015)026
TeV or higher. We can relax the constraints from B s → µ + µ − and B → X s γ by assuming that A t and flavour-violating squark mass insertions are small enough. In addition, the former process can be further suppressed by raising the masses of the CP -odd as well as the CP -even heavier Higgses, and the latter by increasing the squark masses.
Given large bottom and tau Yukawa couplings, large sbottom and stau corrections to the mass of the lighter CP -even Higgs h might undermine the stop contribution that should lift m h up to the measured value around 125 GeV [68][69][70][71][72], as tan β → ∞. Indeed, figure 2 of ref. [73] shows a rapid drop of m h as tan β grows over around 40. This happens due to a cancellation between the tree-level and the ∆ b contributions to the bottom quark pole mass causing a blow-up of the bottom Yukawa coupling. As mentioned in the same reference, this is however not necessarily the case for a high tan β. First of all, it never happens for µ < 0. Even if µ > 0, once we enter the regime where |∆ b | 1, the bottom Yukawa coupling comes back to O(1), thereby suppressing its negative contribution to m h . Therefore the tan β → ∞ scenario allows to choose third generation sfermion masses which lead to a correct Higgs mass.
Owing to decoupling of the multi-TeV extra Higgs masses required by the B-physics constraints, the MSSM Higgs sector is approximately SM-like [74]. Nonetheless it is instructive to qualitatively discuss the most relevant SM-like Higgs decay modes, h → bb, h → τ + τ − , and h → µ + µ − , which might potentially be affected by large bottom and charged lepton Yukawa couplings if the heavy Higgs masses are not sufficiently high to apply the decoupling limit. Their decay rates can be kept close to the SM values in the following way. The effective H 0 u -f -f coupling, for f = b, τ, µ, is generated by the same loop diagrams that generate m f . Therefore, the h-f -f coupling should be SM-like if h consists purely of H 0 u to a sufficient extent. Remember that the H 0 d -f -f coupling, y f , can be comparable to or even larger than y t . The remaining question is then how to maintain the Higgs mixing angle α small enough for h to avoid excessive H 0 d -contamination. It is easy to see that the 2 × 2 tree-level CP -even Higgs mass matrix becomes diagonal as tan β → ∞ (see e.g. eq. (23) of ref. [75]). The one-loop corrections to its off-diagonal elements are suppressed for small A t and A b , which can follow from gauge-mediated SUSY breaking at the same time as the small Bµ, needed for large tan β. Regarding h → γγ, we can use the fact that the stau-loop contribution (see e.g. refs. [76][77][78]) to the dimension-five effective operator decouples faster than the tau self energy, as the stau masses increase.

Dominance transition and sign change
In the Introduction we have shown that if all SUSY masses are equal, a SUSY µ becomes negative in the limit of tan β → ∞. In this case the wino contributions dominate in both a red µ and ∆ red µ , i.e. in the numerator and denominator of eq. (2.5), and they have opposite signs, see eqs. (2.12), (2.13). However, these equations also show that the resulting a SUSY µ will be positive if any of the three bino contributions dominates.
In the following we will delineate regions in the five-dimensional parameter space of the masses M 1,2 , µ, m L,R in which a SUSY µ is positive. We start by noting that the signs  only depend on mass ratios and that each of the three bino contributions can be expected to dominate if a particular mass hierarchy is valid: • TheBμ LμR contributions dominate for M 1 , m L , m R µ: "large µ-limit" • TheBHμ R contributions dominate for M 1 , µ, m R m L : "μ R -dominance" • TheBHμ L contributions dominate for M 1 , µ, m L |M 2 |, m R . These have the smallest numerical prefactors, and we can expect them to dominate only for rather extreme hierarchies. 3 The phenomenological behaviour and importance of theBμ LμR andBHμ R contributions have also been discussed recently in various contexts in refs. [43-47, 54, 65] and refs. [46,47,65,79], respectively. We first focus on cases where these two contributions become dominant. Parameter regions in which theBHμ L contributions dominate will be presented later on. Figure 2 illustrates plots of the individual contributions to a red µ and ∆ red µ for parameters which interpolate between the equal SUSY mass case and the two mass hierarchies for either "large µ-limit", figure 2(a), or "μ R -dominance", figure 2(b). In each hierarchy case the desired bino contributions become the largest ones for sufficiently large mass ratios. We first discuss figure 2(a), corresponding to the "large µ-limit", in detail. It shows the JHEP10(2015)026 contributions as functions of the ratio |M 2 | = µ over M 1 = m L = m R ; the horizontal axis is given by x = log 10 (µ/m R ). We pay particular attention to the dominance transition among the contributions and the sign change of a red µ and ∆ red µ .
• Around x = 0, where all masses are equal, the wino contribution dominates in both a red µ and ∆ red µ . For larger values of x, theBμ LμR contributions increase proportionally to µ whereas the other contributions are suppressed as 1/µ. Hence the latter tend to zero for large x, except for theWHL contribution to ∆ red µ , which tends to a constant because it is also proportional to M 2 and |M 2 | = µ in this plot.
• At x ≈ 0.36, the dominant contribution to a red µ changes from the wino contribution to theBμ LμR contribution. For larger mass ratios, x ≈ 1.00, the same dominance change also happens in ∆ red µ and thus a SUSY µ is always positive for x > 1. • For x < 0, the wino contributions remain dominant.
Similarly, we discuss figure 2(b), which corresponds to "μ R -dominance". It shows the contributions as functions of the ratio m L over M 1 = |M 2 | = µ = m R ; the horizontal axis is given by x = log 10 (m L /m R ). To investigate the sign of a SUSY µ comprehensively we plot it in a more general parameter space, as a function of the two mass ratios used in figures 2(b) and (a), i.e. in the plane m L /m R versus µ/m R . Figure 3 displays the results for four different choices of the remaining mass ratios M 1 /m R and |M 2 |/µ. We begin by discussing figure 3(a). It generalizes figure 2 by assuming M 1 = m R and |M 2 | = µ, so the different regions in figure 3(a) can be fully understood from the previous discussion.
• The large, central red/blue overlap region contains the origin at which all relevant SUSY masses are equal. In this region the wino contributions dominate and a SUSY µ is always negative.
• The two white regions, where a SUSY µ is always positive, correspond to the "large µ-limit" and the "μ R -dominance" region.     In figure 3(a), the gaugino masses |M 2 | and M 1 are equated with µ and m R , respectively. This is a reasonable choice since their precise values are qualitatively unimportant for the appearance of the different regions. The dependence on these two gaugino masses is investigated in figures 3(b)-(d), whose presentation is similar to figure 3(a) but with either M 1 = 10 m R , |M 2 | = 10 µ, or M 1 = 0.1 m R and |M 2 | = 100 µ. Together, these plots cover large, representative parts of the full parameter space of the four relevant mass ratios.
Obviously, all plots can be understood in the same way as figure 3(a). As M 1 becomes heavier than m R in figure 3(b), the bino contributions are suppressed, the wino-dominance area expands, and the positive a SUSY  where M SUSY,min is the smallest among the five relevant mass parameters M 1 , |M 2 |, µ, and m L,R [and thus similar but not equal to the lightest SUSY particle (LSP) mass since the latter arises from diagonalization of mass matrices]. The coefficient C is expected to be of order unity because of radiative muon mass generation. It depends only on the four independent mass ratios (up to terms suppressed by powers of M 2 Z /M 2 SUSY ). Hence for each given set of mass ratios we can pose two questions: 1. What is the value of C?
The answers are given in figure 5 for the special choices M 1 = m R and |M 2 | = µ which have already been used in figure 3(a). The figures also show the behaviour of the muon Yukawa coupling, which is proportional to 1/∆ red µ , see eq. (2.6). The four most interesting regions are the "large µ-limit" and "μ R -dominance" regions with either positive or negative M 2 . These are located in the right end and middle top areas of the two plots and correspond to the white regions of figure 3(a) with a positive a SUSY µ . In detail figure 5 shows the following.
• The first important observation is that in the centre of the "large µ-limit" and "μ Rdominance" regions, C lies between 0.2 and 0.4; the corresponding minimum SUSY mass is around 1 TeV. C does not quite reach unity because of the behaviour of the loop functions in ∆ red µ and a red µ .
• In both regions, the magnitude of C, the minimum SUSY mass, and the muon Yukawa coupling depend on the sign of M 2 . This dependence arises from interference of the dominant contributions with subdominant wino contributions, which can be either constructive in ∆ red µ and destructive in a red µ or vice versa.
• There are small white strips in parameter space with larger C and minimum SUSY masses. In these strips ∆ red µ undergoes a sign change and the Yukawa coupling becomes infinite. Here perturbation theory is not trustworthy since two-loop effects are non-negligible due to the accidental cancellation of the one-loop contributions. We exclude these regions from further discussion.
• Of particular interest are the regions with smaller Yukawa couplings, i.e. the "large µ-limit" region for negative M 2 and the "μ R -dominance" region for positive M 2 .
Here C can still increase up to 0.3, and the minimum SUSY mass reaches up to 1.1 TeV. The opposite sign choices lead to the strips where y µ diverges. In the "large µ-limit" the muon Yukawa coupling is always proportional to 1/µ and therefore its magnitude is comparatively small and decreasing as µ increases. In the regions of "μ R -dominance" |y µ | is generally much larger. Eq. (2.12) allows to deduce that for "μ R -dominance", |y µ | cannot become smaller than 0.78 for negative M 2 and 0.12 for positive M 2 , respectively.  In the white regions the denominator ∆ red µ undergoes a sign change, and the perturbation theory becomes untrustworthy due to large y µ . Benchmark parameter points introduced in table 1 are marked with blue squares. Now, we consider constraints on the relevant parameter space for large a SUSY µ . All previous plots are based on the mass-insertion approximation and depend only on mass ratios. For each given set of mass ratios with positive a SUSY µ there is a value of the overall mass scale M SUSY for which a SUSY µ agrees with experiment. In the following discussions we always fix M SUSY to that particular scale. Figure 6 shows the same parameter space as figure 5, but displays regions in which charged SUSY masses become too small (medium grey), the muon Yukawa coupling is non-perturbative (dark grey) or the electroweak vacuum becomes unstable (hatched). In detail, the constraints are the following:  which are far above the LHC limits. Therefore we show as the medium grey regions in figure 6 an exemplary contour corresponding to a chargino or smuon which is lighter than 100 GeV.
• In the dark grey regions the muon Yukawa coupling violates perturbativity, |y µ | > √ 4π. Generally, we regard the model as a pure low-energy model and do not require the stronger constraint that there be no Landau poles at higher scales; such an analysis has been done in refs. [37,38]. It should be further noted that the observable couplings of the lightest Higgs boson are SM-like up to corrections suppressed by powers of the SUSY scale. This is true independently of the value of the fundamental Yukawa coupling due to decoupling.
• As discussed in refs. [44,[76][77][78]80], the combination y µ is limited by the requirement that the electroweak vacuum is metastable with a lifetime longer than the age of the universe. This constrains in particular the parameter regions with large µ as already stressed in ref. [44]. The regions which are excluded by the fitting formula from this reference applied to the muon sector are shown in figure 6. Since it is possible that y ∼ > y t in our scenario, the O(y n ) counterpart of the standard O(y n t ) correction to JHEP10(2015)026 the effective potential might be significant. To scrutinize the effect of the former, we add to the effective potential the one-loop contributions (see e.g. eq. (1.2) of ref. [81]) in the gaugeless limit, with the renormalization scale Q = √ m L m R . They include the (s)lepton-and the Higgs(ino)-as well as the (s)top-loop contributions. We then evaluate the Euclidean action numerically using the method presented in refs. [82,83]. In principle, the metastability bound on y µ depends on the pseudoscalar Higgs mass M A affecting the one-loop correction to the slepton squared mass. Motivated by the B-physics constraints, we set M A = 50 TeV, a high enough value for the extra Higgs states to decouple. By applying the specified procedure to the muon sector, we find that the regions shown in figure 6 are excluded because of an unstable electroweak vacuum.
According to the results shown in figure 6, the requirement of vacuum metastability poses strong constraints on the parameter space, in particular on the "large µ-limit" region. This region is completely excluded for positive M 2 , while for negative M 2 the ratio µ/m R must not be too large. The "μ R -dominance" region is not constrained by vacuum metastability.
The two different approaches to the vacuum metastability yield similar results but the influence of the additional terms is visible in figure 6(b). The Yukawa coupling of the τ lepton, which is another quantity of interest, can be determined similarly to y µ . There is a strong correlation if the stau masses are assumed to be equal to the smuon masses; without this assumption the two Yukawa couplings are basically independent, and vacuum metastability from the tau sector does not further constrain the five-dimensional parameter space. However, a detailed analysis should take into account higher orders in the large tau Yukawa coupling and is beyond the scope of this paper.
Thus we find that out of the four interesting regions in figure 5, three survive. Now we briefly comment on the comparison between the approximate results for y µ as well as a SUSY µ obtained from eqs. (2.7), (2.9) and their exact results according to the appendix. In the three viable regions described above, the deviation amounts to at most a few percent. Outside these regions the approximation becomes worse for very small SUSY masses and in the small strips with divergent muon Yukawa coupling.

Conclusions
We have studied the muon magnetic moment in the MSSM for tan β → ∞. This is a viable limit in which all down-type masses are generated by loop-induced couplings to the "wrong" Higgs H u . The scenario can also be viewed independently of supersymmetry, as a simplified model which realizes the general idea of radiative muon mass generation.
The SUSY contribution to the muon anomalous magnetic moment is given by eq. (2.5) as the ratio of two one-loop quantities. For this reason the scenario leads to large a µ with TeV-scale new physics masses, and it behaves qualitatively very differently from the standard MSSM with tan β ∼   a SUSY µ depends on the ratios between the five relevant mass parameters; figures 2-4 give a comprehensive explanation of which mass ratios lead to positive a SUSY µ .
The three most promising parameter regions with positive a SUSY µ are the "large µlimit" for negative M 2 and the "μ R -dominance" region for positive or negative M 2 . The "large µ-limit" region with positive M 2 is excluded because of vacuum instability. In all three viable cases, large contributions to a SUSY µ are possible without conflicts with other experimental constraints, see figures 5, 6. The discrepancy (1.1) can be fully explained with a lightest SUSY mass above 1 TeV. If this scenario is correct, the LHC would not be able to find the SUSY particles relevant for a µ . However, section 2.2 shows that the scenario could potentially influence and thus be tested in a multitude of observables in B-physics and Higgs physics.
To conclude, we summarize the behaviour in three simple ways. Table 1 lists several benchmark parameter points with SUSY mass parameters at 1 TeV or higher and gives the resulting a SUSY µ and muon Yukawa coupling. Two points correspond to the "large µ-limit" with negative M 2 and three to the "μ R -dominance" region with positive or negative M 2 . Some points involve only one heavy and four lighter SUSY masses; for these large a SUSY µ is possible for less extreme mass splitting than for the ones with two heavy and three lighter masses. The benchmark points further illustrate that the Yukawa coupling is much larger in the case of "μ R -dominance". Figure 7 shows the results of two scans of the five-dimensional parameter space: plotted are the largest possible values of the minimum SUSY mass M SUSY,min which lead to agreement with eq. (1.1). The left panel shows the results as a function of log 10 (m L /M 1 ) (for positive M 2 ), it corresponds to the y-axis of figure 5(a); the right panel as a function of log 10 (µ/M 1 ) (for negative M 2 ), it corresponds to the x-axis of figure 5(b). The results correspond to the regions presented in figure 5 and are just slightly higher than the values obtained there with fixed mass ratios.
Finally, for the regions of interest, a simple approximation in the style of eq. (1.7) can be derived: for the "large µ-limit" region, we consider the limit |µ| |M 1 | = m L = m R ≡ M SUSY , and for the "μ R -dominance" region we consider the limit m L |µ| = |M 1 | = m R ≡ M SUSY . In both of these limits, we obtain In this appendix we provide the full one-loop results for the muon self energy and the magnetic moment, defined in eqs. (2.2), (2.4). We use the notation of ref. [47], which is similar to the ones from refs. [40,41]. Note that in the standard case of the MSSM with moderate tan β ∼ < 50, the combination m µ tan β frequently appears. An important difference is that here we have to take into account the distinction between the tree-level and the physical muon mass, m tree µ = m µ , and we have to use the identification (2.3), which is exactly valid. Here g 2 denotes the SU(2) gauge coupling. The relevant interaction terms for chargino-sneutrino and neutralino-smuon interactions with muons can be written as L int =χ − i c L iνµ P L + c R iνµ P R µν † µ +χ 0 j n L jμ k P L + n R jμ k P R µμ † k + h. c. , (A.2)

JHEP10(2015)026
where the coupling coefficients are defined as The unitary matrices U, V, N diagonalize the chargino and neutralino mass matrices; Uμ is used to diagonalize the smuon mass matrix. The resulting smuon mass eigenvalues are given by A term containing the trilinear coupling vanishes because of v d = 0, and in the last term we have used eq. (A.1). The occurring SUSY parameters are the higgsino mass parameter µ, the gaugino masses M 1 as well as M 2 and the soft mass parameters for the secondgeneration slepton doublet and singlet m L and m R , respectively. With this notation, the MSSM one-loop contributions to the on-shell muon self-energy from chargino-sneutrino or neutralino-smuon loops can be expressed as Σ MSSM µ = Σχ containing the Euler-Mascheroni constant γ E . Owing to the unitarity of the mixing matrices, in eqs. (A.5) the UV divergences as well as the dependence on µ R vanish. Because of v d = 0, the tree-level mass of the muon (as well as the other charged leptons and down-type quarks) equals zero. Therefore, the entire muon mass has to be generated at the loop level, which means m µ = −Σ MSSM µ . (A.8) For a given parameter point, this equation can be used to determine the muon Yukawa coupling y µ . Note that the self energy is approximately linear in y µ , as suggested by the definition of ∆ red µ in eq. (2.2). But the Yukawa coupling also enters nonlinearly, through quadratic terms in n L jμ k (n R jμ k ) * , and via the smuon masses (A.4). Therefore, it has to be determined numerically in an exact calculation. As shown in section 2, however, the nonlinear terms vanish in the mass-insertion approximation.
The SUSY one-loop diagrams contributing to the anomalous magnetic moment of the muon are the same as for the muon self-energy, except for the fact that now a photon couples to the charged SUSY particle in the loop, i.e. the chargino or the smuon. Again neglecting terms suppressed by tan β, the result reads a SUSY  [40] as normalized such that F C 2 (1) = F N 2 (1) = 1. The formulas look identical to the respective terms in the standard case with moderate tan β, but like for the self energy, there are nonlinear terms in y µ , which usually are neglected because of the m µ /M SUSY suppression. In our case, these terms are only suppressed as M W,Z /M SUSY ; hence they are much more significant in an exact calculation. Nevertheless, like for the self energy, in the massinsertion approximation, a SUSY µ is linear in the Yukawa coupling.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.