Improved determination of the Higgs mass in the MSSM with heavy superpartners

We present several advances in the effective field theory calculation of the Higgs mass in MSSM scenarios with heavy superparticles. In particular, we compute the dominant two-loop threshold corrections to the quartic Higgs coupling for generic values of the relevant SUSY-breaking parameters, including all contributions controlled by the strong gauge coupling and by the third-family Yukawa couplings. We also study the effects of a representative subset of dimension-six operators in the effective theory valid below the SUSY scale. Our results will allow for an improved determination of the Higgs mass and of the associated theoretical uncertainty.


Introduction
At the price of doubling the particle content of the Standard Model (SM), supersymmetry (SUSY) provides elegant solutions to several open issues, including the stability of the electroweak (EW) scale, the nature of dark matter and the possibility of embedding the SM in a grand-unified gauge theory. Common features of supersymmetric extensions of the SM are an extended Higgs sector and the existence of tree-level relations between the quartic Higgs couplings and the other couplings of the considered model, which translate into predictions for the Higgs-boson masses. When radiative corrections are included, those predictions are sensitive to the whole particle spectrum of the model, and can be used to constrain its parameter space even before the discovery of SUSY particles.
For the MSSM, both the discovery in 2012 [40,41] of a SM-like Higgs boson with mass about 125 GeV [42] and the negative results of the searches for stop squarks at the LHC [43][44][45][46][47][48] favor scenarios with a SUSY mass scale M S in the TeV range. In particular, the observed value of the Higgs mass requires the radiative correction to the squared-mass parameter, m 2 h , to be at least as large as its tree-level value: if the stops are heavy enough, this can be realized via the dominant top/stop contributions, which are controlled by the top Yukawa coupling, g t ∼ O (1), and are enhanced by logarithms of the ratio between the stop and top masses. A further increase in m 2 h can be obtained if the left-right stop mixing parameter X t is about twice the average stop mass. Roughly speaking, for tan β large enough to almost saturate the treelevel bound on the lightest-scalar mass, m h ≈ 125 GeV requires the average stop mass to be somewhere around 1 TeV for the "maximal" (i.e., most favorable) value of X t , and above 10 TeV for vanishing X t . However, when the SUSY scale is significantly larger than the EW scale, fixed-order calculations of m h such as the ones implemented in the codes mentioned above may become inadequate, because radiative corrections of order n in the loop expansion contain terms enhanced by as much as ln n (M S /m t )-where we take the top mass as a proxy for the EW scale. Indeed, a possible symptom of such heavy-SUSY malaise is the fact that, in scenarios with TeV-scale stop masses and large stop mixing, the spread in the predictions of those codes for m h exceeds the theoretical accuracy of their (largely equivalent) two-loop calculations, which was estimated in the early 2000s to be about 3 GeV [49,50] in what were then considered natural regions of the MSSM parameter space.
In the presence of a significant hierarchy between the SUSY scale and the EW scale, the computation of the Higgs mass needs to be reorganized in an effective field theory (EFT) approach: the heavy particles are integrated out at the scale M S , where they only affect the matching conditions for the couplings of the EFT valid below M S ; the appropriate renormalization-group equations (RGEs) are then used to evolve those couplings between the SUSY scale and the EW scale, where the running couplings are related to physical observables such as the Higgs-boson mass and the masses of gauge bosons and fermions. In this approach, the computation is free of large logarithmic terms both at the SUSY scale and at the EW scale, while the effect of those terms is accounted for to all orders in the loop expansion by the evolution of the couplings between the two scales. More precisely, large corrections can be resummed to the (next-to) nleading-logarithmic (N n LL) order by means of n-loop calculations at the SUSY and EW scales combined with (n +1)loop RGEs. On the other hand, the common procedure of matching the MSSM to a renormalizable EFT -such as the plain SM -in the unbroken phase of the EW symmetry amounts to neglecting corrections suppressed by powers of v 2 /M 2 S , where we denote by v the vev of a SM-like Higgs scalar. Those corrections can in fact be mapped to the effect of non-renormalizable, higher-dimensional operators in the EFT Lagrangian.
The EFT approach to the computation of the MSSM Higgs mass dates back to the early 1990s [51][52][53]. Over the years, it has also been exploited to determine analytically the coefficients of the logarithmic terms in m 2 h at one [54], two [55][56][57][58] and even three [49,59] loops, by solving perturbatively the appropriate systems of boundary conditions and RGEs. However, when the focus was on "natural" scenarios with SUSY masses of a few hundred GeV, the omission of O(v 2 /M 2 S ) terms limited the accuracy of the EFT approach, and the effect of the resummation of logarithmic corrections was not deemed important enough to justify abandoning the fixed-order calculations of the Higgs mass in favor of a complicated EFT set-up with higher-dimensional operators. 2 More recently, an interest in "unnatural" scenarios such as split SUSY [61,62] and high-scale SUSY (see, e.g., Ref. [63]), and then the LHC results pushing the expectations for the SUSY scale into the TeV range, have brought the EFT approach back into fashion. On the one hand, in Ref. [64] the authors of FeynHiggs combined the fixedorder calculation of m h implemented in their code with a resummation of the LL and NLL terms controlled exclusively by g t and by the strong gauge coupling g 3 . On the other hand, three papers [65][66][67] presented updates of the traditional EFT calculation: the use of the state-of-the-art results collected in Ref. [68] for the SM part (i.e., three-loop RGEs and two-loop EW-scale matching conditions), together with the full one-loop and partial two-loop matching conditions at the SUSY scale, allow for a full NLL and partial NNLL resummation of the logarithmic corrections. 3 Several public codes for the EFT calculation of the Higgs mass in the MSSM with heavy SUSY have also been released: SusyHD [69], based on Ref. [67]; MhEFT [70], based on Refs. [65,71] and covering as well scenarios with a light two-Higgs-doublet model (THDM); HSSUSY [72,73], a module of FlexibleSUSY [74] with the same essential features as the original SusyHD; FlexibleEFTHiggs [72,73], which combines a full one-loop computation of m h with a LL resummation of the logarithmic corrections; finally, an EFT approach similar to the one of Ref. [73] was recently implemented in SPheno/SARAH [75].
In MSSM scenarios with stop masses of several TeV, where the effects of O(v 2 /M 2 S ) can be safely neglected, the theoretical uncertainty of the EFT prediction for the Higgs mass stems from missing terms of higher orders in the loop expansion, both in the calculation of the matching conditions at the SUSY scale and in the SM part of the calculation. In Refs. [66,67] such uncertainty was estimated to be at most 1 GeV in a simplified MSSM scenario with degenerate SUSY masses of 10 TeV, tan β = 20 and vanishing X t , where m h ≈ 123.5 GeV. In such scenario, the prediction for m h of the "hybrid" (i.e., fixed-order+partial NLL) calculation of Ref. [64] was about 3 GeV higher, well outside the theoretical uncertainty of the EFT result. In Refs. [67,71] it was suggested that most of the discrepancy came from the determination of the coupling g t used in the resummation procedure, for which Ref. [64] omitted one-loop EW and two-loop QCD effects, consistently with the accuracy of the m h calculation in that paper. Those effects were later included in FeynHiggs, which now also allows for a full NLL and partial NNLL resummation of the logarithmic corrections [76]. In the simplified MSSM scenario mentioned above, the refinements in the resummation procedure of FeynHiggs reduce the discrepancy with the EFT prediction for m h to a few hundred MeV.
As mentioned earlier, MSSM scenarios with stop masses below a couple of TeV and large stop mixing -which are definitely more interesting from the point of view of LHC phenomenology -suffer from even larger spreads in the predictions of different codes for m h . For example, in a benchmark point with degenerate SUSY masses of 1 TeV, tan β = 20, and X t varied so as to maximize m h , the EFT calculation finds m max h ≈ 123 GeV, whereas SoftSusy, SuSpect and SPheno -which implement the same corrections to the Higgs masses, but differ in the determination of the running couplings -find m max h ≈ 124.5 − 126.5 GeV, and the latest version of FeynHiggs [77] finds m max h ≈ 126 − 128 GeV (depending on the code's settings). However, in this case the comparison between the EFT prediction for m h and the various fixed-order (or hybrid) predictions is less straightforward than in scenarios with multi-TeV stop masses, because there is no obvious argument to favor one calculational approach over the others: the O(v 2 /M 2 S ) terms might or might not be negligible, and the logarithmic corrections might or might not be important enough to mandate their resummation. For all approaches, this unsatisfactory situation points to two urgent needs: first, to improve the calculation of m h with the inclusion of higher-order effects; second, to provide a better estimate of the theoretical uncertainty, tailored to the "difficult" region of the parameter space with stop masses about 1-2 TeV.
In this paper we take several steps towards an improved EFT determination of the Higgs mass in the MSSM with heavy superpartners. In particular, in Sect. 2 we compute the two-loop, O(g 6 t ) contribution to the SUSY-scale matching condition for the quartic Higgs coupling -which was previously known only in simplified scenarios [21,65,67] allowing for generic values of all the relevant SUSY-breaking parameters. We also include the two-loop contributions controlled by the bottom and tau Yukawa couplings, addressing some subtleties related to the presence of potentially large tan β-enhanced corrections. Our new results bring the matching condition for the quartic Higgs coupling to the same level, in terms of an expansion in coupling constants, as the two-loop Higgs-mass calculations in SoftSusy, SuSpect and SPheno. In Sect. 3 we study instead the effects of a representative subset of dimension-six operators in the EFT. We obtain both an improvement in our prediction for m h in scenarios with stop masses about 1-2 TeV and a morerealistic estimate of the theoretical uncertainty associated to O(v 2 /M 2 S ) effects. The results presented in this paper have been implemented in modified versions of the codes SusyHD [69] and HSSUSY [72]. All the analytic formulas that proved too lengthy to be printed here are available upon request in electronic form.

Two-loop matching of the quartic Higgs coupling
In this section we describe our calculation of the two-loop matching condition for the quartic Higgs coupling. We consider a set-up in which all SUSY particles as well as a linear combination of the two Higgs doublets of the MSSM are integrated out at a common renormalization scale Q ≈ M S , so that the EFT valid below the matching scale is just the SM. Using the conventions outlined in Sect. 2 of Ref. [66], the two-loop matching condition for the quartic coupling of the SM-like Higgs doublet H takes the form where g and g are the EW gauge couplings, β can be interpreted as the angle that rotates the two original MSSM doublets into a light doublet H and a massive doublet A, and λ n is the n-loop threshold correction to the quartic coupling arising from integrating out the heavy particles at the scale M S . The contributions to λ 1 controlled by the EW gauge couplings and by the top Yukawa coupling, for generic values of all SUSY parameters, were given in Ref. [66], completing and correcting earlier results of Refs. [78,79]. For completeness, we report in the appendix the full result for the one-loop contributions of heavy scalars, including also terms controlled by the bottom and tau Yukawa couplings. However, the only one-loop contributions relevant to our computation of the two-loop threshold correction, where we will consider the "gaugeless" limit g = g = 0, are those proportional to the fourth power of a third-family Yukawa coupling, which read where byĝ f we denote SM-like Yukawa couplings, 4 related to their MSSM counterpartsŷ f byĝ t =ŷ t sin β ,ĝ b = y b cos β andĝ τ =ŷ τ cos β . Moreover, for each fermion species f : N f c is the number of colors; (mf L , mf R ) are the soft SUSY-breaking sfermion masses, i.e. (m Q 3 , m U 3 ), (m Q 3 , m D 3 ) and (m L 3 , m E 3 ) for stops, sbottoms and staus, respectively; is the trilinear soft SUSY-breaking Higgs-sfermion coupling, μ is the Higgs-mass term in the MSSM superpotential, r t = cot β and r b = r τ = tan β ; x f = mf L /mf R ; finally, the loop functions F 1 and F 2 are defined in Appendix A of Ref. [66].
For what concerns the two-loop threshold correction λ 2 , simplified results for the O(g 4 t g 2 3 ) and O(g 6 t ) contributions, valid in the limit m Q 3 = m U 3 = m A = mg ≡ M S (where m A is the mass of the heavy Higgs doublet and mg is the gluino mass), were made available as far back as in Ref. [21]. Among the recent EFT analyses, Refs. [66,67] obtained formulas for the O(g 4 t g 2 3 ) contributions valid for arbitrary values of all the relevant SUSYbreaking parameters. The O(g 6 t ) contributions, on the other hand, were neglected in Ref. [66], while they were included in Refs. [65,67] only through simplified formulas derived from those of Ref. [21]. In this paper we extend the calculations of Refs. [66,67] to obtain all contributions to λ 2 controlled only by the third-family Yukawa couplings, again for arbitrary values of all the relevant SUSY-breaking parameters. Besides improving our knowledge of the O(g 6 t ) contributions from two-loop diagrams involving stops, this allows us to properly account for sbottom and stau contributions that can become relevant at large values of tan β. We also discuss how to obtain the O(g 4 b g 2 3 ) contributions from the known results for the O(g 4 t g 2 3 ) ones. Altogether, our results amount to a complete determination of λ 2 in the limit of vanishing EW gauge (and first-two-generation Yukawa) couplings.

Outline of the calculation
The two-loop, Yukawa-induced threshold correction to the quartic Higgs coupling λ at the matching scale Q can be expressed as where by V 2 , heavy we denote the contribution to the MSSM scalar potential from two-loop diagrams involving sfermions that interact with themselves, with Higgs doublets or with matter fermions and higgsinos only through the thirdfamily Yukawa couplings, as well as from two-loop diagrams involving only the heavy Higgs doublet and matter fermions. The last two terms in Eq. (3) contain additional two-loop contributions that will be described below. In the following we will focus on the contributions to λ 2 that involve the top and bottom Yukawa couplings, and comment only briefly on the inclusion of the contributions that involve the tau Yukawa coupling, which are in general much smaller. However, we stress that the results that we implemented in SusyHD [69] and HSSUSY [72] (and that we make available upon request) do include the tau-Yukawa contributions through two loops.
In the gaugeless limit adopted in our calculation, the fielddependent mass spectrum of the particles that enter the relevant two-loop diagrams simplifies considerably: we can approximate the masses of the lightest Higgs scalar and of the would-be Goldstone bosons to zero, and the masses of all components (scalar, pseudoscalar and charged) of the heavy Higgs doublet to m 2 A ; the charged and neutral components of the two higgsino doublets combine into Dirac spinors with degenerate mass eigenvalues |μ| 2 ; the tree-level mixing angle in the CP-even sector is just α = β − π/2. For the contributions to V 2 , heavy that involve the top and bottom Yukawa couplings, we adapt the results used for the effective-potential calculation of the MSSM Higgs masses in Ref. [25]. 5 To compute the fourth derivative of the effective potential entering Eq. (3) we follow the approach outlined in section 2.3 of Ref. [66]: we express the stop and sbottom masses and mixing angles as functions of field-dependent top and bottom masses, m t =ĝ t |H | and m b =ĝ b |H |, and obtain where the term in the last line is obtained from the terms in the first three lines by swapping top and bottom, and we used the shortcuts The derivatives of the field-dependent stop and sbottom parameters and the limit of vanishing top and bottom masses in Eq. (4) are obtained as described in Ref. [66]. As in the case of the O(g 4 t g 2 3 ) contributions, we find that the fourth derivative of the two-loop effective potential contains terms proportional to ln(m 2 q /Q 2 ), which would diverge for vanishing quark masses but cancel out against similar terms in the contribution denoted as λ shift, f in Eq. (3). Indeed, above the matching scale the one-loop contribution to the quartic Higgs coupling from box diagrams with a top or bottom quark, is expressed in terms of the MSSM couplingsĝ q , whereas below the matching scale the same contribution is expressed in terms of the SM couplings g q . To properly compute the two-loop, Yukawa-only part of the matching condition for the quartic Higgs coupling, we must re-express the MSSM couplings entering δλ g 4 q , q above the matching scale (including those implicit in m q ) according toĝ q → g q (1 + g Y q ), where g Y q denotes the terms controlled by the Yukawa couplings in the threshold correction to g q . In particular, we find for the top and bottom Yukawa couplings where the last term on the right-hand side of each equation reads, in a notation analogous to the one of Eq. (2), and corresponds to the threshold correction to the light-Higgs WFR arising from squark loops. The loop functions F 5 , F 6 and F 9 are defined in appendix A of Ref. [66]. We also remark that Eqs. (7)-(9) assume that the angle β entering the relations between the SM-like couplingsĝ q and their MSSM counterpartsŷ q is renormalized as described in Sect. 2.2 of Ref. [66], removing entirely the contributions of the off-diagonal WFR of the Higgs doublets. Combining the effects of the shifts in the Yukawa couplings with the renormalization of the Higgs fields (keeping into account also the field-dependent quark masses in the logarithms) we obtain the total contribution to λ 2 arising from the quark-box diagrams of Eq. (6), which cancels the logarithmic dependence on the quark masses of the derivatives of V 2 , heavy . We checked that the contributions in Eq. (4) that involve more than two derivatives of the two-loop effective potential cancel out completely against the shift of the corresponding contributions in the one-loop part -namely, the non-logarithmic term in the righthand side of Eq. (6) -so that the final result for λ 2 can be related to the two-loop correction to the light-Higgs mass. This is the same "decoupling" property found in Ref. [66] for the O(g 4 t g 2 3 ) part of λ 2 . Finally, it can be inferred from Eqs. (7)-(10) that the contribution of δ Zq H cancels out of λ shift, f .
The last term in Eq. (3), λ shift,f , arises from shifts in the sfermion contribution to the one-loop matching condition for the quartic Higgs coupling, Eq. (2). In particular, it contains terms arising from the WFR of the Higgs fields, which are not captured by the derivatives of V 2 , heavy , plus additional contributions that arise if we express the one-loop threshold correction in Eq. (2) in terms of the SM Yukawa couplings, g q , instead of the MSSM ones,ĝ q . We remark here that, while the shift of the Yukawa couplings in the quark-box diagrams of Eq. (6) is required for a consistent two-loop matching of the quartic Higgs coupling, an analogous shift in the squark contribution of Eq. (2) is to some extent a matter of choice. In Refs. [66,67] the top Yukawa coupling entering the one-loop part of the threshold correction to the quartic Higgs coupling was interpreted as the SM one. Applying that choice to both the top and the bottom Yukawa couplings, we would find where again the contributions of the WFR of the Higgs fields cancel out against analogous terms in the shifts of the Yukawa couplings. After including in λ 2 the shifts in Eqs. (10) and (11), we checked that, in the limit of part of λ 2 coincides with the simplified result given in Eq. (21) of Ref. [67].
On the other hand, it is well known [80][81][82] that the relation between the bottom Yukawa coupling of the SM and its MSSM counterpart is subject to potentially large corrections enhanced by tan β, which, in the gaugeless limit, arise from diagrams involving either gluino-sbottom or higgsino-stop loops. As discussed, e.g., in Ref. [83], these tan β-enhanced terms can be "resummed" in the DR-renormalized coupling of the MSSM by expressing it aŝ where g b (Q) is the MS-renormalized coupling of the SM, extracted at the EW scale from the bottom mass and evolved up to the matching scale Q with SM RGEs, g Y b is given in Eq. (7), and where C F = 4/3 is a color factor, and we recall that In contrast with our treatment of the top Yukawa coupling, we will therefore choose to interpret the bottom Yukawa coupling entering the one-loop part of the threshold correction to the quartic Higgs coupling as the MSSM one, in order to absorb the tan β-enhanced effects directly in λ 1 . We recall that a similar approach was discussed in Refs. [18,24,25] in the context of the fixed-order calculation of the Higgs masses in the MSSM. With our choice for the bottom Yukawa coupling entering λ 1 , we must omit the term 4 g Y b in the formula for λ shift,f , see Eq. (11), when computing the contributions to λ 2 controlled only by the top and bottom Yukawa couplings. Concerning the O(g 4 b g 2 3 ) contributions, they can be obtained from the O(g 4 t g 2 3 ) contributions computed in Refs. [66,67] via where the notation [ t → b ] in the first line represents the contributions. We note that, in practice, our choice removes from λ 2 potentially large terms characterized by a higher power of tan β than ofĝ b , i.e. terms scaling likeĝ 4 b g 2 3 tan 5 β or likeĝ 4 bĝ 2 t tan 5 β. We now comment on the inclusion of the contributions to λ 2 controlled by the tau Yukawa coupling. The twoloop contributions of O(g 6 τ ), i.e. those involving only the tau Yukawa coupling, do not require a separate calculation, since they can be obtained from the top-only, O(g 6 t ) ones via the 3 and cos β ↔ sin β (see also Ref. [25]). Indeed, as long as we neglect the EW gauge couplings, the threshold correction to the tau Yukawa coupling does not contain any tan β-enhanced terms, and reads where the sfermion contribution to the Higgs WFR, δ Zf H , is obtained including also the stau contribution (with N c = 1) in Eq. (9). We can therefore treat the tau Yukawa coupling in the same way as the top one, expressing the stau contribution to λ 1 in terms of the SM coupling g τ . In addition, "mixed" contributions to the two-loop effective potential controlled by both the tau and the bottom Yukawa couplings arise from diagrams that involve the quartic sbottom-stau coupling; see Appendix B of Ref. [50]. The corresponding contributions to λ 2 can be obtained directly from the derivatives of the effective potential (without additional shifts) with the procedure outlined around Eq. (4), after replacing t → τ in the latter. Finally, the choice of using the MSSM couplinĝ g b in λ 1 spoils the cancellation of Higgs WFR effects in λ shift,f ; see Eq. (11). As a result, when we take into account the O(g 2 τ ) contribution from stau loops in δ Zf H , we find an additional O(g 4 b g 2 τ ) contribution to λ 2 .

Numerical examples
We now provide some illustration of the numerical impact of the newly computed two-loop corrections to the quartic Higgs coupling. To this purpose, we implemented those corrections in modified versions of the codes SusyHD [69] and HSSUSY [72]. All plots presented in this section were One loop To start with, we omit all contributions to λ 2 controlled by the bottom and tau Yukawa couplings, and focus on the effect of extending the contributions controlled by the top Yukawa coupling to generic values of the relevant SUSY-breaking parameters. We consider a scenario in which all SUSY-particle masses are larger than one TeV, but the stop masses are not degenerate. In particular, we take m U 3 = 1.5 TeV and m Q 3 = κ m U 3 , where κ is a scaling parameter that we vary in the range 1 ≤ κ ≤ 4. We also take mg = m A = m U 3 , μ = 4 m U 3 and tan β = 20, and we fix A t via the "maximal" stop mixing condition A t − μ cot β = (6 m Q 3 m U 3 ) 1/2 . For the remaining MSSM parameters, which affect the one-loop part of the calculation, we set all sfermion masses other than those of the stops, as well as the EW gaugino masses, equal to m U 3 , and we take A b = A τ = A t . All of the MSSM parameters listed abovewith the exception of tan β, which is defined as described in Sect. 2.2 of Ref. [66] -are interpreted as DR-renormalized parameters at the scale Q = (m Q 3 m U 3 ) 1/2 .
In Fig. 1 we compare the predictions for m h obtained with the "exact" (i.e., valid for generic SUSY masses) formulas for the top-Yukawa contributions to λ 2 with "approximate" predictions obtained by replacing the scalar and gluino masses of our scenario with the degenerate masses and O(g 6 t ) contributions to λ 2 the simplified formulas given in Refs. [66,67], respectively. In particular, the dotted black line in the left plot of Fig. 1 represents the prediction for m h , as a function of the stop mass ratio κ = m Q 3 /m U 3 , obtained by neglecting all two-loop contributions to the matching of the quartic Higgs coupling, and using the exact results from Refs. [66,67] for the one-loop contributions; the dashed blue line includes also the simpli-  Figure 1 confirms that, as already noticed in Refs. [66,67], the overall effect of the top-Yukawa contributions to λ 2 on the EFT predictions for m h in scenarios with multi-TeV stop masses is rather small, typically less than one GeV. However, the comparison between the dashed and solid lines in the plots of Fig. 1 The meaning of the different curves is explained in the text spectra, the use of simplified formulas with an "average" SUSY mass can lead to a rather poor approximation of the exact results. In particular, the comparison between dashed and solid blue lines shows that by using Eq. (36) of Ref. [66] for the O(g 4 t g 2 3 ) corrections we would significantly overestimate their effect on m h when κ 2 in our scenario. In turn, the dashed and solid red lines show that, by using Eq. (21) of Ref. [67] for the O(g 6 t ) corrections, we could entirely mischaracterize their effect on the Higgs mass: between the point where the solid blue line crosses the solid red one and the point where it crosses the dashed red one, the approximate calculation of the O(g 6 t ) corrections gives a negative shift in m h , while the exact calculation gives a positive shift. We remark, however, that the latter finding depends on the somewhat large value of μ adopted in our scenario: for smaller μ the quality of the approximation for the O(g 6 t ) corrections would improve.
We now turn our attention to the effect of the threshold corrections to the quartic Higgs coupling controlled by the bottom Yukawa coupling. In Fig. 2 we show the EFT prediction for m h as a function of tan β, in a simplified MSSM scenario with all soft SUSY-breaking masses of sfermions and EW gauginos, as well as the heavy Higgs-doublet mass m A , set equal to M S = 1.5 TeV, while the gluino mass is set to mg = 2.5 TeV; the trilinear Higgs-stop coupling A t is fixed by the maximal mixing condition A t − μ cot β = √ 6 M S , and A b = A τ = A t ; finally, we take μ = −1.5 TeV, to enhance the effect of the corrections controlled by the bottom Yukawa coupling. Indeed, negative values of the products μ mg and μA t ensure thatĝ b -which we extract at the matching scale from the SM coupling g b via Eq. (12) -becomes larger for increasing tan β, and possibly hits a pole as the denominator on the right-hand side of Eq. (12) approaches zero. Again, all soft SUSY-breaking parameters as well as μ are renormalized in the DR scheme at the matching scale Q = M S .
The dotted black line in Fig. 2, which shows very little dependence on tan β, represents the prediction for m h obtained by omitting the one-and two-loop corrections to the quartic Higgs coupling controlled by the bottom Yukawa coupling altogether; the dashed black line includes the one-loop O(g 4 b ) contribution to λ 1 , which, as discussed in Sect. 2.1, we express in terms of the MSSM couplingĝ b ; the solid blue line includes also the two-loop O(g 4 b g 2 3 ) contributions to λ 2 ; finally, the solid red line includes also the two-loop contributions to λ 2 . The comparison between the dashed black line and the solid blue and red lines shows that, when expressed in terms of the MSSM couplingĝ b , the O(g 4 b ) contribution to λ 1 already determines the bulk of the dependence of m h on tan β. Indeed, only at rather large tan β, where the dependence becomes steep, can contributions to λ 2 shift the prediction for m h by more than one GeV. Moreover, those corrections partially cancel out for our choice of MSSM parameters.
Finally, we recall that the strong dependence of m h on tan β depicted in Fig. 2 follows from our choice of signs for the products μ mg and μA t . If both of those products were positive instead of negative, the threshold correction g s b + g Y b in Eq. (12) would suppress the MSSM couplingĝ b -as well as the corresponding contributions to the quartic Higgs coupling and, in turn, to m h -at large values of tan β. If the two products had opposite signs, the dependence of m h on tan β would hinge on whether it is g s b or g Y b that prevails in Eq. (12).

On the effects of dimension-six operators
In MSSM scenarios with SUSY masses up to a couple of TeV, the effects suppressed by powers of v 2 /M 2 S -which are not accounted for when the EFT valid below the SUSY scale involves only renormalizable operators -might still be relevant. In the code SusyHD [67,69] the uncertainty of the prediction for the Higgs mass associated to the omission of those effects is obtained by multiplying the contribution to λ 1 from each SUSY particle by a factor 6 (1 ± 2 v 2 /M 2 i ) , where M i is that particle's mass. In a simplified scenario with tan β = 20, degenerate SUSY masses M i ≡ M S and "maximal" X t = √ 6 M S , the uncertainty arising from missing O(v 2 /M 2 S ) effects was thus estimated in Ref. [67] to be about 0.6 GeV for M S = 1 TeV, and to decrease rapidly for larger M S . The total theoretical uncertainty of the EFT prediction for m h , including also the effects of missing higher-order terms in the matching at the SUSY scale and in the SM part of the calculation, was estimated in Ref. [67] to be less than 2 GeV for M S = 1 TeV, where SusyHD finds m h ≈ 123 GeV. As mentioned in Sect. 1, in that scenario the predictions for the Higgs mass of various fixed-order (or hybrid) codes differ form each other by several GeV, and in general lie outside the estimated uncertainty of the EFT result. In this section we aim to improve the EFT calculation of the Higgs mass at moderate values of M S by including some of the most important O(v 2 /M 2 S ) effects, and to appraise the existing estimate of the uncertainty associated to the missing ones.

Outline of the calculation
In the EFT framework, the effects of O(v 2 /M 2 S ) in the predictions for physical observables such as the Higgs mass arise from non-renormalizable, dimension-six effective operators. The most general dimension-six Lagrangian respecting the field content and symmetries of the SM contains a large number of operators; see Refs. [85][86][87][88][89] for recent reviews. In this section we focus on the two operators that induce oneloop corrections to m 2 h proportional to g 2 t m 4 t /M 2 S and twoloop corrections proportional to g 2 t g 2 3 m 4 t /M 2 S , i.e. the terms suppressed by m 2 t /M 2 S in what are usually denoted as oneloop O(α t ) and two-loop O(α t α s ) corrections to the Higgs mass, where α t ≡ g 2 t /(4π) and α s ≡ g 2 3 /(4π). We write the Lagrangian of the SM extended by dimension-six operators as L EFT = L SM − c 6 |H | 6 + c t |H | 2 t R H T q L + h.c. , (16) where q L and t R are third-generation quarks, is the antisymmetric tensor (with 12 = 1) acting on the SU (2) indices, and, to fix our notation, We stress that the choice of considering only the two dimension-six operators shown in Eq. (16) implies that our treatment of the O(v 2 /M 2 S ) effects is by no means complete, even when we restrict the calculation to the "gaugeless" limit g = g = 0. Indeed, to account for the terms proportional to the resummation of the O(v 2 /M 2 S ) logarithmic corrections to m 2 h beyond two loops, even to account only for the effects controlled by the highest powers of g 3 -i.e., the (n + 1)loop terms proportional to g 2 t g 2n 3 m 4 t /M 2 S ln n (M S /m t ) -we should include in Eq. (16) a set of dimension-six operators involving gluons. 8 However, it must be kept in mind that the suppression by a factor m 2 t /M 2 S implies that, for those corrections to be relevant, the argument of the resummed logarithms cannot be too large. As a result, there is no guarantee that the three-loop (and higher) logarithmic effects of O(v 2 /M 2 S ) that we could account for via resummation are more important than other effects that we are neglecting, such as, e.g., nonlogarithmic three-loop corrections unsuppressed by m 2 t /M 2 S . The sure benefits of extending the SM Lagrangian with the two dimension-six operators of Eq. (16) are that i) we include in our calculation of the Higgs mass the O(v 2 /M 2 S ) part of one-and two-loop corrections that are known to be among the most significant ones, and ii) we can exploit our knowledge of the size of those corrections to estimate the theoretical uncertainty associated to other O(v 2 /M 2 S ) effects that we are neglecting.
The boundary conditions on the Wilson coefficients c 6 and c t are obtained by matching the EFT Lagrangian with the full MSSM Lagrangian at a renormalization scale Q ≈ M S . We start by remarking that those two coefficients receive contributions already at the tree level, controlled by the EW gauge couplings and generated when the heavy Higgs doublet -whose mass we denote by m A -is integrated out of the MSSM Lagrangian: However, in the limit of large tan β both contributions scale like 1/ tan 2 β. For tan β 10 , which we require to saturate the tree-level prediction for m h and allow for stop masses around one TeV, the resulting suppression makes the treelevel contributions to c 6 and c t numerically comparable with the one-loop contributions controlled by the EW gauge couplings, which we are not considering in our study. We will therefore omit the tree-level contributions of Eq. (18) altogether in what follows, and we now move on to summarizing our calculation of the one-and two-loop matching conditions relevant to the O(α t ) and O(α t α s ) corrections to the Higgs mass.
Matching of c t : The one-loop matching condition for c t can be derived by equating the expressions for the pole top-quark 8 Focusing on the CP-even operators, those are f abc G a μν G b νρ G c ρμ , |H | 2 G a μν G a μν and t R σ μν T a H T q L G a μν .
mass computed below and above the matching scale: where 1 t (m t ) is the one-loop self energy of the top quark computed with the external momentum p 2 = m 2 t , and v is the Higgs vev in the EFT, whilev = v 2 1 + v 2 2 is the corresponding quantity in the MSSM. We adopt as usual the DR scheme for the MSSM calculation and the MS scheme for the EFT calculation (note, however, that c 1 t is the same in both schemes). We focus here on the O(g 2 3 ) and O(g 3 t g 2 3 ) contributions to the matching conditions for g t and c t , respectively, which are necessary to reproduce the two-loop O(α t α s ) corrections to the Higgs mass. Defininĝ g t (Q) = g t (Q) (1 + g s t ), and considering that the distinction between v andv does not matter at O(g 2 3 ), we can extract g s t and c 1 t from the terms of O(v) and O(v 3 ), respectively, in an expansion of the stop-gluino contribution to the top self energy in powers of v. Starting from Eq. (B2) of Ref. [22] for the unexpanded self energy, we find where the functions F 6 and F 9 can be found in appendix A of Ref. [66], we defined x Q = m Q 3 /mg and x U = m U 3 /mg , and the term in the last line of Eq. (21) is obtained from the terms in the first five lines by swapping x Q and x U . We note that the right-hand side of Eq.
Matching of c 6 : The matching condition for the Wilson coefficient of the operator |H | 6 in Eq. (16) can, in analogy with the calculation of the matching condition for the quartic Higgs coupling described in Sect. 2.1, be obtained from the derivatives with respect to the Higgs field of the sfermion contributions to the effective potential of the MSSM. In particular, the O(g 6 t ) contribution to the one-loop coefficient c 1 6 at the matching scale Q reads c 1 6 (Q) = 1 36 where V 1 ,t is the stop contribution to the Coleman-Weinberg potential of the MSSM, As outlined in Section 2.3 of Ref. [66], the derivatives of V 1 ,t with respect to the Higgs field can easily be computed after expressing the stop masses m 2 t i as functions of the fielddependent top mass m t =ĝ t |H |, leading to where we used for the derivatives of the one-loop potential shortcuts analogous to those defined in Eq. (5) for the derivatives of the two-loop potential. Explicitly, we find where, following the notation of Eq. (2), we defined x t = m Q 3 /m U 3 and X t = X 2 t /(m Q 3 m U 3 ). Eq. (26) agrees with the corresponding results in Refs. [91,92], 9 which employed the method known as "covariant derivative expansion" [93][94][95] to compute the one-loop matching conditions for all bosonic dimension-six operators induced by integrating the squarks out of the MSSM Lagrangian. In the limit of degenerate squark masses m Q 3 = m U 3 = M S , the O(g 6 t ) contribution to c 1 6 reduces to in agreement with the result presented long ago in Ref. [53].
The O(g 6 t g 2 3 ) contribution to the two-loop coefficient c 2 6 at the matching scale Q reads c 2 6 (Q) = 1 36 where V 2 ,t denotes the contribution to the MSSM scalar potential from two-loop diagrams involving the strong gauge interactions of the stop squarks, given e.g. in Eq. (28) of Ref. [66]. The derivatives of V 2 ,t are again obtained from Eq. (25) after expressing the stop masses and mixing angle as functions of the field-dependent top mass. As in the case of the two-loop matching condition for the quartic Higgs coupling discussed in Sect. 2.1, in the derivatives of the two-loop scalar potential we find terms proportional to ln(m 2 t /Q 2 ), which would diverge in the limit of vanishing top mass. Those terms, however, cancel against analogous terms in δc EFT 6 , which represents the one-loop contribution of the dimension-six operators to the Wilson coefficient of |H | 6 as computed in the EFT. In particular, the contribution relevant at O(g 6 t g 2 3 ) arises from a box diagram with a top-quark loop, three regular Yukawa vertices and one dimension-six vertex. We find where for c 1 t (Q) we use the O(g 3 t g 2 3 ) contribution to the matching condition given in Eq. (21). Finally, the third term on the right-hand side of Eq. (28) arises from the fact that, in analogy with our two-loop calculation of the quartic Higgs coupling, we choose to express the one-loop stop contribution to c 1 6 in terms of the MS-renormalized top Yukawa coupling of the EFT, i.e. g t , as opposed to the DR-renormalized coupling of the MSSM, i.e. theĝ t entering Eqs. (25)- (27). 10 The resulting O(g 6 t g 2 3 ) shift in c 2 6 reads δc shift,t 6 = 6 g s t c 1 6 (Q), 10 On the other hand, in the two-loop corrections the distinction between g t andĝ t amounts to a higher-order effect.
where g s t is given in Eq. (20) and c 1 6 (Q) is given in Eq. (26). The analytic formula for c 2 6 (Q) for generic stop and gluino masses is too lengthy to be printed, and we make it available on request in electronic form. For the simplified choice m Q 3 = m U 3 = mg = M S , we obtain We also remark that c 2 6 (Q) contains terms enhanced by powers of the ratios between the gluino mass and the stop masses. In particular, in the simplified scenario where m Q 3 = m U 3 = M S , X t = ± √ 6 M S and mg M S we find where the sign of the first term within round brackets corresponds to the sign of X t . The presence of power-enhanced terms in the heavy-gluino limit is a well-known consequence of the DR renormalization of the parameters in the stop sector, as discussed in Ref. [22] for the fixed-order calculation of the MSSM Higgs masses and in Ref. [67] for the EFT calculation. Those terms would be removed from the twoloop part of c 6 if we interpreted the soft SUSY-breaking stop masses m Q 3 and m U 3 and the stop mixing X t entering the one-loop part as "on-shell"-renormalized parameters.
Comparison with the fixed-order calculation of m 2 h : We now discuss how the inclusion in the EFT Lagrangian of the dimension-six operators shown in Eq. (16) where m 2 h contains the radiative corrections to the tree-level prediction for the Higgs mass, as computed in the EFT. To avoid the occurrence of large logarithms in these corrections, the couplings λ and c 6 in Eq. (33) should be computed at a renormalization scale Q EW of the order of the masses of the particles running in the loops. Focusing on the one-and twoloop terms that account for the desired O(α t ) and O(α t α s ) corrections, we find where the first two lines are the contribution of one-loop diagrams involving top quarks, the third line is the contribution of two-loop diagrams involving top quarks and gluons computed in the MS scheme (the latter was given, e.g., in Ref. [96]), and we can typically take Q EW ≈ m t . Since for the purpose of this calculation the coefficient c t is first generated at one loop, in the first two lines of Eq. (34) we have exploited the relation m t = g t v + c t v 3 and retained 11 only terms linear in c t (note that in those terms, as well as in those of the third line, the difference between m 2 t and g 2 t v 2 amounts to a higher-order effect). Collecting all the terms in Eq. (33) that involve the coefficients of dimension-six operators, we thus find for the one-and two-loop where the one-loop beta function of c 6 entering the squared brackets in Eq. (36) accounts, at the two-loop level, for the fact that in Eq. (33) the coefficient c 6 should be computed at the low scale Q EW . Isolating the relevant terms in the RGEs for the dimension-six operators given in Refs. [97][98][99][100], we have affects all corrections controlled by g t to the quartic Higgs coupling -namely, the threshold corrections at the SUSY scale and the renormalization-group evolution down to the EW scale -as well as the top-quark contributions to m 2 h given in Eq. (34). This results in an "indirect" contribution of a dimension-six operator to the EFT prediction for m 2 h , which combines with the "direct" contributions controlled by c 6 and c t in Eqs. (33) and (34).
In Fig. 3 we show the deviation induced in the EFT prediction for the Higgs mass by the presence of the dimensionsix operators of Eq. (16). The SM parameters used as input for HSSUSY are the same as those listed at the beginning of Sect. 2.2. We consider a simplified MSSM scenario with tan β = 20 and all soft SUSY-breaking masses of sfermions and EW gauginos, as well as the heavy Higgs-doublet mass m A and the higgsino mass μ, set equal to a common SUSY scale M S ; the trilinear Higgs-stop coupling A t is fixed by the maximal mixing condition A t − μ cot β = √ 6 M S , and A b = A τ = A t ; finally, the gluino mass is set to mg = M S in the left plot and to mg = 2 M S in the right plot. We vary the common SUSY scale between M S = 500 GeV and M S = 1 TeV, and interpret the soft SUSY-breaking stop masses and A t as DR-renormalized parameters at the matching scale Q = M S . We remark that, in the considered range of M S , the prediction of HSSUSY for m h (before the introduction of the dimension-six operators) varies between 120.2 and 123 GeV in the left plot and between 118.7 and 121.9 GeV in the right plot, always several GeV below the value measured at the LHC. Therefore, rather than depicting fully realistic scenarios, the figure is meant to illustrate the relative importance of the different effects induced by dimension-six operators, and how those effects get suppressed by an increase in the SUSY scale.
The dashed blue lines in the plots of Fig. 3 represent the inclusion of the sole operator |H | 6 , with the coefficient c 6 computed at one loop and "frozen" at the matching scale Q = M S . This accounts for the O(m 2 t /M 2 S ) part of the O(α t ) corrections to the Higgs mass, as given in Eq. (35). We see that, in these scenarios, the corresponding shift in m h is negative and rather modest, decreasing from about 470 MeV for M S = 500 GeV to about 90 MeV for M S = 1 TeV.
The solid red lines in the plots of Fig. 3 represent instead the inclusion of both of the operators of Eq. (16), with coefficients c 6 and c t computed at two loops and one loop, respectively, and evolved between the scales M S and Q EW with the RGEs of Eq. (37). This accounts also for the O(m 2 t /M 2 S ) part of the O(α t α s ) corrections to the Higgs mass, as given in Eq. (38). In the scenario shown in the left plot, corresponding to mg = M S , the two-loop O(m 2 t /M 2 S ) corrections appear to be rather small, never reaching even ±100 MeV in the considered range of M S . However, we must take into account that the difference between the dashed blue and solid red lines results from the combination of several effects, namely: (i) the c t -induced shift in the value of g t used in the whole calculation, see Eq. (39), whose "indirect" effects on m h we show for illustration as the dot-dashed green lines; (ii) the inclusion of the two-loop part of the matching condition for c 6 at the SUSY scale; (iii) the evolution of c 6 (and of c t ) between the SUSY scale and the EW scale; (iv) the terms controlled by c t in the radiative corrections to the Higgs mass at the EW scale; see Eq. (34). In the scenario of the left plot, the first three of these effects shift m h by several hundred MeV each for M S = 500 GeV, but they undergo significant cancellations, whereas the fourth effect is considerably less important. On the other hand, in the scenario shown in the right plot, corresponding to mg = 2 M S , the "indirect" effects of the shift in g t are reduced due to a smaller value of c t , and the two-loop contribution to the matching of c 6  To assess the relevance of the O(m 2 t /M 2 S ) logarithmic effects beyond two loops, we removed from the EFT prediction for the Higgs mass the higher-order terms that are picked up by solving numerically the RGEs of the Wilson coefficients in Eq. (37). In practice, we compared our results for m h with those obtained by "freezing" c t at the SUSY scale and truncating the evolution of c 6 to the first order in the perturbative expansion -see the terms within square brackets in Eq. (36). We found that these higher-order logarithmic effects are very small in the considered scenarios: even in the one with mg = M S , characterized by a larger value of c t and hence a stronger scale dependence of both c t and c 6 , the resulting shift in m h reaches a maximum of about 20 MeV for M S ≈ 600 GeV, then decreases for larger M S as the suppression by a factor m 2 t /M 2 S begins to prevail over the logarithmic enhancement.
Finally, the dotted black lines in the plots of Fig. 3 represent a naive estimate of the overall size of the O(v 2 /M 2 S ) corrections to the Higgs mass, corresponding to the "EFT uncertainty" implemented in the code SusyHD. Following Ref. [67], we obtain that estimate by multiplying the contribution to λ 1 from each SUSY particle with mass M i by a factor 13   Effects of dimension-six operators on the EFT prediction for the Higgs mass, as a function of the ratio X t /M S , for a common stop mass scale M S = 1 TeV. In the left plot we take mg = 1 TeV, whereas in the right plot we take mg = 2 TeV. The thin vertical lines in the two plots mark the condition X t /M S = √ 6 . The meaning of the different curves and the values of the remaining MSSM parameters are the same as in Fig. 3, as described in the text the two-loop corrections proportional to g 4 t m 4 t /M 2 S -were as large as the ones that we did compute and had the same sign, the estimate of the "EFT uncertainty" implemented in SusyHD would turn out to be sufficiently conservative in the considered scenarios.
It is legitimate to wonder whether the relatively small size of the O(m 2 t /M 2 S ) corrections found in the scenarios of Fig. 3 is just an accident, perhaps related to the choice X t = √ 6 M S made to ensure a near-maximal prediction for the Higgs mass.
To answer this question, in Fig. 4 we show again the deviation induced in the EFT prediction for the Higgs mass by the presence of the dimension-six operators of Eq. (16), this time as a function of the ratio X t /M S . We set M S = 1 TeV, and take all of the remaining MSSM parameters as in the two scenarios of Fig. 3. In particular, we take mg = 1 TeV in the left plot and mg = 2 TeV in the right plot. The thin vertical lines in the two plots of Fig. 4 mark the condition X t /M S = √ 6, i.e. they map the right edge of the corresponding plots in Fig. 3. The meaning of all other lines is the same as in Fig. 3. Figure 4 shows that, for a given value of M S , the impact of the O(m 2 t /M 2 S ) corrections to the Higgs mass can indeed be larger than the one found when X t /M S = √ 6 . This happens in particular for X t ≈ 0, or for values of |X t /M S | larger than √ 6. The figure also shows that for X t ≈ 0 the SusyHD estimate of the "EFT uncertainty" falls short of the computed O(m 2 t /M 2 S ) corrections to the Higgs mass. Indeed, the main contribution to the SusyHD estimate is the one from stops, which -being proportional to the corresponding contribu-tion to λ 1 , see Eq. (2) -is maximized for |X t /M S | = √ 6 and vanishes for X t = 0 (the small non-zero value of the "EFT uncertainty" visible in the plots at X t = 0 is due to the contributions of EW gauginos and higgsinos). In contrast, Eq. (27) shows that the one-loop stop contribution to c 6 does not vanish for X t = 0, yielding a shift in m h of about 130 MeV. However, we must recall that moving away from the "maximal mixing" condition on X t results in a significant decrease in the EFT prediction for the Higgs mass, e.g. for X t = 0 we would find m h ≈ 110.8 GeV. In order to recover a prediction for m h within a few GeV from the observed value, we would need to raise the SUSY scale M S to several TeV, strongly suppressing all effects of dimension-six operators. Therefore, the SusyHD estimate of the "EFT uncertainty" happens to be at its most conservative precisely in the region of the MSSM parameter space where the O(m 2 t /M 2 S ) effects discussed in this section have a chance to be numerically relevant.

Conclusions
If the MSSM is realized in nature, both the measured value of the Higgs mass and the (so far) negative results of the searches for superparticles at the LHC suggest some degree of separation between the SUSY scale M S and the EW scale. In this scenario the MSSM prediction for the Higgs mass is subject to potentially large logarithmic corrections, which can be resummed to all orders in an EFT approach. Over the past few years this has stimulated a considerable amount of activity, aimed, on one hand, at refining the EFT calculation of the MSSM Higgs mass [65][66][67], and, on the other hand, at combining it with the fixed-order calculations implemented in public codes for the determination of the MSSM mass spectrum [64,73,75,76]. Here we contributed to these efforts by providing a complete determination of the two-loop threshold corrections to the quartic Higgs coupling in the limit of vanishing EW gauge (and first-two-generation Yukawa) couplings, for generic values of all the relevant SUSY-breaking parameters. We also studied a class of one-and two-loop corrections to the Higgs mass suppressed by m 2 t /M 2 S , extending the SM Lagrangian with appropriate dimension-six operators. All of our results are available upon request in electronic form, and they were also implemented in modified versions of the codes SusyHD [69] and HSSUSY [72].
The numerical impact of the various corrections computed in this paper turns out to be small, typically below one GeV in regions of the MSSM parameter space where the prediction for the Higgs mass is within a few GeV from the observed value. We stress that this is in fact a desirable feature of the EFT calculation of the Higgs mass: while the logarithmically enhanced corrections are accounted for by the evolution of the parameters between the matching scale and the EW scale, and high-precision calculations at the EW scale can be borrowed from the SM, the small impact of the two-loop corrections computed at the matching scale suggests that the "SUSY uncertainty" associated to uncomputed higher-order terms should be well under control. In principle, the advantages of an EFT approach are less clear-cut when there is only a moderate separation between the SUSY scale and the EW scale, so that the omission of O(v 2 /M 2 S ) effects in the calculation of the Higgs mass is not warranted. However, our study of the dimension-six operators suggests that the naive estimate of the theoretical uncertainty associated to missing O(v 2 /M 2 S ) effects (or "EFT uncertainty") implemented in the code SusyHD is indeed sufficiently conservative in the relevant regions of the MSSM parameter space. The EFT approach also becomes more complicated when some of the new particles are much lighter than the rest. For example, while our results for the two-loop corrections to the quartic Higgs coupling can be directly applied to the standard split-SUSY scenario by taking the limit of vanishing gluino and higgsino masses, scenarios in which both Higgs doublets are light require a dedicated calculation, in which the effective theory valid below the SUSY scale is a THDM (see, e.g., Ref. [71]).
Finally, we recall that the accuracy of the measurement of the Higgs mass at the LHC has already reached the level of a few hundred MeV -i.e., comparable to the effects of the corrections discussed in this paper -and will improve further when more data become available. If SUSY shows up at last, the mass and the couplings of the SM-like Higgs boson will serve as precision observables to constrain MSSM parameters that might not be directly accessible by experiment, especially in scenarios where some of the superparticle masses are in the multi-TeV range. To this purpose, the accuracy of the theoretical predictions will have to match the experimental one, making a full inclusion of two-loop effects in the Higgs-mass calculation unavoidable. Our results should be regarded as necessary steps in that direction. 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 .