Abstract: In this paper, we explore a scenario where the dark matter is a mixture of interacting and non interacting species. Assuming dark matter-photon interactions for the interacting species, we find that the suppression of the matter power spectrum in this scenario can mimic that expected in the case of massive neutrinos. Our numerical studies include present limits from Planck Cosmic Microwave Background data, which render the strength of the dark matter photon interaction unconstrained when the fraction of interacting dark matter is small. Despite the large entangling between mixed dark matter and neutrino masses, we show that future measurements from the Dark Energy Instrument (DESI) could help in establishing the dark matter and the neutrino properties simultaneously, provided that the interaction rate is very close to its current limits and the fraction of interacting dark matter is at least of O (10%). However, for that region of parameter space where a small fraction of interacting DM coincides with a comparatively large interaction rate, our analysis highlights a considerable degeneracy between the mixed dark matter parameters and the neutrino mass scale.

Abstract: Sterile neutrinos with a mass in the eV range have been invoked as a possible explanation of a variety of short baseline (SBL) neutrino oscillation anomalies. However, if one considers neutrino oscillations between active and sterile neutrinos, such neutrinos would have been fully thermalised in the early universe, and would be therefore in strong conflict with cosmological bounds. In this study we first update cosmological bounds on the mass and energy density of eV-scale sterile neutrinos. We then perform an updated study of a previously proposed model in which the sterile neutrino couples to a new light pseudoscalar degree of freedom. Consistently with previous analyses, we find that the model provides a good fit to all cosmological data and allows the high value of H-0 measured in the local universe to be consistent with measurements of the cosmic microwave background. However, new high l polarisation data constrain the sterile neutrino mass to be less than approximately 1 eV in this scenario. Finally, we combine the cosmological bounds on the pseudoscalar model with a Bayesian inference analysis of SBL data and conclude that only a sterile mass in narrow ranges around 1 eV remains consistent with both cosmology and SBL data.

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.