Accurate estimators of power spectra in N-body simulations

Abstract

abridged] A method to rapidly estimate the Fourier power spectrum of a point distribution is presented. This method relies on a Taylor expansion of the trigonometric functions. It yields the Fourier modes from a number of FFTs, which is controlled by the order N of the expansion and by the dimension D of the system. In three dimensions, for the practical value N=3, the number of FFTs required is 20. We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive explicit analytic expression for the spectrum, which allows us to quantify--and correct for--the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N. The only remaining significant source of errors is reduced to the unavoidable cosmic/sample variance due to the finite size of the sample. The analytical calculations are successfully checked against a cosmological N-body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier-Taylor estimator behaves well. We also show how to reach arbitrarily large dynamic range in Fourier space (i.e., high wavenumber), while keeping statistical errors in control, by appropriately "folding" the particle distribution.

arXiv comment18 Pages, 9 Figures. Accepted for publication in MNRAS. The Fourier-Taylor module as well as the associated power spectrum estimator tool we propose is available as an F90 package, POWMES, at http://www.projet-horizon.fr or on request from the authors

References in corpus (25)

• Stable clustering, the halo model and nonlinear cosmological power spectra
• GADGET: A code for collisionless and gasdynamical cosmological simulations
• Analytic model for galaxy and dark matter clustering
• Halo occupation numbers and galaxy bias
• Power Spectrum Analysis of Three-Dimensional Redshift Surveys
• Nonlinear evolution of cosmological power spectra
• How Many Galaxies Fit in a Halo? Constraints on Galaxy Formation Efficiency from Spatial Clustering
• Transients from Initial Conditions in Cosmological Simulations
• Evolution of structure in cold dark matter universes
• Very weak lensing in the CFHTLS Wide: Cosmology from cosmic shear in the linear regime
• Multiscale Gaussian Random Fields for Cosmological Simulations
• Deriving the Nonlinear Cosmological Power Spectrum and Bispectrum from Analytic Dark Matter Halo Profiles and Mass Functions
• Non-Linear Evolution of the Bispectrum of Cosmological Perturbations
• Correcting for the alias effect when measuring the power spectrum using FFT
• The Power Spectrum of Mass Fluctuations Measured from the Lyman-alpha Forest at Redshift z=2.5
• Power Spectrum Correlations Induced by Non-Linear Clustering
• Cosmological Constraints From the 100 Square Degree Weak Lensing Survey
• On measuring the covariance matrix of the nonlinear power spectrum from simulations
• Cosmology and the Halo Occupation Distribution from Small-Scale Galaxy Clustering in the Sloan Digital Sky Survey
• Information content of the non-linear power spectrum: the effect of beat-coupling to large scales
• Point processes and stochastic displacement fields
• Linear perturbative theory of the discrete cosmological N-body problem
• Quantification of discreteness effects in cosmological N-body simulations: I. Initial Conditions
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
• Cosmic Statistics of Statistics: N-point Correlations

Cited by in corpus (81)

• First results from the IllustrisTNG simulations: matter and galaxy clustering
• Cosmology and Fundamental Physics with the Euclid Satellite
• Cosmology and fundamental physics with the Euclid satellite
• The effects of galaxy formation on the matter power spectrum: A challenge for precision cosmology
• The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum
• ECOSMOG: An Efficient Code for Simulating Modified Gravity
• Accurate Estimators of Correlation Functions in Fourier Space
• N-body Simulations for f(R) Gravity using a Self-adaptive Particle-Mesh Code
• Large-scale dark matter simulations
• Exploring the effects of galaxy formation on matter clustering through a library of simulation power spectra
• Matter power spectrum and the challenge of percent accuracy
• Planck early results. VI. The High Frequency Instrument data processing
• Modified Gravity N-body Code Comparison Project
• N-Body Simulations for Coupled Scalar Field Cosmology
• The nonlinear matter and velocity power spectra in f(R) gravity
• Combining perturbation theories with halo models
• Full-Sky Weak Lensing Simulation with 70 Billion Particles
• Redshift-weighted constraints on primordial non-Gaussianity from the clustering of the eBOSS DR14 quasars in Fourier space
• On the road to percent accuracy: nonlinear reaction of the matter power spectrum to dark energy and modified gravity
• Structure Formation in the Symmetron Model
• Exploring Vainshtein mechanism on adaptively refined meshes
• Systematic simulations of modified gravity: symmetron and dilaton models
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• Simulating the quartic Galileon gravity model on adaptively refined meshes
• The CoDECS project: a publicly available suite of cosmological N-body simulations for interacting dark energy models
• Nonlinear Structure Formation with the Environmentally Dependent Dilaton
• Testing the theory of gravity with DESI: estimators, predictions and simulation requirements
• Imprints of Dark Energy on Cosmic Structure Formation I) Realistic Quintessence Models and the Non-Linear Matter Power Spectrum
• Computing General Relativistic effects from Newtonian N-body simulations: Frame dragging in the post-Friedmann approach
• Halo model and halo properties in Galileon gravity cosmologies
• Cosmic Web and Environmental Dependence of Screening: Vainshtein vs. Chameleon
• Higher-order initial conditions for mixed baryon-CDM simulations
• Emulators for the non-linear matter power spectrum beyond $\Lambda$CDM
• Modeling the nonlinear clustering in modified gravity models I: A fitting formula for matter power spectrum of f(R) gravity
• The contributions of matter inside and outside of haloes to the matter power spectrum
• Simultaneous modelling of matter power spectrum and bispectrum in the presence of baryons
• Speeding up $N$-body simulations of modified gravity: Chameleon screening models
• Structure Formation by the Fifth Force III: Segregation of Baryons and Dark Matter
• The effect of baryons on the variance and the skewness of the mass distribution in the universe at small scales
• Unbiased Cosmology Inference from Biased Tracers using the EFT Likelihood
• Structure formation simulations with momentum exchange: alleviating tensions between high-redshift and low-redshift cosmological probes
• Speeding up N-body simulations of modified gravity: Vainshtein screening models
• Simulating Momentum Exchange in the Dark Sector
• The Non-Linear Matter Power Spectrum in Warm Dark Matter Cosmologies
• GRAMSES: a new route to general relativistic $N$-body simulations in cosmology. Part I. Methodology and code description
• Effects of the initial conditions on cosmological $N$-body simulations
• On the Effects of Coupled Scalar Fields on Structure Formation
• Combining perturbation theories with halo models for the matter bispectrum
• Lyman-α Forest and Cosmic Weak Lensing in a Warm Dark Matter Universe
• The effect of photometric redshift uncertainties on galaxy clustering and baryonic acoustic oscillations
• The fully non-linear post-Friedmann frame-dragging vector potential: Magnitude and time evolution from N-body simulations
• Cosmological gravity on all scales II: Model independent modified gravity $N$-body simulations
• Weak lensing predictions for coupled dark energy cosmologies at non-linear scales
• The velocity field in MOND cosmology
• Fitting and forecasting non-linear coupled dark energy
• A Semi-analytic Ray-tracing Algorithm for Weak Lensing
• The Effect of Thermal Velocities on Structure Formation in N-body Simulations of Warm Dark Matter
• Self-Similar Bumps and Wiggles: Isolating the Evolution of the BAO Peak with Power-law Initial Conditions
• Revisiting the matter power spectra in $f(R)$ gravity
• Revisiting the screening mechanism in $f(R)$ gravity
• Decaying warm dark matter and structure formation
• GRAMSES: a new route to general relativistic $N$-body simulations in cosmology. Part II. Initial conditions
• Improved analytical modeling of the non-linear power spectrum in modified gravity cosmologies
• Effective dark matter power spectra in $f(R)$ gravity
• Proca-stinated Cosmology I: A N-body code for the vector Galileon
• Cosmological gravity on all scales III: non-linear matter power spectrum in phenomenological modified gravity
• Stable clustering and the resolution of dissipationless cosmological N-body simulations
• Nonlinear structure formation in Bound Dark Energy
• A scalar field dark matter model and its role in the large scale structure formation in the Universe
• The Impact of Baryons on the Large-Scale Structure of the Universe
• A 1% accurate method to include baryonic effects in galaxy-galaxy lensing models
• Proca-stinated Cosmology II: Matter, Halo, and Lensing Statistics in the vector Galileon
• The Cosmic Mach Number as an Environment Measure for the Underlying Dark Matter Density Field
• Cosmological Simulations of Number Counts
• Estimating power spectrum of discrete cosmic momentum field with fast Fourier transform
• The Effect of Corner Modes in the Initial Conditions of Cosmological Simulations
• An Optimally Weighted Estimator of the Linear Power Spectrum Disentangling the Growth of Density Perturbations Across Galaxy Surveys
• The large-scale monopole of the power spectrum in a Euclid-like survey: wide-angle effects, lensing, and the `finger of the observer'
• The effects on structure of a momentum coupling between dark matter and quintessence
• A neutrino model fit to the CMB power spectrum
• Testing growth rate dependence in cosmological perturbation theory using scale-free models