Fast generation of ensembles of cosmological N-body simulations via mode-resampling

Abstract

We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the large-scale (Gaussian-distributed) Fourier modes in a periodic N-body simulation box in a manner that properly accounts for the nonlinear mode-coupling between large and small scales. We find that our method for adding new large-scale mode realizations recovers the nonlinear power spectrum to sub-percent accuracy on scales larger than about half the Nyquist frequency of the simulation box. Using 20 N-body simulations, we obtain a power spectrum covariance matrix estimate that matches the estimator in Takahashi et al. (2009) (from 5000 simulations) with < 20% errors in all matrix elements. Comparing the rates of convergence, we determine that our algorithm requires ~ 8 times fewer simulations to achieve a given error tolerance in estimates of the power spectrum covariance matrix. The degree of success of our algorithm indicates that we understand the main physical processes that give rise to the correlations in the matter power spectrum. Namely, the large-scale Fourier modes modulate both the degree of structure growth through the variation in the effective local matter density and also the spatial frequency of small-scale perturbations through large-scale displacements. We expect our algorithm to be useful for noise modeling when constraining cosmological parameters from weak lensing (cosmic shear) and galaxy surveys, rescaling summary statistics of N-body simulations for new cosmological parameter values, and any applications where the influence of Fourier modes larger than the simulation size must be accounted for.

13 pages, 11 figures. Matches version accepted for publication in ApJ. New panel in fig. 11 showing the covariance convergence rate with the number of mode resamplings

References in corpus (10)

• Cosmology and the Bispectrum
• The Coyote Universe III: Simulation Suite and Precision Emulator for the Nonlinear Matter Power Spectrum
• Simulations of Wide-Field Weak Lensing Surveys I: Basic Statistics and Non-Gaussian Effects
• Simulations of Baryon Acoustic Oscillations II: Covariance matrix of the matter power spectrum
• Modeling the angular correlation function and its full covariance in Photometric Galaxy Surveys
• Cosmic variance of weak lensing surveys in the non-Gaussian regime
• Shrinkage Estimation of the Power Spectrum Covariance Matrix
• Information content in the halo-model dark-matter power spectrum
• Simulations and cosmological inference: A statistical model for power spectra means and covariances
• Intelligent Design: On the Emulation of Cosmological Simulations

Cited by in corpus (34)

• Weak lensing for precision cosmology
• The Effect of Covariance Estimator Error on Cosmological Parameter Constraints
• Putting the Precision in Precision Cosmology: How accurate should your data covariance matrix be?
• Super-Sample Covariance in Simulations
• Modelling Baryon Acoustic Oscillations with Perturbation Theory and Stochastic Halo Biasing
• Galaxy Bias and non-Linear Structure Formation in General Relativity
• Gaussian covariance matrices for anisotropic galaxy clustering measurements
• Cosmological Structure Formation with Augmented Lagrangian Perturbation Theory
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: theoretical systematics and Baryon Acoustic Oscillations in the galaxy correlation function
• Comparing approximate methods for mock catalogues and covariance matrices I: correlation function
• Super-Sample Signal
• Constraining the halo bispectrum in real and redshift space from perturbation theory and non-linear stochastic bias
• New approach for precise computation of Lyman-alpha forest power spectrum with hydrodynamical simulations
• Perturbation theory approach to predict the covariance matrices of the galaxy power spectrum and bispectrum in redshift space
• Galaxy power-spectrum responses and redshift-space super-sample effect
• Dark Energy Survey Year 1 Results: galaxy mock catalogues for BAO
• Non-linear shrinkage estimation of large-scale structure covariance
• Estimating the power spectrum covariance matrix with fewer mock samples
• Galaxy 2-Point Covariance Matrix Estimation for Next Generation Surveys
• On estimating cosmology-dependent covariance matrices
• Hierarchical probabilistic inference of cosmic shear
• Non-Gaussianity in the Weak Lensing Correlation Function Likelihood -- Implications for Cosmological Parameter Biases
• Response approach to the squeezed-limit bispectrum: application to the correlation of quasar and Lyman-$\alpha$ forest power spectrum
• Extending the halo mass resolution of $N$-body simulations
• The MIP Ensemble Simulation: Local Ensemble Statistics in the Cosmic Web
• Cosmic flows and the expansion of the Local Universe from nonlinear phase-space reconstructions
• Robustness of the covariance matrix for galaxy clustering measurements
• Optimal Weights For Measuring Redshift Space Distortions in Multi-tracer Galaxy Catalogues
• On fast Generation of Cosmological Random Fields
• Compactified Cosmological Simulations of the Infinite Universe
• Keeping It Real: Revisiting a Real-Space Approach to Running Ensembles of Cosmological N-body Simulations
• Cosmological Perturbation Theory as a Tool for Estimating Box-Scale Effects in N-body Simulations
• The Parameter-Level Performance of Covariance Matrix Conditioning in Cosmic Microwave Background Data Analyses
• Improving cosmological covariance matrices with machine learning