Abstract: Using the Fenchel-Eggleston theorem for convex hulls (an extension of the Caratheodory theorem), we prove that any likelihood can be maximized by either a dark matter 1-speed distribution F(v) in Earth's frame or 2-Galactic velocity distribution f(gal) ((u) over right arrow), consisting of a sum of delta functions. The former case applies only to time-averaged rate measurements and the maximum number of delta functions is (N-1), where N is the total number of data entries. The second case applies to any harmonic expansion coefficient of the time-dependent rate and the maximum number of terms is N. Using time-averaged rates, the aforementioned form of F(v) results in a piecewise constant unmodulated halo function (eta) over tilde (BF)-B-0 (v(min)) (which is an integral of the speed distribution) with at most (N-1) downward steps. The authors had previously proven this result for likelihoods comprised of at least one extended likelihood, and found the best-fit halo function to be unique. This uniqueness, however, cannot be guaranteed in the more general analysis applied to arbitrary likelihoods. Thus we introduce a method for determining whether there exists a unique best-fit halo function, and provide a procedure for constructing either a pointwise con fi dence band, if the best-fit halo function is unique, or a degeneracy band, if it is not. Using measurements of modulation amplitudes, the aforementioned form of f(gal) ((u) over right arrow), which is a sum of Galactic streams, yields a periodic time-dependent halo function (eta) over right arrow BF (v(min); t) which at any fixed time is a piecewise constant function of v(min) with at most N downward steps. In this case, we explain how to construct pointwise confidence and degeneracy bands from the time-averaged halo function. Finally, we show that requiring an isotropic Galactic velocity distribution leads to a Galactic speed distribution F(u)that is once again a sum of delta functions, and produces a time-dependent (eta) over tilde BF (v(min); t) function (and a time-averaged (eta) over tilde (0) BF (v(min))) that is piecewise linear, di ff ering significantly from best-fit halo functions obtained without the assumption of isotropy.

Abstract: The μnu SSM provides a solution to the mu-problem of the MSSM and explains the origin of neutrino masses by simply using right-handed neutrino superfields. Given that R-parity is broken in this model, the gravitino is a natural candidate for dark matter since its lifetime becomes much longer than the age of the Universe. We consider the implications of gravitino dark matter in the μnu SSM, analyzing in particular the prospects for detecting gamma rays from decaying gravitinos. If the gravitino explains the whole dark matter component, a gravitino mass larger than 20 GeV is disfavored by the isotropic diffuse photon background measurements. On the other hand, a gravitino with a mass range between 0.1 – 20 GeV gives rise to a signal that might be observed by the FERMI satellite. In this way important regions of the parameter space of the μnu SSM can be checked.

Abstract: The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search. In this work, we evaluate the proposed KDAR neutrino search strategies by realistically modeling both neutrino-nucleus interactions and the response of DUNE. We find that, although reconstruction of the neutrino energy and direction is difficult with current techniques in the relevant energy range, the superb energy resolution, angular resolution, and particle identification offered by DUNE can still permit great signal/background discrimination. Moreover, there are non-standard scenarios in which searches at DUNE for KDAR in the Sun can probe dark matter interactions.

Abstract: 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 fi rmly 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 is still 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 the 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 models with low fine-tuning 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.

Abstract: The energy spectra of particles produced from dark matter (DM) annihilation or decay are one of the fundamental ingredients to calculate the predicted fluxes of cosmic rays and radiation searched for in indirect DM detection. We revisit the calculation of the source spectra for annihilating and decaying DM using the VINCIA shower algorithm in PYTHIA to include QED and QCD final state radiation and diagrams for the EW corrections with massive bosons, not present in the default PYTHIA shower model. We take into account the spin information of the particles during the entire EW shower and the off -shell contributions from massive gauge bosons. Furthermore, we perform a dedicated tuning of the VINCIA and PYTHIA parameters to LEP data on the production of pions, photons, and hyperons at the Z resonance and discuss the underlying uncertainties. To enable the use of our results in DM studies, we provide the tabulated source spectra for the most relevant cosmic messenger particles, namely antiprotons, positrons, gamma rays and the three neutrino flavors, for all the fermionic and bosonic channels and DM masses between 5 GeV and 100 TeV, on github.