Abstract: 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.

Abstract: The elastic scattering between dark matter particles and radiation represents an attractive possibility to solve a number of discrepancies between observations and standard cold dark matter predictions, as the induced collisional damping would imply a suppression of small-scale structures. We consider this scenario and confront it with measurements of the ionization history of the Universe at several redshifts and with recent estimates of the counts of Milky Way satellite galaxies. We derive a conservative upper bound on the dark matter photon elastic scattering cross section of sigma gamma DM < 8 x 10(-10) sigma(T) (m(DM)/GeV) at 95% CL, about one order of magnitude tighter than previous constraints from satellite number counts. Due to the strong degeneracies with astrophysical parameters, the bound on the dark matter-photon scattering cross section derived here is driven by the estimate of the number of Milky Way satellite galaxies. Finally, we also argue that future 21 cm probes could help in disentangling among possible non-cold dark matter candidates, such as interacting and warm dark matter scenarios. Let us emphasize that bounds of similar magnitude to the ones obtained here could be also derived for models with dark matter-neutrino interactions and would be as constraining as the tightest limits on such scenarios.

Abstract: We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce the basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of O(delta t(2)) accuracy like staggered leapfrog and Verlet integration, to Runge-Kutta methods up to O(delta t(4)) accuracy, and the Yoshida and Gauss-Legendre higher-order integrators, accurate up to O(delta t(10)) We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in 3+1 dimensions, including the case of 'self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: i) Interacting scalar fields, ii) Abelian U(1) gauge theories, and iii) Non-Abelian SU(2) gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from O(delta t(2)) up to O(delta t(10)) For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, including when 'self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under SU(2) x U(1). We note that the present manuscript is meant to be part of the theoretical basis for the code CosmoLattice, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, publicly available at http://www.cosrnolattice.net.

Abstract: Around one-third of the point-like sources in the Fermi-LAT catalogues remain as unidentified sources (unIDs) today. Indeed, these unIDs lack a clear, univocal association with a known astrophysical source. If dark matter (DM) is composed of weakly interacting massive particles (WIMPs), there is the exciting possibility that some of these unIDs may actually be DM sources, emitting gamma-rays from WIMPs annihilation. We propose a new approach to solve the standard, machine learning (ML) binary classification problem of disentangling prospective DM sources (simulated data) from astrophysical sources (observed data) among the unIDs of the 4FGL Fermi-LAT catalogue. We artificially build two systematic features for the DM data which are originally inherent to observed data: the detection significance and the uncertainty on the spectral curvature. We do it by sampling from the observed population of unIDs, assuming that the DM distributions would, if any, follow the latter. We consider different ML models: Logistic Regression, Neural Network (NN), Naive Bayes, and Gaussian Process, out of which the best, in terms of classification accuracy, is the NN, achieving around 93 . 3 per cent +/- 0 . 7 per cent performance. Other ML evaluation parameters, such as the True Ne gativ e and True Positive rates, are discussed in our work. Applying the NN to the unIDs sample, we find that the de generac y between some astrophysical and DM sources can be partially solved within this methodology. None the less, we conclude that there are no DM source candidates among the pool of 4FGL Fermi-LAT unIDs.

Abstract: Light primordial black holes (PBHs) with masses smaller than 10(9) g (10(-24) M-circle dot) evaporate before the onset of Big-Bang nucleosynthesis, rendering their detection rather challenging. If efficiently produced, they may have dominated the universe energy density. We study how such an early matter-dominated era can be probed successfully using gravitational waves (GW) emitted by local and global cosmic strings. While previous studies showed that a matter era generates a single-step suppression of the GW spectrum, we instead find a double-step suppression for local-string GW whose spectral shape provides information on the duration of the matter era. The presence of the two steps in the GW spectrum originates from GW being produced through two events separated in time: loop formation and loop decay, taking place either before or after the matter era. The second step – called the knee – is a novel feature which is universal to any early matter-dominated era and is not only specific to PBHs. Detecting GWs from cosmic strings with LISA, ET, or BBO would set constraints on PBHs with masses between 10(6) and 10(9) g for local strings with tension G μ= 10(-11), and PBHs masses between 10(4) and 10(9) g for global strings with symmetry-breaking scale eta = 10(15) GeV. Effects from the spin of PBHs are discussed.