Motivated by the recent galactic center gamma-ray excess identified in the Fermi-LAT data, we perform a detailed study of QCD fragmentation uncertainties in the modeling of the energy spectra of gamma-rays from Dark-Matter (DM) annihilation. When Dark-Matter particles annihilate to coloured final states, either directly or via decays such as $W^{(*)}\to q\bar{q}'$, photons are produced from a complex sequence of shower, hadronisation and hadron decays. In phenomenological studies, their energy spectra are typically computed using Monte Carlo event generators. These results have however intrinsic uncertainties due to the specific model used and the choice of model parameters, which are difficult to asses and which are typically neglected. We derive a new set of hadronisation parameters (tunes) for the \textsc{Pythia~8.2} Monte Carlo generator from a fit to LEP and SLD data at the $Z$ peak. For the first time, we also derive a conservative set of uncertainties on the shower and hadronisation model parameters. Their impact on the gamma-ray energy spectra is evaluated and discussed for a range of DM masses and annihilation channels. The spectra and their uncertainties are also provided in tabulated form for future use. The fragmentation-parameter uncertainties may be useful for collider studies as well.

Constraining the parameters of physical models with $>5-10$ parameters is a widespread problem in fields like particle physics and astronomy. The generation of data to explore this parameter space often requires large amounts of computational resources. The commonly used solution of reducing the number of relevant physical parameters hampers the generality of the results. In this paper we show that this problem can be alleviated by the use of active learning. We illustrate this with examples from high energy physics, a field where simulations are often expensive and parameter spaces are high-dimensional. We show that the active learning techniques query-by-committee and query-by-dropout-committee allow for the identification of model points in interesting regions of high-dimensional parameter spaces (e.g. around decision boundaries). This makes it possible to constrain model parameters more efficiently than is currently done with the most common sampling algorithms and to train better performing machine learning models on the same amount of data. Code implementing the experiments in this paper can be found on GitHub.

We present a new method to interpret the $\gamma$-ray data of our inner Galaxy as measured by the Fermi Large Area Telescope (Fermi LAT). We train and test convolutional neural networks with simulated Fermi-LAT images based on models tuned to real data. We use this method to investigate the origin of an excess emission of GeV $\gamma$-rays seen in previous studies. Interpretations of this excess include $\gamma$ rays created by the annihilation of dark matter particles and $\gamma$ rays originating from a collection of unresolved point sources, such as millisecond pulsars. Our new method allows precise measurements of the contribution and properties of an unresolved population of $\gamma$-ray point sources in the interstellar diffuse emission model.

The Fermi Collaboration has recently updated their analysis of gamma rays from the center of the Galaxy. They reconfirm the presence of an unexplained emission feature which is most prominent in the region of $1-10$ GeV, known as the Galactic Center GeV excess (GCE). Although the GCE is now firmly detected, an interpretation of this emission as a signal of self-annihilating dark matter (DM) particles is not unambiguously possible due to systematic effects in the gamma-ray modeling estimated in the Galactic Plane. In this paper we build a covariance matrix, collecting different systematic uncertainties investigated in the Fermi Collaboration's paper that affect the GCE spectrum. We show that models where part of the GCE is due to annihilating DM can still be consistent with the new data. We also re-evaluate the parameter space regions of the minimal supersymmetric Standard Model (MSSM) that can contribute dominantly to the GCE via neutralino DM annihilation. All recent constraints from DM direct detection experiments such as PICO, LUX, PandaX and Xenon1T, limits on the annihilation cross section from dwarf spheroidal galaxies and Large Hadron Collider limits are considered in this analysis. Due to a slight shift in the energy spectrum of the GC excess with respect to the previous Fermi analysis, and the recent limits from direct detection experiments, we find a slightly shifted parameter region of the MSSM compared to our previous analysis that is consistent with the GCE. Neutralinos with a mass between $85-220$ GeV can describe the excess via annihilation into a pair of $W$-bosons or top quarks. Remarkably, there are low fine-tuning models among the regions that we have found. The complete set of solutions will be probed by upcoming direct detection experiments and with dedicated searches in the upcoming data of the Large Hadron Collider.

Global fit studies performed in the pMSSM and the photon excess signal originating from the Galactic Center seem to suggest compressed electroweak supersymmetric spectra with a $\sim$100 GeV bino-like dark matter particle. We find that these scenarios are not probed by traditional electroweak supersymmetry searches at the LHC. We propose to extend the ATLAS and CMS electroweak supersymmetry searches with an improved strategy for bino-like dark matter, focusing on chargino plus next-to-lightest neutralino production, with a subsequent decay into a tri-lepton final state. We explore the sensitivity for pMSSM scenarios with $\Delta m = m_{\rm NLSP} - m_{\rm LSP} \sim (5 - 50)$ GeV in the $\sqrt{s} = 14$ TeV run of the LHC. Counterintuitively, we find that the requirement of low missing transverse energy increases the sensitivity compared to the current ATLAS and CMS searches. With 300 fb$^{-1}$ of data we expect the LHC experiments to be able to discover these supersymmetric spectra with mass gaps down to $\Delta m \sim 9$ GeV for DM masses between 40 and 140 GeV. We stress the importance of a dedicated search strategy that targets precisely these favored pMSSM spectra.

Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used to search for solutions within the pMSSM. We found solutions that are consistent with astroparticle physics and collider experiments, and provide a fit to the energy spectrum of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with a mass in the range $84-92$~GeV or $87-97$~GeV annihilating into W bosons. A third solutions is found for a neutralino of mass $174-187$~GeV annihilating into top quarks. The best solutions yield a Dark Matter relic density $0.06 < \Omega h^2 <0.13$. These pMSSM solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If the MSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be discovered soon.

We present statistically convergent profile likelihood maps obtained via global fits of a phenomenological Minimal Supersymmetric Standard Model with 15 free parameters (the MSSM-15), based on over 250M points. We derive constraints on the model parameters from direct detection limits on dark matter, the Planck relic density measurement and data from accelerator searches. We provide a detailed analysis of the rich phenomenology of this model, and determine the SUSY mass spectrum and dark matter properties that are preferred by current experimental constraints. We evaluate the impact of the measurement of the anomalous magnetic moment of the muon ($g-2$) on our results, and provide an analysis of scenarios in which the lightest neutralino is a subdominant component of the dark matter. The MSSM-15 parameters are relatively weakly constrained by current data sets, with the exception of the parameters related to dark matter phenomenology ($M_1$, $M_2$, $\mu$), which are restricted to the sub-TeV regime, mainly due to the relic density constraint. The mass of the lightest neutralino is found to be < 1.5 TeV at 99% C.L., but can extend up to 3 TeV when excluding the $g - 2$ constraint from the analysis. Low-mass bino-like neutralinos are strongly favoured, with spin-independent scattering cross-sections extending to very small values, $\sim 10^{-20}$ pb. ATLAS SUSY null searches strongly impact on this mass range, and thus rule out a region of parameter space that is outside the reach of any current or future direct detection experiment. The best-fit point obtained after inclusion of all data corresponds to a squark mass of 2.3 TeV, a gluino mass of 2.1 TeV and a 130 GeV neutralino with a spin-independent cross-section of $2.4 \times 10^{-10}$ pb, which is within the reach of future multi-ton scale direct detection experiments and of the upcoming LHC run at increased centre-of-mass energy.