Abstract: The hypermagnetic helicity density at the electroweak phase transition (EWPT) exceeds many orders of magnitude the galactic magnetic helicity density. Together with previous magnetic helicity evolution calculations after the EWPT and hypermagnetic helicity conversion to the magnetic one at the EWPT, the present calculation completes the description of the evolution of this important topological feature of cosmological magnetic fields. It suggests that if the magnetic field seeding the galactic dynamo has a primordial origin, it should be substantially helical. This should be taken into account in scenarios of galactic magnetic field evolution with a cosmological seed.

Abstract: The observation of space-time variations in fundamental constants would provide strong evidence for the existence of new light degrees of freedom in the theory of Nature. Robustly constraining such scenarios requires exploiting observations that span different scales and probe the state of the Universe at different epochs. In the context of cosmology, both the cosmic microwave background and the Lyman-a forest have proven to be powerful tools capable of constraining variations in electromagnetism, however at the moment there do not exist cosmological probes capable of bridging the gap between recombination and reionization. In the near future, radio telescopes will attempt to measure the 21 cm transition of neutral hydrogen during the epochs of reionization and the cosmic dawn (and potentially the tail end of the dark ages); being inherently sensitive to electromagnetic phenomena, these experiments will offer a unique perspective on space-time variations of the fine-structure constant and the electron mass. We show here that large variations in these fundamental constants would produce features on the 21 cm power spectrum that may be distinguishable from astrophysical uncertainties. Furthermore, we forecast the sensitivity for the Square Kilometer Array, and show that the 21 cm power spectrum may be able to constrain variations at the level of O(10(-3)).

Abstract: Various particle physics models suggest that, besides the (nearly) cold dark matter that accounts for current observations, additional but sub-dominant dark relics might exist. These could be warm, hot, or even contribute as dark radiation. We present here a comprehensive study of two-component dark matter scenarios, where the first component is assumed to be cold, and the second is a non-cold thermal relic. Considering the cases where the non-cold dark matter species could be either a fermion or a boson, we derive consistent upper limits on the non-cold dark relic energy density for a very large range of velocity dispersions, covering the entire range from dark radiation to cold dark matter. To this end, we employ the latest Planck Cosmic Microwave Background data, the recent BOSS DR11 and other Baryon Acoustic Oscillation measurements, and also constraints on the number of Milky Way satellites, the latter of which provides a measure of the suppression of the matter power spectrum at the smallest scales due to the free-streaming of the non-cold dark matter component. We present the results on the fraction f(ncdm) of non-cold dark matter with respect to the total dark matter for different ranges of the non-cold dark matter masses. We find that the 2 sigma limits for non-cold dark matter particles with masses in the range 1-10 keV are f(ncdm) <= 0.29 (0.23) for fermions (bosons), and for masses in the 10-100 keV range they are f(ncdm) <= 0.43 (0.45), respectively.

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 present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, N-eff(SM), that quantifies the cosmological neutrinoto-photon energy densities. The calculation takes into account neutrino flavour oscillations, finite-temperature effects in the quantum electrodynamics plasma to O(e(3)), where e is the elementary electric charge, and a full evaluation of the neutrino-neutrino collision integral. We provide furthermore a detailed assessment of the uncertainties in the benchmark N(eff)(SM )value, through testing the value's dependence on (i) optional approximate modelling of the weak collision integrals, (ii) measurement errors in the physical parameters of the weak sector, and (iii) numerical convergence, particularly in relation to momentum discretisation. Our new, recommended standard-model benchmark is N-eff(SM) 3.0440 +/- 0.0002, where the nominal uncertainty is attributed predominantly to errors incurred in the numerical solution procedure (vertical bar delta N-eff vertical bar similar to 10(-4)), augmented by measurement errors in the solar mixing angle sin(2) theta(12) (vertical bar delta N-eff vertical bar similar to 10(-4)).