Accurate cosmological simulations that include the effect of non-linear matter clustering as well as of massive neutrinos are essential for measuring the neutrino mass scale from upcoming galaxy surveys. Typically, Newtonian simulations are employed and the neutrino distribution is sampled with a large number of particles in order to beat down the shot noise. Here we perform very efficient simulations with light massive neutrinos which require virtually no extra cost over matter-only simulations, and furthermore do not require tracer particles for the neutrinos. Instead, we use a weak-field dictionary based on the recently developed Newtonion motion approach, where Newtonian simulations for matter are paired with a linear relativistic Boltzmann code to allow for an absorption of the neutrino evolution into a time-dependent coordinate transformation. For this, only minimal modifications in existing N-body codes are required, which we have explicitly implemented in $\textit{gevolution}$ and $\textit{gadget-2}$. Our fast method determines the non-linear matter power spectrum to permille-level precision when compared against state-of-the-art simulations that have been performed for $0.1\,{\rm eV} \leq \sum m_\nu \leq 0.3\,$eV.

Reionisation blurring is a non-linear correction to the cosmic microwave background that acts similar to weak gravitational lensing and that can be computed from linear perturbations through a blurring potential. Its impact on the cosmic microwave background is roughly two order of magnitude smaller than that of lensing, in isolation. But the blurring potential is strongly correlated with the lensing potential thereby generating a potentially observable cross-correlation. We compute for the first time the inclusive impact of reionisation blurring on the temperature angular power spectrum and discuss how much it could induce lensing anomalies

We clarify the subtle issue of finding the correct mapping of Newtonian simulations to light-cone observables at very large distance scales. A faithful general-relativistic interpretation specifies a gauge, i.e. a chart that relates the simulation data to points of the space-time manifold. It has already been pointed out that the implicit gauge choice of Newtonian simulations is indeed different from the Poisson gauge that is commonly adopted for relativistic calculations, the difference being most significant at large scales. It is therefore inconsistent, for example, to predict weak-lensing observables from simulations unless this gauge issue is properly accounted for. Using perturbation theory as well as fully relativistic N-body simulations we quantify the systematic error introduced this way, and we discuss several solutions that would render the calculations relativistically self-consistent.

Einstein-Boltzmann Solvers (EBSs) are run on a massive scale by the cosmology community when fitting cosmological models to data. We present a new concept for speeding up such codes with neural networks. The originality of our approach stems from not substituting the whole EBS by a machine learning algorithm, but only its most problematic and least parallelizable step: the integration of perturbation equations over time. This approach offers two significant advantages: the task depends only on a subset of cosmological parameters, and it is blind to the characteristics of the experiment for which the output must be computed (for instance, redshift bins). These allow us to construct a fast and highly re-usable network. In this proof-of-concept paper, we focus on the prediction of CMB source functions, and design our networks according to physical considerations and analytical approximations. This allows us to reduce the depth and training time of the networks compared to a brute-force approach. Indeed, the calculation of the source functions using the networks is fast enough so that it is not a bottleneck in the EBS anymore. Finally, we show that their accuracy is more than sufficient for accurate MCMC parameter inference from Planck data. This paves the way for a new project, CosmicNet, aimed at gradually extending the use and the range of validity of neural networks within EBSs, and saving massive computational time in the context of cosmological parameter extraction.

Observables of cosmic structures are usually not the underlying matter field but biased tracers of matter, such as galaxies or halos. We show how the bias found in Newtonian N-body simulations can be interpreted in terms of the weak-field limit of General Relativity (GR). For this we employ standard Newtonian simulations of cold dark matter and incorporate GR/radiation via a weak-field dictionary that we have recently developed. We find that even when a simple local biasing scheme is employed in the Newtonian simulation, the relativistic bias becomes inherently scale-dependent due to the presence of radiation and GR corrections. This scale-dependence could be in principle observed on large scales in upcoming surveys. As a working example, we apply our methodology to Newtonian simulations for the spherical collapse and recover permille-level agreement between the approaches for extracting the relativistic bias on all considered scales.

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.

Halos and galaxies are tracers of the underlying dark matter structures. While their bias is well understood in the case of a simple Universe composed dominantly of dark matter, the relation becomes more complex in the presence of massive neutrinos. Indeed massive neutrinos introduce rich dynamics in the process of structure formation leading to scale-dependent bias. We study this process from the perspective of general relativity employing a simple spherical collapse model. We find a characteristic signature at the neutrino free-streaming scale in addition to a large-scale feature from general relativity. The scale-dependent halo bias opposes the suppression in the matter distribution due to neutrino free-streaming and leads to corrections of a few percent in the halo power spectrum. It is not only sensitive to the sum of the neutrino-masses, but respond to the individual masses. Accurate models for the neutrino bias are a crucial ingredient for the future data analysis and play an important role in constraining the neutrino masses.

Initial conditions for cosmological N-body simulations are usually calculated by rescaling the present day linear power spectrum obtained from an Einstein-Boltzmann solver to the initial time employing the scale-independent matter growth function. For the baseline Lambda-CDM model, this has been shown to be consistent with General Relativity (GR) even in the presence of relativistic species such as photons. We show that this approach is not feasible in cosmologies with massive neutrinos and present an alternative method employing the Newtonian motion gauge framework.

Anisotropies in the Cosmic Microwave Background (CMB) can be induced during the later stages of cosmic evolution, and in particular during and after the Epoch of Reionisation. Inhomogeneities in the ionised fraction, but also in the baryon density, in the velocity fields and in the gravitational potentials are expected to generate correlated CMB perturbations. We present a complete relativistic treatment of all these effects, up to second order in perturbation theory, that we solve using the numerical Boltzmann code SONG. The physical origin and relevance of all second order terms are carefully discussed. In addition to collisional and gravitational contributions, we identify the diffuse analogue of the blurring and kinetic Sunyaev-Zel'dovich (SZ) effects. Our approach naturally includes the correlations between the imprint from patchy reionisation and the diffuse SZ effects thereby allowing us to derive reliable estimates of the induced temperature and polarisation CMB angular power spectra. In particular, we show that the B-modes generated at intermediate length-scales (l~100) have the same amplitude as the B-modes coming from primordial gravitational waves with a tensor-to-scalar ratio r=10^{-4}.

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.

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.

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 build a statistical description of fermions, taking into account the spin degree of freedom in addition to the momentum of particles, and we detail its use in the context of the kinetic theory of gases of fermions particles. We show that the one-particle distribution function needed to write a Liouville equation is a spinor valued operator. The degrees of freedom of this function are covariantly described by an intensity function and by a polarisation vector which are parallel transported by free streaming. Collisions are described on the microscopic level and lead to a Boltzmann equation for this operator. We apply our formalism to the case of weak interactions, which at low energies can be considered as a contact interaction between fermions, allowing us to discuss the structure of the collision term for a few typical weak-interaction mediated reactions. In particular we find for massive particles that a dipolar distribution of velocities in the interacting species is necessary to generate linear polarisation, as opposed to the case of photons for which linear polarisation is generated from the quadrupolar distribution of velocities.

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.

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.

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 investigate the generation of magnetic fields from non-linear effects around recombination. As tight-coupling is gradually lost when approaching $z\simeq 1100$, the velocity difference between photons and baryons starts to increase, leading to an increasing Compton drag of the photons on the electrons. The protons are then forced to follow the electrons due to the electric field created by the charge displacement; the same field, following Maxwell's laws, eventually induces a magnetic field on cosmological scales. Since scalar perturbations do not generate any magnetic field as they are curl-free, one has to resort to second-order perturbation theory to compute the magnetic field generated by this effect. We reinvestigate this problem numerically using the powerful second-order Boltzmann code SONG. We show that: i) all previous studies do not have a high enough angular resolution to reach a precise and consistent estimation of the magnetic field spectrum; ii) the magnetic field is generated up to $z\simeq 10$; iii) it is in practice impossible to compute the magnetic field with a Boltzmann code for scales smaller than $1\,{\rm Mpc}$. Finally we confirm that for scales of a few ${\rm Mpc}$, this magnetic field is of order $2\times 10^{-29}{\rm G}$, many orders of magnitude smaller than what is currently observed on intergalactic scales.

We develop the transport operator formalism, a new line-of-sight integration framework to calculate the anisotropies of the Cosmic Microwave Background (CMB) at the linear and non-linear level. This formalism utilises a transformation operator that removes all inhomogeneous propagation effects acting on the photon distribution function, thus achieving a split between perturbative collisional effects at recombination and non-perturbative line-of-sight effects at later times. The former can be computed in the framework of standard cosmological perturbation theory with a second-order Boltzmann code such as SONG, while the latter can be treated within a separate perturbative scheme allowing the use of non-linear Newtonian potentials. We thus provide a consistent framework to compute all physical effects contained in the Boltzmann equation and to combine the standard remapping approach with Boltzmann codes at any order in perturbation theory, without assuming that all sources are localised at recombination.

In this note, we study cosmic microwave background (CMB) constraints on primordial non-Gaussianity in DBI galileon models in which an induced gravity term is added to the Dirac-Born-Infeld (DBI) action. In this model, the non-Gaussianity of orthogonal shape can be generated. We provide a relation between theoretical parameters and orthogonal/equilateral non-linear parameters using the Fisher matrix approach for the CMB bispectrum. In doing so, we include the effect of the CMB transfer functions and experimental noise properties by employing the recently developed \textsf{SONG} code. The relation is also shown in the language of effective theory so that it can be applied to general single-field models. Using the bispectrum Fisher matrix and the central values for equilateral and orthogonal non-Gaussianities found by the Planck temperature survey, we provide forecasts on the theoretical parameters of the DBI galileon model. We consider the upcoming Planck polarisation data and the proposed post-Planck experiments COrE and PRISM. We find that Planck polarisation measurements may provide a hint for a non-canonical sound speed at the 68% confidence level. COrE and PRISM will not only confirm a non-canonical sound speed but also exclude the conventional DBI inflation model at more than the 95% and 99% confidence level respectively, assuming that the central values will not change. This indicates that improving constraints on non-Gaussianity further by future CMB experiments is invaluable to constrain the physics of the early universe.

We compute the bispectrum induced in the cosmic microwave background (CMB) temperature and polarisation by the evolution of the primordial density perturbations using the second-order Boltzmann code SONG. We show that adding polarisation increases the signal-to-noise ratio by a factor four with respect to temperature alone and we estimate the observability of this intrinsic bispectrum and the bias it induces on measurements of primordial non-Gaussianity. When including all physical effects except the late-time non-linear evolution, we find for the intrinsic bispectrum a signal-to-noise of $S/N=3.8,\,2.9,\,1.6$ and $0.5$ for, respectively, an ideal experiment with an angular resolution of $\ell_\text{max}=3000$, the proposed CMB surveys PRISM and COrE, and Planck's polarised data; the bulk of this signal comes from the $E$-polarisation and from squeezed configurations. We discuss how CMB lensing is expected to reduce these estimates as it suppresses the bispectrum for squeezed configurations and contributes to the noise in the estimator. We find that the presence of the intrinsic bispectrum will bias a measurement of primordial non-Gaussianity of local type by $f_{NL}^\text{bias}=0.66$ for an ideal experiment with $\ell_\text{max}=3000$. Finally, we verify the robustness of our results by reproducing the analytical approximation for the squeezed-limit bispectrum in the general polarised case.

We estimate the B-polarisation induced in the Cosmic Microwave Background by the non-linear evolution of density perturbations. Using the second-order Boltzmann code SONG, our analysis incorporates, for the first time, all physical effects at recombination. We also include novel contributions from the redshift part of the Boltzmann equation and from the bolometric definition of the temperature in the presence of polarisation. The remaining line-of-sight terms (lensing and time-delay) have previously been studied and must be calculated non-perturbatively. The intrinsic B-mode polarisation is present independent of the initial conditions and might contaminate the signal from primordial gravitational waves. We find this contamination to be comparable to a primordial tensor-to-scalar ratio of $r\simeq10^{-7}$ at the angular scale $\ell\simeq100\,$, where the primordial signal peaks, and $r\simeq 5 \cdot 10^{-5}$ at $\ell\simeq700\,$, where the intrinsic signal peaks. Therefore, we conclude that the intrinsic B-polarisation from second-order effects is not likely to contaminate future searches of primordial gravitational waves.

We compute the spectral distortions of the Cosmic Microwave Background (CMB) polarization induced by non-linear effects in the Compton interactions between CMB photons and cold intergalactic electrons. This signal is of the $y$-type and is dominated by contributions arising from the reionized era. We stress that it is not shadowed by the thermal SZ effect which has no equivalent for polarization. We decompose its angular dependence into $E$- and $B$-modes, and we calculate the corresponding power spectra, both exactly and using a suitable Limber approximation that allows a simpler numerical evaluation. We find that $B$-modes are of the same order of magnitude as $E$-modes. Both spectra are relatively flat, peaking around $\ell=280$, and their overall amplitude is directly related to the optical depth to reionization. Moreover, we find this effect to be one order of magnitude larger than the non-linear kinetic Sunyaev-Zel'dovich effect in galaxy clusters. Finally, we discuss how to improve the detectability of our signal by cross-correlating it with other quantities sourced by the flow of intergalactic electrons.

We develop a new, efficient code for solving the second-order Einstein-Boltzmann equations, and use it to estimate the intrinsic CMB non-Gaussianity arising from the non-linear evolution of density perturbations. The full calculation involves contributions from recombination and less tractable contributions from terms integrated along the line of sight. We investigate the bias that this intrinsic bispectrum implies for searches of primordial non-Gaussianity. We find that the inclusion or omission of certain line of sight terms can make a large impact. When including all physical effects but lensing and time-delay, we find that the local-type f_nl would be biased by f_nl ~ 0.5, below the expected sensitivity of the Planck satellite. The speed of our code allows us to confirm the robustness of our results with respect to a number of numerical parameters.

B-mode polarization of the cosmic background radiation is induced from purely scalar primordial sources at second order in perturbations of the homogeneous, isotropic universe. We calculate the B-mode angular power spectrum C_l^{BB} sourced by the second-order scattering term in the full second-order Boltzmann equations for the polarized radiation phase-space density, which have recently become available. We find that at l\approx 200 the second-order effect is comparable to the first-order effect for a tensor-to-scalar ratio of r=10^{-6}, and to about 2\cdot 10^{-4} at l\approx 1000. It is always negligible relative to the weak-lensing induced contribution.

Within the Closed Time Path (CTP) framework, we derive kinetic equations for particle distribution functions that describe leptogenesis in the presence of several lepton flavours. These flavours have different Standard-Model Yukawa couplings, which induce flavour-sensitive scattering processes and thermal dispersion relations. Kinetic equilibrium, which is rapidly established and maintained via gauge interactions, allows to simplify these equations to kinetic equations for the matrix of lepton charge densities. In performing this simplification, we notice that the rapid flavour-blind gauge interactions damp the flavour oscillations of the leptons. Leptogenesis turns out to be in the parametric regime where the flavour oscillations are overdamped and flavour decoherence is mainly induced by flavour sensitive scatterings. We solve the kinetic equations for the lepton number densities numerically and show that they interpolate between the unflavoured and the fully flavoured regimes within the intermediate parametric region, where neither of these limits is applicable.

Non-gaussianity and B-mode polarization are particularly interesting features of the cosmic microwave background, as -- at least in the standard model of cosmology -- their only sources to first order in cosmological perturbation theory are primordial, possibly generated during inflation. If the primordial sources are small, the question arises how large is the non-gaussianity and B-mode background induced in second-order from the initially gaussian and scalar perturbations. In this paper we derive the Boltzmann hierarchy for the microwave background photon phase-space distributions at second order in cosmological perturbation theory including the complete polarization information, providing the basis for further numerical studies. As an aside we note that the second-order collision term contains new sources of B-mode polarization and that no polarization persists in the tight-coupling limit.