Do Fractals Conﬁrm the General Theory of Relativity?

: The relatively high abundance of fractal properties of complex systems on Earth and in space is considered an argument in support of the general relativity of the geometric theory of gravity. The fractality may be called the fractal symmetry of physical interactions providing self-similarities of complex systems. Fractal symmetry is discrete. A class of geometric solutions of the general relativity equations for a complex scalar ﬁeld is o ﬀ ered. This class allows analogy to spatial fractals in large-scale structures of the universe due to its invariance with respect to the discrete scale transformation of the interval ds ↔ qd (cid:101) s . The method of constructing such solutions is described. As an application, the treatment of spatial variations of the Hubble constant H HST 0 (Riess et al., 2016) is considered. It is noted that the values H HST 0 form an almost fractal set. It has been shown that: a) the variation H HST 0 may be connected with the local gravitational perturbations of the space-time metrics in the vicinity of the galaxies containing Cepheids and supernovae selected for measurements; b) the value of the variation H HST 0 can be a consequence of variations in the space-time metric on the outskirts of the local supercluster, and their self-similarity indicates the fractal distribution of matter in this region. (or galaxies) that were used for the measurement. Variations of the expansion speed relative to the value H Pl 0 indicate the presence of local gravitational perturbations in the vicinity of each source (local changes of the space-time metric relative to the metric of the background reference system). The fractality indicates that the processes generating these gravitational perturbations have an identical nature. Such processes can be gravitational waves generated by the asymmetric collapse of a supernova core or cosmological gravitational perturbations. The fractality of the Hubble constant should also be observed in the fractal cosmological model.


Introduction
For more than a hundred years, on 29 May 1919, the Eddington expedition made observations of a solar eclipse. The purpose of the observation was to test one of the three effects proposed by Einstein to confirm his general relativity (GR)-the curvature of light lines from a distant star in the gravitational field of the Sun. The results of Eddington's observations confirmed the predictions of GR [1]. It was the first experiment to test the main postulate of this theory as a geometrical model of gravitational interaction. Over the past century, many direct experiments have been made. The absolute majority of these experiments confirmed Einstein's theory [2].
However, at present, criticism of GR has renewed, in particular, because of the appearance of, the so-called "dark sector" in cosmology. Astrophysicists and cosmologists introduced the dark matter, and then the dark energy, to explain the characteristics of the rotation curves of galaxies, the stability of galactic disks, the growth of large-scale structures in the early universe, the Hubble parameter (accelerated expansion of space), and the anisotropy of the CMB radiation. Here these problems will not be discussed; but some relevant publications are given in [3]. So far, the "dark sector" methodology could not provide a satisfactory explanation of these problems because of the large variety of dark energy and dark matter models.
Modifications of the GR appeared in response to the "paradoxes of the dark sector". The main disadvantage of these modifications is that they introduce additional fundamental constants of gravitational interaction, which do not allow the use of the effects of the "dark sector". However, these parameters are defined from the same astrophysical observations of galaxies which need to be explained. It does not cause satisfaction. In GR, the fundamental constants measured in terrestrial experiments are used.
Modern information technologies allow the numerical solution of the non-linear equations of GR or its modifications. They are also used in observations, in processing observational data, and in the analysis of these data within a chosen cosmological model. The authors of various cosmological simulations fit their results to observations with good accuracy, which becomes better with an increase in the number of free parameters for baryonic matter (for example, different types of viscosity), for dark matter (from cold to hot), and for dark energy (different variants of scalar field potentials).
In work [3], facts indicating that the dark matter hypothesis is groundless from the point of view of the interpretation of the observations are given. Galaxies and clusters of galaxies are located inside the cosmic web. They are open systems that interact with the circumgalactic medium and the intergalactic medium. The description of the dynamics of such open systems within the framework of GR has not yet been done. Therefore, there is no indisputable evidence that would justify the rejection of GR.
In this article, the relatively high abundance of fractal properties of complex systems on Earth and in space is considered as an argument in support of the GR as a geometric theory of gravity. Geometrically, fractals are self-similar structures. Not always this self-similarity can be observed directly, especially in astrophysical objects and cosmological structures. Meanwhile, a sign of self-similar discrete sets is power dependencies between the physical characteristics of complex systems. For example, let us consider the Tully-Fisher relation L ∼ v 3 in spirals for the luminosity L and the rotation rate of spiral galaxies v. Suppose that the values of the luminosities of galaxies and the rotation rates form geometric progressions (the most obvious example of self-similarity in the state space of the system): The Tully-Fisher relation follows from this formula if we assume ln a ln b = 3. Spiral galaxies are similar physical objects. The power of the Tully-Fisher relation indicates that these objects form a fractal discrete set in the state space of these objects; we observe the projection of this set onto the L − v plane as the Tully-Fisher relation.
The power of the relationship is often met in very different astrophysics objects. Let's give several examples:

•
Main sequence stars obey a mass-luminosity relation, L ∼ m α , α ≈ 3.52; • Period-luminosity relation for pulsating variable stars, L ∼ P α , for classical Cepheids α ≈ 1.15; • Luminosity functions for stars, galaxies, and quasars, Φ ∼ L α , −2.6 < α < −0.6; • Faber-Jackson relation for luminosity and their central stellar velocity dispersion of stars of elliptical galaxies, L ∼ σ 4 ; • Tully-Fisher relation for the luminosity and rotation rate of spiral galaxies, L ∼ v 3 ; • Brightness distribution in the images of bright elliptical galaxies, Power-law spectra of the radio emission of the jets from active galactic nuclei, F ∼ ν −α , 0.6 < α < 0.7 (jets are composed of plasma clouds with a magnetic field); • Baldwin's ratio in the active galactic nucleus for the equivalent width of the emission lines and the luminosity of its galaxy, EW(CaII) ∼ L −α , α > 1; • Red shift distribution of absorption lines in quasar spectra, dN dz ∼ (1 + z) α , 1.67 < α < 2.09; • Spatial correlation function of galaxies in the clusters of galaxies and clusters in superclusters, ξ(r) ∼ r −α , 1.6 < α < 2.2.
The given power relationships are of course correlations. However, new observations only specify them.
The relative high abundance of fractals indicates that perhaps fractality is a fundamental property of physical interactions. The fractality may be called the fractal symmetry of physical interactions providing self-similarity of complex structures. Fractal symmetry is discrete, since a complex system with fractal symmetry consists of self-similar subsystems and the physical parameters of these subsystems form geometric progressions. Fractality can be compared with quantization (discreteness) of the states of the systems in the microcosm.
Can fractal symmetry be a consequence of the geometric model of gravity in GR? This question is logical if we take into account the geometric properties of the phase spaces of the systems, which we call fractals. It is also logical from a historical point of view in the development of science. For example, in the classical work of Minkowski (1909, [4]) it was proposed to back out of the hypothesis of ether, but to combine space and time into a four-dimensional homogeneous and isotropic set of points (the Minkowski world). Minkowski proved his offer relying on postulates of the special theory of relativity and the Lorentz invariance of the Maxwell's equations. According to Minkowski, we discover a homogeneous and isotropic space-time through the invariance of Maxwell's equations with respect to the Lorentz transformations. From the modern point of view, Minkowski considered the Lorentz's symmetry for an electromagnetic field as a consequence of the properties of four-dimensional space-time: the vector potential of the electromagnetic field is also transformed like any vector in four-dimensional space-time. Theoretical cosmology as a science about the structure and evolution of space-time of the universe began with this research by Minkowski. It should be noted that the fundamental hypothesis by Weyl, where he proposed (1918, 1929 [5,6]) to use a scalar field Φ(x k ) (now it is called dilation) with transformation Φ ↔ σ 1/2 Φ and the continuous conformal symmetry of space-time ds 2 ↔ σ x k d s 2 for the electron model. This field should change the Newtonian gravitational interaction of all systems and its consequences would be self-similar properties of these systems. However, the existence of the dilaton field is not detected. Now symmetries play an important role in the development of our ideas about space-time and the physics of matter. The description of the observed properties of matter by means of dynamic equations and symmetries of the interaction fields allows us to understand the nature of the universe at the micro and macro levels. This paradigm inspires cosmologists.
In this paper, a class of geometric solutions of GR is proposed for a complex scalar field. The solutions admit the possibility of spatial fractals in the large-scale structure of the universe [7][8][9]. Quasar groups are an example of such fractals. It was found in [8], that the number of groups N i with size ϑ i depends on the size N i ∼ (ϑ i ) −2.02 .Let's consider the quasar group with number i as the so-called "island system", that has a space-time with an interval ds i . If the intervals of different island systems form a geometric progression ds i+1 = qds i (q is a number, discrete scale transformation), then a distant observer using the metric of the Minkowski's space will observe the geometric progression ϑ i+1 = qϑ i for the angular sizes of the quasar groups. The introduction of a geometric progression for interval means, that we consider the space-time, consists of an ensemble of geometrically self-similar regions.
Below we consider an example of the solution of the GR equations, which allows for the discrete scale transformation of the space-time interval. The complex scalar field is used as a source in the GR equations. It is assumed that the field has a unitary symmetry U (1), which is given in the algebraic form: The complex field is used in the scalar electrodynamics to introduce charges, that are conserved due to symmetry (1). In the Standard Model of the physics of elementary particles SU(n) symmetries are introduced as a generalization of the circle (1) for multidimensional fields.
In the problem considered here, the discrete transform of the field phase ϕ ↔ ϕ + α (α is a constant) without changing its amplitude Ψ (unitary symmetry) corresponds to a field transform ψ ↔ e −iα ψ and a discrete transform of metric tensor g mn (ϕ) ↔ U 0 U 0 g mn ( ϕ) , where U 0 and U 0 are the parameters of the field potentials U = U 0 ψψ * , U = U 0 ψ ψ * , an asterisk means complex conjugation. In this case, the two regions of space-time are geometrically similar, their intervals are connected by a scale transformation ds 2 ↔ U 0 U 0 d s 2 (an analogue of the conformal discrete transformation). An ensemble of such regions can be called a fractal, if ds 2 Note, that in the fractal cosmological model, the simulation of the large-scale structure of the galaxy distribution should require less computational time, because our need is to model the self-similar parts of this structure and not their ensembles as a whole.
The rest of the article is organized as follows. Section 2 presents the method of constructing a geometrical solution of the Lagrange and Einstein equations for a complex scalar field. Section 3 discusses the contradiction in the measurements of the Hubble constant and its interpretation within the framework of the fractal cosmological model.

Exact Geometrical Solution of the Lagrange and Einstein Equations for a Complex Scalar Field with Symmetry U (1)
In the preprint [10], a method is presented for constructing the class of exact solutions discussed here for the GR equations of a complex scalar field. The following describes the revised and expanded version of this method. Let's consider a system with the Hilbert-Einstein's action where the curvature scalar is R, g < 0 is the determinant of the metric tensor g mn , the space-time interval is ds 2 = g mn dx m dx n , the indices run through the values 0, 1, 2, 3, the metric signature is (+ − − −). Here, the physical system of units is maintained to simplify the comparison with the observations. The Lagrangian of a complex scalar field has the form: where h is the Planck constant, U(ψψ * ) is the field potential, which we choose in a simple form U(ψψ * ) = U 0 ψψ * = U 0 Ψ 2 for a field ψ = Ψe iϕ with symmetry (1). From the Lagrange equation we derive an equation for the field phase ϕ(x m ): The goal is to choose the dependence of the Christoffel's symbols Γ l mn and the metric tensor g mn on the derivatives ∂ϕ ∂x m in such away so that Equation (5) for the phase ϕ(x m ) is identically satisfied, and the Einstein equation for the metric tensor g mn can be reduced to an algebraic form for the vector e n = ∂ϕ ∂x n . In this case, the Einstein's equation is invariant relative to the discrete phase transformation ϕ ↔ ϕ + α . Achieving this goal is carried out in two stages.

Selection of the Christoffel's Symbols and Metric Tensor
By direct substitution, we can see that Equation (5) becomes an identity if the following formulas are used: e m e n g mn = U 0 , where b l is a vector. The left-hand sides of formulas (6) and (8) are the tensors, therefore the formulas are covariance invariant. For the covariant conservation of formula (7), it is necessary that the covariant derivatives be equal to zero, e n,m = 0, b l ,m = 0 and We define the metric tensor by the following formula: U 0 e m e n + a m e n + a n e m , where there a m is a vector. From the condition g mn g mn = 4, the equation follows, a m e m = 0. The equation for the covariant derivative g mn,k = 0 is identically satisfied, if e n,m = 0 and a m,n = 0, then Calculation of the Cristoffel's symbols Γ n ml = 1 2 g kn ∂g mk ∂x l + ∂g lk ∂x m − ∂g ml ∂x k and comparison with the definition (7) lead to the following equation: U 0 e n ∂e m ∂x l + a n ∂e m ∂x l + e n ∂a m ∂x l , if the calculations take into account the equalities (12). Equation (13) turns into an identity, if we take into account Equation (12) and b n = U 0 a n .
Formulas (6)- (14) determine the relationships between the vector e m and g mn , Γ l mn , a m , and b l .

Algebraic form of the Einstein's Equations
To determine the vector e n and metric tensor (10), one must use the Einstein's equation R km = 8πG c 4 T km − 1 2 g km T (here we use the notation of the classic textbook [11]). The energy-momentum tensor of the complex field is T km = 2Ψ 2 hc U 0 e k e m . The calculation of the Ricci tensor R mn through the symbols Γ l mn and Equation (12) allows us to obtain the Einstein's equation in the following form: ∂ 2 e m ∂x l ∂x n − ∂ 2 e m ∂x n ∂x l e l + U 0 a l = −ζU 0 (2e m e n + U 0 a m e n + U 0 a n e m ), where ζ = 8πG hc 5 Ψ 2 . The derivatives of the left side of Equation (15) are determined from Equation (9): In this case, Equation (15) is rewritten as follows: ∂e m ∂x k ∂e n ∂x l − ∂e m ∂x n ∂e k ∂x l e k + U 0 a k e l + U 0 a l = −ζ(U 0 ) 2 (2e m e n + U 0 a m e n + U 0 a n e m ).
The Equation (17) has an algebraic equation for a vector e m if we introduce another definition for the tensor ∂e m ∂x n : ∂e m ∂x n = e m γ n + e n γ m , where γ m is a vector. We also use the general form of the solution of Equation (11): where the vector ξ m does not depend on coordinates x k , i.e., ∂ξ m ∂x n and ξ m e m = 3U 0 . As a result, Equation (17) is written in algebraic form: where y = γ k ξ k and z = γ k e k . The vector e m can be found by convolving Equation (20) with the vector e n : Equations (18) and (21) define the class of geometrical solutions for problem (1)-(3) for the metric tensor g mn = 1 U 0 (−2e m e n + e n ξ m + e m ξ n ), and the Cristoffel's symbols if the vectors γ k and ξ k are chosen. Options, that lead to contradictions, should be excluded from this choice: y = 3z, y = 0, z = 0, y = 2z, and γ m γ m = 0. So, the Einstein's equation is reduced to the algebraic form (20). Using it and Equation (18), we find the function e m (x n ).
Convolution of Equation (20) with e m e n , γ m γ n , and ξ m ξ n gives three corresponding equations: 2(y − z) 2 = 3ζU 0 (6U 0 − ξ n ξ n ). is an imaginary number. Then the Equation (18) is reduced to the form (the condition ∂ξ m ∂x n = 0 is used): The Equation (27) includes three unknowns-two functions γ m (x k ), z(x k ), and a constant vector ξ m . Their choice determines the vector e m (x n ). As an example, we present the following choice of equations for functions γ m (x k ), z(x k ): They turn Equation (27) into an identity. The complex solution of Equation (29) has the form: The solution of Equation (28) is more bulky: where dη = ξ n dx n , the functions ξ n ξ n and y are determined from Equations (25)  Note that for the problem in question, Lagrangian (3) is zero, but the field energy density is not zero: The interval of the phase trajectory for the circle (1) is equal to Using formulas (22), one can transform formula (33) to the following form: Formula (34) shows how the local geometry of the phase space (1) is related to the local space-time geometry.
When the phase transformation is ϕ ↔ ϕ + α and the field amplitude is Ψ = Ψ, but the potential parameter is changed U 0 ↔ U 0 , then the Lagrange's equation, the Einstein's equations, covariant vectors e m , ξ m , γ m , the Christoffel's symbols do not change. However, the metric tensor is multiplied by a constant coefficient: As a result, the transformation of the phase of the complex field leads to the compression or stretching of the space-time interval ds 2 ↔ U 0 U 0 d s 2 . This means that the volumes of space-time with fields ψ and ψ are similar geometrically, therefore, the directly observed Newtonian law of law does not change. On the contrary, continuous conformal transformations distort the geometry and are accompanied by a change in the Newtonian law of gravitation.
Thus, if there are two regions of space-time in which complex scalar fields are connected by a phase transformation ϕ ↔ ϕ + α , then the spatial scales in these regions are connected by a transformation dl ↔ (dl) i , respectively. Then a sample of physically similar structures from these regions will form a fractal.

Spatial Variations of the Hubble Constant and Local Gravitational Perturbations
Hubble's article on "relation between distance and radial velocity among extra-galactic nebular" was published 90 years ago [12]. This discovery was the first law of the observational cosmology. The Hubble constant is the most important cosmological parameter. It is used to estimate cosmological distances to galaxies and their clusters. It is used in the theoretical cosmology through the Hubble parameter, H = . a a , where a is the scale factor of the cosmological model and a dot above the letter means the differentiation according to the cosmological time of the model t. The Hubble parameter characterizes the speed of the expansion of the universe space. In a homogeneous and isotropic cosmological model, the Hubble parameter depends only on cosmological time t (at each moment of time, the expansion rate at all points of the 3D-space is identical). Then in the modern era, the Hubble constant H 0 = H(t 0 ) has to be identical in all directions.
The results of the latest local measurements H 0 were very carefully performed in [13,14] using data on Cepheids in 23 galaxies, supernovae SN Ia were observed in 19 galaxies (redshifts 0.01 < z < 0.15, observations by the Hubble Space Telescope). The average value is equal to H HST 0 = (73.48 ± 1.66)km/s/Mpc for the standard cosmological model. Global estimates H 0 are obtained by fitting cosmological parameters of the model to observations of the CMB anisotropy and galaxy clustering. According to the latest data from the Planck collaboration for a spatially flat model, we have H Pl 0 = (67.4 ± 0.5)km/s/Mpc [15]. In recent publications, an unexpectedly large difference between H HST 0 and H Pl 0 is explained by the possibility of the presence of still unclear errors of measurements or hypotheses about the dark energy model (for example, [16]) or about the inhomogeneity of dark matter distribution (for example, [17,18]).
Here we propose a simpler hypothesis to explain the inequality H HST 0 H Pl 0 remaining within the GR framework and without the use of exotic forms of matter. The hypothesis itself appeared due to a careful look at the measurement data H HST 0 in various galaxies [13]. A sample of these data is shown in Table 1. It appears that the measured values H HST 0 form a manifold that is approximately described by the power of low characterizing fractals. It has the following form H HST 0 ≈ q i H * , and q = 1.023. Here H * = 65km/s/Mpc and the values of the exponent i are shown in Table 1. The proximity of value q = 1.023 to 1 may indicate a very small deviation from the isotropic expansion of space.
Let's suppose that this form reflects the real local properties of the space expansion in the vicinity of the supernovae (or galaxies) that were used for the measurement. Variations of the expansion speed relative to the value H Pl 0 indicate the presence of local gravitational perturbations in the vicinity of each source (local changes of the space-time metric relative to the metric of the background reference system). The fractality indicates that the processes generating these gravitational perturbations have an identical nature. Such processes can be gravitational waves generated by the asymmetric collapse of a supernova core or cosmological gravitational perturbations. The fractality of the Hubble constant should also be observed in the fractal cosmological model. Let's consider local gravitational perturbations in the vicinity of the galaxy which hosts Cepheids and supernovae. These sources are used to measure the Hubble constant in [13]. We use the synchronous space-time metric where Latin indices run through values 0, 1, 2, 3, and Greek 1, 2, 3. The coordinates of the background reference system are η and x α , the cosmological time t is related to the conformal time η by the Equation adη = cdt, and h αβ is a local gravitational perturbation. The commoving distance of galaxy r g is measured along the isotropic geodesic of the space-time with metric (1). For the isotropic geodesic ds = 0, then c 2 dt 2 = − g αβ + h αβ k α k β dr 2 . Here the isotropic wave vector is k i = dx i dr , k i k i = 0, and k α = q α a 2 , k α = q α , q α and q α don't depend on t, k α k α = q 2 a 2 . The commoving distance of galaxy r g is equal to  Table 1) was determined from the condition of equality of the distance modules for Cepheids and a supernova in the same galaxy.

Discussion
Above, it is proposed to consider the puzzle of fractals as a consequence of the properties of space-time and the symmetry of the interaction fields of matter. In this case, it is necessary to use the geometric theory of gravity, which is similar to GR. We have described the method of constructing the problems solution (1)-(3) and shown that it had simple and easily verified analytical solutions with geometrical properties. It is important that these solutions for complex scalar fields be generalized to multidimensional fields.
This class of geometric solutions of the general relativity equation allows analogy to spatial fractals in large-scale structures of the universe due to its invariance with respect to discrete scale transformation of interval, ds ↔ qd s . As an application, the treatment of spatial variations of the Hubble constant H HST