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 ΛCDM cosmologies with massless neutrinos to include the effects of massive, but light neutrinos, where the effect of neutrinos on the matter clustering can be adequately modelled using linear perturbation theory for 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.
Read full abstract