# Global analysis of the pMSSM in light of the Fermi GeV excess: prospects for the LHC Run-II and astroparticle experiments

## Abstract

We present a new global fit of the 19-dimensional phenomenological Minimal Supersymmetric Standard Model (pMSSM-19) that comply with all the latest experimental results from dark matter indirect, direct and accelerator dark matter searches. We show that the model provides a satisfactory explanation of the excess of gamma-rays from the Galactic centre observed by the Fermi Large Area Telescope, assuming that it is produced by the annihilation of neutralinos in the Milky Way halo. We identify two regions that pass all the constraints: the first corresponds to neutralinos with a mass GeV annihilating into with a branching ratio of 95% ; the second to heavier neutralinos, with mass GeV annihilating into with a branching ratio of 87%. We show that neutralinos compatible with the Galactic centre GeV excess will soon be within the reach of LHC run-II – notably through searches for charginos and neutralinos, squarks and light smuons – and of Xenon1T, thanks to its unprecedented sensitivity to spin-dependent cross-section off neutrons.

a]Gianfranco Bertone, a]Francesca Calore, b]Sascha Caron, c]Roberto Ruiz, d]Jong Soo Kim, e]Roberto Trotta, a]Christoph Weniger.

[a]GRAPPA, University of Amsterdam, Science Park 904, 1090 GL Amsterdam, Netherlands \affiliation[b]IMAPP, Radboud University Nijmegen, P.O. Box 9010, NL-6500 GL Nijmegen, The Netherlands and Nikhef, Science Park, Amsterdam, The Netherlands \affiliation[c]Instituto de Física Corpuscular, IFIC-UV/CSIC, Valencia, Spain \affiliation[d]Instituto de Física Teórica UAM/CSIC, Madrid, Spain \affiliation[e]Imperial Centre for Inference and Cosmology, Imperial College London, Blackett Laboratory, Prince Consort Road, London SW7 2AZ, UK

## 1 Introduction

The existence of dark matter (DM) is now robustly established [1, 2, 3, 4] and its cosmological abundance measured with high precision [5]. Yet, the fundamental nature of the most abundant matter component in the Universe is unknown. According to the most promising theories, DM is a new fundamental particle. As a consequence, the search for DM is also a search for new physics beyond that of the well-known elementary particles, as laid down in the Standard Model (SM). Weakly interacting massive particles (WIMPs) are the leading DM candidates: they arise in many extensions of the Standard Model of particle physics, and naturally achieve the appropriate relic density through self-annihilation the early Universe. WIMPs can be searched for with three detection strategies: direct detection of the energy recoil of nuclei scattering off DM particles; indirect detection of the final stable products of DM annihilation or decay, as for example gamma rays; and accelerator searches for new particles, in particular at the Large Hadron Collider (LHC).

An excess in the -ray emission from the centre of our Galaxy has been discovered in data from the Large Area telescope (LAT), aboard the Fermi satellite [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. The nature of the so-called Fermi GeV excess remains a mystery. Several explanations have been put forward, the most exciting of which is perhaps DM annihilation in the halo of the Milky Way (see for example [17, 14, 13]). Astrophysical processes have also been suggested: the emission from a population of dim unresolved sources [18, 19, 20, 21], and the inverse Compton emission from a new population of cosmic-ray, either from time-dependent events taking place in the past of the Galaxy [22, 23, 24] or from the high star formation activity in the inner Galaxy [25].

Very recently, two reanalyses of the gamma-ray emission from the inner Galaxy found strong evidence for the Fermi GeV excess being due to hundreds to thousands of dim or unresolved point sources [26, 27], most likely millisecond pulsars [18]. Even more recently, it was shown that part or all of the required millisecond pulsar population could originate from the disruption of globular clusters by tidal forces in the inner Galaxy [28]. Directly detecting some of the millisecond pulsars in the inner Galaxy by radio observations is the next critical step for fully establishing this scenario [26].

It is however striking that the Fermi GeV excess spectrum and spatial distribution are well fitted by what is expected from DM annihilation. The excess could be the first non-gravitational signal of DM particles. It is thus urgent to either corroborate or disprove the DM particle nature of the Fermi GeV excess in the framework of concrete models for physics beyond the SM. Supersymmetry (SUSY) is one of the best-motivated classes of renormalizable extensions of the SM, which can accommodate a stable DM particle together with new degrees of freedom that mediate interactions. The most generic R-parity conserving and phenomenologically viable supersymmetric model is the phenomenological Minimal Supersymmetric Standard Model (pMSSM) [29]. In this work we address the following key question: can models of the pMSSM explain the observed properties of the Fermi GeV excess, while retaining consistency with other experimental data? And if so, what are the detection prospects for future direct detection and collider experiments?

Until now, MSSM scenarios [30] could not reproduce the Fermi GeV excess as observed by [14], as it was impossible to obtain in this framework a light neutralino ( 30 – 40 GeV), as required to fit the Fermi GeV excess spectrum, which could also account for the cosmological DM as measured by Planck [5]. However, [15] demonstrated that higher WIMP masses and annihilation channels different from b-quark pairs can give a good fit to the Fermi excess, owing to the freedom allowed by background model systematics. By fitting the GeV excess data of [15], it has been shown that viable solutions in the MSSM exist [31]. In the context of the pMSSM, [32] demonstrated that a re-assessment of the theoretical uncertainties in the DM signal spectra opens up a new phenomenology at the LHC experiments. Here we present the first systematic study of the pMSSM parameter space through global fitting techniques. This approach exhaustively covers all possible phenomenological signatures, allowing us a complete overview of the viable pMSSM interpretation of the Fermi GeV excess.

The paper is organized as follows: in section 2 we briefly present the well-known framework of the pMSSM and the parameters describing the model. In section 3 we describe the experimental set-up of the global fit and the implementation of the joint likelihood. We present the results of the parameter scan in section 4. We then discuss the implications for future direct and indirect detection experiments in section 5.1 and the prospects for detection at the LHC run II in section 5.2. Finally, we summarize our conclusions in section 6.

## 2 The theoretical framework: the pMSSM

We here study the pMSSM [29], in which the number of free parameters can be reduced to 19, given the present lack of experimental evidence for SUSY and no experimental indication that one requires the full freedom of a 22-dimensional pMSSM at present.

In this model, the lightest supersymmetric particle is the lightest neutralino, , a combination the neutral electroweak gauginos and higgsinos fields. The neutralino is one of the most well-motivated particle DM candidates since it is neutral, stable over cosmological timescales and can lead naturally to the correct DM relic abundance in the early Universe. In what follow, we assume that the neutralino is the particle DM candidate although we do not impose it to fully account for the DM relic abundance as measured by Planck but we allow for subdominant contributions to the DM content (cf. section 3).

We assume first and second generation mass universality, separately in the lepton and quark sectors (table 1). The trilinear couplings of the sfermions enter in the off-diagonal parts of the sfermion mass matrices. Since these entries are proportional to the Yukawa couplings of the respective fermions, we approximate the trilinear couplings associated with the first and second generation fermions to be zero, while the parameters , and represent the third generation trilinear couplings. In our set-up, the Higgs sector is fully described by the ratio of the Higgs vacuum expectation values , the higgsino mass parameter and the mass of the pseudoscalar Higgs , which are more directly related to the phenomenology of the model. This 19-dimensional realization of the pMSSM encapsulates all phenomenologically relevant features of the full model that are of interest for DM and collider experiments. The model parameters are displayed in table 1, along with their prior ranges. All of the input parameters are defined at the SUSY scale .

pMSSM parameters and priors | |||
---|---|---|---|

Flat priors | Log priors | ||

[TeV] | (-5, 5) | ||

[TeV] | (0.1, 5) | ||

[TeV] | (-5, 5) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (0.1,4) | ||

[TeV] | (-10, 10) | ||

[TeV] | (-10,10) | ||

[TeV] | (-10,10) | ||

[TeV] | (-5,5) | ||

[TeV] | (0.01, 5) | ||

[GeV] | [33] (Gaussian prior) | ||

[GeV/cm] | [34] (Gaussian prior) |

## 3 The experimental setup

We implement experimental constraints with a joint likelihood function, whose logarithm takes the following form:

(1) | ||||

where represents the Fermi GeV excess, electroweak precision observables, B and D physics constraints, measurements of the cosmological DM relic density, () direct (indirect) DM detection constraints and () Higgs (sparticles) searches at colliders. We discuss each component in turn:

#### :

For the Fermi GeV excess likelihood we follow the treatment in [15], to account astrophysical uncertainties. In addition, we include uncorrelated 10% uncertainties as DM modeling systematics, following ref. [32]. We marginalize the likelihood function over the uncertainties in the Galactic centre -value, which is assumed to follow a log-normal distribution with mean and variance [17]. When the predicted neutralino relic density, , is smaller than the Planck measurement, , we follow [35] and adopt the so-called “scaling Ansatz”, i.e. we assume that the local ratio of neutralino () to total DM densities () is equal to that for their cosmic abundances:

(2) |

#### :

This implements constraints from Z-pole measurements at LEP [36]. We include the constraint on the effective electroweak mixing angle for leptons , the total width of the Z-boson , the hadronic pole cross-section , as well as the decay width ratios , and the asymmetry parameters , , and , . In addition, we also use the measurement of the mass of the W boson from the LEP experiment [36]. We apply a Gaussian likelihood for all of these quantities, with mean and standard deviation as given in table II of [37].

#### :

The flavor observables related with B and D physics considered are , , , , , , and . For all we apply a Gaussian likelihood and for most of them we use the measurements following table II of [37]. The experimental values assumed for and are and (theoretical uncertainties included [38]) [39] respectively.

#### :

We include the Planck Cosmic Microwave Background data constraint on the DM relic abundance as an upper limit, to allow for the possibility that neutralinos are a sub-dominant DM component. We follow the formalism in the Appendix of [40], using as central value the result from Planck temperature and lensing data [41] with a (fixed) theoretical uncertainty, , to account for the numerical uncertainties entering in the calculation of the relic density.

#### :

#### :

This implements conservative upper limits on the proton spin-dependent cross-section from the IceCube detector in its 79-string configuration [46] (IC-79). Comparable – if slightly weaker – limits have been set for the channel by Super-Kamiokande [47] and ANTARES [48]. The most stringent constraint is for the case where WIMPs annihilate exclusively to pairs. Since the neutrino spectrum generated by bosons is similar to the bosons we apply this constraint whenever the combined branching fraction to and is above 95%. In that case the likelihood is a step function smeared with half a Gaussian (as in Eq. (3.5) of [49]) to account for theoretical and experimental uncertainties that we set to be 50% of the predicted value.

#### :

The likelihood for the Higgs searches has two components: the first implements bounds obtained from Higgs searches at LEP, Tevatron and LHC via HiggsBounds [50], which returns whether a model is excluded or not at the 95% CL. The second component constrains the mass and the production times decay rates of the Higgs-like boson discovered by the LHC experiments ATLAS [51] and CMS [52]. We use HiggsSignals [53] assuming a theoretical uncertainty in the lightest Higgs mass calculation of 2 GeV.

#### :

SUSY searches constraints at LEP and Tevatron follow the likelihood used [49]. We have imposed the strict constraints from a large number of searches for supersymmetry at the LHC. The branching ratios of the sparticles have been calculated with SUSYHIT 1.5 [54]. We have generated the hadronized event samples with Pythia 8.2 [55] and have employed the NNPDF 2.3 parton distribution functions [56]. The generated events are passed on to CheckMATE [57], which is based on the fast detector simulation Delphes 3.10 [58]. CheckMATE tests if the model point in question is excluded or not at confidence level by comparing to current experimental searches at the LHC for supersymmetry in the relevant hadronic and leptonic final states with large missing transverse momentum. We assign a 0 log-likelihood if the point passes all constraints, and exclude it if it fails any of them.

We have only included observables we consider robust in order to be conservative. For instance, we have dropped the electroweak precision observables and in the fit because it is unclear whether the large deviations of 2.5 that are observed with respect to the SM predictions are due to unknown systematic uncertainties or to new physics. The experimental status of the magnetic anomaly of the muon, remains unclear in the face of persistent discrepancies in the determination of the hadronic vacuum-polarization diagram using either or the hadronic decay data. We do not include DM searches in dwarf spheroidal galaxies by the Fermi-LAT [59], as constraints are given only for 100% branching ratio into final states.

Finally, we refer to [37] for details about how the SUSY spectrum and observables are computed.

We use the MultiNest [60] algorithm as implemented in SuperBayeS-v2.0, to perform a global fit of the pMSSM parameter space, including all the data in eq. (1), excepting the SUSY searches at the LHC. This is because the LHC searches evaluation is computationally too expensive to be performed on-the-fly. Our scans were run using both log and flat priors to ensure a complete coverage of the parameter space, gathering samples from likelihood evaluations. Samples have been thinned by a factor of 10, focusing our search to regions of the parameter space that were not clearly ruled out by LHC-Run I constraints. This produced representative samples to which the LHC SUSY searches have been applied. The ensuing samples that pass LHC run I constraints are displayed in figure 1.

## 4 Results

Our global fits identify two distinct viable solutions in the pMSSM parameter
space (figure 1): the first exhibits a WIMP mass of GeV, with the neutralino
annihilating to with a 95% branching ratio. The second solution has a
larger neutralino mass, GeV, and 87% annihilation
final states. The overall best-fit point is in the region has . This is for a fit with 21 free parameters, and 125
Gaussian data points (we do not include limits as their is normalized
to 0 whenever the constraint is satisfied), so we adopt 104 degrees of freedom^{1}

It is important to notice that including theoretical uncertainties in the GC fit is crucial in achieving reasonable -values. The 10% theoretical uncertainty advocated in [32] is a reasonable reflection of the differences in the predicted spectra between current numerical codes. However, in absence of such an uncertainty, the quality of our global fit would degrade to 0.023 and 0.008 for the two best-fit points, respectively.

The contribution to the overall for the two best-fit points from different types of observables are plotted in figure 2. The pulls have been normalized by the number of data points in each group, , to facilitate a visual comparison. We notice that the per data point is distributed fairly evenly across observables. There is a slight preponderance in the contribution coming from the GC Fermi fit (with bins), which is exacerbated if one neglects the 10% theoretical uncertainty in the DM spectra (dashed bars). The contribution to the pull from the LUX likelihood comes almost exclusively from the SD neutron cross-section limit, as the SI constraint is easily satisfied by our best-fit points.

Let us now analyze in more detail both type of solutions.

In the region, model points providing a better fit have a neutralino mostly bino-like ( 80 – 90%) with a similar fraction of both wino and higgsino. Besides we find points in which neutralinos can be dominantly higgsinos with a bino fraction as small as 10%. Those provide a worse fit though, because a large higgsino composition, basically, implies a large annihilation cross section which drops the relic density below the Planck limit leading to a tension with the Fermi GeV excess due to the scaling Ansatz given in eq. (2).

Analogously to Ref. [32], also in the present work, we find a WW solution with a bino-higgsino neutralino composition of about 84 – 92 GeV. However, in our work, this solution has a slightly worse than in Ref. [32]. The reason for this is that we adopt slightly different values of the form-factors for the computation of the SD cross-section. Since such a WW solution is now just below the WW IceCube limits, it is punished a bit in the likelihood.

Typically, since third generation squarks enter into the SI cross section at loop level through LSP-gluon effective interactions [62], their contribution can be comparable to tree-level effective interactions mediated by squarks of the first two families when they are light. In this type of solutions, because current constraints allow sbottoms of a few hundred GeV, their contribution can be sizable and, indeed, cancels out the Higgs exchange contributions. This effect allows to relax the tension with LUX data when Higgs exchange contributions are large.

In terms of the impact of LHC Run-I data, we notice that first and second generation squarks as well as gluinos are decoupled. The stops and the heavier sbottom mass eigenstates are also kinematically inaccessible at the LHC Run-I energies. However, the lighter sbottom eigenstate with a mass around 400 GeV can be produced at the LHC at a considerable rate but the sbottoms evade detection from third generation searches due to complicated cascade decays. The sparticles in the electroweak sector are relatively light. There is a large mass splitting between the SU(2) doublet and singlet sleptons. The SU(2) doublet sleptons are light with masses around 250 GeV and narrowly evade detection in searches for two lepton and large missing transverse momentum final states. The lightest wino-like chargino and the second lightest neutralino escape detection since they are almost mass degenerate with the bino-like neutralino. The production rate of the higgsino eigenstates is too small to yield an observable signal at the LHC run I.

Results about the annihilation cross-section are shown in the top left panel of figure 1. One can see that the points with a better fit exhibit cm/s, consistent with the results found in [15]. We also show the Fermi dwarfs limit for the channel [59], but we emphasize that this limit has not been applied in the fit. In order to compare with the constraint coming from Fermi dwarfs observations, the annihilation cross-section needs to be rescaled by a factor (to account for the possibility of sub-dominant dark matter relic density, which translates according to our Ansatz in a correspondingly reduced local density in the dwarfs), which would suppress the signal well below current and future limits from dwarfs.

Regarding DM direct detection, the top right panel of figure 1 shows the spin-independent (SI) cross-section versus neutralino mass plane. In order to facilitate the visual comparison of our pMSSM models with existing and future limits, we have rescaled the theoretical cross-section by a factor (to account for models where the neutralino does not make up all of the cosmological relic density) and a factor GeV/cm. This accounts for the fact that the local density we have used to predict the number of counts for LUX, , is a nuisance parameter which is generally different from the value assumed by the LUX collaboration in deriving their limit [42], namely 0.3 GeV/cm.

The points that appear above the nominal 95% exclusion limit from LUX cannot be excluded because of a combination of effects: (i) our LUX likelihood is slightly less stringent than what has been published by the LUX collaboration (and depicted in figure 1), (ii) our likelihood function for LUX allows for values of the cross section above the 95% limit to be included (albeit penalized by a smaller likelihood value), and (iii) the global likelihood function can – to an extent – compensate for a poor fit to LUX data by gaining an improvement from other data sets.

Since squarks are heavy, the contribution coming from the exchange of a CP-even Higgs dominates and therefore the SI cross-section scales as . Because, as noticed above, the higgsino fraction in this region is not negligible, it can be large. Indeed, there are model points above the LUX limit allowed by the scaling Ansatz we apply to the local dark matter density and to some extend due to the fact that we vary the local dark matter density. Another interesting feature is that the SI cross-section spans down to pb. That is possible because the heavy CP-even Higgs contribution can be sizable and cancellations with the lightest Higgs channel might occur [63].

In the bottom panels of figure 1 we display the spin dependent (SD) cross-section for scattering off neutrons/protons (left/right panels). While in the case of SI interactions the contributions for proton and neutron are comparable, the SD cross-sections may differ significantly. However, we find a tight correlation between the SD cross-sections for scattering off neutrons and those off protons in our model points. The SD cross-section is dominated by the exchange of a Z boson and therefore the SD cross-section is largely determined by the higgsino content of the neutralino and likewise for the SI cross-section it can be sizable in this region as it can be seen in both panels.

In terms of dark matter direct detection experiments, at present, LUX data represents the strongest constraint on the SD-neutron scattering cross-section because the Xenon contains neutron-odd isotopes therefore we overlay the LUX constraint properly rescaled as for the SI case. For the SD-proton scattering cross-section, IC-79 represents the strongest current constraint for the particular case when the neutralinos annihilate to a WW final state so we show the IC-79 90% CL upper limit [46] , and we have rescaled the value of the SD cross-section by a factor . This assumes that equilibrium between capture and annihilation is reached in the Sun, which is a good assumption for the bulk of the models shown here. In fact, in the region where the neutralino annihilates mainly to and to to a lesser extend, the IC-79 limits apply here and disfavors a large number of model points.

One can see that there are model points above LUX and IC exclusion lines. In those, the higgsino component of the neutralino is dominant over the gaugino one leading to a large SD cross-section for both neutrons and protons. Those points still provide a reasonable fit to the data due to two effects, first because neutralinos with a large higgsino component yields to a relic density sensible below the Planck bound and therefore the scaling Ansatz applies and second because the local dark matter density is a nuisance parameter in our analysis. Beside there is a small fraction of points above the nominal IC-79 limit with a branching ratio to and low enough to evade the stringent IC-79 bound.

In the region, points providing a better fit have a neutralino dominantly bino-like () with a of higgsino. Those points have the characteristic that the neutralinos annihilate to top quark pairs via an exchange of a right-handed stop which is relatively heavy ( 1 TeV). This is possible because on top of the non-helicity suppression the neutralino-stop-top coupling component, which is proportional to the top quark Yukawa coupling, is sizeable due to the non-vanishing higgssino fraction of the neutralino. We also found points in which the right-handed stops are light ( 300 GeV) being almost bino-like (also found in [32]). However those provide a worse fit to the data because higgssino-stop loops have a sizable positive contribution which leads to values above the experimental constraint.

Annihilation into through a t-channel stop exchange requires, for bino-like neutralinos, stops with masses of a few hundred GeV. On the other hand, relatively light stops are not able to lift the tree-level Higgs mass to fit a 125 GeV Higgs. In order to enhance the DM annihilation cross sections to match the Planck measurement, the neutralino coannihilates with sneutrinos which have to be about the same mass as the DM particle ( 200 GeV). This induces the splitting in the left/right sleptons spectrum.

These benchmark points are characterized by left handed slepton next-to-lightest supersymmetric particle with masses above 200 GeV. The higgsinos have masses around 260 GeV. Since the resulting mass differences between the bino-like LSP and the sleptons/higgsinos are small, the final states are rather soft and thus the detection is suppressed in events with dilepton or trilepton final states and large missing transverse energy momentum. The production rate of the sbottom is quite suppressed and hence avoids detection. The spectrum of the remaining supersymmetric particles is decoupled.

The phenomenology of the model points in the region in terms of DM detection is similar to the one. The main difference is in that the IC-79 limits does not apply and therefore larger SD cross-sections are possible. The main difference is that the higgsino composition is not as large as in the WW type of solutions and therefore the SD cross sections for neutron interactions are below the LUX current sensitivity.

Lastly, in figure 3, we show the spectrum of the Fermi GeV excess together with the systematic uncertainties associated with the galactic diffuse emission modeling [15]. We compare the data with the spectra of the pMSSM model points giving the best global -value in the two regions identified in Figure 1.

It is apparent from Figure 3 that the best-fit DM spectra are systematically offset from the mean values of the Galactic center excess spectrum (gray dots and boxes), by about 1 to 2 sigma, and do not provide a good fit to the data at first sight. However, since the systematic astrophysical uncertainties, indicated by the gray boxes () are correlated, this still provides an acceptable fit to the data. To illustrate this point, we show with black dots and error bars the excess spectrum where we moved all data points systematically down, according to the freedom allowed by the covariance (the error bars show now statistical errors only). Together with the 10% uncorrelated systematic modeling uncertainty that we adopted for the DM signal, this provides a reasonable fit to the data, with -values, whereas without the DM signal modeling uncertainties, the -values would be prohibitively small (see Figure 2 above).

## 5 Discovery potential: dark matter detection experiments and LHC run II

### 5.1 Implications for direct and indirect dark matter searches

Generally, spin-dependent and spin-independent scattering cross sections are driven by
the higgsino content of the neutralino.
Therefore, as explained above, the sizable higgsino fraction of points in the two best-fit regions imply large SI cross-sections
that makes the direct detection prospects promising, although the SI cross-section range spans down to
pb due to cancellations with the lightest Higgs.
In figure 1, top right panel, we also display the projected sensitivity limit (defined as the 90% CL exclusion limit) for the Xenon1T experiment and an
hypothetical liquid Xe detector with 500 t yr exposure [61]. The latter experiment essentially
saturates the ultimate detection floor set by coherent neutrino scattering [64].
Xenon1T data will be crucial in discovering or firmly ruling out models belonging to the island and will
probe a significant fraction of the parameter space preferred by the first type of solutions.
In figure 1, bottom left panel, we overlay the projected 90% exclusion limit for Xenon1T [61].
We find that Xenon1T will be able to prove the entirety of the SD neutron scattering cross-section parameter space favoured by our models. ^{2}

The bottom right panel of figure 1 shows that our best-fit points easily evade the constraint set by IC-79 on the SD proton cross-section in the implementation we adopted in this paper. However, an event-level implementation of the likelihood (including the events’ energies [66]) would increase the constraining power of the IC-79 limit, to the point that some of the surviving models could be probed [67].

Finally, as for indirect detection, the preferred parameter space is mostly out of reach even for the future 10-yr Fermi analysis of dwarf spheroidal galaxies [68], where an further improvement of current sensitivities by a factor 2 – 3 can be expected. Note that points with a large annihilation cross section as shown in the top right panel of figure 1 usually correspond to suppressed relic densities, making these points hard to detect due to the factor suppression of their signal.

### 5.2 Prospects for detection at the LHC run II

The models yielding the best -values show several interesting properties relevant for LHC searches. In the following we briefly discuss the discovery prospects depending on the produced SUSY particles. Sparticle mass distributions and MSSM parameters relevant for LHC searches are shown in figure 4.

#### Light squarks:

In the top left panel of Fig. 4 we plot the lightest first/second generation squark mass vs the neutralino mass.
The first and second generation squarks have masses GeV, i.e. above
the usual run I constraints. About of the
models have squark masses below GeV
^{3}

#### Stop:

In the top right panel of Fig. 4 we plot the stop mass vs the neutralino mass. Some models have light stops with masses down to 200 – 300 GeV decaying to chargino and a b-jet. The neutralino has a mass of around 95 GeV. These models are not excluded by current LHC searches [69]. Another interesting region also found in [32] has a stop mass of around 200 – 220 GeV and a mass of the lightest neutralino around 180 GeV. A slight excess in the ATLAS data prevents exclusion with run-1 data in this region [69]. These solutions will likely be tested with early run II analyses. Other solutions yield much heavier stop masses decaying predominantly into the heavier neutralino and chargino states. Dedicated searches for such decays are important.

#### Sbottom:

In the central left panel of Fig. 4 we plot the sbottom mass vs the neutralino mass. Several model points have a sbottom mass as low as 400 GeV. The points are not excluded in our procedure due to multi-step cascade decays involving heavy neutralinos.

Typically, the lighter sbottom state typically has masses around 400 GeV while the lightest neutralino is mostly bino-like with a mass around 90 GeV. The second lightest wino-like neutralino lies around 107 GeV. The two heavier neutralinos are higgsino dominated states and have masses around 250 GeV.

If the sbottoms predominantly decayed into a bottom quark and the lightest neutralino, these benchmark points would clearly be excluded. However, the bottom squark decays into all four neutralino as well as the lighter chargino mass eigenstates with comparable rates. The decay modes of the neutralino and chargino mass eigenstates are relatively complex and hence the limits from simplified sbottom searches do not apply. For instance, the second lightest neutralino eigenstate has large hadronic three body decay modes into the lightest neutralino via off shell Z bosons. Moreover, the second lightest neutralino can radiatively decay into a photon and the lightest neutralino. The corresponding lighter chargino eigentstate has relatively large leptonic three body decay modes via off shell W bosons. Finally, the two heaviest higgsino dominated neutralino mass eigenstates mainly decay into electroweak gauge bosons and the lighter electroweakino states. As a consequence, many events have leptons in the final state which are vetoed in the searches for direct sbottom production. In addition, the higher final state multiplicity tends to soften the net missing transverse momentum distribution compared to scenarios with direct sbottom decays into a bottom quark and the LSP. As a result fewer events pass the selection cuts of the relevant sbottom searches. We explicitly tested those light sbottom scenarios against experimental searches at the LHC with the computer tool CheckMATE and confirmed that those points were allowed.

As seen also for the light squarks about 70% of the models would be excluded with a sbottom limit of GeV.

#### Gluinos:

In the central right panel of Fig. 4 we plot the sbottom mass vs the neutralino mass. Gluinos have masses GeV. An upcoming early run II exclusion on gluinos with a mass up to 2500 GeV would exclude about of these models. Run II searches for first and second generation squarks will thus likely be more constraining.

#### Sleptons:

In the bottom left panel of Fig. 4 we plot the lightest smuon mass vs the neutralino mass. The lightest smuons found in the best fit models have masses GeV in about of the best fit points. This makes searches for smuons in run II very sensitive to these solutions.

#### Chargino/neutralino:

In the bottom right panel of Fig. 4 we plot the chargino mass vs the neutralino mass. Due to the GeV excess likelihood several neutralinos and charginos are typically light, the lightest having a mass fixed to GeV for the solutions and GeV for the solutions. The higgsino component in the and solutions typically involves that is only slightly larger than these mass scales, leading to 2 more neutralinos and the light chargino at masses around 100 GeV () or 200 – 300 GeV (). These states are often mass compressed with the lightest neutralino which makes the solutions evade LHC chargino/neutralino searches so far. Dedicated chargino/neutralino searches will have sensitivity to most models, e.g. by a mono-jet and soft lepton search as proposed also in ref. [32]. The Wino mass scale is quite unconstrained and lies between 100 GeV and 1.5 TeV. The Wino will decay to lighter states yielding final states with Z, W or Higgs bosons.

#### Heavy Higgs:

About of the best models have GeV making searches for heavy Higgs bosons very sensitive. Several chargino and neutralino states are light and have a large coupling to . Consequently heavy Higgs decays to charginos and neutralinos can have huge branching ratios up to competing with top and bottom decays. Dedicated searches for heavy Higgs bosons decaying into final states with with missing transverse momentum would help to constrain these scenarios.

## 6 Conclusions

In this paper we addressed the issue of finding model points in the pMSSM that might explain simultaneously the large set of independent data we gathered from astrophysics, cosmology and high-energy particle physics. We showed that no tension exists between currently available particle physics constraints and the interpretation of the Fermi GeV excess in terms of dark matter annihilation in the framework of the pMSSM. Furthermore, we found evidence for two regions that are able to explain the gamma-ray data, while being consistent with other various experimental constraints: (a) a first region where the neutralino is mostly bino-like and the dominant annihilation channel today is 95% into bosons pairs and (b) a second region where the the annihilation into top-quark pairs dominates and the neutralino is again mainly bino-like. We showed that these models are very appealing since they will be soon in the reach of the next generation of direct detection experiments – Xenon1T will probe the entirety of the best-fit regions thanks to its sensitivity to spin-dependent neutron cross-section – and of the LHC run II, in particular through searches for charginos and neutralinos, squarks and light smuons.

Acknowledgments. We warmly thank T. Bringmann, Marco Selvi, Christopher McCabe, Pat Scott and Tim Stefaniak for useful discussions, Marc Schumann and Marco Selvi for providing us with direct detection limits and William Parker for assistance with some figures. G.B. (P.I.) and F.C. acknowledge support from the European Research Council through the ERC starting grant WIMPs Kairos. R. RdA, is supported by the Ramón y Cajal program of the Spanish MICINN and also thanks the support of the Spanish MICINN’s Consolider-Ingenio 2010 Programme under the grant MULTIDARK CSD2209-00064, the Invisibles European ITN project (FP7-PEOPLE-2011-ITN, PITN-GA-2011-289442-INVISIBLES and the “SOM Sabor y origen de la Materia” (FPA2011-29678) and the “Fenomenologia y Cosmologia de la Fisica mas alla del Modelo Estandar e lmplicaciones Experimentales en la era del LHC” (FPA2010-17747) MEC projects. The work of J. S. Kim has been partially supported by the MINECO, Spain, under contract FPA2013-44773-P; Consolider-Ingenio CPAN CSD2007-00042 and the Spanish MINECO Centro de excelencia Severo Ochoa Program under grant SEV-2012-0249. R.T. acknowledges partial support from an EPSRC “Pathways to Impact” grant. C.W. was supported by the Netherlands Organization for Scientific Research (NWO) through a Vidi grant. We gratefully acknowledge the use of the Cartesius supercomputer (Amsterdam) and of Imperial College London’s HPC. This work was supported by Grant ST/N000838/1 from the Science and Technology Facilities Council.

### Footnotes

- We emphasize that the calculation of the number of actual degrees of freedom is not trivial. One would have to consider the number of active data points, as well as the number of effective parameters in the model. This can only be done properly via extensive Monte Carlo simulations of the data. The simple counting argument we adopt is meant to be representative of what one would get in the simplest scenario.
- We note that the Xenon1T exclusion limits in [61] are obtained by applying a scaling factor derived from the comparison between SI and SD results of Xenon100 and result to be slightly stronger than the limits quoted in ref. [65], which instead adopts a 60 times improvement w.r.to XENON100 (still Xenon1T will be able to entirely probe our identified parameter space).
- The expected reach of the ATLAS and CMS experiments for full MSSM models is not estimated. We compare with mass scales which are likely excluded with HL-LHC data for many MSSM scenarios.

### References

- G. Jungman, M. Kamionkowski, and K. Griest, Supersymmetric dark matter, Phys. Rept. 267 (1996) 195–373, [hep-ph/9506380].
- L. Bergstrom, Nonbaryonic dark matter: Observational evidence and detection methods, Rept. Prog. Phys. 63 (2000) 793, [hep-ph/0002126].
- G. Bertone, D. Hooper, and J. Silk, Particle dark matter: Evidence, candidates and constraints, Phys. Rept. 405 (2005) 279–390, [hep-ph/0404175].
- G. Bertone, ed., Particle Dark Matter: Observations, Models and Searches. Cambridge University Press, 2010.
- The Planck Collaboration, The Scientific Programme of Planck, ArXiv Astrophysics e-prints (Apr., 2006) [astro-ph/0].
- L. Goodenough and D. Hooper, Possible Evidence For Dark Matter Annihilation In The Inner Milky Way From The Fermi Gamma Ray Space Telescope, arXiv:0910.2998.
- Fermi/LAT Collaboration Collaboration, V. Vitale and A. Morselli, Indirect Search for Dark Matter from the center of the Milky Way with the Fermi-Large Area Telescope, arXiv:0912.3828.
- D. Hooper and L. Goodenough, Dark Matter Annihilation in The Galactic Center As Seen by the Fermi Gamma Ray Space Telescope, Phys.Lett. B697 (2011) 412–428, [arXiv:1010.2752].
- D. Hooper and T. Linden, On The Origin Of The Gamma Rays From The Galactic Center, Phys.Rev. D84 (2011) 123005, [arXiv:1110.0006].
- K. N. Abazajian and M. Kaplinghat, Detection of a Gamma-Ray Source in the Galactic Center Consistent with Extended Emission from Dark Matter Annihilation and Concentrated Astrophysical Emission, Phys.Rev. D86 (2012) 083511, [arXiv:1207.6047].
- C. Gordon and O. Macias, Dark Matter and Pulsar Model Constraints from Galactic Center Fermi-LAT Gamma Ray Observations, Phys.Rev. D88 (2013) 083521, [arXiv:1306.5725].
- D. Hooper and T. R. Slatyer, Two Emission Mechanisms in the Fermi Bubbles: A Possible Signal of Annihilating Dark Matter, Phys.Dark Univ. 2 (2013) 118–138, [arXiv:1302.6589].
- K. N. Abazajian, N. Canac, S. Horiuchi, and M. Kaplinghat, Astrophysical and Dark Matter Interpretations of Extended Gamma Ray Emission from the Galactic Center, arXiv:1402.4090.
- T. Daylan, D. P. Finkbeiner, D. Hooper, T. Linden, S. K. N. Portillo, et al., The Characterization of the Gamma-Ray Signal from the Central Milky Way: A Compelling Case for Annihilating Dark Matter, arXiv:1402.6703.
- F. Calore, I. Cholis, and C. Weniger, Background model systematics for the Fermi GeV excess, arXiv:1409.0042.
- Fermi-LAT Collaboration, S. Murgia Talk given at the 2014 Fermi Symposium, Nagoya, Japan, October 20-24 (2014).
- F. Calore, I. Cholis, C. McCabe, and C. Weniger, A Tale of Tails: Dark Matter Interpretations of the Fermi GeV Excess in Light of Background Model Systematics, arXiv:1411.4647.
- K. N. Abazajian, The Consistency of Fermi-LAT Observations of the Galactic Center with a Millisecond Pulsar Population in the Central Stellar Cluster, JCAP 1103 (2011) 010, [arXiv:1011.4275].
- F. Calore, M. Di Mauro, F. Donato, and F. Donato, Diffuse gamma-ray emission from galactic pulsars, Astrophys.J. 796 (2014) 1, [arXiv:1406.2706].
- Q. Yuan and B. Zhang, Millisecond pulsar interpretation of the Galactic center gamma-ray excess, arXiv:1404.2318.
- J. Petrovic, P. D. Serpico, and G. Zaharijas, Millisecond pulsars and the Galactic Center gamma-ray excess: the importance of luminosity function and secondary emission, arXiv:1411.2980.
- E. Carlson and S. Profumo, Cosmic Ray Protons in the Inner Galaxy and the Galactic Center Gamma-Ray Excess, Phys.Rev. D90 (2014) 023015, [arXiv:1405.7685].
- J. Petrovic, P. D. Serpico, and G. Zaharijas, Galactic Center gamma-ray ”excess” from an active past of the Galactic Centre?, JCAP 1410 (2014), no. 10 052, [arXiv:1405.7928].
- I. Cholis, C. Evoli, F. Calore, T. Linden, C. Weniger, et al., The Galactic Center GeV Excess from a Series of Leptonic Cosmic-Ray Outbursts, arXiv:1506.05119.
- D. Gaggero, M. Taoso, A. Urbano, M. Valli, and P. Ullio, Towards a realistic astrophysical interpretation of the Galactic center excess, arXiv:1507.06129.
- R. Bartels, S. Krishnamurthy, and C. Weniger, Strong support for the millisecond pulsar origin of the Galactic center GeV excess, arXiv:1506.05104.
- S. K. Lee, M. Lisanti, B. R. Safdi, T. R. Slatyer, and W. Xue, Evidence for Unresolved Gamma-Ray Point Sources in the Inner Galaxy, arXiv:1506.05124.
- T. D. Brandt and B. Kocsis, Disrupted Globular Clusters Can Explain the Galactic Center Gamma Ray Excess, arXiv:1507.05616.
- MSSM Working Group Collaboration, A. Djouadi et al., The Minimal supersymmetric standard model: Group summary report, hep-ph/9901246.
- M. Cahill-Rowley, J. Gainer, J. Hewett, and T. Rizzo, Towards a Supersymmetric Description of the Fermi Galactic Center Excess, JHEP 1502 (2015) 057, [arXiv:1409.1573].
- P. Agrawal, B. Batell, P. J. Fox, and R. Harnik, WIMPs at the Galactic Center, arXiv:1411.2592.
- A. Achterberg, S. Caron, L. Hendriks, R. Ruiz de Austri, and C. Weniger, A description of the Galactic Center excess in the Minimal Supersymmetric Standard Model, arXiv:1502.05703.
- ATLAS, CDF, CMS, D0 Collaboration, First combination of Tevatron and LHC measurements of the top-quark mass, arXiv:1403.4427.
- M. Pato, L. Baudis, G. Bertone, R. Ruiz de Austri, L. E. Strigari, et al., Complementarity of Dark Matter Direct Detection Targets, Phys.Rev. D83 (2011) 083505, [arXiv:1012.3458].
- G. Bertone, D. G. Cerdeno, M. Fornasa, R. R. de Austri, and R. Trotta, Identification of Dark Matter particles with LHC and direct detection data, Phys.Rev. D82 (2010) 055008, [arXiv:1005.4280].
- ALEPH, DELPHI, L3, OPAL, SLD, LEP Electroweak Working Group, SLD Electroweak Group, SLD Heavy Flavour Group Collaboration, S. Schael et al., Precision electroweak measurements on the resonance, Phys.Rept. 427 (2006) 257–454, [hep-ex/0509008].
- C. Strege, G. Bertone, G. Besjes, S. Caron, R. Ruiz de Austri, et al., Profile likelihood maps of a 15-dimensional MSSM, JHEP 1409 (2014) 081, [arXiv:1405.0622].
- A. Arbey, M. Battaglia, F. Mahmoudi, and D. MartÃnez Santos, Supersymmetry confronts : Present and future status, Phys.Rev. D87 (2013), no. 3 035026, [arXiv:1212.4887].
- CMS, LHCb Collaboration, CMS and L. Collaborations, Combination of results on the rare decays from the CMS and LHCb experiments, .
- G. Bertone, K. Kong, R. R. de Austri, and R. Trotta, Global fits of the Minimal Universal Extra Dimensions scenario, Phys.Rev. D83 (2011) 036008, [arXiv:1010.2023].
- Planck Collaboration, P. Ade et al., Planck 2013 results. XVI. Cosmological parameters, Astron.Astrophys. 571 (2014) A16, [arXiv:1303.5076].
- LUX Collaboration, D. Akerib et al., First results from the LUX dark matter experiment at the Sanford Underground Research Facility, Phys.Rev.Lett. 112 (2014) 091303, [arXiv:1310.8214].
- C. Savage, A. Scaffidi, M. White, and A. G. Williams, LUX likelihood and limits on spin-independent and spin-dependent WIMP couplings with LUXCalc, arXiv:1502.02667.
- QCDSF Collaboration, G. S. Bali et al., Strangeness Contribution to the Proton Spin from Lattice QCD, Phys.Rev.Lett. 108 (2012) 222001, [arXiv:1112.3354].
- P. Junnarkar and A. Walker-Loud, Scalar strange content of the nucleon from lattice QCD, Phys.Rev. D87 (2013) 114510, [arXiv:1301.1114].
- IceCube Collaboration, M. Aartsen et al., Search for dark matter annihilations in the Sun with the 79-string IceCube detector, Phys.Rev.Lett. 110 (2013), no. 13 131302, [arXiv:1212.4097].
- Super-Kamiokande Collaboration, K. Choi et al., Search for neutrinos from annihilation of captured low-mass dark matter particles in the sun by super-kamiokande, Phys. Rev. Lett. 114 (Apr, 2015) 141301.
- ANTARES Collaboration, S. Adrian-Martinez et al., First results on dark matter annihilation in the Sun using the ANTARES neutrino telescope, JCAP 1311 (2013) 032, [arXiv:1302.6516].
- R. R. de Austri, R. Trotta, and L. Roszkowski, A Markov chain Monte Carlo analysis of the CMSSM, JHEP 0605 (2006) 002, [hep-ph/0602028].
- P. Bechtle, O. Brein, S. Heinemeyer, O. StÃ¥l, T. Stefaniak, et al., : Improved Tests of Extended Higgs Sectors against Exclusion Bounds from LEP, the Tevatron and the LHC, Eur.Phys.J. C74 (2014), no. 3 2693, [arXiv:1311.0055].
- ATLAS Collaboration, G. Aad et al., Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC, Phys.Lett. B716 (2012) 1–29, [arXiv:1207.7214].
- CMS Collaboration, S. Chatrchyan et al., Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC, Phys.Lett. B716 (2012) 30–61, [arXiv:1207.7235].
- P. Bechtle, S. Heinemeyer, O. StÃ¥l, T. Stefaniak, and G. Weiglein, : Confronting arbitrary Higgs sectors with measurements at the Tevatron and the LHC, Eur.Phys.J. C74 (2014), no. 2 2711, [arXiv:1305.1933].
- A. Djouadi, M. Muhlleitner, and M. Spira, Decays of supersymmetric particles: The Program SUSY-HIT (SUspect-SdecaY-Hdecay-InTerface), Acta Phys.Polon. B38 (2007) 635–644, [hep-ph/0609292].
- T. SjÃ¶strand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, et al., An Introduction to PYTHIA 8.2, Comput.Phys.Commun. 191 (2015) 159–177, [arXiv:1410.3012].
- R. D. Ball, V. Bertone, S. Carrazza, C. S. Deans, L. Del Debbio, et al., Parton distributions with LHC data, Nucl.Phys. B867 (2013) 244–289, [arXiv:1207.1303].
- M. Drees, H. Dreiner, D. Schmeier, J. Tattersall, and J. S. Kim, CheckMATE: Confronting your Favourite New Physics Model with LHC Data, Comput.Phys.Commun. 187 (2014) 227–265, [arXiv:1312.2591].
- DELPHES 3 Collaboration, J. de Favereau et al., DELPHES 3, A modular framework for fast simulation of a generic collider experiment, JHEP 1402 (2014) 057, [arXiv:1307.6346].
- Fermi-LAT Collaboration, M. Ackermann et al., Searching for Dark Matter Annihilation from Milky Way Dwarf Spheroidal Galaxies with Six Years of Fermi-LAT Data, arXiv:1503.02641.
- F. Feroz, M. Hobson, and M. Bridges, MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics, Mon.Not.Roy.Astron.Soc. 398 (2009) 1601–1614, [arXiv:0809.3437].
- M. Schumann, L. Baudis, L. Bütikofer, A. Kish, and M. Selvi, Dark matter sensitivity of multi-ton liquid xenon detectors, ArXiv e-prints (June, 2015) [arXiv:1506.08309].
- M. Drees and M. Nojiri, Neutralino - nucleon scattering revisited, Phys. Rev. D48 (1993) 3483–3501, [hep-ph/9307208].
- V. Mandic, A. Pierce, P. Gondolo, and H. Murayama, The Lower bound on the neutralino nucleon cross-section, hep-ph/0008022.
- J. Billard, L. Strigari, and E. Figueroa-Feliciano, Solar neutrino physics with low-threshold dark matter detectors, Phys. Rev. D91 (2015), no. 9 095023, [arXiv:1409.0050].
- M. Garny, A. Ibarra, M. Pato, and S. Vogl, On the spin-dependent sensitivity of XENON100, Phys.Rev. D87 (2013), no. 5 056002, [arXiv:1211.4573].
- IceCube Collaboration, P. Scott et al., Use of event-level neutrino telescope data in global fits for theories of new physics, JCAP 1211 (2012) 057, [arXiv:1207.0810].
- IceCube Collaboration, I. Collaboration, , to appear (2015).
- T. Bringmann and C. Weniger, Gamma Ray Signals from Dark Matter: Concepts, Status and Prospects, Phys. Dark Univ. 1 (2012) 194–217, [arXiv:1208.5481].
- ATLAS Collaboration, G. Aad et al., Search for top squark pair production in final states with one isolated lepton, jets, and missing transverse momentum in 8 TeV collisions with the ATLAS detector, JHEP 1411 (2014) 118, [arXiv:1407.0583].