Beyond the Lognormal Approximation: a General Simulation Scheme
arXiv:2006.04653 · doi:10.1093/mnras/staa2294
Abstract
We present a public code to generate random fields with an arbitrary probability distribution function (PDF) and an arbitrary correlation function. The algorithm is cosmology-independent, applicable to any stationary stochastic process over a three dimensional grid. We implement it in the case of the matter density field, showing its benefits over the lognormal approximation, which is often used in cosmology for generation of mock catalogues. We find that the covariance of the power spectrum from the new fast realizations is more accurate than that from a lognormal model. As a proof of concept, we also apply the new simulation scheme to the divergence of the Lagrangian displacement field. We find that information from the correlation function and the PDF of the displacement-divergence provides modest improvement over other standard analytical techniques to describe the particle field in the simulation. This suggests that further progress in this direction should come from multi-scale or non-local properties of the initial matter distribution.
14 pages, 17 figures. The code for ITAM is available at https://github.com/tos-1/ITAM . The code for MUSCLE is at https://github.com/tos-1/MUSCLE . Minor modifications to match the accepted MNRAS version
References in corpus (12)
- The DESI Experiment Part I: Science,Targeting, and Survey Design
- Resumming Cosmological Perturbations via the Lagrangian Picture: One-loop Results in Real Space and in Redshift Space
- Analytic model for the matter power spectrum, its covariance matrix, and baryonic effects
- nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics
- VOBOZ: An Almost-Parameter-Free Halo-Finding Algorithm
- Generating Log-normal Mock Catalog of Galaxies in Redshift Space
- Bayesian analysis of the dynamic cosmic web in the SDSS galaxy survey
- Quantifying distortions of the Lagrangian dark-matter mesh in cosmology
- New fitting formula for cosmic non-linear density distribution
- Precision Prediction of the Log Power Spectrum
- Removable Matter-Power-Spectrum Covariance from Bias Fluctuations
- Statistical connection of peak counts to power spectrum and moments in weak lensing field