Abstract: One of the most promising strategies to identify the nature of dark matter consists in the search for new particles at accelerators and with so-called direct detection experiments. Working within the framework of simplified models, and making use of machine learning tools to speed up statistical inference, we address the question of what we can learn about dark matter from a detection at the LHC and a forthcoming direct detection experiment. We show that with a combination of accelerator and direct detection data, it is possible to identify newly discovered particles as dark matter, by reconstructing their relic density assuming they are weakly interacting massive particles (WIMPs) thermally produced in the early Universe, and demonstrating that it is consistent with the measured dark matter abundance. An inconsistency between these two quantities would instead point either towards additional physics in the dark sector, or towards a non-standard cosmology, with a thermal history substantially different from that of the standard cosmological model.

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: Stringent constraints from direct detection experiments and the Large Hadron Collider motivate us to consider models in which the dark matter does not directly couple to the Standard Model, but that instead annihilates into hidden sector particles which ultimately decay through small couplings to the Standard Model. We calculate the gamma-ray emission generated within the context of several such hidden sector models, including those in which the hidden sector couples to the Standard Model through the vector portal (kinetic mixing with Standard Model hypercharge), through the Higgs portal (mixing with the Standard Model Higgs boson), or both. In each case, we identify broad regions of parameter space in which the observed spectrum and intensity of the Galactic Center gamma-ray excess can easily be accommodated, while providing an acceptable thermal relic abundance and remaining consistent with all current constraints. We also point out that cosmic-ray antiproton measurements could potentially discriminate some hidden sector models from more conventional dark matter scenarios.

Abstract: Nonzero neutrino masses are required by the existence of flavour oscillations, with values of the order of at least 50 meV. We consider the gravitational clustering of relic neutrinos within the Milky Way, and used the N – one-body simulation technique to compute their density enhancement factor in the neighbourhood of the Earth with respect to the average cosmic density. Compared to previous similar studies, we pushed the simulation down to smaller neutrino masses, and included an improved treatment of the baryonic and dark matter distributions in the Milky Way. Our results are important for future experiments aiming at detecting the cosmic neutrino background, such as the Princeton Tritium Observatory for Light, Early-universe, Massive-neutrino Yield (PTOLEMY) proposal. We calculate the impact of neutrino clustering in the Milky Way on the expected event rate for a PTOLEMY-like experiment. We find that the effect of clustering remains negligible for the minimal normal hierarchy scenario, while it enhances the event rate by 10 to 20% (resp. a factor 1.7 to 2.5) for the minimal inverted hierarchy scenario (resp. a degenerate scenario with 150 meV masses). Finally we compute the impact on the event rate of a possible fourth sterile neutrino with a mass of 1.3 eV.

Abstract: After the first four years of data taking, the IceCube neutrino telescope has observed 54 high-energy starting events (HESE) with deposited energies between 20TeV and 2PeV. The background from atmospheric muons and neutrinos is expected to be of about 20 events, all below 100TeV, thus pointing towards the astrophysical origin of about 8 events per year in that data set. However, their precise origin remains unknown. Here, we perform a detailed analysis of this event sample (considering simultaneously the energy, hemisphere and topology of the events) by assuming two contributions for the signal events: an isotropic power-law flux and a flux from decaying heavy dark matter. We fit the mass and lifetime of the dark matter and the normalization and spectral index of an isotropic power-law flux, for various decay channels of dark matter. We find that a significant contribution from dark matter decay is always slightly favored, either to explain the excess below 100TeV, as in the case of decays to quarks or, as in the case of neutrino channels, to explain the three multi-PeV events. Also, we consider the possibility to interpret all the data by dark matter decays only, considering various combinations of two decay channels. We show that the decaying dark matter scenario provides a better fit to HESE data than the isotropic power-law flux.