*Communication* **Cosmological Model with Interconnection between Dark Energy and Matter**

**Gennady S. Bisnovatyi-Kogan**


**Abstract:** It is accepted in the present cosmology model that the scalar field, which is responsible for the inflation stage in the early universe, transforms completely into matter, and the accelerated universe expansion is presently governed by dark energy (DE), whose origin is not connected with the inflationary scalar field. We suppose here that dark matter (DM) has a common origin with a small variable component of dark energy (DEV). We suggest that DE may presently have two components, one of which is the Einstein constant Λ, and another, smaller component DEV (Λ*V*) comes from the remnants of the scalar field responsible for inflation, which gave birth to the origin of presently existing matter. In this note we consider only the stages of the universe expansion after recombination, *z* 1100, when DM was the most abundant component of the matter, therefore we suggest for simplicity that a connection exists between DM and DEV so that the ratio of their densities remains constant over all the stages after recombination, *ρDM* = *αρDEV*, with a constant *α*. One of the problems revealed recently in cosmology is a so-called Hubble tension (HT), which is the difference between values of the present Hubble constant, measured by observation of the universe at redshift *z* - 1, and by observations of a distant universe with CMB fluctuations originated at *z* ∼ 1100. In this paper we suggest that this discrepancy may be explained by deviation of the cosmological expansion from a standard Lambda-CDM model of a flat universe, due to the action of an additional variable component DEV. Taking into account the influence of DEV on the universe's expansion, we find the value of *α* that could remove the HT problem. In order to maintain the almost constant DEV/DM energy density ratio during the time interval at *z* < 1100, we suggest the existence of a wide mass DM particle distribution.

**Keywords:** dark energy; dark matter; Hubble constant

#### **1. Introduction**

During his long and successful scientific career, Yu. N. Gnedin worked also on problems on the border of physics and cosmology, namely on the possibility of direct searching, in astronomical observations, for axions, belonging to the family of Goldston bosons. These tiny particles, introduced by theoretical physicists, have been proposed as candidates for dark matter, whose presence was necessary for interpretation of different astronomical observations. During the years 1992–2009, Yu. N. Gnedin published at least six papers on this intriguing topic [1–6]. The cosmological model considered in this paper demands the presence of very light particles in DM, and axions could be the best candidates for that.

Here we consider a cosmological model that differs slightly from the widely accepted Λ*CDM* model by having two components of the dark energy instead of the usually considered one component represented by Λ. The model considered here is constructed in order to explain the so called "Hubble Tension", which is the observed discrepancy between different averaged "local" measurements of the Hubble constant from one side, and distant measurements based on the analysis of the CMB fluctuations, connected with

**Citation:** Bisnovatyi-Kogan, G.S. Cosmological Model with Interconnection between Dark Energy and Matter. *Universe* **2021** , *7* , 412. https://doi.org/10.3390/ universe7110412

Academic Editors: Galina L. Klimchitskaya, Vladimir M. Mostepanenko and Nazar R. Ikhsanov

Received: 11 October 2021 Accepted: 27 October 2021 Published: 29 October 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

the recombination epoch at redshift *z* ∼ 1100. In order to compare these two measurements, it was necessary to specify the cosmological model for determining the present value of the Hubble constant from distant measurements. We accept that the origin of HT is connected with the not fully appropriate Λ*CDM* model, and suggest its modification removing this discrepancy.

While the present value of the Hubble constant (HC) is one of the most important cosmological parameters, its measurements for many years have been performed by different astronomical groups ("local measurements"). These measurements have been based on several steps in order to have a possibility for independent distance measurements of more and more distant objects with observed runaway velocities. The use of different steps in this sequence by different groups has resulted in finding different values for the HC, from ∼50 km/s/Mps by the Sandage–Tamman group to ∼100 km/s/Mps by the de Vacouleurs group [7]. Over time the local HC measurements at redshift *z* - 1 have been substantially improved due to construction of big telescopes and measurements by the Hubble mission, which permitted the narrowing of this interval of the HC to 72–75 km/s/Mps.

The measurements of the CMB fluctuations by instruments in the satellites WMAP and PANCK gave the possibility for fully independent HC measurements at the recombination epoch. HT appears when you compare the present values of the HC obtained from CMB using a simple Λ*CDM* model with local measurements. It is claimed that the discrepancy in these two values is statistically significant in the range 4.5*σ* to 6.3*σ* [8]; see, nevertheless [9].

It is accepted in the present cosmology model that the scalar field, which is responsible for the inflation stage in the early universe, transforms completely into matter, and the accelerated universe expansion is presently governed by dark energy (DE), whose origin is not connected with the inflationary scalar field. We suppose here that dark matter (DM) has a common origin with a small variable component of dark energy (DEV). We suggest that DE presently may have two components, one of which is the Einstein constant Λ, and another, smaller component DEV (Λ*V*) comes from the remnants of the scalar field responsible for inflation, which gave birth to the origin of presently existing matter. In this note we consider only the stages of the universe's expansion after recombination, *z* 1100, when DM is the most abundant component of the matter; therefore, we suggest for simplicity that a connection exists between DM and DEV so that the ratio of their densities remains constant over all the stages after recombination *ρDM* = *αρDEV*, with a constant *α*.

In this paper we suggest that this discrepancy may be explained by the deviation of the cosmological expansion from a standard Lambda-CDM model of a flat universe, due to the action of an additional variable component DEV. Taking into account the influence of DEV on the universe's expansion, we find the value of *α* that could remove the HT problem. In order to maintain the almost constant DEV/DM energy density ratio during the time interval at *z* < 1100, we suggest the existence of a wide mass DM particle distribution.

#### **2. Universe with Common Origin of DM and DE**

The scalar field with the potential *V*(*φ*), where *φ* is the intensity of the scalar field, is considered the main reason for the inflation [10,11], but see [12]. The equation for the scalar field in the expanding universe is written as [13]

$$
\ddot{\phi} + 3\frac{d}{a}\phi = -\frac{dV}{d\phi}.\tag{1}
$$

Here *a* is a scale factor in the flat expanding universe [7]. The density *ρ<sup>V</sup>* and pressure *PV* of the scalar field <sup>1</sup> are defined as [13]

$$
\rho\_V = \frac{\phi^2}{2} + V, \quad P\_V = \frac{\phi^2}{2} - V. \tag{2}
$$

Consider the universe with the initial scalar field, at initial intensity *φin* and initial potential *Vin*, and at zero derivative *φ*˙ *in* = 0. The derivative of the scalar field intensity is growing on the initial stage of inflation.

Let us suggest that after reaching the relation

$$
\dot{\phi}^2 = \text{2a}V,\tag{3}
$$

it is preserved during further expansion. The kinetic part of the scalar field is transforming into matter, presumably dark matter, and the constant *α* determines the ratio of the the dark energy density, represented by *V*, to the matter density, represented by the kinetic term. As follows from observations, the main part of DE is represented presently by DE, which may be considered as the Einstein constant Λ. At earlier times the input of constant Λ is smaller than the input of Λ*V*, for a wide interval of constant *α* values.

Let us consider an expanding flat universe, described by the Friedmann equation [7]

$$\frac{d^2}{d^2} = \frac{8\pi G}{3}\rho + \frac{\Lambda}{3}.\tag{4}$$

Introduce

$$
\rho\_{\Phi} = V, \quad P\_{\Phi} = -V, \quad \rho\_{m} = \frac{\phi^{2}}{2}, \quad P\_{m} = \beta \frac{\phi^{2}}{2}, \quad \text{with} \quad \rho\_{m} = a\rho\_{\Phi}. \tag{5}
$$

We suggest that only part *β* of the kinetic term makes the input into the pressure of the matter, so it follows from (3) and (5)

$$
\rho = \rho\_{\phi} + \rho\_{m} = (1 + a)V, \quad P = P\_{\phi} + P\_{m} = -(1 - a\beta)V. \tag{6}
$$

The adiabatic condition

$$\frac{d\rho}{d\rho + P} = -\frac{d\mathcal{V}}{\mathcal{V}} = -\mathfrak{Z}\frac{da}{a}, \quad \mathcal{V} \quad \text{is the volume,} \tag{7}$$

may be written as

$$\dot{\rho} = -3\frac{\dot{a}}{a}(\rho\_{\phi} + \rho\_{m} + P\_{\phi} + P\_{m}) = -3\frac{\dot{a}}{a}(\rho\_{m} + P\_{m}) = -3a\frac{1+\beta}{1+a}\frac{\dot{a}}{a}\rho. \tag{8}$$

$$\frac{\rho}{\rho\_\*} = \left(\frac{a}{a\_\*}\right)^{\frac{3a}{2(1+a)}} \quad \text{for} \quad \beta = 0. \tag{9}$$

#### *2.1. A Universe with* Λ *= 0*

Suggest first that cosmological constant Λ = 0, and DE is created only by the part of the scalar field Λ*V*, represented by *V*. The expressions for the total density *ρ*, scaling factor *a*, and Hubble "constant" *H* follow from (4)–(8) as

$$\frac{a}{a\_\*} = (6\pi G\rho\_\* t^2)^{\frac{1+a}{3a(1+\beta)}} \left(\frac{a(1+\beta)}{1+a}\right)^{\frac{2(1+a)}{3a(1+\beta)}} = \left(\frac{\rho\_\*}{\rho}\right)^{\frac{1+a}{3a(1+\beta)}} = \left(\frac{t}{t\_\*}\right)^{\frac{2(1+a)}{3a(1+\beta)}},$$

$$\rho = \left(\frac{1+a}{a(1+\beta)}\right)^2 \frac{1}{6\pi G t^2}, \quad H = \frac{a}{a} = \frac{2(1+a)}{3a(1+\beta)t}.\tag{10}$$

Here *ρ*<sup>∗</sup> = *ρ*(*t*∗), *a*<sup>∗</sup> = *a*(*t*∗), *t*<sup>∗</sup> is an arbitrary time moment. Write the expressions for particular cases. For *β* = 1/3 (radiation dominated universe) it follows from (10)

$$\frac{a}{a\_\*} = (6\pi \mathcal{G}\rho\_\* t^2)^{\frac{1+a}{4a}} \left[\frac{4\alpha}{3(1+a)}\right]^{\frac{1+a}{2a}} = \left(\frac{\rho\_\*}{\rho}\right)^{\frac{1+a}{4a}} = \left(\frac{t}{t\_\*}\right)^{\frac{1+a}{2a}},$$

$$\rho = \left(\frac{3(1+a)}{4a}\right)^2 \frac{1}{6\pi \mathcal{G}t^2}, \quad H = \frac{a}{a} = \frac{1+a}{2at}.\tag{11}$$

For the value of *β* = 0 (dusty universe, *z* < 1100) we have

$$\frac{a}{a\_\*} = \left(6\pi G\rho\_\* t^2\right)^{\frac{1+a}{3a}} \left[\frac{a}{1+a}\right]^{\frac{2(1+a)}{3a}} = \left(\frac{\rho\_\*}{\rho}\right)^{\frac{1+a}{3a}} = \left(\frac{t}{t\_\*}\right)^{\frac{2(1+a)}{3a}},$$

$$\rho = \left(\frac{1+a}{a}\right)^2 \frac{1}{6\pi G t^2}, \quad H = \frac{\dot{a}}{a} = \frac{2(1+a)}{3at}.\tag{12}$$

#### *2.2. A Universe in the Presence of the Cosmological Constant* Λ

Equations (5)–(8) are valid in the presence of Λ. The solution of Equation (4) with nonzero Λ is written in the form

$$\left(\frac{a}{a\_\*}\right)^{\frac{3\alpha(1+\beta)}{2(1+a)}} = \sqrt{\frac{8\pi G\rho\_\*}{\Lambda c^2}} \sinh\left(\sqrt{\frac{\Lambda}{3}} \frac{3\alpha(1+\beta)}{2(1+a)}ct\right) = \sqrt{\frac{\rho\_\*}{\rho}};\tag{13}$$

$$\sqrt{\frac{\Lambda c^2}{8\pi G\rho}} = \sinh\left(\sqrt{\frac{\Lambda}{3}}\frac{3a(1+\beta)}{2(1+a)}ct\right), \quad H = \frac{d}{a} = \sqrt{\frac{\Lambda c^2}{3}}\coth\left(\sqrt{\frac{\Lambda}{3}}\frac{3a(1+\beta)}{2(1+a)}ct\right). \tag{14}$$

For the dusty universe (*β* = 0), after recombination at *z* < 1100 we have

$$\left(\frac{a}{a\_\*}\right)^{\frac{3\alpha}{2(1+a)}} = \sqrt{\frac{8\pi G\rho\_\*}{\Lambda c^2}} \sinh\left(\sqrt{\frac{\Lambda}{3}} \frac{3\alpha}{2(1+a)}ct\right) = \sqrt{\frac{\rho\_\*}{\rho}};\tag{15}$$

$$\sqrt{\frac{\Lambda c^2}{8\pi G\rho}} = \sinh\left(\sqrt{\frac{\Lambda}{3}}\frac{3a}{2(1+a)}ct\right), \quad H = \frac{d}{a} = \sqrt{\frac{\Lambda c^2}{3}}\coth\left(\sqrt{\frac{\Lambda}{3}}\frac{3a}{2(1+a)}ct\right). \tag{16}$$

The dusty universe without DEV is described by relations following from (15) and (16) at *α* → ∞, giving

$$\left(\frac{a}{a\_\*}\right)^{\frac{3}{2}} = \sqrt{\frac{8\pi G\rho\_\*}{\Lambda c^2}} \sinh\left(\sqrt{\frac{\Lambda}{3}}\frac{3}{2}ct\right) = \sqrt{\frac{\rho\_\*}{\rho}};\tag{17}$$

$$\sqrt{\frac{\Lambda c^2}{8\pi G\rho}} = \sinh\left(\sqrt{\frac{\Lambda}{3}}\frac{3}{2}ct\right), \qquad H = \frac{d}{a} = \sqrt{\frac{\Lambda c^2}{3}}\coth\left(\sqrt{\frac{\Lambda}{3}}\frac{3}{2}ct\right). \tag{18}$$

#### **3. Hubble Tension**

Recently, a challenge in cosmology was formulated because of different values, obtained from different experiments, of the Hubble constant at the present epoch. There is a significant discrepancy (tension) between the *Planck* measurement from cosmic microwave background (CMB) anisotropy, where the best-fit model gives [14–16],

$$H\_0^{\rm Pl8} = 67.36 \pm 0.54 \text{ km s}^{-1} \text{Mpc}^{-1},\tag{19}$$

and measurements using type Ia supernovae (SNIa) calibrated with Cepheid distances [17–21],

$$H\_0^{\mathbb{R}19} = 74.03 \pm 1.42 \text{ km s}^{-1} \text{Mpc}^{-1}. \tag{20}$$

Measurements using time delays from lensed quasars [22] gave the value *H*<sup>0</sup> = 73.3+1.7 <sup>−</sup>1.8 km s−<sup>1</sup> Mpc<sup>−</sup>1, while in [23] it was found that

*<sup>H</sup>*<sup>0</sup> <sup>=</sup> 72.4 <sup>±</sup> 1.9 km s−<sup>1</sup> Mpc−<sup>1</sup> using the tip of the red giant branch applied to SNIa, which is independent of the Cepheid distance scale. Analysis of a compilation of these and other recent high- and low-redshift measurements shows [24] that the discrepancy between *Planck* [16] and any three independent late-Universe measurements is between 4*σ* and 6*σ*. Different sophisticated explanations for the appearance of HT have been proposed [25–30] (see also [31–34]) and new experiments have been proposed for checking the reliability of this tension [35] (see also review [36]).

Dark matter (DM) and dark energy (DE) represent about 96% of the universe constituents [14,17,18], but their origin is still not clear. The present value of DE density may be represented by the Einstein cosmological constant Λ [37], but may also be a result of the action of the Higgs-type scalar fields, which are supposed to be the reason for the inflation in the early universe [10] (see also [11,12,38]). The value of the induced Λ, suggested for the inflation, is many orders of magnitude larger than its present value, and no attempts have been made to find a connection between them. The origin of DM is even more vague. There are numerous suggestions for its origin [39–41], but none of these possibilities has been experimentally or observationally confirmed, while many of them have been disproved.

To explain the origin of the Hubble Tension, we introduce a variable part of the cosmological "constant" Λ*V*, proportional to the matter density *ρDM* = *αρDEV*. This part of Λ*<sup>V</sup>* influences the cosmological expansion at large redshifts, where the influence of the real Einstein constant Λ is negligible. The value of Λ*<sup>V</sup>* is represented by a small component of DE, which we define as DEV. We suppose here, without knowledge of the physical properties of DM particles, that there is a wide spectrum of DM particle, which could be produced by DEV until present time.

This seems necessary because at decreasing of the DEV field strength in the expanding universe it would be able to make mutual transformations only with DM particles of decreasing mass. The existence of particles with a very low rest mass (axions [42]) is considered often as a candidate for DM. Note that in the paper of Yu N. Gnedin [5], mutual transformations between axions and electromagnetic photons have been considered, instead of the hypothetical "scalar field" in our model.

#### **4. Removing the Hubble Tension**

We consider a model of the expanding universe after recombination, at *z* ≤ 1100, with a fixed ratio *α* of energy densities between DEV, connected with a scalar field, cosmological constant Λ, and DM. If the mass spectrum of DM particles prolongs to very small masses, then we may expect an almost constant DEV–DM ratio. In the inflation model of the universe, only a scalar field was born at the very beginning, and matter was created in the process of expansion from the dynamic part of the scalar field density.

Here we show that in presence of DEV the Hubble value *H* is decreasing with time slower than without it. This creates a larger present value of *H*0, removing the Hubble tension at *α* ∼ 140. The main idea of removing the tension is the following. The CMB measurements give the value of the Hubble constant *Hr*, at the redshift *z* ∼ 1100, close to the moment of recombination. This value is used for calculation of the present value of *H*0.

For analysis of the Hubble tension it is more convenient to use logarithmic variables, so that from (10), (19) and (20) we have

$$
\log H\_0^{\rm Pl8} = \log 67.36 = 1.83; \quad \log H\_0^{\rm RI9} = \log 74.03 = 1.87;
$$

$$
\log \frac{H\_0^{\rm RI9}}{H\_0^{\rm Pl8}} = \Delta \log H\_0 = 0.04. \tag{21}
$$

The Planck value *H<sup>P</sup> <sup>r</sup>* was measured at the moment of recombination *zr* ≈ 1100, and extrapolated to present time using a dusty flat Friedmann model with a very small input from the cosmological constant Λ, when the values *z*, *ω*, *a* are connected as [7]

$$z + 1 = \frac{\omega}{\omega\_0} = \frac{a\_0}{a}.\tag{22}$$

In the case of an equipartition universe the extrapolation should be performed using Equation (16). Numerical modeling of large-scale structure formation gives preference to the cold dark matter model, corresponding to *Pm* ≈ 0, *β* = 0. In our interpretation the Hubble tension is connected with incorrect extrapolation by Equation (16) without DEV. From the value of Hubble tension in Equation (21) we may estimate *α* from the condition that both measurements are correct, but the reported value *HP*<sup>18</sup> <sup>0</sup> is coming from the incorrect extrapolation, and the actual present epoch value of the Hubble constant is determined by *HR*<sup>19</sup> <sup>0</sup> .

From Equations (16) and (22) we obtain for a dusty universe, at *β* = 0, the following connection between the recombination redshift *zr* ≈ 1100, the present age of the universe *t*<sup>0</sup> and the age of the universe *tr*, corresponding to the recombination, in the form

$$z\_{n1} + 1 = \frac{a\_\*}{a\_{n1}} = \left[ \sqrt{\frac{8\pi G\rho\_\*}{\Lambda c^2}} \sinh\left(\sqrt{\frac{\Lambda}{3}} \frac{3a}{2(1+a)} ct\_{n1}\right) \right]^{-\frac{2(1+a)}{3a}} = \left(\frac{\rho\_\*}{\rho\_{ra}}\right)^{-\frac{(1+a)}{3a}},\tag{23}$$

$$H = \frac{d}{a} = \sqrt{\frac{\Lambda c^2}{3}} \coth\left(\sqrt{\frac{\Lambda}{3}} \frac{3\alpha}{2(1+\alpha)} c t\_{ra}\right). \tag{24}$$

Here indices with *α* indicate the values in the universe with DEV. The index "0" is related to present time at *z*<sup>0</sup> = 0. The index "*r*" is related to the moment of recombination. The values in the universe without DEV do not contain *α* in the indices, and are written as

$$z\_r + 1 = \frac{a\_\*}{a\_r} = \left[\sqrt{\frac{8\pi G\rho\_\*}{\Lambda c^2}} \sinh\left(\sqrt{\frac{\Lambda}{3}} \frac{3}{2} c t\_r\right)\right]^{-\frac{2}{3}} = \left(\frac{\rho\_\*}{\rho\_r}\right)^{-\frac{1}{3}}\,\tag{25}$$

$$H = \frac{d}{a} = \sqrt{\frac{\Lambda c^2}{3}} \coth\left(\sqrt{\frac{\Lambda}{3}} \frac{3}{2} ct\_r\right). \tag{26}$$

The observational data are connected with a redshift, so we find the expression for the time as a function of the redshift in the form

$$ct\_{r\mathfrak{a}} = \frac{2(1+a)}{3a} \sqrt{\frac{3}{\Lambda}} \sinh^{-1} \left[ \sqrt{\frac{\Lambda c^2}{8\pi G\rho\_\*}} (z\_{r\mathfrak{a}} + 1)^{\frac{3a}{2(1+a)}} \right] \tag{27}$$

$$\approx \frac{a+1}{a\sqrt{6\pi G\rho\_\*}} (z\_{r\mathfrak{a}} + 1)^{\frac{3a}{2(1+a)}},$$

$$ct\_{0\mathfrak{a}} = \frac{2(1+a)}{3a} \sqrt{\frac{3}{\Lambda}} \sinh^{-1} \left( \sqrt{\frac{\Lambda c^2}{8\pi G\rho\_\*}} \right). \tag{28}$$

At *t* = *tr<sup>α</sup>* the argument of "sinh" is very small, so only the first term in the expansion remains. Corresponding values for the universe without DEV are written as

$$ct\_r = \frac{2}{3} \sqrt{\frac{3}{\Lambda}} \sinh^{-1} \left[ \sqrt{\frac{\Lambda c^2}{8\pi G \rho\_\*}} (z\_r + 1)^{\frac{3}{2}} \right] \approx \frac{1}{\sqrt{6\pi G \rho\_\*}} (z\_r + 1)^{\frac{3}{2}},\tag{29}$$

$$ct\_{0\alpha} = \frac{2(1+\alpha)}{3\alpha} \sqrt{\frac{3}{\Lambda}} \sinh^{-1} \left( \sqrt{\frac{\Lambda c^2}{8\pi G \rho\_\*}} \right). \tag{30}$$

For *zr<sup>α</sup>* = *zr* ≡ *zrec* we obtain a connection between *tr<sup>α</sup>* and *tr* as

$$t\_{\rm ra} = \frac{a+1}{a} (z\_{\rm rec} + 1)^{\frac{3}{2(1+a)}} t\_r. \tag{31}$$

The ratio of Hubble constants at present time *t*<sup>0</sup> to its value at the recombination time *trα*, in the presence of DEV, using (16), is written as

$$\frac{H\_{0\mu}}{H\_{1\mu}} = \frac{\coth\left(\sqrt{\frac{\Lambda}{3}} \frac{3a}{2(1+a)} c t\_0\right)}{\coth\left(\sqrt{\frac{\Lambda}{3}} \frac{3a c t\_{\text{fl}}}{2(1+a)}\right)} = \frac{\tanh\left(\sqrt{\frac{\Lambda}{3}} \frac{3a c t\_{\text{fl}}}{2(1+a)}\right)}{\tanh\left(\sqrt{\frac{\Lambda}{3}} \frac{3a}{2(1+a)} c t\_0\right)} = \frac{\sqrt{\frac{\Lambda}{3}} \frac{3a}{2(1+a)} c t\_{\text{fl}}}{\tanh\left(\sqrt{\frac{\Lambda}{3}} \frac{3a}{2(1+a)} c t\_0\right)}.\tag{32}$$

The corresponding values for the case without DEV, at *α* → ∞, with recombination time *tr* are written as

$$\frac{H\_0}{H\_r} = \frac{\coth\left(\sqrt{\frac{\Delta}{3}} \frac{3}{2} ct\_0\right)}{\coth\left(\sqrt{\frac{\Delta}{3}} \frac{3ct\_r}{2}\right)} = \frac{\tanh\left(\sqrt{\frac{\Delta}{3}} \frac{3ct\_r}{2}\right)}{\tanh\left(\sqrt{\frac{\Delta}{3}} \frac{3}{2} ct\_0\right)} = \frac{\sqrt{\frac{\Delta}{3}} \frac{3}{2} ct\_r}{\tanh\left(\sqrt{\frac{\Delta}{3}} \frac{3}{2} ct\_0\right)}.\tag{33}$$

Using Equations (31)–(33), we obtain the ratio between the correct value *H*0*α*, identified with the local measurement of *H*, and the value *H*0, obtained by calculations without account of DEV, in the form

$$\frac{H\_{0\alpha}}{H\_0} = \frac{\tanh\left(\sqrt{\frac{\Lambda}{3}}\frac{3}{2}ct\_0\right)}{\tanh\left(\sqrt{\frac{\Lambda}{3}}\frac{3\alpha}{2(1+\alpha)}ct\_0\right)}(z\_{\text{rec}}+1)^{\frac{3}{2(1+\alpha)}}.\tag{34}$$

#### *Numerical Estimations*

Let us consider the following commonly accepted, averaged parameters of the universe, used in relation to the HT explanation.

Local Hubble constant *Hl* = 73 km/s/Mpc.

Distantly measured Hubble constant *Hd* = 67.5 km/s/Mpc.

Total density of the flat universe *<sup>ρ</sup>tot* = <sup>2</sup> · <sup>10</sup>−29*h*<sup>2</sup> = 1.066 · <sup>10</sup>−<sup>29</sup> g/cm3, where h = H/100 km/s/Mpc, with *h* = 0.73.

Distantlymeasured densities of the universe components, *<sup>ρ</sup>*Λ<sup>∗</sup> <sup>=</sup> 0.7*ρtot* <sup>=</sup> 7.5 · <sup>10</sup>−<sup>30</sup> g/cm3; *ρm*<sup>∗</sup> = 0.3*ρtot*.

Locally measured cosmological constant density [18] *ρ*Λ<sup>0</sup> = (0.44 ÷ 0.96) *ρtot*, (2*σ* statistics).

Λ = <sup>8</sup>*πGρ*Λ<sup>∗</sup> *<sup>c</sup>*<sup>2</sup> <sup>=</sup> 1.40 · <sup>10</sup>−<sup>56</sup> cm<sup>−</sup>2.

Average age of the universe *<sup>t</sup>*<sup>0</sup> = 4.35 · <sup>10</sup><sup>17</sup> s [43].

The ratio of two Hubble constants is *HTr* = *Hl Hd* = 1.08.

Identifying *H*0*<sup>α</sup>* ≡ *Hl*, *H*<sup>0</sup> ≡ *Hd*, we obtain from (34)

*H*0*<sup>α</sup> H*<sup>0</sup> = 1.08 = tanh- Λ 3 3 <sup>2</sup> *ct*<sup>0</sup> tanh- Λ 3 3*α* <sup>2</sup>(1+*α*) *ct*<sup>0</sup> (*zrec* <sup>+</sup> <sup>1</sup>) <sup>3</sup> <sup>2</sup>(1+*α*) (35)

$$\eta = \frac{\tanh(1.337)}{\tanh\left(1.337 \frac{\alpha}{1+a}\right)} (1101)^{\frac{3}{2(1+a)}} = \frac{0.87095}{\tanh\left(1.337 \frac{\alpha}{1+a}\right)} (1101)^{\frac{3}{2(1+a)}} \dots$$

Finally we obtain the equation for *α* in the form

$$
\tanh\left(1.337\frac{a}{1+a}\right) = 0.8064(1101)^{\frac{3}{2(1+a)}}, \quad a \approx 140. \tag{36}
$$

Taking *ρmr* = 0.3*ρtot*, we obtain *ρDEV* = 0.0022*ρtot*, and the effective dark energy density at present time should be equal to Ω0*eff* = 0.7022, which is inside the limits of local measurement of Λ0, as mentioned above. The local matter density of the flat universe at present time is Ω0*<sup>m</sup>* = 0.2978, leading to the dark matter density Ω0*dm* = 0.2578, for the baryon density Ω0*<sup>b</sup>* = 0.04.

#### **5. Discussion**

We consider a cosmological model, where the DE input consists of two components: cosmological constant Λ, and small variable part DEV, which is connected with a scalar field, originated in the early universe and creating inflation. The energy density of this variable part is uniquely connected with a matter density, and the matter itself was created during the inflation stage. For simplicity we consider a cosmological model with a linear connection between energy densities as *αρDEV* = *ρDM*. This model could give a solution to the problem of Hubble Tension, which appears by using an inappropriate extrapolation when finding the present value of the Hubble constant from Planck observational data. We have solved the Friedmann equation in the presence of this relation, and have found the value of *α* at which HT disappeared. The present DEV density needed for explanation of HT phenomena is very small relative to the cosmological constant Λ. It influences the cosmological expansion at larger redshifts, where the input of the Einstein constant Λ is small. Presently the situation is opposite, Λ Λ*V*, because decreasing of matter density during cosmological expansion determines the transition from the quasi-Friedmann to quasi-de Sitter stage. The estimation of the density *ρ*Λ*<sup>V</sup>* at the present epoch corresponds to ΩΛ*<sup>V</sup>* ≈ 0.0022, which slightly increases the present dark energy density. In the flat universe it determines decreasing of the dark matter density due to transfer of its energy into DEV in the condition of "equipartition" at constant *α*.

We have used for estimations a constant ratio of *ρ*Λ*<sup>V</sup>* /*ρ<sup>m</sup>* = 1/*α* for the universe expansion after recombination, at *z* < 1100, but deviations from this law should not qualitatively change the conclusion that a relatively small average contribution of the variable ΩΛ*<sup>V</sup>* may explain the difference in Hubble constant measured at local and high-*z* distances.

The present parameters of the LCDM model have been estimated from the analysis of CMB fluctuation measurements in WMAP and PLANCK experiments, having a powerlaw spectrum of adiabatic scalar perturbations. The procedure is based on a search of extremes in the multidimensional parameter space. The presence of HT (if real) adds an additional restriction to this problem. The universe parameters obtained in this process may be changed with this additional restriction. The computations could be performed in the presence of a variable *α*. Decreasing of dark matter leads to decreasing of the field amplitude, which may prevent the energy exchange between DM and DEV in the absence of very light DM particles.

In our model the DM should be represented by wide mass spectrum particles, and not by unique mass CDM particles, which are usually considered now. By analogy with CMB, the lowest mass of DM particles should not presently exceed the value ∼ (Ω*DEV*/Ω*CMB*)1/4 × *kTCMB* ≈ <sup>7</sup> · <sup>10</sup>−<sup>4</sup> eV, to retain the possibility of an almost constant *<sup>α</sup>*.

**Funding:** This work was supported by the Russian Science Foundation (grant No. 18-12-00378).

**Acknowledgments:** The author is very grateful to O. Yu. Tsupko for valuable comments.

**Conflicts of Interest:** The author declares no conflict of interest.

#### **Note**

<sup>1</sup> In most equations below it is taken that *c* = 1.

#### **References**

