#### Abstract

We present the impact of two-loop corrections on the mass of the lightest Higgs boson in the minimal *R*-symmetric supersymmetric standard model (MRSSM). These shift the Higgs boson mass up by typically 5 GeV or more. The dominant corrections arise from strong interactions, and from the gluon and its superpartners, the sgluon and Dirac gluino, and these corrections further increase with large Dirac gluino mass. The two-loop contributions governed purely by Yukawa couplings and the MRSSM , parameters are smaller. We also update our earlier analysis which showed that the MRSSM can accommodate the measured Higgs and *W* boson masses. Including the two-loop corrections increases the parameter space where the theory prediction agrees with the measurement.

#### 1. Introduction

The recent discovery at the LHC of a particle consistent with the long sought Higgs boson seemingly completes the Standard Model (SM). The mass of the particle is measured with an astonishingly high accuracy of GeV [1]. The precise determination of this mass is of paramount importance not only within the context of the Standard Model, but also for finding the path beyond it. In fact, a number of experimental observations suggest that the SM cannot be the ultimate theory and many theoretical scenarios for the beyond SM (BSM) physics have been proposed in past decades. In some models of BSM, in particular in supersymmetric extensions of the SM, the Higgs boson mass can be predicted. However, the current experimental accuracy is far better than theoretical predictions for Higgs boson mass in any given model of BSM physics. From the point of view of theory, the best accuracy has been achieved in the minimal supersymmetric extension of the SM (MSSM), in which the discovery of the Higgs boson and the determination of its mass have given a new impetus to the theoretical efforts. The most recent improvements comprise the inclusion of leading three-loop corrections [2, 3], resummation of leading logarithms beyond the two-loop level [4, 5], inclusion of the external momenta of two-loop self-energies [6, 7], and the evaluation of the -contributions in the complex MSSM [8, 9]. The MSSM two-loop corrections controlled by Yukawa couplings and have been known for quite some time for the real MSSM (see the above references for an overview of the literature).

The absence of any direct signal of supersymmetric particle production at the LHC and the observed Higgs boson mass of ~125 GeV being rather close to the upper value of 135 GeV achievable in the MSSM are a strong motivation to consider nonminimal SUSY scenarios. In fact, nonminimal SUSY models can lift the Higgs boson mass (at the tree-level by new - or -term contributions or at the loop level from additional new states), which makes these models more natural by reducing fine-tuning. They can also weaken SUSY limits either by predicting compressed spectra or by reducing the expected missing transverse energy or by reducing production cross sections. The comparison of the measured Higgs boson mass with the theoretically predicted values in any given model is therefore highly desirable. Although the theoretical calculations for the SM-like Higgs boson mass in such models are less advanced, progress is being made in the development of highly automated tools which greatly facilitate the computations in nonminimal SUSY models: SARAH [10–12] automatically generates spectrum generators similar to SPheno [13, 14]; FlexibleSUSY [15] automatically generates spectrum generators similar to Softsusy [16].

In a recent paper [17] we considered the MRSSM, a highly motivated supersymmetric model with continuous -symmetry [18, 19] distinct from the MSSM. Since -symmetry forbids soft Majorana gaugino masses as well as the higgsino mass term, additional superfields are needed. The MRSSM has been constructed in [20] as a minimal viable model of this type. It contains adjoint chiral superfields with -charge 0 for each gauge sector and two additional Higgs weak iso-doublet superfields with -charge 2. It has been also argued that -symmetry generically forbids large contributions to CP- and flavor-violating observables due to the absence of chirality-changing Dirac gluino couplings [20, 21], relaxing flavor constraints on the sfermion sector, although recently it has been shown that the dramatic chirality-flip suppression of [20] can only work in a limited number of scenarios, and in general a certain correlation between flavor structure of fermion masses and superpartner spectrum is required [22]. Also, Dirac gluinos suppress the production cross section for squarks, making squarks below the TeV scale generically compatible with LHC data. Furthermore, models with -symmetry and/or Dirac gauginos contain promising dark matter candidates [23–25], and the collider physics of the extra, non-MSSM-like states has been studied [26–34].

In [17] the complete next-to-leading order computation and discussion of the lightest Higgs boson and boson masses have been performed. We showed that the model can accommodate measured values of these observables for interesting regions of parameter space with stop masses of order 1 TeV (a similar analysis has been done in [35], where also a welcome reduction of the level of fine-tuning was found). The outcome of the paper was not obvious since in the MRSSM (i) the lightest Higgs boson tree-level mass is typically reduced compared to the MSSM due to mixing with additional scalars, (ii) the stop mixing is absent, and (iii) -symmetry necessarily introduces an scalar triplet, which can increase already at the tree-level. Nevertheless, we identified benchmark points BMP1, BMP2, and BMP3 illustrating different viable parameter regions for , respectively, and also verified that they are not excluded by further experimental constraints from Higgs observables, collider, and low-energy physics.

These promising results motivate a more precise computation of the Higgs boson mass in the MRSSM and a more precise parameter analysis. Technically, this is facilitated by the Mathematica package SARAH, recently updated by providing SPheno routines, which calculate two-loop corrections to the CP-even Higgs scalars masses in the effective potential approximation and the gaugeless limit [36]. This is the level of precision of the established MSSM predictions except for the refinements mentioned above. It is also the level of precision at which the proof [37] applies that the employed regularization by dimensional reduction preserves supersymmetry. First applications of the improved SARAH version to the calculations of the Higgs boson masses in the -parity violating MSSM [38] and next-to-minimal SSM [39] have been published.

Since a judicious choice of the model parameters was needed to meet experimental constraints and an estimate of unknown two-loop contributions was presented, it is of immediate interest to verify our findings at higher precision with the new SARAH version. The aim of the current paper is to calculate two-loop corrections for the Higgs boson mass in the same MRSSM setup as in [17] and present an update of the results obtained there.

The paper is organized as follows. After a short recapitulation of the MRSSM setup in Section 2, we explain in Section 3 our calculation framework and discuss the dependence of two-loop corrections on parameters that entered already at the one-loop level. The dependence on parameters that enter only at the two-loop level is investigated in Section 4. In Section 5 we provide an update to the analysis presented in [17] using the two-loop corrected masses of Higgses, before concluding in Section 6.

#### 2. The MRSSM

The MRSSM has been constructed in [20] as a minimal supersymmetric model with unbroken continuous -symmetry. The superpotential of the model reads as where are the MSSM-like Higgs weak iso-doublets and are the singlet, weak iso-triplet, and -Higgs weak iso-doublets, respectively. The usual MSSM -term is forbidden; instead -terms involving -Higgs fields are allowed. -terms are similar to the usual Yukawa terms, where -Higgs and or play the role of the quark/lepton doublets and singlets.

The usual soft mass terms of the MSSM scalar fields are allowed just like in the MSSM. In contrast, -terms and soft Majorana gaugino masses are forbidden by -symmetry. The fermionic components of the chiral adjoint, for each standard model gauge group , , , respectively, are paired with standard gauginos to build Dirac fermions and the corresponding mass terms. The Dirac gaugino masses generated by -type spurions produce additional terms with the auxiliary -fields in the Lagrangian, which after being eliminated through their equations of motion lead to the appearance of Dirac masses in the scalar sector as well. For our phenomenological studies of two-loop effects we take the soft-breaking scalar mass terms that have been considered in [17]where the holomorphic mass terms for adjoint scalars, which might lead to tachyonic states, have been neglected (see also [40, 41] for discussions that these terms can be subdominant within a broad definition of gauge mediation).

The electroweak symmetry breaking (EWSB) is triggered by nonzero vacuum expectation values of neutral EW scalars, which are parameterized as -Higgs bosons carry -charge 2 and therefore do not develop vacuum expectation values. We stress that in general the mixing of with and leads to a reduction of the lightest Higgs boson mass at the tree-level compared to the MSSM.

#### 3. Higgs Mass Dependence on Superpotential Parameters

We now present the MRSSM Higgs boson mass prediction at the two-loop level. We use the same renormalization scheme as in [17], where all SUSY parameters are defined in the scheme and , , , and are determined by minimizing the effective potential at the two-loop order. The discussion is divided into two parts. In the present section we begin with the one-loop contributions, which are dominated by terms of , where collectively denotes squares of the superpotential couplings and . We then discuss the two-loop contributions of , that is, ones which depend on parameters which already play a role at the one-loop level. In the subsequent section we then discuss those two-loop corrections which involve new parameters.

In the usual MSSM, the one-loop contributions to the Higgs boson mass are dominated by top/stop contributions. In the MRSSM, these contributions are also important, but they are simpler since stop mixing is forbidden by -symmetry (corresponding to the MSSM parameter ). This implies that the top/stop contributions cannot reach values as high as in the MSSM for a given stop mass scale. However, as mentioned above, the MRSSM superpotential contains new terms governed by and which have a Yukawa-like structure. References [17, 35] have given a useful analytical approximation for these contributions. In the limit , , , and large , we get This result shows a behavior proportional to , , and . This is similar to the top/stop contributions as ’s and appear in a similar fashion in superpotential.

We expect therefore that the two-loop result will depend on these model parameters (which already entered at the one-loop level) in a manner similar to the pure top quark/squarks two-loop contributions, that is, similar to the MSSM contributions without stop mixing.

In Figures 1 and 2 the dependence of the lightest Higgs boson mass calculated at tree-, one-, and two-loop levels for two benchmarks BMP1 and BMP3 on different model parameters is shown. All parameters except the ones shown on the horizontal axes are set to the values of the benchmark points defined in [17] (see Table 2). Indeed behavior of the two-loop corrections is very similar to the one of the corresponding one-loop corrections. The numerical impact of the two-loop -contributions is rather small, typically less than GeV, except for very large , where they can reach several GeV. Particularly, the strong dependence for large is already manifest for the tree-level mass; this is due to the mixing with the singlet state already present in the tree-level mass matrix.

One should remember that very large one-loop contributions are required to bring the predicted Higgs boson mass close to the experimental one. In the preferred parameter regions, , are large but still moderate enough not to blow up the two-loop contributions. Although the large values of couplings will lead to Landau poles at scales below the GUT scale; this is not problem for our phenomenological approach, and even in top-down approaches as in [40] perturbativity up to the GUT scale is not required.

Overall, the total two-loop contributions (including the ones to be discussed in the subsequent section) are in the range between and GeV, except in the very large regions. This is in agreement with the estimate given in [17], and it confirms the validity of the perturbative expansion in spite of the large one-loop corrections.

#### 4. QCD Corrections and the Two-Loop Corrected Higgs Boson Mass

At two-loop level the strongly interacting sector and the strong coupling appear directly in the Higgs boson mass predictions. These two-loop corrections involve not only the gluon but also the Dirac gluino and the sgluon, the scalar component of the octet superfield . They can be expected to be sizable, and they depend on the gluino Dirac mass and sgluon soft mass parameters. These parameters already play a role at lower order, appearing in corrections to (through threshold corrections to ), though the influence on, for example, top mass is negligible. The gluino Dirac mass parameter appears not only directly as the gluino mass but via (2), also in couplings and mass terms of sgluons, inducing the mass splitting. In [17] a simplifying assumption was made that masses of the scalar and pseudoscalar components of (complex) sgluon field were equal, since it was unimportant for that analysis between the real and imaginary parts of the sgluon field, . The masses of the scalar sgluons and pseudoscalar sgluons are related by the tree-level formula , where is equal to the soft-breaking parameter [26, 28]. The relevant vertices and Feynman rules are depicted in Figure 3. We assume real , so only the scalar acquires the direct coupling to sfermions proportional to , via (2).

The structure of the strong corrections is thus markedly different from the MSSM case, where only the Majorana gluino and the gluon appear. In the following, we study the magnitude and the behavior of the corrections as a function of the parameters and .

##### 4.1. Analytic Formulas

As in the previous section, we begin with an analytic approximation for the leading contributions of , that is, two-loop strong corrections proportional to . This provides us with qualitative insight and serves as a check of the code. Generally, in the gaugeless limit (in which the two electroweak gauge couplings are neglected), the two-loop corrections from gluinos and sgluons contribute only to the diagonal part of submatrix of the scalar Higgs boson mass matrix. In the MRSSM terms contribute only to element. This already constitutes a difference to the MSSM, where -term violates -symmetry and Peccei-Quinn symmetry leading to couplings of stops to .

Figure 4 shows two-loop diagrams contributing to the Higgs boson mass at that explicitly depend on and/or . These diagrams provide the following contribution to the effective potential: where the functions and are defined in [42]. The effective potential depends on through stop masses, which in the gaugeless limit approach Equation (6) can be obtained from [42] by applying translation rules from real fields to complex ones. Many such rules can be found in [36]; an additional rule needed here for the case of a Lagrangian , where , is .

An important difference to the MSSM is that contributions with fermion mass insertions, corresponding to -type contributions in [42], are not present in the MRSSM. Such contributions vanish due to the lack of - mixing between squarks. Hence the gluino mass appears in a simpler way than in the MSSM. Likewise, the sgluon only enters via the -type diagram of Figure 4. An -type diagram vanishes due to the color structure.

The corresponding two-loop contribution to Higgs boson mass matrix element in zero-momentum approximation is then given by (as pointed out in [36], in SARAH and SPheno the two-loop tadpole contributions are included directly in vacuum minimization condition and not in (8))For large , corrections of order cannot be neglected any more. But since they contribute only to matrix element, their impact on mass of the lightest Higgs, which stems mainly from element, is small. Results of (6) were compared with the results of two-loop routines from the SARAH-generated SPheno module.

##### 4.2. Numerical Analysis

We now turn to the numerical analysis of the complete two-loop corrections to the SM-like Higgs boson mass, using the full evaluation within the framework of SARAH and SPheno. Figures 5(a) and 5(b) focus on the gluino and sgluon mass dependence, which arises mainly from the corrections; they show the two-loop corrections as a function of the gluino mass parameter for two different values of the soft sgluon mass, and 10 TeV for two benchmarks BMP1 and BMP3; other parameters are fixed at benchmark values. For comparison, the two-loop result without the sgluon contribution is shown as well (i.e., without the first diagram of Figure 4). We also plot the MSSM prediction with strong stop mixing and without any sfermion mixing at tree-level.

**(a)**

**(b)**

**(c)**

Figures 5(a) and 5(b) show that the dependence in the MRSSM without sgluon contributions is very similar to the one in the MSSM without stop mixing. The corresponding thin solid red and thin dashed light blue curves in Figure 5 show a characteristic drop for large gluino masses. This is understandable as in the MSSM without sfermion mixing the gluino contribution is precisely the same as in the MRSSM and given by the two corresponding diagrams in Figure 4. The Dirac or Majorana nature of the gluino does not matter since the Dirac partner, the octet superfield , has no direct couplings to quark superfields. A few TeV gluino masses slightly increase the Higgs boson mass, but for larger values of the function becomes negative and drives the correction downwards.

In the full MRSSM calculations, including the sgluon diagrams strongly changes the behavior. Surprisingly, the full MRSSM two-loop contributions resemble the MSSM contributions with large stop mixing. In both cases, large gluino masses strongly enhance the Higgs boson mass, however, for different reasons. In the MSSM the increase can be traced back to the additional -type diagram which is directly proportional to and which vanishes in the limit of no stop mixing. In the MRSSM, on the other hand, the sgluon diagram grows with due to both the sgluon-stop-stop coupling, which scales like , and an increase in the scalar (but not pseudoscalar) sgluon mass. Due to the sgluon contributions the total two-loop contributions to the Higgs boson mass in the MRSSM are larger than the ones in the MSSM. They are further increased by heavy sgluons.

Figure 5(c) compares the numerical impact of individual contributions by successively switching off contributions. It allows us to read off the contributions from sgluon, gluino, and gluon, of , and the remaining two-loop contributions (particularly contributions). The gluon diagrams alone contribute approximately GeV. The negative gluino and the positive sgluon corrections together amount to an additional upward shift of the Higgs boson mass, which can reach several GeV for large Dirac gluino masses. The remaining contributions are far smaller and amount to around GeV for contributions and GeV for the remaining contributions.

#### 5. Update of Benchmarks

In this section we present an update of the analysis of [17], using the more precise evaluation of the Higgs boson mass. Reference [17] studied the mass predictions of and lightest Higgs bosons in the MRSSM and showed that agreement with experimental data is possible, in spite of tree-level shifts from violations of custodial symmetry and from mixing with other Higgs states, respectively.

Table 2 shows benchmark parameter points defined in that reference. They exemplify parameter regions in which and agree with experiment. They are characterized by large , rather light Dirac higgsinos and gauginos, and they have , respectively.

For all three benchmark points the two-loop correction to is around GeV. As discussed in the previous sections, the largest part of this is due to corrections. The MRSSM-specific corrections of are small since the values of , though large, are still not as large as needed to make these corrections dominate; see Figures 1 and 2 for two out of three benchmarks. The magnitude of the total two-loop correction is consistent with the theory error estimate given in [17].

The upward shift of implies that it is easier to obtain agreement with the measured value; that is, smaller values of are sufficient. In Table 3 we provide new, slightly modified benchmark points, whose definitions differ only in the values of . The two-loop Higgs boson mass prediction agrees well with experiment, and the good agreement of with experiment is unchanged. Likewise, both the old and the new set of benchmark points pass checks against HiggsBounds [43–45] and HiggsSignals [46, 47].

In Figure 6 we give an update to some of the subfigures from Figures 4 and 5 of [17]. These show the predictions of and as contour lines in several two-dimensional parameter spaces. The Higgs boson mass is evaluated at the two-loop level. As discussed before, with the exception of the regions of very large , there is a general positive contribution to the lightest Higgs boson mass between 4 and 5 GeV. Accordingly, the contour lines, in particular the central green region in which the Higgs boson mass agrees with experiment, shift to slightly lower values of . Also, the overlap region, where Higgs and boson masses agree with experiment, is enlarged.

**(a)**

**(b)**

#### 6. Conclusions

In this work we have presented the impact of two-loop corrections on the mass of the lightest Higgs boson in the MRSSM. The calculation has been performed using the framework of SARAH in the approximation of the vanishing electroweak gauge couplings and external momenta of the Higgs self-energies. The code has been cross-checked with an analytic calculation of the most important new corrections. We have separately analyzed the impact of contributions involving -couplings, which already appear in the one-loop corrections, and of the strong corrections involving gluon, Dirac gluino, and sgluon exchange.

In the previous work [17] and the present paper we have found that the lightest Higgs boson mass in the MRSSM differs from the one in the usual MSSM in several respects. At tree-level the additional mixing with additional scalar states reduces the MRSSM Higgs mass below the MSSM value. At the one-loop level, the top/stop contributions cannot be as large as in the MSSM, because stop mixing is forbidden by -symmetry. However, the new contributions from the superpotential -terms have a similar structure as the top/stop contributions. If -couplings are similar in magnitude to the top Yukawa coupling, the lightest Higgs boson mass can easily be in the ballpark of the experimentally allowed range.

The two-loop corrections governed by these -couplings, however, amount to only 1 GeV or less in parameter regions in which the Higgs boson mass agrees with experiment. The most important two-loop contributions are the strong corrections of . As we have shown the Dirac gluino and gluon contributions alone are very similar to the MSSM strong contributions for vanishing stop mixing. The inclusion of the sgluons changes the picture. The sgluon contributions are positive and rise with the Dirac gluino mass, such that the total corrections of the MRSSM are larger than the ones of the MSSM, independently of the magnitude of stop mixing.

Overall, the MRSSM two-loop corrections to the lightest Higgs boson mass are typically positive. For example, for the benchmark parameter points proposed in [17], the two-loop corrections to the Higgs boson mass amount to approximately GeV, within the error estimate of that reference. Since perturbation theory shows a converging behavior and since -corrections are subdominant (for less than around 1.2), we estimate the remaining theory uncertainty to be not larger than the one of the MSSM.

The positive two-loop corrections make it easier to achieve agreement between the theory prediction for the lightest Higgs boson mass and the measured value. We have provided an update of the analysis of [17], showing parameter regions of simultaneous agreement of the Higgs and boson mass predictions with experiment. Compared to [17], the allowed parameter regions are slightly larger and located at smaller values of -couplings.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank Kilian Nickel and Florian Staub for communication about SARAH. This work is supported in part by the Polish National Science Centre Grants under OPUS-2012/05/B/ST2/03306, DEC-2012/05/B/ST2/02597, the European Commission through the Contract PITN-GA-2012-316704 (HIGGSTOOLS), the German DAAD PPP Poland Project 56269947 “Dark Matter at Colliders,” and DFG Research Training Group 1504 and the DFG Grant STO 876/4-1.