Interpreting 95 GeV di-photon/ b ¯ b excesses as a lightest Higgs boson of the MRSSM

The Minimal R-symmetric Supersymmetric Standard Model (MRSSM) is a well motivated BSM model which can accommodate the observed 125 GeV Higgs boson in agreement with electroweak precision observables, in particular with the W boson mass and T parameter. In the 2016 paper we showed that the SM-like 125 GeV Higgs state can be also realised as the second-to-lightest scalar of the MRSSM, leaving room for another sub-100 GeV state. Motivated by the recent ATLAS and CMS observation of the di-photon excess at a mass of around 95 GeV we investigate the possibility whether this could be the lightest CP-even MRSSM scalar in a variation of our benchmarks presented in the 2016 work. We show that such a state can also simultaneously explain the excess in the b ¯ b final state observed around the same mass value at LEP. Due to the R-symmetric nature of the model, a light singlet-like Higgs state leads necessarily to a light bino-singlino Dirac dark matter candidate, which can give a correct relic density while evading current experimental bounds. Dark matter and LHC searches place further bounds on this scenario and point to parameter regions which are viable and of interest for the LHC Run III and upcoming dark matter experiments.


Introduction
After the discovery of a Standard Model (SM)-like Higgs boson at the LHC in 2012, it remains an open question whether this is the only fundamental scalar particle or whether it actually heralds a beyond the SM theory (BSM) with an extended scalar sector.Such sectors are notorious in many BSM theories that aim to address at least some unsolved puzzles of the SM, like the hierarchy of scales, neutrino masses, dark matter, baryogenesis etc.
In the absence of any other direct experimental BSM signal (new coloured particles, heavy fermions etc.), a careful study of the Higgs sector properties might be the only way for inferring properties of new physics in coming years.Any BSM scenario must allow to accommodate a SM-like 125 GeV state (h SM 125 ) with properties consistent with current measurements.While the couplings of the h SM 125 state to gauge bosons and third generation fermions and muons have been checked to be SM-like within experimental errors, there is still some room for small deviations.And the h SM 125 self-interactions are still unconstrained.However, an intriguing question is whether there are additional scalar particles, possibly even with smaller mass.Searches for scalar particles below 125 GeV have been performed in the past at LEP and Tevatron and now at the LHC.Results based on Run 1 and the first years of Run 2 data collected by the CMS in di-photon channel showed an excess of 2.9σ (local) for a hypothetical mass of 95.3 GeV [1].Recently CMS confirmed the excess of di-photon events at the 2.9σ level for a mass of 95. 4 GeV based on full Run 2 data set [2].ATLAS search found 1.7σ (local) around 95 GeV [3].It is worth recalling an excess around 95 GeV in the b b channel in LEP searches [4].An analogous small excess has also been seen in the τ τ channel in CMS [5].
These findings have triggered speculations that at least some of those excesses could arise from the production of a new particle (see for example [6,7] and references therein).
In the 2016 paper [8] we considered a possibility of a light (below 125 GeV) singlet-like scalar in the Minimal R-symmetric Supersymmetric model (MRSSM) [9].The MRSSM is an attractive alternative to the MSSM, elegantly solving some of its shortcomings while also significantly altering its phenomenology.For example, the SUSY flavour problem is alleviated as R-symmetry forbids left-right squark mixing which is responsible for many flavour-violating interactions.Scalar quarks can be lighter than in the MSSM as the collider limits for squark masses are weakened [10][11][12][13][14][15] when simultaneously the electroweak and Higgs sectors contain new interactions which can push the SM-like Higgs boson mass up to the observed value for smaller top-squark masses [8,[16][17][18].Those interactions can also contribute to the W-boson mass [8,17,19].The MRSSM also contains various possibilities to explain dark matter [20][21][22] as well as predicts new color-octet scalars [23][24][25][26][27] and Dirac gauginos [28][29][30].The lepton flavour properties have also been analysed in Refs.[31,32], and in relation to the muon g −2 in Ref. [33].
However, it remained an open question whether the aforementioned experimental deviations could be accommodated within the MRSSM.In this work we will show that this is indeed the case.
The paper is structured as follows.In Sec. 2 we recap experimental hints pointing to the existence of a ∼95 GeV scalar resonance.In Sec. 3 we review the basics of the MRSSM and describe the model setup explaining the LHC and LEP excesses -we also propose two benchmark points.In Sec. 4 we confront the benchmark points of Sec. 3 with Higgs physics, dark matter and general LHC constraints, before concluding in Sec. 5.

Hints of a new scalar with the mass between 95 and 100 GeV
In this section we briefly summarize the experimental hints for a ∼ 95 GeV scalar resonance.First, there is a long standing LEP anomaly observed in the Zb b final state in the 95 ≲ m b b ≲ 100 invariant mass window [4].This anomaly can be explained for example by a scalar state s with a mass of 98 GeV whose combined production and branching ratio is roughly an order of magnitude smaller than that of a hypothetical SM-like Higgs h SM 98 of the same mass [34] µ Hints of such a state emerged also in the 8 and 13 TeV LHC data [1][2][3].A naive combination of ATLAS and CMS results points to the 3.1σ (local) excess with mass 95.4 GeV and signal strength [35] µ Finally, it is worthwhile mentioning that there is also a hint of an analogous excess in the τ τ channel with µ τ τ = 1.2 ± 0.5 [5].
As the LEP excess is very broad, both µ Zb b and µ γγ can be addressed simultaneously by a single state.Such scenario is possible in many BSM models with an extended Higgs sector.In non-supersymmetric models it has been analysed in a plethora of models (see for example [35][36][37]).In SUSY however, it cannot be realized in its minimal version [38] but can be accommodated in extended SUSY, like the NMSSM (e.g.[7,34]), µνSSM (e.g.[39,40]) or, as we will show, the MRSSM.

A light CP-even scalar in the MRSSM
MRSSM is very well motivated, with phenomenology distinctly different from the MSSM as it contains a continuous, unbroken at the low scale U R (1) R-symmetry under which both superfields and Grassmannian variables θ are charged [41,42].At the level of the Lagrangian, the U R (1) transformation of the coordinate θ → e iφ θ implies that the superpotential W has R-charge +2.Assigning zero R-charges to the SM components of supermultiplets (in analogy to R-parity) forbids then Majorana gaugino and higgsino mass terms.Therefore an enlarged field content is needed to account for non-vanishing of those masses.
In the MRSSM the standard MSSM matter, Higgs and gauge superfields are augmented by the R-charge 0 adjoint chiral superfields Ô, T , Ŝ for each gauge The R-charges of the superfields and the corresponding bosonic and fermionic components.
sector, respectively, and two Higgs iso-doublet superfields Rd,u with R-charge 2 (see Tab. 1 for summary of the MRSSM particle content including the R-charge assignements).The MRSSM superpotential takes the form of where standard notation is used for the MSSM-like fields.This is then supplemented by a D-term SUSY breaking Lagrangian which generates the Dirac mass terms for gauginos.Superfields Ĥu,d , S and T present in Eq. 3 mix to form physical Higgs boson.The mass matrix of the CP-even neutral Higgs bosons in the weak basis (ϕ d , ϕ u , ϕ S , ϕ T ) is then given by [17] with the sub-matrices .
Here m 2 A = 2B µ / sin 2β (with the usual definition of B µ ), we use shorthand notation The vacuum expectation values v d,u,S,T of neutral scalars are assumed real.The MSSM-like M MSSM and the singlet-triplet-like M 22 submatrices mix through the M 21 term.
In the analogy to the MSSM, the SM-like Higgs state is considered to be ϕ u -dominated while the ϕ T is assumed heavy since a large soft breaking mass term m T is a natural and sufficient way to suppress the tree-level triplet contribution to the W boson mass and ρ parameter.This leads to two phenomenologically viable scenarios, with either the SM-like Higgs being the lightest state [17] or the next-to-lightest [8], with a low mass singlet scalar.
To understand how one can realize a scenario with m h 1 ≈ 95 and m h 2 ≈ 125 GeV it is instructive to consider the limit in which not only the triplet decouples (leading to a Higgs sector field content as in the singlet extended 2HDM [43][44][45], albeit with non-trivial relations between couplings) but also in which the MSSM-like pseudoscalar Higgs mass M A and the value of tan β = v u /v d become large (this last approximation will turn out also to be a necessary condition to fit LEP and LHC signal strengths).In this limit it is enough to focus on the 2x2 sub-matrix corresponding to the (ϕ u , ϕ S ) fields only, which reads where ∆m 2 rad denotes expected large quantum corrections to the mass of the SM-like Higgs.A setup with a ∼ 95 GeV scalar therefore requires the matrix element (2,2) small compared to the (1,1) one.This requires and since v u ≈ v = 246 GeV for large tan β, the coupling λ u therefore must be also very small.The off-diagonal matrix element must remain small in order not to disturb the properties of the SM-like Higgs state h 2 .As we will see in the next section, the dark matter data will impose additional constraints among these parameters.
Let us first comment on original benchmark points BMP4-6 of Ref. [8].While the mixing and masses are more-or-less appropriate to fix the LEP excess, none of the points exhibits an enhancement in the γγ channel.Moreover, BMP5 and BMP6 are nowadays excluded by dark matter direct search at LUX-ZEPPELIN [46], while BMP4 is in tension with p-value of 0.083 (as computed by micrOMEGAs).BMP4 and BMP6 are also excluded by stau searches (the most constraining analysis are [47] and [48], respectively) while BMP5 is excluded by jets + MET search [49].
In Tab. 2 we show two representative benchmark points, with decay patterns of the two lightest CP-even Higgses listed in Tab. 3.These benchmarks are characterised by 1• mostly a singlet-like Higgs state h 1 with the mass of 95.4 GeV and a SM-like state h 2 with m h 2 ≈ 125. 25 GeV, • where h 1 fits the LEP excess with (approximating production via partial widths) • and where h 1 is produced at the LHC with • due to small bino-singlino mass parameter M D B both benchmarks feature a light neutralino, which can serve as a dark matter candidate, • however BMP7 and BMP8 differ in how the relic density of DM is achieved, see next section.
Since colour and electrically charged particles are (apart for right handed staus in BMP8) fairly heavy (≳ 0.5 TeV), the µ γγ is not enhanced relative to µ Zb b via contribution from new particles in the loop.Desired values of µ Zb b and µ γγ are rather achieved almost exclusively via Higgs mixing.The ratios of Γ(h 1 → ZZ)/Γ(h SM 95.4 → ZZ), Γ(h 1 → gg)/Γ(h SM 95.4 → gg), and Γ(h 1 → γγ)/Γ(h SM 95.4 → γγ) are around 0.13 and have a simple scaling with the S − H u mixing matrix element.This makes them fully correlated.The non-universality, needed to explain the difference between µ Zb b and µ γγ comes from ratios of branching ratios in Eqs. 9 and 10.Since the total width of scalar states between 90 -130 GeV is dominated by the decay to b b, the BR(h 1 → b b) is fairly insensitive to the mixing effects which cancel between the nominator and denominator in the branching ratio, giving BR(h there is a difference between suppression for the width generated via the H d and H u admixtures).This is no longer the case for the BR(h 1 → γγ) which gets enhanced in relation to BR(h SM 95.4 → γγ) due to a weaker suppression of Γ(h 1 → γγ)/Γ(h SM 95.4 → γγ) in relation to total width of h 1 .This setup is similar to the one described in [7] in the NMSSM.However, the required pattern of masses and mixings is far more difficult to realise in the MRSSM since its constraining nature connects different sectors of the model which means that experimental limits from the dark matter and general collider phenomenology have impact on the Higgs sector.
To understand how the desired pattern of partial widths described above is achieved for parameter points given in Tab. 2 we make a set of 2d scans around BMP7.We focus here on masses of two lightest CP-even Higgses (Fig. 1) and LEP and LHC signal strengths (Fig. 2), relegating the analysis of experimental limits, like dark matter or the consistency of the Higgs sector with experimental data, to Sec. 4. In our setup, where the lightest Higgs is mainly a singlet (with mass given approximately by m 2 S + 4(M D B ) 2 as show in Fig. 1a) and µ u ≳ 0.5 TeV to avoid direct searches of SUSY particles at the LHC (see discussion in Sec.4.3), λ u is constrained to a very small range of λ u ∈ [−0.05, −0.0025] (see Figs. 1d and 1e) and is strongly correlated with µ u (Fig. 1e), as it appears as a λ u µ eff,− u product in the (2, 1) entry of the mass matrix in Eq. 7. The light CP-even state with mass of 95 GeV could be also achieved for a larger value of |λ u | when larger mixing with the SM-like state compensates an increased value of m 2 S + 4(M D B ) 2 , but such points would then be excluded by the Higgs sector analysis in Sec.4.1.Therefore, for a given range of λ u , the parameter µ u influences nontrivially mainly the SM-like Higgs mass via one-loop corrections, where µ u can appear without  the λ u parameter.The situation is somewhat similar for LHC and LEP signal strengths.With a small allowed region of λ u (Fig. 2a) and rather simple functional dependence of µ Zb b (Figs. 2b, 2c and 2d), the non-trivial dependence appears mainly for µ γγ as it is loop induced.For other combinations of parameters appearing in Eq. 7 but not shown in Fig. 1 or Fig. 2 masses or signal strength, respectively, are mostly correlated not leading to any non-trivial relations.We should mention that neither BMP7 nor BMP8 fit the small µ τ τ excess at 95 GeV.This excess cannot be achieved via the mixing while fitting at the same time µ Zb b as in the THDM type II, where h 1 couplings to b b and τ τ are fully correlated (from the point of view of the pure Higgs sector, our setup is analogues to the type II version of S2HDM from Ref. [6]).This can be revisited in the future when FlexibleSUSY will be able to compute loop corrections to the τ τ decay, especially in a scenario with light staus, although expecting such large corrections seems unreasonable.

Confronting with experimental observations
In this section we discuss the phenomenological properties of benchmark points proposed in Sec. 3. In particular, we asses their viability in light of constraints coming from Higgs sector measurements, dark matter properties and general collider phenomenology.The properties of predicted Higgs bosons are checked against experimental data using HiggsTools v1.1.4[65], whereas dark matter properties and LHC constraints are checked using the combination of micrOMEGAs v6.0 [66][67][68] and SModelS v2.3.3 [69][70][71][72], based on SARAH generated CalcHEP model file [58,73].We show a selection of 2-dimensional scans around the BMPs to explain which values of parameters are necessary to avoid current experimental limits.Finally, in Sec.4.3 we comment briefly on prospects for discovering of the MRSSM at future collider experiments.

Higgs sector
The Higgs sector of the model was checked against experimental data using FlexibleSUSY interface to HiggsTools. 2 The validity of any parameter point is assessed by comparing its χ 2 with the minimal χ 2 in the SM.In the SM we find the minimal value χ 2 SM = 151.55 at the Higgs boson mass of 125.25 GeV (all of the numbers use HiggsSignals database v1.1 and assume 3% uncertainty on the Higgs mass prediction in the BSM model).The allowed region in the parameter space at 95% C.L. corresponds to ∆χ 2 ≡ χ 2 − χ 2 SM ≲ 5.99. 3 Tab.4 shows the p-values for both BMPs computed from this ∆χ 2 for 2-degrees of freedom. 4Both points are within 2σ, though more than 1σ away, so if the LHC anomalies persist while the measurements of SM-like Higgs boson improve, this solution might become disfavoured.The points also predict a roughly SM-like h 2 Zγ coupling, which is of no consequence now [77] but might become important once the precision of this measurement increases and statistically significant deviation is observed.
In Fig. 3 we show example regions allowed by measurements of SM-like Higgs properties for selected parameters relevant to properties of the 95 GeV scalar.While properties of the second-to-lightest Higgs boson are determined mostly by other parameters, which were chosen to make it mostly SM-like, the mixing with the singlet can destabilize them.Since some mixing is needed to give the 95 GeV state its couplings to SM particles, there is obviously also an   White regions is where no spectrum could be generated.
upper limit on it.This is seen for example in Fig. 3b where the increase of M D B , which appears in the (2,1) element of the mass matrix in Eq. 7 as a combination √ 2λ u µ eff,− u + g 1 M D B , has to be compensated by a larger value of |λ u |.
Properties of the non-SM Higgs with the mass of 95 GeV were checked using a HiggsBounds component of HiggsTools (database version 1.4).This state is mostly constrained by the wanted 95 GeV LEP anomaly [4], with the ratio of production cross sections to the 95% limit (referred by HiggsBounds as obsRatio) being 0.4 (Tab.4).In Fig. 4 we show example 2d regions around BMP7 allowed by HiggsBounds.When the admixture of SM-like Higgs in the lightest state becomes too large the points start to get excluded by searches of SM or non-SM like Higgses (e.g.[78] or [79]) as the lightest state couples more-and-more strongly to SM particles and the second-to-lightest Higgs starts to deviate from the SM.
Since the remaining charged and CP-even neutral scalar states are fairly heavy, they evade any experimental bounds.The only Higgs-like states which are light are the lightest CP-odd scalars with masses m a = 24.85(BMP7) and 54.20 (BMP8) GeV.Those states are fairly long lived with partial widths 10 −12 − 10 −11 GeV and couple almost exclusively to γγ.Thus they cannot easily be produced at pp and e + e − colliders via the γγ channel, since the production cross-section is extremely small and therefore they easily evade all experimental bounds.

Dark matter sector
As explained in Sec. 3, the light-singlet setup with the bino-singlino mass parameter m D B ∼ 60 GeV in the MRSSM inevitable leads to a light (lighter than around 60 GeV) dark matter candidate (LSP), making the DM relic density and direct detection limits important constraints.The two BMPs in Tab. 2 were selected such as to represent different ways how the correct relic density is achieved.For BMP7 the LSP correct relic density is achieved by the s-channel resonant LSP pair annihilation into Z bosons, while in BMP8 stau is light enough to annihilate LSP pairs via t-channel τR exchange into tau leptons.The direct detection expectations crucially depend on the relation between the µ u parameter and the first generation squark masses.
We compute dark matter observables with micrOMEGAs using SLHA output generated by FlexibleSUSY (FlexibleSUSY output for the MRSSM is fully compatible with SARAH generated CalcHEP/micrOMEGAs model with the caveat that one has to call the output SLHA file SPheno.spc.MRSSM).

Relic density
As we said, the correct relic density in the setup with m χ ≲ m Z can be achieved in two ways: s-channel annihilation via resonant Z (BMP7) or annihilation via t-channel τR (BMP8): BMP7: in this scenario with heavy staus, the correct relic density can only be achieved if m χ 1 ∼ m Z /2.In the MRSSM the Z boson couples to neutralinos only via their (R-)Higgsino admixture.This is controlled by tan β and the relevant µ eff,± i parameter (where the dependence on tan β is negligible so long as tan β is not ∼ 1).With M D B ≲ m Z and the size of the offdiagonal mixing determined by the electroweak scale this puts a strong bounds on the size of µ eff,− u .The smallness of µ eff,− u ≈ µ i + 1 √ 2 λ i v S can be achieved by making both elements small or by relaying on mutual cancellation.Since v S is large in the setup with light singlet, it forces λ u to be relatively small if no artificial cancellation between µ u and λ u v S is enforced.The generation squarks leading to a strong correlation between µ u and m q as seen in Fig. 6b.Even in the case when Higgsino component is not relevant for the relic density, it still plays a role in direct detection, also for µ u > 1 TeV leading to a µ um q correlation also in the case of BMP8.

Direct collider constraints
In Tab. 2 we list masses of selected BSM particles.The two light states in our setup which are in danger of being excluded by direct searches at the LHC are the lightest ρ + chargino, right handed staus and squarks.The remaining, collider relevant states (including MRSSM specific ones like Dirac gluinos or color octet scalars) were chosen to be heavy as their masses do not influence observables we are interested in this work.
Collider limits from direct production of BSM particles at the LHC were checked with SModelS (analysis database v2.3.0,official) using output provided by micrOMEGAs. 5s shown in Fig. 7a, there is a lower limit on µ u of around 500 GeV, making it low enough to allow the values of µ u necessary for a correct DM relic density.In Sec.4.2.2 we have determined the first generation squark masses needed in the case of µ u ∼ 650 GeV.Such masses are currently slightly above the current experimental limit, as shown in Fig. 7b.We note that the limits in Fig. 7b are stronger than actually needed, as we assumed there that all masses of first 2 generation squarks (m Q , m D , m U ) are equal.In reality, as can be seen in Tab. 2, we require that only some of them are small.Finally, as is important in the case of BMP8, right handed staus with 100% branching ratio τR → τ χ 0 are excluded at 95% CL for masses < 89.8 [81].At LHC they are mostly Regions allowed (green) and excluded (red) by direct searches of BSM particles at 95% C.L. as given by SModelS.White regions is where no spectrum could be generated.Since in Fig. 7b and 7c all first and second generation soft squark masses were set equal, m Q;1,2 = m U ;1,2 = m D;1,2 ≡ m q , to facilitate scanning, none of the points corresponds directly to BMP7 or BMP8.unconstrained [82] (which is also the result reported by SModelS).
In the future, light scalars will be searched for not only at the LHC but also at the planned e + e − colliders.Considered e + e − machines like CLIC, ILC, FCCee, CEPC and C 3 would allow to precisely measure properties of the 95 GeV state in the "Higgs-strahlung" process [83][84][85][86][87].The search of light scalar in the "Higgs-strahlung" process was also defined as on of focus topics for the ECFA study on Higgs/Top/EW factories [88].Definite disentanglement of the MRSSM from other models featuring a singlet extended Higgs sector would however require a discovery of some of the MRSSM specific particles like R-Higgses.They would be pair produced at both lepton and hadron colliders decaying into pairs of R-charge 1 SUSY particles, e.g.R − → χχ − or R 0 → lL l * R [89].This would result in a interesting but experimentally very challenging signature with long decay chains and multi-particle final states.If there are other R-charge 1 particles than neutralinos one would expect to observe their direct production first if kinematically allowed.Therefore the other option to narrow down the MRSSM is by focusing on its N = 2 gauge sector by looking for sgluons or the Dirac vs. Majorana nature of neutralinos and gluinos.Sgluons appear in few BSM extensions and are often studied in the context of effective models (see for example [90]).However the distinct feature of the MRSSM is that the lightest sgluon, which decays usually directly to top quarks [26], is a pseudoscalar which could be identified by the top quark decay product distribution.Meanwhile, the Dirac nature of gauginos could be established by looking at their direct production or by their impact on production of other BSM states (via an intermediate s or t channel Dirac gaugino) [14,15,28,30,91].

Summary and conclusions
The Minimal R-symmetric Supersymmetric Standard Model is a highly predictive type of supersymmetric model, featuring an unbroken U R (1) R-symmetry at the electroweak scale.Its constraining nature alters the phenomenology and makes it different from popular SUSY models like the MSSM or the NMSSM, due to a presence of Dirac gauginos, no left-right sfermion mixing, colour octet scalars and other distinct features.In the past it was show that this is reflected in a very distinct Higgs, lepton and quark flavour violation, muon magnetic moment and general collider phenomenologies.
Motivated by emerging hints of a low lying scalar resonance we have demonstrated in this work that the MRSSM can accommodate both of the excesses observed at around 95 GeV at LEP in e + e − → Zb b and LHC in pp → γγ.While the setup of the Higgs sector is at least partially similar to a set of singlet-extended 2HDMs, the non-trivial issue is how to realise the necessary pattern of masses and Higgs mixing within the constrains coming from different sectors of the MRSSM.This is especially obvious in comparison with non-supersymmetric models, where one has a complete freedom in setting masses and mixing of Higgs bosons independently of other parameters.In contrast, in the MRSSM the parameters controlling the Higgs mass matrix are intrinsically connected to dark matter sector.Light CP-even Higgs leads to an even lighter dark matter candidate.Relic density and direct detection experiments determine the MRSSM's bino-singlino mass parameter M D B to be ∼45 GeV and fix a relation between up-(R)Higgsino mass parameter µ u and the first generation squark masses.The µ u and M D B parameters in turn influence the tree-level mixing between the light scalar state and the SM-like Higgs, which implies that a simultaneous successful fit to all these observables is a non-trivial exercise.
In this work we have identified a region of parameter space where the LEP and LHC excesses and dark matter relic density can be simultaneously accommodated while evading current experimental constraints from LHC and dark matter direct detection experiments and predicting SM-like Higgs boson in agreement with current measurements.
Finally, we have used this opportunity to showcase our new interface between FlexibleSUSY and HiggsTools, allowing for a seamless validation of Higgs sectors of a broad class of user defined (supersymmetric and non-supersymmetric) models.

Figure 3 :Figure 4 :
Figure 3: Parameter regions around BMP7 allowed (green) and excluded (red) by SM-like Higgs data at 95% C.L. as reported by HiggsSignals.White regions is where no spectrum could be generated.

Figure 6 :
Figure 6: Parameter regions around BMP7 allowed (green) and excluded (red) by dark matter direct detection experiments at 95% C.L. White regions is where no spectra could be generated.Since in Fig. 6b all first and second generation soft squark masses were set equal, m Q;1,2 = m U ;1,2 = m D;1,2 ≡ m q , to facilitate scanning, none of the points corresponds directly to BMP7 or BMP8.

Figure 7 :
Figure 7:Regions allowed (green) and excluded (red) by direct searches of BSM particles at 95% C.L. as given by SModelS.White regions is where no spectrum could be generated.Since in Fig.7band 7c all first and second generation soft squark masses were set equal,m Q;1,2 = m U ;1,2 = m D;1,2 ≡ m q , tofacilitate scanning, none of the points corresponds directly to BMP7 or BMP8.

Table 2 :
Benchmark points for the scenario discussed here: input parameters, parameters determined via tadpole equations and selected predicted, phenomenologically relevant, pole masses.Dimensionfull parameters are given in GeV or GeV 2 , as appropriate.Input values are listed in the upper part of the table, while derived masses of some light physical states are in the lower part.FlexibleSUSY input/output files (including more significant digits of input parameters) are attached to the arXiv version of this work.

Table 3 :
Widths and branching ratios (in %) of the 95 and 125 GeV Higgs states for benchmark points BMP7 and 8 (see Tab. 2).More complete spectrum files are attached to the arXiv version of this work.

Table 4 :
Summary of experimental constraints passed by BMPs.