Research Article  Open Access
The Higgs Sector of the Minimal SUSY Model
Abstract
I review the Higgs sector of the extension of the minimal supersymmetric standard model (MSSM). I will show that the gauge kinetic mixing plays a crucial role in the Higgs phenomenology. Two light bosons are present, a MSSMlike one and a like one, which mix at one loop solely due to the gauge mixing. After briefly looking at constraints from flavour observables, new decay channels involving righthanded (s)neutrinos are presented. Finally, how model features pertaining to the gauge extension affect the model phenomenology, concerning the existence of RParityconserving minima at loop level and the Higgstodiphoton coupling, will be reviewed.
1. Introduction
The recently discovered Higgs boson is considered as the last missing piece of the standard model (SM) of particle physics. Nonetheless, several firm observations univocally call for its extension, mainly, but not limited to, the presence of dark matter, the neutrino masses and mixing pattern, the stability of the SM vacuum, and the hierarchy problem. Supersymmetry (SUSY) has long been considered as the most appealing framework to extend the SM. Its minimal realisations (MSSM and its constrained versions (for a review, see [1])) start however to feel considerable pressure to accommodate the recent findings, especially the measured Higgs mass of GeV. Despite not in open contrast with the MSSM, the degree of fine tuning required to achieve it is more and more felt as unnatural. In order to alleviate this tension, nonminimal SUSY realisations can be considered. One can either extend the MSSM by the inclusion of extra singlets (e.g., NMSSM [2]) or by extending its gauge group. Concerning the latter, one of the simplest possibilities is to add an additional Abelian gauge group. I will focus here on the presence of group which can be a result of an heterotic string theory (and hence Mtheory) [3–5]. This model, the minimal Parityconserving supersymmetric standard model (BLSSM in short), was proposed in [6, 7] and neutrino masses are obtained via a type I seesaw mechanism. Furthermore, it could help to understand the origin of Parity and its possible spontaneous violation in supersymmetric models [6–8] as well as the mechanism of leptogenesis [9, 10].
It was early pointed out that the presence of two Abelian gauge groups in this model gives rise to kinetic mixing terms of the formwhich are allowed by gauge and Lorentz invariance [11], as and are gaugeinvariant quantities by themselves; see, for example, [12]. Even if these terms are absent at tree level at a particular scale, they will in general be generated by RGE effects [13, 14]. These terms can have a sizable effect on the mass spectrum of this model, as studied in detail in [15], and on the dark matter, where several scenarios would not work if kinetic mixing is neglected, as thoroughly investigated in [16]. In this work, I will review the properties of the Higgs sector of the model. Two light states exist, a MSSMlike boson and a like boson. After reviewing the model, I will show that a large portion of parameter space exists where the SMlike Higgs boson has a mass compatible with its measure, both in a “normal” ( GeV) and in an “inverted” hierarchy ( GeV), also in agreement with bounds from low energy observables and dark matter relic abundance. The phenomenological properties of the two lightest Higgs bosons will be systematically investigated, where once again the gauge mixing is shown to be fundamental. The presence of extra Dterms arising from the new sector, as compared to models based on the SM gauge symmetry, has a large impact on the model phenomenology. They affect both the vacuum structure of the model and the Higgs sector, in particular enhancing the Higgstodiphoton coupling. Both of these issues will be reviewed here, although the latter is disfavoured by recent data [17], to show model features beyond the MSSM.
2. The Model
For a detailed discussion of the masses of all particles as well as of the corresponding oneloop corrections, we refer to [15]. Attention will be paid to the main aspects of the kinetic mixing since it has important consequence for the scalar sector. For the numerical investigations that will be shown, we used the SPheno version [22, 23] created with SARAH [24–28] for the BLSSM. For the standardised model definitions, see [29], while for a review of the model implementation in SARAH, see [30]. This spectrum calculator performs a twoloop RGE evaluation and calculates the mass spectrum at one loop. In addition, it calculates the decay widths and branching ratios (BRs) of all SUSY and Higgs particles as well as low energy observables like . We will discuss the most constrained scenario with a universal scalar mass , a universal gaugino mass , and trilinear softbreaking couplings proportional to the superpotential coupling () at the GUT scale. Other input parameters are , , , , and . They will be defined in the following section. The numerical study here presented has been performed by randomly scanning over the independent input parameters above described via the SSP toolbox [31], while low energy observables such as BR() and BR() have been evaluated with the FlavourKit package [32]. Furthermore, during the scans, all points have been checked with HiggsBounds4.1.1 [33–36], both in the “normal” hierarchy and in the “inverted” hierarchy case.
2.1. Particle Content and Superpotential
The model consists of three generations of matter particles including righthanded neutrinos which can, for example, be embedded in 16plets. Moreover, below the GUT scale, the usual MSSM Higgs doublets are present as well as two fields and responsible for the breaking of . The field is also responsible for generating a Majorana mass term for the righthanded neutrinos and thus we interpret its charge as its lepton number. The same goes for , and we call these fields bileptons since they carry twice the lepton number of (anti)neutrinos. The quantum numbers of the chiral superfields with respect to are summarised in Table 1.

The superpotential is given byand we have the additional soft SUSYbreaking terms:where , are generation indices. Without loss of generality, one can take and to be real. The extended gauge group breaks to as the Higgs fields and bileptons receive vacuum expectation values ():We define in analogy to the ratio of MSSM ().
2.2. Gauge Kinetic Mixing
As already mentioned in the Introduction, the presence of two Abelian gauge groups in combination with the given particle content gives rise to a new effect absent in any model with just one Abelian gauge group: gauge kinetic mixing. This can be seen most easily by inspecting the matrix of the anomalous dimension, which for our model at one loop readswith typical GUT normalisation of the two Abelian gauge groups, that is, for and for [7]. Therefore, even if at the GUT scale the kinetic mixing terms are zero, they are induced via RGE evaluation at lower scales. It turns out that it is more convenient to work with noncanonical covariant derivatives rather than with offdiagonal fieldstrength tensors as in (1). However, both approaches are equivalent [37]. Therefore, in the following, we consider covariant derivatives of the form , where is a vector containing the charges of the field with respect to the two Abelian gauge groups, is the gauge coupling matrixand contains the gauge bosons .
As long as the two Abelian gauge groups are unbroken, we have still the freedom to perform a change of basis by means of a suitable rotation. A convenient choice is the basis where , since in this case only the Higgs doublets contribute to the gauge boson mass matrix of the sector, while the impact of and is only in the offdiagonal elements. Therefore, we choose the following basis at the electroweak scale [38]:
When unification at some large scale ( GeV) is imposed, that is, and , at SUSY scale, we get [15]
2.3. Tadpole Equations
The minimisation of the scalar potential is here described in the socalled tadpole method. We can solve the treelevel tadpole equations arising from the minimum conditions of the vacuum with respect to , , , and . Using and , we obtain, and, thus, we find an approximate relation between and For the numerical results, the oneloop corrected equations are used, which lead to a shift of the solutions in (9)–(12).
2.4. The Scalar Sector
In this model, MSSM complex doublets and bilepton complex singlets are present, yielding even, odd, and charged physical scalars.
Concerning the even scalars, the MSSM and bilepton sectors are almost decoupled, mixing exclusively due to the gauge kinetic mixing. In first approximation, the mass matrix is blockdiagonal and has mass eigenstates that mimic the MSSM case. In practice, it turns out that only two Higgs bosons are light (hereafter called and , one per sector), while the other two are very heavy (above the TeV scale). The lightest scalars are well defined states, being either almost exclusively doubletlike or bileptonlike. It is worth stressing that their mixing is small (see Figure 4) and solely due to the gauge kinetic mixing (see also [39]).
Concerning the physical pseudoscalars and , their masses are given byFor completeness, we note that the mass of charged Higgs boson reads as in the MSSM as
In this model, the odd and charged Higgses are typically very heavy. In (10), we see that, compared to the MSSM, there is a nonnegligible contribution from the gauge kinetic mixing. LHC searches limit and TeV, since [40, 41]at C.L. Notice that recent reanalysis of LEP precision data also constrains TeV at C.L. [42]. A consequence of this strong constraint in the BLSSM is that the first terms in (10)–(12) can be large, pushing for odd and charged Higgs masses in the TeV range.
The very large bound on the mass is in contrast with the nonSUSY version of the model, where the gauge couplings are free parameters and can be much smaller, hence yielding lower mass bounds. The latter need to be evaluated as a function of both gauge couplings [43].
Next, we describe the sneutrino sector, which shows two distinct features compared to the MSSM. Firstly, it gets enlarged by the superpartners of the righthanded neutrinos. Secondly, even more drastically, a splitting between the real and imaginary parts of each sneutrino occurs resulting in twelve states: six scalar sneutrinos and six pseudoscalar ones [44, 45]. The origin of this splitting is the term in the superpotential (see (2)), which is operator after the breaking of . In the case of complex trilinear couplings or terms, a mixing between the scalar and pseudoscalar particles occurs, resulting in 12 mixed states and consequently in a mass matrix.
To gain some feeling for the behaviour of the sneutrino masses, we can consider a simplified setup: neglecting kinetic mixing as well as leftright mixing, the masses of the Rsneutrinos at the SUSY scale can be expressed asIn addition, we treat the parameters , , , , , and as independent. The different effects on the sneutrino masses can easily be understood by inspecting (17). The first two terms give always a positive contribution whereas the third one gives a contribution that can be potentially large which differs in sign between the scalar and pseudoscalar states, therefore inducing a large mass splitting between the states. Further, this contribution can either be positive or negative depending on the sign of . For example, choosing and positive, one finds that the even (odd) sneutrino is the lightest one for (). This is pictorially shown in Figure 1, as a function of the GUTscale input parameter , for a choice of the other parameters. One notices that the even (odd) sneutrino is the lightest one when the GeV Higgs boson is predominantly (). It is worth pointing out here that, as will be described in the following section, when GeV, the nexttolightest Higgs boson can decay into pairs of even sneutrinos, but not into similar channel with odd sneutrinos. Being predominantly a bilepton field, when this decay is open, it saturates its BRs; see Figure 3. Regarding the decay into odd sneutrinos, this channel is accessible (i.e., is light enough) only in the region where is the SMlike Higgs boson, that is, mainly coming from the doublets. In this case, however, this decay channel is mitigated by the small scalar mixing and is not overwhelming (unlike for , now mainly from the bileptons).
Depending on the parameters, either type of sneutrinos can get very light. For the LSP, it can be a suitable dark matter candidate [16] and yield extra fully invisible decay channels to the Higgs bosons, thereby increasing their invisible widths. In the case of the decay into the odd sneutrino, since this can happen mainly for the SMlike Higgs boson, one should account for the constraints on the former [17]. Eventually, the Rsneutrinos could also get tachyonic or develop dangerous Parityviolating . While the first possibility is taken into account in our numerical evaluation by SPheno, and such points are excluded from our scans, the second case will be reviewed in the following subsection.
The last important sector for considerations that will follow is the one of the charged sleptons. See [46] for further details. New SUSY breaking Dterm contributions to the masses appear, which can be parametrised as a function of the mass and of asTheir impact is larger for the sleptons than for the squarks by a factor of due to the different charges (). It is possible to vary the stau mass by GeV with respect to the MSSM case while keeping the impact on the squarks under control. Having different sfermion masses in the BLSSM as compared to the MSSM has a net impact on the Higgs phenomenology, in particular in enhancing the coupling while keeping unaltered the SMlike Higgs coupling to gluons. As described at the end of this review, the new Dterms coming from the sector can further reduce the stau mass entering in the effective interaction (while ensuring a pole mass of GeV, compatible with exclusions) (with pole mass we denote the oneloop corrected mass at , while in the loop, leading to the effective coupling, the running treelevel mass at enters, being the SMlike Higgs boson; i.e., GeV) leading this mechanism to work also in the constrained version of the model. This mechanism has been recently reanalysed also in [47] in the very same model.
2.5. The Issue of RParity Conservation
We have encountered so far several neutral scalar fields which could develop , beside the Higgs bosons. If of fields charged under QCD and electromagnetism are forbidden because the latter are good symmetries, Rsneutrino , which are not by themselves problematic, would unavoidably break RParity. The issue of conserving RParity is of fundamental importance, since this is a builtin symmetry in our model where is gauged. We will therefore restrain ourselves to parameter configurations where the global minimum is RParity conserving.
When all neutral scalar fields are allowed to get , it is not trivial even at the tree level to find which is the deeper global minimum and whether it is of a “good” type, here defined as having the correct broken symmetries and being RParity conserving. One possible way to study this issue is to start from a simplified set of input parameters yielding a correct treelevel global minimum when only the Higgs fields get . and then look for the true global minimum when all other neutral fields (mainly Rsneutrinos) acquire , both at the tree level and at loop level. See [48] for further details.
At the tree level there seems to exist regions where the BLSSM has a stable, RParityconserving global minimum with the correct broken and unbroken gauge groups. For this to happen one needs the Rsneutrino Yukawa coupling to be not so large and the trilinear parameter to be not large compared to the soft scalar mass , as, intuitively, large and can lead to large negative contributions to the potential energy for large values of , as well as reducing the effective Rsneutrino masses, as described above and clear from Figure 1.
It turns out that when loop corrections are taken into account, few points all over such regions of parameters exist where RParity is not preserved anymore, or where or is unbroken. This is apparently due to a very finely tuned breaking of and which often does not survive loop corrections. The reason for this is that, besides the known large contributions of third generation (s)fermions, the additional new particles of the sector also play an important role. As previously described for the charged sleptons sector, new SUSY breaking Dterm contributions to the masses appear; see (18). Since, as shown in (16), the experimental bounds require to be in the multiTeV range, these contributions can be much larger than in the MSSM sector, resulting in the observed importance of the corresponding loop contributions. Furthermore, these contributions are also responsible for the restoration of at the oneloop level.
Ultimately, overall safe regions of parameters cannot be found where the correct vacuum structure can be ensured. At the same time, if naive trends can be spotted for bad points to appear, these have nonetheless to be checked case by case due to the highly nontrivial scalar potential, and it might be possible that neighbour configurations still hold a valid global minimum. We will not check the validity of our scans from the vacuum point of view in the following, being confident that if any point is ruled out, a neighbour one yielding a very similar phenomenology can be found, which is allowed.
3. A Quick Look at Flavour Observables
Before moving to the Higgs phenomenology, we briefly show the impact on the BLSSM model when considering the constraints arising from low energy observables. For a review of the observables as well as for the impact on general SUSY models encompassing a seesaw mechanism, see [49, 50].
We consider here only the two most constraining ones, BR() and BR(). The present exclusions are BR() [18] and BR() [19]. In Figure 2 we plot these branching ratios as a function of the mass of the lightest (in black) and nexttolightest (in red) SMlike neutrino, which display some pattern for evading the bounds. In particular, they are required to be rather light, below eV, while the model, due to the limited scans here performed, seems to prefer configurations with neutrinos heavier than 0.01 eV, hence the preferred region in between. Lighter mass values are nonetheless also allowed.
(a)
(b)
For convenience, the impact of satisfying the earlier bounds will be shown only in the inverted hierarchy case, due to the smaller density of configurations therein. Instead, points not allowed in the normal hierarchy case are automatically dropped.
Regarding the longlasting discrepancy, in the setup investigated here charginos and charged Higgses are too heavy, same for the boson, while the neutralino and sneutrino are too weakly coupled, to give a significant enhancement over the SM prediction.
4. Higgs Phenomenology
We review here the phenomenology of the Higgs sector, showing a first survey of its phenomenological features. First, results when normal hierarchy is imposed are presented. Then, we will show that the inverted hierarchy is also possible on a large portion of the parameter space. Without aim for completeness, the results are here presented as the starting point for a more thorough investigation. Finally, how model features pertaining to the extended gauge sector impinge on the Higgs phenomenology and in particular how the Higgstodiphoton branching ratio can be easily enhanced in this model, despite the experimental data now converging to a more SMlike behaviour than in the recent past, are described.
4.1. Normal Hierarchy
In this subsection we discuss the normal hierarchy case, with the lightest Higgs boson being the SMlike one (i.e., predominantly from the doublets), and a heavier Higgs boson predominantly from the bilepton fields (those carrying number and responsible for the spontaneous breaking). Their mixing is going to be small and solely due to the kinetic mixing.
In Figure 3 we first inspect the heavy Higgs boson branching ratios. Besides the standard decay modes, the decay into a pair of SM Higgs bosons exists, as well as two new characteristic channels of this model, comprising righthanded (s)neutrinos:(1). Its BR can be up to before the top quark threshold and around afterwards.(2). A similar decay channel exists for the boson. The BRs are , up to depending on the heavy Higgs and neutrino masses.(3), where, is the even sneutrino and the LSP, hence providing fully invisible decays of the heavy Higgs. If kinematically open, it saturates the Higgs BRs. Notice that only points with very light even sneutrinos are shown, possible only for very large and negative (see Figure 1).
While the first two channels exist also in the nonSUSY version of the model (however, in the nonSUSY model, the Higgs mixing angle is a free parameter, directly impacting these branching ratios) (see, e.g., [51]), the last one, involving the even sneutrino, is truly new and rather intriguing. This is because the sneutrino is light and it can be a viable LSP candidate if its mass is smaller than , as in this case [16]. It however implies that the heavy Higgs is predominantly bileptonlike, with a light Higgs very much SMlike. This can be seen in Figure 4, where the points with large BR() (in red) have the lowest mixing between and the SM scalar doublet fields, of the order of . It immediately follows that this channel will have very small cross section at the LHC, when considering SMlike Higgs production mechanisms. This is true for all heavy Higgs masses GeV. The GeV Higgs is well SMlike, with tiny reduction of its couplings to the SM particle content. On the other hand, the heavy Higgs is feebly mixed with the doublets, suppressing its interactions with the SM particles and hence its production cross section. This can be seen in Figure 5(a). Considering only the gluonfusion production mechanism, and multiplying it by the relevant BR, we get the cross sections for the choice of channels displayed therein. The most constraining channels, and , are also compared to the exclusions at the LHC for TeV from [20] and [21], respectively. The channels are well below current exclusions, which are hence not shown.
(a)
(b)
We see that all (starting from GeV) the displayed configurations are allowed by the current searches (the exclusions shown by solid curves of the same color as the depicted channel). This is because of the suppression of the heavy Higgs boson cross sections due to the small scalar mixing.
In the lower plot the cross sections for the new channels are displayed. Those pertaining to model configurations for which the heavy Higgs boson decays to the even sneutrino (LSP), yielding a fully invisible decay mode, are displayed in red. Contrary to all other cases, the production of the heavy Higgs for this channel is via vector boson fusion as searched for at the LHC [52]. Typical cross sections range between fb and fb. The channel is shown in blue and it can yield cross sections of fb for GeV. Last is the channel. It can be sizable only for very light masses, fb for GeV, although the further decay chain of the heavy neutrinos has to be accounted for. The latter can give spectacular multileptonic final states of the heavy Higgs boson ( and ) or high jet multiplicity ones (), via and in a 2 : 1 ratio (modulo threshold effects). Further, these decays are typically seesawsuppressed and can therefore give rise to displaced vertices [53].
4.2. Inverted Hierarchy
In this subsection we discuss the inverted hierarchy case, where is the SMlike boson and a lighter Higgs boson exists.
We start once again by presenting the BRs for the nexttolightest Higgs boson in Figure 6. This time, however, this is the SMlike boson, hence predominantly from the doublets. It has the same new channels as the heavy Higgs in the normal hierarchy, the only difference being the odd Rsneutrino instead of the even one. This is simply because the inverted hierarchy can happen only for large positive values, where only the odd Rsneutrino can be light; see Figure 1. The configurations not allowed by the low energy observables or by HiggsBounds are displayed as gray points. We see that may have sizable decays into pairs of the lighter Higgs bosons, yielding jets final states. This decay is still allowed with rates up to few percent. Further, rare decays into pairs of heavy neutrinos are also present, with BRs below the permil level. This channel can give rise to rare multilepton/jets decays for the SMlike Higgs boson, which are searched for at the LHC, even in combination with searches for displaced vertices [54]. The last available channel is the decay into pairs of odd Rsneutrinos. Being the LSP, it will increase the invisible decay width and hence give largerthanexpected widths for the SMlike boson. Its rate is obviously constrained, and a precise evaluation of the allowed range is needed. It however goes beyond the scope of the present review and we postpone it to a future publication.
Regarding the lightest Higgs boson (), this will obviously decay predominantly into pairs of jets, see Figure 7. Notice that due to its large bilepton fraction it can also decay into pairs of very light RH neutrinos, at sizable rates depending on the neutrino masses. As in Figure 6, the nonallowed configurations are displayed as gray points. We see that the pattern of decays is not affected by the inclusion of the constraints, in the sense that this channel stays viable. Once again, the latter will yield multilepton/jet final state, which will be very soft and hence very challenging for the LHC. However, also in this case displaced vertices may appear.
As in the previous section, we show in Figure 8 the mixing between the Higgs mass eigenstates and the doublet fields as a function of the light Higgs mass, to show that is here rather SMlike. Once more, the gray points displayed here are excluded by the low energy observables and by HiggsBounds.
Finally, the production cross sections for the lightest Higgs boson can be evaluated. In Figure 9 we compare the direct production (for the main SM production mechanisms, gluon fusion and vector boson fusion) with the pair production via decays only via gluon fusion, . When the latter channel is kinematically open, that is, GeV, the lightest Higgs boson pair production has cross sections up to pb at the LHC at TeV, and it can give rare , , or () decays of the SMlike Higgs boson. A thorough analysis of the phenomenology of the Higgs sector in the BLSSM for the upcoming LHC run 2, based on the first investigations shown here, will be performed soon.
4.3. Enhancement of the Diphoton Rate
A feature of gaugeextended models is that new SUSYbreaking Dterms arise, which give further contributions to the sparticle masses. In the case of the model under consideration, we showed discussing (18) that these terms can be large and that they bring larger corrections to sleptons than to squarks. We already discussed how the vacuum structure of the BLSSM is affected by this. Here, we discuss the impact of the new Dterms on the Higgs phenomenology, focusing on the Higgstodiphoton decay, despite being disfavoured by most recent data [17], as an illustrative case. See [46] for further details.
To start our discussion let us briefly review the partial decay width of the Higgs boson into two photons within the MSSM and its singlet extensions. This can be written as (see, e.g., [55])corresponding to the contributions from charged SM fermions, bosons, charged Higgs, charginos, charged sleptons, and squarks, respectively. The amplitudes at the lowest order for the spin1, spin, and spin0 particle contributions can be found, for instance, in [55]. denotes the coupling between the Higgs boson and the particle in the loop and is its electric charge. In the SM, the largest contribution is given by the loop, while the toploop leads to a small reduction of the decay rate. In the MSSM, it is possible to get large contributions due to sleptons and squarks, although it is difficult to realise such a scenario in a constrained model with universal sfermion masses [56–58]. In singlet or triplet extension of the MSSM also the chargino and charged Higgs can enhance the loop significantly [59, 60]. However, this is only possible for large singlet couplings which lead to a cutoff well below the GUT scale. In contrast, it is possible to enhance the diphoton ratio in the BLSSM due to light staus even in the case of universal boundary conditions at the GUT scale. We show this by calculating explicitly the contributions of the stau: Here, and represent the Dterm contributions of the left and righthanded stau and we have neglected subleading contributions. Given that , for fixed values of the other parameters, and can be used to enhance the rate by suppressing the denominator.
We turn now to a fully numerical analysis to demonstrate the mechanism to enhance the Higgstodiphoton rate as a feature of the model with an extended gauge sector. This is a result of reducing the stau mass at the Higgs mass scale via extra Dterms as shown discussing (18). We recall here that this mechanism leaves the stop mass and hence, as we will show, the Higgstogluons effective coupling nearly unchanged. In Table 2 we have collected two possible scenarios that provide SMlike Higgs particle in the mass range preferred by LHC results displaying an enhanced diphoton rate. In the first point, the lightest even scalar eigenstate is the SMlike Higgs boson while the light bilepton is roughly twice as heavy. In Figure 10 we show that all the features arise from the extended gauge sector: it is sufficient to change only to obtain an enhanced diphoton signal and the correct dark matter relic density while keeping the mass of the SMlike Higgs nearly unchanged. The dark matter candidate in this scenario is the lightest neutralino that is mostly a bileptino (the superpartner of the bileptons). The correct abundance for is obtained due to a coannihilation with the light stau. In the second point, the SMlike Higgs is accompanied by a light scalar around GeV which couples weakly to the SM gauge bosons, compatibly with the LEP excess [61–63]. In this case, the LSP is a odd sneutrino which annihilates very efficiently due to large . This usually results in a small relic density. To get an abundance which is large enough to explain the dark matter relic, the mass of the sneutrino has to be tuned below [16]. This can be achieved by slightly increasing and by tuning the Majorana Yukawa couplings , which tends to increase the SMlike Higgs mass for the given point. It is worth mentioning that a neutralino LSP with the correct relic density in the stau coannihilation region can also be found in this scenario. Notice that both points yield rates consistent with observations in the channels (measured at the LHC) (being ), as well as an effective Higgstogluon coupling close to 1.

(a)
(b)
(c)
5. Conclusions
In this review I described the extension of the MSSM, focusing in particular on the scalar sector, described in detail. The fundamental role that the gauge kinetic mixing plays in this sector has been underlined.
The comparison to the most constraining low energy observables showed that a preferred region for the light neutrino masses exists to evade these bounds. Then, I presented a first systematic investigation of the phenomenology of the Higgs sector of this model, showing that both the normal hierarchy and the inverted hierarchy of the two lightest Higgs bosons are naturally possible in a large portion of the parameter space. Particular attention has been devoted to analysis of the new decay channels comprising both the even and odd Rsneutrinos, which are a peculiarity of the BLSSM. Based on these first findings, a thorough analysis of the Higgs sector in the BLSSM at the upcoming LHC run 2 will be soon prepared. The fit of the SMlike Higgs boson to the LHC data will also be performed with HiggsSignals [64].
Finally, I described how in the BLSSM model (and in general in gaugeextended MSSM models) the Higgstodiphoton decay can be easily enhanced. Despite being disfavoured by most recent data, this feature is a consequence of the potentially large new SUSYbreaking Dterms arising from the sector. At the same time these terms affect also the vacuum structure of the model, where naive RParityconserving configurations at the tree level could develop deeper RParityviolating global minima, or partially restore the symmetry at one loop. It is however possible to still find RParityconserving global minima on the whole parameter space, which can either accommodate an enhancement of the Higgstodiphoton decay or fit the most recent Higgs data.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The author would like to thank S. Moretti and C. H. ShepherdThemistocleous for helpful discussions in the early stages of this work. He is also really grateful to all his collaborators, and in particular to Florian Staub. He further acknowledges support from the TheorieLHC France initiative of the CNRS/IN2P3 and from the French ANR 12 JS05 002 01 BATS@LHC.
References
 J. Ellis, F. Luo, K. A. Olive, and P. Sandick, “The Higgs mass beyond the CMSSM,” The European Physical Journal C, vol. 73, article 2403, 2013. View at: Publisher Site  Google Scholar
 U. Ellwanger, C. Hugonie, and A. M. Teixeira, “The nexttominimal supersymmetric standard model,” Physics Reports, vol. 496, no. 12, pp. 1–77, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 W. Buchmüller, K. Hamaguchi, O. Lebedev, and M. Ratz, “Supersymmetric Standard Model from the heterotic string (II),” Nuclear Physics B, vol. 785, no. 12, pp. 149–209, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. Ambroso and B. A. Ovrut, “THE BL/electroweak hierarchy in smooth heterotic compactifications,” International Journal of Modern Physics A, vol. 25, no. 13, p. 2631, 2010. View at: Publisher Site  Google Scholar
 M. Ambroso and B. A. Ovrut, “The mass spectra, hierarchy and cosmology of BL MSSM heterotic compactifications,” International Journal of Modern Physics A, vol. 26, no. 9, pp. 1569–1627, 2010. View at: Publisher Site  Google Scholar
 S. Khalil and A. Masiero, “Radiative $BL$ symmetry breaking in supersymmetric models,” Physics Letters B, vol. 665, no. 5, pp. 374–377, 2008. View at: Publisher Site  Google Scholar
 P. Fileviez Perez and S. Spinner, “Fate of $R$ parity,” Physical Review D, vol. 83, no. 3, Article ID 035004, 7 pages, 2011. View at: Publisher Site  Google Scholar
 V. Barger, P. Fileviez Perez, and S. Spinner, “Minimal gauged $U{(1)}_{BL}$ model with spontaneous R parity violation,” Physical Review Letters, vol. 102, Article ID 181802, 2009. View at: Publisher Site  Google Scholar
 J. Pelto, I. Vilja, and H. Virtanen, “Leptogenesis $BL$ in gauged supersymmetry with MSSM Higgs sector,” Physical Review D, vol. 83, Article ID 055001, 2011. View at: Publisher Site  Google Scholar
 K. S. Babu, Y. Meng, and Z. Tavartkiladze, “New ways to leptogenesis with gauged B−L symmetry,” Physics Letters B, vol. 681, no. 1, pp. 37–43, 2009. View at: Publisher Site  Google Scholar
 B. Holdom, “Two U(1)'s and ϵ charge shifts,” Physics Letters B, vol. 166, no. 2, pp. 196–198, 1986. View at: Publisher Site  Google Scholar
 K. S. Babu, C. Kolda, and J. MarchRussell, “Implications of generalized ZZ′ mixing,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 57, no. 11, pp. 6788–6792, 1998. View at: Publisher Site  Google Scholar
 F. del Aguila, G. D. Coughlan, and M. Quirós, “Gauge coupling renormalisation with several U(1) factors,” Nuclear Physics B, vol. 307, no. 3, pp. 633–648, 1988. View at: Publisher Site  Google Scholar
 F. del Aguila, J. Gonzalez, and M. Quiros, “Renormalization group analysis of extended electroweak models from the heterotic string,” Nuclear Physics B, vol. 307, no. 3, pp. 571–632, 1988. View at: Publisher Site  Google Scholar
 B. O'Leary, W. Porod, and F. Staub, “Mass spectrum of the minimal SUSY $BL$ model,” Journal of High Energy Physics, vol. 2012, no. 5, article 042, 2012. View at: Google Scholar
 L. Basso, B. O'Leary, W. Porod, and F. Staub, “Dark matter scenarios in the minimal SUSY $BL$ model,” Journal of High Energy Physics, vol. 2012, no. 9, article 54, 2012. View at: Publisher Site  Google Scholar
 V. Khachatryan, A. M. Sirunyan, A. Tumasyan et al., “Precise determination of the mass of the Higgs boson and tests of compatibility of its couplings with the standard model predictions using proton collisions at 7 and 8 TeV,” The European Physical Journal C, vol. 7, article 212, 2014. View at: Publisher Site  Google Scholar
 J. Adam, X. Bai, A. M. Baldini et al., “New constraint on the existence of the ${\mu}^{+}\to {e}^{+}\gamma $ decay,” Physical Review Letters, vol. 110, Article ID 201801, 2013. View at: Publisher Site  Google Scholar
 U. Bellgardt, G. Otter, R. Eichler et al., “Search for the decay μ→ e^{+}e^{+}e^{−},” Nuclear Physics B, vol. 299, no. 1, pp. 1–6, 1988. View at: Publisher Site  Google Scholar
 CMS Collaboration, “Search for the standard model Higgs boson in the H to WW to lnujj decay channel in pp collisions at the LHC,” Tech. Rep. CMSPASHIG13027, CERN, Geneva, Switzerland, 2012. View at: Google Scholar
 S. Chatrchyan, V. Khachatryan, A. M. Sirunyan et al., “Search for a standardmodellike Higgs boson with a mass in the range 145 to 1000 GeV at the LHC,” The European Physical Journal C, vol. 73, article 2469, 2013. View at: Publisher Site  Google Scholar
 W. Porod, “SPheno, a program for calculating supersymmetric spectra, SUSY particle decays and SUSY particle production at ${e}^{+}{e}^{}$ colliders,” Computer Physics Communications, vol. 153, no. 2, pp. 275–315, 2003. View at: Publisher Site  Google Scholar
 W. Porod and F. Staub, “SPheno 3.1: extensions including flavour, CPphases and models beyond the MSSM,” Computer Physics Communications, vol. 183, no. 11, pp. 2458–2469, 2012. View at: Publisher Site  Google Scholar
 F. Staub, “Sarah,” http://arxiv.org/abs/0806.0538. View at: Google Scholar
 F. Staub, “From superpotential to model files for FeynArts and CalcHep/CompHep,” Computer Physics Communications, vol. 181, no. 6, pp. 1077–1086, 2010. View at: Publisher Site  Google Scholar
 F. Staub, “Automatic calculation of supersymmetric renormalization group equations and loop corrections,” Computer Physics Communications, vol. 182, no. 3, pp. 808–833, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Staub, “SARAH 3.2: dirac gauginos, UFO output, and more,” Computer Physics Communications, vol. 184, no. 7, pp. 1792–1809, 2013. View at: Publisher Site  Google Scholar
 F. Staub, “SARAH 4 : a tool for (not only SUSY) model builders,” Computer Physics Communications, vol. 185, no. 6, pp. 1773–1790, 2014. View at: Publisher Site  Google Scholar
 L. Basso, A. Belyaev, D. Chowdhury et al., “Proposal for generalised supersymmetry les Houches Accord for seesaw models and PDG numbering scheme,” Computer Physics Communications, vol. 184, no. 3, pp. 698–719, 2013. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Staub, “Exploring new models in all detail with SARAH,” http://arxiv.org/abs/1503.04200. View at: Google Scholar
 F. Staub, T. Ohl, W. Porod, and C. Speckner, “A tool box for implementing supersymmetric models,” Computer Physics Communications, vol. 183, no. 10, pp. 2165–2206, 2012. View at: Publisher Site  Google Scholar
 W. Porod, F. Staub, and A. Vicente, “A flavor kit for BSM models,” The European Physical Journal C, vol. 74, article 2992, 2014. View at: Publisher Site  Google Scholar
 P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, and K. E. Williams, “HiggsBounds: confronting arbitrary Higgs sectors with exclusion bounds from LEP and the Tevatron,” Computer Physics Communications, vol. 181, no. 1, pp. 138–167, 2010. View at: Publisher Site  Google Scholar
 P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, and K. E. Williams, “HiggsBounds 2.0.0: confronting neutral and charged Higgs sector predictions with exclusion bounds from LEP and the Tevatron,” Computer Physics Communications, vol. 182, no. 12, pp. 2605–2631, 2011. View at: Publisher Site  Google Scholar
 P. Bechtle, O. Brein, S. Heinemeyer et al., “Recent developments in HiggsBounds and a preview of HiggsSignals,” PoS CHARGED, http://arxiv.org/abs/1301.2345. View at: Google Scholar
 P. Bechtle, O. Brein, S. Heinemeyer et al., “HiggsBounds4 : improved tests of extended Higgs sectors against exclusion bounds from LEP, the Tevatron and the LHC,” The European Physical Journal C, vol. 74, article 2693, 2014. View at: Publisher Site  Google Scholar
 R. M. Fonseca, M. Malinsky, W. Porod, and F. Staub, “Running soft parameters in SUSY models with multiple U(1) gauge factors,” Nuclear Physics B, vol. 854, no. 1, pp. 28–53, 2012. View at: Publisher Site  Google Scholar
 P. H. Chankowski, S. Pokorski, and J. Wagner, “Z′ and the AppelquistCarrazzone decoupling,” The European Physical Journal C—Particles and Fields, vol. 47, no. 1, pp. 187–205, 2006. View at: Publisher Site  Google Scholar
 W. Abdallah, S. Khalil, and S. Moretti, “Double Higgs peak in the minimal SUSY $BL$ model,” Physical Review D, vol. 91, no. 1, Article ID 014001, 6 pages, 2015. View at: Publisher Site  Google Scholar
 G. Aad, B. Abbott, J. Abdallah et al., “Search for highmass dilepton resonances in $pp$ collisions at $\sqrt{s}=8$ TeV with the ATLAS detector,” Physical Review D, vol. 90, Article ID 052005, 2014. View at: Publisher Site  Google Scholar
 V. Khachatryan, A. M. Sirunyan, A. Tumasyan et al., “Search for physics beyond the standard model in dilepton mass spectra in protonproton collisions at $\sqrt{\stackrel{}{s}}=8$ TeV,” Journal of High Energy Physics, vol. 2015, no. 4, article 025, 2015. View at: Publisher Site  Google Scholar
 G. Cacciapaglia, C. Csáki, G. Marandella, and A. Strumia, “The minimal set of electroweak precision parameters,” Physical Review D, vol. 74, no. 3, Article ID 033011, 12 pages, 2006. View at: Publisher Site  Google Scholar
 L. Basso, K. Mimasu, and S. Moretti, “Nonexotic $Z\text{'}$ signals in ${l}^{+}{l}^{}$, $bb$ and $t\stackrel{}{t}$ final states at the LHC,” Journal of High Energy Physics, vol. 2012, no. 11, article 060, 2012. View at: Publisher Site  Google Scholar
 M. Hirsch, H. V. KlapdorKleingrothaus, and S. G. Kovalenko, “BLviolating masses in softly broken supersymmetry,” Physics Letters B, vol. 398, no. 34, pp. 311–314, 1997. View at: Publisher Site  Google Scholar
 Y. Grossman and H. E. Haber, “Sneutrino mixing phenomena,” Physical Review Letters, vol. 78, no. 18, pp. 3438–3441, 1997. View at: Publisher Site  Google Scholar
 L. Basso and F. Staub, “Enhancing $h\to \gamma \gamma $ with staus in supersymmetric models with an extended gauge sector,” Physical Review D, vol. 87, no. 1, Article ID 015011, 6 pages, 2013. View at: Publisher Site  Google Scholar
 A. Hammad, S. Khalil, and S. Moretti, “Higgs boson decays into γγ and Zγ in the MSSM and BLSSM,” http://arxiv.org/abs/1503.05408. View at: Google Scholar
 J. E. CamargoMolina, B. O'Leary, W. Porod, and F. Staub, “Stability of $R$ parity in supersymmetric models extended by $U{(1)}_{BL}$,” Physical Review D, vol. 88, Article ID 015033, 2013. View at: Publisher Site  Google Scholar
 A. Abada, M. E. Krauss, W. Porod, F. Staub, A. Vicente, and C. Weiland, “Lepton flavor violation in lowscale seesaw models: SUSY and nonSUSY contributions,” Journal of High Energy Physics, vol. 2014, no. 11, article 048, 2014. View at: Publisher Site  Google Scholar
 A. Vicente, “Lepton flavor violation beyond the MSSM,” http://arxiv.org/abs/1503.08622. View at: Google Scholar
 L. Basso, S. Moretti, and G. M. Pruna, “Phenomenology of the minimal $BL$ extension of the standard model: the Higgs sector,” Physical Review D, vol. 83, Article ID 055014, 2011. View at: Publisher Site  Google Scholar
 CMS Collaboration, “Search for invisible decays of Higgs bosons in the vector boson fusion production mode,” Tech. Rep. CMSPASHIG14038, CERN, Geneva, Switzerland, 2015, http://cds.cern.ch/record/2007270. View at: Google Scholar
 L. Basso, A. Belyaev, S. Moretti, and C. H. ShepherdThemistocleous, “Phenomenology of the minimal $BL$ extension of the standard model: ${Z}^{\prime}$ and neutrinos,” Physical Review D, vol. 80, no. 5, Article ID 055030, 14 pages, 2009. View at: Publisher Site  Google Scholar
 L. Basso, A. Belyaev, J. Fiaschi, S. Moretti, I. Tomalin, and M. Thomas, In preparation.
 A. Djouadi, “The anatomy of electroweak symmetry breaking: tome I: the Higgs boson in the Standard Model,” Physics Reports, vol. 457, no. 1–4, pp. 1–216, 2008. View at: Publisher Site  Google Scholar
 M. Carena, S. Gori, N. R. Shah, and C. E. M. Wagner, “A 125 GeV SMlike Higgs in the MSSM and the $\gamma \gamma $ rate,” Journal of High Energy Physics, vol. 2012, no. 3, article 014, 2012. View at: Publisher Site  Google Scholar
 U. Ellwanger, “A Higgs boson near 125 GeV with enhanced diphoton signal in the NMSSM,” Journal of High Energy Physics, vol. 2012, no. 3, article 044, 2012. View at: Google Scholar
 R. Benbrik, M. Gomez Bock, S. Heinemeyer, O. Stål, G. Weiglein, and L. Zeune, “Confronting the MSSM and the NMSSM with the discovery of a signal in the two photon channel at the LHC,” European Physical Journal C, vol. 72, article 2171, 2012. View at: Publisher Site  Google Scholar
 K. SchmidtHoberg and F. Staub, “Enhanced $h\to \gamma \gamma $ rate in MSSM singlet extensions,” Journal of High Energy Physics, vol. 2012, no. 10, article 195, 2012. View at: Publisher Site  Google Scholar
 A. Delgado, G. Nardini, and M. Quiros, “Large diphoton Higgs rates from supersymmetric triplets,” Physical Review D, vol. 86, no. 11, Article ID 115010, 7 pages, 2012. View at: Publisher Site  Google Scholar
 ALEPH Collaboration, DELPHI Collaboration, L3 Collaboration, OPAL Collaboration, and The LEP Working Group for Higgs Boson Searches, “Search for the Standard Model Higgs boson at LEP,” Physics Letters B, vol. 565, Article ID 0306033, pp. 61–75, 2003. View at: Publisher Site  Google Scholar
 G. Bélanger, U. Ellwanger, J. F. Gunion, Y. Jiang, S. Kraml, and J. H. Schwarz, “Higgs bosons at 98 and 125 GeV at LEP and the LHC,” Journal of High Energy Physics, vol. 2013, no. 1, article 069, 2013. View at: Publisher Site  Google Scholar
 M. Drees, “Supersymmetric explanation of the excess of Higgslike events at the LHC and at LEP,” Physical Review D, vol. 86, Article ID 115018, 2012. View at: Publisher Site  Google Scholar
 P. Bechtle, S. Heinemeyer, O. Stal, T. Stefaniak, and G. Weiglein, “HiggsSignals: confronting arbitrary Higgs sectors with measurements at the Tevatron and the LHC,” The European Physical Journal C, vol. 74, article 2711, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Lorenzo Basso. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.