Sterile neutrinos with a mass in the eV range have been invoked as a possible explanation of a variety of short baseline (SBL) neutrino oscillation anomalies. However, if one considers neutrino oscillations between active and sterile neutrinos, such neutrinos would have been fully thermalised in the early universe, and would be therefore in strong conflict with cosmological bounds. In this study we first update cosmological bounds on the mass and energy density of eV-scale sterile neutrinos. We then perform an updated study of a previously proposed model in which the sterile neutrino couples to a new light pseudoscalar degree of freedom. Consistently with previous analyses, we find that the model provides a good fit to all cosmological data and allows the high value of $H_0$ measured in the local universe to be consistent with measurements of the cosmic microwave background. However, new high $\ell$ polarisation data constrain the sterile neutrino mass to be less than approximately 1~eV in this scenario. Finally, we combine the cosmological bounds on the pseudoscalar model with a Bayesian inference analysis of SBL data and conclude that only a sterile mass in a narrow range around 1~eV remains consistent with both cosmology and SBL data.

We investigate how sterile neutrinos with a range of masses influence cosmology in MeV-scale reheating temperature scenarios. By computing the production of sterile neutrinos through the combination of mixing and scattering in the early Universe, we find that light sterile neutrinos, with masses and mixings as inferred from short-baseline neutrino oscillation experiments, are consistent with big-bang nucleosynthesis (BBN) and cosmic microwave background (CMB) radiation for the reheating temperature of ${\cal O}(1)$ MeV if the parent particle responsible for reheating decays into electromagnetic components (radiative decay). In contrast, if the parent particle mainly decays into hadrons (hadronic decay), the bound from BBN becomes more stringent. In this case, the existence of the light sterile neutrinos can be cosmologically excluded, depending on the mass and the hadronic branching ratio of the parent particle.

From a theoretical point of view, there is a strong motivation to consider an MeV-scale reheating temperature induced by long-lived massive particles with masses around the weak scale, decaying only through gravitational interaction. In this study, we investigate lower limits on the reheating temperature imposed by big-bang nucleosynthesis assuming both radiative and hadronic decays of such massive particles. For the first time, effects of neutrino self-interactions and oscillations are taken into account in the neutrino thermalization calculations. By requiring consistency between theoretical and observational values of light element abundances, we find that the reheating temperature should conservatively be $T_{\rm RH} \gtrsim 1.8$ MeV in the case of the 100% radiative decay, and $T_{\rm RH} \gtrsim$ 4-5 MeV in the case of the 100% hadronic decays for particle masses in the range of 10 GeV to 100 TeV.

We present $N$-body simulations in which either all, or a fraction of, the cold dark matter decays non-relativistically to a relativistic, non-interacting dark radiation component. All effects from radiation and general relativity are self-consistently included at the level of linear perturbation theory, and our simulation results therefore match those from linear Einstein-Boltzmann codes such as CLASS in the appropriate large-scale limit. We also find that standard, Newtonian $N$-body simulations adequately describe the non-linear evolution at smaller scales ($k \gtrsim 0.1 \, h/{\mathrm{Mpc}}$) in this type of model, provided that the mass of the decaying component is modified correctly, and that the background evolution is correctly treated. That is, for studies of small scales, effects from general relativity and radiation can be safely neglected.

We present $N$-body simulations which are fully compatible with general relativity, with dark energy consistently included at both the background and perturbation level. We test our approach for dark energy parameterised as both a fluid, and using the parameterised post-Friedmann (PPF) formalism. In most cases, dark energy is very smooth relative to dark matter so that its leading effect on structure formation is the change to the background expansion rate. This can be easily incorporated into Newtonian $N$-body simulations by changing the Friedmann equation. However, dark energy perturbations and relativistic corrections can lead to differences relative to Newtonian $N$-body simulations at the tens of percent level for scales $k < (10^{-3} \unicode{x2013} 10^{-2})\,\mathrm{Mpc}^{-1}$, and given the accuracy of upcoming large scale structure surveys such effects must be included. In this paper we will study both effects in detail and highlight the conditions under which they are important. We also show that our $N$-body simulations exactly reproduce the results of the Boltzmann solver CLASS for all scales which remain linear.

Bayesian parameter inference depends on a choice of prior probability distribution for the parameters in question. The prior which makes the posterior distribution maximally sensitive to data is called the Jeffreys prior, and it is completely determined by the response of the likelihood to changes in parameters. Under the assumption that the likelihood is a Gaussian distribution, the Jeffreys prior is a constant, i.e. flat. However, if one parameter is constrained by physical considerations, the Gaussian approximation fails and the flat prior is no longer the Jeffreys prior. In this paper we compute the correct Jeffreys prior for a multivariate normal distribution constrained in one dimension, and we apply it to the sum of neutrino masses $\Sigma m_\nu$ and the tensor-to-scalar ratio $r$. We find that one-dimensional marginalised posteriors for these two parameters change considerably and that the 68% and 95% Bayesian upper limits increase by 9% and 4% respectively for $\Sigma m_\nu$ and 22% and 3% for $r$. Adding the prior to an existing chain can be done as a trivial importance sampling in the final step of the analysis proces.

We show how Newtonian cosmological simulations can be employed to investigate the non-linear evolution of two particle species in a relativistic context. We discuss the application for massive neutrinos and other multi-species systems such as Cold Dark Matter (CDM) plus baryons or Warm Dark Matter (WDM). We propose a method that allows us to perform simulations including massive neutrinos and general relativistic effects at almost the same computational cost as ordinary CDM only N-body simulations, employing tailor-made initial conditions and a dictionary for the interpretation of the simulation output.

Cosmological $N$-body simulations are typically purely run with particles using Newtonian equations of motion. However, such simulations can be made fully consistent with general relativity using a well-defined prescription. Here, we extend the formalism previously developed for $\Lambda$CDM cosmologies with massless neutrinos to include the effects of massive, but light neutrinos. We have implemented the method in two different $N$-body codes, CONCEPT and PKDGRAV, and demonstrate that they produce consistent results. We furthermore show that we can recover all appropriate limits, including the full GR solution in linear perturbation theory at the per mille level of precision.

We introduce the Generalised Lensing and Shear Spectra GLaSS code which is available for download from https://github.com/astro-informatics/GLaSS It is a fast and flexible public code, written in Python, that computes generalized spherical cosmic shear spectra. The commonly used tomographic and spherical Bessel lensing spectra come as built-in run-mode options. GLaSS is integrated into the Cosmosis modular cosmological pipeline package. We outline several computational choices that accelerate the computation of cosmic shear power spectra. Using GLaSS, we test whether the assumption that using the lensing and projection kernels for a spatially-flat universe -- in a universe with a small amount of spatial curvature -- negligibly impacts the lensing spectrum. We refer to this assumption as The Spatially-Flat Universe Approximation, that has been implicitly assumed in all cosmic shear studies to date. We confirm that The Spatially-Flat Universe Approximation has a negligible impact on Stage IV cosmic shear experiments.

Including massive neutrinos in $N$-body simulations is a challenging task due to the large thermal velocities of the neutrinos. In particle based codes this leads to problems of shot-noise due to insufficient sampling of the neutrino momentum distribution function. In this paper we investigate the benefits and drawbacks of a scheme first suggested in a paper by Banerjee et al. in which the initial neutrino distribution is symmetrised in momentum space. We confirm that this method reduce shot-noise significantly, but we also find that it generates some spurious power in the neutrino power spectrum at intermediate and small scales. We speculate that this happens because many neutrinos in the simulation sample the same underlying dark matter structures while moving through the simulation. By carefully tuning the number of directions in momentum space of the initial neutrino distribution we show that some improvements can be made over the case where initial neutrino directions are purely random. At redshifts $z\gtrsim 3$ the method works very well, but at smaller redshifts significant improvements are not possible due to the spurious power generation.

In this paper the non-linear effect of massive neutrinos on cosmological structures is studied in a conceptually new way. We have solved the non-linear continuity and Euler equations for the neutrinos on a grid in real space in $N$-body simulations, and closed the Boltzmann hierarchy at the non-linear Euler equation using the stress and pressure perturbations from linear theory. By comparing with state-of-the art cosmological neutrino simulations, we are able to simulate the non-linear neutrino power spectrum very accurately. This translates into a negligible error in the matter power spectrum, and so our CONCEPT code is ideally suited for extracting the neutrino mass from future high precision non-linear observational probes such as Euclid.

We show how standard Newtonian N-body simulations can be interpreted in terms of the weak-field limit of general relativity by employing the recently developed Newtonian motion gauge. Our framework allows the inclusion of radiation perturbations and the non-linear evolution of matter. We show how to construct the weak-field metric by combining Newtonian simulations with results from Einstein-Boltzmann codes. We discuss observational effects on weak lensing and ray tracing, identifying important relativistic corrections.

We consider the impact of neutrino self-interactions described by an effective four-fermion coupling on cosmological observations. Implementing the exact Boltzmann hierarchy for interacting neutrinos first derived in [arxiv:1409.1577] into the Boltzmann solver CLASS, we perform a detailed numerical analysis of the effects of the interaction on the cosmic microwave background (CMB) anisotropies, and compare our results with known approximations in the literature. While we find good agreement between our exact approach and the relaxation time approximation used in some recent studies, the popular $\left( c_{\text{eff}}^2,c_{\text{vis}}^2 \right)$-parameterisation fails to reproduce the correct scale dependence of the CMB temperature power spectrum. We then proceed to derive constraints on the effective coupling constant $G_{\text{eff}}$ using currently available cosmological data via an MCMC analysis. Interestingly, our results reveal a bimodal posterior distribution, where one mode represents the standard $\Lambda$CDM limit with $G_{\rm eff} \lesssim 10^8 \, G_{\rm F}$, and the other a scenario in which neutrinos self-interact with an effective coupling constant $G_{\rm eff} \simeq 3 \times 10^9 \, G_{\rm F}$.

Newtonian N-body simulations have been employed successfully over the past decades for the simulation of the cosmological large-scale structure. Such simulations usually ignore radiation perturbations (photons and massless neutrinos) and the impact of general relativity (GR) beyond the background expansion. This approximation can be relaxed and we discuss three different approaches that are accurate to leading order in GR. For simulations that start at redshift less than about 100 we find that the presence of early radiation typically leads to percent-level effects on the numerical power spectra at large scales. Our numerical results agree across the three methods, and we conclude that all of the three methods are suitable for simulations in a standard cosmology. Two of the methods modify the N-body evolution directly, while the third method can be applied as a post-processing prescription.

In this paper we present a fast and accurate numerical algorithm for the computation of hyperspherical Bessel functions of large order and real arguments. For the hyperspherical Bessel functions of closed type, no stable algorithm existed so far due to the lack of a backwards recurrence. We solved this problem by establishing a relation to Gegenbauer polynomials. All our algorithms are written in C and are publicly available at Github [https://github.com/lesgourg/class_public]. A Python wrapper is available upon request.

Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This back-scaling method can account for the effect of inhomogeneous residual thermal radiation at early times, which is absent in the Newtonian simulations. We analyse this procedure from a fully relativistic perspective, employing the recently-proposed Newtonian motion gauge framework. We find that N-body simulations for LambdaCDM cosmology starting from back-scaled initial conditions can be self-consistently embedded in a relativistic space-time with first-order metric potentials calculated using a linear Boltzmann code. This space-time coincides with a simple "N-body gauge" for z<50 for all observable modes. Care must be taken, however, when simulating non-standard cosmologies. As an example, we analyse the back-scaling method in a cosmology with decaying dark matter, and show that metric perturbations become large at early times in the back-scaling approach, indicating a breakdown of the perturbative description. We suggest a suitable "forwards approach" for such cases.

We investigate how inflation model selection is affected by the presence of additional free-streaming relativistic degrees of freedom, i.e. dark radiation. We perform a full Bayesian analysis of both inflation parameters and cosmological parameters taking reheating into account self-consistently. We compute the Bayesian evidence for a few representative inflation scenarios in both the standard $\Lambda$CDM model and an extension including dark radiation parametrised by its effective number of relativistic species $N_\mathrm{eff}$. Using a minimal dataset (Planck low-$\ell$ polarisation, temperature power spectrum and lensing reconstruction), we find that the observational status of most inflationary models is unchanged. The exceptions are potentials such as power-law inflation that predict large values for the scalar spectral index that can only be realised when $N_\mathrm{eff}$ is allowed to vary. Adding baryon acoustic oscillations data and the B-mode data from BICEP2/Keck makes power-law inflation disfavoured, while adding local measurements of the Hubble constant $H_0$ makes power-law inflation slightly favoured compared to the best single-field plateau potentials. This illustrates how the dark radiation solution to the $H_0$ tension would have deep consequences for inflation model selection.

Cosmological $N$-body simulations are the standard tool to study the emergence of the observed large-scale structure of the Universe. Such simulations usually solve for the gravitational dynamics of matter within the Newtonian approximation, thus discarding general relativistic effects such as the coupling between matter and radiation ($\equiv$ photons and neutrinos). In this paper we investigate novel hybrid simulations which incorporate interactions between radiation and matter to the leading order in General Relativity, whilst evolving the matter dynamics in full non-linearity according to Newtonian theory. Our hybrid simulations come with a relativistic space-time and make it possible to investigate structure formation in a unified framework. In the present work we focus on simulations initialized at $z=99$, and show that the extracted matter power spectrum receives up to $3\%$ corrections on very large scales through radiation. Our numerical findings compare favourably with linear analytical results from Fidler et al. (2016), from which we deduce that there cannot be any significant non-linear mode-coupling induced through linear radiation corrections.

The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1$\%$ corrections from radiation on large scales if Newtonian simulations are initialised before $z=50$. We provide a relativistic framework in which \emph{unmodified} Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in $N$-body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein--Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the \emph{Newtonian motion} gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

The Short BaseLine (SBL) neutrino oscillation anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such a neutrino is incompatible with cosmological data, in particular observations of the Cosmic Microwave Background (CMB) anisotropies. However, this conclusion can change by invoking new physics. One possibility is to introduce a secret interaction in the sterile neutrino sector mediated by a light pseudoscalar. In this pseudoscalar model, CMB data prefer a sterile neutrino mass that is fully compatible with the mass ranges suggested by SBL anomalies. In addition, this model predicts a value of the Hubble parameter which is completely consistent with local measurements.

We present a fully relativistic calculation of the matter bispectrum at second order in cosmological perturbation theory assuming a Gaussian primordial curvature perturbation. For the first time we perform a full numerical integration of the bispectrum for both baryons and cold dark matter using the second-order Einstein-Boltzmann code, SONG. We review previous analytical results and provide an improved analytic approximation for the second-order kernel in Poisson gauge which incorporates Newtonian nonlinear evolution, relativistic initial conditions, the effect of radiation at early times and the cosmological constant at late times. Our improved kernel provides a percent level fit to the full numerical result at late times for most configurations, including both equilateral shapes and the squeezed limit. We show that baryon acoustic oscillations leave an imprint in the matter bispectrum, making a significant impact on squeezed shapes.

We study a coupled quintessence model with pure momentum exchange and present the effects of such an interaction on the Cosmic Microwave Background (CMB) and matter power spectrum. For a wide range of negative values of the coupling parameter $\beta$ structure growth is suppressed and the model can reconcile the tension between Cosmic Microwave Background observations and structure growth inferred from cluster counts. We find that this model is as good as $\Lambda$CDM for CMB and baryon acoustic oscillation (BAO) data, while the addition of cluster data makes the model strongly preferred, improving the best-fit $\chi^2$-value by more than $16$.

The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an Einstein-Boltzmann solver. We show that in most gauges the initial displacements generated in this way receive a first-order relativistic correction. We define a new gauge, the $N$-body gauge, in which this relativistic correction vanishes and show that a conventional Newtonian $N$-body simulation includes all first-order relativistic contributions (in the absence of radiation) if we identify the coordinates in Newtonian simulations with those in the relativistic $N$-body gauge.

We present a new flexible, fast and accurate way to implement massive neutrinos, warm dark matter and any other non-cold dark matter relics in Boltzmann codes. For whatever analytical or numerical form of the phase-space distribution function, the optimal sampling in momentum space compatible with a given level of accuracy is automatically found by comparing quadrature methods. The perturbation integration is made even faster by switching to an approximate viscous fluid description inside the Hubble radius, which differs from previous approximations discussed in the literature. When adding one massive neutrino to the minimal cosmological model, CLASS becomes just 1.5 times slower, instead of about 5 times in other codes (for fixed accuracy requirements). We illustrate the flexibility of our approach by considering a few examples of standard or non-standard neutrinos, as well as warm dark matter models.

Sterile neutrinos in the electronvolt mass range are hinted at by a number of terrestrial neutrino experiments. However, such neutrinos are highly incompatible with data from the Cosmic Microwave Background and large scale structure. This paper discusses how charging sterile neutrinos under a new pseudoscalar interaction can reconcile eV sterile neutrinos with terrestrial neutrino data. We show that this model can reconcile eV sterile neutrinos in cosmology, providing a fit to all available data which is way better than the standard $\Lambda$CDM model with one additional fully thermalized sterile neutrino. In particular it also prefers a value of the Hubble parameter much closer to the locally measured value.

Sterile neutrinos are thermalised in the early Universe via oscillations with the active neutrinos for certain mixing parameters. The most detailed calculation of this thermalisation process involves the solution of the momentum-dependent quantum kinetic equations, which track the evolution of the neutrino phase space distributions. Until now the collision terms in the quantum kinetic equations have always been approximated using equilibrium distributions, but this approximation has never been checked numerically. In this work we revisit the sterile neutrino thermalisation calculation using the full collision term, and compare the results with various existing approximations in the literature. We find a better agreement than would naively be expected, but also identify some issues with these approximations that have not been appreciated previously. These include an unphysical production of neutrinos via scattering and the importance of redistributing momentum through scattering, as well as details of Pauli blocking. Finally, we devise a new approximation scheme, which improves upon some of the shortcomings of previous schemes.

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. Such sterile neutrinos are incompatible with cosmology because they suppress structure formation unless they can be prevented from thermalising in the early Universe or removed by subsequent decay or annihilation. Here we present a novel scenario in which both sterile neutrinos and dark matter are coupled to a new, light pseudoscalar. This can prevent thermalisation of sterile neutrinos and make dark matter sufficiently self-interacting to have an impact on galactic dynamics and possibly resolve some of the known problems with the standard cold dark matter scenario. Even more importantly it leads to a strongly self-interacting plasma of sterile neutrinos and pseudoscalars at late times and provides an excellent fit to CMB data. The usual cosmological neutrino mass problem is avoided by sterile neutrino annihilation to pseudoscalars. The preferred value of $H_0$ is substantially higher than in standard $\Lambda$CDM and in much better agreement with local measurements.

The existence of a cosmic neutrino background can be probed indirectly by CMB experiments, not only by measuring the background density of radiation in the universe, but also by searching for the typical signatures of the fluctuations of free-streaming species in the temperature and polarisation power spectrum. Previous studies have already proposed a rather generic parametrisation of these fluctuations, that could help to discriminate between the signature of ordinary free-streaming neutrinos, or of more exotic dark radiation models. Current data are compatible with standard values of these parameters, which seems to bring further evidence for the existence of a cosmic neutrino background. In this work, we investigate the robustness of this conclusion under various assumptions. We generalise the definition of an effective sound speed and viscosity speed to the case of massive neutrinos or other dark radiation components experiencing a non-relativistic transition. We show that current bounds on these effective parameters do not vary significantly when considering an arbitrary value of the particle mass, or extended cosmological models with a free effective neutrino number, dynamical dark energy or a running of the primordial spectrum tilt. We conclude that it is possible to make a robust statement about the detection of the cosmic neutrino background by CMB experiments.

We present a new method for calculating CMB anisotropies in a non-flat Friedmann universe, relying on a very stable algorithm for the calculation of hyperspherical Bessel functions, that can be pushed to arbitrary precision levels. We also introduce a new approximation scheme which gradually takes over in the flat space limit, and significant speeds-up calculations. Our method is implemented in the Boltzmann code CLASS. It can be used to benchmark the accuracy of the CAMB code in curved space, which is found to match expectations. For default precision settings, corresponding to 0.1% for scalar temperature spectra and 0.2% for scalar polarisation spectra, our code is two to three times faster, depending on curvature. We also simplify the temperature and polarisation source terms significantly, so the different contributions to the $C_\ell$'s are easy to identify inside the code.

Dark Matter is essential for structure formation in the late Universe so it must be stable on cosmological time scales. But how stable exactly? Only assuming decays into relativistic particles, we report an otherwise model independent bound on the lifetime of Dark Matter using current cosmological data. Since these decays affect only the low-$\ell$ multipoles of the CMB, the Dark Matter lifetime is expected to correlate with the tensor-to-scalar ratio $r$ as well as curvature $\Omega_k$. We consider two models, including $r$ and $r+\Omega_k$ respectively, versus data from Planck, WMAP, WiggleZ and Baryon Acoustic Oscillations, with or without the BICEP2 data (if interpreted in terms of primordial gravitational waves). This results in a lower bound on the lifetime of CDM given by 160Gyr (without BICEP2) or 200Gyr (with BICEP2) at 95% confidence level.

The apparent discrepancy between the value of the tensor-to-scalar ratio reported by the BICEP2 collaboration, $r = 0.20^{+0.07}_{-0.05}$ at 68% CL, and the Planck upper limit, $r < 0.11$ at 95% CL, has attracted a great deal of attention. In this short note, we show that this discrepancy is mainly due to an `apples to oranges' comparison. The result reported by BICEP2 was measured at a pivot scale $k_* = 0.05$ Mpc$^{-1}$, assuming $n_t = 0$, whereas the Planck limit was provided at $k_* = 0.002$ Mpc$^{-1}$, assuming the slow-roll consistency relation $n_t = -r/8$. One should obviously compare the BICEP2 and Planck results under the same circumstances. By imposing $n_t = 0$, the Planck constraint at $k_* = 0.05$ Mpc$^{-1}$ becomes $r < 0.135$ at $95%$ CL, which can be compared directly with the BICEP2 result. Once a plausible dust contribution to the BICEP2 signal is taken into account (DDM2 model), $r$ is reduced to $r = 0.16^{+0.06}_{-0.05}$ and the discrepancy becomes of order $1.3\sigma$ only.

We present a detailed and self-contained analytical derivation of the evolution of sub-horizon cosmological perturbations before decoupling, based on previous work by S. Weinberg. These solutions are valid in the minimal LCDM scenario, to first order in perturbation theory, in the tight-coupling limit and neglecting neutrino shear stress. We compare them to exact numerical solutions computed by a Boltzmann code, and we find the two to be in very good agreement. The analytic solutions show explicitly that CDM and the baryon-photon fluid effectively behave as separate self-gravitating fluids until the epoch of baryon drag. This in turn leads to the surprising conclusion that the CMB is much less sensitive to the clustering properties of minimally coupled Dark Matter models than what would be naively expected.

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. However, such neutrinos seem incompatible with cosmology because they have too large an impact on cosmic structure formation. Here we show that new interactions in the sterile neutrino sector can prevent their production in the early Universe and reconcile short baseline oscillation experiments with cosmology.

This paper presents the linearised Boltzmann equation for photons for scalar, vector and tensor perturbations in flat, open and closed FLRW cosmologies. We show that E- and B-mode polarisation for all types can be computed using only a single hierarchy. This was previously shown explicitly for tensor modes in flat cosmologies but not for vectors, and not for non-flat cosmologies.

While the cosmic baryon asymmetry has been measured at high accuracy to be 6.1*10^-10, a corresponding lepton asymmetry could be as large as 10^-2 if it hides in the neutrino sector. It has been known for some time that a large asymmetry can be generated from a small initial asymmetry given the existence of a sterile neutrino with a mass less than the mass of the active neutrino. While the magnitude of the final lepton asymmetry is deterministic, its sign has been conjectured to be chaotic in nature. This has been proven in the single momentum approximation, also known as the quantum rate equations, but has up to now not been established using the full momentum dependent quantum kinetic equations. Here we investigate this problem by solving the quantum kinetic equations for a system of 1 active and 1 sterile neutrino on an adaptive grid. We show that by increasing the resolution oscillations in the lepton asymmetry are eliminated and that the sign of the lepton asymmetry is in fact deterministic. This paper also serves as a launch paper for the adaptive solver LASAGNA which is available at http://users-phys.au.dk/steen.

Recent cosmological data favour additional relativistic degrees of freedom beyond the three active neutrinos and photons, often referred to as 'dark' radiation. Light sterile neutrinos is one of the prime candidates for such additional radiation. However, constraints on sterile neutrinos based on the current cosmological data have been derived using simplified assumptions about thermalisation of the sterile neutrino at the Big Bang Nucleosynthesis (BBN) epoch. These assumptions are not necessarily justified and here we solve the full quantum kinetic equations in the (1 active + 1 sterile) scenario and derive the number of thermalised species just before BBN begins (T~1MeV) for null (L=0) and large (L=0.01) initial lepton asymmetry and for a range of possible mass-mixing parameters. We find that the full thermalisation assumption during the BBN epoch is justified for initial small lepton asymmetry only. Partial or null thermalisation occurs when the initial lepton asymmetry is large.

Boltzmann codes are used extensively by several groups for constraining cosmological parameters with Cosmic Microwave Background and Large Scale Structure data. This activity is computationally expensive, since a typical project requires from 10'000 to 100'000 Boltzmann code executions. The newly released code CLASS (Cosmic Linear Anisotropy Solving System) incorporates improved approximation schemes leading to a simultaneous gain in speed and precision. We describe here the three approximations used by CLASS for basic LambdaCDM models, namely: a baryon-photon tight-coupling approximation which can be set to first order, second order or to a compromise between the two; an ultra-relativistic fluid approximation which had not been implemented in public distributions before; and finally a radiation streaming approximation taking reionisation into account.

In the last few years there has been some interest in WIMP Dark Matter models featuring a velocity dependent cross section through the Sommerfeld enhancement mechanism, which is a nonrelativistic effect due to massive bosons in the dark sector. In the first part of this article, we find analytic expressions for the boost factor for three different model potentials, the Coulomb potential, the spherical well and the spherical cone well and compare with the numerical solution of the Yukawa potential. We find that the resonance pattern of all the potentials can be cast into the same universal form. In the second part of the article we perform a detailed computation of the Dark Matter relic density for models having Sommerfeld enhancement by solving the Boltzmann equation numerically. We calculate the expected distortions of the CMB blackbody spectrum from WIMP annihilations and compare these to the bounds set by FIRAS. We conclude that only a small part of the parameter space can be ruled out by the FIRAS observations.