We propose several pseudoentropy measures that agree well with the Wehrl entropy, but are significantly faster to compute. All of them are rotationally invariant measures of entanglement very sensitive to non-Gaussianity, anisotropy, and statistical dependence of spherical harmonic coefficients. We provide a simple proof that the projection pseudoentropy converges to the Wehrl entropy with increasing dimensionality of the ancilla projection space. Furthermore, for $l=2$, we show that both the Wehrl entropy and the angular pseudoentropy can be expressed as functions of the squared chordal distance of multipole vectors. We also show that the angular pseudoentropy can distinguish between Gaussian and non-Gaussian temperature fluctuations at large multipoles and henceforth provides a non-brute-force method for identifying non-Gaussianities. This allows us to study possible hints of statistical anisotropy and non-Gaussianity in the CMB up to multipole $l=1000$ using Planck 2015/2018, and WMAP 7-yr data. We find that $l=5$ and $28$ have a large entropy at $2$--$3\sigma$ significance and a slight hint towards a connection of this with the cosmic dipole. On a wider range of large angular scales we do not find indications of isotropy/Gaussianity violation. We also find a small-scale range, $l\in[895,905]$, that is incompatible with the assumptions at about $3\sigma$ level, although how much this significance can be reduced by taking into account the selection effect is left as an open question. We find overall similar results in our analysis of the 2015 and the 2018 data. Finally, we also demonstrate how a range of angular momenta can be studied with the range angular pseudoentropy. Our main purpose is to introduce the methods, analyze their mathematical background, and demonstrate their usage for providing researchers in this field with an additional tool.

Multipole vectors and pseudoentropies provide powerful tools for a numerically fast and vivid investigation of possible statistically anisotropic, respectively non-Gaussian signs in CMB temperature fluctuations. After reviewing and linking these two conceptions we compare their application to data analysis using the Planck 2015 NILC full sky map.

The statistical cosmological principle states that observables on the celestial sphere are sampled from a rotationally invariant distribution. Previously certain large scale anomalies which violate this principle have been found, for example an alignment of the lowest multipoles with the cosmic dipole direction. In this work we continue the search for possible anomalies using multipole vectors which represent a convenient tool for this purpose. In order to study the statistical behavior of multipole vectors, we revisit several construction methods. We investigate all four full-sky foreground-cleaned maps from the Planck 2015 release with respect to four meaningful physical directions using computationally cheap statistics that have a simple geometric interpretation. We find that the full-sky SEVEM map deviates from all the other cleaned maps, as it shows a strong correlation with the Galactic Pole and Galactic Center. The other three maps COMMANDER, NILC and SMICA show a consistent behavior. On the largest angular scales, $l \leq 5$, as well as on intermediate scales, $l = 20, 21, 22, 23, 24$, all of them are unusually correlated with the cosmic dipole direction. These scales coincide with the scales on which the angular power spectrum deviates from the Planck 2015 best-fit {\Lambda}CDM model. In the range $2 \leq l \leq 50$ as a whole there is no unusual behavior visible globally. We do not find abnormal intramultipole correlation, i.e. correlation of multipole vectors inside a given multipole without reference to any outer direction.