In this talk, we discussed the study of the QCD uncertainties on particle spectra from DM annihilation which we studied for the first time in [1]}. We demonstrated that the relative differences between the predictions of different multi-purposes MC event generators (Herwig 7.1.3, Pythia 8.235 and Sherpa 2.2.5) cannot be used to define a conservative estimate of QCD uncertainties particularly in the bulk of the spectra. We studied a complementary approach by using the same modeling paradigm (Pythia8) to define parametric variations taking the default Monash tune as our baseline and performed several retunings using data from LEP. Next, we show quantitatively the impact of the QCD uncertainties on the spectra of gamma-rays from DM annihilation in two benchmark points in the pMSSM. Full data tables which can be used to update those in the PPPC4DMID are public now on Zenodo and can be found in http://doi.org/10.5281/zenodo.3764809.
