L-PICOLA: A parallel code for fast dark matter simulation

Abstract

Robust measurements based on current large-scale structure surveys require precise knowledge of statistical and systematic errors. This can be obtained from large numbers of realistic mock galaxy catalogues that mimic the observed distribution of galaxies within the survey volume. To this end we present a fast, distributed-memory, planar-parallel code, L-PICOLA, which can be used to generate and evolve a set of initial conditions into a dark matter field much faster than a full non-linear N-Body simulation. Additionally, L-PICOLA has the ability to include primordial non-Gaussianity in the simulation and simulate the past lightcone at run-time, with optional replication of the simulation volume. Through comparisons to fully non-linear N-Body simulations we find that our code can reproduce the z = 0 power spectrum and reduced bispectrum of dark matter to within 2% and 5% respectively on all scales of interest to measurements of Baryon Acoustic Oscillations and Redshift Space Distortions, but 3 orders of magnitude faster. The accuracy, speed and scalability of this code, alongside the additional features we have implemented, make it extremely useful for both current and next generation large-scale structure surveys. L-PICOLA is publicly available at https://cullanhowlett.github.io/l-picola.

The redshift-space momentum power spectrum – II. Measuring the growth rate from the combined 2MTF and 6dFGSv surveys

“…As shown in Howlett et al (2015b), the approximate nature of the simulations compared to a fully non-linear Nbody simulation leads to a suppression of non-linear clus-…”

Assessing non-linear models for galaxy clustering III: theoretical accuracy for stage IV surveys

“…Namely, the TNS model [12] with a Lorentzian damping prefactor and an effective field theory of LSS (EFTofLSS) [13,14] inspired model. Both these models attempt to model the quasi non-linear regime of structure formation in different ways and both were shown to do well in reproducing the measured halo power spectrum from COLA [15][16][17][18][19] simulations when f was fixed to its fiducial value. In that work we were concerned with determining the degeneracies of each model's dofs as well as forecasting the model's constraint on f in the context of a future spectroscopic survey.…”

Cosmological N-body simulations with suppressed variance

“…However, the non-linear evolution of intermediate and large scales is accurately captured (Howlett et al 2015;Koda et al 2016), at a fraction of the computational cost of a traditional N-body simulation. This enables the rapid simulation of extremely large volumes, which in turn allows very precise calculations of different statistics that serve as a benchmark for the performance of our method.…”