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.

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.

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.

We have calculated the non-linear effects of generic fermionic and bosonic hot dark matter components in cosmological N-body simulations. For sub-eV masses, the non-linear power spectrum suppression caused by thermal free-streaming resembles the one seen for massive neutrinos, whereas for masses larger than 1eV, the non-linear relative suppression of power is smaller than in linear theory. We furthermore find that in the non-linear regime, one can map fermionic to bosonic models by performing a simple transformation.

We perform a Fourier space decomposition of the dynamics of non-linear cosmological structure formation in LCDM models. From N-body simulations involving only cold dark matter we calculate 3-dimensional non-linear density, velocity divergence and vorticity Fourier realizations, and use these to calculate the fully non-linear mode coupling integrals in the corresponding fluid equations. Our approach allows for a reconstruction of the amount of mode coupling between any two wavenumbers as a function of redshift. With our Fourier decomposition method we identify the transfer of power from larger to smaller scales, the stable clustering regime, the scale where vorticity becomes important, and the suppression of the non-linear divergence power spectrum as compared to linear theory. Our results can be used to improve and calibrate semi-analytical structure formation models.

We study the expected variance of measurements of the Hubble constant, $H_0$, as calculated in either linear perturbation theory or using non-linear velocity power spectra derived from $N$-body simulations. We compare the variance with that obtained by carrying out mock observations in the N-body simulations, and show that the estimator typically used for the local Hubble constant in studies based on perturbation theory is different from the one used in studies based on N-body simulations. The latter gives larger weight to distant sources, which explains why studies based on N-body simulations tend to obtain a smaller variance than that found from studies based on the power spectrum. Although both approaches result in a variance too small to explain the discrepancy between the value of $H_0$ from CMB measurements and the value measured in the local universe, these considerations are important in light of the percent determination of the Hubble constant in the local universe.

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.

We use N-body simulations to find the effect of neutrino masses on halo properties, and investigate how the density profiles of both the neutrino and the dark matter components change as a function of the neutrino mass. We compare our neutrino density profiles with results from the N-one-body method and find good agreement. We also show and explain why the Tremaine-Gunn bound for the neutrinos is not saturated. Finally we study how the halo mass function changes as a function of the neutrino mass and compare our results with the Sheth-Tormen semi-analytic formulae. Our results are important for surveys which aim at probing cosmological parameters using clusters, as well as future experiments aiming at measuring the cosmic neutrino background directly.

The cosmic neutrino background anisotropy is calculated for massive neutrino states by solving the full Boltzmann equation. The effect of weak gravitational lensing, including the Limber approximation, is also derived for massive particles, and subsequently applied to the case of massive neutrinos.

We present the first simulation capable of resolving the structure of neutrino clustering on Mpc scales. The method combines grid- and particle-based methods and achieves very good accuracy on both small and large scales, while keeping CPU consumption under control. Such simulations are not only ideal for calculating the non-linear matter power spectrum but also particularly relevant for studies of how neutrinos cluster in galaxy- or cluster-sized halos. We perform the largest neutrino N-body simulation to date, effectively containing 10 different neutrino hot dark matter components with different thermal properties.

We present a novel, fast and precise method for including the effect of light neutrinos in cosmological N-body simulations. The effect of the neutrino component is included by using the linear theory neutrino perturbations in the calculation of the gravitational potential in the N-body simulation. By comparing this new method with the full non-linear evolution first presented in \cite{Brandbyge1}, where the neutrino component was treated as particles, we find that the new method calculates the matter power spectrum with an accuracy better than 1% for \sum m_\nu \lesssim 0.5 eV at z = 0. This error scales approximately as (\sum m_\nu)^2, making the new linear neutrino method extremely accurate for a total neutrino mass in the range 0.05 - 0.3 eV. At z = 1 the error is below 0.3% for \sum m_\nu \lesssim 0.5 eV and becomes negligible at higher redshifts. This new method is computationally much more efficient than representing the neutrino component by N-body particles.

We have performed detailed studies of non-linear structure formation in cosmological models with light neutrinos. For the first time the effect of neutrino thermal velocities has been included in a consistent way, and the effect on the matter power spectrum is found to be significant. The effect is large enough to be measured in future, high precision surveys. Additionally, we provide a simple but accurate analytic expression for the suppression of fluctuation power due to massive neutrinos. Finally, we describe a simple and fast method for including the effect of massive neutrinos in large-scale N-body simulations which is accurate at the 1% level for \sum m_nu < 0.15 eV.

The extragalactic cosmic gamma-ray background (CGB) is an interesting channel to look for signatures of dark matter annihilation. In particular, besides the imprint in the energy spectrum, peculiar anisotropy patterns are expected compared to the case of a pure astrophysical origin of the CGB. We take into account the uncertainties in the dark matter clustering properties on sub-galactic scales, deriving two possible anisotropy scenarios. A clear dark matter angular signature is achieved when the annihilation signal receives only a moderate contribution from sub-galactic clumps and/or cuspy haloes. Experimentally, if galactic foregrounds systematics are efficiently kept under control, the angular differences are detectable with the forthcoming GLAST observatory, provided that the annihilation signal contributes to the CGB for a fraction >10-20%. If, instead, sub-galactic structures have a more prominent role, the astrophysical and dark matter anisotropies become degenerate, correspondingly diluting the DM signature. As complementary observables we also introduce the cross-correlation between surveys of galaxies and the CGB and the cross-correlation between different energy bands of the CGB and we find that they provide a further sensitive tool to detect the dark matter angular signatures.

We have observed oscillations in the nearby G2 subgiant star beta Hyi using high-precision velocity observations obtained over more than a week with the HARPS and UCLES spectrographs. The oscillation frequencies show a regular comb structure, as expected for solar-like oscillations, but with several l=1 modes being strongly affected by avoided crossings. The data, combined with those we obtained five years earlier, allow us to identify 28 oscillation modes. By scaling the large frequency separation from the Sun, we measure the mean density of beta Hyi to an accuracy of 0.6%. The amplitudes of the oscillations are about 2.5 times solar and the mode lifetime is 2.3 d. A detailed comparison of the mixed l=1 modes with theoretical models should allow a precise estimate of the age of the star.