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: At present, cosmological observations set the most stringent bound on the neutrino mass scale. Within the standard cosmological model (Lambda CDM), the Planck collaboration reports Sigma m(v)< 0.12 eV at 95 % CL. This bound, taken at face value, excludes many neutrino mass models. However, unstable neutrinos, with lifetimes shorter than the age of the universe <tau>(nu) less than or similar to t(U), represent a particle physics avenue to relax this constraint. Motivated by this fact, we present a taxonomy of neutrino decay modes, categorizing them in terms of particle content and final decay products. Taking into account the relevant phenomenological bounds, our analysis shows that 2-body decaying neutrinos into BSM particles are a promising option to relax cosmological neutrino mass bounds. We then build a simple extension of the type I seesaw scenario by adding one sterile state nu (4) and a Goldstone boson phi, in which nu (i)-> nu (4)phi decays can loosen the neutrino mass bounds up to Sigma m(v) similar to 1 eV, without spoiling the light neutrino mass generation mechanism. Remarkably, this is possible for a large range of the right-handed neutrino masses, from the electroweak up to the GUT scale. We successfully implement this idea in the context of minimal neutrino mass models based on a U(1)(mu-tau) flavor symmetry, which are otherwise in tension with the current bound on Sigma m(v).

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.