Pseudo-Cℓ MASTER Algorithm

• Pseudo-Cℓ is a technique that computes unbiased two-point power spectra by correcting for mode mixing caused by survey window functions.
• The method uses a coupling matrix derived from the window function's power spectrum and Wigner-3j symbols to deconvolve the pseudo-spectrum.
• Enhancements like iMASTER incorporate a fiducial theoretical model during bin inversion to eliminate first-order biases in high-precision cosmological analyses.

The pseudo-$C_\ell$ MASTER algorithm is a methodology for extracting unbiased estimates of two-point power spectra (and correlation functions) from cosmological maps measured on the sky, accounting for incomplete sky coverage, complex survey geometries, and variable noise. The approach addresses the linear mixing of power spectrum modes induced by the application of a window function (survey geometry, mask, and/or pixel weighting) and provides a practical, accurate solution for a wide variety of large-scale structure (LSS) and cosmic microwave background (CMB) analyses. Subsequent improvements, notably the $i$Master algorithm, generalize the methodology to achieve unbiasedness in the presence of complex windows, utilizing a fiducial theoretical model within the deconvolution procedure (Singh, 2021).

1. Fundamental Concepts and Definitions

The core object for pseudo-$C_\ell$ estimation is the window-weighted overdensity field: $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$ where $\delta(\mathbf n)$ is the mean-zero overdensity field defined on the sphere (or flat sky) and $w(\mathbf n)$ encodes the pixel-space survey mask and possibly the inverse-noise weighting.

The pseudo-power spectrum is computed via a spherical harmonic expansion: $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$ where $\delta^w_{\ell m}$ are the spherical harmonic coefficients of the weighted field.

The key insight is that the expectation value of $\widetilde{C}_\ell$ is not the true sky power spectrum $C_\ell$, but a linear convolution: $i$0 where $i$1—the "coupling matrix"—is computed solely from the survey window's power spectrum $i$2 (autocorrelation of $i$3) and angular-momentum recoupling (Wigner-3j symbols). This linear relation can be inverted to "unbias" the bandpowers by binning $i$4 and inverting the corresponding binned coupling matrix $i$5: $i$6

2. Derivation and Standard MASTER Implementation

The classical MASTER procedure begins with the explicit windowing of the overdensity field, transformation to harmonic space, and construction of the pseudo-$i$7. By the convolution theorem (on the sphere), this pixel-space multiplication introduces $i$8-space mode mixing governed by the matrix $i$9: $C_\ell$0 where $C_\ell$1 and $C_\ell$2 are the field spins ($C_\ell$3 for scalars, $C_\ell$4 for shear fields), and $C_\ell$5 is the pseudo-spectrum of the window.

Binning in $C_\ell$6 limits the size of $C_\ell$7, making its inversion numerically tractable.

3. Standard Approximations and Sources of Residual Bias

Standard MASTER implementations utilize further approximations:

• Constant-in-bin assumption: The true $C_\ell$8 is assumed nearly constant within each $C_\ell$9-bin.
• Analytic binning operators: Simple forms for binning (e.g., rectangles) are used, with $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$0 and $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$1 approximately satisfying $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$2.
• Neglect of off-diagonal covariance: Non-diagonal elements are often ignored in inversion.

These assumptions lead to a binned coupling: $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$3 However, in practice, the non-constancy of $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$4 and window-induced broad mode-coupling produce a generic residual bias: $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$5 where $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$6 shares the form of $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$7 but is weighted by the true spectrum $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$8. Critically, $$\delta^w(\mathbf{n}) = w(\mathbf{n})\,\delta(\mathbf{n}),$$9, enhancing the impact of small errors $\delta(\mathbf n)$0 at high $\delta(\mathbf n)$1, often yielding $\delta(\mathbf n)$2 biases consistent with the scale of precision targeted in modern surveys.

4. The $\delta(\mathbf n)$3Master Algorithm and Removal of First-Order Bias

The $\delta(\mathbf n)$4Master enhancement leverages a fiducial theoretical model $\delta(\mathbf n)$5 to define a more accurate binning inversion, removing the "constant-in-bin" approximation. The procedure is:

1. Construct a generalized binned coupling matrix:

$\delta(\mathbf n)$6

where $\delta(\mathbf n)$7 bins data and $\delta(\mathbf n)$8 is defined such that $\delta(\mathbf n)$9 for $w(\mathbf n)$0.

2. Invert $w(\mathbf n)$1 to yield unbiased bandpowers:

$w(\mathbf n)$2

By matching the binning operator inversion exactly to the (potentially complex) model within each bin, this approach eliminates first-order residual biases from window complexity. Iterative updates of $w(\mathbf n)$3 are possible if the cosmological posterior deviates significantly from the assumed theory.

5. Configuration Space Formalism and Two-Point Correlation Generalization

The methodology for treating power spectra with complex windows is mirrored in configuration space. Given a two-point function $w(\mathbf n)$4 measured over $w(\mathbf n)$5, the corresponding pseudo-$w(\mathbf n)$6 is: $w(\mathbf n)$7 with $w(\mathbf n)$8 as the Wigner-small-$w(\mathbf n)$9 kernel.

For spin-2 fields (e.g., weak lensing shear), $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$0 correlators are used: $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$1 Pseudo-$$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$2 are formed from $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$3 and separated into $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$4 and $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$5 estimators; window convolution and E/B leakage are treated analogously. Imposing sharp cuts on $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$6 (e.g., restricting to $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$7) convolves the correlation function with an explicitly computable configuration-space kernel.

6. Practical Considerations and Error Propagation

Rigorous practical performance of the (i)MASTER approach requires:

• Window function $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$8 estimation accurate to $$\widetilde{C}_\ell = \frac{1}{2\ell+1} \sum_{m=-\ell}^{+\ell} |\delta^w_{\ell m}|^2,$$9, since the convolution matrix mixes power out to these scales.
• Understanding window errors, which break down into:
  • Additive errors: Spurious contributions to the measured $\delta^w_{\ell m}$0 uncorrelated with cosmological signal, scaling as $\delta^w_{\ell m}$1.
  • Multiplicative errors: Fluctuations of the form $\delta^w_{\ell m}$2, shifting $\delta^w_{\ell m}$3 proportionally to $\delta^w_{\ell m}$4.
• Application across a range of observables, including photometric galaxy maps, spectroscopic samples, intensity-mapping data, weak lensing, and CMB fields.
• For high-precision surveys (Stage-IV), window uncertainties must be propagated analytically, e.g., via marginalization over $\delta^w_{\ell m}$5 modes using priors of the Bridle–Hobson class.

A plausible implication is that, as power spectrum analyses attain sub-percent precision, robust treatment of window and mask uncertainties via approaches such as (i)MASTER becomes non-negotiable for unbiased cosmological inference.

7. Significance, Applicability, and Outlook

The pseudo-$\delta^w_{\ell m}$6 MASTER and $\delta^w_{\ell m}$7Master algorithms define a gold standard for unbiased power spectrum estimation in the presence of realistic survey masks and noise in cosmological analyses. The formalism is universally applicable to any survey geometry or window, extends naturally to both harmonic and configuration space, enables accurate E/B mode separation, and maintains computational tractability (Singh, 2021). The iMaster algorithm, in particular, restores unbiasedness even for highly structured masks or narrow bins at minimal additional computational complexity.

Addressing subtleties—including the propagation of window uncertainties and correcting higher-order bias, especially at high multipole—remains an active area, especially as next-generation surveys demand even greater accuracy in the recovery of cosmological parameters.