Despite the wealth of $Planck$ results, there are difficulties in disentangling the primordial non-Gaussianity of the Cosmic Microwave Background (CMB) from the secondary and the foreground non-Gaussianity (NG). For each of these forms of NG the lack of complete data introduces model-dependencies. Aiming at detecting the NGs of the CMB temperature anisotropy $\delta T$, while paying particular attention to a model-independent quantification of NGs, our analysis is based upon statistical and morphological univariate descriptors, respectively: the probability density function $P(\delta T)$, related to ${\mathrm v}_{0}$, the first Minkowski Functional (MF), and the two other MFs, ${\mathrm v}_{1}$ and ${\mathrm v}_{2}$. From their analytical Gaussian predictions we build the discrepancy functions $\Delta_{k}$ ($k=P,0,1,2$) which are applied to an ensemble of $10^{5}$ CMB realization maps of the $\Lambda$CDM model and to the $Planck$ CMB maps. In our analysis we use general Hermite expansions of the $\Delta_{k}$ up to the $12^{th}$ order, where the coefficients are explicitly given in terms of cumulants. Assuming hierarchical ordering of the cumulants, we obtain the perturbative expansions generalizing the $2^{nd}$ order expansions of Matsubara to arbitrary order in the standard deviation $\sigma_0$ for $P(\delta T)$ and ${\mathrm v}_0$, where the perturbative expansion coefficients are explicitly given in terms of complete Bell polynomials. The comparison of the Hermite expansions and the perturbative expansions is performed for the $\Lambda$CDM map sample and the $Planck$ data. We confirm the weak level of non-Gaussianity ($1$-$2$)$\sigma$ of the foreground corrected masked $Planck$ $2015$ maps.

The Friedmann-Lemaitre-Robertson-Walker (FLRW) metric assumes comoving spatial rigidity of metrical properties. The curvature term in comoving coordinates is environment-independent and cannot evolve. In the standard model, structure formation is interpreted accordingly: structures average out on the chosen metrical background, which remains rigid in comoving coordinates despite nonlinear structure growth. The latter claim needs to be tested, since it is a hypothesis that is not derived using general relativity. We introduce a test of the comoving rigidity assumption by measuring the two-point auto-correlation function on comoving scales---assuming FLRW comoving spatial rigidity---in order to detect shifts in the baryon acoustic oscillation (BAO) peak location for Large Red Galaxy (LRG) pairs of the Sloan Digital Sky Survey Data Release 7. In tangential directions, subsets of pairs overlapping with superclusters or voids show the BAO peak. The tangential BAO peak location for overlap with Nadathur & Hotchkiss superclusters is $4.3\pm1.6$ Mpc/h less than that for LRG pairs unselected for supercluster overlap, and $6.6\pm2.8$ Mpc/h less than that of the complementary pairs. Liivamagi et al. superclusters give corresponding differences of $3.7\pm2.9$ Mpc/h and $6.3\pm2.6$ Mpc/h, respectively. We have found moderately significant evidence (Kolmogorov--Smirnov tests suggest very significant evidence) that the BAO peak location for supercluster-overlapping pairs is compressed by about 6% compared to that of the complementary sample, providing a potential challenge to FLRW models and a benchmark for predictions from models based on an averaging approach that leaves the spatial metric a priori unspecified.

The 3-torus (T^3) FLRW model better fits the nearly zero large-scale auto-correlation of the WMAP CMB sky maps than the infinite flat model. The T^3 model's WMAP parameters imply approximately equal-redshift topological lensing at z \sim 6. We investigate observational strategies for rejecting the T^3 solution or providing candidate topologically lensed galaxy pairs. T^3 holonomies are applied to (i) existing z \sim 6 observations and (ii) simulated observations, creating multiply connected catalogues. Corresponding simply connected catalogues are generated. Each catalogue is analysed using a successive filter method and collecting matched quadruples. Quadruple statistics between the multiply and simply connected catalogues are compared. The expected rejection of the hypothesis, or detection of candidate topologically lensed galaxies, is possible at a significance of 5% for a pair of T^3 axis-centred northern and southern surveys if photometric redshift accuracy is \sigma(\zphot) < 0.01 for a pair of nearly complete 100 deg^2 surveys with a total of > 500 galaxies over 4.3 < z < 6.6, or for a pair of 196 deg^2 surveys with > 400 galaxies and \sigma(\zphot) < 0.02 over 4<z<7. Dropping the maximum time interval in a pair from \Delta t =1 Gyr/h to \Delta t =0.1 Gyr/h yields a requirement of \sigma(\zphot) < 0.005 or \sigma(\zphot) < 0.01, respectively. Millions of z \sim 6 galaxies will be observed over fields of these sizes during the coming decades, implying much stronger constraints. The question is not if the hypothesis will be rejected or confirmed, it is when.