We perform an analysis of the Cosmic Web as a complex network, which is built on a $\Lambda$CDM cosmological simulation. For each of nodes, which are in this case dark matter halos formed in the simulation, we compute 10 network metrics, which characterize the role and position of a node in the network. The relation of these metrics to topological affiliation of the halo, i.e. to the type of large scale structure, which it belongs to, is then investigated. In particular, the correlation coefficients between network metrics and topology classes are computed. We have applied different machine learning methods to test the predictive power of obtained network metrics and to check if one could use network analysis as a tool for establishing topology of the large scale structure of the Universe. Results of such predictions, combined in the confusion matrix, show that it is not possible to give a good prediction of the topology of Cosmic Web (score is $\approx$ 70 $\%$ in average) based only on coordinates and velocities of nodes (halos), yet network metrics can give a hint about the topological landscape of matter distribution.

We investigate the radial density distribution of the dynamical dark energy inside the white dwarfs (WDs) and its possible impact on their intrinsic structure. The minimally-coupled dark energy with barotropic equation of state which has three free parameters (density, equation of state and effective sound speed) is used. We analyse how such dark energy affects the mass-radius relation for the WDs because of its contribution to the joint gravitational potential of the system. For this we use Chandrasekhar model of the WDs, where model parameters are the parameter of the chemical composition and the relativistic parameter. To evaluate the dark energy distribution inside a WD we solve the conservation equation in the spherical static metric. Obtained distribution is used to find the parameters of dark energy for which the deviation from the Chandrasekhar model mass-radius relation become non-negligible. We conclude also, that the absence of observational evidence for existence of WDs with untypical intrinsic structure (mass-radius relation) gives us lower limit for the value of effective sound speed of dark energy $c_s^2 \gtrsim 10^{-4}$ (in units of speed of light).

The recent detection of the 21-cm absorption signal by the EDGES collaboration has been widely used to constrain the basic properties of dark matter particles. However, extracting the parameters of the 21-cm absorption signal relies on a chosen parametrisation of the foreground radio emission. Recently, the new parametrisations of the foreground and systematics have been proposed, showing significant deviations of the 21-cm signal parameters from those assumed by the original EDGES paper. In this paper, we consider this new uncertainty, comparing the observed signal with the predictions of several dark matter models, including the widely-used cold dark matter (CDM) model, 1-3 keV warm dark matter models (WDM), and 7 keV sterile neutrino (SN7) model, capable of producing the reported 3.5 keV line. We show that all these dark matter models cannot be statistically distinguished using the available EDGES data.

Formation of halos in the Dark Ages from initial spherical perturbations is analyzed in a four component Universe (dark matter, dark energy, baryonic matter and radiation) in the approximation of relativistic hydrodynamics. Evolution of density and velocity perturbations of each component is obtained by integration of a system of nine differential equations from $z=10^8$ up to virialization, which is described phenomenologically. It is shown that the number density of dark matter halos with masses $M\sim10^8-10^9\,\mathrm{M_{\odot}}$ virialized at $z\sim10$ is close to the number density of galaxies in comoving coordinates. The dynamical dark energy of classical scalar field type does not significantly influence the evolution of the other components, but dark energy with a small value of effective sound speed can affect the final halo state. Simultaneously, the formation/dissociation of the first molecules have been analyzed in the halos which are forming. The results show that number densities of molecules $\rm{H_2}$ and $\rm{HD}$ at the moment of halo virialization are $\sim10^3$ and $\sim400$ times larger, respectively, than on a uniformly expanding background. It is caused by increased density and rates of reactions at quasilinear and nonlinear evolution stages of density and velocity of the baryonic component of halos. It is shown also that the temperature history of the halo is important for calculating the concentration of molecular ions with low binding energy. So, in a halo with virial temperature $\sim10^5$ K the number density of the molecular ion HeH$^+$ is approximately 100 times smaller than that on the cosmological background.

The formation of large voids in the Cosmic Web from the initial adiabatic cosmological perturbations of space-time metric, density and velocity of matter is investigated in cosmological model with the dynamical dark energy accelerating expansion of the Universe. It is shown that the negative density perturbations with the initial radius of about 50 Mpc in comoving to the cosmological background coordinates and the amplitude corresponding to the r.m.s. temperature fluctuations of the cosmic microwave background lead to the formation of voids with the density contrast up to $-$0.9, maximal peculiar velocity about 400 km/s and the radius close to the initial one. An important feature of voids formation from the analyzed initial amplitudes and profiles is establishing the surrounding overdensity shell. We have shown that the ratio of the peculiar velocity in units of the Hubble flow to the density contrast in the central part of a void does not depend or weakly depends on the distance from the center of the void. It is also shown that this ratio is sensitive to the values of dark energy parameters and can be used to find them based on the observational data on mass density and peculiar velocities of galaxies in the voids.

We analyse the evolution of cosmological perturbations which leads to the formation of large isolated voids in the Universe. We assume that initial perturbations are spherical and all components of the Universe (radiation, matter and dark energy) are continuous media with perfect fluid energy--momentum tensors, which interact only gravitationally. Equations of the evolution of perturbations for every component in the comoving to cosmological background reference frame are obtained from equations of energy and momentum conservation and Einstein's ones and are integrated numerically. Initial conditions are set at the early stage of evolution in the radiation-dominated epoch, when the scale of perturbation is much larger than the particle horizon. Results show how the profiles of density and velocity of matter and dark energy are formed and how they depend on parameters of dark energy and initial conditions. In particular, it is shown that final matter density and velocity amplitudes change within range $\sim$4-7\% when the value of equation-of-state parameter of dark energy $w$ vary in the range from --0.8 to --1.2, and change within $\sim$1\% only when the value of effective sound speed of dark energy vary over all allowable range of its values.

We analyse the evolution of cosmological perturbations which leads to the formation of large voids in the distribution of galaxies. We assume that perturbations are spherical and all components of the Universe - radiation, matter and dark energy - are continuous media with ideal fluid energy-momentum tensors, which interact only gravitationally. Equations of the evolution of perturbations in the comoving to cosmological background reference frame for every component are obtained from equations of conservation and Einstein's ones and are integrated by modified Euler method. Initial conditions are set at the early stage of evolution in the radiation-dominated epoch, when the scale of perturbation is mush larger than the particle horizon. Results show how the profiles of density and velocity of matter in spherical voids with different overdensity shells are formed.

We analyse the evolution of scalar field dark energy in the spherical halos of dark matter at the late stages of formation of gravitationally bound systems in the expanding Universe. The dynamics of quintessential dark energy at the center of dark matter halo strongly depends on the value of effective sound speed $c_s$ (in units of speed of light). If $c_s\sim1$ (classical scalar field) then the dark energy in the gravitationally bound systems is only slightly perturbed and its density is practically the same as in cosmological background. The dark energy with small value of sound speed ($c_s<0.1$), on the contrary, is important dynamical component of halo at all stages of their evolution: linear, non-linear, turnaround, collapse, virialization and later up to current epoch. These properties of dark energy can be used for constraining the value of effective sound speed $c_s$ by comparison the theoretical predictions with observational data related to the large scale gravitationally bound systems.

We investigate the non-linear evolution of spherical density and velocity perturbations of dark matter and dark energy in the expanding Universe. For that we have used the conservation and Einstein equations to describe the evolution of gravitationally coupled inhomogeneities of dark matter, dark energy and radiation from linear stage in the early Universe to the non-linear one at the current epoch. The simple method of the numerical integration of the system of non-linear differential equations for evolution of the central part of halo is proposed. The results are presented for halo of cluster ($k=2$ Mpc$^{-1}$) and supercluster scales ($k=0.2$ Mpc$^{-1}$) and show that quintessential scalar field dark energy with small value of effective speed of sound $c_s<0.1$ can give noticeable impact on the formation of large scale structures in the expanding Universe.

The distribution of dark energy density in the vicinity of compact static objects is analyzed. Dark energy is assumed to be in the form of a scalar field with three parameters: the background density, the equation of state parameter and the effective sound speed. Compact object is assumed to be a homogeneous spherical object of constant radius. We use the solutions of the hydrodynamical equations for dark energy in the gravitational fields of such objects for cases of static distribution of dark energy in the vicinity of star and stationary accretion onto black hole in order to analyze the possibility of constraining of the parameters of dark energy from astrophysical data. We show that dependence of density of dark energy in the vicinity of such object on the effective sound speed, background density and equation of state parameter of dark energy makes it possible to try such tests. Here we exploit the accuracy of determination of masses of Sun and black hole in the center of Milky Way to obtain the lower limit on the effective sound speed of dark energy.

We study the dynamical properties and space distribution of dark energy in the weak and strong gravitational fields caused by inhomogeneities of matter in the static world of galaxies and clusters. We show that the dark energy in the weak gravitational fields of matter density perturbations can condense or dilute, but amplitudes of its perturbations remain very small on all scales. We illustrate also how the "accretion" of the phantom dark energy onto the matter overdensity forms the dark energy underdensity. We analyze the behavior of dark energy in the gravitational fields of stars and black holes with the Schwarzschild metric. It is shown that, in the case of stars, the static solution of the differential equations for energy-momentum conservation exists and describes the distribution of density of dark energy inside and outside of a star. We have found that for stars and galaxies its value differs slightly from the average and is a bit higher for the quintessential scalar field as dark energy and a bit lower for the phantom one. The difference grows with the decrease of the effective sound speed of dark energy and is large in the neighborhood of neutron stars. We obtain and analyze also the solutions of equations that describe the stationary accretion of the dark energy as a test component onto the Schwarzschild black hole. It is shown that the rate of change of mass of the dark energy is positive in the case of quintessential dark energy and is negative in the case of the phantom one.

We analyze the behavior of the scalar field as dark energy of the Universe in a static world of galaxies and clusters of galaxies. We find the analytical solutions of evolution equations of the density and velocity perturbations of dark matter and dark energy, which interact only gravitationally, along with the perturbations of metric in a static world with background Minkowski metric. Using them it was shown that quintessential and phantom dark energy in the static world of galaxies and clusters of galaxies is gravitationally stable and can only oscillate by the influence of self-gravity. In the gravitational field of dark matter perturbations it is able to condense monotonically, but the amplitude of density and velocity perturbations on all scales remains small. It was illustrated also, that the "accretion" of phantom dark energy in the region of dark matter overdensities causes formation of dark energy underdensities - the regions with negative amplitude of density perturbations of dark energy.