Complementarity of direct and indirect searches in the pMSSM

We explore the pMSSM parameter space in view of the constraints from SUSY and monojet searches at the LHC, from Higgs data and flavour physics observables, as well as from dark matter searches. We show that whilst the simplest SUSY scenarios are already ruled out, there are still many possibilities left over in the pMSSM. We discuss the complementarity between different searches and consistency checks which are essential in probing the pMSSM and will be even more important in the near future with the next round of data becoming available.


Introduction
Data from the first run of the LHC and from dark matter (DM) direct detection experiments provide us already with important constraints on the supersymmetric parameter space. We consider the phenomenological MSSM (pMSSM), which is the most general MSSM set-up with R-parity and CP conservation respecting the minimal flavour violation (MFV) paradigm. With its 19 independent parameters, the pMSSM offers sufficient freedom to the masses and couplings to explore the supersymmetric parameter space and the implications of the recent particle and astroparticle physics data in a largely unbiased way. The 19 parameters are scanned over in the ranges given in Table 1 following the methodology described in Refs. [1,2].
In the following, for the relic density, we apply the upper bound of the cold dark matter density measured by Planck [15], allowing for the possibility that dark matter is composed of multiple components or that the early Universe properties differ from that of the standard cosmological model [16,17]. The constraints from flavour physics and low energy data are also applied as described in [18].

SUSY and monojet searches at the LHC
The direct searches for MET signatures with jets and leptons provide important constraints on the strongly interacting SUSY particle masses. However, these searches lose sensitivity in the regions where the mass splittings are small. Such cases often lead to soft jets which are invisible in the detectors. Monojet searches on the other hand are efficient in such situations and can be used to improve the sensitivity [19]. In Fig. 1, the effect of the SUSY (solid line) and monojet (dotted line) search results on the squark and gluino masses is presented. It is remarkable that while the limits on the squark and gluino masses are above 1.5-2 TeV in the constrained or simplified MSSM scenarios, in a more general scenario such as the pMSSM with no ad hoc universality assumption for the SUSY particle masses, this is not any more the case. The figure reveals that many model points with masses below 1 TeV are still allowed and that 30% of the points have still squark masses below 800 GeV and gluino masses below 1.5 TeV.
To investigate further the complementarity between the SUSY and monojet searches we demonstrate in Fig. 2 the fraction of excluded points, for the lightest squark/gluino mass vs. neutralino mass, and for the mass splitting between the lightest squark/gluino and neutralino vs. the neutralino mass. The plots show that squark and gluino masses below 400 GeV are now strongly disfavoured regardless of the lightest neutralino mass, and that the monojet searches are very sensitive in the region of small mass splittings and can improve the reach by more than 100 GeV.

Higgs searches
The discovery of a Higgs boson of ∼126 GeV at the LHC has strong implications for supersymmetry. In particular, it is generally assumed that the discovered Higgs boson corresponds to the lightest MSSM CP-even state. The mass of the lightest CP-even Higgs boson is related in the MSSM to the stop sector parameters, and in particular to the geometric average of the stop masses M S and to the stop mixing parameter X t [20]. In Fig. 3, we show that for M S < 1 TeV, a Higgs mass of 126 GeV can only be reached for M S ∼ √ 6 |X t |, corresponding to the maximal mixing scenario, while other mixings can be allowed for larger M S [21].
In Fig. 4 we present the distribution of the pMSSM points compatible with the h boson mass and the observed yields, in the (Mt 1 , X t ) plane [22,23]. We notice again that small values of mixing parameter |X t | are  clearly disfavoured, and that stop masses as low as 400 GeV are still compatible with the data.
The determination of the Higgs branching fractions to two photons, two W or Z bosons, and to two fermions has strong implications, in particular when complemented by the heavier Higgs state searches [24]. In Fig. 5, we show the constraints on the (M A , tan β) plane from the determination of the light Higgs branching ratios and from searches for heavier states. As can be seen, the combination of the H/A → τ + τ − channel and the mass and measured branching ratios for the lightest h boson exclude the region with M A < 320 GeV for all values of tan β. The ZZ channel, and to a lesser extent the WW channel, close the low M A corner from tan β 2 up to the ττ limit for M A < 230 GeV.
The constraints derived by the study of the Higgs sector are becoming an essential part of the probe of the SUSY parameter space at the LHC.

Flavour physics
In addition to the direct searches, indirect information from the flavour physics observables are of great interest in probing the supersymmetric parameter space. One of the most constraining observables is the rare leptonic decay of B s → µ + µ − , as it can receive supersymmetric contributions scaling with tan 6 β/M 4 A , increasing the branching ratio by orders of magnitude. This branching ratio has been measured by LHCb and CMS [25] with the central value well compatible with the Standard Model prediction, strongly limiting large deviations due to SUSY contributions [26]. The impact of the present and future determinations of BR(B s → µ + µ − ) on the parameters most sensitive to its rate, (M A , tan β) and (M A , Mt 1 ), is shown in Fig. 6, where we give all the valid pMSSM points from our scan (black points), those with 123 < M h < 129 GeV (grey points) and, highlighted in green, those in agreement with the present BR(B s → µ + µ − ) range and the estimated ultimate constraint at 95% C.L. [26]. The constraints from BR(B s → µ + µ − ) affect the same pMSSM region, at large values of tan β and small values of M A , also probed by the H/A → τ + τ − direct Higgs searches. As demonstrated in the previous section, the search for the H/A → τ + τ − decay has already excluded a significant portion of the parameter space where large effects on BR(B s → µ + µ − ) are expected. We also note that the stop sector is further constrained by direct searches in b-jets + MET channels, which disfavour small values of Mt 1 . The figure shows that it is difficult for M A and Mt 1 to be simultaneously light.

Dark matter direct detection
Indirect constraints can also be obtained from dark matter direct detection experiments. They rely on the measurement of the recoil energy of a nucleus with which a dark matter particle may interact. The interpretation of dark matter direct search results in terms of scattering cross sections of dark matter with matter is however subject to strong assumptions on the local density and velocity of dark matter. Other dark matter observables can be considered additionally. The dark matter relic density is sensitive to the annihilation cross section of neutralinos, but also of other SUSY particles which are close in mass to the neutralino. Other constraints can come from indirect detection, which relies on the detection of SM particles produced through the annihilation of dark matter particles. The results suffer strongly from astrophysical assumptions and therefore we do not consider them here.
A well motivated candidate in the MSSM for dark matter is the lightest neutralino. Its scattering with matter is expected to be mediated mainly by a Z boson or a neutral Higgs boson. For this reason, direct detection also probes the (M A , tan β) plane which is also constrained by heavy Higgs searches and flavour physics observables as discussed in the previous sections. This redundancy obtained through searches in multiple channels sensitive in the same regions is essential to fully probe the (M A , tan β) parameter space as it is possible in specific cases to furtively evade one of the constraints, and also for the purpose of consistency checks. The LUX collaboration presented recently stringent limits on the scattering cross sections of dark matter particles with matter [28]. On the other hand, several direct detection experiments had reported excess of events in the light dark matter regions of a few GeV up to tenths of GeV [29][30][31][32]. While it is possible to explain such excesses in the pMSSM [18,33], the LUX limits seem to exclude the regions reported by those experiments.
In Figs. 7 and 8, we present the constraints from LUX on M A and tan β. As can be seen in the figures, important constraints are obtained by dark matter direct detection for large values of tan β. We notice also that small M A < 200 GeV are excluded by these searches independently of any other constraints, and that 90% of the points with M A 350 GeV are ruled out.
In Fig. 9, the spin-independent scattering cross section of a neutralino with matter is shown as a function of the dark matter mass, together with current and prospective limits of direct detection experiments [34]. The blue region represents the spread of our pMSSM points  : Dark matter direct detection spin-independent scattering cross section of a neutralino with a nucleon as a function of the dark matter mass (from [34]). The dotted red line corresponds to the current LUX limits. The solid black line is the projected limit of the future LZ experiment. The black dashed line corresponds to the neutrino background limit. The blue zone shows the spread of our pMSSM points in this plane.
in this plane. The LUX limit excludes about 30% of the model points. This fraction will become larger with the future experiment results becoming available, however, very small scattering cross sections can be obtained in the pMSSM so that even by reaching the neutrino background limit still a non-negligible fraction of the points would survive that could in a large extent be probed at the LHC high luminosity run or at future colliders.

Conclusion
In view of the negative results of the searches in channels with MET conducted by the ATLAS and CMS experiments at the end of the 7 and 8 TeV runs, it is clear that the simplest SUSY benchmarks are now ruled out and it is of great importance to reinterpret the results in more general scenarios with minimal theoretical assumptions. By studying the 19-parameter phenomenological MSSM, we have demonstrated that a wide parameter space of solutions with SUSY masses below 1 TeV, compatible with all the current data still exists. Even at the end of the next LHC run, many of these solutions will not be tested. In such a situation, the only way to point to a specific scenario, or falsify the MSSM, would be to take advantage of interplay between different sectors. Regarding the LHC direct searches, in addition to direct SUSY search results, information from the Higgs sector namely from the mass and coupling measurements as well as information from other-than-SUSY searches such as monoparticle searches can be used to tightly constrain the parameter space. Moreover, indirect search results in particular from the flavour physics sector and rare decays, and also from the dark matter searches will allow us to deeply probe the MSSM parameter space. Such complementarities were discussed in this report. In particular, we showed that the monojet searches are sensitive to the regions with small mass splittings where the SUSY searches are weaker. Also the Higgs searches allow us to explore efficiently the (M A , tan β) plane which is also probed by both data from B s → µ + µ − and from dark matter direct detection. These complementarities will be further exploited with the results of the next LHC runs, future dark matter detection experiments as well as the future Super-B factory.