Minkowski Tensors are tensorial generalizations of the scalar Minkowski Functionals. Due to their tensorial nature they contain additional morphological information of structures, in particular about shape and alignment, in comparison to the scalar Minkowski functionals. They have recently been used [29] to study the statistical isotropy of temperature and E mode data from the Planck satellite. The calculation in [29] relied on stereographic projection of the fields to extract the shape and alignment information. In this work, we calculate Minkowski Tensors directly on the sphere and compute the net alignment in the data, based on a recent work that extends the definition of Minkowski Tensors to random fields on curved spaces. This method circumvents numerical errors that can be introduced by the stereographic projection. We compare the resulting net alignment parameter values obtained from the frequency coadded CMB temperature data cleaned by the SMICA pipeline, to those obtained from simulations that include instrumental beam effects and residual foreground and noise. We find very good agreement between the two within $\approx$ 1$\sigma$ . We further compare the alignments obtained from the beam-convolved CMB maps at individual Planck frequencies to those in the corresponding simulations. We find no significant difference between observed data and simulations across all Planck frequencies, except for the 30 GHz channel. For the 30 GHz channel we find $\approx$ 2$\sigma$ difference between the data and the simulations. This mild disagreement most likely originates from inaccurate estimation of the instrumental beam at 30 GHz.

We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimensional field, and reviewing the definition of Minkowski Tensors. Focusing on the translational invariant statistic $W^{1,1}_{2}$ - a $2 \times 2$ matrix - we calculate its eigenvalues for both the entire excursion set ($\Lambda_{1},\Lambda_{2}$) and for individual connected regions and holes within the set ($\lambda_{1},\lambda_{2}$). The ratio of eigenvalues $\Lambda_{2}/\Lambda_{1}$ informs us of the presence of global anisotropies in the data, and $\langle \lambda_{2}/\lambda_{1} \rangle$ is a measure of the mean shape of peaks and troughs in the density field. We study these quantities for a Gaussian field, then consider how they are modified by the effect of gravitational collapse using the latest Horizon Run 4 cosmological simulation. We find $\Lambda_{1,2}$ are essentially independent of gravitational collapse, as the process maintains statistical isotropy. However, the mean shape of peaks is modified significantly - overdensities become relatively more circular compared to underdensities of the same area. When applying the statistic to a redshift space distorted density field, we find a significant signal in the eigenvalues $\Lambda_{1,2}$, suggesting that they can be used to probe the large-scale velocity field.

We generalize the translation invariant tensor-valued Minkowski Functionals which are defined on two-dimensional flat space to the unit sphere. We apply them to level sets of random fields. The contours enclosing boundaries of level sets of random fields give a spatial distribution of random smooth closed curves. We obtain analytic expressions for the ensemble expectation values for the matrix elements of the tensor-valued Minkowski Functionals for isotropic Gaussian and Rayleigh fields. We elucidate the way in which the elements of the tensor Minkowski Functionals encode information about the nature and statistical isotropy (or departure from isotropy) of the field. We then implement our method to compute the tensor-valued Minkowski Functionals numerically and demonstrate how they encode statistical anisotropy and departure from Gaussianity by applying the method to maps of the Galactic foreground emissions from the PLANCK data.

CMB polarization data is usually analyzed using $E$ and $B$ modes because they are scalars quantities under rotations along the lines of sight and have distinct physical origins. We explore the possibility of using the Stokes parameters $Q$ and $U$ for complementary analysis and consistency checks in the context of searches for non-Gaussianity. We show that the Minkowski Functionals (MFs) of $Q,U$ are invariant under local rotations along the lines of sight even though $Q,U$ are spin-2 variables, for full sky analysis. The invariance does not hold for incomplete sky. For local type primordial non-Gaussianity, when we compare the non-Gaussian deviations of MFs for $Q,U$ to what is obtained for $E$ mode or temperature fluctuations, we find that the amplitude is about an order of magnitude lower and the shapes of the deviations are different. This finding can be useful in distinguishing local type non-Gaussianity from other origins of non-Gaussianity in the observed data. Lastly, we analyze the sensitivity of the amplitudes of the MFs for $Q$, $U$ and the number density of singularities of the total polarization intensity to the tensor-to-scalar ratio, $r$, and find that all of them decrease as $r$ increases.

Tensor Minkowski Functionals (TMFs) are tensorial generalizations of the usual Minkowski Functionals which are scalar quantities. We introduce them here for use in cosmological analysis, in particular to analyze CMB maps. They encapsulate information about the shapes and the orientation of structures. We focus on one of the TMFs, namely $W_2^{1,1}$, which is the generalization of the genus. The ratio of the eigenvalues of the average of $W_2^{1,1}$ over all structures, $\alpha$, encodes the net orientation; and the average of the ratios of the eigenvalues of $W_2^{1,1}$ for each structure, $\beta$, encodes the net anisotropy. We have developed a code that computes $W_2^{1,1}$, and from it $\alpha$ and $\beta$, for a set of structures on the plane. We compute $\alpha$ and $\beta$ as functions of threshold levels for simulated Gaussian and isotropic CMB fields. We obtain $\alpha$ to be one for both temperature and $E$ mode, which means that we recover the input statistical isotropy of density fluctuations in the simulations. The level of net anisotropy of hotspots and coldspots in the CMB fields is quantified by $\beta\sim 0.62$. Then we compute $\alpha$ and $\beta$ for data from PLANCK. We find that the temperature field agrees with the standard LCDM prediction of no net orientation within $3-\sigma$. However, we find that $E$ mode data shows a net orientation that deviates from the theoretical expectation at $14-\sigma$. The possible origin of this deviation may be instrumental effects or other sources and needs to be investigated further. For the net anisotropy we obtain values of $\beta$ for both temperature and $E$ mode that are consistent with the expectations from the standard LCDM simulations. Accurate measurements of $\alpha$ and $\beta$ can be used to test the standard model of cosmology and to search for deviations from it.

We study the signatures of local type primordial non-Gaussianity, parametrized by $\fnl$, of scalar perturbations in CMB polarization using the probability distribution functions, Minkowski Functionals and Betti numbers. We show that the lowest order non-Gaussian deviation of the PDF of the total polarization intensity is at order $(\fnl\sigma)^2$. We calculate the non-Gaussian deviations of Minkowski Functionals and Betti numbers from simulated polarization maps. We find that $E$ mode polarization provides independent and equally strong constraint on $\fnl$ as temperature fluctuations. The non-Gaussian signal in the total polarization intensity, however, is much weaker and has a relatively large cosmic variance and hence may not be useful for detecting local type non-Gaussianity.