Interplay of LHC and dark matter searches in the MSSM 1

PARTICLE PHYSICS RESEARCH ARTICLE Interplay of LHC and dark matter searches in the MSSM 1 Alexandre Arbey 1,2,3 *, and Farvah Mahmoudi 4 Received: 10 April 2015 Accepted: 16 June 2015 Published: 24 July 2015 *Corresponding author: Alexandre Arbey, Université de Lyon, Université Lyon 1, F-69622 Villeurbanne Cedex, France; Centre de Recherche Astrophysique de Lyon, CNRS, UMR 5574, F-69561 Saint-Genis Laval Cedex, France; Ecole Normale Supérieure de Lyon, Lyon Cedex 07, France Email: alexandre.arbey@ens-lyon.fr Reviewing editor: Anna Lipniacka, University of Bergen, Norway Additional information is available at the end of the article Abstract: The Minimal Supersymmetric extension of the Standard Model (MSSM) provides suitable candidates for cold dark matter. We discuss here the constraints from dark matter direct detection and cosmological dark matter density, as well as LHC data from Higgs, SUSY and monojet searches, and flavour physics data, in the context of the phenomenological MSSM (pmssm) with neutralino dark matter. We show that the complementarity of the different sectors is essential to probe the pmssm parameter space. Subjects: Cosmology; Particle & High Energy Physics; Theoretical Physics Keywords: dark matter; LHC; supersymmetry; Higgs 1. Introduction The Minimal Supersymmetric extension of the Standard Model (MSSM) is one of the most studied new physics models. When R-parity is conserved, supersymmetric particles can only be produced in pairs off Standard Model (SM) particles, and the lightest supersymmetric particle (LSP) is stable. If the LSP is neutral and weakly interacting, it can constitute a suitable dark matter candidate. This is the case of the lightest neutralino, which interacts weakly enough to make an ideal dark matter candidate. Studies in the MSSM however suffer from the large number of parameters more than a hundred which makes a systematic analysis impossible. On the other hand, constrained MSSM scenarios, such as the constrained MSSM, have only a handful of parameters, but cover only a small region of the MSSM parameter space. The phenomenological MSSM (pmssm) (Djouadi, XXXX) offers a convenient compromise: it is the most general version of the MSSM with CP and R-parity conservations respecting minimal flavour violation, which has 19 independent parameters. In this work, we will concentrate on the pmssm with neutralino dark matter. We study its 19 parameter space by performing flat scans on all the parameters, as described in Arbey, Battaglia, and Mahmoudi (2012a, 2012b), for masses up to 3.5 TeV. Other analyses in the pmssm can be found for example in AbdusSalam, Allanach, Quevedo, ABOUT THE AUTHORS Alexandre Arbey and Farvah Mahmoudi are associate professors at University of Lyon and members of the CERN Theory Division. A Arbey is a cosmologist and astroparticle physicist, working on the Dark Matter problem, as well as dark energy and Big-Bang nucleosynthesis. F Mahmoudi is a particle physicist, specialist of Higgs, Flavour and Dark Matter physics. Their work consists in combining the constraints from all their sectors of expertise to constrain New Physics models. The project reported in this effort is performed in the context of Supersymmetry, but it can extend to other categories of particle physics scenarios beyond the Standard Model. PUBLIC INTEREST STATEMENT Cosmological and astrophysical observations reveal that dark matter constitutes the major part of matter in the Universe. Yet, the current Standard Model of particle physics does not provide any particle candidate to explain its nature. The supersymmetric extension of the Standard Model suggests that dark matter could be made of neutralinos, which is a neutral and weakly interacting particle not been discovered yet. In this paper, the authors investigate the current constraints on supersymmetric dark matter, by considering simultaneously the results from cosmological observations, dark matter particle detection experiments, as well as LHC searches for new particles. 2015 The Author(s). This open access article is distributed under a Creative Commons Attribution (CC-BY) 4.0 license. Page 1 of 10

Feroz and Hobson (2010), Conley, Gainer, Hewett, Le, and Rizzo (XXXX), Sekmen et al. (2012), Cahill- Rowley, Hewett, Hoeche, Ismail, and Rizzo (2012), CMS Collaboration (XXXXa, XXXXb), Kowalska, Roszkowski, Sessolo, and Trojanowski (2014), Cahill-Rowley et al. (2015). In the following, we first consider LHC constraints from supersymmetry (SUSY) and monojet searches. In Section 3, we study the consequences of the Higgs boson discovery. In Section 4, we account for constraints from flavour physics. The results of dark matter searches are discussed in Section 5. Finally, we consider all the constraints simultaneously to assess their complementarity, before concluding in Section 6. 2. SUSY and monojet searches at the LHC We first consider the constraints from SUSY and monojet searches at the LHC in the context of the pmssm. Since supersymmetric particles are pair produced and the LSP is stable, supersymmetric search channels at the LHC feature elusive neutralinos in the final state. Hence, the experimental signatures consist in final states with missing transverse energy (MET). In this analysis, we consider the jets+met channels, relevant for the search for first and second generation squarks and gluinos (ATLAS Collaboration, XXXXa), as well as b-tagged jets+met channels for sbottom and stop searches (ATLAS Collaboration, XXXXb), and chargino/neutralino searches in two leptons+met (ATLAS Collaboration, XXXXc) and three leptons+met (ATLAS Collaboration, XXXXd) channels. We generate the events with PYTHIA 8.1 (Sjostrand, Mrenna, & Skands, 2008), with the CTEQ6L1 parton distribution functions (Pumplin et al., 2002). Detector simulation is performed with DELPHES 3.0 (Ovyn, Rouby, & Lemaitre, XXXX). Contrary to the previous analyses in the pmssm, here, we systematically use the monojet search results, which are particularly important in the case of compressed spectra. To study the consequences of the monojet searches on the pmssm parameter space, we have designed a specific method to reinterpret the monojet search data in the context of the MSSM, as described in Arbey, Battaglia and Mahmoudi (2014). Monojet channels consist of one hard jet+met in the final state. The results are usually presented in effective or simplified models (ATLAS Collaboration, XXXXe; CMS Collaboration, XXXXc), in which the final state is interpreted as generated by two dark matter particles and one jet, as shown in the left panel of Figure 1. However, within a specific model like the MSSM, there are other means to generate a monojet final state at the LHC. In particular, since the LHC is a hadronic machine, the production of squarks and gluinos which have larger cross sections is more likely. Moreover, in compressed MSSM scenarios, the squarks, gluinos and neutralinos can have similar masses, so that the squark or gluino would mostly decay to soft jet(s) + one neutralino. Due to the cuts applied in the monojet searches, the soft jet(s) will not be observed, and the final state with two squarks/gluinos + one hard jet will be the same as the typical monojet signature. Such a topology is presented in the right panel of Figure 1. For this reason, to assess the implications of the monojet searches in the MSSM, it is important to consider the processes pp (N)LSP + (N)LSP + j, where j is a hard jet and (N)LSP stands for the (next-to-)lsps. We compute the matrix elements corresponding to the full diagrams with MadGraph 5 (Alwall et al., 2014) and generate events using the CTEQ6L1 parton distribution functions. The parton showering is performed with PYTHIA 6 (Sjostrand, Mrenna, & Skands, 2006), and detector simulation with DELPHES 3.0. We consider the ATLAS and CMS monojet results (ATLAS Collaboration, XXXXe; CMS Collaboration, XXXXc) to obtain the exclusion of the pmssm points by the monojet searches. Figure 1. Schematic representation of monojet channels in the effective or simplified models (left panel) and additional diagrams in the MSSM with neutralino dark matter (right panel). Page 2 of 10

In Figure 2, we show the difference of exclusion potential between the case where only neutralino LSPs are considered in the monojet final states (left panel), and the case with in addition squarks and gluinos decaying to soft jets and neutralinos (right panel), for the LHC 14 TeV run with 300 fb 1. As can be seen from the figure, considering the squark and gluino final states would lead to a large increase in the monojet production cross section, enhancing substantially the monojet sensitivity. Figure 2. Monojet production cross section as a function of the neutralino mass, in the case where only neutralino LSPs are considered in the monojet final states (left panel), and in the case with squarks and gluinos decaying to soft jets and neutralinos (right panel), for the LHC 14 TeV run with 300 fb 1. Notes: The blue points satisfy the monojet constraints, while the red points are excluded at 95% C.L. The current exclusion limits by the SUSY and monojet searches are presented in Figure 3 as functions of the lightest squark and gluino masses. As expected the direct SUSY searches strongly probe the small squark and gluino masses, and the monojet searches nicely complement them up to an Figure 3. Fraction of excluded points by the SUSY searches (solid lines) in addition to the monojet searches (dotted lines), as a function of the lightest squark mass (left panel) and of the gluino mass (right panel). Figure 4. Fraction of excluded points in the lightest squark/ gluino vs. neutralino mass plane (left panel) and the lightest squark/gluino mass splitting with the neutralino vs. neutralino mass plane (right panel). Notes: The grey zones correspond to regions where no valid pmssm point is found. The grey line delimits the region where 99% of the points are excluded by the SUSY searches only, and the yellow line the region where the points are excluded by the combination of the SUSY and monojet searches (Arbey et al., 2014). Page 3 of 10

exclusion of close to 100% of the points for masses below 400 GeV. In Figure 4, the excluded points are presented in two-dimensional plots, in the lightest squark/gluino vs. neutralino mass plane and in the lightest squark/gluino mass splitting with the neutralino vs. neutralino mass plane. We see that the SUSY searches can probe squark/gluino masses up to 1.5 TeV when the neutralino is light, whereas the monojet searches improve the constraints on the neutralino mass by 100 GeV for small mass splittings. This shows the important complementarity between the SUSY and monojet searches, in particular when investigating compressed scenarios with small mass splittings. 3. Higgs searches The discovery of a Higgs boson at 125 GeV (Aad et al., 2012; Chatrchyan et al., 2012) has strong implications for the MSSM, as it can be considered to be the lightest MSSM Higgs boson h. In the MSSM, the lightest Higgs mass and couplings are mainly determined by the parameters M A (CP-odd Higgs mass), M S (geometric average of the stop masses), X t = A t μ cot β (stop mixing parameter) and tan β. To obtain a Higgs mass of 125 GeV, it is necessary to have a large M S or a large X t (Arbey, Battaglia, Djouadi, Mahmoudi, & Quevillon, 2012). In our scans, we impose a Higgs mass of 125 ± 3 GeV to account for the experimental and theoretical uncertainties, and consider the measured signal strengths in the γγ, ZZ, W + W, b b and τ + τ channels, and impose the limits from heavier Higgs H A τ + τ searches, as described in Arbey, Battaglia, Djouadi, and Mahmoudi (2012a,2012b). The results are presented in Figure 5 in the (M A, tan β) and (M t 1, X t ) parameter planes. It is clear that large values of M A ( 350 GeV) and X t are favoured by the current measurements. The searches for heavier Higgs states can impose further constraints on the MSSM parameter space. In Figure 6, the constraints obtained from the searches at 8 TeV for H τ + τ, H ZZ and H b b are shown in the left panel. H τ + τ gives the strongest constraints, particularly important at large tan β. The H ZZ channel on the other hand provides constraints in the small tan β region, which are complementary to the constraints from the light Higgs signal strength measurements. In the right panel, extrapolations for the 14 TeV run are provided for the same channels, as well as for H t t. Figure 5. pmssm points in agreement with the Higgs measurements projected on the (M A, tan β) (left panel) and (M t 1, X t ) (right panel) parameter planes. Notes: The black points are in agreement with the Higgs mass and H A τ + τ constraints. The dark and light green points are compatible at 68% and 95% C.L. with the additional h signal strength measurements, respectively (Arbey et al., 2012b). Figure 6. Fraction of excluded points (blue scale) by the heavy Higgs searches for the 8 TeV run (left panel) and extrapolations for the 14 TeV run (right panel). Notes: The lines delimit the excluded regions from different heavy Higgs search channels (Arbey, Battaglia, & Mahmoudi, 2013). Page 4 of 10

Figure 7. Constraints from BR(B s μ + μ ) in the (M A, tan β) (left panel) and (M A, M t 1 ) (right panel) parameter planes. Notes: The black points correspond to the valid pmssm points, the grey points are also compatible with the Higgs mass constraints. The dark green points are in agreement at 95% C.L. with the current BR(B s μ + μ ) constraints, and the light green points show the projection for the ultimate LHCb constraints from BR(B s μ + μ ). For comparison, the solid line in the left panel shows exclusion limits by Higgs searches for H A τ + τ (Arbey, Battaglia, Mahmoudi, & Martinez Santos, 2013). As can be seen, H τ + τ searches will probe a much larger parameter space, while H t t will provide complementary information at small tan β. 4. Flavour physics constraints B physics observables provide also additional constraints on the pmssm parameter region. In particular, the branching ratio of B s μ + μ is of interest since it can receive large contributions from scalar and pseudoscalar operators in the MSSM. Using SuperIso (Mahmoudi, 2008, 2009), we have computed and studied the constraints from this observable on the pmssm. In Figure 7, we present the constraints from BR(B s μ + μ ) in the (M A, tan β) and (M A, M t 1 ) parameters planes. As can be seen, the branching fraction of B s μ + μ is particularly constraining at large tan β and disfavours small M A and M t 1. In this respect, it is complementary to the Higgs searches, and in particular the H A τ + τ channel is sensitive to the same parts of the parameter space as BR(B s μ + μ ). Other important constraints from flavour physics are provided by BR(B X s γ), BR(B X s ll), BR(B τν) and BR(B K ll), which can also set limits on the pmssm parameters (Mahmoudi, Neshatpour, & Virto, 2014). 5. Dark matter searches Dark matter can be considered as currently the only (indirect) evidence for new physics. Cold dark matter is likely to be made of weakly interacting massive particles (WIMPs), so elusive that none has been detected yet. Cold dark matter constitutes about 23% of the total energy density of the Universe. In a scenario like the pmssm with neutralino dark matter, the density of dark matter particles, the so-called relic density, can be computed. This relic density can then be compared to the observed cosmological density of dark matter to set constraints on the MSSM parameters, namely: 0.076 < Ω χ h 2 < 0.163. (1) However, the calculation of the relic density is based on assumptions on the cosmological past of the Universe, in particular before Big-Bang nucleosynthesis where the Universe cannot be observed. Small deviations from these assumptions can lead to large differences (in general increases) in the relic density (Arbey & Mahmoudi, 2008, 2010a; Gelmini, Gondolo, Soldatenko, & Yaguna, 2006; Kamionkowski & Turner, 1990). In addition, dark matter may be made of several components, and the neutralino can therefore constitute only part of it. For these reasons, we prefer to apply a loose relic density constraint: 10 4 < Ω χ h 2 < 0.163. (2) We compute the relic density with SuperIso Relic (Arbey & Mahmoudi, 2010b). In Figure 8, we show the relic density as a function of the neutralino mass, for the different types of neutralinos. While a bino leads in general to a too large relic density because of its suppressed couplings to the Higgs and Z bosons, a higgsino of about 1.5 TeV and a wino of about 2.7 TeV can lead naturally to a relic density compatible with the measured dark matter density. For binos, it is mandatory to have either resonant Page 5 of 10

Figure 8. Relic density as a function of the lightest neutralino mass. Notes: The colours of the points correspond to the neutralino types. The black line gives the central value of the observed dark matter density. annihilation channels or coannihilation channels to get a small enough relic density. Therefore, the relic density constraints have strong consequences on the selection of the neutralino natures. Other important constraints come from the dark matter direct detection experiments such as XENON 100 (Aprile et al., 2012) and LUX (Akerib et al., 2014). These experiments are sensitive to the scattering of WIMPs with nuclei, and can set limits on the WIMP mass as well as their spin-independent scattering cross sections with protons in absence of detection 2. We compute the neutralino scattering cross section with micromegas (Belanger, Boudjema, Pukhov, & Semenov, 2009). In Figure 9, we present the neutralino spin-independent scattering cross section with protons as a function of the neutralino mass, for the different neutralino types. The higgsinos have larger scattering cross sections, and are therefore currently more probed by the XENON and LUX experiments. Winos have smaller cross sections, and will mostly be probed by the LZ experiment. Finally, binos interact only very weakly and are more likely to escape detection by the direct detection experiments. In Figure 10, the exclusion by the LUX experiment is shown in the (M A, tan β) parameter plane. It is clear that dark matter direct detection restricts mostly large tan β and small M A values. One can notice that direct detection is complementary to the heavy Higgs and flavour physics searches since the same part of the parameter space is probed. Since in addition to M A and tan β other MSSM Figure 9. Spin-independent scattering cross section of neutralinos with protons as a function of the lightest neutralino mass. The colour of the points denotes the dominant nature of the lightest neutralino. Notes: The solid lines show the latest XENON and LUX limits, the dashed line represents the prospective LZ limit and the dotted line represents the limit below which neutrinos will act as a large background for the dark matter searches. Page 6 of 10

Figure 10. Constraints from the LUX experiment in the (M A, tan β) parameter plane. Notes: The black points are excluded by the LUX limits, while the green points are allowed. parameters are also relevant for these searches, such redundancies are important to probe deeply the pmssm parameter space. 6. Summary and conclusions Direct SUSY and monojet searches, Higgs and flavour physics data, as well as dark matter constraints provide complementary information on the pmssm parameters. To highlight the interplay between Figure 11. Fraction of pmssm points excluded by the SUSY direct searches (left panel), in addition to the monojet searches (right panel), in the neutralino proton scattering cross section vs. neutralino mass. Notes: The lines show the XENON-100 and LUX dark matter direct detection limits (Arbey et al., 2014). Figure 12. Fraction of pmssm points excluded by the SUSY direct searches, monojet searches and LUX data as a function of neutralino mass (left panel) and mass splitting between the lightest squark or gluino and the lightest neutralino (right panel). Page 7 of 10

the different observables, in this section we pre-impose at the same time the constraints from Higgs physics, flavour physics and relic density. In Figure 11, we show the fraction of excluded pmssm points in the neutralino scattering cross section with matter vs. neutralino mass parameter plane, considering simultaneously constraints from SUSY direct searches and monojet searches. The limits from LUX are also displayed in the same plane. We first see that the LHC searches are complementary to dark matter direct detection, and probe different parts of the parameter space. In particular, LHC searches are sensitive to the strong sector of the MSSM, while direct dark matter detection is more sensitive to the neutralino sector. We also see that the LHC probes regions well below the LUX limits, and that monojet searches nicely increase the LHC discovery potential in the regions where the neutralino scattering cross section with matter is very small. In Figure 12, the fraction of excluded points is shown as a function of the neutralino mass, and of the mass splitting between the lightest squark or gluino and the lightest neutralino. It shows that the SUSY and monojet searches probe more strongly scenarios with light neutralinos, excluding up to 65% of the points for light neutralinos, and that the dark matter direct detection limits can increase this exclusion up to 75%, presenting a nice complementarity between the two sectors. Additionally, the LHC searches are very powerful for mass splittings up to 700 GeV where more than 50% of the points are excluded, while direct detection is less dependent on the mass splittings. To conclude, we have demonstrated the synergy between the different observables of separate sectors, namely the LHC, flavour and dark matter sectors. During the next LHC run, this interplay will become even more important, either in the case no new phenomenon is found to ensure that no region of the parameter space remains uncovered, or in the case new phenomena or particles are discovered to allow for a determination of the underlying scenario and parameters. Acknowledgements This work was done in collaboration with Marco Battaglia to whom the authors are grateful. AA thanks the organizers of the From Higgs to Dark Matter 2014 Geilo Symposium for their invitation and support, and acknowledges partial support from the PEPS Physique Théorique et ses Interfaces of the CNRS. Funding The authors received no direct funding for this research. Author details Alexandre Arbey 1,2,3 E-mail: alexandre.arbey@ens-lyon.fr Farvah Mahmoudi 4 E-mail: mahmoudi@in2p3.fr 1 Université de Lyon, Université Lyon 1, F-69622Villeurbanne Cedex, France. 2 Centre de Recherche Astrophysique de Lyon, CNRS, UMR 5574, F-69561 Saint-Genis Laval Cedex, France. 3 Ecole Normale Supérieure de Lyon, Lyon Cedex 07, France. 4 CERN Theory Division, Physics Department, CH-1211 Geneva 23, Switzerland. Citation information Cite this article as: Interplay of LHC and dark matter searches in the MSSM, Alexandre Arbey, & Farvah Mahmoudi, Cogent Physics (2015), 2: 1064575. Notes 1. This article is associated with the From Higgs to Dark Matter 2014 meeting in Geilo, Norway. 2. It should be noted that several independent dark matter experiments had claimed for possible detection of WIMPs of the order of 10 GeV (Bernabei et al., 2010; Aalseth et al., 2011, Angloher et al., 2012, Ahmed et al., 2011. While the pmssm can provide scenarios compatible with such data (Arbey, Battaglia, & Mahmoudi, 2012, 2013), they are actually disfavoured by the latest LUX results. References Aad, G., et al. (ATLAS Collaboration). (2012). Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC. Physics Letters B, 716, 1 (arxiv:1207.7214 [hep-ex]). Aalseth, C. E., Barbeau, P. S., Colaresi, J., Collar, J. I., Diaz Leon, J., Fast, J. E.,... Yocum, K. M. (2011). Search for an annual modulation in a P-type point contact germanium Dark Matter detector. Physical Review Letters, 107, 141301 (arxiv:1106.0650 [astro-ph.co]). AbdusSalam, S. S., Allanach, B. C., Quevedo, F., Feroz, F., & Hobson, M. (2010). Fitting the phenomenological MSSM. Physical Review D, 81, 095012 (arxiv:0904.2548 [hep-ph]). Ahmed, Z., et al. (CDMS-II Collaboration). (2011). Results from a low-energy analysis of the CDMS II germanium data. Physical Review Letters, 106, 131302 (arxiv:1011.2482 [astro-ph.co]). Akerib, D. S., et al. (LUX Collaboration). (2014). First results from the LUX dark matter experiment at the Sanford Underground Research Facility. Physical Review Letters, 112, 091303 (arxiv:1310.8214 [astro-ph.co]). Alwall, J., Frederix, R., Frixione, S., Hirschi, V., Maltoni, F., Mattelaer, O.,... Zaro, M. (2014). The automated computation of tree-level and next-to-leading order Page 8 of 10

differential cross sections, and their matching to parton shower simulations. Journal of High Energy Physics, 1407, 079 (arxiv:1405.0301 [hep-ph]). Angloher, G., Bauer, M., Bavykina, I., Bento, A., Bucci, C., Ciemniak, C., Deuter, G.,... Zöller, A. (2012). Results from 730 kg days of the CRESST-II Dark Matter search. European Physical Journal C, 72, 1971 (arxiv:1109.0702 [astro-ph.co]). Aprile, E., et al. (XENON100 Collaboration). (2012). Dark Matter results from 225 live days of XENON100 data. Physical Review Letters, 109, 181301 (arxiv:1207.5988 [astro-ph.co]). Arbey, A., Battaglia, M., Djouadi, A., & Mahmoudi, F. (2012). The Higgs sector of the phenomenological MSSM in the light of the Higgs boson discovery. Journal of High Energy Physics, 1209, 107 (arxiv:1207.1348 [hep-ph]). Arbey, A., Battaglia, M., Djouadi, A., Mahmoudi, F., & Quevillon, J. (2012). Implications of a 125 GeV Higgs for supersymmetric models. Physics Letters B, 708, 162 (arxiv:1112.3028 [hep-ph]). Arbey, A., Battaglia, M., Djouadi, A., & Mahmoudi, F. (2013). An update on the constraints on the phenomenological MSSM from the new LHC Higgs results. Physics Letters B, 720, 153 (arxiv:1211.4004 [hep-ph]). Arbey, A., Battaglia, M., Mahmoudi, F., & Martinez Santos, D. (2013). Supersymmetry confronts Bs μ + μ - : Present and future status. Physical Review D, 87(3), 035026 (arxiv:1212.4887 [hep-ph]). Arbey, A., Battaglia, M., & Mahmoudi, F. (2012a). Light neutralino Dark Matter in the pmssm: Implications of LEP, LHC and Dark Matter searches on SUSY particle spectra. European Physical Journal C, 72, 2169 (arxiv:1205.2557 [hep-ph]). Arbey, A., Battaglia, M., & Mahmoudi, F. (2012b). Implications of LHC searches on SUSY particle spectra: The pmssm parameter space with neutralino Dark Matter. European Physical Journal C, 72, 1847 (arxiv:1110.3726 [hep-ph]). Arbey, A., Battaglia, M., & Mahmoudi, F. (2012c). Constraints on the MSSM from the Higgs sector: A pmssm study of Higgs searches, B 0 s μ+ μ - and Dark Matter direct detection. European Physical Journal C, 72, 1906 (arxiv:1112.3032 [hep-ph]). Arbey, A., Battaglia, M., & Mahmoudi, F. (2013a). Supersymmetry with Light Dark Matter confronting the recent CDMS and LHC results. Physical Review D, 88, 095001 (arxiv:1308.2153 [hep-ph]). Arbey, A., Battaglia, M., & Mahmoudi, F. (2013b). Supersymmetric heavy Higgs Bosons at the LHC. Physical Review D, 88(1), 015007 (arxiv:1303.7450 [hep-ph]). Arbey, A., & Mahmoudi, F. (2008). SUSY constraints from relic density: High sensitivity to pre-bbn expansion rate. Physics Letters B, 669, 46 (arxiv:0803.0741 [hep-ph]). Arbey, A., & Mahmoudi, F. (2010a). SUSY constraints, relic density, and very early universe. Journal of High Energy Physics, 1005, 051 (arxiv:0906.0368 [hep-ph]). Arbey, A., & Mahmoudi, F. (2010b). SuperIso relic: A program for calculating relic density and flavor physics observables in supersymmetry. Computer Physics Communications, 181, 1277 (arxiv:0906.0369 [hep-ph]). Arbey, A., Battaglia, M., & Mahmoudi, F. (2014). Combining monojet, supersymmetry, and Dark Matter searches. Physical Review D, 89(7), 077701 (arxiv:1311.7641 [hep-ph]). ATLAS Collaboration. (XXXXa). Search for squarks and gluinos in events with jets and MET. Note ATLAS-CONF-2013-047. ATLAS Collaboration. (XXXXb). Search for direct third generation squark pair production in final states with missing transverse momentum and two b-jets in s = 8 TeV pp collisions with the ATLAS detector. Note ATLAS-CONF-2013-053. ATLAS Collaboration. (XXXXc). Searches for direct slepton-pair and chargino-pair production in final states with two opposite-sign leptons, missing transverse momentum and no jets in 20 fb 1 of pp collisions at s = 8 TeV with the ATLAS detector. Note ATLAS-CONF-2013-049. ATLAS Collaboration. (XXXXd). Search for direct gaugino production in events with three or more leptons and MET. Note ATLAS-CONF-2013-035. ATLAS Collaboration. (XXXXe). Search for pair-produced top squarks decaying into a charm quark and the lightest neutralinos with 20.3 fb -1 of pp collisions at s=8~tev with the ATLAS detector at the LHC. ATLAS-CONF-2013-068. Belanger, G., Boudjema, F., Pukhov, A., & Semenov, A. (2009). Dark matter direct detection rate in a generic model with micromegas 2.2. Computer Physics Communications, 180, 747 (arxiv:0803.2360 [hep-ph]). Bernabei, R., et al. (DAMA and LIBRA Collaborations). (2010). New results from DAMA/LIBRA. European Physical Journal C, 67, 39 (arxiv:1002.1028 [astro-ph.ga]). Cahill-Rowley, M., Cotta, R., Drlica-Wagner, A., Funk, S., Hewett, J., Ismail, A.,... Wood, M. (2015). Complementarity of Dark Matter searches in the phenomenological MSSM. Physical Review D, 91(5), 055011 (arxiv:1405.6716 [hep-ph]). Cahill-Rowley, M. W., Hewett, J. L., Hoeche, S., Ismail, A., & Rizzo, T. G. (2012). The new look pmssm with neutralino and gravitino LSPs. European Physical Journal C, 72, 2156 (arxiv:1206.4321 [hep-ph]). Chatrchyan, S., et al. (CMS Collaboration). (2012). Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC. Physics Letters B, 716, 30 (arxiv:1207.7235 [hep-ex]). CMS Collaboration. (XXXXa). Phenomenological MSSM interpretation of the CMS 2011 5fb-1 results. CMS-PAS-SUS-12-030. CMS Collaboration. (XXXXb). Phenomenological MSSM interpretation of the CMS 7 and 8 TeV results. CMS-PAS-SUS-13-020. CMS Collaboration. (XXXXc). Search for new physics in monojet events in pp collisions at sqrt(s)= 8 TeV. CMS-PAS-EXO-12-048. Conley, J. A., Gainer, J. S., Hewett, J. L., Le, M. P., & Rizzo, T. G. (XXXX). Supersymmetry without prejudice at the 7 TeV LHC. (arxiv:1103.1697[hep-ph]). Djouadi, A., et al. (MSSM Working Group Collaboration). (XXXX). The minimal supersymmetric standard model: Group summary report. hep-ph/9901246. Gelmini, G., Gondolo, P., Soldatenko, A., & Yaguna, C. E. (2006). The effect of a late decaying scalar on the neutralino relic density. Physical Review D, 74, 083514 (hep-ph/0605016). Kamionkowski, M., & Turner, M. S. (1990). Thermal relics: Do we know their abundances? Physical Review D, 42, 3310. Kowalska, K., Roszkowski, L., Sessolo, E. M., & Trojanowski, S. (2014). Low fine tuning in the MSSM with Higgsino Dark Matter and unification constraints. Journal of High Energy Physics, 1404, 166 (arxiv:1402.1328 [hep-ph]). Mahmoudi, F. (2008). SuperIso: A program for calculating the isospin asymmetry of B K* γ in the MSSM. Computer Physics Communications, 178, 745 (arxiv:0710.2067 [hep-ph]). Mahmoudi, F. (2009). SuperIso v2.3: A program for calculating flavor physics observables in supersymmetry. Computer Physics Communications, 180, 1579 (arxiv:0808.3144 [hep-ph]). Mahmoudi, F., Neshatpour, S., & Virto, J. (2014). B K* μ+ μ-optimised observables in the MSSM. European Physical Journal C, 74(6), 2927 (arxiv:1401.2145 [hep-ph]). Ovyn, S., Rouby, X., & Lemaitre, V. (XXXX). DELPHES, a framework for fast simulation of a generic collider experiment. (arxiv:0903.2225 [hep-ph]). Pumplin, J., Stump, D. R., Huston, J., Lai, H. L., Nadolsky, P. M., & Tung, W. K. (2002). New generation of parton distributions with uncertainties from global QCD analysis. Journal of High Energy Physics, 0207, 012 (hep-ph/0201195). Sekmen, S., Kraml, S., Lykken, J., Moortgat, F., Padhi, S., Pape, L.,... Spiropulu, M. (2012). Interpreting LHC SUSY searches in the phenomenological MSSM. Journal of High Energy Physics, 1202, 075 (arxiv:1109.5119 [hep-ph]). Page 9 of 10

Sjostrand, T., Mrenna, S., & Skands, P. Z. (2006). PYTHIA 6.4 physics and manual. Journal of High Energy Physics, 0605, 026 (hep-ph/0603175). Sjostrand, T., Mrenna, S., & Skands, P. Z. (2008). A brief introduction to PYTHIA 8.1. Computer Physics Communications, 178, 852 (arxiv:0710.3820 [hep-ph]). 2015 The Author(s). This open access article is distributed under a Creative Commons Attribution (CC-BY) 4.0 license. You are free to: Share copy and redistribute the material in any medium or format Adapt remix, transform, and build upon the material for any purpose, even commercially. The licensor cannot revoke these freedoms as long as you follow the license terms. Under the following terms: Attribution You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use. No additional restrictions You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits. Cogent Physics (ISSN: 2331-1940) is published by Cogent OA, part of Taylor & Francis Group. Publishing with Cogent OA ensures: Immediate, universal access to your article on publication High visibility and discoverability via the Cogent OA website as well as Taylor & Francis Online Download and citation statistics for your article Rapid online publication Input from, and dialog with, expert editors and editorial boards Retention of full copyright of your article Guaranteed legacy preservation of your article Discounts and waivers for authors in developing regions Submit your manuscript to a Cogent OA journal at www.cogentoa.com Page 10 of 10