We apply the Minkowski tensor statistics to three dimensional Gaussian random fields. Minkowski tensors contain information regarding the orientation and shape of excursion sets, that is not present in the scalar Minkowski functionals. They can be used to quantify globally preferred directions, and additionally provide information on the mean shape of subsets of a field. This makes them ideal statistics to measure the anisotropic signal generated by redshift space distortion in the low redshift matter density field. We review the definition of the Minkowski tensor statistics in three dimensions, focusing on two coordinate invariant quantities $W^{0,2}_{1}$ and $W^{0,2}_{2}$. We calculate the ensemble average of these $3 \times 3$ matrices for an isotropic Gaussian random field, finding that they are proportional to products of the identity matrix and a corresponding scalar Minkowski functional. We show how to numerically reconstruct $W^{0,2}_{1}$ and $W^{0,2}_{2}$ from discretely sampled fields and apply our algorithm to isotropic Gaussian fields generated from a linear $\Lambda$CDM matter power spectrum. We then introduce anisotropy by applying a linear redshift space distortion operator to the matter density field, and find that both $W^{0,2}_{1}$ and $W^{0,2}_{2}$ exhibit a distinct signal characterised by inequality between their diagonal components. We discuss the physical origin of this signal and how it can be used to constrain the redshift space distortion parameter $\Upsilon \equiv f/b$.

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.

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.

We introduce the numbers of hot and cold spots, $n_h$ and $n_c$, of excursion sets of the CMB temperature anisotropy maps as statistical observables that can discriminate different non-Gaussian models. We numerically compute them from simulations of non-Gaussian CMB temperature fluctuation maps. The first kind of non-Gaussian model we study is the local type primordial non-Gaussianity. The second kind of models have some specific form of the probability distribution function from which the temperature fluctuation value at each pixel is drawn, obtained using HEALPIX. We find the characteristic non-Gaussian deviation shapes of $n_h$ and $n_c$, which is distinct for each of the models under consideration. We further demonstrate that $n_h$ and $n_c$ carry additional information compared to the genus, which is just their linear combination, making them valuable additions to the Minkowski Functionals in constraining non-Gaussianity.