Abstract: We revisit a model in which neutrino masses and mixing are described by a two right-handed (RH) neutrino seesaw scenario, implying a strictly hierarchical light neutrino spectrum. A third decoupled RH neutrino, N-DM with mass M-DM, plays the role of cold dark matter (DM) and is produced by the mixing with a source RH neutrino, Ns with mass M-S, induced by Higgs portal interactions. The same interactions are also responsible for N-DM decays. We discuss in detail the constraints coming from DM abundance and stability conditions showing that in the hierarchical case, for M-DM >> M-S, there is an allowed window on M-DM values necessarily implying a contribution, from DM decays, to the high-energy neutrino flux recently detected by IceCube. We also show how the model can explain the matter-antimatter asymmetry of the Universe via leptogenesis in the quasi-degenerate limit. In this case, the DM mass should be within the range 300 GeV less than or similar to M-S < M-DM < 10PeV. We discuss the specific properties of this high-energy neutrino flux and show the predicted event spectrum for two exemplary cases. Although DM decays, with a relatively hard spectrum, cannot account for all the IceCube high-energy data, we illustrate how this extra source of high-energy neutrinos could reasonably explain some potential features in the observed spectrum. In this way, this represents a unified scenario for leptogenesis and DM that could be tested during the next years with more high-energy neutrino events.

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.