A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos

Bayer, Adrian E.; Banerjee, Arka; Feng, Yu

doi:10.1088/1475-7516/2021/01/016

Cited by 32 publications

(31 citation statements)

References 108 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…As for COLA, the order of convergence has, to our knowledge, not been discussed in the literature. The FASTPM approach has been recently extended to include the modelling of massive neutrinos (Bayer et al 2021), and also ported to TENSORFLOW (Modi et al 2020b).…”

Section: Fastpmmentioning

confidence: 99%

“…4.4 for a review of such approaches). Recently, Bayer et al (2021) proposed an extension to the FASTPM algorithm by modifying the kick and drift operator; and Wright et al (2017) extended the COLA algorithm by incorporating the scaledependent effect of neutrinos in the growth factors. On the other hand, an accurate prediction for the impact of neutrinos can also be achieved in post-processing by using cosmology-rescaling methods (Zennaro et al 2019;Angulo and White 2010;Contreras et al 2020c), or as a gauge transformation (Partmann et al 2020).…”

Section: Massive Neutrinosmentioning

confidence: 99%

See 1 more Smart Citation

Large-scale dark matter simulations

Angulo

Hahn

2022

Living Rev Comput Astrophys

105

View full text Add to dashboard Cite

We review the field of collisionless numerical simulations for the large-scale structure of the Universe. We start by providing the main set of equations solved by these simulations and their connection with General Relativity. We then recap the relevant numerical approaches: discretization of the phase-space distribution (focusing on N-body but including alternatives, e.g., Lagrangian submanifold and Schrödinger–Poisson) and the respective techniques for their time evolution and force calculation (direct summation, mesh techniques, and hierarchical tree methods). We pay attention to the creation of initial conditions and the connection with Lagrangian Perturbation Theory. We then discuss the possible alternatives in terms of the micro-physical properties of dark matter (e.g., neutralinos, warm dark matter, QCD axions, Bose–Einstein condensates, and primordial black holes), and extensions to account for multiple fluids (baryons and neutrinos), primordial non-Gaussianity and modified gravity. We continue by discussing challenges involved in achieving highly accurate predictions. A key aspect of cosmological simulations is the connection to cosmological observables, we discuss various techniques in this regard: structure finding, galaxy formation and baryonic modelling, the creation of emulators and light-cones, and the role of machine learning. We finalise with a recount of state-of-the-art large-scale simulations and conclude with an outlook for the next decade.

show abstract

Section: Fastpmmentioning

confidence: 99%

Section: Massive Neutrinosmentioning

confidence: 99%

Large-scale dark matter simulations

Angulo

Hahn

2022

Living Rev Comput Astrophys

105

View full text Add to dashboard Cite

show abstract

“…While the default version of FastPM uses an extra particle species to simulate massive neutrinos (Bayer et al 2021c), which are labelled as ncdm (not-cold dark matter), we use a modified version to treat ncdm same as that in AbacusSummit, i.e., the effect of massive neutrinos only contribute to the Hubble expansion rate but not to the gravitational forces from clustering.…”

Section: Configuration Of Desi-fastpm Simulationsmentioning

confidence: 99%

The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations

Ding,

Chuang,

et al. 2022

Preprint

Self Cite

View full text Add to dashboard Cite

Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise 3D map of our Universe. The survey effective volume reaches ∼ 20 Gpc 3 ℎ −3 . It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. AbacusSummit is a suite of high-resolution dark-matter-only simulations designed for this purpose, with 200 Gpc 3 ℎ −3 (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-𝑁-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variances than the expected DESI statistical variance at the scales we are interested in, e.g., 𝑘 < 0.3 ℎMpc −1 . We further generalize the method for other cosmologies with only one realization in AbacusSummit suite to extend the effective volume ∼ 20 times. In summary, our proposed strategy of combining high fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

show abstract

“…Furthermore, there are even attempts to solve the 6D Vlasov-Poisson equation directly [31]. A very recent development was done in [32] where they successfully realised neutrinos as particles in the quasi-N-body code FastPM. Furthermore, there are several approximative methods, such the cosmological rescaling method [33,34].…”

Section: Neutrino Physicsmentioning

confidence: 99%

A Minimal Model for Massive Neutrinos in Newtonian N-body Simulations

Heuschling¹,

Partmann²,

Fidler³

2022

Preprint

View full text Add to dashboard Cite

We present a novel method for including the impact of massive neutrinos in cold dark matter N-body simulations. Our approach is compatible with widely employed Newtonian N-body codes and relies on only three simple modifications. First, we use commonly employed backscaling initial conditions, based on the cold dark matter plus baryon power spectrum instead of the total matter power spectrum. Second, the accurate Hubble rate is employed in both the backscaling and the evolution of particles in the N-body code. Finally, we shift the final particle positions in a post-processing step to account for the integrated effect of neutrinos on the particles in the simulation. However, we show that the first two modifications already capture most of the relevant neutrino physics for a large range of observationally interesting redshifts and scales. The output of the simulations are the cold dark matter and baryon distributions and can be analysed using standard methods. All modifications are simple to implement and do not generate any computational overhead. By implementing our methods in the N-body codes gadget-4 and gevolution, we show that any state-of-the-art Newtonian N-body code can be utilised out of the box. Our method is also compatible with higher order Lagrangian perturbation theory initial conditions and accurate for masses up to at least m ν = 0.3 eV. Being formulated in relativistic gauge theory, in addition to including the impact of massive neutrinos, our method further includes relativistic corrections relevant on the large scales for free.

show abstract

A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos

Cited by 32 publications

References 108 publications

Large-scale dark matter simulations

Large-scale dark matter simulations

The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations

A Minimal Model for Massive Neutrinos in Newtonian N-body Simulations

Contact Info

Product

Resources

About