Abstract: We study the impact of non-standard momentum distributions of cosmic neutrinos on the anisotropy spectrum of the cosmic microwave background and the matter power spectrum of the large scale structure. We show that the neutrino distribution has almost no unique observable imprint, as it is almost entirely degenerate with the effective number of neutrino flavours, N-eff, and the neutrino mass, m(nu). Performing a Markov chain Monte Carlo analysis with current cosmological data, we demonstrate that the neutrino mass bound heavily depends on the assumed momentum distribution of relic neutrinos. The message of this work is simple and has to our knowledge not been pointed out clearly before: cosmology allows that neutrinos have larger masses if their average momentum is larger than that of a perfectly thermal distribution. Here we provide an example in which the mass limits are relaxed by a factor of two.

Abstract: We present a full Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and Cosmic Microwave Background observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, namely Normal Ordering or Inverted Ordering, with current data. Two possible parametrizations (three neutrino masses versus the lightest neutrino mass plus the two oscillation mass splittings) and priors (linear versus logarithmic) are exhaustively examined. We find that the preference for NO is only driven by neutrino oscillation data. Moreover, the values of the Bayes factor indicate that the evidence for NO is strong only when the scan is performed over the three neutrino masses with logarithmic priors; for every other combination of parameterization and prior, the preference for NO is only weak. As a by-product of our Bayesian analyses, we are able to (a) compare the Bayesian bounds on the neutrino mixing parameters to those obtained by means of frequentist approaches, finding a very good agreement; (b) determine that the lightest neutrino mass plus the two mass splittings parametrization, motivated by the physical observables, is strongly preferred over the three neutrino mass eigenstates scan and (c) find that logarithmic priors guarantee a weakly-to-moderately more efficient sampling of the parameter space. These results establish the optimal strategy to successfully explore the neutrino parameter space, based on the use of the oscillation mass splittings and a logarithmic prior on the lightest neutrino mass, when combining neutrino oscillation data with cosmology and neutrinoless double beta decay. We also show that the limits on the total neutrino mass Sigma m(nu) can change dramatically when moving from one prior to the other. These results have profound implications for future studies on the neutrino mass ordering, as they crucially state the need for self-consistent analyses which explore the best parametrization and priors, without combining results that involve different assumptions.

Abstract: The intersection of the cosmic and neutrino frontiers is a rich field where much discovery space still remains. Neutrinos play a pivotal role in the hot big bang cosmology, influencing the dynamics of the universe over numerous decades in cosmological history. Recent studies have made tremendous progress in understanding some properties of cosmological neutrinos, primarily their energy density. Upcoming cosmological probes will measure the energy density of relativistic particles with higher precision, but could also start probing other properties of the neutrino spectra. When convolved with results from terrestrial experiments, cosmology can become even more acute at probing new physics related to neutrinos or even Beyond the Standard Model (BSM). Any discordance between laboratory and cosmological data sets may reveal new BSM physics and/or suggest alternative models of cosmology. We give examples of the intersection between terrestrial and cosmological probes in the neutrino sector, and briefly discuss the possibilities of what different laboratory experiments may see in conjunction with cosmological observatories.

Abstract: We employ state-of-the art cosmological observables including supernova surveys and BAO information to provide constraints on the mass and mixing angle of a non-resonantly produced sterile neutrino species, showing that cosmology can effectively rule out sterile neutrinos which decay between BBN and the present day. The decoupling of an additional heavy neutrino species can modify the time dependence of the Universe's expansion between BBN and recombination and, in extreme cases, lead to an additional matter-dominated period; while this could naively lead to a younger Universe with a larger Hubble parameter, it could later be compensated by the extra radiation expected in the form of neutrinos from sterile decay. However, recombination-era observables including the Cosmic Microwave Background (CMB), the shift parameter R-CMB and the sound horizon r(s) from Baryon Acoustic Oscillations (BAO) severely constrain this scenario. We self-consistently include the full time-evolution of the coupled sterile neutrino and standard model sectors in an MCMC, showing that if decay occurs after BBN, the sterile neutrino is essentially bounded by the constraint sin(2) theta less than or similar to 0.026(m(s)/eV)(-2).

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.