The MSSM without gluinos; an effective field theory for the stop sector

In this article we study the MSSM with stops and Higgs scalars much lighter than gluinos and squarks of the first two generations. In this setup, one should use an effective field theory with partial supersymmetry in which the gluino and heavy squarks are integrated out in order to connect SUSY parameters (given at a high scale) to observables in the stop sector. In the construction of this effective theory, valid below the gluino mass scale, we take into account O(α3)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O(\alpha _3)$$\end{document} and O(Yt,b2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O(Y_{t,b}^2)$$\end{document} effects and calculate the matching as well as the renormalization group evolution. As a result, the running of the parameters for the stop sector is modified with respect to the full MSSM and SUSY relations between parameters are broken. We show that for some couplings sizable numerical differences exist between the effective field theory approach and the naive calculation based on the MSSM running.


Introduction
There are several theoretical arguments for a light stop in supersymmetric theories. Foremost, in natural supersymmetry (SUSY) light stops are required to cancel the quadratic divergence of the Higgs mass originating from the selfenergy involving a top quark, while the other supersymmetric partners can be much heavier [1,2] due to the smaller couplings to the Higgs. Moreover, the renormalization group equations (RGE) of the minimal supersymmetric standard model (MSSM) generically drive the bilinear mass term parameters of the third generation squarks to lower values (compared to the first two generations) due to their nonnegligible Yukawa couplings [3][4][5][6][7][8]. Although the measured Higgs mass of around 125 GeV [9,10] prefers rather heavy (around the TeV scale) [11][12][13] rather than light stops in the MSSM, this is not necessarily the case in the NMSSM [14], in λSUSY models [15], models with light sneutrinos [16] or in supersymmertic models with additional D-term [17] or F-term [18] contributions to the scalar potential. Also large (or even maximal [19][20][21]) stop mixing angles help to get the right Higgs mass with rather light stops.
LHC searches for top squarks (using simplified models) set a lower bound on its mass of around mt 1 = 300 GeV, which however heavily depends on the neutralino mass. Depending on the stop and the neutralino mass, different decay modes are studied. For the decay channelt 1 → tχ 0 1 [22][23][24], the limits are quite stringent, even though for light neutralinos very light stops cannot be excluded due to the high tt-background [25]. The three-body decayt 1 → W bχ 0 1 was analyzed theoretically in [26] and experimentally in [27]. Finally the decayt 1 → cχ 0 1 and the less important four-body decayt 1 →χ 0 1 d i ff are treated in [28][29][30] and constraints were derived by the ATLAS collaboration from the monojet analysis in [31]. Some bounds can be avoided in kinematic boundary regions or once non-minimal flavor violation is included. However, recently efforts of closing these gaps have been made [32][33][34][35] and stops should in general not be lighter than 300 GeV. Nevertheless, the mass bound for the stop is still weaker than the strong bounds on the squark masses of the first two generations and also on the gluino mass [36,37]. For sbottom quarks LHC searches suggest masses of above 800 GeV [38,39]. The bounds on sparticles with EW interactions only are much less stringent [40][41][42][43][44]. For example, in the case of heavy winos the Higgsino mass parameter μ has only to be larger than 350 GeV [45]. It can be shown, however, that by changing the assumptions on the composition of charginos and neutralinos, collider limits can get even further weakened [46][47][48][49]. For the Higgs bosons, dif-ferent fits [50][51][52][53][54][55][56][57][58] suggest an alignment limit, in which the lightest CP-even Higgs boson takes the role of the SM Higgs. Collider limits on non-SM Higgs bosons for large values of tan β suggest that CP-odd Higgs bosons should be heavier than 800 GeV [59,60].
If the gluino (or the squarks of the first two generations [2,61]) is much heavier than the stops, an effective theory (EFT) with partial SUSY must be constructed in which the gluino (squarks) is integrated out [62][63][64]. Such a hierarchy can for example be achieved for MSSM-like models in a Scherk-Schwarz breaking scenario [65][66][67][68]. The construction of the effective theory for the stop sector is the goal of this article. Assuming a common large mass of order M for the gluino and the squarks of the first two generations, we compute the matching condition between the full MSSM and the effective theory, including one-loop contributions which are enhanced by powers of M. Furthermore, since some supermultiplets are partially integrated out in the effective theory, the supersymmetric relations between gauge/Yukawa couplings, gaugino/Higgsino couplings and four-scalar couplings are broken in the effective theory by radiative corrections. Therefore, these couplings in the effective theory have an independent renormalization group evolution, as discussed in [62,[69][70][71][72][73][74][75][76] mainly for the gauginomatter couplings.
This article is structured as follows: In the next section we establish our effective theory for the stop sector and calculate the matching as well as the running of the relevant parameters at order α 3 and Higgs self-coupling effects). This section is followed by a numerical analysis in Sect. 3. Finally we conclude in Sect. 4.

The effective theory for the stop sector
The aim of this section is to construct the effective theory for the MSSM stop sector, including O(α 3 , Y 2 t,b ) and enhanced effects. As noted before, we assume that the gluino and the squarks of the first two generations are much heavier, with masses of the order M, than the stops, the Higgs scalars and the Higgsinos. The left-handed sbottom is also assumed to be light such that it remains in the effective theory, forming an SU (2) multiplet with the left-handed stop. However, we assume that the right-handed sbottom is heavy, with the mass of the order M. Therefore, we consider the following effective Lagrangian which is valid below the scale M:  (2) indices in the two-squark-two-Higgs interactions are contracted within the parentheses. i, j are the SU (2) indices and the dot denotes the contraction of SU (2) indices as For simplicity, we also assume that the electroweak gauginos and sleptons are heavy. However, since we neglect O(g 2 1 ), O(g 2 2 ) effects in the following, relaxing this assumption would leave our RGEs unchanged. We also ignore the non-holomorphic Higgs-quark couplings which are induced at the loop level [77][78][79][80][81][82][83][84].

Tree-level matching
At the matching scale M the Lagrangian of Eq. (1) has to be compared to the one of the full MSSM (see for example [85][86][87][88]) which originates from the superpotential the soft SUSY breaking terms and the D terms where T A are the generators of SU (3) in the fundamental representation. The matching conditions for the bilinear terms and the trilinear couplings arē The couplings between squarks and Higgs bosons are generated by F-and D-terms in the MSSM Lagrangian. At the scale M, they are given by keeping only Yukawa couplings and g 3 .

One-loop matching
For the matching, we need to include the one-loop effects enhanced by powers of M since their contributions may be comparable to the tree level ones shown in the previous subsection. They can only appear in bilinear and trilinear terms, as seen by dimensional analysis. The bilinear terms receive the following shifts at the matching scale μ R = M: For the trilinear term the shift reads All the other parameters relevant for the stop sector are dimensionless and therefore do not receive any M enhanced corrections.

Renormalization group evolution
The running of the full MSSM parameters [3][4][5][6][7] is known at the two-loop level [8,[89][90][91][92]. Here we give the one-loop beta functions to O(α 3 , Y 2 t,b ) for the parameters of our effective theory in Eq. (1). The corresponding results for the full MSSM are summarized in the appendix. For the strong coupling constant we have (t ≡ log μ R , where μ R denotes the renormalization scale) where the first term on the right hand side is the SM contribution. The effective quark-quark-Higgs Yukawa couplings evolve according to while the evolution of the ones entering the Higgsino-quarksquark vertex is determined by For the Higgs mass parameters we find and for the bilinear squark mass terms The Higgsino mass in the effective theory evolves as and the effective trilinear Hqq coupling as Finally for the quartic H Hqq andqqqq couplings one obtains Note that in all equations above we assumed real parameters. However, all formula can easily be generalized to the complex case by simply replacing a square by the absolute value squared. By integrating these RGEs from M to the stop mass scale mt , we obtain the O(α 3 , Y t,b ) contributions enhanced by log(M/mt ).

Stop masses
In the effective theory, the stop mass matrix in the (t L ,t R ) basis reads by SU (2) gauge symmetry.

Numerical analysis
From the previous analysis, we can see that, by integrating out the gluino and the squarks of the first two generations, parameters which were originally related via SUSY • The top Yukawa coupling Y t In the full MSSM, the Yukawa coupling Y t of the superpotential enters top-top-Higgs, stop-stop-Higgs couplings as well as stop-squark-Higgsino couplings in the same way. However, in the EFT these couplings are independent quantities and they evolve differently below the scale M. This is depicted in Fig. 1, where the evolution of Y t in the naive approach using MSSM RGE is compared to those ofȲ t ,Ȳ q 3L ,Ȳ t R andȲt ≡μ t /μ in the EFT. When the values ofȲ t and Y t are determined at the stop mass scale to give the SM running top mass, their values at the scale M are quite different. Note that these couplings are dimensionless and therefore do not depend on the choice of the parameters formt ,mQ andĀ t . • The quartic coupling of right-handed stops λ 8 In the full MSSM the quartic coupling of right-handed stops λ 8 is given by 1 6 g 2 3 by SUSY relation and evidently also evolves in the same way as 1 6 g 2 3 . However, in the EFT λ 8 andḡ 2 3 follow different RGEs below the scale M, as seen in Fig. 2. The relative difference at the scale mt amounts to roughly 30%. Again, since λ 8 has no mass dimension, its running does not depend onmt ,mQ and A t .
Among the quartic scalar couplings λ 1−8 , the running of λ 8 in the EFT exhibits the largest deviation from the one in the full theory. This is due to symmetry factors, leading to large coefficients of the box diagrams and self-couplings which are responsible for a change in sign on the g 4 3 -dependence. The deviations of the other couplings λ 1−7 from the ones in the full theory are either positive or negative, but are smaller than 20% for our parameter set. We therefore do not show the figures of their runnings here.

Conclusions
In this article, we constructed an effective theory of the stop sector obtained from the full MSSM by integrating out the first and second generation of squarks and the gluino (which we assume to have a common mass of the order M). We computed the matching effects for the dimensionful quantities which are enhanced by powers of M at O(α 3 , Y 2 t,b ). In addition, we obtained the complete O(α 3 , Y 2 t,b ) RGEs of the couplings within the EFT. In the numerical analysis we highlighted that couplings which are related via SUSY identities within the full MSSM have different RGEs within the EFT, which can lead to sizable differences. We illustrated this effect for the top Yukawa couplings and the quartic coupling of right-handed stops, finding differences up to 30% between the EFT and the naive approach. Such deviations could play a role in a future test of the stop-stop or stop-Higgs interactions which also enter the calculation of the Higgs mass.