Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation

Abstract

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization.

arXiv commentUpdated to reflect published version

References in corpus (82)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Simba: Cosmological Simulations with Black Hole Growth and Feedback
• Large-Scale Galaxy Bias
• Transients from Initial Conditions in Cosmological Simulations
• The BAHAMAS project: Calibrated hydrodynamical simulations for large-scale structure cosmology
• Resumming Cosmological Perturbations via the Lagrangian Picture: One-loop Results in Real Space and in Redshift Space
• Towards a realistic population of simulated galaxy groups and clusters
• Self-regulated growth of supermassive black holes by a dual jet/heating AGN feedback mechanism: methods, tests and implications for cosmological simulations
• Solving Large Scale Structure in Ten Easy Steps with COLA
• The MassiveBlack-II Simulation: The Evolution of Halos and Galaxies to z~0
• The Quijote simulations
• The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum
• Uncorrelated Modes of the Nonlinear Power Spectrum
• nbodykit: an open-source, massively parallel toolkit for large-scale structure
• Nonlinear perturbation theory with halo bias and redshift-space distortions via the Lagrangian picture
• Convolution Lagrangian perturbation theory for biased tracers
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• FastPM: a new scheme for fast simulations of dark matter and halos
• The Lagrangian-space Effective Field Theory of Large Scale Structures
• The BACCO Simulation Project: Exploiting the full power of large-scale structure for cosmology
• Cosmological $N$-body simulations with suppressed variance
• DEMNUni: The clustering of large-scale structures in the presence of massive neutrinos
• A new analysis of galaxy 2-point functions in the BOSS survey, including full-shape information and post-reconstruction BAO
• A Lagrangian effective field theory
• Euclid preparation: IX. EuclidEmulator2 -- Power spectrum emulation with massive neutrinos and self-consistent dark energy perturbations
• Matter power spectrum and the challenge of percent accuracy
• Nonlinear Cosmological Power Spectra in Real and Redshift--Space
• The Aemulus Project I: Numerical Simulations for Precision Cosmology
• Redshift-Space Distortions in Lagrangian Perturbation Theory
• Consistent Modeling of Velocity Statistics and Redshift-Space Distortions in One-Loop Perturbation Theory
• The Mira-Titan Universe: Precision Predictions for Dark Energy Surveys
• The Zeldovich approximation
• Galaxy Bias and Primordial Non-Gaussianity
• Towards optimal extraction of cosmological information from nonlinear data
• Cosmological constraints from the tomographic cross-correlation of DESI Luminous Red Galaxies and Planck CMB lensing
• The Gaussian streaming model and Lagrangian effective field theory
• AI-assisted super-resolution cosmological simulations
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• AbacusHOD: A highly efficient extended multi-tracer HOD framework and its application to BOSS and eBOSS data
• Initial Conditions for Accurate N-Body Simulations of Massive Neutrino Cosmologies
• The cosmology dependence of galaxy clustering and lensing from a hybrid $N$-body-perturbation theory model
• SCORCH I: The Galaxy-Halo Connection in the First Billion Years
• UNIT project: Universe $N$-body simulations for the Investigation of Theoretical models from galaxy surveys
• Fast Large Scale Structure Perturbation Theory using 1D FFTs
• Simulations and symmetries
• Lagrangian perturbation theory at one loop order: successes, failures, and improvements
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• The Mildly Non-Linear Regime of Structure Formation
• The clustering of luminous red galaxies at z $\sim$ 0.7 from eBOSS and BOSS data
• Stringent ${\sigma }_8$ constraints from small-scale galaxy clustering using a hybrid MCMC+emulator framework
• Five-percent measurements of the growth rate from simulation-based modelling of redshift-space clustering in BOSS LOWZ
• Statistical properties of paired fixed fields
• Reconstructing large-scale structure with neutral hydrogen surveys
• Inverted initial conditions: exploring the growth of cosmic structure and voids
• The Aemulus Project V: Cosmological constraint from small-scale clustering of BOSS galaxies
• An $n$-th order Lagrangian Forward Model for Large-Scale Structure
• GridSPT: Grid-based calculation for perturbation theory of large-scale structure
• CARPool: fast, accurate computation of large-scale structure statistics by pairing costly and cheap cosmological simulations
• Lagrangian or Eulerian; Real or Fourier? Not All Approaches to Large-Scale Structure Are Created Equal
• Hefty enhancement of cosmological constraints from the DES Y1 data using a Hybrid Effective Field Theory approach to galaxy bias
• Learning effective physical laws for generating cosmological hydrodynamics with Lagrangian Deep Learning
• On the reach of perturbative methods for dark matter density fields
• Priors on red galaxy stochasticity from hybrid effective field theory
• Priors on Lagrangian bias parameters from galaxy formation modelling
• A gradient based method for modeling baryons and matter in halos of fast simulations
• Cosmological Evidence Modelling: a new simulation-based approach to constrain cosmology on non-linear scales
• sCOLA: The N-body COLA Method Extended to the Spatial Domain
• CARPool Covariance: Fast, unbiased covariance estimation for large-scale structure observables
• Modeling Nearest Neighbor distributions of biased tracers using Hybrid Effective Field Theory
• Shell-crossing in a $\Lambda$CDM Universe
• The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations
• Covariance of the matter power spectrum including the survey window function effect: N-body simulations vs. fifth-order perturbation theory on grid
• Suppressing cosmic variance with paired-and-fixed cosmological simulations: average properties and covariances of dark matter clustering statistics
• N-point Statistics of Large-Scale Structure in the Zel'dovich Approximation
• Statistics of biased tracers in variance-suppressed simulations
• High mass and halo resolution from fast low resolution simulations
• Grid-based calculations of redshift-space matter fluctuations from perturbation theory: UV sensitivity and convergence at the field level
• Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates

Cited by in corpus (14)

• Baryon Acoustic Oscillation Theory and Modelling Systematics for the DESI 2024 results
• Aemulus $\nu$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos
• Precision Redshift-Space Galaxy Power Spectra using Zel'dovich Control Variates
• PNG-UNITsims: Halo clustering response to primordial non-Gaussianities as a function of mass
• Euclid preparation. XLI. Galaxy power spectrum modelling in real space
• HOD-Dependent Systematics in Emission Line Galaxies for the DESI 2024 BAO analysis
• Mitigating the noise of DESI mocks using analytic control variates
• HOD-Dependent Systematics for Luminous Red Galaxies in the DESI 2024 BAO Analysis
• On general-relativistic Lagrangian perturbation theory and its non-perturbative generalization
• CSST Cosmological Emulator III: Hybrid Lagrangian Bias Expansion Emulation of Galaxy Clustering
• ForestFlow: predicting the Lyman-$\alpha$ forest clustering from linear to nonlinear scales
• Variance Reduction via Simultaneous Importance Sampling and Control Variates Techniques Using Vegas
• Simulation budgeting for hybrid effective field theories
• Clustering analysis of medium-band selected high-redshift galaxies