The investigation of distant galaxy formation and evolution is a powerful tool to constrain dark matter scenarios, supporting and in some cases surpassing other astrophysical and experimental probes. The recent completion of the Hubble Frontier Fields (HFF) programme combining ultra-deep Hubble Space Telescope observations and the magnification power of gravitational lensing produced by foreground galaxy clusters has enabled the detection of the faintest primordial galaxies ever studied. Here we show how the number density of such primordial galaxies allows to constrain a variety of DM models alternative to CDM. In particular, it provides stringent limits on the mass of thermal WDM candidates, on the parameter space of sterile neutrino production models, and on other DM scenarios featuring particles in the keV mass range which is also supported by recent detections of a 3.5keV X-ray line. These constraints are robust and independent of the baryonic physics modeling of galaxy formation and evolution. Fuzzy DM (ultralight DM particles) results strongly disfavored.

Structure formation at small cosmological scales provides an important frontier for dark matter (DM) research. Scenarios with small DM particle masses, large momenta or hidden interactions tend to suppress the gravitational clustering at small scales. The details of this suppression depend on the DM particle nature, allowing for a direct link between DM models and astrophysical observations. However, most of the astrophysical constraints obtained so far refer to a very specific shape of the power suppression, corresponding to thermal warm dark matter (WDM), i.e., candidates with a Fermi-Dirac or Bose-Einstein momentum distribution. In this work we introduce a new analytical fitting formula for the power spectrum, which is simple yet flexible enough to reproduce the clustering signal of large classes of non-thermal DM models, which are not at all adequately described by the oversimplified notion of WDM. We show that the formula is able to fully cover the parameter space of sterile neutrinos (whether resonantly produced or from particle decay), mixed cold and warm models, fuzzy dark matter, as well as other models suggested by effective theory of structure formation (ETHOS). Based on this fitting formula, we perform a large suite of N-body simulations and we extract important nonlinear statistics, such as the matter power spectrum and the halo mass function. Finally, we present first preliminary astrophysical constraints, based on linear theory, from both the number of Milky Way satellites and the Lyman-$\alpha$ forest. This paper is a first step towards a general and comprehensive modeling of small-scale departures from the standard cold DM model.

We perform a global Bayesian analysis of currently available neutrino data, putting data from oscillation experiments, neutrinoless double beta decay ($0\nu\beta\beta$), and precision cosmology on an equal footing. We evaluate the discovery potential of future $0\nu\beta\beta$ experiments and the Bayes factor of the two possible neutrino mass ordering schemes for different prior choices. We show that the indication for normal ordering is still very mild and does not strongly depend on realistic prior assumptions or different combinations of cosmological data sets. We find a wide range for $0\nu\beta\beta$ discovery potential, depending on the absolute neutrino mass scale, mass ordering and achievable background level.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We show that the measured abundance of ultra-faint lensed galaxies at $z\approx 6$ in the Hubble Frontier Fields (HFF) provides stringent constraints on the parameter space of i) Dark Matter models based on keV sterile neutrinos; ii) the "fuzzy" wavelike Dark Matter models, based on Bose-Einstein condensate of ultra-light particles. For the case of the sterile neutrinos, we consider two production mechanisms: resonant production through the mixing with active neutrinos and the decay of scalar particles. For the former model, we derive constraints for the combination of sterile neutrino mass $m_{\nu}$ and mixing parameter $\sin^2(2\theta)$ which provide the tightest lower bounds on the mixing angle (and hence on the lepton asymmetry) derived so far by methods independent of baryonic physics. For the latter we compute the allowed combinations of the scalar mass, its coupling to the Higgs field, and the Yukawa coupling of the scalar to the sterile neutrinos. We compare our results to independent, existing astrophysical bounds on sterile neutrinos in the same mass range. For the case of "fuzzy" Dark Matter, we show that the observed number density $\approx 1/{\rm Mpc}^3$ of high-redshift galaxies in the HFF sets a lower limit $m_\psi\geq 8\cdot 10^{-22}$ eV (at 3-$\sigma$ confidence level) on the particle mass, a result that strongly disfavors wavelike bosonic Dark Matter as a viable model for structure formation. We discuss the impact on our results of uncertainties due to systematics in the selection of highly magnified, faint galaxies at high redshifts.

We investigate the early Universe production of sterile neutrino Dark Matter by the decays of singlet scalars. All previous studies applied simplifying assumptions and/or studied the process only on the level of number densities, which makes it impossible to give statements about cosmic structure formation. We overcome these issues by dropping all simplifying assumptions (except for one we showed earlier to work perfectly) and by computing the full course of Dark Matter production on the level of non-thermal momentum distribution functions. We are thus in the position to study all aspects of the resulting settings and apply all relevant bounds in a reliable manner. We have a particular focus on how to incorporate bounds from structure formation on the level of the linear power spectrum, since the simplistic estimate using the free-streaming horizon clearly fails for highly non-thermal distributions. Our work comprises the most detailed and comprehensive study of sterile neutrino Dark Matter production by scalar decays presented so far.

The simplest way to create sterile neutrinos in the early Universe is by their admixture to active neutrinos. However, this mechanism, connected to the Dark Matter (DM) problem by Dodelson and Widrow (DW), cannot simultaneously meet the relic abundance constraint as well as bounds from structure formation and X-rays. Nonetheless, unless a symmetry forces active-sterile mixing to vanish exactly, the DW mechanism will unavoidably affect the sterile neutrino DM population created by any other production mechanism. We present a semi-analytic approach to the DW mechanism acting on an arbitrary initial abundance of sterile neutrinos, allowing to combine DW with any other preceding production mechanism in a physical and precise way. While previous analyses usually assumed that the spectra produced by DW and another mechanism can simply be added, we use our semi-analytic results to discuss the validity of this assumption and to quantify its accurateness, thereby also scrutinising the DW spectrum and the derived mass bounds. We then map our results to the case of sterile neutrino DM from the decay of a real SM singlet coupled to the Higgs. Finally, we will investigate aspects of structure formation beyond the usual simple free-streaming estimates in order to judge on the effects of the DW modification on the sterile neutrino DM spectra generated by scalar decay.

We perform a detailed and illustrative study of the production of keV sterile neutrino Dark Matter (DM) by decays of singlet scalars in the early Universe. In the current study we focus on providing a clear and general overview of this production mechanism. For the first time we study all regimes possible on the level of momentum distribution functions, which we obtain by solving a system of Boltzmann equations. These quantities contain the full information about the production process, which allows us to not only track the evolution of the DM generation but to also take into account all bounds related to the spectrum, such as constraints from structure formation or from avoiding too much dark radiation. In particular we show that this simple production mechanism can, depending on the regime, lead to strongly non-thermal DM spectra which may even feature more than one peak in the momentum distribution. These cases could have particularly interesting consequences for cosmological structure formation, as their analysis requires more refined tools than the simplistic estimate using the free-streaming horizon. Here we present the mechanism including all concepts and subtleties involved, for now using the assumption that the effective number of relativistic degrees of freedom is constant during DM production, which is applicable in a significant fraction of the parameter space. This allows us to derive analytical results to back up our detailed numerical computations, thus leading to the most comprehensive picture of keV sterile neutrino DM production by singlet scalar decays that exists up to now.

The search for sharp features in the gamma-ray spectrum is a promising approach to identify a signal from dark matter annihilations over the astrophysical backgrounds. In this paper we investigate the generation of gamma-ray lines and internal bremsstrahlung signals in a toy model where the dark matter particle is a real scalar that couples to a lepton and an exotic fermion via a Yukawa coupling. We show that the Fermi-LAT and H.E.S.S. searches for line-like spectral features severely constrain regions of the parameter space where the scalar dark matter is thermally produced. Finally, we also discuss the complementarity of the searches for sharp spectral features with other indirect dark matter searches, as well as with direct and collider searches.

Dark matter particles captured in the Sun would annihilate producing a neutrino flux that could be detected at the Earth. In some channels, however, the neutrino flux lies in the MeV range and is thus undetectable at IceCube, namely when the dark matter particles annihilate into $e^+e^-$, $\mu^+\mu^-$ or light quarks. On the other hand, the same interaction that mediates the annihilations into light fermions also leads, via higher order effects, to the production of weak gauge bosons (and in the case of quarks also gluons) that generate a high energy neutrino flux potentially observable at IceCube. We consider in this paper tree level annihilations into a fermion-antifermion pair with the associated emission of one gauge boson and one loop annihilations into two gauge bosons, and we calculate the limits on the scattering cross section of dark matter particles with protons in scenarios where the dark matter particle couples to electrons, muons or light quarks from the non-observation of an excess of neutrino events in the direction of the Sun. We find that the limits on the spin-dependent scattering cross section are, for some scenarios, stronger than the limits from direct detection experiments.

We investigate the prospects to observe a high energy neutrino signal from dark matter annihilations in the Sun in scenarios where the dark matter is a Majorana fermion that couples to a quark and a colored scalar via a Yukawa coupling. In this minimal scenario, the dark matter capture and annihilation in the Sun can be studied in a single framework. We find that, for small and moderate mass splitting between the dark matter and the colored scalar, the two-to-three annihilation $q \bar q g$ plays a central role in the calculation of the number of captured dark matter particles. On the other hand, the two-to-three annihilation into $q \bar q Z$ gives, despite its small branching fraction, the largest contribution to the neutrino flux at the Earth at the highest energies. We calculate the limits on the model parameters using IceCube observations of the Sun and we discuss their interplay with the requirement of equilibrium of captures and annihilations in the Sun and with the requirement of thermal dark matter production. We also compare the limits from IceCube to the limits from direct detection, antiproton measurements and collider searches.