Abstract: 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(*) -> qq-', 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 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.

Abstract: The measured branching fractions of B-mesons into leptonic final states derived by the LHCb, Belle and BaBar collaborations hint towards the breakdown of lepton flavour universality. In this work we take at face value the so-called R-D(()*()) observables that are defined as the ratios of neutral B-meson charged-current decays into a D-(*())-meson, a charged lepton and a neutrino final state in the tau and light lepton channels. A well-studied and simple solution to this charged current anomaly is to introduce a scalar leptoquark S-1 that couples to the second and third generation of fermions. We investigate how S-1 can also serve as a mediator between the Standard Model and a dark sector. We study this scenario in detail and estimate the constraints arising from collider searches for leptoquarks, collider searches for missing energy signals, direct detection experiments and the dark matter relic abundance. We stress that the production of a pair of leptoquarks that decays into different final states (i.e. the commonly called “mixed” channels) provides critical information for identifying the underlying dynamics, and we exemplify this by studying the t tau b nu and the resonant S-1 plus missing energy channels. We find that direct detection data provides non-negligible constraints on the leptoquark coupling to the dark sector, which in turn affects the relic abundance. We also show that the correct relic abundance can not only arise via standard freeze-out, but also through conversion-driven freeze-out. We illustrate the rich phenomenology of the model with a few selected benchmark points, providing a broad stroke of the interesting connection between lepton flavour universality violation and dark matter.

Abstract: Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron de-cays, before particles that are stable on astrophysical time scales are produced. Antiprotons produced in this way may leave footprints in experiments such as AMS-02. Several groups have reported an excess of events in the antiproton flux in the rigidity range of 10-20 GV. However, the theoretical modeling of baryon production is not straightforward and relies in part on phenomenological models in Monte Carlo event generators. In this work, we assess the impact of QCD uncertainties on the spectra of antiprotons from dark-matter annihila-tion. As a proof-of-principle, we show that for a two-parameter model that depends only on the thermally-averaged annihilation cross section ((o -v)) and the dark-matter mass (Mx), QCD uncertainties can affect the best-fit mass by up to ti 14% (with large uncertainties for large DM masses), depending on the choice of Mx and the annihilation channel (bb over bar or W+W-), and (o -v) by up to ti 10%. For comparison, changes to the underlying diffusion parameters are found to be within 1%-5%, and the results are also quite resilient to the choice of cosmic-ray propagation model. These findings indicate that QCD uncertainties need to be included in future DM analyses. To facilitate full-fledged analyses, we provide the spectra in tabulated form including QCD uncertainties and code snippets to perform mass interpolations and quick DM fits. The code can be found in this GitHub [1] repository.