IJAA  Vol.11 No.2 , June 2021
Locally Varying Hubble Parameter in Terms of Reduced Friedmann Equation
Show more
Abstract: Assessment of the Hubble parameter as an indicator of the expansion rate of the universe holds a central position in the field of astronomy. From its initial estimate of about 500 km⋅sec-1⋅parsc-1, this value had been steadily amended as the observational tools became more accurate and precise. Despite this, a gap remains between the value of observations relating to local and nonlocal estimations of the Hubble parameter that gave rise to what became known as the Hubble tension. This tension is addressed here while dealing with space fabric as a cosmological fluid that undergoes transition.

1. Introduction

The ΛCDM model which represents standard cosmology had a remarkable success in matching the theoretical predictions of the 6-parameter model to the observational findings.

However, there had been many recent challenges to the integrity of this model in the form of local estimations of the Hubble parameter [1] [2] [3] [4] which have statistically significant difference from those values estimated by non-local measurements, mainly through CMB observations [5] [6].

In addition, the presence of enhanced lensing amplitude whose value exceeds the theoretical prediction ( A L > 1 ) [7] [8] is interpreted as an indication of a closed universe—in contradiction to the basic assumption of the ɅCDM model of a flat cosmology.

Here we will try to address those two questions concerning the Hubble tension and the curvature tension which were tackled previously [9] using degree of freedom approach, but this time the Friedmann equations would be employed.

2. Hubble Parameter-Radius of the Universe (rue) Relationship

The Hubble parameter is related to the universe radius via the scalar factor as follows

r u e = a ( t ) r u e o (1-2)

where r u e , r u e o denote the radius of the universe at any time (t) and at present time (to)

d r u e d t = r u e o d a ( t ) d t (2-2)

H = d a ( t ) d t a ( t ) = r u e o d a ( t ) d t r u e o a ( t ) = d r u e ( t ) d t r u e ( t ) (3-2)

3. Decelerating Pressure-Less Universe

This analysis is based on a model for quanton fields [10], which are space and time varying, when interacting they generate inflationary and gravitational interactions of space fabric which are at origin of dark energy and dark matter.

Here, the standard form of Friedmann equations will be tackled, while later on a slightly modified version will be introduced. Friedman acceleration equation is written as

r u e r u e = 4 π G 3 ( ϱ m + 3 p c 2 ) (1-3)

ϱ e : energy density = K ϱ t 2 (2-3)

r u e = K r u t 2 3 (3-3)

where mass density is defined as

ϱ m = ϱ e 2 c 2 (4-3)

K ϱ , K r u : constants of proportionality [9]. For the case of pressure-less fluid this reduces to

r u e r u e = 4 π G ϱ m 3 (5-3)

Given that

r u e = K r u t 2 3 (6-3)

r u e = d r u e d t = 2 3 K r u t 1 3 (7-3)


ϱ m = k ϱ c 2 t 2 (8-3)

Substituting with this equivalent value for r u e

r u e = 4 π G ϱ m 3 r u e = 4 π G k ϱ 3 c 2 t 2 r u e (9-3)

t = ( r u e k r u ) 3 2 , 1 t 2 = ( k r u r u e ) 3 (10-3)

substituting in (8-3) with value of 1 t 2 and integrating w.r.t r u e

r u e d r u e = 4 π G k ϱ 3 c 2 ( k r u r u e ) 3 r u e d r u e (11-3)

r u e d r u e = 4 π G k ϱ k r u 3 3 c 2 r u e 2 d r u e (12-3)

d ( r u e ) d t d r u e = d r u e d r u e d r u e d t d r u e = d r u e d t d r u e = r u e d r u e = ( r u e ) 2 2 (13-3)

( r u e ) 2 = 8 π G k ϱ k r u 3 3 c 2 r u e + constant (14-3)

( r u e ) 2 ( r u e ) 2 = 8 π G k ϱ k r u 3 3 c 2 r u e 3 + constant (15-3)

Given that k r u 3 r u e 3 = 1 t 2

( r u e ) 2 ( r u e ) 2 = H 2 = 8 π G ϱ m 3 + constant (16-3)

The constant of integration equals k c 2 r u e 2 , which is the usual form of Friedmann equation. If the curvature term is small enough to be neglected in comparison to density term.

( r u e r u e ) 2 = 8 π G ϱ m 3 c 2 (17-3)

which is the first Friedmann equation of Einstein DeSitter model for flat matter dominated universe this time in a closed cosmology model.

Now the Hubble parameter

H = d r u e ( t ) d t r u e ( t ) = 2 3 K r u t 1 3 K r u t 2 3 = 2 3 t (18-3)

Back to the first Friedmann equation and from before

8 π G ϱ m 3 = 4 9 t 2 = H 2 (19-3)

The closed inflationary model of the nature r u e = K r u = t 2 3 of a four dimensional energy density translates into 3 dimensional mass density of the type ϱ m = K a 3 . As for the three dimensional universe, the four degrees of freedom of expanding energy density are the result of adding three degrees of freedom for the qunanton volume in addition to one degree of freedom due to quanton splitting [9], and when considering length, its 4 3 degrees of freedom are the result of adding one degree of freedom of the quanton radius and ( 1 3 ) degrees of freedom as a result of quanton splitting.

In this closed inflationary model, the energy density plays the pivotal role due to its relationship with the energy degree of freedom as this was the major parameter which controlled the inflationary process and contributed to the uniformity and the homogeneity of space fabric as it is today.

Despite being a closed model, it had a clear matter density term dominance over curvature term especially in the earliest part of inflationary history.

The relative ratio between the density term ( 8 π G ϱ m 3 ) and the curvature term ( c 2 r u e 2 ) is as follows

R R = 8 π G ϱ m 3 c 2 r u e 2 = 8 π G K ϱ 3 c 2 t 2 K r u 2 t 4 3 c 2 = 8 π G K ϱ K r u 2 3 c 4 t 2 3 (20-3)

The inflationary process of the universe started with a dominant density term in comparison to the curvature term with a relative ratio RR of almost 1: 2 × 1022 (at t = t Planck ).

Gradually the curvature term assumed greater relative importance, and this gradual loss of dominance for the density term can be understood as energy density expansion is expressed as ϱ e = k ϱ t 2 in other words energy density inversely proportional to true time squared ( t 2 ), while the radius of curvature is put as r u e = k r u t 2 3 as the curvature term ( c 2 r u e 2 ) is inversely proportional to ( t 4 3 ).

The dominance of density term is manifested as a decelerating inflationary regime which weakens as time progresses, and gradually both terms converge. In fact for very small values of (t) the curvature term can be disregarded and this closed inflationary model can be viewed as very similar to flat model.

4. Transition Time and Afterwards

The inflation of the universe was accompanied by a reduction of the energy density which led to a weakness of decelerating effect with regard to the rate of change of its radius ( d r u e d t ). The deceleration came to a halt as the rate d r u e d t approached the value of (c). This happened as the space fabric whose constituent quantons must expand in a degree of freedom manner or in other words d r o d t = c , where r o denotes the outer radius of the universe.

From this point on, the universe radius is expressed in terms of the equivalent sphere whose radius equals r u e (which is also equal to Hubble radius), outer and inner radii r o , r i of a shell.

The transition from sphere to shell shaped universe is the result of a universe volume V u being smaller than volume encompassed by the outer shell V o , since the rate of generation of this volume d V u d t is smaller than the rate of increase of the outer sphere volume d V o d t [9].

The quantons now forth behave as a quantum relativistic (photon like) gas (keeping in mind that photons were defined as relativistic quantons) [10] which obey the relationship

3 p = U V = ϱ e [11] [12] (1-4)

The Friedmann acceleration equation

r u e r u e = 4 π 3 ( ϱ m + 3 p c 2 ) (2-4)

During post transition era, the pressure has an accelerating effect which nullifies the decelerating gravitational effect of density term.

In other words

3 p c 2 = ϱ m (3-4)

and consequently

r u e r u e = zero , or r u e = zero (4-4)

which leads to r u e = constant = c .

Under such conditions the second Hubble equation becomes

H 2 = c 2 r u e 2 (5-4)

r u e : the equivalent radius of the universe.

At present an equivalent universe radius equals 1.36 × 1026 meters and can be expressed as

r u e = c t H (6-4)

tH = Hubble time and this relationship is valid for t > tr. Figure 1 represents the

Figure 1. Representation of Hubble parameter transition from pressure-less fluid of the nature H 0 = 2 3 t to the relativistic fluid of the nature H 0 = c r u e .

transition from a decelerating fluid (the higher curve on the left) to a relativistic fluid (the lower curve on the left) where the transition time is the instance when the Hubble parameters of the two flow regimes were momentarily equal.

5. Local Hubble Parameter

Possible explanations for the inconsistency between the local and the universal values of the Hubble parameter, apart from systematic errors, range from early dark energy [13] [14], acoustic dark energy [15], modified gravitation [16], Dynamic dark energy [17] [18], and neutrino physics [19].

This model proposed an inflationary scenario that transits from sphere to shell shaped universe [9], under such a model the Hubble parameter can have a local value, while its universal value can still be assessed via non local phenomena like CMB.

Previously, the values of the Hubble parameter were obtained in terms of the outer ( r o ), inner radius ( r i ), and their rates of change with time, the same methodology can be elaborated further to allow for the determination of the Hubble parameter at any radius ( r x ).

Remembering that

r o ( t ) = c ( t t r ) + K r u t r 2 3 (1-5)

d r o d t = c (2-5)

t r : transition time

r i ( t ) = r o 3 3 4 π V u 3 (3-5)

d r i d r o = 1 r i 2 ( r o 2 1 4 π c d V u d t ) (4-5)

d r i d t = d r o d t d r i d r o = c d r i d r o (5-5)

Now for the Hubble parameter an any point on the radius r x .

Given that r o > r x > r i and

x = r x r i r o r i (6-5)

r x = x ( r o r i ) + r i = r o ( 1 x ) ( r o r i ) (7-5)

After rearrangement

r x = x r o + ( 1 x ) r i (8-5)

When differentiating w.r.t time

d r x d t = x d r o d t + ( 1 x ) d r i d t (9-5)

And the local Hubble parameter

H ( x ) = d r x d t r x = x d r o d t + ( 1 x ) d r i d t r x (10-5)

H ( x ) = x d r o d t + ( 1 x ) d r i d t x ( r o r i ) + r i = x d r o d t + ( 1 x ) d r i d t x r o + ( 1 x ) r i (11-5)

Figure 2 charts the variation of the Hubble parameter along the radial distance r x .

Figure 2. Representation for the variation of the local value of the Hubble parameter along the radial direction r x , r o > r x > r i .

6. Newtonian Interpretation of the Transition Time

Based on the Newtonian concept of gravitation [20] [21] an illustration can be made about the physical meaning of each of Friedmann equations, given that the mass of the universe M u , the second Newton law is written as

d 2 r u e d t 2 = G M u r u e 2 (1-6)

After multiplying both terms by d r u e d t and integrating, the energy equation is obtained.

1 2 ( d r u e d t ) 2 = G M u r u e + U (2-6)

where the density term represents the gravitational potential energy and the constant (U) represents the total energy of the universe, while the time derivative term represents the kinetic energy.

For the case of the symmetrically expanding universe the mass can be expressed as M u = 4 π ϱ e r u e 3 3 when dividing by r u e 2 .

In pre-transition era the kinetic energy equaled the gravitational decelerating potential, while in post transition universe total energy is represented by the kinetic energy as there was no decelerating potential.

The kinetic energy of the universe is transformed from being equal to the gravitational potential to being equal to the total energy in the universe in the post transition era.

7. The Dual Nature of Density Term

The motivation behind the following analysis is to study how the cosmological Transition took place given that it must have been a smooth one between the two epochs. The condition of such smoothness is that both the density and the curvature terms must momentarily be equal at the instant of transition. Previously, it was shown that the density term had a relative dominance in magnitude over the curvature term. The curvature term approaches the density term but not at any point the universe’s inflationary history did both terms become equal. Under such conditions the cosmological transition took place under a modified gravitation. Here a distinction is made between the energy density spatial expansion (whose pattern is uniform throughout the inflationary history of the universe), and the volumetric expansion which is altered at the transition instance. Previously the universe radius was defined as

r u e = k r u t 2 3 (1-7)

where the degrees of freedom corresponding to ( r u e ) are

D o f r u e = 4 3 (2-7)

This allows universe radius to be put in the form

r u e = k r u t D o f r u e 2 (3-7)

Similarly for the case of spatial expansion which corresponds to the two dimensional expansion of the quanton fields and the spatial radius is represented by r s , the degrees of freedom corresponding to the spatial radius

D o f r s = 4 2 = 2 (4-7)

r s = k r s t D o f r s 2 = k r s t (5-7)

The value of k r s equals the constant (c) as the spatial expansion corresponds a degree of freedom expansion of the quanton energy fields.

The spatial expansion is based on Hubble parameter of the discrete quanton Fields.

r s = c t (6-7)

v = d r s d t = c (7-7)

The spatial expansion based Hubble parameter (Hs) can then be found as follows

H s = d r s d t r s = c c t = 1 t (t: true time) (8-7)

Now the modified Friedmann equation becomes

H s 2 = 8 π G ϱ m 3 = 1 t 2 (9-7)

This equation remains time invariant. (pre- and post-transition), and this is understandable as the expansion nature of quanton fields remains constant irrespective of the inflationary nature of the universe.

The Friedman equation for the volumetric expansion in a pre-transition world now becomes

H v 2 = ( 2 3 t ) 2 = ( 2 3 ) 2 8 π G ϱ m 3 (10-7)

Hv: Hubble parameter for the volumetric expansion and the ratio between the volumetric and spatial expansion Hubble parameters for the pre-transition era

H v H s = t t H = 2 3 (11-7)

This ratio between the volumetric and the spatial expansion Hubble parameters (as well as the ratio between Hubble time and true time) was a constant one during the pre-transition era as it was a case of inflation into a spherically shaped universe.

In pre-transition time the density term played both roles, namely volumetric and spatial expansion, however this duality came to an end with the advent of the transition and afterwards, as the curvature term assumed the role of volumetric expansion while the density term maintained the role of spatial expansion and the ratio between the two in post transition era is no longer a constant one as it was

H v H s = c 3 4 π k v u t 2 3 1 t = t t H (12-7)

= c t 1 3 3 4 π k v u 3 (13-7)

where the universe volume is defined as

V u = k v u t 2 [9] (14-7)

and its equivalent radius

r u e = [ 3 4 π k v u 3 ] t 2 3 (15-7)

and the volumetric Hubble parameter of the post-transition world now becomes

H v = c r u e (16-7)

Previously the value of the vacuum energy density and the total energy content of the universe was provided based on empirical data [22] [23]. Now more refined value for both can be determined based on (9-7), which yields the following values of vacuum energy density and total energy content of the universe respectively.

ϱ e = 3 c 2 8 π G = 8.45 × 10 10 j m 3 (17-7)

E u = 8.9 × 10 69 joules (18-7)

Equation (9-7) allows for the following

3 c 2 = 8 π G k ϱ (19-7)

k ϱ = 3 c 2 8 π G (20-7)

For a closed universe spatially symmetric universe,

V u = E u ϱ e (21-7)

r u e = 3 E u 4 π ϱ e 3 = 3 E u 4 π k ϱ 3 t 2 3 = k r u t 2 3 (22-7)

Substituting with value of k ϱ from (20-7)

k r u = 3 8 π G E u 4 π 3 c 2 3 = 2 G E u c 2 3 (23-7)

r u e = 2 G E u c 2 3 t 2 3 (24-7)

V u = 8 π G E u 3 c 2 t 2 (25-7)

Equations (23, 24, 25-7) are very important as they relate the cosmological parameters in terms of the fundamental physical constants G, c in the absence of the related time constants (k’s).

The concept of uniform density expansion of the form 8 π G ϱ e 3 c 2 = 1 t 2 allows for the development of more refined values of actual energy density while maintaining the symmetries of the quanton physical and cosmological ratios as illustrated in the supportive data.

Figure 3 displays relationship between true time and Hubble time which is based on Equations (10-7) and (16-7) for pre and post transition respectively.

Figure 3. The extrapolated relationship between true time and Hubble time.

Note that the intersection point where the two values are equal marks another possible transition when density term and curvature term are equal in value.

Table 1 offers the modified version of the Friedmann equation for the spatial and volumetric expansion.

Table 1. The density term represents both the volumetric and the spatial expansion in the pre-transition Era while it represents only the spatial expansion in post transition era.

The value of Hubble time is greater than true time up till now, and consequently, the value of the term ( 1 t 2 ) is bound to be greater than the value of the term ( 1 t H 2 ), therefore, the actual mass density is always greater than the critical density.

Based on the previous conclusion the value of the curvature parameter can be obtained

Ω = ϱ m ϱ m c = t H 2 t 2 = 8.45 × 10 10 7.74 × 10 10 = 1.09 (26-7)

Ω k = 1 Ω = 0.09 (27-7)

An explanation of the near zero value of the curvature parameter as a result of the constraining of observational cosmological parameters [24] lies in the fact that this near zero curvature parameter (which is perceived as an evidence of a flat universe) belongs to a closed universe which is expanding at a rate that equals that of an open universe due to geometric degeneracy which allows various closed models to generate the same observational results as flat model, namely

ϱ m c ( perceived critical mass density ) = 3 H v 2 8 π G (28-7)

H v 2 = ( c r u e ) 2 (actual volumetric expansion) (29-7)

8 π G ϱ m c 3 = ( c r u e ) 2 (30-7)

ϱ m c = 8 π G c 2 3 r u e o 2 (31-7)

8. Effect of Presence of Normal Matter

8.1. Departure from Parameter Ratio Symmetry with Time

Previous work [9] the presented cosmological inflationary model which the quanton and the cosmological parameters were determined based on their energy degree of freedom. The relative ratios of those parameters between the present day and the Planck Era ( R R x ) were deduced in the form

R R x = ( t 0 t P ) D o f x 2 = ( 4.37 × 10 17 5.39 × 10 44 ) D o f x 2 = ( 2.84 ) D o f x 2 (1-8)

where R R x is the relative ratio between the value of the parameter (x) at the present day to its value at Planck time given its degrees of freedom (Dof). That model was based on either the absence of the normal matter or considering it to be as diffusive as dark matter. In either case the result was ideal ratios that do not deviate from that of time relative ratio. But in reality the agglomerated, non-diffusive nature of the normal matter reduces its gravitational effect so allowing energy density further expansion and as a result the quanton and the cosmological parameter relative ratios are higher by two or 3 hundredths from that ideal relative ratio.

8.2. Underestimation of the Energy Density

The extent of lensing of observational CMB radiation is affected by the matter distribution, as a result the observed curvature parameter ( Ω k = 0.045 ) is an underestimation of the actual space curvature while the ideal relationship 1 t 2 = 8 π G ϱ m 3 provides an overestimation of the energy density. As a result of the combination of those effects, the actual curvature parameter is somewhere between the observed value and the ideal value.

9. Ethical Statement

The author declares that this work fully complies with the ethical guidelines as had been stated by the journal.

10. Data Availability

The data that support the findings of this study are available from the corresponding author, upon request or alternatively it can be downloaded via the following links:




11. Discussion

Though this model adheres to the main outlines of the GR based standard cosmology but amendments were introduced which can be summarized as follows

1) Distinction is being made between spatial and volumetric expansion of energy density due to the nature of quantized two dimensional quanton space and time varying fields expanding in three dimensional space.

2) In the pre-transition era

a) The Hubble parameter is represented by the modified density term only despite being a closed model (k ≠ zero as density term is greater in magnitude than curvature term).

b) As a result of this distinction between spatial expansion volumetric expansion, a factor of ( 2 3 ) was introduced to density term to account for this difference when calculating the universal volumetric Hubble parameter.

3) In post transition era, introduction of a locally varying Hubble parameter which seems to be a departure from general principle of cosmology which necessitates that inflation is the same everywhere in the universe, but otherwise the reminder of that principle still holds since the degree of freedom expansion of the statistically distributed quanton densities as well as parameter synchronization ensured the space fabric uniformity and homogeneity.

Finally Equation (9-7) provided a constraint that broke the geometric degeneracy and provided a unique closed inflationary model that corresponds to observational values.

12. Conclusions

1) Cosmic transition from decelerating pressure-less fluid to a relativistic fluid allows Hubble parameter to be varied locally, while the universal Hubble parameter can be assessed in terms of curvature term of Friedman equation.

2) Both density and curvature terms carry their own time imprints (true time and Hubble time) which express spatial and volumetric expansion of energy density respectively.

3) Smooth cosmic transition requires that the density and the curvature terms be equal momentarily, this is possible only under modified gravitation.


A preprint of this paper can be found at

Cite this paper: Kassem, A. (2021) Locally Varying Hubble Parameter in Terms of Reduced Friedmann Equation. International Journal of Astronomy and Astrophysics, 11, 175-189. doi: 10.4236/ijaa.2021.112010.

[1]   Riess, A.G., Casertano, S., Yuan, W., Macri, L.M. and Scolnic, D. (2019) Large Magellanic Cloud Cepheid Standards Provide a 1% Foundation for the Determination of the Hubble Constant and Stronger Evidence for Physics beyond ΛCDM. Astrophysical Journal, 876, 85.

[2]   Wong, K.C., Suyu, S.H., Chen, G.F., Rusu, C.L., et al. (2019) H0LiCOW—XIII. A 2.4 per cent Measurement of H0 from Lensed Quasars: 5.3σ Tension between Early- and Late-Universe Probe. Monthly Notices of the Royal Astronomical Society, 498, 1420-1439.

[3]   De Jaeger, T., Stahl, B.E., Zheng, W., Filippenko, A.V., Riess, A.G. and Galbany, L. (2020) A Measurement of the Hubble Constant from Type II Supernovae. Monthly Notices of the Royal Astronomical Society, 496, 3402-3411.

[4]   Huang, C.D., Riess, A.G., Yuan, W., Macri, L.M., Zakamska, N.L., Casertano, S., et al. (2020) Hubble Space Telescope Observations of Mira Variables in the SN Ia Host NGC 1559: An Alternative Candle to Measure the Hubble Constant. The Astrophysical Journal, 889, 5.

[5]   Aghanim, N., Akrami, Y., Ashdown, M., et al. (2020) Planck 2018 Results VI. Cosmological Parameters. Astronomy & Astrophysics, 641, A6.

[6]   Aiola, S., Calabrese, E., Maurin, L., Naess, S., Schmitt, B.L., et al. (2020) The Atacama Cosmology Telescope: DR4 Maps and Cosmological Parameters. Journal of Cosmology and Astroparticle Physics, 12, 047.

[7]   Aghanim, N., Akrami, Y., Ashdown, M., et al. (2020) Planck 2018 Results. V. CMB Power Spectra and Likelihoods. Astronomy & Astrophysical Journal, 641, A5.

[8]   Ade, P.A.R., Aghanim, N., Arnaud, M., et al. (2016) Planck 2015 Results XIII. Cosmological Parameters. Astronomy &Astrophysics, 594, A13.

[9]   Kassem, A. (2020) Symmetries in the Universe, a Quanton Origin. International Astronomy and Astrophysics Research Journal, 2, 52-70.

[10]   Kassem, A. (n.d.) Quanton Based Model of Field Interactions.

[11]   Das, S. and Roychowdhury, D. (2020) Thermodynamics of a Photon Gas with an Invariant Energy Scale. Physical Review D, 81, Article ID: 085039.

[12]   Mansuripur, M. and Han, P. (2017) Thermodynamics of Radiation Pressure and Photon Momentum. Dholakia, K. and Spalding, G.C., Eds., Optical Trapping and Optical Micromanipulation XIV, Proceedings of SPIE Vol. 10347, 1-20.

[13]   Poulin, V., Smith, T.L., Karwal, T. and Kamionkowski, M. (2019) Early Dark Energy Can Resolve the Hubble Tension. Physical Review Letters, 122, 221301.

[14]   Niedermann, F. and Sloth, M.S. (2020) Resolving the Hubble Tension with New Early Dark Energy. Physical Review D, 102, Article ID: 063527.

[15]   Lin, M.X., Benevento, G., Hu, W. and Raveri, M. (2019) Acoustic Dark Energy: Potential Conversion of the Hubble. Physical Review D, 100, 063542.

[16]   Abadi, T. and Kovetz, E.D. (2021) Can Conformally Coupled Modified Gravity Solve the Hubble Tension? Physical Review D, 103, 023530.

[17]   Di Valentino, E., Ferreira, R.Z., Visinelli, L. and Danielsson, U. (2019) Late Time Transitions in the Quintessence Field and the H0 Tension. Physics of the Dark Universe, 26, 100385.

[18]   Pan, S., Yang, W., DiValentino, E., Saridakis, E.N. and Chakraborty, S. (2019) Interacting Scenarios with Dynamical Dark Energy: Observational Constraints and Alleviation of the H0 Tension. Physical Review D, 100, 103520.

[19]   Lancaster, L., Cyr-Racine, F.-Y., Knox, L. and Pan, Z. (2017) A Tale of Two Modes: Neutrino Free-Streaming in the Early Universe. Journal of Cosmology and Astroparticle Physics, 7, 033.

[20]   Tipler, F.J.Rigorous Newtonian cosmologyNewtonian cosmology revisited

[21]   Telkamp, H. (2016) Machian Derivation of the Friedmann Equation. Physical Review D, 94, 043520.

[22]   Oldershaw, R.L. (1989) Towards a Resolution of the Vacuum Energy Density Crisis.

[23]   Valev, D. (2014) Estimations of Total Mass and Energy of the Universe. Physics International, 5, 15-20.

[24]   Efstathiou, G. and Gratton, S. (2020) The Evidence for a Spatially Flat Universe. Monthly Notices of the Royal Astronomical Society: Letters, 496, L91-L95.