Implications of the Higgs discovery for the MSSM

The implications of the discovery of the Higgs boson at the LHC with a mass of approximately 125 GeV are summarised in the context of the minimal supersymmetric extension of the Standard Model, the MSSM. Discussed are the implications from the measured mass and production/decay rates of the observed particle and from the constraints in the search for the heavier Higgs states at the LHC.


Introduction
The historical discovery by ATLAS and CMS of a particle with a mass of approximately 125 GeV [1,2] and properties that are compatible with those of a scalar Higgs boson [3][4][5][6][7][8] has far reaching consequences not only for the Standard Model (SM) of the electroweak and strong interactions, but also for new physics models beyond it. This is particularly true for supersymmetric theories (SUSY) [9][10][11] that are widely considered to be the most attractive extensions of the SM as they naturally protect the Higgs mass against large radiative corrections and stabilise the hierarchy between the electroweak and Planck scales, besides of allowing for the unification of the three gauge coupling constants and providing a good candidate for the dark matter in the universe, the lightest SUSY particle.
In the minimal supersymmetric extension of the SM (MSSM), two Higgs doublet fields H u and H d are required to break the electroweak symmetry, leading to a physical spectrum with five Higgs particles: two CP-even, h and H , a CP-odd, A, and two charged, H ± , states [6,8]. Two parameters are needed to describe the MSSM Higgs sector at the tree level: one Higgs mass, which is generally taken to be that of the pseudoscalar boson M A , and the ratio of vacuum expectation values of the two Higgs fields, tan β = v d /v u , expected to lie in the range 1 tan β 60. The masses of the CP-even h, H and the charged H ± states, as well as the a e-mail: djouadi@th.u-psud.fr mixing angle α in the CP-even sector are uniquely defined in terms of these two inputs at tree level, but this nice property is spoiled at higher orders [12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28].
At high M A values, M A M Z , one is in the so-called decoupling regime [29] in which the neutral CP-even state h is light and has almost exactly the properties of the SM Higgs boson, i.e. its couplings to fermions and gauge bosons are the same as the standard Higgs, while the other CP-even H and the charged H ± bosons become heavy and mass degenerate with the A state, M H ≈ M H ± ≈ M A , and they decouple from the massive gauge bosons. In this regime, the MSSM Higgs sector thus looks almost exactly like the one of the SM with its unique Higgs boson.
There is, however, one major difference between the two cases: while in the SM the Higgs mass is essentially a free parameter (and should simply be smaller than about 1 TeV in order to ensure unitarity in the high-energy scattering of massive gauge bosons), the lightest MSSM CP-even Higgs particle mass is bounded from above and, depending on the SUSY parameters that enter the important quantum corrections, is restricted to M max h ≈ 90-130 GeV. The lower value comes from experimental constraints, in particular Higgs searches at LEP [30,31], while the upper bound assumes a SUSYbreaking scale that is not too high, M S O (1 TeV), in order to avoid too much fine-tuning in the model. Hence, the requirement that the MSSM h boson coincides with the one observed at the LHC, i.e. with M h ≈ 125 GeV and almost SM-like couplings as the LHC data seem to indicate, would place very strong constraints on the MSSM parameters, in particular the SUSY scale M S , through their contributions to the radiative corrections to the Higgs sector. This comes in addition to the limits that have been obtained from the search of the heavier Higgs states at the LHC, as well as from the negative search for supersymmetric particles.
In this review, we summarise the implications of the available LHC Higgs results for the MSSM Higgs sector. We first discuss the consequences of the M h measured value for the various unconstrained (with the many free parameters defined at the weak scale) and constrained (with parameters obeying some universal boundary conditions at the high scale) versions of the MSSM. We then discuss the impact of the measured production and decay rates of the observed particle on the various Higgs couplings and, hence, the MSSM parameters. The impact of the negative search of the heavy H, A and H ± states is summarised. An outlook is given in a concluding section. (1) where we use the short-hand notation s β ≡ sin β etc. and introduce the radiative corrections by a general 2 × 2 matrix ΔM 2 i j . One can then easily derive the neutral CP even Higgs boson masses and the mixing angle α that diagonalises the h and H states, H = cos α H 0 d +sin α H 0 u and h = − sin α H 0 d + cos α H 0 u : The by far leading one-loop radiative corrections to the mass matrix of Eq. (1) are controlled by the top Yukawa coupling, λ t = m t /v sin β with v = 246 GeV, which appears with the fourth power. One obtains a very simple analytical expression for the radiative correction matrix ΔM 2 i j if only this contribution is taken into account [12][13][14] ΔM 2 11 ∼ ΔM 2 12 ∼ 0 , where M S is the geometric average of the two stop masses M S = √ mt 1 mt 2 defined to be the SUSY-breaking scale and X t is the stop mixing parameter given by X t = A t −μ/ tan β with μ the higgsino mass parameter;m t is the running MS top quark mass to account for the leading two-loop QCD corrections in a renormalisation-group improved approach (some refinements can be include as well).
Other soft SUSY-breaking parameters, in particular μ and A b (and in general the corrections controlled by the bottom Yukawa coupling λ b = m b /v cos β, which at large value of μ tan β become relevant) as well as the gaugino mass parameters M 1,2,3 , provide a small but non-negligible correction to ΔM 2 i j and can thus also have an impact on the loop corrections [15,16,[25][26][27][28].
The maximal value of the h mass, M max h is given in the leading one-loop approximation above by (6) and is obtained for the choice of parameters [25][26][27][28]: -a decoupling regime with heavy A states, M A ∼ O(TeV); -large values of the parameter tan β, tan β 10; -heavy stops, i.e. large M S values and we choose in general M S ≤ 3 TeV to avoid a too large fine-tuning [32,33]; -a stop trilinear coupling X t = √ 6M S , the so-called maximal mixing scenario that maximises the stop loops [34].
If the parameters are optimised as above, the maximal h mass value can reach the level of M max h ≈ 130 GeV. An important aspect is that in the decoupling regime M A M Z , the heavier CP-even and the charged Higgs states become almost degenerate in mass with the CP-odd state, M H ≈ M H ± ≈ M A , while the mixing angle α becomes close to α ≈ π 2 − β making the couplings of the light h state to fermions and massive gauge bosons SM-like, and decoupling the H, H ± from the weak bosons as is the case for the state A by virtue of CP invariance.
In this section, we discuss the implications of the measured mass value of the observed Higgs boson at the LHC  that we identify with the lightest state h of the MSSM. We consider the phenomenological MSSM [63] in which the relevant soft SUSY parameters are allowed to vary freely (but with some restrictions) and constrained MSSM scenarios such as the minimal supergravity (mSUGRA) [64][65][66][67], gauge mediated (GMSB) [68][69][70][71][72] and anomaly mediated (AMSB) [73][74][75] supersymmetry-breaking models (for a review, see again Ref. [8]). We also discuss the implications of such an M h value for scenarios in which the supersymmetric spec- < M h < 127 GeV in the [M S , X t ] plane for some selected range of tan β values (right); from Ref. [35] trum is extremely heavy, the so-called split SUSY [76][77][78] or high-scale SUSY models [79,80]. Finally, a new parametrisation of the Higgs sector, which uses the crucial information M h = 125 GeV, is discussed [81].

Implications for the phenomenological MSSM
In an unconstrained MSSM, there is a large number of soft SUSY-breaking parameters, O(100), but analyses can be performed in the so-called "phenomenological MSSM" (pMSSM) [63], in which CP conservation, flavour diagonal sfermion mass and coupling matrices and universality of the first and second sfermion generations are imposed. The pMSSM involves then 22 free parameters in addition to those of the SM: besides tan β and M A , these are the higgsino mass μ, the three gaugino masses M 1,2,3 , the diagonal leftand right-handed sfermion mass parameters mf L ,R and the trilinear sfermion couplings A f .
As discussed above, an estimate of the upper bound on M h can be obtained by including the corrections that involve only the parameters M S and X t . However, to be more precise, one could scan the full pMSSM 22 parameter space in order to include the subleading corrections. To do so, one can use RGE programs such as Suspect [82] which calculate the Higgs and superparticle spectrum in the MSSM including the most up-to-date information [25].
To obtain the value M max h with the full radiative corrections, a large scan of the pMSSM parameters in an uncorrelated way was performed [35,36] in the domains: The results are shown in Fig. 1 What are the implications for the mass of the lightest stop statet 1 ? This is illustrated in the right-hand side of Fig. 1, where the contours are shown in the [M S , X t ] plane in which one obtains 123 < M h < 127 GeV from the pMSSM scan; the regions in which tan β 3, 5 and 60 are highlighted. One sees again that a large part of the parameter space is excluded if the Higgs mass constraint is imposed. In particular, large M S values, in general corresponding to large mt 1 are favoured. However, as M S = √ mt 1 mt 2 , the possibility that mt 1 is of the order of a few 100 GeV is still allowed, provided that stop mixing (leading to a significant mt 1 , mt 2 splitting) is large [36,[57][58][59]. Masses above 1 TeV for the scalar partners of light quarks and for the gluinos are also required by the direct searches of SUSY particles at the LHC [86,87], confirming the need of high M S values. Nevertheless, relatively light stops as well as electroweak sparticles such as sleptons, charginos and neutralinos are still possible allowing for a "natural SUSY" [33] in spite of the value M h ≈ 125 GeV. Nevertheless, the present LHC SUSY searches [86,87] are constraining more and more this natural scenario.

Implications for constrained MSSM scenarios
In constrained MSSM scenarios (cMSSM), the various soft SUSY-breaking parameters obey a number of universal boundary conditions at a high energy scale, thus reducing the number of basic input parameters to a handful. The various soft SUSY-breaking parameters are evolved via the MSSM renormalisation group equations down to the low energy scale M S where the conditions of proper electroweak symmetry breaking (EWSB) are imposed.
Three classes of such models have been widely discussed in the literature. There is first the minimal supergravity (mSUGRA) model [64][65][66][67] in which SUSY-breaking is assumed to occur in a hidden sector which communicates with the visible sector only via flavour-blind gravitational interactions, leading to universal soft breaking terms, namely a common m 1/2 , m 0 , A 0 values for the gaugino masses, sfermion masses and sfermion trilinear couplings. Then come the gauge mediated [68][69][70][71][72] and anomaly mediated [73][74][75] SUSY-breaking (GMSB and AMSB) scenarios in which SUSY-breaking is communicated to the visible sector via, respectively, gauge interactions and a super-Weyl anomaly.
These models are described by tan β, the sign of μ and a few continuous parameters. Besides of allowing for both signs of μ, requiring 1 ≤ tan β ≤ 60 and, to avoid excessive fine-tuning in the EWSB conditions, imposing the bound M S = M EWSB < 3TeV, we adopt the following ranges for the input parameters of these scenarios: Hence, in contrast to the pMSSM, the various parameters which enter the radiative corrections to M h are not all independent in these constrained scenarios, as a consequence of the relations between SUSY-breaking parameters that are set at the high-energy scale and the requirement that electroweak symmetry breaking is triggered radiatively for each set of input parameters. The additional constraints make that it is not possible to freely tune the parameters that enter the Higgs sector to obtain the pMSSM maximal value of M h . In order to obtain even a rough determination of M max h in a given constrained SUSY scenario, it is necessary to scan through the allowed range of values for the basic input parameters.
Using again the program Suspect, a full scan of these scenarios has been performed in Ref. [35] and the results for M max h are shown in the left-hand side of Fig. 2 as a function of tan β, the input parameter that is common to all models, and in the right-hand side of the figure as a function of M S . In the adopted parameter space of the models and with the central values of the SM inputs, the obtained upper h mass value is M max h ≈121 GeV in the AMSB scenario, i.e. much less that 125 GeV, while in the GMSB scenario one has M max h ≈122 GeV (these values are obtained for tan β ≈ 20). Thus, clearly, these two scenarios are disfavoured if the lightest h particle has indeed a mass in the range 123-127 GeV and M S 3 TeV. In mSUGRA, one obtains M max h = 128 GeV and, thus, some parameter space would still survive the M h constraint.
The upper bound on M h in these scenarios can be qualitatively understood by considering in each model the allowed values of the trilinear coupling A t , which essentially determines the stop mixing parameter X t and thus the value of M h for a given scale M S . In GMSB, one has A t ≈ 0 at relatively high scales and its magnitude does not significantly increase in the evolution down to the scale M S ; this implies that we are almost in the no-mixing scenario which gives a low value of M max h as can be seen from Fig. 1. In AMSB, one has a nonzero A t that is fully predicted at any renormalisation scale in terms of the Yukawa and gauge couplings; however, the ratio A t /M S with M S determined from the overall SUSY breaking scale m 3/2 turns out to be rather small, implying again that we are close to the no-mixing scenario. Finally, in the mSUGRA model, since we have allowed A t to vary in a wide range as |A 0 | ≤ 9 TeV, one can get a large A t /M S ratio, which leads to a heavier Higgs particle. However, one cannot easily reach A t values such that X t /M S ≈ √ 6 so that we are not in the maximal-mixing scenario and the higher upper bound on M h in the pMSSM cannot be reached.
In the case of mSUGRA, one can study several interesting special cases: the no-scale scenario with m 0 ≈ A 0 ≈ 0 [88,89], the scenario m 0 ≈ 0 and A 0 ≈ − 1 4 m 1/2 , which approximately corresponds to the constrained nextto-MSSM (cNMSSM) [90,91], A 0 ≈ −m 0 , which corresponds to a very constrained MSSM (VCMSSM) [92], and a non-universal Higgs mass model (NUHM) [93] in which the soft SUSY-breaking scalar mass terms are different for the sfermions and for the two Higgs doublet fields.
In two particular cases, namely the "no-scale" and the "approximate cNMSSM" scenarios, the upper bound on M h is much lower than in the more general mSUGRA case and, in fact, barely reaches M h ≈ 123 GeV. The main reason is that these scenarios involve small values of A 0 at the GUT scale, A 0 ≈ 0 for no-scale and A 0 ≈ − 1 4 m 1/2 for the cNMSSM, which lead to A t values at the weak scale that are too low to generate a significant stop mixing and, hence, one is again close to the no-mixing scenario. Thus, only a very small fraction of the parameter space of these two sub-classes of the mSUGRA model survive if we impose 123 < M h < 127 GeV. These models should thus have a very heavy sfermion spectrum as a value M S 3 TeV is required to increase M max h . In the VCMSSM case, the value M h 125 GeV can be reached as |A 0 | can be large for large m 0 , A 0 ≈ −m 0 , allowing for typical mixing.
Finally, since the NUHM is more general than mSUGRA as we have two more free parameters, the [tan β, M h ] area shown in Fig. 2 is larger than in mSUGRA. However, since we are in the decoupling regime and the value of M A does not matter much (as long as it a larger than a few hundred GeV) and the key weak-scale parameters entering the determination of M h , i.e. tan β, M S and A t are approximately the same in both models, one obtains a bound M max h that is only slightly higher in NUHM compared to the mSUGRA case.
In these constrained scenarios and, in particular in the general mSUGRA model, most of the scanned points giving the appropriate Higgs mass correspond to the decoupling regime of the MSSM Higgs sector and, hence, to an h boson with a SM-Higgs cross section and branching ratios. Furthermore, as the resulting SUSY spectrum for M h = 125 ± 2 GeV is rather heavy in these scenarios (easily evading the LHC limits from direct sparticle searches [86]), one obtains very small contributions to observables like the anomalous muon magnetic moment (g − 2) μ and to B-physics observables such as the rates BR(B s → μ + μ − ) or BR(b → sγ ) [94]. Hence, the resulting spectrum complies with all currently available constraints. In addition, as will be discussed later, the correct cosmological density for the LSP neutralino required by recent measurements [95] can easily be satisfied. The M h value provides thus a unique constraint in this decoupling regime.

Split and high-scale SUSY models
In the preceding discussion, we have always assumed that the SUSY-breaking scale is relatively low, M S 3 TeV, which implies a natural SUSY scenario [33] with supersymmetric and heavier Higgs particles that could be observed at the LHC. However, as already mentioned, this choice is mainly dictated by fine-tuning considerations which are a rather subjective matter as there is no compelling criterion to quantify the acceptable amount of tuning. One could well abandon the SUSY solution to the hierarchy problem and have a very high M S , which implies that, except for the lightest h boson, no other scalar particle is accessible at the LHC or at any foreseen collider. This argument has been advocated to construct the socalled split SUSY scenario [76][77][78] in which the soft SUSYbreaking mass terms for all the scalars of the theory, except for one Higgs doublet, are extremely large, i.e. their common value M S is such that M S 1 TeV (such a situation occurs e.g. in some string motivated models [96]). Instead, the mass parameters for the spin-1 2 particles, the gauginos and the higgsinos, are left in the vicinity of the EWSB scale, allowing for a solution to the dark matter problem and a successful gauge coupling unification, the two other SUSY virtues. The split SUSY models are much more predictive than the usual pMSSM as only a handful parameters are needed to describe the low-energy theory. Besides the common value M S of the soft SUSY-breaking sfermion and one Higgs mass parameters, the basic inputs are essentially the three gaugino masses M 1,2,3 (which can be unified to a common value at M GUT as in mSUGRA), the higgsino parameter μ and tan β. The trilinear couplings A f , which are expected to have values  [35] close to the EWSB scale set by the gaugino/higgsino masses that are much smaller than M S , will play a negligible role.
Concerning the Higgs sector, the main feature of split SUSY is that at the high scale M S , the boundary condition on the quartic Higgs coupling is determined by SUSY: where g and g are the SU(2) and U(1) gauge couplings. Here, tan β is not a parameter of the low-energy effective theory as it enters only the boundary condition above and cannot be interpreted as the ratio of the two Higgs vevs. If the scalars are very heavy, they will lead to radiative corrections in the Higgs sector that are significantly enhanced by large logarithms, log(M S /M EWSB ) where M EWSB ≈ |μ|, M 2 . In order to have reliable predictions, one has to properly decouple the heavy states from the low-energy theory and resum the large logarithmic corrections; in addition, the radiative corrections due to the gauginos and the higgsinos have to be implemented. Following the early work of Ref. [76][77][78], a comprehensive study of the split SUSY spectrum has been performed in Ref. [97]. All the features of the model have been implemented in the code SuSpect [82] upon which the analysis presented in Ref. [35] and summarised here is based.
One can adopt an even more radical attitude than in split SUSY and assume that the gauginos and higgsinos are also very heavy, with a mass close to the scale M S ; this is the case in the so-called high-scale SUSY model [79,80]. Here, one abandons not only the SUSY solution to the fine-tuning problem but also the solution to the dark matter problem by means of the LSP and the successful unification of the gauge couplings. However, there will still be a trace of SUSY at low energy: the matching of the SUSY and low-energy theories is indeed encoded in the Higgs quartic coupling λ of Eq. (8).
Hence, even if broken at very high scales, SUSY would still lead to a "light" Higgs whose mass will give information on M S and tan β.
The treatment of the Higgs sector of the high-scale SUSY scenario is similar to that of split SUSY: one simply needs to decouple the gauginos and higgsinos from the low-energy spectrum (in particular remove their contributions to the renormalisation group evolution of the gauge and Yukawa couplings and to the radiative corrections to M h ) and set their masses to M S . The version of the program Suspect which handles the split SUSY case can be adapted to also cover the Using this tool, a scan in the [tan β, M S ] plane has been performed to determine the value of M h in the split SUSY and high-scale SUSY scenarios; in the former case, M EWSB ≈ √ |M 2 μ| ≈ 246 GeV was chosen for the low scale. The results are shown in Fig. 3, where M h is displayed as a function of M S for selected values of tan β in both split (left plot) and high-scale (right plot) SUSY.
As expected, the maximal M h values are obtained at high tan β and M S values and, at the scale M S ≈ 10 16 GeV at which the couplings g and g approximately unify in the split SUSY scenario, one obtains M h ≈ 160 GeV for the higher tan β = 50 value. Not included is the error bands in the SM inputs that would lead to an uncertainty of about 2 GeV on M h , which is now mainly due to the 1 GeV uncertainty on m t . In addition, the zero-mixing scenario was assumed as the parameter A t is expected to be much smaller than M S ; this approximation might not be valid for M S values below 10 TeV and a maximal mixing A t /M S = √ 6 would increase the Higgs mass value by up to 10 GeV at M S = O(1 TeV) as was discussed earlier for the pMSSM. In the high-scale SUSY scenario, one obtains a value M h ≈ 142 GeV (with again an uncertainty of approximately 2 GeV from the top mass) for high tan β values and at the unification scale M S ≈ 10 14 GeV [79,80]. Much smaller M h values, in the 120 GeV range, can be obtained for lower scales and tan β.
Hence, the requirement that the Higgs mass is in the range 123 M h 127 GeV imposes strong constraints on the parameters of these two models. For this mass range, very large scales are needed for tan β ≈ 1 in the high-scale SUSY scenario, while scales not too far from M S ≈ 10 4 GeV are required at tan β 1 in both the split and the high-scale scenarios. In this case, SUSY should manifest itself at scales much below M GUT if M h ≈ 125 GeV.

Splitting the Higgs and sfermion sectors
In the previous high-scale scenarios, the Higgs mass parameters were assumed to be related to the mass scale of the scalar fermions in such a way that the masses of the heavier Higgs particles are also of the order of the SUSY scale, M A ≈ M S . However, this needs not to be true in general and one can, for instance, have a NUHM-like scenario where the Higgs masses are decoupled from those of the sfermions. If one is primarily concerned with the MSSM Higgs sector, one may be rather conservative and allow any value for M A irrespective of the SUSY-breaking scale M S . This is the quite "model-independent" approach that has been advocated in Refs. [98][99][100][101]: take M A as a free parameter of the pMSSM, with values ranging from O(100 GeV) up to O(M S ), but make no restriction on M S , which can be set to any value, even very high.
An important consequence of this possibility is that it reopens the low tan β region, tan β 3, which was long thought to be forbidden if one requires a SUSY scale M S 1 TeV, as a result of the limit M h 114 GeV from the negative search of a SM-like Higgs boson at LEP [31]. If the SUSY scale is large enough, these small tan β values would become viable again. To estimate the required magnitude of M S , one can still use Suspect in which the possibility M S 1 TeV is implemented [97] with the full set of radiative corrections up to two loops included. In Fig. 4, displayed are the contours in the plane [tan β, M S ] for fixed mass values M h = 120-132 GeV of the observed Higgs state (these include a 3 GeV theoretical uncertainty and also a 3 GeV uncertainty on the top quark mass [84] that is conservatively added linearly in the extreme cases). The maximal mixing X t = √ 6M S scenario is assumed with 1 TeV gaugino/higgsino mass parameters.
One observes that values of tan β ≈ 1 are possible and allow for an acceptable M h value provided the scale M S is large enough. For instance, while one can accommodate a scale M S ≈ 1 TeV with tan β ≈ 5, a large scale M S ≈ 20 TeV is required to obtain tan β ≈ 2; to reach the limit tan β = 1, an order of magnitude increase of M S will be required. Outside the decoupling regime, the obtained M S for a given M h value will be of course larger. For completeness, also shown is the contour for the LEP2 limit M h = 114 GeV which illustrates the fact that tan β ≈ 1 is still allowed provided that M S 20 TeV.

A new parametrisation of the Higgs sector
It was pointed out in Refs. [98,[102][103][104] that when the measured value of the h boson mass M h = 125 GeV is taken into account, the MSSM Higgs sector with solely the dominant radiative corrections included can be again described with only two free parameters such as tan β and M A as was the case at tree level. In other words, the dominant radiative corrections that involve the SUSY parameters are fixed by the value of M h . This observation leads to a rather simple parametrisation of the MSSM Higgs sector.
More specifically, let us assume that in the 2×2 matrix for the radiative corrections to the CP-even Higgs mass matrix Eq. (1), only the leading ΔM 2 22 entry of Eq. (5) that involves the by far dominant stop-top sector contribution is taken into account; this is the so-called approximation and its refinements [15,16,28]. In this ΔM 2 22 ΔM 2 11 , ΔM 2 12 limit, one can simply trade ΔM 2 22 for the by now known h mass value M h = 125 GeV and obtain This was called the habemus MSSM or hMSSM in Ref. [81]. However, this interesting and simplifying feature has to been demonstrated for all MSSM parameters and, in particular, one needs to prove that the impact of the subleading  2 11 and ΔM 2 12 . M A = 300 GeV, tan β = 5 or and 30 are taken and the other parameters are varied as described in the text [81] corrections ΔM 2 11 and ΔM 2 12 is small. To do so, a scan of the pMSSM parameter space using the program SuSpect, in which the full two-loop radiative corrections to the Higgs sector are implemented, has been performed [81]. For a chosen (tan β,M A ) input set, the soft-SUSY parameters that play an important role in the Higgs sector are varied in the following ranges: |μ| ≤ 3 TeV, |A t , A b | ≤ 3M S , 1 TeV ≤ M 3 ≤ 3 TeV and 0.5 TeV ≤ M S ≤ 3 TeV (≈3 TeV is the scale up to which programs such as SuSpect are expected to be reliable). The usual GUT relation between the weak scale gaugino masses 6M 1 = 3M 2 = M 3 has been assumed and A u , A d , A τ = 0 has been set (these last parameters have little impact on the radiative corrections). The MSSM Higgs sector parameters have been computed all across the parameter space, selecting the points which satisfy the constraint 123 ≤ M h ≤ 129 GeV when uncertainties are included. For each of theses points, the Higgs parameters have been compared to those obtained in the simplified MSSM approximation, ΔM 2 11 = ΔM 2 12 = 0, with the lightest Higgs boson mass as input. While the requirement that M h should lie in the range 123-129 GeV has been made, M h was allowed to be different from the one obtained in the "exact" case ΔM 2 11 , ΔM 2 12 = 0. Displayed in Fig. 5 are the differences between the values of the mass M H and the mixing angle α that are obtained when the two possibilities ΔM 2 11 = ΔM 2 12 = 0 and ΔM 2 11 , ΔM 2 12 = 0 are considered. This is shown in the plane [M S , X t ] with X t = A t − μ cot β when all other parameters are scanned as above. The A boson mass was fixed to M A = 300 GeV (a similar result was obtained for M A ≈ 1 TeV) and two representative values tan β = 5 and 30 are used. The conservative approach of plotting only points which maximise these differences has been adopted.
In all cases, the difference between the two M H values is very small (in fact, much smaller than the H boson total decay width Γ H ), less than a few percent, while for α the difference does not exceed ≈0.025 for low values of tan β, but at high tan β values, one can reach the level of ≈0.05 in some rare situations (large values of μ, which enhance the μ tan β contributions). Nevertheless, at high enough tan β, we are far in the decoupling regime already for M A 200 GeV and such a difference does not significantly affect the couplings of the h and H bosons which, phenomenologically, are the main ingredients.
Hence, even when including the full set of radiative corrections, it remains a good approximation to use Eq. (9) to derive the parameters M H and α in terms of the inputs tan β, M A and the measured M h value.
In the case of the charged Higgs boson (whose physics is described by tan β, M H ± and eventually α), the radiative Fig. 6 The SM-like Higgs boson branching ratios in the mass range 120-130 GeV (left) and its production cross sections at proton colliders as a function of the c.m. energy (right) [107] corrections to M H ± are much smaller for large enough M A and one has, at the few percent level in most cases (which is again smaller than the total H ± decay width), In conclusion, this approximation allows to ignore the radiative corrections to the Higgs masses and their complicated dependence on the MSSM parameters and to use a simple formula to derive the other parameters of the Higgs sector, α, M H as well as M H ± . This considerably simplifies phenomenological analyses in the MSSM Higgs sector which up to now rely either on large scans of the parameter space (as in the previous subsections) or resort to benchmark scenarios in which most of the MSSM parameters are fixed (as is the case of Ref. [105] for instance).

Implications of the Higgs production rates
3.1 Light Higgs decay and production at the LHC In many respects, the Higgs particle was born under a very lucky star as the mass value of ≈125 GeV (although too high for a natural SUSY) allows to produce it at the LHC in many redundant channels and to detect it in a variety of decay modes. This allows detailed studies of the Higgs properties as will be discussed in this section.
We start by summarizing the production and decay at the LHC of a light SM-like Higgs particle, which should correspond to the lightest MSSM h boson in the decoupling regime. First, for M h ≈ 125 GeV, the Higgs mainly decays into bb pairs but the decays into W W * and Z Z * final states, before allowing the gauge bosons to decay leptonically W → ν and Z → ( = e, μ), are also significant. The h → τ + τ − channel (as well as the gg and cc decays that are not detectable at the LHC) is also of significance, while the clean loop induced h → γ γ mode can be easily detected albeit its small rates. The very rare h → Z γ and even h → μ + μ − channels should be accessible at the LHC but only with a much larger data sample. This is illustrated in the left-hand side of Fig. 6, where the decay branching fractions of a SM-like Higgs are displayed for the narrow mass range M h = 120-130 GeV On the other hand, many Higgs production processes have significant cross sections as is shown in the right-hand side of Fig. 6, where they are displayed at a proton collider at various past, present and foreseen center of mass energies for a 125 GeV SM-like Higgs boson; the MSTW parton densities [106] have been used.
While the by far dominant gluon fusion mechanism gg → h (ggF) has extremely large rates (≈20 pb at √ s = 7-8 TeV), the subleading channels, i.e. the vector boson fusion (VBF) qq → hqq and the nHiggs-strahlung (HV) qq → hV with V = W, Z mechanisms, have cross sections which should allow for a study of the Higgs particle already at √ s 8 TeV with the amount of integrated luminosity, ≈25 fb −1 , which has been collected by each experiment. The Higgstop associated process pp → tth (ttH) would require higher energy and luminosity.
This pattern already allows the ATLAS and CMS experiments to observe the Higgs boson in several channels and to measure some its couplings in a reasonably accurate way. The channels that have been searched are h → Z Z * → 4 ± , h → W W * → 2 2ν, h → γ γ where the Higgs is mainly produced in ggF with subleading contributions from h j j in the VBF process, h → τ τ where the Higgs is produced in association with one (in ggF) and two (in VBF) jets, and finally h → bb with the Higgs produced in the HV process. One can ignore for the moment the additional search channels h → μμ and h → Z γ for which the sensitivity is still too low with the data collected so far. A convenient way to scrutinise the couplings of the produced h boson is to consider their deviation from the SM expectation. One then considers for a given search channel the signal strength modifier μ, which, with some approximation, can be identified with the Higgs production cross section times decay branching fractions normalised to the SM value. For the h → X X decay channel, one would have in the narrow width approximation, which, from the experimental side would correspond to where N ev X X is the measured number of events in the X X channel, the selection efficiency and L the luminosity.
ATLAS and CMS have provided the signal strengths for the various final states with a luminosity of, respectively, ≈5 fb −1 for the 2011 run at √ s = 7 TeV and ≈20 fb −1 for the 2012 run at √ s = 8 TeV. The constraints given by the two collaborations are shown in Fig. 7.
When the various analysed Higgs search channels are combined, this leads to a global signal strength [108,109] ATLAS : μ tot = 1.30 ± 0.30 CMS : μ tot = 0.87 ± 0.23 (13) which shows a good agreement with the SM expectation.
In fact, when the ATLAS and CMS values are combined, one finds a global signal strength that is very close to unity, implying that the observed Higgs is rather SM-like. Hence, already with the rather limited statistics at hand, the accuracy of the measurements in Eq. (13) is reaching the 20 % level for the ATLAS and CMS collaborations. This is at the same time impressive and worrisome. Indeed, as mentioned earlier the main Higgs production channel is the top and bottom quark loop mediated gluon fusion mechanism and, at √ s = 7 or 8 TeV, the three other mechanisms contribute at a level below 15 % when their rates are added and before kinematical cuts are applied. The majority of the signal events presently observed at the LHC, in particular in the main search channels h → γ γ, h → Z Z * → 4 , h → W W * → 2 2ν and, to a lesser extent, h → τ τ , thus come from the ggF mechanism which is known to be affected by large theoretical uncertainties.
As a matter of fact, although the cross section σ (gg → h) is known up next-to-next-to-leading order (NNLO) in perturbative QCD (and at least at NLO for the electroweak interaction) [110], there is a significant residual scale dependence which points to the possibility that still higher order contributions beyond NNLO cannot be totally excluded. In addition, as the process is of O(α 2 s ) at LO and is initiated by gluons, there are sizable uncertainties due to the gluon parton distribution function (PDF) and the value of the coupling α s . A third source of theoretical uncertainties, the use of an effective field theory (EFT) approach to calculate the radiative corrections beyond the NLO approximation, should in principle also be considered [111,112]. In addition, large uncertainties arise when the gg → h cross section is broken into the jet categories h +0 j, h +1 j and h +2 j [113]. In total, the combined theoretical uncertainty has been estimated to be of order Δ th ≈ ±15 % by the LHC Higgs cross section working group [110] and it would increase up to Δ th ≈ ±20 % if the EFT uncertainty is also included 2 [112].
Hence, the theoretical uncertainty is already at the level of the accuracy of the cross section measured by the ATLAS and CMS collaborations, Eq. (13). Another drawback of the analyses is that they involve strong theoretical assumptions on the total Higgs width since some contributing decay channels not accessible at the LHC are assumed to be SM-like and possible invisible Higgs decays in scenarios beyond the SM are supposed not to occur.
In Ref. [107], following earlier work [114][115][116][117], it has been suggested to consider the decay ratios D X X defined as for a specific production process p = ggF, VBF, VH or all (for inclusive production) and for a given decay channel h → X X when the reference channel h → V V is used. In these ratios, the cross sections σ p ( pp → h) and hence, their significant theoretical uncertainties will cancel out, leaving out only the ratio of decay branching fractions and hence of partial decay widths. These can be obtained with the program HDECAY [118] which includes all higher order effects and are affected by much smaller uncertainties. Thus, the total decay width which includes contributions from channels not under control such as possible invisible Higgs decays, do not appear in the ratios D p X X . Some common experimental systematical uncertainties such as the one from the luminosity measurement and the small uncertainties in the Higgs decay branching ratios also cancel out. We are thus, in principle, left with only with the statistical uncertainty and some (non common) systematical errors. The ratios D X X involve, up to kinematical factors and known radiative corrections, only the ratios |c X | 2 / |c V | 2 of the Higgs reduced couplings to the particles X and V compared to the SM expectation, c X ≡ g h X X /g SM h X X . For the time being, three independent ratios can be considered: D γ γ , D τ τ and D bb . D γ γ is the ratio of the inclusive ATLAS and CMS di-photon and Z Z channels that are largely dominated by the ggF mechanism; D τ τ is the signal strength ratio in the τ τ and W W searches where one selects Higgs production in ggF with an associated jet or in the VBF production mechanism; D bb is the ratio of the h → bb and h → W W decays in hV production for which the sensitivities are currently too low.
In order to test the compatibility of the couplings of the M h = 125 GeV Higgs state with the SM expectation, one can perform a fit based on the χ 2 R function The errors δ(μ X X /μ V V ) are computed assuming no correlations between the different final state searches. The uncertainties on the ratios are derived from the individual errors that are dominated by the experimental uncertainties as one expects that the theoretical uncertainties largely cancel out in the ratios D γ γ , D bb and D τ τ .
For the signal strengths above, the theoretical uncertainties have to be treated as a bias (and not as if they were associated with a statistical distribution) and the fit has to be performed for the two extremal values of the signal strengths: μ i | exp ±δμ i /μ i | th with the theoretical uncertainty δμ i /μ i | th conservatively assumed to be ±20 % for both the gluon and the vector boson fusion mechanisms (because of the contamination due to gg → h + 2 j in the latter case) and ≈5 % for hV associated production.

Fit of the Higgs couplings and their ratios
A large number of analyses of the Higgs couplings from the LHC data have been performed in the SM and its extensions and a partial list is given in Refs. .
In the MSSM, the couplings of the CP-even Higgs particles h and H to gauge bosons and fermions, compared to the SM Higgs couplings, are changed by factors that involve the sine and the cosine of the mixing angles β and α. Outside the decoupling regime where they reach unity, the reduced couplings (i.e. normalised to their SM values) of the lighter h state to third generation t, b, τ fermions and gauge bosons V = W/Z are for instance given by They thus depend not only on the two inputs [tan β, M A ], as occurs at tree level, but, a priori, on the entire MSSM spectrum as a result of the radiative corrections, in the same way as the Higgs masses. In principle, as discussed earlier, knowing tan β and M A and fixing M h to its measured value, the couplings can be determined in general. However, this is true when only the radiative corrections to the Higgs masses are included. Outside the regime in which the pseudoscalar A boson and the supersymmetric particles are very heavy, there are also direct radiative corrections to the Higgs couplings not contained in the mass matrix of Eq. (1) and which can alter this simple picture. First, in the case of b-quarks, additional one-loop vertex corrections modify the tree-level hbb coupling: they grow as m b μ tan β and can be very large at high tan β. The dominant component comes from the SUSY-QCD corrections with sbottom-gluino loops that can be approximated by ) [151]. Outside the decoupling regime the c b coupling reads with tan α → −1/ tan β for M A M Z . A large Δ b would significantly alter the dominant h → bb partial width and affect the branching fractions of all other decay modes.
In addition, the htt coupling is derived indirectly from the gg → h production cross section and the h → γ γ decay branching ratio, two processes that are generated by triangular loops. In the MSSM, these loops involve not only the top quark (and the W boson in the decay h → γ γ ) but also contributions from supersymmetric particles, if not too heavy. In the case of gg → h production, only the contributions of stops is generally important. Including the later and working in the limit M h m t , mt 1,2 , the coupling c t from the ggF process 3 is approximated by [152] which shows that indeed,t contributions can be very large for light stops and for large stop mixing. In the h → γ γ decay rate, because the t,t electric charges are the same, the htt coupling is shifted by the same amount. If one ignores the usually small contributions of the other sparticles (to be discussed in the next subsection), the htt vertex can be simply parametrised by the effective coupling of Eq. (20). We note that the h couplings to τ leptons and c quarks do not receive the direct corrections of Eqs. (19) and (20) and one should still have c c = c 0 t and c τ = c 0 b . However, using c t,b or c 0 t,b in this case has almost no impact in practice as these couplings appear only in the branching ratios for the decays h → cc and τ + τ − , which are small and the direct corrections should not be too large. One can thus, in a first approximation, assume that c c = c t and c τ = c b . Another caveat is due to the invisible Higgs decays which are assumed to be absent and which will be discussed later. Hence, because of the direct corrections, the Higgs couplings cannot be described only by β and α as in Eq. (18). To characterise the Higgs particle at the LHC, it was advocated that at least three independent h couplings should be considered, namely c t , c b and c V = c 0 V [81]. One can thus define the following effective Lagrangian:  and [c t , c b ]. As an example, the fit of the signal strengths and their ratios in the [c t , c b ] plane is displayed in Fig. 9. In this two-dimensional case, the best-fit point is located at c t = 0.88 and c b = 0.97, while c V 1. Note that although for the best-fit point one has c b 1, actually c b 1 in most of the 1σ region.
Using the formulae Eq. (9) for the angle α and using the input M h ≈ 125 GeV, one can make a fit in the plane [tan β, M A ]. This is shown in Fig. 10, where the 68, 95 and 99 % CL contours from the signal strengths and their ratios are displayed when the theory uncertainty is taken as a bias. The best-fit point when the latter uncertainty is set to zero is obtained for the values tan β = 1 and M A = 557 GeV, which implies for the other parameters using

An excess in the γ γ channel?
In the early LHC data, a significant excess in the h → γ γ detection channel was observed, raising the hope that it could be the first signal for physics beyond the SM. This excess has  [81] unfortunately faded away with more statistics and with the full 25 fb −1 data collected at √ s = 7+8 TeV, there is now only a ≈ 2σ excess in ATLAS which measures μ γ γ = 1.6 ± 0.3, while the signal strength measured by the CMS collaboration is μ γ γ = 0.9±0.3, which is SM-like. Nevertheless, it would be interesting to briefly discuss this excess as, besides the fact that it has triggered a vast literature, the h → γ γ channel is the one where new physics, and SUSY in particular, is most likely to manifest itself.
First, it has been realised early that this excess, if not due to a statistical fluctuation, can be easily explained or reduced in the context of the SM by invoking the large theoretical uncertainties that affect the production times decay rate in the dominant channel, gg → h → γ γ . This is shown in Fig. 11, where the ATLAS and CMS ratios R γ γ ≡ μ γ γ and their combination, obtained with the ≈ 10 fb −1 data collected at √ s = 7+8 TeV, is compared to the theory uncertainty bands obtained by the LHC Higgs group [110] and in Ref. [112]. It is clear that including the theory uncertainty as a bias helps to reduce the discrepancy between measurement and expectation and e.g. the excess reduces to 1.3σ from the original 2σ value. Ignoring this option, let us summarise the various possibilities that could explain this excess in the context of the MSSM. Deviations of μ γ γ from the SM value may be due to modifications of either the production cross section or the decay branching fraction or to both. The h decay branching fractions may be modified by a change of the h total decay width. Since the dominant decay mode is h → bb, a change of the effective hbb coupling by the direct vertex correc- Fig. 11 The value of μ γ γ given by the ATLAS and CMS collaborations with the ≈10 fb −1 data collected at √ s = 7 and 8 TeV, as well as their combination, compared to two estimates of the theoretical uncertainty bands; from Ref. [153] tions of Eq. (19) outside the decoupling regime can change all other Higgs rates including h → γ γ . The total width can also be modified by additional decay channels to SUSY particles and the only ones that are allowed by experimental constraints are invisible decays into the LSP that will be discussed later.
Nevertheless, these two possibilities would not only affect the h → γ γ rate but also those of other channels such as h → Z Z where no excess has been observed. It is thus more appropriate to look at deviation in the h → γ γ loop induced decay only. In the MSSM, the h → γ γ process receives contributions from scalar top and bottom quarks, staus and charginos as briefly is summarised below.
-Light stops: as already discussed, the M h = 125 GeV constraint requires large M S = √ mt 1 mt 2 and/or X t values. If mt 1 500 GeV, one should have maximal mixing X t ≈ √ 6M S and, in this case, the ht 1t1 coupling of Eq. (20) is large and leads to a sizeable change of the gg → h → γ γ rate; cf. Fig. 12 (left). However, an enhancement of the h → γ γ rate is over-compensated by a suppression of σ (gg → h) that seems not to occur. μ γ γ is enhanced only in the no-mixing case, X t ≈ 0, which requires very heavy stops which decouple from the amplitude [57- 59,152]. -Light sbottoms: a lightb R state does not conflict with the M h value as its corrections to the mass are small. For mb 1 500 GeV, it contributes to the hgg loop but it reduces the gg → h production rate; Fig. 12 (right). In turn, it has little impact on the h → γ γ rate because of the largely dominating W loop and the small b 1 electric charge. For mb 1 1 TeV, as indicated by direct LHC searches, μ γ γ is unaffected by sbottoms loops [152]. -Light staus: they lead to the largest contributions and have received most of the attention in the literature; see e.g.
Ref. [154][155][156]. For low mτ L/R values, a few 100 GeV, and large mixing X τ = A τ − μ tan β, with tan β ≈ 60 and |μ|=0.5-1 TeV, the lighter stau state has a mass close to the LEP2 bound, mτ 1 ≈ 100 GeV and its coupling to the h boson, g hττ ∝ m τ X τ , is huge. Theτ 1 contribution can hence significantly increase BR(h → γ γ ), up to 50 % [154][155][156], but this occurs only for extreme choices of the parameters. -Light charginos: the hχ + i χ − i couplings are in general small and are maximal when the χ ± i states are almost equal higgsino-wino mixtures. For a mass above 100 GeV and maximal couplings to the h boson, the χ ± 1 contributions to the h → γ γ rate do not exceed the 10-15 % level (with a sign being the same as the sign of μ) [157,158].   [149] Of course, different contributions can sum up resulting in more sizeable shifts. However, a 50 % deviation of the rate is unlikely and occurs only in extreme situations.

Invisible Higgs decays?
Invisible decays can also affect the properties of the observed h particle. In the MSSM, because of the LEP2 constraints, the only possible invisible channel for the h boson is into pairs of the LSP neutralinos, h → χ 0 1 χ 0 1 . BR inv can be important for mχ0 1 < 60GeV and for not too large M 1 and |μ| values which make the LSP a higgsino-gaugino mixture with significant couplings to the h state. Such an LSP would have the relic density Ωh 2 required by the WMAP results [95] since it will annihilate efficiently through the s-channel exchange of the h boson. However, BR inv should be small in this case. This is exemplified in Fig. 13, where log 10 (Ω χ h 2 ) is shown as a function of m χ 0 1 for the pMSSM points that satisfy the LHC Higgs constraints and BR(h → χ 0 1 χ 0 1 ) ≥ 15 %. Only a small area in the region 30 m χ 0 1 60 GeV fulfils these conditions.
The invisible Higgs decay width can be constrained indirectly by a fit of the Higgs couplings and in particular with the signal strength μ Z Z , which is the most accurate one and has the least theoretical ambiguities. Γ inv H enters in the signal strength through the total width Γ tot H calculated with free coefficients c f and c V . The resulting 1σ or 2σ ranges are shown in Fig. 14, where c f is freely varied while c V = 1 [150]. This gives Γ inv H /Γ SM H 50 % at the 95 % CL if the assumption c f = c V = 1 is made.
A more model independent approach would be to perform direct searches for missing transverse energy. These  [150]. The dependence on the theory uncertainties are shown by the black curves that indicate the other possible extreme domains. The direct upper limit on Γ inv H from direct searches at LHC for c V = c f = 1 [159] is also shown have been conducted by ATLAS [159] and CMS [163] in the pp → hV process with V → j j, and in the VBF channel, qq → qq E T / . As an example, we show in Fig. 15 (left) the CMS results for the Higgs cross section times BR inv versus M h when the two channels are combined. For M h ≈ 125 GeV a bound BR inv 50 % is obtained at the 95 % CL.
A more promising search for invisible decays is the monojet channel. In the ggF mode, an additional jet can be emitted at NLO leading to gg → h j final states and, because the QCD corrections are large, σ (H + 1 j) is not much smaller than σ (h + 0 j). The NNLO corrections besides significantly increasing the h + 0 j and h + 1 j rates, lead to h + 2 j events that also occur in VBF and VH. Hence, if the Higgs is coupled to invisible particles, it may recoil against hard QCD radiation, leading to monojet events.
In Refs. [160][161][162], it has been shown that the monojet signature carries a good potential to constrain the invisible decay width of a ≈ 125 GeV Higgs boson. In a model independent fashion, constraints can be placed on the rates Recent monojet searches made by CMS and ATLAS [164,165] are sensitive to R inv close to unity. This is shown in Fig  15 (right) where the best-fit region to the LHC Higgs data is displayed in the Br inv -c gg plane, where c gg is the deviation of σ (gg → h) from the SM expectation [160]. For the SM value c gg = 0, Br inv 20 % is disfavoured at 95 % CL while for c gg > 0, a larger rate is allowed, up to Br inv ∼ 50 %. best-fit regions to the combined LHC Higgs data. The black region is excluded by the monojet constraints while the red region is excluded by the ATLAS Z + E/ T search [159]; from Ref. [160] The Higgs invisible rate and the dark matter detection rate in direct astrophysical searches are correlated in Higgs-portal models. Considering the generic cases of scalar, fermionic and vectorial dark matter particles χ that couple only to the Higgs, one can translate in each case the LHC constraint BR(h → inv.) into a constraint on the Higgs couplings to the χ particles. It turns out that these constraints are competitive with those derived from the bounds on the dark matter scattering cross section on nucleons from the best experiment so far, XENON [95]. This is shown in Fig. 16, where the maximum allowed values of the scattering cross sections are given in the three cases assuming BR inv χ 20 %. The obtained spinindependent rates σ SI χ p are stronger than the direct limit from the XENON100 experiment in the entire M χ 1 2 M h range. In other words, the LHC is currently the most sensitive dark matter detection apparatus, at least in the context of simple Higgs-portal models.

Determination of the Higgs parity
Apart from the measurement of the couplings, one also needs in principle to establish that the observed Higgs state is indeed a CP even scalar particle and hence with J PC = 0 ++ quantum numbers 4 . It is known that the Higgs to vector boson (hV V ) coupling is a possible tool to probe these quantum numbers at the LHC [170,171]. This can be done by studying Fig. 16 Bounds on the spin-independent direct detection cross section σ SI χ p in Higgs-portal models derived for an invisible branching fraction of 20 % (colored lines) for a 125 GeV Higgs. These are compared to the current and future direct bounds from the XENON experiment (black lines). From Ref. [166] various kinematical distributions in the Higgs decay and production processes. One example is the threshold behaviour of the M Z * spectrum in the h → Z Z * → 4 decay channel and another is the azimuthal distribution between the decay planes of the two lepton pairs arising from the Z , Z * bosons from the Higgs decay. These are sensitive to both the spin and the parity of the Higgs.
With the 25 fb −1 data collected so far, the ATLAS and CMS collaborations performed a matrix-element likelihood analysis which exploits the kinematics and Lorenz structure of the h → Z Z * → 4 channel to see whether the angular distributions are more compatible with the 0 + or 0 − hypothesis (as well as the spin-2 possibility) [172,173]. Assuming that it has the same couplings as the SM Higgs boson and Fig. 17 Discrimination between the 0 + and 0 − parity hypotheses for the observed Higgs boson using the kinematics of the h → Z Z * → 4 channel by the ATLAS (left) and CMS (right) collaborations with the data collected at 7+8 TeV [172,173] that it is produced mainly from the dominant ggF process, the observed particle is found to be compatible with a 0 + state and the 0 − possibility is excluded at the 97.8% confidence level or higher; see Fig. 17.
Other useful diagnostics of the CP nature of the Higgs boson that also rely on the different tensorial structure of the hV V coupling can be made in the VBF process. It was known since a long time that in this channel, the distribution in the azimuthal angle between the two jets produced in association with the Higgs discriminates a CP-even from a CP-odd state [174]. This has been extended recently to other observables in this process, like the rapidity separation between the two jets [175][176][177][178].
Recently, the VBF channel pp → H j j has been reanalysed in the presence of an anomalous hV V vertex that parametrises different spin and CP assignments of the produced Higgs boson [178]. The anomalous hV V coupling is introduced by allowing for an effective Lagrangian with higher dimensional operators, which include four momentum terms which are absent in the SM. It was shown that the kinematics of the forward tagging jets in this process is highly sensitive to the structure of the anomalous coupling and that it can effectively discriminate between different assignments for the spin (spin-0 versus spin-2) and the parity (CP-even versus CP-odd) of the produced particle. In particular, it was found that the correlation between the separation in rapidity and the transverse momenta of the scattered quarks, in addition to the already discussed distribution of the azimuthal jet separation, can be significantly altered compared to the SM expectation. This is exemplified in Fig. 18, where the difference in rapidity between tagging jets (Δy j j ) for each of the higher dimensional operators in the hV V couplings is displayed.
These kinematical variables define new corners of the phase space that have not been explored by the experiments at the LHC to probe anomalous hV V couplings and to check the Higgs parity. In addition, some of these observables sig- Fig. 18 Normalised distribution of the difference in rapidity between the scattered jets in VBF for each of the SM and BSM operators (spin-2, CP-even and CP-odd) individually [178] nificantly depend on the c.m. energy and strong constraints on anomalous couplings can be obtained by performing measurements at the LHC with energies of √ s = 8 and 14 TeV. Finally, the associated hV production channel can be used as the invariant mass of the V h system as well as the p T and rapidities of the h and V bosons are also sensitive to anomalous hV V couplings.
Nevertheless, there is a caveat in the analyses relying on the hV V couplings. Since a CP-odd state has no treelevel V V couplings, all the previous processes project out only the CP-even component of the hV V coupling [179][180][181][182] even if the state is a CP-even and -odd mixture. Thus, in the CP studies above, one is simply verifying a posteriori that indeed the CP-even component is projected out.
A better way to measure the parity of the Higgs boson is to study the signal strength in the h → V V channels [150,183]. Indeed, the hV V coupling takes the general form g μν hV V = gives an upper bound on the CP-odd contribution to the hV V coupling.
Using μ V V and the ratios μ γ γ /μ V V and μ τ τ /μ V V as in Eq. (17), it was demonstrated that the particle has indeed a large CP component, 50 % at the 95 % CL, if the Higgs couplings to fermions are SM like. This is shown in Fig. 19, where one sees that the pure CP-odd possibility is excluded at the 3σ level, irrespective of the (mixed CP) Higgs couplings to fermions provided that |c f | 2 = 1.

Implications from heavy Higgs searches
We turn now to the constraints on the MSSM Higgs sector that can be obtained from the search of the heavier H/A and H ± states at the LHC and start with a brief summary of their production and decay properties.

H, A, H ± decays and production at the LHC
The production and decay pattern of the MSSM Higgs bosons crucially depend on tan β. In the decoupling regime that is indicated by the h properties, the heavier CP-even H boson has approximately the same mass as the A state and its interactions are similar. Hence, the MSSM Higgs spectrum will consist of a SM-like Higgs h ≡ H SM and two pseudoscalarlike particles, Φ = H/A. The H ± boson will also be mass degenerate with the Φ states and the intensity of its couplings to fermions will be similar. In the high tan β regime, the couplings of the non-SM like Higgs bosons to b quarks and to τ leptons are so strongly enhanced, and the couplings to top quarks and massive gauge bosons suppressed, that the pattern is rather simple. This is first the case for the decays: the Φ → tt channel and all other decay modes are suppressed to a level where their branching ratios are negligible and the Φ states decay almost exclusively into τ + τ − and bb pairs, with branching ratios of BR(Φ → τ τ ) ≈ 10 % and BR(Φ → bb) ≈ 90 %. The H ± boson decay into τ ν τ final states with a branching fraction of almost 100 % for H ± masses below the tb threshold, M H ± m t + m b , and a branching ratio of only ≈10 % for masses above this threshold while the rate for H ± → tb will be at the ≈90 % level in most cases.
Concerning the production, the strong enhancement of the b-quark couplings at high tan β makes that only two processes are relevant in this case: gg → Φ fusion with the b-loop included and associated production with b-quarks, gg/qq → bb + Φ, which is equivalent to the fusion process bb → Φ when no additional final b-quark is present. All other processes, in particular V Φ, ttΦ and VBF have suppressed rates. In both the bb and the gg fusion cases, as M Φ m b , chiral symmetry holds and the rates are approximately the same for the CP-even H and CP-odd A bosons. While σ (gg → Φ) is known up to NLO in QCD [184], σ (bb → Φ) is instead known up to NNLO [185,186].
The most powerful search channel for the heavier MSSM Higgs particles at the LHC is by far the process pp → gg + bb → Φ → τ + τ − . The precise values of the cross section times branching fraction for this process at the LHC have been updated in Refs. [110,112] and an assessment of the associated theoretical uncertainties has been made. It turns out that, in the production cross section, the total uncertainty from scale variation, the PDFs and α s as well as from the bquark mass are not that small: Δ TH σ ( pp → Φ) × BR(Φ → τ τ ) ≈ ±25 % in the entire M Φ range probed at the LHC at √ s = 8 TeV; Fig. 20.
Besides the QCD uncertainty, three other features could alter the rate σ ( pp → Φ → τ τ ) in the MSSM and they are related to the impact of the SUSY particle contributions: Thus, only in the unlikely cases where the decay H → τ 1τ1 has a branching rate of the order of 50 %, the squark loop contribution to the gg → H process is of the order 50 %, or the Δ b SUSY correction is larger than 100 %, and one can change the pp → Φ → τ τ rate by ≈25 %, which is the level of the QCD uncertainty. One thus expects σ ( pp → Φ) × BR(Φ → τ τ ) to be extremely robust and to depend almost exclusively on M A and tan β.
Finally, for the charged Higgs boson, the dominant search channel is in H ± → τ ν final states with the H ± bosons produced in top quark decays for M H ± m t − m b ≈ 170 GeV, pp → tt with t → H + b → τ νb. This is particularly true at high tan β values when BR(t → H + b) ∝m 2 b tan 2 β is significant. For higher H ± masses, one should rely on the three-body production process pp → tbH ± → tbτ ν, but the rates are presently rather small.
In the low tan β regime, tan β 5, the phenomenology of the heavier A, H, H ± bosons is richer [98,[187][188][189][190]. Starting with the cross sections, we display in Fig. 21 the rates for the relevant production processes at the LHC with √ s = 8 TeV assuming tan β = 2.5. For smaller tan β values, the rates except for pp → H/A + bb are even larger as the H/A + tt and H V V couplings are less suppressed. Fig. 21 The production cross sections of the MSSM heavier neutral Higgs bosons at the LHC at √ s = 8 for tan β = 2.5; only the main production channels are considered [98] Fig. 22 The H/A/H ± decay branching ratios as functions of the Higgs masses for tan β = 2.5 [98] Because of CP invariance which forbids AV V couplings, there is no AV and Aqq processes while the rates for associated tt A and bb A are small because the Att (Abb) couplings are suppressed (not sufficiently enhanced) compared to the SM Higgs. Only the gg → A process with the dominant t and sub-dominant b contributions included provides large rates. The situation is almost the same for H : only gg → H is significant at M H 300 GeV and tan β 5; the VBF and HV modes give additional small contributions for tan β ≈ 1. For H ± , the dominant production channel is again top quark decays, t → H + b for M H ± 170 GeV as for tan β 5, the m t / tan β piece of the H ± tb coupling becomes large; for higher H ± masses, the main process to be considered is gg/qq → H ± tb.
Turning to the H/A/H ± decay pattern, it can be rather involved at low tan β. A summary is as follows for tan β 3; see also Fig. 22 Hence, many decay and production channels need to be considered in this low tan β regime.

Constraints from the LHC Higgs searches
The most efficient channel to probe the heavier MSSM Higgs bosons is by far pp → gg + bb → H/A → τ + τ − . Searches for this process have been performed by ATLAS with ≈5 fb −1 data at the 7 TeV run [191] and by CMS with ≈5+12 fb −1 data at the 7 and 8 TeV runs [192]. Upper limits on the production cross section times decay branching ratio have been set and they can be turned into constraints on the MSSM parameter space.
In Fig. 23, the sensitivity is displayed of the CMS pp → Φ → τ τ analysis with 17 fb −1 of data in the [tan β, M A ] plane. The excluded region, obtained from the observed limit at the 95 % CL, is drawn in blue. The dotted line represents the median expected limit which turns out to be weaker than the observed limit. As can be seen, this constraint is extremely restrictive and, for values M A 250 GeV, it excludes almost the entire intermediate and high tan β regimes, tan β 5. The constraint is less effective for a heavier A boson, but even for M A ≈ 400 GeV the high tan β 10 region is excluded and there is even sensitivity to large values M A ≈ 800 GeV for tan β 50.
There are, however, some caveats to this exclusion limit as discussed previously. The first one is that there is a theoretical uncertainty of order of ±25 % that affects the gg → Φ and bb → Φ production cross sections which, when included, will make the constraint slightly weaker as one then needs to consider the lower value predicted for the production rate. A second caveat is that SUSY effects, direct corrections to the production and H/A decays into sparticles, could alter the rate. However, as previously argued, σ ( pp → Φ) × BR(Φ → τ τ ) is robust against these SUSY effects and the latter will unlikely make a substantial change of the cross section times branching fraction. Finally, the constraint is specifically given in the maximal mixing scenario X t /M S = √ 6 with M S = 1 TeV. The robustness of σ × BR makes that the exclusion limit is actually almost model independent and is valid in far more situations than the "MSSM M max h scenario" quoted there, an assumption that can be removed without any loss.
In fact, the exclusion limit can also be extended to the low tan β region which, in the chosen scenario with M S = 1 TeV, is excluded by the LEP2 limit on the lighter h mass (the green area in the figure) but should resurrect if the SUSY scale is kept as a free parameter. Note also that H/A bosons have also been searched for in the channel gg → bbΦ with Φ → bb (requiring more than three tagged b jets in the final state) but the constraints are much less severe than the ones derived from the τ τ channel [193].

Fig. 24
The H ± limits from ATLAS with √ s = 8 TeV and ≈20 fb −1 data in the channel t → bH + → bτ ν [194] Turning to the H + boson [194,195], the most recent result has been provided by the ATLAS collaboration using the full ≈20 fb −1 data collected at √ s = 8 TeV. The H ± search as been performed using the τ plus jets channel with a hadronically decaying τ lepton in the final state. For M H ± 160 GeV, the results are shown in Fig. 24. Here, the relevant process is top quark decays, t → H + b with the decay H + → τ ν having a branching ratio of almost 100 % at moderate to high tan β. For these high values, the H + tb coupling has a component ∝ m b tan β, which makes BR(t → ! H + b) rather large. Almost the entire tan β 10 region is excluded by the ATLAS analysis.
In addition, the branching fraction for the decay t → bH + is also significant at low tan β values, when the component of the coupling g tbH + ∝m t / tan β becomes dominant. On the other hand, the branching fraction for the decay H ± → τ ν does not become very small as it has competition only from H + → cs, which, even for tan β ≈ 1, does not dominate. Hence, the rates for pp → tt with t → bH + → bτ ν are comparable for tan β ≈ 3 and tan β ≈ 30 and the processes can also probe the low tan β region. This is exemplified in Fig. 24 where one can see that the entire area below tan β ≈ 5 is also excluded. There remains then, for H ± masses close to 90 GeV (where the detection efficiency is lower) and 160 GeV (where one is limited by the phase space), the intermediate area with tan β ≈ 5-10 where the H ± tb coupling is not strongly enhanced.
This ATLAS search has been extended to larger values of M H ± where the charged Higgs is produced in association with top quarks, gb → t H + , but the constraints are poor (only the region tan β 50 is excluded for M H ± = 200-300 GeV) as the cross section for this process is low.
The reopening of the low tan β region allows to consider a plethora of very interesting channels for the heavier Higgs bosons to be also investigated at the LHC: heavier CP-even H decays into massive gauge bosons H → W W, Z Z and Higgs bosons H → hh, CP-odd Higgs decays into a vector and a Higgs boson, A → h Z, CP-even and CP-odd Higgs decays into top quarks, H/A → tt, and even the charged Higgs decay H ± → W h. These final states have been searched for in the context of a heavy SM Higgs boson or for new resonances in some non-SUSY beyond the SM scenarios and the analyses can be adapted to the case of the heavier MSSM Higgs bosons. They would then allow to cover a larger part of the parameter space of the MSSM Higgs sector in a modelindependent way, i.e. without using the information on the scale M S and more generally on the SUSY particle spectrum that appear in the radiative corrections.
In Ref. [98] a preliminary analysis of these channels has been performed using current information given by the ATLAS and CMS collaborations in the context of searches for the SM Higgs boson or other heavy resonances (in particular new Z or Kaluza-Klein gauge bosons that decay into tt pairs). The results are shown in Fig. 25 with an extrapolation to the full 25 fb −1 data of the 7+8 TeV LHC run (it has been assumed that the sensitivity scales simply as the square root of the number of events). The sensitivities from the usual H/A → τ + τ − and t → bH + → bτ ν channels are also shown. The green and red areas correspond to the domains where the H → V V and H/A → tt channels become constraining. The sensitivities in the H → hh and A → h Z modes are given by, respectively, the yellow and brown areas which peak in the mass range M A = 250-350 GeV that is visible at low tan β values. The outcome is impressive. These channels, in particular the H → V V and H/A → tt processes, are very constraining as they cover the entire low tan β area that was previously excluded by the LEP2 bound up to M A ≈ 500 GeV. Even A → h Z and H → hh would be visible at the current LHC in small portions of the parameter space.

Could the observed state be the heavier H boson?
Let us briefly discuss the possibility, raised with the early LHC data, that the observed particle is the heavier MSSM H boson, as advocated for instance in Refs. [60][61][62]105]. The possibility M H ≈ 125 GeV, with H couplings close to those of the SM Higgs, occurs for low values of M A , ≈100-120 GeV, and moderate values of tan β, ≈ 10. In this case, H has approximately SM-like properties, while h has a mass of order 100 GeV or below and suppressed couplings to vector bosons. A dedicated scan for this region of parameter space has been performed [36] and the results were confronted with the measured Higgs mass M h = 123-129 GeV and couplings that comply with the LHC ≈ 10 fb −1 data collected at √ s = 7+8 TeV. Both the signal strengths in the various search channels of the observed Higgs boson and the limits from the pp → τ + τ − channel obtained by the CMS collaboration have been considered.
It was found that among the large flat scan with 10 8 points, only ≈ 2 × 10 −5 of the generated points would remain after imposing these LHC constraints. These points were then excluded by applying the constraints from flavour physics [94] (see also Ref. [196]), mainly the radiative decay b → sγ and dark matter constraints [95] (as they do not satisfy the constraint of 10 −4 < Ωh 2 < 0.155 when accounting for all uncertainties). The updated pp → τ + τ − search performed by CMS with 17 fb −1 data, which excludes all values tan β 5 for M A 250 GeV as shown in Fig. 23, now definitely rules out this scenario. This is exemplified in Fig. 26, where we zoom in the [M A , tan β] plane for low values of the inputs and apply the constraints listed above. The small region in which the H boson was allowed to be the observed state (black points) by the previous H/A → τ + τ − CMS search (dashed blue line) is excluded by the new data (in red). In fact, the latest ATLAS limits from H ± searches given in Fig. 24 also exclude now the possibility M A ≈ 100-120 GeV and, hence, the scenario where H is the observed Higgs state 5 .

Higgs production with SUSY particles
Finally, let us comment on the possibility of the Higgs bosons being produced in processes involving sparticles. First of all, there is the option of Higgs decays into SUSY particles. In the case of the lighter h boson, the only possibility when the LEP2 constraints are taken into account is the decay h → χ 0 1 χ 0 1 , which has been discussed in the context of invisible Higgs decays in Sect. 3.4. In view of the strong LHC limits on squark masses, the only SUSY channels of the heavier H/A/H ± states that might be kinematically open would be the decays into chargino, neutralinos and sleptons. For H/A, these decays have been discussed in the context of the τ searches as they might reduce the H/A → τ τ branching fractions but no specific search for these SUSY final states has been performed so far.
Turning to associated Higgs production with sparticles, the most important process was expected to be pp → t 1t1 + Higgs which could benefit from the possibly large Higgs-stop coupling [199][200][201][202]. The large value of M S and hence the lightest stop mass from the current constraint makes this process unlikely. Another possibility would be associated production with staus where the phase space could be more favourable but the rates are in general much smaller.
The only channel which could lead to a detectable signal with the data collected so far would be Higgs particles from decays of charginos and neutralinos. In particular the decays χ 0 2 → χ 0 1 h, with χ 0 2 directly produced in association with χ ± 1 in the process pp → χ 0 2 χ ± 1 leading a lepton, a Higgs (decaying either into bb or into multi-leptons via h → Z Z * , W W * ) and missing energy [203][204][205][206][207][208]. Fig. 27 The interpretations of the CMS results from the combination of all lepton and E/ T searches with ≈20 fb −1 data collected at √ s = 8 TeV. The expected and observed limits on the pp → χ 0 2 χ ± 1 cross section times the χ 0 2 χ ± 1 → W hχ 0 1 χ 0 1 branching fraction (with the green band is for experimental uncertainties) is compared to the theoretical prediction [209] The CMS collaboration has reported the results for searches of leptons and missing energy with a luminosity of ≈20 fb −1 data collected at √ s = 8 TeV [209]. They set a limit on the cross section times branching ratio for the possible SUSY process pp → χ 0 2 χ ± 1 with χ 0 2 → χ 0 1 h and χ ± 1 → W χ 0 1 . As can be observed from Fig. 27, where the cross section times branching ratio is displayed as a function of the masses m χ ± 1 = m χ 0 2 (with the assumption that the LSP neutralino is very light, m χ 0 1 = 1 GeV), the data show no excess over the SM backgrounds.

What next?
The last 2 years were extremely rich and exciting for particle physics. With the historical discovery at the LHC of a Higgs boson by the ATLAS and CMS collaboration crowned by a Nobel price this fall, and the first probe of its basic properties, they witnessed a giant step in the unravelling of the mechanism that breaks the electroweak symmetry and generates the fundamental particle masses. They promoted the SM as the appropriate theory, up to at least the Fermi energy scale, to describe three of Nature's interactions, the electromagnetic, weak and strong forces, However, it is clear that these 2 years have also led to some frustration, as no signal of physics beyond the SM has emerged from the LHC data. The hope of observing some signs of the new physics models that were put forward to address the hierarchy problem, which is deeply rooted in the Higgs mechanisms, with supersymmetric theories being the most attractive ones, did not materialise.
The discovery of the Higgs boson and the non-observation of new particles has nevertheless far reaching consequences for supersymmetric theories and, in particular, for their simplest low-energy formulation, the MSSM. The mass of approximately 125 GeV of the observed Higgs boson implies that the scale of SUSY-breaking is rather high, at least O(TeV). This is backed up by the limits on the masses of strongly interacting SUSY particles set by the ATLAS and CMS searches, which in most cases exceed the TeV range [86,87]. This implies that if SUSY is indeed behind the stabilisation of the Higgs mass against very high scales that enter via quantum corrections, it is either fine-tuned at the permille level at least or its low-energy manifestation is more complicated than expected.
The production and decay rates of the observed Higgs particles, as well as its spin and parity quantum numbers, as measured by the ATLAS and CMS collaborations with the ≈ 25 fb −1 data collected at √ s = 7 and 8 TeV, indicate that its couplings to fermions and gauge bosons are approximately SM-like. In the context of the MSSM, this implies that we seem to be in the decoupling regime and this 125 GeV particle can be only identified with the lightest h boson, while the other H/A/H ± states must be heavier than approximately the Fermi scale. This last feature is also backed up by the constraints from direct searches of these heavier Higgs states at the LHC.
This drives up to the question that is now very often asked in particle physics (and elsewhere): what to do next? The answer is, for me, obvious: we are only in the beginning of a new era. 6 Indeed, it was expected since a long time that the probing of the EWSB mechanism will be at least a two chapters story. The first one is the search and the observation of a Higgs-like particle that will confirm the scenario of the SM and most of its extensions, that is, a spontaneous symmetry breaking by a scalar field that develops a non-zero vacuum expectation value. This long chapter has just been closed by the ATLAS and CMS collaborations with the spectacular observation of a Higgs boson. This observation opens a second and equally important chapter: the precise determination of the Higgs profile and the unravelling of the EWSB mechanism itself.
A more accurate measurement of the Higgs couplings to fermions and gauge bosons will be mandatory to establish the exact nature of the mechanism and, eventually, to pin down effects of new physics if additional ingredients beyond those of the SM are involved. This is particularly true in weakly interacting theories such as SUSY in which the quantum effects are expected to be small. These measurements could be performed at the upgraded LHC with an energy close to 6 One can rightfully use here the words of Winston Churchill in November 1942 after the battle of El Alamein (which in Arabic literally means "the two flags" but could also mean "the two worlds" or even "the two scientists"!): "Now, this is not the end; it is not even the beginning to the end; but it is, perhaps, the end of the beginning". √ s = 14 TeV, in particular if a very high luminosity, a few ab −1 , is achieved [210][211][212].
At this upgrade, besides improving the measurements performed so far, rare but important channels such as associated Higgs production with top quarks, pp → tth, and Higgs decays into μ + μ − and Z γ states could be probed. Above all, a determination of the self-Higgs coupling could be made by searching for double Higgs production e.g. in the gluon fusion channel gg → hh [213][214][215]; this would be a first step towards the reconstruction of the scalar potential that is responsible of EWSB. A proton collider with an energy √ s = 30 to 100 TeV could do a similar job [212]. In a less near future, a high-energy lepton collider, which is nowadays discussed in various options (ILC, TLEP, CLIC, μ collider) would lead to a more accurate probing of the Higgs properties [216][217][218][219][220][221][222][223][224], promoting the scalar sector of the theory to the high-precision level of the gauge and fermionic sectors achieved by LEP and SLC [30].
Besides the high precision study of the already observed Higgs, one should also continue to search for the heavy states that are predicted by SUSY, not only the superparticles but also the heavier Higgs bosons. The energy upgrade to ≈14 TeV (and eventually beyond) and the planed order of magnitude (or more) increase in luminosity will allow to probe much higher mass scales than presently.
In conclusion, it is not yet time to give up on supersymmetry and on new physics in general but, rather, to work harder to be fully prepared for the more precise and larger data that will be delivered by the upgraded LHC. It will be soon enough to "philosophise" in 2 years from now, when the physics landscape will become clearer.