Compressed Gaussian likelihood for the Planck low-$\ell$ data

Abstract

We present a compressed Gaussian likelihood for the Planck CMB low-$\ell$ E-mode polarization data, constructed from the SRoll2 likelihood which provides the tightest constraint on the reionization optical depth $\ddot{I}„$ to date. The non-Gaussian form of CMB low-$\ell$ TT and EE likelihoods makes them incompatible with Fisher matrix analyses that require an analytic Gaussian $\ddot{I}{‡}^2$, such as the Fisher-bias formalism and Fisher forecasts. We show that the $\ddot{I}{‡}^2$ of an offset log-normal likelihood takes a Gaussian form in the log-transformed power spectrum amplitudes, and can therefore serve as a proxy for the true Gaussian likelihood of this variable in Fisher matrix analyses, without any explicit change of variables. Building on this, we compress the SRoll2 likelihood into a small number of piecewise offset log-normal functions and validate it against the full SRoll2 likelihood via MCMC combined with Planck and ACT DR6 data, finding excellent agreement across all $\hat{I}›$CDM parameters and in extended cosmological models. We further demonstrate that Fisher matrix uncertainty estimates from our compressed likelihood agree well with the full MCMC posteriors. We release our compressed likelihood planck-gaussian-lowl, a lightweight Python package incorporating the compressed low-$\ell$ TT likelihood from previous work, allowing a straightforward incorporation of the Planck CMB low-$\ell$ data into any Gaussian-likelihood-based analysis. The package is publicly available at \href{https://github.com/nanoomlee/planck-gaussian-lowl}{github.com/nanoomlee/planck-gaussian-lowl}.

7+2 pages, 4 figures. Comments are welcome