Perturbed FLRW Metric Explains the Difference in Measurements of the Hubble Constant

1. Introduction

In the $\mathsf{\Lambda }$CDM (Lambda Cold Dark Matter) cosmological model, based on the Friedmann-Lemaître-Robertson-Walker metric, matter moves synchronously with the expanding space. Perturbed metrics, [1,4] are a generalization of the FLRW metric. Deviations from standard cosmology are investigated by considering flat FLRW space-time in Painlev$\acute{e}$–Gullstrand coordinates [5]. With this coordinate choice the space is not expanding, although the distance between galaxies is certainly increasing. In Ref. refcite5,6 is used a distorted stereographic projection of hyperconical universes, which are 4D hypersurfaces embedded into 5D Minkowski spacetime. This model is a perturbed flat FRW metric with linear expansion, that takes the perspective of explaining Hubble tension [8]. In our model Universe is a radially inhomogeneous and described by the example of perturbed FLRW metric.

Determination of the Hubble constant using the cosmic distance ladder [9,10] and gravitational lensing [11,12] yields a difference of about 1/10 of its value. We consider the dependence of the change rate of photon energy on the direction of its motion. It is found applying the principle of the extremum of the photon’s energy integral [13,14,15,16,17] and the relationship between it and the energy of a material particle obtained using Lagrange mechanics [13].

2. The Solution of the Einstein Equations for the Perturbed Metric

A generalization of the flat FLRW metric

$$d{s}^2=c^{2}d{\acute{t}}^2-a{\left(\acute{t}\right)}^{2}d{r}^2-a{\left(\acute{t}\right)}^2{r}^2(d{\theta }^2+{sin}^2\theta d{\phi }^2)$$

(1)

is obtained by transforming the time interval

$$d\acute{t}=dt-f\left(t\right)rdr,$$

(2)

where $f\left(t\right)$ is an additional metric function. The linear element takes the following form:

$$d{s}^2=c^{2}d{t}^2-2cf\left(t\right)rdrdt-(a{\left(t\right)}^2-f{\left(t\right)}^2{r}^2)d{r}^2-a{\left(t\right)}^2{r}^2(d{\theta }^2+\mathrm{si}{\mathrm{n}}^2\theta d{\phi }^2).$$

(3)

For small r, this metric is approximated by the metric of flat expanding spacetime. The components of the Einstein tensor are given by

$$G^{11}={\displaystyle \frac{1}{a^6}}\left\{3a^4{\dot{a}}^2-6a^{3}f\dot{a}+r^2\left[2(a^2{f}^2-a^{3}f\dot{a})\dot{f}+2a{f}^3\dot{a}-4a^2{f^2\dot{a}}^2\right]+\right.$$

$$\left.r^4\left[{f^4\dot{a}}^2+2a{f}^3\dot{a}\dot{f}\right]\right\},$$

(4)

$$G^{12}=G^{21}={\displaystyle \frac{fr}{a^6}}[\left(2a\dot{f}\dot{a}+{f\dot{a}}^2\right)r+\left(2a\dot{f}+2f\dot{a}-3a{\dot{a}}^2\right)a],$$

(5)

$$G^{22}={\displaystyle \frac{1}{a^6}}\left\{-2a^3\ddot{a}-a^2{\dot{a}}^2+2a^2\dot{f}+2af\dot{a}+r^2\left(2af\dot{a}\dot{f}+f^2{\dot{a}}^2\right)\right\},$$

(6)

$$G^{33}=G^{44}{\sin }^2\theta ={\displaystyle \frac{1}{{a^{6}r}^2}}\left\{-2a^3\ddot{a}-a^2{\dot{a}}^2+2a^2\dot{f}+2af\dot{a}+\right.$$

$$\left.r^2\left(a\ddot{a}{f}^2-f^2{\dot{a}}^2+f{a}^2\ddot{f}+a^2{\dot{f}}^2+af\dot{a}\dot{f}\right)\right\},$$

(7)

where ($\dot{}$) denotes differentiation with respect to $x^1=ct$.

The components of the contravariant energy-momentum tensor correspond to the energy density, momentum density, energy flux through unit surfaces per unit time and stresses [18]. The energy-momentum tensor of matter in contravariant form is given by

$$T^{ij}=\left(c^2\mathsf{\rho }+p\right)u^i{u}^j-g^{ij}p.$$

(8)

The non-zero contravariant coefficients for the metric (3) will be as follows:

$$g^{11}={\displaystyle \frac{a^2-f^2{r}^2}{a^2}},g^{12}=g^{21}=-{\displaystyle \frac{fr}{a^2}},g^{22}=-{\displaystyle \frac{1}{a^2}},$$

(9)

$$g^{33}=-{\displaystyle \frac{1}{a^2{r}^2}},g^{44}=-{\displaystyle \frac{1}{a^2{r}^2{sin}^2\theta }}.$$

(10)

In the absence of angular momentum of matter ($u^3{=u}^4=0$), the non-zero components of the energy-momentum tensor (8) are given by

$$T^{11}=\left(c^2\mathsf{\rho }+p\right){\left(u^1\right)}^2-\left(1-{\displaystyle \frac{f^2{r}^2}{a^2}}\right)p,$$

(11)

$$T^{12}=T^{21}=\left(c^2\mathsf{\rho }+p\right)u^1{u}^2+{\displaystyle \frac{fr}{a^2}}p,$$

(12)

$$T^{22}=\left(c^2\mathsf{\rho }+p\right){\left(u^2\right)}^2+{\displaystyle \frac{1}{a^2}}p,$$

(13)

$$T^{33}=T^{44}{\sin }^2\theta ={\displaystyle \frac{1}{a^2{r}^2}}p.$$

(14)

Let’s consider a model with a cosmological constant $\mathsf{\Lambda }=0$. We will seek a solution for a region where the radial velocity $u^2$ is small, and the term of component $T^{22}$ containing its square is negligible and can be disregarded. Next, we will justify this assumption. Also, we assume that the terms of the Einstein tensor components (4)-(7) containing $r^2$ in the numerator are small and can be neglected. Given the satisfaction of this condition and for $g^{11}$, Einstein equations $G^{ij}=\chi T^{ij}+g^{ij}\mathsf{\Lambda }$ with the energy-momentum tensor components (11)-(14), yield

$${\displaystyle \frac{3a{\dot{a}}^2-6f\dot{a}}{a^3}}={\chi c}^2\mathsf{\rho },$$

(15)

$${\displaystyle \frac{2a\dot{f}+2f\dot{a}-3a{\dot{a}}^2}{a^3}}fr=\chi a^2\left(c^2\mathsf{\rho }+p\right)u^1{u}^2+\chi frp,$$

(16)

$${\displaystyle \frac{-2a^2\ddot{a}-a{\dot{a}}^2+2a\dot{f}+2f\dot{a}}{a^3}}=\chi p.$$

(17)

The absence of additional energy and momentum influx implies the equality of the corresponding densities $T^{12},T^{21}$ to zero, i.e., the right-hand side of equation (16). From its left-hand side we obtain

$$2a\dot{f}+2f\dot{a}-3a{\dot{a}}^2=0.$$

(18)

Solving this equation for f under the condition $a\left(t_0\right)=1$ we find

$$f\left(t\right)={\displaystyle \frac{1}{a}}\left(f_0+{\displaystyle \frac{3}{2}}c{\int }_{c{t}_0}^{ct}{\dot{a}}^{2}adt\right)$$

(19)

with $f\left(t_0\right)=f_0$. Substituting this into equation (15) under the condition

$$\mathsf{\rho }={\displaystyle \frac{{\mathsf{\rho }}_0}{{a\left(t\right)}^3}},$$

(20)

with average density of the Universe at the present time ${\mathsf{\rho }}_0=\mathsf{\rho }\left(t_0\right)$ yields

$$3a{\dot{a}}^2-{\displaystyle \frac{6f_0\dot{a}}{a}}-{\displaystyle \frac{9\dot{a}}{a}}c{\int }_{{ct}_0}^{ct}{\dot{a}}^{2}adt={\chi c}^2{\mathsf{\rho }}_0.$$

(21)

Expressing the integral from here and differentiating it with respect to $x^1$, from the obtained equality we find

$$\chi c^2{\mathsf{\rho }}_0\left(1-{\displaystyle \frac{a\ddot{a}}{{\dot{a}}^2}}\right)=-3{\dot{a}}^{2}a+3\ddot{a}{a}^2.$$

(22)

This equation has a solution

$$a\left(t\right)=e^{c{A}_1\left(t-t_0\right)},$$

(23)

where $A_1$ is a constant.

From equations (17) and (18), we obtain the pressure

$$p=-{\displaystyle \frac{2}{c\chi }}{\displaystyle \frac{d}{dt}}\left({\displaystyle \frac{\dot{a}}{a}}\right),$$

(24)

which, due to (23), turns out to be $p=0$. In this case, from (16) under the condition of the absence of energy and momentum influx from the outside (18), we find the radial component of the matter velocity in the comoving reference frame.

$$u^2=0.$$

(25)

Thus, it turns out to be a valid assumption that the radial velocity in the region under consideration is small and may not be taken into account in the component of the energy-momentum tensor (13).

The function (19), with the scale factor of space (23), takes the form

$$f\left(t\right)=e^{-c{A}_1\left(t-t_0\right)}{f}_0+{\displaystyle \frac{1}{2}}{A}_1\left(e^{2c{A}_1\left(t-t_0\right)}-e^{-{cA}_1\left(t-t_0\right)}\right)$$

(26)

and equation (21) yields

$${{3A}_1}^2-6f_0{A}_1={\chi c}^2{\mathsf{\rho }}_0.$$

(27)

3. Change in Photon Energy

Let’s find how the energy of a photon changes in the considered space-time applying the principle of the extremum of the photon’s energy integral [13,14,15,16,17]. The Euler-Lagrange equations for the covariant components of the energy-momentum vector of a light-like particle of unit energy

$${\displaystyle \frac{d{p}_{\lambda }}{d\mu }}-F_{\lambda }=0$$

(28)

with

$$F_{\lambda }={\displaystyle \frac{1}{2u_1{u}^1}}{\displaystyle \frac{\partial g_{ij}}{\partial x^{\lambda }}}{u}^i{u}^j,$$

(29)

corresponding to the energy, take the form

$${\displaystyle \frac{d{p}_1}{d\mu }}=-{\displaystyle \frac{1}{u_{1ph}{u}_{ph}^1}}\left[\dot{f}r{u}_{ph}^1{u}_{ph}^2+(\dot{a}a-\dot{f}f{r}^2{{\left)\right(u}_{ph}^2)}^2+\dot{a}a{\mathsf{\Omega }}^2\right]$$

(30)

under the designation

$${\mathsf{\Omega }}^2=r^2\left[({\left(u_{ph}^3\right)}^2+\mathrm{si}{\mathrm{n}}^2\theta {\left(u_{ph}^4\right)}^2\right].$$

(31)

From the equality of interval (3) to zero, we obtain

$$u_{ph}^2={\displaystyle \frac{-fr{u}_{ph}^1+\sigma \sqrt{a^2{\left(u_{ph}^1\right)}^2-(a^2-f^2{r}^2){\mathsf{\Omega }}^2}}{a^2-f^2{r}^2}},$$

(32)

where $\sigma$ takes on values of ±1 depending on the direction of the light. For small deviations from radial motion:

$${\mathsf{\Omega }}^2<<{{(u}_{ph}^1)}^2$$

(33)

we find

$$u_{ph}^2={\displaystyle \frac{1}{fr+\sigma a}}{u}_{ph}^1-{\displaystyle \frac{\sigma {\mathsf{\Omega }}^2}{2a{u}_{ph}^1}}.$$

(34)

The expression for the time component of the covariant velocity vector is as follows:

$$u_{1ph}=u_{ph}^1-fr{u}_{ph}^2={\displaystyle \frac{a}{a+\sigma fr}}{u}_{ph}^1+{\displaystyle \frac{\sigma fr{\mathsf{\Omega }}^2}{2a{u}_{ph}^1}}.$$

(35)

Further we will consider $fr,\dot{a}/a,\dot{f}/f$ to be small quantities. After substituting the components $u_{ph}^2$ and $u_{1ph}$ into equation (30) without accounting for higher order small quantities, it is rewritten in the form

$${\displaystyle \frac{d{p}_1}{d\mu }}=-{\displaystyle \frac{\sigma \dot{a}+\dot{f}r}{fr+\sigma a}}.$$

(36)

Since the covariant and contravariant temporal components of the generalized momenta of a light-like particle

$$p_{\lambda }={\displaystyle \frac{u_{\lambda }}{u^1{u}_1}},$$

(37)

$$p^{\lambda }={\displaystyle \frac{u^{\lambda }}{u^1{u}_1}}$$

(38)

are related by

$$p^1={\displaystyle \frac{u_{ph}^1}{u_{1ph}}}{p}_1=\left(1+\sigma {\displaystyle \frac{fr}{a}}\right)p_1,$$

(39)

the rate of energy change can be written in the form

$${\displaystyle \frac{d{p}^1}{d\mu }}=\left(1+\sigma {\displaystyle \frac{fr}{a}}\right){\displaystyle \frac{d{p}_1}{d\mu }}+\sigma {\displaystyle \frac{1}{a^2{u}_{ph}^1}}\left[r(\dot{fa}-\dot{af})u_{ph}^1{+fau}_{ph}^2\right].$$

(40)

Substituting $u_{ph}^2$ from equation (34) and $d{p}_1/d\mu$ from equation (36) into this expression, we find

$${\displaystyle \frac{d{p}^1}{d\mu }}=-{\displaystyle \frac{\dot{f}r+\sigma \dot{a}}{\sigma a}}+{\displaystyle \frac{1}{{\sigma a}^2}}\left[r\left(\dot{fa}-\dot{af}\right)+{\displaystyle \frac{fa}{fr+\sigma a}}-{\displaystyle \frac{\sigma {f\mathsf{\Omega }}^2}{2{\left(u_{ph}^1\right)}^2}}\right].$$

(41)

Assuming that the terms in this equation containing r are small in the region under consideration, we can write

$${\displaystyle \frac{d{p}^1}{d\mu }}=-{\displaystyle \frac{\dot{a}}{a}}+{\displaystyle \frac{f}{{\sigma a}^2}}-{\displaystyle \frac{{f\mathsf{\Omega }}^2}{2a^2{\left(u_{ph}^1\right)}^2}}.$$

(42)

In the case of the motion of light emitted towards the observer, we choose $\sigma =-1$. Then, upon substitution of (23) and (26), this expression is transformed to become

$${\displaystyle \frac{d{p}^1}{d\mu }}={-A}_1\left({\displaystyle \frac{3}{2}}-{\displaystyle \frac{1}{2}}{e}^{-3{cA}_1\left(t-t_0\right)}\right)-{f_{0}e}^{-3c{A}_1\left(t-t_0\right)}-$$

$${\displaystyle \frac{{\mathsf{\Omega }}^2}{2{\left(u_{ph}^1\right)}^2}}\left[{\displaystyle \frac{1}{2}}{A}_1\left(1-e^{-3{cA}_1\left(t-t_0\right)}\right)+f_0{e}^{-3c{A}_1\left(t-t_0\right)}\right].$$

(43)

In the region under consideration, the energy of a radially moving photon is given by

$$E_{ph}=h{\mathsf{\rho }}_0\left(1-r{\left.{\displaystyle \frac{d{p}^1}{d\mu }}\right|}_{t=t_0}\right),$$

(44)

where ${\mathsf{\rho }}_0$ - is the value of the photon frequency at the observation point. Since in the cosmological model described by an orthogonal metric, the Hubble constant is expressed as

$$\breve{H}=c{\displaystyle \frac{\dot{a}\left(t_0\right)}{a\left(t_0\right)}},$$

(45)

from equation (27), when comparing it with the Friedman equation

$$\chi c^2\mathsf{\rho }={\displaystyle \frac{3}{a^2}}\left[k+{\left(\dot{a}\right)}^2\right]-\mathsf{\Lambda },$$

(46)

it follows that $A_1$ can be a quantity of the same order as the Hubble parameter. After substituting the function f (26) without small quantities, equation (43) takes the form

$${\left.{\displaystyle \frac{d{p}^1}{d\mu }}\right|}_{t=t_0}={-A}_1{-f}_0-{\displaystyle \frac{f_0{\mathsf{\Omega }}^2}{2{\left(u_{ph}^1\right)}^2}}.$$

(47)

Since in the region under consideration the spacetime metric (3) approximates the Minkowski metric, the affine parameter will be close to $\mu =ct$.

4. Change in the Energy of a Material Particle in the Co-Moving Frame

The Lagrangian of a material particle [13] in space-time (3) is given by

$$L=cm\left\{{\left(u^1\right)}^2-2f\left(t\right)r{u}^1{u}^2-\right.$$

$${\left.[a{\left(t\right)}^2-f{\left(t\right)}^2{r}^2]{\left(u^2\right)}^2-a{\left(t\right)}^2{r}^2({\left(u^3\right)}^2+\mathrm{si}{\mathrm{n}}^2\theta {\left(u^4\right)}^2)\right\}}^{1/2}.$$

(48)

The motion of a material particle in a gravitational field, like that of a photon, in accordance with equation (28) at $\mu =s$ can be represented [10,14,16] as the result of the action of generalized forces

$$F_{\lambda }={\displaystyle \frac{\partial L}{\partial x^{\lambda }}}={\displaystyle \frac{1}{2}}cm{\displaystyle \frac{\partial g_{ij}}{\partial x^{\lambda }}}{u}^i{u}^j,$$

(49)

which, however, do not generally form a first-order tensor. The components of the vector of generalized forces associated with $F_{\lambda }$,

$$F^k={\displaystyle \frac{1}{2}}cm{g}^{k\lambda }{\displaystyle \frac{\partial g_{ij}}{\partial x^{\lambda }}}{u}^i{u}^j$$

(50)

are related to gravitational forces ${\mathrm{F}}^k=c{F}^k$. The equations of motion (28) for the contravariant momentum are transformed into the form

$${\displaystyle \frac{d{\overline{p}}^k}{d\mu }}={\displaystyle \frac{d{p}^k}{d\mu }}+{\displaystyle \frac{d{\overleftrightarrow{p}}^k}{d\mu }}=F^k,$$

(51)

where

$${\displaystyle \frac{d{\overleftrightarrow{p}}^k}{d\mu }}=g^{k\lambda }{\displaystyle \frac{\partial g_{\lambda i}}{\partial x^j}}{u}^j{p}^i$$

(52)

represents the energy imparted to the gravitational field.

Let’s consider a particle of matter in a comoving frame, in which it is motionless (25). The component of the force vector (50) corresponding to energy, $F_m^1$, becomes zero. The rate of change of energy transferred to the gravitational field (52) will be

$${\displaystyle \frac{d{p}_m^1}{ds}}=cm{\displaystyle \frac{f\dot{f}{r}^2}{a^2}}.$$

(53)

Due to (51), the energy of the material particle changes at a rate of

$${\displaystyle \frac{d{p}_m^1}{ds}}=-cm{\displaystyle \frac{f\dot{f}{r}^2}{a^2}}.$$

(54)

The derivative of the function $f\left(t\right)$ has the form

$$\dot{f}={-A_1{f}_{0}e}^{-c{A}_1\left(t-t_0\right)}+A_1^2\left(e^{2c{A}_1\left(t-t_0\right)}+{\displaystyle \frac{1}{2}}{e}^{-c{A}_1\left(t-t_0\right)}\right).$$

(55)

Now we can express the change in energy of the material particle at the present time as

$${\left.{\displaystyle \frac{d{p}_m^1}{ds}}\right|}_{t=t_0}=cm{A}_1\left(f_0+{\displaystyle \frac{3}{2}}{A}_1\right)f_0{r}^2.$$

(56)

In the considered region, the energy of the material particle will be

$$E_m=c\left(cm-r{\left.{\displaystyle \frac{d{p}_m^1}{ds}}\right|}_{t=t_0}\right).$$

(57)

5. Cosmological Redshift in the Radial Motion of Photons

The equation (27) is rewritten in the form

$$3{\left(A_1+f_0\right)}^2-12f_0{A}_1-3{f_0}^2={\chi c}^2{\mathsf{\rho }}_0.$$

(58)

We assume that $f_0$ is of the same order of magnitude as $A_1$ and, consequently, with the Hubble constant. Since in the considered region the quantity $fr$ is small, taking into account the expressions for the rate of change of particle energies (47) and (56), the change in the energy of a material particle will be small in comparison with the change in the equivalent energy of a photon (44) because of

$$\left({\left.{\displaystyle \frac{d{p}_m^1}{ds}}\right|}_{t=t_0}\right)/\left(h{v}_0{\left.{\displaystyle \frac{d{p}^1}{d\mu }}\right|}_{t=t_0}\right)=o\left(fr\right),$$

(59)

where the differentiation parameters are

$$s\approx \mu \approx ct.$$

(60)

The observed redshift of photon energy is considered as a result of the change in the ratio of their energy (44), (47) to the energy of atoms (57). In a small neighborhood without small higher-order quantities due to (59), this will be

$${\displaystyle \frac{E_{ph}}{E_m}}={\displaystyle \frac{h{v}_0}{c^{2}m}}\left[1-r\left(A_1+f_0\right)\right],$$

(61)

from which the value of the Hubble constant at present time

$$H_0=c{(A}_1{+f}_0),$$

(62)

can be deduced, obtained as a result of registering radially moving photons.

In this case, equation (58) takes the form

$$3{\left({\displaystyle \frac{H_0}{c}}\right)}^2-12{\displaystyle \frac{H_0}{c}}{{f_0+9f}_0}^2={\chi c}^2{\mathsf{\rho }}_0.$$

(63)

6. Observed Length Element

The length element

$$d{l}^2=\left({\displaystyle \frac{g_{1p}{g}_{1q}}{g_{11}}}-g_{pq}\right)d{x}^{p}d{x}^q$$

(64)

with $p,q=1,2,3$ in the space-time (3) is given by

$$dl=a\left(t\right)\sqrt{d{r}^2+d{\theta }^2+{sin}^2\theta d{\phi }^2},$$

(65)

that is, the space in its proper frame is homogeneous. However, the observed space in the coordinate frame appears radially inhomogeneous. The radial elements of the length of the metric (3) at constant time

$$d{l}^{\prime }=\sqrt{a{\left(t\right)}^2-f{\left(t\right)}^2{r}^2}dr$$

(66)

is less than obtained from (65) with $d\theta =d\phi =0$ and approach it with decreasing r.

7. Difference in Hubble Constant Measurements

Let’s consider how the energy of a photon will change in the presence of gravitational lensing. We denote the change in the energy of the photon when moving along the path $L_l$, which includes gravitational lensing, as

$$\delta p_l^1={\int }_{L_l}{\displaystyle \frac{d{p}^1}{d\mu }}d\mu$$

(67)

and along the radial path as

$$\delta p_r^1={\int }_{L_r}{\displaystyle \frac{d{p}^1}{d\mu }}d\mu .$$

(68)

We assume that the coordinate distance $\mathsf{\Delta }r$ from the beginning of the paths to the observer is the same. Under the gravitational lensing, neglecting small quantities, the Hubble constant is determined as follows:

$$H_0=-{\displaystyle \frac{c\left(\delta p_l^1-\delta p_r^1\right)}{\delta {\mu }_l-\delta {\mu }_r}},$$

(69)

where $\delta {\mu }_l$ and $\delta {\mu }_r$ are the change in the affine parameter along each of the paths.

The coordinate system is chosen such that when the photon moves along the path $L_l$ the condition $\theta ={\displaystyle \frac{\pi }{2}}$ is satisfied. If we take into account only the influence of the deviation from radial motion, then the difference in the time taken for the photon to travel along both paths will be

$$\delta t={\displaystyle \frac{1}{c}}{\int }_0^{\mathsf{\Delta }r}(1-cos\left(\alpha r\right))dr,$$

(70)

where in the first approximation we consider

$$\cos \left(\left(r\right)\right)=-{\displaystyle \frac{u_{ph}^2}{\sqrt{{\left(u_{ph}^2\right)}^2+{\left(r{u}_{ph}^4\right)}^2}}}.$$

(71)

Due to condition (33), the integral (70) transforms into the form

$$\delta t={\displaystyle \frac{1}{c}}{\int }_0^{\mathsf{\Delta }r}{\displaystyle \frac{{\left(r{u}_{ph}^4\right)}^2}{2{\left(u_{ph}^2\right)}^2}}dr.$$

(72)

In the considered region, we have

$$\delta {\mu }_l-\delta {\mu }_r=c\delta t.$$

(73)

By substituting the corresponding value from (47) into (67) and (68), we find the additional difference in energy change caused by the discrepancy in the rate of its during radial and angular motion

$$\mathsf{\Delta }(\delta p_l^1-\delta p_r^1)=-{\int }_0^{\mathsf{\Delta }r}{\displaystyle \frac{f_0{\left(r{u}_{ph}^4\right)}^2}{2{\left(u_{ph}^1\right)}^2}}dr.$$

(74)

The change in the value of the Hubble constant when light moves along a curved path will be

$$\mathsf{\Delta }{H}_0=-{\displaystyle \frac{c\mathsf{\Delta }\left(\delta p_l^1-\delta p_r^1\right)}{\delta {\mu }_l-\delta {\mu }_r}}.$$

(75)

Since for the considered motion of the photon, the relations

$$u_{ph}^1\approx -a{u}_{ph}^2\approx {-u}_{ph}^2$$

(76)

hold, by substituting the values (73) and (9) into (75), we obtain

$$\mathsf{\Delta }{H}_0={cf}_0.$$

(77)

Let’s estimate the values of $A_1$ and $f_0$ using the results of obtaining the Hubble constant with and without gravitational lensing. The Hubble constant was found through gravitational lensing of the cosmic microwave background: 68.3±1.5 $\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$ [11] and of radiation from the supernovae: ${66.6}_{-3.3}^{+4.1}\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$[12]. The values obtained through the extragalactic distance ladder to supernovae are 73.5±1.1 $\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$[9] and ${75.4}_{-3.7}^{+3.8}$$\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$[10]. These results yield the value $\mathsf{\Delta }{H}_0={cf}_0={-7}_{-2.6}^{+2.4}$$\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$, and due to relation (62) corresponding to the radial motion of the photon we obtain $c{A}_1=81.5\pm 2.4$$\mathrm{km}{\mathrm{s}}^{-1}{\mathrm{Mpc}}^{-1}$.

When the value of $f_0$ is negative, the photon moving towards observer from the star additionally loses energy compared to the motion tangential to the circle from the center of observation or in the opposite direction. The dependence of the spacetime properties on the observer’s position is similar to the manifestation of the spacetime curvature in the non-flat FLRW model. In view of (63), the negative value of the coefficient $f_0$ excludes the possibility of explaining the presence of dark energy using this perturbed metric.

8. Conclusions

In the considered model in region with small perturbations of the FLRW metric, while conserving energy and momentum, this part of the Universe expands exponentially with zero pressure. Solution is obtained that relates the constant of relative expansion of space for a model with a zero cosmological constant. According to this model, in addition to the actual expansion of space, its additional movement in the radial direction is also established. The considered metric can also be applied to a cosmological model with a nonzero $\mathsf{\Lambda }$ -term.

The redshift in the spectra of galaxies is found using Lagrange mechanics in its application to the principle of the extremal integral of photon energy. The light used to obtain redshift through the ladder distance travels along a radial path, unlike gravitational lensing, where its pass has a tangential component. In the considered spacetime, from the observer’s standpoint, a photon loses energy faster when moving radially towards them compared to its motion with a tangential component. Based on the difference in value of the Hubble constant obtained using gravitational lensing and distance ladder the metric coefficients of considered space-time were found in a local region in which their perturbations are small. The value of these coefficients corresponds to the additional movement of space away from the observer.

This perturbation model explains the Hubble tension but does not come close to explaining the existence of dark energy. Certain parameters correspond to a region with low redshift.