REVIEW 2 major objections 3 minor 1 cited by
On the equivalence between galaxy angular correlation function and power spectrum in constraining primordial non-Gaussianity
T0 review · 2 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Scale cuts break the equivalence of galaxy angular power spectrum and correlation function, making small-angle clustering a viable fNL probe.
desk verdict The all-sky scale-cut equivalence argument is clean and worth publishing, but the abstract's partial-sky claim is an extrapolation from fsky=1 forecasts, not a computed result. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing mechanism is the integral constraint: in an all-sky survey, the observed two-point statistics are defined relative to the mean galaxy density estimated from the data itself, which removes the monopole $\ell=0$ (equivalently subtracts $C_{\ell=0}/4\pi$ from the correlation function). This subtraction is a single local operation in harmonic space but a non-local, angle-dependent correction in real space, so the PNG-induced low-multipole power is redistributed across all angular separations of $w(\theta)$. Combined with the scale-dependent bias $\Delta b(k)\propto f_{\rm NL}/k^2$ of Eq. (4), this makes small-angle and partial-sky correlation measurements carry $f_{\rm NL}$ information that harmonic-space scale cuts would discard.
What would settle it
Run a large-volume N-body simulation seeded with a known nonzero $f_{\rm NL}$ (or with halos assigned the scale-dependent bias of Eq. 4), populate it with a high-redshift galaxy-like sample, measure $w(\theta)$ on a partial-sky footprint, and compare the recovered $f_{\rm NL}$ constraint with the paper's Fisher forecast; if the small-angle bins below the nonlinear scale do not retain the predicted $f_{\rm NL}$ signal, or the realized covariance is larger, the central claim would be refuted.
Extended reading notes
Core claim
The central claim is that the observed correlation function and observed power spectrum are equivalent, $\widetilde w(\theta) = \sum_{\ell\ge0} \frac{2\ell+1}{4\pi} \widetilde C_\ell P_\ell(\cos\theta)$ (Eq. 13), once the integral constraint is imposed by dropping the monopole. Yet this equivalence is unstable in practice: the local-type PNG signal in $C_\ell$ is a low-multipole phenomenon, while the same signal enters $w(\theta)$ as an additive, angle-dependent correction from the integral constraint, so it appears at all separations, even below the nonlinear scale. Consequently, when scale cuts remove the largest scales, $C_\ell$ loses most of its $f_{\rm NL}$ sensitivity while $w(\theta)$ keeps it; the paper forecasts that $w(\theta)$ measured up to about 10 degrees for a high-redshift sample ($4<z<7$) yields $\sigma_{f_{\rm NL}}\sim O(1)$ and that combining all bins can reach $\sigma_{f_{\rm NL}}\sim 1$. The paper further notes that the same mechanism should apply to three-dimensional clustering, since any finite survey must estimate its mean density from itself.
Load-bearing premise
The forecast assumes the linear scale-dependent bias model of Eqs. (3)-(5) stays valid down to the small angular scales included, and that the Gaussian, $f_{\rm sky}=1$ covariance is an accurate error model; if nonlinear galaxy bias or unmodeled small-scale systematics contaminate those scales, the projected precision on $f_{\rm NL}$ would not be reached.
Editorial extensions
If this is right
- A high-redshift galaxy sample ($4<z<7$) can reach $\sigma_{f_{\rm NL}}\sim O(1)$ from $w(\theta)$ measured only up to about 10 degrees, nearly matching the full-sky harmonic-space result.
- Combining all six photometric redshift bins, small-angle $w(\theta)$ gives $\sigma_{f_{\rm NL}}\sim 1$, comparable to the $C_\ell$ analysis, when the nonlinear scale is included.
- For partial-sky surveys, real-space analysis is more efficient than harmonic-space analysis, because the PNG signal is not concentrated on the largest scales that are hardest to recover from a cut sky.
- The equivalence between $C_\ell$ and $w(\theta)$ Fisher forecasts is recovered only when $w(\theta)$ reaches scales below the nonlinear scale; cutting $w(\theta)$ at the nonlinear scale loses most of the $f_{\rm NL}$ information.
- The same non-local integral-constraint redistribution should appear in 3D correlation functions of spectroscopic surveys, suggesting small-scale 3D clustering also carries PNG information, with radial modes adding further constraining power.
Reading between the lines
- By the same mechanism, any configuration-space statistic that estimates the mean from the same data will spread a low-wavenumber signal across all scales; this could be exploited for other scale-dependent signatures (e.g., modified gravity or massive neutrinos), not just local PNG.
- The paper's logic implies a testable strategy: instead of discarding small-angle correlation data to avoid nonlinear bias, one could keep it and marginalize over nonlinear nuisance parameters, since exactly those scales carry the PNG signal.
- A natural extension is to replace the $f_{\rm sky}=1$ Gaussian covariance with the true survey window function; the paper approximates partial-sky effects, so a full window-function treatment may shift which angular bins dominate the $f_{\rm NL}$ signal.
- If validated, these forecasts suggest that wide but shallow surveys optimized for photometric galaxies at $z\gtrsim4$ could be competitive with spectroscopic surveys for $f_{\rm NL}$, because the needed angular scales are modest.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies whether the galaxy angular power spectrum C_ℓ and the angular correlation function w(θ) are equivalent for constraining local-type primordial non-Gaussianity (f_NL). It defines the integral-constrained power spectrum ~C_ℓ by subtracting the monopole (Eq. 9) and the corresponding correlation function ~w(θ) by subtracting the full-sky average (Eq. 12), and shows that the two are Legendre transforms of each other (Eq. 13). The authors argue that the f_NL signal in C_ℓ is confined to low multipoles, whereas in w(θ) the integral constraint spreads a constant offset across all angular scales, so that scale cuts break the equivalence. Using Gaussian Fisher forecasts for a six-bin LSST-like sample plus a Lyman-break-galaxy sample at 4<z<7, they report σ_f_NL ≈ 1.5 from w(θ) up to θ_max ≈ 10 deg for the high-redshift sample, and σ_f_NL ≈ 1 for all bins combined. The abstract and conclusion further claim that f_NL can be extracted from w(θ) measured in a survey with partial area coverage.
Significance. The exact derivation of the integral constraint in Eqs. (9)–(13) is clean, self-contained, and mathematically solid; it makes a useful conceptual point about how a low-multipole signal is redistributed in configuration space. The observation that a scale cut in ℓ versus a scale cut in θ breaks the equivalence differently is an original and practically relevant insight, particularly for surveys with limited sky coverage. The Fisher setup is transparent and the forecasts are reproducible from the text. However, the paper's headline claim about partial-sky surveys is not backed by the calculation, which assumes f_sky=1 and a full-sky integral constraint. Since that claim is load-bearing in the abstract and conclusion, the manuscript needs revision.
major comments (2)
- [Abstract, Section IV.A.2, Section V] The abstract and Section V state that PNG information can be extracted from w(θ) measured in a survey with partial area coverage, but the Fisher analysis never models a partial survey. The covariance in Eq. (16) assumes f_sky = 1, and the integral constraint in Eq. (12) subtracts the full-sky monopole C_{ℓ=0}/4π. The variation of θ_max in Fig. 4 and Table I corresponds to discarding pairs at large separations in an all-sky survey, not to restricting the survey footprint. A real partial survey covers a fraction f_sky of the sky; its window function removes or attenuates modes with ℓ ≲ π/θ_survey, not only ℓ = 0, and the Gaussian covariance in Eq. (16) scales roughly as 1/f_sky. For θ_max = 10 deg, f_sky = (1−cos 10°)/2 ≈ 0.0076, so the quoted σ_{f_NL} ≈ 1.5 for the z = [4,7] sample would degrade by about a factor of 11, making the forecast noncompetitive. The closing remark in Section V that window effects are 'rather straightforward' is not sufficient, because the effect is not benign: it suppresses the low-ℓ modes that carry the PNG signal.
- [Section IV.A.3 and Fig. 3] Fig. 3 presents forecasts for w(θ) as θ_min is decreased below the nonlinear scale θ_NL, which is the region to the left of the vertical dashed line in each panel, and the abstract highlights that the PNG signature in w(θ) extends below the nonlinear scale. However, the model in Eqs. (3)–(5) and the covariance in Eq. (16) are linear-theory and Gaussian quantities, and the text states that linear theory is assumed throughout. No nonlinear bias model or nonlinear covariance is provided for the small-angle bins that drive the improved σ_{f_NL} in that region. The forecast numbers in that part of Fig. 3 are therefore not a robust prediction. This does not affect the exact equivalence result in Eqs. (9)–(13), but it does affect the quantitative conclusion that the full f_NL information can be recovered from w(θ) on scales below θ_NL.
minor comments (3)
- [Table I] The columns labeled σ_f_NL (θ_max = [10,180] deg) and σ_f_NL (ℓ_min = [1,30]) are ambiguous: they contain two numbers per redshift bin without a clear indication of which number corresponds to which scale cut. The authors should split these into separate columns or label the pairs explicitly.
- [Notation throughout] The notation for the angular power spectrum and correlation function alternates between C^i_i_{gg,ℓ} and C_{gg,ℓ}, and similarly for w, which makes the text slightly harder to follow. A consistent notation, including the tilde for integral-constrained quantities, would improve readability.
- [Section IV.A.2, Eq. (18)] The shot-noise covariance formula in Eq. (18) is derived for an all-sky survey with N_pair ∝ 4π; for partial sky the pair count should scale with the survey area. This is consistent with the paper's all-sky assumption, but it should be stated explicitly near Eq. (18) so that readers do not apply it to partial-sky cases.
Circularity Check
No significant circularity: the central equivalence is an algebraic identity and the fNL forecasts are parameter-free Fisher projections from an assumed fiducial model.
full rationale
The paper's main derivation, Eqs. (9)-(13), defines the integral-constrained power spectrum and correlation function by subtracting the monopole and then notes that the two are Legendre transforms of each other. This is an identity, and the paper does not dress it up as an empirically fitted prediction. The substantive claim—that scale cuts break the equivalence because the PNG contribution to C_l is confined to low multipoles while the integral-constraint subtraction spreads that signal across all angular separations in w(theta)—follows from the Legendre transform and the k^{-2} behavior of the scale-dependent bias, not from any fitted parameter. The Fisher forecasts use a stated fiducial cosmology, assumed bias and number densities, and Gaussian covariance with fsky=1; no observational data are fitted and then renamed as a prediction. The suggestion that w(theta) from partial sky coverage could be useful is an extrapolation from the full-sky theta_max variation and is explicitly qualified by the remark that survey window effects are not modeled; that is a scope/correctness concern, not circularity. Self-citations (e.g., Kurita & Takada 2023 for window effects, Takada & Jain 2003 for pair counting) appear only as background methodology and are not load-bearing for the paper's core result. The paper is therefore self-contained against the benchmark of its own claims, and no derivation step reduces to its inputs.
Assumptions & free parameters
free parameters (3)
- Linear galaxy bias b1(z)=1+z =
1+z
- Nonlinear-scale threshold epsilon =
0.3
- PNG bias parameter p =
1
assumptions (6)
- standard math Legendre transform relation w(theta) = sum (2l+1)/(4pi) C_l P_l(cos theta) and orthogonality of Legendre polynomials
- domain assumption Local-type PNG model Phi = phi + fNL(phi^2 - <phi^2>)
- domain assumption Scale-dependent bias formula Delta b = 2(b1-p) fNL delta_c / alpha(k)
- domain assumption All-sky survey with fsky=1 and exact integral constraint as monopole subtraction
- domain assumption Gaussian covariance with shot noise, no non-Gaussian or cross-bin covariance
- domain assumption Linear theory plus variance threshold epsilon=0.3 defines the linear regime
Cite this review
Pith. "Pith review of On the equivalence between galaxy angular correlation function and power spectrum in constraining primordial non-Gaussianity." pith.science (2026). https://pith.science/paper/Z4IORIQV
@misc{pith2026250112661,
author = {Pith},
title = {Pith review of: On the equivalence between galaxy angular correlation function and power spectrum in constraining primordial non-Gaussianity},
year = {2026},
howpublished = {\url{https://pith.science/paper/Z4IORIQV}},
note = {Machine review of arXiv:2501.12661}
}
abstract
We investigate the angular power spectrum ($C_\ell)$ and angular correlation function ($w(\theta)$) of galaxy number density field in the presence of the local-type primordial non-Gaussianity (PNG), explicitly accounting for the integral constraint in an all-sky survey. We show that the PNG signature in $C_{\ell}$ is confined to low multipoles in the linear regime, whereas its signature in $w(\theta)$ extends across a wide range of angular scales, including those below the nonlinear scale. Therefore, the equivalence between $C_\ell$ and $w(\theta)$ can be violated when scale cuts of multipoles or angular scales -- for example, to mitigate systematic effects -- are applied in the analysis. Assuming samples of photometric galaxies divided into multiple redshift bins in the range $0<z<7$, we forecast the precision of constraining the PNG parameter ($f_{\rm NL}$) from the hypothetical measurements of $C_\ell$ or $w(\theta)$ assuming different scale cuts in the multipoles or angular scales, respectively. Our results imply that the PNG information can be extracted from $w(\theta)$ on relatively small angular scales such as $\lesssim 10$ degree for a high-redshift galaxy sample or from $w(\theta)$ measured in a survey with partial area coverage.
Figures
Forward citations
Cited by 1 Pith paper
-
Imprints of Large-Scale Structures in the Anisotropies of the Cosmological Gravitational Wave Background
The cosmological gravitational wave background is predicted to cross-correlate with galaxy positions through the late integrated Sachs-Wolfe effect, giving a source-discrimination handle and a forecast improvement of ...
Reference graph
Works this paper leans on
-
[1]
is the linear matter power spectrum today. We denote the linear growth factor that is normalized to the scale fac- tor during the matter-dominated era as D(z), and D0(z) ≡ D(z)/D(z = 0) is the linear growth factor normalized today as D0(z = 0) = 1. Note that the amplitude of the linear mat- ter power spectrum is normalized by the amplitude today, i.e., us...
work page 2015
-
[2]
We con- sider two different setups for the redshift range similar to the Ref
Galaxy redshift distribution We consider an LSST Y1-like survey following the DESC science requirements document (DESC SRD) [27]. We con- sider two different setups for the redshift range similar to the Ref. [28]. First, we adopt the following redshift distribution for the galaxies within redshift range z = 0 − 4: dN dz ∝ z2 exp[−(z/0.26)0.94], (14) norma...
-
[3]
We take analysis setups that would be used in an actual analysis
Covariance The covariance matrix describes the statistical errors of an observable, here the power spectrum or the correlation func- tion for each galaxy sample. We take analysis setups that would be used in an actual analysis. When we search for a PNG signal from data, we would assume a fiducial ΛCDM model with no PNG, i.e. fNL = 0 , to model the covaria...
-
[4]
Fisher matrix Using the covariance matrix, we can compute the Fisher matrix for a hypothetical measurement of C ii gg,ℓ or wii gg in the i-th redshift bin, to forecast the precision of parameter estima- tion: F i αβ ≡ ∂dT ∂pα C−1 ∂d ∂pβ (19) where d is the data vector that is either of the power spec- trum or the correlation function, i.e. C ii gg,ℓ or wi...
- [5]
-
[6]
A. H. Guth and S. Y . Pi, Phys. Rev. Lett. 49, 1110 (1982)
work page 1982
-
[7]
Maldacena, Journal of High Energy Physics 2003, 013 (2003), astro-ph/0210603
J. Maldacena, Journal of High Energy Physics 2003, 013 (2003), astro-ph/0210603
arXiv 2003
-
[8]
N. Bartolo, E. Komatsu, S. Matarrese, and A. Riotto, Phys. Rep. 402, 103 (2004), arXiv:astro-ph/0406398
arXiv 2004
Show all 42 references
-
[9]
Dalal, O
N. Dalal, O. Dor ´e, D. Huterer, and A. Shirokov, Physical Re- view D 77, 123514 (2008), 0710.4560
2008 arXiv
-
[10]
Dor ´e, J
O. Dor ´e, J. Bock, M. Ashby, P. Capak, A. Cooray, R. de Putter, T. Eifler, N. Flagey, Y . Gong, S. Habib, et al., arXiv e-prints arXiv:1412.4872 (2014), 1412.4872
2014 arXiv
-
[11]
Creminelli and M
P. Creminelli and M. Zaldarriaga, JCAP 2004, 006 (2004), astro-ph/0407059
2004 arXiv
-
[12]
Pajer, F
E. Pajer, F. Schmidt, and M. Zaldarriaga, Phys. Rev. D 88, 083502 (2013), 1305.0824
2013 arXiv
- [13]
-
[14]
Desjacques, D
V . Desjacques, D. Jeong, and F. Schmidt, Phys. Rep. 733, 1 (2018), 1611.09787
2018 arXiv
-
[15]
Slosar, C
A. Slosar, C. Hirata, U. Seljak, S. Ho, and N. Padmanabhan, JCAP 8, 31 (2008), 0805.3580
2008 arXiv
-
[16]
S. Ho, A. Cuesta, H.-J. Seo, R. de Putter, A. J. Ross, M. White, N. Padmanabhan, S. Saito, D. J. Schlegel, E. Schlafly, et al., Astrophys. J. 761, 14 (2012), 1201.2137
2012 arXiv
-
[17]
A. J. Ross, W. J. Percival, A. Carnero, G.-b. Zhao, M. Man- era, A. Raccanelli, E. Aubourg, D. Bizyaev, H. Brewington, J. Brinkmann, et al., Mon. Not. Roy. Astron. Soc. 428, 1116 (2013), 1208.1491
2013 arXiv
-
[18]
Rezaie, A
M. Rezaie, A. J. Ross, H.-J. Seo, E.-M. Mueller, W. J. Percival, G. Merz, R. Katebi, R. C. Bunescu, J. Bautista, J. R. Brown- stein, et al., Mon. Not. Roy. Astron. Soc. 506, 3439 (2021), 2106.13724
2021 arXiv
- [19]
-
[20]
M. S. Cagliari, E. Castorina, M. Bonici, and D. Bianchi, arXiv e-prints arXiv:2309.15814 (2023), 2309.15814
2023 arXiv
-
[21]
Riquelme, S
W. Riquelme, S. Avila, J. Garc ´ıa-Bellido, A. Porredon, I. Fer- rero, K. C. Chan, R. Rosenfeld, H. Camacho, A. G. Adame, A. Carnero Rosell, et al., Mon. Not. Roy. Astron. Soc.523, 603 (2023), 2209.07187. 9
2023 arXiv
-
[22]
Rezaie, A
M. Rezaie, A. J. Ross, H.-J. Seo, H. Kong, A. Porredon, L. Samushia, E. Chaussidon, A. Krolewski, A. de Mattia, F. Beutler, et al., arXiv e-prints arXiv:2307.01753 (2023), 2307.01753
2023 arXiv
-
[23]
Chaussidon, C
E. Chaussidon, C. Y `eche, A. de Mattia, C. Payerne, P. McDon- ald, A. J. Ross, S. Ahlen, D. Bianchi, D. Brooks, E. Burtin, et al., arXiv e-prints arXiv:2411.17623 (2024), 2411.17623
2024 arXiv
-
[24]
Biagetti, Galaxies 7, 71 (2019), 1906.12244
M. Biagetti, Galaxies 7, 71 (2019), 1906.12244
2019 arXiv
-
[25]
X. Fang, E. Krause, T. Eifler, and N. MacCrann, JCAP 2020, 010 (2020), 1911.11947
2020 arXiv
-
[26]
T. M. C. Abbott, M. Aguena, A. Alarcon, S. Allam, O. Alves, A. Amon, F. Andrade-Oliveira, J. Annis, S. Avila, D. Bacon, et al., Phys. Rev. D 105, 023520 (2022), 2105.13549
2022 arXiv
-
[27]
Miyatake, S
H. Miyatake, S. Sugiyama, M. Takada, T. Nishimichi, X. Li, M. Shirasaki, S. More, Y . Kobayashi, A. J. Nishizawa, M. M. Rau, et al., Phys. Rev. D 108, 123517 (2023), 2304.00704
2023 arXiv
-
[28]
Terasawa, X
R. Terasawa, X. Li, M. Takada, T. Nishimichi, S. Tanaka, S. Sugiyama, T. Kurita, T. Zhang, M. Shirasaki, R. Takahashi, et al., arXiv e-prints arXiv:2403.20323 (2024), 2403.20323
2024 arXiv
-
[29]
Jeong, F
D. Jeong, F. Schmidt, and C. M. Hirata, Phys. Rev. D 85, 023504 (2012), 1107.5427
2012 arXiv
-
[30]
Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, et al., Astronomy & Astrophysics 594, A13 (2016), 1502.01589
2016 arXiv
-
[31]
T. L. D. E. S. Collaboration, R. Mandelbaum, T. Eifler, R. Hlo ˇzek, T. Collett, E. Gawiser, D. Scolnic, D. Alonso, H. Awan, R. Biswas, et al., The lsst dark energy science collaboration (desc) science requirements document (2021), 1809.01669
2021 arXiv
-
[32]
Schmittfull and U
M. Schmittfull and U. Seljak, Phys. Rev. D 97, 123540 (2018), 1710.09465
2018 arXiv
-
[33]
Harikane, Y
Y . Harikane, Y . Ono, M. Ouchi, C. Liu, M. Sawicki, T. Shibuya, P. S. Behroozi, W. He, K. Shimasaku, S. Arnouts, et al., Astro- phys. J. Suppl. 259, 20 (2022), 2108.01090
2022 arXiv
-
[34]
Crocce, A
M. Crocce, A. Cabr ´e, and E. Gazta ˜naga, Mon. Not. Roy. As- tron. Soc. 414, 329 (2011), 1004.4640
2011 arXiv
-
[35]
Takada and B
M. Takada and B. Jain, Mon. Not. Roy. Astron. Soc. 344, 857 (2003), astro-ph/0304034
2003 arXiv
-
[36]
Joachimi, P
B. Joachimi, P. Schneider, and T. Eifler, Astronomy & Astro- physics 477, 43 (2008), 0708.0387
2008 arXiv
-
[37]
Senatore and M
L. Senatore and M. Zaldarriaga, Journal of Cosmol- ogy and Astroparticle Physics 2015, 013–013 (2015), ISSN 1475-7516, URL http://dx.doi.org/10.1088/ 1475-7516/2015/02/013
2015
-
[38]
G. F. R. Ellis and J. E. Baldwin, Mon. Not. Roy. Astron. Soc. 206, 377 (1984)
1984
-
[39]
Blake and J
C. Blake and J. Wall, Nature (London) 416, 150 (2002), astro- ph/0203385
2002
-
[40]
Y . Itoh, K. Yahata, and M. Takada, Phys. Rev. D 82, 043530 (2010), 0912.1460
2010 arXiv
-
[41]
Takada, R
M. Takada, R. S. Ellis, M. Chiba, J. E. Greene, H. Aihara, N. Arimoto, K. Bundy, J. Cohen, O. Dor ´e, G. Graves, et al., Publ. Astron. Soc. Japan 66, R1 (2014), 1206.0737
2014 arXiv
-
[42]
N. Hand, Y . Feng, F. Beutler, Y . Li, C. Modi, U. Seljak, and Z. Slepian, The Astronomical Journal 156, 160 (2018), ISSN 1538-3881, URL http://dx.doi.org/10.3847/ 1538-3881/aadae0
2018
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.