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 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: We show that the baryon asymmetry of the Universe can be generated by a strongly coupled right handed neutrino condensate which also drives inflation. The resulting model has only a small number of parameters, which completely determine not only the baryon asymmetry of the Universe and the mass of the right handed neutrino but also the inflationary phase. This feature allows us to make predictions that will be tested by current and planned experiments. As compared to the usual approach our dynamical framework is both economical and predictive.

Abstract: We discuss aidnogenesis,(1) i.e. the generation of a dark matter asymmetry, via new sphaleron processes associated to an extra non-abelian gauge symmetry common to both the visible and the dark sectors. Such a theory can naturally produce an abundance of asymmetric dark matter which is of the same size as the lepton and baryon asymmetries, as suggested by the similar sizes of the observed baryonic and dark matter energy content, and provide a definite prediction for the mass of the dark matter particle. We discuss in detail a minimal realization in which the Standard Model is only extended by dark matter fermions which form “dark baryons” through an SU(3) interaction, and a (broken) horizontal symmetry that induces the new sphalerons. The dark matter mass is predicted to be similar to 6GeV, close to the region favored by DAMA and CoGeNT. Furthermore, a remnant of the horizontal symmetry should be broken at a lower scale and can also explain the Tevatron dimuon anomaly.

Abstract: We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations which incorporate cold dark matter (CDM) and neutrinos as independent particle species. Our set of simulations explore the properties of neutrinos in a reference Lambda CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass 10(11) – 10(15) h(-1) M-circle dot, over a redshift range z = 0 – 2. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula, once the neutrino contribution to the total matter is removed. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than similar to 10(13.5) h(-1) M-circle dot. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. The results we obtain are numerically converged in terms of neutrino profiles at the 10% level for scales above similar to 200 h(-1) kpc at z = 0, and are stable with respect to box-size and starting redshift of the simulation. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.