Pith. sign in

REVIEW 2 major objections 3 minor 19 references

Position-space curved-sky anisotropy quadratic estimation

T0 review · 2 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Optimal curved-sky quadratic estimators for CMB anisotropies follow from one covariance-response kernel and one weight formula.

desk verdict Useful methods supplement, but the printed optimal-weight recipe in Eq. (4.4) has a spin-index sign error that makes it unusable as-is. read the letter →

arxiv 1908.02016 v1 pith:ZYLE2MJG submitted 2019-08-06 astro-ph.CO

Julien Carron This is my paper · ORCID
classification astro-ph.CO
keywords cosmicmicrowavebackgroundquadraticestimatorsCMBlensingstatisticalanisotropyspin-weightcorrelationfunctionsgradientandcurlmodescurvedskyWignersmall-dtransforms
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper establishes a general recipe for building optimal quadratic estimators of statistical anisotropies in the cosmic microwave background on the curved sky. Its central claim is that, once the data covariance response $W$ of Eq. (3.3) is known, the optimal joint gradient and curl estimator is a separable product of two position-space, spin-weighted filtered maps with explicit weights, and that responses and Gaussian noise biases follow from compact one-dimensional Wigner $d$-function integrals. This matters because lensing, patchy reionization, cosmic birefringence, point sources, and noise inhomogeneities are all anisotropy probes that can be treated as special cases of the same construction, with no separate derivations needed. The document also shows that the standard temperature, polarization, and minimum-variance lensing estimators are recovered as gradient-mode special cases.

What carries the argument

The load-bearing object is the covariance response kernel $W^{a,st}_{\ell}$ of Eq. (3.3), which encodes how the data covariance changes under a spin-$r$ anisotropy source. The estimator itself is the separable quadratic form of Eq. (1.7), a product of two spin-weighted filtered maps, whose weights are set by Eq. (4.4) as the first Newton-Raphson step of the Gaussian likelihood evaluated at zero anisotropy. Responses and noise covariances are then computed by projecting these position-space products onto gradient and curl modes, which reduces every integral to a one-dimensional Wigner $d$-function transform. This combination of a position-space correlation-function representation and the covariance response is what turns each new anisotropy source into a known optimal estimator.

What would settle it

Simulate Gaussian CMB skies containing a known anisotropy whose covariance response is $W$, apply the optimal weights of Eq. (4.4), and compare the estimator's variance to the inverse of the corresponding Fisher matrix; any significant excess variance would show that the analytic optimality argument misses something. On the noise side, averaging the quadratic estimator over isotropic simulations should reproduce the analytic $N^{(0)}_{\ell}$ of Eqs. (2.4)-(2.5) to numerical precision, so a mismatch would falsify the Gaussian covariance calculations.

Watch

Extended reading notes

Core claim

The core claim is that for a Gaussian CMB whose covariance responds to a spin-$r$ anisotropy through the separable kernel $W$ of Eq. (3.3), the minimum-variance quadratic estimator of the anisotropy's gradient and curl modes is given by Eq. (1.7) with the weights of Eq. (4.4): one leg of the estimator is a delta function in spin, the other leg is proportional to $W$, and both legs act on inverse-variance-filtered maps. The same formalism supplies analytic responses, Eqs. (3.5)-(3.7), and Gaussian noise biases between arbitrary pairs of estimators, Eqs. (2.4)-(2.5), all reduced to one-dimensional Wigner small-$d$ integrals. In the lensing case the construction reproduces the known full-sky lensing estimators, including temperature, polarization, and minimum-variance versions, and extends them to curl modes and to anisotropy sources of arbitrary spin.

Load-bearing premise

The optimality argument assumes that the anisotropy enters only through the separable covariance response $W$ of Eq. (3.3), that the CMB remains Gaussian under the anisotropy, and that the Newton-Raphson step is taken at zero anisotropy; if any of these fails, the proposed weights are approximate rather than optimal.

Editorial extensions

If this is right

  • Any anisotropy whose covariance response has the separable form of Eq. (3.3) immediately gets an optimal joint gradient and curl estimator: supplying $W$ fixes all weights, responses, and leading Gaussian noise biases.
  • The analytic responses of Eqs. (3.5)-(3.7) allow gradient and curl estimates to be normalized independently, so one pipeline can produce unbiased maps for lensing, birefringence, or patchy reionization.
  • The Gaussian noise-bias formulae of Eqs. (2.4)-(2.5) provide the leading $N^{(0)}_{\ell}$ terms between arbitrary pairs of estimators, the quantity needed to debias anisotropy power spectra and cross-spectra.
  • The known temperature, polarization, and minimum-variance lensing estimators are recovered as gradient-mode special cases of the single weight rule of Eq. (4.4), confirming the construction against established results.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct test of the optimality claim would be to feed a numerical covariance response $W$ measured from simulated anisotropic skies, compare the resulting estimator variance with the inverse Fisher matrix, and see whether any excess variance reveals terms beyond Eq. (3.3).
  • The same Newton-Raphson logic could be applied to non-separable covariance responses numerically, which would quantify how much optimality is lost when the separable form is violated.
  • Because the method only assumes Gaussianity under the anisotropy, the estimator construction could be reused for foregrounds or secondary anisotropies whose covariance response is known, not just the CMB examples listed in the paper.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. The paper presents a position-space, spin-weight formalism for quadratic estimators (QEs) of statistical anisotropy of the CMB on the curved sky. It derives analytic expressions for Gaussian noise biases (Section II), estimator responses (Section III), and optimal joint gradient/curl weights from a general covariance response W (Section IV). The document is a supplement to the public Planck 2018 lensing pipeline and applies the formalism to lensing, modulation, patchy reionization, polarization rotation, point sources, and noise inhomogeneities.

Significance. If the formulas are correct, the paper is a useful reference that unifies several anisotropy estimators into one separable QE framework and reduces response and noise computations to one-dimensional Wigner small-d transforms. It also makes concrete contact with known Okamoto-Hu lensing estimators. However, the printed weight formula in Section IV contains a spin-index inconsistency that prevents the central recipe from being used as written; this significantly limits the immediate utility of the paper until corrected.

major comments (2)
  1. [Section IV, Eqs. (4.3)-(4.4)] The printed weight formula is internally inconsistent with the expression it is meant to encode. In Eq. (4.3), the second leg has the spherical harmonic s+r Y_lm, so in the notation of Eq. (1.7) its output spin is t_o = s+r, and the factor is 2 W^{-r,st}_l / t_2. Equation (4.4), however, gives w^{-s+r,t}_l = 2/t_2 W^{-r,-st}_l, i.e. output spin -s+r and a second index -st on W. Since the first leg has output spin -s, the total output spin would be -2s+r rather than r. For r=0, s=2, this would produce a spin -4 object, contradicting the scalar modulation and polarization-rotation examples in Section V. Unless there is an unstated index convention, Eq. (4.4) must be corrected to w^{s+r,t}_l = 2/t_2 W^{-r,st}_l (or the equivalent with clearly defined index symmetries), and the derivation should be checked against the examples.
  2. [Section IV, between Eqs. (4.2) and (4.3)] The transition from the functional derivative in Eq. (4.2) to the explicit weight expression in Eq. (4.3) is not shown; the text states only that performing the derivative using representation (3.3) gives the result. Because the resulting weight formula Eq. (4.4) contains a sign/spin inconsistency, the omitted algebra is not merely a presentation issue. The authors should provide the intermediate spin-weight manipulation for at least one case, or otherwise verify that the printed equations satisfy the spin-sum rule s_o + t_o = r.
minor comments (3)
  1. [Section IV, Eq. (4.4)] The symbols s2 and t2 are defined in the text as "1(s=0) or 2(s≠0)", but this notation is easy to misread as s^2 and t^2; using s_2 and t_2 would be clearer.
  2. [Section I.B, Eq. (1.7)] The same symbol w is used for both legs even though the weights may differ; introducing separate symbols, e.g. w^{(1)} and w^{(2)}, would reduce ambiguity in equations such as Eq. (4.4).
  3. [Sections II and III] The paper states that most formulas follow by applying Eq. (1.10), but no representative derivation is shown. Since this is a methods supplement, adding a short worked derivation of one of the N(0) or response formulas would greatly help readers verify the sign conventions.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the estimator responses and optimal weights are derived from a forward covariance model and a likelihood expansion, not fitted to or renamed from the predicted quantities.

full rationale

The paper is a methods supplement whose central objects are derived rather than fitted. The estimator definition in Eq. (1.7) is generic, and the response calculation in Section III B takes the covariance response W from Eq. (3.3) as an input and derives analytic estimator responses in Eqs. (3.5)-(3.7). The optimal weights in Section IV are obtained from a Newton-Raphson/log-likelihood gradient derivation, explicitly citing the external framework of Hanson and Lewis [10], and the result is benchmarked against the externally established Okamoto-Hu estimators [11]. No fitted parameter is relabeled as a prediction, and no load-bearing step invokes a self-citation as its justification. The only self-reference is to the author's public pipeline repository, which is code rather than a result baked into the derivation. The apparent spin-index inconsistency in Eq. (4.4) noted by a skeptical reader is a possible typographical or convention issue, but it is not a circularity: it concerns internal consistency of a printed formula, not a claim whose output is equivalent to its input. The derivation is self-contained against external benchmarks (Okamoto-Hu lensing estimators), so the appropriate circularity score is 0.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters or invented entities. The paper's results rest on standard angular momentum algebra and standard CMB assumptions (Gaussianity, isotropic noise for analytic N0, separable response models).

assumptions (4)
  • standard math Standard spin-weighted spherical harmonic and Wigner d-matrix identities, including Eq. (1.10), hold.
    Used to derive all response and noise formulas in Sections II and III.
  • domain assumption Inverse-variance filtered CMB maps are statistically isotropic for the analytic N(0) estimates.
    Section II states the analytic N(0) assume isotropy of the filtered maps; this is an approximation, with RDN0 as a more accurate alternative.
  • domain assumption The anisotropy response is fully captured by the separable covariance perturbation of Eq. (3.3) with spin-r weights W^{a,st}_l.
    This parametrization underlies the weight derivation and the examples in Section V.
  • domain assumption Under the anisotropy, the CMB remains Gaussian; the leading Newton-Raphson step gives near-optimal estimates for small anisotropy.
    Section IV, first paragraph, explicitly assumes Gaussianity; lensing generates non-Gaussianity, making the optimality approximate.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Position-space curved-sky anisotropy quadratic estimation." pith.science (2026). https://pith.science/paper/ZYLE2MJG

@misc{pith2026190802016,
  author       = {Pith},
  title        = {Pith review of: Position-space curved-sky anisotropy quadratic estimation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZYLE2MJG}},
  note         = {Machine review of arXiv:1908.02016}
}
read the original abstract

This document supplements the release of the Planck 2018 CMB lensing pipeline, now made publicly available. It collects calculations relevant to curved-sky separable quadratic estimators in the spin-weight, position-space correlation function formalism, including analytic calculations of estimator responses and Gaussian noise biases between arbitrary pairs of quadratic estimators. It also contains the derivation of optimal, joint gradient and curl mode quadratic estimators for parametrized anisotropy of arbitrary spin.

Figures

Figures reproduced from arXiv: 1908.02016 by the authors.

Figure 1
Figure 1. FIG. 1. The geometry and angles in Eq. (1.5), with the local axes in green. Figure originally from Wayne Hu tutorials, [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 1 canonical work pages

  1. [1]

    Planck 2018 results. VIII. Gravitational lensing,

    N. Aghanim et al. (Planck), “Planck 2018 results. VIII. Gravitational lensing,” (2018), arXiv:1807.06210 [astro-ph.CO]

    arXiv 2018

  2. [2]

    B-mode CMB Polarization from Patchy Screening during Reionization,

    Cora Dvorkin, Wayne Hu, and Kendrick M. Smith, “B-mode CMB Polarization from Patchy Screening during Reionization,” Phys. Rev. D79, 107302 (2009), arXiv:0902.4413 [astro-ph.CO]

    arXiv 2009

  3. [3]

    Detection of Gravitational Lensing in the Cosmic Microwave Back- ground,

    Kendrick M. Smith, Oliver Zahn, and Olivier Dore, “Detection of Gravitational Lensing in the Cosmic Microwave Back- ground,” Phys. Rev. D76, 043510 (2007), arXiv:0705.3980 [astro-ph]

    arXiv 2007

  4. [4]

    HEALPix - A Framework for high resolution discretization, and fast analysis of data distributed on the sphere,

    K. M. Gorski, Eric Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, and M. Bartelman, “HEALPix - A Framework for high resolution discretization, and fast analysis of data distributed on the sphere,” Astrophys. J. 622, 759–771 (2005), arXiv:astro-ph/0409513 [astro-ph]

    arXiv 2005

  5. [5]

    Weak gravitational lensing of the cmb,

    Antony Lewis and Anthony Challinor, “Weak gravitational lensing of the cmb,” Phys. Rept. 429, 1–65 (2006), arXiv:astro- ph/0601594 [astro-ph]

    arXiv 2006

  6. [6]

    Fast estimation of polarization power spectra using correlation functions,

    Gayoung Chon, Anthony Challinor, Simon Prunet, Eric Hivon, and Istvan Szapudi, “Fast estimation of polarization power spectra using correlation functions,” Mon. Not. Roy. Astron. Soc. 350, 914 (2004), arXiv:astro-ph/0303414 [astro-ph]

    arXiv 2004

  7. [7]

    Lensed CMB power spectra from all-sky correlation functions,

    Anthony Challinor and Antony Lewis, “Lensed CMB power spectra from all-sky correlation functions,” Phys. Rev. D71, 103010 (2005), arXiv:astro-ph/0502425 [astro-ph]

    arXiv 2005

  8. [8]

    A Measurement of the Cosmic Microwave Background Gravitational Lensing Potential from 100 Square Degrees of SPTpol Data,

    K. T. Story et al. (SPT), “A Measurement of the Cosmic Microwave Background Gravitational Lensing Potential from 100 Square Degrees of SPTpol Data,” Astrophys. J. 810, 50 (2015), arXiv:1412.4760 [astro-ph.CO]

    arXiv 2015

Show all 19 references
  1. [9]

    The Atacama Cosmology Telescope: Two-Season ACTPol Lensing Power Spectrum,

    Blake D. Sherwin et al. (ACT), “The Atacama Cosmology Telescope: Two-Season ACTPol Lensing Power Spectrum,” Submitted to: Phys. Rev. D (2016), arXiv:1611.09753 [astro-ph.CO]

    2016 arXiv

  2. [10]

    Estimators for CMB Statistical Anisotropy,

    Duncan Hanson and Antony Lewis, “Estimators for CMB Statistical Anisotropy,” Phys. Rev. D80, 063004 (2009), arXiv:0908.0963 [astro-ph.CO]

    2009 arXiv

  3. [11]

    CMB lensing reconstruction on the full sky,

    Takemi Okamoto and Wayne Hu, “CMB lensing reconstruction on the full sky,” Phys. Rev.D67, 083002 (2003), arXiv:astro- ph/0301031 [astro-ph]

    2003

  4. [12]

    CMB temperature lensing power recon- struction,

    Duncan Hanson, Anthony Challinor, George Efstathiou, and Pawel Bielewicz, “CMB temperature lensing power recon- struction,” Phys. Rev. D83, 043005 (2011), arXiv:1008.4403 [astro-ph.CO]

    2011 arXiv

  5. [13]

    The shape of the CMB lensing bispectrum,

    Antony Lewis, Anthony Challinor, and Duncan Hanson, “The shape of the CMB lensing bispectrum,” JCAP 1103, 018 (2011), arXiv:1101.2234 [astro-ph.CO]

    2011 arXiv

  6. [14]

    CMB lensing reconstruction biases in cross-correlation with large-scale structure probes,

    Giulio Fabbian, Antony Lewis, and Dominic Beck, “CMB lensing reconstruction biases in cross-correlation with large-scale structure probes,” (2019), arXiv:1906.08760 [astro-ph.CO]

    2019 arXiv

  7. [15]

    Geometry of weak lensing of CMB polarization,

    Anthony Challinor and Gayoung Chon, “Geometry of weak lensing of CMB polarization,” Phys. Rev. D66, 127301 (2002), arXiv:astro-ph/0301064 [astro-ph]

    2002 arXiv

  8. [16]

    CMB lensing and primordial squeezed non-Gaussianity,

    Ruth Pearson, Antony Lewis, and Donough Regan, “CMB lensing and primordial squeezed non-Gaussianity,” JCAP 1203, 011 (2012), arXiv:1201.1010 [astro-ph.CO]. 8

    2012 arXiv

  9. [17]

    Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity,

    P. A. R. Ade et al. (Planck), “Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity,” Astron. Astrophys. 571, A24 (2014), arXiv:1303.5084 [astro-ph.CO]

    2014 arXiv

  10. [18]

    Reconstructing Patchy Reionization from the Cosmic Microwave Background,

    Cora Dvorkin and Kendrick M. Smith, “Reconstructing Patchy Reionization from the Cosmic Microwave Background,” Phys. Rev. D79, 043003 (2009), arXiv:0812.1566 [astro-ph]

    2009 arXiv

  11. [19]

    Extragalactic Foreground Contamination in Temperature-based CMB Lens Reconstruction,

    Stephen J. Osborne, Duncan Hanson, and Olivier Dor´ e, “Extragalactic Foreground Contamination in Temperature-based CMB Lens Reconstruction,” JCAP 1403, 024 (2014), arXiv:1310.7547 [astro-ph.CO]

    2014 arXiv

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.