pith. sign in

arxiv: 2605.30246 · v1 · pith:LYFQAG4Enew · submitted 2026-05-28 · 🌌 astro-ph.CO

Efficient computation of the galaxy angular bispectrum in redshift space

Zucheng Gao , Zvonimir Vlah , Anthony Challinor This is my paper

Pith reviewed 2026-06-29 05:40 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords angular bispectrumredshift spacelarge-scale structureflat-sky approximationunequal-time correlatorstree-level perturbation theoryangular power spectrum
0
0 comments X

The pith

Galaxy angular bispectrum reduces to sum of angular power spectra products

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

The paper derives the tree-level angular bispectrum in redshift space and shows it equals an exact sum of products of two angular power spectra. This identity lets the flat-sky approximation be applied, which cuts computation time sharply while keeping essential line-of-sight modes. Validation finds sub-percent agreement with full-sky results across scales and configurations when narrow Gaussian windows are used. The efficiency gain matters because direct bispectrum calculations have remained too slow for routine use in modeling large-scale structure surveys.

Core claim

The tree-level unequal-time angular bispectrum in real and redshift space can be expressed as a sum of products of two angular power spectra. This full-sky result supports the flat-sky approximation that preserves key line-of-sight mode information while discarding extraneous contributions. The method reaches sub-percent agreement with direct full-sky integration for multipoles above 5 at redshift 1 in equilateral, squeezed, and folded configurations under narrow Gaussian radial windows of width 0.01.

What carries the argument

Decomposition of the bispectrum into a sum of products of two angular power spectra

If this is right

  • The flat-sky method matches full-sky results to sub-percent level for ell greater than or equal to 5 in multiple triangle shapes.
  • The reduction works for both equal-time and unequal-time cases.
  • On small scales the results agree with the Limber approximation where applicable.
  • A public Python implementation is supplied for reuse.

Where Pith is reading between the lines

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

  • Wider windows typical of photometric surveys may need separate accuracy checks.
  • The same power-spectrum-product structure could be tested for higher-order angular statistics.
  • Embedding the code in survey pipelines would quantify real-world speed gains.

Load-bearing premise

Validation holds only for narrow Gaussian radial window functions with redshift width 0.01.

What would settle it

Running the full-sky versus flat-sky comparison with wider windows such as sigma_z = 0.1 and checking whether sub-percent agreement survives for low multipoles.

read the original abstract

Efficient computation of the angular bispectrum is an essential part of modelling large-scale structure observations, but it still remains an extremely challenging task. In this work, we compute the tree-level, unequal-time angular bispectrum in both real and redshift space. By deriving full-sky results, we show that the bispectrum can be expressed as a sum of products of two angular power spectra, enabling the use of our recently developed flat-sky approximation to enhance computational efficiency significantly. This flat-sky formalism preserves key line-of-sight mode information while discarding extraneous full-sky contributions. We validate our approach by comparing it with direct full-sky integration, finding excellent agreement across a wide range of scales and redshifts for all bispectrum configurations. At redshift $z = 1$, we achieve sub-percent agreement (for multipoles $\ell \gtrsim 5$) between full-sky and flat-sky results for equilateral, squeezed, and folded configurations, using narrow Gaussian radial window functions ($\sigma_z = 0.01$) in both equal-time and unequal-time scenarios. On small scales, where direct full-sky integration becomes computationally prohibitive, our results align with the Limber approximation (where applicable), confirming the robustness and accuracy of our implementation. To facilitate future studies, we provide a \texttt{Python} implementation of our results, which is publicly available on \texttt{GitHub}.

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

1 major / 0 minor

Summary. The paper derives the tree-level unequal-time angular bispectrum (real and redshift space) and shows it factors as a sum of products of two angular power spectra. This factoring enables application of a previously developed flat-sky approximation for efficient computation while preserving key line-of-sight information. Results are validated by direct comparison to full-sky integration, achieving sub-percent agreement for ℓ ≳ 5 at z=1 across equilateral, squeezed and folded configurations when using narrow Gaussian radial windows (σ_z=0.01); the implementation is released as open-source Python code.

Significance. If the central derivation and validation hold, the work supplies a practical route to computing angular bispectra at tree level that is substantially faster than direct full-sky integration on small scales. The public code release is a clear strength that supports reproducibility. The approach is particularly relevant for modeling upcoming spectroscopic and photometric surveys where bispectrum measurements are computationally demanding.

major comments (1)
  1. [Abstract and validation results] Abstract and validation section: sub-percent agreement is reported only for narrow Gaussian windows (σ_z=0.01). No tests are presented for broader windows (σ_z ≳ 0.05) typical of photometric surveys, even though the text states that the flat-sky step 'preserves key line-of-sight mode information' only under the narrow-window regime. Because the efficiency gain rests on this approximation, the lack of broader validation is load-bearing for the claimed generality of the method.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive assessment and for identifying the need to clarify the regime of applicability. We address the single major comment below.

read point-by-point responses

  1. Referee: [Abstract and validation results] Abstract and validation section: sub-percent agreement is reported only for narrow Gaussian windows (σ_z=0.01). No tests are presented for broader windows (σ_z ≳ 0.05) typical of photometric surveys, even though the text states that the flat-sky step 'preserves key line-of-sight mode information' only under the narrow-window regime. Because the efficiency gain rests on this approximation, the lack of broader validation is load-bearing for the claimed generality of the method.

    Authors: We agree that the validation and efficiency claims are demonstrated specifically for narrow radial windows (σ_z = 0.01), which are representative of spectroscopic survey bins. The algebraic factoring of the tree-level angular bispectrum into products of angular power spectra is derived without assuming narrow windows and holds for arbitrary radial selection functions. However, the subsequent flat-sky approximation for those power spectra (and the associated preservation of line-of-sight information) is explicitly presented in the manuscript as valid in the narrow-window regime. For broader windows typical of photometric surveys the approximation error grows and direct integration or alternative methods remain necessary. We will revise the abstract and the opening of the validation section to state this scope more explicitly and to note that extension to photometric windows is left for future work. This change removes any implication of unrestricted generality while preserving the paper’s core result. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives the tree-level full-sky angular bispectrum (real and redshift space) and shows via explicit calculation that it factors as a sum of products of two angular power spectra. This mathematical reduction is presented as an independent result, validated by direct comparison to full-sky numerical integration rather than by construction or fit. The efficiency step invokes a prior flat-sky approximation described as 'our recently developed,' but this is applied after the derivation and the combined results are checked against external full-sky benchmarks across configurations; the self-citation is not load-bearing for the core factoring claim itself. No self-definitional loops, fitted inputs renamed as predictions, uniqueness theorems, or ansatzes smuggled via citation appear in the provided derivation steps.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard tree-level perturbation theory in cosmology and the assumption that the radial selection functions can be treated as narrow Gaussians; no new free parameters are introduced beyond those already present in the angular power spectra, and no new entities are postulated.

axioms (2)
  • domain assumption Tree-level perturbation theory suffices for the bispectrum on the scales considered.
    Invoked when stating the computation is performed at tree level.
  • domain assumption The flat-sky approximation preserves the dominant line-of-sight contributions when radial windows are narrow.
    Central to the efficiency claim and validated only for sigma_z = 0.01.

pith-pipeline@v0.9.1-grok · 5776 in / 1489 out tokens · 19342 ms · 2026-06-29T05:40:13.261102+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 · 15 internal anchors

  1. [1]

    Euclid Definition Study Report

    R. Laureijs, J. Amiaux, S. Arduini, J.L. Auguères, J. Brinchmann, R. Cole et al., Euclid Definition Study Report, arXiv e-prints (2011) arXiv:1110.3193 [ 1110.3193]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  2. [2]

    The DESI Experiment Part I: Science,Targeting, and Survey Design

    DESI Collaboration, A. Aghamousa, J. Aguilar, S. Ahlen, S. Alam, L.E. Allen et al., The DESI Experiment Part I: Science,Targeting, and Survey Design, arXiv e-prints (2016) arXiv:1611.00036 [1611.00036]

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  3. [3]

    LSST Dark Energy Science collaboration, The N5K Challenge: Non-Limber Integration for LSST Cosmology, Open Journal of Astrophysics6 (2023) 1 [ 2212.04291]

    work page arXiv 2023

  4. [4]

    Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report

    D. Spergel, N. Gehrels, C. Baltay, D. Bennett, J. Breckinridge, M. Donahue et al., Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report, arXiv e-prints (2015) arXiv:1503.03757 [ 1503.03757]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  5. [5]

    The power spectrum and bispectrum of SDSS DR11 BOSS galaxies I: bias and gravity

    H. Gil-Marín, J. Noreña, L. Verde, W.J. Percival, C. Wagner, M. Manera et al., The power spectrum and bispectrum of SDSS DR11 BOSS galaxies – I. Bias and gravity, Mon. Not. Roy. Astron. Soc.451 (2015) 539 [ 1407.5668]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  6. [6]

    The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the power spectrum and bispectrum of the DR12 BOSS galaxies

    H. Gil-Marín, W.J. Percival, L. Verde, J.R. Brownstein, C.-H. Chuang, F.-S. Kitaura et al., The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the power spectrum and bispectrum of the DR12 BOSS galaxies, Mon. Not. Roy. Astron. Soc. 465 (2017) 1757 [ 1606.00439]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  7. [7]

    Cabass, M.M

    G. Cabass, M.M. Ivanov, O.H.E. Philcox, M. Simonović and M. Zaldarriaga, Constraints on Single-Field Inflation from the BOSS Galaxy Survey, Phys. Rev. Lett.129 (2022) 021301 [2201.07238]. – 26 – 10 11 10 10 10 9 10 8 B 1 2 3 Real-space Flat-sky Limber 10 11 10 10 10 9 10 8 Redshift-space Partial-Limber 101 102 103 10 0 10 Rel. Err. [%] 101 102 103 10 0 10...

    work page arXiv 2022

  8. [8]

    Ivanov, O.H.E

    M.M. Ivanov, O.H.E. Philcox, G. Cabass, T. Nishimichi, M. Simonović and M. Zaldarriaga, Cosmology with the galaxy bispectrum multipoles: Optimal estimation and application to BOSS data, Phys. Rev. D107 (2023) 083515 [ 2302.04414]

    work page arXiv 2023

  9. [9]

    The Galaxy Power Spectrum and Bispectrum in Redshift Space

    V. Desjacques, D. Jeong and F. Schmidt, The Galaxy Power Spectrum and Bispectrum in Redshift Space, JCAP 12 (2018) 035 [ 1806.04015]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  10. [10]

    Pardede, E

    K. Pardede, E. Di Dio and E. Castorina, Wide-angle effects in the galaxy bispectrum, JCAP 09 (2023) 030 [ 2302.12789]

    work page arXiv 2023

  11. [11]

    Benabou, A

    J.N. Benabou, A. Testa, C. Heinrich, H.S. Grasshorn Gebhardt and O. Doré, Galaxy bispectrum in the spherical Fourier-Bessel basis, Phys. Rev. D109 (2024) 103507 [ 2312.15992]

    work page arXiv 2024

  12. [12]

    Beyond the Linear-Order Relativistic Effect in Galaxy Clustering: Second-Order Gauge-Invariant Formalism

    J. Yoo and M. Zaldarriaga, Beyond the Linear-Order Relativistic Effect in Galaxy Clustering: – 27 – Second-Order Gauge-Invariant Formalism, Phys. Rev. D90 (2014) 023513 [ 1406.4140]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  13. [13]

    The bispectrum of relativistic galaxy number counts

    E. Di Dio, R. Durrer, G. Marozzi and F. Montanari, The bispectrum of relativistic galaxy number counts, JCAP 01 (2016) 016 [ 1510.04202]

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  14. [14]

    Montandon, E

    T. Montandon, E. Di Dio, C. Rampf and J. Adamek, Angular bispectrum of matter number counts in cosmic structures, 2501.05422

    work page arXiv

  15. [15]

    Efficient Evaluation of Cosmological Angular Statistics

    V. Assassi, M. Simonović and M. Zaldarriaga, Efficient evaluation of angular power spectra and bispectra, JCAP 11 (2017) 054 [ 1705.05022]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  16. [16]

    The Full-Sky Angular Bispectrum in Redshift Space

    E. Di Dio, R. Durrer, R. Maartens, F. Montanari and O. Umeh, The Full-Sky Angular Bispectrum in Redshift Space, JCAP 04 (2019) 053 [ 1812.09297]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  17. [17]

    Z. Gao, Z. Vlah and A. Challinor, Flat-sky angular power spectra revisited, Journal of Cosmology and Astroparticle Physics2024 (2024) 003

    2024

  18. [18]

    Z. Gao, A. Raccanelli and Z. Vlah, Asymptotic connection between full- and flat-sky angular correlators, Phys. Rev. D108 (2023) 043503 [ 2306.02993]

    work page arXiv 2023

  19. [19]

    Planck collaboration, Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys.641 (2020) A6 [ 1807.06209]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  20. [20]

    On the Statistics of Biased Tracers in the Effective Field Theory of Large Scale Structures

    R. Angulo, M. Fasiello, L. Senatore and Z. Vlah, On the Statistics of Biased Tracers in the Effective Field Theory of Large Scale Structures, JCAP 09 (2015) 029 [ 1503.08826]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  21. [21]

    Fujita, V

    T. Fujita, V. Mauerhofer, L. Senatore, Z. Vlah and R. Angulo, Very Massive Tracers and Higher Derivative Biases, JCAP 01 (2020) 009 [ 1609.00717]

    work page arXiv 2020

  22. [22]

    Raccanelli and Z

    A. Raccanelli and Z. Vlah, Observed power spectrum and frequency-angular power spectrum, Phys. Rev. D 108 (2023) 043537 [ 2306.00808]

    work page arXiv 2023

  23. [23]

    The Intrinsic Alignment of Galaxies and its Impact on Weak Gravitational Lensing in an Era of Precision Cosmology

    M.A. Troxel and M. Ishak, The Intrinsic Alignment of Galaxies and its Impact on Weak Gravitational Lensing in an Era of Precision Cosmology, Phys. Rept. 558 (2014) 1 [ 1407.6990]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  24. [24]

    Euclid collaboration, Euclid preparation. VII. Forecast validation for Euclid cosmological probes, Astron. Astrophys.642 (2020) A191 [ 1910.09273]

    work page arXiv 2020

  25. [25]

    Linear Redshift Distortions: A Review

    A.J.S. Hamilton, Linear redshift distortions: A Review, in Ringberg Workshop on Large Scale Structure, 8, 1997, DOI [ astro-ph/9708102]

    work page internal anchor Pith review Pith/arXiv arXiv 1997

  26. [26]

    Kaiser, Clustering in real space and in redshift space, MNRAS 227 (1987) 1

    N. Kaiser, Clustering in real space and in redshift space, MNRAS 227 (1987) 1

    1987

  27. [27]

    Uncorrelated Modes of the Nonlinear Power Spectrum

    A.J.S. Hamilton, Uncorrelated modes of the nonlinear power spectrum, Mon. Not. Roy. Astron. Soc. 312 (2000) 257 [ astro-ph/9905191]

    work page internal anchor Pith review Pith/arXiv arXiv 2000

  28. [28]

    Bakx, N.E

    T. Bakx, N.E. Chisari and Z. Vlah, COBRA: Optimal Factorization of Cosmological Observables, 2407.04660. – 28 –

    work page arXiv