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.

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.