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: In the framework of a 3+1 scheme with an additional inert state, we consider the thermalisation of sterile neutrinos in the early Universe taking into account the full 4 x 4 mixing matrix. The evolution of the neutrino energy distributions is found solving the momentum-dependent kinetic equations with full diagonal collision terms, as in previous analyses of flavour neutrino decoupling in the standard case. The degree of thermalisation of the sterile state is shown in terms of the effective number of neutrinos, N-eff, and its dependence on the three additional mixing angles (theta(14), theta(24), theta(34)) and on the squared mass difference Delta m(41)(2) is discussed. Our results are relevant for fixing the contribution of a fourth light neutrino species to the cosmological energy density, whose value is very well constrained by the final Planck analysis. For the preferred region of active-sterile mixing parameters from short-baseline neutrino experiments, we find that the fourth state is fully thermalised (N-eff similar or equal to 4).

Abstract: In this paper we analyze a dark matter model inspired by theories with extra dimensions. The dark matter candidate corresponds to the first Kaluza-Klein mode of an real scalar added to the Standard Model. The tower of new particles enriches the calculation of the relic abundance. For large mass splitting, the model converges to the predictions of the inert singlet dark matter model. For nearly degenerate mass spectrum, coannihilations increase the cross-sections used for direct and indirect dark matter searches. Moreover, the Kaluza-Klein zero mode can mix with the SM higgs and further constraints can be applied.

Abstract: The impact of dark matter-neutrino interactions on the measurement of the cosmological parameters has been investigated in the past in the context of massless neutrinos exclusively. Here we revisit the role of a neutrino-dark matter coupling in light of ongoing cosmological tensions by implementing the full Boltzmann hierarchy for three massive neutrinos. Our tightest 95% CL upper limit on the strength of the interactions, parameterized via u(chi) = sigma(0)/sigma(Th) (m(chi)/100GeV)(-1), is u(chi) <= 3.34 . 10(-4), arising from a combination of Planck TTTEEE data, Planck lensing data and SDSS BAO data. This upper bound is, as expected, slightly higher than previous results for interacting massless neutrinos, due to the correction factor associated with neutrino masses. We find that these interactions significantly relax the lower bounds on the value of sigma 8 that is inferred in the context of Lambda CDM from the Planck data, leading to agreement within 1-2 sigma with weak lensing estimates of sigma 8, as those from KiDS1000. However, the presence of these interactions barely affects the value of the Hubble constant H-0.

Abstract: The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate with direct measurements that the non-neutrino backgrounds are below the expected cosmological signal. In this paper we discuss in details the theoretical aspects of the experiment and its physics goals. In particular, we mainly address three issues. First we discuss the sensitivity of PTOLEMY to the standard neutrino mass scale. We then study the perspectives of the experiment to detect the CNB via neutrino capture on tritium as a function of the neutrino mass scale and the energy resolution of the apparatus. Finally, we consider an extra sterile neutrino with mass in the eV range, coupled to the active states via oscillations, which has been advocated in view of neutrino oscillation anomalies. This extra state would contribute to the tritium decay spectrum, and its properties, mass and mixing angle, could be studied by analyzing the features in the beta decay electron spectrum.