A Faster Fourier Transform? Computing Small-Scale Power Spectra and Bispectra for Cosmological Simulations in $\mathcal{O}(N^2)$ Time

Abstract

We present $\mathcal{O}(N^2)$ estimators for the small-scale power spectrum and bispectrum in cosmological simulations. In combination with traditional methods, these allow spectra to be efficiently computed across a vast range of scales, requiring orders of magnitude less computation time than Fast Fourier Transform based approaches alone. These methods are applicable to any tracer; simulation particles, halos or galaxies, and take advantage of the simple geometry of the box and periodicity to remove almost all dependence on large random particle catalogs. By working in configuration-space, both power spectra and bispectra can be computed via a weighted sum of particle pairs up to some radius, which can be reduced at larger $k$, leading to algorithms with decreasing complexity on small scales. These do not suffer from aliasing or shot-noise, allowing spectra to be computed to arbitrarily large wavenumbers. The estimators are rigorously derived and tested against simulations, and their covariances discussed. The accompanying code, HIPSTER, has been publicly released, incorporating these algorithms. Such estimators will be of great use in the analysis of large sets of high-resolution simulations.

arXiv comment31 pages, 6 figures, accepted by MNRAS. Code available at https://HIPSTER.readthedocs.io

References in corpus (73)

• The cosmological simulation code GADGET-2
• Detection of the Baryon Acoustic Peak in the Large-Scale Correlation Function of SDSS Luminous Red Galaxies
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples
• Cosmological Constraints from the SDSS Luminous Red Galaxies
• Astrophysics Source Code Library
• Power Spectrum Analysis of Three-Dimensional Redshift Surveys
• The 2dF Galaxy Redshift Survey: The bias of galaxies and the density of the Universe
• Cosmology from cosmic shear power spectra with Subaru Hyper Suprime-Cam first-year data
• The BAHAMAS project: Calibrated hydrodynamical simulations for large-scale structure cosmology
• The WiggleZ Dark Energy Survey: the growth rate of cosmic structure since redshift z=0.9
• The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues
• The Cosmological Analysis of the SDSS/BOSS data from the Effective Field Theory of Large-Scale Structure
• Cosmological Parameters from the BOSS Galaxy Power Spectrum
• The Quijote simulations
• The Bispectrum as a Signature of Gravitational Instability in Redshift-Space
• Power Spectrum Super-Sample Covariance
• Cosmological constraints from cosmic shear two-point correlation functions with HSC survey first-year data
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in Fourier-space
• Cosmology and the Bispectrum
• nbodykit: an open-source, massively parallel toolkit for large-scale structure
• Power Spectrum Correlations Induced by Non-Linear Clustering
• Accurate Estimators of Correlation Functions in Fourier Space
• The power spectrum and bispectrum of SDSS DR11 BOSS galaxies I: bias and gravity
• Combining Full-Shape and BAO Analyses of Galaxy Power Spectra: A 1.6% CMB-independent constraint on H0
• Fast Estimators for Redshift-Space Clustering
• A Measurement of the Quadrupole Power Spectrum in the Clustering of the 2dF QSO Survey
• Super-Sample Covariance in Simulations
• The Bispectrum of IRAS Galaxies
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: BAO measurement from the LOS-dependent power spectrum of DR12 BOSS galaxies
• Cosmological Parameters and Neutrino Masses from the Final Planck and Full-Shape BOSS Data
• The Weak Lensing Signal and the Clustering of BOSS Galaxies II: Astrophysical and Cosmological Constraints
• Measuring line-of-sight dependent Fourier-space clustering using FFTs
• The Aemulus Project I: Numerical Simulations for Precision Cosmology
• Detection of Baryon Acoustic Oscillation Features in the Large-Scale 3-Point Correlation Function of SDSS BOSS DR12 CMASS Galaxies
• Galaxy Power Spectrum Multipoles Covariance in Perturbation Theory
• Computing the Three-Point Correlation Function of Galaxies in $\mathcal{O}(N^2)$ Time
• The Abacus Cosmos: A Suite of Cosmological N-body Simulations
• Assessment of the Information Content of the Power Spectrum and Bispectrum
• An optimal FFT-based anisotropic power spectrum estimator
• A complete FFT-based decomposition formalism for the redshift-space bispectrum
• The Weak Lensing Signal and the Clustering of BOSS Galaxies I: Measurements
• The power spectrum and bispectrum of SDSS DR11 BOSS galaxies II: cosmological interpretation
• A Detection of the Baryon Acoustic Oscillation Features in the SDSS BOSS DR12 Galaxy Bispectrum
• Accelerating the 2-point and 3-point galaxy correlation functions using Fourier transforms
• Perturbation theory approach to predict the covariance matrices of the galaxy power spectrum and bispectrum in redshift space
• Rapid Separable Analysis of Higher Order Correlators in Large Scale Structure
• A fast estimator for the bispectrum and beyond - A practical method for measuring non-Gaussianity in 21-cm maps
• Modeling the large-scale redshift-space 3-point correlation function of galaxies
• A Practical Computational Method for the Anisotropic Redshift-Space 3-Point Correlation Function
• Efficient Evaluation of Cosmological Angular Statistics
• Testing general relativity with the multipole spectra of the SDSS luminous red galaxies
• Galaxy power-spectrum responses and redshift-space super-sample effect
• Three-Point Statistics from a New Perspective
• The matter power spectrum in redshift space using effective field theory
• Large Covariance Matrices: Smooth Models from the 2-Point Correlation Function
• Fast Estimation of Gravitational and Primordial Bispectra in Large Scale Structures
• RascalC: A Jackknife Approach to Estimating Single and Multi-Tracer Galaxy Covariance Matrices
• Estimating Covariance Matrices for Two- and Three-Point Correlation Function Moments in Arbitrary Survey Geometries
• Computing the Small-Scale Galaxy Power Spectrum and Bispectrum in Configuration-Space
• Disconnected Covariance of 2-point Functions in Large-Scale Structure
• Large Covariance Matrices: Accurate Models Without Mocks
• Developing the 3-Point Correlation Function For the Turbulent Interstellar Medium
• Impact of bias and redshift-space modelling for the halo power spectrum: Testing the effective field theory of large-scale structure
• Tau flavored dark matter and its impact on tau Yukawa coupling
• Measurement of Redshift Space Power Spectrum for BOSS galaxies and the Growth Rate at redshift 0.57
• Simulating the Anisotropic Clustering of Luminous Red Galaxies with Subhalos: A Direct Confrontation with Observation and Cosmological Implications
• The Scaling of the Redshift Power Spectrum: Observations from the Las Campanas Redshift Survey
• A Generalisation For The Infinite Integral Over Three Spherical Bessel Functions
• Rotation method for accelerating multiple-spherical Bessel function integrals against a numerical source function
• Computing 3 point correlation function randoms counts without the randoms catalogue
• Indefinite Integrals of Spherical Bessel Functions
• Advancing the matter bispectrum estimation of large-scale structure: a comparison of dark matter codes
• The Plane Wave Expansion, Infinite Integrals and Identities involving Spherical Bessel Functions

Cited by in corpus (15)

• Large-scale dark matter simulations
• Precision analysis of the redshift-space galaxy bispectrum
• The GIGANTES dataset: precision cosmology from voids in the machine learning era
• The CAMELS project: public data release
• Cosmology Without Window Functions: Cubic Estimators for the Galaxy Bispectrum
• ENCORE: An $\mathcal{O}(N_{\mathrm{g}}^2)$ Estimator for Galaxy $N$-Point Correlation Functions
• Cosmology with cosmic web environments I. Real-space power spectra
• The 3D Lyman-$\alpha$ Forest Power Spectrum from eBOSS DR16
• Modeling the Marked Spectrum of Matter and Biased Tracers in Real- and Redshift-Space
• Beyond Yamamoto: Anisotropic Power Spectra and Correlation Functions with Pairwise Lines-of-Sight
• On the Fast Random Sampling and Other Properties of the Three Point Correlation Function in Galaxy Surveys
• An Exact Integral-to-Sum Relation for Products of Bessel Functions
• Predictive uncertainty on improved astrophysics recovery from multifield cosmology
• The Compressed 3D Lyman-Alpha Forest Bispectrum
• Cosmology with the Redshift-Space Galaxy Bispectrum Monopole at One-Loop Order