Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation
arXiv:2205.15327 · doi:10.1088/1475-7516/2022/09/059
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.
Updated to reflect published version
References in corpus (16)
- Array Programming with NumPy
- The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
- The DESI Experiment Part I: Science,Targeting, and Survey Design
- Simba: Cosmological Simulations with Black Hole Growth and Feedback
- Transients from Initial Conditions in Cosmological Simulations
- Resumming Cosmological Perturbations via the Lagrangian Picture: One-loop Results in Real Space and in Redshift Space
- AbacusSummit: A Massive Set of High-Accuracy, High-Resolution -Body Simulations
- Redshift-Space Distortions in Lagrangian Perturbation Theory
- Consistent Modeling of Velocity Statistics and Redshift-Space Distortions in One-Loop Perturbation Theory
- The Gaussian streaming model and Lagrangian effective field theory
- Simulations and symmetries
- Lagrangian perturbation theory at one loop order: successes, failures, and improvements
- The Aemulus Project V: Cosmological constraint from small-scale clustering of BOSS galaxies
- sCOLA: The N-body COLA Method Extended to the Spatial Domain
- Statistics of biased tracers in variance-suppressed simulations
- 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 : 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- 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