We examine the prospects of the high redshift Long Gamma Ray Bursts (LGRB) Hubble Diagram as a test of the basic cosmological principles. Analysis of the Hubble Diagram allows us to test several fundamental cosmological principles using the directly observed flux-distance-redshift relation. Modern LGRB data, together with the correlation between the spectral peak energy and the isotropic-equivalent radiated energy (the so-called Amati relation) can be used for construction of the Hubble Diagram at the model-independent level. We emphasise observational selection effects, which inevitably exist and distort the theoretically predicted relations. An example is the weak and strong gravitational lensing bias effect for high redshift LGRB in the presence of limited observational sensitivity (Malmquist bias). After bias correction, there is a tendency to vacuum dominated models with $\Omega_\Lambda \rightarrow 0.9$, $\Omega_\text{m} \rightarrow 0.1$. Forthcoming gamma-ray observations by the space THESEUS mission together with ground and space based multimessenger facilities will allow us to improve essentially the restrictions on alternative basic principles of cosmological models.

It is possible to reduce the discrepancy between the local measurement of the cosmological parameter $H_0$ and the value derived from the $Planck$ measurements of the Cosmic Microwave Background (CMB) by considering contamination of the CMB by emission from some medium around distant extragalactic sources, such as extremely cold coarse-grain dust. Though being distant, such a medium would still be in the foreground with respect to the CMB, and, as any other foreground, it would alter the CMB power spectrum. This could contribute to the dispersion of CMB temperature fluctuations. By generating a few random samples of CMB with different dispersions, we have checked that the increased dispersion leads to a smaller estimated value of $H_0$, the rest of the cosmological model parameters remaining fixed. This might explain the reduced value of the $Planck$-derived parameter $H_0$ with respect to the local measurements. The signature of the distant foreground in the CMB traced by SNe was previously reported by the authors of this paper -- we found a correlation between the SN redshifts, $z_{\rm SN}$, and CMB temperature fluctuations at the SNe locations, $T_{\rm SN}$. Here we have used the slopes of the regression lines $T_{\rm SN}\,/\,z_{\rm SN}$ corresponding to different {\it Planck} wave bands in order to estimate the possible temperature of the distant extragalactic medium, which turns out to be very low, about 5\,K. The most likely ingredient of this medium is coarse-grain ($grey$) dust, which is known to be almost undetectable, except for the effect of dimming remote extragalactic sources.

New arguments supporting the reality of large-scale fluctuations in the density of the visible matter in deep galaxy surveys are presented. A statistical analysis of the radial distributions of galaxies in the COSMOS and HDF-N deep fields is presented. Independent spectral and photometric surveys exist for each field, carried out in different wavelength ranges and using different observing methods. Catalogs of photometric redshifts in the optical (COSMOS-Zphot) and infrared (UltraVISTA) were used for the COSMOS field in the redshift interval $0.1 < z < 3.5$, as well as the zCOSMOS (10kZ) spectroscopic survey and the XMM-COSMOS and ALHAMBRA-F4 photometric redshift surveys. The HDFN-Zphot and ALHAMBRA-F5 catalogs of photometric redshifts were used for the HDF-N field. The Pearson correlation coefficient for the fluctuations in the numbers of galaxies obtained for independent surveys of the same deep field reaches $R = 0.70 \pm 0.16$. The presence of this positive correlation supports the reality of fluctuations in the density of visible matter with sizes of up to 1 000 Mpc and amplitudes of up to 20% at redshifts $z \sim 2$. The absence of correlations between the fluctuations in different fields (the correlation coefficient between COSMOS and HDF-N is $R = -0.20 \pm 0.31$) testifies to the independence of structures visible in different directions on the celestial sphere. This also indicates an absence of any influence from universal systematic errors (such as "spectral voids"), which could imitate the detection of correlated structures.

The methods of determining the fractal dimension and irregularity scale in simulated galaxy catalogs and the application of these methods to the data of the 2dF and 6dF catalogs are analyzed. Correlation methods are shown to be correctly applicable to fractal structures only at the scale lengths from several average distances between the galaxies, and up to (10-20)% of the radius of the largest sphere that fits completely inside the sample domain. Earlier the correlation methods were believed to be applicable up to the entire radius of the sphere and the researchers did not take the above restriction into account while finding the scale length corresponding to the transition to a uniform distribution. When an empirical formula is applied for approximating the radial distributions in the samples confined by the limiting apparent magnitude, the deviation of the true radial distribution from the approximating formula (but not the parameters of the best approximation) correlate with fractal dimension. An analysis of the 2dF catalog yields a fractal dimension of 2.20 +-m 0.25 on scale lengths from 2 to 20 Mpc, whereas no conclusive estimates can be derived by applying the conditional density method for larger scales due to the inherent biases of the method. An analysis of the radial distributions of galaxies in the 2dF and 6dF catalogs revealed significant irregularities on scale lengths of up to 70 Mpc. The magnitudes and sizes of these irregularities are consistent with the fractal dimension estimate of D = 2.1-2.4