Cross-spectra likelihood for robust $τ$ constraints from all satellite polarisation data

Abstract

The Thomson scattering optical depth to reionisation, $τ$, one of the six parameters of the $Λ$CDM model, is primarily constrained by the large-scale E-mode polarisation of the Cosmic Microwave Background (CMB). In this work, we present the E-mode Likelihood for Cross-Analysis (elica), a multi-frequency, harmonic-space likelihood that combines all currently available large-scale satellite polarisation data, namely the Planck LFI 70 GHz channel, the Planck HFI 100 and 143 GHz channels processed with the SRoll2 map-making algorithm, and the WMAP Ka, Q, and V bands. The likelihood is built on an extension of the Hamimeche-Lewis formalism to multi-field partial-sky observations. We validate the pipeline using 500 realistic simulations and find that retaining all cross-spectra and the WMAP-LFI auto-spectrum eliminates the significant bias present when all spectra are retained, while preserving comparable uncertainties in the recovered value of $τ$. From the low-$\ell$ E-mode power spectrum alone, we obtain $τ= 0.0575_{-0.0058}^{+0.0048}$ (68% CL). Combining elica with the Planck low-$\ell$ temperature likelihood and the CamSpec high-$\ell$ likelihood, we find $τ= 0.0581_{-0.0059}^{+0.0048}$ and $\ln(10^{10}A_{\mathrm{s}}) = 3.048_{-0.012}^{+0.011}$. Including ACT{} DR6 + Planck CMB lensing and DESI DR2 BAO measurements, we derive an upper bound on the total neutrino mass of $\sum m_ν< 0.069$ eV (95% CL). Our results, obtained through careful cross-validation of all available large-scale polarisation datasets, robustly confirm that the optical depth remains relatively low. This severely constrains the possibility of explaining, or even significantly reducing, the tension between DESI-BAO and CMB observations with a high value of $τ$. The elica likelihood is publicly available.

Figures (1)

• Figure 1: Comparison of data from different channels with the fiducial model. Auto-spectra are shown in the top panel and cross-spectra in the bottom panel. For reference, we also plot the mean spectra obtained from 500 simulations. The error bars represent the standard deviation derived from the same simulations. The fiducial model assumes $\tau = 0.06$ and $A_{\mathrm{s}} = 2.12 \times 10^{-9}$. The grey line corresponds to an alternative model with $\tau = 0.09$ and a rescaled amplitude $A_{\mathrm{s}} = 2.25 \times 10^{-9}$.