Precise prediction for the Higgs-Boson Masses in the μ ν SSM with three right-handed neutrino superfields

The μνSSM is a simple supersymmetric extension of the Standard Model (SM) capable of describing neutrino physics in agreement with experiments. We perform the complete one-loop renormalization of the neutral scalar sector of the μνSSM with three generation of right-handed neutrinos in a mixed on-shell/DR scheme. We calculate the full one-loop corrections to the neutral scalar masses of the μνSSM. The one-loop contributions are supplemented by available MSSM higher-order corrections. We obtain numerical results for a SM-like Higgs-boson mass consistent with experimental bounds, while simultaneously agreeing with neutrino oscillation data. We illustrate the distinct phenomenology of the μνSSM in scenarios in which one or more right-handed sneutrinos are lighter than the SM-like Higgs boson, which might be substantially mixed with them. These scenarios are experimentally accessible, on the one hand, through direct searches of the right-handed sneutrinos decaying into SM particles, and on the other hand, via the measurements of the SM-like Higgs-boson mass and its couplings. In this way the parameter space of the μνSSM can be probed without the need to propose model dependent searches at colliders. Finally, we demonstrate how the μνSSM can simultaneously accommodate two excesses measured at LEP and LHC at ∼ 96 GeV at the 1σ level, while at the same time reproducing neutrino masses and mixings in agreement with neutrino oscillation measurements. ∗email: thomas.biekotter@csic.es †email: Sven.Heinemeyer@cern.ch ‡email: c.munoz@uam.es ar X iv :1 90 6. 06 17 3v 2 [ he pph ] 1 6 A ug 2 01 9


Introduction
The scalar particle at ∼ 125 GeV discovered by the ATLAS [1] and CMS [2] experiments has so far shown to be consistent with the Standard Model (SM) Higgs-boson prediction. The Higgs boson was the last missing piece in the description of electroweak symmetry breaking (EWSB) and the generation of masses of fundamental particles within the SM. The measurement of the mass of this new state already reached a remarkable precision [3]: 1 m H = 125.09 ± 0.21(stat.) ± 0.11(syst.) . (1.1) µνSSM solves the µ-and the ν-problem (neutrino masses) simultaneously without the need to introduce additional energy scales beyond the SUSY-breaking scale. In contrast to the (N)MSSM, R-parity and lepton number are not conserved, leading to a completely different phenomenology, characterized by distinct prompt or displaced decays of the LSP [29][30][31][32]. Although the LSP is not stable anymore, the µνSSM can provide a dark matter candidate with a gravitino with a lifetime longer than the age of the universe [33][34][35][36]. The breaking of R-parity is induced by a neutrino Yukawa term, with the size of the couplings Y ν ij ≤ 10 −6 determined by the electroweak seesaw. Because of the values of Y ν ij , mixings between SM particles and their supersymmetric partners are suppressed. Nevertheless, the additional sources of mixing effects induce a conceptually modified spectrum compared to the MSSM and the NMSSM. The spectrum will be described in detail in Sect. 2.
SUSY relates the quartic couplings of the neutral scalar potential to the gauge couplings of the underlying field theory. Therefore, within SUSY the scalar masses can be predicted in terms of other model parameters, and the precise value of the SM-like Higgs-boson mass is of particular significance. However, the SM-like Higgs-boson mass predictions strongly depend on quantum corrections which can be calculated only to certain order in perturbation theory. Missing higherorder contributions lead to a sizable amount of uncertainty which is usually of a few GeV (see below for details), hence an order of magnitude larger than the experimental uncertainty. This is why a lot of effort is made to predict the Higgs mass to the highest possible precision [37]. We briefly summarize the status of Higgs-mass predictions in the MSSM, the NMSSM, and the µνSSM in the following.
In the MSSM the tree-level mass can be predicted by just two SUSY parameters, i.e., the ratio of the vevs of the Higgs doublets tan β, and either the mass of the CP-odd Higgs boson M A or the mass of the charged Higgs boson M H ± , leading to an upper bound given by the Zboson mass. Large loop corrections are needed to achieve a Higgs-boson mass of ∼ 125 GeV. Beyond the one-loop level, the dominant two-loop corrections of O(α t α s ) [38][39][40][41][42][43], O(α 2 t ) [44,45], O(α b α s ) [46,47] and O(α t α b ) [46] are known (here we use α f = (Y f ) 2 /(4π), with Y f denoting the fermion Yukawa coupling). These corrections, together with a resummation of leading and subleading logarithms from the top/stop sector [48] (see also [49,50] for more details on this type of approach), a resummation of leading contributions from the bottom/sbottom sector [46,47,[51][52][53][54] (see also [55,56]) and momentum-dependent two-loop contributions [57,58] (see also [59]) are included in the public code FeynHiggs [40,48,[60][61][62][63][64][65][66][67]. The most recent version of FeynHiggs contains an improved effective field theory calculation relevant for large SUSY scales [64,66,68]. The complete two-loop QCD contributions in the CP-violating MSSM were calculated in Ref. [69], but not yet included in FeynHiggs. A (nearly) full two-loop effective potential (EP) calculation, including even the leading three-loop corrections, has also been published [70,71], which is, however, not publicly available as a computer code. Furthermore, another leading three-loop calculation of O(α t α 2 s ), depending on the various SUSY mass hierarchies, has been performed [72,73], resulting in the code H3m and is now available as a stand-alone code [74]. It was proven that regularization by dimensional reduction preserves supersymmetry at the required three-loop order [75]. A new calculation of the three-loop contributions of the O(α t α 2 s ) extends the validity of these corrections to the whole parameter space of the CP-conserving MSSM [76]. Most recently, the leading logarithmic terms of the O(α t α 3 s ) have been obtained (see the updated version of the public code Himalaya) [77]. The theoretical uncertainty on the lightest CP-even Higgs-boson mass within the MSSM from unknown higher-order contributions is still at the level of about 2 − 3 GeV for scalar top masses at the TeV-scale, where the actual uncertainty depends on the considered parameter region [62,[78][79][80][81][82].
In the NMSSM the full one-loop calculation including the momentum dependence has been performed in the DR renormalization scheme in Ref. [83,84], or in a mixed on-shell (OS)-DR scheme in Ref. [85][86][87]. Dominant two-loop contributions of O(α t α s , α 2 t ) have been calculated in the leading logarithmic approximation [88,89], and of O(α t α s , α b α s ) in the DR scheme in the EP approach [83]. The two-loop corrections involving only superpotential couplings were given in Ref. [90]. A two-loop calculation of the O(α t α s ) corrections with the top/stop sector renormalized in the OS scheme or in the DR scheme was provided in Ref. [91], while the two-loop corrections of O(α 2 t ) in the CP-violating NMSSM were calculated in a mixed OS-DR scheme [92]. These contributions are implemented in the public code NMSSMCalc. A consistent combination of a full one-loop calculation with all corrections beyond one-loop in the MSSM approximation was given in Ref. [87]. According to a comparison of the various two-loop contributions, at present the theoretical uncertainties from unknown higher-order corrections in the NMSSM are expected to be still larger than for the MSSM [92][93][94].
Beyond the MSSM and the NMSSM, only generic DR-calculations of Higgs-boson mass corrections exist publicly available. An automated calculation of the full one-loop corrections, supplemented by partial two-loop corrections to neutral scalars [90] is implemented in the Mathematica package SARAH [95,96], which can be used to produce a spectrum generator based on the public code SPheno [97]. A hybrid Higgs-boson mass calculation combining effective field theory and fixed-order calculations for a generic class of SUSY models is publicly available in the code FlexibleSUSY [98], also using the expression for the renormalization group equations and fixed-order self-energies as they are calculated by SARAH.
In a previous publication we presented the first calculation of radiative corrections to the neutral scalars in a mixed OS-DR scheme for the µνSSM with only one generation of right-handed neutrinos [99]. We described in detail the renormalization of the scalar potential, including the full one-loop quantum corrections. We consistently combined the full one-loop corrections with the leading MSSM-like two-loop contributions using FeynHiggs. We showed that the contributions from the (s)top-and the (s)bottom-sector are also dominant in the µνSSM, therefore proving that the combination of the one-loop result together with the two-loop contributions from FeynHiggs provides a calculation of the SM-like Higgs-boson mass at a similar accuracy as the NMSSM prediction. In this work, we go one step beyond and extend the calculation to the full µνSSM with three generations of right-handed neutrinos. A striking difference between the one-and the threegeneration case is that in the latter case the neutrino sector can be described in agreement with experimental results without having to rely on the radiative generation of neutrino masses. On account of this, we are able to present benchmark scenarios of the µνSSM accurately accommodating a SM-like Higgs boson at ∼ 125 GeV, as well as correct neutrino mass differences and mixing angles. In addition, we show that it is possible to simultaneously explain two excesses measured at LEP and CMS at a mass of ∼ 96 GeV at the 1σ level.
The paper is organized as follows. In Sect. 2 we describe the model and explain the particle mixings in each sector. In Sect. 3 we give details about the renormalization of the neutral scalar potential at the one-loop level, including the full set of free parameters of the µνSSM. We present the renormalization conditions applied to extract the parameter counterterms, either in the neutral scalar or the neutral fermion sector. In Sect. 4 we explain the extraction of the one-loop corrections to the CP-even scalar masses, based on the renormalization prescription introduced before. We also describe the incorporation of higher-order contributions from the MSSM. In Sect. 5 we discuss a set of benchmark scenarios with several light Higgs bosons. We conclude in Sect. 6.

The model: µνSSM with three generations of right handed neutrinos
The superpotential of the µνSSM with three generations of right-handed neutrinos is written as are the left-chiral quark and lepton superfield doublets, andû c i ,d c i ,ê c i and ν c i are the right-chiral quark and lepton superfields. i and j are family indices running from one to three and a, b = 1, 2 are indices of the fundamental representation of SU(2) with ab the totally antisymmetric tensor and ε 12 = 1. The color indices are not written out. Y u , Y d and Y e are the usual Yukawa couplings also present in the MSSM. The trilinear singlet self couplings κ ijk and the trilinear coupling with the Higgs doublets λ i in the second row are analogues to the couplings of the singlet in the superpotential of the Z 3 -symmetric NMSSM. The µ-term is generated dynamically after the spontaneous EWSB, when the right-handed sneutrinos obtain a vev. The κ-term forbids a global U(1) symmetry avoiding the existence of a Goldstone boson in the CP-odd sector. The remarkable difference to the NMSSM is the additional Yukawa coupling Y ν which induces explicit breaking of R-parity through the λ-and κ-terms, and which justifies the interpretation of the singlet superfields as right-handed neutrino superfields. It should be pointed out that in this case lepton number is not conserved anymore, and also the flavor symmetry in the leptonic sector is broken. A more complete motivation of this superpotential can be found in Refs. [22,24,29].
Working in the framework of low-energy SUSY the corresponding soft SUSY-breaking Lagrangian is given by In the first four lines the fields denote the scalar components of the corresponding superfields. In the last line the fields denote the fermionic superpartners of the gauge bosons. We will neglect flavor mixing in the squark and the quark sector, so the soft masses will be diagonal and we write m 2 where the summation convention on repeated indices is not implied, and the quark Yukawas is a 3-dimensional vector in family space, which is always regarded to be absent in the tree-level Lagrangian of the µνSSM, because it spoils the electroweak seesaw mechanism. We include it here, because the operator is generated at the one-loop level and the parameters m 2 H d L are needed to renormalize the scalar potential. The same reasoning applies for the non-diagonal elements of the soft slepton masses (m 2 L ) i =j and (m 2 ν ) i =j . Theoretically, the absence of soft mass parameters mixing different fields at tree level, (m 2 ) i =j , etc., can be justified by the diagonal structure of the Kähler metric in certain supergravity models, or when the dilaton field is the source of SUSY breaking in string constructions [29]. Notice also that when the down-type Higgs doublet superfield is interpreted as a fourth family of leptons the parameters m 2 H d L can be seen as off-diagonal elements of m 2 L [100]. After the EWSB the neutral scalar fields will acquire a vev. This includes the left-and righthanded sneutrinos, because they are not protected by lepton-number conservation as in the MSSM and the NMSSM. We define the decomposition 3) in such a way that after the EWSB they develop the real vevs which is valid assuming CP-conservation, as we will do throughout this paper. The neutral scalar potential of the µνSSM is given at tree level with all parameters chosen to be real by the soft terms and the F -and D-term contributions of the superpotential. One finds (2.11)

The neutral scalar sector
Using the decomposition from Eqs. (2.4) -(2.6) the linear and bilinear terms in the fields define the tadpoles T ϕ and the scalar CP-even and CP-odd neutral mass matrices m 2 ϕ and m 2 σ after electroweak symmetry breaking, where we collectively denote with the CP-even and CP-odd scalar fields, respectively. The linear terms are only allowed for CP-even fields and given by where the h i are the CP-even scalar fields in the mass eigenstate basis. Without CP-violation in the scalar sector the matrix U H is real. Similarly, for the CP-odd scalars we define the rotation matrix U A that diagonalizes the mass matrix m 2 σ , which includes the neutral Goldstone boson A 1 = G 0 . Because of the smallness of the neutrino Yukawa couplings Y ν ij , which also implies that the left-handed sneutrino vevs v iL have to be small, so that the tadpole coefficients vanish at tree level [24], the mixing of the left-handed sneutrinos with the doublet fields and the singlets will be small.
It is a well known fact that the quantum corrections to the Higgs potential are highly significant in supersymmetric models, see e.g. Refs. [78,101,102] for reviews. As in the NMSSM [18], the upper bound on the lowest Higgs mass squared at tree level is relaxed through additional contributions from the right-handed sneutrinos [24]; Nevertheless, quantum corrections were still shown to contribute significantly especially to the prediction of the SM-like Higgs-boson mass [85,87,94,99,[103][104][105][106][107]. In Ref. [99] we already investigated how important the unique loop corrections of the µνSSM beyond the NMSSM are in realistic scenarios, considering only one generation of right-handed neutrinos, finding only negligible differences compared to the NMSSM-like corrections. This is related to the small size of the neutrino Yukawas Y ν compared to the other couplings in the superpotential. In this paper we go beyond Ref. [99] and investigate the complete µνSSM with three right-handed neutrinos. Thus, genuine effects from the µνSSM are guaranteed to play a role in the prediction of the SM-like Higgs-boson mass just by the presence of additional singlets, whose couplings to the Higgs doublet fields are not suppressed by the size of Y ν . Furthermore, the model can accommodate neutrino data at tree level, so we will be able to describe the phenomenology related to both the scalar and the fermionic sector (and their interplay) more precisely.
If the mixing of the CP-even and -odd sneutrinos with the doublet fields is small, one can obtain approximate analytical expressions for the tree-level masses of the sneutrinos. For the left-handed sneutrinos the dominant terms are proportional to the inverse of their vevs. In particular, for the CP-even left-handed sneutrinos one finds for the diagonal entry of the mass matrix m 2 ϕ , where we defined the effective µ-term as Note that the first term in Eq. (2.23) can usually be neglected as long as v iR v d , v u . Each CP-odd left-handed sneutrino is nearly degenerate with the corresponding CP-even one, though they are slightly lighter due to different contributions proportional to the gauge couplings, For the CP-even right-handed sneutrinos we find, under the assumptions that non-diagonal elements of κ ijk and (m 2 In very good approximation this matrix has the form The eigenvalues of such a matrix are a − b , a − b and a + 2b, and only the mass eigenstate corresponding to the latter eigenvalue mixes with the SM-like Higgs boson [108]. Later we will make use of this fact to simplify the accommodation of SM-like Higgs boson properties, even when all three right-handed neutrinos have masses close to or even below 125 GeV, because two righthanded sneutrinos will be conveniently decoupled from the remaining scalars, and interact very weakly with SM particles. In extensions of the NMSSM with several gauge singlets, this decoupling can lead to practically stable particles. In the µνSSM this is not possible, because the decoupling cannot be exact, even when κ ijk and (m 2 ν R ) ij are diagonal. We stress that the universality of the κ iii is stable with respect to the RGE running, if also the λ i parameters are universal. In this case, one can deduce from the explicit form of the one-loop counterterms δκ ijk shown in the App. B.2.2 that differences in the running are exclusively generated by terms proportional to (Y ν ) 2 , which are negligible in realistic scenarios. For v iR v d , v u Eq. (2.26) further simplifies to while for the CP-odd right-handed sneutrinos one finds Thus, to avoid tachyons both in the CP-even and -odd scalar spectrum, we will follow the sign convention A κ iii < 0, κ iii > 0 and v iR > 0. Before we come to the one-loop renormalization of the neutral scalar potential we briefly describe the other relevant sectors of the µνSSM.

Squark sector
The numerically most important one-loop corrections to the scalar potential are expected from the stop/top-sector, analogous to the (N)MSSM [104][105][106][109][110][111], due to the huge Yukawa coupling of the (scalar) top. The tree-level mass matrices of the squarks differ slightly from the ones in the MSSM. Neglecting flavor mixing in the squark sector, one finds for the up-type squark mass matrix M u i of generation i, It should be noted that in the non-diagonal element explicitly appear the neutrino Yukawa couplings. This term arises in the F-term contributions of the squark potential through the quartic coupling of up-type quarks and one left-handed and the right-handed sneutrinos after EWSB. The mass eigenstates u i1 and u i2 are obtained by the unitary transformation Similarly, for the down-type squarks it is The mass eigenstates d i1 and d i2 are obtained by the unitary transformation (2.37)

Charged scalar sector
Since R-parity, lepton-number conservation and lepton-flavor universality are broken, the six charged left-and right-handed sleptons mix with each other and with the two charged scalars from the Higgs doublets. In the basis C T = (H − d * , H + u , e * iL , e * jR ) we find the following mass terms in the Lagrangian: Assuming CP conservation m 2 H + is a symmetric matrix of dimension 8, The entries are given in App. A.3. The mass matrix is diagonalized by an orthogonal matrix U + : where the diagonal elements of m 2 H + diag are the squared masses of the mass eigenstates which include the charged Goldstone boson H + 1 = G ± .

Charged fermion sector
The charged leptons mix with the charged gauginos and the charged higgsinos. Following the notation of Ref. [29], we write the relevant part of the Lagrangian in terms of two-component spinors ( The 5 × 5 mass matrix m e is defined by The mass matrix is diagonalized by two unitary matrices U e L and U e R , where m diag e contains the masses of the charged fermions in the mass eigenstate base Note that terms mixing the SM leptons and the MSSM-like charginos in Eq. (2.43) are suppressed by the size of the left-handed vevs v iL or the neutrino Yukawa couplings Y ν ij . The smallness of v iL in comparison to the other vevs and M 2 assures the decoupling of the three leptons from the Higgsino and the wino, prohibiting substantial lepton-flavor-universality and lepton-number violation in the charged fermion sector.

Neutral fermion sector
The three left-handed neutrinos and the right-handed neutrinos mix with the neutral Higgsinos and gauginos. Again, following Ref. [29], we write the relevant part of the Lagrangian in terms of two-component spinors ( where m ν is the 10 × 10 symmetric mass matrix. The neutral fermion mass matrix is determined by Because of the Majorana nature of the neutral fermions we can diagonalize m ν with the help of just a single -but complex -unitary matrix U V , with where λ 0 are the two-component spinors in the mass basis. The eigenvalues of the diagonalized mass matrix m diag ν are the masses of the neutral fermions in the mass eigenstate basis. The mass matrix has a seesaw structure, assuring that the three lightest eigenvalues will be very small, so that the mass eigenstates λ 0 1,2,3 can practically be identified with the SM left-handed neutrinos. Components from the MSSM-like neutralinos and the right-handed neutrinos are negligible for the three lightest states. Thus, the left-handed neutrino mixing can in very good approximation (using diagonal Y e ) be expressed in the usual PMNS formalism [112,113] by the three mixing angles θ 12 , θ 13 with i, j = {1, 2, 3}, and we used the short-hand notation c x = cos θ x and s x = sin θ x . In our numerical studies we fitted the experimentally well measured quantities (2.52) We restricted ourselves in the neutrino sector to a tree-level analysis, because the one-loop corrections turn out to be moderate in size (in the normal hierarchy pattern) [28] and can always be compensated by a small shift in the neutrino Yukawa couplings Y ν ij without affecting the conclusions drawn in the scalar sector, in particular for the observables related to the SM-like Higgs boson.
To reduce the parameter space in our analysis, we usually assume the couplings Y ν to be diagonal, as we do for the lepton Yukawa couplings Y e . We emphasize that non-diagonal Y ν ij are not required to reproduce the correct neutrino mixing, because sizable flavor mixing is always present after the diagonalization of m ν , generated by the mixing terms of the left-handed neutrino states with the gauginos, Higgsinos and right-handed neutrinos. Quantitatively, this can be illustrated assuming universal parameters λ : Eq. (2.54) demonstrates that substantial flavor mixing is practically unavoidable in the µνSSM. The first two terms are of particular importance. The first term can be attributed to the mixing with the right-handed neutrinos and higgsinos, and the other terms also include the gaugino mixing. Note that for moderate values of tan β and not too small values of λ the first two terms are the dominant contributions. They contain diagonal and non-diagonal contributions that can easily be adjusted by an appropriate choice of the parameters Y ν i , v iL and the soft gaugino mass parameters M 1 and M 2 . These parameters play only a minor role in the predictions for the SM-like Higgs boson mass and its mixing with the right-handed sneutrinos. Thus, the above mentioned parameters will be used to reproduce neutrino physics in agreement with experimental limits, without having to worry about spoiling the properties of the SM-like Higgs boson.

Renormalization of the Higgs potential at One-Loop
At tree level the part of the Higgs potential relevant for the masses of the scalars is given by the tadpole coefficients in Eqs. (2.13)-(2.16) and the CP-even and CP-odd scalar mass matrix elements in Eqs. (2.17) and (2.18). We want to employ a renormalization procedure as close as possible to the ones used in the (N)MSSM. Therefore, we define in the following subsection certain replacements to obtain a new set of free parameters. The new set of free parameter will permit us to make use of a mixed On-Shell (OS)/DR renormalization scheme. The precise definition of the counterterms of the free parameters will be given in Sect. 3.2. Finally, we describe the renormalization conditions applied on each parameter and the extraction of the counterterms in Sect. 3.3.

Parameter replacements
The vevs of the doublet Higgs fields v u and v d are substituted by the MSSM-like parameters tan β and v according to Note that the definition of v 2 differs from the one in the MSSM by the term v iL v iL . This allows to maintain the relations between v 2 and the gauge-boson masses as they are in the MSSM. Numerically, the difference in the definition of v 2 is negligible. Maintaining the functional form of tan β as it is in the (N)MSSM is convenient to facilitate the comparison of the quantum corrections in the µνSSM and the (N)MSSM, as the one-loop counterterm of tan β is expressed without having to include the counterterms for the left-handed sneutrino vevs [99]. The gauge couplings g 1 and g 2 will be replaced by the gauge-boson masses M W and M Z , This is reasonable, because the gauge-boson masses are well measured physical observables, so we can define them as OS parameters. The explicit dependence of the quantum corrections on the mass counterterm for M 2 W drops out at the one-loop level, but it will contribute implicitly in the definition of the counterterm for v 2 . The scalar soft masses m 2 H d and m 2 Hu and the diagonal elements of the soft slepton mass matrices m 2 L and m 2 ν are replaced by the tadpole coefficient in which they appear. Using the tadpole coefficients as input parameters facilitates the absorption of quantum corrections that would spoil the true vacuum of the potential. Alternatively, one could also trade the vevs for the tadpole coefficients, and keep the soft masses as input parameters. However, it is computationally much more convenient to use the vevs as input and solve the tadpole equations for the squared soft masses, because they appear linearly, while solving the tadpole equations for the vevs, using the soft masses as input, is a complex non-linear problem with multiple solutions. The complete set of independent parameters is summarized in Tab. 1.

Counterterms
The entries of the neutral scalar mass matrices are functions of the independent parameters, and we define their renormalization as  Table 1: Set of independent parameters initially entering the tree-level Higgs potential of the µνSSM in the first row, and final choice of free parameters after the substitutions defined in the text.
The mass counterterms δm 2 ϕ and δm 2 σ enter the renormalized one-loop scalar self-energies. They are given as a linear combination of the counterterms of the independent parameters, which we define as The divergent parts of the counterterms are fixed to cancel the UV divergences. The finite pieces, and thus the meaning of the parameters, have to be fixed by renormalization conditions. We will adopt a mixed renormalization scheme, where tadpoles and gauge boson masses are fixed OS, and the other parameters are fixed in the DR scheme. The exact renormalization conditions will be given in Sect. 3.3. The dependence of the mass counterterms δm 2 ϕ and δm 2 σ on the counterterms of the free parameters is given at the one-loop level by the first order expansion w.r.t. the free parameters, We define the mixing matrices U H and U A to diagonalize the renormalized mass matrices, so they do not have to be renormalized. The expressions for the counterterms of the scalar mass matrices in the mass eigenstate basis are then given by where we emphasize that δm 2 h and δm 2 A are not diagonal, as they would be in a purely OS renormalization procedure which is often used in theories with flavor mixing [114].
The field renormalization required to obtain finite scalar self-energies at arbitrary momentum, is defined by where √ Z and δZ are 8 × 8 dimensional matrices and the equal sign is valid at the one-loop level. In contrast to the (N)MSSM, the field renormalization is not diagonal in the interaction basis. The reason is that the µνSSM explicitly breaks lepton-number conservation and leptonflavor universality, resulting in kinetic mixing terms at one-loop order. 2 For the CP-even and CP-odd neutral scalar fields the definition in Eq. (3.10) implies the following field renormalization in the mass eigenstate basis:

Renormalization conditions
In this section we briefly describe our choice for the renormalization conditions. We start with the OS conditions for the gauge boson mass parameters and the tadpole coefficients followed by our definitions for the DR renormalized parameters, including the field renormalization. All counterterms are extracted diagrammatically by calculating one-loop corrections to linear, bilinear or trilinear terms of the Lagrangian, and identifying the part of the corrections that had to be absorbed individually by the counterterms of the parameters appearing in the tree-level expression of the term. We generated the Feynman diagrams using our FeynArts [116] model file, which was initially created with SARAH version 4.12.0 [95,117]. We modified the model file by hand to be able to use FormCalc [118] for further evaluations and to improve the analytical and numerical evaluation of the rather large expressions. Since the divergent parts of one-loop counterterms can in principle also be derived from the one-loop beta functions, for which generic analytical formulas exist [119][120][121][122][123][124], the diagrammatic calculation of the counterterms was an excellent test for the correctness of our FeynArts model file.
The determination of the counterterms for the set of independent parameters was done in a specific order, because in some cases the definition of the renormalization condition of one counterterm depends on other counterterms, that necessarily had to be determined before. In Fig. 1 we give an overview of the strategy for the extraction of the counterterms. We also highlight in color the sectors of the µνSSM in which the corresponding counterterm was extracted (see caption). The exact definition of the counterterms and their final analytic expressions in terms of UV divergences for the DR counterterms are listed in App. B.
Therein, divergent parts are expressed proporional to ∆, where loop integral are solved in 4 − 2ε dimensions and γ E = 0.5772 . . . is the Euler-Mascharoni constant. Since the field renormalization constants contribute only via divergent parts, they do not contribute to the finite result after canceling divergences in the self-energies. As regularization scheme we choose dimensional reduction [125,126] which was shown to be SUSY conserving at the one-loop level [127]. In contrast to the OS renormalization scheme our field renormalization matrices are hermitian. This holds also true for the field renormalization in the mass eigenstate basis, because as already mentioned the rotations in Eq. (2.19) and Eq. (2.21) diagonalize the renormalized tree-level scalar mass matrices, so Eqs. (3.12) do not introduce non-hermitian parts into the field renormalization that would have to be canceled by a renormalization of the mixing matrices U H and U A themselves.

OS conditions
The SM gauge boson masses are renormalized OS requiring whereΣ T stands for the transverse part of the renormalized gauge boson self-energy. For their mass counterterms these conditions yield Here the Σ T (without the hat) denote the transverse part of the unrenormalized gauge boson self-energies. For the tadpole coefficients T ϕ i the OS conditions read where T (1) ϕ i are the one-loop contributions to the linear terms of the scalar potential, stemming from tadpole diagrams shown in Fig. 2. The tadpole diagrams are calculated in the mass eigenstate basis h. The one-loop tadpole contributions in the interaction basis ϕ are then obtained by the rotation Accordingly we find for the one-loop tadpole counterterms

DR conditions
For practical purposes we decided to renormalize all remaining parameters in the DR scheme, reflecting the fact that there are no physical observables yet that could be directly related to them. The counterterms of each parameter were obtained by calculating the divergent parts of one-loop corrections to different scalar and fermionic two-and three-point functions. We sketch the determination of the counterterms in the (possible) order in which they can be successively derived (see Fig. 1). The general strategy for extracting the counterterms of the free parameters is the following. At first one finds a relatively simple tree-level expression, containing the parameter whose counterterm one wants to extract, and, apart from that, exclusively parameters whose countertems are already known. In our case we used bilinear and trilinear couplings in the neutral scalar and fermionic sector for this. Then one calculates the one-loop corrections to the term by evaluating the corresponding Feynman diagrams. As we only need the divergent parts of the loop corrections for DR conditions, we are able to calculate the diagrams in the gauge eigenstate basis, where the number of diagrams is  drastically reduced. Once the divergent contributions are known, the counterterm to be identified directly follows from the expression of the renormalized Green's functions.
In this chapter we only state the general formulas for the renormalized two-and three-point Green's functions. The exact conditions used to extract each countertem is listed in App. B.2. Therein, we also show the resulting analytic expressions for the counterterms renormalized in the DR scheme. 3 A more detailed discussion in the case of the µνSSM with one generation of righthanded neutrinos can be found in Ref. [99].

Neutral fermion sector
We derived most of the counterterms in the neutral fermion sector. The mass matrix elements of the neutral fermions (m ν ) ij get one-loop corrections via the neutral fermion self-energies χ 0 i χ 0 j , that for Majorana fermions can be decomposed as Defining for the renormalized mass matrix the renormalized scalar part of the self-energies at zero momentum is given bŷ The field renormalization constants can be obtained in the DR scheme by calculating the divergent part of the fermionic piece, The divergent parts of the self-energies of the neutral fermions were calculated diagrammatically in the interaction basis, where diagrams with mass insertions have to be included. If (m ν ) ij containes just a single parameter whose counterterm is unknow, Eq. (3.21) provides a definition for the missing counterterm once the mass counterterm δ (m ν ) ij is expressed in terms of the counterterms of the fundamental parameters. Unfortunately, the right-handed vevs v iR always appear in sums over the family index and never isolated. Therefore we calculated loop corrections to the three elements (m ν ) i,7 = −Y ν ij v jR / √ 2, which provides us with a linear system of three independent equations; that can be solved analytically for the three counterterms δv 2 iR = 2v iR δv iR .

Neutral scalar trilinear couplings
General one-loop scalar three-point functions can be renormalized by wavefunction counterterms and the specific vertex counterterm aŝ where δZ ij are the scalar field renormailzation constants defined in Eq. (3.10), Γ ϕ i ϕ j ϕ k are the treelevel couplings, Γ ϕ i ϕ j ϕ k are the one-loop corrections obtained by evaluating the non-irreducible one-loop three-point diagrams, and δΓ (1) ϕ i ϕ j ϕ k is the coupling counterterm given as a function of the counterterms of the independent parameters. The field renormalization constants are defined as DR-parameters. We calculate the UV-divergent part of the derivative of the scalar CP-even self-energies in the interaction basis and define As before, if at tree level Γ ϕ i ϕ j ϕ k just contains a single parameter whose counterterm is still unknown, the counterterm can be extracted from the divergent part of the one-loop corrections Γ (1) ϕ i ϕ j ϕ k demanding that the renormalized quantity is finite. Similarly to the vevs v iR , for the parameters κ iij it is not possible to find a tree-level expression where each element appears isolated. However, using the renormalized expression in Eq. (3.24) for the vertex we can extract the counterterms for the three subsets (δκ 11j , δκ 22j , δκ 33j ) by renormalizing the subset of vertices (Γ ν R ). Thus, for each subset (κ iij ; i = 1, 2, 3) we get a linear system of three (j = 1, 2, 3) equations to extract the counterterms δκ iij from the condition that the renormalized one-loop three-point function is finite.

Neutral scalar masses
The soft scalar masses appear in the bilinear terms of the Higgs potential. They can be renormalized by calculating radiative corrections to scalar self-energies. Since our final aim is to obtain loop corrections for the CP-even scalars, we used the CP-odd scalar sector to extract the counterterms of the soft masses to have an independent crosscheck of both neutral scalar sectors.
The general form of the renormalized scalar self-energies at the one-loop level iŝ where X = (ϕ, σ) represents either the CP-even or the CP-odd scalar fields and we made use of the fact that the field renormalization constants δZ and the mass matrix m 2 X are real. Demanding that the renormalized self-energiesΣ A i A j are finite in the mass eigenstate basis we can define the divergent parts of the mass counterterms via where the field counterterms in the mass eigenstate basis were defined in Eq. (3.12), and renormalized as DR parameters like the ones for the CP-even scalars (see Eq. (3.25)), and the masses squared m 2 A i are the eigenvalues of the diagonal CP-odd scalar mass matrix m 2 A . In Fig. 3 we show the diagrams that have to be calculated to obtain the quantum corrections to scalar self-energies at the one-loop level in the mass eigenstate basis.
We calculated all diagrams in the 't Hooft-Feynman gauge in which the Goldstone bosons A 1 and H ± 1 and the ghost fields u ± and u Z have the same masses as the corresponding gauge bosons. Calculating the CP-odd self-energies Σ A i A i diagrammatically, we obtain the mass counterterms in mass eigenstate basis through the Eq. (3.28). Now inverting the rotation in Eq. (3.9) yields the mass counterterms for the CP-odd self-energies in the interaction basis, (3.29) Analytically, following the expansion in Eq. (3.8), some of the mass counterterms δm 2 σ depend on the counterterms of the soft mass parameters. We use this dependences to extract the counterterms of (m 2

Loop corrected scalar masses
In the previous section we have defined an OS/DR renormalization scheme for the µνSSM neutral scalar sector. This can be applied (via the FeynArts model file, in which the counterterms are implemented) to any higher-order correction in the µνSSM. As a first application, we evaluate the full one-loop corrections to the CP-even scalar sector in the µνSSM. In the following we emphasize the differences w.r.t the analysis with just one right-handed neutrino from Ref. [99].

Evaluation at the one-loop level
The one-loop renormalized self-energies in the mass eigenstate basis are given bŷ with the field renormalization constants δZ H and the mass counter terms δm 2 h in the mass eigenstate basis defined by the rotations in Eq. (3.12) and Eq. (3.9), respectively. Σ h i h j is the unrenormalized self-energy obtained by calculating the diagrams shown in Fig. 3 with the CP-even states h on the external legs. The self-energies were calculated in the 't Hooft-Feynman gauge, so that gaugefixing terms do not yield counterterm contributions in the Higgs sector at the one-loop level. The loop integrals were regularized using dimensional reduction [125,126] and numerically evaluated for arbitrary real momentum using LoopTools [118]. The contributions from complex values of p 2 were approximated using a Taylor expansion with respect to the imaginary part of p 2 up to first order.
In Eq. (4.1) we already made use of the fact that δZ H is real and symmetric in our renormalization scheme. The mass counterterms are defined as functions of the counterterms of the free parameters following Eq. (3.8) and Eq. (3.9). They contain finite contributions from the tadpole counterterms and from the counterterm for the gauge boson mass M 2 Z . The matrix δm 2 h is real and symmetric.
The renormalized self-energies enter the inverse propagator matrix The loop-corrected scalar masses squared are the zeroes of the determinant of the inverse propagator matrix. The determination of corrected masses has to be done numerically when one wants to account for the momentum-dependence of the renormalized self-energies. This is done by an iterative method that has to be carried out for each of the six squared loop-corrected masses. 4

Inclusion of higher orders
In Eq. (4.2) we did not include the superscript (1) in the self-energies. Restricting the numerical evaluation to a pure one-loop calculation would lead to very large theoretical uncertainties. These can be avoided by the inclusion of corrections beyond the one-loop level. Here we follow the approach of Refs. [87,99] and supplement the µνSSM one-loop results by higher-order corrections in the MSSM limit as provided by FeynHiggs (version 2.13.0) [40,48,[60][61][62][63][64]67]. 5 In this way the leading and subleading two-loop corrections are included, as well as a resummation of large logarithmic terms, see the discussion in Sect. 1, In the partial two-loop contributionsΣ (2 ) h we take over the corrections of O(α s α t , α s α b , α 2 t , α t α b ), assuming that the MSSM-like corrections approximate numerically well the corresponding µνSSM corrections. This assumption is reasonable since the only difference between the squark sector of the µνSSM in comparison to the MSSM are numerically suppressed terms proportional to v jL Y ν jk v kR in the non-diagonal element of the up-type squark mass matrices (see Eq. (2.31)) and proportional to v iL v iL in the diagonal elements of the up-and down-type squark mass matrices (see Eqs. (2.30), (2.32), (2.34) and (2.36)). Furthermore, in Refs. [87,107] the quality of the MSSM approximation was tested in the NMSSM, showing that the genuine NMSSM contributions are in most cases subleading. Our results in Ref. [99] confirm that the same holds true in the µνSSM for the SM-like Higgs-boson mass. The same is expected for the contributions stemming from the resummation of large logarithmic terms given byΣ resum h .

Numerical analysis
In the following we will present several benchmark points (BPs) that illustrate the phenomenology of the scalar sector of the µνSSM. We concentrate on scenarios in which a right-handed sneutrino is mixed with the SM-like Higgs boson. By setting κ ijk = κ δ ij δ jk , as explained in Sect. 2.1, we achieve that only a single right-handed sneutrino substantially mixes with the SM-like Higgs boson. Naturally, the mass scale of the right-handed sneutrinos will then be of the order of the SM-like Higgs boson. However, scenarios in which the decay of the SM-like Higgs boson to two right-handed sneutrinos is kinematically allowed, are experimentally very constrained [108].
In contrast to most of the previous studies of the µνSSM with three generations of righthanded neutrinos [24,25,27,28, 108], we will not always make the simplifying assumption that genuine low-energy µνSSM-parameters have universal values independent of the family index. In Sect. 5.3 we elaborate on the effect of non-universal λ i on the SM-like Higgs-boson mass, while keeping λ λ λ 2 = λ i λ i constant. Since we know from Eq. (2.22) that the tree-level mass of the SM-like Higgs boson strongly depends on λ λ λ 2 , it will be discussed whether the loop corrections increase the dependence on the individual values λ i .
We consider the following experimental constraints on the scenarios presented: -We use the public code HiggsBounds v.5.2.0 [129][130][131][132][133] to determine whether a BP has been excluded by cross section limits from Higgs searches at LEP, LHC or Tevatron. These searches are mostly sensitive to the heavy Higgs and the right-handed sneutrinos, if these are substantially mixed with the SM-like Higgs boson. The production of the left-handed sneutrinos is much smaller at the LHC, and signals from their decay usually demand dedicated searches [30], especially if the left-handed sneutrino is the LSP [29,31].
-The properties of the SM-like Higgs boson, i.e., its mass and signal rates at LHC and Tevatron, are checked using the public code HiggsSignals v.2.2.1 [134][135][136]. Here we assume a theoretical mass uncertainty of 3 GeV. HiggsSignals provides us with a χ 2 -analysis of n obs = 106 observables in the 7+8 TeV data package and n obs = 101 observables in the 13 TeV data package. In our plots we show the reduced χ 2 red = χ 2 /n obs , where a value of χ 2 red = 1 means that on average the signal rates of the SM-like Higgs boson are at the level of the ±1σ range of the measurements.
-The properties of the neutrino sector are in agreement with the measurements of the masssquared differences and the mixing angles obtained from neutrino oscillation experiments. We check that our predictions are within the ±3σ bands published by the NuFit collaboration [ where we considered the normal mass ordering which is now favored by experiments [139]. A genetic algorithm was used to find parameter points that minimize the sum of squared deviations between theoretical prediction and experimental values specified above [140]. Even though the µνSSM allows for flavor-violating decays of leptons, the existing experimental bounds (for instance on µ → eγ) are automatically fulfilled when the constrains on neutrino masses are taken into account [26].
For the necessary input of HiggsBounds and HiggsSignals we compute the decays of the scalars at leading order, but with the loop-corrected mixing matrix elements inserted in the expressions of the scalar couplings. In the limit of vanishing external momentum, which we used in the determination of the mixing matrix elements for the couplings, this method corresponds to include the finite wave-function renormalization factors (Z-factors) for each external scalar [63,141]. For loopinduced decays and off-shell decays to vector bosons we implemented analytic results from the MSSM well known in the literature [142][143][144][145][146], and scaled the expressions with effective couplings defined by the mixing matrix elements and tan β to obtain the result for the scalars in the µνSSM.
For the coupling to b quarks we included the running bottom mass and for the decay to gluons the tan β  Table 2: Low-energy values for the parameters, as defined in the text, of the scan over λ. Dimensionful parameters are given in GeV. The parameters in the last row are fitted to neutrino oscillation data.
running of α s from M Z to the mass of the decaying scalar, and finally add leading higher-order QCD corrections [144,147]. As described in Sect. 2.5 we use the left-handed vevs v iL , the soft gaugino masses M 1 and M 2 , and the neutrino Yukawa couplings Y ν ij to fit the neutrino masses and mixings accurately, making use of the fact that they can be modified without spoiling the properties of the SM-like Higgs boson. Besides for the scenario presented in Sect. 5.4, it will be sufficient to just consider diagonal non-zero elements of Y ν ij . Because we concentrate here on the scalar sector of the µνSSM, and since the fitting has to be done numerically, we do the fitting in our scans just in one particular point for each analysis. By varying a parameter, the prediction for the neutrino properties can be outside the experimentally allowed range in some points. We indicate in our plots when this is the case. Since the neutral fermion mass matrix is of dimension 10, with large hierarchies between the neutrino sector and the remaining part, including one-loop corrections is time-consuming and numerically very challenging. Therefore we stick to a tree-level analysis for the neutrinos. However, we checked for several points that the one-loop corrections are sub-leading and can in principle be compensated by a slight change of the parameters. 6

Scan over λ
The first scenario we are presenting is one with a light right-handed µ-sneutrino that mixes substantially with the SM-like Higgs boson. We show the chosen parameters in Tab. 2. To simplify the notation we define λ = λ i , A λ = A λ i , A ν = A ν ii , κ = κ iii and A κ = A κ iii and vanishing otherwise.

The soft parameters are given in terms of
and m e = m e ii and vanishing otherwise. We vary over the universal parameter λ, while keeping the remaining parameters fixed. For the right-handed e-and τ -sneutrino vevs we chose v 1,3R = 1 TeV, but set a smaller value of v 2R = 765 GeV for the µ-sneutrino vev to decrease the mass of the CP-even µ-sneutrino to the range around the SM-like Higgs-boson mass. The choice to pick ν µR as the light right-handed sneutrino is of no relevance. The large value of κ = 0.5 assures that the other two right-handed sneutrinos will have masses between 300 and 400 GeV, well above 125 GeV. Because the SM-like Higgs boson mass will get additional contributions from the mixing with ν R µR , tan β can be chosen rather low.
As mentioned in the previous subsection, we fit the properties of the neutrinos in just one particular point of the parameter scan. In this scenario, this was done for λ = 0. 168 Tree-level One-loop Two-loop Figure 4: CP-even scalar spectrum in the scan over λ at tree level, one-loop level and partial two-loop level. We show in the brackets the dominant composition of the tree-level, one-loop and two-loop mass eigenstates h (0) , h (1) and h (2) , in the experimentally allowed region of the plot. The desired SM-like Higgs-boson mass is indicated with the horizontal green band, assuming a theory uncertainty of 3 GeV. The red regions are excluded by direct searches for additional scalars. In the yellow region the prediction for the mixing angles of the neutrinos lies outside of the 3σ band of the experimental measurement. On top we show χ 2 red for various Higgs-boson signal strength measurements at LHC. trilinear parameters A ν ii to adjust the masses of the left-handed sneutrinos. For the prediction of the masses of the right-handed sneutrinos and the SM-like Higgs boson, the fitted parameters only play a minor role.
In Fig. 4 we show the resulting spectrum of the light CP-even scalars. The remaining CP-even scalars not shown in the plot have masses above 300 GeV and do not play a role in the following discussion. The dotted lines represent the tree-level masses, the dashed lines the masses including the full one-loop corrections, and the solid lines the one-loop + partial two-loop + resummed (referred to as two-loop in the following) corrected masses, as explained in Sect. 4.2. We mark four regions in the plot which are excluded either by HiggsBounds (red), or by not being in agreement with the neutrino oscillation data (yellow). We stress that region 2 is just excluded for the precise choice of parameters shown in Tab. 2. A new fit of the neutrino properties for each value of λ could easily accommodate predictions for the properties of the neutrinos in agreement with experiments. However, since this would exclusively affect the phenomenology of the heavier left-handed sneutrinos in the scalar sector, we do not apply the fit for each value of λ. This spectrum is characterized by the interplay between the light ν R µR and the SM-like Higgs boson. For small λ the two lightest loop-corrected mass eigenstates h 1 and h 2 have roughly an equal amount of H R u -and ν R νR -admixture (see also Fig. 5). Consequently, region 1 is excluded by direct searches at the LHC, because the diphoton resonance search for a SM-like higgs boson excludes h 1 via its decay to photons [148]. At λ ∼ 0.14 the point is reached where the mass of h 1 drops well below 125 GeV. Thus, beyond that point h 1 can be identified with ν R µR , as the doublet-component of h 1 shrinks to values of roughly ∼ 10%. h 2 , on the other hand, sheds its sneutrino admixture, so that it can be identified as the SM-like Higgs boson, and the large quantum corrections from the 1 . Naturally, this mixing will also affect the SM-like Higgs-boson properties. In this way, scenarios like the one shown here will be tested by experiments in two different and complementary ways, both caused by the mixing of ν R µR and the SM-like Higgs boson: Firstly, direct searches for additional Higgs bosons can be applied to ν R µR , because it is directly coupled to SM particles. Secondly, precise measurements of the SM-like Higgs-boson couplings can detect (or exclude) possible variations from SM predictions. To illustrate the possible modifications, we show in Fig. 6 the effective coupling of the two light CP-even scalar mass eigenstates to up-and down-type quarks normalized to the SM-prediction which in good approximation can be expressed via the loop-corrected mixing matrix elements U H(2) ij and β; In the experimentally allowed region the effective coupling of the SM-like Higgs boson to up-type quarks shows deviations of roughly 10%. This is of the order of precision expected by measurements of the SM Higgs boson couplings at the High-Luminosity LHC [149], and (depending on the centerof-mass energy deployed) an order of magnitude larger than the uncertainty expected for these kind of measurements at a possible future e + e − collider like the ILC [150][151][152]. Comparing to Fig. 4 we can see that the region where the effective couplings are closest to one, meaning equal to the SM prediction, does not coincide with the region where the χ 2 red from HiggsSignals is minimized. This is because the mass of the SM-like Higgs boson is slightly too small in this range of λ, so even including a theoretical uncertainty of 3 GeV some signal strength measurements implemented in HiggsSignals are not accounted for by h 2 and χ 2 red becomes worse.

Scan over v R
In Sect. 5.1 we showed that light right-handed sneutrinos with masses in the vicinity of the SMlike Higgs boson are theoretically possible and can induce measurable modifications of the SM-like Higgs-boson properties. Using data of direct searches and measurements of the couplings of the SM-like Higgs boson, the parameter space of these scenarios can be constrained effectively. In this section we present a scenario that is not excluded by current searches in which all three of the CPeven right-handed sneutrinos will have masses below 125 GeV. We chose the parameters appearing in the mass terms of the ν R iR to be universal, i.e., λ : and A κ := A κ iii . As explained at the beginning of Sect. 5, the universality of κ assures that only one of the ν R iR mixes substantially with the SM-like Higgs boson, while the other two are practically decoupled. This makes it easier to control the total admixture of the doublet components of the ν R iR . The complete set of free parameters is shown in Tab. 3. In this scenario we scan over v R , because they appear linearly in the Majorana-like mass terms of the ν R iR , so it is a convenient parameter to control their masses. Compared to the scan over λ in Sect. 5.1 the overall behavior of the SM-like Higgs boson is aligned more to the SM predictions by decreasing λ. Consequently, because at tree level the additional contribution proportional to λ λ λ 2 is smaller, tan β is larger to increase the quantum corrections to the SM-like Higgs-boson mass. We also decrease κ to make the masses of the ν R iR smaller. As before, the parameters in the last row of Tab. 3 were fitted to accurately predict the left-handed neutrino masses and mixings. The fit was done in the point v R = 1226 GeV, but in this case the neutrino data is described accurately over the whole range of v R at tree level.
In Fig. 7 we show the resulting light CP-even scalar spectrum. In the experimentally allowed region (1213 GeV ≤ v R ≤ 1235 GeV) the lightest mass eigenstates are two almost degenerate righthanded sneutrinos. The third right-handed sneutrino is roughly 20 GeV heavier, and it acquires substantial mixing with the SM-like Higgs boson. Naturally, the ν R iR increase their masses when v R becomes larger, but also the SM-like Higgs boson mass increases, because the mixing with the ν R 3R gives additional contributions.
The scenario is excluded experimentally for very small values of v R , because the two lightest mass eigenstates h 1,2 become lighter than half the mass of the SM-like Higgs boson h 4 , so the decays of h 4 into h 1,2 opens up. Experimental searches for the decay of the SM-like Higgs boson into two lighter scalars that subsequently decay into two b-jets and a pair of µ-leptons [153] exclude region 1 in Fig. 7. These additional decay channels of the SM-like Higgs boson are also the reason why the χ 2 red rapidly increases in region 1, because it suppresses ordinary SM-like decays of h 4 . When v R increases above 1235 GeV further constrains from direct searches for additional Higgs bosons and measurement of the properties of the SM-like Higgs boson become relevant. χ 2 red quickly increases above 2 at v R ∼ 1237 GeV. Already at v R ∼ 1235 GeV the scenario is excluded by LEP searches [154]. Note, that in the red region 2 the mixing of ν R 3R with the SM-like Higgs boson enlarges, while ν R 3R is kinematically still in reach of being produced at LEP via the Higgsstrahlungprocess. Consequently, the channel ee → (h 3 )Z → (bb)Z, where h 3 is identified with ν R 3R , excludes this interval. Interestingly, in the experimentally allowed region, where the mass of ν R 3R is even smaller, LEP data cannot rule out this scenario. The reasons for this is not only the smaller mixing of ν R 3R with the SM-like Higgs boson, but also that in the mass range below 100 GeV LEP saw a slight excess over the SM background (see also Sect. 5.4) [154].
Beyond region 2 the current scenario is experimentally excluded by the measurement of the SMlike Higgs-boson mass in the gray region 3 and by the LHC cross section measurement of the process pp → h 1 → γγ in region 4 [155]. This is because in region 3 the cross-over point is reached, in which the masses of the ν R iR become larger than the SM-like Higgs-boson mass. Through the interference effects the SM-like Higgs-boson mass is pushed to lower values beyond that point. In region 4 the mass eigenstate corresponding to the SM-like Higgs boson is the lightest one at just about 118 GeV. Even though there are two scalars in the mass range of the experimentally measured Higgs-boson mass, there is no contribution to any signal-strength measurement at the LHC, reflected by the fact that the χ 2 red is huge in region 4. The reason is that these states correspond to the practically singlet like right-handed neutrino states. The third right-handed sneutrino carrying the doublet admixture taken from the SM-like Higgs boson has a mass of over 140 GeV. Hence, it also does not contribute to signal-strength measurements of the SM-like Higgs boson.
On a side note we briefly discuss the remaining light scalar h 5 in Fig. 7, which is the lefthanded τ -sneutrino at roughly 235-240 GeV. The fit to the neutrino oscillation data generated a hierarchy between the vevs of the left-handed sneutrinos, with v 3L being the largest. As a result, since dominant tree-level contributions to the ν R iL -masses scale with inverse of v iL (see Eq. (2.23)), the ν R τ L is the lightest CP-even left-handed sneutrino. It is rather invisible to usual searches for additional Higgs bosons at colliders, but a dedicated analysis of LHC data was proposed to search for light left-handed sneutrinos in the framework of the µνSSM [30]. However, the analysis in Ref. [30] concentrated on τ -sneutrinos as the LSP, whereas here there are even lighter SUSY particles in the spectrum. For a detailed discussion of distinct signatures at the LHC related to left-handed sneutrinos within the µνSSM we refer to the literature [29][30][31].

Scan over λ
As already explained at the beginning of this section, there is no theoretical reason to choose the µνSSM-like parameters universal w.r.t. the family index. While the tree-level upper bound on the lightest CP-even scalar mass approximately depends on the term λ λ λ 2 = λ i λ i , and not on the individual λ i , this is not the case for the SM-like Higgs-boson mass, as soon as mixing-effects induced Tree-level One-loop Two-loop Figure 7: Light CP-even scalar spectrum in the scan over v R . Shown are the masses at tree level (dotted), at the one-loop level (dashed) and at the partial two-loop level (solid). We show in the brackets the dominant composition of the loop-corrected mass eigenstates h (2) in the experimentally allowed region of v R . The desired SM-like Higgs-boson mass is indicated with the horizontal green band, assuming a theory uncertainty of 3 GeV. The red regions are excluded by direct searches for additional scalars. In the gray region the SM-like Higgs-boson mass is not predicted accurately. On top we show χ 2 red for various Higgs-boson signal strength measurements at LHC.  Table 4: The same as in Tab. 2 for the scan over λ i while λ λ λ 2 = const.
by the right-handed sneutrinos are considered. This effects can enter at tree level, or via radiative corrections proportional to λ i . These radiative corrections depend on the masses and the mixing of each of the right-handed sneutrino. Since it is not the case that all three ν iR are degenerate, the radiative corrections are expected to depend strongly on the individual values of λ i . Also, when λ λ λ 2 is fixed, the µ-term which is dynamically generated after EWSB and linearly dependent on λ i cannot be constant when the λ i are varied. This can be another source of corrections to the SM-like Higgs-boson mass that explicitly depend on the individual values of the λ i . However, the loop corrections proportional to λ i are an order of magnitude smaller than the ones stemming from the (s)top-sector, partially because quantum contributions to the SM-like Higgs-boson mass at the one-loop level proportional to λ i depend on the singlet-admixture of the SM-like Higgs boson which, in turn, cannot be too large to not spoil the measured signal strengths at the LHC. Nevertheless, we will give here a rough idea of how large the remnant effect of nonuniversal λ i on the SM-like Higgs-boson mass can be while λ λ λ 2 is kept constant. We performed a parameter scan over all possible values of the λ i in a scenario in which ν R µR has a mass between 92 and 115 GeV and mixes substantially with the SM-like Higgs boson. The free parameters were set to the values shown in Tab. 4.
On can see that the scenario is very similar to the one in Sect. 5.1. When the λ i are chosen uniformly we recover the BP at λ = 0.168 in Fig. 4 which lies in the middle of the experimentally allowed parameter region. In Fig. 8 we illustrate the dependence of the SM-like Higgs-boson mass on the individual values of λ i . We show triangle plots [156] with the values of λ 2 i on the axes and with their sum λ λ λ 2 = const. The colors of the points indicate the mass of the SM-like Higgs boson at tree level (top left), at one-loop level (top right), the difference of the SM-like Higgs-boson mass at tree level and one-loop level (bottom left), and the one-loop mass of the right-handed µ-sneutrino (bottom right) which is the lightest CP-even scalar in this scenario. We do not show the two-loop mass of the SM-like Higgs boson, because the two-loop corrections supplemented from FeynHiggs are purely MSSM-like corrections independent of λ i , thus not playing a role in the following discussion. However, the parameters are chosen such that the corrections beyond the one-loop level shift the SM-like Higgs-boson mass into the vicinity of ∼ 125 GeV (see below). In the upper right plot one sees that the one-loop mass of the SM-like Higgs boson is the largest in the central point in which all λ i are equal. The tree-level mass, on the other hand, shows the opposite behavior and is the largest in the corners of the upper left plot in which one of the λ i is practically zero. The one-loop mass varies in the experimentally allowed region by more than 1 GeV. This demonstrates that for an accurate prediction of the SM-like Higgs-boson mass it is crucial to include the independent contributions of all three λ i to the radiative corrections, when mixing effects between the right-handed sneutrinos and the SM-like Higgs boson are sizable.
Note that the variation of the one-loop mass would be even larger if we neglect the experimental h SM (bottom left), and the one-loop mass of ν µR (bottom right) for fixed λ λ λ 2 but varying λ i , indicated by the colors (and sizes for better visibility). Crossed points are excluded by either HiggsBounds, HiggsSignals with χ 2 red ≥ 150/101 = 1.485 for the 13 TeV data set, or because the SM-like Higgs boson mass including the two-loop corrections is smaller than 123 GeV or larger than 127 GeV.
constrains. In this scenario, the main exclusion limit is the requirement to have the SM-like Higgsboson mass above 123 GeV at two-loop. In the corners of the plots, the mass of ν R µR , shown in the lower right plot, increases to values very close to the SM-like Higgs-boson mass. This increases the mixing between both scalars which, in turn, reduces the radiative one-loop corrections to the SM-like Higgs boson. Practically speaking, parts of the loop-corrections "are lost" to ν R µR . This is why the difference between the tree-level and the one-loop mass of the SM-like Higgs boson, shown in the lower left plot of Fig. 8, is the smallest when the mass of ν R µR is the largest. We also emphasize that the variation of the difference between tree-level and one-loop mass, as a measure for the size of genuine one-loop corrections, is more than twice as large as the variation of the one-loop mass of the SM-like Higgs boson. One can see a compensation of lower tree-level mass, but larger loop corrections in the center of the plots, leading to a more stable one-loop mass of the SM-like Higgs boson. This is due to the fact that the loop-corrected mass eigenstate of the SM-like Higgs boson is the physical state, whose properties are constrained by the experimental measurements. Therefore, the singlet-admixture cannot be very large at loop-level, and the parameter dependence of the SM-like Higgs-boson mass induced by this admixture necessarily cannot be too large. The tree-level state, on the other hand, is not physical and can have larger mixing with ν R µR , leading to a stronger dependence on parameters related to the sneutrinos, like λ i in this case. Thus, even though the overall dependence of the SM-like Higgs-boson mass on parameters beyond the MSSM can be very large at tree level, the dependence will be diminished at loop-level in the parameter space that fulfills experimental constraints on the SM-like Higgs-boson properties.

CMS and LEP excess at ∼ 96 GeV
Searches for the SM Higgs boson at LEP can nowadays be used to constrain the parameter space of models with additional scalar particles with masses below the SM-like Higgs-boson mass. Interestingly, an excess over the SM background was observed at a mass of 95-98 GeV in the Higgsstrahlung production channel with an associated decay of the Higgs boson to a pair of b-quarks [154]. Assuming that this excess can be explained by an additional Higgs boson, the signal strength, i.e, the cross section times branching ratio to b-quarks normalized to the SM prediction for a SM-like Higgs boson at the same mass, was extracted in Ref. [157] to be µ LEP = 0.117 ± 0.057 . (5.4) At about the same mass, CMS observed an excess over the SM background in the pp production with associated decay to diphotons in the 8 TeV and the 13 TeV data [158,159]. For the CMS excess, the signal strength was reported in Refs. [159,160] to be In a previous publication, we showed that it is possible to accommodate both excesses simultaneously at the 1σ level in the µνSSM with just one right-handed neutrino [99]. We described how a right-handed sneutrino at ∼ 96 GeV can acquire substantial couplings to SM particles via its mixing with the SM-like Higgs boson. 7 However, we did not include an accurate prediction of the properties of the neutrinos in our analysis in the µνSSM with one generation of right-handed neutrinos, because in that case at least one neutrino mass has to be generated via quantum corrections. On the contrary, in the µνSSM with three generations of right-handed neutrinos, we can describe the mass differences and mixings of the neutrinos at tree level, which is of course much more feasible.  Table 5: Parameters of the scan to fit the LEP and the CMS excesses. Dimensionful parameters are given in GeV. If the family index is omitted the parameter has a universal value independent of the index.
The values of the free parameters to fit the excesses are shown in Tab. 5. The scalar used to fit the excesses is the right-handed τ -sneutrino ν R τ R . This is assured by setting κ 333 κ 111,222 , so that ν R τ R has a smaller mass than the other two right-handed sneutrinos. Then, v 3R and A κ 333 are used to tune the mass of ν R τ R to be at ∼ 95-98 GeV. A sufficiently large mixing of ν R τ R with the SM-like Higgs boson is achieved with a large value of λ 3 ∼ 0.54, while λ 1,2 = 0.01 are very small to avoid that the effective µ-term becomes very large. Alternatively, one could have used smaller values for v 1,2R , but then the other two CP-even sneutrinos ν R e,µ R would have been very light as well, potentially carrying away some of the mixing between the SM-like Higgs boson and the righthanded sneutrinos. 8 Since λ 3 is large, the SM-like Higgs boson receives additional contribution to the tree-level mass. This is why tan β is set to a small value, and, besides A u 3 = −650 GeV, the soft trilinear parameters A u,d,e can be set to zero. Apart from that, tan β shall not be much larger than one to not suppress the coupling of ν R τ R to t-quarks, which scales with the inverse of sin β (see Eq. (5.3)). Finally, we choose a range for A λ 3 in which the mixing between ν R τ R and the SM-like Higgs boson is sufficiently large. We will show results for small ranges of the parameters λ 3 , v 3R , A λ 3 and A κ 333 . While v 3R and A κ 333 are mainly correlated to the mass of ν R τ R , λ 3 and A λ 3 affect the doublet composition of ν 3R . This certainly is not an exhaustive parameter scan covering the complete parameter space, but the scan gives an idea of how the excesses can be accommodated within the µνSSM, and it resembles the solution we found in the µνSSM with just one right-handed neutrino [99].
In Fig. 9 we show the results for the signal strength of the LEP excess µ LEP (top) and of the CMS excess µ CMS (bottom). In both plots the colors of the points indicate the SM-like Higgs-boson mass, while the mass of ν R τ R is shown on the horizontal axis. The signal strengths were calculated in the narrow-width approximation, and the branching ratio and cross section ratios w.r.t. the SM where calculated using the effective coupling approximation as explained in Ref. [99]. One can immediately see that it is rather easy to achieve the experimental value of µ LEP , whereas the largest values for µ CMS reached in our scan are just below the lower limit of µ CMS = 0.4, i.e., 1σ below the central value. This is due to the fact that the main decay channel of ν R τ R is the decay to a pair of bottom quarks, and it is harder to achieve a substantial branching ratio to diphotons required for the CMS excess. Nevertheless, both excesses are fitted at the 1σ level considering the experimental uncertainties, while fitting the neutrino data and being in agreement with the experimental constraints on the SM-like Higgs boson, which we again checked with HiggsSignals assuming a theoretical uncertainty of the SM-like Higgs-boson mass of 3 GeV.
In Fig. 10 we show the correlation of both signal strengths, with the colors encoding the mass of ν R τ R . The strong correlation one can see has its origin in the fact that both signal strengths increase with the amount of doublet-component of ν R τ R . In principle, one could achieve a further enhancement of µ CMS and a suppression of µ LEP by suppressing the down-type doublet component of ν R τ R . Then, the branching ratio to bottom quarks becomes smaller, and the diphoton branching ratio increases because of the smaller total decay width due to the reduction of the decay width to bottom quarks. However, finding such points is difficult, because the dominant terms mixing the right-handed sneutrinos ν R iR with the doublet fields H R d and H R u scale equally with λ i , A λ i , κ ijk and v iR at tree level, as can be seen in Eqs. (A.4) and (A.5). The only difference are the factors v d and v u in each equation, respectively, which cannot be exploited too much, because, as mentioned before, tan β should not be too far from one. This is why an extensive scan of the vast parameter space of the µνSSM would be necessary to find parameter points in which µ CMS is further enhanced without increasing µ LEP too much, which however lies beyond the scope of this paper.
Instead, we will focus on the rest of the spectrum, which heavily depends on the values of the neutrino Yukawa couplings Y ν ij and the vevs of the left-handed sneutrinos v iL , once the remaining parameters are fixed to the values listed in Tab. 5. As an example, we show in Tab. 6 two possible sets of parameters that accommodate accurate neutrino masses and mixings in the parameter scan of this section. In contrast to the other scenarios we presented before, here we will make use of non-zero values of the non-diagonal elements of Y ν ij in one of the BPs. Naturally, this simplifies the accommodation of neutrino properties in agreement with experimental data, because there are six more free parameters that can be adjusted. The price to pay is that there is usually more than one set of parameters of Y ν ij and v iL that give accurate predictions for the neutrino sector. The reason for showing two distinct scenarios is that in the µνSSM the scalar sector is deeply related to the neutrino sector. Thus, different sets of parameters predict fundamentally different scalar spectra, and since there is no theoretical argument that the neutrino Yukawa couplings have to be diagonal, we used the additional freedom to present a point in which, on top of the explanation of the LEP 8.109 · 10 −8 0 0 0 1.154 · 10 −7 0 BP2 7.088 · 10 −8 1.181 · 10 −8 −3.404 · 10 −9 1.902 · 10 −8 1.238 · 10 −7 1.783 BP1 0 0 8.855 · 10 −7 1.890 · 10 −5 2.601 · 10 −4 1.871 · 10 −4 BP2 −2.103 · 10 −9 6.923 · 10 −9 1.383 · 10 −8 1.792 · 10 −6 2.072 · 10 −4 3.673   In both BPs the lightest BSM particle is the right-handed τ -neutrino. This is because κ 333 has to be small to decrease the mass of the corresponding sneutrino ν R τ R . Consequently, also the Majorana mass term for the neutrino will be small and m ν τ R is small. However, its mass is still above half the SM-like Higgs-boson mass, so the decay of the SM-like Higgs boson into ν τ R is forbidden. The striking difference between both BPs is the mass scale of the left-handed τ -sneutrinos and sleptons. In BP2 the Yukawa coupling Y ν 33 is the smallest diagonal element of Y ν , and the vev corresponding to the third family of left-handed sneutrinos v 3L is the largest of the three. This reduces the masses of the left-handed CP-even and CP-odd τ -sneutrino (see Eq. (2.23)) and the τ -slepton to values below the SM-like Higgs-boson mass. In BP1, on the other hand, v 2L is the largest of the three left-handed vevs, and therefore the left-handed µ-sneutrinos and µ-slepton are the lightest left-handed sfermions, although still more than twice as heavy as the SM-like Higgs boson. In this way, the phenomenology of both BPs is distinct, even though the properties of the right-handed τ -sneutrino, which is the particle used to fit the excesses, are not affected by the left-handed sector. This is because its branching ratios are dominantly given by the mixing-effects with the SM-like Higgs boson, which is not suppressed by the small neutrino Yukawa couplings Y ν ij . Nonetheless, BP2 can give rise to additional interesting signal at colliders. A dedicated analysis of the collider phenomenology of light left-handed τ -sneutrinos/sleptons at the LHC can be found in Refs. [29,30], where it was shown that there are no direct bounds from LEP/LHC searches that can be used to set lower limits on the masses of these particles in the framework of the µνSSM. This analysis made use of the fact that the charged τ L can be produced, and thus provide a source for the left-handed sneutrinos, since it decays into lighter ν R τ L or ν I τ L . An important feature is that the subsequent decays of ν R,I τ L can be prompt or displaced. However, in Refs. [29,30] it is assumed that the left-handed sneutrino is the LSP. This is not the case here, since the right-handed τ -neutrino is even lighter. Signals at colliders from the lightest BSM particle ν τ R are not expected, because it is a gauge singlet, thus cannot be produced directly. In principle, it can be produced indirectly via the decay of the sfermions. However, the spectrum is very compressed, such that a pair production of ν τ R from the decays of ν τ R at ∼ 96 GeV or the SM-like Higgs boson at ∼ 125 GeV is kinematically forbidden, and the production of a pair of a right-and a left-handed τ -neutrino is suppressed by the size of Y ν 33 .

Conclusion and Outlook
The µνSSM is a simple SUSY extension of the SM that is capable of describing neutrino physics in agreement with experimental data. As in other SUSY models, higher-order corrections are crucial to reach a theoretical uncertainty at the same level of the (anticipated) experimental accuracy. So far, higher-order corrections in the µνSSM had been restricted to DR calculations, which suffer from the disadvantage that they cannot be directly connected to (possibly future observed) new BSM particles. More recently we had evaluated the corrections to the neutral scalar masses and mixings, but restricting ourselves to one generation of heavy neutrinos [99].
In this paper we have presented the complete one-loop renormalization of the neutral scalar sector of the µνSSM with three generation of right-handed neutrinos in a mixed on-shell/DR scheme. In this way, for the first time, it is possible to evaluate the masses and mixings in the neutral scalar sector with high precision, while simultaneously describe correctly the experimental neutrino data, such as mass differences and mixing angles. An on-shell (OS) renormalization has been chosen for parameters that can be directly identified with (potentially) observable states, whereas a DR renormalization has been chosen for the remaining parameters. We provide details on the full set renormalization parameters, which were implemented into a FeynArts model file (that can be provided by the authors upon request).
We have performed the calculation of the masses of the neutral scalars in the µνSSM with three generations of heavy neutrinos, taking into account the full set of one-loop corrections as described above. These corrections have been supplemented by contributions to the neutral Higgs-boson sector of the MSSM at and beyond the two-loop level as provided by FeynHiggs. These corrections are crucial to obtain a reliable prediction of the mass of an SM-like Higgs boson around ∼ 125 GeV.
The masses of the neutral scalar bosons have been evaluated in a set of scenarios that exemplify the relevant dependences on the underlying µνSSM parameters. The scenarios are in agreement with all available searches for additional Higgs bosons (via the code HiggsBounds), as well with the properties of the SM-like Higgs boson as measured at the LHC (via the code HiggsSignals), while at the same time reproducing correctly the measured values of neutrino mass differences and mixing angles. In a first scenario we varied the assumed to be universal parameter λ. We find that large one-loop corrections for the right-handed sneutrinos arise, that are crucial to accurately account for possible mixing effects between them and the SM-like Higgs boson. In a second scenario we have varied the assumed to be universal vev of the right-handed sneutrinos. Here we find that all three CP-even right-handed sneutrinos can have masses below ∼ 125 GeV, without being excluded by cross-section limits from direct searches for additional Higgs bosons. In our third scenario we deviate from the intergenerational universality assumptions. We take λ λ λ 2 := λ i λ i to be constant, but vary instead the individual λ i . We find that the non-universality of the λ i has an important impact on the predictions of the neutral scalar masses. It has been shown that in the experimentally allowed parameter space the non-universality of the λ i can account for deviations of the SM-like Higgs-boson mass of ∼ 1 GeV, emphasizing the importance to consider the full set of µνSSM parameters at the one-loop level.
As a final example we have discussed how the µνSSM can describe two excesses in the searches for light Higgs bosons in the vicinity of ∼ 96 GeV. These are a ∼ 3 σ excess in the diphoton final state as reported by CMS and a ∼ 2 σ excess in the bb final state as published by LEP. We demonstrated that the µνSSM can account at the ∼ 1 σ level for both excesses, while being in agreement with all available Higgs-boson searches and measurements, as well as the available neutrino data. We are eagerly awaiting updated experimental analyses from ATLAS and CMS to confirm or refute these excesses.
Further explorations of the scalar sector of the µνSSM are necessary to cover the wealth of (possible) phenomenology that this model offers. This includes further studies dropping the (artificial) intergenerational universality assumptions, where in this paper we have taken only the first step. It furthermore includes the evaluation of the charged scalar sector at the one-loop level and beyond, which is expected to be important for the phenomenology at the LHC. Also studies going beyond the LHC searches may become relevant, such as analyses of the possibility of a first order phase transition, leading to gravitational waves created in the early universe. We leave these studies for future work.
where the indices run from 1 to 3. We checked that the coefficients of the divergent part of the field renormalization counterterms are equal to the one-loop anomalous dimensions of the corresponding superfields γ (1) ij , neglecting the terms proportional to the gauge couplings g 1 and g 2 , and divided by the loop factor 16π 2 , i.e., which is the same relation that holds in the (N)MSSM.

B.2.1 Renormalization conditions
In the following we list the renormalization conditions for the counterterms of the parameter used to renormalize the neutral scalar potential:  where δX stands for one of the counterterms just mentioned, and β (1) X is the one-loop coefficient of the beta function of the parameter X, which could be obtained by the help of the mathematica package SARAH [169].
On the contrary, the counterterms of the soft masses stated in Eqs. (B.40)-(B.42) are the ones derived from the one-loop beta function we obtained with SARAH, which were then numerically checked to be equal to the counterterms for (m 2   Table 7: Values for parameters of the standard model in GeV.
The value for v corresponds to a value for the Fermi constant of G F = 1.166 38 × 10 −5 GeV −2 . The values for the gauge boson masses define the cosine of the weak mixing angle to be c w = 0.881 535. Note that since the SM leptons mix with the Higgsinos and gauginos in the µνSSM, the lepton masses are not the real phyiscal input parameters. However, the mixing is tiny, so there will always be three mass eigenstates in the charged fermion sector corresponding to the three standard model leptons, having approximately the masses m e , m µ and m τ . This is why we use the measured values for these masses, such as m τ in Tab. 7, and then calculate the real input parameters, which are the Yukawa couplings