In astrophysics, the Hubble constant H₀ [1] is a parameter to analyze the universe. Nevertheless, it is also one of the lesser-known values.

In 1916, Einstein found the general relativity laws [2]. His equations expect that the universe is either expanding or in a Big Crunch. He could have been the first to predict the universe’s expansion, but influenced by the popular idea, Einstein forced his model to be static with a cosmological constant Λ. In 1922, Friedmann showed from relativity that the universe expands at a calculable rate [3]. In 1927, Georges Lemaitre published independent research [4], giving what is now known as Hubble’s law. In 1929, Hubble discovered the universe’s expansion [1]. Equation (1) gives Hubble’s law, with v being the receding speed in km·s⁻¹, D being the distance between the observed object and the observer, and H₀ being the Hubble constant. He measured about H₀ ≈ 500 km·s⁻¹·MParsec⁻¹. His high value was due to a wrong calibration of the cepheids used to evaluate distances. Hubble’s law was correct, but H₀ was remaining to be found with accuracy.

v = D H 0 (1)

Physicists get H₀ based on far cosmic objects (Cepheids, supernovae, red giants, etc.) or local measurements (CMB, universal gravitational constant G, etc.). Including error margins of published values (see the software in Annex A), H₀ is between 19 to 174 km·s⁻¹·MParsec⁻¹. However, two values are often measured ~69.2 and ~72.1 km·s⁻¹·MParsec⁻¹. An irreconcilable tension between some H₀ values shows up [5]. Even with good accuracies, their error margins do not always overlap. It may let us think that only one of these values is right. No one considered it possible that both values may be in some way correct.

In 2019, we wrote an article [6] explaining what may look like two values for H₀. We calculated the universe age, obtained a result of complex type, and an apparent age of the universe of ~14.4 billion years. The complete equation may be approximated by 1/H₀, giving ~13.56 billion years. We notice that there is a difference of ~4.25% between the approximated and the non-approximated values.

Cosmologists use 1/H₀ to calculate the universe’s age. Thus, if we could measure the apparent age of the universe with no approximation, we would conclude wrongly that the Hubble constant is ~4.25% lower than it should be.

We hypothesize that two values of H₀ are somehow obtained from an approximated and non-approximated equation of the apparent age of the universe. The confusion leads to a tension between two values when there should be only one.

We summarize our cosmological model [6] to get H₀ as a function of α, c, and r_e. We found ways to measure H₀ locally by using the Cosmological Microwave Background (CMB) temperature T and by using the universal gravitational constant G [6] [7]. Based on our model, we found a theoretical equation to calculate H₀ from CODATA values (Committee of Data for Science and Technology) [8].

H 0 = c α 19 β r e ≈ 72.09548580 ( 32 )     km ⋅ s − 1 ⋅ MParsec − 1 where β = 3 − 5 ≈ 0.76 (2)

We want to validate this theoretical value of H₀ and highlight the tension between two measured values of H₀. We list the results of the most recent measures of H₀ and build a graph showing someway the popularity of each H₀ value range.

We list 42 H₀ equations. Certain overcome the difficulties to do experimental measurements. We use one of them as a third measurement of H₀. Our cosmological model shows that H₀ and the speed of light are not constant.

A compact form of notation is used to display tolerances (i.e., 2.734(10) K means 2.734 ± 0.010 K). The CODATA 2014 [8] is used to compare the results of our new equations with the articles published in 2019 and 2020.

Light speed in a vacuum c = 2 99792458   m ⋅ s − 1

Permeability of free space μ 0 = 4 π × 10 − 7     NA − 2

Permittivity of free space ε 0 ≈ 8.854187817 × 10 − 12     F ⋅ m − 1

Universal gravitational constant G ≈ 6.67408 ( 31 ) × 10 − 11     m 3 ⋅ kg − 1 ⋅ s − 2

Electron rest mass m e ≈ 9.10938356 ( 11 ) × 10 − 31   kg

Classical electron radius r e ≈ 2.8179403227 ( 19 ) × 10 − 15   m

Electron charge q e ≈ − 1.6021766208 ( 98 ) × 10 − 19     C

Planck length L p ≈ 1.616229 ( 38 ) × 10 − 35   m

Planck time t p ≈ 5.39116 ( 13 ) × 10 − 44   s

Planck mass m p ≈ 2 .176470 ( 51 ) × 10 − 8     kg

Planck constant h ≈ 6.626070040 ( 81 ) × 10 − 34     J ⋅ s

Fine-structure constant α ≈ 7.2973525664 ( 17 ) × 10 − 3

Boltzmann constant k b ≈ 1.38064852 ( 79 ) × 10 − 23   J ⋅ K − 1

Rydberg constant R ∞ ≈ 10973731.568508 ( 65 )   m − 1

Our theory is based on a cosmological model officially shown in 2019 [6], but it summarizes papers we wrote at www.pragtec.com/physique since 2011. First, we outline some main milestones as we did in 2020 [7].

We calculate the universe’s age with our cosmological model to understand what seems to be two potential values of H₀. We measure H₀ by observing cosmological objects. Universe’s age Δt_u is of complex type and results from the integral of the inverse of the expanding speed of the material universe v_m with the element of distance dr evaluated between the universe’s center of mass (at r = 0) and the apparent material universe radius of curvature at our location r_u.

Δ t u = ∫ r = 0 r = r u 1 v m ( r ) d r = Δ t h u + Δ t 0 h (35)

The Δt_hu is the time elapsed between the horizon (r = r_h) and here (r = r_u):

Δ t h u = ∫ r = r h r = r u 1 v m ( r ) d r ( Realtyperesult ) (36)

The Δt_0h is the elapsed time between r = 0 and the horizon r = r_h:

Δ t 0 h = ∫ r = 0 r = r h 1 v m ( r ) d r ( Imaginarytyperesult ) (37)

At the universe horizon r = r_h, the speed of light is zero. We cannot see beyond the horizon. The delay Δt_hu is the time elapsed between the horizon h and our actual position r_u in the universe. The delay Δt_0h is the time elapsed between the center of mass of the universe and the horizon r_h (given by Equation (38)).

r h = 2 G m u k 2 (38)

Performing the integral calculation of Equation (35), we get Equation (39).

∫ 1 v m ( r ) d r = ( z ( r ) + 2 G ⋅ m u ln ( 2 [ k 2 r + z ( r ) ] ) ) β k 3 where     z ( r ) = k 4 r 2 − 4 G 2 m u 2 (39)

We can decompose the age of the universe Δt_u into two parts, Δt_hu and Δt_0h. The value Δt_hu represents the time elapsed between r = r_h (at the horizon) and our actual position r = r_u in the universe. The value Δt_0h gives the time elapsed between r = 0 (at the Big Bang) and r = r_h (at the horizon).

Δ t u = Δ t h u + Δ t 0 h ≈ ( 9.50 + 10.47 i ) × 10 9     yearswhere     i = − 1 (40)

The imaginary time Δt_0h means that it elapses independently of our time. We cannot see an event between r = 0 and r = r_h, and an observer located between r = 0 and r_h could not see us. The Δt_hu equation is:

Δ t h u = 1 H 0 ( ( ω + 2 ln [ ω + β ( 2 + 5 ) ] − ln ( 4 ) ) 2 22 + 10 5 ) ≈ 7 10 H 0 where ω = β 2 ( 9 + 4 5 ) − 4 (41)

The precise equation for Δt_0h is:

Δ t 0 h = − ( 2 + π ) H 0 22 + 10 5 i ≈ − 77 100 H 0 i (42)

The modulus of the complex age Δt_u gives the universe’s apparent age T_u.

T u = | Δ t u | = | Δ h u + Δ t 0 h | = ( Δ t h u ) 2 + ( Δ t 0 h i ) 2 ≈ 14.14 × 10 9     years (43)

T u = 1 H 0 ( ( ω + 2 ln [ ω + β ( 2 + 5 ) ] − ln ( 4 ) ) 2 22 + 10 5 ) 2 + ( [ − ( 2 + π ) 22 + 10 5 ] i ) 2 ︸ η ≈ 1 ( withabout4 .25%oferror ) (44)

As the square root over the accolade is approximatively equal to 1, we get:

T u ≈ 1 H 0 ≈ 13.56 × 10 9     years (45)

The value of the correction factor between Equations (43) and (45) is η.

η = H 0 T u = H 0 ( Δ t h u ) 2 + ( Δ t 0 h i ) 2 ≈ 1.042516951 (46)

This η explains why scientists currently measure two values of H₀. Scientists can only size the apparent age of the universe with different techniques. They cannot measure the real part and the imaginary part of the universe’s age.

There is no “local” or “far” value of H₀. There is only one H₀. Some techniques give H₀ directly, and others need a correction factor. There is no need for any correction factor when H₀ is calculated from Equation (33), measured with the CMB temperature with Equation (16), or with the universal gravitational constant G with Equation (34). Other techniques may get similar results than Equation (43), and if we impose that value to fit with Equation (45), we get H ′ 0 .

H ′ 0 = 1 | Δ t u | ≈ 1 14.14 × 10 9     years ≈ H 0 η ≈ 69.2 km ⋅ s − 1 ⋅ MParsec − 1 (47)

However, Equation (45) gives the actual H₀ value:

H 0 ≈ 1 13.56 × 10 9     years ≈ 72.1 km ⋅ s − 1 ⋅ MParsec − 1 (48)

If scientists could measure the real part of the universe’s age and associate this value with 1/H₀, they would obtain the following value.

H 0 ≈ 1 9.50 × 10 9     years ≈ 102.94 km ⋅ s − 1 ⋅ MParsec − 1 (49)

If scientists could measure the imaginary part somehow, the association of this value with 1/H₀ (like in Equation (45)) would give the following H₀ value.

H 0 ≈ 1 10.47 × 10 9     years ≈ 93.39 km ⋅ s − 1 ⋅ MParsec − 1 (50)

with different types of experiments to measure the apparent age of the universe, scientists usually get either ~H₀ ≈ 69.2 or ~72.1 km·s⁻¹·Mparsec⁻¹. We assume that all calibration factors are used. New techniques could require other unknown corrective factors that have nothing to do with the related phenomenon.

The articles rarely give enough details to check if the process used needs η. Scientists must verify if the η factor is required for their approach.

In 1929, Hubble made the first observational-based measurements with cepheids and got H₀ ≈ 500 km·s⁻¹·MParsec⁻¹ [1]. Sadly, even with a correct principle, his value is higher than the typical value due to errors in distance calibrations.

Let us validate our theoretical H₀ with an adequate interpretation of 508 measurements found on the Internet. The ends of their tolerance ranges give 1016 values. To find H₀ that has the highest probability to be measured, we compile the number of crossings with the tolerance ranges for each value of H₀. It generates a curve with two tips (Figure 1). The higher it is, the greater the chances are that this value of H₀ may be part of many tolerance ranges among the collected data.

A simple statistical phenomenon may be described with a Gaussian function. For fitting a wavy curve, it is necessary to make the sum of many Gaussians. A simpler model with fewer degrees of liberty must always be privileged.

A curve fit is done by summing different Gaussians (shown in Figure 2). A better gap fitting reduces the risk of finding other results. Thus, we gave a heavier weight (×10) to all data located between 69.2 and 72.1 km·s⁻¹·MParsec⁻¹ (from our theory). We tried with and without this approach, and it gives about the same result. As it improves the gap fitting, we kept this approach.

Each tip in Figure 1 is approximated in Figure 2 with two positive Gaussians. We force these curves to be around two means, even though there are four positive curves. It removes two degrees of liberty. We must add a negative Gaussian to model the gap between the two mean values. We must elaborate on this negative Gaussian. Our theory predicts “two close values” of H₀. On the curve, a deep gap shows up. It is impossible to get such a gap by only adding positives Gaussians which give two little bumps without any gap. To get a real gap, we must add a negative Gaussian. Let us see in Figure 3 what would look like a curve fit without any negative Gaussian. Since the tips are close, they mix up to build only one tip.

The Gaussian sum in Figure 3 peaks around H₀ ≈ 71.11 km·s⁻¹·MParsec⁻¹. The result is not close to our theoretical H₀ ≈ 72.09548580 km·s⁻¹·MParsec⁻¹ (Equation (33)), but it is about what is found if statistics were used through the whole data set, thinking they should see only one tip. Moreover, Jang & Lee showed a similar value of H₀ ≈ 71.17 km·s⁻¹·MParsec⁻¹ (listed in our software in Annex A) that supposedly reduces the tension between the values obtained by cepheids (calibrated on SNe Ia) and CMB.

In Figure 1, we find two groups around H₀ ≈ 69.7 and 71.8 km·s⁻¹·MParsec⁻¹. It is known that there is currently a tension between two groups [5]. A significant gap appears between the two tips. The only way to create such a gap is to

withdraw values nearby a specific value. It would then create a negative Gaussian, such as in Figure 2. It is delicate to debate why some values may have been withdrawn. It could be intentional or not. In the past, it was difficult to see a difference between these groups. Now, the tolerances are small enough to clearly see two groups. With recent growing tensions between these two clans, some may be inclined to shrink or shift some tolerance ranges when it overlaps with neighbor values.

In Figure 4, we apply η to the curves around H₀ ≈ 69.882 km·s⁻¹·MParsec⁻¹. Then, all curves stand around H₀ ≈ 72.36 km·s⁻¹·MParsec⁻¹. Then, with the curves of Figure 4, we build the curve in Figure 5. Figure 6 is a zoom of its tip.

We want to know the precise value of H₀ for which the derivative of the Gaussian summation is 0. It corresponds to the highest probability of getting the true H₀ value. Unfortunately, the derivative of a Gaussian summation is not an easy equation to get in a software. We rather use a numerical technic to get it. In Figure 6, we show a zoom of the quadratic curve fit around the tip value. Using the equation, we take the derivative and find its maximum. The quadratic equation has the following form:

y ( x ) = A x 2 + B x + C Atthetip,theslopeis : y ′ ( x ) = d x d y = 2 A x + B = 0 → x = − B 2 A where x = H 0 , y = numberofcrossingswithdifferenttoleranceranges (51)

H 0 = − B 2 A ≈ 2423.2459592464 2 × 16.8057572117 ≈ 72.0957088907   km ⋅ s − 1 ⋅ MParsec − 1 (52)

This result is well centered on our theoretical value within 3 parts per million. Our approach considers that both clans are someway right. Indeed, their different approaches and results also highlight a new phenomenon. It gives credit to our theory of the universe’s complex age that predicts a few possible fake H₀ values.

We have 508 data. Each has a tolerance range (that may be symmetrical or not) that generates two H₀ values. Therefore, there are a total of i_max = 1016 data at the end. The following equation depicts the statistical error e_t:

e t ≈ ∑ i = 1 i = i max e i 2 i max (53)

We mention that 16 H₀ values in our software in Annex A come from statistics. We kept them since some are mixed up with new valuable data information. So, we modify Equation (53) to remove them to reduce their impact on the total e_t error. We use the following equation where n = 2 × 16 = 32 (each data generates two H₀ values) is the number of elements to exclude from our sample. The total e_t error reduces with the square root of the number of elements included in our sample.

e t ≈ ∑ i = 1 i = i max e i 2 i max ⋅ [ i max i max − n ] (54)

If n = 0, we fall back on Equation (53). With i_max = 1016, Equation (53) gives e_t ≈ ±0.32 km·s⁻¹·MParsec⁻¹. With i_max = 1016 and n = 32, Equation (54) rounds up to e_t ± 0.33 km·s⁻¹·MParsec⁻¹. The impact of these n elements has a very little impact.

H 0 ≈ 72.0957 ± 0.33   km ⋅ s − 1 ⋅ MParsec − 1 (55)

In Annex A, we supply the software used to get this result. All the main steps enumerated in this article are clearly shown. The software uses starting values (found via Excel) to fit the original curve with 5 Gaussian curves (#0 to #4 to use the same numbers as the software). Each Gaussian uses three parameters: μ is the mean value, σ represents the variance, and m is a multiplication factor.

f ( H 0 ) ≈ m σ 2 π e − 1 2 ( H 0 − μ σ ) 2 (56)

Here are the values for the 5 Gaussian curves used to fit the original curve:

For Gaussians #1 and #2, we force the software to use the same mean value. We do the same thing for Gaussians #3 and #4. We also note that the multiplication factor m of Gaussian #0 is negative. With these values, we stopped iterating when the sum of squares of errors was lower than 22000. We see in Figure 1 that the obtained approximated curve fits well the original curve. In our software (Annex A), the iterations start with values close to what they should be.

The specificity of our approach is to say that the two clans are someway right. However, we must apply a correction factor to one of them. Indirectly, it gives credit to a complex universe age that predicts a few possible fake values of H₀.

After reading this article, scientists should continue their work as they were doing, without applying any correction factor to their raw data. The correction factor should only be used on the final Gaussian curve to analyze data.

To avoid repeating everything unnecessarily, we recall different identities that will be used later to determine H₀. Planck units are commonly defined as follows.

Planckmass : m p = h c 2 π G ≈ 2.18 × 10 − 8     kg (58)

Plancktime : t p = h G 2 π c 5 = L p c ≈ 5.91 × 10 − 44     s (59)

Plancklength : L p = h G 2 π c 3 = c t p ≈ 1.61 × 10 − 35     m (60)

PlanckTemperature : T p = h c 3 2 π G k b 2 = m p c 2 k b ≈ 1.42 × 10 32     K (61)

Planckcharge : q p = 2 c h ε 0 = − q e α ≈ 1.88 × 10 − 18     C (62)

The fine-structure constant α is linked to Rydberg constant R_∞ and the electron mass m_e by the following equation:

R ∞ = c m e α 2 2 h (63)

The speed of light c is given as a function of μ₀ and ε₀.

c = 1 μ 0 ε 0 (64)

Associating the mass-energy of a Planck particle with its wave energy and then, using Equations (31), (32), (64), and (62), we get Planck charge q_p defined several ways and as a function of c, G, and h like the other Planck units.

q p = 2 h c μ 0 = 4 π m p L p μ 0 = β G h 2 π μ 0 r e 2 c 4 α 19 (65)

The electron’s charge is determined from the mass of the electron m_e, the classical electron radius r_e, and the vacuum permeability μ₀.

q e = 4 π m e r e μ 0 ≈ − 1.60 × 10 − 19     C (66)

Let us calculate the precise value of the average temperature T of the CMB. We first make equal Equations (16) and (33). Then, we replace G by Equation (31), and we get rid of Planck constant h by its value from Equation (32).

T = m e c 2 k b ( 15 β 6 α 17 π 3 ) 1 / 4 ≈ 2.7367958 ( 16 )   K (67)

This CMB temperature is like Kimura with 2.737 K [23].

For an academic purpose and to show the interdependence of H₀ with the other “constants”, we will enumerate equations using various universe parameters. Some overcome the inherent difficulties in measuring H₀ and show a roundabout way of obtaining an accurate value of it. We also find some others which depend on interesting values, or more precise ones. Using the constants c, k_b, T, m_e, r_e, h, G, μ₀, ε₀, m_u, R_u, R_∞, q_e, q_p, t_p, l_p, T_p, m_p, m_ph, and β , we find many equations.

The H₀ parameter is not constant since 1/H₀ represents an approximation of the apparent universe’s age, and H₀ get smaller over time. Since the universe is old, H₀ changes slowly. If the constancy of all the universe’s parameters is maintained as it is currently done in metrology, the universe’s age and H₀ will seem constant.

Results of 508 different experiences reduce the error by 508^1/2 ≈ 22.5. It may look like a significant number, but it is nothing besides what has been done to measure the electron characteristics accurately. Particle accelerators use millions of electrons at each experiment, and they repeat these many times to find something new. Computers analyze the collisions’ results to make the electron’s characteristics more and more accurate. It is why there is no manner to get better results than that of Equation (33), as it is based on well-known characteristics of the electron. We will see further many other equations that give precise results.

Replacing G by Equation (31) in Equation (16), we get Equation (68).

H 0 = π 3 k b 2 T 2 α 9 8 r e 15 m e c 3 h 3 β 5 (68)

Replacing h in Equation (68) by using Equation (32), we get Equation (69).

H 0 = k b 2 T 2 π 3 α 21 15 r e 2 m e 4 β 5 c 6 (69)

Replacing h in Equation (16) by using Equation (32), we get Equation (70).

H 0 = k b 2 T 2 β 2 π 3 G α 15 r e 3 m e 3 c 8 (70)

Replacing T in Equation (70) by using Equation (67), we get Equation (71).

H 0 = β α 9 G m e r e 3 (71)

Using Equation (32) in Equation (33), we get Equation (72).

H 0 = 2 π m e c 2 α 18 β h (72)

With Equations (63) and (72), we get the most accurate equation.

H 0 = 4 π c R ∞ α 16 β (73)

Using Equation (31) in Equation (73), we get Equation (74).

H 0 = 4 π G m e R ∞ β 3 / 2 c r e α 4 (74)

Equation (74) gives H₀ with G₂₀₁₄ (from CODATA 2014). The measurement of H₀(G₂₀₁₄) ≈ 72.105(36) km·s⁻¹·MParsec⁻¹ is our fourth way to measure H₀.

Using Equation (63) in Equation (74), we get Equation (75). This equation is also a good candidate for measuring H₀ as a function of G.

H 0 = 8 π G h R ∞ 2 β 3 / 2 r e c 2 α 4 (75)

Using Equation (32) in Equation (75), we get Equation (76).

H 0 = 16 π 2 G m e R ∞ 2 β 3 / 2 c α 7 (76)

This equation is another good candidate for measuring H₀ as a function of G. The measure still gives the same result as Equation (74).

We will enumerate other equations without making all the rather fastidious demonstrations. However, all these may be found from previous equations.

H 0 = 2 π m p h c 2 h (77)

H 0 = 2 c 2 π R ∞ α 55 G m e (78)

H 0 = β α 9 m e c r e 3 (79)

H 0 = 8 β G m e ( π R ∞ ) 3 α 9 (80)

H 0 = 8 c β r e ( π R ∞ ) 3 α 29 (81)

H 0 = 2 π c 5 α 57 h G (82)

H 0 = 2 c π β R ∞ α 25 r e (83)

H 0 = 2 π β m e c 3 α 37 h r e (84)

H 0 = 2 π 2 h ( k b 4 T 4 15 m u α 2 c 2 β 4 ) 1 / 3 (85)

H 0 = 2 π k b T p h α 57 (86)

H 0 = c μ 0 q e 2 α 19 β 4 π m e r e 2 (87)

H 0 = q e 2 α 19 β 4 π c ε 0 m e r e 2 (88)

H 0 = q p 2 α 18 β 2 h ε 0 r e (89)

H 0 = c μ 0 q p 2 α 20 β 4 π m e r e 2 (90)

H 0 = q p 2 α 20 β 4 π c ε 0 m e r e 2 (91)

H 0 = c β q p 2 α 19 4 π ε 0 k b T p r e 2 (92)

H 0 = c β q e 2 α 17 4 π ε 0 k b T p r e 2 (93)

H 0 = c q e 2 β 2 4 ε 0 k b T r e 2 ( 15 α 97 π 7 ) 1 / 4 (94)

H 0 = 2 k b T β h ( π 7 α 55 15 ) 1 / 4 (95)

H 0 = k b T β r e m e c ( π 3 α 59 15 ) 1 / 4 (96)

H 0 = 4 π k b T β c μ 0 q e 2 ( π 3 α 59 15 ) 1 / 4 (97)

H 0 = 4 π k b T β c μ 0 q p 2 ( π 3 α 55 15 ) 1 / 4 (98)

H 0 = G k b T c 3 r e 2 ( π 3 15 α 21 ) 1 / 4 (99)

H 0 = 2 G m e k b T h r e c 2 ( π 7 15 α 25 ) 1 / 4 (100)

H 0 = G m p k b T m e c 3 r e 2 ( π 3 α 21 15 β 2 ) 1 / 4 (101)

H 0 = G m u k b T m e c 3 r e 2 ( π 3 α 135 15 β 2 ) 1 / 4 (102)

H 0 = G m e k b T m p c 3 r e 2 ( π 3 β 2 15 α 63 ) 1 / 4 (103)

H 0 = 16 c G k b T ε 0 m e 2 q e 4 ( π 11 15 α 21 ) 1 / 4 (104)

H 0 = β m e 2 8 π 3 G c 3 α 15 h 3 (105)

H 0 = c 3 α 39 G m e β (106)

The last equation measures H₀ from G since all other constants are accurate. Many equations are excellent candidates for measuring H₀ as a function of G or T. These equations could represent valuable tools for cosmologists.

This document gives 42 equations of H₀ as a function of various universe parameters. Since H₀ may be defined using different parameters, we suggest that some of the most critical universe parameters are well linked, as much in the infinitely small as in the infinitely large, and H₀ is part of these.

We want to explain why Hubble parameter H₀ cannot be constant over time. As simple as it is, the reverse of Hubble parameter H₀ is related to the apparent age of the universe (see Equation (45)). Consequently, the H₀ parameter is changing over time. It is, therefore, by abuse of language that we call H₀ the Hubble “constant”. To be more precise, we should say the Hubble “parameter”.

When H₀ is expressed in km·s⁻¹·MParsec⁻¹, the ninth digit after the dot changes every year. It goes completely unnoticed. More than that, even if we could achieve this precision in our measurements of H₀, it would still go unnoticed since we forced c to be constant in 1983. In metrology, scientists choose the speed of light as a standard. Even though c changes every year, if we force it to be constant, we willfully readjust all other constants and units (distance, time, and mass) as a function of c to keep it constant. Then, H₀ looks constant as other parameters.

This article aimed to show that our theoretical value from Equation (33) (giving H₀ ≈ 72.09548580(32) km·s⁻¹·MParsec⁻¹) [6] is the right one, despite a growing tension [5] between values around 69.2 and 72.1 km·s⁻¹·MParsec⁻¹.

With 508 data (from [24] to [310] shown in our software in Annex A), a graph showing the actual tension [5] between two values is shown. We decomposed the curve into Gaussians. A negative one is required to explain the large gap between the two H₀ values, and it is due to withdrawn values. So, we restored them by removing that curve. Then, we applied a η ≈ 1.042516951 correction factor (from our theory) to the curves located at ~H₀ ≈ 69.2 km·s⁻¹·MParsec⁻¹. Our theory highlights a misunderstanding of the link between 1/H₀ and the universe’s apparent age. With the proper correction factor applied, we get a statistical value of H₀ ≈ 72.0957 ± 0.33 km·s⁻¹·MParsec⁻¹, which is close to our theoretical value. Our discovery of the η factor may help to reduce the tension between scientists. Someway we show that even if two H₀ values seem to be commonly found with various techniques, both are accurate if a proper correction factor is used.

With a new cosmological model, we get an apparent age of the universe of about 14.14 billion years. The exact formula is approximated from an elaborate integral result by the well-known 1/H₀ equation that gives 13.56 billion years. Different techniques may lead to either value. It depends if it is an attempt to measure the universe’s age locally or far away. There is no “local” or “distant” value of H₀, as some may pretend [46] [47]. Sticking their measurement of the apparent age of the universe to 1/H₀, most cosmologists get results that stand around 69.2 or 72.1 km·s⁻¹·MParsec⁻¹. Our hypothesis may explain the actual tension [5] relative to these two values. However, there is only one true H₀ value, and the other one is just misinterpreted as being the Hubble constant without quite being so.

Even if many theoretical equations of H₀ are shown in this article, we highlight that we also found a few interesting ways to measure the H₀ accurately using the CMB temperature T and the value of the universal gravitational constant G from CODATA 2014. These results confirm our theoretical value.

H₀ ≈ 72.06(90) km·s⁻¹·MParsec⁻¹ and 71.95(50) km·s⁻¹·MParsec⁻¹ from T,

H₀ ≈ 72.086(36) km·s⁻¹·MParsec⁻¹ and 72.105(36) km·s⁻¹·MParsec⁻¹ from G.

For an academic purpose, we enumerated 42 equations of H₀ using different parameters. These equations showed that H₀ is intricated with all other “constants”. For metrology purposes, the speed of light in a vacuum is forced to be constant to be an unchanging standard. If this situation is considered valid in a metrology context, H₀ should also be considered constant and become part of the CODATA. However, if 1/H₀ represents an approximation of the universe’s age, it would also make sense to say that H₀ is changing over time.

Einstein’s and Schwarzschild’s equations show that massive objects such as the universe influence the speed of light. As the universe expands, its density diminishes, and the local speed of light increases over time.

The fine-structure constant α is unitless and may be described as a ratio where the variation rate at the numerator counterbalances the variation rate at the denominator. Apart from α and β , all “constants” used to describe H₀ in our equations somehow emanate from fundamental units such as the meter, the second, and the kilogram. These units are now defined by the speed of light. As H₀ describes the universe’s age and depends on many unit-dependent “constants” based on c, we should consider c and all universe’s unit-dependent parameters as changing over time. Forcing c to be constant is necessary for metrology purposes, but it is not in the interest of physicists for explaining phenomena. An accurate value of H₀ has a great interest in deepening our understanding of the universe.

The author claims that he has no conflict of interest in connection with the publication of this article.

Mercier, C. (2021) A New Physics Would Explain What Looks Like an Irreconcilable Tension between the Values of Hubble Constants and Allows H₀ to Be Calculated Theoretically Several Ways. Journal of Modern Physics, 12, 1656-1707. https://doi.org/10.4236/jmp.2021.1212098

// This software finds "the best" experimental value of H0 with a set of 508 data

//Compiled on Dev-C++ 5.11 available for free at:

// https://sourceforge.net/projects/orwelldevcpp/

#include

#include

#include

#define printf __mingw_printf

#define nbH0 508 //Number of measurements of H0 analyzed

#define Pi 3.141592654 //Definition of Pi

double Mean[5],Sigma[5],Multiplier[5];//Caracteristics of Gaussian curves

double A,B,C; //Coefficients of the quadratic equation of the tip

double LMSTip; //Least mean square for the approximation tip curve

double BEH; //Best estimate of H0

double H0[2*nbH0-1],NbCrossings[2*nbH0-1];

double nbCrossings[2*nbH0-1]; //Vector of number of crossings

double Accuracy_ppm; //Accuracy of H0 compared to the theoretical value

double TVH = 72.09548580; //Theoretical value of H0 (in km/(s*MParsec)

unsigned int PosTipIndex; //Index corresponding to the tip of nbCrossings array

int n = 20; //Sample before & after PosTipIndex to build the tip equation

double Sqr(double value) {return value*value;} //***Returns the square value

//***Function that returns the square root value

double Sqrt(double Value) {

double D; /*Dummy value*/ double V; /*Returned value*/ int i; /*Counter*/

V=0; D=Value;

for (i=0;i<=50;i++) { //Gives 50 bits of precision

D=-D/2;

if (D<0) while (V*V>Value) V=V+D; else while (V*V

return V;

} //End of Sqrt

double Exp(double Value) { return pow(2.718281828,Value);} //***e^Value

//***This procedure creates a table of 2*nbH increasing values

//***of H0[i] with the tolerances T[i] and the signs +/- Variation[i].

void CreateTableOfCrossingH0Ranges(void) {

double PT[nbH0],NT[nbH0]; //Positive tolerance and negative tolerance

double Variation[2*nbH0]; //Variation from the tolerance range

double DH0,DSgn; //Dummy H0 and Sgn used to put H in increasing order

int DeltaCrossings; //Variation of on the number of crossings

int i,j,k; //Counters

i=0;

//We enumerate all H0 values and their tolerance range found on Internet

//In brackets, we add the bibliographic references

H0[i]=69; PT[i]=+16; NT[i++]=-8; // [24] Abbott et al.

H0[i]=70; PT[i]=+2.7; NT[i++]=-2.7; // [25] Addisson

H0[i]=72.4; PT[i]=+3.9; NT[i++]=-4.8; // [25] Addisson

H0[i]=73.1; PT[i]=+3.3; NT[i++]=-3.9; // [25] Addisson

H0[i]=73.2; PT[i]=+1.3; NT[i++]=-1.3; // [25] Addisson

H0[i]=68.7; PT[i]=+1.3; NT[i++]=-1.3; // [25] Addisson

H0[i]=73.5; PT[i]=+5.3; NT[i++]=-5.3; // [26] Baxter & Sherwin

H0[i]=73.3; PT[i]=+0.7; NT[i++]=-0.7; // [27] Blakeslee et al.

H0[i]=73.78; PT[i]=+0.84; NT[i++]=-0.84; // [28] Bonilla

H0[i]=73.577; PT[i]=+0.106; NT[i++]=-0.106; // [29] Dainotti et al.

H0[i]=73.493; PT[i]=+0.144; NT[i++]=-0.144; // [29] Dainotti et al.

H0[i]=73.222; PT[i]=+0.262; NT[i++]=-0.262; // [29] Dainotti et al.

H0[i]=73.664; PT[i]=+0.223; NT[i++]=-0.223; // [29] Dainotti et al.

H0[i]=73.576; PT[i]=+0.105; NT[i++]=-0.105; // [29] Dainotti et al.

H0[i]=73.513; PT[i]=+0.142; NT[i++]=-0.142; // [29] Dainotti et al.

H0[i]=73.192; PT[i]=+0.265; NT[i++]=-0.265; // [29] Dainotti et al.

H0[i]=73.678; PT[i]=+0.223; NT[i++]=-0.223; // [29] Dainotti et al.

H0[i]=71.8; PT[i]=+3.9; NT[i++]=-3.3; // [30] Denzel et al.

H0[i]=72.94; PT[i]=+0.75; NT[i++]=-0.75; // [31] Di Valentino

H0[i]=72.7; PT[i]=+1.1; NT[i++]=-1.1; // [31] Di Valentino

H0[i]=68.8; PT[i]=+45.7; NT[i++]=-25.5; // [32] Gayathri et al.

H0[i]=62.3; PT[i]=+9.1; NT[i++]=-9.1; // [33] Hagstotz et al.

H0[i]=70.5; PT[i]=+2.37; NT[i++]=-2.37; // [34] Kethan et al.

H0[i]=72.86; PT[i]=+0.036; NT[i++]=-0.036; //Mercier (this document)

H0[i]=72.105; PT[i]=+0.036; NT[i++]=-0.036; //Mercier (this document)

H0[i]=68.3; PT[i]=+4.6; NT[i++]=-4.6; // [35] Mukherjee et al.

H0[i]=70; PT[i]=+0.5; NT[i++]=-0.5; // [36] Park et al.

H0[i]=65.1; PT[i]=+3; NT[i++]=-5.4; // [37] Philcox et al.

H0[i]=65.6; PT[i]=+3.4; NT[i++]=-3.5; // [37] Philcox et al.

H0[i]=70.6; PT[i]=+3.7; NT[i++]=-5; // [37] Philcox et al.

H0[i]=78.3; PT[i]=+2.9; NT[i++]=-2.9; // [38] Qi et al.

H0[i]=73.6; PT[i]=+1.8; NT[i++]=-1.6; // [38] Qi et al.

H0[i]=73; PT[i]=+1.4; NT[i++]=-1.4; // [39] Riess et al.

H0[i]=73.2; PT[i]=+1.3; NT[i++]=-1.3; // [39] Riess et al.

H0[i]=72.1; PT[i]=+2; NT[i++]=-2; // [40] Soltis et al.

H0[i]=69.5; PT[i]=+4; NT[i++]=-4; // [41] Wang & Giannios

H0[i]=71; PT[i]=+20; NT[i++]=-20; // [42] Zhang et al.

H0[i]=67.4; PT[i]=+0.5; NT[i++]=-0.5; // [43] Aghanim et al.

H0[i]=67.73; PT[i]=+0.41; NT[i++]=-0.41; // [44] Benevento

H0[i]=68.22; PT[i]=+0.39; NT[i++]=-0.39; // [44] Benevento

H0[i]=72.5; PT[i]=+1.85; NT[i++]=-1.85; // [44] Benevento

H0[i]=69.17; PT[i]=+1.09; NT[i++]=-1.09; // [44] Benevento

H0[i]=74.5; PT[i]=+5.6; NT[i++]=-6.1; // [45] Birrer et al.

H0[i]=67.4; PT[i]=+4.1; NT[i++]=-3.2; // [45] Birrer et al.

H0[i]=75.35; PT[i]=+1.68; NT[i++]=-1.68; // [46] Camarena & Marra

H0[i]=74; PT[i]=+0.625; NT[i++]=-0.625; // [47] Chang & Zhu

H0[i]=73.8; PT[i]=+6.3; NT[i++]=-5.8; // [48] Coughlin et al.

H0[i]=71.2; PT[i]=+3.2; NT[i++]=-3.1; // [48] Coughlin et al.

H0[i]=72.4; PT[i]=+1.4; NT[i++]=-1.4; // [49] D'Agostino

H0[i]=71.5; PT[i]=+1.3; NT[i++]=-1.3; // [49] D'Agostino

H0[i]=71.54; PT[i]=+1.78; NT[i++]=-1.78; // [50] Dai WM et al.

H0[i]=73.12; PT[i]=+1.14; NT[i++]=-1.14; // [50] Dai WM et al.

H0[i]=66.2; PT[i]=+4.4; NT[i++]=-4.2; // [51] Dietrich et al.

H0[i]=69.9; PT[i]=+0.84; NT[i++]=-0.86; // [52] Gonzalez et al.

H0[i]=71; PT[i]=+4; NT[i++]=-4; // [53] González-Serrena et al.

H0[i]=74.62; PT[i]=+12.35; NT[i++]=-11.34; // [54] Haboury

H0[i]=71.89; PT[i]=+11.02; NT[i++]=-10.17; // [54] Haboury

H0[i]=76.44; PT[i]=+55.76; NT[i++]=-50.17; // [54] Haboury

H0[i]=50.9; PT[i]=+31.1; NT[i++]=-31.9; // [54] Haboury

H0[i]=50.81; PT[i]=+28.19; NT[i++]=-27.81; // [54] Haboury

H0[i]=71; PT[i]=+2; NT[i++]=-3; // [55] Harvey

H0[i]=65.9; PT[i]=+1.5; NT[i++]=-1.5; // [56] Holanda et al.

H0[i]=65.9; PT[i]=+4.4; NT[i++]=-4; // [56] Holanda et al.

H0[i]=64.3; PT[i]=+4.5; NT[i++]=-4.4; // [56] Holanda et al.

H0[i]=66.8; PT[i]=+13.4; NT[i++]=-9.2; // [57] Howlett & Davis

H0[i]=64.8; PT[i]=+7.3; NT[i++]=-7.2; // [57] Howlett & Davis

H0[i]=75.8; PT[i]=+5.2; NT[i++]=-4.9; // [58] Jaeger et al.

H0[i]=65.8; PT[i]=+3.5; NT[i++]=-3.5; // [59] Kim et al.

H0[i]=72.3; PT[i]=+1.4; NT[i++]=-1.4; // [60] Kreisch et al.

H0[i]=71.5; PT[i]=+11.9; NT[i++]=-10.6; // [61] Li & Zhang

H0[i]=74.7; PT[i]=+5.8; NT[i++]=-5.8; // [62] Lombriser

H0[i]=72.06; PT[i]=+0.09; NT[i++]=-0.09; // [7] Mercier

H0[i]=74; PT[i]=+1.6; NT[i++]=-1.6; // [63] Millon et al.

H0[i]=74.2; PT[i]=+1.7; NT[i++]=-1.8; // [63] Millon et al.

H0[i]=50.4; PT[i]=+28.1; NT[i++]=-19.5; // [64] Mukherjee et al.

H0[i]=62.2; PT[i]=+29.5; NT[i++]=-19.7; // [64] Mukherjee et al.

H0[i]=43.1; PT[i]=+24.6; NT[i++]=-11.4; // [64] Mukherjee et al.

H0[i]=67.6; PT[i]=+4.3; NT[i++]=-4.2; // [64] Mukherjee et al.

H0[i]=68.6; PT[i]=+14; NT[i++]=-8.5; // [65] Nicolaou et al.

H0[i]=69.6; PT[i]=+1; NT[i++]=-1.3; // [66] Niedermann and Sloth

H0[i]=71.4; PT[i]=+1; NT[i++]=-1; // [66] Niedermann and Sloth

H0[i]=72; PT[i]=+12; NT[i++]=-8.2; // [67] Palmese et al.

H0[i]=69.03; PT[i]=+0.87; NT[i++]=-0.87; // [68] Pandey et al.

H0[i]=70.6; PT[i]=+1.1; NT[i++]=-1.1; // [68] Pandey et al.

H0[i]=68.44; PT[i]=+0.52; NT[i++]=-0.52; // [68] Pandey et al.

H0[i]=68.1; PT[i]=+0.58; NT[i++]=-0.58; // [68] Pandey et al.

H0[i]=73.9; PT[i]=+3; NT[i++]=-3; // [69] Pesce et al.

H0[i]=68.6; PT[i]=+1.8; NT[i++]=-1.8; // [70] Pogosian et al.

H0[i]=74.03; PT[i]=+1.42; NT[i++]=-1.42; // [71] Rui-Yun et al.

H0[i]=75.1; PT[i]=+2.3; NT[i++]=-2.3; // [72] Schombert et al.

H0[i]=74.2; PT[i]=+2.7; NT[i++]=-3; // [73] Shajib et al.

H0[i]=67.52; PT[i]=+0.96; NT[i++]=-0.95; // [74] Sharov & Sinyakov

H0[i]=70.87; PT[i]=+1.63; NT[i++]=-1.62; // [74] Sharov & Sinyakov

H0[i]=69; PT[i]=+29; NT[i++]=-14; // [75] Vasylyev & Filippenko

H0[i]=67; PT[i]=+41; NT[i++]=-26; // [75] Vasylyev & Filippenko

H0[i]=71; PT[i]=+34; NT[i++]=-30; // [75] Vasylyev & Filippenko

H0[i]=70; PT[i]=+29; NT[i++]=-18; // [75] Vasylyev & Filippenko

H0[i]=72.3; PT[i]=+2.9; NT[i++]=-2.8; // [76] Vogl

H0[i]=75.3; PT[i]=+3; NT[i++]=-2.9; // [77] Wei & Melia

H0[i]=75.3; PT[i]=+1.9; NT[i++]=-1.9; // [77] Wei & Melia

H0[i]=67.9; PT[i]=+1.1; NT[i++]=-1.3; // [78] Wu et al.

H0[i]=72; PT[i]=+2.1; NT[i++]=-2.5; // [78] Wu et al.

H0[i]=73.65; PT[i]=+1.95; NT[i++]=-2.26; // [79] Yang et al.

H0[i]=67.95; PT[i]=+0.78; NT[i++]=-1.03; // [80] Zhang & Huang

H0[i]=69.81; PT[i]=+2.22; NT[i++]=-2.7; // [80] Zhang & Huang

H0[i]=66.75; PT[i]=+3.42; NT[i++]=-4.23; // [80] Zhang & Huang

H0[i]=70.75; PT[i]=+1.55; NT[i++]=-1.55; // [81] Agrawal

H0[i]=73.7; PT[i]=+1.4; NT[i++]=-1.4; // [82] Anderson

H0[i]=72.5; PT[i]=+2.1; NT[i++]=-2.3; // [83] Birrer

H0[i]=67.4; PT[i]=+0.5; NT[i++]=-0.5; // [84] Chang et al.

H0[i]=82.8; PT[i]=+9.4; NT[i++]=-8.3; // [85] Chen et al.

H0[i]=70.1; PT[i]=+5.3; NT[i++]=-4.5; // [85] Chen et al.

H0[i]=77; PT[i]=+4; NT[i++]=-4.6; // [85] Chen et al.

H0[i]=75.6; PT[i]=+3.2; NT[i++]=-3; // [85] Chen et al.

H0[i]=76.8; PT[i]=+2.6; NT[i++]=-2.6; // [85] Chen et al.

H0[i]=75.7; PT[i]=+4.5; NT[i++]=-4.4; // [86] Collett

H0[i]=76.8; PT[i]=+4.2; NT[i++]=-3.8; // [86] Collett

H0[i]=74.2; PT[i]=+3; NT[i++]=-2.9; // [86] Collett

H0[i]=67.6; PT[i]=+1.1; NT[i++]=-1.1; // [87] Cuceu et al.

H0[i]=67.4; PT[i]=+6; NT[i++]=-6.2; // [88] Domínguez

H0[i]=66.6; PT[i]=+1.6; NT[i++]=-1.6; // [88] Domínguez

H0[i]=70.3; PT[i]=+1.36; NT[i++]=-1.35; // [89] Dutta et al.

H0[i]=77; PT[i]=+37; NT[i++]=-18; // [90] Fishbach et al.

H0[i]=76; PT[i]=+19; NT[i++]=-13; // [90] Fishbach et al.

H0[i]=69.8; PT[i]=+0.8; NT[i++]=-0.8; // [91] Freedman et al.

H0[i]=68.09; PT[i]=+0.45; NT[i++]=-0.45; // [92] Guo et al.

H0[i]=69.34; PT[i]=+0.93; NT[i++]=-0.93; // [92] Guo et al.

H0[i]=69.67; PT[i]=+0.95; NT[i++]=-0.94; // [92] Guo et al.

H0[i]=69.36; PT[i]=+0.82; NT[i++]=-0.82; // [92] Guo et al.

H0[i]=69.25; PT[i]=+0.99; NT[i++]=-0.99; // [92] Guo et al.

H0[i]=74; PT[i]=+16; NT[i++]=-8; // [93] Hotokezaka et al.

H0[i]=70.3; PT[i]=+5.3; NT[i++]=-5; // [93] Hotokesaka et al.

H0[i]=82.4; PT[i]=+8.4; NT[i++]=-8.3; // [94] Jee et al.

H0[i]=67; PT[i]=+3; NT[i++]=-3; // [95] Kozmanyan et al.

H0[i]=72.2; PT[i]=+2.1; NT[i++]=-2.1; // [96] Liao et al.

H0[i]=73; PT[i]=+2.8; NT[i++]=-3; // [96] Liao et al.

H0[i]=67.8; PT[i]=+1.3; NT[i++]=-1.3; // [97] MacAulay et al.

H0[i]=67.37; PT[i]=+0.62; NT[i++]=-0.62; // [98] Martinelli

H0[i]=68.8; PT[i]=+1.6; NT[i++]=-1.6; // [98] Martinelli

H0[i]=73.9; PT[i]=+2.5; NT[i++]=-2.5; // [98] Martinelli

H0[i]=67.68; PT[i]=+0.46; NT[i++]=-0.46; // [98] Martinelli

H0[i]=68.4; PT[i]=+1; NT[i++]=-1; // [98] Martinelli

H0[i]=69.2; PT[i]=+1.5; NT[i++]=-1.5; // [98] Martinelli

H0[i]=67.51; PT[i]=+0.61; NT[i++]=-0.61; // [98] Martinelli

H0[i]=68.9; PT[i]=+1.1; NT[i++]=-1.1; // [98] Martinelli

H0[i]=72.1; PT[i]=+2.1; NT[i++]=-1.8; // [98] Martinelli

H0[i]=67.75; PT[i]=+0.46; NT[i++]=-0.46; // [98] Martinelli

H0[i]=68.59; PT[i]=+0.86; NT[i++]=-0.86; // [98] Martinelli

H0[i]=69.6; PT[i]=+1.3; NT[i++]=-1.3; // [98] Martinelli

H0[i]=71.505; PT[i]=+0.03; NT[i++]=-0.03; // [2] Mercier

H0[i]=69; PT[i]=+1.7; NT[i++]=-1.7; // [99] Park & Ratra

H0[i]=69.8; PT[i]=+1.8; NT[i++]=-1.8; // [99] Park & Ratra

H0[i]=68.9; PT[i]=+1.7; NT[i++]=-1.7; // [99] Park & Ratra

H0[i]=70.1; PT[i]=+1.9; NT[i++]=-1.9; // [99] Park & Ratra

H0[i]=68.5; PT[i]=+1.8; NT[i++]=-1.8; // [99] Park & Ratra

H0[i]=69.6; PT[i]=+1.9; NT[i++]=-1.9; // [99] Park & Ratra

H0[i]=72; PT[i]=+1.9; NT[i++]=-1.9; // [100] Reid

H0[i]=73.5; PT[i]=+1.4; NT[i++]=-1.4; // [100] Reid

H0[i]=74.22; PT[i]=+1.82; NT[i++]=-1.82; // [101] Riess et al.

H0[i]=74.03; PT[i]=+1.42; NT[i++]=-1.42; // [101] Riess et al.

H0[i]=72.8; PT[i]=+1.1; NT[i++]=-1.1; // [102] Riess

H0[i]=74.3; PT[i]=+1; NT[i++]=-1; // [102] Riess

H0[i]=71.6; PT[i]=+3.8; NT[i++]=-4.9; // [103] Rusu et al.

H0[i]=67.99; PT[i]=+0.91; NT[i++]=-0.88; // [104] Ryan

H0[i]=68.24; PT[i]=+2.39; NT[i++]=-2.33; // [104] Ryan

H0[i]=66.79; PT[i]=+2.6; NT[i++]=-2.32; // [104] Ryan

H0[i]=66.8; PT[i]=+2.5; NT[i++]=-2.3; // [104] Ryan

H0[i]=66.13; PT[i]=+1.38; NT[i++]=-2.09; // [104] Ryan

H0[i]=67.1; PT[i]=+2.4; NT[i++]=-2.3; // [104] Ryan

H0[i]=68.44; PT[i]=+0.7; NT[i++]=-0.69; // [104] Ryan

H0[i]=69.32; PT[i]=+1.42; NT[i++]=-1.42; // [104] Ryan

H0[i]=68; PT[i]=+2.27; NT[i++]=-1.94; // [104] Ryan

H0[i]=66.6; PT[i]=+2.2; NT[i++]=-1.9; // [104] Ryan

H0[i]=67.19; PT[i]=+1; NT[i++]=-1.6; // [104] Ryan

H0[i]=66.8; PT[i]=+1.8; NT[i++]=-1.7; // [104] Ryan

H0[i]=63.13; PT[i]=+6.48; NT[i++]=-6.48; // [105] Saha & Sahoo

H0[i]=74.2; PT[i]=+2.7; NT[i++]=-3; // [106] Shajib et al.

H0[i]=75; PT[i]=+40; NT[i++]=-32; // [107] Soares-Santos

H0[i]=78; PT[i]=+96; NT[i++]=-24; // [107] Soares-Santos

H0[i]=73.1; PT[i]=+0.7; NT[i++]=-0.7; // [108] Taubenberger et al.

H0[i]=68; PT[i]=+14; NT[i++]=-7; // [109] Tiwari et al.

H0[i]=68; PT[i]=+18; NT[i++]=-8; // [109] Tiwari et al.

H0[i]=73.9; PT[i]=+1; NT[i++]=-1; // [110] Verde et al.

H0[i]=72.5; PT[i]=+1.2; NT[i++]=-1.2; // [110] Verde et al.

H0[i]=73.3; PT[i]=+1.7; NT[i++]=-1.8; // [111] Wong et al.

H0[i]=72.4; PT[i]=+2; NT[i++]=-2; // [112] Yuan et al.

H0[i]=68.36; PT[i]=+0.53; NT[i++]=-0.52; // [113] Zhang & Huang

H0[i]=64.9; PT[i]=+4.6; NT[i++]=-4.3; // [114] Zeng and Yan

H0[i]=67.4; PT[i]=+1.1; NT[i++]=-1.2; // [115] Abbott et al.

H0[i]=69.3; PT[i]=+0.4; NT[i++]=-0.6; // [115] Abbott et al.

H0[i]=73.24; PT[i]=+1.74; NT[i++]=-1.74; // [116] Benetti et al.

H0[i]=72.5; PT[i]=+2.1; NT[i++]=-2.1; // [117] Bolejko

H0[i]=68.1; PT[i]=+2; NT[i++]=-2; // [117] Bolejko

H0[i]=76; PT[i]=+8; NT[i++]=-8; // [118] Braatz

H0[i]=69.3; PT[i]=+4.2; NT[i++]=-4.2; // [118] Braatz

H0[i]=71.9; PT[i]=+7.1; NT[i++]=-7.1; // [119] Cantiello et al.

H0[i]=73.24; PT[i]=+1.74; NT[i++]=-1.74; // [120] Chen

H0[i]=67.4; PT[i]=+0.5; NT[i++]=-0.5; // [120] Chen

H0[i]=73.24; PT[i]=+1.74; NT[i++]=-1.74; // [121] Choudhury & Choubey

H0[i]=72.8; PT[i]=+1.6; NT[i++]=-1.6; // [122] Dhawan et al.

H0[i]=55; PT[i]=+7; NT[i++]=-20;// [123] Di Valentino & Melchiorri

H0[i]=67.06; PT[i]=+1.68; NT[i++]=-1.68; // [124] Gomez-Valent

H0[i]=68.9; PT[i]=+1.96; NT[i++]=-1.96; // [124] Gomez-Valent

H0[i]=68.45; PT[i]=+2; NT[i++]=-2; // [124] Gomez-Valent

H0[i]=73.5; PT[i]=+4.6; NT[i++]=-4.7; // [125] Grillo

H0[i]=72.8; PT[i]=+4.3; NT[i++]=-4.1; // [125] Grillo

H0[i]=69.8; PT[i]=+5.3; NT[i++]=-4.1; // [125] Grillo

H0[i]=70.38; PT[i]=+0.6; NT[i++]=-0.6; // [126] Hoeneisen et al.

H0[i]=71.17; PT[i]=+1.66; NT[i++]=-1.66; // [127] Lee & Jang

H0[i]=73.52; PT[i]=+1.62; NT[i++]=-1.62; // [128] Riess et al.

H0[i]=73.83; PT[i]=+1.48; NT[i++]=-1.48; // [128] Riess et al.

H0[i]=73.48; PT[i]=+1.66; NT[i++]=-1.66; // [129] Riess et al.

H0[i]=74.4; PT[i]=+4.9; NT[i++]=-4.9; // [130] Van Putten

H0[i]=74.5; PT[i]=+7.3; NT[i++]=-7.3; // [130] Van Putten

H0[i]=74.9; PT[i]=+2.6; NT[i++]=-2.6; // [130] Van Putten

H0[i]=66.8; PT[i]=+1.9; NT[i++]=-1.9; // [130] Van Putten

H0[i]=73.75; PT[i]=+1.44; NT[i++]=-1.44; // [130] Van Putten

H0[i]=70; PT[i]=+12; NT[i++]=-8; // [131] Vitale

H0[i]=67; PT[i]=+4; NT[i++]=-4; // [132] Yu et al.

H0[i]=67.498; PT[i]=+7.97; NT[i++]=-3.278; // [133] Zhang

H0[i]=70; PT[i]=+12; NT[i++]=-8; // [134] Abbott

H0[i]=72.5; PT[i]=+2.5; NT[i++]=-8; // [135] Bethapudi & Desai

H0[i]=71.9; PT[i]=+2.4; NT[i++]=-3; // [136] Bonvin et al.

H0[i]=69.2; PT[i]=+1.4; NT[i++]=-2.2; // [136] Bonvin et al.

H0[i]=79; PT[i]=+4.4; NT[i++]=-4.2; // [136] Bonvin et al.

H0[i]=73.75; PT[i]=+2.11; NT[i++]=-2.11; // [137] Cardona

H0[i]=67.81; PT[i]=+0.92; NT[i++]=-0.92; // [137] Cardona

H0[i]=66.93; PT[i]=+0.62; NT[i++]=-0.62; // [137] Cardona

H0[i]=73.46; PT[i]=+1.4; NT[i++]=-1.4; // [137] Cardona

H0[i]=68.3; PT[i]=+2.7; NT[i++]=-2.6; // [138] Chen Yun et al.

H0[i]=68.4; PT[i]=+2.9; NT[i++]=-3.3; // [138] Chen Yun et al.

H0[i]=65; PT[i]=+6.6; NT[i++]=-6.6; // [138] Chen Yun et al.

H0[i]=67.9; PT[i]=+2.4; NT[i++]=-2.4; // [138] Chen Yun et al.

H0[i]=68; PT[i]=+2.8; NT[i++]=-2.8; // [139] Farooq

H0[i]=73.24; PT[i]=+1.74; NT[i++]=-1.74; // [139] Farooq

H0[i]=72.72; PT[i]=+1.67; NT[i++]=-1.67; // [140] Feeney et al.

H0[i]=73.15; PT[i]=+1.78; NT[i++]=-1.78; // [140] Feeney et al.

H0[i]=67.6; PT[i]=+0.7; NT[i++]=-0.6; // [141] Grieb et al.

H0[i]=73; PT[i]=+1.75; NT[i++]=-1.75; // [142] Guo & Zhang

H0[i]=73.24; PT[i]=+1.74; NT[i++]=-1.74; // [143] Hjorth et al.

H0[i]=69.13; PT[i]=+0.24; NT[i++]=-0.24; // [144] Huang and Huang

H0[i]=71.66; PT[i]=+1.8; NT[i++]=-1.8; // [145] Jang & Lee

H0[i]=73.72; PT[i]=+2.03; NT[i++]=-2.03; // [145] Jang & Lee

H0[i]=71.17; PT[i]=+1.66; NT[i++]=-1.66; // [145] Jang & Lee

H0[i]=66.2; PT[i]=+8.9; NT[i++]=-8.9; // [146] Pritychenko

H0[i]=67.2; PT[i]=+6.9; NT[i++]=-6.9; // [146] Pritychenko

H0[i]=69.13; PT[i]=+2.34; NT[i++]=-2.34; // [147] Wang et al.

H0[i]=73.24; PT[i]=+1.74; NT[i++]=-1.74; // [148] Wei & Wu

H0[i]=69.6; PT[i]=+0.7; NT[i++]=-0.7; // [148] Wei & Wu

H0[i]=73.1; PT[i]=+5.7; NT[i++]=-6; // [149] Wong et al.

H0[i]=72.5; PT[i]=+3.1; NT[i++]=-3.1; // [150] Zhang et al.

H0[i]=67.8; PT[i]=+0.9; NT[i++]=-0.9; // [151] Ade et al.

H0[i]=66; PT[i]=+6; NT[i++]=-6; // [152] Gao et al.

H0[i]=70.1; PT[i]=+0.34; NT[i++]=-0.34; // [153] Ichiki et al.

H0[i]=66.5; PT[i]=+1.8; NT[i++]=-1.8; // [154] Ludovic et al.

H0[i]=64.2; PT[i]=+1.9; NT[i++]=-1.9; // [154] Ludovic et al.

H0[i]=91.8; PT[i]=+5.3; NT[i++]=-5.3; // [155] Moresco et al.

H0[i]=72.25; PT[i]=+2.51; NT[i++]=-2.51; // [156] Riess et al.

H0[i]=72.04; PT[i]=+2.67; NT[i++]=-2.67; // [156] Riess et al.

H0[i]=76.18; PT[i]=+2.37; NT[i++]=-2.37; // [156] Riess et al.

H0[i]=74.5; PT[i]=+3.27; NT[i++]=-3.27; // [156] Riess et al.

H0[i]=73.24; PT[i]=+1.74; NT[i++]=-1.74; // [156] Riess et al.

H0[i]=76.2; PT[i]=+3.4; NT[i++]=-3.4; // [157] Tully et al.

H0[i]=75; PT[i]=+2; NT[i++]=-2; // [157] Tully et al.

H0[i]=68.17; PT[i]=+1.55; NT[i++]=-1.56; // [158] Cheng & Qing Guo

H0[i]=68.11; PT[i]=+1.69; NT[i++]=-1.69; // [158] Cheng & Qing Guo

H0[i]=68.11; PT[i]=+0.86; NT[i++]=-0.86; // [158] Cheng & Qing Guo

H0[i]=67.7; PT[i]=+1.1; NT[i++]=-1.1; // [159] Cuesta et al.

H0[i]=69.8; PT[i]=+2.6; NT[i++]=-2.6; // [160] Jang & Lee

H0[i]=72.2; PT[i]=+3.3; NT[i++]=-3.3; // [160] Jang & Lee

H0[i]=68.1; PT[i]=+5.9; NT[i++]=-5.9; // [161] Kumar et al.

H0[i]=73; PT[i]=+26; NT[i++]=-22; // [162] Kuo et al.

H0[i]=70.6; PT[i]=+2.6; NT[i++]=-2.6; // [163] Rigault et al.

H0[i]=68.8; PT[i]=+3.3; NT[i++]=-3.3; // [163] Rigault et al.

H0[i]=67.3; PT[i]=+1.2; NT[i++]=-1.2; // [164] Ade et al.

H0[i]=70.8; PT[i]=+2.4; NT[i++]=-2.4; // [165] Ben-Dayan et al.

H0[i]=69.6; PT[i]=+0.7; NT[i++]=-0.7; // [166] Bennett et al.

H0[i]=64.9; PT[i]=+4.2; NT[i++]=-4.2; // [167] Busti et al.

H0[i]=72.5; PT[i]=+2.5; NT[i++]=-2.5; // [168] Efstathiou

H0[i]=70.6; PT[i]=+3.3; NT[i++]=-3.3; // [168] Efstathiou

H0[i]=74.1; PT[i]=+2.2; NT[i++]=-2.2; // [169] Lima & Cunha

H0[i]=70; PT[i]=+2.2; NT[i++]=-2.2; // [170] Bennett et al.

H0[i]=69.32; PT[i]=+0.8; NT[i++]=-0.8; // [170] Bennett et al.

H0[i]=68; PT[i]=+4.8; NT[i++]=-4.8; // [171] Braatz et al.

H0[i]=68; PT[i]=+2.8; NT[i++]=-2.8; // [172] Farooq & Bathra

H0[i]=73.8; PT[i]=+2.4; NT[i++]=-2.4; // [172] Farooq & Bathra

H0[i]=69.7; PT[i]=+2.4; NT[i++]=-2.4; // [173] Hinshaw et al.

H0[i]=70.4; PT[i]=+2.5; NT[i++]=-2.5; // [173] Hinshaw et al.

H0[i]=69.33; PT[i]=+0.88; NT[i++]=-0.88; // [173] Hinshaw et al.

H0[i]=70.2; PT[i]=+1.4; NT[i++]=-1.4; // [173] Hinshaw et al.

H0[i]=70; PT[i]=+3; NT[i++]=-3; // [174] Humphreys et al.

H0[i]=68; PT[i]=+9; NT[i++]=-9; // [175] Kuo et al.

H0[i]=49.97; PT[i]=+0.19; NT[i++]=-0.19; // [176] Pietrzynski et al.

H0[i]=68.9; PT[i]=+7.1; NT[i++]=-7.1; // [177] Reid et al.

H0[i]=72.1; PT[i]=+3.2; NT[i++]=-2.3; // [178] Salvatelli et al.

H0[i]=74.1; PT[i]=+2.1; NT[i++]=-2.1; // [179] Scowcroft et al.

H0[i]=69; PT[i]=+6; NT[i++]=-6; // [180] Sereno and Pacificz

H0[i]=80; PT[i]=+5.8; NT[i++]=-5.7; // [181] Suyu et al.

H0[i]=75.2; PT[i]=+4.4; NT[i++]=-4.2; // [181] Suyu et al.

H0[i]=73.1; PT[i]=+2.4; NT[i++]=-3.6; // [181] Suyu et al.

H0[i]=74.4; PT[i]=+3; NT[i++]=-3; // [182] Tully et al.

H0[i]=71.3; PT[i]=+2; NT[i++]=-2; // [183] Xia et al.

H0[i]=73.8; PT[i]=+2.4; NT[i++]=-2.4; // [184] Calabrese et al.

H0[i]=68; PT[i]=+2.8; NT[i++]=-2.8; // [184] Calabrese et al.

H0[i]=69.7; PT[i]=+2.5; NT[i++]=-2.5; // [184] Calabrese et al.

H0[i]=74.3; PT[i]=+3.1; NT[i++]=-3.1; // [185] Chavez

H0[i]=67; PT[i]=+3.2; NT[i++]=-3.2; // [186] Colless et al.

H0[i]=74.3; PT[i]=+3; NT[i++]=-3; // [187] Freedman et al.

H0[i]=70.2; PT[i]=+0.14; NT[i++]=-0.14; // [188] Pozzo

H0[i]=75.4; PT[i]=+2.9; NT[i++]=-2.9; // [189] Riess et al.

H0[i]=56; PT[i]=+2; NT[i++]=-2; // [190] Wang

H0[i]=68; PT[i]=+5.5; NT[i++]=-5.5; // [191] Chen & Ratra

H0[i]=67; PT[i]=+3.2; NT[i++]=-3.2; // [192] Beutler et al.

H0[i]=71; PT[i]=+2.5; NT[i++]=-2.5; // [193] Jarosik et al.

H0[i]=70.4; PT[i]=+1.3; NT[i++]=-1.4; // [193] Jarosik et al.

H0[i]=74.8; PT[i]=+3.1; NT[i++]=-3.1; // [194] Riess et al.

H0[i]=74.4; PT[i]=+2.5; NT[i++]=-2.5; // [194] Riess et al.

H0[i]=73.8; PT[i]=+2.4; NT[i++]=-2.4; // [194] Riess et al.

H0[i]=73; PT[i]=+2; NT[i++]=-2; // [195] Freedman & Madore

H0[i]=66; PT[i]=+6; NT[i++]=-4; // [196] Paraficz et Hjorth

H0[i]=76; PT[i]=+3; NT[i++]=-3; // [196] Paraficz et Hjorth

H0[i]=70.6; PT[i]=+3.1; NT[i++]=-3.1; // [197] Suyu et al.

H0[i]=69.7; PT[i]=+4.9; NT[i++]=-5; // [197] Suyu et al.

H0[i]=70.5; PT[i]=+1.3; NT[i++]=-1.3; // [198] Hinshaw et al.

H0[i]=71.9; PT[i]=+2.6; NT[i++]=-2.7; // [198] Hinshaw et al.

H0[i]=70.5; PT[i]=+1.3; NT[i++]=-1.3; // [199] Komatsu et al.

H0[i]=70.4; PT[i]=+1.4; NT[i++]=-1.4; // [199] Komatsu et al.

H0[i]=70.9; PT[i]=+1.3; NT[i++]=-1.3; // [199] Komatsu et al.

H0[i]=70.1; PT[i]=+1.3; NT[i++]=-1.3; // [199] Komatsu et al.

H0[i]=74.2; PT[i]=+3.6; NT[i++]=-3.6; // [200] Riess et al.

H0[i]=84.2; PT[i]=+6; NT[i++]=-6; // [201] Russell

H0[i]=83.4; PT[i]=+8; NT[i++]=-8; // [201] Russell

H0[i]=88; PT[i]=+6; NT[i++]=-6; // [201] Russell

H0[i]=61.7; PT[i]=+1.2; NT[i++]=-1.1; // [202] Leith et al.

H0[i]=67; PT[i]=+13; NT[i++]=-10; // [203] Vuissoz et al.

H0[i]=63; PT[i]=+7; NT[i++]=-3; // [203] Vuissoz et al.

H0[i]=70; PT[i]=+6; NT[i++]=-6; // [204] Oguri

H0[i]=68; PT[i]=+6; NT[i++]=-6; // [204] Oguri

H0[i]=73.5; PT[i]=+3.2; NT[i++]=-3.2; // [205] Spergel et al.

H0[i]=73.2; PT[i]=+3.1; NT[i++]=-3.2; // [205] Spergel et al.

H0[i]=70.4; PT[i]=+1.5; NT[i++]=-1.6; // [205] Spergel et al.

H0[i]=76.9; PT[i]=+3.9; NT[i++]=-3.4; // [206] Bonamente et al.

H0[i]=73.7; PT[i]=+4.6; NT[i++]=-3.8; // [206] Bonamente et al.

H0[i]=77.6; PT[i]=+4.8; NT[i++]=-4.3; // [206] Bonamente et al.

H0[i]=70.8; PT[i]=+1.9; NT[i++]=-1.8; // [207] Hütsi

H0[i]=74.92; PT[i]=+2.28; NT[i++]=-2.28; // [208] Ngeow and Kanbur

H0[i]=74.37; PT[i]=+2.27; NT[i++]=-2.27; // [208] Ngeow and Kanbur

H0[i]=62.3; PT[i]=+1.3; NT[i++]=-1.3; // [209] Sandage et al.

H0[i]=60.9; PT[i]=+1.3; NT[i++]=-1.3; // [209] Sandage et al.

H0[i]=60.7; PT[i]=+1.5; NT[i++]=-1.5; // [209] Sandage et al.

H0[i]=72; PT[i]=+6; NT[i++]=-6; // [210] Wang et al.

H0[i]=73.2; PT[i]=+7; NT[i++]=-7; // [211] Gibson & Brook

H0[i]=75; PT[i]=+7; NT[i++]=-7; // [212] Hamuy

H0[i]=65; PT[i]=+12; NT[i++]=-12; // [212] Hamuy

H0[i]=58; PT[i]=+2; NT[i++]=-2; // [213] Magain

H0[i]=58; PT[i]=+2; NT[i++]=-2; // [214] Olivares et al.

H0[i]=73; PT[i]=+4; NT[i++]=-4; // [215] Riess

H0[i]=69; PT[i]=+8; NT[i++]=-8; // [216] Schmidt et al.

H0[i]=66; PT[i]=+8; NT[i++]=-8; // [217] Stritzinger et al.

H0[i]=78; PT[i]=+9; NT[i++]=-9; // [217] Stritzinger et al.

H0[i]=67; PT[i]=+30; NT[i++]=-18; // [218] Udomprasert et al.

H0[i]=64; PT[i]=+7; NT[i++]=-4; // [219] Boffi & Riess

H0[i]=33; PT[i]=+5; NT[i++]=-5; // [220] Dumin

H0[i]=69; PT[i]=+12; NT[i++]=-12; // [221] Jimenez et al.

H0[i]=75; PT[i]=+7; NT[i++]=-6; // [222] Koopmans

H0[i]=70; PT[i]=+7; NT[i++]=-7; // [223] Mei et al.

H0[i]=68; PT[i]=+6; NT[i++]=-6; // [223] Mei et al.

H0[i]=68; PT[i]=+5; NT[i++]=-5; // [223] Mei et al.

H0[i]=71; PT[i]=+4; NT[i++]=-4; // [223] Mei et al.

H0[i]=77; PT[i]=+19; NT[i++]=-15; // [224] Saunders et al.

H0[i]=85; PT[i]=+20; NT[i++]=-17; // [224] Saunders et al.

H0[i]=72; PT[i]=+5; NT[i++]=-5; // [225] Spergel et al.

H0[i]=71; PT[i]=+4; NT[i++]=-3; // [225] Spergel et al.

H0[i]=63; PT[i]=+2; NT[i++]=-2; // [226] Fassnacht et al.

H0[i]=72; PT[i]=+8; NT[i++]=-8; // [227] Freedman

H0[i]=57; PT[i]=+23; NT[i++]=-16; // [228] Grainge et al.

H0[i]=48; PT[i]=+7; NT[i++]=-4; // [229] Kochanek

H0[i]=71; PT[i]=+6; NT[i++]=-6; // [229] Kochanek

H0[i]=72; PT[i]=+8; NT[i++]=-8; // [229] Kochanek

H0[i]=62; PT[i]=+7; NT[i++]=-7; // [229] Kochanek

H0[i]=75; PT[i]=+8; NT[i++]=-8;// [230] Tikhonov & Galazoutdinova

H0[i]=81; PT[i]=+5; NT[i++]=-5;// [230] Tikhonov & Galazoutdinova

H0[i]=59; PT[i]=+15; NT[i++]=-10; // [231] Treu & Koopmans

H0[i]=71; PT[i]=+2; NT[i++]=-2; // [232] Freedman et al.

H0[i]=71; PT[i]=+3; NT[i++]=-3; // [232] Freedman et al.

H0[i]=70; PT[i]=+5; NT[i++]=-5; // [232] Freedman et al.

H0[i]=72; PT[i]=+9; NT[i++]=-9; // [232] Freedman et al.

H0[i]=82; PT[i]=+6; NT[i++]=-6; // [232] Freedman et al.

H0[i]=72; PT[i]=+8; NT[i++]=-8; // [232] Freedman et al.

H0[i]=65; PT[i]=+5; NT[i++]=-5; // [233] Itoh

H0[i]=76; PT[i]=+1.3; NT[i++]=-1.3; // [234] Jensen et al.

H0[i]=72; PT[i]=+2.3; NT[i++]=-2.3; // [234] Jensen et al.

H0[i]=65; PT[i]=+5; NT[i++]=-5; // [235] Koopmans et al.

H0[i]=71; PT[i]=+8; NT[i++]=-8; // [236] Liu & Graham

H0[i]=64; PT[i]=+14; NT[i++]=-18; // [237] Mason et al.

H0[i]=66; PT[i]=+14; NT[i++]=-11; // [237] Mason et al.

H0[i]=70; PT[i]=+7; NT[i++]=-7; // [238] Mei et al.

H0[i]=69; PT[i]=+4; NT[i++]=-4; // [239] Tonry

H0[i]=71; PT[i]=+6; NT[i++]=-6; // [240] Willick & Puneet

H0[i]=63; PT[i]=+4.3; NT[i++]=-4.3; // [241] Xiao-Feng et al.

H0[i]=69; PT[i]=+4; NT[i++]=-4; // [242] Ferrarese et al.

H0[i]=68; PT[i]=+2; NT[i++]=-2; // [243] Gibson et al.

H0[i]=71; PT[i]=+6; NT[i++]=-6; // [244] Mould et al.

H0[i]=68; PT[i]=+6; NT[i++]=-6; // [244] Mould et al.

H0[i]=71; PT[i]=+4; NT[i++]=-4; // [245] Sakai et al.

H0[i]=77; PT[i]=+7; NT[i++]=-7; // [246] Tikhonov et al.

H0[i]=69; PT[i]=+12; NT[i++]=-19; // [247] Biggs et al.

H0[i]=69; PT[i]=+18; NT[i++]=-12; // [248] Chae KH

H0[i]=74; PT[i]=+18; NT[i++]=-17; // [248] Chae KH

H0[i]=42; PT[i]=+9; NT[i++]=-9; // [249] Collier et al.

H0[i]=73; PT[i]=+6; NT[i++]=-6; // [250] Freedman et al.

H0[i]=64; PT[i]=+8; NT[i++]=-6; // [251] Jha et al.

H0[i]=85; PT[i]=+27; NT[i++]=-23; // [252] Mason & Myers

H0[i]=61; PT[i]=+15; NT[i++]=-14; // [252] Mason & Myers

H0[i]=61; PT[i]=+23; NT[i++]=-21; // [252] Mason & Myers

H0[i]=80; PT[i]=+19; NT[i++]=-17; // [252] Mason & Myers

H0[i]=68; PT[i]=+21; NT[i++]=-19; // [252] Mason & Myers

H0[i]=71; PT[i]=+5; NT[i++]=-5; // [252] Mason & Myers

H0[i]=86; PT[i]=+24; NT[i++]=-24; // [253] Mazumdar & Narasimba

H0[i]=67; PT[i]=+7; NT[i++]=-7; // [254] Tanvir et al.

H0[i]=62.9; PT[i]=+1.6; NT[i++]=-1.6; // [255] Tripp & Branch

H0[i]=62; PT[i]=+2; NT[i++]=-2; // [255] Tripp & Branch

H0[i]=60; PT[i]=+10; NT[i++]=-10; // [256] Branch

H0[i]=66; PT[i]=+15; NT[i++]=-14; // [257] Goicoechea et al.

H0[i]=77; PT[i]=+8; NT[i++]=-8; // [258] Harris et al.

H0[i]=47; PT[i]=+23; NT[i++]=-15; // [259] Hughes & Birkinshaw

H0[i]=82; PT[i]=+8; NT[i++]=-8; // [260] Lauer et al.

H0[i]=89; PT[i]=+10; NT[i++]=-10; // [260] Lauer et al.

H0[i]=65.2; PT[i]=+1.3; NT[i++]=-1.3; // [261] Riess et al.

H0[i]=63.8; PT[i]=+1.3; NT[i++]=-1.3; // [261] Riess et al.

H0[i]=55; PT[i]=+8; NT[i++]=-8; // [262] Tammann & Labhardt

H0[i]=60; PT[i]=+6; NT[i++]=-6; // [263] Tripp

H0[i]=70; PT[i]=+5; NT[i++]=-5; // [264] Giovanelli

H0[i]=76; PT[i]=+8; NT[i++]=-8; // [264] Giovanelli

H0[i]=67; PT[i]=+8; NT[i++]=-8; // [264] Giovanelli

H0[i]=75; PT[i]=+6; NT[i++]=-6; // [265] Gregg

H0[i]=67; PT[i]=+8; NT[i++]=-8; // [266] Hjorth & Tanvir

H0[i]=70; PT[i]=+7; NT[i++]=-7; // [266] Hjorth & Tanvir

H0[i]=60; PT[i]=+40; NT[i++]=-23; // [267] Holzapfel et al.

H0[i]=78; PT[i]=+34; NT[i++]=-28; // [267] Holzapfel et al.

H0[i]=78; PT[i]=+60; NT[i++]=-40; // [267] Holzapfel et al.

H0[i]=58; PT[i]=+10; NT[i++]=-5; // [268] Hoyle et al.

H0[i]=74; PT[i]=+10; NT[i++]=-10; // [269] Schechter

H0[i]=52.5; PT[i]=+2.5; NT[i++]=-2.5; // [270] Sciama

H0[i]=54.8; PT[i]=+0.3; NT[i++]=-0.3; // [270] Sciama

H0[i]=81; PT[i]=+6; NT[i++]=-6; // [271] Tonry et al.

H0[i]=69; PT[i]=+8; NT[i++]=-8; // [272] Amendola

H0[i]=80; PT[i]=+17; NT[i++]=-17; // [272] Amendola

H0[i]=49.5; PT[i]=+4.5; NT[i++]=-4.5; // [273] Biesiada

H0[i]=65; PT[i]=+8; NT[i++]=-8; // [274] Forbes et al.

H0[i]=103; PT[i]=+59; NT[i++]=-28; // [275] Kobayashi

H0[i]=82; PT[i]=+56; NT[i++]=-24; // [275] Kobayashi

H0[i]=60; PT[i]=+24; NT[i++]=-13; // [275] Kobayashi

H0[i]=51; PT[i]=+10; NT[i++]=-7; // [275] Kobayashi

H0[i]=33; PT[i]=+22; NT[i++]=-9; // [275] Kobayashi

H0[i]=74; PT[i]=+26; NT[i++]=-15; // [275] Kobayashi

H0[i]=63; PT[i]=+28; NT[i++]=-15; // [275] Kobayashi

H0[i]=80; PT[i]=+17; NT[i++]=-17; // [276] Mallik

H0[i]=87; PT[i]=+7; NT[i++]=-7; // [276] Mallik

H0[i]=55; PT[i]=+3; NT[i++]=-3; // [277] Schaefer

H0[i]=56; PT[i]=+3; NT[i++]=-3; // [277] Schaefer

H0[i]=82.5; PT[i]=+5.9; NT[i++]=-3; // [278] Grogin & Narayan

H0[i]=82.5; PT[i]=+8.7; NT[i++]=-5.6; // [278] Grogin & Narayan

H0[i]=71; PT[i]=+30; NT[i++]=-25; // [279] Herbig

H0[i]=74.6; PT[i]=+47; NT[i++]=-33; // [280] Holzapfel et al.

H0[i]=38; PT[i]=+18; NT[i++]=-16; // [281] Jones

H0[i]=80; PT[i]=+17; NT[i++]=-17; // [282] Kennicutt Jr et al.

H0[i]=73; PT[i]=+11; NT[i++]=-11; // [283] Mould

H0[i]=81; PT[i]=+11; NT[i++]=-11; // [283] Mould

H0[i]=84; PT[i]=+16; NT[i++]=-16; // [283] Mould

H0[i]=76; PT[i]=+10; NT[i++]=-10; // [283] Mould

H0[i]=82; PT[i]=+11; NT[i++]=-11; // [283] Mould

H0[i]=71; PT[i]=+10; NT[i++]=-10; // [283] Mould

H0[i]=80; PT[i]=+17; NT[i++]=-17; // [283] Mould

H0[i]=80; PT[i]=+17; NT[i++]=-17; // [284] Nakamura & Suto

H0[i]=58; PT[i]=+6; NT[i++]=-6; // [285] Rephaeli

H0[i]=51; PT[i]=+7; NT[i++]=-7; // [286] Schaefer

H0[i]=61; PT[i]=+12; NT[i++]=-12; // [286] Schaefer

H0[i]=26; PT[i]=+5; NT[i++]=-5; // [286] Schaefer

H0[i]=69; PT[i]=+8; NT[i++]=-8; // [287] Tanvir et al.

H0[i]=78; PT[i]=+11; NT[i++]=-11; // [288] Whitmore & Schweizer

H0[i]=65; PT[i]=+25; NT[i++]=-25; // [289] Birkinshaw & Hughes

H0[i]=55; PT[i]=+17; NT[i++]=-17; // [289] Birkinshaw & Hughes

H0[i]=80; PT[i]=+17; NT[i++]=-17; // [290] Freedman

H0[i]=84; PT[i]=+5; NT[i++]=-5; // [291] Lu et al.

H0[i]=73; PT[i]=+6; NT[i++]=-6; // [292] Schmidt & Kirshner

H0[i]=90; PT[i]=+10; NT[i++]=-10; // [293] Tully

H0[i]=43.5; PT[i]=+2.7; NT[i++]=-2.7; // [294] Duemmler

H0[i]=77; PT[i]=+8; NT[i++]=-8; // [295] Lauer & Postman

H0[i]=51; PT[i]=+5; NT[i++]=-5; // [295] Lauer & Postman

H0[i]=75; PT[i]=+30; NT[i++]=-30; // [296] Leibundgut & Pinto

H0[i]=40; PT[i]=+9; NT[i++]=-9; // [297] Birkinshaw

H0[i]=45; PT[i]=+12; NT[i++]=-12; // [297] Birkinshaw

H0[i]=82; PT[i]=+7; NT[i++]=-7; // [298] Tonry

H0[i]=52; PT[i]=+2; NT[i++]=-2; // [299] Sandage & Tammann

H0[i]=45; PT[i]=+3; NT[i++]=-3; // [299] Sandage & Tammann

H0[i]=73; PT[i]=+10; NT[i++]=-10; // [300] Visvanathan

H0[i]=50; PT[i]=+10; NT[i++]=-10; // [301] Sandage & Tammann

H0[i]=52; PT[i]=+2; NT[i++]=-2; // [301] Sandage & Tammann

H0[i]=50; PT[i]=+7; NT[i++]=-7; // [301] Sandage & Tammann

H0[i]=67; PT[i]=+10; NT[i++]=-10; // [302] Dressler

H0[i]=74.3; PT[i]=+11; NT[i++]=-11; // [304] Visvanathan

H0[i]=74.3; PT[i]=+11; NT[i++]=-11; // [305] Hanes

H0[i]=76; PT[i]=+8; NT[i++]=-8; // [306] Bottinelli & Gouguenheim

H0[i]=50.3; PT[i]=+4.3; NT[i++]=-4.3; // [307] Sandage & Tammann

H0[i]=56.9; PT[i]=+3.4; NT[i++]=-3.4; // [308] Sandage & Tammann

H0[i]=57; PT[i]=+6; NT[i++]=-6; // [309] Sandage & Tammann

H0[i]=55.5; PT[i]=+8.7; NT[i]=-8.7; // [310] Sandage & Tammann

//Creates an H0 array that contains all the extremities of the tolerance ranges

for (i=0;i<=nbH0-1;i++) {

H0[i+nbH0]=H0[i]+PT[i]; Variation[i+nbH0]=PT[i];

H0[i]=H0[i]+NT[i]; Variation[i]=NT[i]; }

//Sorts H0 array in ascending order with corresponding Variation of tolerance

for (j=0;j<=2*nbH0-2;j++) {

for (i=j+1;i<=2*nbH0-1;i++) {

if (H0[i]

DH0=H0[j]; DSgn=Variation[j]; H0[j]=H0[i];

Variation[j]=Variation[i]; H0[i]=DH0; Variation[i]=DSgn; } } }

for (i=0;i<=2*nbH0-1;i++) { //Builds the nbCrossings array

if (i==0) { nbCrossings[i]=1; }

else {

if (Variation[i]<0) {nbCrossings[i]=nbCrossings[i-1]+1;}

if (Variation[i]>0) {nbCrossings[i]=nbCrossings[i-1]-1;}

if (H0[i]==H0[i-1]) {

j=i; DeltaCrossings=0;

do {

if (Variation[j]<0) {DeltaCrossings=DeltaCrossings+1;}

if (Variation[j]>0) {DeltaCrossings=DeltaCrossings-1;}

j=j-1; } while (H0[j]==H0[i]);

for(k=i;k>j;k--) {nbCrossings[k]=nbCrossings[j]+DeltaCrossings;}

} } } } //End of CreateTableOfCrossingH0Ranges

//***Function that returns the y coordinate corresponding to x for non

//***centered Gaussian curve

double GaussianCurve(double x, double Mean, double Sigma, double Multiplier) {

double y; //Coordinate y corresponding to x for a non centered Gaussian

y=(Multiplier/(Sigma*Sqrt(2*Pi)))*exp(-0.5*Sqr((x-Mean)/Sigma));

return y; } //End of CreateApproximativeCurve

//***These are the best Gaussians to fit the nbCrossing array as a function of H0***

double FindsGaussianCurvesLS(double Mean[5],double Sigma[5], double Multiplier[5]) {

int i,j; /*Counters*/ double LS = 0; //Least square

double Sum; //Sum of the 5 Gaussien curve for a specific H0 value

for (j=0;j<=2*nbH0-1;j++) {

Sum=0;

for (i=0;i<=4;i++) {

Sum=Sum+GaussianCurve(H0[j],Mean[i],Sigma[i],Multiplier[i]); }

//We give a heavier weight to any error between 69.2 and 72.1 to

// model the gap between these values

if ((H0[i]>=69.2)&&(H0[i]<=72.1)) {

LS=LS+10*(Sqr(nbCrossings[j]-Sum));}

else { LS=LS+Sqr(nbCrossings[j]-Sum); }

} return LS; } //End of FindsGaussianCurvesLS

//***This function finds the best Gaussians to fit the real curve

void FindsBestGaussiansToFitRealCurve(void) {

int i; //Counter

double DMean[5],DSigma[5],DMultiplier[5]; //Dummy arrays

double LS, PLS; //Least Square and Previous Least Square

double StepMean = 0.1, StepSigma = 0.1, StepMultiplier = 0.1; //Variations

int nbMeanNotImproved = 0; //Tells how many times not improved

int nbSigmaNotImproved = 0; //Tells how many times not improved

int nbMultiplierNotImproved = 0; //Tells how many times not improved

//Starting values (approximative values only)

Mean[0]=71; Sigma[0]=1; Multiplier[0]=-280;

Mean[1]=68; Sigma[1]=17; Multiplier[1]=3800;

Mean[2]=Mean[1]; Sigma[2]=3; Multiplier[2]=880;

Mean[3]=73; Sigma[3]=7; Multiplier[3]=1200;

Mean[4]=Mean[3]; Sigma[4]=2; Multiplier[4]=470;

//Fills the 3 dummy arrays DMean, DSigma and DMultiplier

//with the same values than the arrays Mean, Sigma and Multiplier

for(i=0;i<=4;i++) {

DMean[i]=Mean[i]; DSigma[i]=Sigma[i]; DMultiplier[i]=Multiplier[i]; }

//Tries to find the 5 best Gaussians to fit the curve

do {

for (i=0;i<=4;i++) {

//We improve Mean[i], but we force

//Mean[2] = Mean[1] & Mean[4] = Mean[3]

if ((i!=2)&&(i!=4)) {

PLS=FindsGaussianCurvesLS(Mean,Sigma,Multiplier);

DMean[i]=Mean[i]+StepMean;

if (i==1) {DMean[2]=DMean[i];}

if (i==3) {DMean[4]=DMean[i];}

LS=FindsGaussianCurvesLS(DMean,DSigma,DMultiplier);

if (LS

Mean[i]=DMean[i];

if (i==1) {Mean[2]=DMean[i];}

if (i==3) {Mean[4]=DMean[i];}

nbMeanNotImproved=0; }

else {

DMean[i]=Mean[i]-StepMean;

if (i==1) {DMean[2]=DMean[i];}

if (i==3) {DMean[4]=DMean[i];}

LS=FindsGaussianCurvesLS(DMean,DSigma,DMultiplier);

if (LS

Mean[i]=DMean[i];

if (i==1) {Mean[2]=DMean[i];}

if (i==3) {Mean[4]=DMean[i];}

nbMeanNotImproved=0; }

else {

DMean[i]=Mean[i];

if (i==1) {DMean[2]=DMean[i];}

if (i==3) {DMean[4]=DMean[i];}

nbMeanNotImproved++;

if (nbMeanNotImproved>=100) {

nbMeanNotImproved=0; StepMean=StepMean/10; } } } }

//We try to improve Sigma[i]

PLS=FindsGaussianCurvesLS(Mean,Sigma,Multiplier);

DSigma[i]=Sigma[i]+StepSigma;

LS=FindsGaussianCurvesLS(DMean,DSigma,DMultiplier);

if (LS

else {

DSigma[i]=Sigma[i]-StepSigma;

LS=FindsGaussianCurvesLS(DMean,DSigma,DMultiplier);

if (LS

Sigma[i]=DSigma[i]; nbSigmaNotImproved=0; }

else {

DSigma[i]=Sigma[i]; nbSigmaNotImproved++;

if (nbSigmaNotImproved>=100) {

nbSigmaNotImproved=0; StepSigma=StepSigma/10; } } }

//We try to improve Multiplier[i]

PLS=FindsGaussianCurvesLS(Mean,Sigma,Multiplier);

DMultiplier[i]=Multiplier[i]+StepMultiplier;

LS=FindsGaussianCurvesLS(DMean,DSigma,DMultiplier);

if (LS

Multiplier[i]=DMultiplier[i]; nbMultiplierNotImproved=0; }

else {

DMultiplier[i]=Multiplier[i]-StepMultiplier;

LS=FindsGaussianCurvesLS(DMean,DSigma,DMultiplier);

if (LS

Multiplier[i]=DMultiplier[i]; nbMultiplierNotImproved=0; }

else {

DMultiplier[i]=Multiplier[i]; nbMultiplierNotImproved++;

if (nbMultiplierNotImproved>=100) {

nbMultiplierNotImproved=0;

StepMultiplier=StepMultiplier/10; } } } }

LS=FindsGaussianCurvesLS(DMean,DSigma,DMultiplier);

} while (LS>=22000); //Sets a stop point

for(i=0;i<=4;i++) {

printf("\n Mean[%i]=%10lf Sigma[%i]=%10lf Multiplier[%i]=%10lf",

i,Mean[i],i,Sigma[i],i,Multiplier[i]);

} } //End of FindsBestGaussiansToFitRealCurve

//***The function begins by shifting Mean[1] and Mean[2] of the two

//***Gaussian curves that are around H0=69.2km/(s*MParsec) to

//***H0=72.1km/(s*MParsec) with a theoretical factor of 1.042516951. Then

//***the function that modifies the global H0 array builds the summation of

//***the 4 positive Gaussian curves.

void CreatesFinalGaussianCurve(void) {

int i; //Dummy index value

//We shift Mean[1] and Mean[2] with a theoretical factor of 1.042516951

Mean[1]=Mean[1]*1.042516951; Mean[2]=Mean[1];

//We omit i=0 to remove the negative Gaussian curve

for(i=1;i<=nbH0*2-1;i++) {

nbCrossings[i]=GaussianCurve(H0[i],Mean[1],Sigma[1],Multiplier[1]);

nbCrossings[i]=nbCrossings[i]+GaussianCurve(H0[i],Mean[2],Sigma[2],

Multiplier[2]);

nbCrossings[i]=nbCrossings[i]+GaussianCurve(H0[i],Mean[3],Sigma[3],

Multiplier[3]);

nbCrossings[i]=nbCrossings[i]+GaussianCurve(H0[i],Mean[4],Sigma[4],

Multiplier[4]); } } //End of CreatesFinalGaussianCurve

//***Returns the Least Mean Square of the equation DA*x^2+DB*x+DC.

double FindsLMS(double DA, double DB, double DC) {

int i; /*Dummy index value*/ double LMS = 0; //Least Mean Square

for(i=PosTipIndex-n;i<=PosTipIndex+n;i++){

LMS=LMS+Sqr(nbCrossings[i]-(DA*H0[i]*H0[i]+DB*H0[i]+DC));

} return LMS; } //End of FindsLMS

//***Returns A, B, and C coefficients of the quadratic equation of the tip

void ApproximatesTipEquation(void) {

double LMSTipMin; //Reminds the lowest value of least mean square

double PLMSTip, NLMSTip; //LMSTip for a forward and backward step

double DA,DB,DC; //Dummy values of A, B and C coefficients

double StepA,StepB,StepC; //Step variation of the coefficients

double nbCrossingsMax = 0; //Maximum number crossings at the tip

double xa,xb,xc,ya,yb,yc; /*3 coordinates*/ int i; //Dummy index value

for(i=0;i<=nbH0*2-1;i++) { //Finds the index of the approximated tip

if (nbCrossings[i]>nbCrossingsMax) {

nbCrossingsMax=nbCrossings[i]; PosTipIndex=i; } }

i=PosTipIndex; xa=H0[i-n]; xb=H0[i]; xc=H0[i+n];

ya=nbCrossings[i-n]; yb=nbCrossings[i]; yc=nbCrossings[i+n];

//Sets coefficients ABC

A=((yc-ya)/((xc-xa)*(xc-xb)))-((yb-ya)/((xb-xa)*(xc-xb)));

B=((yb-ya)/(xb-xa))-A*(xb+xa); C=ya-A*xa*xa-B*xa; DA=A;DB=B;DC=C;

} //End of ApproximatesTipEquation

//***Function that returns the Best estimate of H0

void BestEstimateOfH0(void) {

CreateTableOfCrossingH0Ranges(); FindsBestGaussiansToFitRealCurve();

CreatesFinalGaussianCurve(); ApproximatesTipEquation();

BEH=-B/(2*A); Accuracy_ppm = (BEH-TVH)/(TVH*1E-6);

LMSTip=FindsLMS(A,B,C); } //End of BesEstimateOfH0

int main(void) {

BestEstimateOfH0();

printf("\n\n Equation of the tip: y = %.10lfx^2+ %.10lfx+ %.10lf",A,B,C);

printf("\n Best estimate of H0 = %.10lf km/(s*MParsec)",BEH);

printf("\n Theoretical H0 = %.10lf km/(s*MParsec)",TVH);

printf("\n Relative accuracy versus theoretical value = %.10lf ppm",

Accuracy_ppm);

getchar(); return 0;

} //End of main