The emergence of large-scale low-dimensional structures in modern scientific literature is associated with fluctuation inhomogeneities in matter, which have passed from the inflationary era to the new cosmological era of flat space. At the same time, the dominant paradigm for more than 100 years was the Friedmann model of homogeneous space (in the early stages, it was transformed into the De Sitter model). However, by now an extremely large amount of observational material has been obtained, indicating that the concept of “homogeneity” of the Universe requires a special interpretation. The presence of (quasi)order in the form of ubiquitous large low-dimensional (coherent) structures instead of the expected homogeneous distribution of matter in the region of the Universe accessible to observation suggests the existence of at least background physical mechanisms that effectively compete over long time intervals with the Friedmann expansion (in terms of the spatial distribution of gravitating matter).

The consideration of cosmological models has historically been associated with the analysis of the properties of metrics, which are solutions of Einstein’s field equations. The distribution of many masses, which creates a fairly arbitrary, and time-varying gravitational field, together with the field equations can be solved analytically only in exceptional cases, and their numerical solution requires colossal labor costs even with the use of the most modern computing technology. However, there are powerful analytical methods for the qualitative analysis of the properties of solutions to systems of such equations, which can be used to study the properties of astrophysical processes, leaving regions with complex geometry or with fast processes for machine modeling. However, the primary issue for the correct writing of the Vlasov-Einstein system of matter-field equations with a self-consistent field at the stage of cosmological structure formation is the identification of the dependencies of the metric coefficients at the current time stage. At present, the prevailing hypothesis about the initial epoch of the expansion of the Universe is the inflationary one, which uses for its explanation a certain formal “inflaton field”, which goes beyond the framework of the classical model ΛCDM. However, based on the concept of a self-consistent field, one can, relying on the ideas of A. A. Vlasov, consider the possibility of realizing the de Sitter stage of expansion as a special form of action of pseudo-forces of metric origin, caused by a special initial distribution of particles on the manifold ${ℳ}_8$ of Laptev supporting elements [1] [2] (a generalization of the E. Cartan space), which is a submanifold of the tangent bundle over the configuration set of the system of particles ${ℳ}_x$; In this article we will limit ourselves to introducing on it the usual Riemannian metric [3] , the operations with which do not require the definition of new concepts.

It can be shown that for long time intervals (including the modern era) the equations with a self-consistent field smoothly transform into the Vlasov-Poisson equations. In this case, their integral Volterra form has, under certain, quite physically obvious and reasonable assumptions, periodic in time solutions [4] . In the region of structure growth, a local change in the metric occurs, caused by this process (as a phase transition of matter in hydrodynamic macroparameters), and it is possible to trace the dynamics of the development of structures from their birth to decay (caused both by the expansion of space and by the increase in internal desynchronization of the constituent parts of macrostructures—which is an essential important part of their genesis: statistical decoherence of multi-time ensembles of particles at large distances).

In this paper, we consider a method for analyzing a system of kinetic Vlasov equations field and equations for the gravitational field, which consistently presents the cosmological aspects of the approaches proposed in the works [5] [6] , and allows us to obtain in a very simple way the general evolution of the observed macrostructures in the Universe, as well as to establish general criteria for their development, using only completely transparent analytical approaches.

We consider the self-consistent system of Vlasov-Einstein equations for a set of massive particles described by the distribution function $F\left(x,v\right)$ ( $x\in X\subseteq {ℝ}_x^4$, $V\subseteq {ℝ}_v^4$):

${\text{div}}_4\left(vF\left(x,v\right)\right)=0,R^{\mu \nu }-\frac{1}{2}{R}^{†}{g}^{\mu \nu }=\kappa T^{\mu \nu },R^{†}\equiv R-2\text{Λ},\kappa =\frac{8\pi G}{c^4},$(1)

where: ${\text{div}}_4\left(\mathrm{...}\right)$—covariant divergence operator:

${\text{div}}_4\left(vF\right)=v^{\mu }{\nabla }_{\mu }F+F{\nabla }_{\mu }{v}^{\mu },{\nabla }_{\mu }F=\frac{\partial F}{\partial x^{\mu }}+\frac{\partial F}{\partial v^{\nu }}\frac{\partial v^{\nu }}{\partial x^{\mu }},\mu ,\nu ,\cdots =\overline{0,3}$(2)

${\nabla }_{\eta }{v}^{\nu }=\frac{\partial v^{\nu }}{\partial x^{\eta }}+{\text{Γ}}_{\eta \zeta }^{\nu }{v}^{\zeta }\equiv w_{\eta }^{\nu }\left[x^{\eta },v^{\nu },f\right],{\nabla }_{\eta }\left(v^{\nu }F\right)=v^{\nu }{\nabla }_{\eta }F+F\left(\frac{\partial v^{\nu }}{\partial x^{\eta }}+{\text{Γ}}_{\eta \zeta }^{\nu }{v}^{\zeta }\right).$(3)

Multiplication by $v^{\eta }\equiv \text{d}{x}^{\eta }/\text{d}\tau$ and convolution of the first Equation (3) by index $\eta$ leads to the classical form of the Euler-Lagrange equations given in manuals on general relativity, however, on the right-hand side there will be a quantity with the dimension of acceleration: $v^{\eta }{w}_{\eta }^{\nu }\equiv {\mathfrak{a}}^{\nu }$. In this case, the system (3) (with respect to 4-velocities) is overdetermined ( $\max \mu \cdot \max \nu >\dim v$) and requires the introduction of additional integrability conditions. As such, we can take the conditions $\hat{\Omega }\left[w\right]\equiv {\nabla }_{\nu }{w}_{\mu }^{\eta }-{\nabla }_{\mu }{w}_{\nu }^{\eta }=0$. Since for (3) we have

${\nabla }_{\nu }{\nabla }_{\mu }{v}^{\eta }-{\nabla }_{\mu }{\nabla }_{\nu }{v}^{\eta }={\nabla }_{\nu }{w}_{\mu }^{\eta }-{\nabla }_{\mu }{w}_{\nu }^{\eta }$, ${\left(v^{\eta }\right)}^{\prime \text{}\prime }{}_{\nu \mu }={\left(v^{\eta }\right)}^{\prime \text{}\prime }{}_{\mu \nu }$, then $\hat{\Omega }\left[w\right]=R_{\alpha \mu \nu }^{\eta }{v}^{\alpha }=0$

(for a flat space the integrability condition is satisfied automatically, because in this case $R_{\alpha \mu \nu }^{\eta }=0$).

In the case of a space with a non-vanishing Riemann tensor, there is a special form of the tensor $w_{\nu }^{\eta }={}_{\left(0\right)}{w}_{\nu }^{\eta }$, which corresponds to the velocity field emanating from one point in the absence of dynamic forces (accelerations) in the system:

${}_{\left(0\right)}{w}_{\nu }^{\eta }=\mathscr{H}\left({\delta }_{\nu }^{\eta }-{v^{\eta }{v}_{\nu }|}_{v=v\left[{}_{\left(0\right)}{w}_{\nu }^{\eta }\right]}/c^2\right)$, $\mathscr{H}=\text{const}$. The latter, according to is a consequence of the energy conservation condition (in terms of momenta $p=p\left(v\right)$) on “the mass hyperboloid ${ℳ}_v$” [7] , and also, according to [5] , the fact that the expression for the acceleration there automatically cancels out:

$\begin{array}{l}{v^{\eta }{\nabla }_{\eta }{v}^{\nu }|}_{v=v\left[{}_{\left(0\right)}{w}_{\nu }^{\eta }\right]}={\left(\frac{\text{d}{v}^{\nu }}{\text{d}\tau }+{\text{Γ}}_{\eta \mu }^{\nu }{v}^{\eta }{v}^{\mu }\right)|}_{v=v\left[{}_{\left(0\right)}{w}_{\nu }^{\eta }\right]}\\ ={v^{\eta }\mathscr{H}\left({\delta }_{\eta }^{\nu }-v_{\eta }{v}^{\nu }/c^2\right)|}_{v=v\left[{}_{\left(0\right)}{w}_{\nu }^{\eta }\right]}=\mathscr{H}\left(v^{\nu }-v^{\nu }\right)=0.\end{array}$

On the other hand, convolution of the fundamental tensor $g_{\nu \mu }$ with the Equation (3) for ${}_{\left(0\right)}{w}_{\nu }^{\eta }$ gives:

${\nabla }_{\eta }c\equiv 0,{\nabla }_{\eta }{e}_{\mu }=\mathscr{H}{c}^{-1}\left(g_{\mu \eta }-e_{\mu }{e}_{\eta }\right),e_{\eta }{e}^{\eta }=1.$(4)

In the geodetic coordinate system ${\mathfrak{K}}_4$, defined by a congruence of geodesics, one of which (parameterized by the natural real parameter $s$) is perpendicular to

the 3-dimensional hypersurface on which the coordinates ${{}_{\left(g\right)}{x}^k|}_{k=1,2,3}$ are given

(the ${\mathfrak{K}}_4$ system describes the entire 4-space), we have a general expression for

the coordinates $s={\mathfrak{X}}^0\equiv x^0$, ${}_{\left(g\right)}{x}^k={\mathfrak{X}}^k\left(x^0,\cdots ,x^3\right)={x^k|}_{k=1,2,3}$ (redesignation of coordinates). The interval in it can be expressed in the form: $\text{d}{s}^2=g_{00}{\left(\text{d}{x}^0\right)}^2-{\mathfrak{G}}_{ik}\text{d}{\mathfrak{X}}^i\text{d}{\mathfrak{X}}^k=g_{\mu \nu }\text{d}{x}^{\mu }\text{d}{x}^{\nu }$, if $g_{\mu \nu }={\mathfrak{G}}_{\mu \nu }+\underset{\text{const}}{\underset{︸}{e_{\mu }{e}_{\nu }}}$, since $e_{\nu }=e^{\nu }=\left(1,0,0,0\right)$ in the geodetic coordinate system, ${g_{\mu 0}|}_{\mu >0}\equiv 0$ ( $g_{00}\equiv 1$).

Thus, from the second equation of the system (4), and the last obtained relation

it follows that ${\nabla }_{\eta }{e}_{\nu }=\mathscr{H}{c}^{-1}{\mathfrak{G}}_{\mu \nu }$. Since in geodetic coordinates $g_{00}=1$, ${g_{0k}|}_{k>0}\equiv 0$, then ${\nabla }_{\eta }{e}_{\nu }=-{\text{Γ}}_{\eta \nu }^{\mu }{e}_{\mu }=\frac{1}{2}\partial g_{\eta \nu }/\partial x^0$.

Consequently, $\partial g_{\eta \nu }/\partial x^0=\partial \mathfrak{G}/\partial x^0$; finally we obtain ${\nabla }_{\eta }{e}_{\nu }=\frac{1}{2}\partial {\mathfrak{G}}_{\eta }{e}_{\nu }/\partial x^0=\left(\mathscr{H}/c\right){\mathfrak{G}}_{\eta }{e}_{\nu }$. Assuming $c,\mathscr{H}=\text{const}$, the last equation can

be considered as an ordinary differential equation, the solution of which is expressed through an arbitrary tensor function of three arguments ${\text{Ξ}}_{\eta \nu }$: ${\mathfrak{G}}_{\eta }{e}_{\nu }=\text{exp}\left(2\mathscr{H}{x}^0/c\right){\text{Ξ}}_{\eta \nu }\left({\mathfrak{X}}^1,{\mathfrak{X}}^2,{\mathfrak{X}}^3\right)$. Interval for the divergent velocity field of the medium

$\text{d}{s}^2={\left(\text{d}{x}^0\right)}^2-\text{exp}\left(2\mathscr{H}{x}^0/c\right)\cdot {\text{Ξ}}_{\eta \nu }\left({\mathfrak{X}}^1,{\mathfrak{X}}^2,{\mathfrak{X}}^3\right)\text{d}{x}^{\eta }\text{d}{x}^{\nu }.$(5)

For Riemannian geometry, the statement about local equivalence in the neighborhood of a chosen point $\mathfrak{O}\in {\mathcal{X}}_2$ is valid on the hypersurface ${\mathcal{X}}_2=\left\{x^0=\text{const}\right\}$ the above metric and the quasi-Euclidean metric of the local neighborhood $O\left(\tilde{\mathfrak{O}}\right)\subset {\mathcal{X}}_2$

$\text{d}{s}^2={\left(\text{d}{x}^0\right)}^2-\underset{k=1,2,3}{{\displaystyle \sum }}\text{d}{\left({\tilde{\mathfrak{X}}}^k\right)}^2,{\tilde{\mathfrak{X}}}^k\equiv \text{exp}\left(\mathscr{H}{x}^0/c\right){\mathfrak{X}}^k$(6)

(the “tilde” symbol will denote coordinate scaling). Spatial coordinates change on manifolds to which the “time/evolutionary” axis is orthogonal in accordance with the law ${\mathfrak{X}}^k\left({\mathfrak{X}}^0\right)=\text{exp}\left(\mathscr{H}{\mathfrak{X}}^0/c\right)\cdot {\mathfrak{X}}^k\left({\mathfrak{X}}^0=0\right)$ with the expansion rate $\text{d}{\mathfrak{X}}^k/\text{d}{\mathfrak{X}}^0=\left(\mathscr{H}{\mathfrak{X}}^0/c\right){\mathfrak{X}}^k$. Obviously, the similarity to the De Sitter metric cannot be just a coincidence.

The condition $\hat{\Omega }\left[w\right]=R_{\alpha \zeta ,\mu \nu }{v}^{\alpha }=0$ leads, after applying the covariant derivative to it and using the Bianchi identity, to a covariant differential equation for the Riemann tensor $v^{\eta }{\nabla }_{\eta }{R}_{\alpha \zeta ,\mu \nu }+2\mathscr{H}{R}_{\alpha \zeta ,\mu \nu }=0$; for the scalar curvature $R$, from here we have an ordinary differential equation (along the evolutionary axis of the geodetic coordinate system, which leads to a formal replacement of $\nabla \to \text{d}/\text{d}{\mathfrak{X}}^0$) for the independent variable ${\mathfrak{X}}^0$), the solution of which is: $R\left({\mathfrak{X}}^0\right)=R\left(0\right)\text{exp}\left(-\mathscr{H}{x}^0/c\right)$. Thus, the primary curvature of the 4–space should decrease with time (if $\mathscr{H}=\text{const}>0$).

It should be noted that the assumptions made above do not include the assumption of isotropy and homogeneity of space. The only important factor that complements the classical consideration of the metric of expanding space-time in general relativity is the assumption of an initially present diverging field of matter velocities at a certain initial moment of consideration (that is, in fact, a specialized Cauchy condition). In this case, matter moves “in a non-force regime” (by inertia), without taking into account the influence of the self-consistent field. The duration of the expansion phase of the need for significant consideration of the energy-momentum tensor can be estimated by the value of the parameter $\mathscr{H}$.

Recent observations of the inhomogeneous structure of the Universe [8] can lead us to a model in which at some point there are many regions with matter in a state with the distribution of particle velocities that were described above. It is quite possible that the synchronization of the mentioned regions is absent, and the competition of the inertial motion of matter particles with their motion in a self-consistent field leads (for sufficiently long times) to the establishment of a regime of effective negative temperature of the system, which leads to the appearance of megastructures such as walls and voids. It should be noted that in this case the function $\mathscr{H}\to \left\{\mathfrak{h}\left(z\right)\right\}$, where $\left\{\mathfrak{h}\left(z\right)\right\}$ is a some set of local indicators of inertial expansion depending on internal parameters (density, temperature, etc.); its may lead to different growth rates of cosmological systems and demonstrate the so-called “Hubble tension” effect (the reasons for this conclusion are given in the next paragraph).

The classical diffusion equation for processes (with characteristic time ${\mathfrak{t}}_0$) in 3-dimensional space can be obtained from the continuity equation

$\partial \rho \left({}_{3}x,t\right)/\partial t+\hat{N}\left[{}_{3}v;f\right]=0,\hat{N}\left[{}_{3}v;f\right]\equiv {\displaystyle \int \text{div}}\left({}_{3}vf\left({}_{3}x,{}_{3}v,t\right)\right){\text{d}}_{3}v,$

when passing to finite differences: $\hat{N}\left[{}_{3}v;f\right]\to {\displaystyle {\sum }_{j=1}^3{\displaystyle \int {}_3{v}_j\left(f\left({}_{3}x+{}_{3}v{\mathfrak{t}}_0+\cdots \right)-f\left({}_{3}x\right)\right)/\left({}_{3}v{\mathfrak{t}}_0+\cdots \right){\text{d}}_{3}v}}$. After the transformations we obtain $\partial \rho /\partial t+\text{div}\left({}_3{v}^{†}+{}_3{V}_D^{†}\right)\rho =\text{Ψ}\left(\rho \right)$, where ${}_3{v}^{†}={\displaystyle \int {}_{3}vf{\text{d}}_{3}v/\rho }$, ${}_3{V}_D^{†}\rho =-D{\nabla }_3\rho$, $D=\left({\mathfrak{t}}_0/2\right){\displaystyle \int {}_3{v}_j^{2}f{\text{d}}_{3}v/\rho }$. For the region of

cosmological macrostructure growth against the background of a homogeneous medium in the inertial regime (if we neglect the presence of force interaction from sources of an external or self-consistent gravitational field, for example, in a channel of matter motion with an ordered average velocity), the covariant diffusion

equation has the form ${\nabla }_4\left(v+V_D\right)\varrho =\text{Ψ}\left(\rho \right)$, where: $\varrho \left(x\right)=\text{d}N/\left({\left|g\right|}^{1/2}\text{d}x\right)$—the particle density in the 4-volume, ${\nabla }_{\mu }{v}^{\nu }=w_{\mu }^{\nu }$, $V_D^{\nu }\varrho =-{\mathcal{D}}^{\mu \nu }{\nabla }_{\nu }\varrho$, ${\mathcal{D}}^{\mu \nu }$—diffusion tensor, $V_D^{\nu }$—diffusion flow velocity, function $\text{Ψ}\left(\varrho \right)={\nu }_{ph}\left(\varrho -{\varrho }_0\right){{\left(1-\varrho /{\varrho }_{\max }\right)}^n|}_{n=1,2}$ displays the phase transition process

“isotropic medium—ordered structure” ( ${\nu }_{ph}$ is phenomenologic coefficient, ${\varrho }_0,{\varrho }_{\max }$ are concentrations of particles in both phases). For us, however, it is essential that the equation for the velocity vector of the medium by multiplying both of its parts by $\text{d}{x}^{\mu }/\text{d}\tau$ and folding by index $\mu$, leads, as previously noted, to the classical form of the Euler-Lagrange equation of geodesics, and on the right-hand side there is a quantity with the dimension of acceleration $v^{\beta }{w}_{\beta }^{\nu }$; this system of 16 equations requires additional restrictions on the Christoffel coefficients, since it is significantly overdetermined ( $\dim v=4<16$). We will restrict ourselves to the form of the tensor already considered earlier when imposing the

necessary additional conditions: $w_{\mu }^{\nu }=\mathfrak{h}{\left({\delta }_{\mu }^{\nu }-v^{\nu }{v}_{\mu }/v^2\right)|}_{v^2=v^{\nu }{v}_{\mu }}$ ( $\mathfrak{h}=\text{const}$, the

diagonal structure of the tensor $w_{\mu }^{\nu }=\text{diag}\left(\mathfrak{h},\cdots ,\mathfrak{h}\right)$). Obviously, in this case for $\mathfrak{h}\ne 0$ we postulate a situation in which there is no rest state (the growth of the structure, ensured by the condition $\nabla v=w_{\nu }^{\nu }>0$, continues indefinitely).

Let us introduce a geodetic coordinate system ${\mathfrak{K}}_4$ (similar to that in item 2). In the congruence of geodesics, we choose a direction with a unit vector $e^{\mu }\left(=\left(1,0^3\right)\right)$ and perform a parametrization on the chosen curve by a real parameter $s$, on the 3-hypersurface drawn normally to the direction $e^{\mu }$ we

introduce the coordinates ${{}_{\left(g\right)}{x}^k|}_{k=1,2,3}$. An arbitrary point in 4-space can now be described in new coordinates $s={\mathfrak{X}}^0$, ${}_{\left(g\right)}{x}^k={\mathfrak{X}}^k\left(x^0,\cdots ,x^3\right)$. The interval element in this case has the form $\text{d}{s}^2=\text{d}{\left({\mathfrak{X}}^0\right)}^2-{\mathfrak{G}}_{ij}\left({\mathfrak{X}}^{\mu }\right)\text{d}{\mathfrak{X}}^i\text{d}{\mathfrak{X}}^j$ ( $\mu =\overline{0,3}$; $i,j=\overline{1,3}$). The metric coefficients can, as before, be written in the form $g_{\mu \nu }={\mathfrak{G}}_{\mu \nu }+e^{\mu }{e}^{\nu }$.

Thus, the covariant derivative

${\nabla }_{\mu }\left({}_{\left(e\right)}{v}_{\nu }\right)=\mathfrak{h}{\left|v\right|}^{-1}\left({\mathfrak{G}}_{\mu \nu }+e_{\mu }{e}_{\nu }-\left({}_{\left(e\right)}{v}_{\mu }\right)\left({}_{\left(e\right)}{v}_{\nu }\right)\right),v_{\left(e\right)}{v}^{\nu }=v^{\nu }.$

Using the coordinate system introduced earlier, solutions for the velocities can

be sought in the form ${}_{\left(e\right)}{v}_{\nu }={}_{\left(e\right)}{v}^{\nu }=e^{\nu }=\left(1,0^3\right)$, so that the equation with a covariant derivative is reduced to the form $\partial {\mathfrak{G}}_{\mu \nu }/\partial {\mathfrak{X}}^0=\left(2\mathfrak{h}/\left|v\right|\right){\mathfrak{G}}_{\mu \nu }$. From here we obtain the general form of the metric tensor $\partial {\mathfrak{G}}_{\mu \nu }={}_{\left(1\right)}\mathfrak{Q}\left({\mathfrak{X}}^0\right)\cdot {}_{\left(2\right)}{\mathfrak{Q}}_{\mu \nu }\left({\mathfrak{X}}^1,{\mathfrak{X}}^2,{\mathfrak{X}}^3\right)$, where ${}_{\left(1\right)}\mathfrak{Q}\left({\mathfrak{X}}^0\right)=\text{exp}\left(\mathfrak{h}{\mathfrak{X}}^0/\left|v\right|\right)$, ${}_{\left(2\right)}{\mathfrak{Q}}_{\mu \nu }\left(\mathrm{...}\right)$ are some admissible functions of their arguments).

Non-zero Christoffel symbols of the second kind for this metric:

${\text{Γ}}_{0j}^i={\delta }_j^i{\left({}_{\left(1\right)}\mathfrak{Q}\right)}^{\prime }{}_{{\mathfrak{X}}^0}/{}_{\left(1\right)}\mathfrak{Q},{\text{Γ}}_{ij}^0={}_{\left(2\right)}{\mathfrak{Q}}_{ij}{}_{\left(1\right)}\mathfrak{Q}{\left({}_{\left(1\right)}\mathfrak{Q}\right)}^{\prime }{}_{{\mathfrak{X}}^0},{\text{Γ}}_{ij}^k={\text{Γ}}_{ij}^k\left({}_{\left(2\right)}{\mathfrak{Q}}_{ij}\right).$

The features of the motion of particles located in the region of a growing large structure can be determined from the Euler-Lagrange equations:

$\frac{\text{d}{v}^i}{\text{d}{\mathfrak{X}}^0}+{\text{Γ}}_{kj}^i{v}^k{v}^j+2\left({\left({}_{\left(1\right)}\mathfrak{Q}\right)}^{\prime }{}_{{\mathfrak{X}}^0}/{}_{\left(1\right)}\mathfrak{Q}\right)v^i{v}^0={\mathfrak{a}}^i.$

For small velocities, the 2nd term on the left-hand side can be neglected (compared to the 3rd), and we have an ordinary inhomogeneous differential equation, the solution of the evolution problem (analogically to crystals or plasmoids as in [5] ) can be written in quadratures:

$v^i\left({\mathfrak{X}}^0\right)~v^i\left({\mathfrak{X}}^0=0\right)\text{exp}\left(-{\displaystyle {\int }_0^{{\mathfrak{X}}^0}{\hat{\omega }}_1}\left({}_{\left(1\right)}\mathfrak{Q}\right)\text{d}{\mathfrak{X}}^0\right)+{\displaystyle {\int }_0^{{\mathfrak{X}}^0}{\hat{\omega }}_2}\left({}_{\left(1\right)}\mathfrak{Q}\left({\mathfrak{X}}_{\text{*}}^0\right),{\mathfrak{a}}^i\left({\mathfrak{X}}_{\text{*}}^0\right)\right)\text{d}{\mathfrak{X}}_{\text{*}}^0,$

${\hat{\omega }}_1\left({}_{\left(1\right)}\mathfrak{Q}\right)=2\left({\left({}_{\left(1\right)}\mathfrak{Q}\right)}^{\prime }{}_{{\mathfrak{X}}^0}/{}_{\left(1\right)}\mathfrak{Q}\right)v^0,$

${\hat{\omega }}_2\left({}_{\left(1\right)}\mathfrak{Q}\left({\mathfrak{X}}_{\text{*}}^0\right)\right)=\text{exp}\left(-{\displaystyle {\int }_0^{{\mathfrak{X}}^0-{\mathfrak{X}}_{\text{*}}^0}{\hat{\omega }}_1\left({}_{\left(1\right)}\mathfrak{Q}\right)}\right){\mathfrak{a}}^i\left({\mathfrak{X}}_{\text{*}}^0\right)\text{d}{\mathfrak{X}}_{\text{*}}^0.$

This gives rise to the effect of a non-force metric impact on particles, leading to their slowing down in the region of growth of the cosmological structure (at

${\left({}_{\left(1\right)}\mathfrak{Q}\right)}^{\prime }{}_{{\mathfrak{X}}^0}/{}_{\left(1\right)}\mathfrak{Q}>0$).

In addition to the simplest spherically symmetric conditions on the velocity, one should also consider filamentary structures, to which a significant number of large-scale cosmological objects observed by astronomers belong. In this case, one should consider the interval element

$\text{d}{s}^2={\left(\text{d}{\mathfrak{X}}^0\right)}^2-{\left(\text{d}{\mathfrak{X}}^3\right)}^2-{\left({}_{\left(I\right)}\mathfrak{Q}\left({\mathfrak{X}}^0\right)\right)}^2\underset{i,j=1}{\overset{2}{{\displaystyle \sum }}}{}_{\left(II\right)}{\mathfrak{Q}}_{ij}\left({\mathfrak{X}}^1,{\mathfrak{X}}^2\right)\text{d}{\mathfrak{X}}^i\text{d}{\mathfrak{X}}^j.$

In this case, a significant number of Christoffel coefficients are non-zero:

${{\text{Γ}}_{j,k}^i|}_{i,j,k=1,2},{\text{Γ}}_{2,4}^2,{{\text{Γ}}_{j,k}^4|}_{j,k=1,2}.$

Their singular points are determined by the conditions: 1) ${}_{\left(I\right)}\mathfrak{Q}=0$, 2) ${}_{\left(II\right)}{\mathfrak{Q}}_{11}\cdot {}_{\left(II\right)}{\mathfrak{Q}}_{22}-{}_{\left(II\right)}{\mathfrak{Q}}_{12}=0$. These same conditions determine the singularities for the elements of the Riemann curvature tensor $R_{12,12}$, $R_{14,14}$, $R_{14,24}$, $R_{24,24}$. Thus, at these points, conditions arise for the growth of the cosmological structure—in particular, the first condition describes the axis of symmetry of the mega-object, the second—the degeneration of the first quadratic form ( $EG-F^2=0$—the simplest case is plane or “Zel’dovich pancake”) of the two-dimensional surface (the stratification of surfaces). Moreover, at a singular point of the 2nd type, the scalar curvature locally changes sign, which indicates the reverse influence of the process of structure emergence on the conditions created there (the transformation of the singularity into the starting point of the process).

We can write out the Euler-Lagrange equations for the new metric to determine the motion of particles in the vicinity of the cosmological structure:

$\frac{\text{d}{v}^i}{\text{d}{\mathfrak{X}}^0}+{\text{Γ}}_{20}^2{v}^2{v}^0+{\text{Γ}}_{10}^1{v}^1{v}^0={\mathfrak{a}}^i,{{\text{Γ}}_{j0}^j|}_{j=1,2}=\frac{1}{2}{\left(\text{ln}\left({}_{\left(I\right)}\mathfrak{Q}\left({\mathfrak{X}}^0\right)\right)\right)}_{{\mathfrak{X}}^0}.$

Its solution will completely coincide with the previously obtained one in terms of structure.

At present, the prevailing concept of describing the early stage of the Universe’s evolution is the inflationary theory. But its implementation requires the fulfillment of very unusual conditions. At the same time, A. A. Vlasov, within the framework of his theory, proposed a simple model of biological and crystalline structures, which can be formulated for cosmology and astrophysical objects, such as topological associations of voids and cosmic walls and the cosmic web. In this paper, we consider the mechanism of inertial motion in a system of many gravitationally bound particles under the condition that special requirements are met for the initial distribution of matter. It is shown that in this case the emerging effect of forceless motion can be used to explain the structure of cosmological objects.