REVIEW 8 cited by
PolyBin3D: A Suite of Optimal and Efficient Power Spectrum and Bispectrum Estimators for Large-Scale Structure
Not yet reviewed by Pith; the record is open.
This paper has not been read by Pith yet. Machine review is queued; the pith claim, tier, and objections will appear here once it completes.
SPECIMEN: schema-true, not a live event
T0 review · schema-true
One-sentence machine reading of the paper's core claim.
pith:XXXXXXXX · record.json · timestamp
abstract
By measuring, modeling and interpreting cosmological datasets, one can place strong constraints on models of the Universe. Central to this effort are summary statistics such as power spectra and bispectra, which condense the high-dimensional data into low-dimensional representations. In this work, we introduce a modern set of estimators for computing such statistics from three-dimensional clustering data, and provide a flexible Python/Cython implementation; PolyBin3D. Working in a maximum-likelihood formalism, we derive general estimators for the two- and three-point functions, which yield unbiased spectra regardless of the survey mask and weighting scheme. These can be directly compared to theory without the need for mask-convolution. Furthermore, we present a numerical scheme for computing the optimal (minimum-variance) estimators for a given survey, which is shown to reduce error-bars on large-scales. Our Python package includes both general "unwindowed" estimators and their idealized equivalents (appropriate for simulations), each of which are efficiently implemented using fast Fourier transforms and Monte Carlo summation tricks, and additionally supports GPU acceleration using JAX. These are extensively validated in this work, with Monte Carlo convergence (relevant for masked data) achieved using only a small number of iterations (typically $<10$ for bispectra). This will allow for fast and unified measurement of two- and three-point functions from current and upcoming survey data.
Forward citations
Cited by 8 Pith papers
-
Control variates from Eulerian and Lagrangian perturbation theory: Application to the bispectrum
A shifted, Zeldovich-resummed control variate with tree-level bispectrum reduces N-body matter-bispectrum variance by up to 10^4 at low k, enabling sub-2% precision from a single 1 (Gpc/h)^3 box.
-
Searching for Inflationary Physics with the CMB Trispectrum: 2. Code & Validation
PolySpec is a public, validated suite of near-optimal quartic estimators that lets the CMB four-point function probe many inflationary models beyond local non-Gaussianity.
-
Searching for Inflationary Physics with the CMB Trispectrum: 1. Primordial Theory & Optimal Estimators
A set of quasi-optimal CMB trispectrum estimators is derived for local, EFT, direction-dependent, spinning-particle, point-source, and lensing templates, enabling first collider searches.
-
Hermes - Towards an Optimal High-Performance Algorithm for Cosmic Statistics of Large Data Sets
Hermes/PyHermes reconstructs catalogues in a scaling-function basis and unifies CIC, 2PCF, 3PCF, marked, and operator-based cosmic statistics as reusable window operations with FFT/MPI/GPU scaling.
-
One-Loop Galaxy Bispectrum: Consistent Theory, Efficient Analysis with COBRA, and Implications for Cosmological Parameters
Using the COBRA basis, the one-loop EFT galaxy bispectrum runs in about one second per cosmology and recovers the input cosmology of PTChallenge simulations up to k_max = 0.15 h/Mpc (monopole) and 0.12 h/Mpc (higher m...
-
Searching for Inflationary Physics with the CMB Trispectrum: 3. Constraints from Planck
A comprehensive Planck PR4 trispectrum analysis finds no primordial non-Gaussianity across 33 templates and sets leading constraints, including tau_NL loc < 1500 at 95% CL.
-
An optimal quadratic estimator for window-free cosmic shear power spectra
A window-free quadratic estimator that is statistically optimal when the true covariance is known improves cosmic shear E-mode precision by 5–15% at ℓ≲500 and suppresses E→B leakage.
-
Inflation, Open Universes, and Dark Energy
A slightly open universe with Ω_k ≈ 3×10^{-3} pulls n_s down to ~0.969, reconciling plateau inflation models with combined CMB and DESI data.
Discussion (0). Continue with ORCID to comment.