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.

The spectrum of the Cosmic Microwave Background follows Planck's black body radiation formula and shows a remarkable constant temperature of T = 2.7. About 380 000 years after the Big Bang at a temperature of T = 3000 Kelvin in the matter dominated era the electrons combine with the protons and 4He and the photons move freely in the neutral universe. So the temperature and distribution of the photons give us information of the universe 380 000 years after the Big Bang. Information about earlier times can, in principle, be derived from the Cosmic Neutrino Background (relic neutrinos). The neutrinos decouple already about 1 second after the Big Bang at a temperature of around 1 MeV or 10^{10} Kelvin. Today their temperature is about 1.95 Kelvin. Registration of these neutrinos is an extremely challenging experimental problem, which can hardly be solved with the present technologies. On the other hand it represents a tempting opportunity to check one of key elements of the Big Bang Cosmology and to probe the early stages of the universe evolution. The search for the Cosmic Neutrino Background with the induced beta decay: relic neutrino + 3H --> 3He + e-, is the topic of this contribution. The signal would show up as a peak in the electron spectrum by an energy with the neutrino mass above the Q value. We discuss the prospects of this approach and argue that it is able to set limits on the Cosmic Neutrino density in our vicinity. We also discuss critically ways to increase with modifications of the present KATRIN spectrometer the Tritium source intensity by a factor 100, which would yield about 170 counts of relic neutrino captures per year. Presently such an increase of the Tritium source intensity seems not to be possible. But one should be able to find an upper limit for the local density of the relic neutrinos in our galaxy.

Theories beyond the standard model include a number of new particles some of which might be light and weakly coupled to ordinary matter. Such particles affect equation of state of nuclear matter and can shift admissible masses of neutron stars to higher values. The internal structure of neutron stars is modified provided the ratio between coupling strength and mass squared of a weakly interacting light boson is above $g^2/\mu^2 \sim 25 ~\mathrm{GeV}^{-2}$. We provide limits on the couplings with the strange sector, which cannot be achieved from laboratory experiments analysis. When the couplings to the first family of quarks is considered the limits imposed by the neutron stars are not more stringent than the existing laboratory ones. The observations on neutron stars give evidence that equation of state of the $\beta$-equilibrated nuclear matter is stiffer than expected from many-body theory of nuclei and nuclear matter. A weakly interacting light vector boson coupled predominantly to the second family of the quarks can produce the required stiffening.

A generalized method to calculate the excitation-energy dependent parity ratio in the nuclear level density is presented, using the assumption of Poisson distributed independent quasi particles combined with BCS occupation numbers. It is found that it is crucial to employ a sufficiently large model space to allow excitations both from low-lying shells and to higher shells beyond a single major shell. Parity ratios are only found to equilibrate above at least 5-10 MeV of excitation energy. Furthermore, an overshooting effect close to major shells is found where the parity opposite to the ground state parity may dominate across a range of several MeV before the parity ratio finally equilibrates. The method is suited for large-scale calculations as needed, for example, in astrophysical applications. Parity distributions were computed for all nuclei from the proton dripline to the neutron dripline and from Ne up to Bi. These results were then used to recalculate astrophysical reaction rates in a Hauser-Feshbach statistical model. Although certain transitions can be considerably enhanced or suppressed, the impact on astrophysically relevant reactions remains limited, mainly due to the thermal population of target states in stellar reaction rates.

We present Dirac-Brueckner-Hartree-Fock calculations for isospin asymmetric nuclear matter which are based on improved approximations schemes. The potential matrix elements have been adapted for isospin asymmetric nuclear matter in order to account for the proton-neutron mass splitting in a more consistent way. The proton properties are particularly sensitive to this adaption and its consequences, whereas the neutron properties remains almost unaffected in neutron rich matter. Although at present full Brueckner calculations are still too complex to apply to finite nuclei, these relativistic Brueckner results can be used as a guidance to construct a density dependent relativistic mean field theory, which can be applied to finite nuclei. It is found that an accurate reproduction of the Dirac-Brueckner-Hartree-Fock equation of state requires a renormalization of these coupling functions.

A new scheme for testing nuclear matter equations of state (EsoS) at high densities using constraints from neutron star phenomenology and a flow data analysis of heavy-ion collisions is suggested. An acceptable EoS shall not allow the direct Urca process to occur in neutron stars with masses below $1.5~M_{\odot}$, and also shall not contradict flow and kaon production data of heavy-ion collisions. Compact star constraints include the mass measurements of 2.1 +/- 0.2 M_sun (1 sigma level) for PSR J0751+1807, of 2.0 +/- 0.1 M_sun from the innermost stable circular orbit for 4U 1636-536, the baryon mass - gravitational mass relationships from Pulsar B in J0737-3039 and the mass-radius relationships from quasiperiodic brightness oscillations in 4U 0614+09 and from the thermal emission of RX J1856-3754. This scheme is applied to a set of relativistic EsoS constrained otherwise from nuclear matter saturation properties with the result that no EoS can satisfy all constraints simultaneously, but those with density-dependent masses and coupling constants appear most promising.

We have examined the importance of momentum dependent induced nucleon currents such as weak-magnetism and pseudoscalar couplings to the amplitude of neutrinoless double beta decay in the mechanisms of light and heavy Majorana neutrino as well as in that of Majoron emission. Such effects are expected to occur in all nuclear models in the direction of reducing the light neutrino matrix elements by about 30%. To test this we have performed a calculation of the nuclear matrix elements of the experimentally interesting nuclei A = 76, 82, 96, 100, 116, 128, 130, 136 and 150 within the pn-RQRPA. We have found that indeed such corrections vary somewhat from nucleus to nucleus, but in all cases they are greater than 25 percent. In the case of heavy neutrino the effect is much larger (a factor of 3). Combining out results with the best presently available experimental limits on the half-life of the neutrinoless double beta decay we have extracted new limits on the effective neutrino mass (light and heavy) and the effective Majoron coupling constant.

We study the pionic decay of a possible dibaryon d' --> N N \pi in the microscopic quark shell model. The initial d' dibaryon wave function (J^P=0^-, T=0) consists of one 1\hbar\omega six-quark shell-model s^5p[51]_X configuration. The most important final six-quark configurations s^6[6]_X, s^4p^2[42]_X and (s^4p^2-s^52s)[6]_X are properly projected onto the NN channel. The final state NN interaction is investigated by means of two phase-equivalent - but off-shell different - potential models. We demonstrate that the decay width \Gamma_d' depends strongly on the short-range behavior of the NN wave function. In addition, the width \Gamma_d' is very sensitive to the mass and size of the d' dibaryon. For dibaryon masses slightly above the experimentally suggested value M_d'=2.065 GeV, we obtain a pionic decay width of \Gamma_d' = 0.18 - 0.32 MeV close to the experimental value \Gamma_d' = 0.5 MeV.

We calculate the mass and structure of a J^P=0^-, T=0 six-quark system using a colored diquark-tetraquark cluster wave function and a nonrelativistic quark model Hamiltonian. The calculated mass is some 350 MeV above the empirical value if the same confinement strength as in the nucleon is used. If the effective two-body confinement strength is weaker in a compound six-quark system than in a single baryon, as expected from a simple harmonic oscillator model, one obtains M_d' = 2092 MeV close to experiment.