pith. sign in

arxiv: 2504.20992 · v2 · pith:IMWSGHOCnew · submitted 2025-04-29 · 🌌 astro-ph.CO · gr-qc

Constraints on primordial non-Gaussianity from Quaia

Pith reviewed 2026-05-22 18:13 UTC · model grok-4.3

classification 🌌 astro-ph.CO gr-qc
keywords primordial non-Gaussianityf_NLquasarsCMB lensingangular clusteringlarge-scale structureQuaia catalogtomographic bins
0
0 comments X

The pith

Quasar clustering combined with CMB lensing cross-correlations constrains local primordial non-Gaussianity to f_NL = -20.5^{+19.0}_{-18.1}

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

The paper analyses the large-scale angular clustering of quasars in the Quaia catalog together with their cross-correlation to maps of CMB lensing convergence. It targets the scale-dependent bias that local-type primordial non-Gaussianity would produce on biased tracers. Using the universality relation with response parameter p_φ=1, the combined auto- and cross-correlation signals in two tomographic redshift bins produce the reported constraint on f_NL. A reader would care because this limit tests early-universe inflation models that differ in their predicted level of non-Gaussianity, and the result is presented as the tightest obtained so far from angular statistics of this kind.

Core claim

Using the universality relation to characterise the response of the quasar overdensity to PNG (p_φ=1), the combination of quasar auto-correlation and its cross-correlation with CMB lensing in two tomographic redshift bins yields f_NL = -20.5^{+19.0}_{-18.1} at 68% C.L. The same data give f_NL = -28.7^{+26.1}_{-24.6} when p_φ=1.6 is assumed instead. CMB lensing cross-correlations alone produce consistent though looser bounds of f_NL = -13.8^{+26.7}_{-25.0} and f_NL = -15.6^{+42.3}_{-34.8} for the two response values. These are stated to be the tightest constraints on f_NL from angular clustering statistics and CMB lensing cross-correlations to date.

What carries the argument

The scale-dependent bias that local primordial non-Gaussianity induces on the overdensity of biased tracers such as quasars, modelled through the universality relation with fixed response parameter p_φ

If this is right

  • Including cross-correlations between the two tomographic redshift bins would further reduce the uncertainty on f_NL.
  • CMB lensing cross-correlations by themselves already produce constraints consistent with the full combination but with larger errors.
  • The same methodology applied to other large-volume tracers would be expected to deliver comparable or tighter limits on local PNG.

Where Pith is reading between the lines

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

  • Repeating the analysis on independent quasar samples at similar redshifts would test whether the central value of f_NL shifts when different selection effects are present.
  • If the response parameter p_φ itself varies with quasar luminosity or redshift, the central value of f_NL would move while the error bars might remain similar.
  • Joint analysis with bispectrum measurements from the CMB could reveal whether the clustering and three-point statistics favour the same value of f_NL.

Load-bearing premise

The main potential sources of systematic contamination in the Quaia sample can be mapped out and their impact mitigated sufficiently to trust the large-scale clustering signal.

What would settle it

A future measurement that finds the large-scale quasar clustering amplitude or its cross-correlation with CMB lensing to differ from the modelled PNG response after all known systematics have been subtracted would invalidate the reported f_NL bounds.

read the original abstract

We analyse the large-scale angular clustering of quasars in the \gaia-\unwise quasar catalog, \quaia, and their cross-correlation with maps of the lensing convergence of the Cosmic Microwave Background (CMB), to constrain the level of primordial non-Gaussianity (PNG). Specifically, we target the scale-dependent bias that would be induced by PNG on biased tracers of the matter inhomogeneities on large scales. The \quaia sample is particularly well suited for this analysis, given the large effective volume covered, and our ability to map out the main potential sources of systematic contamination and mitigate their impact. Using the universality relation to characterise the response of the quasar overdensity to PNG ($p_\phi=1$), we report constraints on the local-type PNG parameter $f_{\rm NL}$ of $f_{\rm NL} =-20.5^{+19.0}_{-18.1}$ (68\% C.L.) by combining the quasar auto-correlation and its cross-correlation with CMB lensing in two tomographic redshift bins (or $f_{\rm NL} =-28.7^{+26.1}_{-24.6}$ if assuming a lower response for quasars, $p_\phi=1.6$). The error on $f_{\rm NL}$ can be further improved if the cross-correlation between the tomographic redshift bins is included. Using the CMB lensing cross-correlations alone, we find $f_{\rm NL} =-13.8^{+26.7}_{-25.0}$ and $f_{\rm NL} = -15.6^{+42.3}_{-34.8}$ for $p_\phi=1$ and $p_\phi=1.6$ respectively. These are the tightest constraints on $f_{\rm NL}$ to date from angular clustering statistics and cross-correlations with CMB lensing.

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 analyzes the large-scale angular clustering of quasars in the Gaia-unWISE Quaia catalog and their cross-correlation with CMB lensing convergence maps to constrain local-type primordial non-Gaussianity. Using two tomographic redshift bins and the universality relation with p_φ=1 (or 1.6), the authors report f_NL = -20.5^{+19.0}_{-18.1} (68% C.L.) from the combination of quasar auto-correlation and CMB-lensing cross-correlation, with additional results from cross-correlations alone and when including inter-bin correlations; they claim these are the tightest constraints on f_NL from angular clustering statistics and CMB lensing cross-correlations to date.

Significance. If the reported constraints hold after rigorous validation, the result provides competitive limits on f_NL that can help discriminate among single-field and multi-field inflationary scenarios. The large effective volume of the Quaia sample and the use of cross-correlations with CMB lensing to reduce certain systematics are positive aspects. The work also explores the sensitivity to the response parameter p_φ, which is a useful check.

major comments (2)
  1. [Abstract and associated methods description] The central claim that the observed large-scale angular power spectra are driven by PNG-induced scale-dependent bias rather than residuals rests on the assertion that main systematics (stellar contamination, dust extinction, selection effects) have been mapped and mitigated. However, without explicit quantification of residual contamination levels in the lowest multipoles of the two tomographic bins (or null tests demonstrating that any residual 1/k^2-like signature is sub-dominant to the statistical uncertainty), the reliability of the combined f_NL = -20.5^{+19.0}_{-18.1} constraint cannot be fully assessed.
  2. [Results section describing the combined constraints] The covariance estimation between the quasar auto-power spectra and the cross-power spectra with CMB lensing, as well as between the two tomographic redshift bins, is not described in sufficient detail to verify that the reported error bars on f_NL properly account for correlations; this directly affects the robustness of the combined constraint and the claim that including inter-bin cross-correlations further improves the error.
minor comments (2)
  1. [Abstract] The abstract states that the error on f_NL can be further improved if the cross-correlation between the tomographic redshift bins is included, but the quantitative improvement and the corresponding constraint value are not provided; adding this would strengthen the presentation.
  2. [Introduction or methods] Notation for the response parameter is given as p_φ in the abstract but should be consistently defined with its physical meaning (e.g., relating to the universality relation for the quasar bias response) at first use in the main text.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for the constructive comments. We address each major comment below and indicate the revisions we will make to improve the presentation and robustness of the analysis.

read point-by-point responses
  1. Referee: [Abstract and associated methods description] The central claim that the observed large-scale angular power spectra are driven by PNG-induced scale-dependent bias rather than residuals rests on the assertion that main systematics (stellar contamination, dust extinction, selection effects) have been mapped and mitigated. However, without explicit quantification of residual contamination levels in the lowest multipoles of the two tomographic bins (or null tests demonstrating that any residual 1/k^2-like signature is sub-dominant to the statistical uncertainty), the reliability of the combined f_NL = -20.5^{+19.0}_{-18.1} constraint cannot be fully assessed.

    Authors: We agree that explicit quantification of residual contamination is important for validating the f_NL constraints. The manuscript already describes the mapping and mitigation of stellar contamination, dust extinction, and selection effects (Sections 3 and 4), including the construction of masks and weighting schemes based on Gaia and unWISE data. However, we acknowledge that the presentation would benefit from additional explicit numbers for residual levels at the lowest multipoles (ℓ ≲ 10) in each tomographic bin and from dedicated null tests isolating any residual scale-dependent signature. We will revise the text and add a new subsection (or appendix) presenting these quantifications and null-test results, confirming that any residual 1/k^2-like contamination remains sub-dominant to the statistical uncertainty. revision: yes

  2. Referee: [Results section describing the combined constraints] The covariance estimation between the quasar auto-power spectra and the cross-power spectra with CMB lensing, as well as between the two tomographic redshift bins, is not described in sufficient detail to verify that the reported error bars on f_NL properly account for correlations; this directly affects the robustness of the combined constraint and the claim that including inter-bin cross-correlations further improves the error.

    Authors: We thank the referee for noting the need for greater detail on the covariance. The covariance matrix is constructed from a hybrid approach combining analytical Gaussian terms with mock-based estimates that explicitly include the cross-covariance between auto- and cross-spectra as well as inter-bin correlations (Section 5.2). To address the comment, we will expand this section with a more complete description of the mock generation, the block structure of the covariance matrix, and validation tests (e.g., comparison with jackknife estimates). This will make clear how the reported uncertainties on f_NL incorporate all relevant correlations and why the inclusion of inter-bin cross-correlations tightens the constraints. revision: yes

Circularity Check

0 steps flagged

No circularity in derivation of f_NL constraints from Quaia clustering

full rationale

The paper obtains f_NL constraints by fitting the observed large-scale angular auto- and cross-power spectra (quasar auto-correlation plus CMB lensing cross-correlation in two tomographic bins) to a theoretical model that includes the PNG-induced scale-dependent bias term. The universality relation (with p_φ=1 or 1.6) is invoked as a standard external theoretical input to translate the bias response into f_NL; it is not defined in terms of the present data or result. No equation reduces the reported constraint to a fitted parameter by construction, no load-bearing step relies on a self-citation chain for uniqueness, and the central result remains a direct data-driven likelihood fit rather than a renaming or tautological re-expression of inputs. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The analysis rests on standard cosmological modeling of PNG effects on biased tracers and the assumption that systematics in the Quaia sample are controllable; no new entities are postulated.

free parameters (1)
  • p_phi
    Response of quasar overdensity to PNG, set to 1 via universality relation or alternatively to 1.6; not fitted from the current data but chosen as modeling assumption.
axioms (1)
  • domain assumption Universality relation characterizing the response of quasar overdensity to primordial non-Gaussianity
    Invoked to set p_phi=1 for the main constraint.

pith-pipeline@v0.9.0 · 5878 in / 1460 out tokens · 109333 ms · 2026-05-22T18:13:31.315309+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 4 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Euclid preparation: Testing multi-field inflation with galaxy power spectrum and bispectrum

    astro-ph.CO 2026-05 conditional novelty 5.0

    Validates redshift-space power spectrum and bispectrum analysis on Abacus-PNG mocks to recover unbiased f_NL constraints for Euclid spectroscopic sample.

  2. Constraining primordial non-Gaussianity from DESI DR1 quasars and Planck PR4 CMB Lensing

    astro-ph.CO 2025-12 unverdicted novelty 5.0

    Cross-correlation of DESI DR1 quasars with Planck PR4 CMB lensing constrains local f_NL to 2^{+28}_{-34} (p=1.6) or 6^{+20}_{-24} (p=1.0), tightening previous limits by 35%.

  3. New constraints on primordial non-Gaussianity from large-scale cross-correlations of CMB lensing and the cosmic infrared background

    astro-ph.CO 2026-05 unverdicted novelty 4.0

    Dust-cleaned CIB and CMB lensing cross-correlations yield f_NL^local = 43 ± 23, tightening constraints on local primordial non-Gaussianity.

  4. Tracing the high-z cosmic web with Quaia: catalogues of voids and clusters in the quasar distribution

    astro-ph.CO 2025-09 unverdicted novelty 4.0

    Using 708,483 quasars at 0.8<z<2.2 from Quaia, the authors identify 12,842 voids and 41,111 clusters whose radii, densities and profiles agree with 50 mock catalogues to within 5-10%, with no ultra-large structures ex...

Reference graph

Works this paper leans on

79 extracted references · 79 canonical work pages · cited by 4 Pith papers · 31 internal anchors

  1. [1]

    Non-Gaussianity from Inflation: Theory and Observations

    N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto,Non-Gaussianity from inflation: theory and observations, Phys. Rep.402 (2004) 103 [astro-ph/0406398]

  2. [2]

    TASI Lectures on Inflation

    D. Baumann,TASI Lectures on Inflation, arXiv e-prints (2009) arXiv:0907.5424 [0907.5424]

  3. [3]

    Planck 2018 results. IX. Constraints on primordial non-Gaussianity

    Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi et al., Planck 2018 results. IX. Constraints on primordial non-Gaussianity, A&A 641 (2020) A9 [1905.05697]

  4. [4]

    G. Jung, M. Citran, B. van Tent, L. Dumilly and N. Aghanim,Constraints on primordial non-Gaussianity from Planck PR4 data, arXiv e-prints (2025) arXiv:2504.00884 [2504.00884]

  5. [5]

    Single field consistency relation for the 3-point function

    P. Creminelli and M. Zaldarriaga,A single-field consistency relation for the three-point function, J. Cosmology Astropart. Phys.2004 (2004) 006 [astro-ph/0407059]

  6. [6]

    The imprints of primordial non-gaussianities on large-scale structure: scale dependent bias and abundance of virialized objects

    N. Dalal, O. Doré, D. Huterer and A. Shirokov,Imprints of primordial non-Gaussianities on large-scale structure: Scale-dependent bias and abundance of virialized objects, Phys. Rev. D77 (2008) 123514 [0710.4560]

  7. [7]

    Castorina, N

    E. Castorina, N. Hand, U. Seljak, F. Beutler, C.-H. Chuang, C. Zhao et al.,Redshift-weighted constraints on primordial non-Gaussianity from the clustering of the eBOSS DR14 quasars in Fourier space, J. Cosmology Astropart. Phys.2019 (2019) 010 [1904.08859]

  8. [8]

    Krolewski, W.J

    A. Krolewski, W.J. Percival, S. Ferraro, E. Chaussidon, M. Rezaie, J.N. Aguilar et al., Constraining primordial non-Gaussianity from DESI quasar targets and Planck CMB lensing, J. Cosmology Astropart. Phys.2024 (2024) 021 [2305.07650]

  9. [9]

    Cagliari, E

    M.S. Cagliari, E. Castorina, M. Bonici and D. Bianchi,Optimal constraints on Primordial non-Gaussianity with the eBOSS DR16 quasars in Fourier space, J. Cosmology Astropart. Phys. 2024 (2024) 036 [2309.15814]

  10. [10]

    McCarthy, M.S

    F. McCarthy, M.S. Madhavacheril and A.S. Maniyar,Constraints on primordial non-Gaussianity from halo bias measured through CMB lensing cross-correlations, Phys. Rev. D108 (2023) 083522 [2210.01049]

  11. [11]

    Radio Galaxy populations and the multi-tracer technique: pushing the limits on primordial non-Gaussianity

    L.D. Ferramacho, M.G. Santos, M.J. Jarvis and S. Camera,Radio galaxy populations and the multitracer technique: pushing the limits on primordial non-Gaussianity, MNRAS 442 (2014) 2511 [1402.2290]

  12. [12]

    Bermejo-Climent, R

    J.R. Bermejo-Climent, R. Demina, A. Krolewski, E. Chaussidon, M. Rezaie, S. Ahlen et al., Constraints on primordial non-Gaussianity from the cross-correlation of DESI Luminous Red Galaxies and P lanck CMB lensing, arXiv e-prints (2024) arXiv:2412.10279 [2412.10279]. 6https://data.cmb-s4.org/planck_pr4.html – 30 –

  13. [13]

    Chaussidon , C

    E. Chaussidon, C. Yèche, A. de Mattia, C. Payerne, P. McDonald, A.J. Ross et al., Constraining primordial non-Gaussianity with DESI 2024 LRG and QSO samples, arXiv e-prints (2024) arXiv:2411.17623 [2411.17623]

  14. [14]

    Storey-Fisher, D.W

    K. Storey-Fisher, D.W. Hogg, H.-W. Rix, A.-C. Eilers, G. Fabbian, M.R. Blanton et al.,Quaia, the Gaia-unWISE Quasar Catalog: An All-sky Spectroscopic Quasar Sample, ApJ 964 (2024) 69 [2306.17749]

  15. [15]

    CMB lensing from Planck PR4 maps

    J. Carron, M. Mirmelstein and A. Lewis,CMB lensing from Planck PR4 maps, J. Cosmology Astropart. Phys. 2022 (2022) 039 [2206.07773]

  16. [16]

    Measuring cosmological parameters with galaxy surveys

    M. Tegmark,Measuring Cosmological Parameters with Galaxy Surveys, Phys. Rev. Lett.79 (1997) 3806 [astro-ph/9706198]

  17. [17]

    Vallenari, A.G.A

    Gaia Collaboration, A. Vallenari, A.G.A. Brown, T. Prusti, J.H.J. de Bruijne, F. Arenou et al., Gaia Data Release 3. Summary of the content and survey properties, A&A 674 (2023) A1 [2208.00211]

  18. [18]

    Meisner, D

    A.M. Meisner, D. Lang, E.F. Schlafly and D.J. Schlegel,unWISE Coadds: The Five-year Data Set, PASP131 (2019) 124504 [1909.05444]

  19. [19]

    A.J. Ross, J. Bautista, R. Tojeiro, S. Alam, S. Bailey, E. Burtin et al.,The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Large-scale structure catalogues for cosmological analysis, MNRAS 498 (2020) 2354 [2007.09000]

  20. [20]

    Alonso, G

    D. Alonso, G. Fabbian, K. Storey-Fisher, A.-C. Eilers, C. García-García, D.W. Hogg et al., Constraining cosmology with the Gaia-unWISE Quasar Catalog and CMB lensing: structure growth, J. Cosmology Astropart. Phys.2023 (2023) 043 [2306.17748]

  21. [21]

    Alonso, O

    D. Alonso, O. Hetmantsev, G. Fabbian, A. Slosar and K. Storey-Fisher,Measurement of the power spectrum turnover scale from the cross-correlation between CMB lensing and Quaia, arXiv e-prints (2024) arXiv:2410.24134 [2410.24134]

  22. [22]

    Chiang,Corrected SFD: A More Accurate Galactic Dust Map with Minimal Extragalactic Contamination, ApJ 958 (2023) 118 [2306.03926]

    Y.-K. Chiang,Corrected SFD: A More Accurate Galactic Dust Map with Minimal Extragalactic Contamination, ApJ 958 (2023) 118 [2306.03926]

  23. [23]

    Williams et al.in preparation.(2025)

    A. Williams et al.in preparation.(2025)

  24. [24]

    Piccirilli, G

    G. Piccirilli, G. Fabbian, D. Alonso, K. Storey-Fisher, J. Carron, A. Lewis et al.,Growth history and quasar bias evolution at z < 3 from Quaia, J. Cosmology Astropart. Phys.2024 (2024) 012 [2402.05761]

  25. [25]

    MASTER of the CMB Anisotropy Power Spectrum: A Fast Method for Statistical Analysis of Large and Complex CMB Data Sets

    E. Hivon, K.M. Górski, C.B. Netterfield, B.P. Crill, S. Prunet and F. Hansen,MASTER of the Cosmic Microwave Background Anisotropy Power Spectrum: A Fast Method for Statistical Analysis of Large and Complex Cosmic Microwave Background Data Sets, ApJ 567 (2002) 2 [astro-ph/0105302]

  26. [26]

    A unified pseudo-$C_\ell$ framework

    D. Alonso, J. Sanchez, A. Slosar and LSST Dark Energy Science Collaboration,A unified pseudo-Cℓ framework, MNRAS 484 (2019) 4127 [1809.09603]

  27. [27]

    Estimating the large-scale angular power spectrum in the presence of systematics: a case study of Sloan Digital Sky Survey quasars

    B. Leistedt, H.V. Peiris, D.J. Mortlock, A. Benoit-Lévy and A. Pontzen,Estimating the large-scale angular power spectrum in the presence of systematics: a case study of Sloan Digital Sky Survey quasars, MNRAS 435 (2013) 1857 [1306.0005]

  28. [28]

    Unbiased pseudo-Cl power spectrum estimation with mode projection

    F. Elsner, B. Leistedt and H.V. Peiris,Unbiased pseudo-Cℓ power spectrum estimation with mode projection, MNRAS 465 (2017) 1847 [1609.03577]

  29. [29]

    Adachi, M.A.O

    Polarbear Collaboration, S. Adachi, M.A.O. Aguilar Faúndez, K. Arnold, C. Baccigalupi, D. Barron et al.,A Measurement of the Degree-scale CMB B-mode Angular Power Spectrum with POLARBEAR, ApJ 897 (2020) 55 [1910.02608]

  30. [30]

    Hervías-Caimapo, K

    C. Hervías-Caimapo, K. Wolz, A. La Posta, S. Azzoni, D. Alonso, K. Arnold et al.,The Simons – 31 – Observatory: Validation of reconstructed power spectra from simulated filtered maps for the Small Aperture Telescope survey, arXiv e-prints (2025) [ 2502.00946]

  31. [31]

    The Atacama Cosmology Telescope: DR6 Power Spectra, Likelihoods and $\Lambda$CDM Parameters

    T. Louis, A. La Posta, Z. Atkins, H.T. Jense, I. Abril-Cabezas, G.E. Addison et al.,The Atacama Cosmology Telescope: DR6 Power Spectra, Likelihoods andΛCDM Parameters, arXiv e-prints (2025) [ 2503.14452]

  32. [32]

    Cornish et al.in preparation.(2025)

    T. Cornish et al.in preparation.(2025)

  33. [33]

    Mittal, O.T

    V. Mittal, O.T. Oayda and G.F. Lewis,The cosmic dipole in the Quaia sample of quasars: a Bayesian analysis, MNRAS 527 (2024) 8497 [2311.14938]

  34. [34]

    How to measure CMB power spectra without losing information

    M. Tegmark,How to measure CMB power spectra without losing information, Phys. Rev. D55 (1997) 5895 [astro-ph/9611174]

  35. [35]

    Quadratic estimator for CMB cross-correlation

    S. Vanneste, S. Henrot-Versillé, T. Louis and M. Tristram,Quadratic estimator for CMB cross-correlation, Phys. Rev. D98 (2018) 103526 [1807.02484]

  36. [36]

    Rybicki and W.H

    G.B. Rybicki and W.H. Press,Interpolation, Realization, and Reconstruction of Noisy, Irregularly Sampled Data, ApJ 398 (1992) 169

  37. [37]

    García-García, D

    C. García-García, D. Alonso and E. Bellini,Disconnected pseudo-Cl covariances for projected large-scale structure data, J. Cosmology Astropart. Phys.2019 (2019) 043 [1906.11765]

  38. [38]

    Clustering of quasars in SDSS-IV eBOSS : study of potential systematics and bias determination

    P. Laurent, S. Eftekharzadeh, J.-M. Le Goff, A. Myers, E. Burtin, M. White et al.,Clustering of quasars in SDSS-IV eBOSS: study of potential systematics and bias determination, J. Cosmology Astropart. Phys.2017 (2017) 017 [1705.04718]

  39. [39]

    Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix

    J. Hartlap, P. Simon and P. Schneider,Why your model parameter confidences might be too optimistic. Unbiased estimation of the inverse covariance matrix, A&A 464 (2007) 399 [astro-ph/0608064]

  40. [40]

    Press and P

    W.H. Press and P. Schechter,Formation of Galaxies and Clusters of Galaxies by Self-Similar Gravitational Condensation, ApJ 187 (1974) 425

  41. [41]

    Constraints on local primordial non-Gaussianity from large scale structure

    A. Slosar, C. Hirata, U. Seljak, S. Ho and N. Padmanabhan,Constraints on local primordial non-Gaussianity from large scale structure, J. Cosmology Astropart. Phys.2008 (2008) 031 [0805.3580]

  42. [42]

    Barreira, G

    A. Barreira, G. Cabass, F. Schmidt, A. Pillepich and D. Nelson,Galaxy bias and primordial non-Gaussianity: insights from galaxy formation simulations with IllustrisTNG, J. Cosmology Astropart. Phys. 2020 (2020) 013 [2006.09368]

  43. [43]

    Barreira,Predictions for local PNG bias in the galaxy power spectrum and bispectrum and the consequences for fN L constraints, J

    A. Barreira,Predictions for local PNG bias in the galaxy power spectrum and bispectrum and the consequences for fN L constraints, J. Cosmology Astropart. Phys.2022 (2022) 033 [2107.06887]

  44. [44]

    What galaxy surveys really measure

    C. Bonvin and R. Durrer,What galaxy surveys really measure, Phys. Rev. D84 (2011) 063505 [1105.5280]

  45. [45]

    The linear power spectrum of observed source number counts

    A. Challinor and A. Lewis,Linear power spectrum of observed source number counts, Phys. Rev. D84 (2011) 043516 [1105.5292]

  46. [46]

    Ultra large-scale cosmology in next-generation experiments with single tracers

    D. Alonso, P. Bull, P.G. Ferreira, R. Maartens and M.G. Santos,Ultra-large-scale Cosmology in Next-generation Experiments with Single Tracers, ApJ 814 (2015) 145 [1505.07596]

  47. [47]

    Constraining ultra large-scale cosmology with multiple tracers in optical and radio surveys

    D. Alonso and P.G. Ferreira,Constraining ultralarge-scale cosmology with multiple tracers in optical and radio surveys, Phys. Rev. D92 (2015) 063525 [1507.03550]

  48. [48]

    Guedezounme, S

    S.L. Guedezounme, S. Jolicoeur and R. Maartens,Primordial non-Gaussianity – the effects of relativistic and wide-angle corrections to the power spectrum, arXiv e-prints (2024) arXiv:2412.06553 [2412.06553]

  49. [49]

    Unequal-Time Correlators for Cosmology

    T.D. Kitching and A.F. Heavens,Unequal-time correlators for cosmology, Phys. Rev. D95 (2017) 063522 [1612.00770]. – 32 –

  50. [50]

    Chisari and A

    N.E. Chisari and A. Pontzen,Unequal time correlators and the Zel’dovich approximation, Phys. Rev. D100 (2019) 023543 [1905.02078]

  51. [51]

    Limber,The Analysis of Counts of the Extragalactic Nebulae in Terms of a Fluctuating Density Field., ApJ 117 (1953) 134

    D.N. Limber,The Analysis of Counts of the Extragalactic Nebulae in Terms of a Fluctuating Density Field., ApJ 117 (1953) 134

  52. [52]

    Extended Limber Approximation

    M. LoVerde and N. Afshordi,Extended Limber approximation, Phys. Rev. D78 (2008) 123506 [0809.5112]

  53. [53]

    Chisari, D

    N.E. Chisari, D. Alonso, E. Krause, C.D. Leonard, P. Bull, J. Neveu et al.,Core Cosmology Library: Precision Cosmological Predictions for LSST, ApJS 242 (2019) 2 [1812.05995]

  54. [54]

    Leonard, T

    C.D. Leonard, T. Ferreira, X. Fang, R. Reischke, N. Schoeneberg, T. Tröster et al.,The N5K Challenge: Non-Limber Integration for LSST Cosmology, The Open Journal of Astrophysics6 (2023) 8 [2212.04291]

  55. [55]

    X. Fang, E. Krause, T. Eifler and N. MacCrann,Beyond Limber: efficient computation of angular power spectra for galaxy clustering and weak lensing, J. Cosmology Astropart. Phys. 2020 (2020) 010 [1911.11947]

  56. [56]

    Likelihood Analysis of CMB Temperature and Polarization Power Spectra

    S. Hamimeche and A. Lewis,Likelihood analysis of CMB temperature and polarization power spectra, Phys. Rev. D77 (2008) 103013 [0801.0554]

  57. [57]

    Planck 2018 results. VI. Cosmological parameters

    Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi et al., Planck 2018 results. VI. Cosmological parameters, A&A 641 (2020) A6 [1807.06209]

  58. [58]

    emcee: The MCMC Hammer

    D. Foreman-Mackey, D.W. Hogg, D. Lang and J. Goodman,emcee: The MCMC Hammer, PASP 125 (2013) 306 [1202.3665]

  59. [59]

    Recovering 3D clustering information with angular correlations

    J. Asorey, M. Crocce, E. Gaztañaga and A. Lewis,Recovering 3D clustering information with angular correlations, MNRAS 427 (2012) 1891 [1207.6487]

  60. [60]

    Mueller, M

    E.-M. Mueller, M. Rezaie, W.J. Percival, A.J. Ross, R. Ruggeri, H.-J. Seo et al.,The clustering of galaxies in the completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Primordial non-Gaussianity in Fourier Space, arXiv e-prints (2021) arXiv:2106.13725 [2106.13725]

  61. [61]

    Cagliari, M

    M.S. Cagliari, M. Barberi-Squarotti, K. Pardede, E. Castorina and G. D’Amico,Bispectrum constraints on Primordial non-Gaussianities with the eBOSS DR16 quasars, arXiv e-prints (2025) arXiv:2502.14758 [2502.14758]

  62. [62]

    Rezaie, A.J

    M. Rezaie, A.J. Ross, H.-J. Seo, H. Kong, A. Porredon, L. Samushia et al.,Local primordial non-Gaussianity from the large-scale clustering of photometric DESI luminous red galaxies, MNRAS 532 (2024) 1902 [2307.01753]

  63. [63]

    Science-driven 3D data compression

    D. Alonso,Science-driven 3D data compression, MNRAS 473 (2018) 4306 [1707.08950]

  64. [64]

    Parameter constraints from cross-correlation of CMB lensing with galaxy clustering

    M. Schmittfull and U. Seljak,Parameter constraints from cross-correlation of CMB lensing with galaxy clustering, Phys. Rev. D97 (2018) 123540 [1710.09465]

  65. [65]

    The bispectrum of galaxies from high-redshift galaxy surveys: primordial non-Gaussianity and non-linear galaxy bias

    E. Sefusatti and E. Komatsu,Bispectrum of galaxies from high-redshift galaxy surveys: Primordial non-Gaussianity and nonlinear galaxy bias, Phys. Rev. D76 (2007) 083004 [0705.0343]

  66. [66]

    Moradinezhad Dizgah, M

    A. Moradinezhad Dizgah, M. Biagetti, E. Sefusatti, V. Desjacques and J. Noreña,Primordial non-Gaussianity from biased tracers: likelihood analysis of real-space power spectrum and bispectrum, J. Cosmology Astropart. Phys.2021 (2021) 015 [2010.14523]

  67. [67]

    D’Amico, M

    G. D’Amico, M. Lewandowski, L. Senatore and P. Zhang,Limits on primordial non-Gaussianities from BOSS galaxy-clustering data, Phys. Rev. D111 (2025) 063514 [2201.11518]. – 33 –

  68. [68]

    Cabass, M.M

    G. Cabass, M.M. Ivanov, O.H.E. Philcox, M. Simonović and M. Zaldarriaga,Constraints on multifield inflation from the BOSS galaxy survey, Phys. Rev. D106 (2022) 043506 [2204.01781]

  69. [69]

    Harscouet, J.A

    L. Harscouet, J.A. Cowell, J. Ereza, D. Alonso, H. Camacho, A. Nicola et al.,Fast Projected Bispectra: the filter-square approach, The Open Journal of Astrophysics8 (2025) 6 [2409.07980]

  70. [70]

    Alonso, J

    D. Alonso, J. Sanchez and A. Slosar,Namaster: A unified pseudo-cℓ framework, Monthly Notices of the Royal Astronomical Society484 (2019) 4127

  71. [71]

    Chisari, D

    N.E. Chisari, D. Alonso, E. Krause et al.,Core cosmology library: Precision cosmological predictions for lsst, The Astrophysical Journal Supplement Series242 (2019) 2

  72. [72]

    Virtanen, R

    P. Virtanen, R. Gommers, T.E. Oliphant et al.,Scipy 1.0: Fundamental algorithms for scientific computing in python, Nature Methods 17 (2020) 261

  73. [73]

    Harris, K.J

    C.R. Harris, K.J. Millman, S.J. van der Walt et al.,Array programming with NumPy, Nature 585 (2020) 357

  74. [74]

    Hunter,Matplotlib: A 2d graphics environment, Computing in Science & Engineering9 (2007) 90

    J.D. Hunter,Matplotlib: A 2d graphics environment, Computing in Science & Engineering9 (2007) 90

  75. [75]

    Foreman-Mackey, D.W

    D. Foreman-Mackey, D.W. Hogg, D. Lang and J. Goodman,emcee: The mcmc hammer, Publications of the Astronomical Society of the Pacific125 (2013) 306

  76. [76]

    Zonca, L

    A. Zonca, L. Singer, D. Lenz, M. Reinecke, C. Rosset, E. Hivon et al.,healpy: equal area pixelization and spherical harmonics transforms for data on the sphere in python, Journal of Open Source Software4 (2019) 1298

  77. [77]

    Vanneste, J.-C

    S. Vanneste, J.-C. Hamilton and J. Grain,xqml: A pure python implementation of the qml method, 2021

  78. [78]

    Aghanim, Y

    Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi et al., Planck 2018 results. VIII. Gravitational lensing, A&A 641 (2020) A8 [1807.06210]

  79. [79]

    Gratton and A

    S. Gratton and A. Challinor,Understanding parameter differences between analyses employing nested data subsets, MNRAS 499 (2020) 3410 [1911.07754]. – 34 – Analysis setting fNL A1 χ2/Ndof PTE Fiducial −48.6+25.5 −28.5 1.02 ± 0.04 54.3/55 0.50 p = 1.6 −71.3+31.7 −38.0 1.03 ± 0.04 50.1/55 0.66 No Cgg ℓ −56.0+26.5 −29.7 1.00 ± 0.05 24.6/30 0.74 No Cgκ ℓ −45.8...