Normalized to: Dakin, J.

**5** article(s) in total. **6** co-authors, from **1** to **4** common article(s). Median position in authors list is **1,0**.

[1] oai:arXiv.org:1904.11773 [pdf] - 1901906

Fully relativistic treatment of decaying cold dark matter in $N$-body
simulations

Comments:

Submitted: **2019-04-26**

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.

[2] oai:arXiv.org:1904.05210 [pdf] - 1945760

Dark energy perturbations in $N$-body simulations

Comments:

Submitted: **2019-04-10**

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.

[3] oai:arXiv.org:1811.00904 [pdf] - 1851577

Fully relativistic treatment of light neutrinos in $N$-body simulations

Comments: 11 pages, 4 figures

Submitted: **2018-11-02**

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.

[4] oai:arXiv.org:1712.03944 [pdf] - 1842294

$\nu$CO$N$CEPT: Cosmological neutrino simulations from the non-linear
Boltzmann hierarchy

Comments: 28 pages, 10 figures

Submitted: **2017-12-11**

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.

[5] oai:arXiv.org:1510.07621 [pdf] - 1300429

CONCEPT - The COsmological $N$-body CodE in PyThon

Comments: 34 pages

Submitted: **2015-10-26**

CO$N$CEPT (COsmological $N$-body CodE in PyThon) is a free and open-source
code for cosmological $N$-body simulations on massively parallel computers with
distributed memory. Collisionless dark matter is the only implemented particle
species. Gravity can be computed using the PP, PM or the P$^{3}$M algorithm.
The goal of CO$N$CEPT is to make it pleasant to work with cosmological $N$-body
simulations - for the cosmologist as well as for the source code developer.
This is the user guide. The source code and additional documentation can be
found at https://github.com/jmd-dk/concept/