pith. machine review for the scientific record. sign in

arxiv: 2605.04864 · v1 · submitted 2026-05-06 · 🌌 astro-ph.CO

Recognition: unknown

Efficient estimators for power spectrum and bispectrum multipole measurements

Yunchen Xie , Ruiyang Zhao , Gan Gu , Xiaoma Wang , Xiaoyong Mu , Yuting Wang , Gong-Bo Zhao , Florian Beutler , John A. Peacock

Authors on Pith no claims yet

Pith reviewed 2026-05-08 17:14 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords power spectrumbispectrummultipolesFFT estimatorslarge-scale structureshot noisespherical harmonicsgalaxy surveys
0
0 comments X

The pith

Symmetries of spherical-harmonic Fourier transforms halve the computational cost for power and bispectrum multipole estimators.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Large-scale galaxy surveys require fast estimators for directional dependence in clustering statistics. The work demonstrates that conjugation and parity symmetries in the transforms of real fields remove redundant parts of the calculation, cutting the cost exactly in half for both power-spectrum and bispectrum multipoles. A further algebraic reduction expresses high even multipoles using only lower-order polynomials with controlled change to the measured quantity. A compressed bispectrum estimator and an analytic shot-noise formula complete the set of optimizations for practical use on large datasets.

Core claim

The paper claims that a suite of FFT-based estimators for power-spectrum and bispectrum multipoles can be made efficient by using exact conjugation and parity symmetries of spherical-harmonic-weighted Fourier transforms of real fields to eliminate redundant magnetic sub-configurations. This reduces computational cost by a factor of 2. High-order even multipoles are rewritten in terms of lower-order Legendre polynomials for additional savings, and a new TripoSH bispectrum estimator is introduced by compressing along an alternative triangle side, together with an analytic treatment of bispectrum shot noise.

What carries the argument

Exact conjugation and parity symmetries in spherical-harmonic-weighted Fourier transforms of real fields that eliminate redundant magnetic sub-configurations.

If this is right

  • High-order even power-spectrum multipoles can be measured using only low-ℓ fields.
  • The TripoSH estimator substantially reduces FFT scaling for quadrupole bispectrum configurations in the large-k-bin limit.
  • Bispectrum shot noise can be subtracted using an analytic integral over triangle-constrained k-space volumes without additional FFTs.
  • The optimizations support an open-source package for efficient clustering measurements in large galaxy surveys.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • These reductions could enable real-time or repeated measurements on the largest upcoming survey volumes.
  • The approach might generalize to other anisotropic statistics involving higher-order correlations.
  • Direct implementation tests would show the exact speedup factor achieved in practice on survey-like data volumes.
  • The small deviation in modified multipoles may require updating covariance matrices when combining with traditional estimators.

Load-bearing premise

Rewriting high-order even multipoles algebraically in terms of lower-order Legendre polynomials produces only a small and controlled deviation from the traditional Yamamoto definition.

What would settle it

Compare the power-spectrum multipoles obtained from the modified estimator against those from the standard Yamamoto estimator on the same set of mock galaxy catalogs to check whether the differences remain small and controlled as claimed.

Figures

Figures reproduced from arXiv: 2605.04864 by Florian Beutler, Gan Gu, Gong-Bo Zhao, John A. Peacock, Ruiyang Zhao, Xiaoma Wang, Xiaoyong Mu, Yunchen Xie, Yuting Wang.

Figure 1
Figure 1. Figure 1: From left to right: the mean ratio, standard deviation ratio, and SNR ratio of the shot-noise-subtracted power spectrum multipoles Pbℓb to the standard ones Pbℓ, computed from 2048 MultiDark-Patchy mocks of the NGC subsample in the redshift slice z ∈ [0.2, 0.5]. are modeled based on effective field theory (EFT) using the IR-resummed Eulerian flavor of Velocileptors code (S.-F. Chen et al. 2020; M. Maus et … view at source ↗
Figure 2
Figure 2. Figure 2: Parameter fitting results using Pb0, Pb2, incorporating Pb4 and its approximation Pb4b. the compression direction to be k3, yielding a new bispectrum estimator denoted as Bb. Bbℓ1ℓ2L (k1, k2) = Hℓ1ℓ2LNℓ1ℓ2L I X m1m2M view at source ↗
Figure 3
Figure 3. Figure 3: Cosmological parameter contours derived from the Fisher matrices of Sugiyama’s Bb202 and our Bb202. The blue dotted (no fill) and gray dashed (no fill) ellipses correspond to constraints from the two quadrupole estimators, respectively. The green-filled ellipse represents the constraints from their shared monopole. The orange- and purple-filled ellipses show the joint constraints obtained by combining the … view at source ↗
Figure 4
Figure 4. Figure 4: Left: Schematic diagram of rigid body rotation of a triangle in Fourier space from its initial position to an arbitrary position. Right: Schematic diagram of spherical harmonic integration on an arbitrary two-dimensional spherical surface in three-dimensional space. this approach still relies on spherical Bessel functions jℓ(kx), whose real-time evaluation remains relatively complex and offers only a limit… view at source ↗
Figure 5
Figure 5. Figure 5: First five orders of gℓ(t) given by Eq. (54) calculated numerically (dashed lines) and analytically (solid lines). 5.4. Derivation of shot noise in three different bispectrum estimators 5.4.1. Scoccimarro estimator Using the special integrals derived in the previous section, we first calculate the relatively simpler shot noise of the Scoccimarro estimator. Following the notation of Eq. (40), we can readily… view at source ↗
Figure 6
Figure 6. Figure 6: Top left: Shot noise term S000|i=j̸=k scaled by k 2 , evaluated from MD-patchy mocks using Sugiyama’s original method, the FFT method, and the analytical method, respectively. The solid lines represent the mean value, while the colored bands indicate the standard deviation across 200 mocks. Top right and bottom left: Same format as top left, but for S202|i=j̸=k and S220|i=j̸=k, respectively. Bottom right: … view at source ↗
Figure 7
Figure 7. Figure 7: Left: Difference between the two approximate methods and the reference method, divided by the standard deviation of the bispectrum multipoles evaluated using the FFT method. Different colors represent different multipoles, and the two methods are distinguished by distinct markers. Right: Comparison of the full bispectrum standard deviations obtained from the three methods view at source ↗
Figure 8
Figure 8. Figure 8: Residuals of the multipole B220 when the shot noise is computed using the Sugiyama/analytical method, relative to the FFT method, measured from 200 MD-patchy mock catalogs. The dashed lines, along with the corresponding coloured shaded regions, represent the mean residual and the standard deviation for the Sugiyama method, while the solid lines with error bars represent the mean residual and the standard d… view at source ↗
Figure 9
Figure 9. Figure 9: presents the mean and standard deviation of the traditional power spectrum multipoles P0, P2, P4, P6, and P8, measured from 2048 MultiDark-Patchy mock catalogs of the BOSS CMASS NGC sample in the redshift slice z ∈ [0.2, 0.5]. These measurements serve as a baseline for validating the modified estimators Pbℓb introduced in Section 3.1. 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0 500 1000 1500 2000 kP`(k) [ h −2Mpc… view at source ↗
Figure 10
Figure 10. Figure 10: Top left: Mean (solid line) and standard deviation (color band) of the compressed bispectrum monopole, along with quadrupoles derived from two distinct compression directions, for fiducial-cosmology mocks. Top right: Mean shape of the measured bispectrum monopole as a function of varying cosmological parameters. Bottom left: Same as top right, but for the Sugiyama quadrupole. Bottom right: Same as top rig… view at source ↗
Figure 11
Figure 11. Figure 11: A schematic diagram illustrating how to derive the total effective volume for a triangular configuration in Fourier space. J. DERIVATION OF Gℓ(T) Integrating and averaging over the spherical surface SP shown in the right panel of view at source ↗
read the original abstract

Large galaxy surveys demand fast and scalable estimators for anisotropic clustering statistics beyond the monopole. We present a suite of efficient FFT-based estimators for power-spectrum and bispectrum multipoles, built upon exact conjugation and parity symmetries of spherical-harmonic--weighted Fourier transforms of real fields. These symmetries eliminate redundant magnetic sub-configurations, thereby reducing the computational cost by a factor of 2. For the Yamamoto power-spectrum multipoles, we further decrease the cost of high-order even multipoles by algebraically expressing ${L}_{2n}$ in terms of lower-order Legendre polynomials, thereby measuring modified high-order multipoles using only low-$\ell$ fields with a small and controlled deviation from the traditional definition. We introduce a new TripoSH bispectrum estimator obtained by compressing the Scoccimarro bispectrum along an alternative triangle side, which substantially reduces the FFT scaling for commonly used quadrupole configurations in the large-$k$-bin limit. We also derive an analytic treatment of bispectrum shot noise by integrating spherical-harmonic kernels over the triangle-constrained $k$-space volumes, avoiding additional FFTs or costly spherical-Bessel evaluations and enabling fast and accurate shot-noise subtraction. Based on these optimizations, we also introduce CosmoNPC, an open-source Python package for large-scale-structure clustering measurements.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a suite of efficient FFT-based estimators for power-spectrum and bispectrum multipoles in large-scale structure analyses. It exploits exact conjugation and parity symmetries of spherical-harmonic-weighted Fourier transforms of real fields to eliminate redundant magnetic sub-configurations and halve the computational cost. For Yamamoto power-spectrum multipoles, high-order even multipoles L_{2n} are algebraically rewritten in terms of lower-order Legendre polynomials, yielding modified estimators with a claimed small controlled deviation from the standard definition. A new TripoSH bispectrum estimator is introduced by compressing the Scoccimarro estimator along an alternative triangle side, reducing FFT scaling in the large-k-bin limit for quadrupole configurations. An analytic treatment of bispectrum shot noise is derived by integrating spherical-harmonic kernels over triangle-constrained k-space volumes. The work also releases the open-source CosmoNPC Python package implementing these methods.

Significance. If the controlled approximations are shown to be sufficiently accurate, the estimators could enable substantially faster multipole measurements for upcoming surveys such as DESI and Euclid while preserving compatibility with standard cosmological modeling. The open-source implementation and analytic shot-noise formula represent concrete strengths that support reproducibility and practical adoption.

major comments (2)
  1. [§3.2] §3.2, Eq. (14): The algebraic rewriting of L_{2n} in terms of lower-order Legendre polynomials is asserted to produce only a small and controlled deviation from the traditional Yamamoto definition. No analytic bound on the k-dependent difference or explicit numerical comparison (e.g., relative error plots across 0.01 < k < 0.2 h/Mpc) is provided; without this, the practical utility for analyses relying on precise multipole ratios remains conditional.
  2. [§5.3] §5.3, Eq. (28): The analytic bispectrum shot-noise formula is derived by integrating kernels over constrained k-space volumes. Its accuracy for all triangle configurations (especially squeezed and folded limits) must be validated against direct numerical integration or Monte Carlo estimates to confirm absence of systematic residuals at the level required for covariance modeling.
minor comments (2)
  1. [Table 1] Table 1: The timing benchmarks should specify the exact k-bin widths and number of realizations used, to allow direct reproduction of the reported speed-up factors.
  2. [§2.1] §2.1: Notation for the spherical-harmonic-weighted fields (e.g., distinction between Y_{ℓm} and the modified weights) could be made more explicit to avoid confusion with standard Yamamoto conventions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and will revise the manuscript to incorporate the suggested validations, thereby strengthening the presentation of the estimators.

read point-by-point responses

  1. Referee: [§3.2] §3.2, Eq. (14): The algebraic rewriting of L_{2n} in terms of lower-order Legendre polynomials is asserted to produce only a small and controlled deviation from the traditional Yamamoto definition. No analytic bound on the k-dependent difference or explicit numerical comparison (e.g., relative error plots across 0.01 < k < 0.2 h/Mpc) is provided; without this, the practical utility for analyses relying on precise multipole ratios remains conditional.

    Authors: We agree that explicit validation strengthens the claim of a controlled deviation. The algebraic relation follows directly from the recurrence properties of Legendre polynomials and the even parity of the multipoles, ensuring the difference is bounded by the width of the k-bin and the smoothness of the underlying power spectrum. In the revised manuscript we will add both the explicit analytic expression for the k-dependent difference and numerical comparisons (relative error plots) between the modified and standard Yamamoto estimators over 0.01 < k < 0.2 h/Mpc, confirming that the fractional deviation remains well below the statistical precision of upcoming surveys. revision: yes

  2. Referee: [§5.3] §5.3, Eq. (28): The analytic bispectrum shot-noise formula is derived by integrating kernels over constrained k-space volumes. Its accuracy for all triangle configurations (especially squeezed and folded limits) must be validated against direct numerical integration or Monte Carlo estimates to confirm absence of systematic residuals at the level required for covariance modeling.

    Authors: The formula is obtained by exact integration of the spherical-harmonic kernels over the triangle-constrained k-space volumes, so it is analytically precise by construction. Nevertheless, to address possible numerical subtleties in implementation and to demonstrate robustness for covariance applications, we will include in the revised manuscript direct comparisons of the analytic expression against both numerical quadrature and Monte Carlo sampling for representative triangle configurations, with particular emphasis on the squeezed and folded limits. These tests will confirm the absence of systematic residuals at the sub-percent level. revision: yes

Circularity Check

0 steps flagged

No significant circularity; efficiency claims rest on standard symmetries and explicit algebraic approximations.

full rationale

The derivation chain relies on exact conjugation/parity properties of spherical-harmonic-weighted Fourier transforms (standard mathematical identities, not self-defined) and algebraic rewriting of even Legendre polynomials L_{2n} into lower-order terms, which the paper explicitly flags as producing a controlled deviation rather than an exact match to the Yamamoto estimator. No equations reduce claimed predictions or efficiencies to fitted parameters, self-citations, or ansatzes imported from prior author work. The bispectrum shot-noise analytic treatment and TripoSH compression follow from direct integration over k-space volumes without circular reduction. The central efficiency factor-of-2 reduction is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard properties of Fourier transforms and spherical harmonics together with algebraic identities; no new free parameters or postulated entities are introduced.

axioms (2)
  • domain assumption Spherical-harmonic-weighted Fourier transforms of real fields obey exact conjugation and parity symmetries that eliminate redundant magnetic sub-configurations.
    Invoked to achieve the factor-of-two cost reduction for both power spectrum and bispectrum estimators.
  • standard math Higher even multipoles L_{2n} can be expressed in terms of lower-order Legendre polynomials with only small controlled deviation from the standard definition.
    Used to lower the FFT cost for high-order even multipoles in the Yamamoto estimator.

pith-pipeline@v0.9.0 · 5553 in / 1429 out tokens · 98507 ms · 2026-05-08T17:14:24.167102+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 26 canonical work pages

  1. [1]

    Abdul Karimet al.[DESI], Phys

    Abdul Karim, M., et al. 2025, Phys. Rev. D, 112, 083515, doi: 10.1103/tr6y-kpc6

    work page doi:10.1103/tr6y-kpc6 2025

  2. [2]

    , keywords =

    Adame, A. G., et al. 2025, JCAP, 09, 008, doi: 10.1088/1475-7516/2025/09/008

    work page doi:10.1088/1475-7516/2025/09/008 2025

  3. [3]

    , keywords =

    Alam, S., Albareti, F. D., Prieto, C. A., et al. 2015, The Astrophysical Journal Supplement Series, 219, 12, doi: 10.1088/0067-0049/219/1/12

    work page doi:10.1088/0067-0049/219/1/12 2015

  4. [4]

    B., Weber, H

    Arfken, G. B., Weber, H. J., & Harris, F. E. 2011, Mathematical methods for physicists: a comprehensive guide (Academic press)

    2011

  5. [5]

    2021, JCAP, 11, 031, doi: 10.1088/1475-7516/2021/11/031

    Beutler, F., & McDonald, P. 2021, JCAP, 11, 031, doi: 10.1088/1475-7516/2021/11/031

    work page doi:10.1088/1475-7516/2021/11/031 2021

  6. [6]

    2014, Monthly Notices of the Royal Astronomical Society, 443, 1065–1089, doi: 10.1093/mnras/stu1051

    Beutler, F., Saito, S., Seo, H.-J., et al. 2014, Monthly Notices of the Royal Astronomical Society, 443, 1065–1089, doi: 10.1093/mnras/stu1051

    work page doi:10.1093/mnras/stu1051 2014

  7. [7]

    Bianchi, D., Gil-Mar´ ın, H., Ruggeri, R., & Percival, W. J. 2015, Monthly Notices of the Royal Astronomical Society: Letters, 453, L11–L15, doi: 10.1093/mnrasl/slv090

    work page doi:10.1093/mnrasl/slv090 2015

  8. [8]

    2020, Journal of Physics A: Mathematical and Theoretical

    Cahn, R., & Slepian, Z. 2020, Journal of Physics A: Mathematical and Theoretical

    2020

  9. [9]

    Consistent Modeling of Velocity Statistics and Redshift-Space Distortions in One-Loop Perturbation Theory

    Chen, S.-F., Vlah, Z., & White, M. 2020, JCAP, 07, 062, doi: 10.1088/1475-7516/2020/07/062

    work page doi:10.1088/1475-7516/2020/07/062 2020

  10. [10]

    The Astrophysical Journal , author =

    Feldman, H. A., Kaiser, N., & Peacock, J. A. 1994, Astrophys. J., 426, 23, doi: 10.1086/174036 Gil-Mar´ ın, H., Nore˜ na, J., Verde, L., et al. 2015, MNRAS, 451, 539, doi: 10.1093/mnras/stv961 Gil-Mar´ ın, H., Percival, W. J., Verde, L., et al. 2017, MNRAS, 465, 1757, doi: 10.1093/mnras/stw2679 Gil-Mar´ ın, H., Nore˜ na, J., Verde, L., et al. 2015, Monthl...

    work page doi:10.1086/174036 1994

  11. [11]

    Constraining M_ with the bispectrum

    Hahn, C., & Villaescusa-Navarro, F. 2021, Journal of Cosmology and Astroparticle Physics, 2021, 029, doi: 10.1088/1475-7516/2021/04/029

    work page doi:10.1088/1475-7516/2021/04/029 2021

  12. [12]

    2020, Journal of Cosmology and Astroparticle Physics, 2020, 040–040, doi: 10.1088/1475-7516/2020/03/040

    Scoccimarro, R. 2020, Journal of Cosmology and Astroparticle Physics, 2020, 040–040, doi: 10.1088/1475-7516/2020/03/040

    work page doi:10.1088/1475-7516/2020/03/040 2020

  13. [13]

    2017, Journal of Cosmology and Astroparticle Physics, 2017, 002, doi: 10.1088/1475-7516/2017/07/002

    Hand, N., Li, Y., Slepian, Z., & Seljak, U. 2017, Journal of Cosmology and Astroparticle Physics, 2017, 002, doi: 10.1088/1475-7516/2017/07/002

    work page doi:10.1088/1475-7516/2017/07/002 2017

  14. [14]

    Astronomy & Astrophysics , author =

    Hartlap, J., Simon, P., & Schneider, P. 2006, Astronomy &amp; Astrophysics, 464, 399–404, doi: 10.1051/0004-6361:20066170

    work page doi:10.1051/0004-6361:20066170 2006

  15. [15]

    2023, Measuringf NL with the SPHEREx Multi-tracer Redshift Space Bispectrum, https://arxiv.org/abs/2311.13082

    Heinrich, C., Dore, O., & Krause, E. 2023, Measuringf NL with the SPHEREx Multi-tracer Redshift Space Bispectrum, https://arxiv.org/abs/2311.13082

    work page arXiv 2023

  16. [16]

    Jing, Y. P. 2005, The Astrophysical Journal, 620, 559, doi: 10.1086/427087

    work page doi:10.1086/427087 2005

  17. [17]

    , year = 2016, volume =

    Kitaura, F.-S., Rodr´ ıguez-Torres, S., Chuang, C.-H., et al. 2016, Monthly Notices of the Royal Astronomical Society, 456, 4156, doi: 10.1093/mnras/stv2826

    work page doi:10.1093/mnras/stv2826 2016

  18. [18]

    An analysis of parameter compression and Full-Modeling techniques with Velocileptors for DESI 2024 and beyond

    Maus, M., et al. 2025, JCAP, 01, 138, doi: 10.1088/1475-7516/2025/01/138

    work page doi:10.1088/1475-7516/2025/01/138 2025

  19. [19]

    Philcox, O. H. E., & Slepian, Z. 2022, Proceedings of the National Academy of Science, 119, e2111366119, doi: 10.1073/pnas.2111366119 Rodr´ ıguez-Torres, S. A., Chuang, C.-H., Prada, F., et al. 2016, Monthly Notices of the Royal Astronomical Society, 460, 1173, doi: 10.1093/mnras/stw1014 33

    work page doi:10.1073/pnas.2111366119 2022

  20. [20]

    2000, ApJ, 544, 597, doi: 10.1086/317248

    Scoccimarro, R. 2000, The Astrophysical Journal, 544, 597, doi: 10.1086/317248

    work page doi:10.1086/317248 2000

  21. [21]

    2015, PhRvD, 92, 083532, doi: 10.1103/PhysRevD.92.083532 Seethapuram Sridhar, S., Williams, W., Price, D., Breen, s., & Ball, L

    Scoccimarro, R. 2015, Phys. Rev. D, 92, 083532, doi: 10.1103/PhysRevD.92.083532

    work page doi:10.1103/physrevd.92.083532 2015

  22. [22]

    N., et al

    Scoccimarro, R., Colombi, S., Fry, J. N., et al. 1998, ApJ, 496, 586, doi: 10.1086/305399

    work page doi:10.1086/305399 1998

  23. [23]

    S., Saito, S., Beutler, F., & Seo, H.-J

    Sugiyama, N. S., Saito, S., Beutler, F., & Seo, H.-J. 2019, MNRAS, 484, 364, doi: 10.1093/mnras/sty3249

    work page doi:10.1093/mnras/sty3249 2019

  24. [24]

    S., Shiraishi, M., & Okumura, T

    Sugiyama, N. S., Shiraishi, M., & Okumura, T. 2018, MNRAS, 473, 2737, doi: 10.1093/mnras/stx2333

    work page doi:10.1093/mnras/stx2333 2018

  25. [25]

    A., Moskalev, A

    Varshalovich, D. A., Moskalev, A. N., & Khersonskii, V. K. 1988, Quantum Theory of Angular Momentum: Irreducible Tensors, Spherical Harmonics, Vector Coupling Coefficients, 3nj Symbols (World Scientific Publishing Company), doi: 10.1142/0270

    work page doi:10.1142/0270 1988

  26. [26]

    2020, The Quijote simulations, Astrophys

    Villaescusa-Navarro, F., Hahn, C., Massara, E., et al. 2020, The Astrophysical Journal Supplement Series, 250, 2, doi: 10.3847/1538-4365/ab9d82

    work page doi:10.3847/1538-4365/ab9d82 2020

  27. [27]

    2024, Communications Physics, 7, doi: 10.1038/s42005-024-01624-7

    Wang, Y., Zhao, G.-B., Koyama, K., et al. 2024, Communications Physics, 7, doi: 10.1038/s42005-024-01624-7

    work page doi:10.1038/s42005-024-01624-7 2024

  28. [28]

    A., & Nishioka, H

    Yamamoto, K., Nakamichi, M., Kamino, A., Bassett, B. A., & Nishioka, H. 2006, Publications of the Astronomical Society of Japan, 58, 93, doi: 10.1093/pasj/58.1.93

    work page doi:10.1093/pasj/58.1.93 2006