Minkowski Functionals of CMB polarisation intensity with Pynkowski: theory and application to Planck data

The analysis of the Cosmic Microwave Backgound (CMB) has yielded very precise information about the composition and evolution of the Universe, especially thanks to the study of its Angular Power Spectrum. However, this quantity is blind to the presence of non--Gaussianities and deviations from statistical isotropy; the study of these effects can be performed with other statistics such as Minkowski Functionals (MFs). These tools have been largely applied to the CMB temperature data without any significant detection of a deviation from Gaussianity and isotropy for this field. In this work, we extend the formalism of MFs to the squared CMB polarisation intensity, $P^2=Q^2+U^2$. We use the Gaussian Kinematic Formula to analytically determine the theoretical predictions of MFs for a Gaussian isotropic field. We then develop a software which computes these quantities on $P^2$ Healpix maps and apply it to simulations in order to verify the robustness of both theory and methodology. Finally, we compute the MFs of Planck $P^2$ maps and compare them with the ones from realistic simulations which include CMB and systematics. We find no significant deviations from primordial Gaussianity or isotropy in Planck CMB polarisation data. However, MFs could play an important role in the analysis of statistical properties of polarised CMB signal from upcoming observations with improved sensitivity. We publicly release the software to compute MFs in arbitrary scalar Healpix maps as a Python package called $\texttt{Pynkowski}$ (https://github.com/javicarron/pynkowski). It is modular, fully documented, and designed to be easy to use in different applications.